Methods and systems for determining manufacturing and operating parameters for a deviated downhole well component
Systems and methods for determining manufacturing or operating parameters for a deviated downhole well component, including a method that includes representing a tubular string as nodes separated by segments, determining transfer matrices for determining an ith node's state vector from an ith−1 node's state vector, and defining initial state vector values for the reference node. The nodes are numerable from 1 to N with an initial, mechanically constrained reference node representable with i=0, and each is associated with a state vector describing a corresponding node position and one or more forces present at said node. The method further includes applying the transfer matrices to obtain each of the state vectors' values, deriving from at least one of the state vectors a parameter value for said component, and specifying a component having said parameter value. The parameter value can include a centralizer or stabilizer composition, manufacturing dimensions, or position.
Latest LANDMARK GRAPHICS CORPORATION Patents:
- Deep learning model with dilation module for fault characterization
- Determining fault surfaces from fault attribute volumes
- Event prediction using state-space mapping during drilling operations
- Operating wellbore equipment using a distributed decision framework
- Managing gas bubble migration in a downhole liquid
This application claims priority to Provisional U.S. Application Ser. No. 61/837,986, titled “Methods and Systems for Modeling a Deviated Downhole Well Component” and filed Jun. 21, 2013 by Robello Samuel and Zhenying Wang, which is incorporated herein by reference.
BACKGROUNDAs world demand for petrochemical products has continued to increase, oil and gas companies have had to expand their exploration and production efforts into developing increasingly deep wells. As a result, the structures constructed to form a well must be capable of operating under larger loads and stresses than ever before. Because failures can have costly consequences, it is important to design all well structures with appropriate safety margins.
One example of such a structure is the well casing. A well casing is a tubular structure generally made of a steel pipe surrounded by a concrete layer that secures the steel pipe to the surrounding formation, thus defining the outside wall of the well. The concrete provides support to the steel pipe, as well as additional isolation layer between the formation and fluids flowing within the casing. In order to determine the correct materials and dimensions for the various casing components, engineers frequently perform computer simulations to model various casing configurations under simulated downhole conditions. The simulations provide the engineer with information regarding the various loads and stresses to which the casing might be subjected, and enable potential designs to be evaluated.
But casing designs are only as good as the underlying simulation model. While simulations of single section casings in vertical wells are generally well understood and produce accurate results, tapered casings, deviated casings and casings with fluid flow restrictions represent indeterminate complex mechanical systems that can be very difficult or impractical to model using existing techniques. While methods do exist wherein these more complex systems are modeled as simpler single-section vertical wells with the results being adjusted to include additional safety margins, such methods can incur a significant risk, given the lack of quantifiable data to support the selected margins.
A better understanding of the various disclosed embodiments can be obtained when the following detailed description is considered in conjunction with the attached drawings, in which:
It should be understood that the drawings and corresponding detailed description do not limit the disclosure, but on the contrary, they provide the foundation for understanding all modifications, equivalents, and alternatives falling within the scope of the appended claims.
DETAILED DESCRIPTIONThe paragraphs that follow describe illustrative systems and methods for determining manufacturing or operating parameters for a deviated downhole well component. Examples are provided within the context of a tapered and deviated well casing that is mechanically constrained at opposite ends. The mechanics of such a casing are explained and illustrated, and matrices are presented that mathematically describe the known and unknown forces acting at various points along the casing. Finally, methods and systems are described that combine the various matrices to compute the forces present along the casing.
The disclosed systems and methods are best understood when described in an illustrative usage context. Accordingly,
Because boreholes are routinely drilled to ten thousand feet or more in depth and can be steered horizontally for perhaps twice that distance, a well casing string is inserted into the borehole and is cemented to the borehole wall to provide support to the borehole and isolation between the formation and the fluids flowing within the well casing string. In the example of
As can be seen in the illustrative example of
In at least some illustrative embodiments, the forces acting on the above-described well casing string are modeled by dividing the system into a series of N subsystems that each only interact with adjacent subsystems, and then determining the forces acting on each subsystem in sequence.
Vi=[uiviαiFxiFhiMi−1]iT (1)
where,
ui is the true vertical depth (TVD) of node i relative to the reference node;
vi is the horizontal distance of node i from the reference node;
αi is the inclination angle of the casing segment at node i;
Fxi is the vertical force present at node i;
Fhi is the horizontal force present at node i; and
Mi is the bending moment present at node i.
The reference node of the illustrative example is located at the casing hanger and is designated as node 0, and the node at the opposite end of the well casing string and furthest away from node 0 is designated as node N.
After dividing the illustrative well casing string into segments, the unknown forces acting on each segment's downhole node are determined by starting at the reference node at one end of the first segment (where all elements of the state vector are known) and determining the state vector for the next node at the opposite end of the segment. The state vector for the downhole node of the first segment is determined by computing a cross product of a transfer matrix and the reference node's state vector. This transfer matrix is defined based upon the known elements of the downhole node's state vector, the previous node's state vector and/or the reference node's state vector. For the first segment the previous node is also the reference node (i.e., i−1=0). In at least some illustrative embodiments, the transfer matrix Ti for a well casing string segment defined between nodes i and i−1 is described using known state vector elements as,
where,
-
- lih is the horizontal distance defined as (vi−vi−1);
- lix is the vertical distance defined as (ui−ui−1);
- Δαi−10 is the change in inclination angle at the ith−1 node defined as (αi−1−α0);
- Δαi0 is the change in inclination angle at the ith node defined as (αi−α0);
- (EI)i is the product of Young's modulus and a moment of inertia of the component at the ith node; and
- (EA)i is the product of Young's modulus and a cross-sectional area of the component at the ith node.
For the first node, i=1 and the state vector V1 is expressed as,
V1=T1×V0 (3).
Equation (3) expresses the state vector V1 as a set of constrained linear equations for the well casing segment that can be solved to determine the unknown forces present at node 1. For each subsequent node i, the cross product of the node's transfer matrix and the reference node's state vector is combined with prior cross products for nodes 1 through i−1 to determine the state vector for node i, thus determining the unknown forces at each node i (i.e., Fxi, Fhi and/or Mi). Once these forces are known, the axial force at a node can be computed as,
Fαi=Fxi×cos(αi)+Fhi×sin(αi) (4).
In at least some illustrative embodiments, a transfer matrix method (TMM) is used to perform the combination of cross products. Using TMM, the combination of cross products producing the ith node's state vector is expressed as,
Vi=(Π1iTi)×V0 (5).
Because each of the prior products of products for each node 1 through i−1 have already been calculated, the product of products does not need to be recalculated for each node. Instead, each node's transfer matrix can be combined with a cumulative transfer matrix, thus avoiding duplicative computations. For example, for a three node string this would be expressed as,
Tacc=T1; V1=Tacc×V0 (6),
Tacc=Tacc×T2; V2=Tacc×V0=(T1×T2)×V0 (7), and
Tacc=Tacc×T3; V3=Tacc×V0=(T1×T2×T3)×V0 (8).
The cross product T1×T2 has already been calculated in equation (7) and saved as cumulative transfer matrix Tacc, which is reused without recalculation in equation (8).
It should be noted that the transfer matrix is not limited to the specific embodiment of equation (2). For example, plugs such as cementing plugs present within the casing string and reductions in the cross-sectional area of the casing such as reduction 15 of
where,
-
- Fpx is a vertical force present on a plug or reduction located at the ith node; and
- Fph is a horizontal force present on the plug or reduction.
Other transfer matrices may include, for example, parameters that describe the load imposed on a casing string by a salt formation (i.e., “salt loading”). A wide variety of transfer matrices suitable for use with the methods described herein will become apparent to those of ordinary skill in the art, and all such variations of transfer matrices are within the scope of the present disclosure.
The algorithmic approach of equations (6) through (8) for computing equation (5) is suitable for implementation by software executing on a computer system such as the illustrative system shown in
Located within processing subsystem 330 of computer system 300 is a display interface 352, a processor 356, a peripheral interface 358, an information storage device 360, a network interface 362 and a memory 370. Bus 364 couples each of these elements to each other and transports their communications. Network interface 362 enables communications with other systems (e.g., via the Internet with a central database server housing additional modeling parameters and suitable for saving the results of the modeling). In accordance with user input received via peripheral interface 358 and program instructions from memory 370 and/or information storage device 360, processor 356 processes input from the user and applies it to the well casing string data to perform the disclosed methods and present the results to the user. Storage device 360 may be implemented using any number of known non-transitory information storage media, including but not limited to magnetic disks, solid-state storage devices and optical storage disks.
Various software modules are shown loaded into memory 370 of
If the segment associated with the current node is one that restricts fluid flow, such as a plug or a casing cross-sectional area reduction (block 406; matrix definition module 376) the transfer matrix is defined in terms of one or more forces present at the restriction (block 408; matrix definition module 376). If the segment associated with the current node is a well casing string, the transfer matrix is defined in terms of the reference node's state vector, the previous node's state vector and/or the known elements of the current node's state vector (block 410; matrix definition module 376). Once the transfer matrix is defined for node i, the cumulative transfer matrix Tacc is either initialized as T1 for i=1, or updated as the cross product of the current node's transfer matrix Ti and Tacc for i>1 (block 412; cross product module 380). The cross product is then utilized to determine the node's full state vector, which is used to determine the unknown element values of the current node's state vector (block 414; TMM module 380). In at least some illustrative embodiments, block 414 is skipped for the first node (i=1) if all the state vector values are already known (e.g., if the values are readily measurable at a starting point at the surface).
If additional well casing string segments remain (block 418; TMM module 380) node index i is incremented and the current node's state vector becomes the prior node's state vector (block 416; TMM module 380). The above-described process is then repeated for the next segment and node along the well casing string (blocks 404 through 418). If there are no additional well casing string segments (block 418; TMM module 380), the previously unknown and now calculated elements of each state vector are subsequently used as a basis for determining casing manufacturing parameters such as dimensions and composition, or as a basis for locating centralizers to position the casing within a borehole (block 420, parameter derivation module 382). The resulting manufacturing parameters or centralizer position are respectively provided as composition and/or dimensional specifications to manufacturing personnel, or as a position for a centralizer or stabilizer to well operators/component installers (block 422, presentation module 384), ending the method (block 424).
In at least some illustrative embodiments, the calculated forces at each node along a casing string are indicated on a graphical representation of the well casing string. Once the forces have been determined, the axial load (e.g., tension) present at a given node can be computed, for example by using equation (4). Casing string parameters such as segment lengths, wall thickness and material compositions may then be determined from the computed axial load. These parameters provide the required casing safety margins at a reduced cost when compared to existing methods that overestimate the required casing string parameters. In other illustrative embodiments, the positions of one or more centralizers (determined based upon the computed and displayed forces) are on the casing string's graphical representation. The calculated forces at each node are used to determine the side force to which a casing segment i (located between two nodes i and i+1) is subjected, for example by using the equation,
SideForcei=√{square root over ((Fx(i+1)−Fxi)2+(Fh(i+1)−Fhi)2)} (9).
The contact points between the casing and borehole wall can be determined from the state vector Vi and the well trajectory. The computed side force and the contact points are then used to determine a suitable centralizer(s) and the best centralizer position(s), for example, at the contact point(s) between the casing segment and the borehole wall.
Numerous other modifications, equivalents, and alternatives, will become apparent to those skilled in the art once the above disclosure is fully appreciated. For example, although the embodiments described use the TMM to determine the transfer matrices, other analytical methods are also suitable to determine transfer matrices used to determine the unknown elements of a node of interest. Further, although the examples provided are applied within the context of static modeling, static snapshots repeatedly performed over time may be combined to provide dynamic as well as near-real-time or real-time modeling of the well casing string (e.g., to predict the loading on a well casing string as a cementing plug progresses down the string). Additionally, although the disclosed embodiments describe modeling well casing strings, any of a wide variety of well components may be modeled to determine the manufacturing and/or operating parameters of said components. These well components include but are not limited to drillstrings, workstrings, production strings and coiled tubing stings. Other well components that restrict fluid flow within a running tubular string (e.g., packers) are also within the scope of the disclosure. It is intended that the following claims be interpreted to embrace all such modifications, equivalents, and alternatives where applicable.
Claims
1. A method for determining manufacturing or operating parameters for a deviated downhole well component, the method comprising:
- representing a tubular string prior to cementing in a borehole as a sequence of nodes separated by tapered segments, said nodes being numerable from i=1 to N with an initial, mechanically constrained reference node located at a casing hanger representable with i=0, a final, mechanically constrained Nth node located at a shoe representable with i=N, and each node being associated with a state vector describing a position of the corresponding node in the borehole and one or more forces present on the tubular string at said corresponding node;
- determining a sequence of transfer matrices enabling the determination of an ith node's state vector from an ith−1 node's state vector;
- defining values of an initial state vector for the reference node;
- applying the transfer matrices to obtain values for each of the state vectors;
- deriving from at least one of the state vectors a parameter value for said component, the parameter value being in a set consisting of a composition, manufacturing dimensions, and a position for a centralizer or stabilizer of the tubular string; and
- specifying a component having said parameter value, wherein said parameter value facilitates manufacture of the component or positioning of the component for cementing in the borehole.
2. The method of claim 1, wherein said specifying includes providing the composition or a dimensional specification to a manufacturer of said component.
3. The method of claim 1, wherein said specifying includes providing the position for the centralizer or stabilizer to an installer of said component.
4. The method of claim 1, wherein each state vector comprises a vertical position ui, a horizontal position vi, and an inclination angle αi, associated with node i; and further comprises a vertical force Fxi, a horizontal force Fhi and a bending moment Mi present at node i; and wherein the state vector is representable as [ui, vi, αi, Fxi, Fhi, Mi, 1]T.
5. The method of claim 4, wherein the transfer matrix for the ith node, said ith node associated with a tubular segment, is representable as: [ 1 0 - l i h ( ( l i h ) 3 6 × ( EI ) i + l i x ( EA ) i ) 0 - ( l i h ) 2 2 × ( EI ) i l i x + Δ α i - 1 0 × l i h 0 1 l i x 0 ( - ( l i x ) 3 6 × ( EI ) i + l i h ( EA ) i ) ( l i x ) 2 2 × ( EI ) i l i h - Δ α i - 1 0 × l i x 0 0 1 - ( l i h ) 2 2 × ( EI ) i ( l i x ) 2 2 × ( EI ) i ( l i x ) 2 + ( l i h ) 2 ( EI ) i Δ α i 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 l i h - l i x 1 0 0 0 0 0 0 0 1 ] i
- where, lih is a horizontal distance defined as (vi−vi−1); lix is a vertical distance defined as (ui−ui−1); Δαi−10 is a change in inclination angle at the ith−1 node defined as (αi−1−α0); Δαi0 is a change in inclination angle at the ith node defined as (αi−α0); (EI)i is a product of Young's modulus and a moment of inertia of the component at the ith node; and (EA)i is a product of Young's modulus and a cross-sectional area of the component at the ith node.
6. The method of claim 4, wherein said deriving includes deriving an axial force Fαi present at the ith node.
7. The method of claim 4, wherein the transfer matrix for the ith node, said ith node associated with a flow restriction, is representable as: [ 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0 F p x 0 0 0 0 0 1 F p h 0 0 0 0 0 0 1 ] i
- where, Fpx is a vertical force present on a plug located at the ith node; and Fph is a horizontal force present on the plug.
8. The method of claim 1, wherein the Nth node is also mechanically constrained.
9. The method of claim 1, wherein the component comprises a running string selected from the group consisting of a well casing string, a drillstring, a production string and a coil tubing.
10. The method of claim 9, wherein the running string comprises a tapered segment, a cross-section size change, a packer or a plug.
11. The method of claim 1, further comprising cementing the component having said parameter value in the borehole.
12. The method of claim 1, further comprising manufacturing the component having said parameter value or positioning the component having said parameter value in the borehole.
13. A system that determines manufacturing and operating parameters for a deviated downhole well component, the system comprising:
- a memory having deviated downhole well component modeling software; and
- one or more processors coupled to the memory, the software causing the one or more processors to: represent a tubular string prior to cementing in a borehole as a sequence of nodes separated by tapered segments, said nodes being numerable from i=1 to N with an initial, mechanically constrained reference node located at a casing hanger representable with i=0, a final, mechanically constrained Nth node located at a shoe representable with i=N, and each node being associated with a state vector describing a position of the corresponding node in the borehole and one or more forces present on the tubular string at said corresponding node; determine a sequence of transfer matrices enabling the determination of an ith node's state vector from an ith−1 node's state vector; define values of an initial state vector for the reference node; apply the transfer matrices to obtain values for each of the state vectors; derive from at least one of the state vectors a parameter value for said component, the parameter value being in a set consisting of a composition, manufacturing dimensions, and a position for a centralizer or stabilizer of the tubular string; and specify a component having said parameter value, wherein said parameter value facilitates manufacture of the component or positioning of the component for cementing in the borehole.
14. The system of claim 13, wherein the one or more processors specify said component at least in part by providing the composition or a dimensional specification to a manufacturer of said component.
15. The system of claim 14, wherein the transfer matrix for the ith node, said ith node associated with a tubular segment, is representable as: [ 1 0 - l i h ( ( l i h ) 3 6 × ( EI ) i + l i x ( EA ) i ) 0 - ( l i h ) 2 2 × ( EI ) i l i x + Δ α i - 1 0 × l i h 0 1 l i x 0 ( - ( l i x ) 3 6 × ( EI ) i + l i h ( EA ) i ) ( l i x ) 2 2 × ( EI ) i l i h - Δ α i - 1 0 × l i x 0 0 1 - ( l i h ) 2 2 × ( EI ) i ( l i x ) 2 2 × ( EI ) i ( l i x ) 2 + ( l i h ) 2 ( EI ) i Δ α i 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 l i h - l i x 1 0 0 0 0 0 0 0 1 ] i
- where, lih is a horizontal distance defined as (vi−vi−1); lix is a vertical distance defined as (ui−ui−1); Δαi−10 is a change in inclination angle at the ith−1 node defined as (αi−1−α0); Δαi0 is a change in inclination angle at the ith node defined as (αi−α0); (EI)i is a product of Young's modulus and a moment of inertia of the component at the ith node; and (EA)i is a product of Young's modulus and a cross-sectional area of the component at the ith node.
16. The system of claim 13, wherein the one or more processors specify said component at least in part by providing the position for the centralizer or stabilizer to an installer of said component.
17. The system of claim 13, wherein each state vector comprises a vertical position ui, a horizontal position vi, and an inclination angle αi, associated with node i; and further comprises a vertical force Fxi, a horizontal force Fhi and a bending moment Mi present at node i; and wherein the state vector is representable as [ui, vi, αi, Fxi, Fhi, Mi, 1]T.
18. The system of claim 17, wherein the one or more processors derive said parameter value at least in part by deriving an axial force Fαi present at the ith node.
19. The system of claim 17, wherein the transfer matrix for the ith node, said ith node associated with a flow restriction, is representable as: [ 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0 F p x 0 0 0 0 0 1 F p h 0 0 0 0 0 0 1 ] i
- where, Fpx is a vertical force present on a plug located at the ith node; and Fph is a horizontal force present on the plug.
20. The system of claim 13, wherein the Nth node is also mechanically constrain ed.
21. The system of claim 13, wherein the component comprises a running string selected from the group consisting of a well casing string, a drillstring, a production string and a coil tubing.
22. The system of claim 21, wherein the running string comprises a tapered segment, a cross-section size change, a packer or a plug.
23. The system of claim 13, further comprising cementing the component having said parameter value in the borehole.
24. The system of claim 13, further comprising manufacturing the component having said parameter value or positioning the component having said parameter value in the borehole.
6785641 | August 31, 2004 | Huang |
6816082 | November 9, 2004 | Laborde |
20070024464 | February 1, 2007 | Lemenager et al. |
20100032165 | February 11, 2010 | Bailey et al. |
20110253447 | October 20, 2011 | Gleitman |
20110320047 | December 29, 2011 | Stone et al. |
20120203515 | August 9, 2012 | Pita |
2014/204521 | December 2014 | WO |
- PCT International Search Report and Written Opinion, dated Apr. 29, 2014, Appl No. PCT/US2014/011998, “Methods and Systems for Determining Manufacturing and Operating Parameters for a Deviated Downhole Well Component,” Filed Jan. 17, 2014, 10 pgs.
- Canadian Application Serial No. 2,913,203; Canadian Examination Report; Apr. 26, 2018, 6 pages.
- “CA Office Action”, Dec. 22, 2016, Appl No. 2,913,203, “Methods and Systems for Determining Manufacturing and Operating Parameters for a Deviated Downhole Well Component,” Filed Jan. 17, 2014.
- “International Preliminary Report on Patentability”, dated Dec. 30, 2015, Appl No. PCT/US2014/011998, “Methods and Systems for Determining Manufacturing and Operating Parameters for a Deviated Downhole Well Component,” Filed Jan. 17, 2014, 6.
- Fazaelizadeh, “Real Time Torque and Drag Analysis during Directional Drilling”, Fazaelizadeh, Mohammad, “Real Time Torque and Drag Analysis during Directional Drilling,” A Thesis Submitted to Faculty of Graduate Studies Partial Fulfilment of the Requirements for the Degree of Doctor of Philosophy, Mar. 2013, Calgary, Alberta, 200.
- Perneder, “A Three-Dimensional Mathematical Model of Directional Drilling”, Perneder, Luc, “A Three-Mathematical Model of Directional Drilling,” a Dissertation Submitted to the Faculty of the Graduate School of the University of Minnesota, Jan. 2013, 143.
- Tveitdal, “Torque & drag analyses of North Sea Wells using new 3D model.”, Tveitdal, Terje, “Torque & drag analyses of North Sea Wells using new 3D model.,” University of Stavanger Faculty of Science and Technology Master's Thesis, Jun. 15, 2011, 75.
Type: Grant
Filed: Jan 17, 2014
Date of Patent: Aug 20, 2019
Patent Publication Number: 20160108704
Assignee: LANDMARK GRAPHICS CORPORATION (Houston, TX)
Inventors: Robello Samuel (Cypress, TX), Zhenying Wang (Missouri City, TX)
Primary Examiner: Manuel L Barbee
Application Number: 14/895,165
International Classification: E21B 41/00 (20060101); E21B 7/04 (20060101); G06F 17/16 (20060101); E21B 17/00 (20060101);