共查询到20条相似文献，搜索用时 31 毫秒
1.
It is important to include the viscous effect in seismic numerical modelling and seismic migration due to the ubiquitous viscosity in an actual subsurface medium. Prestack reverse‐time migration (RTM) is currently one of the most accurate methods for seismic imaging. One of the key steps of RTM is wavefield forward and backward extrapolation and how to solve the wave equation fast and accurately is the essence of this process. In this paper, we apply the time‐space domain dispersion‐relation‐based finite‐difference (FD) method for visco‐acoustic wave numerical modelling. Dispersion analysis and numerical modelling results demonstrate that the time‐space domain FD method has great accuracy and can effectively suppress numerical dispersion. Also, we use the time‐space domain FD method to solve the visco‐acoustic wave equation in wavefield extrapolation of RTM and apply the source‐normalized cross‐correlation imaging condition in migration. Improved imaging has been obtained in both synthetic and real data tests. The migration result of the visco‐acoustic wave RTM is clearer and more accurate than that of acoustic wave RTM. In addition, in the process of wavefield forward and backward extrapolation, we adopt adaptive variable‐length spatial operators to compute spatial derivatives to significantly decrease computing costs without reducing the accuracy of the numerical solution. 相似文献
2.
叠前逆时偏移在理论上是现行偏移方法中最为精确的一种成像方法,其实现过程中的核心步骤之一是波动方程的波场延拓,而波场延拓的本质是求解波动方程,所以精确、快速地求解波动方程对逆时偏移至关重要.本文采用一种基于时空域频散关系的有限差分方法来求解声波方程,分析其频散和稳定性,实现波场数值模拟,并将分析和模拟结果与传统有限差分法进行对比.分析结果和模型数值模拟结果都表明时空域有限差分法模拟精度更高、稳定性更好.将时空域高阶有限差分法应用到叠前逆时偏移波场延拓的方程求解中,然后再利用归一化互相关成像条件成像,理论模型数据偏移处理获得了精度更高的成像.同时,在逆时偏移波场延拓的实现中,采用自适应变长度的空间差分算子求解空间导数的有限差分策略,在不影响数值模拟和成像精度的前提下,有效地提高了计算效率. 相似文献
3.
与地面地震资料相比,VSP资料具有分辨率高、环境噪声小及能更好地反映井旁信息等优点.常规VSP偏移主要对上行反射波进行成像,存在照明度低、成像范围受限等问题.为了增加照明度、拓宽成像范围、提高成像精度,本文采用直达波除外的所有声波波场数据(全波),包括一次反射波、多次反射波等进行叠前逆时偏移成像.针对逆时偏移中的四个关键问题,即波场延拓、吸收边界条件、成像条件及低频噪声的压制,本文分别采用自适应变空间差分算子长度的优化有限差分方法(自适应优化有限差分方法)求解二维声波波动方程以实现高精度、高效率的波场延拓,采用混合吸收边界条件压制因计算区域有限所引起的人工边界反射,采用震源归一化零延迟互相关成像条件进行成像,采用拉普拉斯滤波方法压制逆时偏移中产生的低频噪声.本文对VSP模型数据的逆时偏移成像进行了分析,结果表明:自适应优化有限差分方法比传统有限差分方法具有更高的模拟精度与计算效率,适用于VSP逆时偏移成像;全波场VSP逆时偏移成像比上行波VSP逆时偏移的成像范围大、成像效果好;相对于反褶积成像条件,震源归一化零延迟互相关成像条件具有稳定性好、计算效率高等优点.将本文方法应用于某实际VSP资料的逆时偏移成像,进一步验证了本文方法的正确性和有效性. 相似文献
4.
在实际地球介质中传播的地震波会产生衰减和频散现象,因此其更接近于黏弹性介质,在地震处理中补偿这些黏性影响是十分必要的。基于波动方程的叠前深度偏移中进行吸收衰减补偿更准确,也更有物理意义,直接求解双程波动方程的逆时偏移(RTM)能够成像大倾角复杂构造,具有诸多优势。然而当考虑吸收衰减补偿时通常会产生不稳定现象,大部分研究都是在逆时偏移的波场延拓中进行波数域的低通滤波来解决这个问题。本文采用广义标准线性固体的黏声波动方程进行吸收衰减补偿的Q--RTM方法,通过添加正则化项的方式来稳定延拓过程。添加正则化项本质上是低通滤波,滤波窗口是指数形式,在时空域有明确的形式,可以阻止发生高频不稳定。与直接滤波相比,正则化参数可以是空变的,因此比较适合剧烈变化的区域,我们还发现震源归一化互相关成像条件更适合Q--RTM方法。 相似文献
5.
Reverse-time migration (RTM) is based on seismic numerical modeling algorithms, and the accuracy and efficiency of RTM strongly depend on the algorithm used for numerical solution of wave equations. Finite-difference (FD) methods have been widely used to solve the wave equation in seismic numerical modeling and RTM. In this paper, we derive a series of time–space domain staggered-grid FD coefficients for acoustic vertical transversely isotropic (VTI) equations, and adopt these difference coefficients to solve the equations, then analyze the numerical dispersion and stability, and compare the time–space domain staggered-grid FD method with the conventional method. The numerical analysis results demonstrate that the time–space domain staggered-grid FD method has greater accuracy and better stability than the conventional method under the same discretizations. Moreover, we implement the pre-stack acoustic VTI RTM by the conventional and time–space domain high-order staggered-grid FD methods, respectively. The migration results reveal that the time–space domain staggered-grid FD method can provide clearer and more accurate image with little influence on computational efficiency, and the new FD method can adopt a larger time step to reduce the computation time and preserve the imaging accuracy as well in RTM. Meanwhile, when considering the anisotropy in RTM for the VTI model, the imaging quality of the acoustic VTI RTM is better than that of the acoustic isotropic RTM. 相似文献
6.
In real strata anisotropy and viscosity extensively exists. They degraded waveforms in amplitude, resulting in which reducing of image resolution. To obtain high-precision imaging of deep reservoirs, we extended the separated viscous and anisotropic reverse time migration (RTM) to a stable viscoacoustic anisotropic RTM for vertical transverse isotropic (VTI) media, based on single generalized standard and linear solid (GSLS) media theory.. We used a pseudo-spectral method to develop the numerical simulation. By introducing a regularization operator to eliminate the high-frequency instability problem, we built a stable inverse propagator and achieved viscoacoustic VTI media RTM. High-resolution imaging results were obtained after correcting for the effects of anisotropy and viscosity. Synthetic tests verify the validity and accuracy of algorithm. 相似文献
7.
时间域常Q黏声波方程,由于含分数阶时间导数项,数值求解需要大量内存,计算效率低,不利于地震偏移的实施.通过一系列近似,可将该方程简化为介质频散效应和衰减效应解耦的分数阶拉普拉斯算子黏声波方程,数值求解内存需求少,计算效率高.本文采用交错网格有限差分逼近时间导数,改进的伪谱法计算空间导数,PML吸收边界去除边界反射,对该方程进行数值离散和地震正演模拟,开展地震数据的黏声介质逆时偏移,实现波场逆时延拓过程中同时完成频散校正和衰减补偿.改善深层构造的成像精度,数值结果表明,基于分数阶拉普拉斯算子解耦的黏声介质地震正演模拟与逆时偏移可大幅度提高地震模拟计算效率,偏移剖面明显优于常规声波偏移剖面,极大改善深层构造的成像品质. 相似文献
8.
高斯束逆时偏移是一种兼具计算效率和成像精度的深度域成像方法,能够面向目标成像.地下介质中黏滞性普遍存在,利用传统各向同性或完全弹性的成像方法处理黏滞性探区的数据会降低分辨率,并导致成像位置不准确和振幅欠估计等问题.本文在高斯束逆时偏移的基础上,通过对震源点和检波点处的波场进行衰减补偿,并结合高斯束求解时的角度信息,实现了黏声介质角度域高斯束逆时偏移方法.最后通过模型和实际资料试算对本文方法的正确性和适用性进行了验证.试算结果表明:相比于声波高斯束逆时偏移,本文方法能够对黏滞性引起的吸收衰减进行有效补偿,同时提取的角度域共成像点道集(ADCIGs)不仅可以用于分角度叠加成像压制成像噪声,而且能够为后续的偏移速度分析提供支撑. 相似文献
9.
叠前深度偏移理论及方法一直是地震数据成像中研究的热点问题.业界对单程波叠前深度偏移方法和逆时深度偏移开展了深入的研究,但对双程波方程波场深度延拓理论及成像方法的研究还鲜有报道.本文以地表记录的波场值为基础,利用单程波传播算子估计波场对深度的偏导数,为在深度域求解双程波方程提供充分的边界条件,并提出利用矩阵分解理论实现双程波方程的波场深度外推.通过对强速度变化介质中传播波场的计算,与传统的单程波偏移方法相比,本文提出的偏移方法计算的波场与常规有限差分技术计算的波场相一致,证明了本方法计算的准确性.通过对SEAM模型的成像,在相同的成像参数下,与传统的单程波偏移算法和逆时深度偏移算法方法相比,本文提出的偏移方法能够提供更少的虚假成像和更清晰的成像结果.本文所提偏移算法具有深度偏移和双程波偏移的双重特色,推动和发展了双程波叠前深度偏移的理论和实践. 相似文献
10.
各向异性研究对地下介质精确成像有着重要的意义,在当前计算机硬件迅速发展及宽方位地震数据采集日益普遍的情况下,成像必须考虑介质的各向异性.逆时偏移是基于双程波动方程的较为精确的数值解的成像方法,所以相对于其他地震成像方法,它具有很大的优势,譬如不受反射界面的倾角限制、偏移速度结构合适时能够使回转波及多次波正确成像.在各向同性介质中,可使用标量波方程来模拟波场.而在各向异性介质中,P波和SV波是相互耦合的,即不存在单纯的标量波传播,通常利用能代表耦合波场中P波分量运动学特征的拟声波(qP波)进行偏移成像.本文中,我们推导出了TTI介质下qP波控制方程.该方程可采用显式有限差分格式进行求解.通过声学近似,若沿对称轴方向的剪切波速度为零,对于对称轴方向不变且ε≥δ的模型来说,可得到稳定的数值解.但对于TTI介质来说,由于沿对称轴方向各向异性参数是变化的,声学近似会引起波场传播及数值计算的不稳定.因此,我们提出了正则化有限横波的方法,很好地解决了这一问题.最后,给出了Foothill模型的测试结果及某探区实际资料试算结果,展示了采用这个方程进行复杂TTI模型正演和高质量逆时偏移成像结果,证实了该方法的正确性和实际资料应用中的有效性. 相似文献
11.
近年来,面向实际应用的TI介质准P波正演模拟与逆时偏移成像技术受到空前的关注.基于常规耦合型传播方程的正演模拟方法不仅存在伪横波及频散假象干扰,而且还遭受模型参数限制(η0)和不稳定影响;而纯qP波方程的推导繁琐,且由于方程中包含拟微分算子造成求解难度大且精度有限.为此,本文首先构建了一种适用于任意TI介质的纯qP波传播算子,然后借助Low-rank分解求取该算子中的空间-波数域矩阵,同时引入Cerjan衰减边界条件来压制边界反射干扰,最终实现了一种间接的纯qP波波场外推方案,并将其成功应用于复杂TI介质正演模拟与逆时偏移成像中.通过开展数值模拟,并与其他方法对比表明:①该方法既避免了纯qP波方程的繁琐推导,又克服了耦合型方程对模型参数的限制;②还彻底消除了残余伪横波噪音及数值频散;③且能适应较大时间或空间步长及高频震源,是一种相对准确且稳定的各向异性纵波正演与成像策略. 相似文献
12.
有限差分法广泛应用于地震波场的数值延拓,确定合适的有限差分算子以减小数值频散是有限差分法的一个重要研究内容。近年来为了进一步抑制数值频散和增加时间步长,新的有限差分模板得到了应用,对于此,前人使用泰勒展开方法和最小二乘方法确定有限差分算子系数。本文在以前工作的基础上,使用改进的线性方法确定新模板的有限差分系数,并与传统模板线性方法进行对比;通过频散分析和正演模拟验证出新模板线性方法能够更好地保持频散关系,在相同的精度下效率提高了一倍,从而说明了改进的线性方法的有效性。 相似文献
13.
叠前逆时偏移技术是解决地震成像问题的有力工具,但由于计算量大、成像噪音以及存储量大等原因没有得到广泛的应用.本文给出了逆时偏移的实现过程,分析了高阶有限差分格式的稳定性与频散关系.针对叠前逆时偏移计算量大的问题,使用图形处理器(Graphic Processing Unit,简称GPU)实现算法加速,比传统的CPU计算速度提高了一个数量级.文中对理论模型进行了计算,并与单程波偏移方法做比较,结果表明:叠前逆时偏移有效突破了成像倾角限制,对垂直断层、盐丘空腔内幕等特殊构造成像效果均有显著提高.本文尚未涉及成像噪音去除以及存储量等问题,笔者将另文阐述. 相似文献
14.
叠前逆时偏移是目前成像精度最高的地震偏移方法之一,其实现过程中的一个重要步骤是数值求解全波方程,所以快速有效求解全波方程的数值算法对逆时偏移至关重要. 四阶近似解析辛可分Runge-Kutta (NSPRK) 方法是近年发展的一种具有高效率、高精度的数值求解波动方程的保辛差分方法, 能在粗网格条件下有效压制数值频散, 从而提高计算效率, 节省计算机内存需求量. 本文利用四阶NSPRK方法构造的基本思想,发展了具有六阶空间精度的NSPRK方法,并对新的六阶NSPRK方法进行了详细的稳定性和数值频散分析,以及计算效率比较和波场模拟. 同时将该方法用于声波叠前逆时偏移中, 得到一种时间上保辛、空间具有六阶精度、低数值频散、可应用大步长进行波场延拓并能长时计算的叠前逆时偏移方法,对Sigsbee2B模型进行了偏移成像, 并和四阶NSPRK方法、传统的六阶差分方法、四阶Lax-Wendroff correction (LWC) 方法进行了对比. 数值结果表明, 基于六阶NSPRK方法的叠前逆时偏移能得到更好的成像结果, 是一种优于四阶NSPRK方法、传统的六阶差分方法、四阶LWC叠前逆时偏移的方法, 尤其是在粗网格情况下具有更明显的优越性. 相似文献
15.
各向异性介质纯P波方程完全不受横波的干扰,在一定程度上可以减缓由于介质各向异性引起的数值不稳定,本文推导了具有垂直对称轴的横向各向同性(VTI)介质纯P波一阶速度-应力方程.由于纯P波方程存在一个分数形式的伪微分算子,无法直接采用有限差分法求解.针对该问题,本文采用伪谱法和高阶有限差分法联合求解波动方程,重点分析了混合法求解纯P波一阶速度-应力方程的稳定性问题,并给出了混合法求解纯P波方程的稳定性条件.数值模拟结果表明纯P波方程伪谱法和高阶有限差分混合法能够进行复杂介质的正演模拟,在强变速度、变密度的地球介质中仍然具有较好的稳定性. 相似文献
16.
常规τ值法假设应力松弛时间与应变延迟时间近似相等, 造成了常Q模型拟合精度低. 本文利用精确的广义流变体模型Q值计算公式, 研究改进的τ值法求解常Q模型参数. 根据地震波散射理论, 推导了基于广义流变体模型的黏滞性介质一阶波恩近似方程, 结合位移-速度关系得到了含卷积完全匹配层边界条件的黏滞性介质应力-速度方程的一阶波恩近似表达式. 通过数值实验验证并对比了黏滞性介质中全波波动方程、 一阶波恩近似方程以及单程波波动方程的波场特征, 讨论了基于流变体模型的黏滞性介质一阶波恩近似方程对速度扰动和Q扰动的适应性, 以及对旅行时和振幅精度的影响. 相似文献
17.
本文基于Kjartansson常Q模型理论,推导了常Q衰减介质中黏声波和黏弹性波的速度-应力方程,并采用基于二项式窗函数的优化交错网格有限差分方法进行了数值模拟,同时引入不分裂的复频移卷积完全匹配层(CPML)吸收边界条件,以消除边界反射.使用基于自适应时间步长记忆方法的中心差分近似时间分数阶导数,与常用的短时记忆方法相比,提高了波动方程的离散化精度和计算效率.通过对比均匀模型下声波的数值解与解析解,验证了算法的精确性,并进一步分析了不同品质因子下地震波的频散及衰减特征.对BP盐丘模型的数值模拟结果可以较好地反映本文数值方法对复杂介质的适应性及频散压制效果. 相似文献
18.
逆时偏移被认为是对地下复杂构造进行成像的精确偏移方法,尤其是能够有效地对回转波、绕射波、多次波等各种波动现象进行成像.近几年来随着并行计算机和存储设备的快速发展,逆时偏移方法备受关注.本文采用伪谱法实现了VSP逆时偏移,该方法不仅实施简便,而且计算效率高,精度好.并运用反周期扩展法来消除伪谱法中特殊的周期性边界效应问题.对VSP绕射点模型进行试算,分析了因覆盖次数不足在近井区产生的假象问题.对地堑模型和半圆隆起模型也获得了较好的VSP逆时偏移成像效果.并分别对VSP全波波场及分离出的上行波场进行了逆时偏移成像,可明显发现直达波在炮点和检波点位置处收敛成像,也产生了很强的成像噪声.最后对某地区实际观测的VSP资料进行了逆时偏移成像,并与Kirchhoff法VSP偏移结果和地面地震偏移结果进行了对比,显示了VSP逆时偏移在近井区成像上的优势. 相似文献
19.
随着计算机硬件技术的发展以及高分辨率勘探需求的增加,我们希望能够更准确地模拟地下介质,得到更丰富的地层信息.然而,传统的声学假设并不能描述实际地层所存在各向异性和黏滞性,使得成像分辨率较低.为了实现深部储层的高精度成像,本文同时考虑了介质的各向异性和黏滞性,从TI介质弹性波的基本理论出发,结合各向异性GSLS理论,并通过声学近似方法导出基于GSLS模型的各向异性衰减拟声波方程.数值模拟表明该方程既能准确地描述各向异性介质下的准P波运动学规律,又能体现地层的吸收衰减效应;模型逆时偏移结果表明,在实现成像过程中考虑各向异性和黏滞性的影响,能对高陡构造清晰成像,且剖面振幅相对均衡,分辨率较高. 相似文献
20.
研究不同偏移方法的成像机理是实现复杂构造高精度成像的前提,研究了两类共炮点域的相关型叠前深度偏移成像方法：基于单程波动方程的叠前深度成像方法和基于双程波动方程的叠前深度成像方法,同时对比了它们在计算效率、数据存储量、成像精度、成像机理、速度敏感性等方面的差异及其共性。以复杂构造模型为例,采用了傅里叶有限差分法（FFD）和逆时成像法（RTM）,这两种方法实现了121个共炮点道集的叠前深度偏移成像处理。计算结果表明,当速度准确时,两种深度域波动方程成像方法均可以恢复出各个地质反射界面,其中逆时偏移对陡倾角成像效果显著,当速度存在百分比误差和随机扰动情况时,逆时成像结果要差于单程波方法,因此,逆时偏移方法对速度的敏感性较大,且低频噪声较为严重。 相似文献