CN108204936B - Characterization method of compact reservoir microscopic pore structure - Google Patents
Characterization method of compact reservoir microscopic pore structure Download PDFInfo
- Publication number
- CN108204936B CN108204936B CN201611164619.1A CN201611164619A CN108204936B CN 108204936 B CN108204936 B CN 108204936B CN 201611164619 A CN201611164619 A CN 201611164619A CN 108204936 B CN108204936 B CN 108204936B
- Authority
- CN
- China
- Prior art keywords
- throat
- mercury
- radius
- interval
- mercury saturation
- 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.)
- Active
Links
- 239000011148 porous material Substances 0.000 title claims abstract description 68
- 238000012512 characterization method Methods 0.000 title claims description 10
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 claims abstract description 87
- 229910052753 mercury Inorganic materials 0.000 claims abstract description 87
- 238000002474 experimental method Methods 0.000 claims abstract description 17
- 238000012545 processing Methods 0.000 claims abstract description 16
- 238000000034 method Methods 0.000 claims abstract description 15
- 238000002347 injection Methods 0.000 claims abstract description 4
- 239000007924 injection Substances 0.000 claims abstract description 4
- 239000011435 rock Substances 0.000 claims description 7
- 230000001186 cumulative effect Effects 0.000 claims description 3
- 238000011161 development Methods 0.000 abstract description 8
- 238000011156 evaluation Methods 0.000 abstract description 6
- 230000000694 effects Effects 0.000 abstract description 4
- 238000004891 communication Methods 0.000 description 8
- 238000009825 accumulation Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
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)
- Analysing Materials By The Use Of Radiation (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
The invention provides a method for characterizing a micro-pore structure of a tight reservoir, which comprises the following steps: step 1, carrying out mercury injection experiment data processing; step 2, drawing a histogram of the throat radius and the interval mercury saturation frequency by using the processed data; step 3, selecting the radius R of the peak value of the main peak according to the radius of the throat and the interval mercury saturation frequency histogramf(ii) a Step 4, carrying out mercury intrusion data interception processing; and 5, calculating the configuration coefficient of the throat and the pores, and solving the sum of the product of the radius of the throat and the mercury saturation of the corresponding interval by using mercury intrusion experiment characteristic data to serve as the configuration coefficient of the throat and the pores. The method provides a feasible method for evaluating the micro-pore structure of the compact reservoir, and has great application prospects in compact reservoir development mode selection, development effect evaluation and reservoir evaluation.
Description
Technical Field
The invention relates to the field of evaluation of tight reservoirs, in particular to a method for characterizing a micro-pore structure of a tight reservoir.
Background
The physical property of the compact reservoir is poor and is mainly influenced by the micro-pore structure, the development effect, the development mode and the like of the compact reservoir are greatly influenced by the micro-pore structure, and how to characterize the micro-pore structure characteristics of the compact reservoir is the key point of description of the compact reservoir. The medium-high permeability reservoir has thick throat, large and many pores and high communication degree between the throat and the pores, and the characterization of the micro-pore structure mainly focuses on the size of the throat. However, the compact reservoir has the basic characteristics of few pores, thin throat and poor communication degree between the throat and the pores, and the major influences on the development mode and the development effect are the radius of the throat and the communication degree between the throat and the pores. The representation of the micro-pore structure of the compact reservoir layer not only needs to describe the size and distribution of the throat, but also needs to describe the configuration relationship between the throat and the pores. The invention provides a novel method for describing a compact reservoir micro-pore structure mainly from the aspect of describing the configuration relationship between a throat and pores, and enriches the characterization parameters of the existing micro-pore structure.
Disclosure of Invention
The invention aims to provide a parameter and a method for the representation of a micro-pore structure of a tight reservoir, in particular to describe the configuration relation of a throat and pores.
The object of the invention can be achieved by the following technical measures: method for characterizing the micro-pore structure of a tight reservoir, such tight reservoirThe characterization method of the micro-pore structure comprises the following steps: step 1, carrying out mercury injection experiment data processing; step 2, drawing a histogram of the throat radius and the interval mercury saturation frequency by using the processed data; step 3, selecting the radius R of the peak value of the main peak according to the radius of the throat and the interval mercury saturation frequency histogramf(ii) a Step 4, carrying out mercury intrusion data interception processing; and 5, calculating the configuration coefficient of the throat and the pores, and solving the sum of the product of the radius of the throat and the mercury saturation of the corresponding interval by using mercury intrusion experiment characteristic data to serve as the configuration coefficient of the throat and the pores.
The object of the invention can also be achieved by the following technical measures:
in step 1, based on the original data obtained by the rock core indoor mercury intrusion experiment, obtaining the section mercury saturation corresponding to each throat radius:
ΔSgi=Sgi-Sgi-1
in the formula: Δ SgiInterval mercury saturation,%; sgiCumulative mercury saturation,%, corresponding to the i data points; sgi-1Cumulative mercury saturation,%, Sg corresponding to the i-1 data points0=0;
Thus, the interval mercury saturation frequency is calculated:
in the formula: pSgiInterval mercury saturation frequency,%; Δ SgiInterval mercury saturation,%; sg is the cumulative mercury saturation,%.
In step 2, according to the mercury intrusion experiment processing data obtained in step 1, a frequency histogram is drawn by taking the throat radius as an abscissa and taking the section mercury saturation frequency corresponding to the throat radius as an ordinate.
In step 3, according to the frequency histogram characteristics drawn in step 2, the throat radius corresponding to the main peak value is selected as the radius R of the main peak valuefThe size of the throat radius that best communicates with the pores is characterized.
In step 4, on the basis of the processing data obtained in step 1, data points with the interval mercury saturation frequency of 1% are discarded, and the data points with the interval mercury saturation frequency of less than 1% are discarded to obtain characteristic data representing the micro-pore structure characteristics, so that the influence of small probability data on the statistical result in the statistical process is eliminated, and the wrong judgment on the micro-pore structure characteristics is reduced.
In step 5, calculating the configuration relation coefficient of the micro pore structure throat and the pore of the tight reservoir by using the characteristic data obtained in step 4:
in the formula: pZThe coefficient of throat and pore configuration is mum percent; r isiIs the radius of the throat, mum; Δ SgiInterval mercury saturation,%.
The invention discloses a compact reservoir microscopic pore structure characterization method, relates to the field of reservoir microscopic pore structure characterization, and mainly aims to quantitatively describe the communication degree between a compact reservoir throat and pores. The method comprises the steps of utilizing original data of an indoor rock core mercury intrusion experiment, carrying out data processing twice, obtaining a main peak value radius (R) through a throat radius and interval mercury saturation frequency histogramf) And throat and pore space configuration factor (P)Z). Major peak radius (R)f) Mainly characterizes how large throat radius in the reservoir is best communicated with pores, and the throat and pore configuration coefficient (P)Z) And characterizing the communication degree of the throat and the pore on the whole reservoir. Obtaining interval mercury saturation frequency P corresponding to each throat radius by using original data obtained by indoor mercury intrusion experiments of rock coresSgi(ii) a According to the mercury intrusion experiment processing data, drawing a histogram with the throat radius as the abscissa and the corresponding interval mercury saturation frequency as the ordinate, and selecting the throat radius corresponding to the main peak value as the radius (R) of the main peak valuef) And taking the interval mercury saturation frequency of 1% as a cutoff value, discarding data points with the interval mercury saturation frequency of less than 1%, obtaining characteristic data representing the characteristics of the micro-pore structure, and calculating the coefficient of the relationship between the micro-pore structure throat of the compact reservoir and the pore configuration. The invention is the evaluation of the micro-pore structure of the compact reservoirProvides a feasible method, and has great application prospect in compact reservoir development mode selection, development effect evaluation and reservoir evaluation.
Drawings
FIG. 1 is a schematic diagram of a core mercury intrusion test procedure in an embodiment of the present disclosure;
FIG. 2 is a flow chart of an embodiment of a tight reservoir micro-pore structure characterization method according to an embodiment of the present invention;
FIG. 3 shows a histogram of throat radius and bin mercury saturation frequency in an embodiment of the present invention.
Detailed Description
In order to make the aforementioned and other objects, features and advantages of the present invention comprehensible, preferred embodiments accompanied with figures are described in detail below.
As shown in fig. 2, fig. 2 is a flow chart of the tight reservoir micro-pore structure characterization method of the present invention.
The pressure in the indoor rock core mercury injection experiment process rises continuously, the pressure reflects the size of the throat radius, mercury enters the rock core continuously in the pressure rising process, if the process is differentiated, the pressure is increased by delta P, the section mercury saturation degree delta Sg is increased, the pressure P corresponds to the throat radius, and the section mercury saturation degree delta Sg can be considered as the pore volume communicated with the throat radius. Based on the theory, the big throat radius of the compact reservoir and the good communication between the pore can be evaluated by using mercury intrusion experimental data, and the communication degree between the whole core throat and the pore can be evaluated.
In step 101, data processing is performed to calculate the interval mercury saturation and the interval mercury saturation frequency (interval mercury saturation/total mercury saturation) of each throat radius according to the accumulated mercury saturation obtained from the mercury intrusion experiment. In one embodiment, the section mercury saturation corresponding to each throat radius is obtained by using the original data obtained by the indoor mercury intrusion experiment of the rock core:
ΔSgi=Sgi-Sgi-1
in the formula: Δ SgiInterval mercury saturation,%; sgiAccumulation of points for i data pointsMercury saturation,%; sgi-1Cumulative mercury saturation,%, Sg corresponding to the i-1 data points0=0。
Thus, the interval mercury saturation frequency is calculated:
in the formula: pSgiInterval mercury saturation frequency,%; Δ SgiInterval mercury saturation,%; sg is the cumulative mercury saturation,%. The results are shown in Table 1.
TABLE 1 original mercury intrusion experiment data processing table
The flow proceeds to step 102.
In step 102, a histogram is plotted with the throat radius as the abscissa and the corresponding section mercury saturation frequency as the ordinate, based on the mercury intrusion experimental processing data obtained in step 101, as shown in fig. 1. The flow proceeds to step 103.
In step 103, according to the distribution form of the throat radius mercury saturation frequency histogram, the throat radius corresponding to the peak value of the main peak is selected to represent the size of the throat radius best communicated with the pores, and the throat radius is called as the throat radius of the peak value of the main peak. In one embodiment, according to the frequency histogram feature drawn in step 102, the throat radius corresponding to the main peak is selected as the main peak radius (R)f) The size of the throat radius that best communicates with the pores is characterized, as in figure 1. The flow proceeds to step 104.
In step 104, the mercury intrusion data is intercepted, the data point with the interval mercury saturation frequency of 1% is taken as a cutoff value, and the data point with the interval mercury saturation frequency of less than 1% is discarded, so as to obtain mercury intrusion experimental characteristic data. In one embodiment, the characteristic data representing the characteristic of the micro-pore structure is obtained by discarding data points having a mercury saturation frequency of less than 1% with a mercury saturation frequency of 1% as a cutoff value, as shown in table 2.
TABLE 2 mercury intrusion test characteristic data processing table and calculation results of throat and pore configuration coefficients
The flow proceeds to step 105
In step 105, calculating the configuration coefficient of the throat and the pore, and calculating the sum of the product of the radius of the throat and the mercury saturation of the interval by using the characteristic data to serve as the configuration coefficient of the throat and the pore, so as to represent the communication degree of the throat and the pore of the whole sample. Calculating the configuration relation coefficient of the micro pore structure throat and the pore of the tight reservoir by using the characteristic data obtained in the step 104:
in the formula: pZThe coefficient of throat and pore configuration is mum percent; r isiIs the radius of the throat, mum; Δ SgiInterval mercury saturation,%. The results are shown in Table 2, and the process is complete.
Claims (1)
1. The characterization method of the compact reservoir micro-pore structure is characterized by comprising the following steps:
step 1, carrying out mercury injection experiment data processing;
step 2, drawing a histogram of the throat radius and the interval mercury saturation frequency by using the processed data;
step 3, selecting the radius R of the peak value of the main peak according to the radius of the throat and the interval mercury saturation frequency histogramf;
Step 4, carrying out mercury intrusion data interception processing;
step 5, calculating the configuration coefficient of the throat and the pore, and solving the sum of the product of the radius of the throat and the mercury saturation of the corresponding interval by using mercury intrusion experiment characteristic data as the configuration coefficient of the throat and the pore;
in step 1, based on the original data obtained by the rock core indoor mercury intrusion experiment, obtaining the section mercury saturation corresponding to each throat radius:
ΔSgi=Sgi-Sgi-1
in the formula: Δ SgiInterval mercury saturation,%; sgiCumulative mercury saturation,%, corresponding to the i data points; sgi-1Cumulative mercury saturation,%, Sg corresponding to the i-1 data points0=0;
Thus, the interval mercury saturation frequency is calculated:
in the formula: pSgiInterval mercury saturation frequency,%; Δ SgiInterval mercury saturation,%; sg is cumulative mercury saturation,%;
in step 2, according to the mercury intrusion experiment processing data obtained in step 1, drawing a frequency histogram with the throat radius as the abscissa and the corresponding interval mercury saturation frequency as the ordinate;
in step 3, according to the frequency histogram characteristics drawn in step 2, the throat radius corresponding to the main peak value is selected as the radius R of the main peak valuefCharacterizing the size of the throat radius that best communicates with the pores;
in step 4, on the basis of the processing data obtained in step 1, discarding data points with the interval mercury saturation frequency of less than 1% as a cutoff value to obtain characteristic data representing the micro-pore structure characteristics, so as to eliminate the influence of small probability data on the statistical result in the statistical process and reduce the false judgment on the micro-pore structure characteristics;
in step 5, calculating the configuration relation coefficient of the micro pore structure throat and the pore of the tight reservoir by using the characteristic data obtained in step 4:
in the formula: pZThe coefficient of throat and pore configuration is mum percent; r isiIs the radius of the throat, mum; Δ SgiInterval mercury saturation,%.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611164619.1A CN108204936B (en) | 2016-12-16 | 2016-12-16 | Characterization method of compact reservoir microscopic pore structure |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611164619.1A CN108204936B (en) | 2016-12-16 | 2016-12-16 | Characterization method of compact reservoir microscopic pore structure |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108204936A CN108204936A (en) | 2018-06-26 |
CN108204936B true CN108204936B (en) | 2020-10-13 |
Family
ID=62601351
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611164619.1A Active CN108204936B (en) | 2016-12-16 | 2016-12-16 | Characterization method of compact reservoir microscopic pore structure |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108204936B (en) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103278436A (en) * | 2013-02-01 | 2013-09-04 | 西安石油大学 | Quantitative characterization method of low penetration double-medium sandstone oil reservoir microscopic aperture structure |
CN103670384A (en) * | 2013-11-18 | 2014-03-26 | 中国石油天然气集团公司 | Reservoir classification method and system |
CN103900942A (en) * | 2014-01-10 | 2014-07-02 | 中国石油大学(华东) | Method for continuously characterizing microscopic pore throat parameters of reservoir based on high-pressure mercury injection analysis |
CN104458525A (en) * | 2013-09-13 | 2015-03-25 | 中国石油天然气股份有限公司 | Acquisition and processing method for characterizing characteristics of micro-pore throat by conventional mercury intrusion curve |
CN105334149A (en) * | 2015-11-24 | 2016-02-17 | 中国石油大学(北京) | Micro-pore structure evaluation and reservoir classification method for tight reservoirs |
CN105447298A (en) * | 2014-09-29 | 2016-03-30 | 中国石油化工股份有限公司 | Reservoir stratum analysis method and application thereof |
-
2016
- 2016-12-16 CN CN201611164619.1A patent/CN108204936B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103278436A (en) * | 2013-02-01 | 2013-09-04 | 西安石油大学 | Quantitative characterization method of low penetration double-medium sandstone oil reservoir microscopic aperture structure |
CN104458525A (en) * | 2013-09-13 | 2015-03-25 | 中国石油天然气股份有限公司 | Acquisition and processing method for characterizing characteristics of micro-pore throat by conventional mercury intrusion curve |
CN103670384A (en) * | 2013-11-18 | 2014-03-26 | 中国石油天然气集团公司 | Reservoir classification method and system |
CN103900942A (en) * | 2014-01-10 | 2014-07-02 | 中国石油大学(华东) | Method for continuously characterizing microscopic pore throat parameters of reservoir based on high-pressure mercury injection analysis |
CN105447298A (en) * | 2014-09-29 | 2016-03-30 | 中国石油化工股份有限公司 | Reservoir stratum analysis method and application thereof |
CN105334149A (en) * | 2015-11-24 | 2016-02-17 | 中国石油大学(北京) | Micro-pore structure evaluation and reservoir classification method for tight reservoirs |
Non-Patent Citations (1)
Title |
---|
沧东凹陷孔二段致密砂岩储集层微观孔喉结构特征;程秋萌 等;《新疆石油地质》;20161031;第37卷(第5期);第5-6节,图4 * |
Also Published As
Publication number | Publication date |
---|---|
CN108204936A (en) | 2018-06-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109030311B (en) | Based on nuclear magnetic resonance T2Method for classifying and identifying pore structure of spectrum sensitive parameter | |
CN103278436B (en) | Quantitative characterization method of low penetration double-medium sandstone oil reservoir microscopic aperture structure | |
CN104458525B (en) | Acquisition and processing method for characterizing characteristics of micro-pore throat by conventional mercury intrusion curve | |
CN110927035A (en) | Method for calculating saturation of irreducible tight sandstone bound water | |
CN105781539A (en) | Saturability well logging calculation method of tight oil and gas reservoir | |
CN106525684B (en) | A kind of bearing calibration of the tight sand clancy behavior scale based on pore throat character | |
CN110133035B (en) | Method for quantitatively representing pore throat structure of sandstone reservoir by coupling nuclear magnetic resonance with constant-speed mercury intrusion | |
CN107329182B (en) | Method and device for determining permeability of reservoir | |
CN110702585B (en) | Rock compression coefficient calculation method | |
CN103308433A (en) | Method for analyzing and evaluating tight sandstone reservoir diagenetic facies based on porosity evolution | |
CN106897531B (en) | Quantitative evaluation method for permeability of low-permeability limestone reservoir | |
CN106777515B (en) | Method for analyzing production energy of tight gas well based on rock core experiment data | |
CN108204936B (en) | Characterization method of compact reservoir microscopic pore structure | |
CN103778328A (en) | Reservoir sensitive pore throat extracting method based on standard-deviation analysis | |
CN115078210B (en) | Shale pore structure testing method | |
CN109598245B (en) | Edible oil transverse relaxation attenuation curve signal feature extraction method based on 1D-CNN | |
CN105334293A (en) | Flow state identification method and device for dense gas in multiple media | |
CN108805158B (en) | Compact oil reservoir lithogenous phase dividing method | |
CN113624847A (en) | Method for establishing prediction model of shale hydration damage coefficient and prediction method | |
CN108982326A (en) | A kind of simplification prediction technique of the saturation soil body osmotic coefficient based on fractal theory | |
CN110632677B (en) | Carbonate rock geothermal reservoir classification method | |
CN109211756B (en) | Method for predicting deformation soil air intake value based on VG model | |
CN112284991A (en) | Fractal theory-based porous medium permeability prediction method | |
CN113253354A (en) | Nuclear magnetic T2Method, device, medium and electronic device for determining spectrum cutoff value | |
CN112014906B (en) | Compact reservoir evaluation 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 |