CN103615233A - Method and equipment for measuring axial force and lateral force of casing in three-dimensional borehole - Google Patents

Method and equipment for measuring axial force and lateral force of casing in three-dimensional borehole Download PDF

Info

Publication number
CN103615233A
CN103615233A CN201310629634.9A CN201310629634A CN103615233A CN 103615233 A CN103615233 A CN 103615233A CN 201310629634 A CN201310629634 A CN 201310629634A CN 103615233 A CN103615233 A CN 103615233A
Authority
CN
China
Prior art keywords
unit
described micro
micro unit
axial force
relational expression
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.)
Pending
Application number
CN201310629634.9A
Other languages
Chinese (zh)
Inventor
马振
齐海鹰
曹传文
彭松良
曲绍刚
张成博
刘洪芹
宋阳
钟满发
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Petrochina Co Ltd
Original Assignee
Petrochina Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN201310629634.9A priority Critical patent/CN103615233A/en
Publication of CN103615233A publication Critical patent/CN103615233A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Force Measurement Appropriate To Specific Purposes (AREA)

Abstract

The invention provides a method and equipment for measuring axial force and lateral force of a casing in a three-dimensional borehole, wherein the method comprises the following steps: acquiring track measuring point data of a three-dimensional borehole; selecting a casing string between any two data points from the track measuring point data as a microcell; the method comprises the following steps of collecting curvature, length, effective gravity, moment of inertia of a cross section, elastic modulus, first inclination angle, second inclination angle, first azimuth angle, second azimuth angle and friction coefficient of a well hole of a micro unit; determining the axial force of the second end, the axial force of the first end and the lateral force of the micro unit per unit length; and determining the axial force and the lateral force of the casing in the three-dimensional borehole according to the axial force of the second end of the microcell, the axial force of the first end and the lateral force of the microcell per unit length. The method can be used for casing strength design, has very important reference value for scheme design and adjustment of drilling, well completion, oil production and well repair operation, and is also a basis for evaluating the collapse resistance of the casing.

Description

A kind of method and apparatus of measuring three-dimensional well middle sleeve axial force and lateral force
Technical field
The present invention is about oil exploitation and bore completion technique field, particularly about the designing technique of three-dimensional well middle sleeve intensity, is a kind of method and apparatus of measuring three-dimensional well middle sleeve axial force and lateral force concretely.
Background technology
In the oil production rate of China, the ratio that viscous crude occupies increases year by year.Thickened oil recovery mainly adopts thermodynamic oil extraction process, comprises that pit shaft heating, hot water flooding, hot vapour drive, steam soak and combustion in situ etc., and what wherein really realize large-scale application is steam injection.Yet steam injection, when improving exploitation economic benefit, has also injured sleeve pipe strata structure around, thereby has even caused oil well to be scrapped the application life of having reduced oil well.Sleeve pipe is as the passage on producing zone and ground, and state quality is being related to the life-span of safety in production and Oil/gas Well, also affects drilling cost simultaneously.
At present, in prior art, generally believe that thermal production well sleeve damages too high the causing of axial heat expansion stress, therefore the axial heat expansion stress of usining must not surpass tubing yield limit as design criterion, thereby the axial force of wear sleeve and the prediction of lateral force or assessment, the schematic design of drilling well, completion, oil recovery and workover treatment and adjustment are all had to very important reference value.And the axial force and the lateral force that accurately calculate wear sleeve are the anti-bases of squeezing the ability of ruining of assessment sleeve pipe.
The distortion of directional well, horizontal well middle sleeve, to be subject to force and motion be the complicated mechanics difficult problem of sleeve pipe in long and narrow crooked hole.In prior art, the general just crooked bending with combined axial and lateral load beam theory that adopts of Exact Solution of this difficult problem, according to three moments euqation, write out multivariate linear equations, then with computer simultaneous solution, obtain, but formula, data that this method relates to are a lot of, and formula and method are all very complicated, computational process is loaded down with trivial details tediously long, greatly reduce the determination efficiency of axial force and the lateral force of wear sleeve, and then affected the efficiency to the schematic design of drilling well, completion, oil recovery and workover treatment and adjustment according to the axial force of wear sleeve and lateral force.
Summary of the invention
The problems referred to above that exist in order to overcome prior art, the invention provides a kind of method and apparatus of measuring three-dimensional well middle sleeve axial force and lateral force, utilize finite difference method to propose a kind of approximate, easy scheme, based on hole trajectory data, sleeve pipe is divided into the micro unit section of unit length, by sleeve pipe lowermost end, started to calculate piecemeal from lower to upper, meet petroleum production engineering requirement, not only can be used for casing strength design, to drilling well, completion, the schematic design of oil recovery and workover treatment and adjustment all have very important reference value, still assess the anti-basis of squeezing the ability of ruining of sleeve pipe.
One of object of the present invention is, a kind of method of measuring three-dimensional well middle sleeve axial force and lateral force is provided, and comprising: the track measuring point data that obtains three-dimensional well; From described track measuring point data, choose casing string between any two data points as a micro unit; Gather curvature, length, the effective gravity of described micro unit, the moment of inertia of cross section of described micro unit, the modulus of elasticity of described micro unit; Gather the first hole angle, the second hole angle, first party parallactic angle and the second party parallactic angle of the track measuring point that described micro unit is corresponding, the coefficient of friction resistance of well; The axial force of the second end of determining described micro unit according to the coefficient of friction resistance of the moment of inertia of the curvature of described micro unit, length, effective gravity, cross section, modulus of elasticity and the first hole angle, the second hole angle, first party parallactic angle and second party parallactic angle, well is, the lateral force of the axial force of first end, described micro unit unit length; According to the axial force of the second end of described micro unit, the lateral force of the axial force of first end, described micro unit unit length is determined axial force and the lateral force of three-dimensional well middle sleeve.
One of object of the present invention is, a kind of equipment of measuring three-dimensional well middle sleeve axial force and lateral force is provided, and comprising: track measuring point data acquisition device, for obtaining the track measuring point data of three-dimensional well; Micro unit selecting device, for choosing casing string between any two data points as a micro unit from described track measuring point data; The first harvester, for gathering curvature, length, the effective gravity of described micro unit, the moment of inertia of cross section of described micro unit, the modulus of elasticity of described micro unit; The second harvester, for gathering the first hole angle, the second hole angle, first party parallactic angle and the second party parallactic angle of the track measuring point that described micro unit is corresponding, the coefficient of friction resistance of well; Micro unit is to power determinator, for the axial force of the second end of determining described micro unit according to the coefficient of friction resistance of the moment of inertia of the curvature of described micro unit, length, effective gravity, cross section, modulus of elasticity and the first hole angle, the second hole angle, first party parallactic angle and second party parallactic angle, well, the lateral force of the axial force of first end, described micro unit unit length; Sleeve pipe is to power determinator, for according to the axial force of the second end of described micro unit, the lateral force of the axial force of first end, described micro unit unit length is determined axial force and the lateral force of three-dimensional well middle sleeve.
Beneficial effect of the present invention is, a kind of method and apparatus of measuring three-dimensional well middle sleeve axial force and lateral force is provided, utilize the basic assumption of finite difference method, utilize the method for space coordinate transformation and mathematical analysis, the casing string of take between any two data points is a unit, derived casing string lower unit axial force, unit axial force on tubing string, the axial force that the factors such as tubing string accelerated motion and liquid resistance produce on tubular element, the lateral force that tubing string distortion causes, total lateral force in full-shape plane, total lateral force in binormal direction, the series of parameters equation of the total lateral force of tubular element, utilize finite difference method to propose a kind of approximate, easy scheme, based on hole trajectory data, sleeve pipe is divided into the micro unit section of unit length, by sleeve pipe lowermost end, started to calculate piecemeal from lower to upper, meet petroleum production engineering requirement, not only can be used for casing strength design, to drilling well, completion, the schematic design of oil recovery and workover treatment and adjustment all have very important reference value, still assess the anti-basis of squeezing the ability of ruining of sleeve pipe.
For above and other object of the present invention, feature and advantage can be become apparent, preferred embodiment cited below particularly, and coordinate appended graphicly, be described in detail below.
Accompanying drawing explanation
In order to be illustrated more clearly in the embodiment of the present invention or technical scheme of the prior art, to the accompanying drawing of required use in embodiment or description of the Prior Art be briefly described below, apparently, accompanying drawing in the following describes is only some embodiments of the present invention, for those of ordinary skills, do not paying under the prerequisite of creative work, can also obtain according to these accompanying drawings other accompanying drawing.
A kind of flow chart of measuring the method for three-dimensional well middle sleeve axial force and lateral force that Fig. 1 provides for the embodiment of the present invention;
Fig. 2 is the particular flow sheet of the step S105 in Fig. 1;
Fig. 3 is the particular flow sheet of the step S201 in Fig. 2;
Fig. 4 is the particular flow sheet of the step S202 in Fig. 2;
Fig. 5 is the particular flow sheet of the step S203 in Fig. 2;
A kind of structured flowchart of measuring the equipment of three-dimensional well middle sleeve axial force and lateral force that Fig. 6 provides for the embodiment of the present invention;
Fig. 7 be in a kind of equipment of measuring three-dimensional well middle sleeve axial force and lateral force provided by the invention micro unit to the structured flowchart of power determinator 500;
Fig. 8 is the structured flowchart of the first relational expression determination module 501 in a kind of equipment of measuring three-dimensional well middle sleeve axial force and lateral force provided by the invention;
Fig. 9 is the structured flowchart of the second relational expression determination module 502 in a kind of equipment of measuring three-dimensional well middle sleeve axial force and lateral force provided by the invention;
Figure 10 is the structured flowchart of the 3rd relational expression determination module 503 in a kind of equipment of measuring three-dimensional well middle sleeve axial force and lateral force provided by the invention;
Figure 11 is the schematic diagram of setting up the relational expression of axial load and other factors;
Figure 12 is is that the sleeve pipe that node, sleeve pipe are divided between unit, any two data points is a cell schematics by hole trajectory data point.
The specific embodiment
Below in conjunction with the accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is clearly and completely described, obviously, described embodiment is only the present invention's part embodiment, rather than whole embodiment.Embodiment based in the present invention, those of ordinary skills, not making the every other embodiment obtaining under creative work prerequisite, belong to the scope of protection of the invention.
The distortion of directional well, horizontal well middle sleeve, to be subject to force and motion be the complicated mechanics difficult problem of sleeve pipe in long and narrow crooked hole.The just crooked bending with combined axial and lateral load beam theory of general employing that solves of the prior art, according to three moments euqation, write out multivariate linear equations, then with computer simultaneous solution, obtain, but formula, data that this method relates to are a lot of, and formula and method are all very complicated, computational process is loaded down with trivial details tediously long, greatly reduce the determination efficiency of axial force and the lateral force of wear sleeve, and then affected the efficiency to the schematic design of drilling well, completion, oil recovery and workover treatment and adjustment according to the axial force of wear sleeve and lateral force.
Based on this, a kind of method of measuring three-dimensional well middle sleeve axial force and lateral force that the present invention proposes, the particular flow sheet that Fig. 1 is the method, as shown in Figure 1, described method comprises:
S101: the track measuring point data that obtains three-dimensional well.
S102: choose casing string between any two data points as a micro unit from described track measuring point data.The present invention, in order to set up the universal model that calculates three-dimensional well middle sleeve axial load, first considers a micro unit between two well track measuring points.The hole trajectory data point of take is node, and sleeve pipe is divided into unit, and the sleeve pipe between any two data points is a micro unit, as shown in figure 12.
S103: gather curvature, length, the effective gravity of described micro unit, the moment of inertia of cross section of described micro unit, the modulus of elasticity of described micro unit.In concrete embodiment, can establish L sfor the length of described micro unit, the curvature that K is micro unit, the effective gravity that q is described micro unit, the moment of inertia of the cross section that I is described micro unit, the modulus of elasticity of the steel that E is described micro unit.
S104: gather the first hole angle, the second hole angle, first party parallactic angle and the second party parallactic angle of the track measuring point that described micro unit is corresponding, the coefficient of friction resistance of well.In concrete embodiment, can establish α 1hole angle for the upper extreme point of track measuring point corresponding to described micro unit, is called the first hole angle, α 2hole angle for the lower extreme point of track measuring point corresponding to described micro unit, is called the second hole angle,
Figure BDA0000426754520000041
azimuth for the upper extreme point of described micro unit, is called first party parallactic angle,
Figure BDA0000426754520000042
azimuth for the lower extreme point of micro unit, is called second party parallactic angle, the coefficient of friction resistance that μ is well.
S105: the axial force of the second end of determining described micro unit according to the coefficient of friction resistance of the moment of inertia of the curvature of described micro unit, length, effective gravity, cross section, modulus of elasticity and the first hole angle, the second hole angle, first party parallactic angle and second party parallactic angle, well is, the lateral force of the axial force of first end, described micro unit unit length.Fig. 2 is the particular flow sheet of step S105.
S106: according to the axial force of the second end of described micro unit, the lateral force of the axial force of first end, described micro unit unit length is determined axial force and the lateral force of three-dimensional well middle sleeve.Step S105 calculates after the axial force and lateral force of micro unit, whole three-dimensional well middle sleeve can be considered by several micro units and forms, and axial force corresponding to several micro units and lateral force stack can be determined to axial force and the lateral force of whole sleeve pipe.In concrete embodiment, also can make the length of sleeve pipe micro unit equal this section of casing string length, step S105 calculates the axial force of micro unit and lateral force and the axial force that lateral force is whole three-dimensional well middle sleeve.
Fig. 2 is the particular flow sheet of the step S105 in Fig. 1, and as shown in Figure 2, step S105 specifically comprises:
S201: determine the axial force of the second end of described micro unit, the relational expression of the axial force of the lateral force of unit length and first end is called the first relational expression according to the coefficient of friction resistance of the moment of inertia of the curvature of described micro unit, length, effective gravity, cross section, modulus of elasticity and the first hole angle, the second hole angle, well.Fig. 3 is the particular flow sheet of step S201, and as shown in Figure 3, step S201 specifically comprises:
S301: determine that according to the curvature of described micro unit and length the full-shape of described micro unit changes.The present invention, in order to set up the universal model that calculates three-dimensional well middle sleeve axial load, first considers a micro unit between two well track measuring points, as shown in figure 11, sets up the relational expression of axial load and other factors.For the ease of deriving, the present invention's hypothesis: the curvature of (1) micro unit is constant (can first try to achieve by minimum-curvature method); (2) casing axis and borehole axis overlap, and the curvature of the implicit micro unit of this hypothesis is identical with hole curvature; Well track between (3) two measuring points is positioned at a space plane; (4) flexural deformation of sleeve pipe is still within elastic range.
The hypothesis that is constant according to the curvature of micro unit, can be according to the curvature of the length of micro unit and micro unit, and the full-shape that is calculated micro unit by following formula changes θ:
θ=KL s (1)
Wherein, K is the curvature of micro unit unit; L slength for micro unit.
S302: determine the lateral force that distortion causes according to the modulus of elasticity of the moment of inertia of the cross section of described micro unit, described micro unit and curvature.If F ethe lateral force causing for casing deformation.By following formula, calculated:
F E=11.3EIK 3 (2)
Wherein, I is the moment of inertia of sleeve pipe cross section; E is the modulus of elasticity of steel; K is the curvature of micro unit.
S303: according to the coefficient of friction resistance of the full-shape variation of described micro unit, length, effective gravity, well, distortion cause lateral force and the first hole angle, the second hole angle determines the axial force of the second end of described micro unit, the relational expression of the axial force of the lateral force of unit length and first end is called the first relational expression.
The effective gravity vector of unit length sleeve pipe is:
q → = q k → - - - ( 3 )
If the axial force of the micro unit lower end i.e. axial force of the second end is T 2, the lateral force F of unit length n, the axial force of its upper end is the axial force T of first end 1can be calculated by following formula:
T 1 = T 2 + L s cos θ 2 [ q cos α ‾ + μ ( F E + F n ) ] - - - ( 4 )
α ‾ = ( α 1 + α 2 ) / 2
Wherein, T 1for the axial force of described micro unit first end, T 2for the axial force of described micro unit the second end, L sfor the length of described micro unit, the full-shape that θ is described micro unit changes, the effective gravity that q is described micro unit,
Figure BDA0000426754520000063
α 1for the first hole angle of track measuring point corresponding to described micro unit, α 2for the second hole angle of track measuring point corresponding to described micro unit, the coefficient of friction resistance that μ is well, gets when roofbolt moves upward "+", gets "-", F when roofbolt moves downward efor the lateral force that distortion causes, F nlateral force for described micro unit unit length.
As shown in Figure 2, step S105 also comprises:
S202: determine the axial force of the second end of described micro unit, the relational expression of the total lateral force in the axial force of first end and full-shape plane is called the second relational expression according to the curvature of described micro unit, length, effective gravity, the first hole angle, the second hole angle, first party parallactic angle and second party parallactic angle.Fig. 4 is the particular flow sheet of step S202, and as shown in Figure 4, step S202 specifically comprises:
S401: determine that according to the curvature of described micro unit and length the full-shape of described micro unit changes.The hypothesis that is constant according to the curvature of micro unit, can be according to the curvature of the length of micro unit and micro unit, and the full-shape that is calculated micro unit by formula (1) changes θ.
S402: tangent vector corresponding to the first end points of determining described micro unit according to the first described angle of slope and first party parallactic angle.According to the hypothesis of the axis of micro unit and borehole axis coincidence, the unit tangent vector of micro unit upper extreme point
Figure BDA0000426754520000064
tangent vector corresponding to the first end points that is micro unit can be expressed as by hole angle and the azimuth of corresponding well track measuring point:
τ → 1 = τ 11 i → + τ 12 j → + τ 13 k → - - - ( 5 a )
Figure BDA0000426754520000066
τ 13=cosα 1 (5d)
Wherein, α 1for i.e. the first hole angle of the hole angle of the upper extreme point of micro unit;
Figure BDA0000426754520000068
the first party parallactic angle that for the azimuth of the upper extreme point of micro unit is; First subscript of tangential component represents the serial number of measuring point; Under second, be designated as: " 1 " represents direct north, " 2 " represent due east direction, and " 3 " represent vertical.
S403: tangent vector corresponding to the second end points of determining described micro unit according to the second described angle of slope and second party parallactic angle.According to the hypothesis of the axis of micro unit and borehole axis coincidence, the unit tangent vector of micro unit lower extreme point
Figure BDA0000426754520000071
tangent vector corresponding to the second end points that is micro unit can be expressed as:
τ → 2 = τ 21 i → + τ 22 j → + τ 23 k → - - - ( 6 a )
Figure BDA0000426754520000074
τ 23=cosα 2 (6d)
Wherein, α 2for i.e. the second hole angle of the hole angle of the lower extreme point of micro unit;
Figure BDA0000426754520000075
for the azimuth of the lower extreme point of micro unit is second party parallactic angle.
S404: will obtain the unit binormal vector of described micro unit after tangent vector corresponding to described the first end points tangent vector multiplication cross corresponding with the second end points, unit.The unit binormal vector of micro unit can be by obtaining after the multiplication cross of the tangent vector of two-end-point unit:
m → = 1 sin θ τ → 1 × τ → 2 = m 1 i → + m 2 j → + m 3 k → - - - ( 7 )
Wherein, the sine that the full-shape of micro unit changes is the sine of micro unit two ends unit tangent vector angle, i.e. mould after two unit tangent vector multiplication crosses.
S405: the unit tangent vector of determining described micro unit mid point according to tangent vector corresponding to described the first end points tangent vector corresponding with the second end points.The unit tangent vector of micro unit mid point is:
τ → 0 = τ → 1 + τ → 2 | τ → 1 + τ → 2 | = τ 01 i → + τ 02 j → + τ 03 k → - - - ( 8 )
S406: described unit binormal vector and described unit tangent vector are carried out to multiplication cross, obtain the unit principal normal vector of described micro unit.The unit principal normal vector of micro unit can be by its unit binormal vector and mid point the multiplication cross of unit tangent vector obtain:
n → = m → × τ → 0 = n 1 i → + n 2 j → + n 3 k → - - - ( 9 a )
Wherein,
n 1=m 2τ 03-m 3τ 02 (9b)
n 2=m 3τ 01-m 1τ 03 (9c)
n 3=m 1τ 02-m 2τ 01 (9d)
S407: determine the axial force of the second end of described micro unit, the relational expression of the total lateral force in the axial force of first end and full-shape plane is called the second relational expression according to the full-shape variation of described micro unit, length, effective gravity, unit principal normal vector.
Total lateral force in full-shape plane is:
F ndp = - ( T 1 + T 2 ) sin θ 2 + L s q → · n → - - - ( 10 a )
Also can be expressed as:
F ndp = - ( T 1 + T 2 ) sin θ 2 + n 3 L s q - - ( 10 b )
Wherein, F ndpfor the total lateral force in full-shape plane, T 1for the axial force of described micro unit first end, T 2for the axial force of described micro unit the second end, L sfor the length of described micro unit, the full-shape that θ is described micro unit changes, for the effective gravity vector of described micro unit,
Figure BDA0000426754520000087
for unit principal normal vector.
As shown in Figure 2, step S105 also comprises:
S203: determine the total lateral force in the binormal direction of described micro unit according to the length of described micro unit, effective gravity, the first hole angle, the second hole angle, first party parallactic angle and second party parallactic angle, be called the 3rd relational expression.Fig. 5 is the particular flow sheet of step S203, and as shown in Figure 5, this step specifically comprises:
S501: tangent vector corresponding to the first end points of determining described micro unit according to the first described angle of slope and first party parallactic angle.According to the hypothesis of the axis of micro unit and borehole axis coincidence, the unit tangent vector of micro unit upper extreme point
Figure BDA0000426754520000088
tangent vector corresponding to the first end points that is micro unit can be represented by hole angle and the azimuth of corresponding well track measuring point, if formula (5a) is to as shown in formula (5d).
S502: tangent vector corresponding to the second end points of determining described micro unit according to the second described angle of slope and second party parallactic angle.According to the hypothesis of the axis of micro unit and borehole axis coincidence, the unit tangent vector of micro unit lower extreme point
Figure BDA0000426754520000089
tangent vector corresponding to the second end points that is micro unit can represent, if formula (6a) is to as shown in formula (6d).
S503: will obtain the unit binormal vector of described micro unit after tangent vector corresponding to described the first end points tangent vector multiplication cross corresponding with the second end points, unit.The unit binormal vector of micro unit can be by obtaining after the multiplication cross of the tangent vector of two-end-point unit, as shown in Equation (7).
S504: determine the total lateral force in the binormal direction of described micro unit according to the length of described micro unit, effective gravity, unit binormal vector, be called the 3rd relational expression.Total lateral force in binormal direction is:
F np = L s q → · m → = m 3 q L s - - - ( 11 a )
Wherein, F npfor the total lateral force in binormal direction, L sfor the length of described micro unit,
Figure BDA0000426754520000084
for the effective gravity vector of described micro unit,
Figure BDA0000426754520000085
for unit binormal vector.
Figure BDA0000426754520000091
As shown in Figure 2, step S105 also comprises:
S204: determine the lateral force of the described micro unit unit length in three-dimensional well according to the total lateral force in total lateral force of described full-shape plane, binormal direction, be called the 4th relational expression.In three-dimensional well, the lateral force of the unit length of a micro unit is total lateral force of full-shape plane and the vector of total lateral force of vertical full-shape plane.Because they are mutually vertical, so it is as follows to obtain the design formulas of lateral force of micro unit unit length:
F n = F ndp 2 + F np 2 L s - - - ( 12 )
S205: the axial force of the second end of determining described micro unit according to described the first relational expression, the second relational expression, the 3rd relational expression, the 4th relational expression is, the lateral force of the axial force of first end, described micro unit unit length.From formula (4) and (10b), if calculate axial force, just must first know lateral force, on the other hand, as calculated lateral force, also must first know axial force, therefore, coupling mutually between lateral force and axial force, because their decoupling zero expression formula is very complicated, so the present invention's available iterative method in concrete embodiment solves, and then the axial force of the second end of determining described micro unit is, the lateral force of the axial force of first end, described micro unit unit length.In concrete embodiment, such as solving by following alternative manner.
Calculate axial force, lateral force that each micro unit is corresponding, such as, 100m, 200m, 300m ... place axial force and lateral force be how many (now element length is 100m).Iterative method concrete steps are as follows:
(1) make sleeve pipe micro unit length equal this section of casing string length.
(2) calculate full-shape variation, hole angle variation, azimuthal variation, average hole angle, mean square parallactic angle, unit normal vector component and the unit binormal vector component in the vertical direction in the vertical direction of micro unit, look into the coefficient of friction resistance of getting this position, unit.
(3) make the axial force of micro unit upper end equal the axial force of its lower end.
(4) lateral force of overlapping pipe range by formula (10b), (11a) and (12) unit of account.
(5) by formula (4), calculated the axial force of micro unit upper end.
(6) lateral force of again overlapping pipe range by formula (10b), (11a) and (12) unit of account.
(7) lateral force of the unit pipe range of relatively calculating in (4) step and (6) step, if their difference is less than permissible value, finishes the iteration of this unit; Otherwise return to (4) step.
Be as mentioned above a kind of method of measuring three-dimensional well middle sleeve axial force and lateral force provided by the invention, the method is to utilize the basic assumption of finite difference method, utilize the method for space coordinate transformation and mathematical analysis, the casing string of take between any two data points is a unit, derived casing string lower unit axial force, unit axial force on tubing string, the axial force that the factors such as tubing string accelerated motion and liquid resistance produce on tubular element, the lateral force that tubing string distortion causes, total lateral force in full-shape plane, total lateral force in binormal direction, the series of parameters equation of the total lateral force of tubular element, axial force and the lateral force of three-dimensional well middle sleeve micro unit are calculated, this method not only can be used for casing strength design, to drilling well, completion, the schematic design of oil recovery and workover treatment and adjustment all have very important reference value, still assess the anti-basis of squeezing the ability of ruining of sleeve pipe.
A kind of equipment of measuring three-dimensional well middle sleeve axial force and lateral force that the present invention also proposes, Fig. 6 is for measuring the structured flowchart of the equipment of three-dimensional well middle sleeve axial force and lateral force, and as shown in Figure 6, described equipment comprises:
Track measuring point data acquisition device 100, for obtaining the track measuring point data of three-dimensional well.
Micro unit selecting device 200, for choosing casing string between any two data points as a micro unit from described track measuring point data.The present invention, in order to set up the universal model that calculates three-dimensional well middle sleeve axial load, first considers a micro unit between two well track measuring points.The hole trajectory data point of take is node, and sleeve pipe is divided into unit, and the sleeve pipe between any two data points is a micro unit, as shown in figure 12.
The first harvester 300, for gathering curvature, length, the effective gravity of described micro unit, the moment of inertia of cross section of described micro unit, the modulus of elasticity of described micro unit.In concrete embodiment, can establish L sfor the length of described micro unit, the curvature that K is micro unit, the effective gravity that q is described micro unit, the moment of inertia of the cross section that I is described micro unit, the modulus of elasticity of the steel that E is described micro unit.
The second harvester 400, for gathering the first hole angle, the second hole angle, first party parallactic angle and the second party parallactic angle of the track measuring point that described micro unit is corresponding, the coefficient of friction resistance of well.In concrete embodiment, can establish α 1hole angle for the upper extreme point of track measuring point corresponding to described micro unit, is called the first hole angle, α 2hole angle for the lower extreme point of track measuring point corresponding to described micro unit, is called the second hole angle, azimuth for the upper extreme point of described micro unit, is called first party parallactic angle,
Figure BDA0000426754520000102
azimuth for the lower extreme point of micro unit, is called second party parallactic angle, the coefficient of friction resistance that μ is well.
Micro unit is to power determinator 500, for the axial force of the second end of determining described micro unit according to the coefficient of friction resistance of the moment of inertia of the curvature of described micro unit, length, effective gravity, cross section, modulus of elasticity and the first hole angle, the second hole angle, first party parallactic angle and second party parallactic angle, well, the lateral force of the axial force of first end, described micro unit unit length.Fig. 7 is that micro unit is to the concrete structure block diagram of power determinator.
Sleeve pipe is to power determinator 600, for according to the axial force of the second end of described micro unit, the lateral force of the axial force of first end, described micro unit unit length is determined axial force and the lateral force of three-dimensional well middle sleeve.Micro unit calculates after the axial force and lateral force of micro unit to power determinator, whole three-dimensional well middle sleeve can be considered by several micro units and forms, and axial force corresponding to several micro units and lateral force stack can be determined to axial force and the lateral force of whole sleeve pipe.In concrete embodiment, also can make the length of sleeve pipe micro unit equal this section of casing string length, micro unit calculates the axial force of micro unit and lateral force and the axial force that lateral force is whole three-dimensional well middle sleeve to power determinator.
Fig. 7 is that in a kind of equipment of measuring three-dimensional well middle sleeve axial force and lateral force provided by the invention, micro unit is to the structured flowchart of power determinator 500, and as shown in Figure 7, micro unit specifically comprises to power determinator 500:
The first relational expression determination module 501, for determining the axial force of the second end of described micro unit, the relational expression of the axial force of the lateral force of unit length and first end is called the first relational expression according to the coefficient of friction resistance of the moment of inertia of the curvature of described micro unit, length, effective gravity, cross section, modulus of elasticity and the first hole angle, the second hole angle, well.Fig. 8 is the structured flowchart of the first relational expression determination module 501, and as shown in Figure 8, the first relational expression determination module 501 specifically comprises:
Full-shape changes determining unit 5011, for determining that according to the curvature of described micro unit and length the full-shape of described micro unit changes.The present invention, in order to set up the universal model that calculates three-dimensional well middle sleeve axial load, first considers a micro unit between two well track measuring points, as shown in figure 11, sets up the relational expression of axial load and other factors.For the ease of deriving, the present invention's hypothesis: the curvature of (1) micro unit is constant (can first try to achieve by minimum-curvature method); (2) casing axis and borehole axis overlap, and the curvature of the implicit micro unit of this hypothesis is identical with hole curvature; Well track between (3) two measuring points is positioned at a space plane; (4) flexural deformation of sleeve pipe is still within elastic range.
The hypothesis that is constant according to the curvature of micro unit, can be according to the curvature of the length of micro unit and micro unit, and the full-shape that is calculated micro unit by following formula changes θ:
θ=KL s (1)
Wherein, the curvature that K is micro unit; L slength for micro unit.
Distortion lateral force determining unit 5012, for determining according to the modulus of elasticity of the moment of inertia of the cross section of described micro unit, described micro unit and curvature the lateral force that distortion causes.If F ethe lateral force causing for casing deformation.By following formula, calculated:
F E=11.3EIK 3 (2)
Wherein, I is the moment of inertia of sleeve pipe cross section; E is the modulus of elasticity of steel; K is the curvature of micro unit.
The first relational expression determining unit 5013, for the lateral force and the tangent vector tangent vector corresponding with the second end points corresponding to the first end points that cause according to the coefficient of friction resistance of the full-shape variation of described micro unit, length, effective gravity, well, distortion, determine the axial force of the second end of described micro unit, the relational expression of the axial force of the lateral force of unit length and first end is called the first relational expression.
The effective gravity vector of unit length sleeve pipe is:
q → = q k → - - - ( 3 )
If the axial force of the micro unit lower end i.e. axial force of the second end is T 2, the lateral force F of unit length n, the axial force of its upper end is the axial force T of first end 1can be calculated by following formula:
T 1 = T 2 + L s cos θ 2 [ q cos α ‾ + μ ( F E + F n ) ] - - - ( 4 )
α ‾ = ( α 1 + α 2 ) / 2
Wherein, T 1for the axial force of described micro unit first end, T 2for the axial force of described micro unit the second end, the length that Ls is described micro unit, the full-shape that θ is described micro unit changes, the effective gravity that q is described micro unit, α 1for the first hole angle of track measuring point corresponding to described micro unit, α 2for the second hole angle of track measuring point corresponding to described micro unit, the coefficient of friction resistance that μ is well, gets when roofbolt moves upward "+", gets "-", F when roofbolt moves downward efor the lateral force that distortion causes, F nlateral force for described micro unit unit length.
As shown in Figure 7, micro unit also comprises to power determinator 500:
The second relational expression determination module 502, for determining the axial force of the second end of described micro unit, the relational expression of the total lateral force in the axial force of first end and full-shape plane is called the second relational expression according to the curvature of described micro unit, length, effective gravity, the first hole angle, the second hole angle, first party parallactic angle and second party parallactic angle.Fig. 9 is the structured flowchart of the second relational expression determination module 502, and as shown in Figure 9, the second relational expression determination module 502 specifically comprises:
Full-shape changes determining unit 5021, for determining that according to the curvature of described micro unit and length the full-shape of described micro unit changes.The hypothesis that is constant according to the curvature of micro unit, can be according to the curvature of the length of micro unit and micro unit, and the full-shape that is calculated micro unit by formula (1) changes θ.
The first tangent vector determining unit 5022, for determining the tangent vector corresponding to the first end points of described micro unit according to the first described angle of slope and first party parallactic angle.According to the hypothesis of the axis of micro unit and borehole axis coincidence, the unit tangent vector of micro unit upper extreme point tangent vector corresponding to the first end points that is micro unit can be expressed as by hole angle and the azimuth of corresponding well track measuring point:
τ → 1 = τ 11 i → + τ 12 j → + τ 13 k → - - - ( 5 a )
Figure BDA0000426754520000127
τ 13=cosα 1 (5d)
Wherein, α 1for i.e. the first hole angle of the hole angle of the upper extreme point of micro unit;
Figure BDA0000426754520000128
the first party parallactic angle that for the azimuth of the upper extreme point of micro unit is; First subscript of tangential component represents the serial number of measuring point; Under second, be designated as: " 1 " represents direct north, " 2 " represent due east direction, and " 3 " represent vertical.The second tangent vector determining unit 5023, for determining the tangent vector corresponding to the second end points of described micro unit according to the second described angle of slope and second party parallactic angle.According to the hypothesis of the axis of micro unit and borehole axis coincidence, the unit tangent vector of micro unit lower extreme point
Figure BDA0000426754520000131
tangent vector corresponding to the second end points that is micro unit can be expressed as:
τ → 2 = τ 21 i → + τ 22 j → + τ 23 k → - - - ( 6 a )
Figure BDA0000426754520000133
Figure BDA0000426754520000134
τ 23=cosα 2 (6d)
Wherein, α 2for i.e. the second hole angle of the hole angle of the lower extreme point of micro unit;
Figure BDA0000426754520000135
for the azimuth of the lower extreme point of micro unit is second party parallactic angle.
Binormal vector determining unit 5024, for obtaining the unit binormal vector of described micro unit after tangent vector corresponding to described the first end points tangent vector multiplication cross corresponding with the second end points, unit.The unit binormal vector of micro unit can be by obtaining after the multiplication cross of the tangent vector of two-end-point unit:
m → = 1 sin θ τ → 1 × τ → 2 = m 1 i → + m 2 j → + m 3 k → - - - ( 7 )
Wherein, the sine that the full-shape of micro unit changes is the sine of micro unit two ends unit tangent vector angle, i.e. mould after two unit tangent vector multiplication crosses.
Unit tangent vector determining unit 5025, for determining the unit tangent vector of described micro unit mid point according to tangent vector corresponding to described the first end points tangent vector corresponding with the second end points.The unit tangent vector of micro unit mid point is:
τ → 0 = τ → 1 + τ → 2 | τ → 1 + τ → 2 | = τ 01 i → + τ 02 j → + τ 03 k → - - - ( 8 )
Principal normal vector determining unit 5026, for described unit binormal vector and described unit tangent vector are carried out to multiplication cross, obtains the unit principal normal vector of described micro unit.The unit principal normal vector of micro unit can be by its unit binormal vector and mid point the multiplication cross of unit tangent vector obtain:
n → = m → × τ → 0 = n 1 i → + n 2 j → + n 3 k → - - - ( 9 a )
Wherein,
n 1=m 2τ 03-m 3τ 02 (9b)
n 2=m 3τ 01-m 1τ 03 (9c)
n 3=m 1τ 02-m 2τ 01 (9d)
The second relational expression determining unit 5027, for determining the axial force of the second end of described micro unit, the relational expression of the total lateral force in the axial force of first end and full-shape plane is called the second relational expression according to the full-shape variation of described micro unit, length, effective gravity, unit principal normal vector.
Total lateral force in full-shape plane is:
F ndp = - ( T 1 + T 2 ) sin θ 2 + L s q → · n → - - - ( 10 a )
Also can be expressed as:
F ndp = - ( T 1 + T 2 ) sin θ 2 + n 3 L s q - - ( 10 b )
Wherein, F ndpfor the total lateral force in full-shape plane, T 1for the axial force of described micro unit first end, T 2for the axial force of described micro unit the second end, L sfor the length of described micro unit, the full-shape that θ is described micro unit changes,
Figure BDA0000426754520000143
for the effective gravity vector of described micro unit,
Figure BDA0000426754520000144
for unit principal normal vector.
As shown in Figure 7, micro unit also comprises to power determinator 500:
The 3rd relational expression determination module 503, for determining according to the length of described micro unit, effective gravity, the first hole angle, the second hole angle, first party parallactic angle and second party parallactic angle, the total lateral force in the binormal direction of described micro unit be called the 3rd relational expression.Figure 10 is the structured flowchart of the 3rd relational expression determination module 503, and as shown in Figure 10, the 3rd relational expression determination module 503 specifically comprises:
The first tangent vector determining unit 5031, for determining the tangent vector corresponding to the first end points of described micro unit according to the first described angle of slope and first party parallactic angle.According to the hypothesis of the axis of micro unit and borehole axis coincidence, the unit tangent vector of micro unit upper extreme point
Figure BDA0000426754520000145
tangent vector corresponding to the first end points that is micro unit can be represented by hole angle and the azimuth of corresponding well track measuring point, if formula (5a) is to as shown in formula (5d).
The second tangent vector determining unit 5032, for determining the tangent vector corresponding to the second end points of described micro unit according to the second described angle of slope and second party parallactic angle.According to the hypothesis of the axis of micro unit and borehole axis coincidence, the unit tangent vector of micro unit lower extreme point
Figure BDA0000426754520000146
tangent vector corresponding to the second end points that is micro unit can represent, if formula (6a) is to as shown in formula (6d).
Binormal vector determining unit 5033, for obtaining the unit binormal vector of described micro unit after tangent vector corresponding to described the first end points tangent vector multiplication cross corresponding with the second end points, unit.The unit binormal vector of micro unit can be by obtaining after the multiplication cross of the tangent vector of two-end-point unit, as shown in Equation (7).
The 3rd relational expression determining unit 5034, for determine the total lateral force in the binormal direction of described micro unit according to the length of described micro unit, effective gravity, unit binormal vector, is called the 3rd relational expression.Total lateral force in binormal direction is:
F np = L s q → · m → = m 3 q L s - - - ( 11 a )
Wherein, F npfor the total lateral force in binormal direction, L sfor the length of described micro unit,
Figure BDA0000426754520000148
for the effective gravity vector of described micro unit,
Figure BDA0000426754520000149
for unit binormal vector.
Figure BDA0000426754520000151
As shown in Figure 7, micro unit also comprises to power determinator 500:
The 4th relational expression determination module 504, the lateral force of the described micro unit unit length of three-dimensional well is called the 4th relational expression for determining according to the total lateral force in total lateral force of described full-shape plane, binormal direction.In three-dimensional well, the lateral force of the unit length of a micro unit is total lateral force of full-shape plane and the vector of total lateral force of vertical full-shape plane.Because they are mutually vertical, so it is as follows to obtain the design formulas of lateral force of micro unit unit length:
F n = F ndp 2 + F np 2 L s - - - ( 12 )
Micro unit is measured module 505 to power, for the axial force of the second end of determining described micro unit according to described the first relational expression, the second relational expression, the 3rd relational expression, the 4th relational expression, the lateral force of the axial force of first end, described micro unit unit length.From formula (4) and (10b), if calculate axial force, just must first know lateral force, on the other hand, as calculated lateral force, also must first know axial force, therefore, coupling mutually between lateral force and axial force, because their decoupling zero expression formula is very complicated, so the present invention's available iterative method in concrete embodiment solves, and then the axial force of the second end of determining described micro unit is, the lateral force of the axial force of first end, described micro unit unit length.In concrete embodiment, such as solving by following alternative manner.
Calculate axial force, lateral force that each micro unit is corresponding, such as, 100m, 200m, 300m ... place axial force and lateral force be how many (now element length is 100m).Iterative method concrete steps are as follows:
(1) make sleeve pipe micro unit length equal this section of casing string length.
(2) calculate full-shape variation, hole angle variation, azimuthal variation, average hole angle, mean square parallactic angle, unit normal vector component and the unit binormal vector component in the vertical direction in the vertical direction of micro unit, look into the coefficient of friction resistance of getting this position, unit.
(3) make the axial force of micro unit upper end equal the axial force of its lower end.
(4) lateral force of overlapping pipe range by formula (10b), (11a) and (12) unit of account.
(5) by formula (4), calculated the axial force of micro unit upper end.
(6) lateral force of again overlapping pipe range by formula (10b), (11a) and (12) unit of account.
(7) lateral force of the unit pipe range of relatively calculating in (4) step and (6) step, if their difference is less than permissible value, finishes the iteration of this unit; Otherwise return to (4) step.
Be as mentioned above a kind of equipment of measuring three-dimensional well middle sleeve axial force and lateral force provided by the invention, the present invention is calculated axial force and the lateral force of three-dimensional well middle sleeve micro unit by finite difference method, this programme not only can be used for casing strength design, the schematic design of drilling well, completion, oil recovery and workover treatment and adjustment are all had to very important reference value, still assess the anti-basis of squeezing the ability of ruining of sleeve pipe.
Below in conjunction with specific embodiment, introduce in detail technical scheme of the present invention.The hole trajectory data point of take is node, and sleeve pipe roofbolt is divided into unit, and the sleeve pipe roofbolt between any two data points is a unit, as shown in figure 12.Because whole casing string may be comprised of the sleeve pipe of different model, the upper and lower side position of casing section may, between 2 hole trajectory data points, therefore, need to increase node in casing section boundary position.The calculation process of axial force and the lateral force of one section of micro unit below.
(1) by the depth measurement of casing section top, from well track module, fetch corresponding track data point sequence number (top), hole angle and azimuth; By casing section bottom depth measurement, from well track module, fetch corresponding track data point sequence number (bottom), hole angle and azimuth.
(2) if track data point sequence number corresponding to the casing section top track data point sequence number corresponding with casing section bottom is identical, enters next step, otherwise forward (10) step to.
(3) make sleeve pipe micro unit length equal this section of casing string length.
(4) calculate full-shape variation, hole angle variation, azimuthal variation, average hole angle, mean square parallactic angle, unit normal vector component and the unit binormal vector component in the vertical direction in the vertical direction of micro unit, look into the coefficient of friction resistance of getting this position, unit.
(5) make the axial force of micro unit upper end equal the axial force of its lower end.
(6) by the lateral force of formula (10b), (11a) and (12) unit of account pipe range.
(7) by formula (6), calculated the axial force of micro unit upper end.
(8) again by the lateral force of formula (10b), (11a) and (12) unit of account pipe range.
(9) lateral force of the unit pipe range of relatively calculating in (6) step and (8) step, if their difference is less than permissible value, finishes the iteration of this unit; Otherwise return to (6) step.
(10) this section of tubing string is divided into (bottom-top+1) individual unit, tubular element calculates from (bottom+1) to (top+1) circulation, and cyclic variable is KU, and increment step-length is-1.
(11) if KU equals (bottom+1), tubular element is most on the lower one, and track data point sequence number corresponding to upper end, unit is bottom, and the hole trajectory data of lower end obtains by interpolation; If KU equals (top+1), tubular element is the most top one, and track data point sequence number corresponding to lower end, unit is (top+1), and the hole trajectory data of upper end obtains by interpolation; If KU is between (bottom+1) with (top+1), track data point sequence number corresponding to upper end, unit is (KU-1), and that lower end is corresponding is KU.
(12) return to execution step for (4) to (9) step.
In sum, useful achievement of the present invention is: a kind of method and apparatus of measuring three-dimensional well middle sleeve axial force and lateral force is provided, utilize the basic assumption of finite difference method, utilize the method for space coordinate transformation and mathematical analysis, the casing string of take between any two data points is a unit, derived casing string lower unit axial force, unit axial force on tubing string, the axial force that the factors such as tubing string accelerated motion and liquid resistance produce on tubular element, the lateral force that tubing string distortion causes, total lateral force in full-shape plane, total lateral force in binormal direction, the series of parameters equation of the total lateral force of tubular element, utilize finite difference method to propose a kind of approximate, easy scheme, based on hole trajectory data, sleeve pipe is divided into the micro unit section of unit length, by sleeve pipe lowermost end, started to calculate piecemeal from lower to upper, meet petroleum production engineering requirement, not only can be used for casing strength design, to drilling well, completion, the schematic design of oil recovery and workover treatment and adjustment all have very important reference value, still assess the anti-basis of squeezing the ability of ruining of sleeve pipe.
One of ordinary skill in the art will appreciate that all or part of flow process realizing in above-described embodiment method, can come the hardware that instruction is relevant to complete by computer program, described program can be stored in general computer read/write memory medium, this program, when carrying out, can comprise as the flow process of the embodiment of above-mentioned each side method.Wherein, described storage medium can be magnetic disc, CD, read-only store-memory body (Read-Only Memory, ROM) or random store-memory body (Random Access Memory, RAM) etc.
Those skilled in the art can also recognize that the various functions that the embodiment of the present invention is listed are to realize by hardware or software the designing requirement of depending on specific application and whole system.Those skilled in the art can, for every kind of specific application, can make in all sorts of ways and realize described function, but this realization should not be understood to exceed the scope of embodiment of the present invention protection.
In the present invention, applied specific embodiment principle of the present invention and embodiment are set forth, the explanation of above embodiment is just for helping to understand method of the present invention and core concept thereof; , for one of ordinary skill in the art, according to thought of the present invention, all will change in specific embodiments and applications, in sum, this description should not be construed as limitation of the present invention meanwhile.

Claims (18)

1. a method of measuring three-dimensional well middle sleeve axial force and lateral force, is characterized in that, described method specifically comprises:
Obtain the track measuring point data of three-dimensional well;
From described track measuring point data, choose casing string between any two data points as a micro unit;
Gather curvature, length, the effective gravity of described micro unit, the moment of inertia of cross section of described micro unit, the modulus of elasticity of described micro unit;
Gather the first hole angle, the second hole angle, first party parallactic angle and the second party parallactic angle of the track measuring point that described micro unit is corresponding, the coefficient of friction resistance of well;
The axial force of the second end of determining described micro unit according to the coefficient of friction resistance of the moment of inertia of the curvature of described micro unit, length, effective gravity, cross section, modulus of elasticity and the first hole angle, the second hole angle, first party parallactic angle and second party parallactic angle, well is, the lateral force of the axial force of first end, described micro unit unit length;
According to the axial force of the second end of described micro unit, the lateral force of the axial force of first end, described micro unit unit length is determined axial force and the lateral force of three-dimensional well middle sleeve.
2. method according to claim 1, it is characterized in that, according to the coefficient of friction resistance of the moment of inertia of the curvature of described micro unit, length, effective gravity, cross section, modulus of elasticity and the first hole angle, the second hole angle, first party parallactic angle and second party parallactic angle, well, determine the axial force of the second end of described micro unit, the lateral force of the axial force of first end, described micro unit unit length specifically comprises:
According to the coefficient of friction resistance of the moment of inertia of the curvature of described micro unit, length, effective gravity, cross section, modulus of elasticity and the first hole angle, the second hole angle, well, determine the axial force of the second end of described micro unit, the relational expression of the axial force of the lateral force of unit length and first end is called the first relational expression;
According to the curvature of described micro unit, length, effective gravity, the first hole angle, the second hole angle, first party parallactic angle and second party parallactic angle, determine the axial force of the second end of described micro unit, the relational expression of the total lateral force in the axial force of first end and full-shape plane is called the second relational expression;
According to the length of described micro unit, effective gravity, the first hole angle, the second hole angle, first party parallactic angle and second party parallactic angle, determine the total lateral force in the binormal direction of described micro unit, be called the 3rd relational expression;
According to the total lateral force in total lateral force of described full-shape plane, binormal direction, determine the lateral force of the described micro unit unit length in three-dimensional well, be called the 4th relational expression;
The axial force of the second end of determining described micro unit according to described the first relational expression, the second relational expression, the 3rd relational expression, the 4th relational expression is, the lateral force of the axial force of first end, described micro unit unit length.
3. method according to claim 2, it is characterized in that, according to the coefficient of friction resistance of the moment of inertia of the curvature of described micro unit, length, effective gravity, cross section, modulus of elasticity and the first hole angle, the second hole angle, well, determine the axial force of the second end of described micro unit, the relational expression of the axial force of the lateral force of unit length and first end is called the first relational expression and specifically comprises:
According to the curvature of described micro unit and length, determine that the full-shape of described micro unit changes;
According to the modulus of elasticity of the moment of inertia of the cross section of described micro unit, described micro unit and curvature, determine the lateral force that distortion causes;
According to the coefficient of friction resistance of the full-shape variation of described micro unit, length, effective gravity, well, distortion cause lateral force and the first hole angle, the second hole angle determines the axial force of the second end of described micro unit, the relational expression of the axial force of the lateral force of unit length and first end is called the first relational expression.
4. method according to claim 3, is characterized in that, the first described relational expression is:
T 1 = T 2 + L s cos θ 2 [ q cos α ‾ + μ ( F E + F n ) ]
Wherein, T 1for the axial force of described micro unit first end, T 2for the axial force of described micro unit the second end, L sfor the length of described micro unit, the full-shape that θ is described micro unit changes, the effective gravity that q is described micro unit,
Figure FDA0000426754510000022
α 1for the first hole angle of track measuring point corresponding to described micro unit, α 2for the second hole angle of track measuring point corresponding to described micro unit, the coefficient of friction resistance that μ is well, F efor the lateral force that distortion causes, F nlateral force for described micro unit unit length.
5. according to the method described in claim 2 or 4, it is characterized in that, according to the curvature of described micro unit, length, effective gravity, the first hole angle, the second hole angle, first party parallactic angle and second party parallactic angle, determine the axial force of the second end of described micro unit, the relational expression of the total lateral force in the axial force of first end and full-shape plane is called the second relational expression and specifically comprises:
According to the curvature of described micro unit and length, determine that the full-shape of described micro unit changes;
According to the first described angle of slope and first party parallactic angle, determine the tangent vector corresponding to the first end points of described micro unit;
According to the second described angle of slope and second party parallactic angle, determine the tangent vector corresponding to the second end points of described micro unit;
The unit binormal vector of described micro unit will be obtained after tangent vector corresponding to described the first end points tangent vector multiplication cross corresponding with the second end points, unit;
According to tangent vector corresponding to described the first end points tangent vector corresponding with the second end points, determine the unit tangent vector of described micro unit mid point;
Described unit binormal vector and described unit tangent vector are carried out to multiplication cross, obtain the unit principal normal vector of described micro unit;
According to the full-shape variation of described micro unit, length, effective gravity, unit principal normal vector, determine the axial force of the second end of described micro unit, the relational expression of the total lateral force in the axial force of first end and full-shape plane is called the second relational expression.
6. method according to claim 5, is characterized in that, the second described relational expression is:
F ndp = - ( T 1 + T 2 ) sin θ 2 + L s q → · n →
Wherein, F ndpfor the total lateral force in full-shape plane, T 1for the axial force of described micro unit first end, T 2for the axial force of described micro unit the second end, L sfor the length of described micro unit, the full-shape that θ is described micro unit changes,
Figure FDA0000426754510000032
for the effective gravity vector of described micro unit,
Figure FDA0000426754510000033
for unit principal normal vector.
7. method according to claim 6, it is characterized in that, according to the length of described micro unit, effective gravity, the first hole angle, the second hole angle, first party parallactic angle and second party parallactic angle, determine the total lateral force in the binormal direction of described micro unit, be called the 3rd relational expression and specifically comprise:
According to the first described angle of slope and first party parallactic angle, determine the tangent vector corresponding to the first end points of described micro unit;
According to the second described angle of slope and second party parallactic angle, determine the tangent vector corresponding to the second end points of described micro unit;
The unit binormal vector of described micro unit will be obtained after tangent vector corresponding to described the first end points tangent vector multiplication cross corresponding with the second end points, unit;
According to the length of described micro unit, effective gravity, unit binormal vector, determine the total lateral force in the binormal direction of described micro unit, be called the 3rd relational expression.
8. method according to claim 7, is characterized in that, the 3rd described relational expression is:
F np = L s q → · m →
Wherein, F npfor the total lateral force in binormal direction, L sfor the length of described micro unit,
Figure FDA0000426754510000035
for the effective gravity vector of described micro unit,
Figure FDA0000426754510000036
for unit binormal vector.
9. method according to claim 8, is characterized in that, the 4th described relational expression is:
F n = F ndp 2 + F np 2 L s .
10. an equipment of measuring three-dimensional well middle sleeve axial force and lateral force, is characterized in that, described equipment specifically comprises:
Track measuring point data acquisition device, for obtaining the track measuring point data of three-dimensional well;
Micro unit selecting device, for choosing casing string between any two data points as a micro unit from described track measuring point data;
The first harvester, for gathering curvature, length, the effective gravity of described micro unit, the moment of inertia of cross section of described micro unit, the modulus of elasticity of described micro unit;
The second harvester, for gathering the first hole angle, the second hole angle, first party parallactic angle and the second party parallactic angle of the track measuring point that described micro unit is corresponding, the coefficient of friction resistance of well;
Micro unit is to power determinator, for the axial force of the second end of determining described micro unit according to the coefficient of friction resistance of the moment of inertia of the curvature of described micro unit, length, effective gravity, cross section, modulus of elasticity and the first hole angle, the second hole angle, first party parallactic angle and second party parallactic angle, well, the lateral force of the axial force of first end, described micro unit unit length;
Sleeve pipe is to power determinator, for according to the axial force of the second end of described micro unit, the lateral force of the axial force of first end, described micro unit unit length is determined axial force and the lateral force of three-dimensional well middle sleeve.
11. equipment according to claim 10, is characterized in that, described micro unit specifically comprises to power determinator:
The first relational expression determination module, for determining the axial force of the second end of described micro unit, the relational expression of the axial force of the lateral force of unit length and first end is called the first relational expression according to the coefficient of friction resistance of the moment of inertia of the curvature of described micro unit, length, effective gravity, cross section, modulus of elasticity and the first hole angle, the second hole angle, well;
The second relational expression determination module, for determining the axial force of the second end of described micro unit, the relational expression of the total lateral force in the axial force of first end and full-shape plane is called the second relational expression according to the curvature of described micro unit, length, effective gravity, the first hole angle, the second hole angle, first party parallactic angle and second party parallactic angle;
The 3rd relational expression determination module, for determine the total lateral force in the binormal direction of described micro unit according to the length of described micro unit, effective gravity, the first hole angle, the second hole angle, first party parallactic angle and second party parallactic angle, is called the 3rd relational expression;
The 4th relational expression determination module, the lateral force of the described micro unit unit length of three-dimensional well is called the 4th relational expression for determining according to the total lateral force in total lateral force of described full-shape plane, binormal direction;
Micro unit is measured module to power, for the axial force of the second end of determining described micro unit according to described the first relational expression, the second relational expression, the 3rd relational expression, the 4th relational expression, the lateral force of the axial force of first end, described micro unit unit length.
12. equipment according to claim 11, is characterized in that, the first described relational expression determination module specifically comprises:
Full-shape changes determining unit, for determining that according to the curvature of described micro unit and length the full-shape of described micro unit changes;
Distortion lateral force determining unit, for determining according to the modulus of elasticity of the moment of inertia of the cross section of described micro unit, described micro unit and curvature the lateral force that distortion causes;
The first relational expression determining unit, for according to the coefficient of friction resistance of the full-shape variation of described micro unit, length, effective gravity, well, distortion cause lateral force and the first hole angle, the second hole angle determines the axial force of the second end of described micro unit, the relational expression of the axial force of the lateral force of unit length and first end is called the first relational expression.
13. equipment according to claim 12, is characterized in that, the first described relational expression is:
T 1 = T 2 + L s cos θ 2 [ q cos α ‾ + μ ( F E + F n ) ]
Wherein, T 1for the axial force of described micro unit first end, T 2for the axial force of described micro unit the second end, L sfor the length of described micro unit, the full-shape that θ is described micro unit changes, the effective gravity that q is described micro unit,
Figure FDA0000426754510000052
α 1for the first hole angle of track measuring point corresponding to described micro unit, α 2for the second hole angle of track measuring point corresponding to described micro unit, the coefficient of friction resistance that μ is well, F efor the lateral force that distortion causes, F nlateral force for described micro unit unit length.
14. according to the equipment described in claim 11 or 13, it is characterized in that, the second described relational expression determination module specifically comprises:
Full-shape changes determining unit, for determining that according to the curvature of described micro unit and length the full-shape of described micro unit changes;
The first tangent vector determining unit, for determining the tangent vector corresponding to the first end points of described micro unit according to the first described angle of slope and first party parallactic angle;
The second tangent vector determining unit, for determining the tangent vector corresponding to the second end points of described micro unit according to the second described angle of slope and second party parallactic angle;
Binormal vector determining unit, for obtaining the unit binormal vector of described micro unit after tangent vector corresponding to described the first end points tangent vector multiplication cross corresponding with the second end points, unit;
Unit tangent vector determining unit, for determining the unit tangent vector of described micro unit mid point according to tangent vector corresponding to described the first end points tangent vector corresponding with the second end points;
Principal normal vector determining unit, for described unit binormal vector and described unit tangent vector are carried out to multiplication cross, obtains the unit principal normal vector of described micro unit;
The second relational expression determining unit, for determining the axial force of the second end of described micro unit, the relational expression of the total lateral force in the axial force of first end and full-shape plane is called the second relational expression according to the full-shape variation of described micro unit, length, effective gravity, unit principal normal vector.
15. equipment according to claim 14, is characterized in that, the second described relational expression is:
F ndp = - ( T 1 + T 2 ) sin θ 2 + L s q → · n →
Wherein, F ndpfor the total lateral force in full-shape plane, T 1for the axial force of described micro unit first end, T 2for the axial force of described micro unit the second end, L sfor the length of described micro unit, the full-shape that θ is described micro unit changes,
Figure FDA0000426754510000054
for the effective gravity vector of described micro unit,
Figure FDA0000426754510000055
for unit principal normal vector.
16. equipment according to claim 15, is characterized in that, the 3rd described relational expression determination module specifically comprises:
The first tangent vector determining unit, for determining the tangent vector corresponding to the first end points of described micro unit according to the first described angle of slope and first party parallactic angle;
The second tangent vector determining unit, for determining the tangent vector corresponding to the second end points of described micro unit according to the second described angle of slope and second party parallactic angle;
Binormal vector determining unit, for obtaining the unit binormal vector of described micro unit after tangent vector corresponding to described the first end points tangent vector multiplication cross corresponding with the second end points, unit;
The 3rd relational expression determining unit, for determine the total lateral force in the binormal direction of described micro unit according to the length of described micro unit, effective gravity, unit binormal vector, is called the 3rd relational expression.
17. equipment according to claim 16, is characterized in that, the 3rd described relational expression is:
F np = L s q → · m →
Wherein, F npfor the total lateral force in binormal direction, L sfor the length of described micro unit,
Figure FDA0000426754510000063
for the effective gravity vector of described micro unit,
Figure FDA0000426754510000064
for unit binormal vector.
18. equipment according to claim 17, is characterized in that, the 4th described relational expression is:
F n = F ndp 2 + F np 2 L s .
CN201310629634.9A 2013-11-29 2013-11-29 Method and equipment for measuring axial force and lateral force of casing in three-dimensional borehole Pending CN103615233A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310629634.9A CN103615233A (en) 2013-11-29 2013-11-29 Method and equipment for measuring axial force and lateral force of casing in three-dimensional borehole

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310629634.9A CN103615233A (en) 2013-11-29 2013-11-29 Method and equipment for measuring axial force and lateral force of casing in three-dimensional borehole

Publications (1)

Publication Number Publication Date
CN103615233A true CN103615233A (en) 2014-03-05

Family

ID=50165947

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310629634.9A Pending CN103615233A (en) 2013-11-29 2013-11-29 Method and equipment for measuring axial force and lateral force of casing in three-dimensional borehole

Country Status (1)

Country Link
CN (1) CN103615233A (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104879115A (en) * 2015-05-15 2015-09-02 中国石油天然气集团公司 Method and device for determining underground measurement-while-drilling parameter
CN106548030A (en) * 2016-11-03 2017-03-29 中国石油天然气股份有限公司 Multi-packer composite pipe column stress analysis and calculation method
CN108397130A (en) * 2018-02-02 2018-08-14 中国石油天然气集团有限公司 Boring method and device
CN112100796A (en) * 2019-06-18 2020-12-18 中国石油天然气集团有限公司 Drilling track determination method and device for preventing casing in shale gas well from being damaged

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5464058A (en) * 1993-01-25 1995-11-07 James N. McCoy Method of using a polished rod transducer
RU2221227C1 (en) * 2002-04-17 2004-01-10 Ооо "Нпп "Грант" Strain-measuring force cell for dynamometer testing of well sucker-rod pumps
CN201292838Y (en) * 2008-11-24 2009-08-19 中国石油集团长城钻探工程有限公司 Pipe column state detection device
CN102778311A (en) * 2012-07-23 2012-11-14 北京航空航天大学 Multidimensional radial pressure measuring device of underground sucker rod
CN103616117A (en) * 2013-11-29 2014-03-05 中国石油天然气股份有限公司 Method and equipment for measuring axial force and lateral force of sucker rod string in three-dimensional well

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5464058A (en) * 1993-01-25 1995-11-07 James N. McCoy Method of using a polished rod transducer
RU2221227C1 (en) * 2002-04-17 2004-01-10 Ооо "Нпп "Грант" Strain-measuring force cell for dynamometer testing of well sucker-rod pumps
CN201292838Y (en) * 2008-11-24 2009-08-19 中国石油集团长城钻探工程有限公司 Pipe column state detection device
CN102778311A (en) * 2012-07-23 2012-11-14 北京航空航天大学 Multidimensional radial pressure measuring device of underground sucker rod
CN103616117A (en) * 2013-11-29 2014-03-05 中国石油天然气股份有限公司 Method and equipment for measuring axial force and lateral force of sucker rod string in three-dimensional well

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
覃成锦: "《油气井套管柱载荷分析及优化设计研究》", 《中国博士学位论文全文数据库 工程科技I辑》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104879115A (en) * 2015-05-15 2015-09-02 中国石油天然气集团公司 Method and device for determining underground measurement-while-drilling parameter
CN104879115B (en) * 2015-05-15 2017-09-01 中国石油天然气集团公司 A kind of downhole drill determination method for parameter and device
CN106548030A (en) * 2016-11-03 2017-03-29 中国石油天然气股份有限公司 Multi-packer composite pipe column stress analysis and calculation method
CN106548030B (en) * 2016-11-03 2019-06-11 中国石油天然气股份有限公司 multi-packer composite pipe column stress analysis and calculation method
CN108397130A (en) * 2018-02-02 2018-08-14 中国石油天然气集团有限公司 Boring method and device
CN108397130B (en) * 2018-02-02 2019-11-08 中国石油天然气集团有限公司 Boring method and device
CN112100796A (en) * 2019-06-18 2020-12-18 中国石油天然气集团有限公司 Drilling track determination method and device for preventing casing in shale gas well from being damaged
CN112100796B (en) * 2019-06-18 2024-05-28 中国石油天然气集团有限公司 Drilling track determination method and device for preventing casing in shale gas well from being damaged

Similar Documents

Publication Publication Date Title
CN103616117B (en) Method and equipment for measuring axial force and lateral force of sucker rod string in three-dimensional well
CN103226641B (en) Coupling calculation method of deepwater gas-liquid two-phase flow circulating temperature and pressure
CN102182453B (en) Borehole wall sloughing analysis method
CN103590818B (en) Method and device for semi-analytically determining temperature distribution in shaft
CN103615233A (en) Method and equipment for measuring axial force and lateral force of casing in three-dimensional borehole
CN109033504B (en) Oil-water well casing damage prediction method
CN102682195B (en) Semisubmersible platform transient state bored shaft temperature computation method
CN102434119B (en) Control method of Casson model drilling fluid rheological parameters
CN103726815B (en) A kind of CO 2drive produced well pit shaft fluidised form is determined and parameter optimization method
Shook et al. Use of tracers and temperature to estimate fracture surface area for EGS reservoirs
CN103883267B (en) Arrangement method and device of centralizer for drill rod
CN104389589A (en) Method and system for determining wellbore temperature field distribution based on hollow rod
CN104747152A (en) Heavy oil deposit multi-element hot fluid throughput cross-well gas channeling degree detection method
CN110298057B (en) Supercritical carbon dioxide fracturing fracture extension calculation method
CN106639879A (en) Method for predicting build-up rate of bottom hole assembly
CN105178930A (en) Method and device for calculating steam thermal parameters in ground steam injection pipeline
CN113109162B (en) Rock fracture initiation pressure calculation method based on thermo-fluid-solid coupling
Zhang et al. The enhanced geothermal system heat mining prediction based on fracture propagation simulation of thermo-hydro-mechanical-damage coupling: Insight from the integrated research of heat mining and supercritical CO2 fracturing
CN106018472A (en) Method and device for testing heat conductivity coefficient of heat insulation pipe of steam injection well
CN116167302A (en) Description method of artificial complex cracks in natural gas hydrate yield increase simulation
CN110929447A (en) Shaft temperature field numerical calculation method in thickening acid acidification process
CN105403590A (en) Heat conductivity coefficient testing method and device for heat insulation pipe
CN114282387A (en) DTS-based comprehensive evaluation method for steam injection horizontal well of heavy oil reservoir
CN108090320A (en) A kind of round platform screw type energy pile Heat Transfer Calculation
CN109784003A (en) A kind of optimum organization method based on when repairing horizontal well tubing string

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20140305