Method and device for constructing a medium-filled three-dimensional digital core model

Abstract

The present application relates to oil and gas exploration and development technical field, especially to a kind of medium filling three-dimensional digital core model construction method and device.The method includes, based on two-dimensional core picture, the type, distribution form and filling medium content of core filling medium are obtained;Based on the distribution form and filling medium content of filling medium, the three-dimensional digital core model of medium in pore center filling and the three-dimensional digital core model of medium in pore edge filling are constructed using the three-dimensional digital core model of identifying skeleton and pore;Based on the three-dimensional digital core model of medium in pore center filling and the three-dimensional digital core model of medium in pore edge filling, medium filling three-dimensional digital core model is constructed.The present application studies the influence of different lithology and different filling form on core petrophysical property, so as to improve the ability of identification and quantitative evaluation of filling medium type and content.

Description

Technical Field

[0001] This invention relates to the field of oil and gas exploration and development technology, and in particular to a method and apparatus for constructing a three-dimensional digital core model filled with a medium. Background Technology

[0002] Oil and gas exist in the pore spaces of reservoirs. During reservoir formation and hydrocarbon migration, diagenesis and differences in the thermal history of hydrocarbon generation result in multiple stages of mineral filling in the reservoir pores and different occurrence forms of these filling minerals. Exploration practice shows that in the Sichuan Basin's marine carbonate gas fields and in the Ordos, Junggar, and Songliao Basins, the pore spaces of tight sandstone reservoirs are commonly filled with minerals such as bitumen, quartz, and dolomite. These minerals, acting as filling media, occupy the reservoir space, disrupt the pore structure, reduce reservoir effectiveness, and severely affect reservoir properties and productivity. Therefore, clarifying the impact of filling media on reservoir logging response characteristics and reservoir effectiveness is crucial for stratigraphic evaluation of carbonate and tight sandstone reservoirs. Furthermore, studying the influence of changes in filling media content and morphology on pore structure and properties helps in understanding the evolution of unconventional oil and gas reservoirs in ultra-deep and ultra-old strata, providing assistance for unconventional oil and gas sedimentology research.

[0003] The filling medium occupies part of the reservoir space and disrupts the pore structure of the reservoir, which not only reduces the porosity of the reservoir but also significantly reduces its permeability, seriously affecting the reservoir's physical properties and productivity. Therefore, the identification of the filling medium type and the determination of the filling medium content are crucial for reservoir logging evaluation. However, current research on filling media mainly focuses on the formation mechanism and genetic identification, geochemical characteristics, distribution characteristics, and their impact on reservoir performance. There is relatively little research on the rock physical characteristics of asphalt reservoirs and asphalt logging evaluation methods in asphalt-filled reservoirs. Quartz, calcite, and dolomite filling sandstone result in poor physical properties, low porosity and permeability of tight sandstone, and a complex pore structure, leading to complex rock resistivity characteristics and obvious non-Archite phenomena. As a result, the saturation evaluation of tight sandstone reservoirs remains a major challenge in reservoir logging evaluation both domestically and internationally.

[0004] To address this, researchers employed experimental methods, using core samples before and after media dissolution to study physical properties, density, rock acoustic transit time, resistivity, and nuclear magnetic resonance (NMR) experiments. This research investigated the impact of filling media on the petrophysical characteristics and logging response of carbonate and tight sandstone reservoirs. However, dissolving filling media in the laboratory presents challenges. Furthermore, for tight sandstone reservoirs, low permeability and difficulty in displacement result in non-wetting phase saturation generally below 60% during rock electrical experiments, failing to reflect the electrical properties across the entire water saturation range and thus differing from actual reservoir conditions. The lack of these logging parameters and experimental difficulties in practical applications limit the application of such methods. Numerical simulation based on three-dimensional digital cores integrates theoretical and experimental rock physics, effectively complementing rock physics. Therefore, core numerical simulation technology allows for flexible simulation of changes in filling media type and content, enabling the study of the influence of filling media content, type, and morphology on core pore structure and rock resistivity. Compared to traditional laboratory rock physics experiments, core-based three-dimensional numerical simulations are more efficient, reusable, and facilitate the study of the influence of single microscopic factors.

[0005] Reservoir logging response is directly controlled by the mineral and fluid composition of the reservoir. Previous research generally agreed that the presence of porosity in the filling medium leads to a "three-rise, two-fall" characteristic in the logging response: increased acoustic transit time (AC, DTC), increased resistivity and natural gamma, and decreased neutrons and density. However, under the same lithology, cores with significantly different porous fillings show decreased acoustic transit time, contradicting previous findings. Furthermore, under the same porous filling, the porosity, permeability, and rock physical properties of different cores differ from previous studies. Therefore, it is difficult to determine whether the logging response changes caused by lithofacies or those caused by the filling are stronger, increasing the difficulty of quantitative logging evaluation. Numerical simulation techniques can easily construct three-dimensional digital core models of different lithological fillings, thus laying the foundation for studying the impact of different fillings and the same filling of different lithologies on resistivity. Summary of the Invention

[0006] The purpose of this invention is to provide a method and apparatus for constructing a three-dimensional digital core model with medium filling. This innovative achievement provides a method for constructing a three-dimensional digital core with medium filling and simulating and analyzing its rock physical properties. The aim is to construct three-dimensional digital cores of carbonate rocks with different filling medium contents and forms, and to study the influence of different lithologies and filling forms on the rock physical properties of the cores, thereby improving the ability to identify and quantitatively evaluate the type and content of filling media. This invention contributes to proposing a method for constructing a three-dimensional digital core model with medium filling, providing innovative ideas and theoretical guidance for well logging evaluation of marine carbonate rocks in the Sichuan Basin and tight sandstone reservoirs in the Ordos, Junggar, and Songliao Basins. Simultaneously, the correct identification of the filling medium and the accurate evaluation of the reservoir filling medium content can not only effectively guide the optimization of well logging data acquisition and gas testing, but also provide support for the deployment of gas reservoir development wells. Furthermore, it has significant theoretical and practical implications for finding ancient oil reservoirs, exploring their distribution, scale, and formation mechanisms, and demonstrating the oil and gas exploration potential of the study area. To achieve the above objectives, this invention provides the following technical solution:

[0007] This invention provides a method for constructing a three-dimensional digital core model filled with a medium, the method comprising:

[0008] Based on two-dimensional core images, the type, distribution morphology, and content of the core filling medium were obtained;

[0009] Based on the distribution morphology and content of the filling medium, a three-dimensional digital core model with the medium filling the center of the pores and a three-dimensional digital core model with the medium filling the edge of the pores are constructed using a three-dimensional digital core model that identifies the skeleton and pores.

[0010] Based on the three-dimensional digital core model with the medium filling the pore center and the three-dimensional digital core model with the medium filling the pore edge, a three-dimensional digital core model with medium filling is constructed.

[0011] Furthermore, the three-dimensional digital core model filled with the construction medium at the pore centers includes,

[0012] A three-dimensional digital core model that identifies the skeleton and pores is used as the initial three-dimensional digital core model for filling the pore center of the medium.

[0013] Based on the initial three-dimensional digital core model of the filling medium distribution morphology, filling medium content and filling of the pore center of the medium, a three-dimensional digital core model of the medium filling the pore center is constructed using image morphology method.

[0014] Furthermore, the initial three-dimensional digital core model based on the distribution morphology, content, and filling of the pore centers of the filling medium is used to construct a three-dimensional digital core model with the medium filling the pore centers using image morphology methods, including:

[0015] Based on the initial three-dimensional digital core model filled with the pore center of the medium, the first skeleton set and the first pore set are obtained, and the number of elements in the first pore set is obtained.

[0016] Based on image morphology, an opening operation is performed on spherical structural elements of different diameters in the first pore set to obtain the first opening operation result set, and the number of first pore elements in the first opening operation result set and the proportion of the number of first pore elements in the first pore set are obtained.

[0017] Based on the first set of opening operation results, obtain the difference set of the opening operation results sets of two adjacent diameter sphere structuring elements.

[0018] Perform Euclidean distance calculation on the difference set to obtain the number of pore elements in the corresponding Euclidean distance set and the proportion of the corresponding number of pore elements in the first pore set of the initial three-dimensional digital core model filled with the pore center of the medium.

[0019] The set of opening operation results for spherical structural elements of different diameters, the number of pore elements in the set of opening operation results for each diameter, the proportion of the number of pore elements in the set of opening operation results for each diameter in the pore set of the three-dimensional digital core model filled with medium pores, the difference set of the set of opening operation results for spherical structural elements of two adjacent diameters, the set of Euclidean distances in the difference set of the set of opening operation results for spherical structural elements of adjacent diameters, the proportion of the number of pore elements in each Euclidean distance in the difference set in the pore set of the initial three-dimensional digital core model filled with medium pores, and the content of the filling medium are used to construct a three-dimensional digital core model filled with medium pores.

[0020] Furthermore, the initial three-dimensional digital core model filled with the pore centers of the medium is represented by the first initial three-dimensional digital core set CORE1, wherein,

[0021] CORE1 = MATRIX 中心 ∪PORE 中心 ;

[0022] In the formula, MATRIX 中心 PORE represents the first skeleton set. 中心 Let U denote the first set of pores, and U denote the union of the first set of skeletons, MATRIX. 中心 Represented as:

[0023] MATRIX 中心 ={px 中心 |phase(px 中心 )=0};

[0024] First Pore Set PORE 中心 Represented as:

[0025] PORE 中心 ={px 中心 |phase(px 中心 )=1};

[0026] In the formula, px 中心 This represents an element in the initial three-dimensional digital core model filled with pore centers of the medium, where phase represents the phase state function, phase(px 中心 ) = 0 means px 中心 The element is a skeleton, phase(px) 中心 ) = 1 represents px 中心 The element is porosity.

[0027] Furthermore, the three-dimensional digital core model of the medium filling the pore center includes a pore space model of the medium gradually filling the pore center. The pore space model of the medium gradually filling the pore center is obtained through the three-dimensional digital core assembly CORE of the center-filled core. 中心 Specifically, it means:

[0028] CORE 中心 =MATRIX 中心 ∪PORE 中心 ∪BIT 中心 ;

[0029] BIT 中心 ={px 中心 |phase(px 中心 )=3};

[0030] Among them, CORE 中心 MATRIX represents a three-dimensional digital core assembly filled with media at its center. 中心 Represents the first skeleton set, PORE 中心 Represents the first set of pores, BIT 中心 Represents the centrally filled medium set, phase(px) 中心 ) = 3 represents px 中心 The element is the filling medium, and phase represents the phase function.

[0031] Furthermore,

[0032] If OSB dmax >TBtarget 中心 Then MORPH dmax All the pores in the middle are transformed into the filling medium, and the central filling medium set is BIT. 中心 Represented as:

[0033] BIT中心 =MORPH dmax

[0034] MORPH dmax ={px 中心 |OPEN(px 中心 ,dmax)=1};

[0035] The number of elements in the filling medium set is Onum dmax ;

[0036] In the formula, TBtarget 中心 Indicates the content of the first filling medium; MORPH dmax Represents the set of results of the opening operation on the sphere structuring element with the largest diameter dmax in the first set of pores; Onum dmax OSB represents the number of elements in the set of sphere structure elements satisfying the maximum diameter dmax based on image morphology within the first set of pores; dmax Represents Onum dmax First pore set PORE 中心 The proportion of the number of elements;

[0037] OPEN represents the opening operator, px 中心 This represents the element in the initial three-dimensional digital core model filled with pores in the medium, where dmax represents the maximum diameter of the spherical structural element, and OPEN(px) represents the maximum diameter of the spherical structural element. 中心 dmax) = 1 indicates that px 中心 The result of the opening operation of the sphere structuring element with the maximum diameter dmax.

[0038] Furthermore, if OSB i <TBtarget 中心 <OSB i-1 And DSB1>TBtarget 中心 -OSB i When (i = 2, 3, 4, ..., dmax), then MORPH i All the pores in the medium are transformed into a filling medium.

[0039] The set of filling centers at the pore centers of the medium is represented as follows:

[0040] BIT 中心 =BIT1 中心 ∪BIT2 中心 ;

[0041] BIT1 中心 =MORPH i MORPH i ={px 中心 |OPEN(px中心 ,i)=1};

[0042] BIT2 中心 ={px 中心 Randomly select px 中心 , px 中心 ∈DIST1};

[0043] DMORPH i-1 =MORPH i-1 -MORPH i ;

[0044] DIST1 = {px 中心 |dist(px 中心 ) = 1, and px 中心 ∈DMORPH i-1};

[0045] In the formula, BIT 中心 BIT1 represents the centrally filled medium set. 中心 Indicates center-filled medium set 1, BIT2 中心 Indicates center-filled medium set 2, MORPH i OSB represents the set of results of the opening operation on the sphere structuring element with diameter i in the first set of pores. i MORPH i The elements in the first pore set PORE 中心 The proportion of elements in the total number; MORPH i-1 OSB represents the set of results from the opening operation of sphere structuring elements with diameter i-1 in the first set of pores. i-1 MORPH i-1 The elements in the first pore set PORE 中心 Percentage of elements; px 中心 This represents an element in the initial three-dimensional digital core model filled with pores in the medium; OPEN represents the opening operator; i, i-1 represent the diameters of the spherical structural elements, respectively; OPEN(px 中心 , i) = 1 means px 中心 The result of the opening operation of a sphere structuring element with diameter i is satisfied; DMORPH i-1 DIST1 represents the difference set of the results of the opening operation between sphere structural elements with diameters i-1 and i in the first pore set based on image morphology; DIST1 represents DMORPH. i-1 The set consisting of elements in the set whose Euclidean distance is 1; DSB1 represents the set of elements in DIST1 in the first pore set PORE. 中心 The percentage of elements; dist represents the Euclidean distance operator.

[0046] Furthermore,

[0047] If OSB i <TBtarget<OSB i-1 and When i = 2, 3, 4, ..., dmax),

[0048] The set of filling centers at the pore centers of the medium is represented as follows:

[0049] BIT 中心 =BIT1' 中心 ∪BIT2' 中心 ∪BIT3' 中心 ;

[0050] BIT1' 中心 =MORPH i MORPH i ={px 中心 |OPEN(px 中心 ,i)=1};

[0051]

[0052] BIT3' 中心 ={px 中心 Randomly select px 中心 , px 中心 ∈DIST j+1};

[0053] DMORPH i-1 =MORPH i-1 -MORPH i ;

[0054] DIST j ={px 中心 |dist(px 中心 ) = j, and px 中心 ∈DMORPH i-1};

[0055] DIST j+1 ={px 中心 |dist(px 中心 ) = j + 1, and px 中心 ∈DMORPH i-1};

[0056] In the formula, BIT 中心 BIT1′ represents the set of center-filled media. 中心 This indicates the center-filled medium set 1', BIT2' 中心 This indicates the central filling medium set 2', BIT3' 中心Indicates center-filled medium set 3'; MORPH i MORPH represents the set of results from the opening operation of spherical structuring elements with diameter i in the first set of pores. i-1 DMORPH represents the set of results from the opening operation of spherical structuring elements with diameter i-1 in the first set of pores. i-1 OSB represents the difference set of the results of the opening operation between sphere structuring elements of diameters i-1 and i in the first set of pores, based on image morphology. i MORPH i The elements in the first pore set PORE 中心 Percentage of elements; OSB i-1 MORPH i-1 The elements in the first pore set PORE 中心 Percentage of elements; px 中心 This represents an element in the initial three-dimensional digital core model filled with pores in the medium; OPEN represents the opening operator; i, i-1 represent the diameters of the spherical structural elements, respectively; OPEN(px 中心 , i) = 1 indicates that the center of px satisfies the opening operation result of a sphere with diameter i; DIST j Indicates DMORPH i-1 The set of elements whose Euclidean distance is j; DSB j DIST j The elements in the first pore set PORE 中心 The proportion of elements in the total number; DIST j+1 Indicates DMORPH i-1 The set of elements whose Euclidean distance is j+1; DSB j+1 DIST j+1 The elements in the first pore set PORE 中心 The percentage of elements; dist represents the Euclidean distance operator; dist(px 中心 ) indicates that the element px is calculated. 中心 Euclidean distance; ∪ denotes the union, ∪ j＝1 i DIST j Denotes sets DIST1, DIST2, ..., DIST i The union of .

[0057] Furthermore, a three-dimensional digital core model of the medium filling the pore edges is constructed, including:

[0058] Based on the three-dimensional digital core model that identifies the skeleton and pores, an initial three-dimensional digital core model with pore edge filling is constructed.

[0059] Based on the initial three-dimensional digital core model of the filling medium distribution morphology, filling medium content and filling of the pore edges, a three-dimensional digital core model of the medium filling the pore edges is constructed using the image Euclidean distance map method.

[0060] Furthermore, the method of constructing a three-dimensional digital core model of the medium filling the pore edges using an image Euclidean distance map includes:

[0061] Based on the initial three-dimensional digital core model filled with the pore edges of the medium, the second skeleton set and the second pore set are obtained, and the number of elements in the second pore set is obtained.

[0062] The Euclidean distance method is used to calculate the Euclidean distance from the elements in the second pore set to the second skeleton set.

[0063] Based on the Euclidean distance from the elements in the second pore set to the second skeleton set, the second pore set is divided into multiple second Euclidean distance sets, and the number of elements in each second Euclidean distance set and the proportion of the number of elements in each Euclidean distance set in the second pore set are obtained.

[0064] Based on multiple sets of second Euclidean distances, the proportion of the number of elements in each Euclidean distance set in the second pore set, and the content of the filling medium, a three-dimensional digital core model of the filling at the edge of the medium is constructed.

[0065] Furthermore, the initial three-dimensional digital core model filled with the pore edges of the medium is represented by the second initial three-dimensional digital core set CORE2, wherein,

[0066] CORE2 = MATRIX 边缘 ∪PORE 边缘 ;

[0067] In the formula, MATRIX 边缘 PORE represents the second skeleton set. 边缘 Let U denote the second set of pores, and U denote the union of the second skeleton set, MATRIX. 边缘 Represented as:

[0068] MATRIX 边缘 ={px 边缘 |phase(px 边缘 )=0};

[0069] Second pore set PORE 边缘 Represented as:

[0070] PORE 边缘 ={px 边缘 |phase(px 边缘 )=1};

[0071] In the formula, px 边缘 This represents an element in a three-dimensional digital core set filled with pore edges of the medium, where phase represents the phase state function, phase(px 边缘 ) = 0 means px 边缘 The element is a skeleton, phase(px) 边缘 ) = 1 represents px 边缘 The element is porosity.

[0072] Furthermore, the three-dimensional digital core model with pore edge filling includes a pore space model where the medium gradually fills the pore edges. This pore space model, with the medium gradually filling the pore edges, is connected to the three-dimensional digital core model CORE with center filling. 边缘 Specifically, it means:

[0073] CORE 边缘 =MATRIX 边缘 ∪PORE 边缘 ∪BIT 边缘 ;

[0074] BIT 边缘 ={px 边缘 |phase(px 边缘 )=3};

[0075] Among them, CORE 边缘 This represents a 3D digital core model after edge filling; MATRIX 边缘 Represents the second skeleton set; PORE 边缘 Represents the second set of pores; BIT 边缘 Represents the set of edge-filling media, phase(px) 边缘 ) = 3 indicates that the element is used to fill the medium; px 边缘 This represents an element in the initial three-dimensional digital core model filled with pore edges of the medium; phase represents the phase function.

[0076] Furthermore, if PDSB1 > TBtarget 边缘 Then, nump is randomly selected from the PDIST1 set. 边缘 ×TBtarget 边缘 Each pore element is transformed into a filling medium, and the edge filling medium set BIT in the three-dimensional digital core model of the medium edge filling is established. 边缘 ,

[0077] BIT 边缘 ={px 边缘 Randomly select px 边缘 , px 边缘 ∈PDIST1};

[0078] PDIST1 = {px 边缘 |dist(px 边缘 ) = 1, and px 边缘 ∈PORE2};

[0079] The number of elements in the filling medium set is nump 边缘 ×TBtarget 边缘 ;

[0080] In the formula, BIT 边缘 TBtarget represents the edge-filled media set. 边缘 Indicates the content of the edge-filling medium, nump 边缘 PORE represents the second set of pores filled at the edge. 边缘 The number of elements, px 边缘 This represents an element in the initial three-dimensional digital core model filled with pore edges of the medium. `dist` represents the Euclidean distance operator, and `PDIST1` represents the element in the second pore set `PORE`. 边缘 The second Euclidean distance set consists of elements with a distance of 1 in the middle Euclidean distance set. PDSB1 represents the proportion of the number of elements with a distance of 1 in the second Euclidean distance set in the second pore set.

[0081] Furthermore, when hour,

[0082] The filling medium set is as follows:

[0083] BIT 边缘 =BIT1 边缘 ∪BIT2 边缘 ;

[0084]

[0085] BIT2 边缘 ={px 边缘 Randomly select px 边缘 , px 边缘 ∈PDIST i+1};

[0086] PDIST i+1 ={px 边缘 |dist(px 边缘 ) = i + 1, and px 边缘 ∈PORE2};

[0087] In the formula, BIT 边缘 BIT1 represents the set of edge-filling media. 边缘 BIT2 represents the edge-filling medium set 1. 边缘 Represents edge-filling medium set 2, TBtarget边缘 Indicates the content of edge-filling medium, px 边缘 This represents an element in the initial three-dimensional digital core model filled with pore edges of the medium; dist represents the Euclidean distance operator; PDIST j Indicated in the second pore set PORE 边缘 The second Euclidean distance set consisting of elements with a medium Euclidean distance of j, PDSB j PDIST represents the set of Euclidean distances j. j The proportion of medium elements in the second pore set; PDIST i+1 Indicated in the second pore set PORE 边缘 The second Euclidean distance set consisting of elements with a median Euclidean distance of i+1, PDSB i+1 PDIST represents the set with Euclidean distance i+1. j+1 The proportion of the element quantity in the second pore set; ∪ j＝1 i PDIST j Let PDIST1, PDIST2, ..., PDIST be the set of PDIST. i The union of .

[0088] Furthermore, the construction of a medium-filled three-dimensional digital core model based on a three-dimensional digital core model with the medium filling the pore center and a three-dimensional digital core model with the medium filling the pore edge includes:

[0089] First, based on the three-dimensional digital core model with the medium filling the center of the pores, the center is filled to a predetermined range. Then, based on the three-dimensional digital core model with the medium filling the edge of the pores, the edge filling is completed to obtain the three-dimensional digital core model with the medium filling.

[0090] Furthermore, the construction of a medium-filled three-dimensional digital core model based on a three-dimensional digital core model with the medium filling the pore center and a three-dimensional digital core model with the medium filling the pore edge also includes:

[0091] First, based on the three-dimensional digital core model with the medium filling the pore edges, the edge filling is carried out to a predetermined range. Then, based on the three-dimensional digital core model with the medium filling the pore center, the center filling is completed to obtain the three-dimensional digital core model with the medium filling.

[0092] The present invention also provides a device for constructing a three-dimensional digital core model filled with a medium, the device comprising,

[0093] The acquisition module is used to obtain the type, distribution pattern, and content of the core filling medium based on two-dimensional core images;

[0094] The first construction module is used to construct a three-dimensional digital core model with the medium filling the center of the pores and a three-dimensional digital core model with the medium filling the edge of the pores, based on the distribution pattern and content of the filling medium and using a three-dimensional digital core model that identifies the skeleton and pores.

[0095] The second construction module is used to construct a medium-filled three-dimensional digital core model based on a three-dimensional digital core model with the medium filling the pore center and a three-dimensional digital core model with the medium filling the pore edge.

[0096] The technical effects and advantages of this invention are as follows:

[0097] 1. To address the issue of non-Arch phenomenon in resistivity logging caused by media filling, and the relatively limited research on rock physical properties related to media filling, a three-dimensional digital core model of media filling is constructed to simulate rock physical properties.

[0098] 2. In view of the difficulties in the experimental process of studying the effects of the dissolution of the filling medium on the density, porosity, acoustic waves and resistivity of rocks before and after the dissolution of the filling medium in the laboratory, compared with the traditional laboratory rock physics experiments, the three-dimensional numerical simulation based on the rock core is more efficient, can be reused, and is convenient for studying the influence of single micro-factors. Rock physics numerical simulation is used to calculate the elastic parameters and resistivity properties of rocks.

[0099] 3. To address the difficulty in determining the relative strength of logging response changes caused by lithofacies versus those caused by porous filling media, which complicates quantitative logging evaluation, a three-dimensional digital core model with different lithologies and filling media was constructed.

[0100] 4. For filling media with different filling methods and different filling volumes, based on the initial three-dimensional digital core, a three-dimensional digital core model of the medium with different filling methods and different filling volumes is constructed in the core pore space using methods based on image Euclidean distance map and image morphology.

[0101] Other features and advantages of the invention will be set forth in the following description, and will be apparent in part from the description, or may be learned by practicing the invention. The objects and other advantages of the invention may be realized and obtained by means of the structures pointed out in the description and the drawings. Attached Figure Description

[0102] Figure 1 This is a flowchart of the analysis method in a specific embodiment of the present invention;

[0103] Figure 2 To study the structural diagram of the filling medium in the pore center of the layer, the filling method and filling content of the filling medium in the pore space are obtained based on the diagram;

[0104] Figure 3 To study the structural diagram of the filling medium in the pore edge, the filling method and filling content of the filling medium in the pore space are obtained based on the diagram.

[0105] Figure 4 A visualization of the skeleton and pore sets in a three-dimensional digital core model for identifying the skeleton and pores;

[0106] Figure 5 A perspective view of the pore set distribution in a three-dimensional core model for identifying the skeleton and pores;

[0107] Figure 6 A three-dimensional digital core model diagram after the medium has filled the pore center;

[0108] Figure 7 A perspective view of the distribution of pore set and filling medium set in a three-dimensional core model after the medium has filled the pore centers;

[0109] Figure 8 A perspective view of the distribution of the filling medium in a three-dimensional core model after the medium has filled the pore centers;

[0110] Figure 9 Perspective view of the pore aggregate distribution in a three-dimensional core model after the medium has filled the pore centers;

[0111] Figure 10 The three-dimensional digital core model after the medium is filled in the pore center serves as the initial three-dimensional digital core model diagram for subsequent medium filling at the pore edge;

[0112] Figure 11 The perspective view of the pore set distribution in the initial three-dimensional digital core model after the medium is filled in the pore center, which is the three-dimensional digital core model after the medium is filled in the pore edge;

[0113] Figure 12 A three-dimensional digital core model diagram showing that the medium first fills the pore center and then fills the pore edge;

[0114] Figure 13 A perspective view of the filling medium and pore set distribution of a three-dimensional digital core model in which the medium first fills the pore center and then fills the pore edge;

[0115] Figure 14 A perspective view of the distribution of the filling medium in a three-dimensional digital core model, in which the medium is first filled at the center of the pores and then at the edge of the pores.

[0116] Figure 15 A perspective view of the pore set distribution of a three-dimensional digital core model in which the pore center of the medium is filled first and then the pore edge of the medium is filled;

[0117] Figure 16 A pore network model diagram for identifying the pore set in a three-dimensional digital core model to identify the skeleton and pores;

[0118] Figure 17 A pore network model diagram of the pore set in a three-dimensional digital rock core after filling the pore centers of the medium;

[0119] Figure 18 A pore network model diagram of the pore set in a three-dimensional digital rock core after the pore center of the medium is filled first and then the pore edge of the medium is filled;

[0120] Figure 19 To identify the pore network model of the pore set in the three-dimensional digital rock core, the pore size distribution curve and volume average pore size were obtained. The pore size distribution curve and volume average pore size after filling the pore center of the medium, as well as the pore size distribution curve and volume average pore size after filling the pore center and edge of the Jiezi pores.

[0121] Figure 20 Porosity and resistivity of 3D digital cores for identifying pores and skeleton, as well as resistivity comparison diagrams of the pore center and pore edge after filling. Detailed Implementation

[0122] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0123] To address the shortcomings of existing technologies, this invention discloses a method for constructing a three-dimensional digital core model filled with a medium, the method comprising:

[0124] Step 1: Based on two-dimensional core images, identify the type and distribution of the backfill medium, and determine the content of the backfill medium through image processing.

[0125] Step 2: Based on the distribution morphology and content of the filling medium, a three-dimensional digital core model of the medium filling the pore center is constructed using a three-dimensional digital core model that identifies the skeleton and pores.

[0126] Step 3: Based on the distribution morphology and content of the filling medium, a three-dimensional digital core model of the medium filling the pore edges is constructed using a three-dimensional digital core model that identifies the skeleton and pores.

[0127] Step 4: Based on the three-dimensional digital core model with the medium filling the pore center and the three-dimensional digital core model with the medium filling the pore edge, construct a three-dimensional digital core model with medium filling.

[0128] Combination Figure 1 In a specific embodiment of the present invention, the construction of a three-dimensional digital core model in which the medium fills the pore center in step 2 includes establishing a three-dimensional digital core model that identifies the skeleton and pores by performing fixed threshold segmentation on X-ray CT scan images, specifically:

[0129] A three-dimensional digital core model for identifying the skeleton and pores was established by using X-ray CT scan images for fixed threshold segmentation, or a three-dimensional digital core model for identifying the skeleton and pores was constructed using numerical reconstruction algorithms based on two-dimensional images of the skeleton and pores, serving as the initial three-dimensional digital core model for filling the pore centers of the medium. Based on the distribution morphology, content, and initial three-dimensional digital core model of the filling medium, a three-dimensional digital core model of the medium filling the pore centers was constructed using image morphology methods.

[0130] Furthermore, step S2 specifically includes the following sub-steps:

[0131] S201. Obtain the type of filling medium, the content of filling medium, and the maximum diameter dmax of the core pores identified in step S1;

[0132] S202. Based on the three-dimensional digital core model that identifies the skeleton and pores, the initial three-dimensional digital core model for filling the pore center of the medium is used as the initial three-dimensional digital core model. The phase state function of each element in the initial three-dimensional digital core model filled with the pore center of the medium is obtained, and the first skeleton set and the first pore set of the initial three-dimensional digital core model filled with the pore center of the medium are divided. The first skeleton set represents the set of solid elements in the three-dimensional digital core, and the first pore set represents the set of first pore elements in the three-dimensional digital core. The number of elements in the first pore set is obtained.

[0133] S203. In the first pore set representing the pore space, the image morphology-based method uses spherical structural elements of different diameters to perform opening operations in the pore set. The opening operation results set of spherical structural elements of different diameters is obtained through the opening operation, which is the first opening operation result set. The number of first pore elements in the first opening operation result set and the proportion of the number of first pore elements in the first opening operation result set in the first pore set are obtained.

[0134] S204. Based on the first opening operation result set, obtain the difference set of the opening operation result sets of two adjacent diameter sphere structural elements, and perform Euclidean distance calculation on the difference set. According to the Euclidean distance, divide the difference set into different Euclidean distance sets, and obtain the number of pore elements in the corresponding Euclidean distance set and the proportion of the number of pore elements in the pore set of the initial three-dimensional digital core model filled with the pore center of the medium.

[0135] S205. Based on the first opening operation result set, the number of first pore elements in the first opening operation result set, the proportion of the number of first pore elements in the first opening operation result set in the first pore set, the difference set of the opening operation result sets of two adjacent diameter spherical structural elements in the first opening operation result set, the Euclidean distance set in the difference set, the number of pore elements in the Euclidean distance set and the proportion of the corresponding number of pore elements in the first pore set of the initial three-dimensional digital core model filled with medium pore center, and the content of filling medium, construct a three-dimensional digital core model filled with medium pore center.

[0136] In a specific embodiment of the present invention, in step S202, any element in the three-dimensional digital core model set is represented by px. 中心 The phase function `phase` is used to indicate whether a point is part of the skeleton or a pore, where `phase(px)` represents the phase state. 中心 ) = 0 means px 中心 The element is a skeleton, phase(px) 中心 ) = 1 represents px 中心 The element is a pore, phase(px) 中心 ) = 3 indicates that the element is filling the medium; CORE1 represents the first initial three-dimensional digital core set, used to identify the three-dimensional digital core model of the skeleton and pores; the first initial three-dimensional digital core set CORE1 uses the first skeleton set MATRIX 中心 and the first pore set PORE 中心 The union representation, i.e., CORE1 = MATRIX 中心 ∪PORE 中心 Furthermore, the core skeleton assembly MATRIX 中心 and core pore set PORE 中心 Represented as:

[0137] MATRIX 中心 ={px 中心 |phase(px 中心 )=0};

[0138] PORE 中心 ={px 中心 |phase(px 中心 )=1};

[0139] Calculate the first pore set PORE 中心 The number of elements in nump 中心 .

[0140] In a specific embodiment of the present invention, step S203 mainly involves performing an opening operation on the first pore set of the initial three-dimensional digital core model for spherical structural elements of different diameters, the number of elements in the set, and the proportion of the number of elements in the set to the total number of elements in the pore space set.

[0141] ① Determine the maximum diameter of the core pores as dmax;

[0142] ② Select the first pore set PORE from the set CORE1, which is the initial three-dimensional digital core model filled with the pore centers of the medium. 中心 and the first pore set PORE 中心 The number of elements in nump 中心 .

[0143] ③ Select the diameter of the sphere structural element in the image morphology opening operation as i, and the range of i is from 1 to dmax, i.e. i = 1, 2, ..., dmax.

[0144] ④ For the first pore PORE 中心 Each element in px 中心 A morphologically based method is used to perform an opening operation OPEN(px,i) on spherical structural elements with diameter i. OPEN(px,i) = 0 indicates that the element does not meet the opening operation requirements for spherical structural elements with diameter i in the pore space set, while OPEN(px,i) = 1 indicates that the element meets the opening operation requirements for spherical structural elements with diameter i in the pore space set. This yields the set MORPH of opening operation results for spherical structural elements with diameter i. i This forms the first set of results from the opening operation.

[0145] MORPH i ={px|OPEN(px,i)=1},(i=1,2,…,dmax);

[0146] ⑤ Based on the first luck set MORPH i Calculate the number of elements in a set using the Onum function. i ,(i=1,2,…,dmax);

[0147] ⑥ Based on the first open operation result set MORPH i MORPH i The number of elements in the set (Onum) i ,(i=1,2,…,dmax), calculate MORPHi The number of elements in the set is in the first pore set PORE 中心 The proportion of OSBi in the total.

[0148] OSBi = Onum i / nump 中心 , (i = 1, 2, ..., dmax);

[0149] In a specific embodiment of the present invention, in step S204,

[0150] ① Select the first set of results from the opening operation, MORPH i The number of elements in the set, Onum i and MORPH i The number of elements in the set is in the first pore set PORE 中心 The proportion of OSBi in (i = 1, 2, ..., dmax-1).

[0151] ② Select the set of opening operation results of the sphere structuring element with diameter i+1 in the first opening operation set, MORPH i+1 The number of elements in the set, Onum i+1 and MORPH i+1中心 The number of elements in the set is in the first pore set PORE 中心 The proportion of OSB i+1 , (i = 1, 2, ..., dmax-1).

[0152] ③ Based on the set MORPH in the result of the first opening operation i MORPH i+1 DMORPH is obtained through difference operations. i (i = 1, 2, ..., dmax-1).

[0153] DMORPH i =MORPH i -MORPH i+1 ;

[0154] ④ Calculate the difference set DMORPH i any element in px 中心 To the set MORPH i+1 The Euclidean distance is j = dist(px) 中心 ), (i = 1, 2, ..., dmax-1);

[0155] ⑤ In the difference set DMORPH i In (i = 1, 2, ..., dmax-1), based on the Euclidean distance j = dist(px) 中心(j = 1, 2, ..., n, where n is the largest Euclidean distance in the set) The difference set DMORPH i The set is partitioned into a series of Euclidean distance sets, DIST. j .

[0156] DIST j ={px 中心 |dist(px 中心 ) = j, and px 中心 ∈DMORPH i-1},(j=1,2,...,n);

[0157]

[0158] ⑥ Calculate the difference set DMORPH i Central Euclidean distance set DIST j The number of elements in Dnum j中心 ;

[0159] ⑦ Calculate the difference set DMORPH i Central Euclidean distance set DIST j The number of elements in Dnum j First pore set PORE 中心 The proportion of DSB j ;

[0160] DSB j =Dnum j / nump 中心 j = 1, 2, ..., n;

[0161] In a specific embodiment of the present invention, in step S205, based on the set MORPH... i and the percentage of elements in this set (OSB) i and the difference set of adjacent opening operations, DMORPH i-1 Central Euclidean distance set DIST i The number and proportion of pores in DSB i and filling medium content TBtarget 中心 This allows the medium to gradually fill the pore space from the pore center (i = 2, 3, ..., dmax);

[0162] S2051, if OSB dmax >TBtarget 中心 Then MORPH dmax All pores in the set are transformed into filling media, establishing a media set BIT.

[0163]

[0164] The number of elements in the filling medium set is Onum dmax .

[0165] In the formula, TBtarget 中心 Indicates the content of the first filling medium, MORPH dmax Onum represents the set of elements in the first set of pores that satisfy the result of the opening operation of a sphere with diameter dmax. dmax MORPH dmax Number of elements in the set; OSB dmax Represents Onum dmax First pore set PORE 中心 The proportion of the number of elements;

[0166] Complete the center filling of the medium pores and obtain the center-filled medium assembly (BIT). 中心 The next step is to proceed to step S2053;

[0167] S2052, If OSB i <TBtarget 中心 <OSB i-1 i = 2, 3, 4, ..., d max hour;

[0168] ① First, open the set of structural elements of the sphere with diameter i, MORPH. i All elements in the set are transformed into the filling medium set BIT1. 中心 ,

[0169] That is: BIT1 中心 =MORPH i MORPH i ={px 中心 |OPEN(px 中心 ,i)=1};

[0170] In the formula, BIT1 中心 This indicates the filling medium set 1, MORPH i OSB represents the set of elements in the first set of pores that satisfy the result of the opening operation of a sphere with diameter i. i MORPH i The elements in the first pore set PORE 中心 Percentage of elements; px 中心 OPEN(px) represents the element in the initial three-dimensional digital core model filled with the pore centers of the medium. 中心 ,i)=1 indicates that the center of px satisfies the opening operation result of the sphere structure element with diameter i.

[0171] ②If DSB1>TBtarget 中心 -OSB i (i = 2, 3, 4, ..., d) max Then, randomly select (nump) from the DIST1 set. 中心 ×TBtarget 中心 -OSB i The pores are transformed into filling media, establishing a filling media set BIT2. 中心 ,

[0172] BIT2 中心 ={px 中心 Randomly select px 中心 , px 中心 ∈DIST1};

[0173] DMORPH i-1 =MORPH i-1 -MORPH i ;

[0174] DIST1 = {px 中心 |dist(px 中心 ) = 1, and px 中心 ∈DMORPH i-1};

[0175] Complete the filling of the pore center of the medium and obtain the set of pore center fillings (BIT). The next step is to proceed to step S2053.

[0176] OSB i <TBtarget 中心 <OSB i-1 And DSB1>TBtarget 中心 -OSB i (i = 2, 3, 4, ..., d) max When ), the filling medium set is:

[0177] BIT 中心 =BIT1 中心 ∪BIT2 中心 ;

[0178] BIT1 中心 =MORPH i MORPH i ={px 中心 |OPEN(px 中心 ,i)=1};

[0179] BIT2 中心 ={px 中心 Randomly select px 中心 , px中心 ∈DIST1};

[0180] DMORPH i-1 =MORPH i-1 -MORPH i ;

[0181] DIST1 = {px 中心 |dist(px 中心 ) = 1, and px 中心 ∈DMORPH i-1};

[0182] In the formula, BIT 中心 BIT1 represents the centrally filled medium set. 中心 Indicates center-filled medium set 1, BIT2 中心 Indicates center-filled medium set 2, MORPH i OSB represents the set of elements (i = 2, 3, ..., dmax) in the first set of pores that satisfy the opening operation of a sphere with diameter i. i MORPH i The elements in the first pore set PORE 中心 The proportion of elements in the total number; MORPH i-1 OSB represents the set of elements (i = 2, 3, ..., dmax) resulting from the opening operation of spheres with diameter i-1 in the first pore set. i-1 MORPH i-1 The elements in the first pore set PORE 中心 Percentage of elements; px 中心 This represents an element in the initial three-dimensional digital core model filled with pores in the medium. OPEN represents the opening operator, and i and i-1 represent the diameters of the spherical structural elements, respectively. OPEN(px 中心 ,i)=1 indicates that the center of px satisfies the opening operation result of the sphere structure element with diameter i (i=2,3,…,dmax).

[0183] MORPH i DMORH represents the set of elements in the first set of pores that satisfy the result of the opening operation of a sphere structure element with diameter i; i-1 DIST1 represents the difference set of the results of the opening operation between sphere structural elements with diameters i-1 and i in the first pore set based on image morphology; DIST1 represents DMORPH. i-1 The Euclidean distance in the set is 1, i.e., dist(px) 中心 The set consisting of elements where ) = 1. DSB1 represents the set of elements in DIST1 in the first pore set PORE. 中心The percentage of elements; dist represents the Euclidean distance operator.

[0184] ③If OSB i <TBtarget 中心 <OSB i-1 (i = 2, 3, 4, ..., d) max And OSB i +∑ j＝1 i DSB j <TBtarget 中心 <OSB i +∑ j＝1 i+1 DSB j >TBtarget 中心 -OSB i ;

[0185] First, open the set of structuring elements of the sphere with diameter i, MORPH. i All elements in the set are transformed into the filling medium set BIT1′. 中心 ,

[0186] That is: BIT1' 中心 =MORPH i MORPH i ={px 中心 |OPEN(px 中心 ,i)=1};

[0187] In the formula, BIT1′ 中心 This indicates the filling medium set 1, MORPH i OSB represents the set of elements in the first set of pores that satisfy the result of the opening operation of a sphere with diameter i. i MORPH i The elements in the first pore set PORE 中心 Percentage of elements; px 中心 OPEN(px) represents the element in the initial three-dimensional digital core model filled with the pore centers of the medium. 中心 ,i)=1 indicates that the center of px satisfies the opening operation result of the sphere structure element with diameter i.

[0188] Secondly, the first open operation set MORPH i-1 MORPH i DMORPH difference set i-1 All elements in the set whose Euclidean distance is less than or equal to i fill the medium set BIT2′. 中心 ,Right now:

[0189]

[0190] DIST j ={px 中心 |dist(px 中心 ) = j, and px 中心 ∈DMORPH i-1};

[0191] DMORPH i-1 =MORPH i-1 -MORPH i ;

[0192] In the formula, BIT2′ 中心 ∪ represents the set of centrally filled media 2'; ∪ represents the union of sets 2' and 2''. j＝1 i DIST j Let DIST1, DIST2, ..., DIST be the set. i Union of; MORPH i OSB represents the set of elements in the first set of pores that satisfy the result of the opening operation of a sphere with diameter i; i MORPH i The elements in the first pore set PORE 中心 The proportion of elements in the total number; MORPH i-1 OSB represents the set of elements in the first set of pores that satisfy the opening operation of a sphere with diameter i-1; i-1 MORPH i-1 The elements in the first pore set PORE 中心 Percentage of elements; px 中心 OPEN represents the element in the initial three-dimensional digital core model filled with the pore center of the medium; OPEN represents the opening operator; i, i-1 represent the diameters of the spherical structural element, respectively; OPEN(px 中心 , i) = 1 indicates that the center of px satisfies the opening operation result of a sphere structuring element with diameter i; DMORH i-1 DIST represents the difference set of the results of the opening operation between sphere structural elements with diameters i-1 and i in the first set of pores, based on image morphology; j Indicates DMORPH i-1 The set of elements whose Euclidean distance is j; DSB j DIST j The elements in the first pore set PORE 中心 The percentage of the total number of elements.

[0193] Finally, in the first open set MORPH i-1 MORPH iDMORPH difference set i-1 The set whose Euclidean distance is equal to i+1 is randomly selected (TBtarget) 中心 -OSB i -∑ j＝1 i DSB j )×nump 中心 Each element is transformed into a filling medium set BIT3′ 中心 ,Right now:

[0194] BIT3' 中心 ={px 中心 Randomly select px 中心 , px 中心 ∈DIST j+1};

[0195] DIST j+1 ={px 中心 |dist(px 中心 ) = j + 1, and px 中心 ∈DMORPH i-1};

[0196] DMORPH i-1 =MORPH i-1 -MORPH i ;

[0197] In the formula, BIT3′ 中心 Indicates the central filling medium set 3'; px 中心 This represents an element in the initial three-dimensional digital core model filled with pore centers of the medium; MORPH i OSB represents the set of elements in the first set of pores that satisfy the result of the opening operation of a sphere with diameter i; i MORPH i The elements in the first pore set PORE 中心 The proportion of elements in the total number; MORPH i-1 OSB represents the set of elements in the first set of pores that satisfy the result of the opening operation of a sphere with diameter i-1. i-1 MORPH i-1 The elements in the first pore set PORE 中心 The percentage of elements in the total number; i, i-1 represent the diameters of the spherical structural elements; DMORH i-1 This represents the difference set of the results of the opening operation between sphere structural elements of diameter i-1 and i in the first pore set, based on image morphology. j Indicates DMORPH i-1 The set of elements whose Euclidean distance is j; DSBj DIST j The elements in the first pore set PORE 中心 nump of element count 中心 The proportion in; MORPH i-1 DMORH represents the set of elements in the first set of pores that satisfy the opening operation of spherical structural elements with diameter i-1; i-1 This represents the difference set of the results of the opening operation between sphere structural elements of diameter i-1 and i in the first pore set, based on image morphology. j+1 Indicates DMORPH i-1 The set consisting of elements whose Euclidean distance is j+1. DSB j+1 DIST j+1 The elements in the first pore set PORE 中心 The percentage of elements; dist represents the Euclidean distance operator; dist(px 中心 ) indicates that the element px is calculated. 中心 Euclidean distance.

[0198] That is, when OSB i <TBtarget 中心 <OSB i-1 And OSB i +∑ j＝1 i DSB j <TBtarget 中心 <OSB i +∑ j＝1 i+1 At that time, the central filling medium set is:

[0199] BIT 中心 =BIT1' 中心 UBIT2' 中心 ∪BIT3' 中心 ;

[0200] BIT1' 中心 =MORPH i MORPH i ={px 中心 |OPEN(px 中心 ,i)=1};

[0201]

[0202] BIT3' 中心 ={px 中心 Randomly select px 中心 , px 中心 ∈DIST j+1};

[0203] DMORPH i-1 =MORPH i-1 -MORPH i ;

[0204] DIST j ={px 中心 |dist(px 中心 ) = j, and px 中心 ∈DMORPH i-1};

[0205] DIST j+1 ={px 中心 |dist(px 中心 ) = j + 1, and px 中心 ∈DMORPH i-1}(i=2,3,…,dmax);

[0206] In the formula, BIT1′ 中心 Indicates center-filling medium set 1'; BIT2' 中心 Indicates the central filling medium set 2'; BIT3' 中心 ∪ represents the set of centrally filled media 3'; ∪ represents the union, ∪ j＝1 i DIST j Denotes sets DIST1, DIST2, ..., DIST i Union of; px 中心 This represents an element in the initial three-dimensional digital core model filled with pores in the medium; OPEN(px 中心 , i) = 1 indicates that the center of px satisfies the opening operation result of a sphere structure element with diameter i; MORPH i MORPH i-1 OSB represents the set of elements in the first set of pores that satisfy the opening operation of spheres with diameters of i and i-1. i MORPH i The elements in the first pore set PORE 中心 Percentage of elements; OSB i-1 MORPH i-1 The elements in the first pore set PORE 中心 The proportion of elements in the total number; DMORH i-1 MORPH represents the set of results of opening operations on sphere structural elements with diameters i-1 and i in the first set of pores, based on image morphology. i MORPH i-1 The difference set (i = 2, 3, ..., dmax). DIST jIndicates DMORPH i-1 The set of elements whose Euclidean distance is j; DSB j DIST j The elements in the first pore set PORE 中心 The proportion of elements in the total number; DIST j+1 Indicates DMORPH i-1 The set of elements whose Euclidean distance is j+1; DSB j+1 DIST j+1 The elements in the first pore set PORE 中心 nump of element count 中心 The proportion in; dist represents the Euclidean distance operator; dist(px 中心 ) indicates that the element px is calculated. 中心 Euclidean distance.

[0207] This completes the process of filling the center of the medium pores.

[0208] S2053, Establishing a central filling medium set (BIT) 中心 Subsequently, the three-dimensional digital core assembly CORE 中心 Represented as the skeleton set MATRTIX 中心 First pore set PORE 中心 BIT (Body Filler) 中心 The union of, i.e.:

[0209] CORE 中心 =MATRIX 中心 UPORE 中心 ∪BIT 中心 ;

[0210] In the formula, MATRIX 中心 PORE represents the first skeleton set. 中心 CORE represents the first set of pores. 中心 A three-dimensional digital core set representing a pore space model where the medium gradually fills the pore center, where U represents the union, and the first skeleton set is MATRIX. 中心 Represented as:

[0211] MATRIX 中心 ={px 中心 |phase(px 中心 )=0};

[0212] First Pore Set PORE 中心 Represented as:

[0213] PORE 中心 ={px 中心|phase(px 中心 )=1};

[0214] Central Medium Filling Set (BIT) 中心 Represented as:

[0215] BIT 中心 ={px 中心 |phase(px 中心 )=3};

[0216] In the formula, px 中心 This represents an element in the initial three-dimensional digital core model filled with pore centers of the medium, where phase represents the phase state function, phase(px 中心 ) = 0 means px 中心 The element is a skeleton, phase(px) 中心 ) = 1 represents px 中心 The element is a pore. phase(px) 中心 ) = 3 represents px 中心 The element is the filling medium.

[0217] In a specific embodiment of the present invention, step 3, constructing a three-dimensional digital core model with the medium filling the pore edges, includes:

[0218] Based on the three-dimensional digital core model that identifies the skeleton and pores, an initial three-dimensional digital core model of the medium filling the pore edges is constructed. Based on the distribution morphology of the filling medium, the content of the filling medium, and the initial three-dimensional digital core model of the medium filling the pore edges, a three-dimensional digital core model of the medium filling the pore edges is constructed using the image Euclidean distance map method.

[0219] Specifically, it includes the following sub-steps:

[0220] S301. Obtain the filling medium type and filling medium content TBtarget obtained in step S1. 边缘 ;

[0221] S302. The three-dimensional digital core model for identifying the skeleton and pores is used as the initial three-dimensional digital core model for filling the pore edges of the medium. Based on the phase function of each element in the three-dimensional digital core model, the second skeleton set and the second pore set are divided, and the number of elements in the second pore set is obtained.

[0222] S303. In the identified pore set of the pore space, the Euclidean distance method is used to calculate the Euclidean distance from the element in the second pore set to the second skeleton set; based on the Euclidean distance from the element in the second pore set to the second skeleton set, the second pore set is divided into multiple second Euclidean distance sets, and the number of elements in each second Euclidean distance set and the proportion of the number of elements in each Euclidean distance set in the second pore set are obtained.

[0223] S304. Based on multiple sets of second Euclidean distances, the proportion of the number of elements in each Euclidean distance set in the second pore set, and the content of the filling medium, a three-dimensional digital core model of the filling medium edge is constructed.

[0224] In a specific embodiment of the present invention, in step S302, the three-dimensional digital core model identifying the skeleton and pores is used as the initial three-dimensional digital core model for filling the pore edges of the medium, denoted by the second initial core set CORE2, wherein:

[0225] CORE2 = MATRIX 边缘 ∪PORE 边缘 ;

[0226] In the formula, MATRIX 边缘 PORE represents the second skeleton set. 边缘 Let U denote the second set of pores, and U denote the union of the second skeleton set, MATRIX. 边缘 Represented as:

[0227] MARTIX 边缘 ={px 边缘 |phase(px 边缘 )=0};

[0228] Second pore set PORE 边缘 Represented as:

[0229] PORE 边缘 ={px 边缘 |phase(px 边缘 )=1},

[0230] Second pore set PORE 边缘 The number of elements in nump 边缘 .

[0231] In the formula, px 边缘 This represents the elements in the initial three-dimensional digital core set filled with pore edges of the medium, where phase represents the phase state function, phase(px 边缘 ) = 0 means px 边缘 The element is a skeleton, phase(px) 边缘 ) = 1 represents px边缘 The element is porosity.

[0232] In a specific embodiment of the present invention, in step S303,

[0233] ① The second pore set PORE of the initial three-dimensional digital core model CORE2, which identifies the skeleton and pores, and is used as the pore edge filling model of the medium. 边缘 Based on the calculated Euclidean distance dist, the pore set PORE is... 边缘 Let PDIST be the set of second Euclidean distances with Euclidean distance i. i边缘 The union of, i.e.:

[0234]

[0235] The minimum Euclidean distance is 1, the maximum Euclidean distance is n, and the second pore set PORE 边缘 Mid-range skeleton ensemble MATRIX 边缘 The elements with a Euclidean distance of i form the second Euclidean distance set PDIST. i ;

[0236] ② Fill the pore edges of the medium with the pore set PORE of the initial three-dimensional digital core model CORE2. 边缘 PDIST, a second Euclidean distance set based on Euclidean distance partitioning i Calculate the n sets of second Euclidean distances PDIST in (i = 1, 2, ..., n). i The number of elements in each set is divided into PDnum. i If i = 1, 2, ..., n, then we have

[0237] ③ Based on the filling of the pore edges of the medium with the second pore set PORE of the three-dimensional digital core model CORE2 边缘 The number of elements nump 边缘 With n sets of second Euclidean distances PDIST i Number of elements PDnum i Calculate the corresponding second Euclidean distance set PDIST i Number of elements in PDnum i Second pore set PORE 边缘 The proportion of PDSB in i ;

[0238]

[0239] In a specific embodiment of the present invention, in step S304, the pore set PORE of the initial three-dimensional digital core model CORE2, which is filled with pore edges of the medium, is... 边缘The second Euclidean distance set PDIST i (i = 1, 2, ..., n) and the percentage of elements in this set PDSB i and filling medium content TBtarget 边缘 The medium is gradually filled from the edge with a Euclidean distance of 1 towards the pore space with an increasing Euclidean distance. The specific steps are as follows:

[0240] S3041, If ​​PDSB 1边缘 >TBtarget 边缘 Then, nump is randomly selected from the PDIST1 set. 边缘 ×TBtarget 边缘 The pores are transformed into filling media, and a three-dimensional digital core model of the filling media edge is established, in which the edge filling media set BIT is formed. 边缘 The next step is step S3043.

[0241] BIT 边缘 ={px 边缘 Randomly select px 边缘 , px 边缘 ∈PDIST1};

[0242] PDIST1 = {px 边缘 |dist(px 边缘 ) = 1, and px 边缘 ∈PORE 边缘 The number of elements in the filling medium set is nump 边缘 ×TBtarget 边缘 ;

[0243] In the formula, BIT 边缘 TBtarget represents the edge-filled media set. 边缘 Indicates the content of the edge-filling medium, nump 边缘 PORE represents the second set of pores filled at the edge. 边缘 The number of elements, px 边缘 This represents an element in the initial three-dimensional digital core model filled with pore edges of the medium. `dist` represents the Euclidean distance operator, and `PDIST1` represents the element in the second pore set `PORE`. 边缘 The second Euclidean distance set consists of elements with a distance of 1 in the Euclidean distance set. PDSB1 represents the proportion of the number of elements in the second Euclidean distance set with a distance of 1 in the second pore set.

[0244] S3042. If the pore edges of the medium are filled with the pore set PORE of the three-dimensional digital core model CORE2... 边缘 The set of all second Euclidean distances PDIST with mean Euclidean distance dist≤i jThe percentage of elements in the sequence (j = 1, 2, ..., i) is less than TBtarget. 边缘 The set of all elements PDIST with dist≤i+1 j The sum of the percentages of elements in the sequence (j = 1, 2, ..., i, i+1) is greater than TBtarget. 边缘 hour,

[0245] Right now hour,

[0246] Then ① fill the pore edges of the medium with the pore set PORE of the three-dimensional digital core model CORE2. 边缘 The medium Euclidean distance is the set of all second Euclidean distances less than or equal to i, PDIST. j The elements (j = 1, 2, ..., i) are used as filling media to establish the filling medium set BIT1. 边缘 ;

[0247] Right now

[0248] In the formula, BIT1 边缘 Indicates edge-filling medium set 1, TBtarget 边缘 Indicates the content of edge-filling medium, px 边缘 Represents the second pore set PORE 边缘 In the element, dist represents the Euclidean distance operator, dist(px 边缘 ) represents the pore set PORE 边缘 medium element px 边缘 Euclidean distance to the skeleton set; PDIST j Indicated in the second pore set PORE 边缘 The second Euclidean distance set consisting of elements with a medium Euclidean distance of j, PDSB j This indicates the number of elements with Euclidean distance j in the Euclidean distance set in the second pore set PORE. 边缘 The proportion in. ∪ j＝1 i PDIST j Let PDIST1, PDIST2, ..., PDIST be the set of PDIST. i The union of ; ∪ represents the union of .

[0249] ② In the set PDIST of dist = i+1 i+1 Random selection Each porous element is transformed into a filling medium, namely BIT2. 边缘 ={px 边缘 Randomly select px 边缘 , px 边缘 ∈PDIST i+1};

[0250] ③ The filling medium content is TBtarget 边缘 In a three-dimensional digital core model filled with pore edges, the collection of edge-filling media (BIT) 边缘 For: BIT 边缘 =BIT1 边缘 ∪BIT2 边缘 This completes the filling of the pore center of the medium and obtains the set of BITs of the edge medium filling. 边缘 ;

[0251] when hour,

[0252] The filling medium set is as follows:

[0253] BIT 边缘 =BIT1 边缘 ∪BIT2 边缘 ;

[0254]

[0255] BIT2 边缘 ={px 边缘 Randomly select px 边缘 , px 边缘 ∈PDIST i+1};

[0256] PDIST i+1 ={px 边缘 |dist(px 边缘 ) = i + 1, and px 边缘 ∈PORE 边缘};

[0257] In the formula, BIT 边缘 BIT1 represents the set of edge-filling media. 边缘 BIT2 represents the edge-filling medium set 1. 边缘 Represents edge-filling medium set 2, TBtarget 边缘 Indicates the content of edge-filling medium, px 边缘 PORE represents the second pore set in a three-dimensional digital core model where the pore edges of the medium are filled. 边缘 In the elements, OPEN represents the opening operator, dist represents the Euclidean distance operator, and dist(px) represents the distance operator. 边缘 ) represents the pore set PORE 边缘 medium element px 边缘 Euclidean distance to the skeleton set; DIST j Let j represent the second Euclidean distance set of elements with Euclidean distance j.

[0258] BIT1 边缘 PORE represents the second pore set of the three-dimensional digital core model CORE2, which is filled with pore edges of the medium. 边缘 In the set of Euclidean distances less than or equal to i, the union of the sets of Euclidean distances is BIT2. 边缘 The second pore set PORE represents the second initial three-dimensional digital core model CORE2. 边缘 In the equation, the set of Euclidean distances PDIST with Euclidean distance equal to i+1 is... i+1 The set is composed of randomly selected elements.

[0259] S3043. After establishing the medium-filled assemblies (BITs), the three-dimensional digital core assemblies (COREs) are created. 边缘 It can be expressed as:

[0260] CORE 边缘 =MATRIX 边缘 ∪PORE 边缘 ∪BIT 边缘 ;

[0261] BIT 边缘 ={px 边缘 |phase(px 边缘 )=3};

[0262] Wherein, CORE2 represents the initial three-dimensional digital core set filled with edge media, CORE 边缘 MATRIX represents a three-dimensional digital core assembly after the media edge is filled. 边缘 Represents the second skeleton set, PORE 边缘 Represents the second set of pores, BIT 边缘 Represents the set of edge-filling media, phase(px) 中心 ) = 3 indicates that the element fills the medium, px 边缘 The second pore set PORE in the initial three-dimensional digital model representing the filling of the pore centers of the medium 边缘 The elements in.

[0263] In a specific embodiment of the present invention, in step S4, three-dimensional core models with different filling methods of the medium are constructed by using different combinations of image morphology and image Euclidean distance map methods. Specifically, based on steps S2 and S3, according to the adjustment of the order of steps S2 and S3, a three-dimensional core model is constructed in which the medium first fills the pore edge and then fills the pore center, or a three-dimensional core model is constructed in which the medium first fills the pore center and then fills the pore edge.

[0264] Furthermore, the construction process of the three-dimensional data core model, in which the medium first fills the pore edges and then the pore center, is as follows:

[0265] S4A1. Based on the identification of the skeleton and pores, the initial three-dimensional digital core set CORE is represented as the union of the skeleton set MATRIX and the pore set PORE, that is, CORE = MATRIX∪PORE.

[0266] S4A2. Based on the content of the central filling medium, a three-dimensional data core model of the medium filling the center of the pores is constructed in the current pore set PORE using image morphology method, i.e., CORE=MATRIX∪PORE∪BIT1;

[0267] S4A3. Based on the three-dimensional digital core model of the medium filling the pore center, i.e., CORE=MATRIX∪PORE∪BIT1, according to the content of the edge filling medium, a three-dimensional digital core model of the medium filling the pore center and then filling the pore edge is constructed in the current pore set PORE using the image Euclidean distance method, i.e., CORE=MATRIX∪PORE∪BIT1∪BIT2.

[0268] Furthermore, the construction process of a three-dimensional digital core model, in which the medium first fills the pore edges and then the pore center, is as follows:

[0269] S4B1. Based on the three-dimensional digital core model that identifies the skeleton and pores, the initial three-dimensional digital core set CORE is represented as the union of the skeleton set MATRIX and the pore set PORE, i.e., CORE = MATRIX ∪ PORE.

[0270] S4B2. Based on the content of the edge-filling medium, a three-dimensional digital core model of the medium filling the pore edge is constructed in the current pore set PORE using the image Euclidean distance method, i.e., CORE=MATRIX∪PORE∪BIT2;

[0271] S4B3. Based on the three-dimensional digital core model of the medium filling the pore edges, i.e., CORE=MATRIX∪PORE∪BIT2, according to the amount of medium filling the center, a three-dimensional digital core model of the medium filling the pore center after filling the pore edges is constructed in the current pore set PORE using an image morphology-based method, i.e., CORE=MATRIX∪PORE∪BITI1∪BIT2.

[0272] This invention can also be configured such that step S5, based on the constructed three-dimensional digital core models with different filling methods of the medium, extracts pore structure parameters, simulates rock physical properties, and studies the influence of medium filling on pore space and rock physical properties. This includes the following sub-steps:

[0273] S501, Assignment of resistivity properties of each phase in three-dimensional digital core filled with medium;

[0274] S502. Finite element method simulation calculation of core resistivity;

[0275] S503. Extract the pore network model using the maximum sphere method to obtain pore structure parameters;

[0276] S504. Study the influence of dielectric filling on pore structure parameters and resistivity.

[0277] This innovative achievement provides a method for constructing three-dimensional digital cores with medium filling and simulating and analyzing rock physical properties. Its effects are mainly reflected in the following aspects:

[0278] 1. To address the issue of non-Arch phenomenon in resistivity logging caused by media filling, and the relatively limited research on rock physical properties related to media filling, a three-dimensional digital core model of media filling is constructed to simulate rock physical properties.

[0279] 2. In view of the difficulties in the experimental process of studying the effects of the dissolution of the filling medium on the density, porosity, acoustic waves and resistivity of rocks before and after the dissolution of the filling medium in the laboratory, compared with the traditional laboratory rock physics experiments, the three-dimensional numerical simulation based on the rock core is more efficient, can be reused, and is convenient for studying the influence of single micro-factors. Rock physics numerical simulation is used to calculate the elastic parameters and resistivity properties of rocks.

[0280] 3. To address the difficulty in determining the relative strength of logging response changes caused by lithofacies versus those caused by porous filling media, which complicates quantitative logging evaluation, a three-dimensional digital core model with different lithologies and filling media was constructed.

[0281] 4. For filling media with different filling methods and different filling volumes, based on the initial three-dimensional digital core, a three-dimensional digital core model of the medium with different filling methods and different filling volumes is constructed in the core pore space using methods based on image Euclidean distance map and image morphology.

[0282] The present invention also provides a device for constructing a three-dimensional digital core model filled with a medium, the device comprising,

[0283] The acquisition module is used to obtain the type, distribution pattern, and content of the core filling medium based on two-dimensional core images;

[0284] The first construction module is used to construct a three-dimensional digital core model with the medium filling the center of the pores and a three-dimensional digital core model with the medium filling the edge of the pores, based on the distribution pattern and content of the filling medium and using a three-dimensional digital core model that identifies the skeleton and pores.

[0285] The second construction module is used to construct a medium-filled three-dimensional digital core model based on a three-dimensional digital core model with the medium filling the pore center and a three-dimensional digital core model with the medium filling the pore edge.

[0286] The following are exemplary embodiments of the invention as defined by the claims and their equivalents, taken in conjunction with the accompanying drawings, to aid in a comprehensive understanding. The specific details described herein are to be considered exemplary only and not to limit the scope of the invention. Therefore, those skilled in the art can make various changes and modifications to the embodiments without departing from the scope and spirit of the invention.

[0287] Example 1

[0288] As a preferred embodiment of the present invention, please refer to the appendix to the specification. Figure 1 As shown in the figure, this embodiment discloses a method for constructing three-dimensional core data of medium filling and simulating and analyzing rock physical properties. The method includes the following steps:

[0289] S1. Based on two-dimensional core images, identify the type and distribution of the backfill medium, and determine the content of the backfill medium through image processing.

[0290] S2. Based on the distribution morphology and content of the filling medium, a three-dimensional digital core model with medium filling the center of the pores is constructed using image morphology methods by utilizing a three-dimensional digital core model that identifies the skeleton and pores.

[0291] S3. Based on the distribution morphology and content of the filling medium, a three-dimensional digital core model is constructed by using the three-dimensional digital core model that identifies the skeleton and pores, and the method of image Euclidean distance map is used to construct a three-dimensional digital core model in which the medium fills the pore edges.

[0292] S4. Based on the identification of skeleton and pores, three-dimensional digital core models with different filling methods of the medium are constructed by using different combinations of image morphology and image Euclidean distance map methods.

[0293] S5. Based on a three-dimensional digital core model with medium filling, extract pore structure parameters, simulate rock physical properties, and study the influence of medium filling on pore space and rock physical properties.

[0294] Example 2

[0295] As another preferred embodiment of the present invention, this embodiment is a further detailed description and supplement to the technical solution of the present invention based on the above-described embodiment 1. In this embodiment, Figure 2-9As shown, a three-dimensional digital core model with the medium filling the center is constructed. Based on the distribution morphology and content of the filling medium, a three-dimensional digital core model with the medium filling the pore center is constructed using image morphology methods, utilizing a three-dimensional digital core model that identifies the skeleton and pores. The specific steps are as follows:

[0296] S201. Obtain the filling medium type and filling medium content TBtarget identified in step S1. 中心 And the maximum diameter of the core pores, dmax; Figure 2 The filling medium was identified as pore center filling. Figure 3 The display identifies the filling medium, and the filling type is pore edge filling.

[0297] S202, Three-dimensional digital core model based on identification of skeleton and pores (e.g.) Figure 4 As the initial three-dimensional digital core set filling the pore centers of the medium, the first skeleton set and the first pore set are divided according to the phase state function of each element; Figure 5 Displays the distribution of the first pore set in a three-dimensional digital core model that identifies the skeleton and pores.

[0298] S203. In the first pore set representing the pore space, the image morphology-based method uses spherical structural elements of different diameters to perform opening operations in the pore set. The opening operation results of the spherical structural elements of different diameters are obtained through the opening operation, and the number of elements in the opening operation results set of each diameter spherical structural element and the proportion of the number of elements in the first pore set are obtained.

[0299] S204. Based on the set of opening operation results of spherical structural elements with different diameters, obtain the difference set of the opening operation results set of two adjacent diameter spherical structural elements, and perform Euclidean distance calculation in the difference set set. According to the Euclidean distance, divide the difference set set into different Euclidean distance sets, and obtain the number of pore elements in the corresponding Euclidean distance set and the proportion of the corresponding number of pore elements in the first pore set.

[0300] S205. Based on the first opening operation result set, the number of pore elements in the opening operation result set of sphere structural elements of each diameter, the proportion of the number of pore elements in the opening operation result set of sphere structural elements of each diameter in the first pore set, the difference set of the opening operation result sets of sphere structural elements of two adjacent diameters, the Euclidean distance set in the difference set of the opening operation result sets of sphere structural elements of adjacent diameters, the proportion of the number of pore elements in the Euclidean distance set in the second pore set, and the content of the filling medium, construct a three-dimensional digital core model of the filling of the pore center of the medium. Figure 6 The three-dimensional digital core model shows the filling of the pore center of the medium. Figure 7This displays the distribution of the pore set and the filling medium set in the three-dimensional digital core model after the pore centers of the medium have been filled. Figure 8 This displays the distribution of the filling medium aggregate in the three-dimensional digital core model after the pore centers of the medium have been filled. Figure 9 The distribution of pore sets in a three-dimensional digital core model after the pore centers of the medium are filled.

[0301] In step S202,

[0302] Any element in the 3D digital core model set is represented by px. 中心 The phase function `phase` is used to indicate whether a point is part of the skeleton or a pore, where `phase(px)` represents the phase state. 中心 ) = 0 means px 中心 The element is a skeleton, phase(px) 中心 ) = 1 represents px 中心 The element is a pore, phase(px) 中心 ) = 3 indicates that the element is filling the medium; the three-dimensional digital core model for identifying the skeleton and pores is represented by the set CORE1, which is the first skeleton set MATRIX. 中心 and the first pore set PORE 中心 The union of, i.e., CORE1 = MATRIX 中心 ∪PORE 中心 Furthermore, the core skeleton assembly MATRIX 中心 and core pore set PORE 中心 Represented as:

[0303] MATRIX 中心 ={px 中心 |phase(px 中心 )=0};

[0304] PORE 中心 ={px 中心 |phase(px 中心 )=1};

[0305] Calculate the first pore set PORE 中心 The number of elements in nump 中心 .

[0306] In step S203, the main task is to realize the opening operation set of spherical structural elements of different diameters in the first pore set of the initial three-dimensional digital core model, the number of elements in the set, and the proportion of the number of elements in the set to the total number of elements in the pore space set.

[0307] ① Determine the maximum diameter of the core pores as dmax;

[0308] ② The set of initial three-dimensional digital core models (CORE) filled in the pore centers of the medium. 中心 The first pore set PORE is selected to represent the first pore set. 中心 and the first pore set PORE 中心 The number of elements in nump 中心 .

[0309] ③ Select the diameter of the sphere structural element in the image morphology opening operation as i, and the range of i is from 1 to dmax, i.e. i = 1, 2, ..., dmax;

[0310] ④ For the first pore PORE 中心 Each element in px 中心 A morphologically based method is used to perform an opening operation OPEN(px,i) on spherical structural elements with diameter i. OPEN(px,i) = 0 indicates that the element does not satisfy the opening operation requirement for spherical structural elements with diameter i in the pore space set, while OPEN(px,i) = 1 indicates that the element satisfies the opening operation requirement for spherical structural elements with diameter i in the pore space set. This yields the set MORPH of opening operation results for spherical structural elements with diameter i. i This forms the first set of results from the opening operation.

[0311] MORPH i ={px|OPEN(px,i)=1},(i=1,2,…,dmax);

[0312] ⑤ Based on the first open operation result set MORPH i Calculate the number of elements in a set using the Onum function. i ,(i=1,2,…,dmax);

[0313] ⑥ Based on the first open operation result set MORPH i MORPH i The number of elements in the set (Onum) i ,(i=1,2,…,dmax), calculate MORPH i The number of elements in the set is in the first pore set PORE 中心 The proportion of OSB i .

[0314] OSBi = Onum i / nump 中心 , (i = 1, 2, ..., dmax);

[0315] In step S204,

[0316] ① Select the first set of results from the opening operation, MORPH iThe number of elements in the set, Onum i and MORPH i The number of elements in the set is in the first pore set PORE 中心 The proportion of OSB i , (i = 1, 2, ..., dmax-1).

[0317] ② Select the set of opening operation results of the sphere structuring element with diameter i+1 in the first opening operation set, MORPH i+1 The number of elements in the set, Onum i+1 and MORPH i+1中心 The number of elements in the set is in the first pore set PORE 中心 The proportion of OSB i+1 , (i = 1, 2, ..., dmax-1).

[0318] ③ Based on the set MORPH in the result of the first opening operation i MORPH i+1 DMORPH is obtained through difference operations. i (i = 1, 2, ..., dmax-1).

[0319] DMORPH i =MORPH i -MORPH i+1 ;

[0320] ④ Calculate the difference set DMORPH i any element in px 中心 To the set MORPH i+1 The Euclidean distance is j = dist(px) 中心 ), (i = 1, 2, ..., dmax-1);

[0321] ⑤ In the difference set DMORPH i In (i = 1, 2, ..., dmax-1), based on the Euclidean distance j = dist(px) 中心 (j = 1, 2, ..., n, where n is the largest Euclidean distance in the set) The difference set DMORPH i The set is partitioned into a series of Euclidean distance sets, DIST. j .

[0322] DIST j ={px 中心 |dist(px 中心 ) = j, and px 中心 ∈DMORPH i-1},(j=1,2,...,n);

[0323]

[0324] ⑥ Calculate the difference set DMORPH i Central Euclidean distance set DIST j The number of elements in Dnum j中心 ;

[0325] ⑦ Calculate the difference set DMORPH i Central Euclidean distance set DIST j The number of elements in Dnum j First pore set PORE 中心 The proportion of DSB j ;

[0326] DSB j =Dnum j / nump 中心 j = 1, 2, ..., n;

[0327] In step S205,

[0328] Based on the set MORPH i and the percentage of elements in this set (OSB) i and the difference set of adjacent opening operations, DMORPH i-1 Central Euclidean distance set DIST i The number and proportion of pores in DSB i and filling medium content TBtarget 中心 This allows the medium to gradually fill the pore space from the pore center (i = 2, 3, ..., dmax);

[0329] S2051, if OSB dmax >TBtarget 中心 Then MORPH dmax All pores in the set are transformed into filling media, establishing a media set BIT.

[0330]

[0331] The number of elements in the filling medium set is Onum dmax .

[0332] In the formula, TBtarget 中心 Indicates the content of the first filling medium, MORPH dmax Onum represents the set of elements in the first set of pores that satisfy the result of the opening operation of a sphere with diameter dmax. dmax MORPH dmax Number of elements in the set; OSB dmax Represents Onum dmaxFirst pore set PORE 中心 The proportion of the number of elements;

[0333] Complete the center filling of the medium pores and obtain the center-filled medium set (BIT). 中心 The next step is to proceed to step S2053;

[0334] S2052, If OSB i <TBtarget 中心 <OSB i-1 i = 2, 3, 4, ..., d max hour;

[0335] ① First, open the set of structural elements of the sphere with diameter i, MORPH. i All elements in the set are transformed into the filling medium set BIT1. 中心 ,

[0336] That is: BIT1 中心 =MORPH i MORPH i ={px 中心 |OPEN(px 中心 ,i)=1};

[0337] In the formula, BIT1 中心 Indicates center-filled medium set 1, MORPH i OSB represents the set of elements in the first set of pores that satisfy the result of the opening operation of a sphere with diameter i. i MORPH i The elements in the first pore set PORE 中心 Percentage of elements; px 中心 This represents an element in the initial three-dimensional digital core model filled with pores in the medium; OPEN(px 中心 ,i)=1 indicates that the center of px satisfies the opening operation result of the sphere structure element with diameter i.

[0338] ②If DSB1>TBtarget 中心 -OSB i (i = 2, 3, 4, ..., d) max Then, randomly select (nump) from the DIST1 set. 中心 ×TBtarget 中心 -OSB i ) pores are transformed into filling medium, establishing a central filling medium set 2 (BIT2) 中心 ),

[0339] BIT2 中心 ={px 中心Randomly select px 中心 , px 中心 ∈DIST1};

[0340] DMORPH i-1 =MORPH i-1 -MORPH i ;

[0341] DIST1 = {px 中心 |dist(px 中心 ) = 1, and px 中心 ∈DMORPH i-1};

[0342] Complete the filling of the medium pore center and obtain the collection of medium pore center filling BIT. The next step is to proceed to step S2053.

[0343] OSB i <TBtarget 中心 <OSB i-1 And DSB1>TBtarget 中心 -OSB i (i = 2, 3, ..., d) max When ), the filling medium set is:

[0344] BIT 中心 =BIT1 中心 ∪BIT2 中心 ;

[0345] BIT1 中心 =MORPH i MORPH i ={px 中心 |OPEN(px 中心 ,i)=1};

[0346] BIT2 中心 ={px 中心 Randomly select px 中心 , px 中心 ∈DIST1};

[0347] DMORPH i-1 =MORPH i-1 -MORPH i ;

[0348] DIST1 = {px 中心 |dist(px 中心 ) = 1, and px 中心 ∈DMORPH i-1};

[0349] In the formula, BIT 中心 BIT1 represents the centrally filled medium set. 中心 This indicates the filling medium set 1, BIT2. 中心 Indicates the central filling medium set 2; px 中心 This represents an element in the initial three-dimensional digital core model filled with pores in the medium; OPEN represents the opening operator, and i and i-1 represent the diameters of the spherical structural elements, respectively; OPEN(px 中心 ,i)=1 indicates that the center of px satisfies the opening operation result of a sphere structuring element with diameter i, where i=2,3,…,d max MORPH i OSB represents the set of elements in the first set of pores that satisfy the result of the opening operation of a sphere with diameter i; i MORPH i The elements in the first pore set PORE 中心 The proportion of elements in the total number; MORPH i-1 OSB represents the set of elements resulting from the opening operation of spheres with diameter i-1 in the first pore set; i-1 MORPH i-1 The elements in the first pore set PORE 中心 The percentage of the total number of elements.

[0350] MORPH i DMORH represents the set of elements in the first set of pores that satisfy the result of the opening operation of a sphere structure element with diameter i; i-1 DIST1 represents the difference set of the results of the opening operation between sphere structural elements of diameter i-1 and i in the first pore set based on image morphology. i-1 The Euclidean distance in the set is 1, i.e., dist(px) 中心 The set consisting of elements where ) = 1. DSB1 represents the set of elements in DIST1 in the first pore set PORE. 中心 The percentage of elements; dist represents the Euclidean distance operator.

[0351] ③If OSB i <TBtarget 中心 <OSB i-1 (i = 2, 3, ..., d) max And OSB i +∑ j＝1 i DSB j <TBtarget 中心 <OSB i +∑ j＝1 i+1 DSBj >TBtarget 中心 -OSB i ;

[0352] First, open the set of structuring elements of the sphere with diameter i, MORPH. i All elements in the set are transformed into the filling medium set BIT1′. 中心 ,

[0353] That is: BIT1' 中心 =MORPH i MORPH i ={px 中心 |OPEN(px 中心 ,i)=1};

[0354] In the formula, BIT1′ 中心 Indicates center-filled medium set 1'; MORPH i OSB represents the set of elements in the first set of pores that satisfy the result of the opening operation of a sphere with diameter i; i MORPH i The elements in the first pore set PORE 中心 Percentage of elements; px 中心 This represents an element in the initial three-dimensional digital core model filled with pores in the medium; OPEN(px 中心 ,i)=1 indicates that the center of px satisfies the opening operation result of the sphere structure element with diameter i.

[0355] Secondly, the first open operation set MORPH i-1 MORPH i DMORPH difference set i-1 All elements in the set whose distance is less than or equal to i fill the medium set BIT2′. 中心 ,Right now:

[0356]

[0357] DIST j ={px 中心 |dist(px 中心 ) = j, and px 中心 ∈DMORPH i-1};

[0358] DMORPH i-1 =MORPH i-1 -MORPH i ;

[0359] In the formula, BIT2′ 中心Indicates center-filled medium set 2', MORPH i OSB represents the set of elements in the first set of pores that satisfy the result of the opening operation of a sphere with diameter i. i MORPH i The elements in the first pore set PORE 中心 The proportion of elements in the total number; MORPH i-1 OSB represents the set of elements in the first set of pores that satisfy the result of the opening operation of a sphere with diameter i-1. i-1 MORPH i-1 The elements in the first pore set PORE 中心 Percentage of elements; px 中心 This represents an element in the initial three-dimensional digital core model filled with pores in the medium. OPEN represents the opening operator, and i and i-1 represent the diameters of the spherical structural elements, respectively. OPEN(px 中心 ,i)=1 indicates that the center of px satisfies the opening operation result of a sphere structure element with diameter i. DMORH i-1 This represents the difference set of the results of the opening operation between sphere structural elements of diameter i-1 and i in the first pore set, based on image morphology. j Indicates DMORPH i-1 The set of elements whose Euclidean distance is j. ∪ denotes the union of sets. j＝1 i DIST j Represents sets DIST1, DIST2, ..., DIST i The union of DSB. j DIST j The elements in the first pore set PORE 中心 The proportion of the number of elements;

[0360] Finally, in the first open set MORPH i-1 MORPH i DMORPH difference set i-1 The set whose Euclidean distance is equal to i+1 is randomly selected (TBtarget) 中心 -OSB i -∑ j＝1 i DSB j )×nump 中心 Each element is transformed into a filling medium set BIT3′ 中心 ,Right now:

[0361] BIT3' 中心 ={px 中心 Randomly select px 中心, px 中心 ∈DIST j+1};

[0362] DIST j+1 ={px 中心 |dist(px 中心 ) = j + 1, and px 中心 ∈DMORPH i-1};

[0363] DMORPH i-1 =MORPH i-1 -MORPH i ;

[0364] In the formula, BIT3′ 中心 Indicates center-filled medium set 3; MORPH i MORPH represents the set of elements in the first set of pores that satisfy the result of the opening operation of a sphere with diameter i; i-1 DMORH represents the set of elements in the first set of pores that satisfy the opening operation of spherical structural elements with diameter i-1; i-1 DIST represents the difference set of the results of the opening operation between sphere structural elements with diameters i-1 and i in the first set of pores, based on image morphology; j+1 Indicates DMORPH i-1 The set of elements whose Euclidean distance is j+1; DSB j+1 DIST j+1 The elements in the first pore set PORE 中心 nump of element count 中心 The proportion of DIST; j Indicates DMORPH i-1 The set of elements whose Euclidean distance is j; OSB i MORPH i The elements in the first pore set PORE 中心 Percentage of elements; OSB i-1 MORPH i-1 The elements in the first pore set PORE 中心 Percentage of elements; px 中心 This represents an element in the initial three-dimensional digital core model filled with pore centers of the medium, nump 中心 PORE represents the first set of pores. 中心 Number of elements; DSB j DIST j The elements in the first pore set PORE 中心 nump of element count 中心The proportion of.

[0365] That is, when OSB i <TBtarget 中心 <OSB i-1 And OSB i +∑ j＝1 i DSB j <TBtarget 中心 <OSB i +∑ j＝1 i+1 At that time, the filling medium set is as follows:

[0366] BIT 中心 =BIT1' 中心 ∪BIT2' 中心 ∪BIT3' 中心 ;

[0367] BIT1' 中心 =MORPH i MORPH i ={px 中心 |OPEN(px 中心 ,i)=1};

[0368]

[0369] BIT3' 中心 ={px 中心 Randomly select px 中心 , px 中心 ∈DIST j+1};

[0370] DMORPH i-1 =MORPH i-1 -MORPH i ;

[0371] DIST j ={px 中心 |dist(px 中心 ) = j, and px 中心 ∈DMORPH i-1};

[0372] DIST j+1 ={px 中心 |dist(px 中心 ) = j + 1, and px 中心 ∈DMORPH i-1}(i=2,3,…,dmax);

[0373] In the formula, BIT1′ 中心Indicates the central filling medium set 1'; BIT2' 中心 Indicates the central filling medium set 2'; BIT3' 中心 Indicates center-filled medium set 3'; MORPH i MORPH i-1 OSB represents the set of elements in the first set of pores that satisfy the opening operation of spheres with diameters of i and i-1; i MORPH i The elements in the first pore set PORE 中心 Percentage of elements; OSB i-1 MORPH i-1 The elements in the first pore set PORE 中心 The proportion of elements in the total number; DMORH i-1 MORPH represents the set of results of opening operations on spherical structural elements with diameters i-1 and i in the first set of pores, based on image morphology. i MORPH i-1 The difference set (i = 2, 3, ..., dmax); DIST j Indicates DMORPH i-1 The set of elements whose Euclidean distance is j; DSB j DIST j The elements in the first pore set PORE 中心 The proportion of elements in the total number; DIST j+1 Indicates DMORPH i-1 The set consisting of elements whose Euclidean distance is j+1. DSB j+1 DIST j+1 The elements in the first pore set PORE 中心 nump of element count 中心 The proportion in; dist represents the Euclidean distance operator; dist(px 中心 ) indicates the calculation element px 中心 Euclidean distance; px 中心 OPEN(px) represents an element in the initial three-dimensional digital core model filled with pore centers of the medium. 中心 , i) = 1 indicates that the center of px satisfies the opening operation result of the sphere structuring element with diameter i; ∪ represents the union, ∪ j =1 i DIST j Represents sets DIST1, DIST2, ..., DIST i The union of .

[0374] OPEN represents the opening operator, where i and i-1 represent the diameters of the sphere structuring element, respectively; OPEN(px中心 , i) = 1 indicates that the center of px satisfies the opening operation result of a sphere structure element with diameter i; DSB j DIST j The elements in the first pore set PORE 中心 The percentage of elements in the total number; DSB j+1 DIST j+1 The elements in the first pore set PORE 中心 The percentage of elements; dist represents the Euclidean distance operator.

[0375] This completes the process of filling the center of the medium pores.

[0376] S2053, Establishing a central filling medium set (BIT) 中心 Subsequently, the three-dimensional digital core assembly CORE 中心 Represented as a skeleton set MATRIX 中心 Pore ​​set PORE 中心 BIT (Body Filler) 中心 The union of, i.e.:

[0377] CORE 中心 =MATRIX 中心 ∪PORE 中心 ∪BIT 中心 ;

[0378] In the formula, MATRIX 中心 PORE represents the first skeleton set. 中心 CORE represents the first set of pores. 中心 A three-dimensional digital core set representing a pore space model where the medium gradually fills the pore center, where U represents the union, and the first skeleton set is MATRIX. 中心 Represented as:

[0379] MATRIX 中心 ={px 中心 |phase(px 中心 )=0};

[0380] First Pore Set PORE 中心 Represented as:

[0381] PORE 中心 ={px 中心 |phase(px 中心 )=1};

[0382] Central Medium Filling Set (BIT) 中心 Represented as:

[0383] BIT 中心 ={px中心 |phase(px 中心 )=3};

[0384] In the formula, px 中心 This represents an element in the initial three-dimensional digital core model filled with pore centers of the medium, where phase represents the phase state function, phase(px 中心 ) = 0 means px 中心 The element is a skeleton, phase(px) 中心 ) = 1 represents px 中心 The element is a pore. phase(px) 中心 ) = 3 represents px 中心 The element is the filling medium.

[0385] Example 3

[0386] As another preferred embodiment of the present invention, this embodiment further supplements and elaborates on the technical solution of the present invention based on the above-described Embodiment 1 or Embodiment 2. Based on the distribution morphology and content of the filling medium, a three-dimensional digital core model identifying the skeleton and pores is used, and an image Euclidean distance map method is employed to construct a three-dimensional digital core model of the medium filling at the pore edges. The specific steps are as follows:

[0387] S301. Obtain the filling medium type and filling medium content TBtarget obtained in step S1. 边缘 ;

[0388] S302. Based on the phase state function of each point in the three-dimensional digital core model, the second skeleton and the second pore set of the second initial three-dimensional core are divided, and the number of elements in the second pore set is obtained.

[0389] S303. In the second set of pores in the identified pore space, calculate the Euclidean distance from each pore point to the nearest skeleton point, obtain the Euclidean distance set of the pore space, and obtain the number of pore elements in each Euclidean distance set and the proportion of the number of pore elements in the second pore set.

[0390] S304. Based on the set of Euclidean distances in the pore space of the three-dimensional digital core, the number and proportion of pore elements in each set of Euclidean distances, and the content of the filling medium, a three-dimensional digital core model of the filling medium edge is constructed.

[0391] Further optimized

[0392] In step S302, the three-dimensional digital core model identifying the skeleton and pores serves as the initial three-dimensional digital core model for filling the pore edges of the medium, denoted by the second initial core set CORE2, where:

[0393] CORE2 = MATRIX边缘 ∪PORE 边缘 ;

[0394] In the formula, MATRIX 边缘 PORE represents the second skeleton set. 边缘 Let U denote the second set of pores, and U denote the union of the second skeleton set, MATRIX. 边缘 Represented as:

[0395] MARTIX 边缘 ={px 边缘 |phase(px 边缘 )=0};

[0396] Second pore set PORE 边缘 Represented as:

[0397] PORE 边缘 ={px 边缘 |phase(px 边缘 )=1},

[0398] Second pore set PORE 边缘 The number of elements in nump 边缘 .

[0399] In the formula, px 边缘 This represents the elements in the initial three-dimensional digital core set filled with pore edges of the medium, where phase represents the phase state function, phase(px 边缘 ) = 0 means px 边缘 The element is a skeleton, phase(px) 边缘 ) = 1 represents px 边缘 The element is porosity.

[0400] In step S303, ① the initial three-dimensional digital core model CORE2, which identifies the skeleton and pores, is used as the pore set PORE for filling the pore edges of the medium. 边缘 Based on the calculated Euclidean distance dist, the pore set PORE is... 边缘 Let PDIST be the set of second Euclidean distances with Euclidean distance i. i边缘 The union of, i.e.:

[0401]

[0402] The minimum Euclidean distance is 1, the maximum Euclidean distance is n, and the second pore set PORE 边缘 Mid-range skeleton ensemble MATRIX 边缘 The set of elements with a Euclidean distance of i is PDIST. i ;

[0403] ② The pore set PORE of the initial three-dimensional digital core model CORE2, which is filled at the pore edges of the medium. 边缘 PDIST, a second Euclidean distance set based on Euclidean distance partitioning i Calculate the n Euclidean distance sets PDIST in (i = 1, 2, ..., n). i The number of elements in each set is divided into PDnum. i If i = 1, 2, ..., n, then we have

[0404] ③ The pore set PORE of the initial three-dimensional digital core model CORE2, which is filled with pore edges of the medium. 边缘 The number of elements nump 边缘 With n sets of second Euclidean distances PDIST i Number of elements PDnum i Calculate the corresponding Euclidean distance set PDIST i Number of elements in PDnum i In the pore set PORE 边缘 The proportion of PDSB in i ,

[0405]

[0406] In step S304, the pore set PORE of the initial three-dimensional digital core model CORE2, which is filled with pore edges of the medium, is used as the basis. 边缘 The second Euclidean distance set PDIST i (i = 1, 2, ..., n) and the percentage of elements in this set PDSB i and filling medium content TBtarget 边缘 The medium is gradually filled from the edge with a Euclidean distance of 1 towards the pore space with an increasing Euclidean distance. The specific steps are as follows:

[0407] S3041, If ​​PDSB 1边缘 >TBtarget 边缘 Then, nump is randomly selected from the PDIST1 set. 边缘 ×TBtarget 边缘 The pores are transformed into filling media, and a three-dimensional digital core model of the filling media edge is established, in which the edge filling media set BIT is formed. 边缘 The next step is step S3043.

[0408] BIT 边缘 ={px 边缘 Randomly select px 边缘 , px 边缘 ∈PDIST1};

[0409] PDIST1 = {px 边缘 |dist(px 边缘 ) = 1, and px 边缘 ∈PORE 边缘 The number of elements in the filling medium set is nump 边缘 ×TBtarget 边缘 ;

[0410] In the formula, BIT 边缘 TBtarget represents the set of edge-filled media. 边缘 Indicates the content of the edge-filling medium, nump 边缘 PORE represents the second set of pores filled at the edge. 边缘 The number of elements, px 边缘 This represents an element in the initial three-dimensional digital core model filled with pore edges of the medium. `dist` represents the Euclidean distance operator, and `PDIST1` represents the element in the second pore set `PORE`. 边缘 The set consisting of elements with a distance of 1 in the Euclidean distance set, PDSB1 represents the proportion of the number of elements with a distance of 1 in the Euclidean distance set in the second pore set.

[0411] S3042, If the pore set PORE of the second initial three-dimensional digital core model CORE2 is... 边缘 The set of all second Euclidean distances PDIST with mean Euclidean distance dist≤i j The percentage of elements in the sequence (j = 1, 2, ..., i) is less than TBtarget. 边缘 The set of all second Euclidean distances PDIST for dist≤i+1 j The sum of the percentages of elements in the sequence (j = 1, 2, ..., i, i+1) is greater than TBtarget. 边缘 hour,

[0412] Right now hour,

[0413] Then ① the pore set PORE of the initial three-dimensional digital core model CORE2, which fills the pore edges of the medium, is... 边缘 The medium Euclidean distance is the set of all second Euclidean distances less than or equal to i, PDIST. j The elements (j = 1, 2, ..., i) are used as filling media to establish edge filling medium set 1 (BIT1). 边缘 );

[0414] Right now

[0415] In the formula, BIT1 边缘 Indicates edge-filling medium set 1, TBtarget 边缘 Indicates the content of edge-filling medium, px边缘 Represents the second pore set PORE 边缘 In the element, dist represents the Euclidean distance operator, dist(px 边缘 ) represents the pore set PORE 边缘 medium element px 边缘 Euclidean distance to the skeleton set; PDIST j Indicated in the second pore set PORE 边缘 The second Euclidean distance set, PDSB, is composed of elements with a medium Euclidean distance of j. j This indicates the number of elements with Euclidean distance j in the second Euclidean distance set in the second pore set PORE. 边缘 The proportion in. ∪ j＝1 i PDIST j Let PDIST1, PDIST2, ..., PDIST be the set. i The union of ∪; ∪ means union.

[0416] ② In the set PDIST of dist = i+1 i+1 Random selection Each porous element is transformed into a filling medium, namely BIT2. 边缘 ={px 边缘 Randomly select px 边缘 , px 边缘 ∈PDIST i+1};

[0417] ③ The filling medium content is TBtarget 边缘 In the three-dimensional digital core model of pore edge filling, the edge filling medium assembly BIT 边缘 For: BIT 边缘 =BIT1 边缘 ∪BIT2 边缘 This completes the filling of the pore center of the medium and obtains the edge-filled medium assembly (BIT). 边缘 ;

[0418] when hour,

[0419] The edge filling medium set is as follows:

[0420] BIT 边缘 =BIT1 边缘 ∪BIT2 边缘 ;

[0421]

[0422] BIT2 边缘 ={px 边缘 Randomly select px边缘 , px 边缘 ∈PDIST i+1};

[0423] PDIST i+1 ={px 边缘 |dist(px 边缘 ) = i + 1, and px 边缘 ∈PORE 边缘};

[0424] In the formula, BIT 边缘 Indicates the set of edge-filling media; BIT1 边缘 BIT2 represents the edge-filling medium set 1; 边缘 Indicates edge-filling medium set 2; TBtarget 边缘 Indicates the content of edge-filling medium, px 边缘 This represents the second pore set PORE in the second initial three-dimensional digital core model. 边缘 The elements in the text; OPEN represents the opening operator, dist represents the Euclidean distance operator, dist(px) 边缘 ) represents the set of pores, PORE 边缘 medium element px 边缘 Euclidean distance to the skeleton set; PDIST j Let j represent the second Euclidean distance set of elements with Euclidean distance j.

[0425] Among them, BIT1 边缘 The second pore set PORE represents the second initial three-dimensional digital core model CORE2. 边缘 In the set, the union of the sets of Euclidean distances less than or equal to i; BIT2 边缘 The second pore set PORE represents the second initial three-dimensional digital core model CORE2. 边缘 In the equation, the second Euclidean distance set PDIST is equal to i+1. i+1 The set is formed by randomly selecting elements from the set.

[0426] S3043. After establishing the medium-filled assemblies (BITs), the three-dimensional digital core assemblies (COREs) are created. 边缘 It can be expressed as:

[0427] CORE 边缘 =MATRIX 边缘 ∪PORE 边缘 ∪BIT 边缘 ;

[0428] BIT 边缘 ={px 边缘 |phase(px 边缘)=3};where,

[0429] CORE2 represents the second initial three-dimensional digital core set filled with edge media. 边缘 MATRIX represents a three-dimensional digital core assembly after the media edge is filled. 边缘 Represents the second skeleton set, PORE 边缘 Represents the second set of pores, BIT 边缘 Represents the set of edge-filling media, phase(px) 中心 ) = 3 indicates that the element fills the medium, px 边缘 The second set of pores, PORE, represents the second initial three-dimensional digital model filled with the pore centers of the medium. 边缘 The elements in.

[0430] Example 4

[0431] As another preferred embodiment of the present invention, this embodiment is a further detailed supplement and explanation of the technical solution of the present invention based on the above-described embodiments 1, 2, or 3. In this embodiment, reference is made to the appendix to the specification. Figure 4 To the attached Figure 15 Based on a three-dimensional digital core model that identifies the skeleton and pores, different combinations of image morphology and image Euclidean distance mapping methods are used to construct three-dimensional digital core models with different filling methods of the medium. Specifically, based on steps S2 and S3, by adjusting the order of steps S2 and S3, a three-dimensional digital core model is constructed where the medium first fills the pore edges and then the pore center, or a three-dimensional digital core model is constructed where the medium first fills the pore center and then the pore edges. Specifically, the image morphology method and image Euclidean distance mapping method described in Examples 2 and 3 above are used for edge filling and center filling.

[0432] As one implementation method of this embodiment, the construction process of a three-dimensional data core model in which the medium is first filled at the pore edges and then filled at the pore center is as follows:

[0433] S4A1. Based on the identification of the skeleton and pores, the initial three-dimensional digital core set CORE is represented as the union of the skeleton set MATRIX and the pore set PORE, that is, CORE = MATRIX∪PORE. Figure 4 Displays the skeleton and pore sets of a 3D digital core model that identifies the skeleton and pores. Figure 5 A three-dimensional display of the pore set in a three-dimensional digital core model that identifies the skeleton and pores.

[0434] S4A2. Based on the content of the central filling medium, a three-dimensional data core model of the medium filling the center of the pores is constructed in the current pore set PORE using image morphology method, i.e., CORE=MATRIX∪PORE∪BIT1; Figure 6 The three-dimensional digital core model shows the filling of the pore center of the medium. Figure 7 This displays the distribution of the pore set and the filling medium set in a three-dimensional digital core model after the pore centers of the medium have been filled. Figure 8 This displays the distribution of the filling medium assembly in a three-dimensional digital core model after the pore centers of the medium have been filled. Figure 9 Distribution of pore sets in a three-dimensional digital core model after filling the pore centers of the medium.

[0435] S4A3. Based on the three-dimensional digital core model of the medium filling the pore center, i.e., CORE=MATRIX∪PORE∪BIT1, according to the content of the edge filling medium, a three-dimensional digital core model of the medium filling the pore center and then filling the pore edge is constructed in the current pore set PORE using the image Euclidean distance method, i.e., CORE=MATRIX∪PORE∪BIT∪BIT2. Figure 10 Display of the pore set and skeleton set of the three-dimensional digital core model after the pore centers of the medium have been filled. Figure 11 This displays the set of pores in the three-dimensional digital core model after the pore center of the medium has been filled, followed by pore edge filling. Figure 12 The distribution of pore sets and framework sets in a three-dimensional digital core model, where the medium first fills the pore centers and then the pore edges. Figure 13 The distribution of pore sets and filling medium sets in a three-dimensional digital core model shows that the medium first fills the pore centers and then the pore edges. Figure 14 This shows the distribution of the pore set in a three-dimensional digital core model, where the pore center is filled first, followed by the pore edges. Figure 15 The distribution of the filling medium set in a three-dimensional digital core model, in which the medium first fills the pore center and then the pore edge.

[0436] As another implementation of this embodiment, the construction process of the three-dimensional digital core model in which the medium is first filled at the pore edges and then filled at the pore center is as follows:

[0437] S4B1. Based on the three-dimensional digital core model that identifies the skeleton and pores, the initial three-dimensional digital core set CORE is represented as the union of the skeleton set MATRIX and the pore set PORE, i.e., CORE = MATRIX ∪ PORE.

[0438] S4B2. Based on the content of the edge-filling medium, a three-dimensional digital core model of the medium filling the pore edge is constructed in the current pore set PORE using the image Euclidean distance method, i.e., CORE=MATRIX∪PORE∪BIT2;

[0439] S4B3. Based on the three-dimensional digital core model of the medium filling the pore edges, i.e., CORE=MATRIX∪PORE∪BIT2, according to the amount of medium filling the center, a three-dimensional digital core model of the medium filling the pore center after filling the pore edges is constructed in the current pore set PORE using an image morphology-based method, i.e., CORE=MATRIX∪PORE∪BIT1∪BIT2.

[0440] Example 5

[0441] As another preferred embodiment of the present invention, this embodiment is a further detailed supplement and explanation of the technical solution of the present invention based on the above-described embodiments 1, 2, 3, or 4. In this embodiment, reference is made to the accompanying description. Figure 16 , attached Figure 17 , attached Figure 18 , attached Figure 19 and attached Figure 20 As shown, Figure 16 , Figure 17 , Figure 18 The influence of medium filling on pore structure is illustrated using a pore network model. Figure 19 and Figure 20 This shows the changes in pore structure parameters and resistivity.

[0442] Step S5 further includes the following sub-steps:

[0443] S501, Assignment of resistivity properties of each phase in three-dimensional digital core filled with medium;

[0444] S502. Finite element method simulation calculation of core resistivity;

[0445] S503. Extract the pore network model using the maximum sphere method to obtain pore structure parameters; Figure 16 A pore network model that displays a 3D digital core pore set that identifies pores and the core structure. Figure 17 This displays a pore network model of the 3D digital core pore set after filling the pore centers of the medium in the pore set using a pore and skeleton identification model. Figure 18 This displays a pore network model of the 3D digital core pore set based on the identification of pores and skeleton, after filling the center of the medium pores in the pore set and then filling the edge of the medium pores.

[0446] S504. Study the influence of dielectric filling on pore structure parameters and resistivity. Figure 19 The display shows the changes in pore structure parameters after the medium is filled. It can be seen that the filling of the medium center blocks the large pores, resulting in a decrease in pore radius, while the filling of the medium edge blocks the corners of the pores, resulting in a decrease in small pores and an increase in pore radius. Figure 20 The effect of dielectric filling on resistivity is shown. Filling the pore center reduces porosity by 0.05 and increases resistivity by 100 ohm-meters, while filling the pore edge reduces porosity by 0.015 and increases resistivity by 600 ohm-meters. This indicates that dielectric filling causes a decrease in porosity and an increase in resistivity, and the effect of filling the pore edge on resistivity is more significant.

[0447] Finally, it should be noted that the above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. Although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art can still modify the technical solutions described in the foregoing embodiments or make equivalent substitutions for some of the technical features. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the protection scope of the present invention.

Claims

1. A method for constructing a three-dimensional digital core model filled with a medium, characterized in that, The method includes, Based on two-dimensional core images, the type, distribution morphology, and content of the core filling medium were obtained; Based on the distribution morphology and content of the filling medium, a three-dimensional digital core model with the medium filling the center of the pores and a three-dimensional digital core model with the medium filling the edge of the pores are constructed using a three-dimensional digital core model that identifies the skeleton and pores. Based on three-dimensional digital core models with medium filling the pore center and three-dimensional digital core models with medium filling the pore edge, a three-dimensional digital core model with medium filling is constructed. The three-dimensional digital core model, in which the construction medium fills the pore centers, includes: A three-dimensional digital core model that identifies the skeleton and pores is used as the initial three-dimensional digital core model for filling the pore center of the medium. Based on the initial three-dimensional digital core model filled with the pore center of the medium, the first skeleton set and the first pore set are obtained, and the number of elements in the first pore set is obtained. Based on image morphology, an opening operation is performed on spherical structural elements of different diameters in the first pore set to obtain the first opening operation result set, and the number of first pore elements in the first opening operation result set and the proportion of the number of first pore elements in the first pore set are obtained. Based on the first set of opening operation results, obtain the difference set of the opening operation results sets of two adjacent diameter sphere structure elements; Perform Euclidean distance calculation on the difference set to obtain the number of pore elements in the corresponding Euclidean distance set and the proportion of the corresponding number of pore elements in the first pore set of the initial three-dimensional digital core model filled with the pore center of the medium. Based on the opening operation result sets of spherical structural elements of different diameters, the number of pore elements in each diameter's opening operation result set, the proportion of the number of pore elements in each diameter's opening operation result set in the pore set of the three-dimensional digital core model filled with medium pore centers, the difference set of the opening operation result sets of two adjacent diameter spherical structural elements, the Euclidean distance set in the difference set of the opening operation result sets of adjacent diameter spherical structural elements, the proportion of the number of pore elements in each Euclidean distance in the difference set in the pore set of the initial three-dimensional digital core model filled with medium pore centers, and the filling medium content, a three-dimensional digital core model filled with medium pore centers is constructed; where, if OSB dmax >TBtarget 中心 Then MORPH dmax All the pores in the medium are transformed into a filling medium; If OSB i <TBtarget 中心 <OSB i-1 And DSB1>TBtarget 中心 -OSB i When i = 2, 3, 4, ..., dmax, then MORPH i All the pores in the medium are transformed into a filling medium; If OSB i <TBtarget 中心 <OSB i-1 and , When i = 2, 3, 4, ..., dmax, first open the set of sphere structuring elements with diameter i, MORPH. i All elements in the set are transformed into the filling medium set BIT1′. 中心 Secondly, the first open operation set MORPH... i-1 MORPH i DMORPH difference set i-1 All elements in the set whose distance is less than or equal to i fill the medium set BIT2′. 中心 Finally, in the first open set MORPH i-1 MORPH i DMORPH difference set i-1 The set whose medium-Euclidean distance is equal to i+1 is randomly selected. Each element is transformed into a filling medium set BIT3′ 中心 ,in, TBtarget 中心 Indicates the content of the first filling medium, MORPH dmax OSB represents the set of results from the opening operation of the sphere structuring element with the largest diameter dmax in the first set of pores. dmax Represents Onum dmax First pore set PORE 中心 The proportion of the number of elements, Onum dmax OSB represents the number of elements in the set of sphere structure elements satisfying the maximum diameter dmax based on image morphology within the first set of pores; i MORPH i The elements in the first pore set PORE 中心 Percentage of elements; OSB i-1 MORPH i-1 The elements in the first pore set PORE 中心 The percentage of elements in the total number of elements; DSB1 indicates that the elements in DIST1 are in the first pore set PORE. 中心 The proportion of the number of elements, DIST1 represents DMORPH i-1 A set consisting of elements whose Euclidean distance is 1; DMORPH i-1 This represents the difference set of the sets of sphere structural elements with diameters of i-1 and i within the first set of pores, satisfying the image morphology-based opening operation. (DSB) j DIST j The elements in the first pore set PORE 中心 The proportion of elements in the total number; DIST j+1 Indicates DMORPH i-1 The set consisting of elements whose Euclidean distance is j+1; DSB j+1 DIST j+1 The elements in the first pore set PORE 中心 The percentage of elements; nump 中心 PORE represents the first set of pores. 中心 The number of elements in the middle.

2. The method for constructing a three-dimensional digital core model with medium filling according to claim 1, characterized in that, The initial three-dimensional digital core model filled with the pore centers of the medium is represented by the first initial three-dimensional digital core set CORE1, wherein... In the formula, MATRIX 中心 PORE represents the first skeleton set. 中心 Let U denote the first set of pores, and U denote the union of the first set of skeletons, MATRIX. 中心 Represented as: First Pore Set PORE 中心 Represented as: In the formula, px 中心 This represents an element in the initial three-dimensional digital core model filled with pore centers of the medium, where phase represents the phase state function, phase(px 中心 )=0 means px 中心 The element is a skeleton, phase(px) 中心 )=1 means px 中心 The element is porosity.

3. The method for constructing a three-dimensional digital core model with medium filling according to claim 2, characterized in that, The three-dimensional digital core model of medium filling the pore center includes a pore space model of medium gradually filling the pore center. The pore space model of medium gradually filling the pore center is obtained through the center-filled three-dimensional digital core assembly CORE. 中心 Specifically, it means: Among them, CORE 中心 MATRIX represents a three-dimensional digital core assembly filled with media at its center. 中心 Represents the first skeleton set, PORE 中心 Represents the first set of pores, BIT 中心 Represents the centrally filled medium set, phase(px) 中心 =3 means px 中心 The element is the filling medium, and phase represents the phase function.

4. The method for constructing a three-dimensional digital core model with medium filling according to claim 3, characterized in that, If OSB dmax >TBtarget 中心 Then MORPH dmax All the pores in the middle are transformed into the filling medium, and the central filling medium set is BIT. 中心 Represented as: ， The number of elements in the filling medium set is Onum dmax ; In the formula, TBtarget 中心 Indicates the content of the first filling medium, MORPH dmax Onum represents the set of results of the opening operation on the spherical structuring element with the largest diameter dmax in the first set of pores. dmax OSB represents the number of elements in the set of sphere structure elements satisfying the maximum diameter dmax based on image morphology within the first set of pores; dmax Represents Onum dmax First pore set PORE 中心 The proportion of the number of elements; OPEN represents the opening operator, px 中心 This represents the elements in the initial three-dimensional digital core model filled with pore centers of the medium, where dmax represents the maximum diameter of the spherical structural element, and OPEN (px) represents the maximum diameter of the spherical structural element. 中心 dmax) = 1 means px 中心 The result of the opening operation of the sphere structuring element with the maximum diameter dmax.

5. The method for constructing a three-dimensional digital core model with medium filling according to claim 3, characterized in that, If OSB i <TBtarget 中心 <OSB i-1 And DSB1>TBtarget 中心 -OSB i When i = 2, 3, 4, ..., dmax, then MORPH i All the pores in the medium are transformed into a filling medium. The set of central filling media is represented as follows: In the formula, BIT 中心 BIT1 represents the centrally filled medium set. 中心 Indicates center-filled medium set 1, BIT2 中心 Indicates center-filled medium set 2, MORPH i OSB represents the set of results of the opening operation on the sphere structuring element with diameter i in the first pore set. i MORPH i The elements in the first pore set PORE 中心 The proportion of elements in the total number; MORPH i-1 OSB represents the set of results of the opening operation on sphere structuring elements with diameter i-1 in the first set of pores. i-1 MORPH i-1 The elements in the first pore set PORE 中心 Percentage of elements; px 中心 This represents an element in the initial three-dimensional digital core model filled with pores in the medium. OPEN represents the opening operator, and i and i-1 represent the diameters of the spherical structural elements, respectively. OPEN(px 中心 i) = 1 indicates that px 中心 The result of the opening operation of a sphere structuring element with diameter i is satisfied; DMORPH i-1 DIST1 represents the difference set of the results of the opening operation between sphere structural elements with diameters of i-1 and i in the first set of pores, based on image morphology; DIST1 represents DMORPH. i-1 The set consisting of elements in the set whose Euclidean distance is 1; DSB1 represents the set of elements in DIST1 in the first pore set PORE. 中心 The percentage of elements; dist represents the Euclidean distance operator.

6. The method for constructing a three-dimensional digital core model with medium filling according to claim 3, characterized in that, If OSB i <TBtarget 中心 <OSB i-1 and When i = 2, 3, 4, ..., dmax, The set of central filling media is represented as follows: In the formula, BIT 中心 BIT1 represents the centrally filled medium set. 中心 Indicates center-filled medium set 1', BIT2 中心 Indicates the central filling medium set 2', BIT3 中心 Indicates center-filled medium set 3', MORPH i MORPH represents the set of results of the opening operation on the spherical structuring element with diameter i in the first set of pores. i-1 DMORPH represents the set of results of the opening operation on spherical structuring elements with diameter i-1 in the first set of pores. i-1 OSB represents the difference set of the results of the opening operation between sphere structuring elements of diameters i-1 and i in the first set of pores, based on image morphology. i MORPH i The elements in the first pore set PORE 中心 Percentage of elements in the total number of elements; OSB i-1 MORPH i-1 The elements in the first pore set PORE 中心 Percentage of elements; px 中心 This represents an element in the initial three-dimensional digital core model filled with pores in the medium. OPEN represents the opening operator, and i and i-1 represent the diameters of the spherical structural elements, respectively. OPEN(px 中心 (i) = 1 indicates that the center of px satisfies the opening operation result of a sphere structuring element with diameter i; DIST j Indicates DMORPH i-1 The set of elements whose Euclidean distance is j; DSB j DIST j The elements in the first pore set PORE 中心 The proportion of elements in the total number; DIST j+1 Indicates DMORPH i-1 The set consisting of elements whose Euclidean distance is j+1; DSB j+1 DIST j+1 The elements in the first pore set PORE 中心 The percentage of elements; dist represents the Euclidean distance operator; dist(px 中心 ) indicates the calculation element px 中心 Euclidean distance; ∪ denotes the union, ∪ j=1 i DIST j Let DIST1, DIST2, ..., DIST be the set. i The union of .

7. The method for constructing a three-dimensional digital core model with medium filling according to claim 1, characterized in that, Constructing a three-dimensional digital core model of the medium filling the pore edges includes: Based on the three-dimensional digital core model that identifies the skeleton and pores, an initial three-dimensional digital core model with pore edge filling is constructed. Based on the initial three-dimensional digital core model of the filling medium distribution morphology, filling medium content and filling of the pore edges, a three-dimensional digital core model of the medium filling the pore edges is constructed using the image Euclidean distance map method.

8. The method for constructing a three-dimensional digital core model with medium filling according to claim 7, characterized in that, The method of constructing a three-dimensional digital core model of the medium filling the pore edges using an image Euclidean distance map includes: Based on the initial three-dimensional digital core model filled with the pore edges of the medium, the second skeleton set and the second pore set are obtained, and the number of elements in the second pore set is obtained. The Euclidean distance method is used to calculate the Euclidean distance from the elements in the second pore set to the second skeleton set. Based on the Euclidean distance from the elements in the second pore set to the second skeleton set, the second pore set is divided into multiple second Euclidean distance sets, and the number of elements in each second Euclidean distance set and the proportion of the number of elements in each Euclidean distance set in the second pore set are obtained. Based on multiple sets of second Euclidean distances, the proportion of the number of elements in each Euclidean distance set in the second pore set, and the content of the filling medium, a three-dimensional digital core model of the filling at the edge of the medium is constructed.

9. The method for constructing a three-dimensional digital core model with medium filling according to claim 8, characterized in that, The initial three-dimensional digital core model filled with the pore edges of the medium is represented by the second initial three-dimensional digital core set CORE2, wherein... In the formula, MATRIX 边缘 PORE represents the second skeleton set. 边缘 Let U denote the second set of pores, and U denote the union of the second skeleton set, MATRIX. 边缘 Represented as: Second pore set PORE 边缘 Represented as: In the formula, px 边缘 This represents an element in a three-dimensional digital core set filled with pore edges of the medium, where phase represents the phase state function, phase(px 边缘 )=0 means px 边缘 The element is a skeleton, phase(px) 边缘 )=1 means px 边缘 The element is porosity.

10. The method for constructing a three-dimensional digital core model with medium filling according to claim 8, characterized in that, The three-dimensional digital core model of the medium filling the pore edges includes a pore space model in which the medium gradually fills the pore edges, and a pore space model in which the medium gradually fills the pore edges is connected to the center-filled three-dimensional digital core model CORE. 边缘 Specifically, it means: Among them, CORE 边缘 This represents a 3D digital core model after edge filling; MATRIX 边缘 Represents the second skeleton set; PORE 边缘 Represents the second set of pores; BIT 边缘 Represents the set of edge-filling media, phase(px) 边缘 )=3 indicates that the element is used to fill the medium; px 边缘 The element represents the initial three-dimensional digital core model filled with pore edges of the medium; phase represents the phase function.

11. The method for constructing a three-dimensional digital core model with medium filling according to claim 9, characterized in that, If PDSB1>TBtarget 边缘 Then, nump is randomly selected from the PDIST1 set. 边缘 TBtarget 边缘 Each pore element is transformed into a filling medium, and the edge filling medium set BIT in the three-dimensional digital core model of the medium edge filling is established. 边缘 , ， The number of elements in the filling medium set is nump 边缘 TBtarget 边缘 ; In the formula, BIT 边缘 TBtarget represents the edge-filled media set. 边缘 Indicates the content of the edge-filling medium, nump 边缘 PORE represents the second set of pores filled at the edge. 边缘 The number of elements, in pixels (px) 边缘 This represents an element in the initial three-dimensional digital core model filled with pore edges of the medium; dist represents the Euclidean distance operator; PDIST1 represents the element in the second pore set PORE. 边缘 The second Euclidean distance set consists of elements with a distance of 1 in the middle Euclidean distance set. PDSB1 represents the proportion of the number of elements with a distance of 1 in the second Euclidean distance set in the second pore set.

12. The method for constructing a three-dimensional digital core model with medium filling according to claim 9, characterized in that, when hour, The edge filling medium set is as follows: , ； In the formula, BIT 边缘 BIT1 represents the set of edge-filling media. 边缘 BIT2 represents the edge-filling medium set 1. 边缘 Indicates filling medium set 2, TBtarget 边缘 Indicates the content of edge-filling medium, px 边缘 This represents an element in the initial three-dimensional digital core model filled with pore edges of the medium; `dist` represents the Euclidean distance operator; `PDIST` j Indicated in the second pore set PORE 边缘 The second Euclidean distance set consisting of elements with a medium Euclidean distance of j, PDSB j PDIST represents the set of Euclidean distances j. j The proportion of medium elements in the second pore set, PDIST i+1 Indicated in the second pore set PORE 边缘 The second Euclidean distance set, PDSB, is composed of elements with a median Euclidean distance of i+1. i+1 PDIST represents the set with Euclidean distance i+1. i+1 The proportion of the number of elements in the second pore set, ∪ j=1 i PDIST j Let PDIST1, PDIST2, ..., PDIST be the set. i The union of .

13. The method for constructing a three-dimensional digital core model with medium filling according to claim 1, characterized in that, The construction of a medium-filled three-dimensional digital core model based on a three-dimensional digital core model with the medium filling the pore center and a three-dimensional digital core model with the medium filling the pore edge includes: First, based on the three-dimensional digital core model with the medium filling the center of the pores, the center is filled to a predetermined range. Then, based on the three-dimensional digital core model with the medium filling the edge of the pores, the edge filling is completed to obtain the three-dimensional digital core model with the medium filling.

14. The method for constructing a three-dimensional digital core model with medium filling according to any one of claims 1 or 13, characterized in that, The construction of a medium-filled three-dimensional digital core model based on a three-dimensional digital core model with the medium filling the pore center and a three-dimensional digital core model with the medium filling the pore edge further includes: First, based on the three-dimensional digital core model with the medium filling the pore edges, the edge filling is carried out to a predetermined range. Then, based on the three-dimensional digital core model with the medium filling the pore center, the center filling is completed to obtain the three-dimensional digital core model with the medium filling.

15. A device for constructing a three-dimensional digital core model filled with a medium, characterized in that, The device includes, The acquisition module is used to obtain the type, distribution pattern, and content of the core filling medium based on two-dimensional core images; The first construction module is used to construct, based on the distribution morphology and content of the filling medium, a three-dimensional digital core model that identifies the skeleton and pores, and a three-dimensional digital core model that shows the medium filling the pore center and the medium filling the pore edge. The three-dimensional digital core model, in which the construction medium fills the pore centers, includes: A three-dimensional digital core model that identifies the skeleton and pores is used as the initial three-dimensional digital core model for filling the pore center of the medium. Based on the initial three-dimensional digital core model filled with the pore center of the medium, the first skeleton set and the first pore set are obtained, and the number of elements in the first pore set is obtained. Based on image morphology, an opening operation is performed on spherical structural elements of different diameters in the first pore set to obtain the first opening operation result set, and the number of first pore elements in the first opening operation result set and the proportion of the number of first pore elements in the first pore set are obtained. Based on the first set of opening operation results, obtain the difference set of the opening operation results sets of two adjacent diameter sphere structure elements; Perform Euclidean distance calculation on the difference set to obtain the number of pore elements in the corresponding Euclidean distance set and the proportion of the corresponding number of pore elements in the first pore set of the initial three-dimensional digital core model filled with the pore center of the medium. Based on the opening operation result sets of spherical structural elements of different diameters, the number of pore elements in each diameter's opening operation result set, the proportion of the number of pore elements in each diameter's opening operation result set in the pore set of the three-dimensional digital core model filled with medium pore centers, the difference set of the opening operation result sets of two adjacent diameter spherical structural elements, the Euclidean distance set in the difference set of the opening operation result sets of adjacent diameter spherical structural elements, the proportion of the number of pore elements in each Euclidean distance in the difference set in the pore set of the initial three-dimensional digital core model filled with medium pore centers, and the filling medium content, a three-dimensional digital core model filled with medium pore centers is constructed; where, if OSB dmax >TBtarget 中心 Then MORPH dmax All the pores in the medium are transformed into a filling medium; If OSB i <TBtarget 中心 <OSB i-1 And DSB1>TBtarget 中心 -OSB i When i = 2, 3, 4, ..., dmax, then MORPH i All the pores in the medium are transformed into a filling medium; If OSB i <TBtarget 中心 <OSB i-1 and When i = 2, 3, 4, ..., dmax, first, all elements in the opening operation set MORPHi of the sphere structure elements with diameter i are transformed into the filling medium set BIT1′. 中心 Secondly, the first open operation set MORPH... i-1 The difference set of MORPHi, DMORPH i-1 All elements in the set whose distance is less than or equal to i fill the medium set BIT2′. 中心 Finally, in the first open set MORPH i-1 MORPH i DMORPH difference set i-1 The set whose medium-Euclidean distance is equal to i+1 is randomly selected. Each element is transformed into a filling medium set BIT3′ 中心 , in, TBtarget 中心 Indicates the content of the first filling medium, MORPH dmax OSB represents the set of results from the opening operation of the sphere structuring element with the largest diameter dmax in the first set of pores. dmax Represents Onum dmax First pore set PORE 中心 The proportion of the number of elements, Onum dmax OSB represents the number of elements in the set of sphere structure elements satisfying the maximum diameter dmax based on image morphology within the first set of pores; i MORPH i The elements in the first pore set PORE 中心 Percentage of elements; OSB i-1 MORPH i-1 The elements in the first pore set PORE 中心 The percentage of elements in the total number of elements; DSB1 indicates that the elements in DIST1 are in the first pore set PORE. 中心 The proportion of the number of elements, DIST1 represents DMORPH i-1 A set consisting of elements whose Euclidean distance is 1; DMORPH i-1 This represents the difference set of the sets of sphere structural elements with diameters of i-1 and i within the first set of pores, satisfying the image morphology-based opening operation. (DSB) j DIST j The elements in the first pore set PORE 中心 The proportion of elements in the total number; DIST j+1 Indicates DMORPH i-1 The set consisting of elements whose Euclidean distance is j+1; DSB j+1 DIST j+1 The elements in the first pore set PORE 中心 The percentage of elements; nump 中心 PORE represents the first set of pores. 中心 Number of elements in the middle; The second construction module is used to construct a medium-filled three-dimensional digital core model based on a three-dimensional digital core model with the medium filling the pore center and a three-dimensional digital core model with the medium filling the pore edge.

Images