CN106529092A - Finite element calculating method of shear deformation force of well casing - Google Patents
Finite element calculating method of shear deformation force of well casing Download PDFInfo
- Publication number
- CN106529092A CN106529092A CN201611154853.6A CN201611154853A CN106529092A CN 106529092 A CN106529092 A CN 106529092A CN 201611154853 A CN201611154853 A CN 201611154853A CN 106529092 A CN106529092 A CN 106529092A
- Authority
- CN
- China
- Prior art keywords
- finite element
- shearing
- well casing
- sleeve pipe
- casing
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
Abstract
The invention provides a finite element calculating method of shear deformation force of a well casing. According to the method, the calculation of the shear deformation force of the well casing caused by diastrophism of stratums is simulated based on finite element software. The method comprises 11 steps. The finite element calculating method of the shear deformation force of the well casing has the advantages that through the finite element calculating method of the shear deformation force of the well casing, the morphology of a pipe body of the casing when the casing is subject to the effect of stratum diastrophism force under a well can be accurately simulated, and the size of the biggest shear and diastrophism force the pipe body can bear can also be accurately simulated; through the finite element calculating, the cost of conducting shear tests on casing pipes of well drilling and well cementation can be sharply reduced. The result of the shear force calculated through the finite element calculating method is compared with an object test, and the error is within 5%. The calculated result obtained through the finite element calculating method has an important guiding significance in the design of pipe strings of the casing pipe.
Description
Technical field
The present invention relates to drilling well and oil recovery technique field, are related to a kind of detrusion power FEM calculation side of well casing
Method.
Background technology
During drilling well extraction, stratum after sand production in oil, can be caused loose, form cavity.Cause the loss of propping of subterranean formations
And sliding, horizontal shear action will be produced to well drilling and cementing sleeve pipe, it is easy to cause the failure by shear of sleeve pipe even completely wrong
It is disconnected, drilling well can be made to be forced to stop extraction, cause huge economic loss.
Chinese patent literature number discloses a kind of well drilling and cementing reservoir hunch pilot system for ZL201610333516.7, can
Shearing force is applied with the well cementing casing pipe to different tube diameters, for simulating the shear action of the relative displacement to well cementing casing pipe of rock stratum.
But the document can only obtain the detrusion power of sleeve pipe by the method tested, test expends substantial amounts of physical resources and financial resources, and
The maximum shear load of testing equipment has certain restriction, and the test specification to being carried out also has certain limitation.
The content of the invention
It is an object of the invention to provide a kind of detrusion power finite element method of well casing, can be with quantitative analysis
The maximum shear stress that the lateral displacement change of sleeve pipe and sleeve pipe can bear, can be that prevention casing failure and bending deformation provide reason
By foundation, size of the sleeve pipe in the failure by shear power of down-hole is determined, so as to preferably select suitable sleeve pipe for oil field.For drilling well
Engineering provides technical support.
For achieving the above object, the detrusion power finite element method of the well casing that the present invention is provided includes following
Step：
1) the three-dimensional thin cylinder geometrical model of well casing is set up by finite element software mapping, outside the thin cylinder
Footpath is the outer diameter D of well casing, and wall thickness t of the wall thickness for sleeve pipe, thin cylinder length L are 1000mm；
2) the body overall length symmetric position by finite element software in above-mentioned thin cylinder geometrical model sets up two shearings
The three-dimensional thin cylinder geometrical model of ring, the internal diameter dj of shearing ring are identical with outer tube diameter D, and outer diameter D j of shearing ring is more than body
Outer diameter D, length s of shearing ring is 20～60mm, two shearing rings apart from a be 127mm；
3) described sleeve pipe is carbon steel or low alloy steel material, is covered in finite element software according to the featured configuration of material
The elastic modelling quantity of pipe, stress-strain diagram, Poisson's ratio, the performance of density of material；
4) it is rigid body to arrange two shearing rings, i.e., any deformation occurs in no instance；
5) in finite element software, before calculating, the loading procedure for arranging is needed to be single or multiple steps i.e. analysis step,
The detrusion of well casing is set as quasistatic, the form for arranging analysis step is analyzed for Dynamic Implicit, the time that calculates is 1 second,
Iteration duration 0.01 second；
6) geometrical model of the sleeve pipe three-dimensional thin cylinder to being set up carries out stress and strain model, and the mode of stress and strain model is first
Seed point is set, after carry out the automatic grid division of finite element software, 3 layers of grid are at least divided on wall thickness direction, grid
Form is hexahedral element；
7) radial displacement for constraining the sleeve pipe both ends of the surface is 0；
8) it is surface-to-surface contact to arrange the contact form between the outer surface and the inner surface of two shearing rings of sleeve pipe；
9) to an applying radial displacement in two shearing rings, another barycenter applies position in opposite direction equal in magnitude
Move, shift value is at least 4mm；
10) realize that the computing of the pre-treatment step completed by step 1～9 is solved by the solver of finite element software, meter
The result of calculation will be stored in post processing file；
11) value of the shearing force suffered by body is extracted in the result of calculation obtained from step 10, and is made with shearing
The shear force value curve of the body that the radial displacement of ring increases, is analyzed to the curve, and the value of shearing force is with shearing ring
Displacement increase and speedup slows down, begin with the mutation uprushed with bust afterwards, the maximum before mutation be described sleeve pipe
The shearing force maximum that can be born.
The effect of the present invention is, by the detrusion power finite element method of well casing, accurately can to simulate
Go out the pattern of sleeve pipe body when down-hole is acted on by stratum changing of the relative positions power, and the maximum shear changing of the relative positions power that body can bear
Size.Can be greatly lowered by FEM calculation carries out the cost of well drilling and cementing reservoir hunch test.What the present invention was calculated
Shearing force result is contrasted with actual loading test, and within 5%, the result of calculation of the present invention is designed with error to sleeve pipe pipe string
Important directive significance.
Description of the drawings
Fig. 1 is the three-dimensional finite element model schematic diagram of the detrusion power finite element method of the sleeve pipe of the present invention；
Fig. 2 is the shearing ring schematic diagram of the present invention
Fig. 3 is the sleeve pipe stress and strain model schematic diagram of the present invention；
Fig. 4 is the shearing force born by the sleeve pipe that the detrusion power finite element method of the sleeve pipe of the present invention is extracted
With the curve map of shearing ring alternate displacement.
Fig. 5-1,5-2 are that the result of finite element of the detrusion power finite element method of the sleeve pipe of the present invention is illustrated
Figure.
In figure：
A：Body；B：Shearing ring；C：By the body deformed after shearing force；D：Outer tube diameter；t：Body wall thickness；L：Body
Length；s：The length degree of shearing ring；a：The distance between two shearing rings；Dj：Shearing ring external diameter；dj：Shearing ring internal diameter
Specific embodiment
The detrusion power finite element method of the well casing of the present invention is illustrated with reference to accompanying drawing.
The detrusion power finite element method of the well casing that the present invention is provided is comprised the following steps：
1) the three-dimensional thin cylinder geometrical model of well casing is set up by finite element software mapping, outside the thin cylinder
Footpath is the outer diameter D of well casing, and wall thickness is the wall thickness t of sleeve pipe, and thin cylinder length L is 1000mm, as shown in the A of Fig. 1.
2) the body overall length symmetric position by finite element software in above-mentioned thin cylinder geometrical model sets up two shearings
The three-dimensional thin cylinder geometrical model of ring, the internal diameter dj of shearing ring are identical with outer tube diameter D, and outer diameter D j of shearing ring is more than body
Outer diameter D, length s of shearing ring is 20～60mm, two shearing rings apart from a be 127mm, as shown in the B and Fig. 2 of Fig. 1.
3) described sleeve pipe is carbon steel or low alloy steel material, is covered in finite element software according to the featured configuration of material
The elastic modelling quantity of pipe, stress-strain diagram, Poisson's ratio, the performance of density of material.
4) it is rigid body to arrange two shearing rings, i.e., any deformation occurs in no instance.
5) in finite element software, before calculating, the loading procedure for arranging is needed to be single or multiple steps i.e. analysis step,
The detrusion of well casing is set as quasistatic, the form for arranging analysis step is analyzed for Dynamic Implicit, the time that calculates is 1 second,
Iteration duration 0.01 second.
6) geometrical model of the sleeve pipe three-dimensional thin cylinder to being set up carries out stress and strain model, and the mode of stress and strain model is first
Seed point is set, after carry out the automatic grid division of finite element software, 3 layers of grid are at least divided on wall thickness direction, grid
Form is hexahedral element, as shown in Figure 3.
7) radial displacement for constraining the sleeve pipe both ends of the surface is 0.
8) it is surface-to-surface contact to arrange the contact form between the outer surface and the inner surface of two shearing rings of sleeve pipe；
9) to an applying radial displacement in two shearing rings, another barycenter applies position in opposite direction equal in magnitude
Move, shift value is at least 4mm.
10) this step is solved to the FEM model that 1～9 step is set up.When solving, equation is spatially adopted
After Finite Element Method is discrete, become ODE：
F=M (u)+C (u)+K (u)
The equation, the displacement of any instant, speed, acceleration are solved by Dynamic Implicit analysis NewMark methods used
It is all interrelated, realize solving using iteration and solution simultaneous equations.The result of calculating will be stored in post processing file.
11) from step 10) value of shearing force suffered by body is extracted in the result of calculation that obtains, and make with shearing
The shear force value curve of the body that the radial displacement of ring increases, is analyzed to the curve, and the value of shearing force is with shearing ring
Displacement increase and speedup slows down, begin with the mutation uprushed with bust afterwards, the maximum before mutation be described sleeve pipe
The shearing force maximum that can be born.As shown in figure 4, the shearing force that the maximum before mutation is described sleeve pipe can bear is maximum
Value.Fig. 5-1,5-2 are that body deforms schematic diagram caused by shearing force.
Two concrete calculated examples are given below.
Embodiment 1
Casing gauge especially footpath D=139.7mm, wall thickness t=12.7mm, grade of steel Q125, its elastic modulus E=206GPa, pool
Pine is than μ=0.3, yield strength σS=862～1034MPa, tensile strength sigmab>931MPa.Carried out according to above-mentioned Finite Element
Modeling, tube length 1000mm, two shearing ring width s=40mm, at a distance of a=127mm.Fig. 4 is to be carried by finite element software
The shearing force for taking and the graph of a relation of the distance of the movement of shearing ring, can obtain maximum 1625KN of the shearing force before mutation,
So the value is the maximum shear stress that the sleeve pipe can bear.
Embodiment 2
Casing gauge especially footpath D=146.7mm, wall thickness t=15.9mm, grade of steel TP140V, its elastic modulus E=206GPa,
Poisson's ratio μ=0.3, yield strength σS=985～1080MPa, tensile strength sigmab>1034MPa.According to above-mentioned Finite Element
It is modeled, tube length 1000mm, two shearing ring width s=40mm, at a distance of a=127mm.Can be with by finite element software
Maximum 2405KN of the shearing force before mutation is obtained, so the value is the maximum shear stress that the sleeve pipe can bear.
Embodiment 3
Casing gauge especially footpath D=177.8mm, wall thickness t=10.36mm, grade of steel P110, its elastic modulus E=206GPa, pool
Pine is than μ=0.3, yield strength σS=758～965MPa, tensile strength sigmab>862MPa.Carried out according to above-mentioned Finite Element
Modeling, tube length 1000mm, two shearing ring width s=50mm, at a distance of a=127mm.Can be obtained by finite element software
Maximum 1454KN of the shearing force before mutation, so the value is the maximum shear stress that the sleeve pipe can bear.
Claims (1)
1. the detrusion power finite element method of a kind of well casing, the method are simulated due to ground based on finite element software
The changing of the relative positions of layer causes well casing that the calculating of detrusion occurs, and it is characterized in that：The method is comprised the following steps：
1) the three-dimensional thin cylinder geometrical model of well casing is set up by finite element software mapping, the external diameter of the thin cylinder is
The outer diameter D of well casing, wall thickness t of the wall thickness for sleeve pipe, thin cylinder length L are 1000mm；
2) the body overall length symmetric position by finite element software in above-mentioned thin cylinder geometrical model sets up two shearing rings
Three-dimensional thin cylinder geometrical model, the internal diameter dj of shearing ring are identical with outer tube diameter D, and outer diameter D j of shearing ring is more than outer tube diameter
D, length s of shearing ring is 20～60mm, two shearing rings apart from a be 127mm；
3) described sleeve pipe is carbon steel or low alloy steel material, according to the featured configuration described sleeve pipe of material in finite element software
Elastic modelling quantity, stress-strain diagram, Poisson's ratio, the performance of density of material；
4) it is rigid body to arrange two shearing rings, i.e., any deformation occurs in no instance；
5) in finite element software, before calculating, need the loading procedure for arranging to be single or multiple steps i.e. analysis step, set
The detrusion of well casing is quasistatic, and the form for arranging analysis step is analyzed for Dynamic Implicit, and the time that calculates is 1 second, iteration
Duration 0.01 second；
6) geometrical model of the sleeve pipe three-dimensional thin cylinder to being set up carries out stress and strain model, and the mode of stress and strain model is first to arrange
Seed point, after carry out the automatic grid division of finite element software, 3 layers of grid, the form of grid are at least divided on wall thickness direction
For hexahedral element；
7) radial displacement for constraining the well casing both ends of the surface is 0；
8) it is surface-to-surface contact to arrange the contact form between the outer surface of well casing and the inner surface of two shearing rings；
9) to an applying radial displacement in described two shearing rings, another barycenter applies position in opposite direction equal in magnitude
Move, displacement is at least 4mm；
10) realize that the computing of the pre-treatment step completed by step 1～9 is solved by the solver of finite element software, calculating
As a result will be stored in post processing file；
11) from step 10) value of shearing force suffered by body is extracted in the result of calculation that obtains, and make with shearing ring
The shear force value curve of the body that radial displacement increases, is analyzed to the curve, and the value of shearing force is with the position of shearing ring
Moving increases and speedup slows down, and begins with the mutation uprushed with bust afterwards, and the maximum before mutation can be held for described sleeve pipe
The shearing force maximum received.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611154853.6A CN106529092B (en) | 2016-12-14 | 2016-12-14 | Finite element calculation method for shearing deformation force of well drilling casing |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611154853.6A CN106529092B (en) | 2016-12-14 | 2016-12-14 | Finite element calculation method for shearing deformation force of well drilling casing |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106529092A true CN106529092A (en) | 2017-03-22 |
| CN106529092B CN106529092B (en) | 2020-01-21 |
Family
ID=58339550
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611154853.6A Active CN106529092B (en) | 2016-12-14 | 2016-12-14 | Finite element calculation method for shearing deformation force of well drilling casing |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106529092B (en) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106991235A (en) * | 2017-04-01 | 2017-07-28 | 中国石油天然气集团公司 | Cementing concrete ring integrity analysis Forecasting Methodology and device |
| CN107169240A (en) * | 2017-06-22 | 2017-09-15 | 贵州财经大学 | Returning based on shoulder hole drags computational methods and device |
| CN108279173A (en) * | 2018-01-02 | 2018-07-13 | 中国石油天然气集团公司 | A kind of casing anti-shear performance evaluation method |
| CN110717284A (en) * | 2019-08-16 | 2020-01-21 | 中国石油天然气集团有限公司 | Analysis and test evaluation method for shear-resistant bearing capacity of casing |
| CN111625979A (en) * | 2020-05-28 | 2020-09-04 | 中国船舶工业集团公司第七0八研究所 | Strength checking and loading method for pile leg bolt hole |
| CN113673126A (en) * | 2021-07-28 | 2021-11-19 | 中国石油大学(北京) | Method and device for calculating finite element of annular space with pressure of multilayer casing for well drilling |
| CN113824070A (en) * | 2021-09-22 | 2021-12-21 | 深圳市骏鼎达新材料股份有限公司 | Protective sleeve with self-positioning function |
| CN113836656A (en) * | 2021-09-14 | 2021-12-24 | 临海伟星新型建材有限公司 | Calculation method for reducing amount of fluorine-silicon modified PERT barrier liner pipe by adopting finite element algorithm |
| WO2022227488A1 (en) * | 2021-04-28 | 2022-11-03 | 天津钢管制造有限公司 | Method for obtaining external pressure and collapse resistance ability of double-layer casing |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105760564A (en) * | 2014-12-19 | 2016-07-13 | 中国石油天然气股份有限公司 | Method and device for analyzing oil-string casing failure |
| CN106018125A (en) * | 2016-05-19 | 2016-10-12 | 中国矿业大学 | Shear experiment system for drilling and cementing casings |
-
2016
- 2016-12-14 CN CN201611154853.6A patent/CN106529092B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105760564A (en) * | 2014-12-19 | 2016-07-13 | 中国石油天然气股份有限公司 | Method and device for analyzing oil-string casing failure |
| CN106018125A (en) * | 2016-05-19 | 2016-10-12 | 中国矿业大学 | Shear experiment system for drilling and cementing casings |
Non-Patent Citations (2)
| Title |
|---|
| 于浩 等: "页岩气体积压裂过程中套管失效机理研究", 《中国安全生产科学技术》 * |
| 蒋可 等: "页岩气水平井固井质量对套管损坏的影响", 《天然气工业》 * |
Cited By (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106991235B (en) * | 2017-04-01 | 2020-03-10 | 中国石油天然气集团公司 | Method and device for analyzing and predicting integrity of well cementation cement sheath |
| CN106991235A (en) * | 2017-04-01 | 2017-07-28 | 中国石油天然气集团公司 | Cementing concrete ring integrity analysis Forecasting Methodology and device |
| CN107169240A (en) * | 2017-06-22 | 2017-09-15 | 贵州财经大学 | Returning based on shoulder hole drags computational methods and device |
| CN107169240B (en) * | 2017-06-22 | 2020-09-11 | 贵州财经大学 | Back-dragging calculation method and device based on stepped hole |
| CN108279173A (en) * | 2018-01-02 | 2018-07-13 | 中国石油天然气集团公司 | A kind of casing anti-shear performance evaluation method |
| CN108279173B (en) * | 2018-01-02 | 2020-08-07 | 中国石油天然气集团公司 | Method for evaluating anti-shearing performance of casing |
| CN110717284A (en) * | 2019-08-16 | 2020-01-21 | 中国石油天然气集团有限公司 | Analysis and test evaluation method for shear-resistant bearing capacity of casing |
| CN111625979A (en) * | 2020-05-28 | 2020-09-04 | 中国船舶工业集团公司第七0八研究所 | Strength checking and loading method for pile leg bolt hole |
| CN111625979B (en) * | 2020-05-28 | 2023-08-15 | 中国船舶工业集团公司第七0八研究所 | Strength checking and loading method for pile leg bolt holes |
| WO2022227488A1 (en) * | 2021-04-28 | 2022-11-03 | 天津钢管制造有限公司 | Method for obtaining external pressure and collapse resistance ability of double-layer casing |
| CN113673126A (en) * | 2021-07-28 | 2021-11-19 | 中国石油大学(北京) | Method and device for calculating finite element of annular space with pressure of multilayer casing for well drilling |
| CN113673126B (en) * | 2021-07-28 | 2024-02-13 | 中国石油大学(北京) | Method and device for calculating annular space pressure finite element of multilayer casing for well drilling |
| CN113836656A (en) * | 2021-09-14 | 2021-12-24 | 临海伟星新型建材有限公司 | Calculation method for reducing amount of fluorine-silicon modified PERT barrier liner pipe by adopting finite element algorithm |
| WO2023039932A1 (en) * | 2021-09-14 | 2023-03-23 | 临海伟星新型建材有限公司 | Method for calculating diameter reduction amount of fluorosilicone-modified pert barrier liner pipe by using finite element algorithm |
| CN113824070B (en) * | 2021-09-22 | 2022-11-18 | 深圳市骏鼎达新材料股份有限公司 | Protective sleeve with self-positioning function |
| CN113824070A (en) * | 2021-09-22 | 2021-12-21 | 深圳市骏鼎达新材料股份有限公司 | Protective sleeve with self-positioning function |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106529092B (en) | 2020-01-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106529092A (en) | Finite element calculating method of shear deformation force of well casing | |
| CN105319603B (en) | The Forecasting Methodology of tight sandstone reservoir complexity chicken-wire cracking | |
| CN109271738B (en) | Numerical inversion method for acquiring Weibull distribution parameters of roadway surrounding rock | |
| CN104361211B (en) | The structure and methods for using them of Statistical Damage Constitutive Model For Rock | |
| CN112036098A (en) | Method for simulating hydraulic fracture propagation numerical value of deep oil and gas reservoir | |
| US9836561B2 (en) | Optimizing multistage hydraulic fracturing design based on three-dimensional (3D) continuum damage mechanics | |
| CN113820750B (en) | Method for quantitatively predicting mudstone structural cracks based on elastoplastic mechanics | |
| CN110005407B (en) | Salt cavern gas storage stability evaluation method | |
| CN104963657A (en) | Oil reservoir numerical simulation method and apparatus based on high water-cut period monitoring data constraints | |
| CN109033504B (en) | Oil-water well casing damage prediction method | |
| CN105089595A (en) | Oil reservoir numerical simulation method and device under the action of horizontal fracturing fracture diversion | |
| CN106372297A (en) | Method for determining safe vertical distance between shield and karst cave in sand karst stratum | |
| CN106555575A (en) | A kind of computational methods of the stratum initial cracking pressure of deep layer shale | |
| CN115267905B (en) | Method for predicting crack lost circulation in drilling engineering of complex structural area | |
| CN108952657B (en) | Horizontal well platform fracturing fracture length determination method and device | |
| CN107991188B (en) | Method for predicting hydraulic fracture complexity based on core residual stress level | |
| Yu et al. | On how asymmetric stimulated rock volume in shales may impact casing integrity | |
| Zhang et al. | Strata movement and stress evolution when mining two overlapping panels affected by hard stratum | |
| CN106372295A (en) | Method for determining safe horizontal distance between shield and karst cave in sand karst stratum | |
| CN115508201A (en) | Fracturing string safety checking method | |
| WO2017030838A1 (en) | Mechanisms-based fracture model for geomaterials | |
| Yan et al. | A multiscale analysis of adjacent fault dislocation mechanism induced by tunnel excavation based on continuous-discrete coupling method | |
| CN106844858B (en) | Stratum fracture development zone prediction method and device | |
| CN114458274B (en) | Rock capacity expansion method | |
| Chen et al. | Case study: Treating parameters optimization in casing deformation wells of shale gas in Sichuan Basin |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20190628 Address after: No. 396 Jintang Highway, Dongli District, Tianjin 300301 Applicant after: Tianjin Steel Tube Manufacturing Co., Ltd. Address before: No. 396 Jintang Highway, Dongli District, Tianjin 300301 Applicant before: Tianjin Steel Pipe Group Co., Ltd. |
|
| TA01 | Transfer of patent application right | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |