Research on anisotropy of shale oil reservoir based on rock physics model

Rock physics modeling is implemented for shales in the Luojia area of the Zhanhua topographic depression. In the rock physics model, the clay lamination parameter is introduced into the Backus averaging theory for the description of anisotropy related to the preferred alignment of clay particles, and the Chapman multi-scale fracture theory is used to calculate anisotropy relating to the fracture system. In accordance with geological features of shales in the study area, horizontal fractures are regarded as the dominant factor in the prediction of fracture density and anisotropy parameters for the inversion scheme. Results indicate that the horizontal fracture density obtained has good agreement with horizontal permeability measured from cores, and thus confirms the applicability of the proposed rock physics model and inversion method. Fracture density can thus be regarded as an indicator of reservoir permeability. In addition, the anisotropy parameter of the P-wave is higher than that of the S-wave due to the presence of horizontal fractures. Fracture density has an obvious positive correlation with P-wave anisotropy, and the clay content shows a positive correlation with S-wave anisotropy, which fully shows that fracture density has a negative correlation with clay and quartz contents and a positive relation with carbonate contents.

Anisotropy rock physics model for the Longmaxi shale gas reservoir,Sichuan Basin,China

The preferred orientation of clay minerals dominates the intrinsic anisotropy of shale. We introduce the clay lamination （CL） parameter to the Backus averaging method to describe the intrinsic shale anisotropy induced by the alignment of clay minerals. Then, we perform the inversion of CL and the Thomsen anisotropy parameters. The direct measurement of anisotropy is difficult because of the inability to measure the acoustic velocity in the vertical direction in boreholes and instrument limitations. By introducing the parameter CL, the inversion method provides reasonable estimates of the elastic anisotropy in the Longmaxi shale. The clay content is weakly correlated with the CL parameter. Moreover, the parameter CL is abnormally high at the bottom of the Longmaxi and Wufeng Formations, which are the target reservoirs. Finally, we construct rock physics templates to interpret well logging and reservoir properties.

Analysis of rock physics response of gas-bearing volcanic reservoir based on three-phase poroelastic theory

Unlike previous theories with velocity and/or elastic modulus averaging, we use a three-phase porous rock physics model developed by Santos for analyzing the seismic response of two immiscible fluids in saturated porous media. Considering reservoir reference pressure and coupling drag of two fluids in pores, the effects of frequency, porosity, and gas saturation on the phase velocities of the P-and S-waves are discussed in detail under field conditions. The effects of porosity and gas saturation on Vp/Vs are also provided. The data for our numerical experiments are from a sample of deep volcanic rock from Daqing. The numerical results show that the frequency dispersion effect can be ignored for deep volcanic rocks with low porosity and low permeability. It is concluded that for deep volcanic rocks the effect of gas content in pores on Vp/Vs is negligible but the effect of porosity is significant when there is a certain amount of water contained in the pores. The accurate estimate of lithology and porosity in this case is relatively more important.

Gas sand distribution prediction by prestack elastic inversion based on rock physics modeling and analysis

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.

Anisotropic rock physics models for interpreting pore structures in carbonate reservoirs

We developed an anisotropic effective theoretical model for modeling the elastic behavior of anisotropic carbonate reservoirs by combining the anisotropic self-consistent approximation and differential effective medium models.By analyzing the measured data from carbonate samples in the TL area,a carbonate pore-structure model for estimating the elastic parameters of carbonate rocks is proposed,which is a prerequisite in the analysis of carbonate reservoirs.A workflow for determining elastic properties of carbonate reservoirs is established in terms of the anisotropic effective theoretical model and the pore-structure model.We performed numerical experiments and compared the theoretical prediction and measured data.The result of the comparison suggests that the proposed anisotropic effective theoretical model can account for the relation between velocity and porosity in carbonate reservoirs.The model forms the basis for developing new tools for predicting and evaluating the properties of carbonate reservoirs.

Rock physics inversion based on an optimized MCMC method

Rock physics inversion is to use seismic elastic properties of underground strata for predicting reservoir petrophysical parameters.The Markov chain Monte Carlo(MCMC)algorithm is commonly used to solve rock physics inverse problems.However,all the parameters to be inverted are iterated simultaneously in the conventional MCMC algorithm.What is obtained is an optimal solution of combining the petrophysical parameters with being inverted.This study introduces the alternating direction(AD)method into the MCMC algorithm(i.e.the optimized MCMC algorithm)to ensure that each petrophysical parameter can get the optimal solution and improve the convergence of the inversion.Firstly,the Gassmann equations and Xu-White model are used to model shaly sandstone,and the theoretical relationship between seismic elastic properties and reservoir petrophysical parameters is established.Then,in the framework of Bayesian theory,the optimized MCMC algorithm is used to generate a Markov chain to obtain the optimal solution of each physical parameter to be inverted and obtain the maximum posterior density of the physical parameter.The proposed method is applied to actual logging and seismic data and the results show that the method can obtain more accurate porosity,saturation,and clay volume.

Brittleness index and seismic rock physics model for anisotropic tight-oil sandstone reservoirs

Brittleness analysis becomes important when looking for sweet spots in tightoil sandstone reservoirs. Hence, appropriate indices are required as accurate brittleness evaluation criteria. We construct a seismic rock physics model for tight-oil sandstone reservoirs with vertical fractures. Because of the complexities in lithology and pore structure and the anisotropic characteristics of tight-oil sandstone reservoirs, the proposed model is based on the solid components, pore connectivity, pore type, and fractures to better describe the sandstone reservoir microstructure. Using the model, we analyze the brittleness sensitivity of the elastic parameters in an anisotropic medium and establish a new brittleness index. We show the applicability of the proposed brittleness index for tight-oil sandstone reservoirs by considering the brittleness sensitivity, the rock physics response characteristics, and cross-plots. Compared with conventional brittleness indexes, the new brittleness index has high brittleness sensitivity and it is the highest in oil-bearing brittle zones with relatively high porosity. The results also suggest that the new brittleness index is much more sensitive to elastic properties variations, and thus can presumably better predict the brittleness characteristics of sweet spots in tight-oil sandstone reservoirs.

Prediction of brittleness based on anisotropic rock physics model for kerogen-rich shale

The construction of a shale rock physics model and the selection of an appropriate brittleness index （B/） are two significant steps that can influence the accuracy of brittleness prediction. On one hand, the existing models of kerogen-rich shale are controversial, so a reasonable rock physics model needs to be built. On the other hand, several types of equations already exist for predicting the BI whose feasibility needs to be carefully considered. This study constructed a kerogen-rich rock physics model by performing the self- consistent approximation and the differential effective medium theory to model intercoupled clay and kerogen mixtures. The feasibility of our model was confirmed by comparison with classical models, showing better accuracy. Templates were constructed based on our model to link physical properties and the BL Different equations for the BI had different sensitivities, making them suitable for different types of formations. Equations based on Young＇s Modulus were sensitive to variations in lithology, while those using Lame＇s Coefficients were sensitive to porosity and pore fluids. Physical information must be considered to improve brittleness prediction.

A rock physics model for the characterization of organic-rich shale from elastic properties

Kerogen content and kerogen porosity play a significant role in elastic properties of organic-rich shales. We construct a rock physics model for organic-rich shales to quantify the effect of kerogen content and kerogen porosity using the Kuster and Toks6z theory and the selfconsistent approximation method. Rock physics modeling results show that with the increase of kerogen content and kerogen-related porosity, the velocity and density of shales decrease, and the effect of kerogen porosity becomes more obvious only for higher kerogen content. We also find that the Poisson＇s ratio of the shale is not sensitive to kerogen porosity for the case of gas saturation. Finally, for the seismic reflection responses of an organic-rich shale layer, forward modeling results indicate the fifth type AVO re- sponses which correspond to a negative intercept and a positive gradient. The absolute values of intercept and gradient increase with kerogen content and kerogen porosity, and present predictable variations associated with velocities and density.

The construction of shale rock physics model and brittleness prediction for high-porosity shale gas-bearing reservoir

Due to the huge differences between the unconventional shale and conventional sand reservoirs in many aspects such as the types and the characteristics of minerals,matrix pores and fluids,the construction of shale rock physics model is significant for the exploration and development of shale reservoirs.To make a better characterization of shale gas-bearing reservoirs,we first propose a new but more suitable rock physics model to characterize the reservoirs.We then use a well A to demonstrate the feasibility and reliability of the proposed rock physics model of shale gas-bearing reservoirs.Moreover,we propose a new brittleness indicator for the high-porosity and organic-rich shale gas-bearing reservoirs.Based on the parameter analysis using the constructed rock physics model,we finally compare the new brittleness indicator with the commonly used Young’s modulus in the content of quartz and organic matter,the matrix porosity,and the types of filled fluids.We also propose a new shale brittleness index by integrating the proposed new brittleness indicator and the Poisson’s ratio.Tests on real data sets demonstrate that the new brittleness indicator and index are more sensitive than the commonly used Young’s modulus and brittleness index for the high-porosity and high-brittleness shale gas-bearing reservoirs.

Boundary-refl ected waves and ultrasonic coda waves in rock physics experiments

Ultrasonic coda waves are widely usea to stuay hign-trequency scattering, however, ultrasonic coda waves are strongly affected by interference from by boundary-reflected waves. To understand the effect of boundary-reflected waves, we performed ultrasonic experiments using aluminum and shale samples, and the rotating staggered-mesh finite-difference method to simulate the wavefield. We analyzed the wavefield characteristics at the different receiving points and the interference characteristics of the boundary-reflected waves with the ultrasonic coda wave, and the effect of sample geometry on the ultrasonic coda waves. The increase in the aspect ratio of the samples delays the interference effect of the laterally reflected waves and reduces the effect on the ultrasonic coda waves. The main waves interfering with the ultrasonic coda waves are laterally reflected PP-, PS-, PPP-, and PPS-waves. The scattering and attenuation of the high-frequency energy in actual rocks can weaken the interference of laterally reflected waves with the ultrasonic coda waves.

Construction of complex digital rock physics based on full convolution network

Digital rock physics(DRP)is a paramount technology to improve the economic benefits of oil and gas fields,devise more scientific oil and gas field development plans,and create digital oil and gas fields.Currently,a significant gap is present between DRP theory and practical applications.Conventional digital-core construction focuses only on simple cores,and the recognition and segmentation effect of fractures and pores of complex cores is poor.The identification of rock minerals is inaccurate,which leads to the difference between the digital and actual cores.To promote the application of DRP in developing oil and gas fields,based on the high-precision X-ray computed tomography scanning technology,the U-Net deep learning model of the full convolution neural network is used to segment the pores,fractures,and matrix from the complex rock core with natural fractures innovatively.Simultaneously,the distribution of rock minerals is divided,and the distribution of rock conditions is corrected by X-ray diffraction.A pore—fracture network model is established based on the equivalent radius,which lays the foundation for fluid seepage simulation.Finally,the accuracy of the established a digital core is verified by the porosity measured via nuclear magnetic resonance technology,which is of great significance to the development and application of DRP in oil and gas fields.

Rock physics and seismic reflectivity parameterization and amplitude variation with offsets inversion in terms of total organic carbon indicator

Total organic carbon (TOC) prediction with elastic parameter inversions has been widely used in the identification and evaluation of source rocks. However, the elastic parameters used to predict TOC are not only determined by TOC but also depend on the other physical properties of source rocks. Besides, the TOC prediction with the elastic parameters inversion is an indirect method based on the statistical relationship obtained from well logs and experiment data. Therefore, we propose a rock physics model and define a TOC indicator mainly affected by TOC to predict TOC directly. The proposed rock physics model makes the equivalent elastic moduli of source rocks parameterized by the TOC indicator. Combining the equivalent elastic moduli of source rocks and Gray’s approximation leads to a novel linearized approximation of the P-wave reflection coefficient incorporating the TOC indicator. Model examples illustrate that the novel reflectivity approximation well agrees with the exact Zoeppritz equation until incident angles reach 40°. Convoluting the novel P-wave reflection approximation with seismic wavelets as the forward solver, an AVO inversion method based on the Bayesian theory is proposed to invert the TOC indicator with seismic data. The synthetic examples and field tests validate the feasibility and stability of the proposed AVO inversion approach. Using the inversion results of the TOC indicator, TOC is directly and accurately estimated in the target area.

Acoustic-electrical properties and rock physics models for shale-oil formations:prediction of reservoir properties of interbedded sandstone and shale layers

In recent years,the Yanchang shale-oil formations of the Ordos Basin are rich in reserves with complex lithology and structure characteristics,low porosity and low permeability,and weak anomalies for oil and water discriminations,have been the key targets of unconventional oil/gas resource exploration and development in the relevant areas.The joint acoustic-electrical(AE)properties can be used to interpret reservoir lithology,mineralogy,pore structure,and fluid saturation.To conduct tests of thin section analysis,X-ray diff raction,and ultrasonic and electrical experiments at diff erent pressures and saturation degrees,cores from the shale-oil formations in the Q area of the basin are collected.The variations in AE properties with respect to clay content,porosity,pressure(microfracture),and saturation are analyzed.The experimental results indicate that the rock physics behaviors of sandstones with diff erent clay contents vary significantly.The AE properties of clean sandstones are basically dependent on the microfractures(pressure),while for muddy sandstones,the clay content is an important factor affecting the responses.The target reservoir consists of interbedded sandstone and shale layers.The AE equivalent medium equations and the Gurevich theory are applied to establish the joint models for the diff erent lithologies and simulate the variations in AE properties with respect to fluid type,pore structure,and mineral components.The three-dimensional joint templates of clean and muddy sandstones,as well as shale,are developed based on the elastic and electrical attributes and then calibrated using the experimental and well-log data.The reservoir properties are estimated with the templates and validated by the log data.The results indicate that the joint templates based on lithology characteristics can eff ectively characterize the properties of interbedded sandstone and shale layers.Furthermore,the combined application of AE data provides more beneficial information for the assessment of rock properties,leading to precise estimates that conform with the actual formation conditions.

Seismic rock physics modeling of fractures and fluids in a tight gas sandstone reservoir

Tight sandstone gas reservoirs have characteristics of low porosity and low permeability,complex pore structure.In the paper,we consider the pore and micro-fracture systems in tight sandstone reservoir and perform rock physics modeling of sandstone gas reservoirs based on Chapman.Then we analyze the effects of relaxation time,and difference fluid types on dispersion and attention,and seismic responses.We found that the magnitude of the relaxation time is related to the gas saturation and the viscosity coefficient of the fluid.The magnitude of the relaxation time affects the frequency range of the velocity dispersion and attenuation,and the type of fluid affects the magnitude of the velocity dispersion and attenuation.As the relaxation time decreases,the calculated amplitude deceases.And the amplitude and the waveform of the synthetic seismic record are related to different fluid types and formation thickness.

Simplifying the Integration of Petrophysics and Rock-Physics to Identify Hydrocarbon Bearing Rocks on Seismic

A considerable effort has been made in the literature for quality assurance (QA) and quality checking (QC) of the petrophysical log data for computation of reservoir rock property parameters. Well log data plays an integral role in building a rock physics model for quantitative interpretation (QI) work. A poor-quality rock physics model may lead to significant financial and HSSE implications by drilling wells in undesired locations. Historically, a variety of techniques have been used including histograms and cross plots for reviewing the feasibility of petrophysical logs for QI work. However, no attempt has ever been made to introduce a simplified workflow. This paper serves two-fold. It provides a simplified step by step approach for building a petrophysics/rock physics model. A case study has been presented to compare the synthetic seismogram generated from the simplified workflow with the actual seismic trace at well locations. Secondly, the paper shows how a few key cross plots and rock property parameters provide adequate information to validate petrophysical data, distinguish overburden and reservoir sections, and to help identify fluids and saturation trends within the reservoir sands. In the mentioned case study, the robustness of the simplified rock physics model has helped seismic data to successfully distinguish hydrocarbon bearing reservoir sands from non-reservoir shales.

Facies-constrained prestack seismic probabilistic inversion driven by rock physics

Seismic Rock physics plays a bridge role between the rock moduli and physical properties of the hydrocarbon reservoirs.Prestack seismic inversion is an important method for the quantitative characterization of elasticity,physical properties,lithology and fluid properties of subsurface reservoirs.In this paper,a high order approximation of rock physics model for clastic rocks is established and one seismic AVO reflection equation characterized by the high order approximation(Jacobian and Hessian matrix)of rock moduli is derived.Besides,the contribution of porosity,shale content and fluid saturation to AVO reflectivity is analyzed.The feasibility of the proposed AVO equation is discussed in the direct estimation of rock physical properties.On the basis of this,one probabilistic AVO inversion based on differential evolution-Markov chain Monte Carlo stochastic model is proposed on the premise that the model parameters obey Gaussian mixture probability prior model.The stochastic model has both the global optimization characteristics of the differential evolution algorithm and the uncertainty analysis ability of Markov chain Monte Carlo model.Through the cross parallel of multiple Markov chains,multiple stochastic solutions of the model parameters can be obtained simultaneously,and the posterior probability density distribution of the model parameters can be simulated effectively.The posterior mean is treated as the optimal solution of the model to be inverted.Besides,the variance and confidence interval are utilized to evaluate the uncertainties of the estimated results,so as to realize the simultaneous estimation of reservoir elasticity,physical properties,discrete lithofacies and dry rock skeleton.The validity of the proposed approach is verified by theoretical tests and one real application case in eastern China.

Azimuthally pre-stack seismic inversion for orthorhombic anisotropy driven by rock physics

Based on the long-wavelength approximation,a set of parallel vertical fractures embedded in periodic thin interbeds can be regarded as an equivalent orthorhombic medium. Rock physics is the basis for constructing the relationship between fracture parameters and seismic response. Seismic scattering is an effective way to inverse anisotropic parameters. In this study,we propose a reliable method for predicting the Thomsen's weak anisotropic parameters and fracture weaknesses in an orthorhombic fractured reservoir using azimuthal pre-stack seismic data. First, considering the influence of fluid substitution in mineral matrix, porosity, fractures and anisotropic rocks, we estimate the orthorhombic anisotropic stiffness coefficients by constructing an equivalent rock physics model for fractured rocks. Further, we predict the logging elastic parameters, Thomsen's weak parameters, and fracture weaknesses to provide the initial model constraints for the seismic inversion. Then, we derive the P-wave reflection coefficient equation for the inversion of Thomsen's weak anisotropic parameters and fracture weaknesses.Cauchy-sparse and smoothing-model constraint regularization taken into account in a Bayesian framework, we finally develop a method of amplitude variation with angles of incidence and azimuth(AVAZ) inversion for Thomsen's weak anisotropic parameters and fracture weaknesses, and the model parameters are estimated by using the nonlinear iteratively reweighted least squares(IRLS) strategy. Both synthetic and real examples show that the method can directly estimate the orthorhombic characteristic parameters from the azimuthally pre-stack seismic data, which provides a reliable seismic inversion method for predicting Thomsen's weak anisotropic parameters and fracture weaknesses.

Petrophysical parameters inversion for heavy oil reservoir based on a laboratory-calibrated frequency-variant rock-physics model

Heavy oil has high density and viscosity, and exhibits viscoelasticity. Gassmann's theory is not suitable for materials saturated with viscoelastic fluids. Directly applying such model leads to unreliable results for seismic inversion of heavy oil reservoir. To describe the viscoelastic behavior of heavy oil, we modeled the elastic properties of heavy oil with varying viscosity and frequency using the Cole-Cole-Maxwell (CCM) model. Then, we used a CCoherent Potential Approximation (CPA) instead of the Gassmann equations to account for the fluid effect, by extending the single-phase fluid condition to two-phase fluid (heavy oil and water) condition, so that partial saturation of heavy oil can be considered. This rock physics model establishes the relationship between the elastic modulus of reservoir rock and viscosity, frequency and saturation. The viscosity of the heavy oil and the elastic moduli and porosity of typical reservoir rock samples were measured in laboratory, which were used for calibration of the rock physics model. The well-calibrated frequency-variant CPA model was applied to the prediction of the P- and S-wave velocities in the seismic frequency range (1–100 Hz) and the inversion of petrophysical parameters for a heavy oil reservoir. The pre-stack inversion results of elastic parameters are improved compared with those results using the CPA model in the sonic logging frequency (∼10 kHz), or conventional rock physics model such as the Xu-Payne model. In addition, the inversion of the porosity of the reservoir was conducted with the simulated annealing method, and the result fits reasonably well with the logging curve and depicts the location of the heavy oil reservoir on the time slice. The application of the laboratory-calibrated CPA model provides better results with the velocity dispersion correction, suggesting the important role of accurate frequency dependent rock physics models in the seismic prediction of heavy oil reservoirs.

Quantitative characterization of tight gas sandstone reservoirs using seismic data via an integrated rock-physics-based framework

Seismic characterizing of tight gas sandstone (TGS) reservoirs is essential for identifying promising gas-bearing regions. However, exploring the petrophysical significance of seismic-inverted elastic properties is challenging due to the complex microstructures in TGSs. Meanwhile, interbedded structures of sandstone and mudstone intensify the difficulty in accurately extracting the crucial tight sandstone properties. An integrated rock-physics-based framework is proposed to estimate the reservoir quality of TGSs from seismic data. TGSs with complex pore structures are modeled using the double-porosity model, providing a practical tool to compute rock physics templates for reservoir parameter estimation. The VP/VS ratio is utilized to predict the cumulative thickness of the TGS reservoirs within the target range via the threshold value evaluated from wireline logs for lithology discrimination. This approach also facilitates better capturing the elastic properties of the TGSs for quantitative seismic interpretation. Total porosity is estimated from P-wave impedance using the correlation obtained based on wireline log analysis. After that, the three-dimensional rock-physics templates integrated with the estimated total porosity are constructed to interpret microfracture porosity and gas saturation from velocity ratio and bulk modulus. The integrated framework can optimally estimate the parameters dominating the reservoir quality. The results of the indicator proposed based on the obtained parameters are in good agreement with the gas productions and can be utilized to predict promising TGS reservoirs. Moreover, the results suggest that considering microfracture porosity allows a more accurate prediction of high-quality reservoirs, further validating the applicability of the proposed method in the studied region.