The Dahuiyeh (Zarand) earthquake of 2005 February 22 in central Iran: reactivation of an intramountain reverse fault   总被引：3，自引：0，他引：3  

M. Talebian  J. Biggs  M. Bolourchi  A. Copley  A. Ghassemi  M. Ghorashi  J. Hollingsworth  J. Jackson  E. Nissen  B. Oveisi  B. Parsons  K. Priestley  A. Saiidi 《Geophysical Journal International》2006,164(1):137-148

  相似文献   

Large‐scale poroelastic fractured reservoirs modeling using the fast multipole displacement discontinuity method          下载免费PDF全文

A. Verde  A. Ghassemi 《国际地质力学数值与分析法杂志》2016,40(6):865-886

An effective approach to modeling the geomechanical behavior of the network and its permeability variation is to use a poroelastic displacement discontinuity method (DDM). However, the approach becomes rather computationally intensive for an extensive system of cracks, particularly when considering coupled diffusion/deformation processes. This is because of additional unknowns and the need for time‐marching schemes for the numerical integration. The Fast Multipole Method (FMM) is a technique that can accelerate the solution of large fracture problems with linear complexity with the number of unknowns both in memory and CPU time. Previous works combining DDM and FMM for large‐scale problems have accounted only for elastic rocks, neglecting the fluid leak‐off from the fractures into the matrix and its influence on pore pressure and stress field. In this work we develop an efficient geomechanical model for large‐scale natural fracture networks in poroelastic reservoirs with fracture flow in response to injection and production operations. Accuracy and computational performance of the proposed method with those of conventional poroelastic DDM are compared through several case studies involving up to several tens of thousands of boundary elements. The results show the effectiveness of the FMM approach to successfully evaluate field‐scale problems for the design of exploitation strategies in unconventional geothermal and petroleum reservoirs. An example considering faults reveals the impact of reservoir compartmentalization because of sealing faults for both geomechanical and flow variables under elastic and poroelastic rocks. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

Modeling 3D thermal fracture propagation by transient cooling using virtual multidimensional internal bonds          下载免费PDF全文

K. Huang  A. Ghassemi 《国际地质力学数值与分析法杂志》2016,40(17):2293-2311

Thermal fracturing can play an important role in development of unconventional petroleum and geothermal resources. Thermal fractures can result from the nonlinear deformation of the rock in response to thermal stress related to cold water injection as well as heating. Before the rock reaches the final failure stage, material softening and bulk modulus degradation can cause changes in the thermo‐mechanical properties of the solid. In order to capture this aspect of the rock fracture, a virtual multidimensional internal bond‐based thermo‐mechanical model is derived to track elastic, softening, and the failure stages of the rock in response to the temporal changes of its temperature field. The variations in thermo‐mechanical properties of the rock are derived from a nonlinear constitutive model. To represent the thermo‐mechanical behavior of pre‐existing fractures, the element partition method is employed. Using the model, numerical simulation of 3D thermal fracture propagation in brittle rock is carried out. Results of numerical simulations provide evidence of model verification and illustrate nonlinear thermal response and fracture development in rock under uniform cooling. In addition, fracture coalescence in a cluster of fractures under thermal stress is illustrated, and the process of thermal fracturing from a wellbore is captured. Results underscore the importance of thermal stress in reservoir stimulation and show the effectiveness of the model to predict 3D thermal fracturing. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

Modeling three‐dimensional hydraulic fracture propagation using virtual multidimensional internal bonds     

K. Huang  Z. Zhang  Ahmad Ghassemi 《国际地质力学数值与分析法杂志》2013,37(13):2021-2038

Propagation of fractures, especially those emanating from wellbores and closed natural fractures, often involves Mode I and Mode II, and at times Mode III, posing significant challenges to its numerical simulation. When an embedded inclined fracture is subjected to compression, the fracture edge is constrained by the surrounding materials so that its true propagation pattern cannot be simulated by 2D models. In this article, a virtual multidimensional internal bond (VMIB) model is presented to simulate three‐dimensional (3D) fracture propagation. The VMIB model bridges the processes of macro fracture and micro bond rupture. The macro 3D constitutive relation in VMIB is derived from the 1D bond in the micro scale and is implemented in a 3D finite element method. To represent the contact and friction between fracture surfaces, a 3D element partition method is employed. The model is applied to simulate fracture propagation and coalescence in typical laboratory experiments and is used to analyze the propagation of an embedded fracture. Simulation results for single and multiple fractures illustrate 3D features of the tensile and compressive fracture propagation, especially the propagation of a Mode III fracture. The results match well with the experimental observation, suggesting that the presented method can capture the main features of 3D fracture propagation and coalescence. Moreover, by developing an algorithm for applying pressure on the fracture surfaces, propagation of a natural fracture is also simulated. The result illustrates an interesting and important phenomenon of Mode III fracture propagation, namely the fracture front segmentation. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

Simulation of hydraulic fracture propagation near a natural fracture using virtual multidimensional internal bonds     

Z. Zhang  A. Ghassemi 《国际地质力学数值与分析法杂志》2011,35(4):480-495

A virtual multidimensional internal bond (VMIB) model developed to simulate the propagation of hydraulic fractures using the finite‐element method is formulated within the framework of the virtual internal bond theory (VIB) that considers a solid as randomized material particles in the micro scale, and derives the macro constitutive relation from the cohesive law between the material particles with an implicit fracture criterion. Hydraulic pressure is applied using a new scheme that enables simulation of hydraulically driven cracks. When the model is applied to study hydraulic fracture propagation in the presence of a natural fracture, the results show the method to be very effective. It shows that although the in situ stress ratio is the dominant factor governing the propagation direction, a natural fault can also strongly influence the hydraulic fracture behavior. This influence is conditioned by the shear stiffness of the fault and the distance to the original hydraulic fracture. The model results show that when the fault is strong in shear, its impact on hydraulic fracture trajectory is weak and the hydraulic fracture will likely penetrate the fault. For a weak fault, however, the fracture tends to be arrested at the natural fault. The distance between the fault and the hydraulic fracture is also important; the fault influence increases with decreasing distance. The VMIB does not require selection of a fracture criterion and remeshing when the fracture propagates. Therefore, it is advantageous for modeling fracture initiation and propagation in naturally fractured rock. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

Pore scale study of permeability and tortuosity for flow through particulate media using Lattice Boltzmann method     

Ali Ghassemi  Ali Pak 《国际地质力学数值与分析法杂志》2011,35(8):886-901

In this paper, Lattice Boltzmann method (LBM) has been used to study the effects of permeability and tortuosity on flow through saturated particulate media and identify the relationships between permeability and tortuosity with other parameters such as particles diameter, grain specific surface, and porosity. LBM is a simple kinematic model that can incorporate the essential physics of microscopic and mesoscopic processes involved in flow through granular soils. The obtained results indicate that the 2D LB model, due to its inherent theoretical advantages, is capable of demonstrating that the porosity and specific surface are the most influential parameters in determining the intrinsic permeability of granular media. The obtained results show that particles' size diameter has a two‐fold effect on the coefficient of permeability: one is through specific surface and the other is by tortuosity factor. Numerical study also reveals that tortuosity of granular soils decreases almost linearly with increasing the porosity. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

Mathematical modelling of sea water intrusion,Nauru Island     

F. Ghassemi  A. J. Jakeman  G. Jacobson 《水文研究》1990,4(3):269-281

A two-dimensional finite element model for density dependent groundwater flow was calibrated to simulate sea water intrusion in Nauru Island in the Central Pacific Ocean. Nauru Island occupies an area of 22 km² and supports a population of 8500. The island has been mined for its phosphate deposits and current reserves indicate that the mine has about eight years life remaining. The water supply of the island is about one third dependent on imported water which is also used as ballast on the phosphate ships. Imported water will not be available in the future, and a hydrogeological investigation shows that the island is underlain by a fresh water layer, less than 5 m thick. The freshwater layer overlies a thick transition zone of brackish water which in turn overlies sea water. Simulation of several management options shows that it is possible to substitute current importation of fresh water by careful extraction from the groundwater resources of the island.  相似文献   

Static and dynamic analysis on slope stability using a DFN-DEM approach on the right abutment of the Karun 4 dam     

Kefayati  Sadegh  Baghbanan  Alireza  Torkan  Masoud  Hashemolhosseini  Hamid  Dehnavi  Roohollah Narimani 《地震工程与工程振动(英文版)》2020,19(4):937-951

Uncertainty in input fracture geometric parameters during analysis of the stability of jointed rock slopes is inevitable and therefore the stochastic discrete fracture network (DFN) — distinct element method (DEM) is an efficient modeling tool. In this research, potentially unstable conditions are detected in the right abutment of the Karun 4 dam and downstream of the dam body as a case study. Two extreme states with small and relatively large block sizes are selected and a series of numerical DEM models are generated using a number of validated DFN models. Stability of the rock slope is assessed in both static and dynamic loading states. Based on the design basis earthquake (DBE) and maximum credible earthquake (MCE) expected in the dam site, histories of seismic waves are applied to analyze the stability of the slope in dynamic earthquake conditions. The results indicate that a MCE is likely to trigger sliding of rock blocks on the rock slope major joint. Furthermore, the dynamic analysis also shows a local block failure by the DBE, which can consequently lead to slope instability over the long term. According to the seismic behavior of the two models, larger blocks are prone to greater instability and are less safe against earthquakes.

  相似文献   

Lithospheric structure of NW Iran from P and S receiver functions     

Fataneh Taghizadeh-Farahmand  Forough Sodoudi  Narges Afsari  Mohammad R. Ghassemi 《Journal of Seismology》2010,14(4):823-836

We computed P and S receiver functions to investigate the lithospheric structure beneath the northwest Iran and compute the Vp/Vs ratio within the crust of this seismologically active area. Our results enabled us to map the lateral variations of the Moho as well as those of the lithosphere–asthenosphere boundary (LAB) beneath this region. We selected data from teleseismic events (M_b?>?5.5, epicentral distance between 30° and 95° for P receiver functions and M_b?>?5.7, epicentral distance between 60° and 85° for S receiver functions) recorded from 1995 to 2008 at 8 three-component short-period stations of Tabriz Telemetry Seismic Network. Our results obtained from P receiver functions indicate clear conversions at the Moho boundary. The Moho depth was firstly estimated from the delay time of the Moho converted phase relative to the direct P wave. Then we used the H-Vp/Vs stacking algorithm of Zhu and Kanamori to estimate the crustal thickness and Vp/Vs ratio underneath the stations with clear Moho multiples. We found an average Moho depth of 48 km, which varies between 38.5 and 53 km. The Moho boundary showed a significant deepening towards east and north. This may reveal a crustal thickening towards northeast possibly due to the collision between the Central Iran and South Caspian plates. The obtained average Vp/Vs ratio was estimated to be 1.76, which varies between 1.73 and 1.82. The crustal structure was also determined by modeling of P receiver functions. We obtained a three-layered model for the crust beneath this area. The thickness of the layers is estimated to be 6–11, 18–35, and 38–53 km, respectively. The average of the shear wave velocity was calculated to be 3.4 km/s in the crust and reaches 4.3 km/s below the Moho discontinuity. The crustal thickness values obtained from P receiver functions are in good agreement with those derived by S receiver functions. In addition, clear conversions with negative polarity were observed at ~8.7 s in S receiver functions, which could be related to the conversion at the LAB. This may show a relatively thin continental lithosphere of about 85 km implying that the lithosphere was influenced by various geodynamical reworking processes in the past.  相似文献   

Application of modified qualitative index for surveillance of water-filtration process in turbidity removal by different media   总被引：1，自引：1，他引：0  

G. Badalians Gholikandi  M. Noorisepehr  E. Dehghanifard  A. Koolivand  A. Dehnavi  S. Moalej 《International Journal of Environmental Science and Technology》2012,9(4):691-700

Several media have been used in treatment plants, however, their efficiency for turbidity removal, which is determined by qualitative indices, has been considered. Current qualitative indices such as turbidity and escaping particle number could not completely measure the efficiency of the filtration system; therefore defining new qualitative indices is essential. In this study, the efficiency of two different dual media filters in turbidity removal was compared in different operating condition using qualitative indices. The pilot consisted of a filter column (1-m depth) in which the filter-1 was consisted of a layer of anthracite (450-mm depth) and a layer of silica sand (350-mm depth); and filter-2 had the same media characteristics except for the first layer that was light expanded clay aggregates (LECA). Turbidities of 10, 20, and 30 NTU, coagulant concentrations of 4, 8, and 12?ppm and filtration rates of 10, 15, and 20?m/h were considered as variables. Results showed that the media of filter-2 is a suitable substitute for the media of filter-1 (P value?<?0.05). Turbidity removal efficiencies in different condition were 79.97?±?1.79 to 91.37?±?1.23% for the filter-2 and 75.12?±?2.75 to 86.82?±?1.3% for the filter-1. The LECA layer efficiency in turbidity removal was independent of filtration rates and due to its low head loss; LECA can be used as a proper medium. Results also showed that the particle index was a suitable index as a substitute for turbidity and escaping particle number as indices.  相似文献