The performances of the cement-based materials can be improved by the incorporation of polypropylene fiber, but the damage processes become more complex with different fiber contents at the same time. The acoustic emi...The performances of the cement-based materials can be improved by the incorporation of polypropylene fiber, but the damage processes become more complex with different fiber contents at the same time. The acoustic emission(AE) technology can achieve the global monitoring of internal damage in materials. The evolution process of failure mode and damage degree of polypropylene fiber reinforced mortar and concrete were analyzed by measuring the AE energy, RA value, AF value and b value. It was found that the cement matrix cracked on the initial stage, the cracks further developed on the medium stage and the fibers were pulled out on the last stage. The matrix cracked with minor injury cracks, but the fiber broke with serious damage cracks. The cumulative AE energy was proportional to the polypropylene fiber reinforced concrete and mortar's ductility. The damage mode and damage degree can be judged by identifying the damage stage obtained by the analysis of the AF value.展开更多
This study proposes an elastic finite difference(FD)time domain method with variable grids in three-dimensional cylindrical coordinates.The calculations will diverge and become less accurate by conventional cylindrica...This study proposes an elastic finite difference(FD)time domain method with variable grids in three-dimensional cylindrical coordinates.The calculations will diverge and become less accurate by conventional cylindrical FD as the grid size gradually becomes more extensive with the increasing radius.To prevent grids from being too coarse in far fields,we compensate for the grid cell infl ation by refi ning the grid step in the azimuthal direction.The variable grid FD in the cylindrical coordinate systems has a higher effi ciency in solving acoustic logging while drilling(LWD)problems because the grid boundaries are consistent with those of the drill collar and the borehole.The proposed algorithm saves approximately 94%of the FD grids,80%of the computation time,and memory with a higher calculation accuracy than the FD on rectangular grids for the same models.We also calculate the acoustic LWD responses of the fl uid-fi lled borehole intersecting with fractures.Refl ections are generated at the fractures,which can be equivalent to an additional scattering source.The mode conversions between the collar and the Stoneley waves are revealed.The Stoneley spectra are more sensitive to the fracture.Finally,the logs in a heterogeneous formation with two refl ectors far from the borehole are modeled,and a means of estimating the azimuth of geological interfaces from refl ections is proposed.展开更多
The viscosity of fracturing fluid and in-situ stress difference are the two important factors that affect the hydraulic fracturing pressure and propagation morphology. In this study, raw coal was used to prepare coal ...The viscosity of fracturing fluid and in-situ stress difference are the two important factors that affect the hydraulic fracturing pressure and propagation morphology. In this study, raw coal was used to prepare coal samples for experiments, and clean fracturing fluid samples were prepared using CTAB surfactant. A series of hydraulic fracturing tests were conducted with an in-house developed triaxial hydraulic fracturing simulator and the fracturing process was monitored with an acoustic emission instrument to analyze the influences of fracturing fluid viscosity and horizontal in-situ stress difference on coal fracture propagation. The results show that the number of branched fractures decreased, the fracture pattern became simpler, the fractures width increased obviously, and the distribution of AE event points was concentrated with the increase of the fracturing fluid viscosity or the horizontal in-situ stress difference. The acoustic emission energy decreases with the increase of fracturing fluid viscosity and increases with the increase of horizontal in situ stress difference. The low viscosity clean fracturing fluid has strong elasticity and is easy to be compressed into the tip of fractures, resulting in complex fractures. The high viscosity clean fracturing fluids are the opposite. Our experimental results provide a reference and scientific basis for the design and optimization of field hydraulic fracturing parameters.展开更多
In this paper,we derived the relationships between the travel time difference of sPn and Pn and the local earthquake focal depth.In these equations,the travel time difference of sPn and Pn is not related to the epicen...In this paper,we derived the relationships between the travel time difference of sPn and Pn and the local earthquake focal depth.In these equations,the travel time difference of sPn and Pn is not related to the epicentral distance,but depends only on the regional crustal mode and the focal depth.According to the equations,we provided a simple and accurate method to determine local earthquake focal depth by using the travel time difference between phase sPn and Pn.This method has been used to determine the focal depths of two earthquake of MS6.1 and MS5.6 which occurred at the junction of Panzhihua and Huili,Sichuan on August 30 and 31,2008.The results were compared to those from other sources such as the China Earthquake Networks Center,and the comparison shows that the results are accurate and reliable.展开更多
Buried object which exists at extremely shallow underground can be detected by using an acoustic vibration and a SLDV (scanning laser Doppler vibrometer). Non-contact acoustic exploration can be realized by using ai...Buried object which exists at extremely shallow underground can be detected by using an acoustic vibration and a SLDV (scanning laser Doppler vibrometer). Non-contact acoustic exploration can be realized by using air coupled sound. It was confirmed that discovery and the identification of a buried thing were possible using the OFR method (optimum frequency range method). However, in this method, only the amplitude of the vibration speed spectrum is used and did not use the phase information. Therefore a new imaging method is proposed that used the phase information of the vibration speed spectrum. From the experimental results, the possibility of the outline extraction by using the phase difference between each scan point is confirmed. As a future task, the phase difference image will become clearer by changing the scan point density and long duration time of emitted wave.展开更多
Acoustic-elastic coupled media is often encountered in most marine explorations, and accurate simulation of acoustic-elastic coupled media is of great significance. At present, the study of acoustic-elastic coupled me...Acoustic-elastic coupled media is often encountered in most marine explorations, and accurate simulation of acoustic-elastic coupled media is of great significance. At present, the study of acoustic-elastic coupled media still assumes that the solid of the acoustic-elastic coupled media is isotropic, but this assumption is not in accordance with the actual situation. In this paper, we derive the solid media of acoustic-elastic coupled media from isotropic media to anisotropic media, and propose an acoustic-elastic coupled medium based ontransverse isotropic media with vertical symmetric axes(VTI) to improve the accuracy of forward modeling. Based on the relationship between the Thomsen parameter and the coefficient matrix of the anisotropic elastic wave equation, we transform the Thomson parameter into a velocity model with anisotropic properties. We use a staggered grid finite difference method to simulate the propagation of a wavefield in a three-dimensional acoustic-elastic coupled media. We obtain the snapshots of the wave field when the solid of the acoustic-elastic coupled media is an isotropic medium and a VTI media. When the solid of the acoustic-elastic coupled media is considered VTI media, we can observe the qP wave and qS wave that cannot be observed in the isotropic medium from the wave field snapshot. We can also find that the seismic records obtained by the method we use are more realistic. The algorithm proposed in this paper is of great significance for high-precision ocean numerical simulation.展开更多
Fracture identification is important for the evaluation of carbonate reservoirs. However, conventional logging equipment has small depth of investigation and cannot detect rock fractures more than three meters away fr...Fracture identification is important for the evaluation of carbonate reservoirs. However, conventional logging equipment has small depth of investigation and cannot detect rock fractures more than three meters away from the borehole. Remote acoustic logging uses phase-controlled array-transmitting and long sound probes that increase the depth of investigation. The interpretation of logging data with respect to fractures is typically guided by practical experience rather than theory and is often ambiguous. We use remote acoustic reflection logging data and high-order finite-difference approximations in the forward modeling and prestack reverse-time migration to image fractures. First, we perform forward modeling of the fracture responses as a function of the fracture-borehole wall distance, aperture, and dip angle. Second, we extract the energy intensity within the imaging area to determine whether the fracture can be identified as the formation velocity is varied. Finally, we evaluate the effect of the fracture-borehole distance, fracture aperture, and dip angle on fracture identification.展开更多
The acoustic tools widely used in borehole well logging and being developed in borehole acoustic reflection imaging do not have the function of azimuthal measurement due to a symmetric source, so they can not be used ...The acoustic tools widely used in borehole well logging and being developed in borehole acoustic reflection imaging do not have the function of azimuthal measurement due to a symmetric source, so they can not be used to evaluate the azimuthal character of borehole formation. In this paper, a 3D finite difference method was used to simulate the acoustic fields in a fluid-filled borehole generated by a traditional monopole source and a phased arc array. Acoustic waveforms were presented for both cases. The analysis of the simulated waveforms showed that different from the monopole source, the acoustic energy generated by the phased arc array transmitter mainly radiated to the borehole in a narrow azimuthal range, which was the key technique to implement azimuthal acoustic well logging. Similar to the monopole source, the waveforms generated by the phased arc array in the fluid-filled borehole also contain compressional (P) waves and shear (S) waves refracted along the borehole, which is the theoretical foundation of phased arc array acoustic well logging.展开更多
A finite difference/boundary integral procedure to determine theacoustic reflected pressure from a fluid-loaded bi-laminated plate isdescribed. The bi-laminate is composed of a piezoelectric layer andan elastic layer ...A finite difference/boundary integral procedure to determine theacoustic reflected pressure from a fluid-loaded bi-laminated plate isdescribed. The bi-laminate is composed of a piezoelectric layer andan elastic layer in contact with the fluid, and is held by anacoustically hard baffle. In the numerical model, the fluid pressureat fluid/solid interface is replaced by a continuum of point sourcesweighted by the normal acceleration of the elastic plate, and thegoverning equation system is solved in the solid domain.展开更多
The acoustic reflected pressure from a periodic elastic/piezoelectric laminated plate is studied for the purpose of acoustic reflection control.A finite difference/boundary integral procedure to determine the reflecte...The acoustic reflected pressure from a periodic elastic/piezoelectric laminated plate is studied for the purpose of acoustic reflection control.A finite difference/boundary integral procedure to determine the reflected pressure from the fluid-loaded plate is described.In the numerical model,a Green's function in the form of infinite sum is employed and a boundary integral is performed to replace the fluid pressure at fluid/solid interface by a continuum of point sources weighted by the normal acceleration of the elastic plate.The equation system is then solved only in the solid domain.It is demonstrated that an appropriate applied voltage potential across the piezoelectric layer has the effect of cancelling the fundamental propagating mode,and there,is no reflection for frequencies up to the cut-off frequency of the next propagating mode if the fundamental mode has been eliminated.展开更多
Full-waveform velocity inversion based on the acoustic wave equation in the time domain is investigated in this paper. The inversion is the iterative minimization of the misfit between observed data and synthetic data...Full-waveform velocity inversion based on the acoustic wave equation in the time domain is investigated in this paper. The inversion is the iterative minimization of the misfit between observed data and synthetic data obtained by a numerical solution of the wave equation. Two inversion algorithms in combination with the CG method and the BFGS method are described respectively. Numerical computations for two models including the benchmark Marmousi model with complex structure are implemented. The inversion results show that the BFGS-based algorithm behaves better in inversion than the CG-based algorithm does. Moreover, the good inversion result for Marmousi model with the BFGS-based algorithm suggests the quasi-Newton methods can provide an important tool for large-scale velocity inversion. More computations demonstrate the correctness and effectives of our inversion algorithms and code.展开更多
The inhomogenous ocean waveguide,which leads the amplitude and phase of the signal arriving at a hydrophone array to fluctuate,is one of the causes that make the array gain deviate from its ideal value.The relationshi...The inhomogenous ocean waveguide,which leads the amplitude and phase of the signal arriving at a hydrophone array to fluctuate,is one of the causes that make the array gain deviate from its ideal value.The relationship between the array gain and the fluctuant acoustic channel is studied theoretically.The analytical expression of the array gain is derived via an acoustic channel transfer function on the assumption that the ambient noise field is isotropic.The expression is expanded via the Euler formula to give an insight into the effect of the fluctuant acoustic channel on the array gain.The result demonstrates that the amplitude fluctuation of the acoustic channel transfer functions has a slight effect on the array gain;however,the uniformity of the phase difference between the weighting coefficient and the channel transfer function on all the hydrophones in the array is a major factor that leads the array gain to further deviate from its ideal value.The numerical verification is conducted in the downslope waveguide,in which the gain of a horizontal uniform linear array(HLA)with a wide-aperture operating in the continental slope area is considered.Numerical result is consistent with the theoretical analysis.展开更多
During earthquakes and strong underground explosions it is possible to observe two different effects. The first one is connected with the destruction of media, and this causes acoustic and later hybrid acoustic-electr...During earthquakes and strong underground explosions it is possible to observe two different effects. The first one is connected with the destruction of media, and this causes acoustic and later hybrid acoustic-electromagnetic waves in an epicenter in the atmosphere and in the ionosphere. Another one is connected with cracks in crystals of rocks, which seems more interesting, because it is possible to recognize the natural earthquakes and industrial explosions. In the first effects due to nonlinear elastic modules the acoustic waves move through the lithosphere and transform their spectra from VLF (very low frequencies ~ 1 - 10 kHz) at the depth of about 30 km into the lower part of ELF (extremely low frequencies, ~ 3 Hz - 1 kHz) on the Earth’s surface, then they pass the atmosphere and penetrate into the ionosphere. During the nonlinear acoustic passage through the atmosphere and the ionosphere, the spectrum transforms from ELF into ULF (ultra low frequencies, < 1 Hz) range. In this review article the classification of spectrum and analysis of two cases of the destruction of rocks in the lithosphere is presented. The rocks possess piezoelectric and piezomagnetic properties. In this case the electromagnetic emission is excited by the fracturing in plates of crystals. The difference of emission from piezoelectric and magnetic plates in cases of industrial explosions and natural seismic events including volcanic phenomena gives a possibility to analyze the method of its identification. The consideration is based on the model of the plate of a finite size with an uniformly moving crack.展开更多
In 3D frequency domain seismic forward and inversion calculation,the huge amount of calculation and storage is one of the main factors that restrict the processing speed and calculation efficiency.The frequency domain...In 3D frequency domain seismic forward and inversion calculation,the huge amount of calculation and storage is one of the main factors that restrict the processing speed and calculation efficiency.The frequency domain finite-difference forward simulation algorithm based on the acoustic wave equation establishes a large bandwidth complex matrix according to the discretized acoustic wave equation,and then the frequency domain wave field value is obtained by solving the matrix equation.In this study,the predecessor's optimized five-point method is extended to a 3D seven-point finite-difference scheme,and then a perfectly matched layer absorbing boundary condition(PML)is added to establish the corresponding matrix equation.In order to solve the complex matrix,we transform it to the equivalent real number domain to expand the solvable range of the matrix,and establish two objective functions to transform the matrix solving problem into an optimization problem that can be solved using gradient methods,and then use conjugate gradient algorithm to solve the problem.Previous studies have shown that in the conjugate gradient algorithm,the product of the matrix and the vector is the main factor that affects the calculation efficiency.Therefore,this study proposes a method that transform bandwidth matrix and vector product problem into some equivalent vector and vector product algorithm,thereby reducing the amount of calculation and storage.展开更多
A numerical-analytical solution for seismic and acoustic-gravity waves propagation is applied to a heterogeneous “Earth-Atmosphere” model. Seismic wave propagation in an elastic half-space is described by a system o...A numerical-analytical solution for seismic and acoustic-gravity waves propagation is applied to a heterogeneous “Earth-Atmosphere” model. Seismic wave propagation in an elastic half-space is described by a system of first order dynamic equations of elasticity theory. Propagation of acoustic-gravity waves in the atmosphere is described by the linearized Navier-Stokes equations with the wind. The algorithm proposed is based on the integral Laguerre transform with respect to time, the finite integral Fourier transform along the spatial coordinate with the finite difference solution of the reduced problem.展开更多
基金Funded by National Natural Science Foundation of China(No.51009058)Postdoctoral Science Foundation of China(No.2011M501160)
文摘The performances of the cement-based materials can be improved by the incorporation of polypropylene fiber, but the damage processes become more complex with different fiber contents at the same time. The acoustic emission(AE) technology can achieve the global monitoring of internal damage in materials. The evolution process of failure mode and damage degree of polypropylene fiber reinforced mortar and concrete were analyzed by measuring the AE energy, RA value, AF value and b value. It was found that the cement matrix cracked on the initial stage, the cracks further developed on the medium stage and the fibers were pulled out on the last stage. The matrix cracked with minor injury cracks, but the fiber broke with serious damage cracks. The cumulative AE energy was proportional to the polypropylene fiber reinforced concrete and mortar's ductility. The damage mode and damage degree can be judged by identifying the damage stage obtained by the analysis of the AF value.
基金supported by the National Natural Science Foundation of China(Grant Nos.12174421,11774373,11734017,and 42074215).
文摘This study proposes an elastic finite difference(FD)time domain method with variable grids in three-dimensional cylindrical coordinates.The calculations will diverge and become less accurate by conventional cylindrical FD as the grid size gradually becomes more extensive with the increasing radius.To prevent grids from being too coarse in far fields,we compensate for the grid cell infl ation by refi ning the grid step in the azimuthal direction.The variable grid FD in the cylindrical coordinate systems has a higher effi ciency in solving acoustic logging while drilling(LWD)problems because the grid boundaries are consistent with those of the drill collar and the borehole.The proposed algorithm saves approximately 94%of the FD grids,80%of the computation time,and memory with a higher calculation accuracy than the FD on rectangular grids for the same models.We also calculate the acoustic LWD responses of the fl uid-fi lled borehole intersecting with fractures.Refl ections are generated at the fractures,which can be equivalent to an additional scattering source.The mode conversions between the collar and the Stoneley waves are revealed.The Stoneley spectra are more sensitive to the fracture.Finally,the logs in a heterogeneous formation with two refl ectors far from the borehole are modeled,and a means of estimating the azimuth of geological interfaces from refl ections is proposed.
基金National Natural Science Foundation of China (51974176, 52174194, 51934004)Shandong Provincial Colleges and Universities Youth Innovation and Technology Support Program (2019KJH006)+1 种基金Taishan Scholars Project (TS20190935)Shandong outstanding youth fund (ZR2020JQ22).
文摘The viscosity of fracturing fluid and in-situ stress difference are the two important factors that affect the hydraulic fracturing pressure and propagation morphology. In this study, raw coal was used to prepare coal samples for experiments, and clean fracturing fluid samples were prepared using CTAB surfactant. A series of hydraulic fracturing tests were conducted with an in-house developed triaxial hydraulic fracturing simulator and the fracturing process was monitored with an acoustic emission instrument to analyze the influences of fracturing fluid viscosity and horizontal in-situ stress difference on coal fracture propagation. The results show that the number of branched fractures decreased, the fracture pattern became simpler, the fractures width increased obviously, and the distribution of AE event points was concentrated with the increase of the fracturing fluid viscosity or the horizontal in-situ stress difference. The acoustic emission energy decreases with the increase of fracturing fluid viscosity and increases with the increase of horizontal in situ stress difference. The low viscosity clean fracturing fluid has strong elasticity and is easy to be compressed into the tip of fractures, resulting in complex fractures. The high viscosity clean fracturing fluids are the opposite. Our experimental results provide a reference and scientific basis for the design and optimization of field hydraulic fracturing parameters.
基金funded by the special support projectentitled "Sorting out and processing of seismic data " of central public-interest basic scientific and technological research of Institute of Crustal DynamicsChina Earthquake Administration (ZDJ2007-4)
文摘In this paper,we derived the relationships between the travel time difference of sPn and Pn and the local earthquake focal depth.In these equations,the travel time difference of sPn and Pn is not related to the epicentral distance,but depends only on the regional crustal mode and the focal depth.According to the equations,we provided a simple and accurate method to determine local earthquake focal depth by using the travel time difference between phase sPn and Pn.This method has been used to determine the focal depths of two earthquake of MS6.1 and MS5.6 which occurred at the junction of Panzhihua and Huili,Sichuan on August 30 and 31,2008.The results were compared to those from other sources such as the China Earthquake Networks Center,and the comparison shows that the results are accurate and reliable.
文摘Buried object which exists at extremely shallow underground can be detected by using an acoustic vibration and a SLDV (scanning laser Doppler vibrometer). Non-contact acoustic exploration can be realized by using air coupled sound. It was confirmed that discovery and the identification of a buried thing were possible using the OFR method (optimum frequency range method). However, in this method, only the amplitude of the vibration speed spectrum is used and did not use the phase information. Therefore a new imaging method is proposed that used the phase information of the vibration speed spectrum. From the experimental results, the possibility of the outline extraction by using the phase difference between each scan point is confirmed. As a future task, the phase difference image will become clearer by changing the scan point density and long duration time of emitted wave.
基金Supported by Major Project of National Science and Technology of China(No.2016ZX05026-002-003)National Natural Science Foundation of China(No.41374108)
文摘Acoustic-elastic coupled media is often encountered in most marine explorations, and accurate simulation of acoustic-elastic coupled media is of great significance. At present, the study of acoustic-elastic coupled media still assumes that the solid of the acoustic-elastic coupled media is isotropic, but this assumption is not in accordance with the actual situation. In this paper, we derive the solid media of acoustic-elastic coupled media from isotropic media to anisotropic media, and propose an acoustic-elastic coupled medium based ontransverse isotropic media with vertical symmetric axes(VTI) to improve the accuracy of forward modeling. Based on the relationship between the Thomsen parameter and the coefficient matrix of the anisotropic elastic wave equation, we transform the Thomson parameter into a velocity model with anisotropic properties. We use a staggered grid finite difference method to simulate the propagation of a wavefield in a three-dimensional acoustic-elastic coupled media. We obtain the snapshots of the wave field when the solid of the acoustic-elastic coupled media is an isotropic medium and a VTI media. When the solid of the acoustic-elastic coupled media is considered VTI media, we can observe the qP wave and qS wave that cannot be observed in the isotropic medium from the wave field snapshot. We can also find that the seismic records obtained by the method we use are more realistic. The algorithm proposed in this paper is of great significance for high-precision ocean numerical simulation.
基金supported by National Petroleum Major Project(Grant No.2011ZX05020-008)
文摘Fracture identification is important for the evaluation of carbonate reservoirs. However, conventional logging equipment has small depth of investigation and cannot detect rock fractures more than three meters away from the borehole. Remote acoustic logging uses phase-controlled array-transmitting and long sound probes that increase the depth of investigation. The interpretation of logging data with respect to fractures is typically guided by practical experience rather than theory and is often ambiguous. We use remote acoustic reflection logging data and high-order finite-difference approximations in the forward modeling and prestack reverse-time migration to image fractures. First, we perform forward modeling of the fracture responses as a function of the fracture-borehole wall distance, aperture, and dip angle. Second, we extract the energy intensity within the imaging area to determine whether the fracture can be identified as the formation velocity is varied. Finally, we evaluate the effect of the fracture-borehole distance, fracture aperture, and dip angle on fracture identification.
基金supported by the National Natural ScienceFoundation of China(Grant Nos.10534040,40574049 and 40874097)the Research Fund for the Doctoral Programof Higher Education(Grant No.20070425028)the Foundation of State Key Laboratory of Petroleum Resourceand Prospecting,China University of Petroleum(Grant No.PRPDX2008-08).
文摘The acoustic tools widely used in borehole well logging and being developed in borehole acoustic reflection imaging do not have the function of azimuthal measurement due to a symmetric source, so they can not be used to evaluate the azimuthal character of borehole formation. In this paper, a 3D finite difference method was used to simulate the acoustic fields in a fluid-filled borehole generated by a traditional monopole source and a phased arc array. Acoustic waveforms were presented for both cases. The analysis of the simulated waveforms showed that different from the monopole source, the acoustic energy generated by the phased arc array transmitter mainly radiated to the borehole in a narrow azimuthal range, which was the key technique to implement azimuthal acoustic well logging. Similar to the monopole source, the waveforms generated by the phased arc array in the fluid-filled borehole also contain compressional (P) waves and shear (S) waves refracted along the borehole, which is the theoretical foundation of phased arc array acoustic well logging.
基金the National Natural Science Foundation of China (No.10172039).
文摘A finite difference/boundary integral procedure to determine theacoustic reflected pressure from a fluid-loaded bi-laminated plate isdescribed. The bi-laminate is composed of a piezoelectric layer andan elastic layer in contact with the fluid, and is held by anacoustically hard baffle. In the numerical model, the fluid pressureat fluid/solid interface is replaced by a continuum of point sourcesweighted by the normal acceleration of the elastic plate, and thegoverning equation system is solved in the solid domain.
基金Project supported by the National Natural Science Foundation of China(No.10172039).
文摘The acoustic reflected pressure from a periodic elastic/piezoelectric laminated plate is studied for the purpose of acoustic reflection control.A finite difference/boundary integral procedure to determine the reflected pressure from the fluid-loaded plate is described.In the numerical model,a Green's function in the form of infinite sum is employed and a boundary integral is performed to replace the fluid pressure at fluid/solid interface by a continuum of point sources weighted by the normal acceleration of the elastic plate.The equation system is then solved only in the solid domain.It is demonstrated that an appropriate applied voltage potential across the piezoelectric layer has the effect of cancelling the fundamental propagating mode,and there,is no reflection for frequencies up to the cut-off frequency of the next propagating mode if the fundamental mode has been eliminated.
文摘Full-waveform velocity inversion based on the acoustic wave equation in the time domain is investigated in this paper. The inversion is the iterative minimization of the misfit between observed data and synthetic data obtained by a numerical solution of the wave equation. Two inversion algorithms in combination with the CG method and the BFGS method are described respectively. Numerical computations for two models including the benchmark Marmousi model with complex structure are implemented. The inversion results show that the BFGS-based algorithm behaves better in inversion than the CG-based algorithm does. Moreover, the good inversion result for Marmousi model with the BFGS-based algorithm suggests the quasi-Newton methods can provide an important tool for large-scale velocity inversion. More computations demonstrate the correctness and effectives of our inversion algorithms and code.
基金Project supported by the National Natural Science Foundation of China(Grant No.11534009)
文摘The inhomogenous ocean waveguide,which leads the amplitude and phase of the signal arriving at a hydrophone array to fluctuate,is one of the causes that make the array gain deviate from its ideal value.The relationship between the array gain and the fluctuant acoustic channel is studied theoretically.The analytical expression of the array gain is derived via an acoustic channel transfer function on the assumption that the ambient noise field is isotropic.The expression is expanded via the Euler formula to give an insight into the effect of the fluctuant acoustic channel on the array gain.The result demonstrates that the amplitude fluctuation of the acoustic channel transfer functions has a slight effect on the array gain;however,the uniformity of the phase difference between the weighting coefficient and the channel transfer function on all the hydrophones in the array is a major factor that leads the array gain to further deviate from its ideal value.The numerical verification is conducted in the downslope waveguide,in which the gain of a horizontal uniform linear array(HLA)with a wide-aperture operating in the continental slope area is considered.Numerical result is consistent with the theoretical analysis.
文摘During earthquakes and strong underground explosions it is possible to observe two different effects. The first one is connected with the destruction of media, and this causes acoustic and later hybrid acoustic-electromagnetic waves in an epicenter in the atmosphere and in the ionosphere. Another one is connected with cracks in crystals of rocks, which seems more interesting, because it is possible to recognize the natural earthquakes and industrial explosions. In the first effects due to nonlinear elastic modules the acoustic waves move through the lithosphere and transform their spectra from VLF (very low frequencies ~ 1 - 10 kHz) at the depth of about 30 km into the lower part of ELF (extremely low frequencies, ~ 3 Hz - 1 kHz) on the Earth’s surface, then they pass the atmosphere and penetrate into the ionosphere. During the nonlinear acoustic passage through the atmosphere and the ionosphere, the spectrum transforms from ELF into ULF (ultra low frequencies, < 1 Hz) range. In this review article the classification of spectrum and analysis of two cases of the destruction of rocks in the lithosphere is presented. The rocks possess piezoelectric and piezomagnetic properties. In this case the electromagnetic emission is excited by the fracturing in plates of crystals. The difference of emission from piezoelectric and magnetic plates in cases of industrial explosions and natural seismic events including volcanic phenomena gives a possibility to analyze the method of its identification. The consideration is based on the model of the plate of a finite size with an uniformly moving crack.
基金supported by the National Natural Science Foundation of China(Project U1901602&41790465)Key Special Project for Introduced Talents Team of Southern Marine Science and Engineering Guangdong Laboratory(Guangzhou)(GML2019ZD0203)+2 种基金Shenzhen Key Laboratory of Deep Offshore Oil and Gas Exploration Technology(Grant No.ZDSYS20190902093007855)Shenzhen Science and Technology Program(Grant No.KQTD20170810111725321)the leading talents of Guangdong province program(Grant No.2016LJ06N652).
文摘In 3D frequency domain seismic forward and inversion calculation,the huge amount of calculation and storage is one of the main factors that restrict the processing speed and calculation efficiency.The frequency domain finite-difference forward simulation algorithm based on the acoustic wave equation establishes a large bandwidth complex matrix according to the discretized acoustic wave equation,and then the frequency domain wave field value is obtained by solving the matrix equation.In this study,the predecessor's optimized five-point method is extended to a 3D seven-point finite-difference scheme,and then a perfectly matched layer absorbing boundary condition(PML)is added to establish the corresponding matrix equation.In order to solve the complex matrix,we transform it to the equivalent real number domain to expand the solvable range of the matrix,and establish two objective functions to transform the matrix solving problem into an optimization problem that can be solved using gradient methods,and then use conjugate gradient algorithm to solve the problem.Previous studies have shown that in the conjugate gradient algorithm,the product of the matrix and the vector is the main factor that affects the calculation efficiency.Therefore,this study proposes a method that transform bandwidth matrix and vector product problem into some equivalent vector and vector product algorithm,thereby reducing the amount of calculation and storage.
文摘A numerical-analytical solution for seismic and acoustic-gravity waves propagation is applied to a heterogeneous “Earth-Atmosphere” model. Seismic wave propagation in an elastic half-space is described by a system of first order dynamic equations of elasticity theory. Propagation of acoustic-gravity waves in the atmosphere is described by the linearized Navier-Stokes equations with the wind. The algorithm proposed is based on the integral Laguerre transform with respect to time, the finite integral Fourier transform along the spatial coordinate with the finite difference solution of the reduced problem.