简介：通过地震数据获取裂缝储藏中流体的性质并对流体类型进行识别,是地震勘探岩性反演的重要问题之一。由于地震波的速度、储层的密度等弹性参数对某些流体不具有很强的敏感性,使只依赖振幅信息进行流体识别的传统AVO方法面临困境。作为传统叠前振幅反演的一个拓展,频变AVO（FDAVO）技术进一步考虑了振幅对频率的依赖关系,将这种依赖关系与地下裂缝结构、流体填充对应起来,能带来更丰富的流体信息。利用该技术,本文提出了一种基于地震数据参数化Chapman模型的贝叶斯反演新方法（BIDCMP）,它包含两步算法,即,FDAVO反演储层的非弹性属性和贝叶斯框架下的流体识别。首先,通过匹配观测数据和模型数据,构造差函数反演裂缝储层非弹性参数。随后,在贝叶斯框架下,使用马尔科夫随机场（MRF）作为先验模型,联合多参数场识别流体。本方法在计算过程中,除综合考虑了弹性参数场、测井资料等常规信息外,还特别地加人了第一步中反演得的非弹性参数的约束,从而充分利用了流体粘性差异,最后在最大后验概率（MAP）准则下输出最佳岩性一流体识别结果。分别对合成地震记录和模拟岩性—流体剖面验证本文方法的有效性,结果证明本文方法获得的流体识别结果准确可信。

简介：本文介绍了在B隧道勘查中应用高分辨电磁成像技术探测隐伏构造的一个成功实例.数据采集采用高频电磁成像系统(STRATAGEMEH4,频率范围从1Hz到90kHz).野外采集正交的电磁场分量,从地磁成像剖面中提取目标体有关电磁学信息.为获取高质量的采集数据,野外采用包含天然场源和全张量可控源的混合场源.B隧道坐落在位于中国中部的湖北省的西部,隧道埋深不到200m,但区域地质作用导致它的地质机构十分复杂.第一次勘查过程中,勘探人员误把一个脱落体的露头当成基岩面,施工过程中出现了冒顶.第二次勘探时采用高分辨电磁法和折射地震法,这次勘探找到了隐伏的基岩面和一个隐伏断裂.勘探结果与后来隧道挖掘揭示的构造吻合.

简介：基于Kelvin横向各向同性（KEL-VTI）介质的复物性参数矩阵,在弱各向异性和弱衰减的假设下,提出KEL-VTI介质各向异性复相速度和品质因子的近似解;结合KEL-VTI介质模型,讨论了qP、qSV、qSH波的相位与能量的传播特点;进一步针对淮南煤矿的地质情况建立了典型的KEL-VTI介质模型,数值模拟结果显示PP、PSV波的理论波场与该地区实际采集的三分量地震的纵波和转换波波场吻合程度较好,说明KEL-VTI介质假设是对这种典型煤田地震地质条件的较好近似,有助于多分量地震数据的吸收衰减补偿研究。

简介：本文以中观孔隙结构的White模型为基础，构建了部分饱和孔隙介质模型，利用Biot方程的建立思路和Johnson推导的体变模量，推导了部分饱和孔隙介质中的纵波方程，并以平面波为例，求取了方程的衰减系数，分析了地震勘探频带范围内地震波的衰减特性。结果表明：在部分饱和孔隙介质中，地震波在低频段也会发生明显的衰减和频散现象，频率越大，衰减越大；且第二纵波的衰减比第一纵波更为明显；这一结论弥补了Biot理论在描述地震勘探频带范围内波的衰减现象的不足。文中还研究了孔隙度、饱和度和模型内径尺寸对第纵波衰减特性的影响机理，主要表现在在地震勘探频带范围内，波的衰减随孔隙度的增大而增大，随含油气饱和度的增大而减小，当孔隙内径尺寸小于二分之一外径尺寸时，波的衰减随内径尺寸的增大而增大，当内径尺寸大于二分之一外径尺寸时，波的衰减随内径尺寸增大而减小。

简介：地震正演模拟技术是研究地震波在复杂介质中传播规律的有效途经，尤其在地质构造及其复杂的中国西部地区，其意义更是重大。本文介绍了两种新的正演模拟技术：有限元有限差分方法（FE—FDM）和任意精细积分方法（ADPI），并结合实例分析来验证FE—FDM和ADPI算法的实际效果，结果证明这两种方法能够有效地适用于复杂介质下的地震波传播性质的研究。

简介：白云岩成岩收缩晶间孔、洞、缝与构造网状缝相互沟通可以组成良好的油气储层,但是这种复杂裂缝-孔隙型储层分布随机性强、发育程度和差异性大,储层预测难度大.本文以中国西部酒泉盆地青西油田下白垩统下沟组湖相白云岩裂缝-孔隙型储层为例,提出地球物理综合预测碳酸盐岩裂缝的方法.本文描述了综合地球物理方法预测碳酸盐岩裂缝储层的实施和应用效果.

简介：根据Chapman理论模型,在各向异性介质（如HTI介质）中,当入射角在0-45。范围内,慢横波会发生较大的衰减和频散,且对流体粘度敏感,而P波和快横波则比较小。对于沿裂隙法向传播的慢横波,其振幅受流体影响很大。因此,在P波响应对流体不敏感的情况下,可利用慢横波来获得裂隙型油气藏的流体信息。本文分析了胜利油田垦71地区三维三分量地震数据,检测出的慢横波振幅和旅行时异常与该区的测井资料十分吻合。分析结果还发现,与含油区相比,含水区会产生更高的横波分裂。在含水区,慢横波振幅会产生明显变化,而在含油区则几乎没有变化。

简介：多分量勘探相比常规纵波勘探而言有许多的优势，PP／PS联合AVO分析和反演是强而有效的储层识别方法。在本文中，我们推导了更精确的转换波AVO反射系数公式；并在入射角小于30度时进行了进一步简化，其对于转换波AVO分析和反演而言更加简单有效。基于该近似，我们进行了PP／PS联合AVO反演。实际资料的实例表明，反演得到的纵横波速度比可很好地识别岩性及油气。反演得到的流体因子和泊松比等其他属性在储层处也显示出明显异常，效果显著。

简介：通过常规测井曲线小波变换来评价裂缝参数是一个较新的研究课题。本文针对火山岩裂缝储层,对常规测井曲线做coif5、bior4.4和db5基小波的小波变换,通过变换后的分解信号与成像测井裂缝密度的对比研究,寻找出与裂缝密度相匹配的基小波,再将常规测井曲线在这个基小波下的分解信号与曲线变化率法相结合建立裂缝指示曲线来识别裂缝发育段,最后通过裂缝指示曲线与成像测井的裂缝密度之间的量化关系来比较准确的评价裂缝密度。将裂缝指示曲线法在松辽盆地南部进行了应用,通过裂缝指示曲线求取的裂缝密度与成像测井的裂缝密度之间的相对误差为0.53。

简介：瑞利波勘探主要是建立在弹性介质分层半空间模型基础上的。当实际地层中包含孔隙介质层时，需要将孔隙介质简化为弹性介质来进行分析和处理。这种简化处理究竟对瑞利波勘探造成怎样的影响是本文所研究的问题。本文基于弹性介质与孔隙介质共同构成的分层半空间模型，首先推导了分层半空间中两种介质处于不同相对位置时的瑞利波频散方程，解决了不同阶数矩阵之间的变量传递问题然后，针对传递矩阵法在求解频散函数时可能出现的溢出问题，给出了一种可以有效提高计算范围的解决方案；同时，提出了一套新的数值算法用于复频散方程的快速求解。数值计算结果表明：当孔隙介质位于半空间表面时对低频瑞利波频散特性的影响最为显著，而在其它情况下对瑞利波频散特性的影响相对较小。

简介：非裂化的完全匹配层(NPML)吸收边界条件首先用于有限差分方法模拟固体中弹性波的传播。本文将NPML加以推广用来模拟孔隙介质中弹性波传播。采用空间四阶精度和时间二阶精度中心差分将Biot方程离散化到交错网格。我们选取柱状双层模型来验证NPML有效性。并将NPML应用到储层环境下裸眼井模型,数值结果和离散波数法得到结果进行了对比,符合较好。

简介：针对地面场源在地层介质中所产生瞬变电磁场的数值计算方法和响应分析问题，本文展开研究．解法方面，传统离散镜像法多采用复数运算、数字滤波等方法计算量大，针对这些问题，提出一种改进的离散镜像方法：基于Gaver-Stehfest概率变换算法将电磁场解式实数化，选用Prony方法对目标核函数进行指数级数逼近，根据离散镜像原理和近似系数闭合式求解瞬变电磁场．通过试算均质模型中瞬变电磁场并对比汉克尔变换的数字滤波法所得结果，证明该方法有效且具有较好的精度和适用性．继而基于该方法计算地表磁偶源在典型地电模型中产生的瞬变电磁场，对“地面激发一地层中测量”方式所得感应磁场水平分量响应进行分析并得出结论瞬变场水平分量响应与地电结构、观测时间、空间位置等因素有关，感应磁场水平分量响应反映出涡流场分布及其垂向梯度变化情况，在探测异常体的工作中应尽量选取零偏移距、较大偏移距位置钻孔或用较大观测延时以减小背景场对勘测结果的影响。文中所用离散镜像方法与正演计算结论可为相关研究工作提供参考依据。

简介：本文是三分量地震在陆相薄互层地震勘探中应用成功的首例。大庆长垣喇嘛甸油田处于高含水开发后期,储层预测的主要难题是密井网条件下厚度2m以上砂体的边界识别问题,常规单一纵波地震解释存在多解性,引入横波信息可增加地震解释的可靠性。本文根据纵、横波对气藏的响应特征,利用纵、横波联合振幅属性、分频属性对气藏进行了分析,能够很好地反应气藏的边界,体现了横波受流体影响较小的优势。通过测井曲线岩石物理定量的分析了该地区主要是密度、λρ和μρ三种物性参数能够很好地反映储层岩性变化。在以地震数据分频为核心的高精度地震参考标准层解释基础上利用纵、横波联合反演得出的纵、横波速度比、λρ、μρ等属性与测井解释的储层具有很好的一致性,多波属性λρ沿层切片描述了开发小层的砂体平面展布,与井砂岩图具有较高的一致性,刻画了井间砂体边界变化的细节,并指示了剩余油挖潜的有利区域。

简介：应用多分量地震资料进行成像时通常需要先做波场分离,然后再对分离的波型进行成像。其中,波场分离可以在空间域或波数域实现。然而,由于用交错网格有限差分进行弹性波场数值模拟时,用来进行波数域波场分离的质点振动速度分量定义在不同网格节点上,本文提出了利用波数域插值方法来估算同一网格节点所需质点振动速度值;进而给出了先进行波数域插值后进行波场分离的波数域保幅波场分离方案。数值实验结果表明波数域插值方法具有较高的插值精度且保幅波场分离方法具有较好的保幅性,将本文方法进一步应用于弹性波逆时偏移可以获得保幅性较好的成像结果且对存在一定程度速度误差情况具有较好的适应性。

简介：这篇文章重点研究改进的Gabor小波（improvedGaborwavelet，IGW）变）并讨论了它在地震信号处理和解释中的应用ThispaperintroducesanimprovedGaborwaveletanditscompletetransform，andmainlyanalysestheirpropertiesanddiscussesapplicationsofthesepropertiesinseismicsignalprocessandinterpretation。改进的Gabor小波变换具有以下特性：1）IGWT把时域信号映射到时间一频率域，而传统Gabor小波变换把时间信号映射到时间一尺度域；2）IGWT可用于信号分频，通过固定变换的主频参数dominantfrequency，并变换能提取相应的子带信号，且其主频部分的信息与原信号相应频率部分的信息一致，通过调节变换的分辨率因子，变换能有效控制子带信号的带宽；3）用IGWT和IGWIT构建的滤波器有良好的时一频局部性，在指定时一频范围内能实现针对性滤波。文章用仿真实验和实际用例验证IGWT的这些特性，并在提高地震信号分辨率、地震信号分频和识别小断层等地震信号处理和解释等方面的应用中取得良好效果。

简介：地震波场数值模拟中不可避免地会出现边界反射,一般采用吸收边界条件以压制人工边界反射。目前常用的分裂式完全匹配层（PML）边界条件需要在边界处进行特殊处理,尤其是在三维情况下需要将变量分裂为三个分量,增加了数值模拟的计算时间和内存占用量。与分裂式PML吸收边界条件相比,混合吸收边界条件（HABC）具有易于实现、计算量小和吸收效果好等优点,可以提高三维波动方程数值模拟的计算效率。本文将基于一阶Higdon单程波方程的混合吸收边界条件从二维计算域发展到三维,提出了适用于三维弹性波数值模拟的混合吸收边界条件。均匀模型以及复杂模型的三维数值模拟结果表明,混合吸收边界条件与传统的完全匹配层边界条件相比,具有效率高、吸收效果好的优势。

简介：本文提出了一种基于Kirchhoff积分偏移和逆时偏移的联合速度分析及成像方法,采用剩余曲率分析及层剥离策略进行偏移速度建模。本文方法改善了Kirchhoff积分偏移在复杂构造时计算精度不高和逆时偏移计算效率慢的缺点,兼有计算效率高、成像精度高的优点,并将其应用在反射波法隧道超前预报中。通过模型试算,发现隧道中使用逆时偏移的成像结果在多方面优于Kirchhoff积分偏移的成像结果;通过对实测数据的处理,验证本方法计算效率较高,建立的速度模型合理,成像剖面清晰,结合地质调绘资料可对隧道开挖前方地质构造做出较准确预报。

简介：现阶段的叠前地震反演技术中用于描述反射系数与纵、横波速度和密度之间的关系几乎完全是Zoeppritz方程的近似式，由于这些近似公式在大角度和弹性参数变化剧烈时误差较大，这不仅降低了反演解的精度，而且增加了叠前反演的多解性。本文探索了直接利用Zoeppritz方程求解精确反射系数的理论方法，并基于广义线性反演理论详细推导了基于叠前大角度地震资料的纵、横波速度和密度三参数同步反演算法，同时在反演过程引入正则化约束阻尼因子和共轭梯度算法，有效降低了反演的不适定性和提高了反演收敛性。理论模型试算和实际工区应用表明，本文提出的反演方法能够有效利用大角度（一般入射角〉30°）的叠前地震数据，获得更精确的地震弹性参数反演结果，并且反演结果忠实于地震资料，与井吻合较好。

简介：基于PetrJilek正交各向异性介质中的P-SV波反射系数精确公式,通过近似得到VTI和HTI介质中的P-SV波近似公式,同时结合Kenneth横波弹性阻抗原理,推导建立了VTI和HTI介质中P-SV波的弹性阻抗近似公式;在此基础上构建了弱各向异性介质中的流体识别因子公式,通过理论模型对Castagna和Smith的岩性组合进行了识别;最后详细分析了各向异性参数的取值对P-SV波弹性阻抗及识别因子产生的影响及变化趋势。

简介：阵列声波信号是典型的非线性、非平稳信号，Hilbert～Huang变换（HHT）是处理非平稳信号的一种比较新的时频分析方法。通过对信号进行经验模态分解（EMD）和对瞬时频率的求解，可以获得声波信号的时一频谱。其关键技术就是进行经验模态分解，任何非平稳的信号都可以分解为有限数目并且具有一定物理意义的固有模态函数。EMD方法可以理解为以声波信号极值特征尺度为度量的时频滤波过程。滤波器充分保留了声波信号本身的非线性和非平稳特征，在声波信号的滤波和去噪中具有很大的优势。文中介绍了HHT时频滤波的实现过程，并列举了一些声波测井波列实例，说明了该方法的有效性。