Modelling of Subsidence at the Wairakei Geothermal Field, New Zealand
|
|
- Shavonne Matthews
- 5 years ago
- Views:
Transcription
1 Modelling of Subsidence at the Wairakei Geothermal Field, New Zealand Wesley Koros 1, John O Sullivan 1, Justin Pogacnik 1 and Mike O Sullivan 1 1 Department of Engineering Science, University of Auckland, Private Bag 92019, Auckland 1142, New Zealand w.koros@auckland.ac.nz Keywords: ABAQUS, TOUGH2, Subsidence, Elastoplastic, Wairakei geothermal field. ABSTRACT Significant localised ground subsidence related to geothermal production has been experienced within Wairakei geothermal field since the 1950s. Much of current knowledge of the cause of this subsidence anomaly has come from data gathered recently during an intensive subsidence investigation program. This anomaly is understood to be caused by both a reduction in pore pressure and hydrothermally altered weak formations. In order to quantitatively study the subsidence and predict its future behaviour, a subsidence modelling project is being carried out. A new radially symmetric subsidence model has been created combining geomechanics modelling using the finite element ABAQUS code and mass and heat flow modelling using TOUGH2. The coupled modelling approach calculates the local fluid pressure and temperature and then the related stress state using a nonlinear elasto-plastic constitutive model. The model parameters were determined using field data, laboratory measurements and calibration. The results show that the model is capable of accurately representing the subsidence that has occurred at Wairakei and that it will be a valuable tool for predicting rates of subsidence as a result of future production. 1. INTRODUCTION 1.1 Background Subsidence related to geothermal fluid production has occurred at Wairakei geothermal field in New Zealand, as shown in Figure 1. Damage related to subsidence can result in significant economic loss. At Wairakei field, for instance, subsidence has caused casing damage, tension cracks, ponding of the Wairakei stream, all requiring remedial action. This phenomenon attracted attention as early as the 1970s, after its discovery through re-levelling surveys conducted since the 1950s. Continuous monitoring of existing bench marks (BMs) has involved the use of variety of techniques that include First- Order levelling, a global positioning system (GPS) network (Energy Surveys, 2009) and, recently, Interferometric Synthetic Aperture Radar (InSAR) (Hole et al., 2007). This monitoring has produced important information regarding the magnitude and the rates of subsidence. The subsidence anomaly at Wairakei has evolved with time, as documented in publications by Hatton (1970), Allis (2000), Allis et al. (2009) and Bromley et al. (2015). Generally, subsidence results from compaction of low permeability, compressible layers when their effective stress is gradually increased as a result of pore pressure reduction. An increase in effective stress can be responsible for measurable deformation of reservoir material and this deformation is transmitted to the surface as subsidence. Figure 1: Subsidence bowls, monitoring wells (red) and bench marks (green) at Wairakei-Tauhara fields (modified after Bromley et al., 2010). Figure 2: A map showing Wairakei subsidence bowl, key bench marks (BMs) and wells November 2016 Auckland, New Zealand.
2 Subsidence at Wairakei has complex origins, arising from varying material properties and heterogeneous stratigraphy within the geothermal system. It is understood that ground subsidence at Wairakei is the response of local hydrothermally altered and compressible soft, porous materials to the changes in pressure within the geothermal system. Subsidence models are necessary to help understand the phenomenon. Previously one-dimensional models (Allis and Zhan, 2000; Allis, 2001; Koros et al. 2014, 2015) have been applied to model ground subsidence at Wairakei field. Generally, predictions made with one-dimensional models are limited. For example, they cannot predict the shape and extent of a subsidence bowl and hence the need for 2-D and 3-D multi-dimensional models. The present model is radially symmetric, utilizing r-z coordinates and thus is an approximation to the 3-D situation. Some horizontal detail can be represented in radial symmetry (r-z) and the vertical dimension is fully included. The modelling approach used in this work is based on Biot s (1941) poroelastic theory and Terzaghi s 1923 concept of effective stress (Terzaghi, 1943). Biot s poro-elastic theory can be applied to purely elastic and reversible deformation only. However, some reservoir material at Wairakei field contains yielding components that produce nonlinear elastoplastic behaviour. Therefore the concept of effective stress can be applied more broadly to elasto-plastic solid constitutive laws. 1.2 Purpose and scope The goal of this study is to predict subsidence at Wairakei subsidence bowl with a radial symmetric model (r-z) (shown in Figure 2) using a realistic geological setting and accurate representation of both fluid flow and rock mechanics. We use a loose coupling approach to solve a separate mass and energy transport model from the rock mechanics model. The effects of the model resolution were also investigated by comparing results for the pressure evolution and subsidence from both fine and coarse grids. The reservoir, mass and heat transfer, model was developed to reproduce the observed pressure response to production over the last 50 years. It covers a circle of radius 2970 m and extends from the surface down to a depth of 415 m. The model area includes a number of monitoring wells within and near the subsidence bowl. The reservoir model is implemented using AUTOUGH2 (Yeh et al., 2012) the University of Auckland s version of TOUGH2 (Pruess, 1999). Calibration of the model was carried out by adjusting permeability (k z and k r ), bottom mass flow and the model stratigraphy. This process was aimed at getting the model to match the observed pressure data. For the Wairakei subsidence bowl, sufficient data are available to calibrate the coupled THM (thermo-hydromechanical) radial model discussed here. With appropriate constitutive relationship, the calibrated model was able to capture the observed behavior and could be used to forecast potential future subsidence. The simulated results over a period of 50 years of production are presented and compared with the observed subsidence data. 1.3 Subsidence Mechanisms and History Subsidence in the Wairakei bowl is attributed mainly to a pressure decline in the low pressure shallow steam zone and consequent compaction of intensely altered, highly porous and compressible sediments and tuff materials within the Huka Falls unit. This is discussed in greater detail in the report by Bromley et al. (2010). The investigators agreed that the subsurface decline of pore pressure as result of geothermal production and hydrothermal alteration of weak formations are the main causes of subsidence at Wairakei. These pressure declines associated with geothermal production resulted in an increase in effective stress. Comparison of early levelling surveys of bench marks (BMs) conducted between the 1950s and 1970s at Wairakei field revealed the early development of subsidence. These surveys showed that the greatest subsidence had occurred at the eastern end of the Wairakei production field (Hatton, 1970; Allis, 1982). More re-levelling surveys of existing and newly established bench marks (survey networks), including the BMs shown in Figure 1, have shown continued ground subsidence, with a bowl shaped pattern, at BMs AA13, P128 and A97, as presented in Figure 3. All BMs start from a zero which in most cases was set in the 1950s. The surveys, which were limited in temporal and spatial details, indicated that subsidence continued slowly into the early 1960s. After which, the subsidence magnitude began increasing in the late 1960s and 1970s, as shown in Figure 3. These rapid rates of subsidence were followed by gradually decreasing rates. The subsidence rate at Wairakei subsidence anomaly at BM P128, (near monitoring well WKM 15) had declined to 45mm/yr in the re-levelling surveys and to 70 mm/yr ( ), down from 240mm/yr in the survey. At BM AA13 located at the edge of the bowl near monitoring well WKM 14, the subsidence rates have also declined to 34mm/yr (Energy Surveys, 2009). The bench marks and GPS (global positioning systems) surveys and InSAR analysis have delineated three other localized subsidence bowls within the Tauhara sector of the Wairakei-Tauhara geothermal field (see Figure 1). The geomechanics model was created using ABAQUS (ABAQUS, 2002). One of the key aspects of the rockmechanics model was the use of a nonlinear, elasto-plastic constitutive model for reservoir rock behaviour.the inclusion of plastic deformation has proven to be essential to accurately capture subsidence timing in Wairakei (Koros et al., 2014).The aim of calibration of the ABAQUS model was to match the subsidence data. 2
3 Figure 3: Subsidence trends at benchmarks (BMs) within and outside the Wairakei subsidence bowl from BMs positions are shown in Figure Stratigraphy and rock material properties The stratigraphic sequence observed in WKM15 and WKM14 is a series of sedimentary and volcanic deposits. Geophysical and geological surveys have shown Wairakei field to consist of water-laid sediments and volcanics. They comprise inter-bedded mudstones, siltstones, sandstones and tephra lavas. The volcanic units dominate deeper part of the geothermal system while sediments are more common near the surface. Recently, new volcanic and sediment formations have been identified, key stratigraphic units established and their depths reviewed (Rosenberg et al., 2009; Bignall et al., 2010). The recent exploratory drilling exercise provided knowledge and insights on geological controls on permeability in the reservoir. Wairakei Ignimbrite is covered by the Waiora formation. The Waiora Formation is made up of pyroclastic,volcaniclastic and lava deposits. Overlying the Waiora Formation is the Huka Falls Formation (HFF), which is made up of lacustrine sediments and water deposited tuffs (Grindley, 1965; Rosenberg et al., 2009). The HFF is subdivided into sub-units and is covered by superficial deposits up to the surface. The drilling in the subsidence bowl at WKM15 and outside the bowl at WKM14 revealed varying geology with different characteristics. Geotechnical investigations of compressibility of the cores from these wells provide the parameters which assisted the subsidence modelling discussed here. The impracticality of obtaining compressibility information in the field required the use of laboratory tests on sample cores to provide rock properties. The selected cores were subjected to K o tests whose results were reported in Pender, (2009) and Bromley et al. (2010).The cores from monitoring wells within and outside the subsidence bowl are of varying characteristics between and within different strata. In particular, cores from WKM15 in the upper Waiora Formation (depths: m) were relatively weak. Huka Falls Formation units (depths: m) were weak, with a low stiffness as observed from Point Load Index results and stiffness values from K 0 tests, presented in Figure 10.9 in Bromley, et al. (2010). At the same depth ranges in WKM14, the formations are moderately strong and have a larger stiffness. The variation in strength and stiffness of these formations (from weak to strong and low to high stiffness) is sensitive to the mineralogical composition of each formation. For example, clay mineralogy (smectite) properties (e.g. compressibility) present in the Huka Falls units are affected by thermal and hydrological changes. Strength of materials within the formations is dependent on variables such as clay content and water content. For instance, a material with high smectite content, as evident in WKM15, is more compressible than the rock with a low smectite content observed at WKM14. An increase in water content increases deformability. This is evidenced in Figure 10.9 in Bromley et al. (2010), which shows an increase in moisture content results in a reduction of uniaxial compressive strength and stiffness. 1.5 Approach to subsidence modeling We solve the subsidence problem by coupling a mass and heat transfer problem and a rock mechanics problem. The loosely coupled approach, applied to study subsidence at Wairakei, consists of the following steps: An r-z reservoir model (TOUGH2) is run to determine pore pressure and temperature variation within the reservoir system in response to production. An r-z rock mechanics model (ABAQUS) uses the pore pressure decline and temperature provided by the reservoir model at a sequence of times to calculate stresses and strains in the rock matrix, and subsidence at the surface. The reservoir model (TOUGH2) and geomechanics model (ABAQUS) were linked through a Python-based interface. A detailed explanation of this process is given in Pogacnik, et al. (2015). The ABAQUS subsidence model implemented in this work is based on Biot s (1941) poro-elastic theory. Geothermal fluid production cause significant pore pressure decline within the reservoir system. As a consequence, the effective stress increases and changes the state of stress and deformation in the rock matrix, possibly leading to subsidence at the ground surface. According to both Biot (1941) and Terzaghi (1943), effective stress is responsible for deformation of reservoir rock matrix material and is expressed as: Here, and are total stress tensor, pore fluid pressure and Kronecker delta respectively. The parameter is Biot s coefficient which controls the degree of pore pressure-stress coupling. It takes values between 0 and 1. A realistic subsidence prediction model requires development of a model that takes into account compressibility variation with increase in effective stress. Here we use a nonlinear, elasto-plastic model to simulate subsidence due to geothermal fluid production. Appropriate material properties were obtained from a large dataset of K o tests performed on core samples from a monitoring well within the Wairakei bowl (Pender, 2009). 2. COMPUTER MODEL DESCRIPTION 2.1 Geometry An axisymmetric model was used, hence implying horizontal stratigraphy and radial symmetry of material properties. The model is centred at WKM15 and extends radially for 2970 m towards WKM14 located outside the bowl. The simulations performed use both coarse and fine grids. Our fine grid has 83 x 5 m blocks (415 m) and the coarse grid has 16 x 25 m blocks (400 m). The model grids are shown in Figure 4. The left edge of the model coincides with the central axis of WKM15 (within the Wairakei bowl). The horizontal grid spacing is 25 m at the center of the subsidence bowl and increases logarithmically outward to 2970m. (1) 3
4 The geology of the model domain (shown in Fig. 5) consists of several formations down to ~400 m depth according to the lithostratigraphy provided in Rosenberg et al. (2009). elastic and non-elastic parameters obtained both from in situ measurements, laboratory tests and our previous 1-D subsidence modelling work. Figure 4: (a) Fine grid layout (b) Coarse grid layout. 2.2 TOUGH2 reservoir model No-flow boundary conditions are imposed on the reservoir model at the left (centre of the subsidence bowl) and right (edge of the bowl) sides of the grid. The top (ground surface) boundary is maintained at a constant temperature of 15 C and pressure of 0.1 MPa to represent atmospheric conditions. The base of the reservoir system is at 415 m (400 m for coarse model) where a thick low-permeability Waiora unit begins. Initial temperature at the bottom of the model is 186 C consistent with typical bottom temperature at Eastern borefield wells (Bixley, 2009). At the base of the model mass is extracted to represent the effects of deep production. Part of the model calibration task is to adjust the production rate at the base of the model so that the deep pressure matches the field data. All simulations were done with TOUGH2, with a water-air fluid property module (EOS3). 2.3 ABAQUS rock mechanics model Boundary conditions for the ABAQUS model are: no displacement on the right and bottom boundaries, vertical displacement only on the left boundary and traction-free conditions on the top boundary. Some of the reservoir material at Wairakei shows elasto-plastic behavior and a corresponding constitutive model was used, based on the Model Calibration The in situ geological data along with pore fluid pressure measurements are available. These allowed for a realistic evaluation of geomechanical properties of the formations included in the model, which is essential for a reliable prediction of subsidence. The horizontal stratigraphy was manually adjusted to fit the correct lithostratigraphy, as shown in Figure 5. The fine model with populated geology produced better results compared to coarse model because of refined geology that represented entire geology within the bowl. Permeability values and the mass flows at the base of the model were adjusted to match the pressure decline at three elevations. Manual calibration and the PEST software were used for parameter estimation. PEST calibration software (Doherty, 2004) was integrated with TOUGH2 and ABAQUS simulators to carry out a joint inversion of both the reservoir and rock mechanics models together, after a reasonably good TOUGH2 model had been obtained. This automated calibration process involved searching for a set of model parameters that gives the lowest possible value of objective function. PEST was used to identify permeability and stiffness values of defined rocktypes that give the closest match of the model results to the observed data. Stiffness values calibrated in this study fell within the range of values measured on core samples. Good calibration of the reservoir model is required obtaining a good match to the pressure behavior in the reservoir and consequently obtaining a good match to the subsidence. Reservoir model calibration was achieved by matching the model output to the measured pressure profiles at a few representative points. The pressure observations started in late 1950s and have continued to the present time. The model was calibrated over 50 years ( ). Mass drawn from the bottom of the model and permeabilities were adjusted to match historic pressure measurements (see Figure 6). These were achieved through application of pytough library (Croucher, 2011; Wellmann et al., 2012). The pore-water pressure and temperature results from the reservoir model at five-yearly intervals were used as input for stress-deformation analysis carried out with ABAQUS, and hence the calculation of subsidence. A loosely coupled model was used, i.e., the effect of stress changes on permeability, porosity and hence fluid flow was not considered. The rock-type structures of the two models are shown in Figure 5. As expected some of the detail of the fine model is lost in the coarse model. Decisions were made about the values of material parameters applicable to Wairakei geothermal field. For example, the magnitude of Biot s coefficient depends on the ratio of grain and bulk compressibility and this coefficient has not been determined for Wairakei rocks. However, laboratory measurements have been conducted on other sedimentary core samples to determine its value (Fabre and Gostkiewicz, 1997; Hou et al., 2005; Braun, 2009). Their observations revealed the magnitude of for clay-like material and high porosity sedimentary rocks ranges between Therefore, Biot-Willis coefficient was set to 0.5, hopefully typical of soft material at Wairakei.
5 Elasto-plastic behaviour for the rock matrix was applied, assuming von Mises yield condition applies and isotropic hardening was applied in predicting where yielding occurs. the subdivided Huka Falls, Taupo and Oruanui Formations sedimentary units used in the fine model had to be amalgamated in some cases. It was expected that the use of the fine rather than the coarse grid might have a significant impact on the mass and heat flow and hence subsidence. However the differences between the fine and coarse grids results are surprisingly small. Simulated results of pressure decline at three elevations are shown in Figure 6. These good results were obtained after calibrating the mass withdrawn from the base and the permeabilities. Pressure trends show a rapid decline after 1960 but a slow decline by the 1980s. A lag in the onset of significant pressure decline causes a delay at the start of subsidence. Changes in pressure trends cause the most significant changes in the subsidence values. These pressuretrends have been applied in modelling subsidence at BMs AA13, A97 (outside the bowl) and BM P128 (center of the bowl). The results of mass and heat transfer model are presented in the form of pressure profiles. The pressure profiles generated due to lowering of water-table and decreases in pore-pressure over 50 years of production are shown in Figure 6. These results were imported into ABAQUS for stress-deformation analysis and to predict subsidence. The plots show that the fine model results are better than the coarse model results, i.e. they fit the data better, but they are both quite similar. A limited amount of time was spent on calibrating the coarse model and it probably could be further improved. Figure 5: Grid layout with calibrated geology for (a) Fine model and (b) Coarse model. A more detailed discussion of this constitutive relation is given in Chaboche (2008) and Khan & Huang (1995). In addition, for yielding materials, a critical state parameter, initial and yield stresses and the corresponding plastic strains were defined. For detailed discussion on these parameters see Pogacnik, et al. (2015). 3. DISCUSSION OF RESULTS 3.1 General discussion The modelling results presented in this study are derived from a loose coupling of TOUGH2 and ABAQUS models and involved use of both fine and coarse grids. Good simulation results were obtained, both at the centre of the bowl and some distance away from the centre. 3.2 TOUGH2 Reservoir Model At Wairakei bowl, most of the subsidence occurs within Upper Huka Falls and Upper Waiora Formations. These formations are at elevations , m.a.s.l respectively. The reservoir materials are subdivided into sedimentary units that can be accurately matched to the geological model in the fine grid model. In the coarse grid, 5
6 Figure 6: Pressure trends from (a) fine model and (b) Coarse model at depths 150 m and 250 m (elevations:117 masl and 202 masl respectively) at the Wairakei bowl. Symbols are well data and lines (red, black, green, magenta) are simulated profiles from depth 250 m towards the surface (elevations: 290 masl and 350 masl respectively). Plots of pressure vs. depth, at five-yearly intervals are shown in Fig. 7. The pressures near the key zones for subsidence ( , m.a.s.l) are similar. 3.2 ABAQUS Rock Mechanics Model The geometry for stress-deformation model analysis in ABAQUS was imported from the both fine and coarse reservoir models shown in Figure 4. Material properties used in stress-deformation analysis were determined from a combination of laboratory tests results and in situ values published in geotechnical reports (Bromley et al, 2010; Pender, 2009). The simulation of the loosely-coupled TOUGH2-ABAQUS model was conducted for the period The results of the calculated subsidence, together with observed values at BM AA13 (outside the bowl near WKM 14), A97 (outside the bowl near WK 10) and BM P128 (at the centre of the bowl near WKM 15) for the period 1950s to 1997 are presented in Figure 8. The total observed subsidence from 1950 to 1997 is about 2.07 m and m at BM AA13 and P128 respectively. The vertical deformation in the models matches the data well. The largest subsidence of about m occurred at the centre of the bowl, with decreasing subsidence modelled with increasing radial distance from the centre. Simulated subsidence at 962 m and 1300 m was about 4.5 m and 2.1 m respectively in both coarse and fine models. A comparison between the observed subsidence and the results from the simulations in both model grids shows a good match. The match is somewhat better for the fine model but further calibration could improve the coarse model results. Detailed geology was only available at WKM 15 at the inner boundary of the model and WKM 14 near the outer boundary of the model. The geology in between was assigned by interpolation with some adjustments made during the calibration process. Figure 7: Simulated pressure-elevation trends for (a) Fine model and (b) Coarse model, at 5-yearly time increments between 1950 and
7 Allis, R., Bromley, C. Currie, S.: Update on subsidence at the Wairakei-Tauhara geothermal system, New Zealand. Geothermics, 38, (2009). Allis, R.G.: Review of subsidence at Wairakei field, New Zealand. Geothermics, 29, (2000). Allis, R. G.: Comparison of Subsidence at Wairakei, Broadlands and Kawerau Fields, New Zealand. Proc. 8 th Workshop Geothermal Reservoir Engineering, Stanford University, California (1982). Biot, M.A.: General theory of three-dimensional consolidation. Journal of Applied Physics,12(2), (1941). Figure 8: Subsidence profiles for (a) Fine model and (b) Coarse model for 50 years ( ) of production. Horizontal distance is given as x: at x= 0.0 m (BM P128, near WKM15; at x=962.0 m (BM A97, near WK10, outside the bowl) and at x= m (BM AA13, near WKM14, outside the bowl). 4. CONCLUSION Numerical simulations were conducted, with a radially symmetric model, with loosely coupled mass and heat transfer and rock mechanics to model subsidence at Wairakei geothermal field. Both coarse and fine models were applied. Simulations with the calibrated models gave acceptable results for the 50 years of production. Simulation results revealed the following: Calibration of key parameters, particularly permeability, the mass withdrawal, stiffness and yield parameters led to achieving a good match to the surface displacement. The subsidence model was progressively developed by first calibrating the TOUGH2 reservoir model to match the observed pressure decline and then calibrating the ABAQUS rock mechanics model to match subsidence. The magnitude of subsidence is very high at the centre of the bowl compared that further out. Attention was also focused on constitutive behaviour of the rocks. An elasto-plastic constitutive model was required to effectively represent the behaviour of some of the soft rocks. The application of an axisymmetric model, with fine and coarse grids, to modelling subsidence at Wairakei has shown that it is possible to make good models of subsidence if sufficient data are available. ACKNOWLEDGEMENTS The authors gratefully acknowledge the support by Contact Energy Ltd and University of Auckland. REFERENCES ABAQUS User s Manual: ABAQUS/Standard User s Manual, Volume II, Version 6.3. Hibbit & Sorenson, Inc., Pawtucket, Rhode Island. (2002). Bignall, G.,Milicich, S.D., Ramirez, L.E., Rosenberg, M.D., Kilgour, G.N.,Rae, A.J.: Geology of the Wairakei- Tauhara Geothermal System, New Zealand. Proc. Worlds Geothermal Congress 2010, Bali, Indonesia. 8 pp. (2010). Braun E.: Determination with RACOS of the 3D in situ loadings, the pore pressure effectivenesses and the deformation parameters for several strata in the Salzwedel reservoir area. Report submitted to TU Clausthal, 64 pp. (2009). Bromley, C.J., Currie, S., Jolly, S., Mannington, W.: Subsidence: an Update on New Zealand Geothermal Deformation Observations and Mechanisms. Proc. World Geothermal Congress 2015, Melbourne, Australia. 9 pp. (2015). Bromley, C.J., Currie, S., Ramsay, G., Rosenberg, M., Pender, M., O Sullivan, M., Lynne, B., Lee, S.G., Brockbank, K., Glynn-Morris, T., Mannington, W., Garvey, J.: Tauhara Stage II Geothermal Project: Subsidence Report, GNS Science. (2010). Chaboche, J.L.: A review of some plasticity and viscoplasticity constitutive theories. International Journal of Plasticity, 24, (2008). Croucher, A.E.: PyTOUGH: a Python scripting library for automating TOUGH2 simulations, Proc. 33 rd New Zealand Geothermal Workshop, University of Auckland, Auckland, NZ. (2011). Doherty, J.: PEST, Model-Independent Parameter Estimation, User Manual: 5th Edition, Watermark Numerical Computing, Brisbane, Australia. (2004) Energy surveys,: Regional Precise levelling surveys for Wairakei &Tauhara. Report submitted to Contact Energy. (2009). Fabre, D., Gustkiewicz, J.: Poroelastic properties of limestones and sandstones under hydrostatic conditions. Int. J. Rock Mech. Min. Sci. 34 (1), (1997). Grindley, G.W.: The geology, structure and exploitation of the Wairakei field, Taupo, New Zealand. NZ Geological Survey Bulletin, 75, 131 pp. (1965). Hatton, J. W.: Ground Subsidence of a Geothermal Field during exploitation. Geothermics, 2, (1970). 7
8 Hole, J. K., Bromley, C. J., Stevens, N. F., Wadge, G.: Subsidence in the geothermal fields of the Taupo Volcanic Zone, New Zealand from 1996 to 2005 measured by InSAR. Journal of Volcanology and Geothermal Research, 166, (2007). Hou, Z., Somerville, J., Hutcheon, R.: Experimental investigation of the hydro-mechanical behavior of tight gas formation sandstone and its impacts on the wellbore stability. ISRM-Eurock, Brno, Czech Republic (2005). Khan, A.S., Huang, S.: Continuum Theory of Plasticity. John Wiley & Sons Inc., New York. (2005). Koros, W., O Sullivan, M., Pogacnik, J., O Sullivan, J.: Coupled One-Dimensional Subsidence Modelling Using Elastic-Plastic Rheology. Proc. 36 th New Zealand Geothermal Workshop, Auckland, New Zealand. 9 pp. (2014). Koros, W., O Sullivan, M., Pogacnik, J., O Sullivan, J.: Comparison of Coupled One-Dimensional Subsidence Models. Proc. World Geothermal Congress 2015, Melbourne, Australia. 10 pp. (2015). Pender, M.: Appendix 9: Triaxial testing on core from the Wairakei-Tauhara geothermal field. Interpretive Procedures Report,University of Auckland, Auckland. (2009). Pogacnik, J., Koros, W., O Sullivan, J., O Sullivan, M.: Using ABAQUS to Simulate the Triaxial Stress Strain Behavior of Core Samples from the Wairakei/Tauhara Geothermal Field. Proc. 37 th New Zealand Geothermal Workshop, Taupo, New Zealand. 11 pp. (2015). Pruess, K., Oldenburg, C., Moridis, G.: TOUGH2 User s Guide Version 2.0, University of California, Berkeley, California, USA. (1999). Rosenberg, M.D., Bignall, G., Rae, A.J.: The geological framework of the Wairakei-Tauhara Geothermal System, New Zealand. Geothermics, 38, (2009). Terzaghi, K.: Theoretical Soil Mechanics, Wiley, New York, N. Y., 510 pp. (1943). Wellmann, J.F., Croucher, A.E., and Regenauer-Lieb, K.: Python scripting libraries for subsurface fluid and heat flow simulations with TOUGH2 and SHEMAT. Computers & Geosciences, 43, (2012). Yeh, A., Croucher, A. E., and O Sullivan, M. J.: Recent developments in the AUTOUGH2 simulator. Proceedings of the TOUGH Symposium Lawrence Berkeley National Laboratory, Berkeley, California. 8pp (2012). 8
OVERVIEW OF THE WAIRAKEI-TAUHARA SUBSIDENCE INVESTIGATION PROGRAM
PROCEEDINGS, Thirty-Sixth Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, January 31 - February 2, 2011 SGP-TR-191 OVERVIEW OF THE WAIRAKEI-TAUHARA SUBSIDENCE INVESTIGATION
More informationComparison of Coupled One-Dimensional Subsidence Models
Proceedings World Geothermal Congress 2015 Melbourne, Australia, 19-25 April 2015 Comparison of Coupled One-Dimensional Subsidence Models Wesley Koros 1, John O Sullivan 1, Justin Pogacnik 1, Mike O Sullivan
More informationCHARACTERIZING FEED ZONES IN GEOTHERMAL FIELDS: INTEGRATED LEARNINGS FROM COMPLETION TESTING, IMAGE LOGS AND CONTINUOUS CORE
CHARACTERIZING FEED ZONES IN GEOTHERMAL FIELDS: INTEGRATED LEARNINGS FROM COMPLETION TESTING, IMAGE LOGS AND CONTINUOUS CORE Trystan Glynn-Morris, Katie Mclean and Kerin Brockbank Contact Energy Ltd, Wairakei
More informationNOTICE CONCERNING COPYRIGHT RESTRICTIONS
NOTICE CONCERNING COPYRIGHT RESTRICTIONS This document may contain copyrighted materials. These materials have been made available for use in research, teaching, and private study, but may not be used
More informationUsing Repeat Microgravity Measurements to Track Reinjection in Liquid-Dominated Fields
Proceedings World Geothermal Congress 25 Antalya, Turkey, 24-29 April 25 Using Repeat Microgravity Measurements to Track Reinjection in Liquid-Dominated Fields Trevor M. Hunt Institute of Geological &
More informationTrevor M. Hunta and Warwick M. %slingb
PROCEEDINGS, Nineteenth Workshop on Geothennal Reservoir Engineering Stanford University, Stanford, California, January 16-2, 1994 SGP-TR-147 Estimation of reservoir permeability using gravity change measurements
More informationGEOTECHNICAL INVESTIGATIONS OF SUBSIDENCE IN THE WAIRAKEI- TAUHARA GEOTHERMAL FIELD
GEOTECHNICAL INVESTIGATIONS OF SUBSIDENCE IN THE WAIRAKEI- TAUHARA GEOTHERMAL FIELD Graham Ramsay 1, Trystan Glynn-Morris 2, Michael Pender 3 and Melvyn Griffiths 4 1 Beca Infrastructure Limited P O Box
More informationRESPONSE OF WAIRAKEI GEOTHERMAL RESERVOIR TO 40 YEARS OF PRODUCTION
RESPONSE OF WAIRAKEI GEOTHERMAL RESERVOIR TO 4 YEARS OF PRODUCTION Allan Clotworthy Contact Energy Ltd, Private Bag 21, Taupo, New Zealand Key Words: geothermal, production, reinjection, reservoir, Wairakei
More informationNumerical Simulation of Devolution and Evolution of Steam-Water Two-Phase Zone in a Fractured Geothermal Reservoir at Ogiri, Japan
GRC Transactions, Vol. 37, 2013 Numerical Simulation of Devolution and Evolution of Steam-Water Two-Phase Zone in a Fractured Geothermal Reservoir at Ogiri, Japan Yohei Tateishi 1, Ryuichi Itoi 1, Toshiaki
More informationAdvances in Subsidence Modelling of Exploited Geothermal Fields
Proceedings World Geothermal Congress 2005 Antalya, Turkey, 24-29 April 2005 Advances in Subsidence Modelling of Exploited Geothermal Fields Philip J. White, James V. Lawless, Sergei Terzaghi, Wataru Okada
More informationCOMPUTER MODELLING OF HEAT AND MASS FLOW IN STEAMING GROUND AT KARAPITI THERMAL AREA, NEW ZEALAND
PROCEEDINGS, Twenty-Ninth Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, January 26-28, 2004 SGP-TR-175 COMPUTER MODELLING OF HEAT AND MASS FLOW IN STEAMING GROUND
More informationCase Study: Tauhara New Zealand. Santiago de Chile, May 2014
Case Study: Tauhara New Zealand The Geothermal Institute University of Auckland Bridget Lynne Santiago de Chile, 26-29 May 2014 Bridget Y. Lynne Mick Pender Trystan Glynn-Morris Tauhara Subsidence Case
More informationTu P8 08 Modified Anisotropic Walton Model for Consolidated Siliciclastic Rocks: Case Study of Velocity Anisotropy Modelling in a Barents Sea Well
Tu P8 08 Modified Anisotropic Walton Model for Consolidated Siliciclastic Rocks: Case Study of Velocity Anisotropy Modelling in a Barents Sea Well Y. Zhou (Rock Solid Images), F. Ruiz (Repsol), M. Ellis*
More informationMethods of Interpreting Ground Stress Based on Underground Stress Measurements and Numerical Modelling
University of Wollongong Research Online Coal Operators' Conference Faculty of Engineering and Information Sciences 2006 Methods of Interpreting Ground Stress Based on Underground Stress Measurements and
More informationA Constitutive Framework for the Numerical Analysis of Organic Soils and Directionally Dependent Materials
Dublin, October 2010 A Constitutive Framework for the Numerical Analysis of Organic Soils and Directionally Dependent Materials FracMan Technology Group Dr Mark Cottrell Presentation Outline Some Physical
More informationAPPLICATION OF 1D HYDROMECHANICAL COUPLING IN TOUGH2 TO A DEEP GEOLOGICAL REPOSITORY GLACIATION SCENARIO
PROCEEDINGS, TOUGH Symposium 2015 Lawrence Berkeley National Laboratory, Berkeley, California, September 28-30, 2015 APPLICATION OF 1D HYDROMECHANICAL COUPLING IN TOUGH2 TO A DEEP GEOLOGICAL REPOSITORY
More informationIN-SITU STRESS ESTIMATION IN OFFSHORE EASTERN MEDITERRANEAN WITH FINITE ELEMENT ANALYSIS
8 th GRACM International Congress on Computational Mechanics Volos, 12 July 15 July 2015 IN-SITU STRESS ESTIMATION IN OFFSHORE EASTERN MEDITERRANEAN WITH FINITE ELEMENT ANALYSIS Anna Kyriacou 1, Panos
More informationLand subsidence due to groundwater withdrawal in Hanoi, Vietnam
Land Subsidence (Proceedings of the Fifth International Symposium on Land Subsidence, The Hague, October 1995). 1AHS Publ. no. 234, 1995. 55 Land subsidence due to groundwater withdrawal in Hanoi, Vietnam
More informationSHEAR-SLIP ANALYSIS IN MULTIPHASE FLUID-FLOW RESERVOIR ENGINEERING APPLICATIONS USING TOUGH-FLAC
PROCEEDINGS, TOUGH Symposium 6 Lawrence Berkeley National Laboratory, Berkeley, California, May 5 7, 6 SHEAR-SLIP ANALYSIS IN MULTIPHASE FLUID-FLOW RESERVOIR ENGINEERING APPLICATIONS USING TOUGH-FLAC Jonny
More informationEvaluation of Stress and Strain Induced by the Rock Compaction on a Hydrocarbon Well Completion Using Contact Interfaces with Abaqus
Evaluation of Stress and Strain Induced by the Rock Compaction on a Hydrocarbon Well Completion Using Contact Interfaces with Abaqus G. Capasso and G. Musso * eni E&P, S. Donato Milanese, MI, Italy * now
More informationA Geothermal Reservoir Simulator with AD-GPRS
Proceedings World Geothermal Congress 2015 Melbourne, Australia, 19-25 April 2015 A Geothermal Reservoir Simulator with Zhi Yang Wong, Roland Horne and Denis Voskov Department of Energy Resources Engineering,
More informationInfluences of Geological Depositional Processes and Environments on Geothermal Strata from Wairakei-Tauhara Geothermal Field, New Zealand
Proceedings World Geothermal Congress 2015 Melbourne, Australia, 19-25 April 2015 Influences of Geological Depositional Processes and Environments on Geothermal Strata from Wairakei-Tauhara Geothermal
More informationWellbore stability analysis in porous carbonate rocks using cap models
Wellbore stability analysis in porous carbonate rocks using cap models L. C. Coelho 1, A. C. Soares 2, N. F. F. Ebecken 1, J. L. D. Alves 1 & L. Landau 1 1 COPPE/Federal University of Rio de Janeiro, Brazil
More informationStudy on Numerical Simulation of Steam Huff and Puff Based on Deformable Medium Model
SCIREA Journal of Mine Engineering http://www.scirea.org/journal/mine October 1, 2016 Volume 1, Issue1, October 2016 Study on Numerical Simulation of Steam Huff and Puff Based on Deformable Medium Model
More informationCritical Borehole Orientations Rock Mechanics Aspects
Critical Borehole Orientations Rock Mechanics Aspects By R. BRAUN* Abstract This article discusses rock mechanics aspects of the relationship between borehole stability and borehole orientation. Two kinds
More informationTIME-DEPENDENT BEHAVIOR OF PILE UNDER LATERAL LOAD USING THE BOUNDING SURFACE MODEL
TIME-DEPENDENT BEHAVIOR OF PILE UNDER LATERAL LOAD USING THE BOUNDING SURFACE MODEL Qassun S. Mohammed Shafiqu and Maarib M. Ahmed Al-Sammaraey Department of Civil Engineering, Nahrain University, Iraq
More informationEarth Research 3-Dimensional Geological Modelling of Geothermal Systems in New Zealand a New Visualisation Tool
Proceedings World Geothermal Congress 2010 Bali, Indonesia, 25-29 April 2010 Earth Research 3-Dimensional Geological Modelling of Geothermal Systems in New Zealand a New Visualisation Tool Sarah D. Milicich
More informationMODELING OF GAS MIGRATION THROUGH LOW-PERMEABILITY CLAY USING INFORMATION ON PRESSURE AND DEFORMATION FROM FAST AIR INJECTION TESTS
PROCEEDINGS, TOUGH Symposium 2015 Lawrence Berkeley National Laboratory, Berkeley, California, September 28-30, 2015 MODELING OF GAS MIGRATION THROUGH LOW-PERMEABILITY CLAY USING INFORMATION ON PRESSURE
More information4D stress sensitivity of dry rock frame moduli: constraints from geomechanical integration
Title 4D stress sensitivity of dry rock frame moduli: constraints from geomechanical integration Authors Bloomer, D., Ikon Science Asia Pacific Reynolds, S., Ikon Science Asia Pacific Pavlova, M., Origin
More informationEffect of Hydrothermal Alteration on Rock Properties in the Tauhara Geothermal Field, New Zealand
Proceedings World Geothermal Congress 2015 Melbourne, Australia, 19-25 April 2015 Effect of Hydrothermal Alteration on Rock Properties in the Tauhara Geothermal Field, New Zealand P. Mielke, A. Prieto,
More informationIntegrated Geophysical Model for Suswa Geothermal Prospect using Resistivity, Seismics and Gravity Survey Data in Kenya
Proceedings World Geothermal Congress 2015 Melbourne, Australia, 19-25 April 2015 Integrated Geophysical Model for Suswa Geothermal Prospect using Resistivity, Seismics and Gravity Survey Data in Kenya
More informationUnderstanding hydraulic fracture variability through a penny shaped crack model for pre-rupture faults
Penny shaped crack model for pre-rupture faults Understanding hydraulic fracture variability through a penny shaped crack model for pre-rupture faults David Cho, Gary F. Margrave, Shawn Maxwell and Mark
More informationSubsidence: an Update on New Zealand Geothermal Deformation Observations and Mechanisms
Proceedings World Geothermal Congress 2015 Melbourne, Australia, 19-25 April 2015 Subsidence: an Update on New Zealand Geothermal Deformation Observations and Mechanisms Chris J. Bromley 1, Steve Currie
More informationRock Physics Perturbational Modeling: Carbonate case study, an intracratonic basin Northwest/Saharan Africa
Rock Physics Perturbational Modeling: Carbonate case study, an intracratonic basin Northwest/Saharan Africa Franklin Ruiz, Carlos Cobos, Marcelo Benabentos, Beatriz Chacon, and Roberto Varade, Luis Gairifo,
More informationChapter 1 INTRODUCTION
CONSTITUTIVE MODELS AND CHALK Chapter 1 INTRODUCTION Computing power has advanced significantly in the years since soil mechanics and rock mechanics first became mature fields. Probably the single most
More informationA Closer Look At Hydrothermal Alteration and Fluid-Rock Interaction Using Scanning Electron Microscopy
Proceedings World Geothermal Congress 2015 Melbourne, Australia, 19-25 April 2015 A Closer Look At Hydrothermal Alteration and Fluid-Rock Interaction Using Scanning Electron Microscopy Bridget Y. Lynne
More informationTensor character of pore pressure/stress coupling in reservoir depletion and injection
Tensor character of pore pressure/stress coupling in reservoir depletion and injection Müller, B., Altmann, J.B., Müller, T.M., Weißhardt, A., Shapiro, S., Schilling, F.R., Heidbach, O. Geophysical Institute
More informationFRACTURE REORIENTATION IN HORIZONTAL WELL WITH MULTISTAGE HYDRAULIC FRACTURING
SPE Workshop OILFIELD GEOMECHANICS Slide 1 FRACTURE REORIENTATION IN HORIZONTAL WELL WITH MULTISTAGE HYDRAULIC FRACTURING A. Pimenov, R. Kanevskaya Ltd. BashNIPIneft March 27-28, 2017 Moscow, Russia Slide
More informationProduction-induced stress change in and above a reservoir pierced by two salt domes: A geomechanical model and its applications
Production-induced stress change in and above a reservoir pierced by two salt domes: A geomechanical model and its applications Peter Schutjens, Jeroen Snippe, Hassan Mahani, Jane Turner, Joel Ita and
More informationA 3D MULTI-DISCIPLINARY INTERPRETATION OF THE BASEMENT OF THE TAUPO VOLCANIC ZONE, NEW ZEALAND
PROCEEDINGS, Thirty-Seventh Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, January 30 - February 1, 2012 SGP-TR-194 A 3D MULTI-DISCIPLINARY INTERPRETATION OF THE
More informationAvailable online at ScienceDirect. Procedia Engineering 158 (2016 )
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 158 (2016 ) 344 349 VI ITALIAN ONFERENE OF RESEARHERS IN GEOTEHNIAL ENGINEERING Geotechnical Engineering in Multidisciplinary
More informationLiquefaction potential of Rotorua soils
Pearse-Danker, E. (2013) Liquefaction potential of Rotorua soils Proc. 19 th NZGS Geotechnical Symposium. Ed. CY Chin, Queenstown Liquefaction potential of Rotorua soils E Pearse-Danker Coffey Geotechnics
More informationNATURAL STATE AND FUTURE PRODUCTION MODELLING OF ARJUNO- WELIRANG GEOTHERMAL FIELD, INDONESIA
NATURAL STATE AND FUTURE PRODUCTION MODELLING OF ARJUNO- WELIRANG GEOTHERMAL FIELD, INDONESIA Ardha Wardana 1, John O Sullivan 1 and Michael O Sullivan 1 1 Geothermal Institute and Department of Engineering
More informationINTERPRETATION OF INTERFERENCE EFFECTS IN THREE PRODUCTION WELLS IN THE KAWERAU GEOTHERMAL FIELD, NEW ZEALAND. Lynell Stevens and Kevin J Koorey
PROCEEDINGS, Twenty-First Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, January 22-24, 1996 SGP-TR-151 INTERPRETATION OF INTERFERENCE EFFECTS IN THREE PRODUCTION
More informationand reservoir simulation. Figure 1. Interdependence of experimental, poromechanical, and field study variables. 1 GeoFluids 2020_v2.
UT GeoFluids 2020 Research Plan Peter B. Flemings1, John T. Germaine2, Maria A. Nikolinakou3 1. pflemings@jsg.utexas.edu, 2. john.germaine@tufts.edu. 3. mariakat@mail.utexas.edu Executive summary The UT
More informationNumerical Simulation Study of the Mori Geothermal Field, Japan
Proceedings World Geothermal Congress 2010 Bali, Indonesia, 25-29 April 2010 Numerical Simulation Study of the Mori Geothermal Field, Japan Kazuyoshi Osada 1, Mineyuki Hanano 2, Kei Sato 1, Tatsuya Kajiwara
More informationLandslide FE Stability Analysis
Landslide FE Stability Analysis L. Kellezi Dept. of Geotechnical Engineering, GEO-Danish Geotechnical Institute, Denmark S. Allkja Altea & Geostudio 2000, Albania P. B. Hansen Dept. of Geotechnical Engineering,
More informationIn situ fracturing mechanics stress measurements to improve underground quarry stability analyses
In situ fracturing mechanics stress measurements to improve underground quarry stability analyses Anna M. Ferrero, Maria R. Migliazza, Andrea Segalini University of Parma, Italy Gian P. Giani University
More informationCPT Data Interpretation Theory Manual
CPT Data Interpretation Theory Manual 2016 Rocscience Inc. Table of Contents 1 Introduction... 3 2 Soil Parameter Interpretation... 5 3 Soil Profiling... 11 3.1 Non-Normalized SBT Charts... 11 3.2 Normalized
More informationJohn E. Gale 1 and Eunjeong Seok 2
Field and Laboratory Coupled Fracture Deformation-Pore Pressure-Permeability Experiments That Provide Insight for Depressurization of Fractured Rock Slopes John E. Gale 1 and Eunjeong Seok 2 1 Fracflow
More informationJihoon Kim, George J. Moridis, John Edmiston, Evan S. Um, Ernest Majer. Earth Sciences Division, Lawrence Berkeley National Laboratory 24 Mar.
TOUGH+ROCMECH for the Analysis of coupled Flow, Thermal, Geomechanical and Geophysical Processes Code Description and Applications to Tight/Shale Gas Problems Jihoon Kim, George J. Moridis, John Edmiston,
More informationP16 Gravity Effects of Deformation Zones Induced by Tunnelling in Soft and Stiff Clays
P16 Gravity Effects of Deformation Zones Induced by Tunnelling in Soft and Stiff Clays V. Blecha* (Charles University) & D. Mašín (Charles University) SUMMARY We calculated gravity response of geotechnical
More informationNumerical and Theoretical Study of Plate Load Test to Define Coefficient of Subgrade Reaction
Journal of Geotechnical and Transportation Engineering Volume 1 Issue 2 Numerical and Theoretical Study of Plate Load Test to Define Coefficient of Subgrade Reaction Naeini and Taherabadi Received 9/28/2015
More informationLiquefaction-Induced Lateral Spreading Misko Cubrinovski University of Canterbury, Christchurch, New Zealand
US New Zealand Japan International Workshop Liquefaction-Induced Ground Movements Effects UC Berkeley, California, 2 4 November 2016 Liquefaction-Induced Lateral Spreading Misko Cubrinovski University
More informationNUMERICAL MODELING STUDY OF SIBAYAK GEOTHERMAL RESERVOIR, NORTH SUMATRA, INDONESIA
PROCEEDINGS, Twenty-Sixth Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, January 29-31, 21 SGP-TR-168 NUMERICAL MODELING STUDY OF SIBAYAK GEOTHERMAL RESERVOIR,
More informationTopic 7: Geophysical and Remote Sensing Models. Presenter: Greg Ussher Organisation:SKM
Topic 7: Geophysical and Remote Sensing Models Presenter: Greg Ussher Organisation:SKM Email: gussher@globalskm.com Overview Surface-based imaging of subsurface characteristics That indicate geothermal
More informationExact Solutions of Two-dimensional and Tri-dimensional Consolidation Equations
Exact Solutions of Two-dimensional and Tri-dimensional Consolidation Equations Romolo Di Francesco GEO&GEO Instruments! - research & development Teramo (TE), Italy E-mail: romolo.difrancesco@vodafone.it
More information6298 Stress induced azimuthally anisotropic reservoir - AVO modeling
6298 Stress induced azimuthally anisotropic reservoir - AVO modeling M. Brajanovski* (Curtin University of Technology), B. Gurevich (Curtin University of Technology), D. Nadri (CSIRO) & M. Urosevic (Curtin
More informationAnalyzing effect of fluid flow on surface subsidence
Analyzing effect of fluid flow on surface subsidence in mining area Y. Abousleiman", M. Bai\ H. Zhang', T. Liu" and J.-C. Roegiers* a. School of Engineering and Architecture, The Lebanese American University,
More informationRegional-Scale Salt Tectonics Modelling: Bench-Scale Validation and Extension to Field-Scale Predictions
Software Limited Bench-Scale to Field-Scale Regional-Scale Salt Tectonics Modelling: Bench-Scale Validation and Extension to Field-Scale Predictions Thornton, Dean A., Roberts, Daniel T., Crook, Anthony
More informationWELLBORE STABILITY ANALYSIS IN CHEMICALLY ACTIVE SHALE FORMATIONS
S911 WELLBORE SABILIY ANALYSIS IN CHEMICALLY ACIVE SHALE FORMAIONS by Xiang-Chao SHI *, Xu YANG, Ying-Feng MENG, and Gao LI State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation, Southwest
More informationModule 9 : Foundation on rocks. Content
FOUNDATION ON ROCKS Content 9.1 INTRODUCTION 9.2 FOUNDATION TYPES ON ROCKS 9.3 BEARING CAPCITY- SHALLOW FOUNDATION 9.3.1 Ultimate bearing capacity 9.3.2 Safe bearing pressure 9.3.3 Estimation of bearing
More informationThe evolution of the reservoir stress state throughout the history of production
International Journal of Petroleum and Geoscience Engineering (IJPGE) 2 (2): 181- ISSN 2289-4713 Academic Research Online Publisher Research paper The evolution of the reservoir stress state throughout
More informationStrength, creep and frictional properties of gas shale reservoir rocks
ARMA 1-463 Strength, creep and frictional properties of gas shale reservoir rocks Sone, H. and Zoback, M. D. Stanford University, Stanford, CA, USA Copyright 21 ARMA, American Rock Mechanics Association
More informationMODELING OF ELASTO-PLASTIC MATERIALS IN FINITE ELEMENT METHOD
MODELING OF ELASTO-PLASTIC MATERIALS IN FINITE ELEMENT METHOD Andrzej Skrzat, Rzeszow University of Technology, Powst. Warszawy 8, Rzeszow, Poland Abstract: User-defined material models which can be used
More informationTHE INFLUENCE OF ROOF BOLTS LOCATION ON ITS INTERACTION WITH THE ROCK MASS.
THE INFLUENCE OF ROOF BOLTS LOCATION ON ITS INTERACTION WITH THE ROCK MASS. M. Cała 1, A. Tajduś 1 ABSTRACT This paper examines the influence of roof bolts location on its interaction with rock mass in
More informationractical Geomechanics for Unconventional Resources
P ractical Geomechanics for Unconventional Resources 24-26 October 2012, Calgary, Canada Practical Geomechanics for Unconventional Resources Nowadays, unconventional resources have been brought into the
More informationCrosswell tomography imaging of the permeability structure within a sandstone oil field.
Crosswell tomography imaging of the permeability structure within a sandstone oil field. Tokuo Yamamoto (1), and Junichi Sakakibara (2) (1) University of Miami and Yamamoto Engineering Corporation, (2)
More information11/22/2010. Groundwater in Unconsolidated Deposits. Alluvial (fluvial) deposits. - consist of gravel, sand, silt and clay
Groundwater in Unconsolidated Deposits Alluvial (fluvial) deposits - consist of gravel, sand, silt and clay - laid down by physical processes in rivers and flood plains - major sources for water supplies
More information1. Background. is usually significantly lower than it is in uniaxial tension
NOTES ON QUANTIFYING MODES OF A SECOND- ORDER TENSOR. The mechanical behavior of rocks and rock-like materials (concrete, ceramics, etc.) strongly depends on the loading mode, defined by the values and
More informationRock Mechanics Laboratory Tests for Petroleum Applications. Rob Marsden Reservoir Geomechanics Advisor Gatwick
Rock Mechanics Laboratory Tests for Petroleum Applications Rob Marsden Reservoir Geomechanics Advisor Gatwick Summary A wide range of well established and proven laboratory tests are available for petroleum
More informationModeling pressure response into a fractured zone of Precambrian basement to understand deep induced-earthquake hypocenters from shallow injection
Modeling pressure response into a fractured zone of Precambrian basement to understand deep induced-earthquake hypocenters from shallow injection S. Raziperchikolaee 1 and J. F. Miller 1 Abstract Analysis
More informationBest Practice Reservoir Characterization for the Alberta Oil Sands
Best Practice Reservoir Characterization for the Alberta Oil Sands Jason A. McLennan and Clayton V. Deutsch Centre for Computational Geostatistics (CCG) Department of Civil and Environmental Engineering
More informationReservoir Simulator Compaction Modelling: A Predictor for Accelerated Coupled Rock Mechanics -- Reservoir Simulation
Reservoir Simulator Compaction Modelling: A Predictor for Accelerated Coupled Rock Mechanics -- Reservoir Simulation by Øystein Pettersen Centre for Integrated Petroleum Research, Bergen, Norway ECMOR
More informationMultiphysics modelling of the Mandel-Cryer effect
Int. Jnl. of Multiphysics Volume 10 Number 1 2016 11 Multiphysics modelling of the Mandel-Cryer effect E Holzbecher* German University of Technology in Oman (GUtech) Muscat, Oman ABSTRACT In porous medium
More informationSoil stiffness measured in oedometer tests
Wesley, L. & Pender, M. (2008) Proc. 18 th NZGS Geotechnical Symposium on Soil-Structure Interaction. Ed. CY Chin, Auckland s Laurence Wesley and Michael Pender Faculty of Engineering, University of Auckland,
More informationUSER S MANUAL 1D Seismic Site Response Analysis Example University of California: San Diego August 30, 2017
USER S MANUAL 1D Seismic Site Response Analysis Example http://www.soilquake.net/ucsdsoilmodels/ University of California: San Diego August 30, 2017 Table of Contents USER'S MANUAL TABLE OF CONTENTS Page
More informationNew Plymouth CBD Site Subsoil Class: Results from ground investigation
Alexander, G.J., Chin, C.Y., Kayser, C. & Bradshaw, J. (2017) New Plymouth CBD Site Subsoil Proc. 20 th NZGS Geotechnical Symposium. Eds. GJ Alexander & CY Chin, Napier New Plymouth CBD Site Subsoil Class:
More informationModelling the excavation damaged zone in Callovo-Oxfordian claystone with strain localisation
Modelling the excavation damaged zone in Callovo-Oxfordian claystone with strain localisation B. Pardoen - F. Collin - S. Levasseur - R. Charlier Université de Liège ArGEnCo ALERT Workshop 2012 Aussois,
More informationMODELLING OF CONDITIONS CLOSE TO GEOTHERMAL HEAT SOURCES
MODELLING OF CONDITIONS CLOSE TO GEOTHERMAL HEAT SOURCES Sæunn Halldórsdóttir (ÍSOR), Inga Berre (UiB), Eirik Keilegavlen (UiB) and Guðni Axelsson (ÍSOR) Photo: Auður Agla, ÍSOR Heat transfer in geothermal
More informationRock Cutting Analysis Employing Finite and Discrete Element Methods
Journal of Mechanics Engineering and Automation 6 (2016) 100-108 doi: 10.17265/2159-5275/2016.02.006 D DAVID PUBLISHING Rock Cutting Analysis Employing Finite and Discrete Element Methods Carla Massignani
More information2D Embankment and Slope Analysis (Numerical)
2D Embankment and Slope Analysis (Numerical) Page 1 2D Embankment and Slope Analysis (Numerical) Sunday, August 14, 2011 Reading Assignment Lecture Notes Other Materials FLAC Manual 1. 2. Homework Assignment
More information3D simulations of an injection test done into an unsaturated porous and fractured limestone
3D simulations of an injection test done into an unsaturated porous and fractured limestone A. Thoraval *, Y. Guglielmi, F. Cappa INERIS, Ecole des Mines de Nancy, FRANCE *Corresponding author: Ecole des
More informationON THE FACE STABILITY OF TUNNELS IN WEAK ROCKS
33 rd 33 Annual rd Annual General General Conference conference of the Canadian of the Canadian Society for Society Civil Engineering for Civil Engineering 33 e Congrès général annuel de la Société canadienne
More informationReservoir Geomechanics with ABAQUS
Reservoir Geomechanics with ABAQUS B. Bostrøm and E. Skomedal Statoil ASA, Norway Abstract: The coupled hydro-mechanical behavior of two North Sea high-pressure/hightemperature gas-condensate fields during
More informationIntegrating Geomechanics and Reservoir Characterization Examples from Canadian Shale Plays
Integrating Geomechanics and Reservoir Characterization Examples from Canadian Shale Plays AAPG Geosciences Technology Workshops Geomechanics and Reservoir Characterization of Shale and Carbonates July
More informationReservoir Rock Properties COPYRIGHT. Sources and Seals Porosity and Permeability. This section will cover the following learning objectives:
Learning Objectives Reservoir Rock Properties Core Sources and Seals Porosity and Permeability This section will cover the following learning objectives: Explain why petroleum fluids are found in underground
More informationSTRATIGRAPHY AND HYDROTHERMAL ALTERATION ENCOUNTERED BY MONITOR WELLS COMPLETED AT NGATAMARIKI AND ORAKEI KORAKO IN 2011
STRATIGRAPHY AND HYDROTHERMAL ALTERATION ENCOUNTERED BY MONITOR WELLS COMPLETED AT NGATAMARIKI AND ORAKEI KORAKO IN 2011 Catherine Boseley 1, Greg Bignall 2, Andrew Rae 2, Isabelle Chambefort 2, Brandon
More informationIAEA SAFETY STANDARDS Geotechnical Aspects of Site Evaluation and Foundations in NPPs, NS-G-3.6
IAEA SAFETY STANDARDS Geotechnical Aspects of Site Evaluation and Foundations in NPPs, NS-G-3.6 Regional Workshop on Volcanic, Seismic, and Tsunami Hazard Assessment Related to NPP Siting Activities and
More informationDOWN-HOLE SEISMIC SURVEY AND VERTICAL ELECTRIC SOUNDINGS RABASKA PROJECT, LÉVIS, QUÉBEC. Presented to :
DOWN-HOLE SEISMIC SURVEY AND VERTICAL ELECTRIC SOUNDINGS RABASKA PROJECT, LÉVIS, QUÉBEC Presented to : TERRATECH 455, René-Lévesque Blvd. West Montreal, Québec HZ 1Z3 Presented by : GEOPHYSICS GPR INTERNATIONAL
More informationThe Effect of Stress Arching on the Permeability Sensitive Experiment in the Su Lige Gas Field
The Effect of Stress Arching on the Permeability Sensitive Experiment in the Su Lige Gas Field Fanliao Wang, Xiangfang Li, Gary Couples, Mingchuan Wang, Yiqun Zhang and Jingjing Zhao THE EFFECT OF STRESS
More information8. ELASTIC PROPERTY CORRECTIONS APPLIED TO LEG 154 SEDIMENT, CEARA RISE 1
Shackleton, N.J., Curry, W.., Richter, C., and ralower, T.J. (Eds.), 1997 Proceedings of the Ocean Drilling Program, Scientific Results, Vol. 154 8. ELSTIC PROPERTY CORRECTIONS PPLIED TO LEG 154 SEDIMENT,
More informationTable of Contents Chapter 1 Introduction to Geotechnical Engineering 1.1 Geotechnical Engineering 1.2 The Unique Nature of Soil and Rock Materials
Table of Contents Chapter 1 Introduction to Geotechnical Engineering 1.1 Geotechnical Engineering 1.2 The Unique Nature of Soil and Rock Materials 1.3 Scope of This Book 1.4 Historical Development of Geotechnical
More informationAnalysis of Fracture Propagation under Thermal Stress in Geothermal Reservoirs
Proceedings World Geothermal Congress 2015 Melbourne, Australia, 19-25 April 2015 Analysis of Fracture Propagation under Thermal Stress in Geothermal Reservoirs Ahmad Ghassemi, Sergej Tarasovs Mailing
More informationNUMERICAL ANALYSIS OF DESICCATION, SHRINKAGE AND CRACKING IN LOW PLASTICITY CLAYEY SOILS
IS Numerical - Modelling analysis of Shrinkage of desiccation, Cracking shrinkage, Porous and Media cracking in low plasticity clayey soils XIV International Conference on Computational Plasticity. Fundamentals
More informationAn Updated Numerical Model of Rotorua Geothermal Field
Proceedings World Geothermal Congress 2015 Melbourne, Australia, 19-25 April 2015 An Updated Numerical Model of Rotorua Geothermal Field Thomas M.P. Ratouis, Michael O Sullivan and John O Sullivan Department
More informationNUMERICAL ANALYSIS OF DAMAGE OF RIVER EMBANKMENT ON SOFT SOIL DEPOSIT DUE TO EARTHQUAKES WITH LONG DURATION TIME
Proceedings of the International Symposium on Engineering Lessons Learned from the 2011 Great East Japan Earthquake, March 1-4, 2012, Tokyo, Japan NUMERICAL ANALYSIS OF DAMAGE OF RIVER EMBANKMENT ON SOFT
More informationUSER S MANUAL 1D Seismic Site Response Analysis Example University of California: San Diego August 30, 2017
USER S MANUAL 1D Seismic Site Response Analysis Example http://www.soilquake.net/ucsdsoilmodels/ University of California: San Diego August 30, 2017 Table of Contents USER'S MANUAL TABLE OF CONTENTS Page
More informationGeothermEx, Inc. GEOTHERMAL RESERVOIR ASSESSMENT METHODOLOGY FOR THE SCIENTIFIC OBSERVATION HOLE PROGRAM, KILAUEA EAST RIFT ZONE, HAWAII TASK 1 REPORT
(415) 527 9876 CABLE ADDRESS- GEOTHERMEX TELEX 709152 STEAM UD FAX (415) 527-8164 Geotherm Ex, Inc. RICHMOND. CALIFORNIA 94804-5829 GEOTHERMAL RESERVOIR ASSESSMENT METHODOLOGY FOR THE SCIENTIFIC OBSERVATION
More informationPetroleum Geomechanics for Shale Gas
Petroleum Geomechanics for Shale Gas Prof. Lyesse LALOUI Chair Professor «Gaz Naturel» Ecole Polytechnique Fédérale de Lausanne, EPFL Acknowledgement: V. Favero, A. Ferrari, L. Chao Unconventional Gas
More information