Method and equipment for optimizing hydraulic parameters of deepwater managed pressure drilling in real time

Abstract

The present invention provides a method and equipment for optimizing hydraulic parameters of deepwater managed pressure drilling in real time. The method comprises: acquiring overflow parameters in the current drilling process, performing preprocessing and feature extraction on the overflow parameters, and inputting the overflow parameters into trained support vector machine identification models for overflow judgment. If overflow occurs at the current drilling depth, reducing an opening of a throttle valve and increasing a displacement of a submarine pump, measuring a wellhead back pressure and calculating a bottom hole pressure, and judging whether overflow continues is performed. If overflow continues, mixing high-density drilling fluid with the original drilling fluid, pumping the mixture into a wellbore annulus, and performing the above operations of reducing the throttle valve opening, increasing the displacement of the submarine pump, calculating the bottom hole pressure and judging whether overflow continues is performed until overflow no longer occurs.

Claims

1. A method for optimizing hydraulic parameters of deepwater managed pressure drilling in real time, the method comprising: acquiring overflow parameters in the current drilling process in real time, performing preprocessing and feature extraction on the overflow parameters, and inputting the overflow parameters after preprocessing and feature extraction into trained support vector machine identification models for overflow judgment; and when it is judged that overflow occurs at the current drilling depth, reducing the opening of a throttle valve on a throttle pipeline, increasing a wellhead back pressure and increasing a displacement of a submarine pump and a displacement of drilling fluid; measuring the wellhead back pressure and calculating a bottom hole pressure according to the measured wellhead back pressure; judging whether overflow continues to occur in case that the calculated bottom hole pressure does not fall into a safety window; and under the condition that overflow continues to occur, mixing high-density drilling fluid with the original drilling fluid, pumping the mixture into a wellbore annulus from a drill pipe, and performing the above operations of reducing the opening of the throttle valve, increasing the displacement of the submarine pump, calculating the bottom hole pressure and judging whether overflow continues to occur until overflow no longer occurs, wherein, the trained support vector machine identification models comprise: a flow identification model, a mud pit increment identification model and a standpipe pressure identification model; and the step of acquiring overflow parameters in the current drilling process in real time, performing preprocessing and feature extraction on the overflow parameters, and inputting the overflow parameters after preprocessing and feature extraction into trained support vector machine identification models for overflow judgment comprises: acquiring a flow differential of an inlet and an outlet, a mud pit increment and a standpipe pressure in the current drilling process in real time, performing preprocessing and feature extraction on the flow differential, the mud pit increment and the standpipe pressure, inputting the flow differential, the mud pit increment and the standpipe pressure after preprocessing and feature extraction into corresponding support vector machine identification models for overflow judgment, and processing an overflow probability under each identification model by an information fusion model to judge whether overflow occurs at the current drilling well depth.

2. The method according to claim 1, wherein the step of calculating the bottom hole pressure according to the measured wellhead back pressure comprises: determining flow calculation parameters after overflow of the managed pressure drilling; determining complex fluid components in the overflow state; establishing a wellbore dual-multi model by considering the complex flow in a wellbore in the overflow state; determining a core auxiliary equation and a boundary condition; performing grid partition and numerical discrete on a solution domain of the dual-multi model; and solving the dual-multi model to obtain the bottom hole pressure under the current measured wellhead back pressure.

3. The method according to claim 2, wherein the flow calculation parameters comprise: a wellbore structure, a drilling tool assembly, stratum data, a gas-liquid-solid phase displacement monitored on a drilling platform, a drilling fluid density, a drilling fluid viscosity, the wellhead back pressure, a wellhead temperature and pressure and the current drilling depth of a drill bit; the complex fluid components comprise: drilling fluid, inflow crude oil, stratum water, broken rock debris, hydrate, hydrocarbon gas, CO.sub.2 and H.sub.2S when a hydrate layer is drilled through.

4. The method according to claim 2, wherein the wellbore dual-multi model comprises: a gas phase continuity equation, a liquid phase continuity equation, a solid phase continuity equation, a supercritical phase continuity equation, a momentum equation,_ and an energy equation.

5. Equipment for optimizing hydraulic parameters of deepwater managed pressure drilling in real time, the equipment comprising: a meter configured to measure overflow parameters and a wellhead back pressure in the current drilling process in real time; and a controller configured to perform preprocessing and feature extraction on the overflow parameters, inputting the overflow parameters after preprocessing and feature extraction into trained support vector machine identification models for overflow judgment, when it is judged that overflow occurs at the current drilling depth, perform the following operations: reducing the opening of a throttle valve on a throttle pipeline, increasing a wellhead back pressure and increase a displacement of a submarine pump and a displacement of drilling fluid, calculating a bottom hole pressure according to a collected wellhead back pressure, judging whether overflow continues to occur in case that the calculated bottom hole pressure does not fall into a safety window; and under the condition that overflow continues to occur, mixing high-density drilling fluid with the original drilling fluid, pumping the mixture into a wellbore annulus from a drill pipe, and performing the above operations of reducing the opening of the throttle valve, increasing the displacement of the submarine pump, calculating the bottom hole pressure and judging whether overflow continues to occur until overflow no longer occurs, wherein, the trained support vector machine identification models comprise: a flow identification model, a mud pit increment identification model and a standpipe pressure identification model; the meter is configured to measure a flow differential of an inlet and an outlet, a mud pit increment and a standpipe pressure in the current drilling process in real time; and the controller is configured to perform preprocessing and feature extraction on the flow differential, the mud pit increment and the standpipe pressure, input the flow differential, the mud pit increment and the standpipe pressure after preprocessing and feature extraction into corresponding support vector machine identification models for overflow judgment to obtain an overflow probability under each identification model, and process the overflow probability under each identification model by an information fusion model to judge whether overflow occurs at the current drilling well depth.

6. The equipment according to claim 5, wherein the operation of calculating the bottom hole pressure according to the measured wellhead back pressure comprises: determining flow calculation parameters after overflow of the managed pressure drilling; determining complex fluid components in the overflow state; establishing a wellbore dual-multi model by considering the complex flow in a wellbore in the overflow state; determining a core auxiliary equation and a boundary condition; performing grid partition and numerical discrete on a solution domain of the dual-multi model; and solving the dual-multi model to obtain the bottom hole pressure under the current measured wellhead back pressure.

7. The equipment according to claim 6, wherein the flow calculation parameters comprise: a wellbore structure, a drilling tool assembly, stratum data, a gas-liquid-solid phase displacement monitored on a drilling platform, a drilling fluid density, a drilling fluid viscosity, the wellhead back pressure, a wellhead temperature and pressure and the current drilling depth of a drill bit; the complex fluid components comprise: drilling fluid, inflow crude oil, stratum water, broken rock debris, hydrate, hydrocarbon gas, CO.sub.2 and H.sub.2S when a hydrate layer is drilled through.

8. The equipment according to claim 6, wherein the wellbore dual-multi model comprises: a gas phase continuity equation, a liquid phase continuity equation, a solid phase continuity equation, a supercritical phase continuity equation, a momentum equation, and an energy equation.

Description

BRIEF DESCRIPTION OF DRAWINGS

(1) The accompanying drawings are intended to provide further understanding of the embodiment of the present invention, constitute a part of the specification, and together with the following specific embodiments, are used to explain the embodiments of the present invention, but do not constitute a limitation to the embodiments of the present invention. In the accompanying drawings:

(2) FIG. 1 is a flowchart of real-time optimization of hydraulic parameters of deepwater managed pressure drilling based on a dual-multi model and big data fusion;

(3) FIG. 2 is a flowchart of solving a dual-multi model to obtain a bottom hole pressure according to an embodiment of the present invention;

(4) FIG. 3 is a flowchart of optimization of hydraulic parameters based on a dual-multi model in an overflow state; and

(5) FIG. 4 is a structural schematic diagram of equipment for optimizing hydraulic parameters in real time based on a dual-multi model and big data fusion according to an embodiment of the present invention.

DETAILED DESCRIPTION OF THE EMBODIMENTS

(6) The specific implementations of embodiments of the present invention are described below in detail with reference to the accompanying drawings. It should be understood that the specific embodiments described herein are only used to illustrate and interpret the present invention and are not intended to limit the embodiments of the present invention.

(7) The marine managed pressure drilling technology can meet the requirements of exploration and development of the natural gas field under the complex marine drilling environment. The existing domestic and foreign marine managed pressure drilling technology applications, which are mainly focused on the double-gradient drilling and control mud cap drilling technologies and mainly aim at single-phase flow of the drilling fluid in the wellbore and gas-liquid two-phase flow under the gas injection working condition, have high dependence degree on data of well logging during drilling. Meanwhile, the overflow monitoring method adopted in the drilling site is mainly focused on a threshold method, so the false alarm rate is high. In the aspect of the managed pressure drilling theoretical technology, the existing land managed pressure drilling technology considers the dissolution and precipitation of the acid gas with high CO.sub.2 and H.sub.2S content in the drilling fluid, but ignores the influence of the phase change of the acid gas in the wellbore and the formation of the natural gas hydrate in the high-pressure and low-temperature environment near the seabed mud line on the pressure of the wellbore. Therefore, it is of great significance to realize the early monitoring of overflow in the wellbore by the big data fusion analysis method and calculate the pressure of the wellbore in real time by the wellbore multi-component and multi-phase flow model aiming at the found overflow working condition so as to realize accurate managed pressure drilling and timely discover and handle the abnormal condition.

(8) In view of the problems of the managed pressure drilling theory and the early monitoring of overflow when the deepwater natural gas field is drilled to the reservoir and the hydrate reservoir is penetrated in the drilling process, the present invention provides a method for optimizing hydraulic parameters of deepwater managed pressure drilling in real time based on a dual-multi model and big data fusion. The method specifically comprises:

(9) 1. Constructing Three Pressure Prediction Profiles of a Deepwater Seabed Stratum

(10) The three pressure prediction profiles of the deepwater seabed stratum are constructed according to logging information and adjoining well data before drilling on the platform.

(11) 2. Early Monitoring Overflow of Deepwater Drilling Based on Big Data Fusion

(12) A database is formed by historical drilling data of the current development block and deepwater drilling overflow data in the existing literature, data of the database are subjected to preprocessing and feature extraction and are trained by a support vector machine, a kernel function in the support vector machine is optimized by a particle swarm algorithm at the same time to obtain optimal trained support vector machine models, overflow parameters in the current drilling process are acquired in real time and are input into the trained support vector machine identification models for overflow judgment after preprocessing and feature extraction to obtain an overflow probability under each identification model, finally whether overflow occurs at the current drilling well depth is judged by an information fusion model, overflow handling is performed if overflow is monitored, and normal drilling is performed if overflow is not monitored.

(13) 3. Maintaining Safe Managed Pressure Drilling when Overflow is not Monitored

(14) When overflow at the current drilling depth is not monitored through an overflow risk judgment method, the managed pressure drilling shall be continued in combination with the three pressure profiles of stratum.

(15) 4. Performing Real-Time Simulated Calculation on a Bottom Hole Pressure when Overflow is Monitored

(16) When overflow is monitored to occur at the current drilling depth, the complex flow state in the wellbore is analyzed, a dual-multi model of the wellbore (that is, an eight-component four-phase flow control equation set) is established, and the bottom hole pressure under the current wellhead back pressure is calculated. The step of predicting the bottom hole pressure in real time by the dual-multi model is as follows:

(17) (1) flow calculation parameters after overflow of the managed pressure drilling are determined. The calculation parameters mainly comprise: a wellbore structure, a drilling tool assembly, stratum data, a gas-liquid-solid phase displacement monitored on a drilling platform, a drilling fluid density, a drilling fluid viscosity, a real-time wellhead back pressure, a wellhead temperature and pressure and the current drilling depth of a drill bit.

(18) (2) Complex fluid components in the overflow state are determined. The complex fluid components are focused on eight components, specifically comprising: drilling fluid, inflow crude oil, stratum water, broken rock debris, hydrate, hydrocarbon gas, CO.sub.2 and H.sub.2S when a hydrate layer is drilled through.

(19) (3) A wellbore dual-multi model is established by considering the complex flow in the wellbore in the overflow state. The multiple phases in the wellbore are mainly focused on a gas phase, a liquid phase, a solid phase and a supercritical phase, and the dual-multi model comprises continuity equations of a gas phase (hydrocarbon gas, CO.sub.2 and H.sub.2S invaded in the stratum), a liquid phase (drilling fluid, produced stratum water and crude oil), a solid phase (rock debris and a hydrate phase) and a supercritical phase, a momentum equation and an energy equation.

(20) (4) A core auxiliary equation and a boundary condition are determined. To realize accurate solution of the multi-component and multi-phase flow control equation set established in (3), it is necessary to establish a certain calculation auxiliary equation and determine an initial boundary condition, wherein the core auxiliary equation comprises: a hydrate formation and decomposition equation, a solubility calculation equation of hydrocarbon gases (CH.sub.4, C.sub.2H.sub.6, C.sub.3H.sub.8, etc.) and acid gases (CO.sub.2, H.sub.2S), a supercritical phase discrimination equation, a stratum hydrocarbon gas production equation, a stratum acid gas production equation, etc.

(21) (5) A solution domain of the dual-multi model is subjected to grid partition and numerical discrete. The established dual-multi model is subjected to grid partition of a time domain and a space domain in the wellbore to determine a time step length and a space step length. Meanwhile, the continuity equation, the momentum equation and the energy equation in the dual-multi model are subjected to numerical discrete by a four-point finite difference method.

(22) (6) The bottom hole pressure under the current measured wellhead back pressure is obtained by solving the dual-multi model. According to the obtained three pressure prediction profiles of the stratum, an initial value of the bottom hole pressure is assumed in real time, the dual-multi model is solved to obtain the bottom hole pressure under the current wellhead back pressure value and obtain a multi-phase flow parameter in the wellbore. The multi-phase flow parameter comprises: temperature and pressure distribution in a wellbore annulus and volume fraction of each phase and each component.

(23) 5. Optimizing the hydraulic parameters by combining the three pressure prediction profiles of the stratum and the real-time simulated calculation by the dual-multi model. The safe drilling pressure window at the current well depth is obtained through three pressure prediction profiles of the stratum. The marine drilling platform rapidly adjusts the density of the drilling fluid, the opening of the throttle valve and the displacement of the submarine pump and accurately controls the bottom hole pressure to ensure that the bottom hole pressure is within an appropriate safe window.

(24) 6. Performing Real-Time Overflow Monitoring and Real-Time Calculation by the Dual-Multi Model to Realize Safe and Efficient Drilling

(25) In combination with the early monitoring of overflow and calculation on the bottom hole pressure by the dual-multi model when overflow occurs, the hydraulic parameters are adjusted in real time according to the three pressure profiles of the stratum (an appropriate safe pressure window can be determined according to the three pressure profiles of the stratum), and safe and efficient managed pressure drilling is maintained.

(26) FIG. 1 is a calculation flowchart of a method for optimizing hydraulic parameters of marine managed pressure drilling in real time based on a multi-component and multi-phase flow model. The main implementation steps are as follows:

(27) 1. Construction of Three Pressure Prediction Profiles of a Deepwater Seabed Stratum

(28) The three pressure prediction profiles of the deepwater seabed stratum are constructed according to logging information and adjoining well data before drilling on the platform.

(29) 2. Early Monitoring of Overflow of Deepwater Drilling Based on Big Data Fusion

(30) (1) A database is formed by historical drilling data of the current development block and deepwater drilling overflow data in the existing literature, data of the database are subjected to preprocessing and feature extraction, and an error penalty factor and a nuclear parameter in a support vector machine (SVM) are optimized by a particle swarm algorithm to obtain the optimal trained support vector machine overflow identification modules (a flow identification model, a mud pit increment identification model and a standpipe identification model);

(31) (2) the overflow parameters (a flow differential of an inlet and an outlet, a mud pit increment and a standpipe pressure) in the current drilling process are acquired in real time and are input into the trained support vector machine (SVM) identification models after preprocessing and feature extraction for overflow judgment to obtain an overflow probability under each identification model; and

(32) (3) whether overflow occurs at the current drilling well depth is judged by an information fusion model, overflow handling is performed if overflow is monitored, and normal drilling is performed if overflow is not monitored.

(33) Preprocessing in (1) adopts Fourie transform filtering and noise reduction processing, abnormal points with large fluctuation are removed, and monitoring parameters with small fluctuation are subjected to smoothing processing by a mean filtering method:

(34) $f_n=\frac{1}{N}\underset{k=0}{\overset{N-1}{.Math.}}F\left(k\right)W_N^{-\mathrm{kn}}$
herein, f.sub.n is data after filter change; F(k) is a finite length sequence with a length M; N is an interval length of Fourie transform, wherein N>M; k=0, 1, 2, . . . N−1 is a frequency variable; and W.sub.N is a rotation factor and n is a time variable.

(35) The feature extraction of the obtained data in (1) mainly aims at the representation of variation of each overflow monitoring parameter within a certain time, optimization of the support vector machine (SVM) by the particle swarm algorithm mainly aims at the error penalty factor C and the nuclear parameter σ.sup.2, and a fitness function of the support vector machine is as follows:

(36) $f_{\mathrm{fit}}=\frac{1}{n}\underset{i=1}{\overset{n}{.Math.}}.Math.\overline{y_i}-y_i.Math.$
wherein n is sample capacity, y.sub.i is training set output, y.sub.i is optimization output, and after specified iterations are reached, optimization output of the optimal parameter is stopped to obtain the optimal support vector machine model.

(37) The information fusion model for overflow judgment in (3) is mainly focused on D-S multi-source information. Firstly, according to the overflow probability under each identification model obtained in (2), a normalization constant is calculated:

(38) $K=1-\underset{\mathrm{flow}.Math.\mathrm{no}\mathrm{flow}}{.Math.}\underset{i=1}{\overset{2}{.Math.}}{m}_i\left(A_i\right)$
wherein K is normalization constant; m.sub.i(A.sub.i) is occurrence probability of overflow or non-overflow; in A.sub.i, i=1 represents overflow event, and i=2 represents non-overflow event; flow is overflow state; no flow is non-overflow state.

(39) The obtained overflow occurrence probability M (flow) is:

(40) $\begin{array}{cc}M\left(\mathrm{flow}\right)=\frac{1}{K}\underset{\mathrm{flow}.Math.\mathrm{no}\mathrm{flow}=\left\{\mathrm{flow}\right\}}{.Math.}\underset{i=1}{\overset{2}{.Math.}}{m}_i\left(A_i\right)& \end{array}$
if the probability after fusion is higher than a certain threshold (for example, 0.5), it shows that overflow occurs in the managed pressure drilling process; and if the probability after fusion is lower than the threshold (for example, 0.5), normal managed pressure drilling is performed and no overflow occurs.

(41) 3. Real-Time Simulated Calculation of the Bottom Hole Pressure when Overflow is Monitored

(42) (1) Determination of flow calculation parameters during overflow Multi-phase flow calculation parameters of the deepwater managed pressure drilling are obtained, wherein the calculation parameters mainly comprise: a wellbore structure, a drilling tool assembly, stratum data, a gas-liquid-solid phase displacement in drilling, physical data of drilling fluid, real-time wellhead back pressure, temperature and pressure at the seabed mud line wellhead and the current drilling depth of the drill bit.

(43) (2) Determination of Complex Fluid Components in the Overflow State

(44) When overflow occurs, the fluid in the wellbore are focused on 8 components, specifically comprising: drilling fluid, inflow crude oil, stratum water, broken rock debris, hydrate, hydrocarbon gas, CO.sub.2 and H.sub.2S when a hydrate layer is drilled through.

(45) (3) Establishment of a Dual-Multi Model in the Wellbore by Considering the Complex Fluid Flow in the Wellbore in the Overflow State

(46) For the complex situation of the fluid components and the flow state in the marine managed pressure drilling process, a dual-multi model in the wellbore is established. The “dual-multi” in the model refers to eight-component and four-phase flow, specifically comprising: a gas phase (hydrocarbon gas, CO.sub.2 and H.sub.2S invaded in the stratum), a liquid phase (drilling fluid and produced stratum water), a solid phase (rock debris and a hydrate phase) and a supercritical phase. The dual-multi model contains continuity equations of the four phases, a total momentum equation and an energy equation. Parameters required by calculation comprise: the densities ρ.sub.m, ρ.sub.w, ρ.sub.c, ρ.sub.g, ρ.sub.CO2, ρ.sub.H, ρ.sub.SC, ρ.sub.o and ρ.sub.H2S (kg/m.sup.3) of the drilling fluid, stratum water, rock debris, stratum hydrocarbon gas, CO.sub.2, hydrate phase, supercritical phase, crude oil and H.sub.2S at the local temperature and pressure; the local up-hole velocities v.sub.m, v.sub.w, v.sub.c, v.sub.g, v.sub.CO2, v.sub.H, v.sub.SC, v.sub.o and v.sub.H2S (m/s) of the drilling fluid, stratum water, rock debris, stratum hydrocarbon gas, CO.sub.2, hydrate phase, supercritical phase, crude oil and H.sub.2S; the local volume fractions E.sub.m, E.sub.w, E.sub.c, E.sub.g, E.sub.CO2, E.sub.H, E.sub.SC, E.sub.O and E.sub.H2S (non-dimensional) of the drilling fluid, stratum water, rock debris, stratum hydrocarbon gas, CO.sub.2, hydrate phase, supercritical phase, crude oil and H.sub.2S; the sectional area A (m.sup.2) of the annulus; the mass q.sub.g (kg/s.Math.m.sup.3) of natural gas (CH.sub.4, C.sub.2H.sub.6 and C.sub.3H.sub.8) produced per unit time and per unit thickness; the mass fraction x.sub.g (zero dimension) of the natural gas in the hydrate; the formation/decomposition rate of r.sub.H (kg/s.Math.m) of the hydrate of the natural gas per unit length in the wellbore; the formation/decomposition rate r.sub.sc (kg/s.Math.m) of the supercritical phase per unit length in the wellbore; the mass q.sub.c, q.sub.w, q.sub.CO2, q.sub.H2S, q.sub.o, q.sub.SC (kg/s) of the rock debris, produced water, CO.sub.2, H.sub.2S, crude oil and supercritical phase per unit thickness; the solubility R.sub.i (m.sup.3/m.sup.3) of the natural gas (i=CH.sub.4, C.sub.2H.sub.6 and C.sub.3H.sub.8) in the drilling fluid; the density ρ.sub.gi (kg/m.sup.3) of the natural gas (i=CH.sub.4, C.sub.2H.sub.6 and C.sub.3H.sub.8) in the standard state; the solubility of R.sub.i (m.sup.3/m.sup.3) of the acid gas (i=CO.sub.2 and H.sub.2S) in the drilling fluid; an angle of inclination α (°); pressure P, Pa; coordinates z and m in the flow direction; gravitational acceleration g (m/s.sup.2), and annulus friction F.sub.r (Pa); temperature T (° C.) in the annulus; stratum temperature T.sub.ei (° C.); temperature T.sub.t (° C.) in the drill pipe; the decomposition heat ΔH.sub.H (J/mol) of the hydrate phase, the average mean molecular weight M.sub.H of the hydrate phase, the decomposition heat ΔH.sub.SC (J/mol) of the supercritical phase, and the average mean molecular weight M.sub.SC (kg/mol) of the supercritical phase; the mass flow rate w (kg/s) of the fluid; the specific heat capacity C (J/kg° C.) of the fluid; the volume fraction E (non-dimensional) of the fluid; the volume fraction ρ (kg/m.sup.3) of the fluid; the temperature T.sub.a (° C.) of the annulus fluid; the total heat transfer coefficient U.sub.a (non-dimensional) of the annulus fluid and the stratum; the total heat transfer coefficient Ut (non-dimensional) of the annulus fluid and the drill pipe; the outer diameter r.sub.co (m) of the return pipeline; the inner diameter r.sub.ti (m) of the drill pipe; the heat conductivity coefficient k.sub.e (W/(m.Math.° C.)) of the stratum; the transient heat transfer function T.sub.D (non-dimensional); the critical pressure P.sub.ci(i=CO.sub.2.sub., H.sub.2.sub.S) (MPa) of the acid gases CO.sub.2 and H.sub.2S; the critical temperature T.sub.pci(i=CO.sub.2.sub., H.sub.2.sub.S) (K) of the acid gases CO.sub.2 and H.sub.2S; the contents y.sub.i(i≤C.sub.1.sub.,C.sub.2.sub.,C.sub.3.sub.) (non-dimensional) of C.sub.1, C.sub.2 and C.sub.3 gases in the gas phase; the Henry's constant H.sub.i(i=C.sub.1.sub., C.sub.2.sub., C.sub.3.sup.) (non-dimensional) of the C.sub.1, C.sub.2 and C.sub.3 gases in the liquid phase; the contents y.sub.i(i=CO.sub.2.sub.,H.sub.2.sub.S) (non-dimensional) of the CO.sub.2 and H.sub.2S gases in the gas phase; the Henry's constant H.sub.i(i=CO.sub.2.sub., H.sub.2.sub.S) (non-dimensional) of the CO.sub.2 and H.sub.2S gases in the liquid phase; the bottom hole pressure p.sub.b (MPa) at the time t; and depth h (m).

(47) 1) A Gas Phase Continuity Equation

(48) {circle around (1)} stratum hydrocarbon gas

(49) $\begin{array}{cc}\frac{\partial }{\partial t}\left(A{\rho }_g{E}_g+\frac{\underset{i=C_1,C_2,C_3}{.Math.}{R}_i{\rho }_{\mathrm{gi}}A{E}_o}{B_o}\right)+\frac{\partial }{\partial z}\left(A{\rho }_g{v}_g{E}_g+\frac{\underset{i=C_1,C_2,C_3}{.Math.}{R}_i{\rho }_{\mathrm{gi}}A{E}_o{v}_o}{B_o}\right)=q_g-x_g{r}_H& \end{array}$

(50) {circle around (2)} CO.sub.2 gas

(51) $\begin{array}{cc}\frac{\partial }{\partial t}\left(A{\rho }_{{\mathrm{CO}}_2}{E}_{{\mathrm{CO}}_2}+\frac{R_{{\mathrm{CO}}_2}{\rho }_{{\mathrm{gCO}}_2}A{E}_o}{B_o}\right)+\frac{\partial }{\partial z}\left(A{\rho }_{{\mathrm{CO}}_2}{v}_{{\mathrm{CO}}_2}{E}_{{\mathrm{CO}}_2}+\frac{R_{{\mathrm{CO}}_2}{\rho }_{{\mathrm{gCO}}_2}A{E}_o{v}_o}{B_o}\right)=q_{{\mathrm{CO}}_2}& \end{array}$

(52) {circle around (3)} H.sub.2S gas

(53) $\begin{array}{cc}\frac{\partial }{\partial t}\left(A{\rho }_{H_{2}S}{E}_{H_{2}S}+\frac{R_{H_{2}S}{\rho }_{g{H}_{2}S}A{E}_o}{B_o}\right)+\frac{\partial }{\partial z}\left(A{\rho }_{H_{2}S}{v}_{H_{2}S}{E}_{H_{2}S}+\frac{R_{H_{2}S}{\rho }_{g{H}_{2}S}A{E}_o{v}_o}{B_o}\right)=q_{H_{2}S}& \end{array}$

(54) 2) A Liquid Phase Continuity Equation

(55) {circle around (1)} drilling fluid

(56) $\frac{\partial }{\partial t}\left(A{\rho }_m{E}_m\right)+\frac{\partial }{\partial z}\left(A{\rho }_m{v}_m{E}_m\right)=-\left(1-x_g\right)r_H-r_{SC}$

(57) {circle around (2)} stratum water

(58) $\frac{\partial }{\partial t}\left(A{E}_w{\rho }_w\right)+\frac{\partial }{\partial z}\left(A{E}_w{v}_w{\rho }_w\right)=q_w$

(59) {circle around (3)} crude oil

(60) 0$\begin{array}{cc}\frac{\partial }{\partial t}\left(A{E}_o{\rho }_o-\frac{\underset{i=C_1,C_2,C_3,{\mathrm{CO}}_2,H_{2}S}{.Math.}{R}_i{\rho }_{\mathrm{gi}}A{E}_o}{B_o}\right)+\frac{\partial }{\partial z}\left(A{E}_o{\rho }_o{v}_o-A\frac{\underset{i=C_1,C_2,C_3,{\mathrm{CO}}_2,H_{2}S}{.Math.}{R}_i{\rho }_{\mathrm{gi}}A{E}_o{v}_o}{B_o}\right)=q_0& \end{array}$

(61) 3) A Solid Phase Continuity Equation

(62) {circle around (1)} rock debris

(63) $\frac{\partial }{\partial t}\left(A{\rho }_c{E}_c\right)+\frac{\partial }{\partial z}\left(A{\rho }_c{E}_c{v}_c\right)=q_c$

(64) {circle around (2)} hydrate phase

(65) $\frac{\partial }{\partial t}\left(A{\rho }_H{E}_H\right)+\frac{\partial }{\partial z}\left(A{\rho }_H{v}_H{E}_H\right)=r_H$

(66) 4) A Supercritical Phase Continuity Equation

(67) $\frac{\partial }{\partial t}\left(A{\rho }_{\mathrm{SC}}{E}_{\mathrm{SC}}\right)+\frac{\partial }{\partial z}\left(A{\rho }_{\mathrm{SC}}{v}_{\mathrm{SC}}{E}_{\mathrm{SC}}\right)=q_{\mathrm{SC}}+r_{\mathrm{SC}}$

(68) wherein the volume fraction of all the phases:
E.sub.m+E.sub.c+E.sub.w+E.sub.g+E.sub.o+E.sub.CO.sub.2+E.sub.H.sub.2.sub.S+E.sub.H+E.sub.SC=1

(69) 5) A Momentum Equation

(70) $\frac{\partial }{\partial t}\left[A\underset{i=1}{\overset{4}{.Math.}}\left(\underset{j=1}{\overset{n}{.Math.}}{E}_j{\rho }_j{v}_j\right)\right]+\frac{\partial }{\partial s}\left[A\underset{j=1}{\overset{4}{.Math.}}\left(\underset{j=1}{\overset{n}{.Math.}}{E}_j{\rho }_j{v}_j^2\right)\right]+\mathrm{Ag}\cos \alpha \left[\underset{j=1}{\overset{4}{.Math.}}\left(\underset{j=1}{\overset{n}{.Math.}}{E}_j{\rho }_j\right)\right]+\frac{d\left(Ap\right)}{dz}+\frac{d\left(A{F}_r\right)}{dz}=0$

(71) 6) An Energy Equation

(72) $\frac{\partial }{\partial t}\left[\underset{j=1}{\overset{4}{.Math.}}\left(\underset{j=1}{\overset{n}{.Math.}}A{\rho }_j{E}_j{C}_j{T}_a\right)\right]+\frac{\partial }{\partial s}\left[\underset{j=1}{\overset{4}{.Math.}}\left(\underset{j=1}{\overset{n}{.Math.}}A{\rho }_j{E}_j{w}_j{C}_j{T}_a\right)\right]=2\left[\frac{2\pi r_{co}{U}_a{k}_e}{k_e+r_{co}{U}_a{T}_D}\left(T_{\mathrm{ei}}-T\right)-2\pi r_{\mathrm{ti}}{U}_t\left(T-T_t\right)\right]-\frac{r_H\Delta H_H}{M_H}-\frac{r_{SC}\Delta H_{SC}}{M_{SC}}$

(73) (4) Determination of a Core Auxiliary Equation and a Boundary Condition

(74) 1) A Core Auxiliary Equation

(75) To accurately solve the established dual-multi model, it is necessary to combine with a hydrate formation equation, a gas solubility calculation equation and a supercritical judgment core auxiliary equation for solution, and it is also necessary to combine with a stratum fluid phase discrimination equation, a wellbore friction equation, a flow pattern judgment and a gas-liquid-solid three-phase slip equation and other models.

(76) {circle around (1)} A hydrate formation and decomposition equation:
r.sub.H=(P,T)

(77) {circle around (2)} A gas solubility prediction equation:
R.sub.i(i=C.sub.1.sub.,C.sub.2.sub.,C.sub.3.sub.)=(T,P,y.sub.i(i=C.sub.1.sub.,C.sub.2.sub.,C.sub.3.sub.),H.sub.i(i=C.sub.1.sub.,C.sub.2.sub.,C.sub.3.sub.))
R.sub.i(i=CO.sub.2.sub.,H.sub.2.sub.S)=(T,P,y.sub.i(i=CO.sub.2.sub.,H.sub.2.sub.S),H.sub.i(i=CO.sub.2.sub.,H.sub.2.sub.S))

(78) {circle around (3)} A supercritical phase judgment equation:

(79) $\hspace{1em}\{\begin{array}{c}P\ge P_{\mathrm{pci}\left(i=C{O}_2,H_{2}S\right)}\\ T\ge T_{\mathrm{pci}\left(j=C{O}_2,H_{2}S\right)}\end{array}$

(80) 2) An Initial Boundary Condition

(81) Solution of a stratum temperature field: the wellhead temperature is read through measurement, the stratum temperature, that is, T.sub.h=T.sub.o+ΔTh at the current drilling well depth h is obtained according to the temperature gradient ΔT of the stratum, and the temperature serves as the initial temperature at the time t.

(82) When no overflow occurs, normal drilling is performed:

(83) $E_o\left(z,0\right)=E_w\left(z,0\right)=E_H\left(z,0\right)=E_{H_{2}S}\left(z,0\right)=E_{C{O}_2}\left(z,0\right)=E_g\left(z,0\right)=E_{SC}\left(z,0\right)=0$$E_c\left(z,0\right)=\frac{v_{sc}\left(z,0\right)}{C_c{v}_{sl}\left(z,0\right)+v_{cr}\left(z,0\right)}$$E_m=1-E_c$

(84) wherein v.sub.sc, v.sub.sl and v.sub.cr are the drift velocity (kg/m.sup.3) of rock debris, liquid phase and rock debris settlement; C.sub.c is a velocity distribution coefficient

(85) Under the drilling overflow working condition and the drilling stoppage cyclic working condition, the initial boundary condition is:

(86) {circle around (1)} the drilling overflow working condition

(87) $\hspace{1em}\{\begin{array}{c}p\left(h,t\right)=p_b\\ q_g\left(h,t\right)=q_g\\ q_w\left(h,t\right)=q_w\\ q_c\left(h,t\right)=q_c\\ q_s\left(h,t\right)=q_s\\ q_{SC}\left(h,t\right)=q_{SC}\end{array}$

(88) {circle around (2)} the drilling stoppage cyclic working condition

(89) $\hspace{1em}\{\begin{array}{c}p\left(h,t\right)=p_b\\ q_g\left(h,t\right)=q_g\\ q_w\left(h,t\right)=q_w\\ q_c\left(h,t\right)=0\\ q_s\left(h,t\right)=q_s\\ q_{SC}\left(h,t\right)=q_{SC}\end{array}$

(90) (5) Grid Partition and Numerical Discrete of a Solution Domain of the Dual-Multi Model

(91) 1) To accurately solve the multi-component and multi-phase flow control equation sets established in 3 and 4, it is necessary to perform space domain and time domain grid partition on a definite solution domain. The space grid of the wellbore annulus is subjected to fixed step length division, wherein any grid length: Δz.sub.i=z.sub.i+1−z.sub.i. To track the front of the multi-phase flow in real time, the time step length Δt is obtained according to the relationship between the velocity v.sub.g of free gas and the space grid length Δz.sub.i at this position:

(92) 0$\Delta t=\frac{\Delta z_i}{v_g}$

(93) 2) the Multi-Component and Multi-Phase Flow Control Equation Set in the Wellbore is Subjected to Numerical Discrete

(94) The established dual-multi model (the continuity equations, the momentum equation and the energy equation) is subjected to numerical discrete by a finite difference method. According to the characteristic of the time domain and the space domain in the wellbore, a four-point difference format is adopted. By taking the continuity equation of the rock debris as an example, the four-point difference discrete equation is as follows:

(95) ${\left(A{E}_c{\rho }_c{u}_c\right)}_{j+1}^{n+1}-{\left(A{E}_c{\rho }_c{u}_c\right)}_j^{n+1}=\frac{\Delta z}{2\Delta t}\left[{\left(A{E}_c{\rho }_c\right)}_j^n+{\left(A{E}_c{\rho }_c\right)}_{j+1}^n-{\left(A{E}_c{\rho }_c\right)}_j^{n+1}-{\left({\mathrm{AE}}_c{\rho }_c\right)}_{j+1}^{n+1}\right]+\frac{\Delta z}{2}\left[q_{\mathrm{cj}}^{n+1}+q_{\mathrm{cj}+1}^{n+1}\right]$

(96) (6) Calculation of the Bottom Hole Pressure Under the Wellhead Back Pressure at the Current Time by the Dual-Multi Model

(97) The solution of the dual-multi model is as same as the existing computer solution, as shown in FIG. 2, the marine drilling platform obtains the wellhead back pressure and the calculation parameter at the time n, and the multi-phase flow parameters and the bottom hole pressure in the wellbore at the time n are obtained by solving the dual-multi model. The multi-phase flow parameters comprise temperature and pressure distribution at different positions of a riser and the stratum and the volume fraction and velocity distribution of each phase and each component. If it is necessary to predict the bottom hole pressure at the next time n+1, the calculated multi-phase flow parameters in the wellbore at the time n may serve as the initial condition of the time n+1. The multi-phase flow parameter and the bottom hole pressure at the time n+1 are obtained by solving the dual-multi model.

(98) 5. Realizing Safe Drilling of the Deepwater Managed Pressure Drilling by Combining the Three Pressure Prediction Profiles of the Stratum and the Real-Time Simulated Calculation by the Dual-Multi Model

(99) The opening of the throttle valve and the displacement of the submarine pump are adjusted in real time based on the real-time simulated calculation by the dual-multi model, and managed pressure drilling is continued by combining the method of adjusting the density of the drilling fluid in real time. The specific steps are shown in FIG. 3. When downhole overflow is monitored to occur, the opening of the throttle valve on the throttle pipeline is reduced and the wellhead back pressure is increased; meanwhile, the displacement of the submarine pump is increased and the displacement of the drilling fluid is increased. The bottom hole pressure of the wellhead back pressure at the current time is calculated by the dual-multi model. In combination with the prediction judgment of the three pressure profiles, if overflow continues, drilling fluid with higher density relative to the original drilling fluid is mixed with the original drilling fluid and the mixture is pumped into the wellbore annulus from the drill pipe for drilling; meanwhile, the bottom hole pressure is calculated in real time by the multi-component and multi-phase flow model until the bottom hole pressure falls into an appropriate pressure window, wherein the density of the mixed drilling fluid is determined by the following formula:

(100) ${\rho }_{\mathrm{mix}}=\frac{V_m{\rho }_m+V_h{\rho }_h}{\left(V_m+V_h\right)}$
wherein ρ.sub.mix is the density (g/cm.sup.3) of the mixed drilling fluid; V.sub.m is the volume (cm.sup.3) of the drilling fluid used during drilling in the mud pit; V.sub.h (cm.sup.3) is the volume of the used drilling fluid with high density; ρ.sub.m (g/cm.sup.3) is the density of the drilling fluid during drilling; and ρ.sub.h (g/cm.sup.3) is the concentration of the drilling fluid with high density.

(101) FIG. 4 is a structural schematic diagram of equipment for optimizing hydraulic parameters in real time based on a dual-multi model and big data fusion according to an embodiment of the present invention. As shown in FIG. 4, the equipment comprises: a meter configured to measure overflow parameters and a wellhead back pressure in the current drilling process in real time; and a controller configured to perform preprocessing and feature extraction on the overflow parameters, inputting the overflow parameters after preprocessing and feature extraction into trained support vector machine identification models for overflow judgment, when it is judged that overflow occurs at the current drilling depth, perform the following operations: reducing the opening of a throttle valve on a throttle pipeline, increasing the wellhead back pressure and increasing a displacement of a submarine pump and a displacement of drilling fluid, calculating a bottom hole pressure according to the collected wellhead back pressure, and judging whether overflow continues to occur in case that the calculated bottom hole pressure does not fall into a safety window; and under the condition that overflow continues to occur, mixing high-density drilling fluid with the original drilling fluid, pumping the mixture into a wellbore annulus from a drill pipe and performing the above operations of reducing the opening of the throttle valve, increasing the displacement of the submarine pump, calculating the bottom hole pressure and judging whether overflow continues to occur until overflow no longer occurs.

(102) The components, the performed operations and the relevant benefits of the equipment can be referenced to the description of the method for optimizing the hydraulic parameters of the deepwater managed pressure drilling in real time based on the dual-multi model and big data fusion, which are not elaborated herein.

(103) It should also be noted that the term “comprise”, “include”, or any other variant thereof is intended to cover a non-exclusive inclusion, such that a process, method, product, or device that includes a series of elements includes not only those elements, but also other elements not explicitly listed, or elements that are inherent to such a process, method, product, or device.

(104) Without more restrictions, an element defined by the phrase “comprising a . . . ” does not exclude the presence of another same element in a process, method, product, or device that includes the element.

(105) The above is only an embodiment of the present application and is not intended to limit the present application. For those skilled in the art, the present application may have various modifications and changes. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principle of the present application should be included within the scope of the claims of the present application.