With respect to oceanic fluid dynamics,certain models have appeared,e.g.,an extended time-dependent(3+1)-dimensional shallow water wave equation in an ocean or a river,which we investigate in this paper.Using symbolic...With respect to oceanic fluid dynamics,certain models have appeared,e.g.,an extended time-dependent(3+1)-dimensional shallow water wave equation in an ocean or a river,which we investigate in this paper.Using symbolic computation,we find out,on one hand,a set of bilinear auto-Backlund transformations,which could connect certain solutions of that equation with other solutions of that equation itself,and on the other hand,a set of similarity reductions,which could go from that equation to a known ordinary differential equation.The results in this paper depend on all the oceanic variable coefficients in that equation.展开更多
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.展开更多
We investigate the blow-up effect of solutions for a non-homogeneous wave equation u_(tt)−∆u−∆u_(t)=I_(0+)^(α)(|u|^(p))+ω(x),where p>1,0≤α<1 andω(x)with∫_(R)^(N)ω(x)dx>0.By a way of combining the argum...We investigate the blow-up effect of solutions for a non-homogeneous wave equation u_(tt)−∆u−∆u_(t)=I_(0+)^(α)(|u|^(p))+ω(x),where p>1,0≤α<1 andω(x)with∫_(R)^(N)ω(x)dx>0.By a way of combining the argument by contradiction with the test function techniques,we prove that not only any non-trivial solution blows up in finite time under 0<α<1,N≥1 and p>1,but also any non-trivial solution blows up in finite time underα=0,2≤N≤4 and p being the Strauss exponent.展开更多
Internal multiple interference,affecting both seismic data processing and interpretation,has been observed for long time.Although great progress has been achieved in developing a variety of internal-multiple-eliminati...Internal multiple interference,affecting both seismic data processing and interpretation,has been observed for long time.Although great progress has been achieved in developing a variety of internal-multiple-elimination(IME)methods,how to increase accuracy and reduce cost of IME still poses a significant challenge.A new method is proposed to effectively and efficiently eliminate internal multi-ples,along with its application in internal-multiple-eliminated-migration(IMEM),addressing this issue.This method stems from two-way wave equation depth-extrapolation scheme and associated up/down wavefield separation,which can accomplish depth-extrapolation of both up-going and down-going wavefields simultaneously,and complete internal-multiple-elimination processing,adaptively and effi-ciently.The proposed method has several features:(1)input data is same as that for conventional migration:source signature(used for migration only),macro velocity model,and receiver data,without additional requirements for source/receiver sampling;(2)method is efficient,without need of iterative calculations(which are typically needed for most of IME algorithms);and(3)method is cost effective:IME is completed in the same depth-extrapolation scheme of IMEM,without need of a separate pro-cessing and additional cost.Several synthesized data models are used to test the proposed method:one-dimensional model,horizontal layered model,multi-layer model with one curved layer,and SEG/EAGE Salt model.Additionally,we perform a sensitivity analysis of velocity using smoothed models.This analysis reveals that although the accuracy of velocity measurements impacts our proposed method,it significantly reduces internal multiple false imaging compared to traditional RTM techniques.When applied to actual seismic data from a carbonate reservoir zone,our method demonstrates superior clarity in imaging results,even in the presence of high-velocity carbonate formations,outperforming conven-tional migration methods in deep strata.展开更多
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.展开更多
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.展开更多
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 consider two Cauchy systems of coupled two wave equations in the whole line R under one or two frictional dampings,where the coupling terms are either of order one with respect to the time variable or...In this paper,we consider two Cauchy systems of coupled two wave equations in the whole line R under one or two frictional dampings,where the coupling terms are either of order one with respect to the time variable or of order two with respect to the space variable.We prove some L^(2)(R)-norm decay estimates of solutions and their higher-order derivatives with respect to the space variable,where the decay rates depend on the number of the present frictional dampings,the regularity of the initial data,and some connections between the speeds of wave propagation of the two wave equations.Both our systems are considered under weaker conditions on the coefficients than the ones considered in the literature and they include the case where only one frictional damping is present,so they generate new difficulties and represent new situations that have not been studied earlier.展开更多
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.展开更多
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.展开更多
In this paper,the generalized Boussinesq wave equation u tt-uxx+a(um) xx+buxxxx=0 is investigated by using the bifurcation theory and the method of phase portraits analysis.Under the different parameter conditions,the...In this paper,the generalized Boussinesq wave equation u tt-uxx+a(um) xx+buxxxx=0 is investigated by using the bifurcation theory and the method of phase portraits analysis.Under the different parameter conditions,the exact explicit parametric representations for solitary wave solutions and periodic wave solutions are obtained.展开更多
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.展开更多
In seismic exploration, it is common practice to separate the P-wavefield from the S-wavefield by the elastic wavefield decomposition technique, for imaging purposes. However, it is sometimes difficult to achieve this...In seismic exploration, it is common practice to separate the P-wavefield from the S-wavefield by the elastic wavefield decomposition technique, for imaging purposes. However, it is sometimes difficult to achieve this, especially when the velocity field is complex. A useful approach in multi-component analysis and modeling is to directly solve the elastic wave equations for the pure P- or S-wavefields, referred as the separate elastic wave equa- tions. In this study, we compare two kinds of such wave equations： the first-order （velocity-stress） and the second- order （displacement-stress） separate elastic wave equa- tions, with the first-order （velocity-stress） and the second- order （displacement-stress） full （or mixed） elastic wave equations using a high-order staggered grid finite-differ- ence method. Comparisons are given of wavefield snap- shots, common-source gather seismic sections, and individual synthetic seismogram. The simulation tests show that equivalent results can be obtained, regardless of whether the first-order or second-order separate elastic wave equations are used for obtaining the pure P- or S-wavefield. The stacked pure P- and S-wavefields are equal to the mixed wave fields calculated using the corre- sponding first-order or second-order full elastic wave equations. These mixed equations are computationallyslightly less expensive than solving the separate equations. The attraction of the separate equations is that they achieve separated P- and S-wavefields which can be used to test the efficacy of wave decomposition procedures in multi-com- ponent processing. The second-order separate elastic wave equations are a good choice because they offer information on the pure P-wave or S-wave displacements.展开更多
The solutions to a linear wave equation can satisfy the principle of superposition, i.e., the linear superposition of two or more known solutions is still a solution of the linear wave equation. We show in this articl...The solutions to a linear wave equation can satisfy the principle of superposition, i.e., the linear superposition of two or more known solutions is still a solution of the linear wave equation. We show in this article that many nonlinear wave equations possess exact traveling wave solutions involving hyperbolic, triangle, and exponential functions, and the suitable linear combinations of these known solutions can also constitute linear superposition solutions to some nonlinear wave equations with special structural characteristics. The linear superposition solutions to the generalized KdV equation K（2,2,1）, the Oliver water wave equation, and the k（n, n） equation are given. The structure characteristic of the nonlinear wave equations having linear superposition solutions is analyzed, and the reason why the solutions with the forms of hyperbolic, triangle, and exponential functions can form the linear superposition solutions is also discussed.展开更多
Based on the pure quasi-P wave equation in transverse isotropic media with a vertical symmetry axis （VTI media）, a quasi-P wave equation is obtained in transverse isotropic media with a tilted symmetry axis （TTI me...Based on the pure quasi-P wave equation in transverse isotropic media with a vertical symmetry axis （VTI media）, a quasi-P wave equation is obtained in transverse isotropic media with a tilted symmetry axis （TTI media）. This is achieved using projection transformation, which rotates the direction vector in the coordinate system of observation toward the direction vector for the coordinate system in which the z-component is parallel to the symmetry axis of the TTI media. The equation has a simple form, is easily calculated, is not influenced by the pseudo-shear wave, and can be calculated reliably when δ is greater than ε. The finite difference method is used to solve the equation. In addition, a perfectly matched layer （PML） absorbing boundary condition is obtained for the equation. Theoretical analysis and numerical simulation results with forward modeling prove that the equation can accurately simulate a quasi-P wave in TTI medium.展开更多
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.展开更多
基金financially supported by the Scientific Research Foundation of North China University of Technology(Grant Nos.11005136024XN147-87 and 110051360024XN151-86).
文摘With respect to oceanic fluid dynamics,certain models have appeared,e.g.,an extended time-dependent(3+1)-dimensional shallow water wave equation in an ocean or a river,which we investigate in this paper.Using symbolic computation,we find out,on one hand,a set of bilinear auto-Backlund transformations,which could connect certain solutions of that equation with other solutions of that equation itself,and on the other hand,a set of similarity reductions,which could go from that equation to a known ordinary differential equation.The results in this paper depend on all the oceanic variable coefficients in that equation.
基金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 Natural Science Foundation of China(Grant No.62363005).
文摘We investigate the blow-up effect of solutions for a non-homogeneous wave equation u_(tt)−∆u−∆u_(t)=I_(0+)^(α)(|u|^(p))+ω(x),where p>1,0≤α<1 andω(x)with∫_(R)^(N)ω(x)dx>0.By a way of combining the argument by contradiction with the test function techniques,we prove that not only any non-trivial solution blows up in finite time under 0<α<1,N≥1 and p>1,but also any non-trivial solution blows up in finite time underα=0,2≤N≤4 and p being the Strauss exponent.
基金supported by the National Natural Science Foundation of China(Grant No.42004103)Sichuan Science and Technology Program(2023NSFSC0257)the CNPC Innovation Found(2022DQ02-0306).
文摘Internal multiple interference,affecting both seismic data processing and interpretation,has been observed for long time.Although great progress has been achieved in developing a variety of internal-multiple-elimination(IME)methods,how to increase accuracy and reduce cost of IME still poses a significant challenge.A new method is proposed to effectively and efficiently eliminate internal multi-ples,along with its application in internal-multiple-eliminated-migration(IMEM),addressing this issue.This method stems from two-way wave equation depth-extrapolation scheme and associated up/down wavefield separation,which can accomplish depth-extrapolation of both up-going and down-going wavefields simultaneously,and complete internal-multiple-elimination processing,adaptively and effi-ciently.The proposed method has several features:(1)input data is same as that for conventional migration:source signature(used for migration only),macro velocity model,and receiver data,without additional requirements for source/receiver sampling;(2)method is efficient,without need of iterative calculations(which are typically needed for most of IME algorithms);and(3)method is cost effective:IME is completed in the same depth-extrapolation scheme of IMEM,without need of a separate pro-cessing and additional cost.Several synthesized data models are used to test the proposed method:one-dimensional model,horizontal layered model,multi-layer model with one curved layer,and SEG/EAGE Salt model.Additionally,we perform a sensitivity analysis of velocity using smoothed models.This analysis reveals that although the accuracy of velocity measurements impacts our proposed method,it significantly reduces internal multiple false imaging compared to traditional RTM techniques.When applied to actual seismic data from a carbonate reservoir zone,our method demonstrates superior clarity in imaging results,even in the presence of high-velocity carbonate formations,outperforming conven-tional migration methods in deep strata.
基金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.
文摘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.
基金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 consider two Cauchy systems of coupled two wave equations in the whole line R under one or two frictional dampings,where the coupling terms are either of order one with respect to the time variable or of order two with respect to the space variable.We prove some L^(2)(R)-norm decay estimates of solutions and their higher-order derivatives with respect to the space variable,where the decay rates depend on the number of the present frictional dampings,the regularity of the initial data,and some connections between the speeds of wave propagation of the two wave equations.Both our systems are considered under weaker conditions on the coefficients than the ones considered in the literature and they include the case where only one frictional damping is present,so they generate new difficulties and represent new situations that have not been studied earlier.
基金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.
基金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.
基金Supported by the National Natural Science Foundation of China under Grant No. 10974160the Scientific Research Foundation of the Education Department of Sichuan Province of China under Grant No. 10ZA004
文摘In this paper,the generalized Boussinesq wave equation u tt-uxx+a(um) xx+buxxxx=0 is investigated by using the bifurcation theory and the method of phase portraits analysis.Under the different parameter conditions,the exact explicit parametric representations for solitary wave solutions and periodic wave solutions are obtained.
文摘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.
基金partially supported by China National Major Science and Technology Project (Subproject No:2011ZX05024-001-03)
文摘In seismic exploration, it is common practice to separate the P-wavefield from the S-wavefield by the elastic wavefield decomposition technique, for imaging purposes. However, it is sometimes difficult to achieve this, especially when the velocity field is complex. A useful approach in multi-component analysis and modeling is to directly solve the elastic wave equations for the pure P- or S-wavefields, referred as the separate elastic wave equa- tions. In this study, we compare two kinds of such wave equations： the first-order （velocity-stress） and the second- order （displacement-stress） separate elastic wave equa- tions, with the first-order （velocity-stress） and the second- order （displacement-stress） full （or mixed） elastic wave equations using a high-order staggered grid finite-differ- ence method. Comparisons are given of wavefield snap- shots, common-source gather seismic sections, and individual synthetic seismogram. The simulation tests show that equivalent results can be obtained, regardless of whether the first-order or second-order separate elastic wave equations are used for obtaining the pure P- or S-wavefield. The stacked pure P- and S-wavefields are equal to the mixed wave fields calculated using the corre- sponding first-order or second-order full elastic wave equations. These mixed equations are computationallyslightly less expensive than solving the separate equations. The attraction of the separate equations is that they achieve separated P- and S-wavefields which can be used to test the efficacy of wave decomposition procedures in multi-com- ponent processing. The second-order separate elastic wave equations are a good choice because they offer information on the pure P-wave or S-wave displacements.
文摘The solutions to a linear wave equation can satisfy the principle of superposition, i.e., the linear superposition of two or more known solutions is still a solution of the linear wave equation. We show in this article that many nonlinear wave equations possess exact traveling wave solutions involving hyperbolic, triangle, and exponential functions, and the suitable linear combinations of these known solutions can also constitute linear superposition solutions to some nonlinear wave equations with special structural characteristics. The linear superposition solutions to the generalized KdV equation K（2,2,1）, the Oliver water wave equation, and the k（n, n） equation are given. The structure characteristic of the nonlinear wave equations having linear superposition solutions is analyzed, and the reason why the solutions with the forms of hyperbolic, triangle, and exponential functions can form the linear superposition solutions is also discussed.
基金supported by the National Natural Science Foundation of China(No.41674118)the national science and technology major project(No.2016ZX05027-002)
文摘Based on the pure quasi-P wave equation in transverse isotropic media with a vertical symmetry axis （VTI media）, a quasi-P wave equation is obtained in transverse isotropic media with a tilted symmetry axis （TTI media）. This is achieved using projection transformation, which rotates the direction vector in the coordinate system of observation toward the direction vector for the coordinate system in which the z-component is parallel to the symmetry axis of the TTI media. The equation has a simple form, is easily calculated, is not influenced by the pseudo-shear wave, and can be calculated reliably when δ is greater than ε. The finite difference method is used to solve the equation. In addition, a perfectly matched layer （PML） absorbing boundary condition is obtained for the equation. Theoretical analysis and numerical simulation results with forward modeling prove that the equation can accurately simulate a quasi-P wave in TTI medium.
基金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.
