针对探地雷达(GPR)双参数全波形反演中电导率反演精度差、双参数存在串扰现象、反演计算量大、易陷入局部极值等问题.作者将具有多参数调节功能的L-BFGS算法引入到GPR时间域全波形反演中,它避免了对Hessian矩阵的直接存储与精确求解,减小了存储量和计算量.结合参数调节因子的选取,有效减小了同步反演时介电常数与电导率的串扰影响,在不降低介电常数反演精度的前提下,提高电导率参数的反演精度.通过在反演目标函数中加载改进全变差正则化方法,提高了反演的稳定性,使目标体边缘轮廓更加清晰.首先以简单模型为例,对比了单尺度反演与多尺度串行反演策略的优劣,说明多尺度串行反演有利于逐步搜索全局最优解;而开展参数调节因子的选取实验,说明合适的参数调节因子可以有效改善介质电导率的反演精度;测试了不同正则化的反演效果,表明改进全变差正则化能提高反演稳定性,显著降低模型重构误差.最后,分别对含噪合成数据和实测数据进行了反演测试,说明本文提出的多尺度、双参数反演具有较强的鲁棒性,能提供更丰富的信息约束,重构图像界面清晰、反演效果好.
相似文献探地雷达(GPR)时间域全波形反演计算量巨大,内存要求高,在微机上计算难度大.本文中作者基于GPU并行加速的维度提升反演策略,采用优化的共轭梯度法,避免了Hessian矩阵的计算,在普通微机上实现了时间域全波形二维GPR双参数(介电常数和电导率)快速反演.论文首先推导了二维TM波的时域有限差分法(FDTD)的交错网格离散差分格式及波场更新策略.然后,基于Lagrange乘数法,将约束问题转化为无约束最小问题,构建了共轭梯度法反演目标函数,采用Fletcher-Reeves公式与非精确线搜索Wolfe准则,确保了梯度方向修正因子及迭代步长选取的合理性.而GPU并行计算及维度提升反演策略的应用,数倍地提升了反演速度.最后,开展了3个模型的合成数据的反演实验,分别从观测方式、梯度优化及天线频率等方面,分析了这些因素对雷达全波形反演的影响,说明双参数的反演较单一的介电常数反演,能提供更丰富的信息约束,有效提高模型重建的精度.
相似文献全波形反演可以为叠前深度偏移成像提供更高精度的速度模型,但该方法具有较强的非线性,对初始速度模型的依赖性较强,尤其是在实际应用中,地质条件复杂多变,速度变化不连续,增加了反演非线性程度,常常使反演陷入局部极小值,影响反演的精度.全变差约束在图像去噪领域应用广泛,属于非光滑约束,在去噪过程中能有效的保留图像的不连续界面和边缘信息.本文提出基于Hinge损失函数的垂向全变差约束全波形反演方法,在全变差约束的基础上,利用Hinge损失函数控制模型的更新方向,并使用原-对偶混合梯度算法进行求解,给出这一优化问题的迭代格式,有效提高了对地下不连续界面的重构精度,同时也降低反演对初始速度模型的依赖程度.数值算例证明:与常规全波形反演方法相比,基于全变差约束的全波形反演方法可以有效的重构速度模型中的不连续界面,尤其对高速体边缘的重构效果更明显,但该方法对初始速度模型的依赖性仍然较强;基于Hinge损失函数的垂向全变差约束全波形反演方法降低了对初始速度模型的依赖程度,可以从一个较差的初始模型通过循环迭代的方式最终得到同样精确的速度模型,较好的重构了高速体边缘和不连续界面.
相似文献本文提出解调包络方法来重构地震记录中缺失的低频信号,同时该方法能够降低全波形反演的非线性程度;提出伴随状态震源函数反演方法来得到精确的震源函数,并推导了梯度计算公式;解调包络方法结合低通滤波技术,实现了从低频到高频的多尺度反演策略,有效缓解了全波形反演的周波跳跃问题.数值算例证明了解调包络、伴随状态震源函数反演方法和低通滤波多尺度反演策略的可行性及优越性.震源函数反演精度测试结果表明:即使观测记录在缺失低频信息的情况下,也能反演得到精确的震源函数.缺失低频测试和抗噪能力测试结果表明:即使地震数据中缺失9Hz以下的低频信号或者信噪比极低的情况下,利用反演得到的精确震源函数进行解调包络多尺度全波形反演,同样可以得到高精度的全波形反演结果.与Hilbert包络全波形反演对比结果表明:解调包络在重构低频和降低伴随震源主频方面具有一定优势.
相似文献全波形反演是一种建立高精度速度模型的有力工具,是偏移模式和层析模式的联合.然而,当初始模型较差、数据缺失低频成分和大偏移距数据缺失时,常规波形反演的层析成分更新较弱.因此,反演过程以偏移模式为主,容易导致反演快速陷入局部极小值.本文发展了基于波数域梯度场分解的多尺度波形反演方法(WGDFWI),从梯度场中分离出层析成分,在反演的初期主要依赖层析分量更新背景速度场,为常规全波形反演建立良好的初始模型.首先,基于一种高效的隐式波场分离方法,将梯度场分解为层析成分和偏移成分.然后在层析梯度上应用二维波数域滤波器,以缓解偏移成分泄露的问题,并利用多尺度反演策略,增强反演的稳定性.利用双层模型和Marmousi模型进行试算的结果表明,该方法可以有效重构背景速度模型,为常规波形反演提供良好的初始模型,有效提高反演精度.
相似文献多震源编码技术可以提高全波形反演的计算效率,但同时会引入串扰噪声使反演结果质量降低. 全变分约束可以有效地压制层内噪声并突出模型界面,其与多震源技术的结合,能在大大提高弹性波全波形反演效率的同时提高反演质量. 本文提出了一种高效的动态多震源全波形反演策略,可以在离散串扰噪声的同时保证照明的均匀性. 根据残留串扰噪声的分布特征,构建了与之匹配的基于各向异性全变分约束的弹性波全波形反演方法. 为了减少周期跳跃效应,将基于稀疏约束的低频重构算法应用于弹性波地震记录,给出了利用快速梯度投影算法求解各向异性全变分约束的全波形反演流程. 模型数据测试结果表明本文的方法不仅能有效地抑制多震源方法导致的串扰噪声,同时也能降低观测数据中的噪声对反演结果的影响.
相似文献低频成分缺失和地下速度强烈变化会导致严重的周期跳现象,是地震数据全波形反演的难题.通过对地震数据加时间阻尼和时间积分降主频处理,提出了一种可有效去除周期跳现象的多主频波场时间阻尼全波形反演方法.由浅到深的速度不准确会造成波形走时失配和走时失配的累积.浅部速度的准确反演可有效地减小深部波形走时失配与周期跳现象.对地震数据施加时间阻尼得到时间阻尼数据,利用不同阻尼值的时间阻尼地震数据实现由浅到深的全波形反演.低主频波场的周期跳现象相对高主频波场的要弱.对地震波场进行不同阶的时间积分以得到不同主频的波场,把低主频波场的全波形反演结果作为高主频波场全波形反演的初始模型.应用缺失4 Hz以下频谱成分的二维盐丘模型合成数据验证所提出的全波形反演方法的正确性和有效性,数值试验结果显示多主频波场的时间阻尼全波形反演方法对缺失低频成分地震数据和地下速度强烈变化具有很好的适应性.
相似文献