CN110132818B - Method and system for obtaining permeability of sediments containing natural gas hydrate - Google Patents
Method and system for obtaining permeability of sediments containing natural gas hydrate Download PDFInfo
- Publication number
- CN110132818B CN110132818B CN201910454469.5A CN201910454469A CN110132818B CN 110132818 B CN110132818 B CN 110132818B CN 201910454469 A CN201910454469 A CN 201910454469A CN 110132818 B CN110132818 B CN 110132818B
- Authority
- CN
- China
- Prior art keywords
- pore
- permeability
- hydrate
- throat
- network model
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 230000035699 permeability Effects 0.000 title claims abstract description 96
- 238000000034 method Methods 0.000 title claims abstract description 54
- 239000013049 sediment Substances 0.000 title claims abstract description 31
- NMJORVOYSJLJGU-UHFFFAOYSA-N methane clathrate Chemical compound C.C.C.C.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O NMJORVOYSJLJGU-UHFFFAOYSA-N 0.000 title claims description 23
- 239000011148 porous material Substances 0.000 claims abstract description 174
- 238000004364 calculation method Methods 0.000 claims abstract description 29
- 150000004677 hydrates Chemical class 0.000 claims abstract description 13
- 230000015572 biosynthetic process Effects 0.000 claims description 25
- 239000012530 fluid Substances 0.000 claims description 18
- 238000004590 computer program Methods 0.000 claims description 6
- 238000012512 characterization method Methods 0.000 claims description 5
- 238000004891 communication Methods 0.000 claims description 3
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Natural products C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 abstract description 16
- 239000003345 natural gas Substances 0.000 abstract description 9
- -1 natural gas hydrates Chemical class 0.000 abstract description 9
- 238000011161 development Methods 0.000 abstract description 8
- 238000004088 simulation Methods 0.000 description 12
- 239000007789 gas Substances 0.000 description 6
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 238000011160 research Methods 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 150000001875 compounds Chemical class 0.000 description 2
- 238000000354 decomposition reaction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000001257 hydrogen Substances 0.000 description 2
- 229910052739 hydrogen Inorganic materials 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- OYPRJOBELJOOCE-UHFFFAOYSA-N Calcium Chemical compound [Ca] OYPRJOBELJOOCE-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 229910052791 calcium Inorganic materials 0.000 description 1
- 239000011575 calcium Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000002002 slurry Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N13/00—Investigating surface or boundary effects, e.g. wetting power; Investigating diffusion effects; Analysing materials by determining surface, boundary, or diffusion effects
- G01N13/04—Investigating osmotic effects
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/082—Investigating permeability by forcing a fluid through a sample
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/082—Investigating permeability by forcing a fluid through a sample
- G01N15/0826—Investigating permeability by forcing a fluid through a sample and measuring fluid flow rate, i.e. permeation rate or pressure change
Landscapes
- Chemical & Material Sciences (AREA)
- Physics & Mathematics (AREA)
- Biochemistry (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Dispersion Chemistry (AREA)
- Fluid Mechanics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a method and a system for obtaining the permeability of sediments containing natural gas hydrates. The system comprises: a pore network model building module; a pore-filling hydrate model addition module; a permeability conduction equation establishing module; a permeability calculation module; the method can be made by secondary development on the basis of an open source computing tool OpenPNM of a pore network model. The invention can calculate the permeability when the pore network model contains the hydrate and the hydrate does not completely block the pore throat, provides the selection of the positions of various hydrates and solves the problem that the permeability of the deposit containing the hydrate changes along with the change of the saturation of the hydrate.
Description
Technical Field
The invention belongs to the technical field of hydrate reservoir formation and development, relates to a method for obtaining the permeability of hydrate sediments, and particularly relates to a method and a system for obtaining the permeability of sediments containing natural gas hydrates based on a numerical simulation algorithm of a pore network model.
Background
Natural gas hydrates are a class of clathrate compounds formed by hydrogen bonding water molecules to envelope small gas molecules [ Sloan E D, Koh C a. boundary hydrates of natural gas [ M ]. Third Edition. bocaaraton: CRC Press,2007], when the gas component is methane dominated natural gas, because of its flammability also known as "combustible ice", which is a promising source of energy due to its high reserves, in the last decades the major energy consumer has invested a lot in hydrate exploration and development technology research [ L i J F, Ye J L, Qin X W, et al. the first of moisture organic hydrates in South China Sea [ J ]. China geomology, 2018,1(1): 5-16. deposit, in the development process, physical properties of hydrates including permeability, porosity in moisture in pure Sea [ J ]. gel, 2018,1(1): 5-16. deposit, in the development process, physical properties of hydrates including porosity, porosity in moisture.
The hydrate occupies the space of the pore and the throat when being formed in the pore of the sediment, so that the permeability of the sediment is reduced; upon hydrate decomposition, pore space is released allowing for increased permeability of the deposit. Therefore, the permeability of the hydrate-bearing reservoir will vary with the hydrate saturation (the ratio of the hydrate volume in the sediment pores to the total volume of the sediment pores), and accurate prediction of the change in the permeability of the hydrate-bearing reservoir is a key factor in the successful development of hydrates. Because the development of indoor experiments under the condition of a thermodynamic stability domain of the hydrate is difficult, the numerical simulation is a feasible research method for the influence of the generation and decomposition of the hydrate on the permeability of the sediment.
In the prior art, a Pore Network Model (PNM) described in document [1] flow in Pore medium-Pore Model and porous Science [ J ]. Current Opinion in porosity & Interface Science,2001,6(3): 197) is a mature porous medium multiphase flow permeability simulation method, and is also widely applied to a numerical simulation study of the saturation and permeability properties of hydrates, document [2] (Dai S, Seol Y. water permeability in moisture-bearing characteristics: A Pore-slurry [ J ]. Geopolysical Research L, 2014,41(12): 4176-84, document [3] J, calcium J.C.compressive porosity and hydrate, and is applied to a hydrate formation simulation method of a hydrate, a hydrate formation simulation method, a hydrate formation method, a hydrate simulation method, a hydrate formation method, a hydrate simulation method, a hydrate formation method, a hydrate simulation method, a hydrate formation method, a hydrate simulation method, a hydrate formation method, a hydrate simulation method, a hydrate formation method, a hydrate formation method, a hydrate simulation method, a hydrate formation method, a hydrate formation method, a hydrate formation method, a hydrate formation method, a hydrate formation method, a hydrate formation method, a hydrate formation method, a hydrate formation.
Disclosure of Invention
In order to overcome the defects of the prior art, the invention provides a method and a system for acquiring the permeability of sediments containing natural gas hydrates.
For ease of illustration, the convention herein:
saturation degree of hydrate: the ratio of the hydrate volume in the sediment pores to the total volume of the sediment pores is dimensionless.
Hydrate: short for natural gas hydrate. The gas is a cage-shaped compound formed by connecting water molecules with small gas molecules in an envelope manner by hydrogen bonds, and when the gas component is natural gas mainly containing methane, the gas component is also called 'combustible ice' due to combustibility.
Pore space: non-solid fraction in the sediment.
Pore size: the larger of the independent spaces in the pore space.
A throat: a narrow portion of the pore space connecting the two pores.
And (3) node: the geometric center point of the aperture.
A conduit: the communication part between two adjacent nodes of the pore comprises: a left side node half pore, a throat between two nodes and a right side node half pore.
The steps of calculating the permeability of the deposit by the pore network model method (PNM) are as follows:
1. and establishing a pore network geometric model.
2. The permeability conductivity of the pores and throat in the network model was calculated.
3. All nodes of the network model are conserved in volume.
4. And (5) assigning the inlet and outlet pressure.
5. The pressure at each node is calculated.
6. The inlet flow and outlet flow are calculated.
7. Permeability was calculated using Darcy's (Darcy) law.
On the basis of a pore network model method, a pore hydrate model is added in a pore network geometric model, and a permeability conductivity calculation method of pores and throats in the pore network model under the condition of containing hydrate is given, so that the pore network model method is suitable for permeability calculation of hydrate-containing sediments. In the implementation process of the invention, the pore network model method adopts a pore network model open source computing tool-OpenPNM as a bottom layer, and secondary development is carried out on the basis of the pore network model open source computing tool-OpenPNM.
The technical scheme provided by the invention is as follows:
adding a pore filling type hydrate model in pore throat space of a finished pore network model under the condition that pores of sediments are filled with natural gas hydrate, establishing a permeability conduction equation under the condition, and calculating to obtain the permeability of the natural gas hydrate-containing pore network model; the method comprises the following steps:
A. establishing a pore network model, and adding a pore filling type hydrate model into a pore throat space of the pore network model;
the nodes in the established pore network model represent the central points of pores in the sediment; the conduit in the pore network geometric model represents a communication part between the central points of two adjacent pores in the sediment. The catheter comprises: the left node half pore, the throat between two nodes and the right node half pore in two adjacent pore nodes.
In specific implementation, adding a pore filling type hydrate model in a pore throat space of a pore network model is to add a pore hydrate geometric characteristic characterization array parameter in a geometric class in OpenPNM, and comprises the following steps:
A1. the geometric characteristic characterization array parameters of the added pore hydrate comprise: diameter array h of hydrate generated in pore centeriThe thickness array h of the hydrate generated on the surface of the porepThickness array h of hydrate generated on surface of throatt；
A2. Setting the pore diameter array in the pore network model as D, the throat length array as L, and the throat diameter array as T, when hydrate grows on the pore surface, the pore diameter array is D-hpWhen the hydrate grows on the surface of the throat, the diameter of the throat is recorded as an array T-ht。
B. Deducing a hydrate-containing conduit permeability conduction equation, and calculating to obtain the inter-conduit permeability of the natural gas hydrate-containing pore network model;
in specific implementation, in the calculation class in OpenPNM, the calculation method for the permeability between nodes when the hydrate is contained in the added pore includes the following operations:
B1. for the hydrate generated in the center of the pore, the conduit between the nodes is divided into 5 sections (numbered ④ - ⑧ in figure 2 b), and the form of each section of conduit is changed into a shape capable of simplifying the calculation of permeability by simplifying the geometric form of the conduit (shown as d in figure 2), wherein the sections numbered ⑤ - ⑦ in the figure are cylindrical, the sections ④ and ⑧ are annular columns, and the permeability of the annular section of the annular columns is (the permeability g is)loop) The formula of (1):
wherein D is the outside diameter of the ring section (pore diameter), hiThe diameter in the annular section (hydrate diameter) and L the length of the segment (═ 0.5 h)i) And mu is pore fluid viscosity.
Permeability (g) of cylindrical sectionc) The calculation formula is as follows:
the permeability (g) of the entire catheter is expressed as formula 3 according to the Hagen-Poiseuille equation:
wherein, gzone1～gzone55 segments into which the inter-node conduit is divided;
B2. for hydrates generated on the surfaces of pores and throats, the conduits between nodes are divided into 3 sections (① - ③ in FIG. 2 a), and the 3 sections of conduits are simplified into cylinders (① - ③ in FIG. 2c) by simplifying the conduit shape (FIG. 2 c). The formula for the permeability of the section of the cylinder is given by the formula 2, and only the diameter D in the formula 2 needs to be replaced by D-hp(①, ③ in FIG. 2c) and T-ht(② in FIG. 2c) to obtain the interductal permeability;
C. and calculating the permeability of the sample by adopting a pore network model method.
The fluid flow regime in the internodal conduit is assumed to be laminar. The fluid seepage flow Qij in the throat has the following form:
wherein p is the pressure of the fluid, rij and lij are the radius and length of the throat respectively, mu is the viscosity coefficient of the fluid, and gij is the throat backflow rate.
Assuming the fluid in the formation pore is incompressible, the sum of the flow rates of all throats communicating with it should be zero for a single pore, i.e.:
∑Qijeither case 0 (formula 5)
Therefore, a linear equation system with the pressure as an unknown number can be obtained, and the pressure distribution of the formation pore space can be obtained after the solution. And finally, calculating the inlet flow and the outlet flow of the stratum, and solving the absolute permeability of the stratum by applying Darcy's law.
And calculating the permeability of the natural gas hydrate through the steps.
In specific implementation, the calculation process can execute a calculation method carried by the tool OpenPNM, so that the permeability of the gas hydrate can be obtained.
(II) an acquisition device for the permeability of the sediments containing the natural gas hydrate, which comprises a memory and a processor; the memory for storing a computer program; the processor is used for realizing the method for acquiring the permeability of the natural gas hydrate-containing sediment when the computer program is executed.
(iii) a computer-readable storage medium having stored thereon a computer program which, when executed by a processor, implements the above-described method for obtaining permeability of a gas hydrate-containing deposit.
(IV) the system for acquiring the permeability of the sediments containing the natural gas hydrate is characterized in that the system is developed for the second time on the basis of a pore network model open source computing tool OpenPNM; the system comprises the following modules:
the pore network model establishing module is used for establishing a pore network model of the sediment;
the pore filling type hydrate model adding module is used for adding a pore filling type hydrate model in the pore throat space of the pore network model;
a permeability conduction equation establishing module, wherein the permeability comprises the permeability of an annular section of an annular column shape and the permeability of a cylindrical section;
and the permeability calculation module is used for calculating the permeability of the natural gas hydrate-containing pore network model by using a calculation method for the permeability between nodes under the condition that the pores contain hydrates.
The invention has the beneficial effects that:
the invention provides a method for acquiring the permeability of sediments containing natural gas hydrates. The invention realizes the continuous simulation of the permeability of the water-containing hydrate by adopting an improved pore network model and a simplified formula for the first time. By utilizing the solution provided by the invention, the permeability of the pore network model containing the hydrate and the pore throat which is not completely blocked by the hydrate can be calculated, the selection of the positions where various hydrates exist is provided, and the problem that the permeability of the deposit containing the hydrate changes along with the change of the saturation of the hydrate is solved.
Drawings
FIG. 1 is a diagram of the present invention adding custom hydrate geometry parameters in the geometry class in OpenPNM; the method comprises the following steps: hydrate is generated in a diameter array hi at the center of the pore; hp is a thickness array of the hydrate generated on the surface of the pore; ht is the array of the thickness of the hydrate generated on the surface of the throat.
FIG. 2 is an inter-node catheter computation model established by the present invention;
the device comprises a conduit 3 section, a conduit 5 section, a conduit between nodes, a hydrate section, a throat between nodes, a cylindrical section between nodes and a cylindrical section between nodes, wherein (a) is a conduit model between nodes when the hydrate grows on the surface of a pore, (b) is a conduit model between nodes when the hydrate grows in the center of the pore, (c) is a simplified calculation model of the permeability of the conduit between nodes when the hydrate grows on the surface of the pore, (d) is a calculation model of the permeability of the conduit between nodes when the hydrate grows in the center of the pore, ① - ③ respectively represent a left node semi-pore section, a throat between nodes section and a right node semi-pore section in the conduit 3 section when the hydrate grows on the surface of the pore, and ④ - ⑧ respectively represent a left node cylindrical section, a right node cylindrical.
Fig. 3 is a block diagram of an implementation of the method of the present invention.
Detailed Description
The invention will be further described by way of examples, without in any way limiting the scope of the invention, with reference to the accompanying drawings.
The invention provides a method and a system for acquiring the permeability of sediments containing natural gas hydrates.
FIG. 2 is a calculation model of an inter-node conduit, wherein (a) is an inter-node conduit model when a hydrate grows on a pore surface, (b) is an inter-node conduit model when the hydrate grows at a pore center, (c) is a simplified calculation model of inter-node conduit permeability when the hydrate grows at the pore surface, (d) is an inter-node conduit permeability calculation model when the hydrate grows at the pore center, ① - ③ respectively represent a left-node semi-pore section, an inter-node throat section and a right-node semi-pore section in a conduit 3 section when the hydrate grows at the pore surface, and ④ - ⑧ respectively represent a left-node annular cylindrical section, a left-node cylindrical section, an inter-node throat section, a right-node cylindrical section and a right-node annular cylindrical section in an inter-node conduit 5 section when the hydrate grows at the pore center.
FIG. 3 shows an implementation flow of the method of the present invention, which includes the following steps:
A. adding a pore hydrate geometric feature characterization array in a geometric class in an open source computing tool-OpenPNM of a pore network model, wherein the array comprises the following steps:
A1. diameter array h generated in the center of pore by adding hydrateiAdding hydrate to form an array of thickness h on the surface of the porespAddition of hydrate to form thickness array h on throat surfacet；
A2. Setting the pore diameter array in the pore network model as D, the throat length array as L, and the throat diameter array as T, when hydrate grows on the pore surface, the pore diameter array is D-hpWhen the hydrate grows on the surface of the throat, the diameter of the throat is recorded as an array T-ht。
B. In a calculation class in an open source calculation tool-OpenPNM of a pore network model, adding a calculation method for permeability between nodes under the condition that a pore contains a hydrate:
B1. for hydrates generated in the center of the pore, the internodal conduit is divided into 5 sections (see FIG. 2),
as shown in b of FIG. 2, wherein the segments ⑤ - ⑦ are cylindrical, the segments ④, ⑧ are ring-column-shaped, and the permeability (permeability, g) of the ring-column-shaped ring-shaped cross sectionloop) The calculation formula of (2):
wherein D is the outside diameter of the circular section, D is the inside diameter of the circular section, L is the length of the section, and μ is the viscosity of the pore fluid.
Permeability (g) of cylindrical sectionc) The calculation formula is as follows:
the permeability (g) of the entire catheter is expressed as formula 3 according to the Hagen-Poiseuille equation:
B2. for the hydrates generated on the surfaces of the pores and the throat, the calculation formula of the permeability between the conduits replaces the diameters of the pores and the throat with D-hpAnd T-htAnd executing according to the self-contained computing method of the tool OpenPNM. The main calculation method comprises the following steps of;
the fluid flow regime in the internodal conduit is assumed to be laminar. The fluid seepage flow Qij in the throat has the following form:
wherein p is the pressure of the fluid, rij and lij are the radius and length of the throat respectively, mu is the viscosity coefficient of the fluid, and gij is the throat backflow rate.
Assuming the fluid in the formation pore is incompressible, the sum of the flow rates of all throats communicating with it should be zero for a single pore, i.e.:
∑Qijeither case 0 (formula 5)
Therefore, a linear equation system with the pressure as an unknown number can be obtained, and the pressure distribution of the formation pore space can be obtained after the solution. And finally, calculating the inlet flow and the outlet flow of the stratum, and solving the absolute permeability of the stratum by applying Darcy's law. And calculating the permeability of the natural gas hydrate through the steps.
The following example illustrates the implementation of the method of the present invention by generating a pore network of a cubic wire frame with a three-directional number of 100 × 100 × 100 in cartesian coordinates, setting the pore node spacing to 0.63mm, the pore diameter (D) distribution to a normal distribution, a median of the distribution to 0.42mm, and a standard deviation to 0.1mm, generating the diameters of all 1000000 pores by a random algorithm, the throat diameter to be 0.7 times the diameter of the adjacent smaller pores, and the geometric features of the above pore network model to be recorded in the form of a dictionary.
The hydrate diameter at each step was [0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9] × D (pore diameter) assuming hydrate formation at the center of the pore and setting the number of hydrate growth steps to 10.
The permeability of the model of the pore network containing the hydrate was calculated. And C, calling the inter-node catheter permeability calculation model in the step B when the OpenPNM is used for calculating the permeability. Mixing [0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9]]× D (pore diameter) as the hydrate diameter are sequentially substituted into the calculation to obtain the permeability under different hydrate sizes, and then the hydrate saturation S can be obtained by the following formulah：
In the formula, ∑ VhydrateSum of all hydrate volumes, ∑ VporeIs the sum of the volumes of all pore spaces.
It is noted that the disclosed embodiments are intended to aid in further understanding of the invention, but those skilled in the art will appreciate that: various substitutions and modifications are possible without departing from the spirit and scope of the invention and appended claims. Therefore, the invention should not be limited to the embodiments disclosed, but the scope of the invention is defined by the appended claims.
Claims (4)
1. A method for obtaining the permeability of the sediments containing the natural gas hydrate comprises the steps of establishing a pore network model under the condition that the pores of the sediments are filled with the natural gas hydrate, adding a pore filling type hydrate model into the pore throat space of the pore network model, establishing a permeability conduction equation under the condition, and calculating to obtain the permeability of the pore network model containing the natural gas hydrate;
the method comprises the following steps:
A. establishing a pore network model, and adding a pore filling type hydrate model into a pore throat space of the pore network model; the method comprises the following steps:
A1. the nodes in the established pore network model represent the central points of pores in the sediment; the conduit in the pore network model represents a communication part between the central points of two adjacent pores in the sediment; the catheter comprises: a left node half-pore, a throat between two nodes and a right node half-pore in two adjacent pore nodes;
A2. adding a pore filling type hydrate model in a pore throat space of the pore network model, namely adding geometric characteristic characterization array parameters of the pore hydrate into a geometric class of the model;
A21. the geometric characteristic characterization array parameters of the added pore hydrate comprise: diameter array h of hydrate generated in pore centeriThe thickness array h of the hydrate generated on the surface of the porepThickness array h of hydrate generated on surface of throatt；
A22. Setting the pore diameter array in the pore network model as D, the throat length array as L and the throat diameter array as T, and when hydrate grows on the pore surface, setting the pore diameter array as D-hp(ii) a When the hydrate grows on the surface of the throat, the diameter of the throat is recorded as an array T-ht；
B. Establishing a permeability conduction equation, and calculating the permeability of a pore network model containing the natural gas hydrate by adopting a calculation method of the permeability between nodes under the condition that the pores contain the hydrate; the following operations are performed:
B1. for the hydrate generated in the center of the pore, dividing the conduit between the nodes into 5 sections, including a cylindrical shape and an annular column shape;
by permeability gloopThe permeability of a ring-shaped section having a ring-column shape is represented and calculated by formula 1:
in the formula 1, D is the outer diameter of the annular section, namely the diameter of a pore; h isiThe inner diameter of the annular section is the diameter of the hydrate, L is the length of the section, mu is the viscosity of the pore fluid;
by permeability gcRepresents the cylindrical section permeability, calculated by equation 2:
the permeability g of the internodal catheter is expressed as formula 3:
wherein, gzone1～gzone55 segments into which the inter-node conduit is divided;
B2. for the hydrate generated on the surfaces of the pore and the throat, the diameter of the pore and the diameter of the throat are respectively taken as D-hpAnd T-htThen calculating to obtain the permeability between the catheters;
assuming that the fluid flow state in the conduit between the nodes is laminar flow; the fluid seepage flow Qij in the throat is expressed by equation 4:
wherein p is the pressure of the fluid; rij and lij are the radius and length of the throat, respectively; μ is the viscosity coefficient of the fluid; gij is the throat inverse flow rate;
assuming that the fluid in the formation pores is incompressible, for a single pore, the sum of the flow rates of all throats communicating with it should be zero, expressed as equation 5:
∑Qijeither case 0 (formula 5)
Thereby obtaining a linear equation system with the pressure as an unknown number; solving to obtain the pressure distribution of the formation pores;
calculating the inlet flow and the outlet flow of the stratum, and solving the absolute permeability of the stratum by applying Darcy's law;
and obtaining the permeability of the sediments containing the natural gas hydrate through the steps.
2. The device for acquiring the permeability of the sediments containing the natural gas hydrate comprises a memory and a processor;
the memory for storing a computer program;
the processor, when executing the computer program, is configured to implement the method for obtaining permeability of a gas hydrate-containing deposit according to claim 1.
3. A computer-readable storage medium, characterized in that the storage medium has stored thereon a computer program which, when being executed by a processor, carries out the method for obtaining natural gas hydrate-containing deposit permeability according to claim 1.
4. The system is characterized in that the system is developed for the second time on the basis of a pore network model open source computing tool OpenPNM; the system comprises the following modules:
the pore network model establishing module is used for establishing a pore network model of the sediment;
the pore filling type hydrate model adding module is used for adding a pore filling type hydrate model in the pore throat space of the pore network model;
a permeability conduction equation establishing module, wherein the permeability comprises the permeability of an annular section of an annular column shape and the permeability of a cylindrical section;
and the permeability calculation module is used for calculating the permeability of the natural gas hydrate-containing pore network model by using a calculation method for the permeability between nodes under the condition that the pores contain hydrates.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910454469.5A CN110132818B (en) | 2019-05-28 | 2019-05-28 | Method and system for obtaining permeability of sediments containing natural gas hydrate |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910454469.5A CN110132818B (en) | 2019-05-28 | 2019-05-28 | Method and system for obtaining permeability of sediments containing natural gas hydrate |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110132818A CN110132818A (en) | 2019-08-16 |
| CN110132818B true CN110132818B (en) | 2020-07-28 |
Family
ID=67582616
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910454469.5A Expired - Fee Related CN110132818B (en) | 2019-05-28 | 2019-05-28 | Method and system for obtaining permeability of sediments containing natural gas hydrate |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110132818B (en) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110988307A (en) * | 2019-12-26 | 2020-04-10 | 中国科学院广州能源研究所 | Method for evaluating quality of marine natural gas hydrate reservoir |
| CN111862306B (en) * | 2020-07-03 | 2023-09-19 | 大连理工大学 | CT image-based three-dimensional modeling method for pore filling type hydrate deposit |
| CN112014293A (en) * | 2020-09-07 | 2020-12-01 | 中国石油大学(华东) | Method and device for representing hydrate reservoir seepage capability |
| CN112151125B (en) * | 2020-09-24 | 2022-04-19 | 中国石油大学(华东) | Hydrate occurrence form division method based on pore network model |
| CN112304988B (en) * | 2020-10-29 | 2022-02-01 | 中国石油大学(北京) | Method, device and equipment for determining occurrence state of natural gas hydrate |
| CN112816386B (en) * | 2020-12-31 | 2023-08-18 | 中国石油大学(华东) | Method for measuring permeability of hydrate reservoir in hydrate phase change process |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109632579A (en) * | 2018-12-29 | 2019-04-16 | 西南石油大学 | A kind of shale clay mineral is forced from pipette prediction technique |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6289725B1 (en) * | 1998-04-30 | 2001-09-18 | Bechtel Bwxt Idaho, Llc | Field matric potential sensor |
| CN101846605B (en) * | 2010-03-19 | 2012-07-25 | 中国科学院武汉岩土力学研究所 | Device for testing compressibility of soil containing natural gas hydrate and testing method thereof |
| CN105547964A (en) * | 2016-01-15 | 2016-05-04 | 中国石油大学(华东) | One-way flowing determination method for air-water relative permeability under state of different hydrate saturations |
| CN106780744B (en) * | 2016-12-27 | 2020-04-10 | 中国石油天然气集团公司 | Method for constructing multi-scale three-dimensional digital core by utilizing CT images with different resolutions |
| KR101930068B1 (en) * | 2017-04-03 | 2018-12-17 | 한국기계연구원 | Membrane Filter for Detecting Target Particle, Method for Manufacturing the Same, and Method for Detecting Target Particle Using the Same |
| CN107991212B (en) * | 2017-10-26 | 2020-06-16 | 青岛海洋地质研究所 | Fractal dimension measuring and calculating method for effective pores of hydrate-containing sediments |
| CN108729908B (en) * | 2018-05-21 | 2021-11-09 | 中国石油大学(华东) | Dense oil flow simulation and permeability prediction method based on pore network model |
| CN109557253B (en) * | 2018-11-02 | 2019-07-23 | 广州海洋地质调查局 | A kind of comprehensive hydrate simulation system and its experimental method |
-
2019
- 2019-05-28 CN CN201910454469.5A patent/CN110132818B/en not_active Expired - Fee Related
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109632579A (en) * | 2018-12-29 | 2019-04-16 | 西南石油大学 | A kind of shale clay mineral is forced from pipette prediction technique |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110132818A (en) | 2019-08-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110132818B (en) | Method and system for obtaining permeability of sediments containing natural gas hydrate | |
| Cheng et al. | A rigorous compressible streamline formulation for two-and three-phase black-oil simulation | |
| CN112360422B (en) | Shale gas reservoir fractured horizontal well yield prediction method and system | |
| CN108729908A (en) | A kind of oily flow simulating of densification based on pore network model and Permeability Prediction method | |
| CN106599449A (en) | Well test explanation method for cave volume calculation | |
| CN105574320B (en) | The evaluation method of the effective percolation ability of low permeability sandstone reservoirs | |
| CN109684685A (en) | Production capacity and reservoir method for analyzing stability under a kind of multilateral well hydrate decompression mining conditions | |
| CN107622328A (en) | A kind of shale gas reservoir is segmented more cluster pressure break horizontal well production prediction methods | |
| CN110472372B (en) | Dual-medium-based permeability prediction method and system | |
| CN107132170B (en) | A kind of determination method and apparatus of reservoir stress sensitivity | |
| CN105547964A (en) | One-way flowing determination method for air-water relative permeability under state of different hydrate saturations | |
| CA3108482C (en) | Method, medium, terminal and device for evaluating layered water injection efficiency of oil reservoir | |
| CN113484216B (en) | Method for evaluating water phase flowback rate and reasonable flowback pressure difference of tight sandstone gas reservoir | |
| CN111353205A (en) | Method for calculating stratum pressure and dynamic capacity of water-producing gas well of tight gas reservoir | |
| CN108801888B (en) | Indoor degradation simulation device for organic soil | |
| CN105863625A (en) | Quinary-coefficient classification evaluating method for tight oil reservoir | |
| CN112943176A (en) | Method for measuring gas-oil relative permeability curve for simulating injection and production of oil reservoir type gas storage | |
| CN109357986A (en) | The stifled measuring method for adjusting overall process phase percolation curve of High water cut reservoir long cores | |
| CN110702722B (en) | Based on nuclear magnetic resonance T2Method for calculating fractal dimension of porous medium containing hydrate by spectrum | |
| CN108729912A (en) | Production split method suitable for reservoir numerical simulation | |
| CN115235974A (en) | Hydrate argillaceous silt methane adsorption model construction method and experimental device | |
| CN112151125B (en) | Hydrate occurrence form division method based on pore network model | |
| CN105675441A (en) | Gas-water relative permeability measurement method under radial flow condition at different hydrate saturations | |
| US20230251178A1 (en) | Random particle generation method based on particle size distribution | |
| US20220178252A1 (en) | Method and system for predicting production of fractured horizontal well in shale gas reservoir |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200728 |