Seismic wavefield modeling is important for improving seismic data processing and interpretation. Calculations of wavefield propagation are sometimes not stable when forward modeling of seismic wave uses large time st...Seismic wavefield modeling is important for improving seismic data processing and interpretation. Calculations of wavefield propagation are sometimes not stable when forward modeling of seismic wave uses large time steps for long times. Based on the Hamiltonian expression of the acoustic wave equation, we propose a structure-preserving method for seismic wavefield modeling by applying the symplectic finite-difference method on time grids and the Fourier finite-difference method on space grids to solve the acoustic wave equation. The proposed method is called the symplectic Fourier finite-difference （symplectic FFD） method, and offers high computational accuracy and improves the computational stability. Using acoustic approximation, we extend the method to anisotropic media. We discuss the calculations in the symplectic FFD method for seismic wavefield modeling of isotropic and anisotropic media, and use the BP salt model and BP TTI model to test the proposed method. The numerical examples suggest that the proposed method can be used in seismic modeling of strongly variable velocities, offering high computational accuracy and low numerical dispersion. The symplectic FFD method overcomes the residual qSV wave of seismic modeling in anisotropic media and maintains the stability of the wavefield propagation for large time steps.展开更多
Laplace–Fourier(L-F)domain finite-difference(FD)forward modeling is an important foundation for L-F domain full-waveform inversion(FWI).An optimal modeling method can improve the efficiency and accuracy of FWI.A fl e...Laplace–Fourier(L-F)domain finite-difference(FD)forward modeling is an important foundation for L-F domain full-waveform inversion(FWI).An optimal modeling method can improve the efficiency and accuracy of FWI.A fl exible FD stencil,which requires pairing and centrosymmetricity of the involved gridpoints,is used on the basis of the 2D L-F domain acoustic wave equation.The L-F domain numerical dispersion analysis is then performed by minimizing the phase error of the normalized numerical phase and attenuation propagation velocities to obtain the optimization coefficients.An optimal FD forward modeling method is finally developed for the L-F domain acoustic wave equation and applied to the traditional standard 9-point scheme and 7-and 9-point schemes,where the latter two schemes are used in discontinuous-grid FD modeling.Numerical experiments show that the optimal L-F domain FD modeling method not only has high accuracy but can also be applied to equal and unequal directional sampling intervals and discontinuous-grid FD modeling to reduce computational cost.展开更多
Global acoustic simulations are significant in revealing the internal and physical structure of the Earth.However,due to the limited flexibility of grids and the difficulties in handling boundaries,the traditional fin...Global acoustic simulations are significant in revealing the internal and physical structure of the Earth.However,due to the limited flexibility of grids and the difficulties in handling boundaries,the traditional finite-difference method(FDM)is usually less used in global simulations.Nevertheless,these issues can be well resolved by employing a multi-block structured grid to discretize circular regions.In this paper,we propose an O-H grid approach to partition the circular region and utilize the curvilinear grid finite-difference method(CGFDM)to solve the acoustic wave equation within this circular domain.By appropriately stretching the grid,the interconnections between each grid block are sufficiently smooth for stable information exchange.To verify the efficacy of this method,we conducted three numerical experiments,by comparing results with alternative approaches.Our test results demonstrate good agreement between our findings and the reference solutions.Since the proposed algorithm can effectively solve wave propagation problems in circular regions,it can contribute to 2D global simulation,particularly in interpreting the Earth’s interior.展开更多
全波形反演作为一种高分辨率地下成像技术,在实际应用中仍面临低频数据缺失、周期跳跃以及噪声干扰等一系列挑战,因此提出一种基于Fourier近似Wasserstein距离的隐式全波形反演方法。在目标函数层面,利用Wasserstein距离与Fourier系数...全波形反演作为一种高分辨率地下成像技术,在实际应用中仍面临低频数据缺失、周期跳跃以及噪声干扰等一系列挑战,因此提出一种基于Fourier近似Wasserstein距离的隐式全波形反演方法。在目标函数层面,利用Wasserstein距离与Fourier系数之间的有界关系,构建基于Fourier近似的低频增强型目标函数。通过对地震数据中低波数成分进行加权增强,提高模型背景速度的恢复质量,从而有效缓解周期跳跃问题。在模型重参数化层面,引入基于正弦激活函数的神经网络对速度模型进行隐式表征,利用其连续的高维映射特性替代传统的离散网格表示,在提高模型表征自由度的同时为反演过程引入由参数化带来的平滑约束。选取Marmousi、Marmousi2和SEG Salt 3个典型模型,分别在低频缺失、线性初始模型、含噪场景和强散射介质等条件下评估方法的有效性。结果表明,在初始模型较差、低频缺失、噪声干扰和强散射等复杂条件下,所提方法均表现出优异的鲁棒性与收敛性,能够有效克服周期跳跃并显著提升反演精度。展开更多
Prestack reverse time migration （RTM） is an accurate imaging method ofsubsurface media. The viscoacoustic prestack RTM is of practical significance because itconsiders the viscosity of the subsurface media. One of t...Prestack reverse time migration （RTM） is an accurate imaging method ofsubsurface media. The viscoacoustic prestack RTM is of practical significance because itconsiders the viscosity of the subsurface media. One of the steps of RTM is solving thewave equation and extrapolating the wave field forward and backward; therefore, solvingaccurately and efficiently the wave equation affects the imaging results and the efficiencyof RTM. In this study, we use the optimal time-space domain dispersion high-order finite-difference （FD） method to solve the viscoacoustic wave equation. Dispersion analysis andnumerical simulations show that the optimal time-space domain FD method is more accurateand suppresses the numerical dispersion. We use hybrid absorbing boundary conditions tohandle the boundary reflection. We also use source-normalized cross-correlation imagingconditions for migration and apply Laplace filtering to remove the low-frequency noise.Numerical modeling suggests that the viscoacoustic wave equation RTM has higher imagingresolution than the acoustic wave equation RTM when the viscosity of the subsurface isconsidered. In addition, for the wave field extrapolation, we use the adaptive variable-lengthFD operator to calculate the spatial derivatives and improve the computational efficiencywithout compromising the accuracy of the numerical solution.展开更多
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...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.展开更多
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 ...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.展开更多
During seismic wave propagation on a free surface, a strong material contrast boundary develops in response to interference by P- and S- waves to create a surface-wave phenomenon. To accurately determine the effects o...During seismic wave propagation on a free surface, a strong material contrast boundary develops in response to interference by P- and S- waves to create a surface-wave phenomenon. To accurately determine the effects of this interface on surface-wave propagation, the boundary conditions must be accurately modeled. In this paper, we present a numerical approach based on the dynamic poroelasticity for a space–time-domain staggered-grid finite-difference simulation in porous media that contain a free-surface boundary. We propose a generalized stess mirror formulation of the free-surface boundary for solids and fluids in porous media for the grid mesh on which lays the free-surface plane. Its analog is that used for elastic media, which is suitable for precise and stable Rayleigh-type surface-wave modeling. The results of our analysis of first kind of Rayleigh （R1） waves obtained by this model demonstrate that the discretization of the mesh in a similar way to that for elastic media can realize stable numerical solutions with acceptable precision. We present numerical examples demonstrating the efficiency and accuracy of our proposed method.展开更多
Numerical simulation in the frequency-space domain has inherent advantages, such as： it is possible to simulate wave propagation from multiple sources simultaneously; there are no cumulative errors; only the interest...Numerical simulation in the frequency-space domain has inherent advantages, such as： it is possible to simulate wave propagation from multiple sources simultaneously; there are no cumulative errors; only the interesting frequencies can be selected; and it is more suitable for wave propagation in viscoelastic media. The only obstacle to using the method is the requirement of huge computer storage. We extend the compressed format for storing the coefficient matrix. It can reduce the required computer storage dramatically. We get the optimal coefficients by least-squares method to suppress the numerical dispersion and adopt the perfectly matched layer （PML） boundary conditions to eliminate the artificial boundary reflections. Using larger grid intervals decreases computer storage requirements and provides high computational efficiency. Numerical experiments demonstrate that these means are economic and effective, providing a good basis for elastic wave imaging and inversion.展开更多
The staggered-grid finite-difference （SGFD） method has been widely used in seismic forward modeling. The precision of the forward modeling results directly affects the results of the subsequent seismic inversion and...The staggered-grid finite-difference （SGFD） method has been widely used in seismic forward modeling. The precision of the forward modeling results directly affects the results of the subsequent seismic inversion and migration. Numerical dispersion is one of the problems in this method. The window function method can reduce dispersion by replacing the finite-difference operators with window operators, obtained by truncating the spatial convolution series of the pseudospectral method. Although the window operators have high precision in the low-wavenumber domain, their precision decreases rapidly in the high-wavenumber domain. We develop a least squares optimization method to enhance the precision of operators obtained by the window function method. We transform the SGFD problem into a least squares problem and find the best solution iteratively. The window operator is chosen as the initial value and the optimized domain is set by the error threshold. The conjugate gradient method is also adopted to increase the stability of the solution. Approximation error analysis and numerical simulation results suggest that the proposed method increases the precision of the window function operators and decreases the numerical dispersion.展开更多
The Z-Axis tiPPer eiectromagnetic （ZTEM） technique is based on a frequency-domain airbome electromagnetic system that measures the natural magnetic field. A survey area was divided into several blocks by using the M...The Z-Axis tiPPer eiectromagnetic （ZTEM） technique is based on a frequency-domain airbome electromagnetic system that measures the natural magnetic field. A survey area was divided into several blocks by using the Maxwell＇s equations, and the magnetic components at the center of each edge of the grid cell are evaluated by applying the staggered-grid finite-difference method. The tipper and its divergence are derived to complete the 3D ZTEM forward modeling algorithm. A synthetic model is then used to compare the responses with those of 2D finite-element forward modeling to verify the accuracy of the algorithm. ZTEM offers high horizontal resolution to both simple and complex distributions of conductivity. This work is the theoretical foundation for the interpretation of ZTEM data and the study of 3D ZTEM inversion.展开更多
Numerical simulation in transverse isotropic media with tilted symmetry axis（TTI） using the standard staggered-grid finite-difference scheme（SSG）results in errors caused by averaging or interpolation. In order to ...Numerical simulation in transverse isotropic media with tilted symmetry axis（TTI） using the standard staggered-grid finite-difference scheme（SSG）results in errors caused by averaging or interpolation. In order to eliminate the errors, a method of rotated staggered-grid finite-difference scheme（RSG） is proposed. However, the RSG brings serious numerical dispersion. The compact staggered-grid finite-difference scheme（CSG） is an implicit difference scheme, which use fewer grid points to suppress dispersion more effectively than the SSG. This paper combines the CSG with the RSG to derive a rotated staggered-grid compact finite-difference scheme（RSGC）. The numerical experiments indicate that the RSGC has weaker numerical dispersion and better accuracy than the RSG.展开更多
Aprogram applying an algorithm of finite-difference time-domain method is established that can be used to calculate the resonant frequencies of cavities of arbitrary modes, coupled with the discrete Fourier transform....Aprogram applying an algorithm of finite-difference time-domain method is established that can be used to calculate the resonant frequencies of cavities of arbitrary modes, coupled with the discrete Fourier transform. Several coaxial resonators including the empty coaxial cavity, re-enrant coaxial cavity and partially stepped resonator are studied with this method, especially the spurious mode resonant frequencies of coaxial cavi- ties. The numerical results thus obtained are shown to be in excellent agreement with those obtainable through rigorous theoretical solutions and experiment results.展开更多
To the most of velocity fields, the traveltimes of the first break that seismic waves propagate along rays can be computed on a 2-D or 3-D numerical grid by finite-difference extrapolation. Under ensuring accuracy, t...To the most of velocity fields, the traveltimes of the first break that seismic waves propagate along rays can be computed on a 2-D or 3-D numerical grid by finite-difference extrapolation. Under ensuring accuracy, to improve calculating efficiency and adaptability, the calculation method of first-arrival traveltime of finite-difference is de- rived based on any rectangular grid and a local plane wavefront approximation. In addition, head waves and scat- tering waves are properly treated and shadow and caustic zones cannot be encountered, which appear in traditional ray-tracing. The testes of two simple models and the complex Marmousi model show that the method has higher accuracy and adaptability to complex structure with strong vertical and lateral velocity variation, and Kirchhoff prestack depth migration based on this method can basically achieve the position imaging effects of wave equation prestack depth migration in major structures and targets. Because of not taking account of the later arrivals energy, the effect of its amplitude preservation is worse than that by wave equation method, but its computing efficiency is higher than that by total Green′s function method and wave equation method.展开更多
Conventional finite-difference(FD)methods cannot model acoustic wave propagation beyond Courant-Friedrichs-Lewy(CFL)numbers 0.707 and 0.577 for two-dimensional(2D)and three-dimensional(3D)equal spacing cases,respectiv...Conventional finite-difference(FD)methods cannot model acoustic wave propagation beyond Courant-Friedrichs-Lewy(CFL)numbers 0.707 and 0.577 for two-dimensional(2D)and three-dimensional(3D)equal spacing cases,respectively,thereby limiting time step selection.Based on the definition of temporal and spatial FD operators,we propose a variable-length temporal and spatial operator strategy to model wave propagation beyond those CFL numbers while preserving accuracy.First,to simulate wave propagation beyond the conventional CFL stability limit,the lengths of the temporal operators are modified to exceed the lengths of the spatial operators for high-velocity zones.Second,to preserve the modeling accuracy,the velocity-dependent lengths of the temporal and spatial operators are adaptively varied.The maximum CFL numbers for the proposed method can reach 1.25 and 1.0 in high velocity contrast 2D and 3D simulation examples,respectively.We demonstrate the effectiveness of our method by modeling wave propagation in simple and complex media.展开更多
Although the phase-shift seismic processing method has characteristics of high accuracy, good stability, high efficiency, and high-dip imaging, it is not able to adapt to strong lateral velocity variation. To overcome...Although the phase-shift seismic processing method has characteristics of high accuracy, good stability, high efficiency, and high-dip imaging, it is not able to adapt to strong lateral velocity variation. To overcome this defect, a finite-difference method in the frequency-space domain is introduced in the migration process, because it can adapt to strong lateral velocity variation and the coefficient is optimized by a hybrid genetic and simulated annealing algorithm. The two measures improve the precision of the approximation dispersion equation. Thus, the imaging effect is improved for areas of high-dip structure and strong lateral velocity variation. The migration imaging of a 2-D SEG/EAGE salt dome model proves that a better imaging effect in these areas is achieved by optimized phase-shift migration operator plus a finite-difference method based on a hybrid genetic and simulated annealing algorithm. The method proposed in this paper is better than conventional methods in imaging of areas of high-dip angle and strong lateral velocity variation.展开更多
A finite-difference approach is used to develop a time-dependent mild-slope equation incorporating the effects of bottom dissipation and nonlinearity. The Enler predietor-corrector method and the three-point finite-di...A finite-difference approach is used to develop a time-dependent mild-slope equation incorporating the effects of bottom dissipation and nonlinearity. The Enler predietor-corrector method and the three-point finite-difference method with varying spatial steps are adopted to discretize the time derivatives and the two-dimensional horizontal ones, respectively, thus leading both the time and spatial derivatives to the second-order accuracy. The boundary conditions for the present model are treated on the basis of the general conditions for open and fixed boundaries with an arbitrary reflection coefficient and phase shift. Both the linear and nonlinear versions of the numerical model are applied to the wave propagation and transformation over an elliptic shoal on a sloping beach, respectively, and the linear version is applied to the simulation of wave propagation in a fully open rectangular harbor. From comparison of numerical results with theoretical or experimental ones, it is found that they are in reasonable agreement.展开更多
基金supported by National Natural Science Foundation of China(41504109,41404099)the Natural Science Foundation of Shandong Province(BS2015HZ008)the project of "Distinguished Professor of Jiangsu Province"
文摘Seismic wavefield modeling is important for improving seismic data processing and interpretation. Calculations of wavefield propagation are sometimes not stable when forward modeling of seismic wave uses large time steps for long times. Based on the Hamiltonian expression of the acoustic wave equation, we propose a structure-preserving method for seismic wavefield modeling by applying the symplectic finite-difference method on time grids and the Fourier finite-difference method on space grids to solve the acoustic wave equation. The proposed method is called the symplectic Fourier finite-difference （symplectic FFD） method, and offers high computational accuracy and improves the computational stability. Using acoustic approximation, we extend the method to anisotropic media. We discuss the calculations in the symplectic FFD method for seismic wavefield modeling of isotropic and anisotropic media, and use the BP salt model and BP TTI model to test the proposed method. The numerical examples suggest that the proposed method can be used in seismic modeling of strongly variable velocities, offering high computational accuracy and low numerical dispersion. The symplectic FFD method overcomes the residual qSV wave of seismic modeling in anisotropic media and maintains the stability of the wavefield propagation for large time steps.
基金National Natural Science Foundation of China(no.41604037)Natural Science Foundation of Hubei Province(no.2022CFB125)+2 种基金Open Fund of Key Laboratory of Exploration Technologies for Oil and Gas Resources(Yangtze University)Ministry of Education(no.K2021-09)College Students'Innovation and Entrepreneurship Training Program(no.2019053)。
文摘Laplace–Fourier(L-F)domain finite-difference(FD)forward modeling is an important foundation for L-F domain full-waveform inversion(FWI).An optimal modeling method can improve the efficiency and accuracy of FWI.A fl exible FD stencil,which requires pairing and centrosymmetricity of the involved gridpoints,is used on the basis of the 2D L-F domain acoustic wave equation.The L-F domain numerical dispersion analysis is then performed by minimizing the phase error of the normalized numerical phase and attenuation propagation velocities to obtain the optimization coefficients.An optimal FD forward modeling method is finally developed for the L-F domain acoustic wave equation and applied to the traditional standard 9-point scheme and 7-and 9-point schemes,where the latter two schemes are used in discontinuous-grid FD modeling.Numerical experiments show that the optimal L-F domain FD modeling method not only has high accuracy but can also be applied to equal and unequal directional sampling intervals and discontinuous-grid FD modeling to reduce computational cost.
基金supported by National Key Research and Development Program of China(No.2022YFF0800602)Guangdong Provincial Key Laboratory of Geophysical High-resolution Imaging Technology(No.2022B1212010002)Shenzhen Science and Technology Program(No.KQTD20170810111725321).
文摘Global acoustic simulations are significant in revealing the internal and physical structure of the Earth.However,due to the limited flexibility of grids and the difficulties in handling boundaries,the traditional finite-difference method(FDM)is usually less used in global simulations.Nevertheless,these issues can be well resolved by employing a multi-block structured grid to discretize circular regions.In this paper,we propose an O-H grid approach to partition the circular region and utilize the curvilinear grid finite-difference method(CGFDM)to solve the acoustic wave equation within this circular domain.By appropriately stretching the grid,the interconnections between each grid block are sufficiently smooth for stable information exchange.To verify the efficacy of this method,we conducted three numerical experiments,by comparing results with alternative approaches.Our test results demonstrate good agreement between our findings and the reference solutions.Since the proposed algorithm can effectively solve wave propagation problems in circular regions,it can contribute to 2D global simulation,particularly in interpreting the Earth’s interior.
文摘全波形反演作为一种高分辨率地下成像技术,在实际应用中仍面临低频数据缺失、周期跳跃以及噪声干扰等一系列挑战,因此提出一种基于Fourier近似Wasserstein距离的隐式全波形反演方法。在目标函数层面,利用Wasserstein距离与Fourier系数之间的有界关系,构建基于Fourier近似的低频增强型目标函数。通过对地震数据中低波数成分进行加权增强,提高模型背景速度的恢复质量,从而有效缓解周期跳跃问题。在模型重参数化层面,引入基于正弦激活函数的神经网络对速度模型进行隐式表征,利用其连续的高维映射特性替代传统的离散网格表示,在提高模型表征自由度的同时为反演过程引入由参数化带来的平滑约束。选取Marmousi、Marmousi2和SEG Salt 3个典型模型,分别在低频缺失、线性初始模型、含噪场景和强散射介质等条件下评估方法的有效性。结果表明,在初始模型较差、低频缺失、噪声干扰和强散射等复杂条件下,所提方法均表现出优异的鲁棒性与收敛性,能够有效克服周期跳跃并显著提升反演精度。
基金This research was supported by the National Nature Science Foundation of China （No. 41074100） and the Program for NewCentury Excellent Talents in the University of the Ministry of Education of China （No. NCET- 10-0812）.
文摘Prestack reverse time migration （RTM） is an accurate imaging method ofsubsurface media. The viscoacoustic prestack RTM is of practical significance because itconsiders the viscosity of the subsurface media. One of the steps of RTM is solving thewave equation and extrapolating the wave field forward and backward; therefore, solvingaccurately and efficiently the wave equation affects the imaging results and the efficiencyof RTM. In this study, we use the optimal time-space domain dispersion high-order finite-difference （FD） method to solve the viscoacoustic wave equation. Dispersion analysis andnumerical simulations show that the optimal time-space domain FD method is more accurateand suppresses the numerical dispersion. We use hybrid absorbing boundary conditions tohandle the boundary reflection. We also use source-normalized cross-correlation imagingconditions for migration and apply Laplace filtering to remove the low-frequency noise.Numerical modeling suggests that the viscoacoustic wave equation RTM has higher imagingresolution than the acoustic wave equation RTM when the viscosity of the subsurface isconsidered. In addition, for the wave field extrapolation, we use the adaptive variable-lengthFD operator to calculate the spatial derivatives and improve the computational efficiencywithout compromising the accuracy of the numerical solution.
基金supported by the National Natural Science Foundation of China(No.41474110)Shell Ph.D. Scholarship to support excellence in geophysical research
文摘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.
文摘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.
基金sponsed by National Natural Science Foundation of China(NSFC,Grant No.41304077)the Natural Basic Research Program of China(the“973 Project,”Grant No.2013CB733303)Postdoctoral Science Foundation of China(Grant No.2014T70740)
文摘During seismic wave propagation on a free surface, a strong material contrast boundary develops in response to interference by P- and S- waves to create a surface-wave phenomenon. To accurately determine the effects of this interface on surface-wave propagation, the boundary conditions must be accurately modeled. In this paper, we present a numerical approach based on the dynamic poroelasticity for a space–time-domain staggered-grid finite-difference simulation in porous media that contain a free-surface boundary. We propose a generalized stess mirror formulation of the free-surface boundary for solids and fluids in porous media for the grid mesh on which lays the free-surface plane. Its analog is that used for elastic media, which is suitable for precise and stable Rayleigh-type surface-wave modeling. The results of our analysis of first kind of Rayleigh （R1） waves obtained by this model demonstrate that the discretization of the mesh in a similar way to that for elastic media can realize stable numerical solutions with acceptable precision. We present numerical examples demonstrating the efficiency and accuracy of our proposed method.
基金supported by the 863 Program (Grant no.2006AA09Z323)the 973 Program (Grant No.2006CB202402)
文摘Numerical simulation in the frequency-space domain has inherent advantages, such as： it is possible to simulate wave propagation from multiple sources simultaneously; there are no cumulative errors; only the interesting frequencies can be selected; and it is more suitable for wave propagation in viscoelastic media. The only obstacle to using the method is the requirement of huge computer storage. We extend the compressed format for storing the coefficient matrix. It can reduce the required computer storage dramatically. We get the optimal coefficients by least-squares method to suppress the numerical dispersion and adopt the perfectly matched layer （PML） boundary conditions to eliminate the artificial boundary reflections. Using larger grid intervals decreases computer storage requirements and provides high computational efficiency. Numerical experiments demonstrate that these means are economic and effective, providing a good basis for elastic wave imaging and inversion.
基金jointly supported by the NSF(No.41720104006)the Strategic Priority Research Program of the Chinese Academy of Sciences(A)(No.XDA14010303)+2 种基金the National Oil and Gas Project(Nos.2016ZX05002-005-007HZ and 2016ZX05014-001-008HZ)the Shandong Innovation Project(No.2017CXGC1602)the Qingdao Innovation Project(Nos.16-5-1-40-jch and 17CX05011)
文摘The staggered-grid finite-difference （SGFD） method has been widely used in seismic forward modeling. The precision of the forward modeling results directly affects the results of the subsequent seismic inversion and migration. Numerical dispersion is one of the problems in this method. The window function method can reduce dispersion by replacing the finite-difference operators with window operators, obtained by truncating the spatial convolution series of the pseudospectral method. Although the window operators have high precision in the low-wavenumber domain, their precision decreases rapidly in the high-wavenumber domain. We develop a least squares optimization method to enhance the precision of operators obtained by the window function method. We transform the SGFD problem into a least squares problem and find the best solution iteratively. The window operator is chosen as the initial value and the optimized domain is set by the error threshold. The conjugate gradient method is also adopted to increase the stability of the solution. Approximation error analysis and numerical simulation results suggest that the proposed method increases the precision of the window function operators and decreases the numerical dispersion.
基金supported by the Natural Science Foundation of China(No.41374078)Geological Survey Projects of Ministry of Land and Resources of China(No.12120113086100 and 12120113101300)
文摘The Z-Axis tiPPer eiectromagnetic （ZTEM） technique is based on a frequency-domain airbome electromagnetic system that measures the natural magnetic field. A survey area was divided into several blocks by using the Maxwell＇s equations, and the magnetic components at the center of each edge of the grid cell are evaluated by applying the staggered-grid finite-difference method. The tipper and its divergence are derived to complete the 3D ZTEM forward modeling algorithm. A synthetic model is then used to compare the responses with those of 2D finite-element forward modeling to verify the accuracy of the algorithm. ZTEM offers high horizontal resolution to both simple and complex distributions of conductivity. This work is the theoretical foundation for the interpretation of ZTEM data and the study of 3D ZTEM inversion.
文摘Numerical simulation in transverse isotropic media with tilted symmetry axis（TTI） using the standard staggered-grid finite-difference scheme（SSG）results in errors caused by averaging or interpolation. In order to eliminate the errors, a method of rotated staggered-grid finite-difference scheme（RSG） is proposed. However, the RSG brings serious numerical dispersion. The compact staggered-grid finite-difference scheme（CSG） is an implicit difference scheme, which use fewer grid points to suppress dispersion more effectively than the SSG. This paper combines the CSG with the RSG to derive a rotated staggered-grid compact finite-difference scheme（RSGC）. The numerical experiments indicate that the RSGC has weaker numerical dispersion and better accuracy than the RSG.
文摘Aprogram applying an algorithm of finite-difference time-domain method is established that can be used to calculate the resonant frequencies of cavities of arbitrary modes, coupled with the discrete Fourier transform. Several coaxial resonators including the empty coaxial cavity, re-enrant coaxial cavity and partially stepped resonator are studied with this method, especially the spurious mode resonant frequencies of coaxial cavi- ties. The numerical results thus obtained are shown to be in excellent agreement with those obtainable through rigorous theoretical solutions and experiment results.
基金National Natural Science Foundation of China (49894190-024) and Geophysical Prospecting Key Laboratory Foundation of China National Petroleum Corporation.
文摘To the most of velocity fields, the traveltimes of the first break that seismic waves propagate along rays can be computed on a 2-D or 3-D numerical grid by finite-difference extrapolation. Under ensuring accuracy, to improve calculating efficiency and adaptability, the calculation method of first-arrival traveltime of finite-difference is de- rived based on any rectangular grid and a local plane wavefront approximation. In addition, head waves and scat- tering waves are properly treated and shadow and caustic zones cannot be encountered, which appear in traditional ray-tracing. The testes of two simple models and the complex Marmousi model show that the method has higher accuracy and adaptability to complex structure with strong vertical and lateral velocity variation, and Kirchhoff prestack depth migration based on this method can basically achieve the position imaging effects of wave equation prestack depth migration in major structures and targets. Because of not taking account of the later arrivals energy, the effect of its amplitude preservation is worse than that by wave equation method, but its computing efficiency is higher than that by total Green′s function method and wave equation method.
基金the National Natural Science Foundation of China(No.41874144)the Research Foundation of China University of PetroleumBeijing at Karamay(RCYJ2018A-01-001).
文摘Conventional finite-difference(FD)methods cannot model acoustic wave propagation beyond Courant-Friedrichs-Lewy(CFL)numbers 0.707 and 0.577 for two-dimensional(2D)and three-dimensional(3D)equal spacing cases,respectively,thereby limiting time step selection.Based on the definition of temporal and spatial FD operators,we propose a variable-length temporal and spatial operator strategy to model wave propagation beyond those CFL numbers while preserving accuracy.First,to simulate wave propagation beyond the conventional CFL stability limit,the lengths of the temporal operators are modified to exceed the lengths of the spatial operators for high-velocity zones.Second,to preserve the modeling accuracy,the velocity-dependent lengths of the temporal and spatial operators are adaptively varied.The maximum CFL numbers for the proposed method can reach 1.25 and 1.0 in high velocity contrast 2D and 3D simulation examples,respectively.We demonstrate the effectiveness of our method by modeling wave propagation in simple and complex media.
基金the Open Fund(PLC201104)of State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation (Chengdu University of Technology)the National Natural Science Foundation of China(No.61072073)the Key Project of Education Commission of Sichuan Province(No.10ZA072)
文摘Although the phase-shift seismic processing method has characteristics of high accuracy, good stability, high efficiency, and high-dip imaging, it is not able to adapt to strong lateral velocity variation. To overcome this defect, a finite-difference method in the frequency-space domain is introduced in the migration process, because it can adapt to strong lateral velocity variation and the coefficient is optimized by a hybrid genetic and simulated annealing algorithm. The two measures improve the precision of the approximation dispersion equation. Thus, the imaging effect is improved for areas of high-dip structure and strong lateral velocity variation. The migration imaging of a 2-D SEG/EAGE salt dome model proves that a better imaging effect in these areas is achieved by optimized phase-shift migration operator plus a finite-difference method based on a hybrid genetic and simulated annealing algorithm. The method proposed in this paper is better than conventional methods in imaging of areas of high-dip angle and strong lateral velocity variation.
基金This work wasjointlysupported by the National Natural Science Foundation of China(Grant No.40106008) the National Natural Science Fundfor Distinguished Young Scholars(Grant No.40225014)
文摘A finite-difference approach is used to develop a time-dependent mild-slope equation incorporating the effects of bottom dissipation and nonlinearity. The Enler predietor-corrector method and the three-point finite-difference method with varying spatial steps are adopted to discretize the time derivatives and the two-dimensional horizontal ones, respectively, thus leading both the time and spatial derivatives to the second-order accuracy. The boundary conditions for the present model are treated on the basis of the general conditions for open and fixed boundaries with an arbitrary reflection coefficient and phase shift. Both the linear and nonlinear versions of the numerical model are applied to the wave propagation and transformation over an elliptic shoal on a sloping beach, respectively, and the linear version is applied to the simulation of wave propagation in a fully open rectangular harbor. From comparison of numerical results with theoretical or experimental ones, it is found that they are in reasonable agreement.
