CN116246007A - Microscopic residual oil type conversion quantitative analysis method based on CT technology - Google Patents
Microscopic residual oil type conversion quantitative analysis method based on CT technology Download PDFInfo
- Publication number
- CN116246007A CN116246007A CN202111493785.7A CN202111493785A CN116246007A CN 116246007 A CN116246007 A CN 116246007A CN 202111493785 A CN202111493785 A CN 202111493785A CN 116246007 A CN116246007 A CN 116246007A
- Authority
- CN
- China
- Prior art keywords
- residual oil
- clusters
- oil
- microscopic
- type
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 116
- 238000006243 chemical reaction Methods 0.000 title claims abstract description 96
- 238000004445 quantitative analysis Methods 0.000 title claims abstract description 35
- 238000005516 engineering process Methods 0.000 title abstract description 9
- 238000006073 displacement reaction Methods 0.000 claims abstract description 158
- 230000008569 process Effects 0.000 claims abstract description 49
- 238000002591 computed tomography Methods 0.000 claims abstract description 41
- 239000011148 porous material Substances 0.000 claims abstract description 40
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 58
- 230000011218 segmentation Effects 0.000 claims description 20
- 238000012545 processing Methods 0.000 claims description 19
- 230000015654 memory Effects 0.000 claims description 15
- 229920006395 saturated elastomer Polymers 0.000 claims description 15
- 230000009466 transformation Effects 0.000 claims description 15
- 238000000926 separation method Methods 0.000 claims description 12
- 238000003860 storage Methods 0.000 claims description 11
- 239000000126 substance Substances 0.000 claims description 11
- 239000003795 chemical substances by application Substances 0.000 claims description 10
- 230000009467 reduction Effects 0.000 claims description 8
- 238000007781 pre-processing Methods 0.000 claims description 6
- 238000004590 computer program Methods 0.000 claims description 3
- 230000009286 beneficial effect Effects 0.000 abstract description 4
- 239000003921 oil Substances 0.000 description 326
- 238000004364 calculation method Methods 0.000 description 9
- 239000003086 colorant Substances 0.000 description 6
- 230000006870 function Effects 0.000 description 5
- 230000006798 recombination Effects 0.000 description 5
- 238000005215 recombination Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 4
- 230000018109 developmental process Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000009826 distribution Methods 0.000 description 4
- 239000011435 rock Substances 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 230000002776 aggregation Effects 0.000 description 3
- 238000004220 aggregation Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- NLKNQRATVPKPDG-UHFFFAOYSA-M potassium iodide Chemical compound [K+].[I-] NLKNQRATVPKPDG-UHFFFAOYSA-M 0.000 description 3
- 238000004040 coloring Methods 0.000 description 2
- 238000007405 data analysis Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000003709 image segmentation Methods 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 238000002347 injection Methods 0.000 description 2
- 239000007924 injection Substances 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- 238000005065 mining Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 238000012800 visualization Methods 0.000 description 2
- 238000005481 NMR spectroscopy Methods 0.000 description 1
- 241001085205 Prenanthella exigua Species 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004624 confocal microscopy Methods 0.000 description 1
- 239000010779 crude oil Substances 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005530 etching Methods 0.000 description 1
- 238000012921 fluorescence analysis Methods 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 238000002955 isolation Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 229920000642 polymer Polymers 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
- 239000011206 ternary composite Substances 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
- 239000011800 void material Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T15/00—3D [Three Dimensional] image rendering
- G06T15/50—Lighting effects
- G06T15/80—Shading
- G06T15/87—Gouraud shading
-
- 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/11—Region-based segmentation
-
- 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
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Graphics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
The invention discloses a quantitative analysis method for conversion of microscopic residual oil types based on CT technology. The method comprises the following steps: collecting a core sample of a region to be detected, simulating a displacement process for the core sample, and collecting two-dimensional CT scanning images corresponding to the core sample in at least two different displacement states in the displacement process; respectively recombining two-dimensional CT scanning images under each displacement state to obtain an oil phase three-dimensional model under the displacement state; according to the three-dimensional shape factor G of the residual oil clusters in the oil phase three-dimensional model and the Euler number E of the pore space N For the residual oil clustersDividing microscopic occurrence form types of the equipment; and obtaining a microscopic residual oil type conversion three-dimensional model according to the conversion relation of the residual oil clusters. And (3) converting the microscopic residual oil type into a three-dimensional model, and calculating the conversion rate of the residual oil clusters and the actual utilization ratio. The invention has more accurate type division of the residual oil, more accurate quantitative analysis and is beneficial to improving the utilization degree of microscopic residual oil.
Description
Technical Field
The invention relates to the field of quantitative analysis of conversion of types of micro residual oil, in particular to a quantitative analysis method of conversion of types of micro residual oil based on a CT technology.
Background
In the current exploitation work, the utilization of the micro-residual oil is difficult, especially the utilization of the micro-residual oil of a medium-high permeability oil reservoir is difficult, and one of the main reasons is that the type of the residual oil is not accurately divided, and an accurate quantitative analysis method for converting the micro-residual oil type is lacking. There are two main methods for analyzing the microscopic residual oil at present. One is a direct method, i.e. a technology for representing the distribution rule of residual oil in different development stages based on imaging technology, such as microscopic glass etching model, fluorescence analysis method, nuclear magnetic resonance and the like. The other is an indirect method, namely, a mathematical simulation model is established, and the model is characterized by a computer. The method has the advantages of being free from interference of external conditions, such as geological complexity, experimental conditions, sample representativeness and the like, and being capable of presenting three-dimensional space residual oil evolution rules.
Disclosure of Invention
The inventor finds that in the prior art, the direct method has the defect that the existing state and distribution rule of the residual oil on the two-dimensional plane can only be simulated, and the crude oil utilization ratio and the residual oil distribution characteristic in the three-dimensional space cannot be presented. The indirect method has the defects that the seepage rule of the oil-water phase in the reservoir is complex, and the difficulty of accurately establishing a mathematical model is high. No matter how the direct method and the indirect method currently relate to the conversion of the types of the micro residual oil, the division accuracy of the types of the micro residual oil is insufficient, the method for definitely determining the conversion process between the micro residual oil of different types at different stages in the oil reservoir development process is lacking, the quantitative conversion analysis is accurate, the utilization capacity of different displacement mediums on various micro residual oils is difficult to clear, the quantitative analysis method for lack of accuracy on the occurrence state, the utilization degree and the type conversion of the micro residual oil at different mining stages (a water injection stage, a chemical driving stage and a subsequent water driving stage) is difficult to clear, and the problems of low utilization degree and the like on the micro residual oil are caused. In order to at least partially solve the technical problems existing in the prior art, the inventor makes the invention, and through a specific embodiment, a method and a device for quantitative analysis of conversion of micro residual oil types based on CT technology are provided.
In a first aspect, an embodiment of the present invention provides a method for generating a three-dimensional model for converting a micro residual oil type, including:
collecting a core sample of a region to be detected, simulating a displacement process for the core sample, and collecting two-dimensional CT scanning images corresponding to the core sample in at least two different displacement states in the displacement process;
respectively recombining two-dimensional CT scanning images under each displacement state to obtain an oil phase three-dimensional model under the displacement state;
aiming at the oil phase three-dimensional model in the displacement state, according to the three-dimensional shape factor G of the residual oil clusters and the Euler number E of the pore space in the oil phase three-dimensional model N Dividing microscopic occurrence form types of the residual oil clusters, and distinguishing and marking pixel points of the residual oil clusters under various microscopic occurrence form types;
according to the transformation relation between the residual oil clusters of the preset microscopic occurrence form types and the residual oil clusters of all the residual oil clusters under two different displacement states, determining the pixel points which belong to the transformation part, the net non-use part and the actual use part and correspond to the pixel points of the residual oil clusters of the preset microscopic occurrence form types, and performing distinguishing marking to obtain a three-dimensional transformation model of the microscopic residual oil types.
In some alternative embodiments, the displacement state includes: a pre-displacement state, an in-displacement state, and a post-displacement state; the pre-displacement state is a saturated oil state.
In some optional embodiments, after collecting two-dimensional CT scan images corresponding to the core sample in at least two different displacement states during the displacement process, the method further includes: preprocessing a two-dimensional CT scanning image; the pretreatment comprises the following steps: segmentation processing, noise reduction and gray level binarization; the segmentation process refers to dividing different substances on a two-dimensional CT scanning image, firstly segmenting pores, then segmenting pores entering oil, and then segmenting pores swept by a displacement agent.
In some optional embodiments, for the two-dimensional CT scan image under each displacement state, oil-water segmentation is performed on the two-dimensional CT scan image, an oil-phase two-dimensional model is extracted from the two-dimensional CT scan image after oil-water segmentation, and according to the oil-phase two-dimensional model, an oil-phase three-dimensional model under the displacement state is recombined.
In some alternative embodiments, the oil-water separation includes: the pixel value intervals of the oil and the water are respectively divided by a threshold segmentation method.
In some alternative embodiments, the three-dimensional form factor G of the remaining oil clusters is calculated with the euler number E of the pore space N According to G and E N The microscopic occurrence forms of the residual oil clusters are divided into different types in different intervals, and pixel points of the residual oil clusters with various microscopic occurrence form types are marked and named.
In some optional embodiments, for a preset microscopic occurrence pattern type, the pixel points of the remaining oil clusters of the preset microscopic occurrence pattern type in the back-drive state intersect with the pixel points of the remaining oil clusters of all types except the preset microscopic occurrence pattern type in the front-drive state, so as to obtain pixel points corresponding to the remaining oil clusters of all types converted into the preset microscopic occurrence pattern type from the front-drive state to the back-drive state;
subtracting pixel points corresponding to the residual oil clusters of the preset microscopic occurrence form type from pixel points of the residual oil clusters of the preset microscopic occurrence form type in the back-driving state from pixel points corresponding to the residual oil clusters of all the residual types converted from the front-driving state to the back-driving state, so as to obtain pixel points corresponding to the net unused part of the residual oil clusters of the preset microscopic occurrence form type in the back-driving state;
Respectively intersecting the pixel points of the preset microcosmic occurrence form type residual oil clusters in the front displacement state with the pixel points of the residual oil clusters of all types remained in the rear displacement state to obtain pixel points corresponding to the residual oil clusters of all types converted from the preset microcosmic occurrence form type residual oil clusters in the front displacement state to the rear displacement state;
and subtracting the pixel points corresponding to the net unused part of the residual oil cluster of the preset microscopic occurrence form type in the back-driving state from the pixel points corresponding to the residual oil clusters of the preset microscopic occurrence form type in the front-driving state, and subtracting the pixel points corresponding to the actual unused part of the residual oil cluster of the preset microscopic occurrence form type converted into the residual oil clusters of all types in the back-driving state from the front-driving state to the back-driving state to obtain the pixel points corresponding to the actual unused part of the residual oil cluster of the preset microscopic occurrence form type in the front-driving state to the back-driving state.
In a second aspect, an embodiment of the present invention provides a quantitative analysis method for conversion of a micro-residual oil type, including:
Establishing a microscopic residual oil type conversion three-dimensional model;
the three-dimensional model for converting the micro residual oil type is obtained by a generation method of the three-dimensional model for converting the micro residual oil type;
and (3) converting the microscopic residual oil type into a three-dimensional model, and calculating the residual oil cluster conversion rate and the actual utilization ratio of the preset microscopic occurrence form type.
In some optional embodiments, counting the number of pixels corresponding to the conversion of the preset microscopic occurrence pattern type remaining oil clusters from the previous displacement state to the subsequent displacement state into the remaining all types of remaining oil clusters; counting the number of pixels of the residual oil clusters of the preset microscopic occurrence form type in a front displacement state; the number of pixels corresponding to the residual oil clusters of the preset microscopic occurrence form type is converted from the front displacement state to the rear displacement state, and the conversion rate of the residual oil clusters of the preset microscopic occurrence form type is obtained by dividing the number of pixels corresponding to the residual oil clusters of the preset microscopic occurrence form type in the front displacement state by the number of pixels corresponding to the residual oil clusters of the preset microscopic occurrence form type, wherein the conversion rate is the conversion rate of the residual oil clusters of the preset microscopic occurrence form type;
In some optional embodiments, counting the number of pixels corresponding to the actual usage portion of the remaining oil clusters of the preset microscopic occurrence pattern type from the previous displacement state to the subsequent displacement state; and dividing the number of pixels corresponding to the actual utilization part of the residual oil clusters of the preset microscopic occurrence form type from the front displacement state to the rear displacement state by the number of pixels of the residual oil clusters of the preset microscopic occurrence form type in the front displacement state to obtain the actual utilization ratio of the residual oil clusters of the preset microscopic occurrence form type in the front displacement state.
In a third aspect, an embodiment of the present invention provides a device for generating a three-dimensional model for converting a type of micro-residual oil, including:
the oil phase three-dimensional model building module is used for collecting a core sample of a region to be detected, simulating a displacement process for the core sample, and collecting two-dimensional CT scanning images corresponding to the core sample in two different displacement states in the displacement process; respectively recombining two-dimensional CT scanning images under each displacement state to obtain an oil phase three-dimensional model under the displacement state;
the residual oil cluster type marking module is used for aiming at the oil phase three-dimensional model in the displacement state, and according to the three-dimensional shape factor G of the residual oil clusters in the oil phase three-dimensional model and the Euler number E of the pore space N Dividing microscopic occurrence form types of the residual oil clusters, and distinguishing and marking pixel points of the residual oil clusters under various microscopic occurrence form types;
the transformation relation determining module is used for determining pixel points belonging to a transformation part, a net unused part and an actual unused part in the pixel points of the residual oil clusters of the preset microscopic occurrence form type according to transformation relations between the residual oil clusters of the preset microscopic occurrence form type and the residual oil clusters of all the residual oil types in two different displacement states and performing distinguishing marks to obtain a microscopic residual oil transformation three-dimensional model.
In a fourth aspect, an embodiment of the present invention provides a quantitative analysis device for conversion of a type of micro-residual oil, including:
the conversion three-dimensional model building module is used for building a conversion three-dimensional model of the type of the micro residual oil; the three-dimensional model for converting the micro residual oil type is obtained by a generation method of the three-dimensional model for converting the micro residual oil type;
the quantitative analysis module is used for converting the microscopic residual oil type into a three-dimensional model and calculating the residual oil cluster conversion rate and the actual utilization ratio of the preset microscopic occurrence form type.
In a fifth aspect, an embodiment of the present invention provides an electronic device, including: the system comprises a memory, a processor and a computer program stored on the memory and capable of running on the processor, wherein the processor realizes a generation method of a three-dimensional model for converting the type of the micro-residual oil or a quantitative analysis method for converting the type of the micro-residual oil when executing the program.
In a sixth aspect, an embodiment of the present invention provides a computer storage medium, where computer executable instructions are stored, where the computer executable instructions when executed by a processor implement a method for generating a three-dimensional model of conversion of the aforementioned micro-residual oil or a method for quantitative analysis of conversion of the aforementioned micro-residual oil type.
The technical scheme provided by the embodiment of the invention has the beneficial effects that at least:
the embodiment of the invention provides a generation method of a microscopic residual oil type conversion three-dimensional model, which comprises the steps of preprocessing a CT scanning image, including segmentation processing, noise reduction, gray level binarization and the like, so that different substances in the image are more accurately divided, the influence of non-characteristic factors such as image noise and the like is removed, the obtained image is more accurate, and the degree of automation of image processing is improved by combining imageJ with Avizo, so that the image processing is more accurate and rapid, and quantitative analysis based on the image is more accurate. Euler number E of three-dimensional shape factor G and pore space N By combining, the type of the micro-residual oil is divided according to the basis, so that the division of the micro-residual oil type is more accurate. Euler number E of three-dimensional shape factor G and pore space N The combination is used as the dividing basis, and is more reasonable and accurate than the independent dependence on the three-dimensional shape factor G as the dividing basis. According to the requirement of quantitative analysis of fineness, core samples can be acquired time-by-time for multiple times in the displacement process so as to analyze and obtain the type conversion condition and the use condition of the micro residual oil under the action of different types of displacement agents at different times in the displacement process, thereby reasonably adjusting the displacement agents on site, selecting the specific types of displacement agents at specific times and being beneficial to improving the use degree of the micro residual oil.
The second embodiment of the invention provides a quantitative analysis method for conversion of microscopic residual oil types, which is characterized in that a three-dimensional model is constructed, the areas where certain types of residual oil are located under different displacement states are marked, and the pixel points of the areas are calculated, so that quantitative analysis can be accurately and intuitively performed on the conversion process of the residual oil, various parameters such as the conversion rate of the residual oil, the actual utilization ratio and the like are conveniently calculated, reference is conveniently provided for on-site workers, and the displacement agent is reasonably adjusted, thereby improving the utilization degree of the residual oil. Because only the core sample is subjected to simulated displacement, the cost is low compared with the displacement experiment of the actual exploitation process. And because the two-dimensional image processing, three-dimensional image recombination, image segmentation, pixel point calculation and other processes are completed by the assistance of computer software, the invention has high degree of automation, is accurate, efficient, quick and visual, and further reduces the cost.
Additional features and advantages of the invention will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention. The objectives and other advantages of the invention will be realized and attained by the structure particularly pointed out in the written description and claims thereof as well as the appended drawings.
The technical scheme of the invention is further described in detail through the drawings and the embodiments.
Drawings
The accompanying drawings are included to provide a further understanding of the invention and are incorporated in and constitute a part of this specification, illustrate the invention and together with the embodiments of the invention, serve to explain the invention. In the drawings:
FIG. 1 is a flow chart of a method for generating a three-dimensional model of conversion of micro-residual oil types in an embodiment of the invention;
FIG. 2 is an original two-dimensional CT scan image in accordance with an embodiment of the present invention;
FIG. 3 is a pre-processed two-dimensional CT scan image according to an embodiment of the present invention;
FIG. 4 is a two-dimensional scan image of oil before water separation in accordance with an embodiment of the present invention;
FIG. 5 is a gray scale view of a two-dimensional scanned image after oil water separation in an embodiment of the present invention;
FIG. 6 is a three-dimensional image of a two-dimensional CT processed image reconstruction in accordance with an embodiment of the present invention;
FIG. 7 is a diagram of a microscopic occurrence pattern classification signature of the remaining oil in an embodiment of the present invention;
FIG. 8 is a three-dimensional image of a digital core marked with sets in accordance with one embodiment of the present invention;
FIG. 9 is a schematic view of the source of cluster residual oil in a water flooding process according to an embodiment of the present invention;
FIG. 10 is a flow chart of a quantitative analysis method for conversion of micro-residual oil types in an embodiment of the invention;
FIG. 11 is a schematic diagram of a three-dimensional model for converting a type of oil using micro-surplus according to an embodiment of the invention;
FIG. 12 is a block diagram of a generation apparatus for converting a micro-residual oil type into a three-dimensional model in accordance with an embodiment of the present invention;
FIG. 13 is a block diagram of a quantitative analysis device for conversion of micro-residual oil types in an embodiment of the invention.
Detailed Description
Exemplary embodiments of the present disclosure will be described in more detail below with reference to the accompanying drawings. While exemplary embodiments of the present disclosure are shown in the drawings, it should be understood that the present disclosure may be embodied in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the disclosure to those skilled in the art.
In order to solve the problems that in the prior art, the division accuracy of the types of the micro-residual oil is insufficient, a method for definitely determining the conversion process between the micro-residual oil of different types at different stages in the oil reservoir development process is insufficient, the quantitative conversion analysis is accurate, the utilization capability of different displacement mediums on various micro-residual oils is difficult to clear, and the quantitative analysis method for the occurrence state, the utilization degree and the type conversion of the micro-residual oil at different mining stages (a water injection stage, a chemical driving stage and a subsequent water driving stage) is lack, so that the utilization degree of the micro-residual oil is low and the like is caused.
Example 1
The first embodiment of the invention provides a method for generating a three-dimensional model of conversion of micro residual oil types, which is shown by referring to fig. 1, and comprises the following steps:
s1, collecting a core sample of a region to be detected, simulating a displacement process for the core sample, and collecting two-dimensional CT scanning images corresponding to the core sample in at least two different displacement states in the displacement process.
Because the digital core technology is adopted, only the core sample is subjected to simulated displacement, so the cost is very low compared with the displacement experiment in the actual exploitation process.
Wherein the displacement state comprises: a pre-displacement state, an in-displacement state, and a post-displacement state; the pre-displacement state is a saturated oil state.
The state during the displacement may include a plurality of states depending on one or more parameters, for example, a state after displacing for N minutes or hours, a state after displacing for 2N minutes or hours, a state after displacing for n+2 minutes or hours, or the like, and may be divided into a state of displacing for M milliliters of oil depending on the volume of the displaced oil, a state of displacing for 2 milliliters of oil, or the like, and a state of displacing for N minutes or hours after displacing for M milliliters of oil, or the like, depending on the displacement time parameter. Wherein the end of the displacement means that the displacement is ended when the water content of the produced liquid reaches more than 95%.
The displacement mode comprises water flooding, polymer flooding, ternary composite flooding and the like.
The core sample is shot or scanned by using an experimental instrument (such as a high-power optical microscope or an X-ray CT scanner) to obtain a large number of core two-dimensional CT scanning images, and then the two-dimensional images are overlapped and heavy to form a three-dimensional digital core through a modeling program or software, wherein the three-dimensional digital core mainly comprises a sequential slice imaging method, a laser scanning confocal microscopy method and an X-ray CT scanning method. Referring to fig. 2, fig. 2 is an original two-dimensional CT scan image.
And performing image preprocessing on the two-dimensional CT scanning image. The preprocessing comprises segmentation processing, noise reduction, gray level binarization and the like.
The segmentation process refers to dividing different substances on a two-dimensional CT scanning image, firstly segmenting pores, then segmenting pores entering oil, and then segmenting pores swept by a displacement agent. The rock core projection image directly obtained by CT scanning expresses the change of the density of the rock core composition substances through the change of gray values, and the pores in the rock core can not be directly distinguished from the framework. Therefore, it is necessary to divide the different components in a two-dimensional CT scan image for further utilization of the image. The gray value change regions of two different substances have crossing regions, the pore positions of which can be identified by naked eyes, but some boundary regions cannot be divided by naked eyes, for example, the boundary portions of the pores and the framework are fuzzy and difficult to distinguish, so that the image needs to be divided so as to reasonably divide the pores and the framework. For example, the liquid in the core comprises water and oil, after potassium iodide is added into the water, the water-containing part presents bright white, so that the water can be conveniently identified, if a certain pore is smaller, the water entering the pore is less, the brightness of the pore is darker, if a certain pore is larger and the oil is more, the gray value of the image is larger due to the large density of the oil, and the color of the image is closer to black.
Image noise reduction is the process of removing pixel noise from an image. There are a number of noise reduction algorithms available. Through proper noise reduction processing, the integrity (i.e. main characteristics) of the original information can be maintained as much as possible, and useless information in the image can be removed.
Gray level binarization refers to that a gray level image is selected through a proper threshold value to obtain a binarized image which can still reflect the whole or partial characteristics of the image, namely, the whole image is obviously black-white, so that different substances on the image are more obviously distinguished.
In one embodiment, imageJ software may be optionally used to perform segmentation assistance, noise reduction, gray scale binarization, and the like on the two-dimensional CT scan image. Besides basic image operation such as coloring, gray level adjustment and the like, the ImageJ software can also create statistical information according to user-defined parameters, such as image region and pixel statistics, spacing and angle calculation, can create a histogram and a section view, and perform fourier transformation and the like. Therefore, the two-dimensional CT scanning image can be processed more rapidly and accurately through the assistance of imageJ software. Referring to fig. 3, fig. 3 is a preprocessed two-dimensional CT scan image.
And S2, respectively recombining two-dimensional CT scanning images in each displacement state to obtain an oil phase three-dimensional model in the displacement state. The detailed steps are as follows:
and respectively carrying out oil-water segmentation on the two-dimensional CT scanning images in each displacement state, extracting an oil-phase two-dimensional model from the oil-water segmented two-dimensional CT scanning images, and recombining an oil-phase three-dimensional model in the displacement state according to the oil-phase two-dimensional model.
The oil-water separation refers to the separation of pixel value intervals of oil and water by a threshold separation method, so as to achieve the purpose of oil-water separation, wherein the separation effect can be more obvious through coloring treatment, for example, referring to fig. 4 and 5, fig. 4 is a two-dimensional scanning image before oil-water separation, fig. 5 is a gray scale image of a two-dimensional scanning image after oil-water separation, in a colorful two-dimensional scanning image after oil-water separation, an image area where oil is located is set to be red, an image area where water is located is set to be blue, and an image area where a skeleton is located is set to be black.
The recombination means to recombine the two-dimensional image after the segmentation processing into a three-dimensional image. For example, in this embodiment, 1440 processed two-dimensional CT images are reconstructed into a three-dimensional image, and the three-dimensional image shown in fig. 6 is referred to, so as to achieve the purpose of extracting three-dimensional pores or oil-water distribution.
In the process of segmentation and recombination, imageJ and Avizo are combined for use, so that the degree of automation of image processing is improved, and the image processing is more accurate and rapid. The Avizo is software for visualization and analysis of scientific and industrial data, and has the functions of three-dimensional reconstruction, three-dimensional image data rendering, visualization of flow simulation results inside a three-dimensional model through an advanced vector field, calculation and quantification of data such as density, distance, area, volume and the like through a statistics module, different marks made on the basis of individual pixel allocation to distinguish different structures to generate the three-dimensional model, further data analysis and the like. Therefore, through Avizo software, the three-dimensional recombination can be realized, the conversion process of the micro residual oil can be displayed, different substances such as pores, frameworks and different types of micro residual oil can be marked, and the calculation and quantitative analysis of data such as different substance densities, distances, areas, volumes, custom units and the like of the pores, the frameworks and the different types of micro residual oil can be performed through a statistics module. Therefore, with the aid of Avizo software, the degree of automation of image processing and data analysis is improved, and quantitative analysis on the conversion of the micro residual oil is more accurate and rapid.
The accuracy of oil-water segmentation is directly related to image preprocessing, wherein the pixel noise has a great influence on oil-water segmentation.
Extracting an oil phase two-dimensional model by using the two-dimensional image after oil-water segmentation, recombining the oil phase two-dimensional models corresponding to different displacement states into an oil phase three-dimensional model, and according to the three-dimensional shape factor G and the Euler number E of the pore space N And respectively dividing types of microscopic occurrence forms of the oil in the oil phase three-dimensional model, marking different colors for pixels of the oil with different microscopic occurrence form types, and marking the pixels of the oil with the same microscopic occurrence form type with the same color.
And (3) extracting an oil phase two-dimensional image from the two-dimensional CT processed image after oil-water segmentation under different displacement states, and establishing an oil phase two-dimensional model. For example, the portions other than the oil phase in the two-dimensional image are set to be colorless or white, the obtained image is an extracted oil phase two-dimensional image, and an oil phase two-dimensional model is created from these extracted oil phase two-dimensional images. And then reconstructing an oil phase three-dimensional image according to the oil phase two-dimensional model. Referring to fig. 5, fig. 5 is a three-dimensional image of an oil phase reconstructed from two-dimensional images of the oil phase.
Step S3, aiming at the oil phase three-dimensional model in the displacement state, according to the three-dimensional shape factor G of the residual oil clusters and the Euler number E of the pore space in the oil phase three-dimensional model N Dividing microscopic occurrence form types of the residual oil clusters, and distinguishing and marking pixel points of the residual oil clusters under various microscopic occurrence form types;
the residual oil enrichment mode on the pore scale is various and is influenced by the microscopic reservoir rock pore space structure type. Euler number E combining three-dimensional form factor G with pore space N The microscopic occurrence forms of the residual oil can be quantitatively distinguished. The three-dimensional shape factor represents the similarity degree of oil drops and spheres, the larger the shape factor is, the higher the similarity degree is, and for one sphere, the shape factor reaches the maximum value of 1, and the calculation formula is as follows:
wherein G is a three-dimensional form factor; pi is a circumference rate constant, and V is the volume of the residual oil clusters in cubic meters; s is the surface area of the remaining oil clusters, in square meters. The remaining oil clusters are a spatially aggregated morphology of the remaining oil.
The image topology parameter is a generic term for various feature parameters used for shape matching and object recognition in the field of image processing. Euler number (Euler number) is a space topology description mode, also called Euler feature number (Euler-Poincare Characteristic), and is divided into a three-dimensional surface and a three-dimensional body, and can be kept unchanged after operations such as image translation and rotation, so that the Euler number (Euler number) is widely applied to image processing aspects such as geological sandstone analysis and shadow detection, and the formula is as follows:
E N ＝b 0 -b 1 +b 2
Wherein E is N Euler number as pore space; b 0 The number of the connectors is the number; b 1 The number of holes refers to the maximum number of cuts without breaking, also called tunnel; b 2 Representing the number of holes, also known as cavity.
Calculating the Euler number E of the three-dimensional shape factor G and the pore space of the residual oil clusters N According to G and E N The microscopic occurrence forms of the residual oil clusters are divided into different types in different intervals, and pixel points of the residual oil clusters with various microscopic occurrence form types are marked and named.
Euler number E combining three-dimensional form factor G with void space in this patent N The microscopic appearance of the remaining oil can be quantitatively divided into five types including cluster, column, porous, membranous and isolated. Wherein, the cluster judgment basis is G < 0.1; the columnar judgment basis is that G is more than or equal to 0.1 and less than 0.3; the porous judgment basis is that G is more than or equal to 0.3 and less than 0.7 and E N < 1; the film-shaped judgment basis is G which is more than or equal to 0.3<0.7 and E N 1 or more; the determination basis of the isolation state is G is more than or equal to 0.7.
Euler number E of three-dimensional shape factor G and pore space N By combining, as a basis, the microscopic occurrence form of the oil in the three-dimensional model is divided into types, namely the microscopic residual oil is divided into types, so that the division of the microscopic residual oil types is more accurate. Euler number E of three-dimensional shape factor G and pore space N The combination is used as the dividing basis, and is more reasonable and accurate than the independent dependence on the three-dimensional shape factor G as the dividing basis. And marking different colors on the pixels of the oil with different microscopic occurrence forms, and marking the pixels of the oil with the same microscopic occurrence form with the same color.
For example, referring to fig. 7, a basis for quantitatively distinguishing the microscopic occurrence forms of the residual oil is input into the Avizo software or other software capable of realizing the similar functions, the Avizo classifies the residual oil according to the rule, marks different colors for distinguishing, and in the color version of fig. 7, different types of residual oil clusters are in different colors. In fig. 7, the darker colored portion represents that this type of micro-residual oil is denser in this portion.
And S4, determining pixel points which belong to a conversion part, a net non-use part and an actual use part and correspond to the pixel points in the residual oil cluster pixel points of the preset microscopic occurrence form type according to the conversion relation between the residual oil clusters of the preset microscopic occurrence form type and the residual oil clusters of all the residual oil types in two different displacement states, and performing distinguishing marking to obtain a microscopic residual oil type conversion three-dimensional model.
The process for determining the pixel points corresponding to the conversion part, the net unused part and the actual unused part in the residual oil cluster pixel points of the preset microscopic occurrence form type comprises the following steps of;
For a preset microscopic occurrence form type, intersecting pixel points of the residual oil clusters of the preset microscopic occurrence form type in a back-drive state with pixel points of the residual oil clusters of all types except the preset microscopic occurrence form type in a front-drive state to obtain pixel points corresponding to the residual oil clusters of all types converted into the preset microscopic occurrence form type from the front-drive state to the back-drive state;
subtracting pixel points corresponding to the residual oil clusters of the preset microscopic occurrence form type from pixel points of the residual oil clusters of the preset microscopic occurrence form type in the back-driving state from pixel points corresponding to the residual oil clusters of all the residual types converted from the front-driving state to the back-driving state, so as to obtain pixel points corresponding to the net unused part of the residual oil clusters of the preset microscopic occurrence form type in the back-driving state;
respectively intersecting the pixel points of the preset microcosmic occurrence form type residual oil clusters in the front displacement state with the pixel points of the residual oil clusters of all types remained in the rear displacement state to obtain pixel points corresponding to the residual oil clusters of all types converted from the preset microcosmic occurrence form type residual oil clusters in the front displacement state to the rear displacement state;
And subtracting the pixel points corresponding to the net unused part of the residual oil cluster of the preset microscopic occurrence form type in the back-driving state from the pixel points of the residual oil cluster of the preset microscopic occurrence form type in the front-driving state, and subtracting the pixel points corresponding to the residual oil clusters of all the remaining types from the residual oil cluster of the preset microscopic occurrence form type in the front-driving state to the residual oil cluster of all the remaining types in the back-driving state, so as to obtain the pixel points corresponding to the actual used part of the residual oil cluster of the preset microscopic occurrence form type from the residual oil cluster of the preset microscopic occurrence form type in the front-driving state to the residual oil cluster of the preset microscopic occurrence form type in the back-driving state.
In order to make the above-mentioned process of determining the pixel points corresponding to the conversion portion, the net unused portion, the actual used portion, etc. clearer, the preset microscopic occurrence form type is exemplified by cluster, and the displacement mode is exemplified by water driving. On the premise of defining the mutual conversion of the cluster residual oil types, firstly, the direction of the cluster residual oil in different stages (namely, different displacement states) comprises a mobilized part and an unused part, wherein the mobilized part comprises a produced part and a type conversion (cluster conversion into other types) part, and the unused part comprises a net unused part and a type conversion (other types are converted into cluster) part. The cluster residual oil can be divided into a direction and a source direction in the displacement process. In this direction, the cluster-shaped residual oil can be divided into the extracted residual oil and the residual oil converted from other types, and in this direction, the cluster-shaped residual oil can be divided into the residual oil which is not used and the residual oil converted from other types.
For the unused part of the cluster-shaped residual oil, the cluster-shaped residual oil in the water driving state (which can be the state in the process of driving or the state after the completion of driving) is intersected with the rest four types of residual oil except the cluster-shaped residual oil in the saturated oil state, the cluster-shaped residual oil is converted from the rest four types of residual oil in the water driving process, as shown in the following formula,
a (driven non-cluster) =a (driven cluster) Σa (saturated non-cluster)
The cluster-driven oil is a cluster-shaped residual oil set under a water-driven state (can be a middle state of displacement or after the end of displacement), wherein A (saturated non-cluster) is a cluster-shaped residual oil set which is formed by converting the four types of residual oil into the cluster-shaped residual oil set under the water-driven state, and the cluster-shaped residual oil set is formed by converting the four types of residual oil into the cluster-shaped residual oil set under the saturated oil state.
The cluster residual oil obtained by making a difference between the cluster residual oil and the part of the residual oil in the water driving state is a net unused part of the actual cluster residual oil, as shown in the following formula,
a (driving net motionless) =a (driving cluster) -a (driving non-cluster)
Wherein, A (driving cluster) is a cluster-shaped residual oil set under a water driving state (can be a middle driving state or can be after the driving is finished), A (driving non-cluster-rotating cluster) is a cluster-shaped residual oil set converted from the rest four types of residual oil under the water driving state, and A (driving cluster is clean) is a set of clean and unused parts of actual cluster-shaped residual oil under the water driving state.
For the cluster-like surplus oil-utilizing portion, it is possible to subdivide into an actual-utilizing portion and a surplus oil-converting portion. And (3) respectively carrying out intersection operation on the cluster-shaped residual oil in the saturated oil state and four types of residual oil except for the cluster-shaped residual oil in the water flooding state, wherein the obtained residual oil is a part converted from the cluster-shaped residual oil into the four types of residual oil in the water flooding process, namely a cluster-shaped residual oil conversion part. As shown in the following formula,
a (driven to non-cluster) =a (saturated cluster) Σa (driven to non-cluster)
Wherein, A (cluster saturation) is the aggregation of cluster-shaped residual oil in a saturated oil state, A (cluster non-driving) is the aggregation of four types of residual oil except for the cluster-shaped residual oil in a water driving state, and A (cluster driving to non-cluster driving) is the aggregation of the parts of the four types of residual oil converted from the cluster-shaped residual oil into the rest in the water driving state.
Subtracting the net unused part of the cluster-shaped residual oil and the cluster-shaped residual oil conversion part of the cluster-shaped residual oil in the water flooding state from the cluster-shaped residual oil in the saturated oil state, namely the actual using capacity of water flooding development on the residual oil, namely the actual using part of the cluster-shaped residual oil; as shown in the following formula,
a (cluster actual) =a (saturated cluster) -a (drive net actual) -a (drive to non-cluster)
Wherein, A (cluster actual) is the set of the actual use part of the cluster residual oil in the water driving state, A (saturated cluster) is the set of the cluster residual oil in the saturated oil state, A (driving cluster net non-movement) is the set of the net non-use part of the actual cluster residual oil in the water driving state, and A (driving cluster to non-cluster) is the set of the rest four types of residual oil parts converted from the cluster residual oil in the water driving state.
And marking the sets respectively to obtain the micro residual oil type conversion three-dimensional model. Referring to fig. 8, fig. 8 is a three-dimensional image of the digital core after marking each set. For example, in the Avizo software or other software capable of realizing similar functions, the above parts are subjected to set operation, and different colors of pixel marks corresponding to each set are distinguished, so as to obtain a micro residual oil type conversion three-dimensional model. Referring to fig. 9, a, b, and c show the directions of the cluster residual oil in the water driving process, and d, e, and f show the sources of the cluster residual oil in the water driving process. In the color version of fig. 9, the remaining oils of each type appear different colors.
According to the requirement of quantitative analysis of fineness, core samples can be acquired time-by-time for multiple times in the displacement process so as to analyze and obtain the type conversion condition and the use condition of the micro residual oil under the action of different types of displacement agents at different times in the displacement process, thereby reasonably adjusting the displacement agents on site, selecting the specific types of displacement agents at specific times and being beneficial to improving the use degree of the micro residual oil.
Example two
The second embodiment of the invention provides a quantitative analysis method for conversion of micro residual oil types, which is shown by referring to fig. 10, and comprises the following steps:
S5, establishing a microscopic residual oil type conversion three-dimensional model;
for a specific step of creating a three-dimensional model of the conversion of the micro-residual oil type, reference is made to embodiment one.
And S6, converting the three-dimensional model by using the microscopic residual oil type, and calculating the residual oil cluster conversion rate and the actual utilization ratio of the preset microscopic occurrence form type. Referring to fig. 11, a process in which the remaining oil cluster conversion rate of the preset microscopic occurrence morphology type is calculated is as follows:
counting the number of pixel points corresponding to the conversion of the preset microscopic occurrence form type residual oil clusters from the front displacement state to the rear displacement state into the residual oil clusters of all the types;
counting the number of pixels of the residual oil clusters of the preset microscopic occurrence form type in a front displacement state;
the number of pixels corresponding to the residual oil clusters of the preset microscopic occurrence form type is converted from the front displacement state to the rear displacement state, and the conversion rate of the residual oil clusters of the preset microscopic occurrence form type is obtained by dividing the number of pixels corresponding to the residual oil clusters of the preset microscopic occurrence form type in the front displacement state by the number of pixels corresponding to the residual oil clusters of the preset microscopic occurrence form type, wherein the conversion rate is the conversion rate of the residual oil clusters of the preset microscopic occurrence form type;
Taking the cluster residual oil conversion rate as an example, the calculation formula of the cluster residual oil conversion rate is as follows,
the actual utilization ratio of the residual oil clusters of the preset microscopic occurrence form type is calculated as follows:
counting the number of pixel points corresponding to the actual use part of the residual oil clusters of the preset microscopic occurrence form type from the front displacement state to the rear displacement state;
and dividing the number of pixels corresponding to the actual utilization part of the residual oil clusters of the preset microscopic occurrence form type from the front displacement state to the rear displacement state by the number of pixels of the residual oil clusters of the preset microscopic occurrence form type in the front displacement state to obtain the actual utilization ratio of the residual oil clusters of the preset microscopic occurrence form type in the front displacement state.
The calculation formula of the actual utilization ratio of the cluster residual oil is as follows,
similarly, other parameters may be designed and the result calculated based on the number of pixels in the set marked in the previous step. For example, a class of microscopic occurrence pattern residual oil extraction portions is calculated, taking a cluster-like residual oil extraction portion as an example, a calculation formula is,
a (cluster drive) =a (saturated cluster) -a (net drive) a (cluster drive to non-cluster drive)
Wherein, A (cluster driving) is the collection of the extraction parts of cluster residual oil in the water driving state, and other conforming meanings are the same as the previous description.
By using the microscopic residual oil conversion three-dimensional model, the areas where certain types of residual oil are located under different displacement states are found, and the pixel points of the areas are calculated, so that quantitative analysis can be accurately and intuitively carried out on the residual oil conversion process, various parameters such as residual oil conversion rate, actual utilization ratio and the like are conveniently calculated, reference is conveniently provided for on-site workers, and the utilization degree of the residual oil is improved. Because only the core sample is subjected to simulated displacement, the cost is low compared with the displacement experiment of the actual exploitation process. And because the two-dimensional image processing, three-dimensional image recombination, image segmentation, pixel point calculation and other processes are completed by the assistance of computer software, the invention has high degree of automation, is accurate, efficient, quick and visual, and further reduces the cost.
Example III
The third embodiment of the invention provides a device for generating a three-dimensional model for converting a micro residual oil type, which has a flow shown in fig. 12 and comprises:
the oil phase three-dimensional model building module 101 is used for collecting a core sample of a region to be detected, simulating a displacement process for the core sample, and collecting two-dimensional CT scanning images corresponding to the core sample in two different displacement states in the displacement process; respectively recombining two-dimensional CT scanning images under each displacement state to obtain an oil phase three-dimensional model under the displacement state;
The remaining oil cluster type marking module 102 is configured to, for the three-dimensional oil phase model in the displacement state, determine, according to the three-dimensional shape factor G of the remaining oil clusters in the three-dimensional oil phase model and the euler number E of the pore space N Dividing microscopic occurrence form types of the residual oil clusters, and distinguishing and marking pixel points of the residual oil clusters under various microscopic occurrence form types;
the conversion relation determining module 103 is configured to determine, according to the conversion relation between the remaining oil clusters of the preset microscopic occurrence form type and the remaining oil clusters of all types in the two different displacement states, pixels belonging to the conversion part, the net non-use part, and the pixel corresponding to the actual use part in the pixel points of the remaining oil clusters of the preset microscopic occurrence form type, and perform differential marking to obtain a microscopic remaining oil conversion three-dimensional model.
Example IV
The fourth embodiment of the invention provides a quantitative analysis device for conversion of micro residual oil types, the flow of which is shown in fig. 13, comprising:
the conversion three-dimensional model building module 111 is used for building a conversion three-dimensional model of the micro residual oil type; the three-dimensional model for converting the micro residual oil type is obtained by a generation method of the three-dimensional model for converting the micro residual oil type;
The quantitative analysis module 112 is configured to apply the microscopic remaining oil type conversion three-dimensional model, and calculate a remaining oil cluster conversion rate and an actual usage duty ratio of a preset microscopic occurrence form type.
An embodiment of the present invention provides an electronic device, including: the system comprises a memory, a processor and a computer program stored on the memory and capable of running on the processor, wherein the processor realizes a generation method of a three-dimensional model for converting the type of the micro-residual oil or a quantitative analysis method for converting the type of the micro-residual oil when executing the program.
The embodiment of the invention provides a computer storage medium, which is characterized in that computer executable instructions are stored in the computer storage medium, and the computer executable instructions realize a generation method of a three-dimensional model for converting micro residual oil or a quantitative analysis method for converting the micro residual oil type when being executed by a processor.
Unless specifically stated otherwise, terms such as processing, computing, calculating, determining, displaying, or the like, may refer to an action and/or process of one or more processing or computing systems, or similar devices, that manipulates and transforms data represented as physical (e.g., electronic) quantities within the processing system's registers or memories into other data similarly represented as physical quantities within the processing system's memories, registers or other such information storage, transmission or display devices. Information and signals may be represented using any of a variety of different technologies and techniques. For example, data, instructions, commands, information, signals, bits, symbols, and chips that may be referenced throughout the above description may be represented by voltages, currents, electromagnetic waves, magnetic fields or particles, optical fields or particles, or any combination thereof.
It should be understood that the specific order or hierarchy of steps in the processes disclosed are examples of exemplary approaches. Based on design preferences, it is understood that the specific order or hierarchy of steps in the processes may be rearranged without departing from the scope of the present disclosure. The accompanying method claims present elements of the various steps in a sample order, and are not meant to be limited to the specific order or hierarchy presented.
In the foregoing detailed description, various features are grouped together in a single embodiment for the purpose of streamlining the disclosure. This method of disclosure is not to be interpreted as reflecting an intention that the claimed embodiments of the subject matter require more features than are expressly recited in each claim. Rather, as the following claims reflect, invention lies in less than all features of a single disclosed embodiment. Thus the following claims are hereby expressly incorporated into this detailed description, with each claim standing on its own as a separate preferred embodiment of this invention.
Those of skill would further appreciate that the various illustrative logical blocks, modules, circuits, and algorithm steps described in connection with the embodiments disclosed herein may be implemented as electronic hardware, computer software, or combinations of both. To clearly illustrate this interchangeability of hardware and software, various illustrative components, blocks, modules, circuits, and steps have been described above generally in terms of their functionality. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the overall system. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present disclosure.
The steps of a method or algorithm described in connection with the embodiments disclosed herein may be embodied directly in hardware, in a software module executed by a processor, or in a combination of the two. A software module may reside in RAM memory, flash memory, ROM memory, EPROM memory, EEPROM memory, registers, hard disk, a removable disk, a CD-ROM, or any other form of storage medium known in the art. An exemplary storage medium is coupled to the processor such the processor can read information from, and write information to, the storage medium. In the alternative, the storage medium may be integral to the processor. The processor and the storage medium may reside in an ASIC. The ASIC may reside in a user terminal. The processor and the storage medium may reside as discrete components in a user terminal.
For a software implementation, the techniques described herein may be implemented with modules (e.g., procedures, functions, and so on) that perform the functions described herein. These software codes may be stored in memory units and executed by processors. The memory unit may be implemented within the processor or external to the processor, in which case it can be communicatively coupled to the processor via various means as is known in the art.
The foregoing description includes examples of one or more embodiments. It is, of course, not possible to describe every conceivable combination of components or methodologies for purposes of describing the aforementioned embodiments, but one of ordinary skill in the art may recognize that many further combinations and permutations of various embodiments are possible. Accordingly, the embodiments described herein are intended to embrace all such alterations, modifications and variations that fall within the scope of the appended claims. Furthermore, as used in the specification or claims, the term "comprising" is intended to be inclusive in a manner similar to the term "comprising," as interpreted when employed as a transitional word in a claim. Furthermore, any use of the term "or" in the specification of the claims is intended to mean "non-exclusive or".
Claims (14)
1. The method for generating the three-dimensional model by converting the type of the micro residual oil is characterized by comprising the following steps of:
collecting a core sample of a region to be detected, simulating a displacement process for the core sample, and collecting two-dimensional CT scanning images corresponding to the core sample in at least two different displacement states in the displacement process;
respectively recombining two-dimensional CT scanning images under each displacement state to obtain an oil phase three-dimensional model under the displacement state;
Aiming at the oil phase three-dimensional model in the displacement state, according to the three-dimensional shape factor G of the residual oil clusters and the Euler number E of the pore space in the oil phase three-dimensional model N Dividing microscopic occurrence form types of the residual oil clusters, and distinguishing and marking pixel points of the residual oil clusters under various microscopic occurrence form types;
according to the transformation relation between the residual oil clusters of the preset microscopic occurrence form types and the residual oil clusters of all the residual oil clusters under two different displacement states, determining the pixel points which belong to the transformation part, the net non-use part and the actual use part and correspond to the pixel points of the residual oil clusters of the preset microscopic occurrence form types, and performing distinguishing marking to obtain a three-dimensional transformation model of the microscopic residual oil types.
2. The method of claim 1, wherein the displacement state comprises:
a pre-displacement state, an in-displacement state, and a post-displacement state; the pre-displacement state is a saturated oil state.
3. The method of claim 1, wherein after collecting two-dimensional CT scan images corresponding to the core sample in at least two different displacement states during the displacement, further comprising:
preprocessing a two-dimensional CT scanning image;
The pretreatment comprises the following steps: segmentation processing, noise reduction and gray level binarization;
the segmentation process refers to dividing different substances on a two-dimensional CT scanning image, firstly segmenting pores, then segmenting pores entering oil, and then segmenting pores swept by a displacement agent.
4. The method of claim 1, wherein the reconstructing the two-dimensional CT scan image for each displacement state to obtain the three-dimensional model of the oil phase for the displacement state comprises:
and respectively carrying out oil-water segmentation on the two-dimensional CT scanning images in each displacement state, extracting an oil-phase two-dimensional model from the oil-water segmented two-dimensional CT scanning images, and recombining an oil-phase three-dimensional model in the displacement state according to the oil-phase two-dimensional model.
5. The method of claim 4, wherein the oil-water separation comprises:
the pixel value intervals of the oil and the water are respectively divided by a threshold segmentation method.
6. The method of claim 1, wherein the three-dimensional form factor G and euler number E of pore space are based on the remaining oil clusters in the three-dimensional model of the oil phase N Dividing microscopic occurrence form types of the residual oil clusters, including:
Calculating the Euler number E of the three-dimensional shape factor G and the pore space of the residual oil clusters N According to G and E N Microcosmic of remaining oil clusters in different regionsThe occurrence forms are divided into different types, and pixel points of residual oil clusters of various microscopic occurrence form types are marked and named.
7. The method of claim 1, wherein determining, from the conversion relationships between the remaining oil clusters of the preset microscopic occurrence pattern type and the remaining oil clusters of all types in the two different displacement states, the pixel points corresponding to the conversion part, the net non-use part, and the actual use part among the pixel points of the remaining oil clusters of the preset microscopic occurrence pattern type includes:
for a preset microscopic occurrence form type, intersecting pixel points of the residual oil clusters of the preset microscopic occurrence form type in a back-drive state with pixel points of the residual oil clusters of all types except the preset microscopic occurrence form type in a front-drive state to obtain pixel points corresponding to the residual oil clusters of all types converted into the preset microscopic occurrence form type from the front-drive state to the back-drive state;
Subtracting pixel points corresponding to the residual oil clusters of the preset microscopic occurrence form type from pixel points of the residual oil clusters of the preset microscopic occurrence form type in the back-driving state from pixel points corresponding to the residual oil clusters of all the residual types converted from the front-driving state to the back-driving state, so as to obtain pixel points corresponding to the net unused part of the residual oil clusters of the preset microscopic occurrence form type in the back-driving state;
respectively intersecting the pixel points of the preset microcosmic occurrence form type residual oil clusters in the front displacement state with the pixel points of the residual oil clusters of all types remained in the rear displacement state to obtain pixel points corresponding to the residual oil clusters of all types converted from the preset microcosmic occurrence form type residual oil clusters in the front displacement state to the rear displacement state;
and subtracting the pixel points corresponding to the net unused part of the residual oil cluster of the preset microscopic occurrence form type in the back-driving state from the pixel points corresponding to the residual oil clusters of the preset microscopic occurrence form type in the front-driving state, and subtracting the pixel points corresponding to the actual unused part of the residual oil cluster of the preset microscopic occurrence form type converted into the residual oil clusters of all types in the back-driving state from the front-driving state to the back-driving state to obtain the pixel points corresponding to the actual unused part of the residual oil cluster of the preset microscopic occurrence form type in the front-driving state to the back-driving state.
8. A quantitative analysis method for conversion of a microscopic residual oil type, comprising:
establishing a microscopic residual oil type conversion three-dimensional model;
the three-dimensional model for converting the micro residual oil type is obtained by the generation method of the three-dimensional model for converting the micro residual oil type according to any one of claims 1 to 7;
and (3) converting the microscopic residual oil type into a three-dimensional model, and calculating the residual oil cluster conversion rate and the actual utilization ratio of the preset microscopic occurrence form type.
9. The method of claim 8, wherein the calculating the residual oil cluster conversion rate of the preset microscopic occurrence pattern type by using the microscopic residual oil type conversion three-dimensional model comprises:
counting the number of pixel points corresponding to the conversion of the preset microscopic occurrence form type residual oil clusters from the front displacement state to the rear displacement state into the residual oil clusters of all the types;
counting the number of pixels of the residual oil clusters of the preset microscopic occurrence form type in a front displacement state;
and dividing the number of pixels corresponding to the residual oil clusters of the preset microscopic occurrence form type from the front displacement state to the rear displacement state by the number of pixels corresponding to the residual oil clusters of all the remaining types, so as to obtain the conversion ratio of the residual oil clusters of the preset microscopic occurrence form type to the residual oil, namely the conversion rate of the residual oil clusters of the preset microscopic occurrence form type.
10. The method of claim 8, wherein the converting the microscopic residual oil type into the three-dimensional model, calculating the actual usage duty ratio of the residual oil clusters of the preset microscopic occurrence pattern type, comprises:
counting the number of pixel points corresponding to the actual use part of the residual oil clusters of the preset microscopic occurrence form type from the front displacement state to the rear displacement state;
and dividing the number of pixels corresponding to the actual utilization part of the residual oil clusters of the preset microscopic occurrence form type from the front displacement state to the rear displacement state by the number of pixels of the residual oil clusters of the preset microscopic occurrence form type in the front displacement state to obtain the actual utilization ratio of the residual oil clusters of the preset microscopic occurrence form type in the front displacement state.
11. A device for generating a three-dimensional model of conversion of a type of micro-residual oil, comprising:
the oil phase three-dimensional model building module is used for collecting a core sample of a region to be detected, simulating a displacement process for the core sample, and collecting two-dimensional CT scanning images corresponding to the core sample in two different displacement states in the displacement process; respectively recombining two-dimensional CT scanning images under each displacement state to obtain an oil phase three-dimensional model under the displacement state;
The residual oil cluster type marking module is used for aiming at the oil phase three-dimensional model in the displacement state, and according to the three-dimensional shape factor G of the residual oil clusters in the oil phase three-dimensional model and the Euler number E of the pore space N Dividing microscopic occurrence form types of the residual oil clusters, and distinguishing and marking pixel points of the residual oil clusters under various microscopic occurrence form types;
the transformation relation determining module is used for determining pixel points belonging to a transformation part, a net unused part and an actual unused part in the pixel points of the residual oil clusters of the preset microscopic occurrence form type according to transformation relations between the residual oil clusters of the preset microscopic occurrence form type and the residual oil clusters of all the residual oil types in two different displacement states and performing distinguishing marks to obtain a microscopic residual oil transformation three-dimensional model.
12. A microscopic remaining oil type conversion quantitative analysis device, characterized by comprising:
the conversion three-dimensional model building module is used for building a conversion three-dimensional model of the type of the micro residual oil; the three-dimensional model for converting the micro residual oil type is obtained by the generation method of the three-dimensional model for converting the micro residual oil type according to any one of claims 1 to 7;
the quantitative analysis module is used for converting the microscopic residual oil type into a three-dimensional model and calculating the residual oil cluster conversion rate and the actual utilization ratio of the preset microscopic occurrence form type.
13. An electronic device, comprising: a memory, a processor and a computer program stored on the memory and executable on the processor, the processor implementing the method for generating a three-dimensional model of conversion of a micro-residual oil type according to any one of claims 1 to 7 or the method for quantitative analysis of conversion of a micro-residual oil type according to any one of claims 8 to 10 when the program is executed.
14. A computer storage medium having stored therein computer executable instructions which when executed by a processor implement the method of generating a three-dimensional model of conversion of micro-residual oil according to any one of claims 1 to 7 or the quantitative analysis method of type conversion of micro-residual oil according to any one of claims 8 to 10.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111493785.7A CN116246007A (en) | 2021-12-08 | 2021-12-08 | Microscopic residual oil type conversion quantitative analysis method based on CT technology |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111493785.7A CN116246007A (en) | 2021-12-08 | 2021-12-08 | Microscopic residual oil type conversion quantitative analysis method based on CT technology |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN116246007A true CN116246007A (en) | 2023-06-09 |
Family
ID=86633663
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111493785.7A Pending CN116246007A (en) | 2021-12-08 | 2021-12-08 | Microscopic residual oil type conversion quantitative analysis method based on CT technology |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116246007A (en) |
-
2021
- 2021-12-08 CN CN202111493785.7A patent/CN116246007A/en active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106248559B (en) | A kind of five sorting technique of leucocyte based on deep learning | |
| CN104237103B (en) | A kind of pore communication quantitatively characterizing method and device | |
| Baveye et al. | Observer-dependent variability of the thresholding step in the quantitative analysis of soil images and X-ray microtomography data | |
| HORGAN | Mathematical morphology for analysing soil structure from images | |
| CN108876901A (en) | A kind of digital cores method for reconstructing based on two dimensional image and multi-point statistic | |
| Bauer et al. | From computed microtomography images to resistivity index calculations of heterogeneous carbonates using a dual-porosity pore-network approach: Influence of percolation on the electrical transport properties | |
| CN105445160B (en) | The void characteristics and its extracting method of a kind of asphalt | |
| CN109285222A (en) | The building of organic shale high-resolution digital rock core and analysis method | |
| Duan et al. | Image partitioning into convex polygons | |
| CN108734714B (en) | Method for analyzing carbonate rock structure based on Matlab | |
| CN112051200B (en) | Quantitative evaluation method and device for pore structure of tight sandstone reservoir | |
| CN109961446A (en) | CT/MR three-dimensional image segmentation processing method, device, equipment and medium | |
| Friese et al. | Analysis of tomographic mineralogical data using YaDiV—Overview and practical case study | |
| CN113222992B (en) | Crack characteristic characterization method and system based on multi-fractal spectrum | |
| CN111739149B (en) | Oil-water distribution continuity restoration method for rock CT scanning image | |
| CN112163621A (en) | Compact sandstone reservoir pore structure classification and characterization method based on micro ct technology | |
| CN115546605A (en) | Training method and device based on image labeling and segmentation model | |
| CN109472790A (en) | A kind of machine components defect inspection method and system | |
| CN112767266B (en) | Deep learning-oriented microscopic endoscope image data enhancement method | |
| Kang et al. | Construction of complex digital rock physics based on full convolution network | |
| CN113570652A (en) | Quantitative analysis method for sandstone reservoir mineral intercrystalline pores based on SEM image | |
| CN116246007A (en) | Microscopic residual oil type conversion quantitative analysis method based on CT technology | |
| CN110660051A (en) | Tensor voting processing method based on navigation pyramid | |
| CN113724223B (en) | YOLOv3 data set production method and system based on optical microscope | |
| CN115147539A (en) | Indication-kriging-based pore structure residual oil characterization method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |