Microcrack Porosity Estimation Based on Rock Physics Templates: A Case Study in Sichuan Basin, China

: Low porosity-permeability structures and microcracks, where gas is produced, are the main characteristics of tight sandstone gas reservoirs in the Sichuan Basin, China. In this work, an analysis of amplitude variation with offset (AVO) is performed. Based on the experimental and log data, sensitivity analysis is performed to sort out the rock physics attributes sensitive to microcrack and total porosities. The Biot–Rayleigh poroelasticity theory describes the complexity of the rock and yields the seismic properties, such as Poisson’s ratio and P-wave impedance, which are used to build rock-physics templates calibrated with ultrasonic data at varying effective pressures. The templates are then applied to seismic data of the Xujiahe formation to estimate the total and microcrack porosities, indicating that the results are consistent with actual gas production reports.


Introduction
The development of society has led to a significant increase in the demand of oil and gas resources, and the depletion of conventional petroleum resources made the exploration and extraction of unconventional hydrocarbon resources necessary [1,2]. Tight sandstone reservoirs are widely distributed and account for a high production of China's total natural gas [3,4]. Large-scale gas fields with great potential have been discovered in more than ten basins, including Sichuan, Ordos, Songliao, Tuha, and Junggar [5].
Deep burial generally leads to high mechanical compaction and cementation of a sandstone reservoir, so that tight and heterogeneous reservoir rocks are generated, where microcracks are developed [6]. These microcracks affect the rock elastic properties and control the fluid flow, providing channels for hydrocarbon migration [7][8][9]. Their identification is a key factor for evaluating reservoir quality [10,11].
The study of microcracks in tight sandstones has become an important topic of rockphysics studies [12]. Hudson [13] established a relation between microcrack density and rock elastic properties, while Smith [14] stated that in low porosity rocks, the effect of cracks on seismic velocity can be important as mineral composition, total porosity, and fluid type. Yoon [15] established a relation between aspect ratio and microcrack density. Cheng [16] studied the effect of effective pressure and fluids. Zhang et al. [17,18] proposed a differential poroelastic model to describe wave propagation and dissipation in fluid-saturated rocks which contain inclusions at multiple scales.
The use of RPT (rock physics template; the list of abbreviations is given in Abbreviations) is widespread: Avseth and Ødegaard [19] linked the reservoir characteristics with elastic attributes, whereas Avseth et al. [20] used the Biot-Gassmann theory to discriminate between oil and gas in reservoirs. RPT was also used by Xin and Han [21] to estimate reservoir lithology and fluids. Michel [22] predicted porosity and oil saturation in shales, based on RPTs built with P-wave impedance (I P ) and the ratio between the P-and S-wave velocities (V P /V S ). Gupta et al. [23] identified thin sandy layers saturated with oil, and Ba et al. [24] estimated porosity and gas saturation. Carcione and Avseth [25] evaluated organic content, hydrocarbon saturation, and in-situ pressure in source rocks, and Liu et al. [26] proposed a combination of fluid indicators. Pang et al. [27] predicted microcrack properties in tight sandstone reservoirs, based on P-wave impedance and attenuation.
The Biot-Rayleigh (BR) theory was introduced by Ba et al. [28,29] for describing wave propagation characteristics in rocks with multi-phase pore structures. This work is mainly based on BR theory. We consider ten wells in the proposed area and perform a sensitivity analysis by using the ultrasonic and well-log data. On the basis of Poisson's ratio (ν) and P-wave impedance, we establish multi-scale RPTs for tight sandstone gas reservoirs, which relate the total and microcrack porosities to the seismic properties, which are obtained by means of seismic inversion.

Geology
The West Sichuan depression, bounded by the Longmenshan thrust and Longquanshan uplift belts, is a key area of natural gas exploration in China [30]. The maximum buried depth of the Xujiahe Formation exceeds 4.5 km, and the average thickness of the target layer is approximately 120 m. The work area was subjected to multi-stage Indosinian, Yanshan, and Himalayan tectonic movements, which resulted in diverse structural fractures that are controlled by complex diagenesis processes. A mechanical compaction is an aspect of diagenesis that destroys the primary pores. The sedimentary facies are mainly about the braided river delta facies. The channels overlap with each other from different periods, forming a longitudinal superposition and plane contiguous sand bodies [31].

Reservoir Characteristics
The tight sandstones are mainly composed of litharenites, with small amounts of lithic quartz arenite and feldspathic litharenite. In terms of mineral content, the average amounts of quartz and feldspar are 69% and 8.2%, respectively [32]. The effect of compaction on the reservoirs is more significant than that of cementation. Figure 1 shows that the compaction process is responsible for the developed geological characteristics such as complex pore structures and network microcracks. The average secondary, primary, and microcrack porosities are 2.9%, 1.1%, and 0.3%, respectively. The grain size varies from fine to medium, and the sorting property varies from poor to good. Moreover, the reservoir exhibits low porosity, low matrix permeability, and small pore-throat radii [31]. The microcracks maintain the reservoir connectivity and facilitate the gas production [33].

AVO Characteristics
Amplitude variations with offset (AVO) discriminate lithologies and help to detect hydrocarbons [34]. Rutherford and Williams [35] categorized the AVO responses of mudstone/gas-bearing sandstone interface into the three types. Subsequently, Castagna and Swan [36] added a fourth AVO response type. We analyze the AVO characteristics based on the profiles of Well P, shown in Figure 2, where the target layer is indicated with dashed red lines. Table 1 shows the seismic properties of the sandstone model of the Xujiahe Formation based on the log data. To obtain the P-and S-wave velocities and density, a 60 m interval is selected. The model data are collected from the interval above the top interface of the target layer that has a relatively stable P-wave velocity. For the target layer, we consider average values.

AVO Characteristics
Amplitude variations with offset (AVO) discriminate lithologies and help to detect hydrocarbons [34]. Rutherford and Williams [35] categorized the AVO responses of mudstone/gas-bearing sandstone interface into the three types. Subsequently, Castagna and Swan [36] added a fourth AVO response type. We analyze the AVO characteristics based on the profiles of Well P, shown in Figure 2, where the target layer is indicated with dashed red lines. Table 1 shows the seismic properties of the sandstone model of the Xujiahe Formation based on the log data. To obtain the P-and S-wave velocities and density, a 60 m interval is selected. The model data are collected from the interval above the top interface of the target layer that has a relatively stable P-wave velocity. For the target layer, we consider average values. Figure 3 illustrates the AVO of the top (red color) and bottom (blue color) interfaces. The characteristics of the curves agree with the fact that that the P-and S-wave impedances of the target layer are smaller than those of the overburden mudstone. The top and bottom responses are type IV and I AVO curves, respectively.    Figure 3 illustrates the AVO of the top (red color) and bottom (blue color) interfaces. The characteristics of the curves agree with the fact that that the P-and S-wave impedances of the target layer are smaller than those of the overburden mudstone. The top and bottom responses are type IV and I AVO curves, respectively.   Figure 4 shows an idealization of the tight sandstone of the Xujiahe Formation, which is characterized by a complex pore structure and microcrack network. The total porosity is the sum of intergranular and microcrack porosity.   Figure 4 shows an idealization of the tight sandstone of the Xujiahe Formation, which is characterized by a complex pore structure and microcrack network. The total porosity is the sum of intergranular and microcrack porosity.  Figure 4 shows an idealization of the tight sandstone of the Xujiahe Formation, which is characterized by a complex pore structure and microcrack network. The total porosity is the sum of intergranular and microcrack porosity.  The flowchart for the prediction of the reservoir properties is given in Figure 5. The microcrack porosity is estimated by using the quantitative relation between the elastic attributes and reservoir properties, based on a rock-physics model and aided by geological, log data, and ultrasonic data.

Theory and Flowchart
The flowchart for the prediction of the reservoir properties is given in Figure 5. The microcrack porosity is estimated by using the quantitative relation between the elastic attributes and reservoir properties, based on a rock-physics model and aided by geological, log data, and ultrasonic data. The main minerals are quartz, feldspar, and clay. The total porosity is less than 12%. The procedure to establish the rock-physics model is as follows.
(1) The Voigt-Reuss-Hill equation [37][38][39] is used to compute the elastic modulus of mineral mixture VRH M according to the mineral composition: The main minerals are quartz, feldspar, and clay. The total porosity is less than 12%. The procedure to establish the rock-physics model is as follows.
(1) The Voigt-Reuss-Hill equation [37][38][39] is used to compute the elastic modulus of mineral mixture M VRH according to the mineral composition: where f i and M i denote the volume fraction and elastic modulus, respectively, of the i-th component.
(2) According to the pore structure shown in Figure 4, we use the differential equivalent medium (DEM) theory [40] to add spherical pores and oblate cracks, whose aspect ratio are 1 and 0.0005, respectively, to add pores and microcracks into the matrix, and obtain the bulk and shear moduli of the rock skeleton (starred quantities), The initial conditions are K * (0) = K 1 and µ * (0) = µ 1 , where K 1 and µ 1 are the bulk and shear moduli, respectively, of the initial main material (phase 1); K 2 and µ 2 are the bulk and shear moduli, respectively, of the inclusion that is gradually added (phase 2); y is the volume content of phase 2; and P ( * 2) and Q ( * 2) [41] are related to the shape of the inclusion.
(3) The equations proposed by Batzle and Wang [42] are used to compute the bulk modulus, viscosity, and density of pore fluids at different temperatures and pressures.
(4) The BR theory is used to compute the wave response. Appendix A shows the computation of the complex wave number k. The complex velocity is v = ω/k, where ω is the angular frequency, and the P-wave velocity and quality factor are [7] where "Re" and "Im" take real and imaginary parts, respectively.
We consider the inclusion radius of 100 µm, and the bulk and shear moduli of the matrix are 33 GPa and 45 GPa, respectively. Gas has a bulk modulus of 0.02 GPa, a density of 0.089 g/cm 3 , and viscosity of 0.016 × 10 −3 Pa·s. The crack porosity is set to 0.2%, and the total porosity varies as shown in Figure 6, where the P-wave velocity and dissipation factor are plotted as a function of frequency. Increasing porosity implies increasing velocity dispersion, attenuation, and relaxation frequency.
where "Re" and "Im" take real and imaginary parts, respectively.
We consider the inclusion radius of 100 μm, and the bulk and shear moduli of the matrix are 33 GPa and 45 GPa, respectively. Gas has a bulk modulus of 0.02 GPa, a density of 0.089 g/cm 3 , and viscosity of 0.016 × 10 −3 Pa·s. The crack porosity is set to 0.2%, and the total porosity varies as shown in Figure 6, where the P-wave velocity and dissipation factor are plotted as a function of frequency. Increasing porosity implies increasing velocity dispersion, attenuation, and relaxation frequency.  Figure 7 shows the P-and S-wave velocities as a function of the total porosity for various microcrack porosities at 1 MHz. As can be seen, the velocities decrease with increasing total and microcrack porosities, as expected.  Figure 7 shows the P-and S-wave velocities as a function of the total porosity for various microcrack porosities at 1 MHz. As can be seen, the velocities decrease with increasing total and microcrack porosities, as expected.

Ultrasonic Experiments and Sensitivity Analysis
To investigate the effects of the microcracks, ultrasonic experiments at 1 MHz were performed. We select a sample with a porosity of 4.39% and a grain bulk modulus of 39 MPa. The grain, dry-rock, and wet-rock densities are 2.691, 2.573, and 2.62 g/cm 3 , respectively, where the sample is saturated with water. The experimental setup proposed by Guo et al. [43] is used to measure the velocities at 20 °C, with the ultrasonic pulse method. The sample is sealed with a rubber sleeve and placed in the vessel. The pore pressure and temperature are fixed. Effective pressures of 5, 10, 15, 20, 25, 30, and 35 MPa are applied to the sample, in both the gas-saturated and water-saturated cases, and the velocities are measured.
As shown in Figure 8a, the measured P-and S-wave velocities increase with effective pressure. The increase of the confining pressure at a fixed pore pressure will lead to the gradual closure of internal microcracks (especially at low effective pressures), so as to stiffen the rock skeleton and increase its elastic moduli (and the wave velocities).

Ultrasonic Experiments and Sensitivity Analysis
To investigate the effects of the microcracks, ultrasonic experiments at 1 MHz were performed. We select a sample with a porosity of 4.39% and a grain bulk modulus of 39 MPa. The grain, dry-rock, and wet-rock densities are 2.691, 2.573, and 2.62 g/cm 3 , respectively, where the sample is saturated with water. The experimental setup proposed by Guo et al. [43] is used to measure the velocities at 20 • C, with the ultrasonic pulse method. The sample is sealed with a rubber sleeve and placed in the vessel. The pore pressure and temperature are fixed. Effective pressures of 5, 10, 15, 20, 25, 30, and 35 MPa are applied to the sample, in both the gas-saturated and water-saturated cases, and the velocities are measured.
As shown in Figure 8a, the measured P-and S-wave velocities increase with effective pressure. The increase of the confining pressure at a fixed pore pressure will lead to the gradual closure of internal microcracks (especially at low effective pressures), so as to stiffen the rock skeleton and increase its elastic moduli (and the wave velocities). temperature are fixed. Effective pressures of 5, 10, 15, 20, 25, 30, and 35 MPa are applied to the sample, in both the gas-saturated and water-saturated cases, and the velocities are measured.
As shown in Figure 8a, the measured P-and S-wave velocities increase with effective pressure. The increase of the confining pressure at a fixed pore pressure will lead to the gradual closure of internal microcracks (especially at low effective pressures), so as to stiffen the rock skeleton and increase its elastic moduli (and the wave velocities).  Figure 8b shows the velocity ratio. In the full gas saturation case, P S V V increases with increasing effective pressure and the opposite behavior occurs for full water saturation case.
In addition, we measured the total porosity in the range 5-35MPa. Microcracks with small aspect ratios tend to close first with increasing effective pressure and the relation between porosity and pressure changes from exponential to linear [44].  Figure 8b shows the velocity ratio. In the full gas saturation case, V P /V S increases with increasing effective pressure and the opposite behavior occurs for full water saturation case.
In addition, we measured the total porosity in the range 5-35 MPa. Microcracks with small aspect ratios tend to close first with increasing effective pressure and the relation between porosity and pressure changes from exponential to linear [44].
Experimental measurements can be performed to predict microcrack porosity based on the relation between porosity and effective pressure [45][46][47]. Stiff porosity can be obtained through a linear extrapolation of this relation and the microcrack porosity can be estimated as the difference between the total and stiff porosities. Figure 9 shows the different porosities as a function of the effective pressure. Experimental measurements can be performed to predict microcrack porosity based on the relation between porosity and effective pressure [45][46][47]. Stiff porosity can be obtained through a linear extrapolation of this relation and the microcrack porosity can be estimated as the difference between the total and stiff porosities. Figure 9 shows the different porosities as a function of the effective pressure. The sensitivity analysis method like Sobol's indices can be used to analyze the interactions between the parameters, which can be applied with the given forward modeling equation or method/procedure [48]. In this work, sensitivity analysis to the total and mi- Figure 9. Total, stiff, and microcrack porosities as a function of the effective pressure. The sensitivity analysis method like Sobol's indices can be used to analyze the interactions between the parameters, which can be applied with the given forward modeling equation or method/procedure [48]. In this work, sensitivity analysis to the total and microcrack porosities as inputs is performed based on the actually measured data from the experiments and well logs as puts, without considering the detailed modeling equations (or procedures). Based on the observed experimental or log data, the fluid sensitivity indicator (FSI) [49][50][51][52] has been proposed to analyze how the rock elastic properties are affected by the different fluid saturation statuses. Whereas, in this study, a similar method is adopted to analyze the relative variations of rock physics attributes with respect to the changes in total/microcrack porosity, and it is defined as where A m and A represent the value of rock physics attribute at the minimal total (or microcrack) porosity and the average value of rock physics attribute within the considered range, respectively. Eleven attributes are considered, including density and the P-and S-wave velocities. The other eight are V P /V S , P-wave impedance (I P ), S-wave impedance (I S ), Poisson's ratio (ν), shear modulus (µ), the first Lamé constant (λ), λρ, and V P /ρ. Figure 10 shows the sensitivity indices to microcrack porosity at ultrasonic frequencies, where we observe that λρ, λ, ν, µ, and I P are the most sensitive ones.   The sensitivity analysis to total porosity is performed for the sonic log data, whe we consider the gas-saturated layers of Well P. Figure 11 shows the results, where λρ λ , ν , μ , and P I are the most sensitive. In Figure 10, P V ρ , P V , and P I have similar value while that of P I is slightly higher than the values of P V ρ and V However, it is obvious that in Figure 11 the value of P I is higher than the values P V ρ and P V . According to Figures 10 and 11, ν and P I are sensitive to both p rosities and can be considered to build the rock-physics templates. The sensitivity analysis to total porosity is performed for the sonic log data, where we consider the gas-saturated layers of Well P. Figure 11 shows the results, where λρ, λ, ν, µ, and I P are the most sensitive. In Figure 10, V P /ρ, V P , and I P have a similar value while that of I P is slightly higher than the values of V P /ρ and V P . However, it is obvious that in Figure 11 the value of I P is higher than the values of V P /ρ and V P . According to Figures 10 and 11, ν and I P are sensitive to both porosities and can be considered to build the rock-physics templates. λ , ν , μ , and P I are the most sensitive. In Figure 10, P V ρ , P V , and P I have similar value while that of P I is slightly higher than the values of P V ρ and V However, it is obvious that in Figure 11 the value of P I is higher than the values P V ρ and P V . According to Figures 10 and 11, ν and P I are sensitive to both p rosities and can be considered to build the rock-physics templates. Figure 11. Sensitivity indices varying with the total porosity, based on data from Well P. Figure 12 shows the RPTs with respect to total and microcrack porosities at ultrason and log frequencies, where the black and red curves isolines of total and microcrack p rosities, respectively.  Figure 12 shows the RPTs with respect to total and microcrack porosities at ultrasonic and log frequencies, where the black and red curves isolines of total and microcrack porosities, respectively.

Calibration
Ultrasonic and log data are used to calibrate the RPTs. Figure 13 displays an RPT at 1 MHz and the inclusion radius of 50 μm, where the color of the scatters indicates microcrack porosity, and the maximum value is 0.54%. The Poisson's ratio gradually increases with decreasing microcrack porosity; the plot shows that the template is in agreement with the ultrasonic data.

Calibration
Ultrasonic and log data are used to calibrate the RPTs. Figure 13 displays an RPT at 1 MHz and the inclusion radius of 50 µm, where the color of the scatters indicates microcrack porosity, and the maximum value is 0.54%. The Poisson's ratio gradually increases with decreasing microcrack porosity; the plot shows that the template is in agreement with the ultrasonic data. Figure 14 shows an RPT at 10 kHz and the inclusion radius of 200 µm, where the color bar indicates total porosity (the data are from the log data of well P). Most of the data have a porosity greater than 4%. With the increasing porosity, Poisson's ratio increases and P-wave impedance decreases, and the agreement between data and theory is good.

Calibration
Ultrasonic and log data are used to calibrate the RPTs. Figure 13 displays an RPT at 1 MHz and the inclusion radius of 50 μm, where the color of the scatters indicates microcrack porosity, and the maximum value is 0.54%. The Poisson's ratio gradually increases with decreasing microcrack porosity; the plot shows that the template is in agreement with the ultrasonic data.  data have a porosity greater than 4%. With the increasing porosity, Poisson's ratio increases and P-wave impedance decreases, and the agreement between data and theory is good. Figure 14. Log data and corresponding RPT. The color bar indicates the total porosity.

Microcrack Estimation
Based on approximated reflection coefficients [53], a pre-stack AVA inversion [54][55][56] is adopted to obtain the Poisson ratio and P-wave impedance and build the templates. Then, the microcrack porosity is estimated, with data points outside the template consid-

Microcrack Estimation
Based on approximated reflection coefficients [53], a pre-stack AVA inversion [54][55][56] is adopted to obtain the Poisson ratio and P-wave impedance and build the templates. Then, the microcrack porosity is estimated, with data points outside the template considered non-reservoir. Figure 15 shows the P-wave impedance and Poisson's ratio obtained from a seismic survey line. Well M shows high P-wave impedances at the top and low P-wave impedances at the bottom of the target layer, and high Poisson's ratio in the middle. On the other hand, Well F shows the low impedance and Poisson's ratio. Figure 16 shows the total and microcrack porosities predicted for the survey. Both porosities are high at Well M and low at Well F. The gas production reports of wells M and F are 100.5 × 10 4 m 3 per day and 4.4 × 10 4 m 3 per day, respectively, indicating that the predictions are consistent with the reports.   Figures 17 and 18 show results for another seismic survey line. Wells C and P show high P-wave impedances in the middle section, while Poisson's ratio is low throughout the target layer around Well C. This well exhibits low total and microcrack porosities (see Figure 18). Although Well P shows a higher total porosity, the microcrack porosity is low. Well J shows the best potential for gas production with high total and microcrack porosi-   Figures 17 and 18 show results for another seismic survey line. Wells C and P show high P-wave impedances in the middle section, while Poisson's ratio is low throughout the target layer around Well C. This well exhibits low total and microcrack porosities (see Figure 18). Although Well P shows a higher total porosity, the microcrack porosity is low. Well J shows the best potential for gas production with high total and microcrack porosi-  Figures 17 and 18 show results for another seismic survey line. Wells C and P show high P-wave impedances in the middle section, while Poisson's ratio is low throughout the target layer around Well C. This well exhibits low total and microcrack porosities (see Figure 18). Although Well P shows a higher total porosity, the microcrack porosity is low. Well J shows the best potential for gas production with high total and microcrack porosities. The reports show that Well C produced almost no gas with non-industrial gas production flow, while Wells P and J produce 2.14 × 10 4 and 6.4 × 10 4 m 3 gas per day, respectively. The predictions of the porosities are consistent with the production status.  3D seismic slices of the work area of ~471 km 2 are produced from the RPTs. According to the daily gas production rate, ten wells are classified into three categories, namely, extremely low (less than   3D seismic slices of the work area of ~471 km 2 are produced from the RPTs. According to the daily gas production rate, ten wells are classified into three categories, namely, extremely low (less than  3D seismic slices of the work area of~471 km 2 are produced from the RPTs. According to the daily gas production rate, ten wells are classified into three categories, namely, extremely low (less than 1 × 10 4 m 3 ), low (1 × 10 4 − 7 × 10 4 m 3 ) and high (higher than 7 × 10 4 m 3 ) gas production wells. Figures 19 and 20 show the corresponding maps of total and microcrack porosities, respectively. Energies 2021, 14, x FOR PEER REVIEW 15 of 20 Figure 19. 3D slice of reservoir total porosity. Black, white, and red circles indicate extremely low, low, and high gas production wells, respectively. The predictions show that Well A is located in an area with extremely low total and microcrack porosities. The total porosity of Wells C and D are higher than that of Well A, but their microcrack porosities are low. Wells F, G, J, and P are located in an area with low Figure 19. 3D slice of reservoir total porosity. Black, white, and red circles indicate extremely low, low, and high gas production wells, respectively.  The predictions show that Well A is located in an area with extremely low total and microcrack porosities. The total porosity of Wells C and D are higher than that of Well A, but their microcrack porosities are low. Wells F, G, J, and P are located in an area with low The predictions show that Well A is located in an area with extremely low total and microcrack porosities. The total porosity of Wells C and D are higher than that of Well A, but their microcrack porosities are low. Wells F, G, J, and P are located in an area with low total and microcrack porosities, whereas Wells E, N, and M show high total and microcrack porosities. According to the production test report, Well A is a water well with no gas production, and Wells C and D produce almost no gas. The three wells are classified as the lowest production. Wells P, F, G, and J produce 2.14 × 10 4 , 4.4 × 10 4 , 5.76 × 10 4 , and 6.4 × 10 4 m 3 gas per day, respectively, classified as low production. Wells E, N and M produce 17.66 × 10 4 , 25 × 10 4 , and 100.5 × 10 4 m 3 gas per day, respectively, and are classified as high production. The predictions basically agree with the reports, so that the RPTs can discriminate between low and high gas production reservoirs.

Conclusions
This study uses rock-physics templates to estimate the total and microcrack porosities of tight sandstone gas reservoirs. The models are based on the poroelasticity theory to obtain the elastic attributes at different frequency bands, namely, how Poisson's ratio, velocity, and other relevant attributes are affected by those porosities. Ultrasonic experiments are performed under full water and gas saturations, showing that the wave velocities increase with effective pressure. When the rock is fully saturated with gas, the ratio between the Pand S-wave velocities increases with effective pressure, while for water, an opposite trend is observed.
A sensitivity analysis shows that P-wave impedance and Poisson's ratio are most sensitive to both porosities. The attributes are considered to build the templates which are calibrated with experimental and log data. Finally, the templates are applied to seismic data to estimate the porosities, indicating that the predictions are consistent with the gas production reports.
At present, the tight sandstone gas reservoirs of Sichuan Basin have contributed to the major proportion of the total gas production output of the petroleum industries in China. However, the tight sandstone reservoirs generally exhibit the characteristics of low porosity, low permeability, high heterogeneity, and deep burial. There are difficulties in the precise prediction and detailed description of high-quality reservoirs based on the traditional geophysical exploration methods. In this work, we propose a multiscale rock physics modeling method by incorporating the effects of microcrack porosity on rock physics attributes. A workflow is presented to establish the relationships between rock physical properties and seismic wave responses, based on which the prediction of microcrack porosity is achieved in the 3D seismic survey work area. The application of the method is effective as a case study.
The proposed method can be applied in rock physics modeling based on the elastic/seismic responses of reservoir rock. It may be affected by the problem of multiplicity of solutions in the parameter inversion and reservoir prediction. Further studies may incorporate the wave attenuation in rock physics modeling, so that more attributes will be available for validating and calibrating the templates, and the influence by the multiplicity of solutions can be alleviated. Furthermore, this paper is focused on the studies on tight gas sandstone reservoirs, and for the different types of complex reservoirs, the approach should be extended and adjusted to achieve effective applications. Specifically, for those tight/shale oil reservoirs, the viscoelasticity characteristics of fluid have to be analyzed and incorporated in rock physics modeling, so that some reasonable results can be obtained.
For future studies, the different sensitivity analysis methods should be included to analyze the interactions between different attributes with respect to rock basic properties. Furthermore, the multi-phases of different types of microcracks, mesopores and fractures need be considered to establish a more realistic rock model, with which the in situ reservoirs can be better described and simulated, in combination with the latest techniques of digital rock physics. (

Conflicts of Interest:
The authors declare no conflict of interest.

AVO Amplitude variation with offset BR
Biot-Rayleigh DEM Differential equivalent medium RPT Rock physics template FSI Fluid sensitivity indicator