Abstract
Carbonate reservoirs exhibit complex pore structure, which significantly affects the elastic properties and seismic response, as well as the prediction of physical parameters. As one of the main factors impacting fluid prediction, pore structure parameter directly involves in few inversion methods. In order to directly predict pore structure parameter in inversion, a novel quantitative reflection coefficient formula is proposed, that integrate Russell's poroelasticity theory with Sun's petrophysical model. This formula separates fluid bulk modulus from porosity and pore structure parameter, allowing for accurate determination of pore-fluid distribution through Bayesian framework. Both theoretical model analysis and multi-component digital core experiments of carbonates validate the importance of pore structure parameter on fluid identification. The practical application of carbonate reservoirs in Sichuan Basin demonstrates that the proposed fluid factor, eliminating the prediction illusion caused by heterogeneity in porosity and pore structure parameter within strata, provides more precise and reliable predictions compared to the Russell fluid factor. Furthermore, the similarity between the Russell fluid factor obtained directly from the Russell approximation and the Russell fluid factor calculated indirectly from the proposed method confirms the stability and accuracy of the new reflection coefficient formula.
1 Introduction
With the advancement of exploration and development technology, the targets of exploration and development are becoming increasingly intricate. There has been a decline in reserves and grades of both oil and gas, with a growing proportion of low porosity and low permeability. In terms of carbonate reservoirs, predicting reservoirs is much challenging due to their deep burial depth and weak seismic responses. Moreover, they exhibit strong heterogeneity and significant variations in physical properties and thickness over short distances, as well as complex pore types. All these characteristics contribute to heightened difficulty in identifying oil and gas within them. Therefore, one of the current challenges in geophysics lies in effectively identifying fluids within target reservoirs with complex porous media.
The effect of porosity and pore types on the effective elastic properties of reservoirs is very critical in hydrocarbon prediction. A study conducted by ) showed that the variations in the pore aspect ratio can lead to changes in wave velocity exceeding 2000 m/s or more. Currently, significant advancements have been achieved in the research of pore structure parameters. For example, ) and ) examined the influence of two- and three-dimensional pore shape on elastic properties, respectively. ) derived approximate analytical expressions for the elastic parameters of dry/saturated fractured rock based on DEM theory. ) proposed a new petrophysical modeling method using the Gassmann equation and Eshelby–Walsh ellipsoidal fracture theory to extract a parameter characterizing variations in the pore structure, thereby demonstrating its profound influence on the elastic properties of the rock. ) provided an analytical expression for the initial minimum aspect ratio of soft pores based on poroelasticity theory and extended the “squirt flow” model based on the characteristic aspect ratio by adding soft pores iteratively to analyze the influence of complex pore distribution on the squirt flow and its possible velocity dispersion characteristics. ) established a pore-scale numerical simulation method for elastic wave propagation in porous media based on the extraction of the pore structure parameter from digital core images in view of the complex pore structure and significant heterogeneity of carbonate rocks, verified its reliability with natural core data, and quantitatively analyzed the effect of carbonate pore structures on the propagation velocity and the scattering attenuation of elastic waves. ) proposed a unified expression by consolidating various petrophysical models for porous media, including Pride, Geertsma, and Keys-Xu, aiming to establish a more universal application across different rock types. Subsequently, ) further emphasized the significance of the pore aspect ratio as a fundamental benchmark due to the ambiguous physical definition and inconsistent quantitative representation of the pore structure parameter. These studies collectively indicate that the current research on the pore structure parameter primarily focuses on small-scale cores and logs but lacks applicability for fluid prediction at the seismic scale.
In the past few decades, there has been significant advancement in the fluid identification technology of seismic reservoirs, from a “qualitative” approach based on seismic amplitude anomaly in the 1980s to a “quantitative” method relying on the fluid factor in the present stage. With the development of a prestack seismic inversion, the fluid sensitive terms have evolved from assessing the relative variation in elastic parameters to presently recognizing physical parameters with distinct petrophysical significance. ) and ) both proposed the construction methods of the fluid factor for porous fluid-saturated rocks. ) and ) used the Biot–Gassmann theory to refine the Aki–Richards approximation under saturated fluid conditions and introduced the fluid factor, which can be directly involved in seismic fluid detection. ) derived a seismic reflection coefficient formula that incorporates porosity by combining the Russell approximation with the Nur model. ) then extended this formula to PP–PS joint inversions, which further improves the fluid prediction in heterogeneous reservoirs but does not consider the pore structure parameter. ) and ) established a direct relationship among fluid factors and P- and S-wave moduli based on the petrophysical model of porous elastic media to obtain a new Zoeppritz approximation formula based on the P- and S-wave moduli, which circumvents the issue related to accurately determining density fluid factor calculation and has been successfully applied in an exploration area in eastern China. ) proposed a novel fracture fluid factor that can simultaneously detect fracture development and fluid properties by combining P-wave anisotropic fracture prediction and Russell fluid factor into the Cartesian coordinate system in order to address the challenge of fluid identification in anisotropy, which achieved promising application in igneous areas. ) employed sequential Gaussian simulation and Metropolis sampling algorithm based on Bayesian’s theoretical framework to directly estimate the Russell fluid factor, which enhanced the accuracy of fluid factor identification. Although substantial developments have been made in recent years regarding the algorithms of fluid prediction, the issue of fluid prediction under non-homogeneous conditions still faces challenges such as strong multi-solution and less precise prediction. This can be attributed to the omission of pore structure parameters in the widely used Gassmann fluid equation for bidirectional media containing fluid, making it complicated and challenging to accurately quantify the impact of pore structure. ) and ) introduced the squirt flow model into an inversion, along with the relevant parameters of the pore structure, making an initial step toward integrating the pore structure parameter in seismic reservoir prediction. ) also proposed a new method combining the pore structure with the Bayesian non-linear simultaneous inversion of physical parameters, further improving the importance and involvement of the pore structure in the reservoir prediction. However, the aforementioned methods including pore structures can only be obtained in an indirect way. In other words, traditional physical parameters are first calculated, and then, specific petrophysical theories are employed to derive the pore structure and other fluid parameters. This approach does not directly incorporate the pore structure parameter as the factor in the inversion process of seismic reflection coefficient. Therefore, it is imperative to further study the inversion method of the pore structure from indirect to direct inversions in order to streamline the procedure, while improving the inversion accuracy.
In order to mitigate the impact of complex pore structure on fluid prediction in carbonate reservoirs, a novel inversion method is proposed that integrates Sun’s petrophysical model and Russell’s poroelasticity theory to derive a decoupled reflection coefficient formula for fluid bulk modulus , porosity , and pore structure parameter . A two-layer theoretical geological model is used for forward modeling simulation to verify the accuracy of the newly derived formula and assess the impact of porosity and the pore structure parameter on the seismic response. Then, an improved decorrelation method is employed to eliminate correlations between multiple parameters, ensuring stability and robustness during the inversion process. In the practical application in a specific exploration area within the Sichuan Basin, this study first demonstrates through digital core experiments that the pore structure parameter has great influence on elastic moduli. Furthermore, the results based on the actual data show that the fluid bulk modulus as the proposed fluid indicator eliminates the interference of the pore structure ( and ) and can provide more accurate and reliable fluid distribution than Russell’s fluid factor .
2 Methods
Establishing a direct mathematical relationship between theoretical or empirical petrophysical models and expressions of seismic reflection coefficients is the key to quantitative characterization of fluid prediction involving the pore structure.
2.1 Poroelasticity theory
First, one of the widely used quantitative formulas of seismic response for fluid identification is Biot–Gassmann theory (), which elucidates the pore/fluid interaction in homogeneous porous media and defines the functional relationship between seismic velocity and fluid. Its corresponding formula is shown as follows:where and represent the bulk moduli for the saturated and dry porous rocks, and denote the shear moduli for the saturated and dry porous rocks, represents the Biot coefficient, and signifies the modulus of pressure that drives water into strata without changing its volume. ) derived the Gassmann formula based on the pore-elasticity theory and found that the fluid/porosity term was related to multiple petrophysical parameters, including porosity, which could be quantitatively characterized by saturated elastic parameters:where and denote the bulk moduli of the fluid and saturated rock frame and represents the porosity. Therefore, formula (2) can be reformulated as follows:
In addition, by substituting the fluid factor into Aki–Richard approximation, Russell proposed a formula for the seismic response that can directly obtain by inversion () as follows:where denotes the incidence angle, signifies the density, and and represent the velocity ratios of P and S wave in saturated and dry rocks, respectively.
2.2 Reflection coefficient formula for decoupling porosity and pore structure parameter
Second, in terms of petrophysics, ) proposed a petrophysical model based on Biot theory, introducing the pore structure parameter to effectively depict the effect of the pore shape on seismic velocity. Given its inclusion of both the pore structure parameter and porosity , Sun’s model exhibits superior applicability for characterizing reservoirs with a complex pore structure compared to conventional petrophysical models:
It is worth noting that the pore structure parameter represents the pore shape, particularly the aspect ratio , which does not fall within a numerical range from 0 to 1. Through digital core simulation, , clarified their non-linear relationship: when and when .
According to the derivation by ), it can be concluded that when is extremely small (). Additionally, in tight sandstone or carbonate with low porosity, exhibits a significant influence on the seismic response, which greatly contributes to heterogeneity. Therefore, formula (6) can be simplified as follows:
It can be found by comparing formula 9, 4that
The exponential term in formula 10 poses challenges for numerical analysis, and therefore, it can be expanded by the Taylor series as follows:
The selection of different orders in approximation (12) can be tailored to meet different accuracy requirements. To achieve the decoupling of from and by direct inversion, formula 12 in the first-order approximation is substituted into formula 5 to obtain a novel fluid prediction formula as follows:where represents the dry rock matrix term () and denotes the structural density term.
3 Accuracy analysis
For verifying the rationality of the proposed reflection coefficient in formula 13 for decoupling the pore structure parameter and porosity simultaneously, a two-layer geological model based on the logging data from carbonate reservoirs is designed to analyze its accuracy. The model parameters are shown in Table 1. The exact Zoeppritz formula, Aki–Richards approximation, and the proposed formula are used to calculate the reflection coefficient, and their comparison is shown in Figure 1. The reflection coefficient obtained from the new formula exhibits good agreement with that from the exact Zoeppritz formula for incidence angle below 35°, while maintaining similar accuracy to the A&K approximation for incidence angle above 35°. Therefore, Figure 1 shows that there is no loss of computational accuracy in the proposed formula of separating the pore structure despite the increase in input parameters.
TABLE 1
| (GPA) | (GPA) | (m/s) | (m/s) | (kg/m3) | (%) | (−) | |
|---|---|---|---|---|---|---|---|
| Layer 1 | 42.5 | 13.5 | 4,000 | 2,350 | 2,450 | 6 | 3 |
| Layer 2 | 18 | 14.5 | 4,175 | 2,420 | 2,460 | 8 | 6.5 |
Parameters of the two-layer strata model.
kf is fluid bulk modulus, μ is shear modulus, vp is P-wave velocity, vs is S-wave velocity, ρ is density, φ is porosity, γ is pore structure parameter.
FIGURE 1
4 Importance of the pore structure
It should be noted that there are five terms to be solved in the proposed approximation (13), including porosity and the pore structure parameter . Since these two parameters are not commonly used in a seismic inversion, it is necessary to confirm their considerable influence on the AVO reflection coefficient.
A two-layer model with the same basic elastic parameters (; ; ) is built to reveal the effect of the pore structure on the reflection coefficient by varying only or , as shown in Figure 2. Figure 2A shows that when is fixed at 5 and changes from 1% to 10%, the reflection coefficient exhibits significant changes. Similarly, in Figure 2B, when is fixed at 6% and changes from 2 to 8, corresponding to the transition from micro-fractures to vugular pores in pore types, noticeable variations in the reflection coefficient are observed. The forward simulation results indicate that both and have a considerable impact on the seismic response, akin to conventional elastic parameters. Upon comparing Figure 2A with (b), it is observed that at low values of , the seismic response of shows similarities with that of . The consistent variation trend of reflection coefficient in both figures suggests that an increase in the complexity of the pore structure leads to a more pronounced seismic response. Therefore, the influence of on the variation of seismic reflection should not be neglected, necessitating the consideration of both and as pivotal factors in fluid prediction.
FIGURE 2
5 Inversion for decoupling the pore structure
5.1 Bayesian framework
The solution of the proposed method, as indicated by formula 13, involves five parameters to be solved for. This places higher demands on the stability and accuracy of the inversion compared to the traditional AVO method, which typically require only three parameters. The proposed formula requires the association of (where takes an integer greater than or equal to five) equations, with each equation corresponding to an incidence angle. Assuming that each seismic trace has incidence angles and sampling points, let , , , , and ; consequently, formula 13 can be rewritten as follows:where , , , , , and denotes a unit matrix of .
For simplicity of expression, formula 14 can be reformulated as follows:where , , and .
The solution of formula 15 provides the five parameters; however, it suffers from severe ill-posedness in practical applications. Therefore, based on the Bayesian framework, the issue of well-posedness can be effectively improved by introducing the prior information of the model parameters into the regularization of inversion. Thus, the posterior probability density distribution of model parameters iswhere signifies the prior distribution of the model parameters, represents the likelihood function that characterizes the noise distribution of seismic gathers, and denotes a constant that can be ignored if the posterior probability distribution function remains unchanged.
The prestack seismic gathers in actual data contain a certain level of noise. Assuming that the noise is uncorrelated with a mean value of 0, following the likelihood function of the Gaussian distribution, thus, the likelihood function of noise can be expressed as follows:where and represent the covariance and the mean square error of P-wave noise, respectively.
The prior distribution of the model parameters can be either univariate or multivariate. Considering that five involved parameters, which have certain correlation, need to be obtained by the proposed method, the selection of multivariate distribution can effectively mitigate the ill-conditioning issue. While Gaussian distribution only provides consistent weighting coefficients, which may impact the sparsity of inversion results, Cauchy distribution can yield non-consistent weighting coefficients with sparsity effects and greater geological significance. Therefore, the multivariate Cauchy distribution is adopted as the prior distribution for the model parameters, with its specific formula given as follows:where signifies the correlation matrix, which can be obtained by maximum expectation estimation. denotes the matrix of , and its expression is as follows:
By substituting formulas 17, 18 into Formula 16, the posterior probability density of model parameters can be obtained as follows:
After the algebraic transformation of Formula 20, the objective function should be as follows:where determines the sparsity of inversion results, with larger values leading to sparser outcomes. Conversely, if is set too small, it may result in a distorted inversion.
5.2 Decorrelation method
The statistical correlation between elastic parameters is widely recognized as a major factor contributing to the instability of the prestack inversion. The proposed method involves five parameters, out of which three parameters (fluid bulk modulus , porosity , and pore structure parameter ) are further decomposed by the fluid factor , resulting in a closer relationship among them, thus leading to a substantial increase in inversion instability. Evidently, relying solely on the prior probability distribution scheme is insufficient for eliminating the parameter correlations. To address this issue, an improved decorrelation method () is adopted to obtain parameters with the lowest possible correlation, which uses the variance matrix and two linear transformations based on Chen’s previous work () to convert relevant data into “white data.”
The decorrelation method is illustrated using a three-parameter dataset as an example, since it involves five parameters in the proposed inversion, which is not convenient for figuring. As shown in Figure 3A, the sequences x, y, and z exhibit significant correlation, and their covariance matrix is expressed as follows:where , , and represent the variance of , , and , respectively, and , , and denote the covariance of the three parameters. The singular value decomposition of the covariance matrix iswhere signifies the eigenvector matrix and denotes the eigenvalue matrix. Upon applying to , the resulting sequences are shown in Figure 3B. However, despite this transformation, the newly obtained data still have a weak correlation. To achieve further decorrelation, we apply to , resulting in . Figure 3C shows that can be referred to as “white data,” where are almost completely unrelated.
FIGURE 3
According to the aforementioned method, one of the covariance matrices of , , , , and can be decomposed as
The covariance matrix can be extended to two sparse eigenvector matrix and eigenvalues matrix by considering time samples. Therefore, the transformation observed in Eq. 21 can be represented as follows:where and . The reflection coefficients of , , , , and can be obtained by using the iterative reweighed least-squares to solve the objective function (25), followed by obtaining the five parameters through the trace integral.
6 Synthetic data test
The proposed inversion of decoupling porosity and pore structure parameter involves a total of five parameters, and thus, it is imperative to verify its feasibility and anti-noise capabilities. A set of well data from carbonates in the actual working area of the Sichuan Basin was selected for testing purposes. First, prior to inversion, the measured data underwent the Backus averaging process () which transformed the data from log scale to seismic scale, followed by time–depth conversion to convert it from depth domain to time domain. Second, the reflection coefficients of the well data at various sampling times and angles (5°–60°) are obtained by forward simulation based on the exact Zoeppritz equation, subsequently convolved with a 30 Hz Ricker wavelet to generate synthetic seismic data, as shown in Figure 4A. Third, Gaussian random noise with the signal-to-noise ratio (S/N) of 2 was added to the synthetic record, as shown in Figure 4B.
FIGURE 4
Figures 5A, B show the inversion results without and with noise, respectively, displaying the three direct inversion parameters of fluid bulk modulus , which characterize reservoirs, porosity , and pore structure parameter . Additionally, the Russell fluid factor is calculated indirectly using formula 10 or (12) following the proposed inversion. In the absence of noise, both the direct and indirect inversions yield consistent results with well curves. However, due to a weaker seismic response than the fluid bulk modulus , porosity , and pore structure parameter exhibit slightly lower resolution. At S/N = 2, the direct inversion results (, , and ) show a slight decrease in accuracy compared to noise-free conditions, while errors became more apparent in indirect inversion’s due to error accumulation from indirect calculation. Nevertheless, even with added noise during the inversion process, the inversion results still maintain a similar trend as observed in the well data, indicating a good stability of the proposed method for practical applications.
FIGURE 5
7 Actual data application
The well data from an exploration area in the Sichuan Basin is utilized to validate the feasibility and applicability of the proposed prestack inversion technique for decoupling both pore structure parameter and porosity simultaneously. The target reservoir comprises dolomitic gas-bearing carbonate with low porosity but an intricate pore structure, which constitutes one of the influential factors affecting fluid prediction.
7.1 Digital core analysis
Before fluid detection, component analysis and simulation of different pore shapes were conducted on a total of seven cores from three wells in various formations to construct multi-component digital core models for evaluating the influence of the pore structure in this working area. Figure 6A shows a gray-scale slice obtained through a CT scan depicting the actual core. Based on the analysis of the actual core, the multi-component digital core was constructed, as shown in Figure 6B. The digital core comprised two dominant minerals: dolomite (light-colored) and quartz (dark-colored), and the dark-colored component also containing 7% kaolinite. As numerous previous studies have demonstrated the effect of porosity on elastic parameters, such as those by ) and ), the effect of the pore structure parameter on elastic parameters within different components in the digital cores is solely focused on being discussed. The relationship between pore structure parameter and saturated bulk modulus or shear modulus of the digital core is shown in Figures 7A, B. The of each mineral component is greatly affected by when , with the of dolomite experiencing a change of over 40% and quartz approximately 32%, in particular. However, the effect of on is negligible and can be disregarded, as all three mineral components vary by less than 10% when . However, in the simulation of filling oil or gas, the measured fluid bulk modulus remains constant regardless of variations in the three mineral components or changes in the pore structure parameter . Therefore, if fluid identification does not take the impact of the pore structure parameter into account, its prediction results will exhibit significant multi-solution possibilities and uncertainties.
FIGURE 6
FIGURE 7
7.2 Fluid prediction
Reservoir analysis is initially performed on the well data, using the same well in the anti-noise experiment, as shown in Figure 8. Among the three logged carbonate gas-bearing reservoirs (highlighted by red and blue boxes in Figure 8), Russell fluid factor exhibited conspicuous manifestation solely within the second reservoir, and fluid bulk modulus demonstrated consistent characteristics in all the three reservoirs. According to Eqs 3, 10, Russell fluid factor serves as an integrated indicator of fluid within the saturated porous rock, whereas fluid bulk modulus acts as a fluid indicator independent of porosity and pore structure parameter . Therefore, by analyzing the curves of and , it becomes discernible as to which factor affects the prediction outcomes: The red-boxed reservoir is primarily influenced by , whereas the blue-boxed one is affected by . Notably, the third reservoir represents a thinly interbedded formation controlled by both and , thereby indicating its highly heterogeneous nature.
FIGURE 8
Second, both the proposed inversion and Russell’s inversion are calculated for fluid prediction. Additionally, indirectly obtained from the results of the proposed inversion by formula 10 is also obtained for comparison, as shown in Figure 9. The outcomes of fluid bulk modulus (Figure 9A), porosity (Figure 9B), and pore structure parameter (Figure 9C) closely match the data from the test well. Moreover, the profile of direct (Figure 9E) coincides with that of indirect (Figure 9D), indicating the stability and reliability of the proposed inversion. The result of (Figure 9A), which is decoupled from and , exhibits a more pronounced distribution pattern in reservoirs compared to that of . Notably, between them, there exists a significant disparity in fluid prediction for the third reservoir—the thin interbedded gas-bearing formation. By magnifying the red box near the well in Figures 9A, E,as shown in Figure 10, it can be observed that (Figure 10A) displays better consistency with verified gas production tests within the thin interbedded gas reservoirs than does (Figure 10B), which shows that after eliminating the influence of and through inversion, the prediction of accurately reflects reservoir distribution.
FIGURE 9
FIGURE 10
8 Conclusion
The pore structure parameter , like porosity , is one of the dominant factors affecting fluid prediction in heterogeneous reservoirs. By combining Russell’s poroelasticity theory and Sun’s petrophysical model, a new reflection coefficient formula is proposed, which decouples fluid bulk modulus from porosity and pore structure parameter , thus eliminating the interference of the pore-related factors on reservoir prediction. Theoretical experiments demonstrate that the proposed inversion retains comparable computational accuracy to the conventional A&K approximation and has good anti-noise ability under the condition of five involved parameters. Both theoretical model and digital core analyses reveal that the pore structure parameter can have a significant impact on the prediction of heterogeneous reservoirs with low porosity. In the example from the Sichuan Basin, compared with Russell fluid factor , fluid bulk modulus from the proposed inversion method provides more accurate and distinct reservoir distribution in carbonate gas-bearing strata. However, there are still two prerequisites for the application of the proposed inversion. First, as a less commonly used parameter, the pore structure parameter needs to be obtained from sufficient logging data in the working area under generally industrial process. The method used to obtain the pore structure parameter in this paper is the empirical method proposed by , , which is still inadequate in terms of efficiency and verifiability. Second, in order to ensure the robustness of the proposed inversion method with five parameters, high-quality seismic data are also indispensable.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding author.
Author contributions
RH: methodology and writing-original draft. EW: conceptualization and project administration. WW: methodology and visualization. GY: data curation and formal analysis. XZ: validation and writing-review and editing. WZ: investigation and formal analysis. DH: conceptualization and investigation.
Funding
This work was funded by the Key Technology Research Project of CNPC, grant number 2021ZG03, and the Science and Technology Project of Petrochina, grant number 2022DJ8004.
Conflict of interest
Authors RH, EW, WW, GY, XZ, WZ, and DH were employed by PetroChina.
Publisher’s note
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.
References
1
BackusG.GilbertF. (1968). The resolving power of gross Earth data. Geophys. J. Int.16 (2), 169–205. 10.1111/j.1365-246x.1968.tb00216.x
2
BerrymanJ. G. (1999). Origin of Gassmann’s equations. GEOPHYSICS64 (5), 1627–1629. 10.1190/1.1444667
3
BiotM. A. (1941). General theory of three-dimensional consolidation. J. Appl. Phys.12 (2), 155–164. 10.1063/1.1712886
4
ChenJ. J.YinX. Y. (2007). Three-parameter AVO waveform inversion based on Bayesian theorem. Chin. J. Geophys.50 (4), 1251–1260. 10.1007/s11442-007-0020-2
5
DengJ.ZhouH.WangH.ZhaoJ.WangS. (2015). The influence of pore structure in reservoir sandstone on dispersion properties of elastic waves. Chin. J. Geophys.58 (9), 3389–3400. 10.6038/cjg20150931
6
DuB.YangW.ZhangJ.YongX.GaoJ.LiH. (2019). Matrix-fluid decoupling-based joint PP-PS wave seismic inversion for fluid identification. GEOPHYSICS84 (3), R477–R487. 10.1190/geo2017-0376.1
7
FanX.ZhangG.ZhangJ. (2019). Prediction method of pore structure parameters of tight sandstone. Seg. Tech. Program Expand. Abstr.89, 3678–3682. 10.1190/segam2019-3199837.1
8
GassmannF. (1951). Über die elastizitat poroser medien: vierteljahrschrift der naturforschenden. Gesellshaft96, 1–23.
9
HanD.BatzleM. L. (2004). Gassmann’s equation and fluid‐saturation effects on seismic velocities. GEOPHYSICS69 (2), 398–405. 10.1190/1.1707059
10
HeX.HeZ.WangX.XiongX.JiangL. (2012). Rock skeleton models and seismic porosity inversion. Appl. Geophys.9 (3), 349–358. 10.1007/s11770-012-0345-1
11
HeX.LinK.ZhangZ.HeZ. (2018). Quantitative comparison of the influence of porosity and pore structure on reservoir characteristics. Geophys. Prospect. Petroleum57 (2), 179–185.
12
JiangL.WenX.ZhouD.HeZ.HeX. (2012). The constructing of pore structure factor in carbonate rocks and the inversion of reservoir parameters. Appl. Geophys.9 (2), 223–232. 10.1007/s11770-012-0333-5
13
KachanovM.TsukrovI.ShafiroB. (1994). Effective moduli of solids with cavities of various shapes. Appl. Mech. Rev.47 (1S), S151–S174. 10.1115/1.3122810
14
KriefM.GaratJ.StellingwerffJ.VentreJ. (1990). A petrophysical interpretation using the velocities of P and S waves. Log. Anal.31, 355–369.
15
LiH.ZhangJ.PanH.GaoQ. (2021). Nonlinear simultaneous inversion of pore structure and physical parameters based on elastic impedance. Sci. China Earth Sci.64 (6), 977–991. 10.1007/s11430-020-9773-8
16
RussellB. H.GrayD.HampsonD. P. (2011). Linearized AVO and poroelasticity. GEOPHYSICS76 (3), C19–C29. 10.1190/1.3555082
17
RussellB. H.HedlinK.HiltermanF. J.LinesL. R. (2003). Fluid‐property discrimination with AVO: A biot-gassmann perspective. GEOPHYSICS68 (1), 29–39. 10.1190/1.1543192
18
SunR.YinX.WangB.ZhangG. (2016). A direct estimation method for the Russell fluid factor based on stochastic seismic inversion. Chin. J. Geophys.59 (3), 1143–1150. 10.6038/cjg20160334
19
SunW.HeZ.LiY.ZhangF.ZhouY. (2015). A new factor of fluid-filled fractures and its application. Chin. J. Geophys.58 (7), 2536–2545. 10.6038/cjg20150728
20
SunY. (2000). Core-log-seismic integration in hemipelagic marine sediments on the eastern flank of the Juan de Fuca Ridge. Proc. Ocean Drill. Program, Sci. Results168, 21–23. 10.2973/odp.proc.sr.168.009.2000
21
SunY.GoldbergD. (1997). Effects of aspect ratio on wave velocities in fractured rocks. Seg. Tech. Program Expand. Abstr.67, 925–928. 10.1190/1.1886170
22
WangF.BianH.ZhangY.DuanC.ChenG. (2016). Hilbert-Huang transform combined with smoothed pseudo Wigner-Ville time-frequency distribution to identify reservoir fluid properties. Geophys. Prospect. Petroleum55 (6), 851–860. 10.3969/j.issn.1000-1441.2016.06.010
23
WangW.LiG.YuX.YangW.WangW.FengL.et al (2017). A meta-analysis of adenosine A2A receptor antagonists on levodopa-induced dyskinesia in vivo. Seg. Tech. Program Expand. Abstr.87, 702–706. 10.3389/fneur.2017.00702
24
YinX.ZhangS. (2014). Bayesian inversion for effective pore-fluid bulk modulus based on fluid-matrix decoupled amplitude variation with offset approximation. GEOPHYSICS79 (5), R221–R232. 10.1190/geo2013-0372.1
25
ZhangT.SunY.TianJ.LuL.QinF.ZhaoX.et al (2018b). Two-parameter prestack seismic inversion of porosity and pore-structure parameter of fractured carbonate reservoirs: part 2 - applications. Interpretation6 (4), SM9–SM17. 10.1190/int-2018-0019.1
26
ZhangT.SunY. (2018a). Two-parameter prestack seismic inversion of porosity and pore structure parameter of fractured carbonate reservoirs: part 1 — methods. Interpretation6, 1–SM8. 10.1190/int-2017-0219.1
27
ZhaoJ.PanJ.HuY.LiJ.LiuX.LiC.et al (2021a). Digital rock physics-based studies on effect of pore types on elastic properties of carbonate reservoir Part 1: imaging processing and elastic modelling. Chin. J. Geophys.64 (2), 656–669. 10.6038/cjg2021O0228
28
ZhaoJ.PanJ.HuY.LiJ.LiuX.LiC.et al (2021b). Digital rock physics-based studies on effect of pore types on elastic properties of carbonate reservoir Part 2: pore structure factor characterization and inversion of reservoir. Chin. J. Geophys.64 (2), 670–671. 10.23736/S0026-4806.20.06558-1
29
ZimmermanR. W. (1986). Compressibility of two-dimensional cavities of various shapes. J. Appl. Mech.53 (3), 500–504. 10.1115/1.3171802
30
ZongZ.WangY. (2019). Amplitude variation with angle inversion for fluid discrimination with the consideration of squirt flow. Seg. Tech. Program Expand. Abstr.89, 734–738. 10.1190/segam2019-3216431.1
31
ZongZ.YinX.WuG. (2012). AVO inversion and poroelasticity with P- and S-wave moduli. GEOPHYSICS77 (6), N17–N24. 10.1190/geo2011-0214.1
32
ZongZ.YinX.WuG.WuZ. (2015). Elastic inverse scattering for fluid variation with time-lapse seismic data. GEOPHYSICS80 (2), WA61–WA67. 10.1190/geo2014-0183.1
Summary
Keywords
fluid prediction, pore structure parameter, fluid factor, porosity, carbonate
Citation
He R, Wang E, Wang W, Yan G, Zheng X, Zhao W and He D (2023) Bayesian fluid prediction by decoupling both pore structure parameter and porosity. Front. Earth Sci. 11:1269597. doi: 10.3389/feart.2023.1269597
Received
30 July 2023
Accepted
29 August 2023
Published
15 September 2023
Volume
11 - 2023
Edited by
Xintong Dong, Jilin University, China
Reviewed by
Xiangbo Gong, Jilin University, China
Shiqi Dong, Changjiang River Scientific Research Institute (CRSRI), China
Updates
Copyright
© 2023 He, Wang, Wang, Yan, Zheng, Zhao and He.
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: Run He, he_run732@126.com
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.