CN108195735A - Capillary pressure curve classification method - Google Patents
Capillary pressure curve classification method Download PDFInfo
- Publication number
- CN108195735A CN108195735A CN201711296710.3A CN201711296710A CN108195735A CN 108195735 A CN108195735 A CN 108195735A CN 201711296710 A CN201711296710 A CN 201711296710A CN 108195735 A CN108195735 A CN 108195735A
- Authority
- CN
- China
- Prior art keywords
- curves
- capillary pressure
- curve
- pressure curve
- pressure
- 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
Landscapes
- Chemical & Material Sciences (AREA)
- Dispersion Chemistry (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Measuring Volume Flow (AREA)
- Earth Drilling (AREA)
Abstract
The invention discloses a capillary pressure curve classification method, which comprises the steps of fitting a capillary pressure curve, and determining the number N of pore systems and displacement pressure P reflected by the capillary pressure curvedAnd according to PdSorting the curves by size; then comparing the similarity of any two curves to ensure that the two curves have the same N value, calculating the absolute value of the saturation difference corresponding to the two capillary pressure curves one by one, and if the saturation difference of all the points is less than the allowable error, dividing the two curves into the same category. The invention ensures the similarity of the overall shapes of the curves, can adjust the classification precision according to the requirements and has larger flexibility.
Description
Technical field
The present invention relates to a kind of capillary pressure curve sorting techniques, quantitatively divide for the capillary pressure curve of pore media
Analysis belongs to reservoir rock pore structure evaluation field.
Background technology
For reservoir rock pore configuration research, capillary pressure experiment is to study the basis of throat size, different hairs
Pipe pressure tracing pattern feature reflects the different pore structure characteristic of rock.In evaluating reservoir, usually according to capillary pressure
The morphological feature of curve classifies to reservoir, so as to improve the accuracy of evaluating reservoir.Therefore, capillary pressure curve is carried out
Accurately classification is the basis of reservoir classification and evaluation.
Traditional capillary pressure curve sorting technique is nearly all to utilize Pd（Replacement pressure）、R35（Mercury saturation degree 35% is right
The throat radius answered）、R50（50% corresponding throat radius of mercury saturation degree）And Swir（Irreducible water saturation）Wait curvilinear characteristics parameter
One or more of, classified by the method for cluster to it.These sorting techniques are there are clearly disadvantageous part, usually
The larger curve of morphological differences is divided into same class, especially the application effect in the stronger carbonate rock of anisotropism compared with
Difference.Such as Fig. 1（a）In two capillary pressure curves significantly there is different morphological features, but R35 and R50 values of two curves
It is all very close, if classified according only to R35 or R50 to it, the two can be divided into same class；Equally in Fig. 1（b）
In, according only to PdAnd SwirAlso the curve that two can be differed greatly is divided into same class.Even with Pd, R35, R50 and SwirDeng more
A parameter carries out compressive classification, it is also possible to which differ greatly two curves are divided into same class, such as Fig. 1（c）It is shown.
It can be seen that the prior art classifies to the classification of capillary pressure curve only in accordance with characteristic parameter, do not account for
The similitude of curve configuration, except accuracy, rapidity and the flexibility of low capillary pressure curve classification.
Invention content
It is an object of the invention to overcome the above problem of the existing technology, a kind of capillary pressure curve classification side is provided
Method.Present invention ensures that the similitude of curve configuration, and the precision of classification can be adjusted as needed, have larger flexible
Property.
To achieve the above object, the technical solution adopted by the present invention is as follows：
A kind of capillary pressure curve sorting technique, it is characterised in that：Capillary pressure curve is fitted first, determines hollow billet pressure
The pore system quantity N and replacement pressure P that force curve is reflectedd, and according to PdSize curve is ranked up；Then compare
The similitude of arbitrary two curves, it is ensured that two curves have identical N values, and pressure spot calculates two capillary pressure curves one by one
The absolute value of corresponding saturation degree difference, if the saturation degree difference of all the points is both less than allowable error, by two curves
It is divided into same category.
In the method, utilize（1）Formula is fitted all original capillary pressure curves one by one, determines each curve
The pore system quantity N and replacement pressure P reflectedd；
（1）
In formula：
Pc- note mercury pressure（Bar or Mpa）；
Sw- water saturation（Decimal）；
The quantity of N-pore system（Dimensionless）；
A, the coefficient of b, c, d-needs fitting（Dimensionless）；
tanh- hyperbolic tangent function.
In the method, by replacement pressure PdPrinciple from high to low, is ranked up capillary pressure curve, compiles successively
Number, and the category attribute of all curves is initialized as 0.
In the method, fluid saturation allowable error Err is determined, i.e., two similar curves are corresponding to arbitrary pressure spot
Fluid saturation absolute value of the difference maximum value.
The maximum absolute value value is bigger, and classification is thicker, and obtained class categories are fewer；Maximum absolute value value is smaller,
Classification is finer, and class categories are more.
The maximum absolute value value is determined according to the requirement of the quantity and nicety of grading of capillary pressure curve, if classification
It is relatively thick, then the value between 0.1 ~ 0.2；Such as sophisticated category, then the value between 0.05 ~ 0.1.
It is described that two curves are divided into same category, compared by loop iteration, allow each curve all with it is other
All curves carry out a similarity system design, if two curves have identical N values, and two songs at each pressure spot
The absolute value of fluid saturation difference corresponding to line is both less than allowable error Err, then two curves is divided into same category.
It is described two curves are divided into same category to be specially：In arbitrary pressure spot PciLocate, corresponding to two curves
Saturation degree is respectively SwA and SwB, and absolute difference is | SwB-SwA |, if all met at all pressure spots | SwB-SwA |<
Err then illustrates that two tracing patterns are similar.
Advantage using the present invention is：
1st, the present invention classifies according to the quantity of pore system, it is ensured that the capillary pressure with different aperture system quantity is bent
Line belongs to different classifications.
2nd, the present invention is by comparing the saturation difference corresponding to each pressure spot, it is ensured that the phase of curve configuration
Like property.
3rd, conventional method is unable to control the precision of classification, and method proposed by the present invention can adjust the essence of classification as needed
Degree has larger flexibility.
To sum up, method proposed by the present invention can carry out fast automatic classification to a large amount of capillary pressure curves, suitable for any
The capillary pressure curve of form has very big flexibility, is reservoir classification and evaluation, rock type division, pore configuration research
It lays a good foundation.
Description of the drawings
Fig. 1 a-c are existing general characteristics parametric classification schematic diagram；
Fig. 2 is the corresponding saturation degree mathematic interpolation schematic diagram of arbitrary pressure spot on capillary pressure curve.
Specific embodiment
Embodiment 1
A kind of capillary pressure curve sorting technique, is first fitted capillary pressure curve, determines that capillary pressure curve institute is anti-
The pore system quantity N and replacement pressure P reflectedd, and according to PdSize curve is ranked up；Then more arbitrary two songs
The similitude of line, it is ensured that two curves have identical N values, and pressure spot calculates satisfying corresponding to two capillary pressure curves one by one
With the absolute value of degree difference, if the saturation degree difference of all the points is both less than allowable error, two curves are divided into same
Classification.
In the method, utilize（1）Formula is fitted all original capillary pressure curves one by one, determines each curve
The pore system quantity N and replacement pressure P reflectedd；
（1）
In formula：
Pc- note mercury pressure（Bar or Mpa）；
Sw- water saturation（Decimal）；
The quantity of N-pore system（Dimensionless）；
A, the coefficient of b, c, d-needs fitting（Dimensionless）；
tanh- hyperbolic tangent function.
In the method, by replacement pressure PdPrinciple from high to low, is ranked up capillary pressure curve, compiles successively
Number, and the category attribute of all curves is initialized as 0.
In the method, fluid saturation allowable error Err is determined, i.e., two similar curves are corresponding to arbitrary pressure spot
Fluid saturation absolute value of the difference maximum value.
The maximum absolute value value is bigger, and classification is thicker, and obtained class categories are fewer；Maximum absolute value value is smaller,
Classification is finer, and class categories are more.
The maximum absolute value value is determined according to the requirement of the quantity and nicety of grading of capillary pressure curve, if classification
It is relatively thick, then the value between 0.1 ~ 0.2；Such as sophisticated category, then the value between 0.05 ~ 0.1.
It is described that two curves are divided into same category, compared by loop iteration, allow each curve all with it is other
All curves carry out a similarity system design, if two curves have identical N values, and two songs at each pressure spot
The absolute value of fluid saturation difference corresponding to line is both less than allowable error Err, then two curves is divided into same category.
It is described two curves are divided into same category to be specially：In arbitrary pressure spot PciLocate, corresponding to two curves
Saturation degree is respectively SwA and SwB, and absolute difference is | SwB-SwA |, if all met at all pressure spots | SwB-SwA |<
Err then illustrates that two tracing patterns are similar.
Embodiment 2
Capillary pressure curve is fluid saturation（Wetting phase or non-wetted phase）Pressure-dependent curve, the hole knot with rock
Structure feature is closely related.Two similar capillary pressure curves, the fluid under the conditions of the uniform pressure corresponding to two curves are satisfied
Should be very close with degree, illustrate that two curves do not have similitude, and reflect difference if fluid saturation differs greatly
Pore structure characteristic.Capillary pressure curve automatic classification method proposed by the present invention, the method recycled by double iterative, allows
Each curve all carries out a similarity system design with other all curves, and ensures that same class curve has identical N values.From
After the completion of dynamic classification, each curve has the category attribute value of an assigned integer（GP）, identical property value then belongs to
Same category.
Based on this principle, the automatic classification method of capillary pressure curve proposed by the present invention specifically includes following four step
Suddenly：
The first step utilizes（1）Formula is fitted all original capillary pressure curves one by one, determines what each curve was reflected
Pore system quantity（N）And replacement pressure（Pd）；
（1）
In formula：
Pc- note mercury pressure（Bar or Mpa）；
Sw- water saturation（Decimal）；
The quantity of N-pore system（Dimensionless）；
A, the coefficient of b, c, d-needs fitting（Dimensionless）；
tanh- hyperbolic tangent function.
Second step, by replacement pressure PdPrinciple from high to low, is ranked up capillary pressure curve, number consecutively, and
The category attribute of all curves is initialized as 0.
Third walks, and determines fluid saturation allowable error Err, i.e., two similar curves are in the stream corresponding to arbitrary pressure spot
The absolute value of the difference maximum value of body saturation degree.The value is bigger, and classification is thicker, and obtained class categories are fewer；The value is smaller, point
Class is finer, and class categories are more.Usually can the value be determined according to the requirement of the quantity and nicety of grading of capillary pressure curve,
If classification it is thicker, can between 0.1 ~ 0.2 value；Such as require sophisticated category, then can between 0.05 ~ 0.1 value.
4th step, loop iteration compare, and each curve is allowed all to carry out a similarity system design with other all curves.Such as
Two curves of fruit have identical N values, and the fluid saturation difference at each pressure spot corresponding to two curves is exhausted
Allowable error Err is both less than to value, then two curves are divided into same category.As shown in Fig. 2, in pressure spot PciPlace, two
Saturation degree corresponding to curve is respectively SwA and SwB, and absolute difference is | SwB-SwA |.It is if all full at all pressure spots
Foot | SwB-SwA |<Err then illustrates that two tracing patterns are similar.
Claims (8)
1. a kind of capillary pressure curve sorting technique, it is characterised in that：Capillary pressure curve is fitted first, determines hollow billet
The pore system quantity N and replacement pressure P that pressure curve is reflectedd, and according to PdSize curve is ranked up；Then compare
The similitude of more arbitrary two curves, it is ensured that two curves have identical N values, and pressure spot calculates two capillary pressure songs one by one
The absolute value of saturation degree difference corresponding to line, if the saturation degree difference of all the points is both less than allowable error, by two songs
Line is divided into same category.
2. capillary pressure curve sorting technique according to claim 1, it is characterised in that：In the method, utilize（1）Formula
All original capillary pressure curves are fitted one by one, determine that pore system quantity N that each curve reflected and row drive
Pressure Pd；
（1）
In formula：
Pc- note mercury pressure（Bar or Mpa）；
Sw- water saturation（Decimal）；
The quantity of N-pore system（Dimensionless）；
A, the coefficient of b, c, d-needs fitting（Dimensionless）；
tanh- hyperbolic tangent function.
3. capillary pressure curve sorting technique according to claim 2, it is characterised in that：In the method, pressure is driven by row
Power PdPrinciple from high to low, is ranked up capillary pressure curve, number consecutively, and the category attribute of all curves is initial
Turn to 0.
4. capillary pressure curve sorting technique according to claim 3, it is characterised in that：In the method, fluid is determined
Saturation degree allowable error Err, i.e. two similar curves the fluid saturation corresponding to arbitrary pressure spot absolute value of the difference most
Big value.
5. capillary pressure curve sorting technique according to claim 4, it is characterised in that：The maximum absolute value value is got over
Greatly, classify thicker, obtained class categories are fewer；Maximum absolute value value is smaller, and classification is finer, and class categories are more.
6. capillary pressure curve sorting technique according to claim 5, it is characterised in that：The maximum absolute value value according to
The requirement of the quantity and nicety of grading of capillary pressure curve determines, if classifying thicker, the value between 0.1 ~ 0.2；Such as
Sophisticated category, the then value between 0.05 ~ 0.1.
7. capillary pressure curve sorting technique according to claim 6, it is characterised in that：It is described to be divided into two curves
In same category, compared by loop iteration, each curve is allowed all to carry out a similarity system design with other all curves, such as
Two curves of fruit have identical N values, and the fluid saturation difference at each pressure spot corresponding to two curves is exhausted
Allowable error Err is both less than to value, then two curves are divided into same category.
8. capillary pressure curve sorting technique according to claim 7, it is characterised in that：It is described to be divided into two curves
Same category is specially：In arbitrary pressure spot PciPlace, the saturation degree corresponding to two curves are respectively SwA and SwB, difference
Absolute value is | SwB-SwA |, if all met at all pressure spots | SwB-SwA |<Err then illustrates that two tracing patterns are similar.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711296710.3A CN108195735B (en) | 2017-12-08 | 2017-12-08 | Capillary pressure curve classification method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711296710.3A CN108195735B (en) | 2017-12-08 | 2017-12-08 | Capillary pressure curve classification method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108195735A true CN108195735A (en) | 2018-06-22 |
| CN108195735B CN108195735B (en) | 2020-05-22 |
Family
ID=62573651
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711296710.3A Active CN108195735B (en) | 2017-12-08 | 2017-12-08 | Capillary pressure curve classification method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108195735B (en) |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070113619A1 (en) * | 2004-10-01 | 2007-05-24 | Newsham Kent E | Electronic humidity chamber for vapor desorption to determine high capillary pressures |
| CN102141637A (en) * | 2010-01-28 | 2011-08-03 | 中国石油天然气股份有限公司 | A Method for Continuous and Quantitative Evaluation of Reservoir Pore Structure Using Nuclear Magnetic Resonance Logging Data |
| CN103267721A (en) * | 2013-05-03 | 2013-08-28 | 中国石油天然气集团公司 | Method for evaluating water containing characteristic and occurrence state of compact sandstone storage layer aperture |
| CN103675945A (en) * | 2013-12-17 | 2014-03-26 | 中国石油天然气股份有限公司 | Method and equipment for measuring saturation of hole type reservoir |
| CN104101905A (en) * | 2013-04-11 | 2014-10-15 | 中国石油天然气集团公司 | Reservoir classification method based on rock electricity parameters |
| CN104179498A (en) * | 2013-05-22 | 2014-12-03 | 中国石油化工股份有限公司 | Zone-dividing saturation simulation method |
| CN104573198A (en) * | 2014-12-23 | 2015-04-29 | 长江大学 | Method for reconstructing digital rock core and pore network model based on random fractal theory |
| CN105190366A (en) * | 2013-03-08 | 2015-12-23 | 沙特阿拉伯石油公司 | Determining continuous capillary pressure curves for subsurface earth formations using saturation and NMR log data |
| CN105353419A (en) * | 2015-11-25 | 2016-02-24 | 中国地质大学（北京） | Nuclear magnetism capillary pressure curve construction method based on reservoir stratum classification |
| CN107153733A (en) * | 2017-05-09 | 2017-09-12 | 中国石油集团川庆钻探工程有限公司 | Capillary pressure curve fitting method for rock with complex pore structure |
-
2017
- 2017-12-08 CN CN201711296710.3A patent/CN108195735B/en active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070113619A1 (en) * | 2004-10-01 | 2007-05-24 | Newsham Kent E | Electronic humidity chamber for vapor desorption to determine high capillary pressures |
| CN102141637A (en) * | 2010-01-28 | 2011-08-03 | 中国石油天然气股份有限公司 | A Method for Continuous and Quantitative Evaluation of Reservoir Pore Structure Using Nuclear Magnetic Resonance Logging Data |
| CN105190366A (en) * | 2013-03-08 | 2015-12-23 | 沙特阿拉伯石油公司 | Determining continuous capillary pressure curves for subsurface earth formations using saturation and NMR log data |
| CN104101905A (en) * | 2013-04-11 | 2014-10-15 | 中国石油天然气集团公司 | Reservoir classification method based on rock electricity parameters |
| CN103267721A (en) * | 2013-05-03 | 2013-08-28 | 中国石油天然气集团公司 | Method for evaluating water containing characteristic and occurrence state of compact sandstone storage layer aperture |
| CN104179498A (en) * | 2013-05-22 | 2014-12-03 | 中国石油化工股份有限公司 | Zone-dividing saturation simulation method |
| CN103675945A (en) * | 2013-12-17 | 2014-03-26 | 中国石油天然气股份有限公司 | Method and equipment for measuring saturation of hole type reservoir |
| CN104573198A (en) * | 2014-12-23 | 2015-04-29 | 长江大学 | Method for reconstructing digital rock core and pore network model based on random fractal theory |
| CN105353419A (en) * | 2015-11-25 | 2016-02-24 | 中国地质大学（北京） | Nuclear magnetism capillary pressure curve construction method based on reservoir stratum classification |
| CN107153733A (en) * | 2017-05-09 | 2017-09-12 | 中国石油集团川庆钻探工程有限公司 | Capillary pressure curve fitting method for rock with complex pore structure |
Non-Patent Citations (1)
| Title |
|---|
| 汪娟 等: "伊拉克中部Ａ油田上白垩统Khasib组储层流动单元研究", 《地质科技情报》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108195735B (en) | 2020-05-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Huang et al. | Quantitative methods of evaluating image segmentation | |
| Azzalini et al. | Clustering via nonparametric density estimation: The R package pdfCluster | |
| CN103810723B (en) | Method for tracking target based on interframe constraint super-pixel coding | |
| CN102930301B (en) | Image classification method based on characteristic weight learning and nuclear sparse representation | |
| US20090097733A1 (en) | Method and apparatus for clustering and visualization of multicolor cytometry data | |
| CN102663100A (en) | Two-stage hybrid particle swarm optimization clustering method | |
| CN103984953A (en) | Cityscape image semantic segmentation method based on multi-feature fusion and Boosting decision forest | |
| CN113505225B (en) | A small-sample medical relationship classification method based on multi-layer attention mechanism | |
| CN105574547A (en) | Integrated learning method and device adapted to weight of dynamically adjustable base classifier | |
| CN101984340A (en) | Method for identifying flow type of soft grain two-phase turbulence based on artificial intelligence | |
| CN103737877B (en) | Calculating method of plastic injection moulding technological window | |
| CN104463208A (en) | Multi-view semi-supervised collaboration classification algorithm with combination of agreement and disagreement label rules | |
| CN102163239A (en) | Classifier integration method based on floating classification threshold | |
| CN107220977A (en) | The image partition method of Validity Index based on fuzzy clustering | |
| CN102842043B (en) | Particle swarm classifying method based on automatic clustering | |
| CN103246808A (en) | Low-permeability sandstone water-containing gas reservoir classification method | |
| CN102693361A (en) | Drawing method for trend curve of large dataset | |
| WO2019026523A1 (en) | Method for calculating clustering evaluation value, and method for determining number of clusters | |
| CN109902692A (en) | An Image Classification Method Based on Local Region Depth Feature Coding | |
| CN106203515A (en) | Method for selecting high-dimensional small sample data features by multi-criterion fusion | |
| CN108195735A (en) | Capillary pressure curve classification method | |
| Kane et al. | Determining the number of clusters for a k-means clustering algorithm | |
| CN105894035A (en) | SAR image classification method based on SAR-SIFT and DBN | |
| CN109977797B (en) | Optimization method of first-order target detector based on sorting loss function | |
| CN103310452B (en) | The image partition method that a kind of weight automatically selects |
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 |