Uncertainty of volume fraction in bimrock using the scan-line method and its application in the estimation of deformability parameters

The volume fraction within a bimrock or bimsoil is an essential parameter that is useful for estimating the engineering properties of heterogeneous geomaterials. This paper presents analytical and numerical solutions to quantify the uncertainty of volume fraction measurements in bimrock/bimsoil using a scan-line method. The analytical solutions for the mean and variance of volume fraction estimates are based on a representative volume element model. The numerical solution is obtained through simulations of scan-line measurements. This work also employs physical tests using CT scan images from artificial bimrock/bimsoil to validate these solutions. The results demonstrate that the uncertainties of the volume fraction depend on the magnitude of the volume fraction of the blocks, the diameter of the blocks, and the length of the scan line. The proposed analytical and numerical solutions are compared with existing physical experimental tests and analytical solutions. An illustrative example to demonstrate the estimation of the uncertainty of volume fraction using the scan-line measurement is present. Finally, an example application of the volume fraction characterization in the geological engineering, in terms of Young’s modulus estimation and characterization, is provided.

The V f of interest could be the portion of discrete-phase inclusions (e.g., certain minerals, blocks, and particulates) or continuous-phase matrix in a composite material. In this study, the volume fraction of a block, V b , is defined as the ratio of the block volume over the total volume, which is the same as the volumetric block proportion, VBP, established in some literature (Lindquist 1994;Lindquist and Goodman 1994;Kahraman and Alber 2006, 2009Kahraman et al. 2015;Kalender et al. 2014;Medley 1994Medley , 1997Medley , 2001Medley and Goodman 1994;Medley and Sanz Rehermann 2004;Sonmez et al. 2004Sonmez et al. , 2006aSonmez et al. , b, 2016Tsesarsky et al. 2014;Tien et al. 2010).
Measurements of V f can be made through stereology and can be obtained via 0D, 1D, 2D, or 3D measurements depending on the dimension of the sampling windows used (Russ and Dehoff 1999). That is, 0D measurements are sampled by points, 1D by lines, 2D by areas, and 3D by volumes. Thus, "volume" is used in a generalized sense, with the resulting "volume" fraction being the ratio of the phase "volume" over the total volume sampled. That is, the volume fraction is defined as the point fraction, P f , for 0D, the line fraction, L f , for 1D, the area fraction, A f , for 2D, and the volume fraction, V f , for 3D. The Delesse principle (1847) shows that all these different dimensional measurements can lead to the same results, namely: It is difficult to determine the true volume fraction for bimrocks/bimsoils, especially for large scales. The common practice is to employ either 1D or 2D measurements, or both (Sonmez et al. 2004). Many investigators resolved to use aerial fractions of outcrops (2D surfaces), augmented with field and laboratory tests to obtain the geometric information of bimrock/bimsoil from outcrops (Medley 1994;Xu et al. 2008Xu et al. , 2011Coli et al. 2012;Kahraman et al. 2015;Ymeti et al. 2017;Liu et al. 2018;Meng et al. 2018;Yang et al. 2019). However, it is not always possible to obtain areal fraction of blocks through image analysis if the blocks are dyed by matrix or the chromatism between blocks and matrix is small. As illustrated in Fig.  1, 1D measurement is perhaps the easiest and most efficient way to obtain an aerial fraction of block (A b ) in both laboratory and field scenarios. The 1D measurement has also been used for the estimation of area fraction A f of a surface (i.e., a cross-section, an image), as in the estimation of lake area from remote sensing images (Stein and Yifru 2010) and air voids from concrete surface (ASTM 2012).
Further, according to the Delesse principle (1847), the expected value of A b is equal to V b if the sampling area is sufficiently large such that it can represent the geometry of true bimrocks/bimsoils. In this paper, we assume that scan-line measurements of the outcrops, images, or surfaces can reveal the 3D geometry of bimrocks/bimsoils. Accordingly, this paper uses V b to express the areal fraction of block.
To avoid sampling bias, 1D lines sampling should be applied systematically. This is readily achieved using equally spaced, straight parallel lines (Fig. 1) that capture spatial variability within a given system. This systematic 1D line measurement is referred to herein as the scan-line method. To apply the scan-line method, the length of the sampling lines has to be selected, which may be estimated based on the confidence intervals. Quantification of the uncertainty in V b is thus the first step to apply the scanline method.
Various approaches have been proposed in the past to guide the selection of the required sampling length. For instance, some rules of thumb have been proposed. Holmes (1921) recommended that the total length of the lines be 100 times the average diameter of the grains being measured. Krumbein and Pettijohn (1938) considered that a total traverse length equal to 1000 times the largest particle would give a fair degree of accuracy. ASTM C457-11 (2012) requires the line length to be 1000 to 3000 times the maximum diameter of the void. Analytical solutions provide improved recommendations, but so far are limited to the cases where V f is small. Hilliard and Cahn (1961) presented an analytical solution, which is only valid when the V f of a void is less than 5%, as follows: where E(•) represents the expected value, CV(•) represents the coefficient of variance, L indicates the total length of the scan line, and L c is the length of the interception when the scan line passes through a void. If the particle shape is circular, Eq.
(2) can be simplified as follows: Medley and Goodman (1994) developed empirical solutions. They applied hand-tracing on a photograph of mélange cores to evaluate the convergence of V b with respect to the scan-line sampling lengths. Medley (1994) further prepared a series of bimrock samples when investigating the uncertainty associated with volume fraction estimation, with the known V b varying within the range of 13 to 55%. Medley further developed a chart to determine the uncertainty factor for various scan-line lengths and V b (Medley 1997(Medley , 2001.
In the following, the concept of "represented volume element," or RVE, is introduced first. A new analytical solution on V b measurement uncertainty that is based on RVE is then presented. This is followed a detailed description of the verification efforts. The verification is carried out first by numerical simulation and further by applying to controlled physical tests. A new empirical equation is derived from the numerical simulation results. Finally, an example is presented for estimation of the uncertainty of volume fraction; the results are further used to estimate the deformability parameters, including the shear modulus and Young's modulus.

Representative volume element
The measurement of the volume fraction of inclusions within a rock matrix is used as an example for introducing the RVE concept. Consider the case that circular inclusions of the same size are randomly and uniformly distributed within a 2D matrix, RVE can be thought of a square element with a circular inclusion that has a volume fraction within the element the same as V b to be measured. A randomly drawn scan line, as shown in Fig. 2, would repeatedly hit and miss the inclusions. If the RVE is superimposed and aligned with the nearby circular inclusion that intersects with the scan line, as depicted in Fig. 2, it becomes clear that estimates of V b using scan lines can conceptually be obtained by simply studying a single RVE. When the blocks are uniformly and randomly distributed within a matrix, a straight line of length L becomes N line segments randomly drawn inside the RVE, where NL s equals L and L s is the width of the RVE. From this, a theoretical derivation makes it possible to estimate V b using random variables.
It is important to note that RVE, the natural of the phase, inclusion determines how high V b can go. For instance, in a 2D setting with uniform circular inclusions the maximum V b that can be attained is 78.5% (= πr 2 / (2R) 2 ) when the circular inclusions circumscribe the boundary of RVE. While in a 3D setting, since a sphere circumscribes to the boundary of a cube RVE can work with V b only up 52.4% (= 4/3πr 3 /(2R) 3 ). But in reality, these limits can be broken when the inclusions have different sizes. To model that, the inclusions should be composed of polydisperse circle or spheres within an RVE. Such formulations are presented also.

Mathematical derivation
The analytical solution is first derived with the following considerations (assumptions): (i) The inclusions are circular; (ii) the scan lines are assumed to be vertical within a RVE; (iii) the inclusions distribution within the matrix follows isotropic, uniform, and random (called IUR) conditions; (iv) the RVE is square with a width L s , contains a circular inclusion of diameter D; (v) the value of V b is in the range of 0 to π/4." It thus gives the expected V b From the RVE model (Fig. 3), a scan line passing through x would intercept the inclusion with an intercept length, L c , as follows, The relationship of L c with V b is simply which further gives The variance of linear fraction of block (L b ) in an RVE, V REV (L b ) can be expressed by And Eq. (8) can also be expressed by L c and L s as shown as follows: And Eq. (9) can be rewritten as Substituting Eq. (7) into Eq. (10), Eq. (10) can be represented using L c and L b , as follows: As shown in Fig. 2, each scan line can be viewed as a summation process. Thus, the statistics of its L b will converge to a Gaussian distribution, which can be explained through the central limit theorem (CLT) (Wackerly et al. 2008). The mean value and uncertainty of L b estimation at various sampling sizes (N = L/L s ) can be determined as follows: where V(•) represents variance. By substituting Eq. (6) and Eq. (11) into Eq. (12), the following is obtained: Equation (13) can determine the uncertainty of the linear fraction of a block by accessing the length of the scan line, the intercept lengths, and the expected linear fraction of the block, which can be obtained from a scan-line measurement.
Although a block diameter cannot be measured directly using the scan-line method, D can be calculated using Eq. (7). Hence, Eq. (13) can be rewritten as Equation (14) is the same as the analytical solution proposed by Tien et al. (2010), which can estimate the uncertainty of a scan-line measurement under monodisperse block inclusion.  (14) cannot be used under this condition. To tackle this issue, this study employs the concept of the "total sum of squares" (TSS) (Wackerly et al. 2008), which includes the "sum of squares for treatment" (SST, the first term of Eq. (15)), and the "sum of squares for error" (SSE, the second term of Eq. (15)), where SST can be viewed as the measurement uncertainty from the block size variability and SSE can be viewed as the measurement uncertainties from the spatial probability of each block size: A polydisperse inclusion ( Fig. 4a) can be reassembled as shown in Fig. 4b, with each block size having the same A f , indicating there are no A f variations between each block size. Thus, the SST is equal to zero, and the TSS is equal to the SSE. Hence, Eq. (15) can be rewritten as: Then, the total number of measurements in Eq. (16) is divided, which can be rewritten as follows:  where f(d) represents the possibility density function of intercepts from the ith block size, which can also be represented as the block proportion, Ci = ith block area / total block area.
Substituting Eq. (14) into Eq. (17), and allowing the TSS to have an equivalent diameter, D e , the following is obtained: Equation (18) can also be expressed in series as follows: Substituting Eq. (7) into Eq. (18), the resulting equation can be rewritten as: Assuming that for each block E(L c ) = a i E(L c ) eq , where a i is a multiple factor based on the equivalent block size and E(L c ) eq is an equivalent intercept mean, Eq. (15) can then be rewritten as The results of Eq. (21) are the same as those of Eq. (7), which means Eq. (13) can be used to determine the uncertainty of L b in polydisperse inclusions. In addition, Eq. (14) should be rewritten as follows: For a reliability design or other in situ investigations, most engineers apply CV ¼ ffiffiffiffiffiffiffiffi ffi V : ð Þ p =E ⋅ ð Þ to determine the uncertainty of the measurement. Equations (13) and (22) can therefore also be presented as a CV, which are expressed as Equations (23) and (24) present an analytical solution for the uncertainty of the volume fraction measurements using scan-line method. This analytical solution is applicable when

Numerical simulations
This paper presents a MATLAB code that can address a volume fraction measurement in a bimrock/bimsoil using a scan line. A 2D domain with uniform and random circular block dispersions in a square domain was generated. A series of scan lines were placed on the feature. A series of L b were obtained. The mean and standard deviation of L b were calculated.

Periodic boundary
While generating the 2D domain, it is impossible to avoid wall (or edge) effects (Bentz and Garboczi 1999), which result in non-uniform dispersion near the edges. To reduce the effect of the edge, four strips (2D in width) were trimmed from the four edges of an extended square (the generated domain) (Tien et al. 2010). However, this approach may result in V b scatter. To tackle this issue, a periodic boundary is applied in this study. If a block intersects one edge of the boundaries, the outside part of the block will be mapped onto its opposite side, as shown in Fig. 5. The advantages of the periodic boundary method are to eliminate the wall effect and to generate features with the precisely desired volume fraction.  71, 9.43, 18.9, 28.3, 37.8, 47.1, 56.5 30, 60, 120, 240 10:1:1, 3:2:1, 1:1:1, 1:1:10 10 cm, 5 cm, 2.5 cm a Block proportion C 1 :C 2 :C 3 of block diameters equal to 10 cm, 5 cm, and 2.5 cm, respectively.

Sizes of domain and block
The size of the domain should be sufficiently large to accommodate a large number of blocks embedded into the matrix to meet the conditions of being statistically isotropic and homogeneous. In this study, the size of the domain is 400 cm × 400 cm. The diameters of the blocks are 2.5 cm, 5 cm, and 10 cm.

Simulation procedures
The numerical simulation procedures are as follows: I. Define a square domain of size B × B II. Determine V b within this domain III. Place circular blocks randomly, and do not allow any overlap between blocks IV. Generate a large number of scan lines (say 10,000) in Fig. 6. Then, calculate linear fractions of blocks, L b , for each line Repeat steps (1) through (4) 100 times for each V b ; thus, this work regenerates a heterogeneous rock mass 100 times for each V b and calculates the total mean L b (100 × 10,000) and CV.

Parametric studies
Parametric studies were performed to evaluate the effects on the uncertainty of volume fraction estimation. Several  Medley's (1997) physical model tests parameters, block sizes, relative lengths of scan lines, L/D e , and volume fractions of blocks, V b , were selected. Three block sizes were used: D 1 = 10 cm, D 2 = 5 cm, and D 3 = 2.5 cm, including monodisperse and polydisperse types. Four different relative lengths of scan lines (L/D e ) were studied, including L/D e = 30, L/D e = 60, L/D e = 120, and L/D e = 240. The range of V b is from 3.14 to 70.1%. The parameters for the simulated bimrock/bimsoil are as shown in Table 1.

Simulation results
Parametric study for the relative length of a scan line The simulation results are shown in Fig. 7. Cases of lower and higher V b are given as examples. It can be clearly seen that V b affects the uncertainty and distribution of L b even under the same dimensionless length of a scan line (or sampling size). For lower V b (= 0.0506) (shown in Fig. 7a as a lower dimensionless length of a scan line, L/ D e = 30), the distribution of L b shows a non-normal or a non-log-normal distribution. As L/D e increases, the distribution becomes more similar to a normal distribution, as shown in Fig. 7a-d. For higher V b (= 0.506) for L/D e = 30 (shown in Fig. 7a), the distribution of L b has an approximately normal distribution. The CV of L b is inversely proportional to the square root of L/D e , which can be predicted by the analytical solution (Eq. (19)) and explained based on the CLT.
Equivalent diameter of the blocks For polydisperse cases, the equivalent diameter of the blocks can be calculated using Eq. (19). For example, the diameters of blocks D 1 , D 2 , and D 3 are 10 cm, 5 cm, and 2.5 cm, respectively, and the ratio of block proportions is C 1 :C 2 :C 3 = 6:4:2. The equivalent diameter of the block is D e ¼ 10 Â 6=12 þ 5 Â 4=12 þ 2:5 Â 2=12 ¼ 7:08 m: There is a correlation between the measurement uncertainty and block sizes, as shown in Eq. (19). From previous parametric studies, the uncertainty is proportional to the square root of the sampling size. Therefore, CV(L b ) multiplied by (L/D e ) could be a constant at a particular V b . Other similar  (Tien et al. 2010(Tien et al. , 2011(Tien et al. , 2012(Tien et al. , 2015 also concluded this phenomenon. A parametric study of block sizes is shown in Fig. 8. All the uncertainties of V b estimations for the various block sizes are located on a curved line. This means the numerical simulations for various block sizes followed Eq. (19).
Regression curve for numerical solutions Based on the results of the numerical simulations, a relationship between CV, V b , and L/D e was obtained via a regression procedure (R 2 = 0.99): A comparison of the analytical and numerical solutions is shown in Fig. 8.

Validation: physical model test and CT scan images
The results of the physical model developed by Medley (1997) and the CT scan images are adapted to validate the analytical and numerical solutions.

Physical model
The physical bimrock models were composed of plaster of Paris matrices in which were embedded ellipsoidal blocks to  Experimental data Hilliard and Cahn (1961) This study

Analytical solution
Numerical solution Medley's (1997) physical investigate the uncertainty of a V b estimate when using a scan line. The V b of physical bimrock models ranged from 13 to 55%. The particle size distribution followed a specific fractal dimension (Medley 1997), which was in agreement with the particle size distribution (PSD) of mélanges from San Francisco. However, the PSD was for 3D bimrock features. Hence, obtaining the D e from a 3D PSD is necessary. If a surface cuts though a 3D monodisperse block inclusion, then D e can be calculated from Eq. (19) as follows: in which D 3D is the block diameter in a 3D feature.
If a surface passes though the 3D polydisperse inclusion, then, by substituting Eq. (26) into Eq. (18), and using a PSD for the expression, the following is obtained: where f(d PSD ) represents a PSD function. Medley's (1997) physical model included five different elliptical block sizes, whose minor, middle, and major semi  . 13 Comparison of the CVs of block volume fractions: experimental data, analytical solutions, and numerical solution axes were respectively 3, 3, and 6 mm; 6, 6, and 12 mm; 12, 12, and 24 mm; 24, 24, and 48 mm; and 48, 48, and 96 mm, respectively. If the orientation of an elliptical block follows uniformly random conditions, then the 1D measurement results (mean and uncertainty) will be close to the spherical block dispersion (Tien et al. 2011). The elliptical blocks can be represented as an effective sphere, which has the same volume as the ellipsoids. A PSD of the blocks in Medley's physical model is shown in Fig. 9. After substituting the PSD (Fig. 9) into Eq. (28), the equivalent block diameter, D e , is equal to 32.5 mm. The detailed calculations are shown in Table 2.
Furthermore, this work compares the proposed analytical solution with Eq. (24), the numerical simulation, and Medley's model. The results indicate that the proposed solution is in good agreement with Medley's physical model tests, as shown in Fig. 10 and Table 3.

CT scan images
This study uses a CT scan to obtain internal cross-section images from Iceland sand and crushed rock.

Iceland sand
This sample consists of Iceland sand and an epoxy matrix (Gräser, personal communication). The resolution of case I is 669 pixels × 669 pixels, corresponding to a real size of 6.21 mm × 6.21 mm. The original CT scan images are in grayscale. To separate the block and matrix more clearly, the grayscale images are enhanced to binary images through image processing (shown in Fig. 11). A scanline analysis code (SAC) was developed to detect which blocks would be intersected by scan-lines under binary images. The number of scan lines in the SAC depends on the image resolution, i.e., if there are 669 pixels in the horizontal direction, then there will be 669 scan lines in the vertical direction. The SAC records the "linear fraction of block" from each scan line and the "intercept length" from each intercepted block.
In case I, the scan-line length, L, is 6.21 mm, and 669 scan lines are generated in the SAC. The SAC also recorded 12,511 intercept lengths, which can be substituted into Eq. (21), and D e = 0.366 mm. A mean of L b = 55.4% and a standard deviation of L b = 7.83% were obtained.
Comparisons of Hilliard and Cahn's (1961) analytical solution (Eqs. (2) or (3)), the proposed analytical solution, and the numerical solution are shown in Table 3. They indicate that the proposed analytical and numerical solutions agree well with statistical results from case I.

Crushed rock
This sample consists of crushed rock (from Kaoping River, Taiwan) and an epoxy matrix, identified as case II (Fig. 12a). The image size is 483 pixels × 483 pixels, and its physical size is 1.95 cm × 1.95 cm. Both cases I and II are high V b (more than 50%) images. To obtain images (cases III-VII) with low V b , a computer program was used to delete blocks in case II randomly, as shown in Fig.  12b-f.
Comparisons of the CVs of the block volume fraction for the CT scan images, Hilliard and Cahn's (1961) analytical solution, and the present study are made in Table 3 and Fig. 13. It can be clearly seen that the proposed analytical solution and the numerical solution compare very well with the measured data. However, Hilliard and Cahn's (1961) solution will overestimate in higher V b conditions.

Volume fraction estimation
The location shown in the illustration of the scan-line method is at + 17 km, No. 114 Provincial Highway, Taoyuan County, Taiwan. The geology of this location is  Table 5 The elastic parameters of laterite-gravel formation (Lin 1986) Fig. 1b (called case VIII). The boundary of the block and matrix is 32 mm (Church et al. 1987). The intercept lengths of the blocks are shown in Table 4. From this table, L b can be obtained as follows: After substituting the L c data into Eq. (21), D e can be determined using: Then, by substituting L b , D e , and L into Eq. (24), the following is obtained: According to a numerical simulation of this case, the histogram of the linear fractions of the rock is approximately a normal distribution, as shown in Fig. 14. The standard deviation is equal to 0.0254. The 95% confidence interval of L b is 0.338 ± 2 ⋅ 0.0254. Based on the data of the scan-line field measurement, the equivalent Fig. 16 The shear modulus of laterite-gravel formation vs volume fraction of block Fig. 15 The bulk modulus of laterite-gravel formation vs volume fraction of block diameter of the blocks is 6.29 cm. According to the Delesse principle (1847), the expected V b is 0.338 with a 95% confidence level in the interval of 0.287 to 0.389. In most geotechnical parameter investigations, CV < 10% corresponds to "low variability" (Phoon et al. 1995). However, if the acceptable CV is more tolerant, e.g., CV = 20%, then the requirement is L = 1.55 m.

Estimation of deformability parameters (shear modulus and Young's modulus)
It is generally difficult to determine the deformability of laterite-gravel by conventional laboratory and field tests. The deformability of laterite-gravel may be calculated based on micromechanics, for example, self-consistent scheme (Hershey 1954;Hill 1965) and differential scheme (McLaughlin 1977). The required input parameters for micromechanics include bulk moduli (K) and shear moduli (G) of laterite and gravel (Lin 1986) (as shown in Table 5) and V f of gravel (denoted herein as V b ). The elastic ratio of gravel and laterite, G of gravel over G of laterite (G G /G L ) or K of gravel over K of laterite (K G /K L ), are about 1100. The elastic modulus ratio (Young's modulus over uniaxial compressive strength (E G /UCS G ) is about 500. Based on the data from the scan-line measurement of case VIII (in section 4.1), the mean and CV of V b are 0.338 and 0.0751, respectively.
In this paper, we utilize differential scheme (Eq. (32) to Eq. (33)) assuming V b as a normal random variable and using Monte Carlo approach. where Substitute K L (which is equal to K when V b = 0), G L (which is equal to G when V b = 0),K G , G G , and the normal random variable V b into the differential scheme (McLaughlin 1977). Then, using the 4th order Runge-Kutta method (Butcher 1996) to solve the partial differential equations of the differential scheme, the bulk modulus,K, and shear modulus,G, of laterite-gravel formation under a specific V b value can then be obtained. The simulation results of bulk modulus and shear modulus are shown in Figs. 15 and 16, respectively. The result shows that the influence of stiffness of gravel is insignificant on the overall deformability of laterite-gravel mixture. In addition, Young's modulus (E) and Poisson ratio (υ) can be also obtained by substituting bulk modulus and shear modulus into Eq. (34) and Eq. (35).
Shear modulus and Young's modulus are important design parameters in geological engineering. The histograms of shear modulus and Young's modulus obtained in this study are shown in Figs. 17 and 18, respectively. These histograms indicate that both design parameters follow approximately normal distribution. The mean and CV of shear modulus and Young's modulus are shown in Table 6, which can provide an estimate of the 95% confidence interval of shear modulus in the range of 22.4 to It should be noted that the deformability of bimrocks/ bimsoils often exhibits spatial and tempo variations (Chen et al. 2019). These variations may be induced by the geometry of bimrocks/bimsoils such as the volume fraction of block. With the proposed approach, the variation of the deformability parameters can be estimated, and the results may be used, for example, for reliability-based design of engineering systems in the bimrocks/bimsoils (Chen et al. 2019).

Summary and conclusions
This paper presented analytical and numerical solutions to quantify the uncertainty of volume fraction measurements in bimrock/bimsoil using a scan-line method. The proposed method had been validated using published physical experimental data and analytical solutions. As the volume fraction estimated with the scan-line measurements is linked to deformability parameters, such as modulus of elasticity and shear modulus, an estimate of these parameters along with their uncertainty characteristics can readily be made. Such information enables the deterministic, as well as the reliability-based, analysis and design of engineering systems in the bimrocks/bimsoils. I. This paper proposed an analytical solution to quantify the uncertainty of volume fraction estimation in a bimrock/bimsoil using a scan-line method. A geometry model based on representative volume elements was adapted to derive the uncertainty of the volume fraction of a block. The solution is a function of the volume fraction of the block (V b ), the diameter of the block, and the length of the scan line. The applicable range of V b is between 0 and π/4. II. For the numerical approach, a 2D domain was generated, which used a "periodic boundary" to avoid the "wall effect." The results of the numerical simulations agreed well with the analytical solution.
III. Experimental results of the physical model and CT scan images were adapted to validate the analytical and numerical solutions. The proposed analytical and numerical solutions showed a good agreement with the experiments. IV. The proposed analytical solution for estimation of the volume fraction is further extended into a deterministic micromechanical model to estimate the mean and variation of bimrock's mechanical properties such as the shear modulus and Young's modulus. The variation of these deformability parameters is also characterized.