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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
Fractured reservoirs always show anisotropic amplitude features,i.e.the reflection amplitude of seismic waves varies with offset and azimuth （AVOZ）.A noise attenuation fracture inversion algorithm is presented for f...Fractured reservoirs always show anisotropic amplitude features,i.e.the reflection amplitude of seismic waves varies with offset and azimuth （AVOZ）.A noise attenuation fracture inversion algorithm is presented for fracture detection based on P-wave AVOZ.The conventional inversion method always fails when applied to limited azimuth data because of the existence of noise.In our inversion algorithm,special attention is paid to suppressing the noise during inversion,to overcome the limitation of the conventional inversion method on limited azimuth data.Numerical models are employed to illustrate the effectiveness of the method.The inversion algorithm is then applied to Tazhong 45 area field data which is acquired under limited azimuth distribution.Compared with cores and fullbore formation microimage （FMI）,the inverted results （fracture density and orientation） are reasonable,suggesting that the inversion algorithm is feasible for fracture prediction in the Tarim Basin.展开更多
We present systematic investigations on the physics,detection performance and inversion of logging-while-drilling extradeep azimuthal resistivity measurements(EDARM).First,the definitions of EDRAM measurements are dis...We present systematic investigations on the physics,detection performance and inversion of logging-while-drilling extradeep azimuthal resistivity measurements(EDARM).First,the definitions of EDRAM measurements are discussed,followed by the derivation of the attenuation and phase-shift geometrical factors to illustrate the relative contributions of formation units to the observed signals.Then,a new definition of detection depth,which considers the uncertainty of inversion results caused by the data noise,is proposed to quantify the detection capability of ED ARM.Finally,the B ayesian theory associated with Markov chain Monte Carlo sampling is introduced for fast processing of EDARM data.Numerical results show that ED ARM is capable of detecting the azimuth and distance of remote bed boundaries,and the detection capability increases with increasing spacing and resistivity contrast.The EDARM tool can accommodate a large range of formation resistivity and is able to provide the resistivity anisotropy at arbitrary relative dipping angles.In addition,multiple bed boundaries and reservoir images near the borehole are readily obtained by using the Bayesian inversion.展开更多
In this article, a method of high order, staggered grid, finite difference was used to model snapshots of seismic wave propagation and multi-azimuth surface seismic records in viscoelastic cracked monoclinic media. Th...In this article, a method of high order, staggered grid, finite difference was used to model snapshots of seismic wave propagation and multi-azimuth surface seismic records in viscoelastic cracked monoclinic media. The modeling results clearly show the anisotropy caused by the crack fill material and the viscoelastic wave field through the monoclinic media. Summarizing the regularity of seismic records for different azimuths, this study is a significant attempt to understand the wave fields propagated in real media and to propose a more rational theoretical model for the future.展开更多
Seismic anisotropy provides important constraints on deformation patterns of Earth＇s material. Rayleigh wave dispersion data with azimuthal anisotropy can be used to invert for depth-dependent shear wavespeed azimuth...Seismic anisotropy provides important constraints on deformation patterns of Earth＇s material. Rayleigh wave dispersion data with azimuthal anisotropy can be used to invert for depth-dependent shear wavespeed azimuthal anisotropy, therefore reflecting depth-varying deformation patterns in the crust and upper mantle. In this study, we propose a two-step method that uses the Neighborhood Algorithm （NA） for the point-wise inversion of depth-dependent shear wavespeeds and azimuthal anisotropy from Rayleigh wave azimuthally anisotropic dispersion data. The first step employs the NA to estimate depth- dependent Vsv （or the elastic parameter L） as well as their uncertainties from the isotropic part Rayleigh wave dispersion data. In the second step, we first adopt a difference scheme to compute approximate Rayleigh-wave phase velocity sensitivity kernels to azimuthally anisotropic parameters with respect to the velocity model obtained in the first step. Then we perform the NA to estimate the azimuthally anisotropic parameters Gc/L and Gs/L at depths separately from the corresponding cosine and sine terms of the azimuthally anisotropic dispersion data. Finally, we compute the depth-dependent magnitude and fast polariza- tion azimuth of shear wavespeed azimuthal anisotropy. The use of the global search NA and Bayesian analysis allows for more reliable estimates of depth-dependent shear wavespeeds and azimuthal anisotropy as well as their uncertainties.We illustrate the inversion method using the azimuthally anisotropic dispersion data in SE Tibet, where we find apparent changes of fast axes of shear wavespeed azimuthal anisotropy between the crust and uppermost mantle.展开更多
In this study,the cylindrical finite-volume method(FVM)is advanced for the efficient and high-precision simulation of the logging while drilling(LWD)orthogonal azimuth electromagnetic tool(OAEMT)response in a three-di...In this study,the cylindrical finite-volume method(FVM)is advanced for the efficient and high-precision simulation of the logging while drilling(LWD)orthogonal azimuth electromagnetic tool(OAEMT)response in a three-dimensional(3 D)anisotropic formation.To overcome the ill-condition and convergence problems arising from the low induction number,Maxwell’s equations are reformulated into a mixed Helmholtz equation for the coupled potentials in a cylindrical coordinate system.The electrical fi eld continuation method is applied to approximate the perfectly electrical conducting(PEC)boundary condition,to improve the discretization accuracy of the Helmholtz equation on the surface of metal mandrels.On the base,the 3 D FVM on Lebedev’s staggered grids in the cylindrical coordinates is employed to discretize the mixed equations to ensure good conformity with typical well-logging tool geometries.The equivalent conductivity in a non-uniform element is determined by a standardization technique.The direct solver,PARDISO,is applied to efficiently solve the sparse linear equation systems for the multi-transmitter problem.To reduce the number of calls to PARDISO,the whole computational domain is divided into small windows that contain multiple measuring points.The electromagnetic(EM)solutions produced by all the transmitters per window are simultaneously solved because the discrete matrix,relevant to all the transmitters in the same window,is changed.Finally,the 3 D FVM is validated against the numerical mode matching method(NMM),and the characteristics of both the coaxial and coplanar responses of the EM field tool are investigated using the numerical results.展开更多
Azimuthal electromagnetic(EM)logging while drilling(LWD)has been extensively used in high-angle and horizontal(HA/HZ)wells.However,due to the effects of formation anisotropy,accurate geosteering decision and formation...Azimuthal electromagnetic(EM)logging while drilling(LWD)has been extensively used in high-angle and horizontal(HA/HZ)wells.However,due to the effects of formation anisotropy,accurate geosteering decision and formation evaluations have become increasingly difficult.To quantitatively analyze the effect of anisotropy on tool responses and data processing,this paper investigates the sensitivity of EM LWD measurements to electric anisotropy and inversion accuracy via forward modeling and inversion.First,a sensitivity factor is defined to quantitatively analyze the sensitivity of the magnetic field components and synthetic signals to electric anisotropy.Then,azimuthal EM LWD responses in anisotropic layered formations are simulated,and the sensitivities to formation parameters for compensated and uncompensated tool configurations are comparatively analyzed.Finally,we discuss the effects of the inversion model on bed boundary inversion in anisotropic formations.Numerical simulation and inversion results show that azimuthal EM LWD can be significantly affected by electric anisotropy.Fortunately,by using a symmetrical compensation configuration,the sensitivity of the geosignals to electric anisotropy can be suppressed,and the boundary detection capability can be further enhanced.Anisotropy normally gives rise to separated resistivity curves and abnormal"horns";moreover,complicated nonlinear distortion can also arise in geosignals as the tool approaches a bed boundary.If anisotropy effects are ignored in the inversion process,the estimated bed boundary and formation resistivity are usually unreliable,which may mislead geosteering decisions.展开更多
If the thickness of coal seams and the lithology of both roofs and floors of coal seams have not changed at all or only a little, then it is thought that the elastic anisotropy of coal seams depends mainly on fracture...If the thickness of coal seams and the lithology of both roofs and floors of coal seams have not changed at all or only a little, then it is thought that the elastic anisotropy of coal seams depends mainly on fractures and obeys the horizontally symmetric model of an azimuth anisotropy. For a fixed offset, the amplitude A of the reflection P-wave and the cosine of 2φ has an approximately linear relation, (φ is the source-detector azimuth with respect to the fracture strike. Based on this relationship, many things can be done, such as the extraction of macro bins, the correction of residual normal moveout, the formation of azimuth gather, the transformation and normalization of azimuth gathers and the extraction of reflection wave amplitudes of coal seams. The least squares method was used to inverse theoretically the direction and density of fractures of coal seams. The result is in good agreement with the regional geological structure, indicating that the azimuth anisotropic analysis of the P-wave is feasible in evaluating the density and direction of fractures in coal seams.展开更多
Acquisition footprint is a new concept to describe the seismic noise in three-dimensional seismic exploration and it is closely related to geometry and observation shuttering.At present,the study on acquisition footpr...Acquisition footprint is a new concept to describe the seismic noise in three-dimensional seismic exploration and it is closely related to geometry and observation shuttering.At present,the study on acquisition footprints has become a hot spot.In partnership with the Dagang Oilfield,we used the channel sand body seismic physical model to study the characteristics of wide/narrow azimuth acquisition footprints and analyzed and compared the two types of footprints and their effects on target imaging.In addition,the footprints caused by data processing of the normal moveout offset(NMO)stretching aberration were discussed.These footprints are located only in the shallow or middle layer in the time slice,and possibly affect the imaging of shallow target layers,and have no influence on deep target imaging.Seismic physical modeling has its advantages in the study of acquisition footprints.展开更多
基金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.
基金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 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.
基金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 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.
基金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.
文摘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.
文摘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.
文摘Fractured reservoirs always show anisotropic amplitude features,i.e.the reflection amplitude of seismic waves varies with offset and azimuth （AVOZ）.A noise attenuation fracture inversion algorithm is presented for fracture detection based on P-wave AVOZ.The conventional inversion method always fails when applied to limited azimuth data because of the existence of noise.In our inversion algorithm,special attention is paid to suppressing the noise during inversion,to overcome the limitation of the conventional inversion method on limited azimuth data.Numerical models are employed to illustrate the effectiveness of the method.The inversion algorithm is then applied to Tazhong 45 area field data which is acquired under limited azimuth distribution.Compared with cores and fullbore formation microimage （FMI）,the inverted results （fracture density and orientation） are reasonable,suggesting that the inversion algorithm is feasible for fracture prediction in the Tarim Basin.
基金co-funded by Chinese Postdoctoral Science Foundation(2018M640663)the National Natural Science Foundation of China(41474100,41574118,41674131)National Science and Technology Major Project of the Ministry of Science and Technology of China(2017ZX05009-001)
文摘We present systematic investigations on the physics,detection performance and inversion of logging-while-drilling extradeep azimuthal resistivity measurements(EDARM).First,the definitions of EDRAM measurements are discussed,followed by the derivation of the attenuation and phase-shift geometrical factors to illustrate the relative contributions of formation units to the observed signals.Then,a new definition of detection depth,which considers the uncertainty of inversion results caused by the data noise,is proposed to quantify the detection capability of ED ARM.Finally,the B ayesian theory associated with Markov chain Monte Carlo sampling is introduced for fast processing of EDARM data.Numerical results show that ED ARM is capable of detecting the azimuth and distance of remote bed boundaries,and the detection capability increases with increasing spacing and resistivity contrast.The EDARM tool can accommodate a large range of formation resistivity and is able to provide the resistivity anisotropy at arbitrary relative dipping angles.In addition,multiple bed boundaries and reservoir images near the borehole are readily obtained by using the Bayesian inversion.
文摘In this article, a method of high order, staggered grid, finite difference was used to model snapshots of seismic wave propagation and multi-azimuth surface seismic records in viscoelastic cracked monoclinic media. The modeling results clearly show the anisotropy caused by the crack fill material and the viscoelastic wave field through the monoclinic media. Summarizing the regularity of seismic records for different azimuths, this study is a significant attempt to understand the wave fields propagated in real media and to propose a more rational theoretical model for the future.
基金supported by National Natural Science Foundation of China (Grant number 41222028)China National Special Fund for Earthquake Scientific Research in Public Interest (Grant number 201508008)the Fundamental Research Funds for the Central Universities (WK2080000053)
文摘Seismic anisotropy provides important constraints on deformation patterns of Earth＇s material. Rayleigh wave dispersion data with azimuthal anisotropy can be used to invert for depth-dependent shear wavespeed azimuthal anisotropy, therefore reflecting depth-varying deformation patterns in the crust and upper mantle. In this study, we propose a two-step method that uses the Neighborhood Algorithm （NA） for the point-wise inversion of depth-dependent shear wavespeeds and azimuthal anisotropy from Rayleigh wave azimuthally anisotropic dispersion data. The first step employs the NA to estimate depth- dependent Vsv （or the elastic parameter L） as well as their uncertainties from the isotropic part Rayleigh wave dispersion data. In the second step, we first adopt a difference scheme to compute approximate Rayleigh-wave phase velocity sensitivity kernels to azimuthally anisotropic parameters with respect to the velocity model obtained in the first step. Then we perform the NA to estimate the azimuthally anisotropic parameters Gc/L and Gs/L at depths separately from the corresponding cosine and sine terms of the azimuthally anisotropic dispersion data. Finally, we compute the depth-dependent magnitude and fast polariza- tion azimuth of shear wavespeed azimuthal anisotropy. The use of the global search NA and Bayesian analysis allows for more reliable estimates of depth-dependent shear wavespeeds and azimuthal anisotropy as well as their uncertainties.We illustrate the inversion method using the azimuthally anisotropic dispersion data in SE Tibet, where we find apparent changes of fast axes of shear wavespeed azimuthal anisotropy between the crust and uppermost mantle.
基金supported jointly by Strategic Pilot Science and Technology Project of Chinese Academy of Sciences (No. XDA14020102)National key research and development plan (No. 2017YFC0601805)+5 种基金National Natural Science Foundation of China (No. 41574110)Youth Foundation of Hebei Educational Committee (No. QN2018217)Hebei Higher Education Teaching Reform Research and Practice(No. 2018GJJG328)Zhangjiakou science and technology bureau(No. 1821011B)Doctoral Fund of Hebei Institute of Architecture and Civil Engineering (No. B-201606)Academic Team Innovation Ability Improvement Project of Hebei Institute of Architecture and Civil Engineering(TD202011)。
文摘In this study,the cylindrical finite-volume method(FVM)is advanced for the efficient and high-precision simulation of the logging while drilling(LWD)orthogonal azimuth electromagnetic tool(OAEMT)response in a three-dimensional(3 D)anisotropic formation.To overcome the ill-condition and convergence problems arising from the low induction number,Maxwell’s equations are reformulated into a mixed Helmholtz equation for the coupled potentials in a cylindrical coordinate system.The electrical fi eld continuation method is applied to approximate the perfectly electrical conducting(PEC)boundary condition,to improve the discretization accuracy of the Helmholtz equation on the surface of metal mandrels.On the base,the 3 D FVM on Lebedev’s staggered grids in the cylindrical coordinates is employed to discretize the mixed equations to ensure good conformity with typical well-logging tool geometries.The equivalent conductivity in a non-uniform element is determined by a standardization technique.The direct solver,PARDISO,is applied to efficiently solve the sparse linear equation systems for the multi-transmitter problem.To reduce the number of calls to PARDISO,the whole computational domain is divided into small windows that contain multiple measuring points.The electromagnetic(EM)solutions produced by all the transmitters per window are simultaneously solved because the discrete matrix,relevant to all the transmitters in the same window,is changed.Finally,the 3 D FVM is validated against the numerical mode matching method(NMM),and the characteristics of both the coaxial and coplanar responses of the EM field tool are investigated using the numerical results.
基金supported by the National Natural Science Foundation of China(No.41674131,No.41974146,and No.41904109)the Shandong Province Postdoctoral Innovation Projects(sdbh20180025)the Fundamental Research Funds for the Central Universities(No.17CX06041)。
文摘Azimuthal electromagnetic(EM)logging while drilling(LWD)has been extensively used in high-angle and horizontal(HA/HZ)wells.However,due to the effects of formation anisotropy,accurate geosteering decision and formation evaluations have become increasingly difficult.To quantitatively analyze the effect of anisotropy on tool responses and data processing,this paper investigates the sensitivity of EM LWD measurements to electric anisotropy and inversion accuracy via forward modeling and inversion.First,a sensitivity factor is defined to quantitatively analyze the sensitivity of the magnetic field components and synthetic signals to electric anisotropy.Then,azimuthal EM LWD responses in anisotropic layered formations are simulated,and the sensitivities to formation parameters for compensated and uncompensated tool configurations are comparatively analyzed.Finally,we discuss the effects of the inversion model on bed boundary inversion in anisotropic formations.Numerical simulation and inversion results show that azimuthal EM LWD can be significantly affected by electric anisotropy.Fortunately,by using a symmetrical compensation configuration,the sensitivity of the geosignals to electric anisotropy can be suppressed,and the boundary detection capability can be further enhanced.Anisotropy normally gives rise to separated resistivity curves and abnormal"horns";moreover,complicated nonlinear distortion can also arise in geosignals as the tool approaches a bed boundary.If anisotropy effects are ignored in the inversion process,the estimated bed boundary and formation resistivity are usually unreliable,which may mislead geosteering decisions.
基金Projects 40574058 supported by the National Natural Science Foundation of China, 2005cb221500 the National Key Basic Research and Development(973) Program of China and 03(2007) the Scientific and Technological Project about Geology and Mineral Resources of Henan Land Resources Department
文摘If the thickness of coal seams and the lithology of both roofs and floors of coal seams have not changed at all or only a little, then it is thought that the elastic anisotropy of coal seams depends mainly on fractures and obeys the horizontally symmetric model of an azimuth anisotropy. For a fixed offset, the amplitude A of the reflection P-wave and the cosine of 2φ has an approximately linear relation, (φ is the source-detector azimuth with respect to the fracture strike. Based on this relationship, many things can be done, such as the extraction of macro bins, the correction of residual normal moveout, the formation of azimuth gather, the transformation and normalization of azimuth gathers and the extraction of reflection wave amplitudes of coal seams. The least squares method was used to inverse theoretically the direction and density of fractures of coal seams. The result is in good agreement with the regional geological structure, indicating that the azimuth anisotropic analysis of the P-wave is feasible in evaluating the density and direction of fractures in coal seams.
基金Heterogeneous formation geophysical response characteristics (973 Program, subject number: 2007CB209601)Continental reservoir seismic physical model (CNPC Fundamental Research Projects, subject number: 06A10102)
文摘Acquisition footprint is a new concept to describe the seismic noise in three-dimensional seismic exploration and it is closely related to geometry and observation shuttering.At present,the study on acquisition footprints has become a hot spot.In partnership with the Dagang Oilfield,we used the channel sand body seismic physical model to study the characteristics of wide/narrow azimuth acquisition footprints and analyzed and compared the two types of footprints and their effects on target imaging.In addition,the footprints caused by data processing of the normal moveout offset(NMO)stretching aberration were discussed.These footprints are located only in the shallow or middle layer in the time slice,and possibly affect the imaging of shallow target layers,and have no influence on deep target imaging.Seismic physical modeling has its advantages in the study of acquisition footprints.
