CN110853138B - Construction method of dual-medium carbonate pore-crack dual-network model - Google Patents
Construction method of dual-medium carbonate pore-crack dual-network model Download PDFInfo
- Publication number
- CN110853138B CN110853138B CN201911150076.1A CN201911150076A CN110853138B CN 110853138 B CN110853138 B CN 110853138B CN 201911150076 A CN201911150076 A CN 201911150076A CN 110853138 B CN110853138 B CN 110853138B
- Authority
- CN
- China
- Prior art keywords
- network model
- crack
- pore
- dual
- space
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 title claims abstract description 39
- 238000010276 construction Methods 0.000 title claims abstract description 15
- 239000011148 porous material Substances 0.000 claims abstract description 74
- 238000000034 method Methods 0.000 claims abstract description 40
- 230000009977 dual effect Effects 0.000 claims abstract description 26
- 239000011435 rock Substances 0.000 claims abstract description 17
- 238000004458 analytical method Methods 0.000 claims abstract description 10
- 238000004891 communication Methods 0.000 claims abstract description 10
- 230000000877 morphologic effect Effects 0.000 claims abstract description 8
- 238000003709 image segmentation Methods 0.000 claims abstract description 5
- 238000000926 separation method Methods 0.000 claims abstract description 4
- 238000012545 processing Methods 0.000 claims abstract description 3
- 238000000605 extraction Methods 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 7
- 238000004422 calculation algorithm Methods 0.000 claims description 6
- 238000002591 computed tomography Methods 0.000 claims description 6
- 238000009826 distribution Methods 0.000 claims description 6
- 230000008602 contraction Effects 0.000 claims description 3
- 238000007781 pre-processing Methods 0.000 claims description 3
- 125000005587 carbonate group Chemical group 0.000 claims description 2
- 230000007480 spreading Effects 0.000 claims description 2
- 238000003892 spreading Methods 0.000 claims description 2
- 238000011161 development Methods 0.000 abstract description 7
- 230000018109 developmental process Effects 0.000 abstract description 7
- 239000012530 fluid Substances 0.000 description 8
- 239000011800 void material Substances 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000005457 optimization Methods 0.000 description 5
- 238000011160 research Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000012512 characterization method Methods 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 238000005325 percolation Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000010339 dilation Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 230000002035 prolonged effect Effects 0.000 description 1
- 238000007670 refining Methods 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
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
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
-
- G—PHYSICS
- 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/60—Analysis of geometric attributes
- G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2200/00—Indexing scheme for image data processing or generation, in general
- G06T2200/08—Indexing scheme for image data processing or generation, in general involving all processing steps from image acquisition to 3D model generation
-
- 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 Vision & Pattern Recognition (AREA)
- Chemical & Material Sciences (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- Computer Graphics (AREA)
- Dispersion Chemistry (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Geophysics And Detection Of Objects (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
The invention relates to the technical field of oil gas development, in particular to a method for constructing a dual-medium carbonate pore-crack dual-network model. The method comprises the following steps: scanning the dual-medium carbonate rock sample to obtain a three-dimensional image; performing image segmentation processing on the three-dimensional image to extract a space to be identified; performing morphological analysis on the space to be identified, and completing separation of a pore space and a crack space according to morphological analysis results; respectively constructing a network model of the separated pore space and the separated crack space by adopting a sphere filling method to obtain a pore network model and a crack network model; and adding a communication throat between the pore network model and the crack network model to complete the construction of the pore-crack dual network model. The invention can accurately construct the pore network model and the crack network model when constructing the dual-medium carbonate network model, and simultaneously combine the pore network model and the crack network model into the pore-crack dual-network model.
Description
Technical Field
The invention relates to the technical field of oil gas development, in particular to a method for constructing a dual-medium carbonate pore-crack dual-network model.
Background
With the continuous enhancement of the exploration and development of China, the types of oil and gas reservoirs encountered at present are more and more complex, and particularly the exploration reserves of pore-fracture type carbonate reservoirs are continuously increased. The pore-crack type carbonate gas reservoir is different from the conventional sandstone gas reservoir, and a crack and pore dual-medium system is formed in the pore-crack type carbonate gas reservoir due to the existence of a large number of natural cracks, so that the seepage characteristics and the fluid exchange relationship are more complex, and the difficulty in developing a development scheme is greater.
For the research of the seepage characteristics of carbonate rock with cracks and pores simultaneously developed, the common practice is to construct a rock mass system with the cracks and the common pore medium, and solve and calculate the flow characteristics of coupling fluid in different pore and crack mediums by using a finite element method. Because the operation process of the method is complex and a large enough model needs to be constructed for the core sample with strong heterogeneity, the research period is greatly prolonged, and the method cannot be popularized and used on a large scale.
The pore network model has become an important tool for the characterization of the microscopic pore structure and the prediction of the microscopic seepage characteristics of the core since the introduction of the core analysis field in the 50 th century of 20 th year. The traditional pore network model construction method is mainly used for porous media, and due to continuity of crack space and wide three-dimensional expansion, the conventional modeling method can lose part of information so that the true three-dimensional topological structure of the crack cannot be accurately captured. In addition, for dual-medium carbonate rock with pores and cracks which develop simultaneously, fluid exchange exists between the pores and the cracks, and a certain error obviously exists in the simultaneous characterization of the pores and the cracks by using a set of network models, so that the later research and analysis are influenced.
To sum up, when constructing a dual-medium carbonate network model, two problems need to be solved simultaneously: firstly, realizing the accurate construction of a fracture network model; secondly, after a network model is respectively constructed for two mediums of the pore crack, the two mediums are fused into a set of pore-crack dual network model.
Disclosure of Invention
Aiming at the defects existing in the prior art, the invention provides a dual-medium carbonate pore-crack dual-network model construction method, which can accurately construct a pore network model and a crack network model when constructing a dual-medium carbonate network model in application, and simultaneously fuse the pore network model and the crack network model into a pore-crack dual-network model.
The technical scheme adopted by the invention is as follows:
the construction method of the dual-medium carbonate pore-crack dual-network model comprises the following steps:
step one: collecting a dual-medium carbonate sample, and scanning the dual-medium carbonate sample to obtain a nondestructive three-dimensional image;
step two: performing image segmentation processing on the three-dimensional image, and extracting a space to be identified, which contains a pore space and a crack space, from the three-dimensional image;
step three: performing morphological analysis on the space to be identified, and completing separation of a pore space and a crack space according to morphological analysis results;
step four: respectively constructing a network model of the separated pore space and the separated crack space by adopting a sphere filling method to obtain a pore network model and a crack network model;
step five: and adding a communication throat between the pore network model and the crack network model to complete the construction of the pore-crack dual network model.
In the first step, a CT scanner is used to perform CT scanning on the dual-medium carbonate sample to obtain a nondestructive three-dimensional image.
As a preferred aspect of the above technical solution, the specific process for obtaining a three-dimensional image of a dual medium carbonate sample includes: firstly, scanning a dual-medium carbonate rock sample by using a CT scanner to obtain projection data; then, reconstructing the projection data into a core gray image by using an image reconstruction method; and finally, preprocessing the gray level image to obtain a three-dimensional image of the dual-medium carbonate rock sample.
As a preferred aspect of the above-mentioned aspect, in the third step, the process of performing morphological analysis on the void space and the crack space includes:
s31, measuring and obtaining three-dimensional length and three-dimensional minimum width data of a pore space and a crack space;
s32, calculating to obtain the ratio of the three-dimensional length to the three-dimensional minimum width;
s33, setting a threshold value, comparing the calculated ratio with the threshold value, and judging the crack space or the pore space according to the comparison result.
As a preferable mode of the above technical solution, in the fourth step, the sphere filling method includes: the three-dimensional space spreading and topological structure is depicted by adopting a central line extraction method, then an inscribed sphere is established by taking a central line voxel as a sphere center and taking a space boundary point as a tangent point, finally, a rod body is established to connect all continuous inscribed spheres in a centering manner, the opening degree of the space is represented by the interconnected inscribed spheres, and the connectivity and seepage channel information of the connected rod body represent the space.
As the optimization of the technical scheme, the specific steps of constructing the network model for the crack space are as follows:
s41, center line extraction: extracting continuous voxel chains, wherein the voxel chains are crack center lines of a crack space;
s42, determining an expandable space by using each voxel on a central line as a basic point and adopting an expansion algorithm, after finding the range of the double-medium carbonate skeleton voxels nearest to each voxel, identifying an inscribed sphere corresponding to the basic point by adopting a contraction algorithm, and finally calculating to obtain an upper limit and a lower limit of the inscribed sphere radius;
s43, after a series of inscribed spheres taking a central line voxel as a sphere center are extracted, optimizing the extracted inscribed spheres to obtain a residual inscribed sphere set;
and S43, sequentially centering and connecting the rest inscribed sphere sets through the rod bodies to form a crack space for a network model.
As a preferred embodiment of the present invention, in step S43, the process of optimizing the extracted inscribed sphere includes: firstly, according to the square size of the radius of each voxel corresponding to the inscribed sphere, each voxel is represented in a linked list form; and searching the linked list for the inscribed sphere corresponding to the voxel with the inclusion relationship, deleting the included inscribed sphere with the inclusion relationship, and finally leaving the rest inscribed sphere set which is not included.
As a preferred embodiment of the above technical solution, in the fifth step, the specific step of constructing the pore-crack dual network model includes:
s51, marking contact surfaces of a pore space and a crack space, finding out respective maximum inscribed spheres in a pore network model and a crack network model corresponding to the contact surfaces, and numbering the maximum inscribed spheres;
s52, obtaining a throat radius and coordination number distribution relation communicated between the pore network model and the crack network model according to the number and the area of the contact surfaces and the corresponding number of the maximum inscribed sphere;
and S53, connecting the largest inscribed spheres in the pore network model and the crack network model corresponding to the contact surface through the throat according to the throat radius and coordination number distribution relation to form a pore-crack dual network model.
The beneficial effects of the invention are as follows:
1. the invention can accurately construct the pore network model and the crack network model when constructing the dual-medium carbonate network model, and simultaneously combine the pore network model and the crack network model into a set of pore-crack dual-network model.
2. The invention utilizes a method of combining central line extraction with maximum inscribed sphere to establish a fracture network model for finely representing the real fracture space topological structure in the rock, and the model can be directly used for parameter representation and seepage characteristic simulation of the fracture space.
3. The pore-crack dual network model constructed by the invention can be regarded as a quasi-static model, the fluid flow in the model is controlled by the capillary pressure, and the problem of singular points caused by the flow in cracks when the physical properties of the fluid change can be ignored. By utilizing the invasion-percolation theory, the relative permeability prediction of each phase of fluid can be realized when multiphase fluid flows simultaneously during dual-medium carbonate rock seepage simulation and oilfield development.
4. The pore-crack dual network model constructed by the invention can be regarded as a quasi-static model, the calculation speed is greatly improved compared with the traditional method for solving and calculating the flow characteristics of coupling fluid in different media of pores and cracks by using a finite element method, the method can be directly applied to development and production of oil and gas fields, and mechanism recognition and data support are provided for the design of an oil and gas field development scheme.
Drawings
In order to more clearly illustrate the embodiments of the invention or the technical solutions in the prior art, the drawings that are required in the embodiments or the description of the prior art will be briefly described, it being obvious that the drawings in the following description are only some embodiments of the invention, and that other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a schematic flow chart of the steps of the present invention;
FIG. 2 is a micrometer CT scan image of dual media carbonate;
FIG. 3 is a diagram of a pore network model, a fracture network model and a constituent pore-fracture dual network model of dual medium carbonate;
FIG. 4 is a schematic drawing of extraction of void space and crack space from CT images;
FIG. 5 is a schematic view of a three-dimensional length search angle in example 2;
FIG. 6 is a schematic illustration of the void space and fracture space separated;
FIG. 7 is a schematic view of the centerline extraction in example 3;
FIG. 8 is a schematic diagram of the maximum inscribed sphere and percolation path partitioning in example 3 showing the crack space;
FIG. 9 is a schematic diagram of a method for characterizing a split network model in example 3;
FIG. 10 is a schematic illustration of the pore and crack interface markings of example 4;
FIG. 11 is a schematic diagram of the embodiment 4 before the throat is connected between the pores and the cracks;
FIG. 12 is a schematic representation of example 4 after the throat is connected between the pores and the cracks.
Detailed Description
The invention is further described with reference to the drawings and specific examples. It should be noted that the description of these examples is for aiding in understanding the present invention, but is not intended to limit the present invention. Specific structural and functional details disclosed herein are merely representative of example embodiments of the invention. This invention may, however, be embodied in many alternate forms and should not be construed as limited to the embodiments set forth herein.
It should be appreciated that the terms first, second, etc. are used merely to distinguish between descriptions and should not be construed as indicating or implying relative importance. Although the terms first, second, etc. may be used herein to describe various elements, these elements should not be limited by these terms. These terms are only used to distinguish one element from another element. For example, a first element could be termed a second element, and, similarly, a second element could be termed a first element, without departing from the scope of example embodiments of the present invention.
It should be understood that the term "and/or" is merely an association relationship describing the associated object, and means that three relationships may exist, for example, a and/or B may mean: the terms "/and" herein describe another associative object relationship, indicating that there may be two relationships, e.g., a/and B, may indicate that: the character "/" herein generally indicates that the associated object is an "or" relationship.
It should be understood that in the description of the present invention, the terms "upper", "vertical", "inner", "outer", etc. indicate an azimuth or a positional relationship in which the inventive product is conventionally put in use, or an azimuth or a positional relationship that are conventionally understood by those skilled in the art, are merely for convenience of describing the present invention and simplifying the description, and do not indicate or imply that the apparatus or elements referred to must have a specific azimuth, be configured and operated in a specific azimuth, and thus should not be construed as limiting the present invention.
It will be understood that when an element is referred to as being "connected," "connected," or "coupled" to another element, it can be directly connected or coupled to the other element or intervening elements may be present. In contrast, when an element is referred to as being "directly connected" or "directly coupled" to another element, there are no intervening elements present. Other words used to describe relationships between elements (e.g., "between … …" pair "directly between … …", "adjacent" pair "directly adjacent", etc.) should be interpreted in a similar manner.
In the description of the present invention, it should also be noted that, unless explicitly specified and limited otherwise, the terms "disposed," "mounted," and "connected" are to be construed broadly, and may be, for example, fixedly connected, detachably connected, or integrally connected; can be mechanically or electrically connected; can be directly connected or indirectly connected through an intermediate medium, and can be communication between two elements. The specific meaning of the above terms in the present invention will be understood in specific cases by those of ordinary skill in the art.
The terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of example embodiments of the invention. As used herein, the singular forms "a", "an" and "the" are intended to include the plural forms as well, unless the context clearly indicates to the contrary. It will be further understood that the terms "comprises," "comprising," "includes," "including" and/or "including," when used herein, specify the presence of stated features, integers, steps, operations, elements, and/or components, and do not preclude the presence or addition of one or more other features, numbers, steps, operations, elements, components, and/or groups thereof.
It should also be noted that in some alternative implementations, the functions/acts noted may occur out of the order noted in the figures. For example, two figures shown in succession may in fact be executed substantially concurrently or the figures may sometimes be executed in the reverse order, depending upon the functionality/acts involved.
In the following description, specific details are provided to provide a thorough understanding of example embodiments. However, it will be understood by those of ordinary skill in the art that the example embodiments may be practiced without these specific details. For example, a system may be shown in block diagrams in order to avoid obscuring the examples with unnecessary detail. In other embodiments, well-known processes, structures, and techniques may not be shown in unnecessary detail in order to avoid obscuring the example embodiments.
Example 1:
the embodiment provides a dual-medium carbonate pore-crack dual-network model construction method, which is as shown in fig. 1:
the method comprises the following steps:
step one: and carrying out CT scanning on the dual-medium carbonate rock sample by using the micron CT to obtain a nondestructive real sample three-dimensional image.
Step two: and (3) extracting pore spaces and crack spaces from the three-dimensional image obtained in the step one by utilizing an image segmentation technology. Because the pores and the cracks are extracted as the same phase, the pores and the cracks need to be further separated, and the aim of respectively constructing the network model is fulfilled.
Step three: and (3) carrying out morphological analysis on the pore and crack space extracted in the step two, calculating the three-dimensional length and the three-dimensional minimum width information of each pore/crack, and dividing the pore and the crack through the ratio of the two.
Step four: firstly, a central line extraction method is utilized to truly draw out three-dimensional space spread and topology structures of the crack in the crack space separated in the step three, then, a continuous expansion searching method is adopted by taking each central line voxel as a sphere center, the crack boundary is ascertained, the mutually connected spheres are used for representing the real opening of the crack, stick-shaped construction is used for connecting the spheres, and connectivity and seepage channel information of the crack are represented; and simultaneously, constructing a network model for the pore space separated in the step three by using the sphere filling method to obtain a pore network model.
Step five: and (3) after the pore network model and the crack network model are respectively obtained according to the step four and the step five, constructing the pore-crack dual network model by adding a communication throat between the maximum spheres representing the pores and the cracks.
As shown in fig. 2, the original dual medium carbonate micron CT scan image is shown in fig. 3, which is the pore network model, fracture network model and final pore-fracture dual network model finally obtained by the present method.
The method for establishing the digital rock core by using the gray level image obtained by the micron CT scanning mainly comprises the following three steps: 1. scanning the processed rock sample by using CT equipment to obtain projection data of the rock sample; 2. reconstructing the obtained projection data into a core gray image by using an image reconstruction method, as shown in a part a in fig. 4; 3. after preprocessing the gray level image, the void and crack space in the gray level image is separated from the rock skeleton by using an image segmentation method, as shown in part b in fig. 4.
Example 2:
as an optimization of the above embodiment, the void space and the crack space can be distinguished by calculating the ratio of the three-dimensional length and the three-dimensional minimum width of the void and the crack, and selecting a suitable threshold value for effective separation. The three-dimensional length measurement is to place each pore/crack in a coordinate system as shown in fig. 5, take 18 degrees as a step length along the theta and phi directions, search from 0 degrees to 162 degrees, record the distance between the head and tail ends of each step, wherein the maximum distance is the three-dimensional length. Two planes which are parallel to each other and are in contact with the edges of the pores/cracks and are made in the three-dimensional length direction rotate by taking the three-dimensional length direction as an axis, and the minimum distance between the two planes is the three-dimensional minimum width. Because the crack is in a flat form, the ratio of the three-dimensional length to the minimum width is far greater than that of the pore, so that the crack and the pore can be distinguished by the size of the ratio. The result of distinguishing the pores and cracks of fig. 4 using the ratio of three-dimensional length to minimum width is shown in fig. 6.
Example 3:
as an optimization of the above embodiment, when the network model is constructed for the fracture space separated in the third step, the following operations are sequentially adopted to obtain a fracture network model:
1. extracting a central line: firstly, giving a value to each voxel of a fracture digital core model, wherein the value is the distance from the voxel to the nearest rock skeleton; and then continuously refining the digital core model of the fracture according to the distance map until a continuous voxel chain is left, wherein the voxel chain is the extracted fracture center line as shown in fig. 7.
2. With each centerline voxel on the centerline as a base point, a dilation algorithm is first used to determine the possible search space. After finding the nearest skeleton voxel range from each voxel, adopting a contraction algorithm to identify the inscribed sphere corresponding to the center voxel; and finally, calculating to obtain the upper limit value and the lower limit value of the radius of the inscribed sphere.
3. After a series of inscribed spheres centered on the centerline voxels are extracted, the extracted inscribed spheres need to be optimized because some spheres are a subset of other spheres. The optimization method is characterized in that each voxel is represented in a linked list mode according to the square size of the inscribed sphere radius corresponding to the voxel. When the loop calculation is to r2=0, the operation is stopped. For an inscribed sphere with a sphere center of p0= (x 0, y0, z 0) and a radius of R0, searching an inscribed sphere corresponding to a voxel p meeting the condition dist2 (p, p 0) < R02 in the generated linked list, and judging whether an inclusion relationship exists between the inscribed sphere and the appointed sphere. The found inscribed sphere with the containing relation is deleted from the linked list, and finally the largest inscribed sphere set which is not contained by other spheres is left in the linked list.
4. Further, the largest inscribed ball set is converted into clusters connected with each other, as shown in fig. 8, the main balls in each cluster define the size of the crack space in the area where the cluster is located, the nodes connected between different clusters define the size of the seepage channels of the crack space, after one seepage channel is found, the links of the channel to two main balls can be obtained, and the links are composed of the seepage channels and all inscribed balls between the main balls. As shown in fig. 9, these links also actually form the skeleton of the overall fracture network model.
Example 4:
as an optimization for the above embodiment, after the fracture network model and the pore network model are respectively constructed, the following operation flows are sequentially adopted, and a communication throat is added between the largest inscribed spheres representing the pores and the fractures, so that the construction of the pore-fracture dual network model is realized:
1. after the pore space and the crack space are obtained, as shown in fig. 10, the contact surface of the pore space and the crack space, which are in contact with each other, is marked, and the largest inscribed sphere number in the pore network model and the crack network model corresponding to the contact surface is found.
2. According to the number N and the area S of each contact surface and the corresponding number I of the largest inscribed sphere, establishing a throat radius and coordination number distribution relation of communication between the pores and the cracks.
3. As shown in fig. 11 to 12, according to the distribution relation of the throat radius and coordination number, the pore and the maximum sphere of the crack corresponding to the contact surface are selected to connect the two throats with proper radius and number, so as to construct the pore-crack dual network model.
The method accurately describes the three-dimensional space topological structure of the crack by a method of extracting the central line and adding the maximum inscription sphere, finely characterizes the information of the space spread, the opening degree, the connectivity and the like of the crack, and solves the problem of accurately constructing the crack network model when constructing the dual-medium carbonate network model; in addition, the method builds a communication throat between the pore and the crack based on the real communication surface and the communication relation of the pore and the crack, and solves the problem that two mediums of the pore and the crack are respectively built into a network model and then are fused into a set of pore-crack dual network model.
The invention is not limited to the alternative embodiments described above, but any person may derive other various forms of products in the light of the present invention. The above detailed description should not be construed as limiting the scope of the invention, which is defined in the claims and the description may be used to interpret the claims.
Claims (4)
1. The method for constructing the dual-medium carbonate pore-crack dual-network model is characterized by comprising the following steps of:
step one: collecting a dual-medium carbonate sample, and scanning the dual-medium carbonate sample to obtain a nondestructive three-dimensional image;
step two: performing image segmentation processing on the three-dimensional image, and extracting a space to be identified, which contains a pore space and a crack space, from the three-dimensional image;
step three: performing morphological analysis on the space to be identified, and completing separation of the pore space and the crack space according to morphological analysis results, wherein the method comprises the following steps: s31, measuring and obtaining three-dimensional length and three-dimensional minimum width data of a pore space and a crack space; s32, calculating to obtain the ratio of the three-dimensional length to the three-dimensional minimum width; s33, setting a threshold value, comparing the calculated ratio with the threshold value, and judging a crack space or a pore space according to a comparison result;
step four: respectively constructing a network model for the separated pore space and the separated crack space by adopting a sphere filling method to obtain a pore network model and a crack network model, wherein the sphere filling method comprises the steps of adopting a central line extraction method to draw three-dimensional space spreading and topological structures, then taking a central line voxel as a sphere center, taking a space boundary point as a tangent point to establish an inscribed sphere, finally establishing a rod body to connect each continuous inscribed sphere in a centering manner, representing the opening degree of the space by the interconnected inscribed spheres, and representing the connectivity of the space and the seepage channel information by the connected rod body; the specific steps of constructing the network model for the crack space are as follows: s41, center line extraction: extracting continuous voxel chains, wherein the voxel chains are crack center lines of a crack space; s42, determining an expandable space by using each voxel on a central line as a basic point and adopting an expansion algorithm, after finding the range of the double-medium carbonate skeleton voxels nearest to each voxel, identifying an inscribed sphere corresponding to the basic point by adopting a contraction algorithm, and finally calculating to obtain an upper limit and a lower limit of the inscribed sphere radius; s43, after a series of inscribed spheres taking a central line voxel as a sphere center are extracted, optimizing the extracted inscribed spheres to obtain a residual inscribed sphere set; s43, sequentially centering and connecting the rest inscribed sphere sets through the rod bodies to form a crack space for a network model;
step five: a communication throat is added between the pore network model and the crack network model, so that the construction of a pore-crack dual network model is completed; the specific steps for constructing the pore-crack dual network model comprise: s51, marking contact surfaces of a pore space and a crack space, finding out respective maximum inscribed spheres in a pore network model and a crack network model corresponding to the contact surfaces, and numbering the maximum inscribed spheres; s52, obtaining a throat radius and coordination number distribution relation communicated between the pore network model and the crack network model according to the number and the area of the contact surfaces and the corresponding number of the maximum inscribed sphere; and S53, connecting the largest inscribed spheres in the pore network model and the crack network model corresponding to the contact surface through the throat according to the throat radius and coordination number distribution relation to form a pore-crack dual network model.
2. The dual medium carbonate pore-fracture dual network model construction method according to claim 1, characterized in that: in the first step, a CT scanner is adopted to carry out CT scanning on a dual-medium carbonate rock sample, and a nondestructive three-dimensional image is obtained.
3. The dual medium carbonate pore-fracture dual network model construction method according to claim 2, characterized in that: the specific process for acquiring the three-dimensional image of the dual-medium carbonate sample comprises the following steps: firstly, scanning a dual-medium carbonate rock sample by using a CT scanner to obtain projection data; then, reconstructing the projection data into a core gray image by using an image reconstruction method; and finally, preprocessing the gray level image to obtain a three-dimensional image of the dual-medium carbonate rock sample.
4. The dual medium carbonate pore-fracture dual network model construction method according to claim 1, characterized in that: in step S43, the process of optimizing the extracted inscribed sphere includes: firstly, according to the square size of the radius of each voxel corresponding to the inscribed sphere, each voxel is represented in a linked list form; and searching the linked list for the inscribed sphere corresponding to the voxel with the inclusion relationship, deleting the included inscribed sphere with the inclusion relationship, and finally leaving the rest inscribed sphere set which is not included.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911150076.1A CN110853138B (en) | 2019-11-21 | 2019-11-21 | Construction method of dual-medium carbonate pore-crack dual-network model |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911150076.1A CN110853138B (en) | 2019-11-21 | 2019-11-21 | Construction method of dual-medium carbonate pore-crack dual-network model |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110853138A CN110853138A (en) | 2020-02-28 |
CN110853138B true CN110853138B (en) | 2023-08-18 |
Family
ID=69603447
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911150076.1A Active CN110853138B (en) | 2019-11-21 | 2019-11-21 | Construction method of dual-medium carbonate pore-crack dual-network model |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110853138B (en) |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111754500B (en) * | 2020-06-30 | 2023-06-27 | 中国科学院地质与地球物理研究所 | Rock fracture network topology structure describing system |
CN112082917B (en) * | 2020-08-03 | 2021-04-30 | 西南石油大学 | Gas-water unsteady two-phase seepage simulation method based on dynamic network simulation |
CN112098293B (en) * | 2020-08-03 | 2021-06-01 | 西南石油大学 | Unsteady gas-water two-phase seepage simulation method based on pore fracture dual-medium gas reservoir |
CN112233166A (en) * | 2020-09-11 | 2021-01-15 | 安徽理工大学 | Pore size distribution evaluation method based on porous medium three-dimensional pore space image |
CN111932642B (en) * | 2020-09-27 | 2021-01-01 | 湖南大学 | Method, device and equipment for measuring and calculating volume of structural crack and storage medium |
CN112132966A (en) * | 2020-09-29 | 2020-12-25 | 成都理工大学 | Shale fracture network connectivity characterization method based on topological structure |
CN112504928B (en) * | 2020-10-14 | 2022-11-04 | 中国石油天然气股份有限公司 | Method and device for determining connectivity of fractures in reservoir rock |
CN113160107A (en) * | 2020-11-03 | 2021-07-23 | 清能艾科(深圳)能源技术有限公司 | Rock core micro-crack extraction method and device, electronic equipment and storage medium |
CN113310877B (en) * | 2021-06-07 | 2022-12-06 | 中国石油大学(华东) | Method for constructing multi-scale rock pore network model with differentiated regional structure |
CN113405966B (en) * | 2021-06-08 | 2022-08-23 | 浙江广天构件集团股份有限公司 | Method for calculating pore size distribution of cement-based material particle accumulation system |
CN114220094B (en) * | 2022-02-22 | 2022-04-26 | 中国科学院地质与地球物理研究所 | Pore and crack identification method and system based on two-dimensional rock core scanning image |
US11952891B2 (en) | 2022-08-22 | 2024-04-09 | Saudi Arabian Oil Company | Systems and method for constraining 3D fracture model properties using X-ray micro-computed tomography of core plugs for naturally fractured reservoirs |
CN116152259B (en) * | 2023-04-23 | 2023-07-04 | 西南石油大学 | Reservoir permeability calculation method based on graph neural network |
CN117745979B (en) * | 2024-02-21 | 2024-05-03 | 山东科技大学 | Three-dimensional fracture-pore coupling network simulation generation method and system |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2001073476A1 (en) * | 2000-03-27 | 2001-10-04 | Ortoleva Peter J | Method for simulation of enhanced fracture detection in sedimentary basins |
CN105260543A (en) * | 2015-10-19 | 2016-01-20 | 中国石油天然气股份有限公司 | Double-pore model-based multi-medium oil gas flow simulation method and device |
CN109033672A (en) * | 2018-08-09 | 2018-12-18 | 中国石油天然气股份有限公司 | The dynamic crack determination method and device of pore constriction network model in displacement simulation |
CN110126058A (en) * | 2019-05-20 | 2019-08-16 | 重庆大学 | A kind of rock sample preparation method based on CT visualization and 3D printing |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101929973B (en) * | 2009-06-22 | 2012-10-17 | 中国石油天然气股份有限公司 | Quantitative calculation method for hydrocarbon saturation of fractured reservoir |
-
2019
- 2019-11-21 CN CN201911150076.1A patent/CN110853138B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2001073476A1 (en) * | 2000-03-27 | 2001-10-04 | Ortoleva Peter J | Method for simulation of enhanced fracture detection in sedimentary basins |
CN105260543A (en) * | 2015-10-19 | 2016-01-20 | 中国石油天然气股份有限公司 | Double-pore model-based multi-medium oil gas flow simulation method and device |
CN109033672A (en) * | 2018-08-09 | 2018-12-18 | 中国石油天然气股份有限公司 | The dynamic crack determination method and device of pore constriction network model in displacement simulation |
CN110126058A (en) * | 2019-05-20 | 2019-08-16 | 重庆大学 | A kind of rock sample preparation method based on CT visualization and 3D printing |
Non-Patent Citations (1)
Title |
---|
基于孔与裂隙网络模型的平行微裂隙对驱油的影响规律研究;刘海娇 等;力学学报;第50卷(第04期);890-898 * |
Also Published As
Publication number | Publication date |
---|---|
CN110853138A (en) | 2020-02-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110853138B (en) | Construction method of dual-medium carbonate pore-crack dual-network model | |
CN108763711B (en) | Permeability prediction method based on rock core scanning image block numerical simulation | |
CN101556703B (en) | Method for establishing network model based on serial section image | |
CN107817199A (en) | A kind of construction method of tight sand multi-scale porosity model and application | |
CN113609696B (en) | Multi-scale multi-component digital core construction method and system based on image fusion | |
CN109242985B (en) | Method for determining key parameters of pore structure from three-dimensional image | |
CN105551004A (en) | Core CT image processing-based remaining oil micro-occurrence representing method | |
CN113963130B (en) | Construction method of fracture network model for rock core fracture | |
CN106324002A (en) | Carbonatite pore structure characterization method based on rock classification and multi-scale digital cores | |
CN104885124A (en) | Method for producing a three-dimensional characteristic model of a porous material sample for analysis of permeability characteristics | |
CN109900617B (en) | Method for calculating permeability curve of fractured formation based on acoustoelectric imaging log | |
CN110320137A (en) | A kind of Multiscale Fusion method based on digital cores | |
CN116152259B (en) | Reservoir permeability calculation method based on graph neural network | |
CN111739149B (en) | Oil-water distribution continuity restoration method for rock CT scanning image | |
CN113310877B (en) | Method for constructing multi-scale rock pore network model with differentiated regional structure | |
CN112903555A (en) | Porous medium permeability calculation method and device considering pore anisotropy | |
CN115309846B (en) | Road network structure identification method based on parallel coefficients | |
US11867869B2 (en) | Multiple porosity micromodel | |
CN115221792A (en) | Natural gas hydrate type dividing device and method based on machine learning | |
CN114705606A (en) | Blocking method of key seepage nodes in rock based on networked analysis | |
CN114663640A (en) | Submarine geographic entity demarcation and classification method based on landform and structural characteristics | |
CN110441204B (en) | Digital core simulation-based compact reservoir fracturing fluid damage digital evaluation method | |
CN113870200B (en) | Method for quantitatively analyzing pore coordination number in sandstone and carbonate reservoir | |
CN114299098A (en) | Rapid rock permeability analysis method based on throat characteristics of scanned image | |
CN117853766B (en) | Tunnel fracture coplanarity matching method and system based on tunnel face and borehole image |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |