Journal of Engineering Geology

---

Hoek and Brown suggested a method to estimate the strength and deformation modulus parameters of rock masses. The method was then widely used in rock engineering designs. In such designs, the mean values of Hoek and Brown parameters are often used which are not proper values due to the variability of rock mass properties within a great range of values. In such cases, probability analysis of rock mass properties is highly important. The geological strength index is one of the most important parameters in Hoek and Brown equations. Determination of this parameter includes greater uncertainties than determining other parameters. In this paper, based on the results of rock mechanical tests carried out on rock samples of Gol-Gohar iron ore mine, and the required field surveys, the sensitivity of rock mass geomechanical properties on the type of the statistical distribution function of the geological strength index in statistical analysis of these parameters using Monte Carlo simulation method was investigated. The results showed that the sensitivity of Hoek and Brown equations to determine different rock mass geomechanical parameters varies as the type of the statistical distribution function of the geological strength index changes. The sensitivity of geomechanical parameters such as internal friction angle, cohesion, total strength and rock mass modulus on the type of the statistical distribution function of the geological strength index is much less than parameters such as uniaxial compressive strength and tension strength of rock mass. The greatest variations based on changes of the type of the statistical distribution function of the geological strength index are less than 5% for the internal friction angle, cohesion and total strength, less than 10% for the modulus, and less than 25% for the uniaxial compressive strength and tension strength.

---

Marly rocks of Abtalkh formation were classified by Q, RMR, RSR and RQD rock mass classification systems using 222 meters logs from exploratory boreholes in Doosti dam site. The results show that the RMR is the most suitable method for classification of studied rock masses and has highest correlation coefficient with RQD. The validity of different Q-RMR equations was studied using error ratio (ER). Cameron et al. (1981) and Morno (1982) equations have lowest ER and highest validity for studied marlstones. Bieniawski (1989) and Cameron (1981) relationships are lower and higher limits of equations for marly rocks respectively. 

---

Combination of Adoptive Network based Fuzzy Inference System (ANFIS) and subtractive clustering (SC) has been used for estimation of deformation modulus (E_m) and rock mass strength (UCS_m) considering depth of measurement. To do this, learning of the ANFIS based subtractive clustering (ANFISBSC) was performed firstly on 125 measurements of 9 variables such as rock mass strength (UCS_m), deformation modulus (E_m), depth, spacing, persistence, aperture, intact rock strength (UCS_i), geomechanical rating (RMR) and elastic modulus (E_i). Then, at second phase, testing the trained ANFISBSC structure has been perfomed on 40 data measurements. Therefore, predictive rock mass models have been developed for 2-6 variables where model complexity influences the estimation accuracy. Results of multivariate simulation of rock mass for estimating UCS_m and E_m have shown that accuracy of the ANFISBSC method increases coincident with development of model from 2 variables to 6 variables. According to the results, 3-variable model of ANFISBSC method has general estimation of both UCS_m and E_m corresponding with 20% to 30% error while the results of multivariate analysis are successfully improved by 6-variable model with error of less than 3%. Also, dip of the fitted line on data point of measured and estimated UCS_m and E_m for 6-variable model approaches about 1 respect to 0.94 for 3- variable model. Therefore, it can be concluded that 6-variable model of ANFISBSC gives reasonable prediction of UCS_m and E_m.

---

Introduction
Safety and sustainability of infrastructures which were placed in or on rock mass mainly control by geometrically size distribution and physical and mechanical characteristics of rock blocks that is created by intersection of discontinuities. hence identification of rock blocks has a key role in mechanical analysis and hydraulic behaviour of jointed rock mass. Detection process of blocks have many applications in rock mechanic which could be referred to their use in the numerical methods like discrete element method or in analysis of continuous deformation of discontinuities. As pioneer researchers, Goodman and Shi, Warburton and Heliot could be known as leaders in the field of diagnosis of rock mass blocks. Warburton provides a method based on geometric parameters of rock mass and developed a software based on it. Warburton in his work assumed discontinuities as parallel and infinite. In the earlier works, discontinuities were considered as infinite panes. So, just convex blocks were distinguishable. Concave blocks were diagnosis in more detailed researches that is created by finite discontinuities. Basically, methods based on finite planes was classified into two branches. Aforementioned branches were based on blocks detection based on topology concepts and assemble of block elements. Lin at al. presented detection method that assumed discontinuities as finite planes and worked based on topology theory. This method could realize convex and concave blocks of rock mass. Ikegawa and Hudson, Jing presented the similar methods using more accurate process. Sharma et al. presented an equation for calculating the volume of rock blocks in their work. Ferreira provided a method based on graph theory which is better than other method considering time and complicity. Based on this method, firstly vertices were detected in two dimensions and then created a graph based in vertices and edges which in next step constitute polygons that are form in two-dimension blocks. In the present research, it is developed high-speed algorithms to identify the blocks. New method was developed in MATLAB software that by assuming infinite discontinuities and inclusion of a set of joints. we have identified created blocks and calculated their volume and at last block volume histogram were draw that paves the way to obtain their distribution function.
Material and methods
Infinite planes are used to simulate of discontinuities.in this study, each discontinuity is represented by a plane in a three-dimensional Euclidean space. To identify the block, a certain volume of rock mass space should be considered as study region. The studied volume is called domain. By the intersection of discontinuity planes in space, rocky blocks are created in the domain. First, vertices should be recognized at first as first step in block detection. Then, edges are diagnoses and after that it's time to specify the polygons and finally, polyhedron or blocks are obtained by joining edges together. Each vertex in space is created by the intersection of three nonparallel planes. In fact, the vertex is the interface of three planes in the Euclidean space. The next element in the block metric process is the diagnosis of the edges or the blocks' edges. All edges are sections on the lines which created by the intersection of the planes in space. first the parallel vector of all the lines resulting from the intersection of the pair of planes is obtained.
After detection of edges, it’s time to identify polygons that form key element of blocks. Each polygon of a block is formed from their constituent unit. In this step, polygons belong to each discontinuity plane is identified separately. Some edges are determined that are start from the end of selected edge between other edges. In this state, if there is just one edge, that edge is record as the next edge of first polygon. If there is more than one edge from the edge of the selected edges, the angle is calculated between each possible of end edge with the selected edge.
In the next step, it’s time to diagnosis polyhedrons that have created by discontinuities intersection. In the previous step, possible polygons were obtained for each discontinuity. In this stage, it is used the principle which is designed this algorithm that two polygons that formed a block have a common edge. So, the first polygon of first discontinuity is consider as first polygon of first block to recognize block.
Results and discussion
According to the developed algorithm, MATLAB software was used to model the discontinuities. The computational and graphic capabilities of this software have created a lot of attractions for most researchers to use its potential. The strengths of this software are high computing power with its graphical accuracy. The code developed in MATLAB is called RockBlock2 that is designed using a graphical user interface (GUI) to make it easy to use. To illustrate how the program works, there are 29 discontinuities given to the program. The program first takes the dip and dip direction of discontinuities along with the desired point on it and calculates the parameters that make up the equation of discontinuity planes.
Input data is stored in a separate Excel file that was previously introduced to the program. In the next step, the program attempts to identify the vertices. The program stores the coordinates of each corner, with the assignment of a number to it, in the matrix of the corners, which is in fact the Excel file that was previously introduced to the program to use in the next steps, after recognizing vertices on the area.
Identifying the edges is the next step that the program done. At this stage, the program begins to identify each single edge using the data from the previous step that means the coordinates of the corners and the algorithm defined.
The coordinates of the beginning and end of each edge along with its number are stored and maintained in the edge matrix in the Excel file format. In the stage of identifying the polygons, the polygons are formed by joining the edges together. This matrix is a special matrix that its matrix matrices are matrix itself. The matrix of polygons is a row matrix; whose number is the number of discontinuities. Because, as it mentioned in the chapter of the algorithm, the polygons are found by separation of discontinuities. Therefore, each column of the polygons matrix is consisting of faces that are on a certain discontinuity.
The next step begins the process of identifying the blocks, or the same polygons by the program. At this step, the program starts the identification process using the features found in the previous step and the algorithm defined for it. At this stage, the identified blocks are stored in the blocks matrix. By identifying blocks, the program calculates the volume of each block and finally draw its volume histogram. In fact, a volume histogram is presented to illustrate how the block volume is distributed. Obtaining the distribution of blocks or, in other words, achieving a block probability distribution function is an essential step in the behavior of rock mass. Because one of the most important consequences of the presence of discontinuities is the fragmentation of the rock material under the block intervals. By having the block distribution function, it is possible to produce a blockbuster method using random methods, such as Monte Carlo, and to analyze it in various and arbitrary modes.
Conclusion
To identify and study the rocky blocks created by discontinuities, a hierarchical algorithm was designed and developed in MATLAB software. This algorithm identifies and records blocks, consisting of blocks, edges, and facets of the blocks forming components, including stone blocks. This algorithm, which is written for user-friendly ease with the use of graphical coding capabilities, shows a very fast performance using the parallel computing power of MATLAB software. The developed code calculates the dip and dip direction of discontinuities using the geometric properties, and calculates the blocks created in three dimensions and calculates their volume. This histogram code displays the calculated volumes.
The results show that the developed code with its fast performance, while identifying the blocks, calculates and records their volumes without errors. The ability to display the step-by-step process of identifying blocks is one of the clear features of this code. Information about edge is also records and is available for auxiliary applications. Histogram of block volume is one of the most important results of the developed code, which can have different applications.
Identification of created rocky blocks is used both in the stability analysis and rock mass simulations such as Discrete Fracture Network modeling. Determination of block volume distribution function which is done using histogram is one of the most important uncertainties in three-dimensional rock masses behavior that can play a key role in optimizing the design of structures involved in rock mass. Therefore, considering the key role of blocks volume, identifying and calculating block volumes and, consequently, plotting their histogram and determining the distribution function governing them, has a key role in the static and dynamic analysis of rock base structures. ./files/site1/files/123/1Extended_Abstract.pdf

---

A jointed rock mass presents a more complex design problem than the other rock masses. The complexity arises from the number (greater than two) of joint sets which define the degree of discontinuity of medium.  The condition that arises in these types of rock masses is the generation of discrete rock blocks, of various geometries. They defined by the natural fracture surfaces and the excavation surface. Stability problems in blocky jointed rock are generally associated with gravity falls of blocks from the roof and sidewalls. Whereas for block defined in the crown of tunnel,the requirement is to examine the potential for displacement of each block under the influence of the surface tractions arising from the local stress field and the gravitational load, in this paper various types of wedge formation in the crown of tunnel due to intersection of joint sets with various dip were examined. The state of stability of the wedge was then assessed through the factor of safety against roof failure. Following that the formed wedges in New York city and Washington D.C tunnels crown were investigated with limiting equilibrium analytical method and by use of Hoek and Brown failure criterion. The obtained results from analytical method corresponded with field observation.
 

---