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.展开更多
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.展开更多
The elliptic azimuthal anisotropy coefficient(v_(2))of the identified particles at midrapidity(|η|<0.8)was investigated in p-Pb collisions at√s_(NN)=5.02 TeV using a multi-phase transport model(AMPT).The calculat...The elliptic azimuthal anisotropy coefficient(v_(2))of the identified particles at midrapidity(|η|<0.8)was investigated in p-Pb collisions at√s_(NN)=5.02 TeV using a multi-phase transport model(AMPT).The calculations of differential v_(2)based on the advanced flow extraction method of light flavor hadrons(pions,kaons,protons,andΛ)in small collision systems were extended to a wider transverse momentum(p_(T))range of up to 8 GeV/c for the first time.The string-melting version of the AMPT model provides a good description of the measured p_(T)-differential v_(2)of the mesons but exhibits a slight deviation from the baryon v_(2).In addition,we observed the features of mass ordering at low p_(T)and the approximate number-of-constituentquark(NCQ)scaling at intermediate p_(T).Moreover,we demonstrate that hadronic rescattering does not have a significant impact on v_(2)in p-Pb collisions for different centrality selections,whereas partonic scattering dominates in generating the elliptic anisotropy of the final particles.This study provides further insight into the origin of collective-like behavior in small collision systems and has referential value for future measurements of azimuthal anisotropy.展开更多
Cement density monitoring plays a vital role in evaluating the quality of cementing projects,which is of great significance to the development of oil and gas.However,the presence of inhomogeneous cement distribution a...Cement density monitoring plays a vital role in evaluating the quality of cementing projects,which is of great significance to the development of oil and gas.However,the presence of inhomogeneous cement distribution and casing eccentricity in horizontal wells often complicates the accurate evaluation of cement azimuthal density.In this regard,this paper proposes an algorithm to calculate the cement azimuthal density in horizontal wells using a multi-detector gamma-ray detection system.The spatial dynamic response functions are simulated to obtain the influence of cement density on gamma-ray counts by the perturbation theory,and the contribution of cement density in six sectors to the gamma-ray recorded by different detectors is obtained by integrating the spatial dynamic response functions.Combined with the relationship between gamma-ray counts and cement density,a multi-parameter calculation equation system is established,and the regularized Newton iteration method is employed to invert casing eccentricity and cement azimuthal density.This approach ensures the stability of the inversion process while simultaneously achieving an accuracy of 0.05 g/cm^(3) for the cement azimuthal density.This accuracy level is ten times higher compared to density accuracy calculated using calibration equations.Overall,this algorithm enhances the accuracy of cement azimuthal density evaluation,provides valuable technical support for the monitoring of cement azimuthal density in the oil and gas industry.展开更多
Imaging detection is an important means to obtain target information.The traditional imaging detection technology mainly collects the intensity information and spectral information of the target to realize the classif...Imaging detection is an important means to obtain target information.The traditional imaging detection technology mainly collects the intensity information and spectral information of the target to realize the classification of the target.In practical applications,due to the mixed scenario,it is difficult to meet the needs of target recognition.Compared with intensity detection,the method of polarization detection can effectively enhance the accuracy of ground object target recognition(such as the camouflage target).In this paper,the reflection mechanism of the target surface is studied from the microscopic point of view,and the polarization characteristic model is established to express the relationship between the polarization state of the reflected signal and the target surface parameters.The polarization characteristic test experiment is carried out,and the target surface parameters are retrieved using the experimental data.The results show that the degree of polarization(DOP)is closely related to the detection zenith angle and azimuth angle.The(DOP)of the target is the smallest in the direction of light source incidence and the largest in the direction of specular reflection.Different materials have different polarization characteristics.By comparing their DOP,target classification can be achieved.展开更多
The article presents an experimental study on the flow of an eutectic gallium alloy in a cylindrical cell,which is placed in an alternating magnetic field.The magnetic field is generated by a coil connected to an alte...The article presents an experimental study on the flow of an eutectic gallium alloy in a cylindrical cell,which is placed in an alternating magnetic field.The magnetic field is generated by a coil connected to an alternating current source.The coil is located at a fixed height in such a way that its plane is perpendicular to the gravity vector,which in turn is parallel to the axis of the cylinder.The position of the cylinder can vary in height with respect to the coil.The forced flow of the considered electrically conductive liquid is generated due to the action of the localized electromagnetic force.It is assumed that under the action of the alternating magnetic field,the liquid is heated uniformly,and the resulting heat is quickly absorbed by the forced flow,so that liquid free convection can be neglected.The experiment is carried out using an ultrasonic Doppler anemometer.One transducer is installed in the axially located cylinder sluice and the other transducer is placed in the near-wall region.According to the results,a velocity profile,corresponding to a two-tori flow pattern can be hardly obtained in the low frequency range of the power supply.However,this is possible in the high frequency range.The average velocity profiles depend essentially on the location of the coil relative to the cell.The spectral analysis of velocity signals shows that the amplitude of the velocity pulsations is comparable to the average value of the flow velocity.Such experimental results and their verification through comparison with numerical calculations are intended to support the development of new methods for reducing the intensity of vortex flows during the electromagnetic separation of impurities through an electromagnetic induction mechanism(able to produce an electromotive force that displaces particles).展开更多
The Binchuan Basin in northwest Yunnan, southwest China, is a rift basin developed at the intersection of the Red River Fault and Chenghai Fault, where historical earthquakes have occurred. Understanding the fine velo...The Binchuan Basin in northwest Yunnan, southwest China, is a rift basin developed at the intersection of the Red River Fault and Chenghai Fault, where historical earthquakes have occurred. Understanding the fine velocity structure of the shallow crust in this region can help improve earthquake location accuracy and our understanding of the relationship between fault zone structures and fault slip behaviors. Using the continuous waveform data recorded by 381 dense array stations in 2017, we obtained 7 915 Rayleigh-wave phase velocity dispersion curves in the period band of 0.2–6 s from ambient noise cross-correlation functions after rigorous data processing and quality control. We determined 3D isotropic and azimuthally anisotropic shear wave velocity models at depths above 6 km in the shallow crust based on the direct surface wave azimuthal anisotropic tomography method. The isotropic model reveals a strong correspondence between the S-wave velocity structure at depths of 0–1 km and the regional topography and lithology. The Binchuan depocenter, Zhoucheng depocenter, Xiangyun Basin, and Xihai Rift Basin are primarily composed of Quaternary deposits, which show low-velocity anomalies, while the regions with the Paleozoic shale, limestone, and basalt exhibit high-velocity anomalies. The nearly N–S orientation of fast directions from azimuthal anisotropy models are mainly controlled by the active Binchuan Fault with N–S strike as well as the NNW-oriented primary compressive stress.展开更多
When inverting the S-wave velocity and azimuthal anisotropy from ambient noise data, it is always to obtain the partial overlapped inversion results in contiguous different regions. Merging different data to achieve a...When inverting the S-wave velocity and azimuthal anisotropy from ambient noise data, it is always to obtain the partial overlapped inversion results in contiguous different regions. Merging different data to achieve a consistent model becomes an essential requirement. Based on the S-wave velocity and azimuthal anisotropy obtained from different contiguous regions, this paper introduces three kinds of methods for merging data. For data from different regions with partial overlapping areas, the merged results could be calculated by direct average weighting(DAW), linear dynamic weighting(LDW), and Gaussian function weighting(GFW), respectively. Data tests demonstrate that the LDW and GFW methods can effectively merge data by reasonably allocating data weights to capitalize on the data quality advantages in each zone. In particular, they can resolve the data smoothness at the boundaries of data areas, resulting in a consistent data model in larger regions. This paper presents the effective methods and valuable experiences that can be referred to as advancing data merging technology.展开更多
The existing research results show that a fixed single station must conduct three consecutive frequency shift measurements and obtain the target’s moving speed by constructing two frequency difference equations. This...The existing research results show that a fixed single station must conduct three consecutive frequency shift measurements and obtain the target’s moving speed by constructing two frequency difference equations. This article proposes a new method that requires only two consecutive measurements. While using the azimuth measurement to obtain the angular difference between two radial distances, it also conducts two consecutive Doppler frequency shift measurements at the same target azimuth. On the basis of this measurement, a frequency difference equation is first constructed and solved jointly with the Doppler frequency shift equation. By eliminating the velocity variable and using the measured angular difference to obtain the target’s lead angle, the target’s velocity can be solved by using the Doppler frequency shift equation again. The new method avoids the condition that the target must move equidistantly, which not only provides an achievable method for engineering applications but also lays a good foundation for further exploring the use of steady-state signals to achieve passive positioning.展开更多
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.展开更多
基金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.
基金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.
基金This work was supported by the Key Laboratory of Quark and Lepton Physics(MOE)in Central China Normal University(Nos.QLPL2022P01,QLPL202106)Natural Science Foundation of Hubei Provincial Education Department(No.Q20131603)+2 种基金National key research,development program of China(No.2018YFE0104700)National Natural Science Foundation of China(No.12175085)Fundamental research funds for the Central Universities(No.CCNU220N003).
文摘The elliptic azimuthal anisotropy coefficient(v_(2))of the identified particles at midrapidity(|η|<0.8)was investigated in p-Pb collisions at√s_(NN)=5.02 TeV using a multi-phase transport model(AMPT).The calculations of differential v_(2)based on the advanced flow extraction method of light flavor hadrons(pions,kaons,protons,andΛ)in small collision systems were extended to a wider transverse momentum(p_(T))range of up to 8 GeV/c for the first time.The string-melting version of the AMPT model provides a good description of the measured p_(T)-differential v_(2)of the mesons but exhibits a slight deviation from the baryon v_(2).In addition,we observed the features of mass ordering at low p_(T)and the approximate number-of-constituentquark(NCQ)scaling at intermediate p_(T).Moreover,we demonstrate that hadronic rescattering does not have a significant impact on v_(2)in p-Pb collisions for different centrality selections,whereas partonic scattering dominates in generating the elliptic anisotropy of the final particles.This study provides further insight into the origin of collective-like behavior in small collision systems and has referential value for future measurements of azimuthal anisotropy.
基金The authors would like to acknowledge the support of the National Natural Science Foundation of China(41974127,42174147).References。
文摘Cement density monitoring plays a vital role in evaluating the quality of cementing projects,which is of great significance to the development of oil and gas.However,the presence of inhomogeneous cement distribution and casing eccentricity in horizontal wells often complicates the accurate evaluation of cement azimuthal density.In this regard,this paper proposes an algorithm to calculate the cement azimuthal density in horizontal wells using a multi-detector gamma-ray detection system.The spatial dynamic response functions are simulated to obtain the influence of cement density on gamma-ray counts by the perturbation theory,and the contribution of cement density in six sectors to the gamma-ray recorded by different detectors is obtained by integrating the spatial dynamic response functions.Combined with the relationship between gamma-ray counts and cement density,a multi-parameter calculation equation system is established,and the regularized Newton iteration method is employed to invert casing eccentricity and cement azimuthal density.This approach ensures the stability of the inversion process while simultaneously achieving an accuracy of 0.05 g/cm^(3) for the cement azimuthal density.This accuracy level is ten times higher compared to density accuracy calculated using calibration equations.Overall,this algorithm enhances the accuracy of cement azimuthal density evaluation,provides valuable technical support for the monitoring of cement azimuthal density in the oil and gas industry.
基金supported by the National Key Laboratory of Electromagnetic Space Security(JCKY2023230C009).
文摘Imaging detection is an important means to obtain target information.The traditional imaging detection technology mainly collects the intensity information and spectral information of the target to realize the classification of the target.In practical applications,due to the mixed scenario,it is difficult to meet the needs of target recognition.Compared with intensity detection,the method of polarization detection can effectively enhance the accuracy of ground object target recognition(such as the camouflage target).In this paper,the reflection mechanism of the target surface is studied from the microscopic point of view,and the polarization characteristic model is established to express the relationship between the polarization state of the reflected signal and the target surface parameters.The polarization characteristic test experiment is carried out,and the target surface parameters are retrieved using the experimental data.The results show that the degree of polarization(DOP)is closely related to the detection zenith angle and azimuth angle.The(DOP)of the target is the smallest in the direction of light source incidence and the largest in the direction of specular reflection.Different materials have different polarization characteristics.By comparing their DOP,target classification can be achieved.
基金supported by Russian Science Foundation Grant RSF-22-19-20106。
文摘The article presents an experimental study on the flow of an eutectic gallium alloy in a cylindrical cell,which is placed in an alternating magnetic field.The magnetic field is generated by a coil connected to an alternating current source.The coil is located at a fixed height in such a way that its plane is perpendicular to the gravity vector,which in turn is parallel to the axis of the cylinder.The position of the cylinder can vary in height with respect to the coil.The forced flow of the considered electrically conductive liquid is generated due to the action of the localized electromagnetic force.It is assumed that under the action of the alternating magnetic field,the liquid is heated uniformly,and the resulting heat is quickly absorbed by the forced flow,so that liquid free convection can be neglected.The experiment is carried out using an ultrasonic Doppler anemometer.One transducer is installed in the axially located cylinder sluice and the other transducer is placed in the near-wall region.According to the results,a velocity profile,corresponding to a two-tori flow pattern can be hardly obtained in the low frequency range of the power supply.However,this is possible in the high frequency range.The average velocity profiles depend essentially on the location of the coil relative to the cell.The spectral analysis of velocity signals shows that the amplitude of the velocity pulsations is comparable to the average value of the flow velocity.Such experimental results and their verification through comparison with numerical calculations are intended to support the development of new methods for reducing the intensity of vortex flows during the electromagnetic separation of impurities through an electromagnetic induction mechanism(able to produce an electromotive force that displaces particles).
基金supported by the National Key R&D Program of China (2021YFC3000704)the Special Fund of Key Laboratory of Earthquake Prediction,CEA (2021IEF0103)the National Natural Science Foundation of China (42125401)。
文摘The Binchuan Basin in northwest Yunnan, southwest China, is a rift basin developed at the intersection of the Red River Fault and Chenghai Fault, where historical earthquakes have occurred. Understanding the fine velocity structure of the shallow crust in this region can help improve earthquake location accuracy and our understanding of the relationship between fault zone structures and fault slip behaviors. Using the continuous waveform data recorded by 381 dense array stations in 2017, we obtained 7 915 Rayleigh-wave phase velocity dispersion curves in the period band of 0.2–6 s from ambient noise cross-correlation functions after rigorous data processing and quality control. We determined 3D isotropic and azimuthally anisotropic shear wave velocity models at depths above 6 km in the shallow crust based on the direct surface wave azimuthal anisotropic tomography method. The isotropic model reveals a strong correspondence between the S-wave velocity structure at depths of 0–1 km and the regional topography and lithology. The Binchuan depocenter, Zhoucheng depocenter, Xiangyun Basin, and Xihai Rift Basin are primarily composed of Quaternary deposits, which show low-velocity anomalies, while the regions with the Paleozoic shale, limestone, and basalt exhibit high-velocity anomalies. The nearly N–S orientation of fast directions from azimuthal anisotropy models are mainly controlled by the active Binchuan Fault with N–S strike as well as the NNW-oriented primary compressive stress.
基金supported by the National Natural Science Foundation of China (Project 42330311)the Central Public-interest Scientific Institution Basal Research Fund (No. 2021IEF0103)the National Key R&D Project of China (2017YFC1500304)。
文摘When inverting the S-wave velocity and azimuthal anisotropy from ambient noise data, it is always to obtain the partial overlapped inversion results in contiguous different regions. Merging different data to achieve a consistent model becomes an essential requirement. Based on the S-wave velocity and azimuthal anisotropy obtained from different contiguous regions, this paper introduces three kinds of methods for merging data. For data from different regions with partial overlapping areas, the merged results could be calculated by direct average weighting(DAW), linear dynamic weighting(LDW), and Gaussian function weighting(GFW), respectively. Data tests demonstrate that the LDW and GFW methods can effectively merge data by reasonably allocating data weights to capitalize on the data quality advantages in each zone. In particular, they can resolve the data smoothness at the boundaries of data areas, resulting in a consistent data model in larger regions. This paper presents the effective methods and valuable experiences that can be referred to as advancing data merging technology.
文摘The existing research results show that a fixed single station must conduct three consecutive frequency shift measurements and obtain the target’s moving speed by constructing two frequency difference equations. This article proposes a new method that requires only two consecutive measurements. While using the azimuth measurement to obtain the angular difference between two radial distances, it also conducts two consecutive Doppler frequency shift measurements at the same target azimuth. On the basis of this measurement, a frequency difference equation is first constructed and solved jointly with the Doppler frequency shift equation. By eliminating the velocity variable and using the measured angular difference to obtain the target’s lead angle, the target’s velocity can be solved by using the Doppler frequency shift equation again. The new method avoids the condition that the target must move equidistantly, which not only provides an achievable method for engineering applications but also lays a good foundation for further exploring the use of steady-state signals to achieve passive positioning.
基金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.
