The three parameters of P-wave velocity, S-wave velocity, and density have remarkable differences in conventional prestack inversion accuracy, so study of the consistency inversion of the ＂three parameters＂ is very ...The three parameters of P-wave velocity, S-wave velocity, and density have remarkable differences in conventional prestack inversion accuracy, so study of the consistency inversion of the ＂three parameters＂ is very important. In this paper, we present a new inversion algorithm and approach based on the in-depth analysis of the causes in their accuracy differences. With this new method, the inversion accuracy of the three parameters is improved synchronously by reasonable approximations and mutual constraint among the parameters. Theoretical model calculations and actual data applications with this method indicate that the three elastic parameters all have high inversion accuracy and maintain consistency, which also coincides with the theoretical model and actual data. This method has good application prospects.展开更多
By substituting rock skeleton modulus expressions into Gassmann approximate fluid equation, we obtain a seismic porosity inversion equation. However, conventional rock skeleton models and their expressions are quite d...By substituting rock skeleton modulus expressions into Gassmann approximate fluid equation, we obtain a seismic porosity inversion equation. However, conventional rock skeleton models and their expressions are quite different from each other, resuling in different seismic porosity inversion equations, potentially leading to difficulties in correctly applying them and evaluating their results. In response to this, a uniform relation with two adjusting parameters suitable for all rock skeleton models is established from an analysis and comparison of various conventional rock skeleton models and their expressions including the Eshelby-Walsh, Pride, Geertsma, Nur, Keys-Xu, and Krief models. By giving the two adjusting parameters specific values, different rock skeleton models with specific physical characteristics can be generated. This allows us to select the most appropriate rock skeleton model based on geological and geophysical conditions, and to develop more wise seismic porosity inversion. As an example of using this method for hydrocarbon prediction and fluid identification, we apply this improved porosity inversion, associated with rock physical data and well log data, to the ZJ basin. Research shows that the existence of an abundant hydrocarbon reservoir is dependent on a moderate porosity range, which means we can use the results of seismic porosity inversion to identify oil reservoirs and dry or water-saturated reservoirs. The seismic inversion results are closely correspond to well log porosity curves in the ZJ area, indicating that the uniform relations and inversion methods proposed in this paper are reliable and effective.展开更多
Seismic inversion is one of the most widely used technologies for reservoir prediction. Many good results have been obtained but sometimes it fails to differentiate the lithologies and identify the fluids. However, se...Seismic inversion is one of the most widely used technologies for reservoir prediction. Many good results have been obtained but sometimes it fails to differentiate the lithologies and identify the fluids. However, seismic prestack elastic inversion based on rock physics modeling and analysis introduced in this paper is a significant method that can help seismic inversion and interpretation reach a new quantitative （or semi-quantitative） level from traditional qualitative interpretation. By doing rock physics modeling and forward perturbation analysis, we can quantitatively analyze the essential relationships between rock properties and seismic responses and try to find the sensitive elastic properties to the lithology, porosity, fluid type, and reservoir saturation. Finally, standard rock physics templates （RPT） can be built for specific reservoirs to guide seismic inversion interpretation results for reservoir characterization and fluids identification purpose. The gas sand distribution results of the case study in this paper proves that this method has unparalleled advantages over traditional post-stack methods, by which we can perform reservoir characterization and seismic data interpretation more quantitatively and efficiently.展开更多
Stochastic seismic inversion is the combination of geostatistics and seismic inversion technology which integrates information from seismic records, well logs, and geostatistics into a posterior probability density fu...Stochastic seismic inversion is the combination of geostatistics and seismic inversion technology which integrates information from seismic records, well logs, and geostatistics into a posterior probability density function （PDF） of subsurface models. The Markov chain Monte Carlo （MCMC） method is used to sample the posterior PDF and the subsurface model characteristics can be inferred by analyzing a set of the posterior PDF samples. In this paper, we first introduce the stochastic seismic inversion theory, discuss and analyze the four key parameters： seismic data signal-to-noise ratio （S/N）, variogram, the posterior PDF sample number, and well density, and propose the optimum selection of these parameters. The analysis results show that seismic data S/N adjusts the compromise between the influence of the seismic data and geostatistics on the inversion results, the variogram controls the smoothness of the inversion results, the posterior PDF sample number determines the reliability of the statistical characteristics derived from the samples, and well density influences the inversion uncertainty. Finally, the comparison between the stochastic seismic inversion and the deterministic model based seismic inversion indicates that the stochastic seismic inversion can provide more reliable information of the subsurface character.展开更多
What are the anomalous seismic reflection bodies at depths of over 6000m？Are they reefs or igneous rock？This is a difficult problem for seismic techniques,but the GMES technique can handle it .The GMES technique is ...What are the anomalous seismic reflection bodies at depths of over 6000m？Are they reefs or igneous rock？This is a difficult problem for seismic techniques,but the GMES technique can handle it .The GMES technique is a joint exploration technique combining gravity,magnetic,electrical,and seismic techniques.The specific procedure is to conduct a 2D interface-constrained CEMP inversion using 2D seismic and log data followed by a property parameter inversion of the anomalous bodics using gravity and seismic data by the stripping technique.We then estimate the physical properties ofthe anomalous bodies,such as density,susceptibility,resistivity,velocity,and etc.to deduce the geological features of the bodies and provide a basis for drilling decisions.The work in the TZ area reported in this paper shows the applicability of the technique.展开更多
On April 20, 2013, a magnitude M s 7.0 earthquake occurred in Lushan, Sichuan Province, China, and caused heavy casualties and economic losses. Based on the local broadband waveforms in Sichuan and adjacent provinces ...On April 20, 2013, a magnitude M s 7.0 earthquake occurred in Lushan, Sichuan Province, China, and caused heavy casualties and economic losses. Based on the local broadband waveforms in Sichuan and adjacent provinces regional networks and teleseismic broadband records from IRIS stations, the focal mechanism and the focal depth are determined by the CAP (Cut And Paste) and its upgraded methods, CAPtele and CAPjoint, respectively. The results show that the focal mechanisms and depth from different methods are steady, and the best double couple solution derived from the joint inversion is 210°, 44°, and 91° for strike, dip, and rake angles respectively for one nodal plane and 29°, 46°, and 89° for another with 16 km focal depth and M w 6.66 moment magnitude. In order to verify the reliability of the results, a number of tests are performed based on local seismograms with different velocity models. They indicate that there is about 10 degree's fluctuation in focal mechanisms and about 2 km variation in focal depth with a complex velocity structure. Furthermore, inverted by re-sampling the teleseismic waveforms on the basis of epicentral distance, the solutions are consistent with each other, which manifests that the teleseismic records are effective for constraining source parameters of the Lushan earthquake.展开更多
Elastic wave inverse scattering theory plays an important role in parameters estimation of heterogeneous media.Combining inverse scattering theory,perturbation theory and stationary phase approximation,we derive the P...Elastic wave inverse scattering theory plays an important role in parameters estimation of heterogeneous media.Combining inverse scattering theory,perturbation theory and stationary phase approximation,we derive the P-wave seismic scattering coefficient equation in terms of fluid factor,shear modulus and density of background homogeneous media and perturbation media.With this equation as forward solver,a pre-stack seismic Bayesian inversion method is proposed to estimate the fluid factor of heterogeneous media.In this method,Cauchy distribution is utilized to the ratios of fluid factors,shear moduli and densities of perturbation media and background homogeneous media,respectively.Gaussian distribution is utilized to the likelihood function.The introduction of constraints from initial smooth models enhances the stability of the estimation of model parameters.Model test and real data example demonstrate that the proposed method is able to estimate the fluid factor of heterogeneous media from pre-stack seismic data directly and reasonably.展开更多
基金sponsored by the National Major Program (No. 2011ZX05006-006)the 973 Program of China (No. 2011CB201104)Technical Research of Elastic Flooding Boundary and Well Network Optimization at the Development Late Stage of Low Permeable Oil Field (No. 2011ZX05009)
文摘The three parameters of P-wave velocity, S-wave velocity, and density have remarkable differences in conventional prestack inversion accuracy, so study of the consistency inversion of the ＂three parameters＂ is very important. In this paper, we present a new inversion algorithm and approach based on the in-depth analysis of the causes in their accuracy differences. With this new method, the inversion accuracy of the three parameters is improved synchronously by reasonable approximations and mutual constraint among the parameters. Theoretical model calculations and actual data applications with this method indicate that the three elastic parameters all have high inversion accuracy and maintain consistency, which also coincides with the theoretical model and actual data. This method has good application prospects.
基金supported by the National Nature Science Foundation of China(Grant No.41174114)Important National Science and Technology Specific Projects(Grant No.2011ZX05025-005-010)
文摘By substituting rock skeleton modulus expressions into Gassmann approximate fluid equation, we obtain a seismic porosity inversion equation. However, conventional rock skeleton models and their expressions are quite different from each other, resuling in different seismic porosity inversion equations, potentially leading to difficulties in correctly applying them and evaluating their results. In response to this, a uniform relation with two adjusting parameters suitable for all rock skeleton models is established from an analysis and comparison of various conventional rock skeleton models and their expressions including the Eshelby-Walsh, Pride, Geertsma, Nur, Keys-Xu, and Krief models. By giving the two adjusting parameters specific values, different rock skeleton models with specific physical characteristics can be generated. This allows us to select the most appropriate rock skeleton model based on geological and geophysical conditions, and to develop more wise seismic porosity inversion. As an example of using this method for hydrocarbon prediction and fluid identification, we apply this improved porosity inversion, associated with rock physical data and well log data, to the ZJ basin. Research shows that the existence of an abundant hydrocarbon reservoir is dependent on a moderate porosity range, which means we can use the results of seismic porosity inversion to identify oil reservoirs and dry or water-saturated reservoirs. The seismic inversion results are closely correspond to well log porosity curves in the ZJ area, indicating that the uniform relations and inversion methods proposed in this paper are reliable and effective.
文摘Seismic inversion is one of the most widely used technologies for reservoir prediction. Many good results have been obtained but sometimes it fails to differentiate the lithologies and identify the fluids. However, seismic prestack elastic inversion based on rock physics modeling and analysis introduced in this paper is a significant method that can help seismic inversion and interpretation reach a new quantitative （or semi-quantitative） level from traditional qualitative interpretation. By doing rock physics modeling and forward perturbation analysis, we can quantitatively analyze the essential relationships between rock properties and seismic responses and try to find the sensitive elastic properties to the lithology, porosity, fluid type, and reservoir saturation. Finally, standard rock physics templates （RPT） can be built for specific reservoirs to guide seismic inversion interpretation results for reservoir characterization and fluids identification purpose. The gas sand distribution results of the case study in this paper proves that this method has unparalleled advantages over traditional post-stack methods, by which we can perform reservoir characterization and seismic data interpretation more quantitatively and efficiently.
基金supported by the National Major Science and Technology Project of China on Development of Big Oil-Gas Fields and Coalbed Methane (No. 2008ZX05010-002)
文摘Stochastic seismic inversion is the combination of geostatistics and seismic inversion technology which integrates information from seismic records, well logs, and geostatistics into a posterior probability density function （PDF） of subsurface models. The Markov chain Monte Carlo （MCMC） method is used to sample the posterior PDF and the subsurface model characteristics can be inferred by analyzing a set of the posterior PDF samples. In this paper, we first introduce the stochastic seismic inversion theory, discuss and analyze the four key parameters： seismic data signal-to-noise ratio （S/N）, variogram, the posterior PDF sample number, and well density, and propose the optimum selection of these parameters. The analysis results show that seismic data S/N adjusts the compromise between the influence of the seismic data and geostatistics on the inversion results, the variogram controls the smoothness of the inversion results, the posterior PDF sample number determines the reliability of the statistical characteristics derived from the samples, and well density influences the inversion uncertainty. Finally, the comparison between the stochastic seismic inversion and the deterministic model based seismic inversion indicates that the stochastic seismic inversion can provide more reliable information of the subsurface character.
文摘What are the anomalous seismic reflection bodies at depths of over 6000m？Are they reefs or igneous rock？This is a difficult problem for seismic techniques,but the GMES technique can handle it .The GMES technique is a joint exploration technique combining gravity,magnetic,electrical,and seismic techniques.The specific procedure is to conduct a 2D interface-constrained CEMP inversion using 2D seismic and log data followed by a property parameter inversion of the anomalous bodics using gravity and seismic data by the stripping technique.We then estimate the physical properties ofthe anomalous bodies,such as density,susceptibility,resistivity,velocity,and etc.to deduce the geological features of the bodies and provide a basis for drilling decisions.The work in the TZ area reported in this paper shows the applicability of the technique.
基金supported by National Natural Science Foundation of China (Grant Nos.41174086,41074052 and 41021003)Special Project Seismic Commonweal research (Grant No.201308013)Key Development Program of Chinese Academy of Science (Grant No.KZZD-EW-TZ-05)
文摘On April 20, 2013, a magnitude M s 7.0 earthquake occurred in Lushan, Sichuan Province, China, and caused heavy casualties and economic losses. Based on the local broadband waveforms in Sichuan and adjacent provinces regional networks and teleseismic broadband records from IRIS stations, the focal mechanism and the focal depth are determined by the CAP (Cut And Paste) and its upgraded methods, CAPtele and CAPjoint, respectively. The results show that the focal mechanisms and depth from different methods are steady, and the best double couple solution derived from the joint inversion is 210°, 44°, and 91° for strike, dip, and rake angles respectively for one nodal plane and 29°, 46°, and 89° for another with 16 km focal depth and M w 6.66 moment magnitude. In order to verify the reliability of the results, a number of tests are performed based on local seismograms with different velocity models. They indicate that there is about 10 degree's fluctuation in focal mechanisms and about 2 km variation in focal depth with a complex velocity structure. Furthermore, inverted by re-sampling the teleseismic waveforms on the basis of epicentral distance, the solutions are consistent with each other, which manifests that the teleseismic records are effective for constraining source parameters of the Lushan earthquake.
基金supported by the National Basic Research Program of China(Grant No.2013CB228604)the National Grand Project for Science and Technology(Grant Nos.2011ZX05030-004-002,2011ZX05019-003&2011ZX05006-002)
文摘Elastic wave inverse scattering theory plays an important role in parameters estimation of heterogeneous media.Combining inverse scattering theory,perturbation theory and stationary phase approximation,we derive the P-wave seismic scattering coefficient equation in terms of fluid factor,shear modulus and density of background homogeneous media and perturbation media.With this equation as forward solver,a pre-stack seismic Bayesian inversion method is proposed to estimate the fluid factor of heterogeneous media.In this method,Cauchy distribution is utilized to the ratios of fluid factors,shear moduli and densities of perturbation media and background homogeneous media,respectively.Gaussian distribution is utilized to the likelihood function.The introduction of constraints from initial smooth models enhances the stability of the estimation of model parameters.Model test and real data example demonstrate that the proposed method is able to estimate the fluid factor of heterogeneous media from pre-stack seismic data directly and reasonably.
