Next Article in Journal
Optimum Organic Rankine Cycle Design for the Application in a CHP Unit Feeding a District Heating Network
Previous Article in Journal
An MPC Scheme with Enhanced Active Voltage Vector Region for V2G Inverter
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

High-Order AVO Inversion for Effective Pore-Fluid Bulk Modulus Based on Series Reversion and Bayesian Theory

1
SINOPEC Exploration Production Research Institute, Haidian District, Beijing 100083, China
2
College of Geophysics, China University of Petroleum—Beijing, Changping District, Beijing 102249, China
3
College of Petroleum, China University of Petroleum—Beijing at Karamay, Karamay, Xinjiang 834000, China
*
Author to whom correspondence should be addressed.
Energies 2020, 13(6), 1313; https://doi.org/10.3390/en13061313
Submission received: 31 January 2020 / Revised: 4 March 2020 / Accepted: 6 March 2020 / Published: 12 March 2020
(This article belongs to the Section H: Geo-Energy)

Abstract

:
Pore-fluid identification is one of the key technologies in seismic exploration. Fluid indicators play important roles in pore-fluid identification. For sandstone reservoirs, the effective pore-fluid bulk modulus is more susceptible to pore-fluid than other fluid indicators. AVO (amplitude variation with offset) inversion is an effective way to obtain fluid indicators from seismic data directly. Nevertheless, current methods lack a high-order AVO equation for a direct, effective pore-fluid bulk modulus inversion. Therefore, based on the Zoeppritz equations and Biot–Gassmann theory, we derived a high-order P-wave AVO approximation for an effective pore-fluid bulk modulus. Series reversion and Bayesian theory were introduced to establish a direct non-linear P-wave AVO inversion method. By adopting this method, the effective pore-fluid bulk modulus, porosity, and density can be inverted directly from seismic data. Numerical simulation results demonstrate the precision of our proposed method. Model and field data evaluations show that our method is stable and feasible.

1. Introduction

Pore-fluid identification remains a key part of hydrocarbon reservoir exploration. Fluid indicators are crucial to distinguishing pore-fluid. The first fluid indicator is defined as the weighting difference of P- and S-wave velocities (Smith and Gidlow, 1987) [1]. Hence, Fatti et al. (1994) proposed the weighting difference of P- and S-wave reflection coefficients as a fluid indicator [2]. As the recognition technology developed, the definition of the fluid indicator expanded into the impedance domain. Goodway et al. (1997) showed the advantages of Lamé parameters in fluid discrimination while lithologic parameters remained unchanged [3]; however, the identification accuracy of Lamé parameters will decrease if lithologic parameters change. Quakenbush et al. (2006) defined a fluid indicator with a bulk modulus, which was more sensitive to pore-fluid and more independent of a rock matrix than the Lamé parameters [4]. Russell et al. (2011) generalized the fluid indicator to a Gassmann fluid factor and enhanced the ability to distinguish pore-fluid [5]. Most conventional fluid indicators are a combination of density and velocity. However, these fluid indicators can often provide uncertain results when determining fluids because their sensitivity is conditional on the mixed effects of pore-fluid and rock matrices (Zhang et al., 2010) [6]. To improve the accuracy of fluid identification, Yin and Zhang (2014) presented an effective pore-fluid bulk modulus as a fluid indicator. It has been proven that the effective pore-fluid bulk modulus is advantageous in identifying reservoirs because it is only related to pore-fluid [7].
AVO inversion is a useful tool to extract fluid indicators from seismic data directly. AVO inversion methods can be divided into linear and nonlinear methods. The former uses linear AVO approximations to establish inversion equations, while the latter uses exact expressions or high order AVO approximations. Yin and Zhang (2014) deduced a linear four-term AVO approximation and proposed a Bayesian AVO inversion method for an effective pore-fluid bulk modulus. By using this method, the effective pore-fluid bulk modulus could be inverted from seismic data directly [7]. However, the four terms of this approximation are not completely independent, and this method needs four angle gathers. Yin et al. (2018) put forward a new linear three-term AVO approximation and enhanced the stability of inversion [8]. Nevertheless, we must analyze the exponential relationship between porosity and shear modulus when we invert field data by employing this approximation. Linear approximations have high precision under small incident angles and parameter differences at both sides of an interface; however, they will suffer a loss of accuracy with increasing incident angles and parameter differences. At present, neither exact nor nonlinear equations exist for an effective pore-fluid bulk modulus.
In this paper, an effective pore-fluid bulk modulus is considered as a fluid indicator. Based on Zoeppritz equations and the Biot–Gassmann theory, we first deduce a high-order AVO approximation for the effective pore-fluid bulk modulus. Consequently, a direct nonlinear AVO inversion method is formulated for the fluid indicator in a series reversion and Bayesian theory framework. Model and field data tests follow.

2. Theory and Method

2.1. Sensitivity Analysis of the Effective Pore-Fluid Bulk Modulus

Various fluid indicators have been proposed to identify pore-fluid. In this section, we analyze the sensitivity of several commonly used fluid indicators to detect pore-fluid changes in sandstone reservoirs. By using the three types of sandstone model displayed in Table 1 [9], the P-wave modulus, Gassmann fluid factor, and effective pore-fluid bulk modulus are computed. The calculation results are shown in Table 2. To analyze the sensitivity of these three fluid indicators, we define the relative variation as an identified parameter and visualize this identified parameter in the bar graph (Figure 1). The larger the relative variation value is, the more sensitive the fluid indicator is to the change of pore-fluid. The relative variations of the effective pore-fluid bulk modulus are the largest in all of the three pictures in Figure 1. This demonstrates that the effective pore-fluid bulk modulus has more advantages in identifying pore-fluid than the other two fluid indicators did.

2.2. P-Wave High-Order AVO Approximations for the Effective Pore-Fluid Bulk Modulus

An effective pore-fluid bulk modulus has shown a clear advantage as a fluid indicator for identifying pore-fluid. AVO inversion is a useful technology for extracting an effective pore-fluid bulk modulus from seismic data directly, though a high-order approximation for inversion is nonexistent. Here, a new estimation is deduced based on the Zoeppritz equations and Biot–Gassmann theory. The Zoeppritz equations are the foundation of AVO methodology [10]. Based on the relationships among seismic amplitudes, this equation can be expressed as (Keys, 1989) [11]
[ X 1 B 2 X 2 C X 1 D 2 X 2 1 X 2 B X 1 C 2 X 2 D X 2 B 2 X 1 X 2 B ( 1 2 B 2 X 2 ) 2 A D 2 X 1 C 2 X 2 A D ( 1 2 D 2 X 2 ) ( 1 2 B 2 X 2 ) 2 B 2 X 1 B 2 X 2 A C ( 1 2 D 2 X 2 ) 2 A D 2 X 1 D 2 X 2 ] [ R P P R P S T P P T P S ] = [ X 1 X 2 2 B 2 X 1 X 2 1 2 B 2 X 2 ] ,
where X represents sin θ , with θ as the incident angle; R P P and R P S are P- and S-wave reflection coefficients, respectively; and T P P and T P S are P- and S-wave transmission coefficients, respectively.
For the incident P-wave case, A , B , C , and D are defined as (Kim and Innanen, 2013) [12]
A = ( 1 Δ ρ ρ ¯ ) 1 ,
B = γ s a t 1 ,
C = { ( 1 Δ ρ ρ ¯ ) [ γ d r y 2 γ s a t 2 ( 1 Δ μ μ ¯ ) 1 + ( 1 γ d r y 2 γ s a t 2 ) ( 1 Δ f f ¯ ) 1 ] } 1 2 ,
D = [ γ s a t 2 ( 1 Δ μ μ ¯ ) 1 ( 1 Δ f f ¯ ) ] 1 2 ,
where γ d r y and γ s a t denote the velocity ratios of P- and S-wave in dry rock and saturated rock, respectively. Δ f / f ¯ , Δ μ / μ ¯ and Δ ρ / ρ ¯ are calculated by
Δ f f ¯ = 2 ( f i + 1 f i ) ( f i + 1 + f i ) ,
Δ μ μ ¯ = 2 ( μ i + 1 μ i ) ( μ i + 1 + μ i ) ,
Δ ρ ρ ¯ = 2 ( ρ i + 1 ρ i ) ( ρ i + 1 + ρ i ) ,
where f and μ represent the Gassmann fluid factor and shear modulus, respectively, which can be calculated by (Russell et al., 2011) [5]
f = ρ α 2 γ d 2 μ ,
μ = ρ β 2 ,
where ρ is density, and i and i + 1 are the upper and lower layers of the interface, respectively.
According to the Gassmann equation, the Gassmann fluid factor can be derived as (Han and Batzle, 2003) [13]
f = G K f ,
where K f is the effective pore-fluid bulk modulus, and
G = ( 1 K n ) 2 ϕ ,
K n = K d r y K s ,
where ϕ , K s , and K d r y stand for the porosity of porous rock, the bulk modulus of solid grain, and the bulk modulus of dry porous rock, respectively.
Based on the critical porosity model, Nur et al. (1998) provide the following equations [14]:
K d r y = K s ( 1 ϕ ϕ c ) ,
μ d r y = μ s ( 1 ϕ ϕ c ) ,
where ϕ c , μ s , and μ d r y denote critical porosity, the shear modulus of solid grain, and the shear modulus of dry porous rock, respectively.
For the same reservoir, we assume that critical porosity is constant. Replacing μ with μ d r y , and substituting Equations (5) and (7a,b) into Equation (2a–d), we have
A = ( 1 Δ ρ ρ ¯ ) 1 ,
B = γ s a t 1 ,
C = { ( 1 Δ ρ ρ ¯ ) [ ( γ d r y 2 γ s a t 2 ) ( 1 γ ϕ Δ ϕ ϕ ¯ ) 1 + ( 1 γ d r y 2 γ s a t 2 ) ( 1 Δ ϕ ϕ ¯ Δ K f K ¯ f ) 1 ] } 1 2 ,
D = [ γ s a t 2 ( 1 γ ϕ Δ ϕ ϕ ¯ ) 1 ( 1 Δ ρ ρ ¯ ) ] 1 2 .
By substituting Equation (8a–d) into (1) and applying the Taylor expansion to Δ K f / K ¯ f , Δ ϕ / ϕ ¯ and Δ ρ / ρ ¯ , the P-wave reflection coefficient is deduced as follows:
R p p = R p p ( 1 ) + R p p ( 2 ) + R p p ( 3 ) + ,
where, R p p ( i )     ( i = 1 , 2 , 3 , ) are the different order terms of P-wave reflection coefficients, involving the effective pore-fluid bulk modulus, porosity, and density. The ith-order approximation is written as R p p = R p p ( 1 )    + + R p p ( i )     ( i = 1 , 2 , 3 , ) . Using the method of numerical analysis, we study the effect of the expansion order number on the accuracy of the deduced approximation. For small incident angle and parameter differences across the interface, the relative errors of different order approximations are all less than 1%. With an increase of the incident angle and parameter differences across the interface, the accuracy of the first- and second-order approximations decreases heavily. The relative error of the third-order approximation increases but is still less than 1%. When the expansion order number is greater than 3, the accuracy of the deduced approximation increases slightly, but the calculation efficiency decreases heavily. We considered both accuracy and efficiency, and decided to use the third-order approximation, which is expressed as
R p p = R p p ( 1 ) + R p p ( 2 ) + R p p ( 3 ) ,
R p p ( 1 ) = W 11 Δ K f K ¯ f + W 12 Δ ϕ ϕ ¯ + W 13 Δ Δ ρ ρ ¯ ,
R p p ( 2 ) = W 21 ( Δ K f K ¯ f ) 2 + W 22 ( Δ ϕ ϕ ¯ ) 2 + W 23 ( Δ ρ ρ ¯ ) 2 + W 24 ( Δ K f K ¯ f ) ( Δ ϕ ϕ ¯ ) + W 25 ( Δ K f K ¯ f ) ( Δ ρ ρ ¯ ) + W 26 ( Δ ϕ ϕ ¯ ) ( Δ ρ ρ ¯ ) ,
R p p ( 3 ) = W 30 ( Δ K f K ¯ f ) 3 + W 31 ( Δ ϕ ϕ ¯ ) 3 + W 32 ( Δ ρ ρ ¯ ) 3 + W 33 ( Δ K f K ¯ f ) 2 ( Δ ϕ ϕ ¯ ) + W 34 ( Δ K f K ¯ f ) 2 ( Δ ρ ρ ¯ ) + W 35 ( Δ ϕ ϕ ¯ ) 2 ( Δ K f K ¯ f ) + W 36 ( Δ ϕ ϕ ¯ ) 2 ( Δ ρ ρ ¯ ) + W 37 ( Δ ρ ρ ¯ ) 2 ( Δ K f K ¯ f ) + W 38 ( Δ ρ ρ ¯ ) 2 ( Δ ϕ ϕ ¯ ) + W 39 ( Δ K f K ¯ f ) ( Δ ϕ ϕ ¯ ) ( Δ ρ ρ ¯ ) .
Using the sand model in Table 3 (Yin and Zhang, 2014) [7], we compare the accuracy of these approximations. P-wave reflection coefficients are computed through the Zoeppritz equations [10], Aki–Richards approximation (Aki and Richards, 1980) [15], and different order approximations, respectively. As shown in Figure 2a, the curves of the Aki–Richards approximation and different order approximations all are close to the Zoeppritz curve. Small elastic parameter perturbations across reflectors are caused only by pore-fluid (Model 1). In Figure 2b,c, note that the AVO curves computed from the Aki–Richards approximation and different order approximations are close to the Zoeppritz curve for incident angles less than 30°. However, for larger incident angles the curve of the deduced third-order approximation is the closest to the Zoeppritz curve. This is because the elastic parameter perturbations across reflectors are caused by only porosity or both pore-fluid and porosity (Models 2 and 3).
Next, we analyze the information contained in reflections and the stability of inversion for the first-order approximation. If the second- and higher-order terms of Equation (9) are neglected, it becomes a linear approximation. The P-wave reflection coefficient is written as
R p p R p p ( 1 ) = [ ( 1 4 1 4 γ d r y 2 γ s a t 2 ) sec 2 θ ] Δ K f K ¯ f + [ 2 1 γ s a t 2 γ ϕ sin 2 θ + ( 1 4 1 4 γ d r y 2 γ s a t 2 + 1 4 γ d r y 2 γ s a t 2 γ ϕ ) sec 2 θ ] Δ ϕ ϕ ¯ + ( 1 4 sec 2 θ + 1 2 ) Δ ρ ρ ¯ .
In this equation, the effective pore-fluid bulk modulus, porosity, and density have different proportions of information in their reflection coefficients. As a result, the stability of inversion for these three elastic parameters is different. Equation (11) can be expressed as
R p p = [ ( 1 4 1 4 γ d r y 2 γ s a t 2 ) sec 2 θ [ 2 1 γ s a t 2 γ ϕ sin 2 θ + ( 1 4 1 4 γ d r y 2 γ s a t 2 + 1 4 γ d r y 2 γ s a t 2 γ ϕ ) sec 2 θ ]   1 4 sec 2 θ + 1 2 ] T [ Δ K f K ¯ f Δ ϕ ϕ ¯ Δ ρ ρ ¯ ] ,
where the first term on the right of Equation (12) is the mapping operator between elastic parameters and reflection coefficients.
Nicolao et al. (1993) proposed that the eigenvalues of this mapping operator matrix can be used to measure the information of elastic parameters [16]. The eigenvalues are proportional to the information content. The direction cosine of the eigenvector of the mapping operator matrix can reflect the difficulty of inverting these elastic parameters. The greater the absolute value of the direction cosine is, the easier the inversion is. Moreover, when the direction cosine equals zero, the inversion is unviable. By using the first model in Table 3, the eigenvalues corresponding to the effective pore-fluid bulk modulus, porosity, and density are shown in Figure 3a. When the incident angle is smaller than 30°, the energy of the first eigenvalue (blue) is larger than that of the second eigenvalue (red) and the third eigenvalue (black). Consequently, we only need to analyze the eigenvector of the first eigenvalue to get a different degree in inverting these elastic parameters. As displayed in Figure 3b, the direction cosine of the effective pore-fluid bulk modulus, porosity, and density is far away from zero. This means that these three elastic parameters can be inverted by Equation (11).
The stability of inversion is very important. Therefore, the method proposed by Yin et al. (2018) is used to analyze the stability of Equation (11). Yin et al. (2018) deduced a three-term linear approximation including the effective pore-fluid bulk modulus [8], an approximation of which can be displayed as
R p p = [ ( 1 4 1 4 γ d r y 2 γ s a t 2 ) sec 2 θ ] Δ K f K ¯ f + [ ( γ d r y 2 4 γ s a t 2 2 γ s a t 2 sin 2 θ ) + 1 r ( sec 2 θ 4 γ d r y 2 2 γ s a t 2 sec 2 θ + 2 γ s a t 2 sin 2 θ ) ] Δ f m f ¯ m + ( 1 4 sec 2 θ + 1 2 ) Δ ρ ρ ¯ ,
where f m = ϕ μ , r is the fit coefficient of ϕ and μ according to field data.
The logarithm of the condition number of weighting coefficients can measure the stability of the inversion equation. The smaller the logarithm is, the more stable the inversion is. We compute the values of Equations (11) and (13), the results of which are shown in Figure 4. With the increase of incident angle, both of these two conditional numbers decrease. However, the conditional number of Equation (13) is larger than that of Equation (11). This indicates that the stability of inverting the effective pore-fluid bulk modulus, porosity, and density increases with an increase in the incident angle, and that the stability of Equation (11) is better than that of Equation (13).
The weighting coefficients before elastic parameters in Equation (11) include γ d r y 2 / γ s a t 2 . γ s a t is the ratio of P- to S-wave velocities in saturated rock and can be calculated by logging data, while γ d r y is the ratio of P- to S-wave velocities in dry rock. Russell et al. (2003, 2006, 2011) [5,17,18] estimated the value of γ d r y 2 through theoretical analysis and practical calculations, and proposed a possible range value for γ d r y 2 .
Here, we analyze the effect of γ d r y 2 / γ s a t 2 on the weighting coefficients of the effective pore-fluid bulk modulus, porosity, and density. According to Russell et al. (2011) [5], we assume that γ s a t 2 is 4 and γ d r y 2 is 2.33, 2.50, and 3.00, respectively. As shown in Figure 5, the weighting coefficient curves of these three elastic parameters are computed with an incident angle from 0° to 45°. With the variation of γ d r y 2 , the pore-fluid bulk modulus (Figure 5a) and porosity (Figure 5b) curves change, while the density curve (Figure 5c) is unaffected. Therefore, γ d r y 2 / γ s a t 2 influences only the effective pore-fluid bulk modulus and porosity, while the density is unchanged.
In this paper, we use the third-order approximation to establish a nonlinear AVO inversion method. This approximation includes the first-, second-, and third-order terms:
R p p R p p ( 1 ) + R p p ( 2 ) + R p p ( 3 ) .

2.3. Effective Pore-Fluid Bulk Modulus Nonlinear AVO Inversion Method

AVO inversion includes linear and nonlinear approaches. By adopting the third-order approximation, the AVO inversion method is strictly nonlinear. Because the efficiency of the nonlinear inversion method is pretty low in practice, series reversion is applied. Series reversion is an effective approach to transforming a nonlinear relationship into a linear relationship (Frank and Tat, 2001) [19]. For example, for the expression
y = x 1 x = i = 1 n x i ,
the solution is
x = y / ( 1 + y ) .
By using series reversion to solve Equation (15), x can be expanded in series as
y = i = 1 n x i ,
where x i is the ith-order item of y .
Substituting Equation (17) into (15), we can express x i as
x i = ( 1 ) i + 1 y i ,        i = 1 , 2 , , n .
By adopting Equation (18) into (17), x can be calculated as
x = x 1 + x 2 + x 3 + = y y 2 + y 3 + = y / ( 1 y ) .
In this way, the P-wave reflection coefficient, the effective pore-fluid bulk modulus, porosity, and density can be expanded in series. There are some tests for the effect of the expansion item number on accuracy and efficiency. When these parameters are expanded into three items, this method has advantages in both accuracy and efficiency. These parameters can be expressed as
R p p = R p p 1 + R p p 2 + R p p 3 ,
K f = K f 1 + K f 2 + K f 3 ,
ϕ = ϕ 1 + ϕ 2 + ϕ 3   ,
ρ = ρ 1 + ρ 2 + ρ 3   .
Combining Equation (20a–d) with (10a–d), we obtain
R p p 3 = R p p ( 1 ) ( K f 3 , ϕ 3 , ρ 3 ) + R p p ( 2 ) ( K f 2 , ϕ 2 , ρ 2 ) + R p p ( 3 ) ( K f 2 , ϕ 2 , ρ 2 ) + R p p ( 3 ) ( K f 1 , ϕ 1 , ρ 1 ) ,
R p p 2 = R p p ( 1 ) ( K f 2 , ϕ 2 , ρ 2 ) + R p p ( 2 ) ( K f 1 , ϕ 1 , ρ 1 ) ,
R p p 1 = R p p ( 1 ) ( K f 1 , ϕ 1 , ρ 1 ) .
Additionally, using the approach of Zhang and Weglein (2009) [20], we can calculate R p p 1 , R p p 2 , and R p p 3 through R p p in four steps. First, according to Equation (21c), K f 1 , ϕ 1 , and ρ 1 are inverted by the Bayesian theory. Second, K f 2 , ϕ 2 , and ρ 2 are inverted via the least square method based on Equation (21b). Third, K f 3 , ϕ 3 , and ρ 3 are inverted using the last step from Equation (21a). Finally, K f , ϕ , and ρ are computed using Equation (20a–d).
In the first three steps, a linear inversion AVO method is functional. The AVO inversion is based on the convolution model, which is written as
d = W R = W G m ,
where, d , W , m , and G are the seismic data, seismic wavelet, elastic parameters, and mapping operator between d and m , respectively. Because G is usually an irreversible matrix, its conjugate matrix is commonly adopted to solve Equation (22). In general, the objective function of AVO inversion takes the form of an L2 norm, which is expressed as
J min = 1 2 ( G m d ) 2 d t .
AVO inversion is an ill-conditioned and ill-posed problem. The single point estimation method is commonly used to obtain optimal solutions for the objective function. However, it is unable to evaluate solutions. Thus, this paper introduces the Bayesian inversion theory. The mathematical expression of the Bayesian theory can be expressed as
P ( m | d ) = P ( d | m ) P ( m ) / P ( d ) ,
where P ( m | d ) denotes the posterior probability function, P ( d | m ) is the likelihood function, P ( m ) represents the prior distribution function (on which the results of this theory mainly depend), and P ( d ) is marginal distribution, which is always a constant. It is assumed that the prior information complies with multivariate Cauchy distribution, which is expressed as
P ( m ) = P o m exp [ 2 i = 1 N ln ( 1 + m T Φ i m ) ] ,
P o m = 1 π 2 N | Ψ | N / 2 ,
Φ = ( D i ) T Ψ 1 D i ,
D n l i = { 1 , n = 1 , l = i 1 , n = 2 , l = i + N 1 , n = 3 , l = i + 2 N 0 , e l s e ,
where N is the number of samples, and Ψ is a 3 × 3 scale matrix.
The likelihood function obeys multivariate Gaussian distribution as
P ( d | m ) = 1 ( 2 π ) N 2 det C d exp { 1 2 ( d G m ) T C d 1 ( d G m ) } } ,
where C d denotes the noise covariance matrix.
The posterior probability distribution of parameters to be solved can be indicated as
P ( m | d ) = { 1 2 ( d - G m ) T C d 1 ( d - G m ) 2 i = 1 N ln ( 1 + m T Φ i m ) } .
In terms of Equations (23) and (27), the objective function becomes
J ( m ) = 1 2 ( d - G m ) T C d 1 ( d - G m ) + 2 i = 1 N ln ( 1 + m T Φ i m ) .
By taking the derivative of Equation (28) with respect to m , m is given.

3. Results

In this section, a logging model and field data are used to validate the advantages of our new inversion method.

3.1. Model Data Test

A logging model from sandstone reservoirs was used to verify the feasibility of our inversion method. The synthetic seismic record generated by the Ricker wavelet with a dominant frequency of 25 Hz is shown in Figure 6a (without noise) and Figure 6b (signal-to-noise ratio equals 2). Figure 7a displays the inverted results of the synthetic gather without noise. The logging data, elastic parameters inverted by the previous linear method (Yin et al., 2018), and proposed nonlinear method in this paper are shown using black, green, and red lines. The errors between the logging data and the inverted results are displayed in Figure 7b. The inversion results of the effective pore-fluid bulk modulus, porosity, and density had the same trend as the logging data, while the curves of the proposed nonlinear method were closer to the logging data curves. Figure 7c displays the values of L2 norm and shows that these two methods could obtain a satisfactory result within five iterations. The final value of the L2 norm in the proposed method case was smaller than that in the method of Yin et al. (2018). Furthermore, we added noise to the synthetic angle gathers, as shown in Figure 6b. The inverted elastic parameters of the synthetic angle gathers with and without noise are displayed in Figure 8a. Figure 8b shows the errors of logging data and the inverted results. The error difference was relatively larger in the noise case than in the noise-free case. However, the error between inverted elastic parameters and logging data was not so severe. As can be seen in Figure 8c, the curves of the L2 norm were close. Using the synthetic gathers without noise, we compared the accuracy and efficiency of the proposed nonlinear method and the traditional fully nonlinear method (Zhou et al., 2016) [21]. Figure 9a shows the inverted results of these two methods. The inverted values were both close to those of the logging data. The curves of error are displayed in Figure 9b, without much difference between the curves. However, the traditional fully nonlinear method took about 15 times as long as the proposed method.

3.2. Field Data Test

Additionally, a real 2D dataset was used to measure the feasibility of the proposed method, with mainly sandstone reservoirs as the target layers. Figure 10 is the post-stack seismic section, with gas layers at 1.601 and 1.662 s, and water layers at 1.579, 1.620, 1.628, and 1.642 s, according to the logging interpretation. Based on the logging data, we calculated the effective pore-fluid bulk modulus to test its sensitivity to pore-fluid. The cross-plot of the pore-fluid bulk modulus and porosity is shown in Figure 11. The different pore-fluid types can be clearly identified by the pore-fluid bulk modulus. We transformed the pre-stack gathers from the offset domain to the angle domain and extracted angle gathers from 3° to 38°. A wavelet was statistically determined from the field seismic data for inversion. Finally, the proposed inversion method was applied to obtain the elastic parameters.
The inverted effective pore-fluid bulk modulus and density are displayed in Figure 12. In Figure 12a, the values of purple area range from −1.5 to −0.5 GPa, which means gas sand, while the red color ranges from −4 to −3 GPa, representing water sand. There are two purple areas matching the gas sandstone reservoirs in the profile of the effective pore-fluid bulk modulus. Four red layers exist and all correspond to the positions of water sandstone reservoirs. As shown in Figure 12b, the gas sand in purple has an approximate value of 2.30 g/cm3. The value of the water sand in blue is around 2.38 g/cm3. For the density profile, there are two purple layers and one blue layer. However, the first purple layer at 1.579 s is not gas sand but water sand. The bottom three layers cannot be identified on the profile. Note that the effective pore-fluid bulk modulus showed a better correlation with fluid interpretation than density did.

4. Discussion

In this study, we addressed three topics. First, we favored the effective pore-fluid bulk modulus as a fluid indicator to identify pore-fluid in a sandstone reservoir. A data analysis model followed and verifies the high sensitivity to pore-fluid change by the effective pore-fluid bulk modulus (Table 2 and Figure 1). With this in mind, a new high-order AVO approximation was deduced based on the Zoeppritz equations and Biot–Gassmann theory. Appendix A includes the detailed formula derivations. Figure 2, Figure 3 and Figure 4 show the high accuracy and stability of this new approximation. Second, we introduce the series reversion and Bayesian theory to build a nonlinear AVO inversion method. The results of logging data tests showed that the proposed nonlinear inversion method had higher accuracy than the linear method (Figure 7). The stability and noise-insensitivity of the proposed method were high because of the Bayesian framework (Figure 8). The proposed method showed high accuracy and efficiency by comparison with the traditional fully nonlinear inversion method (Figure 9). The field data test showed the practicability of the proposed method and the strong identification ability of the effective pore-fluid bulk modulus (Figure 11). Third, we briefly illustrated the limitation and expectations of this work. Compared with other fluid indicators, the effective pore-fluid bulk modulus was more sensitive to pore-fluid change. However, this methodology was only applicable for sandstone reservoirs. Moreover, fluid indicators should be selected according to the actual reservoir conditions. An assumption that the lithology of the target reservoir was constant was made during the derivation of the new approximation. This presumption is inconsistent with real conditions, and will produce errors when applying the new approximation into field data. A new study direction would be to address this assumption. Finally, we use an AVO inversion method to obtain the effective pore-fluid bulk modulus from seismic data. This method relies heavily on the initial model, which is difficult to build satisfactorily for field data inversion. Deep learning would be a promising alternative for the future inverted method.

5. Conclusions

An effective pore-fluid bulk modulus is regarded as a fluid indicator that can identify pore-fluid in sandstone. Based on the Zoeppritz equations and Biot–Gassmann theory, we first derived a high-order approximation that estimated the effective pore-fluid bulk modulus directly. By introducing the series reversion and Bayesian theory, a nonlinear AVO inversion method was proposed to invert the effective pore-fluid bulk modulus. Model results demonstrate that the high-order AVO approximation had higher precision than linear approximations and could be applied to AVO inversion. The results of the field data test show the stability, effectiveness, and feasibility of our method in identifying pore-fluid in a sandstone reservoir.

Author Contributions

Conceptualization, L.S., Y.S., D.C., and Y.L.; Data curation, L.S., D.C., and Y.S.; Funding acquisition, L.S., Y.L., and J.L.; Investigation, L.S. and Y.L.; Methodology, Y.S., D.C., and Y.L.; Resources, L.S. and Y.L.; Supervision, J.L.; Writing—original draft, L.S. and Y.S.; Writing—review & editing, Y.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Important National Science & Technology Specific Project of China under contract 2016ZX05002005-003.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

The contents of this appendix show the weighting coefficients before Δ K f / K ¯ f , Δ ϕ / ϕ ¯ , and Δ ρ / ρ ¯ of R P P ( 1 ) , R P P ( 2 ) and R P P ( 3 ) .

Appendix A.1. The Weighting Coefficients before Δ K f / K ¯ f , Δ ϕ / ϕ ¯ , and Δ ρ / ρ ¯ of R P P ( 1 )

R p p ( 1 ) = W 11 Δ K f K ¯ f + W 12 Δ ϕ ϕ ¯ + W 13 Δ Δ ρ ρ ¯   ,
where
W 11 = 1 4 ( 1 γ d r y 2 γ s a t 2 ) sec 2 ( θ ) ,
W 12 = 1 4 ( 1 γ d r y 2 γ s a t 2 ) sec 2 ( θ ) + [ 1 4 γ d r y 2 γ s a t 2 sec 2 ( θ ) 2 γ s a t 2 sin 2 ( θ ) ] γ ϕ ,
W 13 = 1 2 1 4 sec 2 ( θ ) .

Appendix A.2. The Weighting Coefficients before Δ K f / K ¯ f , Δ ϕ / ϕ ¯ , and Δ ρ / ρ ¯ of R P P ( 2 )

R p p ( 2 ) = W 21 ( Δ K f K ¯ f ) 2 + W 22 Δ ϕ ϕ ¯ 2 + W 23 Δ ρ ρ ¯ 2 + W 24 Δ K f K ¯ f Δ ϕ ϕ ¯ + W 25 Δ K f K ¯ f Δ ρ ρ ¯ + W 26 Δ ϕ ϕ ¯ Δ ρ ρ ¯ ,
where
W 21 = 1 8 ( 1 γ d r y 2 γ s a t 2 ) [ sin 2 ( θ ) + γ d r y 2 γ s a t 2 ] ,
W 22 = 1 8 ( 1 γ d r y 2 γ s a t 2 ) [ sin 2 ( θ ) + γ d r y 2 γ s a t 2 ] + { 1 8 γ d r y 2 γ s a t 2 [ sec 2 ( θ ) γ d r y 2 γ s a t 2 ] 1 γ s a t 2 ( 1 1 γ s a t ) sin 2 ( θ ) } γ ϕ 2 1 4 γ d r y 2 γ s a t 2 ( 1 γ d r y 2 γ s a t 2 ) γ ϕ ,
W 23 = 1 8 ( 1 2 γ s a t ) sin 2 ( θ ) ,
W 24 = 1 4 ( 1 γ d r y 2 γ s a t 2 ) [ sin 2 ( θ ) + γ d r y 2 γ s a t 2 ] 1 4 γ d r y 2 γ s a t 2 ( 1 γ d r y 2 γ s a t 2 ) γ ϕ ,
W 25 = 1 4 ( 1 γ d r y 2 γ s a t 2 ) sin 2 ( θ ) ,
W 26 = 1 4 ( 1 γ d r y 2 γ s a t 2 ) sin 2 ( θ ) + 1 γ s a t 2 ( 1 γ d r y 2 4 ) sin 2 ( θ ) γ ϕ .

Appendix A.3. The Weighting Coefficients before Δ K f / K ¯ f , Δ ϕ / ϕ ¯ , and Δ ρ / ρ ¯ of R P P ( 3 )

R p p ( 3 ) = W 30 Δ K f K ¯ f 3 + W 31 Δ ϕ ϕ ¯ 3 + W 32 Δ ρ ρ ¯ 3 + W 33 Δ K f K ¯ f 2 Δ ϕ ϕ ¯ + W 34 Δ K f K ¯ f 2 Δ ρ ρ ¯ + W 35 Δ ϕ ϕ ¯ 2 Δ K f K ¯ f + W 36 Δ ϕ ϕ ¯ 2 Δ ρ ρ ¯ + W 37 Δ ρ ρ ¯ 2 Δ K f K ¯ f + W 38 Δ ρ ρ ¯ 2 Δ ϕ ϕ ¯ + W 39 Δ K f K ¯ f Δ ϕ ϕ ¯ Δ ρ ρ ¯ ,
where
W 30 = 3 64 ( 1 γ d r y 4 γ s a t 4 ) sin 2 ( θ ) 3 64 γ d r y 2 γ s a t 2 ( 1 + 1 3 sin 2 ( θ ) ) + 1 64 + 7 64 γ d r y 4 γ s a t 4 5 64 γ d r y 6 γ s a t 6 ( 1 1 5 sin 2 ( θ ) ) ,
W 31 = [ 3 64 ( 1 γ d r y 4 γ s a t 4 ) sin 2 ( θ ) 3 64 γ d r y 2 γ s a t 2 ( 1 + 1 3 sin 2 ( θ ) ) 1 64 + 7 64 γ d r y 4 γ s a t 4 5 64 γ d r y 6 γ s a t 6 ( 1 1 5 sin 2 ( θ ) ) ] + [ 5 64 γ d r y 6 γ s a t 6 ( 1 1 5 sin 2 ( θ ) ) 1 2 1 γ s a t 2 ( 1 + γ d r y 2 γ s a t 2 ) sin 2 ( θ ) + 1 16 γ d r y 2 γ s a t 2 sec 2 ( θ ) 1 8 γ d r y 4 γ s a t 4 ( 1 1 γ s a t 2 sin 2 ( θ ) ) + 3 4 1 γ s a t 3 sin 2 ( θ ) ] γ ϕ 3 + [ 1 8 1 γ s a t 2 sin 2 ( θ ) 3 64 γ d r y 2 γ s a t 2 ( 1 + γ d r y 4 γ s a t 4 ) sin 2 ( θ ) 1 4 γ d r y 2 γ s a t 2 ( 1 1 2 γ d r y 2 γ s a t 2 ) sin 2 ( θ ) 11 32 γ d r y 4 γ s a t 4 ( 1 3 11 sin 2 ( θ ) ) + 7 64 γ d r y 2 γ s a t 2 ( 1 + 15 7 γ d r y 4 γ s a t 4 ) ] γ + [ ( 1 4 γ d r y 2 γ s a t 2 1 2 1 γ s a t 3 ) ( 1 γ d r y 2 γ s a t 2 ) sin 2 ( θ ) 15 64 γ d r y 6 γ s a t 6 ( 1 1 5 sin 2 ( θ ) ) + 23 64 γ d r y 4 γ s a t 4 ( 1 3 23 sin 2 ( θ ) ) 1 8 γ d r y 2 γ s a t 2 ] γ ϕ 2 ,
W 32 = 1 64 3 64 sin 2 ( θ ) + 1 16 1 γ s a t sin 2 ( θ ) ,
W 33 = [ 9 64 ( 1 γ d r y 4 γ s a t 4 ) sin 2 ( θ ) 9 64 γ d r y 2 γ s a t 2 ( 1 + 1 3 sin 2 ( θ ) ) + 3 64 + 21 64 γ d r y 4 γ s a t 4 15 64 γ d r y 6 γ s a t 6 ( 1 1 5 sin 2 ( θ ) ) ] + [ 1 8 1 γ s a t 2 sin 2 ( θ ) 3 64 γ d r y 2 γ s a t 2 ( 1 + γ d r y 4 γ s a t 4 ) sin 2 ( θ ) 1 4 γ d r y 2 γ s a t 2 ( 1 1 2 γ d r y 2 γ s a t 2 ) sin 2 ( θ ) 11 32 γ d r y 4 γ s a t 4 ( 1 3 11 sin 2 ( θ ) ) + 7 64 γ d r y 2 γ s a t 2 ( 1 + 15 7 γ d r y 4 γ s a t 4 ) ] γ ϕ ,
W 34 = 1 32 γ d r y 2 γ s a t 2 ( 1 + 5 sin 2 ( θ ) ) 1 64 γ d r y 4 γ s a t 4 sec 2 ( θ ) 1 64 9 64 sin 2 ( θ ) ,
W 35 = [ 9 64 ( 1 γ d r y 4 γ s a t 4 ) sin 2 ( θ ) 9 64 γ d r y 2 γ s a t 2 ( 1 + 1 3 sin 2 ( θ ) ) + 3 64 + 21 64 γ d r y 4 γ s a t 4 15 64 γ d r y 6 γ s a t 6 ( 1 1 5 sin 2 ( θ ) ) ] + 2 [ 1 8 1 γ s a t 2 sin 2 ( θ ) 3 64 γ d r y 2 γ s a t 2 ( 1 + γ d r y 4 γ s a t 4 ) sin 2 ( θ ) 1 4 γ d r y 2 γ s a t 2 ( 1 1 2 γ d r y 2 γ s a t 2 ) sin 2 ( θ ) 11 32 γ d r y 4 γ s a t 4 ( 1 3 11 sin 2 ( θ ) ) + 7 64 γ d r y 2 γ s a t 2 ( 1 + 15 7 γ d r y 4 γ s a t 4 ) ] γ ϕ + [ ( 1 4 γ d r y 2 γ s a t 2 1 2 1 γ s a t 3 ) ( 1 γ d r y 2 γ s a t 2 ) sin 2 ( θ ) 15 64 γ d r y 6 γ s a t 6 ( 1 1 5 sin 2 ( θ ) ) + 23 64 γ d r y 4 γ s a t 4 ( 1 3 23 sin 2 ( θ ) ) 1 8 γ d r y 2 γ s a t 2 ] γ ϕ 2 ,
W 36 = 1 32 γ d r y 2 γ s a t 2 ( 1 + 5 sin 2 ( θ ) ) 1 64 γ d r y 4 γ s a t 4 sec 2 ( θ ) 1 64 9 64 sin 2 ( θ ) + [ 3 4 1 γ s a t 2 ( 1 1 γ s a t ) sin 2 ( θ ) 1 8 γ d r y 2 γ s a t 2 ( 1 2 γ s a t 2 ) sin 2 ( θ ) 1 64 γ d r y 4 γ s a t 4 sec 2 ( θ ) ] γ ϕ 2 + ( 1 γ d r y 2 γ s a t 2 ) ( 1 4 1 γ s a t 2 sin 2 ( θ ) 1 32 γ d r y 2 γ s a t 2 sec 2 ( θ ) ) γ ϕ ,
W 37 = ( 1 γ d r y 2 γ s a t 2 ) ( 9 64 sin 2 ( θ ) 1 8 1 γ s a t sin 2 ( θ ) 1 64 ) ,
W 38 = ( 1 γ d r y 2 γ s a t 2 ) ( 9 64 sin 2 ( θ ) 1 8 1 γ s a t sin 2 ( θ ) 1 64 ) + [ 1 64 γ d r y 2 γ s a t 2 ( 1 9 sin 2 ( θ ) ) 1 16 1 γ s a t ( 1 + 2 γ d r y 2 γ s a t 2 ) sin 2 ( θ ) 1 8 1 γ s a t 2 sin 2 ( θ ) ] γ ϕ ,
W 39 = 1 16 γ d r y 2 γ s a t 2 ( 1 + 5 sin 2 ( θ ) ) 1 32 γ d r y 4 γ s a t 2 sec 2 ( θ ) 1 32 9 32 sin 2 ( θ ) + ( 1 γ d r y 2 γ s a t 2 ) ( 1 4 1 γ s a t 2 sin 2 ( θ ) 1 32 γ d r y 2 γ s a t 2 sec 2 ( θ ) ) γ ϕ .

References

  1. Smith, G.C.; Gidlow, P.M. Weighting stacking for rock property estimation and detection of gas. Geophys. Prospect. 1987, 35, 993–1014. [Google Scholar] [CrossRef]
  2. Fatti, J.L.; Smith, G.C.; Vail, P.J.; Strauss, P.J.; Levitt, P.R. Detection of gas in sandstone reservoirs using AVO analysis: A 3-D seismic case history using the Geostack technique. Geophysics 1994, 59, 1362–1376. [Google Scholar] [CrossRef]
  3. Goodway, B.; Chen, T.; Downton, J. Improved AVO Fluid Detection and Lithology Discrimination Using Lame Petrophysical Parameters: “λρ”, “μρ” and “λ/μ Fluid Stack” from P and S Inversions; SEG Expanded Abstracts: Dallas, TX, USA, 1997; pp. 183–186. [Google Scholar]
  4. Quakenbush, M.; Shang, B.; Tuttle, C. Poisson impedance. Lead. Edge 2006, 25, 128–138. [Google Scholar] [CrossRef]
  5. Russell, B.H.; Gray, D.; Hampson, D.P. Linearized AVO and poroelasticity. Geophysics 2011, 76, C19–C29. [Google Scholar] [CrossRef]
  6. Zhang, S.; Yin, X.; Zhang, F. Quasi Fluid Modulus for Delicate Lithology and Fluid Discrimination; SEG Expanded Abstracts: Denver, CO, USA, 2010; pp. 404–408. [Google Scholar]
  7. Yin, X.; Zhang, S. Bayesian inversion for effective pore-fluid bulk modulus based on fluid-matrix decoupled amplitude variation with offset approximation. Geophysics 2014, 79, R221–R232. [Google Scholar] [CrossRef]
  8. Yin, X.; Wang, H.; Cao, D.; Zhou, Q.; Sun, W. Three-term AVO approximation of Kf-fm and prestack seismic inversion for deep reservoirs. Oil Geophys. Prospect. 2018, 53, 129–135. [Google Scholar]
  9. Hilterman, F.J. Seismic Amplitude Interpretation; Society of Exploration Geophysicists: Tulsa, OK, USA, 2001; Volume 4. [Google Scholar]
  10. Zoeppritz, K. Erdbebenwellen VIII B: Über die Reflexion and Durchgang seismischer Wellen durch Unstetigkeitsflächen. Gott. Nachr. 1919, 1, 66–84. [Google Scholar]
  11. Keys, R.G. Polarity reversals in reflections from layered media. Geophysics 1989, 54, 900–905. [Google Scholar] [CrossRef]
  12. Kim, S.; Innanen, K. Linear and nonlinear poroelastic AVO. CREWES Res. Rep. 2013, 25, 1–28. [Google Scholar]
  13. Han, D.; Batzle, M. Gain Function and Hydrocarbon Indicators; SEG Expanded Abstracts: Dallas, TX, USA, 2003; pp. 1695–1698. [Google Scholar]
  14. Nur, A.; Mavko, G.; Dvorkin, J.; Galmudi, D. Critical porosity: A key to relating physical properties to porosity in rocks. Lead. Edge 1998, 17, 357–362. [Google Scholar] [CrossRef]
  15. Aki, K.; Richards, P.G. Quantitative Seismology; W. H. Freeman and Company: San Francisco, CA, USA, 1980. [Google Scholar]
  16. De Nicolao, A.; Drufuca, G.; Rocca, F. Eigenvalues and eigenvectors of linearized elastic inversion. Geophysics 1993, 58, 670–679. [Google Scholar] [CrossRef]
  17. Russell, B.H.; Hedlin, K.; Hilterman, F.J. Fluid-property discrimination with AVO: A Biot-Gassmann perspective. Geophysics 2003, 68, 29–39. [Google Scholar] [CrossRef] [Green Version]
  18. Russell, B.H.; Gray, D.; Hampson, D.P.; Lines, L.R. Linearized AVO and poroelasticity. CREWES Res. Rep. 2006, 18, 1–25. [Google Scholar] [CrossRef]
  19. Frank, H.W.; Tat, S.Y. New applications of non-linear chirp scaling in SAR data processing. IEEE Trans. Geosci. Remote Sens. 2001, 39, 946–953. [Google Scholar]
  20. Zhang, H.; Weglein, A.B. Direct nonlinear inversion of multiparameter 1D elastic media using the inverse scattering series. Geophysics 2009, 74, WCD15–WCD27. [Google Scholar] [CrossRef]
  21. Zhou, L.; Li, J.Y.; Chen, X.H. Nonlinear three-term AVO inversion based on exact Zoeppritz equation. Chin. J. Geophys. 2016, 59, 2663–2673. [Google Scholar]
Figure 1. The relative variations of fluid indicators in sandstone models: (a) Type 1; (b) Type 2; (c) Type 3. M, F, and Kf represent the P-wave modulus, Gassmann fluid factor, and the effective pore-fluid bulk modulus.
Figure 1. The relative variations of fluid indicators in sandstone models: (a) Type 1; (b) Type 2; (c) Type 3. M, F, and Kf represent the P-wave modulus, Gassmann fluid factor, and the effective pore-fluid bulk modulus.
Energies 13 01313 g001
Figure 2. P-wave reflection coefficients computed by different equations: (a) Model 1; (b) Model 2; (c) Model 3.
Figure 2. P-wave reflection coefficients computed by different equations: (a) Model 1; (b) Model 2; (c) Model 3.
Energies 13 01313 g002
Figure 3. (a) Comparisons of three eigenvalues vs the maximum incident angle in the data domain; (b) Comparisons of directional cosines of the first eigenvector vs the maximum incident angle in the model domain.
Figure 3. (a) Comparisons of three eigenvalues vs the maximum incident angle in the data domain; (b) Comparisons of directional cosines of the first eigenvector vs the maximum incident angle in the model domain.
Energies 13 01313 g003
Figure 4. Analysis of the condition number of the coefficient matrix in Equations (11) and (13).
Figure 4. Analysis of the condition number of the coefficient matrix in Equations (11) and (13).
Energies 13 01313 g004
Figure 5. Weighting coefficients of γ d r y 2 = 2.33 , γ d r y 2 = 2.50 , and γ d r y 2 = 3.00 . (a) The effective pore-fluid bulk modulus; (b) porosity; (c) density.
Figure 5. Weighting coefficients of γ d r y 2 = 2.33 , γ d r y 2 = 2.50 , and γ d r y 2 = 3.00 . (a) The effective pore-fluid bulk modulus; (b) porosity; (c) density.
Energies 13 01313 g005
Figure 6. P-wave synthetic angle gathers: (a) without noise; (b) SNR = 2.
Figure 6. P-wave synthetic angle gathers: (a) without noise; (b) SNR = 2.
Energies 13 01313 g006
Figure 7. (a) Inversion results of the effective pore-fluid bulk modulus, porosity, and density by the linear and proposed methods; (b) Error; (c) The values of the L2 norm.
Figure 7. (a) Inversion results of the effective pore-fluid bulk modulus, porosity, and density by the linear and proposed methods; (b) Error; (c) The values of the L2 norm.
Energies 13 01313 g007
Figure 8. (a) Inversion results of the effective pore-fluid bulk modulus, porosity, and density based on the synthetic gathers with and without noise; (b) Error; (c) The values of the L2 norm.
Figure 8. (a) Inversion results of the effective pore-fluid bulk modulus, porosity, and density based on the synthetic gathers with and without noise; (b) Error; (c) The values of the L2 norm.
Energies 13 01313 g008
Figure 9. (a) Inversion results of the effective pore-fluid bulk modulus, porosity, and density by the fully nonlinear and proposed methods; (b) Error.
Figure 9. (a) Inversion results of the effective pore-fluid bulk modulus, porosity, and density by the fully nonlinear and proposed methods; (b) Error.
Energies 13 01313 g009
Figure 10. Field data seismic section.
Figure 10. Field data seismic section.
Energies 13 01313 g010
Figure 11. Cross-plot of the effective pore-fluid bulk modulus and porosity.
Figure 11. Cross-plot of the effective pore-fluid bulk modulus and porosity.
Energies 13 01313 g011
Figure 12. Inverted parameters sections: (a) the effective pore-fluid bulk modulus section; (b) the density section. Purple, black, and green indicate gas sand, water sand, and the surrounding rocks, respectively.
Figure 12. Inverted parameters sections: (a) the effective pore-fluid bulk modulus section; (b) the density section. Purple, black, and green indicate gas sand, water sand, and the surrounding rocks, respectively.
Energies 13 01313 g012aEnergies 13 01313 g012b
Table 1. Three types of sandstone models by Hilterman (2001) [9].
Table 1. Three types of sandstone models by Hilterman (2001) [9].
TypePore-FluidP-Wave Velocity (m/s)S-Wave Velocity (m/s)Density (g/cm3)Porosity (%)
1Water411524532.3215
Gas405025262.2120
2Water304815952.2315
Gas278116652.0820
3Water21348602.1115
Gas15439011.8820
Table 2. The calculation results of fluid indicators in sandstone models.
Table 2. The calculation results of fluid indicators in sandstone models.
TypePore-FluidP-wave Modulus (GPa)Gassmann Fluid Factor (GPa)Effective Pore-Fluid Bulk Modulus (GPa)
1Water39.29−16.56−17.66
Gas36.25−20.16−12.90
2Water20.72−19.75−21.07
Gas16.09−6.98−5.58
3Water9.613.373.59
Gas4.48−1.63−1.70
Table 3. Parameters of the sand model (Yin and Zhang, 2014) [7].
Table 3. Parameters of the sand model (Yin and Zhang, 2014) [7].
Pore-Fluid K f   ( GPa ) ρ (g/cm3) V p   ( m / s ) V s   ( m / s ) ϕ (%) ϕ c   ( % )
Model 1Gas0.101.99192012302540
Water2.382.26259011502540
Model 2Water2.382.26259011502540
Water2.382.34358020702540
Model 3Gas0.102.12334021702540
Water2.382.26259011502540

Share and Cite

MDPI and ACS Style

Shi, L.; Sun, Y.; Liu, Y.; Cova, D.; Liu, J. High-Order AVO Inversion for Effective Pore-Fluid Bulk Modulus Based on Series Reversion and Bayesian Theory. Energies 2020, 13, 1313. https://doi.org/10.3390/en13061313

AMA Style

Shi L, Sun Y, Liu Y, Cova D, Liu J. High-Order AVO Inversion for Effective Pore-Fluid Bulk Modulus Based on Series Reversion and Bayesian Theory. Energies. 2020; 13(6):1313. https://doi.org/10.3390/en13061313

Chicago/Turabian Style

Shi, Lei, Yuhang Sun, Yang Liu, David Cova, and Junzhou Liu. 2020. "High-Order AVO Inversion for Effective Pore-Fluid Bulk Modulus Based on Series Reversion and Bayesian Theory" Energies 13, no. 6: 1313. https://doi.org/10.3390/en13061313

APA Style

Shi, L., Sun, Y., Liu, Y., Cova, D., & Liu, J. (2020). High-Order AVO Inversion for Effective Pore-Fluid Bulk Modulus Based on Series Reversion and Bayesian Theory. Energies, 13(6), 1313. https://doi.org/10.3390/en13061313

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop