Login| Sign Up| Help| Contact|

Patent Searching and Data


Title:
BEHIND PIPE EVALUATION OF CUT AND PULL TENSION PREDICTION IN WELL ABANDONMENT AND INTERVENTION OPERATIONS
Document Type and Number:
WIPO Patent Application WO/2016/028537
Kind Code:
A1
Abstract:
Casing removal operations may utilize an analytical model to the estimate the total force to (1) break bonding between a cut casing and the annular materials contained between the cut casing and the surrounding borehole and (2) pull the cut casing upward through the borehole while taking into account frictional forces inside and outside the cut casing, shear forces, and the like. Such models may further predict the casing overpull and total tension needed to extract the cut casing from the borehole. Using these models in conjunction with the specifications of the casing retrieval equipment may allow for determining the cutting depth for producing the cut casing with higher accuracy than previously available.

Inventors:
SUPARMAN FNU (US)
FOX PHILIP EDMUND (US)
FRISCH GARY JAMES (US)
Application Number:
PCT/US2015/044508
Publication Date:
February 25, 2016
Filing Date:
August 10, 2015
Export Citation:
Click for automatic bibliography generation   Help
Assignee:
HALLIBURTON ENERGY SERVICES INC (US)
International Classes:
G05B19/02; E21B29/00; G06F9/455
Foreign References:
US20110048701A12011-03-03
US20120242342A12012-09-27
US4538684A1985-09-03
US20090266544A12009-10-29
US20020060076A12002-05-23
Other References:
See also references of EP 3152386A4
Attorney, Agent or Firm:
KAISER, Iona et al. (500 North Capitol Street N.W, Washington District of Columbia, US)
Download PDF:
Claims:
CLAIMS

The invention claimed is:

1. A method comprising :

measuring at least one characteristic of fluids and solids disposed in an annulus defined between a casing and a wellbore with at least one of: a cement bond log tool, an ultrasonic device, or a neutron sensor, and thereby providing an as-logged annulus;

modeling a liquid-filled annulus based on the as-logged annulus without solids therein;

modeling a plurality of cut casings at different cutting depths;

calculating a total force for the as-logged annulus {Ftot s ,g) for extracting each of the plurality of cut casings;

calculating a total force for the liquid-filled annulus (Ftot q) for extracting each of the plurality of cut casings; and

determining a preferred cutting depth based on the Ftot S log for each of the plurality of cut casings, the Ftot q for each of the plurality of cut casings, and specifications of casing retrieval equipment.

2. The method of claim 1 further comprising :

performing a casing removal operation at the preferred cutting depth with the casing retrieval equipment.

3. The method of claim 1, wherein Ftotias log is calculated according to Equation 1 :

Ptot.as log ~ EWCS + FFin as iog + FF0Ut as iog + FS0Ut as iog ~^ F shear Equation 1 where: EWCS is a total effective casing weight of the cut casing;

PPin,as iog is a total fluid frictional force between a fluid contained within the cut casing and an inner surface of the cut casing in the as- logged annulus for a total length of the cut casing;

PPout.as iog is a total fluid frictional force between a fluid contained within the annulus and an outer surface of the cut casing in the as- logged annulus for the total length of the cut casing;

PSout.as iog is a total solid frictional force between solids contained within the annulus and the outer surface of the cut casing in the as-logged annulus for the total length of the cut casing; P shear is a total shear force to break the bond between a cement sheath in the annulus and the cut casing in the as-logged annulus;

PPin,iiq is a total fluid frictional force between a fluid contained within the cut casing and the inner surface of the cut casing in the liquid- filled annulus for the total length of the cut casing; and

PPout.uq is a total fluid frictional force between a fluid contained within the annulus and the outer surface of the cut casing in the liquid- filled annulus for the total length of the cut casing; and

4. The method of claim 1 or 3, wherein Ftot q is calculated according to Equation 2:

Pfot.liq ~ EWCS + FFinMq + FF0ut,iiq Equation 2 where: EWCS is a total effective casing weight of the cut casing;

PPin.asiog is a total fluid frictional force between a fluid contained within the cut casing and an inner surface of the cut casing in the as- logged annulus for a total length of the cut casing;

PPout.asiog is a total fluid frictional force between a fluid contained within the annulus and an outer surface of the cut casing in the as- logged annulus for the total length of the cut casing;

PSout,asiog is a total solid frictional force between solids contained within the annulus and the outer surface of the cut casing in the as-logged annulus for the total length of the cut casing;

P shear is a total shear force to break the bond between a cement sheath in the annulus and the cut casing in the as-logged annulus;

PPin.uq 'S a total fluid frictional force between a fluid contained within the cut casing and the inner surface of the cut casing in the liquid- filled annulus for the total length of the cut casing; and

PPout.uq 'S a total fluid frictional force between a fluid contained within the annulus and the outer surface of the cut casing in the liquid- filled annulus for the total length of the cut casing; and

5. The method of claim 4 further comprising:

calculating a drag force {Fdrag) according to Equation 4

Pdrag ~ P tot, as log ~ Ptot.liq ~ P^ ' out, as log PS out, as log P shear ~ PPout.liq Equation 4, wherein determining the preferred cutting depth is further based on

6. The method of claim 5, wherein determining the preferred cutting depth comprises: calculating a threshold force {Fthresh) based on the specifications of casing retrieval equipment and the preferred cutting depth provides for Ftot as log < Fthresh and/or Fdrag + Ftot liq < Fthresh.

7. The method of claim 6, wherein determining the preferred cutting depth further comprises: calculating a tension for the as-logged annulus (Ttot,as iog) using Ftot S log i a tension for the liquid-filled annulus {TtotMq ) using Ftot,iiq> a tension for the drag force {Tdrag ) using Fdrag i and a threshold tension (Tthresh) based on the specifications of casing retrieval equipment, wherein the preferred cutting depth provides for Ttot as iog < Tthresh and/or Tdrag + Ttot iq

Tthresh ·

8. The method of claim 6, wherein determining the preferred cutting depth further comprises: calculating an overpull for the as-logged annulus (OPtot,as iog) using Ftot S log i an overpull for the liquid-filled annulus {OPtot liq ) using Ftot,iiq> a n overpull for the drag force {OPdrag) using Fdrag i and a threshold overpull (OPthresh) based on the specifications of casing retrieval equipment, wherein the preferred cutting depth provides for OPtot as iog < OPthresh and/or

OPdrag + OPtot,liq OPthresh -

9. A method comprising :

measuring at least one characteristic of fluids and solids disposed in an annulus formed by a casing and a wellbore with at least one of: a cement bond log, an ultrasonic device, or a neutron sensor, thereby providing an as- logged annulus;

modeling a liquid-filled annulus based on the as-logged annulus without solids therein;

modeling a plurality of cut casings at different cutting depths;

calculating a drag force (Fdrag ) for extracting each of the plurality of cut casings based on the as-logged annulus and the liquid-filled annulus;

determining an preferred cutting depth based on the Fdrag for each of the plurality of cut casings and specifications of casing retrieval equipment.

10. The method of claim 9 further comprising :

performing a casing removal operation at the preferred cutting depth with the casing retrieval equipment.

11. The method of claim 9, wherein Fdrag is calculated according to Equation 4 :

Fdrag ~ FF out.as log FS Out.as log F shear ~ FFout.liq Equation 4 where : FFout as iog is a total flu id frictional force between a fluid contained within the annulus and an outer surface of the cut casing in the as- logged annu lus for the total length of the cut casing;

FSout,as iog is a total solid frictional force between solids contained within the annulus and an inner surface of the cut casing in the as-logged annulus for the total length of the cut casing ;

F shear is a total shear force to break the bond between a cement sheath in the annulus and the cut casing in the as-logged annu lus; and

FF0Utiuq is a total fluid frictional force between a fluid contained within the annu lus and the outer surface of the cut casing in the liqu id- filled annu lus for the total length of the cut casing .

12. The method of claim 9, wherei n determining the preferred cutting depth comprises : calcu lating a threshold force {Fthresh) based on the specifications of casing retrieval equ ipment, wherein the preferred cutting depth provides for Fdrag < Fthresh .

13. The method of claim 9, wherei n determining the preferred cutting depth comprises : calculating a tension {Tdrag) using Fdrag and a threshold tension (Tthresh) based on the specifications of casing retrieval equ ipment, wherein the preferred cutting depth provides for Tdrag < Tthresh .

14. The method of claim 9, wherei n determining the preferred cutting depth comprises : calculating an overpu ll {OPdrag) using Fdrag and a threshold overpull (OPthresh) based on the specifications of casing retrieval equipment, and wherein the preferred cutting depth provides for OPdrag < OPthresh .

15. A system comprising :

a conveyance extendable into a wellbore lined at least partially with a casing, wherein the conveyance is coupled to at least one of: a cement bond log tool, an ultrasonic device, or a neutron sensor, and wherein fluids and solids are contained in an annulus defined between the casing and the wellbore; and a control system that includes a non-transitory mediu m readable by a processor and storing instructions that when executed by the processor cause the control system to : measure at least one characteristic of the fluids and solids the at least one of: the cement bond log tool, the ultrasonic device, or the neutron sensor, and thereby providing an as-logged annulus;

model a liquid-filled annulus based on the as-logged annulus without solids therein;

model a plurality of cut casings at different cutting depths; calculate a total force for the as-logged annulus {Ftot,asiog) for extracting each of the plurality of cut casings;

calculate a total force for the liquid-filled annulus {FtotMq) for extracting each of the plurality of cut casings; and

determine a preferred cutting depth based on the Ftot Slog for each of the plurality of cut casings, the Ftot q for each of the plurality of cut casings, and specifications of casing retrieval equipment.

16. The system of claim 15, wherein Ftot S log is calculated according to Equation 1:

Ptot.as log ~ EWCS + FFin as iog + FF0Ut as iog + FS0Ut as iog ~^ F shear Equation 1 where: EWCS is a total effective casing weight of the cut casing;

PPin,asiog is a total fluid frictional force between a fluid contained within the cut casing and an inner surface of the cut casing in the as- logged annulus for a total length of the cut casing;

PPout.asiog is a total fluid frictional force between a fluid contained within the annulus and an outer surface of the cut casing in the as- logged annulus for the total length of the cut casing;

PSout,asiog is a total solid frictional force between solids contained within the annulus and the outer surface of the cut casing in the as-logged annulus for the total length of the cut casing;

Pshear is a total shear force to break the bond between a cement sheath in the annulus and the cut casing in the as-logged annulus;

PPin.uq 'S a total fluid frictional force between a fluid contained within the cut casing and the inner surface of the cut casing in the liquid- filled annulus for the total length of the cut casing; and

PPout.uq 'S a total fluid frictional force between a fluid contained within the annulus and the outer surface of the cut casing in the liquid- filled annulus for the total length of the cut casing; and

17. The system of claim 15 or 16, FtotMq is calculated according to Equation 2:

Pfot.liq ~ EWCS + FFinMq + FF0UtMq Equation 2 where: EWCS is a total effective casing weight of the cut casing;

PPin,asiog is a total fluid frictional force between a fluid contained within the cut casing and an inner surface of the cut casing in the as- logged annulus for a total length of the cut casing;

PPout.asiog is a total fluid frictional force between a fluid contained within the annulus and an outer surface of the cut casing in the as- logged annulus for the total length of the cut casing;

PSout,asiog is a total solid frictional force between solids contained within the annulus and the outer surface of the cut casing in the as-logged annulus for the total length of the cut casing;

Pshear is a total shear force to break the bond between a cement sheath in the annulus and the cut casing in the as-logged annulus;

PPin.uq 'S a total fluid frictional force between a fluid contained within the cut casing and the inner surface of the cut casing in the liquid- filled annulus for the total length of the cut casing; and

PPout.uq 'S a total fluid frictional force between a fluid contained within the annulus and the outer surface of the cut casing in the liquid- filled annulus for the total length of the cut casing; and

18. A system comprising:

a conveyance extendable into a wellbore lined at least partially with a casing, wherein the conveyance is coupled to at least one of: a cement bond log tool, an ultrasonic device, or a neutron sensor, and wherein fluids and solids are contained in an annulus defined between the casing and the wellbore; and a control system that includes a non-transitory medium readable by a processor and storing instructions that when executed by the processor cause the control system to:

measure at least one characteristic of fluids and solids disposed in an annulus formed by a casing and a wellbore with at least one of: a cement bond log, an ultrasonic device, or a neutron sensor, thereby providing an as-logged annulus; model a liquid-filled annu lus based on the as-logged annulus without solids therein ;

model a plu rality of cut casings at different cutting depths; calcu late a drag force {Fdrag) for extracting each of the plurality of cut casings based on the as-logged annu lus and the liquid-filled annulus;

determine an preferred cutting depth based on the Fdrag for each of the plurality of cut casings and specifications of casing retrieval equipment.

19. The system of claim 18, wherein Fdrag is calcu lated according to Equation 4 :

Fdrag ~ FF out.as log FS Out.as log F shear ~ FFout.liq Equation 4 where : FFoutiOLS log is a total flu id frictional force between a fluid contained within the annulus and an outer surface of the cut casing in the as- logged annu lus for the total length of the cut casing;

FSout S log is a total solid frictional force between solids contained within the annulus and an inner surface of the cut casing in the as-logged annulus for the total length of the cut casing ;

F shear is a total shear force to break the bond between a cement sheath in the annulus and the cut casing in the as-logged annu lus; and

FF0Utiuq is a total fluid frictional force between a fluid contained within the annu lus and the outer surface of the cut casing in the liqu id- filled annu lus for the total length of the cut casing .

Description:
BEHIND PIPE EVALUATION OF CUT AND PULL TENSION PREDICTION IN WELL ABANDONMENT AND INTERVENTION OPERATIONS

BACKGROUND

[0001] The exemplary embodiments described herein relate to casing removal operations.

[0002] The plug and abandonment (PA) of wells in offshore operations have two very important aspects addressed in the present disclosure. The critical issues associated with PA activities in offshore areas, especially in deep water, are safety and ecological environment protection. PA operations require cost and technologically effective methods. Determination of cutting depths and prediction or estimation of the forces required to pull and remove a portion of the casing lining a wellbore after cutting is becoming a very important set of information. The existing standard method to define free casing (also referred to herein as "pipe"), or pipe that is not bonded, is typically performed using a cement bond log (CBL) and often in conjunction with an ultrasonic device.

[0003] The evaluation of casing bond integrity is typically looked at for distinction between free pipe versus the extent of bonding of the pipe to determine a target cutting depth. The previous methods qualitatively define the least amount of force to pull the portion of the casing with minimized risk of having some stuck point above cutting depth. The determination of the depth to cut and pull a casing sometimes involves a trial and error process, which often leads to increased costs and time required to complete the project. The industry needs more accurate and predictable methods to determine optimal cutting depths.

BRIEF DESCRIPTION OF THE DRAWINGS

[0004] The following figures are included to illustrate certain aspects of the embodiments, and should not be viewed as exclusive embodiments. The subject matter disclosed is capable of considerable modifications, alterations, combinations, and equivalents in form and function, as will occur to those skilled in the art and having the benefit of this disclosure.

[0005] FIG. 1 illustrates a casing disposed in a wellbore. [0006] FIG. 2 illustrates a casing disposed in a wellbore where a cut line of the casing is located above a cement sheath.

[0007] FIG. 3 is a schematic process by which the analytical methods described herein may be executed .

[0008] FIG. 4 illustrates an exemplary plot of F tot S log and F totMq as a fu nction of wellbore depth.

[0009] FIG. 5 illustrates an exemplary plot of F drag as a function of wellbore depth.

[0010] FIG. 6 illustrates an exemplary system su itable for performing the analyses and methods of the present application.

[0011] FIGS. 7 and 8 illustrate crossplots of the predicted and actual data for the overpull and tension, respectively.

[0012] FIGS. 9 and 10 illustrate the calculated overpu ll and tension as a function of depth using the F totMq and the F tot S log , respectively.

DETAILED DESCRIPTION

[0013] The exemplary embodiments described herein relate to casing removal operations and, more specifically, to analytical models that estimate the total force to ( 1) break bonding between a cut casing and the annular materials contained between the cut casing and the su rrounding borehole and (2) pull the cut casing upward through the borehole. As used herein, the term "cut casing" refers to the portion of the casing from the cutting depth to the top of the casing (i. e. , the uphole end of the casing) . As used herein, the term "annular materials" refers to materials (e.g. , solids, liquids, gases, and combinations thereof) disposed in an annular space and does not imply a specific shape or configuration of the materials.

[0014] The analytical model of the present application uses measu rements from ( 1) nuclear sensors, (2) a cement bond log (CBL) tool, (3) ultrasonic devices, or (4) a combination thereof to evaluate and determine characteristics of annu lar materials disposed in the annu lus defined between casing and a wellbore wall . The measurements from these sensors and tools may provide information about the com position of the annular materials and the bonding characteristics of the annular materials and the casing outer wall . Exemplary annular materials may include, but are not limited to, solid cement, settled drilling fluid solids, mud, water, hydrates, light fluid, gas, and combinations thereof.

[0015] Each annular material contributes to the frictional forces experienced by the casing outer wall during removal of the cut casing from the borehole. By using the bonding information from CBL and/or ultrasonic measurements and the type of materials in the annulus as determined using the nuclear sensors, the shear strength of the annulus contents may be quantitatively determined, which allows for predicting the casing overpull (i.e. , amount of force needed to pull the cut casing upward and out of the borehole) and total tension needed to extract the cut casing from the borehole. The casing overpull and total tension for extraction of the cut casing may be used to determine, with higher accuracy than previously available, a cutting depth for producing the portion of the casing to be extracted.

[0016] Implementation of the analytical model described herein may be particularly useful for offshore casing removal operations and for preparing side tracking wells. Additionally, the analytical model and methods described herein may be applied to onshore shale formations for detecting shallow gas leaks to surface or fluid movement to shallow reservoirs including surface water layers. Further, the analytical model and methods described herein may be useful for identifying gas bubbling/flowing behind the casing, which is useful in preventing or mitigating environmental issues.

[0017] In some embodiments, the annular materials between the casing and the borehole may be analyzed with one or more of: a CBL tool, an ultrasonic device (e.g. , a circumferential acoustic scanning tool), or a nuclear sensor (e.g. , a neutron density sensor, a gamma-density sensor, a thermal neutron sensor, a dual-space neutron tool, a pulsed neutron sensor, or any combination thereof). For example, one or more CBL tools may be used for the methods and analyses described herein. In another example, one or more CBL tools in combination with one or more nuclear sensors may be used for the methods and analyses described herein. In yet another example, one or more ultrasonic devices in combination with one or more nuclear sensors may be used for the methods and analyses described herein. In another example, one or more nuclear sensors may be used for the methods and analyses described herein. [0018] The measurements may be performed by any suitable technique including wireline logging . In some instances where more than one measurement is performed, the CBL tool, ultrasonic device, and nuclear sensor may be on the same downhole tool or string of downhole tools. In alternate embodiments, the foregoing sensors and tools may be included on separate downhole tools, where the measurements are performed during different trips through the borehole.

[0019] The measu rements described herein in addition to the casing specifications, the fluid properties, and a borehole trajectory survey may be input into an analytical model of the present disclosu re to yield a suggested cutting depth (i. e. , a depth or location in a wellbore to cut the casing so as to produce the cut casing for extraction) . Then, in some instances, the casing may be cut in a location based on the suggested cutting depth, and the cut casing produced may be extracted and retrieved from the borehole.

[0020] FIG. 1 illustrates a casing 102 disposed in a wellbore 100. The casing 102 and wellbore 100 create an annu lus 104 with annu lar materials disposed therein . As illustrated, a portion of the annulus 104 may include solid cement 106 and another portion of the annu lus 104 may be filled with an outer fluid 108 (e.g. , settled drilling flu id solids, mud, water, hydrates, light fluid, gas, or combinations thereof). Fu rther, the casing 102 contains an inner fl uid 110. For clarity, the adjectives "outer" and "inner" refer to exterior and interior locations, respectively, with respect to a casing . FIG. 1 fu rther illustrates a cut line 112 at a specified cutting depth such that a cut casing 114 is in contact with the solid cement 106, the outer fluid 108, and the inner fluid 110.

[0021 ] FIG. 2 illustrates a casing 202 disposed in a wellbore 200 where a cut line 212 of the casing 202 is located above a cement sheath 206. Accordingly, a cut casing 214 in FIG. 2 is in contact with only an inner fluid 210 in the casing 202 and an outer fluid 208 (e.g. , settled drilling flu id solids, mud, water, hydrates, light fluid, gas, or combinations thereof) in an annulus 204 formed by the wellbore 200 and the casing 202.

[0022] The analytical model of the present disclosure estimates the total frictional force (F tot ) for extracting the cut casing from the wellbore for the two scenarios illustrated in FIGS. 1 and 2. First, ( 1) nuclear sensors, (2) a CBL tool, (3) u ltrasonic devices, or (4) a combination thereof evaluate the location and composition of annu lar materials. Then, the F tot may be calculated for a one or more cutting depths for the two scenarios: ( 1) an as-logged annulus where the composition and location of the annular materials from the foregoing analyses and the cutting depth is at a location along the solid cement as illustrated in FIG. 1 {F tot as log ) and (2) a liquid-filled annulus at the same cutting depth but assuming that all annular materials down to the cutting depth are outer fluids as illustrated in FIG. 2 (F tot liq ) . Therefore, the two scenarios calculate the F tot for the two extremes where the F tot S log is the greatest expected force to extract the cut casing from the well bore and the F tot q is the least expected force, which may each calculated according to Equations 1 and 2, respectively.

Ftot.as log ~ FW CS + FFi n as i 0 g + FF 0Ut as i 0 g + FS 0Ut as i 0 g + Fshear Equation 1 where: EW CS is the total effective casing weight (Equation 3);

PPin,as iog is a total fluid frictional force between a fluid contained within the cut casing 114 (i.e. , the inner fluid 110) and an inner surface 116 of the cut casing 114 in the as-logged annulus 104 for a total length of the cut casing 114;

FFout.as iog is a total fluid frictional force between a fluid contained within the annulus 104 (i.e., the outer fluid 108) and an outer surface 118 of the cut casing 114 in the as-logged annulus 104 for the total length of the cut casing 114;

FS out S log is a total solid frictional force between solids contained within the annulus 104 (e.g., settled drilling fluid solids, the cement sheath 106, and the like) and the outer surface 118 of the cut casing 114 in the as-logged annulus for the total length of the cut casing 114;

F shear is a total shear force to break the bond between the cement sheath 106 in the annulus 104 and the cut casing 114 in the as- logged annulus 104;

Ffot.liq ~ EW CS + FF inMq + FF 0UtMq Equation 2 where: FF inMq is a total fluid frictional force between a fluid contained within the cut casing 214 (i.e. , the inner fluid 210) and the inner surface 216 of the cut casing 214 in the liquid-filled annulus for the total length of the cut casing 214; and

PPout.uq is a total fluid frictional force between a fluid contained within the annu lus 204 (i. e. , the outer fluid 208) and the outer surface 218 of the cut casing 214 in the liquid-filled annulus for the total length of the cut casing 214;

EW CS = GW CS * BF Equation 3 where : GW CS is gross weight of the cut casing ; and

BF is the is the buoyance factor of the cut casing

[0023] Some embodiments may involve modeling the liqu id-filled annulus based on the as-logged annulus without solids therein. Further, some embodiments may involve modeling a plurality of cut casings at different cutting depths and calculating F tot S log and F totMq for each of the modeled cut casing . Described fu rther below, the F tot S log and F totMq for each of the modeled cut casing along with specifications of the casing retrieval equipment may be used to determine a preferred cutting depth.

[0024] Referring again to the foregoing equations, Equation 2 for FtotMq does not include FS 0Ut or F shear terms because in the liqu id-filled annulus model, solids are not present along the cut casing 214.

[0025] The difference between F tot S log and F totMq is the total drag force (F drag ) for extracting the cut casing from the wellbore and is given in Equation 4. Because the frictional forces from the inner fluids 110,210 are the same (i. e. , FF in S log = FF inMq ), the terms FF in S log and FF inMq and each of the EW CS are not present in the expanded portion of Equation 4.

^drag ~ F tot,as log ~ F tot,liq — FF out,as i og + FS 0Ut CLS i og + F sf t e∞ - FF 0Ut n q

Equation 4

[0026] Each of FF in as log , FF in Uq , FF out S log , FS out as log , F shear , and FF out,uq ma y be calculated according to Equations 5-9. FFin.as log = FF in q = ∑ * (pi * g * z) * (π * d t * Az) Equation 5 where: a fl i is the friction coefficient of the inner fluid;

i is the density of the inner fluid;

g is the acceleration of gravity;

z is the depth of the cut line;

d t is the inner diameter of the cut casing; and

Az is the length of the cut casing

FF 0Ut a S log = ∑ a fl}0 * (p 0 * g * z) * (1 - S cv ) * (π * d 0 * Δζ) Equation 6 where: a fl 0 is the friction coefficient of the outer fluid;

Po is the density of the outer fluid;

d 0 is the outer diameter of the cut casing; and

S cv is the fraction Az of the cut casing contacting the outer solid (i. e. , the solid in the annulus)

FS 0Ut a S log = ∑ a so i * ( o * g * z) * (S cv ) * (π * d 0 * Δζ) Equation 7 where: a soi is the friction coefficient of the outer solid

Fsnear = ∑f «sh * S„ * 5/ * Δζ Equation 8 where: a sft is the shear coefficient for the annular material; and

SI is the shear index as a function of impedance, bulk density of the annular material, and the thickness of the annular material

FFout.uq = ∑ a fi,o * (Po * g * z) * (π * d 0 * Az) Equation 9

[0027] Some embodiments may involve modeling the liquid-filled annulus based on the as-logged annulus without solids therein. Further, some embodiments may involve modeling a plurality of cut casings at different cutting depths and calculating F drag based on the as-logged annulus and the liquid-filled annulus for each of the modeled cut casing. Described further below, F drag for each of the modeled cut casing along with specifications of the casing retrieval equipment may be used to determine a preferred cutting depth.

[0028] In the conditions where wells are constructed with multiple casings concentrically-arranged within one another, the above equations and analytical models described herein may be updated to accou nt for the contents of a second annulus defined between concentric casings and, if applicable, any deeper annu lar regions.

[0029] The foregoing description of the analytical models disclosed herein are for vertical portions of a wellbore, but may be adapted for deviated wellbore. To do so, the friction force for deviated portions of the wellbore may be projected in two components : ( 1) the horizontal component and (2) the vertical component. The total friction force will be the resultant of those two components.

[0030] FIG. 3 illustrates a process by which the analytical methods described herein may be executed . Inputs may include the casing specifications {e.g. , d i t z, Δζ, and d 0 ), the fluid properties {e.g. , p 0 and p t ), the well trajectory survey, and the sensor measu rements {e.g. , the S cv may be measured by the ultrasonic device). Regarding the various frictional coefficients, each may be estimated based on the compositions to which the particular frictional coefficient relates. Fu rther, in some instances, data from previous cut and pull operations may be used to u pdate the various frictional coefficients using error minimization for a best fit cu rve character.

[0031] The analytical model may produce one or more outputs. In some instances, the analytical model may output one or more of the frictional forces described above {F drag , F tot S log , or F totMq ) . In some instances, the analytical model may use the frictional forces to derive the tension and/or overpull associated with the cut and pu ll operation. Tension and overpull may be calcu lated by a variety of known methods that take into account the frictional forces {F drag , F tot as log , or F tot Uq ) .

[0032] In some embodiments, the methods described herein may execute the analytical model using a variety of cutting depths. Then, the output may be the calculated frictional forces, tension, and/or overpull as a function of depth. For example, FIG. 4 illustrates an exemplary plot of F tot as log and F totMq as a function of wellbore depth. As illustrated, the two frictional forces are essentially equal until reaching the top of the cement (horizontal dotted like), then the F tot S log increases at a greater rate than the F totMq because the all-liquid scenario has less frictional force with no solid frictional or shear components. FIG. 5 illustrates the corresponding F drag , which as described above, is the difference between the F tot S log and F totMq . Similarly, the method may involve plotting tension and/or overpull as a function of depth for one or more of the frictional forces.

[0033] In cut and pull operations, it is typically desired to cut at the deepest depth possible taking into consideration the specifications for the casing retrieval equipment (e.g. , the hoisting apparatus, conveyance, platform, and the like) and safety (e.g., operati ng at no more than 10% less of the maximum casing retrieval equipment specifications) . In some instances, a threshold value for the frictional force (F thresh ) based on the rig capacity for maximum pu ll within safety limits, the tension, and/or the overpu ll may be determined based on the safety considerations and the casing retrieval equ ipment specifications. The cutting depth may then be determined such that F tot S log < F thresh and/or Fdrag + Ftot.uq < F thresh - Similarly, the cutting depth may be determined such that the threshold tension and/or overpu ll is between a minimum tension and/or overpull is based on the F tot}Uq and a maximum tension and/or overpu ll is based on the F tot S log . That is, in some embodiments, determining the preferred cutting depth may involve calculating a tension for the as-logged annulus (T tot S i og ) using F tot S log , a tension for the liqu id-filled annulus (T totMq ) using F totMq , a tension for the drag force (T drag ) using F drag , and/or a threshold tension (T thresh ) based on the specifications of casing retrieval equipment, wherein the preferred cutting depth provides for T tot S log < T thresh and/or T drag + T tot q < T thresh . In some instances, determining the preferred cutting depth may involve calculating calcu lating an overpull for the as-logged annulus (OP tot S log ) using F tot S log , an overpull for the liquid-filled annu lus (OP tot liq ) using F tot liq , an overpu ll for the drag force (OP drag ) using F drag , and/or a threshold overpu ll (OP thresh ) based on the specifications of casing retrieval equ ipment, wherein the preferred cutting depth provides for OP tot S log < OP thresh and/or OP drag + OP tot q < OP thresh .

[0034] In some embodiments, the mud may be circulated within the wellbore after cutting and before pulling the casing . This may decrease the total tension requ ired for pu lling the casing because some of the solid may be lifted to the surface during the mud circulation .

[0035] FIG. 6 provides an exemplary illustration of a system 610 suitable for performing the analyses and methods described herein, according to at least some embodiments. While this example is illustrated as an offshore system 610, those skilled in the art will recognize the applicability and corresponding modification for onshore systems, without departing from the scope of the disclosure.

[0036] As illustrated, a semi-su bmersible platform 612 is centered over a submerged oil and gas formation 614 located below sea floor 616. A subsea conduit 618 extends from a deck 620 of the platform 612 to a wellhead installation 622 including su bsea blow-out preventers 624. The platform 612 has a hoisting apparatus 626 and a derrick 628 for conveying a tool 632 coupled to a conveyance 630 through a wellbore 636 penetrating the su bterranean formation 614. Although depicted as a wireline, conveyance 630 may be, but is not limited to, a wireline, a coiled tubing, a slickline, a cable, a downhole tractor, and the like.

[0037] Annular materials 638 are disposed between the casing 634 and the wellbore 636. As illustrated, the annular materials 638 transition from a liquid 640 to cement 644 at a point along the wellbore 636. The tool 632 may be used to measu re a property of the annular materials 638 along a length or portion of the wellbore 636. The measu rements may be conveyed to a control system 648 and used as inputs for the analytical model described herein for determining a cutting line. Then, the casing 634 may be cut at or nea r the cutting line and pulled to the platform 612 (or the su rface for a land-based system) .

[0038] The control system(s) 648 and corresponding computer hardware used to implement the various illustrative blocks, modules, elements, components, methods, and algorithms described herein can include a processor configu red to execute one or more sequences of instructions, programming stances, or code stored on a non-transitory, computer-readable mediu m . The processor can be, for example, a general pu rpose microprocessor, a microcontroller, a digital signal processor, an application specific integrated circuit, a field programmable gate array, a programmable logic device, a controller, a state machine, a gated logic, discrete hardware components, an artificial neural network, or any like suitable entity that can perform calculations or other manipulations of data. In some embodiments, computer hardware can further include elements such as, for example, a memory (e.g., random access memory (RAM), flash memory, read only memory (ROM), programmable read only memory (PROM), erasable programmable read only memory (EPROM)), registers, hard disks, removable disks, CD-ROMS, DVDs, or any other like suitable storage device or mediu m .

[0039] Executable sequences described herein can be implemented with one or more sequences of code contained in a memory. In some embodiments, such code can be read into the memory from another machine- readable mediu m. Execution of the sequences of instructions contained in the memory can cause a processor to perform the process steps described herein. One or more processors in a multi-processing arrangement can also be employed to execute instruction sequences in the memory. In addition, hardwired circu itry can be used in place of or in combination with software instructions to implement various embodiments described herein. Thus, the present embodiments are not limited to any specific combination of hardware and/or software.

[0040] As used herein, a machine-readable mediu m will refer to any medium that directly or indirectly provides instructions to a processor for execution . A machine-readable mediu m can take on many forms including, for example, non-volatile media, volatile media, and transmission media. Nonvolatile media can include, for example, optical and magnetic disks. Volatile media can include, for example, dynamic memory. Transmission media can include, for example, coaxial cables, wire, fiber optics, and wires that form a bus. Common forms of machine-readable media can include, for example, floppy disks, flexible disks, hard disks, magnetic tapes, other like magnetic media, CD- ROMs, DVDs, other like optical media, punch cards, paper tapes and like physical media with patterned holes, RAM, ROM, PROM, EPROM and flash EPROM .

[0041 ] Embodiments described herein may include Embodiments A, B, C, or D .

[0042] Embodiment A is a method that includes measu ring at least one characteristic of fluids and solids disposed in an annulus defined between a casing and a wellbore with at least one of: a cement bond log tool, an ultrasonic device, or a neutron sensor, and thereby providing an as-logged annulus; modeling a liqu id-filled annu lus based on the as-logged annu lus without solids therein ; modeling a plurality of cut casings at different cutting depths; calculating a total force for the as-logged annu lus {F tot S i og ) for extracting each of the plurality of cut casings; calculating a total force for the liqu id-filled annulus {F tot q ) for extracting each of the plurality of cut casings; and determining a preferred cutting depth based on the F totias log for each of the plurality of cut casings, the F totMq for each of the plurality of cut casings, and specifications of casing retrieval equipment.

[0043] Embodiment B is system that includes a conveyance extendable into a wellbore lined at least partially with a casing, wherein the conveyance is coupled to at least one of: a cement bond log tool, an ultrasonic device, or a neutron sensor, and wherein fluids and solids are contained in an annulus defined between the casing and the wellbore; and a control system that includes a non-transitory medium readable by a processor and storing instructions that when executed by the processor cause the control system to: measure at least one characteristic of the fluids and solids the at least one of: the cement bond log tool, the ultrasonic device, or the neutron sensor, and thereby providing an as-logged annulus; model a liquid-filled annulus based on the as-logged annulus without solids therein; model a plurality of cut casings at different cutting depths; calculate a total force for the as-logged annulus ( F tot,as iog) for extracting each of the plurality of cut casings; calculate a total force for the liquid-filled annulus (F tot q ) for extracting each of the plurality of cut casings; and determine a preferred cutting depth based on the F tot S log for each of the plurality of cut casings, the F tot q for each of the plurality of cut casings, and specifications of casing retrieval equipment.

[0044] Embodiments A and B may further include one or more of the following : Element 1 : the method further including performing (or the instruction further causing the control system to perform) a casing removal operation at the preferred cutting depth with the casing retrieval equipment; Element 2: wherein F t ot,as iog is calculated according to Equation 1; Element 3 : wherein F totMq is calculated according to Equation 2; Element 4 : the method further including calculating (or the instruction further causing the control system to calculate) a drag force (F drag ) according to Equation 4, wherein determining the preferred cutting depth is further based on F drag ; Element 5 : the method further including Element 4 and calculating (or the instruction further causing the control system to perform Element 4 and calculate) a threshold force (F thresh ) based on the specifications of casing retrieval equipment, and wherein the preferred cutting depth provides for F tot S log < F thresh and/or F drag + F totMq < F thresh ; Element 6: the method further including Element 5 and calculating (or the instruction fu rther causing the control system to perform Element 5 and calcu late) a tension for the as-logged annulus {T to t,as io g ) using F tot S log i a tension for the liqu id-filled annulus (TtotMq) u si n 9 F tot,uq> a tension for the drag force {T drag ) using F drag i and a threshold tension (T thresh ) based on the specifications of casing retrieval equipment, wherein the preferred cutting depth provides for T totiOS log < T thresh and/or T drag + T tot liq < T thresh ; and Element 7 : the method further including Element 5 and calculating (or the instruction fu rther causing the control system to perform Element 5 and calculate) an overpull for the as-logged annulus (OP t ot,as iog) using F tot S log i an overpull for the liquid-filled annu lus {OP tot liq ) using F tot,uq> a n overpull for the drag force {OP drag ) using F drag , and a threshold overpull (OP thresh ) based on the specifications of casing retrieval equ ipment, wherein the preferred cutting depth provides for OP tot S i og < OP thresh and/or OP d rag + OPtot.iiq < 0P t h res h - Exemplary combinations of elements may include, but are not limited to, Element 1 in combination with Elements 2 and/or 3 and optionally fu rther in combination with at least one of Elements 4-7; Elements 2 and 3 in combination and optionally further in combination with at least one of Elements 4-7; Element 1 in combination with at least one of Elements 4-7; and at least two of Elements 4-7 in combination.

[0045] Embodiment C is a method that includes measuring at least one characteristic of fluids and solids disposed in an annu lus formed by a casing and a wellbore with at least one of: a cement bond log, an ultrasonic device, or a neutron sensor, thereby providing an as-logged annulus; modeling a liqu id-filled annulus based on the as-logged annulus without solids therein ; modeling a plu rality of cut casings at different cutting depths; calcu lating a drag force ( F d rag) for extracting each of the plurality of cut casings based on the as-logged annulus and the liquid-filled annulus; determining an preferred cutting depth based on the F drag for each of the plurality of cut casings and specifications of casing retrieval equipment.

[0046] Embodiment D is a system that includes a conveyance extendable into a wellbore lined at least partially with a casing, wherein the conveyance is coupled to at least one of: a cement bond log tool, an ultrasonic device, or a neutron sensor, and wherein fluids and solids are contained in an annulus defined between the casing and the wellbore; and a control system that includes a non-transitory medium readable by a processor and storing instructions that when executed by the processor cause the control system to: measure at least one characteristic of fluids and solids disposed in an annulus formed by a casing and a wellbore with at least one of: a cement bond log, an ultrasonic device, or a neutron sensor, thereby providing an as-logged annulus; model a liquid-filled annulus based on the as-logged annulus without solids therein; model a plurality of cut casings at different cutting depths; calculate a drag force (F drag ) for extracting each of the plurality of cut casings based on the as-logged annulus and the liquid-filled annulus; and determine an preferred cutting depth based on the F drag for each of the plurality of cut casings and specifications of casing retrieval equipment.

[0047] Embodiments C and D may further include one or more of the following : Element 1 ; Element 8: wherein F drag is calculated according to Equation 4; Element 9 : wherein the method further comprises calculating (or the instruction further causing the control system to calculate) a threshold force (Ft h res h ) based on the specifications of casing retrieval equipment, wherein the preferred cutting depth provides for F drag < F thresh ; Element 10 : wherein the method further comprises calculating (or the instruction further causing the control system to calculate) a tension {T drag ) using F drag and a threshold tension (Tt h res h ) based on the specifications of casing retrieval equipment, wherein the preferred cutting depth provides for T drag < T thresh ; and Element 11 : wherein the method further comprises calculating (or the instruction further causing the control system to calculate) an overpull {OP drag ) using F drag and a threshold overpull (OP thresh ) based on the specifications of casing retrieval equipment, and wherein the preferred cutting depth provides for OP drag < OP thresh . Exemplary combinations of elements may include, but are not limited to, Element 1 in combination with Element 8 and optionally further in combination with at least one of Elements 9- 11; Element 1 in combination with at least one of Elements 9- 11; Element 8 in combination with at least one of Elements 9-11; and two or more of 9-11 in combination.

[0048] Unless otherwise indicated, all numbers expressing quantities of ingredients, properties such as molecular weight, reaction conditions, and so forth used in the present specification and associated claims are to be understood as being modified in all instances by the term "about." Accordingly, unless indicated to the contrary, the numerical parameters set forth in the following specification and attached claims are approximations that may vary depending upon the desired properties sought to be obtained by the embodiments of the present disclosure. At the very least, and not as an attempt to limit the application of the doctrine of equivalents to the scope of the claim, each numerical parameter should at least be construed in light of the number of reported significant digits and by applying ordinary rounding techniques.

[0049] One or more illustrative embodiments incorporating the embodiments disclosed herein are presented herein. Not all features of a physical implementation are described or shown in this application for the sake of clarity. It is understood that in the development of a physical embodiment incorporating the embodiments of the present disclosure, numerous implementation-specific decisions must be made to achieve the developer's goals, such as compliance with system-related, business-related, government- related and other constraints, which vary by implementation and from time to time. While a developer's efforts might be time-consuming, such efforts would be, nevertheless, a routine undertaking for those of ordinary skill the art and having benefit of this disclosure.

[0050] While compositions and methods are described herein in terms of "comprising" various components or steps, the compositions and methods can also "consist essentially of" or "consist of" the various components and steps.

[0051] To facilitate a better understanding of the embodiments of the present disclosure, the following examples of preferred or representative embodiments are given. In no way should the following examples be read to limit, or to define, the scope of the invention.

EXAMPLES

[0052] Linear Programing from MATLAB® (a high-level language and interactive environment for numerical computation, visualization, and programming, available from MathWorks) was used to solve the modeling logic in order to find the friction coefficients for FF in S log i FF inMq i FF out S log i FS out S log i Fs h ear, and FF 0UtMq described above in Equations 5-9. The basic equation was written in the form of matrix as Equation 10, where Y is the matrix output, A is the matrix input, and x is the friction coefficient that is to be optimized. Y = A - x Equation 10

[0053] When optimizing x and fl 0 , each of x and fl 0 were varied as a function of the inner fluid viscosity and the outer fluid viscosity, respectively. Similarly, a sol also can be separated individually for multiple solids due to mud settling in combination with other solid types such as cement. The a sh was likewise individually separated for solid settling and other modeled solid types such as cement. All the friction coefficients were optimized by linear programing using MATLAB® software with actual existing data from the cut and pull operations in a historical history matching manner.

[0054] Comparison of actual overpull and total tension applied when cutting and pulling casing for abandonment (7 wells) compared with the prediction calculations from the analytical model with the friction coefficients derived as described above. FIGS. 7 and 8 are crossplots of the predicted and actual data for the overpull and tension, respectively. These crossplots illustrate that the analytical model estimates overpull and total tension values accurately, which is evidenced by the high values for R 2 .

[0055] Using the derived frictional coefficients, the overpull and tension was calculated for a liquid-filled annulus and an as-logged annulus. FIGS. 9 and 10 provide the calculated overpull and tension as a function of depth using the F totMq (liquid-filled annulus) and the F tot S log (as-logged annulus), respectively. This illustrates the graphical outputs from the analytical model described herein that may be used in determining a cutting depth.

[0056] Therefore, the present invention is well adapted to attain the ends and advantages mentioned as well as those that are inherent therein. The particular embodiments disclosed above are illustrative only, as the present invention may be modified and practiced in different but equivalent manners apparent to those skilled in the art having the benefit of the teachings herein. Furthermore, no limitations are intended to the details of construction or design herein shown, other than as described in the claims below. It is therefore evident that the particular illustrative embodiments disclosed above may be altered, combined, or modified and all such variations are considered within the scope and spirit of the present invention. The invention illustratively disclosed herein suitably may be practiced in the absence of any element that is not specifically disclosed herein and/or any optional element disclosed herein. While compositions and methods are described in terms of "comprising," "containing," or "including" various components or steps, the compositions and methods can also "consist essentially of" or "consist of" the various components and steps. All numbers and ranges disclosed above may vary by some amount. Whenever a numerical range with a lower limit and an upper limit is disclosed, any number and any included range falling within the range is specifically disclosed. In particular, every range of values (of the form, "from about a to about b," or, equivalently, "from approximately a to b," or, equivalently, "from approximately a-b") disclosed herein is to be understood to set forth every number and range encompassed within the broader range of values. Also, the terms in the claims have their plain, ordinary meaning unless otherwise explicitly and clearly defined by the patentee. Moreover, the indefinite articles "a" or "an," as used in the claims, are defined herein to mean one or more than one of the element that it introduces.