Rock physical characteristics of deep dolomite under complex geological conditions:A case study of 4th Member of Sinian Dengying Formation in the Sichuan Basin,China

The deep-ultra deep carbonate reservoir in China,commonly subjected to modification of multi-stage diagenesis,has extremely high heterogeneity.Conventional rock physics analysis cannot accurately identify the elastic responses of reservoir.Here,the rock physics properties of the dolomite from the 4th Member of the Sinian Dengying Formation are experimentally measured,and the change law of rock physics characteristics is investigated within the framework of the diagenetic processes by the analysis of the elastic and petrologic characteristics,pore structure,and sedimentary environments.The results show that the differentiated diagenesis results in different pore structure characteristics and microtexture characteristics of the rock.The microbial dolomite of the algal mound-grain beach facies is subjected to the contemporaneous microbial dolomitization and seepage-reflux dolomitization,penecontemporaneous selective dissolution,burial dolomitization,and hydrothermal dolomitization.The resultant crystalline dolomite is found with one main type of the dolomite crystal contact boundaries,and the dissolution pore is extensive development.The siliceous,muddy,and limy dolomite of the interbeach sea environment mainly experiences the weak capillary concentration dolomitization,intensive mechanical compaction-induced densification,and burial dolomitization.Such crystalline dolomite is observed with four types of contact boundaries,namely the dolomite contact,clay contact,quartz contact,and calcite contact boundaries,and porosity mostly attributed to residual primary inter-granular or crystalline pores.The samples with the same crystal boundary condition have consistent correlations between the compressional-and shear-wave velocities,and between the compressional-wave velocity and the velocity ratio.Additionally,the variation of the acoustic velocity with effective pressure and the intensity of pore-scale fluid-related dispersion are controlled by the differentiation of pore structure types of the samples.The varied effects of soft pores like micro-cracks on the compressional-and shearwave velocity causes considerable changes in the relationships between the compressional-and shearwave velocities,compressional-wave velocity and velocity ratio,and porosity and acoustic velocity.This research is an attempt to demonstrate a novel method for investigating the rock physics variation of rock during the geological process,and the obtained findings can provide the rock physics basis for seismic prediction of the characteristics of deep carbonate reservoirs.

Sweet spot prediction in tight sandstone reservoir based on well-bore rock physical simulation

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.

Establishment and application of the anisotropic shale-rock physical model in the observation coordinate system

No shale-rock physical model has been established in the observation coordinate system.To this end,this paper carried out anisotropic wave velocity tests on shale rock and compared the Thomsen,Daley,and Berryman solutions to characterize anisotropic acoustic wave velocity.Finally,the Daley solution was selected.Based on basic rock physical models,such as SCA and DEM methods,and combined with the Daley solution,an anisotropic shale-rock physical model was established in the observation coordinate system and applied in Well B1 in the Luzhou area,Sichuan Basin.Our research conclusions were as follows:1.for the samples from the same core,the P-wave velocities in three directions were in the order VP11>VP45>VP33,shear wave velocity VS11 was the largest,but VS33 and VS45 did not follow the law of Vs33>Vs45 for some samples;2.the Daley solution,which not only considers the accuracy requirements but also has a complete expression of P-,SV-,and SH-waves,is most suitable for characterization of anisotropic wave velocity in this study area;3.the rock physical model constructed in the observation coordinate system has high accuracy,in which the absolute value of the relative error of the P-wave slowness was between 0%and 5.05%(0.55%on average),and that of shear-wave slowness was between 0%and 6.05%(0.59%on average);4.the acoustic waves recorded in Well B1 in the observation coordinate system were very different from those in the constitutive coordinate system.The relative difference of the P-wave was between 6.76%and 30.84%(14.68%on average),and that of the S-wave was between 7.00%and 23.44%(13.99%on average).The acoustic slowness measured in the observation coordinate system,such as in a deviated well or a horizontal well section,must be converted to the constitutive coordinate system before it can be used in subsequent engineering applications;5.the anisotropic shale-rock physical model built in the observation coordinate system proposed in this paper can provide basic data and guidance for subsequent pore pressure prediction,geomechanical modeling,and fracturing stimulation design for deviated and horizontal wells.

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.

A seismic elastic moduli module for the measurements of low-frequency wave dispersion and attenuation of fluid-saturated rocks under different pressures

Knowledge about the seismic elastic modulus dispersion,and associated attenuation,in fluid-saturated rocks is essential for better interpretation of seismic observations taken as part of hydrocarbon identification and time-lapse seismic surveillance of both conventional and unconventional reservoir and overburden performances.A Seismic Elastic Moduli Module has been developed,based on the forced-oscillations method,to experimentally investigate the frequency dependence of Young's modulus and Poisson's ratio,as well as the inferred attenuation,of cylindrical samples under different confining pressure conditions.Calibration with three standard samples showed that the measured elastic moduli were consistent with the published data,indicating that the new apparatus can operate reliably over a wide frequency range of f∈[1-2000,10^(6)]Hz.The Young's modulus and Poisson's ratio of the shale and the tight sandstone samples were measured under axial stress oscillations to assess the frequency-and pressure-dependent effects.Under dry condition,both samples appear to be nearly frequency independent,with weak pressure dependence for the shale and significant pressure dependence for the sandstone.In particular,it was found that the tight sandstone with complex pore microstructure exhibited apparent dispersion and attenuation under brine or glycerin saturation conditions,the levels of which were strongly influenced by the increased effective pressure.In addition,the measured Young's moduli results were compared with the theoretical predictions from a scaled poroelastic model with a reasonably good agreement,revealing that the combined fluid flow mechanisms at both mesoscopic and microscopic scales possibly responsible for the measured dispersion.

A rock-physical modeling method for carbonate reservoirs at seismic scale

Strong heterogeneity and complex pore systems of carbonate reservoir rock make its rock physics model building and fluid substitution difficult and complex. However, rock physics models connect reservoir parameters with seismic parameters and fluid substitution is the most effective tool for reservoir prediction and quantitative characterization. On the basis of analyzing complex carbonate reservoir pore structures and heterogeneity at seismic scale, we use the gridding method to divide carbonate rock into homogeneous blocks with independent rock parameters and calculate the elastic moduli of dry rock units step by step using different rock physics models based on pore origin and structural feature. Then, the elastic moduli of rocks saturated with different fluids are obtained using fluid substitution based on different pore connectivity. Based on the calculated elastic moduli of rock units, the Hashin-Shtrikman-Walpole elastic boundary theory is adopted to calculate the carbonate elastic parameters at seismic scale. The calculation and analysis of carbonate models with different combinations of pore types demonstrate the effects of pore type on rock elastic parameters. The simulated result is consistent with our knowledge of real data.

Seismic rock physical modelling for gas hydrate-bearing sediments

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.

Numerical simulation of non-Archie electrophysical property of saturated rock with lattice Boltzmann method

The electrophysical property of saturated rocks is very important for reservoir identification and evaluation. In this paper, the lattice Boltzmann method （LBM） was used to study the electrophysical property of rock saturated with fluid because of its advantages over conventional numerical approaches in handling complex pore geometry and boundary conditions. The digital core model was constructed through the accumulation of matrix grains based on their radius distribution obtained by the measurements of core samples. The flow of electrical current through the core model saturated with oil and water was simulated on the mesoscopic scale to reveal the non-Archie relationship between resistivity index and water saturation （I-Sw）. The results from LBM simulation and laboratory measurements demonstrated that the I-Sw relation in the range of low water saturation was generally not a straight line in the log-log coordinates as described by the Archie equation. We thus developed a new equation based on numerical simulation and physical experiments. This new equation was used to fit the data from laboratory core measurements and previously published data. Determination of fluid saturation and reservoir evaluation could be significantly improved by using the new equation.

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.

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.

The constructing of pore structure factor in carbonate rocks and the inversion of reservoir parameters

With a more complex pore structure system compared with clastic rocks, carbonate rocks have not yet been well described by existing conventional rock physical models concerning the pore structure vagary as well as the influence on elastic rock properties. We start with a discussion and an analysis about carbonate rock pore structure utilizing rock slices. Then, given appropriate assumptions, we introduce a new approach to modeling carbonate rocks and construct a pore structure algorithm to identify pore structure mutation with a basis on the Gassmann equation and the Eshelby-Walsh ellipsoid inclusion crack theory. Finally, we compute a single well＇s porosity using this new approach with full wave log data and make a comparison with the predicted result of traditional method and simultaneously invert for reservoir parameters. The study results reveal that the rock pore structure can significantly influence the rocks＇ elastic properties and the predicted porosity error of the new modeling approach is merely 0.74%. Therefore, the approach we introduce can effectively decrease the predicted error of reservoir parameters.

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.

Calculations of rock matrix modulus based on a linear regression relation

The rock matrix bulk modulus or its inverse, the compressive coefficient, is an important input parameter for fluid substitution by the Biot-Gassmann equation in reservoir prediction. However, it is not easy to accurately estimate the bulk modulus by using conventional methods. In this paper, we present a new linear regression equation for calculating the parameter. In order to get this equation, we first derive a simplified Gassmann equation by using a reasonable assumption in which the compressive coefficient of the saturated pore fluid is much greater than the rock matrix, and, second, we use the Eshelby- Walsh relation to replace the equivalent modulus of a dry rock in the Gassmann equation. Results from the rock physics analysis of rock sample from a carbonate area show that rock matrix compressive coefficients calculated with water-saturated and dry rock samples using the linear regression method are very close （their error is less than 1%）. This means the new method is accurate and reliable.

Rock skeleton models and seismic porosity inversion

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.

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.

Effects of micrite microtextures on the elastic and petrophysical properties of carbonate reservoirs

Apparent differences in sedimentation and diagenesis exist between carbonate reservoirs in different areas and affect their petrophysical and elastic properties.To elucidate the relevant mechanism,we study and analyze the characteristics of rock microstructure and elastic properties of carbonates and their variation regularity using 89 carbonate samples from the different areas The results show that the overall variation regularities of the physical and elastic properties of the carbonate rocks are controlled by the microtextures of the microcrystalline calcite,whereas the traditional classification of rock-and pore-structures is no longer applicable.The micrite microtextures can be divided,with respect to their morphological features,into porous micrite,compact micrite,and tight micrite.As the micrites evolves from the first to the last type,crystal boundaries are observed with increasingly close coalescence,the micritic intercrystalline porosity and pore-throat radius gradually decrease;meanwhile,the rigidity of the calcite microcrystalline particle boundary and elastic homogeneity are enhanced.As a result,the seismic elastic characteristics,such as permeability and velocity of samples,show a general trend of decreasing with the increase of porosity.For low-porosity rock samples(φ<5%)dominated by tight micrite,the micritic pores have limited contributions to porosity and permeability and the micrite elastic properties are similar to those of the rock matrix.In such cases,the macroscopic physical and elastic properties are more susceptible to the formation of cracks and dissolution pores,but these features are controlled by the pore structure.The pore aspect ratio can be used as a good indication of pore types.The bulk modulus aspect ratio for dissolution pores is greater than 0.2,whereas that of the intergranular pores ranges from 0.1 to 0.2.The porous and compact micrites are observed to have a bulk modulus aspect ratio less than 0.1,whereas the ratio of the tight micrite approaches 0.2。

Permeability Estimation of Shale Oil Reservoir with Laboratory-derived Data: A Case Study of the Chang 7 Member in Ordos Basin

The shale oil reservoir within the Yanchang Formations of Ordos Basin harbors substantial oil and gas resources and has recently emerged as the primary focus of unconventional oil and gas exploration and development.Due to its complex pore and throat structure,pronounced heterogeneity,and tight reservoir characteristics,the techniques for conventional oil and gas exploration and production face challenges in comprehensive implementation,also indicating that as a vital parameter for evaluating the physical properties of a reservoir,permeability cannot be effectively estimated.This study selects 21 tight sandstone samples from the Q area within the shale oil formations of Ordos Basin.We systematically conduct the experiments to measure porosity,permeability,ultrasonic wave velocities,and resistivity at varying confining pressures.Results reveal that these measurements exhibit nonlinear changes in response to effective pressure.By using these experimental data and effective medium model,empirical relationships between P-and S-wave velocities,permeability and resistivity and effective pressure are established at logging and seismic scales.Furthermore,relationships between P-wave impedance and permeability,and resistivity and permeability are determined.A comparison between the predicted permeability and logging data demonstrates that the impedance–permeability relationship yields better results in contrast to those of resistivity–permeability relationship.These relationships are further applied to the seismic interpretation of shale oil reservoir in the target layer,enabling the permeability profile predictions based on inverse P-wave impedance.The predicted results are evaluated with actual production data,revealing a better agreement between predicted results and logging data and productivity.