Equivalent staggered-grid(ESG) as a new family of schemes has been utilized in seismic modeling,imaging,and inversion.Traditionally,the Taylor series expansion is often applied to calculate finite-difference(FD) coeff...Equivalent staggered-grid(ESG) as a new family of schemes has been utilized in seismic modeling,imaging,and inversion.Traditionally,the Taylor series expansion is often applied to calculate finite-difference(FD) coefficients on spatial derivatives,but the simulation results suffer serious numerical dispersion on a large frequency zone.We develop an optimized equivalent staggered-grid(OESG) FD method that can simultaneously suppress temporal and spatial dispersion for solving the second-order system of the 3 D elastic wave equation.On the one hand,we consider the coupling relations between wave speeds and spatial derivatives in the elastic wave equation and give three sets of FD coefficients with respect to the P-wave,S-wave,and converted-wave(C-wave) terms.On the other hand,a novel plane wave solution for the 3 D elastic wave equation is derived from the matrix decomposition method to construct the time-space dispersion relations.FD coefficients of the OESG method can be acquired by solving the new dispersion equations based on the Newton iteration method.Finally,we construct a new objective function to analyze P-wave,S-wave,and C-wave dispersion concerning frequencies.The dispersion analyses show that the presented method produces less modeling errors than the traditional ESG method.The synthetic examples demonstrate the effectiveness and superiority of the presented method.展开更多
The nonlinear interactions of waves with a double-peaked power spectrum have been studied in shallow water. The starting point is the prototypical equation for nonlinear unidirectional waves in shallow water, i.e. the...The nonlinear interactions of waves with a double-peaked power spectrum have been studied in shallow water. The starting point is the prototypical equation for nonlinear unidirectional waves in shallow water, i.e. the Korteweg de Vries equation. By means of a multiple-scale technique two defocusing coupled Nonlinear SchrCMinger equations are derived. It is found analytically that plane wave solutions of such a system are unstable for small perturbations, showing that the existence of a new energy exchange mechanism which can influence the behavior of ocean waves in shallow water.展开更多
In this paper we are going to derive two numerical methods for solving the coupled nonlinear Schrodinger-Boussinesq equation. The first method is a nonlinear implicit scheme of second order accuracy in both directions...In this paper we are going to derive two numerical methods for solving the coupled nonlinear Schrodinger-Boussinesq equation. The first method is a nonlinear implicit scheme of second order accuracy in both directions time and space;the scheme is unconditionally stable. The second scheme is a nonlinear implicit scheme of second order accuracy in time and fourth order accuracy in space direction. A generalized method is also derived where the previous schemes can be obtained by some special values of l. The proposed methods will produced a coupled nonlinear tridiagonal system which can be solved by fixed point method. The exact solutions and the conserved quantities for two different tests are used to display the robustness of the proposed schemes.展开更多
文摘Equivalent staggered-grid(ESG) as a new family of schemes has been utilized in seismic modeling,imaging,and inversion.Traditionally,the Taylor series expansion is often applied to calculate finite-difference(FD) coefficients on spatial derivatives,but the simulation results suffer serious numerical dispersion on a large frequency zone.We develop an optimized equivalent staggered-grid(OESG) FD method that can simultaneously suppress temporal and spatial dispersion for solving the second-order system of the 3 D elastic wave equation.On the one hand,we consider the coupling relations between wave speeds and spatial derivatives in the elastic wave equation and give three sets of FD coefficients with respect to the P-wave,S-wave,and converted-wave(C-wave) terms.On the other hand,a novel plane wave solution for the 3 D elastic wave equation is derived from the matrix decomposition method to construct the time-space dispersion relations.FD coefficients of the OESG method can be acquired by solving the new dispersion equations based on the Newton iteration method.Finally,we construct a new objective function to analyze P-wave,S-wave,and C-wave dispersion concerning frequencies.The dispersion analyses show that the presented method produces less modeling errors than the traditional ESG method.The synthetic examples demonstrate the effectiveness and superiority of the presented method.
基金The paper was financially supporrted by the National Natural Science Foundation of China (Grant No40476062)Foundation of Hebei University of Science and Technology
文摘The nonlinear interactions of waves with a double-peaked power spectrum have been studied in shallow water. The starting point is the prototypical equation for nonlinear unidirectional waves in shallow water, i.e. the Korteweg de Vries equation. By means of a multiple-scale technique two defocusing coupled Nonlinear SchrCMinger equations are derived. It is found analytically that plane wave solutions of such a system are unstable for small perturbations, showing that the existence of a new energy exchange mechanism which can influence the behavior of ocean waves in shallow water.
文摘In this paper we are going to derive two numerical methods for solving the coupled nonlinear Schrodinger-Boussinesq equation. The first method is a nonlinear implicit scheme of second order accuracy in both directions time and space;the scheme is unconditionally stable. The second scheme is a nonlinear implicit scheme of second order accuracy in time and fourth order accuracy in space direction. A generalized method is also derived where the previous schemes can be obtained by some special values of l. The proposed methods will produced a coupled nonlinear tridiagonal system which can be solved by fixed point method. The exact solutions and the conserved quantities for two different tests are used to display the robustness of the proposed schemes.
