CN111028094B - Equivalent TOC calculation method based on shale scanning electron microscope image extraction - Google Patents
Equivalent TOC calculation method based on shale scanning electron microscope image extraction Download PDFInfo
- Publication number
- CN111028094B CN111028094B CN201911112480.XA CN201911112480A CN111028094B CN 111028094 B CN111028094 B CN 111028094B CN 201911112480 A CN201911112480 A CN 201911112480A CN 111028094 B CN111028094 B CN 111028094B
- Authority
- CN
- China
- Prior art keywords
- shale
- electron microscope
- equivalent
- sample
- image
- 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
- 238000001000 micrograph Methods 0.000 title claims abstract description 38
- 238000000605 extraction Methods 0.000 title claims abstract description 17
- 238000004364 calculation method Methods 0.000 title claims abstract description 10
- 238000000034 method Methods 0.000 claims abstract description 42
- 239000005416 organic matter Substances 0.000 claims abstract description 38
- 230000011218 segmentation Effects 0.000 claims abstract description 12
- 238000003384 imaging method Methods 0.000 claims abstract description 11
- 238000012545 processing Methods 0.000 claims abstract description 11
- 230000008569 process Effects 0.000 claims abstract description 8
- 238000007781 pre-processing Methods 0.000 claims abstract description 7
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims abstract description 5
- 229910052799 carbon Inorganic materials 0.000 claims abstract description 5
- 239000011028 pyrite Substances 0.000 claims description 40
- NIFIFKQPDTWWGU-UHFFFAOYSA-N pyrite Chemical compound [Fe+2].[S-][S-] NIFIFKQPDTWWGU-UHFFFAOYSA-N 0.000 claims description 40
- 229910052683 pyrite Inorganic materials 0.000 claims description 40
- 239000011148 porous material Substances 0.000 claims description 22
- 238000001914 filtration Methods 0.000 claims description 13
- 230000005484 gravity Effects 0.000 claims description 6
- 238000005498 polishing Methods 0.000 claims description 6
- 238000000227 grinding Methods 0.000 claims description 4
- 235000002566 Capsicum Nutrition 0.000 claims description 3
- 239000006002 Pepper Substances 0.000 claims description 3
- 235000016761 Piper aduncum Nutrition 0.000 claims description 3
- 235000017804 Piper guineense Nutrition 0.000 claims description 3
- 235000008184 Piper nigrum Nutrition 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 3
- 238000011049 filling Methods 0.000 claims description 3
- 150000003839 salts Chemical class 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 230000001174 ascending effect Effects 0.000 claims description 2
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 244000203593 Piper nigrum Species 0.000 claims 1
- 230000000694 effects Effects 0.000 claims 1
- 238000011160 research Methods 0.000 abstract description 7
- 238000012360 testing method Methods 0.000 abstract description 7
- 238000005259 measurement Methods 0.000 abstract description 5
- 238000003672 processing method Methods 0.000 abstract description 4
- 239000007789 gas Substances 0.000 description 6
- 230000009286 beneficial effect Effects 0.000 description 4
- 239000000126 substance Substances 0.000 description 4
- 230000008602 contraction Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000005192 partition Methods 0.000 description 3
- 241000722363 Piper Species 0.000 description 2
- -1 argon ion Chemical class 0.000 description 2
- 238000002485 combustion reaction Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000013384 organic framework Substances 0.000 description 2
- 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
- 238000012512 characterization method Methods 0.000 description 1
- 239000002734 clay mineral Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
- 238000010972 statistical evaluation Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/136—Segmentation; Edge detection involving thresholding
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/187—Segmentation; Edge detection involving region growing; involving region merging; involving connected component labelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20024—Filtering details
- G06T2207/20032—Median filtering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30181—Earth observation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Business, Economics & Management (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Marine Sciences & Fisheries (AREA)
- Mining & Mineral Resources (AREA)
- Agronomy & Crop Science (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- Animal Husbandry (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Geometry (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
The embodiment of the invention discloses an equivalent TOC (total organic carbon) calculation method based on shale scanning electron microscope image extraction, which comprises the following steps: processing a sample, scanning and imaging, and preprocessing an electron microscope image; threshold segmentation, namely dividing the elements of the shale sample electron microscope image into four components; identifying and correcting the white edge, and extracting the content of an organic matter framework; the method has the advantages that the shale organic matter skeleton content is evaluated by a digital image processing method, the equivalent TOC content of the shale is calculated, the shale reservoir gas storage capacity and the shale reservoir transport capacity can be accurately evaluated, meanwhile, the method is higher in universality, shale samples of different types and different electron microscope scanning imaging qualities can be researched, the research is facilitated, in addition, the method is short in measurement period, simple in processing process and capable of greatly reducing the test cost.
Description
Technical Field
The embodiment of the invention relates to the technical field of petroleum and natural gas geological exploration, in particular to an equivalent TOC (total organic carbon) calculation method based on shale scanning electron microscope image extraction.
Background
Shale gas has a huge reserve as a clean unconventional energy source. The extraction of the shale pore structure characteristics is the basis for carrying out accurate reserve evaluation and pore network migration capability calculation and is also the key point of shale gas exploration and development research. The shale has complex content of components and developed pore structure, and the main components of the shale are organic matters, inorganic matters and pyrite. Organic pores and organic cracks develop in the organic matters; inorganic pores and inorganic cracks develop in the inorganic substance.
The TOC content of shale is typically obtained by combustion testing and well logging. These are TOC analyses based on field tests, which are not only highly inaccurate, but also, more importantly, have a high cost and a long measurement period. Shale structure qualitative observation and quantitative characterization based on scanning electron microscope images are hot spots which are formed in recent years. The observation of a scanning electron microscope also reaches enough precision, the shale pore structure can be directly subjected to statistical evaluation, and information contained in the image can be effectively quantitatively extracted based on a reasonable image processing technology.
However, the existing method for extracting the TOC content of the shale has the following defects:
(1) the existing extraction method for the TOC content of shale has inaccurate extraction result and insufficient effective information mining, cannot accurately evaluate the gas storage capacity and the transportation capacity of a shale reservoir, and brings great inaccuracy to research work;
(2) the existing shale TOC content extraction method is low in applicability, cannot be used for researching shale samples of different types and different electron microscope scanning imaging qualities, and is not beneficial to research;
(3) the existing method for extracting the TOC content of the shale has the disadvantages of long measurement period, complex treatment process and high cost, and causes resource waste.
Disclosure of Invention
Therefore, the equivalent TOC calculation method based on shale scanning electron microscope image extraction is provided by the embodiment of the invention, the method adopts a digital image processing method to evaluate the shale organic matter skeleton content, and calculates the equivalent TOC content of shale, so that the shale reservoir gas storage capacity and transport capacity can be accurately evaluated, meanwhile, the method has higher universality, can be used for researching shale samples of different types and different electron microscope scanning imaging qualities, and is beneficial to research, in addition, the method has the advantages of short measurement period, simple processing process and lower cost, can greatly save manpower and material resources, and can effectively solve the problems in the prior art.
In order to achieve the above object, the embodiments of the present invention provide the following technical solutions: an equivalent TOC calculation method based on shale scanning electron microscope image extraction comprises the following steps:
s100, processing a sample, scanning and imaging, and preprocessing an electron microscope image;
s200, performing threshold segmentation, and dividing elements of the shale sample electron microscope image into four components;
s300, identifying and correcting the 'white edge', and extracting the content of an organic matter framework;
s400, converting the two-dimensional plane parameters into equivalent three-dimensional parameters by using ascending dimension conversion, and calculating the equivalent TOC content of the shale sample.
Further, in step S100, the method specifically includes:
s101, selecting a sample, and performing sample cutting, surface grinding and polishing treatment to obtain a sample slice with a size required by the observation of a scanning electron microscope;
s102, obtaining a series of scanning electron microscope images of the shale sample with a representative volume scale through scanning electron microscope observation;
and S103, preprocessing the image, cutting off the black edge of the image, and filtering the electron microscope image.
Further, in step S103, the method specifically includes:
extracting four vertex angles of a rectangular region with the image gray value not being 0, reserving the region in the rectangle, and cutting out the region outside the rectangle;
and eliminating the non-smoothness of the shape scanning of the cut image by mean filtering, and then eliminating the salt and pepper noise of the image by median filtering.
Further, the step of performing mean filtering on the removed image specifically includes: generating a 3 × 3 (or 5 × 5) template, calculating an average value of 8 (or 24) pixels (excluding the target point itself) around the target point, replacing the pixel value of the target point with the average value, and then sequentially calculating the average of all pixel points of the image to obtain a filtered image;
the median filtering performed on the removed image specifically includes: and generating a 3 × 3 (or 5 × 5) template, solving a median value of 8 (or 24) pixels (excluding the target point) around the target point, replacing the pixel value of the target point with the median value, and then solving the median values of all pixel points of the image in sequence to obtain the filtered image.
Further, in step S200, the method specifically includes:
s201, selecting a representative image in the cut images, namely an image containing more pores, organic matters, inorganic matters and pyrite;
s202, respectively adjusting and distinguishing three threshold values of the four components by using a binarization method, namely a pore threshold value, a segmentation threshold value of organic matters and inorganic matters and a segmentation threshold value of inorganic matters and pyrite;
and S203, after the three thresholds are obtained, performing threshold segmentation on the scanning electron microscope image of the sample, and dividing elements of each image into four components, namely holes, organic matters (including holes), inorganic matters (including holes) and pyrite.
Further, in step S300, the identification and correction of the white edge specifically includes:
s301, sequentially marking each communicated pyrite region in the sample image, counting the area of each pyrite region, and extracting the boundary coordinate, the barycentric coordinate, the long axis length, the short axis length and the equivalent radius of each pyrite region;
s302, calculating the ratio of the long axis to the short axis of each pyrite block;
s303, counting the average gray value of 9 pixel points in the 3 multiplied by 3 area size near the center of gravity point;
s304, if the ratio of the long axis to the short axis of the pyrite block is larger than 10 or the average gray value of 9 pixel points near the gravity center is smaller than the pyrite threshold, determining that the pyrite block is actually the interference of the 'white edge' at the periphery of the hole seam, and judging the real composition of the 'white edge'.
Further, when judging the real composition of the "white border":
finding out a minimum rectangle containing the block according to the boundary coordinates, and giving out coordinates of four vertexes of the rectangle;
expanding the equivalent radius size outwards from the extracted minimum rectangle to obtain a new rectangle as the neighborhood of the hole;
when a certain outwards expanded neighborhood edge exceeds the boundary, inwards shrinking according to the step length of 3 pixels until the shrinking of the boundary-free range is not exceeded for the first time, and the shrinking processes of all edges are independent;
counting the neighborhoods of the communicating blocks, and respectively calculating the content of organic matters and inorganic matters in each neighborhood, particularly the content of the organic matters and the inorganic matters in the region contacted with the pyrite region;
judging the real composition of the white edge according to the relative content of the organic matter component and the inorganic matter component in the neighborhood, if the organic matter content is high, considering the white edge as the organic matter, and filling the white edge area as the gray value of the organic matter; if the inorganic matter content is high, the real component of the white edge is considered to be inorganic matter, and the white edge area is filled with the gray value of the inorganic matter.
Further, in step S300, the specific steps of extracting the content of the organic framework are as follows:
s305, extracting and counting the number of pixels of which all gray values are greater than a pore threshold value and less than a basic threshold value in the corrected shale sample electron microscope image;
s306, counting the number N of the pixels in the S305 o Multiplying the actual area value of the single pixel representation to obtain the area S occupied by organic matters of the shale sample o ;
S o =N o ×S single
S307, counting the total pixel number N of the sample s Multiplying the single pixel representation actual area value to serve as the shale sample actual area S s 。
S s =N s ×S single
Further, in step S400, the transformation specifically includes:
the area parameter is converted according to the following formula:
wherein s is p Is a two-dimensional area, S v Is the three-dimensional equivalent area.
Further, in step S400, when calculating the equivalent TOC content of the shale sample:
s401, calculating the three-dimensional equivalent organic matter framework area of the organic matter framework, and dividing the three-dimensional equivalent organic matter framework area by the total area of the three-dimensional equivalent sample to obtain the three-dimensional equivalent organic matter framework content gamma o ;
γ o =S v,o /S v,s =s p,o /s p,s =N o /N s
S402, calculating the three-dimensional equivalent TOC content of the shale sample according to the following formula:
where ρ is k For the density of kerogen, 1.05g/cm is generally taken 3 ,ρ r The density of the shale is generally 2.55g/cm 3 ,ω k The carbon content of the kerogen is 70-85% generally.
The embodiment of the invention has the following advantages:
(1) according to the method, the content of the shale organic matter framework is evaluated by adopting a digital image processing method, and the equivalent TOC content of the shale is calculated, so that the gas storage capacity and the transportation capacity of a shale reservoir can be accurately evaluated;
(2) the method has higher universality, can be used for researching shale samples of different types and different electron microscope scanning imaging qualities, and is beneficial to research;
(3) the invention has short measuring period and simple processing process, and can greatly reduce the testing cost.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below. It should be apparent that the drawings in the following description are merely exemplary, and that other embodiments can be derived from the drawings provided by those of ordinary skill in the art without inventive effort.
FIG. 1 is a schematic view of the overall flow structure of the present invention;
Detailed Description
The present invention is described in terms of particular embodiments, other advantages and features of the invention will become apparent to those skilled in the art from the following disclosure, and it is to be understood that the described embodiments are merely exemplary of the invention and that it is not intended to limit the invention to the particular embodiments disclosed. All other embodiments, which can be obtained by a person skilled in the art without making any creative effort based on the embodiments in the present invention, belong to the protection scope of the present invention.
As shown in fig. 1, the invention provides an equivalent TOC calculation method based on shale scanning electron microscope image extraction, which comprises the following steps:
and S100, processing a sample, scanning and imaging, and preprocessing an electron microscope image.
In step S100, the method specifically includes:
and S101, selecting a sample, and performing sample cutting, surface grinding and polishing treatment to obtain a sample slice with the size required by the observation of a scanning electron microscope.
And S102, observing through a scanning electron microscope to obtain a series of scanning electron microscope images of the shale sample with a representative volume scale.
Cutting the obtained shale core into blocks of 1cm multiplied by 0.5cm, selecting a surface of 1cm multiplied by 0.5cm as an observation surface, grinding, and putting into argon ion polishing equipment for polishing. After polishing, the sample is placed in an observation bin of a scanning electron microscope device for scanning electron microscope observation, and a (REV) scanning electron microscope image reaching a representative size is obtained, wherein the size of the image in this embodiment is 398.5 μm × 398.5 μm.
And S103, preprocessing the image, cutting off the black edge of the image, and filtering the electron microscope image.
In step S103, the method specifically includes:
extracting four vertex angles of a rectangular region with the image gray value not being 0, reserving the region in the rectangle, and cutting out the region outside the rectangle; and eliminating the non-smoothness of the shape scanning of the cut image by mean filtering, and then eliminating the salt and pepper noise of the image by median filtering.
The mean filtering of the clipped image specifically includes: and generating a 3 × 3 (or 5 × 5) template, calculating an average value of 8 (or 24) pixels (excluding the target point itself) around the target point, replacing the pixel value of the target point with the average value, and then sequentially calculating the average value of all pixel points of the image to obtain the filtered image.
The median filtering of the clipped image specifically includes: and generating a 3 × 3 (or 5 × 5) template, solving a median value of 8 (or 24) pixels (excluding the target point) around the target point, replacing the pixel value of the target point with the median value, and then solving the median values of all pixel points of the image in sequence to obtain the filtered image.
And S200, performing threshold segmentation, and dividing elements of the shale sample electron microscope image into four components.
In step S200, the method specifically includes:
step S201, selecting a representative image in the cut images, namely an image containing more pores, organic matters, inorganic matters and pyrite;
step S202, three threshold values for distinguishing the four components, namely a pore threshold value, a partition threshold value of organic matters and inorganic matters and a partition threshold value of inorganic matters and pyrite, are respectively adjusted by a binarization method.
Adjusting and extracting the pore threshold, taking the gray critical values with different sizes as the threshold, and judging the quality of the gray value according to two standards: first, whether the extracted wells coincide with those seen by the naked eye; and whether the obtained face porosity is close to the measured porosity or not. The best gray value is selected as the final division threshold, and the final value of the pore threshold in this embodiment is 70.
When the organic matter and inorganic matter partition threshold is selected, the grayscale critical values with different sizes larger than the pore threshold are taken for extraction test, the grayscale value of the extracted organic matter region closest to the result recognized by naked eyes is the basic threshold for distinguishing organic matter from inorganic matter, and in this embodiment, the final value of the basic threshold is 135.
When the inorganic substance and pyrite segmentation threshold is selected, the grayscale critical values with different sizes larger than the basic threshold are taken for extraction test, and the grayscale value with the closest result of the extracted pyrite region and the result identified by naked eyes is used as the pyrite threshold for distinguishing the inorganic substance and the pyrite, and the final value of the pyrite threshold in the embodiment is 205.
The main components in the shale sample in this example are organic pores, inorganic pores, organic fissures, inorganic fissures, organic matter, inorganic matter, pyrite, and other clay minerals, where the gray values of the organic pores, inorganic pores, organic fissures, and inorganic fissures are the smallest (black) and cannot be distinguished by a single gray value, the gray values of the organic matter are larger (black gray), the gray values of the inorganic matter are larger (gray), and the gray values of the pyrite are the largest (off-white or white).
And S203, after the three thresholds are obtained, performing threshold segmentation on the scanning electron microscope image of the sample, dividing elements of each image into four components of holes, organic matters (including holes), inorganic matters (including holes) and pyrite, sequentially marking the connected blocks by using a bwleael method, and marking all pixels in each connected block as the same number.
And S300, identifying and correcting the white edges, and extracting the content of the organic matter skeleton.
In step S300, the identification and correction of the white edge specifically includes:
and S301, sequentially marking each communicated pyrite region in the sample image, counting the area of each pyrite region, and extracting the boundary coordinate, the barycentric coordinate, the long axis length, the short axis length and the equivalent radius of each pyrite region.
Step S302, a ratio of a long axis to a short axis of each pyrite block is calculated.
Step S303, counting the average gray value of 9 pixel points in the area size of 3 multiplied by 3 near the center of gravity.
And S304, judging, if the ratio of the long axis to the short axis of the pyrite block is more than 10, or the average gray value of 9 pixel points near the gravity center is less than the pyrite threshold, determining that the pyrite block is substantially the interference of the 'white edge' at the periphery of the hole seam, and judging the real composition of the 'white edge'.
Judging the real composition of the white edge:
finding out the minimum rectangle containing the block according to the boundary coordinates, and giving out coordinates of four vertexes of the rectangle; expanding the equivalent radius size outwards from the extracted minimum rectangle to obtain a new rectangle as the neighborhood of the hole; when a certain outward expansion neighborhood edge meets the exceeding boundary, the contraction is inwards performed according to the step length of 3 pixels until the first time that the contraction does not exceed the boundary range is met, and the contraction processes of all edges are independent.
Counting the neighborhoods of the communicating blocks, and respectively calculating the content of organic matters and inorganic matters in each neighborhood, particularly the content of the organic matters and the inorganic matters in the region contacted with the pyrite region; judging the real composition of the white edge according to the relative content of the organic matter component and the inorganic matter component in the neighborhood, if the organic matter content is high, considering the white edge as the organic matter, and filling the white edge area as the gray value of the organic matter; if the inorganic matter content is high, the real component of the "white edge" is considered to be inorganic matter, and the "white edge" region is filled with the gray value of the inorganic matter.
In step S300, the specific steps of extracting the organic framework content are as follows:
and S305, extracting and counting the number of pixels of which all gray values are greater than a hole threshold value and less than a basic threshold value in the sample image in the corrected shale sample electron microscope image.
Step 306, counting the number N of the pixels counted in the step S305 o Multiplying the actual area value of the single pixel representation to obtain the area S occupied by organic matters of the shale sample o ;
S o =N o ×S single
Step S307, counting the total pixel number N of the sample s Multiplying the actual area value of the single pixel representation to be used as the actual area S of the shale sample s 。
S s =N s ×S single
In step S400, the transformation specifically includes:
the area parameter is converted as follows:
wherein s is p Is a two-dimensional area, S v Is the three-dimensional equivalent area.
In step S400, when calculating the equivalent TOC content of the shale sample:
step S401, calculating the three-dimensional equivalent organic matter framework area of the organic matter framework, and dividing the three-dimensional equivalent organic matter framework area by the total area of the three-dimensional equivalent sample to obtain the three-dimensional equivalent organic matter framework content gamma o ;
γ o =S v,o /S v,s =s p,o /s p,s =N o /N s
Step S402, calculating the three-dimensional equivalent TOC content of the shale sample according to the following formula:
where ρ is k The density of kerogen is generally 1.05g/cm 3 ,ρ r Is density of shaleDegree, typically 2.55g/cm 3 ,ω k The carbon content of the kerogen is 70-85% generally.
The calculated equivalent TOC values for the shale samples were between 2.83% and 3.44%, while the TOC of the same shale samples as measured by combustion was 4.05%. Therefore, the TOC predicted value based on the scanning electron microscope image is equivalent to the TOC value of the shale sample measured by the physicochemical technology, namely the TOC predicted value and the TOC value of the shale sample have equivalent consistency on the TOC content, and the TOC predicted value and the TOC value can be contrastively analyzed with certain accuracy.
The method has the following advantages: according to the method, the content of the shale organic matter framework is evaluated by adopting a digital image processing method, and the equivalent TOC content of the shale is calculated, so that the gas storage capacity and the transportation capacity of a shale reservoir can be accurately evaluated; the method has higher universality, can be used for researching shale samples of different types and different electron microscope scanning imaging qualities, and is beneficial to research, the processing accuracy rate of the shale samples with better electron microscope scanning imaging can reach 95 percent by the method, and the processing accuracy rate of the shale samples with poorer imaging quality can also reach more than 85 percent; the method has the advantages of short measurement period, simple processing process, low cost, great saving of manpower and material resources, and good reliability, and the difference of the test accuracy in comparison with the traditional physical and chemical technology is within 10%.
Although the invention has been described in detail above with reference to a general description and specific examples, it will be apparent to one skilled in the art that modifications or improvements may be made thereto based on the invention. Accordingly, such modifications and improvements are intended to be within the scope of the invention as claimed.
Claims (8)
1. An equivalent TOC calculation method based on shale scanning electron microscope image extraction is characterized by comprising the following steps:
s100, processing a sample, scanning and imaging, and preprocessing an electron microscope image;
s200, performing threshold segmentation, and dividing elements of the shale sample electron microscope image into four components;
s300, identifying and correcting a 'white edge', and extracting the content of an organic matter framework, wherein the 'white edge' is the phenomenon that the gray value of the edge of a pore is abnormally increased due to the edge enhancement effect of a scanning electron microscope;
s400, converting the two-dimensional plane parameters into equivalent three-dimensional parameters by using ascending dimension conversion, and calculating the equivalent TOC content of the shale sample;
calculating the three-dimensional equivalent organic matter framework area of the organic matter framework, and dividing the three-dimensional equivalent organic matter framework area by the total three-dimensional equivalent sample area to obtain the three-dimensional equivalent organic matter framework content gamma o ;
γ o =S v,o /S v,s =s p,o /s p,s =N o /N s ;
Wherein N is o Is the number of pixels, N s Is the total number of pixels, s p,o Representing the two-dimensional area, s, corresponding to No pixels p,s Representing the two-dimensional area, S, corresponding to Ns total pixels of the sample v,o Represents the three-dimensional equivalent area, S, corresponding to No pixels v,s Representing the three-dimensional equivalent area corresponding to Ns total pixels of the sample;
calculating the three-dimensional equivalent TOC content of the shale sample according to the following formula:
where ρ is k Taking the density of kerogen to be 1.05g/cm 3 ,ρ r Taking the density of shale as 2.55g/cm 3 ,ω k 70-85% of kerogen as carbon content.
2. The method for calculating the equivalent TOC extracted based on the shale scanning electron microscope image according to claim 1, wherein in step S100, the method specifically comprises:
s101, selecting a sample, and performing sample cutting, surface grinding and polishing treatment to obtain a sample slice with a size required by the observation of a scanning electron microscope;
s102, obtaining a series of scanning electron microscope images of the shale sample with a representative volume scale through scanning electron microscope observation;
and S103, preprocessing the image, cutting off the black edge of the image, and filtering the electron microscope image.
3. The method for calculating the equivalent TOC extracted based on the shale scanning electron microscope image according to claim 2, wherein in step S103, the method specifically comprises:
extracting four vertex angles of a rectangular region with the image gray value not being 0, reserving the region in the rectangle, and cutting out the region outside the rectangle;
and eliminating the non-smoothness of the shape scanning of the cut image by mean filtering, and then eliminating the salt and pepper noise of the image by median filtering.
4. The method for calculating equivalent TOC extracted based on shale scanning electron microscope images as claimed in claim 1, wherein in step S200, the method specifically comprises:
s201, selecting a representative image in the cut image, namely an image containing more pores, organic matters, inorganic matters and pyrite;
s202, respectively adjusting and distinguishing three threshold values of the four components by using a binarization method, namely a pore threshold value, a segmentation threshold value of organic matters and inorganic matters and a segmentation threshold value of inorganic matters and pyrite;
and S203, after the three thresholds are obtained, performing threshold segmentation on the scanning electron microscope images of the samples, and dividing elements of each image into four components, namely holes, organic matters containing the holes, inorganic matters containing the holes and pyrite.
5. The method for calculating equivalent TOC based on shale scanning electron microscope image extraction according to claim 1, wherein in step S300, the identification and correction of the white edge specifically comprises:
s301, sequentially marking each communicated pyrite region in the sample image, counting the area of each pyrite region, and extracting the boundary coordinate, the barycentric coordinate, the long axis length, the short axis length and the equivalent radius of each pyrite region;
s302, calculating the ratio of the long axis to the short axis of each pyrite block;
s303, counting the average gray value of 9 pixel points in the 3 multiplied by 3 area size near the center of gravity point;
s304, if the ratio of the long axis to the short axis of the pyrite block is larger than 10 or the average gray value of 9 pixel points near the gravity center is smaller than a pyrite threshold, determining that the pyrite block is actually white edge interference at the periphery of the hole, and judging the real composition of the white edge.
6. The equivalent TOC calculation method based on shale scanning electron microscope image extraction according to claim 5, characterized in that when the real composition of "white edge" is judged:
finding out the minimum rectangle containing the block according to the boundary coordinates of the 'white mark' and giving out coordinates of four vertexes of the rectangle;
expanding the equivalent radius size outwards from the extracted minimum rectangle to obtain a new rectangle as the neighborhood of the hole;
when a certain outwards expanded neighborhood edge exceeds the boundary, inwards shrinking according to the step length of 3 pixels until the shrinking is stopped when the first requirement that the boundary range is not exceeded, and the shrinking processes of all edges are mutually independent;
counting the neighborhoods of the communicating blocks, and respectively calculating the content of organic matters and inorganic matters in the neighborhoods;
judging the real composition of the white edge according to the relative content of the organic matter component and the inorganic matter component in the neighborhood, if the organic matter content is high, considering the white edge as the organic matter, and filling the white edge area as the gray value of the organic matter; if the inorganic matter content is high, the real component of the white edge is considered to be inorganic matter, and the white edge area is filled with the gray value of the inorganic matter.
7. The method for calculating equivalent TOC based on shale scanning electron microscope image extraction according to claim 1, wherein in step S300, the specific steps of extracting organic matter skeleton content are as follows:
s305, extracting and counting the number of pixels of which all gray values are greater than a pore threshold value and less than a basic threshold value in the corrected shale sample electron microscope image;
s306, counting the number N of the pixels in the S305 o Multiplying the actual area value of the single pixel representation to obtain the area S occupied by organic matters of the shale sample o ;
S o =N o ×S single
S307, counting the total pixel number N of the sample s Multiplying the actual area value of the single pixel representation to be used as the actual area S of the shale sample s ;
S s =N s ×S single 。
8. The method for calculating equivalent TOC extracted based on shale scanning electron microscope images according to claim 1, wherein in step S400, the transformation specifically comprises:
the area parameter is converted as follows:
wherein s is p Is a two-dimensional area, S v Is the three-dimensional equivalent area.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911112480.XA CN111028094B (en) | 2019-11-14 | 2019-11-14 | Equivalent TOC calculation method based on shale scanning electron microscope image extraction |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911112480.XA CN111028094B (en) | 2019-11-14 | 2019-11-14 | Equivalent TOC calculation method based on shale scanning electron microscope image extraction |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111028094A CN111028094A (en) | 2020-04-17 |
CN111028094B true CN111028094B (en) | 2022-09-20 |
Family
ID=70201496
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911112480.XA Active CN111028094B (en) | 2019-11-14 | 2019-11-14 | Equivalent TOC calculation method based on shale scanning electron microscope image extraction |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111028094B (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111879804B (en) * | 2020-07-03 | 2021-10-08 | 中国石油大学(北京) | Method, device and system for determining organic carbon content of shale |
CN112348831B (en) * | 2020-11-05 | 2022-11-11 | 中国石油大学(华东) | Shale SEM image segmentation method based on machine learning |
CN114627083B (en) * | 2022-03-16 | 2023-11-03 | 贝光科技(苏州)有限公司 | Shale pore type identification method based on secondary electron signal image |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016007170A1 (en) * | 2014-07-11 | 2016-01-14 | Halliburton Energy Services, Inc. | Imaging a porous rock sample using a nanoparticle suspension |
CN106547028A (en) * | 2015-09-21 | 2017-03-29 | 中国石油化工股份有限公司 | The method and apparatus of prediction shale reservoir TOC |
US10190998B1 (en) * | 2018-08-29 | 2019-01-29 | Research Institute Of Petroleum Exploration & Development, Dagang Oil Field Of Cnpc | Method and device for evaluating and predicting a shale oil enrichment areas of fault lacustrine basins |
-
2019
- 2019-11-14 CN CN201911112480.XA patent/CN111028094B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016007170A1 (en) * | 2014-07-11 | 2016-01-14 | Halliburton Energy Services, Inc. | Imaging a porous rock sample using a nanoparticle suspension |
CN106547028A (en) * | 2015-09-21 | 2017-03-29 | 中国石油化工股份有限公司 | The method and apparatus of prediction shale reservoir TOC |
US10190998B1 (en) * | 2018-08-29 | 2019-01-29 | Research Institute Of Petroleum Exploration & Development, Dagang Oil Field Of Cnpc | Method and device for evaluating and predicting a shale oil enrichment areas of fault lacustrine basins |
Non-Patent Citations (4)
Title |
---|
"Porosity estimation in kerogen-bearing shale gas reservoirs";Hongyn Yu;《Journal of Natural Gas Science and Engineer》;20180430;全文 * |
"基于图像分析技术的页岩微观孔隙特征定性及定量表征";孙寅森;《地球科学进展》;20160731;全文 * |
"泥页岩有机显微组分的扫描电镜形貌特征及识别方法";焦淑静;《电子显微学报》;20180430;全文 * |
"浅谈页岩气测井有机碳含量的计算方法";钟雷文;《能源技术 能源与环境》;20170630;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN111028094A (en) | 2020-04-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111028094B (en) | Equivalent TOC calculation method based on shale scanning electron microscope image extraction | |
WO2019134252A1 (en) | Method and device for automated portrayal and accurate measurement of width of structural crack | |
CN110910411B (en) | Shale crack automatic extraction method with size self-adaption function | |
CN107923240B (en) | Method for determining porosity associated with organic matter in a well or formation | |
US9558588B2 (en) | Method for building a 3D model of a rock sample | |
CN110910444B (en) | RES scale shale equivalent three-dimensional pore parameter rapid extraction method | |
CN104237103A (en) | Quantitative characterization method and quantitative characterization device for pore connectivity | |
CN112819809B (en) | Method for quantifying morphology of mineral particles in rock | |
CN110414341A (en) | A kind of Blasting In The Open Mining heap measurement statistical method based on unmanned air vehicle technique | |
Idowu et al. | Wettability analysis using micro-CT, FESEM and QEMSCAN, and its applications to digital rock physics | |
Payton et al. | The influence of grain shape and size on the relationship between porosity and permeability in sandstone: A digital approach | |
Carbonneau | The threshold effect of image resolution on image‐based automated grain size mapping in fluvial environments | |
CN113298795A (en) | Shale oil and gas reservoir pore fine classification extraction method and device | |
CN116452613B (en) | Crack contour extraction method in geological survey | |
CN115829990B (en) | Natural crack identification method based on imaging logging image processing | |
Lyu et al. | An image processing approach to measuring the sphericity and roundness of fracturing proppants | |
CN106524998A (en) | Method for measuring tunnel linear exposed structural surface based on three-dimensional laser scanning technique | |
CN103761525B (en) | Micro-CT digital image processing-based component detection method of cathode carbon block for aluminium electrolysis | |
CN113177949B (en) | Large-size rock particle feature recognition method and device | |
Saadatfar et al. | Rock fabric and texture from digital core analysis | |
CN108765381B (en) | Hydraulic engineering crack extraction method and system and information processing terminal | |
CN115047013A (en) | Shale rock processing method and device | |
Ehsan et al. | Particle size analysis by image processing technique | |
CN115410049B (en) | Classification evaluation method and device for rock erosion degree | |
Jiu et al. | A method to extract the content, radius and specific surface area of maceral compositions in coal reservoirs based on image modeling |
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 |