The Tan-Lu Fault Zone is a large NNE-trending fault zone that has a substantial effect on the development of eastern China and its earthquake disaster prevention efforts. Aiming at the azimuthally anisotropic structur...The Tan-Lu Fault Zone is a large NNE-trending fault zone that has a substantial effect on the development of eastern China and its earthquake disaster prevention efforts. Aiming at the azimuthally anisotropic structure in the upper crust and seismogenic tectonics in the Hefei segment of this fault, we collected phase velocity dispersion data of fundamental mode Rayleigh waves from ambient noise cross-correlation functions of ~400 temporal seismographs in an area of approximately 80 × 70 km along the fault zone. The period band of the dispersion data was ~0.5–10 s. We inverted for the upper crustal three-dimensional(3-D) shear velocity model with azimuthal anisotropy from the surface to 10 km depth by using a 3-D direct azimuthal anisotropy inversion method. The inversion result shows the spatial distribution characteristics of the tectonic units in the upper crust. Additionally, the deformation of the Tan-Lu Fault Zone and its conjugated fault systems could be inferred from the anisotropy model. In particular, the faults that have remained active from the early and middle Pleistocene control the anisotropic characteristics of the upper crustal structure in this area. The direction of fast axes near the fault zone area in the upper crust is consistent with the strike of the faults, whereas for the region far away from the fault zone, the direction of fast axes is consistent with the direction of the regional principal stress caused by plate movement. Combined with the azimuthal anisotropy models in the deep crust and uppermost mantle from the surface wave and Pn wave, the different anisotropic patterns caused by the Tan-Lu Fault Zone and its conjugated fault system nearby are shown in the upper and lower crust. Furthermore,by using the double-difference method, we relocated the Lujiang earthquake series, which contained 32 earthquakes with a depth shallower than 10 km. Both the Vs model and earthquake relocation results indicate that earthquakes mostly occurred in the vicinity of structural boundaries with fractured media, with high-level development of cracks and small-scale faults jammed between more rigid areas.展开更多
Determining the orientation of in-situ stresses is crucial for various geoscience and engineering appli-cations.Conventional methods for estimating these stress orientations often depend on focal mechanism solutions(F...Determining the orientation of in-situ stresses is crucial for various geoscience and engineering appli-cations.Conventional methods for estimating these stress orientations often depend on focal mechanism solutions(FMSs)derived from earthquake data and formation micro-imager(FMI)data from well logs.However,these techniques can be costly,depth-inaccurate,and may lack spatial coverage.To address this issue,we introduce the use of three-dimensional(3D)seismic data(active sources)as a lateral constraint to approximate the 3D stress orientation field.Recognizing that both stress and fracture patterns are closely related to seismic velocity anisotropy,we derive the orientation of azimuthal anisotropy from multi-azimuth 3D seismic data to compensate for the lack of spatial stress orientation information.We apply our proposed workflow to a case study in the Weiyuan area of the Sichuan Basin,China,a region targeted for shale gas production.By integrating diverse datasets,including 3D seismic,earthquakes,and well logs,we develop a comprehensive 3D model of in-situ stress(orientations and magnitudes).Our results demonstrate that the estimated anisotropy orientations from 3D seismic data are consistent with the direction of maximum horizontal principal stress(SHmax)obtained from FMIs.We analyzed 12 earthquakes(magnitude>3)recorded between 2016 and 2020 for their FMSs and compressional axis(P-axis)orientations.The derived SHmax direction from our 3D stress model is 110°ES(East-South),which shows excellent agreement with the FMSs(within 3.96°).This close alignment validates the reliability and precision of our integrated method for predicting 3D SHmax orientations.展开更多
The characterization of subsurface formations via the analysis of near-wellbore velocity profiles represents a crucial method in geophysical exploration.This technique enables the evaluation of key parameters,includin...The characterization of subsurface formations via the analysis of near-wellbore velocity profiles represents a crucial method in geophysical exploration.This technique enables the evaluation of key parameters,including rock brittleness,wellbore stability,fracturing effects,and invasion extent,thereby enhancing comprehension of formation structures and informing exploration strategies.However,traditional near-wellbore formation velocity imaging methods exhibit two principal limitations.First,these methods lack azimuthal sensitivity,yielding results averaged across all directions.Second,they are computationally intensive and impractical for well-site environments.To address these drawbacks,we developed a rapid 3D velocity imaging method for array acoustic logging instru ments equipped with azimuthal receivers,capable of producing 3D imaging results efficiently.The workflow entails the following steps:(1)Band-pass filtering of logging data to mitigate scattered wave interference caused by formation heterogeneity near the wellbore;(2)combination of receivers with varying detection ranges in each direction to derive radial velocity sequences,followed by integration of ray-tracing theory to obtain 2D velocity distributions;and(3)synthesis of final 3D velocity imaging results via interpolation of these 2D datasets.In the velocity sequence extraction process,we significantly reduced the computational load by employing an adaptive time window,ensuring rapid and stable application in well-site settings.We utilized the finite difference method to construct well models with heteroge neous formations.The compressional and shear wave 3D velocity imaging results derived from synthetic data correlated with the model,demonstrating the azimuthal sensitivity of our proposed method.Furthermore,we applied this method to a well in West China,successfully identifying the azimuth of nearwellbore anisotropy.展开更多
Reservoirs with a group of vertical fractures in a vertical transversely isotropic(VTI)background are considered as orthorhombic(ORT)medium.However,fracture detection in ORT medium using seismic inversion methods rema...Reservoirs with a group of vertical fractures in a vertical transversely isotropic(VTI)background are considered as orthorhombic(ORT)medium.However,fracture detection in ORT medium using seismic inversion methods remains challenging,as it requires the estimation of more than eight parameters.Assuming the reservoir to be a weakly anisotropic ORT medium with small contrasts in the background elastic parameters,a new azimuthal elastic impedance equation was first derived using parameter combinations and mathematical approximations.This equation exhibited almost the same accuracy as the original equation and contained only six model parameters:the compression modulus,anisotropic shear modulus,anisotropic compression modulus,density,normal fracture weakness,and tangential fracture weakness.Subsequently,a stepwise inversion method using second-order derivatives of the elastic impedance was developed to estimate these parameters.Moreover,the Thomsen anisotropy parameter,epsilon,was estimated from the inversion results using the ratio of the anisotropic compression modulus to the compression modulus.Synthetic examples with moderate noise and field data examples confirm the feasibility and effectiveness of the inversion method.The proposed method exhibited accuracy similar to that of previous inversion strategies and could predict richer vertical fracture information.Ultimately,the method was applied to a three-dimensional work area,and the predictions were consistent with logging and geological a priori information,confirming the effectiveness of this method.Summarily,the proposed stepwise inversion method can alleviate the uncertainty of multi-parameter inversion in ORT medium,thereby improving the reliability of fracture detection.展开更多
Range-azimuth imaging of ground targets via frequency-modulated continuous wave(FMCW)radar is crucial for effective target detection.However,when the pitch of the moving array constructed during motion exceeds the phy...Range-azimuth imaging of ground targets via frequency-modulated continuous wave(FMCW)radar is crucial for effective target detection.However,when the pitch of the moving array constructed during motion exceeds the physical array aperture,azimuth ambiguity occurs,making range-azimuth imaging on a moving platform challenging.To address this issue,we theoretically analyze azimuth ambiguity generation in sparse motion arrays and propose a dual-aperture adaptive processing(DAAP)method for suppressing azimuth ambiguity.This method combines spatial multiple-input multiple-output(MIMO)arrays with sparse motion arrays to achieve high-resolution range-azimuth imaging.In addition,an adaptive QR decomposition denoising method for sparse array signals based on iterative low-rank matrix approximation(LRMA)and regularized QR is proposed to preprocess sparse motion array signals.Simulations and experiments show that on a two-transmitter-four-receiver array,the signal-to-noise ratio(SNR)of the sparse motion array signal after noise suppression via adaptive QR decomposition can exceed 0 dB,and the azimuth ambiguity signal ratio(AASR)can be reduced to below-20 dB.展开更多
A P-band polarimetric synthetic aperture radar(PolSAR)sensor has deep penetration ability into and through the vegetation canopies in forested environments.Thus,the sensor is of great potential to accurately assess fo...A P-band polarimetric synthetic aperture radar(PolSAR)sensor has deep penetration ability into and through the vegetation canopies in forested environments.Thus,the sensor is of great potential to accurately assess forest parameters such as coverage,stand density,and tree height.Unfortunately,the radar backscatter from complex terrain can adversely impact the backscatter from trees or forests,and forest parameters assessed can be erroneous.Thus,reducing the topographic impact is an urgent must.In this study,a topographic compensation algorithm has been studied.To assess the algorithm’s validity and effectiveness,we applied it to P-band PolSAR datasets in four forested areas in the US.Trees in the forest stands have diverse species,and the topographic conditions of the terrain differ.Significant topographic impact on the P-band PolSAR data exists before the topographic compensation algorithm.After the algorithm,the impact decreases noticeably qualitatively and quantitatively.The algorithm is valid and effective in reducing the topographic influence on the PolSAR data and,consequently,provides a better chance of retrieving accurate forest parameters.展开更多
The purpose of this study is to analyze the impact of upper ocean dynamics on velocity bunching,represented by az-imuthal cutoff wavelength(i.e.,sea surface wind,wave,and current).In this study,over 1400 dual-polarize...The purpose of this study is to analyze the impact of upper ocean dynamics on velocity bunching,represented by az-imuthal cutoff wavelength(i.e.,sea surface wind,wave,and current).In this study,over 1400 dual-polarized(vertical-vertical(VV)and vertical-horizontal(VH))Sentinel-1(S-1)synthetic aperture radar(SAR)images collected in tropical cyclones(TC)are utilized.These images are combined with wind and rain observations from the stepped-frequency microwave radiometer(SFMR),wave simu-lations conducted using a third-generation numerical wave model,WAVEWATCH-III(WW3),and SAR-derived wind information collected from CyclObs winds.The WW3-simulated significant wave height(SWH)is validated against measurements from HY-2B altimeter taken in August and September 2021,yielding a root mean square error(RMSE)of 0.48 m and a correlation coefficient(COR)of 0.88.The SAR-based azimuthal cutoff wavelengths in VV polarization,which quantitatively represent the effect of velocity bunching,are compared with theoretical values calculated using WW3-simulated SWH.A notable relationship is observed between the difference in azimuthal cutoff wavelength and SAR-derived wind speed and WW3-simulated SWH.Analysis results show that the correlation between SAR-based azimuthal cutoff wavelength and SWH is stronger than that with wind and current.Finally,a machine learning algorithm is used to develop an algorithm aimed at simulating the azimuthal cutoff wavelength in TCs,including wind,wave,and incidence angle.This method yields an RMSE of 8.90 m,a COR of 0.91,and a scatter index of 0.04 for VV-polar-ization SAR.展开更多
During daylight laser polarization sensing of high-level clouds(HLCs),the lidar receiving system generates a signal caused by not only backscattered laser radiation,but also scattered solar radiation,the intensity and...During daylight laser polarization sensing of high-level clouds(HLCs),the lidar receiving system generates a signal caused by not only backscattered laser radiation,but also scattered solar radiation,the intensity and polarization of which depends on the Sun’s location.If a cloud contains spatially oriented ice particles,then it becomes anisotropic,that is,the coefficients of directional light scattering of such a cloud depend on the Sun’s zenith and azimuth angles.In this work,the possibility of using the effect of anisotropic scattering of solar radiation on the predictive ability of machine learning algorithms in solving the problem of predicting the HLC backscattering phase matrix(BSPM)was evaluated.The hypothesis that solar radiation scattered on HLCs has no effect on the BSPM elements of such clouds determined with a polarization lidar was tested.The operation of two algorithms for predicting the BSPM elements is evaluated.To train the first one,meteorological data were used as input parameters;for the second algorithm,the azi-muthal and zenith angles of the Sun’s position were added to the meteorological parameters.It is shown that there is no significant improvement in the predictive ability of the algorithm.展开更多
We propose a method for mOdeling azimuthal AVO responses from a fractured i reflector. The method calculates the integrated reflected wavetrains, and the wavetrains contain elastodynamic information including the con...We propose a method for mOdeling azimuthal AVO responses from a fractured i reflector. The method calculates the integrated reflected wavetrains, and the wavetrains contain elastodynamic information including the contrast in impedance and anelasticity i across interfaces, the intemal anisotropic propagation, the dispersion and attenuation along i the wave path, and tuning and interference. The results suggest that for large angles of incidence, the velocity dispersion and attenuation increase the amplitudes of PP waves from the top and decrease those from the bottom. For azimuthal responses at specific angles of incidence, the reflected wavetrains of PP waves tend to have longer duration with increasing azimuth. In contrast, model-converted PSV and PSH reflections show stable azimuthal features and are less affected by the reflector thickness. The amplitudes of PSV reflections increase with increasing azimuth; moreover, the waves have no reflection energy at 0° and 90° azimuth and maximum amplitude at 45° azimuth.展开更多
The effect of the wide and narrow azimuth 3D observation systems on seismic imaging precision is becoming a hot area for studies of high precision 3D seismic acquisition methods in recent years. In this paper we utili...The effect of the wide and narrow azimuth 3D observation systems on seismic imaging precision is becoming a hot area for studies of high precision 3D seismic acquisition methods in recent years. In this paper we utilize 3D physical seismic modeling experiments. A 3D channel sand body physical seismic model is constructed and two acquisition systems are designed with wide azimuth （16 lines） and narrow azimuth （8 lines） to model 3D seismic data acquisition and processing seismic work flows. From analysis of migrated time slice data with high quality and small size, we conclude that when the overlying layers are smooth and lateral velocities have little change, both wide and narrow azimuth observation systems in 3D acquisition can be used for obtaining high precision imaging and equivalent resolution of the channel sand body.展开更多
We developed a novel cement evaluation logging tool,named the azimuthally acoustic bond tool（AABT）,which uses a phased-arc array transmitter with azimuthal detection capability.We combined numerical simulations and ...We developed a novel cement evaluation logging tool,named the azimuthally acoustic bond tool（AABT）,which uses a phased-arc array transmitter with azimuthal detection capability.We combined numerical simulations and field tests to verify the AABT tool.The numerical simulation results showed that the radiation direction of the subarray corresponding to the maximum amplitude of the first arrival matches the azimuth of the channeling when it is behind the casing.With larger channeling size in the circumferential direction,the amplitude difference of the casing wave at different azimuths becomes more evident.The test results showed that the AABT can accurately locate the casing collars and evaluate the cement bond quality with azimuthal resolution at the casing-cement interface,and can visualize the size,depth,and azimuth of channeling.In the case of good casingcement bonding,the AABT can further evaluate the cement bond quality at the cementformation interface with azimuthal resolution by using the amplitude map and the velocity of the formation wave.展开更多
This study focuses on the factors that may affect the feasibility of performing elliptical anisotropy analysis on azimuthal PP and PS-wave data in HTI media, with the aim of using the modeling results as guidance in r...This study focuses on the factors that may affect the feasibility of performing elliptical anisotropy analysis on azimuthal PP and PS-wave data in HTI media, with the aim of using the modeling results as guidance in real seismic data application. Our results reveal that there is an offset limitation for both PP- and PS-waves in elliptical anisotropy fitting, and that PS-waves show a wider applicable offset range and larger observable azimuthal anisotropy than PP-waves. The major axis of the elliptical fit to the amplitudes of the R-component is perpendicular to the fracture strike, which is opposite to that in PP-wave analysis. The azimuthal interval travel time of PS-waves shows a nearly elliptical distribution and the major axis of the fit ellipse is perpendicular to the fracture strike, which is same as that in PP-wave analysis. For data within the applicable offset range, the anisotropic magnitude obtained from amplitude and travel time attributes of PP- and PS-waves exhibits a dependence on fracture density, and the major to minor axis ratio of the fit ellipse may be used to infer the relative distribution of fracture densities.展开更多
According to the characteristics of gravity passive navigation, this paper presents a novel gravity passive navigation system （GPNS）, which consists of the rate azimuth platform （RAP）, gravity sensor, digitally st...According to the characteristics of gravity passive navigation, this paper presents a novel gravity passive navigation system （GPNS）, which consists of the rate azimuth platform （RAP）, gravity sensor, digitally stored gravity maps, depth sensor and relative log. The algorithm of rate azimuth platform inertial navigation system, error state-space equations, measurement equations and GPNS optimal filter are described. In view of the measurements made by an onboard gravity sensor the Eotvos effect is introduced in the gravity measurement equation of a GPNS optimal filter. A GPNS is studied with the Matlab/Simulink tools; simulation results demonstrate that a GPNS has small errors in platform attitude and position. Because the inertial navigation platform is the rate azimuth platform in the GPNS and gravity sensor is mounted on the rate azimuth platform, the cost of the GPNS is lower than existing GPNS＇s and according to the above results the GPNS meets the need to maintain accuracy navigation for underwater vehicles over long intervals.展开更多
The cost of the gravity passive inertial navigation system will be lower witha rate azimuth platform and gravity sensor constituting a gravity measurement and navigationsystem. According to the system performance char...The cost of the gravity passive inertial navigation system will be lower witha rate azimuth platform and gravity sensor constituting a gravity measurement and navigationsystem. According to the system performance characteristics, we study the rate azimuth platforminertial navigation system (RAPINS), give the system navigation algorithm, error equations of theattitude, velocity and position of the rate azimuth platform, and random error models of theaccelerometer and gyro. Using the MATLAB/Simulink tools, we study the RAPINS and RAPINS withvelocity damping. Simulation results demonstrate that the RAPINS with velocity damping has smallerrors in platform attitude and position and satisfies gravity measurement and navigationrequirement.展开更多
Because the vector of angular momentum of ESG (electrostatically suspended gyroscope) maintains the fixed direction in inertial space, it may be regarded as a fixed star. The astronavigation aigorithm is used to estim...Because the vector of angular momentum of ESG (electrostatically suspended gyroscope) maintains the fixed direction in inertial space, it may be regarded as a fixed star. The astronavigation aigorithm is used to estimate the azimuth angle and the gyro constant drift in the paper. The relative errors which affect the estimation accuracy of the azimuth angle are the analysed.展开更多
The recently built China Digital Seismic Network consists of the China National Digital Seismic Network （CNDSN）, 31 regional seismic networks and several small aperture arrays with more than 1 000 stations including...The recently built China Digital Seismic Network consists of the China National Digital Seismic Network （CNDSN）, 31 regional seismic networks and several small aperture arrays with more than 1 000 stations including 850＋ broadband stations. It forms a gigantic seismic array that provides an unprecedented opportunity to study the Earth＇s deep interior besides its routine task of seismic monitoring. Many modern seismic studies rely on rotation of vertical and horizontal components in order to separate different types of seismic waves. Knowledge of the orientations of the two horizontal components thus is important to perform a correction rotation. We analyzed particle motions of teleseismic P waves recorded by the network and used them to estimate the northcomponent azimuth of each station. An SNR-weighted-multi-event method was introduced to obtain component azimuths that best explain the P-wave particle motions of all the events recorded at a station. The method provides robust estimates including a measurement error calculated from background noise levels. We found that about one third of the stations have some sort of problems, including misorientation of the two horizontal components, mislabeling and polarity reversal in one or more components. These problems need to be taken into account for any rotation based seismic studies.展开更多
Quantitative inversion of fracture weakness plays an important role in fracture prediction.Considering reservoirs with a set of vertical fractures as horizontal transversely isotropic media,the logarithmic normalized ...Quantitative inversion of fracture weakness plays an important role in fracture prediction.Considering reservoirs with a set of vertical fractures as horizontal transversely isotropic media,the logarithmic normalized azimuthal elastic impedance(EI)is rewritten in terms of Fourier coefficients(FCs),the 90°ambiguity in the azimuth estimation of the symmetry axis is resolved by judging the sign of the second FC,and we choose the FCs with the highest sensitivity to fracture weakness and present a feasible inversion workflow for fracture weakness,which involves:(1)the inversion for azimuthal EI datasets from observed azimuthal angle gathers;(2)the prediction for the second FCs and azimuth of the symmetry axis from the estimated azimuthal EI datasets;and(3)the estimation of fracture weakness combining the extracted second FCs and azimuth of the symmetry axis iteratively,which is constrained utilizing the Cauchy sparse regularization and the low-frequency regularization in a Bayesian framework.Tests on synthetic and field data demonstrate that the 90°ambiguity in the azimuth estimation of the symmetry axis has been removed,and reliable fracture weakness can be obtained when the estimated azimuth of the symmetry axis deviates less than 30°,which can guide the prediction of fractured reservoirs.展开更多
A type of specific fractured hydrocarbon reservoir, a shale fractured reservoir, exists in the Shengli oilfield. Due to very small porosity of this type, low sensitivity to the variation of petrochemical property para...A type of specific fractured hydrocarbon reservoir, a shale fractured reservoir, exists in the Shengli oilfield. Due to very small porosity of this type, low sensitivity to the variation of petrochemical property parameters, and strong anisotropy, it is very difficult to explore for them. So far, there is no set of mature methods for recognition of direction, distribution, and density of the fractures by an integrated analysis of geologic, geophysical, well log, drilling data, and etc. This paper presents a new method for acoustic impedance variation with azimuth （IPVA）, based on existing fracture detection methods. Seismic acquisition, processing, and recognition techniques were developed for detecting directional vertical fractures using multi-azimuth P wave data in combination with the seismic and geological features of shale fractures in the Luojia area. The IPVA research is carried out for recognizing the distribution, strike, and density of fractures based on the study of velocity variation with azimuth （VVA） and amplitude variation with azimuth （AVA） for full azimuth P wave data at different CMP positions. Through practical application in the Luojia area, primary results have been obtained which verifies that the IPVA method provides good potential for quantitative detection of parallel, high angle, shale fractures.展开更多
Over 800 accelerograms recorded by 272 ground-level stations during the Wenchuan earthquake are used to analyze the influence of rupture distance, local site conditions and azimuth on peak ground accelerations （PGAs...Over 800 accelerograms recorded by 272 ground-level stations during the Wenchuan earthquake are used to analyze the influence of rupture distance, local site conditions and azimuth on peak ground accelerations （PGAs）. To achieve a better understanding of the characteristics of ground motions, the spatial distributions of the EW, NS and UD components of PGAs are obtained. Comparisons between the EW and NS components, the fault-normal and fault-parallel components, and the vertical and horizontal components of PGAs are performed, and the regression formula of the vertical-to-horizontal ratio of PGAs is developed. The attenuation relationship of peak horizontal accelerations （PHAs） is compared with several contemporary attenuation relationships. In addition, an analysis of residuals is conducted to identify the potential effects of rupture distance, azimuth and site conditions on the observed values of PHAs. The analysis focuses on medium-hard soil site conditions, as they provided most of the data used in this study.展开更多
基金financially supported by the National Key Research and Development Program of China (2022YFC3005600)the Foundation of the Anhui Educational Commission (2023AH051198)+1 种基金the National Natural Science Foundation of China (42125401 and 42104063)the Joint Open Fund of Mengcheng National Geophysical Observatory (MENGO-202201)。
文摘The Tan-Lu Fault Zone is a large NNE-trending fault zone that has a substantial effect on the development of eastern China and its earthquake disaster prevention efforts. Aiming at the azimuthally anisotropic structure in the upper crust and seismogenic tectonics in the Hefei segment of this fault, we collected phase velocity dispersion data of fundamental mode Rayleigh waves from ambient noise cross-correlation functions of ~400 temporal seismographs in an area of approximately 80 × 70 km along the fault zone. The period band of the dispersion data was ~0.5–10 s. We inverted for the upper crustal three-dimensional(3-D) shear velocity model with azimuthal anisotropy from the surface to 10 km depth by using a 3-D direct azimuthal anisotropy inversion method. The inversion result shows the spatial distribution characteristics of the tectonic units in the upper crust. Additionally, the deformation of the Tan-Lu Fault Zone and its conjugated fault systems could be inferred from the anisotropy model. In particular, the faults that have remained active from the early and middle Pleistocene control the anisotropic characteristics of the upper crustal structure in this area. The direction of fast axes near the fault zone area in the upper crust is consistent with the strike of the faults, whereas for the region far away from the fault zone, the direction of fast axes is consistent with the direction of the regional principal stress caused by plate movement. Combined with the azimuthal anisotropy models in the deep crust and uppermost mantle from the surface wave and Pn wave, the different anisotropic patterns caused by the Tan-Lu Fault Zone and its conjugated fault system nearby are shown in the upper and lower crust. Furthermore,by using the double-difference method, we relocated the Lujiang earthquake series, which contained 32 earthquakes with a depth shallower than 10 km. Both the Vs model and earthquake relocation results indicate that earthquakes mostly occurred in the vicinity of structural boundaries with fractured media, with high-level development of cracks and small-scale faults jammed between more rigid areas.
基金supported by the National Key R&D Program of China(Grant No.2020YFA0710604)NSFC(Grant No.42374064).
文摘Determining the orientation of in-situ stresses is crucial for various geoscience and engineering appli-cations.Conventional methods for estimating these stress orientations often depend on focal mechanism solutions(FMSs)derived from earthquake data and formation micro-imager(FMI)data from well logs.However,these techniques can be costly,depth-inaccurate,and may lack spatial coverage.To address this issue,we introduce the use of three-dimensional(3D)seismic data(active sources)as a lateral constraint to approximate the 3D stress orientation field.Recognizing that both stress and fracture patterns are closely related to seismic velocity anisotropy,we derive the orientation of azimuthal anisotropy from multi-azimuth 3D seismic data to compensate for the lack of spatial stress orientation information.We apply our proposed workflow to a case study in the Weiyuan area of the Sichuan Basin,China,a region targeted for shale gas production.By integrating diverse datasets,including 3D seismic,earthquakes,and well logs,we develop a comprehensive 3D model of in-situ stress(orientations and magnitudes).Our results demonstrate that the estimated anisotropy orientations from 3D seismic data are consistent with the direction of maximum horizontal principal stress(SHmax)obtained from FMIs.We analyzed 12 earthquakes(magnitude>3)recorded between 2016 and 2020 for their FMSs and compressional axis(P-axis)orientations.The derived SHmax direction from our 3D stress model is 110°ES(East-South),which shows excellent agreement with the FMSs(within 3.96°).This close alignment validates the reliability and precision of our integrated method for predicting 3D SHmax orientations.
基金supported by the National Natural Science Foundation of China(Grant Nos.42174129 and 41374143)the project of"Vice President of Science and Technology"of Changping District,Beijing(Grant No.202502003019)。
文摘The characterization of subsurface formations via the analysis of near-wellbore velocity profiles represents a crucial method in geophysical exploration.This technique enables the evaluation of key parameters,including rock brittleness,wellbore stability,fracturing effects,and invasion extent,thereby enhancing comprehension of formation structures and informing exploration strategies.However,traditional near-wellbore formation velocity imaging methods exhibit two principal limitations.First,these methods lack azimuthal sensitivity,yielding results averaged across all directions.Second,they are computationally intensive and impractical for well-site environments.To address these drawbacks,we developed a rapid 3D velocity imaging method for array acoustic logging instru ments equipped with azimuthal receivers,capable of producing 3D imaging results efficiently.The workflow entails the following steps:(1)Band-pass filtering of logging data to mitigate scattered wave interference caused by formation heterogeneity near the wellbore;(2)combination of receivers with varying detection ranges in each direction to derive radial velocity sequences,followed by integration of ray-tracing theory to obtain 2D velocity distributions;and(3)synthesis of final 3D velocity imaging results via interpolation of these 2D datasets.In the velocity sequence extraction process,we significantly reduced the computational load by employing an adaptive time window,ensuring rapid and stable application in well-site settings.We utilized the finite difference method to construct well models with heteroge neous formations.The compressional and shear wave 3D velocity imaging results derived from synthetic data correlated with the model,demonstrating the azimuthal sensitivity of our proposed method.Furthermore,we applied this method to a well in West China,successfully identifying the azimuth of nearwellbore anisotropy.
基金sponsorship of the National Natural Science Foundation of China(42430809,42274157,42030103,42404132)the Fund of State Key Laboratory of Deep Oil and Gas,China University of Petroleum(East China)(SKLDOG2024-ZYTS-02)+5 种基金the Postdoctoral Fellowship Program of CPSF(GZB20240850)the Postdoctoral Project of Qingdao(QDBSH20240102082)the Fundamental Research Funds for the Central Universities(24CX07004A,24CX06036A)the CNPC Innovation Fund(2024DQ02-0505,2024DQ02-0136)the Innovation fund project for graduate student of China University of Petroleum(East China)the Fundamental Research Funds for the Central Universities(24CX04002A).
文摘Reservoirs with a group of vertical fractures in a vertical transversely isotropic(VTI)background are considered as orthorhombic(ORT)medium.However,fracture detection in ORT medium using seismic inversion methods remains challenging,as it requires the estimation of more than eight parameters.Assuming the reservoir to be a weakly anisotropic ORT medium with small contrasts in the background elastic parameters,a new azimuthal elastic impedance equation was first derived using parameter combinations and mathematical approximations.This equation exhibited almost the same accuracy as the original equation and contained only six model parameters:the compression modulus,anisotropic shear modulus,anisotropic compression modulus,density,normal fracture weakness,and tangential fracture weakness.Subsequently,a stepwise inversion method using second-order derivatives of the elastic impedance was developed to estimate these parameters.Moreover,the Thomsen anisotropy parameter,epsilon,was estimated from the inversion results using the ratio of the anisotropic compression modulus to the compression modulus.Synthetic examples with moderate noise and field data examples confirm the feasibility and effectiveness of the inversion method.The proposed method exhibited accuracy similar to that of previous inversion strategies and could predict richer vertical fracture information.Ultimately,the method was applied to a three-dimensional work area,and the predictions were consistent with logging and geological a priori information,confirming the effectiveness of this method.Summarily,the proposed stepwise inversion method can alleviate the uncertainty of multi-parameter inversion in ORT medium,thereby improving the reliability of fracture detection.
基金supported by the National Natural Science Foundation of China under Grant 62301051.
文摘Range-azimuth imaging of ground targets via frequency-modulated continuous wave(FMCW)radar is crucial for effective target detection.However,when the pitch of the moving array constructed during motion exceeds the physical array aperture,azimuth ambiguity occurs,making range-azimuth imaging on a moving platform challenging.To address this issue,we theoretically analyze azimuth ambiguity generation in sparse motion arrays and propose a dual-aperture adaptive processing(DAAP)method for suppressing azimuth ambiguity.This method combines spatial multiple-input multiple-output(MIMO)arrays with sparse motion arrays to achieve high-resolution range-azimuth imaging.In addition,an adaptive QR decomposition denoising method for sparse array signals based on iterative low-rank matrix approximation(LRMA)and regularized QR is proposed to preprocess sparse motion array signals.Simulations and experiments show that on a two-transmitter-four-receiver array,the signal-to-noise ratio(SNR)of the sparse motion array signal after noise suppression via adaptive QR decomposition can exceed 0 dB,and the azimuth ambiguity signal ratio(AASR)can be reduced to below-20 dB.
基金supported by the National Natural Science Foundation of China under Grants No.41771401 and No.42350710201.
文摘A P-band polarimetric synthetic aperture radar(PolSAR)sensor has deep penetration ability into and through the vegetation canopies in forested environments.Thus,the sensor is of great potential to accurately assess forest parameters such as coverage,stand density,and tree height.Unfortunately,the radar backscatter from complex terrain can adversely impact the backscatter from trees or forests,and forest parameters assessed can be erroneous.Thus,reducing the topographic impact is an urgent must.In this study,a topographic compensation algorithm has been studied.To assess the algorithm’s validity and effectiveness,we applied it to P-band PolSAR datasets in four forested areas in the US.Trees in the forest stands have diverse species,and the topographic conditions of the terrain differ.Significant topographic impact on the P-band PolSAR data exists before the topographic compensation algorithm.After the algorithm,the impact decreases noticeably qualitatively and quantitatively.The algorithm is valid and effective in reducing the topographic influence on the PolSAR data and,consequently,provides a better chance of retrieving accurate forest parameters.
基金supported by the National Natural Science Foundation of China(Nos.42076238,42376174)the Natural Science Foundation of Shanghai(No.23ZR1426900).
文摘The purpose of this study is to analyze the impact of upper ocean dynamics on velocity bunching,represented by az-imuthal cutoff wavelength(i.e.,sea surface wind,wave,and current).In this study,over 1400 dual-polarized(vertical-vertical(VV)and vertical-horizontal(VH))Sentinel-1(S-1)synthetic aperture radar(SAR)images collected in tropical cyclones(TC)are utilized.These images are combined with wind and rain observations from the stepped-frequency microwave radiometer(SFMR),wave simu-lations conducted using a third-generation numerical wave model,WAVEWATCH-III(WW3),and SAR-derived wind information collected from CyclObs winds.The WW3-simulated significant wave height(SWH)is validated against measurements from HY-2B altimeter taken in August and September 2021,yielding a root mean square error(RMSE)of 0.48 m and a correlation coefficient(COR)of 0.88.The SAR-based azimuthal cutoff wavelengths in VV polarization,which quantitatively represent the effect of velocity bunching,are compared with theoretical values calculated using WW3-simulated SWH.A notable relationship is observed between the difference in azimuthal cutoff wavelength and SAR-derived wind speed and WW3-simulated SWH.Analysis results show that the correlation between SAR-based azimuthal cutoff wavelength and SWH is stronger than that with wind and current.Finally,a machine learning algorithm is used to develop an algorithm aimed at simulating the azimuthal cutoff wavelength in TCs,including wind,wave,and incidence angle.This method yields an RMSE of 8.90 m,a COR of 0.91,and a scatter index of 0.04 for VV-polar-ization SAR.
基金supported by the Government of the Russian Federation grant number 075-15-2025-009 of 28 February 2025 and by the Russian Science Foundation,Grant No.24-72-10127.
文摘During daylight laser polarization sensing of high-level clouds(HLCs),the lidar receiving system generates a signal caused by not only backscattered laser radiation,but also scattered solar radiation,the intensity and polarization of which depends on the Sun’s location.If a cloud contains spatially oriented ice particles,then it becomes anisotropic,that is,the coefficients of directional light scattering of such a cloud depend on the Sun’s zenith and azimuth angles.In this work,the possibility of using the effect of anisotropic scattering of solar radiation on the predictive ability of machine learning algorithms in solving the problem of predicting the HLC backscattering phase matrix(BSPM)was evaluated.The hypothesis that solar radiation scattered on HLCs has no effect on the BSPM elements of such clouds determined with a polarization lidar was tested.The operation of two algorithms for predicting the BSPM elements is evaluated.To train the first one,meteorological data were used as input parameters;for the second algorithm,the azi-muthal and zenith angles of the Sun’s position were added to the meteorological parameters.It is shown that there is no significant improvement in the predictive ability of the algorithm.
基金sponsored by the National Natural Science Foundation of China(under Grant Nos.41404090 and U1262208
文摘We propose a method for mOdeling azimuthal AVO responses from a fractured i reflector. The method calculates the integrated reflected wavetrains, and the wavetrains contain elastodynamic information including the contrast in impedance and anelasticity i across interfaces, the intemal anisotropic propagation, the dispersion and attenuation along i the wave path, and tuning and interference. The results suggest that for large angles of incidence, the velocity dispersion and attenuation increase the amplitudes of PP waves from the top and decrease those from the bottom. For azimuthal responses at specific angles of incidence, the reflected wavetrains of PP waves tend to have longer duration with increasing azimuth. In contrast, model-converted PSV and PSH reflections show stable azimuthal features and are less affected by the reflector thickness. The amplitudes of PSV reflections increase with increasing azimuth; moreover, the waves have no reflection energy at 0° and 90° azimuth and maximum amplitude at 45° azimuth.
基金supported by the National Basic Research Program (the 973 Program, No. 2007CB209601).
文摘The effect of the wide and narrow azimuth 3D observation systems on seismic imaging precision is becoming a hot area for studies of high precision 3D seismic acquisition methods in recent years. In this paper we utilize 3D physical seismic modeling experiments. A 3D channel sand body physical seismic model is constructed and two acquisition systems are designed with wide azimuth （16 lines） and narrow azimuth （8 lines） to model 3D seismic data acquisition and processing seismic work flows. From analysis of migrated time slice data with high quality and small size, we conclude that when the overlying layers are smooth and lateral velocities have little change, both wide and narrow azimuth observation systems in 3D acquisition can be used for obtaining high precision imaging and equivalent resolution of the channel sand body.
基金supported by the National Natural Science Foundation of China(Nos.11204380,11374371,61102102,and11134011)National Science and Technology Major Project(No.2011ZX05020-009)+1 种基金China National Petroleum Corporation(Nos.2014B-4011,2014D-4105,and 2014A-3912)PetroChina Innovation Foundation(No.2014D-5006-0307)
文摘We developed a novel cement evaluation logging tool,named the azimuthally acoustic bond tool（AABT）,which uses a phased-arc array transmitter with azimuthal detection capability.We combined numerical simulations and field tests to verify the AABT tool.The numerical simulation results showed that the radiation direction of the subarray corresponding to the maximum amplitude of the first arrival matches the azimuth of the channeling when it is behind the casing.With larger channeling size in the circumferential direction,the amplitude difference of the casing wave at different azimuths becomes more evident.The test results showed that the AABT can accurately locate the casing collars and evaluate the cement bond quality with azimuthal resolution at the casing-cement interface,and can visualize the size,depth,and azimuth of channeling.In the case of good casingcement bonding,the AABT can further evaluate the cement bond quality at the cementformation interface with azimuthal resolution by using the amplitude map and the velocity of the formation wave.
文摘This study focuses on the factors that may affect the feasibility of performing elliptical anisotropy analysis on azimuthal PP and PS-wave data in HTI media, with the aim of using the modeling results as guidance in real seismic data application. Our results reveal that there is an offset limitation for both PP- and PS-waves in elliptical anisotropy fitting, and that PS-waves show a wider applicable offset range and larger observable azimuthal anisotropy than PP-waves. The major axis of the elliptical fit to the amplitudes of the R-component is perpendicular to the fracture strike, which is opposite to that in PP-wave analysis. The azimuthal interval travel time of PS-waves shows a nearly elliptical distribution and the major axis of the fit ellipse is perpendicular to the fracture strike, which is same as that in PP-wave analysis. For data within the applicable offset range, the anisotropic magnitude obtained from amplitude and travel time attributes of PP- and PS-waves exhibits a dependence on fracture density, and the major to minor axis ratio of the fit ellipse may be used to infer the relative distribution of fracture densities.
文摘According to the characteristics of gravity passive navigation, this paper presents a novel gravity passive navigation system （GPNS）, which consists of the rate azimuth platform （RAP）, gravity sensor, digitally stored gravity maps, depth sensor and relative log. The algorithm of rate azimuth platform inertial navigation system, error state-space equations, measurement equations and GPNS optimal filter are described. In view of the measurements made by an onboard gravity sensor the Eotvos effect is introduced in the gravity measurement equation of a GPNS optimal filter. A GPNS is studied with the Matlab/Simulink tools; simulation results demonstrate that a GPNS has small errors in platform attitude and position. Because the inertial navigation platform is the rate azimuth platform in the GPNS and gravity sensor is mounted on the rate azimuth platform, the cost of the GPNS is lower than existing GPNS＇s and according to the above results the GPNS meets the need to maintain accuracy navigation for underwater vehicles over long intervals.
文摘The cost of the gravity passive inertial navigation system will be lower witha rate azimuth platform and gravity sensor constituting a gravity measurement and navigationsystem. According to the system performance characteristics, we study the rate azimuth platforminertial navigation system (RAPINS), give the system navigation algorithm, error equations of theattitude, velocity and position of the rate azimuth platform, and random error models of theaccelerometer and gyro. Using the MATLAB/Simulink tools, we study the RAPINS and RAPINS withvelocity damping. Simulation results demonstrate that the RAPINS with velocity damping has smallerrors in platform attitude and position and satisfies gravity measurement and navigationrequirement.
文摘Because the vector of angular momentum of ESG (electrostatically suspended gyroscope) maintains the fixed direction in inertial space, it may be regarded as a fixed star. The astronavigation aigorithm is used to estimate the azimuth angle and the gyro constant drift in the paper. The relative errors which affect the estimation accuracy of the azimuth angle are the analysed.
基金supported by NSF grant EAR-063566(F.N.)National Natural Science Foundation of China grant 40774042(J.L.)
文摘The recently built China Digital Seismic Network consists of the China National Digital Seismic Network （CNDSN）, 31 regional seismic networks and several small aperture arrays with more than 1 000 stations including 850＋ broadband stations. It forms a gigantic seismic array that provides an unprecedented opportunity to study the Earth＇s deep interior besides its routine task of seismic monitoring. Many modern seismic studies rely on rotation of vertical and horizontal components in order to separate different types of seismic waves. Knowledge of the orientations of the two horizontal components thus is important to perform a correction rotation. We analyzed particle motions of teleseismic P waves recorded by the network and used them to estimate the northcomponent azimuth of each station. An SNR-weighted-multi-event method was introduced to obtain component azimuths that best explain the P-wave particle motions of all the events recorded at a station. The method provides robust estimates including a measurement error calculated from background noise levels. We found that about one third of the stations have some sort of problems, including misorientation of the two horizontal components, mislabeling and polarity reversal in one or more components. These problems need to be taken into account for any rotation based seismic studies.
基金the sponsorship of the National Natural Science Foundation of China(41674130)National Grand Project for Science and Technology(2016ZX05002-005)for funding this research.
文摘Quantitative inversion of fracture weakness plays an important role in fracture prediction.Considering reservoirs with a set of vertical fractures as horizontal transversely isotropic media,the logarithmic normalized azimuthal elastic impedance(EI)is rewritten in terms of Fourier coefficients(FCs),the 90°ambiguity in the azimuth estimation of the symmetry axis is resolved by judging the sign of the second FC,and we choose the FCs with the highest sensitivity to fracture weakness and present a feasible inversion workflow for fracture weakness,which involves:(1)the inversion for azimuthal EI datasets from observed azimuthal angle gathers;(2)the prediction for the second FCs and azimuth of the symmetry axis from the estimated azimuthal EI datasets;and(3)the estimation of fracture weakness combining the extracted second FCs and azimuth of the symmetry axis iteratively,which is constrained utilizing the Cauchy sparse regularization and the low-frequency regularization in a Bayesian framework.Tests on synthetic and field data demonstrate that the 90°ambiguity in the azimuth estimation of the symmetry axis has been removed,and reliable fracture weakness can be obtained when the estimated azimuth of the symmetry axis deviates less than 30°,which can guide the prediction of fractured reservoirs.
文摘A type of specific fractured hydrocarbon reservoir, a shale fractured reservoir, exists in the Shengli oilfield. Due to very small porosity of this type, low sensitivity to the variation of petrochemical property parameters, and strong anisotropy, it is very difficult to explore for them. So far, there is no set of mature methods for recognition of direction, distribution, and density of the fractures by an integrated analysis of geologic, geophysical, well log, drilling data, and etc. This paper presents a new method for acoustic impedance variation with azimuth （IPVA）, based on existing fracture detection methods. Seismic acquisition, processing, and recognition techniques were developed for detecting directional vertical fractures using multi-azimuth P wave data in combination with the seismic and geological features of shale fractures in the Luojia area. The IPVA research is carried out for recognizing the distribution, strike, and density of fractures based on the study of velocity variation with azimuth （VVA） and amplitude variation with azimuth （AVA） for full azimuth P wave data at different CMP positions. Through practical application in the Luojia area, primary results have been obtained which verifies that the IPVA method provides good potential for quantitative detection of parallel, high angle, shale fractures.
基金National Natural Science Foundation of China Under Grant No. 90715038, 50878199 and 50808166National Basic Research Program of China Under Grant No. 2007CB714200
文摘Over 800 accelerograms recorded by 272 ground-level stations during the Wenchuan earthquake are used to analyze the influence of rupture distance, local site conditions and azimuth on peak ground accelerations （PGAs）. To achieve a better understanding of the characteristics of ground motions, the spatial distributions of the EW, NS and UD components of PGAs are obtained. Comparisons between the EW and NS components, the fault-normal and fault-parallel components, and the vertical and horizontal components of PGAs are performed, and the regression formula of the vertical-to-horizontal ratio of PGAs is developed. The attenuation relationship of peak horizontal accelerations （PHAs） is compared with several contemporary attenuation relationships. In addition, an analysis of residuals is conducted to identify the potential effects of rupture distance, azimuth and site conditions on the observed values of PHAs. The analysis focuses on medium-hard soil site conditions, as they provided most of the data used in this study.
