This paper presents an analogical study between electromagnetic and elastic wave fields,with a one-to-one correspondence principle established regarding the basic wave equations,the physical quantities and the differe...This paper presents an analogical study between electromagnetic and elastic wave fields,with a one-to-one correspondence principle established regarding the basic wave equations,the physical quantities and the differential operations.Using the electromagnetic-to-elastic substitution,the analogous relations of the conservation laws of energy and momentum are investigated between these two physical fields.Moreover,the energy-based and momentum-based reciprocity theorems for an elastic wave are also derived in the time-harmonic state,which describe the interaction between two elastic wave systems from the perspectives of energy and momentum,respectively.The theoretical results obtained in this analysis can not only improve our understanding of the similarities of these two linear systems,but also find potential applications in relevant fields such as medical imaging,non-destructive evaluation,acoustic microscopy,seismology and exploratory geophysics.展开更多
Based on an algebraically Rossby solitary waves evolution model,namely an extended(2+1)-dimensional Boussinesq equation,we firstly introduced a special transformation and utilized the Hirota method,which enable us to ...Based on an algebraically Rossby solitary waves evolution model,namely an extended(2+1)-dimensional Boussinesq equation,we firstly introduced a special transformation and utilized the Hirota method,which enable us to obtain multi-complexiton solutions and explore the interaction among the solutions.These wave functions are then employed to infer the influence of background flow on the propagation of Rossby waves,as well as the characteristics of propagation in multi-wave running processes.Additionally,we generated stereogram drawings and projection figures to visually represent these solutions.The dynamical behavior of these solutions is thoroughly examined through analytical and graphical analyses.Furthermore,we investigated the influence of the generalized beta effect and the Coriolis parameter on the evolution of Rossby waves.展开更多
In the generalized continuum mechanics(GCM)theory framework,asymmetric wave equations encompass the characteristic scale parameters of the medium,accounting for microstructure interactions.This study integrates two th...In the generalized continuum mechanics(GCM)theory framework,asymmetric wave equations encompass the characteristic scale parameters of the medium,accounting for microstructure interactions.This study integrates two theoretical branches of the GCM,the modified couple stress theory(M-CST)and the one-parameter second-strain-gradient theory,to form a novel asymmetric wave equation in a unified framework.Numerical modeling of the asymmetric wave equation in a unified framework accurately describes subsurface structures with vital implications for subsequent seismic wave inversion and imaging endeavors.However,employing finite-difference(FD)methods for numerical modeling may introduce numerical dispersion,adversely affecting the accuracy of numerical modeling.The design of an optimal FD operator is crucial for enhancing the accuracy of numerical modeling and emphasizing the scale effects.Therefore,this study devises a hybrid scheme called the dung beetle optimization(DBO)algorithm with a simulated annealing(SA)algorithm,denoted as the SA-based hybrid DBO(SDBO)algorithm.An FD operator optimization method under the SDBO algorithm was developed and applied to the numerical modeling of asymmetric wave equations in a unified framework.Integrating the DBO and SA algorithms mitigates the risk of convergence to a local extreme.The numerical dispersion outcomes underscore that the proposed SDBO algorithm yields FD operators with precision errors constrained to 0.5‱while encompassing a broader spectrum coverage.This result confirms the efficacy of the SDBO algorithm.Ultimately,the numerical modeling results demonstrate that the new FD method based on the SDBO algorithm effectively suppresses numerical dispersion and enhances the accuracy of elastic wave numerical modeling,thereby accentuating scale effects.This result is significant for extracting wavefield perturbations induced by complex microstructures in the medium and the analysis of scale effects.展开更多
In this paper,we establish global classical solutions of semilinear wave equations with small compact supported initial data posed on the product space R^(3)×T.The semilinear nonlinearity is assumed to be of the ...In this paper,we establish global classical solutions of semilinear wave equations with small compact supported initial data posed on the product space R^(3)×T.The semilinear nonlinearity is assumed to be of the cubic form.The main ingredient here is the establishment of the L^(2)-L^(∞)decay estimates and the energy estimates for the linear problem,which are adapted to the wave equation on the product space.The proof is based on the Fourier mode decomposition of the solution with respect to the periodic direction,the scaling technique,and the combination of the decay estimates and the energy estimates.展开更多
This paper studies the(2+1)-dimensional Hirota-Satsuma-Ito equation.Based on an associated Hirota bilinear form,lump-type solution,two types of interaction solutions,and breather wave solution of the(2+1)-dimensional ...This paper studies the(2+1)-dimensional Hirota-Satsuma-Ito equation.Based on an associated Hirota bilinear form,lump-type solution,two types of interaction solutions,and breather wave solution of the(2+1)-dimensional Hirota-Satsuma-Ito equation are obtained,which are all related to the seed solution of the equation.It is interesting that the rogue wave is aroused by the interaction between one-lump soliton and a pair of resonance stripe solitons,and the fusion and fission phenomena are also found in the interaction between lump solitons and one-stripe soliton.Furthermore,the breather wave solution is also obtained by reducing the two-soliton solutions.The trajectory and period of the one-order breather wave are analyzed.The corresponding dynamical characteristics are demonstrated by the graphs.展开更多
For the (2 + 1)-dimensional nonlinear dispersive Boussinesq equation, by using the bifurcation theory of planar dynamical systems to study its corresponding traveling wave system, the bifurcations and phase portraits ...For the (2 + 1)-dimensional nonlinear dispersive Boussinesq equation, by using the bifurcation theory of planar dynamical systems to study its corresponding traveling wave system, the bifurcations and phase portraits of the regular system are obtained. Under different parametric conditions, various sufficient conditions to guarantee the existence of analytical and non-analytical solutions of the singular system are given by using singular traveling wave theory. For certain special cases, some explicit and exact parametric representations of traveling wave solutions are derived such as analytical periodic waves and non-analytical periodic cusp waves. Further, two-dimensional wave plots of analytical periodic solutions and non-analytical periodic cusp wave solutions are drawn to visualize the dynamics of the equation.展开更多
In this paper, we study the propagation and its failure to propagate (pinning) of a travelling wave in a Nagumo type equation, an equation that describes impulse propagation in nerve axons that also models population ...In this paper, we study the propagation and its failure to propagate (pinning) of a travelling wave in a Nagumo type equation, an equation that describes impulse propagation in nerve axons that also models population growth with Allee effect. An analytical solution is derived for the traveling wave and the work is extended to a discrete formulation with a piecewise linear reaction function. We propose an operator splitting numerical scheme to solve the equation and demonstrate that the wave either propagates or gets pinned based on how the spatial mesh is chosen.展开更多
A periodically homoclinic solution and some rogue wave solutions of (1+1)-dimensional Boussinesq equation are obtained via the limit behavior of parameters and different polynomial functions. Besides, the mathematics ...A periodically homoclinic solution and some rogue wave solutions of (1+1)-dimensional Boussinesq equation are obtained via the limit behavior of parameters and different polynomial functions. Besides, the mathematics reasons for different spatiotemporal structures of rogue waves are analyzed using the extreme value theory of the two-variables function. The diversity of spatiotemporal structures not only depends on the disturbance parameter u0 </sub>but also has a relationship with the other parameters c<sub>0</sub>, α, β.展开更多
In this paper, we consider the generalized Korteweg-de-Vries (KdV) equations which are remarkable models of the water waves mechanics, the shallow water waves, the quantum mechanics, the ion acoustic waves in plasma, ...In this paper, we consider the generalized Korteweg-de-Vries (KdV) equations which are remarkable models of the water waves mechanics, the shallow water waves, the quantum mechanics, the ion acoustic waves in plasma, the electro-hydro-dynamical model for local electric field, signal processing waves through optical fibers, etc. We determine the useful and further general exact traveling wave solutions of the above mentioned NLDEs by applying the exp(−τ(ξ))-expansion method by aid of traveling wave transformations. Furthermore, we explain the physical significance of the obtained solutions of its definite values of the involved parameters with graphic representations in order to know the physical phenomena. Finally, we show that the exp(−τ(ξ))-expansion method is convenient, powerful, straightforward and provide more general solutions and can be helping to examine vast amount of travelling wave solutions to the other different kinds of NLDEs.展开更多
In this paper, we study the existence of standing waves for the nonlinear Schrödinger equation with combined power-type nonlinearities and a partial harmonic potential. In the L<sup>2</sup>-supercriti...In this paper, we study the existence of standing waves for the nonlinear Schrödinger equation with combined power-type nonlinearities and a partial harmonic potential. In the L<sup>2</sup>-supercritical case, we obtain the existence and stability of standing waves. Our results are complements to the results of Carles and Il’yasov’s artical, where orbital stability of standing waves have been studied for the 2D Schrödinger equation with combined nonlinearities and harmonic potential.展开更多
By employing the Hirota’s bilinear method and different test functions, the breather solutions of HSI equation with different structures are obtained based on symbolic calculation with perturbation parameters. Some n...By employing the Hirota’s bilinear method and different test functions, the breather solutions of HSI equation with different structures are obtained based on symbolic calculation with perturbation parameters. Some new lump solitons are found in the process of studying the degradation behavior of breather solutions. The interaction between lump solution and soliton solution is constructed in the form of lump solution, and the motion trajectory of lump is obtained. In addition, the theorem of lump solitons and N-solitons superposition is given and proved. The superposition formula of lump is derived from the theorem, and its spatial evolution behavior is given.展开更多
Traditional pre-stack depth migration can only provide subsurface structural information. However, simple structure information is insufficient for petroleum exploration which also needs amplitude information proporti...Traditional pre-stack depth migration can only provide subsurface structural information. However, simple structure information is insufficient for petroleum exploration which also needs amplitude information proportional to reflection coefficients. In recent years, pre-stack depth migration algorithms which preserve amplitudes and based on the one- way wave equation have been developed. Using the method in the shot domain requires a deconvolution imaging condition which produces some instability in areas with complicated structure and dramatic lateral variation in velocity. Depth migration with preserved amplitude based on the angle domain can overcome the instability of the one-way wave migration imaging condition with preserved amplitude. It can also offer provide velocity analysis in the angle domain of common imaging point gathers. In this paper, based on the foundation of the one-way wave continuation operator with preserved amplitude, we realized the preserved amplitude prestack depth migration in the angle domain. Models and real data validate the accuracy of the method.展开更多
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.展开更多
Based on arbitrarily wide-angle wave equations,a reverse-time propagation scheme is developed by substituting the partial derivatives of depth and time with central differences. The partial derivative of horizontal di...Based on arbitrarily wide-angle wave equations,a reverse-time propagation scheme is developed by substituting the partial derivatives of depth and time with central differences. The partial derivative of horizontal direction is replaced with high order difference. The imaging condition is computed by solving the eikonal equations. On the basis of above techniques,a prestack reverse-time depth migration algorithm is developed. The processing exam-ples of synthetic data show that the method can remove unwanted internal reflections and decrease the migration noise. The method also has the advantage of fidelity and is applicable of dip angle reflector imaging.展开更多
The perfectly matched layer (PML) is a highly efficient absorbing boundary condition used for the numerical modeling of seismic wave equation. The article focuses on the application of this technique to finite-eleme...The perfectly matched layer (PML) is a highly efficient absorbing boundary condition used for the numerical modeling of seismic wave equation. The article focuses on the application of this technique to finite-element time-domain numerical modeling of elastic wave equation. However, the finite-element time-domain scheme is based on the second- order wave equation in displacement formulation. Thus, the first-order PML in velocity-stress formulation cannot be directly applied to this scheme. In this article, we derive the finite- element matrix equations of second-order PML in displacement formulation, and accomplish the implementation of PML in finite-element time-domain modeling of elastic wave equation. The PML has an approximate zero reflection coefficients for bulk and surface waves in the finite-element modeling of P-SV and SH wave propagation in the 2D homogeneous elastic media. The numerical experiments using a two-layer model with irregular topography validate the efficiency of PML in the modeling of seismic wave propagation in geological models with complex structures and heterogeneous media.展开更多
In this paper the decay of global solutions to some nonlinear dissipative wave equations are discussed, which based on the method of prior estimate technique and a differenece inequality.
This paper deals with the initial-boundary value mixed problems for nonlinear wave equations. By introducing the 'blowing-up facts K(u,u_i)', We may discuss the blowing up behaviours of solutions in finite tim...This paper deals with the initial-boundary value mixed problems for nonlinear wave equations. By introducing the 'blowing-up facts K(u,u_i)', We may discuss the blowing up behaviours of solutions in finite time to the mixed problems with respect to Neumann boundary and Dirichlet boundary for various nonlinear conditions and initial value conditions which usually meet.展开更多
The existence and the nonexistence,the uniqueness and the energy decay estimate of solution for the fourth-order nonlinear wave equation utt+αΔ2 u-bΔut-βΔu+ut|ut|^r+g(u)=0 in Ω×(0,∞) are studied w...The existence and the nonexistence,the uniqueness and the energy decay estimate of solution for the fourth-order nonlinear wave equation utt+αΔ2 u-bΔut-βΔu+ut|ut|^r+g(u)=0 in Ω×(0,∞) are studied with the boundary condition u=(u)/(υ)=0 onΩ and the initial condition u(x,0)=u0(x),ut(x,0)=u1(x,0) in bounded domain ΩR^n ,n≥1.The energy decay rate of the global solution is estimated by the multiplier method.The blow-up result of the solution in finite time is established by the ideal of a potential well theory,and the existence of the solution is gotten by the Galekin approximation method.展开更多
The Finite Difference (FD) method is an important method for seismic numerical simulations. It helps us understand regular patterns in seismic wave propagation, analyze seismic attributes, and interpret seismic data...The Finite Difference (FD) method is an important method for seismic numerical simulations. It helps us understand regular patterns in seismic wave propagation, analyze seismic attributes, and interpret seismic data. However, because of its discretization, the FD method is only stable under certain conditions. The Arbitrary Difference Precise Integration (ADPI) method is based on the FD method and adopts an integration scheme in the time domain and an arbitrary difference scheme in the space domain. Therefore, the ADPI method is a semi-analytical method. In this paper, we deduce the formula for the ADPI method based on the 3D elastic equation and improve its stability. In forward modeling cases, the ADPI method was implemented in 2D and 3D elastic wave equation forward modeling. Results show that the travel time of the reflected seismic wave is accurate. Compared with the acoustic wave field, the elastic wave field contains more wave types, including PS- and PP- reflected waves, transmitted waves, and diffracted waves, which is important to interpretation of seismic data. The method can be easily applied to elastic wave equation numerical simulations for eoloical models.展开更多
基金funded by the National Natural Science Foundation of China(Grant No.12404507)the Natural Science Research of Jiangsu Higher Education Institutions of China(Grant No.24KJB140013)the Scientific Startup Foundation of Nanjing Normal University(Grant No.184080H201B49).
文摘This paper presents an analogical study between electromagnetic and elastic wave fields,with a one-to-one correspondence principle established regarding the basic wave equations,the physical quantities and the differential operations.Using the electromagnetic-to-elastic substitution,the analogous relations of the conservation laws of energy and momentum are investigated between these two physical fields.Moreover,the energy-based and momentum-based reciprocity theorems for an elastic wave are also derived in the time-harmonic state,which describe the interaction between two elastic wave systems from the perspectives of energy and momentum,respectively.The theoretical results obtained in this analysis can not only improve our understanding of the similarities of these two linear systems,but also find potential applications in relevant fields such as medical imaging,non-destructive evaluation,acoustic microscopy,seismology and exploratory geophysics.
基金Supported by the National Natural Science Foundation of China(No.32360249)the Natural Science Foundation of Inner Mongolia Autonomous Region of China(No.2022QN01003)+2 种基金the University Scientific Research Project of Inner Mongolia Autonomous Region of China(No.NJZY22484)the Scientific Research Improvement Project of Youth Teachers of Inner Mongolia Autonomous Region of China(No.BR230161)the Inner Mongolia Agricultural University Basic Discipline Scientific Research Launch Fund(No.JC2020003)。
文摘Based on an algebraically Rossby solitary waves evolution model,namely an extended(2+1)-dimensional Boussinesq equation,we firstly introduced a special transformation and utilized the Hirota method,which enable us to obtain multi-complexiton solutions and explore the interaction among the solutions.These wave functions are then employed to infer the influence of background flow on the propagation of Rossby waves,as well as the characteristics of propagation in multi-wave running processes.Additionally,we generated stereogram drawings and projection figures to visually represent these solutions.The dynamical behavior of these solutions is thoroughly examined through analytical and graphical analyses.Furthermore,we investigated the influence of the generalized beta effect and the Coriolis parameter on the evolution of Rossby waves.
基金supported by project XJZ2023050044,A2309002 and XJZ2023070052.
文摘In the generalized continuum mechanics(GCM)theory framework,asymmetric wave equations encompass the characteristic scale parameters of the medium,accounting for microstructure interactions.This study integrates two theoretical branches of the GCM,the modified couple stress theory(M-CST)and the one-parameter second-strain-gradient theory,to form a novel asymmetric wave equation in a unified framework.Numerical modeling of the asymmetric wave equation in a unified framework accurately describes subsurface structures with vital implications for subsequent seismic wave inversion and imaging endeavors.However,employing finite-difference(FD)methods for numerical modeling may introduce numerical dispersion,adversely affecting the accuracy of numerical modeling.The design of an optimal FD operator is crucial for enhancing the accuracy of numerical modeling and emphasizing the scale effects.Therefore,this study devises a hybrid scheme called the dung beetle optimization(DBO)algorithm with a simulated annealing(SA)algorithm,denoted as the SA-based hybrid DBO(SDBO)algorithm.An FD operator optimization method under the SDBO algorithm was developed and applied to the numerical modeling of asymmetric wave equations in a unified framework.Integrating the DBO and SA algorithms mitigates the risk of convergence to a local extreme.The numerical dispersion outcomes underscore that the proposed SDBO algorithm yields FD operators with precision errors constrained to 0.5‱while encompassing a broader spectrum coverage.This result confirms the efficacy of the SDBO algorithm.Ultimately,the numerical modeling results demonstrate that the new FD method based on the SDBO algorithm effectively suppresses numerical dispersion and enhances the accuracy of elastic wave numerical modeling,thereby accentuating scale effects.This result is significant for extracting wavefield perturbations induced by complex microstructures in the medium and the analysis of scale effects.
文摘In this paper,we establish global classical solutions of semilinear wave equations with small compact supported initial data posed on the product space R^(3)×T.The semilinear nonlinearity is assumed to be of the cubic form.The main ingredient here is the establishment of the L^(2)-L^(∞)decay estimates and the energy estimates for the linear problem,which are adapted to the wave equation on the product space.The proof is based on the Fourier mode decomposition of the solution with respect to the periodic direction,the scaling technique,and the combination of the decay estimates and the energy estimates.
基金Project supported by the National Natural Science Foundation of China (Grant Nos.12275172 and 11905124)。
文摘This paper studies the(2+1)-dimensional Hirota-Satsuma-Ito equation.Based on an associated Hirota bilinear form,lump-type solution,two types of interaction solutions,and breather wave solution of the(2+1)-dimensional Hirota-Satsuma-Ito equation are obtained,which are all related to the seed solution of the equation.It is interesting that the rogue wave is aroused by the interaction between one-lump soliton and a pair of resonance stripe solitons,and the fusion and fission phenomena are also found in the interaction between lump solitons and one-stripe soliton.Furthermore,the breather wave solution is also obtained by reducing the two-soliton solutions.The trajectory and period of the one-order breather wave are analyzed.The corresponding dynamical characteristics are demonstrated by the graphs.
文摘For the (2 + 1)-dimensional nonlinear dispersive Boussinesq equation, by using the bifurcation theory of planar dynamical systems to study its corresponding traveling wave system, the bifurcations and phase portraits of the regular system are obtained. Under different parametric conditions, various sufficient conditions to guarantee the existence of analytical and non-analytical solutions of the singular system are given by using singular traveling wave theory. For certain special cases, some explicit and exact parametric representations of traveling wave solutions are derived such as analytical periodic waves and non-analytical periodic cusp waves. Further, two-dimensional wave plots of analytical periodic solutions and non-analytical periodic cusp wave solutions are drawn to visualize the dynamics of the equation.
文摘In this paper, we study the propagation and its failure to propagate (pinning) of a travelling wave in a Nagumo type equation, an equation that describes impulse propagation in nerve axons that also models population growth with Allee effect. An analytical solution is derived for the traveling wave and the work is extended to a discrete formulation with a piecewise linear reaction function. We propose an operator splitting numerical scheme to solve the equation and demonstrate that the wave either propagates or gets pinned based on how the spatial mesh is chosen.
文摘A periodically homoclinic solution and some rogue wave solutions of (1+1)-dimensional Boussinesq equation are obtained via the limit behavior of parameters and different polynomial functions. Besides, the mathematics reasons for different spatiotemporal structures of rogue waves are analyzed using the extreme value theory of the two-variables function. The diversity of spatiotemporal structures not only depends on the disturbance parameter u0 </sub>but also has a relationship with the other parameters c<sub>0</sub>, α, β.
文摘In this paper, we consider the generalized Korteweg-de-Vries (KdV) equations which are remarkable models of the water waves mechanics, the shallow water waves, the quantum mechanics, the ion acoustic waves in plasma, the electro-hydro-dynamical model for local electric field, signal processing waves through optical fibers, etc. We determine the useful and further general exact traveling wave solutions of the above mentioned NLDEs by applying the exp(−τ(ξ))-expansion method by aid of traveling wave transformations. Furthermore, we explain the physical significance of the obtained solutions of its definite values of the involved parameters with graphic representations in order to know the physical phenomena. Finally, we show that the exp(−τ(ξ))-expansion method is convenient, powerful, straightforward and provide more general solutions and can be helping to examine vast amount of travelling wave solutions to the other different kinds of NLDEs.
文摘In this paper, we study the existence of standing waves for the nonlinear Schrödinger equation with combined power-type nonlinearities and a partial harmonic potential. In the L<sup>2</sup>-supercritical case, we obtain the existence and stability of standing waves. Our results are complements to the results of Carles and Il’yasov’s artical, where orbital stability of standing waves have been studied for the 2D Schrödinger equation with combined nonlinearities and harmonic potential.
文摘By employing the Hirota’s bilinear method and different test functions, the breather solutions of HSI equation with different structures are obtained based on symbolic calculation with perturbation parameters. Some new lump solitons are found in the process of studying the degradation behavior of breather solutions. The interaction between lump solution and soliton solution is constructed in the form of lump solution, and the motion trajectory of lump is obtained. In addition, the theorem of lump solitons and N-solitons superposition is given and proved. The superposition formula of lump is derived from the theorem, and its spatial evolution behavior is given.
基金supported by the National 863 Program(Grant No.2006AA06Z206)the National 973 Program(Grant No.2007CB209605)CNPC geophysical laboratories and Ph.D innovative funding in China University of Petroleum(East China)
文摘Traditional pre-stack depth migration can only provide subsurface structural information. However, simple structure information is insufficient for petroleum exploration which also needs amplitude information proportional to reflection coefficients. In recent years, pre-stack depth migration algorithms which preserve amplitudes and based on the one- way wave equation have been developed. Using the method in the shot domain requires a deconvolution imaging condition which produces some instability in areas with complicated structure and dramatic lateral variation in velocity. Depth migration with preserved amplitude based on the angle domain can overcome the instability of the one-way wave migration imaging condition with preserved amplitude. It can also offer provide velocity analysis in the angle domain of common imaging point gathers. In this paper, based on the foundation of the one-way wave continuation operator with preserved amplitude, we realized the preserved amplitude prestack depth migration in the angle domain. Models and real data validate the accuracy of the method.
基金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.
文摘Based on arbitrarily wide-angle wave equations,a reverse-time propagation scheme is developed by substituting the partial derivatives of depth and time with central differences. The partial derivative of horizontal direction is replaced with high order difference. The imaging condition is computed by solving the eikonal equations. On the basis of above techniques,a prestack reverse-time depth migration algorithm is developed. The processing exam-ples of synthetic data show that the method can remove unwanted internal reflections and decrease the migration noise. The method also has the advantage of fidelity and is applicable of dip angle reflector imaging.
基金sponsored by the National Natural Science Foundation of China Research(Grant No.41274138)the Science Foundation of China University of Petroleum(Beijing)(No.KYJJ2012-05-02)
文摘The perfectly matched layer (PML) is a highly efficient absorbing boundary condition used for the numerical modeling of seismic wave equation. The article focuses on the application of this technique to finite-element time-domain numerical modeling of elastic wave equation. However, the finite-element time-domain scheme is based on the second- order wave equation in displacement formulation. Thus, the first-order PML in velocity-stress formulation cannot be directly applied to this scheme. In this article, we derive the finite- element matrix equations of second-order PML in displacement formulation, and accomplish the implementation of PML in finite-element time-domain modeling of elastic wave equation. The PML has an approximate zero reflection coefficients for bulk and surface waves in the finite-element modeling of P-SV and SH wave propagation in the 2D homogeneous elastic media. The numerical experiments using a two-layer model with irregular topography validate the efficiency of PML in the modeling of seismic wave propagation in geological models with complex structures and heterogeneous media.
文摘In this paper the decay of global solutions to some nonlinear dissipative wave equations are discussed, which based on the method of prior estimate technique and a differenece inequality.
文摘This paper deals with the initial-boundary value mixed problems for nonlinear wave equations. By introducing the 'blowing-up facts K(u,u_i)', We may discuss the blowing up behaviours of solutions in finite time to the mixed problems with respect to Neumann boundary and Dirichlet boundary for various nonlinear conditions and initial value conditions which usually meet.
文摘The existence and the nonexistence,the uniqueness and the energy decay estimate of solution for the fourth-order nonlinear wave equation utt+αΔ2 u-bΔut-βΔu+ut|ut|^r+g(u)=0 in Ω×(0,∞) are studied with the boundary condition u=(u)/(υ)=0 onΩ and the initial condition u(x,0)=u0(x),ut(x,0)=u1(x,0) in bounded domain ΩR^n ,n≥1.The energy decay rate of the global solution is estimated by the multiplier method.The blow-up result of the solution in finite time is established by the ideal of a potential well theory,and the existence of the solution is gotten by the Galekin approximation method.
基金supported by the National Science and Technology Major Project of China(Grant No. 2011ZX05004-003,2011ZX05014-006-006)the National Key Basic Research Program of China(Grant No. 2013CB228602)the Natural Science Foundation of China(Grant No. 40974066)
文摘The Finite Difference (FD) method is an important method for seismic numerical simulations. It helps us understand regular patterns in seismic wave propagation, analyze seismic attributes, and interpret seismic data. However, because of its discretization, the FD method is only stable under certain conditions. The Arbitrary Difference Precise Integration (ADPI) method is based on the FD method and adopts an integration scheme in the time domain and an arbitrary difference scheme in the space domain. Therefore, the ADPI method is a semi-analytical method. In this paper, we deduce the formula for the ADPI method based on the 3D elastic equation and improve its stability. In forward modeling cases, the ADPI method was implemented in 2D and 3D elastic wave equation forward modeling. Results show that the travel time of the reflected seismic wave is accurate. Compared with the acoustic wave field, the elastic wave field contains more wave types, including PS- and PP- reflected waves, transmitted waves, and diffracted waves, which is important to interpretation of seismic data. The method can be easily applied to elastic wave equation numerical simulations for eoloical models.