A coal rock pore network calibration method, device and equipment

Abstract

The specification provides a coal rock pore network calibration method, device and equipment. The method comprises: obtaining an initial fracture system of a target coal rock, an initial matrix pore of the target coal rock, a preset main road equivalent rule and a preset secondary road equivalent rule; according to the preset main road equivalent rule, the initial fracture system is equivalent to a main road, and according to the preset secondary road equivalent rule, the initial matrix pore is equivalent to a secondary road; extracting micron-scale characteristic parameters in the main road and nanometer-scale characteristic parameters in the secondary road; according to the micron-scale characteristic parameters and the nanometer-scale characteristic parameters, a target coal rock pore network is generated; and the target coal rock pore network is calibrated as a whole to obtain a calibration result. Based on the above method, effective calibration of the target coal rock pore network can be realized.

Description

Technical Field

[0001] This specification relates to the field of coalbed methane development technology, and in particular to a method, apparatus and equipment for calibrating coal and rock pore networks. Background Technology

[0002] The coal pore network includes matrix pores and fracture systems. Reconstructing and calibrating the coal pore network is of great significance for the exploration, development, and exploitability evaluation of coalbed methane.

[0003] However, based on existing methods, it is impossible to obtain the true coal and rock pore network and to effectively calibrate the coal and rock pore network, thus making it impossible to accurately determine the effective flow space of the target coal and rock, thereby reducing the efficiency of coalbed methane development.

[0004] There is currently no effective solution to the above problems. Summary of the Invention

[0005] This specification provides a method, apparatus, and equipment for calibrating coal and rock pore networks, which can achieve effective calibration of coal and rock pore networks.

[0006] On the one hand, the embodiments of this specification provide a method for calibrating coal and rock pore networks, including:

[0007] Obtain the initial fracture system, initial matrix porosity, preset main road equivalent rules, and preset secondary road equivalent rules of the target coal and rock;

[0008] According to the preset main road equivalence rule, the initial fracture system is equivalent to a main road, and according to the preset secondary road equivalence rule, the initial matrix pores are equivalent to secondary roads.

[0009] Extract the micron-scale feature parameters of the main road and the nano-scale feature parameters of the secondary road;

[0010] Based on the micron-scale characteristic parameters and the nano-scale characteristic parameters, a target coal and rock pore network is generated;

[0011] The target coal and rock pore network was calibrated as a whole, and the calibration results were obtained.

[0012] Furthermore, the micrometer-scale characteristic parameters include at least one of the following: fracture system porosity, face cleavage spacing, face cleavage extension length, face cleavage aperture, face cleavage non-penetration probability, face cleavage roughness, end cleavage spacing, end cleavage aperture, end cleavage non-penetration probability, and end cleavage roughness; the nanometer-scale characteristic parameters include at least one of the following: matrix porosity and anisotropy characterization parameters.

[0013] Further, generating the target coal and rock pore network based on the micrometer-scale characteristic parameters and the nanometer-scale characteristic parameters includes:

[0014] Based on the micron-scale characteristic parameters, the growth space of the facet cleavage and the distribution characteristics of the facet cleavage in the growth space are determined;

[0015] The end-cut growth space is determined based on the end-cut spacing and the end-cut opening.

[0016] Based on the non-penetration probability of the end cut and the roughness of the end cut, the distribution characteristics of the end cut in the end cut growth space are determined;

[0017] Based on the distribution characteristics of the face cleavage growth space and the distribution characteristics of the end cleavage growth space, a reconstructed coal and rock fracture system is generated;

[0018] Based on the reconstructed coal and rock fracture system and the nanoscale characteristic parameters, a target coal and rock pore network is generated.

[0019] Further, the step of generating the target coal pore network based on the reconstructed coal fracture system and the nanoscale characteristic parameters includes:

[0020] Determine the porosity of the reconstructed coal-rock fracture system;

[0021] The porosity of the reconstructed coal-rock fracture system is compared with the porosity of the fracture system by performing a difference calculation.

[0022] When the result of the difference processing is less than a preset threshold, the reconstructed coal and rock fracture system is taken as the target coal and rock fracture system;

[0023] Based on the target coal and rock fracture system and the nanoscale characteristic parameters, a target coal and rock pore network is generated.

[0024] Further, generating the target coal pore network based on the target coal fracture system and the nanoscale characteristic parameters includes:

[0025] The target coal and rock fracture system is divided into blocks to obtain a generated space; wherein, the generated space includes multiple spatial coordinate points;

[0026] Determine the state of each spatial coordinate point among the plurality of spatial coordinate points, so as to obtain the reconstructed coal and rock matrix porosity based on the state;

[0027] Determine the porosity of the reconstructed coal and rock matrix pores;

[0028] The porosity of the reconstructed coal matrix porosity and the matrix porosity in the nanoscale characteristic parameters are subjected to difference processing. When the result of the difference processing is less than a preset threshold, the reconstructed coal matrix porosity is taken as the target coal matrix porosity.

[0029] Based on the target coal and rock fracture system and the target coal and rock matrix pores, a target coal and rock pore network is generated.

[0030] Further, determining the state of each spatial coordinate point among the plurality of spatial coordinate points includes:

[0031] Acquire nanoscale binarized scanning images;

[0032] By traversing the nanoscale binary scan image, multi-point neighborhood templates corresponding to the anisotropic characterization parameters are obtained;

[0033] Based on the multi-point neighborhood template, the state of each spatial coordinate point among the multiple spatial coordinate points is determined one by one.

[0034] Furthermore, the overall calibration of the target coal and rock pore network to obtain calibration results includes:

[0035] The target coal and rock pore network is calibrated as a whole according to the following formula:

[0036]

[0037] Among them, D CB Let r be the box dimension of the target coal and rock. k Let be the radius of the circle. For a radius of r k The minimum number of circles (spheres) required to cover the target coal and rock pore network;

[0038] The box dimension of the target coal and rock is used as the calibration result.

[0039] On the other hand, embodiments of this specification also provide a calibration device for coal and rock pore networks, comprising:

[0040] The acquisition module is used to acquire the initial fracture system of the target coal and rock, the initial matrix porosity of the target coal and rock, the preset main road equivalent rules, and the preset secondary road equivalent rules;

[0041] An equivalent module is used to convert the initial fracture system into a main road according to the preset main road equivalence rule, and to convert the initial matrix pores into a secondary road according to the preset secondary road equivalence rule.

[0042] The extraction module is used to extract micron-scale feature parameters from the main road and nanon-scale feature parameters from the secondary road;

[0043] The generation module is used to generate a target coal and rock pore network based on the micron-scale characteristic parameters and the nano-scale characteristic parameters.

[0044] The calibration module is used to perform overall calibration of the target coal and rock pore network and obtain calibration results.

[0045] In another aspect, this application also provides a calibration device for coal and rock pore networks, including a processor and a memory for storing processor-executable instructions, wherein the processor executes the instructions to implement the calibration method for coal and rock pore networks in the above embodiments.

[0046] In another aspect, this application also provides a computer-readable storage medium storing computer instructions thereon, wherein the computer-readable storage medium executes the instructions to implement the calibration method for coal and rock pore networks in the above embodiments.

[0047] This specification provides a method, apparatus, and equipment for calibrating a coal and rock pore network. First, it obtains the initial fracture system, initial matrix porosity, preset primary path equivalence rules, and preset secondary path equivalence rules of the target coal and rock. Second, it equivalences the initial fracture system to primary paths according to the preset primary path equivalence rules, and the initial matrix porosity to secondary paths according to the preset secondary path equivalence rules. Further, it extracts micron-scale feature parameters from the primary paths and nanon-scale feature parameters from the secondary paths. Finally, it generates the target coal and rock pore network based on the micron-scale and nanon-scale feature parameters. The target coal and rock pore network is then calibrated as a whole to obtain the calibration result. By equivalencing the initial fracture system to primary paths according to the preset primary path equivalence rules and the initial matrix porosity to secondary paths according to the preset secondary path equivalence rules, a realistic target coal and rock pore network can be obtained, achieving effective calibration of the coal and rock pore network and laying the foundation for subsequent flow simulation. Attached Figure Description

[0048] To more clearly illustrate the embodiments of this specification, the accompanying drawings used in the embodiments will be briefly introduced below. The drawings described below are only some embodiments recorded in this specification. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0049] Figure 1 This is a schematic flowchart of a method for calibrating a coal and rock pore network according to an embodiment of the present invention;

[0050] Figure 2 These are real coal samples from the embodiments of this invention;

[0051] Figure 3 This is the target coal and rock two-dimensional fracture system in the embodiments of the present invention;

[0052] Figure 4 This is the target coal and rock three-dimensional fracture system in the embodiments of the present invention;

[0053] Figure 5 This is the lateral anisotropy feature in the embodiments of the present invention;

[0054] Figure 6 This is the longitudinal anisotropy feature in the embodiments of the present invention;

[0055] Figure 7 This is the homogeneous connectivity feature in the embodiments of the present invention;

[0056] Figure 8 This is a two-dimensional two-point neighborhood template in the embodiments of the present invention;

[0057] Figure 9 This is the target coal and rock two-dimensional pore network in the embodiments of the present invention;

[0058] Figure 10 This is the three-dimensional pore network of the target coal and rock in the embodiments of the present invention;

[0059] Figure 11 This is a schematic diagram of the structural composition of the calibration device for coal and rock pore networks in an embodiment of the present invention;

[0060] Figure 12 This is a schematic diagram of the structural composition of a computer device provided in an embodiment of the present invention. Detailed Implementation

[0061] To enable those skilled in the art to better understand the technical solutions in this specification, the technical solutions in the embodiments of this specification will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of this specification, and not all embodiments. Based on the embodiments in this specification, all other embodiments obtained by those skilled in the art without creative effort should fall within the scope of protection of this specification.

[0062] Considering that the coal pore network includes matrix pores and fracture systems, obtaining and calibrating the coal pore network is of great significance for the exploration, development, and exploitability evaluation of coalbed methane. Current methods cannot obtain a true coal pore network and cannot effectively calibrate it, thus failing to accurately determine the effective flow space of the target coal and reducing the efficiency of coalbed methane development.

[0063] To address the aforementioned issues, this application proposes a method, apparatus, and equipment for calibrating coal and rock pore networks, thereby achieving effective calibration of coal and rock pore networks and improving coalbed methane development efficiency.

[0064] Based on the above approach, this specification proposes a calibration method for coal and rock pore networks. First, it obtains the initial fracture system, initial matrix porosity, preset main pathway equivalence rules, and preset secondary pathway equivalence rules of the target coal and rock. Second, according to the preset main pathway equivalence rules, the initial fracture system is equivalent to a main pathway, and according to the preset secondary pathway equivalence rules, the initial matrix porosity is equivalent to a secondary pathway. Further, it extracts micron-scale feature parameters from the main pathways and nano-scale feature parameters from the secondary pathways. Based on the micron-scale and nano-scale feature parameters, a target coal and rock pore network is generated. Finally, the target coal and rock pore network is calibrated as a whole to obtain the calibration result. Based on the calibration result, the effective flow space of the target coal and rock is determined. (See reference...) Figure 1 As shown in the embodiments of this specification, a method is provided. In specific implementation, the method may include the following.

[0065] S101: Obtain the initial fracture system of the target coal and rock, the initial matrix porosity of the target coal and rock, the preset main road equivalent rules, and the preset secondary road equivalent rules.

[0066] S102: The initial fracture system is equivalent to a main road according to the preset main road equivalence rule, and the initial matrix pores are equivalent to secondary roads according to the preset secondary road equivalence rule.

[0067] S103: Extract the micrometer-scale feature parameters of the main road and the nanometer-scale feature parameters of the secondary road.

[0068] S104: Generate the target coal and rock pore network based on the micron-scale characteristic parameters and the nano-scale characteristic parameters.

[0069] S105: Perform overall calibration of the target coal and rock pore network to obtain calibration results.

[0070] In some embodiments, the initial fracture system can be the coal-rock fracture system to be reconstructed, and the initial matrix porosity can be the coal-rock matrix porosity to be reconstructed. The coal-rock fracture system and the coal-rock matrix porosity constitute a coal-rock pore network. The preset main road equivalence rule can be used to equate the initial fracture system to a main road, and the preset secondary road equivalence rule can be used to equate the initial matrix porosity to a secondary road. The main road can be analogous to the main road in a highway network, and the secondary road can be analogous to the secondary road in a highway network. By equating the initial fracture system to a main road and the initial matrix porosity to a secondary road, the actual characteristics of coal and rock can be better reflected, laying the foundation for the subsequent accurate reconstruction of the coal-rock fracture system.

[0071] In some embodiments, before obtaining the initial fracture system of the target coal and rock, the initial matrix porosity of the target coal and rock, the preset main road equivalent rules, and the preset secondary road equivalent rules, the specific implementation may further include:

[0072] S1: Acquire micron-scale scanning images and nano-scale scanning images of the target coal and rock;

[0073] S2: Binarize the micron-scale scanning image and the nano-scale scanning image of the target coal and rock respectively to obtain the micron-scale binarized scanning image and the nano-scale binarized scanning image.

[0074] S3: Extract feature parameters from the micrometer-scale binarized scan images and the nanometer-scale binarized scan images respectively to obtain micrometer-scale feature parameters and nanometer-scale feature parameters;

[0075] S4: Generate the initial fracture system of the target coal and rock based on the micron-scale characteristic parameters, and generate the initial matrix pores of the target coal and rock based on the nano-scale characteristic parameters.

[0076] In some embodiments, micron-scale and nanoscale scan images of the target coal and rock can be acquired using different image scanning devices. An optimal threshold T for the aforementioned micron-scale and nanoscale scan images can be calculated using a threshold segmentation algorithm. Based on this optimal threshold, the micron-scale and nanoscale scan images are binarized. The binarization process can be expressed by the following formula:

[0077]

[0078] Among them, f s (x,y) represents the grayscale value of the scanned image, f(x,y) represents the grayscale value of the binarized image, and T represents the optimal threshold.

[0079] In some embodiments, feature parameters are extracted from the micrometer-scale binarized scan images and the nanometer-scale binarized scan images, respectively, to obtain micrometer-scale feature parameters and nanometer-scale feature parameters. The micrometer-scale feature parameters may include micrometer-scale features such as porosity, roughness, cleavage spacing, cleavage extension length, cleavage aperture, cleavage non-penetration probability, end-cleavage spacing, end-cleavage aperture, and end-cleavage non-penetration probability. The nanometer-scale feature parameters may include nanometer-scale matrix porosity and anisotropy characterization parameters.

[0080] In some embodiments, based on the extracted feature parameters of the binarized scan image, micron-scale feature parameters of the main road and nano-scale feature parameters of the secondary road can be further extracted. The micron-scale feature parameters of the main road include at least one of the following: porosity of the fracture system, cleavage spacing, cleavage extension length, cleavage aperture, cleavage non-penetration probability, cleavage roughness, end cleavage spacing, end cleavage aperture, end cleavage non-penetration probability, and end cleavage roughness. The nano-scale feature parameters of the secondary road include at least one of the following: matrix porosity and anisotropic characterization parameters.

[0081] In some embodiments, the generation of the target coal and rock pore network based on the micrometer-scale characteristic parameters and the nanometer-scale characteristic parameters may, in specific implementation, include:

[0082] S1: Based on the micrometer-scale characteristic parameters, determine the face cleavage growth space and the distribution characteristics of face cleavage in the face cleavage growth space;

[0083] S2: Determine the end cut growth space based on the end cut spacing and the end cut opening;

[0084] S3: Determine the distribution characteristics of the end cut in the end cut growth space based on the end cut non-penetration probability and the end cut roughness;

[0085] S4: Based on the distribution characteristics of the face cleavage growth space and the distribution characteristics of the end cleavage growth space, a reconstructed coal and rock fracture system is generated;

[0086] S5: Generate the target coal pore network based on the reconstructed coal fracture system and the nanoscale characteristic parameters.

[0087] In some embodiments, determining the facet growth space and the distribution characteristics of facets in the facet growth space based on the micrometer-scale characteristic parameters may, in specific implementations, include:

[0088] S1: Determine the surface cleavage growth space based on the surface cleavage spacing, the surface cleavage extension length, and the surface cleavage opening.

[0089] S2: Determine the distribution characteristics of the face cut in the face cut growth space based on the face cut non-penetration probability and the face cut roughness.

[0090] In some embodiments, when generating the growth space of surface cleavage, a certain random deviation can be added to the surface cleavage spacing, surface cleavage extension length, and surface cleavage aperture. When generating the growth space of end cleavage, a certain random deviation can be added to the end cleavage spacing and end cleavage aperture. By adding random deviations, the coal-rock fracture system can be reconstructed. It should be noted that each time the above reconstruction process is repeated, the obtained coal-rock fracture system is different.

[0091] In some embodiments, the generation of the target coal pore network based on the reconstructed coal fracture system and the nanoscale characteristic parameters may, in specific implementation, include:

[0092] S1: Determine the porosity of the reconstructed coal-rock fracture system;

[0093] S2: Perform difference processing on the porosity of the reconstructed coal-rock fracture system and the porosity of the fracture system;

[0094] S3: When the result of the difference processing is less than a preset threshold, the reconstructed coal and rock fracture system is taken as the target coal and rock fracture system;

[0095] S4: Generate a target coal and rock pore network based on the target coal and rock fracture system and the nanoscale characteristic parameters.

[0096] In some embodiments, the porosity of the reconstructed coal-rock fracture system is compared with the porosity of the fracture system to obtain an error value between the porosity of the reconstructed coal-rock fracture system and the porosity of the fracture system. A certain condition can be set; when the porosity of the reconstructed coal-rock fracture system meets this condition, the iteration stops, the final reconstructed coal-rock fracture system is output, and this final reconstructed coal-rock fracture system is used as the target coal-rock fracture system. The condition can be comparing the error value with a preset threshold. If the error value is less than the preset threshold, the reconstructed coal-rock fracture system is directly output; if it is greater than the preset threshold, the iteration continues to re-acquire the reconstructed coal-rock fracture system until the porosity of the re-acquired reconstructed coal-rock fracture system is less than the preset threshold, at which point the iteration stops. Of course, the conditions set are not limited to the examples above. Those skilled in the art may make other changes based on the technical essence of the embodiments in this specification, but as long as the achieved function and effect are the same as or similar to the embodiments in this specification, they should all be covered within the protection scope of the embodiments in this specification.

[0097] In some embodiments, the generation of the target coal pore network based on the target coal fracture system and the nanoscale characteristic parameters may, in specific implementation, include:

[0098] S1: The target coal and rock fracture system is divided into blocks to obtain a generated space; wherein, the generated space includes multiple spatial coordinate points;

[0099] S2: Determine the state of each spatial coordinate point among the plurality of spatial coordinate points, so as to obtain the reconstructed coal and rock matrix porosity based on the state;

[0100] S3: Determine the porosity of the reconstructed coal and rock matrix pores;

[0101] S4: Perform a difference processing on the porosity of the reconstructed coal matrix porosity and the matrix porosity in the nanoscale characteristic parameters. When the result of the difference processing is less than a preset threshold, the reconstructed coal matrix porosity is taken as the target coal matrix porosity.

[0102] S5: Generate a target coal and rock pore network based on the target coal and rock fracture system and the target coal and rock matrix pores.

[0103] In some embodiments, the porosity of the reconstructed coal matrix pores and the matrix pore porosity in the nanoscale characteristic parameters are differentially processed to obtain the error value between the porosity of the reconstructed coal matrix pores and the matrix pore porosity in the nanoscale characteristic parameters. A certain condition can be set, and when the porosity of the reconstructed coal matrix pores meets the condition, the iteration stops, the final reconstructed coal matrix pores are output, and the final reconstructed coal matrix pores are used as the target coal matrix pores. The above condition can be that the error value is compared with a preset threshold. If it is less than the preset threshold, the reconstructed coal matrix pores are directly output. If it is greater than the preset threshold, the iteration continues to re-acquire the reconstructed coal matrix pores until the porosity of the re-acquired reconstructed coal matrix pores is less than the preset threshold, at which point the iteration stops.

[0104] In some embodiments, the determination of the state of each spatial coordinate point among the plurality of spatial coordinate points may, in specific implementation, include:

[0105] S1: Acquire nanoscale binarized scanning images;

[0106] S2: Traverse the nanoscale binary scan image to obtain the multi-point neighborhood templates corresponding to the anisotropic characterization parameters;

[0107] S3: Based on the multi-point neighborhood template, determine the state of each spatial coordinate point among the multiple spatial coordinate points one by one.

[0108] In some embodiments, the target coal and rock fracture system is a matrix composed of 0 and 1. Therefore, determining the state of each spatial coordinate point, that is, determining whether the state of each spatial coordinate point is 0 or 1, if it is 0, it can be said that the spatial coordinate point is the matrix, and if it is 1, it can be said that the spatial coordinate point is the pore. By determining whether the spatial coordinate point is the pore or the matrix, the foundation can be laid for obtaining the reconstructed coal and rock matrix pores.

[0109] In some embodiments, after the overall calibration of the target coal and rock pore network is performed and the calibration results are obtained, the specific implementation may further include: determining the effective flow space of the target coal and rock based on the calibration results.

[0110] In some embodiments, the above-described overall calibration of the target coal and rock pore network to obtain calibration results may, in specific implementation, include:

[0111] S1: The target coal and rock pore network is calibrated as a whole according to the following formula:

[0112]

[0113] Among them, D CB Let r be the box dimension of the target coal and rock. k N is the radius of the circle (sphere). rk For a radius of r k The minimum number of circles (spheres) required to cover the target coal and rock pore network;

[0114] S2: Use the box dimension of the target coal and rock as the calibration result.

[0115] In some embodiments, during implementation, a set of radii r of circles (spheres) are set. k You can then obtain a corresponding set. ln1 / r k The x-axis is N. rk Using the vertical axis as the ordinate and representing it in a logarithmic coordinate system, the slope of the obtained curve is the box dimension of the target coal and rock.

[0116] In some embodiments, by obtaining the box dimension, the overall pore network of the target coal and rock can be effectively calibrated. Effective calibration can lay the foundation for subsequent accurate and rapid simulation of the effective flow space of the target coal and rock, thereby improving the efficiency of coalbed methane development.

[0117] In a specific scenario example, the calibration method for coal and rock pore networks provided in the embodiments of this specification can be applied to achieve effective calibration of coal and rock pore networks. Before calibration, a target coal and rock pore network needs to be generated. The target coal and rock pore network is composed of a target coal and rock fracture system and a target coal and rock matrix pore. The following describes how to reconstruct the target coal and rock fracture system and the target coal and rock matrix pore in a specific implementation scenario.

[0118] To obtain the target coal and rock fracture system, the following methods are mainly adopted:

[0119] First, the growth space of the face cleavage is determined based on the spacing, extension length, and opening of the face cleavage, with a certain random deviation applied. Then, the growth direction and actual extension length of the face cleavage are calculated based on the non-penetration probability, and the face cleavage distribution is generated by combining this with the face cleavage roughness. Next, the growth space of the end cleavage is determined based on the spacing and opening of the end cleavage, with a certain random deviation applied. Finally, the growth direction and extension length of the end cleavage are calculated based on the non-penetration probability, and the end cleavage distribution is generated by combining this with the end cleavage roughness.

[0120] Secondly, a reconstructed coal-rock fracture system is generated based on the distribution of surface cleavage and end cleavage, and the porosity of the reconstructed coal-rock fracture system is calculated. The porosity of the reconstructed coal-rock fracture system is then subtracted from the porosity of the fracture system in the micron-scale characteristic parameters of the main road to obtain the error value between the porosity of the reconstructed coal-rock fracture system and the porosity of the fracture system. It is then determined whether the error value between the porosity of the reconstructed coal-rock fracture system and the porosity of the fracture system is less than a preset threshold. If it is less than the preset threshold, the reconstructed coal-rock fracture system is used as the target coal-rock fracture system, i.e., the target coal-rock fracture system is output. If it is not less than the preset threshold, the process needs to be repeated to obtain the distribution of surface cleavage and end cleavage, and the reconstructed coal-rock fracture system needs to be regenerated until the error value between the porosity of the reconstructed coal-rock fracture system and the porosity of the fracture system is less than the preset threshold, thus obtaining the target coal-rock fracture system.

[0121] By acquiring the target coal and rock fracture system, the coal and rock pore network can be more accurately reflected. Figure 2 A real coal sample was shown. Figure 3 The target coal rock two-dimensional fracture system is shown (at this point, matrix pores have not yet been added).

[0122] It should be noted that, for reference Figure 4 As shown, Figure 4 The target coal and rock three-dimensional fracture system is shown. The target coal and rock fracture system obtained by the method in the above embodiments can be a two-dimensional fracture system or a three-dimensional fracture system. This specification does not make a specific limitation on this.

[0123] To obtain the target matrix porosity, the following methods are mainly adopted:

[0124] First, the target coal and rock fracture system obtained above is divided into blocks to obtain a generation space; the generation space contains multiple spatial coordinate points. For example, a 50×50 generation space contains spatial coordinate points (1,1), (1,2), ... (50,50). The block division process can be understood as dividing the target coal and rock fracture system into blocks and defining the blocks as the generation space of the target coal and rock matrix pores.

[0125] Secondly, by traversing the aforementioned nanoscale binarized scanning images, multi-point neighborhood templates corresponding to the anisotropic characterization parameters are obtained. These anisotropic characterization parameters, in the two-dimensional case, can be categorized into lateral, longitudinal, and homogeneous connectivity features. (See [reference needed]). Figures 5 to 7 As shown, see reference Figures 5 to 7 The transverse anisotropy, longitudinal anisotropy, and homogeneous connectivity characteristics are shown respectively. Anisotropy can be represented using neighborhood templates. Based on the acquired multi-point neighborhood templates, the state of each spatial coordinate point is determined one by one. Since the target coal and rock fracture system is a matrix composed of 0s and 1s, determining the state of each spatial coordinate point involves determining whether its state is 0 or 1. If it is 0, it indicates that the spatial coordinate point is the matrix; if it is 1, it indicates that the spatial coordinate point is a pore.

[0126] Further reading Figure 8 As shown, taking a two-dimensional neighborhood of two points as an example, when point (i, j) is 0 or 1, there are four possibilities for point (i, j+1) to be 0 or 1: when point (i, j) is 0, point (i, j+1) is 0 or 1; when point (i, j) is 1, point (i, j+1) is 0 or 1. By traversing all points in the generation space, the probability of each point is obtained. For example, the probability that point (i, j) is 0 and point (i, j+1) is 0 can be obtained. For instance, if the probability that point (1, 1) is 1 (pore) and point (1, 2) is 1 (pore) is 0.1, a random number is generated between 0 and 1. If the random number is less than 0.1, point (1, 2) is considered 1 (pore); otherwise, it is 0 (matrix). After obtaining these probabilities, the pores of the coal and rock matrix can be reconstructed.

[0127] Finally, the porosity of the reconstructed coal matrix is ​​calculated. The difference between the porosity of the reconstructed coal matrix and the porosity of the matrix in the nanoscale characteristic parameters is processed to obtain the error value between the porosity of the reconstructed coal matrix and the porosity of the matrix. It is then determined whether the error value between the porosity of the reconstructed coal matrix and the porosity of the matrix is ​​less than a preset threshold. If it is less than the preset threshold, the reconstructed coal matrix porosity is taken as the target coal matrix porosity. If it is not less than the preset threshold, the iteration is repeated. That is, the state of each point in the generation space is judged one by one according to the obtained multi-point neighborhood template, so as to re-obtain the reconstructed coal matrix porosity. The iteration stops when the re-obtained reconstructed coal matrix porosity is less than the preset threshold, and the target coal matrix porosity is finally output.

[0128] It should be noted that, for reference Figure 9 , Figure 10 As shown, Figure 9 The two-dimensional pore network of the target coal and rock is shown. Figure 10 The three-dimensional pore network of the target coal and rock is shown. The pore network of the target coal and rock obtained by the method in the above embodiments can be a two-dimensional pore network or a three-dimensional pore network. This specification does not make a specific limitation on this.

[0129] The above method will be described below with reference to a specific embodiment. However, it is worth noting that this specific embodiment is only for better illustration of this application and does not constitute an improper limitation of this application.

[0130] Before implementation, firstly, micron-scale and nano-scale scan images of the target coal and rock are acquired; secondly, the micron-scale and nano-scale scan images of the target coal and rock are binarized to obtain micron-scale binarized scan images and nano-scale binarized scan images; further, feature parameters are extracted from the micron-scale and nano-scale binarized scan images to obtain micron-scale feature parameters and nano-scale feature parameters; finally, based on the micron-scale feature parameters, the initial fracture system of the target coal and rock is generated, and based on the nano-scale feature parameters, the initial matrix pores of the target coal and rock are generated.

[0131] In practice, the initial fracture system of the target coal and rock, the initial matrix porosity of the target coal and rock, the preset main road equivalent rules, and the preset secondary road equivalent rules are first obtained.

[0132] Secondly, the initial fracture system is equivalent to a main road according to the preset main road equivalence rule, and the initial matrix pores are equivalent to secondary roads according to the preset secondary road equivalence rule; then, the micron-scale feature parameters in the main road and the nano-scale feature parameters in the secondary road are extracted.

[0133] Furthermore, based on the micrometer-scale characteristic parameters of face cleavage spacing, face cleavage extension length, and face cleavage aperture, and with the application of a certain random bias, the growth space of face cleavage is determined; based on the face cleavage non-penetration probability and face cleavage roughness, the distribution characteristics of face cleavage in the face cleavage growth space are determined; based on the nanometer-scale characteristic parameters of end cleavage spacing and end cleavage aperture, and with the application of a certain random bias, the growth space of end cleavage is determined; based on the end cleavage non-penetration probability and end cleavage roughness, the distribution characteristics of end cleavage in the end cleavage growth space are determined; based on the distribution characteristics of face cleavage in the face cleavage growth space and end cleavage in the end cleavage growth space, a reconstructed coal-rock fracture system is generated. The porosity of the reconstructed coal-rock fracture system is calculated and then compared with the porosity of the fracture system in the micrometer-scale characteristic parameters. If the error value between the reconstructed coal-rock fracture system porosity and the fracture system porosity in the micrometer-scale characteristic parameters is less than a preset threshold, the reconstructed coal-rock fracture system is output and used as the target coal-rock fracture system. If the error value is greater than or equal to the preset threshold, the reconstructed coal-rock fracture system is obtained again according to the above steps until the error value between the porosity of the reconstructed coal-rock fracture system and the fracture system porosity in the micrometer-scale characteristic parameters is less than the preset threshold. Then the iteration stops and the reconstructed coal-rock fracture system is output as the target coal-rock fracture system.

[0134] Furthermore, the acquired target coal and rock fracture system is divided into blocks to obtain a generation space; the generation space includes multiple spatial coordinate points; the acquired nanoscale binarized scan images are traversed to obtain multi-point neighborhood templates corresponding to anisotropic characterization parameters; based on the multi-point neighborhood templates, the state of each spatial coordinate point among the multiple spatial coordinate points is determined one by one, the state including 0 or 1 (matrix or pore); based on the determined state of each spatial coordinate point, the reconstructed coal and rock matrix pores are obtained and the porosity of the reconstructed coal and rock matrix pores is calculated; the porosity of the reconstructed coal and rock matrix pores is then determined. If the error value between the reconstructed coal matrix porosity and the matrix porosity in the nanoscale feature parameters is less than a preset threshold, then the reconstructed coal matrix porosity is output and used as the target coal matrix porosity. If the error value is greater than or equal to the preset threshold, then the reconstructed coal matrix porosity is obtained again according to the above steps until the error value between the porosity of the reconstructed coal matrix porosity and the matrix porosity in the nanoscale feature parameters is less than the preset threshold. Then the iteration stops and the reconstructed coal matrix porosity is output as the target coal matrix porosity.

[0135] Finally, based on the target coal and rock fracture system and the target coal and rock matrix pores, a target coal and rock pore network is generated. Then, the box dimension from the fractal dimension is used to calibrate the target coal and rock pore network as a whole, obtaining a specific value for the box dimension. After obtaining the specific value of the box dimension, the effective flow space of the target coal and rock can be determined, thereby effectively improving the efficiency of coalbed methane development.

[0136] Although this specification provides the following examples or appendices Figure 11 The methods, steps, or apparatus structures shown may include more or fewer combined operational steps or module units based on conventional or non-inventive methods. In steps or structures where there is no logically necessary causal relationship, the execution order of these steps or the module structure of the apparatus is not limited to the execution order or module structure shown in the embodiments or drawings of this specification. When the methods or module structures described are applied in actual devices, servers, or terminal products, they can be executed sequentially or in parallel according to the methods or module structures shown in the embodiments or drawings (e.g., in parallel processor or multi-threaded processing environments, or even distributed processing or server cluster implementation environments).

[0137] Based on the above-described method for calibrating coal and rock pore networks, this specification also provides an embodiment of a calibration device for coal and rock pore networks. For example... Figure 11 As shown, the device may specifically include the following modules:

[0138] The acquisition module 1101 is used to acquire the initial fracture system of the target coal and rock, the initial matrix porosity of the target coal and rock, the preset main road equivalent rules and the preset secondary road equivalent rules;

[0139] The equivalent module 1102 is used to convert the initial fracture system into a main road according to the preset main road equivalent rule, and to convert the initial matrix pores into a secondary road according to the preset secondary road equivalent rule.

[0140] Extraction module 1103 is used to extract micron-scale feature parameters in the main road and nano-scale feature parameters in the secondary road;

[0141] The generation module 1104 is used to generate a target coal and rock pore network based on the micron-scale characteristic parameters and the nano-scale characteristic parameters.

[0142] The calibration module 1105 is used to perform overall calibration of the target coal and rock pore network and obtain calibration results.

[0143] In some embodiments, the micrometer-scale characteristic parameters in the acquisition module 1101 include at least one of the following: fracture system porosity, face cleavage spacing, face cleavage extension length, face cleavage aperture, face cleavage non-penetration probability, face cleavage roughness, end cleavage spacing, end cleavage aperture, end cleavage non-penetration probability, and end cleavage roughness; the nanometer-scale characteristic parameters include at least one of the following: matrix porosity and anisotropy characterization parameters.

[0144] In some embodiments, the generation module 1104 can be specifically used to: determine the face cleavage growth space and the distribution characteristics of face cleavage in the face cleavage growth space based on the micrometer-scale characteristic parameters; determine the end cleavage growth space based on the end cleavage spacing and the end cleavage aperture; determine the distribution characteristics of end cleavage in the end cleavage growth space based on the end cleavage non-penetration probability and the end cleavage roughness; generate a reconstructed coal and rock fracture system based on the distribution characteristics of the face cleavage growth space and the distribution characteristics of the end cleavage growth space; and generate a target coal and rock pore network based on the reconstructed coal and rock fracture system and the nanometer-scale characteristic parameters.

[0145] In some embodiments, the generation module 1104 may further be used to determine the porosity of the reconstructed coal-rock fracture system; perform difference processing on the porosity of the reconstructed coal-rock fracture system and the porosity of the fracture system; when the result of the difference processing is less than a preset threshold, the reconstructed coal-rock fracture system is used as the target coal-rock fracture system; and generate a target coal-rock pore network based on the target coal-rock fracture system and the nanoscale characteristic parameters.

[0146] In some embodiments, the generation module 1104 can further be used to perform block division processing on the target coal and rock fracture system to obtain a generation space; wherein, the generation space includes multiple spatial coordinate points; determine the state of each spatial coordinate point among the multiple spatial coordinate points, so as to obtain the reconstructed coal and rock matrix porosity according to the state; determine the porosity of the reconstructed coal and rock matrix porosity; perform difference processing on the porosity of the reconstructed coal and rock matrix porosity and the matrix porosity in the nanoscale characteristic parameters, and when the result of the difference processing is less than a preset threshold, the reconstructed coal and rock matrix porosity is taken as the target coal and rock matrix porosity; generate a target coal and rock pore network according to the target coal and rock fracture system and the target coal and rock matrix porosity.

[0147] In some embodiments, the generation module 1104 may also be used to acquire a nanoscale binarized scan image; traverse the nanoscale binarized scan image to acquire a multi-point neighborhood template corresponding to the anisotropic characterization parameters; and determine the state of each spatial coordinate point among the plurality of spatial coordinate points according to the multi-point neighborhood template.

[0148] In some embodiments, the calibration module 1105 can be specifically used to perform overall calibration of the target coal and rock pore network according to the following formula:

[0149]

[0150] Among them, D CB Let r be the box dimension of the target coal and rock. k Let be the radius of the circle. For a radius of rk The minimum number of circles required to cover the pore network of the target coal and rock; the box dimension of the target coal and rock is used as the calibration result.

[0151] It should be noted that the units, devices, or modules described in the above embodiments can be implemented by computer chips or physical entities, or by products with certain functions. For ease of description, the above devices are described by dividing them into various modules according to their functions. Of course, in implementing this specification, the functions of each module can be implemented in one or more software and / or hardware, or the module that implements the same function can be implemented by a combination of multiple sub-modules or sub-units, etc. The device embodiments described above are merely illustrative. For example, the division of units is only a logical functional division, and there may be other division methods in actual implementation. For example, multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Furthermore, the coupling or direct coupling or communication connection between the devices or units shown or discussed can be through some interfaces, and the indirect coupling or communication connection between devices or units can be electrical, mechanical, or other forms.

[0152] As can be seen from the above, the calibration device for coal and rock pore networks provided in the embodiments of this specification can accurately reflect the geometric morphological characteristics of coal and rock pore networks, laying the foundation for subsequent flow simulation and improving the development efficiency of coalbed methane resources.

[0153] This specification also provides an electronic device for calibrating coal and rock pore networks, including a processor and a memory for storing processor-executable instructions. Specifically, the processor can perform the following steps according to the instructions: acquiring the initial fracture system of the target coal and rock, the initial matrix pores of the target coal and rock, a preset main road equivalence rule, and a preset secondary road equivalence rule; equipping the initial fracture system as a main road according to the preset main road equivalence rule, and equipping the initial matrix pores as secondary roads according to the preset secondary road equivalence rule; extracting micron-scale feature parameters from the main roads and nano-scale feature parameters from the secondary roads; generating the target coal and rock pore network based on the micron-scale feature parameters and the nano-scale feature parameters; and performing overall calibration of the target coal and rock pore network to obtain the calibration result.

[0154] To execute the above instructions more accurately, please refer to... Figure 12 As shown in the embodiments of this specification, another specific electronic device is also provided, wherein the electronic device includes a network communication port 1201, a processor 1202 and a memory 1203, and the above structures are connected by internal cables so that the various structures can perform specific data interaction.

[0155] Specifically, the network communication port 1201 can be used to acquire the initial fracture system of the target coal and rock, the initial matrix porosity of the target coal and rock, the preset main road equivalent rules, and the preset secondary road equivalent rules.

[0156] The processor 1202 can be specifically used to: equip the initial fracture system as a primary road according to the preset primary road equivalence rule; equip the initial matrix pores as secondary roads according to the preset secondary road equivalence rule; extract micron-scale feature parameters from the primary roads and nano-scale feature parameters from the secondary roads; generate a target coal and rock pore network based on the micron-scale feature parameters and the nano-scale feature parameters; and perform overall calibration on the target coal and rock pore network to obtain calibration results.

[0157] The memory 1203 can be used to store the corresponding instruction program.

[0158] In this embodiment, the network communication port 1201 can be a virtual port bound to different communication protocols, thereby enabling the sending or receiving of different data. For example, the network communication port can be a port responsible for web data communication, a port responsible for FTP data communication, or a port responsible for email data communication. Furthermore, the network communication port can also be a physical communication interface or communication chip. For example, it can be a wireless mobile network communication chip, such as GSM or CDMA; it can also be a Wi-Fi chip; or it can be a Bluetooth chip.

[0159] In this embodiment, the processor 1202 can be implemented in any suitable manner. For example, the processor can take the form of a microprocessor or processor and a computer-readable medium storing computer-readable program code (e.g., software or firmware) executable by the (micro)processor, logic gates, switches, application-specific integrated circuits (ASICs), programmable logic controllers, and embedded microcontrollers, etc. This specification is not limiting.

[0160] In this embodiment, the memory 1203 may include multiple layers. In a digital system, anything that can store binary data can be a memory. In an integrated circuit, a circuit with storage function but no physical form is also called a memory, such as RAM, FIFO, etc. In a system, a storage device with a physical form is also called a memory, such as a memory stick, TF card, etc.

[0161] This specification also provides a computer storage medium based on the above method, wherein the computer storage medium stores computer program instructions, which, when executed, perform the following: acquiring the initial fracture system of the target coal and rock, the initial matrix porosity of the target coal and rock, a preset main road equivalence rule, and a preset secondary road equivalence rule; equipping the initial fracture system as a main road according to the preset main road equivalence rule, and equipping the initial matrix porosity as a secondary road according to the preset secondary road equivalence rule; extracting micron-scale feature parameters from the main road and nanon-scale feature parameters from the secondary road; generating a target coal and rock pore network based on the micron-scale feature parameters and the nanon-scale feature parameters; and performing overall calibration on the target coal and rock pore network to obtain calibration results.

[0162] In this embodiment, the storage medium includes, but is not limited to, Random Access Memory (RAM), Read-Only Memory (ROM), cache, hard disk drive (HDD), or memory card. The memory can be used to store computer program instructions. The network communication unit can be an interface configured according to standards specified in the communication protocol for network connection communication.

[0163] While this specification provides the steps of operation for the methods described in the embodiments or flowcharts, more or fewer steps may be included based on conventional or non-inventive means. The order of steps listed in the embodiments is merely one possible order of execution among many steps and does not represent the only possible order. In actual device or client product execution, the methods shown in the embodiments or drawings may be executed sequentially or in parallel (e.g., in a parallel processor or multi-threaded processing environment, or even a distributed data processing environment). The terms "comprising," "including," or any other variations thereof are intended to cover a non-exclusive inclusion, such that a process, method, product, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, product, or apparatus. Without further limitations, the presence of other identical or equivalent elements in a process, method, product, or apparatus that includes said elements is not excluded. The terms "first," "second," etc., are used to denote names and do not indicate any particular order.

[0164] Those skilled in the art will also know that, besides implementing the controller using purely computer-readable program code, the same functions can be achieved by logically programming the method steps, making the controller function as logic gates, switches, application-specific integrated circuits (ASICs), programmable logic controllers (PLCs), and embedded microcontrollers. Therefore, such a controller can be considered a hardware component, and the devices within it used to implement various functions can also be considered structures within that hardware component. Alternatively, the devices used to implement various functions can be considered as both software modules implementing the method and structures within a hardware component.

[0165] This specification can be described in the general context of computer-executable instructions that are executed by a computer, such as program modules. Generally, program modules include routines, programs, objects, components, data structures, classes, etc., that perform a specific task or implement a specific abstract data type. This specification can also be practiced in distributed computing environments, where tasks are performed by remote processing devices connected via a communication network. In distributed computing environments, program modules can reside in local and remote computer storage media, including storage devices.

[0166] As can be seen from the above description of the embodiments, those skilled in the art can clearly understand that this specification can be implemented by means of software plus necessary general-purpose hardware platforms. Based on this understanding, the technical solutions of this specification can essentially be embodied in the form of a software product. This computer software product can be stored in a storage medium, such as ROM / RAM, magnetic disk, optical disk, etc., and includes several instructions to cause a computer device (which may be a personal computer, mobile terminal, server, or network device, etc.) to execute the methods described in the various embodiments or some parts of the embodiments of this specification.

[0167] The various embodiments in this specification are described in a progressive manner. Similar or identical parts between embodiments can be referred to interchangeably. Each embodiment focuses on its differences from other embodiments. This specification can be used in numerous general-purpose or special-purpose computer system environments or configurations. Examples include: personal computers, server computers, handheld or portable devices, tablet devices, multiprocessor systems, microprocessor-based systems, set-top boxes, programmable electronic devices, network PCs, minicomputers, mainframe computers, and distributed computing environments including any of the above systems or devices, etc.

[0168] Although this specification has been described by way of examples, those skilled in the art will recognize that many variations of this specification are possible without departing from its spirit, and it is intended that the appended claims cover such variations without departing from the spirit of this specification.

Claims

1. A method for calibrating coal and rock pore networks, characterized in that, include: Obtain the initial fracture system, initial matrix porosity, preset main road equivalent rules, and preset secondary road equivalent rules of the target coal and rock; According to the preset main road equivalence rule, the initial fracture system is equivalent to a main road, and according to the preset secondary road equivalence rule, the initial matrix pores are equivalent to secondary roads. Extract the micron-scale feature parameters of the main road and the nano-scale feature parameters of the secondary road; Based on the micron-scale characteristic parameters and the nano-scale characteristic parameters, a target coal and rock pore network is generated; The target coal and rock pore network was calibrated as a whole, and the calibration results were obtained. The step of generating the target coal and rock pore network based on the micrometer-scale characteristic parameters and the nanometer-scale characteristic parameters includes: Based on the micron-scale characteristic parameters, the growth space of the facet cleavage and the distribution characteristics of the facet cleavage in the growth space are determined; The end cleavage growth space is determined based on the end cleavage spacing and end cleavage opening in the micrometer-scale characteristic parameters. Based on the end cut non-penetration probability and end cut roughness in the micrometer-scale characteristic parameters, the distribution characteristics of the end cut in the end cut growth space are determined; Based on the distribution characteristics of the face cleavage growth space and the distribution characteristics of the end cleavage growth space, a reconstructed coal and rock fracture system is generated; Based on the reconstructed coal and rock fracture system and the nanoscale characteristic parameters, a target coal and rock pore network is generated. The step of generating the target coal pore network based on the reconstructed coal fracture system and the nanoscale characteristic parameters includes: Determine the porosity of the reconstructed coal-rock fracture system; The porosity of the reconstructed coal-rock fracture system is compared with the porosity of the fracture system in the micron-scale characteristic parameters. When the result of the difference processing is less than a preset threshold, the reconstructed coal and rock fracture system is taken as the target coal and rock fracture system; The target coal and rock fracture system is divided into blocks to obtain a generated space; wherein, the generated space includes multiple spatial coordinate points; Determine the state of each spatial coordinate point among the plurality of spatial coordinate points, so as to obtain the reconstructed coal and rock matrix porosity based on the state; Determine the porosity of the reconstructed coal and rock matrix pores; The porosity of the reconstructed coal matrix porosity and the matrix porosity in the nanoscale characteristic parameters are subjected to difference processing. When the result of the difference processing is less than a preset threshold, the reconstructed coal matrix porosity is taken as the target coal matrix porosity. Based on the target coal and rock fracture system and the target coal and rock matrix pores, a target coal and rock pore network is generated.

2. The method according to claim 1, characterized in that, The micrometer-scale characteristic parameters include at least one of the following: fracture system porosity, face cleavage spacing, face cleavage extension length, face cleavage aperture, face cleavage non-penetration probability, face cleavage roughness, end cleavage spacing, end cleavage aperture, end cleavage non-penetration probability, and end cleavage roughness; the nanometer-scale characteristic parameters include at least one of the following: matrix porosity and anisotropy characterization parameters.

3. The method according to claim 1, characterized in that, Determining the state of each spatial coordinate point among the plurality of spatial coordinate points includes: Acquire nanoscale binarized scanning images; By traversing the nanoscale binary scan image, multi-point neighborhood templates corresponding to the anisotropic characterization parameters are obtained; Based on the multi-point neighborhood template, the state of each spatial coordinate point among the multiple spatial coordinate points is determined one by one.

4. The method according to claim 1, characterized in that, The target coal and rock pore network was calibrated as a whole, and the calibration results were obtained, including: The target coal and rock pore network is calibrated as a whole according to the following formula: in, Let the box dimension be the target coal and rock. Let be the radius of the circle. For the radius is The minimum number of circles required to cover the target coal and rock pore network; The box dimension of the target coal and rock is used as the calibration result.

5. A calibration device for coal and rock pore networks, characterized in that, include: The acquisition module is used to acquire the initial fracture system of the target coal and rock, the initial matrix porosity of the target coal and rock, the preset main road equivalent rules, and the preset secondary road equivalent rules; An equivalent module is used to convert the initial fracture system into a main road according to the preset main road equivalence rule, and to convert the initial matrix pores into a secondary road according to the preset secondary road equivalence rule. The extraction module is used to extract micron-scale feature parameters from the main road and nanon-scale feature parameters from the secondary road; The generation module is used to generate a target coal and rock pore network based on the micron-scale characteristic parameters and the nano-scale characteristic parameters. The calibration module is used to perform overall calibration of the target coal and rock pore network and obtain calibration results; The step of generating the target coal and rock pore network based on the micrometer-scale characteristic parameters and the nanometer-scale characteristic parameters includes: Based on the micron-scale characteristic parameters, the growth space of the facet cleavage and the distribution characteristics of the facet cleavage in the growth space are determined; The end cleavage growth space is determined based on the end cleavage spacing and end cleavage opening in the micrometer-scale characteristic parameters. Based on the end cut non-penetration probability and end cut roughness in the micrometer-scale characteristic parameters, the distribution characteristics of the end cut in the end cut growth space are determined; Based on the distribution characteristics of the face cleavage growth space and the distribution characteristics of the end cleavage growth space, a reconstructed coal and rock fracture system is generated; Based on the reconstructed coal and rock fracture system and the nanoscale characteristic parameters, a target coal and rock pore network is generated. The step of generating the target coal pore network based on the reconstructed coal fracture system and the nanoscale characteristic parameters includes: Determine the porosity of the reconstructed coal-rock fracture system; The porosity of the reconstructed coal-rock fracture system is compared with the porosity of the fracture system in the micron-scale characteristic parameters. When the result of the difference processing is less than a preset threshold, the reconstructed coal and rock fracture system is taken as the target coal and rock fracture system; The target coal and rock fracture system is divided into blocks to obtain a generated space; wherein, the generated space includes multiple spatial coordinate points; Determine the state of each spatial coordinate point among the plurality of spatial coordinate points, so as to obtain the reconstructed coal and rock matrix porosity based on the state; Determine the porosity of the reconstructed coal and rock matrix pores; The porosity of the reconstructed coal matrix porosity and the matrix porosity in the nanoscale characteristic parameters are subjected to difference processing. When the result of the difference processing is less than a preset threshold, the reconstructed coal matrix porosity is taken as the target coal matrix porosity. Based on the target coal and rock fracture system and the target coal and rock matrix pores, a target coal and rock pore network is generated.

6. A calibration device for coal and rock pore networks, characterized in that, It includes a processor and a memory for storing processor-executable instructions, wherein the processor, when executing the instructions, implements the steps of the method according to any one of claims 1-4.

7. A computer-readable storage medium, characterized in that, It stores computer instructions that, when executed by a processor, implement the steps of the method according to any one of claims 1 to 4.

Images