Journal of Engineering Geology

---

Evaluation of ground response is one of the most important issues that should be considered in seismic geotechnical engineering field. Alongside the earthquake path associated to regional soil, generally earth movement in places with soft soil is greater than the movement in places with harder soil. This paper is focused to identify local soil condition of Ardekan city which influences on earthquake wave shaking. Therefore after drilling boreholes, implementing geotechnical investigations and down hole geophysical tests, the soil layer characteristics and thicknesses have been obtained. Then shear wave velocity along with soil density have been determined. With using these data it is developed a shaking geotechnical models for different city regions. Finally the ground movement parameters have been determined by   the available data obtained such as density, wave velocity along with using the equivalent linear method employing EERA program. This work was prepared for the return period of 75, 475 and2475 years. It is found that northwest region of city has the most amplification in comparison to other regions.

---

The choice of a suitable bearing capacity of soil becomes the most important issue to be considered in any project. This paper describes analytical investigation conducted to evaluate the ultimate bearing capacity of adjacent footings in various spacings of footings. Bearing capacity of adjacent footings is determined based on virtual retaining wall method by applying equilibrium between active and passive forces. Results indicate the ultimate bearing capacity of each foundation changes due to the interference effect of the failure surface in the soil and it depends on footings spacing. In the present study, effect of soil type, depth of adjacent footings and reinforced soil is investigated on bearing capacity of adjacent footings. Results indicate closely spaced footings, can decrease or increase bearing capacity of adjacent footings with respect to single footings. Also, reinforced soil increases bearing capacity of closely spaced footings with respect to single footings on unreinforced soil, it depends on footings spacing. Finally, the predicted results are compared with those reported from experiments, analytical and numerical results performed by others, indicating an acceptable agreement.
 

---

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.

---

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

---

Introduction
The effect of surface geology on seismic movement is known and acceptable and this effects can consider important factor in movement resulting from earthquake. studying intensity and dispersal of recent decade earthquake destruction indicated importance of construction effects and surface geology conditions more than ever. From view point of  engineering, earthquake importance is in light of effects that these earthquakes is created in construct such as dams, powerhouses, bridges, residential areas and industrial installations that in most cases, this constructs are building not only on rack mass surface but on earth surface, e.g . alluvial layers placed on bed stone. The effect of soil layers on earthquake waves is result of complex processes that this effect can exist as dynamic support under stability soil conditions that is called as intensification from it.
Material and methods
There are multiple methods in order to determining effect of construct and affecting in on earth potent movement features, that among them are covered less-cost numeral methods and more site output and beacuse of reason are using from these methods in order to analyzing respond to earthquake vibrations. In this paper try to using data resulting from drining 5 boreholes on Tabas city construction are studying the effect of and also comparing numeral methods of analyzing building such as equivalent linear and non-linear analysis for earthquake return periods of 75, 475, 2475 using NERA and EERA softwares.
Results
Taken together reinforcement rate and also maximum velocity in earth surface can explain that Dehshak region and Tabas center areas include more intensification conditions. On the other hand, south zone of Sarasyab sector and then Imamzade area include higher solidity and least intensification. Based on done studies are suggested to guided urban development programm more toward Hossein – ebne - Mousakazem Emamzadeh in order to exert from more suitable buildings. Also, regarding to EERA high-leval evalution and non-linear soil bahavior for earthquake with 2475 high return period is used from NERA software for analyzing construction effect to obtain maxium more realistic surface velocity ../files/site1/files/121/Rajabzadehi_Abstract(1).pdf

---

In urban areas, it is essential to protect the existing adjacent structures and underground facilities from the damage due to tunneling. In order to minimize the risk, a tunnel engineer needs to be able to make reliable prediction of ground deformations induced by tunneling. Numerous investigations have been conducted in recent years to predict the settlement associated with tunneling; the selection of appropriate method depends on the complexity of the problems. This research intends to develop a method based on Artificial Neural Network (ANN) for the prediction of tunnelling-induced surface settlement. Surface settlements above a tunnel due to tunnel construction are predicted with the help of input variables that have direct physical significance. The data used in running the network models have been obtained from line 2 of Mashhad subway tunnel project. In order to predict the tunnelling-induced surface settlement, a Multi-Layer Perceptron (MLP) analysis is used. A three-layer, feed-forward, back-propagation neural network, with a topology of 7-24-1 was found to be optimum. For optimum ANN architecture, the correlation factor and the minimum of Mean Squared Error are 0.963 and 2.41E-04, respectively. The results showed that an appropriately trained neural network could reliably predict tunnelling-induced surface settlement.

---

In recent years, with the growing use of the nailing method for stabilizing excavation walls, there has been a need for a comprehensive investigation of the behavior of this method. In the  previous studies, the behavior of nailed walls has been investigated in static and dynamic states and under different conditions. However, due to the different feature of near-field ground motions, it is  necessary to study the effect of these motions on the behavior of the nailed walls. Near-fault ground motion is significantly affected by the earthquake record direction and the rupture mechanism. So, in this study, to compare the effects of near-field and far-field ground motions, a two-dimensional (2D) soil- nailed wall was considered. PLAXIS 2D was used for the modeling of the soil-nailed wall system. An excavation with a dimension of 10 meters in height was taken into the account. In this study, 10 records (Five fault-normal near-field ground motion records and five far-field ground motion records), were recorded  on the rock and  applied to the model. These ground motion records were derived from the near-fault ground motion record set used by Baker. These records were scaled to the Peak Ground Acceleration (PGA) of 0.35g and then applied to the bottom of the finite element models. Mohr-Coulomb model was then used to describe the soil behavior, and Elasto-plastic model was employed for the nails. A damping ratio of 0.05 was considered at the fundamental periods of the soil layer. The results showed that the  generated values of bending moment, shear force and axial force in nails under the effect of the near-fault ground motions were  more than those in the far-ault ground motions. These values were  almost equal to 23% for the maximum bending moment, 30% for the  shear force,  and 22% for the axial force. The created displacement under the effect of near-fault ground motions was  more than that in the far-fault since a higher energy was  applied to the model in the near-field ground motions during a short time (pulse-like ground motions). In contrast, in the far-fault ground motions, due to the more uniform distribution of energy during the record, such pulse-like displacements were not observed in the system response. Increasing in nail length and soil densification, decreases the displacement of the soil-nailed wall but does not change the general behavior of the soil under the effect of near-field ground motions. Based on the obtained results, for a constant PGA, there were  positive correlations between the values of the  maximum displacement on the top of the wall and  the PGV values of near-fault ground motion records. However, the mentioned correlations were  not observed in the case of far-fault ground motions.

---

Introduction
Rock mass deformation modulus is one of  the major parameters has to be considered in the design phase of arch dams. Due to filling and discharging of reservoir and corresponding loading and unloading on the dam abutments, irreversible deformation takes place within the rock mass and consequently, increases the potential of creating a separation between dam body and abutments. Therefore, the rock mass modulus must be more than an alowable value in order to prevent arch dam failure. Regarding small core samples and lack of joints and other similar discontinuities in samples, the determined modulus through performing laboratory tests is higher than those obtained through in-situ tests. The available technique to estimate the rock mass deformation modulus is divided into two classes as direct and indirect methods. In direct methods, the rock mass deformation modulus is measured via performing in-situ tests such as plate loading test while it is estimated through empirical equations using rock mass classification and laboratory test results in indirect methods. These equations are developed based on regression analysis between the rock mass modulus calculated via in-situ tests, the rock mass classification and laboratory test results. Although application of these equations is simple and cost-effective, the results are doubtful and cannot be used in the design phase of arch dam due to the heterogeneous nature of rock mass and rock type variability. The numbers of micro-cracks which are developed after gallery excavation using drilling and blasting technique are more close to the loading plate. Thus, calculated modulus in these points is lower than reality. The displacement in the points far from loading plate was near to zero while the transmitted load which is calculated applying ASTM D4394 standard is more than reality in small galleries. Consequently, the calculated modulus was extremely larger than real values and sometimes even more than intact value. The empirical equations are site dependent and they are just applicable in sites with similar geotechnical condition. It is obvious that in-situ tests, such as plate loading, are the appropriate method in order to determine the modulus of deformation, however, due to some simplification in the data processing such as semi-infinite boundary condition, the application of numerical simulation as a data processing tool is more appropriate. In this research, the Beheshtabad dam was introduced and the geology characteristics of dam site were investigated. Applying direct and indirect methods, the rock mass modulus of dam abutments is calculated.
Material and Methods
The dam site is placed approximately at a distance of 2.7 km from the intersection of Koohrange and Beheshtabad river. In accordance with geological studies, the rocks in the site could be categorized in four units combined of Dolomite, Dolomitic Limestone, Limestone, Marl and Marly Limestone. Applying empirical equation the rock mass modulus of dam abutments is evaluated based on the laboratory test results and rock mass engineering classification systems. In addition, ASTM D4394 is applied to investigate the results of ten plate loading tests which are executed in the right and left abutments. To interpret the plate loading test results in the right abutment, a three-dimensional Fast Lagrange Analysis of Continuum (FLAC3D) model is developed.
Result and Discussion
To process the numerical simulation results, back analysis as a data processing tool is used. In this approach, the input parameters of numerical model will be changed in the way that the measured quantities by extensometers at the monitoring points are almost equal with the computed ones via numerical model at the corresponding points. Based on the sensitivity analysis carried out on the Mohr-Coulomb failure criterion parameters, the friction coefficient and cohesion variation do not affect the displacements calculated via numerical simulation as the more portion of gallery displacements are elastic. The error function is minimum when the rock mass modulus is 12 GPa and the horizontal to vertical stress ratio (K0) is equal to 0.5. The evaluated rock mass modulus based on the numerical simulation is two times lower than corresponding one evaluated applying empirical equation as a result of empirical equation uncertainty. Consideration of stress decrement under loading plate shows lower level of stress decrement under loading plate in ASTM D4394 compared to numerical simulation. This is why, the rock mass modulus, calculated based on ASTM D4394, increases dramatically by getting distance from the loading plate. 
Conclusion
The empirical methods estimating the modulus of deformation based on rock mass classification systems tend to evaluate large value of modulus especially for the weak massive rocks.
As a result of galleries dimensions and semi-infinite boundary condition assumed in ASTM D4394, the calculated rock mass modulus increases dramatically by getting distance from loading plate. Therefore, the numerical simulation was applied to process the plate loading test results. A new normalized error function was developed based on measured displacements and the rock mass modulus in the right abutment was determined 12 GPa which is very lower than the calculated value using ASTM D 4394. Also, as a result of numerical simulation, the rock mass is uniform. The stress increment perpendicular to the loading plate was calculated applying numerical simulation which is 0-90 percent lower than those suggested by ASTM D 4394. 
 

---

Existence of sodium elements in fine-grained of some soils causes dispersive phenomenon in them.  Failure to accurately identify dispersive clays leads to damage because dispersive clay soil particles disperse under certain conditions and wash away quickly. This research assesses dispersive degree of outcropped soils in southeast of Yazd. Finally, the modification of soil dispersive potential was investigated by using nanocellulose.  After performing a series of physical, chemical and mechanical tests, characterization/ specification of the studied soils were determined. Then dispersive degree was specified by conducting chemical, pin hole, crumb and double hydrometer tests. Finally, soil dispersivity stabilization was performed using sample preparation with 0.5, 1 and 1.5% nanocellulose. The results showed that the studied soils have moderate dispersive in borehole A and extreme to slightly dispersive in borehole B. Therefore, it can be concluded that the closer we get to the center of the plain, the greater the dispersibility. The results of the dispersibility stabilization soil tests indicate that the optimum moisture content and dry specific gravity increases and decreases. Uniaxial strength and CBR increases. Also, it shows that the increase in nanocellulose has a positive effect on the modified samples and improves the soil dispersibility in this area../files/site1/files/152/%D8%B2%D8%A7%D8%B1%D8%B9.pdf
 

---

In geotechnical engineering, rock mechanics and engineering geology, depending on the project design, uniaxial strength and static Youngchr('39')s modulus of rocks are of vital importance. The direct determination of the aforementioned parameters in the laboratory, however, requires intact and high-quality cores and preparation of their specimens have some limitations. Moreover, performing these tests is time-consuming and costly. Therefore, in this study, it was tried to precisely predict the desirable parameters using physical characteristics and ultrasonic tests. To do so, two methods, i.e. principal components regression and support vector regression, were employed. The parameters used in modelling included density, P- wave velocity, dynamic Poisson’s ratio and porosity. Accordingly, the experimental results conducted on 115 limestone rock samples, including uniaxial compressive and ultrasonic tests, were used and the desired parameters in the modelling were extracted using the laboratory results. By means of correlation coefficient (R²), normalized mean square error (NMSE) and Mean absolute error (MAE), the developed models were validated and their accuracy were evaluated. The obtained results showed that both methods could estimate the target parameters with high accuracy. In support vector regression, Particle Swarm Optimization method was used for determining optimal values of box constraint mode and epsilon mode, and the modelling was conducted using four kernel functions, including linear, quadratic, cubic and Gaussian. Here, the quadratic kernel function yielded the best result for UCS and cubic kernel function yielded the best result for E_s. In addition, comparing the results of the principal components regression and the support vector regression indicated that the latter outperformed the former.

---