Two-dimensional digital core pore filling method and device, electronic equipment and medium

Abstract

The invention discloses a two-dimensional digital core pore filling method and device, electronic equipment and a medium. The method comprises the following steps: performing cutting and data preprocessing on three-dimensional digital core data to obtain a plurality of two-dimensional digital core images as preparation data; obtaining a binary segmentation image according to the preparation data, and taking the binary segmentation image as use data including an initial core image and porosity; performing connected component analysis on the initial rock core image, obtaining a marked image, and calculating the number of connected regions; carrying out random sorting and random selection on the connected regions to obtain a connected region array; and filling the pores based on the initial core image, the connected region array and the filling proportion. On the basis of an image processing algorithm, the digital core pore filling method is provided by combining multiple algorithms, and the pore filling proportion and the filling mode can be manually controlled.

Description

Technical Field

[0001] The present invention relates to the fields of digital rock physics and image processing, and more specifically, to a two-dimensional digital core pore filling method, device, electronic equipment and medium. Background Art

[0002] Rock filling is a common phenomenon in geological processes, which occurs when cracks and pores in the stratum are filled with various substances. In oil and gas exploration, rock filling may bring a series of problems, affecting the effectiveness and efficiency of exploring and developing oil and gas resources. Traditional rock physics experimental methods are difficult to study specific filling materials and filling types. Digital core technology is an effective method to study rock filling characteristics. Digital core technology is a method to digitally represent rock samples under certain geological conditions through mathematical modeling methods, computer image processing and high-resolution CT imaging, and to obtain core properties by using numerical simulation and other methods to equate physical parameters to the constructed digital core model. Digital core technology can more accurately calculate rock electrical properties, elastic properties, seepage properties and other property parameters, providing a new choice for rock feature analysis. Image processing and morphological methods have a wide range of applications and good application effects. Combining them with digital rock physics is an effective way to process cores. Therefore, it is of practical significance to deal with rock karst and filling problems based on image morphological methods.

[0003] There are few image processing algorithms in the field of digital rock physics, most of which are grayscale segmentation, edge extraction and other algorithms. The data obtained are mostly binary segmentation images, and then numerical simulation is carried out to calculate the attribute parameters of digital core acoustics, seepage and other aspects. The above methods are difficult to meet the actual needs of filling and modeling karst holes. In karst and reservoirs with strong filling effects, there is a problem of deviation between actual drilling and demonstration. The drilling results reveal that the reservoir is highly heterogeneous, severely filled, and the prediction consistency rate is low. It is urgent to carry out rock physical testing and modeling related research with different filling methods and filling types. Traditional rock physics methods are difficult to study partial areas of rocks, and it is difficult to explore the impact of local factors on the overall rock properties.

[0004] There is still a need to develop a two-dimensional digital core pore filling method.

[0005] The information disclosed in the background technology section of the present invention is only intended to deepen the understanding of the general background technology of the present invention, and should not be regarded as acknowledging or suggesting in any form that the information constitutes the prior art already known to those skilled in the art. Summary of the invention

[0006] The present invention proposes a two-dimensional digital core pore filling method, device, electronic equipment and medium. Based on an image processing algorithm, a plurality of algorithms are combined to propose a digital core pore filling method, which can artificially control the pore filling ratio and filling method.

[0007] In a first aspect, an embodiment of the present disclosure provides a two-dimensional digital core pore filling method, comprising:

[0008] Cut and preprocess the three-dimensional digital core data to obtain multiple two-dimensional digital core images as preliminary data;

[0009] Acquire a binary segmentation image according to the prepared data as use data, including an initial core image and porosity;

[0010] Performing connected component analysis on the initial core image, obtaining a labeled image, and calculating the number of connected regions;

[0011] Randomly sort and randomly select the connected areas to obtain the connected area array A random ;

[0012] The pores are filled based on the initial core image, the connected region array and the filling ratio.

[0013] As a specific implementation manner of the embodiment of the present disclosure, the preliminary data is segmented by a threshold segmentation method to obtain the binary segmented image.

[0014] As a specific implementation of the embodiment of the present disclosure, the porosity is obtained by the following steps:

[0015] Segmentation is performed at each pixel position (i, j) in the initial core image, where the grayscale value of the pore pixel is set to 0 and the total pore pixel is S 1 , the skeleton pixel gray value is set to 2, and the total skeleton pixel value is recorded as S 2 , and then calculate the porosity.

[0016] As a specific implementation of the embodiment of the present disclosure, the porosity is:

[0017] P porosity =S 1 / (S 1 +S 2 )

[0018] Among them, P porosity is the porosity.

[0019] As a specific implementation of the embodiment of the present disclosure, the number of connected areas is:

[0020] Nfeatures =label(I initial =0)

[0021] Among them, N features is the number of connected regions, I initial is the initial core image, label() is the labeling function, I initial =0 indicates the area with 0 pixels, that is, the pore area.

[0022] As a specific implementation of the embodiment of the present disclosure, the connected region array is:

[0023] A random =random_permutatoin(S size )

[0024] Among them, A random is the connected region array, random_permutatoin() is the random sorting function, S size is the size of the connected area.

[0025] As a specific implementation method of the embodiment of the present disclosure, the pores are filled:

[0026] I final =Fill_Porosity(I initial , A random , P fill )

[0027] Among them, P fill is the filling ratio, I final For the filled core, Fill_Porosity() is the digital core image filling function.

[0028] In a second aspect, the present disclosure also provides a two-dimensional digital core pore filling device, including:

[0029] The data preparation module performs cutting and data preprocessing on the three-dimensional digital core data to obtain multiple two-dimensional digital core images as preparation data;

[0030] A binary segmentation module, which obtains a binary segmentation image according to the prepared data as use data, including an initial core image and porosity;

[0031] A connectivity analysis module performs connectivity component analysis on the initial core image, obtains a labeled image, and calculates the number of connected regions;

[0032] The sorting and selection module randomly sorts and randomly selects the connected areas to obtain an array of connected areas;

[0033] A filling module fills the pores based on the initial core image, the connected region array and the filling ratio.

[0034] As a specific implementation manner of the embodiment of the present disclosure, the preliminary data is segmented by a threshold segmentation method to obtain the binary segmented image.

[0035] As a specific implementation of the embodiment of the present disclosure, the porosity is obtained by the following steps:

[0036] Segmentation is performed at each pixel position (i, j) in the initial core image, where the grayscale value of the pore pixel is set to 0 and the total pore pixel is S 1 , the skeleton pixel gray value is set to 2, and the total skeleton pixel value is recorded as S 2 , and then calculate the porosity.

[0037] As a specific implementation of the embodiment of the present disclosure, the porosity is:

[0038] P porosity =S 1 / (S 1 +S 2 )

[0039] Among them, P porosity is the porosity.

[0040] As a specific implementation of the embodiment of the present disclosure, the number of connected areas is:

[0041] N features =label(I initial =0)

[0042] Among them, N features is the number of connected regions, I initial is the initial core image, label() is the labeling function, I initial =0 indicates the area with 0 pixels, that is, the pore area.

[0043] As a specific implementation of the embodiment of the present disclosure, the connected region array is:

[0044] A random =random_permutatoin(S size )

[0045] Among them, A random is the connected region array, random_permutatoin() is the random sorting function, S size is the size of the connected area.

[0046] As a specific implementation method of the embodiment of the present disclosure, the pores are filled:

[0047] I final =Fill_Porosity(I initial , A random , P fill )

[0048] Among them, P fill is the filling ratio, I final For the filled core, Fill_Porosity() is the digital core image filling function.

[0049] In a third aspect, an embodiment of the present disclosure further provides an electronic device, the electronic device comprising:

[0050] A memory storing executable instructions;

[0051] A processor runs the executable instructions in the memory to implement the two-dimensional digital core pore filling method.

[0052] In a fourth aspect, an embodiment of the present disclosure further provides a computer-readable storage medium, which stores a computer program, and when the computer program is executed by a processor, the two-dimensional digital core pore filling method is implemented.

[0053] Its beneficial effects are:

[0054] The present invention sets a certain filling ratio to fill the pores, generates a digital core image with continuous or discontinuous fillers, and constructs a pore filling method by combining a variety of image morphology algorithms. This can more realistically simulate the pore structure of underground rock formations, and provides an effective method for exploring the impact of filling types on reservoirs, providing strong support for geological research and resource development.

[0055] The methods and apparatus of the present invention have other features and advantages that will be apparent from, or will be described in detail in, the accompanying drawings and subsequent detailed descriptions incorporated herein, which together serve to explain the specific principles of the present invention. BRIEF DESCRIPTION OF THE DRAWINGS

[0056] The above and other objects, features and advantages of the present invention will become more apparent through a more detailed description of exemplary embodiments of the present invention in conjunction with the accompanying drawings, wherein like reference numerals generally represent like components throughout the exemplary embodiments of the present invention.

[0057] Figure 1A flow chart showing the steps of a two-dimensional digital core pore filling method according to an embodiment of the present invention.

[0058] Figure 2 A schematic diagram showing digital core data according to an embodiment of the present invention is shown.

[0059] Figure 3a and Figure 3b Schematic diagrams showing random filling effects of original data and core data according to an embodiment of the present invention are shown respectively.

[0060] Figure 4a and Figure 4b Schematic diagrams showing the continuous filling effect of raw data and core data according to an embodiment of the present invention are shown respectively.

[0061] Figure 5 A block diagram of a two-dimensional digital core pore filling device according to an embodiment of the present invention is shown.

[0062] Description of reference numerals:

[0063] 201. Data preparation module; 202. Binary segmentation module; 203. Connectivity analysis module; 204. Sorting and selection module; 205. Filling module. DETAILED DESCRIPTION

[0064] The preferred embodiments of the present invention will be described in more detail below. Although the preferred embodiments of the present invention are described below, it should be understood that the present invention can be implemented in various forms and should not be limited to the embodiments set forth herein.

[0065] To facilitate understanding of the solutions and effects of the embodiments of the present invention, six specific application examples are given below. Those skilled in the art should understand that the examples are only for facilitating understanding of the present invention, and any specific details thereof are not intended to limit the present invention in any way.

[0066] Example 1

[0067] Figure 1 A flow chart showing the steps of a two-dimensional digital core pore filling method according to an embodiment of the present invention.

[0068] like Figure 1As shown, the two-dimensional digital core pore filling method includes: step 101, cutting and data preprocessing for three-dimensional digital core data, and obtaining multiple two-dimensional digital core images as preliminary data; step 102, obtaining a binary segmentation image according to the preliminary data as use data, including an initial core image and porosity; step 103, performing connected component analysis on the initial core image, obtaining a labeled image, and calculating the number of connected areas; step 104, randomly sorting and randomly selecting the connected areas to obtain an array of connected areas; step 105, filling the pores based on the initial core image, the connected area array and the filling ratio.

[0069] In one example, the prepared data is segmented by a threshold segmentation method to obtain a binary segmentation image.

[0070] In one example, the porosity is obtained by the following steps:

[0071] Segmentation is performed at each pixel position (i, j) in the initial core image, where the grayscale value of the pore pixel is set to 0 and the total pore pixel is S 1 , the skeleton pixel gray value is set to 2, and the total skeleton pixel value is recorded as S 2 , and then calculate the porosity.

[0072] In one example, the porosity is:

[0073] P porosity =S 1 / (S 1 +S 2 )

[0074] Among them, P porosity is the porosity.

[0075] In one example, the number of connected regions is:

[0076] N features =label(I initial =0)

[0077] Among them, N features is the number of connected regions, I initial is the initial core image, label() is the labeling function, I initial =0 indicates the area with 0 pixels, that is, the pore area.

[0078] In one example, the connected components array is:

[0079] A random =random_permutatoin(S size )

[0080] Among them, Arandom is the connected region array, random_permutatoin() is the random sorting function, S size is the size of the connected area.

[0081] In one example, the pores are filled:

[0082] I final =Fill_Porosity(I initial , A random , P fill )

[0083] Among them, P fill is the filling ratio, I final For the filled core, Fill_Porosity() is the digital core image filling function.

[0084] Specifically, the purpose of the present invention is to address the difficulty of conducting research on partial areas of rocks in conventional rock physics experiments, and therefore it is difficult to analyze the impact of the filling area on the overall properties of the core. With the support of high-precision X-ray CT scanning equipment, a random and continuously assembled pore filling method is provided around digital core technology. The filling method of the present invention acts on binary segmentation grayscale images commonly used in the digital core analysis process, combines and optimizes random selection algorithms, connected component analysis methods, random sorting and screening methods, grayscale value filling methods, and other methods to achieve a digital core model construction method with controllable filling ratios and filling types.

[0085] According to the three-dimensional digital core data obtained by high-precision CT scanning, cutting and data preprocessing are performed to obtain a large number of two-dimensional digital core images as preliminary data; the threshold segmentation method is used to accurately segment the two-dimensional digital core images to obtain binary segmentation images as the use data of the algorithm. The porosity is recorded as P porosity , the initial core image is denoted as I initial , in I initial Each pixel position (i, j) in is segmented, where the grayscale value of the pore pixel is set to 0 and the total pore pixel is S 1 , the skeleton pixel gray value is set to 2, and the total skeleton pixel value is recorded as S 2 The porosity calculation formula is:

[0086] P porosity =S 1 / (S 1 +S 2 )

[0087] The three-dimensional digital core data is relatively large and has a high demand for physical memory, and the numerical simulation calculations for the three-dimensional model require extremely large computing resources; the two-dimensional digital core properties can effectively reflect the underground medium conditions to a certain extent, and have certain practical significance, and numerical simulation is a process based on segmented images. Therefore, the use of two-dimensional segmented images as experimental and application data in the present invention is reasonable.

[0088] A pore filling algorithm based on random continuous components was constructed and tested using the prepared data. The working principle and implementation steps of the algorithm are as follows:

[0089] Connected component analysis is performed on the initial image. The initial core image is denoted as I initial , the marked image is I label , the number of connected regions is N features , using connected component analysis to label I initial The pore-connected area in the image is generated to generate the labeled image I label And calculate the number of connected regions N features . Among them I label , N features =label(I initial =0).

[0090] Randomly sort and select the connected areas. size , randomly selected connected area array A random , randomly sort the connected area size array to obtain a randomly selected connected area array. random =random_permutatoin(S size ).

[0091] Fill the pores. Note the filling ratio P fill , the core after filling is I final , select an area from the randomly connected area and randomly fill some pixels in the selected area until the predetermined filling ratio P is reached fill , where I final =Fill_Porosity(I initial , A random , P fill ).

[0092] Example 2

[0093] The present invention also provides a two-dimensional digital core pore filling device, comprising:

[0094] The data preparation module performs cutting and data preprocessing on the three-dimensional digital core data to obtain multiple two-dimensional digital core images as preparation data;

[0095] A binary segmentation module, which obtains a binary segmentation image based on the prepared data as the use data, including an initial core image and porosity;

[0096] The connectivity analysis module performs connected component analysis on the initial core image, obtains the labeled image, and calculates the number of connected regions;

[0097] The sorting and selection module randomly sorts and randomly selects the connected areas to obtain an array of connected areas;

[0098] The filling module fills the pores based on the initial core image, the connected region array and the filling ratio.

[0099] In one example, the prepared data is segmented by a threshold segmentation method to obtain a binary segmentation image.

[0100] In one example, the porosity is obtained by the following steps:

[0101] Segmentation is performed at each pixel position (i, j) in the initial core image, where the grayscale value of the pore pixel is set to 0 and the total pore pixel is S 1 , the skeleton pixel gray value is set to 2, and the total skeleton pixel value is recorded as S 2 , and then calculate the porosity.

[0102] In one example, the porosity is:

[0103] P porosity =S 1 / (S 1 +S 2 )

[0104] Among them, P porosity is the porosity.

[0105] In one example, the number of connected regions is:

[0106] N features =label(I initial =0)

[0107] Among them, N features is the number of connected regions, I initial is the initial core image, label() is the labeling function, I initial =0 indicates the area with 0 pixels, that is, the pore area.

[0108] In one example, the connected components array is:

[0109] A random =random_permutatoin(S size)

[0110] Among them, A random is the connected region array, random_permutatoin() is the random sorting function, S size is the size of the connected area.

[0111] In one example, the pores are filled:

[0112] I final =Fill_Porosity(I initial , A random , P fill )

[0113] Among them, P fill is the filling ratio, I final For the filled core, Fill_Porosity() is the digital core image filling function.

[0114] Specifically, the three-dimensional digital core data obtained by high-precision CT scanning is cut and preprocessed to obtain a large number of two-dimensional digital core images as preliminary data; the threshold segmentation method is used to accurately segment the two-dimensional digital core images to obtain binary segmentation images as the use data of the algorithm. The porosity is recorded as P porosity , the initial core image is denoted as I initial , in I initial Each pixel position (i, j) in is segmented, where the grayscale value of the pore pixel is set to 0 and the total pore pixel is S 1 , the skeleton pixel gray value is set to 2, and the total skeleton pixel value is recorded as S 2 The porosity calculation formula is:

[0115] P porosity =S 1 / (S 1 +S 2 )

[0116] The three-dimensional digital core data is relatively large and has a high demand for physical memory, and the numerical simulation calculations for the three-dimensional model require extremely large computing resources; the two-dimensional digital core properties can effectively reflect the underground medium conditions to a certain extent, and have certain practical significance, and numerical simulation is a process based on segmented images. Therefore, the use of two-dimensional segmented images as experimental and application data in the present invention is reasonable.

[0117] A pore filling algorithm based on random continuous components was constructed and tested using the prepared data. The working principle and implementation steps of the algorithm are as follows:

[0118] Connected component analysis is performed on the initial image. The initial core image is denoted as Iinitial , the marked image is I label , the number of connected regions is N features , using connected component analysis to label I initial The pore-connected area in the image is generated to generate the labeled image I label And calculate the number of connected regions N features . Among them I label , N features =label(I initial =0).

[0119] Randomly sort and select the connected areas. size , randomly selected connected area array A random , randomly sort the connected area size array to obtain a randomly selected connected area array. random =random_permutatoin(S size ).

[0120] Fill the pores. Note the filling ratio P fill , the core after filling is I final , select an area from the randomly connected area and randomly fill some pixels in the selected area until the predetermined filling ratio P is reached fill , where I final =Fill_Porosity(I initial , A random , P fill ).

[0121] Example 3

[0122] The present invention selects a real drilled porous carbonate rock core, digitizes it through a high-resolution CT scanning method, and then produces test data to demonstrate the effect of the pore filling algorithm.

[0123] Figure 2 A schematic diagram showing digital core data according to an embodiment of the present invention is shown.

[0124] like Figure 2 As shown in the figure, the data is obtained by high-precision CT scanning and reconstruction of porous carbonate rocks, and two segmented images (pores and matrix) are obtained through multiple image preprocessing methods such as threshold segmentation.

[0125] High-precision CT scanning and reconstruction are performed on porous carbonate cores to obtain a three-dimensional digital core model. The obtained two-dimensional digital core image is cut and preprocessed to obtain a qualified two-dimensional digital core image as test data.

[0126] Figure 3a and Figure 3b Schematic diagrams of random filling effects of original data and core data according to an embodiment of the present invention are respectively shown, wherein the pore area (black part) is covered by the random filling area (grey part).

[0127] Figure 4a and Figure 4b Schematic diagrams of continuous filling effects of original data and core data according to an embodiment of the present invention are respectively shown, wherein the pore area (black part) is covered by the continuous filling area (grey part).

[0128] This method is used to test the test data, and the test results are as follows Figure 3a , Figure 3b , Figure 4a , Figure 4b As shown, Figure 4a , Figure 4b The pore interconnected areas are completely filled until the preset filling ratio is reached. The filling ratio and filling type requirements can be set artificially, and the accuracy meets the needs of actual work.

[0129] Example 4

[0130] Figure 5 A block diagram of a two-dimensional digital core pore filling device according to an embodiment of the present invention is shown.

[0131] like Figure 5 As shown, the two-dimensional digital core pore filling device comprises:

[0132] The data preparation module 201 performs cutting and data preprocessing on the three-dimensional digital core data to obtain a plurality of two-dimensional digital core images as preparation data;

[0133] A binary segmentation module 202, which obtains a binary segmentation image according to the prepared data as the use data, including an initial core image and porosity;

[0134] Connectivity analysis module 203, performs connected component analysis on the initial core image, obtains the labeled image, and calculates the number of connected regions;

[0135] The sorting and selecting module 204 randomly sorts and randomly selects the connected regions to obtain a connected region array;

[0136] The filling module 205 fills the pores based on the initial core image, the connected region array and the filling ratio.

[0137] As an optional solution, the prepared data is segmented by a threshold segmentation method to obtain a binary segmentation image.

[0138] Alternatively, obtain the porosity by:

[0139] Segmentation is performed at each pixel position (i, j) in the initial core image, where the grayscale value of the pore pixel is set to 0 and the total pore pixel is S 1 , the skeleton pixel gray value is set to 2, and the total skeleton pixel value is recorded as S 2 , and then calculate the porosity.

[0140] As an option, the porosity is:

[0141] P porosity =S 1 / (S 1 +S 2 )

[0142] Among them, P porosity is the porosity.

[0143] As an option, the number of connected regions is:

[0144] N features =label(I initial =0)

[0145] Among them, N features is the number of connected regions, I initial is the initial core image, label() is the labeling function, I initial =0 indicates the area with 0 pixels, that is, the pore area.

[0146] As an alternative, the connected area array is:

[0147] A random =random_permutatoin(S size )

[0148] Among them, A random is the connected region array, random_permutatoin() is the random sorting function, S size is the size of the connected area.

[0149] As an option, fill the pores:

[0150] I final =Fill_Porosity(I initial , A random , P fill )

[0151] Among them, P fill is the filling ratio, I finalFor the filled core, Fill_Porosity() is the digital core image filling function.

[0152] Example 5

[0153] This embodiment provides an electronic device, which includes: a memory storing executable instructions; and a processor, which runs the executable instructions in the memory to implement the above-mentioned two-dimensional digital core pore filling method.

[0154] An electronic device according to an embodiment of the present disclosure includes a memory and a processor.

[0155] The memory is used to store non-temporary computer-readable instructions. Specifically, the memory may include one or more computer program products, which may include various forms of computer-readable storage media, such as volatile memory and / or non-volatile memory. The volatile memory may, for example, include random access memory (RAM) and / or cache memory (cache), etc. The non-volatile memory may, for example, include read-only memory (ROM), hard disk, flash memory, etc.

[0156] The processor may be a central processing unit (CPU) or other forms of processing units having data processing capabilities and / or instruction execution capabilities, and may control other components in the electronic device to perform desired functions. In one embodiment of the present disclosure, the processor is used to run the computer-readable instructions stored in the memory.

[0157] Those skilled in the art should be able to understand that in order to solve the technical problem of how to obtain a good user experience, the present embodiment may also include well-known structures such as a communication bus and an interface, and these well-known structures should also be included in the protection scope of the present disclosure.

[0158] For detailed description of this embodiment, reference may be made to the corresponding descriptions in the aforementioned embodiments, which will not be repeated here.

[0159] Example 6

[0160] This embodiment provides a computer-readable storage medium, which stores a computer program. When the computer program is executed by a processor, the two-dimensional digital core pore filling method is implemented.

[0161] According to the computer-readable storage medium of the embodiment of the present disclosure, non-transitory computer-readable instructions are stored thereon. When the non-transitory computer-readable instructions are executed by a processor, all or part of the steps of the above-mentioned methods of each embodiment of the present disclosure are executed.

[0162] The above-mentioned computer-readable storage media include, but are not limited to: optical storage media (e.g., CD-ROM and DVD), magneto-optical storage media (e.g., MO), magnetic storage media (e.g., magnetic tape or mobile hard disk), media with built-in rewritable non-volatile memory (e.g., memory card) and media with built-in ROM (e.g., ROM box).

[0163] Those skilled in the art should understand that the purpose of the above description of the embodiments of the present invention is only to exemplarily illustrate the beneficial effects of the embodiments of the present invention, and is not intended to limit the embodiments of the present invention to any given examples.

[0164] The embodiments of the present invention have been described above, and the above description is exemplary, not exhaustive, and is not limited to the disclosed embodiments. Many modifications and changes will be apparent to those skilled in the art without departing from the scope and spirit of the described embodiments.

Claims

1. A two-dimensional digital core pore filling method, It is characterized in that include: Cut and preprocess the three-dimensional digital core data to obtain multiple two-dimensional digital core images as preliminary data; Acquire a binary segmentation image according to the prepared data as use data, including an initial core image and porosity; Performing connected component analysis on the initial core image, obtaining a labeled image, and calculating the number of connected regions; Randomly sort and randomly select the connected areas to obtain the connected area array A random ; The pores are filled based on the initial core image, the connected region array and the filling ratio.

2. The two-dimensional digital core pore filling method according to claim 1, in, The prepared data is segmented by a threshold segmentation method to obtain the binary segmentation image.

3. The two-dimensional digital core pore filling method according to claim 1, in, The porosity is obtained by the following steps: Segmentation is performed at each pixel position (i, j) in the initial core image, where the grayscale value of the pore pixel is set to 0 and the total pore pixel is S 1 , the skeleton pixel gray value is set to 2, and the total skeleton pixel value is recorded as S 2 , and then calculate the porosity.

4. The two-dimensional digital core pore filling method according to claim 3, in, The porosity is: P porosity =S 1 / (S 1 +S 2 ) Among them, P porosity is the porosity.

5. The two-dimensional digital core pore filling method according to claim 1, in, The number of connected regions is: N features =label(I initial =0) Among them, N features is the number of connected regions, I initial is the initial core image, label() is the labeling function, I initial =0 indicates the area with 0 pixels, that is, the pore area.

6. The two-dimensional digital core pore filling method according to claim 1, in, The connected area array is: A random =random_permutatoin(S size ) Among them, A random is the connected region array, random_permutatoin() is the random sorting function, S size is the size of the connected area.

7. The two-dimensional digital core pore filling method according to claim 1, in, Filling the pores: I final =Fill_Porosity(I initial ,A random ,P fill ) Among them, P fill is the filling ratio, I final For the filled core, Fill_Porosity() is the digital core image filling function.

8. A two-dimensional digital core pore filling device, It is characterized in that include: The data preparation module performs cutting and data preprocessing on the three-dimensional digital core data to obtain multiple two-dimensional digital core images as preparation data; A binary segmentation module, which obtains a binary segmentation image according to the prepared data as use data, including an initial core image and porosity; A connectivity analysis module performs connectivity component analysis on the initial core image, obtains a labeled image, and calculates the number of connected regions; The sorting and selection module randomly sorts and randomly selects the connected areas to obtain an array of connected areas; A filling module fills the pores based on the initial core image, the connected region array and the filling ratio.

9. An electronic device, It is characterized in that The electronic device comprises: A memory storing executable instructions; A processor, wherein the processor runs the executable instructions in the memory to implement the two-dimensional digital core pore filling method according to any one of claims 1 to 7.

10. A computer-readable storage medium, It is characterized in that The computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the two-dimensional digital core pore filling method according to any one of claims 1 to 7 is implemented.

Images