稀疏脉冲法
包括最大似然反褶积、L1范数反褶积、最小熵反褶积、最大熵反褶积、同态反褶积等，稀疏脉冲反演是基于脉冲反褶积基础上的递推反演方法,其基本假设是地层的强反射系数是稀疏分布的。从地震道中根据稀疏的原则提取反射系数,与子波褶积后生成合成地震记录;利用合成地震记录与原始地震道残差的大小修改参与褶积的反射系数个数,再作合成地震记录;如此迭代,最终得到一个能最佳逼近原始地震道的反射系数序列。该方法适用于井数较少的地区,其主要优点是能够获得宽频带的反射系数,较好地解决地震反演的多解性问题,从而使反演结果更趋于真实。
约束稀疏脉冲反演采用一个快速的趋势约束脉冲反演算法,用解释层位和井约束控制波阻抗的趋势和幅值范围,脉冲算法产生了宽带结果,恢复了缺失的低频和高频成分;同时,再加入根据井的波阻抗的趋势约束。约束稀疏脉冲反演最小误差函数是J=∑(ri)p+λq∑(di-si)q++α2∑(ti-Zi)2(1)
式中:ri为样点的反射系数;zi为样点的波阻抗;di是原始地震道;si是合成地震道;Zi介于井约束的最大和最小波阻抗之间;ti是用户提供的波阻抗趋势;α为趋势最小匹配加权因子;p,q为L模因子;i是地震道样点序号;λ为数据不匹配加权因子。
如果从最大似然反褶积中求反射系数r(t),则在上述过程中为了得到可靠的反射系数估计值,可以单独输入波阻抗信息作为约束条件,从而求得最合理的波阻抗模型
Z(t)=Z(t-1)（1+r(t)）/（1-r(t)） (2)
稀疏脉冲法假设反射系数是稀疏的、离散的，利用测井资料可以得到井旁道的准确反射系数，通过上述反褶积方法，在测井资料、地质模型的约束下，逐道递推子波、反射系数，从而反演出波阻抗、速度等数据。
常规递推法与稀疏脉冲反演法主要是利用反褶积方法来恢复反射系数序列，由经过标定的反射系数序列递推出相对波阻抗，然后加上从声波测井和地质模型中得到的低频分量，最终得到反演波阻抗。这两类方法的主要缺陷是选择可靠低频信息较为困难，由反射系数递推波阻抗过程中误差积累快，当反射系数存在较大误差时，递推出来的波阻抗剖面会面貌全非。
2012_WPI_23_时间空间域全波形反演 精品资料
仅供学习与交流，如有侵权请联系网站删除 谢谢2 时间-空间域全波形反演
任浩然，王华忠
同济大学海洋与地球科学学院波现象与反演成像研究组，上海，200092
摘 要 地震波反演成像的理论研究地震波场在地球介质中的传播现象，从地震数据直接反求物性参数。而全波形反演是利用全波场进行地震反演的方法。在反演框架中，正演方法是反演成败的基础。基于二次型泛函，非线性寻优可以采用梯度导引类的解法，也可以采用牛顿类的解法。我们主要讨论梯度导引类方法求解时间-空间域全波形反演的方案。基于模型试验，讨论了全波形反演方法的优缺点，并就全波形反演的实用化提出了相应策略。
[关键词]：地震波反演，全波形反演，特征波场，非线性反演
1 引言
从Tarantola(1984)建立以波动方程为基础的地震反演方法以后的近三十年来，地震波反演方法已经建立了较为完善的理论体系。这些方法可以分别从不同的模型参数化、不同的正演方法和不同的反演策略来区分。其中，介质模型的参数化与正演方法相互对应。对模型有什么样的认识，才能建立什么样的正演方程,有什么样的正演方程才会有什么样的反演策略。
逆散射层析技术基于波场的线性近似，将地下的成像点看作是互不相干的散射点。传统上的散射层析只针对均匀背景模型（Devaney，1982；Wu和
Toksöz，1987）。后来发展的广义散射层析基于De Wolf近似(De Wolf，1985)，这种对散射波场进行多次前向散射和单次后向散射的近似使得广义散射层析在非均匀背景上进行反演(朱小三，2010)。基于旅行时的层析反演方法可以分为：利用初至波或折射波的近地表速度反演(刘玉柱，2007)，利用反射波的深层速度反演（Billette和Lambare，1998），利用井间透射波旅行时反演井旁精细速度反演（Harris等，1995）等。基于射线理论的各种反演方法本质上还是模拟波场的旅行时信息，因此可以反演影响旅行时的速度等参数。然而，勘探工业发展的现实表明，仅仅研究速度参数是不够的。后期的地震解释要求偏移成像是一个保振幅的，即按照偏移或者反演估计出的反射系数随角度的变化应当能够反应真实的地下物性。这就要求对模型的参数化要更加贴近物精品资料
- 1 - 波动方程反演问题的一种新的逼近方法
波动方程反演问题是一个既重要又有挑战性的数学问题，它在物理学，工程学，医学，经济学等多个领域有广泛的应用。其本质是从测量参数求解未知函数，这个过程称为反演。由于传统反演方法的计算复杂度和模型误差问题，迫切需要一种新的反演方法来提高反演效率和精度。
本文提出了一种新的反演方法，用于求解波动方程反演问题。我们将波动方程模型化为一个非线性最优化问题，分别使用快收敛的误差函数法和混合粒子群优化算法来求解未知参数。在实验室数据的基础上，我们将反演方法应用于几种典型的波动方程模型，并进行了详细的分析。结果表明，与传统方法相比，本文提出的反演方法具有更高的反演准确性和更快的收敛速度。
1.动方程反演问题
波动方程反演问题是指从测量参数求解未知函数，这类问题在多个领域（如物理学，工程学，医学，经济学等）有着重要的应用。传统的反演方法包括拟牛顿迭代法，反演积分法，局部最小二乘法，最小范数法和最小距离法等，但这些方法存在模型误差问题，计算复杂度大等缺点。
为了提高计算效率和反演精度，人们提出了一种新的反演方法，即数值优化。在数值优化方法中，最常用的是最优化算法，它将反演问题转化为一个最优控制问题。最优化算法可以有效地求解复杂的优化问题，其中包括粒子群优化算法，蚁群算法，模拟退火算法，遗传 - 2 - 算法，混合粒子群优化算法等。
2.的反演方法
本文提出一种新的反演方法，通过将波动方程模型化为一个非线性最优化问题来求解未知参数。我们将波动方程模型的反演问题转换为一个约束最优化问题，将实测的数据作为目标函数，引入一种收敛快的误差函数作为最优化函数。该误差函数可以快速将优化过程聚焦于最优解，并且不会收敛于局部最优解。
同时，为了提高反演精度，我们将混合粒子群优化算法应用于解决未知参数问题。该算法结合了粒子群优化算法和基于模拟退火的优化算法，可以有效搜索全局最优解，而且算法收敛速度快、效率高、可靠性高。
2021年2月第56卷第1期
•综合研究•文章编号：1000-7210(2021)01-0164-08
基于Zoeppritz方程的叠前和叠后
混合多参数非线性地震反演
张凌远®张宏兵尚作萍®严立志® 任权®
(①河海大学地球科学与工程学院，江苏南京211100;②河海大学力学与材料学院.江苏南京211100)
张凌远，张宏兵，尚作萍，严立志，任权.基于Zoeppritz方程的叠前和叠后混合多参数非线性地霣反演.石油地球 物理勘探，2021,56 (lh 丨 <>4-171.
摘要目前通过Shuey近似式或Gray近似式及基追踪理论反演泊松比，或依据叠前反演获得的纵、横波速度计 算纵横波速度比，但不可避免地降低了参数反演精度。为此，尝试通过精确的Zoeppritz方程直接反演纵横波速 度比或泊松比。基于垂直入射反射系数和精确Zo印pritz方程的褶积模型，结合边界保护正则化建立叠前和叠 后地震反演的目标函数，以降低反演不适定性的不利影响，并使用快速模拟退火算法进行全局非线性优化。由 于多参数之间的量级差异以及多参数同步随机搜索增大了不确定性，在反演时通过在改进Fatti方程中引入密 度与纵波速度、横波速度与纵波速度两个拟合式提高反演精度与稳定性。实际资料反演结果表明：精确Zoep- prhz方程直接反演纵横波速度比的精度明显高于间接反演，各参数反演结果与测丼数据基本吻合；叠后反演获 得的纵波速度和密度精度较高，叠前反演获得的纵横波速度比的精度高于叠后反演，18°〜24°角度道集的反演 效果好于3°〜9°和33°〜39°角度道集。
关键词叠前反演混合反演 Zoeppritz方程边界保护正则化纵横波速度比气层识别
中图分类号：P631 文献标识码：A doi： 10. 13810/j. cnki. issn. 1000-7210. 2021.01. 019
〇引言
最近十几年，基于Zoeppritz方程的叠前反演可
