MEASURING FORMATION POROSITY AND PERMEABILITY

Abstract

Values for porosity and permeability of core samples in a borehole are estimated by generating radial waves with an acoustic source in fluid around the core sample, and measuring pressure in the fluid. Moreover, the acoustic source operates at frequency close to a resonant frequency of the core sample. After the acoustic source no longer operates at the resonant frequency, pressure in the fluid attenuates over time. The pressure attenuation is recorded by the pressure measurements, along with the pressure in the fluid at the first harmonic (spectral component). The pressure attenuation and spectral component each are dependent on porosity and permeability of the core sample. Thus values for the porosity and permeability are determined based on the arithmetic relationships between pressure attenuation and the spectral component and porosity and permeability.

Claims

1. A method of estimating characteristics of a core sample in a borehole comprising (a) measuring values of pressure in a fluid surrounding the core sample that fluctuate in response to radial waves that propagate through the fluid and the core sample; (b) identifying an attenuation of pressure in the fluid over time based on step (a); (c) identifying a pressure in the fluid at the frequency of a resonant peak; and (d) estimating at least one of porosity or permeability of the core sample based on the values identified in steps (b) and (c).

2. The method of claim 1, wherein the radial waves are at a frequency at about a resonant frequency of the core sample.

3. The method of claim 2, wherein the step of identifying an attenuation of pressure in the fluid takes place after radial waves at the resonant frequency of the core sample no longer propagate through the fluid.

4. The method of claim 1, further comprising providing an acoustic source for generating the radial waves.

5. The method of claim 1, wherein the acoustic source operates at about the first resonant frequency of the core sample.

6. The method of claim 1, wherein the resonant peak of step (c) is at the first resonant peak.

7. The method of claim 1, further comprising obtaining the core sample with a coring tool, and wherein step (a) takes place in the coring tool.

8. A method of estimating characteristics of a core sample in a borehole comprising: (a) obtaining the core sample with a coring tool in the borehole and immersing the core sample with a liquid; (b) generating radial waves in the liquid with a transmitter that operates at a resonant frequency of the core sample; (c) identifying an attenuation of peaks of pressure in the fluid over time that are from the radial waves that propagate through the core sample; (d) identifying a pressure in the fluid at the frequency of a resonant peak that is from at least one of the radial waves; and (e) estimating at least one of porosity or permeability of the core sample.

9. The method of claim 8, wherein the step of estimating at least one of porosity or permeability of the core sample is based on the steps of identifying of steps (b) and (c).

10. The method of claim 8, further comprising the generation of radial waves at the resonant frequency is terminated prior to steps (c) and (d).

11. The method of claim 8, wherein the resonant frequency of the core sample is estimated.

12. A system for estimating characteristics of a core sample in a borehole comprising a coring bit having an axial bore in which the core sample is selectively disposed; a tubular transducer that circumscribes the axial bore and selectively oscillates at about the resonant frequency of the core sample; an acoustic receiver disposed in the axial bore and immersed in liquid that contacts the core sample, and that selectively receives radial waves generated by the transducer that reflect from the core sample; and a processor for estimating a value for at least one of a porosity or permeability of the core sample based on a measurement of an attenuation of pressure of the radial waves received by the receiver and a pressure measured in the fluid at a first resonant frequency of the core sample.

13. The system of claim 12, further comprising springs for suspending the receiver in the fluid.

14. The system of claim 12, further comprising a housing on which the coring bit is attached and which defines a coring tool, and a drill string attached to an end of the coring tool for rotating the coring bit.

15. The system of claim 12, wherein the liquid is disposed in the axial bore when the coring bit is at surface.

Description

BRIEF DESCRIPTION OF DRAWINGS

[0007] Some of the features and benefits of the present invention having, been stated, others will become apparent as the description proceeds when taken in conjunction with the accompanying drawings, in which:

[0008] FIG. 1 is a side partial sectional view of an example of a coring tool deployed in a wellbore and in accordance with the present invention.

[0009] FIG. 2 is a side sectional view of an example of a portion of the coring tool of FIG. 1 obtaining core samples in a wellbore and in accordance with the present invention.

[0010] FIG. 3 is a side sectional perspective view of an embodiment of acoustically measuring a core sample in the coring tool of FIG. 2 and in accordance with the present invention.

[0011] FIG. 4 is an example plot of pressure wave attenuation of an acoustic signal adjacent a core sample and in accordance with the present invention.

[0012] FIG. 5 is an example plot of eigenwaves of an acoustic signal adjacent a core sample and in accordance with the present invention.

[0013] While the invention will be described in connection with the preferred embodiments, it will be understood that it is not intended to limit the invention to that embodiment. On the contrary, it is intended to cover all alternatives, modifications, and equivalents, as may be included within the spirit and scope of the invention as defined by the appended claims.

DETAILED DESCRIPTION OF INVENTION

[0014] The method and system of the present disclosure will now be described more fully hereinafter with reference to the accompanying drawings in which embodiments are shown. The method and system of the present disclosure may be in many different forms and should not be construed as limited to the illustrated embodiments set forth herein, rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey its scope to those skilled in the art. Like numbers refer to like elements throughout

[0015] It is to be further understood that the scope of the present disclosure is not limited to the exact details of construction, operation, exact materials, or embodiments shown and described, as modifications and equivalents will be apparent to one skilled in the art. In the drawings and specification, there have been disclosed illustrative embodiments and, although specific tetras are employed, they are used in a generic and descriptive sense only and not for the purpose of limitation.

[0016] Shown in side sectional view in FIG. 1 is an example of a coring tool 10 disposed in a borehole 12 that intersects a subterranean formation 14. The coring tool 10 is deployed on a lower end of a drill string 16. A rotary table 18 couples with the drill string 16 on surface and rotates the drill string 16 so the coring tool 10 can bore into the formation 14 and obtain core samples from the formation 14. The rotary table 18 is included as part of a drilling rig 20 on the surface shown having a derrick 21 for supporting the rotary table 18 and drill string 16.

[0017] Referring now to FIG. 2, shown in side sectional view is a detailed portion of the coring tool 10 in the borehole 12 illustrating that a bit 22 is included on a lower end of the coring tool 10 for excavating in the formation 14. An elongate bore 24 extends axially through the bit 22, and in which a core sample 26 is received after the bit 22 excavates the rock surrounding the core 26. Drilling fluid or mud can be circulated downhole through the coring tool 10 within an annulus 28 shown extending axially through the bit 22 and spaced radially outward from the bore 24. In addition to cooling the bit 22 during excavating, the drilling fluid washes cuttings 30 back uphole that were removed from the formation 14 by the cutting action of the bit 22.

[0018] Optional bearings 32 are shown in the annulus 28 and between an inner body 34 and outer tubing 35. Both the inner body 34 and outer tubing 35 are tubular members, where the inner body 34 is coaxially and rotationally set within the outer tubing 35. The inner body 34 contains the core sample 26, and the outer tubing 35 defines an outer surface of the annulus 28 and includes cutting elements on its outer surface for cutting through the formation 14.

[0019] Also disposed in the annulus 28 is a transmitter 36 for generating an acoustic signal; in one example the signal includes a radial wave R (FIG. 3) that contacts, passes through, and reflects from the sample 26. A receiver 38, which in one example can be a hydrophone, is illustrated spaced axially away from the transmitter 36. Examples exist wherein the transmitter 36 and receiver 38 are integrated into a single unit. In an example, the inner body 34 is filled with a selected known fluid before being deployed downhole. Optionally, an elastic liner (not shown) can be provided on an inner surface of the inner body 34 for preventing damage to the sample core 26 and to prevent formation fluids from escaping the sample core 26 or from within the bore 24.

[0020] Schematically illustrated in side sectional perspective view in FIG. 3 are radial acoustic waves R being generated in fluid F that surrounds the core sample 26, and where reflections of the radial waves R from the core sample 26 are recorded by the receiver 38. As illustrated in the example of FIG. 3, the transmitter 36 is a tubular member that is selectively oscillated to generate the radial acoustic waves R in the fluid F. Example materials for the transmitter 36 include piezoelectric and/or magnetostrictive materials, in an example the transmitter 36 is oscillated by applying electricity to the transmitter 36.

[0021] Also shown in FIG. 3 is that the hydrophone 38 is suspended in the annulus 28 and supported by springs 42. In a non-limiting example, the transmitter 36 is operated at the resonant frequency of the radial eigenwaves of the core sample 26; which in one embodiment is at a frequency that is within about +/5% of the first resonant frequency. Further in this example, the receiver 38 records the pressure of the radial acoustic wave R. Where in one example, the pressure of the radial acoustic wave is the increase in pressure due to the compressional waves.

[0022] Referring back to FIG. 2, in one non-limiting example of operation, the coring tool 10 is rotated in the borehole 12 so that the bit 22 excavates an annular cut into the formation 14 thereby forming the core sample 26 that extends into the bore 24. The transmitter 36 is actuated for a period of time to generate radial waves R (FIG. 3) shown in the fluid F propagating towards the core sample 26. A portion of the radial waves R reflect from the core sample 26 and are measured by the receiver 38; as indicated above the receiver 38 records the pressure of the radial acoustic wave R reflected from the core sample 26.

[0023] Referring now to FIG. 4, shown is a plot 44 that represents pressure in the fluid F measured by the receiver 38, where the recorded pressure is a result of the acoustic signal in the fluid F generated by the transmitter 36. Shown in FIG. 4 is that the peaks 46 of the measured pressure decay over time. The decay of the pressure peaks 46 are modeled by plot 48, which is represented by the expression Ae.sup.(t). In one embodiment, the measured pressure decay over time takes place after the transmitter 36 no longer operates at the resonant frequency. It has been found that the attenuation factor of the pressure waves (the envelope of the maximal values of pressure) is also determined by permeability k and porosity of the formation, where:

=(k,)Equation (1).

[0024] Shown in FIG. 5 is a plot 50 that represents a spectrum of sample eigenwaves.

[0025] Analyzing the behavior of spectral component of pressure at the harmonic source at the minimal resonant frequency can yield useful information. For example, it has been found that the value of pressure at the first maximum (pressure at the frequency of the first resonant peak) is dependent upon formation permeability k and porosity co, where:

p.sub.max=(k,)Equation (2).

[0026] Equations (1) and (2) at given acoustic velocities in fluid and the formation, at corresponding porosities, constitute a set of equations for finding permeability and porosity. Measuring p.sub.max and and solving sets (1), (2), values can be obtained for permeability k and porosity . Provided below are how the functional dependences (k, ) and p.sub.max=(k, ) were theoretically obtained.

[0027] Functions (k, ), (k, ) are found from the theory describing the radial waves in the borehole with saturated porous core sample. The acoustic equations of the two-velocity medium are a linearized version of the following equations, where Equations (3)-(7) describe the propagation of waves in the saturated porous core sample:

[00001]$\begin{array}{cc}.Math.\frac{}{t}+\mathrm{div}.Math..Math.j=0,,& \mathrm{Equation}.Math..Math.\left(3\right)\\ .Math.\frac{g_{\mathrm{ik}}}{t}+g_{\mathrm{jk}}.Math..Math.{}_i.Math.u_j+g_{\mathrm{ij}}.Math..Math.{}_k.Math.u_j+u_j.Math.{}_j.Math.g_{\mathrm{ik}}=0,,& \mathrm{Equation}.Math..Math.\left(4\right)\\ .Math.\frac{j_i}{t}+{}_k.Math.\left({}_s.Math.u_i.Math.u_k+{}_l.Math.v_i.Math.v_k+p.Math..Math.{}_{\mathrm{ik}}+g_{\mathrm{jk}}.Math.h_{\mathrm{ij}}+{}_{\mathrm{ik}}\right)=0,,& \mathrm{Equation}.Math..Math.\left(5\right)\\ \frac{v_i}{t}+\left(v,\right).Math.v_i=-\frac{1}{}.Math.{}_i.Math.p+\frac{{}_s}{2.Math.}.Math.{}_i.Math.{\left(u-v\right)}^2-\frac{1}{2.Math.}.Math.h_{\mathrm{jk}}.Math.{}_i.Math.g_{\mathrm{jk}}+\frac{1}{{}_l}.Math.f_i,,\mathrm{and}& \mathrm{Equation}.Math..Math.\left(6\right)\\ .Math.\frac{S}{t}+\mathrm{div}\left(S.Math.\frac{j}{}+\frac{q}{T}\right)=\frac{R}{T},.& \mathrm{Equation}.Math..Math.\left(7\right)\end{array}$

[0028] The condition of the positive dissipative function R is as follows:

[00002]$\begin{array}{cc}-R=f.Math.\frac{1}{{}_l}.Math.\left(j-.Math..Math.u\right)+q.Math.\frac{T}{T}+A_{\mathrm{ik}}\left({}_i.Math.v_k+{}_k.Math.v_i-\frac{2}{3}.Math.{}_{\mathrm{ik}}.Math..Math.\mathrm{div}.Math..Math.v\right)+a.Math..Math.\mathrm{div}.Math..Math.v,& \mathrm{Equation}.Math..Math.\left(8\right)\end{array}$

[0029] and it determines the form of the dissipative flows:

[00003]$\begin{array}{cc}.Math.q={}_{12}\left(.Math..Math.u-j\right)+.Math.\frac{1}{T}.Math.T,,& \mathrm{Equation}.Math..Math.\left(9\right)\\ .Math.f+{}_k.Math.{}_{\mathrm{ik}}=b\left(.Math..Math.u-j\right)+{}_l.Math.{}_{12}.Math.\frac{1}{T}.Math.T,,\mathrm{and}& \mathrm{Equation}.Math..Math.\left(10\right)\\ {}_{\mathrm{ik}}=A_{\mathrm{ik}}+a.Math..Math.{}_{\mathrm{ik}},A_{\mathrm{ik}}=-\left({}_i.Math.v_k+{}_k.Math.v_i-\frac{2}{3}.Math.{}_{\mathrm{ik}}.Math..Math.\mathrm{div}.Math..Math.v\right),a=-.Math..Math.\mathrm{div}.Math..Math.v..& \mathrm{Equation}.Math..Math.\left(11\right)\end{array}$

[0030] The equation of the state of the saturated porous medium is selected in the following form:

[00004]$\begin{array}{cc}p=p_0-\frac{1}{2}.Math.{\mathrm{Kg}}_{\mathrm{ll}}+\frac{1}{{}_0}.Math.\left(K+\right).Math.+{}_p.Math.S+\frac{1}{2}.Math.{{}_s\left(u-v\right)}^2,,& \mathrm{Equation}.Math..Math.\left(12\right)\\ .Math.T=T_0+{}_s.Math.S+{}_p.Math.\frac{1}{{}_0}.Math.,,\mathrm{and}& \mathrm{Equation}.Math..Math.\left(13\right)\\ h_{\mathrm{ik}}=h_{0,\mathrm{ik}}+.Math..Math.g_{\mathrm{ik}}+\frac{1}{2}.Math..Math..Math.g_{\mathrm{ll}}.Math.{}_{\mathrm{ik}}-\frac{}{{}_0}.Math.{\mathrm{Kg}}_{\mathrm{ll}}.Math.{}_{\mathrm{ik}}-\frac{1}{2}.Math.g^{\mathrm{ik}}.Math.{{}_s\left(u-v\right)}^2..& \mathrm{Equation}.Math..Math.\left(14\right)\end{array}$

[0031] In the formulae above, .sub.l, .sub.s are partial densities of the saturating fluid and porous matrix, respectively, =.sub.l+.sub.s is density of the saturated medium; u, v are velocities of the porous matrix and saturating fluid; j=.sub.su+.sub.lv is the momentum; h.sub.ik is the stress tensor, g.sub.ik is the metric tensor; p is pressure, T is temperature, S is entropy; , and are elastic moduli of the saturated porous medium; .sub.s, .sub.p are thermodynamic parameters of the medium. The dissipative coefficient of the interphase friction b may be found from the ratio b=/(k), where is dynamic viscosity of saturating fluid, k is permeability of the porous medium.

[0032] The linearized isothermal version of Equations (3)-(6) enables one to find the motion equations of the porous matrix u and saturating fluid v in the acoustic approximation:

[00005]$\begin{array}{cc}\frac{{}^2.Math.u}{t^2}-c_t^2.Math..Math..Math.u-a_1.Math.\mathrm{div}.Math..Math.u+a_2.Math.\mathrm{div}.Math..Math.v+\frac{{}_l}{{}_s}.Math.b\left(\frac{u}{t}-\frac{v}{t}\right)=0,,\mathrm{and}& \mathrm{Equation}.Math..Math.\left(15\right)\\ .Math.\frac{{}^2.Math.v}{t^2}+a_3.Math.\mathrm{div}.Math..Math.u-a_4.Math.\mathrm{div}.Math..Math.v-b\left(\frac{u}{t}-\frac{v}{t}\right)=0..& \mathrm{Equation}.Math..Math.\left(16\right)\end{array}$

[0033] Coefficients a.sub.i (1=1, . . . , 4) are determined by the elastic moduli of the saturated porous medium , and :

[00006]$\begin{array}{cc}{a}_1=\frac{1}{{}_s}.Math.\left(\frac{{}_s^2}{{}^2}.Math.+\frac{{}_l^2}{{}^2}.Math.K+\frac{1}{3}.Math.\right),a_2=\frac{{}_l}{{}_s}.Math.\left(\frac{{}_l}{{}^2}.Math.K-\frac{{}_s}{{}^2}.Math.\right),,\mathrm{and}& \mathrm{Equation}.Math..Math.\left(17\right)\\ .Math.a_3=\frac{{}_k}{{}^2}.Math.K-\frac{{}_s}{{}^2}.Math.,a_4=\frac{{}_l}{{}^2}.Math.K+\frac{{}_l}{{}^2}.Math.,.& \mathrm{Equation}.Math..Math.\left(18\right)\end{array}$

[0034] The three moduli K, and , where K=+2/3, may be found from three acoustic velocities c.sub.p1, c.sub.p2, c.sub.i of the ideal hydrodynamic approximation, in accordance with the formulae given below:

[00007]$\begin{array}{cc}={}_s.Math.c_t^2,=\frac{1}{2}.Math.\frac{{}_s}{{}_l}.Math.\left(.Math..Math.c_{p.Math..Math.1}^2+.Math..Math.c_{p.Math..Math.2}^2-4.Math.{}_l.Math.c_t^2-\sqrt{{\left(.Math..Math.c_{p.Math..Math.1}^2-.Math..Math.c_{p.Math..Math.2}^2\right)}^2-\frac{64}{9}.Math.{}_s.Math.{}_l.Math.c_t^4}\right),,\mathrm{and}& \mathrm{Equation}.Math..Math.\left(19\right)\\ =\frac{1}{2}.Math.\left(.Math..Math.c_{p.Math..Math.1}^2+.Math..Math.c_{p.Math..Math.2}^2-\frac{8}{3}.Math.{}_s.Math.c_t^2+\sqrt{{\left(.Math..Math.c_{p.Math..Math.1}^2-.Math..Math.c_{p.Math..Math.2}^2\right)}^2-\frac{64}{9}.Math.{}_s.Math.{}_l.Math.c_t^4}\right)..& \mathrm{Equation}.Math..Math.\left(20\right)\end{array}$

[0035] The acoustic field in borehole fluid is described by the following equation:

[00008]$\begin{array}{cc}\frac{{}^2.Math.v}{t^2}-c_{p.Math..Math.0}^2.Math.\mathrm{div}.Math..Math.v=0,.& \mathrm{Equation}.Math..Math.\left(21\right)\end{array}$

[0036] Where c.sub.p0 is the velocity of sound in borehole fluid.

[0037] The following boundary conditions are shown in FIG. 3: (1) On the surface of the acoustic transmitter, the source with the amplitude p.sup.(0)=p.sub.0(t) is given; (2) At the boundary between the borehole and saturated porous sample, the continuity condition is applied to the mass flow, stress tensor .sub.rr, and partial pressures: (1)u.sub.r.sup.(1)+v.sub.r.sup.(1)=v.sub.r.sup.(0), .sub.rr.sup.(1)=.sub.rr.sup.(0), p.sup.(1)/.sup.(1)=p.sup.(0)/.sub.(0), where is porosity. The latter condition follows from a more general continuity condition of the partial stress tensor in fluid .sub.rr.sup.l(1)=.sub.rr.sup.(0), where .sub.rr.sup.s(1)+.sub.rr.sup.l(1).

[0038] Time derivatives of stress tensor components and pressure can be found via the displacement velocities:

[00009]$\begin{array}{cc}.Math.-\frac{p^{\left(1\right)}}{t}=d_1.Math..Math.\mathrm{div}.Math..Math.u^{\left(1\right)}+d_2.Math..Math.\mathrm{div}.Math..Math.v^{\left(1\right)},,& \mathrm{Equation}.Math..Math.\left(22\right)\\ .Math.\frac{{}_{\mathrm{rr}}^{\left(1\right)}}{t}=d_3.Math..Math.\mathrm{div}.Math..Math.u^{\left(1\right)}+d_4.Math..Math.\mathrm{div}.Math..Math.v^{\left(1\right)}+2.Math.d_5.Math..Math.u_r^{},,\mathrm{and}& \mathrm{Equation}.Math..Math.\left(23\right)\\ .Math.\frac{{}_{\mathrm{rr}}^{\left(0\right)}}{t}=-\frac{p^{\left(0\right)}}{t}=d_0.Math..Math.\mathrm{div}.Math..Math.v^{\left(0\right)},.& \mathrm{Equation}.Math..Math.\left(24\right)\\ \mathrm{Where}.Math.\text{:}.Math..Math.d_1=\left({}_s^{\left(1\right)}.Math.-{}_l^{\left(1\right)}.Math.K\right).Math.\text{/}.Math.{}^{\left(1\right)},d_2=\left({}_l^{\left(1\right)}.Math.+{}_l^{\left(1\right)}.Math.K\right).Math.\text{/}.Math.{}^{\left(1\right)},& \mathrm{Equation}.Math..Math.\left(25\right)\\ .Math.d_3={}_s^{\left(1\right)}.Math..Math.\text{/}.Math.{}^{\left(1\right)}-2.Math..Math.\text{/}.Math.3,d_4={}_l^{\left(1\right)}.Math..Math.\text{/}.Math.{}^{\left(1\right)},,\mathrm{and}& \mathrm{Equation}.Math..Math.\left(26\right)\\ .Math.{}_5=,{}_0={}^{\left(0\right)}.Math.c_{p.Math..Math.0}^2..& \mathrm{Equation}.Math..Math.\left(27\right)\end{array}$

[0039] The set of equations presented with corresponding boundary conditions was solved numerically. The result is the functional dependences (k, ), (k, ).

[0040] Referring back to FIG. 1, a controller 52 is shown in communication with the coring tool 10, wherein the controller 52 can be housed within the coring tool 10 or on surface. In an example, the controller 52 receives the values of pressure measured in the coring tool 10. Yet further optionally, the controller 52 is equipped with a processor for calculating the values of porosity and/or permeability as described above. The controller 52 may include an information handling system (IHS), where the IHS optionally includes one or more of a processor, memory accessible by the processor, nonvolatile storage area accessible by the processor, and logics for performing each of the steps above described.

[0041] The present invention described herein, therefore, is well adapted to carry out the objects and attain the ends and advantages mentioned, as well as others inherent therein. While a presently preferred embodiment of the invention has been given for purposes of disclosure, numerous changes exist in the details of procedures for accomplishing the desired results. These and other similar modifications will readily suggest themselves to those skilled in the art, and are intended to be encompassed within the spirit of the present invention disclosed herein and the scope of the appended claims.