Three-dimensional acoustic wave equation modeling based on the optimal finite-difference scheme

Generally, FD coefficients can be obtained by using Taylor series expansion （TE） or optimization methods to minimize the dispersion error. However, the TE-based FD method only achieves high modeling precision over a limited range of wavenumbers, and produces large numerical dispersion beyond this range. The optimal FD scheme based on least squares （LS） can guarantee high precision over a larger range of wavenumbers and obtain the best optimization solution at small computational cost. We extend the LS-based optimal FD scheme from two-dimensional （2D） forward modeling to three-dimensional （3D） and develop a 3D acoustic optimal FD method with high efficiency, wide range of high accuracy and adaptability to parallel computing. Dispersion analysis and forward modeling demonstrate that the developed FD method suppresses numerical dispersion. Finally, we use the developed FD method to source wavefield extrapolation and receiver wavefield extrapolation in 3D RTM. To decrease the computation time and storage requirements, the 3D RTM is implemented by combining the efficient boundary storage with checkpointing strategies on GPU. 3D RTM imaging results suggest that the 3D optimal FD method has higher precision than conventional methods.

Finite-difference numerical modeling with even-order accuracy in two-phase anisotropic media

To improve the accuracy of the conventional finite-difference method, finitedifference numerical modeling methods of any even-order accuracy are recommended. We introduce any even-order accuracy difference schemes of any-order derivatives derived from Taylor series expansion. Then, a finite-difference numerical modeling method with any evenorder accuracy is utilized to simulate seismic wave propagation in two-phase anisotropic media. Results indicate that modeling accuracy improves with the increase of difference accuracy order number. It is essential to find the optimal order number, grid size, and time step to balance modeling precision and computational complexity. Four kinds of waves, static mode in the source point, SV wave cusps, reflection and transmission waves are observed in two-phase anisotropic media through modeling.

A Hybrid Dung Beetle Optimization Algorithm with Simulated Annealing for the Numerical Modeling of Asymmetric Wave Equations

In the generalized continuum mechanics(GCM)theory framework,asymmetric wave equations encompass the characteristic scale parameters of the medium,accounting for microstructure interactions.This study integrates two theoretical branches of the GCM,the modified couple stress theory(M-CST)and the one-parameter second-strain-gradient theory,to form a novel asymmetric wave equation in a unified framework.Numerical modeling of the asymmetric wave equation in a unified framework accurately describes subsurface structures with vital implications for subsequent seismic wave inversion and imaging endeavors.However,employing finite-difference(FD)methods for numerical modeling may introduce numerical dispersion,adversely affecting the accuracy of numerical modeling.The design of an optimal FD operator is crucial for enhancing the accuracy of numerical modeling and emphasizing the scale effects.Therefore,this study devises a hybrid scheme called the dung beetle optimization(DBO)algorithm with a simulated annealing(SA)algorithm,denoted as the SA-based hybrid DBO(SDBO)algorithm.An FD operator optimization method under the SDBO algorithm was developed and applied to the numerical modeling of asymmetric wave equations in a unified framework.Integrating the DBO and SA algorithms mitigates the risk of convergence to a local extreme.The numerical dispersion outcomes underscore that the proposed SDBO algorithm yields FD operators with precision errors constrained to 0.5‱while encompassing a broader spectrum coverage.This result confirms the efficacy of the SDBO algorithm.Ultimately,the numerical modeling results demonstrate that the new FD method based on the SDBO algorithm effectively suppresses numerical dispersion and enhances the accuracy of elastic wave numerical modeling,thereby accentuating scale effects.This result is significant for extracting wavefield perturbations induced by complex microstructures in the medium and the analysis of scale effects.

3-D acoustic wave equation forward modeling with topography

In order to model the seismic wave field with surface topography, we present a method of transforming curved grids into rectangular grids in two different coordinate systems. Then the 3D wave equation in the transformed coordinate system is derived. The wave field is modeled using the finite-difference method in the transformed coordinate system. The model calculation shows that this method is able to model the seismic wave field with fluctuating surface topography and achieve good results. Finally, the energy curves of the direct and reflected waves are analyzed to show that surface topography has a great influence on the seismic wave＇s dynamic properties.

Frequency domain wave equation forward modeling using gaussian elimination with static pivoting

Frequency domain wave equation forward modeling is a problem of solving large scale linear sparse systems which is often subject to the limits of computational efficiency and memory storage. Conventional Gaussian elimination cannot resolve the parallel computation of huge data. Therefore, we use the Gaussian elimination with static pivoting （GESP） method for sparse matrix decomposition and multi-source finite-difference modeling. The GESP method does not only improve the computational efficiency but also benefit the distributed parallel computation of matrix decomposition within a single frequency point. We test the proposed method using the classic Marmousi model. Both the single-frequency wave field and time domain seismic section show that the proposed method improves the simulation accuracy and computational efficiency and saves and makes full use of memory. This method can lay the basis for waveform inversion.

Mathematical Wave Functions and 3D Finite Element Modelling of the Electron and Positron

The wave/particle duality of particles in Physics is well known. Particles have properties that uniquely characterize them from one another, such as mass, charge and spin. Charged particles have associated Electric and Magnetic fields. Also, every moving particle has a De Broglie wavelength determined by its mass and velocity. This paper shows that all of these properties of a particle can be derived from a single wave function equation for that particle. Wave functions for the Electron and the Positron are presented and principles are provided that can be used to calculate the wave functions of all the fundamental particles in Physics. Fundamental particles such as electrons and positrons are considered to be point particles in the Standard Model of Physics and are not considered to have a structure. This paper demonstrates that they do indeed have structure and that this structure extends into the space around the particle’s center (in fact, they have infinite extent), but with rapidly diminishing energy density with the distance from that center. The particles are formed from Electromagnetic standing waves, which are stable solutions to the Schrödinger and Classical wave equations. This stable structure therefore accounts for both the wave and particle nature of these particles. In fact, all of their properties such as mass, spin and electric charge, can be accounted for from this structure. These particle properties appear to originate from a single point at the center of the wave function structure, in the same sort of way that the Shell theorem of gravity causes the gravity of a body to appear to all originate from a central point. This paper represents the first two fully characterized fundamental particles, with a complete description of their structure and properties, built up from the underlying Electromagnetic waves that comprise these and all fundamental particles.

Hybrid absorbing boundary condition for threedimensional elastic wave modeling

Edge reflections are inevitable in numerical modeling of seismic wavefields, and they are usually attenuated by absorbing boundary conditions. However, the commonly used perfectly matched layer （PML） boundary condition requires special treatment for the absorbing zone, and in three-dimensional （3D） modeling, it has to split each variable into three corresponding variables, which increases the computing time and memory storage. In contrast, the hybrid absorbing boundary condition （HABC） has the advantages such as ease of implementation, less computation time, and near-perfect absorption; it is thus able to enhance the computational efficiency of 3D elastic wave modeling. In this study, a HABC is developed from two-dimensional （2D） modeling into 3D modeling based on the I st Higdon one way wave equations, and a HABC is proposed that is suitable for a 3D elastic wave numerical simulation. Numerical simulation results for a homogenous model and a complex model indicate that the proposed HABC method is more effective and has better absorption than the traditional PML method.

Research on the Seismic Wave Field of Karst Cavern Reservoirs near Deep Carbonate Weathered Crusts

Fracture and cavern hydrocarbon reservoirs in carbonates are an important pool type worldwide. The karst cavern reservoirs are easiest to identify on seismic reflection data. The prediction, exploration, and development of this type of reservoir require theoretical research on seismic wave fields reflected from complex inhomogeneous media. We compute synthetic seismic sections for fluidfilled cavern reservoirs of various heights and widths using random media models and inhomogeneous media elastic wave equations. Results indicate that even caverns significantly smaller than 1/ 4 wavelength are detectible on conventional band-width seismic sections as diffractions migrated into bead-type events. Diffraction amplitude is a function of cavern height and width. We introduce a width-amplitude factor which can be used to calculate the diffraction amplitude of a cavern with a limited width from the diffraction amplitude computed for an infinitely wide cavern.

Second-generation wavelet finite element based on the lifting scheme for GPR simulation

Ground-penetrating radar(GPR)is a highly efficient,fast and non-destructive exploration method for shallow surfaces.High-precision numerical simulation method is employed to improve the interpretation precision of detection.Second-generation wavelet finite element is introduced into the forward modeling of the GPR.As the finite element basis function,the second-generation wavelet scaling function constructed by the scheme is characterized as having multiple scales and resolutions.The function can change the analytical scale arbitrarily according to actual needs.We can adopt a small analysis scale at a large gradient to improve the precision of analysis while adopting a large analytical scale at a small gradient to improve the efficiency of analysis.This approach is beneficial to capture the local mutation characteristics of the solution and improve the resolution without changing mesh subdivision to realize the efficient solution of the forward GPR problem.The algorithm is applied to the numerical simulation of line current radiation source and tunnel non-dense lining model with analytical solutions.Result show that the solution results of the secondgeneration wavelet finite element are in agreement with the analytical solutions and the conventional finite element solutions,thereby verifying the accuracy of the second-generation wavelet finite element algorithm.Furthermore,the second-generation wavelet finite element algorithm can change the analysis scale arbitrarily according to the actual problem without subdividing grids again.The adaptive algorithm is superior to traditional scheme in grid refinement and basis function order increase,which makes this algorithm suitable for solving complex GPR forward-modeling problems with large gradient and singularity.

Daubechies小波有限元求解GPR波动方程

基于可分离小波理论,由一维Daubechies尺度函数的张量积构造二维Daubechies小波基,并将它作为GPR波动方程求解的插值函数,导出了二维Daubechies小波有限元GPR方程离散格式;通过引入转换矩阵,实现小波系数空间与雷达场值之间转换.引入自由度凝聚技术,有效解决了小波有限元求解中小波单元内部自由度过多的问题,节约了计算量并方便与传统有限元法耦合.然后,详细阐述了Daubechies小波有限元联系系数计算方法,有效解决了小波有限元求解偏微分方程的难点与核心问题.最后,以两个典型GPR模型为例,对比了Daubechies小波有限元与传统有限元的雷达正演剖面图与单道波形图,结果表明:在相同的剖分方式及节点数目条件下,Daubechies小波有限元的紧支性与正交性一定程度上提高了求解效率,它与有限元法求解结果能较好地吻合,验证了Daubechies小波有限元算法的正确性.

非零偏移距DSR叠前深度偏移

提出了非零偏移距域双平方根算子的叠前深度偏移法,该方法根据不同速度模型,通过精确正演模型数据,利用双平方根方程进行偏移,实现了偏移成像和动力学特性研究。研究认为:采用中点-半偏移距域的双平方根算子的叠前深度偏移方法,能够适应横向速度变化,具有精度高、无频散、背景噪声弱、耗时短等优点。对岩性油气藏勘探有重要意义。

基于Tesseral2D的水下砂体地震正演计算

基于水下砂层与围岩的波阻抗差异,应用地震波数值模拟软件Tesseral2D建立含有水、粉细砂、砾砂和砾岩的起伏地层模型,在不同的道间距、最小偏移距、子波频率以及不同岩体波速条件下分别对该模型进行正演计算。一般情况下,在震源频率为400 Hz、最小偏移距为5 m、道间距为1 m或2 m时,地震波响应明显、同相轴清晰、干扰波较少,水下地层界面反映良好。进一步结合南京长江第二大桥桥址区的地质资料,建立含有水、淤泥质粉质黏土、粉细砂、砾砂、砂砾卵石和砾岩的水下砂层模型,应用Tesseral2D软件对该模型进行正演研究。结果表明,在震源频率为400 Hz、道间距为2 m以及最小偏移距为5 m时,地震响应能很好地反映水下各岩层界面,特别是能较明显地圈定水下砂层的厚度和分布范围,为实际水下砂体的地震勘探提供理论依据。

Methods for wave equation prestack depth migration and numerical experiments

In this paper the methods of wave theory based prestack depth migration and their implementation are studied. Using the splitting of wave operator, the wavefield extrapolation equations are deduced and the numerical schemes are presented. The numerical tests for SEG/EAEG model with MPI are performed on the PC-cluster. The numerical results show that the methods of single-shot (common-shot) migration and synthesized-shot migration are of practical values and can be applied to field data processing of 3D prestack depth migration.

SeisMod^(TM)波动方程正演模拟软件在圈闭识别中的应用

准噶尔盆地石西凸起侏罗系八道湾组具很大勘探潜力,现有地震资料对低幅度构造背斜圈闭不能进行有效识别.为满足圈闭识别需要,在不同解释方案基础上,利用SeisModTM波动方程正演模拟软件构建3种地质模型,正演模拟出不同地层构造模式下地震响应特征.通过对正演模拟结果对比分析,为背斜圈闭解释方案提供理论依据,对该区侏罗系八道湾组油气勘探提供有效技术手段.

An unsplit complex frequency-shifted perfectly matched layer for second-order acoustic wave equations

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.

基于有限差分低秩分解策略的黏声各向异性纯qP波正演模拟方法

地下介质的黏弹性和各向异性特征会导致地震波出现相位频散和振幅衰减,如果在地震数据处理中忽略这些影响,那么成像结果会出现同相轴畸变和偏移假象。常规黏声TTI介质拟声波方程可用于模拟地震波在黏声各向异性介质中的传播特征,但存在伪横波干扰和受模型参数限制(各向异性参数ε必须大于δ)等问题。为了解决上述问题,将基于声学近似的纯qP波频散关系与常Q衰减模型相结合,推导了一种黏声TTI介质纯qP波方程。该方程包含解耦的相位频散和振幅衰减项,便于实现衰减补偿逆时偏移。基于新推导的方程,发展了有限差分低秩分解求解策略,实现了黏声TTI介质纯qP波正演模拟。数值模拟结果表明,该方程克服了黏声TTI介质拟声波方程的局限,可较为准确且稳定地模拟地震波在黏声各向异性介质中的传播过程。同时,该有限差分低秩分解求解策略继承了有限差分求解方法高效的优势,相比于传统低秩分解法具有更高的计算效率。

矩形/长方体网格单元时间高阶有限差分波场数值模拟

有限差分法在求解不同类型地震波动方程中广泛应用,提高有限差分在时间域和空间域的地震波场数值模拟精度十分重要。为了提高现有时间—空间域高阶差分方法的精度和灵活性,提出了基于矩形/长方体网格单元的改进时间—空间高阶有限差分方案,分别求解二维和三维声波方程。改进差分法主要通过联合旁轴网格节点和传统轴向节点来共同完成偏导数近似,进一步结合平面波理论和泰勒级数展开求取改进模板中的高阶差分系数。数值精度分析表明,改进差分方案相较于传统方法具有更高的数值精度和稳定性。多个模型算例表明,在相同参数条件下,改进差分方法在兼顾空间模拟精度的前提下能够提高时间模拟精度。此外,改进模板沿不同空间方向可以选用不同的网格步长,有效增加了波场模拟的灵活性,能够为后续的地震成像和反演提供有效波场延拓工具。

A three-dimensional coupled numerical model of nonlinear waves in a harbor

A 3-D time-domain numerical coupled model for nonlinear waves acting on a ship in a harbor has been developed in the present study.The whole domain is divided into the inner domain and the outer domain.The inner domain is the area around the ship,where the flow is expressed by the Laplace equation and numerically solved by the finite element method.The other area is the outer domain,where the flow is described by the higher-order Boussinesq equations and numerically solved by the finite difference method.The matching conditions on the interfaces between the inner domain and the outer domain,the procedure of coupled solution,the length of common domain and the mesh generation in the inner domain are discussed in detail.The other coupled model with the flow in the inner domain governed by the simplified linear Euler equations and relevant physical experiment are adopted to validate the present coupled model,and it is shown that the numerical results of the present model agree with the experimental data,so the present model can be used for the study on the effect of nonlinear waves acting on a fixed ship in a large area and provide a reference for the time-domain simulation of nonlinear wave forces on an arbitrary object in a large harbor and the 3-D district computation in the future.

基于广义坐标变换的起伏地表频率域波动方程正演模拟与逆时偏移

高效精确的波动方程正演模拟是地震偏移成像和参数反演的核心,频率域的有限差分方法在正演模拟具备一定优势,包括能够在多震源模拟中计算效率更高,在选择网格间距时更具有灵活性等.然而,这种方法在处理起伏地表问题上较为困难.为了解决起伏地表情况下的频率域波动方程正演模拟与逆时偏移问题,我们提出通过建立广义坐标系,利用传统笛卡尔坐标系到起伏地表坐标系的映射关系,在起伏地表条件下高精度地模拟地震波传播波场.在此基础上,进一步讨论了实现频率域逆时偏移的可行性.并使用均匀模型和Foothills模型的数值测试来验证本文方法基本流程的有效性和计算效率.结果表明,在进行复杂地表场景的频率域地震成像和反演时,本文方法具有作为正演建模求解器的潜力.