Geofluids
Volume 2021 (2021), Article ID 5517489, 11 pages
https://doi.org/10.1155/2021/5517489

Research Article

Quantitative Analysis of the Topological Structure of Rock Pore Network

Dayu Ye,1,2,3 Guannan Liu,2,3 Ning Luo,1,3 Feng Gao,1,3,4 Xinmin Zhu,1,3 and Fengtian Yue3,4

¹State Key Laboratory for Geomechanics and Deep Underground Engineering, China University of Mining and Technology, Xuzhou 221116, China
²Key Laboratory of Deep Earth Science and Engineering (Sichuan University), Ministry of Education, Chengdu, 610065, China
³Mechanics and Civil Engineering Institute, China University of Mining and Technology, Xuzhou 221116, China
⁴Laboratory of Mine Cooling and Coal-Heat Integrated Exploitation, China University of Mining and Technology, Xuzhou 221116, China

Correspondence should be addressed to Guannan Liu; guannanliu@cumt.edu.cn

Received 9 January 2021; Revised 18 April 2021; Accepted 27 May 2021; Published 14 June 2021

Academic Editor: Jian-Chun Xu

Copyright © 2021 Dayu Ye et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

As the most significant nonlinear reservoir, the rocks have complex structural characteristic. The pore structure of the rock is varied in shape and complex in connectivity. However, the prevailing methods for characterising the microstructure of rocks, such as the coordination number method and fractal theory, are still difficult to quantify the structural properties. In this study, based on the CT-scan method and a new complex network theory, the topological characteristics of rocks such as seepage path selection, degree of pore aggregation, pore importance, and pore module structure are analysed. The results show that the scale-free network model is more reliable in characterising the rock pore network than previously published structural models, and a small number of pores are the “key” to the seepage process. Besides, we proposed a new method to quantify the importance of rock pores and present the distribution characteristics and connectivity laws of the rock-pore network. This provides a new method to study the seepage process of the nonlinear reservoirs.

1. Introduction

The highly complex network of pore throats is the main conduit for CBM migration during the extraction [1–4]. Therefore, analysing the topological characteristics of rock pore networks, including morphology, connectivity, and seepage path selection, is one of the main approaches to unravelling the mechanism of hydrocarbon transport and to enhancing the recovery [5–8].

In recently published studies, many scholars have studied the pore network structure and seepage properties of rocks at the microscale [9–12]. Advances in CT and scanning electron microscopy (SEM) have enabled the application of the rock digital core. Using focused ion beam scanning electron microscope (FIB-SEM), scanning electron microscope (SEM), and computed tomography (CT), Leu et al. studied the pore size distribution, orientation, and mineralogical characteristics of shale and mudstone [13]. The digital core of shale with different porosity was extracted by Kelly et al. using FIB-SEM, and the coordination-number method is used to explore the network connectivity of shale pores [14]. Based on digital cores of shale, Yang et al. proposed a superposition algorithm capable of analysing both organic and inorganic pores of the matrix [15].

Furthermore, analyses related to the structural characteristics of rocks have been progressively recognized, and methods for the construction of rock pore-network models have been widely studied. Yang et al. proposed a pore network model to analyse the evolution of shale oil permeability under different organic matter contents [16]. Hajizadeh et al. propose a stochastic multidimensional numerical algorithm based on continuous two-dimensional multipoint statistical simulations coupled with multiscale conditional data using 2D images of the pore structure [17]. Civan et al. propose an improved shale gas transport model that can be used to determine the diffusion coefficient and permeability of shales [18]. Additionally, many published models have investigated the seepage properties of pore rocks which do not satisfy the Darcy’s law of permeability [19, 20]. However, the throat connecting the rock pores is irregular and unevenly distributed. Besides, current methods for characterising rock structure mainly adopt the Bieuler number or macroscopic averaging perspectives. The rock pore structure is highly complex, and the conventional network models still have major difficulties in characterising pore connectivity and other topologies [21–23].

With the advent of the small-world model [24] and scale-free model [25], the theoretical framework of modern complex networks was initially formed. In the past 10 years, as a new interdisciplinary field, complex network theory has penetrated into many disciplines from life sciences to physics [26, 27]. Complex network theory reasonably characterises the structural complexity of network systems, the pore complexity, and the interaction of various complexity factors [28]. At present, studies of rock structures are mainly based on basic topological parameters such as macroscopic averaging or allotment numbers [29–32]. However, it is still difficult to characterise the connectivity of pore networks in a reasonably reliable way. Based on complex network theory, this paper presents a comprehensive analysis of the distribution patterns of rock networks such as average path length, degree distribution, the importance of pores in the seepage process, and the characteristics of pore aggregation. And the results of this study provide a rational approach to reveal the micromechanics of porous rocks.

2. Basic Structure of Rock Pore Network

In this section, the digital cores used for our structural analysis were constructed by the Blunt team at Imperial College [33]. Six different rock samples were selected, and the pore structures were extracted from the rock samples. The three-dimensional pore networks of the rock samples are shown in Figure 1 [34–36].

(a) Berea sandstone

(b) Carbonate C1

(c) Carbonate C2

(d) Sand pack F1

(e) Sandstone S8

(f) Synthetic silica A1

Figure 1: 3D pore network of different rocks (bounding box volume: 300 μm by 300 μm by 300 μm).

Isolated pores in the rock networks do not affect the seepage process, and they are therefore removed during the construction of three-dimensional pore networks. Based on complex network theory, we obtained the basic data of the network structures for the rocks above, as shown in Table 1.

Table 1: Basic information of pore network.

Rock type
Berea sandstone
Carbonate C1
Carbonate C2
Sand pack F1
Sandstone S8
Synthetic silica A1

Porosity
19.6%
23.3%
16.8%
33%
34%
42.9%

Node number
6057
3962
1933
974
1266
1620

Edge number
12098
6700
1932
2651
3833
1619

Average degree
3.995
3.382
1.999
5.444
6.055
1.999

Network density
0.001
0.001
0.001
0.006
0.005
0.001

Network diameter
33
31
137
16
17
48

Average path length
12.445
9.981
45.957
6.399
6.452
22.397

Modularity
0.868
0.919
0.953
0.726
0.745
0.949

Average clustering coefficient
0.296
0.408
0
0.317
0.425
0

Meanwhile, based on complex network theory, we analysed the structural characteristics of the rock pore networks [37]. The degree of a pore is defined as the number of the pore connected with other pores:

where is the average degree of the network, is the number of throats, and is the number of pores in the network.

In order to calculate the transmission efficiency of the rock network, the network diameter D and the average path length are used, which represent the maximum value of the distance between any two pores in the rock network and the average of the distance between any two pores in the network:

where is the value of the shortest path length.

Besides, the average path length can be expressed as:

Theoretical experience suggests that the network diameter and the average path length should be proportional to network size and porosity. However, the pattern in the actual rock samples (Table 1) does not behave in line with theoretical experience. The above experience applies to the homogeneous rock network models, while in the real rock models, the network diameter and average path length could be influenced by a variety of factors such as the number of throats, pore size, and the length of throats.

We applied network density to quantify the denseness of rock pores. The pore network density is the ratio of the number of throats to its upper limit:

In order to analyse whether the network can be divided into different modules, we calculated the modularity of the pore network:

where is the proportion of the throats in the module and the module to the whole rock throats, . And a module is a group of pores which is physically or functionally connected, working together to perform a relatively independent function.

Besides, we applied clustering coefficient to characterise the aggregation of rock pores. The clustering coefficient can be expressed as:

where is the degree of pore , i.e., there are pores connected to pore . If there are interconnected throats among these pores, the number of interconnected throats among the pores is denoted by .

3. Sandstone Pore Network Topology Analysis

In order to analyse the connectivity of the pore network, we calculated the degree distribution for different rock samples, as shown in Figure 2. The degree distributions of the rock networks fit well with the BA scale-free model. There are no eigenvalues for the rock pore networks, and the degree of the network for different samples is mainly less than 5, which can be regarded as the “degree threshold.” When the degree exceeds the “degree threshold,” the number of pores is drastically reduced. This inhomogeneity is probably responsible for the characteristics: a small number of pores (with large degree values) play a decisive role in the seepage process throughout the rock network. Besides, in our previous published research, we have concluded that the scale-free network model is more reliable in characterising the degree distribution of rocks than random networks and other mainstream network methods, as shown in Figure 3 [38].

Figure 2: Degree distribution of three-dimensional pore networks.

Figure 3: Comparison of different models with real-rock distribution characteristics (note: the degree distribution of regular model is the points of , 3, 4, and 6).

Where BA (Barabasi and Albert) scale-free network [25] was proposed in order to explain the generation mechanism of the power-law distribution. As a model for describing heterogeneous networks, the degree distribution function of BA scale-free network is

where is the number of new nodes connected to existing nodes during network expansion, namely,

Among them, is power exponent. Besides, the degree distribution of ER random model satisfies [39]:

And the construction algorithm of WS small world model is as follows [24]:

(1)Start Modeling with a Ring of Regular Networks. The network contains nodes, and each node connects to the nodes closest to it. Where (2)Randomize Reconnection. Reconnect each edge in the network randomly with probability , that is, one endpoint of the edge remains unchanged, while the other endpoint is taken as a randomly selected node in the network. It stipulates that there can be at most one edge between any two different nodes, and each node cannot have an edge connected to itself

In order to measure the importance of pores in seepage networks, eigenvector centrality is used. For two pores and in a rock network, when and are connected, otherwise, . Then, the function of eigenvector centrality of pore can be expressed as:

where and are the eigenvector centrality of pores and , is the set of pores connected to pore , is the total number of pores in the rock, and is the scaling constant. Denoting the adjacency matrix of the network by , Eq. (10) can be expressed in matrix form as:

In this way, solving for eigenvector centrality can be transformed into a problem of computing eigenvalues and vectors or solving a system of equations. According to the Perron-Frobenius theorem, the vector can only be the eigenvector corresponding to the largest eigenvalue of the adjacency matrix . Therefore, the eigenvector centrality of each pore in the rock network can be obtained from the eigenvector corresponding to the largest eigenvalue derived by the curtain iteration method.

The eigenvector centralities of the pore networks are shown in Figure 4, where the abscissa of the figure is the pore Id. For the different rock pore networks, the eigenvector centrality of most pores is less than 0.6. Only a small proportion of the pores have a large centrality. The importance of a pore depends on both the number and the importance of its connected pores. The eigenvector centrality method defines the centrality of a pore as a function of the centrality of the pores it connects. Thus, as the eigenvector centrality increases, the connectivity of the pores in the pore network subsequently increases. In rock seepage network, if the well-connected pores are blocked with rock debris or sealed by the crust, it may result in a significant reduction in the overall connectivity of the rock network. We therefore speculate that pores with large eigenvector centrality have a more significant contribution to rock permeability.

(a) Berea sandstone

(b) Carbonate C1

(c) Carbonate C2

(d) Sand pack F1

(e) Sandstone S8

(f) Synthetic silica A1

Figure 4: Eigenvector centrality distribution of pores.

In order to describe the structure conveniently, the Fruchterman-Reingold (FR) method is used to assign a uniform circular surface. The FR algorithm considers all nodes as electrons, and each node receives two forces: (1) the Coulomb force of the other nodes (repulsion) and (2) the Hooke force of the edge-to-point (gravity). Then, with the interaction of the forces, the whole layout will eventually be called a state of equilibrium [40]. The pore centrality of the networks is shown in Figure 5, where the number is the Id of the pores, and the size of the pores expresses the centrality.

(a) Berea sandstone

(b) Carbonate C1

(c) Carbonate C2

(d) Sand pack

(e) Sandstone S8

(f) Synthetic silica A1

Figure 5: Centrality of rock pores.

The clustering coefficients of the pores for different rock samples are shown in Figure 6. It can be concluded from Figure 6 that the increase of clustering coefficients was initially slow and then became evident. The clustering coefficients for the different networks are distributed between 0 and 0.6. And in conjunction with the conclusion of the analysis in Table 1, we concluded that the connectivity of real rocks is always less than that of the fully-coupled network with a clustering value of 1. As the distance between the pores increases, the possibility of the connection is greatly reduced.

(a) Berea sandstone

(b) Carbonate C1

(c) Carbonate C2

(d) Sand pack

(e) Sandstone S8

(f) Synthetic silica A1

Figure 6: Clustering coefficient of pores.

Previously published studies have typically considered the pore network of rocks to be randomly distributed. However, as shown in Figures 2 and 3, the scale-free network model which satisfies the power-law distribution is more suitable for characterising the degree distribution. Therefore, we simultaneously calculated the modularity of different pore networks to compare the differences between the real rock and the stochastic model in the pore distribution characteristics. The modular characteristics of the different rock samples are shown in Figure 7, where pores in the same module are represented by the same colour and size.

(a) Berea sandstone

(b) Carbonate C1

(c) Carbonate C2

(d) Sand pack

(e) Sandstone S8

(f) Synthetic silica A1

Figure 7: Modularity structure of rock pore network.

It can be concluded from Figure 7 that the pore distributions of the different rocks have obvious modular characteristics. Besides, the modularity of different rock samples ranged from 0.8 to 1 (Table 1), which also demonstrates the modular characteristic of the pore network. We therefore conclude that the pores of the rocks can be divided into groups based on their characteristic and function. In addition, the pore connectivity is proportional to the number of throat connecting the pores. Therefore, we predict that the pore groups with high connectivity have a tendency to the high permeability.

4. Conclusions

In this paper, we proposed a new method to analyse the distribution characteristics of rock-pore structures and the connectivity laws of the rock-pore network. The pore networks were studied by four parameters of the complex theory: degree distribution, eigenvector centrality, clustering coefficient, and the modularity. The main conclusions are shown as follows:

(1)Compared with the random model, the BA scale-free network model is more suitable for charactering the rock-pore networks. And a small number of pores (with large degree values) play a decisive role in the seepage process throughout the rock network(2)The importance of a pore depends on both the number and the importance of its connected pores. Accordingly, we propose a new method for quantifying the importance of rock pores. If the well-connected pores are blocked with rock debris or sealed by the crust, it may result in a significant reduction in the overall connectivity of the rock network(3)The connectivity of real rocks is always less than that of the fully-coupled network. As the distance between the pores increases, the possibility of the connection is greatly reduced. Besides, the rock-pore networks have obvious modular characteristics. The networks can be divided into several groups according to the connectivity. Pore connectivity is proportional to the number of throat connecting the pores. We therefore predict that the pore groups with high connectivity have the high permeability

Data Availability

All data included in this study are available upon request by contact with the corresponding author. The digital cores of Blunt can be found in the website of Imperial College London:http://www.imperial.ac.uk/earth-science/research/research-groups/pore-scale-modelling/micro-ct-images-and-networks/

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the open foundation of Key Laboratory of Deep Earth Science and Engineering (Sichuan University), Ministry of Education No. DESE202103; the National Key Research and Development, the Program of China No. 2020YFA0711800; and the National Natural Science Foundation of China No. 12072363.

References

M. J. Blunt, “Flow in porous media - pore-network models and multiphase flow,” Current Opinion in Colloid & Interface Science, vol. 6, no. 3, pp. 197–207, 2001.

P. E. Øren and S. Bakke, “Reconstruction of Berea sandstone and pore-scale modelling of wettability effects,” Journal of Petroleum Science & Engineering, vol. 39, no. 3-4, pp. 177–199, 2003.

M. J. Blunt, B. Bijeljic, H. Dong et al., “Pore-scale imaging and modelling,” Advances in Water Resources, vol. 51, pp. 197–216, 2013.

G. N. Liu, D. Y. Ye, F. Gao, and J. S. Liu, “A dual fractal poroelastic model for characterizing fluid flow in fractured coal masses,” Geofluids, vol. 2020, Article ID 2787903, 13 pages, 2020.

Y. Ju, Y. Yang, Z. Song, and W. Xu, “A statistical model for porous structure of rocks,” Sci. China Ser. E, vol. 51, no. 11, pp. 2040–2058, 2008.

P. E. Øren and S. Bakke, “Process based reconstruction of sandstones and prediction of transport properties,” Transport in Porous Media, vol. 46, no. 2/3, pp. 311–343, 2002.

X. C. Zhao, The Study of Reconstruction Method of Digital Core and Pore Network Model, China Petroleum University, East China, 2009.

W. Li, J. S. Liu, J. Zeng et al., “A fully coupled multidomain and multiphysics model for evaluation of shale gas extraction,” Fuel, vol. 278, p. 118214, 2020.

G. LIU, B. YU, F. GAO, D. YE, and F. YUE, “Analysis of permeability evolution characteristics based on dual fractal coupling model for coal seam,” Fractals, vol. 28, no. 7, article 2050133, 2020.

G. N. Liu, J. S. Liu, L. Liu, D. Ye, and F. Gao, “A fractal approach to fully-couple coal deformation and gas flow,” Fuel, vol. 240, no. 1, pp. 219–236, 2019.

J. Zeng, J. Liu, W. Li, Y. K. Leong, D. Elsworth, and J. Guo, “Evolution of shale permeability under the influence of gas diffusion from the fracture wall into the matrix,” Energy & Fuels, vol. 34, no. 4, pp. 4393–4406, 2020.

K. Li, S. Kong, P. Xia, and X. Wang, “Microstructural characterisation of organic matter pores in coal-measure shale,” Advances in Geo-Energy Research, vol. 4, no. 4, pp. 372–391, 2020.

L. Leu, A. Georgiadis, M. J. Blunt et al., “Multiscale description of shale pore systems by scanning SAXS and WAXS microscopy,” Energy & Fuels, vol. 30, no. 12, pp. 10282–10297, 2016.

S. Kelly, H. El-Sobky, C. Torres-Verdín, and M. T. Balhoff, “Assessing the utility of FIB-SEM images for shale digital rock physics,” Advances in Water Resources, vol. 18, pp. 1–14, 2015.

Y. F. Yang, J. Yao, C. C. Wang et al., “New pore space characterization method of shale matrix formation by considering organic and inorganic pores,” Journal of Natural Gas Science and Engineering, vol. 27, pp. 496–503, 2015.

Y. F. Yang, K. Wang, L. Zhang, H. Sun, K. Zhang, and J. Ma, “Pore-scale simulation of shale oil flow based on pore network model,” Fuel, vol. 251, pp. 683–692, 2019.

A. Hajizadeh, A. Safekordi, and F. A. Farhadpour, “A multiple-point statistics algorithm for 3D pore space reconstruction from 2D images,” Advances in Water Resources, vol. 34, no. 10, pp. 1256–1267, 2011.

F. Civan, C. S. Rai, and C. H. Sondergeld, “Shale-gas permeability and diffusivity inferred by improved formulation of relevant retention and transport mechanisms,” Transport in Porous Media, vol. 86, no. 3, pp. 925–944, 2011.

W. Li, J. Liu, J. Zeng, Y. K. Leong, and D. Elsworth, “A fully coupled multidomain and multiphysics model for shale gas production,” in Presented at the 5th ISRM Young Scholars' Symposium on Rock Mechanics and International Symposium on Rock Engineering for Innovative Future, 1-4 December, Okinawa, Japan, 2019.

X. Qin, Y. Zhou, and A. P. Sasmito, “An effective thermal conductivity model for fractal porous media with rough surfaces,” Advances in Geo-Energy Research, vol. 3, no. 2, pp. 149–155, 2019.

L. A. Dillard and M. J. Blunt, “Development of a pore network simulation model to study nonaqueous phase liquid dissolution,” Water Resources Research, vol. 36, no. 2, pp. 439–454, 2000.

M. J. Blunt, M. D. Jackson, M. Piri, and P. H. Valvatne, “Detailed physics, predictive capabilities and macroscopic consequences for pore-network models of multiphase flow,” Advances in Water Resources, vol. 25, no. 8-12, pp. 1069–1089, 2002.

H. Q. Zhao, I. F. Macdonald, and M. J. Kwiecien, “Multi-orientation scanning: a necessity in the identification of pore necks in porous media by 3-D computer reconstruction from serial section data,” Journal of Colloid and Interface Science, vol. 162, no. 2, pp. 390–401, 1994.

D. J. Watts and S. H. Strogatz, “Collective dynamics of ‘small-world’ networks,” Nature, vol. 393, no. 6684, pp. 440–442, 1998.

A.-L. Barabasiand and R. Albert, “Emergence of scaling in random networks,” Science, vol. 286, no. 5439, pp. 509–512, 1999.

A. Sanyal, B. R. Lajoie, G. Jain, and J. Dekker, “The long-range interaction landscape of gene promoters,” Nature, vol. 389, pp. 109–113, 2012.

Z. Jiang, M. Liang, Q. Li, and D. Guo, “Optimal dynamic bandwidth allocation for complex networks,” Phys. A Stat. Mech. Appl., vol. 392, no. 5, pp. 1256–1262, 2013.

A. L. Barabasi, Linked: The New Science of Networks, Persus Publishing, Massachusetts, 2002.

J. YAO, C. C. WANG, D. Y. FAN, Z. Q. HUANG, H. SUN, and Y. F. YANG, “The construction method and microscopic flow simulation of carbonate dual pore network model,” Scientia Sinica Physica, Mechanica & Astronomica, vol. 43, no. 7, pp. 896–902, 2013.

Z. Zhang, Y. Shi, H. Li, and W. Jin, “Experimental study on the pore structure characteristics of tight sandstone reservoirs in Upper Triassic Ordos Basin China,” Energy Exploration & Exploitation, vol. 34, no. 3, pp. 418–439, 2016.

J. Zeng, X. Wang, J. Guo, F. Zeng, and F. Zeng, “Modeling of heterogeneous reservoirs with damaged hydraulic fractures,” Journal of Hydrology, vol. 574, pp. 774–793, 2019.

J. Zeng, X. Wang, J. Guo, F. Zeng, and Q. Zhang, “Composite linear flow model for multi-fractured horizontal wells in tight sand reservoirs with the threshold pressure gradient,” Journal of Petroleum Science and Engineering, vol. 165, pp. 890–912, 2018.

Micro-CT Images and Networks, http://www.imperial.ac.uk/earth-science/research/research-groups/pore-scale-modelling/micro-ct-images-and-networks/(accessed on 14 March 2021).

P. Iassonov, T. Gebrenegus, and M. Tuller, “Segmentation of X-ray computed tomography images of porous materials: a crucial step for characterization and quantitative analysis of pore structures,” Water Resources Research, vol. 45, no. 9, article W09415, 2009.

W. L. Zhu, G. A. Gaetani, F. Fusseis, L. G. J. Montesi, and F. de Carlo, “Microtomography of partially molten rocks: three-dimensional melt distribution in mantle peridotite,” Science, vol. 332, no. 6025, pp. 88–91, 2011.

F. Zhang, C. Hui, and J. S. Terblanche, “An interaction switch predicts the nested architecture of mutualistic networks,” Ecology Letters, vol. 14, no. 8, pp. 797–803, 2011.

G. Liu, J. He, R. Li et al., “Topological and dynamic complexity of rock masses based on GIS and complex networks,” Physica A: Statistical Mechanics and its Applications, vol. 512, pp. 1240–1248, 2018.

D. Ye, G. Liu, F. Gao, X. Zhu, and Y. Hu, “A study on the structure of rock engineering coatings based on complex network theory,” Coatings, vol. 10, no. 12, p. 1152, 2020.

B. Bollobas, Random Graphs, New York Academic Press, Washington, DC, USA, 2nd edition, 2011.

T. M. J. Fruchterman and E. M. Reingold, “Graph drawing by force-directed placement,” Software Practice & Experience, vol. 21, no. 11, pp. 1129–1164, 1991.