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

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 α ₁hole angle for the upper extreme point of track measuring point corresponding to described micro unit, is called the first hole angle, α ₂hole 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,

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 ³ (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:

$\stackrel{\→}{q}=q\stackrel{\→}{k}---\left(3\right)$

If the axial force of the micro unit lower end i.e. axial force of the second end is T ₂, the lateral force F of unit length _n, the axial force of its upper end is the axial force T of first end ₁can be calculated by following formula:

${\mathrm{<figure-callout\; id="1"\; label="T"\; filenames="HDA0000426754530000112.png"\; state="\{\{state\}\}">T</figure-callout>}}_1=T_2+\frac{L_s}{\cos \frac{\mathrm{\&theta;}}{2}}[q\cos \stackrel{\&OverBar;}{\mathrm{\&alpha;}}+\mathrm{\&mu;}(F_E+F_n)]---\left(4\right)$

$\stackrel{\&OverBar;}{\mathrm{\&alpha;}}=({\mathrm{\&alpha;}}_1+{\mathrm{\&alpha;}}_2)/2$

Wherein, T ₁for the axial force of described micro unit first end, T ₂for 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,

α ₁for the first hole angle of track measuring point corresponding to described micro unit, α ₂for 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

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:

${\stackrel{\&RightArrow;}{\mathrm{\&tau;}}}_1={\mathrm{\&tau;}}_{11}\stackrel{\&RightArrow;}{i}+{\mathrm{\&tau;}}_{12}\stackrel{\&RightArrow;}{j}+{\mathrm{\&tau;}}_{13}\stackrel{\&RightArrow;}{k}---\left(5a\right)$

τ ₁₃=cosα ₁ (5d)

Wherein, α ₁for i.e. the first hole angle of the hole angle of the upper extreme point of micro unit;

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

tangent vector corresponding to the second end points that is micro unit can be expressed as:

${\stackrel{\&RightArrow;}{\mathrm{\&tau;}}}_2={\mathrm{\&tau;}}_{21}\stackrel{\&RightArrow;}{i}+{\mathrm{\&tau;}}_{22}\stackrel{\&RightArrow;}{j}+{\mathrm{\&tau;}}_{23}\stackrel{\&RightArrow;}{k}---\left(6a\right)$

τ ₂₃=cosα ₂ (6d)

Wherein, α ₂for i.e. the second hole angle of the hole angle of the lower extreme point of micro unit;

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:

$\stackrel{\&RightArrow;}{m}=\frac{1}{\sin \mathrm{\&theta;}}{\stackrel{\&RightArrow;}{\mathrm{\&tau;}}}_1\&times;{\stackrel{\&RightArrow;}{\mathrm{\&tau;}}}_2=m_1\stackrel{\&RightArrow;}{i}+m_2\stackrel{\&RightArrow;}{j}+m_3\stackrel{\&RightArrow;}{k}---\left(7\right)$

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:

${\stackrel{\&RightArrow;}{\mathrm{\&tau;}}}_0=\frac{{\stackrel{\&RightArrow;}{\mathrm{\&tau;}}}_1+{\stackrel{\&RightArrow;}{\mathrm{\&tau;}}}_2}{|{\stackrel{\&RightArrow;}{\mathrm{\&tau;}}}_1+{\stackrel{\&RightArrow;}{\mathrm{\&tau;}}}_2|}={\mathrm{\&tau;}}_{01}\stackrel{\&RightArrow;}{i}+{\mathrm{\&tau;}}_{02}\stackrel{\&RightArrow;}{j}+{\mathrm{\&tau;}}_{03}\stackrel{\&RightArrow;}{k}---\left(8\right)$

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:

$\stackrel{\&RightArrow;}{n}=\stackrel{\&RightArrow;}{m}\&times;{\stackrel{\&RightArrow;}{\mathrm{\&tau;}}}_0=n_1\stackrel{\&RightArrow;}{i}+n_2\stackrel{\&RightArrow;}{j}+n_3\stackrel{\&RightArrow;}{k}---\left(9a\right)$

Wherein,

n ₁=m ₂τ ₀₃-m ₃τ ₀₂ (9b)

n ₂=m ₃τ ₀₁-m ₁τ ₀₃ (9c)

n ₃=m ₁τ ₀₂-m ₂τ ₀₁ (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_{\mathrm{ndp}}=-({\mathrm{<figure-callout\; id="1"\; label="T"\; filenames="HDA0000426754530000112.png"\; state="\{\{state\}\}">T</figure-callout>}}_1+T_2)\sin \frac{\mathrm{\&theta;}}{2}+L_s\stackrel{\&RightArrow;}{q}\&CenterDot;\stackrel{\&RightArrow;}{n}---\left(10a\right)$

Also can be expressed as:

$F_{\mathrm{ndp}}=-({\mathrm{<figure-callout\; id="1"\; label="T"\; filenames="HDA0000426754530000112.png"\; state="\{\{state\}\}">T</figure-callout>}}_1+T_2)\sin \frac{\mathrm{\&theta;}}{2}+n_3{L}_{s}q--\left(10b\right)$

Wherein, F _ndpfor the total lateral force in full-shape plane, T ₁for the axial force of described micro unit first end, T ₂for 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,

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

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

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_{\mathrm{np}}=L_s\stackrel{\&RightArrow;}{q}\&CenterDot;\stackrel{\&RightArrow;}{m}=m_{3}q{L}_s---\left(11a\right)$

Wherein, F _npfor the total lateral force in binormal direction, L _sfor the length of described micro unit,

for the effective gravity vector of described micro unit,

for unit binormal vector.

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=\frac{\sqrt{F_{\mathrm{ndp}}^2+F_{\mathrm{np}}^2}}{L_s}---\left(12\right)$

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 α ₁hole angle for the upper extreme point of track measuring point corresponding to described micro unit, is called the first hole angle, α ₂hole 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,

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 ³ (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:

$\stackrel{\&RightArrow;}{q}=q\stackrel{\&RightArrow;}{k}---\left(3\right)$

If the axial force of the micro unit lower end i.e. axial force of the second end is T ₂, the lateral force F of unit length _n, the axial force of its upper end is the axial force T of first end ₁can be calculated by following formula:

${\mathrm{<figure-callout\; id="1"\; label="T"\; filenames="HDA0000426754530000112.png"\; state="\{\{state\}\}">T</figure-callout>}}_1=T_2+\frac{L_s}{\cos \frac{\mathrm{\&theta;}}{2}}[q\cos \stackrel{\&OverBar;}{\mathrm{\&alpha;}}+\mathrm{\&mu;}(F_E+F_n)]---\left(4\right)$

$\stackrel{\&OverBar;}{\mathrm{\&alpha;}}=({\mathrm{\&alpha;}}_1+{\mathrm{\&alpha;}}_2)/2$

Wherein, T ₁for the axial force of described micro unit first end, T ₂for 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, α ₁for the first hole angle of track measuring point corresponding to described micro unit, α ₂for 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:

${\stackrel{\&RightArrow;}{\mathrm{\&tau;}}}_1={\mathrm{\&tau;}}_{11}\stackrel{\&RightArrow;}{i}+{\mathrm{\&tau;}}_{12}\stackrel{\&RightArrow;}{j}+{\mathrm{\&tau;}}_{13}\stackrel{\&RightArrow;}{k}---\left(5a\right)$

τ ₁₃=cosα ₁ (5d)

Wherein, α ₁for i.e. the first hole angle of the hole angle of the upper extreme point of micro unit;

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

tangent vector corresponding to the second end points that is micro unit can be expressed as:

${\stackrel{\&RightArrow;}{\mathrm{\&tau;}}}_2={\mathrm{\&tau;}}_{21}\stackrel{\&RightArrow;}{i}+{\mathrm{\&tau;}}_{22}\stackrel{\&RightArrow;}{j}+{\mathrm{\&tau;}}_{23}\stackrel{\&RightArrow;}{k}---\left(6a\right)$

τ ₂₃=cosα ₂ (6d)

Wherein, α ₂for i.e. the second hole angle of the hole angle of the lower extreme point of micro unit;

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:

$\stackrel{\&RightArrow;}{m}=\frac{1}{\sin \mathrm{\&theta;}}{\stackrel{\&RightArrow;}{\mathrm{\&tau;}}}_1\&times;{\stackrel{\&RightArrow;}{\mathrm{\&tau;}}}_2=m_1\stackrel{\&RightArrow;}{i}+m_2\stackrel{\&RightArrow;}{j}+m_3\stackrel{\&RightArrow;}{k}---\left(7\right)$

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:

${\stackrel{\&RightArrow;}{\mathrm{\&tau;}}}_0=\frac{{\stackrel{\&RightArrow;}{\mathrm{\&tau;}}}_1+{\stackrel{\&RightArrow;}{\mathrm{\&tau;}}}_2}{|{\stackrel{\&RightArrow;}{\mathrm{\&tau;}}}_1+{\stackrel{\&RightArrow;}{\mathrm{\&tau;}}}_2|}={\mathrm{\&tau;}}_{01}\stackrel{\&RightArrow;}{i}+{\mathrm{\&tau;}}_{02}\stackrel{\&RightArrow;}{j}+{\mathrm{\&tau;}}_{03}\stackrel{\&RightArrow;}{k}---\left(8\right)$

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:

$\stackrel{\&RightArrow;}{n}=\stackrel{\&RightArrow;}{m}\&times;{\stackrel{\&RightArrow;}{\mathrm{\&tau;}}}_0=n_1\stackrel{\&RightArrow;}{i}+n_2\stackrel{\&RightArrow;}{j}+n_3\stackrel{\&RightArrow;}{k}---\left(9a\right)$

Wherein,

n ₁=m ₂τ ₀₃-m ₃τ ₀₂ (9b)

n ₂=m ₃τ ₀₁-m ₁τ ₀₃ (9c)

n ₃=m ₁τ ₀₂-m ₂τ ₀₁ (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_{\mathrm{ndp}}=-({\mathrm{<figure-callout\; id="1"\; label="T"\; filenames="HDA0000426754530000112.png"\; state="\{\{state\}\}">T</figure-callout>}}_1+T_2)\sin \frac{\mathrm{\&theta;}}{2}+L_s\stackrel{\&RightArrow;}{q}\&CenterDot;\stackrel{\&RightArrow;}{n}---\left(10a\right)$

Also can be expressed as:

$F_{\mathrm{ndp}}=-({\mathrm{<figure-callout\; id="1"\; label="T"\; filenames="HDA0000426754530000112.png"\; state="\{\{state\}\}">T</figure-callout>}}_1+T_2)\sin \frac{\mathrm{\&theta;}}{2}+n_3{L}_{s}q--\left(10b\right)$

Wherein, F _ndpfor the total lateral force in full-shape plane, T ₁for the axial force of described micro unit first end, T ₂for 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,

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

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

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_{\mathrm{np}}=L_s\stackrel{\&RightArrow;}{q}\&CenterDot;\stackrel{\&RightArrow;}{m}=m_{3}q{L}_s---\left(11a\right)$

Wherein, F _npfor the total lateral force in binormal direction, L _sfor the length of described micro unit,

for the effective gravity vector of described micro unit,

for unit binormal vector.

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=\frac{\sqrt{F_{\mathrm{ndp}}^2+F_{\mathrm{np}}^2}}{L_s}---\left(12\right)$

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+\frac{L_s}{\cos \frac{\mathrm{\&theta;}}{2}}[q\cos \stackrel{\&OverBar;}{\mathrm{\&alpha;}}+\mathrm{\&mu;}(F_E+F_n)]$

Wherein, T ₁for the axial force of described micro unit first end, T ₂for 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,

α ₁for the first hole angle of track measuring point corresponding to described micro unit, α ₂for 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_{\mathrm{ndp}}=-(T_1+T_2)\sin \frac{\mathrm{\&theta;}}{2}+L_s\stackrel{\&RightArrow;}{q}\&CenterDot;\stackrel{\&RightArrow;}{n}$

Wherein, F _ndpfor the total lateral force in full-shape plane, T ₁for the axial force of described micro unit first end, T ₂for 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,

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_{\mathrm{np}}=L_s\stackrel{\&RightArrow;}{q}\&CenterDot;\stackrel{\&RightArrow;}{m}$

Wherein, F _npfor the total lateral force in binormal direction, L _sfor the length of described micro unit,

for the effective gravity vector of described micro unit,

for unit binormal vector.

9. method according to claim 8, is characterized in that, the 4th described relational expression is:

$F_n=\frac{\sqrt{F_{\mathrm{ndp}}^2+F_{\mathrm{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+\frac{L_s}{\cos \frac{\mathrm{\&theta;}}{2}}[q\cos \stackrel{\&OverBar;}{\mathrm{\&alpha;}}+\mathrm{\&mu;}(F_E+F_n)]$

Wherein, T ₁for the axial force of described micro unit first end, T ₂for 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,

α ₁for the first hole angle of track measuring point corresponding to described micro unit, α ₂for 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_{\mathrm{ndp}}=-(T_1+T_2)\sin \frac{\mathrm{\&theta;}}{2}+L_s\stackrel{\&RightArrow;}{q}\&CenterDot;\stackrel{\&RightArrow;}{n}$

Wherein, F _ndpfor the total lateral force in full-shape plane, T ₁for the axial force of described micro unit first end, T ₂for 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,

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_{\mathrm{np}}=L_s\stackrel{\&RightArrow;}{q}\&CenterDot;\stackrel{\&RightArrow;}{m}$

Wherein, F _npfor the total lateral force in binormal direction, L _sfor the length of described micro unit,

for the effective gravity vector of described micro unit,

for unit binormal vector.

18. equipment according to claim 17, is characterized in that, the 4th described relational expression is:

$F_n=\frac{\sqrt{F_{\mathrm{ndp}}^2+F_{\mathrm{np}}^2}}{L_s}.$

Images