Journal of Engineering Geology

---

Desiccation cracking commonly occurring in compacted clayey soils typically used as landfill liners can result in poor hydraulic performance of the liner. In this research, a simplified image processing technique was developed in order to characterize desiccation cracking intensity in compacted clayey soils. Three pairs of compacted clayey soils were studied in a relatively large scale experiment to evaluate the effect of geotextile cover on desiccation cracking under real-time atmospheric conditions. Digital images were taken from the surface of soils at certain time intervals for 10 months and were analyzed to determine crack intensity factor (CIF). The key parameter in identification of cracks as accurately as possible was found to be sensitivity. Calibration process was based on using %20 of the images with different crack intensities whose crack dimensions and therefore CIF values have been already measured to compare to program output. A calibration coefficient for sensitivity was accordingly determined based on the average difference between the sensitivity introduced by the program and the actual sensitivity calculated based on an overlaying process. Result of verification of this methodology indicated that it can be reliably used to determine CIF of compacted clay soils in a simple yet accurate manner.

---

The Keshvari watershed is located at south east of Khorramabad city in Lorestan province. This area is one part of the folded Zagros zone based on structural geology classification. By consider the type of geological formations, topographic conditions and its area, this watershed is very unstable and capable for occurring landslide. In this study, artificial neural network (ANN) with structure of multi-layer percepteron and Back Propagation learning algorithm used for zonation of landslide risk. The results of ANN showed the final structure of 9-11-1 for zonation of landslide risk in Keshvari watershed. According this zonation, 23.81, 7.53, 6.49, 18.68 and 43.47 percent of area are located in very low, low, moderate, high and very high risk classes, respectively.

---

Determination of uniaxial compressive strength (UCS) of intact rock is an important mechanical parameter required for many engineering projects. In some engineering projects, for example, well drilling has been accomplished for petroleum. The requirement of deep well to take samples to obtain rock core sample for determination of UCS is a difficult task. On the other hand, determination of this parameter is essential in order to analyze well wall stability and well development program. Therefore, the idea of using drilling cuttings is proposed for determination of UCS. In this paper, in order to develop relationship between UCS and single compressive strength (SCS) 7 block sample of microcrystalline limestone from Asmari formation were used. Then UCS test was performed and uniaxial compressive strength was determined. Next, these samples were crushed and 420 single particles were prepared. Then SCS for each particle was determined. Since the shape of particles affects particle strength, shape of particles was modified.  The total particles used for determination of SCS were spherical. In order to study the effect size of particle, particles with diameters 2, 3 and 4 millimeters were prepared and the SCS for each particle has been determined. With the increase of diameter of particles, the SCS has been increased too. In order to eliminate the effect of size of particles, it is defined variable size and strength and proposed chart between them. Coefficient of correlation between SCS and UCS is more than 0.91 which indicates a high correlation between them.

---

One of the geotechnical hazards in the tunnels under high overburden and high in situ stresses is the phenomenon of rock burst. Rock burst is a typical geologic phenomenon caused by excavation in rock masses. In this phenomenon, because of stress released and explosion in rock masses, they are broken as large and small pieces and are distributed, so that leads to damage of peoples or equipments. Therefore, familiar with this phenomenon and its mechanism of occurrence, is need to analyze this issue. The second part of water supply Karaj-Tehran tunnel with a length of 14 km and about 4.5 m diameter is located in Tehran province. Rock burst analysis has been carried out in the tunnel from kilometer 6 to 9.5 that is critical section because of high overburden (up to 800 m) and presence of faults and crushed zones. In this paper, for predicting rock burst in the critical section of second part of Karaj-Tehran tunnel, four criteria including, Strain energy, Rock brittleness, Seismic energy and Tangential stress criterion are used. Analysis results show that units with high overburden have high possibility of rock burst. 

---

The study area is located in Jiroft district, Iran, and is a part of Sahand-Bazman volcanic zone. There are various landslide factors and the importance of each factor are identified qualitatively, based on previous studies and regional specifications. Three landslides were recognized in the study area using direct method (field work) and aerial photographs interpretation. One of these landslides is located in the vicinity of Mohammad Abad of Maskoon Village. The aim of this study is landslide hazard mapping using two integration methods that includes Fuzzy Logic and Hybrid Fuzzy-Weight of Evidence (Hybrid F-W of E). The obtained results of maps from both methods, show a good agreement especially in introducing  high hazard regions. The hybrid method is based on the occurred landslide points and is more rigorous, so hazard regions delineated by this method occupy smaller areas than the areas introduced by fuzzy model. Therefore, hazard maps resulted from Hybrid and Fuzzy methods, can be considered as minimum and maximum limits of landslide hazard in the area, respectively. 

---

In this study using finite element procedure was used to simulate the dynamic behavior of reinforced soil walls, to evaluate their dynamic response on all types of deformation modes, different mechanisms of failure detection and identification of parameters in each of the modes and the mechanisms. Detailed numerical modeling, behavioral models and materials were described and Dynamic response of the physical model has been validated experimentally. Parametric study has been of the wall height of 5 meters by the effective parameters such as hardness, length to height ratio, the vertical reinforcement, wall height, and acceleration inputs. Three modes of deformation were observed. The study showed that occur bulging deformation mode while the use of flexible reinforcement and occur overturning deformation mode while the use of stiffness reinforcement. Stiffness reinforcements have the most effective in changing the type of deformation. Length to height ratio of reinforcements has the minimum effective in changing the type of deformation.

---

Various types of numerical analyses such as   Finite Element Method, Boundary Element Method and Distinct Element Method, are used in rock mechanics and in engineering practices for designing rock structures such as tunnels, underground caverns, slopes, dam foundations and so on. In this paper, the results of back analysis of Koohin tunnel which is located in the first section of Qazvin-Rasht railway have been presented. The main purpose of this paper is to perform the back analysis of the mentioned tunnel with the use of numerical models. For modeling the tunnel, two different sections of 30+150 km and 30+900 km are analyzed with FLAC 2D software.  To perform back analysis the suitable interval of geomechanical parameters according to the tests which were performed on the core drillings has been determined. With the use of direct method in back analysis, the errors of models have been corrected in several steps and finally the geomechanical parameters in 30+150 km station (Elastic Modulus = 0.3 GPa, Cohesion = 0.21 MPa and Internal Angle of Friction = 34^°) and in 30+900 km station (Elastic Modulus = 0.3 GPa, Cohesion = 0.21 MPa and Internal Angle of Friction = 35°) have been achieved. The geomechanical parameters which obtained from back analysis are completely in the chosen interval and compliance with the results of tests which performed on core drillings. On the basis of  geomechanical parameters obtained from back analysis with the parameters which used in the design of the tunnel, the tunnel design and the structure method were confirmed.  

---

Marl rocks are from weak rocks which cause some problems due to high swelling and efflorescence capability, low resistance and durability in construction of engineering structures. Creation of these problems is due to inadequate recognition of engineering geology properties of these rocks. Hence, in this research for determination of the physical and mechanical parameters of marl rocks of Safa dam site, Aterberg limits, density, porosity and moisture percent, uniaxial compressive strength (UCS), direct shear, swelling and three axial compressive strength tests and X-ray diffraction (XRD) and X-ray fluorescence (XRF) analyses were performed on available samples from excavated bores. Results of experiments indicate that the type and contents of minerals and percentage of calcium carbonate in marl rocks, especially clay minerals are effective factors on engineering geology properties of these rocks

---

Determination of stress in earthfill dams is one of the most important parameters in dam safety studies. Stress monitoring can be done using total pressure cells which are typically installed during construction. The cell is installed with its sensitive surface in direct contact with the soil to measure total stress of soil and in combination with piezometers to measure pore-water pressure acting in the soil mass. Total pressure cells needs to be installed with care to get reasonable measurements. However, measurements are often incompatible with the theoretical predictions such that pressure cell results usually have some inaccuracies. There are several parameters effecting pressure cell errors. However, in the present paper it is only focused on the height of embankment and the duration of dam construction. For this purpose, a case study, namely Alborz embankment dam located in northern part of Iran was studied. It is an earth dam with clay core with a height of 78 m. Using the monitoring data and considering the effect of embankment height and construction period parameters, a model is presented to predict the pressure cells error with Gene Expression Programming (GEP) procedure by GeneXProTools 4.0 software. The computed coefficient of correlation (R²) for the proposed model is 0.98 showing a good agreement with the monitoring data. The obtained results indicate that the ratio of height difference to time difference for Alborz dam has a significant role in dam pressure cells errors

---

Comprehensive laboratory tests were performed to assess fatigue behavior of Alvand monzogranite rock subjected to uniaxial cyclic loading. A series of static loading tests was done to obtain the required data for the fatigue tests. Three maximum load levels (85, 90, 95% uniaxial compressive strength (&sigmac)) at amplitudes 70% were used with 1Hz cyclic loading frequency. The results indicated that maximum stress level significantly influenced fatigue behavior of this rock. It was found that fatigue life decreases in a power function with increasing maximum stress level. Accumulative fatigue damage process shows three stages of behavior including crack initiation phase, uniform velocity phase and acceleration phase. Fatigue damage process were analyzed according to axial and lateral maximum and minimum strain, tangent and second modulus, toughness and hysteresis energy in both loading and unloading conditions. Among these parameters, lateral strain, axial strain and second modulus show the best three-stage fatigue damage behavior. Also, it should be noted that most of the cracks generated in parallel to loading direction and lateral strain are affected by more than axial strain.  

---

Today the effects of grouting are usually confirmed by the results of permeability tests but this method is not enough to show the changes in mechanical properties of rock masses. Although many investigators use the in situ tests for evaluation of rock mass mechanical property improvement. But this tests are time consuming and expensive. Grouting reduces the permeability and improves the condition of joints and ultimately increases the rate of rock mass classification in rock engineering. So with measurement of rock mass quality index values (Q-value) in cores obtained from grouted boreholes, the efficiency and success in improving the mechanical properties of rock mass can be showed. This paper for first time introduces Q-logging as a simple method to assess the impact of grouting in improvement of the rock mass quality. Here in, the results of Q-Logging in trial injection panels in the Bakhtiary, Bazoft and Khersan II dams have been examined. The deformation modulus were calculated from the Q-Logging for before and after of grouting. Results show that there is a good agreement between calculated rock mass parameters based on the Q-Logging method and the measured from in-situ test in the studied site. This agreement confirms the efficiency and applicability of the Q- Logging method for assessment of grouting success as well as the estimation of the rock mass parameters in grouted areas. Also it has been shown that the deformation modulus in weak rock mass with low quality has been more improved than rock mass with beater quality.

---

One of the tests that is used for the characterization of soil abrasivity, is LCPC test. LCPC test device is designed for measuring the wear particles as small as pebbles (4-6.3 mm). In this study, some of the most important abrasive minerals were collected from different parts of Iran for analyzing the effect of the geological parameters on the ability of abrasive minerals. Firstly, amount of index minerals abrasivity is measured according to three standards of AFNOR P18-553, AFNOR P18-579 and AFNOR P18-560 that are the preparation of samples for testing, procedure of laboratory tests and analysis of grain size with laboratory sieves. The effect of geological parameters affecting the wear rate of the sample, including five parameters of shape, size, angularity and saturation rate of the environment, has been studied. The effectiveness of these parameters on the abrasivity of samples are studied according to NF ISO 5725 relating to usage of statistics, the accuracy of test method, the repeatability and the ability to reproduce a standard way of testing within laboratory (based on classification index X 06-041). Finally, after ensuring significant effect of these parameters on the abrasivity of minerals by help of SPSS, abrasivity rates for types of minerals that have the hardness below 7 in the Mohs hardness scale, have been predicted.

---

One way of reduction of leakage from beneath of earth dams is using of one contact clay layer with very low permeability and intermediate to high plastisity and connectig it to core of dam. Since, most of fine-grained soil in environtment of dam have low plastisity and preparing it from another place is not economic, use of bentonite in order to improvement of engineering characteristic of borrowed clay is suitable way.
In this search effect of bentonite on geotechnical properties of fine-grained soils with low plastisity are evaluated. Results of this research show that hydraulic conductivity, consolidation coefficient, dry density, colifornia bearing ratio (CBR),.....are decreased with increase in bentonite content but optimum moisture , Aterberg limits , cohession and so on are increased with bentonite addition. Finally, with analysis of obtained result, optimum percent of bentonite is offered in order to improve of engineering properties of used clay in contact region.

---

The in situ measurement of discontinuity geometry of rock mass exposures is a time consuming and sometimes hazardous process. Moreover, a large proportion of the exposure is often inaccessible. Thus, a fast and safe tool is required in order to acquire the information which characterizes the geological/structural regime. Digital image processing techniques provide the necessary tools for realizing this goal. This paper presents a methodology for automated discontinuity trace detection in digital images of rock mass exposures. In this study at first based on difference in gray level discontinuities with the face, fracture traces detected in images of rock face. Then some parameters of discontinuities geometry such as spacing, linear joint density, persistence, trace angle of joints and value of RQD are obtained. The Automated discontinuity geometry analysis system including: 1- Providing a digital image from rock face 2- The pre-processing on the images 3- Detection of edge or joint traces by the canny detector 4- Description of the edges using line detector by the Hough transform 5- The joint sets estimation using fuzzy methods and 6- Description the rock mass geometry properties.

---

Direct observation and experience of past earthquakes together with modeling carried out by researchers, has shown that ground motion acceleration and frequency is affected by the nonlinear behavior of site soil. In the process of assessing the seismic response of structures and lifelines, it is essential to understand the nonlinear behavior of the soil and how it can affect the results. In this paper, the nonlinear behavior of Urmia's subsurface soil is studied by performing one dimensional nonlinear site response analysis in time domain. Artificial acceleration time histories that were synthesized based on the result of seismic hazard analysis, conducted over three return periods, are used as input motion. Spectral acceleration at the ground surface is compared with those calculated for seismic bedrock, and spectral acceleration amplification curves are obtained. These curves show that, the amplification is greater in the central and eastern regions of the city than those for other regions of the city because of a deeper soil profile. The results show that the maximum amplification for higher return period is smaller because of greater soil nonlinear behavior

---

One of the most useful procedures in soil stabilization is lime. Soil improvement using lime is a quick and simple approach which could be included in large and small projects. The objective of soil ‘Improvement’ with quicklime is to achieve an immediate reaction, which significantly strengthens the soil due to the removal of moisture and a chemical change in clays. In order to do a parametric study on the influence of the lime on shear strength preparing the samples is important. In this paper, in addition to considering a method of samples preparation, the effect of lime content, water content and processing time on the shear strength of clay using direct shear test is investigated. The results indicate that the method of samples preparation is effective and is identified that there is an optimum lime and moisture content which maximize shear strength.

---

The objective of this paper is to investigate the bearing capacity of strip foundations next to sand slope. A series of laboratory model tests has been carried out and a new correlation coefficient to estimate the bearing capacity of shallow foundations near slopes is presented. The sand layers were prepared in a steel test tank with inside dimensions 500 ´ 200´ 250 mm. After vertical loading, the applied load and displacement of foundation were recorded and stress-settlement curve is drawn. Finally, the load at which the shear failure of the soil occurs is recorded as ultimate bearing capacity of foundation. The test sand used in this study was Babolsar sand with relative density of 50%. The relative performance of different distance of foundation from the edge of slope and inclination angle of slope are compared using same quantity of soil properties in each test. The results indicate that with increasing distance from the edge of the slope, bearing capacity increases linearly. Also with increasing slope angle, the bearing capacity has declined linearly

---

The present study aims to employ intelligent methods to predict shear wave velocity (V_s) in limestone. Shear wave velocity is one of the most important rock dynamic parameters. Direct determination of this parameter takes time, cost and requires accuracy as well. On the other hand, there is no precise equation for indirect determination. This research attempts to provide some simulations to predict V_s using the information obtained several dams located in Iran, using different approaches, including adaptive neuro-fuzzy inference system (ANFIS) and gene expression programming (GEP). 136 datasets were utilized for modeling and 34 datasets were used for evaluating its performance. Parameters such as Compressional wave velocity (V_p), density (g) and porosity (n) were considered as input parameters. The values of R² and RMSE were 0.958 and 113.620 for ANFIS, where they were 0.928 and 110.006 for GEP respectively. With respect to the accuracy of the intelligent methods, they can be recommended for future studies

---

This paper describes triaxial compression tests conducted to determine the effect of fiber inclusion on stiffness and deformation characteristics of sand-gravel mixtures. Tested soil was a mixture of Babolsar sand from the shores of the Caspian Sea and Karaj River gravel. Portland cement was used as the cementing agent and fibers 12mm in length and 0.023mm in diameter at 0%, 0.5% and 1.0% were added to the mixtures. Triaxial tests were performed on saturated samples in consolidated drained and undrained conditions at confining pressures of 100, 200 and 300 kPa. Deviatoric stress-axial strain, volumetric strain-axial strain, pore pressure-axial strain curves with deformation and stiffness characteristics were investigated. Tests results show that fiber addition increased peak and residual shear strength of the soil. Fiber addition resulted in an increase of the maximum positive and negative volumetric strains. In undrained condition, fiber inclusion caused increase in initial positive pore pressure and final suction. It has also been observed that fibers decreased initial tangent stiffness of the cemented sand-gravel mixture.

---

In contrast with Mobility Factor (MF) and Risk Assessment Code (RAC) indices, I_R attributes a risk share to metal species bound to reducible and oxidizable phases which are totally neglected in both of the two above-mentioned indices. In other words, besides the absolutely mobile fractions, the potentially mobile ones are also regarded in risk evaluation process elaborated by I_R. The different structure of the newly-developed index may prevent risk level underestimation especially in case where a remarkable percent of bulk concentration is accumulated within reducible and oxidizable phases. The independency of the index value to the bulk concentration makes it possible to discuss the potential risk in different levels of bulk concentration. Furthermore, the index capability in indication of risky pollution, regardless of the pollution source type, may prevent the probable misleading caused by distinct separation of bulk concentration into geopogenic and anthropogenic portion