Aspects of the disclosure relate to gamma-gamma measurements used to characterize properties of downhole structures. More specifically, aspects of the disclosure relate to methods and apparatus that may be used to determine casing thickness, collar presence, slot presence, packer integrity, density of the material in the annulus and other features used in recovery of hydrocarbons from geological stratum.
Gamma-gamma measurements may be used in cased-hole environment to estimate formation density. In order to achieve correct results, completion properties should be known and a backscatter detector is commonly used to assess apparent casing thickness. Casing collars are assessed during such investigations. Conventionally, short and long-spacing detectors are then used in combination to estimate formation density, assuming the cement density is known and constant. Cement evaluation through Gamma-Gamma measurements has also been independently addressed by conventional technologies.
The estimation of borehole diameter via a single long spacing detector , and directly using correlation between counting rates and hole diameter for controlled annulus and formation properties are also conventionally known.
While the conventional methods and apparatus described above are helpful in determining specific instances and features, there is a need to provide an apparatus and method that are easier to operate than conventional apparatus and methods.
There is a further need to provide apparatus and methods that do not have the drawbacks of these conventional systems, such as being cumbersome to calculate desired features, while providing a more accurate methodology to achieve better results than are possible with conventional technologies.
There is a still further need to reduce economic costs associated with operations and apparatus described above with conventional tools such that operators may more efficiently characterize needed information.
So that the manner in which the above recited features of the present disclosure can be understood in detail, a more particular description of the disclosure, briefly summarized below, may be had by reference to embodiments, some of which are illustrated in the drawings. It is to be noted that the drawings illustrate only typical embodiments of this disclosure and are therefore not to be considered limiting of its scope, for the disclosure may admit to other equally effective embodiments without specific recitation. Accordingly, the following summary provides just a few aspects of the description and should not be used to limit the described embodiments to a single concept.
In one example embodiment, a method of characterizing at least one downhole feature is disclosed. The method may comprise using a gamma-gamma tool, performing at least one scan to characterize casing properties and an annulus density indicator and correcting apparent densities from casing and annulus density effect.
In another example embodiment, a method is disclosed. The method may comprise producing at least one pulse of energy from a downhole tool, the at least one pulse of energy directed into a casing and annulus of a wellbore. The method may further comprise receiving the at least one pulse of energy at the downhole tool. The method may also comprise characterizing casing properties and an annulus density indicator from the received pulse of energy. The method may also comprise correcting apparent densities from casing and annulus density effect.
In another example embodiment, an apparatus is disclosed. The apparatus may be comprised of a downhole tool configured to transmit and receive gamma radiation into a casing and surrounding annulus of a wellbore, the downhole tool having at least three detectors to receive the gamma radiation. The apparatus may further be comprised of a computing apparatus configured to receive data from the downhole tool, wherein the computing apparatus is configured to perform a correction of data received based upon a combination of a first shallow measurement and a second deeper measurement.
So that the manner in which the above recited features of the present disclosure can be understood in detail, a more particular description of the disclosure, briefly summarized above, may be had by reference to embodiments, some of which are illustrated in the drawings. It is to be noted, however, that the appended drawings illustrate only typical embodiments of this disclosure and are therefore not to be considered limiting of its scope, for the disclosure may admit to other equally effective embodiments.
To facilitate understanding, identical reference numerals have been used, where possible, to designate identical elements that are common to the figures (“FIGS”). It is contemplated that elements disclosed in one embodiment may be beneficially utilized on other embodiments without specific recitation.
In the following, reference is made to embodiments of the disclosure. It should be understood, however, that the disclosure is not limited to specific described embodiments. Instead, any combination of the following features and elements, whether related to different embodiments or not, is contemplated to implement and practice the disclosure. Furthermore, although embodiments of the disclosure may achieve advantages over other possible solutions and/or over the prior art, whether or not a particular advantage is achieved by a given embodiment is not limiting of the disclosure. Thus, the following aspects, features, embodiments and advantages are merely illustrative and are not considered elements or limitations of the claims except where explicitly recited in a claim. Likewise, reference to “the disclosure” shall not be construed as a generalization of inventive subject matter disclosed herein and should not be considered to be an element or limitation of the claims except where explicitly recited in a claim.
Although the terms first, second, third, etc., may be used herein to describe various elements, components, regions, layers and/or sections, these elements, components, regions, layers and/or sections should not be limited by these terms. These terms may be only used to distinguish one element, component, region, layer or section from another region, layer or section. Terms such as “first”, “second” and other numerical terms, when used herein, do not imply a sequence or order unless clearly indicated by the context. Thus, a first element, component, region, layer or section discussed herein could be termed a second element, component, region, layer or section without departing from the teachings of the example embodiments.
When an element or layer is referred to as being “on,” “engaged to,” “connected to,” or “coupled to” another element or layer, it may be directly on, engaged, connected, coupled to the other element or layer, or interleaving elements or layers may be present. In contrast, when an element is referred to as being “directly on,” “directly engaged to,” “directly connected to,” or “directly coupled to” another element or layer, there may be no interleaving elements or layers present. Other words used to describe the relationship between elements should be interpreted in a like fashion. As used herein, the term “and/or” includes any and all combinations of one or more of the associated listed terms.
Some embodiments will now be described with reference to the figures. Like elements in the various figures will be referenced with like numbers for consistency. In the following description, numerous details are set forth to provide an understanding of various embodiments and/or features. It will be understood, however, by those skilled in the art, that some embodiments may be practiced without many of these details, and that numerous variations or modifications from the described embodiments are possible. As used herein, the terms “above” and “below”, “up” and “down”, “upper” and “lower”, “upwardly” and “downwardly”, and other like terms indicating relative positions above or below a given point are used in this description to more clearly describe certain embodiments.
Aspects of the disclosure relate to apparatus and methods for characterizing liner and annulus properties. Such characterization may be done through use of a gamma-gamma tool. In such tools, the transport of gamma ray radiation from a source to a detector in an energy range below 1 MeV is dominated by two types of interactions with the atoms: the scattering off the atom electrons and the photoelectric absorption. The cross sections of both interactions depend on the incident gamma ray energy, weakly for the scattering effect but strongly for the photoelectric effect. This latter effect dominates at low energy, below some 100 s keV or so for common atoms.
The gamma ray reaching a detector is stacked into an energy spectrum. Contiguous bins can also be grouped into so-called energy windows. The gamma-gamma tool used in aspects of the current disclosure may possess any number of detectors. In example embodiments provided herein, three (3) detectors are used. For clarity, the three detectors used are defined as Backscatter (BS), Short Spacing (SS), and Long Spacing (LS), and wherein 3 to 4 energy windows are produced per detector.
There is no exact formulation for the dependence of each energy window to formation electronic density and photoelectric factor. Common models use simple phenomenological formulations that combine an attenuation term from bulk gamma ray flux attenuation that decreases exponentially with the density and Pe and a scattering term proportional to electronic density. The parameters of these models depend on the detector and on the selected energy windows. Detectors that are located close to the source will have a larger contribution from the scattering term, as gamma rays need to scatter back to the tool while long-spaced detectors will have a dominant attenuation term.
In terms of depth of investigation, in addition to the effect of source-to-detector spacing, there is an additional effect linked to the energy window. High energy windows tend to have a shallower depth of investigation (“DOI”) than low energy windows, as high energy gamma ray reaching the detector can only have endured a limited number of scattering and hence could not travel far.
Because the photoelectric factor has a depressing action on low energy windows, removing part of the gamma rays that could have been detected, the DOI of the photoelectric factor is somehow shallower than that of density, and does not depend on spacing in the “attenuation” transport regime.
This difference in depth of investigation between energy windows is implicitly accounted for in processing workflows based on algorithmic inversion of a forward model, when the forward model is defined on a database with radial information (for example mudcakes of different thicknesses).
An example of difference in DOI between different detectors and energy windows is provided in
Aspects of the current disclosure characterize a cased-hole completion using the difference in radial sensitivity of different energy windows and detectors. Aspects of the disclosure provide the following new features:
BS energy windows show different DOI. The 50%-DOI is typically within the 0.3-0.7-in range (for OH-measurements). When casing is present, the shallower measurements are mostly sensitive to the casing, with some residual sensitivity to annulus properties. The measurement with deeper DOI will be more sensitive to the annulus and possibly formation properties and by combining the different energy windows, we can extract information related to casing and annulus properties.
For that purpose, an apparent casing thickness for each energy window of the BS detector is estimated.
The red lines correspond to an example of apparent thickness estimator. It is clear that depending on the true casing thickness, the apparent casing thickness can be strongly overestimated for the low energy windows, with a sensitivity to both annulus and formation properties. The general use of this type of direct estimator can be misleading, but in controlled conditions, with little variation around a “nominal” point, they bring forward useful pieces of information.
In embodiments, a fit-for-conditions solution is used, optimizing the workflow around expected operational conditions. A different objective may be targeted, which is the characterization of completion (casing and annulus), and not characterization of the formation itself.
In the following, it is assumed, for illustration, that the workflow is optimized for a 7-in casing or line, with possibly varying annulus densities, inside a 8.75-in borehole. It is assumed that the open hole formation density is known, or at least have some correlated indicator, that is defined as HDEN.
The first step is to remove the contribution from the formation to the apparent casing thickness estimated described in
Once the residual formation density effects are removed from the counting rates, a casing estimator can be computed. They are now functions of true casing thickness and annulus density. An example of equation that can be used to correlate the counting rates in a given energy window of a backscatter measurement is as follows:
CSTKBSW=α log(β log(BSW)+γ) EQ1
The projection onto the casing absorption axis following the annulus density axis constitutes a measurement called casing absorption. It is zeroed at expected nominal casing absorption point to provide an index that varies around zero.
The projection onto the annulus density axis following the casing absorption axis constitutes a measurement called annulus density, which is scaled to represent density variations from a gas density to a high cement density.
The same process is applied to short spacing detector to obtain two complementary estimations that may be used for quality control, as short spacing is less sensitive to casing and annulus properties than back scatter.
Resulting indicators of casing absorption and annulus density are shown in
A classification flag, see 210, is created based on casing absorption and annulus density indicators. It classifies different depth intervals with the following labels (given as illustration, but not restricted to):
Using an open-hole database, and following existing methods, apparent monosensor densities can be derived from SS and LS in all conditions, including cased hole conditions. Of course, in the latter conditions, these apparent densities are affected by the presence of casing and annulus, and not only by formation properties.
Once the casing properties and the annulus density indicator are estimated in step 1, it is possible to correct the apparent SS and LS densities from the effect of casing and from the effect of an annulus of varying density, according to the previously derived annulus density indicator curves, but of “nominal” thickness. This correction can be calibrated on MCNP simulation and on log data where OH reference density is known. Example of linear calibration is given in
For each detector, a combination of apparent single energy window densities is created to optimize the robustness of the answer. This combination results in a single corrected apparent density for each detector, namely RHSA_CH and RHLA_CH.
These CH apparent densities (RHSA_CH or RHLA_CH) should overlap with the OH formation density (DEN_OH), provided the annulus thickness is nominal. Deviations from nominal annulus thickness can be parameterized by a radial response factor J_RHXA defined as:
This expression is similar to that of the radial response function J defined in the introduction, with the exception that the formation density is replaced by fixed 2.5 g/cc in the denominator. This substitution aims at stabilizing JCH in front of low-density coal layers.
The quantification of actual annulus thickness relies on a transform of J_RHXA into thickness excess to which is added the nominal annulus thickness to provide ATHK_RHXA, as illustrated in
ATHK_RHXA=ATHK_NOM+α(ANN_DEN)×atan(J_RHXA), EQ. 3
Where ATHK_NOM is the monila annulus thickness, and α(ANN_DEN) is a function fitted on a MCNP or controlled log database. Resulting logs are illustrated in
In the presence of casing and annulus, the sensitivity of the different energy windows to formation properties is significantly affected. In particular, the shallow BS measurements quickly lose their sensitivity to formation properties.
To capture the sensitivity of the different BS energy windows to completion properties, transforms are defined between counting rates and apparent casing thickness; these are represented by the black lines in
In
This method may be used to separate and quantify the effects of annulus density and casing thickness. The uncertainty associated with this method, however, was considered unacceptably high due to the scarcity of the cased-hole database and the absence of numerical (MCNP®) modeling to support the exact response. In one embodiment, a more robust analysis is presented. By projecting onto the pink, 0.317-in casing line along the water-annulus casing thickness axis, an indicator is obtained of annulus density, referred to as ANNULUS_DEN. In addition, by projecting on the casing thickness axis along the pink line direction, an estimation of casing thickness is obtained, from which we subtract the nominal casing thickness of 0.317 in to end up with an indicator called CASING_ABSORB.
In another embodiment, interpretation of residual variations as a function of annulus thickness and estimation of such is further provided. For that purpose, an approach based on the radial response function formalism as shown in Equation 4.
In this embodiment, a density correction workflow is provided for a three-step preparatory analysis. First, the radial response function behind a 0.3-in casing must be evaluated; second, the effect on this response of varying annulus density must be accounted for; and last, the response must be zeroed around a nominal 0.875-in annulus thickness to follow the corrected density definition.
To address the first step, numerical simulation results are used that were obtained during the derivation of the density cased-hole algorithm. The annulus material is cement, cement thickness ranges from 0 to 2 in, and formation density ranges from 1.7 to 3.03 g/cc. The modeled response functions take the simple form of a hyperbolic tangent with a single coefficient λ, such that
where ha is the annulus thickness, and λ is a function of the casing thickness CSTK. For each casing thickness, a specific transform calibrated for the no-cement case allows accounting for pure casing effect and hence estimating the residual radial response to cement thickness. For energy window 3 from the short and long spacing detectors, we obtain the results listed in Table 1.
At this stage, there is no estimate of the effect of varying annulus density on the radial response function, so we will use the cement-type response and acknowledge that the annulus thickness estimation is not fully quantitative but instead captures its qualitative dynamic range.
Using Equation 5, for the casing-only-corrected density ρcas,
ρcas=ρf+J(ha)(ρa−ρf) (EQ . 6)
The casing and nominal cement thickness corrected density ρCH becomes
where J0=J(ho) with ho the nominal thickness. By using a Taylor expansion around the nominal thickness, the change in annulus thickness radius may be estimated through
For short spacing and long spacing measurements, using the values listed in Table 1 for 0.3 in casing and cement, the geometrical terms have values CSS3=0.52 and CLS3=1.78. While it is understood that these values are indicative only, they still serve as guidance for the actual implementation. The annulus thickness obtained from the long spacing detector is called ATHK_RHLA in the following section.
In front of coal layers, the formation density ρf can be as low as or even lower than the annulus density, and the computation proposed in Equation 7 may become inaccurate. To avoid presenting confusing results in such cases, the computation is limited to noncoal intervals. Similarly, intervals classified as collar or packer are excluded from the computation.
In one example embodiment, a method of characterizing at least one downhole feature is disclosed. The method may comprise using a gamma-gamma tool, performing at least one scan to characterize casing properties and an annulus density indicator and correcting apparent densities from casing and annulus density effect.
In another example embodiment, the method may be performed wherein the correcting apparent densities include correcting an apparent short-spaced density and long-spaced density.
In another example embodiment, the method may be performed wherein the correcting apparent densities further comprise assuming a nominal annulus thickness.
In another example embodiment, the method may be performed wherein the correcting apparent densities include estimating an annulus thickness from combining short spaced and long spaced densities with a OH density.
In another example embodiment, the method may be performed wherein the performing the at least one scan includes creating a classification flag.
In another example embodiment, the method may be performed wherein the creating the classification flag includes at least one of the following: water in the annulus, cement in the annulus, heavy cement in the annulus, presence of a liner slot.
In another example embodiment, the method may be performed wherein the classification flag includes at least one of the following: a casing collar, a casing de-centralizer and a packer.
In another example embodiment, a method is disclosed. The method may comprise producing at least one pulse of energy from a downhole tool, the at least one pulse of energy directed into a casing and annulus of a wellbore. The method may further comprise receiving the at least one pulse of energy at the downhole tool. The method may also comprise characterizing casing properties and an annulus density indicator from the received pulse of energy. The method may also comprise correcting apparent densities from casing and annulus density effect.
In another example embodiment, the method may be performed wherein the correcting apparent densities include correcting apparent short-spaced densities.
In another example embodiment, the method may be performed wherein the correcting apparent densities include correcting apparent long-spaced densities.
In another example embodiment, the method may be performed wherein the downhole tool produces gamma radiation.
In another example embodiment, the method may be performed wherein the downhole tool has three detectors.
In another example embodiment, the method may be performed wherein the characterizing casing properties and an annulus density indicator from the received pulse of energy further comprise creating a classification flag.
In another example embodiment, the method may be performed wherein the creating the classification flag includes at least one of the following: water in the annulus, cement in the annulus, heavy cement in the annulus, presence of a liner slot.
In another example embodiment, the method may be performed wherein the classification flag includes at least one of the following: a casing collar, a casing de-centralizer and a packer.
In another example embodiment, an apparatus is disclosed. The apparatus may be comprised of a downhole tool configured to transmit and receive gamma radiation into a casing and surrounding annulus of a wellbore, the downhole tool having at least three detectors to receive the gamma radiation. The apparatus may further be comprised of a computing apparatus configured to receive data from the downhole tool, wherein the computing apparatus is configured to perform a correction of data received based upon a combination of a first shallow measurement a second deeper measurement.
In another example embodiment, the apparatus may be comprised wherein the downhole tool is configured to transmit energy in an energy range below 1 MeV.
In another example embodiment, the apparatus may be comprised wherein the computing apparatus and the downhole tool are connected through a wire.
The foregoing description of the embodiments has been provided for purposes of illustration and description. It is not intended to be exhaustive or to limit the disclosure. Individual elements or features of a particular embodiment are generally not limited to that particular embodiment, but, where applicable, are interchangeable and can be used in a selected embodiment, even if not specifically shown or described. The same may be varied in many ways. Such variations are not to be regarded as a departure from the disclosure, and all such modifications are intended to be included within the scope of the disclosure.
While embodiments have been described herein, those skilled in the art, having benefit of this disclosure, will appreciate that other embodiments are envisioned that do not depart from the inventive scope. Accordingly, the scope of the present claims or any subsequent claims shall not be unduly limited by the description of the embodiments described herein.
This Application claims priority to US Provisional Patent Application serial no. 63/164205, which was filed on Mar. 22, 2021. The entirety of which is incorporated herein.
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/US2022/020568 | 3/16/2022 | WO |
Number | Date | Country | |
---|---|---|---|
63164205 | Mar 2021 | US |