CN107132171A - A kind of method for determining compact reservoir pore-size distribution based on pressure mercury N2 adsorption translocation data - Google Patents
A kind of method for determining compact reservoir pore-size distribution based on pressure mercury N2 adsorption translocation data Download PDFInfo
- Publication number
- CN107132171A CN107132171A CN201710238534.1A CN201710238534A CN107132171A CN 107132171 A CN107132171 A CN 107132171A CN 201710238534 A CN201710238534 A CN 201710238534A CN 107132171 A CN107132171 A CN 107132171A
- Authority
- CN
- China
- Prior art keywords
- pore
- aperture
- frequency
- formula
- represent
- 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
- 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 method for determining compact reservoir pore-size distribution based on pressure mercury N2 adsorption translocation data, it obtains the pore-size distribution situation in the reservoir first with N2 adsorption method, and N2 adsorption institute gaging hole larynx is distributed as the distribution situation of mesopore or micropore;Secondly using pressing mercury method to obtain the pore-size distribution situation in the reservoir, mercury injection method institute gaging hole larynx is distributed as mesopore or macropore;The pore-size distribution result for reusing the two measurement carries out throat distribution situation reallocation, asks for each interval aperture proportion, and then obtain the pore-size distribution situation of gamut.The problem of present invention solves different stage throat distribution in fine and close oil and gas reservoir, realizes the global analysis of pore throat character distribution characteristics in reservoir, so as to improve the globality and accuracy of compact reservoir micropore structure sign.
Description
Technical field
The invention belongs to petroleum natural gas exploration field, it is related to the research method of reservoir full aperture distribution, is specifically
A kind of method for determining compact reservoir pore-size distribution based on pressure mercury-N2 adsorption translocation data.
Background technology
With continuing to develop for oil-gas exploration, it is more next that the observation of the microcosmic reservoir pore throat character feature of compact reservoir describes method
Abundanter, from traditional Mercury injection combination Microscopic observation, to CT three-dimensional pore space analogue techniques, the research of reservoir micro throat is constantly
It is combined with mathematics and computer, and tends to quantification.
In the research method of existing reservoir pores distribution, because different measuring method principles are different, and experiment condition
Limitation, the pore-size distribution interval range of every kind of experimental method measurement can all be not quite similar.Mercury injection method is under the conditions of different pressures
By measuring the intrusion that mercury is pressed into rock sample hole, so as 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 its pore radius scope measured, which is generally in mesopore or micropore, mercury injection method, many
Tiny hole is no matter great pressure also can not be pressed into mercury solution, thus can only measure the aperture point of mesopore and macropore
Cloth (in the world, is typically divided into three classes aperture by size:Aperture≤2nm is micropore, and aperture is in 2-50nm scopes are
Hole, aperture >=50nm is macropore).And the distribution of rock sample microscopic aperture is unlikely to be pore size distribution or mesopore micropore point in single macropore
Cloth, two kinds of measuring methods are unable to be visually observed that the pore-size distribution situation in gamut, individually a kind of method measurement meeting
Cause measurement range narrow, be unfavorable for comprehensively researching and analysing reservoir micropore structure, particularly in compact reservoir
In research process, micropore macropore no less important should more increase the research to pore throat character holistic distribution characteristics.
The content of the invention
The purpose of the present invention is to overcome reservoir pores narrow distribution range present in above-mentioned prior art, can not intuitively observe
The problem of pore-size distribution situation there is provided a kind of based on pressure mercury-N2 adsorption translocation data method for determining compact reservoir pore-size distribution,
The pore size distribution range for pressing mercury and N2 adsorption method translocation to obtain is unified, currently known reservoir pressure mercury and N2 adsorption it can provide in fact
It is final to determine compact reservoir gamut pore-size distribution feature by analyzing the pore-size distribution feature of the two in the case of material, with behaviour
Control is simple, it is visual by force with accuracy it is high the characteristics of.
The technical scheme is that:It is a kind of to determine compact reservoir pore-size distribution based on pressure mercury-N2 adsorption translocation data
Method, comprises the following steps:
Step (1), determines the N2 adsorption parameter for the pore diameter range that nitrogen adsorption method can be detected, uses n1Represent that N2 adsorption method is total to
The aperture number measured, RiRepresent the size in the aperture, αiRepresent the frequency shared by the aperture, XiRepresent the aperture institute
The volume accounted for, w1Represent the total pore volume that N2 adsorption method is measured, i=1,2 ... n1。
Step (2), determines the pressure mercury parameter for the pore diameter range that mercury injection method can be detected, uses n2Represent that mercury method measures aperture altogether
Number, rjRepresent the size in the aperture, βjRepresent the frequency shared by the aperture, YjRepresent volume, w shared by the aperture2Table
Show the total pore volume that pressure mercury method is measured, j=1,2 ..., n2。
Step (3), due to nitrogen adsorption method, two methods of mercury injection method measure scope some be overlapping, and overlapping
The volume of part pore throat is certain, then has Xi=Yj, i.e.,:
ω1×αi=ω2×βj... ... ... ... ... formula is 1.
Therefore the variation tendency of the throat distribution for the lap that two methods are measured is identical.By nitrogen adsorption method and
Two groups of throat distribution parameters that mercury injection method is measured can directly read the number in overlapping aperture, and nitrogen adsorption method and mercury injection method are represented with x
The number in the overlapping aperture measured.
Step (4), determines the frequency α in integration step (1)iWith the frequency β in step (2)jRequired mapping ratio.If λ
The mapping ratio needed for integration process is represented, then also as mercury injection method measures large pore to λ and nitrogen adsorption method measures small aperture
Overall prior probability ratio, 2. λ size determined by formula:
1. 2. simultaneous formula obtain with formula:
λtRepresent to calculate resulting 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. obtained altogether by formula, by x prior probability obtained by this x overlapping interval than logical
Cross formula and 4. determine average prior probability ratio:
Step (5), determines the frequency α in integration step (1)iWith the frequency β in step (2)jFrequency afterwards.Represent whole
The frequency in each aperture after conjunction, the then frequency in each aperture after integratingValue determined by following process:
I, when pore size the range of micropores of nitrogen absorption measurement and when the aperture of lap (do not include), frequencyBy
5. formula determines:
II, when the mesopore range for the lap that pore size is measured in two methods, frequency6. determined by formula:
(wherein, work as i=n1During-x+1, j=1 ... ... works as i=n1When, j=x)
III, (aperture for not including lap), frequency when pore size is in the macropore scope that mercury injection method is measuredBy
7. formula determines:
Step (6), what is obtained by step (5)SummationUnitization conversion is carried out, namely is so that each scope
Pore-size distribution frequency sum be changed into the ratio that 100, Ψ represents unitization, 8. Ψ size is determined by formula:
9. determine to finally give the difference measured comprising two kinds of experimental methods by formula using formula 8. middle Ψ size
The overall throat distribution frequency of pore diameter range, αkFinal throat distribution frequency is represented,
Pore-size distribution block diagram is made using the data finally given, the aperture in gamut just can be visually observed that
Distribution situation, is conducive to the research to pore throat character holistic distribution characteristics.
It is preferred that in step (1), foundation《Rock specific surface and pore-size distribution determine static nitrogen adsorption capacity method》SY/T
Flow as defined in 6154-1995 standards carries out N2 adsorption experiment to rock core, determines the hole for the pore diameter range that nitrogen adsorption method can be detected
Footpath distributed constant.
It is preferred that in step (2), foundation《Rock capillary pressure curve is determined》SY/T 5346-2005 normal process is to rock
The heart carries out pressure mercury experiment, determines the pressure mercury parameter for the pore diameter range that mercury injection method can be detected.
Beneficial effects of the present invention:Determine to cause based on pressure mercury-N2 adsorption translocation data there is provided one kind in the embodiment of the present invention
The method of close reservoir pores distribution, tests the small-bore measured distribution by N2 adsorption and intraventricular pressure mercury tests 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, because a certain overlapping interval volume that two methods are measured is identical, therefore can obtain pressure mercury by calculating
The macrovoid cumulative volume that method measurement the is obtained prior probability ratio overall with the fine pore that nitrogen adsorption method measurement is obtained, wherein priori
Likelihood ratio is tried to achieve respectively using the related data in each overlapping aperture in overlapping range, its average value is then taken again, most
Afterwards using required prior probability ratio, the pore-size distribution frequency that two kinds of measuring methods are obtained is integrated, and the aperture after integration point
The progress of cloth frequency is unitization, finally gives the graph of pore diameter distribution of each scope.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, so as to improve
Globality and accuracy that compact reservoir micropore structure is characterized.
Brief description of the drawings
Fig. 1 is flow chart of the method for the present invention;
Fig. 2 is the nitrogen adsorption method graph of pore diameter distribution that the present invention is provided;
Fig. 3 is the mercury injection method graph of pore diameter distribution that the present invention is provided;
Fig. 4 is the pressure mercury and nitrogen adsorption method graph of pore diameter distribution that the present invention is provided;
Fig. 5 is the pressure finally given mercury-N2 adsorption translocation gamut pore-size distribution that the present invention is provided.
Embodiment
Below in conjunction with the accompanying drawings, the specific embodiment of the present invention is described in detail, it is to be understood that of the invention
Protection domain do not limited by embodiment.
Referring to Fig. 1, compact reservoir aperture is determined based on pressure mercury-N2 adsorption translocation data the embodiments of the invention provide one kind
The method of distribution, comprises the following steps:
Step (1), determines the N2 adsorption parameter for the pore diameter range that nitrogen adsorption method can be detected, uses n1Represent that N2 adsorption method is total to
The aperture number measured, RiRepresent the size in the aperture, αiRepresent the frequency shared by the aperture, XiRepresent the aperture institute
The volume accounted for, w1Represent the total pore volume that N2 adsorption method is measured, i=1,2 ... n1。
Step (2), determines the pressure mercury parameter for the pore diameter range that mercury injection method can be detected, uses n2Represent that mercury method measures aperture altogether
Number, rjRepresent the size in the aperture, βjRepresent the frequency shared by the aperture, YjRepresent volume, w shared by the aperture2Table
Show the total pore volume that pressure mercury method is measured, j=1,2 ..., n2。
Step (3), due to nitrogen adsorption method, two methods of mercury injection method measure scope some be overlapping, and overlapping
The volume of part pore throat is certain, then has Xi=Yj, i.e.,:
ω1×αi=ω2×βj... ... ... ... ... formula is 1.
Therefore the variation tendency of the throat distribution for the lap that two methods are measured is identical.By nitrogen adsorption method and
Two groups of throat distribution parameters that mercury injection method is measured can directly read the number in overlapping aperture, and nitrogen adsorption method and mercury injection method are represented with x
The number in the overlapping aperture measured.
Step (4), determines the frequency α in integration step (1)iWith the frequency β in step (2)jRequired mapping ratio.If λ
The mapping ratio needed for integration process is represented, then also as mercury injection method measures large pore to λ and nitrogen adsorption method measures small aperture
Overall prior probability ratio, 2. λ size determined by formula:
1. 2. simultaneous formula obtain with formula:
λtRepresent to calculate resulting 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. obtained altogether by formula, by x prior probability obtained by this x overlapping interval than logical
Cross formula and 4. determine average prior probability ratio:
Step (5), determines the frequency α in integration step (1)iWith the frequency β in step (2)jFrequency afterwards.Represent whole
The frequency in each aperture after conjunction, the then frequency in each aperture after integratingValue determined by following process:
I, when pore size the range of micropores of nitrogen absorption measurement and when the aperture of lap (do not include), frequencyBy
5. formula determines:
II, when the mesopore range for the lap that pore size is measured in two methods, frequency6. determined by formula:
(wherein, work as i=n1During-x+1, j=1 ... ... works as i=n1When, j=x)
III, (aperture for not including lap), frequency when pore size is in the macropore scope that mercury injection method is measuredBy
7. formula determines:
Step (6), what is obtained by step (5)SummationUnitization conversion is carried out, namely is so that each scope
Pore-size distribution frequency sum be changed into the ratio that 100, Ψ represents unitization, 8. Ψ size is determined by formula:
9. determine to finally give the difference measured comprising two kinds of experimental methods by formula using formula 8. middle Ψ size
The overall throat distribution frequency of pore diameter range, αkFinal throat distribution frequency is represented,
Pore-size distribution block diagram is made using the data finally given, the aperture in gamut just can be visually observed that
Distribution situation, is conducive to the research to pore throat character holistic distribution characteristics.
Further, in step (1), foundation《Rock specific surface and pore-size distribution determine static nitrogen adsorption capacity method》SY/T
Flow as defined in 6154-1995 standards carries out N2 adsorption experiment to rock core, determines the hole for the pore diameter range that nitrogen adsorption method can be detected
Footpath distributed constant.
Further, in step (2), foundation《Rock capillary pressure curve is determined》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 be detected.
Definitely to illustrate actual effect of the present invention, now choose Song-liao basin oil field tight sand sample and (bury
Deep 2342.5m) it is used as the research object of this embodiment.
Step one, N2 adsorption experimental method pore diameter range is determined:Choose representative core sample, foundation《Rock
Static nitrogen adsorption capacity method SY/T 6154-1995 are determined than surface and pore-size distribution》Flow as defined in standard carries out nitrogen to rock core
Adsorption experiment obtains the N2 adsorption parameter of mesopore and micropore, n1Represent the aperture number that N2 adsorption method is measured altogether, n1=19,
RiRepresent the size in the aperture, αiRepresent the frequency shared by the aperture, XiRepresent the volume shared by the aperture, w1Represent
The total pore volume that N2 adsorption method is measured, i=1,2 ... 19, specific experiment the results are shown in Table 1, and graph of pore diameter distribution is made by table 1
(Fig. 2).
The nitrogen adsorption method pore size distribution data of table 1
Sequence 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 |
Sequence 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 is determined:Foundation《Rock capillary pressure curve is determined》SY/T 5346-
2005 normal process carry out the pressure mercury parameter that the experiment of rock core mercury injection method obtains macropore and mesopore, and aperture number n is measured altogether2=11,
The size r in the aperture measuredjAnd shared frequency βjIt is shown in Table 2, YjVolume shared by expression the aperture, j=1,2 ..., 11,
w2Represent the total pore volume that pressure mercury method is measured.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).
The mercury injection method pore size distribution data of table 2
Step 3, has Fig. 4 analyses to understand that the interval pore-size distribution trend of the lap that two methods are measured is identical, only
It is that proportion is different in respective measurement range, actual is with identical pore throat space.According in the data of table 1 and table 2
Data determine that nitrogen adsorption method is 6, the pore size distribution data of lap with the overlapping number for the pore diameter range that mercury injection method is determined
It is shown in Table 3.
The lap pore size distribution data of table 3
Aperture/nm | 6.3 | 10 | 16 | 25 | 40 | 63 |
Frequency αi/ % | 8.322 | 11.967 | 9.913 | 16.516 | 25.333 | 11.408 |
Frequency βj/ % | 3.219 | 5.926 | 6.659 | 8.34 | 8.852 | 10.531 |
Step 4, determines the frequency α in integration step (1)iWith the frequency β in step (2)jRequired mapping ratio.If λ
The mapping ratio needed for integration process is represented, 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 w that λ size is measured by the pressure mercury method in step (2)2Size and step
(1) the total pore volume w that the N2 adsorption method in is measured1Size determine, i.e.,:
Again because when positioned at a certain overlapping interval, this overlapping interval volume and mercury injection method of determination of nitrogen adsorption are determined
The volume of this overlapping interval be equal, then have:
Xi=Yj, i.e. ω1×αi=ω2×βj... ... ... ... formula is 2.
By formula, 1. 2. simultaneous is obtained with formula:
3. formula is the evaluation formula of prior probability ratio, λtRepresent to calculate resulting total hole by t-th of lap
Gap volume prior probability ratio, as i=14, j=1, t=1, pore diameter range are located at the 1st overlapping interval, by that analogy, i=19
When, j=6, t=6, pore diameter range are 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, 6 prior probabilities, respectively λ are obtained1=2.585, λ2=2.019, λ3=
1.488, λ4=1.980, λ5=2.862, λ6=1.083, pass through formula by 6 prior probabilities obtained by this 6 overlapping intervals
4. average prior probability ratio is calculated:
Obtained average prior probability ratio is:
Step 5, integrates N2 adsorption and tests the frequency of the micropore measured and mesopore with pressing the mesopore and macropore that mercury experiment is measured
Frequency.The frequency in each aperture after integrating is represented,Value determined by following process:
I, when aperture is in 0.015-6.3nm scopes, the frequency in each aperture of this scope5. determined by formula:
II, when aperture is in 6.3-63nm scopes, the frequency in each aperture of this scope6. determined by formula:
III, when aperture is in 63-630nm scopes, the frequency in each aperture of this scope7. 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 conversion is carried out, namely is so that each scope
Pore-size distribution frequency sum should be 100.Unitization ratio is Ψ, then has
It is final to calculate the different pore size scope throat distribution frequency α for obtaining measuring comprising two kinds of experimental methodsk, then haveThe different pore size scope that two methods are measured can now be obtained
Pore-size distribution frequency αk, as shown in table 5.
The final pore size distribution data of table 5
The a wide range of interior different interval pore-size distribution frequencies finally obtained, are organized into figure by this composite result, that is, obtain
The a wide range of pore-size distribution situation that two methods are measured.It thus can no longer be bound to the limit that experimental measurement method is brought
System, obtains nanoscale to micron-sized reservoir pores distribution situation, such as Fig. 5.
In summary, it is provided in an embodiment of the present invention to determine compact reservoir aperture point based on pressure mercury-N2 adsorption translocation data
The method of cloth, tests the small-bore measured distribution by N2 adsorption and intraventricular pressure mercury tests the large aperture measured and is 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
The a certain overlapping interval volume measured in two methods is identical, therefore can obtain what mercury injection method measurement was obtained by calculating
The macrovoid cumulative volume prior probability ratio overall with the fine pore that nitrogen adsorption method measurement is obtained, wherein prior probability ratio are to utilize
The related data in each overlapping aperture is tried to achieve respectively in overlapping range, its average value is then taken again, finally utilizes required elder generation
Probability ratio is tested, the pore-size distribution frequency that two kinds of measuring methods are obtained is integrated, and the pore-size distribution frequency after integration carries out list
Positionization, finally gives the graph of pore diameter distribution of each scope.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
Globality and accuracy that view hole gap structure is characterized.
Disclosed above is only several specific embodiments of the present 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 change 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, it is characterised in that including such as
Lower step:
Step (1), determines the N2 adsorption parameter for the pore diameter range that nitrogen adsorption method can be detected, uses n1Represent that N2 adsorption method is measured altogether
Aperture number, RiRepresent the size in the aperture, αiRepresent the frequency shared by the aperture, XiRepresent shared by the aperture
Volume, w1Represent the total pore volume that N2 adsorption method is measured, i=1,2 ... n1;
Step (2), determines the pressure mercury parameter for the pore diameter range that mercury injection method can be detected, uses n2Represent that mercury method measures aperture number altogether,
rjRepresent the size in the aperture, βjRepresent the frequency shared by the aperture, YjRepresent volume, w shared by the aperture2Represent pressure
The total pore volume that mercury method is measured, j=1,2 ..., n2;
Step (3), the scope measured due to two methods of nitrogen adsorption method, mercury injection method some be overlapping, and lap
The volume of pore throat is certain, then has Xi=Yj, i.e.,:
ω1×αi=ω2×βj... ... ... ... ... formula is 1.
The two groups of throat distribution parameters measured by nitrogen adsorption method and mercury injection method can directly read the number in overlapping aperture, be represented with x
The number in the overlapping aperture that nitrogen adsorption method and mercury injection method are measured;
Step (4), determines the frequency α in integration step (1)iWith the frequency β in step (2)jRequired mapping ratio;If λ is represented
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 ratio of body prior probability, λ size is determined by formula:
1. 2. simultaneous formula obtain with formula:
λtRepresent to calculate resulting overall prior probability ratio by t-th of lap;Wherein i=1,2 ... n1;J=1,2 ...,
n2;T=1,2 ..., x;
3. obtain x prior probability ratio altogether by formula, formula is passed through by x prior probability ratio obtained by this x overlapping interval
4. average prior probability ratio is determined:
Step (5), determines the frequency α in integration step (1)iWith the frequency β in step (2)jFrequency afterwards;Represent after integrating
Each aperture frequency, then integrate after each aperture frequencyValue determined by following process:
I, when pore size in the range of micropores of nitrogen absorption measurement and when, frequency5. determined by formula:
II, when the mesopore range for the lap that pore size is measured in two methods, frequency6. determined by formula:
Wherein, i=n is worked as1During-x+1, j=1,2 ..., work as i=n1When, j=x;
III, when pore size is in the macropore scope that mercury injection method is measured, frequency7. determined by formula:
Wherein, j=x+1, x+2 ..., n2;
Step (6), what is obtained by step (5)SummationUnitization conversion is carried out, namely is so that the hole of each scope
Footpath distribution frequency sum is changed into the ratio that 100, Ψ represents unitization, and 8. Ψ size is determined by formula:
9. determine to finally give the different pore size measured comprising two kinds of experimental methods by formula using formula 8. middle Ψ size
The overall throat distribution frequency of scope, αkFinal throat distribution frequency is represented,
Wherein k=1,2 ... n1+n2-x;
Pore-size distribution block diagram is made using the data finally given, the pore-size distribution in gamut just can be visually observed that
Situation.
2. a kind of method for determining compact reservoir pore-size distribution based on pressure mercury-N2 adsorption translocation data as claimed in claim 1,
Characterized in that, in step (1), foundation《Rock specific surface and pore-size distribution determine static nitrogen adsorption capacity method》SY/T 6154-
Flow 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 be detected
Parameter.
3. a kind of method for determining compact reservoir pore-size distribution based on pressure mercury-N2 adsorption translocation data as claimed in claim 1,
Characterized in that, in step (2), foundation《Rock capillary pressure curve is determined》SY/T 5346-2005 normal process enters to rock core
Row pressure mercury experiment, determines the pressure mercury parameter for the pore diameter range that mercury injection method can be detected.
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 true CN107132171A (en) | 2017-09-05 |
CN107132171B 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) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108268712A (en) * | 2018-01-05 | 2018-07-10 | 中国石油大学(北京) | The method and apparatus that nuclear magnetic resonance determines pore media capillary pressure |
CN109632594A (en) * | 2017-10-27 | 2019-04-16 | 东北石油大学 | A method of based on high-pressure mercury Multi-scale characterization compact reservoir pore throat characteristic |
CN111141650A (en) * | 2019-12-11 | 2020-05-12 | 中国地质大学(武汉) | 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 |
CN112067531A (en) * | 2020-10-15 | 2020-12-11 | 西安特种设备检验检测院 | Porous material pore size distribution testing method |
CN112326523A (en) * | 2020-09-18 | 2021-02-05 | 中国石油天然气股份有限公司 | Method and device for acquiring pore throat characteristic parameters, computer equipment and storage medium |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105424580A (en) * | 2016-01-14 | 2016-03-23 | 太原理工大学 | Coal full-pore-diameter measuring and pore shape semi-quantifying method |
CN105486621A (en) * | 2015-12-29 | 2016-04-13 | 中国石油天然气股份有限公司 | Test method for shale pore size distribution |
EP3135735A1 (en) * | 2014-04-22 | 2017-03-01 | Tokai Carbon Co., Ltd. | Carbon black and rubber composition |
CN106525691A (en) * | 2016-12-09 | 2017-03-22 | 河南理工大学 | Method for determining full-pore-diameter pore structure of coal through multi-data fusion |
US20170226031A1 (en) * | 2014-08-05 | 2017-08-10 | IFP Energies Nouvelles | Method for separating meta-xylene using a zeolitic adsorbent with a large external surface area |
-
2017
- 2017-04-13 CN CN201710238534.1A patent/CN107132171B/en not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3135735A1 (en) * | 2014-04-22 | 2017-03-01 | Tokai Carbon Co., Ltd. | Carbon black and rubber composition |
US20170226031A1 (en) * | 2014-08-05 | 2017-08-10 | IFP Energies Nouvelles | Method for separating meta-xylene using a zeolitic adsorbent with a large external surface area |
CN105486621A (en) * | 2015-12-29 | 2016-04-13 | 中国石油天然气股份有限公司 | Test method for shale pore size distribution |
CN105424580A (en) * | 2016-01-14 | 2016-03-23 | 太原理工大学 | Coal full-pore-diameter measuring and pore shape semi-quantifying method |
CN106525691A (en) * | 2016-12-09 | 2017-03-22 | 河南理工大学 | Method for determining full-pore-diameter pore structure of coal through multi-data fusion |
Non-Patent Citations (1)
Title |
---|
张云峰等: "基于氮吸附和压汞数据确定致密储层孔径分布", 《深圳大学学报理工版》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109632594A (en) * | 2017-10-27 | 2019-04-16 | 东北石油大学 | A method of based on high-pressure mercury Multi-scale characterization compact reservoir pore throat characteristic |
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 |
CN108268712A (en) * | 2018-01-05 | 2018-07-10 | 中国石油大学(北京) | The method and apparatus that nuclear magnetic resonance determines pore media capillary pressure |
CN108268712B (en) * | 2018-01-05 | 2021-05-11 | 中国石油大学(北京) | Method and device for determining capillary pressure of pore medium by nuclear magnetic resonance |
CN111141650A (en) * | 2019-12-11 | 2020-05-12 | 中国地质大学(武汉) | 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 |
CN112067531A (en) * | 2020-10-15 | 2020-12-11 | 西安特种设备检验检测院 | Porous material pore size distribution testing method |
Also Published As
Publication number | Publication date |
---|---|
CN107132171B (en) | 2019-10-25 |
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 | |
CN109932297B (en) | Method for calculating permeability of tight sandstone reservoir | |
Clerke | Permeability, Relative Permeability, Microscopic Displacement Efficiency and Pore Geometry of M_1 Bimodal Pore Systems in Arab D Limestone | |
CN101725344B (en) | Method for determining litho-electric parameters | |
CN110296931B (en) | Characterization method and system for oil-water relative permeability information of tight sandstone | |
CN109632594B (en) | Method for characterizing pore throat characteristics of tight reservoir on basis of high-pressure mercury intrusion multi-scale | |
CN104819923A (en) | Low-permeability sandstone reservoir pore structure quantitative inversion method based on nuclear magnetic resonance | |
CN103512838A (en) | Method for determining contribution of pores with different apertures in shale reservoir stratum to porosity | |
CN104374685B (en) | A kind of method that digital cores three hole component asks for reservoir permeability | |
CN104777181A (en) | Method for determining nuclear magnetic resonance T2 cutoff value and fluid saturation of tight oil, and apparatus thereof | |
CN105445442A (en) | Method for describing pore feature of shale and enrichment relationship between shale and shale oil | |
CN110487693A (en) | A kind of method of determining mud shale different type porosity | |
CN105447298A (en) | Reservoir stratum analysis method and application thereof | |
CN110161071A (en) | A method of evaluation compact reservoir movable fluid Minimum throat radius | |
CN110700820A (en) | Dessert classification method for compact oil reservoir in northern part of Songliao basin | |
CN107153733A (en) | A kind of capillary pressure curve approximating method of Complicated Pore Structures rock | |
CN106897531A (en) | A kind of method for quantitatively evaluating of hyposmosis limestone reservoir permeability | |
CN109001093A (en) | A kind of quantitatively characterizing method of shale reservoir hole system and the method for evaluating shale reservoir pore character | |
CN104237107A (en) | Method and system for interpreting apparent permeability of low-permeability reservoir in formation | |
Zhu et al. | Multi-scale characterization of organic matter pore space in deep marine shale combined with mathematical morphology | |
CN110909466B (en) | Method for improving calculation accuracy of average capillary pressure curve of differential pore sorting reservoir | |
CN107843531A (en) | The nuclear magnetic resonance parameter characterizing method of the low viscous oil oil reservoir pore structure of hyposmosis | |
CN108920824B (en) | Yield decrement analysis method for narrow river channel reservoir | |
CN106370577A (en) | Fracture cave oil reservoir simulation test method |
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 |