应用时间外推或深度外推法进行地震波形的非线性反演

地震波形非线性反演的经典目标是要获取地球模型,以便对初始条件和给定震源最佳拟合观测地震波形。通常利用迭代法解这个问题:每次迭代都要牵涉给定震源在某介质中解波动方程,该波动方程的解在时间上是反向的,作为源有一定的误差;于是可得到两波场在空间每一点的关系。     

地震波形非线性反演的替代方法

大型非线性反演问题的一个普通例子是地震波形反演，通常解这类问题的方法涉及到求出观测与理论预测值之间某一失配函数的极小值。随着问题规模的增大，需要反演大矩阵的一些方法变得很麻烦，要执行反演并避免稳定性方面的问题就需要相当大的储存量和计算工作。因此，不需要任何大型矩阵反演的方法是很受欢迎的。目前，下降型算法正被广泛使用着。通常，在每次迭代时，由失配函数的梯度得出下降方向，并且根据这个方向或者还有以前的     

时间域地震全波形反演方法进展

地震全波形反演(Full Waveform Inversion)利用地震波场的运动学和动力学信息建立地层速度模型,能够精确刻画模型细节,成为近年来的研究热点.本文对日渐成熟的时间域全波形反演方法的研究进行了综述,着重评述了地震波场模拟方法、震源子波的处理方式、目标泛函的建立方式、优化反演方法、梯度预处理方式等关键环节的进展,指出了全波形反演方法向分频分尺度、全局搜索与局部搜索联合的优化反演、时间域正演和频率域反演的混合域反演等多种技术有机联合的方向发展.     

时间域平面波全波形反演

全波形反演同时利用地震波的振幅和相位信息,基于最小二乘思想反演介质参数,具有高分辨率、高精度等优点.但是,在地震数据的处理过程中运用全波形反演的运算成本非常高.为了解决这个问题,本文通过线性时移函数将多个炮记录压缩为若干个平面波记录,实现了一种基于平面波编码的时间域多尺度全波形反演方法.通过对Marmousi模型的反演测试以及与常规全波形反演对比分析,验证了本文方法的较高的运算效率和较强的抗噪声能力.     

反射地震数据的逐层波形反演

本文针对层状介质并结合梯度法波形反演,提出逐层波形反演的方法. 首先给出介质扰动响应的概念,并在此基础上分析了梯度法波形反演方法. 波形反演实质上是将实测地震记录和预测地震记录的波形残差信息转化为实际地质模型和预测地质模型的模型残差信息. 波形反演的优点是利用大量振幅相位信息得到高分辨率的反演结果, 其缺点是运行耗时大;当初始模型和实际模型相差较大时,迭代算法容易陷入局部极小点,这是因为目标函数和初始模型同实际模型间的差异是非线性的关系. 逐层波形反演方法是使自上而下每一层的目标函数最小,这样总的目标函数也是最小的. 利用二分法速度扫描确定每一层速度不仅提高了运算速度也避免了迭代算法陷入局部极小点的问题. 结合介质扰动响应和目标函数值变化可以更为准确迅速地确定每一层速度和该层界面位置.     

地震道非线性最优化反演

在非线性最优化理论的基础上提出了地震道反演的非线性最优化方法,推导出地震道非线性最优化反演所需的梯度向量、Hasse矩阵,并对测井约束反演、反演稳定性等问题进行了探讨.理论模型和实际资料的反演结果表明,这一方法是地震道反演的一种有效手段.     

一种时间域全波形反演中的地震子波估计方法

全波形反演具有高精度成像能力,然而由理论走向实际应用还存在很多问题.全波形反演通过匹配波形来更新模型,数据的波形与地震子波有直接的关系.本文介绍了一种针对时间域波形反演的子波估计方法,并将其应用到全波形反演中.由于在频率域子波反演可表示为一个线性优化问题,因此本文先给定一个试探子波,在时间域通过有限差分法正演得到地震波场,并在频率域与观测数据比较分析,反演得到预测子波.此外,本文还简要地介绍了如何使用伴随状态法计算全波形反演的梯度.数值实验证明,本文方法反演得到地震子波与真实子波具有很好的吻合度,并在全波形反演中取得了不错的效果.相比不依赖于子波的方法,提前进行子波估计可提高反演效率.此外,本文方法适用于全波形反演中常使用的多尺度反演策略.     

时移地震资料贝叶斯AVO波形反演

针对时移地震差异数据,给出了一种基于贝叶斯理论的AVO波形反演方法.该方法可以利用时移地震差异数据同时反演出纵波阻抗、横波阻抗和密度的变化.利用时移地震资料进行反演,由于采集和处理过程中存在一定的差异,不同年份地震资料在非注采过程影响区域也会存在一定的变化,而该变化会导致反演结果在非注采区域有较大的变化.针对这一问题,本文采用贝叶斯理论框架,将待求的纵横波阻抗、密度变化的先验信息和包含在地震数据中的信息结合起来,对于纵横波阻抗和密度变化,假设其服从Gauss分布,并以时移地震分别反演的结果作为其期望,同时,为了更好地表征储层属性变化,提高分辨率和抑制非注采区域弹性参数的变化,假设弹性参数变化的导数服从改进的Cauchy分布.数值模拟试验和实际资料处理结果皆表明,本文提出的反演方法能够有效地抑制假象,突出储层性质的变化,得到高分辨率的弹性参数变化信息,为研究储层属性的变化和优化开采方案提供更多的有效的信息.     

小波变换多尺度地震波形反演

完全非线性地震波形反演问题是石油地球物理勘探领域中一个非常重要而又难度很大的问题。本文提出了多尺度地震波形反演的小波变换方法,对于一维非线性地震波形反演问题,将该方法和已有的简单迭代法及多重网格法作了比较,数值实验结果表明,本方法效果较好。     

时移地震资料贝叶斯AVO波形反演

针对时移地震差异数据,给出了一种基于贝叶斯理论的AVO波形反演方法.该方法可以利用时移地震差异数据同时反演出纵波阻抗、横波阻抗和密度的变化.利用时移地震资料进行反演,由于采集和处理过程中存在一定的差异,不同年份地震资料在非注采过程影响区域也会存在一定的变化,而该变化会导致反演结果在非注采区域有较大的变化.针对这一问题,本文采用贝叶斯理论框架,将待求的纵横波阻抗、密度变化的先验信息和包含在地震数据中的信息结合起来,对于纵横波阻抗和密度变化,假设其服从Gauss分布,并以时移地震分别反演的结果作为其期望,同时,为了更好地表征储层属性变化,提高分辨率和抑制非注采区域弹性参数的变化,假设弹性参数变化的导数服从改进的Cauchy分布.数值模拟试验和实际资料处理结果皆表明,本文提出的反演方法能够有效地抑制假象,突出储层性质的变化,得到高分辨率的弹性参数变化信息,为研究储层属性的变化和优化开采方案提供更多的有效的信息.     

地震全波形反演方法研究综述

近年来,随着计算机硬件水平的提高,地震全波形反演技术研究快速发展,并有效的推进了油气勘探.本综述主要对当前地震全波形反演主要存在反演非唯一性、噪声敏感性、初始模型强依赖、易陷入局部极值、计算量大等问题进行调研,重点介绍了地震全波形反演方法在时间域,频率域和Laplace域内的各种改进和优化的策略,为全波形反演方法研究提供参考.     

小波变换多尺度地震波形反演

完全非线性地震波形反演问题是石油地球物理勘探领域中一个非常重要而又难度很大的问题。本文提出了多尺度地震波形反演的小波变换方法,对于一维非线性地震波形反演问题,将该方法和已有的简单迭代法及多重网格法作了比较,数值实验结果表明,本方法效果较好。     

地震波形多尺度反演的一点讨论

基于在反演过程中对初始模型依赖性强、易陷入局部极值等问题,本文引入小波分析,提出多尺度地震波形反演方法,从而将参数反演问题转化到小波域中重要系数优化问题。利用多尺度之间的内在联系及小波域中重要系数的稀疏性,有将改进了局部极值、计算量等问题。并对几种多尺度反演策略进行了比较讨论。基于波动方程正演及褶积模型的两种反演方法的数值实便结果显示了本方法良好的效果。     

高铁地震信号直达波波形反演建模

高铁地震信号由高速行驶的高铁列车与铁轨挤压形变后产生,具有频带宽、低频强、重复性好等优点.如何利用高铁地震信号进行地下介质重构,尚处于前期探索阶段.全波形反演建模技术是目前地球物理领域建模精度较高的方法之一,但其需要低频数据才能保证结果的稳定收敛,因此利用全波形反演对高铁地震信号进行近地表建模及属性变化监测和灾害预警具有独特优势.本文通过对高铁地震信号的分析,研究了高铁地震的理论震源信号和基于桥墩的离散震源信号,并通过模型测试实现了利用高铁地震信号的全波形反演建模及属性变化监测.

     

地震波形反演的稀疏约束正则化方法

本文考虑地震波形反演问题.为了克服传统的Tikhonov正则化方法过度光滑的弊端,引入了非线性稀疏约束正则化方法,并采用对偶方法求解稀疏约束泛函的极小点.基于二维声波方程波形反演问题进行了数值模拟,针对不同模型对稀疏约束正则化方法进行了测试.结果表明,稀疏约束正则化方法对不连续介质模型的介质边缘具有良好的识别能力.     

测定断层面解的地震波形反演方法

测定地震的断层面解,是地震学研究的重要任务。近年来的研究结果表明,用综合地震图模拟远场体波波形是测定地震断层面解的有效方法。远场体波地球介质响应近似地为一单位函数,因此远场体波综合地震图计算十分方便。该方法的缺点是只适用于大地震(震级M>6.5)。由于地震震源机制的影响,走滑断层的地震在远场(震中距离     

时间二阶积分波场的全波形反演

通过对波场的时间二阶积分运算以增强地震数据中的低频成分,提出了一种可有效减小对初始速度模型依赖性的地震数据全波形反演方法-时间二阶积分波场的全波形反演方法.根据散射理论中的散射波场传播方程,推导出时间二阶积分散射波场的传播方程,再利用一阶Born近似对时间二阶积分散射波场传播方程进行线性化.在时间二阶积分散射波场传播方程的基础上,利用散射波场反演地下散射源分布,再利用波场模拟的方法构建地下入射波场,然后根据时间二阶积分散射波场线性传播方程中散射波场与入射波场、速度扰动间的线性关系,应用类似偏移成像的公式得到速度扰动的估计,以此建立时间二阶积分波场的全波形迭代反演方法.最后把时间二阶积分波场的全波形反演结果作为常规全波形反演的初始模型可有效地减小地震波场全波形反演对初始模型的依赖性.应用于Marmousi模型的全频带合成数据和缺失4 Hz以下频谱成分的缺低频合成数据验证所提出的全波形反演方法的正确性和有效性,数值试验显示缺失4 Hz以下频谱成分数据的反演结果与全频带数据的反演结果没有明显差异.     

远震体波波形反演获得的震源时间函数

地中海东部大地震成核问题的一个有趣的现象，是它们的震源时间函数的急始特性，使用同一反演方法确定了发生在世界其他一些地区破坏了地表的浅源大地震的震源时间函数都呈现出类似的特性，这些地震的断层滑动类型包括：正断层、走滑断层和冲断层型。这些急始型震源时间函数可能反映了破坏障碍体的困难性，即在破裂成核的动力断裂过程完成后，断层的局部地区仍未破坏，也可能反映了通常发生在大地震破裂期间的某些特殊的物理过程（如蠕变不稳定性）。发生在类似地质构造环境中，且震源机制也类似的地震，反映出共同的模式。     

利用地震波形反演研究震源破裂时空过程

本文对近10多年来利用地震波形反演推断震源破裂时空过程这一领域的研究进展作一简单介绍．论述了理论地震图的计算、反演确定震源参数、以及解释反演结果三方面的发展,并对反演结果进行了归纳总结．研究表明,在本文所提到的地震中,多数大地震的震中附近只有相对小的滑动,而在多数中等地震的震中或其附近有较大的滑动．本文认为,在进一步提高反演解答的稳定性和可信程度的基础上,在大地构造背景下研究震源的性质,对一批不同震级、发生在不同构造区的地震进行系统研究是有意义的．     

利用地震波形反演研究震源破裂时空过程

本文对近10多年来利用地震波形反演推断震源破裂时空过程这一领域的研究进展作一简单介绍．论述了理论地震图的计算、反演确定震源参数、以及解释反演结果三方面的发展,并对反演结果进行了归纳总结．研究表明,在本文所提到的地震中,多数大地震的震中附近只有相对小的滑动,而在多数中等地震的震中或其附近有较大的滑动．本文认为,在进一步提高反演解答的稳定性和可信程度的基础上,在大地构造背景下研究震源的性质,对一批不同震级、发生在不同构造区的地震进行系统研究是有意义的．