基于新的惩罚因子算法的波场重构反演

全波形反演是一种高精度的反演方法,其目标函数是一个强非线性函数,易受局部极值影响,而且反演过程计算量较大.波场重构反演是近几年提出的一种改进的全波形反演理论.该反演方法通过将波动方程作为惩罚项引入到目标函数中,通过拓宽解的寻找空间减弱了局部极小值的影响,而且反演过程不需要计算伴随波场,提高了计算效率.但该反演方法一直缺少准确的惩罚因子算法,直接影响到该方法的准确度.本文将波场重构反演拓展到时间域并利用梯度法进行波场重构.频率域的惩罚因子用来加强波动方程的约束,而时间域惩罚因子表现为调节模拟波场和实际波场的权重因子.为此,我们根据约束优化理论,在波动方程准确以及重构波场与反演参数解耦的假设下,提出以波动方程为目标函数的新的惩罚因子算法.根据波形反演在应用时普遍存在的噪音干扰、子波错误和低频信息缺失的情况下,应用部分Sigsbee2A模型合成数据对本文提出的算法进行实验.数值实验结果表明：基于新的惩罚因子算法,在其他信息不准确的情况下,波场重构反演可以给出高精度的反演结果.     

时间域波场重构反演（英文）

波场重构反演是一种改进的全波形反演理论。该反演方法通过将波动方程引入目标函数中拓宽了解的寻找空间,通过重构真实波场来计算模型梯度,大大提高了计算效率的同时还减弱了局部极小值的影响。但目前该理论基本在频率域进行,而频率域反演对计算内存的需求太高,并且很难应用到实际生产中。因此,本文将波场重构反演拓展到时间域,推导了时间域波场重构的增广方程,结合模型试算结果对波场重构的模型梯度进行了修改。数值实验表明,时间域波场重构反演准确性较高并且对低频信息具有良好的重建能力。     

基于对数目标函数的跨孔雷达频域波形反演

波形反演在探地雷达领域的应用已有十余年历史,但绝大部分算例属于时间域波形反演.频率域波形反演由于能够灵活地选择迭代频率并可以使用不同类型的目标函数,因而更加多样化.本文的频率域波形反演基于时间域有限差分(FDTD)法,采用对数目标函数,可在每一次迭代过程中同时或者单独反演介电常数和电导率.文中详细推导了频率域波形反演的理论公式,给出对数目标函数下的梯度表达式,并使用离散傅氏变换(DFT)实现数据的时频变换,能够有效地减少大模型反演的内存需求.在后向残场源的时频域转换过程中,提出仅使用以当前频点为中心的一个窄带数据,可以消除高频无用信号的干扰,获得可靠的反演结果.为加速收敛,采用每迭代十次则反演频率跳跃一定频带宽度的反演策略.实验证明适当的频率跳跃能够在不降低分辨率的基础上有效地提高反演效率.通过两组不同情形下合成数据反演的分析对比,证明基于对数目标函数的波形反演结果准确可靠.最后,将该方法应用到一组实际数据,得到较好的反演结果.     

复波场的实部和虚部对同时源波形反演的影响

高精度的成像需要准确的速度场信息,波形反演被认为是目前具有最高分辨率的速度反演方法之一.计算效率是目前全波形反演需要考虑的一个主要问题.为了很好地解决计算效率的问题,本文引入了一种高效的无串扰同时源反演方法,并详细介绍了其原理与计算流程.同时,基于此多震源同时反演方法,本文拓展出实波场反演、虚波场反演及复波场反演的三种反演策略,进而分析复波场的实部和虚部对全波形反演的影响.实验表明,相比于复波场全波形反演,无论是实波场反演还是虚波场反演,反演分辨率有所降低;相比于实波场反演,虚波场反演对初始模型的依赖性较小,目标函数的非线性较弱;最后,通过使用组合反演策略,即初始阶段采用虚波场反演,中后期阶段采用复波场反演,不仅可以降低反演的非线性,而且能够保证高精度建模.     

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

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

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

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

基于多网格的频率域全波形反演(英文)

频率域全波形反演虽然克服了时间方向上的局部极小值问题,但是地下介质的复杂性使其在空间域仍然存在局部极小值缺陷。在优化梯度法基础上,本文采用预条件双共轭梯度稳定算法和多重网格方法计算反演中的波场传播和目标函数的梯度,在保证计算速度的同时,减小计算机内存的消耗。频率域波形反演和多重网格的多尺度性质有效改善问题极小值缺陷,加快反演的收敛速度。以局部非均匀的三孔模型和Marmousi模型的数值模拟结果验证了该算法的有效性。     

简化混合域全波形反演多GPU加速策略

全波形反演利用地震记录中的振幅、走时和相位等信息,通过拟合实际地震记录和计算波场来定量提取地下介质的弹性参数,进而为勘探地震成像、速度建模以及大尺度构造演化分析等提供可靠依据.但全波形反演计算量巨大,特别是应用于三维大区块叠前数据时,生产成本仍然很高.本文介绍并比较了时间域和频率域的全波形反演方法,综合两者的优点,最终采用混合域的反演算法,并且在此基础上做了进一步的简化以提高计算效率.针对全波形反演方法应用于大规模叠前数据时易陷入局部极小值的问题,我们提出对模型数据进行分割,同时在数个小模型内进行梯度搜索,然后对比各个局域的梯度,最终找出合适的全局下降方向,以克服局部极小的隐患.该方法能够充分利用GPU的硬件特性.在GPU环境下实现本文所提出的简化混合域全波形反演算法.数值计算实例体现出新方法具有良好的计算效率、反演精度和算法可扩展性.     

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

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

多主频波场的时间阻尼全波形反演

低频成分缺失和地下速度强烈变化会导致严重的周期跳现象,是地震数据全波形反演的难题.通过对地震数据加时间阻尼和时间积分降主频处理,提出了一种可有效去除周期跳现象的多主频波场时间阻尼全波形反演方法.由浅到深的速度不准确会造成波形走时失配和走时失配的累积.浅部速度的准确反演可有效地减小深部波形走时失配与周期跳现象.对地震数据施加时间阻尼得到时间阻尼数据,利用不同阻尼值的时间阻尼地震数据实现由浅到深的全波形反演.低主频波场的周期跳现象相对高主频波场的要弱.对地震波场进行不同阶的时间积分以得到不同主频的波场,把低主频波场的全波形反演结果作为高主频波场全波形反演的初始模型.应用缺失4 Hz以下频谱成分的二维盐丘模型合成数据验证所提出的全波形反演方法的正确性和有效性,数值试验结果显示多主频波场的时间阻尼全波形反演方法对缺失低频成分地震数据和地下速度强烈变化具有很好的适应性.