Severe vibration of underground structures may be induced under blast loads. According to the characteristics of the explosion-induced ground shock wave, a new-type damper, inverse control magneto-rheological（MR） da...Severe vibration of underground structures may be induced under blast loads. According to the characteristics of the explosion-induced ground shock wave, a new-type damper, inverse control magneto-rheological（MR） damper was designed to control the vibration, The high-frequency performance test of the MR damper was carried out on the small shaking table. It is shown that the performance can be modeled by use of the modified Bouc-Wen model, and the Parameters of the model keep stable in the range of 15--50 Hz.展开更多
This study is an exploratory analysis of applying natural language processing techniques such as Term Frequency-Inverse Document Frequency and Sentiment Analysis on Twitter data. The uniqueness of this work is establi...This study is an exploratory analysis of applying natural language processing techniques such as Term Frequency-Inverse Document Frequency and Sentiment Analysis on Twitter data. The uniqueness of this work is established by determining the overall sentiment of a politician’s tweets based on TF-IDF values of terms used in their published tweets. By calculating the TF-IDF value of terms from the corpus, this work displays the correlation between TF-IDF score and polarity. The results of this work show that calculating the TF-IDF score of the corpus allows for a more accurate representation of the overall polarity since terms are given a weight based on their uniqueness and relevance rather than just the frequency at which they appear in the corpus.展开更多
Recently,the metasurfaces for independently controlling the wavefront and amplitude of two orthogonal circularly polarized electromagnetic(EM)waves have been demonstrated to open a way toward spin-multiplexing compact...Recently,the metasurfaces for independently controlling the wavefront and amplitude of two orthogonal circularly polarized electromagnetic(EM)waves have been demonstrated to open a way toward spin-multiplexing compact metadevices.However,these metasurfaces are mostly restricted to a single operation frequency band.The main challenge to achieving multiple frequency manipulations stems from the complicated and time-consuming design caused by multifrequency cross talk.To solve this problem,we propose a deep-learning-assisted inverse design method for designing a dual-spin/frequency metasurface with flexible multiplexing of off-axis vortices.By analyzing the cross talk between different spin/frequency channels based on the deep-learning method,we established the internal mapping relationship between the physical parameters of a meta-atom and its phase responses in multichannels,realizing the rapid inverse design of the spin/frequency multiplexing EM device.As a proof of concept,we demonstrated in the microwave region a dual-frequency arbitrary spin-to-orbit angular momentum converter,a dual-frequency off-axis vector vortex multiplexer,and a large-capacity(16-channel)vortex beam generator.The proposed method may provide a compact and efficient platform for the multiplexing of vortices,which may further stimulate their applications in wireless communication and quantum information science.展开更多
Three-dimensional gravity inversion based on the mass property model is very popular in recent years. The time and efficiency of inversion algorithms is relative to the magnitude of the target mesh. One approach is to...Three-dimensional gravity inversion based on the mass property model is very popular in recent years. The time and efficiency of inversion algorithms is relative to the magnitude of the target mesh. One approach is to search over the entire solution space for a more refined result. However, the inversion will be difficult with the increased parameters in the large search space and the number of computations increases exponentially. ｜n this paper, we propose a novel approach based on the frequency characteristics of the density distribution over the mesh. The purposes of our study are to reduce the parameters of three- dimensional gravity inversion and to lighten the image quality of the inversion result. The results show that the new method can expedite the inversion processing and get a better geological interpretation than tradition methods.展开更多
Although full waveform inversion in the frequency domain can overcome the local minima problem in the time direction, such problem still exists in the space direction because of the media subsurface complexity. Based ...Although full waveform inversion in the frequency domain can overcome the local minima problem in the time direction, such problem still exists in the space direction because of the media subsurface complexity. Based on the optimal steep descent methods, we present an algorithm which combines the preconditioned bi-conjugated gradient stable method and the multi-grid method to compute the wave propagation and the gradient space. The multiple scale prosperity of the waveform inversion and the multi-grid method can overcome the inverse problems local minima defect and accelerate convergence. The local inhomogeneous three-hole model simulated results and the Marmousi model certify the algorithm effectiveness.展开更多
We studied finite-element-method-based two-dimensional frequency-domain acoustic FWI under rugged topography conditions. The exponential attenuation boundary condition suitable for rugged topography is proposed to sol...We studied finite-element-method-based two-dimensional frequency-domain acoustic FWI under rugged topography conditions. The exponential attenuation boundary condition suitable for rugged topography is proposed to solve the cutoff botmdary problem as well as to consider the requirement of using the same subdivision grid in joint multifrequency inversion. The proposed method introduces the attenuation factor, and by adjusting it, acoustic waves are sufficiently attenuated in the attenuation layer to minimize the cutoff boundary effect. Based on the law of exponential attenuation, expressions for computing the attenuation factor and the thickness of attenuation layers are derived for different frequencies. In multifrequency-domain FWI, the conjugate gradient method is used to solve equations in the Gauss-Newton algorithm and thus minimize the computation cost in calculating the Hessian matrix. In addition, the effect of initial model selection and frequency combination on FWI is analyzed. Examples using numerical simulations and FWI calculations are used to verify the efficiency of the proposed method.展开更多
The complexity of an elastic wavefield increases the nonlinearity of inversion, To some extent, multiscale inversion decreases the nonlinearity of inversion and prevents it from falling into local extremes. A multisca...The complexity of an elastic wavefield increases the nonlinearity of inversion, To some extent, multiscale inversion decreases the nonlinearity of inversion and prevents it from falling into local extremes. A multiscale strategy based on the simultaneous use of frequency groups and layer stripping method based on damped wave field improves the stability of inversion. A dual-level parallel algorithm is then used to decrease the computational cost and improve practicability. The seismic wave modeling of a single frequency and inversion in a frequency group are computed in parallel by multiple nodes based on multifrontal massively parallel sparse direct solver and MPI. Numerical tests using an overthrust model show that the proposed inversion algorithm can effectively improve the stability and accuracy of inversion by selecting the appropriate inversion frequency and damping factor in low- frequency seismic data.展开更多
As a high quality seismic imaging method, full waveform inversion （FWI） can accurately reconstruct the physical parameter model for the subsurface medium. However, application of the FWI in seismic data processing i...As a high quality seismic imaging method, full waveform inversion （FWI） can accurately reconstruct the physical parameter model for the subsurface medium. However, application of the FWI in seismic data processing is computationally expensive, especially for the three-dimension complex medium inversion. Introducing blended source technology into the frequency-domain FWI can greatly reduce the computational burden and improve the efficiency of the inversion. However, this method has two issues： first, crosstalk noise is caused by interference between the sources involved in the encoding, resulting in an inversion result with some artifacts; second, it is more sensitive to ambient noise compared to conventional FWI, therefore noisy data results in a poor inversion. This paper introduces a frequency-group encoding method to suppress crosstalk noise, and presents a frequency- domain auto-adapting FWI based on source-encoding technology. The conventional FWI method and source-encoding based FWI method are combined using an auto-adapting mechanism. This improvement can both guarantee the quality of the inversion result and maximize the inversion efficiency.展开更多
Amplitude variations with offset or incident angle （AVO/AVA） inversion are typically combined with statistical methods, such as Bayesian inference or deterministic inversion. We propose a joint elastic inversion met...Amplitude variations with offset or incident angle （AVO/AVA） inversion are typically combined with statistical methods, such as Bayesian inference or deterministic inversion. We propose a joint elastic inversion method in the time and frequency domain based on Bayesian inversion theory to improve the resolution of the estimated P- and S-wave velocities and density. We initially construct the objective function using Bayesian inference by combining seismic data in the time and frequency domain. We use Cauchy and Gaussian probability distribution density functions to obtain the prior information for the model parameters and the likelihood function, respectively. We estimate the elastic parameters by solving the initial objective function with added model constraints to improve the inversion robustness. The results of the synthetic data suggest that the frequency spectra of the estimated parameters are wider than those obtained with conventional elastic inversion in the time domain. In addition, the proposed inversion approach offers stronger antinoising compared to the inversion approach in the frequency domain. Furthermore, results from synthetic examples with added Gaussian noise demonstrate the robustness of the proposed approach. From the real data, we infer that more model parameter details can be reproduced with the proposed joint elastic inversion.展开更多
Short period surface waves generated by a local earthquake recorded by broadband seismometers at distances of about 186 to 778 km from the earthquake’s epicenter located in Cameroon (Central Africa) were processed fo...Short period surface waves generated by a local earthquake recorded by broadband seismometers at distances of about 186 to 778 km from the earthquake’s epicenter located in Cameroon (Central Africa) were processed for group velocity maps and dispersion waveforms using the frequency time analysis (FTAN) method. The resulting group velocity fundamental modes of the extracted Rayleigh and Love waves were used for a joint amplitude spectral and P polarity inversion using moment tensor inversion. The corresponding group velocity dispersion curves, the residual as a function of depth, the amplitude spectra and the moment tensor solutions of the regions from the epicenter to the different stations up to a depth of about 10 km were obtained.展开更多
Tunnel seismic detection methods are effective for obtaining the geological structure around the tunnel face,which is critical for safe construction and disaster mitigation in tunnel engineering.However,there is often...Tunnel seismic detection methods are effective for obtaining the geological structure around the tunnel face,which is critical for safe construction and disaster mitigation in tunnel engineering.However,there is often a lack of accuracy in the acquired geological information and physical properties ahead of the tunnel face in the current tunnel seismic detection methods.Thus,we apply a frequency-domain acoustic full-waveform inversion(FWI)method to obtain high-resolution results for the tunnel structure.We discuss the influence of the frequency group selection strategy and the tunnel observation system settings regarding the inversion results and determine the structural imaging and physical property parameter inversion of abnormal geological bodies ahead of the tunnel face.Based on the conventional strategies of frequency-domain acoustic FWI,we propose a frequency group selection strategy that combines a low-frequency selection covering the vertical wavenumber and a high-frequency selection of antialiasing.This strategy can effectively obtain the spatial structure and physical parameters of the geology ahead of the tunnel face and improve the inversion resolution.In addition,by linearly increasing the side length of the tunnel observation system,we share the influence of the length of the two sides of the observation systems of different tunnels on the inversion results.We found out that the inversion results are the best when the side length is approximately five times the width of the tunnel face,and the influence of increasing the side observation length beyond this range on the inversion results can be ignored.Finally,based on this approach,we invert for the complex multi-stratum model,and an accurate structure and physical property parameters of the complex stratum ahead of the tunnel face are obtained,which verifies the feasibility of the proposed method.展开更多
We propose a new 3D inversion scheme to invert the near-and transition-zone data of CSAMT with topography accurately.In this new method,the earth was discretized into unstructured tetrahedra to fit the ragged topograp...We propose a new 3D inversion scheme to invert the near-and transition-zone data of CSAMT with topography accurately.In this new method,the earth was discretized into unstructured tetrahedra to fit the ragged topography and the vector fi nite-element method was adopted to obtain precise responses and good sensitivity.To simulate the attitude and shape of the transmitter,we divided a long-grounded transmitter into dipoles and integrated these dipoles to obtain good responses in the near-and transition-fi eld zones.Next,we designed an L2 norm-based objective functional and applied a standard quasi-Newton method as the optimization method to solve the inverse problem and guarantee steady convergence.We tested our 3D inversion method first on synthetic data and then on a field dataset acquired from select sites near Changbai Mountain,China.In both tests,the new inversion algorithm achieved excellent fitting between the predicted and observed data,even in near-and transition-fi eld zones,and the inversion results agreed well with the true model.These fi ndings reveal that the proposed algorithm is eff ective for 3D inversion of CSAMT data.展开更多
The controlled source extremely low frequency （CSELF） electromagnetic method is characterized by extremely long and powerful sources and a huge measurement range. Its electromagnetic field can therefore be affected ...The controlled source extremely low frequency （CSELF） electromagnetic method is characterized by extremely long and powerful sources and a huge measurement range. Its electromagnetic field can therefore be affected by the ionosphere and displacement current. Research on 3D forward modeling and inversion of CSELF electromagnetic data is currently in its infancy. This paper makes exploratory attempts to firstly calculate the 1D extremely low frequency electromagnetic field under ionosphere-air-earth coupling circumstances, and secondly analyze the propagation characteristics of the background electromagnetic field. The 3D staggered-grid finite difference scheme for solving for the secondary electric field is adopted and incorporated with the 1D modeling algorithm to complete 3D forward modeling. Considering that surveys can be carried out in the near field and transition zone for lower frequencies, the 3D Limited-memory Broyden-Fletcher-Goldfarb-Shanno （LBFGS） inversion of CSELF electromagnetic data is presented （in which the sources, or primary fields, are included）, with the aim of directly inverting the impedance data, regardless of where it is acquired. Derivation of the objective functional gradient is the core component in the inversion. Synthetic tests indicate that the well-chosen approximation to the Hessian can significantly speed up the inversion. The model responses corresponding to the coexistence of conductive and resistive blocks show that the off-diagonal components of tensor impedance are much more sensitive to the resistivity variation than the diagonal components. In comparison with conventional scalar inversion, tensor inversion is superior in the recoveries of electric anomalies and background resistivity.展开更多
Finding materials with specific properties is a hot topic in materials science.Traditional materials design relies on empirical and trial-and-error methods,requiring extensive experiments and time,resulting in high co...Finding materials with specific properties is a hot topic in materials science.Traditional materials design relies on empirical and trial-and-error methods,requiring extensive experiments and time,resulting in high costs.With the development of physics,statistics,computer science,and other fields,machine learning offers opportunities for systematically discovering new materials.Especially through machine learning-based inverse design,machine learning algorithms analyze the mapping relationships between materials and their properties to find materials with desired properties.This paper first outlines the basic concepts of materials inverse design and the challenges faced by machine learning-based approaches to materials inverse design.Then,three main inverse design methods—exploration-based,model-based,and optimization-based—are analyzed in the context of different application scenarios.Finally,the applications of inverse design methods in alloys,optical materials,and acoustic materials are elaborated on,and the prospects for materials inverse design are discussed.The authors hope to accelerate the discovery of new materials and provide new possibilities for advancing materials science and innovative design methods.展开更多
Multi-wave exploration is an effective means for improving precision in the exploration and development of complex oil and gas reservoirs that are dense and have low permeability. However, convened wave data is charac...Multi-wave exploration is an effective means for improving precision in the exploration and development of complex oil and gas reservoirs that are dense and have low permeability. However, convened wave data is characterized by a low signal-to-noise ratio and low resolution, because the conventional deconvolution technology is easily affected by the frequency range limits, and there is limited scope for improving its resolution. The spectral inversion techniques is used to identify λ/8 thin layers and its breakthrough regarding band range limits has greatly improved the seismic resolution. The difficulty associated with this technology is how to use the stable inversion algorithm to obtain a high-precision reflection coefficient, and then to use this reflection coefficient to reconstruct broadband data for processing. In this paper, we focus on how to improve the vertical resolution of the converted PS-wave for multi-wave data processing. Based on previous research, we propose a least squares inversion algorithm with a total variation constraint, in which we uses the total variance as a priori information to solve under-determined problems, thereby improving the accuracy and stability of the inversion. Here, we simulate the Gaussian fitting amplitude spectrum to obtain broadband wavelet data, which we then process to obtain a higher resolution converted wave. We successfully apply the proposed inversion technology in the processing of high-resolution data from the Penglai region to obtain higher resolution convened wave data, which we then verify in a theoretical test. Improving the resolution of converted PS-wave data will provide more accurate data for subsequent velocity inversion and the extraction of reservoir reflection information.展开更多
Fluid and effective fracture identification in reservoirs is a crucial part of reservoir prediction.The frequency-dependent AVO inversion algorithms have proven to be effective for identifying fluid through its disper...Fluid and effective fracture identification in reservoirs is a crucial part of reservoir prediction.The frequency-dependent AVO inversion algorithms have proven to be effective for identifying fluid through its dispersion property.However,the conventional frequency-dependent AVO inversion algorithms based on Smith&Gidlow and Aki&Richards approximations do not consider the acquisition azimuth of seismic data and neglect the effect of seismic anisotropic dispersion in the actual medium.The aligned fractures in the subsurface medium induce anisotropy.The seismic anisotropy should be considered while accounting for the seismic dispersion properties through fluid-saturated fractured reservoirs.Anisotropy in such reservoirs is frequency-related due to wave-induced fluid-flow(WIFF)between interconnected fractures and pores.It can be used to identify fluid and effective fractures(fluid-saturated)by using azimuthal seismic data via anisotropic dispersion properties.In this paper,based on Rüger’s equation,we derived an analytical expression in the frequency domain for the frequencydependent AVOAz inversion in terms of fracture orientation,dispersion gradient of isotropic background rock,anisotropic dispersion gradient,and the dispersion at a normal incident angle.The frequency-dependent AVOAz equation utilizes azimuthal seismic data and considers the effect of both isotropic and anisotropic dispersion.Reassigned Gabor Transform(RGT)is used to achieve highresolution frequency division data.We then propose the frequency-dependent AVOAz inversion method to identify fluid and characterize effective fractures in fractured porous reservoirs.Through application to high-qualified seismic data of dolomite and carbonate reservoirs,the results show that the method is useful for identifying fluid and effective fractures in fluid-saturated fractured rocks.展开更多
Frequency-dependent amplitude versus offset(FAVO)inversion is a popular method to estimate the frequency-dependent elastic parameters by using amplitude and frequency information of pre-stack seismic data to guide flu...Frequency-dependent amplitude versus offset(FAVO)inversion is a popular method to estimate the frequency-dependent elastic parameters by using amplitude and frequency information of pre-stack seismic data to guide fluid identification.Current frequency-dependent AVO inversion methods are mainly based on elastic theory without the consideration of the viscoelasticity of oil/gas.A fluid discrimination approach is proposed in this study by incorporating the viscoelasticity and relevant FAVO inversion.Based on viscoelastic and rock physics theories,a frequency-dependent viscoelastic solid-liquid decoupling fluid factor is initially constructed,and its sensitivity in fluid discrimination is compared with other conventional fluid factors.Furthermore,a novel reflectivity equation is derived in terms of the viscoelastic solid-liquid decoupling fluid factor.Due to the introduction of viscoelastic theory,the proposed reflectivity is related to frequency,which is more widely applicable than the traditional elastic reflectivity equation directly derived from the elastic reflectivity equation on frequency.Finally,a pragmatic frequency-dependent inversion method is introduced to verify the feasibility of the equation for frequency-dependent viscoelastic solid-liquid decoupling fluid factor prediction.Synthetic and field data examples demonstrate the feasibility and stability of the proposed approach in fluid discrimination.展开更多
The nearly analytic discrete(NAD)method is a kind of finite difference method with advantages of high accuracy and stability.Previous studies have investigated the NAD method for simulating wave propagation in the tim...The nearly analytic discrete(NAD)method is a kind of finite difference method with advantages of high accuracy and stability.Previous studies have investigated the NAD method for simulating wave propagation in the time-domain.This study applies the NAD method to solving three-dimensional(3D)acoustic wave equations in the frequency-domain.This forward modeling approach is then used as the“engine”for implementing 3D frequency-domain full waveform inversion(FWI).In the numerical modeling experiments,synthetic examples are first given to show the superiority of the NAD method in forward modeling compared with traditional finite difference methods.Synthetic 3D frequency-domain FWI experiments are then carried out to examine the effectiveness of the proposed methods.The inversion results show that the NAD method is more suitable than traditional methods,in terms of computational cost and stability,for 3D frequency-domain FWI,and represents an effective approach for inversion of subsurface model structures.展开更多
To tackle the issue of notch frequency and center frequency drift of the L(0,1)mode guided wave in ultra⁃sonic guided wave⁃based stress monitoring of prestressed steel strands,a method using higher⁃order mode plateau ...To tackle the issue of notch frequency and center frequency drift of the L(0,1)mode guided wave in ultra⁃sonic guided wave⁃based stress monitoring of prestressed steel strands,a method using higher⁃order mode plateau fre⁃quencies is adopted.First,the correlation between group velocity peaks and phase velocities at these plateau frequen⁃cies is analyzed.This analysis establishes a quantitative rela⁃tionship between phase velocity and stress in the steel strand,providing a theoretical foundation for stress monitor⁃ing.Then the two⁃dimensional Fourier transform is em⁃ployed to separate wave modes.Dynamic programming techniques are applied in the frequency⁃velocity domain to extract higher⁃order modes.By identifying the group veloc⁃ity peaks of these separated higher⁃order modes,the plateau frequencies of guided waves are determined,enabling indi⁃rect measurement of stress in the steel strand.To validate this method,finite element simulations are conducted under three scenarios.Results show that the higher⁃order modes of transient signals from three different positions can be ac⁃curately extracted,leading to successful cable stress moni⁃toring.This approach effectively circumvents the issue of guided wave frequency drift and improves stress monitoring accuracy.Consequently,it significantly improves the appli⁃cation of ultrasonic guided wave technology in structural health monitoring.展开更多
基金Supported by National Nature Fund and National Civil-Defense Office
文摘Severe vibration of underground structures may be induced under blast loads. According to the characteristics of the explosion-induced ground shock wave, a new-type damper, inverse control magneto-rheological（MR） damper was designed to control the vibration, The high-frequency performance test of the MR damper was carried out on the small shaking table. It is shown that the performance can be modeled by use of the modified Bouc-Wen model, and the Parameters of the model keep stable in the range of 15--50 Hz.
文摘This study is an exploratory analysis of applying natural language processing techniques such as Term Frequency-Inverse Document Frequency and Sentiment Analysis on Twitter data. The uniqueness of this work is established by determining the overall sentiment of a politician’s tweets based on TF-IDF values of terms used in their published tweets. By calculating the TF-IDF value of terms from the corpus, this work displays the correlation between TF-IDF score and polarity. The results of this work show that calculating the TF-IDF score of the corpus allows for a more accurate representation of the overall polarity since terms are given a weight based on their uniqueness and relevance rather than just the frequency at which they appear in the corpus.
基金supported by the National Natural Science Foundation of China(Grant Nos.62271243 and 62071215)the National Key Research and Development Program of China(Grant No.2017YFA0700201)+1 种基金the Joint Fund of Ministry of Education for Equipment Pre-research(Grant No.8091B032112)the Priority Academic Program Development of Jiangsu Higher Education Institutions,the Fundamental Research Funds for the Central Universities,and Jiangsu Provincial Key Laboratory of Advanced Manipulating Technique of Electromagnetic Wave.
文摘Recently,the metasurfaces for independently controlling the wavefront and amplitude of two orthogonal circularly polarized electromagnetic(EM)waves have been demonstrated to open a way toward spin-multiplexing compact metadevices.However,these metasurfaces are mostly restricted to a single operation frequency band.The main challenge to achieving multiple frequency manipulations stems from the complicated and time-consuming design caused by multifrequency cross talk.To solve this problem,we propose a deep-learning-assisted inverse design method for designing a dual-spin/frequency metasurface with flexible multiplexing of off-axis vortices.By analyzing the cross talk between different spin/frequency channels based on the deep-learning method,we established the internal mapping relationship between the physical parameters of a meta-atom and its phase responses in multichannels,realizing the rapid inverse design of the spin/frequency multiplexing EM device.As a proof of concept,we demonstrated in the microwave region a dual-frequency arbitrary spin-to-orbit angular momentum converter,a dual-frequency off-axis vector vortex multiplexer,and a large-capacity(16-channel)vortex beam generator.The proposed method may provide a compact and efficient platform for the multiplexing of vortices,which may further stimulate their applications in wireless communication and quantum information science.
基金supported by the Key Project Fund of the Chinese Academy of Sciences under grant number (kzcx2-yw-203-01)the Major State Basic Research Development Program of China(973 Program,Grant No.2007CB41170404)
文摘Three-dimensional gravity inversion based on the mass property model is very popular in recent years. The time and efficiency of inversion algorithms is relative to the magnitude of the target mesh. One approach is to search over the entire solution space for a more refined result. However, the inversion will be difficult with the increased parameters in the large search space and the number of computations increases exponentially. ｜n this paper, we propose a novel approach based on the frequency characteristics of the density distribution over the mesh. The purposes of our study are to reduce the parameters of three- dimensional gravity inversion and to lighten the image quality of the inversion result. The results show that the new method can expedite the inversion processing and get a better geological interpretation than tradition methods.
基金supported by the China State Key Science and Technology Project on Marine Carbonate Reservoir Characterization (No. 2011ZX05004-003)the Basic Research Programs of CNPC during the 12th Five-Year Plan Period (NO.2011A-3603)+1 种基金the Natural Science Foundation of China (No.41104066)the RIPED Young Professional Innovation Fund (NO.2010-13-16-02, 2010-A-26-02)
文摘Although full waveform inversion in the frequency domain can overcome the local minima problem in the time direction, such problem still exists in the space direction because of the media subsurface complexity. Based on the optimal steep descent methods, we present an algorithm which combines the preconditioned bi-conjugated gradient stable method and the multi-grid method to compute the wave propagation and the gradient space. The multiple scale prosperity of the waveform inversion and the multi-grid method can overcome the inverse problems local minima defect and accelerate convergence. The local inhomogeneous three-hole model simulated results and the Marmousi model certify the algorithm effectiveness.
基金financially supported by the National High Technology Research and Development Program of China(No.2012AA09A20105)the National Science Foundation Network(No.41574127)
文摘We studied finite-element-method-based two-dimensional frequency-domain acoustic FWI under rugged topography conditions. The exponential attenuation boundary condition suitable for rugged topography is proposed to solve the cutoff botmdary problem as well as to consider the requirement of using the same subdivision grid in joint multifrequency inversion. The proposed method introduces the attenuation factor, and by adjusting it, acoustic waves are sufficiently attenuated in the attenuation layer to minimize the cutoff boundary effect. Based on the law of exponential attenuation, expressions for computing the attenuation factor and the thickness of attenuation layers are derived for different frequencies. In multifrequency-domain FWI, the conjugate gradient method is used to solve equations in the Gauss-Newton algorithm and thus minimize the computation cost in calculating the Hessian matrix. In addition, the effect of initial model selection and frequency combination on FWI is analyzed. Examples using numerical simulations and FWI calculations are used to verify the efficiency of the proposed method.
基金supported by the Natural Science Foundation of China(No.41374122)
文摘The complexity of an elastic wavefield increases the nonlinearity of inversion, To some extent, multiscale inversion decreases the nonlinearity of inversion and prevents it from falling into local extremes. A multiscale strategy based on the simultaneous use of frequency groups and layer stripping method based on damped wave field improves the stability of inversion. A dual-level parallel algorithm is then used to decrease the computational cost and improve practicability. The seismic wave modeling of a single frequency and inversion in a frequency group are computed in parallel by multiple nodes based on multifrontal massively parallel sparse direct solver and MPI. Numerical tests using an overthrust model show that the proposed inversion algorithm can effectively improve the stability and accuracy of inversion by selecting the appropriate inversion frequency and damping factor in low- frequency seismic data.
基金financially supported by the National Natural Science Foundation of China(No.41074075/D0409)the National Science and Technology Major Project(No.2011ZX05025-001-04)
文摘As a high quality seismic imaging method, full waveform inversion （FWI） can accurately reconstruct the physical parameter model for the subsurface medium. However, application of the FWI in seismic data processing is computationally expensive, especially for the three-dimension complex medium inversion. Introducing blended source technology into the frequency-domain FWI can greatly reduce the computational burden and improve the efficiency of the inversion. However, this method has two issues： first, crosstalk noise is caused by interference between the sources involved in the encoding, resulting in an inversion result with some artifacts; second, it is more sensitive to ambient noise compared to conventional FWI, therefore noisy data results in a poor inversion. This paper introduces a frequency-group encoding method to suppress crosstalk noise, and presents a frequency- domain auto-adapting FWI based on source-encoding technology. The conventional FWI method and source-encoding based FWI method are combined using an auto-adapting mechanism. This improvement can both guarantee the quality of the inversion result and maximize the inversion efficiency.
基金supported by the National Nature Science Foundation Project(Nos.41604101 and U1562215)the National Grand Project for Science and Technology(No.2016ZX05024-004)+2 种基金the Natural Science Foundation of Shandong(No.BS2014NJ005)Science Foundation from SINOPEC Key Laboratory of Geophysics(No.33550006-15-FW2099-0027)the Fundamental Research Funds for the Central Universities
文摘Amplitude variations with offset or incident angle （AVO/AVA） inversion are typically combined with statistical methods, such as Bayesian inference or deterministic inversion. We propose a joint elastic inversion method in the time and frequency domain based on Bayesian inversion theory to improve the resolution of the estimated P- and S-wave velocities and density. We initially construct the objective function using Bayesian inference by combining seismic data in the time and frequency domain. We use Cauchy and Gaussian probability distribution density functions to obtain the prior information for the model parameters and the likelihood function, respectively. We estimate the elastic parameters by solving the initial objective function with added model constraints to improve the inversion robustness. The results of the synthetic data suggest that the frequency spectra of the estimated parameters are wider than those obtained with conventional elastic inversion in the time domain. In addition, the proposed inversion approach offers stronger antinoising compared to the inversion approach in the frequency domain. Furthermore, results from synthetic examples with added Gaussian noise demonstrate the robustness of the proposed approach. From the real data, we infer that more model parameter details can be reproduced with the proposed joint elastic inversion.
文摘Short period surface waves generated by a local earthquake recorded by broadband seismometers at distances of about 186 to 778 km from the earthquake’s epicenter located in Cameroon (Central Africa) were processed for group velocity maps and dispersion waveforms using the frequency time analysis (FTAN) method. The resulting group velocity fundamental modes of the extracted Rayleigh and Love waves were used for a joint amplitude spectral and P polarity inversion using moment tensor inversion. The corresponding group velocity dispersion curves, the residual as a function of depth, the amplitude spectra and the moment tensor solutions of the regions from the epicenter to the different stations up to a depth of about 10 km were obtained.
基金supported by the National Natural Science Foundation of China(41704146)the Fundamental Research Funds for National Universities,China University of Geosciences(Wuhan)(CUGL180816)。
文摘Tunnel seismic detection methods are effective for obtaining the geological structure around the tunnel face,which is critical for safe construction and disaster mitigation in tunnel engineering.However,there is often a lack of accuracy in the acquired geological information and physical properties ahead of the tunnel face in the current tunnel seismic detection methods.Thus,we apply a frequency-domain acoustic full-waveform inversion(FWI)method to obtain high-resolution results for the tunnel structure.We discuss the influence of the frequency group selection strategy and the tunnel observation system settings regarding the inversion results and determine the structural imaging and physical property parameter inversion of abnormal geological bodies ahead of the tunnel face.Based on the conventional strategies of frequency-domain acoustic FWI,we propose a frequency group selection strategy that combines a low-frequency selection covering the vertical wavenumber and a high-frequency selection of antialiasing.This strategy can effectively obtain the spatial structure and physical parameters of the geology ahead of the tunnel face and improve the inversion resolution.In addition,by linearly increasing the side length of the tunnel observation system,we share the influence of the length of the two sides of the observation systems of different tunnels on the inversion results.We found out that the inversion results are the best when the side length is approximately five times the width of the tunnel face,and the influence of increasing the side observation length beyond this range on the inversion results can be ignored.Finally,based on this approach,we invert for the complex multi-stratum model,and an accurate structure and physical property parameters of the complex stratum ahead of the tunnel face are obtained,which verifies the feasibility of the proposed method.
基金the Strategic Priority Research Program of the Chinese Academy of Sciences(No.XDA14020102)the National Natural Science Foundation of China(Nos.41774125,41530320,41904104)+1 种基金the Key National Research Project of China(No.2018YFC0603300)the S&T Program of Beijing(No.Z181100005718001).
文摘We propose a new 3D inversion scheme to invert the near-and transition-zone data of CSAMT with topography accurately.In this new method,the earth was discretized into unstructured tetrahedra to fit the ragged topography and the vector fi nite-element method was adopted to obtain precise responses and good sensitivity.To simulate the attitude and shape of the transmitter,we divided a long-grounded transmitter into dipoles and integrated these dipoles to obtain good responses in the near-and transition-fi eld zones.Next,we designed an L2 norm-based objective functional and applied a standard quasi-Newton method as the optimization method to solve the inverse problem and guarantee steady convergence.We tested our 3D inversion method first on synthetic data and then on a field dataset acquired from select sites near Changbai Mountain,China.In both tests,the new inversion algorithm achieved excellent fitting between the predicted and observed data,even in near-and transition-fi eld zones,and the inversion results agreed well with the true model.These fi ndings reveal that the proposed algorithm is eff ective for 3D inversion of CSAMT data.
基金supported by the National Natural Science Foundation of China(No.41374078)geological survey project of Ministry of Land and Resources of China(No.12120113101300)
文摘The controlled source extremely low frequency （CSELF） electromagnetic method is characterized by extremely long and powerful sources and a huge measurement range. Its electromagnetic field can therefore be affected by the ionosphere and displacement current. Research on 3D forward modeling and inversion of CSELF electromagnetic data is currently in its infancy. This paper makes exploratory attempts to firstly calculate the 1D extremely low frequency electromagnetic field under ionosphere-air-earth coupling circumstances, and secondly analyze the propagation characteristics of the background electromagnetic field. The 3D staggered-grid finite difference scheme for solving for the secondary electric field is adopted and incorporated with the 1D modeling algorithm to complete 3D forward modeling. Considering that surveys can be carried out in the near field and transition zone for lower frequencies, the 3D Limited-memory Broyden-Fletcher-Goldfarb-Shanno （LBFGS） inversion of CSELF electromagnetic data is presented （in which the sources, or primary fields, are included）, with the aim of directly inverting the impedance data, regardless of where it is acquired. Derivation of the objective functional gradient is the core component in the inversion. Synthetic tests indicate that the well-chosen approximation to the Hessian can significantly speed up the inversion. The model responses corresponding to the coexistence of conductive and resistive blocks show that the off-diagonal components of tensor impedance are much more sensitive to the resistivity variation than the diagonal components. In comparison with conventional scalar inversion, tensor inversion is superior in the recoveries of electric anomalies and background resistivity.
基金funded by theNationalNatural Science Foundation of China(52061020)Major Science and Technology Projects in Yunnan Province(202302AG050009)Yunnan Fundamental Research Projects(202301AV070003).
文摘Finding materials with specific properties is a hot topic in materials science.Traditional materials design relies on empirical and trial-and-error methods,requiring extensive experiments and time,resulting in high costs.With the development of physics,statistics,computer science,and other fields,machine learning offers opportunities for systematically discovering new materials.Especially through machine learning-based inverse design,machine learning algorithms analyze the mapping relationships between materials and their properties to find materials with desired properties.This paper first outlines the basic concepts of materials inverse design and the challenges faced by machine learning-based approaches to materials inverse design.Then,three main inverse design methods—exploration-based,model-based,and optimization-based—are analyzed in the context of different application scenarios.Finally,the applications of inverse design methods in alloys,optical materials,and acoustic materials are elaborated on,and the prospects for materials inverse design are discussed.The authors hope to accelerate the discovery of new materials and provide new possibilities for advancing materials science and innovative design methods.
基金supported by the China National Petroleum Corporation Scientific research and technology development project(Nos.2013E-38-08)
文摘Multi-wave exploration is an effective means for improving precision in the exploration and development of complex oil and gas reservoirs that are dense and have low permeability. However, convened wave data is characterized by a low signal-to-noise ratio and low resolution, because the conventional deconvolution technology is easily affected by the frequency range limits, and there is limited scope for improving its resolution. The spectral inversion techniques is used to identify λ/8 thin layers and its breakthrough regarding band range limits has greatly improved the seismic resolution. The difficulty associated with this technology is how to use the stable inversion algorithm to obtain a high-precision reflection coefficient, and then to use this reflection coefficient to reconstruct broadband data for processing. In this paper, we focus on how to improve the vertical resolution of the converted PS-wave for multi-wave data processing. Based on previous research, we propose a least squares inversion algorithm with a total variation constraint, in which we uses the total variance as a priori information to solve under-determined problems, thereby improving the accuracy and stability of the inversion. Here, we simulate the Gaussian fitting amplitude spectrum to obtain broadband wavelet data, which we then process to obtain a higher resolution converted wave. We successfully apply the proposed inversion technology in the processing of high-resolution data from the Penglai region to obtain higher resolution convened wave data, which we then verify in a theoretical test. Improving the resolution of converted PS-wave data will provide more accurate data for subsequent velocity inversion and the extraction of reservoir reflection information.
基金supported by the National Major Science and Technology Project of China(2016ZX05004003)the National Natural Science Foundation of China(41574103,41974120,U20B2015)Open Fund of State Key Laboratory of Coal Resources and Safe Mining(Grant No.SKLCRSM19KFA08)。
文摘Fluid and effective fracture identification in reservoirs is a crucial part of reservoir prediction.The frequency-dependent AVO inversion algorithms have proven to be effective for identifying fluid through its dispersion property.However,the conventional frequency-dependent AVO inversion algorithms based on Smith&Gidlow and Aki&Richards approximations do not consider the acquisition azimuth of seismic data and neglect the effect of seismic anisotropic dispersion in the actual medium.The aligned fractures in the subsurface medium induce anisotropy.The seismic anisotropy should be considered while accounting for the seismic dispersion properties through fluid-saturated fractured reservoirs.Anisotropy in such reservoirs is frequency-related due to wave-induced fluid-flow(WIFF)between interconnected fractures and pores.It can be used to identify fluid and effective fractures(fluid-saturated)by using azimuthal seismic data via anisotropic dispersion properties.In this paper,based on Rüger’s equation,we derived an analytical expression in the frequency domain for the frequencydependent AVOAz inversion in terms of fracture orientation,dispersion gradient of isotropic background rock,anisotropic dispersion gradient,and the dispersion at a normal incident angle.The frequency-dependent AVOAz equation utilizes azimuthal seismic data and considers the effect of both isotropic and anisotropic dispersion.Reassigned Gabor Transform(RGT)is used to achieve highresolution frequency division data.We then propose the frequency-dependent AVOAz inversion method to identify fluid and characterize effective fractures in fractured porous reservoirs.Through application to high-qualified seismic data of dolomite and carbonate reservoirs,the results show that the method is useful for identifying fluid and effective fractures in fluid-saturated fractured rocks.
基金the sponsorship of National Natural Science Foundation of China(41974119,U1762103)Science Foundation from Innovation and Technology Support Program for Young Scientists in Colleges of Shandong province and Ministry of Science and Technology of China(2020RA2C620131)。
文摘Frequency-dependent amplitude versus offset(FAVO)inversion is a popular method to estimate the frequency-dependent elastic parameters by using amplitude and frequency information of pre-stack seismic data to guide fluid identification.Current frequency-dependent AVO inversion methods are mainly based on elastic theory without the consideration of the viscoelasticity of oil/gas.A fluid discrimination approach is proposed in this study by incorporating the viscoelasticity and relevant FAVO inversion.Based on viscoelastic and rock physics theories,a frequency-dependent viscoelastic solid-liquid decoupling fluid factor is initially constructed,and its sensitivity in fluid discrimination is compared with other conventional fluid factors.Furthermore,a novel reflectivity equation is derived in terms of the viscoelastic solid-liquid decoupling fluid factor.Due to the introduction of viscoelastic theory,the proposed reflectivity is related to frequency,which is more widely applicable than the traditional elastic reflectivity equation directly derived from the elastic reflectivity equation on frequency.Finally,a pragmatic frequency-dependent inversion method is introduced to verify the feasibility of the equation for frequency-dependent viscoelastic solid-liquid decoupling fluid factor prediction.Synthetic and field data examples demonstrate the feasibility and stability of the proposed approach in fluid discrimination.
基金supported by the Joint Fund of Seismological Science(Grant No.U1839206)the National R&D Program on Monitoring,Early Warning and Prevention of Major Natural Disaster(Grant No.2017YFC1500301)+2 种基金supported by IGGCAS Research Start-up Funds(Grant No.E0515402)National Natural Science Foundation of China(Grant No.E1115401)supported by National Natural Science Foundation of China(Grant No.11971258).
文摘The nearly analytic discrete(NAD)method is a kind of finite difference method with advantages of high accuracy and stability.Previous studies have investigated the NAD method for simulating wave propagation in the time-domain.This study applies the NAD method to solving three-dimensional(3D)acoustic wave equations in the frequency-domain.This forward modeling approach is then used as the“engine”for implementing 3D frequency-domain full waveform inversion(FWI).In the numerical modeling experiments,synthetic examples are first given to show the superiority of the NAD method in forward modeling compared with traditional finite difference methods.Synthetic 3D frequency-domain FWI experiments are then carried out to examine the effectiveness of the proposed methods.The inversion results show that the NAD method is more suitable than traditional methods,in terms of computational cost and stability,for 3D frequency-domain FWI,and represents an effective approach for inversion of subsurface model structures.
基金The National Natural Science Foundation of China(No.52278303).
文摘To tackle the issue of notch frequency and center frequency drift of the L(0,1)mode guided wave in ultra⁃sonic guided wave⁃based stress monitoring of prestressed steel strands,a method using higher⁃order mode plateau fre⁃quencies is adopted.First,the correlation between group velocity peaks and phase velocities at these plateau frequen⁃cies is analyzed.This analysis establishes a quantitative rela⁃tionship between phase velocity and stress in the steel strand,providing a theoretical foundation for stress monitor⁃ing.Then the two⁃dimensional Fourier transform is em⁃ployed to separate wave modes.Dynamic programming techniques are applied in the frequency⁃velocity domain to extract higher⁃order modes.By identifying the group veloc⁃ity peaks of these separated higher⁃order modes,the plateau frequencies of guided waves are determined,enabling indi⁃rect measurement of stress in the steel strand.To validate this method,finite element simulations are conducted under three scenarios.Results show that the higher⁃order modes of transient signals from three different positions can be ac⁃curately extracted,leading to successful cable stress moni⁃toring.This approach effectively circumvents the issue of guided wave frequency drift and improves stress monitoring accuracy.Consequently,it significantly improves the appli⁃cation of ultrasonic guided wave technology in structural health monitoring.
