Multiscale and Multiphysical Approaches to Fluids Flow in Unconventional Reservoirs
View this Special IssueResearch Article  Open Access
Hang Bian, Yuxuan Xia, Cheng Lu, Xuwen Qin, Qingbang Meng, Hongfeng Lu, "Pore Structure Fractal Characterization and Permeability Simulation of Natural Gas Hydrate Reservoir Based on CT Images", Geofluids, vol. 2020, Article ID 6934691, 9 pages, 2020. https://doi.org/10.1155/2020/6934691
Pore Structure Fractal Characterization and Permeability Simulation of Natural Gas Hydrate Reservoir Based on CT Images
Abstract
The gaswater twophase seepage process is complex during the depressurization process of natural gas hydrate in a clayey silt reservoir in the South China Sea, the transport mechanism of which has not been clarified as it is affected by the pore structure. In this study, we select six clayey silt samples formed after the dissociation of natural gas hydrate in the South China Sea, employing CT scanning technology to observe the pore structure of clayey silt porous media directly. The original CT scanning images are further processed to get the binarized images of the samples, which can be used for simulation of the porosity and absolute permeability. Based on the fractal geometry theory, pore structures of the samples are quantitatively characterized from the aspect of pore distribution, heterogeneity, and anisotropy (represented by three main fractal geometric parameters: fractal dimension, lacunarity, and succolarity, respectively). As a comparison, the binarized CT images of two conventional sandstone cores are simulated with the same parameters. The results show that the correlation between porosity and permeability of the hydrate samples is poor, while there is a strong correlation among the succolarity and the permeability. Fractal dimension (represents complexity) of clayey silt samples is higher compared with conventional sandstone cores. Lacunarity explains the difference in permeability among samples from the perspective of pore throat diameter and connectivity. Succolarity indicates the extent to which the fluid in the pore is permeable, which can be used to characterize the anisotropy of pore structures. Therefore, these three fractal parameters clarify the relationship between the microstructure and macroscopic physical properties of clayey silt porous media.
1. Introduction
Natural gas hydrate refers to one kind of icelike crystalline compound formed by a methanebased hydrocarbon gas and water in a stable domain controlled by a certain temperature and pressure condition [1]. In the field of exploration and development of the natural gas hydrate, especially in hydrate classification, extraction, dissociation, experimental simulation, productivity evaluation method, etc., the work has become a research hotspot pursued by scientific research institutions, experts, and scholars [2–8]. The first offshore natural gas hydrate production test was conducted by China Geological Survey at Shenhu Area in South China Sea in 2017 [9]. During the depressurization process of production, the solid hydrate is decomposed into methane gas and water. The twophase flow passes through the formation and then flows into the wellbore. In the initial stage, the flow rate is stable. However, as production progresses, the flow rate is going to be greatly reduced, resulting in a significant decrease in production. The abovementioned multiphase flow within pore space is extremely complicated, and its seepage mechanism is not yet clear. It is speculated to be related to pore structure and mineral composition of the hydrate reservoir.
The differences of porosity, connectivity, and topology of pore space always determine the macroscopic seepage properties [10, 11]. The pore space of the reservoir core can be obtained by multilaboratory test methods [12–15]. The obtained pore space images and parameters can also lay the foundation for calculation of permeability [16–18]. Many numerical simulation methods [19–21] were used to obtain permeability of the reservoir, which greatly promotes understanding of the relationship between reservoir porosity and permeability, and the formulation of reservoir development programs.
Fractal geometry theory, as the powerful tool for quantification of complex space, has been used by many scholars in the study of characterization of the pore structure and modelling of transport phenomena in conventional and unconventional reservoirs. Yu and Cheng [22] established the fractal permeability model for the bidispersed porous media and introduced the algorithm of pore space fractal dimension and tortuosity fractal dimension based on the boxcounting method. Cai et al. [23] established a fractal permeability model for creeping microstructures under different pressure conditions based on the data of water flooding experiment combined with CT scanning [24] for clayey silt sediments. Liu et al. [25] used fractal dimensions to establish a fractal theorybased relative permeability model of hydratebearing sediments. N'Diaye et al. [26] computed fractal dimension, lacunarity, and succolarity to provide a precise description of porosity and pore characteristics. Xia et al. [27] proposed that succolarity can be used to predict the permeability of low permeable sandstone reservoirs. Therefore, fractal geometry theory has been proved to be a good method to characterize the reservoir microstructure and to study the macroscopic differences in physical properties of reservoirs. In the study of clayey silt porous media, the macroscopic changes of reservoir physical properties can also be explained from the perspective of microstructure based on fractal geometry theory. However, few studies have used fractal quantitative parameters to characterize the microporous structure of clayey silt porous media and further used to study the change of permeability.
In this study, pore structure images of six hydrate reservoir samples were firstly acquired by microCT scanning. Then, the porosity, permeability, and fractal parameters were obtained based on the binarized pore images. Finally, how microporous structures affect the permeability of clayey silt reservoir samples was analyzed according to the relationship between permeability and pore structure parameters.
2. Materials and Facility
2.1. Hydrate Reservoir Samples
Six clayey silt samples are taken from the field of the first offshore natural gas hydrate production test, which is located at Shenhu Area in the South China Sea [9]. There are eight natural gas hydrate deposits, among which the W1117 deposit (Figure 1) is selected as the optimal target for the production test. The lithology of the reservoir is clayey silt, of which the clay mineral content is high (26%30%), mainly composed of montmorillonite and illite. The average median diameter of the sample particle is 12 μm. The average values of effective porosity and hydrate saturation are 33% and 31%, respectively. Pore space of the reservoir is filled with three phases: solid hydrate, free gas, and water. Detailed geological, ore body, seismic, and logging information of reservoir can be found in literature [9].
2.2. CT (Computed Tomography) Scanner
Six clayey silt samples taken from the hydrate reservoir are subjected to CT scanning by a MicroSCT200 micron CT scanner. The scanning voltage and power are 95 kV and 10 W, respectively. The scanning mode is focused scanning, and all clayey silt samples are 5 mm in diameter and 5 mm in length, so the pixel number and pixel resolution for CT images of all clayey silt samples are and ~1 μm, separately.
3. Methods
After the original CT scanning images obtained, a series of preprocessing including extraction, filtering, and binarization [28] is conducted to acquire the binarized pore space structure images. Based on the binarized images, a variety of data such as sample porosity, permeability, and fractal structural parameters can be acquired and then can be used to construct a relational model between pore structure and physical properties.
3.1. Porosity Calculation
Porosity refers to the ratio of the pore volume to the total volume of the rock. It is an important parameter to quantitatively describe the rock reservoir capacity. The expression is shown below: where is the pore volume, is the total volume of the rock, and is the total volume of the skeleton.
In microscopic sample images, the pore volume and the rock skeleton volume correspond to the pixel volume of different color regions. Since the size represented by each pixel is uniform in all images of the same resolution, Equation (1) can be transformed to where is the number of pore pixels, is the total number of pixels of the sample image, and is the number of sample skeleton pixels. The resulting binarized image is converted into an array matrix as displayed in Figure 2(b), and the sample image is represented by the pixel value of each pixel in the image. A sample image of a cube corresponds to a threedimensional array matrix. By counting the number “0” and “1” in the matrix, the number of pore pixels and rock skeleton pixels can be obtained, respectively.
(a)
(b)
3.2. Fractal Structural Parameters Calculation
In the fractal geometry theory, fractal dimension, lacunarity, and succolarity are the main parameters of fractal analysis. Fractal dimension is used as the basic parameter of fractal quantitative characterization. Among plentiful research results, the boxcounting method is a widely used and accurate method by which the fractal dimension of pore distribution can be obtained based on the structural images of the tested sample. In this study, all fractal dimensions are calculated based on the boxcounting method as well. In the boxcounting method, the image is covered by boxes with scale , and the number of the boxes covered with the pores is recorded. The expression is as follows:
Therefore, the negative value of the slope of the linear fit for the boxnumber and scale of double logarithmic plots in the selfsimilar interval is the pore space fractal dimension .
Lacunarity analysis is a popular analytical tool of fractal microstructure in recent years. It can quantitatively analyze the degree of fractal clustering and characterize the changes of natural images and can consequently distinguish different structures with the same fractal dimension. Lacunarity characterizes the distribution of the observed objects in the image and measures the frequency and magnitude of observed objects. The larger the lacunarity, the more intensive the distribution of the observed object.
The gliding boxcounting method is a general measurement method of lacunarity. The value of lacunarity can be defined as a statistical moment function when divided by the squared value of [29, 30]:
Equation (5) is used to normalize lacunarity, which can exclude the influence of porosity.
Succolarity is one of the crucial parameters of fractal geometry, which indicates the ability of fluid flowing inside the medium, that is, the extent to which the fluid in the pore is permeable or flowing. It can be used to measure the connectivity of the patterns or structures in different directions.
The algorithm of succolarity is also based on the boxcounting method, which evaluates the fluid flow ability in different flow directions by applying a virtual pressure field. The value of succolarity can be characterized as the pressure coverage value divided by the maximum possible coverage value [31]:
3.3. Permeability Calculation
In order to study the relationship between complex microstructure and macroscopic seepage property of the gas hydrate reservoir, we selected the absolute permeability, which is associated with the sample structure, as the research object. Absolute permeability of the rock core is the permeability measured when only one fluid (single phase) flows through the pores of the rock, the fluid does not react with the rock, and the flow of fluid is based on Darcy’s law of linear percolation. When a singlephase fluid with viscosity flows through a porous medium whose crosssection is , length is , and pressure difference is , the amount of fluid flowing through the rock pores within a unit time is
Based on Equation (7), the absolute permeability of the six samples was simulated using the finite volume method (FVM) [32]. Avizo® 9 software is employed to convert the six binary graphic models into FVM models for the absolute permeability simulation. Simulation work is done with the Avizo XLab Hydro Extension module which assumes that the fluid was laminar and the viscosity was 0.001 Pa·s. The velocitypressure coupling solver is used, and the absolute convergence factor is . The fluid flow governing equation is the NavierStokes (NS) equation with the pressure gradient along the direction of calculated permeability. Fluid circulates only in the core, and the core boundary is noncirculating. After completion of the iterations, the flow at the outlet can be obtained. Then, the permeability of each rock sample can be calculated by Darcy’s law.
4. Results and Discussion
Figure 3 (G1G6) shows original grayscale images of the six samples. The pores of the hydrate reservoir samples in Shenhu Area possess ringlike boundaries, which are derived from microorganisms such as diatoms. After a series of preprocessing procedures conducted, the binarized pore structure images of the six samples are obtained, which are expressed in Figure 3 (B1B6).
In order to better illustrate the special physical properties of hydrate reservoir samples, we selected two conventional sandstone cores for comparison. The pixel number and pixel resolution for CT images of sandstone cores shown in Figure 4 are 512^{3} and ~1 μm. According to Figures 3 and 4, unlike the irregular intergranular pores of conventional sandstone reservoirs, the pore morphology of the hydrate reservoir mostly turns out to be spherical and uniform in the image. The reason is that the hydrate reservoirs are mainly composed of fine silt particles with good sorting, so that the intergranular pores are consequently regular. In addition, the hydrate reservoir is greatly affected by microbial processes during diagenesis, which leads to the existence of spheroidal bioclastics and interdental pores.
(a)
(b)
4.1. Porosity and Permeability
After the binarized images shown in Figures 3 and 4 obtained, porosity and permeability can be calculated based on these images. Figure 5 illustrates the pressure field distributed on direction of permeability simulation of the six hydrate samples. Table 1 shows the porosity and absolute permeability of hydrate Samples 16 and conventional sandstone Cores a and b. The ranges of porosity and absolute permeability of eight samples and cores are 0.17~0.39 and 18.48~1497.73 mD, respectively. The simulation results on Core 1 , Core 3 , and Core 6 express that the pressure field and the velocity field are opposite in direction, which results in the inability to calculate the permeability. And the cause of this condition is presumed to be the absence of connected pores on these directions. According to the data in Table 1, the hydrate sample exhibits higher porosity and lower permeability compared with the conventional sandstone core, and the porosity and permeability are poorly correlated (as shown in Figure 6, the correlation coefficient is only 0.4669). That is, the porosity of Sample 1, 2, and 3 is 0.22, 0.26, and 0.17, respectively, but the absolute permeability is much smaller than that of sandstone Core a, whose porosity is 0.18. The porosity of hydrate Sample 4 and Sample 5 has reached a high level at about 0.3, but the permeability is only an average of 200 mD. When the porosity of conventional sandstone Core b is greater than 0.3, the permeability can reach more than 1000 mD. As for Sample 6, although the porosity is also approximately 0.2, its permeability is more than twice of that of Sample 1, 2, and 3. In summary, the permeability of the hydrate sample is lower than that of the conventional sandstone core and the distribution is complex, having poor correlation with the porosity. It is necessary to carry out further analysis for the microscopic pore structure of the sample by using fractal theory to find out the seepage mechanism.

4.2. Fractal Dimension and Permeability
The 3D fractal dimensions of hydrate samples and sandstone cores are calculated by the boxcounting method based on Equation (3). The range of 3D fractal dimension of eight samples and cores is 2.70~2.85. The fractal dimension characterizes a complex feature of the overall pore distribution, which is generally affected by many factors, especially porosity. The larger the porosity, the larger the fractal dimension, and the more complex the pore space distribution develops. From this perspective, in the case of low porosity, hydrate samples exhibit a high fractal dimension compared with conventional sandstone cores ( and ), which proves that the distribution of pore space of the hydrate sample is more complicated than that of the conventional sandstone cores (Table 2).

4.3. Lacunarity and Permeability
In order to observe lacunarity (heterogeneity) of different samples and cores more intuitively, medians of the normalized values at different scales are listed in Table 3. The range of lacunarity midvalue of eight samples and cores is 0.168~0.447. Among the hydrate samples with smaller porosity, it can be seen that though the porosities of Sample 1, 2, 3, and 6 are similar, permeability of Sample 6 is almost double of that of Sample 1, 2, and 3. The results indicate that with the heterogeneity of sample increases, the permeability decreases. The reason is that when porosity is small, permeability of the hydrate sample with uniform pore distribution depends on the connectivity of pore throats. That is, the stronger the pore distribution heterogeneity of the hydrate sample, the worse the connectivity of the pore throats, and the smaller the permeability in consequence. Therefore, for a hydrate sample with a small porosity, the smaller the calculated lacunarity (the weaker the heterogeneity), the easier it is to exhibit a larger permeability. That is, the stronger the pore distribution heterogeneity of the hydrate sample, the worse the pore throat connectivity, and the smaller the permeability consequently. Therefore, for a hydrate sample with small porosity, the smaller the calculated lacunarity (the weaker the heterogeneity), the more likely it is to possess higher permeability.

On the other hand, hydrate samples with larger porosity and sandstone cores show that the stronger the heterogeneity, the higher the permeability. Namely, Core a is more heterogeneous in the case of possessing smaller porosity, and its permeability is basically not lower than Samples 4 and 5; Sample 5 has a lower porosity than Sample 4, but it is highly heterogeneous, and the permeabilities in two directions are greater than Sample 4. The reason is that in the hydrate samples with larger porosity and sandstone cores, the distribution of pore space has already formed pore throats for fluid flow. Currently, the greater the heterogeneity of the pore distribution, the larger the diameter of the pore throats, which consequently leads to the higher permeability. Therefore, in hydrate samples with large porosity and all sandstone cores, the larger the calculated lacunarity (the stronger the heterogeneity), the greater the permeability.
4.4. Succolarity and Permeability
Table 4 shows the calculated value of the 3D succolarity on six directions of six hydrate samples. The range of the results is 0.0024~0.3312. It can be seen that, except for Sample 3, hydrate samples have little difference in succolarity values on all directions. Since succolarity can characterize the physical significance of the anisotropic characteristics of core pore structure distribution, the hydrate sample pore structure has no obvious anisotropic characteristics. For Sample 3, the extremely small succolarity value on direction also indicates that there is no connected pores in this direction, which explains why the permeability cannot be calculated.

Figure 7 is a schematic diagram showing the fitting of succolarity and permeability on different positive directions of six hydrate samples (except Sample 1 , Sample 3 , and Sample 6 , on which direction the permeability cannot be calculated). It can be seen that, succolarity values of the hydrate samples show a good correlation with the permeabilities (correlation coefficient is 0.663). The cause of the incomplete fitting is that the pore radius of Sample 5 is larger (Figure 5), which results in a higher permeability, however the effect of pore radius has not been considered by succolarity. Even so, compared with porosity, succolarity can be used to fit the permeability and establish a prediction model, which gives a certain theoretical guidance to the changing rule of hydrate reservoir permeability.
5. Conclusions
In this study, six natural gas hydrate samples are conducted with CT scanning to get the pore structure, which are used to calculate the porosity, permeability, and fractal parameters. The results show that the fractal parameter can quantitatively characterize the pore structure and analyze the change of permeability of the hydrate reservoir. The following conclusions can be derived: (a)Compared with the conventional sandstone cores, the hydrate samples show higher porosity and lower permeability, and the correlation between porosity and permeability is poor(b)From the comparison results of the fractal dimension, hydrate samples possess higher fractal dimension than conventional sandstone cores and the structural development is more complicated(c)From the results of lacunarity calculation, it turns out that for hydrate samples with smaller porosity, with the lacunarity of pore distribution (the weaker the heterogeneity) decreases, the permeability increases; for the hydrate samples with larger porosity and sandstone cores, with the lacunarity of pore distribution (the stronger the heterogeneity) increases, the permeability increases(d)According to calculation results of succolarity, hydrate samples basically show nonanisotropic development characteristics. It is clear that the correlation between the succolarity value and the permeability value is high. Therefore, a predictive model of permeability can be built to provide theoretical support for hydrate reservoir development
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This research was funded by the National Natural Science Foundation of China (No. 51991365), China Geological Survey Project (No. DD20190232), Key Special Project for Introduced Talents Team of Southern Marine Science and Engineering Guangdong Laboratory (Guangzhou) (GML2019ZD0105), and Key Program of Marine Economy Development (Six Marine Industries) Special Foundation of Department of Natural Resources of Guangdong Province [2020]047.
References
 E. D. Sloan Jr., “Fundamental principles and applications of natural gas hydrates,” Nature, vol. 426, no. 6964, pp. 353–359, 2003. View at: Publisher Site  Google Scholar
 J. Cai, Y. Xia, S. Xu, and H. Tian, “Advances in multiphase seepage characteristics of natural gas hydrate sediments,” Chinese Journal of Theoretical and Applied Mechanics, vol. 52, pp. 208–223, 2020. View at: Google Scholar
 S. Fan, C. Yu, X. Lang, Y. Wang, and J. Chen, “Micronanoscale studies on occurrence and gas production and storage technology of marine gas hydrates,” Earth Science, vol. 43, pp. 1542–1548, 2018. View at: Google Scholar
 X. Ge, J. Liu, Y. Fan, D. Xing, S. Deng, and J. Cai, “Laboratory investigation into the formation and dissociation process of gas hydrate by lowfield NMR technique,” Journal of Geophysical Research: Solid Earth, vol. 123, no. 5, pp. 3339–3346, 2018. View at: Publisher Site  Google Scholar
 M. R. Islam, “A new recovery technique for gas production from Alaskan gas hydrates,” Journal of Petroleum Science and Engineering, vol. 11, no. 4, pp. 267–281, 1994. View at: Publisher Site  Google Scholar
 K. A. Kvenvolden, “Gas hydrates—geological perspective and global change,” Reviews of Geophysics, vol. 31, no. 2, pp. 173–187, 1993. View at: Publisher Site  Google Scholar
 X. Liu, H. Liu, L. Xing, Y. Yin, and J. Wang, “Seismic lowfrequency shadow beneath gas hydrate in the Shenhu area based on the stereoscopic observation system,” Journal of Earth Science, vol. 29, no. 3, pp. 669–678, 2018. View at: Publisher Site  Google Scholar
 L. Zhang, C. Zhang, K. Zhang et al., “Porescale investigation of methane hydrate dissociation using the lattice Boltzmann method,” Water Resources Research, vol. 55, no. 11, pp. 8422–8444, 2019. View at: Publisher Site  Google Scholar
 J. Li, J. L. Ye, X. W. Qin et al., “The first offshore natural gas hydrate production test in South China Sea,” China Geology, vol. 1, no. 1, pp. 5–16, 2018. View at: Publisher Site  Google Scholar
 J. Cai, E. Perfect, C. L. Cheng, and X. Hu, “Generalized modeling of spontaneous imbibition based on HagenPoiseuille flow in tortuous capillaries with variably shaped apertures,” Langmuir, vol. 30, no. 18, pp. 5142–5151, 2014. View at: Publisher Site  Google Scholar
 J. Qiao, J. Zeng, S. Jiang, and Y. Wang, “Impacts of sedimentology and diagenesis on pore structure and reservoir quality in tight oil sandstone reservoirs: implications for macroscopic and microscopic heterogeneities,” Marine and Petroleum Geology, vol. 111, pp. 279–300, 2020. View at: Publisher Site  Google Scholar
 B. Hazra, D. A. Wood, V. Vishal, A. K. Varma, D. Sakha, and A. K. Singh, “Porosity controls and fractal disposition of organicrich Permian shales using lowpressure adsorption techniques,” Fuel, vol. 220, pp. 837–848, 2018. View at: Publisher Site  Google Scholar
 Y. Yao and D. Liu, “Comparison of lowfield NMR and mercury intrusion porosimetry in characterizing pore size distributions of coals,” Fuel, vol. 95, pp. 152–158, 2012. View at: Publisher Site  Google Scholar
 X. Zhao, M. J. Blunt, and J. Yao, “Porescale modeling: effects of wettability on waterflood oil recovery,” Journal of Petroleum Science and Engineering, vol. 71, no. 34, pp. 169–178, 2010. View at: Publisher Site  Google Scholar
 C. Zou, D. Dong, S. Wang et al., “Geological characteristics and resource potential of shale gas in China,” Petroleum Exploration and Development, vol. 37, no. 6, pp. 641–653, 2010. View at: Publisher Site  Google Scholar
 P. C. Carman, “Fluid flow through granular beds,” Transactions Institute of Chemical Engineers, vol. 15, pp. 150–167, 1997. View at: Google Scholar
 P. C. Carman, “The determination of the specific surfaces of powders I,” Journal of the Society of the Chemical Industrialists, vol. 57, pp. 225–234, 1938. View at: Google Scholar
 H. A. Nooruddin and M. E. Hossain, “Modified Kozeny–Carmen correlation for enhanced hydraulic flow unit characterization,” Journal of Petroleum Science and Engineering, vol. 80, no. 1, pp. 107–115, 2011. 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
 C. Manwart, U. Aaltosalmi, A. Koponen, R. Hilfer, and J. Timonen, “LatticeBoltzmann and finitedifference simulations for the permeability for threedimensional porous media,” Physical Review E, vol. 66, no. 1, pp. 16702–16713, 2002. View at: Publisher Site  Google Scholar
 L. Zhang, W. Jing, Y. Yang et al., “The investigation of permeability calculation using digital core simulation technology,” Energies, vol. 12, no. 17, article 3273, 2019. View at: Publisher Site  Google Scholar
 B. Yu and P. Cheng, “A fractal permeability model for bidispersed porous media,” International Journal of Heat and Mass Transfer, vol. 45, no. 14, pp. 2983–2993, 2002. View at: Publisher Site  Google Scholar
 J. Cai, Y. Xia, C. Lu, H. Bian, and S. Zou, “Creeping microstructure and fractal permeability model of natural gas hydrate reservoir,” Marine and Petroleum Geology, vol. 115, article 104282, 2020. View at: Publisher Site  Google Scholar
 C. Lu, Y. Xia, X. Sun et al., “Permeability evolution at various pressure gradients in natural gas hydrate reservoir at the Shenhu area in the South China Sea,” Energies, vol. 12, no. 19, article 3688, 2019. View at: Publisher Site  Google Scholar
 L. Liu, S. Dai, F. Ning, J. Cai, C. Liu, and N. Wu, “Fractal characteristics of unsaturated sands − implications to relative permeability in hydratebearing sediments,” Journal of Natural Gas Science and Engineering, vol. 66, pp. 11–17, 2019. View at: Publisher Site  Google Scholar
 M. N'Diaye, C. Degeratu, J. M. Bouler, and D. Chappard, “Biomaterial porosity determined by fractal dimensions, succolarity and lacunarity on microcomputed tomographic images,” Materials Science and Engineering: C, vol. 33, no. 4, pp. 2025–2030, 2013. View at: Publisher Site  Google Scholar
 Y. Xia, J. Cai, E. Perfect, W. Wei, Q. Zhang, and Q. Meng, “Fractal dimension, lacunarity and succolarity analyses on CT images of reservoir rocks for permeability prediction,” Journal of Hydrology, vol. 579, article 124198, 2019. View at: Publisher Site  Google Scholar
 H. Dong and M. J. Blunt, “Porenetwork extraction from microcomputerizedtomography images,” Physical Review E, vol. 80, no. 3, article 036307, 2009. View at: Publisher Site  Google Scholar
 C. Allain and M. Cloitre, “Characterizing the lacunarity of random and deterministic fractal sets,” Physical Review A, vol. 44, no. 6, pp. 3552–3558, 1991. View at: Publisher Site  Google Scholar
 A. Roy, E. Perfect, W. M. Dunne, N. Odling, and J.W. Kim, “Lacunarity analysis of fracture networks: evidence for scaledependent clustering,” Journal of Structural Geology, vol. 32, no. 10, pp. 1444–1449, 2010. View at: Publisher Site  Google Scholar
 R. H. C. de Melo and A. Conci, “How succolarity could be used as another fractal measure in image analysis,” Telecommunication Systems, vol. 52, no. 3, pp. 1643–1655, 2013. View at: Publisher Site  Google Scholar
 N. Ahmad, J. Rappaz, J.L. Desbiolles et al., “Numerical simulation of macrosegregation: a comparison between finite volume method and finite element method predictions and a confrontation with experiments,” Metallurgical and Materials Transactions A, vol. 29, no. 2, pp. 617–630, 1998. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2020 Hang Bian 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.