To fully extract and mine the multi-scale features of reservoirs and geologic structures in time/depth and space dimensions, a new 3D multi-scale volumetric curvature （MSVC） methodology is presented in this paper. W...To fully extract and mine the multi-scale features of reservoirs and geologic structures in time/depth and space dimensions, a new 3D multi-scale volumetric curvature （MSVC） methodology is presented in this paper. We also propose a fast algorithm for computing 3D volumetric curvature. In comparison to conventional volumetric curvature attributes, its main improvements and key algorithms introduce multi-frequency components expansion in time-frequency domain and the corresponding multi-scale adaptive differential operator in the wavenumber domain, into the volumetric curvature calculation. This methodology can simultaneously depict seismic multi-scale features in both time and space. Additionally, we use data fusion of volumetric curvatures at various scales to take full advantage of the geologic features and anomalies extracted by curvature measurements at different scales. The 3D MSVC can highlight geologic anomalies and reduce noise at the same time. Thus, it improves the interpretation efficiency of curvature attributes analysis. The 3D MSVC is applied to both land and marine 3D seismic data. The results demonstrate that it can indicate the spatial distribution of reservoirs, detect faults and fracture zones, and identify their multi-scale properties.展开更多
A new approach for unparallel trajectory bistatic spotlight SAR imaging is proposed. The approach utilizes the concept of instantaneous Doppler wavenumber and introduces two variants, the sum-range and subtraction-ran...A new approach for unparallel trajectory bistatic spotlight SAR imaging is proposed. The approach utilizes the concept of instantaneous Doppler wavenumber and introduces two variants, the sum-range and subtraction-range, to develop the 2D frequency analytical formula. Based on the assumption of plane wavefront, the transmitting and receiving Doppler are separated and formulated via series reversion. And frequency scaling is applied to focus image. The algorithm is with high computational efficiency, and provides well focus for limited scene imaging. Simulation result confirms the validity of the approach.展开更多
Monitoring of the South China Sea is always one of the focuses in the field of ocean engineering for its particular geographic position. The modal wave number tomography is proposed for monitoring the front and numeri...Monitoring of the South China Sea is always one of the focuses in the field of ocean engineering for its particular geographic position. The modal wave number tomography is proposed for monitoring the front and numerical simulation is performed for the front of the South China Sea. With the empirical orthonormal function (EOF) applied to reduce the parameter search space, the perturbation inversion method is used for inversing sound speed profiles. The 2-D ocean environment used for numerical simulation is selected from the ocean area, located in 20°N, 118°E～20°N, 125°E, near the Luzon Strait in the South China Sea. The ocean environment sound speed distribution in the ocean area under study is obtained from the assimilation of multi-source remote sensing data. The numerical simulation shows that the modal wave number tomography can inverse the average sound speed profile, therefore,it can be used to monitor ocean internal structures such as ocean fronts and eddies which affect sound speed distribution.展开更多
An improved circular synthetic aperture radar（CSAR） imaging algorithm of omega-k（ω-k） type mainly for reconstructing an image on a cylindrical surface is proposed.In the typical CSAR ω-k algorithm,the rage traje...An improved circular synthetic aperture radar（CSAR） imaging algorithm of omega-k（ω-k） type mainly for reconstructing an image on a cylindrical surface is proposed.In the typical CSAR ω-k algorithm,the rage trajectory is approximated by Taylor series expansion to the quadratic terms,which limits the valid synthetic aperture length and the angular reconstruction range severely.Based on the model of the CSAR echo signal,the proposed algorithm directly transforms the signal to the two-dimensional（2D） wavenumber domain,not using approximation processing to the range trajectory.Based on form of the signal spectrum in the wavenumber domain,the formula for the wavenumber domain interpolation of the w-k algorithm is deduced,and the wavenumber spectrum of the reference point used for bulk compression is obtained from numerical method.The improved CSAR ω-k imaging algorithm increases the valid synthetic aperture length and the angular area greatly and hence improves the angular resolution of the cylindrical imaging.Additionally,the proposed algorithm can be repeated on different cylindrical surfaces to achieve three dimensional（3D） image reconstruction.The 3D spatial resolution of the CSAR system is discussed,and the simulation results validate the correctness of the analysis and the feasibility of the algorithm.展开更多
A correlation equation between the UV absorption wavenumbers of 1,4-disubstituted benzenes and the excited-state substituent constant was obtained. For 80 sorts of 1,4- disubstituted benzenes, the correlation coeffici...A correlation equation between the UV absorption wavenumbers of 1,4-disubstituted benzenes and the excited-state substituent constant was obtained. For 80 sorts of 1,4- disubstituted benzenes, the correlation coefficient was 0.9805, and the standard deviation was only 672.27 cm^-1. The results imply that the excited-state substituent constant can be used productively for research on UV energy of 1,4-disubstituted benzenes. The present method provides a new avenue to study the UV absorption spectra of aromatic systems with the excited-state substituent constant, and it is helpful to understand the effect of substituent electrostatic effects on the chemical and physical properties of conjugated compounds with multiple substituents in excited state.展开更多
In this paper, the analytical representations of four wave source functions in high-frequency spectrum range are given on the basis of ocean wave theory and dimensional analysis, and the perturbation method is used to...In this paper, the analytical representations of four wave source functions in high-frequency spectrum range are given on the basis of ocean wave theory and dimensional analysis, and the perturbation method is used to solve the governing equations of ocean wave high-frequency spectrum on the basis of the temporally stationary and locally homogeneous scale relations of microscale wave. The microscale ocean wavenumber spectrum correct to the second order has an explicit structure, its first order part represents the equilibrium between dif- ferent source functions, and its second order part represents the contribution of microscale wave propagation.展开更多
This paper utilizes meridional and zonal velocity data measured by chaff rockets at a high resolution of 25 m in the 67—96 km region during the MAC/SINE campaign to study the char-acteristics of the spectra of the me...This paper utilizes meridional and zonal velocity data measured by chaff rockets at a high resolution of 25 m in the 67—96 km region during the MAC/SINE campaign to study the char-acteristics of the spectra of the meridional and zonal velocity fluctuations. Results suggest that the mean vertical wavenumber spectra have a forked point near vertical wavenumber k = 0.001 m-1. In the vertical wavenumber region smaller than the forked point, the spectra of the meridional and zonal velocities have approximately the same spectral slope, amplitude, and characteristic vertical wavelength, and agree well with the linear saturated spectral model. In the vertical wavenumber region larger than the forked point, however, the meridional and zonal spectra are inconsistent not only in both slope and amplitude, but also in spectral structure. The differences between the observed meridional and zonal spectra provide observational evidence of azimuthal anisotropy of the velocity fluctuation field in the large vertical wavenumber region in mean sense.展开更多
An exact solution based on the wavenumber integration method is proposed and implemented in a numerical model for the acoustic field in a Pekeris waveguide excited by either a point source in cylindrical geometry or a...An exact solution based on the wavenumber integration method is proposed and implemented in a numerical model for the acoustic field in a Pekeris waveguide excited by either a point source in cylindrical geometry or a line source in plane geometry. Besides, an unconditionally stable numerical solution is also presented, which entirely resolves the stability problem in previous methods. Generally the branch line integral contributes to the total field only at short ranges, and hence is usually ignored in traditional normal mode models. However, for the special case where a mode lies near the branch cut, the branch line integral can contribute to the total field significantly at all ranges. The wavenumber integration method is well-suited for such problems. Numerical results are also provided, which show that the present model can serve as a benchmark for sound propagation in a Pekeris waveguide.展开更多
A correlation analysis is performed to investigate the relationship between El Nino-Southern Oscillation （ENSO） and the Antarctic oscillation （AAO） at the quasi-quadrennial （QQ） timescale.It is found that the co...A correlation analysis is performed to investigate the relationship between El Nino-Southern Oscillation （ENSO） and the Antarctic oscillation （AAO） at the quasi-quadrennial （QQ） timescale.It is found that the cold tongue index （CTI） and the AAO index （AAOI） are negatively correlated with about a 7-month lead-time,while they are positively correlated with about a 15-month lag-time.To further explore this relationship,complex empirical orthogonal function analysis is employed in the QQ sea level pressure （SLP） anomalies from 1951 to 2002.The results indicate that,during the ENSO cycle,there exists one kind of global tropical wave of wavenumber 1 （GTW1） propagating eastward.With the traveling of GTW1,the tropical SLP anomaly tends to intrude into the southern mid-latitudes.Accordingly,three strong signals travel synchronously along the circumSouth-Pacific path,and a relatively weak signal extends eastward and poleward over the South Ocean in the Atlantic-Indian Ocean sector.Following the propagation of these signals,the AAO phase tends to be reversed progressively.As a result,there exists an evident lead-lag correlation between CTI and AAOI.It can be concluded that ENSO plays a key role in the phase transition of AAO at the QQ timescale.It is also noticed that this regular relationship is only evident in the canonical ENSO events,for which sea surface temperature （SST） anomalies extend westward from the tropical eastern Pacific.On the other hand,the similar relationships are not found among those atypical ENSO events for which SST anomalies spread eastward from the central Pacific,such as the 1982-1983 ENSO event.展开更多
It is difficult to quantitatively detect defects by using the time domain or frequency domain features of Lamb wave signals due to their dispersion and multimodal characteristics.Therefore,it is important to discover ...It is difficult to quantitatively detect defects by using the time domain or frequency domain features of Lamb wave signals due to their dispersion and multimodal characteristics.Therefore,it is important to discover an intrinsical parameter of Lamb waves that could be used as a damage sensitive feature.In this paper,quantitative defect detection in aluminium plates is carried out by means of wavenumber analysis approach.The wavenumber of excited Lamb wave mode is a fixed value,given a frequency,a thickness and material properties of the target plate.When Lamb waves propagate to the structural discontinuity,new wavenumber components are created by abrupt wavefield change.The new wavenumber components can be identified in the frequency-wavenumber domain.To estimate spatially dependent wavenumber values,a short-space two-dimensional Fourier transform(FT)method is presented for processing wavefield data of Lamb waves.The results can be used to determine the location,size and depth of rectangular notch.The analysis techniques are demonstrated using simulation examples of an aluminium plate with a rectangular notch.Then,the wavenumber analysis method is applied to simulation data that are obtained through a range of notch depths and widths.The results are analyzed and rules of the technique with regards to estimating notch depth are determined.Based on simulation results,guidelines for using the technique are developed.Finally,experimental wavefield data are obtained in aluminium plates with rectangular notches by a full noncontact transceiving method,i.e.,laser-laser method.Band-pass filtering combined with continuous wavelet transform is used to extract a certain frequency component from the full laser-induced wavefield with wide band.Shortspace two-dimensional FT method is used for further processing full wavefield data at a certain frequency to estimate spatially dependent wavenumber values.The consistency of simulation and experimental results shows the effectiveness of proposed wavenumber method for quantitative展开更多
The reliable estimation of the wavenumber space(k-space)of the plates remains a longterm concern for acoustic modeling and structural dynamic behavior characterization.Most current analyses of wavenumber identificatio...The reliable estimation of the wavenumber space(k-space)of the plates remains a longterm concern for acoustic modeling and structural dynamic behavior characterization.Most current analyses of wavenumber identification methods are based on the deterministic hypothesis.To this end,an inverse method is proposed for identifying wave propagation characteristics of twodimensional structures under stochastic conditions,such as wavenumber space,dispersion curves,and band gaps.The proposed method is developed based on an algebraic identification scheme in the polar coordinate system framework,thus named Algebraic K-Space Identification(AKSI)technique.Additionally,a model order estimation strategy and a wavenumber filter are proposed to ensure that AKSI is successfully applied.The main benefit of AKSI is that it is a reliable and fast method under four stochastic conditions:(A)High level of signal noise;(B)Small perturbation caused by uncertainties in measurement points’coordinates;(C)Non-periodic sampling;(D)Unknown structural periodicity.To validate the proposed method,we numerically benchmark AKSI and three other inverse methods to extract dispersion curves on three plates under stochastic conditions.One experiment is then performed on an isotropic steel plate.These investigations demonstrate that AKSI is a good in-situ k-space estimator under stochastic conditions.展开更多
Magnetic field gradient tensor technique provides abundant data for delicate inversion of subsurface magnetic susceptibility distribution. Large scale magnetic data inversion imaging requires high speed and accuracy f...Magnetic field gradient tensor technique provides abundant data for delicate inversion of subsurface magnetic susceptibility distribution. Large scale magnetic data inversion imaging requires high speed and accuracy for forward modeling. For arbitrarily distributed susceptibility data on an undulated surface, we propose a fast 3D forward modeling method in the wavenumber domain based on(1) the wavenumber-domain expression of the prism combination model and the Gauss–FFT algorithm and(2) cubic spline interpolation. We apply the proposed 3D forward modeling method to synthetic data and use weighting coefficients in the wavenumber domain to improve the modeling for multiple observation surfaces, and also demonstrate the accuracy and efficiency of the proposed method.展开更多
It is well known that quantitative estimation of slip distributions on fault plane is one of the most important issues for earthquake source inversion related to the fault rupture process. The characteristics of slip ...It is well known that quantitative estimation of slip distributions on fault plane is one of the most important issues for earthquake source inversion related to the fault rupture process. The characteristics of slip distribution on the main fault play a fundamental role to control strong ground motion pattern. A large amount of works have also suggested that variable slip models inverted from longer period ground motion recordings are relevant for the prediction of higher frequency ground motions. Zhang et al. （Chin J Geophys 56：1412-1417, 2013） and Wang et al. （Chin J Geophys 56：1408-1411,2013） published their source inversions for the fault rupturing process soon after the April 20, 2013 Lushan earthquake in Sichuan, China. In this study, first, we synthesize two forward source slip models： the value of maximum slip, fault dimension, size, and dimension of major asperities, and comer wave number obtained from Wang＇s model is adopted to constrain the gen- eration of k-2 model and crack model. Next, both inverted and synthetic slip models are employed to simulate the ground motions for the Lushan earthquake based on the stochastic finite-fault method. In addition, for a comparison purpose, a stochastic slip model and another k-2 model （k 2 model II） with 2 times value of comer wave number of the original k-2 model （k 2 model I） are also employed for simulation for Lushan event. The simulated results characterized by Modified Mer- calli Intensity （MMI） show that the source slip models based on the inverted and synthetic slip distributions could capture many basic features associated with the ground motion patterns. Moreover, the simulated MMI distributions reflect the rupture directivity effect and the influence of the shallow velocity structure well. On the other hand, the simulated MMI bystochastic slip model and k 2 model II is apparently higher than observed intensity. By contrast, our simulation results show that the higher frequency ground motion is sensitive to the degree of sli展开更多
The current local wavenumber methods for the interpretation of magnetic anomalies compute the locations of geological bodies by solving complex matrices. Presently, such methods require to know the structural index, w...The current local wavenumber methods for the interpretation of magnetic anomalies compute the locations of geological bodies by solving complex matrices. Presently, such methods require to know the structural index, which is a parameter that represents the source type. The structural index is hard to know in real data; consequently, the precision of current methods is low. We present the fast local wavenumber （FLW） method, and define the squared sum of the horizontal and vertical local wavenumbers as the cumulative local wavenumber. The FLW method is the linear combination of the umulative local wavenumberand other wavenumbers, and is used to compute the locations and structural index of the source without a priori information and matrix solution. We apply the FLW method to synthetic magnetic anomalies, and the results suggest that the FLW method is insensitive to background and oblique magnetization. Next, we apply the FLW method to real magnetic data to obtain the location and structural index of the source.展开更多
A new finite difference scheme-SCSD scheme has been proposed based on CD (Central Difference)scheme and SUD (Secondr-order Upwind Difference) scheme. Its basic feature is controIlable convectivesfability and always se...A new finite difference scheme-SCSD scheme has been proposed based on CD (Central Difference)scheme and SUD (Secondr-order Upwind Difference) scheme. Its basic feature is controIlable convectivesfability and always second-order accuracy (Stability-Controllable Second-order Difference ). It hasbeen proven that this scheme is convective-stable if the grid Peclet number .The advanage of this new scheme has been discussed based on the modified wavenumber analysis byusing Fourier transform. This scheme has been applied to the 2-D incompressible convective-diffusiveequation and 2-D compressible Euler equation, and corresponding finite difference equations have beenderived. Numerical examples of 1-D Burgers equation and 2-D transport equation have been presentedto show its good accuracy and controllable convective stability展开更多
基金supported by the National Natural Science Foundation of China (No. 41004054) Research Fund for the Doctoral Program of Higher Education of China (No. 20105122120002)Natural Science Key Project, Sichuan Provincial Department of Education (No. 092A011)
文摘To fully extract and mine the multi-scale features of reservoirs and geologic structures in time/depth and space dimensions, a new 3D multi-scale volumetric curvature （MSVC） methodology is presented in this paper. We also propose a fast algorithm for computing 3D volumetric curvature. In comparison to conventional volumetric curvature attributes, its main improvements and key algorithms introduce multi-frequency components expansion in time-frequency domain and the corresponding multi-scale adaptive differential operator in the wavenumber domain, into the volumetric curvature calculation. This methodology can simultaneously depict seismic multi-scale features in both time and space. Additionally, we use data fusion of volumetric curvatures at various scales to take full advantage of the geologic features and anomalies extracted by curvature measurements at different scales. The 3D MSVC can highlight geologic anomalies and reduce noise at the same time. Thus, it improves the interpretation efficiency of curvature attributes analysis. The 3D MSVC is applied to both land and marine 3D seismic data. The results demonstrate that it can indicate the spatial distribution of reservoirs, detect faults and fracture zones, and identify their multi-scale properties.
基金Supported by the National Natural Science Foundation of China (Grant No. 60502044)Program for New Century Excellent Talents inUniversity (Grant No. NCET-06-0861)
文摘A new approach for unparallel trajectory bistatic spotlight SAR imaging is proposed. The approach utilizes the concept of instantaneous Doppler wavenumber and introduces two variants, the sum-range and subtraction-range, to develop the 2D frequency analytical formula. Based on the assumption of plane wavefront, the transmitting and receiving Doppler are separated and formulated via series reversion. And frequency scaling is applied to focus image. The algorithm is with high computational efficiency, and provides well focus for limited scene imaging. Simulation result confirms the validity of the approach.
文摘Monitoring of the South China Sea is always one of the focuses in the field of ocean engineering for its particular geographic position. The modal wave number tomography is proposed for monitoring the front and numerical simulation is performed for the front of the South China Sea. With the empirical orthonormal function (EOF) applied to reduce the parameter search space, the perturbation inversion method is used for inversing sound speed profiles. The 2-D ocean environment used for numerical simulation is selected from the ocean area, located in 20°N, 118°E～20°N, 125°E, near the Luzon Strait in the South China Sea. The ocean environment sound speed distribution in the ocean area under study is obtained from the assimilation of multi-source remote sensing data. The numerical simulation shows that the modal wave number tomography can inverse the average sound speed profile, therefore,it can be used to monitor ocean internal structures such as ocean fronts and eddies which affect sound speed distribution.
文摘An improved circular synthetic aperture radar（CSAR） imaging algorithm of omega-k（ω-k） type mainly for reconstructing an image on a cylindrical surface is proposed.In the typical CSAR ω-k algorithm,the rage trajectory is approximated by Taylor series expansion to the quadratic terms,which limits the valid synthetic aperture length and the angular reconstruction range severely.Based on the model of the CSAR echo signal,the proposed algorithm directly transforms the signal to the two-dimensional（2D） wavenumber domain,not using approximation processing to the range trajectory.Based on form of the signal spectrum in the wavenumber domain,the formula for the wavenumber domain interpolation of the w-k algorithm is deduced,and the wavenumber spectrum of the reference point used for bulk compression is obtained from numerical method.The improved CSAR ω-k imaging algorithm increases the valid synthetic aperture length and the angular area greatly and hence improves the angular resolution of the cylindrical imaging.Additionally,the proposed algorithm can be repeated on different cylindrical surfaces to achieve three dimensional（3D） image reconstruction.The 3D spatial resolution of the CSAR system is discussed,and the simulation results validate the correctness of the analysis and the feasibility of the algorithm.
基金ACKNOWLEDGMENTS This work was supported by the National Natural Science Foundation of China （No.20772028 and No.20472019）, the Natural Science Foundation of Hunan Province （No.06JJ2002）, and the Scientific Research Fund of Hunan Provincial Education Commission.
文摘A correlation equation between the UV absorption wavenumbers of 1,4-disubstituted benzenes and the excited-state substituent constant was obtained. For 80 sorts of 1,4- disubstituted benzenes, the correlation coefficient was 0.9805, and the standard deviation was only 672.27 cm^-1. The results imply that the excited-state substituent constant can be used productively for research on UV energy of 1,4-disubstituted benzenes. The present method provides a new avenue to study the UV absorption spectra of aromatic systems with the excited-state substituent constant, and it is helpful to understand the effect of substituent electrostatic effects on the chemical and physical properties of conjugated compounds with multiple substituents in excited state.
文摘In this paper, the analytical representations of four wave source functions in high-frequency spectrum range are given on the basis of ocean wave theory and dimensional analysis, and the perturbation method is used to solve the governing equations of ocean wave high-frequency spectrum on the basis of the temporally stationary and locally homogeneous scale relations of microscale wave. The microscale ocean wavenumber spectrum correct to the second order has an explicit structure, its first order part represents the equilibrium between dif- ferent source functions, and its second order part represents the contribution of microscale wave propagation.
基金Acknowledgements This work was supported by the National Natural Science Foundation of China (Grant Nos. 49974038 and 49990450) and by the National Basic Research Project of China (Grant No. G2000078407).
文摘This paper utilizes meridional and zonal velocity data measured by chaff rockets at a high resolution of 25 m in the 67—96 km region during the MAC/SINE campaign to study the char-acteristics of the spectra of the meridional and zonal velocity fluctuations. Results suggest that the mean vertical wavenumber spectra have a forked point near vertical wavenumber k = 0.001 m-1. In the vertical wavenumber region smaller than the forked point, the spectra of the meridional and zonal velocities have approximately the same spectral slope, amplitude, and characteristic vertical wavelength, and agree well with the linear saturated spectral model. In the vertical wavenumber region larger than the forked point, however, the meridional and zonal spectra are inconsistent not only in both slope and amplitude, but also in spectral structure. The differences between the observed meridional and zonal spectra provide observational evidence of azimuthal anisotropy of the velocity fluctuation field in the large vertical wavenumber region in mean sense.
基金supported by the National Natural Science Foundation of China(Grant No.11125420)the Knowledge Innovation Program of the Chinese Academy of Sciences+2 种基金the China Postdoctoral Science Foundation(Grant No.2014M561882)the Doctoral Fund of Shandong ProvinceChina(Grant No.BS2012HZ015)
文摘An exact solution based on the wavenumber integration method is proposed and implemented in a numerical model for the acoustic field in a Pekeris waveguide excited by either a point source in cylindrical geometry or a line source in plane geometry. Besides, an unconditionally stable numerical solution is also presented, which entirely resolves the stability problem in previous methods. Generally the branch line integral contributes to the total field only at short ranges, and hence is usually ignored in traditional normal mode models. However, for the special case where a mode lies near the branch cut, the branch line integral can contribute to the total field significantly at all ranges. The wavenumber integration method is well-suited for such problems. Numerical results are also provided, which show that the present model can serve as a benchmark for sound propagation in a Pekeris waveguide.
基金The Major National Scientific Research Project on Global Change under contract No.2010CB951901the National Science Foundation of China under contract No.40821092Special Fund for Public Welfare Industry (Meteorology) under contract No.GYHY200906018
文摘A correlation analysis is performed to investigate the relationship between El Nino-Southern Oscillation （ENSO） and the Antarctic oscillation （AAO） at the quasi-quadrennial （QQ） timescale.It is found that the cold tongue index （CTI） and the AAO index （AAOI） are negatively correlated with about a 7-month lead-time,while they are positively correlated with about a 15-month lag-time.To further explore this relationship,complex empirical orthogonal function analysis is employed in the QQ sea level pressure （SLP） anomalies from 1951 to 2002.The results indicate that,during the ENSO cycle,there exists one kind of global tropical wave of wavenumber 1 （GTW1） propagating eastward.With the traveling of GTW1,the tropical SLP anomaly tends to intrude into the southern mid-latitudes.Accordingly,three strong signals travel synchronously along the circumSouth-Pacific path,and a relatively weak signal extends eastward and poleward over the South Ocean in the Atlantic-Indian Ocean sector.Following the propagation of these signals,the AAO phase tends to be reversed progressively.As a result,there exists an evident lead-lag correlation between CTI and AAOI.It can be concluded that ENSO plays a key role in the phase transition of AAO at the QQ timescale.It is also noticed that this regular relationship is only evident in the canonical ENSO events,for which sea surface temperature （SST） anomalies extend westward from the tropical eastern Pacific.On the other hand,the similar relationships are not found among those atypical ENSO events for which SST anomalies spread eastward from the central Pacific,such as the 1982-1983 ENSO event.
基金supported by the National Natural Science Foundation of China(Nos.51475012,11772014,and 11272021)
文摘It is difficult to quantitatively detect defects by using the time domain or frequency domain features of Lamb wave signals due to their dispersion and multimodal characteristics.Therefore,it is important to discover an intrinsical parameter of Lamb waves that could be used as a damage sensitive feature.In this paper,quantitative defect detection in aluminium plates is carried out by means of wavenumber analysis approach.The wavenumber of excited Lamb wave mode is a fixed value,given a frequency,a thickness and material properties of the target plate.When Lamb waves propagate to the structural discontinuity,new wavenumber components are created by abrupt wavefield change.The new wavenumber components can be identified in the frequency-wavenumber domain.To estimate spatially dependent wavenumber values,a short-space two-dimensional Fourier transform(FT)method is presented for processing wavefield data of Lamb waves.The results can be used to determine the location,size and depth of rectangular notch.The analysis techniques are demonstrated using simulation examples of an aluminium plate with a rectangular notch.Then,the wavenumber analysis method is applied to simulation data that are obtained through a range of notch depths and widths.The results are analyzed and rules of the technique with regards to estimating notch depth are determined.Based on simulation results,guidelines for using the technique are developed.Finally,experimental wavefield data are obtained in aluminium plates with rectangular notches by a full noncontact transceiving method,i.e.,laser-laser method.Band-pass filtering combined with continuous wavelet transform is used to extract a certain frequency component from the full laser-induced wavefield with wide band.Shortspace two-dimensional FT method is used for further processing full wavefield data at a certain frequency to estimate spatially dependent wavenumber values.The consistency of simulation and experimental results shows the effectiveness of proposed wavenumber method for quantitative
基金supported by the Lyon Acoustics Center of Lyon University,Francefunded by the China Scholarship Council(CSC)。
文摘The reliable estimation of the wavenumber space(k-space)of the plates remains a longterm concern for acoustic modeling and structural dynamic behavior characterization.Most current analyses of wavenumber identification methods are based on the deterministic hypothesis.To this end,an inverse method is proposed for identifying wave propagation characteristics of twodimensional structures under stochastic conditions,such as wavenumber space,dispersion curves,and band gaps.The proposed method is developed based on an algebraic identification scheme in the polar coordinate system framework,thus named Algebraic K-Space Identification(AKSI)technique.Additionally,a model order estimation strategy and a wavenumber filter are proposed to ensure that AKSI is successfully applied.The main benefit of AKSI is that it is a reliable and fast method under four stochastic conditions:(A)High level of signal noise;(B)Small perturbation caused by uncertainties in measurement points’coordinates;(C)Non-periodic sampling;(D)Unknown structural periodicity.To validate the proposed method,we numerically benchmark AKSI and three other inverse methods to extract dispersion curves on three plates under stochastic conditions.One experiment is then performed on an isotropic steel plate.These investigations demonstrate that AKSI is a good in-situ k-space estimator under stochastic conditions.
基金supported by the National Special Plan for the 13th Five-Year Plan of China(No.2017YFC0602204-10)Independent Exploration of the Innovation Project for Graduate Students at Central South University(No.2017zzts176)+3 种基金National Natural Science Foundation of China(Nos.41574127,41404106,and 41674075)Postdoctoral Fund Projects of China(No.2017M622608)National Key R&D Program of China(No.2018YFC0603602)Natural Science Youth Fund Project of the Hunan Province,China(No.2018JJ3642)
文摘Magnetic field gradient tensor technique provides abundant data for delicate inversion of subsurface magnetic susceptibility distribution. Large scale magnetic data inversion imaging requires high speed and accuracy for forward modeling. For arbitrarily distributed susceptibility data on an undulated surface, we propose a fast 3D forward modeling method in the wavenumber domain based on(1) the wavenumber-domain expression of the prism combination model and the Gauss–FFT algorithm and(2) cubic spline interpolation. We apply the proposed 3D forward modeling method to synthetic data and use weighting coefficients in the wavenumber domain to improve the modeling for multiple observation surfaces, and also demonstrate the accuracy and efficiency of the proposed method.
基金supported by the CAS/SAFEA International Partnership Program for Creative Research Teams(KZZD-EW-TZ-19)
文摘It is well known that quantitative estimation of slip distributions on fault plane is one of the most important issues for earthquake source inversion related to the fault rupture process. The characteristics of slip distribution on the main fault play a fundamental role to control strong ground motion pattern. A large amount of works have also suggested that variable slip models inverted from longer period ground motion recordings are relevant for the prediction of higher frequency ground motions. Zhang et al. （Chin J Geophys 56：1412-1417, 2013） and Wang et al. （Chin J Geophys 56：1408-1411,2013） published their source inversions for the fault rupturing process soon after the April 20, 2013 Lushan earthquake in Sichuan, China. In this study, first, we synthesize two forward source slip models： the value of maximum slip, fault dimension, size, and dimension of major asperities, and comer wave number obtained from Wang＇s model is adopted to constrain the gen- eration of k-2 model and crack model. Next, both inverted and synthetic slip models are employed to simulate the ground motions for the Lushan earthquake based on the stochastic finite-fault method. In addition, for a comparison purpose, a stochastic slip model and another k-2 model （k 2 model II） with 2 times value of comer wave number of the original k-2 model （k 2 model I） are also employed for simulation for Lushan event. The simulated results characterized by Modified Mer- calli Intensity （MMI） show that the source slip models based on the inverted and synthetic slip distributions could capture many basic features associated with the ground motion patterns. Moreover, the simulated MMI distributions reflect the rupture directivity effect and the influence of the shallow velocity structure well. On the other hand, the simulated MMI bystochastic slip model and k 2 model II is apparently higher than observed intensity. By contrast, our simulation results show that the higher frequency ground motion is sensitive to the degree of sli
基金This work was supported by the National Key Research and Development Program of China （Nos. 2017YFC0601305, 2017YFC0602203, and 2017YFC0601606）, National Science and Technology Major Project task （No. 2016ZX05027-002-03）, National Natural Science Foundation of China （No. 41604098）, and State Key Program of National Natural Science of China （No. 41430322）.
文摘The current local wavenumber methods for the interpretation of magnetic anomalies compute the locations of geological bodies by solving complex matrices. Presently, such methods require to know the structural index, which is a parameter that represents the source type. The structural index is hard to know in real data; consequently, the precision of current methods is low. We present the fast local wavenumber （FLW） method, and define the squared sum of the horizontal and vertical local wavenumbers as the cumulative local wavenumber. The FLW method is the linear combination of the umulative local wavenumberand other wavenumbers, and is used to compute the locations and structural index of the source without a priori information and matrix solution. We apply the FLW method to synthetic magnetic anomalies, and the results suggest that the FLW method is insensitive to background and oblique magnetization. Next, we apply the FLW method to real magnetic data to obtain the location and structural index of the source.
文摘A new finite difference scheme-SCSD scheme has been proposed based on CD (Central Difference)scheme and SUD (Secondr-order Upwind Difference) scheme. Its basic feature is controIlable convectivesfability and always second-order accuracy (Stability-Controllable Second-order Difference ). It hasbeen proven that this scheme is convective-stable if the grid Peclet number .The advanage of this new scheme has been discussed based on the modified wavenumber analysis byusing Fourier transform. This scheme has been applied to the 2-D incompressible convective-diffusiveequation and 2-D compressible Euler equation, and corresponding finite difference equations have beenderived. Numerical examples of 1-D Burgers equation and 2-D transport equation have been presentedto show its good accuracy and controllable convective stability
