VSP直达波和反射波波动方程走时联合反演

地面资料全波形反演采用低波数回折波敏感核和高波数偏移等时线恢复地下模型的长波长和短波长分量.当回折波的穿透深度有限时,很容易陷入局部极值.波动方程反射波走时反演采用反射波的透射敏感核更新速度模型.当浅层速度存在较大误差时,仍无法得到满意的反演结果.VSP资料中直达波是一种透射波,穿透深度大,可用于井旁背景速度场建模.本文发展了针对VSP观测方式的波动方程走时联合反演(直达波和反射波)方法.采用全局优化参数反演法实现VSP上下行波和纵横波分离,构建了基于下行直达波和上行反射波走时残差的混合目标函数,推导了相应的伴随状态方程和梯度公式,给出了背景速度和反射系数分步反演流程.理论和实际VSP资料应用表明：直达波和反射波走时联合反演可以得到运动学特征精确的速度模型,为后续全波形反演提供可靠的初始模型;联合反演比单独直达波/反射波走时反演的精度高,比直达波和反射波走时级联反演的耗时少.     

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

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

各向异性TI介质qP反射波走时层析成像

地震走时层析成像是反演地层各向异性参数分布的有效方法,但是关于地震各向异性介质走时层析成像的研究并不多,其技术远远没有达到成熟的阶段.在野外数据采集时,地表反射波观测方式相对井间和垂直地震剖面观测方式的成本更低,利用qP反射波走时反演各向异性参数具有更加广泛的实用价值.本文实现的TI介质地震走时层析成像方法结合了TI介质反射波射线追踪算法、走时扰动方程和非线性共轭梯度算法,它可以对任意强度的TI介质模型进行反演,文中尝试利用qP反射波走时重建TI介质模型的参数图像.利用qP反射波对层状介质模型和块状异常体模型进行走时反演,由于qP波相速度对弹性模量参数和Thomsen参数的偏微分不同,所以可以分别反演弹性模量参数和Thomsen参数.数值模拟结果表明:利用qP反射波可以反演出TI介质模型的弹性模量参数与Thomsen参数,不同模型的走时迭代反演达到了较好的收敛效果,与各向同性介质走时反演结果相比较,各向异性介质走时反演结果具有较好的识别能力.     

垂向非均匀介质的走时反演

本文提出了一种利用反射波走时曲线计算垂向非均匀介质速度和反射界面深度的方法。当在地球表面没有获得来自地下某一深度范围内介质的任何信息时,可以认为这一深度范围内的地震波速度具有连续性。利用来自其底部反射界面的反射波走时曲线,可以计算出该深度范围的地震波速度结构。对三种模型进行了理论计算,所得反演解与真实值较为一致,其中计算出的反射界面深度最为精确。利用本文提出的方法可以计算两个相邻反射面之间的垂向非均匀速度结构,如低速层等。     

利用三维弯曲界面的反演方法重建唐山震区莫霍界面三维构造形态

研究利用反射波的走时反演三维弯曲界面和介质层速度的计算方法．各层界面利用分段非完全三次多项式描述．正问题采用一种快速的三维射线追踪方法．反演过程采用变阻尼最小二乘法．数值模拟结果表明，解很快收敛到真模型．处理了通过唐山震区的实测资料，重建了该震区莫霍界面三维构造形态，并揭示了该区域构造与地震活动的关系．      

三维复杂地壳结构非线性走时反演

中国大陆中西部乃至全球造山带普遍具有复杂地壳结构.随着矿产资源勘探和深部探测研究的深入,探测造山带及盆山耦合区下方地壳精细结构正逐渐成为当前面临的巨大挑战.人工源深地震测深方法正越来越清晰地揭示出不同构造域地壳速度结构的基本特征,然而传统的层状结构模型参数化方法难以准确描述复杂地质模型,通常情况下多忽略速度结构的精细间断面且采用层边界平滑处理,难以满足地壳精细结构成像的发展要求.针对上述困难,本文采用最近发展的块状结构建模方案构建三维复杂地壳模型,基于逐段迭代射线追踪正演走时计算方法,推导了走时对三角形界面深度以及网格速度的偏导数,开展了非线性共轭梯度走时反演方法研究.发展了利用直达波和反射波等多震相走时数据对界面深度和网格速度的多参数联合反演方法,并引人不同种类震相数据的权系数和不同类型参数偏导数归一化的方法.数值算例表明,基于块状结构的非线性共轭梯度走时反演方法适用于复杂地壳结构模型,在利用人工源走时数据反演复杂地壳精细结构领域具有良好的应用前景.     

高斯波包反射走时速度反演方法

扰动高斯波包理论指出,在Gabor域描述模型的扰动成分,且入射波场为短时宽带信号时,扰动波场可在时间域通过高斯波包算子描述.在此基础上通过拟合反射波的走时,提出一种速度反演方法.反射波走时残差利用地震道局部波形的互相关函数表示,以走时残差的二范数作为目标函数,优化目标函数实现对速度场的反演.基于一阶Born近似,利用扰动高斯波包理论推导出目标函数对速度场的梯度是本文理论部分的核心内容.梯度包括两部分:正传的背景波场与反传的扰动高斯波包之间的互相关,反传的背景波场和正传的扰动高斯波包之间的互相关.梯度表达式中背景波场和扰动波场均利用高斯波包算子模拟.计算梯度的具体算法中,如何模拟扰动波场,以及如何计算反射波的走时残差是两个要点,文中对此做了详细的讨论.数值实验进一步阐述了反演的实现策略,实验结果表明高斯波包反射走时速度反演方法和实现策略有效可行,并得到了理想的反演结果.     

波动方程初至波多信息联合反演方法

地震初至波中包含着丰富的近地表速度结构信息,如何分阶段、分尺度地利用这些信息进行近地表速度建模是地震勘探中的一个关键问题.在速度反演的不同阶段,综合利用初至波中的不同信息(如走时、包络和波形等)进行联合反演,可以有效地降低反演对初始模型的依赖程度,提高近地表速度模型的反演精度.为此,本文提出了一种统一基于波动方程正演引擎的初至波多信息联合反演方法,该方法同时匹配观测和模拟的初至波走时、包络和波形信息.在不同的反演阶段选择不同的权重因子调节不同信息的权重,这样不仅降低了反演对初始速度模型的依赖,而且自然地实现了多尺度反演.在每轮反演迭代中,一次正演模拟的波场同时应用于初至波走时、包络和波形匹配,无需额外的射线追踪.同时,联合反演方法在一定程度上缓解了串联反演中目标函数漂移问题,提高了近地表速度建模的精度.     

地震立体层析成像的实现方法及效果分析

立体层析成像是一种新的地震反射波层析成像方法,能为叠前深度偏移提供较为精确的宏观速度模型。本文研究了立体层析成像的实现方法,包括斜率与走时数据的拾取、离散速度模型构建和初始化、射线参数的确定、斜率和走时及射线计算以及反演问题解法等,建立立体层析成像的算法流程。并通过对Marmousi模型试验,对立体层析成像运行所需的主要参数,如初始速度模型、拾取数据量、离散网格尺寸、速度平滑权重等进行测试和分析,总结这些不同参数对立体层析反演结果的影响规律,用以指导生产实践。     

数据域特征波反射走时反演

地震波的运动学信息(走时、斜率等)通常用于宏观速度建模.针对走时反演方法,一个基本问题是走时拾取或反射时差的估计.对于成像域反演方法,可以通过成像道集的剩余深度差近似计算反射波时差.在数据域中,反射地震观测数据是有限频带信号,如果不能准确地确定子波的起跳时间,难以精确地确定反射波的到达时间.另一方面,如果缺乏关于模型的先验信息,则很难精确测量自地下同一个反射界面的观测数据同相轴和模拟数据同相轴之间的时差.针对走时定义及时差测量问题,首先从叠前地震数据的稀疏表达出发,利用特征波场分解方法,提取反射子波并估计局部平面波的入射和出射射线参数.进一步,为了实现自动和稳定的走时拾取,用震相的包络极值对应的时间定义反射波的到达时,实现了立体数据中间的自动生成.理论上讲,利用包络极值定义的走时大于真实的反射波走时,除非观测信号具有无限带宽(即delta脉冲).然而,走时反演的目的是估计中-大尺度的背景速度结构,因此走时误差导致的速度误差仍然在可以接受的误差范围内.利用局部化传播算子及特征波聚焦成像条件将特征波数据直接投影到地下虚拟反射点,提出了一种新的反射时差估计方法.既避免了周期跳跃现象以及串层等可能性,又消除了振幅因素对时差测量的影响.最后,在上述工作基础之上,提出了一种基于特征波场分解的新型全自动反射走时反演方法(CWRTI).通过对泛函梯度的线性化近似,并用全变差正则化方法提取梯度的低波数部分,实现了背景速度迭代反演.在理论上,无需长偏移距观测数据或低频信息、对初始模型依赖性低且计算效率高,可以为后续的全波形反演提供可靠的初始速度模型.理论和实际资料的测试结果证明了本文方法的有效性.     

二维复杂层状介质中地震多波走时联合反演成像

采用新近提出的多次波射线追踪正演算法,结合共轭梯度法求解带约束的阻尼最小二乘最优化反演问题,分析讨论了利用多震相走时资料进行联合反演成像的方法及技术.考虑到不同震相走时的拾取误差不同,反演算法中引入了不同震相种类数据的权系数; 由于同时反演速度模型和反射界面起伏中不同模型参数变化对走时影响程度的不同, Jacobi偏导矩阵元素中引入了不同参数的归一化因子; 另外,为了克服射线密度过大(或过小)区域速度模型的过度(或欠)更新问题,反演算法中引入了等权射线密度的概念.几种数值模拟实例表明(含噪声敏感性试验): 多波走时的联合或同时反演成像技术是一种提高走时成像空间分辨率,进而降低重建模型失真度行之有效的方法.     

远震反演中新的地壳异常体改正方法

近年来，远震走时层析成像方法有了长足的发展，在地下结构反演的研究中获得了众多成果。针对射线在台站下方覆盖率较差而导致远震反演方法对地壳速度约束不足的问题，本文提出了一种新的地壳异常体改正方法——对同一台站的相对走时残差进行求和平均去均值，消除了地壳中复杂的速度异常体对上地幔速度结构反演的干扰，反演结果具有更优的相对走时残差分布，同时反演得到的速度模型具有更小的数据方差。     

基于多震相联合反演江苏地区一维速度模型研究

利用地震走时数据,采用联合反演方法获取了江苏地区的一维P波速度模型。与仅采用初至波走时的传统天然地震走时获取方法相比,该方法充分利用了大量存在的续至波参与反演,能有效改进中下地壳的反演能力。针对地震震相目录中常存在震相标识错误的问题,采用的自动判别筛选震相方法能最大限度提高数据走时的精度,可以对不同震相进行有效区分。与其他常用一维速度模型相比,本文反演的模型对Pg、Pn震相走时拟合效果最佳,残差最小。当所用走时数据拥有较高定位精度时,该反演方法能为研究区三维速度结构成像和地震定位提供较可靠的一维速度模型。     

基于振幅加权的拖缆与海底地震资料联合波形反演方法

南黄海新近系沉积层与其下方中-古生代海相地层之间速度差异大,存在强波阻抗界面,导致接收的浅部和中深部地层的反射波振幅差异极大.一般全波形反演方法未考虑不同界面反射波振幅的差异,对南黄海浅部沉积层速度的反演效果好,而对中深部地层速度的反演效果差.为此,本文提出一种基于振幅加权的拖缆与海底地震资料联合波形反演方法.该方法对中深部弱反射波数据施加较大的权重,对浅部强反射波数据施加较小的权重,从而均衡浅、中和深部地层反射波的振幅,使得自浅至深所有地层速度均能得到较好反演效果.同时,该方法利用海底地震资料和拖缆地震记录积分的低频成分反演低波数背景速度场,利用拖缆资料的高频成分提高反演模型的分辨率,充分发挥了拖缆和海底地震资料的优势互补作用,获得了比仅使用单一资料更好的反演结果.根据地震剖面和钻井资料,建立了南黄海中部隆起的速度模型.利用该理论模型的合成地震记录对本文方法进行了测试和分析,验证了本文方法的正确性和有效性.     

考虑水层鸣震的海底多分量数据转换波走时反演方法

海底多分量地震数据在完成上-下行P/S波分离之后,基本消除了检波器端的多次波信号,但仍受到一些源端多次波的干扰.在浅水(水深在500 m以内)、软海底环境下,经海底与海面各反射一次的水层鸣震因能量较强,对反射波偏移成像与速度分析造成极大的困扰.就横波速度分析而言,含源端水层鸣震的PPS波对一次转换PS波形成主导性干扰.本文从走时反演的视角,基于声波方程建立描述PS波和PPS波的Born模拟算子,提出同时拟合这两种S波震相的反射走时反演方法及相应的梯度预条件技术.合成数据实验表明,该方法可克服源端强水层鸣震的影响,改善转换波走时拟合条件,提高横波速度反演精度.     

横向均匀介质中地震走时的正演和反演

横向均匀介质可由若干水平分层足够好地逼近,本文引入零厚度层来描述成层介质,从而建立了回折波、反射波、首波和直达波统一的时距关系以及地震走时随模型参数变化的公式。本文详细讨论了水平分层介质里地震走时曲线的特征,合成地震走时的试验表明了由地震走时,用阻尼最小二乘法可以成功地反演速度模型参数,所反演的模型含有速度梯度层,也可以含有低速层。     

基于声波测井资料的波动方程走时层析反演速度建模

地震勘探目标逐渐由浅层(<2000 m)转向中深层(2000～3500 m),高精度速度建模是实现精确地震成像的关键手段之一.目前,初至波或早至波速度层析建模方法已较为成熟,然而对深层速度建模更为重要的反射波层析建模方法依然面临走时拾取困难及计算量大等难题.针对当前中深层走时层析速度建模方法面临的难题,本文提出了一种基于声波测井资料的波动方程全波走时层析速度建模方法及策略.该方法首先通过早至波层析反演获得浅层及大尺度速度信息;其次,通过提取初次成像剖面的构造倾角信息,利用声波测井资料进行稀疏反演,进一步获得更高精度的速度模型;最后,以此为初始模型,进行全波层析反演获得最终偏移速度模型.数值结果表明,本文所提出的速度建模方法及策略不仅可以避免传统反射波层析建模方法因采用Born近似反偏移所带来的大计算量问题,且大幅提高了中深层速度建模的精度,进而能有效地提高中深层复杂构造区域的地震成像质量.     

平移初至点走时动校正反演(英文)

浅层长排列信号在动校正后出现动校拉伸,本文从走时公式的角度深入研究其产生机理,,论证了常规走时公式不能准确描述同一反射子波中所有采样点的走时,建立在子波中心点走时基础上的传统速度分析不能得到准确的动校正速度。通过理论推导重构了同一子波周期中各信号点的走时公式,进而提出基于初至点走时的动校正反演方法,其反演精度相对于传统方法得到了大幅提升。最终采用初至点走时对长排列记录进行整体搬家取得了无拉伸校平的结果,验证了所重构公式的正确性。     

龙门山断裂带中段深地震宽角反射/折射剖面研究

射线追踪方法是研究地震波在横向非均匀地壳介质中传播的重要方法.本文推导了理论走时对网格化节点速度的偏导数公式，提出了针对深地震宽角反射/折射剖面数据反演的联合迭代法，并使用该方法对横跨龙门山断裂带中段的一条深地震宽角反射/折射剖面进行了反演和解释.首先，对每一炮的观测数据进行一维反演，在此基础上插值出一个粗略的二维速度模型；然后，使用射线追踪方法计算理论走时，再根据理论走时与观测走时的拟合程度对二维模型进行调整，以获得更加接近实际的二维速度模型；最后，利用联合迭代法对观测走时进行反演，经反复迭代使所有接收点理论走时与实测走时的残差平方和最小，最终获取该剖面的二维地壳速度结构.反演结果表明:测线东段的沉积盖层明显厚于中段褶皱带和西部高原，中部褶皱带部分地区出现基岩裸露；构造转换带两侧的地层分界面近于水平层状分布，其西侧的中、下地壳内各存在一个层间速度间断面；构造转换带内存在薄厚不等的低速层，自西向东有增厚趋势.此外，龙门山断裂带的3条主断裂向下深切结晶基底，这是由于西部松潘—甘孜地块自西向东运动，受到刚性扬子地块的阻挡，沿铲式断裂向上爬升所致；而在断层上盘距地表约15 km深处出现的最大剪应力极值区，正是发生汶川M_S8.0地震的震源位置.      

三维复杂层状介质中多震相走时联合反演成像

采用新近提出的改进型不规则最短路径多次波射线追踪正演技术,结合共轭梯度法求解带约束的阻尼最小二乘最优化反演问题,讨论了三维复杂层状模型中利用多震相走时资料进行联合反演成像的技术方法.考虑到不同震相种类走时的拾取误差不同,反演算法中引入了不同震相种类数据的权系数;另外,考虑到同时反演速度模型和反射界面起伏中不同参数变化对...