This is a case study of the application of pre-stack inverted elastic parameters to tight-sand reservoir prediction. With the development of oil and gas exploration, pre-stack data and inversion results are increasing...This is a case study of the application of pre-stack inverted elastic parameters to tight-sand reservoir prediction. With the development of oil and gas exploration, pre-stack data and inversion results are increasingly used for production objectives. The pre-stack seismic property studies include not only amplitude verse offset （AVO） but also the characteristics of other elastic property changes. In this paper, we analyze the elastic property parameters characteristics of gas- and wet-sands using data from four gas-sand core types. We found that some special elastic property parameters or combinations can be used to identify gas sands from water saturated sand. Thus, we can do reservoir interpretation and description using different elastic property data from the pre-stack seismic inversion processing. The pre- stack inversion method is based on the simplified Aki-Richard linear equation. The initial model can be generated from well log data and seismic and geologic interpreted horizons in the study area. The input seismic data is angle gathers generated from the common reflection gathers used in pre-stack time or depth migration. The inversion results are elastic property parameters or their combinations. We use a field data example to examine which elastic property parameters or combinations of parameters can most easily discriminate gas sands from background geology and which are most sensitive to pore-fluid content. Comparing the inversion results to well data, we found that it is useful to predict gas reservoirs using λ, λρ, λ/μ, and K/μ properties, which indicate the gas characteristics in the study reservoir.展开更多
In this article,based on the acoustic measurements of core samples obtained from the low to medium porosity and permeability reservoirs in the WXS Depression,the densities and P and S wave velocities of these core sam...In this article,based on the acoustic measurements of core samples obtained from the low to medium porosity and permeability reservoirs in the WXS Depression,the densities and P and S wave velocities of these core samples were obtained.Then based on these data,a series of elastic parameters were computed.From the basic theory and previous pore fluid research results,we derived a new fluid identification factor（F）.Using the relative variations,Ag/w and Ao/w,of the elastic parameters between gas and water saturated samples and between oil and water saturated samples,λρ,σHSFIF,Kρ,λρ-2μρ,and F as quantitative indicators,we evaluate the sensitivity of the different fluid identification factors to identify reservoir fluids and validate the effects by crossplots.These confirm that the new fluid identification factor（F） is more sensitive for distinguishing oil and water than the traditional method and is more favorable for fliud identification in low to medium porosity and permeability reservoirs.展开更多
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.展开更多
Cauchy priori distribution-based Bayesian AVO reflectivity inversion may lead to sparse estimates that are sensitive to large reflectivities. For the inversion, the computation of the covariance matrix and regularized...Cauchy priori distribution-based Bayesian AVO reflectivity inversion may lead to sparse estimates that are sensitive to large reflectivities. For the inversion, the computation of the covariance matrix and regularized terms requires prior estimation of model parameters, which makes the iterative inversion weakly nonlinear. At the same time, the relations among the model parameters are assumed linear. Furthermore, the reflectivities, the results of the inversion, or the elastic parameters with cumulative error recovered by integrating reflectivities are not well suited for detecting hydrocarbons and fuids. In contrast, in Bayesian linear AVO inversion, the elastic parameters can be directly extracted from prestack seismic data without linear assumptions for the model parameters. Considering the advantages of the abovementioned methods, the Bayesian AVO reflectivity inversion process is modified and Cauchy distribution is explored as a prior probability distribution and the time-variant covariance is also considered. Finally, we propose a new method for the weakly nonlinear AVO waveform inversion. Furthermore, the linear assumptions are abandoned and elastic parameters, such as P-wave velocity, S-wave velocity, and density, can be directly recovered from seismic data especially for interfaces with large reflectivities. Numerical analysis demonstrates that all the elastic parameters can be estimated from prestack seismic data even when the signal-to-noise ratio of the seismic data is low.展开更多
Inversion of Young’s modulus,Poisson’s ratio and density from pre-stack seismic data has been proved to be feasible and effective.However,the existing methods do not take full advantage of the prior information.With...Inversion of Young’s modulus,Poisson’s ratio and density from pre-stack seismic data has been proved to be feasible and effective.However,the existing methods do not take full advantage of the prior information.Without considering the lateral continuity of the inversion results,these methods need to invert the reflectivity first.In this paper,we propose multi-gather simultaneous inversion for pre-stack seismic data.Meanwhile,the total variation（TV）regularization,L1 norm regularization and initial model constraint are used.In order to solve the objective function contains L1norm,TV norm and L2 norm,we develop an algorithm based on split Bregman iteration.The main advantages of our method are as follows：（1）The elastic parameters are calculated directly from objective function rather than from their reflectivity,therefore the stability and accuracy of the inversion process can be ensured.（2）The inversion results are more in accordance with the prior geological information.（3）The lateral continuity of the inversion results are improved.The proposed method is illustrated by theoretical model data and experimented with a 2-D field data.展开更多
Up to now, the primary method for studying critical porosity and porous media are experimental measurements and data analysis. There are few references on how to numerically calculate porosity at the critical point, p...Up to now, the primary method for studying critical porosity and porous media are experimental measurements and data analysis. There are few references on how to numerically calculate porosity at the critical point, pore fluid-related parameters, or framework-related parameters. So in this article, we provide a method for calculating these elastic parameters and use this method to analyze gas-bearing samples. We first derive three linear equations for numerical calculations. They are the equation of density p versus porosity Ф, density times the square of compressional wave velocity p Vp^2 versus porosity, and density times the square of shear wave velocity pVs^2 versus porosity. Here porosity is viewed as an independent variable and the other parameters are dependent variables. We elaborate on the calculation steps and provide some notes. Then we use our method to analyze gas-bearing sandstone samples. In the calculations, density and P- and S-velocities are input data and we calculate eleven relative parameters for porous fluid, framework, and critical point. In the end, by comparing our results with experiment measurements, we prove the viability of the method.展开更多
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.展开更多
For the complicated reservoir description of the GD oilfield, P-wave and S-wave elastic impedance inversion was carried out using pre-stack seismic data to accurately identify the lithology of the reservoir. The joint...For the complicated reservoir description of the GD oilfield, P-wave and S-wave elastic impedance inversion was carried out using pre-stack seismic data to accurately identify the lithology of the reservoir. The joint inversion was performed using three or more partial stacks to overcome the singularity of post-stack seismic inversion that can not satisfy the requirements of complex reservoir description and to avoid the instability of the inversion result caused by low signal-noise ratio in the pre-stack gather. The basic theory of prestack elastic impedance inversion is briefly described in this paper and, using real data of the GD oilfield, the key steps of angle gather wavelet extraction, horizon calibration, S-wave velocity prediction, and elastic parameter extraction were analyzed and studied. The comprehensive interpretation of multiple elastic parameters determined from log analysis is a key to improving the effect ofprestack seismic inversion.展开更多
When there exists anisotropy in underground media elastic parameters of the observed coordinate possibly do not coincide with that of the natural coordinate. According to the theory that the density of potential energ...When there exists anisotropy in underground media elastic parameters of the observed coordinate possibly do not coincide with that of the natural coordinate. According to the theory that the density of potential energy, dissipating energy is independent of the coordinate, the relationship of elastic parameters between two coordinates is derived for two-phase anisotropic media. Then, pseudospectral method to solve wave equations of two-phase anisotropic media is derived. At last, we use this method to simulate wave propagation in two-phase anisotropic media four types of waves are observed in the snapshots, i.e., fast P wave and slow P wave, fast S wave and slow S wave. Shear wave splitting, SV wave cusps and elastic wave reflection and transmission are also observed.展开更多
To study the influence of temperature and water content on ultrasonic wave velocity and to establish the relationship between ultrasonic wave velocity and frozen silty clay strength, ultrasonic tests were conducted to...To study the influence of temperature and water content on ultrasonic wave velocity and to establish the relationship between ultrasonic wave velocity and frozen silty clay strength, ultrasonic tests were conducted to frozen silty clay by using RSM-SY5（T） nonmetal supersonic test meter, and the tensile strength and compressive strength of silty clay were measured under various negative temperatures. Test and analysis results indicate that, ultrasonic wave velocity rapidly changes in the temperature range of-1 ℃ to -5 ℃. Ultrasonic wave velocity increased with an increase of water content until the water content reached the critical water content, while decreased with an increase of water content after the water content exceeded the critical water content. This study showed that there was strong positive correlation between the ul- trasonic wave velocity and the frozen soil strength. As ultrasonic wave velocity increased, either tensile strength or com- pressive strength increased. Based on the experimental data, the relationship between ultrasonic wave velocity and frozen silty clay strength was obtained through regression analysis. It was found that the ultrasonic test technique can be used to test frozen soils and lay the foundation for the determination of frozen soil strength.展开更多
The Ordovician carbonate rocks of the Yingshan formation in the Tarim Basin have a complex pore structure owing to diagenetic and secondary structures. Seismic elastic parameters（e.g., wave velocity） depend on poros...The Ordovician carbonate rocks of the Yingshan formation in the Tarim Basin have a complex pore structure owing to diagenetic and secondary structures. Seismic elastic parameters（e.g., wave velocity） depend on porosity and pore structure. We estimated the average specific surface, average pore-throat radius, pore roundness, and average aspect ratio of carbonate rocks from the Tazhong area. High P-wave velocity samples have small average specific surface, small average pore-throat radius, and large average aspect ratio. Differences in the pore structure of dense carbonate samples lead to fluid-related velocity variability. However, the relation between velocity dispersion and average specific surface, or the average aspect ratio, is not linear. For large or small average specific surface, the pore structure of the rock samples becomes uniform, which weakens squirt fl ow and minimizes the residuals of ultrasonic data and predictions with the Gassmann equation. When rigid dissolved（casting mold） pores coexist with less rigid microcracks, there are significant P-wave velocity differences between measurements and predictions.展开更多
To establish the relationship among reservoir characteristics and rock physical parameters,we construct the well-bore rock physical models firstly,considering the influence factors,such as mineral composition,shale co...To establish the relationship among reservoir characteristics and rock physical parameters,we construct the well-bore rock physical models firstly,considering the influence factors,such as mineral composition,shale content,porosity,fluid type and saturation.Then with analyzing the change rules of elastic parameters along with the above influence factors and the cross-plots among elastic parameters,the sensitive elastic parameters of tight sandstone reservoir are determined,and the rock physics template of sweet spot is constructed to guide pre-stack seismic inversion.The results show that velocity ratio and Poisson impedance are the most sensitive elastic parameters to indicate the lithologic and gas-bearing properties of sweet spot in tight sandstone reservoir.The high-quality sweet spot is characterized by the lower velocity ratio and Poisson impedance.Finally,the actual seismic data are selected to predict the sweet spots in tight sandstone gas reservoirs,so as to verify the validity of the rock physical simulation results.The significant consistency between the relative logging curves and inversion results in different wells implies that the utilization of well-bore rock physical simulation can guide the prediction of sweet spot in tight sandstone gas reservoirs.展开更多
Amplitude variations with offset or incident angle （AVO/AVA） inversion are typically combined with statistical methods, such as Bayesian inference or deterministic inversion. We propose a joint elastic inversion met...Amplitude variations with offset or incident angle （AVO/AVA） inversion are typically combined with statistical methods, such as Bayesian inference or deterministic inversion. We propose a joint elastic inversion method in the time and frequency domain based on Bayesian inversion theory to improve the resolution of the estimated P- and S-wave velocities and density. We initially construct the objective function using Bayesian inference by combining seismic data in the time and frequency domain. We use Cauchy and Gaussian probability distribution density functions to obtain the prior information for the model parameters and the likelihood function, respectively. We estimate the elastic parameters by solving the initial objective function with added model constraints to improve the inversion robustness. The results of the synthetic data suggest that the frequency spectra of the estimated parameters are wider than those obtained with conventional elastic inversion in the time domain. In addition, the proposed inversion approach offers stronger antinoising compared to the inversion approach in the frequency domain. Furthermore, results from synthetic examples with added Gaussian noise demonstrate the robustness of the proposed approach. From the real data, we infer that more model parameter details can be reproduced with the proposed joint elastic inversion.展开更多
Shale mechanical properties are important for shale gas production,but the magnitudes are difficult to estimate,standard size cores are hard to sample,and secondary interstice generation is inevitable.This paper propo...Shale mechanical properties are important for shale gas production,but the magnitudes are difficult to estimate,standard size cores are hard to sample,and secondary interstice generation is inevitable.This paper proposes a method for determining shale macroscale modulus,which is determined at a hierarchy of scales from the nano-to macro-scales.Microscale measurements are upscaled to estimate the corresponding magnitudes at the macroscale.A case study is conducted with Silurian shale samples,using the hierarchy scales,gridding nanoindentation,atomic force microscopy(AFM),mineral liberation analysis(MLA),X-ray diffraction(XRD),and uniaxial compression tests.The mineral compositions are analyzed using MLA and XRD,and the shale composition is described in terms of clay minerals,organic matter,and siliceous and carbonate contents.The variation in the Young’s modulus is analyzed based on the recorded indentation depth curves and modulus distributions.The nanoindentation and AFM results are upscaled to the centimeter scale through the Mori-Tanaka method.The upscaled results exhibit satisfactory fitting with the conventional uniaxial compression results,although the fitting of the upscaled AFM results is better than nanoindentation.The proposed approach can be applied to promptly and comprehensively predict the shale mechanical parameters during shale gas exploration.展开更多
Triaxial cyclic loading tests have been performed to assess the influence of plastic deformation on inelastic deformational properties of anisotropic argillite with bedding planes which is regarded as a kind of transv...Triaxial cyclic loading tests have been performed to assess the influence of plastic deformation on inelastic deformational properties of anisotropic argillite with bedding planes which is regarded as a kind of transversely isotropic media.Considering argillite's anisotropy and inelastic deformational properties,theoretical formulae for calculating oriented elastic parameters were deduced by the unloading curves,which can be better fitted for the description of its elasticity than loading curves.Test results indicate that with the growth of accumulated plastic,strain,the apparent elastic modulus of argillite decreases in a form of exponential decay function,whereas the apparent Poisson ratio increase in a form of power equation.A ratio of unloading recoverable strain to the total strain increment occurred during a loading cycle is defined to illustrate the characteristic relations between anisotropic coupled elasto-plastic deformation and plastic strain.It is significant to observe that high stress level and plastic history have an inhibiting effect on argillite anisotropy.展开更多
The mechanical, electronic and magnetic properties of non-magnetic MgTe and ferro-magnetic （FM） Mgo.75 TM025 Te （TM = Fe, Co, Ni） in the zinc-blende phase are studied by ab-initio calculations for the first time. ...The mechanical, electronic and magnetic properties of non-magnetic MgTe and ferro-magnetic （FM） Mgo.75 TM025 Te （TM = Fe, Co, Ni） in the zinc-blende phase are studied by ab-initio calculations for the first time. We use the generalized gradient approximation functional for computing the structural stability, and mechanical properties, while the modified Becke and Johnson local （spin） density approximation （mBJLDA） is utilized for determining the electronic and magnetic properties. By comparing the energies of non-magnetic and FM calculations, we find that the compounds are stable in the FM phase, which is confirmed by their structural stabilities in terms of enthalpy of formation. Detailed descriptions of elastic properties of Mgo.75TMo.25Te alloys in the FM phase are also presented. For electronic properties, the spin- polarized electronic band structures and density of states are computed, showing that these compounds are direct bandgap materials with strong hybridizations of TM 3d states and Te p states. Further, the ferromagnetism is discussed in terms of the Zener free electron model, RKKY model and double exchange model. The charge density contours in the （110） plane are calculated to study bonding properties. The spin exchange splitting and crystal field splitting energies are also calculated. The distribution of electron spin density is employed in computing the magnetic moments appearing at the magnetic sites （Fe, Co, Ni）, as well as at the non-magnetic sites （Mg, Te）. It is found that the p-d hybridization causes not only magnetic moments on the magnetic sites but also induces negligibly small magnetic moments at the non-magnetic sites.展开更多
Inversion methods for conventional reservoirs cannot be used in tight reservoirs because of small differences in wave impedance between tight reservoirs and surrounding rocks. Tight reservoirs in westem China are char...Inversion methods for conventional reservoirs cannot be used in tight reservoirs because of small differences in wave impedance between tight reservoirs and surrounding rocks. Tight reservoirs in westem China are characterized by strong heterogeneity; thus, it is difficult to predict favorable zones using conventional reservoir inversion methods. In this study, we propose an inversion method based on the Xu-White model. First, we modify the test parameters of the Xu-White model until model P- and S-wave velocities and density values agree with well-logging data. Then, we calculate P- and S-wave velocities, density, and Poisson ratio for different lithologies and oil saturation and construct different geological models for interbedded sand and shale. Subsequently, we use approximations to Zoeppritz equations to perform prestack forward modeling and analyze the response characteristics of prestack gathers for different physical properties, lithology, and oil saturation. We analyze the response of thirteen types of elastic parameters to porosity, clay content, and oil saturation. After the optimization of the real prestack gathers in the Z area of western China, we select an elastic parameter that is the most sensitive to lithology and oil saturation to predict the distribution of oil-producing reservoirs with high accuracy.展开更多
The differential cross sections of elastic and inelastic scattering of ~3He ions on ^（13）C and ^（14）C have been studied at an energy of 37.9 MeV with a double folding model based on M3Y-Reid effective nucleon-nucl...The differential cross sections of elastic and inelastic scattering of ~3He ions on ^（13）C and ^（14）C have been studied at an energy of 37.9 MeV with a double folding model based on M3Y-Reid effective nucleon-nucleon interaction.The resulted parameters have been used for the standard Distorted Wave Born Approximation calculations of angular distributions corresponding to different excitations levels of ^（13）C and ^（14）C and deformation parameters have been deduced.展开更多
There are ambiguities and uncertainties in the recognition of gas hydrate seismic reflections and in quantitative predictions of physical information of natural gas hydrate reservoirs from seismic data. Rock physical ...There are ambiguities and uncertainties in the recognition of gas hydrate seismic reflections and in quantitative predictions of physical information of natural gas hydrate reservoirs from seismic data. Rock physical modelling is a bridge that transforms the seismic information of geophysical observations into physical information, but traditional rock physics models lack descriptions of reservoir micro-structures and pore-filling materials. Considering the mineral compositions and pore microstructures of gas hydrates, we built rock physical models for load-bearing and pore-filling gas hydrate-bearing sediments,describe the mineral compositions, pore connectivity and pore shape using effective media theory, calculated the shear properties of pore-filling gas hydrates using Patchy saturation theory and Generalized Gassmann theory, and then revealed the quantitative relation between the elastic parameters and physical parameters for gas hydrate-bearing sediments. The numerical modelling results have shown that the ratios of P-wave and S-wave velocities decrease with hydrate saturation, the P-wave and S-wave velocities of load-bearing gas hydrate-bearing sediments are more sensitive to hydrate saturation, sensitivity is higher with narrower pores, and the ratios of the P-wave and S-wave velocities of pore-filling gas hydrate-bearing sediments are more sensitive to shear properties of hydrates at higher hydrate saturations. Theoretical analysis and practical application results showed that the rock physical models in this paper can be used to calculate the quantitative relation between macro elastic properties and micro physical properties of gas hydrate-bearing sediments, offer shear velocity information lacking in well logging, determine elastic parameters that have more effective indicating abilities, obtain physical parameters such as hydrate saturation and pore aspect ratios, and provide a theoretical basis and practical guidance for gas hydrate quantitative predictions.展开更多
Tight zones of the gas bearing Kangan and Dalan formations of the South Pars gas field contain a considerable amount of unswept gas due to their low porosity, low permeability and isolated pore types. The current stud...Tight zones of the gas bearing Kangan and Dalan formations of the South Pars gas field contain a considerable amount of unswept gas due to their low porosity, low permeability and isolated pore types. The current study, integrates core data, rock elastic properties and 3D seismic attributes to delineate fight and low-reservoir-quality zones of the South Pars gas field. In the first step, the dynamic reservoir geomechanical parameters were calculated based on empirical relationships from well log data. The log-derived elastic moduli were validated with the available laboratory measurements of core data. Cross plots between estimated porosity and elastic parameters based on Young＇s modulus indicate that low porosity zone coincide with high values of Young＇s module. The results were validated with petro- graphic studies of the available thin sections. The core samples with low porosity and permeability are correlated with strong rocks with tight matrix frameworks and high elastic values. Subsequently, rock elastic properties including Young＇s modulus and Poisson＇s ratio along with porosity were estimated by using neural networks from a collection of 3D post-stack seismic attributes, such as acoustic impedance （ALl）, instantaneous phase of AI and apparent polarity. Distinguishing low reservoir quality areas in pay zones with unswept gas is then facilitated by locating low porosity and high elastic modulus values. An- hydrite zones are identified and eliminated as non-pay zones due to their characterization of zero porosi- ty and high Young modulus values. The methodology described has applications for unconventional re- servoirs more generally, because it is able to distinguish low porosity and permeability zones that are po- tentially productive from those unprospective zones with negligible reservoir quality.展开更多
基金supported by the National Basic Priorities Program "973" Project (Grant No.2007CB209600)China Postdoctoral Science Foundation Funded Project
文摘This is a case study of the application of pre-stack inverted elastic parameters to tight-sand reservoir prediction. With the development of oil and gas exploration, pre-stack data and inversion results are increasingly used for production objectives. The pre-stack seismic property studies include not only amplitude verse offset （AVO） but also the characteristics of other elastic property changes. In this paper, we analyze the elastic property parameters characteristics of gas- and wet-sands using data from four gas-sand core types. We found that some special elastic property parameters or combinations can be used to identify gas sands from water saturated sand. Thus, we can do reservoir interpretation and description using different elastic property data from the pre-stack seismic inversion processing. The pre- stack inversion method is based on the simplified Aki-Richard linear equation. The initial model can be generated from well log data and seismic and geologic interpreted horizons in the study area. The input seismic data is angle gathers generated from the common reflection gathers used in pre-stack time or depth migration. The inversion results are elastic property parameters or their combinations. We use a field data example to examine which elastic property parameters or combinations of parameters can most easily discriminate gas sands from background geology and which are most sensitive to pore-fluid content. Comparing the inversion results to well data, we found that it is useful to predict gas reservoirs using λ, λρ, λ/μ, and K/μ properties, which indicate the gas characteristics in the study reservoir.
基金supported by the the Key Project of Chinese Ministry of Education (Grant No.109035)the National Natural Science Foundation Key Project (Grant No.40830423)Key Projects of Students Extra-curricular Science and Technology Research Program of Schlumberger (Grant No.SLBX0908)
文摘In this article,based on the acoustic measurements of core samples obtained from the low to medium porosity and permeability reservoirs in the WXS Depression,the densities and P and S wave velocities of these core samples were obtained.Then based on these data,a series of elastic parameters were computed.From the basic theory and previous pore fluid research results,we derived a new fluid identification factor（F）.Using the relative variations,Ag/w and Ao/w,of the elastic parameters between gas and water saturated samples and between oil and water saturated samples,λρ,σHSFIF,Kρ,λρ-2μρ,and F as quantitative indicators,we evaluate the sensitivity of the different fluid identification factors to identify reservoir fluids and validate the effects by crossplots.These confirm that the new fluid identification factor（F） is more sensitive for distinguishing oil and water than the traditional method and is more favorable for fliud identification in low to medium porosity and permeability reservoirs.
基金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 High-Tech Research and Development Program of China(863 Program)(No.2008AA093001)
文摘Cauchy priori distribution-based Bayesian AVO reflectivity inversion may lead to sparse estimates that are sensitive to large reflectivities. For the inversion, the computation of the covariance matrix and regularized terms requires prior estimation of model parameters, which makes the iterative inversion weakly nonlinear. At the same time, the relations among the model parameters are assumed linear. Furthermore, the reflectivities, the results of the inversion, or the elastic parameters with cumulative error recovered by integrating reflectivities are not well suited for detecting hydrocarbons and fuids. In contrast, in Bayesian linear AVO inversion, the elastic parameters can be directly extracted from prestack seismic data without linear assumptions for the model parameters. Considering the advantages of the abovementioned methods, the Bayesian AVO reflectivity inversion process is modified and Cauchy distribution is explored as a prior probability distribution and the time-variant covariance is also considered. Finally, we propose a new method for the weakly nonlinear AVO waveform inversion. Furthermore, the linear assumptions are abandoned and elastic parameters, such as P-wave velocity, S-wave velocity, and density, can be directly recovered from seismic data especially for interfaces with large reflectivities. Numerical analysis demonstrates that all the elastic parameters can be estimated from prestack seismic data even when the signal-to-noise ratio of the seismic data is low.
基金supported by the National Natural Science Foundation of China (Nos.61775030,61571096,41301460,61362018,and 41274127)the key projects of Hunan Provincial Department of Education (No.16A174)
文摘Inversion of Young’s modulus,Poisson’s ratio and density from pre-stack seismic data has been proved to be feasible and effective.However,the existing methods do not take full advantage of the prior information.Without considering the lateral continuity of the inversion results,these methods need to invert the reflectivity first.In this paper,we propose multi-gather simultaneous inversion for pre-stack seismic data.Meanwhile,the total variation（TV）regularization,L1 norm regularization and initial model constraint are used.In order to solve the objective function contains L1norm,TV norm and L2 norm,we develop an algorithm based on split Bregman iteration.The main advantages of our method are as follows：（1）The elastic parameters are calculated directly from objective function rather than from their reflectivity,therefore the stability and accuracy of the inversion process can be ensured.（2）The inversion results are more in accordance with the prior geological information.（3）The lateral continuity of the inversion results are improved.The proposed method is illustrated by theoretical model data and experimented with a 2-D field data.
基金supported by the National Natural Science Foundation of China (Grant No.40874052)the Key Laboratory of Geo-detection (China University of Geosciences,Beijing),Ministry of Education
文摘Up to now, the primary method for studying critical porosity and porous media are experimental measurements and data analysis. There are few references on how to numerically calculate porosity at the critical point, pore fluid-related parameters, or framework-related parameters. So in this article, we provide a method for calculating these elastic parameters and use this method to analyze gas-bearing samples. We first derive three linear equations for numerical calculations. They are the equation of density p versus porosity Ф, density times the square of compressional wave velocity p Vp^2 versus porosity, and density times the square of shear wave velocity pVs^2 versus porosity. Here porosity is viewed as an independent variable and the other parameters are dependent variables. We elaborate on the calculation steps and provide some notes. Then we use our method to analyze gas-bearing sandstone samples. In the calculations, density and P- and S-velocities are input data and we calculate eleven relative parameters for porous fluid, framework, and critical point. In the end, by comparing our results with experiment measurements, we prove the viability of the method.
文摘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.
文摘For the complicated reservoir description of the GD oilfield, P-wave and S-wave elastic impedance inversion was carried out using pre-stack seismic data to accurately identify the lithology of the reservoir. The joint inversion was performed using three or more partial stacks to overcome the singularity of post-stack seismic inversion that can not satisfy the requirements of complex reservoir description and to avoid the instability of the inversion result caused by low signal-noise ratio in the pre-stack gather. The basic theory of prestack elastic impedance inversion is briefly described in this paper and, using real data of the GD oilfield, the key steps of angle gather wavelet extraction, horizon calibration, S-wave velocity prediction, and elastic parameter extraction were analyzed and studied. The comprehensive interpretation of multiple elastic parameters determined from log analysis is a key to improving the effect ofprestack seismic inversion.
文摘When there exists anisotropy in underground media elastic parameters of the observed coordinate possibly do not coincide with that of the natural coordinate. According to the theory that the density of potential energy, dissipating energy is independent of the coordinate, the relationship of elastic parameters between two coordinates is derived for two-phase anisotropic media. Then, pseudospectral method to solve wave equations of two-phase anisotropic media is derived. At last, we use this method to simulate wave propagation in two-phase anisotropic media four types of waves are observed in the snapshots, i.e., fast P wave and slow P wave, fast S wave and slow S wave. Shear wave splitting, SV wave cusps and elastic wave reflection and transmission are also observed.
基金supported by the National Natural Science Foundation of China (No.41271080 and No.41230630)the Western Project Program of the Chinese Academy of Sciences(KZCX2-XB3-19)the open fund of Qinghai Research and Observation Base, Key Laboratory of Highway Construction and Maintenance Technology in Permafrost Region Ministry of Transport, PRC (2012-12-4)
文摘To study the influence of temperature and water content on ultrasonic wave velocity and to establish the relationship between ultrasonic wave velocity and frozen silty clay strength, ultrasonic tests were conducted to frozen silty clay by using RSM-SY5（T） nonmetal supersonic test meter, and the tensile strength and compressive strength of silty clay were measured under various negative temperatures. Test and analysis results indicate that, ultrasonic wave velocity rapidly changes in the temperature range of-1 ℃ to -5 ℃. Ultrasonic wave velocity increased with an increase of water content until the water content reached the critical water content, while decreased with an increase of water content after the water content exceeded the critical water content. This study showed that there was strong positive correlation between the ul- trasonic wave velocity and the frozen soil strength. As ultrasonic wave velocity increased, either tensile strength or com- pressive strength increased. Based on the experimental data, the relationship between ultrasonic wave velocity and frozen silty clay strength was obtained through regression analysis. It was found that the ultrasonic test technique can be used to test frozen soils and lay the foundation for the determination of frozen soil strength.
基金supported by the Natural Science Foundation of China(No.41274138)
文摘The Ordovician carbonate rocks of the Yingshan formation in the Tarim Basin have a complex pore structure owing to diagenetic and secondary structures. Seismic elastic parameters（e.g., wave velocity） depend on porosity and pore structure. We estimated the average specific surface, average pore-throat radius, pore roundness, and average aspect ratio of carbonate rocks from the Tazhong area. High P-wave velocity samples have small average specific surface, small average pore-throat radius, and large average aspect ratio. Differences in the pore structure of dense carbonate samples lead to fluid-related velocity variability. However, the relation between velocity dispersion and average specific surface, or the average aspect ratio, is not linear. For large or small average specific surface, the pore structure of the rock samples becomes uniform, which weakens squirt fl ow and minimizes the residuals of ultrasonic data and predictions with the Gassmann equation. When rigid dissolved（casting mold） pores coexist with less rigid microcracks, there are significant P-wave velocity differences between measurements and predictions.
基金supported by the National Key R&D Program of China(Grant No.2018YFC1405900)the Major Projects of National Science and Technology(Grant Nos.2016ZX05011-002,2016ZX05027-002-005)+3 种基金the National Natural Science Foundation of China(Grant No.41806073)the Natural Science Foundation of Shandong Province(Grant No.ZR2017BD014)Shandong Provincial Key Laboratory of Depositional Mineralization and Sedimentary Minerals,Shandong University of Science and Technology(Grant No.DMSM2017042)the Fundamental Research Funds for the Central Universities(Grant No.201964016)
文摘To establish the relationship among reservoir characteristics and rock physical parameters,we construct the well-bore rock physical models firstly,considering the influence factors,such as mineral composition,shale content,porosity,fluid type and saturation.Then with analyzing the change rules of elastic parameters along with the above influence factors and the cross-plots among elastic parameters,the sensitive elastic parameters of tight sandstone reservoir are determined,and the rock physics template of sweet spot is constructed to guide pre-stack seismic inversion.The results show that velocity ratio and Poisson impedance are the most sensitive elastic parameters to indicate the lithologic and gas-bearing properties of sweet spot in tight sandstone reservoir.The high-quality sweet spot is characterized by the lower velocity ratio and Poisson impedance.Finally,the actual seismic data are selected to predict the sweet spots in tight sandstone gas reservoirs,so as to verify the validity of the rock physical simulation results.The significant consistency between the relative logging curves and inversion results in different wells implies that the utilization of well-bore rock physical simulation can guide the prediction of sweet spot in tight sandstone gas reservoirs.
基金supported by the National Nature Science Foundation Project(Nos.41604101 and U1562215)the National Grand Project for Science and Technology(No.2016ZX05024-004)+2 种基金the Natural Science Foundation of Shandong(No.BS2014NJ005)Science Foundation from SINOPEC Key Laboratory of Geophysics(No.33550006-15-FW2099-0027)the Fundamental Research Funds for the Central Universities
文摘Amplitude variations with offset or incident angle （AVO/AVA） inversion are typically combined with statistical methods, such as Bayesian inference or deterministic inversion. We propose a joint elastic inversion method in the time and frequency domain based on Bayesian inversion theory to improve the resolution of the estimated P- and S-wave velocities and density. We initially construct the objective function using Bayesian inference by combining seismic data in the time and frequency domain. We use Cauchy and Gaussian probability distribution density functions to obtain the prior information for the model parameters and the likelihood function, respectively. We estimate the elastic parameters by solving the initial objective function with added model constraints to improve the inversion robustness. The results of the synthetic data suggest that the frequency spectra of the estimated parameters are wider than those obtained with conventional elastic inversion in the time domain. In addition, the proposed inversion approach offers stronger antinoising compared to the inversion approach in the frequency domain. Furthermore, results from synthetic examples with added Gaussian noise demonstrate the robustness of the proposed approach. From the real data, we infer that more model parameter details can be reproduced with the proposed joint elastic inversion.
基金supported by the National Natural Science Foundation of China(Grant No.42072194,U1910205)the Fundamental Research Funds for the Central Universities(800015Z1190,2021YJSDC02)
文摘Shale mechanical properties are important for shale gas production,but the magnitudes are difficult to estimate,standard size cores are hard to sample,and secondary interstice generation is inevitable.This paper proposes a method for determining shale macroscale modulus,which is determined at a hierarchy of scales from the nano-to macro-scales.Microscale measurements are upscaled to estimate the corresponding magnitudes at the macroscale.A case study is conducted with Silurian shale samples,using the hierarchy scales,gridding nanoindentation,atomic force microscopy(AFM),mineral liberation analysis(MLA),X-ray diffraction(XRD),and uniaxial compression tests.The mineral compositions are analyzed using MLA and XRD,and the shale composition is described in terms of clay minerals,organic matter,and siliceous and carbonate contents.The variation in the Young’s modulus is analyzed based on the recorded indentation depth curves and modulus distributions.The nanoindentation and AFM results are upscaled to the centimeter scale through the Mori-Tanaka method.The upscaled results exhibit satisfactory fitting with the conventional uniaxial compression results,although the fitting of the upscaled AFM results is better than nanoindentation.The proposed approach can be applied to promptly and comprehensively predict the shale mechanical parameters during shale gas exploration.
基金Program(2011CB710601) supported by National Basic Research Program of ChinaProject(50925933) supported by National Natural Science Foundation of China+1 种基金Project(2008BAB29B03) supported by National Key Technology Research and Development Program of ChinaProject(2010-122-011) supported by Guizhou Provincial Department of Transportation,China
文摘Triaxial cyclic loading tests have been performed to assess the influence of plastic deformation on inelastic deformational properties of anisotropic argillite with bedding planes which is regarded as a kind of transversely isotropic media.Considering argillite's anisotropy and inelastic deformational properties,theoretical formulae for calculating oriented elastic parameters were deduced by the unloading curves,which can be better fitted for the description of its elasticity than loading curves.Test results indicate that with the growth of accumulated plastic,strain,the apparent elastic modulus of argillite decreases in a form of exponential decay function,whereas the apparent Poisson ratio increase in a form of power equation.A ratio of unloading recoverable strain to the total strain increment occurred during a loading cycle is defined to illustrate the characteristic relations between anisotropic coupled elasto-plastic deformation and plastic strain.It is significant to observe that high stress level and plastic history have an inhibiting effect on argillite anisotropy.
基金the Deanship of Scientific Research at King Saud University for funding this Prolific Research Group (PRG-1436-26)
文摘The mechanical, electronic and magnetic properties of non-magnetic MgTe and ferro-magnetic （FM） Mgo.75 TM025 Te （TM = Fe, Co, Ni） in the zinc-blende phase are studied by ab-initio calculations for the first time. We use the generalized gradient approximation functional for computing the structural stability, and mechanical properties, while the modified Becke and Johnson local （spin） density approximation （mBJLDA） is utilized for determining the electronic and magnetic properties. By comparing the energies of non-magnetic and FM calculations, we find that the compounds are stable in the FM phase, which is confirmed by their structural stabilities in terms of enthalpy of formation. Detailed descriptions of elastic properties of Mgo.75TMo.25Te alloys in the FM phase are also presented. For electronic properties, the spin- polarized electronic band structures and density of states are computed, showing that these compounds are direct bandgap materials with strong hybridizations of TM 3d states and Te p states. Further, the ferromagnetism is discussed in terms of the Zener free electron model, RKKY model and double exchange model. The charge density contours in the （110） plane are calculated to study bonding properties. The spin exchange splitting and crystal field splitting energies are also calculated. The distribution of electron spin density is employed in computing the magnetic moments appearing at the magnetic sites （Fe, Co, Ni）, as well as at the non-magnetic sites （Mg, Te）. It is found that the p-d hybridization causes not only magnetic moments on the magnetic sites but also induces negligibly small magnetic moments at the non-magnetic sites.
基金sponsored by The Major Science and Technology Projects of China National Petroleum Corporation(No.2011E-0303)
文摘Inversion methods for conventional reservoirs cannot be used in tight reservoirs because of small differences in wave impedance between tight reservoirs and surrounding rocks. Tight reservoirs in westem China are characterized by strong heterogeneity; thus, it is difficult to predict favorable zones using conventional reservoir inversion methods. In this study, we propose an inversion method based on the Xu-White model. First, we modify the test parameters of the Xu-White model until model P- and S-wave velocities and density values agree with well-logging data. Then, we calculate P- and S-wave velocities, density, and Poisson ratio for different lithologies and oil saturation and construct different geological models for interbedded sand and shale. Subsequently, we use approximations to Zoeppritz equations to perform prestack forward modeling and analyze the response characteristics of prestack gathers for different physical properties, lithology, and oil saturation. We analyze the response of thirteen types of elastic parameters to porosity, clay content, and oil saturation. After the optimization of the real prestack gathers in the Z area of western China, we select an elastic parameter that is the most sensitive to lithology and oil saturation to predict the distribution of oil-producing reservoirs with high accuracy.
文摘The differential cross sections of elastic and inelastic scattering of ~3He ions on ^（13）C and ^（14）C have been studied at an energy of 37.9 MeV with a double folding model based on M3Y-Reid effective nucleon-nucleon interaction.The resulted parameters have been used for the standard Distorted Wave Born Approximation calculations of angular distributions corresponding to different excitations levels of ^（13）C and ^（14）C and deformation parameters have been deduced.
基金supported by the National Natural Science Foundation of China (Grant No. 41706042)the China Postdoctoral Science Foundation (Grant No. 2015M582060)+2 种基金the Special Fund for Land & Resources Scientific Research in the Public Interest (Grant No. 201511037)the National Key Research and Development Program (Grant No. 2017YFC0307400)the Foundation of Key Laboratory of Submarine Geosciences (Grant No. KLSG1603)
文摘There are ambiguities and uncertainties in the recognition of gas hydrate seismic reflections and in quantitative predictions of physical information of natural gas hydrate reservoirs from seismic data. Rock physical modelling is a bridge that transforms the seismic information of geophysical observations into physical information, but traditional rock physics models lack descriptions of reservoir micro-structures and pore-filling materials. Considering the mineral compositions and pore microstructures of gas hydrates, we built rock physical models for load-bearing and pore-filling gas hydrate-bearing sediments,describe the mineral compositions, pore connectivity and pore shape using effective media theory, calculated the shear properties of pore-filling gas hydrates using Patchy saturation theory and Generalized Gassmann theory, and then revealed the quantitative relation between the elastic parameters and physical parameters for gas hydrate-bearing sediments. The numerical modelling results have shown that the ratios of P-wave and S-wave velocities decrease with hydrate saturation, the P-wave and S-wave velocities of load-bearing gas hydrate-bearing sediments are more sensitive to hydrate saturation, sensitivity is higher with narrower pores, and the ratios of the P-wave and S-wave velocities of pore-filling gas hydrate-bearing sediments are more sensitive to shear properties of hydrates at higher hydrate saturations. Theoretical analysis and practical application results showed that the rock physical models in this paper can be used to calculate the quantitative relation between macro elastic properties and micro physical properties of gas hydrate-bearing sediments, offer shear velocity information lacking in well logging, determine elastic parameters that have more effective indicating abilities, obtain physical parameters such as hydrate saturation and pore aspect ratios, and provide a theoretical basis and practical guidance for gas hydrate quantitative predictions.
文摘Tight zones of the gas bearing Kangan and Dalan formations of the South Pars gas field contain a considerable amount of unswept gas due to their low porosity, low permeability and isolated pore types. The current study, integrates core data, rock elastic properties and 3D seismic attributes to delineate fight and low-reservoir-quality zones of the South Pars gas field. In the first step, the dynamic reservoir geomechanical parameters were calculated based on empirical relationships from well log data. The log-derived elastic moduli were validated with the available laboratory measurements of core data. Cross plots between estimated porosity and elastic parameters based on Young＇s modulus indicate that low porosity zone coincide with high values of Young＇s module. The results were validated with petro- graphic studies of the available thin sections. The core samples with low porosity and permeability are correlated with strong rocks with tight matrix frameworks and high elastic values. Subsequently, rock elastic properties including Young＇s modulus and Poisson＇s ratio along with porosity were estimated by using neural networks from a collection of 3D post-stack seismic attributes, such as acoustic impedance （ALl）, instantaneous phase of AI and apparent polarity. Distinguishing low reservoir quality areas in pay zones with unswept gas is then facilitated by locating low porosity and high elastic modulus values. An- hydrite zones are identified and eliminated as non-pay zones due to their characterization of zero porosi- ty and high Young modulus values. The methodology described has applications for unconventional re- servoirs more generally, because it is able to distinguish low porosity and permeability zones that are po- tentially productive from those unprospective zones with negligible reservoir quality.
