CN109635410A - A kind of permeability discrete element simulation method based on pore system - Google Patents
A kind of permeability discrete element simulation method based on pore system Download PDFInfo
- Publication number
- CN109635410A CN109635410A CN201811487984.5A CN201811487984A CN109635410A CN 109635410 A CN109635410 A CN 109635410A CN 201811487984 A CN201811487984 A CN 201811487984A CN 109635410 A CN109635410 A CN 109635410A
- Authority
- CN
- China
- Prior art keywords
- permeability
- discrete element
- aperture
- discrete
- pore
- 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.)
- Granted
Links
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
Abstract
A kind of permeability discrete element simulation method based on pore system provided by the invention, comprising the following steps: the mean permeability and gas adsorption-desorption curve of reservoir locating for research object in measurement research area；According to the resulting gas absorption-desorption curve of step 1, porosity type is determined；Porosity type obtained in step 2 is simplified, discretization aperture model is obtained；Calculate the pore-size distribution elementary statistics parameter of discretization aperture model obtained in step 3；Discretization aperture module converter obtained in step 3 is permeability discrete element by the pore-size distribution elementary statistics parameter of the discretization aperture model obtained according to step 4；The permeability discrete element expansion that step 5 obtains is embedded in reservoir locating for research object, the permeability discrete element analysis of reservoir locating for research object in the research area is obtained；This method can preferably meet oil gas/geothermal exploration exploitation, particularly oil gas/geothermal exploitation actual demand in reservoir of the permeability there are significant difference.
Description
Technical field
It is the invention belongs to the numerical simulation technology in geological exploration and development technology, in particular to a kind of based on pore system
Permeability discrete element simulation method.
Background technique
Permeability is one of oil gas, underground heat, hydrogeological study important parameter, is either still existed in oil-gas migration
Hot-fluid transmitting aspect all plays a crucial role.Therefore, Permeability Distribution is Geologic modeling, is especially being related to fluid fortune
The fundamental of the Geologic modeling of aspect is moved, accurate setting becomes preferred oil gas exploration Favorable Areas, determines next step oil gas/ground
One important step of hot Distribution Characteristics etc..Only effectively determine the permeability spread rule of main oil and gas reservoir or geothermal reservoir
Rule, could establish reliable geologic assessment model, provide correct technological guidance for the evaluating reservoir in later period and production decision.
In existing method for numerical simulation, the most commonly used is actual measurement permeability mean value is based on, by reservoir permeability in model
It is set as single fixed value, and carries out numerical simulation work on this basis；But this is larger with practical geology situation gap, simulation
As a result it also greatly differs from each other with real background.Some scholars use random approach, assign reservoir each unit lattice with random number permeability
Value；But random distribution characteristic does not simultaneously have practically matter meaning, and the calculated numerical result geologic interpretation of institute is poor.
With being constantly progressive for pore-size distribution research means, can be known by means such as gas absorption-detachment assays
The pore system type and pore-size distribution feature of reservoir, and pore system further can be determined as both sides open-type, ink
Bottle-type, parallel slits type or unilateral slit-type.But pore-size distribution feature not yet applies to Permeability Distribution research, does not exist
It is applied in high resolution numerical simulation.
Therefore, it is necessary to form a set of permeability being distributed with reference to aperture distribution characteristics, more effective characterization reservoir permeability
Discrete element simulation new method, to meet the actual demand of the exploration early periods such as oil gas, underground heat.
This method proposes permeability discrete element simulation method based on pore system type for the first time at home, is embedded in aperture
Data have effectively depicted Permeability Parameters distribution characteristics, can be on the basis of gas absorption-detachment assays data further
Deepen the reservoir permeability regularity of distribution, improves the precision of numerical simulation result.
Summary of the invention
The permeability discrete element simulation method based on pore system that the purpose of the present invention is to provide a kind of, solves existing
Geologic modeling when foundation penetration value it is inaccurate, cause geologic interpretation poor.
In order to achieve the above object, the technical solution adopted by the present invention is that:
A kind of permeability discrete element simulation method based on pore system provided by the invention, comprising the following steps:
Step 1, the mean permeability and gas adsorption-desorption curve of reservoir locating for research object in area are studied in measurement；
Step 2, according to the resulting gas absorption-desorption curve of step 1, porosity type is determined；
Step 3, porosity type obtained in step 2 is simplified, obtains discretization aperture model；
Step 4, the pore-size distribution elementary statistics parameter of discretization aperture model obtained in step 3 is calculated；
Step 5, the pore-size distribution elementary statistics parameter of the discretization aperture model obtained according to step 4 will be in step 3
Resulting discretization aperture module converter is permeability discrete element；
Step 6, the permeability discrete element expansion that step 5 obtains is embedded in reservoir locating for research object, obtains this and grinds
Study carefully the permeability discrete element analysis of reservoir locating for research object in area.
Preferably, firstly, acquiring m sample of stratum different location locating for research object, different depth in the research area,
And each sample permeability is measured, it is recorded as k respectively1、k2、k3……km；
Secondly, calculating the permeability mean value of the correspondence reservoir using formula (1)
Preferably, it is by the specific method that original aperture model conversation is discretization aperture model:
Firstly, the different pore size aspect of model according to obtained in step 2, draws original aperture model；
Secondly, set the total length of original aperture model as d, unit nm；The original aperture model includes n discrete tuples
Unit number is closed, then length corresponding to each discrete element assembled unit is
Then, according to the following average pore size r for calculating i-th of discrete element assembled uniti:
Wherein, SiFor the practical area covered of i-th of discrete unit internal orifice dimension；
Finally, according to the aperture r of above-mentioned resulting various discrete member assembled uniti, splicing drafting discretization aperture combination
Model.
Preferably, the specific of the pore-size distribution elementary statistics parameter of discretization aperture model obtained in step 3 is calculated
Method is:
It include n discrete element assembled unit in discrete element aperture model obtained in setting procedure 2, wherein various discrete
Aperture in first assembled unit is respectively r1、r2、r3……rn, unit nm, then principal statistical parameter may be expressed as:
Wherein,Refer to the aperture mean value of the discrete element aperture model, unit nm；σ (r) refers to the discrete element aperture model
The standard deviation of pore radius, unit nm；Cv () refers to the coefficient of variation of the discrete element aperture model.
Preferably, the combination of permeability discrete element analysis is converted by the combination of discretization aperture model obtained in step 3
Specific method is:
Firstly, including a permeability discrete unit of n ' in setting permeability discrete element analysis combination, then the permeability is discrete
The permeability of each permeability discrete unit is represented sequentially as in meta-model combinationUnit mD；
Wherein, the mean permeability of permeability discrete element analysis combinationFollowing relationship should be met:
The coefficient of permeability variation of permeability discrete element analysis combinationFollowing relationship should be met:
The standard deviation of permeability discrete element analysis combinationThen further indicate that are as follows:
Secondly, supplementing specific conditional relationship according to different pore structures, wherein it is assumed that discrete element combination includes 2
Permeability discrete unit, then:
If pore system is both sides open-type, relational expression should be met:
If pore system is ink bottle-type, relational expression should be met:
If pore system is parallel slits type, relational expression should be met:
If pore system is unilateral slit-type, relational expression should be met:
It is assumed that discrete element combination includes 2+1 permeability discrete unit, then:
If pore system is both sides open-type, relational expression should be met:
If pore system is ink bottle-type, relational expression should be met:
If pore system is parallel slits type, relational expression should be met:
If pore system is unilateral slit-type, relational expression should be met:
Finally, then being solved in discrete element combination corresponding to each discrete element assembled unit by simultaneous formula (5)-(11)
Permeability value
Preferably, the permeability discrete element expansion that step 5 obtains is embedded in using horizontal repetition or longitudinal repetitive mode
In reservoir locating for research object.
Compared with prior art, the beneficial effects of the present invention are:
A kind of permeability discrete element simulation method based on pore system provided by the invention, gas absorption-desorption is real
It tests parameter to combine with permeability measured data, it is special that reservoir permeability distribution can be more reliably characterized according to pore system type
Sign can preferably meet oil gas/geothermal exploration exploitation, particularly oil gas/ground in reservoir of the permeability there are significant difference
The actual demand of heat exploitation.
Detailed description of the invention
Fig. 1 is the permeability discrete element simulation flow chart in the present invention；
Fig. 2 is the pore system type schematic diagram in the present invention；
Fig. 3 is that the hole discrete element in the present invention is converted into permeability discrete element schematic diagram.
Specific embodiment
With reference to the accompanying drawing, the present invention is described in more detail.
As shown in Figure 1, a kind of permeability discrete element simulation method based on pore system provided by the invention, including it is following
Step:
Step 1: determine research each stratum in area analog type, wherein the analog type include fining simulation and
Simplify simulation:
Firstly, determining the research object in research area, and the analog type on stratum locating for the research object is determined as essence
Refinement simulation；And the analog type on other stratum of research area is determined as to simplify simulation；
Step 2: measurement reservoir mean permeability:
Firstly, acquiring multiple samples of stratum different location, different depth locating for the research object, and measure each sample infiltration
Saturating rate, is recorded as k respectively1、k2、k3……km；
Secondly, calculating the permeability mean value of the correspondence reservoir using formula (1)
Step 3: obtaining gas absorption-desorption curve, porosity type is determined:
Firstly, sample in acquisition reservoir, and " rock specific surface and pore-size distribution measure static nitrogen adsorption capacity according to standard
Method (SY/T6154-1995) " carries out gas and (is often referred to N2、CO2Deng) adsorption-desorption experiment, the absorption-for obtaining counter sample is de-
Attached curve；
Secondly, according to the adsorption-desorption curve of acquisition and combining purely related to applied chemistry international federation (IUPAC)
Research achievement judges that the adsorption-desorption curve belongs to H1、H2、H3Or H4Type；
Finally, determining the porosity type of the reservoir, affiliated porosity type packet according to the affiliated type of adsorption-desorption curve
Include both sides open-type, ink bottle-type, parallel slits type and unilateral slit-type (Fig. 2).
Step 4: simplifying porosity type, discretization characterizes pore-size distribution:
In order to which influence of the pore system to reservoir permeability can be characterized in discrete element analysis, need first by original aperture model
It is converted into discretization aperture model, specifically:
Firstly, the different pore size aspect of model according to described in pure and applied chemistry international federation (IUPAC), draws
Original aperture model (Fig. 3), aperture model total length are d, unit nm；Secondly, set the total length of original aperture model as d,
Unit nm；The original aperture model includes n discrete element assembled unit number, then length corresponding to each discrete element assembled unit is
Then, according to the following average pore size r for calculating i-th of discrete element assembled uniti:
Wherein, SiFor i-th of practical area covered of discrete element assembled unit internal orifice dimension；
Finally, according to the aperture r of above-mentioned resulting various discrete member assembled uniti, splicing drafting discretization aperture combination
Model (Fig. 3)；
Step 5: calculated hole diameters are distributed elementary statistics parameter: in order to which pore-size distribution feature is reflected into permeability feature
On, need calculated hole diameters to be distributed principal statistical parameter first.It is assumed that aperture discrete element analysis group is amounted to including n aperture discrete element
Assembled unit, the aperture in each aperture discrete element assembled unit is respectively r1、r2、r3……rn, unit nm.Then principal statistical
Learning parameter may be expressed as:
In above formula,Refer to the aperture mean value in a discrete element combination, can reflect that reservoir average pore size is horizontal in model,
Unit nm；σ (r) refers to the standard deviation of pore radius, reflects the amplitude of variation in aperture, unit nm；And cv () refers to the variation in discrete element combination
Coefficient can reflect the dispersion degree in aperture under the premise of nondimensional.Above-mentioned each index can be the foundation of permeability discrete element
Reference frame is provided.
Step 6: converting permeability discrete element for aperture discrete element: it is assumed that each permeability discrete element group is amounted to comprising n
Permeability discrete unit, the permeability of permeability discrete unit is represented sequentially as in combinationDeng unit
mD.For guarantee permeability discrete element combination mean permeability it is consistent with the overall permeability of reservoir in model, then should meet with
Lower relational expression:
In formula,Represent the permeability mean value in discrete element combination, unit mD.
Meanwhile in order to enable permeability discrete element dispersion degree and pore-size distribution dispersion degree it is almost the same, then permeate
Rate discrete element should also meet formula:
In formula,Represent the coefficient of permeability variation in discrete element combination, dimensionless.Wherein, permeability discrete element
Standard deviationIt can further indicate that are as follows:
Unit mD.
In addition to this, since different pore systems has different pore-size distribution features, different holes need to be directed to
Gap structure supplements specific conditional relationship.It is assumed that discrete element combination includes 2 permeability discrete units, if pore system is both sides
Open-type (H1Type), then it should meet relational expression:
If pore system is ink bottle-type (H2Type), then it should meet relational expression:
If pore system is parallel slits type (H3Type), then it needs to meet relational expression:
If pore system is unilateral slit-type (H4Type), then it should ensure that:
It is assumed that discrete element combination includes 2+1 permeability discrete unit, if pore system is both sides open-type (H1Type), then
Relational expression should be met:
If pore system is ink bottle-type (H2Type), then it should meet relational expression:
If pore system is parallel slits type (H3Type), then it needs to meet relational expression:
If pore system is unilateral slit-type (H4Type), then it should ensure that:
By simultaneous formula (5)-(15), then permeability value corresponding to each discrete unit in discrete element combination can be solved Deng (Fig. 3).
Step 7: extrapolation permeability discrete element is to entire simulation system: as involved in the combination of permeability discrete element from
Throwaway member is limited, and discrete unit need to be extrapolated to entire reservoir model.Therefore, can be repeated by level, it is longitudinal the modes such as repeat,
The permeability discrete unit that can reflect pore system aperture spread is extended into entire fluid reservoir, to be embedded in entire simulation system
Among；The entire simulation system is the TOUGH groupware, which is situated between for simulating one-dimensional, the porous crack of two and three dimensions
In matter, the Numerical Simulation Program of fluid flowing, the heat transfer and Contaminants Transport of multiphase flow, multicomponent and Fei Dengwen.
In order to keep those skilled in the art discrete to the permeability in the invention patent based on pore system
First analogy method has understanding more profound, below will be by taking the area A as an example, referring to flow chart (Fig. 1), to the permeability in the present invention
Discrete element simulation method is described in further detail.
Step 1: geological system and corresponding reservoir that determination need to simulate:, first should be true before carrying out numerical simulation work
Surely the research area that need to simulate, research object and corresponding reservoir.The area A is the shale gas block of Chongqing southeast, studies area
Interior main force's shale gas exploration and development series of strata are Lower Cambrian Series Niutitang group, buried depth about 3000-4000m.According to B in the area A
Well shale gas is drilled well and is disclosed, and Niutitang group integral thickness about 125m, bottom target zone reservoir corresponds to depth 3802-3842m, excellent
Matter shale thickness about 40m, lithology with dark-grey-grey black mud stone, grey matter mud stone, containing silty, black silicalite and greyish black
Color contains based on siliceous mudstone, whole to have many characteristics, such as that the content of organic matter is abundant, evolution grade is higher, hole infiltration condition is extremely low.Therefore,
Can determine substantially: the area A is research area, and the High-quality Reservoir of Niutitang group bottom about 40m is the emphasis shale gas for needing fine analog
Reservoir, other stratum can then simplify processing up and down.
Step 2: measurement reservoir mean permeability: choosing 29 blocks of Niutitang group shale that B well 3802-3842m corresponds to depth
Sample carries out permeability analysis test, and corresponding target zone mean permeability is calculated using formula (1)Since research area is without other research floor positions, there is no need to calculate the flat of other shale gas reservoirs
Equal permeability.
1 B well 3802-3842m core permeability table of table
Step 3: obtaining gas absorption-desorption curve, porosity type is determined: choosing B well 3802-3842m sample and carry out N2
Adsorption-desorption experiment, Yi Faxian: the N of corresponding depth shale samples2Adsorption curve is almost overlapped with desorption curve, and adsorb-de-
Attached curve is whole more straight, according to IUPAC research achievement, indicates that adsorption-desorption curve belongs to H4Type, corresponding hole class
Type is unilateral slit-type.
Step 4: simplifying porosity type, discretization characterizes pore-size distribution: according to step 3, it may be determined that the area A Niutitang group
The porosity type of target zone is unilateral slit-type.Meanwhile according to research area's pore analysis, main micropore size is about 5nm.It is wrapping
In discrete element combination containing 3 discrete units, the aperture that can successively set discrete element is respectively r1=1nm, r2=5nm and r3=
9nm, it is suitable for the combinations of the discrete element aperture model of numerical simulation.
Step 5: calculated hole diameters be distributed elementary statistics parameter: according to it is above-mentioned include 3 discrete units aperture combination,
It can calculated hole diameters mean valueAre as follows:
The standard deviation of pore radius σ (r) are as follows:
Meanwhile being based on aperture mean valueIt, can calculated hole diameters coefficient of variation cv (r) with the standard deviation of pore radius σ (r), it may be assumed that
Above-mentioned each index can provide reference frame for the foundation of permeability discrete element.
Step 6: converting permeability discrete element for aperture discrete element: in the discrete element combination comprising 3 discrete units
In, it is assumed that the permeability of each discrete unit is followed successively by φ1、φ2、φ3Deng unit mD.Permeability to guarantee discrete element combination is equal
Value is consistent with overall permeability in model reservoir, should meet relational expression:
Meanwhile in order to enable permeability discrete element dispersion degree and pore-size distribution dispersion degree it is almost the same, then permeate
Rate discrete element should also meet formula:
In formula, cv (φ) represents the coefficient of permeability variation in discrete element combination, dimensionless.Wherein, permeability discrete element
Standard deviation sigma (φ) can further indicate that are as follows:
Unit mD.In addition to this, due to the N for the southeast Niutitang group target zone that changes2Adsorption-desorption curve belongs to H4Class
Type, corresponding pore system is unilateral slit-type, therefore need to meet relational expression:
By each formula of simultaneous step 6, can solve:
I.e. discrete element combination in each discrete unit permeability successively are as follows: 0.00051mD, 0.00257mD and
0.00463mD。
Step 7: extrapolation permeability discrete element is to entire simulation system: generally speaking, research area stratum is in horizontal distribution,
Therefore discrete element combination can be extended into the entire area A Niutitang group target zone, and be embedded in entire by horizontal duplicate mode
In discrete element analysis.Remaining stratum then can be by uniforming permeability Unified Characterization.
It will be understood by those of skill in the art that the accurate measurement of pore-size distribution and Permeability Parameters is permeability discrete element
The important prerequisite of simulation, and the representativeness of actual sample, reliability can influence the accuracy of later period numerical simulation.Therefore, it is
Guarantee that the calculated result of this method can effectively be compared with practical geological condition, it is necessary to screen to have before numerical simulation
Representational sample, penetration value analog result ability precision with higher.
Claims (6)
1. a kind of permeability discrete element simulation method based on pore system, which comprises the following steps:
Step 1, the mean permeability and gas adsorption-desorption curve of reservoir locating for research object in area are studied in measurement；
Step 2, according to the resulting gas absorption-desorption curve of step 1, porosity type is determined；
Step 3, porosity type obtained in step 2 is simplified, obtains discretization aperture model；
Step 4, the pore-size distribution elementary statistics parameter of discretization aperture model obtained in step 3 is calculated；
Step 5, the pore-size distribution elementary statistics parameter of the discretization aperture model obtained according to step 4 by step 3 gained
Discretization aperture module converter be permeability discrete element；
Step 6, the permeability discrete element expansion that step 5 obtains is embedded in reservoir locating for research object, obtains the research area
The permeability discrete element analysis of reservoir locating for interior research object.
2. a kind of permeability discrete element simulation method based on pore system according to claim 1, which is characterized in that step
In rapid 1, firstly, acquiring m sample of stratum different location locating for research object, different depth in the research area, and measure each
Sample permeability, is recorded as k respectively1、k2、k3......km；
Secondly, calculating the permeability mean value of the correspondence reservoir using formula (1)
3. a kind of permeability discrete element simulation method based on pore system according to claim 1, which is characterized in that step
In rapid 3, it is by the specific method that original aperture model conversation is discretization aperture model:
Firstly, the different pore size aspect of model according to obtained in step 2, draws original aperture model；
Secondly, set the total length of original aperture model as d, unit nm；The original aperture model includes that n discrete element combination is single
First number, then length corresponding to each discrete element assembled unit be
Then, according to the following average pore size r for calculating i-th of discrete element assembled uniti:
Wherein, SiFor the practical area covered of i-th of discrete unit internal orifice dimension；
Finally, according to the aperture r of above-mentioned resulting various discrete member assembled uniti, splicing drafting discretization aperture combination model.
4. a kind of permeability discrete element simulation method based on pore system according to claim 1, which is characterized in that step
In rapid 4, the specific method for calculating the pore-size distribution elementary statistics parameter of discretization aperture model obtained in step 3 is:
It include n discrete element assembled unit in discrete element aperture model obtained in setting procedure 2, wherein various discrete tuple
The aperture closed in unit is respectively r1、r2、r3……rn, unit nm, then principal statistical parameter may be expressed as:
Wherein,Refer to the aperture mean value of the discrete element aperture model, unit nm；σ (r) refers to the hole of the discrete element aperture model
Diameter standard deviation, unit nm；Cv (r) refers to the coefficient of variation of the discrete element aperture model.
5. a kind of permeability discrete element simulation method based on pore system according to claim 1, which is characterized in that step
In rapid 5, it converts the combination of discretization aperture model obtained in step 3 to the specific method of permeability discrete element analysis combination
It is:
Firstly, including a permeability discrete unit of n ', then the permeability discrete element mould in setting permeability discrete element analysis combination
The permeability of each permeability discrete unit is represented sequentially as in type combinationUnit mD；
Wherein, the mean permeability of permeability discrete element analysis combinationFollowing relationship should be met:
The coefficient of permeability variation of permeability discrete element analysis combinationFollowing relationship should be met:
The standard deviation of permeability discrete element analysis combinationThen further indicate that are as follows:
Secondly, supplementing specific conditional relationship according to different pore structures, wherein it is assumed that discrete element combination includes 2n infiltration
Rate discrete unit, then:
If pore system is both sides open-type, relational expression should be met:
If pore system is ink bottle-type, relational expression should be met:
If pore system is parallel slits type, relational expression should be met:
If pore system is unilateral slit-type, relational expression should be met:
It is assumed that discrete element combination includes 2n+1 permeability discrete unit, then:
If pore system is both sides open-type, relational expression should be met:
If pore system is ink bottle-type, relational expression should be met:
If pore system is parallel slits type, relational expression should be met:
If pore system is unilateral slit-type, relational expression should be met:
Finally, then solving infiltration corresponding to each discrete element assembled unit in discrete element combination by simultaneous formula (5)-(11)
Saturating rate value
6. a kind of permeability discrete element simulation method based on pore system according to claim 1, which is characterized in that step
In rapid 6, the permeability discrete element expansion that step 5 obtains is embedded in by research object institute using horizontal repetition or longitudinal repetitive mode
Locate in reservoir.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811487984.5A CN109635410B (en) | 2018-12-06 | 2018-12-06 | Permeability discrete element simulation method based on pore system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811487984.5A CN109635410B (en) | 2018-12-06 | 2018-12-06 | Permeability discrete element simulation method based on pore system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109635410A true CN109635410A (en) | 2019-04-16 |
| CN109635410B CN109635410B (en) | 2021-08-17 |
Family
ID=66071664
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811487984.5A Active CN109635410B (en) | 2018-12-06 | 2018-12-06 | Permeability discrete element simulation method based on pore system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109635410B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110021220A (en) * | 2019-04-28 | 2019-07-16 | 中国华能集团清洁能源技术研究院有限公司 | A kind of building block system geothermal tail water recharge analysis system and application method |
| CN113189305A (en) * | 2021-04-29 | 2021-07-30 | 中国华能集团清洁能源技术研究院有限公司 | Simulation method and system for porous permeable rock |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2090907A1 (en) * | 2008-02-14 | 2009-08-19 | Exxonmobil Upstream Research Company | Method for determining the properties of hydrocarbon reservoirs from geophysical data |
| CN102339339A (en) * | 2010-07-23 | 2012-02-01 | 中国石油化工股份有限公司 | Method for analyzing remaining oil distribution of fractured-vuggy reservoir |
| CN103339488A (en) * | 2011-01-27 | 2013-10-02 | 普拉德研究及开发股份有限公司 | Gas sorption analysis of unconventional rock samples |
| CN104990851A (en) * | 2015-06-23 | 2015-10-21 | 西南石油大学 | Novel experimental research method for shale sensitivity |
| US20150355158A1 (en) * | 2014-06-05 | 2015-12-10 | Geocosm, LLC | Predicting sediment and sedimentary rock properties |
| US20170053046A1 (en) * | 2015-08-17 | 2017-02-23 | Saudi Arabian Oil Company | Capillary pressure analysis for petrophysical statistical modeling |
| CN107885893A (en) * | 2016-09-29 | 2018-04-06 | 中国石油化工股份有限公司 | The construction method of carbon dioxide flooding reservoir heterogeneity is described |
| CN108843300A (en) * | 2018-06-25 | 2018-11-20 | 中国石油天然气股份有限公司 | It is a kind of complexity porous media in determine mainstream channel type method and device |
-
2018
- 2018-12-06 CN CN201811487984.5A patent/CN109635410B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2090907A1 (en) * | 2008-02-14 | 2009-08-19 | Exxonmobil Upstream Research Company | Method for determining the properties of hydrocarbon reservoirs from geophysical data |
| CN102339339A (en) * | 2010-07-23 | 2012-02-01 | 中国石油化工股份有限公司 | Method for analyzing remaining oil distribution of fractured-vuggy reservoir |
| CN103339488A (en) * | 2011-01-27 | 2013-10-02 | 普拉德研究及开发股份有限公司 | Gas sorption analysis of unconventional rock samples |
| US20150355158A1 (en) * | 2014-06-05 | 2015-12-10 | Geocosm, LLC | Predicting sediment and sedimentary rock properties |
| CN104990851A (en) * | 2015-06-23 | 2015-10-21 | 西南石油大学 | Novel experimental research method for shale sensitivity |
| US20170053046A1 (en) * | 2015-08-17 | 2017-02-23 | Saudi Arabian Oil Company | Capillary pressure analysis for petrophysical statistical modeling |
| CN107885893A (en) * | 2016-09-29 | 2018-04-06 | 中国石油化工股份有限公司 | The construction method of carbon dioxide flooding reservoir heterogeneity is described |
| CN108843300A (en) * | 2018-06-25 | 2018-11-20 | 中国石油天然气股份有限公司 | It is a kind of complexity porous media in determine mainstream channel type method and device |
Non-Patent Citations (4)
| Title |
|---|
| H. DONG: "Characterization and comparison of capillary pore structures of digital cement pastes", 《MATERIALS AND STRUCTURES》 * |
| LAWRENCE M. ANOVITZ: "Characterization and Analysis of Porosity and Pore Structures", 《REVIEWS IN MINERALOGY & GEOCHEMISTRY》 * |
| 冯其红: "基于嵌入离散裂缝的页岩气藏视渗透率模型", 《地球科学》 * |
| 杨峰: "页岩储层微观孔隙结构特征", 《石油学报》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110021220A (en) * | 2019-04-28 | 2019-07-16 | 中国华能集团清洁能源技术研究院有限公司 | A kind of building block system geothermal tail water recharge analysis system and application method |
| CN110021220B (en) * | 2019-04-28 | 2024-03-26 | 中国华能集团清洁能源技术研究院有限公司 | Building block type geothermal tail water recharging analysis system and application method |
| CN113189305A (en) * | 2021-04-29 | 2021-07-30 | 中国华能集团清洁能源技术研究院有限公司 | Simulation method and system for porous permeable rock |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109635410B (en) | 2021-08-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101487390B (en) | Archie mode method for confirming initial oil saturation of oil layer | |
| CN106468172A (en) | A kind of Oil in Super-low Permeability sandstone oil reservoir low-resistance reservoir log interpretation method | |
| CN105426612A (en) | Stratum component optimization determining method and device | |
| CN105651966A (en) | Shale oil and gas high-quality reservoir stratum evaluation method and parameter determination method | |
| CN103776981B (en) | A kind of new karst stage time method of discrimination | |
| CN103510946A (en) | Method for evaluating reservoir fluid property through gas logging data | |
| CN104806232B (en) | A kind of method for determining porosity lower limit of fracture | |
| CN105653815A (en) | Reservoir fluid distribution quantitative interpretation method based on rock physical model theory | |
| CN104199121A (en) | Shale gas pool construction and production favorable area comprehensive determining method | |
| CN103867197A (en) | Complex lithology natural gas reservoir interval transit time discriminating method | |
| CN104747183A (en) | Carbonate reservoir comprehensive classification method | |
| CN107102377B (en) | The method of quantitative forecast tight sand favorable oil/gas exploration area | |
| CN104750884A (en) | Quantitative evaluation method of shale oil and gas enrichment index on the basis of multi-factor nonlinear regression | |
| CN104047598A (en) | Method for predicating productivity of nonhomogeneity ancient karst carbonate reservoir | |
| CN106644736A (en) | Method and device for measuring brittleness index of reservoir rock | |
| CN105134185A (en) | Identification method for reservoir fluid properties | |
| CN102913240B (en) | Reservoir fluid recognizing method | |
| CN108982320A (en) | It is a kind of to carry out Complicated Pore Structures reservoir permeability calculation method using grain size parameter | |
| CN103867198A (en) | Method for distinguishing formation density of carbonatite natural gas reservoir | |
| CN109635410A (en) | A kind of permeability discrete element simulation method based on pore system | |
| CN110700820A (en) | Dessert classification method for compact oil reservoir in northern part of Songliao basin | |
| CN104751002A (en) | Method for determining effective sandstone reservoir | |
| CN104459790A (en) | Oil-gas possibility basin effective reservoir analysis method and device | |
| Bogatkov et al. | Fracture network modeling conditioned to pressure transient and tracer test dynamic data | |
| CN104948150B (en) | A kind of method and apparatus of definite stratum replacement pressure |
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 |