CN114645703B - A computing system for simulating wellbore engineering construction - Google Patents

Abstract

The invention provides a computing system for simulating shaft engineering construction, which comprises a drilling basic data storage module, a data processing module and a data processing module, wherein the drilling basic data storage module is used for storing real-time drilling basic data; the well cleaning analysis calculation module is used for calculating the well cleaning degree of the well sections with different inclinations; the well-entering pipe column mechanical analysis and calculation module is used for calculating the pipe column mechanical; the well track middle target calculation module is used for middle target analysis and calculation; the well cementation construction calculation module is used for analyzing and calculating well cementation parameters; the stuck drill accident analysis and calculation module is used for analyzing and calculating stuck drill parameters; and the overflow well-killing construction calculation module is used for analyzing and calculating well-killing construction parameters. The invention realizes the timely calculation and analysis of the related parameters of the drilling engineering, has important significance for meeting the increasing demands of the drilling engineering, and plays an increasingly important role in improving the drilling benefit, enhancing the safety of the drilling operation and the like.

Description

A computing system for simulating wellbore engineering construction

Technical Field

The present invention relates to the technical field of oil and gas drilling and production engineering, and in particular to a computing system for simulating wellbore engineering construction.

Background Art

With the growing demand for oil in modern society, oil exploration and development has shown global characteristics. Drilling construction areas are spread all over the world and are widely distributed in complex surface and underground geological conditions such as oceans, deserts, swamps, hills, and mountains, which increases the complexity of drilling construction and the uncertainty of the drilling construction process, and puts forward higher requirements for drilling engineering design, risk analysis and control technology. The drilling engineering data analysis of complex wells, drilling risk analysis and control, and real-time optimization of drilling plans all require technical support from drilling engineering software. Drilling engineering software plays an increasingly important role in improving drilling efficiency and enhancing the safety of drilling operations.

For example, a patent document entitled "A drilling multi-parameter calculation system and method for drilling simulation" and published on July 20, 2021 with publication number CN113139704A records a drilling multi-parameter calculation system, including: a main server, which is used to obtain simulation data required for simulation calculations, adjust the calculation order of different drilling parameters, and collect and display the calculation results of each category of drilling parameters; multiple computing networks for corresponding category parameters, and the computing networks have: a first-level server, which is used to evaluate the amount of calculation required to calculate the current category parameters after receiving such calculation instructions, allocate a corresponding number of second-level servers, and divide the received corresponding category original data into tasks, and forward the received calculation results to the main server; multiple second-level servers, which are used to calculate according to the divided original data to obtain corresponding calculation results.

Summary of the invention

The purpose of the present invention is to solve at least one of the above-mentioned deficiencies in the prior art. For example, one of the purposes of the present invention is to provide a computing system for simulating wellbore engineering construction to solve the problem that the basic drilling data is complicated and cumbersome in the actual drilling process, and the drilling data cannot be quickly analyzed, so as to improve the real-time computing efficiency of drilling engineering design, risk analysis and control technology.

In order to achieve the above-mentioned purpose, the present invention provides a computing system for simulating wellbore engineering construction, the computing system includes a basic drilling data storage module, a wellbore cleaning analysis and calculation module, a wellbore string mechanical analysis and calculation module, a wellbore trajectory target calculation module, a cementing construction calculation module, a drill stuck accident analysis and calculation module and an overflow well pressure construction calculation module, wherein the basic drilling data storage module is configured to acquire and store real-time basic drilling data of the target well during the drilling process; the wellbore cleaning analysis and calculation module is connected to the basic drilling data storage module and is configured to perform wellbore cleanliness calculation for well sections with different inclinations; the wellbore string mechanical analysis and calculation module is connected to the basic drilling data storage module and is configured to perform wellbore cleanliness calculation for well sections with different inclinations. The pipe string during the operation is subjected to stress analysis calculation, deformation analysis calculation and safety evaluation calculation; the wellbore trajectory target calculation module is connected to the drilling basic data storage module, and is configured to be able to perform target analysis calculation in combination with directional well trajectory monitoring data and wellbore trajectory three-dimensional graphics; the cementing construction calculation module is connected to the drilling basic data storage module, and is configured to be able to analyze and calculate cementing parameters in combination with real-time drilling basic data; the stuck drill accident analysis calculation module is connected to the drilling basic data storage module, and is configured to be able to analyze and calculate stuck drill parameters in combination with real-time drilling basic data; the overflow well killing construction calculation module is connected to the drilling basic data storage module, and is configured to be able to analyze and calculate well killing construction parameters in combination with real-time drilling basic data.

In an exemplary embodiment of the present invention, the wellbore cleaning analysis and calculation module may include a screening submodule, a first input data extraction submodule, a wellbore cleanliness calculation submodule and a wellbore cleanliness determination submodule, wherein the screening submodule is configured to automatically screen and determine the clean area to which the well section to be calculated belongs according to the well inclination angle of the well section to be calculated, and output the partition result; the first input data extraction submodule is respectively connected to the drilling basic data storage module and the screening submodule, and is configured to extract the first input data required for the calculation of the cleanliness of the well section to be calculated from the drilling basic data storage module based on the partition result; the wellbore cleanliness calculation submodule includes a selection unit, a first clean area calculation unit, a second clean area calculation unit and a third clean area calculation unit, wherein the selection unit is connected to the screening submodule and is configured to control the first clean area calculation unit according to the partition result of the well section to be calculated. One of the first cleaning zone calculation unit, the second cleaning zone calculation unit and the third cleaning zone calculation unit is connected to the first input data extraction submodule, the first cleaning zone calculation unit is configured with a wellbore cleaning calculation model for a small-angle well section, and can calculate and output a first calculation result, the first calculation result includes a drill cuttings particle transmission ratio, the second cleaning zone calculation unit is configured with a wellbore cleaning calculation model for a medium-angle well section, and can calculate and output a second calculation result, the second calculation result includes a ratio of the thickness of the cuttings bed to the wellbore diameter and an annulus stop return speed, the third cleaning zone calculation unit is configured with a wellbore cleaning calculation model for a large-angle well section, and can calculate and output a third calculation result, the third calculation result includes a ratio of the thickness of the cuttings bed to the wellbore diameter; the wellbore cleanliness determination submodule is connected to the wellbore cleanliness calculation submodule, and is configured to analyze and determine the cleaning effect of the well section to be calculated based on the calculation result output by the wellbore condition calculation module.

In an exemplary embodiment of the present invention, the wellbore cleaning analysis and calculation module may also include a wellbore cleaning analysis submodule, which is connected to the wellbore cleanliness calculation submodule and is configured to output the wellbore cleaning analysis curve of the target well to judge the wellbore cleanliness, and includes a small-inclination cuttings transmission ratio curve drawing unit, a medium-inclination critical annulus return velocity and annulus return velocity curve drawing unit, a medium-inclination cuttings bed relative thickness curve drawing unit, a high-inclination critical annulus return velocity and annulus return velocity curve drawing unit, a critical annulus displacement curve drawing unit, and a high-inclination cuttings bed relative thickness curve drawing unit.

In an exemplary embodiment of the present invention, the calculation model for wellbore cleaning in a low-angle well section may be as shown in equations (1) to (3):

In formulas (1) to (3), V _sx is the settling velocity of drill cuttings, m/s; d _s is the equivalent diameter of drill cuttings, cm; ρ _s is the density of drill cuttings, g/cm ³ ; ρ _m is the density of drilling fluid, g/cm ³ ; is the cuttings particle shape coefficient, dimensionless; _Va is the drilling fluid annulus return velocity, m/s; _Qa is the drilling fluid flow rate, L/s; _Dh is the wellbore diameter, mm; _Dp is the drill pipe outer diameter, mm; Rt is the cuttings particle transmission ratio, dimensionless.

The calculation model for wellbore cleaning in the medium-angle well section can be shown as formulas (4) to (7):

h _zx ＝0.015D _h (AV+6.15AV ^0.5 )(1+0.587E)(V _lzx -V _a ) Equation (5)

Hzx＝(hzx/Dh)×100% Formula (6)

In formulas (4) to (7), V _lzx is the critical annular return velocity in the medium-inclined well section, m/s; AV is the apparent viscosity of the drilling fluid, mPa.s; θ is the well inclination angle, degree; h _zx is the thickness of the cuttings bed in the medium-inclined well section, mm; H _zx is the relative thickness of the cuttings bed in the medium-inclined well section, dimensionless; V _p is the annular stop return velocity, m/s; A _bed is the cross-sectional area of the cuttings bed perpendicular to the wellbore axis, mm ² ; C is the cuttings concentration in the cuttings bed, dimensionless; g is the gravitational acceleration, m/s ² ; E is the eccentricity of the drill string, dimensionless; L is the cross-sectional width of the cuttings bed perpendicular to the wellbore axis, mm; η is the friction coefficient between the cuttings bed and the lower wellbore wall, generally taken as 0.2; PV is the plastic viscosity, mPa; YP is the dynamic shear force, Pa; K is the consistency coefficient, Pa·s ⁿ ; n is the fluidity index, dimensionless.

The calculation model for wellbore cleaning in a highly deviated well section can be shown in equations (15) to (20):

H _dx = (h _dx / _D h) × 100% Formula (20)

In formulas (15) to (20), V _Ldx is the critical annular return velocity in the highly deviated well section, m/s; V _sd is the cuttings settling velocity in the highly deviated well section, m/s; V _jx is the mechanical penetration rate, m/h; _Cang is the well inclination correction coefficient, dimensionless; C _size is the cuttings size correction coefficient, dimensionless; C _denF is the drilling fluid density correction coefficient, dimensionless; C _rpm is the drill string speed correction coefficient, dimensionless; Q _Ldx is the critical annular displacement without cuttings bed, L/s; A' _bed is the cuttings bed area, mm ² ; h _dx is the cuttings bed thickness in the highly deviated well section, mm; Q _a is the drilling fluid displacement, L/s; H _dx is the relative thickness of the cuttings bed in the highly deviated well section, dimensionless.

In an exemplary embodiment of the present invention, the wellbore string mechanical analysis and calculation module may include a second input data extraction submodule, a casing force analysis and calculation submodule and a casing strength verification submodule, wherein the second input data extraction submodule is connected to the drilling basic data storage module, and is configured to be able to extract the second input data required for the wellbore string mechanical calculation from the drilling basic data storage module; the casing force analysis and calculation submodule is connected to the second input data extraction submodule, and is configured with a drilling tool force calculation model, which can calculate and output the axial tension, tensile strength, external extrusion pressure, extrusion strength, internal pressure and internal pressure resistance of the casing string in different production periods; the casing strength verification submodule is connected to the second input data extraction submodule, and is configured with a casing strength verification algorithm, which can verify the casing string design data.

In an exemplary embodiment of the present invention, the casing strength verification algorithm may include the following steps: (1) obtaining original design data of the casing string; (2) calculating the effective external extrusion pressure p _ce1 at the casing shoe according to the original design data; (3) selecting the i-th section of casing to be verified according to the load and geometric constraints, 1≤i≤m, when i＝1, the section of casing is located at the bottom of the well, and when i＝m, the section of casing is located at the wellhead; (4) calculating the running depth L _i of the i-th section of casing ; (5) Determine whether the internal pressure resistance and tensile strength of the i-th section of casing meet the verification requirements. If the internal pressure resistance and tensile strength of the i-th section of casing meet the verification requirements, proceed to step (6); otherwise, return to step (2) and select a casing with a higher steel grade or wall thickness to redesign the i-th section; (6) Determine whether the i-th section of casing has reached the wellhead. If the i-th section of casing has not reached the wellhead, return to step (3), set i=i+1, calculate the insertion depth of the next section of casing, and verify the internal pressure resistance and tensile strength; if the i-th section of casing has reached the wellhead, the calculation is terminated and the design results of all casing strings are output.

In an exemplary embodiment of the present invention, the wellbore string mechanical analysis calculation module may also include a mechanical analysis result output submodule, which is connected to the casing stress condition calculation submodule and the casing strength verification submodule, and is configured to graphically output the casing string stress condition and strength verification results.

In an exemplary embodiment of the present invention, the borehole trajectory target calculation module may include a third input data extraction submodule, an anti-collision analysis submodule, a borehole trajectory generation submodule, a borehole trajectory prediction submodule and a borehole trajectory data output submodule, wherein the third input data extraction submodule is connected to the drilling basic data storage module and is configured to be able to extract the third input data required for the borehole trajectory target calculation from the drilling basic data storage module; the anti-collision analysis submodule is connected to the third input data extraction submodule and is configured to be able to calculate the borehole trajectory distance between the target well and the adjacent well; the borehole trajectory generation submodule is connected to the third input data extraction submodule and is configured to be able to automatically generate a three-dimensional borehole trajectory map according to the borehole trajectory data; the borehole trajectory prediction submodule is connected to the borehole trajectory generation submodule and is configured to be able to automatically generate a three-dimensional borehole trajectory prediction map when the well depth, well inclination and azimuth of the expected borehole trajectory are input; the borehole trajectory data output submodule is connected to the borehole trajectory generation submodule and is configured to be able to output the borehole trajectory data.

In an exemplary embodiment of the present invention, the wellbore trajectory target calculation module may also include a trajectory playback submodule, a real-time trajectory submodule and a two-dimensional trajectory map drawing submodule, wherein the trajectory playback submodule is connected to the wellbore trajectory generation submodule, and is configured to automatically simulate and playback a demonstration animation of the wellbore drilling process according to the three-dimensional wellbore trajectory map; the real-time trajectory submodule is connected to the wellbore trajectory generation submodule, and is configured to generate a demonstration animation of the wellbore drilling process in real time according to the three-dimensional wellbore trajectory map; the two-dimensional trajectory map drawing submodule is connected to the wellbore trajectory generation submodule, and is configured to automatically draw the wellbore trajectory as a vertical projection map and a horizontal projection map.

In an exemplary embodiment of the present invention, the cementing construction calculation module may include a fourth input data extraction submodule, a cementing construction parameter calculation submodule and a cementing construction dynamic simulation submodule, wherein the fourth input data extraction submodule is connected to the drilling basic data storage module and is configured to be able to extract the fourth input data required for cementing construction calculation from the drilling basic data storage module; the cementing construction parameter calculation submodule is connected to the fourth input data extraction submodule and is configured with a cementing construction calculation model, which can calculate cementing construction parameters; the cementing construction dynamic simulation submodule is connected to the cementing construction parameter calculation submodule and is configured with a cementing construction numerical model, which can visualize and simulate the cementing construction operation process.

In an exemplary embodiment of the present invention, the drill stuck accident analysis and calculation module may include a fifth input data extraction submodule, a torsion circle number calculation submodule, a stuck point location calculation submodule, a jamming agent calculation submodule, a pressure reduction calculation submodule and a falling time calculation submodule, wherein the fifth input data extraction submodule is connected to the drilling basic data storage module and is configured to be able to extract the fifth input data required for the drill stuck accident analysis and calculation from the drilling basic data storage module; the torsion circle number calculation submodule is connected to the fifth input data extraction submodule and is configured with a torsion circle number calculation model and is able to calculate the torsion circle number; the stuck point location ... The extraction submodule is connected to the fifth input data extraction submodule and is configured with a stuck point location calculation model, which can calculate the stuck point location according to the well type and drilling tool type of the target well; the jamming agent calculation submodule is connected to the fifth input data extraction submodule and is configured with a jamming agent dosage calculation model and a maximum pump pressure calculation model, which can respectively calculate and output the total amount of jamming agent and the maximum pump pressure of the jamming agent; the pressure reduction calculation submodule is connected to the fifth input data extraction submodule and is configured with a U-shaped tube effect pressure reduction calculation model, which can calculate and output the pressure reduction parameters; the falling time calculation submodule is connected to the fifth input data extraction submodule and is configured with a falling time calculation model, which can calculate and output the falling time of the ball body of the ball.

In an exemplary embodiment of the present invention, the overflow well-killing construction calculation module may include a sixth input data extraction submodule, a well-killing construction calculation submodule, a well-killing construction order generation submodule, a well-killing construction dynamic simulation submodule and a well-killing operation drawing submodule, wherein the sixth input data extraction submodule is connected to the drilling basic data storage module and is configured to be able to extract the sixth input data required for the overflow well-killing construction calculation from the drilling basic data storage module; the well-killing construction calculation submodule is connected to the sixth input data extraction submodule and is configured with a well-killing construction calculation model, which can calculate the well-killing construction parameters; the well-killing construction order generation submodule is connected to the well-killing construction calculation submodule and is configured to be able to automatically generate a well-killing construction form according to the calculation results of the well-killing construction parameters; the well-killing construction dynamic simulation submodule is connected to the well-killing construction calculation submodule and is configured with a well-killing construction numerical model, which can visualize and simulate the wellbore well-killing operation process; the well-killing operation drawing submodule is connected to the well-killing construction dynamic simulation submodule and is configured to be able to graphically output the comparison results of real-time well-killing data and theoretical well-killing data.

In an exemplary embodiment of the present invention, the computing system may further include a pressure consumption calculation module, the pressure consumption calculation model is connected to the drilling basic data storage module, and includes a seventh input data extraction submodule, a pressure consumption calculation submodule, a pressure consumption prediction submodule, an actual pressure consumption curve drawing submodule and a pressure consumption prediction curve drawing submodule, wherein the seventh input data extraction submodule is connected to the drilling basic data storage module and is configured to be able to extract the seventh input data required for the pressure consumption calculation from the drilling basic data storage module; the pressure consumption calculation submodule is connected to the seventh input data extraction submodule and is configured to be able to The submodule is capable of automatically calculating and outputting the inner-tube pressure loss, outer-tube pressure loss and total circulating pressure loss of the current downhole drilling tool combination; the actual pressure loss curve drawing submodule is connected to the pressure loss calculation submodule, and is configured to automatically draw the calculation result of the pressure loss calculation submodule as an actual pressure loss curve; the pressure loss prediction submodule is connected to the seventh input data extraction submodule, and is configured to automatically calculate and output the circulating pressure loss prediction result of the section to be drilled under different working conditions; the pressure loss prediction curve drawing submodule is connected to the pressure loss prediction submodule, and is configured to automatically draw the calculation result of the pressure loss prediction curve drawing submodule as a pressure loss prediction curve.

Compared with the prior art, the beneficial effects of the present invention include at least one of the following:

(1) The present invention can quickly analyze and calculate drilling parameters based on basic drilling data in real time, which greatly shortens the calculation time and avoids the errors of manual calculation;

(2) The present invention can utilize the real-time data and dynamic data of the drilling site to perform real-time calculation and analysis on the parameters related to the drilling project, thereby providing data support for safe drilling. At the same time, the drilling project parameters can be optimized based on the analysis results, thereby providing real-time guidance for drilling production.

(3) The present invention realizes timely calculation and analysis of drilling engineering related parameters, which is of great significance to meet the ever-increasing demand for drilling engineering and plays an increasingly important role in improving drilling efficiency and enhancing drilling operation safety;

(4) The present invention is the first in China to realize the integrated application of multiple analysis and calculation software for drilling engineering on the same engineering technology information platform, which can provide unified database support, remote communication support, numerical calculation support, graphic visualization support, etc.

BRIEF DESCRIPTION OF THE DRAWINGS

The above and other objects and/or features of the present invention will become more apparent through the following description in conjunction with the accompanying drawings, in which:

FIG1 shows a schematic diagram of the structure of a computing system for simulating wellbore engineering construction according to an exemplary embodiment of the present invention.

Description of reference numerals:

100-drilling basic data storage module, 200-wellbore cleaning analysis and calculation module, 210-sieve submodule, 220-first input data extraction submodule, 230-wellbore cleanliness calculation submodule, 231-selection unit, 232-first clean area calculation unit, 233-second clean area calculation unit, 234-third clean area calculation unit, 240-wellbore cleanliness determination submodule, 250-wellbore cleaning analysis submodule, 300-wellbore string mechanical analysis and calculation module, 310-second input data extraction submodule, 320-casing force analysis and calculation submodule, 330-casing strength verification submodule, 340-mechanical analysis result output submodule, 400-wellbore trajectory target calculation module, 410-third input data extraction submodule, 420-anti-collision analysis submodule, 430-wellbore trajectory generation submodule, 440-wellbore trajectory prediction submodule, 4 50-wellbore trajectory data output submodule, 460-trajectory playback submodule, 470-real-time trajectory submodule, 480-two-dimensional trajectory map drawing submodule, 500-cementing construction calculation module, 510-fourth input data extraction submodule, 520-cementing construction parameter calculation submodule, 530-cementing construction dynamic simulation submodule, 600-drilling accident analysis calculation module, 610-fifth input data extraction submodule, 620-twist circle calculation submodule, 630-stuck point position calculation submodule, 640-de-stuck agent calculation submodule, 650-pressure reduction calculation submodule, 660-falling time calculation submodule, 700-overflow well killing construction calculation module, 710-sixth input data extraction submodule, 720-well killing construction calculation submodule, 730-well killing construction order generation submodule, 740-well killing construction dynamic simulation submodule, 750-well killing operation drawing submodule.

DETAILED DESCRIPTION

Hereinafter, the computing system for simulating wellbore engineering construction of the present invention will be described in detail with reference to exemplary embodiments and accompanying drawings.

It should be noted that “first”, “second”, “third”, “fourth”, etc. are merely for the convenience of description and distinction and should not be understood as indicating or implying relative importance.

In an exemplary embodiment of the present invention, as shown in FIG1 , a computing system for simulating wellbore engineering construction may include a basic drilling data storage module 100, a wellbore cleaning analysis and calculation module 200, a wellbore string mechanics analysis and calculation module 300, a wellbore trajectory target calculation module 400, a cementing construction calculation module 500, a stuck drill accident analysis and calculation module 600, and an overflow well pressure construction calculation module 700.

The basic drilling data storage module is configured to acquire and store real-time basic drilling data of the target well during the drilling process. Real-time basic drilling data includes wellbore trajectory data, wellbore diameter data, wellbore structure data, drilling tool assembly data, drilling fluid performance data, real-time drilling data, casing data, maximum number of twisting circles of the drilling tool allowed when the hanging casing is inverted, low pump stroke test data, wellhead pressure test data, overflow key data, and K value coefficient data.

The wellbore cleaning analysis and calculation module 200 is connected to the drilling basic data storage module 100 and is configured to be able to perform wellbore cleanliness calculation for well sections with different inclinations.

The mechanical analysis and calculation module 300 of the tubing entering the well is connected to the basic drilling data storage module 100, and is configured to perform force analysis and calculation, deformation analysis and calculation, and safety evaluation calculation on the tubing in different operation processes.

The wellbore trajectory target calculation module 400 is connected to the drilling basic data storage module 100, and is configured to be able to perform target analysis and calculation in combination with the directional well trajectory monitoring data and the wellbore trajectory three-dimensional graphics.

The cementing construction calculation module 500 is connected to the drilling basic data storage module 100 and is configured to be able to analyze and calculate cementing parameters in combination with real-time drilling basic data.

The stuck pipe accident analysis and calculation module 600 is connected to the drilling basic data storage module 100 and is configured to analyze and calculate stuck pipe parameters in combination with real-time drilling basic data.

The overflow well-killing construction calculation module 700 is connected to the drilling basic data storage module 100 and is configured to be able to analyze and calculate well-killing construction parameters in combination with real-time drilling basic data.

In this embodiment, the setting requirements of the wellbore cleaning analysis and calculation module are: combining the real-time well depth, wellbore structure, well inclination, drilling fluid parameters, drilling displacement, and cuttings diameter to analyze the wellbore cleanliness in real time, and referencing the safety critical value to optimize the drilling parameters.

Therefore, as shown in FIG. 1 , the borehole cleaning analysis calculation module 200 may include a screen submodule 210 , a first input data extraction submodule 220 , a borehole cleanliness calculation submodule 230 , and a borehole cleanliness determination submodule 240 .

Among them, the screening submodule 210 is configured to automatically screen and determine the clean area to which the well section to be calculated belongs according to the well inclination angle of the well section to be calculated, and output the partition result. For example, when the well inclination angle is 0-30°, the well section to be calculated is a low-inclination well section, belonging to the first clean area; when the well inclination angle is 30-60°, the well section to be calculated is a medium-inclination well section, belonging to the second clean area; when the well inclination angle is 60-90°, the well section to be calculated is a high-inclination well section, belonging to the third clean area.

The first input data extraction submodule 220 is connected to the drilling basic data storage module 100 and the sieve submodule 210 respectively, and is configured to extract the first input data required for the calculation of the cleaning condition of the well section to be calculated from the drilling basic data storage module based on the partition result. The first input data includes: wellbore trajectory data, wellbore data, wellbore structure data, drilling tool assembly data and drilling fluid performance data. In addition, the first input data extraction submodule can also obtain cuttings parameters through user input.

The wellbore cleanliness calculation submodule 230 may include a selection unit 231 , a first clean area calculation unit 232 , a second clean area calculation unit 233 , and a third clean area calculation unit 234 .

Among them, the selection unit 231 is connected to the screening submodule 210, and is configured to control one of the first clean area calculation unit 232, the second clean area calculation unit 233 and the third clean area calculation unit 234 to be connected to the first input data extraction submodule according to the partition result of the well section to be calculated.

When the partition result output by the sieve submodule 210 is the first clean zone, the selection unit 231 controls the first clean zone calculation unit 232 to connect with the first input data extraction submodule 220. The first clean zone calculation unit is configured with a wellbore cleaning calculation model for a small-angle well section, and can calculate and output a first calculation result, which includes a drill cutting particle transmission ratio. For the first clean zone, the area is a small-angle well section (0° to 30°). Since the well inclination angle is small, a cuttings bed is generally not formed, so the transmission ratio of the drill cutting particles is used as a pre-judgment criterion for wellbore cleaning.

The drilling cuttings particle transmission ratio Rt can be calculated using the wellbore cleaning calculation model for a low-angle well section. If Rt≥0.5, it is determined that the cleanliness of the first cleaning area meets the cleaning effect; if Rt<0.5, it is determined that the cleanliness of the first cleaning area does not meet the cleaning effect, and the drilling fluid performance and/or drilling fluid displacement need to be adjusted. The wellbore cleaning calculation model for a low-angle well section can be shown in the following equations (1) to (3):

In formulas (1) to (3), V _sx is the settling velocity of drill cuttings, m/s; d _s is the equivalent diameter of drill cuttings, cm; ρ _s is the density of drill cuttings, g/cm ³ ; ρ _m is the density of drilling fluid, g/cm ³ ; is the cuttings particle shape coefficient, dimensionless; _Va is the drilling fluid annulus return velocity, m/s; Q is the drilling fluid flow rate, L/s; _Dh is the wellbore diameter, mm; _Dp is the drill pipe outer diameter, mm; Rt is the cuttings particle transmission ratio, dimensionless.

The equivalent diameter _ds of the drill cuttings particles may be 0.5 to 3.0 cm. For example, the equivalent diameters of the drill cuttings particles may be 0.5 cm, 1.0 cm, 1.5 cm, 2.0 cm, 2.5 cm, and 3.0 cm, respectively, for calculation.

Cuttings particle shape factor It can be 0.5 to 1.0. For example, when the drill cuttings particle shape is spherical, the drill cuttings particle shape coefficient Can be 1.0; when the drill cuttings particle shape is ellipsoidal, the drill cuttings particle shape coefficient It can be 0.9; when the shape of the drill cuttings particles is polygonal, the shape coefficient of the drill cuttings particles is It can be 0.8; when the drill cuttings particle shape is long strip, the drill cuttings particle shape coefficient It can be 0.6; when the drill cuttings particle shape is flaky, the drill cuttings particle shape coefficient It can be 0.5.

The borehole diameter D _h can be obtained through caliper logging. For the well section where caliper logging has not been performed, the borehole diameter can be calculated based on the drill bit diameter and an expansion factor of 5% to 10%.

When the partition result output by the sieve submodule 210 is the second clean zone, the selection unit 231 controls the second clean zone calculation unit 233 to connect with the first input data extraction submodule 220. The second clean zone calculation unit is configured with a wellbore cleaning calculation model for a medium-inclination well section, and is capable of calculating and outputting a second calculation result, which includes the ratio of the cuttings bed thickness to the wellbore diameter and the annular stop return speed. For the second clean zone, this area is a medium-inclination well section (30° to 60°), which is prone to forming a cuttings bed. The ratio of the cuttings bed thickness to the wellbore diameter and the annular stop return speed can be used as a pre-judgment criterion for wellbore cleaning.

The calculation model for borehole cleaning in a medium-inclined well section can be used to calculate the relative thickness of the cuttings bed in the medium-inclined well section (that is, the ratio of the cuttings bed thickness to the borehole diameter) H _zx and the annular stop return speed V _p . If H _zx ≤10% and V _p ≥V _a , it is determined that the cleanliness of the second cleaning area meets the cleaning effect; if H _zx >10% or V _p <V _a , it is determined that the cleanliness of the second cleaning area does not meet the cleaning effect, and the drilling fluid performance and/or drilling fluid displacement need to be adjusted. The calculation model for borehole cleaning in a medium-inclined well section can be shown in equations (4) to (7):

h _zx ＝0.015D _h (AV+6.15AV ^0.5 )(1+0.587E)(V _lzx -V _a ) (5)

H _zx = (h _zx /D _h )×100% (6)

In formulas (4) to (7), V _lzx is the critical annular return velocity in the medium-inclined well section, m/s; AV is the apparent viscosity of the drilling fluid, mPa.s; θ is the well inclination angle, degrees; h _zx is the thickness of the cuttings bed in the medium-inclined well section, mm; H _zx is the relative thickness of the cuttings bed in the medium-inclined well section, dimensionless; V _p is the annular stop return velocity, m/s; A _bed is the cross-sectional area of the cuttings bed perpendicular to the wellbore axis, mm ² ; C is the cuttings concentration in the cuttings bed, dimensionless; g is the gravitational acceleration, m/s ² ; E is the eccentricity of the drill string, dimensionless, and the default value is 2100000; L is the cross-sectional width of the cuttings bed perpendicular to the wellbore axis, mm; η is the friction coefficient between the cuttings bed and the lower wellbore wall, generally taken as 0.2; PV is the plastic viscosity, mPa; YP is the dynamic shear force, Pa; K is the consistency coefficient, Pa·s ⁿ ; n is the fluidity index, dimensionless.

Among them, the relevant parameters in the wellbore cleaning calculation model for the medium-angle well section can be calculated using the following equations (8) to (14).

In formulas (8) to (14), L is the cross-sectional width of the cuttings bed perpendicular to the wellbore axis, mm; D _h is the wellbore diameter, mm; h _zx is the thickness of the cuttings bed in the medium-inclination well section, mm; A _bed is the cross-sectional area of the cuttings bed perpendicular to the wellbore axis, mm ² ; AV is the apparent viscosity of the drilling fluid, mPa.s; is the reading of the rotational viscometer at 600 rpm; PV is the plastic viscosity, mPa·s; is the reading of the rotational viscometer at 300 revolutions; YP is the dynamic shear force, Pa; K is the consistency coefficient, Pa·s ⁿ ; n is the fluidity index, dimensionless.

When the partition result output by the sieve submodule 210 is the third clean zone, the selection unit 231 controls the third clean zone calculation unit 234 to connect with the first input data extraction submodule 220. The third clean zone calculation unit is configured with a high-angle well section wellbore cleaning calculation model, which can calculate and output the third calculation result, and the third calculation result includes the ratio of the cuttings bed thickness to the wellbore diameter. For the third clean zone, this area is a high-angle well section (60°~90°), which is a stable area of the cuttings bed. The percentage of the wellbore diameter occupied by the cuttings bed can be used as a pre-judgment criterion for wellbore cleaning.

The calculation model for the wellbore cleaning of the highly deviated well section can be used to calculate the relative thickness H _dx of the cuttings bed in the highly deviated well section. If H _dx ≤10%, it is determined that the cleanliness of the third cleaning area meets the cleaning effect; if H _dx >10%, it is determined that the cleanliness of the third cleaning area does not meet the cleaning effect, and the drilling fluid performance and/or drilling fluid displacement need to be adjusted. The calculation model for the wellbore cleaning of the highly deviated well section is shown in equations (15) to (20):

H _dx = (h _dx /D _h )×100% (20)

In formulas (15) to (20), V _Ldx is the critical annular return velocity in the highly deviated well section, m/s; V _sd is the cuttings settling velocity in the highly deviated well section, m/s; V _jx is the mechanical penetration rate, m/h; _Cang is the well inclination correction coefficient, dimensionless; C _size is the cuttings size correction coefficient, dimensionless; C _denF is the drilling fluid density correction coefficient, dimensionless; C _rpm is the drill string speed correction coefficient, dimensionless; Q _Ldx is the critical annular displacement without cuttings bed, L/s; A' _bed is the cuttings bed area, mm ² ; h _dx is the cuttings bed thickness in the highly deviated well section, mm; Q _a is the drilling fluid displacement, L/s; H _dx is the relative thickness of the cuttings bed in the highly deviated well section, dimensionless.

Among them, the relevant parameters in the wellbore cleaning calculation model for high-angle well sections can be calculated using the following equations (21) to (24).

C _ang ＝0.0342θ-0.000233θ ² -0.213 (21)

C _size =1.286-0.04094448d _s (22)

In formulas (21) to (24), _Cang is the correction coefficient for well inclination, dimensionless; _Csize is the correction coefficient for drill cuttings size, dimensionless; _CdenF is the correction coefficient for drilling fluid density, dimensionless; _Crpm is the correction coefficient for drill string speed, dimensionless; θ is the well inclination, degree; _ds is the equivalent diameter of drill cuttings particles, cm; _ρm is the drilling fluid density, g/ ^cm3 ; and N is the drill string speed, r/min.

The wellbore cleanliness determination submodule 240 is connected to the wellbore cleanliness calculation submodule 230 and is configured to analyze and determine the cleaning effect of the well section to be calculated according to the calculation result output by the wellbore condition calculation module.

In addition, the borehole cleaning analysis calculation module 200 may further include a borehole cleaning analysis submodule 250. The borehole cleaning analysis submodule 250 is connected to the borehole cleanliness calculation submodule 230, and is configured to output a borehole cleaning analysis curve of the target well to determine the borehole cleanliness, and includes a small-inclination cuttings transmission ratio curve drawing unit, a medium-inclination critical annular return velocity and annular return velocity curve drawing unit, a medium-inclination cuttings bed relative thickness curve drawing unit, a high-inclination critical annular return velocity and annular return velocity curve drawing unit, a critical annular displacement curve drawing unit, and a high-inclination cuttings bed relative thickness curve drawing unit.

In this embodiment, the setting requirement of the mechanical analysis and calculation unit of the wellbore string is to analyze the stress and deformation of the drilling string during various working conditions and evaluate the safety. Therefore, as shown in FIG1 , the mechanical analysis and calculation module 300 of the wellbore string may include a second input data extraction submodule 310, a casing force analysis and calculation submodule 320, and a casing strength verification submodule 330.

The second input data extraction submodule 310 is connected to the drilling basic data storage module 100 and is configured to extract the second input data required for the mechanical calculation of the wellbore string from the drilling basic data storage module. The second input data includes: wellbore structure data, drilling tool assembly data, real-time drilling data, drilling fluid performance data and wellbore trajectory data. In addition, the second input data extraction submodule can obtain the drilling tool steel grade and strength data by manual input.

The casing force analysis calculation submodule 320 is connected to the second input data extraction submodule 310 and is configured with a drilling tool force calculation model, which can calculate and output the axial tension, tensile strength, external extrusion pressure, collapse strength, internal pressure and internal pressure resistance of the casing string in different production periods.

It should be noted that the working conditions of the drill string in the well vary with the drilling method (rotary table or downhole power drilling) and the drilling process (such as normal drilling, tripping, etc.). Under different working conditions, the drill string has different working states and is subjected to different forces. These include: tensile stress caused by axial tension, shear stress caused by torque, bending stress caused by drill string bending, additional stress caused by drill string buckling, tensile stress caused by drilling fluid pressure, etc. Analysis shows that the direction of drilling pressure is opposite to that of tension, and the direction of friction is related to the trajectory of drilling tool movement. The dynamic load value is relatively small. For the convenience of analysis, only the axial tension caused by the floating weight of the drilling tool is considered here.

Therefore, the drilling tool force calculation model is shown in equations (25) to (37).

① The axial force generated by the buoyant weight of the drill string.

Where, _Fm is the axial force generated by the buoyant weight, N; _Kb is the buoyancy coefficient of the drilling fluid, dimensionless; _Ldp is the drill pipe length, m; _qdp is the weight of the drill pipe per unit length in the air, kg/m; _Lc is the drill collar length, m; _qc is the weight of the drill collar per unit length in the air, kg/m; _ρs is the density of the drill cuttings particles, g/ ^cm3 ; _ρm is the density of the drilling fluid outside the pipe, g/ ^cm3 .

② External extrusion pressure on the drill string.

P _oc =Hγ _m -(HL)γ _f (27)

γ _m =ρ _m g (28)

γ _f =ρ _f g (29)

In the formula, P _oc is the external extrusion pressure at the calculation point, Pa; H is the depth of the calculation point, m; γ _m is the density of the drilling fluid outside the pipe, N/m ³ ; γ _f is the density of the drilling fluid inside the pipe, N/m ³ ; L is the distance between the liquid in the pipe and the wellhead, m; ρ _m is the density of the drilling fluid outside the pipe, kg/m ³ ; ρ _f is the density of the drilling fluid inside the pipe, kg/m ³ .

③Tensile strength of drill string.

As far as drill pipe is concerned, usually the tensile strength of the joint is much greater than the tensile strength of the drill pipe body, so the tensile strength of the drill pipe can usually be based on the tensile strength of the drill pipe body.

The tensile strength of the drill pipe body is:

P＝ _YmA / ¹⁰³ (30)

Where, P is the minimum tensile strength, KN; Y _m is the drill string yield strength, MPa; A is the cross-sectional area of the drill string body, mm ² .

④ Drill pipe anti-collapse strength.

The minimum collapse strength of the pipe body without axial load and internal pressure is determined by the specified minimum yield pressure of the material and the geometric dimensions of the pipe cross section.

The casing strength verification submodule 330 is connected to the second input data extraction submodule 310 and is configured with a casing strength verification algorithm, and is capable of verifying casing string design data.

It should be noted that the method of casing strength verification is: first design from bottom to top according to the collapse strength, and verify the tensile strength and internal pressure strength at the same time. When the tensile strength or internal pressure strength does not meet the requirements, select a casing one level higher than the previous section, change to tensile strength or internal pressure strength design, and verify the collapse strength until the design requirements are met.

Therefore, the casing strength verification algorithm can include the following steps:

(1) Obtain the original data of casing string design. The original data of casing string design may include steel grade, wall thickness, thread type, segment weight, segment length, tensile coefficient, anti-extrusion coefficient and compressive coefficient.

(2) Calculate the effective external extrusion pressure p _ce1 at the casing shoe (that is, deep at the bottom of the casing) based on the original design data.

Among them, the calculation formula for the effective external extrusion pressure p _ce1 is as follows.

a) Vertical Well

①Surface casing and technical casing

For non-plastic creep formations: P _ce = 0.00981[ρ _m -(1-k _m )ρ _min ]H (31)

For plastic creep formations:

② Production of casing and tail pipe

For non-plastic creep formations: P _ce = 0.00981[ρ _m -(1-k _m )ρ _w ]H (33)

For plastic creep formations:

Wherein, P _ce is the effective external pressure, MPa; ρ _min is the minimum drilling fluid density for the next drilling, g/cm ³ ; ρ _m is the drilling fluid density during cementing, g/cm ³ ; ρ _w is the completion fluid density, g/cm ³ ; υ is the Poisson's coefficient of the formation rock, which is 0.3~0.5 and dimensionless; G _v is the pressure gradient of the overlying rock formation, MPa/m; _km is the hollowing coefficient, which is 0~1, 1 means full hollowing; H is the well depth of the calculation point, m.

b) Directional Well

The effective external stress of a directional well should be calculated by converting the measured well depth of the tortuosity and inclined straight section into the vertical well depth.

(3) The i-th section of casing to be checked is selected according to the load and geometric constraints, 1≤i≤m. When i=1, the section of casing is located at the bottom of the well, and when i=m, the section of casing is located at the wellhead.

(4) Calculate the insertion depth of the i-th section of casing.

It should be noted that the running depth H refers to the well depth, and the running length refers to the difference between two well depths, for example, L ₁ =H ₁ -H ² .

For example, the second section of casing is then selected, and the lowerable depth H ₂ of the second section of casing is calculated, thereby determining the lowering length L ₁ of the first section of casing. Since the lowering length L ₁ of the first section of casing depends on the lowerable depth H ₂ of the second section of casing, the second section of casing should be selected with a lower anti-external collapse strength than the first section of casing.

The insertion depth _H2 of the second section of casing is determined by the following formula.

in, b＝C ₁ C ₂ +2C ₂ C ₃ ;

The running length _L1 of the first section of casing is: _L1 = _H1 - _H2 .

When n>3, the insertion depth _H2 of the second section of casing is determined by the following formula.

in, b＝C ₁ C ₂ +2C ₂ C ₃ ;

Wherein, _H2 is the running depth of the second section of casing, m; _Gce is the effective external pressure gradient of the casing, MPa/m; Sc is the specified anti-collapse coefficient; _Pcon is the anti-collapse strength of the nth section of casing, MPa; _Tyn is the yield strength of the nth section of casing, KN; _kf is the buoyancy coefficient; _qi is the mass per unit length of the i-th section of casing below the design section, Kg/m; _Hi is the running depth of the i-th section of casing, m; _ρmin is the minimum drilling fluid density for the next drilling, g/ ^cm3 ; _km is the hollowing coefficient, ranging from 0 to 1, 1 means full hollowing; _An is the inner cross-sectional area of the nth section of casing, ^mm2 .

(5) Determine whether the internal pressure resistance and tensile strength of the i-th section of the casing meet the verification requirements. If the internal pressure resistance and tensile strength of the i-th section of the casing meet the verification requirements, proceed to step (6); if the internal pressure resistance and tensile strength of the i-th section of the casing do not meet the verification requirements, return to step (2) and select a casing with a higher steel grade or wall thickness to redesign the i-th section.

Taking the first section of casing as an example, for the first section of casing, according to the principle of P _ca1 ≥ Sc × P _ce1 (Sc is the specified anti-extrusion coefficient), select the steel grade and wall thickness of the first section of casing, use the casing strength formula to calculate or find out the casing strength, and list the casing performance parameter table.

1. The method for checking the collapse strength of the top of the first section of casing is as follows:

According to the selected casing, the triaxial collapse strength value of the first section of casing is calculated. If P _ca1 ≥ Sc × P _ce1 , the collapse strength check of the first section of casing is considered normal. Sc is the specified collapse safety factor.

The calculation formula of triaxial collapse strength is:

In the formula, _Pca is the triaxial collapse strength, MPa; _Pco is the collapse strength, MPa; _σa is the axial stress, MPa; _Pi is the liquid column pressure in the pipe, MPa; _Yp is the yield strength of the pipe, MPa.

2. The method for checking the internal pressure strength of the top of the first section casing is as follows:

The triaxial anti-internal pressure strength P _ba1 of the top of the first section of casing is calculated according to the triaxial anti-internal pressure strength formula, and the effective internal pressure P _be1 of the first section of casing is calculated according to the effective internal pressure formula. The safety factor of the first section of casing against internal pressure is: Si ₁ = P _ba1 / P _be1 . If Si ₁ ≥ Si, it is considered that the anti-internal pressure strength of the first section of casing meets the requirements; if Si ₁ < Si, it is considered that the anti-internal pressure strength of the first section of casing does not meet the requirements, and a higher level casing should be selected for tensile design. Si is the specified safety factor against internal pressure.

The calculation formula of triaxial internal compressive strength is:

Wherein, _Pba is the internal pressure resistance, MPa; _Pbo is the internal pressure resistance, MPa; _σa is the axial stress, MPa; _Po is the liquid column pressure outside the tube, MPa; _ro is the outer radius of the casing, mm; _ri is the inner radius of the casing, mm; _Yp is the yield strength of the pipe, MPa.

The calculation formula of effective internal pressure is as follows:

a) Gas Well

①Surface casing and technical casing

The maximum internal pressure at the casing shoe of the surface casing and technical casing can be calculated according to the maximum drilling fluid density used last time, and then the maximum internal pressure and effective internal pressure of the casing at any well depth can be calculated based on the maximum internal pressure.

Maximum internal pressure at the casing shoe: P _bs = 0.00981ρ _max H _s (38)

Maximum internal pressure of casing at any well depth:

Effective internal pressure: P _be =P _bh -0.00981ρ _c H (40

② Production of casing and tail pipe

The production casing and production tail pipe can be considered as being fully filled with natural gas. The maximum internal pressure at any well depth is calculated first, and then the effective internal pressure is calculated based on the maximum internal pressure at any well depth.

Maximum internal pressure of casing at any well depth: P _bh = P _p (41)

Effective internal pressure: P _be =P _bh -0.00981ρ _c H (42)

b) Oil wells

①Surface casing and technical casing

Maximum internal pressure of casing at any well depth:

Effective internal pressure:

② Production of casing and tail pipe

Maximum internal pressure at the casing shoe in production without tubing: P _bs = G _p H _s (45) Maximum internal pressure of casing at any well depth:

Maximum internal pressure produced by oil pipe: P _bh = G _p H _s + 0.00981ρ _w H (47)

Effective internal pressure: P _be =P _bh -0.00981ρ _c H (48)

Wherein, P _bs is the maximum internal pressure at the casing shoe, MPa; H _s is the well depth at the casing shoe, m; G _p is the pressure gradient outside the casing, MPa/m; P _bh is the maximum internal pressure of the casing at any well depth, MPa; ρ _g is the gas with the highest drilling fluid density used during drilling, g/cm ³ ; ρ _w is the completion fluid density, g/cm ³ ; H is the well depth at the calculation point, m; P _be is the effective internal pressure, MPa; ρ _c is the medium in the casing (i.e., oil or natural gas), g/cm ³ ; ρ _max is the maximum drilling fluid density for the next drilling, g/cm ³ .

3. The method for checking the tensile strength of the top of the first section casing is as follows:

The triaxial tensile strength T _a1 of the top of the first section of casing is calculated according to the triaxial tensile strength formula, and the effective tensile force T _e1 of the first section of casing is calculated according to the effective tensile force formula. Then the tensile safety factor of the first section of casing is: St ₁ = T _a1 / T _e1 . If St ₁ ≥ St, it is considered that the tensile strength of the first section of casing meets the requirements; if St ₁ < St, it is considered that the tensile strength of the first section of casing does not meet the requirements. St is the specified safety factor against internal pressure.

The calculation formula for triaxial tensile strength is:

Wherein, _Ta is the triaxial tensile strength, MPa; _Po is the liquid column pressure outside the tube, MPa; _Pi is the liquid column pressure inside the tube, MPa; _Ro is the outer radius of the casing, mm; _Ri is the inner radius of the casing, mm; _To is the tensile strength, MPa.

The calculation formula of effective pulling force is as follows:

a) Vertical Well

Where, _Ten is the effective tension at the top of the nth section of casing, MPa; _Li is the running length of the ith section of casing, m; _ρs is the steel density, g/ ^cm3 ; _ρm is the drilling fluid density, g/ ^cm3 .

b) 2D wellbore

① In the deflection section, the pipe string is in contact with the lower well wall (N>0, N is the number of limited torsion circles)

② In the deflection section, the pipe string is in contact with the upper wellbore wall (N<0)

③Decline well section

④Stable inclination section

T _ei+1 =T _ei +q _ei (cosα+μsinα)(L _i +L _i+1 ) (55)

Wherein, _Tei+1 is the effective tension at the top of the i+1th section casing, MPa; _Tei is the effective tension at the top of the i-th section casing, MPa; _Li is the running length of the i-th section casing, m; Li ₊₁ is the running length of the i+1th section casing, m.

4. The method of redesigning the nth section casing is as follows:

If the tensile strength or internal pressure resistance is not satisfied when the above anti-extrusion design reaches the nth section of casing, a higher level casing should be selected and the tensile strength design of this section of casing should be changed to calculate the lowering length L _on of this section of casing according to the tensile strength of the casing.

The running length L _on of the nth section of casing is calculated by the following formula:

It is also necessary to calculate the running length _Lan of the nth section of casing under triaxial stress, which is calculated as follows:

The triaxial tensile strength _Tan of the nth section of casing is calculated from the triaxial tensile strength. After _Lan is calculated from _Tan and the tensile formula, if Then L _n =L _an ; otherwise repeat the above calculation until Then the casing section's internal pressure resistance and collapse resistance are checked until the designed well depth is met.

(6) Determine whether the i-th section of casing has reached the wellhead. If the i-th section of casing has not reached the wellhead, return to step (3), set i=i+1, calculate the insertion depth of the next section of casing, and check the internal pressure resistance and tensile strength; if the i-th section of casing has reached the wellhead, the calculation is terminated and the design results of all casing strings are output.

In addition, the downhole tubing mechanical analysis calculation module 300 may also include a mechanical analysis result output submodule 340, which is connected to the casing force analysis calculation submodule 320 and the casing strength verification submodule 330, and is configured to graphically output the casing string force condition and strength verification results.

In this embodiment, the setting requirement of the wellbore trajectory target analysis and calculation unit is: to perform target analysis and calculation in real time in combination with the directional well trajectory monitoring data and the wellbore trajectory three-dimensional graphics.

Therefore, as shown in FIG. 1 , the borehole trajectory target calculation module 400 may include a third input data extraction submodule 410 , an anti-collision analysis submodule 420 , a borehole trajectory generation submodule 430 , a borehole trajectory prediction submodule 440 , and a borehole trajectory data output submodule 450 .

The third input data extraction submodule 410 is connected to the drilling basic data storage module 100 and is configured to extract the third input data required for the wellbore trajectory target calculation from the drilling basic data storage module. The third input data includes: wellbore trajectory data.

The anti-collision analysis submodule 420 is connected to the third input data extraction submodule 410 and is configured to calculate the wellbore trajectory distance between the target well and the adjacent well.

The wellbore trajectory generation submodule 430 is connected to the third input data extraction submodule 410 and is configured to automatically generate a three-dimensional wellbore trajectory map according to the wellbore trajectory data.

The borehole trajectory prediction submodule 440 is connected to the borehole trajectory generation submodule 430 and is configured to automatically generate a three-dimensional borehole trajectory prediction map when the well depth, well inclination angle, and azimuth of the predicted borehole trajectory are input.

The wellbore trajectory data output submodule 450 is connected to the wellbore trajectory generation submodule 430 and is configured to output wellbore trajectory data.

The minimum curvature method can be used to calculate the wellbore trajectory, where the calculation formula of the minimum curvature method is shown in equations (57) to (63).

When Δα＝0, , calculate as follows:

Vertical depth increment: ΔH _i = ΔL _i cos α _i (60)

Segment displacement: ΔS _i = ΔL _i sinα _i (61)

N coordinate increment:

E coordinate increment:

Wherein, ΔH _i is the vertical depth increment, m; ΔS _i is the segment displacement, m; ΔN _i is the N coordinate increment, m; ΔE _i is the E coordinate increment, m; α _i is the well inclination angle of the ith measuring point, degrees; Δα is the well inclination angle increment, degrees; is the azimuth of the i-th measuring point, degrees; is the azimuth increment, degrees.

In addition, the borehole trajectory target calculation module 400 may also include a trajectory playback submodule 460, a real-time trajectory submodule 470, and a two-dimensional trajectory drawing submodule 480. Among them, the trajectory playback submodule 460 is connected to the borehole trajectory generation submodule 430, and is configured to automatically simulate and playback the demonstration animation of the drilling borehole generation process according to the three-dimensional borehole trajectory map. The real-time trajectory submodule 470 is connected to the borehole trajectory generation submodule 430, and is configured to generate a demonstration animation of the drilling borehole generation process in real time according to the three-dimensional borehole trajectory map. The two-dimensional trajectory drawing submodule 480 is connected to the borehole trajectory generation submodule 430, and is configured to automatically draw the borehole trajectory into a vertical projection map and a horizontal projection map.

In this embodiment, the setting requirements of the cementing construction calculation module are: to optimize the cementing construction parameters in combination with the well depth, cementing design, and cementing pipe string, to provide a theoretical basis for cementing design, and to guide on-site real-time construction.

Therefore, as shown in FIG. 1 , the cementing construction calculation module 500 may include a fourth input data extraction submodule 510 , a cementing construction parameter calculation submodule 520 , and a cementing construction dynamic simulation submodule 530 .

The fourth input data extraction submodule 510 is connected to the drilling basic data storage module 100 and is configured to extract the fourth input data required for cementing construction calculation from the drilling basic data storage module. The fourth input data includes: casing data, drilling tool assembly data, wellbore structure data and the maximum number of twisting circles of the drilling tool allowed when the hanging casing is inverted. The fourth input data extraction submodule can also obtain the hanging casing steel ball falling time data by manual input.

The cementing construction parameter calculation submodule 520 is connected to the fourth input data extraction submodule 510 and is configured with a cementing construction calculation model, which can calculate cementing construction parameters. The cementing construction parameters may include casing elongation, drilling tool elongation, displacement, maximum allowable casing lowering speed, hanging casing steel ball drop time, and maximum allowable drilling tool torsion number when hanging casing is inverted.

The cementing construction dynamic simulation submodule 530 is connected to the cementing construction parameter calculation submodule 520 and is configured with a cementing construction numerical model, which can visualize and simulate the cementing construction operation process.

In this embodiment, the setting requirements of the stuck drill accident analysis and calculation module are: combining real-time data such as well depth, drilling tool assembly, hydraulic parameters, etc., analyzing and calculating the stuck points, making judgments on the downhole status, and providing a reference basis for the next step of handling the stuck drill or other operations.

Therefore, as shown in FIG. 1 , the stuck drill accident analysis calculation module 600 may include a fifth input data extraction submodule 610 , a twisting circle calculation submodule 620 , a stuck point position calculation submodule 630 , a release agent calculation submodule 640 , a pressure reduction calculation submodule 650 , and a drop time calculation submodule 660 .

The fifth input data extraction submodule 610 is connected to the drilling basic data storage module 100 and is configured to extract the fifth input data required for the stuck pipe accident analysis and calculation from the drilling basic data storage module. The fifth input data includes: drilling tool assembly data, wellbore structure data and K value coefficient data. The fifth input data extraction submodule can also obtain drilling tool strength data by manual input.

The twisting circle number calculation submodule 620 is connected to the fifth input data extraction submodule 610 and is configured with a twisting circle number calculation model, which can calculate the twisting circle number. The twisting circle number calculation steps are as follows:

(1) Select the torsion state of the drill bit and calculate the number of torsion circles. The torsion state of the drill bit includes: the large hook hanging weight is equal to the normal drilling hanging weight, the large hook hanging weight is less than the normal drilling hanging weight, and the large hook hanging weight is less than the normal drilling hanging weight. Each state corresponds to a formula.

a) When the hanging weight of the big hook is equal to the weight of the free drill tool in the drilling fluid, the formula for calculating the limited number of torsion turns is as follows:

Where, N is the number of limited torsion circles; H is the depth of the stuck point, that is, the length of the drill string above the stuck point, m; D is the outer diameter of the drill string above the stuck point, cm; G is the shear coefficient, equal to 8×10 ⁴ MPa; τ is the limited shear force, MPa; σ ₁ is the yield strength of steel, MPa; σ ₂ is the actual tensile stress on the dangerous section of the drill string, MPa; q is the weight of each meter of drill string in drilling fluid, kg/m; F is the cross-sectional area of the pipe body, cm ² ; S is the safety factor, which is taken as 1.5.

b) When the hanging weight of the large hook is less than the weight of the free drill bit in the drilling fluid, the formula (64) can be used to calculate the limited number of torsion circles, but the length of the free drill bit should be calculated based on the length of the drill bit above the neutral point.

c) When the weight of the big hook is greater than the weight of the free drill tool in the drilling fluid, the formula for calculating the limited number of torsion turns is as follows:

Where W is the weight of the large hook, KN; other symbols are the same as those in formula (64). The data listed in the standard table of the number of twisting turns cannot be applied at this time.

(2) Through the drill rod limited torsion circle table, the _σ1 value and q value of drill rods of different steel grades and different outer diameters corresponding to the calculated torsion circle values are inferred and added to the drill rod database.

The stuck point location calculation submodule 630 is connected to the fifth input data extraction submodule 610 and is configured with a stuck point location calculation model, which can calculate the stuck point location according to the well type and drilling tool type of the target well. The stuck point location calculation steps are as follows:

(1) Determine the well type and drilling tool type of the target well through the well basic data, and determine the calculation formula for the stuck point location. Well types include vertical wells and inclined wells (including horizontal wells), and drilling tool types include single drill string and composite drill string.

(2) Use the corresponding calculation formula to determine the location of the card point.

a) If the stuck point is in the vertical well section and there is a single drill string in the well, the stuck point position can be calculated based on the elastic elongation of the drill string under a certain tension. The calculation formula is as follows:

Where, L is the length of the free drill string, m; K = EA, is the calculation coefficient; A is the cross-sectional area of the free drill string, ^cm2 ; E is the elastic coefficient of steel, 2.1 × ¹⁰⁵ MPa; △x is the elongation of the free drill string under the action of △f force, cm; △f is the difference between the two pulling forces on the free drill string that exceed its own hanging weight, KN.

If the influence of joints and thickened parts is taken into account, the following formula is formed:

Where S is the influence coefficient of the drill pipe body, joint and thickened part; Aa is the cross-sectional area of the drill pipe body in cm ² ; Ab is the cross-sectional area of the joint in cm ² ; Ac is the cross-sectional area of the thickened part in cm ² .

b) If the stuck point is in the inclined well section and there is a single drill string in the well, the stuck point should be calculated according to formula (68):

Where S is the influence coefficient of the drill pipe body, joints and thickened parts; μ is the friction coefficient between the drill string and the well wall, which should be input between 2 and 4; Δα is the calculated average value of the well inclination angle up to the nth section, in degrees.

c) If the stuck point is in the vertical well section, and the well is a composite drill string, it needs to be divided into several sections according to the outer diameter and wall thickness of the drill string. The length of each section is L ₁ , L ₂ , L ₃ , etc. Under a certain pulling force △f, calculate the pull-up value of each section △x ₁ , △x ₂ , △x ₃ , etc. Therefore, the calculation method of the stuck point position of the composite drill string is as follows:

First, use formula (69) to calculate the elongation △ _xi of each section of the drilling tool starting from the wellhead. Formula (69) is as follows:

Where △ _xi is the elongation of the i-th drill string under the action of the tension △f, cm; △f is the difference between the two tensions that the composite drill string is subjected to that exceed its own suspended weight, KN; K is the calculation coefficient; _Li is the length of the i-th drill string, m.

Then the measured value △x is compared with the calculated values _∑Δxi and ∑Δxi ₊₁ respectively, i=1,2,3,..., until _∑xi ＜Δx＜∑Δxi ₊₁ , and the stuck point position is determined to be on the i+1th section of the drilling tool.

If the measured value △x is less than the calculated value △x ₁ , the stuck point is on the first section of the drilling tool. The stuck point can be directly calculated using equation (66) or (67). The K value at this time is the K value of the first section of the drilling tool.

If the measured value △x is greater than the calculated value _∑Δxi and less than ∑Δxi ₊₁ , the stuck point is on the i+1th section of the drilling tool, and the stuck point position is calculated using formula (70).

L=∑L _i +K _i+1 (Δx-∑Δx _i )/Δf (70)

Wherein, L is the position of the stuck point, m; _Li is the length of the i-th section of the drill string, m; Ki ₊₁ is the calculation coefficient of the i+1-th section of the drill string; △x is the measured value of the elongation of the composite drill string under the action of the tension △f, m; △ _xi is the elongation of the i-th section of the drill string under the action of the tension △f, m; △f is the difference between the two tensions that the composite drill string is subjected to that exceed its own suspended weight, KN.

d) If the stuck point is in the inclined well section and the well contains a composite drill string, the stuck point position is calculated as follows:

First, calculate the elongation △ _xi of each section of the drilling tool starting from the wellhead, and then compare the measured value △x with the calculated value _∑Δxi , ∑Δxi ₊₁ , i=1,2,3,..., until _∑xi ＜Δx＜∑Δxi ₊₁ , and determine that the stuck point is on the i+1th section of the drilling tool.

If the measured value △x is less than the calculated value △x ₁ , the stuck point is on the first section of the drilling tool. The stuck point can be directly calculated using formula (68). The K value at this time is the K value of the first section of the drilling tool.

If the stuck point is in the i+1th section and is also a deviated section, the position of the stuck point is calculated according to (71);

Where, S _i+1 is the influence coefficient of the drill pipe body, joints and thickened parts of the i+1th section; μ is the friction coefficient between the drill string and the well wall, which should be input between 2 and 4; Δα is the calculated average value of the well inclination angle up to the nth section, in degrees.

The release agent calculation submodule 640 is connected to the fifth input data extraction submodule 610 and is configured with a release agent dosage calculation model and a maximum pump pressure calculation model, which can respectively calculate and output the total release agent dosage and the release agent pumping maximum pump pressure.

The calculation steps for the amount of the deblocking agent are as follows:

First, obtain the drill bit position, stuck point depth, and the inner and outer diameters and lengths of the daily drilling tool assembly data. The drill bit position and stuck point depth are used to calculate the drill string length of the stuck section, and the inner and outer diameters and lengths of the daily drilling tool assembly data are used to calculate the annular capacity and the pipe content of the stuck section.

Then calculate the total amount of de-stuck agent, which is equal to the annular space volume expected to be soaked and the volume inside the drill string. The annular space volume is the annular space volume from the drill bit to the stuck point, and a certain additional amount must be added. Therefore, the total amount of de-stuck agent can be calculated using the following formula:

Wherein, Q is the total amount of de-stuck agent, m ³ ; _Q1 is the annulus capacity of the stuck section, m ³ ; _Q2 is the volume of the pipe inside the stuck section, m ³ ; _Q3 is the reserved displacement volume, m ³ ; K is the accessory coefficient, which is generally taken as 1.2; H is the drill string length of the stuck section, m; D is the drill bit diameter, m; _d1 is the outer diameter of the drill collar or drill pipe, m; _d2 is the inner diameter of the drill collar or drill pipe, m.

If a composite drill string or a stepped wellbore is used, the calculation should be carried out in sections according to the different well diameters and the different inner and outer diameters of the pipe string, and the total amount can be obtained by adding up.

The calculation steps for the maximum pump pressure of the release agent are as follows:

Determine whether the density of the de-stuck agent is less than the density of the drilling fluid. If so, determine the relationship between the total amount of de-stuck agent used and the volume of the drill string. Combined with the daily drill tool assembly data, the height of the liquid column of the de-stuck agent in the drill string under different conditions can be obtained, thereby calculating the maximum pump pressure.

If the density of the de-stuck agent is similar to that of the drilling fluid, the pump pressure will not change significantly and does not need to be calculated. If the density of the de-stuck agent is lower than that of the drilling fluid, there will be a pressure difference during displacement, and the maximum pump pressure is the sum of the circulating pump pressure and the pressure difference of the liquid column inside and outside the pipe. If the total amount of the de-stuck agent is less than the volume of the drill string, the maximum pump pressure is reached when the de-stuck agent is completely pumped into the drill string. If the total amount of the de-stuck agent is greater than the volume of the drill string, the maximum pump pressure is reached when the de-stuck agent reaches the drill bit. The maximum pump pressure can be calculated using the following formula:

P＝p ₁ +p ₂ =P ₁ +0.01(ρ ₁ -ρ ₂ )h (73)

Wherein, P is the maximum pump pressure, MPa; _p1 is the circulating pump pressure, MPa; _p2 is the liquid column pressure difference between the release agent and the drilling fluid, MPa; _ρ1 is the drilling fluid density, g/ ^cm3 ; _ρ2 is the density of the release agent, g/ ^cm3 ; h is the liquid column height of the release agent in the drill string, m.

The pressure reduction calculation submodule 650 is connected to the fifth input data extraction submodule 610 and is configured with a U-tube effect pressure reduction calculation model, which can calculate the output pressure reduction parameters. The calculation steps of the U-tube effect pressure reduction method are as follows:

The pressure difference in the sticking drill is the main cause of the stuck drill. It is believed that reducing the pressure difference or even negative pressure difference can release the stuck drill. This method is subject to the following conditions when it is implemented: the well must have been casingd and have complete well control equipment; there must be no high-pressure layer or collapsed layer in the open hole section; the drill string is not connected to a back pressure valve and can be reversed; there is no possibility of blocking the drill bit water hole.

After low-density drilling fluid is pumped in through reverse circulation and the pressure is released back to achieve internal and external balance, the following relationship is observed when the pressure inside and outside the drill string is balanced:

Liquid column height of low-density drilling fluid _M2 :

H ₁ =(ρ ₁ -ρ ₃ )H/(ρ ₁ -ρ ₂ ) (74)

_M2 is the volume of liquid column in the annulus:

Q ₁ = V ₁ H ₁ /1000 (75)

Height of air column A in drill pipe:

H ₂ =(ρ ₁ -ρ ₃ )H/ρ ₁ (76)

The volume of air in the drill pipe:

_Q2 ＝ _V2H2 / ₁₀₀₀ (77)

The total volume of low-density drilling fluid M2 that needs to be pumped in:

Q ₀ =Q ₁ +Q ₂ (78)

Maximum pressure difference when pumping low-density drilling fluid:

P＝(ρ ₁ -ρ ₂ )H/100 (79)

Wherein, _M1 is the original slurry in the well; _M2 is the low-density drilling fluid replaced; H is the vertical depth of the target layer; _H1 is the vertical height of the liquid column of _M2 drilling fluid in the annulus when the pressure inside and outside the drill string is balanced, m; _ρ1 is the density of _M1 drilling fluid, g/ ^cm3 ; _ρ2 is the density of _M2 drilling fluid, g/ ^cm3 ; _ρ3 is the equivalent density of the formation pressure of the target layer or the equivalent density achieved by design, g/ ^cm3 ; _Q1 is the volume of _M2 drilling fluid when the pressure inside and outside the drill string is balanced after the pressure is backed out, ^m3 ; _V1 is the volume per meter of the annulus, L/m; _V2 is the volume per meter of the drill string, L/m; _H2 is the vertical height of the air column in the drill string when the pressure inside and outside the drill string is balanced, m; _Q2 is the volume of air A when the pressure inside and outside the drill string is balanced after the pressure is backed out, ^m3 ; Q ₀ is the volume of _M2 drilling fluid in the final stage of reverse circulation pumping low-density drilling fluid, ^m3 ; P is the maximum pressure difference when pumping low-density drilling fluid, MPa.

The falling time calculation submodule 660 is connected to the fifth input data extraction submodule 610 and is configured with a falling time calculation model, which can calculate and output the falling time of the pitched ball. The falling time calculation steps of the pitched ball are as follows:

When the ball is only allowed to fall freely in the drill pipe and cannot be pumped, the formula for calculating the falling time is:

In the formula, t is the falling time of the copper ball, s; v is the falling speed of the copper ball, m/s; g is the acceleration of gravity, which is 980cm/ ^s2 ; D is the diameter of the copper ball, _cm ; _ρ1 is the density of the copper ball, g/ ^cm3 ; _ρ2 is the density of the drilling fluid, g/ ^cm3 ; μ is the plastic viscosity of the drilling fluid, 10mPa.s; H is the depth of the ball seat, m; _K2 is the time additional coefficient, which is 1.5~2.0.

In this embodiment, the setting requirements of the overflow well-killing construction calculation module are: combining the real-time drilling tool assembly, wellbore structure, and drilling fluid performance in the wellbore, selecting relevant parameters and well-killing methods, and generating a well-killing construction order with one click.

Therefore, as shown in Figure 1, the overflow well killing construction calculation module 700 may include a sixth input data extraction submodule 710, a well killing construction calculation submodule 720, a well killing construction order generation submodule 730, a well killing construction dynamic simulation submodule 740 and a well killing operation drawing submodule 750.

The sixth input data extraction submodule 710 is connected to the drilling basic data storage module 100 and is configured to extract the sixth input data required for the overflow well killing construction calculation from the drilling basic data storage module. The sixth input data includes drilling tool assembly data, wellbore structure data, low pump stroke test data, wellhead pressure test data and overflow key data.

The well killing construction calculation submodule 720 is connected to the sixth input data extraction submodule 710 and is configured with a well killing construction calculation model, and is capable of calculating well killing construction parameters.

The well killing construction form generation submodule 730 is connected to the well killing construction calculation submodule 720 and is configured to automatically generate a well killing construction form according to the calculation results of the well killing construction parameters.

The well killing construction dynamic simulation submodule 740 is connected to the well killing construction calculation submodule 720 and is configured with a well killing construction numerical model, which can visualize and simulate the well killing operation process of the wellbore.

The well killing operation drawing submodule 750 is connected to the well killing construction dynamic simulation submodule 740 and is configured to graphically output the comparison result between the real-time well killing data and the theoretical well killing data.

In this embodiment, the computing system may also include a pressure consumption calculation module, the pressure consumption calculation model is connected to the drilling basic data storage module, and includes a seventh input data extraction submodule, a pressure consumption calculation submodule, a pressure consumption prediction submodule, an actual pressure consumption curve drawing submodule and a pressure consumption prediction curve drawing submodule.

The seventh input data extraction submodule is connected to the drilling basic data storage module and is configured to extract the seventh input data required for pressure loss calculation from the drilling basic data storage module. The seventh input data includes: drilling tool assembly data, wellbore structure data, drilling real-time data and drilling fluid performance data.

The pressure loss calculation submodule is connected to the seventh input data extraction submodule and is configured to automatically calculate and output the inner-tube pressure loss, outer-tube pressure loss and total circulation pressure loss of the current downhole drilling tool assembly.

The actual pressure consumption curve drawing submodule is connected to the pressure consumption calculation submodule, and is configured to automatically draw the calculation result of the pressure consumption calculation submodule as an actual pressure consumption curve.

The pressure loss prediction submodule is connected to the seventh input data extraction submodule and is configured to automatically calculate and output the cyclic pressure loss prediction results of the section to be drilled under different working conditions.

The pressure consumption prediction curve drawing submodule is connected to the pressure consumption prediction submodule and is configured to automatically draw the calculation result of the pressure consumption prediction curve drawing submodule into a pressure consumption prediction curve.

In summary, the beneficial effects of the present invention include at least one of the following:

(1) The present invention can quickly analyze and calculate drilling parameters based on basic drilling data in real time, which greatly shortens the calculation time and avoids the errors of manual calculation;

(2) The present invention can utilize the real-time data and dynamic data of the drilling site to perform real-time calculation and analysis on the parameters related to the drilling project, thereby providing data support for safe drilling. At the same time, the drilling project parameters can be optimized based on the analysis results, thereby providing real-time guidance for drilling production.

(3) The present invention realizes timely calculation and analysis of drilling engineering related parameters, which is of great significance to meet the ever-increasing demand for drilling engineering and plays an increasingly important role in improving drilling efficiency and enhancing drilling operation safety;

(4) The present invention is the first in China to realize the integrated application of multiple analysis and calculation software for drilling engineering on the same engineering technology information platform, which can provide unified database support, remote communication support, numerical calculation support, graphic visualization support, etc.

Although the present invention has been described above in conjunction with the exemplary embodiments and the accompanying drawings, it should be apparent to those skilled in the art that various modifications may be made to the above-described embodiments without departing from the spirit and scope of the claims.

Claims (11)

1. A computing system for simulating wellbore engineering construction, characterized in that the computing system includes a drilling basic data storage module, a wellbore cleaning analysis and calculation module, a wellbore string mechanics analysis and calculation module, a wellbore trajectory target calculation module, a cementing construction calculation module, a stuck drill accident analysis and calculation module, and an overflow well killing construction calculation module, wherein: The drilling basic data storage module is configured to acquire and store real-time drilling basic data of the target well during the drilling process; The wellbore cleaning analysis and calculation module is connected to the drilling basic data storage module and is configured to be able to perform wellbore cleanliness calculation for well sections with different inclinations; The mechanical analysis and calculation module of the wellbore pipe string is connected to the drilling basic data storage module and is configured to perform force analysis and calculation, deformation analysis and calculation and safety evaluation calculation on the pipe string in different operation processes; The wellbore trajectory target calculation module is connected to the drilling basic data storage module and is configured to be able to perform target analysis and calculation in combination with the directional well trajectory monitoring data and the wellbore trajectory three-dimensional graphics; The cementing construction calculation module is connected to the drilling basic data storage module and is configured to be able to analyze and calculate cementing parameters in combination with real-time drilling basic data; The stuck drill accident analysis and calculation module is connected to the drilling basic data storage module and is configured to be able to analyze and calculate stuck drill parameters in combination with real-time drilling basic data; The overflow well-killing construction calculation module is connected to the drilling basic data storage module and is configured to be able to analyze and calculate well-killing construction parameters in combination with real-time drilling basic data; The wellbore cleaning analysis and calculation module includes a screening submodule, a first input data extraction submodule, a wellbore cleanliness calculation submodule and a wellbore cleanliness determination submodule, wherein: The screening submodule is configured to automatically screen and determine the clean zone to which the well section to be calculated belongs according to the well inclination angle of the well section to be calculated, and output the partition result; The first input data extraction submodule is connected to the drilling basic data storage module and the screening submodule respectively, and is configured to extract the first input data required for the calculation of the cleanliness condition of the well section to be calculated from the drilling basic data storage module based on the partition result; The wellbore cleanliness calculation submodule includes a selection unit, a first clean area calculation unit, a second clean area calculation unit and a third clean area calculation unit, wherein: The selection unit is connected to the screening submodule and is configured to control one of the first cleaning area calculation unit, the second cleaning area calculation unit and the third cleaning area calculation unit to be connected to the first input data extraction submodule according to the partition result of the well section to be calculated, The first cleaning zone calculation unit is configured with a wellbore cleaning calculation model for a low-angle well section, and is capable of calculating and outputting a first calculation result, wherein the first calculation result includes a drill cuttings particle transmission ratio, The second cleaning zone calculation unit is configured with a wellbore cleaning calculation model for a medium-angle well section, and is capable of calculating and outputting a second calculation result, wherein the second calculation result includes a ratio of the cuttings bed thickness to the wellbore diameter and an annulus stop return speed. The third cleaning zone calculation unit is configured with a wellbore cleaning calculation model for a high-angle well section, and is capable of calculating and outputting a third calculation result, wherein the third calculation result includes a ratio of the cuttings bed thickness to the wellbore diameter; The wellbore cleanliness determination submodule is connected to the wellbore cleanliness calculation submodule and is configured to analyze and determine the cleaning effect of the well section to be calculated according to the calculation result output by the wellbore condition calculation module; The calculation model of wellbore cleaning in the low-angle well section is shown in equations (1) to (3): In formulas (1) to (3), V _sx is the settling velocity of drill cuttings, m/s; d _s is the equivalent diameter of drill cuttings, cm; ρ _s is the density of drill cuttings, g/cm ³ ; ρ _m is the density of drilling fluid, g/cm ³ ; is the cuttings particle shape factor, dimensionless; _Va is the drilling fluid annular return velocity, m/s; _Qa is the drilling fluid flow rate, L/s; _Dh is the wellbore diameter, mm; _Dp is the drill pipe outer diameter, mm; Rt is the cuttings particle transmission ratio, dimensionless; The calculation model for wellbore cleaning in the medium-angle well section is shown in equations (4) to (7): h _zx = 0.015D _h (AV + 6·15AV ^0.5 )(1 + 0.587E)(V _lzx -V _a ) Formula (5) H _zx =(h _zx /D _h )×100% Formula (6) In formulas (4) to (7), V _lzx is the critical annular return velocity of the medium-inclined well section, m/s; AV is the apparent viscosity of the drilling fluid, mPa.s; θ is the well inclination angle, degree; h _zx is the thickness of the cuttings bed in the medium-inclined well section, mm; H _zx is the relative thickness of the cuttings bed in the medium-inclined well section, dimensionless; V _p is the annular stop return velocity, m/s; A _bed is the cross-sectional area of the cuttings bed perpendicular to the wellbore axis, mm ² ; C is the cuttings concentration in the cuttings bed, dimensionless; g is the gravitational acceleration, m/s ² ; E is the eccentricity of the drill string, dimensionless; L is the cross-sectional width of the cuttings bed perpendicular to the wellbore axis, mm; η is the friction coefficient between the cuttings bed and the lower wellbore wall, which is taken as 0.2; PV is the plastic viscosity, mPa; YP is the dynamic shear force, Pa; K is the consistency coefficient, Pa·s ⁿ ; n is the fluidity index, dimensionless; The calculation model of wellbore cleaning in the highly deviated well section is shown in equations (15) to (20): H _dx = (h _dx / D _h ) × 100% Formula (20) In formulas (15) to (20), V _Ldx is the critical annular return velocity in the highly deviated well section, m/s; V _sd is the cuttings settling velocity in the highly deviated well section, m/s; V _jx is the mechanical penetration rate, m/h; _Cang is the well inclination correction coefficient, dimensionless; C _size is the cuttings size correction coefficient, dimensionless; C _denF is the drilling fluid density correction coefficient, dimensionless; C _rpm is the drill string speed correction coefficient, dimensionless; Q _Ldx is the critical annular displacement without cuttings bed, L/s; A' _bed is the cuttings bed area, mm ² ; h _dx is the cuttings bed thickness in the highly deviated well section, mm; Q _a is the drilling fluid displacement, L/s; H _dx is the relative thickness of the cuttings bed in the highly deviated well section, dimensionless. 2. The computing system for simulating wellbore engineering construction according to claim 1 is characterized in that the borehole cleaning analysis calculation module also includes a borehole cleaning analysis submodule, which is connected to the borehole cleanliness calculation submodule and is configured to output the borehole cleaning analysis curve of the target well to judge the borehole cleanliness, and includes a small-inclination cuttings transmission ratio curve drawing unit, a medium-inclination critical annulus return velocity and annulus return velocity curve drawing unit, a medium-inclination cuttings bed relative thickness curve drawing unit, a high-inclination critical annulus return velocity and annulus return velocity curve drawing unit, a critical annulus displacement curve drawing unit, and a high-inclination cuttings bed relative thickness curve drawing unit. 3. The computing system for simulating wellbore engineering construction according to claim 1, characterized in that the wellbore string mechanical analysis and calculation module comprises a second input data extraction submodule, a casing force analysis and calculation submodule and a casing strength verification submodule, wherein: The second input data extraction submodule is connected to the drilling basic data storage module and is configured to extract the second input data required for the wellbore string mechanics calculation from the drilling basic data storage module; The casing force analysis calculation submodule is connected to the second input data extraction submodule and is configured with a drilling tool force calculation model, which can calculate and output the axial tension, tensile strength, external extrusion pressure, collapse strength, internal pressure and internal pressure resistance of the casing string in different production periods; The casing strength verification submodule is connected to the second input data extraction submodule and is configured with a casing strength verification algorithm, and is capable of verifying casing string design data. 4. The computing system for simulating wellbore engineering construction according to claim 3, wherein the casing strength verification algorithm comprises the following steps: (1) Obtaining original data for casing string design; (2) Calculate the effective external extrusion pressure p _ce1 at the casing shoe based on the original design data; (3) Select the i-th section of casing to be checked according to the load and geometric constraints, 1≤i≤m, when i=1, the section of casing is located at the bottom of the well, when i=m, the section of casing is located at the wellhead; (4) Calculate the running depth _Li of the i-th section of casing; (5) Determine whether the internal pressure resistance and tensile strength of the i-th section of the casing meet the verification requirements. If the internal pressure resistance and tensile strength of the i-th section of the casing meet the verification requirements, proceed to step (6); otherwise, return to step (2) and select a casing with a higher steel grade or wall thickness to redesign the i-th section; (6) Determine whether the i-th section of casing has reached the wellhead. If the i-th section of casing has not reached the wellhead, return to step (3), set i=i+1, calculate the insertion depth of the next section of casing, and check the internal pressure resistance and tensile strength; if the i-th section of casing has reached the wellhead, the calculation is terminated and the design results of all casing strings are output. 5. The computing system for simulating wellbore engineering construction according to claim 3 is characterized in that the mechanical analysis calculation module of the wellbore tubing string also includes a mechanical analysis result output submodule, and the mechanical analysis result output submodule is connected to the casing stress condition calculation submodule and the casing strength verification submodule, and is configured to be able to graphically output the casing string stress condition and strength verification results. 6. The computing system for simulating wellbore engineering construction according to claim 1, characterized in that the wellbore trajectory target calculation module includes a third input data extraction submodule, an anti-collision analysis submodule, a wellbore trajectory generation submodule, a wellbore trajectory prediction submodule and a wellbore trajectory data output submodule, wherein: The third input data extraction submodule is connected to the drilling basic data storage module and is configured to extract the third input data required for the wellbore trajectory target calculation from the drilling basic data storage module; The anti-collision analysis submodule is connected to the third input data extraction submodule and is configured to calculate the wellbore trajectory distance between the target well and the adjacent well; The wellbore trajectory generation submodule is connected to the third input data extraction submodule and is configured to automatically generate a three-dimensional wellbore trajectory map based on the wellbore trajectory data; The wellbore trajectory prediction submodule is connected to the wellbore trajectory generation submodule and is configured to automatically generate a three-dimensional wellbore trajectory prediction map when the well depth, well inclination angle, and azimuth angle of the predicted wellbore trajectory are input; The wellbore trajectory data output submodule is connected to the wellbore trajectory generation submodule and is configured to output wellbore trajectory data. 7. The computing system for simulating wellbore engineering construction according to claim 6, characterized in that the wellbore trajectory on-target computing module further comprises a trajectory playback submodule, a real-time trajectory submodule and a two-dimensional trajectory map drawing submodule, wherein: The trajectory playback submodule is connected to the wellbore trajectory generation submodule and is configured to automatically simulate and playback a demonstration animation of the wellbore drilling generation process according to the three-dimensional wellbore trajectory diagram; The real-time trajectory submodule is connected to the wellbore trajectory generation submodule and is configured to generate a demonstration animation of the wellbore generation process while drilling in real time according to the three-dimensional wellbore trajectory diagram; The two-dimensional trajectory map drawing submodule is connected to the wellbore trajectory generating submodule and is configured to automatically draw the wellbore trajectory into a vertical projection map and a horizontal projection map. 8. The computing system for simulating wellbore engineering construction according to claim 1, characterized in that the cementing construction computing module comprises a fourth input data extraction submodule, a cementing construction parameter computing submodule and a cementing construction dynamic simulation submodule, wherein: The fourth input data extraction submodule is connected to the drilling basic data storage module and is configured to extract the fourth input data required for cementing construction calculation from the drilling basic data storage module; The cementing construction parameter calculation submodule is connected to the fourth input data extraction submodule and is configured with a cementing construction calculation model, capable of calculating cementing construction parameters; The cementing construction dynamic simulation submodule is connected to the cementing construction parameter calculation submodule and is configured with a cementing construction numerical model, which can visualize and simulate the cementing construction operation process. 9. The computing system for simulating wellbore engineering construction according to claim 1, characterized in that the stuck drill accident analysis computing module comprises a fifth input data extraction submodule, a twisting circle calculation submodule, a stuck point location calculation submodule, a de-stuck agent calculation submodule, a pressure reduction calculation submodule and a falling time calculation submodule, wherein: The fifth input data extraction submodule is connected to the drilling basic data storage module and is configured to extract the fifth input data required for the stuck pipe accident analysis and calculation from the drilling basic data storage module; The twisting circle number calculation submodule is connected to the fifth input data extraction submodule and is configured with a twisting circle number calculation model, capable of calculating the twisting circle number; The stuck point location calculation submodule is connected to the fifth input data extraction submodule and is configured with a stuck point location calculation model, capable of calculating the stuck point location according to the well type and drilling tool type of the target well; The jamming agent calculation submodule is connected to the fifth input data extraction submodule and is configured with a jamming agent dosage calculation model and a maximum pump pressure calculation model, which can respectively calculate and output the total jamming agent dosage and the jamming agent pumping maximum pump pressure; The voltage reduction calculation submodule is connected to the fifth input data extraction submodule and is configured with a U-tube effect voltage reduction calculation model, capable of calculating output voltage reduction parameters; The falling time calculation submodule is connected to the fifth input data extraction submodule and is configured with a falling time calculation model, which can calculate and output the falling time of the ball during pitching and holding. 10. The computing system for simulating wellbore engineering construction according to claim 1, characterized in that the overflow well-killing construction computing module comprises a sixth input data extraction submodule, a well-killing construction computing submodule, a well-killing construction order generation submodule, a well-killing construction dynamic simulation submodule and a well-killing operation drawing submodule, wherein: The sixth input data extraction submodule is connected to the drilling basic data storage module and is configured to extract the sixth input data required for overflow well killing construction calculation from the drilling basic data storage module; The well killing construction calculation submodule is connected to the sixth input data extraction submodule and is configured with a well killing construction calculation model, capable of calculating well killing construction parameters; The well killing construction form generation submodule is connected to the well killing construction calculation submodule and is configured to automatically generate a well killing construction form according to the calculation results of the well killing construction parameters; The well killing construction dynamic simulation submodule is connected to the well killing construction calculation submodule and is configured with a well killing construction numerical model, which can visualize and simulate the well killing operation process of the wellbore; The well-killing operation drawing submodule is connected to the well-killing construction dynamic simulation submodule, and is configured to graphically output the comparison result between the real-time well-killing data and the theoretical well-killing data. 11. The computing system for simulating wellbore engineering construction according to claim 1 is characterized in that the computing system also includes a pressure loss calculation module, the pressure loss calculation model is connected to the drilling basic data storage module, and includes a seventh input data extraction submodule, a pressure loss calculation submodule, a pressure loss prediction submodule, an actual pressure loss curve drawing submodule and a pressure loss prediction curve drawing submodule, wherein, The seventh input data extraction submodule is connected to the drilling basic data storage module and is configured to extract the seventh input data required for pressure loss calculation from the drilling basic data storage module; The pressure loss calculation submodule is connected to the seventh input data extraction submodule and is configured to automatically calculate and output the inner-tube pressure loss, outer-tube pressure loss and total circulation pressure loss of the current downhole drilling tool assembly; The actual pressure consumption curve drawing submodule is connected to the pressure consumption calculation submodule and is configured to automatically draw the calculation result of the pressure consumption calculation submodule into an actual pressure consumption curve; The pressure loss prediction submodule is connected to the seventh input data extraction submodule and is configured to automatically calculate and output the cyclic pressure loss prediction results of the section to be drilled under different working conditions; The pressure consumption prediction curve drawing submodule is connected to the pressure consumption prediction submodule and is configured to automatically draw the calculation result of the pressure consumption prediction curve drawing submodule into a pressure consumption prediction curve.