Permeability Calculation of Sand Conglomerate Reservoirs Based on Nuclear Magnetic Resonance (NMR)
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 Kozeny-Carman 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.
The pore structure of reservoir rocks has fractal characteristics . 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. , Zheng et al. , and Bai et al.  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 log-log 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.  have pointed out that sand conglomerate rocks with high heterogeneity tend to exhibit more evident multifractal characteristics or even nonfractal characteristics, on log-log 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.  and Lai et al.  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 . 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 three-dimensional fractal dimensions is set out by Zhang et al. . Therefore, utilizing a single fractal dimension to establish a porosity or permeability model might inevitably constrain subsequent research on permeation performance and reservoir engineering.
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 , some scholars have also done a number of works to widen the application of the model in recent years. Tan et al.  and Chen et al.  applied the IFU model to forecast the permeability with consideration of tortuosity according to the MICP data. Tan et al.  used X-ray 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 Kozeny-Carman equation and relative permeability are derived to broaden the application of the IFU model in reservoir physics.
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. .
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 , 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.
The number of new pores formed in the -th iteration is given by Li et al. :
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 Timur-Coates and SDR models.
The Timur-Coates model  is expressed as follows:
The SDR model  can be expressed as follows: where is the permeability (10-15 m2), 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 Timur-Coates 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 . 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 gravel-cement configuration. By using micro-CT 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”  to convert the T2 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 heat-shrinkable 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 Hagen-Poiseuille equation: where is the flow (cm3/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 tortuosity-porosity analytical model for porous media consisting of square particles introduced by Yu et al. , Yun et al.  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. . 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 
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.
The spectrum of 100% saturated water is the reflection of all pore spaces in the rocks ; 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).
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.
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)
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 three-coefficient 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 Kozeny-Carman Equation
The Kozeny-Carman 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 cross-sectional 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.  and Meng et al. . Thus, permeability of the IFU model, taking the thickness of the boundary layer () into consideration, is given by
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 single-phase, wet-phase, and non-wet-phase 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 .
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.
(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 model-computed 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 Kozeny-Carman 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
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.
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 (2016B-0304).
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 Kozeny-Carman 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 start-up gradient,” Natural Gas Geoscience, vol. 27, no. 1, pp. 142–148, 2016.View at: 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
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.
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 high-resolution CT,” Geological Science and Technology Information, vol. 35, no. 6, pp. 55–62, 2016.View at: 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 bulk-volume irreducible using magnetic resonance,” Logic and Analysis, vol. 39, no. 1, pp. 51–63, 1998.View at: 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
M. J. Yun, M. B. Yu, P. Xu, and J. S. Wu, “Geometrical models for tortuosity of streamlines in three-dimensional 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 X-CT,” 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 pore-throat distribution,” Journal of Shaanxi University of Science & Technology (Natural Science Edition), vol. 6, pp. 100–104, 2014.View at: 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