CN116362031A - Method and device for analyzing permeability evolution law of hydrate reservoir - Google Patents
Method and device for analyzing permeability evolution law of hydrate reservoir Download PDFInfo
- Publication number
- CN116362031A CN116362031A CN202310271933.3A CN202310271933A CN116362031A CN 116362031 A CN116362031 A CN 116362031A CN 202310271933 A CN202310271933 A CN 202310271933A CN 116362031 A CN116362031 A CN 116362031A
- Authority
- CN
- China
- Prior art keywords
- permeability
- parameters
- hydrate
- dimensional geometric
- geometric models
- 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.)
- Pending
Links
- 230000035699 permeability Effects 0.000 title claims abstract description 109
- 238000000034 method Methods 0.000 title claims abstract description 52
- 238000004364 calculation method Methods 0.000 claims abstract description 53
- 239000011159 matrix material Substances 0.000 claims abstract description 53
- 239000011148 porous material Substances 0.000 claims abstract description 53
- 239000002245 particle Substances 0.000 claims abstract description 51
- 238000004088 simulation Methods 0.000 claims abstract description 33
- 150000004677 hydrates Chemical class 0.000 claims abstract description 29
- QNRATNLHPGXHMA-XZHTYLCXSA-N (r)-(6-ethoxyquinolin-4-yl)-[(2s,4s,5r)-5-ethyl-1-azabicyclo[2.2.2]octan-2-yl]methanol;hydrochloride Chemical compound Cl.C([C@H]([C@H](C1)CC)C2)CN1[C@@H]2[C@H](O)C1=CC=NC2=CC=C(OCC)C=C21 QNRATNLHPGXHMA-XZHTYLCXSA-N 0.000 claims description 13
- 238000003860 storage Methods 0.000 claims description 12
- 230000002572 peristaltic effect Effects 0.000 claims description 9
- -1 parallelepiped Chemical compound 0.000 claims description 6
- 238000011160 research Methods 0.000 abstract description 12
- 239000004576 sand Substances 0.000 description 25
- 239000012530 fluid Substances 0.000 description 19
- 238000010586 diagram Methods 0.000 description 11
- 238000004590 computer program Methods 0.000 description 7
- 230000008569 process Effects 0.000 description 7
- 238000009826 distribution Methods 0.000 description 6
- 230000008859 change Effects 0.000 description 5
- 238000013508 migration Methods 0.000 description 5
- 230000005012 migration Effects 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 238000000354 decomposition reaction Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 3
- 238000005192 partition Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 239000004927 clay Substances 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 210000004027 cell Anatomy 0.000 description 1
- 210000005056 cell body Anatomy 0.000 description 1
- 230000008094 contradictory effect Effects 0.000 description 1
- 229910003460 diamond Inorganic materials 0.000 description 1
- 239000010432 diamond Substances 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 238000011065 in-situ storage Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Natural products C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 1
- 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 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 239000013049 sediment Substances 0.000 description 1
- 239000012798 spherical particle Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Abstract
The application discloses a method and a device for analyzing a permeability evolution rule of a hydrate reservoir. The method comprises the following steps: establishing a plurality of three-dimensional geometric models of pore-filling type hydrates with different matrix sizes and different spatial arrangement structures; acquiring parameters of a plurality of three-dimensional geometric models; based on the parameters, respectively carrying out numerical simulation calculation on the plurality of three-dimensional geometric models to obtain a global calculation result; extracting target data from the global calculation result; and determining the evolution rule of the permeability of the hydrate reservoir according to the target data. The method can accurately simulate the influence of the generation of pore filling type hydrates under different particle sizes and matrix space arrangement structures in an actual stratum environment on the permeability. The error caused by the theoretical calculation simplified model is greatly reduced, so that the research result is closer to the actual stratum environment, and the method has important guiding significance on the research of permeability evolution rules for the generation of pore filling type hydrates.
Description
Technical Field
The application relates to the technical field of natural gas hydrates in sea areas, in particular to a method and a device for analyzing a permeability evolution rule of a hydrate reservoir.
Background
The natural gas hydrate has large reserves, belongs to low-carbon clean energy, but can not meet the requirements of commercial exploitation due to the storage property and the gas production property. In situ decomposition of hydrates within a reservoir is accompanied by substantial phase change mass transfer, where permeability is a critical parameter affecting fluid flow, determining hydrate decomposition, migration and production efficiency.
In the prior study on the permeability of a hydrate reservoir, the complexity of a theoretical model excessively depends on a large amount of stratum fine data, so that the difficulty of theoretical calculation is increased. In addition, the existing theoretical model is mostly convenient for research, corresponding simplification treatment is carried out, and dynamic occurrence and evolution migration mechanisms of the hydrate are ignored. Therefore, how to truly reflect the decomposition of hydrate, sand migration and evolution phenomenon of pore structure in the actual stratum environment, and explore the influence of the generation of pore filling type hydrate on the permeability evolution rule, becomes a technical problem to be solved urgently.
Disclosure of Invention
The embodiment of the application aims to provide a method and a device for analyzing the permeability evolution law of a hydrate reservoir, which are used for solving the problem that a theoretical model in the prior art has larger calculation error on the permeability of the hydrate reservoir.
To achieve the above object, a first aspect of the present application provides a method for analyzing a permeability evolution law of a hydrate reservoir, including:
establishing a plurality of three-dimensional geometric models of pore-filling type hydrates with different matrix sizes and different spatial arrangement structures;
acquiring parameters of a plurality of three-dimensional geometric models;
based on the parameters, respectively carrying out numerical simulation calculation on the plurality of three-dimensional geometric models to obtain a global calculation result;
extracting target data from the global calculation result;
and determining the evolution rule of the permeability of the hydrate reservoir according to the target data.
In an embodiment of the present application, creating a plurality of three-dimensional geometric models of pore-containing filled hydrates of different matrix sizes and different spatial arrangements includes:
establishing a plurality of initial three-dimensional geometric models;
and respectively carrying out space structure arrangement on the plurality of initial three-dimensional geometric models according to the particle sizes and the space arrangement structures of the matrix particles of the plurality of initial three-dimensional geometric models so as to obtain a plurality of three-dimensional geometric models.
In an embodiment of the present application, the spatial arrangement structure includes:
cube, rhombic, parallelepiped, and rhombic.
In an embodiment of the present application, the parameters of the plurality of three-dimensional geometric models include:
model parameters, peristaltic flow physical field parameters, and hydrate simulation parameters.
In this embodiment of the present application, performing numerical simulation calculation on the plurality of three-dimensional geometric models to obtain a global calculation result includes:
inputting a plurality of three-dimensional geometric models into finite element software;
acquiring preset corresponding simulation boundary condition parameters;
dividing grids by a plurality of three-dimensional geometric models according to the simulated boundary condition parameters;
and carrying out numerical simulation calculation on the multiple three-dimensional models after grid division to obtain a global calculation result.
In an embodiment of the present application, the method further includes:
and obtaining preset output parameters in finite element software, wherein the output parameters are porosity parameters, hydrate saturation parameters and permeability parameters.
In an embodiment of the present application, the target data includes porosity data, and hydrate saturation data and permeability data, and determining an evolution rule of the permeability of the hydrate reservoir according to the target data includes:
and outputting a relation curve of the saturation and the matrix size of the pore filling type hydrate and the permeability according to the model parameters, the porosity data, the saturation data and the permeability data of the pore filling type hydrate so as to determine the evolution rule of the permeability of the hydrate reservoir.
In embodiments of the present application, the hydrate reservoir permeability satisfies equation (1):
wherein k is permeability, Q is flow, μ is viscosity, L is seepage path length, A is flow channel cross-sectional area, P in For inlet pressure, P out Is the outlet pressure.
A second aspect of the present application provides an apparatus for analyzing the permeability evolution law of a hydrate reservoir, comprising:
a memory configured to store instructions; and
a processor configured to invoke the instructions from the memory and when executing the instructions is capable of implementing a method for analyzing a hydrate reservoir permeability evolution law according to the above.
A third aspect of the present application provides a machine-readable storage medium having stored thereon instructions for causing a machine to perform a method for analysing a hydrate reservoir permeability evolution law according to the above.
Through the technical scheme, a plurality of three-dimensional geometric models of pore-filling type hydrates with different matrix sizes and different space arrangement structures are established; acquiring parameters of a plurality of three-dimensional geometric models; based on the parameters, respectively carrying out numerical simulation calculation on the plurality of three-dimensional geometric models to obtain a global calculation result; extracting target data from the global calculation result; and determining the evolution rule of the permeability of the hydrate reservoir according to the target data. The method can accurately simulate the influence of the generation of pore-containing filling type hydrates under different particle and matrix space arrangement structures in the actual stratum environment on the permeability. The error caused by the theoretical calculation simplified model is greatly reduced, so that the research result is closer to the actual stratum environment, and the method has important guiding significance on the research of permeability evolution rules for the generation of pore filling type hydrates.
Additional features and advantages of embodiments of the present application will be set forth in the detailed description that follows.
Drawings
The accompanying drawings are included to provide a further understanding of embodiments of the present application and are incorporated in and constitute a part of this specification, illustrate embodiments of the present application and together with the description serve to explain, without limitation, the embodiments of the present application. In the drawings:
FIG. 1 schematically illustrates a flow chart of a method for analyzing a hydrate reservoir permeability evolution law according to an embodiment of the present application;
FIG. 2 schematically illustrates a graph of a random distribution model of cubic pore-filled hydrate particles according to an embodiment of the present application;
FIG. 3 schematically illustrates a grid partition of a random distribution of cubic pore-filled hydrate particles according to an embodiment of the present application;
FIG. 4 schematically illustrates an arrangement of porous media fully homogeneous isodiametric sand spheres according to an embodiment of the present application;
FIG. 5 schematically illustrates a pore structure of a porous media fully homogeneous isodiametric sand sphere according to an embodiment of the present application;
FIG. 6 schematically illustrates a grid partition of a porous media fully homogeneous isodiametric sand sphere according to an embodiment of the present application;
FIG. 7 schematically illustrates a graph of output results according to a specific embodiment of the present application;
FIG. 8 schematically illustrates a graph of pore-filled hydrate saturation and matrix size versus permeability according to an embodiment of the present application;
fig. 9 schematically shows a block diagram of an apparatus for analyzing the permeability evolution law of a hydrate reservoir according to an embodiment of the present application.
Detailed Description
For the purposes of making the objects, technical solutions and advantages of the embodiments of the present application more clear, the technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application, and it should be understood that the specific implementations described herein are only for illustrating and explaining the embodiments of the present application, and are not intended to limit the embodiments of the present application. All other embodiments, which can be made by one of ordinary skill in the art without undue burden from the present disclosure, are within the scope of the present application based on the embodiments herein.
It should be noted that, in the embodiment of the present application, directional indications (such as up, down, left, right, front, and rear … …) are referred to, and the directional indications are merely used to explain the relative positional relationship, movement conditions, and the like between the components in a specific posture (as shown in the drawings), and if the specific posture is changed, the directional indications are correspondingly changed.
In addition, if there is a description of "first", "second", etc. in the embodiments of the present application, the description of "first", "second", etc. is for descriptive purposes only and is not to be construed as indicating or implying a relative importance or implicitly indicating the number of technical features indicated. Thus, a feature defining "a first" or "a second" may explicitly or implicitly include at least one such feature. In addition, the technical solutions of the embodiments may be combined with each other, but it is necessary to base that the technical solutions can be realized by those skilled in the art, and when the technical solutions are contradictory or cannot be realized, the combination of the technical solutions should be regarded as not exist and not within the protection scope of the present application.
Fig. 1 schematically shows a flow chart of a method for analyzing the evolution law of permeability of a hydrate reservoir according to an embodiment of the present application. As shown in fig. 1, an embodiment of the present application provides a method for analyzing a hydrate reservoir permeability evolution law, which may include the following steps.
102, acquiring parameters of a plurality of three-dimensional geometric models;
104, extracting target data from the global calculation result;
and 105, determining the evolution rule of the permeability of the hydrate reservoir according to the target data.
In the embodiment of the application, in order to truly reflect the decomposition of the hydrate, the sand migration and the pore structure evolution law in the actual stratum environment, so as to analyze the influence of the generation of the pore filling type hydrate on the permeability, the application provides a method for analyzing the permeability evolution law of the hydrate reservoir. Firstly, a plurality of three-dimensional geometric models of pore-filling type hydrates with different matrix sizes and different spatial arrangement structures can be established through numerical simulation software drawing, the three-dimensional models are composed of sand spheres with completely homogeneous porous media and equal diameters, the matrix sizes are the particle sizes of the sand spheres, and the spatial arrangement structures refer to the spatial arrangement structures of the sand spheres. In one aspect, the spatial arrangement of the hydrates is related to the porosity, which affects the saturation of the hydrates and thus the permeability. On the other hand, the particle size of the sand constituting the hydrate may also affect the permeability. Thus, there is a need to build multiple three-dimensional geometric models of pore-filled hydrates of different matrix sizes and different spatial arrangements.
After the three-dimensional geometric model of the pore filling type hydrate is established, parameter setting can be carried out on the model through simulation software. In embodiments of the present application, the three-dimensional geometric model parameters of the pore-filling type hydrate may include model parameters, peristaltic flow physical field parameters, and hydrate simulation parameters. The model parameters are used for defining a homogeneous spherical sand grain arrangement structure, and mainly comprise the number of particle contact points, the unit volume and the porosity; the laminar flow physical parameters are provided with an inlet speed, an inlet pressure, an outlet speed, an outlet pressure and a non-slip boundary; the hydrate simulation parameters are hydrate matrix radius, density, viscosity, matrix unit side length, hydrate saturation and hydrate particle radius. The grid controlled by the physical field is provided with a boundary layer grid near the wall, and based on the setting of parameters, numerical simulation calculation can be performed on the three-dimensional geometric model to obtain a global calculation result. The data of the evolution of the pore filling type hydrate generated on the permeability, namely the target data, of the internal matrix particles containing the hydrate sediment under different stacking forms can be extracted from the calculation result. According to the extracted target data, a relation graph of the saturation and the permeability of the pore filling type hydrate can be generated, the change trend of the relative and absolute permeability of the pore filling type hydrate can be found according to the graph along with the difference of particle sizes of the particles, and further the evolution condition of the permeability under the conditions of different matrixes and particle sizes and different saturation of the hydrate is analyzed and discussed, so that the evolution rule of the permeability of the hydrate reservoir is found.
Fig. 2 schematically shows a model diagram of a random distribution of cubic pore-filling type hydrate particles according to an embodiment of the present application, and fig. 3 schematically shows a grid-dividing diagram of a random distribution of cubic pore-filling type hydrate particles according to an embodiment of the present application. In practical situations, the occurrence of pore-filling hydrates in a real stratum is randomly distributed, so that a cube model with random distribution as shown in fig. 2 needs to be established to characterize the distribution, and then grid division as shown in fig. 3 is performed for calculation. In the examples of the present application, the ideal state, that is, the fixed position of the center of sphere of the pore-filling type hydrate formation will be described as an example.
Through the technical scheme, a plurality of three-dimensional geometric models of pore-filling type hydrates with different matrix sizes and different space arrangement structures are established; acquiring parameters of a plurality of three-dimensional geometric models; based on the parameters, respectively carrying out numerical simulation calculation on the plurality of three-dimensional geometric models to obtain a global calculation result; extracting target data from the global calculation result; and determining the evolution rule of the permeability of the hydrate reservoir according to the target data. The method can accurately simulate the influence of the generation of pore-containing filling type hydrates under different particle and matrix space arrangement structures in the actual stratum environment on the permeability. The error caused by the theoretical calculation simplified model is greatly reduced, so that the research result is closer to the actual stratum environment, and the method has important guiding significance on the research of permeability evolution rules for the generation of pore filling type hydrates. Fig. 4 schematically shows an arrangement structure of sand spheres of a completely homogeneous isodiametric porous medium according to an embodiment of the present application. As shown in fig. 4, in an embodiment of the present application, creating multiple three-dimensional geometric models of pore-containing filled hydrates of different matrix sizes and different spatial arrangements may include:
establishing a plurality of initial three-dimensional geometric models;
and respectively carrying out space structure arrangement on the plurality of initial three-dimensional geometric models according to the particle sizes and the space arrangement structures of the matrix particles of the plurality of initial three-dimensional geometric models so as to obtain a plurality of three-dimensional geometric models.
In an embodiment of the present application, the spatial arrangement structure may include:
cube, rhombic, parallelepiped, and rhombic.
In the embodiment of the application, since the permeability of the hydrate reservoir is related to the spatial arrangement structure and the matrix size of the pore-containing filled hydrate, the evolution rule of the permeability of the hydrate reservoir can be analyzed by establishing three-dimensional geometric models of different spatial arrangement structures and matrix sizes. Firstly, a plurality of initial three-dimensional geometric models can be established, and then the plurality of initial three-dimensional geometric models are spatially arranged according to the particle sizes and the spatially arranged structures of matrix particles of the plurality of initial three-dimensional geometric models, so that a plurality of three-dimensional geometric models are obtained. In one example, the matrix particles can be classified into four types of medium sand, fine sand, silt and clay according to the size of the matrix particles. Wherein, the grain diameter D of the medium sand is 250.0 mu m less than D and less than or equal to 500.0 mu m; the grain diameter D of the fine sand is 62.5 mu m less than D and less than or equal to 250.0 mu m; the grain diameter D of the silt is more than 2.0 mu m and less than or equal to 62.5 mu m; the particle size D of clay is less than or equal to 2.0 mu m. It should be noted that, the three-dimensional geometric models obtained after the spatial structure arrangement can be divided into four types according to the spatial arrangement structure, and as shown in fig. 4, the three-dimensional geometric models are cubic, oblique cubic, parallelepiped and rhombic in order from left to right. Each spatial arrangement type in turn comprises a plurality of three-dimensional geometric models of different particle sizes.
Fig. 5 schematically illustrates a pore structure of a porous media fully homogeneous isodiametric sand sphere according to an embodiment of the present application. As shown in fig. 5, pore structures of the cube type, the rhombic cube type, the parallelepiped type and the rhombic type are respectively arranged from left to right, and the centers of the hydrate spherical particles are respectively positioned at the centers of the largest pores of each spatial arrangement structure. In fig. 5, the circular hatched portion is a pore-filled hydrate particle, the transparent portion is a pore not occupied by hydrate, and the relatively narrow space in which the sand particles are mutually communicated is pore throat. Wherein, because of the limitation of the geometric relation of the sand particles and the hydrate of the model, when the diameter of the spherical hydrate is equal to the distance between the sand particles, the saturation of the hydrate reaches the maximum. The three-dimensional geometric models of different spatial arrangements have different porosities, and the larger the porosity is, the higher the saturation of the hydrate is. In one example, the maximum saturation of hydrates within the spatial arrangement of the three-dimensional geometric model is 42.73%, 44.33%, 2.52%, 19.76%, respectively. Among the four space arrangement structure types, the connection line of the sphere center of the cube is cube, the included angle of each intersecting straight line is 90 degrees, and the arrangement structure has the largest porosity, so that the arrangement structure is the most ideal space arrangement structure in the actual stratum; the connecting line of the sphere centers of the rhombus is in the shape of a rhombus quadrangular prism, the horizontal planes are arranged in a square shape, and the vertical planes are in the shape of a rhombus with included angles of 60 degrees and 120 degrees; the shape of the connecting line of the sphere centers of the parallelepipeds is an oblique quadrangular prism, the horizontal plane is arranged into a rhombus with an included angle of 60 degrees and 120 degrees, and the vertical direction is arranged into a rhombus with an included angle of 60 degrees and 120 degrees; the shape of the connecting line of the sphere centers of the diamond shapes is a regular rectangular pyramid with equal sides, and the arrangement structure has minimum porosity. In the three-dimensional geometric model, each sand has a fixed and equal number of adjacent contact sand particles, and four types of spatial arrangements are 6, 8, 10 and 12, respectively. According to the configuration of the unit body, the calculation mode of the volume, the porosity and the saturation of the matrix particles satisfies the formula (2):
wherein: r is the radius of matrix particles in the unit body;
initial porosity for the unit cell; v (V) p Is the volume of matrix particles in the cell body; v (V) U Is the volume of the unit body; v (V) H Is the volume of hydrate in the unit body.
In an embodiment of the present application, the parameters of the plurality of three-dimensional geometric models may include:
model parameters, peristaltic flow physical field parameters, and hydrate simulation parameters.
In the embodiment of the application, after the three-dimensional geometric model of the pore filling type hydrate is established, parameter setting can be performed on the model through simulation software. In embodiments of the present application, the three-dimensional geometric model parameters of the pore-filling type hydrate may include model parameters, peristaltic flow physical field parameters, and hydrate simulation parameters. The model parameters are used for defining a homogeneous spherical sand grain arrangement structure, and mainly comprise the number of particle contact points, the unit volume and the porosity; the laminar flow physical parameters are provided with an inlet speed, an inlet pressure, an outlet speed, an outlet pressure and a non-slip boundary; the hydrate simulation parameters are the particle size, density, viscosity, matrix unit side length and the particle size of the hydrate particles.
After the three-dimensional geometric model is established and basic model parameters of the three-dimensional geometric model are set, the physical field parameters of the three-dimensional geometric model can be set continuously. Wherein the physical field parameters may include: inlet velocity, inlet pressure, outlet velocity, and outlet pressure. In practical research, a technician can select proper physical field parameters according to practical requirements. In one example, the physical field parameters are set according to actual engineering study, such as inlet pressure is set to 1MPa, fluid outflow is free, and thus both outlet velocity and outlet pressure are set to free state here. The absolute permeability of matrix particle arrangement is calculated by adopting a liquid permeability measuring method, liquid water is used as the permeability to calculate fluid, and dynamic viscosity is 0.01 Pa.s. In another example, where a physical field is selected, fluid flow should be selected, and the fluid flow regime is peristaltic flow due to a Reynolds number Re well below 1. In the case of selecting a study tree, steady state in the general study should be selected.
Further, the plurality of three-dimensional geometric model parameters further includes hydrate simulation parameters, namely hydrate matrix radius, fluid density, fluid viscosity, matrix unit side length, hydrate particle radius, and hydrate saturation. The parameters can be set according to practical conditions, for example, the radius of the substrate can be set to be 1 multiplied by 10 -6 m, the fluid density can be set to 1000kg/m3, the fluid viscosity can be set to 0.01 Pa.s, and the side length of the matrix unit can be set to 2×10 -6 m, the radius of the hydrate particles can be 7.3205 multiplied by 10 -13 m。
Fig. 6 schematically illustrates a grid partition of a porous medium fully homogeneous isodiametric sand sphere according to an embodiment of the present application. As shown in fig. 6, in the embodiment of the present application, performing numerical simulation calculation on the plurality of three-dimensional geometric models to obtain a global calculation result may include:
inputting a plurality of three-dimensional geometric models into finite element software;
acquiring preset corresponding simulation boundary condition parameters;
dividing grids by a plurality of three-dimensional geometric models according to the simulated boundary condition parameters;
and carrying out numerical simulation calculation on the multiple three-dimensional models after grid division to obtain a global calculation result.
In the embodiment of the application, after the establishment of the plurality of three-dimensional geometric models is completed, numerical simulation calculation can be performed. Firstly, inputting a plurality of three-dimensional geometric models into finite element software, and then obtaining preset corresponding simulation boundary condition parameters to grid the three-dimensional geometric models so as to perform global calculation. The preset corresponding simulated boundary condition parameters comprise a unidirectional peristaltic flow set in the physical field tree, a fluid inlet boundary at the unidirectional peristaltic flow, an inlet pressure and an inlet speed at the inlet boundary, and an outlet boundary set at the physical field window, wherein the inlet boundary and the outlet boundary show symmetry. The boundary is set to limit the fluid with a certain pressure and speed from flowing in from the inlet boundary port, flowing in the pores and pore throats of the matrix, and flowing out from a certain outlet boundary, wherein the freely flowing fluid with a certain speed and pressure flows in the pore throats, and the process of flowing in the pore throats can be regarded as a seepage phenomenon. In practical research, the boundary between fluid inflow and fluid outflow is not limited, and the boundary can be set under the condition of meeting symmetry. Wherein, no slip boundary is set up at the rock matrix grain surface, hydrate grain surface and the upper and lower boundary of the model, and the displacement size.
Further, the simulated boundary condition parameters may also include a grid seed point parameter that is used to indicate the granularity at which the geometric model is grid partitioned. The smaller the preset seed point parameter is, the smaller the granularity of the finite element software for carrying out grid division on the geometric model is, and the higher the precision of finite element calculation on the geometric model is. The finite element software can automatically divide grids of the geometric model according to the set seed point parameters, and the grid division adopts a refined grid controlled by a physical field so as to ensure that a large speed gradient can be analyzed. In addition, the physical field controlled grid also provides boundary layer grids near the walls. Preferably, the simulated boundary condition parameters may also include fluid volume variables, which may be set at the mass properties, which are calculated porosity and permeability parameters.
In an embodiment of the present application, the method may further include:
and obtaining preset output parameters in finite element software, wherein the output parameters are porosity parameters, hydrate saturation parameters and permeability parameters.
In the embodiment of the application, the output parameters can be set as the porosity parameter, the permeability parameter and the hydrate saturation parameter in the input control corresponding to the output manager of the finite element software, wherein the parameters output by the output manager are evolution phenomenon data of the permeability and the different hydrate particle sizes under the simulation of different matrix particle sizes by the finite element software, and the process and result data of the simulated whole permeability phenomenon can be output. For example, the data of absolute permeability and saturation of the hydrate can be output under different sand particle sizes, and the evolution rule of the saturation of the pore-filling type hydrate to the permeability can be obtained through analysis of a data source.
Fig. 7 schematically illustrates a graph of output results according to a specific embodiment of the present application. As shown in fig. 7, if the inlet pressure of the finite element software simulation permeability phenomenon is 1MPa, fig. 7 is the migration phenomenon of the fluid in different matrix particles under the pressure and speed conditions. The direction of the arrow streamline in the figure is the direction of fluid flow, the color of the arrow streamline represents the pressure of the fluid in the pore throat of the matrix, and the color of the section of the model represents the flow velocity of the fluid in the model.
As shown in fig. 7, in the basic unit flow field diagram, the color of the central part of the cubic flow field is lighter, the speed is more obvious, the flow field is wider, and the permeation phenomenon is more obvious. In practical application, corresponding curves can be output according to requirements, so that the evolution condition of the saturation of the pore filling type hydrate on the permeability in the whole simulation process can be analyzed according to the field curves and the permeability data extracted by post-treatment. And extracting the relation between the permeability and the hydrate saturation under different arrangement structures through post-treatment according to the output final result, and analyzing and processing the relation.
Fig. 8 schematically illustrates a graph of pore-filled hydrate saturation and matrix size versus permeability according to an embodiment of the present application. As shown in fig. 8, in an embodiment of the present application, the target data includes porosity data, hydrate saturation data, and permeability data, and determining the evolution law of the permeability of the hydrate reservoir according to the target data may include:
and outputting a relation curve of the pore filling type hydrate saturation and the matrix size and the permeability according to the model parameters, the porosity data, the saturation data and the permeability data of the pore filling type hydrate so as to determine the evolution rule of the permeability of the hydrate reservoir.
In embodiments of the present application, the target data includes porosity data, hydrate saturation data, and permeability data. Since the spatial arrangement of the hydrates is related to the porosity, which affects the saturation of the hydrates and thus the permeability, the saturation data of the pore-filled hydrates can be determined from the porosity data. At the same time, the particle size of the sand particles, i.e. the matrix size, that make up the hydrate can also affect the permeability. In embodiments of the present application, the hydrate reservoir permeability satisfies equation (1):
wherein k is permeability, Q is flow, μ is viscosity, L is seepage path length, A is flow channel cross-sectional area, P in For inlet pressure, P out Is the outlet pressure.
After the permeability data is calculated according to the formula (1), the permeability change under different matrix sizes and different hydrate saturation can be extracted from the permeability data, and when the matrix sizes are different, the permeability of the hydrate can also change correspondingly. For example, in four basic pore structures, namely, cubic, rhombic, parallelepiped, and rhombic, the unit volumes thereof are different depending on the geometrical arrangement. According to the calculation results, the saturation and permeability data of the pore filling type hydrate under the four arrangement structures can be extracted and output. And further quantitatively analyzing the evolution law of the saturation of the pore filling type hydrate on the permeability, namely, combining the saturation data of the pore filling type hydrate, outputting the relation curve of the saturation of the pore filling type hydrate, the size of a matrix and the permeability, so as to determine the evolution law of the permeability of the hydrate reservoir, and further provide a theoretical basis for researching the permeability of the hydrate in the actual stratum environment.
Compared with the method for simulating the permeability evolution phenomenon by using the theoretical model and the test mode, the method can simulate the permeability evolution phenomenon by using the finite element software, can avoid the situation that the theoretical model is too simplified and is difficult to consider the complex stratum environment, can also avoid the situation that the actual stratum environment is difficult to reproduce under the laboratory condition, and causes a great deal of physical and financial resources to be consumed, can simulate the actual stratum condition to a greater extent by researching the algorithm, ensures that the research is more accurate and has guidance, and simultaneously effectively controls the error and the danger coefficient.
Fig. 9 schematically shows a block diagram of an apparatus for analyzing the permeability evolution law of a hydrate reservoir according to an embodiment of the present application. As shown in fig. 9, an embodiment of the present application provides a controller, which may include:
a memory 910 configured to store instructions; and
the processor 920 is configured to call instructions from the memory 910 and when executing the instructions, to implement the method for controlling the boom described above.
Specifically, in embodiments of the present application, the processor 920 may be configured to:
establishing a plurality of three-dimensional geometric models of pore-filling type hydrates with different matrix sizes and different spatial arrangement structures;
acquiring parameters of a plurality of three-dimensional geometric models;
based on the parameters, respectively carrying out numerical simulation calculation on the plurality of three-dimensional geometric models to obtain a global calculation result;
extracting target data from the global calculation result;
and determining the evolution rule of the permeability of the hydrate reservoir according to the target data.
Further, the processor 920 may be further configured to:
establishing a plurality of initial three-dimensional geometric models;
and respectively carrying out space structure arrangement on the plurality of initial three-dimensional geometric models according to the particle sizes and the space arrangement structures of the matrix particles of the plurality of initial three-dimensional geometric models so as to obtain a plurality of three-dimensional geometric models.
In an embodiment of the present application, the spatial arrangement structure includes:
cube, rhombic, parallelepiped, and rhombic.
In an embodiment of the present application, the parameters of the plurality of three-dimensional geometric models include:
model parameters, peristaltic flow physical field parameters, and hydrate simulation parameters.
Further, the processor 920 may be further configured to:
inputting a plurality of three-dimensional geometric models into finite element software;
acquiring preset corresponding simulation boundary condition parameters;
dividing grids by a plurality of three-dimensional geometric models according to the simulated boundary condition parameters;
and carrying out numerical simulation calculation on the multiple three-dimensional models after grid division to obtain a global calculation result.
Further, the processor 920 may be further configured to:
and obtaining preset output parameters in finite element software, wherein the output parameters are porosity parameters, hydrate saturation parameters and permeability parameters.
Further, the processor 920 may be further configured to:
and outputting a relation curve of the pore filling type hydrate saturation and the matrix size and the permeability according to the model parameters, the porosity data, the saturation data and the permeability data of the pore filling type hydrate so as to determine the evolution rule of the permeability of the hydrate reservoir.
In embodiments of the present application, the hydrate reservoir permeability satisfies equation (1):
wherein k is permeability, Q is flow, μ is viscosity, L is seepage path length, A is flow channel cross-sectional area, P in For inlet pressure, P out Is the outlet pressure.
Through the technical scheme, a plurality of three-dimensional geometric models of pore-filling type hydrates with different matrix sizes and different space arrangement structures are established; acquiring parameters of a plurality of three-dimensional geometric models; based on the parameters, respectively carrying out numerical simulation calculation on the plurality of three-dimensional geometric models to obtain a global calculation result; extracting target data from the global calculation result; and determining the evolution rule of the permeability of the hydrate reservoir according to the target data. The method can accurately simulate the influence of the generation of pore-containing filled hydrates under different particle and matrix space arrangement structures in the actual stratum environment on the permeability. The error caused by the theoretical calculation simplified model is greatly reduced, so that the research result is closer to the actual stratum environment, and the method has important guiding significance on the research of permeability evolution rules for the generation of pore filling type hydrates.
Embodiments of the present application also provide a machine-readable storage medium having instructions stored thereon for causing a machine to perform the above-described method for analyzing a hydrate reservoir permeability evolution law.
It will be appreciated by those skilled in the art that embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
In one typical configuration, a computing device includes one or more processors (CPUs), input/output interfaces, network interfaces, and memory.
The memory may include volatile memory in a computer-readable medium, random Access Memory (RAM) and/or nonvolatile memory, etc., such as Read Only Memory (ROM) or flash RAM. Memory is an example of a computer-readable medium.
Computer readable media, including both non-transitory and non-transitory, removable and non-removable media, may implement information storage by any method or technology. The information may be computer readable instructions, data structures, modules of a program, or other data. Examples of storage media for a computer include, but are not limited to, phase change memory (PRAM), static Random Access Memory (SRAM), dynamic Random Access Memory (DRAM), other types of Random Access Memory (RAM), read Only Memory (ROM), electrically Erasable Programmable Read Only Memory (EEPROM), flash memory or other memory technology, compact disc read only memory (CD-ROM), digital Versatile Discs (DVD) or other optical storage, magnetic cassettes, magnetic tape disk storage or other magnetic storage devices, or any other non-transmission medium, which can be used to store information that can be accessed by a computing device. Computer-readable media, as defined herein, does not include transitory computer-readable media (transmission media), such as modulated data signals and carrier waves.
It should also be noted that the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising one … …" does not exclude the presence of other like elements in a process, method, article or apparatus that comprises an element.
The foregoing is merely exemplary of the present application and is not intended to limit the present application. Various modifications and changes may be made to the present application by those skilled in the art. Any modifications, equivalent substitutions, improvements, etc. which are within the spirit and principles of the present application are intended to be included within the scope of the claims of the present application.
Claims (10)
1. A method for analyzing the permeability evolution law of a hydrate reservoir, comprising:
establishing a plurality of three-dimensional geometric models of pore-filling type hydrates with different matrix sizes and different spatial arrangement structures;
acquiring parameters of the plurality of three-dimensional geometric models;
based on the parameters, respectively carrying out numerical simulation calculation on the plurality of three-dimensional geometric models to obtain a global calculation result;
extracting target data from the global calculation result;
and determining the evolution rule of the permeability of the hydrate reservoir according to the target data.
2. The method of claim 1, wherein the creating a plurality of three-dimensional geometric models of pore-containing filled hydrates of different matrix sizes and different spatial arrangements comprises:
establishing a plurality of initial three-dimensional geometric models;
and respectively carrying out space structure arrangement on the plurality of initial three-dimensional geometric models according to the particle sizes and the space arrangement structures of the matrix particles of the plurality of initial three-dimensional geometric models so as to obtain the plurality of three-dimensional geometric models.
3. The method of claim 1, wherein the spatial arrangement comprises:
cube, rhombic, parallelepiped, and rhombic.
4. The method of claim 1, wherein the parameters of the plurality of three-dimensional geometric models comprise:
model parameters, peristaltic flow physical field parameters, and hydrate simulation parameters.
5. The method of claim 1, wherein performing numerical simulation calculations on the plurality of three-dimensional geometric models to obtain global calculation results comprises:
inputting the plurality of three-dimensional geometric models into finite element software;
acquiring preset corresponding simulation boundary condition parameters;
respectively meshing the plurality of three-dimensional geometric models according to the simulated boundary condition parameters;
and carrying out numerical simulation calculation on the multiple three-dimensional models after grid division to obtain a global calculation result.
6. The method according to claim 1, wherein the method further comprises:
and obtaining preset output parameters in finite element software, wherein the output parameters are porosity parameters, hydrate saturation parameters and permeability parameters.
7. The method of claim 1, wherein the target data comprises porosity data, hydrate saturation data, and permeability data, and wherein determining an evolution law of hydrate reservoir permeability from the target data comprises:
and outputting a relation curve of the pore filling type hydrate saturation and the matrix size and the permeability according to the model parameters, the porosity data, the saturation data and the permeability data of the pore filling type hydrate so as to determine the evolution rule of the permeability of the hydrate reservoir.
8. The method of claim 7, wherein the hydrate reservoir permeability satisfies equation (1):
wherein k is permeability, Q is flow, μ is viscosity, L is seepage path length, A is flow channel cross-sectional area, P in For inlet pressure, P out Is the outlet pressure.
9. An apparatus for analyzing the permeability evolution law of a hydrate reservoir, comprising:
a memory configured to store instructions; and
a processor configured to invoke the instructions from the memory and when executing the instructions is capable of implementing the method for analysing the evolution law of permeability of a hydrate reservoir according to any one of claims 1 to 8.
10. A machine-readable storage medium having instructions stored thereon for causing a machine to perform the method for analyzing the permeability evolution law of a hydrate reservoir according to any one of claims 1 to 8.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310271933.3A CN116362031A (en) | 2023-03-16 | 2023-03-16 | Method and device for analyzing permeability evolution law of hydrate reservoir |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310271933.3A CN116362031A (en) | 2023-03-16 | 2023-03-16 | Method and device for analyzing permeability evolution law of hydrate reservoir |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN116362031A true CN116362031A (en) | 2023-06-30 |
Family
ID=86941237
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310271933.3A Pending CN116362031A (en) | 2023-03-16 | 2023-03-16 | Method and device for analyzing permeability evolution law of hydrate reservoir |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116362031A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116611371A (en) * | 2023-07-20 | 2023-08-18 | 浙江晶鸿精密机械制造有限公司 | Method for improving porous bearing performance and porous bearing |
-
2023
- 2023-03-16 CN CN202310271933.3A patent/CN116362031A/en active Pending
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116611371A (en) * | 2023-07-20 | 2023-08-18 | 浙江晶鸿精密机械制造有限公司 | Method for improving porous bearing performance and porous bearing |
| CN116611371B (en) * | 2023-07-20 | 2023-10-20 | 浙江晶鸿精密机械制造有限公司 | Method for improving porous bearing performance and porous bearing |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10671775B2 (en) | Sediment transport simulation with parameterized templates for depth profiling | |
| CN106837315B (en) | Method for representing coupling effect of fractured carbonate rock matrix and fractures | |
| CN106960070A (en) | A kind of seepage simulation method that coal body is reconstructed based on finite element discretization member CT | |
| WO2011163166A2 (en) | Machine, program product, and computer-implemented method to simulate reservoirs as 2.5d unstructured grids | |
| CN116362031A (en) | Method and device for analyzing permeability evolution law of hydrate reservoir | |
| US9589081B2 (en) | Karstification simulation | |
| Hilpert et al. | Calibration of a pore-network model by a pore-morphological analysis | |
| CN105260543A (en) | Double-pore model-based multi-medium oil gas flow simulation method and device | |
| CN106442271A (en) | Rock core permeability simulating method and rock core permeability simulating device | |
| CN108961409B (en) | Method for constructing 3D printing physical model based on oil reservoir three-dimensional geologic body | |
| CA2865866A1 (en) | Method for efficient dynamic gridding | |
| CN113361161B (en) | Shale oil-water two-phase flow simulation method and device based on digital core and storage medium | |
| CN110990926B (en) | Urban surface building hydrodynamic simulation method based on area correction rate | |
| Kerimov et al. | Mechanical trapping of particles in granular media | |
| CA2803144C (en) | Machine, program product, and computer-implemented method to simulate reservoirs as 2.5d unstructured grids | |
| CN107451383B (en) | Calibration method for initial bed sand gradation of plane two-dimensional water sand mathematical model | |
| CN111677486A (en) | Compact oil carbon dioxide huff and puff simulation method and device and storage medium | |
| CN114153007B (en) | Modeling method for water lock effect damage oil-gas layer, damage degree spatial-temporal evolution 4D quantitative and intelligent diagnosis method and system thereof | |
| CN115579074A (en) | Rock surface wettability prediction method and system | |
| Oxtoby et al. | A finite-volume solver for two-fluid flow in heterogeneous porous media based on OpenFOAM R G | |
| RU2626974C2 (en) | Method for determining a karstic region | |
| CN109063348B (en) | Displacement simulation method and device based on pore throat network model | |
| CN113281825A (en) | Rock physical model construction method and device | |
| Bak et al. | Porcolation: an invasion percolation model for mercury porosimetry | |
| CN115755610B (en) | Water injection throughput development numerical simulation system |
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 |