You must log in to edit PetroWiki. Help with editing

Content of PetroWiki is intended for personal use only and to supplement, not replace, engineering judgment. SPE disclaims any and all liability for your use of such content. More information

# Difference between revisions of "Borehole gravimetry"

Line 1: | Line 1: | ||

− | Borehole gravity was pioneered by Smith<ref name="r1" /> and then applied to problems of reservoir evaluation by McCulloh ''et al.''<ref name="r2" /> The borehole gravity meter or gravimeter responds to variations in density. Modern instruments sense a rock volume that is approximately the same as that investigated by deep [[ | + | Borehole gravity was pioneered by Smith<ref name="r1">Smith, N.L. 1950. The Case for Gravity Data From Boreholes. Geophysics 15 (4): 605-636. http://dx.doi.org/ 10.1190/1.1437623.</ref> and then applied to problems of reservoir evaluation by McCulloh ''et al.''<ref name="r2">McCulloh, T.H., Kandle, G.R., and Schoellhamer, J.E. 1968. Application of Gravity Measurements in Wells to Problems of Reservoir Evaluation. Trans., Soc. of Professional Well Log Analysts 9th Annual Logging Symposium, New Orleans, paper O.</ref> The borehole gravity meter or gravimeter responds to variations in density. Modern instruments sense a rock volume that is approximately the same as that investigated by deep [[Resistivity_and_spontaneous_(SP)_logging|resistivity tools]]. Unlike the shallower-sensing [[Density_logging|density log]], the borehole gravimeter is insensitive to wellbore conditions such as rugosity and the presence of casing. Its principal applications are: |

− | |||

− | |||

− | == Theory == | + | *Through-casing time-lapse monitoring of saturations/fluid contacts in gas reservoirs |

+ | *Downhole calibration of surface geophysical mapping of geological structures | ||

+ | |||

+ | == Theory == | ||

Any two masses, ''m''<sub>''i''</sub> and ''m''<sub>''j''</sub>, separated by a distance, ''r'', experience a gravitational force of attraction, ''f'', which is expressed as | Any two masses, ''m''<sub>''i''</sub> and ''m''<sub>''j''</sub>, separated by a distance, ''r'', experience a gravitational force of attraction, ''f'', which is expressed as | ||

− | [[File:Vol5 page 0409 eq 001.png]]....................(1) | + | [[File:Vol5 page 0409 eq 001.png|RTENOTITLE]]....................(1) |

where ''G'' is the universal gravitational constant (6.6726 × 10 –8 cgs units). More specifically, a mass, ''m'', on the surface of the Earth would experience a gravitational force given by | where ''G'' is the universal gravitational constant (6.6726 × 10 –8 cgs units). More specifically, a mass, ''m'', on the surface of the Earth would experience a gravitational force given by | ||

− | [[File:Vol5 page 0409 eq 002.png]]....................(2) | + | [[File:Vol5 page 0409 eq 002.png|RTENOTITLE]]....................(2) |

where ''M'' is the mass of the Earth, ''R'' is its radius, and an acceleration due to gravity, ''g'', is given by | where ''M'' is the mass of the Earth, ''R'' is its radius, and an acceleration due to gravity, ''g'', is given by | ||

− | [[File:Vol5 page 0409 eq 003.png]]....................(3) | + | [[File:Vol5 page 0409 eq 003.png|RTENOTITLE]]....................(3) |

− | Because the Earth is a rotating oblate spheroid, the quantity g at mean sea level varies with latitude, and it must be corrected for tidal effects. The unit of g is the Gal [1 cm/s<sup>2</sup>]. Surface gravity surveys use the milliGal as the preferred unit. Borehole gravity surveys often use the microGal. The acceleration due to gravity, or just "gravity," is measured with a gravimeter. | + | Because the Earth is a rotating oblate spheroid, the quantity g at mean sea level varies with latitude, and it must be corrected for tidal effects. The unit of g is the Gal [1 cm/s<sup>2</sup>]. Surface gravity surveys use the milliGal as the preferred unit. Borehole gravity surveys often use the microGal. The acceleration due to gravity, or just "gravity," is measured with a gravimeter. |

− | ==Measurement== | + | == Measurement == |

− | A borehole gravimeter follows the same principles of operation as a surface gravity meter. It is essentially a very sensitive spring balance. The weight of a horizontal hinged beam with a small mass attached to its free end is balanced by a combination of the tension in a compensating spring and an electrostatic force ('''Fig.1'''). When the acceleration due to gravity increases, the weight of the beam increases, and there is a greater tension in the spring. The spring tension is directly related to the acceleration due to gravity. It is controlled by an adjusting screw for which the number of turns is calibrated in gravitational units. The electrostatic force brings the beam to a horizontal position for reading purposes. It, too, is calibrated in gravitational units. The gravity reading is the difference between the spring tension and the electrostatic force: a tidal gravity correction has to be applied. In this way, differences in gravity can be measured between two places (e.g., between two depth locations in a borehole). Downhole measurements are made by occupying selected gravity stations. They are not continuous measurements with depth. | + | A borehole gravimeter follows the same principles of operation as a surface gravity meter. It is essentially a very sensitive spring balance. The weight of a horizontal hinged beam with a small mass attached to its free end is balanced by a combination of the tension in a compensating spring and an electrostatic force ('''Fig.1'''). When the acceleration due to gravity increases, the weight of the beam increases, and there is a greater tension in the spring. The spring tension is directly related to the acceleration due to gravity. It is controlled by an adjusting screw for which the number of turns is calibrated in gravitational units. The electrostatic force brings the beam to a horizontal position for reading purposes. It, too, is calibrated in gravitational units. The gravity reading is the difference between the spring tension and the electrostatic force: a tidal gravity correction has to be applied. In this way, differences in gravity can be measured between two places (e.g., between two depth locations in a borehole). Downhole measurements are made by occupying selected gravity stations. They are not continuous measurements with depth. |

<gallery widths="300px" heights="200px"> | <gallery widths="300px" heights="200px"> | ||

Line 27: | Line 28: | ||

</gallery> | </gallery> | ||

− | It can be shown<ref name="r3" /> that the difference in gravity, Δ''g'' (mGal), between two locations at the top and bottom of an infinite horizontal reservoir layer penetrated by a vertical well is related to the density, ρ (g/cm<sup>3</sup>), and thickness, ''h''<sub>''m''</sub> (m), of that layer by the expression: | + | It can be shown<ref name="r3">Alixant, J.-L. and Mann, E. 1995. In-Situ Residual Oil Saturation to Gas from Time-Lapse Borehole Gravity. Presented at the SPE Annual Technical Conference and Exhibition, Dallas, Texas, 22-25 October 1995. SPE-30609-MS. http://dx.doi.org/10.2118/30609-MS.</ref> that the difference in gravity, Δ''g'' (mGal), between two locations at the top and bottom of an infinite horizontal reservoir layer penetrated by a vertical well is related to the density, ρ (g/cm<sup>3</sup>), and thickness, ''h''<sub>''m''</sub> (m), of that layer by the expression: |

− | [[File:Vol5 page 0409 eq 004.png]]....................(4) | + | [[File:Vol5 page 0409 eq 004.png|RTENOTITLE]]....................(4) |

where ''F'' is the vertical gradient of gravity (mGal/m), and ''G'' is in cgs units. '''Eq.4''' can be solved for the layer density so that | where ''F'' is the vertical gradient of gravity (mGal/m), and ''G'' is in cgs units. '''Eq.4''' can be solved for the layer density so that | ||

− | [[File:Vol5 page 0409 eq 005.png]]....................(5) | + | [[File:Vol5 page 0409 eq 005.png|RTENOTITLE]]....................(5) |

The gradient, ''F'' (mGal/m), is a function of latitude, ''λ'' (degrees), and elevation, ''h'' (m), as per the Intl. Gravity Formula of 1967, as follows: | The gradient, ''F'' (mGal/m), is a function of latitude, ''λ'' (degrees), and elevation, ''h'' (m), as per the Intl. Gravity Formula of 1967, as follows: | ||

− | [[File:Vol5 page 0410 eq 001.png]]....................(6) | + | [[File:Vol5 page 0410 eq 001.png|RTENOTITLE]]....................(6) |

By substituting for ''F'' in '''Eq.5''', we have | By substituting for ''F'' in '''Eq.5''', we have | ||

− | [[File:Vol5 page 0410 eq 002.png]]....................(7) | + | [[File:Vol5 page 0410 eq 002.png|RTENOTITLE]]....................(7) |

− | where h¯ is the mean elevation of the layer (m), and Δ''g''/''h''<sub>''m''</sub> is in mGal/m. '''Eq.7''' is that most commonly used for deriving density from borehole gravity measurements. If the borehole is deviated at an angle, ''α'', the measured depth interval has to be converted to a true stratigraphic thickness using a specific form of '''Eq.5''' for zero dip. Corrections are needed where the model of an infinite layer breaks down because of the presence of structural discontinuities away from the wellbore. Modern borehole gravimeters can detect gravity differences of a few microGals. | + | where h¯ is the mean elevation of the layer (m), and Δ''g''/''h''<sub>''m''</sub> is in mGal/m. '''Eq.7''' is that most commonly used for deriving density from borehole gravity measurements. If the borehole is deviated at an angle, ''α'', the measured depth interval has to be converted to a true stratigraphic thickness using a specific form of '''Eq.5''' for zero dip. Corrections are needed where the model of an infinite layer breaks down because of the presence of structural discontinuities away from the wellbore. Modern borehole gravimeters can detect gravity differences of a few microGals. |

− | The borehole gravity meter delivers an interval density. It is the only tool that can furnish through-casing density. Where the layer is heterogeneous, the computed density is an average or apparent density. The error in density is a function of the layer thickness. With a LaCoste-Romberg borehole gravimeter, a single measurement of gravity above and below a layer of thickness 6.6 ft [2 m] should result in an error in apparent density of approximately ± 0.025 g/cm<sup>3</sup>. This expected error can be reduced through repeat measurements and by selecting a larger depth interval. Turning this around, the spatial resolution of a borehole gravimeter is governed by the accuracy to which density is required. For example, an accuracy of ±0.01 g/cm<sup>3</sup> would be achieved through three measurements of gravity at the top and at the base of the target layer, provided that the latter is at least 9 ft [2.7 m] thick. | + | The borehole gravity meter delivers an interval density. It is the only tool that can furnish through-casing density. Where the layer is heterogeneous, the computed density is an average or apparent density. The error in density is a function of the layer thickness. With a LaCoste-Romberg borehole gravimeter, a single measurement of gravity above and below a layer of thickness 6.6 ft [2 m] should result in an error in apparent density of approximately ± 0.025 g/cm<sup>3</sup>. This expected error can be reduced through repeat measurements and by selecting a larger depth interval. Turning this around, the spatial resolution of a borehole gravimeter is governed by the accuracy to which density is required. For example, an accuracy of ±0.01 g/cm<sup>3</sup> would be achieved through three measurements of gravity at the top and at the base of the target layer, provided that the latter is at least 9 ft [2.7 m] thick. |

− | Borehole gravity tools have different sizes for different hole conditions. For example, the EDCON tools range in diameter from 3.875 in. [98 mm] for low-temperature (110°C), low-pressure (8,000 psi) applications to 5.25 in. [133 mm] for high-temperature (204°C), high-pressure (20,000 psi) applications. The temperature range can be extended to 260°C with special ring seals. Because of the tool size, there are limits on the deviation of boreholes in which it can be deployed. The measurement stations are located relative to other logging runs by using the gamma log and the [[ | + | Borehole gravity tools have different sizes for different hole conditions. For example, the EDCON tools range in diameter from 3.875 in. [98 mm] for low-temperature (110°C), low-pressure (8,000 psi) applications to 5.25 in. [133 mm] for high-temperature (204°C), high-pressure (20,000 psi) applications. The temperature range can be extended to 260°C with special ring seals. Because of the tool size, there are limits on the deviation of boreholes in which it can be deployed. The measurement stations are located relative to other logging runs by using the gamma log and the [[Casing_collar_locator|casing collar locator]](CCL). The depth of investigation within a homogeneous layer is governed by the contrast between the mud filtrate and formation fluids. It is typically more than 23 ft [7 m]. A larger station spacing, ''h''<sub>''m''</sub>, will not increase this range. It will merely reduce the ability of the tool to see near-well density anomalies. Like surface gravity meters, the tool suffers from drift (of the spring tension), which makes accurate calibration difficult. |

− | ==Application== | + | == Application == |

− | Key thrusts in reservoir evaluation have been the sensing of vuggy, fractured, and heterogeneous reservoirs, in which a deep-sensing porosity measurement is needed to complement the volumetric-sensing capability of deep resistivity logs.<ref name="r4" /><ref name="r5" /> More recent applications have been directed at through-casing monitoring of gas saturations.<ref name="r6" /><ref name="r7" /> In this respect, it is noteworthy that the larger volumes sensed by the borehole gravity meter, relative to conventional density logs, are more closely associated with the simulator grid scale. Moreover, time-lapse gravity measurements are not degraded by structural anomalies. | + | Key thrusts in reservoir evaluation have been the sensing of vuggy, fractured, and heterogeneous reservoirs, in which a deep-sensing porosity measurement is needed to complement the volumetric-sensing capability of deep resistivity logs.<ref name="r4">Jageler, A.H. 1976. Improved Hydrocarbon Reservoir Evaluation Through Use of Borehole-Gravimeter Data. J Pet Technol 28 (6): 709-718. SPE-5511-PA. http://dx.doi.org/10.2118/5511-PA.</ref><ref name="r5">Gournay, L.S. and Lyle, W.D. 1984. Determination of Hydrocarbon Saturation and Porosity Using a Combination Borehole Gravimeter (BHGM) and Deep Investigating Electric Log. Presented at the SPWLA 25th Annual Logging Symposium, 1984. SPWLA-1984-WW.</ref> More recent applications have been directed at through-casing monitoring of gas saturations.<ref name="r6">Popta, J.V., Heywood, J.M.T., Adams, S.J. et al. 1990. Use of Borehole Gravimetry for Reservoir Characterisation and Fluid Saturation Monitoring. Presented at the European Petroleum Conference, The Hague, Netherlands, 21-24 October 1990. SPE-20896-MS. http://dx.doi.org/10.2118/20896-MS.</ref><ref name="r7">Brady, J.L., Wolcott, D.S., and Aiken, C.L.V. 1993. Gravity Methods: Useful Techniques for Reservoir Surveillanc. Presented at the SPE Western Regional Meeting, Anchorage, Alaska, 26-28 May 1993. SPE-26095-MS. http://dx.doi.org/10.2118/26095-MS.</ref> In this respect, it is noteworthy that the larger volumes sensed by the borehole gravity meter, relative to conventional density logs, are more closely associated with the simulator grid scale. Moreover, time-lapse gravity measurements are not degraded by structural anomalies. |

− | As an example, time-lapse borehole gravimetry has been used to determine the residual oil saturation to gas within the oil rim of the onshore Rabi field in Gabon.<ref name="r3" /> The reservoir comprises clean, coarse-grained sands with high-salinity formation water. The required accuracy for residual oil saturation was ± 10 saturation units. A baseline gravity survey was run over an undepleted oil-bearing interval near the gas/oil contact (GOC) (Run 1). As the reservoir is depleted, the GOC moves down and the oil saturation decreases toward its residual value. A second gravity survey (Run 2) allowed the change in gas saturation, Δ''S''<sub>''g''</sub>, to be calculated from the change in measured density, Δ''ρ''<sub>''b''</sub>, porosity, ''ϕ'', and the densities of oil, ''ρ''<sub>''o''</sub>, and gas, ''ρ''<sub>''g''</sub>. | + | As an example, time-lapse borehole gravimetry has been used to determine the residual oil saturation to gas within the oil rim of the onshore Rabi field in Gabon.<ref name="r3">Alixant, J.-L. and Mann, E. 1995. In-Situ Residual Oil Saturation to Gas from Time-Lapse Borehole Gravity. Presented at the SPE Annual Technical Conference and Exhibition, Dallas, Texas, 22-25 October 1995. SPE-30609-MS. http://dx.doi.org/10.2118/30609-MS.</ref> The reservoir comprises clean, coarse-grained sands with high-salinity formation water. The required accuracy for residual oil saturation was ± 10 saturation units. A baseline gravity survey was run over an undepleted oil-bearing interval near the gas/oil contact (GOC) (Run 1). As the reservoir is depleted, the GOC moves down and the oil saturation decreases toward its residual value. A second gravity survey (Run 2) allowed the change in gas saturation, Δ''S''<sub>''g''</sub>, to be calculated from the change in measured density, Δ''ρ''<sub>''b''</sub>, porosity, ''ϕ'', and the densities of oil, ''ρ''<sub>''o''</sub>, and gas, ''ρ''<sub>''g''</sub>. |

− | [[File:Vol5 page 0411 eq 001.png]]....................(8) | + | [[File:Vol5 page 0411 eq 001.png|RTENOTITLE]]....................(8) |

− | Once Δ''S''<sub>''g''</sub> was known, the oil saturation could be calculated, assuming no change in connate-water saturation. | + | Once Δ''S''<sub>''g''</sub> was known, the oil saturation could be calculated, assuming no change in connate-water saturation. |

− | Three surveys were undertaken twelve months apart ('''Fig.2'''). All measurements were made in a data-dedicated borehole. Gravity was measured four to six times at each station, with station intervals as low as 3.3 ft [1 m]. Stations were reoccupied with a shuttle-based system for enhanced depth control. The overall accuracy of the density difference in '''Eq.8''' was 0.015 g/cm<sup>3</sup>. This accuracy corresponds to an accuracy of 0.7 μGal on the station-specific readings and an accuracy of 1.0 μGal on the gravity difference. The residual oil saturation was determined as 15±10 saturation units. Of this uncertainty, eight saturation units could be ascribed to the borehole-gravity measurements and two saturation units to uncertainties in porosity and connate-water saturation. This study set new objectives and standards for borehole gravimetry. | + | Three surveys were undertaken twelve months apart ('''Fig.2'''). All measurements were made in a data-dedicated borehole. Gravity was measured four to six times at each station, with station intervals as low as 3.3 ft [1 m]. Stations were reoccupied with a shuttle-based system for enhanced depth control. The overall accuracy of the density difference in '''Eq.8''' was 0.015 g/cm<sup>3</sup>. This accuracy corresponds to an accuracy of 0.7 μGal on the station-specific readings and an accuracy of 1.0 μGal on the gravity difference. The residual oil saturation was determined as 15±10 saturation units. Of this uncertainty, eight saturation units could be ascribed to the borehole-gravity measurements and two saturation units to uncertainties in porosity and connate-water saturation. This study set new objectives and standards for borehole gravimetry. |

<gallery widths="300px" heights="200px"> | <gallery widths="300px" heights="200px"> | ||

Line 66: | Line 67: | ||

== Nomenclature == | == Nomenclature == | ||

+ | |||

{| | {| | ||

− | |||

− | |||

− | |||

|- | |- | ||

− | |''d'' | + | | ''c'' |

− | |= | + | | = |

− | |gravimeter design constant, L, m | + | | gravimeter design constant, L, m |

− | |- | + | |- |

− | |''f'' | + | | ''d'' |

− | |= | + | | = |

− | |gravitational force of attraction, mL/t<sup>2</sup>, N | + | | gravimeter design constant, L, m |

+ | |- | ||

+ | | ''f'' | ||

+ | | = | ||

+ | | gravitational force of attraction, mL/t<sup>2</sup>, N | ||

|- | |- | ||

− | |''F''<sub>''r''</sub> | + | | ''F''<sub>''r''</sub> |

− | |= | + | | = |

− | |ratio factor for the minimum curvature method | + | | ratio factor for the minimum curvature method |

|- | |- | ||

− | |''F'' | + | | ''F'' |

− | |= | + | | = |

− | |vertical gradient of gravity, 1/t<sup>2</sup>, mGal/m | + | | vertical gradient of gravity, 1/t<sup>2</sup>, mGal/m |

|- | |- | ||

− | |''g'' | + | | ''g'' |

− | |= | + | | = |

− | |acceleration due to gravity, L/t<sup>2</sup>, mGal | + | | acceleration due to gravity, L/t<sup>2</sup>, mGal |

|- | |- | ||

− | |''G'' | + | | ''G'' |

− | |= | + | | = |

− | |universal gravitational constant, L<sup>3</sup>/mt<sup>2</sup>, dyne cm<sup>2</sup>/g<sup>2</sup> | + | | universal gravitational constant, L<sup>3</sup>/mt<sup>2</sup>, dyne cm<sup>2</sup>/g<sup>2</sup> |

|- | |- | ||

− | |''h'' | + | | ''h'' |

− | |= | + | | = |

− | |elevation, L, m | + | | elevation, L, m |

|- | |- | ||

− | |''h''<sub>''m''</sub> | + | | ''h''<sub>''m''</sub> |

− | |= | + | | = |

− | |measured bed thickness, L, m | + | | measured bed thickness, L, m |

|- | |- | ||

− | |''m''<sub>''i''</sub>, ''m''<sub>''j''</sub> | + | | ''m''<sub>''i''</sub>, ''m''<sub>''j''</sub> |

− | |= | + | | = |

− | |any two gravitationally attracting masses, m, kg | + | | any two gravitationally attracting masses, m, kg |

|- | |- | ||

− | |''m'' | + | | ''m'' |

− | |= | + | | = |

− | |mass, m, kg | + | | mass, m, kg |

|- | |- | ||

− | |''M'' | + | | ''M'' |

− | |= | + | | = |

− | |mass of the Earth, m, kg | + | | mass of the Earth, m, kg |

|- | |- | ||

− | |''r'' | + | | ''r'' |

− | |= | + | | = |

− | |distance between two gravitationally attracting masses, L, m | + | | distance between two gravitationally attracting masses, L, m |

|- | |- | ||

− | |''R'' | + | | ''R'' |

− | |= | + | | = |

− | |Earth’s radius, L, km | + | | Earth’s radius, L, km |

|- | |- | ||

− | |[[File:Vol5 page 0410 inline 001.png]] | + | | [[File:Vol5 page 0410 inline 001.png|RTENOTITLE]] |

− | |= | + | | = |

− | |mean elevation of layer, L, m | + | | mean elevation of layer, L, m |

|- | |- | ||

− | |''T'' | + | | ''T'' |

− | |= | + | | = |

− | |tension, mL/t<sup>2</sup>, N | + | | tension, mL/t<sup>2</sup>, N |

|- | |- | ||

− | |''φ'' | + | | ''φ'' |

− | |= | + | | = |

− | |dip, –, degrees | + | | dip, –, degrees |

|- | |- | ||

− | |Δ''g'' | + | | Δ''g'' |

− | |= | + | | = |

− | |difference in gravity between two locations at the top and bottom of an infinite horizontal reservoir layer penetrated by a vertical well, L/t<sup>2</sup>, Gal | + | | difference in gravity between two locations at the top and bottom of an infinite horizontal reservoir layer penetrated by a vertical well, L/t<sup>2</sup>, Gal |

|- | |- | ||

− | |Δ''ρ''<sub>''b''</sub> | + | | Δ''ρ''<sub>''b''</sub> |

− | |= | + | | = |

− | |change in measured density, m/L<sup>3</sup>, g/cm<sup>3</sup> | + | | change in measured density, m/L<sup>3</sup>, g/cm<sup>3</sup> |

|- | |- | ||

− | |''λ'' | + | | ''λ'' |

− | |= | + | | = |

− | |latitude, –, degrees | + | | latitude, –, degrees |

|- | |- | ||

− | |''ρ'' | + | | ''ρ'' |

− | |= | + | | = |

− | |density, m/L<sup>3</sup>, g/cm<sup>3</sup> | + | | density, m/L<sup>3</sup>, g/cm<sup>3</sup> |

|- | |- | ||

− | |''ρ''<sub>''g''</sub> | + | | ''ρ''<sub>''g''</sub> |

− | |= | + | | = |

− | |gas density, m/L<sup>3</sup>, g/cm<sup>3</sup> | + | | gas density, m/L<sup>3</sup>, g/cm<sup>3</sup> |

|- | |- | ||

− | |''ρ''<sub>''o''</sub> | + | | ''ρ''<sub>''o''</sub> |

− | |= | + | | = |

− | |oil density, m/L<sup>3</sup>, g/cm<sup>3</sup> | + | | oil density, m/L<sup>3</sup>, g/cm<sup>3</sup> |

|- | |- | ||

− | |''ϕ'' | + | | ''ϕ'' |

− | |= | + | | = |

− | |porosity | + | | porosity |

|- | |- | ||

− | |''χ'' | + | | ''χ'' |

− | |= | + | | = |

− | |magnetic susceptibility of the rock | + | | magnetic susceptibility of the rock |

|} | |} | ||

− | ==References== | + | == References == |

− | <references | + | <references /> |

− | |||

− | + | == Noteworthy papers in OnePetro == | |

− | + | Use this section to list papers in OnePetro that a reader who wants to learn more should definitely read | |

− | + | == External links == | |

− | + | Use this section to provide links to relevant material on websites other than PetroWiki and OnePetro | |

− | |||

− | |||

− | |||

− | |||

− | |||

− | |||

− | |||

− | |||

− | Use this section to | ||

− | == | + | == See also == |

− | |||

− | + | [[Types_of_logs|Types of logs]] | |

− | [[ | + | [[Borehole_imaging|Borehole imaging]] |

− | [[ | + | [[PEH:Specialized_Well-Logging_Topics]] |

− | [[ | + | [[Category:5.6 Formation evaluation and management]] |

## Latest revision as of 13:28, 12 June 2015

Borehole gravity was pioneered by Smith^{[1]} and then applied to problems of reservoir evaluation by McCulloh *et al.*^{[2]} The borehole gravity meter or gravimeter responds to variations in density. Modern instruments sense a rock volume that is approximately the same as that investigated by deep resistivity tools. Unlike the shallower-sensing density log, the borehole gravimeter is insensitive to wellbore conditions such as rugosity and the presence of casing. Its principal applications are:

- Through-casing time-lapse monitoring of saturations/fluid contacts in gas reservoirs
- Downhole calibration of surface geophysical mapping of geological structures

## Theory

Any two masses, *m*_{i} and *m*_{j}, separated by a distance, *r*, experience a gravitational force of attraction, *f*, which is expressed as

where *G* is the universal gravitational constant (6.6726 × 10 –8 cgs units). More specifically, a mass, *m*, on the surface of the Earth would experience a gravitational force given by

where *M* is the mass of the Earth, *R* is its radius, and an acceleration due to gravity, *g*, is given by

Because the Earth is a rotating oblate spheroid, the quantity g at mean sea level varies with latitude, and it must be corrected for tidal effects. The unit of g is the Gal [1 cm/s^{2}]. Surface gravity surveys use the milliGal as the preferred unit. Borehole gravity surveys often use the microGal. The acceleration due to gravity, or just "gravity," is measured with a gravimeter.

## Measurement

A borehole gravimeter follows the same principles of operation as a surface gravity meter. It is essentially a very sensitive spring balance. The weight of a horizontal hinged beam with a small mass attached to its free end is balanced by a combination of the tension in a compensating spring and an electrostatic force (**Fig.1**). When the acceleration due to gravity increases, the weight of the beam increases, and there is a greater tension in the spring. The spring tension is directly related to the acceleration due to gravity. It is controlled by an adjusting screw for which the number of turns is calibrated in gravitational units. The electrostatic force brings the beam to a horizontal position for reading purposes. It, too, is calibrated in gravitational units. The gravity reading is the difference between the spring tension and the electrostatic force: a tidal gravity correction has to be applied. In this way, differences in gravity can be measured between two places (e.g., between two depth locations in a borehole). Downhole measurements are made by occupying selected gravity stations. They are not continuous measurements with depth.

It can be shown^{[3]} that the difference in gravity, Δ*g* (mGal), between two locations at the top and bottom of an infinite horizontal reservoir layer penetrated by a vertical well is related to the density, ρ (g/cm^{3}), and thickness, *h*_{m} (m), of that layer by the expression:

where *F* is the vertical gradient of gravity (mGal/m), and *G* is in cgs units. **Eq.4** can be solved for the layer density so that

The gradient, *F* (mGal/m), is a function of latitude, *λ* (degrees), and elevation, *h* (m), as per the Intl. Gravity Formula of 1967, as follows:

By substituting for *F* in **Eq.5**, we have

where h¯ is the mean elevation of the layer (m), and Δ*g*/*h*_{m} is in mGal/m. **Eq.7** is that most commonly used for deriving density from borehole gravity measurements. If the borehole is deviated at an angle, *α*, the measured depth interval has to be converted to a true stratigraphic thickness using a specific form of **Eq.5** for zero dip. Corrections are needed where the model of an infinite layer breaks down because of the presence of structural discontinuities away from the wellbore. Modern borehole gravimeters can detect gravity differences of a few microGals.

The borehole gravity meter delivers an interval density. It is the only tool that can furnish through-casing density. Where the layer is heterogeneous, the computed density is an average or apparent density. The error in density is a function of the layer thickness. With a LaCoste-Romberg borehole gravimeter, a single measurement of gravity above and below a layer of thickness 6.6 ft [2 m] should result in an error in apparent density of approximately ± 0.025 g/cm^{3}. This expected error can be reduced through repeat measurements and by selecting a larger depth interval. Turning this around, the spatial resolution of a borehole gravimeter is governed by the accuracy to which density is required. For example, an accuracy of ±0.01 g/cm^{3} would be achieved through three measurements of gravity at the top and at the base of the target layer, provided that the latter is at least 9 ft [2.7 m] thick.

Borehole gravity tools have different sizes for different hole conditions. For example, the EDCON tools range in diameter from 3.875 in. [98 mm] for low-temperature (110°C), low-pressure (8,000 psi) applications to 5.25 in. [133 mm] for high-temperature (204°C), high-pressure (20,000 psi) applications. The temperature range can be extended to 260°C with special ring seals. Because of the tool size, there are limits on the deviation of boreholes in which it can be deployed. The measurement stations are located relative to other logging runs by using the gamma log and the casing collar locator(CCL). The depth of investigation within a homogeneous layer is governed by the contrast between the mud filtrate and formation fluids. It is typically more than 23 ft [7 m]. A larger station spacing, *h*_{m}, will not increase this range. It will merely reduce the ability of the tool to see near-well density anomalies. Like surface gravity meters, the tool suffers from drift (of the spring tension), which makes accurate calibration difficult.

## Application

Key thrusts in reservoir evaluation have been the sensing of vuggy, fractured, and heterogeneous reservoirs, in which a deep-sensing porosity measurement is needed to complement the volumetric-sensing capability of deep resistivity logs.^{[4]}^{[5]} More recent applications have been directed at through-casing monitoring of gas saturations.^{[6]}^{[7]} In this respect, it is noteworthy that the larger volumes sensed by the borehole gravity meter, relative to conventional density logs, are more closely associated with the simulator grid scale. Moreover, time-lapse gravity measurements are not degraded by structural anomalies.

As an example, time-lapse borehole gravimetry has been used to determine the residual oil saturation to gas within the oil rim of the onshore Rabi field in Gabon.^{[3]} The reservoir comprises clean, coarse-grained sands with high-salinity formation water. The required accuracy for residual oil saturation was ± 10 saturation units. A baseline gravity survey was run over an undepleted oil-bearing interval near the gas/oil contact (GOC) (Run 1). As the reservoir is depleted, the GOC moves down and the oil saturation decreases toward its residual value. A second gravity survey (Run 2) allowed the change in gas saturation, Δ*S*_{g}, to be calculated from the change in measured density, Δ*ρ*_{b}, porosity, *ϕ*, and the densities of oil, *ρ*_{o}, and gas, *ρ*_{g}.

Once Δ*S*_{g} was known, the oil saturation could be calculated, assuming no change in connate-water saturation.

Three surveys were undertaken twelve months apart (**Fig.2**). All measurements were made in a data-dedicated borehole. Gravity was measured four to six times at each station, with station intervals as low as 3.3 ft [1 m]. Stations were reoccupied with a shuttle-based system for enhanced depth control. The overall accuracy of the density difference in **Eq.8** was 0.015 g/cm^{3}. This accuracy corresponds to an accuracy of 0.7 μGal on the station-specific readings and an accuracy of 1.0 μGal on the gravity difference. The residual oil saturation was determined as 15±10 saturation units. Of this uncertainty, eight saturation units could be ascribed to the borehole-gravity measurements and two saturation units to uncertainties in porosity and connate-water saturation. This study set new objectives and standards for borehole gravimetry.

**Fig.2 – Time-laps borehole gravimeter data around the GOC in the Rabi field, Gabon. Note the excellent repeatability of the 1995 shuttle data at approximately 1099 m measured depth. The data recorded between 1098.5 and 1101 m measured det were used to assess the remaining oil saturation to gas.**^{[3]}

## Nomenclature

## References

- ↑ Smith, N.L. 1950. The Case for Gravity Data From Boreholes. Geophysics 15 (4): 605-636. http://dx.doi.org/ 10.1190/1.1437623.
- ↑ McCulloh, T.H., Kandle, G.R., and Schoellhamer, J.E. 1968. Application of Gravity Measurements in Wells to Problems of Reservoir Evaluation. Trans., Soc. of Professional Well Log Analysts 9th Annual Logging Symposium, New Orleans, paper O.
- ↑
^{3.0}^{3.1}^{3.2}Alixant, J.-L. and Mann, E. 1995. In-Situ Residual Oil Saturation to Gas from Time-Lapse Borehole Gravity. Presented at the SPE Annual Technical Conference and Exhibition, Dallas, Texas, 22-25 October 1995. SPE-30609-MS. http://dx.doi.org/10.2118/30609-MS. - ↑ Jageler, A.H. 1976. Improved Hydrocarbon Reservoir Evaluation Through Use of Borehole-Gravimeter Data. J Pet Technol 28 (6): 709-718. SPE-5511-PA. http://dx.doi.org/10.2118/5511-PA.
- ↑ Gournay, L.S. and Lyle, W.D. 1984. Determination of Hydrocarbon Saturation and Porosity Using a Combination Borehole Gravimeter (BHGM) and Deep Investigating Electric Log. Presented at the SPWLA 25th Annual Logging Symposium, 1984. SPWLA-1984-WW.
- ↑ Popta, J.V., Heywood, J.M.T., Adams, S.J. et al. 1990. Use of Borehole Gravimetry for Reservoir Characterisation and Fluid Saturation Monitoring. Presented at the European Petroleum Conference, The Hague, Netherlands, 21-24 October 1990. SPE-20896-MS. http://dx.doi.org/10.2118/20896-MS.
- ↑ Brady, J.L., Wolcott, D.S., and Aiken, C.L.V. 1993. Gravity Methods: Useful Techniques for Reservoir Surveillanc. Presented at the SPE Western Regional Meeting, Anchorage, Alaska, 26-28 May 1993. SPE-26095-MS. http://dx.doi.org/10.2118/26095-MS.

## Noteworthy papers in OnePetro

Use this section to list papers in OnePetro that a reader who wants to learn more should definitely read

## External links

Use this section to provide links to relevant material on websites other than PetroWiki and OnePetro