CA2625178A1 - Methods and computer-readable media for determining design parameters to prevent tubing buckling in deviated wellbores - Google Patents
Methods and computer-readable media for determining design parameters to prevent tubing buckling in deviated wellbores Download PDFInfo
- Publication number
- CA2625178A1 CA2625178A1 CA002625178A CA2625178A CA2625178A1 CA 2625178 A1 CA2625178 A1 CA 2625178A1 CA 002625178 A CA002625178 A CA 002625178A CA 2625178 A CA2625178 A CA 2625178A CA 2625178 A1 CA2625178 A1 CA 2625178A1
- Authority
- CA
- Canada
- Prior art keywords
- tubing
- parameter
- predicting
- buckling
- deviated wellbore
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 23
- 238000013461 design Methods 0.000 title claims abstract description 17
- 238000005452 bending Methods 0.000 claims abstract description 44
- 239000003129 oil well Substances 0.000 claims abstract description 16
- 238000006073 displacement reaction Methods 0.000 claims description 6
- 238000005096 rolling process Methods 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 3
- 241000307294 Diaphus theta Species 0.000 claims 1
- 238000003860 storage Methods 0.000 description 17
- 238000004364 calculation method Methods 0.000 description 9
- 238000012545 processing Methods 0.000 description 7
- 238000004590 computer program Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 3
- 101150109831 SIN4 gene Proteins 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000005553 drilling Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B17/00—Drilling rods or pipes; Flexible drill strings; Kellies; Drill collars; Sucker rods; Cables; Casings; Tubings
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/02—Subsoil filtering
- E21B43/10—Setting of casings, screens, liners or the like in wells
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B41/00—Equipment or details not covered by groups E21B15/00 - E21B40/00
Abstract
Methods and computer-readable media are provided for determining design parameters for oil well casing and tubing to prevent buckling in deviated wellbores. Well parameter data including tubing size, tubing weight, well depth, and well geometry is obtained and may be utilized to calculate parameters for predicting the movement of tubing near a packer or centralizer in the deviated wellbore based on the received well parameter data, predicting a total bending moment near the packer or centralizer, predicting a maximum bending stress near the packer or centralizer based on the total bending moment, and predicting the minimum axial force necessary to initiate buckling due to friction, and predicting the onset of buckling for the connection of tubing of different sizes. After the parameters have been calculated, they may be utilized in the design of the oil well casing and tubing to prevent buckling in the deviated wellbore.
Description
METHODS AND COMPUTER-READABLE MEDIA FOR DETERMINING
DESIGN PARAMETERS TO PREVENT TUBING BUCKLING IN
DEVIATED WELLBORES
This application is being filed on 03 October 2006, as a PCT International Patent application in the name of Landmark Graphics Corporation, a U.S.
national corporation, application for the designation of all countries except the US, and Robert F. Mitchell, a U.S. citizen, applicant for the designation of the US
only, and claims priority to U.S. Provisional Patent Application Serial No. 60/723,513, filed October 4, 2005, and U.S. Utility Patent Application No. 11/274,637, filed Noveinber 15, 2005. Both of the aforementioned patent applications are assigned to the same assignee as this application and are expressly incorporated herein by reference.
TECHNICAL FIELD
The present invention is related to the analysis of oil well casing and pipe or tubing buckling caused by critical loading in a wellbore. More particularly, the present inven.tion is related to the accurate determination of critical loading parameters in the design of oil well tubing to prevent buckling in deviated wellbores.
BACKGROUND
In an oil well, casing is typically installed to withstand various pressures which may be present in an open hole or wellbore and to stabilize the pipes or tubing used for drilling. Typically, casing hangs straight down in vertical wells or lies on the low side of the hole in deviated wells. During drilling operations, thermal or pressure loads within a wellbore may produce compressive loads which, if sufficiently high, will cause the initial well configuration to become unstable.
However, since the tubing is confined within the casing (or alternatively an open hole), the tubing can deform into anotller stable configuration, which may be a helical or coil shape in a vertical well or a lateral "S" shaped configuration in a deviated well. The change to the new configurations caused by the deformed tubing is known as "buckling."
In tubing and casing design, the accurate analysis of buckling is important for several reasons. First, buckling generates bending stresses not present in the original configuration. If the stresses in the original (i.e., "unbuckled") configuration were near yield, additional stress could produce failure in the tubing, including permanent plastic deformation called "corkscrewing." Second, buckling causes movement in oil well tubing. That is, buckled tubing (which is coiled) is shorter than straight tubing, and this is an important consideration if the tubing is not fixed. Third, tubing buckling causes the relief of compressive axial loads when the casing surrounding the tubing is fixed.
Previously, models have been developed for analyzing buckling in wells, however, these models suffer from several drawbacks when applied to deviated wells. One drawback with previous models is that tubing bending stress due to buckling will be overestimated for deviated wells. Another drawback with previous models as applied to deviated wells is that they over predict tubing movement.
Still another drawback with previous models is that tubing compliance is overestimated, which may greatly underestimate the axial loads able to be withstood by the surrounding casing. It is with respect to these considerations and others that the various embodiments of the present invention have been made.
SUMMARY
Illustrative embodiments of the present invention address these issues and others by providing a method of determining design parameters for oil well casing and tubing to prevent buckling in a deviated wellbore. According to the method, well parameter data is received which may include tubing size, tubing weight, well deptll, and well geometry. The method further includes calculating a first parameter for predicting the movement of tubing near at least one boundary condition in the deviated wellbore based on the received well parameter data. The boundary condition may be a packer installed in the deviated wellbore, a centralizer installed in the deviated wellbore, or both.
The method further includes calculating a second parameter for predicting a total bending moment near the at least one boundary condition, calculating a third parameter for predicting a maximum bending stress near the at least one boundary condition in the deviated wellbore based on the total bending moment, and calculating a fourth parameter for predicting the minimum axial force necessary to initiate buckling due to friction, based on the received well parameter data.
The method may fiurther include calculating a fifth parameter for predicting the onset of buckling for the connection of tubing of different sizes (i.e., tapered strings) based on the received well parameter data. After the first, second, third, fourth, and fifth parameters have been calculated, they may be utilized in the design of the oil well casing and tubing to prevent buckling in the deviated wellbore.
Other illustrative embodiments of the invention may also be implemented in a computer system or as an article of manufacture such as a computer program product or computer readable media. The computer program product may be a computer storage media readable by a computer system and encoding a computer program of instructions for executing a computer process. The computer program product may also be a propagated signal on a carrier readable by a computing system and encoding a computer program of instructions for executing a computer process.
These and various other features, as well as advantages, which characterize the present invention, will be apparent from a reading of the following detailed description and a review of the associated drawings.
DESCRIPTION OF THE DRAWINGS
FIGURE 1 shows a typical computer system operating environment for illustrative embodiments of the present invention.
FIGURE 2 shows logical operations performed by an illustrative embodiment for calculating a parameter for predicting the movement of tubing near at least one boundary condition in a deviated wellbore.
FIGURE 3 shows logical operations performed by an illustrative embodiment for calculating parameters for predicting total bending moments and maximum bending stresses near at least one boundary condition in a deviated wellbore.
FIGURE 4 shows logical operations performed by an illustrative embodiment for calculating parameters for predicting minimum axial forces necessary to initiate buclcling due to friction in a deviated welibore.
FIGURE 5 shows logical operations performed by an illustrative embodiment for calculating a parameter for predicting the onset of buckling for the connection of tubing of different sizes in a deviated wellbore.
DESIGN PARAMETERS TO PREVENT TUBING BUCKLING IN
DEVIATED WELLBORES
This application is being filed on 03 October 2006, as a PCT International Patent application in the name of Landmark Graphics Corporation, a U.S.
national corporation, application for the designation of all countries except the US, and Robert F. Mitchell, a U.S. citizen, applicant for the designation of the US
only, and claims priority to U.S. Provisional Patent Application Serial No. 60/723,513, filed October 4, 2005, and U.S. Utility Patent Application No. 11/274,637, filed Noveinber 15, 2005. Both of the aforementioned patent applications are assigned to the same assignee as this application and are expressly incorporated herein by reference.
TECHNICAL FIELD
The present invention is related to the analysis of oil well casing and pipe or tubing buckling caused by critical loading in a wellbore. More particularly, the present inven.tion is related to the accurate determination of critical loading parameters in the design of oil well tubing to prevent buckling in deviated wellbores.
BACKGROUND
In an oil well, casing is typically installed to withstand various pressures which may be present in an open hole or wellbore and to stabilize the pipes or tubing used for drilling. Typically, casing hangs straight down in vertical wells or lies on the low side of the hole in deviated wells. During drilling operations, thermal or pressure loads within a wellbore may produce compressive loads which, if sufficiently high, will cause the initial well configuration to become unstable.
However, since the tubing is confined within the casing (or alternatively an open hole), the tubing can deform into anotller stable configuration, which may be a helical or coil shape in a vertical well or a lateral "S" shaped configuration in a deviated well. The change to the new configurations caused by the deformed tubing is known as "buckling."
In tubing and casing design, the accurate analysis of buckling is important for several reasons. First, buckling generates bending stresses not present in the original configuration. If the stresses in the original (i.e., "unbuckled") configuration were near yield, additional stress could produce failure in the tubing, including permanent plastic deformation called "corkscrewing." Second, buckling causes movement in oil well tubing. That is, buckled tubing (which is coiled) is shorter than straight tubing, and this is an important consideration if the tubing is not fixed. Third, tubing buckling causes the relief of compressive axial loads when the casing surrounding the tubing is fixed.
Previously, models have been developed for analyzing buckling in wells, however, these models suffer from several drawbacks when applied to deviated wells. One drawback with previous models is that tubing bending stress due to buckling will be overestimated for deviated wells. Another drawback with previous models as applied to deviated wells is that they over predict tubing movement.
Still another drawback with previous models is that tubing compliance is overestimated, which may greatly underestimate the axial loads able to be withstood by the surrounding casing. It is with respect to these considerations and others that the various embodiments of the present invention have been made.
SUMMARY
Illustrative embodiments of the present invention address these issues and others by providing a method of determining design parameters for oil well casing and tubing to prevent buckling in a deviated wellbore. According to the method, well parameter data is received which may include tubing size, tubing weight, well deptll, and well geometry. The method further includes calculating a first parameter for predicting the movement of tubing near at least one boundary condition in the deviated wellbore based on the received well parameter data. The boundary condition may be a packer installed in the deviated wellbore, a centralizer installed in the deviated wellbore, or both.
The method further includes calculating a second parameter for predicting a total bending moment near the at least one boundary condition, calculating a third parameter for predicting a maximum bending stress near the at least one boundary condition in the deviated wellbore based on the total bending moment, and calculating a fourth parameter for predicting the minimum axial force necessary to initiate buckling due to friction, based on the received well parameter data.
The method may fiurther include calculating a fifth parameter for predicting the onset of buckling for the connection of tubing of different sizes (i.e., tapered strings) based on the received well parameter data. After the first, second, third, fourth, and fifth parameters have been calculated, they may be utilized in the design of the oil well casing and tubing to prevent buckling in the deviated wellbore.
Other illustrative embodiments of the invention may also be implemented in a computer system or as an article of manufacture such as a computer program product or computer readable media. The computer program product may be a computer storage media readable by a computer system and encoding a computer program of instructions for executing a computer process. The computer program product may also be a propagated signal on a carrier readable by a computing system and encoding a computer program of instructions for executing a computer process.
These and various other features, as well as advantages, which characterize the present invention, will be apparent from a reading of the following detailed description and a review of the associated drawings.
DESCRIPTION OF THE DRAWINGS
FIGURE 1 shows a typical computer system operating environment for illustrative embodiments of the present invention.
FIGURE 2 shows logical operations performed by an illustrative embodiment for calculating a parameter for predicting the movement of tubing near at least one boundary condition in a deviated wellbore.
FIGURE 3 shows logical operations performed by an illustrative embodiment for calculating parameters for predicting total bending moments and maximum bending stresses near at least one boundary condition in a deviated wellbore.
FIGURE 4 shows logical operations performed by an illustrative embodiment for calculating parameters for predicting minimum axial forces necessary to initiate buclcling due to friction in a deviated welibore.
FIGURE 5 shows logical operations performed by an illustrative embodiment for calculating a parameter for predicting the onset of buckling for the connection of tubing of different sizes in a deviated wellbore.
DETAILED DESCRIPTION
Illustrative embodiments of the present invention provide for determining design parameters for oil well casing and tubing to prevent buckling in a deviated wellbore. Referring now to the drawings, in which like numerals represent like elements, various aspects of the present invention will be described. In particular, FIGURE 1 and the corresponding discussion are intended to provide a brief, general description of a suitable computing environment in which embodiments of the invention may be iinplemented. While the invention will be described in the general context of program modules that execute in conjunction with program modules that run on an operating system on a personal computer, those skilled in the art will recognize that the invention may also be implemented in combination with other types of computer systems and program modules.
Generally, program modules include routines, programs, components, data structures, and other types of structures that perform particular tasks or implement particular abstract data types. Moreover, those skilled in the art will appreciate that the invention may be practiced with other computer systein configurations, including hand-held devices, multiprocessor systems, inicroprocessor-based or programmable consumer electronics, minicomputers, mainframe computers, and the like. The invention may also be practiced in distributed computing environments where tasks are performed by remote processing devices that are linked through a communications network. In a distributed computing environment, program modules may be located in both local and remote memory storage devices.
Referring now to FIGURE 1, an illustrative computer architecture for a computer 2 utilized in the various einbodiments of the invention will be described.
The computer architecture shown in FIGURE 1 illustrates a conventional desktop or laptop computer, including a central processing unit 5 ("CPU"), a system memory 7, including a random access memory 9 ("RAM") and a read-only memory ("ROM") 11, and a system bus 12 that couples the memory to the CPU 5. A basic input/output system containing the basic routines that help to transfer information between elements within the computer, such as during startup, is stored in the ROM
11. The computer 2 further includes a mass storage device 14 for storing an operating system 16, application programs 26, and seismic data 28, which will be described in greater detail below.
The mass storage device 14 is connected to the CPU 5 through a mass storage controller (not shown) connected to the bus 12. The mass storage device 14 and its associated computer readable media provide non-volatile storage for the computer 2. Although the description of computer readable media contained herein refers to a mass storage device, such as a hard disk or CD-ROM drive, it should be appreciated by those skilled in the art that computer readable media can be any available media that can be accessed by the computer 2.
By way of example, and not limitation, computer readable media may comprise computer storage media and communication media. Computer storage media includes volatile and non-volatile, removable and non-removable media implemented in any method or technology for storage of information such as computer readable instructions, data structures, program modules or other data.
Computer storage media includes, but is not limited to, RAM, ROM, EPROM, EEPROM, flash memory or other solid state memory technology, CD-ROM, digital versatile disks ("DVD"), or other optical storage, rnagnetic cassettes, magnetic tape, magnetic disk storage or other magnetic storage devices, or any other medium which can be used to store the desired information and which can be accessed by the computer 2.
The computer 2 may also include an input/output controller 22 for receiving and processing input from a number of other devices, including a keyboard, mouse, or electronic stylus (not shown in FIGURE 1). Similarly, an input/output controller 22 may provide output to display screen 24, a printer, or other type of output device.
As mentioned briefly above, a number of program modules and data files may be stored in the mass storage device 14 and RAM 9 of the computer 2, including an operating system 16 suitable for controlling the operation of a personal computer. The computer 2 is also capable of executing one or more application programs. In particular, the computer 2 is operative to execute casing and tubing design application program 26. According to the various illustrative embodiments of the invention, the casing and tubing design application program 26 (hereinafter referred to as "the application program 26") comprises program modules for performing various "buclcling" calculations used in the design of oil well casing and tubing. The data files stored in the mass storage device 14 may include well parameter data 28. The well parameter data 28 may include, but is not limited to, well tubing size (e.g., the inside and outside dimensions of the well tubing), tubing weight, well depth, well geometry (e.g., whether a well is vertical, horizontal, or otherwise deviated), radial clearance (i.e., the maximum distance tubing may move from the center of the wellbore or casing until it touches the wall of the wellbore or casing that it is confined by), the moment of inertia for the tubing, the temperature of the tubing in a wellbore, the current pressure in the wellbore, and whether the wellbore contains a packer or centralizer. As is known to those skilled in the art, packers are devices for holding tubing in a wellbore when the tubing is run from the surface. Packers provide a pressure seal for the wellbore and prevent fluids from mixing down hole. Centralizers are mechanical devices (i.e., collars) which are used to position casing concentrically in a wellbore and prevent the casing from lying eccentrically against the wellbore wall.
As will be described in greater detail below, the well parameter data is utilized by the application program 26 to perform buckling calculations for designing oil well casing and tubing. According to one embodiment of the invention, the application program 26 may comprise the WELLCAT application program marketed by LANDMARK GRAPHICS CORPORATION of Houston, Texas. It should be appreciated, however, that the various aspects of the invention described herein may be utilized with other application programs from other manufacturers. Additional details regarding the various calculations performed by the application program 26 will be provided below with respect to FIGURES 2-5.
Referring now to FIGURES 2-5, illustrative logical operations or routines will be described illustrating a process for determining design parameters for oil well casing and tubing to prevent buckling in a deviated wellbore. When reading the discussion of the illustrative routines presented herein, it should be appreciated that the logical operations of various embodiments of the present invention are implemented (1) as a sequence of computer implemented acts or program modules running on a computing system and/or (2) as interconnected machine logic circuits or circuit modules within the computing system. The implementation is a matter of choice dependent on the performance requirements of the computing system implementing the invention. Accordingly, the logical operations illustrated in FIGURES 2-5, and making up illustrative embodiments of the present invention described herein are referred to variously as operations, structural devices, acts or modules. It will be recognized by one skilled in the art that these operations, structural devices, acts and modules may be implemented in software, in firmware, in special purpose digital logic, and any combination thereof without deviating from the spirit and scope of the present invention as recited within the claims attached hereto.
In the following discussion of FIGURES 2-5, a number of formulae utilized by the application program 26 to calculate parameters for predicting various buckling conditions will be described using the following nomenclature:
E = Young's modulus, psi F = axial buckling force P = buckling force, lbf G= pipe shear modulus p = the pitch of a helix, L, ft.
I = moment of inertia of tubing, L4, in4 J= polar moment of inertia of tubing, L4, in4 El = the bending stiffiiess of tubing M = total bending moment, ft-lbf.
M; = bending moment in i direction, ft-lbf.
rc = tubing-casing radial clearance, L, in.
rp = tubing-casing radius, L, in.
do = tubing outside diameter, L, in.
s = measured depth, L, ft.
wc = contact load between a wellbore and tubing wbp = the buoyant weight of the tubing n, = the vertical component of the normal to the wellbore trajectory bZ = the vertical component to the binormal to the wellbore trajectory x = wellbore curvature T = term in contact force equation, dimensionless ul, u2 = tubing displacements, L, in.
wõ = the contact load between the tubing and casing, lbf/ft.
a= coefficient in solutions, L-1, ft"1 (3 = coefficient in solutions, L'1, ft-1 S, = parameters in beam-column equations ( is also the dynamic coefficient of friction in buckling criterion with friction equations) Oso, As1= beam-column solution lengths, L, ft.
E, so, s1 slopes in beam-column solutions, dimensionless 0= angle between the pipe center location and an x coordinate 01= angle in beam-colunui solution, radians ~=dimensionless length =as subscript o indicates initial conditions Referring now to FIGURE 2, an illustrative routine 200 performed by a processing device, such as the CPU 5 of the computer of FIG. 1 will be described for calculating a parameter for predicting the movement of tubing near at least one boundary condition in a deviated wellbore, according to one embodiment of the invention. As defined herein and in the appended claims, a "boundary condition"
Illustrative embodiments of the present invention provide for determining design parameters for oil well casing and tubing to prevent buckling in a deviated wellbore. Referring now to the drawings, in which like numerals represent like elements, various aspects of the present invention will be described. In particular, FIGURE 1 and the corresponding discussion are intended to provide a brief, general description of a suitable computing environment in which embodiments of the invention may be iinplemented. While the invention will be described in the general context of program modules that execute in conjunction with program modules that run on an operating system on a personal computer, those skilled in the art will recognize that the invention may also be implemented in combination with other types of computer systems and program modules.
Generally, program modules include routines, programs, components, data structures, and other types of structures that perform particular tasks or implement particular abstract data types. Moreover, those skilled in the art will appreciate that the invention may be practiced with other computer systein configurations, including hand-held devices, multiprocessor systems, inicroprocessor-based or programmable consumer electronics, minicomputers, mainframe computers, and the like. The invention may also be practiced in distributed computing environments where tasks are performed by remote processing devices that are linked through a communications network. In a distributed computing environment, program modules may be located in both local and remote memory storage devices.
Referring now to FIGURE 1, an illustrative computer architecture for a computer 2 utilized in the various einbodiments of the invention will be described.
The computer architecture shown in FIGURE 1 illustrates a conventional desktop or laptop computer, including a central processing unit 5 ("CPU"), a system memory 7, including a random access memory 9 ("RAM") and a read-only memory ("ROM") 11, and a system bus 12 that couples the memory to the CPU 5. A basic input/output system containing the basic routines that help to transfer information between elements within the computer, such as during startup, is stored in the ROM
11. The computer 2 further includes a mass storage device 14 for storing an operating system 16, application programs 26, and seismic data 28, which will be described in greater detail below.
The mass storage device 14 is connected to the CPU 5 through a mass storage controller (not shown) connected to the bus 12. The mass storage device 14 and its associated computer readable media provide non-volatile storage for the computer 2. Although the description of computer readable media contained herein refers to a mass storage device, such as a hard disk or CD-ROM drive, it should be appreciated by those skilled in the art that computer readable media can be any available media that can be accessed by the computer 2.
By way of example, and not limitation, computer readable media may comprise computer storage media and communication media. Computer storage media includes volatile and non-volatile, removable and non-removable media implemented in any method or technology for storage of information such as computer readable instructions, data structures, program modules or other data.
Computer storage media includes, but is not limited to, RAM, ROM, EPROM, EEPROM, flash memory or other solid state memory technology, CD-ROM, digital versatile disks ("DVD"), or other optical storage, rnagnetic cassettes, magnetic tape, magnetic disk storage or other magnetic storage devices, or any other medium which can be used to store the desired information and which can be accessed by the computer 2.
The computer 2 may also include an input/output controller 22 for receiving and processing input from a number of other devices, including a keyboard, mouse, or electronic stylus (not shown in FIGURE 1). Similarly, an input/output controller 22 may provide output to display screen 24, a printer, or other type of output device.
As mentioned briefly above, a number of program modules and data files may be stored in the mass storage device 14 and RAM 9 of the computer 2, including an operating system 16 suitable for controlling the operation of a personal computer. The computer 2 is also capable of executing one or more application programs. In particular, the computer 2 is operative to execute casing and tubing design application program 26. According to the various illustrative embodiments of the invention, the casing and tubing design application program 26 (hereinafter referred to as "the application program 26") comprises program modules for performing various "buclcling" calculations used in the design of oil well casing and tubing. The data files stored in the mass storage device 14 may include well parameter data 28. The well parameter data 28 may include, but is not limited to, well tubing size (e.g., the inside and outside dimensions of the well tubing), tubing weight, well depth, well geometry (e.g., whether a well is vertical, horizontal, or otherwise deviated), radial clearance (i.e., the maximum distance tubing may move from the center of the wellbore or casing until it touches the wall of the wellbore or casing that it is confined by), the moment of inertia for the tubing, the temperature of the tubing in a wellbore, the current pressure in the wellbore, and whether the wellbore contains a packer or centralizer. As is known to those skilled in the art, packers are devices for holding tubing in a wellbore when the tubing is run from the surface. Packers provide a pressure seal for the wellbore and prevent fluids from mixing down hole. Centralizers are mechanical devices (i.e., collars) which are used to position casing concentrically in a wellbore and prevent the casing from lying eccentrically against the wellbore wall.
As will be described in greater detail below, the well parameter data is utilized by the application program 26 to perform buckling calculations for designing oil well casing and tubing. According to one embodiment of the invention, the application program 26 may comprise the WELLCAT application program marketed by LANDMARK GRAPHICS CORPORATION of Houston, Texas. It should be appreciated, however, that the various aspects of the invention described herein may be utilized with other application programs from other manufacturers. Additional details regarding the various calculations performed by the application program 26 will be provided below with respect to FIGURES 2-5.
Referring now to FIGURES 2-5, illustrative logical operations or routines will be described illustrating a process for determining design parameters for oil well casing and tubing to prevent buckling in a deviated wellbore. When reading the discussion of the illustrative routines presented herein, it should be appreciated that the logical operations of various embodiments of the present invention are implemented (1) as a sequence of computer implemented acts or program modules running on a computing system and/or (2) as interconnected machine logic circuits or circuit modules within the computing system. The implementation is a matter of choice dependent on the performance requirements of the computing system implementing the invention. Accordingly, the logical operations illustrated in FIGURES 2-5, and making up illustrative embodiments of the present invention described herein are referred to variously as operations, structural devices, acts or modules. It will be recognized by one skilled in the art that these operations, structural devices, acts and modules may be implemented in software, in firmware, in special purpose digital logic, and any combination thereof without deviating from the spirit and scope of the present invention as recited within the claims attached hereto.
In the following discussion of FIGURES 2-5, a number of formulae utilized by the application program 26 to calculate parameters for predicting various buckling conditions will be described using the following nomenclature:
E = Young's modulus, psi F = axial buckling force P = buckling force, lbf G= pipe shear modulus p = the pitch of a helix, L, ft.
I = moment of inertia of tubing, L4, in4 J= polar moment of inertia of tubing, L4, in4 El = the bending stiffiiess of tubing M = total bending moment, ft-lbf.
M; = bending moment in i direction, ft-lbf.
rc = tubing-casing radial clearance, L, in.
rp = tubing-casing radius, L, in.
do = tubing outside diameter, L, in.
s = measured depth, L, ft.
wc = contact load between a wellbore and tubing wbp = the buoyant weight of the tubing n, = the vertical component of the normal to the wellbore trajectory bZ = the vertical component to the binormal to the wellbore trajectory x = wellbore curvature T = term in contact force equation, dimensionless ul, u2 = tubing displacements, L, in.
wõ = the contact load between the tubing and casing, lbf/ft.
a= coefficient in solutions, L-1, ft"1 (3 = coefficient in solutions, L'1, ft-1 S, = parameters in beam-column equations ( is also the dynamic coefficient of friction in buckling criterion with friction equations) Oso, As1= beam-column solution lengths, L, ft.
E, so, s1 slopes in beam-column solutions, dimensionless 0= angle between the pipe center location and an x coordinate 01= angle in beam-colunui solution, radians ~=dimensionless length =as subscript o indicates initial conditions Referring now to FIGURE 2, an illustrative routine 200 performed by a processing device, such as the CPU 5 of the computer of FIG. 1 will be described for calculating a parameter for predicting the movement of tubing near at least one boundary condition in a deviated wellbore, according to one embodiment of the invention. As defined herein and in the appended claims, a "boundary condition"
may comprise either a packer or a centralizer installed in a deviated wellbore. The routine 200 begins at operation 210 where the application program 26 receives the well parameter data 28 by retrieving it from the mass storage device 14. As discussed above with respect to FIGURE 1, the well parameter data 28 may include a number of measurements including well tubing size (e.g., the inside and outside dimensions of the well tubing), tubing weight, well depth, well geometry (e.g., whether a well is vertical, horizontal, or 6therwise deviated), radial clearance (i.e., the maximum distance tubing may move from the center of the wellbore or casing until it touches the wall of the wellbore or casing that it is confined by), the moment of inertia for the tubing, the temperature of the tubing in a wellbore, the current pressure in the wellbore, and whether the wellbore contains a packer or centralizer.
It will be appreciated that the well parameter data 28 may also be manually inputted directly into the application program 26 by a user.
The routine 200 then continues from operation 210 at operation 220 where the application program 26 calculates a paraineter for predicting the movement (i.e., displacement) of tubing near a packer in the deviated wellbore when the tubing starts to buckle. In particular, the application program 26 calculates a "beam-column"
solution. As is known to those skilled in the art, a beam-column is a structural member that is subjected to simultaneous axial and transverse loads (i.e., compression and bending). For the packer boundary condition, the application program 26 performs an analysis to calculate a beam-column solution to buckling equations which brings the tubing from a centralized position, tangent to the wellbore, to a point tangent to the wellbore wall. The application program 26 utilizes the following equations to satisfy these conditions:
u1b =[sinsin+(1-cos~o)(cosS
U2b =E[(1-cos~o)(sin4 -4)+(sin4o -4o)(cos4 -1)]/8 S = 40 sin~a -2(1-cos4o) where ~_ -s E is given by:
/cos0-1 and ~o is approximately 3.84333.
The application program 26 then calculates a solution d0/d~ for the above equations which is:
d0 _ F tanh( ~
A~ +
a~ c~ 5 ) where: ~s - 1.01108. It will be appreciated that the above solution equation may be integrated to give theta:
e(0 =1n cosh(-vry%A~ + ~ S ) cosh(~s) The routine 200 then continues from operation 220 at operation 230 where the application program 26 calculates a parameter for predicting the movement (i.e., displacement) of tubing near a centralizer in the deviated wellbore when the tubing starts to buckle. For the centralizer boundary condition, the application program 26 perfonns an analysis to calculate a beam-column solution to buckling equations which brings the tubing from a centralized position, free to rotate, to a point tangent to the wellbore wall. The application program 26 utilizes the following equations to satisfy these conditions:
Ulb =[(4 -sin~)+(cos~o -1)4]/
U 2b =E [4o(sin4 -4)+ (4a -sin~o)4]/
=4o cos4 o -sin4o 41~1_ 4_ -s 25 where s is given by:
/_sin 4o and 4o is approximately 2.505309.
The application program 26 calculates a solution d6/d4 for the above equations which is:
It will be appreciated that the well parameter data 28 may also be manually inputted directly into the application program 26 by a user.
The routine 200 then continues from operation 210 at operation 220 where the application program 26 calculates a paraineter for predicting the movement (i.e., displacement) of tubing near a packer in the deviated wellbore when the tubing starts to buckle. In particular, the application program 26 calculates a "beam-column"
solution. As is known to those skilled in the art, a beam-column is a structural member that is subjected to simultaneous axial and transverse loads (i.e., compression and bending). For the packer boundary condition, the application program 26 performs an analysis to calculate a beam-column solution to buckling equations which brings the tubing from a centralized position, tangent to the wellbore, to a point tangent to the wellbore wall. The application program 26 utilizes the following equations to satisfy these conditions:
u1b =[sinsin+(1-cos~o)(cosS
U2b =E[(1-cos~o)(sin4 -4)+(sin4o -4o)(cos4 -1)]/8 S = 40 sin~a -2(1-cos4o) where ~_ -s E is given by:
/cos0-1 and ~o is approximately 3.84333.
The application program 26 then calculates a solution d0/d~ for the above equations which is:
d0 _ F tanh( ~
A~ +
a~ c~ 5 ) where: ~s - 1.01108. It will be appreciated that the above solution equation may be integrated to give theta:
e(0 =1n cosh(-vry%A~ + ~ S ) cosh(~s) The routine 200 then continues from operation 220 at operation 230 where the application program 26 calculates a parameter for predicting the movement (i.e., displacement) of tubing near a centralizer in the deviated wellbore when the tubing starts to buckle. For the centralizer boundary condition, the application program 26 perfonns an analysis to calculate a beam-column solution to buckling equations which brings the tubing from a centralized position, free to rotate, to a point tangent to the wellbore wall. The application program 26 utilizes the following equations to satisfy these conditions:
Ulb =[(4 -sin~)+(cos~o -1)4]/
U 2b =E [4o(sin4 -4)+ (4a -sin~o)4]/
=4o cos4 o -sin4o 41~1_ 4_ -s 25 where s is given by:
/_sin 4o and 4o is approximately 2.505309.
The application program 26 calculates a solution d6/d4 for the above equations which is:
d0 F2 tanh( ~
A~ +
a~
where .81965. It will be appreciated that the above solution equation may be integrated to give theta:
cosh(~/A4 +
0(~) =1n cosh(~,,) It will be appreciated by those skilled in the art that the buckling calculations discussed above apply to "near" boundary conditions in a wellbore, contrary to previous buckling models which only applied to "far away" from the boundary conditions.
The routine 200 then continues from operation 230 at operation 240 where the application program 26 generates an output table of the results of the calculations performed in operations 220 and 230. In particular, the results may comprise a table of solutions corresponding to various sizes and weights of tubing, well depths, and axial forces at various well depths. The routine 200 then ends.
Referring now to FIGURE 3, an illustrative routine 300 performed by a processing device, such as the CPU 5 of the computer of FIG. 1 will be described for calculating a parameter for predicting total bending moments and maximum bending stresses near a boundary condition in a deviated wellbore, according to one embodiment of the invention. The routine 300 begins at operation 310 where the application program 26 receives the well .paraineter data 28.
The routine 300 then continues from operation 310 at operation 320 where the application program 26 calculates a parameter for predicting the total bending moment of tubing near a packer and/or centralizer for a beam-column solution by utilizing the following equations:
The bending stresses in the tubing are given by:
z a M =Elyd u' =Frd u' i=12 ' ds2 d42 ~
The total bending moment is therefore calculated as:
M=Fr /(d2ul)2+(dzuz)z ~
It should be understood that in the above equations, r is the radial clearance of the tubing in the packer or centralizer and ul and u2 are measures of the lateral displacement of the tubing in the deviated wellbore.
The routine 300 then continues from operation 320 at operation 330 where the application program 26 calculates a parameter for predicting the total bending moment of tubing near a packer and/or centralizer for a full contact solution (i.e., tubing in contact with the wellbore wall) by utilizing the following equation:
z M=Fr f(d4 +(~z~)z The routine 300 then continues from operation 330 at operation 340 where the application program 26 calculates a parameter for predicting the maximum bending stress for tubing near a packer and/or centralizer by utilizing the following equation:
Md o 6b = 21 It will be appreciated by those skilled in the art that, contrary to previous buckling models, the beam-colunm bending moment may exceed the full contact bending moment in both the packer and the centralizer.
The routine 300 then continues from operation 340 at operation 350 where the application program 26 generates an output table of the results of the calculations performed in operations 320 through 340. In particular, the results may coinprise a table of solutions corresponding to various sizes and weights of tubing, well depths, and axial forces at various well depths. The routine 300 then ends.
Referring now to FIGURE 4, an illustrative routine 400 performed by a processing device, such as the CPU 5 of the computer of FIG. 1 will be described for shows logical operations performed by an illustrative embodiment for calculating parameters for predicting minimum axial forces necessary to initiate buckling due to friction in a deviated wellbore. The routine 400 begins at operation 410 where the application program 26 receives the well parameter data 28.
The routine 400 then continues from operation 410 at operation 420 where the application program 26 calculates a parameter for predicting the minimum axial force to initiate buckling when tubing is rolling in a deviated well. In particular, cylindrical tubing lying on the bottom of a deviated well may be subject to rolling friction. The friction gradually produces a lateral force and a moment that is proportional to the lateral displacement of the tubing. In order to account for rolling friction, the application program 26 calculates a critical buckling parameter F
representing the minimum axial force necessary to allow buckling using the equation:
F _ GJ + 4EIw~
rP ~'~
where wc is given by the equation:
Wc = j(wbPnZ - Fx)2 + (wP b, )2 It should be understood that in cases where the tubing is laying on a flat plane, such as a seabed, the minimum axial force equation reduces to:
F=GJ
r2 P
The routine 400 then continues from operation 420 at operation 430 where the application program 26 calculates a parameter for predicting the minimum axial force to initiate buckling when tubing is rotating in a deviated well. In particular, when tubing is rotating the friction force is constant in the lateral direction relative to the tubing. In order to account for friction caused by rotation, the application program calculates the minimum axial force using the equation:
/4EIw~
F=
~~
where the contact load wc is given by the equation:
_ (Wbpn~ -FK) 2 +(wbpbz )2 w~ ~ 1 + ,u 2 The routine 400 then continues from operation 430 at operation 440 where the application program 26 generates an output table of the results of the calculations performed in operations 420 and 430. In particular, the results may comprise a table of solutions corresponding to various sizes and weights of tubing. The routine then ends.
Referring now to FIGURE 5, an illustrative routine 500 performed by a processing device, such as the CPU 5 of the computer of FIG. 1 will be described for shows logical operations performed by an illustrative embodiment for calculating a parameter for predicting the onset of buckling for the connection of tubing of different sizes (i.e., tapered strings) in a deviated wellbore. The routine 500 begins at operation 410 where the application program 26 receives the well parameter data 28.
The routine 500 then continues from operation 510 at operation 520 where the application program 26 calculates a parameter for predicting the onset of buckling for tapered strings by utilizing the following equations:
V1(s) = ri - 2 (r; ri )[abs - sin(abs)l vz (s) - - r'2', [ 1-cos(ocbs)1 ai s E(0, 27c ab = W
O
Gb Where the subscript b refers to the properties of the beam-column. The " "
term means that the beam-colunm solution can move either to the 0 = 0 (+ solution) or to the 0=7r (- solution). This means that the beam column solution can create either a right hand or left hand helix, depending on which way the solution moves. Assuming that the ith solution satisfies the above equations, the application program 26 further utilizes the following equations:
6l (s) -ln[sech( 2 a1s)]
de' = 2 a, tanh( Z
ds a~s) _ P
at EZIZ
Finally, the application program calculates a solution to the following differential equation for tubing in contact with the wellbore wall, provided r; is less than rj:
dO(s) -52 aj r; sd(Xs - ab ~,, k) ds 2rj 1+E+(1-E)sd2(Xs- b ~,,k) r? -r.2 1-~
' 2', k ,X = 2 aj 1+E
1+E
ri P
ai _ where sd(*,k) is a Jacobi elliptic function with parameter k.
It will be appreciated by those skilled in the art that the above buckling calculations account for tubing with different radial clearances and bending stiffness contrary to previous buckling models which only applied to tubing sections of the same size (i.e., they did not apply to tapered strings).
The routine 500 then continues from operation 520 at operation 530 where the application program 26 generates an output table of the results of the calculations performed in operation 520. In particular, the results may comprise a table of solutions corresponding to various sizes and weights of tubing, well depths, and axial forces at various well depths. The routine 500 then ends.
Based on the foregoing, it should be appreciated that the various embodiments of the invention include methods and computer readable media for determining design parameters for oil well casing and tubing to prevent buckling in deviated wellbores. Although the present invention has been described in connection with various illustrative embodiments, those of ordinary skill in the art will understand that many modifications can be made thereto within the scope of the claims that follow. Accordingly, it is not intended that the scope of the invention in any way be limited by the above description, but instead be determined entirely by reference to the claims that follow.
A~ +
a~
where .81965. It will be appreciated that the above solution equation may be integrated to give theta:
cosh(~/A4 +
0(~) =1n cosh(~,,) It will be appreciated by those skilled in the art that the buckling calculations discussed above apply to "near" boundary conditions in a wellbore, contrary to previous buckling models which only applied to "far away" from the boundary conditions.
The routine 200 then continues from operation 230 at operation 240 where the application program 26 generates an output table of the results of the calculations performed in operations 220 and 230. In particular, the results may comprise a table of solutions corresponding to various sizes and weights of tubing, well depths, and axial forces at various well depths. The routine 200 then ends.
Referring now to FIGURE 3, an illustrative routine 300 performed by a processing device, such as the CPU 5 of the computer of FIG. 1 will be described for calculating a parameter for predicting total bending moments and maximum bending stresses near a boundary condition in a deviated wellbore, according to one embodiment of the invention. The routine 300 begins at operation 310 where the application program 26 receives the well .paraineter data 28.
The routine 300 then continues from operation 310 at operation 320 where the application program 26 calculates a parameter for predicting the total bending moment of tubing near a packer and/or centralizer for a beam-column solution by utilizing the following equations:
The bending stresses in the tubing are given by:
z a M =Elyd u' =Frd u' i=12 ' ds2 d42 ~
The total bending moment is therefore calculated as:
M=Fr /(d2ul)2+(dzuz)z ~
It should be understood that in the above equations, r is the radial clearance of the tubing in the packer or centralizer and ul and u2 are measures of the lateral displacement of the tubing in the deviated wellbore.
The routine 300 then continues from operation 320 at operation 330 where the application program 26 calculates a parameter for predicting the total bending moment of tubing near a packer and/or centralizer for a full contact solution (i.e., tubing in contact with the wellbore wall) by utilizing the following equation:
z M=Fr f(d4 +(~z~)z The routine 300 then continues from operation 330 at operation 340 where the application program 26 calculates a parameter for predicting the maximum bending stress for tubing near a packer and/or centralizer by utilizing the following equation:
Md o 6b = 21 It will be appreciated by those skilled in the art that, contrary to previous buckling models, the beam-colunm bending moment may exceed the full contact bending moment in both the packer and the centralizer.
The routine 300 then continues from operation 340 at operation 350 where the application program 26 generates an output table of the results of the calculations performed in operations 320 through 340. In particular, the results may coinprise a table of solutions corresponding to various sizes and weights of tubing, well depths, and axial forces at various well depths. The routine 300 then ends.
Referring now to FIGURE 4, an illustrative routine 400 performed by a processing device, such as the CPU 5 of the computer of FIG. 1 will be described for shows logical operations performed by an illustrative embodiment for calculating parameters for predicting minimum axial forces necessary to initiate buckling due to friction in a deviated wellbore. The routine 400 begins at operation 410 where the application program 26 receives the well parameter data 28.
The routine 400 then continues from operation 410 at operation 420 where the application program 26 calculates a parameter for predicting the minimum axial force to initiate buckling when tubing is rolling in a deviated well. In particular, cylindrical tubing lying on the bottom of a deviated well may be subject to rolling friction. The friction gradually produces a lateral force and a moment that is proportional to the lateral displacement of the tubing. In order to account for rolling friction, the application program 26 calculates a critical buckling parameter F
representing the minimum axial force necessary to allow buckling using the equation:
F _ GJ + 4EIw~
rP ~'~
where wc is given by the equation:
Wc = j(wbPnZ - Fx)2 + (wP b, )2 It should be understood that in cases where the tubing is laying on a flat plane, such as a seabed, the minimum axial force equation reduces to:
F=GJ
r2 P
The routine 400 then continues from operation 420 at operation 430 where the application program 26 calculates a parameter for predicting the minimum axial force to initiate buckling when tubing is rotating in a deviated well. In particular, when tubing is rotating the friction force is constant in the lateral direction relative to the tubing. In order to account for friction caused by rotation, the application program calculates the minimum axial force using the equation:
/4EIw~
F=
~~
where the contact load wc is given by the equation:
_ (Wbpn~ -FK) 2 +(wbpbz )2 w~ ~ 1 + ,u 2 The routine 400 then continues from operation 430 at operation 440 where the application program 26 generates an output table of the results of the calculations performed in operations 420 and 430. In particular, the results may comprise a table of solutions corresponding to various sizes and weights of tubing. The routine then ends.
Referring now to FIGURE 5, an illustrative routine 500 performed by a processing device, such as the CPU 5 of the computer of FIG. 1 will be described for shows logical operations performed by an illustrative embodiment for calculating a parameter for predicting the onset of buckling for the connection of tubing of different sizes (i.e., tapered strings) in a deviated wellbore. The routine 500 begins at operation 410 where the application program 26 receives the well parameter data 28.
The routine 500 then continues from operation 510 at operation 520 where the application program 26 calculates a parameter for predicting the onset of buckling for tapered strings by utilizing the following equations:
V1(s) = ri - 2 (r; ri )[abs - sin(abs)l vz (s) - - r'2', [ 1-cos(ocbs)1 ai s E(0, 27c ab = W
O
Gb Where the subscript b refers to the properties of the beam-column. The " "
term means that the beam-colunm solution can move either to the 0 = 0 (+ solution) or to the 0=7r (- solution). This means that the beam column solution can create either a right hand or left hand helix, depending on which way the solution moves. Assuming that the ith solution satisfies the above equations, the application program 26 further utilizes the following equations:
6l (s) -ln[sech( 2 a1s)]
de' = 2 a, tanh( Z
ds a~s) _ P
at EZIZ
Finally, the application program calculates a solution to the following differential equation for tubing in contact with the wellbore wall, provided r; is less than rj:
dO(s) -52 aj r; sd(Xs - ab ~,, k) ds 2rj 1+E+(1-E)sd2(Xs- b ~,,k) r? -r.2 1-~
' 2', k ,X = 2 aj 1+E
1+E
ri P
ai _ where sd(*,k) is a Jacobi elliptic function with parameter k.
It will be appreciated by those skilled in the art that the above buckling calculations account for tubing with different radial clearances and bending stiffness contrary to previous buckling models which only applied to tubing sections of the same size (i.e., they did not apply to tapered strings).
The routine 500 then continues from operation 520 at operation 530 where the application program 26 generates an output table of the results of the calculations performed in operation 520. In particular, the results may comprise a table of solutions corresponding to various sizes and weights of tubing, well depths, and axial forces at various well depths. The routine 500 then ends.
Based on the foregoing, it should be appreciated that the various embodiments of the invention include methods and computer readable media for determining design parameters for oil well casing and tubing to prevent buckling in deviated wellbores. Although the present invention has been described in connection with various illustrative embodiments, those of ordinary skill in the art will understand that many modifications can be made thereto within the scope of the claims that follow. Accordingly, it is not intended that the scope of the invention in any way be limited by the above description, but instead be determined entirely by reference to the claims that follow.
Claims (10)
1. A method of determining design parameters for oil well casing and tubing to prevent buckling in a deviated wellbore, comprising:
receiving well parameter data comprising at least one of tubing size, tubing weight, well depth, and well geometry;
calculating a first parameter for predicting the movement of tubing near at least one boundary condition in the deviated wellbore based on the received well parameter data;
calculating a second parameter for predicting a total bending moment near the at least one boundary condition based on the received well parameter data;
calculating a third parameter for predicting a maximum bending stress near the at least one boundary condition in the deviated wellbore based on the total bending moment; and calculating a fourth parameter for predicting the minimum axial force necessary to initiate buckling based on the received well parameter data, wherein the first, second, third, and fourth parameters are utilized in the design of the oil well casing and tubing to prevent buckling in the deviated wellbore.
receiving well parameter data comprising at least one of tubing size, tubing weight, well depth, and well geometry;
calculating a first parameter for predicting the movement of tubing near at least one boundary condition in the deviated wellbore based on the received well parameter data;
calculating a second parameter for predicting a total bending moment near the at least one boundary condition based on the received well parameter data;
calculating a third parameter for predicting a maximum bending stress near the at least one boundary condition in the deviated wellbore based on the total bending moment; and calculating a fourth parameter for predicting the minimum axial force necessary to initiate buckling based on the received well parameter data, wherein the first, second, third, and fourth parameters are utilized in the design of the oil well casing and tubing to prevent buckling in the deviated wellbore.
2. The method of claim 1 further comprising:
calculating a fifth parameter for predicting the onset of buckling for the connection of tubing of different sizes based on the received well parameter data, wherein the fifth parameter is utilized in the design of the oil well casing and tubing to prevent buckling in the deviated wellbore.
calculating a fifth parameter for predicting the onset of buckling for the connection of tubing of different sizes based on the received well parameter data, wherein the fifth parameter is utilized in the design of the oil well casing and tubing to prevent buckling in the deviated wellbore.
3. The method of claim 1, wherein calculating a first parameter for predicting the movement of tubing near at least one boundary condition in the deviated wellbore based on the received well parameter data comprises calculating a parameter for predicting the movement of tubing near a packer in the deviated wellbore using the formula:
where .theta.(.xi.) is a buckling parameter for a beam-column solution for tubing located near the packer in the deviated wellbore;
.DELTA..xi. is the change in dimensionless length associated with the tubing where is given by the relationship:
where s is the measured depth of the tubing;
P is the axial buckling force of the tubing; and EI is the bending stiffness of the tubing; and .PHI.s is a numerical constant.
where .theta.(.xi.) is a buckling parameter for a beam-column solution for tubing located near the packer in the deviated wellbore;
.DELTA..xi. is the change in dimensionless length associated with the tubing where is given by the relationship:
where s is the measured depth of the tubing;
P is the axial buckling force of the tubing; and EI is the bending stiffness of the tubing; and .PHI.s is a numerical constant.
4. The method of claim 1, wherein calculating a first parameter for predicting the movement of tubing near at least one boundary condition in the deviated wellbore based on the received well parameter data comprises calculating a parameter for predicting the movement of tubing near a centralizer in the deviated wellbore using the formula:
where .theta.(.xi.) is a buckling parameter for a beam-column solution for tubing located near the centralizer in the deviated wellbore;
.DELTA..xi. is the change in dimensionless length associated with the tubing where .xi. is given by the relationship:
where s is the measured depth of the tubing;
P is the axial buckling force of the tubing; and EI is the bending stiffness of the tubing; and .phi.c is a numerical constant.
where .theta.(.xi.) is a buckling parameter for a beam-column solution for tubing located near the centralizer in the deviated wellbore;
.DELTA..xi. is the change in dimensionless length associated with the tubing where .xi. is given by the relationship:
where s is the measured depth of the tubing;
P is the axial buckling force of the tubing; and EI is the bending stiffness of the tubing; and .phi.c is a numerical constant.
5. The method of claim 1, wherein calculating a second parameter for predicting a total bending moment near at least one boundary condition in the deviated wellbore based on the received well parameter data comprises calculating a parameter for predicting a total bending moment of the tubing near at least one of a packer and a centralizer in the deviated wellbore using the formula:
where M is the total bending moment in a beam-column solution for the packer or centralizer in the deviated wellbore;
F is the bending stiffness of the tubing;
r is the radial clearance of the tubing in the packer or centralizer; and u1 and u2 are measures of the lateral displacement of the tubing in the deviated wellbore.
where M is the total bending moment in a beam-column solution for the packer or centralizer in the deviated wellbore;
F is the bending stiffness of the tubing;
r is the radial clearance of the tubing in the packer or centralizer; and u1 and u2 are measures of the lateral displacement of the tubing in the deviated wellbore.
6. The method of claim 1, wherein calculating a second parameter for predicting a total bending moment near at least one boundary condition in the deviated wellbore based on the received well parameter data comprises calculating a parameter for predicting a total bending moment near at least one of a packer and a centralizer in the deviated wellbore using the formula:
where M is the total bending moment in a full contact solution for the packer or centralizer in the deviated wellbore;
F is the axial buckling force of the tubing;
r is the radial clearance of the tubing in the packer or centralizer; and .theta. is the angle between a tubing center location and an x coordinate on a coordinate axis from the tubing center location to a point tangent to the wall of the deviated wellbore, wherein x = d.theta./d.xi..
where M is the total bending moment in a full contact solution for the packer or centralizer in the deviated wellbore;
F is the axial buckling force of the tubing;
r is the radial clearance of the tubing in the packer or centralizer; and .theta. is the angle between a tubing center location and an x coordinate on a coordinate axis from the tubing center location to a point tangent to the wall of the deviated wellbore, wherein x = d.theta./d.xi..
7. The method of claim 5, wherein calculating a third parameter for predicting a maximum bending stress near the at least one boundary condition based on the total bending moment comprises calculating a parameter for predicting a maximum bending stress near the at least one of a packer and a centralizer in the deviated wellbore using the formula:
where .sigma.b is the maximum bending stress; and d o is the outside diameter of the tubing.
where .sigma.b is the maximum bending stress; and d o is the outside diameter of the tubing.
8. The method of claim 1, wherein calculating a fourth parameter for predicting the minimum axial force necessary to initiate buckling when the tubing is constrained by friction forces based on the received well parameter data comprises using the formula:
where F is the minimum axial force necessary to initiate buckling in the tubing when the tubing is rolling in the deviated wellbore;
G is the shear modulus of the tubing;
J is the polar moment of inertia of the tubing;
r p is the radius of the tubing;
EI is the bending stiffness of the tubing;
w c is the contact load between the deviated wellbore and the tubing;
and r c is the radial clearance of the tubing;
where F is the minimum axial force necessary to initiate buckling in the tubing when the tubing is rolling in the deviated wellbore;
G is the shear modulus of the tubing;
J is the polar moment of inertia of the tubing;
r p is the radius of the tubing;
EI is the bending stiffness of the tubing;
w c is the contact load between the deviated wellbore and the tubing;
and r c is the radial clearance of the tubing;
9. The method of claim 1, wherein calculating a fourth parameter for predicting the minimum axial force necessary to initiate buckling when the tubing is constrained by friction forces based on the received well parameter data comprises using the formula:
where F is the minimum axial force necessary to initiate buckling in the tubing when the tubing is rotating in the deviated wellbore;
EI is the bending stiffness of the tubing;
r p is the radius of the tubing;
r c is the radial clearance of the tubing; and w c is the contact load between the deviated wellbore and the tubing, wherein w c is given by the relationship:
where w bp is the buoyant weight of the tubing;
n z is the vertical component of the normal to the trajectory of the deviated wellbore;
b z is the vertical component to the binormal to the trajectory of the deviated wellbore;
.KAPPA. is the curvature of the deviated wellbore; and µ is the dynamic coefficient of friction with respect to the tubing in the deviated wellbore.
where F is the minimum axial force necessary to initiate buckling in the tubing when the tubing is rotating in the deviated wellbore;
EI is the bending stiffness of the tubing;
r p is the radius of the tubing;
r c is the radial clearance of the tubing; and w c is the contact load between the deviated wellbore and the tubing, wherein w c is given by the relationship:
where w bp is the buoyant weight of the tubing;
n z is the vertical component of the normal to the trajectory of the deviated wellbore;
b z is the vertical component to the binormal to the trajectory of the deviated wellbore;
.KAPPA. is the curvature of the deviated wellbore; and µ is the dynamic coefficient of friction with respect to the tubing in the deviated wellbore.
10. The method of claim 2, wherein calculating a fifth parameter for predicting the onset of buckling for the connection of tubing of different sizes based on the received well parameter data comprises using the formula:
where is a buckling parameter for a beam-column solution to predict the onset of buckling for a connection of a first tubing and a second tubing;
r i is the radial clearance of the first tubing;
r j is the radial clearance of the second tubing, wherein r i < r j;
, where P is the buckling force associated with the connection of the first tubing and the second tubing and E j I j is the bending stiffness of the second tubing;
, where the subscript b refers to the properties of the beam-column solution, where F is the axial buckling force associated with the connection of the first tubing and the second tubing, and E b l b is the bending stiffness; and where sd(*,k) is a Jacobi elliptic function with parameter k.
where is a buckling parameter for a beam-column solution to predict the onset of buckling for a connection of a first tubing and a second tubing;
r i is the radial clearance of the first tubing;
r j is the radial clearance of the second tubing, wherein r i < r j;
, where P is the buckling force associated with the connection of the first tubing and the second tubing and E j I j is the bending stiffness of the second tubing;
, where the subscript b refers to the properties of the beam-column solution, where F is the axial buckling force associated with the connection of the first tubing and the second tubing, and E b l b is the bending stiffness; and where sd(*,k) is a Jacobi elliptic function with parameter k.
Applications Claiming Priority (5)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US72351305P | 2005-10-04 | 2005-10-04 | |
| US60/723,513 | 2005-10-04 | ||
| US11/274,637 US7412368B2 (en) | 2004-11-15 | 2005-11-15 | Methods and computer-readable media for determining design parameters to prevent tubing buckling in deviated wellbores |
| US11/274,637 | 2005-11-15 | ||
| PCT/US2006/038677 WO2007041594A1 (en) | 2005-10-04 | 2006-10-03 | Methods and computer-readable media for determining design parameters to prevent tubing buckling in deviated wellbores |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CA2625178A1 true CA2625178A1 (en) | 2007-04-12 |
| CA2625178C CA2625178C (en) | 2014-06-03 |
Family
ID=37635618
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CA2625178A Expired - Fee Related CA2625178C (en) | 2005-10-04 | 2006-10-03 | Methods and computer-readable media for determining design parameters to prevent tubing buckling in deviated wellbores |
Country Status (9)
| Country | Link |
|---|---|
| US (1) | US7412368B2 (en) |
| EP (1) | EP1931856B1 (en) |
| AT (1) | ATE465323T1 (en) |
| AU (1) | AU2006299480B2 (en) |
| BR (1) | BRPI0616669A2 (en) |
| CA (1) | CA2625178C (en) |
| DE (1) | DE602006013850D1 (en) |
| NO (1) | NO340815B1 (en) |
| WO (1) | WO2007041594A1 (en) |
Families Citing this family (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2011159307A1 (en) * | 2010-06-18 | 2011-12-22 | Landmark Graphics Corporation | Systems and methods for wellbore optimization |
| CA2831056C (en) * | 2011-06-24 | 2017-08-22 | Landmark Graphics Corporation | Systems and methods for determining the moments and forces of two concentric pipes within a wellbore |
| US9043152B2 (en) | 2011-08-08 | 2015-05-26 | Baker Hughes Incorporated | Realtime dogleg severity prediction |
| US9062540B2 (en) * | 2012-05-11 | 2015-06-23 | Baker Hughes Incorporated | Misalignment compensation for deep reading azimuthal propagation resistivity |
| US9416652B2 (en) | 2013-08-08 | 2016-08-16 | Vetco Gray Inc. | Sensing magnetized portions of a wellhead system to monitor fatigue loading |
| US9845671B2 (en) | 2013-09-16 | 2017-12-19 | Baker Hughes, A Ge Company, Llc | Evaluating a condition of a downhole component of a drillstring |
| NO346971B1 (en) | 2013-09-17 | 2023-03-20 | Halliburton Energy Services Inc | Estimation and Calibration of Downhole Buckling Conditions |
| WO2016118154A1 (en) * | 2015-01-23 | 2016-07-28 | Landmark Graphics Corporation | Determining stresses in a pipe under non-uniform exterior loads |
| BR112018075026A2 (en) * | 2016-07-08 | 2019-03-06 | Halliburton Energy Services Inc. | method, machine readable storage device with instructions stored therein, and, system. |
| CN106503399B (en) * | 2016-11-19 | 2017-09-15 | 东北石油大学 | Peupendicular hole hangs the determination method of tubing string Helical Buckling Critical Load |
| GB2575597B (en) * | 2017-06-16 | 2022-03-23 | Landmark Graphics Corp | Optimized visualization of loads and resistances for wellbore tubular design |
| EP3728791A4 (en) | 2017-12-23 | 2021-09-22 | Noetic Technologies Inc. | System and method for optimizing tubular running operations using real-time measurements and modelling |
| CN109931048A (en) * | 2019-03-27 | 2019-06-25 | 南智（重庆）能源技术有限公司 | Oil/gas well tubing and casing integrality detection method and evaluation system |
| CN111177944B (en) * | 2020-01-09 | 2022-04-08 | 暨南大学 | Deep-sea pipeline buckling propagation pressure calculation method based on plate-shell theory |
| CN112883520A (en) * | 2021-03-29 | 2021-06-01 | 珠海市三鑫科技发展有限公司 | Aluminum component bending analysis method based on direct strength method |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5881311A (en) * | 1996-06-05 | 1999-03-09 | Fastor Technologies, Inc. | Data storage subsystem with block based data management |
| NO304709B1 (en) * | 1997-03-20 | 1999-02-01 | Maritime Well Service As | Device for production tubes |
| US6684952B2 (en) * | 1998-11-19 | 2004-02-03 | Schlumberger Technology Corp. | Inductively coupled method and apparatus of communicating with wellbore equipment |
| US7249044B2 (en) * | 2000-10-05 | 2007-07-24 | I2 Technologies Us, Inc. | Fulfillment management system for managing ATP data in a distributed supply chain environment |
| US6526819B2 (en) * | 2001-02-08 | 2003-03-04 | Weatherford/Lamb, Inc. | Method for analyzing a completion system |
| US8533072B2 (en) * | 2002-07-10 | 2013-09-10 | Sap Aktiengesellschaft | Multi-level glogbal available-to-promise |
| WO2004094768A2 (en) * | 2003-04-23 | 2004-11-04 | Th Hill Associates, Inc. | Drill string design methodology for mitigating fatigue failure |
-
2005
- 2005-11-15 US US11/274,637 patent/US7412368B2/en active Active
-
2006
- 2006-10-03 DE DE602006013850T patent/DE602006013850D1/en active Active
- 2006-10-03 BR BRPI0616669-5A patent/BRPI0616669A2/en not_active Application Discontinuation
- 2006-10-03 WO PCT/US2006/038677 patent/WO2007041594A1/en active Application Filing
- 2006-10-03 EP EP06825412A patent/EP1931856B1/en not_active Not-in-force
- 2006-10-03 AU AU2006299480A patent/AU2006299480B2/en not_active Ceased
- 2006-10-03 CA CA2625178A patent/CA2625178C/en not_active Expired - Fee Related
- 2006-10-03 AT AT06825412T patent/ATE465323T1/en not_active IP Right Cessation
-
2008
- 2008-04-28 NO NO20082001A patent/NO340815B1/en not_active IP Right Cessation
Also Published As
| Publication number | Publication date |
|---|---|
| DE602006013850D1 (en) | 2010-06-02 |
| US7412368B2 (en) | 2008-08-12 |
| CA2625178C (en) | 2014-06-03 |
| US20060106588A1 (en) | 2006-05-18 |
| NO20082001L (en) | 2008-06-26 |
| AU2006299480A1 (en) | 2007-04-12 |
| EP1931856B1 (en) | 2010-04-21 |
| WO2007041594A1 (en) | 2007-04-12 |
| AU2006299480B2 (en) | 2012-05-03 |
| BRPI0616669A2 (en) | 2011-06-28 |
| ATE465323T1 (en) | 2010-05-15 |
| NO340815B1 (en) | 2017-06-19 |
| EP1931856A1 (en) | 2008-06-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CA2625178A1 (en) | Methods and computer-readable media for determining design parameters to prevent tubing buckling in deviated wellbores | |
| Mitchell | Comprehensive analysis of buckling with friction | |
| Wu et al. | Coiled tubing buckling implication in drilling and completing horizontal wells | |
| Cunha | Buckling of tubulars inside wellbores: a review on recent theoretical and experimental works | |
| CN111655969B (en) | System and method for optimizing running operations of a pipe using real-time measurements and modeling | |
| Mitchell et al. | Helical buckling of pipe with connectors and torque | |
| EP2723980B1 (en) | Systems and methods for determining the moments and forces of two concentric pipes within a wellbore | |
| Juvkam-Wold et al. | Casing deflection and centralizer spacing calculations | |
| Tan et al. | A prediction model of casing wear in extended-reach drilling with buckled drillstring | |
| Huang et al. | Local-integral coupling model of tubular strings with connectors and its application in periodic sticking analyses | |
| Menand et al. | Axial force transfer of buckled drill pipe in deviated wells | |
| Mitchell | The twist and shear of helically buckled pipe | |
| US10041344B2 (en) | Determining pressure within a sealed annulus | |
| McSpadden et al. | Development of a stiff-string forces model for coiled tubing | |
| CA2913203C (en) | Methods and systems for determining manufacturing and operating parameters for a deviated downhole well component | |
| Aasen et al. | Buckling models revisited | |
| McSpadden et al. | Friction load redistribution for robust and cost efficient design of extended reach and horizontal well completions | |
| Gu et al. | Analysis of Slack-Off Force Transmitted Downhole in Coiled-Tubing Operations | |
| Zhou et al. | Analysis of leak-off tests in shallow marine sediments | |
| Yang et al. | Couette flow of non-Newtonian power-law fluids in narrow eccentric annuli | |
| Cayeux et al. | Buoyancy Force on a Plain or Perforated Portion of a Pipe | |
| Tang | Mechanical Analysis of a Multi-Test String in High-Temperature and High-Pressure Deep Wells. | |
| Xiao et al. | Spring-supported arch model for predicting hydrostatic collapse strength of flexible riser with layer gap | |
| Yang | Selection of Drillstring Components for Effective Horizontal Drilling | |
| Bybee | Axial-Force Transfer of Buckled Drillpipe in Deviated Wells |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| EEER | Examination request | ||
| MKLA | Lapsed |
Effective date: 20201005 |