简介:这篇文章重点研究改进的Gabor小波(improvedGaborwavelet,IGW)变)并讨论了它在地震信号处理和解释中的应用ThispaperintroducesanimprovedGaborwaveletanditscompletetransform,andmainlyanalysestheirpropertiesanddiscussesapplicationsofthesepropertiesinseismicsignalprocessandinterpretation。改进的Gabor小波变换具有以下特性:1)IGWT把时域信号映射到时间一频率域,而传统Gabor小波变换把时间信号映射到时间一尺度域;2)IGWT可用于信号分频,通过固定变换的主频参数dominantfrequency,并变换能提取相应的子带信号,且其主频部分的信息与原信号相应频率部分的信息一致,通过调节变换的分辨率因子,变换能有效控制子带信号的带宽;3)用IGWT和IGWIT构建的滤波器有良好的时一频局部性,在指定时一频范围内能实现针对性滤波。文章用仿真实验和实际用例验证IGWT的这些特性,并在提高地震信号分辨率、地震信号分频和识别小断层等地震信号处理和解释等方面的应用中取得良好效果。
简介:在地震频段内,中观尺度局域流是引起孔隙岩层速度频散与衰减的主要原因,研究中观尺度孔隙中局域流的过渡频率及其对应的衰减(逆品质因子)峰值频率之间的关系,有助于深入分析孔隙介质内部构造的细节。基于Biot孔隙弹性理论得出时间域内局域流流体通量,引入傅立叶变换得到频域上的流体通量。为避免不同单元体之间局域流的相互影响,建立双相孔隙介质单重单元体模型,并且在流体通量基础上选择含水层的过渡频率来等价局域流的过渡频率。此外,针对不同孔隙介质参数,分析了局域流整体过渡频率及相应的衰减峰值频率的变化特征。为了了解边界条件对计算结果造成的影响,建立了多重单元体模型,引出局部局域流过渡频率及相应的衰减峰值频率。数值模拟结果表明,在低饱和度状态下,两种频率的变化趋势相同,但变化幅度存在较大差异;在高饱和度状态下,两种频率的变化趋势与变化幅度都很接近。
简介:噪声衰减是探地雷达信号处理中的关键问题之一。当探测目标埋藏深度比较浅时,其反射信号与直耦信号和地面回波信号相互重叠,直接影响目标反射波到达时刻的检测及目标的正确定位。针对这个问题,本文提出了一种基于Curvelet变换的噪声衰减方法。通过对理论数值模拟数据和实测数据的处理,以及与平均消去法和二维连续小波该方法处理结果的对比,验证了该方法的可行性和有效性。处理结果显示,该方法不仅可以去除背景噪声、同时可以衰减倾斜相关的相干干扰和数据中的随机噪声。与二维连续小波变换方法相比有更高的计算效率。
简介:噪声衰减是探地雷达信号处理中的关键问题之一。当探测目标埋藏深度比较浅时,其反射信号与直耦信号和地面回波信号相互重叠,直接影响目标反射波到达时刻的检测及目标的正确定位。针对这个问题,本文提出了一种基于Curvelet变换的噪声衰减方法。通过对理论数值模拟数据和实测数据的处理,以及与平均消去法和二维连续小波该方法处理结果的对比,验证了该方法的可行性和有效性。处理结果显示,该方法不仅可以去除背景噪声、同时可以衰减倾斜相关的相干干扰和数据中的随机噪声。与二维连续小波变换方法相比有更高的计算效率。更多还原
简介:希尔伯特—黄变换(Hilbert-Huangtransform,HHT)是一种新的适合非平稳和非线性信号的分析方法,由于地震信号一般呈现出非平稳与非线性特性,因此HHT非常适合地震信号的分析。本文首先介绍了HHT中关于经验模态分解(EmpiricalModeDecomposition,EMD)的实现过程,在此基础上分析了几种基于EMD获得本征模态函数(IntrinsicModeFunctions,IMF)来计算瞬时频率的算法,其中利用了两个采样间隔瞬时频率的平均来计算瞬时频率,较好地反映了地震信号频率成分随时间变化的特征。将该方法应用于四川东北部某地区海相碳酸盐岩地层三维地震叠后偏移数据处理,提取"三瞬"地震属性,与传统的希尔伯特变换提取的"三瞬"地震属性进行对比,结果表明基于HHT的"三瞬"地震属性结果具有更高的分辨率,IMF2的瞬时相位能够较好地刻画台地边缘生物礁相,IMF2的瞬时频率亦具有较好的分带性。将IMF2的"三瞬"地震属性与钻井等资料结合分析,能够更好地识别沉积相的分布。
简介:常规地震勘探大多基于反射波理论,对尺度小、倾角陡、构造复杂的灾害性地质体,很难得到精确的成像结果。本文建立2类典型矿井地质模型,采用有限差分正演模拟算法,研究了矿井地震散射波的特征,将等效偏移距偏移(EOM)和基于等效偏移距干涉叠加偏移的散射成像方法用于矿井地震勘探中,重点分析了散射波成像相比常规反射波地震成像的技术优势。研究表明:1)散射成像方法能够提高覆盖次数,充分利用有效的弱散射信号,对不均匀复杂地质体的散射波的成像效果明显优于反射波叠后偏移;2)对于巷道超前探测,它弥补了基于射线理论的叠加偏移的不足,为提高矿井地震勘探分辨率提供了有效的成像方法。
简介:裂缝识别对碳酸盐岩储层测井评价至关重要,传统的测井仪器由于探测深度浅(小于3米)而无法对井壁外围(大于3米)裂缝发育情况进行评价,远探测声波测井仪器采用相控阵发射、同时加大源距,有效提升了测井仪器的探测深度。但由于缺少正演模拟研究,目前对于裂缝的解释往往是基于经验而缺乏理论依据,导致很多现象难以解释。本文利用高阶有限差分方法对远探测声反射波测井裂缝识别进行了正演模拟及叠前逆时偏移成像研究,首先在理论研究的基础上构建正演模型,重点研究不同裂缝离井壁距离、裂缝张开度和倾角的响应特征;其次在单因素变化基础上提取成像区域的能量强度,分析确定出在实际地层速度有波动变化时远探测声反射波测井方法能够识别裂缝的条件;最后通过对裂缝识别的影响因素定量化分析,确定裂缝识别的最大距离、最小张开度和最小倾角,降低了裂缝识别与评价中的多解性。研究成果对远探测声反射波测井仪器的发展、数据处理方法的改进,以及后续的测井解释工作都有一定的借鉴意义。
简介:震源定位是微地震监测关键技术之一。本文提出用于微地震定位的弹性波和加权弹性波(WEW)干涉成像方法。该方法在保证定位精度的同时,还可避免震源假象。通过各向同性水平层介质状模型的数值试验,初步表明该方法可适应低信噪比微震信号、速度随机扰动、较稀疏的检波器分布等情况,并在速度模型存在一定的系统误差时也仍保持较高的定位精度。由于干涉成像方法不需要进行初至拾取,定位效率相对传统走时方法也得到了提高。采用二维断层模型试算Nnumericalresultsofusingatwo-dimensionalfaultmodel,表明方法还能实现多震源定位,且比逆时成像有更高的定位精度。
简介:为了探明介质的非完全弹性对瑞利型槽波传播的影响以及利用瑞利型槽波品质因子QR反演煤层碳化程度、裂隙、煤层厚度等煤层属性,在本文中,我们将复速度引入瑞利型槽波频散方程,计算了煤层瑞利型槽波的品质因子QR曲线;采用控制变量法,分析了瑞利型槽波品质因子QR随煤层厚度、煤层纵横波品质因子以及围岩纵横波品质因子的变化规律。研究表明:瑞利型槽波品质因子QR曲线与群速度曲线变化趋势一致;埃里相(Airy-phase)频率附近瑞利型槽波的品质因子QR最小,且埃里相频率随煤层厚度增加而减小;瑞利型槽波品质因子QR随煤层横波品质因子QS2的增大而增大。利用瑞利型槽波品质因子QR曲线可以实现对瑞利型槽波的吸收衰减进行补偿、对煤层厚度进行预测以及岩性成像。更多还原
简介:地震波场数值模拟方法对理解和分析地震波的传播规律具有着重要的意义。弹性波动方程能够模拟地下介质的实际情况,为偏移和成像提供有效的依据。在弹性波波场数值模拟中,旋转交错网格数值模拟(RSM)修改了标准交错网格数值模拟(SSM)方法,将同类的参数定义在同样的节点上,拓宽了稳定性条件的约束,但在低速区会出现较严重的频散。变阶数差分方法是自适应空间算子长度方法的一种变化和推广。它以理论频散误差研究为基础,结合实际波场传播的情况进行误差计算,对不同速度匹配不同的差分阶数。本文研究了变阶数旋转交错网格数值模拟(VRSM),即是籽变阶数方法应用到RSM中,它可以很好地解决RSM在低速区域的数值频散问题,以及减少不必要的时间损耗;同时讨论了旋转交错网格的理论频散特性,并基于波场分离的方法分析了实际波场传播的频散误差,将原方法的应用范围由声波推广到剪切波,由理论值推广到时变值。在数值模拟试验中,VRSM将被应用于水平层模型和Overthrust模型。通过阶数分配以及相应波场传播效果和计算时间的分析,验证了该方法应用于复杂介质波场模拟中的实用性和有效性。实验的结果表明VRSM能够合理分配不同速度所对应的差分阶数,能保证计算的精确性,并合理控制计算的时间。
简介:传统的AVO正演只考虑了单。界面的反射系数对地震波波场振幅的影响,忽略了地震波在介质中传播的各种传播效应。通过引入地震波在介质中传播的几何扩散、吸收衰减以及透射损失等传播效应,提出了基于射线理论的水平层状介质多波保幅AVO正演方法。推导了水平层状介质多波几何扩散校正公式,来描述多波在介质中传播的几何扩散效应。通过直接引入复旅行时,而无需借助复速度,建立了复旅行时与品质因子的关系,来描述粘弹介质的吸收衰减。直接求解Zoeppritz方程计算多波的透射系数,用于描述多波在介质中传播时的透射损失。数值计算表明,几何扩散、吸收衰减以及透射损失对多波振幅的影响是随偏移距变化而变化的,多波保幅AVO正演需要考虑波传播效应对反射波振幅的改造。
简介:常规多波联合反演采用Zoeppritz方程的近似式构建正演方程,反演过程中需要假定背景纵横波速度比为常数,其反演精度不高,稳定性不好。本文提出了一种基于精确Zoeppritz方程的多波联合反演方法,结合贝叶斯方法进行广义线性反演。本方法基于精确Zoeppritz方程构建正演方程,避免了近似式反演在大角度时引起的误差;利用贝叶斯方法引入模型参数的先验分布信息,作为反演的正则化项,降低了反演的不适定性;反演目标函数中引入低频软约束,稳定了反演低频结果,提高了反演的鲁棒性;在求解反演目标函数时,利用快速算法,降低了反演的运算量。经过模型试算,证明了该方法的优越性和抗噪性;并在实际资料的应用中证明了该方法的实用性和有效性。
简介:分数阶S变换(FRST)具有较强的时频聚集性。利用FRST处理地震数据,通过合适的分数阶参数将频率轴旋转到适当位置,即可实现目标地质特征信息的最佳识别。由于不同的地震信号的最优分数阶参数可能不同,因而对整体的分数阶参数的最优估计不利于对多道地震数据的处理。本文首先利用FRST分离出共频率数据体,并利用共频率数据体进行了低频伴影分析,然后提出FRST和盲分离结合的方法,不需要对地震数据的最优分数阶参数进行估计,即可提取识别有效地质特征信息的独立频谱,提高对地震数据的解释效率。仿真实验表明在分数阶时频域内此方法能有效分离出独立的频率信息。将该方法用于实际的地震数据,并与已知井信息进行比对,验证了其有效性。
简介:针对转换波各向异性叠前时间偏移实践中的四参数速度模型估计中参数间相互影响、难于准确的确定,本文将简化的转换波两参数动校正方程从叠加速度分析拓展到叠前时间偏移速度模型修正,形成了基于两参数动校正方程的新的四参数法转换波叠前时间偏移速度修正新方法和流程。在该方法和流程中,先进行转换波两参数叠加速度分析,获得转换波叠前时间偏移初始的速度和各向异性参数,然后通过分析偏移后的共成像点道集中的剩余时差修正速度和各向异性参数。其中叠前时间偏移速度模型的垂直速度比仍是需要利用纵波和转换波叠加剖面,层位标定后,利用相关的方法获得。有效速度比需结合纵波速度分析,利用Thomsen(1999)公式计算得到初始模型,最终依靠百分比扫描偏移处理获得最终的有效速度比模型。该方法简化了转换波高质量成像的速度估计方法,减小了多个参数估计的不确定性,在实际应用中也取得了较好的成像效果。