This study investigates seismic interferometry in which the Green＇s function is estimated between two receiv- ers by cross-correlation and integration over sources. For smoothly varying source strengths, the dominant...This study investigates seismic interferometry in which the Green＇s function is estimated between two receiv- ers by cross-correlation and integration over sources. For smoothly varying source strengths, the dominant contributions of the correlation integral come from the stationary phase directions in the forward and backward directions from the alignment of the two receivers. Gaussian beams can be used to evaluate the correlation integral and concentrate the amplitudes in a vicinity of the stationary phase regions instead of completely relying on phase interference. Several numerical examples are shown to illustrate how this process works. The use of Gaussian beams for the evaluation of the correlation integral results in stable estimates, and also provides physical insight into the estimation of the Green＇s function based on seismic interferometry.展开更多
Studying the inner structure of intraplate earthquakes originating in aseismic areas,which are poorly covered by seismic networks or as historical earthquakes is usually the only way to get knowledge about their sourc...Studying the inner structure of intraplate earthquakes originating in aseismic areas,which are poorly covered by seismic networks or as historical earthquakes is usually the only way to get knowledge about their source mechanisms,which is partially essential for a deeper understanding of intraplate geodynamics.The epicentral zones of earthquakes are situated in hard-toreach areas,so,using active seismic methods for such purposes is unreasonable or even impossible because of high cost and logistical difficulties.We propose a novel approach that combines diverse passive seismic methods,which allows us to get sufficient information about geological environment structure for such task solutions.As an example,we investigated the inner structure of platform earthquake epicentral zone originated up north of Russia.We used four passive seismic methods:microseismic sounding method,passive seismic interferometry,HVSR method,and microseismic activity method.We show that passive seismic data,recoded in the same installation and processed by these different methods,can provide sufficient information about structure of studied environment,needed to explain source mechanism.In sum,the hypocenter zone is presented by intersection of vertical faults and a lateral fractured zone in the middle crust.Results were confirmed by comparison with results by active seismic methods.展开更多
In common practice in the oil fields,the injection of water and gas into reservoirs is a crucial technique to increase production.The control of the waterflooding front in oil/gas exploitation is a matter of great con...In common practice in the oil fields,the injection of water and gas into reservoirs is a crucial technique to increase production.The control of the waterflooding front in oil/gas exploitation is a matter of great concern to reservoir engineers.Monitoring the waterflooding front in oil/gas wells plays a very important role in adjusting the well network and later in production,taking advantage of the remaining oil po-tential and ultimately achieving great success in improving the recovery rate.For a long time,micro-seismic monitoring,numerical simulation,four-dimensional seismic and other methods have been widely used in waterflooding front monitoring.However,reconciling their reliability and cost poses a significant challenge.In order to achieve real-time,reliable and cost-effective monitoring,we propose an innovative method for waterflooding front monitoring through the similarity analysis of passive source time-lapse seismic images.Typically,passive source seismic data collected from oil fields have extremely low signal-to-noise ratio(SNR),which poses a serious problem for obtaining structural images.The proposed method aims to visualize and analyze underground changes by highlighting time-lapse images and provide a strategy for underground monitoring using long-term passive source data under low SNR conditions.First,we verify the feasibility of the proposed method by designing a theoretical model.Then,we conduct an analysis of the correlation coefficient(similarity)on the passive source time-lapse seismic imaging results to enhance the image differences and identify the simulated waterflooding fronts.Finally,the proposed method is applied to the actual waterflooding front monitoring tasks in Shengli Oilfield,China.The research findings indicate that the monitoring results are consistent with the actual devel-opment conditions,which in turn demonstrates that the proposed method has great potential for practical application and is very suitable for monitoring common development tasks in oil fields.展开更多
A case study of seismic interferometry applied to a small microseismic monitoring network is here presented.The main objectives of this study are(i)to quantify the lateral variability of shear-wave ve-locities in the ...A case study of seismic interferometry applied to a small microseismic monitoring network is here presented.The main objectives of this study are(i)to quantify the lateral variability of shear-wave ve-locities in the studied area,and(ii)to investigate the bias produced by noise directionality and non-stationarity in the velocity estimate.Despite the limited number of stations and the short-period char-acter of the seismic sensors,the empirical Green's functions were retrieved for all station pairs using two years of passive data.Both group and phase velocities were derived,the former using the widespread frequency-time analysis,the latter through the analysis of the real part of the cross-spectra.The main advantage of combining these two methods is a more accurate identification of higher modes,resulting in a reduction of ambiguity during picking and data interpretation.Surface wave tomography was run to obtain the spatial distribution of group and phase velocities for the same wavelengths.The low standard deviation of the results suggests that the sparse character of the network does not limit the applicability of the method,for this specific case.The obtained maps highlight the presence of a lower velocity area that extends from the centre of the network towards southeast.Group and phase velocity dispersion curves have been jointly inverted to retrieve as many shear-wave velocity profiles as selected station pairs.While the average model can be used for a more accurate location of the local natural seismicity,the associated standard deviations give us an indication of the lateral heterogeneity of seismic velocities as a function of depth.Finally,the same velocity analysis was repeated for different time windows in order to quantify the error associated to variations in the noise field.Errors as large as 4%have been found,related to the unfavorable orientation of the receiver pairs with respect to strongly directional noise sources,and to the very short time widows.It was shown that using a one-year time window these errors arereduced to 0.3%.展开更多
The thickness and upper densification structure of an ice sheet are important parameters for dynamic ice sheet modeling and glacier mass balance studies.Seismic ambient noise methods,such as the horizontal-to-vertical...The thickness and upper densification structure of an ice sheet are important parameters for dynamic ice sheet modeling and glacier mass balance studies.Seismic ambient noise methods,such as the horizontal-to-vertical spectral ratio(H/V)method and ambient noise cross-correlation method,are becoming increasingly popular in glacier structure investigations.During China's 39th expedition to Antarctica,seismic ambient noise experiments were conducted to investigate the structure of the ice sheet at Kunlun Station,Dome A,using a seismic nodal system.We obtained a broad band(0.1–10 Hz)H/V curve with a 1-hour noise record from a seismic node.In addition,we extracted the Rayleigh wave dispersion curve with 5-day noise cross-correlation functions from a linear dense seismic array.Three clear peaks were observed in the H/V curve—a lower peak at~0.17 Hz and two higher frequency peaks at~3 Hz and~6 Hz.We inverted the ice sheet thickness using the lower frequency portion of the H/V curve and inverted the upper structure of the ice sheet using the higher frequency portion of the H/V curve jointly with the dispersion curve.Our estimations from ambient noise observations were consistent with those derived from the BedMachine ice sheet thickness dataset and the density profile determined by ground-penetrating radar investigations at the same site.展开更多
Passive surface wave imaging has been a powerful tool for near-surface characterization in urban areas,which extracts surface wave signals from ambient seismic noise and then estimates subsurface shear wave velocity b...Passive surface wave imaging has been a powerful tool for near-surface characterization in urban areas,which extracts surface wave signals from ambient seismic noise and then estimates subsurface shear wave velocity by inversion of the measured phase velocity.The high-frequency(approximately>1 Hz)seismic noise fields in urban environments are dominantly induced by human activities such as the vehicle traffic.Traffic seismic sources are nonrandomly distributed in time and space.Applying standard interferometric techniques to recordings from these nonrandom noise sources makes the Green’s function liable to estimation errors.We analyze the influence of using nonrandom traffic seismic sources for surface wave imaging.With nonrandom traffic seismic sources in time,spurious signals are generated in the cross-correlation function.With nonrandom traffic seismic sources in space,surface-wave phase velocities could be overestimated in the dispersion measurement.We provide an overview of solutions for surface-wave imaging with nonrandom traffic seismic sources in time and space,aiming to improve the retrieval of high-frequency surface waves and achieve reliable results from ultrashort(tens of seconds)observations for near-surface characterization.展开更多
Dense-array ambient noise tomography is a powerful tool for achieving high-resolution subsurface imag-ing,significantly impacting geohazard prevention and control.Conventional dense-array studies,how-ever,require simu...Dense-array ambient noise tomography is a powerful tool for achieving high-resolution subsurface imag-ing,significantly impacting geohazard prevention and control.Conventional dense-array studies,how-ever,require simultaneous observations of numerous stations for extensive coverage.To conduct a comprehensive karst feature investigation with limited stations,we designed a new synchronous-asyn-chronous observation system that facilitates dense array observations.We conducted two rounds of asynchronous observations,each lasting approximately 24 h,in combination with synchronous backbone stations.We achieved wide-ranging coverage of the study area utilizing 197 nodal receivers,with an average station spacing of 7 m.The beamforming results revealed distinct variations in the noise source distributions between day and night.We estimated the source strength in the stationary phase zone and used a weighting scheme for stacking the cross-correlation functions(C ^(1) functions)to suppress the influ-ence of nonuniform noise source distributions.The weights were derived from the similarity coefficients between multicomponent C^(1)functions related to Rayleigh waves.We employed the cross-correlation of C ^(1) functions(C^(2)methods)to obtain the empirical Green’s functions between asynchronous stations.To eliminate artifacts in C ^(2) functions from higher-mode surface waves in C^(1)functions,we filtered the C^(1)functions on the basis of different particle motions linked to multimode Rayleigh waves.The dispersion measurements of Rayleigh waves obtained from both the C^(1)and C^(2)functions were utilized in surface wave tomography.The inverted three-dimensional(3D)shear-wave(S-wave)velocity model reveals two significant low-velocity zones at depths ranging from 40 to 60 m,which align well with the karst caves found in the drilling data.The method of short-term synchronous-asynchronous ambient noise tomography shows promise as a cost-effective and efficient approach for urban geohazard investigations.展开更多
Multiply scattered waves are sensitive to media changes owing to the effect of repeated sampling,superposition and amplification. Based on this characteristic,small-medium changes could be detected by using coda wave ...Multiply scattered waves are sensitive to media changes owing to the effect of repeated sampling,superposition and amplification. Based on this characteristic,small-medium changes could be detected by using coda wave interferometry. In recent years,coda wave interferometry has been widely used in estimating velocity variation with high precision in areas such as seismology and non-destructive testing. This paper systematically presents the principle and research status of coda wave interferometry,and especial focus is placed on the research of media velocity variations by using repeating earthquakes,artificial sources,and ambient noise. Applications of coda wave interferometry can contribute to the more subtle understanding of dynamic evolution process in the medium.展开更多
In conventional seismic exploration,complex overburden prevents us from imaging the deeper underground structure,because near-surface velocity model cannot be efficiently simplified and accurately established. The vir...In conventional seismic exploration,complex overburden prevents us from imaging the deeper underground structure,because near-surface velocity model cannot be efficiently simplified and accurately established. The virtual source method provides a promising solution,which images below a complex overburden without the knowledge of overburden velocities and near-surface changes. In this paper,the authors simulate the forward modeling under complex overburden with sources on the surface and receivers in the horizontal well,and reset the real sources into the horizontal well using the virtual source method,which subtly avoid the distortion effects of the complex overburden. Finally the underground structure can be imaged using Kirchhoff migration,so the feasibility and advantages of the virtual source method are verified.展开更多
The extraction of the Earth＇s Green＇s function from field fluctuations is a rapidly growing area of research. The principle of Green＇s function extraction is often related to the requirement of equipartitioning, wh...The extraction of the Earth＇s Green＇s function from field fluctuations is a rapidly growing area of research. The principle of Green＇s function extraction is often related to the requirement of equipartitioning, which stipulates that the energy of field fluctuations is distributed evenly in some sense. We show the meaning of equipartitioning for a variety of different formulations for Green＇s function retrieval. We show that equipartitioning is not a sufficient condition, and provide several examples that illustrate this point. We discuss the implications of lack of equipartitioning for various schemes for the reconstruction of the Green＇s function in seismology. The theory for Green＇s function extraction is usually based on a statistical theory that relies on ensemble averages. Since there is only one Earth, one usually replaces the ensemble average with a time average. We show that such a replacement only makes sense when attenuation is taken into account, and show how the theory for Green＇s function extraction for oscillating systems can be extended to incorporate attenuation.展开更多
基金supported by U.S. National Science Foundation EAR06-35611U.S. Air Force contract FA8718-08-C-002the members of the Geo-Mathematical Imaging Group (GMIG) at Purdue University
文摘This study investigates seismic interferometry in which the Green＇s function is estimated between two receiv- ers by cross-correlation and integration over sources. For smoothly varying source strengths, the dominant contributions of the correlation integral come from the stationary phase directions in the forward and backward directions from the alignment of the two receivers. Gaussian beams can be used to evaluate the correlation integral and concentrate the amplitudes in a vicinity of the stationary phase regions instead of completely relying on phase interference. Several numerical examples are shown to illustrate how this process works. The use of Gaussian beams for the evaluation of the correlation integral results in stable estimates, and also provides physical insight into the estimation of the Green＇s function based on seismic interferometry.
基金funded by the Russian Federation Ministry of Science and Higher Education according to the research project FUUW-2025-0033 of the FECIAR UrB RASthe Arkhangelsk branch of the Territorial Fund of Geological Information for the North-Western Federal District for providing geological and geophysical reports for the objects under consideration。
文摘Studying the inner structure of intraplate earthquakes originating in aseismic areas,which are poorly covered by seismic networks or as historical earthquakes is usually the only way to get knowledge about their source mechanisms,which is partially essential for a deeper understanding of intraplate geodynamics.The epicentral zones of earthquakes are situated in hard-toreach areas,so,using active seismic methods for such purposes is unreasonable or even impossible because of high cost and logistical difficulties.We propose a novel approach that combines diverse passive seismic methods,which allows us to get sufficient information about geological environment structure for such task solutions.As an example,we investigated the inner structure of platform earthquake epicentral zone originated up north of Russia.We used four passive seismic methods:microseismic sounding method,passive seismic interferometry,HVSR method,and microseismic activity method.We show that passive seismic data,recoded in the same installation and processed by these different methods,can provide sufficient information about structure of studied environment,needed to explain source mechanism.In sum,the hypocenter zone is presented by intersection of vertical faults and a lateral fractured zone in the middle crust.Results were confirmed by comparison with results by active seismic methods.
基金supported by the CNPC-SWPU Innovation Alliance Technology Cooperation Project(2020CX020000)the National Natural Science Foundation of China(42022028)+1 种基金the Natural Science Foundation of Sichuan Province(24NSFSC0808)the China Scholarship Council(202306440144)。
文摘In common practice in the oil fields,the injection of water and gas into reservoirs is a crucial technique to increase production.The control of the waterflooding front in oil/gas exploitation is a matter of great concern to reservoir engineers.Monitoring the waterflooding front in oil/gas wells plays a very important role in adjusting the well network and later in production,taking advantage of the remaining oil po-tential and ultimately achieving great success in improving the recovery rate.For a long time,micro-seismic monitoring,numerical simulation,four-dimensional seismic and other methods have been widely used in waterflooding front monitoring.However,reconciling their reliability and cost poses a significant challenge.In order to achieve real-time,reliable and cost-effective monitoring,we propose an innovative method for waterflooding front monitoring through the similarity analysis of passive source time-lapse seismic images.Typically,passive source seismic data collected from oil fields have extremely low signal-to-noise ratio(SNR),which poses a serious problem for obtaining structural images.The proposed method aims to visualize and analyze underground changes by highlighting time-lapse images and provide a strategy for underground monitoring using long-term passive source data under low SNR conditions.First,we verify the feasibility of the proposed method by designing a theoretical model.Then,we conduct an analysis of the correlation coefficient(similarity)on the passive source time-lapse seismic imaging results to enhance the image differences and identify the simulated waterflooding fronts.Finally,the proposed method is applied to the actual waterflooding front monitoring tasks in Shengli Oilfield,China.The research findings indicate that the monitoring results are consistent with the actual devel-opment conditions,which in turn demonstrates that the proposed method has great potential for practical application and is very suitable for monitoring common development tasks in oil fields.
基金This study is part of the project 2021RUAPON-REACT EU-Finanziamento PON“Ricerca e Innovazione”20142020,grant n.19-G-12543-2,funded by the Italian Ministry of University and Research(MUR)This study was developed in the frame of“The Geosciences for Sustainable Development”project(Budget Minis-tero dell'Universita e della Ricerca-Dipartimenti di Eccellenza 2023-2027,code n.C93C23002690001).
文摘A case study of seismic interferometry applied to a small microseismic monitoring network is here presented.The main objectives of this study are(i)to quantify the lateral variability of shear-wave ve-locities in the studied area,and(ii)to investigate the bias produced by noise directionality and non-stationarity in the velocity estimate.Despite the limited number of stations and the short-period char-acter of the seismic sensors,the empirical Green's functions were retrieved for all station pairs using two years of passive data.Both group and phase velocities were derived,the former using the widespread frequency-time analysis,the latter through the analysis of the real part of the cross-spectra.The main advantage of combining these two methods is a more accurate identification of higher modes,resulting in a reduction of ambiguity during picking and data interpretation.Surface wave tomography was run to obtain the spatial distribution of group and phase velocities for the same wavelengths.The low standard deviation of the results suggests that the sparse character of the network does not limit the applicability of the method,for this specific case.The obtained maps highlight the presence of a lower velocity area that extends from the centre of the network towards southeast.Group and phase velocity dispersion curves have been jointly inverted to retrieve as many shear-wave velocity profiles as selected station pairs.While the average model can be used for a more accurate location of the local natural seismicity,the associated standard deviations give us an indication of the lateral heterogeneity of seismic velocities as a function of depth.Finally,the same velocity analysis was repeated for different time windows in order to quantify the error associated to variations in the noise field.Errors as large as 4%have been found,related to the unfavorable orientation of the receiver pairs with respect to strongly directional noise sources,and to the very short time widows.It was shown that using a one-year time window these errors arereduced to 0.3%.
基金supported by the National Natural Science Foundation of China(NSFC),grant number 42076234the National Key Research and Development Program of China,grant number 2021YFC2801404。
文摘The thickness and upper densification structure of an ice sheet are important parameters for dynamic ice sheet modeling and glacier mass balance studies.Seismic ambient noise methods,such as the horizontal-to-vertical spectral ratio(H/V)method and ambient noise cross-correlation method,are becoming increasingly popular in glacier structure investigations.During China's 39th expedition to Antarctica,seismic ambient noise experiments were conducted to investigate the structure of the ice sheet at Kunlun Station,Dome A,using a seismic nodal system.We obtained a broad band(0.1–10 Hz)H/V curve with a 1-hour noise record from a seismic node.In addition,we extracted the Rayleigh wave dispersion curve with 5-day noise cross-correlation functions from a linear dense seismic array.Three clear peaks were observed in the H/V curve—a lower peak at~0.17 Hz and two higher frequency peaks at~3 Hz and~6 Hz.We inverted the ice sheet thickness using the lower frequency portion of the H/V curve and inverted the upper structure of the ice sheet using the higher frequency portion of the H/V curve jointly with the dispersion curve.Our estimations from ambient noise observations were consistent with those derived from the BedMachine ice sheet thickness dataset and the density profile determined by ground-penetrating radar investigations at the same site.
基金supported by the National Natural Science Foundation of China(No.42374177)。
文摘Passive surface wave imaging has been a powerful tool for near-surface characterization in urban areas,which extracts surface wave signals from ambient seismic noise and then estimates subsurface shear wave velocity by inversion of the measured phase velocity.The high-frequency(approximately>1 Hz)seismic noise fields in urban environments are dominantly induced by human activities such as the vehicle traffic.Traffic seismic sources are nonrandomly distributed in time and space.Applying standard interferometric techniques to recordings from these nonrandom noise sources makes the Green’s function liable to estimation errors.We analyze the influence of using nonrandom traffic seismic sources for surface wave imaging.With nonrandom traffic seismic sources in time,spurious signals are generated in the cross-correlation function.With nonrandom traffic seismic sources in space,surface-wave phase velocities could be overestimated in the dispersion measurement.We provide an overview of solutions for surface-wave imaging with nonrandom traffic seismic sources in time and space,aiming to improve the retrieval of high-frequency surface waves and achieve reliable results from ultrashort(tens of seconds)observations for near-surface characterization.
基金supported by the National Natural Science Foundation of China(41830103)the Project of Nanjing Center of China Geological Survey(DD20190281).
文摘Dense-array ambient noise tomography is a powerful tool for achieving high-resolution subsurface imag-ing,significantly impacting geohazard prevention and control.Conventional dense-array studies,how-ever,require simultaneous observations of numerous stations for extensive coverage.To conduct a comprehensive karst feature investigation with limited stations,we designed a new synchronous-asyn-chronous observation system that facilitates dense array observations.We conducted two rounds of asynchronous observations,each lasting approximately 24 h,in combination with synchronous backbone stations.We achieved wide-ranging coverage of the study area utilizing 197 nodal receivers,with an average station spacing of 7 m.The beamforming results revealed distinct variations in the noise source distributions between day and night.We estimated the source strength in the stationary phase zone and used a weighting scheme for stacking the cross-correlation functions(C ^(1) functions)to suppress the influ-ence of nonuniform noise source distributions.The weights were derived from the similarity coefficients between multicomponent C^(1)functions related to Rayleigh waves.We employed the cross-correlation of C ^(1) functions(C^(2)methods)to obtain the empirical Green’s functions between asynchronous stations.To eliminate artifacts in C ^(2) functions from higher-mode surface waves in C^(1)functions,we filtered the C^(1)functions on the basis of different particle motions linked to multimode Rayleigh waves.The dispersion measurements of Rayleigh waves obtained from both the C^(1)and C^(2)functions were utilized in surface wave tomography.The inverted three-dimensional(3D)shear-wave(S-wave)velocity model reveals two significant low-velocity zones at depths ranging from 40 to 60 m,which align well with the karst caves found in the drilling data.The method of short-term synchronous-asynchronous ambient noise tomography shows promise as a cost-effective and efficient approach for urban geohazard investigations.
基金sponsored by the Task of Department of Earthquake Monitoring and Prediction,CEAthe Special Fund of Earthquake Scientific Research of China Earthquake Administration(201208003)
文摘Multiply scattered waves are sensitive to media changes owing to the effect of repeated sampling,superposition and amplification. Based on this characteristic,small-medium changes could be detected by using coda wave interferometry. In recent years,coda wave interferometry has been widely used in estimating velocity variation with high precision in areas such as seismology and non-destructive testing. This paper systematically presents the principle and research status of coda wave interferometry,and especial focus is placed on the research of media velocity variations by using repeating earthquakes,artificial sources,and ambient noise. Applications of coda wave interferometry can contribute to the more subtle understanding of dynamic evolution process in the medium.
文摘In conventional seismic exploration,complex overburden prevents us from imaging the deeper underground structure,because near-surface velocity model cannot be efficiently simplified and accurately established. The virtual source method provides a promising solution,which images below a complex overburden without the knowledge of overburden velocities and near-surface changes. In this paper,the authors simulate the forward modeling under complex overburden with sources on the surface and receivers in the horizontal well,and reset the real sources into the horizontal well using the virtual source method,which subtly avoid the distortion effects of the complex overburden. Finally the underground structure can be imaged using Kirchhoff migration,so the feasibility and advantages of the virtual source method are verified.
文摘The extraction of the Earth＇s Green＇s function from field fluctuations is a rapidly growing area of research. The principle of Green＇s function extraction is often related to the requirement of equipartitioning, which stipulates that the energy of field fluctuations is distributed evenly in some sense. We show the meaning of equipartitioning for a variety of different formulations for Green＇s function retrieval. We show that equipartitioning is not a sufficient condition, and provide several examples that illustrate this point. We discuss the implications of lack of equipartitioning for various schemes for the reconstruction of the Green＇s function in seismology. The theory for Green＇s function extraction is usually based on a statistical theory that relies on ensemble averages. Since there is only one Earth, one usually replaces the ensemble average with a time average. We show that such a replacement only makes sense when attenuation is taken into account, and show how the theory for Green＇s function extraction for oscillating systems can be extended to incorporate attenuation.
