Research Article  Open Access
Dayu Ye, Guannan Liu, Ning Luo, Feng Gao, Xinmin Zhu, Fengtian Yue, "Quantitative Analysis of the Topological Structure of Rock Pore Network", Geofluids, vol. 2021, Article ID 5517489, 11 pages, 2021. https://doi.org/10.1155/2021/5517489
Quantitative Analysis of the Topological Structure of Rock Pore Network
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 CTscan 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 scalefree 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 rockpore 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 (FIBSEM), 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 FIBSEM, and the coordinationnumber 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 porenetwork 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 twodimensional 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 smallworld model [24] and scalefree 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 threedimensional 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
Isolated pores in the rock networks do not affect the seepage process, and they are therefore removed during the construction of threedimensional 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.

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 scalefree 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 scalefree 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].
Where BA (Barabasi and Albert) scalefree network [25] was proposed in order to explain the generation mechanism of the powerlaw distribution. As a model for describing heterogeneous networks, the degree distribution function of BA scalefree network iswhere 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 PerronFrobenius 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 wellconnected 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
In order to describe the structure conveniently, the FruchtermanReingold (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 edgetopoint (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
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 fullycoupled 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
Previously published studies have typically considered the pore network of rocks to be randomly distributed. However, as shown in Figures 2 and 3, the scalefree network model which satisfies the powerlaw 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
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 rockpore structures and the connectivity laws of the rockpore 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 scalefree network model is more suitable for charactering the rockpore 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 wellconnected 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 fullycoupled network. As the distance between the pores increases, the possibility of the connection is greatly reduced. Besides, the rockpore 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/earthscience/research/researchgroups/porescalemodelling/microctimagesandnetworks/
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  porenetwork models and multiphase flow,” Current Opinion in Colloid & Interface Science, vol. 6, no. 3, pp. 197–207, 2001. View at: Publisher Site  Google Scholar
 P. E. Øren and S. Bakke, “Reconstruction of Berea sandstone and porescale modelling of wettability effects,” Journal of Petroleum Science & Engineering, vol. 39, no. 34, pp. 177–199, 2003. View at: Publisher Site  Google Scholar
 M. J. Blunt, B. Bijeljic, H. Dong et al., “Porescale imaging and modelling,” Advances in Water Resources, vol. 51, pp. 197–216, 2013. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 G. N. Liu, J. S. Liu, L. Liu, D. Ye, and F. Gao, “A fractal approach to fullycouple coal deformation and gas flow,” Fuel, vol. 240, no. 1, pp. 219–236, 2019. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 K. Li, S. Kong, P. Xia, and X. Wang, “Microstructural characterisation of organic matter pores in coalmeasure shale,” Advances in GeoEnergy Research, vol. 4, no. 4, pp. 372–391, 2020. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 S. Kelly, H. ElSobky, C. TorresVerdín, and M. T. Balhoff, “Assessing the utility of FIBSEM images for shale digital rock physics,” Advances in Water Resources, vol. 18, pp. 1–14, 2015. View at: Google Scholar
 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. View at: Publisher Site  Google Scholar
 Y. F. Yang, K. Wang, L. Zhang, H. Sun, K. Zhang, and J. Ma, “Porescale simulation of shale oil flow based on pore network model,” Fuel, vol. 251, pp. 683–692, 2019. View at: Publisher Site  Google Scholar
 A. Hajizadeh, A. Safekordi, and F. A. Farhadpour, “A multiplepoint statistics algorithm for 3D pore space reconstruction from 2D images,” Advances in Water Resources, vol. 34, no. 10, pp. 1256–1267, 2011. View at: Publisher Site  Google Scholar
 F. Civan, C. S. Rai, and C. H. Sondergeld, “Shalegas 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. View at: Publisher Site  Google Scholar
 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, 14 December, Okinawa, Japan, 2019. View at: Google Scholar
 X. Qin, Y. Zhou, and A. P. Sasmito, “An effective thermal conductivity model for fractal porous media with rough surfaces,” Advances in GeoEnergy Research, vol. 3, no. 2, pp. 149–155, 2019. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 M. J. Blunt, M. D. Jackson, M. Piri, and P. H. Valvatne, “Detailed physics, predictive capabilities and macroscopic consequences for porenetwork models of multiphase flow,” Advances in Water Resources, vol. 25, no. 812, pp. 1069–1089, 2002. View at: Publisher Site  Google Scholar
 H. Q. Zhao, I. F. Macdonald, and M. J. Kwiecien, “Multiorientation scanning: a necessity in the identification of pore necks in porous media by 3D computer reconstruction from serial section data,” Journal of Colloid and Interface Science, vol. 162, no. 2, pp. 390–401, 1994. View at: Publisher Site  Google Scholar
 D. J. Watts and S. H. Strogatz, “Collective dynamics of ‘smallworld’ networks,” Nature, vol. 393, no. 6684, pp. 440–442, 1998. View at: Publisher Site  Google Scholar
 A.L. Barabasiand and R. Albert, “Emergence of scaling in random networks,” Science, vol. 286, no. 5439, pp. 509–512, 1999. View at: Publisher Site  Google Scholar
 A. Sanyal, B. R. Lajoie, G. Jain, and J. Dekker, “The longrange interaction landscape of gene promoters,” Nature, vol. 389, pp. 109–113, 2012. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 J. Zeng, X. Wang, J. Guo, F. Zeng, and Q. Zhang, “Composite linear flow model for multifractured horizontal wells in tight sand reservoirs with the threshold pressure gradient,” Journal of Petroleum Science and Engineering, vol. 165, pp. 890–912, 2018. View at: Publisher Site  Google Scholar
 MicroCT Images and Networks, http://www.imperial.ac.uk/earthscience/research/researchgroups/porescalemodelling/microctimagesandnetworks/(accessed on 14 March 2021).
 P. Iassonov, T. Gebrenegus, and M. Tuller, “Segmentation of Xray 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. View at: Google Scholar
 W. L. Zhu, G. A. Gaetani, F. Fusseis, L. G. J. Montesi, and F. de Carlo, “Microtomography of partially molten rocks: threedimensional melt distribution in mantle peridotite,” Science, vol. 332, no. 6025, pp. 88–91, 2011. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 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 forcedirected placement,” Software Practice & Experience, vol. 21, no. 11, pp. 1129–1164, 1991. View at: Publisher Site  Google Scholar
Copyright
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.