CN116386778B - Construction method, device, terminal and storage medium of pore adsorption density model - Google Patents
Construction method, device, terminal and storage medium of pore adsorption density model Download PDFInfo
- Publication number
- CN116386778B CN116386778B CN202310280912.8A CN202310280912A CN116386778B CN 116386778 B CN116386778 B CN 116386778B CN 202310280912 A CN202310280912 A CN 202310280912A CN 116386778 B CN116386778 B CN 116386778B
- Authority
- CN
- China
- Prior art keywords
- pore
- adsorption density
- wall material
- effective
- different
- 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.)
- Active
Links
- 239000011148 porous material Substances 0.000 title claims abstract description 457
- 238000001179 sorption measurement Methods 0.000 title claims abstract description 178
- 238000003860 storage Methods 0.000 title claims abstract description 15
- 238000010276 construction Methods 0.000 title claims description 12
- 239000000463 material Substances 0.000 claims abstract description 178
- 238000000034 method Methods 0.000 claims abstract description 40
- 238000013507 mapping Methods 0.000 claims abstract description 28
- 238000004364 calculation method Methods 0.000 claims abstract description 15
- 239000001307 helium Substances 0.000 claims abstract description 13
- 229910052734 helium Inorganic materials 0.000 claims abstract description 13
- SWQJXJOGLNCZEY-UHFFFAOYSA-N helium atom Chemical compound [He] SWQJXJOGLNCZEY-UHFFFAOYSA-N 0.000 claims abstract description 13
- 238000012900 molecular simulation Methods 0.000 claims abstract description 13
- 238000005516 engineering process Methods 0.000 claims abstract description 11
- 238000004590 computer program Methods 0.000 claims description 21
- 239000005416 organic matter Substances 0.000 claims description 18
- 238000011010 flushing procedure Methods 0.000 claims description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 16
- 239000007789 gas Substances 0.000 description 13
- 229910052500 inorganic mineral Inorganic materials 0.000 description 8
- 239000011707 mineral Substances 0.000 description 8
- 238000010586 diagram Methods 0.000 description 7
- 230000008569 process Effects 0.000 description 7
- 238000004088 simulation Methods 0.000 description 7
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 6
- 230000006870 function Effects 0.000 description 6
- 229910021389 graphene Inorganic materials 0.000 description 6
- 230000003993 interaction Effects 0.000 description 6
- 238000005381 potential energy Methods 0.000 description 6
- 238000011161 development Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 229910021532 Calcite Inorganic materials 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 239000012530 fluid Substances 0.000 description 3
- 229910052900 illite Inorganic materials 0.000 description 3
- 238000013113 molecular simulation experiment Methods 0.000 description 3
- VGIBGUSAECPPNB-UHFFFAOYSA-L nonaaluminum;magnesium;tripotassium;1,3-dioxido-2,4,5-trioxa-1,3-disilabicyclo[1.1.1]pentane;iron(2+);oxygen(2-);fluoride;hydroxide Chemical compound [OH-].[O-2].[O-2].[O-2].[O-2].[O-2].[F-].[Mg+2].[Al+3].[Al+3].[Al+3].[Al+3].[Al+3].[Al+3].[Al+3].[Al+3].[Al+3].[K+].[K+].[K+].[Fe+2].O1[Si]2([O-])O[Si]1([O-])O2.O1[Si]2([O-])O[Si]1([O-])O2.O1[Si]2([O-])O[Si]1([O-])O2.O1[Si]2([O-])O[Si]1([O-])O2.O1[Si]2([O-])O[Si]1([O-])O2.O1[Si]2([O-])O[Si]1([O-])O2.O1[Si]2([O-])O[Si]1([O-])O2 VGIBGUSAECPPNB-UHFFFAOYSA-L 0.000 description 3
- 239000010453 quartz Substances 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 3
- 239000000126 substance Substances 0.000 description 3
- 238000000342 Monte Carlo simulation Methods 0.000 description 2
- 239000004927 clay Substances 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000013508 migration Methods 0.000 description 2
- 230000005012 migration Effects 0.000 description 2
- 238000000329 molecular dynamics simulation Methods 0.000 description 2
- 229910052683 pyrite Inorganic materials 0.000 description 2
- 239000011028 pyrite Substances 0.000 description 2
- NIFIFKQPDTWWGU-UHFFFAOYSA-N pyrite Chemical compound [Fe+2].[S-][S-] NIFIFKQPDTWWGU-UHFFFAOYSA-N 0.000 description 2
- 229910004298 SiO 2 Inorganic materials 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000004090 dissolution Methods 0.000 description 1
- 210000000003 hoof Anatomy 0.000 description 1
- 238000011065 in-situ storage Methods 0.000 description 1
- 150000002500 ions Chemical class 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 229910044991 metal oxide Inorganic materials 0.000 description 1
- 150000004706 metal oxides Chemical class 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000000302 molecular modelling Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000011368 organic material Substances 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C60/00—Computational materials science, i.e. ICT specially adapted for investigating the physical or chemical properties of materials or phenomena associated with their design, synthesis, processing, characterisation or utilisation
-
- 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
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
Abstract
The application provides a method and a device for constructing a pore adsorption density model, a terminal and a storage medium. The method comprises the following steps: calculating corresponding pore adsorption density according to defined pore diameters of different pore wall materials by utilizing a molecular simulation technology, and defining the pore adsorption density as a first pore adsorption density; obtaining effective apertures obtained by calibrating different pore wall materials based on standard helium, and obtaining the mapping relation between the defined aperture and the effective aperture of each pore wall material; according to the defined pore diameter, the first pore adsorption density and the effective pore diameter corresponding to each pore wall material, obtaining the pore adsorption density of the corresponding pore wall material, and defining the pore adsorption density as a second pore adsorption density; and constructing pore adsorption density models of different pore wall materials according to the effective pore diameters and the second pore adsorption densities corresponding to the different pore wall materials. The application can realize the calculation of the pore adsorption density corresponding to different pore wall materials.
Description
Technical Field
The application relates to the technical field of geological reservoirs, in particular to a method, a device, a terminal and a storage medium for constructing a pore adsorption density model.
Background
Shale gas in unconventional oil gas is mainly stored in low-pore ultralow-permeability organic matter-rich shale in adsorption, free and dissolution states, and the in-situ storage mode of self-generation, self-storage and self-preservation with short migration distance and large-scale aggregation is an important characteristic of the shale gas. The pores and cracks in shale are used as main reservoir and migration space in shale reservoirs, and are always the research focus in shale gas exploration and development. As the principal force of shale gas exploration and development at present, the pores and cracks of the sea-phase shale are various in types and different in pore size, and obvious pores or cracks with different pore wall materials or crack wall materials exist at the edges of organic matters or minerals and clay. The influence of pores or cracks with different pore diameters or crack widths on the adsorption capacity of shale gas is explored, the influence of the pores or cracks on the adsorption capacity of shale gas under different pore wall or crack wall materials is beneficial to deep understanding of the gas-containing mechanism of the shale gas, and the exploration and development progress of shale gas resources/reserves in China are quickened by scientific evaluation.
At present, research on development type, influence and distribution rule of shale pore gaps in sea, sea-land transition or land is mature. Numerous studies have also been made on the impact of pores or cracks on methane adsorption, including physical experiments and molecular modeling methods. However, physical experiments are performed on a certain volume of sample, and the obtained result represents the gas content of the whole sample. In shale reservoirs, different pores and cracks of the wall or pore wall materials are often observed at the edges of organic matters or minerals and clay, and research on methane adsorption amount of the pores and cracks due to the difference of the pore wall or the wall materials is lacking.
Disclosure of Invention
The application provides a method, a device, a terminal and a storage medium for constructing a pore adsorption density model, which are used for solving the problem that the prior art lacks in calculating the pore adsorption density of a mixed pore wall.
In a first aspect, the present application provides a method for constructing a pore adsorption density model, including:
calculating corresponding pore adsorption density according to defined pore diameters of different pore wall materials by utilizing a molecular simulation technology, and defining the pore adsorption density as first pore adsorption density, wherein the defined pore diameters are preset pore diameters corresponding to the different pore wall materials;
obtaining effective apertures obtained by calibrating different pore wall materials based on standard helium, and obtaining the mapping relation between the defined aperture and the effective aperture of each pore wall material, wherein the effective apertures are real pore apertures corresponding to the different pore wall materials;
according to the defined pore diameter, the first pore adsorption density and the effective pore diameter corresponding to each pore wall material, obtaining the pore adsorption density of the corresponding pore wall material, and defining the pore adsorption density as a second pore adsorption density;
and constructing a pore adsorption density model of different pore wall materials according to the effective pore diameter and the second pore adsorption density corresponding to the different pore wall materials.
In a second aspect, the present application provides a device for constructing a pore adsorption density model, including:
the first calculation module is used for calculating corresponding pore adsorption density according to defined pore diameters of different pore wall materials by utilizing a molecular simulation technology, defining the pore adsorption density as first pore adsorption density, wherein the defined pore diameters are preset pore diameters corresponding to the different pore wall materials;
the mapping module is used for obtaining effective apertures obtained by calibrating different pore wall materials based on standard helium, and obtaining the mapping relation between the defined aperture and the effective aperture of each pore wall material, wherein the effective aperture is a real pore aperture corresponding to the different pore wall materials;
the second calculation module is used for obtaining the pore adsorption density of the corresponding pore wall material according to the defined pore diameter, the first pore adsorption density and the effective pore diameter corresponding to each pore wall material, and defining the pore adsorption density as a second pore adsorption density;
the construction module is used for constructing pore adsorption density models of different pore wall materials according to the effective pore diameters and the second pore adsorption densities corresponding to the different pore wall materials.
In a third aspect, the present application provides a terminal comprising a memory, a processor and a computer program stored in the memory and executable on the processor, the processor implementing the steps of the method according to the first aspect or any one of the possible implementations of the first aspect when the computer program is executed.
In a fourth aspect, the present application provides a computer readable storage medium storing a computer program which when executed by a processor implements the steps of the method of the first aspect or any one of the possible implementations of the first aspect.
The application provides a construction method, a device, a terminal and a storage medium of a pore adsorption density model, which are characterized in that a molecular simulation technology is utilized to calculate first pore adsorption densities corresponding to defined pore diameters of different pore wall materials, then the mapping relation between the defined pore diameters and the effective pore diameters is determined through the obtained effective pore diameters, so that effective pore diameters and second pore adsorption densities corresponding to the different pore wall materials are obtained, and the pore adsorption density model of the different pore wall materials is constructed according to the effective pore diameters and the second pore adsorption densities. According to the method, the effective aperture can be determined more quickly and more conveniently according to the mapping relation between the defined aperture and the effective aperture, and the calculation of the pore adsorption density corresponding to different pore wall materials can be realized according to the pore adsorption density model of the different pore wall materials, wherein the pore adsorption density model has universality.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings that are needed in the embodiments or the description of the prior art will be briefly described below, it being obvious that the drawings in the following description are only some embodiments of the present application, and that other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flow chart of an implementation of a method for constructing a pore adsorption density model provided by an embodiment of the present application;
FIG. 2 is a graph showing the distribution of methane molecules in different pore wall materials according to an embodiment of the present application;
FIG. 3 is a schematic diagram of a pore adsorption density model for different pore wall materials provided by an embodiment of the present application;
FIG. 4 is a schematic diagram of pore size measurement of pyrite-organic matter provided in an embodiment of the present application;
FIG. 5 is a schematic diagram of pore diameter measurement of quartz-organic material according to an embodiment of the present application;
FIG. 6 is a schematic structural diagram of a device for constructing a pore adsorption density model according to an embodiment of the present application;
fig. 7 is a schematic diagram of a terminal according to an embodiment of the present application.
Detailed Description
In the following description, for purposes of explanation and not limitation, specific details are set forth such as the particular system architecture, techniques, etc., in order to provide a thorough understanding of the embodiments of the present application. It will be apparent, however, to one skilled in the art that the present application may be practiced in other embodiments that depart from these specific details. In other instances, detailed descriptions of well-known systems, devices, circuits, and methods are omitted so as not to obscure the description of the present application with unnecessary detail.
For the purpose of making the objects, technical solutions and advantages of the present application more apparent, the following description will be made by way of specific embodiments with reference to the accompanying drawings.
Fig. 1 is a flowchart of an implementation method of a pore adsorption density model according to an embodiment of the present application, which is described in detail below:
in step 101, a molecular simulation technique is used to calculate a corresponding pore adsorption density according to defined pore diameters of different pore wall materials, and the pore adsorption density is defined as a first pore adsorption density, and the defined pore diameter is a preset pore diameter corresponding to the different pore wall materials.
In the embodiment of the application, according to the defined pore diameters corresponding to different pore wall materials, a molecular simulation technology is adopted to evaluate and calculate different pore gap structure types, so as to obtain the first pore adsorption density.
The pore wall is the wall of the pore and the slit, so that the method is not only suitable for calculating the adsorption density of the pore, but also suitable for calculating the adsorption density of the slit, namely suitable for calculating the adsorption density of the pore.
In one possible implementation, the pore wall material may include illite-organic matter, quartz-organic matter, pyrite-organic matter, and calcite-organic matter.
In the embodiment of the application, graphene is used as a boundary material of pores and cracks instead of organic matters. Wherein the unit cell formula of illite is Si 2 AlO 5 (OH) quartz has a unit cell formula of SiO 2 The unit cell of pyrite is FeS 2 The unit cell of calcite is CaCO 3 The parameters of the various minerals are shown in table 1.
TABLE 1 mineral parameters
Wherein Graphene is Graphene, quartz is Quartz, calcite is Calcite, pyrrite is Pyrite, and Illite is Illite.
In the embodiment of the application, a molecular simulation experiment is carried out by utilizing a molecular simulation technology, wherein the molecular simulation experiment adopts a giant regular Monte Carlo method (Grand Canonical Monte Carlo, GCMC) to evaluate different pore structure types, and the methane adsorption densities of pore seams of which the pore wall materials are illite-organic matters, quartz-organic matters, pyrite-organic matters and calcite-organic matters are respectively simulated, so that theoretical support and evaluation parameters are provided for reservoir evaluation.
Specifically, the interaction between the fluid (gas) along the Z direction and the pore wall in the simulation is described by Steele potential in the simulation, as shown in formula (1):
wherein Δ= 0.335nm, ρ w ＝114nm -3 ，ε w ＝28K，σ w ＝0.3345nm，σ wf ＝(σ w +σ f )/2。
In the module simulation experiment, graphene is generally used as a boundary material of pores and cracks instead of organic matters. CLAYFF is a generic force field suitable for use in simulated hydrate and metal oxide systems where all atoms translate freely as point charges.
The trap pf force field is used to perform non-binding interactions based on pairwise additive Lennard Jones (LJ) potential energy as shown in equation (2):
U(r ij )＝4ε ij [(σ ij /r ij ) 12 -(σ ij /r ij ) 6 ] (2)
wherein r is ij Is the distance epsilon ij Is potential energy well depth sigma ij Is the potential energy.
The interactions between atoms i and j are calculated from the Lorentz-Berthelot combination and the interactions between minerals and methane are calculated from the coulomb potential as shown in equation (3):
E coul ＝e 2 /4лε o ∑ i≠j q i q j /r ij (3)
wherein E is coul Is coulomb potential energy, q i And q j Is a partial charge, e is the charge of an electron, ε o Is the dielectric constant of the vacuum medium, namely epsilon o ＝8.854187817*10 -12 The parameters of the LJ potential energy for F/m, specific minerals and graphene are shown in Table 2.
TABLE 2 LJ potential energy parameters for minerals and graphene
Short-range LJ interactions are truncated at 1.2nm in molecular dynamics simulations, long-range classical interactions and plate geometry are calculated from three-dimensional Ewald summation. In the simulation process, the pressure applied to the outer surface of the reservoir in the GCMC simulation frame is 30MPa, and the temperature is set to 130 ℃ so as to accord with the environment where most of the sea shale is located. At a given temperature and pressure, the fluid in the aperture and the fluid in the external reservoir are in chemical equilibrium. Photographs were taken of Visual Molecular Dynamics (VMD) in CPK model, with the size of each atom in the snapshot being smaller than the actual LJ size. At the beginning of the modeling, methane molecules are randomly distributed in the pores or gaps, and the distribution of methane molecules is shown in fig. 2. Adsorption was then carried out until equilibrium, and the molecular density was measured. During GCMC simulation, different ions in minerals are fixed, C1 molecules are randomly and equiprobability inserted or separated under the influence of chemical potential, 15 ten thousand MC cycles are contained in each adsorption molecule for achieving balance in simulation, and 50 ten thousand MC cycles are contained in each adsorption molecule for sampling density distribution. In addition, the chemical potential was obtained by the Widom insert method using monte carlo simulation and NVT ensemble in a tank, and the density of the mixture at a given pressure and temperature was calculated according to PREOS. And calculating the pore adsorption density of different mixed pore walls according to the defined pores under the given mixed pore walls.
In the embodiment of the application, the molecular simulation technology is utilized, the pore adsorption density of different pore wall materials corresponding to the defined pore diameters can be calculated according to the defined pore diameters corresponding to the different pore wall materials, and a data reference basis is provided for the subsequent effective pore diameter calculation and the pore adsorption density corresponding to the effective pore diameter calculation.
In step 102, obtaining effective apertures obtained by calibrating different pore wall materials based on standard helium, and obtaining a mapping relation between defined apertures and effective apertures of the pore wall materials, wherein the effective apertures are real pore apertures corresponding to the different pore wall materials.
In the embodiment of the application, the effective pore diameters corresponding to different pore wall materials, namely the real pore diameters corresponding to different pore wall materials, are obtained, and the mapping relation between the defined pore diameters corresponding to each pore wall material and the effective pore diameters is determined by utilizing the defined pore diameters corresponding to different pore wall materials preset in the step 101.
In one possible implementation, step 102 may include:
for any hole wall material, standard helium is flushed into the holes of the hole wall material, and the holes corresponding to the hole wall material are calibrated to obtain the effective hole diameter of the hole wall material;
and determining the mapping relation between the defined aperture and the effective aperture corresponding to the same pore wall material.
In the embodiment of the application, pores of different pore wall materials are measured by using standard gas helium, the effective pore diameter corresponding to each pore wall material is determined according to the amount of the standard helium which is flushed into the pores of each pore wall material, and then the defined pore diameter and the effective pore diameter corresponding to the same pore wall material are calculated to obtain the mapping relation between the defined pore diameter and the effective pore diameter corresponding to the same pore wall material.
In one possible implementation manner, after determining the mapping relationship between the defined pore diameter and the effective pore diameter corresponding to the same pore wall material, the method may further include:
and fitting the defined aperture and the effective aperture corresponding to any hole wall material based on a linear formula to obtain a defined-effective aperture model of the hole wall material.
For any hole wall material, after the mapping relation between the defined hole diameter and the effective hole diameter is determined, fitting the defined hole diameter and the effective hole diameter corresponding to the hole wall material according to Origin software by utilizing a linear formula, and determining a definition-effective hole diameter model corresponding to different hole wall materials according to fitting results, wherein the specific definition-effective hole diameter model is shown in a table 3, x is the defined hole diameter, and Y is the effective hole diameter.
TABLE 3 definition of different pore wall materials-effective pore size model
In the embodiment of the application, because the defined aperture is the aperture input in the software, and is not the real adsorption aperture of different pore wall materials in the experimental process, if the effective aperture is not calculated, the adsorption density calculation can generate deviation, so that a definition-effective aperture model corresponding to different pore wall materials is established, the methane adsorption density error caused by the unreal adsorption aperture of the defined aperture in the molecular simulation experiment can be corrected, the accuracy of the adsorption density calculation result is improved, and when the effective aperture corresponding to the pore wall material is required to be determined later, the effective aperture of the pore wall material can be directly obtained by utilizing the definition-effective aperture model according to the defined aperture corresponding to the pore wall material, and the experimental acquisition is not required.
In step 103, according to the defined pore diameter, the first pore adsorption density and the effective pore diameter corresponding to each pore wall material, the pore adsorption density of the corresponding pore wall material is obtained, and the pore adsorption density is defined as the second pore adsorption density.
According to the defined pore diameters corresponding to the different pore wall materials in the step 101 and the first pore adsorption density corresponding to the defined pore diameters, and according to the effective pore diameters corresponding to the defined pore diameters in the step 102, the true adsorption density of the pores of the different pore wall materials, namely the second pore adsorption density corresponding to the pore wall materials, is calculated.
In one possible implementation, step 103 may include:
inputting the defined pore diameter, the first pore adsorption density and the effective pore diameter of any pore wall material into an equal proportion formula, and calculating the second pore adsorption density of the pore wall material.
In the embodiment of the present application, for any pore wall material, the corresponding second pore adsorption density is calculated for the defined pore diameter, the first pore adsorption density and the effective pore diameter of the pore wall material by using an equal proportion formula, for example, assuming that the defined pore diameter of a certain pore wall material is y 1 The corresponding first pore adsorption density is ρ 1 The effective aperture is y 2 The second pore adsorption density ρ is then calculated according to an equal proportion formula 2 The calculation formula of (a) is as formula (4):
the second pore adsorption density can be calculated according to formula (4)
In step 104, a pore adsorption density model of the different pore wall materials is constructed according to the effective pore diameters and the second pore adsorption densities corresponding to the different pore wall materials.
And (3) establishing a pore adsorption density model of the different pore wall materials according to the effective pore diameters corresponding to the different pore wall materials in the step (102) and the second pore adsorption density corresponding to the effective pore diameters in the step (103). The pore adsorption density model is mainly used for calculating the adsorption density of pores of different mixed pore wall materials.
In one possible implementation, step 104 may include:
and fitting the effective pore diameter and the second pore adsorption density corresponding to any pore wall material based on a Boltzmann formula to obtain a pore adsorption density model of the pore wall material.
Wherein Boltzmann' S entropic equation, i.e. s= kln Ω. Where k is a boltzmann constant, S is a macroscopic coefficient entropy, and is a measure of the degree of disorder of molecular motion or arrangement, Ω is a possible microscopic number, and the larger Ω is, the more disordered the system is. From this, the microscopic meaning of entropy can be seen: entropy is a measure of the disorder of thermal motion of molecules within a system.
In the embodiment of the application, for any pore wall material, along with the increase of the effective pore diameter corresponding to the pore wall material, the second pore adsorption density corresponding to the effective pore diameter is obtained, and according to a large number of effective pore diameters and second pore adsorption densities of the pore wall material, the effective pore diameters and the second pore adsorption densities of the pore wall material are fitted through a boltzmann formula, so that a pore adsorption density model of the pore wall material is established, wherein in the process of establishing the pore adsorption density model, the effective pore diameters and the second pore adsorption densities are calculated according to defined pore diameters for different pore wall materials, including but not limited to illite-organic matters, quartz-organic matters, pyrite-organic matters and calcite-organic matters, so that the established pore adsorption density model is suitable for calculating the pore adsorption densities corresponding to different mixed pore wall materials. A specific model of the pore adsorption density of the different pore wall materials is shown in fig. 3.
The application provides a construction method of a pore adsorption density model, which utilizes a molecular simulation technology to calculate first pore adsorption densities corresponding to defined pore diameters of different pore wall materials, and then determines a mapping relation between the defined pore diameters and the effective pore diameters through the obtained effective pore diameters, so as to obtain effective pore diameters and second pore adsorption densities corresponding to the different pore wall materials, and constructs the pore adsorption density model of the different pore wall materials according to the effective pore diameters and the second pore adsorption densities. According to the method, the effective aperture can be determined more quickly and more conveniently according to the mapping relation between the defined aperture and the effective aperture, and the calculation of the pore adsorption density corresponding to different pore wall materials can be realized according to the pore adsorption density model of the different pore wall materials, wherein the pore adsorption density model has universality.
The method for constructing the pore adsorption density model described above will be described below by way of one embodiment example.
Firstly, pore diameter measurement is carried out on sea shale in a cow hoof pond group under different mixed pore wall materials, so that pyrite-organic matter effective pore diameter is 0.012 mu m, quartz-organic matter effective pore diameter is 0.021 mu m, and specific measurement results are shown in fig. 4 and 5.
Then, the effective pore size x=0.012 μm was brought into the pyrite-organic matter adsorption model in fig. 3, i.e., y=5.48+23421.67/(1+e) x+14.03/1.68 ) The adsorption density is 5.484g/nm 3 . The effective pore size x=0.021 μm was brought into the adsorption model of the quartz-organic matter in fig. 3, i.e. y=5.46+13942.7/(1+e) x+15.28/1.9 ) The adsorption density is 5.46g/nm 3 。
It should be understood that the sequence number of each step in the foregoing embodiment does not mean that the execution sequence of each process should be determined by the function and the internal logic, and should not limit the implementation process of the embodiment of the present application.
The following are device embodiments of the application, for details not described in detail therein, reference may be made to the corresponding method embodiments described above.
Fig. 6 is a schematic structural diagram of a device for constructing a pore adsorption density model according to an embodiment of the present application, and for convenience of explanation, only the portions relevant to the embodiment of the present application are shown, which is described in detail below:
as shown in fig. 6, the pore adsorption density model constructing apparatus 6 includes:
the first calculation module 61 is configured to calculate a corresponding pore adsorption density according to defined pore diameters of different pore wall materials by using a molecular simulation technique, define the pore adsorption density as a first pore adsorption density, and define a pore diameter as a preset pore diameter corresponding to the different pore wall materials;
the mapping module 62 is configured to obtain effective apertures obtained by calibrating different pore wall materials based on standard helium, and obtain a mapping relationship between a defined aperture and an effective aperture of each pore wall material, where the effective aperture is a real pore aperture corresponding to the different pore wall materials;
a second calculation module 63, configured to obtain a pore adsorption density of the pore wall material according to the defined pore diameter, the first pore adsorption density and the effective pore diameter corresponding to each pore wall material, and define the pore adsorption density as a second pore adsorption density;
the construction module 64 is configured to construct a pore adsorption density model of different pore wall materials according to the effective pore diameters and the second pore adsorption densities corresponding to the different pore wall materials.
The application provides a construction device of a pore adsorption density model, which utilizes a molecular simulation technology to calculate first pore adsorption densities corresponding to defined pore diameters of different pore wall materials, and then determines a mapping relation between the defined pore diameters and the effective pore diameters through the obtained effective pore diameters, so as to obtain effective pore diameters and second pore adsorption densities corresponding to the different pore wall materials, and constructs the pore adsorption density model of the different pore wall materials according to the effective pore diameters and the second pore adsorption densities. According to the method, the effective aperture can be determined more quickly and more conveniently according to the mapping relation between the defined aperture and the effective aperture, and the calculation of the pore adsorption density corresponding to different pore wall materials can be realized according to the pore adsorption density model of the different pore wall materials, wherein the pore adsorption density model has universality.
In one possible implementation manner, the mapping module may specifically include:
the calibration module is used for aiming at any hole wall material, flushing standard helium into the holes of the hole wall material, and calibrating the holes corresponding to the hole wall material to obtain the effective hole diameter of the hole wall material;
the determining module is used for determining the mapping relation between the defined aperture and the effective aperture corresponding to the same pore wall material.
In one possible implementation, the determining module may be further configured to:
and fitting the defined aperture and the effective aperture corresponding to any hole wall material based on a linear formula to obtain a defined-effective aperture model of the hole wall material.
In one possible implementation, the second computing module may be specifically configured to:
inputting the defined pore diameter, the first pore adsorption density and the effective pore diameter of any pore wall material into an equal proportion formula, and calculating the second pore adsorption density of the pore wall material.
In one possible implementation, the building block may specifically be configured to:
and fitting the effective pore diameter and the second pore adsorption density corresponding to any pore wall material based on a Boltzmann formula to obtain a pore adsorption density model of the pore wall material.
In one possible implementation, the pore wall material may include illite-organic matter, quartz-organic matter, pyrite-organic matter, and calcite-organic matter.
Fig. 7 is a schematic diagram of a terminal according to an embodiment of the present application. As shown in fig. 7, the terminal 7 of this embodiment includes: a processor 70, a memory 71, and a computer program 72 stored in the memory 71 and executable on the processor 70. The processor 70, when executing the computer program 72, implements the steps of the above-described embodiments of the method for constructing the pore adsorption density model, such as steps 101 to 104 shown in fig. 1. Alternatively, the processor 70, when executing the computer program 72, performs the functions of the modules of the apparatus embodiments described above, such as the functions of the modules 61 through 64 shown in fig. 6.
By way of example, the computer program 72 may be partitioned into one or more modules that are stored in the memory 71 and executed by the processor 70 to complete the present application. The one or more modules/units may be a series of computer program instruction segments capable of performing the specified functions, which instruction segments are used for describing the execution of the computer program 72 in the terminal 7. For example, the computer program 72 may be partitioned into modules 61 through 64 shown in fig. 6.
The terminal 7 may be a computing device such as a desktop computer, a notebook computer, a palm computer, a cloud server, etc. The terminal 7 may include, but is not limited to, a processor 70, a memory 71. It will be appreciated by those skilled in the art that fig. 7 is merely an example of the terminal 7 and is not limiting of the terminal 7, and may include more or fewer components than shown, or may combine some components, or different components, e.g., the terminal may further include input and output devices, network access devices, buses, etc.
The processor 70 may be a central processing unit (Central Processing Unit, CPU), or may be another general purpose processor, a digital signal processor (Digital Signal Processor, DSP), an application specific integrated circuit (Application Specific Integrated Circuit, ASIC), a Field-programmable gate array (Field-Programmable Gate Array, FPGA) or other programmable logic device, a discrete gate or transistor logic device, a discrete hardware component, or the like. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like.
The memory 71 may be an internal storage unit of the terminal 7, such as a hard disk or a memory of the terminal 7. The memory 71 may be an external storage device of the terminal 7, such as a plug-in hard disk, a Smart Media Card (SMC), a Secure Digital (SD) Card, a Flash memory Card (Flash Card) or the like, which are provided on the terminal 7. Further, the memory 71 may also include both an internal storage unit and an external storage device of the terminal 7. The memory 71 is used for storing the computer program as well as other programs and data required by the terminal. The memory 71 may also be used for temporarily storing data that has been output or is to be output.
It will be apparent to those skilled in the art that, for convenience and brevity of description, only the above-described division of the functional units and modules is illustrated, and in practical application, the above-described functional distribution may be performed by different functional units and modules according to needs, i.e. the internal structure of the apparatus is divided into different functional units or modules to perform all or part of the above-described functions. The functional units and modules in the embodiment may be integrated in one processing unit, or each unit may exist alone physically, or two or more units may be integrated in one unit, where the integrated units may be implemented in a form of hardware or a form of a software functional unit. In addition, the specific names of the functional units and modules are only for distinguishing from each other, and are not used for limiting the protection scope of the present application. The specific working process of the units and modules in the above system may refer to the corresponding process in the foregoing method embodiment, which is not described herein again.
In the foregoing embodiments, the descriptions of the embodiments are emphasized, and in part, not described or illustrated in any particular embodiment, reference is made to the related descriptions of other embodiments.
Those of ordinary skill in the art will appreciate that the various illustrative elements and algorithm steps described in connection with the embodiments disclosed herein may be implemented as electronic hardware, or combinations of computer software and electronic hardware. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the solution. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present application.
In the embodiments provided in the present application, it should be understood that the disclosed apparatus/terminal and method may be implemented in other manners. For example, the apparatus/terminal embodiments described above are merely illustrative, e.g., the division of the modules or units is merely a logical function division, and there may be additional divisions when actually implemented, e.g., multiple units or components may be combined or integrated into another system, or some features may be omitted or not performed. Alternatively, the coupling or direct coupling or communication connection shown or discussed may be an indirect coupling or communication connection via interfaces, devices or units, which may be in electrical, mechanical or other forms.
The units described as separate units may or may not be physically separate, and units shown as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.
In addition, each functional unit in the embodiments of the present application may be integrated in one processing unit, or each unit may exist alone physically, or two or more units may be integrated in one unit. The integrated units may be implemented in hardware or in software functional units.
The integrated modules/units, if implemented in the form of software functional units and sold or used as stand-alone products, may be stored in a computer readable storage medium. Based on such understanding, the present application may also be implemented by implementing all or part of the flow in the method of the above embodiment, or by instructing the relevant hardware by a computer program, where the computer program may be stored in a computer readable storage medium, and the computer program may be executed by a processor, where the steps of the method embodiment of constructing each pore adsorption density model are implemented. Wherein the computer program comprises computer program code which may be in source code form, object code form, executable file or some intermediate form etc. The computer readable medium may include: any entity or device capable of carrying the computer program code, a recording medium, a U disk, a removable hard disk, a magnetic disk, an optical disk, a computer Memory, a Read-Only Memory (ROM), a random access Memory (Random Access Memory, RAM), an electrical carrier signal, a telecommunications signal, a software distribution medium, and so forth. It should be noted that the computer readable medium may include content that is subject to appropriate increases and decreases as required by jurisdictions in which such content is subject to legislation and patent practice, such as in certain jurisdictions in which such content is not included as electrical carrier signals and telecommunication signals.
The above embodiments are only for illustrating the technical solution of the present application, and not for limiting the same; although the application has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical scheme described in the foregoing embodiments can be modified or some technical features thereof can be replaced by equivalents; such modifications and substitutions do not depart from the spirit and scope of the technical solutions of the embodiments of the present application, and are intended to be included in the scope of the present application.
Claims (6)
1. The construction method of the pore adsorption density model is characterized by comprising the following steps of:
calculating corresponding pore adsorption density according to defined pore diameters of different pore wall materials by utilizing a molecular simulation technology, and defining the pore adsorption density as first pore adsorption density, wherein the defined pore diameters are preset pore diameters corresponding to the different pore wall materials;
obtaining effective apertures obtained by calibrating different pore wall materials based on standard helium, and obtaining the mapping relation between the defined aperture and the effective aperture of each pore wall material, wherein the effective apertures are real pore apertures corresponding to the different pore wall materials;
according to the defined pore diameter, the first pore adsorption density and the effective pore diameter corresponding to each pore wall material, obtaining the pore adsorption density of the corresponding pore wall material, and defining the pore adsorption density as a second pore adsorption density;
constructing pore adsorption density models of different pore wall materials according to the effective pore diameters and the second pore adsorption densities corresponding to the different pore wall materials;
the method for obtaining the mapping relation between the defined aperture and the effective aperture of each pore wall material comprises the following steps:
aiming at any hole wall material, the standard helium is flushed into the holes of the hole wall material, and the holes corresponding to the hole wall material are calibrated to obtain the effective hole diameter of the hole wall material;
determining the mapping relation between the defined aperture and the effective aperture corresponding to the same pore wall material;
the step of obtaining the pore adsorption density of the pore wall material according to the defined pore diameter, the first pore adsorption density and the effective pore diameter corresponding to each pore wall material, and defining the pore adsorption density as a second pore adsorption density, comprising:
inputting a defined pore diameter, a first pore adsorption density and an effective pore diameter of any pore wall material into an equal proportion formula, and calculating a second pore adsorption density of the pore wall material;
the construction of the pore adsorption density model of different pore wall materials according to the effective pore diameter and the second pore adsorption density corresponding to the different pore wall materials comprises the following steps:
and fitting the effective pore diameter and the second pore adsorption density corresponding to any pore wall material based on a Boltzmann formula to obtain a pore adsorption density model of the pore wall material.
2. The method for constructing a pore adsorption density model according to claim 1, wherein after determining the mapping relationship between the defined pore diameter and the effective pore diameter corresponding to the same pore wall material, the method further comprises:
and fitting the defined aperture and the effective aperture corresponding to any hole wall material based on a linear formula to obtain a defined-effective aperture model of the hole wall material.
3. The method of claim 1, wherein the pore wall material comprises illite-organic matter, quartz-organic matter, pyrite-organic matter, and calcite-organic matter.
4. A pore adsorption density model construction apparatus, comprising:
the first calculation module is used for calculating corresponding pore adsorption density according to defined pore diameters of different pore wall materials by utilizing a molecular simulation technology, defining the pore adsorption density as first pore adsorption density, wherein the defined pore diameters are preset pore diameters corresponding to the different pore wall materials;
the mapping module is used for obtaining effective apertures obtained by calibrating different pore wall materials based on standard helium, and obtaining the mapping relation between the defined aperture and the effective aperture of each pore wall material, wherein the effective aperture is a real pore aperture corresponding to the different pore wall materials;
the second calculation module is used for obtaining the pore adsorption density of the corresponding pore wall material according to the defined pore diameter, the first pore adsorption density and the effective pore diameter corresponding to each pore wall material, and defining the pore adsorption density as a second pore adsorption density;
the construction module is used for constructing pore adsorption density models of different pore wall materials according to the effective pore diameters and the second pore adsorption densities corresponding to the different pore wall materials;
wherein the mapping module comprises:
the calibration module is used for aiming at any hole wall material, flushing the standard helium into the holes of the hole wall material, and calibrating the holes corresponding to the hole wall material to obtain the effective hole diameter of the hole wall material;
the determining module is used for determining the mapping relation between the defined aperture and the effective aperture corresponding to the same pore wall material;
the second computing module includes:
inputting a defined pore diameter, a first pore adsorption density and an effective pore diameter of any pore wall material into an equal proportion formula, and calculating a second pore adsorption density of the pore wall material;
the construction module comprises:
and fitting the effective pore diameter and the second pore adsorption density corresponding to any pore wall material based on a Boltzmann formula to obtain a pore adsorption density model of the pore wall material.
5. A terminal comprising a memory, a processor and a computer program stored in the memory and executable on the processor, characterized in that the processor, when executing the computer program, realizes the steps of the method for constructing a pore adsorption density model according to any one of the preceding claims 1 to 3.
6. A computer-readable storage medium storing a computer program, characterized in that the computer program, when executed by a processor, implements the steps of the method of constructing a pore adsorption density model according to any one of claims 1 to 3.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310280912.8A CN116386778B (en) | 2023-03-21 | 2023-03-21 | Construction method, device, terminal and storage medium of pore adsorption density model |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310280912.8A CN116386778B (en) | 2023-03-21 | 2023-03-21 | Construction method, device, terminal and storage medium of pore adsorption density model |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116386778A CN116386778A (en) | 2023-07-04 |
| CN116386778B true CN116386778B (en) | 2023-09-29 |
Family
ID=86960785
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310280912.8A Active CN116386778B (en) | 2023-03-21 | 2023-03-21 | Construction method, device, terminal and storage medium of pore adsorption density model |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116386778B (en) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113705868A (en) * | 2021-08-16 | 2021-11-26 | 长江大学 | Shale gas state quantitative prediction method, system, equipment and storage medium |
| CN115408950A (en) * | 2021-05-26 | 2022-11-29 | 中国石油化工股份有限公司 | Shale gas flow simulation method, storage medium and equipment |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111912957B (en) * | 2020-08-24 | 2021-06-01 | 东北石油大学 | Shale oil quantity detection method and system based on occurrence state |
-
2023
- 2023-03-21 CN CN202310280912.8A patent/CN116386778B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115408950A (en) * | 2021-05-26 | 2022-11-29 | 中国石油化工股份有限公司 | Shale gas flow simulation method, storage medium and equipment |
| CN113705868A (en) * | 2021-08-16 | 2021-11-26 | 长江大学 | Shale gas state quantitative prediction method, system, equipment and storage medium |
Non-Patent Citations (1)
| Title |
|---|
| 利用GCMC分子模拟技术研究页岩气的吸附行为和机理;卢双舫;沈博健;许晨曦;陈国辉;刘可禹;薛庆忠;方志雄;何希鹏;;地球科学(05);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116386778A (en) | 2023-07-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Millan et al. | On the fractal scaling of soil data. Particle-size distributions | |
| CN111624147B (en) | Relative permeability measuring method and device for rock core | |
| Deuflhard et al. | Computational Molecular Dynamics: Challenges, Methods, Ideas: Proceeding of the 2nd International Symposium on Algorithms for Macromolecular Modelling, Berlin, May 21–24, 1997 | |
| Maier et al. | Simulation of flow through bead packs using the lattice Boltzmann method | |
| Maryška et al. | Numerical simulation of fracture flow with a mixed-hybrid FEM stochastic discrete fracture network model | |
| Li et al. | Stochastic analysis of unsaturated flow with probabilistic collocation method | |
| Matthews et al. | Measurement and simulation of the effect of compaction on the pore structure and saturated hydraulic conductivity of grassland and arable soil | |
| WO2017016655A1 (en) | Method of predicting parameters of a geological formation | |
| CN109471166A (en) | A kind of deep carbonate reservoirs shear wave prediction technique based on porosity type inverting | |
| EP3096252A2 (en) | Adaptive multiscale multi-fidelity reservoir simulation | |
| CN104101904B (en) | A kind of quick method for asking for formation shear speed | |
| Wheeler et al. | A multiscale preconditioner for stochastic mortar mixed finite elements | |
| CN105425280A (en) | Prediction method for mineral modulus and pore structure | |
| Shu et al. | The existence of solutions for impulsive fractional partial neutral differential equations | |
| Bayat et al. | Mathematical models for soil particle‐size distribution and their overall and fraction‐wise fitting to measurements | |
| Masihi et al. | Fast estimation of connectivity in fractured reservoirs using percolation theory | |
| CN103527172A (en) | Method for calculating water saturation through variable rock electivity coupling index | |
| CN116386778B (en) | Construction method, device, terminal and storage medium of pore adsorption density model | |
| CN113189674B (en) | Saturation estimation method and system for natural gas hydrate | |
| Kostorz et al. | An efficient and robust method for parameterized nonintrusive reduced‐order modeling | |
| Asensio et al. | Coupling stabilized finite element methods with finite difference time integration for advection–diffusion–reaction problems | |
| Hunt | Scale-dependent hydraulic conductivity in anisotropic media from dimensional cross-over | |
| Xie et al. | Variable dimension fractal-based conversion method between the nuclear magnetic resonance T2 spectrum and capillary pressure curve | |
| Mohd Aji et al. | Molecular simulation study on methane adsorption in amorphous shale structure | |
| CN115408950A (en) | Shale gas flow simulation method, storage medium and equipment |
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 |