CN109447944A - Lithofacies identification method and system for carbonate rock - Google Patents
Lithofacies identification method and system for carbonate rock Download PDFInfo
- Publication number
- CN109447944A CN109447944A CN201811106779.XA CN201811106779A CN109447944A CN 109447944 A CN109447944 A CN 109447944A CN 201811106779 A CN201811106779 A CN 201811106779A CN 109447944 A CN109447944 A CN 109447944A
- Authority
- CN
- China
- Prior art keywords
- rock
- carbonate rock
- image
- lithofacies
- carbonate
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 239000011435 rock Substances 0.000 title claims abstract description 272
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 title claims abstract description 179
- 238000000034 method Methods 0.000 title claims abstract description 31
- 239000002245 particle Substances 0.000 claims description 16
- 238000004590 computer program Methods 0.000 claims description 9
- 208000035126 Facies Diseases 0.000 claims description 5
- 108010001267 Protein Subunits Proteins 0.000 claims description 2
- 238000010276 construction Methods 0.000 claims 2
- 239000004575 stone Substances 0.000 claims 2
- 238000004364 calculation method Methods 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 19
- 238000003384 imaging method Methods 0.000 description 4
- 229910000514 dolomite Inorganic materials 0.000 description 3
- 239000010459 dolomite Substances 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 239000013078 crystal Substances 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 239000000843 powder Substances 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 238000004090 dissolution Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012706 support-vector machine Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0004—Industrial image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/90—Determination of colour characteristics
Landscapes
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Quality & Reliability (AREA)
- Image Analysis (AREA)
Abstract
The invention provides a lithofacies identification method and system for carbonate rocks. The lithofacies identification method of the carbonate rock comprises the following steps: acquiring a color image of the carbonate rock; obtaining a gray level image of the carbonate rock according to the color image of the carbonate rock; obtaining a binary image of the carbonate rock according to the gray level image of the carbonate rock; determining black pixels in the binary image according to the pixel values of the binary image; obtaining the rock structure number according to the average value of the longitudinal number of the black pixels in the binary image and the average value of the transverse number of the black pixels in the binary image; the lithofacies of the carbonate rock can be identified according to the rock structure number, the lithofacies of the carbonate rock can be accurately identified, and the identification working efficiency is improved.
Description
Technical field
The present invention relates to geological exploration fields, and in particular, to a kind of Lithofacies Identification method and system of carbonate rock.
Background technique
The lithofacies Logging Identification Method of carbonate rock mainly has two major classes method: first is that the identification based on conventional logging data
Method, including qualitative cross plot method of identification, quantitative mathematical statistics method such as neural network, support vector machine and fuzzy theory
The methods of.Second is that using new Logging Technology data, such as electric imaging logging identifies lithofacies, Electrical imaging chart method etc..
The above method identifies that lithofacies, greatest problem are that conventional logging data resolution is lower using conventional logging data,
For thin layer, thin interbed can not be identified, and human factor is larger, and identification working efficiency is low.
Summary of the invention
The main purpose of the embodiment of the present invention is to provide a kind of Lithofacies Identification method and system of carbonate rock, with accurate
It identifies the lithofacies of carbonate rock, improves identification working efficiency.
To achieve the goals above, the embodiment of the present invention provides a kind of Lithofacies Identification method of carbonate rock, comprising:
Obtain the color image of carbonate rock；
The gray level image of carbonate rock is obtained according to the color image of carbonate rock；
The bianry image of carbonate rock is obtained according to the gray level image of carbonate rock；
The black picture element in bianry image is determined according to the pixel value of bianry image；
According to black picture element, the average value and black picture element of longitudinal number are lateral in bianry image in bianry image
The average value of number obtains rock structure number；
The lithofacies of carbonate rock are identified according to rock structure number.
The embodiment of the present invention also provides a kind of Lithofacies Identification system of carbonate rock, comprising:
Acquiring unit, for obtaining the color image of carbonate rock；
Gray level image unit obtains the gray level image of carbonate rock for the color image according to carbonate rock；
Bianry image unit obtains the bianry image of carbonate rock for the gray level image according to carbonate rock；
Black picture element determination unit determines the black picture element in bianry image for the pixel value according to bianry image；
Rock structure counting unit, for the average value and black according to black picture element longitudinal number in bianry image
The average value of pixel lateral number in bianry image, obtains rock structure number；
Facies unit, for identifying the lithofacies of carbonate rock according to rock structure number.
The embodiment of the present invention also provides a kind of computer equipment, including memory, processor and storage are on a memory simultaneously
The computer program that can be run on a processor, processor perform the steps of when executing computer program
Obtain the color image of carbonate rock；
The gray level image of carbonate rock is obtained according to the color image of carbonate rock；
The bianry image of carbonate rock is obtained according to the gray level image of carbonate rock；
The black picture element in bianry image is determined according to the pixel value of bianry image；
According to black picture element, the average value and black picture element of longitudinal number are lateral in bianry image in bianry image
The average value of number obtains rock structure number；
The lithofacies of carbonate rock are identified according to rock structure number.
The embodiment of the present invention also provides a kind of computer readable storage medium, is stored thereon with computer program, computer
It is performed the steps of when program is executed by processor
Obtain the color image of carbonate rock；
The gray level image of carbonate rock is obtained according to the color image of carbonate rock；
The bianry image of carbonate rock is obtained according to the gray level image of carbonate rock；
The black picture element in bianry image is determined according to the pixel value of bianry image；
According to black picture element, the average value and black picture element of longitudinal number are lateral in bianry image in bianry image
The average value of number obtains rock structure number；
The lithofacies of carbonate rock are identified according to rock structure number.
The Lithofacies Identification method and system of the carbonate rock of the embodiment of the present invention first obtain the color image of carbonate rock, root
The gray level image of carbonate rock is obtained according to the color image of carbonate rock, obtains carbonate further according to the gray level image of carbonate rock
Then the bianry image of rock determines the black picture element in bianry image according to the pixel value of bianry image, then according to black picture
The average value of element average value and black picture element lateral number in bianry image of longitudinal number in bianry image, obtains
Rock structure number finally identifies the lithofacies of carbonate rock according to rock structure number, can accurately identify the lithofacies of carbonate rock, mention
Height identification working efficiency.
Detailed description of the invention
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, embodiment will be described below
Needed in attached drawing be briefly described, it should be apparent that, the accompanying drawings in the following description is only of the invention some
Embodiment for those of ordinary skill in the art without creative efforts, can also be attached according to these
Figure obtains other attached drawings.
Fig. 1 is the flow chart of the Lithofacies Identification method of carbonate rock in the embodiment of the present invention；
Fig. 2 is the gray level image schematic diagram of carbonate rock in the embodiment of the present invention；
Fig. 3 is the bianry image schematic diagrames of carbonate rock in the embodiment of the present invention；
Fig. 4 a be in the embodiment of the present invention Mx13 well in the gray level image of 4605m-4608m；
Fig. 4 b be in the embodiment of the present invention Mx13 well in the bianry image of 4605m-4608m；
Fig. 5 a be in the embodiment of the present invention Mx13 well in the gray level image of 4592m-4595m；
Fig. 5 b be in the embodiment of the present invention Mx13 well in the bianry image of 4592m-4595m；
Fig. 6 A is the schematic diagram of micrite cloud rock in the embodiment of the present invention；
Fig. 6 B is the gray level image schematic diagram of micrite cloud rock in the embodiment of the present invention；
Fig. 6 C is the bianry image schematic diagrames of micrite cloud rock in the embodiment of the present invention；
Fig. 7 A is the schematic diagram of particle cloud rock in the embodiment of the present invention；
Fig. 7 B is the gray level image schematic diagram of particle cloud rock in the embodiment of the present invention；
Fig. 7 C is the bianry image schematic diagrames of particle cloud rock in the embodiment of the present invention；
Fig. 8 be in the embodiment of the present invention Mx13 well in the Lithofacies Identification result schematic diagram of 4575m-4625m；
Fig. 9 is the structural block diagram of the Lithofacies Identification system of carbonate rock in the embodiment of the present invention.
Specific embodiment
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete
Site preparation description, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on
Embodiment in the present invention, it is obtained by those of ordinary skill in the art without making creative efforts every other
Embodiment shall fall within the protection scope of the present invention.
The lithofacies of carbonate rock can not be accurately identified in view of the prior art, working efficiency is low, and the embodiment of the present invention provides one
The Lithofacies Identification method of kind carbonate rock improves identification working efficiency to accurately identify the lithofacies of carbonate rock.Below in conjunction with attached
The present invention is described in detail for figure.
Fig. 1 is the flow chart of the Lithofacies Identification method of carbonate rock in the embodiment of the present invention.Fig. 2 is in the embodiment of the present invention
The gray level image schematic diagram of carbonate rock.Fig. 3 is the bianry image schematic diagrames of carbonate rock in the embodiment of the present invention.Extremely such as Fig. 1
Shown in Fig. 3, the Lithofacies Identification method of carbonate rock may include:
S101: the color image of carbonate rock is obtained.
S102: the gray level image of carbonate rock is obtained according to the color image of carbonate rock.
S103: the bianry image of carbonate rock is obtained according to the gray level image of carbonate rock.
S104: the black picture element in bianry image is determined according to the pixel value of bianry image.
S105: according to the average value and black picture element of black picture element longitudinal number in bianry image in bianry image
The average value of middle transverse direction number, obtains rock structure number.
S106: the lithofacies of carbonate rock are identified according to rock structure number.
The executing subject of the Lithofacies Identification method of carbonate rock shown in FIG. 1 can be computer.Process as shown in Figure 1
It is found that the Lithofacies Identification method of the carbonate rock of the embodiment of the present invention first obtains the color image of carbonate rock, according to carbonate
The color image of rock obtains the gray level image of carbonate rock, obtains the two-value of carbonate rock further according to the gray level image of carbonate rock
Then image determines the black picture element in bianry image according to the pixel value of bianry image, then according to black picture element in two-value
The average value of average value and the black picture element lateral number in bianry image of longitudinal number, obtains rock structure in image
Number finally identifies the lithofacies of carbonate rock according to rock structure number, can accurately identify the lithofacies of carbonate rock, improves identification work
Make efficiency.
In one embodiment, S102 is specifically included: obtaining carbonate according to the color image of carbonate rock by following formula
The gray level image of rock:
Gray=R*0.299+G*0.587+B*0.114；
Wherein, Gray is gray value, and R is red value, and G is green value, and B is blue valve.
In one embodiment, S103 is specifically included: judging whether the gray value in each pixel of gray level image is greater than first
Preset value；When the gray value in each pixel of gray level image is greater than the first preset value, gray value is adjusted to 0, pixel value
It is 0；When the gray value in each pixel of gray level image is less than or equal to the first preset value, gray value is adjusted to 255, as
Element value is 1.
Wherein, the gray value of gray level image is commonly divided into 0 to 255, totally 256 ranks, wherein 0 most dark (completely black), 255
Most bright (Quan Bai), centre are the grey of transition.First preset value is 100.
In one embodiment, S104 is specifically included: whether the pixel value for judging pixel in bianry image is 0.When pixel value is 0
When, pixel is black picture element, low-resistivity region is represented, such as (such as Fig. 7 C such as bed boundary (as shown in Figure 6 C) or dissolution pore, cracks
It is shown).When pixel value is 1, pixel is white pixel, represents high resistivity portion, such as rock matrix.
In one embodiment, rock structure number is obtained by following formula:
Wherein, ARFNFMIFor rock structure number, xpixavgFor the average value of black picture element lateral number in bianry image,
ypixavgFor the average value of black picture element longitudinal number in bianry image.
In one embodiment, S106 is specifically included: judging whether rock structure number is greater than the second preset value；When rock structure number
When greater than the second preset value, the lithofacies of carbonate rock are micrite cloud lithofacies；When rock structure number is less than or equal to the second preset value
When, the lithofacies of carbonate rock are particle cloud lithofacies.
Wherein, the second preset value is 1.Micrite cloud lithofacies have Lamellar character, including micrite cloud rock, shale micrite cloud rock and
Mud powder crystal cloud rock etc.；Particle cloud lithofacies have blocky feature, including particle cloud rock, fine grain cloud rock, powder crystal cloud rock and residual particles
Yun Yan etc., predominantly coarse grained dolomite.
One of specific embodiment of the invention is as follows:
By taking the Sichuan Basin area Mo Xi Temple of the Dragon King group Mx13 well as an example, practical situations of the invention are elaborated:
1, Basic Geological situation and regional parameters situation:
The Sichuan Basin area Mo Xi Temple of the Dragon King group stratum is a set of Marine source rock, and lithology is mainly dolomite,
A small amount of argillaceous dolomite.Well logging Image-forming instrument (FMI, Formation are scanned using Schlumberger formation resistivity
MicroScanner Image) electric imaging logging instrument measurement electric imaging logging data, 75% or more measurement data coverage rate, completely
Sufficient application conditions.
2, step is realized:
Step 1: obtaining the color image of carbonate rock, the gray scale of carbonate rock is obtained according to the color image of carbonate rock
Image.
Step 2: judging whether the gray value in each pixel of gray level image is greater than 100；When each picture of gray level image
When gray value in element is greater than 100, gray value is adjusted to 0；When the gray value in each pixel of gray level image is less than or waits
When 100, gray value is adjusted to 255.Fig. 4 a be in the embodiment of the present invention Mx13 well in the gray level image of 4605m-4608m.
Fig. 4 b be in the embodiment of the present invention Mx13 well in the bianry image of 4605m-4608m.Fig. 5 a is Mx13 well in the embodiment of the present invention
In the gray level image of 4592m-4595m.Fig. 5 b be in the embodiment of the present invention Mx13 well in the bianry image of 4592m-4595m.Such as
Shown in Fig. 4 a to Fig. 5 b, the digital image processing functions that can use matlab software convert the gray level image of carbonate rock to
Gray level image specifically can be converted bianry image using rgb2gray function by the bianry image of carbonate rock.
Step 3: whether the pixel value for judging pixel in bianry image is 0.When pixel value is 0, pixel is black picture element,
It can determine the black picture element in bianry image.
Step 4: according to the average value and black picture element of black picture element longitudinal number in bianry image in bianry image
The average value of middle transverse direction number, obtains rock structure number.
Fig. 6 A is the schematic diagram of micrite cloud rock in the embodiment of the present invention.Fig. 6 B is the ash of micrite cloud rock in the embodiment of the present invention
Spend image schematic diagram.Fig. 6 C is the bianry image schematic diagrames of micrite cloud rock in the embodiment of the present invention.As shown in Fig. 6 A to Fig. 6 C, mud
Brilliant cloud lithofacies have thin interbed shape feature, and the average value of the lateral number of black picture element is 87 on bianry image, and longitudinal
The average value of number is 27, and rock structure number is greater than 1, is 3.119.Fig. 7 A is the signal of particle cloud rock in the embodiment of the present invention
Figure.Fig. 7 B is the gray level image schematic diagram of particle cloud rock in the embodiment of the present invention.Fig. 7 C is particle cloud rock in the embodiment of the present invention
Bianry image schematic diagrames.As shown in Fig. 7 A to Fig. 7 C, particle cloud lithofacies are blocky feature, the black picture element on bianry image
The average value of lateral number be 20, the average value of longitudinal number is 34, and rock structure number is 0.5893 less than 1.
Step 5: judging whether rock structure number is greater than 1；When rock structure number is greater than 1, the lithofacies of carbonate rock are mud
Brilliant Yun Yanxiang；When rock structure number is less than or equal to 1, the lithofacies of carbonate rock are particle cloud lithofacies.
Fig. 8 be in the embodiment of the present invention Mx13 well in the Lithofacies Identification result schematic diagram of 4575m-4625m.As shown in figure 8,
The 1st in Fig. 8 is depth track, and the 2nd to the 4th is Logging Curves road, wherein the 2nd is natural gamma GR and hole diameter
Curve CAL, the 3rd is neutron, density harmony radio frequency channel, and the 4th is depth resistivity RT and flushed zone resistivity RXO, and the 5th is rock
Heart lithofacies histogram, the 6th is the explanatory note of rock core lithofacies, that is, to the 5th explanation, and the 7th is rock structure number,
For numberical range from 0 to 8, intermediate vertical line is 1, and being greater than 1 is micrite cloud lithofacies, and being less than or equal to 1 is particle cloud lithofacies.8th
For the explanation to Lithofacies Types.
To sum up, the Lithofacies Identification method of the carbonate rock of the embodiment of the present invention first obtains the color image of carbonate rock, root
The gray level image of carbonate rock is obtained according to the color image of carbonate rock, obtains carbonate further according to the gray level image of carbonate rock
Then the bianry image of rock determines the black picture element in bianry image according to the pixel value of bianry image, then according to black picture
The average value of element average value and black picture element lateral number in bianry image of longitudinal number in bianry image, obtains
Rock structure number finally identifies the lithofacies of carbonate rock according to rock structure number, can accurately identify the lithofacies of carbonate rock, mention
Height identification working efficiency.
Based on the same inventive concept, the embodiment of the invention also provides a kind of Lithofacies Identification system of carbonate rock, due to
The principle that the system solves the problems, such as is similar to the Lithofacies Identification method of carbonate rock, therefore the implementation of the system may refer to method
Implementation, overlaps will not be repeated.
Fig. 9 is the structural block diagram of the Lithofacies Identification system of carbonate rock in the embodiment of the present invention.As shown in figure 9, carbonate
The Lithofacies Identification system of rock includes:
Acquiring unit, for obtaining the color image of carbonate rock；
Gray level image unit obtains the gray level image of carbonate rock for the color image according to carbonate rock；
Bianry image unit obtains the bianry image of carbonate rock for the gray level image according to carbonate rock；
Black picture element determination unit determines the black picture element in bianry image for the pixel value according to bianry image；
Rock structure counting unit, for the average value and black according to black picture element longitudinal number in bianry image
The average value of pixel lateral number in bianry image, obtains rock structure number；
Facies unit, for identifying the lithofacies of carbonate rock according to rock structure number.
In a kind of wherein embodiment, bianry image unit is specifically included:
Whether judgment sub-unit, the gray value in each pixel for judging gray level image are greater than the first preset value；
The first adjustment subelement, for gray value to be adjusted to 0；
Second adjustment subelement, for gray value to be adjusted to 255.
In a kind of wherein embodiment, rock structure number is obtained by following formula:
Wherein, ARFNFMIFor rock structure number, xpixavgFor the average value of black picture element lateral number in bianry image,
ypixavgFor the average value of black picture element longitudinal number in bianry image.
In a kind of wherein embodiment, facies unit is specifically used for:
Judge whether rock structure number is greater than the second preset value；
When rock structure number is greater than the second preset value, the lithofacies of carbonate rock are micrite cloud lithofacies；
When rock structure number is less than or equal to the second preset value, the lithofacies of carbonate rock are particle cloud lithofacies.
To sum up, the Lithofacies Identification system of the carbonate rock of the embodiment of the present invention first obtains the color image of carbonate rock, root
The gray level image of carbonate rock is obtained according to the color image of carbonate rock, obtains carbonate further according to the gray level image of carbonate rock
Then the bianry image of rock determines the black picture element in bianry image according to the pixel value of bianry image, then according to black picture
The average value of element average value and black picture element lateral number in bianry image of longitudinal number in bianry image, obtains
Rock structure number finally identifies the lithofacies of carbonate rock according to rock structure number, can accurately identify the lithofacies of carbonate rock, mention
Height identification working efficiency.
The embodiment of the invention also provides a kind of computer equipments, including memory, processor and storage are on a memory
And the computer program that can be run on a processor, processor perform the steps of when executing computer program
Obtain the color image of carbonate rock；
The gray level image of carbonate rock is obtained according to the color image of carbonate rock；
The bianry image of carbonate rock is obtained according to the gray level image of carbonate rock；
The black picture element in bianry image is determined according to the pixel value of bianry image；
According to black picture element, the average value and black picture element of longitudinal number are lateral in bianry image in bianry image
The average value of number obtains rock structure number；
The lithofacies of carbonate rock are identified according to rock structure number.
To sum up, the computer equipment of the embodiment of the present invention first obtains the color image of carbonate rock, according to carbonate rock
Color image obtains the gray level image of carbonate rock, obtains the binary map of carbonate rock further according to the gray level image of carbonate rock
Then picture determines the black picture element in bianry image according to the pixel value of bianry image, then according to black picture element in binary map
The average value of average value and the black picture element lateral number in bianry image of longitudinal number, obtains rock structure number as in,
The lithofacies that carbonate rock is finally identified according to rock structure number, can accurately identify the lithofacies of carbonate rock, improve identification work
Efficiency.
The embodiment of the invention also provides a kind of computer readable storage mediums, are stored thereon with computer program, calculate
Machine program performs the steps of when being executed by processor
Obtain the color image of carbonate rock；
The gray level image of carbonate rock is obtained according to the color image of carbonate rock；
The bianry image of carbonate rock is obtained according to the gray level image of carbonate rock；
The black picture element in bianry image is determined according to the pixel value of bianry image；
According to black picture element, the average value and black picture element of longitudinal number are lateral in bianry image in bianry image
The average value of number obtains rock structure number；
The lithofacies of carbonate rock are identified according to rock structure number.
To sum up, the computer readable storage medium of the embodiment of the present invention first obtains the color image of carbonate rock, according to carbon
The color image of Carbonate Rocks obtains the gray level image of carbonate rock, obtains carbonate rock further according to the gray level image of carbonate rock
Then bianry image determines the black picture element in bianry image according to the pixel value of bianry image, exists then according to black picture element
The average value of average value and the black picture element lateral number in bianry image of longitudinal number, obtains rock in bianry image
Number is constructed, the lithofacies of carbonate rock are finally identified according to rock structure number, the lithofacies of carbonate rock can be accurately identified, improves and knows
Other working efficiency.
Particular embodiments described above has carried out further in detail the purpose of the present invention, technical scheme and beneficial effects
Describe in detail it is bright, it should be understood that the above is only a specific embodiment of the present invention, the guarantor being not intended to limit the present invention
Range is protected, all within the spirits and principles of the present invention, any modification, equivalent substitution, improvement and etc. done should be included in this
Within the protection scope of invention.
Claims (10)
1. a kind of Lithofacies Identification method of carbonate rock characterized by comprising
Obtain the color image of carbonate rock；
The gray level image of carbonate rock is obtained according to the color image of the carbonate rock；
The bianry image of carbonate rock is obtained according to the gray level image of the carbonate rock；
The black picture element in the bianry image is determined according to the pixel value of the bianry image；
According to the average value and the black picture element of the black picture element longitudinal number in the bianry image described two
It is worth the average value of lateral number in image, obtains rock structure number；
The lithofacies of carbonate rock are identified according to the rock structure number.
2. the Lithofacies Identification method of carbonate rock according to claim 1, which is characterized in that according to the carbonate rock
Gray level image obtains the bianry image of carbonate rock, specifically includes:
Judge whether the gray value in each pixel of the gray level image is greater than the first preset value；
When the gray value in each pixel of the gray level image is greater than first preset value, the gray value is adjusted to
0；
When the gray value in each pixel of the gray level image is less than or equal to first preset value, by the gray value
It is adjusted to 255.
3. the Lithofacies Identification method of carbonate rock according to claim 1, which is characterized in that obtain rock by following formula
Stone construction number:
Wherein, ARFNFMIFor rock structure number, xpixavgFor the average value of black picture element lateral number in the bianry image,
ypixavgFor the average value of black picture element longitudinal number in the bianry image.
4. the Lithofacies Identification method of carbonate rock according to claim 1, which is characterized in that according to the rock structure number
The lithofacies for identifying carbonate rock, specifically include:
Judge whether the rock structure number is greater than the second preset value；
When the rock structure number is greater than second preset value, the lithofacies of the carbonate rock are micrite cloud lithofacies；
When the rock structure number is less than or equal to second preset value, the lithofacies of the carbonate rock are particle cloud rock
Phase.
5. a kind of Lithofacies Identification system of carbonate rock characterized by comprising
Acquiring unit, for obtaining the color image of carbonate rock；
Gray level image unit, for obtaining the gray level image of carbonate rock according to the color image of the carbonate rock；
Bianry image unit, for obtaining the bianry image of carbonate rock according to the gray level image of the carbonate rock；
Black picture element determination unit, for determining the black picture in the bianry image according to the pixel value of the bianry image
Element；
Rock structure counting unit, for the average value according to the black picture element longitudinal number in the bianry image, and
The average value of the black picture element lateral number in the bianry image, obtains rock structure number；
Facies unit, for identifying the lithofacies of carbonate rock according to the rock structure number.
6. the Lithofacies Identification system of carbonate rock according to claim 5, which is characterized in that the bianry image unit tool
Body includes:
Whether judgment sub-unit, the gray value in each pixel for judging the gray level image are greater than the first preset value；
The first adjustment subelement, for the gray value to be adjusted to 0；
Second adjustment subelement, for the gray value to be adjusted to 255.
7. the Lithofacies Identification system of carbonate rock according to claim 5, which is characterized in that obtain rock by following formula
Stone construction number:
Wherein, ARFNFMIFor rock structure number, xpixavgFor the average value of black picture element lateral number in the bianry image,
ypixavgFor the average value of black picture element longitudinal number in the bianry image.
8. the Lithofacies Identification system of carbonate rock according to claim 5, which is characterized in that the facies unit is specifically used
In:
Judge whether the rock structure number is greater than the second preset value；
When the rock structure number is greater than second preset value, the lithofacies of the carbonate rock are micrite cloud lithofacies；
When the rock structure number is less than or equal to second preset value, the lithofacies of the carbonate rock are particle cloud rock
Phase.
9. a kind of computer equipment including memory, processor and stores the meter that can be run on a memory and on a processor
Calculation machine program, which is characterized in that the processor performs the steps of when executing the computer program
Obtain the color image of carbonate rock；
The gray level image of carbonate rock is obtained according to the color image of the carbonate rock；
The bianry image of carbonate rock is obtained according to the gray level image of the carbonate rock；
The black picture element in the bianry image is determined according to the pixel value of the bianry image；
According to the average value and the black picture element of the black picture element longitudinal number in the bianry image described two
It is worth the average value of lateral number in image, obtains rock structure number；
The lithofacies of carbonate rock are identified according to the rock structure number.
10. a kind of computer readable storage medium, is stored thereon with computer program, which is characterized in that the computer program
It is performed the steps of when being executed by processor
Obtain the color image of carbonate rock；
The gray level image of carbonate rock is obtained according to the color image of the carbonate rock；
The bianry image of carbonate rock is obtained according to the gray level image of the carbonate rock；
The black picture element in the bianry image is determined according to the pixel value of the bianry image；
According to the average value and the black picture element of the black picture element longitudinal number in the bianry image described two
It is worth the average value of lateral number in image, obtains rock structure number；
The lithofacies of carbonate rock are identified according to the rock structure number.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811106779.XA CN109447944B (en) | 2018-09-21 | 2018-09-21 | Lithofacies identification method and system for carbonate rock |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811106779.XA CN109447944B (en) | 2018-09-21 | 2018-09-21 | Lithofacies identification method and system for carbonate rock |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109447944A true CN109447944A (en) | 2019-03-08 |
| CN109447944B CN109447944B (en) | 2020-08-11 |
Family
ID=65530649
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811106779.XA Active CN109447944B (en) | 2018-09-21 | 2018-09-21 | Lithofacies identification method and system for carbonate rock |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109447944B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110469327A (en) * | 2019-08-02 | 2019-11-19 | 中国石油天然气集团有限公司 | A kind of approximant three-dimensional evaluation method of carbonate reservoir well logging of five steps |
Citations (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1563891A (en) * | 2004-04-20 | 2005-01-12 | 长安大学 | System and method for discriminating road gap |
| CN101950359A (en) * | 2010-10-08 | 2011-01-19 | 郝红卫 | Method for recognizing rock type |
| CN103257081A (en) * | 2013-04-28 | 2013-08-21 | 北京大学 | Method and device for recovering oil and gas reservoir rock mechanics underground in-situ model |
| CN103267678A (en) * | 2013-04-28 | 2013-08-28 | 北京大学 | Synchronous measurement method and apparatus for dynamic and static rock mechanical parameters |
| CN103278614A (en) * | 2013-04-28 | 2013-09-04 | 北京大学 | Method and device for correcting dynamic and static rock mechanical parameters |
| CN103325117A (en) * | 2013-06-17 | 2013-09-25 | 中国石油天然气集团公司 | Rock core picture processing method and system based on MATLAB |
| CN103698811A (en) * | 2013-12-24 | 2014-04-02 | 中国石油天然气股份有限公司 | Carbonate rock structure component well logging quantitative identification method and application thereof |
| CN104535475A (en) * | 2015-01-08 | 2015-04-22 | 中国石油大学(北京) | Carbonate rock microstructure determination method and device |
| US20150317052A1 (en) * | 2014-05-02 | 2015-11-05 | Lhb Soluções Em Informações E Métodos Ltda. | System and process to describe vertical sequences of rocks using gestures |
| CN105431750A (en) * | 2013-06-27 | 2016-03-23 | 地球物理集团公司 | Fracability measurement method and system |
| CN105911603A (en) * | 2016-05-04 | 2016-08-31 | 湖南科技大学 | Natural electric field based four-dimensional geophysical prospecting method |
| CN107133630A (en) * | 2016-02-29 | 2017-09-05 | 中国石油化工股份有限公司 | A kind of method that carbonate porosity type is judged based on scan image |
| CN107239783A (en) * | 2017-06-13 | 2017-10-10 | 中国矿业大学(北京) | Coal-rock identification method based on extension local binary patterns and regression analysis |
| CN107966546A (en) * | 2017-11-21 | 2018-04-27 | 西南石油大学 | A kind of shale lithofacies plane distribution preparation method and shale exploration system |
| CN108072916A (en) * | 2016-11-16 | 2018-05-25 | 中国石油化工股份有限公司 | Lacustrine facies mud shale lithofacies Logging Identification Method based on genetic analysis |
| CN108072748A (en) * | 2016-11-07 | 2018-05-25 | 中国石油天然气股份有限公司 | Carbonate rock lithofacies characterization method based on dominant lithology |
| CN108121013A (en) * | 2017-11-13 | 2018-06-05 | 中国石油天然气股份有限公司 | Method for identifying microbial carbonate rock lithofacies |
-
2018
- 2018-09-21 CN CN201811106779.XA patent/CN109447944B/en active Active
Patent Citations (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1563891A (en) * | 2004-04-20 | 2005-01-12 | 长安大学 | System and method for discriminating road gap |
| CN101950359A (en) * | 2010-10-08 | 2011-01-19 | 郝红卫 | Method for recognizing rock type |
| CN103257081A (en) * | 2013-04-28 | 2013-08-21 | 北京大学 | Method and device for recovering oil and gas reservoir rock mechanics underground in-situ model |
| CN103267678A (en) * | 2013-04-28 | 2013-08-28 | 北京大学 | Synchronous measurement method and apparatus for dynamic and static rock mechanical parameters |
| CN103278614A (en) * | 2013-04-28 | 2013-09-04 | 北京大学 | Method and device for correcting dynamic and static rock mechanical parameters |
| CN103325117A (en) * | 2013-06-17 | 2013-09-25 | 中国石油天然气集团公司 | Rock core picture processing method and system based on MATLAB |
| CN103325117B (en) * | 2013-06-17 | 2016-08-10 | 中国石油天然气集团公司 | A kind of core image processing method based on MATLAB and system |
| CN105431750A (en) * | 2013-06-27 | 2016-03-23 | 地球物理集团公司 | Fracability measurement method and system |
| CN103698811A (en) * | 2013-12-24 | 2014-04-02 | 中国石油天然气股份有限公司 | Carbonate rock structure component well logging quantitative identification method and application thereof |
| US20150317052A1 (en) * | 2014-05-02 | 2015-11-05 | Lhb Soluções Em Informações E Métodos Ltda. | System and process to describe vertical sequences of rocks using gestures |
| CN104535475A (en) * | 2015-01-08 | 2015-04-22 | 中国石油大学(北京) | Carbonate rock microstructure determination method and device |
| CN107133630A (en) * | 2016-02-29 | 2017-09-05 | 中国石油化工股份有限公司 | A kind of method that carbonate porosity type is judged based on scan image |
| CN105911603A (en) * | 2016-05-04 | 2016-08-31 | 湖南科技大学 | Natural electric field based four-dimensional geophysical prospecting method |
| CN108072748A (en) * | 2016-11-07 | 2018-05-25 | 中国石油天然气股份有限公司 | Carbonate rock lithofacies characterization method based on dominant lithology |
| CN108072916A (en) * | 2016-11-16 | 2018-05-25 | 中国石油化工股份有限公司 | Lacustrine facies mud shale lithofacies Logging Identification Method based on genetic analysis |
| CN107239783A (en) * | 2017-06-13 | 2017-10-10 | 中国矿业大学(北京) | Coal-rock identification method based on extension local binary patterns and regression analysis |
| CN108121013A (en) * | 2017-11-13 | 2018-06-05 | 中国石油天然气股份有限公司 | Method for identifying microbial carbonate rock lithofacies |
| CN107966546A (en) * | 2017-11-21 | 2018-04-27 | 西南石油大学 | A kind of shale lithofacies plane distribution preparation method and shale exploration system |
Non-Patent Citations (5)
| Title |
|---|
| LEE JUNG MO 等: "Seismic anisotropy caused by rock fabric", 《GEOPHYS. J. INT》 * |
| 于炳松 等: "《岩石学》", 30 October 2017 * |
| 冯国庆 等: "利用多点地质统计学方法模拟岩相分布", 《西安石油大学学报（自然科学版）》 * |
| 李昌 等: "FMI测井技术在川东地区碳酸盐岩溶孔溶洞型储层评价中的应用", 《海相油气地质》 * |
| 李昌 等: "视岩石结构数技术在测井识别碳酸盐岩岩相中的应用", 《油气地球物理》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110469327A (en) * | 2019-08-02 | 2019-11-19 | 中国石油天然气集团有限公司 | A kind of approximant three-dimensional evaluation method of carbonate reservoir well logging of five steps |
| CN110469327B (en) * | 2019-08-02 | 2022-12-02 | 中国石油天然气集团有限公司 | Five-step approximation type carbonate reservoir well logging three-dimensional evaluation method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109447944B (en) | 2020-08-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| RU2543698C1 (en) | Analysis of petrographic images for detection of capillary pressure in porous media | |
| US20180274364A1 (en) | Method and Device for Determining Karst Development Degree of Reservoir, Computer Readable Storage Medium and Device | |
| MX2007000380A (en) | Computer-based method for while-drilling modeling and visualization of layered subterranean earth formations. | |
| CN108363110B (en) | Spectral analysis method for calculating shale reservoir mineral content and brittleness index through imaging logging | |
| CN110208859A (en) | Oil-base mud well crack quantitative parameter intelligence computation method based on ultrasonic imaging | |
| CN102011583A (en) | Method for identifying reservoir by combining electric imaging with reef geologic model | |
| CN108109157A (en) | A kind of rock mass analysis and assessment method based on digital panorama borehole image | |
| Yamada et al. | Revisiting porosity analysis from electrical borehole images: integration of advanced texture and porosity analysis | |
| CN110837114B (en) | Method and device for identifying rough-surface pyroclastic rock and electronic equipment | |
| CN110147778A (en) | Rare Earth Mine exploits recognition methods, device, equipment and storage medium | |
| Zhang et al. | Research on coal and rock type recognition based on mechanical vision | |
| US10152783B2 (en) | Method and system for two dimensional sediment particle shape characterization based on borehole image | |
| CN109447944A (en) | Lithofacies identification method and system for carbonate rock | |
| US9563602B2 (en) | Method of analyzing 3D geological structure using structure index | |
| Linek et al. | Rock classification based on resistivity patterns in electrical borehole wall images | |
| CN111766637B (en) | Lithology quantitative spectrum method for identifying lithology of tight reservoir | |
| Xian et al. | Fracture Characterisation Through Image Log Conductive Feature Extraction | |
| CN105425296B (en) | Geologic body recognition methods and device | |
| Blondel et al. | Textural analysis and structure-tracking for geological mapping: applications to sonar images from Endeavour Segment, Juan de Fuca Ridge | |
| CN111311559A (en) | Non-deviation crack characterization method | |
| CN109948291A (en) | A kind of discontinuous boundary line direction-adaptive recognition methods of sand body | |
| You et al. | Automatic facies classification from acoustic image logs using deep neural networks | |
| CN106780332A (en) | Full hole well logging video generation device | |
| CN108828687A (en) | A kind of calculation of permeability based on Electrical imaging Areal porosity | |
| CN104062681A (en) | Seismic horizon tracking preprocessing method based on fractional derivative |
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 |