CN105551018A - Extraction method facing each component of carbonate complex reservoir digital core - Google Patents
Extraction method facing each component of carbonate complex reservoir digital core Download PDFInfo
- Publication number
- CN105551018A CN105551018A CN201510875673.6A CN201510875673A CN105551018A CN 105551018 A CN105551018 A CN 105551018A CN 201510875673 A CN201510875673 A CN 201510875673A CN 105551018 A CN105551018 A CN 105551018A
- Authority
- CN
- China
- Prior art keywords
- image
- carbonatite
- threshold
- class
- component
- 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
Abstract
The present invention discloses an extraction method facing each component of a carbonate complex reservoir digital core. The method provided by the invention comprises: performing pretreatment of electrical imaging logging information through loading electrical imaging logging information, and obtaining an electrical imaging logging static image; generating a full-hole electrical imaging image through adoption of an image repairing method; converting the full-hole electrical imaging image to a gray image, obtaining a grouping threshold through adoption of a maximum variance method between a foreground and a background on the gray image, and automatically extracting the components of the carbonate complex reservoir according to the grouping threshold; arranging a threshold regulation sliding bar, and performing interactive and fine extraction of each component of three types of carbonate reservoirs through regulating and dividing the dual threshold. Through adoption of the maximum variance method between the foreground and the background, the components of the carbonate reservoir is automatically extracted, so that the accurate extraction of aperture part, the communication part and the skeleton part of the reservoir is realized. The method provided by the invention is high in precision and high in reliability, the construction of a three-dimensional digital core is achieved so as to allow a digital core model to be furthest the same as an actual rock.
Description
Technical field
The present invention relates to logging technology data processing field, be specifically related to a kind of each component extracting process towards carbonatite complicated reservoirs digital cores.
Background technology
Electric imaging logging data longitudinal frame is high, can clearly reflect well pass the various reservoir characteristics on stratum, extracted by each component of reservoir of research based on Image Logging Data, in conjunction with rock core casting body flake data, carry out 3-dimensional digital rock core structure, make digital cores model farthest identical with the rock of reality, for the reservoir that calculates to a nicety rock physical property from analyze different rock physicses respond between internal relation Research foundation is provided, each component (the hole of carbonate reservoir, connected component and skeleton) extract be 3-dimensional digital rock core build basis, very large on the modeling accuracy impact of 3-dimensional digital rock core.At present, mainly utilize electric imaging logging data, adopt image partition method to extract each component of reservoir.The result precision that it obtains is lower, and digital cores model has larger difference with actual rock, and therefore, anxious a kind of precision to be studied is high, and each component extracting process of carbonate reservoir that the digital cores model obtained is farthest identical with the rock of reality.
Summary of the invention
In view of this, be necessary to provide a kind of precision high, and each component extracting process of carbonate reservoir that the digital cores model obtained is farthest identical with the rock of reality.
Towards each component extracting process of carbonatite complicated reservoirs digital cores, it comprises the following steps:
S1, loading electric imaging logging data, carry out pre-service by electric imaging logging data, obtain electric imaging logging still image;
S2, employing image repair method generate full hole Electrical imaging image;
S3, full hole Electrical imaging image is converted to gray level image, component gray level image adopting varimax between two class carry out carbonate reservoir is extracted automatically;
S4, based on Visualc++6.0 development platform, threshold value is set and regulates slider bar, adopt man-machine interaction mode manual partial to regulate segmentation dual threshold, to the mutual meticulous extraction of each component of carbonatite three class reservoir.
Each component extracting process towards carbonatite complicated reservoirs digital cores of the present invention, it utilizes electric imaging logging data, the component adopting varimax between two class to carry out carbonate reservoir is extracted automatically, realizes accurately dividing the aperture sections of reservoir, connected component and skeleton portion extracting.It is high that the method has precision, the feature that reliability is strong.Solve carbonatite 3-dimensional digital rock core to build, make digital cores model farthest identical with the rock of reality.
Accompanying drawing explanation
Fig. 1 is the FB(flow block) of each component extracting process towards carbonatite complicated reservoirs digital cores of the present invention;
Fig. 2 is the sub-process block diagram of step S3 in Fig. 1;
Fig. 3 is another sub-process block diagram of step S3 in Fig. 1;
Fig. 4 is that full hole Electrical imaging image and each component of carbonatite three class reservoir extract result figure automatically; Wherein, the left side is original Electrical imaging image, and the right is that between two class, varimax carbonatite three class reservoir component extracts result automatically;
Fig. 5 is the mutual meticulous extraction result figure of each component of carbonatite three class reservoir; Wherein, the left side is dual threshold interactive interface, and the right is the mutual meticulous extraction result of carbonatite three class reservoir component.
Embodiment
Clearly understand to make object of the present invention, technical scheme and advantage, below in conjunction with drawings and Examples, the present invention is further elaborated, is to be understood that, specific embodiment described herein only in order to explain the present invention, is not intended to limit the present invention.
The invention provides a kind of each component extracting process towards carbonatite complicated reservoirs digital cores, as shown in Figure 1, it comprises the following steps:
S1, loading electric imaging logging data, electric imaging logging data is carried out pre-service, obtains electric imaging logging still image, concrete, described pre-service comprises pole plate alignment, equalization, generation the dynamic and stalic state image;
S2, employing image repair method generate full hole Electrical imaging image;
S3, full hole Electrical imaging image is converted to gray level image, component gray level image adopting varimax between two class carry out carbonate reservoir is extracted automatically;
S4, based on Visualc++6.0 development platform, threshold value is set and regulates slider bar, adopt man-machine interaction mode manual partial to regulate segmentation dual threshold, to the mutual meticulous extraction of each component of carbonatite three class reservoir.
Wherein, as shown in Figure 2, described step S3 comprises following sub-step:
Varimax between S31, use class, according to the maximum optimal threshold obtaining the classification of carbonatite complicated reservoirs digital cores of variance between prospect and background, becomes foreground image and background image two class by optimal threshold by Iamge Segmentation;
In S32, two classes after singulation, reuse varimax between class and calculate the optimal classification in the subclass of foreground image and background image respectively, and maximum variance between the class of trying to achieve two optimal classifications;
S33, choose the segmentation threshold of corresponding threshold value as carbonatite three class reservoir of maximum variance between two classes respectively.
Described carbonatite three class reservoir comprises aperture sections, connected component and skeleton part, to three class classification problems such as each components of carbonate reservoir, two threshold values are needed to classify, therefore, varimax between two class is proposed, between described pair of class, varimax is for carry out twice maximum variance between clusters separation to gray level image, belongs to one based on varimax between class towards three class problem partitioning algorithms.
Concrete, between two class in varimax, if t is the segmentation threshold of foreground image and background image, foreground image is counted and accounted for image scaled is w
0, average gray is u
0, background image is counted and accounted for image scaled is w
1, average gray is u
1;
The then overall average gray-scale value of described image:
u＝w
0*u
0+w
1*u
1(1)
The variance of described foreground image and background images:
g＝w
0*(u
0-u)*(u
0-u)+w
1(u
1-u)*(u
1-u)(2)
＝w
0*w
1*(u
0-u
1)*(u
0-u
1)
Wherein, when prospect and background inter-class variance maximum time, foreground image and background image difference maximum, corresponding gray scale t be two classes classification optimal threshold.
As shown in Figure 3, the component that between described employing two class, varimax carries out carbonate reservoir automatically extract comprise step by step following:
S31, according to optimal threshold t, carbonatite Electrical imaging Methods in Gray-level Still Images p is divided into p
1and p
2two parts;
In two parts after singulation, to p
1part adopts varimax definite threshold t between class
2, to p
2part adopts varimax definite threshold t between class
3, and according to threshold value t
2by p
1be divided into p
11and p
12, according to threshold value t
3by p
2be divided into p
21and p
22;
From t
2and t
3the maximum corresponding threshold value of both middle inter-class variances is chosen as optimal threshold t in two threshold values
1, by t
1with t, carbonatite Electrical imaging Methods in Gray-level Still Images is divided into three classes, i.e. carbonate reservoir component.
If t
1<t, then according to gray-scale value <t
1, gray-scale value >t, gray-scale value is between t
1and carbonatite Electrical imaging Methods in Gray-level Still Images is divided into three classes by three kinds of threshold conditions between t.
As shown in Figure 4, Fig. 4 is that full hole Electrical imaging image and each component of carbonatite three class reservoir extract result figure automatically; Wherein, the left side is original Electrical imaging image, and the right is that between two class, varimax carbonatite three class reservoir component extracts result automatically.
Further, between employing two class, varimax is determined on the basis of each component threshold value of carbonatite three class reservoir automatically, adopts man-machine interaction mode manual partial to regulate the segmentation dual threshold of each component, the mutual meticulous extraction of each component of carbonatite three class reservoir.
Concrete, after automatically extracting carbonatite three class reservoir component, if be unsatisfied with the result automatically extracted, can man-machine interaction mode be adopted, the segmentation dual threshold that manual adjustments automated manner is determined, to reach optimal segmentation effect.During manual adjustments dual threshold, can the result of the current segmentation of live preview.
As shown in Figure 5, Fig. 5 is the mutual meticulous extraction result figure of each component of carbonatite three class reservoir.The left side is dual threshold interactive interface, and the right is the mutual meticulous extraction result of carbonatite three class reservoir component.As can be seen from Fig. 4 and Fig. 5 comparing result: it is meticulousr that the inventive method finally extracts carbonatite three class reservoir component, more can reflect true carbonate reservoir, and the carbonatite digital cores model built with this is more accurate.
Each component extracting process towards carbonatite complicated reservoirs digital cores of the present invention, it utilizes electric imaging logging data, the component adopting varimax between two class to carry out carbonate reservoir is extracted automatically, realizes accurately dividing the aperture sections of reservoir, connected component and skeleton portion extracting.It is high that the method has precision, the feature that reliability is strong.Solve carbonatite 3-dimensional digital rock core to build, make digital cores model farthest identical with the rock of reality.
The foregoing is only preferred embodiment of the present invention, not in order to limit the present invention, within the spirit and principles in the present invention all, any amendment done, equivalent replacement, improvement etc., all should be included within protection scope of the present invention.
Claims (8)
1. towards each component extracting process of carbonatite complicated reservoirs digital cores, it is characterized in that, described each component extracting process towards carbonatite complicated reservoirs digital cores comprises the following steps:
S1, loading electric imaging logging data, carry out pre-service by electric imaging logging data, obtain electric imaging logging still image;
S2, employing image repair method generate full hole Electrical imaging image;
S3, full hole Electrical imaging image is converted to gray level image, gray level image adopts varimax between two class obtain packet threshold, automatically extracted by the component of packet threshold by carbonate reservoir;
S4, based on Visualc++6.0 development platform, threshold value is set and regulates slider bar, adopt man-machine interaction mode manual partial to regulate segmentation dual threshold, to the mutual meticulous extraction of each component of carbonatite three class reservoir.
2. each component extracting process towards carbonatite complicated reservoirs digital cores according to claim 1, it is characterized in that, described step S3 comprises following sub-step:
S31, use varimax between class, obtain the optimal threshold of carbonatite complicated reservoirs digital cores classification time maximum according to the variance of foreground image and background images, by optimal threshold, Iamge Segmentation is become foreground image and background image two class;
In S32, two classes after singulation, reuse maximum variance between clusters and calculate optimal classification in the subclass of foreground image and background image respectively, and maximum variance between the class of trying to achieve two optimal classifications;
S33, choose the segmentation threshold of the corresponding threshold value of variance that between two classes, maximum variance is maximum as carbonatite three class reservoir respectively.
3. each component extracting process towards carbonatite complicated reservoirs digital cores according to claim 2, is characterized in that,
The overall average gray-scale value of described image:
u＝w
0*u
0+w
1*u
1(1)
The variance of described foreground image and background images:
g＝w
0*(u
0-u)*(u
0-u)+w
1(u
1-u)*(u
1-u)(2)
＝w
0*w
1*(u
0-u
1)*(u
0-u
1)
Wherein: t is the segmentation threshold of foreground image and background image, foreground image is counted and accounted for image scaled is w
0, average gray is u
0, background image is counted and accounted for image scaled is w
1, average gray is u
1.
4. each component extracting process towards carbonatite complicated reservoirs digital cores according to claim 3, is characterized in that,
When prospect and background inter-class variance maximum time, corresponding gray scale t be two classes classification optimal threshold.
5. each component extracting process towards carbonatite complicated reservoirs digital cores according to claim 4, is characterized in that, the component that between described employing two class, varimax carries out carbonate reservoir automatically extract comprise step by step following:
S31, according to optimal threshold t, carbonatite Electrical imaging Methods in Gray-level Still Images p is divided into p
1and p
2two parts;
In S32, two parts after singulation, to p
1part adopts varimax definite threshold t between class
2, to p
2part adopts varimax definite threshold t between class
3, and according to threshold value t
2by p
1be divided into p
11and p
12, according to threshold value t
3by p
2be divided into p
21and p
22;
S33, from t
2and t
3the maximum corresponding threshold value of both middle inter-class variances is chosen as optimal threshold t in two threshold values
1, by t
1with t, carbonatite Electrical imaging Methods in Gray-level Still Images is divided into three classes, i.e. carbonate reservoir component.
6. a kind of each component extracting process towards carbonatite complicated reservoirs digital cores according to claim 5, is characterized in that:
If t
1<t, according to gray-scale value <t
1, gray-scale value >t, gray-scale value is between t
1and carbonatite Electrical imaging Methods in Gray-level Still Images is divided into three classes by three kinds of threshold conditions between t.
7. a kind of each component extracting process towards carbonatite complicated reservoirs digital cores according to claim 1, is characterized in that: carbonatite three class reservoir comprises aperture sections, connected component and skeleton part.
8. a kind of each component extracting process towards carbonatite complicated reservoirs digital cores according to claim 1, is characterized in that: described pre-service comprises pole plate alignment, equalization, generation the dynamic and stalic state image.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510875673.6A CN105551018B (en) | 2015-12-03 | 2015-12-03 | Towards each component extracting method of carbonate rock complicated reservoirs digital cores |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510875673.6A CN105551018B (en) | 2015-12-03 | 2015-12-03 | Towards each component extracting method of carbonate rock complicated reservoirs digital cores |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105551018A true CN105551018A (en) | 2016-05-04 |
| CN105551018B CN105551018B (en) | 2018-04-20 |
Family
ID=55830194
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510875673.6A Active CN105551018B (en) | 2015-12-03 | 2015-12-03 | Towards each component extracting method of carbonate rock complicated reservoirs digital cores |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105551018B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109285222A (en) * | 2018-09-06 | 2019-01-29 | 中国地质大学（北京） | The building of organic shale high-resolution digital rock core and analysis method |
| CN109389128A (en) * | 2018-08-24 | 2019-02-26 | 中国石油天然气股份有限公司 | Electric imaging logging image characteristic automatic extraction method and device |
| CN112414917A (en) * | 2020-11-03 | 2021-02-26 | 西安石油大学 | Method for dividing and characterizing organic pores and inorganic pores of shale oil reservoir |
| CN112881253A (en) * | 2019-11-29 | 2021-06-01 | 华北理工大学 | Method for determining pore area of rock casting body slice picture, face porosity calculation method and application of face porosity calculation method |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6381363B1 (en) * | 1999-03-15 | 2002-04-30 | Grass Valley (U.S.), Inc. | Histogram-based segmentation of images and video via color moments |
| US20070219769A1 (en) * | 2004-06-11 | 2007-09-20 | Otthein Herzog | Method and Device for Segmenting a Digital Representation of Cells |
| CN102194121A (en) * | 2010-03-04 | 2011-09-21 | 天津市天堰医教科技开发有限公司 | Application of improved maximum between-class variance method in tongue crack recognition |
| CN103278864A (en) * | 2013-05-10 | 2013-09-04 | 中国石油天然气股份有限公司 | Method and device for determining geologic feather parameters and distribution of hole seam type reservoir stratum |
| CN104376556A (en) * | 2014-10-31 | 2015-02-25 | 四川大学 | Rock CT image target segmentation method |
| CN104794709A (en) * | 2015-04-10 | 2015-07-22 | 四川大学 | Three-dimensional core image pore and throat segmentation method |
-
2015
- 2015-12-03 CN CN201510875673.6A patent/CN105551018B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6381363B1 (en) * | 1999-03-15 | 2002-04-30 | Grass Valley (U.S.), Inc. | Histogram-based segmentation of images and video via color moments |
| US20070219769A1 (en) * | 2004-06-11 | 2007-09-20 | Otthein Herzog | Method and Device for Segmenting a Digital Representation of Cells |
| CN102194121A (en) * | 2010-03-04 | 2011-09-21 | 天津市天堰医教科技开发有限公司 | Application of improved maximum between-class variance method in tongue crack recognition |
| CN103278864A (en) * | 2013-05-10 | 2013-09-04 | 中国石油天然气股份有限公司 | Method and device for determining geologic feather parameters and distribution of hole seam type reservoir stratum |
| CN104376556A (en) * | 2014-10-31 | 2015-02-25 | 四川大学 | Rock CT image target segmentation method |
| CN104794709A (en) * | 2015-04-10 | 2015-07-22 | 四川大学 | Three-dimensional core image pore and throat segmentation method |
Non-Patent Citations (1)
| Title |
|---|
| 李仪等: "多类类间最大方差算法在图像分割中的应用", 《小型微型计算机系统》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109389128A (en) * | 2018-08-24 | 2019-02-26 | 中国石油天然气股份有限公司 | Electric imaging logging image characteristic automatic extraction method and device |
| US11010629B2 (en) | 2018-08-24 | 2021-05-18 | Petrochina Company Limited | Method for automatically extracting image features of electrical imaging well logging, computer equipment and non-transitory computer readable medium |
| CN109285222A (en) * | 2018-09-06 | 2019-01-29 | 中国地质大学（北京） | The building of organic shale high-resolution digital rock core and analysis method |
| CN112881253A (en) * | 2019-11-29 | 2021-06-01 | 华北理工大学 | Method for determining pore area of rock casting body slice picture, face porosity calculation method and application of face porosity calculation method |
| CN112414917A (en) * | 2020-11-03 | 2021-02-26 | 西安石油大学 | Method for dividing and characterizing organic pores and inorganic pores of shale oil reservoir |
| CN112414917B (en) * | 2020-11-03 | 2023-09-01 | 西安石油大学 | Shale oil reservoir organic pore and inorganic pore dividing and characterizing method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105551018B (en) | 2018-04-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105551018A (en) | Extraction method facing each component of carbonate complex reservoir digital core | |
| CN105678347A (en) | Pedestrian detection method and device | |
| CN106846336A (en) | Extract foreground image, replace the method and device of image background | |
| CN104281858A (en) | Three-dimensional convolutional neutral network training method and video anomalous event detection method and device | |
| CN104517095B (en) | A kind of number of people dividing method based on depth image | |
| CN108596919B (en) | Automatic image segmentation method based on depth map | |
| CN103702098A (en) | In-depth extracting method of three-viewpoint stereoscopic video restrained by time-space domain | |
| CN104200453B (en) | Parallax image correcting method based on image segmentation and credibility | |
| CN110706269B (en) | Binocular vision SLAM-based dynamic scene dense modeling method | |
| EP3779775A1 (en) | Media processing method and related apparatus | |
| Shen et al. | Frustum VoxNet for 3D object detection from RGB-D or Depth images | |
| EP2887314A1 (en) | Video frame conversion method and apparatus | |
| CN104992454A (en) | Regionalized automatic-cluster-change image segmentation method | |
| CN104077727A (en) | Substation electric equipment machine account intelligent type-in method and system | |
| Li et al. | Quality index for view synthesis by measuring instance degradation and global appearance | |
| CN102622768A (en) | Depth-map gaining method of plane videos | |
| CN109049716A (en) | Generation method, device, electronic equipment and the storage medium of 3 D-printing illustraton of model | |
| CN104915952A (en) | Method for extracting local salient objects in depth image based on multi-way tree | |
| CN110378336A (en) | Semantic class mask method, device and the storage medium of target object in training sample | |
| CN103679699A (en) | Stereo matching method based on translation and combined measurement of salient images | |
| CN106952301B (en) | RGB-D image significance calculation method | |
| CN107563275A (en) | A kind of face identification method combined with BIM | |
| CN104050665A (en) | Method and device for estimating foreground dwell time in video image | |
| CN108090436A (en) | A kind of detection method of moving object, system and medium | |
| CN113361567A (en) | Image processing method, image processing device, electronic equipment and storage medium |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB03 | Change of inventor or designer information | ||
| CB03 | Change of inventor or designer information |
Inventor after: Xiao Xiaoling Inventor after: Zhang Xiang Inventor after: Qu Le Inventor after: Du Huanhong Inventor after: Yu Chunhao Inventor before: Zhang Xiang Inventor before: Xiao Xiaoling Inventor before: Qu Le Inventor before: Du Huanhong Inventor before: Yu Chunhao |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |