Journal of Engineering Geology

---

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

---

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

---

The 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 Schmidt hammer provides a quick and inexpensive measure of surface hardness that is widely used for estimating the mechanical properties of rock material such as uniaxial compressive strength and Young’s modulus. On the other hand, Schmidt rebound hardness can be used for a variety of specific applications. In the mining industries, it is used to determine the quality of rock, which is common practice when constructing rock structures such as those found in long wall mining, room and pillar mining, open-pit mining, gate roadways, tunnels, dams, etc. However, a number of issues such as specimen dimensions, water content, hammer type, surface roughness, weathering, testing, data reduction and analysis procedures continue to influence the consistency and reliability of the Schmidt hammer test results. This paper presents: a) a critical review of these basic issues and b) avaluate the effect of temperature, moisture and uniaxial compressive stress on Schmidt hammer hardness. It was found that water content has a significant effect on the Schmidt rebound hardness (SRH) of rocks. So that increase of water content substantially reduced the SRH of samples. Temperature also had a considerable influence on the SRH. However, relationship between SRH decreases with increasing temperature for tested samples were linear. Also tests results showed that uniaxial loading of samples increases the SRH values.

---

In the narrow geosynthetic-reinforced retaining walls a stable rear wall exists in a short distance and so there is no enough space to extend appropriate length of reinforcements. In this case, the probability of overturning of retaining wall increases especially when subjected to earthquake loading. To increase the stability of the wall, reinforcements may be connected to the stable rear surface. Alternative solution is the utilization of full-height cast in-place concrete facing in order to resist the earth pressure by combined actions of reinforcements pullout capacity and facing flexural rigidity. One of the main questions about this type of walls is the portion of earth pressure resisted by the facing. In this study, the seismic earth pressure of narrow geosynthetic-reinforced backfill on rigid facing was evaluated using limit equilibrium approach and horizontal slices method. The critical failure surface was assumed to extend linearly from the wall toe to the rear surface and then moves along the interface of the backfill and rear surface up to the backfill surface. The effects of various parameters such as wall aspect ratio have been investigated. The obtained results show that the applied soil pressure on wall facing will be increased with depth in the upper part of the wall according to the Mononobe-Okabe equation, but its pattern is inversed in the lower part of the wall and it decreases until it reaches to zero at the wall toe. The results of analyses indicate that the attracted soil thrust by the facing increases with lessening of backfill width.

---

Introduction
The prediction of landslide occurrence in a region is very important in reducing the risks and damages caused by this.landslide as a natural disaster in Iran caused a lot of life and financial losses to Iran annually. According to the National Committee on Natural Disaster Reduction of the Ministry of the Interior in 1994, the share of annual damage caused by mass movements in Iran is estimated at 500 billion rials. In the meantime Kurdistan province is the third province in terms of landslide phenomenon after Mazandaran and Golestan. If considering the area is at a higher level. The city of Bijar in this province has a high potential for a wide range of landslides with a combination of mainly mountain topographical factors, lithologic conditions and positioning between two major faults in the region. In this research, using quantitative methods and models on the quantitative  factors of this phenomenon based on the level of information given by past mass movements and influential factors, focusing on artificial neural network method, susceptibility zones were determined by determining the possible risk level.
Knowing such natural events requires proper management of the risks posed by them. On the other hand, artificial neural network as a quantitative model is capable of learning, generalization and decision making, and less need to analyze the accuracy of data in comparison to statistical methods. Map of the susceptibility of the areas to the landslide is an important tool for landuse planning. However, there are many issues in the formation of this phenomenon, which, due to the complexity of the natural processes arising from the relationship between the outcome (dependent variable) and the factors (independent variables), puts into question the general zoning of such areas.
Methodology
Bijar is located in the northeastern part of Kurdistan province, along the longitude of   47 ' 29° to 47 ° 47' east, in latitude 35 ° 35 'to 35' 59 °north. In recent years, the development of the Geographic Information System (GIS) and spatial analysis techniques have improved the risk of indirect zoning. In this regard, artificial neural networks can cover a significant part of these needs.Implementing the neural network model requires learning data. Without learning data, it's virtually impossible to make neural networks. In this paper, learning data shows the occurrence of landslides which have geographical coordinates and were obtained from the Kurdistan Province Natural Resources Organization. In general, learning data in GIS and remote sensing can include data or raster, which in this paper is a point phenomenon and has 144 cases.  However, because of the large extent of the study area and the low number of them, as well as the lack of risk of any landslide zone (from low to very high), the points should be classified as well, and, in terms of numbers, Acceptance. Also, the number of points of relative value In terms of numbers, the conditions are the Normal and the same (that is, the appropriate geographical distribution and distribution in each class) would be more accurate; thus, to create a classifiable spectrum of the AHP Was used. It should be noted that all the maps were standardized in the format and format of the Raster in a matrix (698 rows in 897 columns) identical with a size of 30 * 30 meters. This means that each map has 626,106 pixels of varying value and somewhat similar. In addition, the AHP model was used to categorize the studied area from very desirable (hazardous) to very undesirable (very dangerous) areas. Also, 33 points were added to the learning data on different levels of the map derived from the AHP model. But in order to verify accurately the model, only landslide occurrences were considered.
In order to find out the factors of landslide in Bijar, a map of slope, Aspect, elevation, distance from the fault, distance from the road, distance from the river, Drainage density, lithology and land use using ArcGIS software were prepared and digitized.
After compiling and categorizing these variables, at first, each of the effective criteria in the field was divided into six sub-criteria (land suitability for landslide) from very desirable to very undesirable conditions. The present study utilizes the technique of multi-layer propspert neural networks using post-propagation algorithm (BP). In addition to correcting and editing the layers, the neural network model was implemented using the classification method and applying two types of functions (linear and sigmoid). Then, using the test-error method, the study of the magnitude of the error and the period of the repetition and the change in the number of hidden layers and weights, both functions were performed. Finally, the sigmoid function, which yielded a better result, was selected as the proposed and final function.Order to verify the (accuracy) of the map taken with the existing landslide zones, the final map of the neural network model was again transferred to the ArcGIS software. Finally, the available landscapes on the map resulted from the adaptive neural network model, which, by comparison, gave a percentage and amount Accuracy of each class was achieved.
Result
The input layer were calculated to six classes based on the desirability of mass movements. This decision approach reduces the complexity of the network and improves its performance.
For this purpose. The AHP method was used to define non-slip pixels and range classification.
To implement this method, 9 variables discussed, were scaled up to the most suitable and un suitable range. The final weight of these variables was obtained by using the Thomas saati pair comparison (Table 4), the study area was divided into five categories according to the map for land suitability for landslide hazard. From each class, the 20-pixel from AHP model was selected for network learning in a completely randomized manner. The proposed model is an artificial neural network of MLP multi-layered perceptron with levenberg-marquardt learning algorithm. An early stopping method was used to improve network optimization. Several hidden layers were tested to find the best results. It should be noted that in the structure of all networks, at least the optimal design with the middle one is used, but in their structural composition they are also used with mid-duplex networks. In this paper, the use of tow mid-layers showed better results.  In all Simulations have been made, the mean square error index, as a guide, indicates the network performance in learning the existing model. By changing the number of intermediate neurons and changing the weights as try and error, the most appropriate network model was obtained for the purpose. In this study, the structure of the network with 9 input layers, 2 hidden layers, 1500 repetitions in both functions was accepted as the final structure. The main structure of the neural network with two linear and sigmoid functions was prepared with acceptable error, and the study area was analyzed with a total area of ​​564 km2 with 9 input variables converted into raster data to 30 × 30 pixels. From 564 km2 based on the sigmoid function 61.17% and based on the linear function, 72.76% of the area is unsuitable and very unsuitable in the area where expose to high risk. In both networks, there were very few areas in both optimal and moderate classes (Figures 16 and 17), which indicate the high talent of the area for landslide as a threat. Then, ArcGIS software was used to evaluate the efficiency and accuracy of the model. For this purpose, the point of landslide and zoning maps were combined, compared and anlayzed. The results showed in the sigmoid function 75 items of Landslides were in a very unsuitable range, which included 61% of the total of region.
Conclusion
 In the linear function, approximately 69% of the landslides are in a very unsuitable range and the unsuitable results are about 57%, which results in the success of the model designed in the neural networks (MLP). In the end, the network with sigmoid function is negligibly better than the linear function network.The results show that Bijar and its functions are relatively prone to occurrence of landslides, so that nearly 60% of the city's area is a high risk area with a high risk and only 2% is a low-risk region. The hazardous areas are mainly located around the city of Bijar especially southern and southeast. These areas correspond to high altitudes and maximum fault density and lime lithology with marl (Qom Formation). The model can be very challenging, because of innovative nature of the research, that means need more detailed and comprehensive studies../files/site1/files/121/neiri_Abstract.pdf

---

Introduction
When saturated sandy soils are subjected to seismic loadings, the pore water pressure gradually increases until liquefaction happens and settlement occurs during and after an earthquake. The mentioned problem is attributed to rearrangement of grains and redistribution of voids within the soils. Over the years many methods have been presented to increase liquefaction resistance. However, the main methods utilized in liquefaction mitigation are classified as densification, solidification, drainage and reinforcement techniques. Utilizing scrap tires in soils is a kind of soil reinforcement which has a wide range of application.
Waste material expulsion is one of the environmental problems each country faces. Accumulation of non-degradable polymeric materials in landfills has serious environmental consequences. Efforts to find new ways of soil reinforcement have drawn the attention of researchers towards the use of new recycled materials like scrap tires derivatives. Derivatives of scrap tires have different applications in civil engineering such as reinforcing soft soil, as a drainage layer in landfills and as filler materials.
Material and methods
In this paper a series of 1g shaking table tests were performed to investigate on the effect of tire powders-sand mixture in reducing liquefaction potential, settlements after earthquake and pore water generation. Shaking table is made of Plexiglas with inner dimensions of 200×50×70 cm. At bottom of the container a void chamber is made by using a number 200 sieve so that the saturation process could be done gradually and uniformly. A plastic plate was rigidly fixed at the center of container to separate reinforced and unreinforced samples from each other and waterproofing carefully. Therefore two models (reinforced and unreinforced) can be tested at once with the same input acceleration. An absorbing layer of foam with 2 cm thickness was employed to decrease the effect of boundary conditions in order to avoid a direct confrontation model with a rigid container. Firoozkuh No. 161 sand and tire powders were used for the mixture in reinforced side, and pure sand in unreinforced side. In this study 4 mixture ratio (TC=5%, 10%, 15% and 20%) were done. Both of unreinforced (pure sand) and reinforced (tire powders-sand mixture) models were prepared by wet tamping method, in which soil is mixed with 5% water. Each model was prepared in six layers. The required weight for each layer was considered based on the desired density (relative density is zero) and exact volume of the layer. Each portion was placed into the model container and then tamped to reach desired level. Carbon dioxide (CO₂) was allowed to pass through the specimen at a low pressure in order to replace the air that trapped in the pores of the specimen. Then water was allowed to flow upward through the bottom of the container at low pressures in order to flush out the CO₂ that cause increasing the final degree of saturation. Vibration with approximate uniform amplitude and 2 Hz frequency was applied to the container.
Results and discussion
Results indicate that acceleration within the soil tends to be increased towards the soil surface. On the other hand, after initial liquefaction (that occurred at un-reinforced models), acceleration is decreased due to the increase in excess pore water pressure. Also, it can be seen that the increase in tire powders ratio remarkably reduces the maximum excess pore-water pressure ratio. The settlement of the tire powders-reinforced models is significantly less than the unreinforced models, and with the increase of the tire powder percentage shows a very small increase of volume. The outcomes show that the value of the mean damping ratio versus the shear strain range of 0.01 is increased with the increase in tire powder content. Shear modulus is obtained from the ratio of the difference in maximum and minimum stress and strain developed in desired loop. The maximum of the shear modulus in reinforced models is more than the unreinforced models.
Conclusion
The main aim of the present paper was to investigate the influence of reinforcing a saturated sandy soil with tire powders on the soil dynamic properties and the mitigation of liquefaction potential. The following conclusions were drawn from this research.
- The increase of pore-water pressure leads to a reduction in soil shear stiffness and acceleration amplitude.
- Reinforcing sand with tire powders reduces the excess pore-water pressure ratio because of liquefaction and increases liquefaction resistance. 
- Reinforcing sand with tire powders decreases settlement caused by liquefaction.
- The damping ratio decreases at large shear strain as liquefaction occurs.
- Maximum shear modulus and mean damping ratio of reinforced soil has been increased with increasing tire powder content in the mixture../files/site1/files/123/3BahadoriFarzali.pdf

---

Maximum surface settlement (MSS) is an important parameter for the design and operation of earth pressure balance (EPB) shields that should determine before operate tunneling. Artificial intelligence (AI) methods are accepted as a technology that offers an alternative way to tackle highly complex problems that can’t be modeled in mathematics. They can learn from examples and they are able to handle incomplete data and noisy. The adaptive network–based fuzzy inference system (ANFIS) and hybrid artificial neural network (ANN) with biogeography-based optimization algorithm (ANN-BBO) are kinds of AI systems that were used in this study to build a prediction model for the MSS caused by EPB shield tunneling. Two ANFIS models were implemented, ANFIS-subtractive clustering method (ANFIS-SCM) and ANFIS-fuzzy c–means clustering method (ANFIS-FCM). The estimation abilities offered using three models were presented by using field data of achieved from Bangkok Subway Project in Thailand. In these models, depth, distance from shaft, ground water level from tunnel invert, average face pressure, average penetrate rate, pitching angle, tail void grouting pressure and percent tail void grout filling were utilized as the input parameters, while the MSS was the output parameter. To compare the performance of models for MSS prediction, the coefficient of correlation (R²) and mean square error (MSE) of the models were calculated, indicating the good performance of the ANFIS-SCM model.

---

Introduction
Past research studies have demonstrated that seismic ground motion can vary significantly over distances comparable to the dimensions of long span engineering structures. The accurate determination of earthquake ground motion at the base of long span structures such as dams and bridges whose piers are located on the valleys surface is one of the most important issues in earthquake engineering. In this paper, the spatially variable earthquake ground motions are generated at stations located on the valley slopes, considering the topography effect of a triangular valley. To this end, the simplified geometry of the valley of Masjed Soleyman embankment dam has been used for numerical modeling. The spatially varying ground motions are simulated by using spectral representation method. According to this methodology, the generated time histories are compatible with prescribed response spectra reflecting the wave passage and loss of coherence effects. This method assumes that the response spectrum is identical for all stations i.e., they have the same amplitudes and frequency content. This assumption is not valid for stations located on valley surface in which the amplitude and frequency content of the seismic waves are changed considerably by topography features. It is concluded that the proposed method in this study can lead to artificial spatially variable earthquake ground motions which can be readily reflect the amplification pattern of 2D triangular valleys.
Material and methods
In the first part of this paper, seismic response of a triangular valley is investigated through time history analysis conducted by using FLAC2D computer program. The geometry of the valley analyzed in this paper is chosen close to the valley of the Masjed Soleyman embankment dam. Dynamic analysis is conducted using an artificial earthquake generated by spectral representation method. The material properties are obtained based on the results of a comprehensive study carried out to identify the dynamic characteristics of two large embankment dams in Iran. Spectral amplification functions of seismic waves are calculated by dividing the response spectra of stations located on the slope of the valley to that in base of the valley. These functions are then used as target quantity for generation of spatially variable ground motions at points located on the valley. In this study, spectral representation method, the most widely accepted method for generation of spatially variable ground motions, is developed to take into account the topography effect. According to this methodology, the generated time histories are compatible with prescribed spectral amplification functions reflecting the wave passage and loss of coherence effects. The Harichandran-Vanmarcke coherency model is used to simulate spatially variable seismic ground motions.
Results and discussion
Based on the obtained results the maximum and minimum values of peak acceleration are yielded at the base and at the edge of the valley, respectively. The results indicate considerable increase of the acceleration RMS at points near the edge of the valley. Maximum spectral amplification is also observed at the edge of the valley. For all points located on the valley, the first peak spectral amplification occurred at frequency of 1.15Hz, which can be readily interpreted as the natural frequency of the valley. In order to evaluate the accuracy of the proposed method, the RMS and spectral amplification functions of artificial earthquakes are compared to target quantities. A very good consistency between the spectral amplification of artificial earthquakes and target spectral amplifications was observed in terms of both amplitude and frequency content.
Conclusion
The following conclusions were drawn from this paper.
- Artificial earthquakes generated using proposed method of this paper are in a very good agreement with the amplification pattern of the valley.
- The results of this study can be readily used to investigate the influence of spatial variability of earthquake ground motion on structures like bridges and dams whose supports are located inside the valley and are subjected to multi-support earthquake excitation.
- The proposed method of this paper is not limited only to the valley topography, but it can be effectively used in the generation process of non - uniform artificial earthquakes for stations located on other topography features. The latter can be carried out by establishing the spectral amplification functions of other topography features such as slopes and hills resulted from field or numerical studies.

---

Introduction                               
Generally, in engineering geology physical and mechanical properties of rocks are investigated in macroscopic scale, and less attention is paid to investigate the texture and microstructure developing in rock during deformation. Salt rock, as a best example of ductile rocks, has attracted the attention of many researchers. Compared to silicate rocks, salt rock exhibits extensively ductile behavior at even low temperature and pressure. In micro-tectonics, salt is important, because of it is useful as an analogue material for understanding the microstructural processes and textural development in silicate rocks. Deformed salt rock can display microstructures developed in silicate rocks at high pressures and temperatures. Regarding the similarity between microstructures of salt rock and silicate rocks, investigation of microstructure and deformation mechanism in salt rock can be helpful in understanding the main cause of the squeezing phenomenon in tunnels.
One of the effective factors on squeezing phenomenon is the structures and microstructures of rock. Rock mass classifications that contain rock mass structures are used in the predicting methods. But, so far, no attention has been paid to the role of rock microstructure in predicting the squeezing phenomenon.
This study is aimed to identify deformation mechanisms occurring in microscopic scale in rocks and lead to tunnel convergent in large scale. To achieve this goal, the microstructures in a naturally deformed Late Pre-Cambrian to Early Cambrian Hormuz salt rock from the active Deh Kuyeh salt fountain in Fars province were investigated using Electron Backscatter Diffraction (EBSD).
Materials and Methods
Deh Kuyeh salt diapir was located at about 27 km NE of Lar city. Salt samples were taken from top of the east and west glaciers (S1 and S2) and from the middle part of diapiric stem (sample S3). Raw samples were first cut dry into slabs (approximately 3´2 ´1 cm). Thin sections were prepared following the procedure of Schleder and Urai (2005) and Urai et al. (1987).
Halite crystallographic orientation data were collected using a Zeiss SIGMAVP FEGSEM. EBSD patterns were collected using an accelerating voltage of 30 kV, beam current of ~ 100 nA and a working distance of about 30 mm. Oxford instruments AZTEC software was used for data acquisition. EBSD large step size (50 mm) mapping was used to examine the overall microstructure in each sample. EBSD data were processed using HKL Channel 5 software.
Results and Discussion
All samples showed relatively similar microstructures. Samples comprise a small number of large grains in a matrix of smaller grains. Most grains were irregular in shape with lobate boundaries and internal distortion. Microstructural study revealed that the ductile flow of the salt was accommodated by dislocation creep and dynamic recrystallization. Salt grains show lattice distortion and a prevalence of low-angle boundaries that are evidence for dislocation creep and recovery processes. Misorientation analysis suggests that (110) <110> and (111) <110> slip systems are responsible for crystal plastic deformation of salt grains. Schmid factor analysis showed that stresses acting on inclined directions lead to the maximum activity of these slip systems.
The observed microstructures in the salt are comparable with the microstructures presented for schist samples from Himalaya region. The rock along Himalaya main trusts also showed evidence of dislocation creep and development of crystallographic preferred orientation. Hence, this article suggests that the rock type and its microstructures are the most important factors in occurrence of tunnel convergent.
Conclusions
This article proposes that deformation mechanisms occurring in micro-scale control the rock behavior in large scale. All rocks can behave as a ductile material depending on the temperature and pressure. In intrinsically ductile rocks like salt rock, presence of many active slip systems facilitates rock deformation under lower pressures and temperatures than silicate rocks. High tectonic stresses in shear zones lead to development of a strong shape preferred orientation and crystal preferred orientation in rocks. These microstructures facilitate rock deformation under stresses exiting in tunnels. It can be said that rock type and tectonic history of the area play the most important role in occurrence of squeezing phenomenon. Other factors such as current stress system in the area control deformation speed in tunnel. It seems investigating microstructures of rocks from tunnel route before and after excavation can be effective in identifying places with high possibility of squeezing.

---

Introduction
Soil dynamic properties are used to evaluate the dynamic response of soils at different strain levels in geotechnical engineering. The shear modulus (G) and damping ratio (D) are among the most important dynamic properties of soils. In general, the factors affecting the dynamic behavior of soils are divided into two categories: first; soil type and characteristics such as water content, void ratio and soil plasticity and second; parameters of loads applied on the soil such as the number of loading cycles, loading frequency and loading waveform .Therefore, it is widely believed that the dynamic response of soils partially depends on the characteristics of the load. In this paper, 1-g shaking table tests were employed to investigate the effect of loading waveform and frequency content on dynamic properties of dry sands. The response obtained from soil samples during loading with different frequencies, input accelerations and waveforms were used to generate hysteresis loops of tested samples at different strain amplitudes. Then, hysteresis loops were used to determine the damping ratio and shear modulus at different strain levels. Finally, the effects of loading frequency and waveform on the changes of each parameter (G and D) were investigated.
Material and methods
A hydraulic shaking table with a single degree of freedom, designed and constructed at the Crisis Management Center of Urmia University, was used to conduct the experiments. Firoozkuh No. 161 sand was used in all the experiments. The Firoozkuh sand gradation curve is similar to that of Toyoura sand. In this study, accelerometers were used to measure the acceleration of the input to the sample as well as to record the acceleration caused by the input excitation at different depths of the soil sample. The displacement transducers (LVDT sensors) were also used to measure linear displacement. Each soil sample was constructed using dry Firoozkuh sand poured uniformly into the container from four equal heights of 150 mm to reach a total height of 600 mm. During the compaction process, the accelerometers A₁, A₂, and A₃ were placed at a depth of 150, 300 and 450 mm with respect to the bottom of container. Also, one accelerometer, A₀, was attached rigidly to the container base to measure base acceleration. A displacement transducer (L₁) was placed on the soil surface at a height of 600 mm from the floor of the container to measure the vertical displacement of the surface of the soil. In this study, 42 shaking table tests were performed to study the effect of loading frequency and waveform on dynamic properties of dry sand. The test samples were subjected to rectangular, sinusoidal and triangular loading at frequencies of 0.5 to 9 Hz and at input acceleration of 0.1g and 0.3 g.
Results and discussion
Given the importance of G-γ and D-γ curves in dynamic analyses, the changes in shear modulus with shear strain has been studied. The results show that the shear modulus increases as the frequency increases in all cases, and this increase is observed at lower frequencies and increases with increasing frequency. On the other hand, the shear modulus decreases with increasing shear strain. At a constant testing frequency, soil samples subjected to the rectangular waveform exhibited the largest shear modulus while the samples subjected to the triangular waveform showed the least shear modulus. The shear modulus of the samples under the sinusoidal waveform is barely more than the shear modulus of samples under triangular waveform. Moreover, by increasing the shear strain, the shear modulus values ​​of samples with different waveforms have become closer and the waveform effect is reduced. As for the effect of input acceleration on the shear modulus , increasing the input acceleration increases the shear strain and consequently, decreases the shear modulus in all states (the values ​​of shear modulus in various frequencies and the waveforms under the input acceleration of 0.1 g are larger than the shear modulus values ​​under the input acceleration of 0.3g). In the case of the damping ratio, the results show that, in all cases, damping ratio increases with shear strain. At low strain levels, the damping ratio values at various frequencies and waveforms are low and yet very close. At higher strain levels, the increase in frequency increases the damping ratio. This increase is more significant at higher frequencies. Also, the effect of waveform on the damping ratio is more apparent at larger shear strains, and at such shear strain levels, soil samples under rectangular loading exhibit the largest damping ratio. The damping ratio associated with the sinusoidal and triangular loading are also close to each other and it is a slightly larger for sinusoidal loading. On the other hand, the damping ratio increases with input acceleration. In addition, the effect of increased input acceleration on the increase in the damping ratio is more obvious at higher frequencies mainly due to the increase in shear strain.
Conclusion
In the present study, the effects of loading frequency and waveform on the dynamic properties of dry sand were investigated using shaking table tests. The following conclusions were drawn:

• The shear modulus increases with frequency. The trend is more obvious at larger frequencies. The effect of loading frequency on the damping ratio of the soil at low levels of strain is negligible, and at relatively large strain levels, damping ratio increases with loading frequency.
• Soil samples exhibit the highest shear modulus and damping ratio under rectangular loading. Therefore, in all the tested frequencies and input accelerations, the values of G and D for the rectangular waveforms are greater than those of the sinusoidal and triangular waveforms. The shear modulus and damping ratio for the sinusoidal waveforms are marginally greater than those of triangular waveforms, yet one can consider them practically similar.
• In all cases, the shear strain increased by increasing the amplitude of the input acceleration, and as a result, the shear modulus decreased and the damping ratio increased../files/site1/files/142/2.pdf

---

Introduction
Drill-bit selection is one of the most important aspects of well planning due to the bearing it can have on the overall cost of the well. Bit selection in conventional and slightly inclined wells is a very delicate and complex process. In high angle and horizontal wells it is even more difficult. Historically, drilling engineers have selected bits on the basis of what has been worked well in the area and what has been determined to have the lowest cost run from offset bit records. Often the best bit records were not available for evaluation, because the best bit may not yet have been run, may have been run by a competitor or the engineer was new to the area. As a result the bit program was generally developed by trial and error and at significant additional costs for a large number of wells. In most cases the optimum program was never reached because there was nothing to predict that a bit selection change could further reduce the cost of the well. In this study, an alternative solution approaches using the concept of the power of data mining algorithms to solve the optimum bit program for a given field is proposed.
Material and methods
It has been considered an offset well to be drilled outside the known boundaries of a known field. For this purpose, the seventh well (X-7) of the same field was used as a verification point. The data was trained using the well log and rock bit data of six wells located in the field and the real well log data of well 7 was input as unknown data. These depths are selected based on reported rock bit program. When compared to the real data, it could be observed that the models (adaptive neuro fuzzy inference system, K-nearest neighbors, decision tree, Bayesian classification theory and association rules) estimates the formation hardness accurately. This minor discrepancy was also present with the company’s suggested rock bit program, which was based on the previous wells’ rock bit data.
Results and discussion
In this paper, data mining algorithms for optimum rock bit program estimation is proposed. The accuracy and efficiency of the developed data mining algorithms (adaptive neuro fuzzy inference system, K-nearest neighbors, decision tree, Bayesian classification theory and association rules) that requires sonic and neutron log data input was tested for several real and synthetic cases. In the case of a development? well to be drilled outside the known boundaries of a field the model estimated rock bits with properties that consider the formation hardness correctly but slightly underestimated further rock bit details. The models also produced reasonable rock bit programs for an advance well to be drilled within the known boundaries of a field and a wildcat well drilled in a nearby field with similar rock properties to the training field. Thus it was concluded that the developed adaptive neuro fuzzy inference system is suitable as a front-end system for rock bit selection that could help engineers in decision-making analysis.
Conclusion
Optimum bit selection is one of the important issues in drilling engineering. Usually, optimum bit selection is determined by the lowest cost per foot and is a function of bit cost and performance as well as penetration rate. Conventional optimum rock bit selection program involves development of computer programs created from mathematical models along with information from previously drilled wells in the same area. Based on the data gathered on a daily basis for each well drilled, the optimum drilling program may be modified and revised as unexpected problems arose. The approaches in this study uses the power of data mining algorithms to solve the optimum bit selection problem. In order to achieve this goal, adaptive neuro fuzzy inference system, K-nearest neighbors, decision tree, Bayesian classification theory and association rules were developed by training the models using real rock bit data for several wells in a carbonated field. The training of the basic models involved use of both gamma ray and sonic log data. After that the models were tested using various drilling scenarios in different lithologic units. It was observed that the adaptive neuro fuzzy inference system model has provided satisfactory results.
 
 

---

Introduction
Hydraulic fracturing is one of the most important stimulation methods for oil and gas reservoirs with low permeability. Various factors, such as in-situ stresses, joints and natural fractures of the formation, fluid rheology, mechanical properties of the formation, injection fluid flow rate and perforation operation, effect on the pressure and hydraulic fracture geometry. Many researchers have studied hydraulic fracturing behavior of rocks since decades ago. The researches have showed that hydraulic fracturing operations increase the production of oil wells by up to 30 percent and increase gas wells by 90 percent. Currently, this operation is performed on about 60% of all drilled wells.
Material and methods
In this research, for the experimental investigation of the hydraulic fracturing, considering the reservoir condition, 39 concrete cubic samples with 100 × 100 × 100 mm dimensions and 60 concrete cylindrical samples with a diameter of 54 mm and a height of 110 mm were constructed and the effect of the sample geometry and in-situ stress field on the fracture geometry, breakdown pressure, the pattern of crack propagation and finally the cross fractures in vertical wellbores were investigated.
Results and discussion
In cubic specimens under uniaxial stress, with increasing vertical stress to 8 MPa, first the breakdown pressure has been increased by about 132% and then with increasing vertical stress to 16 MPa, the breakdown pressure has been decreased by about 69%. In cylindrical specimens under uniaxial stress, with increasing vertical stress to 12 MPa, first the breakdown pressure increased by about 113% and then with increasing vertical stress to 16 MPa, the breakdown pressure decreased by about 6%. As the vertical stress increases to a certain limit, the pores and micro-cracks inside the sample close, thus the tensile strength and breakdown pressure increase. In the following, increasing vertical stress causes more small cracks to open and reduces the tensile strength of the rock. In cubic specimens under triaxial stresses, with increasing vertical stress, the breakdown pressure has been increased. Also, in cylindrical specimens under triaxial stresses as the maximum horizontal stress increased, the breakdown pressure increased.
Conclusion
The obtained results demonstrated that increasing the uniaxial stress in the vertical wellbore in both types of studied sample geometry, first the breakdown pressure increases and then from one boundary onwards, with increasing vertical stress, the breakdown pressure decreases. In cubic specimens under triaxial stress, with increasing vertical stress, the breakdown pressure increases. Additionally, in cylindrical specimens under triaxial stresses as the maximum horizontal stress increases, the breakdown pressure increases, so that, in the 8 MPa maximum horizontal stress, the breakdown pressure increases by about 81%../files/site1/files/151/3.pdf
 

---

The soil of the Arvand free zone in the north of Khorramshahr is fine cohesive and cannot be used in earth works. On the other hand, suitable materials for this purpose (coarse-grained soils) are located at the farther distances which a considerable cost requires. In this regard, it is trying to improve the soil with lime and furnace steel slag. This study is focused on improvement of the fine-grained soil by adding various contents of lime and furnace steel slag. For this purpose, after sampling and performance of compaction tests, different amounts of slag (10, 20 and 30% by weight of dry soil) and lime (2, 4 and 6% by weight of dry soil) were added to the soil and after curing for 28 days, the effect of additives on the physical and mechanical properties of soil was investigated by using several tests such as Atterberg limits, compaction, uniaxial compressive strength (UCS) and CBR as soaked and unsoaked. Based on USCS classification the study soil is CL, its plasticity index is about 25% and sulphate ion content is more than 0.5%. Experimental results show that by adding slag and lime at different contents to soil, mechanical properties of soil improve dramatically, so plastic index of soil decreased and UCS and CBR has been increased. Also, the maximum dry unit weight of soil increases and the optimum moisture content decreases. The test results also indicate that the effect of lime on soil is higher than slag and the effect of slag for less than 35% is not considerable, however the test result of unsoaked CBR show that the bearing of soil increase in the more than slag content 20% is significant. According to the previous studies, due to the relatively high sulphate ion content in the soil, the use of lime alone is inappropriate and the slag can only physically improve soil conditions but also chemically prevent the formation of large volume minerals (like Ettringite) by the reaction of lime with soil sulphate ion../files/site1/files/152/%D8%AD%D8%B3%DB%8C%D9%86_%D8%B2%D8%A7%D8%AF%D9%87.pdf

---