CN108829950A - A kind of unconventional reservoir permeability evaluation method based on core image - Google Patents
A kind of unconventional reservoir permeability evaluation method based on core image Download PDFInfo
- Publication number
- CN108829950A CN108829950A CN201810546608.2A CN201810546608A CN108829950A CN 108829950 A CN108829950 A CN 108829950A CN 201810546608 A CN201810546608 A CN 201810546608A CN 108829950 A CN108829950 A CN 108829950A
- Authority
- CN
- China
- Prior art keywords
- sub
- hole
- pore
- component characteristics
- unconventional reservoir
- 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
The present invention provides a kind of unconventional reservoir permeability evaluation method based on core image, includes the following steps：Respective image is obtained after true rock sample is scanned and is identified using rock imaging technique, rock sample in image is divided into corresponding feature；The frequency spectrum of every kind of feature is extracted respectively, is then fitted with mixing fractal method to analysis result, is obtained the fractal parameter of each individual features；Fractal parameter after each fitting is combined, the blowhole model based on mixing fractal theory is obtained；The assessment result of the apparent permeability of current unconventional reservoir can be obtained by calculating after connected probability matrix arrangement to the connected probability between variety classes hole in blowhole model by hole connectivity algorithm.The present invention is from true rock sample data, in conjunction with hole connectivity algorithm, has obtained the connected probability matrix between variety classes hole, and final realize rapidly and accurately assesses unconventional reservoir apparent permeability.
Description
Technical field
The invention belongs to field of petroleum exploitation, more particularly to a kind of can rapidly and accurately assess based on core image
The evaluation method of unconventional reservoir apparent permeability.
Background technique
Shale gas is the natural gas produced from shale reservoir, and ingredient is a kind of important emerging energy based on methane
Source.However successfully exploitation shale gas needs to solve many problems, underlying issue therein first is that rapidly and efficiently assessing page
The apparent permeability of rock.
The many because being known as of shale gas apparent permeability are influenced, mainly there is the following：First, shale is very fine and close, hole
Porosity and permeability are low, and seepage flow space is small, and main pore diameter is in nanometer scale；Second, seepage flow space style is various, connectivity
Poor, scale is different, including the organic hole of nanoscale, the inorganic hole in micron order crack；Third, gas flow and supporting environment are complicated,
(it is present in the form of free gas in blowhole and crack) including free state, ADSORPTION STATE (being adsorbed in blowhole wall surface) and molten
Solve state (being dissolved in kerogen, asphalitine and liquid crude oil) etc..Existing various technical solutions are to unconventional reservoir sample
Calculation amount is huge when computing permeability, needs to establish a kind of new method of rapid evaluation unconventional reservoir sample permeability.
The method in characterization blowhole space has capillary bunch method and interstitial space Reconstruction Method at present, and wherein capillary bunch method is table
Levy one of the classical way in blowhole space.Blowhole is reduced to series of parallel capillary bunch by this method, with research
Flow event of the fluid in high porosity rock flow event.This method model is simple, calculates easy.Develop on this basis
Addition tortuosity, change hollow billet beam diameter and increase a series of equations such as sliding item, to simulate variety classes blowhole
Space and process fluid flow.But the shortcomings that capillary bunch method, is not consider void network structure complicated in rock, does not also examine
Filter hole larynx connectivity, the foundation of model are based on capillary bunch theoretical model, and the pore network composition of this class model has with true rock
Certain gap.
And interstitial space reconstruct rule characterizes blowhole space based on statistics, this method counts the hole in rock first
Gap spread spectrum and ligancy.Further establish the network of a netted unicom pipeline composition, pipeline A venturi, two pipes
The intersection point in road simulates hole, and the aperture that statistics obtains is dispersed into proportion on the intersection point of network, then by the ligancy of statistics
Random blocking pipe is finally completed the reconstruct of interstitial space.The disadvantages of this method is to be reduced to mark by complicated interstitial space
After quasi- square mesh structural model, break up rear random dispersion network access network, eliminate hole and venturi be spatially distributed it is only
Vertical property.In addition, the tortuosity of venturi and the pore-size distribution of venturi itself also do not take in.
Summary of the invention
The object of the present invention is to provide and a kind of unconventional reservoir table can rapidly and accurately be assessed based on core image
See the evaluation method of permeability.
Particularly, the present invention provides a kind of unconventional reservoir permeability evaluation method based on core image, including as follows
Step：
Step 100, respective image is obtained after true rock sample being scanned and identified using rock imaging technique, by image
The organic matter block and ore particles Type division of middle rock sample are component characteristics, by organic matter block, ore particles and rock base
Hole information in matter is divided into pore character, and the other feature in rock sample is divided into supplemental characteristic；
Step 200, the frequency spectrum of every kind of feature is extracted respectively, is then fitted with mixing fractal method to analysis result,
Obtain the fractal parameter of each individual features；
Step 300, the fractal parameter after each fitting is combined, obtains the blowhole mould based on mixing fractal theory
Type；
Step 400, by hole connectivity algorithm to the connected probability between variety classes hole in blowhole model
By being calculated after connected probability matrix arrangement, the assessment result of the apparent permeability of current unconventional reservoir can be obtained.
In an embodiment of the invention, the pore character in the step 100 refers to organic matter block, mine
Organic matter hole in stone particle and the intergranular pore in ore or intragranular hole；The supplemental characteristic refers to natural in the rock sample
Crack or man-made fracture feature.
In an embodiment of the invention, the rock imaging technique includes SEM, FIB-SEM, Micro-CT,
nano-CT。
In an embodiment of the invention, the component characteristics are subdivided into more according to the type difference of ore particles
A corresponding sub- component characteristics；The pore character corresponds to the difference of the sub- component characteristics and is subdivided into multiple corresponding sub-aperture
Gap feature；The supplemental characteristic according to from the relationship of the different sub- component characteristics and to be subdivided into multiple corresponding sub- auxiliary special
Sign.
In an embodiment of the invention, the frequency spectrum of the component characteristics is built upon in wherein each sub- component characteristics
On the equivalent size distribution frequency spectrum of respective subgroup point, the frequency spectrum of the pore character is built upon the equivalent of each sub- pore character
On the spread spectrum in aperture, the frequency spectrum of the supplemental characteristic be built upon fracture spacing represented by sub- supplemental characteristic, seam it is long and
In slit width.
In an embodiment of the invention, every grade of subgroup that the fractal parameter of each sub- component characteristics obtains point
The size distribution frequency spectrum that the equivalent size of feature and its correspondence number collectively form, is equal to the size distribution of the sub- component characteristics
Frequency spectrum；The equivalent diameter for every grade of hole that the fractal parameter of each sub- pore character obtains and its correspondence number collectively form
Size distribution frequency spectrum, be equal to the size distribution frequency spectrum of the sub- pore character；Each sub- supplemental characteristic
The equivalent size for every grade of sub- supplemental characteristic that fractal parameter obtains and its correspond to the crack frequency spectrum that collectively forms of number, and be equal to
The crack frequency spectrum of sub- supplemental characteristic.
In an embodiment of the invention, the mode of connected probability matrix is established such as using interconnected pore probabilistic algorithm
Under：Each subcharacter in three kinds of features is sorted from large to small in the way of equivalent length of side or equivalent diameter, wherein first row
It is classified as number corresponding to equivalent length of side or equivalent diameter for equivalent length of side or equivalent diameter size, second, third column connection is general
Rate, the 4th is classified as the row component type.
In an embodiment of the invention, gathered and operated before sequence.
In an embodiment of the invention, the step of calculating permeability is as follows：
Wherein, the number of iterations i divides shape scale F, hole/block side length λ that i-th iteration generatesi, maximum pore/block
Side length λmax, hole/block number N for generating in i-th iterationi, hole/block number of iteration is not involved in i-th iteration
Ni-soild, A is the sectional area in rock core characterization unit face, CiFor the i-th row block body/pore shape factor, ki(pavg) it is in mean pressure
Power pavgThe permeability of lower i-th row block body/hole.
In an embodiment of the invention, ki(pavg) calculation formula need according to the type of row each in matrix into
Row judgement；
1. setting, component characteristics A, the sequence number consecutively of sub- component characteristics are respectively with A1, A2 ... An, pore character
For B, the sequence number consecutively of sub- pore character is with B1, B2 ... Bn, supplemental characteristic C, sub- supplemental characteristic with C1,
The sequence number consecutively of C2 ... Cn；In component characteristics A, each block for forming sub- component characteristics is porous media；When
In connected probability matrix when the i-th row the 4th column record behavior component characteristics A, ki(pavg) expression formula is：
Wherein, knD(j,pavg) it is average pressure pavgUnder, the permeability value of the sub- pore character of j-th stage；Calculation formula is as follows：
As j=1
As j=2
As j=3
And so on, as j=n
Wherein, λ1,λ2,…,λj,…,λnFor being comprised in porous media representated by the i-th row after being gathered and being operated
Hole equivalent diameter；N1,N2,…,Nj,…,NnTo be comprised in porous Jie representated by the i-th row after being gathered and being operated
The number of hole in matter；Wherein λ1≥λ2≥…≥λj≥…≥λn；
2. when the i-th row the 4th column record behavior is not by poly- hole simultaneously in connected probability matrix, permeability item ki(pavg)
For
3. as the i-th row the 4th column record behavior sub- supplemental characteristic C1 in connected probability matrix, permeability item ki(pavg)
For
4. formula (1) the mesostroma permeability item kmatrixFor
It is above it is various in, μ is gas viscosity, and M is number of moles of gas, and R is ideal gas constant, and T is temperature, ρavgIt is flat
Equal density, pavgFor average pressure, α is tangential momentum adjustment factor,For n kind matrix minerals percentage, k1,
k2,…,knFor the apparent permeability of n kind matrix minerals.
The present invention extracts main constituents, then based on mixing fractal theory, needle from true rock sample data
It to each component, is fitted using fractal parameter, is finally combined each component respectively, obtained based on mixing fractal theory
Blowhole model.Further combined with pore communication algorithm, the connected probability matrix between variety classes hole is obtained, most
A kind of unconventional reservoir permeability evaluation method is constructed eventually, is realized and is rapidly and accurately assessed unconventional reservoir apparent permeability.
Detailed description of the invention
Fig. 1 is one embodiment of the present invention evaluation method flow chart.
Specific embodiment
As shown in Figure 1, unconventional reservoir permeability evaluation method generality disclosed in one embodiment of the present invention includes
Following steps：
Step 100, respective image is obtained after true rock sample being scanned and identified using rock imaging technique, by image
The organic matter block and ore particles Type division of middle rock sample are component characteristics, by organic matter block, ore particles and rock base
Hole information in matter is divided into pore character, and the other feature in rock sample is divided into supplemental characteristic；
Here rock imaging technique includes but is not limited to SEM (scanning electron microscope), FIB-SEM (focused ion
Beam), Micro-CT (microcomputer layer scanning technology), the technologies such as nano-CT.
In above-mentioned each imaging technique, it can analyzed, can also be used after being imaged only with a kind of mode simultaneously
Two or more imaging techniques is combined analysis, because some imaging technique observation scopes are wide, some imaging technique precision are high,
It can integrate to obtain block size distribution and distribution of pores.
Rock can be shale, tight sand or other unconventional reservoir rock samples.
When dividing, by organic matter block, the pyrite block, quartz in rock sample image, feldspar, one kind mine such as montmorillonite
Stone type universal formulation is component characteristics, and can press sub- component characteristics A1, subgroup to the specific ore of each of the component characteristics
Dtex levies A2, the number of sub- component characteristics A3 ... is finely divided class；The holes such as intergranular pore, intragranular hole in image are all drawn
It is divided into pore character, wherein according to organic hole in the distributing position of hole such as organic mass body, hole in pyrite block,
Inorganic hole in matrix etc. is finely divided class by sub- pore character B1, sub- pore character B2, sub- pore character B3 ...；Auxiliary
It is characterized in the other feature in addition to aforementioned two kinds of features, such as crack, density feature, present embodiment is using crack as auxiliary
Feature is helped, crack can be divided into dry and man-made fracture according to type, will according to information such as fracture spacing, seam length, slit widths
Supplemental characteristic be divided into sub- supplemental characteristic C1, sub- supplemental characteristic C2, sub- supplemental characteristic C3 ....
Step 200, the frequency spectrum of every kind of feature is extracted respectively, is then fitted with mixing fractal method to analysis result,
Obtain the fractal parameter of each individual features；
It specifically includes：The scanning result of different scale rock core digital imaging technology based on use extracts component in image
The frequency spectrum of feature, pore character and FRACTURE CHARACTERISTICS, wherein component characteristics are organic plastid block, pyrite block etc. in the type
Size spread spectrum (being the equivalent size frequency spectrum of organic matter block by taking organic matter as an example), pore character is distribution of pores frequency
Spectrum (being organic hole effective aperture frequency spectrum in organic matter block by taking organic hole in organic mass as an example), supplemental characteristic is to split
The fracture spacing of seam stitches slit wide spectrum.
The block size distribution frequency spectrum of component characteristics is fitted using mixing fractal method, obtains fractal parameter AF1,
By fractal parameter AF1 each of obtain sub- component characteristics equivalent size and its correspond to the block size point that number collectively forms
Cloth frequency spectrum, the block size distribution spectral shape are equal to the block size distribution frequency spectrum of respective sub-set dtex sign, i.e. the two phase
It is same or close.
It reuses mixing fractal method to be fitted the pore diameter size distribution frequency spectrum of pore character, obtains a point parameter
The aperture that number BF1, the equivalent diameter of the sub- pore character of every grade for obtaining fractal parameter BF1 and its correspondence number collectively form is big
Small spread spectrum is equal to the size distribution frequency spectrum of corresponding sub- pore character, i.e. the two is identical or close.
Reuse mixing fractal method fracture feature fracture spacing, seam slit wide spectrum be fitted, obtain a point shape
The crack that parameter CF1, the equivalent size of the sub- supplemental characteristic of every grade for obtaining fractal parameter CF1 and its correspondence number collectively form
Frequency spectrum is equal to the crack frequency spectrum of the sub- FRACTURE CHARACTERISTICS, i.e., the two is identical or close.
Step 300, the fractal parameter after each fitting is combined, obtains the blowhole mould based on mixing fractal theory
Type；
In this step, component characteristics centainly include pore character, but pore character not necessarily includes by component characteristics, example
If the block size of organic mass body A1 itself is distributed, contain organic hole B1 inside each organic matter block A1 again, every piece has
Organic hole B1 in machine mass body A1 meets distribution of pores frequency spectrum representated by the above-mentioned organic hole of B1 class.
Step 400, by hole connectivity algorithm to the connected probability between variety classes hole in blowhole model
It is calculated after being arranged in a matrix, the assessment result of the apparent permeability of current unconventional reservoir can be obtained.
The mode for establishing matrix using interconnected pore probabilistic algorithm is as follows：By each subcharacter in three kinds of features according to equivalent
Side length or the mode of equivalent diameter sort from large to small, wherein first is classified as equivalent length of side or equivalent diameter size, and second is classified as
Number corresponding to equivalent length of side or equivalent diameter, third are classified as connected probability, and the 4th is classified as the row type；Wherein equivalent length of side
Or equivalent diameter size：λA1,1>λB2,1>λA1,2.λ represents equivalent length of side or equivalent diameter.λA1,1Represent type A1 in block A
By the equivalent length of side for sorting from large to small rear 1st large sample in block；λB2,1Representing type B2 in hole B, (B2 is not by any in A
Block includes to have neither part nor lot in hole that is poly- and operating, such as inorganic hole) hole in it is macroporous by the after sorting from large to small the 1st
Equivalent diameter；λA1,2By the equivalent length of side for sorting from large to small rear 2nd large sample in the block of type A1 in table block body A.
The present invention extracts main constituents, then based on mixing fractal theory, needle from true rock sample data
Each component is fitted using fractal parameter respectively, is finally combined each component, is obtained based on mixing fractal theory
Blowhole model.Further combined with pore communication algorithm, the connected probability matrix between variety classes hole is obtained, most
A kind of unconventional reservoir permeability evaluation method is constructed eventually, is realized and is rapidly and accurately assessed unconventional reservoir apparent permeability.
2 hole equivalent diameters meet first row from big to small.Using capillary force than size when diameter is equal, non-wetted phase is driven
It is arranged from big to small for capillary force when wetting phase, capillary force is arranged from small to large when wetting phase displacement non-wetted phase.
In step 400, gathered and operated before sequence.
Poly- and operating process specifically is as follows：If in image including organic mass body A1, in organic matter block A1
Hole is organic hole B1, then each organic matter block A1 is taken as porous media, it is organic in each organic matter block A1
Hole B1 meets its fractal cloth, then organic matter block A1 is only taken to participate in sequence in this step, it will not be in organic matter block
Organic hole B1 participate in sequence；If certain in is a kind of, such as inorganic hole B2 is not affiliated with any kind in component characteristics, that
, inorganic hole B2 participation sequence.
In an embodiment of the invention, the step of calculating permeability is as follows：
Wherein, the number of iterations i divides shape scale F, the hole or block side length λ that i-th iteration generatesi, maximum pore or block
Body side length λmax, the hole or block number N that generate in i-th iterationi, the hole or block of iteration are not involved in i-th iteration
Number Ni-soild, A is the sectional area in rock core characterization unit face, CiFor the i-th row block body or the pore shape factor, ki(pavg) it is flat
Equal pressure pavgThe permeability of lower i-th row block body or hole.
λ1,λ2,…,λi,…,λmFor the block or hole equivalent diameter size of the first row the first row in matrix, N1,
N2,…,Ni,…,NmFor the secondary series block or number of apertures in matrix, such as the first row λ1Corresponding block or number of apertures
N1, the i-th row λiCorresponding block or number of apertures Ni。
ki(pavg) calculation formula needs judged according to the type of row each in matrix；
1. setting, component characteristics A, the sequence number consecutively of sub- component characteristics are respectively with A1, A2 ... An, pore character
For B, the sequence number consecutively of sub- pore character is with B1, B2 ... Bn, supplemental characteristic C, sub- supplemental characteristic with C1,
The sequence number consecutively of C2 ... Cn；In component characteristics A, each block for forming sub- component characteristics is porous media；When
In connected probability matrix when the i-th row the 4th column record behavior component characteristics A, ki(pavg) expression formula is：
Wherein, knD(j,pavg) it is average pressure pavgUnder, the permeability value of the sub- pore character of j-th stage；Calculation formula is as follows：
As j=1
As j=2
As j=3
And so on, as j=n
Wherein, λ1,λ2,…,λj,…,λnFor being comprised in porous media representated by the i-th row after being gathered and being operated
Hole equivalent diameter；N1,N2,…,Nj,…,NnTo be comprised in porous Jie representated by the i-th row after being gathered and being operated
The number of hole in matter；Wherein λ1≥λ2≥…≥λj≥…≥λn；
2. when in connected probability matrix the i-th row the 4th column record behavior it is aforementioned poly- and hole when, permeability item ki(pavg)
For
3. as the i-th row the 4th column record behavior sub- supplemental characteristic C1 in connected probability matrix, permeability item ki(pavg)
For
Formula 4. (1) mesostroma permeability item kmatrixFor
It is above it is various in, μ is gas viscosity, and M is number of moles of gas, and R is ideal gas constant, and T is temperature, ρavgIt is flat
Equal density, pavgFor average pressure, α is tangential momentum adjustment factor,For n kind matrix minerals percentage, k1,
k2,…,knFor the gas permeability of n kind matrix minerals.
Embodiment one：
1, Sichuan Basin stone column block Lower Silurian Series Longma small stream group shale samples are taken, are using SEM as rock imaging technique
Example, is scanned sample.
2, scanning result is analyzed, mixing fractal model section is taken to be square, side length is 380 μm.Pyrite characterizes substantially
Unit is single pyrite, the square that organic matter basic unit face is 2.08 μm of side length.
3, the fractal parameter after each fitting is combined, obtains the blowhole model based on mixing fractal theory, led to
It crosses after the connected probability between variety classes hole in blowhole model is arranged in a matrix in hole connectivity algorithm and counts
It calculates, the assessment result of the apparent permeability of current unconventional reservoir can be obtained.
Wherein, organic matter accounts for the 3% of mixing fractal model area, and pyrite has 82 in mixing fractal model plane,
Pyrite average diameter 182nm, microcrack width is between 20nm-85nm, mean breadth 40nm, 5.18 μm of average length, puts down
Equal fracture spacing is 50.3/10000 μm2。
4, it is 99.17nD that the rock sample permeability, which is calculated, in this method, and it is 91.36nD, error that experiment, which measures permeability,
8.55%, computing permeability result is close with true core computing permeability result, and this method is reliable.
So far, although those skilled in the art will appreciate that present invention has been shown and described in detail herein multiple shows
Example property embodiment still without departing from the spirit and scope of the present invention, still can according to the present disclosure directly
Determine or deduce out many other variations or modifications consistent with the principles of the invention.Therefore, the scope of the present invention is understood that and recognizes
It is set to and covers all such other variations or modifications.
Claims (10)
1. a kind of unconventional reservoir permeability evaluation method based on core image, which is characterized in that include the following steps：
Step 100, respective image is obtained after true rock sample being scanned and identified using rock imaging technique, by rock in image
The organic matter block and ore particles Type division of sample are component characteristics, by organic matter block, in ore particles and Rock Matrix
Hole information be divided into pore character, the other feature in rock sample is divided into supplemental characteristic；
Step 200, the frequency spectrum of every kind of feature is extracted respectively, is then fitted, is obtained with mixing fractal method to analysis result
The fractal parameter of each individual features；
Step 300, the fractal parameter after each fitting is combined, obtains the blowhole model based on mixing fractal theory；
Step 400, the connected probability between variety classes hole in blowhole model is pressed by hole connectivity algorithm and is connected
It is calculated after logical probability matrix arrangement, the assessment result of the apparent permeability of current unconventional reservoir can be obtained.
2. unconventional reservoir permeability evaluation method according to claim 1, which is characterized in that
The pore character in the step 100 refers to organic matter block, in organic matter hole and ore in ore particles
Intergranular pore or intragranular hole；The supplemental characteristic refers to intrinsic fracture or man-made fracture feature in the rock sample.
3. unconventional reservoir permeability evaluation method according to claim 1, which is characterized in that
The rock imaging technique includes SEM, FIB-SEM, Micro-CT, nano-CT.
4. unconventional reservoir permeability evaluation method according to claim 3, which is characterized in that
The component characteristics are different according to the type of ore particles and are subdivided into multiple corresponding sub- component characteristics；The hole is special
Sign corresponds to the difference of the sub- component characteristics and is subdivided into multiple corresponding sub- pore characters；The supplemental characteristic according to difference
The relationship of the sub- component characteristics and be subdivided into multiple corresponding sub- supplemental characteristics.
5. unconventional reservoir permeability evaluation method according to claim 4, which is characterized in that
The frequency spectrum of the component characteristics is built upon the equivalent size distribution frequency of respective subgroup point in wherein each sub- component characteristics
In spectrum, the frequency spectrum of the pore character is built upon on the spread spectrum of the effective aperture of each sub- pore character, and the auxiliary is special
The frequency spectrum of sign is built upon in fracture spacing represented by sub- supplemental characteristic, seam length and slit width.
6. unconventional reservoir permeability evaluation method according to claim 5, which is characterized in that each subgroup dtex
The size distribution frequency spectrum that the equivalent size for every grade of sub- component characteristics that the fractal parameter of sign obtains and its correspondence number collectively form,
It is equal to the size distribution frequency spectrum of the sub- component characteristics；Every grade of hole that the fractal parameter of each sub- pore character obtains
The size distribution frequency spectrum that equivalent diameter and its correspondence number collectively form, is equal to the pore size point of the sub- pore character
Cloth frequency spectrum；The equivalent size for every grade of sub- supplemental characteristic that the fractal parameter of each sub- supplemental characteristic obtains and its correspondence number
The crack frequency spectrum collectively formed, with the crack frequency spectrum for being equal to sub- supplemental characteristic.
7. unconventional reservoir permeability evaluation method according to claim 6, which is characterized in that
The mode for establishing connected probability matrix using interconnected pore probabilistic algorithm is as follows：By each subcharacter in three kinds of features according to
The mode of equivalent length of side or equivalent diameter sorts from large to small, wherein and first is classified as equivalent length of side or equivalent diameter size, and second
It is classified as number corresponding to equivalent length of side or equivalent diameter, third column connected probability, the 4th is classified as the row component type.
8. unconventional reservoir permeability evaluation method according to claim 7, which is characterized in that
Gathered before sequence and is operated.
9. unconventional reservoir permeability evaluation method according to claim 8, which is characterized in that
The step of calculating permeability is as follows：
Wherein, the number of iterations i divides shape scale F, hole/block side length λ that i-th iteration generatesi, maximum pore/block side length
λmax, hole/block number N for generating in i-th iterationi, hole/block number of iteration is not involved in i-th iteration
Ni-soild, A is the sectional area in rock core characterization unit face, CiFor the i-th row block body/pore shape factor, ki(pavg) it is in mean pressure
Power pavgThe permeability of lower i-th row block body/hole.
10. unconventional reservoir permeability evaluation method according to claim 9, which is characterized in that
ki(pavg) calculation formula needs judged according to the type of row each in matrix；
1. setting, component characteristics A, the sequence number consecutively of sub- component characteristics are respectively with A1, A2 ... An, pore character B,
The sequence number consecutively of its sub- pore character is with B1, B2 ... Bn, supplemental characteristic C, sub- supplemental characteristic is with C1, C2 ... Cn
Sequence number consecutively；In component characteristics A, each block for forming sub- component characteristics is porous media；Work as connected probability
In matrix when the i-th row the 4th column record behavior component characteristics A, ki(pavg) expression formula is：
Wherein, knD(j,pavg) it is average pressure pavgUnder, the permeability value of the sub- pore character of j-th stage；Calculation formula is as follows：
As j=1
As j=2
As j=3
And so on, as j=n
Wherein, λ1,λ2,…,λj,…,λnFor the hole being comprised in porous media representated by the i-th row after being gathered and being operated
Equivalent diameter；N1,N2,…,Nj,…,NnFor being comprised in porous media representated by the i-th row after being gathered and being operated
The number of hole；Wherein λ1≥λ2≥…≥λj≥…≥λn；
2. when the i-th row the 4th column record behavior is not by poly- hole simultaneously in connected probability matrix, permeability item ki(pavg) be
3. as the i-th row the 4th column record behavior sub- supplemental characteristic C1 in connected probability matrix, permeability item ki(pavg) be
4. formula (1) the mesostroma permeability item kmatrixFor
It is above it is various in, μ is gas viscosity, and M is number of moles of gas, and R is ideal gas constant, and T is temperature, ρavgIt is averagely close
Degree, pavgFor average pressure, α is tangential momentum adjustment factor,For n kind matrix minerals percentage, k1,k2,…,
knFor the apparent permeability of n kind matrix minerals.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810546608.2A CN108829950B (en) | 2018-05-31 | 2018-05-31 | Unconventional reservoir permeability evaluation method based on core image |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810546608.2A CN108829950B (en) | 2018-05-31 | 2018-05-31 | Unconventional reservoir permeability evaluation method based on core image |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108829950A true CN108829950A (en) | 2018-11-16 |
| CN108829950B CN108829950B (en) | 2021-04-27 |
Family
ID=64145245
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810546608.2A Active CN108829950B (en) | 2018-05-31 | 2018-05-31 | Unconventional reservoir permeability evaluation method based on core image |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108829950B (en) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108956415A (en) * | 2018-05-31 | 2018-12-07 | 中国科学院力学研究所 | A kind of permeability saturation curve calculation method for unconventional reservoir sample |
| CN109632604A (en) * | 2019-01-04 | 2019-04-16 | 中国海洋石油集团有限公司 | A kind of pore scale is to rock core scale Relative Permeabilty From Polymer Displacement Experiment method of roughening |
| CN109933939A (en) * | 2019-03-22 | 2019-06-25 | 西南石油大学 | The method for numerical simulation of the unconventional crack initiation of dual media reservoir multiple cracking and extension |
| CN110567858A (en) * | 2019-10-17 | 2019-12-13 | 西南石油大学 | method for predicting shale nano-pore permeability based on fractal theory |
| CN112302606A (en) * | 2020-07-07 | 2021-02-02 | 西南石油大学 | Inversion interpretation method for output profile of low-permeability gas reservoir fractured horizontal well |
| CN112967147A (en) * | 2021-02-04 | 2021-06-15 | 中海石油(中国)有限公司海南分公司 | Bedrock yield contribution rate calculation method considering multi-scale cracks |
| CN113189122A (en) * | 2021-05-13 | 2021-07-30 | 中海石油(中国)有限公司海南分公司 | Perforation damage indoor comprehensive evaluation method |
| CN113297779A (en) * | 2021-06-23 | 2021-08-24 | 中国石油大学(华东) | Shale permeability interpretation method based on dual-medium pore network model |
| CN116342541A (en) * | 2023-03-29 | 2023-06-27 | 中国矿业大学 | Rock-soil body permeability calculation method based on adjacent image pore fusion reconstruction |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100089124A1 (en) * | 2008-09-26 | 2010-04-15 | North Dakota State University | Integrated porous rigid wall and flexible wall permeability test device for soils |
| CN104018829A (en) * | 2014-05-23 | 2014-09-03 | 中国地质大学（北京） | Method for measuring gas and water relative permeability curve through coal-bed gas well production data |
| US8881587B2 (en) * | 2011-01-27 | 2014-11-11 | Schlumberger Technology Corporation | Gas sorption analysis of unconventional rock samples |
| CN105279790A (en) * | 2014-06-13 | 2016-01-27 | 中国石油化工股份有限公司 | Fracture network 3D digital core modeling method |
| US20160337444A1 (en) * | 2015-05-13 | 2016-11-17 | Lunatech, Llc | Social network for electronic vapor device users |
| CN107449706A (en) * | 2017-06-06 | 2017-12-08 | 湖北工业大学 | Deformation soil body saturation, Unsaturated Hydraulic Conductivity Forecasting Methodology based on fractal theory |
-
2018
- 2018-05-31 CN CN201810546608.2A patent/CN108829950B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100089124A1 (en) * | 2008-09-26 | 2010-04-15 | North Dakota State University | Integrated porous rigid wall and flexible wall permeability test device for soils |
| US8881587B2 (en) * | 2011-01-27 | 2014-11-11 | Schlumberger Technology Corporation | Gas sorption analysis of unconventional rock samples |
| CN104018829A (en) * | 2014-05-23 | 2014-09-03 | 中国地质大学（北京） | Method for measuring gas and water relative permeability curve through coal-bed gas well production data |
| CN105279790A (en) * | 2014-06-13 | 2016-01-27 | 中国石油化工股份有限公司 | Fracture network 3D digital core modeling method |
| US20160337444A1 (en) * | 2015-05-13 | 2016-11-17 | Lunatech, Llc | Social network for electronic vapor device users |
| CN107449706A (en) * | 2017-06-06 | 2017-12-08 | 湖北工业大学 | Deformation soil body saturation, Unsaturated Hydraulic Conductivity Forecasting Methodology based on fractal theory |
Non-Patent Citations (3)
| Title |
|---|
| GAOHUI CAO: "《A 3D coupled model of organic matter and inorganic matrix for calculating the permeability of shale》", 《ELSEVIER》 * |
| 李曹雄: "《孔隙网络对页岩储层自吸特征的影响与应用》", 《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》 * |
| 李玉丹 等: "《基于分形理论的页岩基质表现渗透率研究》", 《油气地质和采收率》 * |
Cited By (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108956415B (en) * | 2018-05-31 | 2020-03-17 | 中国科学院力学研究所 | Method for calculating relative permeability curve of unconventional reservoir sample |
| CN108956415A (en) * | 2018-05-31 | 2018-12-07 | 中国科学院力学研究所 | A kind of permeability saturation curve calculation method for unconventional reservoir sample |
| CN109632604A (en) * | 2019-01-04 | 2019-04-16 | 中国海洋石油集团有限公司 | A kind of pore scale is to rock core scale Relative Permeabilty From Polymer Displacement Experiment method of roughening |
| CN109632604B (en) * | 2019-01-04 | 2021-06-15 | 中国海洋石油集团有限公司 | Method for coarsening relative permeability of polymer flooding from pore size to core size |
| CN109933939A (en) * | 2019-03-22 | 2019-06-25 | 西南石油大学 | The method for numerical simulation of the unconventional crack initiation of dual media reservoir multiple cracking and extension |
| CN110567858A (en) * | 2019-10-17 | 2019-12-13 | 西南石油大学 | method for predicting shale nano-pore permeability based on fractal theory |
| CN110567858B (en) * | 2019-10-17 | 2020-11-03 | 西南石油大学 | Method for predicting shale nano-pore permeability based on fractal theory |
| CN112302606B (en) * | 2020-07-07 | 2021-08-24 | 西南石油大学 | Inversion interpretation method for output profile of low-permeability gas reservoir fractured horizontal well |
| CN112302606A (en) * | 2020-07-07 | 2021-02-02 | 西南石油大学 | Inversion interpretation method for output profile of low-permeability gas reservoir fractured horizontal well |
| CN112967147A (en) * | 2021-02-04 | 2021-06-15 | 中海石油(中国)有限公司海南分公司 | Bedrock yield contribution rate calculation method considering multi-scale cracks |
| CN112967147B (en) * | 2021-02-04 | 2024-04-02 | 中海石油(中国)有限公司海南分公司 | Bedrock yield contribution rate calculation method considering multi-scale cracks |
| CN113189122A (en) * | 2021-05-13 | 2021-07-30 | 中海石油(中国)有限公司海南分公司 | Perforation damage indoor comprehensive evaluation method |
| CN113189122B (en) * | 2021-05-13 | 2023-12-22 | 中海石油(中国)有限公司海南分公司 | Perforation damage indoor comprehensive evaluation method |
| CN113297779A (en) * | 2021-06-23 | 2021-08-24 | 中国石油大学(华东) | Shale permeability interpretation method based on dual-medium pore network model |
| CN113297779B (en) * | 2021-06-23 | 2022-07-29 | 中国石油大学(华东) | Shale permeability interpretation method based on dual-medium pore network model |
| CN116342541A (en) * | 2023-03-29 | 2023-06-27 | 中国矿业大学 | Rock-soil body permeability calculation method based on adjacent image pore fusion reconstruction |
| CN116342541B (en) * | 2023-03-29 | 2024-03-22 | 中国矿业大学 | Rock-soil body permeability calculation method based on adjacent image pore fusion reconstruction |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108829950B (en) | 2021-04-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108829950A (en) | A kind of unconventional reservoir permeability evaluation method based on core image | |
| Li et al. | Multi-scale quantitative characterization of 3-D pore-fracture networks in bituminous and anthracite coals using FIB-SEM tomography and X-ray μ-CT | |
| Tang et al. | The effect of the variation in material composition on the heterogeneous pore structure of high-maturity shale of the Silurian Longmaxi formation in the southeastern Sichuan Basin, China | |
| CN103278436B (en) | Quantitative characterization method of low penetration double-medium sandstone oil reservoir microscopic aperture structure | |
| Guo et al. | Quantitative pore characterization and the relationship between pore distributions and organic matter in shale based on Nano-CT image analysis: a case study for a lacustrine shale reservoir in the Triassic Chang 7 member, Ordos Basin, China | |
| CN108956415B (en) | Method for calculating relative permeability curve of unconventional reservoir sample | |
| CN104374682A (en) | Rock core CT (Computed Tomography) scanning analysis method and device | |
| Cao et al. | Classification and controlling factors of organic pores in continental shale gas reservoirs based on laboratory experimental results | |
| CN105551004A (en) | Core CT image processing-based remaining oil micro-occurrence representing method | |
| CN109187305A (en) | A kind of shale structure characterization methods and device | |
| Xie et al. | Quantitative characterization of pore size and structural features in ultra-low permeability reservoirs based on X-ray computed tomography | |
| Cao et al. | Characterization of pore structures and gas transport characteristics of Longmaxi shale | |
| CN109632594B (en) | Method for characterizing pore throat characteristics of tight reservoir on basis of high-pressure mercury intrusion multi-scale | |
| CN110487693A (en) | A kind of method of determining mud shale different type porosity | |
| CN107132171A (en) | A kind of method for determining compact reservoir pore-size distribution based on pressure mercury N2 adsorption translocation data | |
| CN111521539B (en) | Quantitative characterization method for different types of pore structures in compact conglomerate reservoir | |
| CN104833789A (en) | Method for determining air/water relation by using compact sandstone microscopic aperture structure | |
| CN110132816A (en) | Method for analyzing pore structure of organic matter in shale of ancient world | |
| CN111827996B (en) | Multi-parameter comprehensive qualitative compact sandstone reservoir classification method based on mechanical properties | |
| CN109632604B (en) | Method for coarsening relative permeability of polymer flooding from pore size to core size | |
| Cao et al. | Thermally enhanced shale gas recovery: microstructure characteristics of combusted shale | |
| Ghassemzadeh et al. | Pore network simulation of fluid imbibition into paper during coating: II. Characterization of paper's morphology and computation of its effective permeability tensor | |
| Jiang et al. | Investigating the pore structure characteristics and reservoir capacities of lower Jurassic continental shale reservoirs in the northeastern Sichuan Basin, China | |
| Wang et al. | Comparative study of nanoscale pore structure of L ower P alaeozoic marine shales in the M iddle‐U pper Y angtze area, China: I mplications for gas production potential | |
| Li et al. | Detailed characterization of micronano pore structure of tight sandstone reservoir space in three dimensional space: a case study of the Gao 3 and Gao 4 members of Gaotaizi reservoir in the Qijia area of the Songliao basin |
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 |