[0017] The shear wave velocity estimation method based on the inversion of the equivalent elastic modulus of adaptive matrix minerals, the specific implementation steps are as follows:
[0018] (1) Taking a single logging depth point as the research object, according to the logging data (P-wave velocity, density and porosity) at that depth point, the empirical relationship between the longitudinal wave velocity and the shear wave velocity of sand and mudstone is used (Castagna, 1993) to estimate Initial shear wave velocity, and calculate the range, initial value and increment of the Poisson’s ratio of the dry rock and the equivalent bulk modulus of the matrix mineral at the logging depth point.
[0019] The empirical relationship between longitudinal wave velocity and shear wave velocity of sand and mudstone (Castagna, 1993) is shown in formula 1.
[0020] V s = 0.804V p -0.856 (1)
[0021] Where V p Represents longitudinal wave velocity, V s Represents transverse wave velocity.
[0022] Poisson's ratio of dry rock (σ dry The value range of) is set to 0.0~0.4 (the value range of Poisson’s ratio of common sedimentary rocks). The initial value is generally equal to the lower or upper limit of the value range, and the increment is generally between 0.01 and 0.05. between.
[0023] The calculation formula for the lower limit and upper limit of the equivalent bulk modulus range of the matrix mineral is shown in formula 2.
[0024] lower limit: ρ ( V p 2 - 4 3 V s 2 ) / ( 1 - φ ) [ 3 / ( 1 - φ ) ] , Upper limit: [ ρ ( V p 2 - 4 3 V s 2 ) + aφ ] / ( 1 - φ ) [ 3 / ( 1 - φ ) ] - - - ( 2 ) Where a φ Is a correction item related to porosity, a is a constant (generally between 100 and 400), ρ and φ are the density and porosity obtained by logging. The initial value is generally equal to the lower or upper limit of the value range. The increment is generally between 0.0 and 1.0 GPa.
[0025] (2) Two different fluid term calculation methods (based on the Gassmann equation and Gassmann-Boit-Geertsman equation and the method based on the fluid identification factor proposed by Russell et al. in 2003) are used to calculate the fluid term at the same time, and set the above The absolute value of the difference between the two fluid terms is the inversion objective function, and the method of finding the global optimal solution is used to invert the optimal Poisson’s ratio of dry rock and the equivalent bulk modulus of the matrix mineral, when the value of the objective function < Set the threshold (selected according to the calculation accuracy and calculation efficiency, the value is generally between 0 and 0.5), the calculation is stopped, and the equivalent bulk modulus of the rock matrix mineral and the Poisson of dry rock at the current depth point are output ratio.
[0026] Calculation method of fluid term based on Gassmann equation and Gassmann-Boit-Geertsman equation:
[0027] The fluid term referred to in the present invention is the fluid term constituting the longitudinal wave velocity, that is, f in formula 3. Formula 4 is the calculation formula of the fluid term based on Gassmann equation, and Formula 5 is Gassmann-Boit-Geertsman equation, which is a standard one-dimensional quadratic equation with β as a variable. If ρ, V are known p , Σ dry , K 0 , K fl , Φ, β can be solved; then, according to formula 4, the fluid term (f) can be calculated.
[0028] V p = S + f ρ , S = K dry + 4 3 u dry = λ dry + 2 u dry - - - ( 3 )
[0029] f = β 2 M = ( 1 - K dry K 0 ) 2 φ K fl + 1 - φ K 0 - K dry K 0 2 , β = 1 - K dry K 0 , 1 M = β - φ K 0 + φ K fl - - - ( 4 )
[0030] ( Y - 1 ) β 2 + [ Yφ ( K 0 K fl - 1 ) - Y + N K 0 ] β - φ ( Y - N K 0 ) ( K 0 K fl - 1 ) = 0 , Y = 3 ( 1 - σ dry ) ( 1 + σ dry ) , N = ρV P 2 - - - ( 5 )
[0031] In formula 3, formula 4 and formula 5, S represents the framework term, f represents the fluid term, ρ represents the density of fluid-saturated rock (take the logging density value), and K dy , U dry , Λ dry And σ dry Respectively represent the bulk modulus, shear modulus, Lame constant and Poisson’s ratio of dry rock, K 0 Represents the equivalent bulk modulus of the matrix mineral, K fl Represents the bulk modulus of the fluid, and φ represents the porosity.
[0032] The calculation method of fluid term based on the fluid identification factor derived by Russell et al. in 2003:
[0033] ρf = ( Z p 2 - c Z s 2 ) , c = ( V p V s ) dry 2 = 2 ( 1 - σ dry ) 1 - 2 σ dry , Z p =ρV p ,Z s =ρV s (6)
[0034] Ρf in formula 6 is the fluid identification factor, ρ, Z p And Z s They represent the density, P-wave impedance and S-wave impedance of fluid-saturated rocks. If V is known p , V s , Ρ and σ dry , The fluid term (f) can be calculated.
[0035] (3) Based on the Poisson's ratio of the dry rock and the equivalent bulk modulus of the matrix mineral obtained in step 2 inversion, the linear average formula (Equation 7) and Gassmann equation (Equation 8) based on the critical porosity model of Nur (1995) )Calculate the equivalent shear modulus (μ 0 ).
[0036] μ 0 = μ dry ( 1 - φ φ c ) , μ dry = 4 3 K dry ( S - 1 ) , S = 3 ( 1 - σ dry ) 1 + σ dry - - - ( 7 )
[0037] K dry = K sat ( φK 0 K fl + 1 - φ ) - K 0 φK 0 K fl + K sat K 0 - 1 - φ , K sat = ρ ( V p 2 - 4 3 V s 2 ) - - - ( 8 )
[0038] In Equation 7, μ dry And σ dry Represents the shear modulus and Poisson's ratio of dry rock, φ c Represents critical porosity (generally taken as 0.40).
[0039] (4) Based on the equivalent bulk modulus and equivalent shear modulus of the matrix minerals obtained in steps 2 and 3, the simplified Xu-White model (Keys and Xu, 2002) is used to estimate the shear wave velocity. The calculation formula is shown in Equation 9 And formula 10.
[0040] V s = μ d / ρ , μ d =μ 0 (1-φ) q , q = 1 5 X l = s , sh υ l F iijj ( α l ) - - - ( 9 )
[0041] F iijj ( α ) = 2 F 3 + 1 F 4 + F 4 F 5 + F 6 F 7 - F 8 F 9 F 2 F 4 - - - ( 10 - 1 )
[0042] F 2 = 1 + A [ 1 + 3 2 ( g + γ ) - R 2 ( 3 g + 5 γ ) ] + B ( 3 - 4 R ) + A 2 ( A + 3 B ) ( 3 - 4 R ) [ g + γ - R ( g - γ + 2 υ 2 ) ] - - - ( 10 - 2 )
[0043] F 3 = 1 + A 2 [ R ( 2 - γ ) - 1 + α 2 α 2 g ( R - 1 ) ] - - - ( 10 - 3 )
[0044] F 4 = 1 + A 4 [ 3 γ + g - R ( g - γ ) ] - - - ( 10 - 4 )
[0045] F 5 = A [ R ( g + γ - 4 3 ) - g ] + Bγ ( 3 - 4 R ) - - - ( 10 - 5 )
[0046] P 6 =1+A[1+g-R(γ+g)]+B(1-γ)(3-4R) (10-6)
[0047] F 7 = 2 + A 4 [ 9 γ + 3 g - R ( 5 γ + 3 g ) ] + Bγ ( 3 - 4 R ) - - - ( 10 - 7 )
[0048] F 8 = A [ 1 - 2 R + g 2 ( R - 1 ) + υ 2 ( 5 R - 3 ) ] + B ( 1 - γ ) ( 3 - 4 R ) - - - ( 10 - 8 )
[0049] F 9 =A[g(R-1)-Rγ]+Bγ(3-4R) (10-9)
[0050] A=-1, B=0, R = 3 μ 0 3 K 0 + 4 μ 0 , g = α 2 1 - α 2 ( 3 γ - 2 ) , γ = α ( 1 - α 2 ) 3 / 2 [ cos - 1 ( α ) - α 1 - α 2 ] - - - ( 10 - 10 )
[0051] In formula 9 and formula 10, υ l And α l Respectively represent the volume percentage of sandstone and mudstone in the rock matrix (calculated based on the logging mud content) and pore flatness (the flatness of sandstone pores is 0.10, and the flattening of mudstone pores is 0.04, which is similar to the conventional Xu-White model The values are the same), F iiij (α l ) Represents the pore flatness α l Eshelby tensor at time.
[0052] (5) Select the next logging depth point and repeat the calculation in steps 1 to 5; stop the calculation until all the logging depth points in the target depth section are calculated.
[0053] (6) Output the shear wave velocity value estimated at each logging depth point in the target depth section.
[0054] figure 1 with figure 2 It is a calculation example of shear wave velocity estimation based on actual logging data. based on figure 1 The five input parameters shown, including longitudinal wave velocity, density, porosity, water saturation, and shale content, are calculated in strict accordance with the 6 steps of the specific embodiment of the present invention to obtain figure 2 The estimated shear wave velocity curve in, then use the measured shear wave velocity and the estimated shear wave velocity to calculate the absolute value of the error (|estimated shear wave velocity-measured shear wave velocity|) and error percentage The values of the parameters involved in the above calculation are as follows:
[0055] ①The value of the constant a for calculating the equivalent bulk modulus of the matrix mineral: 200;
[0056] ②Increase of equivalent bulk modulus of matrix mineral: 1.0GPa;
[0057] ③The value range of Poisson's ratio of dry rock: 0.0~0.4;
[0058] ④Increase of Poisson's ratio of dry rock: 0.02;
[0059] ⑤ Threshold: 0.5;
[0060] ⑥Sandstone flatness: 0.10, mudstone flatness: 0.04; critical porosity: 0.40.