Multi-component seismic exploration is an important technique in the utilization of P-waves and converted S-waves for oil and gas exploration.It has unique advantages in the structural imaging of gas zones,reservoir p...Multi-component seismic exploration is an important technique in the utilization of P-waves and converted S-waves for oil and gas exploration.It has unique advantages in the structural imaging of gas zones,reservoir prediction,lithology,and gas-water identifi cation,and the development direction and degree of fractures.Multi-component joint inversion is one of the most important steps in multi-component exploration.In this paper,starting from the basic principle of multi-component joint inversion,the diff erences between the method and single P-wave inversion are introduced.Next,the technique is applied to the PLN area of the Sichuan Basin,and the P-wave impedance,S-wave impedance,and density are obtained based on multi-component joint inversion.Through the velocity and lithology,porosity,and gas saturation fi tting formulas,prediction results are calculated,and the results are analyzed.Finally,multi-component joint inversion and single P-wave inversion are compared in eff ective reservoir prediction.The results show that multi-component joint inversion increases the constraints on the inversion conditions,reduces the multi-solution of a single P-wave inversion,and is more objective and reliable for the identification of reservoirs,effectively improving the accuracy of oil and gas reservoir prediction and development.展开更多
With the development of gravity gradient full tensor measurement technique,three-dimensional( 3D) inversion based on gravity gradient tensor can provide more accurate information. But the forward calculation of 3D ful...With the development of gravity gradient full tensor measurement technique,three-dimensional( 3D) inversion based on gravity gradient tensor can provide more accurate information. But the forward calculation of 3D full tensor sensitivity matrix is very time-consuming,which restricts its development and application.According to the symmetry of the kernel function,the authors reconstruct the underground source of geological body to avoid repeat computation of the same value,and work out the corresponding relationship between the response of geological body to the observation point and the response of reconstructed geological body to the observation point. According to the relationship,rapid calculation of full tensor gravity sensitivity matrix can be achieved. The model calculation shows that this method can increase the speed of 30-45 times compared with the traditional calculation method. The sensitivity matrix is applied to the multi-component inversion of gravity gradient. The application of this method on the measured data provides the basis for the promotion of the method.展开更多
Young's modulus and Poisson's ratio are crucial parameters for reservoir characterization and rock brittleness evaluation.Conventional methods often rely on indirect computation or approximations of the Zoeppr...Young's modulus and Poisson's ratio are crucial parameters for reservoir characterization and rock brittleness evaluation.Conventional methods often rely on indirect computation or approximations of the Zoeppritz equations to estimate Young's modulus,which can introduce cumulative errors and reduce the accuracy of inversion results.To address these issues,this paper introduces the analytical solution of the Zoeppritz equation into the inversion process.The equation is re-derived and expressed in terms of Young's modulus,Poisson's ratio,and density.Within the Bayesian framework,we construct an objective function for the joint inversion of PP and PS waves.Traditional gradient-based algorithms often suffer from low precision and the computational complexity.In this study,we address limitations of conventional approaches related to low precision and complicated code by using Circle chaotic mapping,Levy flights,and Gaussian mutation to optimize the quantum particle swarm optimization(QPSO),named improved quantum particle swarm optimization(IQPSO).The IQPSO demonstrates superior global optimization capabilities.We test the proposed inversion method with both synthetic and field data.The test results demonstrate the proposed method's feasibility and effectiveness,indicating an improvement in inversion accuracy over traditional methods.展开更多
The Chilean Pampean flat slab subduction segment is characterized by the nearly horizontal subduction of the Nazca Plate within the depth range of 100-120 km.Numerous seismic tomography studies have been conducted to ...The Chilean Pampean flat slab subduction segment is characterized by the nearly horizontal subduction of the Nazca Plate within the depth range of 100-120 km.Numerous seismic tomography studies have been conducted to investigate its velocity structure;however,they have used only seismic body wave data or surface wave data.As a result,the existing velocity models in the region may have relatively large uncertainties.In this study,we use body wave arrival times from earthquakes occurring in central Chile between 2014 and 2019,as well as Rayleigh wave phase velocity maps at periods of 5-80 s from ambient noise empirical Green’s functions in Chile.By jointly using body wave arrival times and surface wave dispersion data,we refine the VS model and improve earthquake locations in the central Chile subduction zone.Compared with previous velocity models,our velocity model better reveals an eastward-dipping high-velocity plate representing the subducting Nazca Plate,which is 40-50 km thick and is more consistent with the slab thickness estimated by receiver function imaging and thermal modeling.Overall,the intraslab seismicity distribution spatially correlates well with the slab high-velocity anomalies except along the subduction paths of the CopiapóRidge and Juan Fernández Ridge.Additionally,parallel low-velocity stripes are imaged beneath the subducting plate,which are likely associated with the accumulated melts.The joint inversion velocity model also resolves widespread low-velocity anomalies in the crust beneath the Central Volcanic Zone of the central Andes,likely representing crustal magma chambers for various volcanoes.展开更多
Based on the Aki-Richards approximate equations for reflection coefficients and Bayes theorem, we developed an inversion method to estimate P- and S-wave velocity contrasts and density contrast from combined PP and PS...Based on the Aki-Richards approximate equations for reflection coefficients and Bayes theorem, we developed an inversion method to estimate P- and S-wave velocity contrasts and density contrast from combined PP and PS data. This method assumes that the parameters satisfy a normal distribution and introduces the covariance matrix to describe the degree of correlation between the parameters and thus to improve the inversion stability. Then, we suppose that the parameter sequence is subject to the Cauchy distribution and employs another matrix Q to describe the parameter sequence sparseness to improve the inversion result resolution. Tests on both synthetic and real multi-component data prove that this method is valid, efficient, more stable, and more accurate compared to methods using PP data only.展开更多
Multi-component exploration has many advantages over ordinary P-wave exploration. PP/PS joint AVO analysis and inversion are useful and powerful methods to discriminate between reservoir and non-productive lithology. ...Multi-component exploration has many advantages over ordinary P-wave exploration. PP/PS joint AVO analysis and inversion are useful and powerful methods to discriminate between reservoir and non-productive lithology. In this paper, we derive a new PS-wave reflection coefficient approximation equation which is more accurate at larger incidence angles. The equation is simplified for small incidence angles, which makes AVO analysis clearer and easier for angles less than 30 degrees. Based on this approximation, a PP/PS joint inversion is introduced. A real data example shows that oil sands, brine sands and shales can be differentiated based on the P- to S-wave velocity ratio from the PP/PS joint inversion. Fluid factors and Poisson＇s ratio also indicate an anomaly in the target zone at the oil well location.展开更多
Conventional joint PP-PS inversion is based on approximations of the Zoeppritz equations and assumes constant VP/VS;therefore,the inversion precision and stability cannot satisfy current exploration requirements.We pr...Conventional joint PP-PS inversion is based on approximations of the Zoeppritz equations and assumes constant VP/VS;therefore,the inversion precision and stability cannot satisfy current exploration requirements.We propose a joint PP-PS inversion method based on the exact Zoeppritz equations that combines Bayesian statistics and generalized linear inversion.A forward model based on the exact Zoeppritz equations is built to minimize the error of the approximations in the large-angle data,the prior distribution of the model parameters is added as a regularization item to decrease the ill-posed nature of the inversion,low-frequency constraints are introduced to stabilize the low-frequency data and improve robustness,and a fast algorithm is used to solve the objective function while minimizing the computational load.The proposed method has superior antinoising properties and well reproduces real data.展开更多
Based on the synchronous joint gravity and magnetic inversion of single interface by Pilkington and the need of revealing Cenozoic and crystalline basement thickness in the new round of oil-gas exploration, we propose...Based on the synchronous joint gravity and magnetic inversion of single interface by Pilkington and the need of revealing Cenozoic and crystalline basement thickness in the new round of oil-gas exploration, we propose a joint gravity and magnetic inversion methodfor two-layer models by concentrating on the relationship between the change of thicknessI and position of the middle layer and anomaly and discuss the effects of the key parameters. Model tests and application to field data show the validity of this method.展开更多
Multiwave seismic technology promotes the application of joint PP–PS amplitude versus offset (AVO) inversion;however conventional joint PP–PS AVO inversioan is linear based on approximations of the Zoeppritz equatio...Multiwave seismic technology promotes the application of joint PP–PS amplitude versus offset (AVO) inversion;however conventional joint PP–PS AVO inversioan is linear based on approximations of the Zoeppritz equations for multiple iterations. Therefore the inversion results of P-wave, S-wave velocity and density exhibit low precision in the faroffset;thus, the joint PP–PS AVO inversion is nonlinear. Herein, we propose a nonlinear joint inversion method based on exact Zoeppritz equations that combines improved Bayesian inference and a least squares support vector machine (LSSVM) to solve the nonlinear inversion problem. The initial parameters of Bayesian inference are optimized via particle swarm optimization (PSO). In improved Bayesian inference, the optimal parameter of the LSSVM is obtained by maximizing the posterior probability of the hyperparameters, thus improving the learning and generalization abilities of LSSVM. Then, an optimal nonlinear LSSVM model that defi nes the relationship between seismic refl ection amplitude and elastic parameters is established to improve the precision of the joint PP–PS AVO inversion. Further, the nonlinear problem of joint inversion can be solved through a single training of the nonlinear inversion model. The results of the synthetic data suggest that the precision of the estimated parameters is higher than that obtained via Bayesian linear inversion with PP-wave data and via approximations of the Zoeppritz equations. In addition, results using synthetic data with added noise show that the proposed method has superior anti-noising properties. Real-world application shows the feasibility and superiority of the proposed method, as compared with Bayesian linear inversion.展开更多
Accurate determination of seismic velocity of the crust is important for understanding regional tectonics and crustal evolution of the Earth. We propose a stepwise joint linearized inversion method using surface wave ...Accurate determination of seismic velocity of the crust is important for understanding regional tectonics and crustal evolution of the Earth. We propose a stepwise joint linearized inversion method using surface wave dispersion, Rayleigh wave ZH ratio （i.e., ellipticity）, and receiver function data to better resolve 1D crustal shear wave velocity （Vs） structure. Surface wave dispersion and Rayleigh wave ZH ratio data are more sensitive to absolute variations of shear wave speed at depths, but their sensi- tivity kernels to shear wave speeds are different and complimentary. However, receiver function data are more sensitive to sharp velocity contrast （e.g., due to the existence of crustal interfaces） and Vp/Vs ratios. The stepwise inversion method takes advantages of the complementary sensitivities of each dataset to better constrain the Vs model in the crust. We firstly invert surface wave dispersion and ZH ratio data to obtain a 1D smooth absolute vs model and then incorporate receiver function data in the joint inver- sion to obtain a finer Vs model with better constraints on interface structures. Through synthetic tests, Monte Carlo error analyses, and application to real data, we demonstrate that the proposed joint inversion method can resolve robust crustal Vs structures and with little initial model dependency.展开更多
Teleseismic and GPS data were jointly inverted for the rupture process of the 2011 Tohoku earthquake. The inversion results show that it is a bilateral rupture event with an average rupture velocity less than 2.0 km/s...Teleseismic and GPS data were jointly inverted for the rupture process of the 2011 Tohoku earthquake. The inversion results show that it is a bilateral rupture event with an average rupture velocity less than 2.0 km/s along the fault strike direction. The source rupture process consists of three sub-events, the first oc- curred near the hypocenter and the rest two ruptured along the up-dip direction and broke the sea bed, causing a maximum slip of about 30 m. The large-scale sea bed breakage may account for the tremendous tsunami disaster which resulted in most of the death and missing in this mega earthquake.展开更多
A two-dimensional forward and backward algorithm for the controlled-source audio-frequency magnetotelluric （CSAMT） method is developed to invert data in the entire region （near, transition, and far） and deal with ...A two-dimensional forward and backward algorithm for the controlled-source audio-frequency magnetotelluric （CSAMT） method is developed to invert data in the entire region （near, transition, and far） and deal with the effects of artificial sources. First, a regularization factor is introduced in the 2D magnetic inversion, and the magnetic susceptibility is updated in logarithmic form so that the inversion magnetic susceptibility is always positive. Second, the joint inversion of the CSAMT and magnetic methods is completed with the introduction of the cross gradient. By searching for the weight of the cross-gradient term in the objective function, the mutual influence between two different physical properties at different locations are avoided. Model tests show that the joint inversion based on cross-gradient theory offers better results than the single-method inversion. The 2D forward and inverse algorithm for CSAMT with source can effectively deal with artificial sources and ensures the reliability of the final joint inversion algorithm.展开更多
We use interferometric synthetic aperture radar （InSAR） and broadband seismic waveform data to estimate a source model of the 11th July, 2004 M W 6.2 Zhongba earthquake, Tibet of China. This event occurred within th...We use interferometric synthetic aperture radar （InSAR） and broadband seismic waveform data to estimate a source model of the 11th July, 2004 M W 6.2 Zhongba earthquake, Tibet of China. This event occurred within the seismically active zone of southwestern Tibetan Plateau where the east-west extension of the upper crust is observed. Because of limitations in one pair of InSAR data available, there are trade-offs among centroid depth, rupture area and amount of slip. Available seismic data tightly constrain the focal mechanism and centroid depth of the earthquake but not the horizontal location. Together, two complementary data sets can be used to identify the actual fault plane, better constrain the slip model and event location. We first use regional seismic waveform to estimate point source mechanism, then InSAR data is used to obtain better location. Finally, a joint inversion of teleseismic P-waves and InSAR data is performed to obtain a distributed model. Our preferred point source mechanism indicates a seismic moment of ～2.2×10 18 N·m （～M W 6.2）, a fault plane solution of 171° （342 ° ）/42 ° （48 ° ）/-83 ° （-97 ° ）, corresponding to strike/dip/rake, and a depth of 11 km. The fault plane with strike of 171 ？ and dip of 42° is identified as the ruptured fault with the aid of InSAR data. The preferred source model features compact area of slips between depth of 5–11 km and 10 km along strike with maximum slip amplitude of about 1.5 m.展开更多
Considering Zoeppritz equations, reflections of PP and PS are only the function of ratios of density and velocity. So the inversion results will be the same if the ratios are the same but values of density, velocities...Considering Zoeppritz equations, reflections of PP and PS are only the function of ratios of density and velocity. So the inversion results will be the same if the ratios are the same but values of density, velocities of P- wave and S-wave are different without strict constraint. This paper makes efforts to explore nonlinear simultaneous PP and PS inversion with expectation to reduce the ambiguity of AVO analysis by utilizing the redundancy of multi-component AVO measurements. Accurate estimation of ratio parameters depends on independence of input data. There are only two independent AVO attributes for PP reflectivity （i.e. intercept and gradient） and two for PS reflectivity （i.e. pseudo-intercept and pseudo-gradient or extreme amplitude）, respectively. For individual PP and PS inversion, the values of least-squares objective function do not converge around a large neighborhood of chosen true model parameters. Fortunately for joint PP and PS inversion the values of the least-squares objective function show closed contours with single minima. Finally the power function fitting is used to provide a higher precision AVO attributes than traditional polynomial fitting. By using the four independent fitting attributes （two independent attributes for PP and PS respectively）, the inversion of four ratio parameters （velocities and densities） would be estimated with less errors than that in traditional method.展开更多
Joint PP–PS inversion offers better accuracy and resolution than conventional P-wave inversion. P-and S-wave elastic moduli determined through data inversions are key parameters for reservoir evaluation and fluid cha...Joint PP–PS inversion offers better accuracy and resolution than conventional P-wave inversion. P-and S-wave elastic moduli determined through data inversions are key parameters for reservoir evaluation and fluid characterization. In this paper, starting with the exact Zoeppritz equation that relates P-and S-wave moduli, a coefficient that describes the reflections of P-and converted waves is established. This method effectively avoids error introduced by approximations or indirect calculations, thus improving the accuracy of the inversion results. Considering that the inversion problem is ill-posed and that the forward operator is nonlinear, prior constraints on the model parameters and modified low-frequency constraints are also introduced to the objective function to make the problem more tractable. This modified objective function is solved over many iterations to continuously optimize the background values of the velocity ratio, which increases the stability of the inversion process. Tests of various models show that the method effectively improves the accuracy and stability of extracting P and S-wave moduli from underdetermined data. This method can be applied to provide inferences for reservoir exploration and fluid extraction.展开更多
The coal-bearing strata of the deep Upper Paleozoic in the GS Sag have high hydrocarbon potential. Because of the absence of seismic data, we use electromagnetic （MT） and gravity data jointly to delineate the distri...The coal-bearing strata of the deep Upper Paleozoic in the GS Sag have high hydrocarbon potential. Because of the absence of seismic data, we use electromagnetic （MT） and gravity data jointly to delineate the distribution of deep targets based on well logging and geological data. First, a preliminary geological model is established by using three-dimensional （3D） MT inversion results. Second, using the formation density and gravity anomalies, the preliminary geological model is modified by interactive inversion of the gravity data. Then, we conduct MT-constrained inversion based on the modified model to obtain an optimal geological model until the deviations at all stations are minimized. Finally, the geological model and a seismic profile in the middle of the sag is analysed. We determine that the deep reflections of the seismic profile correspond to the Upper Paleozoic that reaches thickness up to 800 m. The processing of field data suggests that the joint MT-gravity modeling and constrained inversion can reduce the multiple solutions for single geophysical data and thus improve the recognition of deep formations. The MT-constrained inversion is consistent with the geological features in the seismic section. This suggests that the joint MT and gravity modeling and constrained inversion can be used to delineate deep targets in similar basins.展开更多
Yushu Ms7.1 earthquake occurred on the Ganzi-Yushu fault zone, across which we carried out a joint relative-gravity and seismic-reflection survey, and then performed a gravity inversion constrained by the seismic-refl...Yushu Ms7.1 earthquake occurred on the Ganzi-Yushu fault zone, across which we carried out a joint relative-gravity and seismic-reflection survey, and then performed a gravity inversion constrained by the seismic-reflection result. Based on the data of complete Bouguer gravity anomaly and seismic reflection, we obtained a layered interface structure in deep crust down to Moho. Our study showed that the inversion could reveal the interfaces of strata along the survey profile and the directions of regional faults in two-dimension. From the characteristics of the observed topography of the Moho basement, we tentatively confirmed that the uplift of eastern edge of Qinghai-Tibet plateau was caused by the subduetion of the Indian plate.展开更多
Joint inversion based on a correlation constraint utilizes a linear correlation function as a structural constraint.The linear correlation function contains a denominator,which may result in a singularity as the objec...Joint inversion based on a correlation constraint utilizes a linear correlation function as a structural constraint.The linear correlation function contains a denominator,which may result in a singularity as the objective function is optimized,leading to an unstable inversion calculation.To improve the robustness of this calculation,this paper proposes a new method in which a sinusoidal correlation function is employed as the structural constraint for joint inversion instead of the conventional linear correlation function.This structural constraint does not contain a denominator,thereby preventing a singularity.Compared with the joint inversion method based on a cross-gradient constraint,the joint inversion method based on a sinusoidal correlation constraint exhibits good performance.An application to actual data demonstrates that this method can process real data.展开更多
Most current prestack AVA joint inversion methods are based on the exact Zoeppritz equation and its various approximations. However, these equations only reflect the relation between reflection coefficients, incidence...Most current prestack AVA joint inversion methods are based on the exact Zoeppritz equation and its various approximations. However, these equations only reflect the relation between reflection coefficients, incidence angles, and elastic parameters on either side of the interface, which means that wave-propagation effects, such as spherical spreading, attenuation, transmission loss, multiples, and event mismatching of P-and S-waves, are not considered and cannot accurately describe the true propagation characteristics of seismic waves. Conventional AVA inversion methods require that these wave-propagation effects have been fully corrected or attenuated before inversion but these requirements can hardly be satisfied in practice. Using a one-dimensional(1 D) earth model, the reflectivity method can simulate the full wavefield response of seismic waves. Therefore, we propose a nonlinear multicomponent prestack AVA joint inversion method based on the vectorized reflectivity method, which uses a fast nondominated sorting genetic algorithm(NSGA II) to optimize the nonlinear multiobjective function to estimate multiple parameters, such as P-wave velocity, S-wave velocity, and density. This approach is robust because it can simultaneously cope with more than one objective function without introducing weight coefficients. Model tests prove the effectiveness of the proposed inversion method. Based on the inversion results, we find that the nonlinear prestack AVA joint inversion using the reflectivity method yields more accurate inversion results than the inversion by using the exact Zoeppritz equation when the wave-propagation effects of transmission loss and internal multiples are not completely corrected.展开更多
A new 3D velocity model of the crust and upper mantle in the southeastern （SE） margin of the Tibetan plateau was obtained by joint inversion of body- and sur- face-wave data. For the body-wave data, we used 7190 eve...A new 3D velocity model of the crust and upper mantle in the southeastern （SE） margin of the Tibetan plateau was obtained by joint inversion of body- and sur- face-wave data. For the body-wave data, we used 7190 events recorded by 102 stations in the SE margin of the Tibetan plateau. The surface-wave data consist of Rayleigh wave phase velocity dispersion curves obtained from ambient noise cross-correlation analysis recorded by a dense array in the SE margin of the Tibetan plateau. The joint inversion clearly improves the Vs model because it is constrained by both data types. The results show that at around 10 km depth there are two low-velocity anomalies embedded within three high-velocity bodies along the Longmenshan fault system. These high-velocity bodies correspond well with the Precambrian massifs, and the two located to the northeast of 2013 Ms 7.0 Lushan earthquake are associated with high fault slip areas during the 2008 Wenchuan earthquake. The aftershock gap between 2013 Lushan earthquake and 2008 Wenchuan earthquake is associated with low-velocity anomalies, which also acts as a barrier zone for ruptures of two earthquakes. Generally large earthquakes （M 〉 5） in the region occurring from 2008 to 2015 are located around the high-velocity zones, indicating that they may act as asperities for these large earthquakes. Joint inversion results also clearly show that there exist low-velocity or weak zones in the mid-lower crust, which are not evenly distributed beneath the SE margin of Tibetan plateau.展开更多
基金This work was supported by“Thirteenth Five-Year”national science and technology major Project(No.2017ZX05018005-004)CNPC fundamental research project(No.2016E-0604)National Natural Science Foundation of China(No.41374111).
文摘Multi-component seismic exploration is an important technique in the utilization of P-waves and converted S-waves for oil and gas exploration.It has unique advantages in the structural imaging of gas zones,reservoir prediction,lithology,and gas-water identifi cation,and the development direction and degree of fractures.Multi-component joint inversion is one of the most important steps in multi-component exploration.In this paper,starting from the basic principle of multi-component joint inversion,the diff erences between the method and single P-wave inversion are introduced.Next,the technique is applied to the PLN area of the Sichuan Basin,and the P-wave impedance,S-wave impedance,and density are obtained based on multi-component joint inversion.Through the velocity and lithology,porosity,and gas saturation fi tting formulas,prediction results are calculated,and the results are analyzed.Finally,multi-component joint inversion and single P-wave inversion are compared in eff ective reservoir prediction.The results show that multi-component joint inversion increases the constraints on the inversion conditions,reduces the multi-solution of a single P-wave inversion,and is more objective and reliable for the identification of reservoirs,effectively improving the accuracy of oil and gas reservoir prediction and development.
基金Support by Project of Geophysical Comprehensive Survey and Information Extraction of Deep Mineral Resources(2016YFC0600505)
文摘With the development of gravity gradient full tensor measurement technique,three-dimensional( 3D) inversion based on gravity gradient tensor can provide more accurate information. But the forward calculation of 3D full tensor sensitivity matrix is very time-consuming,which restricts its development and application.According to the symmetry of the kernel function,the authors reconstruct the underground source of geological body to avoid repeat computation of the same value,and work out the corresponding relationship between the response of geological body to the observation point and the response of reconstructed geological body to the observation point. According to the relationship,rapid calculation of full tensor gravity sensitivity matrix can be achieved. The model calculation shows that this method can increase the speed of 30-45 times compared with the traditional calculation method. The sensitivity matrix is applied to the multi-component inversion of gravity gradient. The application of this method on the measured data provides the basis for the promotion of the method.
基金supported by Fundamental Research Funds for the Central Universities,CHD300102264715National Key Research and Development Program of China under Grant 2021YFA0716902Natural Science Basic Research Program of Shaanxi 2024JCYBMS-199。
文摘Young's modulus and Poisson's ratio are crucial parameters for reservoir characterization and rock brittleness evaluation.Conventional methods often rely on indirect computation or approximations of the Zoeppritz equations to estimate Young's modulus,which can introduce cumulative errors and reduce the accuracy of inversion results.To address these issues,this paper introduces the analytical solution of the Zoeppritz equation into the inversion process.The equation is re-derived and expressed in terms of Young's modulus,Poisson's ratio,and density.Within the Bayesian framework,we construct an objective function for the joint inversion of PP and PS waves.Traditional gradient-based algorithms often suffer from low precision and the computational complexity.In this study,we address limitations of conventional approaches related to low precision and complicated code by using Circle chaotic mapping,Levy flights,and Gaussian mutation to optimize the quantum particle swarm optimization(QPSO),named improved quantum particle swarm optimization(IQPSO).The IQPSO demonstrates superior global optimization capabilities.We test the proposed inversion method with both synthetic and field data.The test results demonstrate the proposed method's feasibility and effectiveness,indicating an improvement in inversion accuracy over traditional methods.
基金funded by the National Natural Science Foundation of China under Grant 42230101the National Agency for Research and Development of Chile(ANID)by Projects AFB180004 and AFB220002the ANID Programa de Cooperación Internacional(PCI)Grant PII-180003.
文摘The Chilean Pampean flat slab subduction segment is characterized by the nearly horizontal subduction of the Nazca Plate within the depth range of 100-120 km.Numerous seismic tomography studies have been conducted to investigate its velocity structure;however,they have used only seismic body wave data or surface wave data.As a result,the existing velocity models in the region may have relatively large uncertainties.In this study,we use body wave arrival times from earthquakes occurring in central Chile between 2014 and 2019,as well as Rayleigh wave phase velocity maps at periods of 5-80 s from ambient noise empirical Green’s functions in Chile.By jointly using body wave arrival times and surface wave dispersion data,we refine the VS model and improve earthquake locations in the central Chile subduction zone.Compared with previous velocity models,our velocity model better reveals an eastward-dipping high-velocity plate representing the subducting Nazca Plate,which is 40-50 km thick and is more consistent with the slab thickness estimated by receiver function imaging and thermal modeling.Overall,the intraslab seismicity distribution spatially correlates well with the slab high-velocity anomalies except along the subduction paths of the CopiapóRidge and Juan Fernández Ridge.Additionally,parallel low-velocity stripes are imaged beneath the subducting plate,which are likely associated with the accumulated melts.The joint inversion velocity model also resolves widespread low-velocity anomalies in the crust beneath the Central Volcanic Zone of the central Andes,likely representing crustal magma chambers for various volcanoes.
基金supported by the China Important National Science & Technology Specific Projects (Grant No. 2011ZX05019-008)the National Natural Science Foundation of China (Grant No. 40839901)
文摘Based on the Aki-Richards approximate equations for reflection coefficients and Bayes theorem, we developed an inversion method to estimate P- and S-wave velocity contrasts and density contrast from combined PP and PS data. This method assumes that the parameters satisfy a normal distribution and introduces the covariance matrix to describe the degree of correlation between the parameters and thus to improve the inversion stability. Then, we suppose that the parameter sequence is subject to the Cauchy distribution and employs another matrix Q to describe the parameter sequence sparseness to improve the inversion result resolution. Tests on both synthetic and real multi-component data prove that this method is valid, efficient, more stable, and more accurate compared to methods using PP data only.
基金supported by the Natural Science Foundation of China (Grant Nos 40974066 and 40821062)the National Basic Research Program of China (Grant No. 2007CB209602)
文摘Multi-component exploration has many advantages over ordinary P-wave exploration. PP/PS joint AVO analysis and inversion are useful and powerful methods to discriminate between reservoir and non-productive lithology. In this paper, we derive a new PS-wave reflection coefficient approximation equation which is more accurate at larger incidence angles. The equation is simplified for small incidence angles, which makes AVO analysis clearer and easier for angles less than 30 degrees. Based on this approximation, a PP/PS joint inversion is introduced. A real data example shows that oil sands, brine sands and shales can be differentiated based on the P- to S-wave velocity ratio from the PP/PS joint inversion. Fluid factors and Poisson＇s ratio also indicate an anomaly in the target zone at the oil well location.
基金supported by the 863 Program of China(No.2013AA064201)
文摘Conventional joint PP-PS inversion is based on approximations of the Zoeppritz equations and assumes constant VP/VS;therefore,the inversion precision and stability cannot satisfy current exploration requirements.We propose a joint PP-PS inversion method based on the exact Zoeppritz equations that combines Bayesian statistics and generalized linear inversion.A forward model based on the exact Zoeppritz equations is built to minimize the error of the approximations in the large-angle data,the prior distribution of the model parameters is added as a regularization item to decrease the ill-posed nature of the inversion,low-frequency constraints are introduced to stabilize the low-frequency data and improve robustness,and a fast algorithm is used to solve the objective function while minimizing the computational load.The proposed method has superior antinoising properties and well reproduces real data.
基金Supported by the National Natural Science Foundation of China(Grant No.40674063)National Hi-tech Research and Development Program of China(863Program)(Grant No.2006AA09Z311)
文摘Based on the synchronous joint gravity and magnetic inversion of single interface by Pilkington and the need of revealing Cenozoic and crystalline basement thickness in the new round of oil-gas exploration, we propose a joint gravity and magnetic inversion methodfor two-layer models by concentrating on the relationship between the change of thicknessI and position of the middle layer and anomaly and discuss the effects of the key parameters. Model tests and application to field data show the validity of this method.
基金supported by the Fundamental Research Funds for the Central Universities of China(No.2652017438)the National Science and Technology Major Project of China(No.2016ZX05003-003)
文摘Multiwave seismic technology promotes the application of joint PP–PS amplitude versus offset (AVO) inversion;however conventional joint PP–PS AVO inversioan is linear based on approximations of the Zoeppritz equations for multiple iterations. Therefore the inversion results of P-wave, S-wave velocity and density exhibit low precision in the faroffset;thus, the joint PP–PS AVO inversion is nonlinear. Herein, we propose a nonlinear joint inversion method based on exact Zoeppritz equations that combines improved Bayesian inference and a least squares support vector machine (LSSVM) to solve the nonlinear inversion problem. The initial parameters of Bayesian inference are optimized via particle swarm optimization (PSO). In improved Bayesian inference, the optimal parameter of the LSSVM is obtained by maximizing the posterior probability of the hyperparameters, thus improving the learning and generalization abilities of LSSVM. Then, an optimal nonlinear LSSVM model that defi nes the relationship between seismic refl ection amplitude and elastic parameters is established to improve the precision of the joint PP–PS AVO inversion. Further, the nonlinear problem of joint inversion can be solved through a single training of the nonlinear inversion model. The results of the synthetic data suggest that the precision of the estimated parameters is higher than that obtained via Bayesian linear inversion with PP-wave data and via approximations of the Zoeppritz equations. In addition, results using synthetic data with added noise show that the proposed method has superior anti-noising properties. Real-world application shows the feasibility and superiority of the proposed method, as compared with Bayesian linear inversion.
基金supported by the National Earthquake Science Experiment in Sichuan and Yunnan Provinces of China(#2016 CESE 0201)National Natural Science Foundation of China(#41574034)China National Special Fund for Earthquake Scientific Research in Public Interest(#201508008)
文摘Accurate determination of seismic velocity of the crust is important for understanding regional tectonics and crustal evolution of the Earth. We propose a stepwise joint linearized inversion method using surface wave dispersion, Rayleigh wave ZH ratio （i.e., ellipticity）, and receiver function data to better resolve 1D crustal shear wave velocity （Vs） structure. Surface wave dispersion and Rayleigh wave ZH ratio data are more sensitive to absolute variations of shear wave speed at depths, but their sensi- tivity kernels to shear wave speeds are different and complimentary. However, receiver function data are more sensitive to sharp velocity contrast （e.g., due to the existence of crustal interfaces） and Vp/Vs ratios. The stepwise inversion method takes advantages of the complementary sensitivities of each dataset to better constrain the Vs model in the crust. We firstly invert surface wave dispersion and ZH ratio data to obtain a 1D smooth absolute vs model and then incorporate receiver function data in the joint inver- sion to obtain a finer Vs model with better constraints on interface structures. Through synthetic tests, Monte Carlo error analyses, and application to real data, we demonstrate that the proposed joint inversion method can resolve robust crustal Vs structures and with little initial model dependency.
基金financially supported by the National Natural Science Foundation of China (Nos. 90915012 and 41090291)the Research Project in Earthquake Science, CEA (No.201108002)
文摘Teleseismic and GPS data were jointly inverted for the rupture process of the 2011 Tohoku earthquake. The inversion results show that it is a bilateral rupture event with an average rupture velocity less than 2.0 km/s along the fault strike direction. The source rupture process consists of three sub-events, the first oc- curred near the hypocenter and the rest two ruptured along the up-dip direction and broke the sea bed, causing a maximum slip of about 30 m. The large-scale sea bed breakage may account for the tremendous tsunami disaster which resulted in most of the death and missing in this mega earthquake.
基金jointly sponsored by the Fundamental Research Funds for the Central Universitiesthe National Natural Science Foundation of China(No.41374078)
文摘A two-dimensional forward and backward algorithm for the controlled-source audio-frequency magnetotelluric （CSAMT） method is developed to invert data in the entire region （near, transition, and far） and deal with the effects of artificial sources. First, a regularization factor is introduced in the 2D magnetic inversion, and the magnetic susceptibility is updated in logarithmic form so that the inversion magnetic susceptibility is always positive. Second, the joint inversion of the CSAMT and magnetic methods is completed with the introduction of the cross gradient. By searching for the weight of the cross-gradient term in the objective function, the mutual influence between two different physical properties at different locations are avoided. Model tests show that the joint inversion based on cross-gradient theory offers better results than the single-method inversion. The 2D forward and inverse algorithm for CSAMT with source can effectively deal with artificial sources and ensures the reliability of the final joint inversion algorithm.
基金study was supported jointly by National Natural Science Foundation of China (Nos.40821160549 and 41074032)CAS Knowledge Innovation Program (No. KZCX2-YW-116-1)Joint Seismological Science Fundation of China(Nos.200808078 and 200708035)
文摘We use interferometric synthetic aperture radar （InSAR） and broadband seismic waveform data to estimate a source model of the 11th July, 2004 M W 6.2 Zhongba earthquake, Tibet of China. This event occurred within the seismically active zone of southwestern Tibetan Plateau where the east-west extension of the upper crust is observed. Because of limitations in one pair of InSAR data available, there are trade-offs among centroid depth, rupture area and amount of slip. Available seismic data tightly constrain the focal mechanism and centroid depth of the earthquake but not the horizontal location. Together, two complementary data sets can be used to identify the actual fault plane, better constrain the slip model and event location. We first use regional seismic waveform to estimate point source mechanism, then InSAR data is used to obtain better location. Finally, a joint inversion of teleseismic P-waves and InSAR data is performed to obtain a distributed model. Our preferred point source mechanism indicates a seismic moment of ～2.2×10 18 N·m （～M W 6.2）, a fault plane solution of 171° （342 ° ）/42 ° （48 ° ）/-83 ° （-97 ° ）, corresponding to strike/dip/rake, and a depth of 11 km. The fault plane with strike of 171 ？ and dip of 42° is identified as the ruptured fault with the aid of InSAR data. The preferred source model features compact area of slips between depth of 5–11 km and 10 km along strike with maximum slip amplitude of about 1.5 m.
基金funded by National 973 Basic Research Developments Program of China (No.2005CB422104)863 National High Technique Research Development Project of China (No.2007AA060505)National Natural Science Foundation of China (No.40839901)
文摘Considering Zoeppritz equations, reflections of PP and PS are only the function of ratios of density and velocity. So the inversion results will be the same if the ratios are the same but values of density, velocities of P- wave and S-wave are different without strict constraint. This paper makes efforts to explore nonlinear simultaneous PP and PS inversion with expectation to reduce the ambiguity of AVO analysis by utilizing the redundancy of multi-component AVO measurements. Accurate estimation of ratio parameters depends on independence of input data. There are only two independent AVO attributes for PP reflectivity （i.e. intercept and gradient） and two for PS reflectivity （i.e. pseudo-intercept and pseudo-gradient or extreme amplitude）, respectively. For individual PP and PS inversion, the values of least-squares objective function do not converge around a large neighborhood of chosen true model parameters. Fortunately for joint PP and PS inversion the values of the least-squares objective function show closed contours with single minima. Finally the power function fitting is used to provide a higher precision AVO attributes than traditional polynomial fitting. By using the four independent fitting attributes （two independent attributes for PP and PS respectively）, the inversion of four ratio parameters （velocities and densities） would be estimated with less errors than that in traditional method.
基金supported by the National Science and Technology Major Project(No.2016ZX05047-002-001)
文摘Joint PP–PS inversion offers better accuracy and resolution than conventional P-wave inversion. P-and S-wave elastic moduli determined through data inversions are key parameters for reservoir evaluation and fluid characterization. In this paper, starting with the exact Zoeppritz equation that relates P-and S-wave moduli, a coefficient that describes the reflections of P-and converted waves is established. This method effectively avoids error introduced by approximations or indirect calculations, thus improving the accuracy of the inversion results. Considering that the inversion problem is ill-posed and that the forward operator is nonlinear, prior constraints on the model parameters and modified low-frequency constraints are also introduced to the objective function to make the problem more tractable. This modified objective function is solved over many iterations to continuously optimize the background values of the velocity ratio, which increases the stability of the inversion process. Tests of various models show that the method effectively improves the accuracy and stability of extracting P and S-wave moduli from underdetermined data. This method can be applied to provide inferences for reservoir exploration and fluid extraction.
基金supported by the National Science and Technology Major Project(No.2016ZX05018006)the National Key Research Development Program(No.2016YFC0601104)the National Natural Science Foundation of China(No.41472136)
文摘The coal-bearing strata of the deep Upper Paleozoic in the GS Sag have high hydrocarbon potential. Because of the absence of seismic data, we use electromagnetic （MT） and gravity data jointly to delineate the distribution of deep targets based on well logging and geological data. First, a preliminary geological model is established by using three-dimensional （3D） MT inversion results. Second, using the formation density and gravity anomalies, the preliminary geological model is modified by interactive inversion of the gravity data. Then, we conduct MT-constrained inversion based on the modified model to obtain an optimal geological model until the deviations at all stations are minimized. Finally, the geological model and a seismic profile in the middle of the sag is analysed. We determine that the deep reflections of the seismic profile correspond to the Upper Paleozoic that reaches thickness up to 800 m. The processing of field data suggests that the joint MT-gravity modeling and constrained inversion can reduce the multiple solutions for single geophysical data and thus improve the recognition of deep formations. The MT-constrained inversion is consistent with the geological features in the seismic section. This suggests that the joint MT and gravity modeling and constrained inversion can be used to delineate deep targets in similar basins.
基金supported by the Key Foundation of Institute of Seismology,China Earthquake Administration( IS200916004)
文摘Yushu Ms7.1 earthquake occurred on the Ganzi-Yushu fault zone, across which we carried out a joint relative-gravity and seismic-reflection survey, and then performed a gravity inversion constrained by the seismic-reflection result. Based on the data of complete Bouguer gravity anomaly and seismic reflection, we obtained a layered interface structure in deep crust down to Moho. Our study showed that the inversion could reveal the interfaces of strata along the survey profile and the directions of regional faults in two-dimension. From the characteristics of the observed topography of the Moho basement, we tentatively confirmed that the uplift of eastern edge of Qinghai-Tibet plateau was caused by the subduetion of the Indian plate.
基金supported by the National Key Research and Development Project of China(No:2017YFC0602201)
文摘Joint inversion based on a correlation constraint utilizes a linear correlation function as a structural constraint.The linear correlation function contains a denominator,which may result in a singularity as the objective function is optimized,leading to an unstable inversion calculation.To improve the robustness of this calculation,this paper proposes a new method in which a sinusoidal correlation function is employed as the structural constraint for joint inversion instead of the conventional linear correlation function.This structural constraint does not contain a denominator,thereby preventing a singularity.Compared with the joint inversion method based on a cross-gradient constraint,the joint inversion method based on a sinusoidal correlation constraint exhibits good performance.An application to actual data demonstrates that this method can process real data.
基金supported by the National Science and Technology Major Project(No.2016ZX05003-003)
文摘Most current prestack AVA joint inversion methods are based on the exact Zoeppritz equation and its various approximations. However, these equations only reflect the relation between reflection coefficients, incidence angles, and elastic parameters on either side of the interface, which means that wave-propagation effects, such as spherical spreading, attenuation, transmission loss, multiples, and event mismatching of P-and S-waves, are not considered and cannot accurately describe the true propagation characteristics of seismic waves. Conventional AVA inversion methods require that these wave-propagation effects have been fully corrected or attenuated before inversion but these requirements can hardly be satisfied in practice. Using a one-dimensional(1 D) earth model, the reflectivity method can simulate the full wavefield response of seismic waves. Therefore, we propose a nonlinear multicomponent prestack AVA joint inversion method based on the vectorized reflectivity method, which uses a fast nondominated sorting genetic algorithm(NSGA II) to optimize the nonlinear multiobjective function to estimate multiple parameters, such as P-wave velocity, S-wave velocity, and density. This approach is robust because it can simultaneously cope with more than one objective function without introducing weight coefficients. Model tests prove the effectiveness of the proposed inversion method. Based on the inversion results, we find that the nonlinear prestack AVA joint inversion using the reflectivity method yields more accurate inversion results than the inversion by using the exact Zoeppritz equation when the wave-propagation effects of transmission loss and internal multiples are not completely corrected.
基金supported by the Natural National Science Foundation of China under grant number 41474039China National Special Fund for Earthquake Scientific Research in Public Interest under grant number 2016 CESE 0201+1 种基金Shanghai Committee of Science and Technology under grant number 14231202600the Fundamental Research Funds for the Central Universities under grant number WK2080000053
文摘A new 3D velocity model of the crust and upper mantle in the southeastern （SE） margin of the Tibetan plateau was obtained by joint inversion of body- and sur- face-wave data. For the body-wave data, we used 7190 events recorded by 102 stations in the SE margin of the Tibetan plateau. The surface-wave data consist of Rayleigh wave phase velocity dispersion curves obtained from ambient noise cross-correlation analysis recorded by a dense array in the SE margin of the Tibetan plateau. The joint inversion clearly improves the Vs model because it is constrained by both data types. The results show that at around 10 km depth there are two low-velocity anomalies embedded within three high-velocity bodies along the Longmenshan fault system. These high-velocity bodies correspond well with the Precambrian massifs, and the two located to the northeast of 2013 Ms 7.0 Lushan earthquake are associated with high fault slip areas during the 2008 Wenchuan earthquake. The aftershock gap between 2013 Lushan earthquake and 2008 Wenchuan earthquake is associated with low-velocity anomalies, which also acts as a barrier zone for ruptures of two earthquakes. Generally large earthquakes （M 〉 5） in the region occurring from 2008 to 2015 are located around the high-velocity zones, indicating that they may act as asperities for these large earthquakes. Joint inversion results also clearly show that there exist low-velocity or weak zones in the mid-lower crust, which are not evenly distributed beneath the SE margin of Tibetan plateau.
