Ocean-Bottom Node(OBN)acquisitions provide both non-converted and converted reflection energy.There is a clear advantage to independently imaging both P-and S-waves,as they provide more information collectively than e...Ocean-Bottom Node(OBN)acquisitions provide both non-converted and converted reflection energy.There is a clear advantage to independently imaging both P-and S-waves,as they provide more information collectively than either does alone.In many conventional converted-wave pre-stack migration algorithms,density is treated as a constant,which is not the real-life case on earth.S-wave velocity and density information is crucial for hydrocarbon detection because it helps in the identification of porefilling fluids.In this paper,we focused on the effect of density on imaging,and developed a method of reverse-time migration(RTM)on converted s-waves of varying densities(VD-RTMCS).Phase correction was required prior to pre-stack migration to avoid constructive interference between data from adjacent sources.Synthetic data examples showed that when density variations were included,image profiles showed advantages in signal-to-noise ratio,vertical resolution and imaging of complex structures.展开更多
Multi-component seismic exploration technology, combining reflected PP- and converted PS-waves, is an effective tool for solving complicated oil and gas exploration problems. The improvement of converted wave resoluti...Multi-component seismic exploration technology, combining reflected PP- and converted PS-waves, is an effective tool for solving complicated oil and gas exploration problems. The improvement of converted wave resolution is one of the key problems. The main factor affecting converted wave resolution is the absorption of seismic waves in overlying strata. In order to remove the effect of absorption on converted waves, inverse Q filtering is used to improve the resolution. In this paper, we present a method to estimate the S-wave Q values from prestack converted wave gathers. Furthermore, we extend a stable and effective poststack inverse Q filtering method to prestack data which uses wave field continuation along the ray path to compensate for attenuation in prestack common shot PP- and PS-waves. The results of theoretical modeling prove that the method of estimating the S-wave Q values has high precision. The results from synthetic and real data prove that the stable inverse Q filtering method can effectively improve the resolution of prestack PP- and PS-waves.展开更多
Based on the empirical Gardner equation describing the relationship between density and compressional wave velocity, the converted wave reflection coefficient extrema attributes for AVO analysis are proposed and the r...Based on the empirical Gardner equation describing the relationship between density and compressional wave velocity, the converted wave reflection coefficient extrema attributes for AVO analysis are proposed and the relations between the extrema position and amplitude, average velocity ratio across the interface, and shear wave reflection coefficient are derived. The extrema position is a monotonically decreasing function of average velocity ratio, and the extrema amplitude is a function of average velocity ratio and shear wave reflection coefficient. For theoretical models, the average velocity ratio and shear wave reflection coefficient are inverted from the extrema position and amplitude obtained from fitting a power function to converted wave AVO curves. Shear wave reflection coefficient sections have clearer physical meaning than conventional converted wave stacked sections and establish the theoretical foundation for geological structural interpretation and event correlation. "The method of inverting average velocity ratio and shear wave reflection coefficient from the extrema position and amplitude obtained from fitting a power function is applied to real CCP gathers. The inverted average velocity ratios are consistent with those computed from compressional and shear wave well logs.展开更多
In this paper, we use the method of pseudo-offset migration (POM) to complete converted wave pre-stack time migration with amplitude-preservation in an anisotropic medium. The method maps the original traces into co...In this paper, we use the method of pseudo-offset migration (POM) to complete converted wave pre-stack time migration with amplitude-preservation in an anisotropic medium. The method maps the original traces into common conversion scatter point (CCSP) gathers directly by POM, which simplifies the conventional processing procedure for converted waves. The POM gather fold and SNR are high, which is favorable for velocity analysis and especially suitable for seismic data with low SNR. We used equivalent anisotropic theory to compute anisotropic parameters. Based on the scattering wave traveltime equation in a VTI medium, the POM pseudo-offset migration in anisotropic media was deduced. By amplitude-preserving POM gather mapping, velocity analysis, stack processing, and so on, the anisotropic migration results were acquired. The forward modeling computation and actual data processing demonstrate the validity of converted wave pre-stack time migration with amplitude-preservation using the anisotropic POM method.展开更多
PS converted-waves (C-waves) have been commonly used to image through gas clouds but the C-wave imaging may also be degraded by the diodic effect introduced by the gas cloud. It may be compensated for using a veloci...PS converted-waves (C-waves) have been commonly used to image through gas clouds but the C-wave imaging may also be degraded by the diodic effect introduced by the gas cloud. It may be compensated for using a velocity perturbation method which decouples the diodic moveout into two parts: the base velocity and the velocity perturbation. Gas clouds are widely distributed in the Sanhu area in the Qaidam basin of northwest China which is rich in natural gas. A land 2D3C seismic dataset is analyzed from the Sanhu area and significant diodic effects are observed in the data which harm the C-wave imaging. The diodic correction is applied to this data and the resultant C-wave imaging and the details of the reservoir structure are significantly improved. The diodic moveout plays an important role in working out the residu~ shear wave statics and the association of diodie correction and shear wave residual statics computation is a key step of C-wave high fidelity imaging in the gas cloud area. Finally, the new process workflow with diodic moveout is given.展开更多
Tracking problem in spherical coordinates with range rate (Doppler) measurements, which would have errors correlated to the range measurement errors, is investigated in this paper. The converted Doppler measurements...Tracking problem in spherical coordinates with range rate (Doppler) measurements, which would have errors correlated to the range measurement errors, is investigated in this paper. The converted Doppler measurements, constructed by the product of the Doppler measurements and range measurements, are used to replace the original Doppler measurements. A de-noising method based on an unbiased Kalman filter (KF) is proposed to reduce the converted Doppler measurement errors before updating the target states for the constant velocity (CV) model. The states from the de-noising filter are then combined with the Cartesian states from the converted measurement Kalman filter (CMKF) to produce final state estimates. The nonlinearity of the de-noising filter states are handled by expanding them around the Cartesian states from the CMKF in a Taylor series up to the second order term. In the mean time, the correlation between the two filters caused by the common range measurements is handled by a mini- mum mean squared error (MMSE) estimation-based method. These result in a new tracking filter, CMDN-EKF2. Monte Carlo simulations demonstrate that the proposed tracking filter can provide efficient and robust performance with a modest computa- tional cost.展开更多
Seismic velocity is important to migration of seismic data, interpretation of lithology and lithofacies as well as prediction of reservoir. The information of shear wave velocity is required to reduce the uncertainty ...Seismic velocity is important to migration of seismic data, interpretation of lithology and lithofacies as well as prediction of reservoir. The information of shear wave velocity is required to reduce the uncertainty for discriminating lithology, identifying fluid type in porous material and calculating gas saturation in reservoir prediction. Based on Zoeppritz equations, a numeral and scanning method was proposed in this paper. Shear wave velocity can be calculated with prestack converted wave data. The effects were demonstrated by inversion of theoretical and real seismic data.展开更多
Conventional land vertical seismic profiling (VSP) exploration usually uses P-wave sources and three-component geophones for receivers, emphasizing P- and converted S-waves. Previous studies show that both dynamite ...Conventional land vertical seismic profiling (VSP) exploration usually uses P-wave sources and three-component geophones for receivers, emphasizing P- and converted S-waves. Previous studies show that both dynamite borehole shots and vertical vibrations from controllable seismic sources at the surface will produce relatively strong pure P-waves and weaker pure S-waves. Interfaces with a large Poisson's ratio difference have a positive influence on the formation of strong transmitted converted S-waves. By a comparative analysis of pure S-waves from sources and converted downgoing S-waves, we believe that the main frequency of pure S-waves is usually lower than pure P-waves while the main frequency of downgoing converted S-waves is close to that of P-waves. We have studied zero-offset and offset VSP data from land P-wave sources. Results show that pure S-waves commonly exist in these data with differences in wave intensity. S-wave velocity can be obtained from the P-wave source zero-offset VSP data. Finally, we discuss the bright future of joint application of VSP P-and S-waves and the full use of S-waves in P-wave source VSP data.展开更多
In transversely isotropic media with a vertical symmetry axis (VTI), the converted-wave (C-wave) moveout over intermediate-to-far offsets is determined by four parameters. These are the C-wave stacking velocity Vc...In transversely isotropic media with a vertical symmetry axis (VTI), the converted-wave (C-wave) moveout over intermediate-to-far offsets is determined by four parameters. These are the C-wave stacking velocity Vc2 , the vertical and effective velocity ratios γ0 and γeff, and the anisotropic parameter χeff. We refer to the four parameters as the C-wave stacking velocity model. The purpose of C-wave velocity analysis is to determine this stacking velocity model. The C-wave stacking velocity model Vc2, γ0, γeff, and χeff can be determined from P-and C-wave reflection moveout data. However, error propagation is a severe problem in C-wave reflection-moveout inversion. The current short-spread stacking velocity as deduced from hyperbolic moveout does not provide sufficient accuracy to yield meaningful inverted values for the anisotropic parameters. The non-hyperbolic moveout over intermediate-offsets (x/z from 1.0 to 1.5) is no longer negligible and can be quantified using a background γ. Non-hyperbolic analysis with a γ correction over the intermediate offsets can yield Vc2 with errors less than 1% for noise free data. The procedure is very robust, allowing initial guesses of γ with up to 20% errors. It is also applicable for vertically inhomogeneous anisotropic media. This improved accuracy makes it possible to estimate anisotropic parameters using 4C seismic data. Two practical work flows are presented for this purpose: the double-scanning flow and the single-scanning flow. Applications to synthetic and real data show that the two flows yield results with similar accuracy but the single-scanning flow is more efficient than the double-scanning flow.展开更多
This paper describes a new algorithm for the 2-D converted-measurement Kalman filter (CMKF) which estimates a target’s Cartesian state given polar position measurements. At each processing index, the new algorithm ch...This paper describes a new algorithm for the 2-D converted-measurement Kalman filter (CMKF) which estimates a target’s Cartesian state given polar position measurements. At each processing index, the new algorithm chooses the more accurate of (1) the sensor’s polar position measurement and (2) the CMKF’s Cartesian position prediction. The new algorithm then computes the raw converted measurement’s error bias and the corresponding debiased converted measurement’s error covariance conditioned on the chosen position estimate. The paper derives explicit expressions for the polar-measurement-conditioned bias and covariance and shows the resulting polar-measurement-conditioned CMKF’s mathematical equivalence with the 2-D modified unbiased CMKF (MUCMKF). The paper also describes a method, based upon the unscented transformation, for approximating the raw converted measurement’s error bias and the debiased converted measurement’s error covariance conditioned on the CMKF’s Cartesian position prediction. Simulation results demonstrate the new CMKF’s improved tracking performance and statistical credibility as compared to those of the 2-D MUCMKF.展开更多
3-D converted-wave data were acquired using digital MEMS (micro-electromechanical system) three component (3C) sensors in the alternating sand and shale sequence in the overburden of the Shengli Ken-71 area. This ...3-D converted-wave data were acquired using digital MEMS (micro-electromechanical system) three component (3C) sensors in the alternating sand and shale sequence in the overburden of the Shengli Ken-71 area. This gives rise to serious non-hyperbolic moveout effects in the converted-wave data due to both the asymmetrical ray path and anisotropic effects. Conventional velocity analysis and moveout correction based on isotropic methods do not flatten reflections events. Here, we use a four-parameter theory to evaluate these effects and process the data. These four parameters include the PS converted wave stacking velocity (Vc2), the vertical velocity ratio (Y0), the effective velocity ratio (Yeff), and the anisotropy parameter (xoff), The method utilizes the moveout information at different offsets to estimate the different parameters and ensures that the events are properly aligned for stacking, As a result, this four-parameter theory leads to an improvement in imaging quality and correlation between the P-waves and converted-waves.展开更多
Studies of converted S-wave data recorded on the ocean bottom seismometer(OBS)allow for the estimation of crustal S-wave velocity,from which is further derived the Vp/Vs ratio to constrain the crustal lithology and ge...Studies of converted S-wave data recorded on the ocean bottom seismometer(OBS)allow for the estimation of crustal S-wave velocity,from which is further derived the Vp/Vs ratio to constrain the crustal lithology and geophysical properties.Constructing a precise S-wave velocity model is important for deep structural research,and inversion of converted S-waves provides a potential solution.However,the inversion of the converted S-wave remains a weakness because of the complexity of the seismic ray path and the inconsistent conversion interface.In this study,we introduced two travel time correction methods for the S-wave velocity inversion and imaged different S-wave velocity structures in accordance with the corresponding corrected S-wave phases using seismic data of profile EW6 in the northeastern South China Sea(SCS).The two inversion models show a similar trend in velocities,and the velocity difference is<0.15 km/s(mostly in the range of 0–0.1 km/s),indicating the accuracy of the two travel time correction methods and the reliability of the inversion results.According to simulations of seismic ray tracing based on different models,the velocity of sediments is the primary influencing factor in ray tracing for S-wave phases.If the sedimentary layer has high velocities,the near offset crustal S-wave refractions cannot be traced.In contrast,the ray tracing of Moho S-wave reflections was not significantly impacted by the velocity of the sediments.The two travel time correction methods have their own advantages,and the application of different approaches is based on additional requirements.These works provide an important reference for future improvements in converted S-wave research.展开更多
At a sampling rate of 100 samples per second,the YRY-4 four-gauge borehole strainmeters(FGBS) are capable of recording transient strains caused by seismic waves such as P and S waves or strain seismograms. At such a...At a sampling rate of 100 samples per second,the YRY-4 four-gauge borehole strainmeters(FGBS) are capable of recording transient strains caused by seismic waves such as P and S waves or strain seismograms. At such a high sampling rate, data from the YRY-4 strainmeters demonstrate fairly satisfactory self-consistency. The strain tensor seismograms demonstrate the senses of motion of P waves, that is, the type of seismic wave travels in the direction of the maximum normal strain change. The observed strain patterns of S waves significantly differ from those of P waves and should contain information about the source mechanism. Spectrum analysis shows that the strain seismograms are consistent with conventional broadband seismograms from the same site.展开更多
The present paper contributes in studying the phase velocities of P- and S-waves in a half space subjected to a compressive initial stress and gravity field. The density and acceleration due to gravity vary quadratica...The present paper contributes in studying the phase velocities of P- and S-waves in a half space subjected to a compressive initial stress and gravity field. The density and acceleration due to gravity vary quadratically along the depth. The dispersion equation is derived in a closed form. It is shown that the phase velocities depend not only on the initial stress, gravity, and direction of propagation but also on the inhomogeneity parameter associated with the density and acceleration due to gravity. Various particular cases are obtained, and the results match with the classical results. Numerical investigations on the phase velocities of P- and S-waves against the wave number are made for various sets of values of the material parameters, and the results are illustrated graphically. The graphical user interface model is developed to generalize the effect.展开更多
Lead is a ubiquitous pollutant and Pb pollution is a global public health problem.Lead has been reported to induce multiple adverse effects,including reproductive toxicity,neurotoxicity,carcinogenicity,nephrotoxicity,...Lead is a ubiquitous pollutant and Pb pollution is a global public health problem.Lead has been reported to induce multiple adverse effects,including reproductive toxicity,neurotoxicity,carcinogenicity,nephrotoxicity,immunotoxicity,and hypertension[1].The traditional medical treatment available for Pb poisoning is chelation,which can save lives in individuals with very high blood Pb levels.The commonly used chelating agents include CaNa2EDTA and meso-2,3-dimercaptosuccinic acid.However,chelation therapy has strong short-term effects on the overall long-term management of Pb exposure.展开更多
The original algorithm for the 2-D debiased converted-measurement Kalman filter (CMKF) specified, with incorrect mathematical justification, a requirement for evaluating the average true bias and covari-ance with the ...The original algorithm for the 2-D debiased converted-measurement Kalman filter (CMKF) specified, with incorrect mathematical justification, a requirement for evaluating the average true bias and covari-ance with the best available polar estimate, rather than exclusively with the polar measurement. Even though this original algorithm yields better tracking performance than the debiased-CMKF algorithm which evaluates the average true bias and covariance exclusively with the polar measurement, this paper shows the specified requirement compromises the statistical consistency between the debiased converted measurement’s error and the average true covariance. To resolve this apparent contradiction, this paper provides the correct empirical explanation for the tracking-performance improvement obtained by the specified requirement.展开更多
An inversion method was applied to crustal earthquakes dataset to find S-wave attenuation characteristics beneath the Eastern Tohoku region of Japan. Accelerograms from 85 shallow crustal earthquakes up to 25 km depth...An inversion method was applied to crustal earthquakes dataset to find S-wave attenuation characteristics beneath the Eastern Tohoku region of Japan. Accelerograms from 85 shallow crustal earthquakes up to 25 km depth and magnitude range between 3.5 and 5.5 were analyzed to estimate the seismic quality factor Qs. A homogeneous attenuation model Qs for the wave propagation path was evaluated from spectral amplitudes, at 24 different frequencies between 0.5 and 20 Hz by using generalized inversion technique. To do this, non-parametric attenuation functions were calculated to observe spectral amplitude decay with hypocentral distance. Then, these functions were parameterized to estimate Qs. It was found that in Eastern Tohoku region, the Qs frequency dependence can be approximated with the function 33 f 1.22 within a frequency range between 0.5 and 20 Hz. However, the frequency dependence of Qs in the frequency range between 0.5 and 6 Hz is best approximated by Qs (f) = 36 f 0.94 showing relatively weaker frequency dependence as compared to the relation Qs (f) = 6 f^ 2.09 for the frequency range between 6 and 15 Hz. These results could be used to estimate source and site parameters for seismic hazard assessment in the region.展开更多
With rapid growth of power demand, transmission capacity is also in urgent need of upgrading. In some cases, converting existing AC transmission lines to DC lines can Improve the transmission capacity and reduce the c...With rapid growth of power demand, transmission capacity is also in urgent need of upgrading. In some cases, converting existing AC transmission lines to DC lines can Improve the transmission capacity and reduce the construction investment. In this paper, the upstream finite element method was expanded to calculate the total electric field of same tower multi-circuit DC lines converted from double-circuit AC lines, and the validity of the algorithm was confirmed by experiments. Taking a DC line converted from a typical same tower 500 kV double-circuit AC transmission line as an example, the surface electric field and the ground total electric field in different pole conductor arrangement schemes were calculated and analyzed, and the critical height of pole conductors for DC lines in residential and non-residential area were determined. Then, the corridor width of DC and AC lines at critical height in residential and non-residential areas before and after AC-DC line transformation were compared. The results indicate that for DC lines converted from common 500 kV double-circuit AC lines, the ground total electric field can meet the requirements of corresponding standard with appropriate pole conductor arrangement schemes.展开更多
The present paper develops an approach of relative orbit determination for satellite formation flight.Inter-satellite measurements by the onboard devices of the satellite were chosen to perform this relative navigatio...The present paper develops an approach of relative orbit determination for satellite formation flight.Inter-satellite measurements by the onboard devices of the satellite were chosen to perform this relative navigation,and the equations of relative motion expressed in the Earth Centered Inertial frame were used to eliminate the assumption of the circular reference orbit.The relative orbit estimation was achieved through a continuous-discrete converted measurement Kalman filter design,in which the measurements were transformed to the inertial frame to avoid the linearization error of the observation equation.In addition,the situation of the coarse measurement period(only microwave radar measurements are available)existing was analyzed.The numerical simulation results verify the validity of the navigation approach,and it has been proved that this approach can be applied to the formation with an elliptical reference orbit.展开更多
The global adoption of Electric Vehicles(EVs)is on the rise due to their advanced features,with projections indicating they will soon dominate the private vehicle market.However,improper management of EV charging can ...The global adoption of Electric Vehicles(EVs)is on the rise due to their advanced features,with projections indicating they will soon dominate the private vehicle market.However,improper management of EV charging can lead to significant issues.This paper reviews the development of high-power,reliable charging solutions by examining the converter topologies used in rectifiers and converters that transfer electricity from the grid to EV batteries.It covers technical details,ongoing developments,and challenges related to these topologies and control strategies.The integration of rapid charging stations has introduced various Power Quality(PQ)issues,such as voltage fluctuations,harmonic distortion,and supra-harmonics,which are discussed in detail.The paper also highlights the benefits of controlled EV charging and discharging,including voltage and frequency regulation,reactive power compensation,and improved power quality.Efficient energy management and control strategies are crucial for optimizing EV battery charging within microgrids to meet increasing demand.Charging stations must adhere to specific converter topologies,control strategies,and industry standards to function correctly.The paper explores microgrid architectures and control strategies that integrate EVs,energy storage units(ESUs),and Renewable Energy Sources(RES)to enhance performance at charging points.It emphasizes the importance of various RES-connected architectures and the latest power converter topologies.Additionally,the paper provides a comparative analysis of microgrid-based charging station architectures,focusing on energy management,control strategies,and charging converter controls.The goal is to offer insights into future research directions in EV charging systems,including architectural considerations,control factors,and their respective advantages and disadvantages.展开更多
基金supported by the National Science and Technology Major Project (No. 2016ZX05027-002)the National Natural Science Foundation of China (No. 41230 318)Qingdao National Laboratory for Marine Science and Technology Innovation Project of Ao-Shan (No. 2015ASKJ03)
文摘Ocean-Bottom Node(OBN)acquisitions provide both non-converted and converted reflection energy.There is a clear advantage to independently imaging both P-and S-waves,as they provide more information collectively than either does alone.In many conventional converted-wave pre-stack migration algorithms,density is treated as a constant,which is not the real-life case on earth.S-wave velocity and density information is crucial for hydrocarbon detection because it helps in the identification of porefilling fluids.In this paper,we focused on the effect of density on imaging,and developed a method of reverse-time migration(RTM)on converted s-waves of varying densities(VD-RTMCS).Phase correction was required prior to pre-stack migration to avoid constructive interference between data from adjacent sources.Synthetic data examples showed that when density variations were included,image profiles showed advantages in signal-to-noise ratio,vertical resolution and imaging of complex structures.
基金supported by the 863 Program(Grant No.2007AA06Z218)
文摘Multi-component seismic exploration technology, combining reflected PP- and converted PS-waves, is an effective tool for solving complicated oil and gas exploration problems. The improvement of converted wave resolution is one of the key problems. The main factor affecting converted wave resolution is the absorption of seismic waves in overlying strata. In order to remove the effect of absorption on converted waves, inverse Q filtering is used to improve the resolution. In this paper, we present a method to estimate the S-wave Q values from prestack converted wave gathers. Furthermore, we extend a stable and effective poststack inverse Q filtering method to prestack data which uses wave field continuation along the ray path to compensate for attenuation in prestack common shot PP- and PS-waves. The results of theoretical modeling prove that the method of estimating the S-wave Q values has high precision. The results from synthetic and real data prove that the stable inverse Q filtering method can effectively improve the resolution of prestack PP- and PS-waves.
基金National 973 Key Basic Research Development Program (No.2005CB422104)SINOPEC's Scientific and Technological Development Program (No.P05063)
文摘Based on the empirical Gardner equation describing the relationship between density and compressional wave velocity, the converted wave reflection coefficient extrema attributes for AVO analysis are proposed and the relations between the extrema position and amplitude, average velocity ratio across the interface, and shear wave reflection coefficient are derived. The extrema position is a monotonically decreasing function of average velocity ratio, and the extrema amplitude is a function of average velocity ratio and shear wave reflection coefficient. For theoretical models, the average velocity ratio and shear wave reflection coefficient are inverted from the extrema position and amplitude obtained from fitting a power function to converted wave AVO curves. Shear wave reflection coefficient sections have clearer physical meaning than conventional converted wave stacked sections and establish the theoretical foundation for geological structural interpretation and event correlation. "The method of inverting average velocity ratio and shear wave reflection coefficient from the extrema position and amplitude obtained from fitting a power function is applied to real CCP gathers. The inverted average velocity ratios are consistent with those computed from compressional and shear wave well logs.
文摘In this paper, we use the method of pseudo-offset migration (POM) to complete converted wave pre-stack time migration with amplitude-preservation in an anisotropic medium. The method maps the original traces into common conversion scatter point (CCSP) gathers directly by POM, which simplifies the conventional processing procedure for converted waves. The POM gather fold and SNR are high, which is favorable for velocity analysis and especially suitable for seismic data with low SNR. We used equivalent anisotropic theory to compute anisotropic parameters. Based on the scattering wave traveltime equation in a VTI medium, the POM pseudo-offset migration in anisotropic media was deduced. By amplitude-preserving POM gather mapping, velocity analysis, stack processing, and so on, the anisotropic migration results were acquired. The forward modeling computation and actual data processing demonstrate the validity of converted wave pre-stack time migration with amplitude-preservation using the anisotropic POM method.
基金supported by the National Natural Science Foundation of China (No. 41074080)the Important National Science & Technology Specific Projects (No. 2011ZX05019-008)
文摘PS converted-waves (C-waves) have been commonly used to image through gas clouds but the C-wave imaging may also be degraded by the diodic effect introduced by the gas cloud. It may be compensated for using a velocity perturbation method which decouples the diodic moveout into two parts: the base velocity and the velocity perturbation. Gas clouds are widely distributed in the Sanhu area in the Qaidam basin of northwest China which is rich in natural gas. A land 2D3C seismic dataset is analyzed from the Sanhu area and significant diodic effects are observed in the data which harm the C-wave imaging. The diodic correction is applied to this data and the resultant C-wave imaging and the details of the reservoir structure are significantly improved. The diodic moveout plays an important role in working out the residu~ shear wave statics and the association of diodie correction and shear wave residual statics computation is a key step of C-wave high fidelity imaging in the gas cloud area. Finally, the new process workflow with diodic moveout is given.
基金Fundamental Research Funds for the Central Universities (HIT.NSRIF.2010097)National Natural Science Foundation of China (61171188, 61201311)
文摘Tracking problem in spherical coordinates with range rate (Doppler) measurements, which would have errors correlated to the range measurement errors, is investigated in this paper. The converted Doppler measurements, constructed by the product of the Doppler measurements and range measurements, are used to replace the original Doppler measurements. A de-noising method based on an unbiased Kalman filter (KF) is proposed to reduce the converted Doppler measurement errors before updating the target states for the constant velocity (CV) model. The states from the de-noising filter are then combined with the Cartesian states from the converted measurement Kalman filter (CMKF) to produce final state estimates. The nonlinearity of the de-noising filter states are handled by expanding them around the Cartesian states from the CMKF in a Taylor series up to the second order term. In the mean time, the correlation between the two filters caused by the common range measurements is handled by a mini- mum mean squared error (MMSE) estimation-based method. These result in a new tracking filter, CMDN-EKF2. Monte Carlo simulations demonstrate that the proposed tracking filter can provide efficient and robust performance with a modest computa- tional cost.
基金National 973 Key Basic Research Development Program (2005CB422104)SINOPEC's Scientific and Techno-logical Development Program(P05063)
文摘Seismic velocity is important to migration of seismic data, interpretation of lithology and lithofacies as well as prediction of reservoir. The information of shear wave velocity is required to reduce the uncertainty for discriminating lithology, identifying fluid type in porous material and calculating gas saturation in reservoir prediction. Based on Zoeppritz equations, a numeral and scanning method was proposed in this paper. Shear wave velocity can be calculated with prestack converted wave data. The effects were demonstrated by inversion of theoretical and real seismic data.
文摘Conventional land vertical seismic profiling (VSP) exploration usually uses P-wave sources and three-component geophones for receivers, emphasizing P- and converted S-waves. Previous studies show that both dynamite borehole shots and vertical vibrations from controllable seismic sources at the surface will produce relatively strong pure P-waves and weaker pure S-waves. Interfaces with a large Poisson's ratio difference have a positive influence on the formation of strong transmitted converted S-waves. By a comparative analysis of pure S-waves from sources and converted downgoing S-waves, we believe that the main frequency of pure S-waves is usually lower than pure P-waves while the main frequency of downgoing converted S-waves is close to that of P-waves. We have studied zero-offset and offset VSP data from land P-wave sources. Results show that pure S-waves commonly exist in these data with differences in wave intensity. S-wave velocity can be obtained from the P-wave source zero-offset VSP data. Finally, we discuss the bright future of joint application of VSP P-and S-waves and the full use of S-waves in P-wave source VSP data.
基金This work is funded by the Edinburgh Anisotropy Project of the British Geological Survey.
文摘In transversely isotropic media with a vertical symmetry axis (VTI), the converted-wave (C-wave) moveout over intermediate-to-far offsets is determined by four parameters. These are the C-wave stacking velocity Vc2 , the vertical and effective velocity ratios γ0 and γeff, and the anisotropic parameter χeff. We refer to the four parameters as the C-wave stacking velocity model. The purpose of C-wave velocity analysis is to determine this stacking velocity model. The C-wave stacking velocity model Vc2, γ0, γeff, and χeff can be determined from P-and C-wave reflection moveout data. However, error propagation is a severe problem in C-wave reflection-moveout inversion. The current short-spread stacking velocity as deduced from hyperbolic moveout does not provide sufficient accuracy to yield meaningful inverted values for the anisotropic parameters. The non-hyperbolic moveout over intermediate-offsets (x/z from 1.0 to 1.5) is no longer negligible and can be quantified using a background γ. Non-hyperbolic analysis with a γ correction over the intermediate offsets can yield Vc2 with errors less than 1% for noise free data. The procedure is very robust, allowing initial guesses of γ with up to 20% errors. It is also applicable for vertically inhomogeneous anisotropic media. This improved accuracy makes it possible to estimate anisotropic parameters using 4C seismic data. Two practical work flows are presented for this purpose: the double-scanning flow and the single-scanning flow. Applications to synthetic and real data show that the two flows yield results with similar accuracy but the single-scanning flow is more efficient than the double-scanning flow.
文摘This paper describes a new algorithm for the 2-D converted-measurement Kalman filter (CMKF) which estimates a target’s Cartesian state given polar position measurements. At each processing index, the new algorithm chooses the more accurate of (1) the sensor’s polar position measurement and (2) the CMKF’s Cartesian position prediction. The new algorithm then computes the raw converted measurement’s error bias and the corresponding debiased converted measurement’s error covariance conditioned on the chosen position estimate. The paper derives explicit expressions for the polar-measurement-conditioned bias and covariance and shows the resulting polar-measurement-conditioned CMKF’s mathematical equivalence with the 2-D modified unbiased CMKF (MUCMKF). The paper also describes a method, based upon the unscented transformation, for approximating the raw converted measurement’s error bias and the debiased converted measurement’s error covariance conditioned on the CMKF’s Cartesian position prediction. Simulation results demonstrate the new CMKF’s improved tracking performance and statistical credibility as compared to those of the 2-D MUCMKF.
文摘3-D converted-wave data were acquired using digital MEMS (micro-electromechanical system) three component (3C) sensors in the alternating sand and shale sequence in the overburden of the Shengli Ken-71 area. This gives rise to serious non-hyperbolic moveout effects in the converted-wave data due to both the asymmetrical ray path and anisotropic effects. Conventional velocity analysis and moveout correction based on isotropic methods do not flatten reflections events. Here, we use a four-parameter theory to evaluate these effects and process the data. These four parameters include the PS converted wave stacking velocity (Vc2), the vertical velocity ratio (Y0), the effective velocity ratio (Yeff), and the anisotropy parameter (xoff), The method utilizes the moveout information at different offsets to estimate the different parameters and ensures that the events are properly aligned for stacking, As a result, this four-parameter theory leads to an improvement in imaging quality and correlation between the P-waves and converted-waves.
基金The National Natural Science Foundation of China under contract Nos 42276062 and 42006071the Seismological Research Foundation for Youths of Guangdong Earthquake Agency under contract No.GDDZY202307+1 种基金the Strategic Priority Research Program of Chinese Academy of Sciences under contract No.XDA22020303the Science and Technology Planning Project of Guangdong Province-Guangdong Collaborative Innovation Center for Earthquake Prevention and Disaster Mitigation Technology under contract No.2018B020207011.
文摘Studies of converted S-wave data recorded on the ocean bottom seismometer(OBS)allow for the estimation of crustal S-wave velocity,from which is further derived the Vp/Vs ratio to constrain the crustal lithology and geophysical properties.Constructing a precise S-wave velocity model is important for deep structural research,and inversion of converted S-waves provides a potential solution.However,the inversion of the converted S-wave remains a weakness because of the complexity of the seismic ray path and the inconsistent conversion interface.In this study,we introduced two travel time correction methods for the S-wave velocity inversion and imaged different S-wave velocity structures in accordance with the corresponding corrected S-wave phases using seismic data of profile EW6 in the northeastern South China Sea(SCS).The two inversion models show a similar trend in velocities,and the velocity difference is<0.15 km/s(mostly in the range of 0–0.1 km/s),indicating the accuracy of the two travel time correction methods and the reliability of the inversion results.According to simulations of seismic ray tracing based on different models,the velocity of sediments is the primary influencing factor in ray tracing for S-wave phases.If the sedimentary layer has high velocities,the near offset crustal S-wave refractions cannot be traced.In contrast,the ray tracing of Moho S-wave reflections was not significantly impacted by the velocity of the sediments.The two travel time correction methods have their own advantages,and the application of different approaches is based on additional requirements.These works provide an important reference for future improvements in converted S-wave research.
基金supported by the Special Fund for Earthquake Research in the Public Interest(No.201108009)
文摘At a sampling rate of 100 samples per second,the YRY-4 four-gauge borehole strainmeters(FGBS) are capable of recording transient strains caused by seismic waves such as P and S waves or strain seismograms. At such a high sampling rate, data from the YRY-4 strainmeters demonstrate fairly satisfactory self-consistency. The strain tensor seismograms demonstrate the senses of motion of P waves, that is, the type of seismic wave travels in the direction of the maximum normal strain change. The observed strain patterns of S waves significantly differ from those of P waves and should contain information about the source mechanism. Spectrum analysis shows that the strain seismograms are consistent with conventional broadband seismograms from the same site.
基金supported by the Research Fellow of Indian School of Mines in Dhanbad (No. 2010DR0016)
文摘The present paper contributes in studying the phase velocities of P- and S-waves in a half space subjected to a compressive initial stress and gravity field. The density and acceleration due to gravity vary quadratically along the depth. The dispersion equation is derived in a closed form. It is shown that the phase velocities depend not only on the initial stress, gravity, and direction of propagation but also on the inhomogeneity parameter associated with the density and acceleration due to gravity. Various particular cases are obtained, and the results match with the classical results. Numerical investigations on the phase velocities of P- and S-waves against the wave number are made for various sets of values of the material parameters, and the results are illustrated graphically. The graphical user interface model is developed to generalize the effect.
基金sponsored by the Central Government Guides Local Scientific and Technological Development Fund Project(YDZX 20201400001443)Shanxi International Science and Technology Cooperation Project(201803D421065)+2 种基金National Natural Science Foundation of China(Grant No.30672621 and 81173473)Taiyuan City Science and Technology Project Special Talents Star Project(120247-08)Basic Research Project of Shanxi Province(202103021223241).
文摘Lead is a ubiquitous pollutant and Pb pollution is a global public health problem.Lead has been reported to induce multiple adverse effects,including reproductive toxicity,neurotoxicity,carcinogenicity,nephrotoxicity,immunotoxicity,and hypertension[1].The traditional medical treatment available for Pb poisoning is chelation,which can save lives in individuals with very high blood Pb levels.The commonly used chelating agents include CaNa2EDTA and meso-2,3-dimercaptosuccinic acid.However,chelation therapy has strong short-term effects on the overall long-term management of Pb exposure.
文摘The original algorithm for the 2-D debiased converted-measurement Kalman filter (CMKF) specified, with incorrect mathematical justification, a requirement for evaluating the average true bias and covari-ance with the best available polar estimate, rather than exclusively with the polar measurement. Even though this original algorithm yields better tracking performance than the debiased-CMKF algorithm which evaluates the average true bias and covariance exclusively with the polar measurement, this paper shows the specified requirement compromises the statistical consistency between the debiased converted measurement’s error and the average true covariance. To resolve this apparent contradiction, this paper provides the correct empirical explanation for the tracking-performance improvement obtained by the specified requirement.
基金a part of author’s M.Sc Research under the project:‘‘Strengthening of Earthquake Engineering Center’’,funded by Higher Education Commission,Government of Pakistan
文摘An inversion method was applied to crustal earthquakes dataset to find S-wave attenuation characteristics beneath the Eastern Tohoku region of Japan. Accelerograms from 85 shallow crustal earthquakes up to 25 km depth and magnitude range between 3.5 and 5.5 were analyzed to estimate the seismic quality factor Qs. A homogeneous attenuation model Qs for the wave propagation path was evaluated from spectral amplitudes, at 24 different frequencies between 0.5 and 20 Hz by using generalized inversion technique. To do this, non-parametric attenuation functions were calculated to observe spectral amplitude decay with hypocentral distance. Then, these functions were parameterized to estimate Qs. It was found that in Eastern Tohoku region, the Qs frequency dependence can be approximated with the function 33 f 1.22 within a frequency range between 0.5 and 20 Hz. However, the frequency dependence of Qs in the frequency range between 0.5 and 6 Hz is best approximated by Qs (f) = 36 f 0.94 showing relatively weaker frequency dependence as compared to the relation Qs (f) = 6 f^ 2.09 for the frequency range between 6 and 15 Hz. These results could be used to estimate source and site parameters for seismic hazard assessment in the region.
文摘With rapid growth of power demand, transmission capacity is also in urgent need of upgrading. In some cases, converting existing AC transmission lines to DC lines can Improve the transmission capacity and reduce the construction investment. In this paper, the upstream finite element method was expanded to calculate the total electric field of same tower multi-circuit DC lines converted from double-circuit AC lines, and the validity of the algorithm was confirmed by experiments. Taking a DC line converted from a typical same tower 500 kV double-circuit AC transmission line as an example, the surface electric field and the ground total electric field in different pole conductor arrangement schemes were calculated and analyzed, and the critical height of pole conductors for DC lines in residential and non-residential area were determined. Then, the corridor width of DC and AC lines at critical height in residential and non-residential areas before and after AC-DC line transformation were compared. The results indicate that for DC lines converted from common 500 kV double-circuit AC lines, the ground total electric field can meet the requirements of corresponding standard with appropriate pole conductor arrangement schemes.
基金Program for Changjiang Scholars and Innovative Research Team in University (PCSIRT)
文摘The present paper develops an approach of relative orbit determination for satellite formation flight.Inter-satellite measurements by the onboard devices of the satellite were chosen to perform this relative navigation,and the equations of relative motion expressed in the Earth Centered Inertial frame were used to eliminate the assumption of the circular reference orbit.The relative orbit estimation was achieved through a continuous-discrete converted measurement Kalman filter design,in which the measurements were transformed to the inertial frame to avoid the linearization error of the observation equation.In addition,the situation of the coarse measurement period(only microwave radar measurements are available)existing was analyzed.The numerical simulation results verify the validity of the navigation approach,and it has been proved that this approach can be applied to the formation with an elliptical reference orbit.
文摘The global adoption of Electric Vehicles(EVs)is on the rise due to their advanced features,with projections indicating they will soon dominate the private vehicle market.However,improper management of EV charging can lead to significant issues.This paper reviews the development of high-power,reliable charging solutions by examining the converter topologies used in rectifiers and converters that transfer electricity from the grid to EV batteries.It covers technical details,ongoing developments,and challenges related to these topologies and control strategies.The integration of rapid charging stations has introduced various Power Quality(PQ)issues,such as voltage fluctuations,harmonic distortion,and supra-harmonics,which are discussed in detail.The paper also highlights the benefits of controlled EV charging and discharging,including voltage and frequency regulation,reactive power compensation,and improved power quality.Efficient energy management and control strategies are crucial for optimizing EV battery charging within microgrids to meet increasing demand.Charging stations must adhere to specific converter topologies,control strategies,and industry standards to function correctly.The paper explores microgrid architectures and control strategies that integrate EVs,energy storage units(ESUs),and Renewable Energy Sources(RES)to enhance performance at charging points.It emphasizes the importance of various RES-connected architectures and the latest power converter topologies.Additionally,the paper provides a comparative analysis of microgrid-based charging station architectures,focusing on energy management,control strategies,and charging converter controls.The goal is to offer insights into future research directions in EV charging systems,including architectural considerations,control factors,and their respective advantages and disadvantages.