基于能量密度的自解耦互相关成像条件

基于能量密度构建的弹性波能量成像条件,可以准确提取弹性波场的能量信息,并可以有效压制背向散射.但是该成像条件得到的成像结果将所有弹性波场能量信息耦合在一起,难以区分纯波模式的能量信息.为此,我们从势能及动能两方面将能量场解耦,得到纯纵波（PP）、纯横波（SS）、转换波（C）能量场.根据能量守恒定律,基于Helmholtz原理对位移场的空间导数进行因式分解,将弹性势能自解耦得到纯纵波、纯横波及转换波势能;引入体应力构建得到P波及S波速度,将弹性动能分解为纯纵波、纯横波和转换波动能.从而,构建基于能量密度自解耦得到纯波能量互相关成像条件.数值实验表明自解耦成像条件可实现弹性波场能量的分解,并得到背向散射压制、振幅有效保持的PP波、SS波和C波成像结果.     

基于保幅波场分离的弹性波逆时偏移方法研究（英文）

欲实现基于弹性波方程的矢量波场逆时偏移纵、横波独立成像,必须在波场延拓过程中实现纵、横波场的分离,散度和旋度算子分离的纵、横波出现振幅与相位的畸变,导致输出成像结果的振幅失真。本文提出一种在弹性波场延拓过程中实现纵、横波保幅分离的方法,在传统的弹性波方程中加入纵波压力、纵波振动速度和横波振动速度方程,实现纵横波的矢量分解,再对分解后的矢量纵波和矢量横波做标量化合成得到保幅分离的纵、横波场,对保幅分离的纵、横波场应用成像条件,然后实现矢量波场逆时偏移的保幅纵横波成像。该方法可以保证分离后纵、横波的振幅与相位不变;同时,分解后的纵波压力和纵波振动速度可用于层间反射噪音压制和横波极性校正,提高多分量地震资料联合逆时偏移的纵、横波成像质量,从而实现保幅弹性波逆时偏移的目的,为叠前深度剖面应用于叠前反演工作奠定基础。     

矢量分离纵横波场的弹性波逆时偏移

弹性波逆时偏移是当前多分量地震资料相对准确的偏移算法,它能够形成多波模式的成像剖面,从而减少纵波勘探的多解性.本文首先依据各向同性介质中矢量分离纵横波场的速度-应力方程组,利用高阶交错网格有限差分数值方法求解弹性波方程,进而构建矢量的纵波和横波波场,不同于散度和旋度算子分离纵横波场的传统方法,文中提出的矢量分离纵横波场方法保持了原始波场的振幅和相位特征.文中也提出将震源归一化的内积成像条件应用于分离后的纯纵波和横波矢量场,由此得到的转换波成像避免了传统弹性波成像方法中出现的极性反转.水平层状和复杂构造模型测试表明,文中提出的基于矢量分离纵横波场的弹性波逆时偏移方法成像精度高,转换波成像PS和SP极性无反转,所形成的多种模式纯波剖面能够准确地对复杂地下构造成像.     

二维TI介质非结构网格弹性波矢量逆时偏移

波场延拓得到的多分量波场中既包含纵波信息也包含横波信息,能否在全波场中实现纵横波的分离对各向同性和各向异性逆时偏移都有非常重要的意义.传统的散度旋度分离只适应于各向同性介质而对各向异性介质却无效.在非规则、非结构网格的弹性波数值模拟方法的基础上,发展了一种适应于各向异性介质的波场分离方法.该方法通过求解Christoffel方程,得到相角和极化角的关系,再利用群角和相角的关系,直接得到群角和极化角的关系.该方法与现存的各向异性波场分离相比,获得的计算效率改进更显著,而且存储量小.用简单各向异性模型和SEG各向异性Hess模型进行测试,都得到了较好的效果,证明了本文方法的有效性.     

等效介质理论中裂缝密度的适用性研究

为了比较裂缝等效介质理论的预测精度,本文利用物理模拟的方法制作了一组9个裂隙密度模型,其中裂隙密度从0增大到12%,裂隙厚度和直径保持不变,分别为0.12mm和3mm.使用0.5MHz的纵横波换能器,采用脉冲透射法得到纵横波的速度和各向异性.然后,将基质的模量、裂缝填充物模量和裂缝参数代入等效介质理论中,计算出纵横波速度和各向异性随着裂隙密度的变化情况.通过纵横波速度的实测值与等效介质理论的定量对比,发现当裂隙密度小于12%时,Hudson理论和线性滑动理论对纵波速度的预测与实测值非常吻合.当裂隙密度小于8%时,Hudson理论和线性滑动理论对快横波速度的预测与实测值非常吻合.当裂隙密度小于9%时,各向异性自洽理论和NIA(非相互作用近似)理论对慢横波速度的预测与实测值非常吻合.通过纵横波速度各向异性实测值与理论计算值定量对比,发现当裂隙密度小于12%时,Hudson理论和线性滑动理论对纵波速度各向异性的预测与实测值更吻合;当裂隙密度小于6%时,各向异性自洽理论和NIA理论对横波速度各向异性的预测与实测值更吻合.     

VTI介质SH-SH波地震反演方法研究

近年来,随着横波可控震源技术的发展,国内外已经实现了纯横波地震勘探.相较于地震纵波,地震横波对横波阻抗、横波速度尤其是各向异性参数变化更为敏感,因此地震横波能够用来更好地估算上述地层参数.VTI(具有垂直对称轴的横向各向同性)是地层介质中广泛存在的一种各向异性形式,对振幅随偏移距变化(AVO)影响显著.本文提出了一种改进的VTI介质SH-SH波反射系数近似公式,新公式具有较高的精度且仅包含两项待求参数：横波阻抗和水平横波(SH波)速度.基于新方程建立了VTI介质SH-SH波反演方法,该方法相比VTI介质的PP波反演不确定性明显下降,同时降低了常规PP波各向异性反演对大角度数据的要求.为了获得独立的横波各向异性参数,进一步地提出了一种基于岩石物理关系的横波各向异性参数估算方法.合成地震数据测试和柴达木盆地九分量地震勘探实际地震数据应用结果表明,新方法能够准确地预测地层的横波阻抗、水平横波(SH波)速度、各向异性参数,为各向异性地层的岩性解释和油气储层预测提供了可靠的解决方案,从而深化了横波地震勘探的应用潜力.     

双相各向异性介质中弹性波传播特征研究

随着地震工程和能源地震勘探的深入发展,人们所遇到的地下介质愈来愈复杂．常规的各向异性介质理论或双相各向同性介质理论难以精确描述含流体的各向异性介质,如裂缝性气藏、含水页岩等．本文以Biot双相各向异性介质理论为基础,利用弹性平面波方程,推导出了任意双相各向异性介质中弹性波的Christoffel方程．根据Christoffel方程,计算并分析了频率对双相横向各向同性介质中弹性波的相速度、衰减、双相振幅比和偏振特征的影响．结果表明,在4类波(快纵波、慢纵波、快横波和慢横波)中,频率对慢纵波影响最大;当耗散很大时,快纵波、快横波和慢横波的流固相振幅比值近似为1．对偏振特征分析的结果表明,在双相各向异性介质中,弹性波的固相位移偏振方向与流相位移偏振方向将不再保持同向或反向,而是呈不同大小的夹角．      

解耦纵横波反射波走时反演

基于常规弹性波动方程的反射波走时反演结合走时和反射波信息可以有效的摄取模型参数中的低波数成分,然而纵横波之间的耦合效应以及纵横波速度对波场的敏感性差异,导致反演的非线性问题增强.为此本文研究了基于解耦波动方程的反射波走时反演,并提出改进的时移互相关目标函数,分别隐式计入射波场快照与反传波场快照的时移量,很大程度的降低了纵波、横波之间的耦合关系,并提高纵横波速度低波数信息的反演质量.最后模型测试证明了本文方法的正确性.     

一种带误差校正的横波速度预测方法

本文提出一种基于纵横波速度关系的带误差校正的横波速度预测方法,它主要针对在进行横波速度预测时所采用纵波模量M近似变换结果与标准的Gassmann关系(K变换)之间存在的误差提出了校正的方法,并给出了完整的基于纵横波速度关系的带误差校正横波预测技术流程,该技术是一个简单实用的横波速度预测工具.     

VTI介质解耦合声波方程数值模拟

Alkhalifah提出的声学近似方程对于各向异性VTI介质地震资料分析处理是一种较为高效的方法,但由于算法本身所产生的伪横波干扰,以及所带来的稳定性问题制约了该技术在实际生产中的应用.本文介绍了纵横波分离的VTI介质解耦合声波近似方程,以解决伪横波干扰问题.并针对数值求解问题,研究了伪解析法在VTI介质声波近似解耦方程数值模拟中的应用.对比数值模拟结果,该方法准确描述了VTI介质qP波波场特征的同时,消除了伪SV波干扰,使得波场传播更稳定,数值频散也更小.     

基于纵横波解耦的三维弹性波逆时偏移

纵横波波场分离是弹性波偏移方法的必要条件,通过纵横波成像的差异可以获取更多地下介质的信息.目前所用的纵横波波场分离方法多采用Helmholtz分解,这样得到的波场不仅物理意义发生了变化,振幅和相位也会发生改变.本文采用纵横波解耦的弹性波方程,将其应用于三维介质,对比分析了纵横波解耦方法相对传统Helmholtz分解方法在相位、振幅上的优势.将该解耦的波场分离方法应用于弹性波逆时偏移,能得到相位、振幅和物理意义不受改变的偏移结果.但是该解耦方法分离得到的纵横波波场均为矢量场,将该波场分离方法用于弹性波逆时偏移,还需要解决矢量场如何得到标量成像结果的问题.本文引入了Poynting矢量,通过Poynting矢量对矢量波场进行标量化,这样就能得到保振幅、相位,且无极性反转的标量PP和PS成像结果.同时针对S波Poynting矢量求取不准确的问题,采用拟S波应力场和S波速度场得到了更加准确的S波Poynting矢量.理论计算证明了本文采用的3D波场解耦的矢量波场分离方法的正确性和引入Poynting矢量对矢量波场进行标量成像的有效性.     

正交各向异性介质中qP波入射的二阶近似反射系数与透射系数

地震波在各向异性介质中以一个准P波(qP)和两个准S波(qS1和qS2)的形式传播.研究三种波的相速度、群速度以及偏振方向等传播性质能够为各向异性介质中的正反演问题提供有效支撑.具有比横向各向同性(TI)介质更一般对称性的正交各向异性介质通常需要9个独立参数对其进行描述,这使得对传播特征的计算更为复杂.当两个准S波速度相近时具有耦合性,从而令慢度的计算产生奇异性.因此,奇异点(慢度面的鞍点和交叉点)附近的反射与透射(R/T)系数的求解不稳定,会导致波场振幅不准确.本文首次通过结合耦合S波射线理论和基于迭代的各向异性相速度与偏振矢量的高阶近似解,得到了适用于正交各向异性介质以qP波入射所产生的二阶R/T系数的计算方法.与基于一阶近似的结果相比,基于二阶近似的方法提高了qP波R/T系数的精度,能得到一阶耦合近似无法表达的准确的qP-qS转换波的R/T系数解,且方法适用于较强的各向异性介质.     

不同尺度裂缝对弹性波速度和各向异性影响的实验研究

裂缝广泛分布于地球介质中并且具有多尺度的特点,裂缝尺度对于油气勘探和开发有着重要的意义.本文制作了一组含不同长度裂缝的人工岩样,其中三块含裂缝岩样中的裂缝直径分别为2 mm、3 mm和4 mm,裂缝的厚度都约为0.06 mm,裂缝密度大致相同(分别为4.8%、4.86%和4.86%).在岩样含水的条件下测试不同方向上的纵横波速度,实验结果表明,虽然三块裂缝岩样中的裂缝密度大致相同,但是含不同直径裂缝岩样的纵横波速度存在差异.在各个方向上,含数量众多的小尺度裂缝的岩样中纵横波速度都明显低于含少量的大尺度裂缝的岩样中纵横波速度.尤其是对纵波速度和SV波速度,在不同尺度裂缝岩样中的差异更明显.在含数量多的小尺度裂缝的岩样中纵波各向异性和横波各向异性最高,而含少量的大尺度的裂缝的岩样中的纵波各向异性和横波各向异性较低.实验测量结果与Hudson理论模型预测结果进行了对比分析,结果发现Hudson理论考虑到了裂缝尺度对纵波速度和纵波各向异性的影响,但是忽略了其对横波速度和横波各向异性的影响.     

基于Lowrank近似的TTI介质一步法qP波正演模拟

近年来,基于Lowrank近似的qP波正演模拟受到广泛关注。通常的Lowrank近似两步法波场外推虽然也能消除伪横波干扰,但是两步法依赖一个实值的相位函数,且在大时间步长不适用。基于Lowrank近似的qP波一步法波场外推方法具有处理复值相位函数的能力,对相位函数进行改进,加入速度梯度项,并成功应用于TTI介质qP波正演模拟。通过试验可知:基于Lowrank近似的一步法波场延拓具有传统两步法的优点,无伪横波干扰,在大时间步长波传播时比两步法更加稳定;各向异性情况下,一步法波场外推适用于任意倾角情况下,波场模拟结果清晰准确,没有发生不稳定。      

部分叠加资料纵横波速度测井约束反演

本文提出一种角度部分叠加资料同时反演纵波和横波速度的方法,以角度部分叠加资料为基础,利用地震波振幅随入射角变化与弹性参数间的数学关系,基于非线性最优化理论,通过将最小平方问题转化为大型带状矩阵的求解问题,采用逐道外推技术依次求得该角度剖面的每个点纵波和横波速度值,利用地震波速度在各向异性介质中与入射角的关系,进一步得到地震波零入射角纵、横波速度,从而可以求得纵横波速度比及泊松比,为地震资料的岩性及含油气性解释提供了丰富的参数信息,为保证反演的稳定性,降低多解性,采用测井、叠后构造层位解释结果的共同约束,并做好叠前去噪以及振幅保持处理.     

VTI介质角度域转换波高斯束偏移成像方法研究

转换波偏移可以利用纵横波波场信息,得到高分辨率的成像结果,从而为油藏描述提供高质量的地震资料.目前的研究主要是利用纵波波场信息进行偏移成像,然而,传统的纵波方法在复杂探区成像时具有一定的局限性.为此,本文在各向异性介质声波射线追踪算法的基础上,推导出各向异性介质转换波射线追踪方程,发展了一种转换波射线追踪算法;并将研究的追踪算法应用到偏移成像中,提出了一种各向异性VTI介质角度域转换波高斯束偏移成像方法.通过各向异性VTI介质断块模型和复杂构造模型试算,说明了本文方法的正确性和有效性.模型试算的结果表明,在考虑地下各向异性时,本文研究的方法具有更好的成像效果,提取的角道集结果可以为偏移速度分析提供依据.     

基于散度和旋度纵横波分离方法的改进

纵、横波的分离是多波多分量地震资料处理中很重要的一步,其分离结果直接影响到后续数据处理的质量.各向同性介质中纵波为无旋场,横波为无散场,因此可以在频率-波数域利用散度和旋度算子对地震记录进行纵、横波分离,但是此处理过程必须知道地表处的纵、横波速度.本文给出了一种估算地表纵、横波速度的方法,可以在纵、横波速度值未知的情况下,将其估算出来.针对弹性波场进行散度和旋度运算时,纵、横波的相位和振幅比发生改变的问题,本文给出了相位和纵、横波振幅比的校正方法.     

各向异性衰减薄层地震响应特征研究

薄储层的叠前地震响应特征研究,特别是针对具有速度各向异性的含流体薄层,对储层描述具有十分重要的意义.文中基于波动方程数值模拟方法,正演得到各向同性弹性、各向同性衰减、速度各向异性、各向异性衰减模型的地震波场,并对比分析了四种模型的纵波(PP)和转换横波(PS)地震反射特征.研究结果表明:在衰减介质背景下,引入各向异性,PP和PS波的反射波振幅较弹性介质均减弱,且衰减因素对薄层振幅的影响强于各向异性.同时,VTI各向异性衰减在单频PS振幅曲线上表现出强差异性,而HTI各向异性衰减则会影响PP和PS波单频振幅曲线的极值点幅值和位置,通过分析单频振幅曲线的极值点振幅和极值点位置对各向异性衰减薄层预测有指导作用,尤其对平时较难分辨的VTI各向异性衰减薄层,单频分析方法的优势更明显.     

库车坳陷高陡构造地区零井源距VSP资料波场分析与识别

库车坳陷属于高陡构造地区,油气资源丰富,勘探潜力巨大.但是该地区地下构造复杂,断裂发育,地层倾角陡,地震波场复杂,地震成像品质较差是制约库车复杂构造研究的主要因素.零井源距VSP资料波场上转换横波很发育,波场记录上转换横波和纵波同时存在造成该地区地面地震波场相当复杂,影响地面地震资料的成像品质.因此为了避免地面地震资料处理过程中各种波场的相互影响,利用VSP波场信息能防止地面地震资料处理过程中可能出现的"陷阱".本文充分发挥VSP波场信息丰富的特点,通过库车坳陷高陡构造地区零井源距VSP波场正演模拟和实际资料波场分析,准确识别纵波和转换横波波场,并利用有效的波场分离手段分离出VSP纵波波场和VSP转换横波波场.通过VSP纵横波波场与实际地面地震资料对比,能很好地识别出了地面地震资料上的纵波和转换横波波场,为库车坳陷高陡构造地区地面地震资料处理提供参考依据.研究表明:库车坳陷高陡构造地区转换横波在地震纵波成像中影响较大,建议库车坳陷高陡构造地区地面地震处理要加强转换横波的识别与应用,进一步加强井地联合处理攻关研究,获得更可靠、更高品质的纵波资料,同时也能得到更为丰富的转换横波资料.     

转换波在储层预测中的研究分析

岩石的各向异性与油气的运移和储集有着密切关系.在油气勘探和油气田开发中,横波对各向异性的敏感性具有重要价值,它可提供一些其他方法无法获取的新信息.纵波在地下传播时,当通过路径的岩石存在各向异性时,会在波阻抗界面产生纵波和横波.这相当于纵波震源同时激发出纵波和横波,利用这些信息就可以对地下岩石的物性进行研究.本文以有限差分为基础,以地下油气田的储藏地质特性为对象,利用弹性波波动方程的传播特性, 研究P_SV转换波,分析其传播特征,用来指导其在储层预测中的应用.通过研究分析可以看出,P_SV转换波通过含油气介质时受影响比较少,能够得到比纵波好的多的成像资料.