Journal of Engineering Geology

---

 Knowing the engineering geological characteristics of carbonate formations is necessary for database. In this research, using petrological study and mechanical tests on 5 types of Ilam-Sarvak formations limestones in Khorramabad city, their engineering geological characteristics were determined and the relationship between physical and mechanical properties have been analyzed. IBM SPSS Statistics (version 19.0) software was used to determine the required relations. The relations have high correlations. Based on the studies on this of thin sections, rocks are characterszed as biomicrite. Limestones of Ilam-Sarvak formations have high hardening and low porosity. These rocks are in medium to high density, very resistant durability index, medium to high UCS and high point load strength category. The rocks are also impermeable. Based on the UCS, modulus ratio of the intact rock, the limestones are CM and CL. According to solubility test, the solution velocity constant was 1.39×10^-6 m/s.

---

The collapsible soils are usually known as soils with open structure and weak bonding between particles. The behavior of such soils is similar to very loose soils. These types of soils, when saturated without any changes in loading or subjecting vibratory loads, experiences huge settlements. The present research deals with investigation on collapsible soils located in the North East of Mashhad. The results of laboratory and in situ tests show that collapsible soil in this region was very sensitive to the increasing of moisture content. This means that an increase in soil moisture content, significant excess settlement occurs during a short time. This indicates that the soil in this region suffers from high potential collapsibility. The huge soil settlement will lead to the stability of existing structures to be at risk. It is, therefore, necessary for the collapsible soil in the region to be improved. For soil improvement, many techniques including moisture mitigation and soil replacement or compaction may be employed. Also stabilization of soil with lime, cement or coarse aggregates are practical methods. Which The results of the present research indicate that stabilization of soil lime is the most appropriate method for increasing bearing capacity of soil and reducing structural settlement.

---

For prediction of exact effect of soil temperature, soil water retention, unsaturated hydraulic conductivity and water flow of soils on seep development, information about soil hydraulic, physical properties soil temperature is necessary. In this study, using geophysical and thermal methods the hydraulic properties and thermal conductivity of unsaturated porous media was estimated. Infiltrated water temperature was 50° C and the temperature data used in inverse solution technique was recorded by 2 thermal sensors during 3255 sec. One of thermal sensor was installed at a depth of 50 cm and the soil hydraulic and thermal parameters were estimated using 48 collected data of the injection water temperature variations. In this research 121 geoelectrical sounding with the Wener, Schlumberger and polar-polar arrays configuration were used. Evaluation and interpretation of the sounding curves show that the water table is located at depth of 75 meters. Top alluvium layer (unsaturated zone) is composed of three layers. In this study, to investigate the distribution and diffusion of pollutants in the unsaturated zone, the brine concentration of 20 g/L was used. Wener and polar-polar arrays were used to detect and image the contaminants plume. The polar-polar configuration shows less spreading of contaminated plume in both vertically and horizontally direction than the Wener array.  Particularly by the Wener configuration the plume depth is more underestimated than the real value. Based on water injection visibility in the vicinity of the trenches and observing the water seepage at the trench face, the real velocity of the water at direction of 35 degrees relative to the horizon was measured equal to 0.159 m/h. Horizontal and vertical contaminant plume expansion coefficient was determined equal to 1.4. Evaluation of thermal and hydraulic properties of unsaturated porous media was done by HYDRUS-2D software. In this study five hydraulic parameters &thetar, &thetas, &alpha, n and Ks and five soil heat transfer parameters including the &lambdaL, &lambdaT, Cn, Co and Cw were included in the parameter estimation process. After implementation of the parameter estimation (adjusted for ​​observed and predicted temperature values) Coefficient of Determination  R² was equal to 0.97 and the objective function value was equal to 11.5. Computational error in the mass balance was about 67 percent. 

---

This paper presents a case study of the instability mechanism, to verify and reinforcement method adopted construct collapsed zone of Sabzkuh water conveyance tunnel in southwest Iran. The instability problems were encountered during tunnel excavation due to the failure, changes in stress field lead to deformation causing dilation and increasing the permeability of sand and gravel layers, local fault gouge zones, landslide and in turn significant reduction in shear strength and collapse in tunnel. IPE Arch Support Technique (IAST) was, used for T1 part of Sabzkuh tunnel zone in order to reinforce the ground around tunnel and to cross the zone falling. In this study, Finite Element Method was employed for the quantitative reinforcement effect with deformation modulus of ground, IPE length and size. As a result, the settlement increases as length increases and decreases with the increase of the deformation modulus of ground and IPE size.  

---

The main objective of this research is to investigate the mechanism of cement-clay-heavy metal contaminant interaction from micro-structural point of view. To achieve this objective series of batch equilibrium and XRD experiments were performed. The results indicate that the addition of 10% cement not only stabilizes the soil, but also at 250 cmol/kg-soil of heavy metals causes 130% increase in heavy metal retention. Furthermore, the XRD analysis shows that in solidified samples with less than 10% cement, the main reason for reduction in peak intensity of clay fraction is due to the presence of heavy metals. However, as the percentage of cement increases, the clay solubility is the main reason for reduction in peak intensity of montmorillonite in XRD test

---

The Dynamic Probing is an important test in site investigation and geotechnical studies. It is used for determination of situation and specification of soil layers when the depth of exploration is moderate. In our country, in the last few years test has been widely used to replace the standard penetration. In this paper, the experiences to use different various types of dynamic probing rigs and the errors that can occur in these experiments are discussed. Then the accuracy of this test is studied to evaluate the specification of fine-grained soils. The most important innovation of this research is the proposition of the new correlations between cone dynamic resistance and undrained shear strength and so compaction percentage in fine-grained soils. The paper encourages the wider application and further development of this test for site investigation in fine-grained soils.

---

The Kandovan village is one of the tourist attractions in East Azarbaijan province of Iran, whose rural houses were excavated within the cone-shaped pyroclastic rocks (in terms of local called keran) several hundred years ago. The present paper discusses the role of engineering geological properties of Kandovan pyroclastic rocks. Kandovan pyroclastic rocks have low resistance against weathering and erosion because their components are plagioclase minerals and pumice fragments with low resistance, welding, sorting and high sphericity and rounding. Although weathering and erosion along existing joints and fractures is the most important causative agent of cone-shaped forms but there is the possibility of further damage of rocks due to continuing these processes. High porosity of rocks has caused that their high capacity for water absorption. High water absorption percent increased sensitivity of rocks against expansion and contraction by freezing-thawing and wetting-drying cycles and low hardness and low their internal strength caused the rocks weathered and disintegrated due to environmental factors. Furthermore, the weak texture of the pyroclastic rocks have caused easy erosion of those by surface waters and wind.

---

Determination of Faults activity rate is among the most important parameters of evaluating faulting hazards. In this paper, active faults on region of Shahid Rajai dam of Sari with radius 100 km based on fractal dimension of faulting and earthquake as well as evaluation of slip rate were classified and those Activity rate were assessed. In order to determine of fractal dimension of faulting and earthquake, Box-counting method was used. For estimate of slip rate beginning the seismicity parameters (a&b) of study area was estimation. Then these parameters were normalized for each fault. Based on the existing relationships and having a&b for each fault, the seismic moment of fault was calculated. Finally according to extant relations for evaluation of seismic moment rate, slip rate of each fault was determined and the faults of study area were classified accordingly. By grading based on fractal dimensions, the faults of North- Alborz, Damghan and Garmsar have been the most active faults in the study area during the last 100 years and according to evaluations of active rate of faulting and earthquake based on fractal dimensions, generally set in category BD and their activity approved. The faults with very low slip rate and with long return period of earthquake, are possible causes of occurrence large earthquakes (856 AD) Gomes and inducement fault namely Damghan Fault is an example of these faults. The faults of Astaneh, Rameh and Cheshmeh-ali have been low slip rate and thereupon be able to develop high seismic moment rate. According to earthquake events in privacy of New faults introduced (Khorram-abad, Majid,…), the high activity rate of these faults seem logical

---

In order to assess mineralogical composition influence on durability and strength of rocks, four samples of hornfelsic rocks were selected from southern and western parts of the city of Hamedan, west of Iran. These rock samples were subjected to mineralogical, physical and mechanical tests in the laboratory. Also, they were evaluated in 15 cycles of slake-durability testing in different pH of sulfuric acid solutions and XRD analysis. Based on the results, the type and amount of minerals, their density and hardness had an influence on the uniaxial compressive strength and the slake durability index of tested rocks. That means, presence of non-dissolved minerals such as graphite in studied rocks, decreases Unconfined Compressing Strength (UCS) but increases the slake durability index. The results of slake-durability test indicated that weight loss of the samples at initial cycles was found to be higher than the end cycles. Also, in these samples, initial minerals in the fresh samples were not exchanged by secondary minerals such as clay minerals. Therefore Hamedan hornfelsic rocks are approximately resistant when were put under accelerated chemical weathering and degradation in the laboratory and natural chemical weathering.

---

In this research attempts were made to estimate the in-situ stresses acting on a hydrocarbon reservoir based on routine activities of acid injection in oil reservoir. It was found that the relation between the re-opening pressure of fracture and principal in-situ stresses can be estimated using rock mechanics equations for the circular underground cavities. An appropriate relation between the maximum and minimum horizontal principal in-situ stresses and reservoir parameters such as permeability, reservoir pressure, Young’s modulus, acid viscosity, injection flow rate and etc., was developed by using the well-known Darcy equations for fluid flow in porous media. Accordingly, knowing the flow rate of acid injection during well operations and some other reservoir parameters, in-situ stresses may be estimated. The method was then successfully applied to a large carbonate reservoir as a case study in south-west of Iran. Maximum and minimum effective horizontal stresses were calculated by employing the presented method. 

---

This paper aims to present an experimental and numerical study on the effect of wetting-drying cycles and plasticity index on the California Bearing Ratio (CBR) of lime stabilized clayey soils. The numerical analysis was carried out based on finite element method for comparison between results of experimental and numerical studies. Three clays with different plasticity indices were mixed with various amounts of hydrated lime and compacted at optimum water content. The CBR tests were conducted to the soils and admixtures after specified curing time and various numbers of wetting-drying cycles. The experimental results indicate that addition of lime content up to 4% causes significant increase in the CBR values. Based on the obtained results the CBR decreases during the wetting phase and increases during the drying phase of each cycle. After 3 cycles the CBR values of lime stabilized clayey soils are increased. Also, for stabilized clays by increasing the plasticity index, the CBR values resulted by increase of lime content are decreased. The comparison between numerical and experimental analyses indicates a good agreement between results.

---

The pseudo-static analysis is one of the conventional methods in embankment dams design and International Commission on Large Dams (ICOLD) suggests using this method before ultimate dynamic analyses. In this research, the static, pseudo-static and dynamic analysis of Masjed Soleyman embankment dam was performed. Using dynamic and pseudo-static analyses results, the safety factor of critical sliding surface was calculated. Permanent displacements of critical sliding surface were evaluated by New mark method and the calculated safety factor was compared. Based on the comparison results in different water levels of the reservoir and by introducing a new equation, the variable horizontal acceleration coefficients in height of the dam body were calculated. Finally, the obtained horizontal acceleration coefficients were compared with the other criteria introduced in different embankment dam's design codes. Totally, the results indicate that the proposed method leads to a larger horizontal acceleration coefficient in higher parts of the dam body.

---

One of the best approaches to reduce transportation problem is to use the underground tunnels. Therefore, Niayesh highway tunnel was performed by the New Austrian Tunnelling Method (NATM) in the northern part of Tehran and it includes north and south tunnels. The excavation of tunnels and other underground structures cause considerable changes in local stress conditions around structures that lead to surface settlement. In this research, surface settlement has been studied for five sections (CS-1 to CS-5) by empirical methods, numerical analysis and actual settlements. For the empirical and numerical methods, O’Reilly and New (1982) method and also finite element method (PLAXIS_2D software) have been used, respectively. On the basis of the obtained results, the numerical method in all sections (except section 3) is in agreement with the actual settlements. While, empirical methods have estimated the settlements more than actual values in those sections.  Also, the achieved results from the aforementioned methods show that the maximum settlement due to tunnel excavation is more than allowable settlement and it is in risk condition

---

 In this research, prediction of concrete strength containing different aggregates using Non-destructive (Ultrasonic) testing through Artificial Neural Networks was carried out. For this purpose, aggregates with different properties were selected from the quarries, and then their destructive and nondestructive properties were obtained in laboratory. The significance of this research, using different aggregates with physical, mechanical and chemical properties also used two different test methods, such as Non-destructive static and dynamic testing, which are respectively uniaxial compressive strength and compressive wave velocity. Thus, this model includes various types of samples and the prediction model includes static and dynamic tests. The results showed that the use of artificial neural networks not only increases the accuracy, but also it reduces costs and time.

---

Buried pipeline system form a key part of global lifeline infrastructure and any significant disruption to the performance of these systems often lead to undesirable impact on regional business, economies or the living condition of citizens. In this paper the response of buried pipelines at fault crossings are studied. A fault movement can be resolved into an axial component, a lateral component in the horizontal plane, and a vertical component. Applying finite element method, the effect of various parameters such as anchored length, internal friction angle of surrounding soil, fault movement and fault crossing angle on the behavior of buried pipeline were studied. Nonlinear behavior for pipe and surrounding soil are assumed using beam-spring model. Results showed an increase in internal friction angle of surrounding soil increases strain and also normalized bending moment and axial force. Comparing bending moment at friction angle of 20^° and 40^° shows about 30% differences with certain crossing angle we can prevent producing large strain and bending moment on pipeline

---

Ground-Penetrating Radar (GPR) is a non-destructive and high-resolution geophysical method which uses high-frequency electromagnetic (EM) wave reflection off buried objects to detect them. In current research this method has been used to identify geometrical parameters of buried cylindrical targets such as tunnel structures. To achieve this aim, relationships between the geometrical parameters of cylindrical targets with the parameters of GPR hyperbolic response have been determined using two intelligent pattern recognition methods known as artificial neural network and template matching. To this goal GPR responses of synthetic cylindrical objects produced by 2D finite-difference method have been used as templates in the neural network and template matching algorithms. The structure of applied neural network has been designed based on extracting discriminant and unique features (eigenvalues and the norm of eigenvalues) from the GPR images and predicting all geometrical parameters of the targets, simultaneously. Also the template matching operation carried out using two diverse similarity approaches, spatial domain convolution and normalized cross correlation in 2D wave number domain. The results of the research show that both two employed intelligent methods can be applied for in situ, real-time, accurate and automatic interpretation of real GPR radargrams, however in general the neural network method has led to less error and better estimation than template matching to predict the geometrical parameters of the cylindrical tar

---

Landslide is one of the mass movement processes that occur in Iran and parts of the world every year. It causes huge human loss and economical damages. In order to check the stability of slopes in Kakasheraf basin, in the first step sliding areas were identified using the aerial photography and field surveys and then distribution map of landslide is provided. The impact of each of these factors which included dip, aspect, altitude, lithology, landuse and distance from the road and drainage are assessed through Arc GIS software merging the effective factors on landslide with the landslide distribution map. Then these factors were prioritized using AHP model. In this study, the fuzzy logic and density area method has been used in the Kakasheraf basin in order to identify landslide hazard zonation. The empirical probability index (EPI) has been used to assess and classify the models outputs in the landslide risk estimation.Results show that the fuzzy logic is more applicable method than density area model for mapping the landslide risk in Kakasheraf basin

---

Mountains are usually formation origin of their neighbor land surface features such as hillsides and plains. some problems and errors may occur in application of RS technique for generation of geology maps and in separation of these units from other similar units. The main objective of this study is to integrate RS and geomorphology approaches for identification of different geomorphology units and finally separation of debris lime stones from massive lime stones in Bahadoran region, Yazd province. For this purpose, a Landsat ETM⁺ image was acquired together with band ratios, principal component analysis and factor analysis approaches to generate lime stone distribution map. Results of this study show that (integration of RS and geomorphology sciences) can better generate the lime stone distribution map compared with the first one

---

 The average of shear wave velocity of the upper 30 m has so far been one of the reliable parameters in seismic site classification in different building codes, despite the numerous weaknesses in the exact explanation of site dynamic characteristics. In this study, an empirical relationship is obtained between the average of the shear wave velocity of the upper 30 m and the average of the shear wave velocity of shallower depths, based on 79 shear wave velocity profiles, in Mashhad. This is followed by the recommendation of proper depths for the dynamic analysis of the site effect based on the information of shear wave velocity profiles and resonance period distribution in the investigated area. The depth of the S-wave velocity profile investigation, required for the analysis of deposit effects has been estimated more than 30 m. whith exception of the southern and western parts of Mashhad (adjacent hillsides). Such depth is estimated as about 80 m for central, eastern, and north-eastern areas, where the resonance period is more than 0.7 s. Therefore, investigation depth of 30 m is only adequate for site classification based on the building codes, and for theoretical analysis deeper studies is needed, in Mashhad

---

Among dynamic parameters, the maximum shear modulus and damping ratio of the materials has the great importance in the dynamic analyses. In this paper, the high frequency GAP-SENSORs (20 kHz) is used to measure shear wave velocity of the sandy and gravelly Dacite materials in the cyclic triaxial equipment. The test procedure is as follows: firstly, a weak impact on the top cap triaxial specimen is induced. The induced low horizontal and vertical displacements time histories (strain in the range 10-6) are recorded using data logger with rate of 100,000 samples per secent for each channel. Then by traveling time and distance between two GAP-SENSORs, wave velocity is measured. Finally, by calculating the poison ratio via horizontal and vertical strains, the maximum shear modulus of the material is calculated. The results indicate that the value of wave velocity is increased by increasing confining pressure, and their values are dependent on particle size. Also, using the reduced vibration amplitude technique at different levels of specimen, damping ratio value is calculated by temping method. By comparing the results of shear modulus and damping ratio using local measurement of axial strain (with GAP-SENSORs in the cyclic triaxial) at small strain level with the results of the proposed method, cofirmed the high precision of the innovative used method.