共查询到20条相似文献，搜索用时 31 毫秒
1.
Jonny Rutqvist Deborah Barr Jens T. Birkholzer Kiyoshi Fujisaki Olaf Kolditz Quan-Sheng Liu Tomoo Fujita Wenqing Wang Cheng-Yuan Zhang 《Environmental Geology》2009,57(6):1347-1360
This paper presents an international, multiple-code, simulation study of coupled thermal, hydrological, and mechanical (THM)
processes and their effect on permeability and fluid flow in fractured rock around heated underground nuclear waste emplacement
drifts. Simulations were conducted considering two types of repository settings (1) open emplacement drifts in relatively
shallow unsaturated volcanic rock, and (2) backfilled emplacement drifts in deeper saturated crystalline rock. The results
showed that for the two assumed repository settings, the dominant mechanism of changes in rock permeability was thermal–mechanically
induced closure (reduced aperture) of vertical fractures, caused by thermal stress resulting from repository-wide heating
of the rock mass. The magnitude of thermal–mechanically induced changes in permeability was more substantial in the case of
an emplacement drift located in a relatively shallow, low-stress environment where the rock is more compliant, allowing more
substantial fracture closure during thermal stressing. However, in both of the assumed repository settings in this study,
the thermal–mechanically induced changes in permeability caused relatively small changes in the flow field, with most changes
occurring in the vicinity of the emplacement drifts. 相似文献
2.
Coupled THM processes in EDZ of crystalline rocks using an elasto-plastic cellular automaton 总被引：5，自引：0，他引：5
Peng-Zhi Pan Xia-Ting Feng Xiao-Hua Huang Qiang Cui Hui Zhou 《Environmental Geology》2009,57(6):1299-1311
This paper aims at a numerical study of coupled thermal, hydrological and mechanical processes in the excavation disturbed
zones (EDZ) around nuclear waste emplacement drifts in fractured crystalline rocks. The study was conducted for two model
domains close to an emplacement tunnel; (1) a near-field domain and (2) a smaller wall-block domain. Goodman element and weak
element were used to represent the fractures in the rock mass and the rock matrix was represented as elasto-visco-plastic
material. Mohr–Coulomb criterion and a non-associated plastic flow rule were adopted to consider the viscoplastic deformation
in the EDZ. A relation between volumetric strain and permeability was established. Using a self-developed EPCA2D code, the elastic, elasto-plastic and creep analyses to study the evolution of stress and deformations, as well as failure
and permeability evolution in the EDZ were conducted. Results indicate a strong impact of fractures, plastic deformation and
time effects on the behavior of EDZ especially the evolution of permeability around the drift. 相似文献
3.
Thermal Constraints on the Emplacement Rate of a Large Intrusive Complex: The Manaslu Leucogranite, Nepal Himalaya 总被引：4，自引：0，他引：4
The emplacement of the Manaslu leucogranite body (Nepal, Himalaya)has been modelled as the accretion of successive sills. Theleucogranite is characterized by isotopic heterogeneities suggestinglimited magma convection, and by a thin (<100 m) upper thermalaureole. These characteristics were used to constrain the maximummagma emplacement rate. Models were tested with sills injectedregularly over the whole duration of emplacement and with twoemplacement sequences separated by a repose period. Additionally,the hypothesis of a tectonic top contact, with unroofing limitingheat transfer during magma emplacement, was evaluated. In thislatter case, the upper limit for the emplacement rate was estimatedat 3·4 mm/year (or 1·5 Myr for 5 km of granite).Geological and thermobarometric data, however, argue againsta major role of fault activity in magma cooling during the leucograniteemplacement. The best model in agreement with available geochronologicaldata suggests an emplacement rate of 1 mm/year for a relativelyshallow level of emplacement (granite top at 10 km), uninterruptedby a long repose period. The thermal aureole temperature andthickness, and the isotopic heterogeneities within the leucogranite,can be explained by the accretion of 20–60 m thick sillsintruded every 20 000–60 000 years over a period of 5Myr. Under such conditions, the thermal effects of granite intrusionon the underlying rocks appear limited and cannot be invokedas a cause for the formation of migmatites. KEY WORDS: granite emplacement; heat transfer modelling; High Himalayan Leucogranite; Manaslu; thermal aureole 相似文献
4.
为研究填充裂隙水流速度对岩体温度和应力的影响，选取甘肃北山地区的花岗岩，制作了稀疏非正交裂隙岩体模
型，采用河砂填充裂隙后进行模型试验；并对模型试验进行离散元数值模拟，分析了模型试验所测岩体温度和热应力与数
值模拟结果的差异和原因。结果显示：填砂裂隙强化了裂隙介质的热导能力，无填充时岩体温度和热应力比填砂时低；模
型试验和离散元模拟均表明，岩体温度和岩体应力随裂隙水流速度增大而减小，但是系统达到稳态所需要的时间变短；模
型试验中斜裂隙水流对温度场起主要作用；由于现有3DEC软件不能考虑水的热物性参数随温度的变化，进而产生自然对
流换热，斜裂隙水流和靠近热源侧的竖裂隙水流对温度场起主要作用。 相似文献
5.
为研究填充裂隙水流速度对岩体温度和应力的影响,选取甘肃北山地区的花岗岩,制作了稀疏非正交裂隙岩体模
型,采用河砂填充裂隙后进行模型试验;并对模型试验进行离散元数值模拟,分析了模型试验所测岩体温度和热应力与数
值模拟结果的差异和原因。结果显示：填砂裂隙强化了裂隙介质的热导能力,无填充时岩体温度和热应力比填砂时低;模
型试验和离散元模拟均表明,岩体温度和岩体应力随裂隙水流速度增大而减小,但是系统达到稳态所需要的时间变短;模
型试验中斜裂隙水流对温度场起主要作用;由于现有3DEC软件不能考虑水的热物性参数随温度的变化,进而产生自然对
流换热,斜裂隙水流和靠近热源侧的竖裂隙水流对温度场起主要作用。 相似文献
6.
The analysis of elasto-plastic boundary value problems using the finite element method involves many discretizations. These lead to the problem of yield surface drift in which the stress state predicted at the end of an elasto-plastic increment of loading does not lie on the current yield surface. As such discrepancies are comulative it is important to ensure that the stresses are corrected back to the yield surface during each increment of loading. In this paper five methods of accounting for this drift are examined. These involve correcting the stresses by projecting back along the plastic flow, the total strain increment and the accumulated effective stress direction. In addition a ‘correct’, method which accounts for the changes in elastic strains which accompany any stress correction is considered. This method is theoretically more sound than the other approximate approaches. All five methods have been used in finite element analyses of the stress changes that occur adjacent to a single pile installed in a uniform deposit of soil on pile loading. The soil was assumed to be normally consolidated and was modelled using a form of modified Cam Clay. Comparison of these results with an analysis, in which yield surface drift was negligible indicated that only the ‘correct’ method and the method involving projecting back along the plastic flow direction give accurate predictions. Substantiai errors occur if the other methods of correcting for yield surface drift are employed. It is recommended that the ‘correct’ method be adopted for finite element calculations. 相似文献
7.
Thermal oil recovery processes involve high pressures and temperatures, leading to large volume changes and induced stresses. These cannot be handled by traditional reservoir simulation because it does not consider coupled geomechanics effects. In this paper we present a fully coupled, thermal half‐space model using a hybrid DDFEM method. A finite element method (FEM) solution is adopted for the reservoir and the surrounding thermally affected zone, and a displacement discontinuity method is used for the surrounding elastic, non‐thermal zone. This approach analyzes stress, pressure, temperature and volume change in the reservoir; it also provides stresses and displacements around the reservoir (including transient ground surface movements) in a natural manner without introducing extra spatial discretization outside the FEM zone. To overcome spurious spatial temperature oscillations in the convection‐dominated thermal advection–diffusion problem, we place the transient problem into an advection–diffusion–reaction problem framework, which is then efficiently addressed by a stabilized finite element approach, the subgrid‐scale/gradient subgrid‐scale method. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
8.
Margaret M. Dekshenieks Eileen E. Hofmann John M. Klinck Eric N. Powell 《Estuaries and Coasts》2000,23(5):593-610
Three models are combined to investigate the effects of changes in environmental conditions on the population structure of the Eastern oyster,Crassostrea virginica. The first model, a time-dependent model of the oyster population as described in Powell et al. (1992, 1994, 1995a,b, 1996, 1997) and Hofmann et al. (1992, 1994, 1995), tracks the distribution, development, spawning, and mortality of sessile oyster populations. The second model, a time-dependent larval growth model as described in Dekshenieks et al. (1993), simulates larval growth and mortality. The final model, a finite element hydrodynamic model, simulates the circulation in Galveston Bay, Texas. The coupled post-settlement-larval model (the oyster model) runs within the finite element grid at locations that include known oyster reef habitats. The oyster model was first forced with 5 yr of mean environmental conditions to provide a reference simulation for Galveston Bay. Additional simulations considered the effects of long-term increases and decreases in freshwater inflow and temperature, as well as decreases in food concentration and total seston on Galveston Bay oyster populations. In general, the simulations show that salinity is the primary environmental factor controling the spatial extent of oyster distribution within the estuary. Results also indicate a need to consider all environmental factors when attempting to predict the response of oyster populations; it is the superposition of a combination of these factors that determines the state of the population. The results from this study allow predictions to be made concerning the effects of environmental change on the status of oyster populations, both within Galveston Bay and within other estuarine systems supporting oyster populations. 相似文献
9.
H.-R. Wenk M. Bortolotti N. Barton E. Oliver D. Brown 《Physics and Chemistry of Minerals》2007,34(9):599-607
Mechanical twinning in polycrystalline quartz was investigated in situ with time-of-flight neutron diffraction and a strain
diffractometer. Dauphiné twinning is highly temperature sensitive. It initiates at a macroscopic differential stress of 50–100 MPa
and, at 500°C, saturates at 400 MPa. From normalized diffraction intensities the patterns of preferred orientation (or texture)
can be inferred. They indicate a partial reversal of twinning during unloading. The remaining twins impose residual stresses
corresponding to elastic strains of 300–400 microstrain. Progressive twinning on loading and reversal during unloading, as
well as the temperature dependence, can be reproduced with finite element model simulations. 相似文献
10.
Shunde Yin 《国际地质力学数值与分析法杂志》2013,37(15):2523-2538
Integration of poromechanics and fracture mechanics plays an important role in understanding a series of thermal fracturing phenomena in subsurface porous media such as cold water flooding for enhanced oil recovery, produced‐water reinjection for waste disposal, cold water injection for geothermal energy extraction, and CO2 injection for geosequestration. Thermal fracturing modeling is important to prevent the potential risks when fractures propagate into undesired zones, and it involves the coupling of heat transfer, mass transport, and stress change as well as the fracture propagation. Analytical method, finite element method, and finite difference method as well as boundary element method have been used to perform the thermal fracturing modeling considering different degrees and combinations of coupling. In this paper, extended finite element method is employed for the thermal fracturing modeling in a fully coupled fashion with remeshing avoided, and the stabilized finite element method is employed to account for the convection‐dominated heat transfer in the fracturing process with numerical oscillation circumvented. With the thermal fracturing model, a hypothetical numerical experiment on cold water injection into a deep warm aquifer is conducted. Results show that parameters such as injection rate, injection temperature, aquifer stiffness, and permeability can affect the fracture development in different ways and extended finite element method and stabilized finite element method provide effective tools for thermal fracturing simulation. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
11.
《Applied Geochemistry》2006,21(2):203-222
Disposal of nuclear waste in deep geological formations is expected to induce thermal fluxes for hundreds of years with maximum temperature reaching about 100–150 °C in the nearfield argillaceous environment. The long-term behavior of clays subjected to such thermal gradients needs to be perfectly understood in safety assessment considerations. In this respect, a Toarcian argillaceous unit thermally disturbed by the intrusion of a 1.1-m wide basaltic dike at the Perthus pass (Herault, France), was studied in detail as a natural analogue. The thermal imprint induced by the dike was evaluated by a mineralogical, chemical and K–Ar study of the <2 μm clay fraction of shale samples collected at increasing distance from the basalt. The data suggest that the mineral composition of the shales was not significantly disturbed when the temperature was below 100–150 °C. Closer to the dike at 150–300 °C, changes such as progressive dissolution of chlorite and kaolinite, increased content of the mixed layers illite–smectite with more illite layers, complete decalcification and subsequent increased content of quartz, were found.At the eastern contact with the dike, the mineral and chemical compositions of both the shales and the basalt suggest water–rock interactions subsequent to the intrusion with precipitation of palagonite and renewed but discrete deposition of carbonate. A pencil cleavage developed in the shales during the dike emplacement probably favored water circulation along the contact. Strontium isotopic data suggest that the fluids of probable meteoric origin, reacted with Bathonian and Bajocian limestones before entering the underlying Toarcian shales.By analogy with deep geological radioactive waste repositories, the results report discrete mineralogical variations of the clays when subjected to temperatures of 100–150 °C that are expected in deep storage conditions. Beyond 150 °C, significant mineralogical changes may alter the physical and chemical properties of the shales, especially of the clay fraction. Also, the development of structural discontinuities in the so-called thermally disturbed zone might be of importance as these discontinuities might become zones for preferential fluid circulation. Finally, the study emphasizes the use of Rb–Sr and K–Ar isotopic systems as tracers of local circulating fluids related to low-grade thermal imprints. 相似文献
12.
Methane hydrate (MH, also called fiery ice) exists in forms of pore filling, cementing and load-bearing skeleton in the methane hydrate bearing sediment (MHBS) and affects its mechanical behavior greatly. To study the changes of macro-scale and micro-scale mechanical behaviors of MHBS during exploitation by thermal recovery and depressurization methods, a novel 2D thermo-hydro-mechanical bonded contact model was proposed and implemented into a platform of distinct element method (DEM), PFC2D. MHBS samples were first biaxially compressed to different deviator stress levels to model different in-situ stress conditions. With the deviator stress maintained at constant, the temperature was then raised to simulate the thermal recovery process or the pore water pressure (i.e. confining pressure for MH bond) was decreased to simulate the depressurization process. DEM simulation results showed that: during exploitation, the axial strain increased with the increase of temperature (in the thermal recovery method) or decrease of pore water pressure (in the depressurization method); sample collapsed during MH dissociation if the deviator stress applied was larger than the compression strength of a pure host sand sample; sample experienced volume contraction but its void ratio was slightly larger than the pure host sand sample at the same axial strain throughout the test. By comparison with the laboratory test results, the new model was validated to be capable of reproducing the exploitation process by thermal recovery and depressurization methods. In addition, some micro-scale parameters, such as contact distribution, bond distribution, and averaged pure rotation rate, were also analyzed to investigate their relationships with the macroscopic responses. 相似文献
13.
Seyed Erfan Saberhosseini Reza Keshavarzi Kaveh Ahangari 《Arabian Journal of Geosciences》2017,10(6):152
One of the crucial consequences of steam assisted gravity drainage (SAGD) process is abnormal reservoir uplifting under thermal steam injection, which can significantly influence the reservoir rock deformation, specifically thin bed reservoirs and causes intensive failures and fractures into the cap rock formations. A thorough understanding of the influences of rock thermo-mechanical properties on reservoir uplifting plays an important role in preventing those aforementioned failures within design and optimization process in SAGD. In addition, coupling of reservoir porous medium and flowing of specific fluid with temperature as an additional degree of freedom with initial pore pressure and in-situ stress condition, are also very challenging parts of geomechanical coupled simulation which would be clearly explained. Thus, a fully coupled thermo-poro-elastic geomechanical model with finite element codes was performed in ABAQUS to investigate the role of rock thermo-mechanical parameters on reservoir vertical uplift during steam injection. It is clearly observed that, any increase in rock thermo-mechanical properties specifically rock’s thermal properties such as specific heat, thermal expansion, and formation’s thermal conductivity, have significant influences on reservoir uplift. So by coupling the temperature as an additional degree of freedom with the coupled pore-fluid stress and diffusion finite element model of SAGD process, the more realistic simulation will be conducted; hence, the errors related to not having heat as an additional degree of freedom will be diminished. In addition, Young’s modulus and specific heat are the rock thermo-mechanical parameters which have the maximum and minimum effects on the reservoir uplift, respectively. 相似文献
14.
15.
16.
The mid-term safety of deep geological nuclear waste repositories is based in part on the presence of a buffer, the main role of which is to isolate the environment from radionuclides. A design evaluation of such a repository is necessary to assess the potential vertical canister movement inside the drift that could reduce the buffer efficiency. A thermo-hydro-mechanical (THM) simulation is performed in a vertical cross-section of the drift. The THM couplings are described, and their influences on the mid-term (300 years) response of the engineered barrier system (EBS) are revealed. This study uses an advanced constitutive model to simulate the THM processes that occur in a specific EBS design case. The near-field simulations of the nuclear waste canister are performed in a two-dimensional finite element configuration that considers the effect of gravity. The focus of this study is on the mechanical behaviour of the buffer, which consists of two different forms of bentonite. Such an approach allows realistic consideration of the effect of the wetting and drying of the buffer material in non-isothermal conditions. Due to a specific design that includes bentonite blocks and pellets, the canister is observed to heave slightly during the re-saturation period, which extends up to 100 years. 相似文献
17.
Obduction emplaces regional-scale fragments of oceanic lithosphere (ophiolites) over continental lithosphere margins of much lower density. For this reason, the mechanisms responsible for obduction remain enigmatic in the framework of plate tectonics. We present two-dimensional (2D) thermo-mechanical models of obduction and investigate the possible dynamics and physical controls of this process. Model geometry and boundary conditions are based on available geological and geochronological data and numerical modeling results are validated against petrological and structural observations of the Oman (Semail) Ophiolite. Our model reproduces the stages of oceanic subduction initiation away from the Arabian margin, the emplacement of the Oman Ophiolite on top of it, and the domal exhumation of the metamorphosed margin through the ophiolitic nappe. A systematic study indicates that 350–400 km of bulk shortening provides the best fit for both maximum pressure–temperature conditions of the metamorphosed margin (1.5–2.5 GPa/450–600 °C) and the dimension of the ophiolitic nappe (~ 170 km width). Our results confirm that a thermal anomaly located close to the Arabian margin (~ 100 km) is needed to initiate obduction. We further suggest that a strong continental basement rheology is a prerequisite for ophiolite emplacement. 相似文献
18.
Based on the Crust2.0 model and the topography data of Chinese continent and its adjacent regions, a three-dimensional finite element model is constructed in terms of the spherical coordinate system. In our numerical model, the average annual ground temperature from 195 meteorological stations and temperature of upper mantle derived from the seismic velocities are adopted as the top and bottom boundary conditions, respectively. The observed thermal conductivity and heat production from P wave velocity based on empirical formula are employed in our numerical model as well. The comparison between the calculated and observed surface heat flow proved that our results are reliable. The temperature beneath the Precambrian cratons is lower than that of other areas for 100–300 °C also. The typical temperature rang at the Moho is estimated to be 800–1000 °C beneath the Tibetan plateau and 500–700 °C beneath the Precambrian cratons (such as Indian plate, Sichuan basin, South China, North China and Tarim), respectively. The thermal state in the eastern part of Sino-Korean craton at the depth deeper than 60 km indicates that it was destructed. The thermal structure in center of Tibetan plateau (especially beneath Qiangtang area) supports the proposed flow of lower crustal or upper mantle material to the east. Generally, the distribution of volcanoes in Chinese continent is consistent with the high temperature areas in the crust or upper mantle. There are many obvious thermal transition zones across the orogenic belts. The thermal transition zone between eastern and western parts in the crust of Chinese continent is consistent with the north–south seismic zone. 相似文献
19.
In this study, short-term surface settlements are predicted for twin tunnels, which are to be excavated in the chainage of
0 + 850 to 0 + 900 m between the Esenler and Kirazlı stations of the Istanbul Metro line, which is 4 km in length. The total
length of the excavation line is 21.2 km between Esenler and Basaksehir. Tunnels are excavated by employing two earth pressure
balance (EPB) tunnel boring machines (TBMs) that have twin tubes of 6.5 m diameter and with 14 m distance from center to center.
The TBM in the right tube follows about 100 m behind the other tube. Segmental lining of 1.4 m length is currently employed
as the final support. Settlement predictions are performed with finite element method by using Plaxis finite element program.
Excavation, ground support and face support steps in FEM analyses are simulated as applied in the field. Predictions are performed
for a typical geological zone, which is considered as critical in terms of surface settlement. Geology in the study area is
composed of fill, very stiff clay, dense sand, very dense sand and hard clay, respectively, starting from the surface. In
addition to finite element modeling, the surface settlements are also predicted by using semi-theoretical (semi-empirical)
and analytical methods. The results indicate that the FE model predicts well the short-term surface settlements for a given
volume loss value. The results of semi-theoretical and analytical methods are found to be in good agreement with the FE model.
The results of predictions are compared and verified by field measurements. It is suggested that grouting of the excavation
void should be performed as fast as possible after excavation of a section as a precaution against surface settlements during
excavation. Face pressure of the TBMs should be closely monitored and adjusted for different zones. 相似文献
20.
The Nisa-Alburquerque granitic batholith (southern Variscan Belt, Iberian Peninsula) has been studied by petrological, structural
and geophysical approaches, obtaining contrasting models for its deep structure and emplacement sequence. In order to test
these models and gain knowledge on the thermal increase induced by the intrusion, we have studied its contact aureole, which
was developed in similar country rock lithologies (mica schists alternating with metasandstones and feldespatic schists) all
along the northern external contact of the batholith. Our results indicate no change in metamorphic grade and some variations
in aureole width, which narrows toward the western sectors of the batholith. Cordierite is the only contact metamorphic mineral
developed together with a high temperature biotite probably related to the granite thermal input. By considering these new
data, together with zircon saturation temperatures within the granite and previous petrological and geophysical studies, we
propose a model in which the feeder zones of the granitic magmas were an eastern main one and a western secondary one. We
have also made comparisons of the metamorphic grade in the country rocks and the xenoliths within the granite. Most of the
xenoliths have the same metamorphic facies as the country rocks (Crd-zone), though some of them contain slightly different
assemblages (And + Crd), which could be explained in different ways: (1) differences in the primary schist compositions, (2)
increased time-span of xenoliths in contact with the melt and (3) xenolith incorporation at slightly higher depths during
final granite ascent. 相似文献