Pipe column erosive wear failure position prediction method and device, electronic equipment and storage medium

Abstract

The invention discloses a tubular column erosive wear failure position prediction method and device, electronic equipment and a storage medium. Comprising the steps that a buckling differential equation and a contact force expression of a target tubular column are obtained, and buckling characteristics of the target tubular column under different loads are determined based on the buckling differential equation and the contact force expression; based on the buckling characteristics of a target tubular column under different loads, establishing a buckling flow channel model of the target tubular column; according to the buckling flow channel model of the target tubular column and the material attribute and fluid attribute of the target tubular column, determining the erosion speed distribution of the target tubular column at the target yield and the target sand production rate; according to the erosion speed distribution, the erosion wear failure position of the target tubular column is predicted, the problems that the erosion research cost is high and consumed time is long under the buckling state of the tubular column can be solved, and technical support is provided for safe design of the tubular column, the service life of the tubular column and quantitative calculation of the residual strength.

Description

Technical Field

[0001] This invention relates to the field of oil production engineering technology, and in particular to a method, device, electronic equipment, and storage medium for predicting the location of tubing erosion and wear failure. Background Technology

[0002] With the continuous development of oil and gas resources, major oilfields in my country are exploring and developing deeper and ultra-deep formations. The lower gas-bearing formations on the southern margin of the Junggar Basin are characterized by deep burial (well depths exceeding 8000m), high formation pressure (maximum formation pressure 171.78MPa), and high temperature (maximum formation temperature 170.13℃). These characteristics inevitably lead to various types of damage to oil and gas tubing strings during downhole operation, threatening the safety of oil and gas well production and operations. Among the more dangerous factors are tubing buckling and erosion. In the high-temperature, high-pressure wells of the southern margin, large amounts of high-temperature fluid flow at high speed from the bottom of the well to the wellhead. During this flow, downhole pressure and temperature surge, making the tubing strings prone to severe buckling. Simultaneously, as the downhole fluid discharge increases, the sand-carrying fluid's prolonged contact and collision with the tubing string causes erosion, thinning of the wall, weakening of strength, and reduction of service life. In severe cases, this can lead to tubing deformation, crushing, rupture, and leakage, causing safety problems. In addition, buckling of the downhole tubing increases the incident angle of high-speed particles carried by the fluid, making the buckled area more susceptible to erosion and wear from particles.

[0003] Currently, comprehensive analyses have been conducted on tubular erosion under various factors. However, for research on the combination of buckling and erosion, existing studies primarily employ experimental methods to test erosion under the buckling morphology of the tubular string, which suffers from high costs and time consumption. Summary of the Invention

[0004] To address the aforementioned problems, the inventors have developed this invention, which, through specific embodiments, provides a method, apparatus, electronic device, and storage medium for predicting the location of tubular erosion and wear failure.

[0005] In a first aspect, embodiments of the present invention provide a method for predicting the location of erosion and wear failure in a tubular string, comprising:

[0006] Obtain the buckling differential equation and contact force expression of the target tubing, and determine the buckling characteristics of the target tubing under different loads based on the buckling differential equation and contact force expression.

[0007] Based on the buckling characteristics of the target tubing under different loads, a buckling flow channel model of the target tubing is established;

[0008] Based on the buckling flow channel model of the target tubing, the material properties of the target tubing, and the fluid properties, determine the erosion velocity distribution of the target tubing at the target production rate and the target sand output rate.

[0009] Based on the erosion velocity distribution, the location of erosion wear failure of the target tubing is predicted.

[0010] Secondly, embodiments of the present invention provide a quantitative safety assessment device for buckling and erosion of high-temperature and high-pressure well tubing, comprising:

[0011] The parameter calculation module is used to obtain the buckling differential equation and contact force expression of the target tubing, and based on the buckling differential equation and contact force expression, determine the buckling characteristics of the target tubing under different loads;

[0012] The model building module is used to build a buckling flow channel model of the target tubing based on the buckling characteristics of the target tubing under different loads.

[0013] The analysis module is used to determine the erosion velocity distribution of the target tubing at the target production rate and target sand output rate based on the buckling flow channel model of the target tubing, the material properties of the target tubing, and the fluid properties of the target tubing.

[0014] The result prediction module is used to predict the location of erosion and wear failure of the target tubing based on the erosion velocity distribution.

[0015] The beneficial effects of the above-described technical solutions provided in the embodiments of the present invention include at least the following:

[0016] By obtaining the buckling differential equation and contact force expression of the target tubing, the buckling characteristics of the target tubing under different loads are determined based on the buckling differential equation and contact force expression. Based on the buckling characteristics of the target tubing under different loads, a buckling flow channel model of the target tubing is established. According to the buckling flow channel model of the target tubing, the material properties and fluid properties of the target tubing, the erosion velocity distribution of the target tubing at the target production rate and target sand output rate is determined. According to the erosion velocity distribution, the erosion wear failure location of the target tubing is predicted. This can solve the problems of high cost and long time consumption in erosion research under the buckling mode of the tubing, and provide technical support for the safe design of tubing, the service life of tubing, and the quantitative calculation of residual strength.

[0017] Other features and advantages of the invention will be set forth in the following description, or may be learned by practicing the invention. The objects and other advantages of the invention may be realized and obtained by means of the structures particularly pointed out in the written description, claims, and drawings.

[0018] The technical solution of the present invention will be further described in detail below with reference to the accompanying drawings and embodiments. Attached Figure Description

[0019] The accompanying drawings are provided to further illustrate the invention and form part of the specification. They are used in conjunction with embodiments of the invention to explain the invention and do not constitute a limitation thereof. In the drawings:

[0020] Figure 1 This is a flowchart of the method in an embodiment of the present invention;

[0021] Figure 2 This is a diagram showing the stress and contact geometry of the tubular column in an embodiment of the present invention;

[0022] Figure 3 This is a diagram showing the geometric relationship between the tubing string and the wellbore along the transverse direction in an embodiment of the present invention.

[0023] Figure 4 This is a buckling flow channel model in an embodiment of the present invention;

[0024] Figure 5 These are velocity field distribution cloud maps within the buckling flow channel model under different working conditions in embodiments of the present invention;

[0025] Figure 6 This illustrates the fluid-gravel motion states at different times in region g of the buckling channel model in this embodiment of the invention.

[0026] Figure 7 This is a cloud map showing the erosion velocity distribution on the inner wall of the buckling flow channel model in this embodiment of the invention.

[0027] Figure 8 This refers to the safety window for erosion rate under different sand output and production rates in the embodiments of the present invention;

[0028] Figure 9 This represents a safety window for the erosion wall thickness and sinking speed under different sand output and production rates in embodiments of the present invention.

[0029] Figure 10 This is a schematic diagram of the structure of an electronic device according to an embodiment of the present invention. Detailed Implementation

[0030] Exemplary embodiments of the present disclosure will now be described in more detail with reference to the accompanying drawings. While exemplary embodiments of the present disclosure are shown in the drawings, it should be understood that the present disclosure may be implemented in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the disclosure to those skilled in the art.

[0031] To address the problems existing in the prior art, embodiments of the present invention provide a method, apparatus, electronic device, and storage medium for predicting the location of tubular erosion and wear failure.

[0032] This invention provides a method for predicting the location of erosion and wear failure in tubing, the process of which is as follows: Figure 1 As shown, it includes the following steps:

[0033] Step S1: Obtain the buckling differential equation and contact force expression of the target tubing, and determine the buckling characteristics of the target tubing under different loads based on the buckling differential equation and contact force expression.

[0034] The buckling characteristics can include the length of the helical buckling segment, the pitch, and the flow channel information. The buckling differential equation is used to describe the load borne by the target tubing.

[0035] Step S2: Based on the buckling characteristics of the target tubing under different loads, establish a buckling flow channel model for the target tubing.

[0036] Step S3: Based on the buckling flow channel model of the target tubing, the material properties of the target tubing, and the fluid properties, determine the erosion velocity distribution of the target tubing at the target production rate and the target sand output rate.

[0037] The target output includes one or more output values, and the target sand output includes one or more sand output values.

[0038] Step S4: Based on the erosion velocity distribution, predict the location of erosion wear failure of the target tubing.

[0039] The method described in this embodiment obtains the buckling differential equation and contact force expression of the target tubing. Based on these equations, it determines the buckling characteristics of the target tubing under different loads. Based on these buckling characteristics, it establishes a buckling flow channel model for the target tubing. According to the buckling flow channel model, the material properties, and fluid properties of the target tubing, it determines the erosion velocity distribution at target production rates and target sand output. Based on this erosion velocity distribution, it predicts the erosion wear failure location of the target tubing. This method solves the problems of high cost and long time consumption in erosion research under the buckling morphology of tubing, providing technical support for tubing safety design, service life, and quantitative calculation of remaining strength.

[0040] In the above method of this embodiment, obtaining the buckling differential equation and contact force expression of the target tubing includes the following steps:

[0041] Obtain the external force vectors and force equilibrium differential equations of the target tubing;

[0042] Based on the external force vector acting on the target tubing, determine the projection of the torque acting on the target tubing in the Oxyz coordinate system;

[0043] Based on the force equilibrium differential equation and the projection of the torque on the target pipe in the Oxyz coordinate system, the deformation control equation of the target pipe is determined.

[0044] Based on the deformation control equation of the target tubing, the buckling differential equation and contact force expression of the target tubing are determined.

[0045] For example, assume the target tubing string is constrained by a straight wellbore. At both ends of the tubing string, it is subjected to axial compressive loads F and torques MT due to the wellhead hangers and downhole packers, as well as weight and wellbore contact forces N. Assume that before buckling, the target tubing string is straight and tightly adhered to the lower wellbore wall. A schematic diagram in the Oxyz coordinate system is shown below. Figure 2 As shown, within the wellbore cross-section, the geometric relationship between the target tubing and the wellbore wall is as follows: Figure 3 As shown.

[0046] The target tubing is subjected to the following external force vector:

[0047] h=(q sinα-N cosθ)iN sinθj-q cosαk,

[0048] In the formula, α is the well inclination angle, °; θ is the angle between the vector and the x direction, °; ​​N is the contact force between the tubing and the wellbore, N; i, j, and k are the direction vectors in the x, y, and z directions in the Oxyz coordinate system, respectively; and q is the weight of the tubing string, N.

[0049] The force equilibrium differential equation for the target tubing is:

[0050]

[0051] The projection of the moment acting on the target tubing in the Oxyz coordinate system is as follows:

[0052]

[0053] In the formula, E is the elastic modulus of the tubular material, Pa; I is the moment of inertia of the tubular cross-section about the neutral axis, m. 4 G is the shear modulus of elasticity of the tubular material, Pa; J is the polar moment of inertia of the cross section, m. 4 γ is the torsion angle, in °.

[0054] Based on the deformation control equations of the bottom drill string assembly derived by Gao Deli et al., the expressions for the complex functions ψ and H are given as follows:

[0055] ψ = r c e iθ =r c (cosθ+isinθ)

[0056] H=(q * sinα-N * cosθ)+N * isinθ.

[0057] In the formula, q * =q / EI,N * =N / EI,

[0058] Therefore, the deformation control equation for the target tubing is:

[0059]

[0060] In the formula,

[0061] Regarding e iθ The derivatives of each order are used to derive and substitute into the deformation control equation of the target tubular column. The real and imaginary parts are then simplified and analyzed to obtain the buckling differential equation and the contact force expression. The buckling differential equation and the contact force expression are as follows:

[0062]

[0063] In the method described in this embodiment, obtaining the buckling differential equation and contact force expression of the target tubing, and determining the buckling characteristics of the target tubing under different loads based on the buckling differential equation and contact force expression, includes the following steps:

[0064] Obtain the buckling differential equation, contact force expression, and geometric model of the target tubing;

[0065] Based on the buckling differential equation and contact force expression of the target tubing, the axial compressive load of the target tubing geometric model is adjusted to determine the buckling characteristics of the target tubing.

[0066] The target tubing geometric model is established based on the material properties, structure, length, and packer setting depth of the target tubing. Based on finite element analysis technology and the target tubing geometric model, the buckling characteristics of the target tubing under different loads are determined according to the production, collapse pressure, internal pressure, and temperature of the target tubing.

[0067] Specifically, the loads borne by the tubing string have been described in detail in the buckling differential equation, mainly including: 1. the lifting force FH at the wellhead, applied by the wellhead equipment; 2. gravity GT, the weight of the tubing string itself; 3. the bottom hole compressive force FB, applied by the packer downhole; 4. the internal pressure and collapse pressure caused by the fluid. The internal pressure is determined by the bottom hole pressure and fluid density, while the external pressure is determined by the annular fluid and the wellhead pressurization value. In the finite element analysis, the bottom axial compressive load generated by the packer is adjusted, and the preset load is gradually increased from 0 kN to determine the buckling characteristics of the target tubing string under different loads.

[0068] In the above method of this embodiment, establishing a buckling flow channel model of the target tubing based on the buckling characteristics of the target tubing under different loads may include the following steps:

[0069] The pipe subjected to complex stress is simplified into a helical cylinder. Based on parameters such as the radius, number of turns, and height of the inner and outer rings of the target pipe, a buckling flow channel model of the target pipe is established.

[0070] Set the fluid properties according to the actual working conditions, including viscosity and density;

[0071] The physical fields of the buckling flow channel model are established, considering interphase coupling and momentum coupling, i.e., the action and reaction forces between phases, following Newton's third law. Simultaneously, each particle in the solid phase is solved individually to obtain detailed particle motion dynamics information, including particle-fluid and particle-wall interactions.

[0072] In the above method of this embodiment, determining the erosion velocity distribution of the target tubing at the target production rate and target sand output rate based on the buckling flow channel model of the target tubing, the material properties of the target tubing, and the fluid properties of the target tubing includes the following steps:

[0073] Based on the buckling flow channel model of the target tubing, the material properties and fluid properties of the target tubing, the motion states of the fluid and sand particles at different times are obtained; based on the motion states of the fluid and sand particles at different times, the erosion velocity distribution of the target tubing at the target output and target sand discharge is determined.

[0074] In the above method of this embodiment, determining the erosion velocity distribution of the target tubing at the target production rate and target sand output rate based on the buckling flow channel model of the target tubing, the material properties of the target tubing, and the fluid properties of the target tubing includes the following steps:

[0075] Using the RNG k-ε turbulence model, based on the buckling channel model of the target tubing, the material properties of the target tubing, and the fluid properties, the erosion velocity distribution of the target tubing at the target production rate and the target sand output rate is determined.

[0076] Among them, the RNG k-ε turbulence model is used to calculate turbulent flow in computational fluid dynamics (CFD). It has a fast convergence speed during calculation and high accuracy in the flow field calculation when there are phenomena such as strong curvature flow, separated flow, and adverse pressure gradient. The RNG k-ε model is a k-epsilon turbulence model based on the renormalization group.

[0077] The transport equations for the RNG k-ε turbulence model are:

[0078]

[0079]

[0080] Among them, G k It is the turbulent kinetic energy generated by the laminar velocity gradient, G b Y is the turbulent kinetic energy generated by buoyancy. M It is a wave generated by diffusion during the transition in compressible turbulence, α k and α ε It is the reciprocal of the effective Prandtl number for turbulence in the k-equation and the ε-equation, S k and S ε It is a user-defined source item, C 1ε C 2ε C 3ε It is a constant, C is obtained through analysis and derivation in the RNG k-ε turbulence model. 1ε =1.42, C 2ε =1.68.

[0081] The scale process is eliminated in RNG by the following turbulent viscosity equation:

[0082]

[0083] C ν ≈100,

[0084] Under the high Reynolds number constraint, we obtain

[0085]

[0086] In the formula, C is derived through RNG theory. μ =0.0845, which is very close to the empirical estimate of 0.09 in the standard k-ε model.

[0087] (1) Turbulent flow is affected by vortices in laminar flow. In CFD, these effects are corrected by modifying the turbulent viscosity. The correction is as follows:

[0088]

[0089] In the formula, μ t0 This is the turbulent viscosity value. Ω is a characteristic vortex number estimated in CFD considering vortices, and α... s It is a constant, depending on whether the fluid flow is a predominantly vortex or moderately vortex flow. For moderately vortex flow, α s =0.07.

[0090] The reciprocal of the Prandtl number in the transmission equation is calculated using the following formula:

[0091]

[0092] When α0 = 1.0, under the high Reynolds number condition, α k and α ε ≈1.393.

[0093] The main difference between the RNG k-ε turbulence model and the standard k-ε model lies in the additional terms in the ε equation:

[0094]

[0095] Among them, eta=Sk / ε, eta0=4.38, and β=0.012.

[0096] In the k-ε turbulence model, the effect of buoyancy on turbulence needs to be considered. The effect of buoyancy is given by the following equation:

[0097]

[0098] In the formula, Pr t It is the Prandtl number of turbulence energy, g i It is the component of gravity in the i-th direction. In the RNG model, Pr t = 1 / α, where β is the coefficient of thermal expansion, defined as:

[0099]

[0100] The degree to which the ε equation is affected by buoyancy depends on the constant C. 3ε It can be calculated using the following formula:

[0101]

[0102] In the formula, v is the velocity component of the fluid parallel to gravity, and u is the velocity component perpendicular to gravity.

[0103] According to fluid mechanics, turbulence near the wall is divided into a viscous sublayer, a buffer zone, and a logarithmic law region. Since the velocity gradient of the boundary layer is large, the boundary layer needs to be set reasonably. That is, the size of the first row of boundary elements should be small enough, and then the thickness adjustment factor and boundary layer stretching factor suitable for the model should be determined.

[0104] Numerical simulations based on the RNG k-ε turbulence model were performed to obtain the motion states of the fluid and sand particles at different times, namely the fluid distribution inside the target pipe and the changes of the sand and gravel group over time. This allowed for the determination of the power source carrying the sand and gravel movement and the concentrated area of ​​collision and erosion between the sand and gravel and the inner wall.

[0105] The method described in this embodiment further includes the following steps:

[0106] The erosion rate at the erosion wear failure location of the target tubing is calculated sequentially at a preset sand output and a preset production rate, and an erosion rate safety window is generated.

[0107] The method described in this embodiment further includes the following steps:

[0108] The erosion wall thickness settling velocity at the erosion wear failure location is calculated sequentially when the target tubing is subjected to a preset sand output and a preset production rate, thereby generating a safe window for the erosion wall thickness settling velocity.

[0109] In the method described in this embodiment, a safety window for erosion velocity and a safety window for erosion wall thickness subsidence velocity are generated, providing a guidance tool for erosion prevention and control under different on-site operating and production conditions.

[0110] Optionally, erosion simulations were performed on tubular strings with different production rates and sand outputs using the same method to obtain variation curves. Then, the Finnie model was used to calculate the location of the most severe buckling section under specific operating conditions. Compared to the Oka, E / CRC, and DNV models, the Finnie model comprehensively considers the inner wall material properties, gravel size, and particle shape, while omitting empirical parameters. The Finnie erosion model measures the incident particles per unit mass (in mm). 3 The volume of surface material removed (in units of / kg) is defined as:

[0111] E(α)=g(α)E 90 ,

[0112]

[0113] Where E(α) is the etching rate, kg / m 2 .s -1 E 90 For vertical incident etching rate, kg / m 2 .s -1 Hv is the hardness of the inner wall of the tubing; v and v' are the velocity and reference velocity, m / s; D and D' are the outer diameter and reference diameter of the tubing, m; k1 and k3 are constants.

[0114]

[0115] Where q1, q2, s1, and s2 are constants. The erosion rate of the oil pipe per second, ΔV. es and the annual wall thickness reduction rate Δd ey It can be represented as:

[0116]

[0117] Where A eff The area of ​​the eroded region is m. 2 ;ve Erosion rate and etching rate, kg / m 2 .s -1 ;ρ t The density of the pipe is kg / m 3 .

[0118] For example, the buckling characteristics of the target tubing under different loads were determined using Abaqus / CAE software. The material properties of the tubing are shown in Table 1. The target tubing depth is 7592m, the packer is set at 7290m, the wellhead oil pressure is 90MPa, the wellhead temperature is 135℃, the annulus protection fluid density is 1.3g / cm3, and the gravel diameter is 1mm with a density of 2750kg / m³. 3 Sand output is 1000 kg / day. Operating condition 1: 620 cubic meters of oil / day - 200,000 cubic meters of gas / day. Operating condition 2: 800 cubic meters of oil / day - 270,000 cubic meters of gas / day.

[0119]

[0120]

[0121] Table 1

[0122] After establishing the geometric model of the target tubing, the axial compressive load of the target tubing geometric model is adjusted based on the buckling differential equation and contact force expression of the target tubing to determine the buckling characteristics of the target tubing under operating conditions 1 and 2. The external pressure can be directly calculated from the density of the annular fluid and the pressurization conditions, while the internal pressure is directly calculated from the bottom hole pressure and the density of the fluid in the tubing.

[0123] A 3D flow channel is created in CAD software, and after meshing, a buckling flow channel model is obtained. The buckling flow channel model is as follows: Figure 4 As shown, the velocity field distribution contour maps of the buckling flow channel model under operating conditions 1 and 2 are as follows: Figure 5 As shown in the figure, the fluid velocity varies depending on the production rate. In the buckling channel model, the fluid velocity distribution is lower near the wall due to the fluid's viscosity. However, due to the complex variations in the buckling channel, the maximum velocity distribution differs across different cross-sections. The velocity distribution at the four different cross-sections clearly deviates from the center.

[0124] Figure 6 The motion state of fluid-gravel in region g of the helical buckling flow channel within 0.02-0.5s is shown. Initially, the gravel moves in the pipe according to the trend of the fluid velocity field. The velocity of the gravel tends to increase as it gets closer to the center of the pipe. When the gravel moves to the first bend, due to the change of the flow channel and the inertia of the gravel itself, most of the gravel collides with the apex of the helical section.

[0125] like Figure 7 As shown, erosion is most severe at the apex of the spiral section under all operating conditions. The maximum erosion velocity in this region under different operating conditions is 8.21 × 10⁻⁶. -8 kg / m 2 .s -1 With 1.48×10 -7 kg / m 2 .s -1 By calculating and analyzing the erosion of buckled tubular columns under different production rates, it was found that the erosion rate increases significantly with increasing production rate, and the severely eroded areas are distributed on the periphery of the helical buckled tubular columns.

[0126] Based on the method described in the above embodiments, erosion simulations were performed on targets at different production rates. As oil production increased, the erosion rate increased significantly. As the sand output Qs increased, the erosion rate of the tubing increased linearly.

[0127] The location of the most severe buckling section was selected for simulation calculation using the erosion model to obtain the erosion wall thickness subsidence rate. Under working conditions 1 and 2, the erosion wall thickness subsidence rate at the location of the most severe buckling section reached 0.16 and 0.29 mm / year, respectively.

[0128] Figure 8 The inner wall erosion rate is a safe window for different sand output and production rates. The inner wall erosion rate changes non-linearly with the sand output and oil production. Figure 9 To establish a safe window for the erosion wall thickness settling velocity under different sand output and production rates, the erosion wall thickness settling velocity exhibits a non-linear relationship with changes in sand output and oil production. In the calculated operating conditions, the wall thickness settling velocity is less than 0.3 mm / year, while when the oil production and sand output reach a certain level, the erosion wall thickness settling velocity can exceed 1 mm / year.

[0129] Those skilled in the art can change the above order without departing from the scope of protection of this disclosure.

[0130] Another embodiment of the present invention provides a quantitative safety assessment device for buckling and erosion of high-temperature and high-pressure well tubing, comprising:

[0131] The parameter calculation module is used to obtain the buckling differential equation and contact force expression of the target tubing, and based on the buckling differential equation and contact force expression, determine the buckling characteristics of the target tubing under different loads;

[0132] The model building module is used to build a buckling flow channel model of the target tubing based on the buckling characteristics of the target tubing under different loads.

[0133] The analysis module is used to determine the erosion velocity distribution of the target tubing at the target production rate and target sand output rate based on the buckling flow channel model of the target tubing, the material properties of the target tubing, and the fluid properties of the target tubing.

[0134] The result prediction module is used to predict the location of erosion and wear failure of the target tubing based on the erosion velocity distribution.

[0135] Regarding the system in the above embodiments, the specific ways in which each module performs operations have been described in detail in the embodiments related to the method, and will not be elaborated here.

[0136] Based on the same inventive concept, embodiments of the present invention also provide an electronic device, the structure of which is as follows: Figure 10 As shown, it includes: a memory, a processor, and a computer program stored in the memory and running on the processor. When the processor executes the computer program, it implements the aforementioned method for predicting the location of tubular erosion and wear failure.

[0137] Based on the same inventive concept, embodiments of the present invention also provide a computer storage medium storing computer-executable instructions, which, when executed by a processor, implement the aforementioned method for predicting the location of tubular erosion and wear failure.

[0138] Any modifications, additions, and equivalent substitutions made within the scope of the principles of this invention shall still fall within the patent coverage of this invention.

Claims

1. A method for predicting the location of erosion and wear failure in tubing, characterized in that, Includes the following steps: Obtain the buckling differential equation and contact force expression of the target tubing, and determine the buckling characteristics of the target tubing under different loads based on the buckling differential equation and contact force expression. Based on the buckling characteristics of the target tubing under different loads, a buckling flow channel model of the target tubing is established; Based on the buckling flow channel model of the target tubing, the material properties of the target tubing, and the fluid properties, determine the erosion velocity distribution of the target tubing at the target production rate and the target sand output rate. Based on the erosion velocity distribution, the location of erosion wear failure of the target tubing is predicted.

2. The method as described in claim 1, characterized in that, The process of obtaining the buckling differential equation and contact force expression of the target tubing, and determining the buckling characteristics of the target tubing under different loads based on the buckling differential equation and contact force expression, includes the following steps: Obtain the buckling differential equation, contact force expression, and geometric model of the target tubing; Based on the buckling differential equation and contact force expression of the target tubing, the axial compressive load of the target tubing geometric model is adjusted to determine the buckling characteristics of the target tubing.

3. The method as described in claim 2, characterized in that, The determination of the erosion velocity distribution of the target tubing at target production rate and target sand output rate, based on the buckling flow channel model, material properties, and fluid properties of the target tubing, includes the following steps: Based on the buckling flow channel model of the target tubing, the material properties and fluid properties of the target tubing, the motion states of the fluid and sand particles at different times are obtained; Based on the motion states of the fluid and sand particles at different times, the erosion velocity distribution of the target tubing at the target production rate and target sand output is determined.

4. The method as described in claim 3, characterized in that, It also includes the following steps: The erosion rate at the erosion wear failure location of the target tubing is calculated sequentially at a preset sand output and a preset production rate, and an erosion rate safety window is generated.

5. The method as described in claim 4, characterized in that, It also includes the following steps: The erosion wall thickness settling velocity at the erosion wear failure location is calculated sequentially when the target tubing is subjected to a preset sand output and a preset production rate, thereby generating a safe window for the erosion wall thickness settling velocity.

6. The method as described in claim 5, characterized in that, The process of obtaining the buckling differential equation and contact force expression of the target tubing includes the following steps: Obtain the external force vectors and force equilibrium differential equations of the target tubing; Based on the external force vector acting on the target tubing, determine the projection of the torque acting on the target tubing in the Oxyz coordinate system; Based on the force equilibrium differential equation and the projection of the torque on the target pipe in the Oxyz coordinate system, the deformation control equation of the target pipe is determined. Based on the deformation control equation of the target tubing, the buckling differential equation and contact force expression of the target tubing are determined.

7. The method as described in claim 6, characterized in that, The determination of the erosion velocity distribution of the target tubing at target production rate and target sand output rate, based on the buckling flow channel model, material properties, and fluid properties of the target tubing, includes the following steps: Using the RNG k-ε turbulence model, based on the buckling channel model of the target tubing, the material properties of the target tubing, and the fluid properties, the erosion velocity distribution of the target tubing at the target production rate and the target sand output rate is determined.

8. A quantitative safety evaluation device for buckling and erosion of high-temperature and high-pressure well tubing strings, characterized in that, include: The parameter calculation module is used to obtain the buckling differential equation and contact force expression of the target tubing, and based on the buckling differential equation and contact force expression, determine the buckling characteristics of the target tubing under different loads; The model building module is used to build a buckling flow channel model of the target tubing based on the buckling characteristics of the target tubing under different loads. The analysis module is used to determine the erosion velocity distribution of the target tubing at the target production rate and target sand output rate based on the buckling flow channel model of the target tubing, the material properties of the target tubing, and the fluid properties of the target tubing. The result prediction module is used to predict the location of erosion and wear failure of the target tubing based on the erosion velocity distribution.

9. An electronic device, characterized in that, include: The memory, the processor, and the computer program stored in the memory and running on the processor, wherein the processor, when executing the computer program, implements the method for predicting the location of tubular erosion and wear failure as described in any one of claims 1 to 7.

10. A computer storage medium, characterized in that, The computer storage medium stores computer-executable instructions, which, when executed, implement the pipe string erosion and wear failure location prediction method according to any one of claims 1 to 7.

Images