Carbonate reservoirs are vital energy storage spaces,including for oil,shale gas,geothermal,and hydrogen energy.Accurate prediction of reservoir characteristics such as permeability and saturated fluid types through n...Carbonate reservoirs are vital energy storage spaces,including for oil,shale gas,geothermal,and hydrogen energy.Accurate prediction of reservoir characteristics such as permeability and saturated fluid types through noninvasive approaches is crucial for optimal storage capability.In this paper,we combine a linear Boolean model and a discrete Fourier transform approach to generate pore‐and fracturepore‐type carbonate rocks.Elastic wave velocity information is necessary to predict permeability in different rock geometry models.Permeability is calculated using the lattice Boltzmann method,and the elastic wave velocity is calculated using a finite element method based on a minimal energy approach.Saturated fluids that contain oil and gas were both considered.Our simulated results reveal that,for pore‐type carbonate,empirical formulas were proposed to estimate permeability through elastic data.However,in fracture‐pore carbonate rocks,the precision of the empirical formula is compromised due to the presence of significant conductive channels within the rock matrix.We also find that using S‐wave velocity and permeability relationships to distinguish oil and gas is better than using P‐wave velocity and permeability relationships under low‐porosity conditions.展开更多
The P-SV wave reflection coefficients in VTI and HTI media were obtained by approximation of the Jilek（2002a and b）equation in orthotropic anisotropic media.An approximate equation for P-SV wave elastic impedance ca...The P-SV wave reflection coefficients in VTI and HTI media were obtained by approximation of the Jilek（2002a and b）equation in orthotropic anisotropic media.An approximate equation for P-SV wave elastic impedance can be derived from the combination of the new coefficients with S-wave elastic impedance（Duffaut et al.,2000）.On this basis, the fluid identification factor in weakly anisotropic media was constructed and used to identify the Castagna and Smith（1994）lithologic combination and achieved good results. Finally,we specifically analyzed the anisotropic parameter impacts P-SV wave elastic impedance and fluid factor trends.展开更多
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.展开更多
Knowledge of the locations of seismic sources is critical for microseismic monitoring. Time-window-based elastic wave interferometric imaging and weighted- elastic-wave （WEW） interferometric imaging are proposed and...Knowledge of the locations of seismic sources is critical for microseismic monitoring. Time-window-based elastic wave interferometric imaging and weighted- elastic-wave （WEW） interferometric imaging are proposed and used to locate modeled microseismic sources. The proposed method improves the precision and eliminates artifacts in location profiles. Numerical experiments based on a horizontally layered isotropic medium have shown that the method offers the following advantages： It can deal with Iow-SNR microseismic data with velocity perturbations as well as relatively sparse receivers and still maintain relatively high precision despite the errors in the velocity model. Furthermore, it is more efficient than conventional traveltime inversion methods because interferometric imaging does not require traveltime picking. Numerical results using a 2D fault model have also suggested that the weighted-elastic-wave interferometric imaging can locate multiple sources with higher location precision than the time-reverse imaging 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.展开更多
In this study,the wave propagation properties of lattice metamaterials with Koch fractal structures are investigated in terms of band structures and directional wave propagation.The analytical models of lattice metama...In this study,the wave propagation properties of lattice metamaterials with Koch fractal structures are investigated in terms of band structures and directional wave propagation.The analytical models of lattice metamaterials are established using the finite element method,and the dispersion relation is solved using the Bloch’s theorem.The band structures of the lattice metamaterials with different numbers of iterations are studied,and the group velocities at a selected frequency are calculated to analyze the directional wave propagation characteristics.Furthermore,dynamic responses of the finite structures are calculated using commercial finite element software to verify the band gaps and directional wave propagation behaviors in the lattice metamaterials.The results show that multiple and low band gaps are present in the lattice materials with various geometric parameters of the Koch fractal,and the position of the lowest band gap decreases as the number of iterations increases.The results indicate the potential applications of lattice metamaterials with Koch fractals for vibration isolation and multi-functional design.展开更多
This article provides the application of the high-order, staggered-grid, finite-difference scheme to model elastic wave propagation in 3-D isotropic media. Here, we use second-order, tempo- ral- and high-order spatial...This article provides the application of the high-order, staggered-grid, finite-difference scheme to model elastic wave propagation in 3-D isotropic media. Here, we use second-order, tempo- ral- and high-order spatial finite-difference formulations with a staggered grid for discretization of the 3-D elastic wave equations of motion. The set of absorbing boundary conditions based on paraxial approximations of 3-D elastic wave equations are applied to the numerical boundaries. The trial re- sults for the salt model show that the numerical dispersion is decreased to a minimum extent, the accuracy high and diffracted waves abundant. It also shows that this method can be used for modeling wave propagation in complex media with the lateral variation of velocity.展开更多
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.展开更多
It is an important subject to probe the structure in the medium by various kinds of detection methods in the geotechnical engineering. Based on the propagation theory of elastic wave in half-space layered medium, the ...It is an important subject to probe the structure in the medium by various kinds of detection methods in the geotechnical engineering. Based on the propagation theory of elastic wave in half-space layered medium, the propagation characteristics of elastic wave in layered medium with different elastic parameters are discussed using dynamic analysis of finite element method. It is known that the S-wave velocity, density and thickness of layer are related to the properties of the elastic wave including waveform characteristics, spectral characteristics and time-frequency characteristics. We pay special attention to the structure with low velocity interlayer. The impact imaging method is applied to the grouting construction of the immersed tube tunnel. Data acquisition and analytical method are introduced in detail. The grouting effects can be qualitatively evaluated by comparing the characteristics of elastic wave before grouting with those after grouting. Finally, a quantitative evaluation is obtained according to the relationship between energy response of elastic wave and impedance ratio.展开更多
Due to their potential properties unlike traditional materials and structures,elastic wave metamaterials have received significant interests in recent years.With the coupling between the acoustic and vibration,their m...Due to their potential properties unlike traditional materials and structures,elastic wave metamaterials have received significant interests in recent years.With the coupling between the acoustic and vibration,their mechanical characteristics can be tuned by the active feedback control system at low frequency ranges in which the traditional passive control is limited.This work illustrates that the superior performances of the effective mass density and sound pressure level(SPL)of an elastic wave metamaterial can be significantly changed by the active control,in which the periodic array of local resonators and orthogonal stiffeners are included.Significantly,based on the locally resonant mechanism,the negative density occurs over a frequency range.Due to the effects of lattice constant,structural damping and other parameters,the SPL with the function of fluid-solid coupling are illustrated and discussed.展开更多
The propagation of elastic waves is studied in a porous solid saturated with two immiscible viscous fluids. The propagation of three longitudinal waves is represented through three scalar potential functions. The lone...The propagation of elastic waves is studied in a porous solid saturated with two immiscible viscous fluids. The propagation of three longitudinal waves is represented through three scalar potential functions. The lone transverse wave is presented by a vector potential function. The displacements of particles in different phases of the aggregate are defined in terms of these potential functions. It is shown that there exist three longitudinal waves and one transverse wave. The phenomena of reflection and refraction due to longitudinal and transverse waves at a plane interface between an elastic solid half-space and a porous solid half-space saturated with two immiscible viscous fluids are investigated. For the presence of viscosity in pore-fluids, the waves refracted to the porous medium attenuate in the direction normal to the interface. The ratios of the amplitudes of the reflected and refracted waves to that of the incident wave are calculated as a non- singular system of linear algebraic equations. These amplitude ratios are used to further calculate the shares of different scattered waves in the energy of the incident wave. The modulus of the amplitude and the energy ratios with the angle of incidence are computed for a particular numerical model. The conservation of the energy across the interface is verified. The effects of variations in non-wet saturation of pores and frequencies on the energy partition are depicted graphically and discussed,展开更多
A Harten-Lax-van Leer-contact （HLLC） approximate Riemann solver is built with elastic waves （HLLCE） for one-dimensional elastic-plastic flows with a hypo- elastic constitutive model and the von Mises＇ yielding cr...A Harten-Lax-van Leer-contact （HLLC） approximate Riemann solver is built with elastic waves （HLLCE） for one-dimensional elastic-plastic flows with a hypo- elastic constitutive model and the von Mises＇ yielding criterion. Based on the HLLCE, a third-order cell-centered Lagrangian scheme is built for one-dimensional elastic-plastic problems. A number of numerical experiments are carried out. The numerical results show that the proposed third-order scheme achieves the desired order of accuracy. The third-order scheme is used to the numerical solution of the problems with elastic shock waves and elastic rarefaction waves. The numerical results are compared with a reference solution and the results obtained by other authors. The comparison shows that the pre- sented high-order scheme is convergent, stable, and essentially non-oscillatory. Moreover, the HLLCE is more efficient than the two-rarefaction Riemann solver with elastic waves （TRRSE）展开更多
Mesh-free finite difference(FD)methods can improve the geometric flexibility of modeling without the need for lattice mapping or complex meshing process.Radial-basisfunction-generated FD is among the most commonly use...Mesh-free finite difference(FD)methods can improve the geometric flexibility of modeling without the need for lattice mapping or complex meshing process.Radial-basisfunction-generated FD is among the most commonly used mesh-free FD methods and can accurately simulate seismic wave propagation in the non-rectangular computational domain.In this paper,we propose a perfectly matched layer(PML)boundary condition for a meshfree FD solution of the elastic wave equation,which can be applied to the boundaries of the non-rectangular velocity model.The performance of the PML is,however,severely reduced for near-grazing incident waves and low-frequency waves.We thus also propose the complexfrequency-shifted PML(CFS-PML)boundary condition for a mesh-free FD solution of the elastic wave equation.For two PML boundary conditions,we derive unsplit time-domain expressions by constructing auxiliary differential equations,both of which require less memory and are easy for programming.Numerical experiments demonstrate that these two PML boundary conditions effectively eliminate artificial boundary reflections in mesh-free FD simulations.When compared with the PML boundary condition,the CFS-PML boundary condition results in better absorption for near-grazing incident waves and evanescent waves.展开更多
Using the active feedback control system on the elastic wave metamaterial,this research concentrates on the sound transmission with the dynamic effective model.The metamaterial is subjected to an incident pressure and...Using the active feedback control system on the elastic wave metamaterial,this research concentrates on the sound transmission with the dynamic effective model.The metamaterial is subjected to an incident pressure and immersed in the external mean flow.The elastic wave metamaterial consists of double plates and the upper and lower four-link mechanisms are attached inside.The vertical resonator is attached by the active feedback control system and connected with two four-link mechanisms.Based on the dynamic equivalent method,the metamaterial is equivalent as a single-layer plate by the dynamic effective parameter.With the coupling between the fluid and structure,the expression of the sound transmission loss(STL)is derived.This research shows the influence of effective mass density on sound transmission properties,and the STL in both modes can be tuned by the acceleration and displacement feedback constants.In addition,the dynamic response and the STL are also changed obviously by different values of structural damping,incident angle(i.e.,the elevation and azimuth angles)and Mach number of the external fluid with the mean flow property.The results for sound transmission by two methods are compared,i.e.,the virtual work principle for double plates and the dynamic equivalent method corresponding to a single one.This paper is expected to be helpful for understanding the sound transmission properties of both pure single-and double-plate models.展开更多
In this work,the three-dimensional(3 D)propagation behaviors in the nonlinear phononic crystal and elastic wave metamaterial with initial stresses are investigated.The analytical solutions of the fundamental wave and ...In this work,the three-dimensional(3 D)propagation behaviors in the nonlinear phononic crystal and elastic wave metamaterial with initial stresses are investigated.The analytical solutions of the fundamental wave and second harmonic with the quasilongitudinal(qP)and quasi-shear(qS_(1) and qS_(2))modes are derived.Based on the transfer and stiffness matrices,band gaps with initial stresses are obtained by the Bloch theorem.The transmission coefficients are calculated to support the band gap property,and the tunability of the nonreciprocal transmission by the initial stress is discussed.This work is expected to provide a way to tune the nonreciprocal transmission with vector characteristics.展开更多
A convolution perfectly matched layer(CPML)can efficiently absorb boundary reflection in numerical simulation.However,the CPML is suitable for the first-order elastic wave equation and is difficult to apply directly t...A convolution perfectly matched layer(CPML)can efficiently absorb boundary reflection in numerical simulation.However,the CPML is suitable for the first-order elastic wave equation and is difficult to apply directly to the second-order elastic wave equation.In view of this,based on the first-order CPML absorbing boundary condition,we propose a new CPML(NCPML)boundary which can be directly applied to the second-order wave equation.We first systematically extend the first-order CPML technique into second-order wave equations,neglecting the space-varying characteristics of the partial damping coefficient in the complex-frequency domain,avoiding the generation of convolution in the time domain.We then transform the technique back to the time domain through the inverse Fourier transform.Numerical simulation indicates that the space-varying characteristics of the attenuation factor have little influence on the absorption effect and increase the memory at the same time.A number of numerical examples show that the NCPML proposed in this study is effective in simulating elastic wave propagation,and this algorithm is more efficient and requires less memory allocation than the conventional PML absorbing boundary.展开更多
Edge reflections are inevitable in numerical modeling of seismic wavefields, and they are usually attenuated by absorbing boundary conditions. However, the commonly used perfectly matched layer （PML） boundary condit...Edge reflections are inevitable in numerical modeling of seismic wavefields, and they are usually attenuated by absorbing boundary conditions. However, the commonly used perfectly matched layer （PML） boundary condition requires special treatment for the absorbing zone, and in three-dimensional （3D） modeling, it has to split each variable into three corresponding variables, which increases the computing time and memory storage. In contrast, the hybrid absorbing boundary condition （HABC） has the advantages such as ease of implementation, less computation time, and near-perfect absorption; it is thus able to enhance the computational efficiency of 3D elastic wave modeling. In this study, a HABC is developed from two-dimensional （2D） modeling into 3D modeling based on the I st Higdon one way wave equations, and a HABC is proposed that is suitable for a 3D elastic wave numerical simulation. Numerical simulation results for a homogenous model and a complex model indicate that the proposed HABC method is more effective and has better absorption than the traditional PML method.展开更多
To understand the evolution of stress-induced elastic wave anisotropy,three triaxial experiments were performed on sandstone specimens with bedding orientations parallel,perpendicular,and oblique to the maximum princi...To understand the evolution of stress-induced elastic wave anisotropy,three triaxial experiments were performed on sandstone specimens with bedding orientations parallel,perpendicular,and oblique to the maximum principal stress.P-wave velocities along 64 different directions on each specimen were monitored frequently to understand the anisotropy change at various stress levels by fitting Thomsen’s anisotropy equation.The results show that the elastic wave anisotropy is very sensitive to mechanical loading.Under hydrostatic loading,the magnitude of anisotropy is reduced in all three specimens.However,under deviatoric stress loading,the evolution of anisotropic characteristics(magnitude and orientation of the symmetry axis)is bedding orientation dependent.Anisotropy reversal occurs in specimens with bedding normal/oblique to the maximum principal stress.P-wave anisotropyε0 is linearly related to volumetric strain Sv and dilatancy,indicating that stress-induced redistribution of microcracks has a significant effect on P-wave velocity anisotropy.The closure of initial cracks and pores aligned in the bedding direction contributes to the decrease of the anisotropy.However,opening of new cracks,aligned in the maximum principal direction,accounts for the increase of the anisotropy.The experimental results provide some insights into the microstructural behavior under loading and provide an experimental basis for seismic data interpretation and parameter selection in engineering applications.展开更多
Small moving vehicles represent an important category of marine engineering tools and devices(equipment)typically used for ocean resource detection and maintenance of marine rights and interests.The lack of efficient ...Small moving vehicles represent an important category of marine engineering tools and devices(equipment)typically used for ocean resource detection and maintenance of marine rights and interests.The lack of efficient power supply modes is one of the technical bottlenecks restricting the effective utilisation of this type of equipment.In this work,the performance characteristics of a new type of elastic-blade/wave-energy converter(EBWEC)and its core energy conversion component(named wave energy absorber)are comprehensively studied.In particular,computational fluid dynamics(CFD)simulations and experiments have been used to analyze the hydrodynamics and performance characteristics of the EBWEC.The pressure cloud diagrams relating to the surface of the elastic blade were obtained through two-way fluid-solid coupling simulations.The influence of blade thickness and relative speed on the performance characteristics of EBWEC was analyzed accordingly.A prototype of the EBWEC and its bucket test platform were also developed.The power characteristics of the EBWEC were analyzed and studied by using the blade thickness and motion cycle as control variables.The present research shows that the EBWEC can effectively overcome the performance disadvantages related to the transmission shaft torque load and power curve fluctuations of rigid blade wave energy converters(RBWEC).展开更多
It is essential to assess the evolution of soil fabric as it has an important role in the mechanical responses of soils during complex loading conditions.This contribution carries out the physical experiments using th...It is essential to assess the evolution of soil fabric as it has an important role in the mechanical responses of soils during complex loading conditions.This contribution carries out the physical experiments using three granular materials in the laboratory.The variations of compression and shear wave velocities(Vp and Vs)are investigated during load-unload cycles under dry and drained conditions.Supplementary discrete element method(DEM)simulations are performed to understand the evolution of soil fabric during the equivalent load-unload cycles using spherical particles.Vp and Vs are not always reversible even though the stress state regains its isotropic condition after unload,indicating that Vp and Vs are governed by not only the stress state but also the fabric change.The variations of Vp/Vs are density-and stress-dependent;a higher level of stress ratio(s01/s03)threshold is observed for denser packings to trigger a significant change in wave velocity ratio(Vp/Vs)for experimental results using spherical glass beads and simulation data using spherical particles.Considering the particle shape,a higher s01/s03 threshold is found for more angular particles than rounded particles.The DEM result reveals that Vp/Vs of spherical particles can be correlated linearly with the evolution of fabric ratio(Fver/Fhor)during loadunload in a pre-peak range under dry and drained conditions.展开更多
基金Chengdu University of Technology Youth Teaching Backbone Project,Grant/Award Number:10912-JXGG2023-09458Sichuan Province Overseas Returnees Science and Technology Excellence Project,Grant/Award Number:10900-23BZ28-02+1 种基金Physics-informed machine learning for coupled modelling of carbon storage,Grant/Award Number:SKLGME021002Chengdu University of Technology Research Start-up Fund,Grant/Award Number:10912-KYQD2022-09458。
文摘Carbonate reservoirs are vital energy storage spaces,including for oil,shale gas,geothermal,and hydrogen energy.Accurate prediction of reservoir characteristics such as permeability and saturated fluid types through noninvasive approaches is crucial for optimal storage capability.In this paper,we combine a linear Boolean model and a discrete Fourier transform approach to generate pore‐and fracturepore‐type carbonate rocks.Elastic wave velocity information is necessary to predict permeability in different rock geometry models.Permeability is calculated using the lattice Boltzmann method,and the elastic wave velocity is calculated using a finite element method based on a minimal energy approach.Saturated fluids that contain oil and gas were both considered.Our simulated results reveal that,for pore‐type carbonate,empirical formulas were proposed to estimate permeability through elastic data.However,in fracture‐pore carbonate rocks,the precision of the empirical formula is compromised due to the presence of significant conductive channels within the rock matrix.We also find that using S‐wave velocity and permeability relationships to distinguish oil and gas is better than using P‐wave velocity and permeability relationships under low‐porosity conditions.
基金sponsored by the National 973 Program(Grant No.2007CB209603)
文摘The P-SV wave reflection coefficients in VTI and HTI media were obtained by approximation of the Jilek（2002a and b）equation in orthotropic anisotropic media.An approximate equation for P-SV wave elastic impedance can be derived from the combination of the new coefficients with S-wave elastic impedance（Duffaut et al.,2000）.On this basis, the fluid identification factor in weakly anisotropic media was constructed and used to identify the Castagna and Smith（1994）lithologic combination and achieved good results. Finally,we specifically analyzed the anisotropic parameter impacts P-SV wave elastic impedance and fluid factor trends.
基金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.
基金supported by the R&D of Key Instruments and Technologies for Deep Resources Prospecting(No.ZDYZ2012-1)National Natural Science Foundation of China(No.11374322)
文摘Knowledge of the locations of seismic sources is critical for microseismic monitoring. Time-window-based elastic wave interferometric imaging and weighted- elastic-wave （WEW） interferometric imaging are proposed and used to locate modeled microseismic sources. The proposed method improves the precision and eliminates artifacts in location profiles. Numerical experiments based on a horizontally layered isotropic medium have shown that the method offers the following advantages： It can deal with Iow-SNR microseismic data with velocity perturbations as well as relatively sparse receivers and still maintain relatively high precision despite the errors in the velocity model. Furthermore, it is more efficient than conventional traveltime inversion methods because interferometric imaging does not require traveltime picking. Numerical results using a 2D fault model have also suggested that the weighted-elastic-wave interferometric imaging can locate multiple sources with higher location precision than the time-reverse imaging 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.
基金Funding for this work has been provided by the National Natural Science Foundation of China(Nos.11872313 and 11502202)National Key R&D Program of China(2017YFB1102801)Fundamental Research Funds for the Central Universities and Seed Foundation of Innovation and Creation for Graduate Students in Northwestern Polytechnical University(CX2020107).
文摘In this study,the wave propagation properties of lattice metamaterials with Koch fractal structures are investigated in terms of band structures and directional wave propagation.The analytical models of lattice metamaterials are established using the finite element method,and the dispersion relation is solved using the Bloch’s theorem.The band structures of the lattice metamaterials with different numbers of iterations are studied,and the group velocities at a selected frequency are calculated to analyze the directional wave propagation characteristics.Furthermore,dynamic responses of the finite structures are calculated using commercial finite element software to verify the band gaps and directional wave propagation behaviors in the lattice metamaterials.The results show that multiple and low band gaps are present in the lattice materials with various geometric parameters of the Koch fractal,and the position of the lowest band gap decreases as the number of iterations increases.The results indicate the potential applications of lattice metamaterials with Koch fractals for vibration isolation and multi-functional design.
文摘This article provides the application of the high-order, staggered-grid, finite-difference scheme to model elastic wave propagation in 3-D isotropic media. Here, we use second-order, tempo- ral- and high-order spatial finite-difference formulations with a staggered grid for discretization of the 3-D elastic wave equations of motion. The set of absorbing boundary conditions based on paraxial approximations of 3-D elastic wave equations are applied to the numerical boundaries. The trial re- sults for the salt model show that the numerical dispersion is decreased to a minimum extent, the accuracy high and diffracted waves abundant. It also shows that this method can be used for modeling wave propagation in complex media with the lateral variation of velocity.
文摘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 National Basic Research Program (973) of China(No.2011CB013505)
文摘It is an important subject to probe the structure in the medium by various kinds of detection methods in the geotechnical engineering. Based on the propagation theory of elastic wave in half-space layered medium, the propagation characteristics of elastic wave in layered medium with different elastic parameters are discussed using dynamic analysis of finite element method. It is known that the S-wave velocity, density and thickness of layer are related to the properties of the elastic wave including waveform characteristics, spectral characteristics and time-frequency characteristics. We pay special attention to the structure with low velocity interlayer. The impact imaging method is applied to the grouting construction of the immersed tube tunnel. Data acquisition and analytical method are introduced in detail. The grouting effects can be qualitatively evaluated by comparing the characteristics of elastic wave before grouting with those after grouting. Finally, a quantitative evaluation is obtained according to the relationship between energy response of elastic wave and impedance ratio.
基金the supports by the National Natural Science Foundation of China(Grants 11922209,11991031 and 12021002)for this research work.
文摘Due to their potential properties unlike traditional materials and structures,elastic wave metamaterials have received significant interests in recent years.With the coupling between the acoustic and vibration,their mechanical characteristics can be tuned by the active feedback control system at low frequency ranges in which the traditional passive control is limited.This work illustrates that the superior performances of the effective mass density and sound pressure level(SPL)of an elastic wave metamaterial can be significantly changed by the active control,in which the periodic array of local resonators and orthogonal stiffeners are included.Significantly,based on the locally resonant mechanism,the negative density occurs over a frequency range.Due to the effects of lattice constant,structural damping and other parameters,the SPL with the function of fluid-solid coupling are illustrated and discussed.
基金Project supported by the Council of Scientific and Industrial Research (CSIR) of New Delhi(Nos. 09/105(0169)/2008-EMR-I and 09/105(0185)/2009-EMR-I)
文摘The propagation of elastic waves is studied in a porous solid saturated with two immiscible viscous fluids. The propagation of three longitudinal waves is represented through three scalar potential functions. The lone transverse wave is presented by a vector potential function. The displacements of particles in different phases of the aggregate are defined in terms of these potential functions. It is shown that there exist three longitudinal waves and one transverse wave. The phenomena of reflection and refraction due to longitudinal and transverse waves at a plane interface between an elastic solid half-space and a porous solid half-space saturated with two immiscible viscous fluids are investigated. For the presence of viscosity in pore-fluids, the waves refracted to the porous medium attenuate in the direction normal to the interface. The ratios of the amplitudes of the reflected and refracted waves to that of the incident wave are calculated as a non- singular system of linear algebraic equations. These amplitude ratios are used to further calculate the shares of different scattered waves in the energy of the incident wave. The modulus of the amplitude and the energy ratios with the angle of incidence are computed for a particular numerical model. The conservation of the energy across the interface is verified. The effects of variations in non-wet saturation of pores and frequencies on the energy partition are depicted graphically and discussed,
基金Project supported by the National Natural Science Foundation of China(Nos.11172050 and11672047)the Science and Technology Foundation of China Academy of Engineering Physics(No.2013A0202011)
文摘A Harten-Lax-van Leer-contact （HLLC） approximate Riemann solver is built with elastic waves （HLLCE） for one-dimensional elastic-plastic flows with a hypo- elastic constitutive model and the von Mises＇ yielding criterion. Based on the HLLCE, a third-order cell-centered Lagrangian scheme is built for one-dimensional elastic-plastic problems. A number of numerical experiments are carried out. The numerical results show that the proposed third-order scheme achieves the desired order of accuracy. The third-order scheme is used to the numerical solution of the problems with elastic shock waves and elastic rarefaction waves. The numerical results are compared with a reference solution and the results obtained by other authors. The comparison shows that the pre- sented high-order scheme is convergent, stable, and essentially non-oscillatory. Moreover, the HLLCE is more efficient than the two-rarefaction Riemann solver with elastic waves （TRRSE）
基金supported by the National Science and Technology Major Project(2016ZX05006-002)the National Natural Science Foundation of China(Nos.41874153,41504097)
文摘Mesh-free finite difference(FD)methods can improve the geometric flexibility of modeling without the need for lattice mapping or complex meshing process.Radial-basisfunction-generated FD is among the most commonly used mesh-free FD methods and can accurately simulate seismic wave propagation in the non-rectangular computational domain.In this paper,we propose a perfectly matched layer(PML)boundary condition for a meshfree FD solution of the elastic wave equation,which can be applied to the boundaries of the non-rectangular velocity model.The performance of the PML is,however,severely reduced for near-grazing incident waves and low-frequency waves.We thus also propose the complexfrequency-shifted PML(CFS-PML)boundary condition for a mesh-free FD solution of the elastic wave equation.For two PML boundary conditions,we derive unsplit time-domain expressions by constructing auxiliary differential equations,both of which require less memory and are easy for programming.Numerical experiments demonstrate that these two PML boundary conditions effectively eliminate artificial boundary reflections in mesh-free FD simulations.When compared with the PML boundary condition,the CFS-PML boundary condition results in better absorption for near-grazing incident waves and evanescent waves.
基金The authors wish to express gratitude for the support provided by the National Natural Science Foundation of China(Grant Nos.11922209,11991031 and 12021002).
文摘Using the active feedback control system on the elastic wave metamaterial,this research concentrates on the sound transmission with the dynamic effective model.The metamaterial is subjected to an incident pressure and immersed in the external mean flow.The elastic wave metamaterial consists of double plates and the upper and lower four-link mechanisms are attached inside.The vertical resonator is attached by the active feedback control system and connected with two four-link mechanisms.Based on the dynamic equivalent method,the metamaterial is equivalent as a single-layer plate by the dynamic effective parameter.With the coupling between the fluid and structure,the expression of the sound transmission loss(STL)is derived.This research shows the influence of effective mass density on sound transmission properties,and the STL in both modes can be tuned by the acceleration and displacement feedback constants.In addition,the dynamic response and the STL are also changed obviously by different values of structural damping,incident angle(i.e.,the elevation and azimuth angles)and Mach number of the external fluid with the mean flow property.The results for sound transmission by two methods are compared,i.e.,the virtual work principle for double plates and the dynamic equivalent method corresponding to a single one.This paper is expected to be helpful for understanding the sound transmission properties of both pure single-and double-plate models.
基金Project supported by the National Natural Science Foundation of China(Nos.11922209,11991031 and 12021002)。
文摘In this work,the three-dimensional(3 D)propagation behaviors in the nonlinear phononic crystal and elastic wave metamaterial with initial stresses are investigated.The analytical solutions of the fundamental wave and second harmonic with the quasilongitudinal(qP)and quasi-shear(qS_(1) and qS_(2))modes are derived.Based on the transfer and stiffness matrices,band gaps with initial stresses are obtained by the Bloch theorem.The transmission coefficients are calculated to support the band gap property,and the tunability of the nonreciprocal transmission by the initial stress is discussed.This work is expected to provide a way to tune the nonreciprocal transmission with vector characteristics.
基金supported by the National Science and Technology Major Special Sub-project of China(No.2016ZX05024-001-008)the National Natural Science Foundation Joint Fund Prcject of China(No.U1562215).
文摘A convolution perfectly matched layer(CPML)can efficiently absorb boundary reflection in numerical simulation.However,the CPML is suitable for the first-order elastic wave equation and is difficult to apply directly to the second-order elastic wave equation.In view of this,based on the first-order CPML absorbing boundary condition,we propose a new CPML(NCPML)boundary which can be directly applied to the second-order wave equation.We first systematically extend the first-order CPML technique into second-order wave equations,neglecting the space-varying characteristics of the partial damping coefficient in the complex-frequency domain,avoiding the generation of convolution in the time domain.We then transform the technique back to the time domain through the inverse Fourier transform.Numerical simulation indicates that the space-varying characteristics of the attenuation factor have little influence on the absorption effect and increase the memory at the same time.A number of numerical examples show that the NCPML proposed in this study is effective in simulating elastic wave propagation,and this algorithm is more efficient and requires less memory allocation than the conventional PML absorbing boundary.
基金supported by the National Natural Science Foundation of China(No.41474110)
文摘Edge reflections are inevitable in numerical modeling of seismic wavefields, and they are usually attenuated by absorbing boundary conditions. However, the commonly used perfectly matched layer （PML） boundary condition requires special treatment for the absorbing zone, and in three-dimensional （3D） modeling, it has to split each variable into three corresponding variables, which increases the computing time and memory storage. In contrast, the hybrid absorbing boundary condition （HABC） has the advantages such as ease of implementation, less computation time, and near-perfect absorption; it is thus able to enhance the computational efficiency of 3D elastic wave modeling. In this study, a HABC is developed from two-dimensional （2D） modeling into 3D modeling based on the I st Higdon one way wave equations, and a HABC is proposed that is suitable for a 3D elastic wave numerical simulation. Numerical simulation results for a homogenous model and a complex model indicate that the proposed HABC method is more effective and has better absorption than the traditional PML method.
基金The research was partially supported by the National Natural Science Foundation of China(Grant Nos.41902297,41872210)the Natural Science Foundation of Hubei Province(Grant No.2018CFB292)Open Research Fund of State Key Laboratory of Geomechanics and Geotechnical Engineering,Institute of Rock and Soil Mechanics,Chinese Academy of Sciences(Grant No.Z017006).
文摘To understand the evolution of stress-induced elastic wave anisotropy,three triaxial experiments were performed on sandstone specimens with bedding orientations parallel,perpendicular,and oblique to the maximum principal stress.P-wave velocities along 64 different directions on each specimen were monitored frequently to understand the anisotropy change at various stress levels by fitting Thomsen’s anisotropy equation.The results show that the elastic wave anisotropy is very sensitive to mechanical loading.Under hydrostatic loading,the magnitude of anisotropy is reduced in all three specimens.However,under deviatoric stress loading,the evolution of anisotropic characteristics(magnitude and orientation of the symmetry axis)is bedding orientation dependent.Anisotropy reversal occurs in specimens with bedding normal/oblique to the maximum principal stress.P-wave anisotropyε0 is linearly related to volumetric strain Sv and dilatancy,indicating that stress-induced redistribution of microcracks has a significant effect on P-wave velocity anisotropy.The closure of initial cracks and pores aligned in the bedding direction contributes to the decrease of the anisotropy.However,opening of new cracks,aligned in the maximum principal direction,accounts for the increase of the anisotropy.The experimental results provide some insights into the microstructural behavior under loading and provide an experimental basis for seismic data interpretation and parameter selection in engineering applications.
基金financially supported by the National Natural Science Foundation of China(Grant Number 51475465)the Hunan Provincial Innovation Foundation for Postgraduate(Grant Number CX2015B014).
文摘Small moving vehicles represent an important category of marine engineering tools and devices(equipment)typically used for ocean resource detection and maintenance of marine rights and interests.The lack of efficient power supply modes is one of the technical bottlenecks restricting the effective utilisation of this type of equipment.In this work,the performance characteristics of a new type of elastic-blade/wave-energy converter(EBWEC)and its core energy conversion component(named wave energy absorber)are comprehensively studied.In particular,computational fluid dynamics(CFD)simulations and experiments have been used to analyze the hydrodynamics and performance characteristics of the EBWEC.The pressure cloud diagrams relating to the surface of the elastic blade were obtained through two-way fluid-solid coupling simulations.The influence of blade thickness and relative speed on the performance characteristics of EBWEC was analyzed accordingly.A prototype of the EBWEC and its bucket test platform were also developed.The power characteristics of the EBWEC were analyzed and studied by using the blade thickness and motion cycle as control variables.The present research shows that the EBWEC can effectively overcome the performance disadvantages related to the transmission shaft torque load and power curve fluctuations of rigid blade wave energy converters(RBWEC).
文摘It is essential to assess the evolution of soil fabric as it has an important role in the mechanical responses of soils during complex loading conditions.This contribution carries out the physical experiments using three granular materials in the laboratory.The variations of compression and shear wave velocities(Vp and Vs)are investigated during load-unload cycles under dry and drained conditions.Supplementary discrete element method(DEM)simulations are performed to understand the evolution of soil fabric during the equivalent load-unload cycles using spherical particles.Vp and Vs are not always reversible even though the stress state regains its isotropic condition after unload,indicating that Vp and Vs are governed by not only the stress state but also the fabric change.The variations of Vp/Vs are density-and stress-dependent;a higher level of stress ratio(s01/s03)threshold is observed for denser packings to trigger a significant change in wave velocity ratio(Vp/Vs)for experimental results using spherical glass beads and simulation data using spherical particles.Considering the particle shape,a higher s01/s03 threshold is found for more angular particles than rounded particles.The DEM result reveals that Vp/Vs of spherical particles can be correlated linearly with the evolution of fabric ratio(Fver/Fhor)during loadunload in a pre-peak range under dry and drained conditions.
