Wavefield separation of multicomponent seismic data to image subsurface structures can be realized in either the space domain or the wavenumber domain. However, as the particle velocity components used in the wavenumb...Wavefield separation of multicomponent seismic data to image subsurface structures can be realized in either the space domain or the wavenumber domain. However, as the particle velocity components used in the wavenumber-domain wavefield separation are not defined at the same grid point with the staggered-grid finite-difference method for elastic wavefield simulation, we propose the wavenumber-domain interpolation method to estimate the required values at the common grid points prior to the wavenumber-domain true-amplitude wavefield separation. Moreover, numerical experiments show that the wavenumber-domain interpolation method has high interpolation accuracy and the trueamplitude wavefield separation method shows good amplitude preservation. The application of the proposed methodology to elastic reverse-time migration can obtain good amplitudepreserved images even in the case of some velocity error.展开更多
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 rectangular notch detection.展开更多
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.展开更多
The wavenumber spectral components WN4 at the mesosphere and low thermosphere(MLT)altitudes(70–10 km)and in the latitude range between±45°are obtained from temperature data(T)observed by the Sounding of the...The wavenumber spectral components WN4 at the mesosphere and low thermosphere(MLT)altitudes(70–10 km)and in the latitude range between±45°are obtained from temperature data(T)observed by the Sounding of the Atmosphere using Broadband Emission Radiometry(SABER)instruments on board the National Aeronautics and Space Administration(NASA)’s Thermosphere–Ionosphere–Mesosphere Energetics and Dynamics(TIMED)spacecraft during the 11-year solar period from 2002 to 2012.We analyze in detail these spectral components WNk and obtain the main properties of their vertical profiles and global structures.We report that all of the wavenumber spectral components WNk occur mainly around 100 km altitude,and that the most prominent component is the wavenumber spectral component WN4 structure.Comparing these long duration temperature data with results of previous investigations,we have found that the yearly variation of spectral component WN4 is similar to that of the eastward propagating non-migrating diurnal tide with zonal wavenumber 3(DE3)at the low latitudes,and to that of the semi-diurnal tide with zonal wavenumber 2(SE2)at the mid-latitudes:the amplitudes of the A4 are larger during boreal summer and autumn at the low-latitudes;at the mid-latitudes the amplitudes have a weak peak in March.In addition,the amplitudes of component WN4 undergo a remarkable short period variation:significant day-to-day variation of the spectral amplitudes A4 occurs primarily in July and September at the low-latitudes.In summary,we conclude that the non-migrating tides DE3 and SE2 are likely to be the origins,at the low-latitudes and the mid-latitudes in the MLT region,respectively,of the observed wavenumber spectral component WN4.展开更多
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.展开更多
基金supported by the National Science Foundation of China(No.41174100)the Large-scale Oil and Gas Field and Coalbed Methane Development Major Projects(No.2011ZX05019-008-08)the China National Petroleum Corporation(No.2014A-3609)
文摘Wavefield separation of multicomponent seismic data to image subsurface structures can be realized in either the space domain or the wavenumber domain. However, as the particle velocity components used in the wavenumber-domain wavefield separation are not defined at the same grid point with the staggered-grid finite-difference method for elastic wavefield simulation, we propose the wavenumber-domain interpolation method to estimate the required values at the common grid points prior to the wavenumber-domain true-amplitude wavefield separation. Moreover, numerical experiments show that the wavenumber-domain interpolation method has high interpolation accuracy and the trueamplitude wavefield separation method shows good amplitude preservation. The application of the proposed methodology to elastic reverse-time migration can obtain good amplitudepreserved images even in the case of some velocity error.
基金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 rectangular notch detection.
基金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 present work is supported by National Science Foundation of China(41604138,41427901,41621063,41474133,41674158,41874179,41322030).
文摘The wavenumber spectral components WN4 at the mesosphere and low thermosphere(MLT)altitudes(70–10 km)and in the latitude range between±45°are obtained from temperature data(T)observed by the Sounding of the Atmosphere using Broadband Emission Radiometry(SABER)instruments on board the National Aeronautics and Space Administration(NASA)’s Thermosphere–Ionosphere–Mesosphere Energetics and Dynamics(TIMED)spacecraft during the 11-year solar period from 2002 to 2012.We analyze in detail these spectral components WNk and obtain the main properties of their vertical profiles and global structures.We report that all of the wavenumber spectral components WNk occur mainly around 100 km altitude,and that the most prominent component is the wavenumber spectral component WN4 structure.Comparing these long duration temperature data with results of previous investigations,we have found that the yearly variation of spectral component WN4 is similar to that of the eastward propagating non-migrating diurnal tide with zonal wavenumber 3(DE3)at the low latitudes,and to that of the semi-diurnal tide with zonal wavenumber 2(SE2)at the mid-latitudes:the amplitudes of the A4 are larger during boreal summer and autumn at the low-latitudes;at the mid-latitudes the amplitudes have a weak peak in March.In addition,the amplitudes of component WN4 undergo a remarkable short period variation:significant day-to-day variation of the spectral amplitudes A4 occurs primarily in July and September at the low-latitudes.In summary,we conclude that the non-migrating tides DE3 and SE2 are likely to be the origins,at the low-latitudes and the mid-latitudes in the MLT region,respectively,of the observed wavenumber spectral component WN4.
基金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.