Experimental Determination of the Biot Elastic Constant: Applications in Formation Evaluation (Sonic Porosity, Rock Strength, Earth Stresses, and Sanding Predictions)

Abstract

Summary We developed an experimental method to obtain the Biot elastic constant of rocks from laboratory dynamic and static measurements. The Biot constant often has been calculated with various empirical equations. The experimental determination of the Biot elastic constant is very important to engineering problems associated with sand control, hydraulic fracturing, wellbore stability, earth stresses, sonic porosity, and estimation of compressional-, P, and shear-, S, wave velocity. Both the dynamic and static moduli of actual reservoir sandstone core samples, jacketed and mounted in a triaxial cell under vacuum, were measured at various confining and overburden stresses. The results obtained show that the Biot constant is a complex function of porosity, permeability, pore-size distribution, and overburden and confining stress, which means that it is not really a constant. Also, the static Biot constant is greater than the dynamic one and their difference increases with increasing overburden stress according to the equation astatic =[1+0.05*(sz)ef]*adynamic (where sz is in Ksi). Moreover, both the experimental static and dynamic Biot constants may be significantly different from values calculated with empirical equations. This study suggests that quantifying the Biot constant in the laboratory may enhance the determination of rock-strength/fracturing, earth stresses, rock subsidence, sanding predictions, P- and S-wave velocities, porosity, and pore fluid from sonic and seismic data. Introduction The Biot1–7 elastic constant, a, of a rock is an important poroelastic parameter that relates stress and pore pressure and describes how compressible the dry skeletal frame is with respect to the solid material composing the dry skeletal frame of the rock. Biot1 measures the ratio of the fluid volume squeezed out to the volume change of the rock if the latter is compressed while allowing the fluid to escape. It is described as Because the petroleum-related rocks are usually saturated, it is important to know how the saturation and pore pressure affect their mechanical and flow properties. Terzaghi's8 effective-stress principle for soils states that we can obtain the effective stress by simply subtracting the fluid pressure from the total stress; i.e., se=st -ap, which means that a=1. This implies that increasing the external stress by some amount produces the same volume change of the porous material as reducing the pore pressure with the same amount. This principle appears to be valid for most properties of soils. However, in petroleum-related rocks, Terzaghi's effective-stress principle may not be valid. Then, a modified effective stress is a function of the Biot constant, a, and given by sef=st -ap. Despite the great significance of a, only a limited amount of laboratory work on its determination has been reported in the literature.9–13 The failure criteria for a saturated rock with a pore pressure are obtained by introducing the effective stress into the dry form of the failure criteria. This means that all rock failure and sand-production prediction models require a known static Biot constant value. So far, researchers, engineers, and geophysicists quite often assume that a=1 (Terzaghi's principle), which is not necessarily true. Alternatively, for the determination of a, they may use various empirical equations.14–17 These equations, however, yield different values that may vary by up to 100% or more depending on the equation used. The primary objective of this study was to determine the Biot elastic constant experimentally, both by dynamic and static measurements, and to establish a correlation between the dynamic and static a. Another objective was to identify any rock properties controlling the Biot elastic constant. Experimental Determination of the Biot Constant In this experimental method, we determined both the dynamic and static moduli of actual reservoir sandstone core samples under high vacuum (<0.15 mbar) and at various confining (s2=s3=sx) and axial (s1=sx) stresses. The vacuum was obtained and maintained in-situ while the rock sample was mounted and tested with a triaxial system. The rock sample is prepared, jacketed, and mounted in the triaxial cell. Then, the cell is closed firmly to prevent leaks and filled with the confining fluid. Vacuum is then pulled out of the sample with a high-power vacuum pump. Once the desired vacuum condition (<0.15 mbar) is established, a multistage triaxial compression test is performed, as discussed in details elsewhere.18 Axial and confining stresses were applied hydraulically. The dynamic and static data were generated at various axial and confining stresses. At each confining-stress stage, several P and S waveforms were recorded as the axial (overburden) stress was increased. The measured P - and S-wave velocities were used to calculate the dynamic Poisson's ratio and the dynamic Bulk, Shear, and Young's moduli of the dry skeletal frame of the rock, Ksk.

Optimizing Drilling Performance by Wellbore Stability and Pore-Pressure Evaluation in Deepwater Exploration

“…Significant effort, therefore, has been put into mechanically-induced instability studies. Wellbore stability is closely related to parameters [1][2][3][4][5][6][7] , such as, rock strength and deformation, in-situ stress regime, porepressure, porosity, lithological composition, clay content, water saturation, type of pore-fluids, temperature, etc. Development of a MEM (Mechanical Earth Model) 7 is essential in incorporating all the above parameters in one consistent model.…”

Experimental Determination of Elastic Anisotropy and Biot's Constant in a Saudi Arabian Reservoir Sandstone

“…Knowing the five stiffness constants from Eqs. 6 through 10, one can evaluate the two Young's moduli and the three Poisson's ratio in a transversely isotropic case (King, 1964(King, , 1969Lo et al, 1986): E V = D / (C 11 2 -C 12 2 ) (11) E H = D / (C 11 C 33 -C 13 2 ) (12) ν 1 = (C 12 C 33 -C 13 2 ) / (C 11 C 33 -C 13 2 ) (13) ν 2 = C 13 (C 11 -C 12 ) / (C 11 C 33 -C 13 2 ) (14) ν 3 = C 13 / (C 11 + C 12 )…”

“…Where D is the determinant, given by: C 11 C 12 C 13 D = C 12 C 11 C 13 (16) C 13 C 13 C 33 ν 1 = is a measure of the effect of the strain in a direction perpendicular to the axis of symmetry on the strain at right angles to it in the same plane. ν 2 = is a measure of the effect of the strain in any direction perpendicular to the axis of symmetry on the strain parallel to the axis of symmetry.…”

“…These correlations, however, yield quite different results. Besides, the variation in Biot's constant with confining pressure and axial stress (overburden and depletion effects) could be significant (Klimentos et al, 1998;Azeemuddin et al, 2001) and may need to be evaluated for each rock under consideration.…”

A New Approach to Borehole Trajectory Optimisation for Increased Hole Stability

“…When the two compressibilities are of the same order of magnitude (when the matrix is weak, and/or in zero or extremely low porosity rocks), it approaches zero. Klimentos et al (1998) showed that the effect of assuming α b equal to unity, due to lack of measured data, could lead to significant errors. When α b is measured, it is usually presented as a single value, which is true for a truly isotropic rock.…”