Determination of Specific Surface of Rock Grains by 2 D Imaging

Specific surface is an important parameter for predicting permeability of porous rocks. Many digital methods have been invented to extract the rock properties via imaging such as Micro-CT. With utilizing 3D volume data, this helps in precise investigation; however, it is neither economically efficient nor can be applied for different situations. In this study, a new approach is developed to estimate rock specific surface using 2D thin section images with micron resolution. One specific conclusion of this study is that there is specific ratio between the specific perimeter of 2D images and the specific surface in the 3D real rock structure. To further investigate this ratio several 3D blocks of rock volume data have been virtually cut in every possible angle and the value of specific perimeter calculated for each obtained 2D thin section. Finally, the predicted value of specific surface for 6 rock types is compared with the real values calculated from the original 3D data. Result indicates acceptable precision of this approach for sandstone rocks whereas not applicable for carbonate rocks.


Introduction
Specific surface area of a porous material is defined as the interstitial surface area of pores per unit of bulk volume or solid volume of the porous material.Specific surface is critical parameter considering fluid permeability of a porous material [1][2][3].Since specific surface is the ratio of area to volume, it has dimension of L −1 .The specific surface of porous granular media is defined as follows [3]: Surface area of grains Volume of rock grains .
Several methods for estimating specific surface of porous materials have been discussed in the literature until now.Collins [1] is one of the first pioneers of using these methods for specific surface measurement.Collins categorized the methods of specific surface determination into three main groups as follows: (A) adsorption measurement of specific chemicals into pore walls, (B) measuring the fluid flow through medium [4], and (C) statistical methods such as stereological, optical, and NMR analysis.
In the adsorption based methods, Methylene Blue Dye, as well as gas adsorption, has been used to determine the surface area of clay minerals for several decades [5].Dogan et al. [6] measured the specific surface of clay minerals by the Brunauer, Emmett, and Teller (BET) method of adsorption of nitrogen gas.Fluid flow methods which are basically an application of Kozeny-Carman equation are widely treated in the literature.This equation is the most famous permeabilityporosity relation, which is widely used in the field of flow in porous media and is the starting point for many other permeability models [7,8].The equation can be written as follows [4]: where  is permeability,  0 is a constant depending on pore shapes,  is porosity, and  is the specific surface area of porous media.Nuclear magnetic resonance (NMR) methods are also utilized widely to characterize rock's pore geometries.Pape et al. [9] developed a fractal analysis method for estimation of specific surface area of rocks using NMR data.The optical method which is based on statistical image processing has been developed by Chalkley [10].The method is based on random tests on two-dimensional images from cross sections of synthetic porous medium [10].In this regard a specified needle with length of  is randomly thrown on the surface of cross section and the number of pores and grains on which the needle is laying should be counted [10].The modified equation which estimates the specific surface of particles using Chalkley's method can be expressed as follows [1]: where  is the number of particle perimeters intersected by the needle, ℎ is the number of pore spaces recumbent under the needle,  is the needle length, and  is the porosity of porous medium.Debye et al. [11] presented a method for estimation of specific surface of heterogeneous solids by measuring the X-ray scattering pattern from different materials.They tried to correlate the porosity of porous structure versus the specific surface of its particles and verified the results by measuring the specific surface using adsorption methods.Berryman [12] modified the method used by Debye et al. [11] by angular averaging between normal vectors of rock grain surface in order to estimate the specific surface value using scattering data of X-ray.Berryman and Blair [13] utilized scanning electron microscopy (SEM) to estimate porous medium porosity and specific surface by applying two-point correlation functions.They compared the results obtained from the Kozeny-Carman equation and experimental permeability measurement to adjust the two-point correlation function for finding the specific surface.Berryman and Blair [14] extended their work to consider tortuosity and electrical formation factor to calculate specific surface using data available for Berea sand stone.Furthermore, they discussed how to choose the optimum magnification of imaging in order to estimate the specific surface from rock cross sections.
Yu et al. [15] developed a fractal two-phase porous media model and estimated the specific surface areas in two and three dimensions.They expressed that the specific surface of porous media is a function of the fractal dimensions and microstructural parameters.In addition, they found that the 2D specific surface (specific perimeter) is in relation to the 3D specific surface, which is dependent on the grain size and fractal dimension.
X-ray microtomographic imaging and visualization of rock material at the pore scale can give an important insight for understanding rock structure's properties.3D images allow one to map in detail the pore and grain structure and interconnectivity of the porous medium [16].Okabe and Blunt [17] utilized multiple-point statistics based on twodimensional thin sections as training images to generate geologically realistic 3D pore-space representations.As they stated, the thin section images can provide multiple-point statistics, which describe the statistical relation between multiple spatial locations and use the probability of occurrence of particular patterns.Using the generated 3D structure they predicted the absolute permeability of the sandstone sample using the effective medium theory and lattice Boltzmann simulation method.
Okabe and Blunt [18] used the approximated value of specific surface from 2D images to reconstruct 3D images of rock.In this regard, they tried to generate a 3D volume with the same specific surface as the estimated value from 2D images.A method of estimating the specific surface from 2D images has been previously developed by Adler et al. [19].
Politis et al. [20] presented a process-based and stochastic reconstruction method of porous media.Their method uses binary SEM images of rock to obtain a reconstruction method based on specific surface measurements.Blunt et al. [21] presented a vast field of experiences in the concept of pore scale imaging and modeling by providing a high resolution data bank of Micro-CT images of rock pieces.Mostaghimi et al. [22] utilized this data bank to present a suitable size of representative elementary volume (REV) for estimation of geometry-based properties such as porosity and specific surface area.
One main challenge in the application of multiple-point statistics to the simulation of three-dimensional (3D) blocks is the lack of suitable 3D training images.Comunian et al. [23] compared three methods to solve this challenge by attending (2D) training images: (a) the method based on aggregation of probabilities, (b) the method based on merging the lists obtained using the impala algorithm from diverse 2D training images and creating a list of compatible data events that is then used for the multiple-point statistics simulation, and (c) the method based on sequential simulations of 2D slices constrained by the conditioning data computed at the previous simulation steps.
In order to estimate the specific surface value of rock grains, cross-sectional images of different rocks have been utilized by coupling image analysis methods and statistical calculations.Initially, the specific surface of different rock structures has been calculated by computational methods with high precision using Micro-CT images provided by Blunt et al. [21].Sahimi [24] stated that photomicrographs of polished sections of a sample porous medium with sufficient contrast between the pores and the matrix can be used to determine the specific surface of rock grains.Sahimi [24] stated that from the relation between 2D (surface) measurements and the properties of the 3D system, an estimate of specific surface can be obtained.The main objective in this study is to correlate the three-dimensional value of grain specific surface versus 2D characteristics of rock sections.The study here also presents a straightforward relationship between specific surface and the ratio of particle perimeter by its area in 2D images.

Methodology
One of the main challenges in extracting data from 2D images of porous medium is how to section the medium.This can be normally done by various angles and positions while the perpendicular orientation is more common.The true values of porosity and specific surface may not be visualized in sectional images because of limited field of view in two dimensions.In this work, the key parameter in estimation of specific surface is the ratio of grain particle perimeter to the area In order to examine the specific surface of grains in more detail, each 3D sample of rock has been gridded to smaller cubic parts.Calculation of specific surface from 3D images in each lattice requires counting the pixels of the image.This however depends on the image resolution.Estimation of specific perimeter requires a specified pattern in sectioning.Lock et al. [25] utilized an averaging approach for the range of 180 degrees of cross section angle in order to find mean pore size of rock images.A similar approach of cross sectioning is used in this work to consider various probabilities and complexities of grain structure.Figure 1 illustrates the pattern of cross sections in the plane of and -axis named as   direction.The black and white parts lying under the cross lines present pores and grains of the rock, respectively.Each cross section indicates a 2D image which was analyzed by computer coding to find the specific perimeter.Next step is to average the data to upscale the results.Averaging specific perimeter in   direction can be written as where ⟨  ⟩ is normalized specific perimeter and   (  ) is a function that gives back the single value of specific perimeter for a cross section with an angle   .Here, ⟨  ⟩ is the average of specific perimeter for all possible angles of cross sections.Figure 2 presents a 3D view of one grid from a binary image of Berea sandstone as one of the rock types investigated in the current work.Binary data are provided from Blunt et al. 's database [21].Consequently, an averaged value for specific perimeter in each arbitrary section can be obtained as The plane that contains each of   ,   , and   is clarified in Figure 3.According to the integral definition, by simplifying the denominator, (6) changes to Although specific perimeter (  ) is a function of cross sectioning an angle, its dependency to sectioning angle cannot be formulized easily.Due to the unstructured and random shape of rock grains, the value of specific perimeter fluctuates randomly versus angle of cross section.Figure 4 shows the variations of specific perimeter for 6 different samples of Berea sandstone versus 90 degrees of cross section rotation.Apparently, no specific trend can be detected for describing the function's behavior.Also, there is a similar situation for one of the rock samples (Berea sandstone) in the different plane of angles labelled as   ,   , and   .Figure 5 presents the changes of specific perimeter for different angles of cross section from 0 to 180 degrees for this sample.As clarified in Figure 1, in order to have an overall view of rock specific perimeter, all possible angles of cross section should be observed and recorded for each rock piece.As illustrated in Figure 5, the value of specific perimeter at 0 ∘ is equal to its value at 180 ∘ , because of a complete rotation of the sectioning line around the sample.
Due to random behavior of specific perimeter for different angles (as visible in Figure 5), a numerical integration method (trapezoidal rule method) is used to solve (7).By estimation of the integrals using this method, the value of normalized specific perimeter has been obtained for each sample grid of rock.This calculation of integral is done using the developed computational code via MATLAB software for each angle in each grid of different rock types that are provided in Figure 6.
In the next step, it has been found that there may be a meaningful relationship between specific perimeter and specific surface of a porous material.Considering the rock grains as perfect spheres, the ratio between the specific surface and normalized specific perimeter ⟨  ⟩ can be stated as Here,  is the sphere radius and   is Geometric Coefficient.This ratio can be used for representing the ideal porous medium with considering the fact that real grain structures do not fully resemble spherical structures.

Results and Discussions
Binary 2D images of different studied rocks are provided in Figure 6.The black area represents the pore space in each image.The heterogeneity of carbonate samples is also visible in these images.The length at the bottom of each image shows the size of the square side.
In order to find the ratio between specific surface and normalized specific perimeter of porous materials, six different rock types are investigated.Berea sandstone as a well-known rock type shows almost uniform scattered data in estimation of Geometric Coefficient (  ). Figure 7 shows that the value of Geometric Coefficient varies versus normalized specific perimeter of 2D images for more than 500 samples.It can be concluded that most of the data points are almost centralized around an average value.Here, the value is about 1.37 for   .
The average   and standard deviation for its range of changes have been calculated for six different rock types using the algorithm presented earlier in this paper.Table 1 presents the detailed results of a running model on several samples.The average and standard deviation for each recorded parameter are presented.The average of Geometric Coefficient for different samples does not vary significantly.This shows that there is a reliable ratio between the specific surface and specific perimeter of rock grains.Although the carbonate samples are significantly heterogeneous and the standard deviation of their parameters has relatively large values,   does not change significantly.
As a key conclusion, the results of applying the present algorithm on 2184 rock samples from six different rock types leads one to declare an averaged value for Geometric Coefficient about 1.35.Size of each cubic sample is fifty pixels which show the real size by multiplying image resolution.The same procedure used in the Figure 7 sample has been applied to the other five rock types and the values of Geometric Coefficient have been plotted.Figure 8 illustrates the scattering of Geometric Coefficient for two types of sandstones and one sand pack, the specifications of which are available in Table 1.The ratio between specific surface and specific perimeter (Geometric Coefficient) for all three groups of samples is slightly changing around the value 1.4.
Figure 9 shows the severe fluctuation of Geometric Coefficient for two carbonate rocks.Apparently, it can be concluded that there is not any specific pattern in scattering data of those samples.Additionally, increasing the size of samples (representative elementary volume) does not indicate a symmetric pattern of scattering.
In order to evaluate the precision of the estimated specific surface, these values are plotted versus the measured specific surfaces for all six rock types.In this regard, hundred random cubes are chosen among each rock body and their specific surfaces have been estimated by multiplying the specific perimeter to 1.35.Figure 10 illustrates how data are scattered around the real values.The Coefficient of Determination which is also known as -Squared indicates the precision of prediction.Here the value is about 0.63 which is reasonable considering the heterogeneity of samples and variety of rock types.As it is visible, the carbonate samples show more diversity from the unit slope line.This indicates that the presented model cannot be used comprehensively to predict carbonates specific surface while it seems to be applicable for sandstones.

Conclusions
In this study, a Geometric Coefficient has been found to estimate the specific surface area () of rock grains from 2D images.Normalized specific perimeter ⟨  ⟩ is the averaged ratio between grains perimeter to their area in the images.
Based on analysis in this study in order to estimate specific surface area of rock grains one can rely on calculating and averaging the specific perimeter from several 2D images and then multiply the normalized value by 1.35.These 2D images may be obtained from thin section microscopy or any other photographical method which scans only the outer surface of the porous media.The results indicate that this approach is relatively robust for estimation of specific surface of sandstones.Using this approach is not recommended for carbonate rocks, and therefore future work is recommended.
Finally, the approach in this study can provide some preliminary observations of grain structures for modelling purposes when 3D data is not available.

Nomenclature
0 : Kozeny constant   : Geometric Coefficient   : Specific perimeter   : Cross section angle in the plane perpendicular to axis    : Cross section angle in the plane perpendicular to axis    : Cross section angle in the plane perpendicular to axis  ⟨  ⟩: Normalized (averaged) specific perimeter ℎ: Hits, number of pore spaces recumbent under the needle : Needle length : Specific surface : Cuts, number of particle perimeters intersect by needle : Permeability : Sphere radius : Porosity.

Figure 1 :
Figure 1: Schematic pattern of 180 degrees of cross sectioning in direction of   .

Figure 3 : 6 Figure 4 :
Figure 3: Clarification of angles and their plates for integral on cross sections.

Figure 5 :
Figure 5: The fluctuation of specific perimeter for different planes of angles   ,   , and   in Berea sandstone.

Figure 6 :
Figure 6: Images of used rock types and textures.

Figure 7 :Figure 8 :
Figure 7: Plot of Geometric Coefficient versus normalized specific perimeter for Berea sandstone.

Table 1 :
Specific surface and perimeter and Geometric Coefficient (  ) calculated by model.