CN111624146A - Method for quantitatively characterizing full-scale pore size distribution characteristics of compact reservoir - Google Patents
Method for quantitatively characterizing full-scale pore size distribution characteristics of compact reservoir Download PDFInfo
- Publication number
- CN111624146A CN111624146A CN201910151245.7A CN201910151245A CN111624146A CN 111624146 A CN111624146 A CN 111624146A CN 201910151245 A CN201910151245 A CN 201910151245A CN 111624146 A CN111624146 A CN 111624146A
- Authority
- CN
- China
- Prior art keywords
- pore
- reservoir
- pore size
- target
- relaxation time
- 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.)
- Pending
Links
- 239000011148 porous material Substances 0.000 title claims abstract description 209
- 238000009826 distribution Methods 0.000 title claims abstract description 102
- 238000000034 method Methods 0.000 title claims abstract description 44
- 230000001186 cumulative effect Effects 0.000 claims abstract description 39
- 238000012512 characterization method Methods 0.000 claims abstract description 25
- 238000005516 engineering process Methods 0.000 claims abstract description 14
- 238000012417 linear regression Methods 0.000 claims abstract description 8
- 238000001225 nuclear magnetic resonance method Methods 0.000 claims abstract description 4
- 238000011161 development Methods 0.000 claims description 12
- 230000018109 developmental process Effects 0.000 claims description 12
- 238000003384 imaging method Methods 0.000 claims description 8
- 230000011218 segmentation Effects 0.000 claims description 8
- 230000000007 visual effect Effects 0.000 claims description 8
- 238000006243 chemical reaction Methods 0.000 abstract description 5
- 238000005481 NMR spectroscopy Methods 0.000 description 19
- 239000012530 fluid Substances 0.000 description 9
- 238000012360 testing method Methods 0.000 description 8
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 description 6
- 229910052753 mercury Inorganic materials 0.000 description 6
- 238000001228 spectrum Methods 0.000 description 6
- 239000007789 gas Substances 0.000 description 5
- 239000011435 rock Substances 0.000 description 5
- 238000009792 diffusion process Methods 0.000 description 4
- 229910001919 chlorite Inorganic materials 0.000 description 3
- 229910052619 chlorite group Inorganic materials 0.000 description 3
- QBWCMBCROVPCKQ-UHFFFAOYSA-N chlorous acid Chemical compound OCl=O QBWCMBCROVPCKQ-UHFFFAOYSA-N 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000002347 injection Methods 0.000 description 3
- 239000007924 injection Substances 0.000 description 3
- 229910052500 inorganic mineral Inorganic materials 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 239000011707 mineral Substances 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 229920006395 saturated elastomer Polymers 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000010835 comparative analysis Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000005755 formation reaction Methods 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 241000950275 Epilachna deltoides Species 0.000 description 1
- 235000019738 Limestone Nutrition 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- XKRFYHLGVUSROY-UHFFFAOYSA-N argon Substances [Ar] XKRFYHLGVUSROY-UHFFFAOYSA-N 0.000 description 1
- 229910052786 argon Inorganic materials 0.000 description 1
- -1 argon ion Chemical class 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 239000002734 clay mineral Substances 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000003628 erosive effect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000000445 field-emission scanning electron microscopy Methods 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 125000004435 hydrogen atom Chemical class [H]* 0.000 description 1
- 239000006028 limestone Substances 0.000 description 1
- 238000013332 literature search Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000001000 micrograph Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005498 polishing Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000004626 scanning electron microscopy Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000001179 sorption measurement Methods 0.000 description 1
- 238000010998 test 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
-
- 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
- G01N2015/0813—Measuring intrusion, e.g. of mercury
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)
- Geophysics And Detection Of Objects (AREA)
Abstract
The application provides a method for quantitatively characterizing full-scale pore size distribution characteristics of a compact reservoir, which comprises the following steps: s1, obtaining a relaxation time cumulative distribution curve of the compact reservoir sample by a nuclear magnetic resonance method; s2, acquiring a pore diameter cumulative distribution curve of the reservoir sample by a field emission scanning electron microscope full-scale pore diameter quantitative characterization technology; and S3, acquiring the target function relation between the relaxation time and the pore size distribution of the reservoir sample by utilizing the relaxation time cumulative distribution curve and the pore size cumulative distribution curve and combining a linear regression method. By the method, the conversion from the relaxation time distribution of compact reservoir gaps to the pore size distribution can be realized, and the full-scale characterization of the compact reservoir pore size is quantitatively characterized.
Description
Technical Field
The invention relates to the technical field of oil and gas exploration, in particular to a method for quantitatively characterizing full-scale pore size distribution characteristics of a compact reservoir.
Background
Modern unconventional oil and gas exploration practices reveal that oil and gas accumulation and distribution in tight sandstone reservoirs are very different from oil and gas occurrence in conventional sandstone reservoirs. In tight sandstone reservoirs, hydrocarbons are mainly collected and distributed in micro-nano pore throat network systems therein. Compared with a millimeter-micron pore throat network system of a conventional sandstone reservoir, the micron-nanometer pore throat network system has the advantages of more complex structure, smaller pore diameter and stronger heterogeneity, and particularly, the full-scale pore diameter distribution characteristics are difficult to describe and need to be realized by establishing a targeted method system.
The existing compact sandstone reservoir pore diameter full-scale quantitative characterization method mostly adopts a method of combining a nuclear magnetic resonance technology and a high-pressure mercury injection technology. In the method, the main measurement range of the nuclear magnetic resonance technology is 5nm to 100 mu m, the detection range of the high-pressure mercury injection technology is 10nm to 100 mu m, and the method combines the throat distribution result obtained in the high-pressure mercury injection experiment with the T of the nuclear magnetic resonance2And carrying out comparative analysis on the relaxation time result, and converting the relaxation time of nuclear magnetic resonance into the pore size distribution of the compact reservoir through formula deduction so as to realize the quantitative characterization of the pore size. However, the method has a more obvious problem that the result tested by the high-pressure mercury intrusion test is the throat distribution characteristic of the compact reservoir, the result tested by the nuclear magnetic resonance test is the pore development characteristic of the compact reservoir, the two are two main components in the microstructure of the reservoir respectively, and the two should not be confused and mutually converted.
Through patent and literature search, the currently disclosed aperture full-scale characterization method is greatly different from the technology provided by the invention, and the technical categories related to the aperture full-scale characterization method provided by the invention and the differences of the technical categories and the technology provided by the invention are described as follows.
The shale pore size distribution test method with the publication number of CN105486621B is mainly characterized in that a gas adsorption test is adopted to represent the pore size distribution characteristics of 2nm-200nm, a high-pressure mercury intrusion test result is adopted, and a Washburn equation is utilized to calculate the residual pore size distribution characteristics. The method is a pore size characterization method designed for shale reservoirs, and various test results are spliced. Shale reservoir and methodThe pore size distribution of the tight sandstone reservoir is obviously different, the invention aims at the pore size characterization of the tight sandstone, adopts the conversion based on the relaxation time of the nuclear magnetic resonance T2, and directly converts T2The spectrum is converted into a pore size distribution curve, and the problem of splicing does not exist. Literature nuclear magnetic resonance T2The research of the spectrum conversion pore radius distribution method is mainly realized by comparing nuclear magnetic resonance T2The pore radius is converted by the spectrum and the pore throat distribution curve measured by a high-pressure mercury intrusion method, and the pore radius is mixed with the throat radius. According to the invention, the nuclear magnetic resonance and scanning electron microscope image processing are both characterized by the pore radius, so that the problem that the pore and the throat do not correspond is effectively avoided.
Therefore, a calculation method for the quantitative distribution characterization of the pore diameter under the full scale of the tight sandstone reservoir is needed to be provided, so that the problem that the distribution of the pore diameter in the tight reservoir is not clear under the full scale due to strong heterogeneity is solved, and the understanding of the quantitative characterization of the pore diameter in the tight reservoir can be solved.
Disclosure of Invention
Aiming at the problems in the prior art, the nuclear magnetic resonance technology and the field emission scanning electron microscope are discovered to have relatively consistent detection ranges through comparative analysis, the main measurement ranges are all about 5nm to 100 mu m, and the pore spaces in all scales of the compact sandstone reservoir can be basically covered, so that the nuclear magnetic resonance technology and the field emission scanning electron microscope are in contact, and the quantitative conversion representation of the relaxation time distribution to the pore size distribution of the nuclear magnetic resonance is realized.
The method for quantitatively characterizing the full-scale pore size distribution characteristics of the compact reservoir provided by the invention comprises the following steps: s1, obtaining a relaxation time cumulative distribution curve of the compact reservoir sample by a nuclear magnetic resonance method; s2, acquiring a pore diameter cumulative distribution curve of the reservoir sample by a field emission scanning electron microscope full-scale pore diameter quantitative characterization technology; and S3, acquiring the target function relation between the relaxation time and the pore size distribution of the reservoir sample by utilizing the relaxation time cumulative distribution curve and the pore size cumulative distribution curve and combining a linear regression method. By the method, the conversion from the relaxation time distribution of compact reservoir gaps to the pore size distribution can be realized, and the full-scale characterization of the compact reservoir pore size is quantitatively characterized.
In a preferred embodiment, step S2 includes: s21, obtaining a target continuous image representing full-scale reservoir pore development characteristics by a plurality of sub-images through an image splicing method; s22, separating the target continuous image by a threshold segmentation method to obtain pore binary images of different types of pores; s23, quantitatively counting the reservoir pore parameters in the pore binary image to obtain the pore diameter cumulative distribution curve. Through the implementation mode, the distribution characteristics of different types of pores of the compact reservoir under the full pore size can be identified, so that the understanding of the pore-throat structure development in the compact reservoir is increased, and the area percentages and the pore size distribution of different sizes and different types of pores are extracted from the two-dimensional image.
In a preferred embodiment, step S21 includes: s11, determining a target view capable of effectively representing full-scale reservoir pore development characteristics through a field emission scanning electron microscope; s12, dividing the target view into a plurality of sub-view grids; s13, respectively carrying out high-resolution imaging on each sub-view grid in the sub-view grids to obtain a plurality of sub-images; and S14, splicing the multiple sub-images in sequence under the premise of ensuring the integrity to obtain the target continuous image.
In a preferred embodiment, step S22 includes: s221, determining a gray threshold; s222, segmenting pixels of which the gray values are greater than or equal to the gray threshold value in the target continuous image into backgrounds; s223, dividing pixels with the gray value smaller than the gray threshold value in the target continuous image into pores to obtain a pore binary image; and S224, repeating the steps S221 to S223 to obtain the pore binary images of different types of pores.
In a preferred embodiment, in step S3, the objective function relationship is expressed by the following equation:
wherein r iscPore size, T, of the reservoir sample2For the relaxation time of the reservoir sample, n and C are the power exponent and constant, respectively, to be determined.
In a preferred embodiment, step S3 includes: s31, determining a relaxation time value and an aperture value corresponding to the same ordinate value in the relaxation time cumulative distribution curve and the aperture cumulative distribution curve; s32, determining the relaxation time value and the aperture value as an abscissa and an ordinate respectively to obtain a target value pair; s33, repeating the steps S31-S32 to obtain at least one target value pair; and S34, substituting the at least one target value pair into the target function relational expression, and solving n and C by using a linear regression method.
In a preferred embodiment, the gray threshold is 30.
In a preferred embodiment, the reservoir pore parameters include pore size, number of pores, and areal porosity.
In a preferred embodiment, the target field of view has a dimension greater than 10 times the reservoir pore size.
In a preferred embodiment, the size of the target field of view is greater than 1 mm.
In a preferred embodiment, the high resolution is 8.0 × 105ppi。
The method for quantitatively characterizing the full-scale pore size distribution of the compact reservoir provided by the invention can convert the nuclear magnetic resonance relaxation time distribution of the pores of the reservoir into the pore size distribution, realize the quantitative characterization of the pore size distribution of the compact reservoir under the full scale, solve the problem of ambiguous distribution of the pore size in the compact reservoir under the full scale due to strong heterogeneity, and also solve the recognition of the quantitative characterization of the pore size in the compact reservoir.
The features mentioned above can be combined in various suitable ways or replaced by equivalent features as long as the object of the invention is achieved.
Drawings
The invention will be described in more detail hereinafter on the basis of embodiments and with reference to the accompanying drawings. Wherein:
FIG. 1 shows a schematic flow diagram of a method for quantitative characterization of tight reservoir full-scale pore size distribution according to an embodiment of the present invention;
FIG. 2 shows a schematic flow diagram of a full-scale aperture quantitative characterization method according to an embodiment of the present invention;
FIG. 3 shows a schematic flow chart of a method of acquiring successive images of a target according to an embodiment of the invention;
FIG. 4 shows a schematic diagram of a plurality of sub-field meshing according to an embodiment of the invention;
FIGS. 5 and 6 show a grayscale image before segmentation of chlorite intergranular pores and a binarized image after segmentation, respectively, according to an embodiment of the present invention;
FIG. 7 shows an image of a two-dimensional target sequential image reflecting a full-scale aperture of a reservoir in accordance with an embodiment of the invention;
FIG. 8 shows a comparison of compact reservoir pore relaxation time cumulative distribution curves and pore size cumulative distribution curves in accordance with an embodiment of the present invention.
In the drawings, like parts are provided with like reference numerals. The drawings are not to scale.
Detailed Description
The invention will be further explained with reference to the drawings.
Fig. 1 is a schematic flow chart 100 of a method for quantitatively characterizing the full-scale pore size of a tight reservoir provided by the present invention. As shown in fig. 1, the method 100 includes:
s110, acquiring a relaxation time cumulative distribution curve of the compact reservoir sample by a nuclear magnetic resonance method;
s120, acquiring a pore diameter cumulative distribution curve of the reservoir sample by a field emission scanning electron microscope full-scale pore diameter quantitative characterization technology;
and S130, acquiring a target function relation between the relaxation time and the pore size distribution of the reservoir sample by utilizing the relaxation time cumulative distribution curve and the pore size cumulative distribution curve and combining a linear regression method.
Specifically, in S110, according to the selected target layer compact sandstone sample, performing nuclear magnetic resonance test on T2Measuring relaxation time, wherein the magnitude of the relaxation time depends on the strength of the acting force of the fluid molecules on the solid surface, and accumulating the relaxation time as an abscissaAnd taking the distribution value as a vertical coordinate to obtain a full-scale pore relaxation time cumulative distribution curve of the compact reservoir sample.
For a pure material sample (such as pure water) in which the surrounding environment and interaction of each hydrogen nucleus are the same, the relaxation of the atomic nuclei in the individual pores can be expressed by a relaxation time through a diffusion surface relaxation model of nuclear magnetic resonance, when T is2Can be expressed as:
the first term on the right side of equation (1) is the bulk relaxation term, the magnitude of which depends on the properties of the saturated fluid. Volume relaxation T2BIt is meant that when a fluid is present in the larger pores (believed to be independent of the pore space), the fluid itself undergoes a decay, also referred to as free relaxation. For a porous medium sample such as reservoir rock, which has a very complex mineral composition and pore structure, a fluid is present in the porous medium, and is divided and surrounded by a plurality of interfaces, the shape and size of the pores are not uniform, and the pore space is not restricted, so that the volume relaxation has no obvious correlation with the pore surface, but is related to the whole test system, such as the temperature, rock wettability, fluid viscosity and the like under the environment where the test system is located, wherein the most important parameters are the properties of the fluid in the pores.
In nmr applications, the relaxation of the fluid is generally negligible, since the relaxation strength of the rock surface is much greater than the free relaxation strength of the fluid.
The third term on the right side of the formula (1) is a diffusion relaxation term, and the diffusion relaxation term can be removed according to an experimental technology of nuclear magnetic resonance diffusion removal measurement. After removing the first term and the third term on the right side of the equation, equation (1) can be simplified as follows:
in the formula, ρ2The surface relaxation strength is controlled by mineral compositionAnd pore surface properties; S/V is the specific surface of an individual pore, the size of which is inversely proportional to the pore radius. It can be seen that when the pore structure is more complex, the pores are smaller and the specific surface area is larger, the influence of the pore surface interaction is stronger, and T2The shorter the time.
As can be seen from formula (2), the observed relaxation time T2Associated with the specific surface of the porous medium. If the pore structure is simplified into a spherical or columnar shape, the relationship between the specific surface area of the pores and the pore diameter is changed toThen there are:
in the formula, FsIs the shape factor of the pore (for spherical pores, F)s3; for columnar voids, Fs2), dimensionless; r iscPore radius, μm.
The actual formation has complicated pore structure, T2The distribution is in a power function relationship with the pore radius:
in the formula: n is a power exponent and is dimensionless.
Because the rho cannot be actually measured by the current equipment method2And FsTherefore, it is impossible to use the formula (4) for converting the nuclear magnetic resonance T2The distribution is transformed with the pore radius distribution curve.
When the core is 100% saturated water, the nuclear magnetic resonance T of the core is utilized2The spectra allow the evaluation of the pore size and its corresponding pore volume distribution. For the pore size distribution characterization of the tight reservoir, the pore size distribution characterization can be obtained by using a field emission scanning electron microscopy, namely, the step S120 provided by the invention: acquiring a pore size cumulative distribution curve of the reservoir sample by using a full-scale quantitative characterization technology of a field emission scanning electron microscope, as shown in fig. 2, the step S120 may specifically include the following steps:
s121, obtaining a target continuous image representing full-scale reservoir pore development characteristics by using a plurality of sub-images through an image splicing method;
s122, separating the target continuous image by adopting a threshold segmentation method to obtain pore binary images of different types of pores;
s123, quantitatively counting reservoir pore parameters in the pore binary image to obtain the pore diameter cumulative distribution curve.
Briefly, pores and throats as low as 5-10 nm in compact sandstone are obtained, and image splicing continuous imaging can be realized through software, namely, a plurality of images with overlapped parts are spliced into a large-scale seamless high-resolution image. After an image capable of effectively standardizing the full pore development characteristics of the reservoir is obtained, the area percentages and pore size distribution of pores of different sizes and different types are extracted from the two-dimensional image.
Herein, the different types of pores include intergranular pores, erosion pores, clay mineral intergranular pores, and the like.
Specifically, as shown in fig. 3, in step S121, a target continuous image capable of characterizing pore development of a pore size full-scale reservoir is obtained, wherein:
in S1211, a target field of view capable of effectively characterizing full-scale reservoir pore development features is determined by scanning electron microscopy. Specifically, under a FIB-SEM double-beam scanning electron microscope, finding a visual field of the tight sandstone sample, wherein as mentioned above, the sizes of pores and throats in the target visual field can be as low as 5-10 nm, and the target visual field of the tight sandstone should be larger than 1mm, so that the tight sandstone sample under the visual field comprises pores of various types and sizes, and the tight sandstone sample can be truly and completely characterized.
As shown in fig. 4, in S1212, the target view selected in S1211 is divided into a plurality of sub-view grids under the FIB-SEM. For example, the number of the sub-view grids is 25 × 25, wherein the sizes of the sub-view grids may be the same or different, and the invention is not limited herein.
Next, each of the plurality of sub-field meshes is subjected to high resolution imaging in S1213. In operation, the FIB-SEM may be focused on a target sub-field grid and then set to a high resolution at which the target sub-field grid is imaged so that the smaller pores therein can also be delineated. It should be understood that during the imaging process, the range of imaging should be larger than the range of the target sub-field grid to enable alignment for complete imaging, avoid omission, and ensure accuracy of the result. And (4) carrying out high-resolution imaging (with consistent resolution) on each sub-view grid one by one according to the steps to obtain a plurality of sub-images reflecting the pore characteristics of the sub-view grids.
It can be seen that the multiple sub-images are small-scale and high-resolution images, and the edges of two adjacent sub-images have an overlapping region. Optionally, the multiple sub-images may also be acquired by other means, such as images obtained at different times, at different viewing angles, or by different sensors.
Preferably, high resolution refers to 8.0 × 105ppi。
In S1214, the multiple sub-images are sequentially re-stitched by the MAPS1.0 software to obtain a target continuous image, that is, multiple images with overlapping regions are stitched into a large seamless high-resolution image, the image can reflect the full-scale pores of the tight sandstone sample at high resolution, and the two-dimensional target continuous image is a grayscale image, and the pores can be identified by naked eyes. It should be understood here that, in the stitching, it should be premised on ensuring the integrity, that is, it needs to be relied on according to the arrangement order of the multiple sub-views in the original target view, so as to ensure the integrity of the finally obtained two-dimensional target continuous image.
After the image capable of effectively representing the full-scale pore development characteristics of the reservoir is obtained, in order to identify and represent different types of pores, a reasonable segmentation algorithm is selected to distinguish the pores from the skeleton, namely, the area percentages and pore size distribution of the pores with different scales and different types are extracted from the gray level image. In the embodiment of the present invention, identification and characterization are performed in S122, that is, a pore binary image of different types of pores is obtained by separating from the target continuous image by using a threshold segmentation method. Specifically, step S120 includes the steps of:
and S1221, determining a gray threshold value.
S1222, dividing the pixels with the gray value larger than or equal to the gray threshold value in the target continuous image into the background;
s1223, dividing pixels with the gray values smaller than the gray threshold value in the target continuous image into pores to obtain a pore binary image;
it is assumed that the gray values of adjacent pixels within the object (aperture) or background (skeleton) in the gray-scale image are similar, while the gray values of the pixels between the object and the background are different. Based onProfessional image processing software, taking a chlorite intercrystalline pore image (fig. 5) as an example, finds that the gray values of the pore regions are mainly distributed in [0, 60 ] through manual adjustment of a segmentation threshold value and visual observation of the pore extraction process]That is, when the division threshold is about 30, the division of substantially all pores can be realized, and the result is shown in fig. 6.
And S1224, repeating the steps S1221 to S1223 to obtain the pore binary images of different types of pores.
And S123, after obtaining the pore binary images of different types of pores of the compact sandstone sample through S122, characterizing the pore parameters. Specifically, based on the binary image, parameters such as the size, the number, the surface porosity and the like of different types of pores can be quantitatively counted by using software, so that the close combination of qualitative pore description and quantitative pore research is realized. And taking the pore radius as a horizontal coordinate and the cumulative distribution value as a vertical coordinate to obtain the full-scale pore radius cumulative distribution curve of the compact reservoir sample.
In the embodiment of the invention, the pore space of the full-scale reservoir is quantitatively represented by taking the tight sandstone as an example, however, the method is not only suitable for the tight sandstone, but also suitable for other rock formations, such as shale, tight limestone and other unconventional reservoirs.
The inventors have made attempts in the practice of a dense oil exploration work using the above method. Taking an Erdos basin H221 well 36# sample as an example, the distribution of chlorite mineral intercrystalline pores is considered, a large-area argon ion polishing sample is selected, an analysis area is divided into 625 grids with 25 parts of transverse and longitudinal directions, the horizontal range of each grid is 50 microns, the minimum pore can be identified to be 58nm, all effective pores can be basically observed, the horizontal range of the whole area is 1.19mm (as shown in figure 6), the lower limit of the large visual field size of the tight sandstone is also exceeded, and the development and distribution characteristics of reservoir pores can be effectively represented.
As shown in table 1 below, a total of 298703 effective pores with an areal porosity of 4.71% were identified by the method of the present invention, with the smallest pores having a pore diameter of 58nm and the largest pores having a pore diameter of 15364 nm. The distribution of the pore diameter is in a multi-modal distribution, and the pore distribution is the largest at 900-1000 nm.
TABLE 1 statistical table for pore characterization of samples H221-36
In S130, the relaxation time cumulative distribution curve obtained in S110 and the pore size cumulative distribution curve obtained in S120 are used in combination with a linear regression method to obtain an objective function relationship between the relaxation time and the pore size distribution of the reservoir sample, so that the nmr cumulative distribution curve and the pore radius cumulative distribution curve can be converted into each other.
the equation (5) is the objective function relationship between relaxation time and pore radius, where T2As pore relaxation time, rcN and C are the power index and constant, respectively, to be determined for the pore radius.
It can be seen that there is a clear similarity in morphology between the pore radius cumulative distribution curve and the relaxation time cumulative distribution curve, indicating that they are comparable. Specifically, step S130 may include the steps of:
s131, determining a relaxation time value T corresponding to the same ordinate value in the relaxation time cumulative distribution curve and the pore size cumulative distribution curve20And the pore size rc0;
S132, converting the relaxation time value T20And the aperture value rc0Respectively determined as abscissa and ordinate to obtain target value pair (T)20,rc0);
S133, repeating the steps S131 to S132 to obtain at least one target value pair (T)2,rc);
S134, the at least one target value pair (T) is processed2,rc) Substituting the obtained product into the objective function relational expression, and solving n and C by using a linear regression method.
Selecting the relaxation time value with the same accumulated content and the corresponding aperture value through the pore radius accumulated distribution curve and the relaxation time accumulated distribution curve, and substituting the values into a formula (5) to obtain the values of C and n, thereby realizing the T of the rock with 100 percent of saturated water2And directly converting the distribution curve into a pore size distribution curve to obtain the quantitative pore size representation of the compact sandstone sample.
Taking the E.deltoid basin H221 well 35# sample as an example, let T2Converting the spectrum and the pore size distribution into cumulative distribution curves, and selecting a series of T with the same cumulative content2The values and corresponding aperture values (FIG. 8) are substituted into equation (5) and C and n are iteratively solved by regression to achieve T2The spectra are converted to pore size distributions. On the basis of the method, the average value of C of a compact reservoir in the 7 sections of the research area is 181, the average value of n is 0.8248, and then the objective function relation between the pore radius and the relaxation time of the sample is as follows:
rc=181T2 1.2124
the method for quantitatively characterizing the full-scale pore size distribution of the compact reservoir provided by the invention can convert the nuclear magnetic resonance relaxation time distribution of the pores of the reservoir into the pore size distribution, realize the quantitative characterization of the pore size distribution of the compact reservoir under the full scale, solve the problem of ambiguous distribution of the pore size in the compact reservoir under the full scale due to strong heterogeneity, and also solve the recognition of the quantitative characterization of the pore size in the compact reservoir.
In the description of the present invention, it is to be understood that the terms "upper", "lower", "bottom", "top", "front", "rear", "inner", "outer", "left", "right", and the like, indicate orientations or positional relationships based on the orientations or positional relationships shown in the drawings, are only for convenience in describing the present invention and simplifying the description, and do not indicate or imply that the device or element being referred to must have a particular orientation, be constructed in a particular orientation, and be operated, and thus, should not be construed as limiting the present invention.
Although the invention herein has been described with reference to particular embodiments, it is to be understood that these embodiments are merely illustrative of the principles and applications of the present invention. It is therefore to be understood that numerous modifications may be made to the illustrative embodiments and that other arrangements may be devised without departing from the spirit and scope of the present invention as defined by the appended claims. It should be understood that various dependent claims and features described herein may be combined in ways other than those described in the original claims. It is also to be understood that features described in connection with a single embodiment may be used in other such embodiments.
Claims (10)
1. A method for quantitatively characterizing the full-scale pore size distribution of a tight reservoir is characterized by comprising the following steps:
s1, obtaining a relaxation time cumulative distribution curve of the compact reservoir sample by a nuclear magnetic resonance method;
s2, acquiring a pore diameter cumulative distribution curve of the reservoir sample by a field emission scanning electron microscope full-scale pore diameter quantitative characterization technology;
and S3, acquiring the target function relation between the relaxation time and the pore size distribution of the reservoir sample by utilizing the relaxation time cumulative distribution curve and the pore size cumulative distribution curve and combining a linear regression method.
2. The method according to claim 1, wherein step S2 includes:
s21, obtaining a target continuous image representing full-scale reservoir pore development characteristics by a plurality of sub-images through an image splicing method;
s22, separating the target continuous image by adopting a threshold segmentation method to obtain pore binary images of different types of pores;
s23, quantitatively counting reservoir pore parameters in the pore binary image to obtain the pore diameter cumulative distribution curve.
3. The method according to claim 2, wherein step S21 includes:
s11, determining a target view capable of effectively representing full-scale reservoir pore development characteristics through a field emission scanning electron microscope;
s12, dividing the target visual area into a plurality of sub visual area grids;
s13, respectively carrying out high-resolution imaging on each sub-view grid in the sub-view grids to obtain a plurality of sub-images;
and S14, splicing the sub-images in sequence on the premise of ensuring the completeness to obtain the target continuous image.
4. The method according to claim 2, wherein step S22 includes:
s221, determining a gray threshold;
s222, segmenting pixels of which the gray values are larger than or equal to the gray threshold value in the target continuous image into backgrounds;
s223, dividing pixels with the gray values smaller than the gray threshold value in the target continuous image into pores to obtain a pore binary image;
s224, repeating the steps S21-S23 to obtain the pore binary images of different types of pores.
5. The method according to claim 1, wherein in step S3, the objective function relationship is expressed by the following equation:
wherein r iscIs the pore size, T, of the reservoir sample2N and C are the power exponent and constant, respectively, to be determined for the relaxation time of the reservoir sample.
6. The method of claim 5, wherein step S3 includes:
s31, determining a relaxation time value and an aperture value corresponding to the same ordinate value in the relaxation time cumulative distribution curve and the aperture cumulative distribution curve;
s32, determining the relaxation time value and the aperture value as an abscissa and an ordinate respectively to obtain a target value pair;
s33, repeating the steps S31-S32 to obtain at least one target value pair;
and S34, substituting the at least one target value pair into the target function relational expression, and solving n and C by using a linear regression method.
7. The method of claim 4, wherein the grayscale threshold is 30.
8. The method of claim 2, wherein the reservoir pore parameters comprise pore size, pore number, and areal porosity.
9. The method of claim 3, wherein the target field of view has a size greater than 10 times a reservoir pore size.
10. The method of claim 9, wherein the size of the target field of view is greater than 1 mm.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910151245.7A CN111624146A (en) | 2019-02-28 | 2019-02-28 | Method for quantitatively characterizing full-scale pore size distribution characteristics of compact reservoir |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910151245.7A CN111624146A (en) | 2019-02-28 | 2019-02-28 | Method for quantitatively characterizing full-scale pore size distribution characteristics of compact reservoir |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111624146A true CN111624146A (en) | 2020-09-04 |
Family
ID=72271626
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910151245.7A Pending CN111624146A (en) | 2019-02-28 | 2019-02-28 | Method for quantitatively characterizing full-scale pore size distribution characteristics of compact reservoir |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111624146A (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112284999A (en) * | 2020-10-26 | 2021-01-29 | 中国石油大学(华东) | Sandstone pore size distribution determination method and application thereof |
CN114428043A (en) * | 2020-09-18 | 2022-05-03 | 中国石油化工股份有限公司 | Porous medium pore size distribution characterization method and electronic equipment |
CN114577673A (en) * | 2022-01-18 | 2022-06-03 | 重庆大学 | NMR characterization method for pore wettability of equal-particle-size coal powder |
CN115078210A (en) * | 2022-05-05 | 2022-09-20 | 中国科学院地质与地球物理研究所 | Shale pore structure testing method |
CN116840124A (en) * | 2023-06-21 | 2023-10-03 | 西南石油大学 | Method for quantitatively representing shale reservoir core hole seam configuration |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105424576A (en) * | 2015-11-12 | 2016-03-23 | 中国科学院地质与地球物理研究所 | Shale porosity distribution uniformity evaluation method |
CN106053315A (en) * | 2016-06-06 | 2016-10-26 | 中国石油大学(华东) | Method for classifying bioclastic limestone reservoir pore structures |
CN108169099A (en) * | 2018-01-17 | 2018-06-15 | 西南石油大学 | A kind of shale gas RESERVOIR PORE STRUCTURE quantitative calculation method based on nuclear magnetic resonance |
CN109283114A (en) * | 2018-11-07 | 2019-01-29 | 西安石油大学 | A kind of quantitatively characterizing method of tight sand oil reservoir full size pore throat |
-
2019
- 2019-02-28 CN CN201910151245.7A patent/CN111624146A/en active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105424576A (en) * | 2015-11-12 | 2016-03-23 | 中国科学院地质与地球物理研究所 | Shale porosity distribution uniformity evaluation method |
CN106053315A (en) * | 2016-06-06 | 2016-10-26 | 中国石油大学(华东) | Method for classifying bioclastic limestone reservoir pore structures |
CN108169099A (en) * | 2018-01-17 | 2018-06-15 | 西南石油大学 | A kind of shale gas RESERVOIR PORE STRUCTURE quantitative calculation method based on nuclear magnetic resonance |
CN109283114A (en) * | 2018-11-07 | 2019-01-29 | 西安石油大学 | A kind of quantitatively characterizing method of tight sand oil reservoir full size pore throat |
Non-Patent Citations (8)
Title |
---|
HUCHENG DENG 等: "Improved pore-structure characterization in shale formations with FESEM technique", 《JOURNAL OF NATURAL GAS SCIENCE AND ENGINEERING》, vol. 35, 25 August 2016 (2016-08-25), pages 309 - 319, XP029777159, DOI: 10.1016/j.jngse.2016.08.063 * |
孙小龙 等: "基于核磁共振标定的高压压汞孔喉分布定量评价方法", 《岩矿测试》 * |
孙小龙 等: "基于核磁共振标定的高压压汞孔喉分布定量评价方法", 《岩矿测试》, vol. 36, no. 6, 30 November 2017 (2017-11-30), pages 601 - 607 * |
李海波 等: "核磁共振T2 谱换算孔隙半径分布方法研究", 《波谱学杂志》 * |
李海波 等: "核磁共振T2 谱换算孔隙半径分布方法研究", 《波谱学杂志》, vol. 25, no. 2, 30 June 2008 (2008-06-30), pages 273 - 280 * |
李爱芬 等: "核磁共振研究致密砂岩孔隙结构的方法及应用", 《中国石油大学学报(自然科学版)》 * |
李爱芬 等: "核磁共振研究致密砂岩孔隙结构的方法及应用", 《中国石油大学学报(自然科学版)》, vol. 39, no. 6, 31 December 2015 (2015-12-31), pages 93 - 95 * |
赵婕: "图像特征提取与语义分析", vol. 1, 重庆大学出版社, pages: 80 - 82 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114428043A (en) * | 2020-09-18 | 2022-05-03 | 中国石油化工股份有限公司 | Porous medium pore size distribution characterization method and electronic equipment |
CN112284999A (en) * | 2020-10-26 | 2021-01-29 | 中国石油大学(华东) | Sandstone pore size distribution determination method and application thereof |
CN112284999B (en) * | 2020-10-26 | 2022-06-17 | 中国石油大学(华东) | Sandstone pore size distribution determination method and application thereof |
CN114577673A (en) * | 2022-01-18 | 2022-06-03 | 重庆大学 | NMR characterization method for pore wettability of equal-particle-size coal powder |
CN114577673B (en) * | 2022-01-18 | 2024-02-02 | 重庆大学 | Pore wettability NMR characterization method of equal-particle-size pulverized coal |
CN115078210A (en) * | 2022-05-05 | 2022-09-20 | 中国科学院地质与地球物理研究所 | Shale pore structure testing method |
CN116840124A (en) * | 2023-06-21 | 2023-10-03 | 西南石油大学 | Method for quantitatively representing shale reservoir core hole seam configuration |
CN116840124B (en) * | 2023-06-21 | 2023-12-22 | 西南石油大学 | Method for quantitatively representing shale reservoir core hole seam configuration |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111624146A (en) | Method for quantitatively characterizing full-scale pore size distribution characteristics of compact reservoir | |
Song et al. | Pore structure characterization and permeability prediction of coal samples based on SEM images | |
CN105352873B (en) | The characterizing method of shale pore structure | |
US8909508B2 (en) | Petrographic image analysis for determining capillary pressure in porous media | |
Rabbani et al. | Comparing three image processing algorithms to estimate the grain-size distribution of porous rocks from binary 2D images and sensitivity analysis of the grain overlapping degree | |
Randle | Electron backscatter diffraction: Strategies for reliable data acquisition and processing | |
Cnudde et al. | 3D characterization of sandstone by means of X-ray computed tomography | |
Cerezo et al. | Aspects of the observation of clusters in the 3‐dimensional atom probe | |
CN110415167A (en) | A kind of rough surface crack generation method and pilot system based on Digital image technology | |
Ohm et al. | Size distribution of coarse-grained soil by sedimaging | |
Muszynski et al. | Particle shape estimates of uniform sands: visual and automated methods comparison | |
CN109444015A (en) | More kens, it is multiple dimensioned under shale reservoir microcellular system identification method | |
CN113158490B (en) | Method for establishing tight sandstone reservoir permeability calculation model | |
Lindqvist et al. | Image analysis applied to engineering geology, a literature review | |
Pal et al. | Porosity estimation by digital image analysis | |
Saur et al. | First investigation of quartz and calcite shape fabrics in strained shales by means of X-ray tomography | |
Germinario et al. | Textural and mineralogical analysis of volcanic rocks by µ-XRF mapping | |
CN113191330A (en) | Method for identifying area growth pores by fusing secondary electrons and back scattering electron images | |
Block et al. | Multifractal analysis of the microdistribution of elements in sedimentary structures using images from scanning electron microscopy and energy dispersive X ray spectrometry | |
Ghous et al. | 3D imaging of reservoir core at multiple scales: Correlations to petrophysical properties and pore scale fluid distributions | |
CN111626975A (en) | Method for quantitatively representing full-scale reservoir pores | |
WO2019074491A1 (en) | Method for patient stratification using image analysis parameter distribution functions | |
CN114518309A (en) | Method and device for effectively three-dimensionally representing interlayer pores in shale oil reservoir | |
Ghalib et al. | Soil stratigraphy delineation by VisCPT | |
EL‐ELLA et al. | Discrimination between depositional environments using grain‐size analyses |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200904 |