有限元素法全倾角波动方程偏移   总被引：2，自引：0，他引：2       下载免费PDF全文

王尚旭  牟永光 《地球物理学报》1989,32(3):308-318

本文试图从两个方面提高波动方程偏移的效果:1.使用全方向波动方程进行偏移,使之适用于任意倾角界面的反射;2.使用有限元素法提高微分方程数值解法的精度。  相似文献   

VSP时空域高角度单程波方程偏移及其应用研究          下载免费PDF全文

范廷恩  黄旭日  马淑芳 《地球物理学报》2016,59(9):3459-3469

偏移成像是VSP数据处理中的一个重要环节,常规的VSP成像方法通常利用VSP-CDP转换或Kirchhoff偏移,均存在保幅性差及成像精度低等问题,而波动方程叠前深度偏移被认为是对地下复杂构造进行成像的精确偏移方法.任意广角波动方程作为一种高精度的空间域单程波波动方程,同时由于只含有二阶偏导数项,易于数值实现,与其他单程波波动方程相比,具有更大的成像倾角,因此是偏移成像的有力工具之一.本文将AWWE推广应用到VSP数据成像中,实现了VSP时空域高角度单程波方程偏移.首先从三维标量任意广角波动方程出发,推导了完全匹配层吸收边界条件,在基本不增加计算量的前提下有效地压制了边界反射成像噪音,同时利用非线性反演算法优选参考速度来提高平方根算子的近似程度,从而提高高角度地层的成像精度.模型数值模拟实验验证了该方法的有效性,同时表明该方法在陡倾角构造情况下能取得很好的成像效果.最后对某地区实际观测的VSP资料进行了偏移成像,并与地面地震偏移结果进行了对比,显示出VSP波动方程偏移在成像分辨率上的优势.  相似文献   

深水崎岖海底地震数据成像方法与应用   总被引：6，自引：0，他引：6       下载免费PDF全文

常旭  刘伊克  杜向东  李绪宣 《地球物理学报》2008,51(1):228-234

本文针对深水崎岖海底地震数据深部成像困难的问题,讨论了波动方程叠前深度域保幅偏移的基本理论,利用数值模型验证了波动方程叠前深度域保幅偏移算法的保持振幅特性,给出了该方法与非保幅偏移方法对于复杂构造成像精度的对比,证明了保幅偏移方法可以提高复杂构造的成像精度.本文还将波动方程叠前深度域保幅偏移算法应用到实际资料的处理中,处理结果表明该方法能够有效消除崎岖海底对深部地层的影响.  相似文献   

波动方程偏移速度建模:一种直接反演方法   总被引：1，自引：1，他引：0       下载免费PDF全文

张江杰  张剑锋 《地球物理学报》2011,54(3):835-844

基于炮域波动方程叠前深度偏移所产出的角道集,发展了不需要多次迭代的波动方程偏移速度建模方法.文中分析了炮域波动方程偏移生成角道集的方法,给出了非均匀介质中角道集同相轴曲率与速度误差的定量关系,发展了基于同相轴曲率的速度模型直接更新算法.在初始模型较合理的情况下,应用一次炮域波动方程偏移计算,即可得到较准确的速度模型.这...  相似文献   

标量波动方程全倾角有限差分偏移          下载免费PDF全文

马在田  纪少游 《地球物理学报》1988,31(6):678-686

传统的有限差分法偏移是建立在单程波波动方程的不同阶数的近似式基础上的。因此,各阶近似式的偏移方法都存在一个倾角极限问题。克服倾角极限和提高极限倾角是八十年代以来有限差分偏移技术的研究目标。在这方面取得了显著的进展。本文从标量波动方程出发,通过函数替代,将它化为一个等价的方程组,用有限差分法解此方程组进行波场外推,实现全倾角偏移成象的目的。与国外已有的使用波动方程全式进行偏移的方法相比较,除方法原理本身不同外,该方法具有稳定性好,计算工作量较少的优点。  相似文献   

共接收点倾斜叠加波动方程偏移          下载免费PDF全文

甘永忠  熊维纲 《地球物理学报》1991,34(3):361-368

共接收点倾斜叠加波动方程偏移,本质上是一种叠前偏移方法.每给定一个斜率P,对经过叠前（动校正前）常规处理的地震记录中的各共接收点道集,沿直线t=τ+px进行倾斜叠加,就形成一个共接收点倾斜叠加剖面.对之进行波动方程偏移,该偏移剖面将代表地下真实构造.对一系列的p,我们可以得到一系列这样的偏移剖面.对它们作共接收点叠加,偏移叠加剖面的信噪比将超过水平叠加剖面.本文导出了在均匀、水平层状及非均匀介质条件下的共接收点倾斜叠加波动方程偏移算法.  相似文献   

波动方程偏移成像阴影的照明补偿   总被引：6，自引：6，他引：6       下载免费PDF全文

陈生昌  马在田  Wu Ru-Shan 《地球物理学报》2007,50(3):844-850

受地下复杂构造和地震数据采集系统的影响，使地震波对地下目标的照明出现不均匀性，地震采集系统难以有效地获取地下某些目标的反射信息，进而使数据偏移成像在这些目标体上出现成像阴影. 根据波场和Green函数的窗口Fourier框架展开，利用角度域波动方程偏移成像和波动方程照明分析，并结合波动方程反演理论，提出一种角度域波动方程偏移成像阴影照明补偿方法. 这种补偿方法能同时考虑地震数据采集系统和波场传播路径对偏移成像的影响，消除复杂构造区的偏移成像阴影，改进波动方程叠前深度偏移成像在复杂构造区的成像效果.  相似文献   

稀疏采样下陡角度构造的波动方程深度偏移成像   总被引：12，自引：1，他引：11       下载免费PDF全文

刘礼农  陈树民  张剑锋 《地球物理学报》2006,49(5):1452-1459

针对稀疏采样,特别是crossline方向稀疏采样的三维地震数据,提出了一个基于波动方程方法的陡角度构造深度偏移成像方法.文中针对一类混合域（空间和波数）的波动方程偏移方法展开研究.通过建立反假频的单程波算子和通过（形式上）填充空道重建理想采样的地震数据,有效地实现了陡倾角地层的准确成像.这一算法策略解决了波动方程偏移方法在实际应用中的一个关键问题,从而更好地发挥了波动方程偏移方法在复杂构造成像上的优势.二维理论数据（Marmousi模型）和三维实际地震资料成像结果表明本文方法是有效的.这一方法对现行的三维地震数据处理有重要的现实意义.  相似文献   

高阶方程偏移的分裂算法   总被引：23，自引：6，他引：23       下载免费PDF全文

马在田 《地球物理学报》1983,26(4):377-389

对于用高阶偏微分方程进行地震偏移,本文提出一种分裂算法。有限差分法目前只能解波动方程的二阶和三阶的近似方程的偏移问题,更高阶的近似方程还没有一种有效的解法。利用高阶方程进行偏移会提高偏移的精度和效果。本文就是为了解决这个问题所进行的探索。 本文阐述了分裂算法的原理,并给出了应用该算法解决偏移问题的实例。  相似文献   

任意广角波动方程深度偏移参数优化分析     

孙歧峰  白清云 《CT理论与应用研究》2011,20(1):29-38

强横向变速和陡倾角介质的成像问题是地震偏移成像的难点.本文在频率-空间域利用有限差分高阶分裂法求解任意广角波动方程实现了地震偏移成像.试验表明,参考速度的选取是影响成像精度的关键因素.脉冲响应测试表明通过优化参数任意广角波动方程能够以较低阶的方程获得较大的偏移角度.Marmousi模型数据叠前深度偏移试验表明,该方法能...  相似文献   

基于波场垂向和水平方向外推的高陡构造偏移          下载免费PDF全文

陈生昌  张博 《地球物理学报》2012,55(4):1300-1306

常规的单程波波动方程偏移成像方法对大角度的高陡构造偏移成像存在内在的限制.根据波动方程在各个空间方向的数学特性和高陡构造反射地震波的传播特征,通过把地震波分解为垂向的上下行波、水平方向的前后行波和左右行波,提出基于波场垂向外推和水平方向外推相结合的单程波波动方程高陡构造偏移成像方法,即用波场垂向外推的单程波波动方程偏移成像方法解决中低角度平缓构造的偏移成像,用波场水平方向外推的单程波波动方程偏移成像方法解决中高角度陡倾构造的偏移成像.这种基于波场垂向和水平方向外推相结合的高陡构造偏移成像方法是常规单程波波动方程叠前深度偏移成像方法的补充和改进,它相对基于全波方程的逆时偏移具有计算效率上的优势.  相似文献   

再论地震数据偏移成像   总被引：1，自引：1，他引：0       下载免费PDF全文

陈生昌  周华敏 《地球物理学报》2016,59(2):643-654

利用地震波正向传播方程对属于波形线性反演问题近似求解方法的地震数据偏移成像进行重新推导,得到了适合散射地震数据的散射偏移成像方法和适合反射地震数据的反射偏移成像方法.以地震波传播的散射理论为出发点,首先根据描述一次散射波正向传播的线性方程研究建立散射地震数据的偏移成像方法理论;利用高频近似对产生散射波场的地下速度扰动函数的空间变化进行近似,推导出地下反射率函数,再由散射波传播方程推导出基于反射率函数的反射波传播方程,然后根据描述一次反射波正向传播的线性方程研究建立反射地震数据的偏移成像方法理论.本文指出和修正了Claerbout偏移成像方法中的不足,提出的地震数据偏移成像方法是对当前偏移成像方法理论的完善,使反射地震数据偏移成像具有了更坚实的数学物理理论基础,得到的偏移成像结果相位正确、位置准确、分辨率提高.  相似文献   

衰减雷达波有限元偏移   总被引：27，自引：5，他引：22       下载免费PDF全文

底青云  许琨  王妙月 《地球物理学报》2000,43(2):257-263

高频雷达波在地球介质中有较强的衰减,反演中不可忽略．为此文中首先给出了含衰减项的雷达波的有限元方程及其偏移理论．用有限差分法或有限元法可正演合成雷达波资料,加入一定的扰动后用含衰减项的雷达波有限元方程做偏移,实例结果表明,考虑衰减项的偏移结果能使界面更好地归位,这为提高探地雷达地质解释的分辨率提供了可能性,为逐渐地实现符合雷达波自身动力学特点的处理系统奠定了基础．  相似文献   

井巷工程反射波超前探测逆时偏移成像   总被引：1，自引：0，他引：1       下载免费PDF全文

鲁光银  熊瑛  朱自强 《地球物理学进展》2009,24(6):2308-2315

探讨了井巷隧道工程反射波法超前地质预报成像问题.从一阶速度－应力弹性波方程出发,推导了二维各向同性介质情况下弹性波逆时传播的高阶差分格式,实现了弹性波在数值空间中的逆时延拓.构建反射波法隧道超前探测中的断层、软弱夹层等介质模型,以反射波探测的正演记录作为初始条件,并从程函方程出发,采用逆时差分格式求取介质模型网格空间中各点的直达波旅行时作为弹性波逆时偏移的成像条件,实现多波多分量资料的逆时偏移.偏移结果表明,逆时偏移能够使隧道壁接收到的波场准确归位,提高隧道反射波超前探测的资料处理的精度.  相似文献   

双复杂条件下带误差补偿的频率空间域有限差分法叠前深度偏移   总被引：4，自引：2，他引：2       下载免费PDF全文

叶月明  李振春  仝兆岐  张辉  徐秀刚  李栋 《地球物理学进展》2008,23(1):136-145

油气勘探的重点正转向复杂地表条件和复杂地质条件的区域.双复杂条件下的叠前深度偏移方法是解决复杂地表条件和复杂地质构造成像的有效手段,基于"逐步累加"的"直接下延"法和"波场上延"法都是解决复杂地表成像的有效手段.波动方程的频率空间域有限差分深度偏移对介质速度横向变化有较强的适应性,适宜于复杂构造的偏移成像.然而,频率-空间域有限差分法求解时引入的误差影响了成像的质量,我们用了带误差补偿的频率-空间域有限差分偏移改善了起伏地表条件下的频空域有限差分偏移质量,对模型和实际资料进行了试算,得到了较好的成像效果.误差补偿可以在若干个外推步长上进行,通过对比分析,我们发现,相对于基于起伏地表的傅立叶有限差分法偏移来说,该方法在改善起伏地表条件下偏移成像质量的同时,也具有较高的运算效率.  相似文献   

共炮检距道集波动方程保幅叠前深度偏移方法   总被引：9，自引：9，他引：0       下载免费PDF全文

刘定进  印兴耀 《地球物理学进展》2007,22(2):492-501

本文提出了一种基于双平方根算子的共炮检距道集波动方程保幅叠前深度偏移方法,将振幅误差补偿作为偏移的一部分与“运动学偏移”一起在偏移过程中实现.其基本内容包括：(1)从保幅的单平方根算子方程出发,推导出由双平方根算子定义的保幅单程波方程;(2)根据地震波摄动理论把速度场分裂为层内常速背景和变速扰动,分别在频率－波数域和频率－空间域求得波场深度延拓的偏移时移量及振幅校正系数,从而得到最终的DSR保幅波场延拓算子;(3)在高频假设条件下,把DSR保幅波场延拓公式中的积分运算进行稳相近似,得到保幅波场延拓的相移公式.理论分析和模型数值试验表明,该方法不但可以使散射能量聚焦、归位,提高成像精度;而且可以输出正确反映地下反射系数的振幅信息,为后续的地震属性分析（如AVO/AVA）提供更真实的地震信息.  相似文献   

Kirchhoff偏移法在探地雷达正演图像处理中的应用   总被引：6，自引：6，他引：0       下载免费PDF全文

戴前伟  冯德山  何继善 《地球物理学进展》2005,20(3):849-853

本文首先通过对波动方程的分析,得出了声波波动方程和雷达波波动方程形式一致性,从而说明了把广泛应用在地震数据处理中的偏移技术引入到GPR资料处理中的可行性;其后说明了时域有限差分法（FDTD）的原理,并用它合成了几种常见的雷达正演剖面;最后利用Kirchhoff积分偏移法对正演所得的雷达剖面进行偏移处理,通过对比偏移处理前后的雷达正演剖面,可知Kirchhoff积分偏移法能使雷达正演剖面中的反射波的归位,绕射波收敛,从而大大地提高了雷达正演剖面的分辨率,更好地指导GPR剖面的地质解释和验证偏移方法的有效性.  相似文献   

Viscoacoustic prestack reverse time migration based on the optimal time-space domain high-order finite-difference method   总被引：3，自引：1，他引：2  

Yan Zhao  Yang Liu  Zhi-Ming Ren 《应用地球物理》2014,11(1):50-62

Prestack reverse time migration （RTM） is an accurate imaging method ofsubsurface media. The viscoacoustic prestack RTM is of practical significance because itconsiders the viscosity of the subsurface media. One of the steps of RTM is solving thewave equation and extrapolating the wave field forward and backward; therefore, solvingaccurately and efficiently the wave equation affects the imaging results and the efficiencyof RTM. In this study, we use the optimal time-space domain dispersion high-order finite-difference （FD） method to solve the viscoacoustic wave equation. Dispersion analysis andnumerical simulations show that the optimal time-space domain FD method is more accurateand suppresses the numerical dispersion. We use hybrid absorbing boundary conditions tohandle the boundary reflection. We also use source-normalized cross-correlation imagingconditions for migration and apply Laplace filtering to remove the low-frequency noise.Numerical modeling suggests that the viscoacoustic wave equation RTM has higher imagingresolution than the acoustic wave equation RTM when the viscosity of the subsurface isconsidered. In addition, for the wave field extrapolation, we use the adaptive variable-lengthFD operator to calculate the spatial derivatives and improve the computational efficiencywithout compromising the accuracy of the numerical solution.  相似文献