Background
Natural gas hydrate is a solid compound with the appearance similar to ice, and is composed of low molecular weightA cage-like structure of a molecular gas (mainly hydrocarbon molecules, such as methane, ethane, and the like, and small molecular gases such as carbon dioxide, hydrogen sulfide, and the like) and water molecules under low-temperature and high-pressure conditions. Natural gas hydrates formed mainly from methane gas are dominant in nature and are generally called as combustible ice because of their appearance similar to ice. Methane hydrate is mainly stored in submarine deep water land slope environment and land permafrost region. The natural gas hydrate can be released to 164-180 m under the standard state3And methane gas of 0.87m3The water of (2). According to conservative estimation, the content of the natural gas hydrate in nature is 21 multiplied by 10m3This is almost twice the known fossil energy on earth, and is considered an ideal alternative to the fossil energy in the 21 st century.
When the concentration of the dissolved methane in the pore water in the stable zone of the hydrate in the marine environment exceeds the solubility of the methane hydrate, the dissolved methane can be crystallized to form the hydrate, and a hydrate layer trap is formed along with the increase of the content of the hydrate, and a free gas layer develops below the hydrate layer trap. The free gas under the hydrate layer under the specified conditions leaks up the channels into the sea floor and forms pits, authigenic carbonates, biocenotes, bubble plumes on the sea floor, such as the ruser hillock, blake sea-platform, etc., the north congo land slope, the norwegian open sea, and the south china sea.
The free gas under the hydrate layer is closed by the capillary force of the hydrate layer, the gas pressure of the interface of the hydrate and the free gas is increased along with the gas accumulation and the gas layer thickness increase, when the gas overpressure overcomes the defect that the capillary is closed, the gas is excited to leak, the overpressure gas pushes pore water to be discharged upwards, pit is formed on the seabed, and the pit depth reflects the fluid destruction strength and the overpressure amplitude of the free gas layer.
Disclosure of Invention
The invention aims to provide a method for estimating the depth of a pit in a marine hydrate area and evaluating a free gas layer, which is favorable for predicting the leakage and disaster of submarine methane on the one hand and evaluating the amount of submarine hydrocarbon resources on the other hand.
In order to achieve the purpose, the technical scheme of the invention is as follows:
a method for marine hydrate zone pit depth estimation and free gas layer evaluation, comprising:
acquiring a gas leakage process parameter, inputting the parameter into a fluid leakage model in the pore medium, and outputting the speed of the leakage fluid by the fluid leakage model in the pore medium;
inputting the speed of the leaked fluid into a pit depth model to obtain the pit depth;
and (4) inverting the methane concentration in the free gas layer under the hydrate through the pit depth.
Further, the fluid leakage model in the pore medium is:
in the above equation, Δ P is the total driving force of fluid transport and is the sum of the pressure drops applied to the air flow column and the water flow column (Δ P)g+ΔPw) Equal to the overpressure of the gas reservoir; ρ is the fluid density, d is the thickness of the free gas layer, μ is the fluid viscosity, V is the fluid velocity, k is the permeability of the deposit, k isrgAnd krwThe relative permeability of the deposit pore and water, hgAnd hwThe height of the air flow column and the water flow column respectively.
Further, the leakage fluid velocity is:
in the formula (2) in an approximate relationship
Knowing the fluid transport velocity as a function of the height of the gas flow column (h)
g=h-h
w) Is increased.
Further, the inputting the seepage fluid speed into the pockmark depth model to obtain the pockmark depth comprises:
fluid leakage velocity satisfaction
Forming a pit on the seabed after the quicksand sediment is removed by the ocean current, wherein the depth of the pit is determined by the bottom boundary of the quicksand sediment body; in formula (2)
By replacing the fluid velocity and by the pit depth h
pmBy replacing the height h of the water column
wObtaining a pit depth model:
in the formula (3)
The method is simplified as follows:
the density and viscosity of the fluid are constant under the condition of meeting the temperature and pressure.
Further, the inversion of the methane concentration in the hydrate lower free gas layer by pit depth comprises:
from the methane leak-off and pit depth models, pit formation is known to be related to gas overpressure, hydrate layer bottom boundary depth and deposit density:
and (3) obtaining a thickness equation of the hydrate closed underlying free gas reservoir through equation deformation, wherein the thickness of the gas layer depends on the depth of the seabed pockmark and the density of fluid in the gas reservoir:
the average fluid density in the gas reservoir is related to the fluid saturation and the overburden pressure; the calculation equation for predicting the fluid density in the gas reservoir is obtained by the deformation of the formula (6):
under a compliant environment, the density of water is assumed to be constant and the density of gas is approximated as a function of hydrostatic pressure, so the gas saturation calculation equation is expressed as:
equation (8) is used to estimate gas saturation, where the depth of the sea bottom pit and the thickness of the gas layer under the containment layer are known from ocean sonic surveys.
Compared with the prior art, the invention has the beneficial effects that:
according to the method, the seabed methane leakage condition is predicted through a fluid leakage model in a pore medium, the leakage condition is estimated according to the leakage condition, and finally the methane concentration in a hydrate underlying free gas layer is inverted through the pit depth. The method is beneficial to predicting the leakage and disaster of the submarine methane on one hand and evaluating the amount of submarine hydrocarbon resources on the other hand.
Detailed Description
Example (b):
the technical solution of the present invention is further described below with reference to the accompanying drawings and examples.
Referring to fig. 1, the method for estimating the depth of a marine hydrate pit and evaluating a free gas layer according to the embodiment mainly includes the following steps:
101. acquiring a gas leakage process parameter, inputting the parameter into a fluid leakage model in the pore medium, and outputting the speed of the leakage fluid by the fluid leakage model in the pore medium;
102. inputting the speed of the leaked fluid into a pit depth model to obtain the pit depth;
103. and (4) inverting the methane concentration in the free gas layer under the hydrate through the pit depth.
Therefore, the method can predict the submarine methane leakage condition through the fluid leakage model in the pore medium, then estimate the leakage condition according to the leakage condition, and finally invert the methane concentration in the hydrate underlying free gas layer through the pit depth. The method is beneficial to predicting the leakage and disaster of the submarine methane on one hand and evaluating the amount of submarine hydrocarbon resources on the other hand.
In the gas leakage process, the gas column and the water column both flow under the driving of the overpressure of the hydrate-enclosed free gas, and the total driving force of the fluid migration is equal to the gas overpressure (rho)w-ρg) And gd. The air flow column grows and pushes the overlying pore water in the leak path upward at the same velocity. Assuming that the hydrate stability zone is a homogeneous body, the fluid leakage rate in the leakage channel is the same, and therefore the fluid leakage model in the pore medium is:
in the above equation, Δ P is the total driving force of fluid transport and is the sum of the pressure drops applied to the air flow column and the water flow column (Δ P)g+ΔPw) Equal to the overpressure of the gas reservoir. ρ is the fluid density, d is the thickness of the free gas layer, μ is the fluid viscosity, μgAnd muwViscosity of gas and water, respectively, V fluid velocity, k permeability of the deposit, krgAnd krwThe relative permeability of the deposit pore and water, hgAnd hwThe height of the air flow column and the water flow column respectively. Since the gas saturation in the gas flow column is assumed to beThe water saturation in the water flow column is 1, and the relative permeability of water and water is 1. Fluid leak migration speed:
in equation (2) in an approximate relationship
The fluid migration velocity can be known as the height (h) of the air flow column
g=h-h
w) Is increased.
The resistance of the single length fluid in the pore medium increases with the increase of the height of the airflow column, that is, the pore fluid pressure of the sediment lattice is gradually increased, when the fluid resistance exceeds the static rock pressure of the corresponding sediment body, the corresponding sediment layer is fluidized to become flowing sand, that is, the speed of the leaked fluid must meet the requirement of the flowing sand
After the quicksand sediment is removed by sea current, a pit is formed on the seabed, and the pit depth is determined by the bottom boundary of the quicksand sediment body. In equation (2)
By replacing the fluid velocity and by the pit depth h
pmBy replacing the height h of the water column
wObtaining a pit depth model:
in equation (3)
The equation can be simplified as:
the density and viscosity of the fluid are constant under certain temperature and pressure conditions.
Through the methane leakage and pit depth models, it can be known that pit formation is related to gas overpressure, hydrate layer bottom boundary depth and deposit density.
Through equation deformation, a thickness equation of the hydrate closed underlying free gas reservoir can be obtained. The gas layer thickness depends on the depth of the subsea pit and the density of the fluid in the gas reservoir.
The average fluid density in a gas reservoir is related to the fluid saturation and the overlying pressure. A calculation equation for predicting the density of the fluid in the gas reservoir can be obtained by a variation of equation (6).
Under certain circumstances, the density of water is considered constant and the density of gas can be approximated as a function of hydrostatic pressure, and therefore the gas saturation calculation equation can be expressed as:
equation (8) can be used to estimate gas saturation, where the depth of the sea bottom pockmark and the thickness of the gas layer under the seal layer can be known from the ocean acoustic survey, and the related data is relatively common.
FIG. 2 is a schematic diagram showing the principle of free gas leakage and sea-bottom pit formation under a hydrate, wherein Z is the depth under the sea bottom, and h is the hydrate stable zone thickness (or hydrate seal depth). The procedure is as follows.
t0The free gas is trapped under the hydrate layer.
t1The gas begins to leak through the seal.
t2The height of the gas column is increased to push water flow to be discharged outwards, the height of the water flow column is correspondingly shortened, and the fluid migration speed is increased continuously.
t3At the moment, the pore pressure in the water-containing flow sediment exceeds the static rock pressure, so pits appear on the sea bottom to form a single air flow channel.
t4When the natural gas in the free gas reservoir is gradually emptied, the overpressure of the pores disappears, and the airflow column in the fluid channel gradually degrades.
t5At that time, the airflow column completely disappears, a gas chimney is left on the seabed, the hydrate sealing effect is recovered, and new gas accumulation is started.
The depth of the sea bottom pockmark can be calculated using equation (4). The unknown variables are free gas layer thickness and methane seal depth.
The methane saturation in the hydrate underlying free gas reservoir can be calculated using equation (8). Unknown variables were free gas layer thickness, pit depth and seal depth.
The above embodiments are only for illustrating the technical concept and features of the present invention, and the purpose thereof is to enable those skilled in the art to understand the contents of the present invention and implement the present invention accordingly, and not to limit the protection scope of the present invention accordingly. All equivalent changes or modifications made in accordance with the spirit of the present disclosure are intended to be covered by the scope of the present disclosure.