SIMULATION METHOD AND SYSTEM OF JOINT EXPLOITATION OF NATURAL GAS HYDRATE, SHALLOW GAS AND DEEP-SEATED GAS

Abstract

Provided is a simulation method and system of joint exploitation of natural gas hydrate, shallow gas and deep-seated gas. The method includes: constructing a simulation model of joint exploitation of a shallow gas layer and a hydrate layer, including setting and meshing the shallow gas layer and the hydrate layer, constructing the simulation model of joint exploitation by setting a formation parameter and a production parameter, and solving the model to acquire productivity data; constructing a simulation model of joint exploitation of the shallow gas layer, the deep-seated gas layer and the hydrate layer, including setting and meshing the shallow gas layer, the deep-seated gas layer and the hydrate layer, constructing a simulation model of joint exploitation of the shallow gas layer, the deep-seated gas layer and the hydrate layer by setting the formation parameter and the production parameter, and solving the model to acquire productivity data.

Claims

1. A simulation method of joint exploitation of natural gas hydrate, shallow gas and deep-seated gas, comprising: Step S1, constructing a simulation model of joint exploitation of a shallow gas layer and a hydrate layer, comprising setting and meshing the shallow gas layer and the hydrate layer, setting a geological parameter, a production parameter and a well control parameter, constructing the simulation model of joint exploitation, and solving the model to acquire productivity data; Step S2, on basis of Step S1, adding a deep-seated gas layer below the hydrate layer, meshing the deep-seated gas layer, setting the geological parameter and the production parameter, constructing a simulation model of joint exploitation of the shallow gas layer, the deep-seated gas layer and the hydrate layer, and solving the simulation model of joint exploitation of the shallow gas layer, the deep-seated gas layer and the hydrate layer to acquire productivity data; wherein in Step S1, a hydrate decomposition kinetic equation is:
CH.sub.4.Math.N.sub.hH.sub.2OCH.sub.4(g)+N.sub.hH.sub.2O.sub.(1); wherein N.sub.h is a number of water molecules bound by hydrate; the hydrate layer follows a mass conservation equation and an energy conservation equation, and a system conservation relationship is expressed as: $\frac{M^{}}{t}+.Math.F^{}=q^{};$ where t is time in unit of s; is a component identifier, and in the system conservation relationship, represents a hydrate component h, a methane component m, a water component w or energy e; M.sup. is a sum of all components of in unit of kg/m.sup.3 or J/m.sup.3; F.sup. is a flowable component of in unit of kg/(m.sup.2.Math.s); q.sup. is a source and sink of in unit of kg/(m.sup.3.Math.s) or J/(m.sup.3.Math.s); the mass conservation equation of the hydrate component is: $M^h=S_H{}_H;$ where M.sup.h is a sum of a mass of the hydrate component in unit of kg/m.sup.3; S.sub.H is a saturation of a hydrate phase; .sub.H is a density of the hydrate phase in unit of kg/m.sup.3; the mass conservation equation of the methane component is: $M^m=S_A{}_A{X}_A^m+S_G{}_G{X}_G^m+S_H{}_H{X}_H^m;$ $F^m=X_A^m{F}_A+X_G^m{F}_G;$ $q^m=X_{q,A}^m{q}_A+X_{q,G}^m{q}_G;$ where M.sup.m is a sum of a mass of the methane component in unit of kg/m.sup.3; is a porosity of a reservoir; S.sub.A is a saturation of a water phase; S.sub.G is a saturation of a gas phase; .sub.A is a density of the water phase in unit of kg/m.sup.3; .sub.G is a density of the gas phase in unit of kg/m.sup.3; X.sub.A.sup.m is a ratio of the mass of the methane component to a mass of the water phase; X.sub.G.sup.m is a ratio of the mass of the methane component to a mass of the gas phase; X.sub.H.sup.m is a ratio of the mass of the methane component to a mass of the hydrate phase; F.sub.A is a mass flow of the water phase passing through per unit cross-sectional area in unit of kg/(m.sup.2.Math.s); F.sub.G is a mass flow of the gas phase passing through per unit cross-sectional area in unit of kg/(m.sup.2.Math.s); F.sup.m is a mass flow of the methane component passing through per unit cross-sectional area in unit of kg/(m.sup.2.Math.s); q.sup.m is a source and sink of the methane component in unit of kg/(m.sup.3.Math.s); q.sub.A is a source and sink of the water phase in unit of kg/(m.sup.3.Math.s); q.sub.G is a source and sink of the gas phase in unit of kg/(m.sup.3.Math.s); X.sub.q,A.sup.m is a ratio of the mass of the methane component to a mass of the source and sink of the water phase; X.sub.q,G.sup.m is a ratio of the mass of the methane component to a mass of the source and sink of the gas phase; the mass conservation equation of the water component is: $M^w=S_A{}_A{X}_A^w+S_G{}_G{X}_G^w+S_H{}_H{X}_H^w;$ $F^w=X_A^w{F}_A+X_G^w{F}_G;$ $q^w=X_{q,A}^w{q}_A+X_{q,G}^w{q}_G;$ where M.sup.w is a sum of a mass of the water component in unit of kg/m.sup.3; X.sub.A.sup.w is a ratio of the mass of the water component to the mass of the water phase; X.sub.G.sup.w is a ratio of the mass of the water component to the mass of the gas phase; X.sub.H.sup.w is a ratio of the mass of the water component to the mass of the hydrate phase; F.sup.w is a mass flow of the water component passing through per unit cross-sectional area in unit of kg/(m.sup.2.Math.s); q.sup.w is a source and sink of the water component in unit of kg/(m.sup.3.Math.s); X.sub.q,A.sup.w is a ratio of the mass of the water component to the mass of the source and sink of the water phase; X.sub.q,G.sup.w is a ratio of the mass of the water component to a mass in the source and sink of the gas phase; the energy conservation equation is: $M^e=\left(1-\right){}_R{H}_R+\underset{=A,G,H}{.Math.}{s}_{}{}_{}{H}_{}+{}_H{s}_H{H}_0;$ $F^e=-K_{c}T+\underset{=A,G}{.Math.}{H}_{}{F}_{};$ $q^e=\underset{=A,G}{.Math.}{H}_{}{q}_{};$ where M.sup.e is a sum of energy in unit of J/m.sup.3; .sub.R is a density of rock in unit of kg/m.sup.3; H.sub.R is an enthalpy of rock in unit of J/kg; s.sub. is a saturation of phase; H.sub. is an enthalpy of phase in unit of J/kg; s.sub.H is a change value of a hydrate saturation in current time step; H.sub.0 is decomposition/formation enthalpy of hydrate in unit of J/kg; F.sup.e is energy flow rate in unit of J/(m.sup.2.Math.s); K.sub.c is a comprehensive thermal conductivity of a system in unit of W/(m.Math.K); F.sub. is a mass flow of phase passing through per unit cross-sectional area in unit of kg/(m.sup.2.Math.s); q.sup.e is a source and sink of energy in unit of J/(m.sup.3.Math.s).

2. The simulation method according to claim 1, wherein in Step S2, the shallow gas layer and the deep-seated gas layer follow the mass conservation equation and the energy conservation equation, and the system conservation relationship is expressed as: $\frac{M^x}{t}+.Math.F^x=q^x;$ where x is a gas component g or energy e; M.sup.x is a sum of all components of x in unit of kg/m.sup.3 or J/m.sup.3; F.sup.x is a flowable component of x in unit of kg/(m.sup.2.Math.s); q.sup.x is a source and sink of x in unit of kg/(m.sup.3.Math.s) or J/(m.sup.3.Math.s); when a mass of the gas component is conserved, $M^g=S_G{}_G{X}_G;$ $F^g=X_G{F}_G;$ $q^g=X_{q,G}{q}_G;$ where M.sup.g is a sum of the mass of the gas component in unit of kg/m.sup.3; X.sub.G is a ratio of the mass of the gas component to the mass of the gas phase; F.sup.g is a mass flow of the gas component passing through per unit cross-sectional area in unit of kg/(m.sup.2.Math.s); q.sup.g is a source and sink of the gas component in unit of kg/(m.sup.3.Math.s); X.sub.q,G is a ratio of the mass of the gas component to the mass of the source and sink of the gas phase; when the energy is conserved, $M^e=\left(1-\right){}_R{H}_R+s_G{}_G{H}_G;$ $F^e=-K_{c}T+H_G{F}_G;$ $q^e=H_G{q}_G;$ where s.sub.G is a saturation of the gas phase; H.sub.G is an enthalpy of the gas phase in unit of J/kg; F.sub.G is the mass flow of the gas phase passing through per unit cross-sectional area in unit of kg/(m.sup.2.Math.s).

3. The simulation method according to claim 1, wherein in Steps S1 and S2, a gas seepage law is as follows: introducing a deviation factor into a gas state equation: $\mathrm{pV}=\mathrm{nZRT};$ where p is a gas absolute pressure in unit of Mpa; T is a gas absolute temperature in unit of K; V is a gas volume in unit of m.sup.3; n is an amount of gas substance in unit of mol; R is a gas constant, which is 8.31410.sup.3 Mpa/(mol.Math.K); Z is a deviation factor of a natural gas which is dimensionless; gas seepage is similar to liquid seepage, when gas is in a laminar flow state, Darcy's seepage law is used to describe a flow state, and for a homogeneous formation and in three-dimensional seepage space, a generalized Darcy's law is: $v=-\frac{K}{}p;$ a flow velocity of a gas flow in the three-dimensional space is expressed as: $v_x=-\frac{K}{}\frac{p}{x};$ $v_y=-\frac{K}{}\frac{p}{x};$ $v_z=-\frac{K}{}\left(\frac{p}{z}+g\right);$ where v is a gas seepage velocity in unit of m/s; K is a formation permeability in unit of D; is a gas viscosity in unit of mPa.Math.s; g is a gravity acceleration in unit of g/cm.sup.3; x, y and z are space coordinate axes; when the gas seepage velocity increases beyond a predetermined value, there is a nonlinear relationship between a seepage velocity and a pressure gradient, which does not satisfy the Darcy's seepage law; in horizontal direction, when there is turbulence and inertia resistance in a process of gas seepage, a nonlinear quadratic equation of motion satisfying dynamic law of the natural gas is: $\frac{\mathrm{dp}}{\mathrm{dx}}=-{\left(\frac{}{K}v+v^2\right)}_{};$ where is a characteristic parameter of a pore structure which influences the turbulence and the inertial resistance; when the seepage velocity increases to deviate from the Darcy's law, the Darcy's law becomes: $v=-\frac{K}{}\frac{\mathrm{dp}}{\mathrm{dx}};$ $\mathrm{wherein}$ $=\frac{1}{\left(1+\mathrm{Kv}/\right)};$ wherein is a turbulence correction coefficient; in process of joint exploitation of hydrate, shallow gas and deep-seated gas, calculation equations of a water phase relative permeability k.sub.rw, a gas phase relative permeability k.sub.rg and a capillary pressure P.sub.c are respectively: $k_{\mathrm{rw}}=k_{\mathrm{rw}0}\stackrel{\_}{S_w}{}^{1/2}{\left[1-{\left(1-\stackrel{\_}{S_w}{}^{1/m}\right)}^m\right]}^2;$ $k_{\mathrm{rg}}=k_{\mathrm{rg}0}\stackrel{\_}{S_g}{}^{1/2}{\left[1-\stackrel{\_}{S_{\mathrm{wh}}}{}^{1/m}\right]}^m;$ $P_c={P_{c0}\left[\stackrel{\_}{S_w}{}^{1/m}-1\right]}^{1-m};$ $\mathrm{where}$ $\stackrel{\_}{S_w}=\left(S_w-S_{\mathrm{wr}}\right)/\left(1-S_{\mathrm{wr}}-S_{\mathrm{gr}}\right);$ $\stackrel{\_}{S_{\mathrm{wh}}}=\left(S_w+S_h-S_{\mathrm{wr}}\right)/\left(1-S_{\mathrm{wr}}-S_{\mathrm{gr}}\right);$ where S.sub.wr is a bound water saturation, S.sub.gr is a residual gas saturation, k.sub.rw0 and k.sub.rg0 are endpoint values of permeability, P.sub.c0 is an endpoint value of the capillary pressure, and m is a van Genuchten parameter.

4. A simulation system of joint exploitation of natural gas hydrate, shallow gas and deep-seated gas, comprising: a simulation module of joint exploitation of natural gas hydrate and shallow gas, which is configured to set and mesh the shallow gas layer and the hydrate layer, set a geological parameter, a production parameter and a well control parameter, construct a simulation model of joint exploitation, and solve the model to acquire productivity data; a simulation module of joint exploitation of natural gas hydrate, shallow gas and deep-seated gas, which is configured to, on basis of the simulation module of joint exploitation of natural gas hydrate and shallow gas, add a deep-seated gas layer below the hydrate layer, mesh the deep-seated gas layer, set the geological parameter and the production parameter, construct a simulation model of joint exploitation of the shallow gas layer, the deep-seated gas layer and the hydrate layer, and solve the simulation model of joint exploitation of the shallow gas layer, the deep-seated gas layer and the hydrate layer to acquire productivity data.

5. A computer device, comprising a processor and a memory for storing instructions executable by the processor, wherein the processor, when executing the instructions, implements steps of the method according to claim 1.

6. The computer device according to claim 5, wherein in Step S2, the shallow gas layer and the deep-seated gas layer follow the mass conservation equation and the energy conservation equation, and the system conservation relationship is expressed as: $\frac{M^x}{t}+.Math.F^x=q^x;$ where x is a gas component g or energy e; M.sup.x is a sum of all components of x in unit of kg/m.sup.3 or J/m.sup.3; F.sup.x is a flowable component of x in unit of kg/(m.sup.2.Math.s); q.sup.x is a source and sink of x in unit of kg/(m.sup.3.Math.s) or J/(m.sup.3.Math.s); when a mass of the gas component is conserved, $M^g=S_G{}_G{X}_G;$ $F^g=X_G{F}_G;$ $q^g=X_{q,G}{q}_G;$ where M.sup.g is a sum of the mass of the gas component in unit of kg/m.sup.3; X.sub.G is a ratio of the mass of the gas component to the mass of the gas phase; F.sup.g is a mass flow of the gas component passing through per unit cross-sectional area in unit of kg/(m.sup.2.Math.s); q.sup.g is a source and sink of the gas component in unit of kg/(m.sup.3.Math.s); X.sub.q,G is a ratio of the mass of the gas component to the mass of the source and sink of the gas phase; when the energy is conserved, $M^e=\left(1-\right){}_R{H}_R+s_G{}_G{H}_G;$ $F^e=-K_{c}T+H_G{F}_G;$ $q^e=X_G{F}_G;$ where s.sub.G is a saturation of the gas phase; H.sub.G is an enthalpy of the gas phase in unit of J/kg; F.sub.G is the mass flow of the gas phase passing through per unit cross-sectional area in unit of kg/(m.sup.2.Math.s).

7. The computer device according to claim 5, wherein in Steps S1 and S2, a gas seepage law is as follows: introducing a deviation factor into a gas state equation: $\mathrm{pV}=\mathrm{nZRT};$ where p is a gas absolute pressure in unit of Mpa; T is a gas absolute temperature in unit of K; V is a gas volume in unit of m.sup.3; n is an amount of gas substance in unit of mol; R is a gas constant, which is 8.31410.sup.3 Mpa/(mol.Math.K); Z is a deviation factor of a natural gas which is dimensionless; gas seepage is similar to liquid seepage, when gas is in a laminar flow state, Darcy's seepage law is used to describe a flow state, and for a homogeneous formation and in three-dimensional seepage space, a generalized Darcy's law is: $v=-\frac{K}{}p;$ a flow velocity of a gas flow in the three-dimensional space is expressed as: $v_x=-\frac{K}{}\frac{p}{x};$ $v_y=-\frac{K}{}\frac{p}{y};$ $v_z=-\frac{K}{}\left(\frac{p}{z}+g\right);$ where v is a gas seepage velocity in unit of m/s; K is a formation permeability in unit of D; is a gas viscosity in unit of mPa s; g is a gravity acceleration in unit of g/cm.sup.3; x, y and z are space coordinate axes; when the gas seepage velocity increases beyond a predetermined value, there is a nonlinear relationship between a seepage velocity and a pressure gradient, which does not satisfy the Darcy's seepage law; in horizontal direction, when there is turbulence and inertia resistance in a process of gas seepage, a nonlinear quadratic equation of motion satisfying dynamic law of the natural gas is: $\frac{\mathrm{dp}}{\mathrm{dx}}=-\left(\frac{}{K}v+v^2\right);$ where is a characteristic parameter of a pore structure which influences the turbulence and the inertial resistance; when the seepage velocity increases to deviate from the Darcy's law, the Darcy's law becomes: $v=-\frac{K}{}\frac{\mathrm{dp}}{\mathrm{dx}};$ $\mathrm{wherein}$ $=\frac{1}{\left(1+\mathrm{Kv}/\right)};$ wherein is a turbulence correction coefficient; in process of joint exploitation of hydrate, shallow gas and deep-seated gas, calculation equations of a water phase relative permeability k.sub.rw, a gas phase relative permeability k.sub.rg and a capillary pressure P.sub.c are respectively: $k_{\mathrm{rw}}=k_{\mathrm{rw}0}{{\stackrel{\_}{S_w}}^{1/2}\left[1-{\left(1-{\stackrel{\_}{S_w}}^{1/m}\right)}^m\right]}^2;$ $k_{\mathrm{rg}}=k_{\mathrm{rg}0}{{\stackrel{\_}{S_g}}^{1/2}\left[1-{\stackrel{\_}{S_{\mathrm{wh}}}}^{1/m}\right]}^m;$ $P_c={P_{c0}\left[{\stackrel{\_}{S_w}}^{1/m}-1\right]}^{1-m};$ $\mathrm{where}$ $\stackrel{\_}{S_w}=\left(S_w-S_{\mathrm{wr}}\right)/\left(1-S_{\mathrm{wr}}-S_{\mathrm{gr}}\right);$ $\stackrel{\_}{S_{\mathrm{wh}}}=\left(S_w+S_h-S_{\mathrm{wr}}\right)/\left(1-S_{\mathrm{wr}}-S_{\mathrm{gr}}\right);$ where S.sub.wr is a bound water saturation, S.sub.gr is a residual gas saturation, k.sub.rw0 and k.sub.rg0 are endpoint values of permeability, P.sub.c0 is an endpoint value of the capillary pressure, and m is a van Genuchten parameter.

8. A non-transitory computer-readable storage medium, having computer instructions stored thereon, wherein the instructions, when executed, implements steps of the method according to claim 1.

9. The non-transitory computer-readable storage medium according to claim 8, wherein in Step S2, the shallow gas layer and the deep-seated gas layer follow the mass conservation equation and the energy conservation equation, and the system conservation relationship is expressed as: $\frac{M^x}{t}+.Math.F^x=q^x;$ where x is a gas component g or energy e; M.sup.x is a sum of all components of x in unit of kg/m.sup.3 or J/m.sup.3; F.sup.x is a flowable component of x in unit of kg/(m.sup.2.Math.s); q.sup.x is a source and sink of x in unit of kg/(m.sup.3.Math.s) or J/(m.sup.3.Math.s); when a mass of the gas component is conserved, $M^g=S_G{}_G{X}_G;$ $F^g=X_G{F}_G;$ $q^g=X_{q,G}{q}_G;$ where M.sup.g is a sum of the mass of the gas component in unit of kg/m.sup.3; X.sub.G is a ratio of the mass of the gas component to the mass of the gas phase; F.sup.g is a mass flow of the gas component passing through per unit cross-sectional area in unit of kg/(m.sup.2.Math.s); q.sup.g is a source and sink of the gas component in unit of kg/(m.sup.3.Math.s); X.sub.q,G is a ratio of the mass of the gas component to the mass of the source and sink of the gas phase; when the energy is conserved, $M^e=\left(1-\right){}_R{H}_R+s_G{}_G{H}_G;$ $F^e=-K_{c}T+H_G{F}_G;$ $q^e=H_G{q}_G;$ where s.sub.G is a saturation of the gas phase; H.sub.G is an enthalpy of the gas phase in unit of J/kg; F.sub.G is the mass flow of the gas phase passing through per unit cross-sectional area in unit of kg/(m.sup.2.Math.s).

10. The non-transitory computer-readable storage medium according to claim 8, wherein in Steps S1 and S2, a gas seepage law is as follows: introducing a deviation factor into a gas state equation: $\mathrm{pV}=\mathrm{nZRT};$ where p is a gas absolute pressure in unit of Mpa; T is a gas absolute temperature in unit of K; V is a gas volume in unit of m.sup.3; n is an amount of gas substance in unit of mol; R is a gas constant, which is 8.31410.sup.3 Mpa/(mol.Math.K); Z is a deviation factor of a natural gas which is dimensionless; gas seepage is similar to liquid seepage, when gas is in a laminar flow state, Darcy's seepage law is used to describe a flow state, and for a homogeneous formation and in three-dimensional seepage space, a generalized Darcy's law is: $v=-\frac{K}{}p;$ a flow velocity of a gas flow in the three-dimensional space is expressed as: $v_x=-\frac{K}{}\frac{p}{x};$ $v_y=-\frac{K}{}\frac{p}{y};$ $v_z=-\frac{K}{}\left(\frac{p}{z}+g\right);$ where v is a gas seepage velocity in unit of m/s; K is a formation permeability in unit of D; is a gas viscosity in unit of mPa.Math.s; g is a gravity acceleration in unit of g/cm.sup.3; x, y and z are space coordinate axes; when the gas seepage velocity increases beyond a predetermined value, there is a nonlinear relationship between a seepage velocity and a pressure gradient, which does not satisfy the Darcy's seepage law; in horizontal direction, when there is turbulence and inertia resistance in a process of gas seepage, a nonlinear quadratic equation of motion satisfying dynamic law of the natural gas is: $\frac{\mathrm{dp}}{\mathrm{dx}}=-\left(\frac{}{K}v+v^2\right);$ where is a characteristic parameter of a pore structure which influences the turbulence and the inertial resistance; when the seepage velocity increases to deviate from the Darcy's law, the Darcy's law becomes: $v=-\frac{K}{}\frac{\mathrm{dp}}{\mathrm{dx}};$ $\mathrm{wherein}$ $=\frac{1}{\left(1+\mathrm{Kv}/\right)};$ wherein is a turbulence correction coefficient; in process of joint exploitation of hydrate, shallow gas and deep-seated gas, calculation equations of a water phase relative permeability k.sub.rw, a gas phase relative permeability k.sub.rg and a capillary pressure P.sub.c are respectively: $k_{\mathrm{rw}}=k_{\mathrm{rw}0}{{\stackrel{\_}{S_w}}^{1/2}\left[1-{\left(1-{\stackrel{\_}{S_w}}^{1/m}\right)}^m\right]}^2;$ $k_{\mathrm{rg}}=k_{\mathrm{rg}0}{{\stackrel{\_}{S_g}}^{1/2}\left[1-{\stackrel{\_}{S_{\mathrm{wh}}}}^{1/m}\right]}^m;$ $P_c={P_{c0}\left[{\stackrel{\_}{S_w}}^{1/m}-1\right]}^{1-m};$ $\mathrm{where}$ $\stackrel{\_}{S_w}=\left(S_w-S_{\mathrm{wr}}\right)/\left(1-S_{\mathrm{wr}}-S_{\mathrm{gr}}\right);$ $\stackrel{\_}{S_{\mathrm{wh}}}=\left(S_w+S_h-S_{\mathrm{wr}}\right)/\left(1-S_{\mathrm{wr}}-S_{\mathrm{gr}}\right);$ where S.sub.wr is a bound water saturation, S.sub.gr is a residual gas saturation, k.sub.rw0 and k.sub.rg0 are endpoint values of permeability, P.sub.c0 is an endpoint value of the capillary pressure, and m is a van Genuchten parameter.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0049] FIG. 1 is a flow chart of a simulation method and system of joint exploitation of natural gas hydrate, shallow gas and deep-seated gas and a system thereof according to an embodiment of the present disclosure.

[0050] FIG. 2 is a schematic diagram of a simulation model of joint exploitation of a shallow gas layer and a hydrate layer constructed according to an embodiment of the present disclosure.

[0051] FIG. 3 is a graph showing the law that the cumulative gas production generated by the joint exploitation of natural gas hydrate and shallow gas and the gas production rate change with time at the bottom hole flowing pressure of 3.75 MPa according to an embodiment of the present disclosure.

[0052] FIG. 4 is a schematic diagram of a simulation model of joint exploitation of a shallow gas layer, a deep-seated gas layer and a hydrate layer constructed according to an embodiment of the present disclosure.

[0053] FIG. 5 is a graph showing the law that the cumulative gas production generated by the joint exploitation of shallow gas, deep-seated gas and natural gas hydrate and the gas production rate change with time at the bottom hole flowing pressure of 3.75 MPa according to an embodiment of the present disclosure.

DETAILED DESCRIPTION OF THE EMBODIMENTS

[0054] The technical schemes in the embodiments of the present disclosure will be clearly and completely described with reference to the drawings in the embodiments of the present disclosure hereinafter. Obviously, the described embodiments are only some embodiments of the present disclosure, rather than all of the embodiments. Based on the embodiment of the present disclosure, all other embodiments obtained by those skilled in the art without creative labor fall within the scope of protection of the present disclosure.

[0055] As shown in FIG. 1, the simulation method of joint exploitation of natural gas hydrate, shallow gas and deep-seated gas includes the following steps.

[0056] S1, constructing a simulation model of joint exploitation of a shallow gas layer and a hydrate layer, includes: setting and meshing the shallow gas layer and the hydrate layer setting a geological parameter, a production parameter, a well control parameter, etc.; constructing the simulation model of joint exploitation, and solving the model to acquire productivity data. A hydrate decomposition kinetic equation is:

CH.sub.4.Math.N.sub.hH.sub.2OCH.sub.4(g)+N.sub.hH.sub.2O.sub.(1); [0057] where N.sub.h is the number of water molecules bound by hydrate.

[0058] The hydrate layer follows a mass conservation equation and an energy conservation equation, and a system conservation relationship is expressed as:

[00015]$\frac{M^{}}{t}+.Math.F^{}=q^{};$ [0059] where t is time in unit of s; is a component identifier, and in the system conservation relationship, it represents a hydrate component h, a methane component m, a water component w or energy e; M.sup. is the sum of all components of in unit of kg/m.sup.3 or J/m.sup.3; F.sup. is a flowable component of in unit of kg/(m.sup.2.Math.s); q.sup. is the source and sink of in unit of kg/(m.sup.3.Math.s) or J/(m.sup.3.Math.s).

[0060] The mass conservation equation of the hydrate component is:

[00016]$M^h=S_H{}_H;$ [0061] where M.sup.h is the sum of the mass of the hydrate component in unit of kg/m.sup.3; S.sub.H is the saturation of a hydrate phase; .sub.H is the density of the hydrate phase in unit of kg/m.sup.3.

[0062] The mass conservation equation of the methane component is:

[00017]$M^m=S_A{}_A{X}_A^m+S_G{}_G{X}_G^m+S_H{}_H{X}_H^m;$$F^m=X_A^m{F}_A+X_G^m{F}_G;$$q^m=X_{q,A}^m{q}_A+X_{q,G}^m{q}_G;$ [0063] where M.sup.m is the sum of the mass of the methane component in unit of kg/m.sup.3; is a porosity of a reservoir; S.sub.A is the saturation of the water phase; S.sub.G is the saturation of the gas phase; .sub.A is the density of the water phase in unit of kg/m.sup.3; .sub.G is the density of the gas phase in unit of kg/m.sup.3; X.sub.A.sup.m is the ratio of the mass of the methane component to the mass of the water phase; X.sub.G.sup.m is the ratio of the mass of the methane component to the mass of the gas phase; X.sub.H.sup.m is the ratio of the mass of the methane component to the mass of the hydrate phase; F.sub.A is a mass flow of the water phase passing through per unit cross-sectional area in unit of kg/(m.sup.2.Math.s); F.sub.G is the mass flow of the gas phase passing through per unit cross-sectional area in unit of kg/(m.sup.2.Math.s); F.sup.m is the mass flow of the methane component passing through per unit cross-sectional area in unit of kg/(m.sup.2.Math.s); q.sup.m is a source and sink of the methane component in unit of kg/(m.sup.3.Math.s); q.sub.A is the source and sink of the water phase in unit of kg/(m.sup.3.Math.s); q.sub.G is the source and sink of the gas phase in unit of kg/(m.sup.3.Math.s); X.sub.q,A.sup.m is the ratio of the mass of the methane component to the mass of the source and sink of the water phase; X.sub.q,G.sup.m is the ratio of the mass of the methane component to the mass of the source and sink of the gas phase.

[0064] The mass conservation equation of the water component is:

[00018]$M^w=S_A{}_A{X}_A^w+S_G{}_G{X}_G^w+S_H{}_H{X}_H^w;$$F^w=X_A^w{F}_A+X_G^w{F}_G;$$q^w=X_{q,A}^w{q}_A+X_{q,G}^w{q}_G;$ [0065] where M.sup.w is the sum of the mass of the water component in unit of kg/m.sup.3; X.sub.A.sup.w is the ratio of the mass of the water component to the mass of the water phase; X.sub.G.sup.w is the ratio of the mass of the water component to the mass of the gas phase; X.sub.H.sup.w is the ratio of the mass of the water component to the mass of the hydrate phase; F.sup.w is the mass flow of the water component passing through per unit cross-sectional area in unit of kg/(m.sup.2.Math.s); q.sup.w is the source and sink of the water component in unit of kg/(m.sup.3.Math.s); X.sub.q,A.sup.w, is the ratio of the mass of the water component to the mass of the source and sink of the water phase; X.sub.q,G.sup.w is the ratio of the mass of the water component to the mass in the source and sink of the gas phase.

[0066] The energy conservation equation is:

[00019]$M^e=\left(1-\right){}_R{H}_R+\underset{=A,G,H}{.Math.}{s}_{}{}_{}{H}_{}+{}_H{s}_H{H}_0;$$F^e=-K_{c}T+\underset{=A,G}{.Math.}{H}_{}{F}_{};$$q^e=\underset{=A,G}{.Math.}{H}_{}{q}_{};$ [0067] where M.sup.e is the sum of energy in unit of J/m.sup.3; .sub.R is the density of rock in unit of kg/m.sup.3; H.sub.R is the enthalpy of rock in unit of J/kg; s.sub. is the saturation of phase; H.sub. is the enthalpy of phase in unit of J/kg; s.sub.H is the change value of the hydrate saturation in the current time step; H.sub.0 is the decomposition/formation enthalpy of hydrate in unit of J/kg; F.sup.e is the energy flow rate in unit of J/(m.sup.2.Math.s); K.sub.c is a comprehensive thermal conductivity of the system in unit of W/(m.Math.K); F.sub. is the mass flow of phase passing through per unit cross-sectional area in unit of kg/(m.sup.2.Math.s); q.sup.e is the source and sink of energy in unit of J/(m.sup.3.Math.s).

[0068] S2, on the basis of Step S1, a deep-seated gas layer is added below the hydrate layer, the deep-seated gas layer is meshed, the geological parameter, and the production parameter, etc. are set, a simulation model of joint exploitation of the shallow gas layer, the deep-seated gas layer and the hydrate layer is constructed, and the simulation model of joint exploitation of the shallow gas layer, the deep-seated gas layer and the hydrate layer is solved to acquire productivity data.

[0069] The shallow gas layer and the deep-seated gas layer follow the mass conservation equation and the energy conservation equation, and the system conservation relationship is expressed as:

[00020]$\frac{M^x}{t}+.Math.F^x=q^x;$ [0070] where x is the gas component g or energy e; M.sup.x is the sum of all components of x in unit of kg/m.sup.3 or J/m.sup.3; F.sup.x is the flowable component of x in unit of kg/(m.sup.2.Math.s); q.sup.x is the source and sink of x in unit of kg/(m.sup.3.Math.s) or J/(m.sup.3.Math.s).

[0071] When the mass of the gas component is conserved,

[00021]$M^g=S_G{}_G{X}_G;$$F^g=X_G{F}_G;$$q^g=X_{q,G}{q}_G;$ [0072] where M.sup.g is the sum of the mass of the gas component in unit of kg/m.sup.3; X.sub.G is the ratio of the mass of the gas component to the mass of the gas phase; F.sup.g is the mass flow of the gas component passing through per unit cross-sectional area in unit of kg/(m.sup.2.Math.s); q.sup.g is the source and sink of the gas component in unit of kg/(m.sup.3.Math.s); X.sub.q,G is the ratio of the mass of the gas component to the mass of the source and sink of the gas phase.

[0073] When the energy is conserved:

[00022]$M^e=\left(1-\right){}_R{H}_R+s_G{}_G{H}_G;$$F^e=-K_{c}T+H_G{F}_G;$$q^e=H_G{q}_G;$ [0074] where s.sub.G is the saturation of the gas phase; H.sub.G is the enthalpy of the gas phase in unit of J/kg; F.sub.G is the mass flow of the gas phase passing through per unit cross-sectional area in unit of kg/(m.sup.2.Math.s).

[0075] The temperature of the deep-seated gas layer is low, and a compressibility of an actual gas is quite different from that of an ideal gas, so a deviation factor is introduced into a gas state equation:

[00023]$\mathrm{pV}=\mathrm{nZRT};$ [0076] where p is a gas absolute pressure in unit of Mpa; T is a gas absolute temperature in unit of K; V is a gas volume in unit of m.sup.3; n is an amount of gas substance in unit of mol; R is a gas constant, which is 8.31410.sup.3 Mpa/(mol.Math.K); Z is a deviation factor of gas which is dimensionless.

[0077] Gas seepage is similar to liquid seepage, when gas is in a laminar state, Darcy's seepage law is used to describe the flow state, and for a homogeneous formation and in the three-dimensional seepage space, a generalized Darcy's law is:

[00024]$v=-\frac{K}{}p.$

[0078] The flow velocity of gas flow in the three-dimensional space is expressed as:

[00025]$v_x=-\frac{K}{}\frac{p}{x};$$v_y=-\frac{K}{}\frac{p}{y};$$v_z=-\frac{K}{}\left(\frac{p}{z}+g\right);$ [0079] where v is a gas seepage velocity in unit of m/s; K is a formation permeability in unit of D; is the gas viscosity in unit of mPa.Math.s; g is the gravity acceleration in unit of g/cm.sup.3; x, y and z are space coordinate axes.

[0080] When the gas seepage velocity increases to a certain extent, influence of a turbulence and an inertia will become more and more obvious, there is a nonlinear relationship between a seepage velocity and a pressure gradient, which does not satisfy the Darcy's seepage law; in horizontal direction, when there is turbulence and inertia resistance in the process of gas seepage, a nonlinear quadratic equation of motion satisfying dynamic law of the natural gas is:

[00026]$\frac{\mathrm{dp}}{\mathrm{dx}}=-\left(\frac{}{K}v+v^2\right);$ [0081] where is a characteristic parameter of a pore structure which influences the turbulence and the inertial resistance; [0082] a first term at the right end of the above equation is viscosity resistance which is proportional to the seepage velocity, and a second term is inertia resistance which is proportional to the square of the seepage velocity, thus when the seepage velocity is small, the influence of the inertia resistance is not taken into account; when the seepage velocity is large, the influence of the turbulence and the inertia is gradually obvious, and when the seepage velocity increases to deviate from the Darcy's law, the inertia will play a leading role; therefore, the above equation is a generalized equation of motion, so that Darcy's law becomes:

[00027]$v=-\frac{K}{}\frac{\mathrm{dp}}{\mathrm{dx}};$$\mathrm{where}$$=\frac{1}{\left(1+\mathrm{Kv}/\right)};$ [0083] where is a turbulence correction coefficient.

[0084] In the process of joint exploitation of hydrate, shallow gas and deep-seated gas, the calculation equation of a water phase relative permeability k.sub.rw, a gas phase relative permeability k.sub.rg and the capillary pressure P.sub.c are, respectively,

[00028]$k_{\mathrm{rw}}=k_{\mathrm{rw}0}{{\stackrel{\_}{S_w}}^{1/2}\left[1-{\left(1-{\stackrel{\_}{S_w}}^{1/m}\right)}^m\right]}^2;$$k_{\mathrm{rg}}=k_{\mathrm{rg}0}{{\stackrel{\_}{S_g}}^{1/2}\left[1-{\stackrel{\_}{S_{\mathrm{wh}}}}^{1/m}\right]}^m;$$P_c={P_{c0}\left[{\stackrel{\_}{S_w}}^{1/m}-1\right]}^{1-m};$$\mathrm{where}$$\stackrel{\_}{S_w}=\left(S_w-S_{\mathrm{wr}}\right)/\left(1-S_{\mathrm{wr}}-S_{\mathrm{gr}}\right);$$\stackrel{\_}{S_{\mathrm{wh}}}=\left(S_w+S_h-S_{\mathrm{wr}}\right)/\left(1-S_{\mathrm{wr}}-S_{\mathrm{gr}}\right);$ [0085] where S.sub.wr is a bound water saturation, S.sub.gr is a residual gas saturation, k.sub.rw0 and k.sub.rg0 are endpoint values of permeability, P.sub.c0 is an endpoint value of the capillary pressure, and m is a van Genuchten parameter.

[0086] The embodiment of the present disclosure further provides a simulation system of joint exploitation of natural gas hydrate, shallow gas and deep-seated gas, including a simulation module of joint exploitation of natural gas hydrate and shallow gas and a simulation module of joint exploitation of natural gas hydrate, shallow gas and deep-seated gas.

[0087] The simulation module of joint exploitation of natural gas hydrate and shallow gas, which is configured to set and mesh the shallow gas layer and the hydrate layer, set a geological parameter, a production parameter and a well control parameter, construct the simulation model of joint exploitation, and solve the model to acquire productivity data.

[0088] The simulation module of joint exploitation of natural gas hydrate, shallow gas and deep-seated gas, which is configured to, on the basis of the simulation module of joint exploitation of natural gas hydrate and shallow gas, add a deep-seated gas layer below the hydrate layer, mesh the deep-seated gas layer, set the geological parameter and the production parameter, construct a simulation model of joint exploitation of the shallow gas layer, the deep-seated gas layer and the hydrate layer, and solve the simulation model of joint exploitation of the shallow gas layer, the deep-seated gas layer and the hydrate layer to acquire productivity data.

[0089] The embodiment of the present disclosure further provides a computer device, including a processor and a memory for storing instructions executable by the processor, where the processor, when executing the instructions, implements steps of the method described above.

[0090] The embodiment of the present disclosure further provides a non-transitory computer-readable storage medium having computer instructions stored thereon, where the instructions, when executed, implements steps of the method described above. Those skilled in the art can understand that all or part of the processes in the method of implementing the above embodiments can be completed by instructing related hardware through a computer program. The computer program can be stored in a non-transitory computer-readable storage medium. When executed, the computer program can include the processes of the embodiments of each method described above. Any reference to the memory, the storage, the database or other media used in various embodiments provided by the present disclosure may include at least one of a non-volatile memory and a volatile memory. The non-volatile memory may include a Read-Only Memory (ROM), a magnetic tape, a floppy disk, a flash memory or an optical memory, etc. The volatile memory may include a Random Access Memory (RAM) or an external cache. By way of illustration and not limitation, RAM can be in various forms, such as a Static Random Access Memory (SRAM) or a Dynamic Random Access Memory (DRAM).

[0091] The technical scheme of the present disclosure will be further described with reference to the embodiments.

[0092] As shown in FIG. 2, the embodiment of the present disclosure only takes into account a shallow gas layer and a natural gas hydrate layer. The upper capping formation and the lower capping formation of the model are impermeable. The data refer to the logging data from China Geological Survey on the pilot production of natural gas hydrate in Shenhu waters of South China Sea in 2017. The depth of the seabed level from the sea level is 1200 m, the pressure at the seabed is 11760 kPa, and the seabed temperature is 3.8 C. The density of seabed sediments is 2600 kg/m.sup.3, and the specific heat capacity of sediments is 1000 J/(kg.Math. C.). At the same time, the depth of the natural gas hydrate layer from the sea surface is set to 1530 m, and the depth of shallow gas layer from the sea surface is set to 1700 m.

[0093] The model is 560 m560 m285 m in size, 560 m in length and width, and 285 m in thickness. The model is divided into 56 grids in XY direction and 35 grids in the vertical direction. The model consists of four layers, namely, an upper capping formation and a lower capping formation, a hydrate layer and a shallow gas layer. The upper capping formation and the lower capping formation are 40 m thick and divided into four grids. The hydrate layer with a thickness of 15 m, is divided into five grids. The shallow gas layer with a thickness of 35 m, is divided into five grids. Other geological types with a thickness of 170 m between the gas hydrate layer and the shallow gas layer are set as invalid grids, and the influence of formation parameters on the results is ignored. A vertical well in the center of the model is depressurized and drilled from the top of the shallow gas layer to the bottom of the hydrate layer. The geological model divided by the formation thickness is shown in FIG. 2, in which the area between the hydrate layer and the shallow gas layer is an invalid grid area.

[0094] The hydrate layer with a thickness of 15 m, includes water, natural gas hydrate and methane gas, with an effective porosity of 33%, a hydrate saturation of 31% and a permeability of 1.5 mD. The shallow gas layer has a thickness of 35 m, a porosity of 25% and a permeability of 6 mD. The reservoir contains a methane phase and a gas phase.

[0095] FIG. 3 is a graph showing the law that the cumulative gas production generated by the joint exploitation of natural gas hydrate, shallow gas and deep-seated gas and the gas production rate change with time at the bottom hole flowing pressure of 3.75 MPa. From FIG. 3, it can be seen that the gas production rate is a curve of first rising, then rapidly decreasing, and then gently decreasing, while the cumulative gas production is a curve of first rapidly increasing and then gently increasing. Therefore, the production process can be divided into two stages. (1) At a stage of mass production of free gas, a very high gas production rate is obtained in the initial period of production and then the gas production rate decreases rapidly, and the maximum gas production rate reaches 1.8310.sup.6 m.sup.3/d. At this stage, there is a large amount of free gas in the shallow gas layer, so the gas production rate and the gas production at the initial period suddenly increase. After that, due to the mass production of free gas in the reservoir, the reservoir pressure decreases and the hydrate starts to decompose. The produced gas includes two parts: free gas and gas produced by hydrate decomposition. The produced gas at this stage is mainly free gas, which can achieve a great gas production rate at the initial stage of exploitation. However, the gas production rate drops rapidly due to the large decrease of free gas in the subsequent period, and the cumulative gas production starts to gradually become flat. (2) At a hydrate decomposition and gas production stage, the hydrate starts to decompose to produce decomposed gas, and the gas production rate is mainly influenced by the decomposition rate of the natural gas hydrate in the reservoir. Due to the large amount of exploitation of shallow gas, the gas pressure drops rapidly, thereby leading to the rapid increase of difference between the bottom hole pressure and the formation pressure, which promotes the decomposition of hydrate and accelerates the decomposition rate of hydrate. However, due to the endothermic decomposition of hydrate, the temperature of the hydrate layer decreases. Without the supplement of external heat, the decomposition rate of the hydrate gradually decreases, and the gas production rate has been declining. However, the decline rate is slow, which basically keeps a stable value of about 3.0310+m.sup.3/d, thus it has good commercial exploitation value.

[0096] As shown in FIG. 4, the embodiment of the present disclosure takes into account the shallow gas layer, the natural gas hydrate layer and the deep-seated gas layer. The upper capping formation and the lower capping formation of the model are water-saturated layers. The depth of the natural gas hydrate layer from the sea surface is set as 1530 m, the depth of the shallow gas layer from the sea surface is set as 1700 m, and the depth of the deep-seated gas layer from the sea surface is set as 4500 m.

[0097] The model is 560 m560 m3077 m in size, 560 m in length and width, and 3077 m in thickness. The model is divided into 56 grids in XY direction and 131 grids in the vertical direction. The model consists of five layers, namely, an upper capping formation and a lower capping formation, a shallow gas layer, a hydrate layer and a deep-seated gas layer. The upper capping formation and the lower capping formation are 40 m thick and divided into four grids. The hydrate layer with a thickness 15 m is divided into five grids. The shallow gas layer with a thickness of 35 m is divided into five grids. The deep-seated gas layer with a thickness 27 m is divided into three grids. Other geological types with a thickness of 155 m between the shallow gas layer and the natural gas hydrate layer and other geological types with a thickness of 2765 m between the natural gas hydrate layer and the deep-seated gas layer are set as invalid grids, and the influence of formation parameters on the results is ignored. A vertical well in the center of the model is depressurized and drilled from the top of the shallow gas layer to the bottom of the deep-seated gas layer. The geological model divided by the formation thickness is shown in FIG. 4, in which the area between the hydrate layer and the shallow gas layer and the area between the shallow gas layer and the deep-seated gas layer are invalid grid areas.

[0098] The hydrate layer with a thickness of 15 m, includes water and natural gas hydrate and methane gas, in which the effective porosity is 33%, the hydrate saturation is 31%, and the permeability is 1.5 mD. The thickness of the shallow gas layer is 35 m, the porosity is 25%, the permeability is 6 mD. The reservoir contains methane gas. The deep-seated gas layer with a thickness of 27 m, includes methane gas, in which the effective porosity is 32%, and the permeability is 7.4 mD.

[0099] FIG. 5 is a graph showing the law that the cumulative gas production generated by the multi-gas joint exploitation model and the gas production rate change with time at the bottom hole flowing pressure of 3.75 MPa. From FIG. 5, it can be seen that the gas production rate is a curve of first rising, then rapidly decreasing, and then gently decreasing, while the cumulative gas production is a curve of first rapidly increasing and then gently increasing. Therefore, the production process can be divided into two stages. (1) At a stage of mass production of free gas, a very high gas production rate is obtained in the initial period of production and then the gas production rate decreases rapidly, and the maximum gas production rate reaches 2.8210.sup.7 m.sup.3/d. At this stage, there is a large amount of free gas in the shallow gas layer and the deep-seated gas layer, so the gas production rate and the gas production at the initial period suddenly increase. The produced gas mainly comes from free gas in the reservoir and methane gas converted from hydrate. The produced gas at this stage is mainly free gas, which can achieve a great gas production rate at the early stage of exploitation. However, the gas production rate drops rapidly due to the large decrease of free gas in the subsequent period, and the cumulative gas production starts to gradually become flat. (2) At a hydrate decomposition and gas production stage, the hydrate starts to decompose to produce methane gas, and the gas production rate is mainly influenced by the decomposition rate of the natural gas hydrate in the reservoir. Due to the large amount of exploitation of shallow gas and deep-seated gas, the gas pressure drops rapidly, thereby leading to the rapid increase of the difference between the bottom hole pressure and the formation pressure, which promotes the decomposition of hydrate and accelerates the decomposition rate of hydrate. Therefore, in order to study the influence of hydrate decomposition on gas well productivity, it is necessary to accurately measure the law of gas production changing with time under different conditions. However, due to the endothermic decomposition of hydrate, the temperature of the hydrate layer decreases. Without the supplement of external heat, the decomposition rate of the hydrate gradually decreases, and the gas production rate has been declining. However, the decline rate is slow, which basically keeps a stable value of about 1.8110.sup.4 m.sup.3/d.

[0100] Of course, the above description is not a limitation of the present disclosure, and the present disclosure is not limited to the above examples. Changes, modifications, additions or substitutions made by those skilled in the art within the essential scope of the present disclosure should also belong to the scope of protection of the present disclosure.