Abstract
Tight-sandstone reservoirs have a complex pore structure with microcracks and intergranular pores, which have a significant impact on the seismic properties. We have performed ultrasonic measurements at different confining pressures for 15 tight-gas sandstone samples of the Xujiahe formation in Western Sichuan Basin, and have available well-log and seismic data of this area. The aim of this work is to estimate the porosity and crack properties for variable pressure conditions. The EIAS (equivalent inclusion-average stress) model is adopted to compute the high- and low-frequency bulk and shear moduli as a function of crack aspect ratio and (soft) and (stiff) porosities. Then, we use the EIAS-Zener anelastic model to obtain the wave properties as a function of frequency, and compare results with those of the constant Q (Kjartansson) one for verification of the robustness of the approach. The corresponding P-wave impedance, density and phase velocity ratio (VP/VS) are computed in order to built 3D rock‐physics templates (RPTs) at the ultrasonic, well-log and seismic frequency bands. The methodology is applied to a survey line crossing two wells, which together with the laboratory experiments, provide calibration suitable data. The estimated stiff porosity and crack porosity and density are consistent with the available data and actual production records, indicating that 3D RPTs provide a useful interpretation tool in seismic exploration and prospect evaluation.
Introduction
Unconventional gas resources offer significant gas production growth potential. Tight-gas sandstones represent an important part of the unconventional production and abundant reserves are yet to be developed (). There are tight-gas sandstone reservoirs almost in all oil and gas bearing areas with significant gas reserves (Zhu et al., 2008). These reservoirs have sandstones with low porosity, permeability and gas saturation and abundant microcracks (, ; ; ). An effective reservoir characterization is the key to their successful exploration (; ). The reservoir contains a large number of natural fractures, which provide main seepage channels and effective storage space, and favor single-well productivity (; Xiao et al., 2019). Therefore, an accurate identification of the pore properties is essential.
In this case, rock physics is required to relate the wave properties to the composition and microstructure of these rocks, which have been affected by complex diagenetic and sedimentary processes during their formation (Vernik and Kachanov, 2010; Zhang et al., 2019; ). Walsh (1965) developed some of the first techniques to study the effects of cracks and confining pressure on the dry-rock moduli. Thomsen (1995) considered aligned circular cracks and analyzed the effects of the associated anisotropy on the elastic properties, and quantified how microcracks affect the wave velocities in tight-gas sandstone reservoirs. Tang (2011) proposed an elastic-wave model to describe the effects of pores and cracks by extending the Biot and BISQ theories. Double-porosity models that consider two phases with different compressibilities are relatively new and useful approaches to model wave propagation in heterogeneous media (; ; ; ).
To model the relation between pore geometry and elasticity attributes, we used the EIAS (equivalent inclusion-average stress) model (), which is consistent with Gassmann equation at low frequencies and with the Hashin-Shtrikman bounds when applied to two-phase systems. Moreover, we used the EIAS-Zener model that generalizes the EIAS model to the full frequency range, by incorporating the Zener (standard linear solid) model (). The crack porosity, average crack aspect ratio and quality factors and phase velocities were obtained as a function of frequency by fitting experimental data at the unrelaxed and relaxed states.
RPTs establish a link between the reservoir properties (e.g., porosity, fluid saturation, clay content, etc) and the elastic properties, such as velocity, density, impedance, and wet-rock stiffness moduli (; ; ; Tan et al., 2020). used the P-wave quality factor and impedance to build RPTs, and shown that the estimated porosity and saturation are consistent with actual gas production results and well-log data. used RPTs based on VP/VS and acoustic impedance applied to carbonate samples.
We used VP/VS, impedance and density to build 3D RPTs, based on the EIAS model and the Zener frequency kernels, to characterize the relation between elastic properties and reservoir attributes, i.e., equant porosity and crack porosity and aspect ratio (). The templates were calibrated with ultrasonic, log and seismic data. Thereafter, a quantitative interpretation of those attributes can be performed by superposing the seismic-inversion data on the templates.
Laboratory Experiments
The Study Area
The reservoir in the Xujiahe formation of the western Sichuan Basin (China) belongs to the braided river delta sedimentary system. The main mineral components are quartz, feldspar, rock debris, mica, and heavy minerals. The size of the mineral particles is mainly medium to fine, most sandstones are sorted well, and the crushed particles are medium to poorly rounded. These sandstone reservoirs, of low porosity and permeability, have a high crack density and a heterogeneous distribution of pore fluids, due to long-term diagenesis.
Rock Specimens and Data
Fifteen samples were collected at approximately 2 km depth. They have low porosity and permeability, and are mainly composed of quartz, feldspar and rare clay. Moreover, there are calcite and quartz cementation in intergranular pores, and cracks. The ultrasonic experiments measurements (samples TSA1-TSA15) follow the experimental set-up of , which is composed of a high pressure vessel and several units to control the temperature and the confining and pore pressures, and perform acoustic wave testing. The acquisition rate for ultrasonic waveforms (1 MHz) is 50 M/s and the time resolution is 0.02 μs.
The tests were performed at zero pore pressure and 22°C, whereas the confining pressure increases from 1 to 35 MPa (1, 5, 10, 15, 20, 25, 30 and 35 MPa). The differential pressure (Pd) is defined here as the difference between the confining and pore pressures. Ultrasonic compressional and shear-wave velocities (VP and VS) at full-gas (nitrogen) were obtained by picking the first arrivals of waveforms with 1 MHz dominant frequency. These samples have low porosity and permeability, with a maximum of 13.91% and 1.37 mD, respectively. The properties are given in Table 1. Figure 1 shows a crossplot of VP and VS for varying differential pressure. As can be seen, the velocities increase with pressure, as expected. The correlation of a linear fitting between VP and VS is excellent, the R-square is 0.956, and the slope of VP vs. VS is 1.6.
TABLE 1
| Sample | Total porosity | Permeability (mD) | Mineral density (g cm−3) |
|---|---|---|---|
| TSA1 | 0.0243 | <0.001 | 2.695 |
| TSA2 | 0.0303 | <0.001 | 2.687 |
| TSA3 | 0.0346 | 0.005 | 2.694 |
| TSA4 | 0.0388 | <0.001 | 2.670 |
| TSA5 | 0.0419 | <0.001 | 2.676 |
| TSA6 | 0.0439 | 0.004 | 2.691 |
| TSA7 | 0.0626 | 0.046 | 2.672 |
| TSA8 | 0.0633 | 0.047 | 2.665 |
| TSA9 | 0.0757 | 0.0225 | 2.679 |
| TSA10 | 0.0855 | 0.082 | 2.670 |
| TSA11 | 0.0865 | 0.028 | 2.668 |
| TSA12 | 0.0897 | 0.140 | 2.662 |
| TSA13 | 0.1326 | 1.21 | 2.659 |
| TSA14 | 0.1335 | 1.320 | 2.653 |
| TSA15 | 0.1391 | 1.370 | 2.660 |
Tight-gas sandstones properties.
FIGURE 1
Figure 2 shows ultrasonic velocities as a function of differential pressure at full gas saturation. Velocities are non-linear at low differential pressure, becoming almost linear above 25 MPa, because cracks with lower aspect ratio close first and when the differential pressure exceeds approximately 25 MPa, the stiffer pores are also affected. The velocities increase faster for samples TSA4, TSA5 and TSA8, that can be due to a higher crack density, resulting in a higher sensitivity to pressure. Similar results have been observed previously by Yin et al. (2017), , , , Sun and Goldberg (1997) with different types of sandstones.
FIGURE 2
Determination of the Experimental Crack Porosity
Cracks with lower aspect ratio close first when the differential pressure increases, leading to higher velocities and decreasing crack density. We use sample TSA13 as an example to estimate crack porosity and the pressure dependency of the total, stiff (equant) and crack porosity (Figure 3) (The former is the sum of the last two porosities). The total porosity as a function of pressure is obtained by measurement, whereas the stiff porosity by a linear extrapolation of the high-pressure trend determined between 30 and 35 MPa. The similar method was used by Yin et al. (2017) and . The crack porosity is estimated by subtracting the stiff one from the total porosity. Figure 4 shows the crack porosity as a function of the differential pressure. When the differential pressure is less than 25 MPa, the crack porosity decreases rapidly and non-linearly. As the pressure increases, the crack porosity exhibits linear characteristics until it approaches zero at high differential pressures.
FIGURE 3
FIGURE 4
Rock-Physics Model
The EIAS Model
The microscopic pore structure is important because it affects the rock skeleton properties and the fluid distribution. The EIAS model can effectively address the relation between pore geometry and the elastic attributes, and it is consistent with the Hashin-Shtrikman bounds for two-phase systems (; ). Assuming a background medium containing spherical pores and cracks, the resulting expressions for the effective high-frequency elastic moduli of the saturated rock arewhere ϕ denotes the rock total porosity, K0 and G0 are the bulk and shear moduli of the mineral mixture, respectively, Kf is the fluid bulk modulus and and are intermediate variables, which can be expressed as (; ; ; Sun et al., 2020), where a is the crack aspect ratio, c is the soft-pore (crack) fraction and 1 - c is the stiff-pore fraction; the crack density is . The coefficients P1 and Q1 correspond to spherical pores and P2 and Q2 are approximations for penny-shaped cracks.
The effective low-frequency moduli (when complete fluid-pressure communication occurs) arewhereand
The detailed derivation of these equations can be found in . The model considers the interactions between cracks and has no restrictions on the crack density. In order to generalize the model to all frequencies, we combine the EIAS model with the Zener (standard linear solid) and Kjartansson kernels, to describe the physics with two very dissimilar approaches and check the robustness of the methodology.
EIAS-Zener Model
The Zener mechanical model can be used to described the frequency dependence of dispersion and attenuation (e.g., ; ). This model satisfies the Kramers-Kronig relations (). The minimum quality factors of the bulk and shear relaxation peaks are given byrespectively, whereas the bulk and shear complex moduli arewhere , f0 is the frequency of the relaxation peak, such that for , and for , .
The phase velocity and quality factor of the body waves areandrespectively, where denotes or , being the frequency-dependent complex P- and S-wave velocities
(), respectively, where ρ is the mass density.
EIAS-Kjartansson Model
presented a constant-Q model specified by the phase velocity at a reference frequency f1 and the value of Q. To combine this model with the EIAS model, we consider f1 = 50 Hz, which is the average frequency of the event corresponding to the target layer, i.e., rad/s. The phase velocity at low frequencies is
For f2 = 1 MHz (the ultrasonic experimental frequency), rad/s, and the phase velocity is
The phase velocity as a function of frequency is (, Eq. 2.214),as an approximation for Q>>1. We obtain after some calculations,Defining,the complex modulus is (, Eq. 2.212)
Therefore, the complex velocity and quality factor arerespectively.
Example
The previous models are used to relate the crack properties (aspect ratio and porosity) and the frequency-dependent elastic moduli, phase velocities and quality factors. The purpose of using two dissimilar models is to verify the robustness of the approach, i.e., if the two models give similar results, we can rely on the methodology. Let us consider sample TSA13 (porosity is 13.26%). We assume that the mineral bulk and shear modulus are 39 and 36 GPa, respectively (), the fluid properties are obtained from . We consider a Zener kernel with f0 = 10 KHz and the parameters of the constant-Q kernel assumed in the previous section. Figure 5 shows the phase velocities and dissipation factors as a function of frequency, where it is clear both kernels honor the velocity dispersion between the seismic and laboratory frequency bands. Therefore, they can be used to build multi-scale 3D RPTs in the ultrasonic, logging and seismic frequency bands.
FIGURE 5
Thereafter, we assume a crack porosity of 0.2% and crack aspect ratios of 0.001, 0.0014, 0.0018, 0.0022, 0.0026. Velocities and dissipation factors as a function of frequency are shown in Figures 6A,B. When the crack porosity is a constant and the crack aspect ratio increases, the rock becomes stiffer and the P-wave velocity increases, but attenuation decreases. If we take a crack aspect ratio of 0.002 and crack porosities of 0.2, 0.3, 0.4, 0.5, 0.6%, the results are shown in Figures 6C,D, where the trend opposite to the previous one.
FIGURE 6
Multi-Scale 3D RPTs
Building the RPTs
By setting the crack aspect ratio, crack porosity and total porosity as variables, the corresponding P-wave impedance, density and VP/VS value are calculated by using the EIAS-Zener model. Then, we build the 3D RPTs at ultrasonic (1 MHz) well-log (10 kHz) and seismic (50 Hz) frequencies, and the seismic template is used to estimate the microstructure of the sandstone reservoirs. Similar templates have been used by to estimate the microfracture porosity of deep carbonate reservoirs. Figures 7A,B show the 3D RPTs at ultrasonic and seismic frequencies, respectively (the well-log template is similar to the seismic one).
FIGURE 7
As the stiff porosity increases, the P-wave impedance and mass density decrease, but the variation in VP/VS is small. With the increase of crack porosity, the P-wave impedance decreases and VP/VS increases, but the density hardly changes. As the crack aspect ratio increases, the P-wave impedance increases, VP/VS decreases. This behavior is similar at all frequencies.
Calibration at Ultrasonic Frequencies
The templates are calibrated with the crack porosity obtained in Section Determination of the Experimental Crack Porosity. Figure 8 superposes the data to the template. The porosity (2.43–13.91%) of samples is basically consistent with the corresponding values in the templates, and the agreement is also good for the crack porosity. What’s more, the average error between the predicted porosities (based on ultrasonic wave attributes) and the measured one is 0.00455. In addition, as this increases, the P-wave impedance decreases.
FIGURE 8
Calibration at Well-Log and Seismic Frequencies
Figure 9 shows the calibration of the 3D RPTs at log and seismic frequencies, where the data is extracted from well B. The template at seismic frequencies is calibrated with inversion results (P-wave impedance, density and VP/VS) of pre-stack seismic data that has been verified, whole the porosity is taken from the logs. Similar methods have been used by . The velocities and density at seismic frequency band are extracted from the seismic traces around the borehole by using a three-term inversion method (). We use the partial-stack method to obtain the seismic data based on pre-stack angle gathers, which improves the signal-noize ratio and suppresses random noises. In addition, the velocities and porosity at sonic-log frequency band are extracted from the well-log data by using the velocity-time semblance method, which avoids picking arrival times and can accurately determine the phase velocities for different modes (Wang et al., 2020), and the density is obtained from density logging. Figure 9 shows that the porosity is basically consistent with the corresponding values in the templates. Figure 10 shows a crossplot of VP and VS of well B, where the velocities increase with decreasing porosity, as expected. The correlation of a linear fitting between VP and VS is excellent, and the slope of VP vs. VS is 1.12.
FIGURE 9
FIGURE 10
The P-wave impedance, density and VP/VS decreases with increasing porosity, and the latter decreases with increasing crack aspect ratio. The porosity is low (less than 7%), which is consistent with the geological characteristics of the reservoirs.
Estimation of the Microstructure from Seismic Data
The predictions of reservoir stiff porosity and crack porosity and density are performed on a 2D seismic line crossing wells A and B. P-wave impedance, density and VP/VS are overlapped onto the 3D RPTs, and the microstructural properties are estimated.
Figure 11A shows the porosity from well log A and B are superposed, showing a good agreement. By taking well B as an example, the average error between the predicted porosities (based on seismic wave attributes) and the measured porosities is 0.00923, which is higher than the one at ultrasonic frequency, because the resolution of the seismic data is low. Figures 11B,C show the crack porosity and density, respectively. It shows that well B has a higher crack porosity and density. The horizontal sections are created from target layer (depth 4,787 m) of the study area in the Xujiahe formation of the western Sichuan Basin (China). Figures 12, 13 show horizontal section of those attributes, which are simulated by using EIAS-Zener and EIAS-Kjartansson models, respectively. The results show that the two models give similar results, thus indicating that the methodology is working. The total porosity and crack porosity and density effectively serve to identify high-quality reservoirs. Well B is located in an area with high porosity and crack content, where the target layer has a good pore connectivity, storage capability, and natural gas accumulation conditions. The gas production report of well B indicates 5148929 m3 per day, while well A yields 57559 m3 per day. The results agree with the actual production of the two wells. Therefore, the methodology is reliable, and is expected to provide an important basis for further rock-physics studies and field applications of seismic exploration.
FIGURE 11
FIGURE 12
FIGURE 13
Conclusion
Rock-physics templates constitute an easy-to-use “tool-box” for pore microstructure properties interpretation from seismic data. We have implemented the EIAS-Zener model to describe the crack properties and elastic moduli of tight-gas sandstone reservoirs as a function of frequency from the seismic to the ultrasonic bands. Since P-wave impedance, mass density and VP/VS are sensitive to the pore and crack properties, we built the templates to estimate these attributes from seismic data. First, we calibrate the templates with ultrasonic, well-log and seismic data. Then, we estimate the total porosity and crack porosity and density from inversion of seismic data. The results indicate that these estimated attributes are consistent with the actual gas production reports. Thus, these templates can effectively be used to estimate the microstructure characteristics of tight-gas sandstone reservoirs based on seismic data.
Statements
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author contributions
WC: Methodology, writing original manuscript JB: Funding support, data and discussion, writing original manuscript, supervision JC: Model building, revising manuscript MP: Seismic data analysis CW: Ultrasonic experimental data analysis.
Funding
The authors are grateful to the support of the Jiangsu Innovation and Entrepreneurship Plan, Jiangsu Province Science Fund for Distinguished Young Scholars (BK20200021), National Natural Science Foundation of China (41974123), China Postdoctoral Science Foundation (2021M690863), Fundamental Research Funds for the Central Universities (B210202112) and the research funds from SINOPEC Key Laboratory of Geophysics.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
1
AkiK.RichardsP. G. (1980). Quantitative Seismology, Theory and Methods. San Francisco: W.H. Freeman and Company.
2
AnjosS. M. C.De RosL. F.SilvaC. M. A. (2003). Chlorite Authigenesis and Porosity Preservation in the Upper Cretaceous marine Sandstones of the Santos Basin, Offshore Eastern Brazil. Int. Assoc. Sedimentologists34, 289–316. 10.1002/9781444304336.ch13
3
BaJ.XuW.FuL.-Y.CarcioneJ. M.ZhangL. (2017). Rock Anelasticity Due to Patchy Saturation and Fabric Heterogeneity: A Double Double-Porosity Model of Wave Propagation. J. Geophys. Res. Solid Earth122 (3), 1949–1976. 10.1002/2016JB013882
4
BaJ.ZhaoJ.CarcioneJ. M.HuangX. (2016). Compressional Wave Dispersion Due to Rock Matrix Stiffening by clay Squirt Flow. Geophys. Res. Lett.43, 6186–6195. 10.1002/2016GL069312
5
BatzleM.WangZ. (1992). Seismic Properties of Pore Fluids. Geophysics57 (11), 1396–1408. 10.1190/1.1443207
6
BerrymanJ. G. (1995). Mixture Theories for Rock Propertie. A Handbook of Physical Constants. Washington, D C: American Geophysical Union, 205–228. 10.1029/RF003p0205
7
CarcioneJ. M.AvsethP. (2015). Rock-physics Templates for clay-rich Source Rocks. Geophysics80 (5), D481–D500. 10.1190/geo2014-0510.1
8
CarcioneJ. M.CavalliniF.BaJ.ChengW.QadrouhA. (2018). On the Kramers-Kronig Relations. Rheologica Acta58 (1-2), 21–28. Available at: https://link.springer.com/article/10.1007%2Fs00397–018-1119–3. 10.1007/s00397-018-1119-3
9
CarcioneJ. M.FarinaB.PolettoF.QadrouhA. N.ChengW. (2020). Seismic Attenuation in Partially Molten Rocks. Phys. Earth Planet. Interiors309, 106568. 10.1016/j.pepi.2020.106568
10
CarcioneJ. M.GeiD.PicottiS.MicheliniA. (2012). Cross-hole Electromagnetic and Seismic Modeling for CO2 Detection and Monitoring in a saline Aquifer. J. Pet. Sci. Eng.100, 162–172. 10.1016/j.petrol.2012.03.018
11
CarcioneJ. M. (2014). Wave Fields in Real Media. Theory and Numerical Simulation of Wave Propagation in Anisotropic, Anelastic, Porous and Electromagnetic media. Third edition. Elsevierextended and revised.
12
ChengW.BaJ.FuL.-Y.LebedevM. (2019). Wave-velocity Dispersion and Rock Microstructure. J. Pet. Sci. Eng.183, 106466. 10.1016/j.petrol.2019.106466
13
ChengW.CarcioneJ. M.QadrouhA. N.AlajmiM.BaJ. (2020). Rock Anelasticity, Pore Geometry and the Biot-Gardner Effect. Rock Mech. Rock Eng.53, 3969–3981. 10.1007/s00603-020-02155-7
14
DavidE. C.FortinJ.SchubnelA.GuéguenY.ZimmermanR. W. (2013). Laboratory Measurements of Low- and High-Frequency Elastic Moduli in Fontainebleau sandstone. Geophysics78 (5), D369–D379. 10.1190/geo2013-0070.1
15
DengJ.ZhouH.WangH.ZhaoJ.WangS. (2015). The Influence of Pore Structure in Reservoir sandstone on Dispersion Properties of Elastic Waves. Chin. J. Geophysics-Chinese58 (9), 3389–3400. 10.6038/cjg20150931
16
DvorkinJ.WallsJ.DavalosG. (2020). Velocity-Porosity-Mineralogy Model for Unconventional Shale and its Applications to Digital Rock Physics. Front. Earth Sci.8. 10.3389/feart.2020.613716
17
EndresA. L.KnightR. J. (1997). Incorporating Pore Geometry and Fluid Pressure Communication into Modeling the Elastic Behavior of Porous Rocks. Geophysics62 (1), 106–117. 10.1190/1.1444110
18
FuB.-Y.GuoJ.FuL.-Y.GlubokovskikhS.GalvinR. J.GurevichB. (2018). Seismic Dispersion and Attenuation in Saturated Porous Rock with Aligned Slit Cracks. J. Geophys. Res. Solid Earth123, 6890–6910. 10.1029/2018JB015918
19
GegenhuberN.PuposJ. (2015). Rock Physics Template from Laboratory Data for Carbonates. J. Appl. Geophys.114, 12–18. 10.1016/j.jappgeo.2015.01.005
20
GuanD. S.NiuJ. Y. (1995). Unconventional Oil and Gas Geology in China. Beijing: Petroleum Industry Press, 60–85.
21
GuoM.-Q.FuL.-Y.BaJ. (2009). Comparison of Stress-Associated Coda Attenuation and Intrinsic Attenuation from Ultrasonic Measurements. Geophys. J. Int.178, 447–456. 10.1111/j.1365-246X.2009.04159.x
22
GuoM. Q.BaJ.MaR. P.ChenT. S.ZhangL.PangM. Q.et al (2018). P-wave Velocity Dispersion and Attenuation in Fluid-Saturated Tight Sandstones: Characteristics Analysis Based on a Double Double-Porosity Structure Model Description. Chin. J. Geophys.61 (3), 1053–1068. 10.6038/cji2018L0678
23
HaoZ.BaJ.ZhangL.ZengQ.JiangLiuR. J.QianW.et al (2016). Rock Physics Inversion Workflow on Reservoir Parameters: A Case Study of Seismic Hydrocarbon Detection in Large-Area Tight Dolomite Reservoirs. J. Seismic Exploration25, 561–588.
24
HarmelenA. V.WeijermarsR. (2018). Complex Analytical Solutions for Flow in Hydraulically Fractured Hydrocarbon Reservoirs with and without Natural Fractures. Appl. Math. Model.56, 137–157. 10.1016/j.apm.2017.11.027
25
KhlaifatA.QutobH.BarakatN. (2011). Tight Gas Sands Development Is Critical to Future World Energy Resources. Society of Petroleum Engineers, 1–12. 10.2118/142049-MS
26
KjartanssonE. (1979). ConstantQ-wave Propagation and Attenuation. J. Geophys. Res.84, 4737–4748. 10.1029/jb084ib09p04737
27
LiH. B.ZhangJ. J.CaiS. J.PanH. J. (2019). 3D Rock Physics Template for Reservoir with Complex Pore Structure. Chin. J. Geophys. (in Chinese)62, 2711–2723. 10.6038/cjg2019K0672
28
LiuL.TangD.WoY.LiuL.SunW. (2019). Favorable Area Prediction of Tight sandstone: A Case Study of the He8 Formation in the Kangning Area, Eastern Ordos Basin, China. J. Pet. Sci. Eng.175, 430–443. 10.1016/j.petrol.2018.12.069
29
MavkoG.MukerjiT.DvorkinJ. (2009). The Rock Physics Handbook. Cambridge Univ. Press. 10.1017/cbo9780511626753
30
PangM.BaJ.CarcioneJ. M.PicottiS.ZhouJ.JiangR. (2019). Estimation of Porosity and Fluid Saturation in Carbonates from Rock-Physics Templates Based on Seismic Q. Geophysics84 (6), M25–M36. 10.1190/geo2019-0031.1
31
PangM.BaJ.CarcioneJ.VesnaverA.MaR.ChenT. (2020b). Analysis of Attenuation Rock-Physics Template of Tight Sandstones: Reservoir Microcrack Prediction. Chin. J. Geophys.63, 4205–4219.
32
PangM.BaJ.FuL. Y.CarcioneJ. M.MarkusU.ZhangL. (2020a). Estimation of Microfracture Porosity in Deep Carbonate Reservoirs Based on 3d Rock-Physics Templates. Interpretation8 (4), 1–41. 10.1190/int-2019-0258.1
33
PervukhinaM.GurevichB.DewhurstD. N.SigginsA. F. (2010). Applicability of Velocity-Stress Relationships Based on the Dual Porosity Concept to Isotropic Porous Rocks. Geophys. J. Int.181, 1473–1479. 10.1111/j.1365-246X.2010.04535.x
34
PrideS. R.BerrymanJ. G.HarrisJ. M. (2004). Seismic Attenuation Due to Wave-Induced Flow. J. Geophys. Res.109 (B1). 10.1029/2003jb002639
35
SmithT. M.SayersC. M.SondergeldC. H. (2009). Rock Properties in Low-Porosity/low-Permeability Sandstones. The Leading Edge28 (1), 48–59. 10.1190/1.3064146
36
StrokerT. M.HarrisN. B.Crawford ElliottW.Marion WamplerJ. (2013). Diagenesis of a Tight Gas Sand Reservoir: Upper Cretaceous Mesaverde Group, Piceance Basin, Colorado. Mar. Pet. Geology.40, 48–68. 10.1016/j.marpetgeo.2012.08.003
37
SunY.CarcioneJ. M.GurevichB. (2020). Squirt-flow Seismic Dispersion Models: A Comparison. Geophys. J. Internat.222, 2068–2082. 10.1093/gji/ggaa274
38
SunY. F.GoldbergD. (1997). Estimation of Aspect-Ratio Changes with Pressure from Seismic Velocities. Geol. Soc. Lond. Spec. Publications122, 131–139. 10.1144/GSL.SP.1997.122.01.10
39
TanW.BaJ.MüllerT.FangG.ZhaoH. (2020). Rock Physics Model of Tight Oil Siltstone for Seismic Prediction of Brittleness. Geophys. Prospecting68 (5), 1554–1574. 10.1111/1365-2478.12938
40
TangX. M. (2011). A Unified Theory for Elastic Wave Propagation through Porous media Containing Cracks—An Extension of Biot’s Poroelastic Wave Theory. Sci. China Earth Sci. (In Chinese)41 (6), 784–795.
41
ThomsenL. (1995). Elastic Anisotropy Due to Aligned Cracks in Porous Rock1. Geophys. Prospect43, 805–829. 10.1111/j.1365-2478.1995.tb00282.x
42
VernikL.KachanovM. (2010). Modeling Elastic Properties of Siliciclastic Rocks. Geophysics75 (6), E171–E182. 10.1190/1.3494031
43
WalshJ. B. (1965). The Effect of Cracks on the Compressibility of Rock. J. Geophys. Res.70, 381–389. 10.1029/JZ070i002p00381
44
WangH.ToksözM. N.FehlerM. C. (2020). Borehole Acoustic Logging – Theory and Methods. Petroleum Engineering. 10.1007/978-3-030-51423-5
45
XiaoZ.DingW.LiuJ.TianM.YinS.ZhouX.et al (2019). A Fracture Identification Method for Low-Permeability sandstone Based on R/S Analysis and the Finite Difference Method: A Case Study from the Chang 6 Reservoir in Huaqing Oilfield, Ordos Basin. J. Pet. Sci. Eng.174, 1169–1178. 10.1016/j.petrol.2018.12.017
46
YinH.ZhaoJ.TangG.ZhaoL.MaX.WangS. (2017). Pressure and Fluid Effect on Frequency-dependent Elastic Moduli in Fully Saturated Tight sandstone. J. Geophys. Res. Solid Earth122, 8925–8942. 10.1002/2017JB014244
47
ZhangL.BaJ.FuL.CarcioneJ. M.CaoC. (2019). Estimation of Pore Microstructure by Using the Static and Dynamic Moduli. Int. J. Rock Mech. Mining Sci.113, 24–30. 10.1016/j.ijrmms.2018.11.005
48
ZhuR.ZouC.ZhangN.WangX.ChengR.LiuL.et al (2008). Diagenetic Fluids Evolution and Genetic Mechanism of Tight sandstone Gas Reservoirs in Upper Triassic Xujiahe Formation in Sichuan Basin, China. Sci. China Ser. D-earth Sci.51 (9), 1340–1353. 10.1007/s11430-008-0102-8
Summary
Keywords
tight-gas sandstone, rock-physics templates, porosity, crack porosity, EIAS-Zener model, attenuation
Citation
Cheng W, Ba J, Carcione JM, Pang M and Wu C (2021) Estimation of the Pore Microstructure of Tight-Gas Sandstone Reservoirs with Seismic Data. Front. Earth Sci. 9:646372. doi: 10.3389/feart.2021.646372
Received
26 December 2020
Accepted
07 June 2021
Published
29 June 2021
Volume
9 - 2021
Edited by
Xin Li, China National Offshore Oil Corporation, China
Reviewed by
Luca De Siena, Johannes Gutenberg University Mainz, Germany
David Iacopini, University of Naples Federico II, Italy
Tingting Gu, Purdue University, United States
Updates
Copyright
© 2021 Cheng, Ba, Carcione, Pang and Wu.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jing Ba, jingba@188.com
This article was submitted to Solid Earth Geophysics, a section of the journal Frontiers in Earth Science
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.