共查询到20条相似文献，搜索用时 781 毫秒
1.
2.
双相介质地震波场数值模拟的迭积微分算子及其PML边界条件 总被引：1，自引：1，他引：0
将基于计算数学中Forsyte 广义正交多项式的迭积微分算子引入到地震波动方程的一阶速度--应力方程的空间微分运算中去,并采用时间错格有限差分算子替代传统的差分算子以匹配高精度的空间迭积微分算子,从而发展一种全新的地震波场正演模拟方法,来解决复杂非均匀介质模型中的波场传播问题.为了大幅衰减人工边界引起的反射,本文将完全匹配层(Perfectly Matched Layer,PML)吸收边界条件引入到所构建的方法中,以解决迭积微分算子法的边界问题.以二维波动方程为例,用迭积微分算子法实现了双相介质的地震波场正演模拟,模拟结果表明,双相介质模型较好地解释了含流体孔隙特性.同时也表明迭积微分算子法是一种非常实用、有效的数值模拟方法. 相似文献
3.
本文针对流固边界耦合介质提出了一种高效、稳定的正演数值模拟方法. 首先,从一阶位移-应力弹性波方程出发,基于海底流固边界的位移和应力的连续性条件,采用三次样条海底界面定量表征方法,推导出不规则海底界面下流固边界耦合介质中的地震波波动方程;其次,通过空间微分的高阶差分格式提高数值模拟的空间精度,并结合已推导的地震波波动方程,将四阶时间微分转换至高阶空间微分,进一步提高了数值模拟的时间精度;最后,在与标量波波动方程数值模拟结果对比分析的基础上,分别利用简单的水平层状模型和复杂海底模型,验证和讨论了本文提出的流固边界耦合介质高阶有限差分地震波正演模拟方法的有效性和准确性. 相似文献
4.
早期的褶积微分算子法都是基于正反傅立叶变换而实现的,其精度比四阶有限差分稍高。本文将计算数学中的Forsyte广义正交多项式微分算子与褶积算子相结合,构建了一个新的快速、高精度褶积微分算子,其计算结果非常接近实验函数微分的精确值,精度与l6阶有限差分相当。粘弹性波动方程更真实地描述了实际地下介质中弹性波的传播规律及其波场特征。本文以二维粘弹性波动方程为例,推导了粘弹性介质波动方程的离散格式,用迭积微分算子法实现了粘弹性介质的地震波场正演模拟,并对其波传播特征进行了分析。计算结果表明该算法能正确模拟粘弹性介质中的地震波,正确地反映粘弹性介质中波场的传播规律。 相似文献
5.
从各向同性介质中波场数值模拟的褶积微分算子法出发,推导出了各向异性双相介质中波场传播数值计算的褶积新算法.将常见的二阶微分Biot波动方程用等效的一阶速度—应力双曲方程表示,其中未知的波场向量包括固相和流体的速度分量和应力分量,由此对方程的时间项使用交错网格差分方法计算,而对空间项则采用褶积微分算法进行求解.对各向异性双相介质在单层介质模型和双层介质模型中的波场特征进行了研究.研究的结果显示,在两层介质分界面上当地震波产生反射时能观测到两类纵波和横波,并且在衰减系数大的介质里慢纵波很难见到. 相似文献
6.
7.
《地球物理学进展》2016,(2)
地震波场数值模拟不仅是研究复杂地区地震资料采集、处理和解释的有效辅助手段,而且是研究地球深部精细构造和地球深部探测的有效工具.射线理论和波动方程理论是地震波场数值模拟的理论基础.射线理论主要刻画地震波在介质中传播的走时场、地震射线等运动学属性;波动方程理论通过求解波动方程来描述地震波在介质中传播的弹性动力学响应(能量衰减、相位特征、偏振属性、以及全波形等).基于波动方程理论的波场数值模拟由于能够引入丰富的波场信息,使得人们对不同介质中地震波的传播过程有了较全面的了解.本文以二维层状均匀介质模型为例,通过射线追踪法和交错网格有限差分法模拟得到的波场快照图、单炮地震记录剖面、合成理论地震图的分析比较,不但对地震波在各向同性层状均匀介质中的传播规律和特点有了深刻的理解和认识,同时又可以相互验证两种不同方法的正确性和有效性. 相似文献
8.
本文提出了一种基于非结构化网格的海洋电磁有限单元正演算法.为了回避场源奇异性,文中选用二次场算法,将背景电阻率设置为水平层状且各向异性,场源在水平层状各向异性介质中所激发的一次场通过汉克尔积分得到.基于Coulomb规范得到二次矢量位和标量位所满足的Maxwell方程组,通过Galerkin加权余量法形成大型稀疏有限元方程,采用不完全LU分解(ILU)预条件因子的quasi-minimum residual(QMR)迭代解法对有限元方程进行求解得到二次矢量位和标量位;进而,利用滑动平均方法得到二次矢量位和标量位在空间的导数,由此得到二次电磁场;通过一维模型对算法的可靠性进行验证,与此同时,针对实际复杂海洋电磁模型,比较有限元模拟结果与积分方程模拟结果,进一步验证算法精度.若干计算结果均表明,文中算法具有良好的通用性,适用于井中电磁、航空电磁,环境地球物理等非均匀且各向异性介质中的电磁感应基础研究. 相似文献
9.
对横向变速介质中振幅保真的共炮检距波场延拓方法进行了研究．在横向线性变速条件下,通过引入新的偏微分方程,给出了该介质中振幅保真的波场延拓方程．这一波场延拓方程是均匀介质中振幅保真延拓方程在横向线性变速介质中的推广．首先由线性变速介质中地震波射线的走时函数,证明了该延拓方程的运动学等价方程对横向线性变速介质中的所有炮检距和各种反射层倾角都适用．然后,通过对模型数据的计算并与传统的ＤＭＯ方法比较,证明了这一波场延拓方程能很好地保持波前的振幅特性,特别是能保持焦散点的波前的特性． 相似文献
10.
三维声波方程相比二维声波方程能够更好的模拟三维空间的地震波传播,模拟标量近似下的弹性波在三维复杂介质的传播过程.基于非规则网格的正演模拟方法的格子法可以处理很好的刻画起伏地表、速度间断面等复杂构造,但是这类方法需要大量的几何描述来描述网格.本文提出了三维六面体双重网格的格子法来模拟声波方程,一方面该方法继承了格子法能够灵活处理自由表面和速度间断面的特性.另一方面,该方法通过双重网格的实现极大的减少了几何描述文件的大小,可以最大的实现GPU加速,实现粗粒度并行,在节省了几何描述空间的同时达到了很高的加速比. 相似文献
11.
与其他偏移方法相比,逆时偏移基于精确的波动方程而不是对其近似,用时间外推来代替深度外推.因此,它具有良好的精度,不受地下构造倾角和介质横向速度变化的限制.激发时间成像条件的求取是叠前逆时偏移的难点之一,本文采用求解程函方程的方法得到地下各点的初至波走时,以此作为叠前逆时偏移的成像条件.基于任意矩形网格和局部平面波前近似的有限差分初至波走时计算方法精度较高并适用于强纵横向变速的复杂介质.试算结果表明,在复杂介质模型中利用叠前逆时深度偏移收到了很好的成像效果. 相似文献
12.
13.
This paper discusses Born/Rytov approximation tomographic velocity inversion methods constrained by the Fresnel zone. Calculations of the sensitivity kernel function and traveltime residuals are critical in tomographic velocity inversion. Based on the Born/Rytov approximation of the frequency-domain wave equation, we derive the traveltime sensitivity kernels of the wave equation on the band-limited wave field and simultaneously obtain the traveltime residuals based on the Rytov approximation. In contrast to single-ray tomography, the modified velocity inversion method improves the inversion stability. Tests of the near-surface velocity model and field data prove that the proposed method has higher accuracy and Computational efficiency than ray theory tomography and full waveform inversion methods. 相似文献
14.
间断有限元方法(Discontinuous Galerkin method,简称DGM)在求解地震波动方程时具有低数值频散、网格剖分灵活等优点,因此,为适应数值模拟对模拟精度和复杂地质结构的要求,本文提出一种新的加权Runge-Kutta间断有限元(weighted Runge-Kutta discontinuous Galerkin,简称WRKDG)方法,用于求解三维D′Alembert介质中声波方程.本文不仅详细推导了其数值格式,特别地,根据常微分方程理论给出了满足数值稳定性条件的一般经验公式,并首次对该方法的数值频散和耗散进行了深入分析,且考虑了耗散参数对结果的影响.同时,我们也对该方法进行了精度测试,并分析了3D情形下WRKDG方法的并行加速比,结果表明3D WRKDG方法具有良好的并行性.最后,我们给出了包含均匀模型、非规则几何模型以及非均匀Marmousi模型在内的数值模拟算例.结果表明,该方法不仅计算准确,能与解析解很好地吻合,且能有效模拟包含球体在内的非规则模型及非均匀Marmousi模型中的衰减声波波场.数值模拟实验进一步验证了WRKDG方法在求解三维D′Alembert介质中声波方程时的正确性和有效性,并获得了对这种强衰减介质中波传播特征的规律性新认识. 相似文献
15.
Tariq Alkhalifah 《Geophysical Prospecting》2002,50(4):373-382
A linearized eikonal equation is developed for transversely isotropic (TI) media with a vertical symmetry axis (VTI). It is linear with respect to perturbations in the horizontal velocity or the anisotropy parameter η. An iterative linearization of the eikonal equation is used as the basis for an algorithm of finite-difference traveltime computations. A practical implementation of this iterative technique is to start with a background model that consists of an elliptically anisotropic, inhomogeneous medium, since traveltimes for this type of medium can be calculated efficiently using eikonal solvers, such as the fast marching method. This constrains the perturbation to changes in the anisotropy parameter η (the parameter most responsible for imaging improvements in anisotropic media). The iterative implementation includes repetitive calculation of η from traveltimes, which is then used to evaluate the perturbation needed for the next round of traveltime calculations using the linearized eikonal equation. Unlike isotropic media, interpolation is needed to estimate η in areas where the traveltime field is independent of η, such as areas where the wave propagates vertically.
Typically, two to three iterations can give sufficient accuracy in traveltimes for imaging applications. The cost of each iteration is slightly less than the cost of a typical eikonal solver. However, this method will ultimately provide traveltime solutions for VTI media. The main limitation of the method is that some smoothness of the medium is required for the iterative implementation to work, especially since we evaluate derivatives of the traveltime field as part of the iterative approach. If a single perturbation is sufficient for the traveltime calculation, which may be the case for weak anisotropy, no smoothness of the medium is necessary. Numerical tests demonstrate the robustness and efficiency of this approach. 相似文献
Typically, two to three iterations can give sufficient accuracy in traveltimes for imaging applications. The cost of each iteration is slightly less than the cost of a typical eikonal solver. However, this method will ultimately provide traveltime solutions for VTI media. The main limitation of the method is that some smoothness of the medium is required for the iterative implementation to work, especially since we evaluate derivatives of the traveltime field as part of the iterative approach. If a single perturbation is sufficient for the traveltime calculation, which may be the case for weak anisotropy, no smoothness of the medium is necessary. Numerical tests demonstrate the robustness and efficiency of this approach. 相似文献
16.
E. BRÜCKL 《Geophysical Prospecting》1987,35(9):973-992
We consider multiply covered traveltimes of first or later arrivals which are gathered along a refraction seismic profile. The two-dimensional distribution of these traveltimes above a coordinate frame generated by the shotpoint axis and the geophone axis or by the common midpoint axis and the offset axis is named a traveltime field. The application of the principle of reciprocity to the traveltime field implies that for each traveltime value with a negative offset there is a corresponding equal value with positive offset. In appendix A procedures are demonstrated which minimize the observational errors of traveltimes inherent in particular traveltime branches or complete common shotpoint sections. The application of the principle of parallelism to an area of the traveltime field associated with a particular refractor can be formulated as a partial differential equation corresponding to the type of the vibrating string. The solution of this equation signifies that the two-dimensional distribution of these traveltimes may be generated by the sum of two one-dimensional functions which depend on the shotpoint coordinate and the geophone coordinate. Physically, these two functions may be interpreted as the mean traveltime branches of the reverse and the normal shot. In appendix B procedures are described which compute these two functions from real traveltime observations by a least-squares fit. The application of these regressed traveltime field data to known time-to-depth conversion methods is straightforward and more accurate and flexible than the use of individual traveltime branches. The wavefront method, the plus-minus method, the generalized reciprocal method and a ray tracing method are considered in detail. A field example demonstrates the adjustment of regressed traveltime fields to observed traveltime data. A time-to-depth conversion is also demonstrated applying a ray tracing method. 相似文献
17.
The influence of the mesoscale eddy on seismic wave propagation and seismic imaging in deep sea is investigated. Based on fundamental fluid equations, an appropriate partial differential equation is derived for the acoustic pressure field in water with eddies, including current effects. Seismic wavefields and synthetic seismograms in the center of the eddy are simulated. Numerical experiments demonstrate that velocity variations caused by the eddy can lead to traveltime perturbations. Further, in seismic images, the reflectors below the water layer are positioned incorrectly due to the perturbation of the eddy and this image perturbation depends linearly on the migration velocity of the layer below the corresponding reflector. The zero-offset seismic profiling throughout the affected area of the eddy shows that the maximum traveltime perturbation appears at the center of the eddy and the structure of horizontal reflectors below the water layer are distorted. 相似文献
18.
层析反演是速度建模中最重要的方法之一,结合偏移成像在成像域进行走时层析速度反演是当前比较成熟有效且广泛应用的技术.本文从高斯束偏移成像条件出发,在波动方程的一阶Born近似和Rytov近似下,推导了成像域走时扰动与速度扰动的线性关系,建立了成像域走时层析方程及其显式表达的层析核函数.该核函数的本质是有限频层析核函数,利用该核函数替换常规射线层析核函数可以明显提高层析反演精度.该核函数的计算关键是背景波场格林函数的计算,本文利用高斯束传播算子计算格林函数进而得到走时层析核函数,实现方式灵活高效且计算精度较高.基于高斯束传播算子的偏移成像与层析成像相结合进行深度域建模迭代,体现了速度建模与偏移成像一体化的思想.数值计算及实际数据应用证明了基于高斯束传播算子的成像域走时层析方法的有效性. 相似文献
19.
The complex‐valued first‐arrival traveltime can be used to describe the properties of both velocity and attenuation as seismic waves propagate in attenuative elastic media. The real part of the complex‐valued traveltime corresponds to phase arrival and the imaginary part is associated with the amplitude decay due to energy absorption. The eikonal equation for attenuative vertical transversely isotropic media discretized with rectangular grids has been proven effective and precise to calculate the complex‐valued traveltime, but less accurate and efficient for irregular models. By using the perturbation method, the complex‐valued eikonal equation can be decomposed into two real‐valued equations, namely the zeroth‐ and first‐order traveltime governing equations. Here, we first present the topography‐dependent zeroth‐ and first‐order governing equations for attenuative VTI media, which are obtained by using the coordinate transformation from the Cartesian coordinates to the curvilinear coordinates. Then, we apply the Lax–Friedrichs sweeping method for solving the topography‐dependent traveltime governing equations in order to approximate the viscosity solutions, namely the real and imaginary parts of the complex‐valued traveltime. Several numerical tests demonstrate that the proposed scheme is efficient and accurate in calculating the complex‐valued P‐wave first‐arrival traveltime in attenuative VTI media with an irregular surface. 相似文献
20.
最短路径射线追踪算法,用预先设置的网络节点的连线表示地震波传播路径,当网络节点稀疏时,获得的射线路径呈之字形,计算的走时比实际走时系统偏大. 本文在波前扩展和反向确定射线路径的过程中,在每个矩形单元内,通过对某边界上的已知走时节点的走时进行线性插值,并利用Fermat原理即时求出从该边界到达其他边界节点的最小走时及其子震源位置和射线路径,发展了相应的动态网络算法. 从而克服了最短路径射线追踪算法的缺陷,大大提高了最小走时和射线路径的计算精度. 相似文献