Pore fractal dimension determination method, device, equipment, medium and program product

Abstract

The present application provides a kind of pore fractal dimension determination method, device, equipment, medium and program product, method includes: based on the microscopic image of shale surface random area, determine pore size distribution data;Based on pore size distribution data and cumulative size distribution algorithm, determine target pore fractal dimension;To target pore fractal dimension is representative test, in the case where test passes, determine target pore fractal dimension as shale pore fractal dimension.The present application determines shale pore size distribution data by analyzing the microscopic image of shale surface random area.Based on shale pore size distribution data and cumulative size distribution algorithm, determine target pore fractal dimension, and the representative test is carried out to the target pore fractal dimension obtained.It solves the technical problem that shale pore fractal dimension cannot be accurately obtained by the scheme of only randomly selecting part of the visual field of shale cross section for microscopic analysis when calculating pore fractal dimension by using the existing pore cumulative size distribution algorithm.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of image data processing, and in particular to a pore fractal dimension determination method, device, equipment, medium and program product. BACKGROUND

[0002] Quantitative research on shale pores by microscopic observation and in-depth analysis of the fractal characteristics of shale pores can provide important reference for efficient exploration and development of shale oil and gas. Based on microscopic observation technology, there are currently algorithms that can perform fractal analysis on the pores of two-dimensional sections of shale samples. The pore cumulative size distribution algorithm requires obtaining the total number of pores in all sections of the shale and pore size distribution data to calculate the fractal dimension of the shale pores. However, detailed microscopic analysis is time-consuming and laborious, and random selection of shale sections for microscopic analysis has limited error in small-scale pore fractal analysis, but the error increases in large-scale pore fractal analysis. Therefore, how to accurately obtain the fractal dimension of shale pores has become a technical problem to be solved. SUMMARY

[0003] The present application provides a pore fractal dimension determination method, device, equipment, medium and program product to solve the defect of inaccurate pore fractal dimension calculation in the existing pore fractal dimension analysis method, and realizes accurate calculation of the fractal dimension of shale pores.

[0004] The present application provides a pore fractal dimension determination method, comprising the following steps.

[0005] Based on the microscopic image of a random area on the surface of the shale, determine the pore size distribution data;

[0006] Based on the pore size distribution data and the cumulative size distribution algorithm, determine the target pore fractal dimension;

[0007] Representative test the target pore fractal dimension, and if the test is passed, determine the target pore fractal dimension as the fractal dimension of the shale pores.

[0008] According to the pore fractal dimension determination method provided by the present application, the determination of the target pore fractal dimension based on the pore size distribution data and the cumulative size distribution algorithm comprises:

[0009] Determine the target variable and the maximum pore diameter of the image based on the pore size distribution data;

[0010] Determine the target pore diameter interval based on the target variable and the maximum pore diameter of the image, and determine the first pore diameter interval based on the target variable;

[0011] Determine the target pore fractal dimension according to the target pore diameter interval, the first pore diameter interval and the cumulative size distribution algorithm.

[0012] The method for determining the pore fractal dimension according to the present application comprises the following steps:

[0013] Adding the fractal dimension assumption value into the adjustment term of the cumulative size distribution algorithm based on the target pore size interval and the first pore size interval;

[0014] Determining the fractal dimension response value corresponding to the fractal dimension assumption value through regression analysis;

[0015] Determining the target pore fractal dimension based on the two-dimensional variable relationship graph drawn based on the fractal dimension assumption value and the fractal dimension response value.

[0016] The method for determining the pore fractal dimension according to the present application comprises the following steps of determining the target pore fractal dimension based on the two-dimensional variable relationship graph drawn based on the fractal dimension assumption value and the fractal dimension response value:

[0017] Determining the target curve corresponding to the fractal dimension assumption value and the fractal dimension response value, and the target straight line where the fractal dimension assumption value is equal to the fractal dimension response value;

[0018] Determining the target value range corresponding to the fractal dimension assumption value;

[0019] Determining the target pore fractal dimension based on the target curve, the target straight line and the target value range.

[0020] The method for determining the pore fractal dimension according to the present application comprises the following steps of performing representative test on the target pore fractal dimension, and determining the target pore fractal dimension as the shale pore fractal dimension in the case where the test is passed:

[0021] Determining the cumulative observation area based on the microscopic image of the random area of the shale surface;

[0022] Determining the target two-dimensional variable relationship graph of the target pore fractal dimension and the cumulative observation area;

[0023] Performing representative test on the target pore fractal dimension based on the cumulative observation area and the target two-dimensional variable relationship graph;

[0024] Determining the target pore fractal dimension as the shale pore fractal dimension in the case where it is determined based on the cumulative observation area and the target two-dimensional variable relationship graph that the test is passed.

[0025] According to the pore fractal dimension determination method provided by the application, the target pore fractal dimension-target two-dimensional variable relationship graph is determined, and then the following steps are included:

[0026] The fluctuation value of the target pore fractal dimension is determined based on the target two-dimensional variable relationship graph.

[0027] The fluctuation value less than the target threshold in the continuous interval is determined as a target fluctuation value.

[0028] In the case that the cumulative observation area corresponding to the target fluctuation value is greater than a first threshold, it is determined that the representative test of the target pore fractal dimension is passed.

[0029] The application further provides a pore fractal dimension determination device, including the following modules:

[0030] A pore size distribution data determination module is configured to determine pore size distribution data based on a microscopic image of a random area on a shale surface.

[0031] A target pore fractal dimension determination module is configured to determine a target pore fractal dimension based on the pore size distribution data and a cumulative size distribution algorithm.

[0032] A representative test module is configured to perform a representative test on the target pore fractal dimension, and in the case that the test is passed, the target pore fractal dimension is determined as a shale pore fractal dimension.

[0033] The application further provides an electronic device including a memory, a processor, and a computer program stored in the memory and executable on the processor, and the processor implements the pore fractal dimension determination method according to any one of the above when executing the program.

[0034] The application further provides a non-transitory computer readable storage medium having a computer program stored thereon, and the computer program is executable on a processor to implement the pore fractal dimension determination method according to any one of the above.

[0035] The application further provides a computer program product including a computer program, and the computer program is executable on a processor to implement the pore fractal dimension determination method according to any one of the above.

[0036] The pore fractal dimension determination method, device, equipment, medium and program product provided by the present application determine shale pore size distribution data by analyzing microscopic images of random areas on a shale surface. Based on the shale pore size distribution data and a cumulative size distribution algorithm, a target pore fractal dimension is determined, and the obtained target pore fractal dimension is subjected to a representative test. The technical problem that the shale pore fractal dimension cannot be accurately obtained by only randomly selecting a part of the field of view of a shale cross section for microscopic analysis in the prior art method for calculating the pore fractal dimension is solved. BRIEF DESCRIPTION OF DRAWINGS

[0037] In order to more clearly illustrate the technical solutions in the present application or the prior art, the following will briefly introduce the drawings needed to be used in the embodiments or the prior art description. Obviously, the drawings in the following description are some embodiments of the present application, and other drawings can be obtained by those skilled in the art without any creative effort on the basis of these drawings.

[0038] Figure 1 is one of the flowcharts of the pore fractal dimension determination method provided by the present application.

[0039] Figure 2 is the second flowchart of the pore fractal dimension determination method provided by the present application.

[0040] Figure 3 is a structural schematic diagram of the pore fractal dimension determination device provided by the present application.

[0041] Figure 4 is a structural schematic diagram of the electronic device provided by the present application. DETAILED DESCRIPTION

[0042] In order to make the objects, technical solutions and advantages of the present application clearer, the technical solutions in the present application will be described clearly and completely below in combination with the drawings in the present application. Obviously, the described embodiments are some embodiments of the present application, but not all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without any creative effort fall within the protection scope of the present application.

[0043] The pore fractal dimension determination method, device, equipment, medium and program product of the present application will be described below in combination with Figures 1-4

[0044] Figure 1 is one of the flowcharts of the pore fractal dimension determination method provided by the present application, as shown in Figure 1 The method comprises the following steps:

[0045] ​Step 100, determining pore size distribution data based on a microscopic image of a random area of a shale surface;

[0046] Specifically, for the problem of large error of pore fractal dimension calculation results existing in the present random selection of shale interface microscopic analysis, the present application improves the accuracy of shale pore fractal dimension calculation through a new pore cumulative size distribution algorithm. First, a microscopic image of a random area of a shale surface is obtained, and the number of pores is extracted to determine the pore size distribution data.

[0047] The specific process is as follows:

[0048] Step 1, mechanically polishing and argon ion polishing the shale surface respectively to form a highly flat plane in the non-pore area of the shale surface to meet the requirements of pore microscopic observation.

[0049] Step 2, microscopic imaging of the random field of the shale surface is performed through microscopic observation, and the cumulative microscopic observation area is recorded; the microscopic observation includes but is not limited to scanning electron microscopy and helium ion microscopy.

[0050] Step 3, binarizing the microscopic image, extracting the pore area, pore number and pore size distribution data in the microscopic image based on the gray scale; the pore size refers to a parameter reflecting the pore size, including but not limited to the equivalent circle diameter and the Feret diameter.

[0051] Step 200, determining the target pore fractal dimension based on the pore size distribution data and the cumulative size distribution algorithm;

[0052] The specific process of calculating the pore fractal dimension through the new cumulative size distribution algorithm is as follows:

[0053] Step 4, determining the maximum pore diameter in the microscopic image, the target pore diameter interval, and the number of pores in the target pore diameter interval;

[0054] Step 5, determining the first pore diameter interval based on the target pore diameter interval and the number of pores in the target pore diameter interval, and the number of pores in the target pore diameter interval is greater than a certain threshold value;

[0055] Step 6, based on the target pore diameter interval, the number of pores in the target pore diameter interval, and the first pore diameter interval, adding a series of fractal dimension assumption values in the adjustment term of the general equation of the cumulative size distribution, and obtaining the fractal dimension response value corresponding to the fractal dimension assumption value through regression analysis;

[0056] Step 7, plot the curve corresponding to the fractal dimension hypothesis value and the fractal dimension response value in the two-dimensional variable relationship diagram, confirm that the curve and the straight line (indicating that the fractal dimension hypothesis value is equal to the fractal dimension response value) have only one intersection point in a certain range (for example, a range greater than 1 and less than 2) of the fractal dimension hypothesis value, and the fractal dimension hypothesis value (or the fractal dimension response value) corresponding to the intersection point is the actual calculation value of the fractal dimension, that is, the target pore fractal dimension in the embodiment.

[0057] Step 300, a representative test is performed on the target pore fractal dimension, and in the case that the test passes, the target pore fractal dimension is determined as the shale pore fractal dimension.

[0058] The present application also provides a method for performing a representative test on a fractal dimension, and the specific process is as follows:

[0059] Step 8, determine the cumulative observation area of microscopic observation, and the two-dimensional variable relationship diagram of the pore fractal dimension and the cumulative observation area;

[0060] Step 9, based on the two-dimensional variable relationship diagram of the pore fractal dimension and the cumulative observation area, determine the fluctuation value of the pore fractal dimension;

[0061] Step 10, in the continuous interval, if the fluctuation value is less than the target threshold value, the fluctuation value is determined as the target fluctuation value, and in the case that the cumulative observation area corresponding to the target fluctuation value is not less than the first threshold value, it is determined that the pore fractal dimension test passes. If the test does not pass, return to step 2.

[0062] The present embodiment determines the shale pore size distribution data by analyzing the microscopic image of the random area of the shale surface. Based on the shale pore size distribution data and the cumulative size distribution algorithm, the target pore fractal dimension is determined, and the obtained target pore fractal dimension is subjected to a representative test. The technical problem that the shale pore fractal dimension cannot be accurately obtained by only randomly selecting a part of the visual field of the shale section for microscopic analysis in the prior art is solved.

[0063] Figure 2 It is a flowchart of the pore fractal dimension determination method provided by the present application, as shown in Figure 2 The method can further include:

[0064] Step 210, determining a target variable and an image maximum pore diameter based on the pore size distribution data;

[0065] Step 220, determining a target pore diameter interval based on the target variable and the image maximum pore diameter, and determining a first pore diameter interval based on the target variable;

[0066] Step 230, determining the target pore fractal dimension according to the target pore size interval, the first pore size interval, and the cumulative size distribution algorithm.

[0067] The target pore size interval in step 4 above refers to a pore size range greater than or equal to the variable and less than the maximum pore size obtained in the micrograph ; the number of pores in the target pore size interval refers to the number of pores in the target pore size interval .

[0068] The first pore size interval in step 5 above refers to the pore size interval corresponding to the value range of the variable To ensure calculation accuracy, the variable when taking values needs to satisfy greater than the first preset value, which can be 100.

[0069] According to the target pore size interval, the first pore size interval, and the cumulative size distribution algorithm obtained above, the target pore fractal dimension is further determined.

[0070] In this embodiment, the target pore fractal dimension is calculated by determining the target pore size interval and the first pore size interval from the pore size distribution data, and a new cumulative size distribution algorithm.

[0071] In one embodiment, the pore fractal dimension determination method provided by the embodiments of the present application can further include:

[0072] Step 231, adding a fractal dimension assumption value in the adjustment term of the cumulative size distribution algorithm based on the target pore size interval and the first pore size interval;

[0073] Step 232, determining a fractal dimension response value corresponding to the fractal dimension assumption value through regression analysis;

[0074] Step 233, determining the target pore fractal dimension based on a two-dimensional variable relationship graph drawn based on the fractal dimension assumption value and the fractal dimension response value.

[0075] In step 6 above, a series of fractal dimension assumption values are added in the adjustment term of the general equation of the cumulative size distribution based on the target pore size interval, the number of pores in the target pore size interval, and the first pore size interval. The calculation method and principle of obtaining the fractal dimension response value corresponding to the fractal dimension assumption value through regression analysis are as follows:

[0076] According to the existing cumulative size distribution algorithm, the total number of pores with a pore size greater than in the entire cross section of the shale is related to the maximum pore size in the entire cross section and the fractal dimension There is a relationship as shown in Equation 1 below.

[0077] ; (1)

[0078] In Step 2 above, the microscopic image is only imaging part of the area of the shale section, and thus the number of pores obtained by shooting is obviously less than the number of pores in the entire section of the shale sample. In addition, due to the limitation of the field size of the microscopic imaging, the pores with a pore size greater than the field size cannot be directly shot, and the maximum pore size in the shooting area must be less than . Therefore, the total number of pores in the entire section of the shale sample with a pore size greater than or equal to can be represented as two parts as shown in Equation 2 below.

[0079] ; (2)

[0080] In Equation 2 above, is the number of pores in the entire section of the shale sample with a pore size greater than or equal to ; and is the number of pores in the entire section of the shale with a pore size greater than or equal to and less than . Wherein, is represented as shown in Equation 3 below.

[0081] ; (3)

[0082] Wherein, is the area of the entire section of the shale sample; is the area of the microscopic observation area; is the number of pores in the microscopic observation area with a pore size greater than or equal to and less than , and according to Equations 1, 2 and 3 above, Equation 4 below can be derived.

[0083] ; (4)

[0084] Equation 1 above is a special case of Equation 4 when , , and , and thus Equation 4 is a general equation of the traditional cumulative size distribution algorithm, which is applicable to randomly selecting part of the field for pore microscopic analysis, from having to perform pore microscopic analysis on the entire section of the shale sample. When a specific analysis process is used to study a specific shale sample, the , , , and​ are fixed values, so the above formula 4 can be further expressed as the following formula 5.

[0085] ; (5)

[0086] wherein, is a constant; is an adjustment term. is expressed as the following formula 6.

[0087] ; (6)

[0088] Fractal dimension The theoretical value range of the value is greater than 1 and less than 2. A series of fractal dimension hypothesis values are substituted for in the above formula 6, and are associated with the above formula 5, to obtain the following formula 7.

[0089] ; (7)

[0090] In the above formula 7, is a fractal dimension response value. Based on the pore quantity and its size distribution data obtained in step 3, regression analysis is performed with the independent variable and the dependent variable , and the reciprocal of the regression slope is the fractal dimension response value corresponding to each fractal dimension hypothesis value .

[0091] The embodiment determines the target pore fractal dimension by drawing a two-dimensional variable relationship graph based on the fractal dimension hypothesis value and the fractal dimension response value.

[0092] In one embodiment, the pore fractal dimension determination method provided by the embodiment of the application can further include:

[0093] Step 233-1, determining a target curve corresponding to the fractal dimension hypothesis value and the fractal dimension response value, and a target straight line where the fractal dimension hypothesis value is equal to the fractal dimension response value;

[0094] Step 233-2, determining a target value range corresponding to the fractal dimension hypothesis value;

[0095] Step 233-3, determining a target pore fractal dimension based on the target curve, the target straight line, and the target value range.

[0096] Specifically, the step 7 above draws a curve corresponding to the fractal dimension hypothesis value and the fractal dimension response value in the two-dimensional variable relationship diagram, confirms that the curve intersects with a straight line where the fractal dimension hypothesis value is equal to the fractal dimension response value, and there is only one intersection point in the range where the fractal dimension hypothesis value is greater than 1 and less than 2, and the fractal dimension hypothesis value corresponding to the intersection point is equal to the fractal dimension response value, that is, the actual calculation value of the fractal dimension.

[0097] The principle is as follows:

[0098] In the two-dimensional variable relationship diagram of the fractal dimension hypothesis value and the corresponding fractal dimension response value , the curve is drawn by using the fractal dimension hypothesis value and the fractal dimension response value collected in the step 6 above, and it is checked that the curve has only one intersection point with the straight line represented by in the range where If there is only one intersection point, according to the formula 5, the formula 6 and the formula 7 above, it is known that at this time, that is, the fractal dimension hypothesis value or the fractal dimension response value corresponding to the intersection point must be the calculated fractal dimension value , that is, the target pore fractal dimension in the embodiment.

[0099] The embodiment determines the target pore fractal dimension through the target curve, the target straight line and the target value range.

[0100] In one embodiment, the pore fractal dimension determination method provided by the embodiment of the present application can further include:

[0101] Step 310, determining an accumulated observation area based on a microscopic image of a random area on a shale surface;

[0102] Step 320, determining a target two-dimensional variable relationship diagram of the target pore fractal dimension and the accumulated observation area;

[0103] Step 330, performing a representative test on the target pore fractal dimension based on the accumulated observation area and the target two-dimensional variable relationship diagram;

[0104] Step 340, in a case where it is determined that the test passes based on the accumulated observation area and the target two-dimensional variable relationship diagram, determining that the target pore fractal dimension is a shale pore fractal dimension.

[0105] The pore fractal dimension determination method provided by the embodiment of the present application can further include:

[0106] Step 350, determining a fluctuation value of the target pore fractal dimension based on the target two-dimensional variable relationship diagram;

[0107] Step 360, determining the fluctuation value less than the target threshold value in the continuous interval as a target fluctuation value;

[0108] Step 370, in the case that the cumulative observation area corresponding to the target fluctuation value is greater than a first threshold value, determining that the representative test of the target pore fractal dimension is passed.

[0109] Specifically, in the step 10, the fluctuation value less than the target threshold value in the continuous interval is determined as the target fluctuation value, and in the case that the cumulative observation area corresponding to the target fluctuation value is not less than the first threshold value, it is determined that the representative test of the target pore fractal dimension is passed. In an available embodiment, the target threshold value can be 0.1, and the first threshold value can be .

[0110] The embodiment further improves the accuracy of the obtained pore fractal dimension by performing the representative test on the target pore fractal dimension.

[0111] The pore fractal dimension determination device provided by the present application is described below, and the pore fractal dimension determination device described below can be correspondingly referred to the pore fractal dimension determination method described above.

[0112] Please refer to Figure 3 , the present application also provides a pore fractal dimension determination device, comprising:

[0113] The pore size distribution data determination module 301 is configured to determine pore size distribution data based on a microscopic image of a random area on a shale surface.

[0114] The target pore fractal dimension determination module 302 is configured to determine a target pore fractal dimension based on the pore size distribution data and a cumulative size distribution algorithm.

[0115] The representative test module 303 is configured to perform a representative test on the target pore fractal dimension, and in the case that the test is passed, determine that the target pore fractal dimension is a shale pore fractal dimension.

[0116] Optionally, the target pore fractal dimension determination module comprises:

[0117] The first determination unit is configured to determine a target variable and an image maximum pore diameter based on the pore size distribution data.

[0118] The first pore diameter interval determination unit is configured to determine a target pore diameter interval based on the target variable and the image maximum pore diameter, and determine a first pore diameter interval based on the target variable.

[0119] A target pore fractal dimension determination unit is configured to determine a target pore fractal dimension according to the target pore size interval, the first pore size interval, and a cumulative size distribution algorithm.

[0120] Optionally, the target pore fractal dimension determination unit comprises:

[0121] A fractal dimension hypothesis value adding unit is configured to add a fractal dimension hypothesis value in an adjustment term of the cumulative size distribution algorithm based on the target pore size interval and the first pore size interval.

[0122] A fractal dimension response value determination unit is configured to determine a fractal dimension response value corresponding to the fractal dimension hypothesis value through regression analysis.

[0123] A second determination unit is configured to determine a target pore fractal dimension based on a two-dimensional variable relationship graph drawn based on the fractal dimension hypothesis value and the fractal dimension response value.

[0124] Optionally, the second determination unit comprises:

[0125] A third determination unit is configured to determine a target curve corresponding to the fractal dimension hypothesis value and the fractal dimension response value, and a target straight line where the fractal dimension hypothesis value is equal to the fractal dimension response value.

[0126] A target value range determination unit is configured to determine a target value range corresponding to the fractal dimension hypothesis value.

[0127] A fourth determination unit is configured to determine a target pore fractal dimension based on the target curve, the target straight line, and the target value range.

[0128] Optionally, the representative test on the target pore fractal dimension, in a case where the test passes, determining that the target pore fractal dimension is a shale pore fractal dimension comprises:

[0129] A cumulative observation area determination unit is configured to determine a cumulative observation area based on a microscopic image of a random area of a shale surface.

[0130] A target two-dimensional variable relationship graph determination unit is configured to determine a target two-dimensional variable relationship graph of the target pore fractal dimension and the cumulative observation area.

[0131] A representative test unit is configured to perform a representative test on the target pore fractal dimension based on the cumulative observation area and the target two-dimensional variable relationship graph.

[0132] A shale pore fractal dimension determination unit is configured to determine that the target pore fractal dimension is a shale pore fractal dimension in a case where it is determined that the test passes based on the cumulative observation area and the target two-dimensional variable relationship graph.

[0133] Optionally, the determining the target pore fractal dimension and the target two-dimensional variable relationship graph of the cumulative observation area includes:

[0134] A fluctuation value determination module is configured to determine a fluctuation value of the target pore fractal dimension based on the target two-dimensional variable relationship graph.

[0135] A target fluctuation value determination module is configured to determine the fluctuation value less than a target threshold value in a continuous interval as a target fluctuation value.

[0136] A representative test passing determination module is configured to determine that the representative test of the target pore fractal dimension is passed in a case where the cumulative observation area corresponding to the target fluctuation value is greater than a first threshold value.

[0137] Figure 4 An example of an entity structure diagram of an electronic device is shown in FIG. 1. Figure 4 As shown in FIG. 1, the electronic device can include a processor 410, a communications interface 420, a memory 430, and a communications bus 440, wherein the processor 410, the communications interface 420, and the memory 430 can complete mutual communication through the communications bus 440. The processor 410 can invoke a logical instruction in the memory 430 to execute a pore fractal dimension determination method, which includes: determining pore size distribution data based on a microscopic image of a random area of a shale surface; determining a target pore fractal dimension based on the pore size distribution data and a cumulative size distribution algorithm; and performing a representative test on the target pore fractal dimension, and determining the target pore fractal dimension as a shale pore fractal dimension in a case where the test is passed.

[0138] In addition, the logical instruction in the memory 430 described above can be implemented in the form of a software function unit and sold or used as an independent product, which can be stored in a computer-readable storage medium. Based on such understanding, the technical solutions of the present application essentially or the part that contributes to the prior art or part of the technical solutions can be embodied in the form of a software product, and the computer software product is stored in a storage medium, including a plurality of instructions to make a computer device (which can be a personal computer, a server, or a network device, etc.) execute all or part of the steps of the methods described in the embodiments of the present application. The foregoing storage medium includes: a U disk, a mobile hard disk, a read-only memory (ROM, Read-Only Memory), a random access memory (RAM, Random Access Memory), a magnetic disk or an optical disk, and various media that can store program codes.

[0139] In another aspect, the present application also provides a computer program product comprising a computer program, which can be stored on a non-transitory computer readable storage medium, and the computer program is executable by a processor to enable the computer to perform the method for determining the pore fractal dimension provided by any of the above methods, which comprises: determining pore size distribution data based on microscopic images of random areas of a shale surface; determining a target pore fractal dimension based on the pore size distribution data and a cumulative size distribution algorithm; and performing a representativeness test on the target pore fractal dimension, and in the case that the test is passed, determining the target pore fractal dimension as the pore fractal dimension of the shale.

[0140] In another aspect, the present application also provides a non-transitory computer readable storage medium having a computer program stored thereon, which is executable by a processor to implement the method for determining the pore fractal dimension provided by any of the above methods, which comprises: determining pore size distribution data based on microscopic images of random areas of a shale surface; determining a target pore fractal dimension based on the pore size distribution data and a cumulative size distribution algorithm; and performing a representativeness test on the target pore fractal dimension, and in the case that the test is passed, determining the target pore fractal dimension as the pore fractal dimension of the shale.

[0141] The apparatus embodiments described above are merely illustrative, wherein the units shown as separate components can or can not be physically separate, and the components shown as units can or can not be physical units, i.e., can be located in one place, or can be distributed on multiple network units. Part or all of the modules can be selected according to actual needs to achieve the purpose of the present embodiment. Those skilled in the art can understand and implement without creative labor.

[0142] From the above description of the embodiments, those skilled in the art can clearly understand that the embodiments can be implemented by means of software plus necessary general hardware platforms, and of course can also be implemented by hardware. Based on such understanding, the above technical solutions, essentially or in other words, the part that contributes to the prior art, can be embodied in the form of a software product, which can be stored in a computer readable storage medium, such as a ROM / RAM, a magnetic disk, an optical disk, etc., and includes a number of instructions to make a computer device (which can be a personal computer, a server, or a network device, etc.) execute the methods described in each embodiment or some parts of the embodiments.

[0143] It should be pointed out finally that the above embodiments are only used to illustrate the technical solutions of the present application, but not to limit the same; and although the present application has been described in detail with reference to the foregoing embodiments, it should be appreciated by those skilled in the art that the technical solutions recorded in the foregoing embodiments can be modified, or some technical features thereof can be replaced equivalently; and these modifications or replacements do not make the essence of the corresponding technical solutions deviate from the spirit and scope of the technical solutions of the embodiments of the present application.

Claims

1. A method for determining the fractal dimension of pores, characterized in that, include: Based on microscopic images of random regions on the shale surface, determine the pore size distribution data; Based on the pore size distribution data and the cumulative size distribution algorithm, the target pore fractal dimension is determined; The target pore fractal dimension is subjected to a representativeness test. If the test passes, the target pore fractal dimension is determined to be the shale pore fractal dimension. The determination of the target pore fractal dimension based on the pore size distribution data and the cumulative size distribution algorithm includes: The target variable and the maximum aperture of the image are determined based on the pore size distribution data; The target aperture range is determined based on the target variable and the maximum aperture of the image, and the first aperture range is determined based on the target variable; The target pore fractal dimension is determined based on the target pore size range, the first pore size range, and the cumulative size distribution algorithm. The step of determining the target pore fractal dimension based on the target pore size range, the first pore size range, and the cumulative size distribution algorithm includes: Based on the target aperture range and the first aperture range, a fractal dimension assumption value is added to the adjustment term of the cumulative size distribution algorithm; The fractal dimension response value corresponding to the assumed fractal dimension value is determined by regression analysis. The target pore fractal dimension is determined by plotting a two-dimensional variable relationship diagram based on the fractal dimension assumption and the fractal dimension response value. The process of performing a representativeness test on the target pore fractal dimension, and determining the target pore fractal dimension as the shale pore fractal dimension if the test passes, includes: The cumulative observed area was determined based on microscopic images of random regions on the shale surface. Determine the relationship between the target pore fractal dimension and the cumulative observed area in a two-dimensional variable diagram. Based on the relationship between the cumulative observed area and the target two-dimensional variable, the representativeness test of the target pore fractal dimension is performed. If the test is passed based on the relationship between the cumulative observed area and the target two-dimensional variable, the target pore fractal dimension is determined to be the shale pore fractal dimension. The process of determining the target two-dimensional variable relationship between the target pore fractal dimension and the cumulative observed area then includes: The fluctuation value of the target pore fractal dimension is determined based on the target two-dimensional variable relationship diagram; The fluctuation value that is less than the target threshold within a continuous interval is determined to be the target fluctuation value; If the cumulative observed area corresponding to the target fluctuation value is greater than the first threshold, the representativeness test of the target pore fractal dimension is passed.

2. The method for determining the fractal dimension of pores according to claim 1, characterized in that, The determination of the target pore fractal dimension by drawing a two-dimensional variable relationship diagram based on the fractal dimension assumption and the fractal dimension response value includes: Determine the target curves corresponding to the fractal dimension assumption value and the fractal dimension response value, and the target straight line where the fractal dimension assumption value is equal to the fractal dimension response value; Determine the target range of values ​​corresponding to the assumed fractal dimension; The target pore fractal dimension is determined based on the target curve, the target straight line, and the target value range.

3. A device for determining the fractal dimension of pores, characterized in that, include: The pore size distribution data determination module is used to determine pore size distribution data based on microscopic images of random areas on the shale surface; The target pore fractal dimension determination module is used to determine the target pore fractal dimension based on the pore size distribution data and the cumulative size distribution algorithm. The representativeness test module is used to perform a representativeness test on the target pore fractal dimension. If the test passes, the target pore fractal dimension is determined to be the shale pore fractal dimension. The determination of the target pore fractal dimension based on the pore size distribution data and the cumulative size distribution algorithm includes: The target variable and the maximum aperture of the image are determined based on the pore size distribution data; The target aperture range is determined based on the target variable and the maximum aperture of the image, and the first aperture range is determined based on the target variable; The target pore fractal dimension is determined based on the target pore size range, the first pore size range, and the cumulative size distribution algorithm. The step of determining the target pore fractal dimension based on the target pore size range, the first pore size range, and the cumulative size distribution algorithm includes: Based on the target aperture range and the first aperture range, a fractal dimension assumption value is added to the adjustment term of the cumulative size distribution algorithm; The fractal dimension response value corresponding to the assumed fractal dimension value is determined by regression analysis. The target pore fractal dimension is determined by plotting a two-dimensional variable relationship diagram based on the fractal dimension assumption and the fractal dimension response value. The process of performing a representativeness test on the target pore fractal dimension, and determining the target pore fractal dimension as the shale pore fractal dimension if the test passes, includes: The cumulative observed area was determined based on microscopic images of random regions on the shale surface. Determine the relationship between the target pore fractal dimension and the cumulative observed area in a two-dimensional variable diagram. Based on the relationship between the cumulative observed area and the target two-dimensional variable, the representativeness test of the target pore fractal dimension is performed. If the test is passed based on the relationship between the cumulative observed area and the target two-dimensional variable, the target pore fractal dimension is determined to be the shale pore fractal dimension. The process of determining the target two-dimensional variable relationship between the target pore fractal dimension and the cumulative observed area then includes: The fluctuation value of the target pore fractal dimension is determined based on the target two-dimensional variable relationship diagram; The fluctuation value that is less than the target threshold within a continuous interval is determined to be the target fluctuation value; If the cumulative observed area corresponding to the target fluctuation value is greater than the first threshold, the representativeness test of the target pore fractal dimension is passed.

4. An electronic device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the processor executes the program, it implements the method for determining the pore fractal dimension as described in any one of claims 1 to 2.

5. A non-transitory computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by the processor, it implements the method for determining the pore fractal dimension as described in any one of claims 1 to 2.

6. A computer program product, comprising a computer program, characterized in that, When the computer program is executed by the processor, it implements the method for determining the pore fractal dimension as described in any one of claims 1 to 2.

Images