Journal of Engineering Geology

---

Tunnel boring machines (TBM) are widely used in excavating urban tunnels. These kinds of machines have different types based on supporting faces and tunnel walls. One type of these machines, is the Earth Pressure Balance (EPB) type that was used in excavating the Line 1 Tunnel of Tabriz Metro. Different parameters such as geological conditions, rock mass properties, dip and machine specifications affect the efficiency of the machine. One method of predicting the efficiency of these machines is to estimate their penetration rates. In this study the value of TBM penetration rates are predicted by an artificial neural network. Predicting of this parameter is so effective for conducting in high risk regions by understanding the time of facing to these regions. The main result of this study is to forecast the penetration rate with an acceptable accuracy and to determine the effective parameters through sensitivity analysis measured by an artificial neural network.

---

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.

---

With respect to development of underground structures and their high construction costs in intra- and inter-municipal transportation, it is necessary to study the stability of such structures. In this research, tunnel stability of Mashhad Metro line 2 with 17 km length is studied. First, the type of sediments and geotechnical properties in Mashhad Plain are investigated. The SPT profiles were prepared using Rock Work 2006 software. The soil classification tests and XRD results show that the soils in this line are mostly clay such as Illite and Kaolinite types. Afterwards, because most of the soils in this line are characterized as fine grained, the ground settlement using PLAXIS V8 software was performed. According to the numerical modeling and the depth of tunnel, the optimum depth for tunnel was determined.

---

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

---

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.

---

Khorasan Razavi province is one of the areas with many chains of Qanat and Mashhad city in the center of this province has developed on areas with Qanat in the west direction. Loads caused by tall buildings and structures that built on old Qanats can make Qanats unstable and consequently Qanat collapse can lead to ground settlement. This paper deals with identifying the exact locations of Qanat chains by aerial photographs and evaluating the main factors that cause Qanat collapsing. Moreover the stability of Qanat was evaluated in numerical modeling by Plaxis software. The geotechnical data, Qanat depth, Qanat lining system and vertical load was used in modeling and after that, the extension of plastic zone around the Qanat underground tunnel was analyzed. The results of this research show that the rate of plastic zone extension and the influence of lining in Qanat stability decreases by increasing in the Qanat depth. As it is estimated while Qanat depth increases as much as one meter, it can tolerate more pressure as much as loads induced by a one-storey building.

---

The Minab (Esteghlal) dam site is located in east of Minab city in Hormozgan province. The Minab active fault cross the reservoir of dam and have an important role in leakage from the reservoir. The joint study of area in ten stations around the reservoir of dam display the four main joint sets. For assessment of leakage of reservoir, the permeability of rock masses outcrops in the reservoir is estimated by hydraulic conductivity HC experimental model. For this, the RQD, GSI and other characteristics of rock mass around the reservoir were measured in field studies. The results show that the permeability of embankments changes from 9.14×10-6 up to 2.02×10-5 m/s. Also the water lost for three different condition of minimum, mean and maximum level of water table is about 0.14, 0.20 and 0.29 m3/s. The results indicate that the discontinuities with trend of east-west and northeast-southwest and also shear fault zone of Minab have main effects in leakage of reservoir.

---

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.

---

Unpredicted and excessive tool wear is a major challenge in some tunnel projects. Due to more need for maintenance, abrasivity problems increese time and cost of excavation activity. A new laboratory test machine has been built in Engineering Geological Laboratory of Ferdowsi University of Mashhad, to obtain a proper view about relation between performance parameters of a TBM and the soil abrasion. In this paper the results of abrasion tests carried out on silica samples, as an abundant hard mineral, by the machin, are presented. Hence 36 tests have been performed on coarse silica sand samples, with various amount of deadweight and rotation speed at different times. The results prove direct relation between time, surcharge and rotation speed versus tool wear. For example a linear relation between tool wear and surcharge or rotation speed was observed. Besides a logarithmic relation was achieved for time effect

---

 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

---

Kerman city and its surrounding towns in terms of morphology, consist of a flat alluvial plain of fine silt and clay materials. These sediments have very gentle slope, and constitute the city's main infrastructure. Natural soils generally are structured by passing times due to the influence of environmental factors. Geological factors in Kerman alluviums have caused fine structure after deposition. This paper studies briefly sedimentary basin Kerman, mineralogical and geotechnical properties of the sediments of these areas.semi-qualitative analysis of samples has shown that the mineral deposits in Kerman are mainly illite, chlorite, smectite and calcite. then In order to evaluate the effect of depositional environments and geological history of the engineering properties of the sediments of the Kerman city, A large number of triaxial tests on reconstituted soil moisture greater than LL and different confining pressures, consolidated drained and consolidated undrained is done . In order to evaluate the applicability of soil behavior in normal conditions and the reconstructed curves and stress - strain sensitive soils and soil structure compared with standard sensitivity and then their resistance have been investigated. The results of triaxial tests can be used to interpret the depositional environments and geological history. Comparing the curves of stress - strain in natural and reconstituted samples indicates that in many cases the behavior of intact and reconstituted soils were similar and cementation and soil structure have not been much development. Also Comparison of electron microscopy images of reconstituted and intact samples, not random arrangement of particular structure and soil compaction within the city limits have confirmed. therefore, soils of Kerman are relatively similar by the influence of depositional environments and geological history. these soils have a lot of structure and cementation and are generally compact and strengthening.

---

In tunnelling in soil mass, in groundwater existing mode, liquefaction, elastic displacements and settlement in soils upon the tunnel, are the risks may attack the excavated underground space stability. In this case study that were performed on second line of Mashhad city subway route, information catched from Standard Penetration Test, in situ and laboratorial tests, were used to optimum numerical values search for soil engineering parameters that could optimize the TBM stationing level. In order to this goal attaining, intelligent, numerical and probabilistic methods were used and the reliability of intelligent and numerical methods with the Safety Factors of tunnel stability, investigated simultaneously. The results were denoting the accordance of intelligent models such as Genetic Algorithm (GA) and Multi objective Genetic Algorithm with Finite Element model's output. So these models could be complement of each others in planning and designing of tunnels and using of them advised in tunneling and excavations.

---

./files/site1/files/6.pdfExtended Abstract
(Paper pages 277-298)
Introduction
Preparation of uniform and repeatable reconstituted sand specimens of required density is a prerequisite for obtaining reliable results from experimental studies. Among different methods of reconstituted specimens, sand pluviation technique is widely adopted by researchers because of its unique advantage. In this study, a new curtain traveling rainer (CTR) is developed for large model sand bed preparation in experimental studies. CTR is a simple and low-cost system which is worked on the principle of air pluviation of sand. It provides specimens with wide range of relative density of sand bed (viz, 30%-90%) and very high degree of spatial uniformity of density distribution while reducing the time of preparation the specimens. A series of laboratory tests is carried out in order to study the performance of the proposed system and the effect of the curtain speed, curtain width, height of the fall and flow rate on the relative density and uniformity of sand specimens. For the sand used in the present study, it was observed that the relative density increases with an increase in the curtain speed and height of the fall. Furthermore, increasing the curtain width results in reducing the relative density.
Material and Methods
The calibration of geotechnical in situ tests in granular soil requires the preparation of large, uniform, and replicable specimens of a desired density. When preparing calibration chamber specimens, the adoption of techniques such as chemical impregnation and the freezing method in order to obtain undisturbed granular specimens becomes unfeasible due to the technical limitations and relatively high expense of these techniques. Conversely, the pluviation method has been widely used because of its ability to simulate the depositional mechanism of the soil and because of its applicability to a wide range of specimens, from small specimens for triaxial tests to large specimens for calibration chamber tests. 
soil sample was chosen from natural white-yellow silica sand mines of Firouzkooh, which can be categorized as poorly graded sand (SP) based on unified classification system (USCS). In order to control sand flow rate at the end of the pluviation path, a series of plates is designed to have rectangular openings with a width range of 2 to 4.5 mm. It should be noted that a sand reservoir is included in the transformer, which enables a uniform sand flow over the rectangular opening. the raining height is set at 100 mm to 500 mm with 100 mm steps. For more accuracy an extra test with  height fall equal 150 mm is performed.
A contiguous system of wheel-rail is operated for effective transmission of the traveling funnel over the entire sample surface in the circular container. The device’s jacking system and its related components are the cause of several limitations, which lead to the implementation of two-joint methods for keeping the sand rain height constant during the pluviation and sample preparation process.
Results and discussion
The CTR system comprises sand transfer compartments from the main hopper to the sample container and a rectangular opening at the bottom of the hopper, which controls the pluviation flow rate. The main concept driver of this research is to reproduce large samples in the most efficient time.  In order to recognize the uniformity of the reconstituted specimen in the vertical and horizontal directions, the variation of density is evaluated by placing 20 cylindrical molds within the specimen.
In this paper the effect of deposition intensity and the effect of height and flow rate on the sample relative density are evaluated.
Calibration of the sample preparatory device is very important in order to produce optional and repeatable samples with a specified relative density, in experimental studies and laboratory models. According to the test results, the effects of drop height and flow rate are investigated. Calibration graphs are presented in Figs. 1 and Fig. 2 for the proposed system in the case of 2.5 cm and 5 cm layer thickness.



Conclusion
This paper aims to extend the existing apparatus to achieve consistent low and high relative density sand samples. The preparation of low relative density samples is particularly important in liquefaction studies in geotechnical earthquake engineering. The comprehensive design and calibration of the CTR system can be concluded in the following points. The proposed method can easily be deployed to produce any arbitrary sample with a wide range of relative densities. Increasing the flow rate given a constant drop height leads to decrease in the relative density and is independent of layer thickness. Keeping constant the drop height and flow rate, higher relative densities can be achieved by increasing the curtain traveling velocity. There is a direct relationship between drop height and relative density. The  results  give  some  information  about  the  deposition  process  and  in particular about the terminal falling height. It can be henceforth stated that the performance of the proposed system is reliable and very acceptable due to high uniformity across the entire sample. 

 

---

---

Introduction
Determination of in situ stress-direction and magnitude are prerequisite for any oil well drilling and oil field development such as hydraulic fracturing. One of the simplest and most widely used methods is called borehole breakout analysis. Breakouts are compression fractures made in the direction of minimum horizontal in situ stress (Sh), if drilling mud pressure be lower than optimum mud pressure. Some borehole imaging logs such as FMI, FMS and UBI are appropriate tools for wellbore fracture detection. These fractures are distinguished in the logs as dark and symmetrical points (or lines) on both sides of the well and are used as an indicator for in situ stress studies. The size and shape of these fractures are strongly depend on the magnitude of the in situ stress. Therefore, many researchers suggested that by analyzing the geometric shape of the borehole breakout is an appropriate technique for estimation of in situ stress components. .... ./files/site1/files/0Extended_Abstract7.pdf
 

---

Introduction
Hazardous waste (solid, liquid or contained gases) is a waste with properties that make it potentially dangerous or harmful to human health or the environment. Site selection and suitable conditions for hazardous wastes landfill is considered as the final stage of waste management that they have high sensivity. The purpose of this study is to identify prone areas to hazardous waste landfill for Chaharmahal and Bakhtiari province using geographic information systems (GIS) as an important tool for the analysis of potential sites and the Analytical Hierarchy Process (AHP) and to provide solutions to optimize the positioning is executed. Firstly, criteria and limitations of environmental, economic, social and physical were determined, then layers of the criteria in GIS were prepared. In this study, the inappropriate areas were first removed from the model, and the suitability of remaining regions as a categorize criterion considered. Categorize criteria for paired comparison using AHP as an efficient tool for determining the relative weight parameters are used to measure and rank the expert choice application imposed. Then the implement paired comparison of the relative weights of the criteria and sub-criteria and criteria for each category were determined. After calculating the net weight and normal weight, normal weight based on standard maps in the GIS environment has been classifieds. Finally, by combining maps and applying criteria FA map, the final map was extracted.
Material and methods
The purpose of this research is to identify and prioritize appropriate areas of special waste disposal using multi-criteria decision-making methods. In order to locate using the GIS, first, identifying, evaluating and selecting criteria and constraints for the construction of landfill, in order to reduce the economic, environmental, and health costs. In the multi-criteria evaluation method, criteria are the basis of decision making, so that a set of criteria is combined and combined to achieve a single combination. In this paper, a two-stage process was used to locate the landfill site. In the first stage, which is recognized as the identification stage of prohibited areas according to different criteria, the study area is divided into two appropriate and inappropriate classes that will be eliminated as prohibited areas for the construction of landfills. In the second stage, the various factors are ranked and weighted according to the relative importance and, finally, places that receive the appropriate points are introduced as areas susceptible to the dumping of special wastes. In order to obtain the digital data of the criteria in the GIS environment from the digital elevation map (DEM), the specifications of the piezometric wells information are available from the regional water organization of the province. The available data such as geological map of the province at a scale of 1: 250,000, satellite images of Landsat and map of land suitability of the province, rainfall data of the synoptic stations of the province and the data of the Environmental Protection Agency were used. In general, the following steps have been taken in the process of locating:
- Identification of effective locating factors (limitations and factors)
- Digitizing and providing the required layers of information using the GIS package
- Identify and eliminate prohibited and inappropriate areas for landfill construction
- Classification and weighting of the factors and layers of information sought
- Integration of layers and the provision of a mapped rate and talent to determine the appropriate areas.
Results and discussion
1. Set limits
In this study, in order to select suitable sites for landfill particular, the criteria and limitations were determined. The information layers for each of the criteria were provided in the GIS environment.
2. Classification and weighting criteria
In the second stage, which is the stage of weighting and rating, of 14 effective criteria were used in site selection. AHP is one of the most efficient techniques, multi-criteria decision. This method is based on comparing factors and to study various scenarios to give managers and decision makers. This technique is one of the most comprehensive system designed for decision-making with multiple criteria.
3. Editor hierarchy to locate
Hierarchical structure is a graphical representation of a real complex problem, which mainly target the problem and at the next general criteria, sub-criteria and options are the way in AHP is used to calculate points based on comparison test.
4. Shipping Weight Matrix Binary comparison and decision-making
After compiling a hierarchical structure, the next step is to evaluate the elements by comparing the test. In general, if the number of options and criteria respectively m and n are then paired comparison matrix of options for comparison matrix m × m and n × n matrix will be a couple of criteria.
5. After weighing and preparing the normal weight of the options, the normalized weights in the GIS environment were added to the criteria map and the Raster and Weighted layers of each criterion were prepared. Due to the wide area of the studied area, the size of each pixel was 50 * 50 m. Then, using the Raster Module, the Criterion Map was combined and a zoning map was prepared for the special waste disposal site.
Conclusions
In the present study, according to various criteria influencing the Hierarchical Analysis Process for prioritizing the criteria in decision making, based on the results, the talent map of the area was prepared for special waste dumping, in which according to the final score of the layers, the area was classified into four appropriate, relatively suitable, relatively inappropriate and inappropriate classes. Suitable areas were 12.64%, relatively fairly 32.47%, relatively inappropriate 30.43%, and inappropriate zones 9.58% of the area of the talent map were included.
./files/site1/files/123/5Extended_Abstract.pdf

---

Introduction
The problematic collapsible soils are deposits with wind origin that constitute about 10% of the total area of ​​the earth. Several countries, including China, Russia, the United States, France, Germany, New Zealand, and Argentina have vast areas of collapsible soils. These deposits usually form a semi-stable honeycomb structure and are highly susceptible to sudden changes in the volume reduction due to becoming humid. Collapsibility and other related issues such as different subsidences, land cracks and landfalls seriously damage the infrastructures constructed on these soils.
 By the growing rate of urbanization in different parts of the world, the probability of construction on these soils and consequently water availability for these soils will increase; as a result, humidity increases and the collapse of these soils may occur. Therefore, studying the behavior of these types of soils is very important. Over the past six decades, many researchers have studied the collapse mechanism of collapsible soils due to becoming humid. Discussions on this subject are summarized in three categories: traditional methods, soil structure studies, and soil mechanics-based methods. In the present work, collapsibility and its controlling factors in the soils of Kerman city are investigated.
 
Material and methods
To determine engineering properties of Kerman deposits in this research, the geotechnical information was gathered and 50 core samples were extracted from different parts of the city. The sampling points were selected such that they could have a high overlap. X-ray diffraction (XRD) was applied to determine the mineral type and soil structures while scanning electron microscopy (SEM) was used to study grain arrangement.
Results and discussion
Geotechnical characteristics of the samples collected from Kerman plain deposits include their physical and mechanical properties. Based on the obtained results, this fine-grained sediment generally includes two CL and CL-ML groups. The mineralogy studies of Kerman city soils show that the minerals in these deposits are mainly illite, chlorite, illite-smectite, calcite, quartz, and gypsum. In order to study the collapsibility level of the soils in Kerman through the field studies, samples were taken from different parts of the city and the tests were carried out to determine the physical properties, collapsibility index, and structural studies. Through the SEM analyses, samples related to Haft Bagh area, Motahhari Town, and Pedar Town revealed an open structure and intergranular pores and thus a high level of collapsibility. On the other hand, in the majority of samples taken from the central part of the city, such as Esteghlal Street, Azadi Square, Bahmaniyar Street, and Hafez Street, the soil aggregates generally have corner-to-corner connectivity, with no specific order in their structure, and the arrangement of the particles is random and irregular. The orientation of the particles mostly shows no sharp pattern. In addition to soil particles, they have shown random and disorientated cavities with small sizes, suggesting the density and compactness of the soil indicating a low to moderate collapsibility. In some areas (e.g., Pedar Township and Motahhari Township), crystalline salt and gypsum crystals are clearly seen. It is expected that by increasing the amount of water, these salts dissolve and their effects can be observed as dissolution cavities.
 The dissolution of soluble crystals can also reduce the strength of the soil structure and ultimately lead to soil degradation. Calcite crystals are also found in some places in the form of calcite cement among the grains, sometimes as single crystals, and sometimes as lime nodules within the soils of Kerman city. Among the stated criteria in this research, Denisov, Holtz, and Hill criteria, the Russian regulations and ASTM standards were employed to assess the potential of the studied soil collapsing. Based on the criterion of the construction regulations of Russia, it was found that the deposits of the city of Kerman are mainly collapsible (L>-0.1).                     
Moreover, based on the Denisov criterion (if e/eL>1.5 the soil is non-collapsible, if it is between 0.75 and 1.5, the soil is prone to collapsing, and if it is between 0.5 and 0.75, the soil is severely collapsible), soils of Kerman are within the range of collapse-prone soils. Finally, based on the ASTM criterion, in some areas of the city like Motahhari Town, Pedar Town, and Haft bagh, soils show a high collapsibility. In comparison, in the central parts of the city, the values of this criterion vary between 0.15 and 11, suggesting the presence of soils with a moderate collapsibility. Comparing the results obtained using these criteria it is seen that areas with a collapsible behavior are relatively similar collapsibility results are obtained.
Conclusion
Based on the achieved results, fine-grained sediments of Kerman city are mainly composed of CL and CL-ML groups. Mineralogy results indicate that the minerals in these deposits are mainly illite, chlorite, illite-smectite, calcite, quartz, and gypsum. SEM results for the central part of Kerman city confirm the compressed and densely compact form of soil particles. The results obtained, using the construction regulations of Russia show that the soils in the study area are collapsible. According to the Denisov criterion, they were found to be prone to collapse. Finally, based on the ASTM results for the central parts of the city, soils exhibit a low to moderate collapsibility. However, in some areas of the city, such as Motahhari and Haft bagh, soils show a complete collapsibility behavior../files/site1/files/131/1Extended_Abstract.pdf
 

---

Introduction
Mashhad city, the second largest metropolis of Iran, is located in an arid and semi-arid region. Overexploitation of groundwater in Mashhad plain has caused up to 22.5-meter drop in the groundwater level from 1984 to 2013. The groundwater depletion in the unconsolidated aquifer has resulted in subsidence and cracks on the land surface. To determine the land subsidence rate map and the reasons for hot spot subsidence, the latest Envisat images of the ESA Space Agency's Archive for Mashhad plain were used. leveling and GPS data were combined with the radar interferometry results and the annual subsidence rate maps with high precision were obtained. Finally, the geology and soil texture maps of study area are compared to the land subsidence map.
Methods and results
To assess the land subsidence in Mashhad plain three methods of leveling, GPS and Insar are used. Leveling data are available in three profile of of Mashhad-Quchan (BCBD), Mashhad-Kalat (BDBE) and Mashhad-Sarakhs (BEBN) in two time interval of 1994-2003. The highest rates of subsidence in the BCBD, BDBE and BEBN lines are 7, 3.5 and 8.1 cm/year, respectively. Six permanent GPS stations have been installed in Mashhad plain, among them, NFRD, GOLM and TOUS have recorded the land subsidence, with the highest annual rate of 21.2 cm/year at TOUS Station. The third method applied to assess the history of land subsidence was InSAR radar interferometry which provided the extent and pattern of subsidence in all of the study area. For this, 23 images of the Envisat ASAR are processed during the 05/24/2010 to 06/30/2003 time period. The highest subsidence rate estimated by this method was 32 cm/year in the northwest of Mashhad. In general, two subsidence bowls, in the northwest and south east of Mashhad city are identified. Fig. 1 shows the annual subsidence rate map in Mashhad plain. Using the root-mean-square error (RMSE), the accuracy of the InSAR method was verified with GPS and leveling data.
Discussion
The rate and distribution of land subsidence in Mashhad plain are affected by geological factors such as soil texture, deposit thickness, geological structures and groundwater drawdown. The geological and geophysical studies and exploratory drilling results in the Mashhad Plain indicate that the bedrock morphology is very rough. The bedrock outcrops in some places while in some other places covered by more than 300 meters alluvial deposits. Generally, by distance from the mountain, alluvium thickness and as a result the likelihood of subsidence would be increased. Mashhad plain is surrounded by the active and quaternary faults in the north and south edges. In the north of Mashhad plain Marly bedrock is uplifted by Tous fault and outcropped in the north of fault. In the south of Mashhad two normal faults have resulted to the increase of alluvium thickness in south and central of Mashhad plain. The change of river pathway also let to deposition of a sequence of the fine-grained and coarse-grained soils in central of plain between Toos and southern branch of South Mashhad fault (F2).
used to draw the cross section
In order to evaluate the subsurface conditions and its effect on the land subsidence, the soil texture are studied using the deep water wells and piezometers log (Figure 2). Fig. 3 shows the longitudinal section (northwest to southeast) of the area. As it can be observed, the soil texture includes of alternation of fine and coarse grains layers (Figs. 4). In this condition, sandy soils help to shortening the drain path of clayey layers and leads to acceleration of the consolidation. The average rate of annual subsidence in the area is 14 cm for one meter of drop in the groundwater level.
Nowadays, in the urban area, due to the urban sewage waters, there is a rising of groundwater level.  Therefore, no land subsidence has occurred in the central parts of the city. It is expected by completion of urban sewage network about 62 million cubic meters of sewage water will be eliminated from the aquifer recharge, which will cause a notable drop in the groundwater level and prominent land subsidence in specific area of the city. Considering the geological conditions and the operation of the existing faults, it is expected that in the case of groundwater drop, no significant subsidence will occur in south of the F2 fault, due to the decrease in the alluvium thickness and to the coarse texture of the soil. But in the northern and northeastern parts of the city, which are located between F2 and the Tous faults, high rate of land subsidence is expected.
Figure 4: The cross section of soil texture and the annual average rate of land subsidence and groundwater level drop
Conclusions
Using the radar interferometry processing, the highest annual rate of subsidence in Mashhad plain is about 32 cm/year. Land subsidence in Mashhad plain has an increasing trend and the geological conditions have a critical role in the subsidence rate and its pattern. Generally, soil texture near the mountain area in South is coarse and grain size decreases toward the center of the plain. But because the outcrop of Marly formation in the north slopes, soil texture is mainly fine grains. In the center of Mashhad plain soil texture constituted of fine and coarse grains which are converted together as inter fingering facieses, which have a critical role in decreasing of the consolidation time and increasing the land subsidence rate. It is predicted by complimenting of the urban wastewater network, the groundwater level will be dropped in the city area and the northwest and southeast subsidence ellipsoids which already can be seen will be connected together. Therefore, the area between F2 and Toos faults, will be shown the highest rate of subsidence, due to high thickness and fine-grained soil texture../files/site1/files/133/5Extended_Abstracts.pdf

---

Introduction
CO₂ injection in deep geological formations, such as depleted oil and gas reservoirs, in addition to the environmental benefits, is one of the effective method for enhanced oil recovery (EOR) as tertiary EOR. Presence of reservoirs with a pressure drop which require injection of gas in the southwest of Iran and having the technical and environmental effects of CO₂ injection have created a huge potential for CO₂ injection to EOR in this region. In the first step, to perform CO₂-EOR, the geomechanical assessment is needed to find out pore pressure, in-situ stress magnitudes and orientations and fractures and faults conditions. In this paper, the initial in-situ pore pressure is predicted using modified Eaton method for 47 wells in the length of the study field and calibrated using repeat formation test and mud pressure data. In-situ stress was obtained by the poroelastic method for 47 wells in the length of the study field and calibrated using leak off test and extended leak off test. Then, the orientation of in-situ stresses is obtained based on image logs. Hydraulical and mechanical activities of fractures and faults were performed by critically-stressed-fault hypothesis
Material and Methods
In this paper, the initial pore pressure is calculated using modified Eaton method and other corrections that are proposed by Azadpour et al. (2015). The estimated initial pore pressure is validated using mud weight pressure (Pmw) and repeat formation tester (RFT) data. In-situ stresses are composed of three orthogonal principal stresses, vertical stress (S_V), maximum horizontal stress (S_H), and minimum horizontal stress (S_h) with specific magnitude and orientations. The magnitude of S_V is calculated by integration of rock densities from the surface to the depth of interest. The poroelastic horizontal strain model is used to determine the magnitudes of the S_H and S_h. Then, the estimated minimum horizontal stress from poroelastic horizontal strain model is validated against direct measurements of LOT and XLOT tests. The orientation of breakouts was determined based on compressively stressed zones observed in the UBI log and using Caliper and Bit Size (BS) logs. The hole elongates perpendicular to the S_H and breakouts develop at the azimuth of S_h. Fractures and faults reactivation analyses are very important, they can potentially propagate upwards into the lower caprock and further through the upper caprock due to CO₂ injection. Fractures and faults identification were performed based on image logs. Based on performed seismic interpretations by NISOC (National Iranian South Oil Company), 15 faults have been detected in the field. Fractures and faults conductivity and activity in the current stress filed affect on fluid flow and mechanical stability or instability of the CO₂ injection site. Critically stressed fault hypothesis, introduced by Barton et al. (1995), states that in a formation with fractures and faults at different angles to the current stress field, the conductivity of fluids through their apertures are controlled by the interplay of principal stress orientations and fracture or fault directions. Hence, conductive and critically stressed fractures and faults in the current stress field were evaluated using critically stressed fault hypothesis. Fractures and faults are plotted in normalized 3D Mohr diagrams (normalized by the vertical stress), therefore conductive and critically stressed fractures and faults were determined.
Results and discussions
The maximum distribution of initial pore pressure was 20-25 MPa in the field and the average of initial pore pressure was 25 MPa in the field. Unlike the World Stress Map, the stress regime is normal in the reservoir. Because the Kazeroon fault and Dezful Embayment act as a strike-slip tensional basin, resulting in the subsidence of Dezful compared with other regions. The frequency distribution of calculated in-situ stress in 47 studied wells in the length of the field has been presented. The maximum frequency distribution of S_V, S_H and S_h were between 60-70, 50-60 and 30-40 MPa, respectively. A large amount of fracturing is observed in 20-25 m below the caprock. Based on the continuity of their low amplitude traces on the acoustic amplitude image of UBI, fractures are classified into 4 classes: discontinuous-open, continuous-open, possible-open and closed fractures. OBMI and UBI image logs processing were performed in 7 wells. As can be seen from the image log, and caliper analysis the most dominant strike of S_H around the well is 27^◦ and S_h strike is 117◦. These have two dominant orientation, some faults are along the strike of the Zagros fold-thrust belt (NW-SE) and the others are perpendicular to the Zagros fold-thrust belt strike (NE-SW).
Based on the normalized 3D Mohr diagrams it is clear that the fractures and faults that are oriented to the S_H will be the most permeable, because the faults and fractures experience the least amount of stresses in the direction of S_H and they have minimum resistance to flow in this direction, therefore will have relatively high permeability. Also, results showed the faults number 15, 6, 10 and 2 will be the most dangerous faults during CO₂ injection.
 
 

---

Introduction
 Analysis of time, location and magnitude of foreshocks and aftershocks has been one of the most important cases for experts in various scientific fields such as: seismology, structural engineering and crisis management, and other interrelated fields. Since this analysis and the result of studies on seismotectonic and cases of earthquakes help us identify the foreshocks and aftershocks with the goal of decreasing losses and nervious stress of the injured community in quake-stricken areas and skilled crisis management. The cause fault of earthquake plays the important role in foreshocks and aftershocks of the earthquake. So, study on fault behaviour is a suitable method for analyzing and routing the basic parameters of foreshocks and aftershocks. Also, foreshocks and aftershock are important parts of any earthquake in a seismic area. The analysis of the basic parameters of the foreshocks is one of the most practical researches for reducing the risk of earthquakes. The identification of behavioral pattern of foreshocks can help researchers detect the active fault conditions for the occurrence of earthquakes in different areas. The present study is concerned with the study of behavioral patterns earthquakes, foreshocks, and aftershock of Zohan earthquake. Experience of large aftershocks in different parts of the world indicates that, following earthquakes and depending on seismic-tectonic conditions, large aftershocks are likely to occur in the earthquake-effected zone, which will aggravate the damage caused by earthquakes (Omi et al., 2013). The main factor contributing to the worsening of damage caused by aftershocks is the performance of structures that are weakened but not destroyed by main earthquakes and are, thus, highly likely to be destroyed by large aftershocks (Saket and Fatemi Aghda, 2006).
Material and methods
The present paper makes use of data collected in a real earthquake and similar expriences in other earthquakes for presenting a practical pattern for predicting primary earthquake patterns, determining the location, magnitude, and time of aftershocks. The target of this case is decreasing the effects of earthquake. To this end, we used the results from studies on basic parameters of foreshocks and aftershocks of Zohan earthquake, and 2012 earthquake in South Khorasan province. The rationale for selecting the aforementioned studies is: location of event, the Zohan earthquake, had been identificated as an area with high risk for the occurrence of earthquakes, although there has been no wide-scale earthquake in this area in the last two decades. These conditions are important causes for more concentrated studies on this area because there is a high chance for wide-scale earthquakes striking this area.
Result and Discussion
In this part of research, we conduct a study on the location, magnitude and depth of foeshocks. Some of the world-wide research suggested that these data can help to predict the time of  mainshocks. Studies conducted on the variations of frequency in foreshocks can follow this goal.
In this paper, the available statistical data such as periodical variations of seismicity in the weeks leading up to the main shock can be used as a tool for estimating the approximate time of a future important earthquake. The weekly variations of seismicity before Zohan earthquake indicate a relative increase and then decrease within a 100 km radius around the epicenter of the main shock.
 
Table 1: Variations of frequency of foreshock based magnitude before Zohan earthquake

Week before main shock	Frequency of foreshock in the Radius of 100Km from main shock	Frequency of foreshock(with M>2.5) in the Radius of 100Km from main shock
6	0	0
5	1	1
4	1	0
3	2	0
2	5	3
1	2	0

 
Studies on numerous earthquakes in Iran and other regions in the world show that the distribution of aftershocks can be related to fault type or the direction of principal stress (Saket and Fatemi Aghda, 2006) and (King et al., 1994). Whereas maximum Coulomb stress change is related to maximum principal stress in earthquakes, the concentration of aftershocks can coincide with the direction of maximum principal stress (σ₁) of the causative fault in mainshock. Considering the direction of maximum principal stress and its adaptation to the scattering of aftershocks, the above hypothesis is confirmed.
Also studies on frequency changes and seismic quiescence of small aftershocks help us in predicting future aftershocks. The results the of presented research by Itawa (2008) on the World earthquake catalogue suggest that seismic quiescence theory is true for different regions of the world. Based on the results of the  study mentioned above, this case can be used as a tool for predicting large aftershocks in Zohan earthquake.

Fig 1. Adaptation of direction of maximum principal stress with scatering of the aftershocks of Zohan earthquake. a: direction of maximum principal stress (σ₁) of the causative fault in mainshock. b- scaterring of the aftershocks
Table 2: Seismic sequience versus magnitude of aftershocks

Row	Seismic Quiescence for aftershocks	Aftershock Magnitude	Data and Time of aftershocks
1	13	3.0	2012/12/05 17:21:03
2	36	3.4	2012/12/05 17:57:03
3	161	3.1	2012/12/05 20:38:09
4	3906	3.9	2012/12/08 13:44:19

In addation, frequency of aftershocks and certain time distance (seismic quiescence) between their can use precursors for detecting the time of large aftershocks. The relevant analysis in this study showed that methods such as: time series beside seismic quiescence can help in conducting a more accurate time forecast of large aftershocks.
Conclusion

• The results of this research suggest that we can identify some of the charactristics of the main shock by focusing on location, magnitude and depth of foeshocks.
• In Zohan earthquake, the direction of maximum principal stress is adpated to the scattering of aftershocks, and this case suggests that there is a specific relationship between them.
• The relevant analysis in this study showed that the methods such as: time series beside seismic quiescence can help conduct a more accurate time forecast of large aftershocks../files/site1/files/144/saket.pdf