Method and device for determining pressure of drill tube in high-temperature and high-pressure well

Abstract

The invention discloses a method and a device for determining the pressure of a high-temperature and high-pressure well tripping borehole, and relates to the technical field of petroleum and natural gas drilling, and the method comprises the following steps: acquiring parameter data in a tripping operation process; grid nodes are divided according to the preset node number and the preset point well depth, and the infinitesimal segment length is obtained; determining hole drift angles at the grid nodes according to the inclinometry data; and according to the mass conservation equation, the momentum conservation equation, the energy conservation equation, the length of the infinitesimal section, the hole drift angle at the grid node and the parameter data, determining the pressure of the tripping well cylinder under temperature and pressure coupling. Therefore, the effect of high temperature and high pressure on the property of the drilling fluid is considered, the pressure calculation precision of the tripping well is improved, the accuracy of shaft pressure prediction is improved, the requirement for accurate regulation and control of the pressure of the tripping well of a high-temperature and high-pressure well is effectively met, and safe and efficient tripping construction operation is guaranteed. The method can improve the pressure calculation precision of the tripping well and improve the prediction accuracy of the pressure of the shaft.

Description

Technical Field

[0001] The present invention relates to the technical field of oil and gas drilling, and in particular to a method and a device for determining the wellbore pressure during tripping of a high-temperature and high-pressure well. Background Art

[0002] This section is intended to provide a background or context to the embodiments of the invention recited in the claims. No admission is made that the description herein is prior art by inclusion in this section.

[0003] The deep formation has prominent characteristics of high temperature and high pressure, and the safe density window of drilling fluid is narrow. With the increasing exploration and development of deep oil and gas resources, the requirements for wellbore pressure prediction during tripping and drilling operations are becoming increasingly higher. Accurately predicting the wellbore pressure during tripping and drilling is of practical significance for ensuring drilling safety and improving operation efficiency.

[0004] The existing wellbore pressure calculation for tripping and drilling does not consider the interaction between wellbore temperature, wellbore pressure and drilling fluid properties. For deep and ultra-deep wells, ignoring the effects of high temperature and high pressure will lead to large deviations in wellbore pressure prediction, causing wellbore pressure imbalance, leading to complex situations such as well leakage and overflow. The existing methods for determining wellbore pressure for tripping and drilling have low accuracy, resulting in greater safety risks in drilling operations. Summary of the invention

[0005] The embodiment of the present invention provides a method for determining the wellbore pressure of a high-temperature and high-pressure well during tripping, which is used to improve the accuracy of wellbore pressure prediction, effectively meet the accurate regulation of the wellbore pressure of a high-temperature and high-pressure well during tripping, and ensure safe and efficient tripping operations. The method includes:

[0006] Obtaining parameter data during the tripping operation; wherein the parameter data at least includes: inclination measurement data and well depth at a preset point;

[0007] According to the preset number of nodes and the preset well depth, the grid nodes are divided and the length of the microelement segment is obtained;

[0008] Determine the well inclination angle at the grid node according to the inclination survey data;

[0009] The tripping wellbore pressure under temperature-pressure coupling is determined based on the mass conservation equation, momentum conservation equation, energy conservation equation, microelement segment length, well inclination angle at grid nodes, and parameter data.

[0010] The embodiment of the present invention further provides a device for determining the wellbore pressure during tripping of a high-temperature and high-pressure well, which is used to improve the accuracy of wellbore pressure prediction, effectively meet the accurate regulation of the wellbore pressure during tripping of a high-temperature and high-pressure well, and ensure safe and efficient tripping operations. The device includes:

[0011] The parameter acquisition module is used to acquire parameter data during the tripping operation; wherein the parameter data at least includes: inclination measurement data and well depth at a preset point;

[0012] The first processing module is used to divide the grid nodes and obtain the micro-element segment length according to the preset node number and the preset point well depth;

[0013] The second processing module is used to determine the well inclination angle at the grid node according to the inclination measurement data;

[0014] The third processing module is used to determine the tripping wellbore pressure under temperature-pressure coupling according to the mass conservation equation, momentum conservation equation, energy conservation equation, microelement segment length, well inclination angle at grid nodes, and parameter data.

[0015] An embodiment of the present invention also provides a computer device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the above-mentioned method for determining the wellbore pressure during drilling of a high-temperature and high-pressure well when executing the computer program.

[0016] An embodiment of the present invention further provides a computer-readable storage medium storing a computer program, which, when executed by a processor, implements the above-mentioned method for determining the wellbore pressure during tripping and drilling of a high-temperature and high-pressure well.

[0017] An embodiment of the present invention further provides a computer program product, which includes a computer program. When the computer program is executed by a processor, the method for determining the wellbore pressure during tripping and drilling of a high-temperature and high-pressure well is implemented.

[0018] In the embodiment of the present invention, parameter data in the process of drilling operation is obtained; wherein the parameter data at least includes: inclination data and well depth of preset points; grid nodes are divided according to the number of preset nodes and the well depth of preset points, and the length of micro-element segments is obtained; the well inclination angle at the grid node is determined according to the inclination data; the wellbore pressure of drilling under temperature-pressure coupling is determined according to the mass conservation equation, momentum conservation equation, and energy conservation equation, micro-element segment length, well inclination angle at the grid node, and parameter data. In this way, the effect of high temperature and high pressure on the properties of drilling fluid is taken into account, the calculation accuracy of wellbore pressure of drilling is improved, the accuracy of wellbore pressure prediction is improved, and the accurate regulation of wellbore pressure of drilling in high temperature and high pressure wells is effectively met, ensuring safe and efficient drilling operation. BRIEF DESCRIPTION OF THE DRAWINGS

[0019] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the following briefly introduces the drawings required for use in the embodiments or the prior art descriptions. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work. In the drawings:

[0020] Figure 1 A flow chart of a method for determining wellbore pressure during tripping of a high-temperature and high-pressure well provided in an embodiment of the present invention;

[0021] Figure 2 A flow chart of a method for dividing grid nodes and obtaining microelement segment length according to a preset number of nodes and a preset well depth provided in an embodiment of the present invention;

[0022] Figure 3 A flow chart of a method for determining the tripping wellbore pressure under temperature-pressure coupling according to a mass conservation equation, a momentum conservation equation, an energy conservation equation, a microelement segment length, a well inclination angle at a grid node, and parameter data provided in an embodiment of the present invention;

[0023] Figure 4 A flow chart of a method for calculating the drilling fluid density, the drilling fluid yield stress, the drilling fluid plastic viscosity, the drilling fluid specific heat capacity and the formation thermal conductivity at a preset pressure value and a preset temperature value provided in an embodiment of the present invention;

[0024] Figure 5 A schematic diagram of a device for determining the wellbore pressure during tripping of a high-temperature and high-pressure well provided in an embodiment of the present invention;

[0025] Figure 6 A schematic diagram of a computer device provided in an embodiment of the present invention. DETAILED DESCRIPTION

[0026] To make the purpose, technical solution and advantages of the embodiments of the present invention more clear, the embodiments of the present invention are further described in detail below in conjunction with the accompanying drawings. Here, the exemplary embodiments of the present invention and their descriptions are used to explain the present invention, but are not intended to limit the present invention.

[0027] The acquisition, storage, use, and processing of data in the technical solution of this application comply with the relevant provisions of national laws and regulations.

[0028] The term "and / or" herein only describes an association relationship, indicating that three relationships may exist. For example, A and / or B may represent the following three situations: A exists alone, A and B exist at the same time, and B exists alone. In addition, the term "at least one" herein represents any combination of at least two of any one or more of a plurality of. For example, including at least one of A, B, and C may represent including any one or more elements selected from the set consisting of A, B, and C.

[0029] In the description of this specification, the terms "include", "including", "have", "contain", etc. are all open terms, which mean including but not limited to. The descriptions with reference to the terms "one embodiment", "a specific embodiment", "some embodiments", "for example", etc. mean that the specific features, structures or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the present application. In this specification, the schematic representations of the above terms do not necessarily refer to the same embodiment or example. Moreover, the specific features, structures or characteristics described may be combined in a suitable manner in any one or more embodiments or examples. The order of steps involved in each embodiment is used to schematically illustrate the implementation of the present application, and the order of steps is not limited and can be appropriately adjusted as needed.

[0030] Research has found that deep formations have prominent high-temperature and high-pressure characteristics and a narrow window of safe density for drilling fluids. With the increasing efforts in the exploration and development of deep oil and gas resources, the requirements for predicting wellbore pressure during tripping operations are becoming increasingly high. Accurately predicting wellbore pressure during tripping operations is of practical significance for ensuring drilling safety and improving operational efficiency.

[0031] The existing wellbore pressure calculation for tripping and drilling does not consider the interaction between wellbore temperature, wellbore pressure and drilling fluid properties. For deep and ultra-deep wells, ignoring the effects of high temperature and high pressure will lead to large deviations in wellbore pressure prediction, causing wellbore pressure imbalance, leading to complex situations such as well leakage and overflow. The existing methods for determining wellbore pressure for tripping and drilling have low accuracy, resulting in greater safety risks in drilling operations.

[0032] Regarding the above studies, Figure 1 As shown, an embodiment of the present invention provides a method for determining the wellbore pressure during tripping of a high-temperature and high-pressure well, comprising:

[0033] S101: Acquire parameter data during the tripping operation; wherein the parameter data at least includes: inclination measurement data and well depth at a preset point;

[0034] S102: Divide the grid nodes according to the preset number of nodes and the preset well depth, and obtain the micro-element segment length;

[0035] S103: determining the well inclination angle at the grid node according to the inclination survey data;

[0036] S104: Determine the tripping wellbore pressure under temperature-pressure coupling according to the mass conservation equation, momentum conservation equation, energy conservation equation, microelement segment length, well inclination angle at grid nodes, and parameter data.

[0037] In the embodiment of the present invention, parameter data in the process of drilling operation is obtained; wherein the parameter data at least includes: inclination data and well depth of preset points; according to the number of preset nodes and the well depth of preset points, the grid nodes are divided, and the length of the micro-element segment is obtained; the well inclination angle at the grid node is determined according to the inclination data; according to the mass conservation equation, the momentum conservation equation, and the energy conservation equation, the micro-element segment length, the well inclination angle at the grid node, and the parameter data, the wellbore pressure of drilling under temperature-pressure coupling is determined. In this way, the effect of high temperature and high pressure on the properties of drilling fluid is taken into account, the calculation accuracy of the wellbore pressure of drilling is improved, the accuracy of wellbore pressure prediction is improved, and the accurate regulation of wellbore pressure of drilling in high temperature and high pressure wells is effectively met, ensuring safe and efficient drilling operation.

[0038] The method for determining the wellbore pressure during tripping and drilling of the above-mentioned high-temperature and high-pressure well is described in detail below.

[0039] Regarding the above S102, if Figure 2 As shown, a flow chart of a method for dividing grid nodes and obtaining microelement segment length according to a preset number of nodes and a preset well depth provided by an embodiment of the present invention includes:

[0040] S201: According to the preset number of nodes and the total well depth, the grid nodes are divided using an equal spacing or a variable spacing method to obtain an initial grid node matrix.

[0041] In one embodiment of the present invention, when equal-interval division is adopted, the average well depth between adjacent preset grid nodes is first calculated according to the preset node quantity and the total well depth, and then the initial grid node matrix is ​​obtained according to the average well depth of each preset grid node.

[0042] For example, if the well depth is 3000m and the number of preset nodes is 5, the average well depth of adjacent preset grid nodes is 600, and the initial grid node matrix is ​​[600, 1200, 1800, 2400, 3000].

[0043] In another embodiment of the present invention, when variable spacing division is adopted, the interval of the variable spacing is set in advance according to the preset number of nodes and the total well depth. For example, the spacing can be set to increase or decrease. Then, the initial grid node matrix is ​​obtained according to the preset variable spacing interval, the preset number of nodes, and the total well depth.

[0044] For example, if the well depth is 3000m, the number of preset nodes is 5, and the spacing between adjacent preset grid nodes is preset to increase, the initial grid node matrix is ​​[600, 1000, 1600, 2400, 3000].

[0045] S202: Add the preset well depth to the initial grid node matrix to obtain a final divided grid node matrix; wherein each element value in the grid node matrix represents the well depth corresponding to a grid node.

[0046] The preset well depths include, for example, the well depth at the casing shoe, the well depth at the drill string diameter change point, and the well depth at the abnormal formation pressure.

[0047] Specifically, according to the well depth of the preset point and the well depth of the preset grid node, the well depth of the preset point is added to the initial grid node matrix to obtain the final divided grid node matrix; wherein each element value in the grid node matrix represents the well depth corresponding to a grid node.

[0048] For example, the initial grid node matrix is ​​[600, 1000, 1600, 2400, 3000], and the preset well depths include: 700m at the casing shoe, 1050m at the drill string diameter change, and 2600m at the abnormal formation pressure. The final divided grid node matrix is ​​[600, 700, 1000, 1050, 1600, 2400, 2600, 3000].

[0049] S203: Determine the length of the microelement segment according to the preset number of nodes and the total well depth.

[0050] In one embodiment of the present invention, the well depth between preset grid nodes is determined according to the total well depth and the number of preset nodes, and the well depth between the preset grid nodes is used as the micro-element segment length.

[0051] With respect to the above S103, determining the well inclination angle at the grid node according to the inclination measurement data, for example, includes: calculating the well inclination angle at the grid node by using a linear interpolation method according to the inclination measurement data.

[0052] In one embodiment of the present invention, the inclination survey data includes, for example: the well depth and the well inclination angle at the inclination survey point above the grid node, and the well depth and the well inclination angle at the inclination survey point below the grid node.

[0053] Specifically, the well inclination angle at the grid node is calculated using the linear interpolation method according to the inclination survey data, including: using the well depth and well inclination angle at the inclination survey point above the grid node, and the well depth and well inclination angle at the inclination survey point below the grid node, to calculate the well inclination angle at the grid node:

[0054]

[0055] Where θ is the well inclination angle at the grid node; D m is the well depth at the grid node; θ1 is the well inclination angle at the inclination measurement point above the grid node; D m,1 is the well depth at the inclination point above the grid node; θ2 is the well inclination angle at the inclination point below the grid node; D m,2 It is the well depth at the inclinometer point below the grid node.

[0056] For the above S104, the parameter data also includes, for example: the total heat transfer coefficient between cement ring, casing and annular fluid, the total heat transfer coefficient between drill pipe and fluid in the pipe, dimensionless temperature, drilling fluid mass flow rate, borehole radius, drill pipe inner radius, test temperature value, test pressure value, annular drilling fluid flow rate, inclination data, drilling fluid density at test temperature value and test pressure value, drilling fluid yield stress at test temperature value and test pressure value, drilling fluid plastic viscosity at test temperature value and test pressure value, drilling fluid specific heat capacity at test temperature value and test pressure value, inlet drilling fluid temperature, formation specific heat capacity, formation thermal conductivity at test temperature value and test pressure value, surface formation temperature, geothermal gradient and total well depth, preset point well depth, inlet drilling fluid temperature, drilling fluid flow rate, annular cross-sectional area. Figure 3 As shown, a flow chart of a method for determining the tripping wellbore pressure under temperature-pressure coupling according to the mass conservation equation, momentum conservation equation, energy conservation equation, microelement segment length, well inclination angle at the grid node, and parameter data provided by an embodiment of the present invention, comprising: for each grid node, executing the following process:

[0057] S301: According to the preset pressure value, preset temperature value and parameter data of the grid node, the drilling fluid density at the preset pressure value and preset temperature value, the drilling fluid yield stress at the preset pressure value and preset temperature value, the drilling fluid plastic viscosity at the preset pressure value and preset temperature value, the drilling fluid specific heat capacity at the preset pressure value and preset temperature value and the formation thermal conductivity at the preset pressure value and preset temperature value are calculated.

[0058] like Figure 4 As shown, the drilling fluid density at a preset pressure value and a preset temperature value, the drilling fluid yield stress at a preset pressure value and a preset temperature value, the drilling fluid plastic viscosity at a preset pressure value and a preset temperature value, the drilling fluid specific heat capacity at a preset pressure value and a preset temperature value, and the formation thermal conductivity at a preset pressure value and a preset temperature value are calculated, including:

[0059] S401: Calculate the density of the drilling fluid at the preset temperature value and the preset pressure value according to the test temperature value, the test pressure value, the drilling fluid density at the test temperature value and the test pressure value, the preset temperature, and the preset pressure.

[0060] In one embodiment of the present invention, the drilling fluid density at the preset temperature value and the preset pressure value is calculated according to the test temperature value, the test pressure value, the drilling fluid density at the test temperature value and the test pressure value, the preset temperature, and the preset pressure, including:

[0061] The following formula is used to calculate the drilling fluid density at the preset temperature value and the preset pressure value based on the test temperature value, the test pressure value, the drilling fluid density at the test temperature value and the test pressure value, the preset temperature, and the preset pressure:

[0062] ρ=ρ0+[A1(T-T0)+A2(P-P0)]

[0063] Among them, ρ is the drilling fluid density at the preset temperature value and the preset pressure value, ρ0 is the drilling fluid density at the test temperature value and the test pressure value, T0 is the test temperature value, P0 is the test pressure value, T is the preset temperature value, P is the preset pressure value, and A1 and A2 are drilling fluid characteristic parameters.

[0064] Specifically, the drilling fluid characteristic parameters are constant values ​​determined historically based on the drilling fluid characteristics. In one embodiment of the present invention, A1=-0.71603, A2=4.12172×10 -7 .

[0065] S402: Calculate the yield stress of the drilling fluid at a preset temperature value and a preset pressure value according to the test temperature value, the test pressure value, and the yield stress of the drilling fluid at the test temperature value and the test pressure value.

[0066] In one embodiment of the present invention, the yield stress of the drilling fluid at a preset temperature value and a preset pressure value is calculated according to the test temperature value, the test pressure value, and the yield stress of the drilling fluid at the test temperature value and the test pressure value, including:

[0067] The following formula is used to calculate the yield stress of the drilling fluid at the preset temperature value and the preset pressure value based on the test temperature value, the test pressure value, and the yield stress of the drilling fluid at the test temperature value and the test pressure value:

[0068]

[0069] Among them, τ y is the yield stress of the drilling fluid at the preset temperature and pressure, τ y,0 is the yield stress of drilling fluid under test temperature and test pressure, and α1, β1, γ1, δ1, and ε1 are the characteristic parameters of drilling fluid.

[0070] In one embodiment of the present invention, α1 = -1.01229 × 10 -2 , β1=-1.03560×10 -4,γ1=5.08647×10 -3 ,δ1=-5.92252×10 -5 , ε1=1.28569×10 -4 .

[0071] S403: Calculate the plastic viscosity of the drilling fluid at the preset temperature value and the preset pressure value according to the preset test temperature value, the test pressure value, and the plastic viscosity of the drilling fluid at the test temperature value and the test pressure value.

[0072] In one embodiment of the present invention, according to the preset test temperature value, the test pressure value, and the plastic viscosity of the drilling fluid at the test temperature value and the test pressure value, the plastic viscosity of the drilling fluid at the preset temperature value and the preset pressure value is calculated, including:

[0073] The following formula is used to calculate the plastic viscosity of the drilling fluid at the preset temperature value and the preset pressure value according to the preset test temperature value, the test pressure value, and the plastic viscosity of the drilling fluid at the test temperature value and the test pressure value:

[0074]

[0075] Among them, μ p is the plastic viscosity of the drilling fluid at the preset temperature and pressure, μ p,0 is the plastic viscosity of the drilling fluid at the test temperature and test pressure values, and α2, β2, γ2, δ2, and ε2 are the characteristic parameters of the drilling fluid.

[0076] In one embodiment of the present invention, α2 = -3.26140 × 10 -2 , β2=2.38849×10 -4 ,γ2=1.98767×10 -2 ,δ2=4.75823×10 -5 ,ε2=-2.44399×10 -4 .

[0077] S404: Calculate the specific heat capacity of the drilling fluid at the preset temperature value and the preset pressure value according to the preset temperature value.

[0078] In one embodiment of the present invention, according to the preset temperature value, calculating the specific heat capacity of the drilling fluid at the preset temperature value and the preset pressure value includes:

[0079] The following formula is used to calculate the specific heat capacity of the drilling fluid at the preset temperature value and the preset pressure value according to the preset temperature value:

[0080] c fl =a+bT+cT -2

[0081] Among them, cfl is the specific heat capacity of the drilling fluid at the preset temperature value T, and a, b, c are the characteristic parameters of the drilling fluid.

[0082] In one embodiment of the present invention, a=4114.83672, b=1.12275, c=-14.55477.

[0083] S405: According to the preset temperature value, calculating the thermal conductivity of the formation at the preset temperature value.

[0084] In one embodiment of the present invention, according to the preset temperature value, calculating the formation thermal conductivity at the preset temperature value includes:

[0085] The following formula is used to calculate the formation thermal conductivity at the preset temperature value according to the preset temperature value:

[0086] k e =k e,0 (1+βT)

[0087] Among them, k e is the thermal conductivity of the formation at the preset temperature T, k e,0 is the thermal conductivity of the formation at 0℃, and β3 is the formation characteristic parameter.

[0088] Specifically, the formation characteristic parameter is a constant value determined historically based on the formation characteristics.

[0089] S302: Calculate the calculated temperature value of the grid node according to the energy conservation equation, the specific heat capacity of the drilling fluid at a preset pressure value and a preset temperature value, the thermal conductivity of the formation at a preset pressure value and a preset temperature value, the total heat transfer coefficient of the cement ring, the casing and the annulus fluid, the total heat transfer coefficient of the drill pipe and the fluid in the pipe, the dimensionless temperature, the drilling fluid mass flow rate, the wellbore radius, the inner radius of the drill pipe, the inlet drilling fluid temperature, the surface formation temperature, the geothermal gradient, and the total well depth.

[0090] In one embodiment of the present invention, according to the energy conservation equation, the specific heat capacity of the drilling fluid at a preset pressure value and a preset temperature value, the thermal conductivity of the formation at a preset pressure value and a preset temperature value, the total heat transfer coefficient of the cement sheath, the casing and the annular fluid, the total heat transfer coefficient of the drill pipe and the fluid in the pipe, the dimensionless temperature, the drilling fluid mass flow rate, the borehole radius, the inner radius of the drill pipe, the inlet drilling fluid temperature, the surface formation temperature, the geothermal gradient, and the total well depth, the calculated temperature value of the grid node is calculated, including: using the energy conservation equation, the specific heat capacity of the drilling fluid at a preset pressure value and a preset temperature value, and the thermal conductivity of the formation at a preset pressure value and a preset temperature value to obtain the following calculated temperature value calculation formula, and according to the total heat transfer coefficient of the cement sheath, the casing and the annular fluid, the total heat transfer coefficient of the drill pipe and the fluid in the pipe, the dimensionless temperature, the drilling fluid mass flow rate, the borehole radius, the inner radius of the drill pipe, the inlet drilling fluid temperature, the surface formation temperature, the geothermal gradient, and the total well depth, the following calculated temperature value calculation formula is used to determine the calculated temperature value:

[0091]

[0092] Among them, T cal To calculate the temperature value, the boundary condition is T t | z＝0 =T ti , T ti is the inlet drilling fluid temperature, T es is the surface temperature, g G is the geothermal gradient, L is the total well depth c fl is the specific heat capacity of the drilling fluid at the preset pressure and temperature, k e is the formation thermal conductivity at the test temperature and pressure, U a is the total heat transfer coefficient between cement sheath, casing and annulus fluid, U t is the total heat transfer coefficient between the drill pipe and the fluid in the pipe, T D is the dimensionless temperature, w is the drilling fluid mass flow rate, r w is the wellbore radius, r t is the inner radius of the drill pipe, and π is the pi.

[0093] Specifically, the energy conservation equation of drilling fluid is:

[0094]

[0095] Among them, T cal To calculate the temperature value (℃), T t is the drilling fluid temperature in the pipe (℃), T ei is the formation temperature (℃), c flis the specific heat capacity of the drilling fluid at the preset pressure and temperature (J / (kg·℃)), c e is the formation specific heat capacity (J / (kg·℃)), k e is the formation thermal conductivity (W / (m·℃)) at the test temperature and pressure, U a is the total heat transfer coefficient between cement sheath, casing and annulus fluid (W / (m 2 ℃)), U t is the total heat transfer coefficient between the drill pipe and the fluid in the pipe (W / (m 2 ℃)), T D is the dimensionless temperature, t D is the dimensionless cycle time, t is the cycle time (h), w is the drilling fluid mass flow rate (kg / s), ρ e is the formation density (kg / m 3 ), z is the well depth (m), r w is the wellbore radius (m), r t is the inner radius of the drill pipe (m), and π is the pi.

[0096] According to the boundary conditions for the above energy conservation equation: T t | z＝0 =R ti , After analysis, the above temperature calculation formula is obtained.

[0097] S303: Determine a temperature error value according to a preset temperature value and a calculated temperature value.

[0098] In one embodiment of the present invention, the temperature error value is determined by using the following formula according to the preset temperature value and the calculated temperature value:

[0099]

[0100] Among them, R cal To calculate the temperature value, T ass is the preset temperature value, and m1 is the temperature error value.

[0101] S304: When the temperature error value is less than the first error threshold, the calculated pressure value of the grid node is calculated according to the mass conservation equation, the momentum conservation equation, the drilling fluid density at the preset temperature value and the preset pressure value, the drilling fluid yield stress at the preset pressure value and the preset temperature value, the drilling fluid plastic viscosity at the preset pressure value and the preset temperature value, the annular drilling fluid flow rate, the annular cross-sectional area, the well inclination angle at the grid node, and the microelement segment length.

[0102] In one embodiment of the present invention, according to the mass conservation equation and the momentum conservation equation, the drilling fluid density at a preset temperature and a preset pressure value, the drilling fluid yield stress at a preset pressure and a preset temperature value, the drilling fluid plastic viscosity at a preset pressure and a preset temperature value, the annular drilling fluid flow rate, the annular cross-sectional area, the well inclination angle at the grid node, and the micro-element segment length, the calculated pressure value of the grid node is calculated, including: the fluctuating pressure per unit length of the annular drilling fluid obtained according to the drilling fluid density at a preset temperature and a preset pressure value, the drilling fluid yield stress at a preset pressure and a preset temperature value, and the drilling fluid plastic viscosity at a preset pressure and a preset temperature value; according to the mass conservation equation and the momentum conservation equation, the following calculated pressure value calculation formula is obtained; according to the drilling fluid density at a preset temperature and a preset pressure value, the annular drilling fluid flow rate, the annular cross-sectional area, the well inclination angle at the grid node, and the micro-element segment length, the calculated pressure value is determined using the following calculated pressure value calculation formula:

[0103]

[0104] Among them, P cal To calculate the pressure value, i represents the i-th grid node, ρ is the drilling fluid density at the preset temperature and pressure, V a is the annular drilling fluid velocity, A a is the cross-sectional area of ​​the annulus (m 2 ), P a To calculate the pressure value, P Surge,a is the fluctuating pressure per unit length of the annular drilling fluid, g is the gravitational acceleration, θ is the well inclination angle at the grid node, and Δz is the length of the microelement segment.

[0105] Specifically, the mass conservation equation for drilling fluid flow is:

[0106]

[0107] Where ρ is the drilling fluid density (kg / m 3 ), V a is the annular drilling fluid velocity (m / s), A a is the cross-sectional area of ​​the annulus (m 2 ), z is the well depth (m)

[0108] The momentum conservation equation for drilling fluid flow is:

[0109]

[0110] Among them, P cal To calculate the pressure value (Pa), P Surge,a is the fluctuating pressure per unit length of the annular drilling fluid (Pa / m), g is the acceleration due to gravity (9.81 m / s2 ), θ is the well inclination angle (°).

[0111] The above mass conservation equation and momentum conservation equation are discretely solved by finite difference method to obtain the calculation formula of pressure value:

[0112]

[0113] In addition, the fluctuating pressure per unit length of the annular drilling fluid obtained according to the drilling fluid density at a preset temperature value and a preset pressure value, the drilling fluid yield stress at a preset pressure value and a preset temperature value, and the drilling fluid plastic viscosity at a preset pressure value and a preset temperature value, for example, includes:

[0114] (1) Fluid flow pattern judgment

[0115] When Re a >Re a,c When Re a ≤Re a,c When , the annular flow is laminar.

[0116]

[0117] Where: Re a is the annular fluid Reynolds number, D w is the inner diameter of the wellbore (m), D d,o is the outer diameter of the drill pipe (m), ρ is the density of the drilling fluid at the preset pressure and temperature (kg / m 3 ), μ p It is the plastic viscosity of the drilling fluid at the preset pressure value and the preset temperature value (Pa·s).

[0118]

[0119] Where: Re a,c is the critical Reynolds number of the annular fluid, τ y is the yield stress of drilling fluid at preset pressure and temperature (Pa), τ w is the yield stress of drilling fluid on the outer wall of drill pipe (Pa), and He is the Hedstrom number.

[0120]

[0121]

[0122] (2) Calculation of laminar flow pressure (Re a ≤Re a,c )

[0123]

[0124] Where: P Surge,a It is the fluctuation pressure per unit length of the annular drilling fluid (Pa / m).

[0125] (3) Calculation of turbulent fluctuation pressure (Re a >Re a,c )

[0126]

[0127] S305: Determine a pressure error value according to a preset pressure value and a calculated pressure value.

[0128] In one embodiment of the present invention, the pressure error value is determined by using the following formula according to the preset pressure value and the calculated pressure value:

[0129]

[0130] Among them, P cal To calculate the pressure value, P ass is the preset pressure value, and m2 is the pressure error value.

[0131] S306: When the pressure error value is less than the second error threshold, the calculated temperature value is determined as the actual pressure value of the grid node, and the process ends.

[0132] S307: When the pressure error value is not less than the second error threshold, the calculated pressure value is used as the preset pressure value, and the process returns to S304.

[0133] S308: When the temperature error value is not less than the first error threshold, return to 301.

[0134] The present invention also provides a device for determining the wellbore pressure of a high-temperature and high-pressure well in an embodiment, as described in the following embodiment. Since the principle of the device to solve the problem is similar to the method for determining the wellbore pressure of a high-temperature and high-pressure well in an embodiment, the implementation of the device can refer to the implementation of the method for determining the wellbore pressure of a high-temperature and high-pressure well in an embodiment, and the repeated parts will not be repeated.

[0135] like Figure 5 FIG. 1 is a schematic diagram of a device for determining the wellbore pressure during tripping of a high-temperature and high-pressure well provided by an embodiment of the present invention, comprising:

[0136] The parameter acquisition module 501 is used to acquire parameter data during the tripping operation; wherein the parameter data at least includes: inclination measurement data and preset well depth;

[0137] The first processing module 502 is used to divide the grid nodes according to the preset number of nodes and the preset well depth, and obtain the micro-element segment length;

[0138] The second processing module 503 is used to determine the well inclination angle at the grid node according to the inclination measurement data;

[0139] The third processing module 504 is used to determine the tripping wellbore pressure under temperature-pressure coupling according to the mass conservation equation, momentum conservation equation, energy conservation equation, microelement segment length, well inclination angle at grid nodes, and parameter data.

[0140] In a possible implementation, the first processing module is specifically used to divide the grid nodes by an equal spacing or variable spacing method according to a preset number of nodes and a total well depth to obtain an initial grid node matrix;

[0141] The well depth of the preset point is added to the initial grid node matrix to obtain the final divided grid node matrix; wherein each element value in the grid node matrix represents the well depth corresponding to a grid node;

[0142] The length of the microelement segment is determined based on the preset number of nodes and the total well depth.

[0143] In a possible implementation, the second processing module is specifically configured to calculate the well inclination angle at the grid node by using a linear interpolation method according to the inclination survey data.

[0144] In a possible implementation, the inclination survey data includes: the well depth and well inclination angle at the inclination survey point above the grid node, and the well depth and well inclination angle at the inclination survey point below the grid node;

[0145] The second processing module is specifically used to calculate the well inclination angle at the grid node by using the well depth and well inclination angle at the inclination measuring point above the grid node and the well depth and well inclination angle at the inclination measuring point below the grid node using the following formula:

[0146]

[0147] Where θ is the well inclination angle at the grid node; D m is the well depth at the grid node; D1 is the well inclination angle at the inclination measurement point above the grid node; D m,1 is the well depth at the inclination point above the grid node; θ2 is the well inclination angle at the inclination point below the grid node; D m,2 It is the well depth at the inclination measurement point below the grid node.

[0148] In a possible implementation, the parameter data also includes: the total heat transfer coefficient of the cement ring, the casing and the annular fluid, the total heat transfer coefficient of the drill pipe and the fluid in the pipe, the dimensionless temperature, the drilling fluid mass flow rate, the wellbore radius, the inner radius of the drill pipe, the test temperature value, the test pressure value, the annular drilling fluid flow rate, the inclination data, the drilling fluid density at the test temperature value and the test pressure value, the drilling fluid yield stress at the test temperature value and the test pressure value, the drilling fluid plastic viscosity at the test temperature value and the test pressure value, the drilling fluid specific heat capacity at the test temperature value and the test pressure value, the inlet drilling fluid temperature, the formation specific heat capacity, the formation thermal conductivity at the test temperature value and the test pressure value, the surface formation temperature, the geothermal gradient and the total well depth, the preset point well depth, the inlet drilling fluid temperature, the drilling fluid flow rate, and the annular cross-sectional area;

[0149] The third processing module is specifically configured to execute the following process for each grid node:

[0150] According to the preset pressure value, preset temperature value and parameter data of the grid node, the drilling fluid density at the preset pressure value and preset temperature value, the drilling fluid yield stress at the preset pressure value and preset temperature value, the drilling fluid plastic viscosity at the preset pressure value and preset temperature value, the drilling fluid specific heat capacity at the preset pressure value and preset temperature value and the formation thermal conductivity at the preset pressure value and preset temperature value are calculated;

[0151] Calculate the calculated temperature value of the grid node according to the energy conservation equation, the specific heat capacity of the drilling fluid at a preset pressure value and a preset temperature value, the thermal conductivity of the formation at a preset pressure value and a preset temperature value, the total heat transfer coefficient of the cement sheath, the casing and the annulus fluid, the total heat transfer coefficient of the drill pipe and the fluid in the pipe, the dimensionless temperature, the drilling fluid mass flow rate, the wellbore radius, the inner radius of the drill pipe, the inlet drilling fluid temperature, the surface formation temperature, the geothermal gradient, and the total well depth;

[0152] Determine the temperature error value according to the preset temperature value and the calculated temperature value;

[0153] When the temperature error value is less than the first error threshold, the calculated pressure value of the grid node is calculated according to the mass conservation equation, the momentum conservation equation, the drilling fluid density at the preset temperature value and the preset pressure value, the drilling fluid yield stress at the preset pressure value and the preset temperature value, the drilling fluid plastic viscosity at the preset pressure value and the preset temperature value, the annular drilling fluid flow rate, the annular cross-sectional area, the well inclination angle at the grid node, and the microelement segment length;

[0154] Determine a pressure error value according to a preset pressure value and a calculated pressure value;

[0155] When the pressure error value is less than a second error threshold, determining the calculated temperature value as the actual pressure value of the grid node;

[0156] When the pressure error value is not less than the second error threshold, the calculated pressure value is used as the preset pressure value, and the step of returning to calculate the calculated pressure value of the grid node according to the mass conservation equation, the momentum conservation equation, the preset temperature value, and the drilling fluid density at the preset pressure value is returned;

[0157] When the temperature error value is not less than the first error threshold, return to the step of calculating the drilling fluid density at the preset pressure value and the preset temperature value, the drilling fluid yield stress at the preset pressure value and the preset temperature value, the drilling fluid plastic viscosity at the preset pressure value and the preset temperature value, the drilling fluid specific heat capacity at the preset pressure value and the preset temperature value, and the formation thermal conductivity at the preset pressure value and the preset temperature value according to the preset pressure value, the preset temperature value and the preset temperature value of the grid node.

[0158] In a possible implementation, the third processing module is specifically used to calculate the drilling fluid density at the preset temperature value and the preset pressure value according to the test temperature value, the test pressure value, the drilling fluid density at the test temperature value and the test pressure value, the preset temperature, and the preset pressure;

[0159] Calculate the yield stress of the drilling fluid at a preset temperature value and a preset pressure value according to the test temperature value, the test pressure value, and the yield stress of the drilling fluid at the test temperature value and the test pressure value;

[0160] Calculating the plastic viscosity of the drilling fluid at the preset temperature value and the preset pressure value according to the preset test temperature value, the test pressure value, and the plastic viscosity of the drilling fluid at the test temperature value and the test pressure value;

[0161] According to the preset temperature value, calculating the specific heat capacity of the drilling fluid at the preset temperature value and the preset pressure value;

[0162] According to the preset temperature value, the thermal conductivity of the formation at the preset temperature value is calculated.

[0163] In a possible implementation, the third processing module is specifically configured to calculate the drilling fluid density using the following formula:

[0164] ρ=ρ0+[A19T-T0)+A2(P-P0)]

[0165] Among them, ρ is the drilling fluid density at the preset temperature value and the preset pressure value, ρ0 is the drilling fluid density at the test temperature value and the test pressure value, T0 is the test temperature value, P0 is the test pressure value, T is the preset temperature value, P is the preset pressure value, and A1 and A2 are drilling fluid characteristic parameters.

[0166] In a possible implementation, the third processing module is specifically configured to calculate the drilling fluid yield stress using the following formula:

[0167]

[0168] Among them, τ y is the yield stress of the drilling fluid at the preset temperature and pressure, τ y,0 is the yield stress of drilling fluid under test temperature and test pressure, and α1, β1, γ1, δ1, and ε1 are the characteristic parameters of drilling fluid.

[0169] In a possible implementation, the third processing module is specifically used to calculate the plastic viscosity of the drilling fluid at a preset temperature value and a preset pressure value according to the preset test temperature value, the test pressure value, and the plastic viscosity of the drilling fluid at the test temperature value and the test pressure value using the following formula:

[0170]

[0171] Among them, μ p is the plastic viscosity of the drilling fluid at the preset temperature and pressure, μ p,0 is the plastic viscosity of the drilling fluid at the test temperature and test pressure values, and α2, β2, γ2, δ2, and ε2 are the characteristic parameters of the drilling fluid.

[0172] In a possible implementation, the third processing module is specifically configured to calculate the specific heat capacity of the drilling fluid at a preset temperature value and a preset pressure value according to the preset temperature value using the following formula:

[0173] c fl =a+bT+cT -2

[0174] Among them, c fl is the specific heat capacity of the drilling fluid at the preset temperature value T, and a, b, c are the characteristic parameters of the drilling fluid.

[0175] In a possible implementation, the third processing module is specifically configured to calculate the formation thermal conductivity at a preset temperature value according to the preset temperature value using the following formula:

[0176] k e =k e,0 (1+βT)

[0177] Among them, k e is the thermal conductivity of the formation at the preset temperature T, k e,0 is the thermal conductivity of the formation at 0℃, and β3 is the formation characteristic parameter.

[0178] In a possible implementation, the third processing module is specifically used to use the energy conservation equation, the specific heat capacity of the drilling fluid at a preset pressure value and a preset temperature value, and the thermal conductivity of the formation at a preset pressure value and a preset temperature value to obtain the following calculated temperature value calculation formula, and the calculated temperature value is determined by using the following calculated temperature value calculation formula according to the total heat transfer coefficient of the cement ring, the casing, and the annular fluid, the total heat transfer coefficient of the drill pipe and the fluid in the pipe, the dimensionless temperature, the drilling fluid mass flow rate, the wellbore radius, the inner radius of the drill pipe, the inlet drilling fluid temperature, the surface formation temperature, the geothermal gradient, and the total well depth:

[0179]

[0180] Among them, T cal To calculate the temperature value, the boundary condition is T t | z＝0 =T ti , T ti is the inlet drilling fluid temperature, T es is the surface temperature, g G is the geothermal gradient, L is the total well depth, c fl is the specific heat capacity of the drilling fluid at the preset pressure and temperature, k e is the formation thermal conductivity at the test temperature and pressure, U a is the total heat transfer coefficient between cement sheath, casing and annulus fluid, U t is the total heat transfer coefficient between the drill pipe and the fluid in the pipe, T D is the dimensionless temperature, w is the drilling fluid mass flow rate, r w is the wellbore radius, r t is the inner radius of the drill pipe, and π is the pi.

[0181] In a possible implementation manner, the third processing module is specifically configured to determine the temperature error value using the following formula according to the preset temperature value and the calculated temperature value:

[0182]

[0183] Among them, T cal To calculate the temperature value, T ass is the preset temperature value, and m1 is the temperature error value.

[0184] In a possible implementation, the third processing module is specifically used to obtain the fluctuating pressure per unit length of the annular drilling fluid according to the drilling fluid density at a preset temperature value and a preset pressure value, the drilling fluid yield stress at a preset pressure value and a preset temperature value, and the drilling fluid plastic viscosity at a preset pressure value and a preset temperature value; according to the mass conservation equation and the momentum conservation equation, the following calculated pressure value calculation formula is obtained; according to the drilling fluid density at a preset temperature value and a preset pressure value, the annular drilling fluid flow rate, the annular cross-sectional area, the well inclination angle at the grid node, and the microelement segment length, the calculated pressure value is determined using the following calculated pressure value calculation formula:

[0185]

[0186] Among them, P cal To calculate the pressure value, i represents the i-th grid node, ρ is the drilling fluid density at the preset temperature and pressure, V a is the annular drilling fluid velocity, A a is the cross-sectional area of ​​the annulus (m 2 ), P a To calculate the pressure value, P Surge,a is the fluctuating pressure per unit length of the annular drilling fluid, g is the gravitational acceleration, θ is the well inclination angle at the grid node, and Δz is the length of the microelement segment.

[0187] In a possible implementation manner, the third processing module is specifically configured to determine the pressure error value using the following formula according to the preset pressure value and the calculated pressure value:

[0188]

[0189] Among them, P cal To calculate the pressure value, P ass is the preset pressure value, and m2 is the pressure error value.

[0190] Based on the above invention concept, Figure 6 As shown, the present invention also proposes a computer device 600, including a memory 610, a processor 620, and a computer program 630 stored in the memory 610 and executable on the processor 620, wherein the processor 620 implements the aforementioned method for determining the wellbore pressure during drilling of a high-temperature and high-pressure well when executing the computer program 630.

[0191] An embodiment of the present invention further provides a computer-readable storage medium storing a computer program, which, when executed by a processor, implements the above-mentioned method for determining the wellbore pressure during tripping and drilling of a high-temperature and high-pressure well.

[0192] An embodiment of the present invention further provides a computer program product, which includes a computer program. When the computer program is executed by a processor, the method for determining the wellbore pressure during tripping and drilling of a high-temperature and high-pressure well is implemented.

[0193] In the embodiment of the present invention, parameter data in the process of drilling operation is obtained; wherein the parameter data at least includes: inclination data and well depth of preset points; grid nodes are divided according to the number of preset nodes and the well depth of preset points, and the length of micro-element segments is obtained; the well inclination angle at the grid node is determined according to the inclination data; the wellbore pressure of drilling under temperature-pressure coupling is determined according to the mass conservation equation, momentum conservation equation, and energy conservation equation, micro-element segment length, well inclination angle at the grid node, and parameter data. In this way, the effect of high temperature and high pressure on the properties of drilling fluid is taken into account, the calculation accuracy of wellbore pressure of drilling is improved, the accuracy of wellbore pressure prediction is improved, and the accurate regulation of wellbore pressure of drilling in high temperature and high pressure wells is effectively met, ensuring safe and efficient drilling operation. .

[0194] Those skilled in the art will appreciate that embodiments of the present invention may be provided as methods, systems, or computer program products. Therefore, the present invention may take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Moreover, the present invention may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0195] The present invention is described with reference to flowcharts and / or block diagrams of methods, devices (systems), and computer program products according to embodiments of the present invention. It should be understood that each process and / or block in the flowchart and / or block diagram, as well as the combination of processes and / or blocks in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0196] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0197] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process in the computer or other programmable device. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0198] The specific embodiments described above further illustrate the objectives, technical solutions and beneficial effects of the present invention in detail. It should be understood that the above description is only a specific embodiment of the present invention and is not intended to limit the scope of protection of the present invention. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present invention should be included in the scope of protection of the present invention.

Claims

1. A method for determining the wellbore pressure during tripping of a high-temperature and high-pressure well, characterized in that: include: Obtaining parameter data during the tripping operation; wherein the parameter data at least includes: inclination measurement data and well depth at a preset point; According to the preset number of nodes and the preset well depth, the grid nodes are divided and the length of the microelement segment is obtained; Determine the well inclination angle at the grid node according to the inclination survey data; The tripping wellbore pressure under temperature-pressure coupling is determined based on the mass conservation equation, momentum conservation equation, energy conservation equation, microelement segment length, well inclination angle at grid nodes, and parameter data.

2. The method for determining the wellbore pressure during tripping of a high-temperature and high-pressure well according to claim 1, characterized in that: According to the preset number of nodes and the preset well depth, the grid nodes are divided and the length of the microelement segment is obtained, including: According to the preset number of nodes and the total well depth, the grid nodes are divided by equal spacing or variable spacing to obtain an initial grid node matrix; The well depth of the preset point is added to the initial grid node matrix to obtain the final divided grid node matrix; wherein each element value in the grid node matrix represents the well depth corresponding to a grid node; The length of the microelement segment is determined based on the preset number of nodes and the total well depth.

3. The method for determining the wellbore pressure during tripping of a high-temperature and high-pressure well according to claim 1, characterized in that: Preset point well depth, including: The well depth at the casing shoe, the well depth at the drill string diameter change point, and the well depth at the abnormal formation pressure.

4. The method for determining the wellbore pressure during tripping of a high-temperature and high-pressure well according to claim 1, characterized in that: Determine the well inclination angle at the grid node based on the inclination survey data, including: Based on the inclination data, the linear interpolation method is used to calculate the well inclination angle at the grid node.

5. The method for determining the wellbore pressure during tripping of a high-temperature and high-pressure well according to claim 3, characterized in that: The inclination survey data include: the well depth and well inclination angle at the inclination survey point above the grid node, and the well depth and well inclination angle at the inclination survey point below the grid node; According to the inclination data, the linear interpolation method is used to calculate the well inclination angle at the grid node, including: The well depth and well inclination angle at the inclination point above the grid node, and the well depth and well inclination angle at the inclination point below the grid node are calculated using the following formula: Where θ is the well inclination angle at the grid node; D m is the well depth at the grid node; θ1 is the well inclination angle at the inclination measurement point above the grid node; D m,1 is the well depth at the inclination point above the grid node; θ2 is the well inclination angle at the inclination point below the grid node; D m,2 It is the well depth at the inclination measurement point below the grid node.

6. The method for determining the wellbore pressure during tripping of a high-temperature and high-pressure well according to claim 1, characterized in that: The parameter data also include: total heat transfer coefficient of cement sheath, casing and annular fluid, total heat transfer coefficient of drill pipe and fluid in pipe, dimensionless temperature, drilling fluid mass flow rate, borehole radius, inner radius of drill pipe, test temperature value, test pressure value, annular drilling fluid flow rate, inclination data, drilling fluid density at test temperature value and test pressure value, drilling fluid yield stress at test temperature value and test pressure value, drilling fluid plastic viscosity at test temperature value and test pressure value, drilling fluid specific heat capacity at test temperature value and test pressure value, inlet drilling fluid temperature, formation specific heat capacity, formation thermal conductivity at test temperature value and test pressure value, surface formation temperature, geothermal gradient and total well depth, preset point well depth, inlet drilling fluid temperature, drilling fluid flow rate, annular cross-sectional area; According to the mass conservation equation, momentum conservation equation, energy conservation equation, microelement length, well inclination angle at grid node, and parameter data, the tripping wellbore pressure under temperature-pressure coupling is determined, including: For each grid node, perform the following process: According to the preset pressure value, preset temperature value and parameter data of the grid node, the drilling fluid density at the preset pressure value and preset temperature value, the drilling fluid yield stress at the preset pressure value and preset temperature value, the drilling fluid plastic viscosity at the preset pressure value and preset temperature value, the drilling fluid specific heat capacity at the preset pressure value and preset temperature value and the formation thermal conductivity at the preset pressure value and preset temperature value are calculated; Calculate the calculated temperature value of the grid node according to the energy conservation equation, the specific heat capacity of the drilling fluid at a preset pressure value and a preset temperature value, the thermal conductivity of the formation at a preset pressure value and a preset temperature value, the total heat transfer coefficient of the cement sheath, the casing and the annulus fluid, the total heat transfer coefficient of the drill pipe and the fluid in the pipe, the dimensionless temperature, the drilling fluid mass flow rate, the wellbore radius, the inner radius of the drill pipe, the inlet drilling fluid temperature, the surface formation temperature, the geothermal gradient, and the total well depth; Determine the temperature error value according to the preset temperature value and the calculated temperature value; When the temperature error value is less than the first error threshold, the calculated pressure value of the grid node is calculated according to the mass conservation equation, the momentum conservation equation, the drilling fluid density at the preset temperature value and the preset pressure value, the drilling fluid yield stress at the preset pressure value and the preset temperature value, the drilling fluid plastic viscosity at the preset pressure value and the preset temperature value, the annular drilling fluid flow rate, the annular cross-sectional area, the well inclination angle at the grid node, and the microelement segment length; Determine a pressure error value according to a preset pressure value and a calculated pressure value; When the pressure error value is less than a second error threshold, determining the calculated temperature value as the actual pressure value of the grid node; When the pressure error value is not less than the second error threshold, the calculated pressure value is used as the preset pressure value, and the step of returning to calculate the calculated pressure value of the grid node according to the mass conservation equation, the momentum conservation equation, the preset temperature value, and the drilling fluid density at the preset pressure value is returned; When the temperature error value is not less than the first error threshold, return to the step of calculating the drilling fluid density at the preset pressure value and the preset temperature value, the drilling fluid yield stress at the preset pressure value and the preset temperature value, the drilling fluid plastic viscosity at the preset pressure value and the preset temperature value, the drilling fluid specific heat capacity at the preset pressure value and the preset temperature value, and the formation thermal conductivity at the preset pressure value and the preset temperature value according to the preset pressure value, the preset temperature value and the preset temperature value of the grid node.

7. The method for determining the wellbore pressure during tripping of a high-temperature and high-pressure well according to claim 6, characterized in that: Calculating the drilling fluid density at a preset pressure value and a preset temperature value, the drilling fluid yield stress at a preset pressure value and a preset temperature value, the drilling fluid plastic viscosity at a preset pressure value and a preset temperature value, the drilling fluid specific heat capacity at a preset pressure value and a preset temperature value, and the formation thermal conductivity at a preset pressure value and a preset temperature value, including: Calculate the density of the drilling fluid at the preset temperature value and the preset pressure value according to the test temperature value, the test pressure value, the density of the drilling fluid at the test temperature value and the test pressure value, the preset temperature, and the preset pressure; Calculate the yield stress of the drilling fluid at a preset temperature value and a preset pressure value according to the test temperature value, the test pressure value, and the yield stress of the drilling fluid at the test temperature value and the test pressure value; Calculating the plastic viscosity of the drilling fluid at the preset temperature value and the preset pressure value according to the preset test temperature value, the test pressure value, and the plastic viscosity of the drilling fluid at the test temperature value and the test pressure value; According to the preset temperature value, calculating the specific heat capacity of the drilling fluid at the preset temperature value and the preset pressure value; According to the preset temperature value, the thermal conductivity of the formation at the preset temperature value is calculated.

8. The method for determining the wellbore pressure during tripping of a high-temperature and high-pressure well according to claim 7, characterized in that: According to the test temperature value, the test pressure value, the drilling fluid density at the test temperature value and the test pressure value, the preset temperature, and the preset pressure, the drilling fluid density at the preset temperature value and the preset pressure value is calculated, including: The following formula is used to calculate the drilling fluid density at the preset temperature value and the preset pressure value based on the test temperature value, the test pressure value, the drilling fluid density at the test temperature value and the test pressure value, the preset temperature, and the preset pressure: ρ=ρ0+[A1(T-T0)+A2(P-P0)] Among them, ρ is the drilling fluid density at the preset temperature value and the preset pressure value, ρ0 is the drilling fluid density at the test temperature value and the test pressure value, T0 is the test temperature value, P0 is the test pressure value, T is the preset temperature value, P is the preset pressure value, and A1 and A2 are drilling fluid characteristic parameters.

9. The method for determining the wellbore pressure during tripping of a high-temperature and high-pressure well according to claim 7, characterized in that: According to the test temperature value, the test pressure value, and the yield stress of the drilling fluid under the test temperature value and the test pressure value, the yield stress of the drilling fluid under the preset temperature value and the preset pressure value is calculated, including: The following formula is used to calculate the yield stress of the drilling fluid at the preset temperature value and the preset pressure value based on the test temperature value, the test pressure value, and the yield stress of the drilling fluid at the test temperature value and the test pressure value: Among them, τ y is the yield stress of the drilling fluid at the preset temperature and pressure, τ y,0 is the yield stress of drilling fluid under test temperature and test pressure, and α1, β1, γ1, δ1, and ε1 are the characteristic parameters of drilling fluid.

10. The method for determining the wellbore pressure during tripping of a high-temperature and high-pressure well according to claim 7, characterized in that: According to the preset test temperature value, the test pressure value, and the plastic viscosity of the drilling fluid at the test temperature value and the test pressure value, the plastic viscosity of the drilling fluid at the preset temperature value and the preset pressure value is calculated, including: The following formula is used to calculate the plastic viscosity of the drilling fluid at the preset temperature value and the preset pressure value according to the preset test temperature value, the test pressure value, and the plastic viscosity of the drilling fluid at the test temperature value and the test pressure value: Among them, μ p is the plastic viscosity of the drilling fluid at the preset temperature and pressure, μ p,0 is the plastic viscosity of the drilling fluid at the test temperature and test pressure values, and α2, β2, γ2, δ2, and ε2 are the characteristic parameters of the drilling fluid.

11. The method for determining the wellbore pressure during tripping of a high-temperature and high-pressure well according to claim 7, characterized in that: According to the preset temperature value, the specific heat capacity of the drilling fluid at the preset temperature value and the preset pressure value is calculated, including: The following formula is used to calculate the specific heat capacity of the drilling fluid at the preset temperature value and the preset pressure value according to the preset temperature value: c fl =a+bT+cT -2 Among them, c fl is the specific heat capacity of the drilling fluid at the preset temperature value T, and a, b, c are the characteristic parameters of the drilling fluid.

12. The method for determining the wellbore pressure during tripping of a high-temperature and high-pressure well according to claim 7, characterized in that: According to the preset temperature value, the formation thermal conductivity under the preset temperature value is calculated, including: The following formula is used to calculate the formation thermal conductivity at the preset temperature value according to the preset temperature value: k e =k e,0 (1+βT) Among them, k e is the thermal conductivity of the formation at the preset temperature T, k e,0 is the thermal conductivity of the formation at 0℃, and β3 is the formation characteristic parameter.

13. The method for determining the wellbore pressure during tripping of a high-temperature and high-pressure well according to claim 6, characterized in that: According to the energy conservation equation, the specific heat capacity of the drilling fluid at the preset pressure value and the preset temperature value, the thermal conductivity of the formation at the preset pressure value and the preset temperature value, the total heat transfer coefficient of the cement sheath, the casing and the annulus fluid, the total heat transfer coefficient of the drill pipe and the fluid in the pipe, the dimensionless temperature, the drilling fluid mass flow rate, the wellbore radius, the inner radius of the drill pipe, the inlet drilling fluid temperature, the surface formation temperature, the geothermal gradient, and the total well depth, the calculated temperature value of the grid node is calculated, including: The energy conservation equation, the specific heat capacity of the drilling fluid at the preset pressure value and the preset temperature value, and the thermal conductivity of the formation at the preset pressure value and the preset temperature value are used to obtain the following calculation temperature value calculation formula. The calculation temperature value is determined by the following calculation temperature value calculation formula according to the total heat transfer coefficient of the cement sheath, the casing and the annulus fluid, the total heat transfer coefficient of the drill pipe and the fluid in the pipe, the dimensionless temperature, the drilling fluid mass flow rate, the wellbore radius, the inner radius of the drill pipe, the inlet drilling fluid temperature, the surface formation temperature, the geothermal gradient, and the total well depth: Among them, T cal To calculate the temperature value, the boundary condition is T t | z＝0 =T ti , T ti is the inlet drilling fluid temperature, T es is the surface temperature, g G is the geothermal gradient, L is the total well depth, c fl is the specific heat capacity of the drilling fluid at the preset pressure and temperature, k e is the formation thermal conductivity at the test temperature and pressure, U a is the total heat transfer coefficient between cement sheath, casing and annulus fluid, U t is the total heat transfer coefficient between the drill pipe and the fluid in the pipe, T D is the dimensionless temperature, w is the drilling fluid mass flow rate, r w is the wellbore radius, r t is the inner radius of the drill pipe, and π is the pi.

14. The method for determining the wellbore pressure during tripping of a high-temperature and high-pressure well according to claim 6, characterized in that: The temperature error value is determined according to the preset temperature value and the calculated temperature value, including: The temperature error value is determined by the following formula based on the preset temperature value and the calculated temperature value: Among them, T cal To calculate the temperature value, T ass is the preset temperature value, and m1 is the temperature error value.

15. The method for determining the wellbore pressure during tripping of a high-temperature and high-pressure well according to claim 6, characterized in that: According to the mass conservation equation and momentum conservation equation, the drilling fluid density at the preset temperature value and the preset pressure value, the drilling fluid yield stress at the preset pressure value and the preset temperature value, the drilling fluid plastic viscosity at the preset pressure value and the preset temperature value, the annular drilling fluid flow rate, the annular cross-sectional area, the well inclination angle at the grid node, and the microelement segment length, the calculated pressure value of the grid node is calculated, including: Fluctuating pressure per unit length of the annular drilling fluid obtained according to the drilling fluid density at a preset temperature value and a preset pressure value, the drilling fluid yield stress at a preset pressure value and a preset temperature value, and the drilling fluid plastic viscosity at a preset pressure value and a preset temperature value; According to the mass conservation equation and momentum conservation equation, the following calculation formula for the calculated pressure value is obtained. According to the drilling fluid density, annular drilling fluid velocity, annular cross-sectional area, well inclination angle at the grid node, and micro-element length at the preset temperature and preset pressure values, the calculated pressure value is determined using the following calculation formula for the calculated pressure value: Among them, P cal To calculate the pressure value, i represents the i-th grid node, ρ is the drilling fluid density at the preset temperature and pressure, V a is the annular drilling fluid velocity, A a is the cross-sectional area of ​​the annulus, P a To calculate the pressure value, P Surge,a is the fluctuating pressure per unit length of the annular drilling fluid, g is the gravitational acceleration, θ is the well inclination angle at the grid node, and Δz is the length of the microelement segment.

16. The method for determining the wellbore pressure during tripping of a high-temperature and high-pressure well according to claim 6, characterized in that: The pressure error value is determined according to the preset pressure value and the calculated pressure value, including: The pressure error value is determined based on the preset pressure value and the calculated pressure value using the following formula: Among them, P cal To calculate the pressure value, P ass is the preset pressure value, and m2 is the pressure error value.

17. A device for determining the wellbore pressure during tripping of a high-temperature and high-pressure well, characterized in that: include: The parameter acquisition module is used to acquire parameter data during the tripping operation; wherein the parameter data at least includes: inclination measurement data and well depth at a preset point; The first processing module is used to divide the grid nodes and obtain the micro-element segment length according to the preset node number and the preset point well depth; The second processing module is used to determine the well inclination angle at the grid node according to the inclination measurement data; The third processing module is used to determine the tripping wellbore pressure under temperature-pressure coupling according to the mass conservation equation, momentum conservation equation, energy conservation equation, microelement segment length, well inclination angle at grid nodes, and parameter data.

18. The apparatus for determining the wellbore pressure during tripping of a high-temperature and high-pressure well according to claim 17, characterized in that: The first processing module is specifically used to divide the grid nodes by equal spacing or variable spacing according to the preset number of nodes and the total well depth to obtain an initial grid node matrix; The well depth of the preset point is added to the initial grid node matrix to obtain the final divided grid node matrix; wherein each element value in the grid node matrix represents the well depth corresponding to a grid node; The length of the microelement segment is determined based on the preset number of nodes and the total well depth.

19. The apparatus for determining the wellbore pressure during tripping of a high-temperature and high-pressure well according to claim 17, characterized in that: The second processing module is specifically used to calculate the well inclination angle at the grid node by using the linear interpolation method according to the inclination survey data.

20. The apparatus for determining the wellbore pressure during tripping of a high-temperature and high-pressure well according to claim 19, characterized in that: The inclination survey data include: the well depth and well inclination angle at the inclination survey point above the grid node, and the well depth and well inclination angle at the inclination survey point below the grid node; The second processing module is specifically used to calculate the well inclination angle at the grid node by using the well depth and well inclination angle at the inclination measuring point above the grid node and the well depth and well inclination angle at the inclination measuring point below the grid node using the following formula: Where θ is the well inclination angle at the grid node; D m is the well depth at the grid node; θ1 is the well inclination angle at the inclination measurement point above the grid node; D m,1 is the well depth at the inclination point above the grid node; θ2 is the well inclination angle at the inclination point below the grid node; D m,2 It is the well depth at the inclination measurement point below the grid node.

21. The apparatus for determining the wellbore pressure during tripping of a high-temperature and high-pressure well according to claim 17, characterized in that: The parameter data also include: total heat transfer coefficient of cement sheath, casing and annular fluid, total heat transfer coefficient of drill pipe and fluid in pipe, dimensionless temperature, drilling fluid mass flow rate, borehole radius, inner radius of drill pipe, test temperature value, test pressure value, annular drilling fluid flow rate, inclination data, drilling fluid density at test temperature value and test pressure value, drilling fluid yield stress at test temperature value and test pressure value, drilling fluid plastic viscosity at test temperature value and test pressure value, drilling fluid specific heat capacity at test temperature value and test pressure value, inlet drilling fluid temperature, formation specific heat capacity, formation thermal conductivity at test temperature value and test pressure value, surface formation temperature, geothermal gradient and total well depth, preset point well depth, inlet drilling fluid temperature, drilling fluid flow rate, annular cross-sectional area; The third processing module is specifically configured to execute the following process for each grid node: According to the preset pressure value, preset temperature value and parameter data of the grid node, the drilling fluid density at the preset pressure value and preset temperature value, the drilling fluid yield stress at the preset pressure value and preset temperature value, the drilling fluid plastic viscosity at the preset pressure value and preset temperature value, the drilling fluid specific heat capacity at the preset pressure value and preset temperature value and the formation thermal conductivity at the preset pressure value and preset temperature value are calculated; Calculate the calculated temperature value of the grid node according to the energy conservation equation, the specific heat capacity of the drilling fluid at a preset pressure value and a preset temperature value, the thermal conductivity of the formation at a preset pressure value and a preset temperature value, the total heat transfer coefficient of the cement sheath, the casing and the annulus fluid, the total heat transfer coefficient of the drill pipe and the fluid in the pipe, the dimensionless temperature, the drilling fluid mass flow rate, the wellbore radius, the inner radius of the drill pipe, the inlet drilling fluid temperature, the surface formation temperature, the geothermal gradient, and the total well depth; Determine the temperature error value according to the preset temperature value and the calculated temperature value; When the temperature error value is less than the first error threshold, the calculated pressure value of the grid node is calculated according to the mass conservation equation, the momentum conservation equation, the drilling fluid density at the preset temperature value and the preset pressure value, the drilling fluid yield stress at the preset pressure value and the preset temperature value, the drilling fluid plastic viscosity at the preset pressure value and the preset temperature value, the annular drilling fluid flow rate, the annular cross-sectional area, the well inclination angle at the grid node, and the microelement segment length; Determine a pressure error value according to a preset pressure value and a calculated pressure value; When the pressure error value is less than a second error threshold, determining the calculated temperature value as the actual pressure value of the grid node; When the pressure error value is not less than the second error threshold, the calculated pressure value is used as the preset pressure value, and the step of returning to calculate the calculated pressure value of the grid node according to the mass conservation equation, the momentum conservation equation, the preset temperature value, and the drilling fluid density at the preset pressure value is returned; When the temperature error value is not less than the first error threshold, return to the step of calculating the drilling fluid density at the preset pressure value and the preset temperature value, the drilling fluid yield stress at the preset pressure value and the preset temperature value, the drilling fluid plastic viscosity at the preset pressure value and the preset temperature value, the drilling fluid specific heat capacity at the preset pressure value and the preset temperature value, and the formation thermal conductivity at the preset pressure value and the preset temperature value according to the preset pressure value, the preset temperature value and the preset temperature value of the grid node.

22. The apparatus for determining the wellbore pressure during tripping of a high-temperature and high-pressure well according to claim 21, characterized in that: The third processing module is specifically used to calculate the drilling fluid density at the preset temperature value and the preset pressure value according to the test temperature value, the test pressure value, the drilling fluid density at the test temperature value and the test pressure value, the preset temperature, and the preset pressure; Calculate the yield stress of the drilling fluid at a preset temperature value and a preset pressure value according to the test temperature value, the test pressure value, and the yield stress of the drilling fluid at the test temperature value and the test pressure value; Calculating the plastic viscosity of the drilling fluid at the preset temperature value and the preset pressure value according to the preset test temperature value, the test pressure value, and the plastic viscosity of the drilling fluid at the test temperature value and the test pressure value; According to the preset temperature value, calculating the specific heat capacity of the drilling fluid at the preset temperature value and the preset pressure value; According to the preset temperature value, the thermal conductivity of the formation at the preset temperature value is calculated.

23. The apparatus for determining the wellbore pressure during tripping of a high-temperature and high-pressure well according to claim 22, characterized in that: The third processing module is specifically used to calculate the drilling fluid density using the following formula: ρ=ρ0+[A1(T-T0)+A2(P-P0)] Among them, ρ is the drilling fluid density at the preset temperature value and the preset pressure value, ρ0 is the drilling fluid density at the test temperature value and the test pressure value, T0 is the test temperature value, P0 is the test pressure value, T is the preset temperature value, P is the preset pressure value, and A1 and A2 are drilling fluid characteristic parameters.

24. The apparatus for determining the wellbore pressure during tripping of a high-temperature and high-pressure well according to claim 22, wherein: The third processing module is specifically used to calculate the drilling fluid yield stress using the following formula: Among them, τ y is the yield stress of the drilling fluid at the preset temperature and pressure, τ y,0 is the yield stress of drilling fluid under test temperature and test pressure, and α1, β1, γ1, δ1, and ε1 are the characteristic parameters of drilling fluid.

25. The apparatus for determining the wellbore pressure during tripping of a high-temperature and high-pressure well according to claim 22, wherein: The third processing module is specifically used to calculate the plastic viscosity of the drilling fluid at a preset temperature value and a preset pressure value according to the preset test temperature value, the test pressure value, and the plastic viscosity of the drilling fluid at the test temperature value and the test pressure value using the following formula: Among them, μ p is the plastic viscosity of the drilling fluid at the preset temperature and pressure, μ p,0 is the plastic viscosity of the drilling fluid at the test temperature and test pressure values, and α2, β2, γ2, δ2, and ε2 are the characteristic parameters of the drilling fluid.

26. The apparatus for determining the wellbore pressure during tripping of a high-temperature and high-pressure well according to claim 22, wherein: The third processing module is specifically used to calculate the specific heat capacity of the drilling fluid at a preset temperature value and a preset pressure value according to the preset temperature value using the following formula: c fl =a+bT+cT -2 Among them, c fl is the specific heat capacity of the drilling fluid at the preset temperature value T, and a, b, c are the characteristic parameters of the drilling fluid.

27. The apparatus for determining the wellbore pressure during tripping of a high-temperature and high-pressure well as claimed in claim 22, characterized in that: The third processing module is specifically used to calculate the formation thermal conductivity at the preset temperature value according to the preset temperature value using the following formula: k e =k e,0 (1+βT) Among them, k e is the thermal conductivity of the formation at the preset temperature T, k e,0 is the thermal conductivity of the formation at 0℃, and β3 is the formation characteristic parameter.

28. A computer device comprising a memory, a processor and a computer program stored in the memory and executable on the processor, characterized in that: When the processor executes the computer program, the method according to any one of claims 1 to 16 is implemented.

29. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the method according to any one of claims 1 to 16 is implemented.

30. A computer program product, characterized in that The computer program product comprises a computer program, and when the computer program is executed by a processor, the method according to any one of claims 1 to 16 is implemented.

Images