共查询到19条相似文献，搜索用时 421 毫秒
1.
2.
传统的高阶有限差分波动方程数值模拟方法采用高阶差分算子近似空间偏导数,能有效抑制空间频散.然而,传统的有限差分法仅采用二阶差分算子近似时间偏导数,这使得地震波场沿时间外推的精度较低.当采用较大的时间采样间隔,传统的有限差分法模拟波场会出现明显的时间频散,甚至不稳定.本文基于新的差分结构和中心网格剖分,发展了一种空间任意偶数阶精度、时间四阶和六阶精度的时空域有限差分方法.基于对离散后的频散关系进行泰勒展开,本文推导了时空域高阶有限差分算子的差分系数.相速度分析表明时间四阶、六阶精度的差分方法能显著地减小传统时间二阶精度差分方法的时间频散.在相同的精度下与传统差分法比较,本文发展的时间四阶、六阶有限差分方法的计算效率比传统方法高.均匀和非匀均介质中的波场数值模拟实验进一步证实本文研究的时空高阶有限差分方法的优越性. 相似文献
3.
拟合检验计算中固体潮理论值的一阶微商值最优计算方法讨论 总被引：1，自引：0，他引：1
本文综合介绍了固体潮理论值的一阶微商值的四种计算方法:分波法、一阶差分法、天顶距微分法和插值微分法。文中给出了各种方法对微商值及其拟合检验的对比计算结果,讨论了一阶差分法的差分间隔和插值微分法的插值间隔。结果表明,在固体潮资料的拟合检验计算中,固体潮理论值的一阶微商的最优计算方法是天顶距微分法和插值微分法,差分法中,△T的取值不宜小于1秒。 相似文献
4.
声波方程数值模拟已广泛应用于理论地震计算,同时构成了地震逆时偏移成像技术的基础.对于有限差分法而言,在满足一定的稳定性条件时,普遍存在着因网格化而形成的数值频散效应.如何有效地缓解或压制数值频散是有限差分方法研究的关键所在.为精确求解空间偏导数,相继发展了高阶差分格式优化方法和伪谱方法.近期,为更好地缓解数值频散,提出了时间-空间域有限差分方法,该方法采用了泰勒展开近似方法来确定有限差分格式系数,因而只能保证在一定的小范围内很好的拟合波场传播规律.为进一步压制数值频散效应,本文引入了时间-空间域特定波数点满足频散关系的方法,根据震源、波速和网格间距确定波数范围,同时考虑了多个传播角度,然后建立方程确定了相应的有限差分格式系数,使得差分系数能在更大范围符合波场传播规律.通过频散分析和正演模拟,验证了本文方法的有效性. 相似文献
5.
声波方程数值模拟已广泛应用于理论地震计算,同时构成了地震逆时偏移成像技术的基础.对于有限差分法而言,在满足一定的稳定性条件时,普遍存在着因网格化而形成的数值频散效应.如何有效地缓解或压制数值频散是有限差分方法研究的关键所在.有限差分格式分为显式有限差分和隐式有限差分.隐式有限差分能够进一步压制数值频散效应.因此本文提出了给定频率范围满足时间-空间域隐式有限差分频散关系的方法,并根据震源频率、波速和网格间距确定波数范围,在此基础上建立方程确定了相应的隐式有限差分系数,使得差分系数能在更大频率范围符合波场传播规律.通过频散分析和正演模拟,验证了本文方法的有效性. 相似文献
6.
有限差分法求解Helmholtz方程,依赖于两点:1差分格式的构造;2高效的求解算法.本文采用平均导数法离散Helmholtz方程.该差分格式有三点好处:1能适用于横纵不等间距采样;2在完全匹配层区域(PML),差分方程与微分方程逐点相容;3能将一个波长内的采样点数减少至少于4.求解离散的Helmholtz方程的算法一般分为直接法和迭代算法.直接法由于内存需求太大而无法适用于大规模问题;基于Krylov子空间的迭代方法结合多重网格预条件算法是一种快速高效求解方法,然而对于横纵不等间距采样(在多重网格中称为各向异性问题),经典的多重网格方法失效.本文分析了经典多重网格的三个重要组成部分:完全加权限制算子,点松弛技术以及双线性延拓算子,进而采用了半粗化技术代替全粗化技术,线松弛技术代替点松弛技术以及依赖差分算子的延拓算子代替双线性延拓算子,使得各向异性问题变得收敛;而且对于非均匀介质中-低频率的迭代问题,我们获得了较为满意的收敛速度. 相似文献
7.
有限差分方法是波场数值模拟的一个重要方法,但常规的有限差分法本身存在着数值频散问题，会降低波场模拟的精度与分辨率，为了克服常规差分算子的数值频散，本文采用25点优化差分算子，再根据最优化理论求取的优化系数，建立了频率空间域中弹性波波动方程的差分格式；为了消除边界反射，引入最佳匹配层，构造了各向同性介质中弹性波方程在不同边界和角点处的边界条件. 最后由弹性波波动方程和边界条件，通过频率域有限差分法，分别利用不同震源对弹性波在均匀各向同性介质、层状介质及凹陷模型中的传播过程进行了数值正演模拟，得到了单频波波场、时间切片和共炮点道集，为下一步的研究工作（如成像、反演）提供了研究基础. 相似文献
8.
《应用地球物理》2017,(2)
地震波场正演模拟是地震资料处理、解释中最为重要的技术之一。地震波场正演模拟在大时间步长、长时程的波场延拓中,存在计算不稳定的问题。本文基于声波方程的Hamilton表述,在波动方程求解中用辛差分格式进行时间网格离散,用傅里叶有限差分进行空间网格离散,提出一种新的保结构地震波场正演模拟方法一辛格式傅里叶有限差分法,在保证计算精度的同时提高计算的稳定性。利用声学近似处理空间-波数混合域的积分算子,将该方法推广至各向异性介质。给出各向同性和各向异性条件下的地震正演模拟的计算流程,并将本文方法用于BP盐丘、BP TTI等模型的波场正演模拟。数值算例表明本文开发的方法适用于速度变化剧烈的复杂介质地震波场正演模拟,计算精度高,数值频散小,在各向异性介质正演中能够有效避免qSV波残余,在大时间步长的迭代计算中稳定性好。本文为在辛算法的框架下实现高精度地震正演模拟提供了一种新的选择。 相似文献
9.
由于时间域内粘弹性介质的本构方程是一种卷积积分形式,无法将它直接离散化数值求解.本文采用GSLS模型逼近谐振Q模型介质的粘弹性;推导了粘弹性介质中实现纵波和横波分解的等价波动方程.同时给出了等价方程的完全匹配吸收边界(PML)条件公式及相应的交错网格任意偶数阶精度有限差分格式.最后应用交错网格高阶有限差分法,求解等价波动方程.实验显示GSLS模型逼近精度高,吸收边界效果好,能够实现纵、横波的完全分离,可以得到高精度的波场快照和合成记录;并且波场快照和合成记录能较好的反映谐振Q模型介质的粘弹性特征.结果证明GSLS模型能够精确地逼近谐振Q模型的粘弹性. 相似文献
10.
波动方程有限差分法正演模拟,对认识地震波传播规律、进行地震属性研究、地震资料地质解释、储层评价等,均具有重要的理论和实际意义.但有限差分法本身固有存在着数值频散问题,数值频散在正演模拟中是一种严重的干扰,会降低波场模拟的精度与分辨率.针对TI介质波场模拟的交错网格有限差分方法,本文从空间网格离散、时间网格离散和算子近似等三个方面对其产生的数值频散进行了分析,并结合其他学者的研究成果给出了TI介质波场模拟中压制数值频散的方法与策略：在已知介质频散关系时,对差分算子可实施算子校正;通过提高差分方程的阶数来提高波场模拟精度;采用流体力学中守恒式方程的通量校正传输方法来压制波场模拟中的数值频散;在实际正演模拟时,采用交错网格高阶有限差分方程,不仅在空间上采用高阶差分,而且在时间上也要采用高阶差分,否则只在单一方向上（空间或时间）提高方程的阶数对压制数值频散也不会取得理想的效果. 相似文献
11.
叠前逆时偏移是目前成像精度最高的地震偏移方法之一,其实现过程中的一个重要步骤是数值求解全波方程,所以快速有效求解全波方程的数值算法对逆时偏移至关重要. 四阶近似解析辛可分Runge-Kutta (NSPRK) 方法是近年发展的一种具有高效率、高精度的数值求解波动方程的保辛差分方法, 能在粗网格条件下有效压制数值频散, 从而提高计算效率, 节省计算机内存需求量. 本文利用四阶NSPRK方法构造的基本思想,发展了具有六阶空间精度的NSPRK方法,并对新的六阶NSPRK方法进行了详细的稳定性和数值频散分析,以及计算效率比较和波场模拟. 同时将该方法用于声波叠前逆时偏移中, 得到一种时间上保辛、空间具有六阶精度、低数值频散、可应用大步长进行波场延拓并能长时计算的叠前逆时偏移方法,对Sigsbee2B模型进行了偏移成像, 并和四阶NSPRK方法、传统的六阶差分方法、四阶Lax-Wendroff correction (LWC) 方法进行了对比. 数值结果表明, 基于六阶NSPRK方法的叠前逆时偏移能得到更好的成像结果, 是一种优于四阶NSPRK方法、传统的六阶差分方法、四阶LWC叠前逆时偏移的方法, 尤其是在粗网格情况下具有更明显的优越性. 相似文献
12.
有限差分法广泛应用于地震波场的数值延拓,确定合适的有限差分算子以减小数值频散是有限差分法的一个重要研究内容。近年来为了进一步抑制数值频散和增加时间步长,新的有限差分模板得到了应用,对于此,前人使用泰勒展开方法和最小二乘方法确定有限差分算子系数。本文在以前工作的基础上,使用改进的线性方法确定新模板的有限差分系数,并与传统模板线性方法进行对比;通过频散分析和正演模拟验证出新模板线性方法能够更好地保持频散关系,在相同的精度下效率提高了一倍,从而说明了改进的线性方法的有效性。 相似文献
13.
一般偏移算法是用反演算子通过解析方法求解.最小二乘偏移方法采用另一种思路,即采用数值方法,通过解一个线性离散反问题来索求解.这样我们试着寻找一个模型匹配地震数据并能表现出其某些特点来得到偏移图像.最小二乘法能减少偏移赝像,得到更精确的偏移效果.Kirchhoff算子在最小二乘偏移方法中应用较广,但需要较多的迭代次数,而且具有Kirchhoff偏移的缺点.本文把最小二乘方法运用到基于波长延拓的波动方程偏移方法中,为提高最小二乘偏移的效率,可采用效率较高的正传播算子和反传播算子.我们利用效率较高,能适应剧烈横向变速的傅立叶有限差分正传播和反传播算子来做叠后最小二乘偏移.数值实例表明,通过少数的共轭梯度法迭代,就能得到与真实模型差别很微小的偏移效果.对于傅式变换我们采用了数值软件FFTW,其变换速度比常规FFT算法一般要快六倍以上,进一步提高了效率.本文算法很容易在并行机上实现,这些特点在处理大型数据时大有裨益. 相似文献
14.
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. 相似文献
15.
地震波场数值模拟是理解地震波在地下介质中的传播特点，帮助解释观测数据的有效手段，而提高计算精度和运算效率是所有波场数值模拟方法研究所追求的目标.有限差分技术是求解波动方程计算效率最高、应用最为广泛的方法之一.但传统的有限差分技术计算过程中的数值频散问题影响了该技术的计算精度与计算效率.本文通过交错网格高阶有限差分技术与通量校正传输方法（Flux|corrected transport method，FCT）相结合, 对横向各向同性介质（Transverse isotropic medium，TI）一阶速度|应力弹性波动方程组进行了数值求解研究.波场快照数值模拟结果表明，本文研究的数值模拟方法与波动方程二阶有限差分方法、交错网格四阶有限差分方法相比，在压制网格数值频散方面有明显的优势，计算精度提高，而且可以利用较大的空间步长，提高计算效率. 相似文献
16.
J. GAZDAG 《Geophysical Prospecting》1980,28(1):60-70
A stacked seismic section represents a wave-field recorded at regularly spaced points on the surface. The seismic migration process transforms this recorded data into a reflectivity display. In recent years, Jon F. Claerbout and his co-workers developed migration techniques based on the numerical approximation of the wave equation by finite difference methods. This paper describes an alternative method, termed ASD (for Accurate Space Derivative), and its application to the wave equation migration problem. In this approach to the numerical solution of partial differential equations, partial derivatives are computed by finite Fourier transform methods. This migration method can accommodate media with vertical as well as horizontal velocity variations. 相似文献
17.
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. 相似文献
18.
研究不同偏移方法的成像机理是实现复杂构造高精度成像的前提,研究了两类共炮点域的相关型叠前深度偏移成像方法：基于单程波动方程的叠前深度成像方法和基于双程波动方程的叠前深度成像方法,同时对比了它们在计算效率、数据存储量、成像精度、成像机理、速度敏感性等方面的差异及其共性。以复杂构造模型为例,采用了傅里叶有限差分法（FFD）和逆时成像法（RTM）,这两种方法实现了121个共炮点道集的叠前深度偏移成像处理。计算结果表明,当速度准确时,两种深度域波动方程成像方法均可以恢复出各个地质反射界面,其中逆时偏移对陡倾角成像效果显著,当速度存在百分比误差和随机扰动情况时,逆时成像结果要差于单程波方法,因此,逆时偏移方法对速度的敏感性较大,且低频噪声较为严重。 相似文献
19.
叠前逆时偏移在理论上是现行偏移方法中最为精确的一种成像方法,其实现过程中的核心步骤之一是波动方程的波场延拓,而波场延拓的本质是求解波动方程,所以精确、快速地求解波动方程对逆时偏移至关重要.本文采用一种基于时空域频散关系的有限差分方法来求解声波方程,分析其频散和稳定性,实现波场数值模拟,并将分析和模拟结果与传统有限差分法进行对比.分析结果和模型数值模拟结果都表明时空域有限差分法模拟精度更高、稳定性更好.将时空域高阶有限差分法应用到叠前逆时偏移波场延拓的方程求解中,然后再利用归一化互相关成像条件成像,理论模型数据偏移处理获得了精度更高的成像.同时,在逆时偏移波场延拓的实现中,采用自适应变长度的空间差分算子求解空间导数的有限差分策略,在不影响数值模拟和成像精度的前提下,有效地提高了计算效率. 相似文献