Research Article  Open Access
You Zhou, Songtao Wu, Zhiping Li, Rukai Zhu, Shuyun Xie, Xiufen Zhai, Lei Lei, "Permeability Calculation of Sand Conglomerate Reservoirs Based on Nuclear Magnetic Resonance (NMR)", Geofluids, vol. 2020, Article ID 1385469, 14 pages, 2020. https://doi.org/10.1155/2020/1385469
Permeability Calculation of Sand Conglomerate Reservoirs Based on Nuclear Magnetic Resonance (NMR)
Abstract
The concept of an intermingled fractal unit (IFU) model was first proposed by Atzeni and Pia in 2008, and their model has since been successfully applied to predict thermal conductivity, electrical conductivity, and the mechanical properties of porous media materials. This paper, based on the Pia IFU model, fits the pore size distribution spectrum to quantitatively characterize the Triassic Karamay Formation conglomerate reservoirs in the Mahu region, in the Junggar Basin of Northwest China, and makes permeability predictions using the free fluid spectrum according to the nuclear magnetic resonance (NMR) experimental data. The results show that the accuracy of the IFU model is significantly higher than that of the classic Coates and SDR models for conglomerate reservoirs with complex pore structures, indicating that this is an effective method to calculate permeability based on NMR. In addition, preliminary discussions are entered into regarding the intermingled fractal expression of the KozenyCarman equation and the relative permeability, in order to widen the application of the IFU model in reservoir physics. The derived expressions appear complicated in form but are straightforward to calculate and apply using computer programming since their iteration parameters are definite. The findings set out in this paper provide a valuable reference for further research of the IFU model in reservoir physics.
1. Introduction
The pore structure of reservoir rocks has fractal characteristics [1]. Because of this, a number of scholars have recently adopted fractal theory to establish the porosity or the permeability model of reservoir rocks and have achieved some productive results. For instance, Li et al. [2], Zheng et al. [3], and Bai et al. [4] have successfully applied fractal theory to build models for forecasting the permeability of tight sandstone reservoirs. With the development of fractal theory, however, it has become apparent that the theory has limitations in characterizing a real reservoir space. When computing fractal dimensions using capillary data, for example, it has been noted that in a loglog coordinate system, the to plot sometimes exhibits multifractal characteristics for some kinds of rocks, instead of being linearly fitted, which significantly influences the accuracy of computation [5–9]. Meanwhile, Sima et al. [10] have pointed out that sand conglomerate rocks with high heterogeneity tend to exhibit more evident multifractal characteristics or even nonfractal characteristics, on loglog coordinates, illustrating that the probability of multifractal characteristics is greater in the conglomerate samples. This section analyzes the fractal feature gained from the capillary buddle subsection model and draws similar conclusions for sand conglomerate rocks (Figure 1). Figures 1(a)–1(c) show that the fractal characteristics of the capillary pressure curve of the sand conglomerate reservoir generally need to be fitted with two or more straight lines, which indicates that different intervals of pore size have different fractal dimensions. Figure 1(d) shows that the fractal dimensions of some large intervals of pore size, calculated by the graphical method, are greater than 3, and this phenomenon is interpreted as being attributable to nonfractal characteristics.
(a) Well X003
(b) Well M152
(c) Well H3
(d) Well H3
To address this problem , Li et al. [11] and Lai et al. [12] proposed to run a weighted average of the mercury saturation of each fractal segment to obtain a weighted fractal dimension. However, this method neglects the detailed features of the pore space and hence only reflects the overall fractal characteristics of rock pores. Moreover, although the fractal dimension is a highly significant parameter in fractal theory, it has limitations in characterizing a real rock pore space. This is because porous media with different pore structures may have the same or similar fractal dimension [13]. Figure 2 shows that although the pore structures of the four samples are different, their fractal dimensions tend to be similar. The specific calculation process for threedimensional fractal dimensions is set out by Zhang et al. [15]. Therefore, utilizing a single fractal dimension to establish a porosity or permeability model might inevitably constrain subsequent research on permeation performance and reservoir engineering.
(a)
(b)
(c)
(d)
Given the limitations of the single fractal dimension, Italian scholar Pia et al. [16–23] introduced an intermingled fractal unit (IFU) model, which has proven to be useful in simulating porous media. Pia and his research team applied the model to predict the thermal conductivity, electrical conductivity, mechanical properties, and permeability of porous media materials. Intuitively, the IFU model can be considered to be several superimposed fractal units. Each of the small fractal units in the IFU model is a modified Sierpiński carpet, with each fractal unit having different iteration parameters. Thus, the IFU model can be regarded as the Sierpiński Carpet modified by altering its iteration rule. The advantage of this model is that it can present detailed characterizations of internal spaces of complicated objects.
At present, applications of the IFU model have relatively little engagement in the field of petroleum engineering. Since Pia first used the IFU model to predict permeability based on mercury intrusion capillary pressure (MICP) data [20], some scholars have also done a number of works to widen the application of the model in recent years. Tan et al. [24] and Chen et al. [25] applied the IFU model to forecast the permeability with consideration of tortuosity according to the MICP data. Tan et al. [26] used Xray computed tomography (CT) and scanning electron microscopy (SEM) data to make permeability prediction based on the model. Li et al. [27, 28] applied a 2D and 3D IFU model to predict permeability using SEM images for shales. Objectively, at present, the applications of the IFU model in reservoir physics basically focus on the prediction of permeability, and the IFU model based on NMR data has not been reported.
This paper, based on the Pia IFU model, constructs a model to fit the pore size distribution spectrum according to NMR measurements and makes permeability predictions using the free fluid spectrum. Meanwhile, some key issues in the process of modeling using NMR data are also pointed out. Also, new intermingled fractal expressions of the KozenyCarman equation and relative permeability are derived to broaden the application of the IFU model in reservoir physics.
2. Methodology
2.1. Introduction to the Pia IFU Model
This section sets out some core steps for the derivation process of iterative formulas for the IFU model. For a detailed explanation of the entire process, see Pia et al. [20].
Without derogation from general principles, a specific IFU is used as an example to examine the relationships between its iteration parameters (Figure 3).
Let be the maximum pore side length of this fractal unit. The side length of the fractal unit will then be . is the number of the child squares removed from each parent square in every iteration, and is the ratio factor (a multiple of equal divisions of the side length in each iteration). Thus, and . is the number of the child squares that do not take part in the iteration. Thus, . The value is the iteration time.
According to fractal geometric theory [29], the fractal dimension is the value of the number of squares that participate in the next iteration divided by the ratio factor :
For the Sierpiński carpet fractal model, the fractal dimension of the 3D spatial model is the value of the 2D plane model plus 1.
To produce a Sierpiński carpet similar to the unit shown in Figure 3, the side length of the new pores formed in the th iteration is defined by Pia et al. [20] as follows:
The number of new pores formed in the th iteration is given by Li et al. [27]:
For ease of theoretical calculation, it is assumed that the real percolation section of every square pore is its maximum incircle and that the diameter, position, and quantity of pores will be the same on any rock section when the 2D IFU model is extended to 3D. If is the length of the model and the total iteration time for the fractal unit is , then the generated total pore volume is given as:
The IFU model can then be utilized to reproduce the pore size distribution spectrum and porosity. The procedure used to reproduce porous microstructures is shown in Figure 4. For a detailed explanation of the process, see Pia’s series of papers.
2.2. Permeability Calculation Using Traditional Classic Models Based on NMR Experiments
NMR experiments can quantitatively characterize the pore structures of rocks and can be applied directly to calculate permeability. The current classic models for calculating permeability using NMR data are the TimurCoates and SDR models.
The TimurCoates model [30] is expressed as follows:
The SDR model [31] can be expressed as follows: where is the permeability (10^{15} m^{2}), is the effective porosity (%), is the geometric mean of (ms), FFI is the saturation of movable fluid (%), BVI is the saturation of immovable fluid (%), and , , , , , and are the corresponding parameters of each model, which are determined by the results of actual core data used to fit experimental permeability. In the Coates model, in order to facilitate the calculation, the values of and are sometimes directly regarded as 4 and 2, respectively, as well as the values of and in the SRD model [32, 33].
In recent years, some scholars have pointed out that the classic TimurCoates and SDR models have certain limitations for calculating the permeability of complex reservoirs. For example, the classic model has good applicability for the calculation of permeability in medium and high permeability reservoirs. However, for reservoirs with low porosity and low permeability, it is difficult to correlate the calculation results with experimental data [34]. This paper applied the classical model to predict the permeability of sand conglomerate reservoirs in the Mahu rim region of the Junggar Basin, and the calculation results exhibit significant discrepancies from the experimental results (see Section 3.2). This is because the sand conglomerate rocks have extremely high heterogeneity due to their complex gravelcement configuration. By using microCT scanning technology on some typical sand conglomerate rock samples, it is indicated that the sand conglomerate rocks in the research area have complex and diverse microscopic pore structures (Figure 5).
Meanwhile, the actual calculations in the classic models are complicated, resulting in large volumes of calculations for the process of converting the NMR spectrum into a pseudocapillary pressure curve using the nonlinear method. Moreover, it is found that the conversion coefficients and conversion indices calculated by the classic model vary greatly, even for rock samples in the same study area. The average value of the region does not satisfy the requirement for calculation accuracy, and the nonlinear fitting method does not explain the theoretical relationship between spectral relaxation time and real pore sizes. Because of this, this paper utilizes the “centrifugation method” [35] to convert the T_{2} spectrum into a corresponding pore size distribution for subsequent permeability prediction based on the IFU model.
There is an issue worth pointing out here. Since the sand conglomerates are easy to split in the process of centrifugal experiment, the cores are wrapped with a special heatshrinkable film in the experiments. Meanwhile, in order to determine the centrifugal force suitable for the samples in the study area, five representative samples were selected to test six different centrifugal forces of 100 psi, 150 psi, 200 psi, 250 psi, 300 psi, and 350 psi. By analyzing the descend range of irreducible water saturation before and after centrifugation, it was found that until the centrifugal force increased from 300 psi to 350 psi, the reduction of irreducible water saturation of five samples with different physical properties tended to be stable. Therefore, we select a centrifugal force of 300 psi to calibrate a cutoff value, and then, the corresponding lower limit of the movable fluid throat radius is 0.0704 μm. The experiment adopts a GeoSpec2 NMR core analyzer produced by Oxford Instruments Magnetic Resonance in England. The working frequency is 2 MHz, the temperature of the instrument is 35°C, the room temperature is 25°C, and the humidity is 50%60%.
2.3. Permeability Calculation Using the IFU Model Based on NMR Experiments
The expression of permeability based on the IFU model is derived in this section.
The flow of a single pore formed after the th iteration is given by the HagenPoiseuille equation: where is the flow (cm^{3}/s), is the pressure difference between the ends of the capillary tube (atm), is the viscosity of fluid (mPa·s), is the characteristic length of the capillary tube (cm), and is the tortuosity (dimensionless parameter). The ratio of the capillary tube’s real length to the characteristic length is expressed:
Based on the average tortuosityporosity analytical model for porous media consisting of square particles introduced by Yu et al. [36], Yun et al. [37] proposed an enhanced expression of the average tortuosity of porous media consisting of circular particles, expressed as follows: where is the porosity. Thus, the total flow () of the pores formed from the first to the th iteration is given as:
According to the equivalent percolation principle,
Finally, the expression for computing permeability can be derived as follows:
Thus, the permeability of a single fractal unit is given in formula (12), and the formula is similar in form to the expression in Pia et al. [20]. If the IFU model contains several fractal units, formulas (9) and (10) can be used to calculate its tortuosity and flow rates, respectively. The total permeability of the IFU model is easily obtained according to the equivalent seepage principle.
3. Results and Discussion
3.1. Conversion of the Spectrum to Pore Size Distribution
The basic principle of the “centrifugation method” is that the pore size , corresponding to the cutoff value determined by the centrifugal experiment, is a constant, and thence, the pore size , corresponding to the th relaxation time , is given as [38]
During the NMR experiment, the centrifugal pressure for testing is 300 psi, so the calculated pore size is 0.0704 μm, corresponding to . By applying formula (13), the pore size corresponding to different relaxation times can be calculated, and the pore size distribution based on spectra can then be constructed as shown in Figure 6.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
The spectrum of 100% saturated water is the reflection of all pore spaces in the rocks [39]; i.e., it contains connected pores, dead pores, and micropores storing bound water. Therefore, for permeability calculations based on NMR data, the free fluid spectrum should be used rather than the 100% saturated water spectrum. This is because only connected pores with percolation ability contribute to permeability in the pore space, and thus, only the free fluid spectrum has good correspondence with the real seepage space. The acquisition method for this value is to subtract the centrifugal spectrum signal from the saturated water spectrum signal, and the remainder is the free fluid spectrum (Figure 7).
(a)
(b)
Also, an abnormality persists for some samples, in that the centrifugal spectrum in some samples is higher than the saturation spectrum (Figures 8(a) and 8(b)), resulting in the appearance of a “negative accumulation area” when the free fluid accumulation curve is obtained (Figures 8(c) and 8(d)). These anomalies in the micropore intervals are difficult to avoid, because they are caused by system errors in the instrumentation. Therefore, when using the IFU model to fit the pore size distribution spectrum of the free fluid, it is suggested that the model is fitted to the highest point of the cumulative curve to circumvent system errors.
(a)
(b)
(c)
(d)
To summarize, the procedures of conversion of the spectrum to pore size distribution are as follows: (1)Subtract the centrifugal spectrum from the 100% saturated water spectrum to obtain the free fluid spectrum(2)Eliminate the “negative accumulation area” to obtain an effective free fluid accumulation curve(3)Use formula (13) to convert the effective free fluid spectrum to the pore diameter cumulative distribution curve
3.2. Permeability Prediction
The steps for computation of permeability with the IFU model can be summarized as follows: (1)Select the fractal units with different parameters for building the IFU model (Table 1), and change the iteration parameters to assure a consistency between the pore diameter distribution spectra of the structural model and the experimental data. Figure 9 reveals fitting results for some representative samples(2)Based on the completion of the above fitting process, formula (12) can be applied in the IFU model to calculate permeability (Table 2)

(a)
(b)
(c)
(d)

In order to make a better comparison between the computed results of the IFU model and the classical models, we subdivide the two classical models into the following four types according to whether the undetermined parameters of the models are constant or not. The relevant parameters of the experiment and models are shown in Tables 3 and 4.


The parameters in Table 3 are substituted into each model in Table 4, and the best model coefficient is determined by programming using the method of threecoefficient regression analysis. It can be seen that in Coates model 1 and SDR model 1, except that the coefficient is a variable, the index is also not a constant. Actually, when the regression statistical method is used to determine the coefficients and exponents in the equation, due to the lack of data points, the variation range of model parameters is closely related to the selection of initial values, which requires some prior knowledge. Therefore, a more feasible method is to take logarithms on both sides of the equation and then calculate each undetermined coefficient by regression statistical analysis.
Table 5 and Figure 10 show that the classical models cannot accurately calculate the permeability of the sand conglomerate reservoirs in the study area. The average relative error of the results calculated using the IFU model is 23.1%, within a certain tolerance for error, indicating that the IFU model is reliable in simulating the movable fluid pore space based on NMR data.

Also, it is noted that the permeability calculation results, based on the IFU model, are mostly smaller than the experimental ones. This is mainly due to the fact that, for reservoirs with complex pore structures such as sand conglomerate reservoirs, it is difficult to replicate the real conditions of a bound water state in centrifugal experiments. Centrifugal parameter settings often have strong regional characteristics. However, in fact, in order to truly replicate a bound water state and obtain the free fluid spectrum, the centrifugal force and centrifugation time must be markedly different for each sample.
Therefore, it is not realistic to increase the centrifugal force and the centrifugation time without reference to other parameters, because “grain shattering” occurs easily during the centrifugation procedure for sand conglomerate rocks. Excessive centrifugal force will damage the samples and then affect the experimental results. Therefore, since it is difficult to accurately measure the content of the movable fluid in centrifugal experiments, the computational accuracy of the IFU model based on NMR data is inevitably somewhat compromised.
3.3. Other Applications of the IFU Model in Reservoir Physics
3.3.1. Intermingled Fractal Expression of the KozenyCarman Equation
The KozenyCarman equation is a semiempirical formula widely used in calculating permeability, whereas the determination of the KC constant is highly empirical. In fact, the KC constant is closely related to the microscopic pore structure and is not a fixed value. Since the IFU model can accurately quantify specific pore throat information in the porous medium, the KC constant, with a practical physical manifestation, can be obtained using the IFU.
In order to facilitate the presentation of the formula derivation process, this section takes the example of an IFU model with only one fractal unit. Taking tortuosity into consideration, the total internal surface area () of the pores produced after the th iteration is straightforwardly obtained according to equation (3):
Thus, the specific surface () of the IFU model is
From equation (4), the porosity () of the IFU model can be rendered as where is the crosssectional area of the model. Note that for some tight reservoirs with extralow porosity and permeability, the pore sizes converted using the NMR free fluid spectrum are the effective flow radius, without considering boundary layer effects. However, for pore sizes converted by mercury injection capillary pressure (MICP) measurements, it is necessary to deduct the thickness of the boundary layer to obtain an effective flow radius, since mercury does not produce a boundary layer in MICP experiments. For quantitative calculation of the thickness of the boundary layer, please refer to the references made by Tian et al. [40] and Meng et al. [41]. Thus, permeability of the IFU model, taking the thickness of the boundary layer () into consideration, is given by
According to the classical KozenyCarman equation, where is the KC constant. By substituting formulas (15)–(17) in (18), the expression of the KC constant is obtained:
It can be seen that the KC constant is not an empirical constant without precise physical meaning but is controlled by factors such as iteration parameters, tortuosity, and the thickness of the boundary layer in the IFU model. Once these parameters are given, the value of KC for each sample can be calculated simply by programming.
3.3.2. Intermingled Fractal Expression of Relative Permeability
Based on the IFU model, the expression of relative permeability can also be derived. Similarly to Section 3.3.1, the IFU model with only one fractal unit is taken as an example in this case. The assumptions for the model are as follows: (a)The reservoir is composed of unequal diameter capillaries with fractal features(b)The fluid flow conforms to Darcy’s law, and the wet phase is distributed in the capillaries whose pore diameters are smaller than the critical capillary radius () and the nonwet phase is distributed in the capillaries whose pore diameters are larger than the critical capillary radius. Therefore, only one fluid exists in a single capillary(c)The boundary layer fluid adheres to the inside of the tube wall in the form of a bound water film(d)The fluid is incompressible, the viscosity is constant, and the percolation process is isothermal
Thus, considering the thickness of the boundary layer, formula (12) can be rewritten as
According to the definition of saturation,
The effective permeability can then be rewritten as
It is indicated from the previous study that, as the iteration time () increases, the diameters of the generated pores reduce. Therefore, if is the critical iteration time, then the nonwet phase exists in the pores of and the wet phase exists in the pores of .
The effective permeability of the wet phase is then given by
The effective permeability of the nonwet phase is given by
The relative permeability of the wet phase can then be derived as follows:
And the relative permeability of the nonwet phase is where , , and are the tortuosities of 100% saturated singlephase, wetphase, and nonwetphase fluids, respectively, and they are determined by formula (9) according to their respective saturations. Also, it is found that the forms of the formulas are similar to those expressed in the classical Burdine model [42].
Although these derived formulas seem complicated in form, in fact, the IFU model can give detailed quantitative information about each pore throat in the porous medium by fitting the pore size distribution spectrum, which makes the calculations of these parameters (permeability, KC constant, relative permeability, etc.) rapid and convenient by computer programming.
4. Conclusions
(1)The classical Coates and SDR models cannot fulfill the precision requirements of permeability calculations for sand conglomerate reservoirs in the study area, and the calculation accuracy of results based on the IFU model is significantly improved. Also, the modelcomputed results are slightly smaller than the experimental ones, because it is difficult to accurately measure the content of the movable fluid in the centrifugal experiment(2)The new intermingled fractal expressions of the KozenyCarman equation and relative permeability were derived based on the IFU model, and the physical meanings of the parameters were explained. The expressions appear complicated in form but are easy to calculate and apply using computer programming, since their iteration parameters are defined(3)The authors believe that the application of the IFU model in reservoir physics is much more than that, and further exploration of possible improvements and applications of the model is necessary and meaningful
Data Availability
All data generated or analyzed during this study are included in this article, and all data included in this study are available upon request by contact with the corresponding author.
Conflicts of Interest
The authors declare that there is no conflict of interest regarding the publication of this paper.
Acknowledgments
This research was financially supported by the National Key Basic Research and Development Program (973 Program), China (2014CB239000), the China National Science and Technology Major Projects (2017ZX05001), and the PetroChina Science and Technology Major Projects (2016B0304).
References
 C. Z. He and M. Q. Hua, “Fractal geometry description of reservoir pore structure,” Oil and Gas Geology, vol. 19, no. 1, pp. 15–23, 1998. View at: Google Scholar
 L. R. Li, S. Y. Yuan, and Y. L. Hu, “A new model for describing the relationship between the permeability and the porosity of fractal porous media,” Journal of Xi'an Shiyou University (Natural Science Edition), vol. 25, no. 3, pp. 49–51, 2010. View at: Google Scholar
 B. Zheng and J. H. Li, “A new fractal permeability model for porous media based on KozenyCarman equation,” Natural Gas Geoscience, vol. 26, no. 1, pp. 193–198, 2015. View at: Google Scholar
 R. T. Bai, Z. P. Li, J. X. Nan, F. P. Lai, H. Li, and Q. Wei, “The fractal permeability model in tight sand reservoir accounts for startup gradient,” Natural Gas Geoscience, vol. 27, no. 1, pp. 142–148, 2016. View at: Google Scholar
 K. Li, “Analytical derivation of brooks–Corey type capillary pressure models using fractal geometry and evaluation of rock heterogeneity,” Journal of Petroleum Science and Engineering, vol. 73, no. 12, pp. 20–26, 2010. View at: Publisher Site  Google Scholar
 H. Y. Liu, Z. Y. Tian, and Z. Y. Xu, “Quantitative evaluation of carbonate reservoir pore structure based on fractal characteristics,” Lithologic Reservoirs, vol. 29, no. 5, pp. 97–105, 2017. View at: Google Scholar
 F. Wang, L. Jiao, Z. Liu, X. Tan, C. Wang, and J. Gao, “Fractal analysis of pore structures in low permeability sandstones using mercury intrusion porosimetry,” Journal of Porous Media, vol. 21, no. 11, pp. 1097–1119, 2018. View at: Publisher Site  Google Scholar
 F. Wang, K. Yang, and J. Cai, “Fractal characterization of tight oil reservoir pore structure using nuclear magnetic resonance and mercury intrusion porosimetry,” Fractals, vol. 26, no. 2, article 1840017, 2018. View at: Publisher Site  Google Scholar
 Q. H. Zhao, Z. P. Li, F. P. Lai, and Y. Du, “Fractal characteristics of rock pores in shale gas reservoir,” Journal of Liaoning Technical University (Natural Science Edition), vol. 35, no. 9, pp. 914–919, 2016. View at: Google Scholar
 L. Q. Sima, G. D. Yang, F. Wu, L. Wang, and F. Meng, “Fractal feature about the pore structure and controlling factor in tight glutenite reservoir in Baikouquan formation of Mahu depression in Junggar Basin,” Well Logging Technology, vol. 40, no. 5, pp. 609–616, 2016. View at: Google Scholar
 D. Y. Li, S. B. Zang, X. J. Ren, H. J. Hui, H. Wang, and F. F. Miao, “Study on pore structure of low permeability reservoirs with fractal theory,” Liaoning Chemical Industry, vol. 39, no. 7, pp. 723–726, 2010. View at: Google Scholar
 J. Lai, G. W. Wang, Y. Q. Zheng, W. Y. Li, and C. Cai, “Method for calculating the fractal dimension of the pore structure of low permeability reservoirs: a case study on the Xujiahe formation reservoir in Central Sichuan basin,” Journal of Northeast Petroleum University, vol. 37, no. 1, pp. 1–7, 2013. View at: Google Scholar
 J. Z. Zhang, Fractal (Second Edition), Tsinghua University Press, Beijing, 2011.
 Y. Zhou, S. Wu, Z. Li et al., “Multifractal study of threedimensional pore structure of sandconglomerate reservoir based on CT images,” Energy & Fuels, vol. 32, no. 4, pp. 4797–4807, 2018. View at: Publisher Site  Google Scholar
 T. F. Zhang, S. Y. Xie, Z. Y. Bao, Y. Wang, and L. Y. Pan, “Fractal and multifractal research on pore system for porous dolomite reservoirs based on highresolution CT,” Geological Science and Technology Information, vol. 35, no. 6, pp. 55–62, 2016. View at: Google Scholar
 G. Pia, L. Casnedi, M. Ionta, and U. Sanna, “On the elastic deformation properties of porous ceramic materials obtained by poreforming agent method,” Ceramics International, vol. 41, no. 9, pp. 11097–11105, 2015. View at: Publisher Site  Google Scholar
 G. Pia, L. Casnedi, and U. Sanna, “Porosity and pore size distribution influence on thermal conductivity of yttriastabilized zirconia: experimental findings and model predictions,” Ceramics International, vol. 42, no. 5, pp. 5802–5809, 2016. View at: Publisher Site  Google Scholar
 G. Pia and L. Casnedi, “Heat transfer in high porous alumina: experimental data interpretation by different modelling approaches,” Ceramics International, vol. 43, no. 12, pp. 9184–9190, 2017. View at: Publisher Site  Google Scholar
 G. Pia and U. Sanna, “A geometrical fractal model for the porosity and thermal conductivity of insulating concrete,” Construction and Building Materials, vol. 44, pp. 551–556, 2013. View at: Publisher Site  Google Scholar
 G. Pia and U. Sanna, “An intermingled fractal units model to evaluate pore size distribution influence on thermal conductivity values in porous materials,” Applied Thermal Engineering, vol. 65, no. 12, pp. 330–336, 2014. View at: Publisher Site  Google Scholar
 G. Pia and U. Sanna, “An intermingled fractal units model and method to predict permeability in porous rock,” International Journal of Engineering Science, vol. 75, pp. 31–39, 2014. View at: Publisher Site  Google Scholar
 G. Pia, C. Siligardi, L. Casnedi, and U. Sanna, “Pore size distribution and porosity influence on sorptivity of ceramic tiles: from experimental data to fractal modelling,” Ceramics International, vol. 42, no. 8, pp. 9583–9590, 2016. View at: Publisher Site  Google Scholar
 G. Pia, “High porous yttriastabilized zirconia with aligned pore channels: morphology directionality influence on heat transfer,” Ceramics International, vol. 42, no. 10, pp. 11674–11681, 2016. View at: Publisher Site  Google Scholar
 X. H. Tan, J. Y. Liu, X. P. Li, L. H. Zhang, and J. Cai, “A simulation method for permeability of porous media based on multiple fractal model,” International Journal of Engineering Science, vol. 95, pp. 76–84, 2015. View at: Publisher Site  Google Scholar
 Z. L. Chen, N. T. Wang, L. Sun, X. H. Tan, and S. Deng, “Prediction method for permeability of porous media with tortuosity effect based on an intermingled fractal units model,” International Journal of Engineering Science, vol. 121, pp. 83–90, 2017. View at: Publisher Site  Google Scholar
 X.H. Tan, X.P. Li, J.Y. Liu, L.H. Zhang, and J. Cai, “Fractal analysis of stress sensitivity of permeability in porous media,” Fractals, vol. 23, no. 2, article 1550001, 2015. View at: Publisher Site  Google Scholar
 C. Li, M. Lin, L. Ji, W. Jiang, and G. Cao, “Investigation of intermingled fractal model for organicrich shale,” Energy & Fuels, vol. 31, no. 9, pp. 8896–8909, 2017. View at: Publisher Site  Google Scholar
 C. Li, M. Lin, L. Ji, W. Jiang, and G. Cao, “Rapid evaluation of the permeability of organicrich shale using the 3D intermingledfractal model,” SPE Journal, vol. 23, no. 6, pp. 2175–2187, 2018. View at: Publisher Site  Google Scholar
 B. B. Mandelbrot, The Fractal Geometry of Nature, Freeman, New York, 1983.
 G. R. Coates, J. Galford, and D. Mardon, “A new characterization of bulkvolume irreducible using magnetic resonance,” Logic and Analysis, vol. 39, no. 1, pp. 51–63, 1998. View at: Google Scholar
 W. E. Kenyon, P. I. Day, C. Straley, and J. F. Willemsen, “A threepart study of NMR longitudinal relaxation properties of watersaturated sandstones,” SPE Formation Evaluation, vol. 3, no. 3, pp. 622–636, 2013. View at: Publisher Site  Google Scholar
 Z. G. Cheng, S. C. Luo, Z. W. Du, S. Chang, G. L. Li, and H. Li, “The method to calculate tight sandstone reservoir permeability using pore throat characteristic parameters,” Well Logging Technology, vol. 38, no. 2, pp. 185–189, 2014. View at: Google Scholar
 L. Z. Zou and Q. H. Liu, “Analysis of permeability model of nuclear magnetic resonance logging,” World Well Logging Technology, vol. 183, no. 3, pp. 27–31, 2011. View at: Google Scholar
 M. S. Zhou, J. Yin, J. L. Pan, S. L. Ren, S. X. Yang, and B. Z. Liu, “Method for effectiveness evaluation of low permeability reservoir based on NMR logging pore,” World Well Logging Technology, vol. 38, no. 4, pp. 452–457, 2014. View at: Google Scholar
 Y. Yao, D. Liu, Y. Che, D. Tang, S. Tang, and W. Huang, “Petrophysical characterization of coals by lowfield nuclear magnetic resonance (NMR),” Fuel, vol. 89, no. 7, pp. 1371–1380, 2010. View at: Publisher Site  Google Scholar
 Y. BoMing and L. JianHua, “A geometry model for tortuosity of flow path in porous media,” Chinese Physics Letters, vol. 21, no. 8, pp. 1569–1571, 2004. View at: Publisher Site  Google Scholar
 M. J. Yun, M. B. Yu, P. Xu, and J. S. Wu, “Geometrical models for tortuosity of streamlines in threedimensional porous media,” Canadian Journal of Chemical Engineering, vol. 84, no. 3, pp. 301–309, 2006. View at: Google Scholar
 Y. B. Yao, D. M. Liu, Y. D. Cai, and J. Q. Li, “Fine quantitative characterization of pores and fractures of coal based on NMR and XCT,” Scientia Sinica Terrae, vol. 40, no. 11, pp. 1598–1607, 2010. View at: Google Scholar
 C. X. Ning, Z. X. Jiang, Z. Y. Gao et al., “Quantitative evaluation of pore connectivity with nuclear magnetic resonance and high pressure mercury injection: a case study of the lower section of Es3 in Zhanhua sag,” Journal of China University of Mining and Technology, vol. 46, no. 3, pp. 578–585, 2017. View at: Google Scholar
 X. F. Tian, L. S. Cheng, X. L. Li et al., “A new method to calculate relative permeability considering the effect of porethroat distribution,” Journal of Shaanxi University of Science & Technology (Natural Science Edition), vol. 6, pp. 100–104, 2014. View at: Google Scholar
 Q. Meng, H. Liu, and J. Wang, “A critical review on fundamental mechanisms of spontaneous imbibition and the impact of boundary condition, fluid viscosity and wettability,” Advances in GeoEnergy Research, vol. 1, no. 1, pp. 1–17, 2017. View at: Publisher Site  Google Scholar
 N. T. Burdine, “Relative permeability calculation from size distribution data,” Transactions of AIME, vol. 198, pp. 71–78, 1953. View at: Google Scholar
Copyright
Copyright © 2020 You Zhou 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.