CN107132171B - A method of compact reservoir pore-size distribution is determined based on pressure mercury-N2 adsorption translocation data - Google Patents
A method of compact reservoir pore-size distribution is determined based on pressure mercury-N2 adsorption translocation data Download PDFInfo
- Publication number
- CN107132171B CN107132171B CN201710238534.1A CN201710238534A CN107132171B CN 107132171 B CN107132171 B CN 107132171B CN 201710238534 A CN201710238534 A CN 201710238534A CN 107132171 B CN107132171 B CN 107132171B
- Authority
- CN
- China
- Prior art keywords
- pore
- aperture
- frequency
- adsorption
- formula
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume, or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume, or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
- G01N15/0886—Mercury porosimetry
Abstract
The invention discloses a kind of methods for determining compact reservoir pore-size distribution based on pressure mercury-N2 adsorption translocation data, obtain the pore-size distribution situation in the reservoir first with N2 adsorption method, N2 adsorption institute gaging hole larynx is distributed as the distribution situation of mesoporous or micropore；Secondly the pore-size distribution situation in the reservoir is obtained using pressure mercury method, mercury injection method institute gaging hole larynx is distributed as mesoporous or macropore；The reallocation of throat distribution situation is carried out using the pore-size distribution result of the two measurement again, seeks each section aperture proportion, and then obtain the pore-size distribution situation of gamut.The present invention solves the problems, such as different stage throat distribution in fine and close oil and gas reservoir, realizes the global analysis of pore throat character distribution characteristics in reservoir, to improve the globality and accuracy of compact reservoir micropore structure characterization.
Description
Technical field
The invention belongs to petroleum natural gas exploration fields, are related to the research method of reservoir full aperture distribution, specifically
A method of compact reservoir pore-size distribution is determined based on pressure mercury-N2 adsorption translocation data.
Background technique
With the continuous development of oil-gas exploration, it is more next that the observation of the microcosmic reservoir pore throat character feature of compact reservoir describes method
It is abundanter, it is combined under the microscope from traditional Mercury injection, arrives CT three-dimensional pore space analogue technique, the research of reservoir micro throat is continuous
It is combined with mathematics and computer, and tends to quantification.
In the research method of existing reservoir pores distribution, since different measurement method principles are different and experiment condition
The pore-size distribution interval range of limitation, the measurement of every kind of experimental method can all be not quite similar.Mercury injection method is under the conditions of different pressures
It is pressed into the intrusion in rock sample hole by measurement mercury, to calculate the distribution situation in aperture, N2 adsorption method is to measure not
With liquid nitrogen under pressure in the adsorbance of sample interior pore throat, pore-size distribution is calculated according to different aperture model.Liquid nitrogen is thinner
Adsorption capacity is stronger in small hole, and the pore radius range measured is generally mesoporous or micropore, many of mercury injection method
Tiny hole is that no matter great pressure can not also be such that mercury solution is pressed into, thus can only measure the aperture point of mesoporous and macropore
(in the world, generally aperture is divided into three classes by size: aperture≤2nm is micropore to cloth, and aperture is in 2-50nm range is
Hole, aperture >=50nm are macropore).And rock sample microscopic aperture is distributed pore size distribution or mesoporous micropore point in the macropore for being unlikely to be single
Cloth, two kinds of measurement methods are unable to be visually observed that the pore-size distribution situation in gamut, and individually a kind of method measures meeting
Cause measurement range narrow, is unfavorable for comprehensively researching and analysing reservoir micropore structure, especially in compact reservoir
In research process, micropore macropore no less important should more increase the research to pore throat character holistic distribution characteristics.
Summary of the invention
The purpose of the present invention is overcome above-mentioned reservoir pores narrow distribution range existing in the prior art, can not intuitively observe
The problem of pore-size distribution situation, provides a kind of method for determining compact reservoir pore-size distribution based on pressure mercury-N2 adsorption translocation data,
The pore size distribution range for pressing mercury and N2 adsorption method translocation to obtain is unified, it can be realized and press mercury and N2 adsorption money in known reservoir
It is final to determine compact reservoir gamut pore-size distribution feature by the pore-size distribution feature of both analyses in the case of material, there is behaviour
Control is simple, the visual feature high with accuracy by force.
The technical scheme is that a kind of determine compact reservoir pore-size distribution based on pressure mercury-N2 adsorption translocation data
Method includes the following steps:
Step (1) determines the N2 adsorption parameter for the pore diameter range that nitrogen adsorption method can detect, uses n1Indicate that N2 adsorption method is total
The aperture number measured, RiIndicate the size in i-th of aperture, αiIndicate frequency shared by i-th of aperture, XiIndicate i-th of aperture
Shared volume, V1Indicate the total pore volume that N2 adsorption method measures, i=1,2 ... n1。
Step (2) determines the pressure mercury parameter for the pore diameter range that mercury injection method can detect, uses n2Indicate that mercury method measures aperture altogether
Number, rjIndicate the size in j-th of aperture, βjIndicate frequency shared by j-th of aperture, YjIndicate volume shared by j-th of aperture,
V2Indicate the total pore volume that pressure mercury method measures, j=1,2 ..., n2。
Step (3), due to nitrogen adsorption method, two methods of mercury injection method measure range some be overlapping, and be overlapped
The volume of part pore throat be it is certain, then have Xi=Yj, it may be assumed that
V1×αi=V2×βj... ... ... ... ... formula is 1.
Therefore the variation tendency of the throat distribution for the lap that two methods measure is identical.By nitrogen adsorption method and
Two groups of throat distribution parameters that mercury injection method measures directly read the number in overlapping aperture, indicate that nitrogen adsorption method and mercury injection method are surveyed with x
The number in the overlapping aperture obtained.
Step (4) determines the frequency α in integration step (1)iWith the frequency β in step (2)jRequired mapping ratio.If λ
Mapping ratio needed for indicating integration process, then also as mercury injection method measures large pore to λ and nitrogen adsorption method measures small aperture
Overall prior probability ratio, 2. the size of λ determines by formula:
2. 1. simultaneous formula is obtained with formula:
λtIt indicates to calculate obtained overall prior probability ratio by t-th of lap；Wherein i=1,2 ... n1；J=1,
2,…,n2；T=1,2 ..., x；
X prior probability ratio is 3. found out altogether by formula, by the obtained x prior probability of this x overlapping interval than logical
It crosses formula and 4. determines average prior probability ratio:
Step (5) determines the frequency α in integration step (1)iWith the frequency β in step (2)jFrequency afterwards.Indicate whole
The frequency in each aperture after conjunction, the then frequency in each aperture after integratingValue pass through following process determine:
I, when pore size in the range of micropores of nitrogen absorption measurement and when (aperture not comprising lap), frequencyBy
5. formula determines:
II, when mesopore range of the pore size in the lap that two methods measure, frequencyIt is 6. determined by formula:
(wherein, work as i=n1When-x+1, j=1 ... ... works as i=n1When, j=x)
III, when pore size in the macropore range that mercury injection method measures (aperture not comprising lap), frequencyBy
7. formula determines:
Micropore size≤the 2nm, the pore diameter range of mesoporous are 2-50nm, macropore diameter >=50nm,
Step (6), what is obtained by step (5)SummationUnitization transformation is carried out, that is to say so that each range
The sum of pore-size distribution frequency become the ratio that 100, Ψ indicates unitization, 8. the size of Ψ is determined by formula:
Using formula 8. in Ψ size by formula 9. determination finally obtain the difference measured comprising two kinds of experimental methods
The whole throat distribution frequency of pore diameter range, αkIndicate final throat distribution frequency,
Pore-size distribution histogram is made using finally obtained data, can be visually observed that the aperture in gamut
Distribution situation is conducive to the research to pore throat character holistic distribution characteristics.
Preferably, in step (1), according to " rock specific surface and pore-size distribution measure static N2 adsorption volumetric method " SY/T
Process as defined in 6154-1995 standard carries out N2 adsorption experiment to rock core, determines the hole for the pore diameter range that nitrogen adsorption method can detect
Diameter distribution parameter.
Preferably, in step (2), according to " measurement of rock capillary pressure curve " SY/T 5346-2005 normal process to rock
The heart carries out pressure mercury experiment, determines the pressure mercury parameter for the pore diameter range that mercury injection method can detect.
Beneficial effects of the present invention: in the embodiment of the present invention, a kind of determined based on pressure mercury-N2 adsorption translocation data is provided and is caused
The method of close reservoir pores distribution, the small-bore distribution measured by N2 adsorption experiment and intraventricular pressure mercury test the large aperture measured
Two kinds of data are distributed, by the data of two kinds of pore size distribution ranges of different sizes, pass through two groups of data statistics, two groups of data overlaps
Pore diameter range, due to a certain overlapping interval volume that two methods measure be it is identical, can pass through calculate obtain pressure mercury
The prior probability ratio of macrovoid total volume that method measurement obtains and the fine pore totality that nitrogen adsorption method measurement obtains, wherein priori
Likelihood ratio is acquired respectively using the related data in overlapping each in overlapping range aperture, then takes its average value again, most
Required prior probability ratio is utilized afterwards, and the pore-size distribution frequency that two kinds of measurement methods are obtained is integrated, and the aperture after integration point
The progress of cloth frequency is unitization, finally obtains the graph of pore diameter distribution of each range.Therefore, the present invention can solve fine and close oil and gas reservoir
The pore-size distribution problem of middle different stage, realizes the global analysis of pore throat character feature in fine and close oil and gas reservoir, to improve
The globality and accuracy of compact reservoir micropore structure characterization.
Detailed description of the invention
Fig. 1 is flow chart of the method for the present invention；
Fig. 2 is nitrogen adsorption method graph of pore diameter distribution provided by the invention；
Fig. 3 is mercury injection method graph of pore diameter distribution provided by the invention；
Fig. 4 is pressure mercury provided by the invention and nitrogen adsorption method graph of pore diameter distribution；
Fig. 5 is finally obtained pressure mercury provided by the invention-N2 adsorption translocation gamut pore-size distribution.
Specific embodiment
With reference to the accompanying drawing, the specific embodiment of the present invention is described in detail, it is to be understood that of the invention
Protection scope be not limited by the specific implementation.
Referring to Fig. 1, compact reservoir aperture is determined based on pressure mercury-N2 adsorption translocation data the embodiment of the invention provides a kind of
The method of distribution, includes the following steps:
Step (1) determines the N2 adsorption parameter for the pore diameter range that nitrogen adsorption method can detect, uses n1Indicate that N2 adsorption method is total
The aperture number measured, RiIndicate the size in i-th of aperture, αiIndicate frequency shared by i-th of aperture, XiIndicate i-th of aperture
Shared volume, V1Indicate the total pore volume that N2 adsorption method measures, i=1,2 ... n1。
Step (2) determines the pressure mercury parameter for the pore diameter range that mercury injection method can detect, uses n2Indicate that mercury method measures aperture altogether
Number, rjIndicate the size in j-th of aperture, βjIndicate frequency shared by j-th of aperture, YjIndicate volume shared by j-th of aperture,
V2Indicate the total pore volume that pressure mercury method measures, j=1,2 ..., n2。
Step (3), due to nitrogen adsorption method, two methods of mercury injection method measure range some be overlapping, and be overlapped
The volume of part pore throat be it is certain, then have Xi=Yj, it may be assumed that
V1×αi=V2×βj... ... ... ... ... formula is 1.
Therefore the variation tendency of the throat distribution for the lap that two methods measure is identical.By nitrogen adsorption method and
Two groups of throat distribution parameters that mercury injection method measures directly read the number in overlapping aperture, indicate that nitrogen adsorption method and mercury injection method are surveyed with x
The number in the overlapping aperture obtained.
Step (4) determines the frequency α in integration step (1)iWith the frequency β in step (2)jRequired mapping ratio.If λ
Mapping ratio needed for indicating integration process, then also as mercury injection method measures large pore to λ and nitrogen adsorption method measures small aperture
Overall prior probability ratio, 2. the size of λ determines by formula:
2. 1. simultaneous formula is obtained with formula:
λtIt indicates to calculate obtained overall prior probability ratio by t-th of lap；Wherein i=1,2 ... n1；J=1,
2,…,n2；T=1,2 ..., x；
X prior probability ratio is 3. found out altogether by formula, by the obtained x prior probability of this x overlapping interval than logical
It crosses formula and 4. determines average prior probability ratio:
Step (5) determines the frequency α in integration step (1)iWith the frequency β in step (2)jFrequency afterwards.Indicate whole
The frequency in each aperture after conjunction, the then frequency in each aperture after integratingValue pass through following process determine:
I, when pore size in the range of micropores of nitrogen absorption measurement and when (aperture not comprising lap), frequencyBy
5. formula determines:
II, when mesopore range of the pore size in the lap that two methods measure, frequencyIt is 6. determined by formula:
(wherein, work as i=n1When-x+1, j=1 ... ... works as i=n1When, j=x)
III, when pore size in the macropore range that mercury injection method measures (aperture not comprising lap), frequencyBy
7. formula determines:
Micropore size≤the 2nm, the pore diameter range of mesoporous are 2-50nm, macropore diameter >=50nm,
Step (6), what is obtained by step (5)SummationUnitization transformation is carried out, that is to say so that each range
The sum of pore-size distribution frequency become the ratio that 100, Ψ indicates unitization, 8. the size of Ψ is determined by formula:
Using formula 8. in Ψ size by formula 9. determination finally obtain the difference measured comprising two kinds of experimental methods
The whole throat distribution frequency of pore diameter range, αkIndicate final throat distribution frequency,
Pore-size distribution histogram is made using finally obtained data, can be visually observed that the aperture in gamut
Distribution situation is conducive to the research to pore throat character holistic distribution characteristics.
Further, in step (1), according to " rock specific surface and pore-size distribution measure static N2 adsorption volumetric method " SY/T
Process as defined in 6154-1995 standard carries out N2 adsorption experiment to rock core, determines the hole for the pore diameter range that nitrogen adsorption method can detect
Diameter distribution parameter.
Further, in step (2), according to " measurement of rock capillary pressure curve " SY/T 5346-2005 normal process pair
Rock core carries out pressure mercury experiment, determines the pressure mercury parameter for the pore diameter range that mercury injection method can detect.
Definitely to illustrate actual effect of the present invention, now chooses Song-liao basin oil field tight sand sample and (bury
Deep 2342.5m) research object as this embodiment.
Step 1, the measurement of N2 adsorption experimental method pore diameter range: choosing representative core sample, according to " rock
Specific surface and pore-size distribution measure static state N2 adsorption volumetric method SY/T 6154-1995 " process as defined in standard is to rock core progress nitrogen
Adsorption experiment obtains the N2 adsorption parameter of mesoporous and micropore, n1Indicate the aperture number that N2 adsorption method measures altogether, n1=19,
RiIndicate the size in i-th of aperture, αiIndicate frequency shared by i-th of aperture, XiIndicate volume shared by i-th of aperture, V1Table
Show the total pore volume that N2 adsorption method measures, i=1,2 ... 19, specific experiment the results are shown in Table 1, and pore-size distribution is made by table 1
Scheme (Fig. 2).
1 nitrogen adsorption method pore size distribution data of table
| Serial number | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 |
| N2 adsorption aperture Ri/nm | 0.015 | 0.025 | 0.04 | 0.065 | 0.1 | 0.16 | 0.25 | 0.4 | 0.63 | 1 |
| Frequency αi/ % | 0.002 | 0.003 | 0.006 | 0.008 | 0.014 | 0.02 | 0.034 | 0.052 | 0.874 | 3.067 |
| Serial number | 11 | 12 | 13 | 14 | 15 | 16 | 17 | 18 | 19 | 20 |
| N2 adsorption aperture Ri/nm | 1.6 | 2.5 | 4 | 6.3 | 10 | 16 | 25 | 40 | 63 | |
| Frequency αi/ % | 3.004 | 3.771 | 5.684 | 8.322 | 11.976 | 9.913 | 16.516 | 25.333 | 11.408 |
Step 2, pressure mercury experimental method pore diameter range measurement: according to " measurement of rock capillary pressure curve " SY/T 5346-
2005 normal process carry out the pressure mercury parameter that the experiment of rock core mercury injection method obtains macropore and mesoporous, measure aperture number n altogether2=11,
The size r in the aperture measuredjAnd shared frequency βjIt is shown in Table 2, YjIndicate volume shared by j-th of aperture, j=1,2 ..., 11,
V2Indicate the total pore volume that pressure mercury method measures.Table 2 is organized into graph of pore diameter distribution (Fig. 3), and by two groups of numbers of Tables 1 and 2
According to being organized into a graph of pore diameter distribution (Fig. 4).
2 mercury injection method pore size distribution data of table
Step 3 has Fig. 4 to analyze it is found that the pore-size distribution trend in the lap section that two methods measure is identical, only
It is that proportion is different in respective measurement range, practical is pore throat space having the same.According in 1 data of table and table 2
Data determine that nitrogen adsorption method and the overlapping number for the pore diameter range that mercury injection method measures are 6, the pore size distribution data of lap
It is shown in Table 3.
3 lap pore size distribution data of table
Step 4 determines the frequency α in integration step (1)iWith the frequency β in step (2)jRequired mapping ratio.If λ
Mapping ratio needed for indicating integration process, then also as mercury injection method measures large pore to λ and nitrogen adsorption method measures small aperture
Overall prior probability ratio, the total pore volume V that the size of λ is measured by the pressure mercury method in step (2)2Size and step
(1) the total pore volume V that the N2 adsorption method in measures2Size determine, it may be assumed that
Again because when being located at a certain overlapping interval, this overlapping interval volume and mercury injection method measurement of determination of nitrogen adsorption
This overlapping interval volume be it is equal, then have:
Xi=Yj, i.e. V1×αi=V2×βj... ... ... ... formula is 2.
By formula, 1. 2. simultaneous is obtained with formula:3. formula is priori
The evaluation formula of probability ratio, λtIt indicates to calculate obtained total pore size volume prior probability ratio by t-th of lap, works as i
When=14, j=1, t=1, pore diameter range are located at the 1st overlapping interval, and so on, when i=19, j=6, t=6, aperture model
It encloses and is located at the 6th overlapping interval.I=14,15 ..., 19；J=1,2 ..., 6；T=1,2 ..., 6.
3. the data of table 3 are substituted into formula, find out 6 prior probabilities, respectively λ1=2.585, λ2=2.019, λ3=
1.488 λ4=1.980, λ5=2.862, λ6=1.083, formula is passed through by obtained 6 prior probabilities of this 6 overlapping intervals
4. calculating average prior probability ratio:
Obtained average prior probability ratio are as follows:
Step 5, the frequency and pressure mercury for integrating micropore and mesoporous that N2 adsorption experiment measures test the mesoporous and macropore measured
Frequency.The frequency in each aperture after indicating integration,Value pass through following process determine:
I, when aperture is in 0.015-6.3nm range, the frequency in each aperture of this rangeIt is 5. determined by formula:
II, when aperture is in 6.3-63nm range, the frequency in each aperture of this rangeIt is 6. determined by formula:
III, when aperture is in 63-630nm range, the frequency in each aperture of this rangeIt is 7. determined by formula:
Pore size distribution data after integration is as shown in table 4 below.
Table 4 integrates back aperture distributed data
Step 6 is obtained by above stepSummationUnitization transformation is carried out, that is to say so that each range
The sum of pore-size distribution frequency should be 100.Unitization ratio is Ψ, then hasThe different pore size range throat distribution frequency measured comprising two kinds of experimental methods is finally calculated
αk, then have(k=1,2,3 ..., 24, j=1,2,3 ..., 24), the different holes that available two methods measure at this time
The pore-size distribution frequency α of diameter rangek, as shown in table 5.
The final pore size distribution data of table 5
The pore-size distribution frequency in a wide range of interior different sections finally obtained, is organized into figure for this composite result to get arriving
The a wide range of pore-size distribution situation that two methods measure.It thus can no longer be bound to limit brought by experimental measurement method
System, obtains nanoscale to micron-sized reservoir pores distribution situation, such as Fig. 5.
In conclusion provided in an embodiment of the present invention determine compact reservoir aperture point based on pressure mercury-N2 adsorption translocation data
The method of cloth, the large aperture that the small-bore distribution and the experiment of intraventricular pressure mercury measured by N2 adsorption experiment measures are distributed two kinds of numbers
According to, by the data of two kinds of pore size distribution ranges of different sizes, by the pore diameter range of two groups of data statistics, two groups of data overlaps, by
It is identical in a certain overlapping interval volume that two methods measure, therefore can obtains mercury injection method measurement by calculating and obtain
The prior probability ratio for the fine pore totality that macrovoid total volume and nitrogen adsorption method measurement obtain, wherein prior probability ratio is to utilize
The related data in each overlapping aperture acquires respectively in overlapping range, then takes its average value again, finally utilizes required elder generation
Probability ratio is tested, the pore-size distribution frequency that two kinds of measurement methods are obtained is integrated, and the pore-size distribution frequency after integration carries out list
Positionization finally obtains the graph of pore diameter distribution of each range.Therefore, the present invention can solve different stage in fine and close oil and gas reservoir
Pore-size distribution problem realizes the global analysis of pore throat character feature in fine and close oil and gas reservoir, so that it is micro- to improve compact reservoir
The globality and accuracy of view hole gap structure characterization.
Disclosed above is only several specific embodiments of the invention, and still, the embodiment of the present invention is not limited to this, is appointed
What what those skilled in the art can think variation should all fall into protection scope of the present invention.
Claims (3)
1. a kind of method for determining compact reservoir pore-size distribution based on pressure mercury-N2 adsorption translocation data, which is characterized in that including such as
Lower step:
Step (1) determines the N2 adsorption parameter for the pore diameter range that nitrogen adsorption method can detect, uses n1Indicate that N2 adsorption method measures altogether
Aperture number, RiIndicate the size in i-th of aperture, αiIndicate frequency shared by i-th of aperture, XiIt indicates shared by i-th of aperture
Volume, V1Indicate the total pore volume that N2 adsorption method measures, i=1,2 ... n1；
Step (2) determines the pressure mercury parameter for the pore diameter range that mercury injection method can detect, uses n2Indicate that mercury method measures aperture number altogether,
rjIndicate the size in j-th of aperture, βjIndicate frequency shared by j-th of aperture, YjIndicate volume shared by j-th of aperture, V2It indicates
The total pore volume that pressure mercury method measures, j=1,2 ..., n2；
Step (3), the range measured due to two methods of nitrogen adsorption method, mercury injection method some be overlapping, and lap
The volume of pore throat be it is certain, then have Xi=Yj, it may be assumed that
V1×αi=V2×βj... ... ... ... ... formula is 1.
The two groups of throat distribution parameters measured by nitrogen adsorption method and mercury injection method directly read the number in overlapping aperture, indicate nitrogen with x
The number in the overlapping aperture that absorption method and mercury injection method measure；
Step (4) determines the frequency α in integration step (1)iWith the frequency β in step (2)jRequired mapping ratio；If λ is indicated
Mapping ratio needed for integration process, then λ is also that mercury injection method measures large pore and measures the total of small aperture with nitrogen adsorption method
2. the size of the ratio of body prior probability, λ is determined by formula:
2. 1. simultaneous formula is obtained with formula:
λtIt indicates to calculate obtained overall prior probability ratio by t-th of lap；Wherein i=1,2 ... n1；J=1,2 ...,
n2；T=1,2 ..., x；
It 3. finds out x prior probability ratio altogether by formula, formula is passed through by the obtained x prior probability ratio of this x overlapping interval
4. determining average prior probability ratio:
Step (5) determines the frequency α in integration step (1)iWith the frequency β in step (2)jFrequency afterwards；After indicating integration
Each aperture frequency, then integrate after each aperture frequencyValue pass through following process determine:
I, when pore size in the range of micropores of nitrogen absorption measurement and when, frequencyIt is 5. determined by formula:
II, when mesopore range of the pore size in the lap that two methods measure, frequencyIt is 6. determined by formula:
Wherein, work as i=n1When-x+1, j=1,2 ..., work as i=n1When, j=x；
III, when pore size is in the macropore range that mercury injection method measures, frequencyIt is 7. determined by formula:
Wherein, j=x+1, x+2 ..., n2；
Micropore size≤the 2nm, the pore diameter range of mesoporous are 2-50nm, macropore diameter >=50nm,
Step (6), what is obtained by step (5)SummationUnitization transformation is carried out, that is to say the hole so that each range
The sum of diameter distribution frequency becomes the ratio that 100, Ψ indicates unitization, 8. the size of Ψ is determined by formula:
Using formula 8. in Ψ size by formula 9. determination finally obtain the different pore size measured comprising two kinds of experimental methods
The whole throat distribution frequency of range, αkIndicate final throat distribution frequency,
Wherein k=1,2 ... n1+n2-x；
Pore-size distribution histogram is made using finally obtained data, the pore-size distribution that can be visually observed that in gamut
Situation.
2. a kind of method for determining compact reservoir pore-size distribution based on pressure mercury-N2 adsorption translocation data as described in claim 1,
It is characterized in that, in step (1), according to " rock specific surface and pore-size distribution measure static N2 adsorption volumetric method " SY/T 6154-
Process as defined in 1995 standards carries out N2 adsorption experiment to rock core, determines the pore-size distribution for the pore diameter range that nitrogen adsorption method can detect
Parameter.
3. a kind of method for determining compact reservoir pore-size distribution based on pressure mercury-N2 adsorption translocation data as described in claim 1,
It is characterized in that, in step (2), according to " measurement of rock capillary pressure curve " SY/T 5346-2005 normal process to rock core into
Row pressure mercury experiment, determines the pressure mercury parameter for the pore diameter range that mercury injection method can detect.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710238534.1A CN107132171B (en) | 2017-04-13 | 2017-04-13 | A method of compact reservoir pore-size distribution is determined based on pressure mercury-N2 adsorption translocation data |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710238534.1A CN107132171B (en) | 2017-04-13 | 2017-04-13 | A method of compact reservoir pore-size distribution is determined based on pressure mercury-N2 adsorption translocation data |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107132171A CN107132171A (en) | 2017-09-05 |
| CN107132171B true CN107132171B (en) | 2019-10-25 |
Family
ID=59716345
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710238534.1A Expired - Fee Related CN107132171B (en) | 2017-04-13 | 2017-04-13 | A method of compact reservoir pore-size distribution is determined based on pressure mercury-N2 adsorption translocation data |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107132171B (en) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109632594B (en) * | 2017-10-27 | 2022-03-18 | 东北石油大学 | Method for characterizing pore throat characteristics of tight reservoir on basis of high-pressure mercury intrusion multi-scale |
| CN108268712B (en) * | 2018-01-05 | 2021-05-11 | 中国石油大学(北京) | Method and device for determining capillary pressure of pore medium by nuclear magnetic resonance |
| CN111141650B (en) * | 2019-12-11 | 2021-03-30 | 中国地质大学（武汉） | Method for representing occurrence state of water in shale by utilizing frozen high-pressure mercury press |
| CN111175214A (en) * | 2020-01-17 | 2020-05-19 | 西安石油大学 | Method for representing full size of pore diameter of unconventional tight reservoir |
| CN111398122A (en) * | 2020-04-03 | 2020-07-10 | 中国矿业大学 | Comprehensive characterization method for heterogeneity characteristics of full-scale pore structure of shale |
| CN112326523A (en) * | 2020-09-18 | 2021-02-05 | 中国石油天然气股份有限公司 | Method and device for acquiring pore throat characteristic parameters, computer equipment and storage medium |
| CN112067531B (en) * | 2020-10-15 | 2021-08-24 | 西安特种设备检验检测院 | Porous material pore size distribution testing method |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP6233710B2 (en) * | 2014-04-22 | 2017-11-22 | 東海カーボン株式会社 | Carbon black and rubber composition |
| FR3024662B1 (en) * | 2014-08-05 | 2016-12-09 | Ifp Energies Now | METHOD FOR SEPARATING META-XYLENE USING EXTERNAL HIGH SURFACE ZEOLITHIC ADSORBENT |
| CN105486621B (en) * | 2015-12-29 | 2018-05-04 | 中国石油天然气股份有限公司 | A kind of shale pore-size distribution test method |
| CN105424580B (en) * | 2016-01-14 | 2018-10-02 | 太原理工大学 | A kind of coal full aperture measures and its hole shape semidefinite quantization method |
| CN106525691A (en) * | 2016-12-09 | 2017-03-22 | 河南理工大学 | Method for determining full-pore-diameter pore structure of coal through multi-data fusion |
-
2017
- 2017-04-13 CN CN201710238534.1A patent/CN107132171B/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| CN107132171A (en) | 2017-09-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107132171B (en) | A method of compact reservoir pore-size distribution is determined based on pressure mercury-N2 adsorption translocation data | |
| CN103278436B (en) | Quantitative characterization method of low penetration double-medium sandstone oil reservoir microscopic aperture structure | |
| Xie et al. | Fractal and multifractal analysis of carbonate pore-scale digital images of petroleum reservoirs | |
| CN101725344B (en) | Method for determining litho-electric parameters | |
| Clerke | Permeability, Relative Permeability, Microscopic Displacement Efficiency and Pore Geometry of M_1 Bimodal Pore Systems in Arab D Limestone | |
| CN104777181B (en) | Method for determining nuclear magnetic resonance T2 cutoff value and fluid saturation of tight oil, and apparatus thereof | |
| CN109932297B (en) | Method for calculating permeability of tight sandstone reservoir | |
| CN104619952A (en) | Digital rock analysis systems and methods with reliable multiphase permeability determination | |
| CN104819923A (en) | Low-permeability sandstone reservoir pore structure quantitative inversion method based on nuclear magnetic resonance | |
| CN105445442B (en) | Portray mud shale pore character and its method with shale oil enrichment relationship | |
| CN108798657B (en) | Logging interpretation method based on drilling fluid logging parameter gas measurement value | |
| CN112284999B (en) | Sandstone pore size distribution determination method and application thereof | |
| CN105447298A (en) | Reservoir stratum analysis method and application thereof | |
| CN110700820A (en) | Dessert classification method for compact oil reservoir in northern part of Songliao basin | |
| CN110702580A (en) | Compact sandstone reservoir pore throat heterogeneity characterization method based on information entropy | |
| CN109001093A (en) | A kind of quantitatively characterizing method of shale reservoir hole system and the method for evaluating shale reservoir pore character | |
| CN110632677B (en) | Carbonate rock geothermal reservoir classification method | |
| Lyu et al. | The effect of particle size on the interpretation of pore structure of shale by N2 adsorption | |
| CN104991280A (en) | Perforated stratum selection method and apparatus based on microscopic pore throat structure index | |
| CN107843531A (en) | The nuclear magnetic resonance parameter characterizing method of the low viscous oil oil reservoir pore structure of hyposmosis | |
| CN106321087B (en) | A method of obtaining lithostratigraphy factor | |
| CN104781648A (en) | Reservoir pore structure classification method based on bulb model | |
| Zhu et al. | Multi-scale characterization of organic matter pore space in deep marine shale combined with mathematical morphology | |
| CN110909466A (en) | Method for improving calculation precision of differential pore sorting reservoir average capillary pressure curve | |
| Song et al. | Pore fractal characteristics of lignite at different temperatures based on mercury intrusion test |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20191025 Termination date: 20200413 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |