The pseudospectral method has been applied to the simulation of seismic wave propagation in 2-D global Earth model. When a whole Earth model is considered, the center of the Earth is included in the model and then sin...The pseudospectral method has been applied to the simulation of seismic wave propagation in 2-D global Earth model. When a whole Earth model is considered, the center of the Earth is included in the model and then singularity arises at the center of the Earth where r=0 since the 1/r term appears in the wave equations. In this paper, we extended the global seismic wavefield simulation algorithm for regular grid mesh to staggered grid configuration and developed a scheme to solve the numerical problems associated with the above singularity for a 2-D global Earth model defined on staggered grid using pseudospectral method. This scheme uses a coordinate transformation at the center of the model, in which the field variables at the center are calculated in Cartesian coordinates from the values on the grids around the center. It allows wave propagation through the center and hence the wavefield at the center can be stably calculated. Validity and accuracy of the scheme was tested by compared with the discrete wavenumber method. This scheme could also be suitable for other numerical methods or models parameterized in cylindrical or spherical coordinates when singularity arises at the center of the model.展开更多
Elastic waves in the seabed generated by low-frequency noise radiating from ships are known as ship seismic waves and can be used to detect and identify ships. To obtain the propagation characteristics of ship seismic...Elastic waves in the seabed generated by low-frequency noise radiating from ships are known as ship seismic waves and can be used to detect and identify ships. To obtain the propagation characteristics of ship seismic waves, an algorithm for calculating Seismic waves at the seafloor is presented based on the staggered-grid finite difference method. The accuracy of the algorithm was tested by comparison with analytical solutions. Numerical simulation of seismic waves generated by a low-frequency point sotmd source in a typical shallow sea environment was carried out. Using various source frequencies and locations in the numerical simulation, we show that the seismic waves in the near field are composed mostly of transmitted S-waves and interface waves while transmitted P-waves are weak near the seafloor. However, in the far field, the wave components of the seismic wave are mainly normal modes and interface waves, with the latter being relatively strong in the waveforms, As the source frequency decreases, the normal modes become smaller and the interface waves dominate the time series of the seismic waves.展开更多
Strong ground motion of an earthquake is simulated by using both staggered grid finite difference method (FDM) and stochastic method, respectively. The acceleration time histories obtained from the both ways and their...Strong ground motion of an earthquake is simulated by using both staggered grid finite difference method (FDM) and stochastic method, respectively. The acceleration time histories obtained from the both ways and their response spectra are compared. The result demonstrates that the former is adequate to simulate the low-frequency seismic wave; the latter is adequate to simulate the high-frequency seismic wave. Moreover, the result obtained from FDM can better reflect basin effects.展开更多
A finite volume method based unstructured grid is presented to solve the two dimensional viscous and incompressible flow.The method is based on the pressure-correction concept and solved by using a semi-staggered grid...A finite volume method based unstructured grid is presented to solve the two dimensional viscous and incompressible flow.The method is based on the pressure-correction concept and solved by using a semi-staggered grid technique.The computational procedure can handle cells of arbitrary shapes,although solutions presented in this paper were only involved with triangular and quadrilateral cells.The pressure or pressure-correction value was stored on the vertex of cells.The mass conservation equation was discretized on the dual cells surrounding the vertex of primary cells,while the velocity components and other scale variables were saved on the central of primary cells.Since the semi-staggered arrangement can’t guarantee a strong coupling relationship between pressure and velocity,thus a weak coupling relationship leads to the oscillations for pressure and velocity.In order to eliminate such an oscillation,a special interpolation scheme was used to construct the pressure-correction equation.Computational results of several viscous flow problems show good agreement with the analytical or numerical results in previous literature.This semi-staggered grid method can be applied to arbitrary shape elements,while it has the most efficiency for triangular cells.展开更多
A new 3D finite-difference(FD) method of spatially asymmetric staggered grids was presented to simulate elastic wave propagation in topographic structures.The method approximated the first-order elastic wave equations...A new 3D finite-difference(FD) method of spatially asymmetric staggered grids was presented to simulate elastic wave propagation in topographic structures.The method approximated the first-order elastic wave equations by irregular grids finite difference operator with second-order time precise and fourth-order spatial precise.Additional introduced finite difference formula solved the asymmetric problem arisen in non-uniform staggered grid scheme.The method had no interpolation between the fine and coarse grids.All grids were computed at the same spatial iteration.Complicated geometrical structures like rough submarine interface,fault and nonplanar interfaces were treated with fine irregular grids.Theoretical analysis and numerical simulations show that this method saves considerable memory and computing time,at the same time,has satisfactory stability and accuracy.展开更多
In this paper, the cell face velocities in the discretization of the continu- ity equation, the momentum equation, and the scalar equation of a non-staggered grid system are calculated and discussed. Both the momentum...In this paper, the cell face velocities in the discretization of the continu- ity equation, the momentum equation, and the scalar equation of a non-staggered grid system are calculated and discussed. Both the momentum interpolation and the linear interpolation are adopted to evaluate the coefficients in the discretized momentum and scalar equations. Their performances are compared. When the linear interpolation is used to calculate the coefficients, the mass residual term in the coefficients must be dropped to maintain the accuracy and convergence rate of the solution.展开更多
The numerical dispersion and computational cost are high for conventional Taylor series expansion staggered-grid finite-difference forward modeling owing to the high frequency of the wavelets and the large grid interv...The numerical dispersion and computational cost are high for conventional Taylor series expansion staggered-grid finite-difference forward modeling owing to the high frequency of the wavelets and the large grid intervals. In this study, the cosine-modulated binomial window function （CMBWF）-based staggered-grid finite-difference method is proposed. Two new parameters, the modulated time and modulated range are used in the new window function and by adjusting these two parameters we obtain different characteristics of the main and side lobes of the amplitude response. Numerical dispersion analysis and elastic wavefield forward modeling suggests that the CMBWF method is more precise and less computationally costly than the conventional Taylor series expansion staggered-grid finite-difference method.展开更多
Numerical simulations of a seismic wavefield are important to analyze seismic wave propagation. Elastic-wave equations are used in data simulation for modeling migration and imaging. In elastic wavefield numerical mod...Numerical simulations of a seismic wavefield are important to analyze seismic wave propagation. Elastic-wave equations are used in data simulation for modeling migration and imaging. In elastic wavefield numerical modeling, the rotated staggered-grid method （RSM） is a modification of the standard staggered-grid method （SSM）. The variable-order method is based on the method of variable-length spatial operators and wavefield propagation, and it calculates the real dispersion error by adapting different finite-difference orders to different velocities. In this study, the variable-order rotated staggered-grid method （VRSM） is developed after applying the variable-order method to RSM to solve the numerical dispersion problem of RSM in low-velocity regions and reduce the computation cost. Moreover, based on theoretical dispersion and the real dispersion error of wave propagation calculated with the wave separation method, the application of the original method is extended from acoustic to shear waves, and the calculation is modified from theoretical to time-varying values. A layered model and an overthrust model are used to demonstrate the applicability of VRSM. We also evaluate the order distribution, wave propagation, and computation time. The results suggest that the VRSM order distribution is reasonable and VRSM produces high-precision results with a minimal computation cost.展开更多
A finite-volume formulation is proposed to solve the three-dimensional,non-hydrostatic Navier-Stokes equations on an unstructured,staggered,z-lever grid,with the goal of simulating non-hydrostatic processes in the fre...A finite-volume formulation is proposed to solve the three-dimensional,non-hydrostatic Navier-Stokes equations on an unstructured,staggered,z-lever grid,with the goal of simulating non-hydrostatic processes in the free-surface flows.The advection and diffusion terms in the momentum equation are discretized explicitly with the Eulerian scheme,which has the attractive property of being conservative.An integral method of the top-layer pressure is applied to account for the full effects of non-hydrostatic pressure at the free-surface layer.It is shown that the results obtained with a small number of vertical layers（e.g.,2-3 layers） are in good agreements with experimental data or analytical solutions,demonstrating the efficiency and accuracy of the model in simulating a range of free-surface flow problems including wave motion and tide-induced motion.展开更多
This paper presents a new horizontal staggered grid(LE grid),which defines h at a gridpoint,and both u and v at the same mid-gridpoint along the x and y directions.A general method is used to deduce the dispersion rel...This paper presents a new horizontal staggered grid(LE grid),which defines h at a gridpoint,and both u and v at the same mid-gridpoint along the x and y directions.A general method is used to deduce the dispersion relationships of describing inertia gravity waves on LE grid and Arakawa A―E grids,which are then compared with the analytical solution(AS)in re-solved or under-resolved cases,using two-order central difference or four-order compact differ-ence scheme from the frequency and group velocity.Results show that in both resolved and under-resolved cases,no matter whether two-order central difference or four-order compact dif-ference scheme is used,the frequency and group velocity discrete errors on LE grid in describing inertia gravity waves are smaller than those of Arakawa A―E grids.At the same time,it is only on LE or Arakawa grid C that the employment of a compact difference scheme of higher difference precision can improve their accuracy in describing inertia gravity waves.However,as for the other four grids(Arakawa A,B,D and E),when the difference precision increases,the accuracy of simulating inertia gravity waves decreases.展开更多
Reverse time migration and full waveform inversion involve the crosscorrelation of two wavefields,propagated in the forward-and reverse-time directions,respectively.As a result,the forward-propagated wavefield needs t...Reverse time migration and full waveform inversion involve the crosscorrelation of two wavefields,propagated in the forward-and reverse-time directions,respectively.As a result,the forward-propagated wavefield needs to be stored,and then accessed to compute the correlation with the backward-propagated wavefield.Boundary-value methods reconstruct the source wavefield using saved boundary wavefields and can significantly reduce the storage requirements.However,the existing boundary-value methods are based on the explicit finite-difference(FD)approximations of the spatial derivatives.Implicit FD methods exhibit greater accuracy and thus allow for a smaller operator length.We develop two(an accuracy-preserving and a memory-efficient)wavefield reconstruction schemes based on an implicit staggered-grid FD(SFD)operator.The former uses boundary wavefields at M layers of grid points and the spatial derivatives of wavefields at one layer of grid points to reconstruct the source wavefield for a(2M+2)th-order implicit SFD operator.The latter applies boundary wavefields at N layers of grid points,a linear combination of wavefields at M–N layers of grid points,and the spatial derivatives of wavefields at one layer of grid points to reconstruct the source wavefield(0≤N<M).The required memory of accuracy-preserving and memory-efficient schemes is(M+1)/M and(N+2)/M times,respectively,that of the explicit reconstruction scheme.Numerical results reveal that the accuracy-preserving scheme can achieve accurate reconstruction at the cost of storage.The memory-efficient scheme with N=2 can obtain plausible reconstructed wavefields and images,and the storage amount is 4/(M+1)of the accuracy-preserving scheme.展开更多
Seismic wave modeling is a cornerstone of geophysical data acquisition, processing, and interpretation, for which finite-difference methods are often applied. In this paper, we extend the velocity- pressure formulatio...Seismic wave modeling is a cornerstone of geophysical data acquisition, processing, and interpretation, for which finite-difference methods are often applied. In this paper, we extend the velocity- pressure formulation of the acoustic wave equation to marine seismic modeling using the staggered-grid finite-difference method. The scheme is developed using a fourth-order spatial and a second-order temporal operator. Then, we define a stability coefficient (SC) and calculate its maximum value under the stability condition. Based on the dispersion relationship, we conduct a detailed dispersion analysis for submarine sediments in terms of the phase and group velocity over a range of angles, stability coefficients, and orders. We also compare the numerical solution with the exact solution for a P-wave line source in a homogeneous submarine model. Additionally, the numerical results determined by a Marmousi2 model with a rugged seafloor indicate that this method is sufficient for modeling complex submarine structures.展开更多
In this paper,a fairly simple 3D immersed interface method based on the CG-Uzawa type method and the level set representation of the interface is employed for solving three-dimensional Stokes flow with singular forces...In this paper,a fairly simple 3D immersed interface method based on the CG-Uzawa type method and the level set representation of the interface is employed for solving three-dimensional Stokes flow with singular forces along the interface.The method is to apply the Taylor’s expansions only along the normal direction and incorporate the jump conditions up to the second normal derivatives into the finite difference schemes.A second order geometric iteration algorithm is employed for computing orthogonal projections on the surface with third-order accuracy.The Stokes equations are discretized involving the correction terms on staggered grids and then solved by the conjugate gradient Uzawa type method.The major advantages of the present method are the special simplicity,the ability in handling the Dirichlet boundary conditions,and no need of the pressure boundary condition.The method can also preserve the volume conservation and the discrete divergence free condition very well.The numerical results show that the proposed method is second order accurate and efficient.展开更多
A depth averaged nonlinear k ε model for turbulent flows in complex geometries has been developed in a boundary fitted coordinate system. The SIMPLEC procedure is used to develop an economical discrete method for ...A depth averaged nonlinear k ε model for turbulent flows in complex geometries has been developed in a boundary fitted coordinate system. The SIMPLEC procedure is used to develop an economical discrete method for staggered grids to analyze flows in a 90° bend. This paper describes how to change a program in rectangular coordinate into a boundary fitted coordinate. The results compare well with experimental data for flow in a meandering channel showing the efficiency of the model and the discrete method.展开更多
In modelling elastic wave propagation in a porous medium, when the ratio between the fluid viscosity and the medium permeability is comparatively large, the stiffness problem of Blot＇s poroelastic equations will be e...In modelling elastic wave propagation in a porous medium, when the ratio between the fluid viscosity and the medium permeability is comparatively large, the stiffness problem of Blot＇s poroelastic equations will be encountered. In the paper, a partition method is developed to solve the stiffness problem with a staggered high-order finite-difference. The method splits the Biot equations into two systems. One is stiff, and solved analytically, the other is nonstiff, and solved numerically by using a high-order staggered-grid finite-difference scheme. The time step is determined by the staggered finite-difference algorithm in solving the nonstiff equations, thus a coarse time step may be employed. Therefore, the computation efficiency and computational stability are improved greatly. Also a perfect by matched layer technology is used in the split method as absorbing boundary conditions. The numerical results are compared with the analytical results and those obtained from the conventional staggered-grid finite-difference method in a homogeneous model, respectively. They are in good agreement with each other. Finally, a slightly more complex model is investigated and compared with related equivalent model to illustrate the good performance of the staggered-grid finite-difference scheme in the partition method.展开更多
We investigated the effect of microscopic distribution modes of hydrates in porous sediments, and the saturation of hydrates and free gas on the elastic properties of saturated sediments. We simulated the propagation ...We investigated the effect of microscopic distribution modes of hydrates in porous sediments, and the saturation of hydrates and free gas on the elastic properties of saturated sediments. We simulated the propagation of seismic waves in gas hydrate-bearing sediments beneath the seafloor, and obtained the common receiver gathers of compressional waves（P-waves） and shear waves（S-waves）. The numerical results suggest that the interface between sediments containing gas hydrates and free gas produces a large-amplitude bottomsimulating reflector. The analysis of multicomponent common receiver data suggests that ocean-bottom seismometers receive the converted waves of upgoing P- and S-waves, which increases the complexity of the wavefield record.展开更多
Biot-flow and squirt-flow are the two most important fluid flow mechanisms in porous media containing fluids. Based on the BISQ （Biot-Squirt） model where the two mechanisms are treated simultaneously, the elastic wa...Biot-flow and squirt-flow are the two most important fluid flow mechanisms in porous media containing fluids. Based on the BISQ （Biot-Squirt） model where the two mechanisms are treated simultaneously, the elastic wave-field simulation in the porous medium is limited to two-dimensions and two-components （2D2C） or two-dimensions and three-components （2D3C）. There is no previous report on wave simulation in three- dimensions and three-components. Only through three dimensional numerical simulations can we have an overall understanding of wave field coupling relations and the spatial distribution characteristics between the solid and fluid phases in the dual-phase anisotropic medium. In this paper, based on the BISQ equation, we present elastic wave propagation in a three dimensional dual-phase anisotropic medium simulated by the staggered-grid high-order finite-difference method. We analyze the resulting wave fields and show that the results are an improvement.展开更多
In this paper, first we calculate finite-difference coefficients of implicit finite- difference methods （IFDM） for the first and second-order derivatives on normal grids and first- order derivatives on staggered gri...In this paper, first we calculate finite-difference coefficients of implicit finite- difference methods （IFDM） for the first and second-order derivatives on normal grids and first- order derivatives on staggered grids and find that small coefficients of high-order IFDMs exist. Dispersion analysis demonstrates that omitting these small coefficients can retain approximately the same order accuracy but greatly reduce computational costs. Then, we introduce a mirrorimage symmetric boundary condition to improve IFDMs accuracy and stability and adopt the hybrid absorbing boundary condition （ABC） to reduce unwanted reflections from the model boundary. Last, we give elastic wave modeling examples for homogeneous and heterogeneous models to demonstrate the advantages of the proposed scheme.展开更多
A three-dimensional （3D） predictor-corrector finite difference method for standing wave is developed. It is applied to solve the 3D nonlinear potential flow equa- tions with a free surface. The 3D irregular tank is ...A three-dimensional （3D） predictor-corrector finite difference method for standing wave is developed. It is applied to solve the 3D nonlinear potential flow equa- tions with a free surface. The 3D irregular tank is mapped onto a fixed cubic tank through the proper coordinate transform schemes. The cubic tank is distributed by the staggered meshgrid, and the staggered meshgrid is used to denote the variables of the flow field. The predictor-corrector finite difference method is given to develop the difference equa- tions of the dynamic boundary equation and kinematic boundary equation. Experimental results show that, using the finite difference method of the predictor-corrector scheme, the numerical solutions agree well with the published results. The wave profiles of the standing wave with different amplitudes and wave lengths are studied. The numerical solutions are also analyzed and presented graphically.展开更多
Based on the two-dimensional(2D)three-component first-order velocity-stress equation,the high order staggered mesh finite difference numerical simulation method was used to simulate the elastic and viscoelastic tilted...Based on the two-dimensional(2D)three-component first-order velocity-stress equation,the high order staggered mesh finite difference numerical simulation method was used to simulate the elastic and viscoelastic tilted transversely isotropic(TTI)media.The perfect matched layer(PML)absorption boundary condition was selected to eliminate the boundary effect.The results show that:(①)Under the condition of fixed elastic parameters of elastic TTI medium,when the polarization angle and azimuth are 60°and 45°respectively,the degree of shear wave splitting is significantly greater than the angle of 0°;②The influence of viscoelasticity on TTI medium is mainly reflected in the amplitude.If the quality factor decreases,the attenuation of the seismic wave amplitude increases,causing the waveform to become wider and distorted.If the quality factor increases,the viscoelastic medium becomes closer to elastic medium;③For TTI medium with different polarization angle and azimuth angle in the upper and lower layers,the shear wave can multiple splits at the interface of medium.The symmetry of seismograms is affected by the polarization angle and azimuth angle of TTI medium;④Viscoelasticity has a great influence on reflected wave,transmitted wave and converted wave in the low-velocity model.When the viscoelasticity is strong,the weaker waves may not be shown.展开更多
基金supported by the National Natural Science Foundation of China under grant Nos.40474012,40874020 and 40821062
文摘The pseudospectral method has been applied to the simulation of seismic wave propagation in 2-D global Earth model. When a whole Earth model is considered, the center of the Earth is included in the model and then singularity arises at the center of the Earth where r=0 since the 1/r term appears in the wave equations. In this paper, we extended the global seismic wavefield simulation algorithm for regular grid mesh to staggered grid configuration and developed a scheme to solve the numerical problems associated with the above singularity for a 2-D global Earth model defined on staggered grid using pseudospectral method. This scheme uses a coordinate transformation at the center of the model, in which the field variables at the center are calculated in Cartesian coordinates from the values on the grids around the center. It allows wave propagation through the center and hence the wavefield at the center can be stably calculated. Validity and accuracy of the scheme was tested by compared with the discrete wavenumber method. This scheme could also be suitable for other numerical methods or models parameterized in cylindrical or spherical coordinates when singularity arises at the center of the model.
基金Supported by the National Natural Science Foundation of China(Nos.51179195,51679248)the National Defense Foundation of China(No.513030203-02)
文摘Elastic waves in the seabed generated by low-frequency noise radiating from ships are known as ship seismic waves and can be used to detect and identify ships. To obtain the propagation characteristics of ship seismic waves, an algorithm for calculating Seismic waves at the seafloor is presented based on the staggered-grid finite difference method. The accuracy of the algorithm was tested by comparison with analytical solutions. Numerical simulation of seismic waves generated by a low-frequency point sotmd source in a typical shallow sea environment was carried out. Using various source frequencies and locations in the numerical simulation, we show that the seismic waves in the near field are composed mostly of transmitted S-waves and interface waves while transmitted P-waves are weak near the seafloor. However, in the far field, the wave components of the seismic wave are mainly normal modes and interface waves, with the latter being relatively strong in the waveforms, As the source frequency decreases, the normal modes become smaller and the interface waves dominate the time series of the seismic waves.
基金National Natural Science Foundation of China (5048003) and DAAD of Munich University, Germany.
文摘Strong ground motion of an earthquake is simulated by using both staggered grid finite difference method (FDM) and stochastic method, respectively. The acceleration time histories obtained from the both ways and their response spectra are compared. The result demonstrates that the former is adequate to simulate the low-frequency seismic wave; the latter is adequate to simulate the high-frequency seismic wave. Moreover, the result obtained from FDM can better reflect basin effects.
基金supported by the Natural Science Foundation of China(No.51676208)the Fundamental Research Funds for the Central Universities(No.18CX07012A and No.19CX05002A)support from the Major Program of the Natural Science Foundation of Shandong Province(No.ZR2019ZD11).
文摘A finite volume method based unstructured grid is presented to solve the two dimensional viscous and incompressible flow.The method is based on the pressure-correction concept and solved by using a semi-staggered grid technique.The computational procedure can handle cells of arbitrary shapes,although solutions presented in this paper were only involved with triangular and quadrilateral cells.The pressure or pressure-correction value was stored on the vertex of cells.The mass conservation equation was discretized on the dual cells surrounding the vertex of primary cells,while the velocity components and other scale variables were saved on the central of primary cells.Since the semi-staggered arrangement can’t guarantee a strong coupling relationship between pressure and velocity,thus a weak coupling relationship leads to the oscillations for pressure and velocity.In order to eliminate such an oscillation,a special interpolation scheme was used to construct the pressure-correction equation.Computational results of several viscous flow problems show good agreement with the analytical or numerical results in previous literature.This semi-staggered grid method can be applied to arbitrary shape elements,while it has the most efficiency for triangular cells.
文摘A new 3D finite-difference(FD) method of spatially asymmetric staggered grids was presented to simulate elastic wave propagation in topographic structures.The method approximated the first-order elastic wave equations by irregular grids finite difference operator with second-order time precise and fourth-order spatial precise.Additional introduced finite difference formula solved the asymmetric problem arisen in non-uniform staggered grid scheme.The method had no interpolation between the fine and coarse grids.All grids were computed at the same spatial iteration.Complicated geometrical structures like rough submarine interface,fault and nonplanar interfaces were treated with fine irregular grids.Theoretical analysis and numerical simulations show that this method saves considerable memory and computing time,at the same time,has satisfactory stability and accuracy.
基金Project supported by the National Natural Science Foundation of China (Nos. 51176204 and 51134006)
文摘In this paper, the cell face velocities in the discretization of the continu- ity equation, the momentum equation, and the scalar equation of a non-staggered grid system are calculated and discussed. Both the momentum interpolation and the linear interpolation are adopted to evaluate the coefficients in the discretized momentum and scalar equations. Their performances are compared. When the linear interpolation is used to calculate the coefficients, the mass residual term in the coefficients must be dropped to maintain the accuracy and convergence rate of the solution.
基金supported by the National Major Research Equipment Development Projects(No.ZDYZ2012-1-02-04)the National Natural Science Foundation of China(No.41474106)
文摘The numerical dispersion and computational cost are high for conventional Taylor series expansion staggered-grid finite-difference forward modeling owing to the high frequency of the wavelets and the large grid intervals. In this study, the cosine-modulated binomial window function （CMBWF）-based staggered-grid finite-difference method is proposed. Two new parameters, the modulated time and modulated range are used in the new window function and by adjusting these two parameters we obtain different characteristics of the main and side lobes of the amplitude response. Numerical dispersion analysis and elastic wavefield forward modeling suggests that the CMBWF method is more precise and less computationally costly than the conventional Taylor series expansion staggered-grid finite-difference method.
基金supported by the National Science and Technology Major Project of China(No.2011ZX05004-003)the National Basic Research Program of China(No.2013CB228602)the National High Tech Research Program of China(No.2013AA064202)
文摘Numerical simulations of a seismic wavefield are important to analyze seismic wave propagation. Elastic-wave equations are used in data simulation for modeling migration and imaging. In elastic wavefield numerical modeling, the rotated staggered-grid method （RSM） is a modification of the standard staggered-grid method （SSM）. The variable-order method is based on the method of variable-length spatial operators and wavefield propagation, and it calculates the real dispersion error by adapting different finite-difference orders to different velocities. In this study, the variable-order rotated staggered-grid method （VRSM） is developed after applying the variable-order method to RSM to solve the numerical dispersion problem of RSM in low-velocity regions and reduce the computation cost. Moreover, based on theoretical dispersion and the real dispersion error of wave propagation calculated with the wave separation method, the application of the original method is extended from acoustic to shear waves, and the calculation is modified from theoretical to time-varying values. A layered model and an overthrust model are used to demonstrate the applicability of VRSM. We also evaluate the order distribution, wave propagation, and computation time. The results suggest that the VRSM order distribution is reasonable and VRSM produces high-precision results with a minimal computation cost.
文摘A finite-volume formulation is proposed to solve the three-dimensional,non-hydrostatic Navier-Stokes equations on an unstructured,staggered,z-lever grid,with the goal of simulating non-hydrostatic processes in the free-surface flows.The advection and diffusion terms in the momentum equation are discretized explicitly with the Eulerian scheme,which has the attractive property of being conservative.An integral method of the top-layer pressure is applied to account for the full effects of non-hydrostatic pressure at the free-surface layer.It is shown that the results obtained with a small number of vertical layers（e.g.,2-3 layers） are in good agreements with experimental data or analytical solutions,demonstrating the efficiency and accuracy of the model in simulating a range of free-surface flow problems including wave motion and tide-induced motion.
基金supported by LASG,Institute of Atmospheric Physics,the ChineseAcademy of Sciences.
文摘This paper presents a new horizontal staggered grid(LE grid),which defines h at a gridpoint,and both u and v at the same mid-gridpoint along the x and y directions.A general method is used to deduce the dispersion relationships of describing inertia gravity waves on LE grid and Arakawa A―E grids,which are then compared with the analytical solution(AS)in re-solved or under-resolved cases,using two-order central difference or four-order compact differ-ence scheme from the frequency and group velocity.Results show that in both resolved and under-resolved cases,no matter whether two-order central difference or four-order compact dif-ference scheme is used,the frequency and group velocity discrete errors on LE grid in describing inertia gravity waves are smaller than those of Arakawa A―E grids.At the same time,it is only on LE or Arakawa grid C that the employment of a compact difference scheme of higher difference precision can improve their accuracy in describing inertia gravity waves.However,as for the other four grids(Arakawa A,B,D and E),when the difference precision increases,the accuracy of simulating inertia gravity waves decreases.
基金partially supported by National Key R&D Program of China(2021YFA0716902)the National Natural Science Foundation of China(42174156)the Fundamental Research Funds for the Central Universities,CHD(300102261107)。
文摘Reverse time migration and full waveform inversion involve the crosscorrelation of two wavefields,propagated in the forward-and reverse-time directions,respectively.As a result,the forward-propagated wavefield needs to be stored,and then accessed to compute the correlation with the backward-propagated wavefield.Boundary-value methods reconstruct the source wavefield using saved boundary wavefields and can significantly reduce the storage requirements.However,the existing boundary-value methods are based on the explicit finite-difference(FD)approximations of the spatial derivatives.Implicit FD methods exhibit greater accuracy and thus allow for a smaller operator length.We develop two(an accuracy-preserving and a memory-efficient)wavefield reconstruction schemes based on an implicit staggered-grid FD(SFD)operator.The former uses boundary wavefields at M layers of grid points and the spatial derivatives of wavefields at one layer of grid points to reconstruct the source wavefield for a(2M+2)th-order implicit SFD operator.The latter applies boundary wavefields at N layers of grid points,a linear combination of wavefields at M–N layers of grid points,and the spatial derivatives of wavefields at one layer of grid points to reconstruct the source wavefield(0≤N<M).The required memory of accuracy-preserving and memory-efficient schemes is(M+1)/M and(N+2)/M times,respectively,that of the explicit reconstruction scheme.Numerical results reveal that the accuracy-preserving scheme can achieve accurate reconstruction at the cost of storage.The memory-efficient scheme with N=2 can obtain plausible reconstructed wavefields and images,and the storage amount is 4/(M+1)of the accuracy-preserving scheme.
基金Supported by the National Natural Science Foundation of China(Nos. 41206043, 40930845)the Open Foundation of Key Laboratory of Marine Geology and Environment of Chinese Academy of Sciences(No. MGE2011KG07)+1 种基金the Knowledge Innovation Program of Chinese Academy of Sciences (No. KZCX2-YW-229)the National Basic Research Program of China (973 Program) (No. 2009CB219505)
文摘Seismic wave modeling is a cornerstone of geophysical data acquisition, processing, and interpretation, for which finite-difference methods are often applied. In this paper, we extend the velocity- pressure formulation of the acoustic wave equation to marine seismic modeling using the staggered-grid finite-difference method. The scheme is developed using a fourth-order spatial and a second-order temporal operator. Then, we define a stability coefficient (SC) and calculate its maximum value under the stability condition. Based on the dispersion relationship, we conduct a detailed dispersion analysis for submarine sediments in terms of the phase and group velocity over a range of angles, stability coefficients, and orders. We also compare the numerical solution with the exact solution for a P-wave line source in a homogeneous submarine model. Additionally, the numerical results determined by a Marmousi2 model with a rugged seafloor indicate that this method is sufficient for modeling complex submarine structures.
基金supported by the Special Project on High-performance Computing under the National Key R&D Program(No.2016YFB0200604)National Natural Science Foundation of China(11971502,11571385)Guangdong Natural Science Foundation(2017A030313017).
文摘In this paper,a fairly simple 3D immersed interface method based on the CG-Uzawa type method and the level set representation of the interface is employed for solving three-dimensional Stokes flow with singular forces along the interface.The method is to apply the Taylor’s expansions only along the normal direction and incorporate the jump conditions up to the second normal derivatives into the finite difference schemes.A second order geometric iteration algorithm is employed for computing orthogonal projections on the surface with third-order accuracy.The Stokes equations are discretized involving the correction terms on staggered grids and then solved by the conjugate gradient Uzawa type method.The major advantages of the present method are the special simplicity,the ability in handling the Dirichlet boundary conditions,and no need of the pressure boundary condition.The method can also preserve the volume conservation and the discrete divergence free condition very well.The numerical results show that the proposed method is second order accurate and efficient.
文摘A depth averaged nonlinear k ε model for turbulent flows in complex geometries has been developed in a boundary fitted coordinate system. The SIMPLEC procedure is used to develop an economical discrete method for staggered grids to analyze flows in a 90° bend. This paper describes how to change a program in rectangular coordinate into a boundary fitted coordinate. The results compare well with experimental data for flow in a meandering channel showing the efficiency of the model and the discrete method.
基金Project supported by the ＂100 Talents Project＂ of the Chinese Academy of Sciences and the Major Program of the National Natural Science Foundation of China （Grant No 10534040）.
文摘In modelling elastic wave propagation in a porous medium, when the ratio between the fluid viscosity and the medium permeability is comparatively large, the stiffness problem of Blot＇s poroelastic equations will be encountered. In the paper, a partition method is developed to solve the stiffness problem with a staggered high-order finite-difference. The method splits the Biot equations into two systems. One is stiff, and solved analytically, the other is nonstiff, and solved numerically by using a high-order staggered-grid finite-difference scheme. The time step is determined by the staggered finite-difference algorithm in solving the nonstiff equations, thus a coarse time step may be employed. Therefore, the computation efficiency and computational stability are improved greatly. Also a perfect by matched layer technology is used in the split method as absorbing boundary conditions. The numerical results are compared with the analytical results and those obtained from the conventional staggered-grid finite-difference method in a homogeneous model, respectively. They are in good agreement with each other. Finally, a slightly more complex model is investigated and compared with related equivalent model to illustrate the good performance of the staggered-grid finite-difference scheme in the partition method.
基金supported by the National Natural Science Foundation of China(No.41174087,41204089)the National Oil and Gas Major Project(No.2011ZX05005-005)
文摘We investigated the effect of microscopic distribution modes of hydrates in porous sediments, and the saturation of hydrates and free gas on the elastic properties of saturated sediments. We simulated the propagation of seismic waves in gas hydrate-bearing sediments beneath the seafloor, and obtained the common receiver gathers of compressional waves（P-waves） and shear waves（S-waves）. The numerical results suggest that the interface between sediments containing gas hydrates and free gas produces a large-amplitude bottomsimulating reflector. The analysis of multicomponent common receiver data suggests that ocean-bottom seismometers receive the converted waves of upgoing P- and S-waves, which increases the complexity of the wavefield record.
基金National Natural Science Foundation (Project number 40604013).
文摘Biot-flow and squirt-flow are the two most important fluid flow mechanisms in porous media containing fluids. Based on the BISQ （Biot-Squirt） model where the two mechanisms are treated simultaneously, the elastic wave-field simulation in the porous medium is limited to two-dimensions and two-components （2D2C） or two-dimensions and three-components （2D3C）. There is no previous report on wave simulation in three- dimensions and three-components. Only through three dimensional numerical simulations can we have an overall understanding of wave field coupling relations and the spatial distribution characteristics between the solid and fluid phases in the dual-phase anisotropic medium. In this paper, based on the BISQ equation, we present elastic wave propagation in a three dimensional dual-phase anisotropic medium simulated by the staggered-grid high-order finite-difference method. We analyze the resulting wave fields and show that the results are an improvement.
基金supported by the National Natural Science Foundation of China(NSFC)(Grant No. 41074100)the Program for New Century Excellent Talents in University of Ministry of Education of China(Grant No. NCET-10-0812)
文摘In this paper, first we calculate finite-difference coefficients of implicit finite- difference methods （IFDM） for the first and second-order derivatives on normal grids and first- order derivatives on staggered grids and find that small coefficients of high-order IFDMs exist. Dispersion analysis demonstrates that omitting these small coefficients can retain approximately the same order accuracy but greatly reduce computational costs. Then, we introduce a mirrorimage symmetric boundary condition to improve IFDMs accuracy and stability and adopt the hybrid absorbing boundary condition （ABC） to reduce unwanted reflections from the model boundary. Last, we give elastic wave modeling examples for homogeneous and heterogeneous models to demonstrate the advantages of the proposed scheme.
基金supported by the Yunnan Provincial Applied Basic Research Program of China(No. KKSY201207019)
文摘A three-dimensional （3D） predictor-corrector finite difference method for standing wave is developed. It is applied to solve the 3D nonlinear potential flow equa- tions with a free surface. The 3D irregular tank is mapped onto a fixed cubic tank through the proper coordinate transform schemes. The cubic tank is distributed by the staggered meshgrid, and the staggered meshgrid is used to denote the variables of the flow field. The predictor-corrector finite difference method is given to develop the difference equa- tions of the dynamic boundary equation and kinematic boundary equation. Experimental results show that, using the finite difference method of the predictor-corrector scheme, the numerical solutions agree well with the published results. The wave profiles of the standing wave with different amplitudes and wave lengths are studied. The numerical solutions are also analyzed and presented graphically.
基金the National Natural Science Foundation of China(Nos.41974048,41574078,41604039,41604102)the Guangxi Natural Science Foundation of China(Nos.2018GXNSFAA138059,2016GXNSFBA380082 and 2018GXNSFBA050005)+1 种基金Guangxi Science and Technology Base and Talent Project(Gui Kc AD19110057)Guangxi High School Junior Teachers Foundation Funding for capacity improvement projects(No.2019KY0264).
文摘Based on the two-dimensional(2D)three-component first-order velocity-stress equation,the high order staggered mesh finite difference numerical simulation method was used to simulate the elastic and viscoelastic tilted transversely isotropic(TTI)media.The perfect matched layer(PML)absorption boundary condition was selected to eliminate the boundary effect.The results show that:(①)Under the condition of fixed elastic parameters of elastic TTI medium,when the polarization angle and azimuth are 60°and 45°respectively,the degree of shear wave splitting is significantly greater than the angle of 0°;②The influence of viscoelasticity on TTI medium is mainly reflected in the amplitude.If the quality factor decreases,the attenuation of the seismic wave amplitude increases,causing the waveform to become wider and distorted.If the quality factor increases,the viscoelastic medium becomes closer to elastic medium;③For TTI medium with different polarization angle and azimuth angle in the upper and lower layers,the shear wave can multiple splits at the interface of medium.The symmetry of seismograms is affected by the polarization angle and azimuth angle of TTI medium;④Viscoelasticity has a great influence on reflected wave,transmitted wave and converted wave in the low-velocity model.When the viscoelasticity is strong,the weaker waves may not be shown.
