共查询到18条相似文献，搜索用时 138 毫秒
1.
偏移成像是VSP数据处理中的一个重要环节,常规的VSP成像方法通常利用VSP-CDP转换或Kirchhoff偏移,均存在保幅性差及成像精度低等问题,而波动方程叠前深度偏移被认为是对地下复杂构造进行成像的精确偏移方法.任意广角波动方程作为一种高精度的空间域单程波波动方程,同时由于只含有二阶偏导数项,易于数值实现,与其他单程波波动方程相比,具有更大的成像倾角,因此是偏移成像的有力工具之一.本文将AWWE推广应用到VSP数据成像中,实现了VSP时空域高角度单程波方程偏移.首先从三维标量任意广角波动方程出发,推导了完全匹配层吸收边界条件,在基本不增加计算量的前提下有效地压制了边界反射成像噪音,同时利用非线性反演算法优选参考速度来提高平方根算子的近似程度,从而提高高角度地层的成像精度.模型数值模拟实验验证了该方法的有效性,同时表明该方法在陡倾角构造情况下能取得很好的成像效果.最后对某地区实际观测的VSP资料进行了偏移成像,并与地面地震偏移结果进行了对比,显示出VSP波动方程偏移在成像分辨率上的优势. 相似文献
2.
基于任意广角波动方程的频率-空间域深度偏移方法研究 总被引：1，自引：1，他引：0
孙歧峰 《CT理论与应用研究》2010,19(2):17-24
传统的单程波动方程偏移算法对大倾角成像困难,本文基于时空域任意广角单程声波方程,通过对时间变量进行傅立叶变换得到其频率-空间域形式,利用有限差分进行离散化,设计并实现了频率-空间域有限差分叠后偏移成像算法。模型试算表明,该方法能够通过参数优化使得低阶数算子适应较大倾角地层的偏移成像。 相似文献
3.
常规的单程波波动方程偏移成像方法对大角度的高陡构造偏移成像存在内在的限制.根据波动方程在各个空间方向的数学特性和高陡构造反射地震波的传播特征,通过把地震波分解为垂向的上下行波、水平方向的前后行波和左右行波,提出基于波场垂向外推和水平方向外推相结合的单程波波动方程高陡构造偏移成像方法,即用波场垂向外推的单程波波动方程偏移成像方法解决中低角度平缓构造的偏移成像,用波场水平方向外推的单程波波动方程偏移成像方法解决中高角度陡倾构造的偏移成像.这种基于波场垂向和水平方向外推相结合的高陡构造偏移成像方法是常规单程波波动方程叠前深度偏移成像方法的补充和改进,它相对基于全波方程的逆时偏移具有计算效率上的优势. 相似文献
4.
Steeply-dipping structures migration based on the wavefield vertical- and horizontal-extrapolation 
下载免费PDF全文
下载免费PDF全文 常规的单程波波动方程偏移成像方法对大角度的高陡构造偏移成像存在内在的限制.根据波动方程在各个空间方向的数学特性和高陡构造反射地震波的传播特征,通过把地震波分解为垂向的上下行波、水平方向的前后行波和左右行波,提出基于波场垂向外推和水平方向外推相结合的单程波波动方程高陡构造偏移成像方法,即用波场垂向外推的单程波波动方程偏移成像方法解决中低角度平缓构造的偏移成像,用波场水平方向外推的单程波波动方程偏移成像方法解决中高角度陡倾构造的偏移成像.这种基于波场垂向和水平方向外推相结合的高陡构造偏移成像方法是常规单程波波动方程叠前深度偏移成像方法的补充和改进,它相对基于全波方程的逆时偏移具有计算效率上的优势. 相似文献
5.
叠前逆时偏移在理论上是现行偏移方法中最为精确的一种成像方法,其实现过程中的核心步骤之一是波动方程的波场延拓,而波场延拓的本质是求解波动方程,所以精确、快速地求解波动方程对逆时偏移至关重要.本文采用一种基于时空域频散关系的有限差分方法来求解声波方程,分析其频散和稳定性,实现波场数值模拟,并将分析和模拟结果与传统有限差分法进行对比.分析结果和模型数值模拟结果都表明时空域有限差分法模拟精度更高、稳定性更好.将时空域高阶有限差分法应用到叠前逆时偏移波场延拓的方程求解中,然后再利用归一化互相关成像条件成像,理论模型数据偏移处理获得了精度更高的成像.同时,在逆时偏移波场延拓的实现中,采用自适应变长度的空间差分算子求解空间导数的有限差分策略,在不影响数值模拟和成像精度的前提下,有效地提高了计算效率. 相似文献
6.
叠前深度偏移理论及方法一直是地震数据成像中研究的热点问题.业界对单程波叠前深度偏移方法和逆时深度偏移开展了深入的研究,但对双程波方程波场深度延拓理论及成像方法的研究还鲜有报道.本文以地表记录的波场值为基础,利用单程波传播算子估计波场对深度的偏导数,为在深度域求解双程波方程提供充分的边界条件,并提出利用矩阵分解理论实现双程波方程的波场深度外推.通过对强速度变化介质中传播波场的计算,与传统的单程波偏移方法相比,本文提出的偏移方法计算的波场与常规有限差分技术计算的波场相一致,证明了本方法计算的准确性.通过对SEAM模型的成像,在相同的成像参数下,与传统的单程波偏移算法和逆时深度偏移算法方法相比,本文提出的偏移方法能够提供更少的虚假成像和更清晰的成像结果.本文所提偏移算法具有深度偏移和双程波偏移的双重特色,推动和发展了双程波叠前深度偏移的理论和实践. 相似文献
7.
8.
逆时偏移成像具有高分辨率的特点,在频率域实施能提升偏移成像效率.数值求解波动方程是逆时偏移成像的关键环节,数值算法的优劣直接决定着偏移过程的计算效率与成像结果的质量.本文在频率域近似解析离散化(NAD)方法的基础上通过优化网格差分模板系数构造改进NAD方法,保证数值离散精度的同时可缓解阻抗矩阵的病态程度,以此提高频率域正演模拟的计算效率并将该方法运用于频率域声波逆时偏移.通过分别对凹陷模型、Sigsbee2B模型以及Marmousi模型进行偏移成像且与四阶标准NAD方法和四阶普通有限差分方法进行比较,改进NAD方法在压制数值频散与提高成像效率方面的优势得到了充分体现. 相似文献
9.
复杂地表的单程波动方程地震叠前正演 总被引：4，自引：0，他引：4
作者基于数学检波器和等时叠加原理,实现了复杂地表的单程波动方程地震叠前正演模拟。该方法采用虚拟的数学检波器接收地下的反射地震信号,灵活地将接收点布置在地表的任何地方,从而满足地表起伏的要求。此外,根据等时叠加原理, 该方法采用单程波动方程进行波场延拓和成像,计算简单快速。通过复杂正断层的数值模拟,得到了高信噪比的共炮集地震记录,并采用适用于起伏地形的深度偏移方法对该共炮集地震记录进行了叠前深度偏移,较好地实现了地震波的偏移归位,从而证明了这里提出的起伏地表的单程波动方程地震叠前正演方法是正确和有效的。 相似文献
10.
地震资料分辨率降低,得不到深层介质的精确信息实际上是由于大地吸收效应的影响.同时与双程波动方程相比单程波动方程避免了多次波的干扰并且计算效率高、占用内存少.本文首先基于开尔芬粘弹性介质模型将品质因子与单程波分步傅立叶法波场延拓算子相结合,实现了粘弹性介质波场延拓,从而将单程波弹性介质波场延拓推广到了粘弹性介质.然后在定位原理,数学检波器原理以及等时叠加原理的基础之上实现了粘弹性介质非零偏移距叠前正演模拟.最后将数值模拟得到的正演记录进行弹性偏移和粘弹性偏移并进行对比分析.通过数值算例可以看出,粘弹性介质叠前正演深层的反射波能量减弱,同相轴变粗,频带变窄,主频减小,分辨率降低;粘弹性偏移不但实现了振幅的恢复,而且同时偏移剖面的垂向空间分辨率也得到了提高. 相似文献
11.
ADI plus interpolation: accurate finite-difference solution to 3D paraxial wave equation 总被引：2，自引：0，他引：2
Yanghua Wang 《Geophysical Prospecting》2001,49(5):547-556
An accurate finite-difference solution is developed for the paraxial wave equation in 3D seismic migration. The conventional alternating-direction-implicit (ADI) scheme used in migration causes errors, because the variables in the migration problem are complex-valued, not real-valued, and the imaginary part of the higher-order spatial derivatives cannot be ignored. The accuracy of the 3D paraxial extrapolator is preserved by (i) retaining these higher-order terms so that it does not produce the apparent azimuthal anisotropy in conventional migration, and (ii) filtering the non-physical evanescent waves during the downward extrapolation. The implementation of the accurate solution consists of two steps: firstly, the application of ADI to solve two tridiagonal systems sequentially, and secondly, an interpolation between the extrapolated wavefields of successive extrapolation levels. The method is computationally efficient as it uses the ADI scheme and, in addition, couples the correction for azimuthal anisotropy and the suppression of evanescent waves into a single operation, the interpolation step. 相似文献
12.
常规完全匹配吸收边界(PML)对以近掠射角入射到界面上的波以及低频波、损耗波都会产生虚假边界反射.基于递归积分的不分裂复频移PML算法,利用复频移拉伸函数,极大地改善了PML边界条件的性能,我们进一步推导出基于递归积分的不分裂复频移PML弹性波方程交错网格高阶差分法,对长条形介质模型进行数值模拟,与常规PML算法进行比较说明该算法对以掠射角入射到PML界面的波以及PML层内损耗波的吸收效果. 相似文献
13.
A staggered-grid high-order finite-difference modeling for elastic wave field in arbitrary tilt anisotropic media 总被引：1，自引：0，他引：1
Introduction The real Earth usually presents anisotropy. Therefore, it is of theoretical and practical sig- nificance for many fields as oil and gas, seismic exploration and production, earthquake prediction, detection of deep structure and so on to study on seismic wave theory, numerical simulation method and its applications in the anisotropic media (Crampin, 1981, 1984; Crampin et al, 1986; Hudson et al, 1996; Liu et al, 1997; Thomsen, 1986, 1995; TENG et al, 1992; HE and ZHANG, 1996)… 相似文献
14.
就大型近场波动的高效数值模拟而言,稳定实现高阶人工边界是一个尚未圆满解决的问题.本文针对使用多次透射公式的SH波动集中质量有限元模拟,依据GKS定理的群速度解释,进一步阐明了人工边界与内域离散格式耦合所导致高频失稳的机理,即两者支持群速度指向内域的外行高频平面谐波,波动能量自发地从人工边界进入內域,从而导致失稳,而这类谐波是由集中质量有限元离散引入的.本文提出了消除此种耦合失稳的一种方法:通过修改有限元刚度阵来改变内域离散格式,并保证修改格式的精度不低于原有格式的精度.理论分析和数值实验表明此法能稳定实现透射边界.本文研究结果具有推广应用前景. 相似文献
15.
基于精细积分法的三维弹性波数值模拟(英文) 总被引：1，自引：0，他引：1
波动方程有限差分法是地震数值模拟中的一种重要的方法,对理解和分析地震传播规律、分析地震属性和解释地震资料有着非常重要的意义。但是有限差分法由于其离散化的思想,产生了不稳定性。精细积分法在有限差分法的基础上,在时间域采用解析解的表达形式,在空间域保留任意差分格式,发展成为半解析的数值方法。本文结合并发展了以往学者的成果,推导了任意精细积分法的三维弹性波正演模拟计算公式,并对其稳定性进行了数值分析。在计算实例中,实现了精细积分法二维和三维弹性波模型的地震正演模拟,对计算结果的分析表明,精细积分法反射信号走时准确,稳定性好,弹性波场相较于声波波场,弹性波波场成分更为丰富,包含了更多波型成分(PP-和PS-反射波、透射波和绕射波),这对实际地震资料的解释和储层分析有重要的意义。实践证明,该方法可直接应用到弹性波的地质模型的数值模拟中。 相似文献
16.
《中国科学:地球科学(英文版)》2021,(6)
The perfectly matched layer(PML) boundary condition has been proven to be effective for attenuating reflections from model boundaries during wavefield simulation. As such, it has been widely used in time-domain finite-difference wavefield simulations. The conventional PML has poor performance for near grazing incident waves and low-frequency reflections. To overcome these limitations, a more complex frequency-shifted stretch(CSF) function is introduced, which is known as the CFSPML boundary condition and can be implemented in the time domain by a recursive convolution technique(CPML). When implementing the PML technique to second-order wave equations, all the existing methods involve adding auxiliary terms and rewriting the wave equations into new second-order partial differential equations that can be simulated by the finite-difference scheme, which may affect the efficiency of numerical simulation. In this paper, we propose a relatively simple and efficient approach to implement CPML for the second-order equation system, which solves the original wave equations numerically in the stretched coordinate. The spatial derivatives in the stretched coordinate are computed by adding a correction term to the regular derivatives. Once the first-order spatial derivatives are computed, we computed the second-order spatial derivatives in a similar way; therefore, we refer to the method as two-step CPML(TS-CPML). We apply the method to the second-order acoustic wave equation and a coupled second-order pseudo-acoustic TTI wave equation. Our simulations indicate that amplitudes of reflected waves are only about half of those computed with the traditional CPML method, suggesting that the proposed approach has computational advantages and therefore can be widely used for forwarding modeling and seismic imaging. 相似文献
17.
在数值模拟中,隐式有限差分具有较高的精度和稳定性.然而,传统隐式有限差分算法大多由于需要求解大型矩阵方程而存在计算效率偏低的局限性.本文针对一阶速度-应力弹性波方程,构建了一种优化隐式交错网格有限差分格式,然后将改进格式由时间-空间域转换为时间-波数域,利用二范数原理建立目标函数,再利用模拟退火法求取优化系数.通过对均匀模型以及复杂介质模型进行一阶速度-应力弹性波方程数值模拟所得单炮记录、波场快照分析表明:这种优化隐式交错网格差分算法与传统的几种显式和隐式交错网格有限差分算法相比不但降低了计算量,而且能有效的压制网格频散,使弹性波数值模拟的精度得到有效的提高. 相似文献
18.
将弹性波方程变换至Hamilton体系,构造适用于弹性波模拟的高效显式二阶辛Runge-Kutta-Nystrm(RKN)格式,运用根数理论得到此格式的阶条件方程组.通过给定系数的限定条件,得到方程的对称解.为了使时间离散误差达到极小,提出数值频率与真实频率比较,通过Taylor展开,得到关于辛系数的限定方程,求解方程组得到最小频散辛RKN格式.对比分析时间演进方程的稳定性,得到使库朗数达到极大值的限定方程,求解方程组得到最稳定辛RKN格式.发现此两种格式为同一格式.新得到的辛RKN格式不依赖于空间离散方法,为了对比的需要,选取有限差分法进行空间离散.在频散、稳定性分析中,与常见辛格式对比,从理论上分析了本文提出的格式在数值频散压制、稳定性提升等方面的优势,数值实验进一步证实了理论分析的正确性. 相似文献