CROSS REFERENCE TO RELATED APPLICATIONS [0001] This application claims priority to U.S. Patent Application No. 63/036,804 filed June 9, 2020 titled “METHOD FOR ESTIMATING BIOTS COEFFICIENT FROM ADVANCED DOWNHOLE LOGGING MEASUREMENTS”, the disclosure of which is incorporated herein by reference in its entirety.

BACKGROUND

[0002] Porous formations, such as shale formations, are increasingly important for oil and gas exploration and production. Accordingly, there is a need to understand the various stresses in porous formations and to develop industrial applications based on these stresses. For example, drilling and hydraulic fracturing operations require an understanding of the stresses in porous formations.

SUMMARY

[0003] Applicant recognized the above problems and developed methods and systems disclosed herein. Embodiments are disclosed to fracture porous formations based on the effective stress of the formations. For example, methods include the steps of determining a static bulk modulus of dry rock or rock frame in a porous formation from downhole data and determining a static bulk modulus of grain minerals in the porous formation from downhole data. In embodiments Biot’s coefficient for the porous formation is determined based on the static bulk modulus of dry rock or rock frame and the static bulk modulus of grain minerals. The effective stress of the porous formation can be determined based on Biot’s coefficient for the porous formation. The effective stress can be determined by subtracting the product of Biot’s coefficient and the pore pressure from the total stress for the formation.

[0004] In embodiments, a system can include one or more hydraulic fracturing pumps to pump fluid into a porous formation at a specified pumping pressure. The system can include at least one processor and computer memory with instructions to cause the system to carry out the steps of the methods disclosed herein. When the effective stress of a porous formation is determined, the effective stress can be used to set the pumping pressure for the one or more hydraulic fracturing pumps. The effective stress of the formation can also be used in other geomechanical applications, including, for example, performing wellbore stability analyses and reservoir integrity assessments.

BRIEF DESCRIPTION OF THE DRAWINGS

[0005] Various embodiments in accordance with the present disclosure will be described with reference to the drawings, in which:

[0006] FIG. 1 illustrates an example porous formation;

[0007] FIG 2. illustrates an example method for determining a rock frame bulk modulus;

[0008] FIG 3 illustrates an example method for determining a rock frame bulk modulus;

[0009] FIG. 4 illustrates an example method for determining a grain bulk modulus;

[0010] FIG. 5 illustrates an example method for determining a grain bulk modulus; and

[0011] FIG. 6 illustrates an example method for hydraulic fracturing.

DETAILED DESCRIPTION

[0012] In the following description, various embodiments will be illustrated by way of example and not by way of limitation. References to various embodiments in this disclosure are not necessarily to the same embodiment, and such references mean at least one. While specific implementations and other details are discussed, it is to be understood that this is done for illustrative purposes only. A person skilled in the relevant art will recognize that other components and configurations may be used without departing from the scope and spirit of the claimed subject matter.

[0013] Hydraulic fracturing (sometimes referred to as “fracing” or “fracking”) is a method to extract hydrocarbons from a porous formation. The effectiveness of hydraulic fracturing depends in part on certain properties of the porous formation. For example, the failure behavior of a formation is based at least in part on the effective stress within the formation. Thus, determining the effective stress accurately is an important part of hydraulic fracturing design. Accurate determination of effective stress also plays an important role in wellbore stability analysis and reservoir integrity assessment during drilling, completion, and production.

[0014] One method of determining effective stress is based on other properties, namely the overburden stress and pore pressure of a formation. Biot’s coefficient is used to determine the effective stress based on the overburden stress and pore pressure. However, others have had difficulty accurately determining Biot’s coefficient. As a result, effective stress is not known with accuracy, which causes further problems with hydraulic fracturing design and other applications as described above. An accurate determination of Biot’ s coefficient is useful for these applications and other geomechanical applications.

[0015] Improved drilling and hydraulic fracturing methods are disclosed. Methods are disclosed for logging downhole measurements and determining Biot’s coefficient from the measurements. Downhole logging measurements include, for example, acoustic wave velocities, mineralogy, and nuclear magnetic resonance (NMR) logging.

[0016] In embodiments, acoustic velocity measurements can be used to determine the static bulk modulus of a rock frame or dry rock through frequency transfer functions and static and dynamic bulk moduli correlation functions. The acoustic velocity measurements can first be used to determine the dynamic bulk modulus. The dynamic values can then be corrected for frequency if measured by ultrasonic or seismic devices or logs. The resulting values can then be converted to static values.

[0017] In embodiments, mineral ogical and NMR measurements are used to determine the grain bulk modulus. The mineralogical measurements can be used to determine the volumetric amount of each mineral, and the NMR measurements can be used to determine diagenetic clay characteristics.

[0018] In Equation 1, Biot’s coefficient (α) is employed to determine the effective stress (s') from the externally applied total stress (a) and pore pressure (P _p ):,

[0019] In Equation 1, K _frmS represents the static bulk modulus of the rock frame or dry rock and K _{g S} represents the static bulk modulus of grains, or forming minerals, from which the rock is composed. Symbol ‘S’ in Equation 1 denotes “static.”

[0020] Figure 1 illustrates an exemplary porous formation 100. The externally applied total stress (s) 102 is shared between the rock frame and pore fluids in subsurface formations. The external total stress 102 can represent any applied exterior stress. When a fluid fills the pores, or voids, in a porous formation 100, the pressure in the pores can be measured. The pore pressure (P _p ) 106 in Figure 1 represents this pressure.

[0021] The effective stress (σ'' 104, not externally applied total stress (a) 102, governs the failure behavior of subsurface formations. The effective stress 104 represents the stress acting on the grain and plays an important role in many engineering applications, including wellbore stability analysis, hydraulic fracturing design, and reservoir integrity assessment during drilling, completion, injection, and production. Therefore, the determination of Biot’s coefficient has multiple geomechanical applications.

[0022] Methodologies disclosed herein allow for the determination of Biot’s coefficient in the field using downhole logging measurements, which improves the design and performance of engineering operations, such as drilling and hydraulic fracturing. The determination can be made in near real time, and adjustments to drilling, completion, and stimulation operations can occur in the field.

[0023] Returning to Equation 1, Biot’s coefficient can be determined by first determining two moduli: the static bulk modulus of the rock frame {K _frmS ) and the static bulk modulus of grain minerals (K _{g s} ). These two values can be measured through hydrostatic compression tests on the rock sample in the laboratory: rock frame bulk modulus can be measured by performing a hydrostatic compression test on a jacketed sample with zero pore pressure, and grain mineral bulk modulus can be measured by performing a hydrostatic compression test on an unjacketed sample, where the pore and confining pressure are equal. Values from hydrostatic testing are called static bulk moduli.

[0024] According to embodiments, the static bulk moduli can be determined using downhole measurements in the field, without lab testing. For example, downhole acoustic velocity measurements provide a way to determine the dynamic bulk modulus (K _{DH ,d Sat} ) of a rock at fully fluid- saturated condition (denoted as Sat) as provided in Equation 2:

[0025] In Equation 2, V _{P Sat} and V _{S Sat} respectively represent the compressional and shear wave velocities respectively, and p _B ,sat represents the bulk density of a fully fluid- saturated rock,. Symbol d denotes dynamic. Formation rock can be saturated with various types of fluids ( e.g water, oil and water, gas and water; or gas, oil, and water) at downhole in-situ conditions of high temperature and pressures, resulting in the frequency dispersion effect on downhole wave velocity measurements. In embodiments, Equation 2 is applied to the dry frame because rock mechanical properties are intrinsic properties and should not be influenced by fluids in the void. Fluid stripping can be performed and the dry frame’s acoustic velocity can be used for the calculations. [0026] According to embodiments, the downhole dynamic bulk modulus (K _{DH d Sat} ) needs frequency correction and a transform applied to the rock frame bulk modulus (K _{frm s} ) for the determination of Biot’s coefficient. In embodiments, transform functions are used to convert downhole dynamic bulk modulus to static rock frame bulk modulus. Figure 2 illustrates such an exemplary method.

[0027] Figure 2 illustrates an exemplary method 200 for determining the static bulk modulus of rock frame. As illustrated in Figure 2, the dynamic bulk moduli of the rock frame is measured 202. The function K _{Lab d Sat} = f(freq, fluid) is built based on the previously determined dynamic bulk moduli 204. This function is then transformed to downhole function K _{DH d Sat} = f(freq, fluid) 206. Extrapolating the downhole function yields the downhole function at low frequency 208. The dynamic bulk modulus K _{DH d jrm} of rock frame at dry conditions can then be derived 210. The method further includes measuring the static bulk modulus and dynamic bulk modulus of the rock frame 212 and building the correlation between the static bulk modulus and dynamic bulk modulus of the rock frame 214. The static bulk modulus of the rock frame K^ _{rm S} can be determined based on the correlation between K _{Lab djrm} and K _frmS built from laboratory measurements and the determined dynamic bulk modulus 216. [0028] Figure 3 illustrates another exemplary method 300 for determining the static bulk modulus of rock frame. First, the dynamic bulk moduli of the rock frame K _{Lab d Sat} is measured under various frequencies 302. A function K _{Lab d Sat} = / (freq, fluid) is built from laboratory dynamic bulk modulus measurements on fully saturated rocks to form the relationship between the dynamic bulk modulus and frequencies applied 304. Then, the function K _{Lab dSat} = f (freq, fluid) is transformed or adjusted to obtain the downhole function K _{DH d Sat} = f (freq, fluid) based on downhole dynamic bulk modulus measurements, which corrects the effect of other factors such as temperature and pressure on the bulk modulus 306.

[0029] Dynamic bulk modulus K _{DH d Sat} at low frequency (LF) is extrapolated from the function K _DH,d,sat = f(frecq, fluid) 308. Then, the dynamic bulk modulus K _{DH djrm} of rock frame or dry rock can be derived using Gassmann fluid substitution model or other available models 310. The dynamic bulk modulus K _{DH djrm} of rock frame or dry rock is assumed to be the same as laboratory dynamic bulk modulus K _{Lab djrm.} The dynamic bulk modulus K _{Lab djrm} and static bulk modulus Kf _rm, s °f the rock frame are measured, for example based on hydrostatic and ultrasonic tests 312. The correlation between K _{Lab djrm} and K _frmS is built from these measurements 314. The static bulk modulus K _{frm S} is obtained from the correlation 316.

[0030] Figure 4 illustrates an exemplary method 400 to determine the grain bulk modulus K _{g S.} Exemplary methods include the step of determining an initial mineral fraction f contributing to the grain bulk modulus 402 and measuring the volumetric fraction of the various minerals 404. The static grain mineral bulk modulus can be determined from based on these inputs 406. In addition, the static grain bulk modulus can be measured 408. The difference between the calculated static grain mineral bulk modulus and the measured static grain bulk modulus can be minimized 410, as discussed in more detail below. A determination can be made if this difference is minimized 412. If so, the grain bulk modulus can be derived 416. If not, the mineral fraction f _t can be adjusted 414. The static grain mineral bulk modulus can be redetermined 406 and the difference between the determined static grain mineral bulk modulus and the measured static grain bulk modulus attempted to be minimized again 410. This loop can be repeated until a minimum is achieved and the grain bulk modulus derived 416.

[0031] Figure 5 illustrates another exemplary method 500 to determine the grain bulk modulus. In embodiments, the initial mineral fraction f _t contributing to the grain bulk modulus is determined based on NMR or other sources 502. The minerology and volumetric fraction V _L or weight fraction of minerals in the formation can be measured with X-ray diffraction XRD analysis on samples from the formation. In other words, the grain bulk modulus may have as an input the mineral ogical measurements or XRD analysis on rock materials (core samples or drill cuttings for example) and their diagenetic clay characteristics. Such embodiments represent an improvement, among other reasons, due to the difficulties measuring and calculating the grain bulk modulus, particularly for certain rock mineralogical compositions.

[0032] The grain bulk modulus of a mixture of different minerals is often difficult to determine. Such a determination can depend on the volumetric fraction of each mineral component. For example, the Voigt (V) and Reuss (R) methods give the theoretical maximum (K _{g S} ) and minimum grain bulk moduli of a mixture of n minerals as

Equation 3

Equation 4 where i = 1,2, ··· , n are the minerals composed of rock; V _L and K _t are the volumetric fractions and bulk modulus of the i ^th mineral, respectively. Symbol ‘s’ denotes static bulk modulus. [0033] In embodiments, the grain bulk modulus is determined based on the characteristics of each mineral in the rock as: where K _{g s} , _mea , _j is the measured static grain bulk moduli of /th core sample, and f, is the fraction (0 to 1.0) of mineral i which contributes to the grain bulk modulus of a mineral mixture. Clay minerals can be diagenetic and detrital (or structural). Their contribution to grain bulk moduli may be different, represented by the fraction of each mineral f _i . The fraction of diagenetic minerals can be derived from NMR logging measurements and neutron-gamma spectroscopy downhole measurements (or XRD analysis).

[0034] The static grain mineral bulk modulus K _g,s,caI can be determined based on the mineral fraction and volumetric fraction 506, for example using Equation 5 above. The static grain bulk modulus K _g , _{S mea} can be based on hydrostatic tests 508. The difference between thee measured and calculated values, respectively K _{g s mea} and K _{g s cai} can be minimized 510, for example using Equation 6 above. If confirmed to be minimized 512, then the grain bulk modulus K _{g S} has been determined based on downhole mineralogy measurements 516. If not minimized, then the mineral fraction can be adjusted 514, and the loop repeated as described above until the difference between K _g,s, mea and K _g , _s , _cal is minimized.

[0035] The disclosed methods provide accurate determinations of the static bulk modulus of the rock frame or dry rock K^ _{rm S} and the static bulk modulus of grains K _{g S} . Biot’s coefficient can be derived from these determinations, for example by using Equation 1 above. Biot’s coefficient can be used in numerous geomechanical applications. [0036] Figure 6 illustrates an exemplary hydraulic fracturing method 600. In the method, the static bulk modulus of a rock frame K _frmS is determined 602, for example using the methods described herein. The static bulk modulus of grains K _{g ,s} is determined 604, for example using the methods described herein. Biot’s coefficient is determined based on these determinations 606, for the formation can then be determined 608, for example using Equation 1. In other words, the product of Biot’s coefficient and the pore pressure is subtracted from the total or overburden stress to yield the effective stress of the formation.

[0037] The effective stress yields important information about the failure behavior of the formation. For instance, the pressure at which fluid must be pumped into the formation is based on the effective stress. Therefore, the formation can be fractured based on the effective stress of the formation 610 determined as described herein. For example, hydraulic fracturing pumps can be set at specified pressure settings based on the effective stress of the formation, thereby improving hydraulic fracturing operations. Other geomechanical applications are also improved with improved information about the effective stress of the formation. For example, determining the effective stress using the methods disclosed herein can be used to analyze wellbore stability or assess reservoir integrity during drilling, completion, and production.

[0038] Various steps discussed herein may be carried out as part of a computer implemented method. In embodiments, computer memory may include instructions that, when executed by at least one processor, carry out the steps of methods disclosed herein. Exemplary systems may include such memory and at least one processor. Exemplary systems may also include such memory and at least one processor and one or more hydraulic fracturing pumps. In embodiments, pressure settings for the one or more hydraulic fracturing pumps are set based on the effective pressure of a formation, which is determined based on methods disclosed herein.

[0039] Although the technology herein has been described with reference to particular embodiments, it is to be understood that these embodiments are merely illustrative of the principles and applications of the present technology. It is therefore to be understood that numerous modifications may be made to the illustrative embodiments and that other arrangements may be devised without departing from the spirit and scope of the present technology as defined by the appended claims.