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.  相似文献   

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.  相似文献   

Embedment effects of flexible foundations on control of structures     

Fereidoun Amini  Masoud Shadlou 《Soil Dynamics and Earthquake Engineering》2011

  相似文献   

Crustal structure of Northwest Zagros (Kermanshah) and Central Iran (Yazd and Isfahan) using teleseismic Ps converted phases     

Narges Afsari  Forogh Sodoudi  Fataneh Taghizadeh Farahmand  Mohammad Reza Ghassemi 《Journal of Seismology》2011,15(2):341-353

Receiver functions are widely employed to detect P-to-S converted waves and are especially useful to image seismic discontinuities in the crust. In this study we used the P receiver function technique to investigate the velocity structure of the crust beneath the Northwest Zagros and Central Iran and map out the lateral variation of the Moho boundary within this area. Our dataset includes teleseismic data (M _b ≥ 5.5, epicentral distance from 30° to 95°) recorded at 12 three-component short-period stations of Kermanshah, Isfahan and Yazd telemetry seismic networks. Our results obtained from P receiver functions indicate clear Ps conversions at the Moho boundary. The Moho depths were firstly estimated from the delay time of the Moho converted phase relative to the direct P wave beneath each network. Then, we used the P receiver function inversion to find the properties of the Moho discontinuity such as depth and velocity contrast. Our results obtained from PRF are in good agreement with those obtained from the P receiver function modeling. We found an average Moho depth of about 42 km beneath the Northwest Zagros increasing toward the Sanandaj-Sirjan Metamorphic Zone and reaches 51 km, where two crusts (Zagros and Central Iran) are assumed to be superposed. The Moho depth decreases toward the Urmieh-Dokhtar Cenozoic volcanic belt and reaches 43 km beneath this area. We found a relatively flat Moho beneath the Central Iran where, the average crustal thickness is about 42 km. Our P receiver function modeling revealed a shear wave velocity of 3.6 km/s in the crust of Northwest Zagros and Central Iran increasing to 4.5 km/s beneath the Moho boundary. The average shear wave velocity in the crust of UDMA as SSZ is 3.6 km/s, which reaches to 4.0 km/s while in SSZ increases to 4.3 km/s beneath the Moho.  相似文献