共查询到18条相似文献，搜索用时 467 毫秒
1.
构建近地表横波速度模型是煤田多分量地震资料处理的重要环节。相较于面波多道分析法,全波形反演在构建近地表横波模型中具有更高的分辨率。然而,在基于梯度的全波形反演中,由于地震记录频带有限、波场的非均匀覆盖以及双重散射等原因导致梯度算子不随深度的增加而缩放,模型深部参数得不到明显更新。目标函数的Hessian算子包含曲率信息,可清晰预测梯度算子中的焦散现象及双重散射产生的伪影,因此,逆Hessian算子则可作为反卷积算子实现对梯度的预处理,加强对模型深部的照明能力。然而Hessian算子具有巨大维度,对其显式计算十分困难。基于此,借鉴逆散射理论的思想,给出勒夫波全波形反演目标函数的拟Hessian算子的表达式,并提出一种梯度预处理的全波形反演方法。将该方法分别应用于断层模型、凹陷模型以及起伏界面模型的重构试验,反演结果表明:与传统的共轭梯度全波形反演方法相比,基于拟Hessian算子的预处理共轭梯度方法可加快收敛速度,提升成像质量。 相似文献
2.
浅地表低速异常体(如空洞)的精确探测一直是地球物理领域具有重要价值与挑战性的研究课题,对城市灾害评估与复杂条件下煤田地震勘探的浅地表精细建模具有重要意义。面波全波形反演适合浅地表高精度建模,然而在实际应用中仍存在诸多亟待解决的难题。针对面波全波形反演方法中多参数串扰、实际数据预处理、震源子波估计等关键问题,开发了一套完整的瑞雷波多参数全波形反演方法流程,以实现对浅地表空洞的精确探测。该方法中模型的横波速度、纵波速度和密度随反演进程均被同步更新,减弱后两个参数偏离真实值对横波速度反演精度的消极影响。采用伴随状态法构建的拟海森算子对梯度进行预处理以压制地表伪影、增强波场照明,提高对小尺度异常的表征能力。通过褶积因子消除波场正演与实际数据采集的维数差,实现3D波场到2D波场的转换。采用校正滤波方法估计震源子波,并在迭代进程中进行动态估计,以减弱特定参数模型不准确的影响。同时,该方法采用多尺度反演策略,减轻由低速异常引起的目标函数非凸性,提高反演稳定性。合成数据和实际案例测试结果表明,瑞雷波多参数全波形反演方法得到的横波、纵波速度模型、密度模型具有基本一致性,其中横波速度模型准确度最高。实测数... 相似文献
3.
4.
层析成像中的波形反演 总被引：1，自引：0，他引：1
在波动方程有限差分数值模拟基础上，用波形记录残差平方并建立目标函数。对目标函数进行梯度优化，我们得出了介质速度的迭代公式。通过计算目标函数梯度我们可以对介质进行速度校正，从而获得反演解。经数值试验证明这是一种很有效的反演方法。 相似文献
5.
将OBFGS(Online-BFGS)优化算法应用到全波形反演中,OBFGS算法是一种不依赖Wolfe准则求取步长的优化算法。基于OBFGS算法的特点,提出了自适应步长策略,可以根据当前目标函数的变化情况自适应调整步长。基于自适应步长的OBFGS算法不要求在调整步长时计算正演和梯度,因此可以加快计算速率。通过对Marmousi型进行反演的精度分析表明,基于自适应步长OBFGS算法全波形反演提高了反演精度。反演效率分析结果表明,基于自适应步长的OBFGS算法的反演耗时较短,具有更快的收敛速度,提高了反演的效率。 相似文献
6.
7.
全波形反演是勘探地球物理领域兴起的核心技术之一,不但可以构建地下速度结构,也能够反演衰减参数(品质因子Q)模型,有助于识别地下介质类型和构造(如流体和煤层陷落柱),对煤和油气等自然资源的勘探和开发有重要意义。参数串扰是黏弹性全波形反演的关键难点,受速度误差影响,反演的Q模型会包含非常强的串扰噪声。针对该问题,提出了基于多目标函数的黏弹性全波形反演理论与方法,首先通过旅行时反演速度结构,再通过中心频率目标函数反演Q模型,最终使用波形差目标函数同时反演速度和Q模型。由于中心频率主要受衰减影响,因此,可有效减弱速度误差对Q反演的影响。最后,通过数值模拟验证了算法可有效地反演速度和Q模型。 相似文献
8.
《物探化探计算技术》2016,(5)
常规的全波形反演对初始模型和低频数据的依赖性较强,反演精确度经常会受到"跳周"现象的严重影响。将小波域最小平方滤波引入到全波形反演中,并利用小波变换的多尺度特性,有效地提高了反演过程的稳定性,减小反演受到"跳周"现象影响的可能性。小波域最小平方滤波具有更高的精度和更好的性能,利用该算法调整模拟波场的相位,从而改善模拟波场和观测波场之间的相位差异,并由此构造一个新的目标函数。同时利用小波变换的多尺度特性将地震数据分解为不同频带数据,实现多尺度反演。数值模拟实验结果表明,基于小波域最小平方滤波的多尺度自适应全波形反演,对初始模型和低频数据的依赖程度降低,较常规全波形反演具有更好的稳定性,受到"跳周"现象影响的可能性大大降低。 相似文献
9.
二维大地电磁正则化共轭梯度法反演算法 总被引：2，自引：0，他引：2
针对大地电磁二维反演中目标函数收敛速度慢而且解的稳定性较差等问题,提出了大地电磁数据的正则化共轭梯度法反演算法(Regularized Conjugate Gradient Algorithm,RCGA)。此算法在构建目标函数时引入正则化的思想,利用共轭梯度法求解最优化问题。在每次迭代过程中根据目标的收敛情况更新正则化因子,有效地解决了迭代时目标函数发散的问题。最后对两个典型模型进行了试算,对其收敛速度及解的稳定性进行了分析,并与传统反演方法的计算结果进行了对比,结果都表明RCGA具有收敛速度快、精度高、结果稳定等优点。 相似文献
10.
全波形反演不仅利用相位和振幅信息,还利用波形的细节变化,具有刻画模型精确细节的能力.在对稀疏矩阵直接LU分解求解的基础上,采用梯度预处理方法对声波介质速度模型进行了反射波全波形反演.采用误差反向传播算法计算目标函数梯度以及伪Hessian矩阵的对角线元素来做梯度预处理.数值模型的实验结果表明,利用有效的频率段便能反演出分辨率较高的速度结构,用低频反演出的结果作为高频反演的初始模型,减少了解的非唯一性.二维高斯光滑初始模型提供了有利的低频信息,得到较好的反演结果.伪Hessian矩阵的预处理吸收了高斯牛顿法的二次收敛优势,在不增加计算量的前提下,加快收敛速度. 相似文献
11.
频率-空间域正演模拟是频率域及Laplace-Fourier域全波形反演的基础,起伏地表条件下波形反演算法的关键是正演算法中考虑起伏地表的影响。基于带PML吸收边界的声波波动方程,在已有最优9点有限差分正演算法的基础上构建了起伏地表条件下频率-空间域正演算法。通过应用变网格技术,进一步提高算法的计算效率、降低内存开销,使得大规模起伏地表模型的频率域正反演问题成为可能。理论分析及数值测试表明：通过对近地表区域进行局部网格加密,可有效地压制由于矩形网格离散引起的角点散射;结合变网格技术可较易获得5倍以上计算效率的提高及内存占用的降低,且随着模型尺度的增加及地表起伏高程差的减小,倍数将显著增加;在细网格与粗网格交界处产生的虚假反射振幅幅值控制在原始波场的2%以内,满足地震波场正反演的需求。 相似文献
12.
13.
全波形反演可以利用叠前地震波场的运动学和动力学信息重建地下速度结构,具有揭示复杂地质背景下构造与岩性细节信息的潜力。然而,庞大的计算量和存储空间需求,限制了全波形反演的发展。在频率多尺度全波形反演中将L BFGS数值优化算法与同时激发震源技术相结合的方法来改善这一现状。首先,对Marmousi模型进行了速度反演：在计算过程中明显发现对计算机内存的占用减少,最终反演结果与实际Marmousi模型的拟合误差为0.095 9,较小;采用10个频带单炮震源正演384炮所需时间约为32 640 s,而采用同时激发震源（384炮）正演一次所需时间仅约为700 s。然后,基于高速楔形体模型进行了抗噪能力研究：原始含噪地震记录信噪比为11.147 3 dB;对反演得到的速度模型进行正演,其地震记录信噪比为22.251 8 dB。最后,基于逆冲断层模型进行了反演速度扰动能力研究,反演得到的最终模型很清晰,与具有速度扰动特性的实际模型非常接近,拟合误差仅为0.036 0。数值模拟试验结果表明：此方法反演精度高,内存开销较小,能够显著提高计算效率,并且具有良好的抗噪能力,能够反演出具有速度扰动特性的介质。 相似文献
14.
似海底反射层的速度异常是识别天然气水合物的重要标志,这里提出了一种针对天然气水合物似海底反射层的全波形反演方法。这种方法分为全局搜索与局部搜索二部份:首先使用遗传算法进行旅行时,反演得到背景速度模型;然后用其作为初始模型,使用共轭梯度算法进行全波形反演。通过对含噪数据的数值试验,算法表现出了较高的稳定性,并确定了进行全波形反演的遗传算子。将这种波形反演方法应用于我国南海北部海域的天然气水合物研究,反演得到了分辨率高于常规速度分析的似海底反射层速度结构,并识别出似海底反射层的速度异常。利用纵波速度反演的结果,计算出沉积物中游离气的含量,认为BSR下方的低速层可以解释为含至少1%游离气的薄层。并分析了研究区内甲烷气的来源,认为该区域游离气兼有生物气和热解气。 相似文献
15.
16.
全波形反演方法利用叠前地震波场的运动学和动力学信息重建地下速度结构,具有揭示复杂地质背景下构造与岩性细节信息的潜力。然而,巨大的计算量是阻碍其发展的一个瓶颈问题。为此,研究者们提出了震源编码技术来减少计算量,但是此方法在模型更新过程中会引进随机串扰噪声,降低反演结果准确性。所以,在保证计算精度的情况下,本文提出了采用逐减随机震源采样的方法来高效计算全波形反演问题。笔者将此方法应用于频率域二维黏滞声波波动方程全波形反演,开始了在频率域进行随机震源采样类方法的研究,计算过程中共使用了依次增大的8个频率段;并应用Overthrust模型来验证此类随机震源采样法的正确性。实验结果表明：基于逐减随机震源采样法的反演结果与实际Overthrust模型的拟合误差为0.065 65,而应用基于全部震源的全波形反演方法得到的反演结果与实际Overthrust模型的拟合误差为0.064 64,两者差别不大;但计算用时由740 min减少到291.2 min,即计算效率提高了2.54倍。为了更好地确定方法的有效性,将其应用于Marmousi模型进行试算。模型试算结果表明：基于逐减随机震源和基于全部震源得到的反演结果与实际Marmousi模型的拟合误差分别为0.080 12和0.078 97,相差不大;但计算用时由1 218.9 min减少到274.4 min,计算效率提高了4.44倍。综上,在保证反演精度的情况下,基于逐减随机震源采样法的频率域全波形反演方法大大减少了计算量,具有不可替代的计算优势,并且没有引进随机串扰噪声。 相似文献
17.
由于裂缝性油气藏具有突出的资源潜力和经济效益,利用地震方法对裂缝储层进行精细的定量描述逐渐成为勘探地球物理的关键任务之一。为了克服以往数据驱动类反演方法无法直接获得裂缝参数、而基于静态等效介质模型驱动的反演方法无法描述孔隙内部结构和流体信息的缺点,笔者提出一种基于动态等效介质模型的储层定量描述新方法。该方法通过频变AVO(amplitude variation with offset)理论建立目标函数并使用全局最优化算法反演裂缝参数。一维和二维模型测试证实,由于充分利用了反射系数频变响应对裂缝密度和时间尺度因子的敏感性,反演方法可以对裂缝储层实现有效描述。 相似文献
18.
利用一维正则化反演进行大地电磁测深数据拟二维反演解释 总被引：1，自引：0，他引：1
大地电磁测深的反演问题是不适定的,其反演结果不稳定,且具有非唯一性。通过在目标函数中采用正则化方法,可以使得不适定反演问题具有稳定的反演结果,并改善解的稳定性和非唯一性问题。为了提高野外大地电磁测深数据的处理效率和初步解释的精度,提出了大地电磁测深数据的一维正则化反演进行拟二维反演解释方法。这里所述的大地电磁测深一维反演解释,与以往的解释方法不同,其思路首先用Bostick反演的深度来控制层参数,使反演计算的模型参数仅存在电阻率;最后采用阻尼高斯-牛顿算法进行反演计算,并将Bostick反演结果作为反演计算的初始模型。通过模型试算,结果表明其处理速度快、解释直观,对野外大地电磁测深数据进行初步反演解释是可行的。 相似文献