Probabilistic characterization of Brazilian tensile strength of rock from both direct and indirect measurements

The determination of Brazilian tensile strength (BTS) is an essential step in the analysis and design of mining structures at a particular site. The estimated design parameters are affected by the inherent variability of BTS, measurement errors arising from laboratory testing, and transformation uncertainty associated with the empirical model linking BTS to other rock properties when it is indirectly estimated. These uncertainties are usually lumped together as the total variability of BTS. However, it is the inherent variability resulting from natural geological factors, not the total variability, that directly affects the actual response of rock structures. Hence, there is a need for proper characterization of the inherent variability of BTS while the measurement errors and transformation uncertainty are explicitly incorporated. This paper develops a Bayesian approach which uses sequential updating, that is multi-input oriented for probabilistic characterization of the inherent variability of BTS of rock. The proposed approach systematically combines previous engineering experience and site information from both direct BTS data and data from indirect tests like point load test to inversely infer the inherent variability of BTS. The proposed approach quantitatively accounts for the effects of measurement errors and transformation uncertainty on the characterization of the inherent variability of BTS. The proposed approach is illustrated and validated using real-life data and simulated data. The result shows that the proposed approach provides a proper characterization of the inherent variability of BTS based on available information from multiple sources. Sensitivity studies are also performed to explore the effects of measurement errors on the performance of the proposed approach.


Introduction
The determination of rock properties is indispensable in mining engineering, providing necessary input for analyses and designs of mining and geotechnical structures.Rock parameters vary from location to location within a rock deposit/site and among different rock sites.This is because rocks are formed from various geological processes and these processes affect the mineral constituents, packing, and properties of rocks.This variability (often referred to as inherent variability), which is natural and arises as a result of geological processes leading to the formation of rocks is one of the uncertainties in rock parameters (e.g., Baecher and Christian 2003;Aladejare 2016;Aladejare and Wang 2017a).During the determination of characteristic values of rock parameters, there are also measurement errors which arise as a result of imperfections.The imperfections may include but are not limited to malfunctioning of test equipment and operator and procedural errors, especially when there is a deviation from the recommendations of the International Society of Rock Mechanics (ISRM) (Ulusay and Hudson 2007) or rock testing and parameter determination.There is also transformation uncertainty, which is introduced when a rock parameter (e.g., Brazilian tensile strength BTS) for mining design and analyses are estimated from indirect measurements (e.g., point load index Is (50) values) through the use of a transformation model (e.g., regression/empirical model linking Is (50) values to BTS).Using transformation models is very prominent in mining engineering because of 678 Page 2 of 10 the difficulties in performing laboratory testing for determination of some rock parameters, especially those requiring special sample preparation (e.g., Shalabi et al. 2007;Diamantis et al. 2011;Kurtuluş et al. 2010;Sharma et al. 2011;Heidari et al. 2012;Kahraman et al. 2012;Khandelwal 2013;Kallu and Roghanchi 2015;Fakir et al. 2017).
The inherent variability which is the actual variability of a rock site is usually lumped together with the measurement errors and transformation uncertainty as the total variability of estimated rock parameters from site investigation (Phoon and Kulhawy 1999).However, using total variability instead of actual variability leads to overestimation of uncertainties in rock parameters.This might lead to design errors in which the expected performance or response of rock structures varies significantly from the actual performance or response at rock sites (Christian and Baecher, 2011;Aladejare and Wang 2017b).Hence, there is a need for rational quantification of various uncertainties in estimated rock/geotechnical parameters in order to properly characterizes the inherent variability of rock properties at a specific site (Wang et al. 2016).This will help to ensure the estimation of the actual performance or response of the rock structure at a site.Performing probabilistic characterization of BTS inherent viability remains a challenge because attention has not been given to BTS unlike uniaxial compressive strength, Young's modulus and the deformation modulus of rocks that have been researched (e.g., Wang and Aladejare 2015;2016a;Aladejare and Wang 2019).Furthermore, there is difficulty in separating inherent variability from measurement errors for direct measurements (e.g., BTS values from Bralizian tensile tests) and transformation uncertainty arising when indirect measurements (e.g., Is (50) values) are used to estimate BTS values.This problem is further compounded by the lack of a sufficient number of site-specific test results that can ensure meaningful statistical analysis for mining designs.It is timeconsuming and costly to prepare rock samples that meet the ISRM standards for most rock strength tests.In addition, for instance, it may be difficult if not impossible to extract standard cores from weak or highly fractured, thinly bedded, and/or block-in matrix rocks (Gokceoglu 2002;Aksoy et al. 2010).
Bayesian methodology provides a logical route to formulate the quantification of inherent variability of rock parameters as an inverse analysis problem.Bayesian approaches provide the platform to quantify the inherent variability of rock parameters from data observed through direct and/ or indirect measurements.Bayesian approaches have been used to determine rock and geotechnical parameters as well as solve other mining engineering problems by integrating information obtained from different sources during site investigation (e.g., Wang et al. 2010;Müller et al. 2014;Wang andAladejare 2015, 2016b;Wang et al. 2015;Cao et al. 2016;Ng et al. 2017;Aladejare and Wang 2019;Aladejare et al. 2021;2022a, b).These previous studies used multisource information that includes prior knowledge available about a rock parameter before a rock project or laboratory analysis, site test data from indirect measurements, and transformation models relating indirect measurements to rock parameter of interest.Aladejare and Wang (2019) likened the process to data fusion, in which information from two or more data sources are combined to provide a more accurate estimation than any of the data sources individually.However, it is not common or reported in literature that indirect testing procedures (e.g., point load tests) are combined with direct testing (e.g., Brazilian tensile tests) to improve the characterization of BTS and help to reduce the effect of limited data often obtained for BTS.Ignoring direct measurements of BTS may not lead to a full depiction of the information on BTS at a site because direct measurement of rock parameters like BTS can provide vital information during the updating of the statistics and probability distributions of BTS from multiple sources of information.The statistics and probability distributions of rock parameters are helpful in understanding the updating processes and for checking the quality of information from different sources (i.e., both direct and indirect test data).It is therefore important to rationally incorporate direct and indirect measurements into the probabilistic characterization of inherent variability of rock parameters.There is no research that has performed probabilistic characterization of BTS inherent variability at a rock site through a systematic combination of results of indirect measurements like Is (50) values for which transformation uncertainty is unavoidable when it is used to estimate BTS and direct measurements of BTS which are associated with measurement errors.Logical depiction and understanding of the effects of measurement errors and transformation uncertainty are critical to the probabilistic characterization of BTS.
The study develops a Bayesian approach that is based on a sequential updating process to perform probabilistic characterization of the inherent variability of BTS of rock at a rock site by integrating prior knowledge with site information from both direct (BTS from Brazilian tensile tests) and indirect (e.g., Is (50) values from point load tests) measurements.Aladejare and Wang (2019) developed a Bayesian sequential updating approach that systematically combines values of rock mass rating (RMR) and tunneling quality index (Q) for the probabilistic characterization of rock mass deformation modulus.Aladejare et al. (2021) developed another Bayesian sequential updating approach that systematically combines values of three different punch tests (i.e., Schmidt rebound hardness (SRH), block punch index (BPI), and point load strength (Is (50) )) for the probabilistic characterization of uniaxial compressive strength (UCS) of rock.However, both cases (i.e., Aladejare and Wang (2019) and Aladejare et al. (2021)) used indirect test data without taking the direct measurements into consideration.The proposed Bayesian approach sequentially makes use of a limited number of Is (50) values from point load tests and BTS values from Brazilian tensile tests at different updating levels.The approach quantitatively accounts for measurement errors and transformation uncertainty during the probabilistic characterization of the inherent variability of BTS at a rock site.The paper starts with the probabilistic modeling of the inherent variability of BTS at a rock site, followed by the development of the Bayesian approach.The proposed approach is then illustrated and validated using real-life and simulated data.In addition, sensitivity studies are performed to explore the effects of measurement errors on the probabilistic characterization of the inherent variability of BTS at a specific rock site.

Probabilistic modeling of inherent variability of Brazilian tensile strength
For the probabilistic modeling of the inherent variability of Brazilian tensile strength in a rock formation or deposit, this study uses a normal random variable to represent the BTS.Many studies in the literature have modeled rock properties as normal random variables (e.g., Sari and Karpuz 2006;Wang and Aladejare 2015;2016a).Therefore, let B and B denote the mean and standard deviation of BTS, respectively.Then, BTS can be expressed as (Cao et al. 2016;Wang and Aladejare 2016b): where Z = standard normal random variable.Using Eq. ( 1), the inherent variability of BTS in a rock formation or deposit is quantified by B and B .In the next section, the Bayesian approach which is based on multi-level updating is presented.The approach updates the knowledge of B and B for the probabilistic characterization of the inherent variability of BTS at a particular rock site based on prior knowledge and information from both Is (50) data and direct BTS measurements from Brazilian tensile strength tests.

Bayesian framework
Test results from different laboratory experiments might not be obtained simultaneously during laboratory investigation into the properties of rock.For example, in this study, Is (50) values are measured in the laboratory through point load tests while direct BTS measurements are obtained from Brazilian tensile tests in the laboratory.It is of practical interest to note (1) BTS = B + B Z the statistics (e.g., B and B ) of BTS according to indirect testing data (e.g., Is (50) values) from the point load tests and how these statistics change as BTS values from Brazilian tensile tests are taken into account.Under the proposed Bayesian framework, the updating of B and B is partitioned into two levels and is sequentially performed to see the evolution of the statistics of B and B at each updating level.Note that updating of B and B performed at one level using both Is (50) from indirect point load tests and BTS from direct Brazilian tensile tests (e.g., Chen and Gilbert 2017;Cao et al. 2016) provides the same probabilistic characterization results as those obtained from the Bayesian approach that is based on sequential updating, but it cannot reflect the evolution of the statistics of BTS as direct BTS measurements are considered in addition to Is (50) values.In this study, Is (50) values are used at the first updating level to update the knowledge of B and B for the probabilistic characterization of the inherent variability of BTS.Direct BTS measurements are then added at the second updating level to further update knowledge of B and B of BTS.
For the two updating levels of the proposed Bayesian approach, the posterior distributions of B and B are used to reflect their updated knowledge (e.g., Cao et al. 2016;Aladejare et al. 2021;2022a)  Equation (4) quantifies the combined or updated information on B and B from the prior knowledge, point load Is (50) , data (i.e., DATA I ), and BTS values (i.e., DATA II ).

Likelihood function for indirect measurements
Point load test measurements, Is * (50) , can be related to BTS * values from Brazilian tensile strength tests by regression model from Sulukcu and Ulusay (2001) as given in Eq. ( 5).
Equation ( 5) above can be transformed to incorporate the transformation uncertainty which is introduced to the estimation of parameters when regression models are used for parameter estimations.The transformation model of Eq. ( 5) is given as where a = 0.4348 ; and T is a normal random variable that has a mean of zero and a standard deviation of T = 2.2700 and it represents the transformation uncertainty associated with Eq. ( 6).
In Eq. ( 6), BTS * represents the measured value, rather than the actual value of BTS that is modeled by Eq. ( 1) of Brazilian tensile strength obtained from Brazilian tensile strength tests because regression models are usually developed by fitting the empirical relationship between rock parameters based on their respective measurements.Therefore, BTS * is associated with measurement errors which arise during Brazilian tensile tests and can be expressed as a summation of the measurement error M and the actual value (i.e., BTS) of Brazilian tensile strength, as follows: where M is a normal random variable.The mean and standard deviation of M denoted by M and M , respectively, are taken as 1 and 0.1 MPa.These adopted statistics for M are consistent with the typical values reported in the literature for most geotechnical parameters (e.g., Phoon and Kulhawy 1999;Shen et al. 2018).The measurement errors are incorporated in this study because of its effects on the characterization of inherent variability of rock/geotechnical parameters.
Combining Eqs. ( 6) and ( 7 From Eq. ( 7), BTS * is a normal random variable with a mean of B + M and a standard deviation of At a given rock formation or site, the statistics of BTS (i.e., B and B ) is used to determine the probability density function (PDF) of BTS * .The BTS is independent of Is (50) data (i.e., DATA I ) for a given set of B and B .There- fore, the P DATA II | B , B , DATA I can be simplified as P DATA II | B , B in Bayesian framework and expressed as: The updated information in terms of the posterior distribution of B and B is obtained based on prior information about BTS, point load test data, and Brazilian tensile test data using Eqs.(4), (9), and (10).Thus, the posterior distribution quantifies the integrated information from the three sources.Using the Eqs.for a specific rock site, equivalent samples of BTS can be generated for probabilistic characterization of the inherent variability of BTS. (9)

Sample simulation from the Bayesian approach
The PDFs of BTS in the updating levels I and II based on the updated information on B and B from Eqs. ( 2) and (3) can be expressed as Eqs.( 11) and ( 12), respectively.
where P BTS| B , B is a normal PDF for a given set of B and B to reflect the normal distribution assumed for BTS; P B , B |DATA I and P B , B |DATA II , DATA I are the pos- terior distribution at updating levels 1 and II, respectively.A large number of BTS samples is generated from Eqs. ( 11) and ( 12) to numerically quantify the updated PDFs of BTS in updating levels I and II, respectively, using Markov chain Monte Carlo (MCMC) simulation.This is consistent with the methodology in many previous studies for solving complicated PDFs such as those presented in Eqs. ( 11) and ( 12) (e.g., Beck and Au 2002;Zhang et al. 2012;Wang and Cao 2013;Juang et al. 2013;Peng et al. 2014;Wang and Aladejare 2016a).Based on the simulated samples of BTS in each updating level, conventional statistical analyses are performed on the samples to evaluate the respective statistics and PDFs of BTS in each updating level.The statistical estimates of BTS summarize the updated knowledge on their respective inherent variability at the site, where Is (50) data and direct BTS measurements from Brazilian tensile tests are obtained.The proposed Bayesian approach provides insights into the evolution of the statistics and probability distribution of BTS as direct and indirect measurements are used.Understanding such evolution is helpful in assessing the quality of information from point load tests and Brazilian tensile tests.It also helps to discover the consistency or otherwise of the information from both tests.

Illustration of the proposed Bayesian approach using real rock site
For illustration, the proposed approach was applied to probabilistically characterize the inherent variability of Brazilian tensile strength at a real rock site in India (Mishra and Basu 2012).Table 1 shows the values of Is * (50) and BTS * obtained for the site.
For the purpose of illustration, 10 data of BTS * were randomly drawn and combined with the 19 Is * (50) as input in the proposed approach from the data in Table 1.Note that it is often the case that for most mining engineering (11) projects, there is a possibility to perform more point load measurements compared to Brazillian tensile strength tests.In addition, a joint uniform distribution with B min = 4.39 MPa; B max = 26.05MPa; B min = 0.91 MPa; and B max = 22.87 MPa was used as the prior distribution in this example, which is consistent with the typical range of BTS of igneous rocks (Aladejare and Wang 2017a).Using the prior knowledge, the 10 randomly drawn data of BTS * and the 19 Is * (50) data from Table 1, a total of 20,000 equivalent samples of BTS were generated in updating levels I and II, respectively, to probabilistically characterize the inherent variability of Brazilian tensile strength at the rock site.Then, the equivalent samples of BTS are compared with the 20 BTS * values from the site.To avoid confusion in this study, only the samples generated in updating level II are shown in this example.

Characterization of the inherent variability of Brazilian tensile strength
Figure 1 shows the 20,000 samples of BTS generated by MCMC simulation in updating level II, based on which a histogram of BTS was constructed.The histogram is presented in Fig. 2, and it includes the 95% confidence interval (i.e., [10 MPa,22 MPa]) ranging from 2.5 to 97.5% quantiles of BTS and estimates of the mean and standard deviation of BTS samples (i.e., mean = 15.24MPa and standard deviation = 2.74 MPa). Figure 3 shows the PDF of BTS estimated from the histogram by a solid line, and the PDF of BTS calculated from direct numerical integration by Eq. ( 12), which is shown by open squares.The two PDFs estimated from the simulated samples of BTS and the direct numerical integration agree well with each other.This validates the Bayesian approach and indicates that the simulated samples in Fig. 1 can sufficiently depict the PDF of BTS after the multi-level updating.By integrating the prior information, Is * (50) data and BTS * data in the proposed approach, the statistics and PDF of BTS satisfactorily characterize the inherent variability of BTS of the granite at the rock site.In addition, Fig. 4 shows the estimated mean and standard deviation of BTS in the two updating levels.It is shown that incorporating BTS * values in updating level II, the mean and standard deviation values of BTS improve and become more consistent with the site BTS.This indicates that integrating the information provided by Is * (50) data and BTS * data, the simulated samples have reduced uncertainty.In addition, the proposed approach is useful in assessing the type of information available during rock parameter characterization and how incorporating information from multiple sources can impact site characterization in terms of reduction in uncertainty.

Effect of assumed measurement error for test data
In this study, the measurement errors M in direct BTS measurements from Brazilian tensile tests are modeled as normal random variables with a mean of M = 1 MPa and M = 0.1 MPa.However, correct characterization of measurement errors for specific rock parameters is difficult because of the possibility of the introduction of different sources of errors  during laboratory testing of rock samples.In addition, it is also difficult to delineate the measurement errors from inherent variability because both uncertainties are usually combined into the scatterness of the rock test data.It is difficult to separate these uncertainties when there is a limited number of laboratory test results, which is often the case during laboratory investigations for the characterization of rock parameters.Assumptions can be made for measurement errors based on engineering judgements, previous experience, or typical values available in the literature.The magnitude of measurement errors assumed represents the extent of belief in the laboratory procedure and testing for characterization of rock parameters, with larger measurement errors indicating a lower degree of belief in laboratory procedures and vice versa.To explore the effects of the assumed M on the probabilistic characterization of the inherent variability of rock parameters, this subsection performs a sensitivity study by assuming four more different values (i.e., 0 MPa, 0.05 MPa, 0.2 MPa, and 0.3 MPa) of M .The four different values in addition to M = 0.1 MPa used in the previ- ous section makes a total of five different values of M .For each assumed M , the proposed Bayesian approach is used to characterize the inherent variability of BTS and estimate the mean and standard deviation of BTS. Figure 5 shows the variation of the estimated mean and standard deviation of BTS against the assumed M by open circles and squares, respectively.
It can be noted that the estimated mean of BTS remains almost unchanged as the value of M increases.This is likely due to the fact that the mean of measurement error remains constant during the characterizations.However, as the assumed M increases, the estimated standard deviation of BTS decreases considerably, which means that the estimated inherent variability of BTS decreases in this example because M remains almost constant.For a given set of site-specific test data (e.g., direct BTS measurements from Brazilian tensile tests), the observed data scatterness is due to inherent variability and measurement errors.When the assumed M is high, more parts of the scatterness of the observed BTS data are attributed to measurement errors during the interpretation of the test data.As can be noted in Fig. 5, the proposed Bayesian approach reflects the degrees of belief in laboratory procedures in a logical way.The implication of the results in Fig. 5 is that when test data are considered perfectly accurate (i.e., M = 0 MPa), the scatterness in measured data during direct measurements is totally attributed to the inherent variability, leading to an overestimation of inherent variability in the BTS.In reality, an increase in M (which indicates a decrease in the degree of belief in the laboratory procedure and testing) means that the magnitude of measurement errors in the scatterness of the measured data is greater, and the estimated inherent variability of BTS decreases as shown in this example.The proposed Bayesian approach is useful in understanding the effect of the degree of belief, which reflects the assumed measurement errors in the laboratory procedure and testing for the characterization of rock parameters.This will be helpful for mining engineers and practitioners because of the difficulty in quantifying the measurement errors involved during the laboratory determination of rock parameters for a specific site.

Simulated rock site
In mining engineering, rock parameters are unavoidably associated with inherent variability.This variability is not known with certainty and is only inferred from prior knowledge like engineering judgement and measured values of rock parameters from direct and indirect laboratory tests.Therefore, a sensitivity study is performed in this section to further validate the proposed Bayesian approach.To perform the sensitivity study, Is *   (50) and BTS * values simulated from a simulated rock site with known inherent variability of BTS are used.In the simulated rock site, the inherent variability of BTS is modeled by a normal random variable with mean B = 12 MPa and B = 3 MPa.BTS * and Is * (50) data are simulated using Eqs.( 7) and ( 8), respectively, while the measurement error is modeled by a normal random variable with M = 1 MPa and M = 0.1 MPa.Ten sets of { Is * (50) , BTS * } values are simulated for different numbers of input data.In this study, four different numbers of input data for both indirect and direct measurements are considered (i.e., n = 5, 10, 20, and 30).The idea of using different numbers of input data is to explore the effect of the number of input data on the output of the proposed Bayesian approach.For a given n value, the 10 sets of { Is * (50) , BTS * } values simu- lated were used together with the prior knowledge in the previous section to simulate 20,000 samples of BTS using the proposed Bayesian approach.Using the 20,000 samples for each case, the estimated mean and standard deviation of BTS were obtained.The procedures were performed 30 times for each n value to obtain 30 values of estimated mean and standard deviation of BTS, and the calculations of 95% confidence intervals using conventional statistical analyses.

Results of sensitivity analysis
The results of the sensitivity analysis performed are plotted in Figs. 6 and 7.In both figures, open circles are used to plot the estimated mean and standard deviation of BTS versus the n value of indirect and direct laboratory measurements used in the sensitivity study.The figures also include the 95% confidence intervals of the estimated mean and standard deviation of BTS represented by error bars.For comparison, the figures also include the true values of B and B which were used to simulate the input data, and represented by dashed lines.As the n value increases from 5 to 30, the shaded part in the middle of the error bars gradually approaches the dashed lines and then plots directly on the dashed lines.This indicates that the estimated mean and standard deviation of BTS are identical to their respective true values that were used to generate indirect and direct input data.At small n values, the values of the estimated mean and standard deviation are scattered to some extent.The scatterness is largely due to statistical uncertainties in the data, which arise as a result of relatively sparse data at such n values.
However, such scatterness gradually disappears as the n value increases, which indicates that the scatterness of the estimated mean and standard deviation of BTS reduces as the number of input data in the proposed approach increases.The estimated inherent variability of BTS obtained from the proposed Bayesian approach becomes more accurate as the number of input data from laboratory investigations used in the approach increases.It can be concluded that the proposed Bayesian approach provides reasonable estimates of the mean and standard deviation of BTS from the integration of both prior knowledge and information from indirect and direct measurements from laboratory tests.

Summary and conclusions
This study proposed a Bayesian approach that integrates engineering experience (i.e., prior knowledge) and measurements from point load tests (i.e., indirect test source in this study) and Brazilian tensile tests (i.e., direct test source in this study) in a logical manner for the probabilistic characterization of the inherent variability of BTS at a rock site.The proposed approach was illustrated and validated using data from a real-life rock site and a simulated rock site.In addition, a sensitivity study was performed to explore the effects of assumed measurement errors on the interpretation of test data during laboratory investigation of rock samples for characterization of rock parameters.The magnitude of measurement errors depicts the degrees of belief in the testing procedures, with increased measurement errors indicating somewhat lack of confidence in the laboratory testing procedures that produced test measurements.The results showed that the proposed Bayesian approach provides a systematic route to incorporate measurement errors into the characterization of site-specific inherent variability of rock parameters in a rational manner.With limited information from indirect and direct test measurements together with prior knowledge, the approach satisfactorily characterized the BTS in both sites.The results of the sensitivity analysis showed that the assumed measurement errors have effects on the inherent variability estimated for a given set of data from laboratory testing.As the number of input data increases, the proposed Bayesian approach reduces the statistical uncertainties and the effects of measurement errors contained in the input data.The insight from the proposed approach will be helpful during the characterization of rock parameters because it makes it logically possible to quantify the inherent variability of rock parameters without lumping it with measurement errors.With the limited number of data often obtained from laboratory testing and investigations into rock parameters, it is impractical to quantify the measurement errors and their effects on the characterization of rock parameters.This further reinforces the usefulness of the proposed approach because such effects are systematically quantified and incorporated into the interpretation of test results for the characterization of rock parameters through the Bayesian approach.
of k I different Is (50) values measured for a rock formation or deposit or site;DATA II = {BTS * 1 , BTS * 2 , … , BTS * k II} is a set of BTS * of k II different BTS values obtained by per- forming Brazilian tensile tests on rock samples recovered from a rock formation, deposi,t or site; the superscripted asterisk (*) indicates the measured value; P B , B is the prior distribution that reflects the information on B and B in the absence of any site-specific test data and is taken as the joint uniform distribution of B and B in this study, which are completely defined by their respective maximum (i.e., B max and B max ) and minimum (i.e., B min and B min ) possible values (Aladejare and Idris 2020); P B , B |DATA I is the posterior distribution of B and B in the first updating level (i.e., updating level I), and it is taken as the prior distribution of B and B in the second updating level (i.e., updating level II); P DATA I | B , B and P DATA II | B , B , DATA I are likelihood functions of the two updating levels that quantify the information on B and B provided by DATA I and DATA II , respectively; and C I and C II are normalizing constants that are inde- pendent of B and B .Substituting Eq. (2) into Eq.(3) gives (2) P B , B |DATA I = C I P B , B P DATA I | B , B (3) P B , B |DATA I , DATA II = C II P B , B |DATA I P DATA II | B , B , DATA I 678 Page 4 of 10

Fig. 5
Fig. 5 Effects of assumed measurement errors on estimated inherent variability of BTS

Fig. 6 Fig. 7
Fig. 6 Mean values of BTS estimated using different numbers of simulated test data B , B |DATA I , DATA II = C I C II P B , B P DATA I | B , B P DATA II | B , B , DATA I + a B Z + a M + T explicitly considers the measurement errors and transformation uncertainty for inferring the inherent variability of BTS based on Is * (50) data.Since the inherent variability, measurement errors, and transformation uncertainty are from different sources, they are assumed to be mutually independent.Therefore, Is * (50) can be taken as a normal random variable with a mean of aμ B + aμ Updating level II uses the direct BTS measurements (i.e., DATA II = { BTS * 1 , BTS * 2, …., BTS * k II }) from Brazilian tensile tests as input in the Bayesian framework.The likelihood function P DATA II | B , B , DATA I of updating level II is used to quantify the information provided by DATA II .
) gives Equation (8) presents a likelihood model that relates B and B to site-specific Is * (50) data.The likelihood model (4) P

Table 1
Mishra and Basu 2012)ts of granite from Malanjkhand Copper Project in the state of Madhya Pradesh, India (fromMishra and Basu 2012)