CN108204936B - Characterization method of compact reservoir microscopic pore structure - Google Patents

Characterization method of compact reservoir microscopic pore structure Download PDF

Info

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
Application number
CN201611164619.1A
Other languages
Chinese (zh)
Other versions
CN108204936A (en
Inventor
李志鹏
刘显太
杨勇
郭迎春
吕广忠
陈利
卜丽侠
何伟
刘峰
孙玉花
高劲松
王玮
王瑞
张书凡
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Chemical Corp
Exploration and Development Research Institute of Sinopec Shengli Oilfield Co
Original Assignee
China Petroleum and Chemical Corp
Exploration and Development Research Institute of Sinopec Shengli Oilfield Co
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Petroleum and Chemical Corp, Exploration and Development Research Institute of Sinopec Shengli Oilfield Co filed Critical China Petroleum and Chemical Corp
Priority to CN201611164619.1A priority Critical patent/CN108204936B/en
Publication of CN108204936A publication Critical patent/CN108204936A/en
Application granted granted Critical
Publication of CN108204936B publication Critical patent/CN108204936B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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

Characterization method of compact reservoir microscopic pore structure
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:
Figure BDA0001183023720000021
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:
Figure BDA0001183023720000022
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:
Figure BDA0001183023720000041
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
Figure BDA0001183023720000042
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
Figure BDA0001183023720000051
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:
Figure BDA0001183023720000052
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:
Figure FDA0002418645310000011
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:
Figure FDA0002418645310000021
in the formula: pZThe coefficient of throat and pore configuration is mum percent; r isiIs the radius of the throat, mum; Δ SgiInterval mercury saturation,%.
CN201611164619.1A 2016-12-16 2016-12-16 Characterization method of compact reservoir microscopic pore structure Active CN108204936B (en)

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)

* Cited by examiner, † Cited by third party
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

Patent Citations (6)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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