Science.gov

Sample records for reservoir simulation problems

  1. Experiences with linear solvers for oil reservoir simulation problems

    SciTech Connect

    Joubert, W.; Janardhan, R.; Biswas, D.; Carey, G.

    1996-12-31

    This talk will focus on practical experiences with iterative linear solver algorithms used in conjunction with Amoco Production Company`s Falcon oil reservoir simulation code. The goal of this study is to determine the best linear solver algorithms for these types of problems. The results of numerical experiments will be presented.

  2. Geothermal reservoir simulation

    NASA Technical Reports Server (NTRS)

    Mercer, J. W., Jr.; Faust, C.; Pinder, G. F.

    1974-01-01

    The prediction of long-term geothermal reservoir performance and the environmental impact of exploiting this resource are two important problems associated with the utilization of geothermal energy for power production. Our research effort addresses these problems through numerical simulation. Computer codes based on the solution of partial-differential equations using finite-element techniques are being prepared to simulate multiphase energy transport, energy transport in fractured porous reservoirs, well bore phenomena, and subsidence.

  3. Interactive reservoir simulation

    SciTech Connect

    Regtien, J.M.M. Por, G.J.A.; Stiphout, M.T. van; Vlugt, F.F. van der

    1995-12-31

    Shell`s new Modular Reservoir Simulator (MoReS) has been equipped with a comprehensive and versatile user interface called FrontEnd. Apart from providing a user-friendly environment for interactive reservoir simulation, FrontEnd serves a software platform for other dynamic simulation and reservoir-engineering applications. It offers to all supported applications a common user interface, enables the re-use of code and reduces overall maintenance and support costs associated with the embedded applications. Because of its features, FrontEnd facilitates the transfer of research results in the form of operational software to end users. When coupled with MoReS, FrontEnd can be used for pre- and post-processing and interactive simulation. The pre-processing options allow data to be inputted by means of various OSF/Motif widgets containing a spreadsheet, text editors, dialogues and graphical input. The display of the input data as well as the post-processing of all simulation results is made possible by a variety of user-defined plot of tabular (e.g. timestep summary) and array (simulation grid) data. During a simulation user-defined plots can be displayed and edited, allowing a close inspection of the results as they are being calculated. FrontEnd has been equipped with a powerful input command language, which gives the batch user as much flexibility and control over the input as the interactive user.

  4. Trends in reservoir simulation

    SciTech Connect

    Nolen, J.S.

    1995-06-01

    The future of reservoir simulation is driven by two different and, on the surface, paradoxical trends. On the one hand, the user base is on average becoming less experienced, and on the other, increasingly complex models are being built to honor the advances in reservoir-description technology. The job of the software development community is to create software that satisfies both the ease-of-use needs of the novice and the accuracy needs of the integrated geoscience team. One of the near-term effects of these demands will be to improve the capabilities and quality of the fully integrated geoscience work-station. This will include the need for implementation of industry-wide data standards. Reservoir simulators will need to incorporate increasing amounts of interactivity and built-in expertise. Accuracy of results will be improved by increased use of unstructured grids, including automatic gridding software with dynamic capabilities. Additional research will focus on complex wells, including both in-flow performance and wellbore hydraulics. Finally, grid size will continue to escalate in step with advances in hardware and software. The growth of grid size will be mitigated by substantial efforts in upscaling, but ultimately parallel computing must provide the mechanism for continued growth.

  5. Models for naturally fractured, carbonate reservoir simulations

    SciTech Connect

    Tuncay, K.; Park, A.; Ozkan, G.; Zhan, X.; Ortoleva, P.; Hoak, T.; Sundberg, K.

    1998-12-31

    This report outlines the need for new tools for the simulation of fractured carbonate reservoirs. Several problems are identified that call for the development of new reservoir simulation physical models and numerical techniques. These include: karst and vuggy media wherein Darcy`s and traditional multi-phase flow laws do not apply; the need for predicting the preproduction state of fracturing and stress so that the later response of effective stress-dependent reservoirs can be predicted; and methods for predicting the fracturing and collapse of vuggy and karst reservoirs in response to draw-down pressure created during production. Specific research directions for addressing each problem are outlined and preliminary results are noted.

  6. Integration of reservoir simulation and geomechanics

    NASA Astrophysics Data System (ADS)

    Zhao, Nan

    Fluid production from tight and shale gas formations has increased significantly, and this unconventional portfolio of low-permeability reservoirs accounts for more than half of the gas produced in the United States. Stimulation and hydraulic fracturing are critical in making these systems productive, and hence it is important to understand the mechanics of the reservoir. When modeling fractured reservoirs using discrete-fracture network representation, the geomechanical effects are expected to have a significant impact on important reservoir characteristics. It has become more accepted that fracture growth, particularly in naturally fractured reservoirs with extremely low permeability, cannot be reliably represented by conventional planar representations. Characterizing the evolution of multiple, nonplanar, interconnected and possibly nonvertical hydraulic fractures requires hydraulic and mechanical characterization of the matrix, as well as existing latent or healed fracture networks. To solve these challenging problems, a reservoir simulator (Advanced Reactive Transport Simulator (ARTS)) capable of performing unconventional reservoir simulation is developed in this research work. A geomechanical model has been incorporated into the simulation framework with various coupling schemes and this model is used to understand the geomechanical effects in unconventional oil and gas recovery. This development allows ARTS to accept geomechanical information from external geomechanical simulators (soft coupling) or the solution of the geomechanical coupled problem (hard coupling). An iterative solution method of the flow and geomechanical equations has been used in implementing the hard coupling scheme. The hard coupling schemes were verified using one-dimensional and two-dimensional analytical solutions. The new reservoir simulator is applied to learn the influence of geomechanical impact on unconventional oil and gas production in a number of practical recovery scenarios

  7. Petroleum reservoir data for testing simulation models

    SciTech Connect

    Lloyd, J.M.; Harrison, W.

    1980-09-01

    This report consists of reservoir pressure and production data for 25 petroleum reservoirs. Included are 5 data sets for single-phase (liquid) reservoirs, 1 data set for a single-phase (liquid) reservoir with pressure maintenance, 13 data sets for two-phase (liquid/gas) reservoirs and 6 for two-phase reservoirs with pressure maintenance. Also given are ancillary data for each reservoir that could be of value in the development and validation of simulation models. A bibliography is included that lists the publications from which the data were obtained.

  8. Next generation oil reservoir simulations

    SciTech Connect

    Joubert, W.

    1996-04-01

    This paper describes a collaborative effort between Amoco Production Company, Los Alamos National Laboratory and Cray Research Inc. to develop a next-generation massively parallel oil reservoir simulation code. The simulator, code-named Falcon, enables highly detailed simulations to be performed on a range of platforms such as the Cray T3D and T3E. The code is currently being used by Amoco to perform a sophisticated field study using multiple geostatistical realizations on a scale of 2-5 million grid blocks and 1000-2000 wells. In this paper we discuss the nature of this collaborative effort, the software design and engineering aspects of the code, parallelization experiences, and performance studies. The code will be marketed to the oil industry by a third-party independent software vendor in mid-1996.

  9. Multigrid methods with applications to reservoir simulation

    SciTech Connect

    Xiao, Shengyou

    1994-05-01

    Multigrid methods are studied for solving elliptic partial differential equations. Focus is on parallel multigrid methods and their use for reservoir simulation. Multicolor Fourier analysis is used to analyze the behavior of standard multigrid methods for problems in one and two dimensions. Relation between multicolor and standard Fourier analysis is established. Multiple coarse grid methods for solving model problems in 1 and 2 dimensions are considered; at each coarse grid level we use more than one coarse grid to improve convergence. For a given Dirichlet problem, a related extended problem is first constructed; a purification procedure can be used to obtain Moore-Penrose solutions of the singular systems encountered. For solving anisotropic equations, semicoarsening and line smoothing techniques are used with multiple coarse grid methods to improve convergence. Two-level convergence factors are estimated using multicolor. In the case where each operator has the same stencil on each grid point on one level, exact multilevel convergence factors can be obtained. For solving partial differential equations with discontinuous coefficients, interpolation and restriction operators should include information about the equation coefficients. Matrix-dependent interpolation and restriction operators based on the Schur complement can be used in nonsymmetric cases. A semicoarsening multigrid solver with these operators is used in UTCOMP, a 3-D, multiphase, multicomponent, compositional reservoir simulator. The numerical experiments are carried out on different computing systems. Results indicate that the multigrid methods are promising.

  10. Reservoir Simulations of Low-Temperature Geothermal Reservoirs

    NASA Astrophysics Data System (ADS)

    Bedre, Madhur Ganesh

    The eastern United States generally has lower temperature gradients than the western United States. However, West Virginia, in particular, has higher temperature gradients compared to other eastern states. A recent study at Southern Methodist University by Blackwell et al. has shown the presence of a hot spot in the eastern part of West Virginia with temperatures reaching 150°C at a depth of between 4.5 and 5 km. This thesis work examines similar reservoirs at a depth of around 5 km resembling the geology of West Virginia, USA. The temperature gradients used are in accordance with the SMU study. In order to assess the effects of geothermal reservoir conditions on the lifetime of a low-temperature geothermal system, a sensitivity analysis study was performed on following seven natural and human-controlled parameters within a geothermal reservoir: reservoir temperature, injection fluid temperature, injection flow rate, porosity, rock thermal conductivity, water loss (%) and well spacing. This sensitivity analysis is completed by using ‘One factor at a time method (OFAT)’ and ‘Plackett-Burman design’ methods. The data used for this study was obtained by carrying out the reservoir simulations using TOUGH2 simulator. The second part of this work is to create a database of thermal potential and time-dependant reservoir conditions for low-temperature geothermal reservoirs by studying a number of possible scenarios. Variations in the parameters identified in sensitivity analysis study are used to expand the scope of database. Main results include the thermal potential of reservoir, pressure and temperature profile of the reservoir over its operational life (30 years for this study), the plant capacity and required pumping power. The results of this database will help the supply curves calculations for low-temperature geothermal reservoirs in the United States, which is the long term goal of the work being done by the geothermal research group under Dr. Anderson at

  11. Large eddy simulation of a pumped- storage reservoir

    NASA Astrophysics Data System (ADS)

    Launay, Marina; Leite Ribeiro, Marcelo; Roman, Federico; Armenio, Vincenzo

    2016-04-01

    The last decades have seen an increasing number of pumped-storage hydropower projects all over the world. Pumped-storage schemes move water between two reservoirs located at different elevations to store energy and to generate electricity following the electricity demand. Thus the reservoirs can be subject to important water level variations occurring at the daily scale. These new cycles leads to changes in the hydraulic behaviour of the reservoirs. Sediment dynamics and sediment budgets are modified, sometimes inducing problems of erosion and deposition within the reservoirs. With the development of computer performances, the use of numerical techniques has become popular for the study of environmental processes. Among numerical techniques, Large Eddy Simulation (LES) has arisen as an alternative tool for problems characterized by complex physics and geometries. This work uses the LES-COAST Code, a LES model under development in the framework of the Seditrans Project, for the simulation of an Upper Alpine Reservoir of a pumped-storage scheme. Simulations consider the filling (pump mode) and emptying (turbine mode) of the reservoir. The hydraulic results give a better understanding of the processes occurring within the reservoir. They are considered for an assessment of the sediment transport processes and of their consequences.

  12. A general formulation for compositional reservoir simulation

    SciTech Connect

    Rodriguez, F.; Guzman, J.; Galindo-Nava, A.

    1994-12-31

    In this paper the authors present a general formulation to solve the non-linear difference equations that arise in compositional reservoir simulation. The general approach here presented is based on newton`s method and provides a systematic approach to generate several formulations to solve the compositional problem, each possessing a different degree of implicitness and stability characteristics. The Fully-Implicit method is at the higher end of the implicitness spectrum while the IMPECS method, implicit in pressure-explicit in composition and saturation, is at the lower end. They show that all methods may be obtained as particular cases of the fully-implicit method. Regarding the matrix problem, all methods have a similar matrix structure; the composition of the Jacobian matrix is however unique in each case, being in some instances amenable to reductions for optimal solution of the matrix problem. Based on this, a different approach to derive IMPECS type methods is proposed; in this case, the whole set of 2nc + 6 equations, that apply in each gridblock, is reduced to a single pressure equation through matrix reduction operations; this provides a more stable numerical scheme, compared to other published IMPCS methods, in which the subset of thermodynamic equilibrium equations is arbitrarily decoupled form the set of gridblock equations to perform such reduction. The authors discuss how the general formulation here presented can be used to formulate and construct an adaptive-implicit compositional simulators. They also present results on the numerical performance of FI, IMPSEC and IMPECS methods on some test problems.

  13. A finite element simulation system in reservoir engineering

    SciTech Connect

    Gu, Xiaozhong

    1996-03-01

    Reservoir engineering is performed to predict the future performance of a reservoir based on its current state and past performance and to explore other methods for increasing the recovery of hydrocarbons from a reservoir. Reservoir simulations are routinely used for these purposes. A reservoir simulator is a sophisticated computer program which solves a system of partial differential equations describing multiphase fluid flow (oil, water, and gas) in a porous reservoir rock. This document describes the use of a reservoir simulator version of BOAST which was developed by the National Institute for Petroleum and Energy Research in July, 1991.

  14. Geomechanically Coupled Simulation of Flow in Fractured Reservoirs

    NASA Astrophysics Data System (ADS)

    Barton, C.; Moos, D.; Hartley, L.; Baxter, S.; Foulquier, L.; Holl, H.; Hogarth, R.

    2012-12-01

    Capturing the necessary and sufficient detail of reservoir hydraulics to accurately evaluate reservoir behavior remains a significant challenge to the exploitation and management of fracture-dominated geothermal reservoirs. In these low matrix permeability reservoirs, stimulation response is controlled largely by the properties of natural and induced fracture networks, which are in turn controlled by the in situ stresses, the fracture distribution and connectivity and the hydraulic behavior of the fractures. This complex interaction of fracture flow systems with the present-day stress field compounds the problem of developing an effective and efficient simulation to characterize, model and predict fractured reservoir performance. We discuss here a case study of the integration of geological, geophysical, geomechanical, and reservoir engineering data to characterize the in situ stresses, the natural fracture network and the controls on fracture permeability in geothermal reservoirs. A 3D geomechanical reservoir model includes constraints on stress magnitudes and orientations, and constraints on mechanical rock properties and the fractures themselves. Such a model is essential to understanding reservoir response to stimulation and production in low matrix permeability, fracture-dominated reservoirs. The geomechanical model for this study was developed using petrophysical, drilling, and wellbore image data along with direct well test measurements and was mapped to a 3D structural grid to facilitate coupled simulation of the fractured reservoir. Wellbore image and stimulation test data were used along with microseismic data acquired during the test to determine the reservoir fracture architecture and to provide control points for a realistic inter-connected discrete fracture network. As most fractures are stress-sensitive, their hydraulic conductivities will change with changes in bottomhole flowing and reservoir pressures, causing variations in production profiles

  15. Rabi multi-sector reservoir simulation model

    SciTech Connect

    Bruijnzeels, C.; O`Halloran, C.

    1995-12-31

    To ensure optimum ultimate recovery of the 46 meter thick oil rim of the Rabi Field in Gabon, a full field simulation model was required. Due to it`s size and complexity, with local cusping, coning and geological circumstances dominating individual well behavior, a single full field model would be too large for existing hardware. A method was developed to simulate the full field with 5 separate sector models, whilst allowing the development in one sector model to have an effect on the boundary conditions of another sector. In this manner, the 13 x 4.5 km field could be simulated with a horizontal well spacing down to 175 meter. This paper focuses on the method used to attach single 3-phase tank cells to a sector simulation grid in order to represent non-simulated parts of the field. It also describes the history matching methodology and how to run a multisector model in forecasting mode. This method can be used for any reservoir, where size and complexity require large reservoir simulation models that normally could not be modeled within the constraints of available computer facilities. Detailed studies can be conducted on specific parts of a field, whilst allowing for dynamic flow and pressure effects caused by the rest of the field.

  16. A flexible gridding scheme for reservoir simulation

    SciTech Connect

    Verma, S.

    1995-12-31

    This paper describes a new control volume based finite difference scheme for petroleum reservoir simulation which can be used with unstructured grids. The numerical scheme to model fluid flow is shown to be easily used for Voronoi grids in 2D. It can also be used with certain geometrical limitations for 3D Voronoi grids. The scheme can be used without any significant limitations for triangle or tetrahedron based grids where control volumes are constructed around their vertices. It assumes uniform properties inside such control volumes. Full, anisotropic and asymmetric permeability tensor can be easily handled with the proposed method. The permeability tensor can vary from block to block. Thus it will be of great value in modeling fluid flow in reservoirs where principal directions of permeability varies between beds or within a bed. The paper also presents an analysis of some of the published flexible gridding schemes which use a control volume type algebraic approximation and demonstrate the advantages of the method presented here. The technique for grid construction is also discussed. Test results presented here demonstrate the need for proper representation of reservoir geometry to predict the correct flow behavior. The gridding scheme described in this paper achieves that purpose.

  17. Iterative Schemes for Time Parallelization with Application to Reservoir Simulation

    SciTech Connect

    Garrido, I; Fladmark, G E; Espedal, M S; Lee, B

    2005-04-18

    Parallel methods are usually not applied to the time domain because of the inherit sequentialness of time evolution. But for many evolutionary problems, computer simulation can benefit substantially from time parallelization methods. In this paper, they present several such algorithms that actually exploit the sequential nature of time evolution through a predictor-corrector procedure. This sequentialness ensures convergence of a parallel predictor-corrector scheme within a fixed number of iterations. The performance of these novel algorithms, which are derived from the classical alternating Schwarz method, are illustrated through several numerical examples using the reservoir simulator Athena.

  18. ADVANCED TECHNIQUES FOR RESERVOIR SIMULATION AND MODELING OF NONCONVENTIONAL WELLS

    SciTech Connect

    Louis J. Durlofsky; Khalid Aziz

    2004-08-20

    Nonconventional wells, which include horizontal, deviated, multilateral and ''smart'' wells, offer great potential for the efficient management of oil and gas reservoirs. These wells are able to contact larger regions of the reservoir than conventional wells and can also be used to target isolated hydrocarbon accumulations. The use of nonconventional wells instrumented with downhole inflow control devices allows for even greater flexibility in production. Because nonconventional wells can be very expensive to drill, complete and instrument, it is important to be able to optimize their deployment, which requires the accurate prediction of their performance. However, predictions of nonconventional well performance are often inaccurate. This is likely due to inadequacies in some of the reservoir engineering and reservoir simulation tools used to model and optimize nonconventional well performance. A number of new issues arise in the modeling and optimization of nonconventional wells. For example, the optimal use of downhole inflow control devices has not been addressed for practical problems. In addition, the impact of geological and engineering uncertainty (e.g., valve reliability) has not been previously considered. In order to model and optimize nonconventional wells in different settings, it is essential that the tools be implemented into a general reservoir simulator. This simulator must be sufficiently general and robust and must in addition be linked to a sophisticated well model. Our research under this five year project addressed all of the key areas indicated above. The overall project was divided into three main categories: (1) advanced reservoir simulation techniques for modeling nonconventional wells; (2) improved techniques for computing well productivity (for use in reservoir engineering calculations) and for coupling the well to the simulator (which includes the accurate calculation of well index and the modeling of multiphase flow in the wellbore

  19. Geothermal reservoir engineering computer code comparison and validation calculations using MUSHRM and CHARGR geothermal reservoir simulators

    SciTech Connect

    Pritchett, J.W.

    1980-11-01

    The essential features of the reservoir codes CHARGR and MUSHRM are described. Solutions obtained for the problem set posed by DOE are presented. CHARGR was used for all six problems; MUSHRM was used for one. These problems are: the 1-D Avdonin solution, the 1-D well test analysis, 2-D flow to a well in fracture/block media, expanding two-phase system with drainage, flow in a 2-D areal reservoir, and flow in a 3-D reservoir. Results for the last problem using both codes are compared. (MHR)

  20. Simulation of California's Major Reservoirs Outflow Using Data Mining Technique

    NASA Astrophysics Data System (ADS)

    Yang, T.; Gao, X.; Sorooshian, S.

    2014-12-01

    The reservoir's outflow is controlled by reservoir operators, which is different from the upstream inflow. The outflow is more important than the reservoir's inflow for the downstream water users. In order to simulate the complicated reservoir operation and extract the outflow decision making patterns for California's 12 major reservoirs, we build a data-driven, computer-based ("artificial intelligent") reservoir decision making tool, using decision regression and classification tree approach. This is a well-developed statistical and graphical modeling methodology in the field of data mining. A shuffled cross validation approach is also employed to extract the outflow decision making patterns and rules based on the selected decision variables (inflow amount, precipitation, timing, water type year etc.). To show the accuracy of the model, a verification study is carried out comparing the model-generated outflow decisions ("artificial intelligent" decisions) with that made by reservoir operators (human decisions). The simulation results show that the machine-generated outflow decisions are very similar to the real reservoir operators' decisions. This conclusion is based on statistical evaluations using the Nash-Sutcliffe test. The proposed model is able to detect the most influential variables and their weights when the reservoir operators make an outflow decision. While the proposed approach was firstly applied and tested on California's 12 major reservoirs, the method is universally adaptable to other reservoir systems.

  1. A structurally complex and dynamic reservoir description for reservoir simulation, Kuparuk River Field, Alaska

    SciTech Connect

    Walsh, T.P.; Leander, M.H.; Wilcox, T.C.

    1995-08-01

    The Kupanuk River Field is a structurally complex giant oil field adjacent to the Prudhoe Bay Field on Alaska`s North Slope. Oil is reservoired within two Early Cretaceous shallow marine sandstone formations, separated stratigraphically by an erosionally truncated marine silt/shale. Subjected to several phases of tectonism, this highly compartmentalized reservoir has been developed on regular 160 acre direct line drive patterns. An integrated team of geoscientists and engineers from BP Exploration (Alaska) Inc. and ARCO Alaska Inc. is presently quantifying the benefits of infill drilling at Kuparuk, and identifying the best locations for well placement. The two primary reservoir characteristics believed to impact the effectiveness of infill drilling are large-scale reservoir heterogeneity, and reservoir comparmentation due to faulting. Multiple thin pay zones within the two reservoir intervals are isolated laterally by faults with magnitudes greater than pay zone thickness. A process and tools designed to construct and maintain a structurally complex reservoir description, shared by the geoscientists and reservoir engineers, are described. Cross-discipline integration is aided by the use of Tech*Logic`s IREX 3-D reservoir modeling and visualization application. The unique architecture of the IREX model allows for representation of very complex structural geometries, and facilitates iteration between reservoir description and simulation, along the seismic to simulation continuum. Modifications to the reservoir description are guided by well-level history matching within the constraints of all available geoscience information. The techniques described will be of particular interest to those working on reservoir description and simulation of structurally complex fields.

  2. Feasibility study of sedimentary enhanced geothermal systems using reservoir simulation

    NASA Astrophysics Data System (ADS)

    Cho, Jae Kyoung

    The objective of this research is to evaluate the preliminary feasibility of commercial geothermal projects, from a sedimentary reservoir with low permeability that requires productivity enhancement, using numerical reservoir simulation. The performance of a sedimentary geothermal reservoir is investigated in terms of reservoir hydraulics and thermal evolution. To build a reliable benchmark for simulation study, validation of the numerical reservoir model with respect to an analytical model is presented, and the process to achieve an acceptable match between the numerical and analytical solutions is described. The analytical model used in this study is based on the work of Gringarten (1978), which consists of a conceptual geothermal reservoir, considering an injection and production well doublet in a homogeneous porous media. A commercial thermal reservoir simulator (STARS from Computer Modeling Group, CMG) is used in this work for numerical modeling. In order to reproduce the analytical model results, the numerical simulation model is modified to include the same assumptions of the analytical model. Simulation model parameters that make the numerical results deviate from the analytical solution, such as the grid block size, time step and no-flow boundary are identified and investigated. An analytical tracer test model proposed by Shook (2000) is numerically modeled. This model allows us to predict the time when the temperature of the produced water decreases by capturing a tracer component at production well. Reservoir simulation models with different porosity and permeability distribution are tested to see the effects of reservoir inhomogeneity and anisotropy. In particular, premature thermal breakthrough due to the presence of high permeability streak in a reservoir model is simulated. In an effort to apply the knowledge we obtained from the analytical solutions, the effects of reservoir rock and water properties, as a function of pressure and temperature, are

  3. Cooperative Learning in Reservoir Simulation Classes: Overcoming Disparate Entry Skills

    ERIC Educational Resources Information Center

    Awang, Mariyamni

    2006-01-01

    Reservoir simulation is one of the core courses in the petroleum engineering curriculum and it requires knowledge and skills in three major disciplines, namely programming, numerical methods and reservoir engineering. However, there were often gaps in the students' readiness to undertake the course, even after completing the necessary…

  4. Estimation of discontinuous coefficients in parabolic systems: Applications to reservoir simulation

    NASA Technical Reports Server (NTRS)

    Lamm, P. D.

    1984-01-01

    Spline based techniques for estimating spatially varying parameters that appear in parabolic distributed systems (typical of those found in reservoir simulation problems) are presented. The problem of determining discontinuous coefficients, estimating both the functional shape and points of discontinuity for such parameters is discussed. Convergence results and a summary of numerical performance of the resulting algorithms are given.

  5. New simulators quantify slimhole effects on reservoir performance

    SciTech Connect

    Azari, M.; Soliman, M.

    1995-03-01

    According to a recent slimhole survey, 60% of the slimhole wells drilled in 1992 were considered an economic success. The remaining 40% of wells drilled, however, either had not been in service long enough for their success to be determined or had failed to provide the expected production results. To more accurately predict production results before choosing to drill a slimhole well, drilling operators must answer the following questions: (1) How would the rate performance and cumulative recovery of a slimhole well compare to a standard well bore completion? (2) How would the reduced tubing ID of a slim hole affect well bore hydraulics? (3) What effect would hole diameter have on the well bore storage period? (4) Will a slimhole well more likely experience serious coning or cresting problems? Reservoir engineers can use several powerful simulators and software programs to determine well performance.

  6. Self-priming hemodynamic reservoir and inline flow meter for a cardiopulmonary bypass simulation.

    PubMed

    Raasch, David; Austin, Jon; Tallman, Richard

    2010-06-01

    Simulator exercises are used at Midwestern University to augment academic and laboratory training toward consolidating particular skills, increasing situation awareness, and preparing the student for practice within the team environment of an operating room. This paper describes an enhanced cardiopulmonary bypass simulator consisting of a self-priming hemodynamic reservoir that includes an inline flow meter. A typical cardiopulmonary bypass adult perfusion circuit was assembled using a roller pump console and integrated oxygenator/heat exchanger/reservoir and primed with 2 liters of water. For patient simulation, a soft-sided reservoir bag was mounted onto an inclined platform. A 1-liter soft-sided bag was placed just above the reservoir, providing an overflow reservoir. The priming line extended to the head of the mannequin. The arterial, venous, and suction lines extended through the open chest. The primed perfusion circuit was connected to ports on the filled reservoir bag. To test the patient simulation, the arterial pump output was adjusted to flow rates ranging from 1-7 liters per minute, with a complete interruption (to zero flow) between each test run. An inline flow meter was added to the bypass circuit and an analog to digital converter board was used to pass flow data into the computer-based simulation program. The use of an inclined hemodynamic reservoir bag proved to be self-priming and functional without problems over a wide range of flows tested. By including a reservoir with the mannequin, plus processing and displaying real-time flow data using the CPB-Sim simulation program, a higher fidelity and more realistic simulation experience was created. PMID:20648900

  7. Optimality conditions for a two-stage reservoir operation problem

    NASA Astrophysics Data System (ADS)

    Zhao, Jianshi; Cai, Ximing; Wang, Zhongjing

    2011-08-01

    This paper discusses the optimality conditions for standard operation policy (SOP) and hedging rule (HR) for a two-stage reservoir operation problem using a consistent theoretical framework. The effects of three typical constraints, i.e., mass balance, nonnegative release, and storage constraints under both certain and uncertain conditions are analyzed. When all nonnegative constraints and storage constraints are unbinding, HR results in optimal reservoir operation following the marginal benefit (MB) principle (the MB is equal over current and future stages. However, if any of those constraints are binding, SOP results in the optimal solution, except in some special cases which need to carry over water in the current stage to the future stage, when extreme drought is certain and a higher marginal utility exists for the second stage. Furthermore, uncertainty complicates the effects of the various constraints. A higher uncertainty level in the future makes HR more favorable as water needs to be reserved to defend against the risk caused by uncertainty. Using the derived optimality conditions, an algorithm for solving a numerical model is developed and tested with the Miyun Reservoir in China.

  8. The sequential method for the black-oil reservoir simulation on unstructured grids

    NASA Astrophysics Data System (ADS)

    Li, Baoyan; Chen, Zhangxin; Huan, Guanren

    2003-11-01

    This paper presents new results for applying the sequential solution method to the black-oil reservoir simulation with unstructured grids. The fully implicit solution method has been successfully applied to reservoir simulation with unstructured grids. However, the complexity of the fully implicit method and the irregularity of the grids result in a very complicated structure of linear equation systems (LESs) and in high computational cost to solve them. To tackle this problem, the sequential method is applied to reduce the size of the LESs. To deal with instable problems caused by the low implicit degree of this method, some practical techniques are introduced to control convergence of Newton-Raphson's iterations which are exploited in the linearization of the governing equations of the black-oil model. These techniques are tested with the benchmark problem of the ninth comparative solution project (CSP) organized by the society of petroleum engineers (SPE) and applied to field-scale models of both saturated and undersaturated reservoirs. The simulation results show that the sequential method uses as little as 20.01% of the memory for solving the LESs and 23.89% of the total computational time of the fully implicit method to reach the same precision for the undersaturated reservoirs, when the same iteration control parameters are used for both solution methods. However, for the saturated reservoirs, the sequential method must use stricter iteration control parameters to reach the same precision as the fully implicit method.

  9. Hydrodynamic modeling of petroleum reservoirs using simulator MUFITS

    NASA Astrophysics Data System (ADS)

    Afanasyev, Andrey

    2015-04-01

    MUFITS is new noncommercial software for numerical modeling of subsurface processes in various applications (www.mufits.imec.msu.ru). To this point, the simulator was used for modeling nonisothermal flows in geothermal reservoirs and for modeling underground carbon dioxide storage. In this work, we present recent extension of the code to petroleum reservoirs. The simulator can be applied in conventional black oil modeling, but it also utilizes a more complicated models for volatile oil and gas condensate reservoirs as well as for oil rim fields. We give a brief overview of the code by providing the description of internal representation of reservoir models, which are constructed of grid blocks, interfaces, stock tanks as well as of pipe segments and pipe junctions for modeling wells and surface networks. For conventional black oil approach, we present the simulation results for SPE comparative tests. We propose an accelerated compositional modeling method for sub- and supercritical flows subjected to various phase equilibria, particularly to three-phase equilibria of vapour-liquid-liquid type. The method is based on the calculation of the thermodynamic potential of reservoir fluid as a function of pressure, total enthalpy and total composition and storing its values as a spline table, which is used in hydrodynamic simulation for accelerated PVT properties prediction. We provide the description of both the spline calculation procedure and the flashing algorithm. We evaluate the thermodynamic potential for a mixture of two pseudo-components modeling the heavy and light hydrocarbon fractions. We develop a technique for converting black oil PVT tables to the potential, which can be used for in-situ hydrocarbons multiphase equilibria prediction under sub- and supercritical conditions, particularly, in gas condensate and volatile oil reservoirs. We simulate recovery from a reservoir subject to near-critical initial conditions for hydrocarbon mixture. We acknowledge

  10. A chemical EOR benchmark study of different reservoir simulators

    NASA Astrophysics Data System (ADS)

    Goudarzi, Ali; Delshad, Mojdeh; Sepehrnoori, Kamy

    2016-09-01

    Interest in chemical EOR processes has intensified in recent years due to the advancements in chemical formulations and injection techniques. Injecting Polymer (P), surfactant/polymer (SP), and alkaline/surfactant/polymer (ASP) are techniques for improving sweep and displacement efficiencies with the aim of improving oil production in both secondary and tertiary floods. There has been great interest in chemical flooding recently for different challenging situations. These include high temperature reservoirs, formations with extreme salinity and hardness, naturally fractured carbonates, and sandstone reservoirs with heavy and viscous crude oils. More oil reservoirs are reaching maturity where secondary polymer floods and tertiary surfactant methods have become increasingly important. This significance has added to the industry's interest in using reservoir simulators as tools for reservoir evaluation and management to minimize costs and increase the process efficiency. Reservoir simulators with special features are needed to represent coupled chemical and physical processes present in chemical EOR processes. The simulators need to be first validated against well controlled lab and pilot scale experiments to reliably predict the full field implementations. The available data from laboratory scale include 1) phase behavior and rheological data; and 2) results of secondary and tertiary coreflood experiments for P, SP, and ASP floods under reservoir conditions, i.e. chemical retentions, pressure drop, and oil recovery. Data collected from corefloods are used as benchmark tests comparing numerical reservoir simulators with chemical EOR modeling capabilities such as STARS of CMG, ECLIPSE-100 of Schlumberger, REVEAL of Petroleum Experts. The research UTCHEM simulator from The University of Texas at Austin is also included since it has been the benchmark for chemical flooding simulation for over 25 years. The results of this benchmark comparison will be utilized to improve

  11. A CUDA based parallel multi-phase oil reservoir simulator

    NASA Astrophysics Data System (ADS)

    Zaza, Ayham; Awotunde, Abeeb A.; Fairag, Faisal A.; Al-Mouhamed, Mayez A.

    2016-09-01

    Forward Reservoir Simulation (FRS) is a challenging process that models fluid flow and mass transfer in porous media to draw conclusions about the behavior of certain flow variables and well responses. Besides the operational cost associated with matrix assembly, FRS repeatedly solves huge and computationally expensive sparse, ill-conditioned and unsymmetrical linear system. Moreover, as the computation for practical reservoir dimensions lasts for long times, speeding up the process by taking advantage of parallel platforms is indispensable. By considering the state of art advances in massively parallel computing and the accompanying parallel architecture, this work aims primarily at developing a CUDA-based parallel simulator for oil reservoir. In addition to the initial reported 33 times speed gain compared to the serial version, running experiments showed that BiCGSTAB is a stable and fast solver which could be incorporated in such simulations instead of the more expensive, storage demanding and usually utilized GMRES.

  12. A simulation-based reservoir management program

    SciTech Connect

    Voskanian, M.M.; Kendall, R.P.; Whitney, E.M.; Coombs, S.; Paul, R.G.; Ershaghi, I.

    1996-05-01

    There are more than 5,200 independent oil and gas producers operating in the US today (based on current IPAA membership figures). These companies are playing an increasingly important role in production of hydrocarbons in California and elsewhere in the US. Pacific Operators Offshore, Inc., in a historic collaboration with its government royalty owners, the California State Lands Commission and the Minerals Management Service of the US Department of Interior, is attempting to redevelop the Carpinteria Offshore Field after two-and-a-half decades of production and partial abandonment by a previous operator. This paper will describe a project which focuses on the distribution of advanced reservoir management technologies (geological, petrophysical, and engineering) to independent producers like Pacific Operators Offshore, Inc. The evolving information highway, specifically the World Wide Web (WWW), serves as the distribution medium. The project to be described in this paper is an example of the implementation of a reservoir management tool which is supported by distributed databases, incorporates a shared computing environment, and integrates stochastic, geological, and engineering modeling.

  13. Large-scale three-dimensional geothermal reservoir simulation on PCs

    SciTech Connect

    Antunez, E.; Moridis, G.; Pruess, K.

    1994-01-01

    TOUGH2, Lawrence Berkeley Laboratory`s general purpose simulator for mass and heat flow and transport was enhanced with the addition of a set of preconditioned conjugate gradient solvers and ported to a PC. The code was applied to a number of large 3-D geothermal reservoir problems with up to 10,000 grid blocks. Four test problems were investigated. The first two involved a single-phase liquid system, and a two-phase system with regular Cartesian grids. The last two involved a two-phase field problem with irregular gridding with production from and injection into a single porosity reservoir, and a fractured reservoir. The code modifications to TOUGH2 and its setup in the PC environment are described. Algorithms suitable for solving large matrices that are generally non-symmetric and non-positive definite are reviewed. Computational work per time step and CPU time requirements are reported as function of problem size. The excessive execution time and storage requirements of the direct solver in TOUGH2 limits the size of manageable 3-D reservoir problems to a few hundred grid blocks. The conjugate gradient solvers significantly reduced the execution time and storage requirements making possible the execution of considerably larger problems (10,000 + grid blocks). It is concluded that the current PCs provide an economical platform for running large-scale geothermal field simulations that just a few years ago could only be executed on mainframe computers.

  14. Large-scale three-dimensional geothermal reservoir simulation on PCs

    SciTech Connect

    Antunez, Emilio; Moridis, George; Pruess, Karsten

    1994-01-20

    TOUGH2, Lawrence Berkeley Laboratory's general purpose simulator for mass and heat flow and transport was enhanced with the addition of a set of preconditioned conjugate gradient solvers and ported to a PC. The code was applied to a number of large 3-D geothermal reservoir problems with up to 10,000 grid blocks. Four test problems were investigated. The first two involved a single-phase liquid system, and a two-phase system with regular Cartesian grids. The last two involved a two-phase field problem with irregular gridding with production from and injection into a single porosity reservoir, and a fractured reservoir. The code modifications to TOUGH2 and its setup in the PC environment are described. Algorithms suitable for solving large matrices that are generally non-symmetric and non-positive definite are reviewed. Computational work per time step and CPU time requirements are reported as function of problem size. The excessive execution time and storage requirements of the direct solver in TOUGH2 limits the size of manageable 3-D reservoir problems to a few hundred grid blocks. The conjugate gradient solvers significantly reduced the execution time and storage requirements making possible the execution of considerably larger problems (10,000+ grid blocks). It is concluded that the current PCs provide an economical platform for running large-scale geothermal field simulations that just a few years ago could only be executed on mainframe computers.

  15. Large-scale three-dimensional geothermal reservoir simulation on PC's

    NASA Astrophysics Data System (ADS)

    Antunez, E.; Moridis, G.; Pruess, K.

    1994-01-01

    TOUGH2, Lawrence Berkeley Laboratory's general purpose simulator for mass and heat flow and transport, was enhanced with the addition of a set of preconditioned conjugate gradient solvers and ported to a PC. The code was applied to a number of large three dimensional (3-D) geothermal reservoir problems with up to 10,000 grid blocks. Four test problems were investigated. The first two involved a single-phase liquid system, and a two-phase system with regular Cartesian grids. The last two involved a two-phase field problem with irregular gridding with production from and injection into a single porosity reservoir, and a fractured reservoir. The code modifications to TOUGH2 and its setup in the PC environment are described. Algorithms suitable for solving large matrices that are generally nonsymmetric and nonpositive definite are reviewed. Computational work per time step and CPU time requirements are reported as function of problem size. The excessive execution time and storage requirements of the direct solver in TOUGH2 limits the size of manageable 3-D reservoir problems to a few hundred grid blocks. The conjugate gradient solvers significantly reduced the execution time and storage requirements making possible the execution of considerably larger problems (10,000+ grid blocks). It is concluded that the current PC's provide an economical platform for running large-scale geothermal field simulations that just a few years ago could only be executed on mainframe computers.

  16. Optimality Conditions for A Two-Stage Reservoir Operation Problem

    NASA Astrophysics Data System (ADS)

    Zhao, J.; Cai, X.; Wang, Z.

    2010-12-01

    This paper discusses the optimality conditions for standard operation policy (SOP) and hedging rule (HR) for a two-stage reservoir operation problem within a consistent theoretical framework. The effects of three typical constraints, which are mass balance, non-negative release and storage constraints under both certain and uncertain conditions have been analyzed. When all non-negative constraints and storage constraints are non-binding, HR results in optimal reservoir operation following the marginal benefit (MB) principle (the MB is equal over the two stages); while if any of the non-negative release or storage constraints is binding, in general SOP results in the optimal solution except two special cases. One of them is a complement of the traditional SOP/HR curve, which happens while the capacity constraint is binding; the other is a special hedging rule, which should be employed to carry over all water in the current stage to the future, when extreme drought is certain and higher marginal utility exists for the second stage. Furthermore, uncertainty complicates the effects of the various constraints but in general higher uncertainty level in the future makes HR a more favorable since water needs to be reserved to defense the risk caused by the uncertainty. Using the derived optimality conditions, an algorithm for solving the model numerically has been developed and tested with hypothetical examples.

  17. Application of advanced reservoir characterization, simulation, and production optimization strategies to maximize recovery in slope and basin clastic reservoirs, West Texas (Delaware Basin). Technical progress report

    SciTech Connect

    Dutton, S.P.

    1996-04-30

    The objective of this project is to demonstrate that detailed reservoir characterization of slope and basin clastic reservoirs in sandstones of the Delaware Mountain Group in the Delaware Basin of West Texas and New Mexico is a cost effective way to recover a higher percentage of the original oil in place through strategic placement of infill wells and geologically based field development. Project objectives are divided into two major phases. The objectives of the reservoir characterization phase of the project are to provide a detailed understanding of the architecture and heterogeneity of two fields, the Ford Geraldine unit and Ford West field, which produce from the Bell Canyon and Cherry Canyon Formations, respectively, of the Delaware Mountain Group and to compare Bell Canyon and Cherry Canyon reservoirs. Reservoir characterization will utilize 3-D seismic data, high-resolution sequence stratigraphy, subsurface field studies, outcrop characterization, and other techniques. One the reservoir-characterization study of both field is completed, a pilot area of approximately 1 mi{sup 2} in one of the fields will be chosen for reservoir simulation. The objectives of the implementation phase of the project are to: (1) apply the knowledge gained from reservoir characterization and simulation studies to increase recovery from the pilot area; (2) demonstrate that economically significant unrecovered oil remains in geologically resolvable untapped compartments; and (3) test the accuracy of reservoir characterization and flow simulation as predictive tools in resource preservation of mature fields. A geologically designed, enhanced recovery program (CO{sub 2} flood, waterflood, or polymer flood) and well-completion program will be developed, and one to three infill well will be drilled and cored. Technical progress is summarized for: geophysical characterization; reservoir characterization; outcrop characterization; and producibility problem characterization.

  18. An adaptive nonlinear solution scheme for reservoir simulation

    SciTech Connect

    Lett, G.S.

    1996-12-31

    Numerical reservoir simulation involves solving large, nonlinear systems of PDE with strongly discontinuous coefficients. Because of the large demands on computer memory and CPU, most users must perform simulations on very coarse grids. The average properties of the fluids and rocks must be estimated on these grids. These coarse grid {open_quotes}effective{close_quotes} properties are costly to determine, and risky to use, since their optimal values depend on the fluid flow being simulated. Thus, they must be found by trial-and-error techniques, and the more coarse the grid, the poorer the results. This paper describes a numerical reservoir simulator which accepts fine scale properties and automatically generates multiple levels of coarse grid rock and fluid properties. The fine grid properties and the coarse grid simulation results are used to estimate discretization errors with multilevel error expansions. These expansions are local, and identify areas requiring local grid refinement. These refinements are added adoptively by the simulator, and the resulting composite grid equations are solved by a nonlinear Fast Adaptive Composite (FAC) Grid method, with a damped Newton algorithm being used on each local grid. The nonsymmetric linear system of equations resulting from Newton`s method are in turn solved by a preconditioned Conjugate Gradients-like algorithm. The scheme is demonstrated by performing fine and coarse grid simulations of several multiphase reservoirs from around the world.

  19. Reservoir simulation in a North Sea reservoir experiencing significant compaction drive

    SciTech Connect

    Cook, C.C.; Jewell, S.

    1995-12-31

    The Valhall field in the Norwegian North Sea is a high porosity chalk reservoir undergoing primary pressure depletion. Over the last ten years there have been a number of computer modeling studies of the field which have all assumed an original oil-in-place of approximately 2,000 MMSTB (318.0{times}10{sup 6}m{sup 3}) to the present due to the addition of wells and the optimization of completion techniques. However, the single most important and unique feature influencing Valhall long term production performance is reservoir rock compaction. This paper describes the mathematical model used to simulate reservoir performance in a compacting reservoir with specific discussion regarding the proportion of oil produced by each physical recovery process. An understanding of the recovery mechanisms and their relative importance is critical for the successful management of the field. This paper also presents an alternative method for evaluating the various recovery processes using a simple solution to the material balance equation. This is used to substantiate the magnitude of the various recovery mechanisms identified in the simulation model.

  20. Multiscale characterization of porous media properties for hydrocarbon reservoir simulation

    NASA Astrophysics Data System (ADS)

    Neeman, Henry; Lao, Hio-Wai; Simpson, Dale; Papavassiliou, Dimitrios V.

    2001-07-01

    Fluid flow through porous materials is critical for understanding and predicting the behavior of systems as diverse in function and scale as hydrocarbon reservoirs, aquifers, filters, membrane separators and even catalytic converters. Recently, there have been calls to incorporate more physics in oil reservoir simulations, as well as to enhance computational capability through the use of High Performance Computing (HPC), in order to improve reservoir management. Accurate prediction of reservoir behavior depends on the physical properties of not only the fluid but also the underlying rock formation. Contemporary approaches to solving these flows involve simulation of only a single physical scale. We are currently developing HiMuST (Hierarchical Multiscale Simulator Technology), an integrated multiscale simulation system for flow through heterogeneous porous materials. HiMuST uses a hierarchy of simulation codes to address the issue of rock property characterization at the pore scale and can self-adjust according to available input data. At the microscopic scale, HiMuST employs the Lattice Boltzmann Method, based on magnetic resonance digitizations of actual rock samples. At the mesoscopic scale, a stochastic model represents a pore network as a randomly generated skeleton of cylindrical pipes, based on physical characteristics determined by the microscopic simulation. We present computational and computer science issues involved in the HPC implementation of the codes and in integrating them into a seamless simulation system. Issues such as portability, scalability, efficiency and extensibility of the final product are also discussed, as well as the numerical methods implemented at each scale. Example simulation results are presented.

  1. On-line Optimization-Based Simulators for Fractured and Non-fractured Reservoirs

    SciTech Connect

    Milind D. Deo

    2005-08-31

    Oil field development is a multi-million dollar business. Reservoir simulation is often used to guide the field management and development process. Reservoir characterization and geologic modeling tools have become increasingly sophisticated. As a result the geologic models produced are complex. Most reservoirs are fractured to a certain extent. The new geologic characterization methods are making it possible to map features such as faults and fractures, field-wide. Significant progress has been made in being able to predict properties of the faults and of the fractured zones. Traditionally, finite difference methods have been employed in discretizing the domains created by geologic means. For complex geometries, finite-element methods of discretization may be more suitable. Since reservoir simulation is a mature science, some of the advances in numerical methods (linear, nonlinear solvers and parallel computing) have not been fully realized in the implementation of most of the simulators. The purpose of this project was to address some of these issues. {sm_bullet} One of the goals of this project was to develop a series of finite-element simulators to handle problems of complex geometry, including systems containing faults and fractures. {sm_bullet} The idea was to incorporate the most modern computing tools; use of modular object-oriented computer languages, the most sophisticated linear and nonlinear solvers, parallel computing methods and good visualization tools. {sm_bullet} One of the tasks of the project was also to demonstrate the construction of fractures and faults in a reservoir using the available data and to assign properties to these features. {sm_bullet} Once the reservoir model is in place, it is desirable to find the operating conditions, which would provide the best reservoir performance. This can be accomplished by utilization optimization tools and coupling them with reservoir simulation. Optimization-based reservoir simulation was one of the

  2. Multigrid Methods for Fully Implicit Oil Reservoir Simulation

    NASA Technical Reports Server (NTRS)

    Molenaar, J.

    1996-01-01

    In this paper we consider the simultaneous flow of oil and water in reservoir rock. This displacement process is modeled by two basic equations: the material balance or continuity equations and the equation of motion (Darcy's law). For the numerical solution of this system of nonlinear partial differential equations there are two approaches: the fully implicit or simultaneous solution method and the sequential solution method. In the sequential solution method the system of partial differential equations is manipulated to give an elliptic pressure equation and a hyperbolic (or parabolic) saturation equation. In the IMPES approach the pressure equation is first solved, using values for the saturation from the previous time level. Next the saturations are updated by some explicit time stepping method; this implies that the method is only conditionally stable. For the numerical solution of the linear, elliptic pressure equation multigrid methods have become an accepted technique. On the other hand, the fully implicit method is unconditionally stable, but it has the disadvantage that in every time step a large system of nonlinear algebraic equations has to be solved. The most time-consuming part of any fully implicit reservoir simulator is the solution of this large system of equations. Usually this is done by Newton's method. The resulting systems of linear equations are then either solved by a direct method or by some conjugate gradient type method. In this paper we consider the possibility of applying multigrid methods for the iterative solution of the systems of nonlinear equations. There are two ways of using multigrid for this job: either we use a nonlinear multigrid method or we use a linear multigrid method to deal with the linear systems that arise in Newton's method. So far only a few authors have reported on the use of multigrid methods for fully implicit simulations. Two-level FAS algorithm is presented for the black-oil equations, and linear multigrid for

  3. The Interfacial Interaction Problem in Complex Multiple Porosity Fractured Reservoirs

    NASA Astrophysics Data System (ADS)

    Suarez-Arriaga, Mario-Cesar

    2003-04-01

    Many productive reservoirs (oil, gas, water, geothermal) are associated to natural fracturing. Fault zones and fractures act as open networks for fluid and energy flow from depth. Their petrophysical parameters are heterogeneous and randomly distributed, conforming extremely complex natural systems. Here, the simultaneous heat and mass flows are coupled to the deformation of thermoporoelastic rocks. The system's volume is divided into N interacting continua, each one occupying a region of space Vn wrapped by a surface Sn (n=1,N). The mass flow is represented by: ∂/∂t ∫ Vn ρf φdV + ∫ Sn F⃗M ṡ n⃗dS = ∫ Vn qMdV (3) Taking into account a non-isothermal process the coupled equation of energy is: ∂/∂t ∫ Vn [φρf hf + (1 - φ)ρrhr]dV + ∫ Sn F⃗E ṡ n⃗dS = ∫ Vn qEdV (4) Where t means time, φ is porosity, ρf, ρr are fluid and rock densities, F⃗M and F⃗E are total mass and energy flows, qM and qE are volumetric mass and energy extracted or injected into Vn, hf and hr are specific enthalpies for fluid and rock respectively. Rock deformation is coupled through the equation: ∇⃗ ṡ (ρf/μK ṡ ∇⃗pφ)Vn = φ (Dtρf + ρf/VφdVφ/dt)Vn (5) K is the absolute permeability tensor, μ means dynamic fluid viscosity, Dt is a total derivative, pφ is pore pressure and Vφ is the volume of pores in Vn. The N media interact with each other, every one has its own parameters and its own interporosity flow. Modelling these coupled phenomena requires to average highly contrasting physical properties, independently of the method used in the solution of equations. A lot of attention has been devoted to develop realistic numerical models to describe flows in reservoirs under exploitation. But to the best of our knowledge very little attention has been focused on the problem of interfacial interaction and averaging petrophysical parameters in multiple porosity reservoirs.

  4. Reservoir simulation improves implementation of Midway Sunset steamflood

    SciTech Connect

    Ellison, T.K.; Clayton, C.A.

    1995-12-31

    Thermal reservoir simulation was utilized to understand, make development recommendations, and project the performance of the Monarch C steamflood in a portion of Mobil`s South Midway Sunset field. The Monarch, a thick sequence of complex turbidite deposition, is characterized by extreme geological heterogeneity (lithofacies-controlled permeability and saturation variation, and mudstone barrier layers). Steamflood performance in the Monarch is related directly to the reservoir quality, and the path of steam flow is significantly influenced by the numerous laterally extensive mudstone barriers. The fine grain clay-bearing sediments were deposited on the anticline, distal from the source, whereas the coarser grain sediments, with little clay, were deposited on the more proximal syncline and steep dip areas. Consequently, steamflood performance improves relative to the crest since reservoir quality improves (including oil saturation), clay content decreases, and structure (dip) becomes more pronounced.

  5. High performance computing for domestic petroleum reservoir simulation

    SciTech Connect

    Zyvoloski, G.; Auer, L.; Dendy, J.

    1996-06-01

    This is the final report of a two-year, Laboratory-Directed Research and Development (LDRD) project at the Los Alamos National Laboratory. High-performance computing offers the prospect of greatly increasing the resolution at which petroleum reservoirs can be represented in simulation models. The increases in resolution can be achieved through large increases in computational speed and memory, if machine architecture and numerical methods for solution of the multiphase flow equations can be used to advantage. Perhaps more importantly, the increased speed and size of today`s computers make it possible to add physical processes to simulation codes that heretofore were too expensive in terms of computer time and memory to be practical. These factors combine to allow the development of new, more accurate methods for optimizing petroleum reservoir production.

  6. NFFLOW: A reservoir simulator incorporating explicit fractures (SPE 153890)

    SciTech Connect

    Boyle, E.J.; Sams, W.N.

    2012-01-01

    NFFLOW is a research code that quickly and inexpensively simulates flow in moderately fractured reservoirs. It explicitly recognizes fractures separately from rock matrix. In NFFLOW fracture flow is proportional to the pressure gradient along the fracture, and flow in the rock matrix is determined by Darcy’s Law. The two flow mechanisms are coupled through the pressure gradient between a fracture and its adjacent rock matrix. Presented is a promising change to NFFLOW that allows for flow across a rock matrix block.

  7. Reservoir characterization with sequential Gaussian simulation constrained by diffraction tomography

    SciTech Connect

    Lo, T.W.; Bermawi, A.

    1994-12-31

    A geostatistical approach for reservoir characterization that honors both surface seismic data and wireline data is described. It first computes a velocity profiles with seismic diffraction tomography, then, performs kriging with an external drift and sequential Gaussian simulation using the velocity profiles as soft data and the sonic logs as hard data. The product is a velocity profile with a resolution as high as that of the smoothed sonic logs, showing lateral velocity variations constrained by surface seismic data.

  8. Quantum simulations of physics problems

    NASA Astrophysics Data System (ADS)

    Somma, Rolando D.; Ortiz, Gerardo; Knill, Emanuel H.; Gubernatis, James

    2003-08-01

    If a large Quantum Computer (QC) existed today, what type of physical problems could we efficiently simulate on it that we could not simulate on a classical Turing machine? In this paper we argue that a QC could solve some relevant physical "questions" more efficiently. The existence of one-to-one mappings between different algebras of observables or between different Hilbert spaces allow us to represent and imitate any physical system by any other one (e.g., a bosonic system by a spin-1/2 system). We explain how these mappings can be performed showing quantum networks useful for the efficient evaluation of some physical properties, such as correlation functions and energy spectra.

  9. Reservoir simulation studies: Wairakei Geothermal Field, New Zealand. Final report

    SciTech Connect

    Pritchett, J.W.; Rice, L.F.; Garg, S.K.

    1980-01-01

    Numerical reservoir simulation techniques were used to perform a history-match of the Wairakei geothermal system in New Zealand. First, a one-dimensional (vertical) model was chosen; realistic stratigraphy was incorporated and the known production history was imposed. The effects of surface and deep recharge were included. Good matches were obtained, both for the reservoir pressure decline history and changes in average discharge enthalpy with time. Next, multidimensional effects were incorporated by treating with a two-dimensional vertical section. Again, good history matches were obtained, although computed late-time discharge enthalpies were slightly high. It is believed that this disparity arises from inherently three-dimensional effects. Predictive calculations using the two-dimensional model suggest that continued future production will cause little additional reservoir pressure drop, but that thermal degradation will occur. Finally, ground subsidence data at Wairakei was examined. It was concluded that traditional elastic pore-collapse models based on classical soil-mechanics concepts are inadequate to explain the observed surface deformation. It is speculated that the measured subsidence may be due to structural effects such as aseismic slippage of a buried reservoir boundary fault.

  10. Galerkin finite-element simulation of a geothermal reservoir

    USGS Publications Warehouse

    Mercer, J.W., Jr.; Pinder, G.F.

    1973-01-01

    The equations describing fluid flow and energy transport in a porous medium can be used to formulate a mathematical model capable of simulating the transient response of a hot-water geothermal reservoir. The resulting equations can be solved accurately and efficiently using a numerical scheme which combines the finite element approach with the Galerkin method of approximation. Application of this numerical model to the Wairakei geothermal field demonstrates that hot-water geothermal fields can be simulated using numerical techniques currently available and under development. ?? 1973.

  11. New Simulator for Non-Equilibrium Modeling of Hydrate Reservoirs

    NASA Astrophysics Data System (ADS)

    Kvamme, B.; Qorbani Nashaqi, K.; Jemai, K.; Vafaei, M.

    2014-12-01

    Due to Gibbs phase rule and combination of first and second law of thermodynamics, hydrate in nature cannot be in equilibrium since they come from different parent phases. In this system hydrate formation and dissociation is affected by local variables such as pressure, temperature and composition with mass and energy transport restrictions. Available simulators have attempted to model hydrate phase transition as an equilibrium reaction. Although those which treated the processes of formation and dissociation as kinetics used model of Kim and Bishnoi based on laboratory PVT experiment, and consequently hard to accept up scaling to real reservoirs condition. Additionally, they merely check equilibrium in terms of pressure and temperature projections and disregard thermodynamic requirements for equilibrium especially along axes of concentrations in phases. Non-equilibrium analysis of hydrate involves putting aside all the phase transitions which are not possible and use kinetic evaluation to measure phase transitions progress in each grid block for each time step. This procedure is Similar to geochemical reservoir simulators logic. As a result RetrasoCodeBright has been chosen as hydrate reservoir simulator and our work involves extension of this code. RetrasoCodeBright (RCB) is able to handle competing processes of formation and dissociation of hydrates as pseudo reactions at each node and each time step according to the temperature, pressure and concentration. Hydrates can therefore be implemented into the structure as pseudo minerals, with appropriate kinetic models. In order to implement competing nature of phase transition kinetics of hydrate formation, we use classical nucleation theory based on Kvamme et al. as a simplified model inside RCB and use advanced theories to fit parameters for the model (PFT). Hydrate formation and dissociation can directly be observed through porosity changes in the specific areas of the porous media. In this work which is in

  12. Generation of saturation functions for simulation models of carbonate reservoirs

    NASA Astrophysics Data System (ADS)

    Huang, Qingfeng

    A rock type is the unit of rock deposited under similar conditions, which went through similar diagenetic processes, producing analogous rock fabric, with distinct set of pore types, and pore throat size distribution, having specific range of porosity and permeability. Rock typing can generally be used as a guide to assign petrophysical characteristics to different zones for detailed reservoir characterization, modeling and simulation, which provide valid frames for reservoir development. It is often assumed that conventional rock types are capable of assigning multiphase flow characteristics, such as capillary pressure and relative permeability to the cells of dynamic simulation models. However, these conventional rock types, or static reservoir rock types (SRRT) fail to capture the actual variability of capillary pressure and relative permeability, due to lack of representation of wettability difference at different elevation above the free water level (FWL) in carbonate reservoirs, especially in the highly heterogeneous reservoirs. This should be resolved through dynamic reservoir rock types (DRRT), in which wettability effect is imposed on the SRRTs to generate saturation functions for simulation models. This research studies Ghedan's comprehensive DRRT model7, and proposes a modified Ghedan's model. First, the defined static rock types are sub-divided into sub-static rock types based on porosity frequency. Second, three curve-fitting programs are coded to generate the related saturation-height functions. These are the modified Ghedan-Okuyiga equation, Cuddy function and Power Law function. Developed from Ghedan-Okuyiga function113, the recommended modified Ghedan-Okuyiga function has been proposed with saturation and implicit porosity as a function of height above FWL in the transition zone. Third, each sub-static rock type is divided into a number of DRRTs by determining the capillary pressure and relative permeability curves in the oil zone from gas

  13. Scale Model Simulation of Enhanced Geothermal Reservoir Creation

    NASA Astrophysics Data System (ADS)

    Gutierrez, M.; Frash, L.; Hampton, J.

    2012-12-01

    Geothermal energy technology has successfully provided a means of generating stable base load electricity for many years. However, implementation has been spatially limited to limited availability of high quality traditional hydro-thermal resources possessing the combination of a shallow high heat flow anomaly and an aquifer with sufficient permeability and continuous fluid recharge. Enhanced Geothermal Systems (EGS) has been proposed as a potential solution to enable additional energy production from the non-conventional hydro-thermal resources. Hydraulic fracturing is considered the primary means of creating functional EGS reservoirs at sites where the permeability of the rock is too limited to allow cost effective heat recovery. EGS reservoir creation requires improved fracturing methodology, rheologically controllable fracturing fluids, and temperature hardened proppants. Although large fracture volumes (several cubic km) have been created in the field, circulating fluid through these full volumes and maintaining fracture volumes have proven difficult. Stimulation technology and methodology as used in the oil and gas industry for sedimentary formations are well developed; however, they have not sufficiently been demonstrated for EGS reservoir creation. Insufficient data and measurements under geothermal conditions make it difficult to directly translate experience from the oil and gas industries to EGS applications. To demonstrate the feasibility of EGS reservoir creation and subsequent geothermal energy production, and to improve the understanding of hydraulic and propping in EGS reservoirs, a heated true-triaxial load cell with a high pressure fluid injection system was developed to simulate an EGS system from stimulation to production. This apparatus is capable of loading a 30x30x30 cubic cm rock sample with independent principal stresses up to 13 MPa while simultaneously providing heating up to 180 degree C. Multiple orientated boreholes of 5 to 10 mm

  14. Geothermal reservoir engineering computer code comparison and validation using the GEONZ simulator program

    SciTech Connect

    Horne, R.N.; Ogbe, D.O.; Temeng, K.; Ramey, H.J. Jr.

    1980-11-14

    It was originally proposed to use the GEOTHERM geothermal simulator program to prepare solutions to the first five of the six Department of Energy computer code comparison sets. Valid solutions were anticipated in all of the five problems attempted, but it was expected that problems 3 and 4 might present some difficulties. A more recent version of the program, called GEONZ became available and was used successfully on problems 3 and 4. The new program, GEONZ, had additional capabilities that enabled it to handle both superheated steam and counterflows of steam and water. The choice of the GEONZ code is discussed, followed by an in-depth description of the solutions obtained for problems 1 through 5. The problem statements are included as Appendix A. The five problems are: 1-D Avdonin Solution, 1-D well test analysis, 2-D flow in fracture/block medium, 2-D phase system with drainage, and flow in a 2-D areal reservoir. (MHR)

  15. Characterization of a complex fluvial-deltaic reservoir for simulation

    SciTech Connect

    Begg, S.H.; Kay, A.; Gustason, G.

    1996-09-01

    Lithotype is shown to be the main geological control on the spatial distribution of reservoir flow properties. The distribution of lithotypes is in turn controlled by lithofacies. Using core data, assemblages of lithotypes were grouped into major facies associations (MFA`s). Log signatures were used to pick MFA`s in uncored wells to provide conditioning data for a stochastic description of their interwell distribution using the sequential indicator simulation (SIS) technique. Deterministic correlation was not enforced. Characteristic object-models of the distribution of lithotypes within each MFA were then generated and converted to fine-scale poro-perm models using core-lug data. These models were upscaled to yield characteristic poro-perm distributions at the reservoir simulation scale. Finally, the MFA model was used as a template to distribute the upscaled poro-perm within the simulation model. The model history-matched rapidly and accurately, even though the wells were conditioned only to the MFA`s and not to foot-by-foot data. Multiple MFA realizations were generated to provide some understanding of the uncertainty in thickness and other rock properties between wells.

  16. Simulating Oil Production from Fractured/Faulted Basement Reservoirs

    NASA Astrophysics Data System (ADS)

    Wang, H.; Forster, C.; Fu, Y.; Huang, C.; Yang, Y.; Deo, M.

    2006-12-01

    A fully-implicit, three-dimensional (3D), three-phase, discrete fault/fracture, black oil simulator provides new insight and understanding of oil production from reservoirs in fractured, low-permeability basement rocks. Results obtained with a controlled volume finite element (CVFE) method compare favorably to those obtained using both single- and dual-porosity finite difference methods (e.g., ECLIPSE). A regularized network of 30 orthogonal faults within a 1000 by 1000 by 200 feet model domain is used to compare the simulator results and to explore the implications of grid sensitivity. In this simple reservoir, cumulative oil recoveries over 900 days of production are similar for CVFE, single-porosity and dual-porosity approaches. CVFE is used to simulate a complex network of intersecting faults that mimic a more realistic basement reservoir with the same fault surface area and fault volume as the regularized network. Cumulative oil production at 900 days is about 3% lower than for the regularized network. The CVFE method provides a much improved ability to represent complex fracture/fault geometries and spatial variations in the internal properties of faults. CVFE simulations of the realistic network illustrate the possible consequences of uncertainty in knowing fracture/fault properties (e.g., porosity, permeability, thickness, dip orientation, connectivity and flow transmissibility). For example, introducing spatial variability in permeability within the fault planes (using spatially randomized patterns of 10, 100 and 1000 md), while retaining a constant geometric mean permeability of 100 md, yields enhanced oil production due to the high-permeability pathways. A 50:50 mix of 10 and 1000 md elements yields 36%OOIP while a 33:33:33 mix of 10, 100 and 1000 md yields 24%OOIP. These results are 26% and 14% greater, respectively, than that obtained for the uniform 100 md case (11%OOIP). This inherent variability, combined with uncertainty in knowing the detailed

  17. CO2/ brine substitution experiments at simulated reservoir conditions

    NASA Astrophysics Data System (ADS)

    Kummerow, Juliane; Spangenberg, Erik

    2015-04-01

    Capillary properties of rocks affect the mobility of fluids in a reservoir. Therefore, the understanding of the capillary pressure behaviour is essential to assess the long-term behaviour of CO2 reservoirs. Beyond this, a calibration of the petrophysical properties on water saturation of reservoir rocks at simulated in situ conditions is crucial for a proper interpretation of field monitoring data. We present a set-up, which allows for the combined measurements of capillary pressure, electric resistivity, and elastic wave velocities under controlled reservoir conditions (pconf = 400 bar, ppore = 180 bar, T = 65 ° C) at different brine-CO2 saturations. The capillary properties of the samples are measured using the micropore membrane technique. The sample is jacketed with a Viton tube (thickness = 4 mm) and placed between two current electrode endcaps, which as well contain pore fluid ports and ultrasonic P and S wave transducers. Between the sample and the lower endcap the hydrophilic semi-permeable micro-pore membrane (pore size = 100 nm) is integrated. It is embedded into filter papers to establish a good capillary contact and to protect the highly sensitive membrane against mechanical damage under load. Two high-precision syringe pumps are used to displace a quantified volume of brine by CO2 and determine the corresponding sample saturation. The fluid displacement induces a pressure gradient along the sample, which corresponds to the capillary pressure at a particular sample saturation. It is measured with a differential pressure sensor in the range between 0 - 0.2 MPa. Drainage and imbibition cycles are performed to provide information on the efficiency of capillary trapping and to get a calibration of the petrophysical parameters of the sample.

  18. Improving multi-objective reservoir operation optimization with sensitivity-informed problem decomposition

    NASA Astrophysics Data System (ADS)

    Chu, J. G.; Zhang, C.; Fu, G. T.; Li, Y.; Zhou, H. C.

    2015-04-01

    This study investigates the effectiveness of a sensitivity-informed method for multi-objective operation of reservoir systems, which uses global sensitivity analysis as a screening tool to reduce the computational demands. Sobol's method is used to screen insensitive decision variables and guide the formulation of the optimization problems with a significantly reduced number of decision variables. This sensitivity-informed problem decomposition dramatically reduces the computational demands required for attaining high quality approximations of optimal tradeoff relationships between conflicting design objectives. The search results obtained from the reduced complexity multi-objective reservoir operation problems are then used to pre-condition the full search of the original optimization problem. In two case studies, the Dahuofang reservoir and the inter-basin multi-reservoir system in Liaoning province, China, sensitivity analysis results show that reservoir performance is strongly controlled by a small proportion of decision variables. Sensitivity-informed problem decomposition and pre-conditioning are evaluated in their ability to improve the efficiency and effectiveness of multi-objective evolutionary optimization. Overall, this study illustrates the efficiency and effectiveness of the sensitivity-informed method and the use of global sensitivity analysis to inform problem decomposition when solving the complex multi-objective reservoir operation problems.

  19. Improved storage efficiency through geologic modeling and reservoir simulation

    SciTech Connect

    Ammer, J.R.; Mroz, T.H.; Covatch, G.L.

    1997-11-01

    The US Department of Energy (DOE), through partnerships with industry, is demonstrating the importance of geologic modeling and reservoir simulation for optimizing the development and operation of gas storage fields. The geologic modeling and reservoir simulation study for the Natural Fuel Gas Supply Corporation CRADA was completed in September 1995. The results of this study were presented at the 1995 Society of Petroleum Engineers` (SPE) Eastern Regional Meeting. Although there has been no field verification of the modeling results, the study has shown the potential advantages and cost savings opportunities of using horizontal wells for storage enhancement. The geologic modeling for the Equitrans` CRADA was completed in September 1995 and was also presented at the 1995 SPE Eastern Regional Meeting. The reservoir modeling of past field performance was completed in November 1996 and prediction runs are currently being made to investigate the potential of offering either a 10 day or 30 day peaking service in addition to the existing 110 day base load service. Initial results have shown that peaking services can be provided through remediation of well damage and by drilling either several new vertical wells or one new horizontal well. The geologic modeling for the Northern Indiana Public Service Company CRADA was completed in November 1996 with a horizontal well being completed in January 1997. Based on well test results, the well will significantly enhance gas deliverability from the field and will allow the utilization of gas from an area of the storage field that was not accessible from their existing vertical wells. Results are presented from these three case studies.

  20. Application of parallel computing techniques to a large-scale reservoir simulation

    SciTech Connect

    Zhang, Keni; Wu, Yu-Shu; Ding, Chris; Pruess, Karsten

    2001-02-01

    Even with the continual advances made in both computational algorithms and computer hardware used in reservoir modeling studies, large-scale simulation of fluid and heat flow in heterogeneous reservoirs remains a challenge. The problem commonly arises from intensive computational requirement for detailed modeling investigations of real-world reservoirs. This paper presents the application of a massive parallel-computing version of the TOUGH2 code developed for performing large-scale field simulations. As an application example, the parallelized TOUGH2 code is applied to develop a three-dimensional unsaturated-zone numerical model simulating flow of moisture, gas, and heat in the unsaturated zone of Yucca Mountain, Nevada, a potential repository for high-level radioactive waste. The modeling approach employs refined spatial discretization to represent the heterogeneous fractured tuffs of the system, using more than a million 3-D gridblocks. The problem of two-phase flow and heat transfer within the model domain leads to a total of 3,226,566 linear equations to be solved per Newton iteration. The simulation is conducted on a Cray T3E-900, a distributed-memory massively parallel computer. Simulation results indicate that the parallel computing technique, as implemented in the TOUGH2 code, is very efficient. The reliability and accuracy of the model results have been demonstrated by comparing them to those of small-scale (coarse-grid) models. These comparisons show that simulation results obtained with the refined grid provide more detailed predictions of the future flow conditions at the site, aiding in the assessment of proposed repository performance.

  1. Cooperative Learning in Reservoir Simulation Classes: Overcoming Disparate Entry Skills

    NASA Astrophysics Data System (ADS)

    Awang, Mariyamni

    2006-10-01

    Reservoir simulation is one of the core courses in the petroleum engineering curriculum and it requires knowledge and skills in three major disciplines, namely programming, numerical methods and reservoir engineering. However, there were often gaps in the students' readiness to undertake the course, even after completing the necessary requirements. The disparate levels of competency of the good and poor students made it difficult to target a certain level. Cooperative learning in the form of projects and peer teaching was designed to address the major concern of disparate entry skills, and at the same time the method used should also succeed in keeping students interest in class, developing communication skills and improving self-learning. Slower and weaker students were expected to benefit from being taught by good students, who were better prepared, and good students would gain deeper comprehension of the subject matter. From evaluations, the approach was considered successful since the overall passing rate was greater than 95% compared to previous years of around 70-80%. It had also succeeded in improving the learning environment in class. Future simulation classes will continue to use the cooperative approach with minor adjustments.

  2. Deriving multiple near-optimal solutions to deterministic reservoir operation problems

    NASA Astrophysics Data System (ADS)

    Liu, Pan; Cai, Ximing; Guo, Shenglian

    2011-08-01

    Even deterministic reservoir operation problems with a single objective function may have multiple near-optimal solutions (MNOS) whose objective values are equal or sufficiently close to the optimum. MNOS is valuable for practical reservoir operation decisions because having a set of alternatives from which to choose allows reservoir operators to explore multiple options whereas the traditional algorithm that produces a single optimum does not offer them this flexibility. This paper presents three methods: the near-shortest paths (NSP) method, the genetic algorithm (GA) method, and the Markov chain Monte Carlo (MCMC) method, to explore the MNOS. These methods, all of which require a long computation time, find MNOS using different approaches. To reduce the computation time, a narrower subspace, namely a near-optimal space (NOSP, described by the maximum and minimum bounds of MNOS) is derived. By confining the MNOS search within the NOSP, the computation time of the three methods is reduced. The proposed methods are validated with a test function before they are examined with case studies of both a single reservoir (the Three Gorges Reservoir in China) and a multireservoir system (the Qing River Cascade Reservoirs in China). It is found that MNOS exists for the deterministic reservoir operation problems. When comparing the three methods, the NSP method is unsuitable for large-scale problems but provides a benchmark to which solutions of small- and medium-scale problems can be compared. The GA method can produce some MNOS but is not very efficient in terms of the computation time. Finally, the MCMC method performs best in terms of goodness-of-fit to the benchmark and computation time, since it yields a wide variety of MNOS based on all retained intermediate results as potential MNOS. Two case studies demonstrate that the MNOS identified in this study are useful for real-world reservoir operation, such as the identification of important operation time periods and

  3. Foam flooding reservoir simulation algorithm improvement and application

    NASA Astrophysics Data System (ADS)

    Wang, Yining; Wu, Xiaodong; Wang, Ruihe; Lai, Fengpeng; Zhang, Hanhan

    2014-05-01

    As one of the important enhanced oil recovery (EOR) technologies, Foam flooding is being used more and more widely in the oil field development. In order to describe and predict foam flooding, experts at domestic and abroad have established a number of mathematical models of foam flooding (mechanism, empirical and semi-empirical models). Empirical models require less data and apply conveniently, but the accuracy is not enough. The aggregate equilibrium model can describe foam generation, burst and coalescence by mechanism studying, but it is very difficult to accurately describe. The research considers the effects of critical water saturation, critical concentration of foaming agent and critical oil saturation on the sealing ability of foam and considers the effect of oil saturation on the resistance factor for obtaining the gas phase relative permeability and the results were amended by laboratory test, so the accuracy rate is higher. Through the reservoir development concepts simulation and field practical application, the calculation is more accurate and higher.

  4. Homogenization and simulation for compositional flow in naturally fractured reservoirs

    NASA Astrophysics Data System (ADS)

    Chen, Zhangxin

    2007-02-01

    A dual porosity model of multidimensional, multicomponent, multiphase flow in naturally fractured reservoirs is derived by the mathematical theory of homogenization. A fully compositional model is considered where there are N chemical components, each of which may exist in any or all of the three phases: gas, oil, and water. Special attention is paid to developing a general approach to incorporating gravitational forces, pressure gradient effects, and effects of mass transfer between phases. In particular, general equations for the interactions between matrix and fracture systems are obtained under homogenization by a careful scaling of these effects. Using this dual porosity compositional model, numerical experiments are reported for the benchmark problems of the sixth comparative solution project organized by the society of petroleum engineers.

  5. Integration of geology and reservoir engineering to produce reservoir simulation model at Cabin Creek Field, Cedar Creek Anticline, Montana

    SciTech Connect

    Pieterson, R.; DiMarco, M.J.; Sodersten, S.S. )

    1996-01-01

    Because of its mature stage of development, a key aspect of continued economic development of the Cedar Creek Anticline (CCA), Montana (STOOIP > 2 billion barrels with over 360 MMstb produced) is the Identification of remaining pods of high S[sub o] within the original field boundaries. Present economic conditions make it essential to select drillsites with high probabilities of success and high prognoses flow rates in these remaking high S. area. Integration of a well-constrained geologic model and reservoir simulation pinpointed remaining pods of oil in a 3-m thick, subvertically fractured, dolomitic limestone reservoir of the Carboniferous Mission Canyon Formation in the Cabin Creek Field of the CCA This resulted in a successful high-flow-rate horizontal well (initial rate >800 BOPD) whose oil production was accurately predicted by a 3-D reservoir simulation. The model has 53,750 gridblocks each of which Is 60 by 60 m. The effect of the natural-fracture network was constrained with the k[sub v]/k[sub h] (vertical to horizontal permeability ratio). The simulation covered a 40-yr. production period. Gross production was Input as a constraint; oil and water rates were matched. Adjustments to absolute permeability, aquifer volume and relative water permeability were required to obtain a match between observed and simulated production rates. The model was fine tuned by matching the production of individual wells in areas with a high remaining S[sub o]. This project demonstrated that (1) interplay of geology and reservoir engineering provided a better reservoir model than could have been done individually, (2) simulation work identified horizontal drilling and recompletion candidates, with one successful horizontal well completed to date, and (3) use of the reservoir simulator for field-scale modeling In conjunction with a well-refined geologic synthesis can successfully pinpoint undeveloped reserves at CCA.

  6. Integration of geology and reservoir engineering to produce reservoir simulation model at Cabin Creek Field, Cedar Creek Anticline, Montana

    SciTech Connect

    Pieterson, R.; DiMarco, M.J.; Sodersten, S.S.

    1996-12-31

    Because of its mature stage of development, a key aspect of continued economic development of the Cedar Creek Anticline (CCA), Montana (STOOIP > 2 billion barrels with over 360 MMstb produced) is the Identification of remaining pods of high S{sub o} within the original field boundaries. Present economic conditions make it essential to select drillsites with high probabilities of success and high prognoses flow rates in these remaking high S. area. Integration of a well-constrained geologic model and reservoir simulation pinpointed remaining pods of oil in a 3-m thick, subvertically fractured, dolomitic limestone reservoir of the Carboniferous Mission Canyon Formation in the Cabin Creek Field of the CCA This resulted in a successful high-flow-rate horizontal well (initial rate >800 BOPD) whose oil production was accurately predicted by a 3-D reservoir simulation. The model has 53,750 gridblocks each of which Is 60 by 60 m. The effect of the natural-fracture network was constrained with the k{sub v}/k{sub h} (vertical to horizontal permeability ratio). The simulation covered a 40-yr. production period. Gross production was Input as a constraint; oil and water rates were matched. Adjustments to absolute permeability, aquifer volume and relative water permeability were required to obtain a match between observed and simulated production rates. The model was fine tuned by matching the production of individual wells in areas with a high remaining S{sub o}. This project demonstrated that (1) interplay of geology and reservoir engineering provided a better reservoir model than could have been done individually, (2) simulation work identified horizontal drilling and recompletion candidates, with one successful horizontal well completed to date, and (3) use of the reservoir simulator for field-scale modeling In conjunction with a well-refined geologic synthesis can successfully pinpoint undeveloped reserves at CCA.

  7. Development of a multiple application reservoir simulator for use on a vector computer

    SciTech Connect

    Kendall, R.P.; Morrell, G.O.; Peaceman, D.W.; Silliman, W.J.; Watts, J.W.

    1983-01-01

    This paper describes a new reservoir simulator possessing four distinctive features. It uses a recently developed formulation that is efficient in both black-oil and compositional modeling. It possesses a design that permits it to solve efficiently very large problems on a vector computer while being reasonably efficient on smaller problems and on scalar machines. It uses a novel approach to data management that permits writing most of the program as if it rar in core despite the fact that, in general, data is continually being temporarily stored on and retrieved from disk. Its well management package uses a flexible data structure and logic sequence that permits modeling of complex well and surface facility networks. The paper presents an overview of the simulator, discussing each of the features mentioned above. It also briefly discusses experience to date in using the new program and the conclusions reached based on this experience.

  8. Nonlinear Flow Numerical Simulation of an Ultra-Low Permeability Reservoir

    NASA Astrophysics Data System (ADS)

    Yu, Rong-Ze; Lei, Qun; Yang, Zheng-Ming; Bian, Ya-Nan

    2010-07-01

    A nonlinear flow mathematical model is established and the grid equation is deduced. A nonlinear flow reservoir numerical simulation program is compiled. The permeability loss coefficient is used to describe the permeability loss. A pilot calculation is made on the basis of actual field data, which reflects the reservoir development characteristics. The numerical simulation program based on nonlinear flow can anticipate the dynamic characteristics of the ultra-low permeability reservoir exploitation more exactly.

  9. Predicting the natural state of fractured carbonate reservoirs: An Andector Field, West Texas test of a 3-D RTM simulator

    SciTech Connect

    Tuncay, K.; Romer, S.; Ortoleva, P.; Hoak, T.; Sundberg, K.

    1998-12-31

    The power of the reaction, transport, mechanical (RTM) modeling approach is that it directly uses the laws of geochemistry and geophysics to extrapolate fracture and other characteristics from the borehole or surface to the reservoir interior. The objectives of this facet of the project were to refine and test the viability of the basin/reservoir forward modeling approach to address fractured reservoir in E and P problems. The study attempts to resolve the following issues: role of fracturing and timing on present day location and characteristics; clarifying the roles and interplay of flexure dynamics, changing rock rheological properties, fluid pressuring and tectonic/thermal histories on present day reservoir location and characteristics; and test the integrated RTM modeling/geological data approach on a carbonate reservoir. Sedimentary, thermal and tectonic data from Andector Field, West Texas, were used as input to the RTM basin/reservoir simulator to predict its preproduction state. The results were compared with data from producing reservoirs to test the RTM modeling approach. The effects of production on the state of the field are discussed in a companion report. The authors draw the following conclusions: RTM modeling is an important new tool in fractured reservoir E and P analysis; the strong coupling of RTM processes and the geometric and tensorial complexity of fluid flow and stresses require the type of fully coupled, 3-D RTM model for fracture analysis as pioneered in this project; flexure analysis cannot predict key aspects of fractured reservoir location and characteristics; fracture history over the lifetime of a basin is required to understand the timing of petroleum expulsion and migration and the retention properties of putative reservoirs.

  10. Chemical stimulation of gas condensate reservoirs: An experimental and simulation study

    NASA Astrophysics Data System (ADS)

    Kumar, Viren

    Well productivity in gas condensate reservoirs is reduced by condensate banking when the bottom hole flowing pressure drops below the dewpoint pressure. Several methods have been proposed to restore gas production rates after a decline due to condensate blocking. Gas injection, hydraulic fracturing, horizontal wells and methanol injection have been tried with limited success. These methods of well stimulation either offer only temporary productivity restoration or are applicable only in some situations. Wettability alteration of the rock in the near well bore region is an economic and efficient method for the enhancement of gas-well deliverability. Altering the wettability of porous media from strongly water-wet or oil-wet to intermediate-wet decreases the residual liquid saturations and results in an increase in the relative permeability to gas. Such treatments also increase the mobility and recovery of condensate from the reservoir. This study validates the above hypothesis and provides a simple and cost-efficient solution to the condensate blocking problem. Screening studies were carried out to identify the chemicals based on structure, solubility and reactivity at reservoir temperature and pressure. Experiments were performed to evaluate these chemicals to improve gas and condensate relative permeabilities. The improvement in relative permeability after chemical treatment was quantified by performing high pressure and high temperature coreflood experiments in Berea sandstone, Texas Cream limestone and reservoir cores using synthetic gas mixtures at reservoir conditions. Experiments were done at high flow rates and for long time periods to evaluate the durability of the treatment. Single well simulation studies were conducted to demonstrate the performance of the chemical treatment in the field. The experimental relative permeability data was modeled using a trapping number dependent relative permeability model and incorporated in the simulations. Effect of

  11. Advanced Techniques for Reservoir Simulation and Modeling of Non-Conventional Wells

    SciTech Connect

    Durlofsky, Louis J.

    2000-08-28

    This project targets the development of (1) advanced reservoir simulation techniques for modeling non-conventional wells; (2) improved techniques for computing well productivity (for use in reservoir engineering calculations) and well index (for use in simulation models), including the effects of wellbore flow; and (3) accurate approaches to account for heterogeneity in the near-well region.

  12. Simulated Three-Point Problems.

    ERIC Educational Resources Information Center

    Leyden, Michael B.

    1979-01-01

    The concept of sloping bedrock strata is portrayed by simple construction of a cardboard model. By use of wires and graph paper, students simulate the drilling of wells and use standard mathematical operations to determine strike and dip of the model stratum. (RE)

  13. Numerical simulation of water injection into vapor-dominated reservoirs

    SciTech Connect

    Pruess, K.

    1995-01-01

    Water injection into vapor-dominated reservoirs is a means of condensate disposal, as well as a reservoir management tool for enhancing energy recovery and reservoir life. We review different approaches to modeling the complex fluid and heat flow processes during injection into vapor-dominated systems. Vapor pressure lowering, grid orientation effects, and physical dispersion of injection plumes from reservoir heterogeneity are important considerations for a realistic modeling of injection effects. An example of detailed three-dimensional modeling of injection experiments at The Geysers is given.

  14. On an iterative ensemble smoother and its application to a reservoir facies estimation problem

    NASA Astrophysics Data System (ADS)

    Luo, Xiaodong; Chen, Yan; Valestrand, Randi; Stordal, Andreas; Lorentzen, Rolf; Nævdal, Geir

    2014-05-01

    For data assimilation problems there are different ways in utilizing the available observations. While certain data assimilation algorithms, for instance, the ensemble Kalman filter (EnKF, see, for examples, Aanonsen et al., 2009; Evensen, 2006) assimilate the observations sequentially in time, other data assimilation algorithms may instead collect the observations at different time instants and assimilate them simultaneously. In general such algorithms can be classified as smoothers. In this aspect, the ensemble smoother (ES, see, for example, Evensen and van Leeuwen, 2000) can be considered as an smoother counterpart of the EnKF. The EnKF has been widely used for reservoir data assimilation (history matching) problems since its introduction to the community of petroleum engineering (Nævdal et al., 2002). The applications of the ES to reservoir data assimilation problems are also investigated recently (see, for example, Skjervheim and Evensen, 2011). Compared to the EnKF, the ES has certain technical advantages, including, for instance, avoiding the restarts associated with each update step in the EnKF and also having fewer variables to update, which may result in a significant reduction in simulation time, while providing similar assimilation results to those obtained by the EnKF (Skjervheim and Evensen, 2011). To further improve the performance of the ES, some iterative ensemble smoothers are suggested in the literature, in which the iterations are carried out in the forms of certain iterative optimization algorithms, e.g., the Gaussian-Newton (Chen and Oliver, 2012) or the Levenberg-Marquardt method (Chen and Oliver, 2013; Emerick and Reynolds, 2012), or in the context of adaptive Gaussian mixture (AGM, see Stordal and Lorentzen, 2013). In Emerick and Reynolds (2012) the iteration formula is derived based on the idea that, for linear observations, the final results of the iterative ES should be equal to the estimate of the EnKF. In Chen and Oliver (2013), the

  15. Geothermal reservoir simulation to enhance confidence in predictions for nuclear waste disposal

    SciTech Connect

    Kneafsey, Timothy J.; Pruess, Karsten; O'Sullivan, Michael J.; Bodvarsson, Gudmundur S.

    2002-06-15

    Numerical simulation of geothermal reservoirs is useful and necessary in understanding and evaluating reservoir structure and behavior, designing field development, and predicting performance. Models vary in complexity depending on processes considered, heterogeneity, data availability, and study objectives. They are evaluated using computer codes written and tested to study single and multiphase flow and transport under nonisothermal conditions. Many flow and heat transfer processes modeled in geothermal reservoirs are expected to occur in anthropogenic thermal (AT) systems created by geologic disposal of heat-generating nuclear waste. We examine and compare geothermal systems and the AT system expected at Yucca Mountain, Nevada, and their modeling. Time frames and spatial scales are similar in both systems, but increased precision is necessary for modeling the AT system, because flow through specific repository locations will affect long-term ability radionuclide retention. Geothermal modeling experience has generated a methodology, used in the AT modeling for Yucca Mountain, yielding good predictive results if sufficient reliable data are available and an experienced modeler is involved. Codes used in geothermal and AT modeling have been tested extensively and successfully on a variety of analytical and laboratory problems.

  16. Long-term Reservoir Routing Simulations Using Data-Driven Approaches

    NASA Astrophysics Data System (ADS)

    Ashouri, H.; Chowdhary, H.; Chinnayakanahalli, K.; Dodov, B.

    2015-12-01

    Flood is a highly complex natural hazard that accounts for major losses to human societies worldwide. Dams built with the aim of mitigating the flood risk significantly modify river flow regimes but unavailability and/or inaccessibility of proper information about reservoir operational rules impose a big hurdle to global flood modeling. This is specifically critical for flood-prone regions where lack of proper representation of reservoir operation can lead to significant under- or overestimation of the flood magnitude, risk, and losses. With the availability of longer in-situ observational data records, as well as advancements in satellite altimetry techniques for measuring reservoir levels, operational rules can be indirectly deduced. In this study, the observed reservoir levels as well as the historical and forecast time series of inflows are incorporated into a stochastic autoregressive moving average statistical modeling scheme to simulate the releases from the dam at each time step. The resulting operational rule curve is used in a reservoir simulation model to simulate the outflows from the reservoirs. The efficiency of the model is examined for three case studies in the United States, including John Martin Reservoir (CO), Coralville Lake (IA, and specifically for the devastating 2008 flood in the state), and Boca Reservoir (CA). Statistical measures are derived and tested to evaluate the accuracy of the simulated hydrographs against USGS streamflow gauge observations. The results prove the capability of the developed model in simulating reasonably accurate outflows from dams and will be presented at the meeting.

  17. Moomba Lower Daralingie Beds (LDB) gas storage project: Reservoir management using a novel numerical simulation technique

    SciTech Connect

    Jamal, F.G.

    1994-12-31

    Engineers managing underground gas storage projects are often faced with challenges involving gas migration, inventory variance, gas quality and inventory-pressures. This paper discusses a unique underground gas storage project where sales gas and ethane are stored in two different but communicating regions of the same reservoir. A commercially available reservoir simulator was used to model the fluid flow behavior in this reservoir, hence, providing a tool for better management and use of the existing gas storage facilities.

  18. Nonlinear filtering in oil/gas reservoir simulation: filter design

    SciTech Connect

    Arnold, E.M.; Voss, D.A.; Mayer, D.W.

    1980-10-01

    In order to provide an additional mode of utility to the USGS reservoir model VARGOW, a nonlinear filter was designed and incorporated into the system. As a result, optimal (in the least squares sense) estimates of reservoir pressure, liquid mass, and gas cap plus free gas mass are obtained from an input of reservoir initial condition estimates and pressure history. These optimal estimates are provided continuously for each time after the initial time, and the input pressure history is allowed to be corrupted by measurement error. Preliminary testing of the VARGOW filter was begun and the results show promise. Synthetic data which could be readily manipulated during testing was used in tracking tests. The results were positive when the initial estimates of the reservoir initial conditions were reasonably close. Further testing is necessary to investigate the filter performance with real reservoir data.

  19. Fishery management problems and possibilities on large southeastern reservoirs

    USGS Publications Warehouse

    Parsons, John W.

    1958-01-01

    In recognition of these problems, the development and application of sound management procedures may be accomplished by controlling species composition and availability of fish through water level control, timber clearing, application of selective toxicants, commercial fishing, introduction of new fish species, and the management of tailwaters and tributaries. Extended research and interagency cooperation are necessary to properly develop and apply sound management. Promotion of angling and regulation of the fishery may best be realized by providing adequate fishing facilities, elimination of certain restricted areas, and in some cases revision of laws and regulations. Biologists must not only meet the present demands for improved sport fishing but must also balance fish yield by increasing the use of food fish by sport and commercial operations.

  20. The Simulation of Inflow Discharge and Suspended Sediment Transport Rate for a Reservoir

    NASA Astrophysics Data System (ADS)

    Wu, Ching-Hsien; Chen, Ching-Nuo; Tsai, Chih-Heng; Tsai, Chang-Tai

    2010-05-01

    The major functions of a reservoir include flood-protection, public water-supply, irrigation, hydropower and tourism. Consequently, these functions can provide great contributions for economic development. Therefore, important issues associated with reservoir watersheds such as soil erosion and deposition must be carefully studied in order to enhance watershed management. Accurate and timely estimation of peak flow discharges into a reservoir is very crucial for flood protection strategies and the general safety of the reservoir. In this study, GIS is applied to a physiographical soil erosion-deposition model, using rainfall data as the primary input, to simulate both stream flow hydrographs and sediment transport into a reservoir. In this study, the Shihmen Reservoir watershed is used as an illustrative example. The Shihmen Reservoir, a multiple-purpose reservoir for irrigation, hydropower, public water-supply, flood-protection, and tourism, is located on the mid-stream reach of the Tahan River. High-concentration sediment-containing flood flow events during storm periods are primarily responsible for the turbidity and sedimentation in the Shihmen Reservoir. Therefore, if incoming stream flow and sediment concentration hydrographs can be determined promptly during storm periods, high-concentration turbid water can be diverted from the reservoir and low-concentration water, during the recession, can be stored. Consequently, the deposition of the reservoir can be reduced extending the life of the reservoir. Therefore, the purpose of this project is to establish a physiographic soil erosion-deposition model for the Shihmen Reservoir watershed. Using rainfall data as input, we can simulate the hydrographs of both stream flow and sediment concentration entering the reservoir. The results obtained from our model can be used as a reference to aid in the operation of the Shihmen Reservoir concerning deposition prevention. The numerical studies show that the peak flow and

  1. An Exact Solution to the Draining Reservoir Problem of the Incompressible and Non-Viscous Liquid

    ERIC Educational Resources Information Center

    Hong, Seok-In

    2009-01-01

    The exact expressions for the drain time and the height, velocity and acceleration of the free surface are found for the draining reservoir problem of the incompressible and non-viscous liquid. Contrary to the conventional approximate results, they correctly describe the initial time dependence of the liquid velocity and acceleration. Torricelli's…

  2. Universal approximators for multi-objective direct policy search in water reservoir management problems: a comparative analysis

    NASA Astrophysics Data System (ADS)

    Giuliani, Matteo; Mason, Emanuele; Castelletti, Andrea; Pianosi, Francesca

    2014-05-01

    The optimal operation of water resources systems is a wide and challenging problem due to non-linearities in the model and the objectives, high dimensional state-control space, and strong uncertainties in the hydroclimatic regimes. The application of classical optimization techniques (e.g., SDP, Q-learning, gradient descent-based algorithms) is strongly limited by the dimensionality of the system and by the presence of multiple, conflicting objectives. This study presents a novel approach which combines Direct Policy Search (DPS) and Multi-Objective Evolutionary Algorithms (MOEAs) to solve high-dimensional state and control space problems involving multiple objectives. DPS, also known as parameterization-simulation-optimization in the water resources literature, is a simulation-based approach where the reservoir operating policy is first parameterized within a given family of functions and, then, the parameters optimized with respect to the objectives of the management problem. The selection of a suitable class of functions to which the operating policy belong to is a key step, as it might restrict the search for the optimal policy to a subspace of the decision space that does not include the optimal solution. In the water reservoir literature, a number of classes have been proposed. However, many of these rules are based largely on empirical or experimental successes and they were designed mostly via simulation and for single-purpose reservoirs. In a multi-objective context similar rules can not easily inferred from the experience and the use of universal function approximators is generally preferred. In this work, we comparatively analyze two among the most common universal approximators: artificial neural networks (ANN) and radial basis functions (RBF) under different problem settings to estimate their scalability and flexibility in dealing with more and more complex problems. The multi-purpose HoaBinh water reservoir in Vietnam, accounting for hydropower

  3. G-2 and G-3 reservoirs, Delta South field, Nigeria - 2. Simulation of water injection

    SciTech Connect

    Thakur, G.C.; Stanat, P.L.; Aruna, M.; Ajayi, S.A.; Poston, S.

    1982-01-01

    A description is given of a two-dimensional, three-phase, black-oil simulation of the G-2 and G-3 reservoirs in the Delta South field offshore Nigeria. The purpose of these studies was to investigate, from an engineering standpoint, various operating schemes for optimizing the oil recovery from each of these highly gravity-segregated reservoirs. 4 refs.

  4. Coupling of geochemical and multiphase flow processes for validation of the MUFITS reservoir simulator against TOUGHREACT

    NASA Astrophysics Data System (ADS)

    De Lucia, Marco; Kempka, Thomas; Afanasyev, Andrey; Melnik, Oleg; Kühn, Michael

    2016-04-01

    Coupled reactive transport simulations, especially in heterogeneous settings considering multiphase flow, are extremely time consuming and suffer from significant numerical issues compared to purely hydrodynamic simulations. This represents a major hurdle in the assessment of geological subsurface utilization, since it constrains the practical application of reactive transport modelling to coarse spatial discretization or oversimplified geological settings. In order to overcome such limitations, De Lucia et al. [1] developed and validated a one-way coupling approach between geochemistry and hydrodynamics, which is particularly well suited for CO2 storage simulations, while being of general validity. In the present study, the models used for the validation of the one-way coupling approach introduced by De Lucia et al. (2015), and originally performed with the TOUGHREACT simulator, are transferred to and benchmarked against the multiphase reservoir simulator MUFITS [2]. The geological model is loosely inspired by an existing CO2 storage site. Its grid comprises 2,950 elements enclosed in a single layer, but reflecting a realistic three-dimensional anticline geometry. For the purpose of this comparison, homogeneous and heterogeneous scenarios in terms of porosity and permeability were investigated. In both cases, the results of the MUFITS simulator are in excellent agreement with those produced with the fully-coupled TOUGHREACT simulator, while profiting from significantly higher computational performance. This study demonstrates how a computationally efficient simulator such as MUFITS can be successfully included in a coupled process simulation framework, and also suggests ameliorations and specific strategies for the coupling of chemical processes with hydrodynamics and heat transport, aiming at tackling geoscientific problems beyond the storage of CO2. References [1] De Lucia, M., Kempka, T., and Kühn, M. A coupling alternative to reactive transport simulations

  5. Using microstructure observations to quantify fracture properties and improve reservoir simulations. Final report, September 1998

    SciTech Connect

    Laubach, S.E.; Marrett, R.; Rossen, W.; Olson, J.; Lake, L.; Ortega, O.; Gu, Y.; Reed, R.

    1999-01-01

    The research for this project provides new technology to understand and successfully characterize, predict, and simulate reservoir-scale fractures. Such fractures have worldwide importance because of their influence on successful extraction of resources. The scope of this project includes creation and testing of new methods to measure, interpret, and simulate reservoir fractures that overcome the challenge of inadequate sampling. The key to these methods is the use of microstructures as guides to the attributes of the large fractures that control reservoir behavior. One accomplishment of the project research is a demonstration that these microstructures can be reliably and inexpensively sampled. Specific goals of this project were to: create and test new methods of measuring attributes of reservoir-scale fractures, particularly as fluid conduits, and test the methods on samples from reservoirs; extrapolate structural attributes to the reservoir scale through rigorous mathematical techniques and help build accurate and useful 3-D models of the interwell region; and design new ways to incorporate geological and geophysical information into reservoir simulation and verify the accuracy by comparison with production data. New analytical methods developed in the project are leading to a more realistic characterization of fractured reservoir rocks. Testing diagnostic and predictive approaches was an integral part of the research, and several tests were successfully completed.

  6. Simulations of Flow Circulations and Atrazine Concentrations in a Midwest U.S. Reservoir

    NASA Astrophysics Data System (ADS)

    Zhao, Xianggui; Gu, Roy R.; Guo, Chuling; Wang, Kui; Li, Shijie

    Atrazine is the most commonly used herbicide in the spring for pre-emergent weed control in the corn cropping area in the Midwestern United States. A frequent high level of herbicide concentrations in reservoirs is a great concern for public health and aquatic ecosystems. In this study, a two-dimensional hydrodynamics and toxic contaminant transport model was applied to Saylorville Reservoir, Iowa, USA. The model simulates physical, chemical, and biological processes and predicts unsteady vertical and longitudinal distributions of a toxic chemical. Model results were validated by measured temperatures and atrazine concentrations. Simulated flow velocities, water temperatures, and chemical concentrations demonstrated that the spatial variation of atrazine concentrations was largely affected by seasonal flow circulation patterns in the reservoir. In particular, the simulated fate and transport of atrazine showed the effect of flow circulation on spatial distribution of atrazine during summer months as the river flow formed an underflow within the reservoir and resulted in greater concentrations near the surface of the reservoir. Atrazine concentrations in the reservoir peaked around the end of May and early June. A thorough understanding of the fate and transport of atrazine in the reservoir can assist in developing operation and pollution prevention strategies with respect to timing, amount, and depth of withdrawal. The responses of atrazine transport to various boundary conditions provide useful information in assessing environmental impact of alternative upstream watershed management practices on the quality of reservoir water.

  7. Simulating the effects of adsorption and capillary forces in geothermal reservoirs

    SciTech Connect

    Sta. Maria, Roman B.; Pingol, Alponso S.

    1996-01-24

    Until recently, geothermal reservoir simulators use flat interface thermodynamics to determine the thermodynamic state of the reservoir. Development of new simulators and the modification of existing ones has now incorporated the physics of curved interface thermodynamics. These simulators account for the effects of sorption and capillary forces. The simulators GSS and TETRAD were used to simulate the performance of a hypothetical vapordominated geothermal reservoir. GSS is a simulator specifically developed to account for adsorption by using adsorption isotherms. On the other hand, TETRAD is a commercial simulator that was modified to account for vapor pressure lowering by using capillary pressure relations. GSS and TETRAD yielded similar results. Thus, the two formulations being used to account for curved interface thermodynamics are practically equivalent. Areas for improvement of both GSS and TETRAD were identified. The hysteresis and temperature dependence of sorption and capillary properties are issues that are needed to be addressed.

  8. Analysis of formation pressure test results in the Mount Elbert methane hydrate reservoir through numerical simulation

    USGS Publications Warehouse

    Kurihara, M.; Sato, A.; Funatsu, K.; Ouchi, H.; Masuda, Y.; Narita, H.; Collett, T.S.

    2011-01-01

    Targeting the methane hydrate (MH) bearing units C and D at the Mount Elbert prospect on the Alaska North Slope, four MDT (Modular Dynamic Formation Tester) tests were conducted in February 2007. The C2 MDT test was selected for history matching simulation in the MH Simulator Code Comparison Study. Through history matching simulation, the physical and chemical properties of the unit C were adjusted, which suggested the most likely reservoir properties of this unit. Based on these properties thus tuned, the numerical models replicating "Mount Elbert C2 zone like reservoir" "PBU L-Pad like reservoir" and "PBU L-Pad down dip like reservoir" were constructed. The long term production performances of wells in these reservoirs were then forecasted assuming the MH dissociation and production by the methods of depressurization, combination of depressurization and wellbore heating, and hot water huff and puff. The predicted cumulative gas production ranges from 2.16??106m3/well to 8.22??108m3/well depending mainly on the initial temperature of the reservoir and on the production method.This paper describes the details of modeling and history matching simulation. This paper also presents the results of the examinations on the effects of reservoir properties on MH dissociation and production performances under the application of the depressurization and thermal methods. ?? 2010 Elsevier Ltd.

  9. Solid hydrocarbon: a migration-of-fines problem in carbonate reservoirs

    SciTech Connect

    Lomando, A.J.

    1986-05-01

    The most familiar example of a migration-of-fines problem is authigenic kaolinite, which can detach, migrate through a pore system, and bridge pore throats, thus reducing permeability. under certain conditions, a similar problem is caused by solid hydrocarbon, independent of a mode of origin, which has precipitated in carbonate pore systems. Cores from several reservoirs in the Lower Cretaceous of east Texas were used as the data base in this study. Three morphotypes of solid hydrocarbon have been identified from thin-section and scanning electron microscope observations: droplets, peanut brittle, and carpets. Droplets are small, individual, rounded particles scattered on pore walls. Peanut brittle ranges from a continuous to discontinuous thin coating with random rounded lumps that probably have droplet precursors. Carpets are thick, continuous coatings and, at the extreme, can effectively occlude whole pores. Initially, solid hydrocarbon reduces permeability without necessarily decreasing porosity significantly. Likewise, solid hydrocarbon cannot be detected directly from wireline logs. Acidizing to enhance communication to the well bore is a common completion procedure in limestone and calcareous sandstone reservoirs. In reservoirs containing solid hydrocarbon, acid etches the substrate and releases solid hydrocarbon, which migrates in the pore system and bridges pore throats. Differential well-bore pressure also may cause solid hydrocarbon to migrate. Therefore, wettability, which controls hydrocarbon adhesion to the pore walls, and the dominant morphotype are important factors in the extent of reservoir damage.

  10. Vortex formation in coalescence of droplets with a reservoir using molecular dynamics simulations.

    PubMed

    Taherian, Fereshte; Marcon, Valentina; Bonaccurso, Elmar; van der Vegt, Nico F A

    2016-10-01

    The flow patterns generated by the coalescence of aqueous ethanol droplets with a water reservoir are investigated using molecular dynamics simulations. The influence of surface tension gradient, which leads to the spreading of the droplet along the liquid-vapor interface of the reservoir, is studied by changing the ethanol concentration of the droplet. The internal circulation (vortex strength) of the droplet and the reservoir are analyzed separately. Simulation results reveal the formation of swirling flows within the droplet at early times when the radius of the coalescence neck due to the capillary forces increases rapidly with time. The vortex strength is found to be higher at lower concentrations of ethanol (higher liquid-vapor surface tension of the droplet), where the driving force for the contact line movement (capillary force) is stronger. The circulation diminishes by moving the center of mass of the droplet toward the reservoir. The lower surface tension of the droplet compared to the reservoir leads to surface tension gradient driven flow, which transports the droplet molecules along the liquid-vapor interface of the reservoir. Such a flow motion results in the generation of convective flows in the underlying water, which forms swirling flows within the reservoir. Therefore, the vortex strength of the reservoir is higher at higher ethanol concentrations of the droplet. The reservoir circulation decays to zero as soon as the ethanol concentration becomes homogeneous along the interface of the pool. The time evolution of circulation within the droplet and the reservoir are correlated with the center of mass motion of the droplet toward the surface, the time variation of the precursor film radius and the dynamic surface tension of the reservoir. PMID:27388133

  11. Integration of seismic methods with reservoir simulation, Pikes Peak heavy oil field, Saskatchewan

    NASA Astrophysics Data System (ADS)

    Zou, Ying

    The Pikes Peak heavy oil field has been operated by Husky Energy Ltd since 1981. Steam injection has been successfully employed to increase production. Efforts in geophysics and reservoir engineering have been made to improve interpretations in the mapping of reservoir conditions. This dissertation developed tools and a working flow for integrating the analysis of time-lapse seismic surveys with reservoir simulation, and applied them to the Pikes Peak field. Two time-lapse 2D seismic lines acquired in February 1991 and March 2000 in the eastern part of the field were carefully processed to produce wavelet and structure matched final sections. Reservoir simulation based on the field reservoir production history was carried out. It provided independent complementary information for the time-lapse seismic analysis. A rock physics procedure based on Gassmann's equation and Batzle and Wang's empirical relationship successfully linked the reservoir engineering to the seismic method. Based on the resultant seismic models, synthetic seismic sections were generated as the analogy of field seismic sections. The integrated interpretation for the Pikes Peak reservoir drew the following conclusions: The areas with a gas saturation difference, between two compared time steps, have seismic differences. Thicker gas zones correspond with large reflectivity changes on the top of the reservoir and larger traveltime delays in the seismic section. The thin gas zones only induce large reflectivity changes on the top of the reservoir, and do not have large time delays below the reservoir zone. High temperature regions also correlate with areas having large seismic energy differences. High temperature with thick gas (steam and methane) zones may be evidence for steam existence. The seismic differences at locations far from the production zone are due to the lower pressure that causes solution gas to evolve from the oil. Pressure changes propagate much faster (˜20 m in one month) than

  12. Student Ecosystems Problem Solving Using Computer Simulation.

    ERIC Educational Resources Information Center

    Howse, Melissa A.

    The purpose of this study was to determine the procedural knowledge brought to, and created within, a pond ecology simulation by students. Environmental Decision Making (EDM) is an ecosystems modeling tool that allows users to pose their own problems and seek satisfying solutions. Of specific interest was the performance of biology majors who had…

  13. Simulation and resistivity modeling of a geothermal reservoir with waters of different salinity

    SciTech Connect

    Pruess, K.; Wilt, M.; Bodvarsson, G.S.; Goldstein, N.E.

    1982-10-01

    Apparent resistivities measured by means of repetitive dipole-dipole surveys show significant changes within the Cerro Prieto reservoir. The changes are attributed to production and natural recharge. To better understand the observed geophysical phenomena a simple reservoir simulation study combined with the appropriate DC resistivity calculations to determine the expected magnitude of apparent resistivity change. We consider production from a liquid-dominated reservoir with dimensions and parameters of the Cerro Prieto A reservoir and assume lateral and vertical recharge of colder and less saline waters. Based on rather schematic one- and two-dimensional reservoir simulations, we calculate changes in formation resistivity which we then transform into changes in apparent resistivity that would be observed at the surface. Simulated changes in apparent resistivities over the production zone show increases of 10 to 20% over a 3 year period at the current rate of fluid extraction. Changes of this magnitude are not only within our ability to discern using proper field techniques, but are consistent in magnitude with some of the observed effects. However, the patterns of apparent resistivity changes in the simulated dipole-dipole pseudosection only partially resemble the observed field data. This is explained by the fact that the actual fluid recharge into the A reservoir is more complicated than assumed in our simple, schematic recharge models.

  14. Optimization of coalbed-methane-reservoir exploration and development strategies through integration of simulation and economics

    SciTech Connect

    Clarkson, C.R.; McGovern, J.M.

    2005-12-15

    The unique properties and complex characteristics of coalbed methane (CBM) reservoirs, and their consequent operating strategies, call for an integrated approach to be used to explore for and develop coal plays and prospects economically. An integrated approach involves the use of sophisticated reservoir, wellbore, and facilities modeling combined with economics and decision-making criteria. A new CBM prospecting tool (CPT) was generated by combining single-well (multilayered) reservoir simulators with a gridded reservoir model, Monte Carlo (MC) simulation, and economic modules. The multilayered reservoir model is divided into pods, representing relatively uniform reservoir properties, and a 'type well' is created for each pod. At every MC iteration, type-well forecasts are generated for the pods and are coupled with economic modules. A set of decision criteria contingent upon economic outcomes and reservoir characteristics is used to advance prospect exploration from the initial exploration well to the pilot and development stages. A novel approach has been used to determine the optimal well spacing should prospect development be contemplated. CPT model outcomes include a distribution of after-tax net present value (ATNPV), mean ATNPV (expected value), chance of economic success (Pe), distribution of type-well and pod gas and water production, reserves, peak gas volume, and capita. Example application of CPT to a hypothetical prospect is provided. An integrated approach also has been used to assist with production optimization of developed reservoirs. For example, an infill-well locating tool (ILT) has been constructed to provide a quick-look evaluation of infill locations in a developed reservoir. An application of ILT to a CBM reservoir is provided.

  15. Comparisons of Simulated Hydrodynamics and Water Quality for Projected Demands in 2046, Pueblo Reservoir, Southeastern Colorado

    USGS Publications Warehouse

    Ortiz, Roderick F.; Galloway, Joel M.; Miller, Lisa D.; Mau, David P.

    2008-01-01

    Pueblo Reservoir is one of southeastern Colorado's most valuable water resources. The reservoir provides irrigation, municipal, and industrial water to various entities throughout the region. The reservoir also provides flood control, recreational activities, sport fishing, and wildlife enhancement to the region. The Bureau of Reclamation is working to meet its goal to issue a Final Environmental Impact Statement (EIS) on the Southern Delivery System project (SDS). SDS is a regional water-delivery project that has been proposed to provide a safe, reliable, and sustainable water supply through the foreseeable future (2046) for Colorado Springs, Fountain, Security, and Pueblo West. Discussions with the Bureau of Reclamation and the U.S. Geological Survey led to a cooperative agreement to simulate the hydrodynamics and water quality of Pueblo Reservoir. This work has been completed and described in a previously published report, U.S. Geological Survey Scientific Investigations Report 2008-5056. Additionally, there was a need to make comparisons of simulated hydrodynamics and water quality for projected demands associated with the various EIS alternatives and plans by Pueblo West to discharge treated water into the reservoir. Plans by Pueblo West are fully independent of the SDS project. This report compares simulated hydrodynamics and water quality for projected demands in Pueblo Reservoir resulting from changes in inflow and water quality entering the reservoir, and from changes to withdrawals from the reservoir as projected for the year 2046. Four of the seven EIS alternatives were selected for scenario simulations. The four U.S. Geological Survey simulation scenarios were the No Action scenario (EIS Alternative 1), the Downstream Diversion scenario (EIS Alternative 2), the Upstream Return-Flow scenario (EIS Alternative 4), and the Upstream Diversion scenario (EIS Alternative 7). Additionally, the results of an Existing Conditions scenario (water years 2000 through

  16. Some mismatches occurred when simulating fractured reservoirs as homogeneous porous media

    SciTech Connect

    Mario Cesar Suarez Arriaga; Fernando Samaniego V.; Fernando Rodriguez

    1996-01-24

    The understanding of transport processes that occur in naturally fractured geothermal systems is far from being complete. Often, evaluation and numerical simulations of fractured geothermal reservoirs, are carried out by assuming equivalent porous media and homogeneous petrophysical properties within big matrix blocks. The purpose of this paper, is to present a comparison between results obtained from numerical studies of a naturally fractured reservoir treated as a simple porous medium and the simulation of some real aspects of the fractured reservoir. A general conclusion outlines the great practical importance of considering even approximately, the true nature of such systems. Our results show that the homogeneous simplified evaluation of the energy resource in a fractured system, could result in unrealistic estimates of the reservoir capacity to generate electricity.

  17. Mathematical simulation of gas-liquid mixture flow in a reservoir and a wellbore with allowance for the dynamical interactions in the reservoir-well system

    NASA Astrophysics Data System (ADS)

    Abbasov, E. M.; Feyzullayev, Kh. A.

    2016-01-01

    Fluid dynamic processes related to mature oil field development are simulated by applying a numerical algorithm based on the gas-liquid mixture flow equations in a reservoir and a wellbore with allowance for the dynamical interaction in the reservoir-well system. Numerical experiments are performed in which well production characteristics are determined from wellhead parameters.

  18. Particle-based simulation of hydraulic fracture and fluid/heat flow in geothermal reservoirs

    NASA Astrophysics Data System (ADS)

    Mora, Peter; Wang, Yucang; Alonso-Marroquin, Fernando

    2013-06-01

    Realizing the potential of geothermal energy as a cheap, green, sustainable resource to provide for the planet's future energy demands that a key geophysical problem be solved first: how to develop and maintain a network of multiple fluid flow pathways for the time required to deplete the heat within a given region. We present the key components for micro-scale particle-based numerical modeling of hydraulic fracture, and fluid and heat flow in geothermal reservoirs. They are based on the latest developments of ESyS-Particle - the coupling of the Lattice Solid Model (LSM) to simulate the nonlinear dynamics of complex solids with the Lattice Boltzmann Method (LBM) applied to the nonlinear dynamics of coupled fluid and heat flow in the complex solid-fluid system. The coupled LSM/LBM can be used to simulate development of fracture systems in discontinuous media, elastic stress release, fluid injection and the consequent slip at joint surfaces, and hydraulic fracturing; heat exchange between hot rocks and water within flow pathways created through hydraulic fracturing; and fluid flow through complex, narrow, compact and gouge-or powder-filled fracture and joint systems. We demonstrate the coupled LSM/LBM to simulate the fundamental processes listed above, which are all components for the generation and sustainability of the hot-fractured rock geothermal energy fracture systems required to exploit this new green-energy resource.

  19. Advanced Techniques for Reservoir Simulation and Modeling of Non-Conventional Wells

    SciTech Connect

    Durlofsky, Louis J.; Aziz, Khalid

    2001-08-23

    Research results for the second year of this project on the development of improved modeling techniques for non-conventional (e.g., horizontal, deviated or multilateral) wells were presented. The overall program entails the development of enhanced well modeling and general simulation capabilities. A general formulation for black-oil and compositional reservoir simulation was presented.

  20. 3-D-numerical approach to simulate an avalanche impact into a reservoir

    NASA Astrophysics Data System (ADS)

    Gabl, R.; Seibl, J.; Gems, B.; Aufleger, M.

    2015-06-01

    The impact of an avalanche into a reservoir induces an impulse wave, which poses a threat to population and infrastructure. For a good approximation of the generated wave height and length as well as the resulting outflow volume over structures and dams, formulas, which base on different simplifying assumptions, can be used. Further project-specific investigations by means of a scale model test or numerical simulations are advisable for complex reservoirs as well as the inclusion of hydraulic structures such as spillways. The paper presents a new approach for a 3-D-numerical simulation of an avalanche impact into a reservoir. In this model concept the energy and mass of the avalanche are represented by accelerated water on the real hill slope. Instead of snow, only water and air are used to simulate the moving avalanche with the software FLOW-3D. A significant advantage of this assumption is the self-adaptation of the model avalanche onto the terrain. In order to reach good comparability of the results with existing research at the ETH Zürich, a simplified reservoir geometry is investigated. Thus, a reference case has been analysed including a variation of three geometry parameters (still water depth in the reservoir, freeboard of the dam and reservoir width).

  1. Application of advanced reservoir characterization, simulation and production optimization strategies to maximize recovery in slope and basin clastic reservoirs, West Texas (Delaware Basin). Annual report

    SciTech Connect

    Dutton, S.P.; Asquith, G.B.; Barton, M.D.; Cole, A.G.; Gogas, J.; Malik, M.A.; Clift, S.J.; Guzman, J.I.

    1997-11-01

    The objective of this project is to demonstrate that detailed reservoir characterization of slope and basin clastic reservoirs in sandstones of the Delaware Mountain Group in the Delaware Basin of West Texas and New Mexico is a cost-effective way to recover a higher percentage of the original oil in place through strategic placement of infill wells and geologically based field development. This project involves reservoir characterization of two Late Permian slope and basin clastic reservoirs in the Delaware Basin, West Texas, followed by a field demonstration in one of the fields. The fields being investigated are Geraldine Ford and Ford West fields in Reeves and Culberson Counties, Texas. Project objectives are divided into two major phases, reservoir characterization and implementation. The objectives of the reservoir characterization phase of the project were to provide a detailed understanding of the architecture and heterogeneity of the two fields, the Ford Geraldine unit and Ford West field. Reservoir characterization utilized 3-D seismic data, high-resolution sequence stratigraphy, subsurface field studies, outcrop characterization, and other techniques. Once reservoir characterized was completed, a pilot area of approximately 1 mi{sup 2} at the northern end of the Ford Geraldine unit was chosen for reservoir simulation. This report summarizes the results of the second year of reservoir characterization.

  2. Simulation of irreversible rock compaction effects on geopressured reservoir response: Topical report

    SciTech Connect

    Riney, T.D.

    1986-12-01

    A series of calculations are presented which quantitatively demonstrate the effects of nonlinear stress-deformation properties on the behavior of geopressured reservoirs. The range of stress-deformation parameters considered is based on information available from laboratory rock mechanics tests performed at the University of Texas at Austin and at Terra Tek, Inc. on cores recovered from geopressured wells. The effects of irreversible formation rock compaction, associated permeability reduction, and repetitive load/unload cycling are considered. The formation rock and geopressured brine properties are incorporated into an existing reservoir simulator using a bilinear model for the irreversible compaction process. Pressure drawdown and buildup testing of a well producing from the geopressured formation is simulated for a suite of calculations covering the range of formation parameters. The results are presented and discussed in terms of the inference (e.g., permeability and reservoir volume) that would be drawn from the simulated test data by an analyst using conventional methods.

  3. Problems related to water quality and algal control in Lopez Reservoir, San Luis Obispo County, California

    USGS Publications Warehouse

    Fuller, Richard H.; Averett, Robert C.; Hines, Walter G.

    1975-01-01

    A study to determine the present enrichment status of Liopez Reservoir in San Luis Obispo county, California, and to evaluate copper sulfate algal treatment found that stratification in the reservoir regulates nutrient release and that algal control has been ineffective. Nuisance algal blooms, particularly from March to June, have been a problem in the warm multipurpose reservoir since it was initially filled following intense storms in 1968-69. The cyanophyte Anabaena unispora has been dominant; cospecies are the diatoms Stephanodiscus astraea and Cyclotella operculata, and the chlorophytes Pediastrum deplex and Sphaerocystis schroeteri. During an A. unispora bloom in May 1972 the total lake surface cell count was nearly 100,000 cells/ml. Thermal stratification from late spring through autumn results in oxygen deficiency in the hypolimnion and metalimnion caused by bacterial oxidation of organic detritus. The anaerobic conditions favor chemical reduction of organic matter, which constitute 10-14% of the sediment. As algae die, sink to the bottom, and decompose, nutrients are released to the hypolimnion , and with the autumn overturn are spread to the epilimnion. Algal blooms not only hamper recreation, but through depletion of dissolved oxygen in the epilimnion may have caused periodic fishkills. Copper sulfate mixed with sodium citrate and applied at 1.10-1.73 lbs/acre has not significantly reduced algal growth; a method for determining correct dosage is presented. (Lynch-Wisconsin)

  4. Radioactive Sediment Transport on Ogaki Dam Reservoir in Fukushima Evacuated Zone: Numerical Simulation Studies by 2-D River Simulation Code

    NASA Astrophysics Data System (ADS)

    Yamada, Susumu; Kitamura, Akihiro; Kurikami, Hiroshi; Machida, Masahiko

    2015-04-01

    Fukushima Daiichi Nuclear Power Plant (FDNPP) accident on March 2011 released significant quantities of radionuclides to atmosphere. The most significant nuclide is radioactive cesium isotopes. Therefore, the movement of the cesium is one of the critical issues for the environmental assessment. Since the cesium is strongly sorbed by soil particles, the cesium transport can be regarded as the sediment transport which is mainly brought about by the aquatic system such as a river and a lake. In this research, our target is the sediment transport on Ogaki dam reservoir which is located in about 16 km northwest from FDNPP. The reservoir is one of the principal irrigation dam reservoirs in Fukushima Prefecture and its upstream river basin was heavily contaminated by radioactivity. We simulate the sediment transport on the reservoir using 2-D river simulation code named Nays2D originally developed by Shimizu et al. (The latest version of Nays2D is available as a code included in iRIC (http://i-ric.org/en/), which is a river flow and riverbed variation analysis software package). In general, a 2-D simulation code requires a huge amount of calculation time. Therefore, we parallelize the code and execute it on a parallel computer. We examine the relationship between the behavior of the sediment transport and the height of the reservoir exit. The simulation result shows that almost all the sand that enter into the reservoir deposit close to the entrance of the reservoir for any height of the exit. The amounts of silt depositing within the reservoir slightly increase by raising the height of the exit. However, that of the clay dramatically increases. Especially, more than half of the clay deposits, if the exit is sufficiently high. These results demonstrate that the water level of the reservoir has a strong influence on the amount of the clay discharged from the reservoir. As a result, we conclude that the tuning of the water level has a possibility for controlling the

  5. Simulation of Hydrodynamics at Stratified Reservoirs Using a Staged Modeling Approach

    SciTech Connect

    Khangaonkar, Tarang P.; Yang, Zhaoqing; Paik, Joongcheol; Sotiropoulos, Fotis

    2008-10-01

    Hydropower reservoirs impounded by high-head dams exhibit complex circulation that confuses the downstream migrating salmon and limits successful collection and passage of fish. Fish passage engineers attempt to modify the hydrothermal behavior at reservoirs through structural and operational modifications and often use hydrodynamic simulations to guide their actions. Simulation of key hydrothermal processes such as (a) development of a stable two-layer stratified system, (b) density-driven currents over a reservoir length scale, and (c) discharge hydraulics near the power generation and fish collection intakes requires highly specialized models applied at differing temporal and spatial scales. A staged modeling approach is presented that uses external coupling of models at varying temporal scales and spatial resolution to simulate the entire hydraulic regime from the mouth of the reservoir at the upstream end to the discharge at the dam. The staged modeling approach is illustrated using a case study where structural modifications were evaluated to improve reservoir stratification and density-driven currents. The model results provided input and valuable insight in the development of a new structure design and configuration for effective fish collection near the forebay of a high-head dam.

  6. Application of Advanced Reservoir Characterization, Simulation, and Production Optimization Strategies to Maximize Recovery in Slope, and Basin Clastic Reservoirs, West Texas (Delaware Basin)

    SciTech Connect

    Shirley P. Dutton

    1997-04-30

    The objective of this project is to demonstrate that detailed reservoir characterization of slope and basin clastic reservoirs in sandstones of the Delaware Mountain Group in the Delaware Basin of West Texas and New Mexico is a cost effective way to recover a higher percentage of the original oil in place through strategic placement of infill wells and geologically based field development. Project objectives are divided into two major phases. The objectives of the reservoir characterization phase of the project are to provide a detailed understanding of the architecture and heterogeneity of two fields, the Ford Geraldine unit and Ford West field, which produce from the Bell Canyon and Cherry Canyon Formations, respectively, of the Delaware Mountain Group and to compare Bell Canyon and Cherry Canyon reservoirs. Reservoir characterization will utilize 3-D seismic data, high-resolution sequence stratigraphy, subsurface field studies, outcrop characterization, and other techniques. Once the reservoir-characterization study of both fields is completed, a pilot area of approximately 1 mi 2 in one of the fields will be chosen for reservoir simulation. The objectives of the implementation phase of the project are to (1) apply the knowledge gained from reservoir characterization and simulation studies to increase recovery from the pilot area, (2) demonstrate that economically significant unrecovered oil remains in geologically resolvable untapped compartments, and (3) test the accuracy of reservoir characterization and flow simulation as predictive tools in resource preservation of mature fields. A geologically designed, enhanced-recovery program (CO 2 flood, waterflood, or polymer flood) and well-completion program will be developed, and one to three infill wells will be drilled and cored. Through technology transfer workshops and other presentations, the knowledge gained in the comparative study of these two fields can then be applied to increase production from the more

  7. Application of Advanced Reservoir Characterization, Simulation, and Production Optimization Strategies to Maximize Recovery in Slope and Basin Clastic Reservoirs, West Texas (Delaware Basin).

    SciTech Connect

    Dutton, S.P.

    1997-10-30

    The objective of this project is to demonstrate that detailed reservoir characterization of slope and basin clastic reservoirs in sandstones of the Delaware Mountain Group in the Delaware Basin of West Texas and New Mexico is a cost effective way to recover a higher percentage of the original oil in place through strategic placement of infill wells and geologically based field development. Project objectives are divided into two major phases. The objectives of the reservoir characterization phase of the project are to provide a detailed understanding of the architecture and heterogeneity of two fields, the Ford Geraldine unit and Ford West field, which produce from the Bell Canyon and Cherry Canyon Formations, respectively, of the Delaware Mountain Group and to compare Bell Canyon and Cherry Canyon reservoirs. Reservoir characterization will utilize 3-D seismic data, high-resolution sequence stratigraphy, subsurface field studies, outcrop characterization, and other techniques. Once the reservoir- characterization study of both fields is completed, a pilot area of approximately 1 mi{sup 2} in one of the fields will be chosen for reservoir simulation. The objectives of the implementation phase of the project are to (1) apply the knowledge gained from reservoir characterization and simulation studies to increase recovery from the pilot area, (2) demonstrate that economically significant unrecovered oil remains in geologically resolvable untapped compartments, and (3) test the accuracy of reservoir characterization and flow simulation as predictive tools in resource preservation of mature fields. A geologically designed, enhanced-recovery program (CO{sub 2} flood, water flood, or polymer flood) and well-completion program will be developed, and one to three infill wells will be drilled and cored. Through technology transfer workshops and other present at ions, the knowledge gained in the comparative study of these two fields can then be applied to increase product ion

  8. Problem reporting management system performance simulation

    NASA Technical Reports Server (NTRS)

    Vannatta, David S.

    1993-01-01

    This paper proposes the Problem Reporting Management System (PRMS) model as an effective discrete simulation tool that determines the risks involved during the development phase of a Trouble Tracking Reporting Data Base replacement system. The model considers the type of equipment and networks which will be used in the replacement system as well as varying user loads, size of the database, and expected operational availability. The paper discusses the dynamics, stability, and application of the PRMS and addresses suggested concepts to enhance the service performance and enrich them.

  9. Motor operated valves problems tests and simulations

    SciTech Connect

    Pinier, D.; Haas, J.L.

    1996-12-01

    An analysis of the two refusals of operation of the EAS recirculation shutoff valves enabled two distinct problems to be identified on the motorized valves: the calculation methods for the operating torques of valves in use in the power plants are not conservative enough, which results in the misadjustement of the torque limiters installed on their motorizations, the second problem concerns the pressure locking phenomenon: a number of valves may entrap a pressure exceeding the in-line pressure between the disks, which may cause a jamming of the valve. EDF has made the following approach to settle the first problem: determination of the friction coefficients and the efficiency of the valve and its actuator through general and specific tests and models, definition of a new calculation method. In order to solve the second problem, EDF has made the following operations: identification of the valves whose technology enables the pressure to be entrapped: the tests and numerical simulations carried out in the Research and Development Division confirm the possibility of a {open_quotes}boiler{close_quotes} effect: determination of the necessary modifications: development and testing of anti-boiler effect systems.

  10. Evaluation of Gas Production Potential of Hydrate Deposits in Alaska North Slope using Reservoir Simulations

    NASA Astrophysics Data System (ADS)

    Nandanwar, M.; Anderson, B. J.

    2015-12-01

    Over the past few decades, the recognition of the importance of gas hydrates as a potential energy resource has led to more and more exploration of gas hydrate as unconventional source of energy. In 2002, U.S. Geological Survey (USGS) started an assessment to conduct a geology-based analysis of the occurrences of gas hydrates within northern Alaska. As a result of this assessment, many potential gas hydrate prospects were identified in the eastern National Petroleum Reserve Alaska (NPRA) region of Alaska North Slope (ANS) with total gas in-place of about 2 trillion cubic feet. In absence of any field test, reservoir simulation is a powerful tool to predict the behavior of the hydrate reservoir and the amount of gas that can be technically recovered using best suitable gas recovery technique. This work focuses on the advanced evaluation of the gas production potential of hydrate accumulation in Sunlight Peak - one of the promising hydrate fields in eastern NPRA region using reservoir simulations approach, as a part of the USGS gas hydrate development Life Cycle Assessment program. The main objective of this work is to develop a field scale reservoir model that fully describes the production design and the response of hydrate field. Due to the insufficient data available for this field, the distribution of the reservoir properties (such as porosity, permeability and hydrate saturation) are approximated by correlating the data from Mount Elbert hydrate field to obtain a fully heterogeneous 3D reservoir model. CMG STARS is used as a simulation tool to model multiphase, multicomponent fluid flow and heat transfer in which an equilibrium model of hydrate dissociation was used. Production of the gas from the reservoir is carried out for a period of 30 years using depressurization gas recovery technique. The results in terms of gas and water rate profiles are obtained and the response of the reservoir to pressure and temperature changes due to depressurization and hydrate

  11. Simulation studies to evaluate the effect of fracture closure on the performance of fractured reservoirs; Final report

    SciTech Connect

    Howrie, I.; Dauben, D.

    1994-03-01

    A three-year research program to evaluate the effect of fracture closure on the recovery of oil and gas from naturally fractured reservoirs has been completed. The overall objectives of the study were to: (1) evaluate the reservoir conditions for which fracture closure is significant, and (2) evaluate innovative fluid injection techniques capable of maintaining pressure within the reservoir. The evaluations of reservoir performance were made by a modern dual porosity simulator, TETRAD. This simulator treats both porosity and permeability as functions of pore pressure. The Austin Chalk in the Pearsall Field in of South Texas was selected as the prototype fractured reservoir for this work. During the first year, simulations of vertical and horizontal well performance were made assuming that fracture permeability was insensitive to pressure change. Sensitivity runs indicated that the simulator was predicting the effects of critical reservoir parameters in a logical and consistent manner. The results confirmed that horizontal wells could increase both rate of oil recovery and total oil recovery from naturally fractured reservoirs. In the second year, the performance of the same vertical and horizontal wells was reevaluated with fracture permeability treated as a function of reservoir pressure. To investigate sensitivity to in situ stress, differing loading conditions were assumed. Simulated natural depletions confirm that pressure sensitive fractures degrade well performance. The severity of degradation worsens when the initial reservoir pressure approaches the average stress condition of the reservoir, such as occurs in over pressured reservoirs. Simulations with water injection indicate that degradation of permeability can be counteracted when reservoir pressure is maintained and oil recovery can be increased when reservoir properties are favorable.

  12. Application of a new scale up methodology to the simulation of displacement processes in heterogeneous reservoirs

    SciTech Connect

    Durlofsky, L.J.; Milliken, W.J.; Dehghani, K.; Jones, R.C.

    1994-12-31

    A general method for the scale up of highly detailed, heterogeneous reservoir cross sections is presented and applied to the simulation of several recovery processes in a variety of geologic settings. The scale up technique proceeds by first identifying portions of the fine scale reservoir description which could potentially lead to high fluid velocities, typically regions of connected, high permeability. These regions are then modeled in detail while the remainder of the domain is coarsened using a general numerical technique for the calculation of effective permeability. The overall scale up method is applied to the cross sectional simulation of three actual fields. Waterflood, steamflood and miscible flood recovery processes are considered. In all these cases, the scale up technique is shown to give coarsened reservoir descriptions which provide simulation results in very good agreement with those of the detailed reservoir descriptions. For these simulations, speedups in computation times, for the coarsened models relative to their fine grid counterparts, range from a factor of 10 to a factor of 200.

  13. A combination of streamtube and geostatical simulation methodologies for the study of large oil reservoirs

    SciTech Connect

    Chakravarty, A.; Emanuel, A.S.; Bernath, J.A.

    1997-08-01

    The application of streamtube models for reservoir simulation has an extensive history in the oil industry. Although these models are strictly applicable only to fields under voidage balance, they have proved to be useful in a large number of fields provided that there is no solution gas evolution and production. These models combine the benefit of very fast computational time with the practical ability to model a large reservoir over the course of its history. These models do not, however, directly incorporate the detailed geological information that recent experience has taught is important. This paper presents a technique for mapping the saturation information contained in a history matched streamtube model onto a detailed geostatistically derived finite difference grid. With this technique, the saturation information in a streamtube model, data that is actually statistical in nature, can be identified with actual physical locations in a field and a picture of the remaining oil saturation can be determined. Alternatively, the streamtube model can be used to simulate the early development history of a field and the saturation data then used to initialize detailed late time finite difference models. The proposed method is presented through an example application to the Ninian reservoir. This reservoir, located in the North Sea (UK), is a heterogeneous sandstone characterized by a line drive waterflood, with about 160 wells, and a 16 year history. The reservoir was satisfactorily history matched and mapped for remaining oil saturation. A comparison to 3-D seismic survey and recently drilled wells have provided preliminary verification.

  14. Quantum simulation of dissipative processes without reservoir engineering

    SciTech Connect

    Di Candia, R.; Pedernales, J. S.; del Campo, A.; Solano, E.; Casanova, J.

    2015-05-29

    We present a quantum algorithm to simulate general finite dimensional Lindblad master equations without the requirement of engineering the system-environment interactions. The proposed method is able to simulate both Markovian and non-Markovian quantum dynamics. It consists in the quantum computation of the dissipative corrections to the unitary evolution of the system of interest, via the reconstruction of the response functions associated with the Lindblad operators. Our approach is equally applicable to dynamics generated by effectively non-Hermitian Hamiltonians. We confirm the quality of our method providing specific error bounds that quantify its accuracy.

  15. Quantum Simulation of Dissipative Processes without Reservoir Engineering

    PubMed Central

    Di Candia, R.; Pedernales, J. S.; del Campo, A.; Solano, E.; Casanova, J.

    2015-01-01

    We present a quantum algorithm to simulate general finite dimensional Lindblad master equations without the requirement of engineering the system-environment interactions. The proposed method is able to simulate both Markovian and non-Markovian quantum dynamics. It consists in the quantum computation of the dissipative corrections to the unitary evolution of the system of interest, via the reconstruction of the response functions associated with the Lindblad operators. Our approach is equally applicable to dynamics generated by effectively non-Hermitian Hamiltonians. We confirm the quality of our method providing specific error bounds that quantify its accuracy. PMID:26024437

  16. Solving search problems by strongly simulating quantum circuits

    PubMed Central

    Johnson, T. H.; Biamonte, J. D.; Clark, S. R.; Jaksch, D.

    2013-01-01

    Simulating quantum circuits using classical computers lets us analyse the inner workings of quantum algorithms. The most complete type of simulation, strong simulation, is believed to be generally inefficient. Nevertheless, several efficient strong simulation techniques are known for restricted families of quantum circuits and we develop an additional technique in this article. Further, we show that strong simulation algorithms perform another fundamental task: solving search problems. Efficient strong simulation techniques allow solutions to a class of search problems to be counted and found efficiently. This enhances the utility of strong simulation methods, known or yet to be discovered, and extends the class of search problems known to be efficiently simulable. Relating strong simulation to search problems also bounds the computational power of efficiently strongly simulable circuits; if they could solve all problems in P this would imply that all problems in NP and #P could be solved in polynomial time. PMID:23390585

  17. Deposition and simulation of sediment transport in the Lower Susquehanna River reservoir system

    USGS Publications Warehouse

    Hainly, R.A.; Reed, L.A.; Flippo, H.N., Jr.; Barton, G.J.

    1995-01-01

    The Susquehanna River drains 27,510 square miles in New York, Pennsylvania, and Maryland and is the largest tributary to the Chesapeake Bay. Three large hydroelectric dams are located on the river, Safe Harbor (Lake Clarke) and Holtwood (Lake Aldred) in southern Pennsylvania, and Conowingo (Conowingo Reservoir) in northern Maryland. About 259 million tons of sediment have been deposited in the three reservoirs. Lake Clarke contains about 90.7 million tons of sediment, Lake Aldred contains about 13.6 million tons, and Conowingo Reservoir contains about 155 million tons. An estimated 64.8 million tons of sand, 19.7 million tons of coal, 112 million tons of silt, and 63.3 million tons of clay are deposited in the three reservoirs. Deposition in the reservoirs is variable and ranges from 0 to 30 feet. Chemical analyses of sediment core samples indicate that the three reservoirs combined contain about 814,000 tons of organic nitrogen, 98,900 tons of ammonia as nitrogen, 226,000 tons of phosphorus, 5,610,000 1tons of iron, 2,250,000 tons of aluminum, and about 409,000 tons of manganese. Historical data indicate that Lake Clarke and Lake Aldred have reached equilibrium, and that they no longer store sediment. A comparison of cross-sectional data from Lake Clarke and Lake Aldred with data from Conowingo Reservoir indicates that Conowingo Reservoir will reach equilibrium within the next 20 to 30 years. As the Conowingo Reservoir fills with sediment and approaches equilibrium, the amount of sediment transported to the Chesapeake Bay will increase. The most notable increases will take place when very high flows scour the deposited sediment. Sediment transport through the reservoir system was simulated with the U.S. Army Corps of Engineers' HEC-6 computer model. The model was calibrated with monthly sediment loads for calendar year 1987. Calibration runs with options set for maximum trap efficiency and a "natural" particle-size distribution resulted in an overall computed trap

  18. Optimization or Simulation? Comparison of approaches to reservoir operation on the Senegal River

    NASA Astrophysics Data System (ADS)

    Raso, Luciano; Bader, Jean-Claude; Pouget, Jean-Christophe; Malaterre, Pierre-Olivier

    2015-04-01

    Design of reservoir operation rules follows, traditionally, two approaches: optimization and simulation. In simulation, the analyst hypothesizes operation rules, and selects them by what-if analysis based on effects of model simulations on different objectives indicators. In optimization, the analyst selects operational objective indicators, finding operation rules as an output. Optimization rules guarantee optimality, but they often require further model simplification, and can be hard to communicate. Selecting the most proper approach depends on the system under analysis, and the analyst expertise and objectives. We present advantage and disadvantages of both approaches, and we test them for the Manantali reservoir operation rule design, on the Senegal River, West Africa. We compare their performance in attaining the system objectives. Objective indicators are defined a-priori, in order to quantify the system performance. Results from this application are not universally generalizable to the entire class, but they allow us to draw conclusions on this system, and to give further information on their application.

  19. Mathematical Simulation of Contaminant Flow in Closed Reservoir

    NASA Astrophysics Data System (ADS)

    Agranat, V. M.; Goudov, A. M.; Perminov, V. A.

    2016-01-01

    A mathematical model of the propagation in flooded mine lightweight contaminant due to allocation of groundwater is considered. Mathematical model was based on an analysis of experimental data and using concept and methods from reactive media mechanics. The boundary-value problem is solved numerically using the finite volume method. The distribution of fields of velocities and concentration of impurity particles in a flooded mine have been obtained at different times. These results can be used to analyze mining water treatment process due to environment and evaluate its further possible improvements.

  20. OGS#PETSc approach for robust and efficient simulations of strongly coupled hydrothermal processes in EGS reservoirs

    NASA Astrophysics Data System (ADS)

    Watanabe, Norihiro; Blucher, Guido; Cacace, Mauro; Kolditz, Olaf

    2016-04-01

    A robust and computationally efficient solution is important for 3D modelling of EGS reservoirs. This is particularly the case when the reservoir model includes hydraulic conduits such as induced or natural fractures, fault zones, and wellbore open-hole sections. The existence of such hydraulic conduits results in heterogeneous flow fields and in a strengthened coupling between fluid flow and heat transport processes via temperature dependent fluid properties (e.g. density and viscosity). A commonly employed partitioned solution (or operator-splitting solution) may not robustly work for such strongly coupled problems its applicability being limited by small time step sizes (e.g. 5-10 days) whereas the processes have to be simulated for 10-100 years. To overcome this limitation, an alternative approach is desired which can guarantee a robust solution of the coupled problem with minor constraints on time step sizes. In this work, we present a Newton-Raphson based monolithic coupling approach implemented in the OpenGeoSys simulator (OGS) combined with the Portable, Extensible Toolkit for Scientific Computation (PETSc) library. The PETSc library is used for both linear and nonlinear solvers as well as MPI-based parallel computations. The suggested method has been tested by application to the 3D reservoir site of Groß Schönebeck, in northern Germany. Results show that the exact Newton-Raphson approach can also be limited to small time step sizes (e.g. one day) due to slight oscillations in the temperature field. The usage of a line search technique and modification of the Jacobian matrix were necessary to achieve robust convergence of the nonlinear solution. For the studied example, the proposed monolithic approach worked even with a very large time step size of 3.5 years.

  1. Application of Advanced Reservoir Characterization, Simulation, and Production Optimization Strategies to Maximize Recovery in Slope and Basin Clastic Reservoirs, West Texas (Delaware Basin)

    SciTech Connect

    Shirley P. Dutton

    1997-07-30

    The objective of this project is to demonstrate that detailed reservoir characterization of slope and basin clastic reservoirs in sandstones of the Delaware Mountain Group in the Delaware Basin of West Texas and New Mexico is a cost effective way to recover a higher percentage of the original oil in place through strategic placement of infill wells and geologically based field development. Project objectives are divided into two major phases. The objectives of the reservoir characterization phase of the project are to provide a detailed understanding of the architecture and heterogeneity of two fields, the Ford Geraldine unit and Ford West field, which produce from the Bell Canyon and Cherry Canyon Formations, respectively, of the Delaware Mountain Group and to compare Bell Canyon and Cherry Canyon reservoirs. Reservoir characterization will utilize 3-D seismic data, high-resolution sequence stratigraphy, subsurface field studies, outcrop characterization, and other techniques. Once the reservoir-characterization study of both fields is completed, a pilot area of approximately 1 mi 2 in one of the fields will be chosen for reservoir simulation.

  2. Reservoir Simulation on the Cerro Prieto Geothermal Field: A Continuing Study

    SciTech Connect

    Castaneda, M.; Marquez, R.; Arellano, V.; Esquer, C.A.

    1983-12-15

    The Cerro Prieto geothermal field is a liquid-dominated geothermal reservoir of complex geological and hydrological structure. It is located at the southern end of the Salton-Mexicali trough which includes other geothermal anomalies as Heber and East Mesa. Although in 1973, the initial power plant installed capacity was 75 MW of electrical power, this amount increased to 180 MW in 1981 as field development continued. It is expected to have a generating capacity of 620 MW by the end of 1985, when two new plants will be completely in operation. Questions about field deliverability, reservoir life and ultimate recovery related to planned installations are being presently asked. Numerical modeling studies can give very valuable answers to these questions, even at the early stages in the development of a field. An effort to simulate the Cerro Prieto geothermal reservoir has been undergoing for almost two years. A joint project among Comision Federal de Electricidad (CFE), Instituto de Investigaciones Electricas (IIE) and Intercomp of Houstin, Texas, was created to perform reservoir engineering and simulation studies on this field. The final project objective is tosimulate the behavior of the old field region when production from additional wells located in the undeveloped field zones will be used for feeding the new power plants.

  3. Geochemical simulations on CO2-fluid-rock interactions in EGS reservoirs

    NASA Astrophysics Data System (ADS)

    Pan, F.; McPherson, B. J.; Lichtner, P. C.; Kaszuba, J. P.; Lo Re, C.; Karra, S.; Lu, C.; Xu, T.

    2012-12-01

    Supercritical CO2 has been suggested as a heat transmission fluid in Enhanced Geothermal Systems (EGS) reservoirs to improve energy extraction. Understanding the geochemical processes of CO2-fluid-rock interactions in EGS reservoirs is significant important to investigate the performance of energy extraction with CO2 instead of water as a working fluid, carbon sequestration and risk assessment. The objectives of this study: (1) to calibrate and evaluate the kinetic rate constants and specific reactive surface areas of minerals based on the batch experimental data conducted by other researchers (collaborators Kaszuba and Lo Ré at the University of Wyoming); (2) to investigate the effects of CO2-fluid-rock geochemical interactions on the energy extraction efficiency, carbon sequestration, and risk assessment. A series of laboratory experiments were conducted (Lo Ré et al., 2012) to investigate the geochemical reactions among water, fractured granite rocks, and injected supercritical CO2 at elevated temperatures of 250 oC, and pressures of 250-450 bars. The batch simulations were firstly conducted to mimic the laboratory experiments with the calibration of mineral reactive surface areas using TOUGHREACT model and parameter estimation software (PEST). Then, we performed 2-D geochemical modeling to simulate the chemical interactions among CO2, fluids, and rocks at high temperatures and pressures of EGS reservoirs. We further investigated the effects of fluid-rock interactions on the energy extraction, carbon sequestration, and risk assessment with CO2 as a heat transmission fluid instead of water for EGS reservoirs. Results of carbonate mineral precipitations suggested that the CO2 as a working fluid instead of water was favorable for EGS reservoirs on the CO2 sequestration. Our simulations also suggested that the energy extraction could be enhanced using CO2 as the transmission fluid compared to water.

  4. Simulating California reservoir operation using the classification and regression-tree algorithm combined with a shuffled cross-validation scheme

    NASA Astrophysics Data System (ADS)

    Yang, Tiantian; Gao, Xiaogang; Sorooshian, Soroosh; Li, Xin

    2016-03-01

    The controlled outflows from a reservoir or dam are highly dependent on the decisions made by the reservoir operators, instead of a natural hydrological process. Difference exists between the natural upstream inflows to reservoirs and the controlled outflows from reservoirs that supply the downstream users. With the decision maker's awareness of changing climate, reservoir management requires adaptable means to incorporate more information into decision making, such as water delivery requirement, environmental constraints, dry/wet conditions, etc. In this paper, a robust reservoir outflow simulation model is presented, which incorporates one of the well-developed data-mining models (Classification and Regression Tree) to predict the complicated human-controlled reservoir outflows and extract the reservoir operation patterns. A shuffled cross-validation approach is further implemented to improve CART's predictive performance. An application study of nine major reservoirs in California is carried out. Results produced by the enhanced CART, original CART, and random forest are compared with observation. The statistical measurements show that the enhanced CART and random forest overperform the CART control run in general, and the enhanced CART algorithm gives a better predictive performance over random forest in simulating the peak flows. The results also show that the proposed model is able to consistently and reasonably predict the expert release decisions. Experiments indicate that the release operation in the Oroville Lake is significantly dominated by SWP allocation amount and reservoirs with low elevation are more sensitive to inflow amount than others.

  5. A methodology for incorporating geomechanically-based fault damage zones models into reservoir simulation

    NASA Astrophysics Data System (ADS)

    Paul, Pijush Kanti

    In the fault damage zone modeling study for a field in the Timor Sea, I present a methodology to incorporate geomechanically-based fault damage zones into reservoir simulation. In the studied field, production history suggests that the mismatch between actual production and model prediction is due to preferential fluid flow through the damage zones associated with the reservoir scale faults, which is not included in the baseline petrophysical model. I analyzed well data to estimate stress heterogeneity and fracture distributions in the reservoir. Image logs show that stress orientations are homogenous at the field scale with a strike-slip/normal faulting stress regime and maximum horizontal stress oriented in NE-SW direction. Observed fracture zones in wells are mostly associated with well scale fault and bed boundaries. These zones do not show any anomalies in production logs or well test data, because most of the fractures are not optimally oriented to the present day stress state, and matrix permeability is high enough to mask any small anomalies from the fracture zones. However, I found that fracture density increases towards the reservoir scale faults, indicating high fracture density zones or damage zones close to these faults, which is consistent with the preferred flow direction indicated by interference and tracer test done between the wells. It is well known from geologic studies that there is a concentration of secondary fractures and faults in a damage zone adjacent to larger faults. Because there is usually inadequate data to incorporate damage zone fractures and faults into reservoir simulation models, in this study I utilized the principles of dynamic rupture propagation from earthquake seismology to predict the nature of fractured/damage zones associated with reservoir scale faults. The implemented workflow can be used to more routinely incorporate damage zones into reservoir simulation models. Applying this methodology to a real reservoir utilizing

  6. Numerical Simulation of Injectivity Effects of Mineral Scaling and Clay Swelling in a Fractured Geothermal Reservoir

    SciTech Connect

    Xu, Tianfu; Pruess, Karsten

    2004-05-10

    A major concern in the development of hot dry rock (HDR) and hot fractured rock (HFR) reservoirs is achieving and maintaining adequate injectivity, while avoiding the development of preferential short-circuiting flow paths such as those caused by thermally-induced stress cracking. Past analyses of HDR and HFR reservoirs have tended to focus primarily on the coupling between hydrology (flow), heat transfer, and rock mechanics. Recent studies suggest that rock-fluid interactions and associated mineral dissolution and precipitation effects could have a major impact on the long-term performance of HFR reservoirs. The present paper uses recent European studies as a starting point to explore chemically-induced effects of fluid circulation in HFR systems. We examine ways in which the chemical composition of reinjected waters can be modified to improve reservoir performance by maintaining or even enhancing injectivity. Chemical manipulations considered here include pH modification and dilution with fresh water. We performed coupled thermo-hydrologic-chemical simulations in which the fractured medium was represented by a one-dimensional MINC model (multiple interacting continua), using the non-isothermal multi-phase reactive geochemical transport code TOUGHREACT. Results indicate that modifying the injection water chemistry can enhance mineral dissolution and reduce clay swelling. Chemical interactions between rocks and fluids will change a HFR reservoir over time, with some changes favorable and others not. A detailed, quantitative understanding of processes and mechanisms can suggest chemical methods for reservoir management, which may be employed to improve the performance of the geothermal system.

  7. Analysis of numerical simulations and influencing factors of seasonal manganese pollution in reservoirs.

    PubMed

    Peng, Hui; Zheng, Xilai; Chen, Lei; Wei, Yang

    2016-07-01

    Seasonal manganese pollution has become an increasingly pressing water quality issue for water supply reservoirs in recent years. Manganese is a redox-sensitive element and is released from sediment under anoxic conditions near the sediment-water interface during summer and autumn, when water temperature stratification occurs. The reservoir water temperature and water dynamic conditions directly influence the formation of manganese pollution. Numerical models are useful tools to quantitatively evaluate manganese pollution and its influencing factors. This paper presents a reservoir manganese pollution model by adding a manganese biogeochemical module to a water quality model-CE-QUAL-W2. The model is applied to the Wangjuan reservoir (Qingdao, China), which experiences manganese pollution during summer and autumn. Field data are used to verify the model, and the results show that the model can reproduce the main features of the thermal stratification and manganese distribution. The model is used to evaluate the manganese pollution process and its four influencing factors, including air temperature, water level, wind speed, and wind directions, through different simulation scenarios. The results show that all four factors can influence manganese pollution. High air temperature, high water level, and low wind speed aggravate manganese pollution, while low air temperature, low water level, and high wind speed reduce manganese pollution. Wind that travels in the opposite direction of the flow aggravates manganese pollution, while wind in the same direction as the flow reduces manganese pollution. This study provides useful information to improve our understanding of seasonal manganese pollution in reservoirs, which is important for reservoir manganese pollution warnings and control. PMID:27068892

  8. Large-scale three-dimensional geothermal reservoir simulation on small computer systems

    SciTech Connect

    Antunez, E.; Moridis, G.; Pruess, K.

    1995-05-01

    The performance of TOUGH2, Lawrence Berkeley Laboratory`s general purpose simulator for mass and heat flow and transport enhanced with the addition of a set of preconditioned conjugate gradient solvers, was tested on three PCs (486-33, 486-66, Pentium-90), a MacIntosh Quadra 800, and a workstation IBM RISC 6000. A two-phase, single porosity, 3-D geothermal reservoir model with 1,411 irregular grid blocks, with production from and injection into the reservoir was used as the test model. The code modifications to TOUGH2 and its setup in each machine environment are described. Computational work per time step and CPU time requirements are reported for each of the machines used. It is concluded that the current PCs provide the best price/performance platform for running large-scale geothermal field simulations that just a few years ago could only be executed on mainframe computers and high-end workstations.

  9. Some Current Problems in Simulator Design, Testing and Use.

    ERIC Educational Resources Information Center

    Caro, Paul W.

    Concerned with the general problem of the effectiveness of simulator training, this report reflects information developed during the conduct of aircraft simulator training research projects sponsored by the Air Force, Army, Navy, and Coast Guard. Problems are identified related to simulator design, testing, and use, all of which impact upon…

  10. Faults simulations for three-dimensional reservoir-geomechanical models with the extended finite element method

    NASA Astrophysics Data System (ADS)

    Prévost, Jean H.; Sukumar, N.

    2016-01-01

    Faults are geological entities with thicknesses several orders of magnitude smaller than the grid blocks typically used to discretize reservoir and/or over-under-burden geological formations. Introducing faults in a complex reservoir and/or geomechanical mesh therefore poses significant meshing difficulties. In this paper, we consider the strong-coupling of solid displacement and fluid pressure in a three-dimensional poro-mechanical (reservoir-geomechanical) model. We introduce faults in the mesh without meshing them explicitly, by using the extended finite element method (X-FEM) in which the nodes whose basis function support intersects the fault are enriched within the framework of partition of unity. For the geomechanics, the fault is treated as an internal displacement discontinuity that allows slipping to occur using a Mohr-Coulomb type criterion. For the reservoir, the fault is either an internal fluid flow conduit that allows fluid flow in the fault as well as to enter/leave the fault or is a barrier to flow (sealing fault). For internal fluid flow conduits, the continuous fluid pressure approximation admits a discontinuity in its normal derivative across the fault, whereas for an impermeable fault, the pressure approximation is discontinuous across the fault. Equal-order displacement and pressure approximations are used. Two- and three-dimensional benchmark computations are presented to verify the accuracy of the approach, and simulations are presented that reveal the influence of the rate of loading on the activation of faults.

  11. Simulation studies to evaluate the effect of fracture closure on the performance of naturally fractured reservoirs. Annual report

    SciTech Connect

    Not Available

    1991-10-01

    The first of a three-year research program to evaluate the effect of fracture closure on the recovery of oil and gas from naturally fractured reservoirs has been completed. The objectives of the study are to (1) evaluate the reservoir conditions where fracture closure is significant, and (2) evaluate innovative fluid injection techniques capable of maintaining pressure within the reservoir. Simulation studies were conducted with a dual porosity simulator capable of simulating the performance of vertical and horizontal wells. Each simulator was initialized using properties typical of the Austin Chalk reservoir in Pearsall Field, Texas. Simulations of both vertical and horizontal well performance were made assuming that fracture permeability was insensitive to pressure change. Sensitivity runs indicate that the simulator is predicting the effects of critical reservoir parameters in a logical and consistent manner. The results to-date confirm that horizontal wells can increase both oil recovery rate and total oil recovery from naturally fractured reservoirs. The year one simulation results will provide the baseline for the ongoing study which will evaluate the performance degradation caused by the sensitivity of fracture permeability to pressure change, and investigate fluid injection pressure maintenance as a means to improve oil recovery performance. The study is likely to conclude that fracture closure decreases oil recovery and that pressure support achieved through fluid injection could be beneficial in improving recovery.

  12. Computer simulation of reservoir depletion and oil flow from the Macondo well following the Deepwater Horizon blowout

    USGS Publications Warehouse

    Hsieh, Paul

    2010-01-01

    This report describes the application of a computer model to simulate reservoir depletion and oil flow from the Macondo well following the Deepwater Horizon blowout. Reservoir and fluid data used for model development are based on (1) information released in BP's investigation report of the incident, (2) information provided by BP personnel during meetings in Houston, Texas, and (3) calibration by history matching to shut-in pressures measured in the capping stack during the Well Integrity Test. The model is able to closely match the measured shut-in pressures. In the simulation of the 86-day period from the blowout to shut in, the simulated reservoir pressure at the well face declines from the initial reservoir pressure of 11,850 pounds per square inch (psi) to 9,400 psi. After shut in, the simulated reservoir pressure recovers to a final value of 10,300 psi. The pressure does not recover back to the initial pressure owing to reservoir depletion caused by 86 days of oil discharge. The simulated oil flow rate declines from 63,600 stock tank barrels per day just after the Deepwater Horizon blowout to 52,600 stock tank barrels per day just prior to shut in. The simulated total volume of oil discharged is 4.92 million stock tank barrels. The overall uncertainty in the simulated flow rates and total volume of oil discharged is estimated to be + or - 10 percent.

  13. Simulation of extreme reservoir level distribution with the SCHADEX method (EXTRAFLO project)

    NASA Astrophysics Data System (ADS)

    Paquet, Emmanuel; Penot, David; Garavaglia, Federico

    2013-04-01

    The standard practice for the design of dam spillways structures and gates is to consider the maximum reservoir level reached for a given hydrologic scenario. This scenario has several components: peak discharge, flood volumes on different durations, discharge gradients etc. Within a probabilistic analysis framework, several scenarios can be associated with different return times, although a reference return level (e.g. 1000 years) is often prescribed by the local regulation rules or usual practice. Using continuous simulation method for extreme flood estimation is a convenient solution to provide a great variety of hydrological scenarios to feed a hydraulic model of dam operation: flood hydrographs are explicitly simulated by a rainfall-runoff model fed by a stochastic rainfall generator. The maximum reservoir level reached will be conditioned by the scale and the dynamics of the generated hydrograph, by the filling of the reservoir prior to the flood, and by the dam gates and spillway operation during the event. The simulation of a great number of floods will allow building a probabilistic distribution of maximum reservoir levels. A design value can be chosen at a definite return level. An alternative approach is proposed here, based on the SCHADEX method for extreme flood estimation, proposed by Paquet et al. (2006, 2013). SCHADEX is a so-called "semi-continuous" stochastic simulation method in that flood events are simulated on an event basis and are superimposed on a continuous simulation of the catchment saturation hazard using rainfall-runoff modelling. The SCHADEX process works at the study time-step (e.g. daily), and the peak flow distribution is deduced from the simulated daily flow distribution by a peak-to-volume ratio. A reference hydrograph relevant for extreme floods is proposed. In the standard version of the method, both the peak-to-volume and the reference hydrograph are constant. An enhancement of this method is presented, with variable peak

  14. Simulated effects of proposed reservoir-development alternatives on streamflow quantity in the White River, Colorado and Utah

    USGS Publications Warehouse

    Kuhn, Gerhard; Ellis, S.R.

    1984-01-01

    Numerous reservoirs have been proposed for the White River basin in Colorado and Utah, primarily to provide water for oil-shale development. A multireservoir-flow model was used to simulate the effects of streamflow withdrawal at four of the proposed reservoirs using historical streamflow data from the 1932-81 water years. The proposed reservoirs considered in the study were Avery, Powell Park, Taylor Draw, and White River Reservoirs; construction of Taylor Draw Dam was completed during the study. Annual streamflow depletions from the White River ranging from about 93,000 to 226,000 acre-feet were simulated for the 50 year period. Simulated streamflow throughout the year generally became smaller and more constant as streamflow throughout the year generally became smaller and more constant as streamflow depletion increased. Minimum streamflow requirements would not have been met for a maximum of 13 years and water-use requirements associated with the proposed reservoirs would not have been met for a maximum of 3 years. The current water-use pattern, which depletes about 40,000 acre-feet per year and is dominated by irrigation of hay meadows and pastureland, was maintained in the simulation. Relations between reservoir active capacity and yield applicable to the White River also were developed. These relations show that reservoir storage of about 400,000 acre-feet is the maximum practicable for the White River. (USGS)

  15. Getting it Together: An Environmental Simulation Problem

    ERIC Educational Resources Information Center

    Hackett, Jay K.

    1975-01-01

    Describes an environmental simulation exercise conducted at an outdoor camp in Colorado. The two simulations involved the physical, social, political and environmental factors which might influence construction of an airport in the Rocky Mountain National Park and factors associated with a nuclear generating plant. (BR)

  16. Measuring the instructional validity of clinical simulation problems.

    PubMed

    Feinstein, E; Gustavson, L P; Levine, H G

    1983-03-01

    Written clinical simulation problems in two formats--forced-choiced and essay--were used to test junior and senior medical students at the conclusion of their pediatric rotations. A comparison was made in the performance of students with varying levels of clinical experience. There seemed to be no consistent pattern of improvement with increased instructional time on the forced-choice management problems. Junior students did show improvement over time on the essay management problems, which also seemed to reflect the increased problem-solving and organizational skills of seniors in comparison to juniors. Correlations across problems and correlations between the clinical problem test and other measurement techniques were very weak, partly due to the low sampling reliability of clinical simulation problems. In this study, clinical simulation problems failed to demonstrate responsiveness to development and maturation in the problem-solving approach to patient care. PMID:10259952

  17. Characterization and simulation of an exhumed fractured petroleum reservoir. Final report, March 18, 1996--September 30, 1998

    SciTech Connect

    Forster, C.B.; Nielson, D.L.; Deo, M.

    1998-12-01

    An exhumed fractured reservoir located near Alligator Ridge in central Nevada provides the basis for developing and testing different approaches for simulating fractured petroleum reservoirs. The fractured analog reservoir comprises a 90 m thickness of silty limestone and shaly interbeds within the Devonian Pilot Shale. A period of regional compression followed by ongoing basin and range extension has created faults and fractures that, in tern, have controlled the migration of both oil and gold ore-forming fluids. Open pit gold mines provide access for observing oil seepage, collecting the detailed fracture data needed to map variations in fracture intensity near faults, build discrete fracture network models and create equivalent permeability structures. Fault trace patterns mapped at the ground surface provide a foundation for creating synthetic fault trace maps using a stochastic procedure conditioned by the outcrop data. Conventional simulations of petroleum production from a 900 by 900 m sub-domain within the reservoir analog illustrate the possible influence of faults and fractures on production. The consequences of incorporating the impact of different stress states (e.g., extension, compression or lithostatic) are also explored. Simulating multiphase fluid flow using a discrete fracture, finite element simulator illustrates how faults acting as conduits might be poorly represented by the upscaling procedures used to assign equivalent permeability values within reservoir models. The parallelized reservoir simulators developed during this project provide a vehicle to evaluate when it might be necessary to incorporate very fine scale grid networks in conventional reservoir simulators or to use finely gridded discrete fracture reservoir simulators.

  18. Feasibility of steam injection process in a thin, low-permeability heavy oil reservoir of Arkansas -- a numerical simulation study

    SciTech Connect

    Sarkar, A.K.; Sarathi, P.S.

    1993-12-01

    This report details the findings of an in-depth study undertaken to assess the viability of the steam injection process in the heavy oil bearing Nacatoch sands of Arkansas. Published screening criteria and DOE`s steamflood predictive models were utilized to screen and select reservoirs for further scrutiny. Although, several prospects satisfied the steam injection screening criteria, only a single candidate was selected for detailed simulation studies. The selection was based on the availability of needed data for simulation and the uniqueness of the reservoir. The reservoir investigated is a shallow, thin, low-permeability reservoir with low initial oil saturation and has an underlying water sand. The study showed that the reservoir will respond favorably to steamdrive, but not to cyclic steaming. Steam stimulation, however, is necessary to improve steam injectivity during subsequent steamdrive. Further, in such marginal heavy oil reservoirs (i.e., reservoir characterized by thin pay zone and low initial oil saturation) conventional steamdrive (i.e., steam injection using vertical wells) is unlikely to be economical, and nonconventional methods must be utilized. It was found that the use of horizontal injectors and horizontal producers significantly improved the recovery and oil-steam ratio and improved the economics. It is recommended that the applicability of horizontal steam injection technology in this reservoir be further investigated.

  19. Numerical simulations of CO2 -assisted gas production from hydrate reservoirs

    NASA Astrophysics Data System (ADS)

    Sridhara, P.; Anderson, B. J.; Myshakin, E. M.

    2015-12-01

    A series of experimental studies over the last decade have reviewed the feasibility of using CO2 or CO2+N2 gas mixtures to recover CH4 gas from hydrates deposits. That technique would serve the dual purpose of CO2 sequestration and production of CH4 while maintaining the geo-mechanical stability of the reservoir. In order to analyze CH4 production process by means of CO2 or CO2+N2 injection into gas hydrate reservoirs, a new simulation tool, Mix3HydrateResSim (Mix3HRS)[1], was previously developed to account for the complex thermodynamics of multi-component hydrate phase and to predict the process of CH4 substitution by CO2 (and N2) in the hydrate lattice. In this work, Mix3HRS is used to simulate the CO2 injection into a Class 2 hydrate accumulation characterized by a mobile aqueous phase underneath a hydrate bearing sediment. That type of hydrate reservoir is broadly confirmed in permafrost and along seashore. The production technique implies a two-stage approach using a two-well design, one for an injector and one for a producer. First, the CO2 is injected into the mobile aqueous phase to convert it into immobile CO2 hydrate and to initiate CH4 release from gas hydrate across the hydrate-water boundary (generally designating the onset of a hydrate stability zone). Second, CH4 hydrate decomposition is induced by the depressurization method at a producer to estimate gas production potential over 30 years. The conversion of the free water phase into the CO2 hydrate significantly reduces competitive water production in the second stage, thereby improving the methane gas production. A base case using only the depressurization stage is conducted to compare with enhanced gas production predicted by the CO2-assisted technique. The approach also offers a possibility to permanently store carbon dioxide in the underground formation to greater extent comparing to a direct injection of CO2 into gas hydrate sediment. Numerical models are based on the hydrate formations at the

  20. The big fat LARS - a LArge Reservoir Simulator for hydrate formation and gas production

    NASA Astrophysics Data System (ADS)

    Beeskow-Strauch, Bettina; Spangenberg, Erik; Schicks, Judith M.; Giese, Ronny; Luzi-Helbing, Manja; Priegnitz, Mike; Klump, Jens; Thaler, Jan; Abendroth, Sven

    2013-04-01

    Simulating natural scenarios on lab scale is a common technique to gain insight into geological processes with moderate effort and expenses. Due to the remote occurrence of gas hydrates, their behavior in sedimentary deposits is largely investigated on experimental set ups in the laboratory. In the framework of the submarine gas hydrate research project (SUGAR) a large reservoir simulator (LARS) with an internal volume of 425 liter has been designed, built and tested. To our knowledge this is presently a word-wide unique set up. Because of its large volume it is suitable for pilot plant scale tests on hydrate behavior in sediments. That includes not only the option of systematic tests on gas hydrate formation in various sedimentary settings but also the possibility to mimic scenarios for the hydrate decomposition and subsequent natural gas extraction. Based on these experimental results various numerical simulations can be realized. Here, we present the design and the experimental set up of LARS. The prerequisites for the simulation of a natural gas hydrate reservoir are porous sediments, methane, water, low temperature and high pressure. The reservoir is supplied by methane-saturated and pre-cooled water. For its preparation an external gas-water mixing stage is available. The methane-loaded water is continuously flushed into LARS as finely dispersed fluid via bottom-and-top-located sparger. The LARS is equipped with a mantle cooling system and can be kept at a chosen set temperature. The temperature distribution is monitored at 14 reasonable locations throughout the reservoir by Pt100 sensors. Pressure needs are realized using syringe pump stands. A tomographic system, consisting of a 375-electrode-configuration is attached to the mantle for the monitoring of hydrate distribution throughout the entire reservoir volume. Two sets of tubular polydimethylsiloxan-membranes are applied to determine gas-water ratio within the reservoir using the effect of permeability

  1. Geology and Petrophysical Characterization of the Ferron Sandstone for 3-D Simulation of a Fluvial-Deltaic Reservoir

    SciTech Connect

    Ann Mattson; Craig B. Forster; Paul B. Anderson; Steve H. Snelgrove; Thomas C. Chidsey, Jr.

    1997-05-20

    The objective of this project is to develop a comprehensive, interdisciplinary, and quantitative characterization of a fluvial-deltaic reservoir which will allow realistic inter-well and reservoir-scale modeling to be constructed for improved oil-field development in similar reservoirs world-wide. The geological and petrophysical properties of the Cretaceous Ferron Sandstone in east-central Utah will be quantitatively determined. Both new and existing data will be integrated into a three-dimensional representation of spatial variations in porosity, storativity, and tensorial rock permeability at a scale appropriate for inter-well to regional-scale reservoir simulation. Four activities continued this quarter as part of the geological and petrophysical characterization of the fluvial-deltaic Ferron Sandstone in the Ivie Creek case-study area: (1) regional stratigraphic interpretation, (2) case-study evaluation, (3) reservoir modeling, and (4) technology transfer.

  2. Prediction of effects of hydraulic fracturing using reservoir and well flow simulation

    SciTech Connect

    Mineyuki Hanano; Tayuki Kondo

    1992-01-01

    This paper presents a method to predict and evaluate effects of hydraulic fracturing jobs by using reservoir and well flow numerical simulation. The concept of the method i5 that steam production rate at the operating well head pressure is predicted with different fracture conditions which would be attained by the hydraulic fracturing jobs. Then, the effects of the hydraulic fracturing is evaluated by comparing the predicted steam production rate and that before the hydraulic fracturing. This course of analysis will suggest how large fracture should be created by the fracturing job to attain large enough increase in steam production at the operating condition and the best scheme of the hydraulic fracturing job.

  3. Parallel, Multigrid Finite Element Simulator for Fractured/Faulted and Other Complex Reservoirs based on Common Component Architecture (CCA)

    SciTech Connect

    Milind Deo; Chung-Kan Huang; Huabing Wang

    2008-08-31

    Black-oil, compositional and thermal simulators have been developed to address different physical processes in reservoir simulation. A number of different types of discretization methods have also been proposed to address issues related to representing the complex reservoir geometry. These methods are more significant for fractured reservoirs where the geometry can be particularly challenging. In this project, a general modular framework for reservoir simulation was developed, wherein the physical models were efficiently decoupled from the discretization methods. This made it possible to couple any discretization method with different physical models. Oil characterization methods are becoming increasingly sophisticated, and it is possible to construct geologically constrained models of faulted/fractured reservoirs. Discrete Fracture Network (DFN) simulation provides the option of performing multiphase calculations on spatially explicit, geologically feasible fracture sets. Multiphase DFN simulations of and sensitivity studies on a wide variety of fracture networks created using fracture creation/simulation programs was undertaken in the first part of this project. This involved creating interfaces to seamlessly convert the fracture characterization information into simulator input, grid the complex geometry, perform the simulations, and analyze and visualize results. Benchmarking and comparison with conventional simulators was also a component of this work. After demonstration of the fact that multiphase simulations can be carried out on complex fracture networks, quantitative effects of the heterogeneity of fracture properties were evaluated. Reservoirs are populated with fractures of several different scales and properties. A multiscale fracture modeling study was undertaken and the effects of heterogeneity and storage on water displacement dynamics in fractured basements were investigated. In gravity-dominated systems, more oil could be recovered at a given pore

  4. Making Problem-Solving Simulations More Realistic.

    ERIC Educational Resources Information Center

    Cotton, Samuel E.

    2002-01-01

    Many problem-solving activities include mathematical principles but students do not use them during the design and experimentation phases before creating a prototype or product. Restricting the amount and/or type of materials available to students will require them to calculate and requisition the materials needed. (JOW)

  5. Three-dimensional geomechanical simulation of reservoir compaction and implications for well failures in the Belridge diatomite

    SciTech Connect

    Fredrich, J.T.; Argueello, J.G.; Thorne, B.J.; Wawersik, W.R.

    1996-12-31

    This paper describes an integrated geomechanics analysis of well casing damage induced by compaction of the diatomite reservoir at the Belridge Field, California. Historical data from the five field operators were compiled and analyzed to determine correlations between production, injection, subsidence, and well failures. The results of this analysis were used to develop a three-dimensional geomechanical model of South Belridge, Section 33 to examine the diatomite reservoir and overburden response to production and injection at the interwell scale and to evaluate potential well failure mechanisms. The time-dependent reservoir pressure field was derived from a three-dimensional finite difference reservoir simulation and used as input to three-dimensional non-linear finite element geomechanical simulations. The reservoir simulation included approximately 200 wells and covered 18 years of production and injection. The geomechanical simulation contained 437,100 nodes and 374,130 elements with the overburden and reservoir discretized into 13 layers with independent material properties. The results reveal the evolution of the subsurface stress and displacement fields with production and injection and suggest strategies for reducing the occurrence of well casing damage.

  6. Three-dimensional geomechanical simulation of reservoir compaction and implications for well failures in the Belridge diatomite

    SciTech Connect

    Fredrich, J.T.; Argueello, J.G.; Thorne, B.J.; Wawersik, W.R. |

    1996-11-01

    This paper describes an integrated geomechanics analysis of well casing damage induced by compaction of the diatomite reservoir at the Belridge Field, California. Historical data from the five field operators were compiled and analyzed to determine correlations between production, injection, subsidence, and well failures. The results of this analysis were used to develop a three-dimensional geomechanical model of South Belridge, Section 33 to examine the diatomite reservoir and overburden response to production and injection at the interwell scale and to evaluate potential well failure mechanisms. The time-dependent reservoir pressure field was derived from a three-dimensional finite difference reservoir simulation and used as input to three-dimensional non-linear finite element geomechanical simulations. The reservoir simulation included -200 wells and covered 18 years of production and injection. The geomechanical simulation contained 437,100 nodes and 374,130 elements with the overburden and reservoir discretized into 13 layers with independent material properties. The results reveal the evolution of the subsurface stress and displacement fields with production and injection and suggest strategies for reducing the occurrence of well casing damage.

  7. Integrated 3-D Ground-Penetrating Radar, Outcrop, and Borehole Data Applied to Reservoir Characterization and Flow Simulation

    SciTech Connect

    George McMechan; Rucsandra Corbeanu; Craig Forster; Kristian Soegaard; Xiaoxian Zeng; Carlos Aiken; Robert Szerbiak; Janok Bhattacharya; Michael Wizevich; Xueming Xu; Stephen Snelgrove; Karen Roche; Siang Joo Lim; Djuro Navakovic; Christopher White; Laura Crossey; Deming Wang; John Thurmond; William Hammon III; Mamadou BAlde; Ari Menitove

    2001-08-31

    OAK-B135 (IPLD Cleared) Existing reservoir models are based on 2-D outcrop studies; 3-D aspects are inferred from correlation between wells, and so are inadequately constrained for reservoir simulations. To overcome these deficiencies, we initiated a multidimensional characterization of reservoir analogs in the Cretaceous Ferron Sandstone in Utah. The study was conducted at two sites (Corbula Gulch and Coyote Basin); results from both sites are contained in this report. Detailed sedimentary facies maps of cliff faces define the geometry and distribution of potential reservoir flow units, barriers and baffles at the outcrop. High resolution 2-D and 3-D ground-penetrating radar (GPR) images extend these reservoir characteristics into 3-D, to allow development of realistic 3-D reservoir models. Models use geometric information from the mapping and the GPR data, petrophysical data from surface and cliff-face outcrops, lab analyses of outcrop and core samples, and petrography. The measurements are all integrated into a single coordinate system using GPS and laser mapping of the main sedimentological features and boundaries.The final step is analysis of results of 3-D fluid flow modeling to demonstrate applicability of our reservoir analog studies to well siting and reservoir engineering for maximization of hydrocarbon production. The main goals of the project are achieved. These are the construction of a deterministic 3-D reservoir analog model from a variety of geophysical and geologic measurements at the field sites, integrating these into comprehensive petrophysical models, and flow simulations through these models. This unique approach represents a significant advance in characterization and use of reservoir analogs.

  8. Intergrated 3-D Ground-Penetrating Radar,Outcrop,and Boreholoe Data Applied to Reservoir Characterization and Flow Simulation.

    SciTech Connect

    McMechan et al.

    2001-08-31

    Existing reservoir models are based on 2-D outcrop;3-D aspects are inferred from correlation between wells,and so are inadequately constrained for reservoir simulations. To overcome these deficiencies, we initiated a multidimensional characterization of reservoir analogs in the Cretaceous Ferron Sandstone in Utah.The study was conducted at two sites(Corbula Gulch Coyote Basin); results from both sites are contained in this report. Detailed sedimentary facies maps of cliff faces define the geometry and distribution of potential reservoir flow units, barriers and baffles at the outcrop. High resolution 2-D and 3-D ground penetrating radar(GPR) images extend these reservoir characteristics into 3-D to allow development of realistic 3-D reservoir models. Models use geometric information from the mapping and the GPR data, petrophysical data from surface and cliff-face outcrops, lab analyses of outcrop and core samples, and petrography. The measurements are all integrated into a single coordinate system using GPS and laser mapping of the main sedimentologic features and boundaries. The final step is analysis of results of 3-D fluid flow modeling to demonstrate applicability of our reservoir analog studies to well siting and reservoir engineering for maximization of hydrocarbon production. The main goals of this project are achieved. These are the construction of a deterministic 3-D reservoir analog model from a variety of geophysical and geologic measurements at the field sites, integrating these into comprehensive petrophysical models, and flow simulation through these models. This unique approach represents a significant advance in characterization and use of reservoir analogs. To data,the team has presented five papers at GSA and AAPG meetings produced a technical manual, and completed 15 technical papers. The latter are the main content of this final report. In addition,the project became part of 5 PhD dissertations, 3 MS theses,and two senior undergraduate research

  9. Numerical Simulation of Subsurface Transport and Groundwater Impacts from Hydraulic Fracturing of Tight/Shale Gas Reservoirs

    NASA Astrophysics Data System (ADS)

    Reagan, M. T.; Moridis, G. J.; Keen, N. D.

    2014-12-01

    The use of reservoir stimulation techniques, such as hydraulic fracturing, has grown tremendously over the last decade, and concerns have arisen that reservoir stimulation creates environmental threats through the creation of permeable pathways that could connect the stimulated reservoir to shallower groundwater aquifers. This study investigates, by numerical simulation, gas and water transport between a deeper tight-gas reservoir and a shallower overlying groundwater aquifer following hydraulic fracturing operations, assuming that the formation of a connecting pathway has already occurred. We focus on two general transport scenarios: 1) communication between the reservoir and aquifer via a connecting fracture or fault and 2) communication via a deteriorated, preexisting nearby well. The simulations explore a range of permeabilities and geometries over time scales, and evaluate the mechanisms and factors that could lead to the escape of gas or reservoir fluid and the contamination of groundwater resources. We also examine the effects of overpressured reservoirs, and explore long-term transport processes as part of a continuing study. We conclude that the key factors driving short-term transport of gas include high permeability for the connecting pathway and the overall volume of the connecting feature. Gas production from the reservoir via a horizontal well is likely to mitigate release through the reduction of available free gas and the lowering of reservoir pressure. We also find that fractured tight-gas reservoirs are unlikely to act as a continuing source of large volumes of migrating gas, and incidents of gas escape are likely to be limited in duration and scope. Reliable field and laboratory data must be acquired to constrain the factors and determine the likelihood of these outcomes.

  10. Simulation studies to evaluate the effect of fracture closure on the performance of naturally fractured reservoirs. Annual report

    SciTech Connect

    Not Available

    1992-11-01

    The second year of this three-year research program to evaluate the effect of fracture closure on the recovery of oil and gas from naturally fractured reservoirs has been completed. The overall objectives of the study are to: (1) evaluate the reservoir conditions where fracture closure is significant, and (2) evaluate innovative fluid injection techniques capable of maintaining pressure within the reservoir. Simulation studies have been conducted with a dual porosity simulator capable of simulating the performance of vertical and horizontal wells. Each simulation model has been initialized with properties typical of the Austin Chalk reservoir in Pearsall Field, Texas. During year one, simulations of both vertical and horizontal well performance were made assuming that fracture permeability was insensitive to pressure charge. The results confirmed that horizontal wells could increase both rate of oil recovery and total oil recovery from naturally fractured reservoirs. During the second year the performances of the same vertical and horizontal wells were evaluated with the assumption that fracture permeability was a function of reservoir pressure. This required repetition of most of the natural depletion cases simulated in year one while invoking the pressure-sensitive fracture permeability option. To investigate sensitivity to in situ stress, two stress conditions were simulated for each primary variable. The water injection cases, begun in year one, were extended to include most of the reservoir parameters investigated for natural depletion, including fracture permeability as a function of net stress and the use of horizontal wells. The results thus far confirm that pressure-sensitive fractures degrade well performance and that the degradation is reduced by water injection pressure maintenance. Furthermore, oil recovery can be significantly increased by water injection pressure maintenance.

  11. Numerical simulation of groundwater movement and managed aquifer recharge from Sand Hollow Reservoir, Hurricane Bench area, Washington County, Utah

    USGS Publications Warehouse

    Marston, Thomas M.; Heilweil, Victor M.

    2012-01-01

    The Hurricane Bench area of Washington County, Utah, is a 70 square-mile area extending south from the Virgin River and encompassing Sand Hollow basin. Sand Hollow Reservoir, located on Hurricane Bench, was completed in March 2002 and is operated primarily as a managed aquifer recharge project by the Washington County Water Conservancy District. The reservoir is situated on a thick sequence of the Navajo Sandstone and Kayenta Formation. Total recharge to the underlying Navajo aquifer from the reservoir was about 86,000 acre-feet from 2002 to 2009. Natural recharge as infiltration of precipitation was approximately 2,100 acre-feet per year for the same period. Discharge occurs as seepage to the Virgin River, municipal and irrigation well withdrawals, and seepage to drains at the base of reservoir dams. Within the Hurricane Bench area, unconfined groundwater-flow conditions generally exist throughout the Navajo Sandstone. Navajo Sandstone hydraulic-conductivity values from regional aquifer testing range from 0.8 to 32 feet per day. The large variability in hydraulic conductivity is attributed to bedrock fractures that trend north-northeast across the study area. A numerical groundwater-flow model was developed to simulate groundwater movement in the Hurricane Bench area and to simulate the movement of managed aquifer recharge from Sand Hollow Reservoir through the groundwater system. The model was calibrated to combined steady- and transient-state conditions. The steady-state portion of the simulation was developed and calibrated by using hydrologic data that represented average conditions for 1975. The transient-state portion of the simulation was developed and calibrated by using hydrologic data collected from 1976 to 2009. Areally, the model grid was 98 rows by 76 columns with a variable cell size ranging from about 1.5 to 25 acres. Smaller cells were used to represent the reservoir to accurately simulate the reservoir bathymetry and nearby monitoring wells

  12. Three-dimensional numerical reservoir simulation of the EGS Demonstration Project at The Geysers geothermal field

    NASA Astrophysics Data System (ADS)

    Borgia, Andrea; Rutqvist, Jonny; Oldenburg, Curt M.; Hutchings, Lawrence; Garcia, Julio; Walters, Mark; Hartline, Craig; Jeanne, Pierre; Dobson, Patrick; Boyle, Katie

    2013-04-01

    -isothermal porous media numerical flow simulator in order to model the evolution and injection-related operational dynamics of The Geysers geothermal field. At the bottom of the domain in the felsite, we impose a constant temperature, constant saturation, low-permeability boundary. Laterally we set no-flow boundaries (no mass or heat flow), while at the top we use a fully aqueous-phase-saturated constant atmospheric pressure boundary condition. We compute initial conditions for two different conceptual models. The first conceptual model has two phases (gas and aqueous) with decreasing proportions of gas from the steam zone downward; the second model has dry steam all the way from the steam zone to the bottom. The first may be more similar to a pre-exploitation condition, before production reduced pressure and dried out the system, while the second is calibrated to the pressure and temperature actually measured in the reservoir today. Our preliminary results are in reasonable agreement with the pressure monitoring at Prati State 31. These results will be used in hydrogeomechanical modeling to plan, design, and validate the effects of injection in the system.

  13. Simulation studies of air transport operational problems

    NASA Technical Reports Server (NTRS)

    Lauber, J. K.; Billings, C. E.; Stevenson, J. E.; Ruffell-Smith, H. P.; Cooper, G. E.

    1976-01-01

    An experimental evaluation of the monitored approach procedure for conducting low visibility instrument approaches is described. Four airline crews each flew 16 approaches using the monitored procedure and 16 using a modified standard procedure in a DC-10 simulator under various conditions of visibility, wind shear and turbulence, and radar vectoring scenarios. In terms of system measures of aircrew performance, no major differences were found. Pilot opinion data indicate that there are some desirable characteristics of the monitored procedure, particularly with reference to the increased role of the flight engineer in conducting low visibility approaches. Rationale for developing approach procedures is discussed.

  14. Finite element method for simulating coupled thermo-hydro-mechanical processes in discretely fractured porous media and application to enhanced geothermal reservoir analysis

    NASA Astrophysics Data System (ADS)

    Watanabe, N.; Wong, L.; Bloecher, G.; Cacace, M.; Kolditz, O.

    2012-12-01

    We present our recent development of the finite element method (FEM) for simulating coupled thermo-hydro-mechanical (THM) processes in discretely fractured porous media and an application to geothermal reservoir modeling for the research test site Gross Schoenebeck in Germany operated by the GFZ German Research Centre for Geosciences. Numerical analysis of multi-physics problems in fractured rocks is important for various geotechnical applications. In particular for enhanced geothermal reservoirs where induced fractures and possibly natural fault systems dominate the system behavior, explicit modeling of those characteristic fractures (i.e. discrete fracture models) is essential to get more detailed understanding of in-situ processes and reliable estimations of heat extraction from those deep reservoirs. However, as fractures are mechanical discontinuities, it is difficult to solve the problems using continuity based numerical methods such as the FEM. Currently, equivalent porous medium or multiple continuum model approaches are often only the way to model fractured rocks with the FEM. The authors have recently developed lower-dimensional interface elements (LIEs) for modeling mechanics-involved coupled processes with pre-existing fractures (Watanabe et al. 2012 IJNME). The method does not require any double nodes unlike conventional interface elements. Moreover, for coupled problems, the approach allows for the use of a single mesh for both mechanical and other related processes such as flow and transport. All the code developments have been carried out within the scientific open source project OpenGeoSys (www.opengeosys.net) (Kolditz et al. 2012 EES). Using both traditional and new simulation techniques, a geothermal reservoir model for the research test site Gross Schoenebeck has been developed. Unstructured meshing of the complex faulted reservoir including both rock matrix and fracture elements has been conducted using recently developed automatic

  15. Simulating Thermal-Hydrologic-Mechanical-Chemical Evolution Surrounding Fluid Injection in a Fractured Porous Geothermal Reservoir

    NASA Astrophysics Data System (ADS)

    Taron, J.; Min, K.; Elsworth, D.

    2006-12-01

    Computational analysis is conducted on the coupled thermal-hydrologic-mechanical-chemical (THMC) behavior of a stimulated EGS geothermal reservoir. Numerical analyses utilize a newly developed simulator capable of examining THMC processes in fractured porous geologic media. The simulator links the thermal-hydrologic- chemical (THC) computational code TOUGHREACT with the mechanical (M) capability of FLAC3D, where the response of pore fluid pressure to mechanical disturbance is treated as an undrained system and mineral precipitation/dissolution generates porosity and permeability change within each dual-permeability continuum. Non-linear permeability response to thermal-hydrologic-mechanical (THM) mechanisms is accommodated via embryonic mechanical and transport constitutive laws, and is considered to act in union with permeability changes associated with the removal or addition of minerals within the system. This construct is applied to the geometry of an injector-withdrawal doublet within the Coso Geothermal field, where in situ stress conditions, thermal state, and mineralogical composition at 3000m depth are extracted from recorded field data. Initial results for feasible parametric settings show that permeability reduction in the vicinity of a cool (80°C) injection well may be significant, within an order of magnitude, and accompanied by large (MPa) changes in the stress field throughout the reservoir for imposed boundary conditions of constant stress.

  16. 3D hybrid tectono-stochastic modeling of naturally fractured reservoir: Application of finite element method and stochastic simulation technique

    NASA Astrophysics Data System (ADS)

    Gholizadeh Doonechaly, N.; Rahman, S. S.

    2012-05-01

    Simulation of naturally fractured reservoirs offers significant challenges due to the lack of a methodology that can utilize field data. To date several methods have been proposed by authors to characterize naturally fractured reservoirs. Among them is the unfolding/folding method which offers some degree of accuracy in estimating the probability of the existence of fractures in a reservoir. Also there are statistical approaches which integrate all levels of field data to simulate the fracture network. This approach, however, is dependent on the availability of data sources, such as seismic attributes, core descriptions, well logs, etc. which often make it difficult to obtain field wide. In this study a hybrid tectono-stochastic simulation is proposed to characterize a naturally fractured reservoir. A finite element based model is used to simulate the tectonic event of folding and unfolding of a geological structure. A nested neuro-stochastic technique is used to develop the inter-relationship between the data and at the same time it utilizes the sequential Gaussian approach to analyze field data along with fracture probability data. This approach has the ability to overcome commonly experienced discontinuity of the data in both horizontal and vertical directions. This hybrid technique is used to generate a discrete fracture network of a specific Australian gas reservoir, Palm Valley in the Northern Territory. Results of this study have significant benefit in accurately describing fluid flow simulation and well placement for maximal hydrocarbon recovery.

  17. Role of reservoir simulation in development and management of complexly-faulted, multiple-reservoir Dulang field, offshore Malaysia: Holistic strategies

    SciTech Connect

    Sonrexa, K.; Aziz, A.; Solomon, G.J.; Bandal, M.; Embong, M.K.; Wahir, H.

    1995-10-01

    The Dulang field, discovered in 1981, is a major oil filed located offshore Malaysia in the Malay Basin. The Dulang Unit Area constitutes the central part of this exceedingly heterogeneous field. The Unit Area consists of 19 stacked shaly sandstone reservoirs which are divided into about 90 compartments with multiple fluid contacts owing to severe faulting. Current estimated put the Original-Oil-In-Place (OOIP) in the neighborhood of 700 million stock tank barrels (MMSTB). Production commenced in March 1991 and the current production is more than 50,000 barrels of oil per day (BOPD). In addition to other more conventional means, reservoir simulation has been employed form the very start as a vital component of the overall strategy to develop and manage this challenging field. More than 10 modeling studies have been completed by Petronas Carigali Sdn. Bhd. (Carigali) at various times during the short life of this field thus far. To add to that, Esso Production Malaysia Inc. (EPMI) has simultaneously conducted a number of independent studies. These studies have dealt with undersaturated compartments as well as those with small and large gas caps. They have paved the way for improved reservoir characterization, optimum development planning and prudent production practices. This paper discusses the modeling approaches and highlights the crucial role these studies have played on an ongoing basis in the development and management of the complexly-faulted, multi-reservoir Dulang Unit Area.

  18. Geothermal reservoir engineering computer code comparison and validation

    SciTech Connect

    Faust, C.R.; Mercer, J.W.; Miller, W.J.

    1980-11-12

    The results of computer simulations for a set of six problems typical of geothermal reservoir engineering applications are presented. These results are compared to those obtained by others using similar geothermal reservoir simulators on the same problem set. The purpose of this code comparison is to check the performance of participating codes on a set of typical reservoir problems. The results provide a measure of the validity and appropriateness of the simulators in terms of major assumptions, governing equations, numerical accuracy, and computational procedures. A description is given of the general reservoir simulator - its major assumptions, mathematical formulation, and numerical techniques. Following the description of the model is the presentation of the results for the six problems. Included with the results for each problem is a discussion of the results; problem descriptions and result tabulations are included in appendixes. Each of the six problems specified in the contract was successfully simulated. (MHR)

  19. Efficient algorithms for distributed simulation and related problems

    SciTech Connect

    Kumar, D.

    1987-01-01

    This thesis presents efficient algorithms for distributed simulation, and for the related problems of termination detection and sequential simulation. Distributed simulation algorithms applicable to the simulation of special classes of systems, such that almost no overhead messages are required are presented. By contrast, previous distributed simulation algorithms, although applicable to the general class of any discrete-event system, usually require too many overhead messages. First, a simple distributed simulation algorithm is defined with nearly zero overhead messages for simulating feedforward systems. An approximate method is developed to predict its performance in simulating a class of feedforward-queuing networks. Performance of the scheme is evaluated in simulating specific subclasses of these queuing networks. It is shown that the scheme offers a high performance for serial-parallel networks. Next, another distributed simulation scheme is defined for a class of distributed systems whose topologies may have cycles. One important problem in devising distributed simulation algorithms is that of efficient detection of termination. With this in mind, a class of termination-detection algorithms using markers is devised. Finally, a new sequential simulation algorithm is developed, based on a distributed one. This algorithm often reduces the event-list manipulations of traditional-event list-driven simulation.

  20. Hydrothermal simulation of a fractured carbonate reservoir in southern Italy and automated detections of optimal positions for geothermal doublet installations

    NASA Astrophysics Data System (ADS)

    Niederau, Jan; Gomez, Sergio; Ebigbo, Anozie; Inversi, Barbara; Marquart, Gabriele; Scrocca, Davide

    2015-04-01

    In this work, we present the results of hydrothermal simulations for assessing the geothermal potential of a fractured carbonate reservoir in Campania (Guardia Lombardi). Local surface heat flows of up to 90 mW/m² suggest that this area is a potential medium-enthalpy geothermal reservoir. The targeted reservoir rocks are fractured shallow-water carbonates (Jurassic to Cretaceous) of the Apulia Platform. During the Apennine orogeny, those carbonates were affected by at least two tectonic phases: Thrust-related folding of the carbonate platform due to compression followed by extension which caused major normal faulting. Based on seismic interpretation, a discretized structural model is set up, comprising the reservoir unit and the overlying sedimentary cover. The model comprises an area of 42 km × 28 km and extends to a depth of about six kilometers. Results of calibrated hydrothermal reservoir simulations suggest that free convection occurs in some parts of the reservoir. For assessing optimal locations for potential hydrothermal doublet systems, a tool was developed which uses the results of the reservoir simulationsin combination with predefined constraints. Those constraints or minimum requirements consider: a) minimum temperature for operating the doublet system, b) minimum matrix permeability allowing for a pumping rate of 40 L/s, and c) social constraints (location of cities or conservation areas, where the construction of a potential geothermal energy plant would be problematic). The optimization tool ranks possible doublet system locations by evaluating an objective function for the minimum requirements. Those locations are further used to extract smaller models from the big reservoir model and simulate the operation of a hypothetical geothermal doublet system. By assessing the optimized results, an optimal location of a geothermal energy plant would produce water with a temperature of 163 °C from a depth of almost 4 km.

  1. A virtual company concept for reservoir management

    SciTech Connect

    Martin, F.D.; Kendall, R.P.; Whitney, E.M.

    1998-12-31

    This paper describes how reservoir management problems were pursued with a virtual company concept via the Internet and World Wide Web. The focus of the paper is on the implementation of virtual asset management teams that were assembled with small independent oil companies. The paper highlights the mechanics of how the virtual team transferred data and interpretations, evaluated geological models of complex reservoirs, and used results of simulation studies to analyze various reservoir management strategies.

  2. Problems Caused by Microbes and Treatment Strategies Monitoring and Preventing Reservoir Souring Using Molecular Microbiological Methods (MMM)

    NASA Astrophysics Data System (ADS)

    Gittel, Antje

    The injection of seawater during the process of secondary oil recovery in offshore oilfields supplies huge amounts of sulphate to the prokaryotic reservoir communities. Together with the presence of oil organics and their degradation products as electron donors, this facilitates the enrichment and growth of sulphate-reducing prokaryotes (SRP) in the reservoir, as well as in pipings and top-side installations (Sunde and Torsvik, 2005; Vance and Thrasher, 2005). The activity of SRP causes severe economic problems due to the reactivity and toxicity of the produced hydrogen sulphide (H2S), one of the major problems being reservoir souring. Besides the use of broad-spectrum biocides or inhibitors for sulphate reduction, the addition of nitrate effectively decreased the net production of H2S in model column studies (Myhr et al., 2002; Hubert et al., 2005; Dunsmore et al., 2006) and field trials (Telang et al., 1997; Bødtker et al., 2008). The mechanisms by which nitrate addition might affect souring control are (i) the stimulation of heterotrophic nitrate-reducing bacteria (hNRB) that outcompete SRP for electron donors, (ii) the activity of nitrate-reducing, sulphide-oxidising bacteria (NR-SOB), and (iii) the inhibition of SRP by the production of nitrite and nitrous oxides (Sunde and Torsvik, 2005; Hubert and Voordouw, 2007).

  3. Mineral trapping of CO2 in operated geothermal reservoirs - Numerical simulations on various scales

    NASA Astrophysics Data System (ADS)

    Kühn, Michael; Stanjek, Helge; Peiffer, Stefan; Clauser, Christoph

    2013-04-01

    A novel approach to store CO2 not only by hydrodynamic trapping within a reservoir, but to convert dissolved CO2 into the geochemically more stable form of calcite in a reaction with calcium obtained from dissolution of sulphates and alkalinity from feldspars or fly ashes is described here. The presentation gives answers to the key questions: • Where are potential geothermal reservoirs with anhydrite abundant? • Does the transfer of anhydrite into calcite work at all and what are the reaction rates? • What are probable alkalinity sources and how fast are they available? Numerical simulation is a means to quantify the entire process of CO2 storage and to deepen the understanding of the detailed chemical processes. We performed numerical simulations on multiple scales. The relevant scales reach from the micro or thin section scale (ca. 1 cm) to the reservoir scale (ca. 10 km). The idea is to provide constraints for smaller scale models from the larger scale and derive functionality from smaller scale models of processes which cannot be resolved in larger scale models, due to restrictions of discretization of the applied numerical mesh. With regard to the 3 questions above we can conclude that the combination of CO2 storage and geothermal energy production is generally feasible because candidate sites are available, anhydrite is transformable into calcite and alkalinity can be provided by fly ashes (Back et al. 2010) or even in-situ (Kühn and Clauser 2006). Based on our laboratory investigations and numerical studies we are able to estimate the storage potential for mineral trapping of CO2 in geothermal reservoirs (Kühn et al. 2009). On the one hand the maximum is unfortunately less than a million tons over the life time of a geothermal heating plant. On the other hand significant storage capacities are available in geological formations with regard to hydrodynamic trapping for millions of tonnes of carbon dioxide. This is why under the current circumstances

  4. Development and use of a 3D simulator in evaluating reservoir development options: Fateh Mishrif, offshore Dubai

    SciTech Connect

    Crick, J.; Singh, S.P.

    1985-01-01

    Fateh Mishrif, one of the largest reservoir offshore Dubai, is a thick carbonate reefal development of which the centre has been completely eroded. During the past fifteen years production has been through natural depletion, peripheral water injection and gas lift. This paper describes the construction of a 3D reservoir model for the Fateh Mishrif reservoir and its application in evaluating further reservoir developments plans. Significant reservoir factors which influence fluid flow are a field wide tar-mat, faulting and the central pinch out. The initial fluids in place description is affected by two oil-water contacts across the field. The model simulates these features accurately. The prediction of future oil rates is based on an efficient gas lift optimization technique on a well by well basis which also provides a forecast of future gas lift requirements for the reservoir. The constraints on total available gas for gas lift and water injection rates, both field wide or regionally, can be modelled adequately. The model has been used for evaluating the expansion of the existing waterflood and future infill drilling potential. An important result of this study has been that the history matching phase has highlighted the need for further data gathering in specific parts of the reservoir.

  5. Parallel computation for reservoir thermal simulation: An overlapping domain decomposition approach

    NASA Astrophysics Data System (ADS)

    Wang, Zhongxiao

    2005-11-01

    In this dissertation, we are involved in parallel computing for the thermal simulation of multicomponent, multiphase fluid flow in petroleum reservoirs. We report the development and applications of such a simulator. Unlike many efforts made to parallelize locally the solver of a linear equations system which affects the performance the most, this research takes a global parallelization strategy by decomposing the computational domain into smaller subdomains. This dissertation addresses the domain decomposition techniques and, based on the comparison, adopts an overlapping domain decomposition method. This global parallelization method hands over each subdomain to a single processor of the parallel computer to process. Communication is required when handling overlapping regions between subdomains. For this purpose, MPI (message passing interface) is used for data communication and communication control. A physical and mathematical model is introduced for the reservoir thermal simulation. Numerical tests on two sets of industrial data of practical oilfields indicate that this model and the parallel implementation match the history data accurately. Therefore, we expect to use both the model and the parallel code to predict oil production and guide the design, implementation and real-time fine tuning of new well operating schemes. A new adaptive mechanism to synchronize processes on different processors has been introduced, which not only ensures the computational accuracy but also improves the time performance. To accelerate the convergence rate of iterative solution of the large linear equations systems derived from the discretization of governing equations of our physical and mathematical model in space and time, we adopt the ORTHOMIN method in conjunction with an incomplete LU factorization preconditioning technique. Important improvements have been made in both ORTHOMIN method and incomplete LU factorization in order to enhance time performance without affecting

  6. Analysis of induced seismicity and heat transfer in geothermal reservoirs by coupled simulation

    NASA Astrophysics Data System (ADS)

    Gan, Quan

    . Conversely, at high relative non-dimensional flow rates the propagating pressure pulse is larger and migrates more quickly through the reservoir but the thermal drawdown is uniform across the reservoir and without the presence of a distinct thermal front, and less capable of triggering late-stage seismicity. In Chapter 2 we develop a dimensionless model to predict the thermal drawdown response, and quantify the relationship between the timing and magnitude of late stage seismic event and the induced thermal stress from thermal drawdown. We evaluate the uniformity of thermal drawdown as a function of a dimensionless flow rate QD that scales with fracture spacing s( m), injection rate q (kg/s), and the distance between the injector and the target point L* ( Qd ∝ qs2 / L*). By assuming the dominant heat transfer by heat conduction within the fractured medium, this model is either capable to predict the timing of induced seismicity by the thermal stress by the analytical formula. Due to the significant influence of fracture network geometry in heat transfer and induced seismicity, a discrete fracture network model is developed (Chapter 3) to couple stress and fluid flow in a discontinuous fractured mass represented as a continuum by coupling the continuum simulator TF_FLAC 3D with cell-by-cell discontinuum laws for deformation and flow. Both equivalent medium crack and permeability tensor approaches are employed to characterize preexisting discrete fractures. The evolution of fracture permeability accommodates stress-dependent aperture under different stress states, including normal closure, shear dilation, and for fracture walls out of contact under tensile loading. This discrete fracture network model is applied (Chapter 4) in a generic reservoir with an initial permeability in the range of 10-17 to 10-16 m2, fracture density of ~0.09m -1 and fractures oriented such that either none, one, or both sets of fractures are critically stressed. For a given reservoir with a pre

  7. Image analysis for Validation of Simulations of Fluid Mix Problem

    SciTech Connect

    Kamath, C; Miller, P

    2007-01-10

    As computer simulations gain acceptance for the modeling of complex physical phenomena, there is an increasing need to validate these simulation codes by comparing them to experiments. Currently, this is done qualitatively, using a visual approach. This is obviously very subjective and more quantitative metrics are needed, especially to identify simulations which are closer to experiments than other simulations. In this paper, we show how image processing techniques can be effectively used in such comparisons. Using an example from the problem of mixing of two fluids, we show that we can quantitatively compare experimental and simulation images by extracting higher level features to characterize the objects in the images.

  8. The Diffusion Simulator - Teaching Geomorphic and Geologic Problems Visually.

    ERIC Educational Resources Information Center

    Gilbert, R.

    1979-01-01

    Describes a simple hydraulic simulator based on more complex models long used by engineers to develop approximate solutions. It allows students to visualize non-steady transfer, to apply a model to solve a problem, and to compare experimentally simulated information with calculated values. (Author/MA)

  9. THE ITA PROBLEM: A Ready-to-Use Simulation

    ERIC Educational Resources Information Center

    Halleck, Gene B.

    2008-01-01

    This simulation probes what is known as the "foreign" teaching assistant problem. The "problem" can be found at large state universities in the United States where international graduate students are required to earn their scholarships by teaching undergraduate courses and arises because of a combination of issues, including the oral proficiency…

  10. A Problem-Solving Simulation for Classical Genetics.

    ERIC Educational Resources Information Center

    Stewart, Jim

    1996-01-01

    Presents a genetics simulation that uses a noncomputer version of BioQUEST's Genetics Construction Kit. Provides students with opportunities for solving realistic genetics problems and allows them to experience how classical geneticists think. Discusses the problem's search space, thinking qualitatively, hypothesis generation and testing, thinking…

  11. A Simulated Research Problem for Undergraduate Metamorphic Petrology.

    ERIC Educational Resources Information Center

    Amenta, Roddy V.

    1984-01-01

    Presents a laboratory problem in metamorphic petrology designed to simulate a research experience. The problem deals with data on scales ranging from a geologic map to hand specimens to thin sections. Student analysis includes identifying metamorphic index minerals, locating their isograds on the map, and determining the folding sequence. (BC)

  12. Multiscale Snow/Icemelt Discharge Simulations into Alpine Reservoirs: adding Glacier Dynamics to a Hydrological Model

    NASA Astrophysics Data System (ADS)

    Schueller, Felix; Förster, Kristian; Hanzer, Florian; Huttenlau, Matthias; Marzeion, Ben; Strasser, Ulrich; Achleitner, Stefan; Kirnbauer, Robert

    2015-04-01

    Glacier and snow runoff in high alpine regions is an essential process in hydrological research for its high relevance on lower altitude areas and hydro-power generation. MUSICALS II (Multiscale Snow/Icemelt Discharge Simulations into Alpine Reservoirs) seeks to identify and quantify water availability and runoff in alpine headwater catchments. The focus is on future changes due to glacier retreat, altering the multi-day and seasonal runoff available for hydropower operations. Our aim is to investigate and improve runoff forecasts by coupling the semi-distributed hydrological model HQSim with a simple glacier evolution model. The glacier model MMBM (Marzeion Mass Balance Model) with its statistical nature allows for fast modelling of the dynamical properties of glaciers. We present the design of the coupled hydrological application for different hydro power headwater catchments in Tyrol. The capabilities of the glacier model to simulate the selected glaciers is shown. Simulated discharge with the original and the coupled model are compared to downstream gauge measurements. Using the multi-objective optimization algorithm AMALGAM (A Multi-ALgorithm, Genetically Adaptive Multiobjective model), we optimize the glacier module parameters fully automatically. The results show the improvements in runoff modelling for past periods, when altering of glaciated catchment parts is considered. This indicates consideration of this process is mandatory for simulating future developments.

  13. Non-equilibrium simulation of CH4 production through the depressurization method from gas hydrate reservoirs

    NASA Astrophysics Data System (ADS)

    Qorbani, Khadijeh; Kvamme, Bjørn

    2016-04-01

    Natural gas hydrates (NGHs) in nature are formed from various hydrate formers (i.e. aqueous, gas, and adsorbed phases). As a result, due to Gibbs phase rule and the combined first and second laws of thermodynamics CH4-hydrate cannot reach thermodynamic equilibrium in real reservoir conditions. CH4 is the dominant component in NGH reservoirs. It is formed as a result of biogenic degradation of biological material in the upper few hundred meters of subsurface. It has been estimated that the amount of fuel-gas reserve in NGHs exceed the total amount of fossil fuel explored until today. Thus, these reservoirs have the potential to satisfy the energy requirements of the future. However, released CH4 from dissociated NGHs could find its way to the atmosphere and it is a far more aggressive greenhouse gas than CO2, even though its life-time is shorter. Lack of reliable field data makes it difficult to predict the production potential, as well as safety of CH4 production from NGHs. Computer simulations can be used as a tool to investigate CH4 production through different scenarios. Most hydrate simulators within academia and industry treat hydrate phase transitions as an equilibrium process and those which employ the kinetic approach utilize simple laboratory data in their models. Furthermore, it is typical to utilize a limited thermodynamic description where only temperature and pressure projections are considered. Another widely used simplification is to assume only a single route for the hydrate phase transitions. The non-equilibrium nature of hydrate indicates a need for proper kinetic models to describe hydrate dissociation and reformation in the reservoir with respect to thermodynamics variables, CH4 mole-fraction, pressure and temperature. The RetrasoCodeBright (RCB) hydrate simulator has previously been extended to model CH4-hydrate dissociation towards CH4 gas and water. CH4-hydrate is added to the RCB data-base as a pseudo mineral. Phase transitions are treated

  14. Coal Chemistry for Mechanical Engineers: From Macromolecular Thermodynamics to Reservoir Simulation

    SciTech Connect

    Romanov, V.

    2007-05-01

    In pilot trials and commercial scale field demonstrations of CO2 storage in coal seams, quite often unexpected problems with coal swelling around injector and reducing injection efficiency (e.g., Allison unit in the San Juan Basin, RECOPOL in Poland, Hokkaido project in Japan, etc.) can stall or even terminate the site development. To avoid the costly mistakes with the prospective site evaluation, the state of the art in reservoir modeling needs to be improved by taking into account coal properties at the macromolecular level. The current models are based on the rock mechanics, which ignores decades of experimental and theoretical studies of interaction between coal and injected fluids. A pseudopolymer approach is introduced to the modelers as a viable alternative, especially, at medium to high fluid pressures. Further, it is discussed how the thermodynamics of CO2 dissolution in the macromolecular network of the coal matrix can be incorporated into geomechanical models.

  15. Use of TOUGHREACT to Simulate Effects of Fluid Chemistry onInjectivity in Fractured Geothermal Reservoirs with High Ionic StrengthFluids

    SciTech Connect

    Xu, Tianfu; Zhang, Guoxiang; Pruess, Karsten

    2005-02-09

    Recent studies suggest that mineral dissolution/precipitation and clay swelling effects could have a major impact on the performance of hot dry rock (HDR) and hot fractured rock (HFR) reservoirs. A major concern is achieving and maintaining adequate injectivity, while avoiding the development of preferential short-circuiting flow paths. A Pitzer ionic interaction model has been introduced into the publicly available TOUGHREACT code for solving non-isothermal multi-phase reactive geochemical transport problems under conditions of high ionic strength, expected in typical HDR and HFR systems. To explore chemically-induced effects of fluid circulation in these systems, we examine ways in which the chemical composition of reinjected waters can be modified to improve reservoir performance. We performed a number of coupled thermo-hydrologic-chemical simulations in which the fractured medium was represented by a one-dimensional MINC model (multiple interacting continua). Results obtained with the Pitzer activity coefficient model were compared with those using an extended Debye-Hueckel equation. Our simulations show that non-ideal activity effects can be significant even at modest ionic strength, and can have major impacts on permeability evolution in injection-production systems. Alteration of injection water chemistry, for example by dilution with fresh water, can greatly alter precipitation and dissolution effects, and can offer a powerful tool for operating hot dry rock and hot fractured rock reservoirs in a sustainable manner.

  16. Reservoir Modeling by Data Integration via Intermediate Spaces and Artificial Intelligence Tools in MPS Simulation Frameworks

    SciTech Connect

    Ahmadi, Rouhollah; Khamehchi, Ehsan

    2013-12-15

    Conditioning stochastic simulations are very important in many geostatistical applications that call for the introduction of nonlinear and multiple-point data in reservoir modeling. Here, a new methodology is proposed for the incorporation of different data types into multiple-point statistics (MPS) simulation frameworks. Unlike the previous techniques that call for an approximate forward model (filter) for integration of secondary data into geologically constructed models, the proposed approach develops an intermediate space where all the primary and secondary data are easily mapped onto. Definition of the intermediate space, as may be achieved via application of artificial intelligence tools like neural networks and fuzzy inference systems, eliminates the need for using filters as in previous techniques. The applicability of the proposed approach in conditioning MPS simulations to static and geologic data is verified by modeling a real example of discrete fracture networks using conventional well-log data. The training patterns are well reproduced in the realizations, while the model is also consistent with the map of secondary data.

  17. Effects of urban flood-detention reservoirs on peak discharges and flood discharges and flood frequencies, and simulation of flood-detention reservoir outflow hydrographs in two watersheds in Albany, Georgia

    USGS Publications Warehouse

    Hess, G.W.; Inman, E.J.

    1994-01-01

    This report describes the effects of flood-detention reservoirs on downstream peak discharges of two urban tributaries to Kinchafoonee Creek (tributaries 1 and 2) in Albany, Georgia and presents simulated flood-detention reservoir outflow hydrographs. Rainfall-runoff data were collected for six years at two stations in these two urban watersheds. Tributary number 1 basin has a drainage area of 0.12 square miles, contains 23.8 percent impervious area, and contains two detention reservoirs. Tributary number 2 basin has a drainage area of 0.09 square miles, contains 12.9 percent impervious area, and has one detention reservoir. The Distributed Routing Rainfall-Runoff Model (DR3M) was calibrated using rainfall-runoff data collected during 1987- 92 at each station. DR3M was then used to simulate long-term (1906-33, 1941-73) peak discharges for these stations for conditions ranging from the existing condition with all detention reservoirs in place to the condition of no detention reservoirs. Flood-frequency relations based on the long-term peak discharges were developed for each simulation by fitting the logarithms of the annual peak discharge data to a Pearson type III distri- bution curve. The effect of detention reservoirs on peak discharge data to a Pearson type III distributio curve. The effect of detention reservoirs on peak discharges was determined by comparison of simulated flood-frequency peak discharges for conditions with and without the detention reservoirs. The comparisons indicated that the removal of flood-detention reservoirs from the tributary number 1 basin would increase the 10-, 50-, and 100-year peak discharges by 164 to 204 percent. Removal of the reservoir from tributary number 2 basin would increase these discharges by about 145 percent.

  18. Simulating reservoir releases to mitigate climate impacts on fish sustainability below Shasta Lake using stochastic and mechanistic modeling approaches

    NASA Astrophysics Data System (ADS)

    Sapin, J. R.; Saito, L.; Rajagopalan, B.; Caldwell, R. J.

    2013-12-01

    Preservation of the Chinook salmon fishery on the Sacramento River in California has been a major concern since the winter-run Chinook was listed as threatened in 1989. The construction of Shasta Dam and Reservoir in 1945 prevented the salmon from reaching their native cold-water spawning habitat, resulting in severe population declines. The temperature control device (TCD) installed at Shasta Dam in 1997 provides increased capabilities of supplying cold-water habitat downstream of the dam to stimulate salmon spawning. However, increased air temperatures due to climate change could make it more difficult to meet downstream temperature targets with the TCD. By coupling stochastic hydroclimatology generation with two-dimensional hydrodynamic modeling of the reservoir we can simulate TCD operations under extreme climate conditions. This is accomplished by stochastically generating climate and inflow scenarios (created with historical data from NOAA, USGS and USBR) as input into a CE-QUAL-W2 model of the reservoir that can simulate TCD operations. Simulations will investigate if selective withdrawal from multiple gates of the TCD are capable of meeting temperature targets downstream of the dam under extreme hydroclimatic conditions. Moreover, our non-parametric methods for stochastically generating climate and inflow scenarios are capable of producing statistically representative years of extreme wet or extreme dry conditions beyond what is seen in the historical record. This allows us to simulate TCD operations for unprecedented hydroclimatic conditions with implications for climate changes in the watershed. Preliminary results of temperature outputs from simulations of TCD operations under extreme climate conditions with CE-QUAL-W2 will be presented. The conditions chosen for simulation are grounded to real-world managerial concerns by utilizing collaborative workshops with reservoir managers to establish which hydroclimatic scenarios would be of most concern for

  19. Simulating the gas hydrate production test at Mallik using the pilot scale pressure reservoir LARS

    NASA Astrophysics Data System (ADS)

    Heeschen, Katja; Spangenberg, Erik; Schicks, Judith M.; Priegnitz, Mike; Giese, Ronny; Luzi-Helbing, Manja

    2014-05-01

    LARS, the LArge Reservoir Simulator, allows for one of the few pilot scale simulations of gas hydrate formation and dissociation under controlled conditions with a high resolution sensor network to enable the detection of spatial variations. It was designed and built within the German project SUGAR (submarine gas hydrate reservoirs) for sediment samples with a diameter of 0.45 m and a length of 1.3 m. During the project, LARS already served for a number of experiments simulating the production of gas from hydrate-bearing sediments using thermal stimulation and/or depressurization. The latest test simulated the methane production test from gas hydrate-bearing sediments at the Mallik test site, Canada, in 2008 (Uddin et al., 2011). Thus, the starting conditions of 11.5 MPa and 11°C and environmental parameters were set to fit the Mallik test site. The experimental gas hydrate saturation of 90% of the total pore volume (70 l) was slightly higher than volumes found in gas hydrate-bearing formations in the field (70 - 80%). However, the resulting permeability of a few millidarcy was comparable. The depressurization driven gas production at Mallik was conducted in three steps at 7.0 MPa - 5.0 MPa - 4.2 MPa all of which were used in the laboratory experiments. In the lab the pressure was controlled using a back pressure regulator while the confining pressure was stable. All but one of the 12 temperature sensors showed a rapid decrease in temperature throughout the sediment sample, which accompanied the pressure changes as a result of gas hydrate dissociation. During step 1 and 2 they continued up to the point where gas hydrate stability was regained. The pressure decreases and gas hydrate dissociation led to highly variable two phase fluid flow throughout the duration of the simulated production test. The flow rates were measured continuously (gas) and discontinuously (liquid), respectively. Next to being discussed here, both rates were used to verify a model of gas

  20. Impact of integrated 3D reservoir modeling/flow simulation on development of deepwater sands, Mars Field, Gulf of Mexico

    SciTech Connect

    Lerch, C.S.; Johnston, J.R.; Juedeman, M.E.

    1996-12-31

    Prospect Mars is a major Gulf of Mexico deep water oil discovery made under joint partnership between operator Shell Offshore Inc. and partner British Petroleum Inc. The discovery lies in 3000 feet of water, located 130 miles southeast of New Orleans, Louisiana. The field was discovered in 1989 and to date 14 significant reservoir intervals from 10,000 to 19000 feet below sea level have been penetrated. Estimated recoverable reserves for the first phase of field development are in excess of 500 MMBE and development plans include installation of a 24 slot tension leg platform and two subsea wells, with first production in mid-1996. Over a two year period a comprehensive effort was directed at creating a new set of reservoir models utilizing an integrated software package developed at Shell E&P Technology Co. This package is able to incorporate pertinent geological, geophysical, and petrophysical data into 3-D reservoir models which can be used to: (1) estimate reservoir quantity, quality, and continuity, (2) predict reservoir production performance, (3) select development well locations, and (4) facilitate reserve estimation. This software allows interpretations from 3-D seismic, well control, and analog outcrops to be effectively integrated and passed to the reservoir model for flow simulation. This integrated effort at modeling ensured a more realistic reservoir picture upon which to base field development. Almost all the development wells pre-drilled prior to platform installation have been affected or designed based on these reservoir models and well results have been used to keep the models updated and evergreen.

  1. Impact of integrated 3D reservoir modeling/flow simulation on development of deepwater sands, Mars Field, Gulf of Mexico

    SciTech Connect

    Lerch, C.S.; Johnston, J.R.; Juedeman, M.E. )

    1996-01-01

    Prospect Mars is a major Gulf of Mexico deep water oil discovery made under joint partnership between operator Shell Offshore Inc. and partner British Petroleum Inc. The discovery lies in 3000 feet of water, located 130 miles southeast of New Orleans, Louisiana. The field was discovered in 1989 and to date 14 significant reservoir intervals from 10,000 to 19000 feet below sea level have been penetrated. Estimated recoverable reserves for the first phase of field development are in excess of 500 MMBE and development plans include installation of a 24 slot tension leg platform and two subsea wells, with first production in mid-1996. Over a two year period a comprehensive effort was directed at creating a new set of reservoir models utilizing an integrated software package developed at Shell E P Technology Co. This package is able to incorporate pertinent geological, geophysical, and petrophysical data into 3-D reservoir models which can be used to: (1) estimate reservoir quantity, quality, and continuity, (2) predict reservoir production performance, (3) select development well locations, and (4) facilitate reserve estimation. This software allows interpretations from 3-D seismic, well control, and analog outcrops to be effectively integrated and passed to the reservoir model for flow simulation. This integrated effort at modeling ensured a more realistic reservoir picture upon which to base field development. Almost all the development wells pre-drilled prior to platform installation have been affected or designed based on these reservoir models and well results have been used to keep the models updated and evergreen.

  2. Conformity assessment for seismic monitoring and reservoir simulation at the Ketzin pilot site - how much conformity can be reached?

    NASA Astrophysics Data System (ADS)

    Lüth, Stefan; Ivanova, Alexandra; Kempka, Thomas

    2016-04-01

    The EU CCS Directive defines three high-level criteria which have to be fulfilled by a site operator in the post closure phase of a storage site before liability can be transferred to the public after site closure. One of these high-level requirements is "Demonstrating conformity between observed and simulated plume behaviour". The observed plume behaviour is derived from geophysical and/or geochemical monitoring. Repeated 3D seismic observations have proven to provide the most comprehensive image of a CO2 plume in various projects such as Sleipner, Weyburn, or Ketzin. The simulated plume behaviour is derived from reservoir simulation using a model calibrated with monitoring results. Plume observations using any monitoring method are always affected by limited resolution and detection ability, and reservoir simulations will only be able to provide an approximated representation of the occurring reservoir processes. Therefore, full conformity between observed and simulated plume behaviour is difficult to achieve, if it is at all. It is therefore of crucial importance for each storage site to understand to what degree conformity can be achieved under realistic conditions, comprising noise affected monitoring data and reservoir models based on geological uncertainties. We applied performance criteria (plume footprint area, lateral migration distance, plume volume, and similarity index) for a comparison between monitoring results (4D seismic measurements) and reservoir simulations, considering a range of seismic amplitude values as noise threshold and a range of minimum thickness of the simulated CO2 plume. Relating the performance criteria to the noise and thickness threshold values allows assessing the quality of conformance between simulated and observed behaviour of a CO2 plume. The Ketzin site is provided with a comprehensive monitoring data set and a history-matched reservoir model. Considering the relatively high noise level, which is inherent for land

  3. Eos modeling and reservoir simulation study of bakken gas injection improved oil recovery in the elm coulee field, Montana

    NASA Astrophysics Data System (ADS)

    Pu, Wanli

    The Bakken Formation in the Williston Basin is one of the most productive liquid-rich unconventional plays. The Bakken Formation is divided into three members, and the Middle Bakken Member is the primary target for horizontal wellbore landing and hydraulic fracturing because of its better rock properties. Even with this new technology, the primary recovery factor is believed to be only around 10%. This study is to evaluate various gas injection EOR methods to try to improve on that low recovery factor of 10%. In this study, the Elm Coulee Oil Field in the Williston Basin was selected as the area of interest. Static reservoir models featuring the rock property heterogeneity of the Middle Bakken Member were built, and fluid property models were built based on Bakken reservoir fluid sample PVT data. By employing both compositional model simulation and Todd-Longstaff solvent model simulation methods, miscible gas injections were simulated and the simulations speculated that oil recovery increased by 10% to 20% of OOIP in 30 years. The compositional simulations yielded lower oil recovery compared to the solvent model simulations. Compared to the homogeneous model, the reservoir model featuring rock property heterogeneity in the vertical direction resulted in slightly better oil recovery, but with earlier CO2 break-through and larger CO2 production, suggesting that rock property heterogeneity is an important property for modeling because it has a big effect on the simulation results. Long hydraulic fractures shortened CO2 break-through time greatly and increased CO 2 production. Water-alternating-gas injection schemes and injection-alternating-shut-in schemes can provide more options for gas injection EOR projects, especially for gas production management. Compared to CO2 injection, separator gas injection yielded slightly better oil recovery, meaning separator gas could be a good candidate for gas injection EOR; lean gas generated the worst results. Reservoir

  4. Drought Prediction Technique For Reservoir Operation

    NASA Astrophysics Data System (ADS)

    Khayyat Kholghi, M.

    After three years critical drought period, Iran is facing a great problem on water re- sources. In this situation, several reservoir dams in this country can play an efficient role for solving an important part of this problem. One of the essential parameters in reservoir simulation is the inflow discharge. Predicting the likelihood of a drought markedly enhance the efficiency of reservoir operation. This study applies the kriging method to predict inflows to Dez reservoir dam in Southern Iran. Thirty-eight years of the one, ten, fifteen and monthly period reservoir inflow data were used in this study. A subsequent reservoir operation simulation is employed to determine the drought lead-time. An efficient reservoir operational strategy can be established with the aid of this time and the probability of successful drought prediction. According to simu- lation results, the kriging approach in terms of the number droughts being correctly predicted the drought lead-time and the probability of success of drought prediction. Keywords: Drought-prediction-kriging-reservoir operation-simulation-Dez reservoir dam-Iran

  5. Numerical simulation of the environmental impact of hydraulic fracturing of tight/shale gas reservoirs on near-surface groundwater: Background, base cases, shallow reservoirs, short-term gas, and water transport

    NASA Astrophysics Data System (ADS)

    Reagan, Matthew T.; Moridis, George J.; Keen, Noel D.; Johnson, Jeffrey N.

    2015-04-01

    Hydrocarbon production from unconventional resources and the use of reservoir stimulation techniques, such as hydraulic fracturing, has grown explosively over the last decade. However, concerns have arisen that reservoir stimulation creates significant environmental threats through the creation of permeable pathways connecting the stimulated reservoir with shallower freshwater aquifers, thus resulting in the contamination of potable groundwater by escaping hydrocarbons or other reservoir fluids. This study investigates, by numerical simulation, gas and water transport between a shallow tight-gas reservoir and a shallower overlying freshwater aquifer following hydraulic fracturing operations, if such a connecting pathway has been created. We focus on two general failure scenarios: (1) communication between the reservoir and aquifer via a connecting fracture or fault and (2) communication via a deteriorated, preexisting nearby well. We conclude that the key factors driving short-term transport of gas include high permeability for the connecting pathway and the overall volume of the connecting feature. Production from the reservoir is likely to mitigate release through reduction of available free gas and lowering of reservoir pressure, and not producing may increase the potential for release. We also find that hydrostatic tight-gas reservoirs are unlikely to act as a continuing source of migrating gas, as gas contained within the newly formed hydraulic fracture is the primary source for potential contamination. Such incidents of gas escape are likely to be limited in duration and scope for hydrostatic reservoirs. Reliable field and laboratory data must be acquired to constrain the factors and determine the likelihood of these outcomes.

  6. Numerical simulation of the environmental impact of hydraulic fracturing of tight/shale gas reservoirs on near-surface groundwater: Background, base cases, shallow reservoirs, short-term gas, and water transport

    SciTech Connect

    Reagan, Matthew T.; Moridis, George J.; Keen, Noel D.; Johnson, Jeffrey N.

    2015-04-18

    Hydrocarbon production from unconventional resources and the use of reservoir stimulation techniques, such as hydraulic fracturing, has grown explosively over the last decade. However, concerns have arisen that reservoir stimulation creates significant environmental threats through the creation of permeable pathways connecting the stimulated reservoir with shallower freshwater aquifers, thus resulting in the contamination of potable groundwater by escaping hydrocarbons or other reservoir fluids. This study investigates, by numerical simulation, gas and water transport between a shallow tight-gas reservoir and a shallower overlying freshwater aquifer following hydraulic fracturing operations, if such a connecting pathway has been created. We focus on two general failure scenarios: (1) communication between the reservoir and aquifer via a connecting fracture or fault and (2) communication via a deteriorated, preexisting nearby well. We conclude that the key factors driving short-term transport of gas include high permeability for the connecting pathway and the overall volume of the connecting feature. Production from the reservoir is likely to mitigate release through reduction of available free gas and lowering of reservoir pressure, and not producing may increase the potential for release. We also find that hydrostatic tight-gas reservoirs are unlikely to act as a continuing source of migrating gas, as gas contained within the newly formed hydraulic fracture is the primary source for potential contamination. Such incidents of gas escape are likely to be limited in duration and scope for hydrostatic reservoirs. Reliable field and laboratory data must be acquired to constrain the factors and determine the likelihood of these outcomes.

  7. Simulation-based inexact chance-constrained nonlinear programming for eutrophication management in the Xiangxi Bay of Three Gorges Reservoir.

    PubMed

    Huang, Y L; Huang, G H; Liu, D F; Zhu, H; Sun, W

    2012-10-15

    Although integrated simulation and optimization approaches under stochastic uncertainty have been applied to eutrophication management problems, few studies are reported in eutrophication control planning where multiple formats of uncertainties and nonlinearities are addressed in forms of intervals and probabilistic distributions within an integrated framework. Since the impounding of Three Gorges Reservoir (TGR), China in 2003, the hydraulic conditions and aquatic environment of the Xiangxi Bay (XXB) have changed significantly. The resulting emergence of eutrophication and algal blooms leads to its deteriorated water quality. The XXB becomes an ideal case study area. Thus, a simulation-based inexact chance-constrained nonlinear programming (SICNP) model is developed and applied to eutrophication control planning in the XXB of the TGR under uncertainties. In the SICNP, the wastewater treatment costs for removing total phosphorus (TP) are set as the objective function; effluent discharge standards, stream water quality standards and eutrophication control standards are considered in the constraints; a steady-state simulation model for phosphorus transport and fate is embedded in the environmental standards constraints; the interval programming and chance-constrained approaches are integrated to provide interval decision variables but also the associated risk levels in violating the system constraints. The model results indicate that changes in the violating level (q) will result in different strategy distributions at spatial and temporal scales; the optimal value of cost objective is from [2.74, 13.41] million RMB to [2.25, 13.08] million RMB when q equals from 0.01 to 0.25; the required TP treatment efficiency for the Baisha plant is the most stringent, which is followed by the Xiakou Town and the Zhaojun Town, while the requirement for the Pingyikou cement plant is the least stringent. The model results are useful for making optimal policies on eutrophication

  8. Simulated Last Glacial Maximum Δ14CATM and the Deep Glacial Ocean Reservoir

    NASA Astrophysics Data System (ADS)

    Mariotti, V.; Paillard, D.; Roche, D. M.; Bouttes, N.; Bopp, L.

    2012-12-01

    Δ14Catm has been estimated at 420 ± 80‰ (INTCAL09) during the Last Glacial Maximum (LGM) compared to preindustrial times (0‰), but mechanisms explaining this difference are not yet resolved. Δ14Catm is a function of cosmogenic production in high atmosphere and of carbon cycling in the Earth system (through carbon exchange with the superficial reservoirs, ocean and continental biosphere). 10Be-based reconstructions show a contribution of the cosmogenic production term of only 200 ± 200‰ at the LGM. The remaining 220‰ of Δ14Catm variation between the LGM and preindustrial times have thus to be explained by changes in the carbon cycle. Recently, Bouttes et al. (2010) proposed to explain most of the difference in atmospheric pCO2 between glacial and interglacial times by brine-induced ocean stratification in the Southern Ocean. This mechanism involves the formation of very saline water masses that can store Dissolved Inorganic Carbon (DIC) in the deep ocean. During glacial times, the sinking of brines is enhanced and more DIC is stored in the deep ocean, lowering atmospheric pCO2. Such an isolated ocean reservoir would be characterized by a low Δ14C signature. Evidence of such 14C-depleted deep waters during the LGM has recently been found in the Southern Ocean (Skinner et al., 2010). The degassing of this carbon with low Δ14C would then reduce Δ14Catm throughout the deglaciation. We have further developed the CLIMBER-2 model to include a cosmogenic production of 14C as well as an interactive atmospheric 14C reservoir. We investigate the role of both sinking of brines and cosmogenic production, alongside iron and vertical diffusion mechanisms to explain changes in Δ14Catm during the last deglaciation. In our simulations, not only the sinking of brine mechanism is consistent with past Δ14C data but also it explains most of the differences in atmospheric pCO2 and Δ14C between LGM and preindustrial times.

  9. Simulation and Analysis of Converging Shock Wave Test Problems

    SciTech Connect

    Ramsey, Scott D.; Shashkov, Mikhail J.

    2012-06-21

    Results and analysis pertaining to the simulation of the Guderley converging shock wave test problem (and associated code verification hydrodynamics test problems involving converging shock waves) in the LANL ASC radiation-hydrodynamics code xRAGE are presented. One-dimensional (1D) spherical and two-dimensional (2D) axi-symmetric geometric setups are utilized and evaluated in this study, as is an instantiation of the xRAGE adaptive mesh refinement capability. For the 2D simulations, a 'Surrogate Guderley' test problem is developed and used to obviate subtleties inherent to the true Guderley solution's initialization on a square grid, while still maintaining a high degree of fidelity to the original problem, and minimally straining the general credibility of associated analysis and conclusions.

  10. Snow-melt Runoff Simulation for Dam Reservoir in the Heavy Snow Region

    NASA Astrophysics Data System (ADS)

    Sato, Yoshinobu; Murai, Akihiro; Sumi, Tetsuya

    2014-05-01

    Increases in land surface temperature will have a significant affect on the hydrological cycle, particularly in regions where the available water resources are mainly dominated by the melting snow or ice. Thus, to clarify the impact of climate change on river discharge in cold and mountainous region is becoming one of the urgent issues for policy making and planning for the integrated river water management under the inevitable warming climate. However, in order to study climate change impacts on water resources in the heavy snow region, snow-melt runoff simulation for dam reservoir should be improved. Because, the available meteorological data for runoff simulation is quite limited, especially in a mountainous regions in Japan. In this study, we analyzed the inflow into the Okutadami Dam in the Agano River basin located in the northern mountainous region in Japan by a distributed hydrological model (Hydro-BEAM: Hydrological river Basin Environment Assessment Model). The Okutadami dam has an important role as one of the largest hydro-power generation dam in Japan. The result of our initial simulation underestimated the inflow significantly, especially in snow melting season, because of small input precipitation. We firstly modified the input precipitation by the JMA (Japan Meteorological Agency)'s climatic value 2010 (monthly 1km2 averaged mesh based precipitation dataset during the period from 1981 to 2010). Due to the modification, simulated annual mean river discharge (water balance) was improved significantly. Secondly, we modified the threshold temperature which distinguishes rainfall and snowfall improved the reproducibility slightly. Lastly, we modified the monthly discharge variation (seasonal change pattern) in snow melting period by considering the effect of heat supply by rainfall on snow surface layer. Consequently, we found that the calculated inflow to the Okutadami dam agreed well the observation. These methods will contribute to clarify the

  11. Numerical simulation of pore size dependent anhydrite precipitation in geothermal reservoirs

    NASA Astrophysics Data System (ADS)

    Mürmann, Mario; Kühn, Michael; Pape, Hansgeorg; Clauser, Christoph

    2013-04-01

    Porosity and permeability of reservoirs are key parameters for an economical use of hot water from geothermal installations and can be significantly reduced by precipitation of minerals, such as anhydrite. The borehole Allermöhe 1 near Hamburg (Germany) represents a failed attempt of geothermal heat mining due to anhydrite precipitation (Baermann et al. 2000). For a risk assessment of future boreholes it is essential to understand how and when anhydrite cementation occurred under reservoir conditions. From core samples of the Allermöhe borehole it was determined that anhydrite precipitation took place in regions of relatively high porosity while regions of low porosity remained uncemented (Wagner et al. 2005). These findings correspond to the fact that e.g. halite precipitation in porous media is found only in relatively large pores (Putnis and Mauthe 2001). This study and others underline that pore size controls crystallization and that it is therefore necessary to establish a relation between pore size and nucleation. The work presented here is based on investigations of Emmanuel and Berkowitz (2007) who present such a relation by applying a thermodynamic approach. However this approach cannot explain the heterogeneous precipitation observed in the Allermöhe core samples. We chose an advanced approach by considering electric system properties resulting in another relation between pore size and crystallization. It is well known that a high fluid supersaturation can be maintained in porous rocks (Putnis and Mauthe 2001). This clearly indicates that a supersaturation threshold exists exceeding thermodynamic equilibrium considerably. In order to quantify spatially heterogeneous anhydrite cementation a theoretical approach was chosen which considered the electric interaction between surface charges of the matrix and calcium and sulphate ions in the fluid. This approach was implemented into the numerical code SHEMAT (Clauser 2003) and used to simulate anhydrite

  12. Simulation of commercial scale CO2 injection into a fracture reservoir

    NASA Astrophysics Data System (ADS)

    Li, Y.; Li, S.; Zhang, Y.

    2011-12-01

    Geologic Carbon Sequestration is a proposed means to reduce atmospheric concentration of carbon dioxide (CO2). At Teapot Dome, Wyoming, CO2 will be injected into the Tensleep Formation, a depleted oil reservoir characterized with significant heterogeneity including facies, faults, and fractures. We've collected geological and engineering characterization data of the entire Teapot Dome field, including core data, well logs, seismic data, and production records. All data were screened for accuracy, before subsets of the data are used to build a geologic reservoir model. A formation structural model is created first by interpreting faults and stratigraphy from 3D seismic data. Formation MicroImager logs are analyzed for fracture characteristics to generate in-situ fracture intensity at wells, which is subsequently interpolated throughout the model with kriging. Based on the intensity and a set of geometric fracture parameters (constrained by outcrop and core measurements), a 3D stochastic Discrete Fracture Network (DFN) is created. One realization of the DFN is upscaled to a dual-porosity simulation model using a variant of the Oda's method. This method upscales the fracture network to an effective grid-block fracture permeability. To account for fluid transfer from matrix to fracture, a sigma factor is computed using average fracture spacings within the grid block. Matrix porosity is populated in the model by kriging interpolation of well-log-derived values. With the dual-porosity model, CO2 is injected near the crest of the dome (the proposed location) at a rate of 2.6 Mt/year for 50 years, with an injection bottomhole pressure set at 160% hydrostatic pressure. Boundary of the model is open except along one bounding fault assumed sealed. Results of the simulation suggest that provided that fluid pressure buildup is not an issue (simulated pressure buildup near the fault is minor), the Tensleep Formation at Teapot Dome can sustain commercial-scale injection over time

  13. Noble gas and hydrocarbon tracers in multiphase unconventional hydrocarbon systems: Toward integrated advanced reservoir simulators

    NASA Astrophysics Data System (ADS)

    Darrah, T.; Moortgat, J.; Poreda, R. J.; Muehlenbachs, K.; Whyte, C. J.

    2015-12-01

    Although hydrocarbon production from unconventional energy resources has increased dramatically in the last decade, total unconventional oil and gas recovery from black shales is still less than 25% and 9% of the totals in place, respectively. Further, the majority of increased hydrocarbon production results from increasing the lengths of laterals, the number of hydraulic fracturing stages, and the volume of consumptive water usage. These strategies all reduce the economic efficiency of hydrocarbon extraction. The poor recovery statistics result from an insufficient understanding of some of the key physical processes in complex, organic-rich, low porosity formations (e.g., phase behavior, fluid-rock interactions, and flow mechanisms at nano-scale confinement and the role of natural fractures and faults as conduits for flow). Noble gases and other hydrocarbon tracers are capably of recording subsurface fluid-rock interactions on a variety of geological scales (micro-, meso-, to macro-scale) and provide analogs for the movement of hydrocarbons in the subsurface. As such geochemical data enrich the input for the numerical modeling of multi-phase (e.g., oil, gas, and brine) fluid flow in highly heterogeneous, low permeability formations Herein we will present a combination of noble gas (He, Ne, Ar, Kr, and Xe abundances and isotope ratios) and molecular and isotopic hydrocarbon data from a geographically and geologically diverse set of unconventional hydrocarbon reservoirs in North America. Specifically, we will include data from the Marcellus, Utica, Barnett, Eagle Ford, formations and the Illinois basin. Our presentation will include geochemical and geological interpretation and our perspective on the first steps toward building an advanced reservoir simulator for tracer transport in multicomponent multiphase compositional flow (presented separately, in Moortgat et al., 2015).

  14. A comparative study of three simulation optimization algorithms for solving high dimensional multi-objective optimization problems in water resources

    NASA Astrophysics Data System (ADS)

    Schütze, Niels; Wöhling, Thomas; de Play, Michael

    2010-05-01

    Some real-world optimization problems in water resources have a high-dimensional space of decision variables and more than one objective function. In this work, we compare three general-purpose, multi-objective simulation optimization algorithms, namely NSGA-II, AMALGAM, and CMA-ES-MO when solving three real case Multi-objective Optimization Problems (MOPs): (i) a high-dimensional soil hydraulic parameter estimation problem; (ii) a multipurpose multi-reservoir operation problem; and (iii) a scheduling problem in deficit irrigation. We analyze the behaviour of the three algorithms on these test problems considering their formulations ranging from 40 up to 120 decision variables and 2 to 4 objectives. The computational effort required by each algorithm in order to reach the true Pareto front is also analyzed.

  15. Numerical simulations of the Macondo well blowout reveal strong control of oil flow by reservoir permeability and exsolution of gas.

    PubMed

    Oldenburg, Curtis M; Freifeld, Barry M; Pruess, Karsten; Pan, Lehua; Finsterle, Stefan; Moridis, George J

    2012-12-11

    In response to the urgent need for estimates of the oil and gas flow rate from the Macondo well MC252-1 blowout, we assembled a small team and carried out oil and gas flow simulations using the TOUGH2 codes over two weeks in mid-2010. The conceptual model included the oil reservoir and the well with a top boundary condition located at the bottom of the blowout preventer. We developed a fluid properties module (Eoil) applicable to a simple two-phase and two-component oil-gas system. The flow of oil and gas was simulated using T2Well, a coupled reservoir-wellbore flow model, along with iTOUGH2 for sensitivity analysis and uncertainty quantification. The most likely oil flow rate estimated from simulations based on the data available in early June 2010 was about 100,000 bbl/d (barrels per day) with a corresponding gas flow rate of 300 MMscf/d (million standard cubic feet per day) assuming the well was open to the reservoir over 30 m of thickness. A Monte Carlo analysis of reservoir and fluid properties provided an uncertainty distribution with a long tail extending down to 60,000 bbl/d of oil (170 MMscf/d of gas). The flow rate was most strongly sensitive to reservoir permeability. Conceptual model uncertainty was also significant, particularly with regard to the length of the well that was open to the reservoir. For fluid-entry interval length of 1.5 m, the oil flow rate was about 56,000 bbl/d. Sensitivity analyses showed that flow rate was not very sensitive to pressure-drop across the blowout preventer due to the interplay between gas exsolution and oil flow rate. PMID:21730177

  16. Numerical simulations of the Macondo well blowout reveal strong control of oil flow by reservoir permeability and exsolution of gas

    PubMed Central

    Oldenburg, Curtis M.; Freifeld, Barry M.; Pruess, Karsten; Pan, Lehua; Finsterle, Stefan; Moridis, George J.

    2012-01-01

    In response to the urgent need for estimates of the oil and gas flow rate from the Macondo well MC252-1 blowout, we assembled a small team and carried out oil and gas flow simulations using the TOUGH2 codes over two weeks in mid-2010. The conceptual model included the oil reservoir and the well with a top boundary condition located at the bottom of the blowout preventer. We developed a fluid properties module (Eoil) applicable to a simple two-phase and two-component oil-gas system. The flow of oil and gas was simulated using T2Well, a coupled reservoir-wellbore flow model, along with iTOUGH2 for sensitivity analysis and uncertainty quantification. The most likely oil flow rate estimated from simulations based on the data available in early June 2010 was about 100,000 bbl/d (barrels per day) with a corresponding gas flow rate of 300 MMscf/d (million standard cubic feet per day) assuming the well was open to the reservoir over 30 m of thickness. A Monte Carlo analysis of reservoir and fluid properties provided an uncertainty distribution with a long tail extending down to 60,000 bbl/d of oil (170 MMscf/d of gas). The flow rate was most strongly sensitive to reservoir permeability. Conceptual model uncertainty was also significant, particularly with regard to the length of the well that was open to the reservoir. For fluid-entry interval length of 1.5 m, the oil flow rate was about 56,000 bbl/d. Sensitivity analyses showed that flow rate was not very sensitive to pressure-drop across the blowout preventer due to the interplay between gas exsolution and oil flow rate. PMID:21730177

  17. Scale bridging in molecular simulation. Recurrent problems and current options

    NASA Astrophysics Data System (ADS)

    Hartmann, Carsten; Delle Site, Luigi

    2015-09-01

    Multiscale and multiphysics approaches have become an integral part of the molecular modeling and simulation toolbox and are used to attack various real-world problems that would be out of reach without these techniques. This special topics issue is devoted to a critical appraisal of some of the most popular scale bridging techniques for molecular simulation. It features regular articles and a "Discussion and Debate" section, in which experts in the field discuss specific articles and general aspects of scale bridging techniques.

  18. [Simulation of nitrogen and phosphorus loss in Siling Reservoir watershed with AnnAGNPS].

    PubMed

    Bian, Jin-yun; Wang, Fei-er; Yang, Jia; Yu, Jie; Lou, Li-ping; Yu, Dan-ping

    2012-08-01

    By using annual agricultural non-point source model (AnnAGNPS), this study simulated the export loading of nitrogen and phosphorus in Siling Reservoir watershed in Tiaoxi Basin, and integrated with the simulation results, the spatial distribution characteristics of non-point source pollution in the watershed was analyzed. The result showed that the export loading of nitrogen and phosphorus had similar characteristics: in the study area, the export loading of nutrients were higher in southern and western regions and lower in northern and eastern regions. Forest land mainly made up of bamboo was the main export source of nitrogen and phosphorus loading with the contribution above 90% of nutrient load of whole watershed. Three fertilization practices such as no fertilizer (CK), site-specific nutrient management (SSNM) and farmers' fertilizaction practice (FFP) were used in the scenario analysis. The scenario analysis showed that to a certain degree, SSNM could reduce the nitrogen and phosphorus loss. Comparing with FFP, the reduction of SSNM in dissolved nitrogen (DN), particle nitrogen (PN), dissolved phosphorus (DP) and particle phosphorus (PP) was 8.17%, 4.33%, 9.08% and 1.02%, respectively. PMID:23213887

  19. Numerical Simulations and Tracer Studies as a Tool to Support Water Circulation Modeling in Breeding Reservoirs

    NASA Astrophysics Data System (ADS)

    Zima, Piotr

    2014-12-01

    The article presents a proposal of a method for computer-aided design and analysis of breeding reservoirs in zoos and aquariums. The method applied involves the use of computer simulations of water circulation in breeding pools. A mathematical model of a pool was developed, and a tracer study was carried out. A simplified model of two-dimensional flow in the form of a biharmonic equation for the stream function (converted into components of the velocity vector) was adopted to describe the flow field. This equation, supplemented by appropriate boundary conditions, was solved numerically by the finite difference method. Next, a tracer migration equation was solved, which was a two-dimensional advection-dispersion equation describing the unsteady transport of a non-active, permanent solute. In order to obtain a proper solution, a tracer study (with rhodamine WT as a tracer) was conducted in situ. The results of these measurements were compared with numerical solutions obtained. The results of numerical simulations made it possible to reconstruct water circulation in the breading pool and to identify still water zones, where water circulation was impeded.

  20. Improving design factors of air diffuser systems based on field conditions of dam reservoirs: CFD simulation approach.

    PubMed

    Shin, Sangmin; Lee, Seungjae; Lee, Sangeun; Yum, Kyungtaek; Park, Heekyung

    2012-01-01

    This study aims to improve the design factors of air diffuser systems that have been analyzed in laboratory experiments, with consideration of the field conditions of dam reservoirs. In this study, the destratification number (D(N)), destratification radius, and efficiency are considered as design factors. The computational fluid dynamics (CFD) simulation experiment is performed in diverse field conditions in order to analyze these factors. The results illustrate the wider range of D(N) values in field conditions and the relationship of the destratification radius and efficiency to D(N). The results can lead to better performance of air diffuser systems and water quality management in dam reservoir sites. PMID:22678200

  1. Geothermal reservoir engineering code: comparison and validation

    SciTech Connect

    Not Available

    1981-02-27

    INTERCOMP has simulated six geothermal reservoir problems. INTERCOMP's geothermal reservoir model was used for all problems. No modifications were made to this model except to provide tabular output of the simulation results in the units used in RFP No. DE-RP03-80SF-10844. No difficulty was encountered in performing the problems described herein, although setting up the boundary and grid conditions exactly as specified were sometimes awkward, and minor modifications to the grid system were necessitated. The results of each problem are presented in tabular and (for many) graphical form.

  2. Development of a compositional model fully coupled with geomechanics and its application to tight oil reservoir simulation

    NASA Astrophysics Data System (ADS)

    Xiong, Yi

    Tight oil reservoirs have received great attention in recent years as unconventional and promising petroleum resources; they are reshaping the U.S. crude oil market due to their substantial production. However, fluid flow behaviors in tight oil reservoirs are not well studied or understood due to the complexities in the physics involved. Specific characteristics of tight oil reservoirs, such as nano-pore scale and strong stress-dependency result in complex porous medium fluid flow behaviors. Recent field observations and laboratory experiments indicate that large effects of pore confinement and rock compaction have non-negligible impacts on the production performance of tight oil reservoirs. On the other hand, there are approximations or limitations for modeling tight oil reservoirs under the effects of pore confinement and rock compaction with current reservoir simulation techniques. Thus this dissertation aims to develop a compositional model coupled with geomechanics with capabilities to model and understand the complex fluid flow behaviors of multiphase, multi-component fluids in tight oil reservoirs. MSFLOW_COM (Multiphase Subsurface FLOW COMpositional model) has been developed with the capability to model the effects of pore confinement and rock compaction for multiphase fluid flow in tight oil reservoirs. The pore confinement effect is represented by the effect of capillary pressure on vapor-liquid equilibrium (VLE), and modeled with the VLE calculation method in MSFLOW_COM. The fully coupled geomechanical model is developed from the linear elastic theory for a poro-elastic system and formulated in terms of the mean stress. Rock compaction is then described using stress-dependent rock properties, especially stress-dependent permeability. Thus MSFLOW_COM has the capabilities to model the complex fluid flow behaviors of tight oil reservoirs, fully coupled with geomechanics. In addition, MSFLOW_COM is validated against laboratory experimental data, analytical

  3. Geological and Petrophysical Characterization of the Ferron Sandstone for 3-D Simulation of a Fluvial-Deltaic Reservoir.

    SciTech Connect

    Allison, M.L.

    1997-07-01

    The objective of this project is to develop a comprehensive, interdisciplinary, and quantitative characterization of a fluvial- deltaic reservoir which will allow realistic inter-well and reservoir-scale modeling to be constructed for improved oil-field development in similar reservoirs world-wide. The geological and petrophysical properties of the Cretaceous Ferron Sandstone in east-central Utah will be quantitatively determined. Both new and existing data will be integrated into a three-dimensional representation of spatial variations in porosity, storativity, and tensorial rock permeability at a scale appropriate for inter-well to regional-scale reservoir simulation. Results could improve reservoir management through proper infill and extension drilling strategies, reduction of economic risks, increased recovery from existing oil fields, and more reliable reserve calculations. Transfer of the project results to the petroleum industry is an integral component of the project. Two activities continued this quarter as part of the geological and petrophysical characterization of the fluvial-deltaic Ferron Sandstone: (1) evaluation of the Ivie Creek case-study area and (2) technology transfer. The Ivie Creek case-study evaluation work during the quarter focused on the two parasequence sets, the Kf-1 and Kf-2, in the lower Ferron Sandstone. This work included: (1) clinoform characterization, (2) parasequence characterization from elevation and isopach maps, and (3) three-dimensional facies modeling. Scaled photomosaic panels from the Ivie Creek amphitheater (south-facing outcrop belt) and Quitchupah Canyon (Fig. 1) provide a deterministic framework for two apparent-dip cross sections. These panels along with other photomosaic coverage and data from five drill holes, ten stratigraphic sections, and 22 permeability transacts (Fig. 1), acquired during two field seasons, provided the necessary information for this geologic evaluation and creation of the models to be used

  4. Geological and Petrophysical Characterization of the Ferron Sandstone for 3-D Simulation of a Fluvial-Deltaic Reservoir

    SciTech Connect

    M. Lee Allison

    1997-03-01

    The objective of this project is to develop a comprehensive, interdisciplinary, and quantitative characterization of a fluvial-deltaic reser v oir which will allow realistic inter-well and reservoir-scale modeling to be constructed for improved oil-field development in similiar reservoirs world-wide. The geological and petrophysical properties of the Cretaceous Ferron Sandstone in east-central Utah will be quantitatively determined . Both new and existing data will be integrated into a three-dimensional representation of spatial variations in porosity, storativity, and tensorial rock permeability at a scale appropriate for inter-well to regional-scale reservoir simulation. Results could improve reservoir management through proper infill and extension drilling strategies, reduction of economic risks, increased recovery from existing oil fields, and more reliable reserve calculations . Transfer of the project results to the petroleum industry is an integral component of the project. Four activities continued this quarter as part of the geological and petrophysical characterization of the fluvial-deltaic Ferron Sandstone in the Ivie Creek case-study area: (1) geostatistics, (2) field description of clinoform bounding surfaces, (3) reservoir modeling, and (4) technology transfer.

  5. Application of advanced reservoir characterization, simulation, and production optimization strategies to maximize recovery in slope and basin clastic reservoirs, West Texas (Delaware Basin). Quarterly report, October 1 - December 31, 1996

    SciTech Connect

    Dutton, S.P.

    1997-01-01

    The objective of this project is to demonstrate that detailed reservoir characterization of slope and basin clastic reservoirs in sandstones of the Delaware Mountain Group in the Delaware Basin of West Texas and New Mexico is a cost effective way to recover a higher percentage of the original oil in place through strategic placement of infill wells and geologically based field development. Project objectives are divided into two major phases. The objectives of the reservoir characterization phase of the project are to provide a detailed understanding of the architecture and heterogeneity of two fields, the Ford Geraldine unit and Ford West field, which produce from the Bell Canyon and Cherry Canyon Formations, respectively, of the Delaware Mountain Group and to compare Bell Canyon and Cherry Canyon reservoirs. Reservoir characterization will utilize 3-D seismic data, high-resolution sequence stratigraphy, subsurface field studies, outcrop characterization, and other techniques. Once the reservoir-characterization study of both fields is completed, a pilot area of approximately 1 mi{sup 2} in one of the fields will be chosen for reservoir simulation. The objectives of the implementation phase of the project are to (1) apply the knowledge gained from reservoir characterization and simulation studies to increase recovery from the pilot area, (2) demonstrate that economically significant unrecovered oil remains in geologically resolvable untapped compartments, and (3) test the accuracy of reservoir characterization and flow simulation as predictive tools in resource preservation of mature fields. A geologically designed, enhanced-recovery program (CO{sub 2} flood, waterflood, or polymer flood) and well-completion program will be developed, and one to three infill wells will be drilled and cored. Technical progress is summarized for: geophysical characterization; reservoir characterization; outcrop characterization; and recovery technology identification and analysis.

  6. Application of advanced reservoir characterization, simulation, and production optimization strategies to maximize recovery in slope and basin clastic reservoirs, West Texas (Delaware Basin). Quarterly report, April 1,1996 - June 30, 1996

    SciTech Connect

    Dutton, S.P.

    1996-07-01

    The objective of this project is to demonstrate that detailed reservoir characterization of slope and basin clastic reservoirs in sandstones of the Delaware Mountain Group in the Delaware Basin of West Texas and New Mexico is a cost effective way to recover a higher percentage of the original oil in place through strategic placement of infill wells and geologically based field development. Project objectives are divided into two major phases. The objectives of the reservoir characterization phase of the project are to provide a detailed understanding of the architecture and heterogeneity of two fields, the Ford Geraldine unit and Ford West field, which produce from the Bell Canyon and Cherry Canyon Formations, respectively, of the Delaware Mountain Group and to compare Bell Canyon and Cherry Canyon reservoirs. Reservoir characterization will utilize 3-D seismic data, high-resolution sequence stratigraphy, subsurface field studies, outcrop characterization, and other techniques. Once the reservoir- characterization study of both fields is completed, a pilot area of approximately 1 mi{sup 2} in one of the fields will be chosen for reservoir simulation. The objectives of the implementation phase of the project are to (1) apply the knowledge gained from reservoir characterization and simulation studies to increase recovery from the pilot area, (2) demonstrate that economically significant unrecovered oil remains in geologically resolvable untapped compartments, and (3) test the accuracy of reservoir characterization and flow simulation as predictive tools in resource preservation of mature fields. A geologically designed, enhanced-recovery program (CO{sub 2} flood, waterflood, or polymer flood) and well-completion program will be developed, and one to three infill wells will be drilled and cored. Progress to date is summarized for reservoir characterization.

  7. Simulation of spatial and temporal distributions of non-point source pollution load in the Three Gorges Reservoir Region.

    PubMed

    Shen, Zhenyao; Qiu, Jiali; Hong, Qian; Chen, Lei

    2014-09-15

    Non-point source (NPS) pollution has become the largest threat to water quality in recent years. Major pollutants, particularly from agricultural activities, which include nitrogen, phosphorus and sediment that have been released into aquatic environments, have caused a range of problems in the Three Gorges Reservoir Region (TGRR), China. It is necessary to identify the spatial and temporal distributions of NPS pollutants and the highly polluted areas for the purpose of watershed management. In this study, the NPS pollutant load was simulated using the Soil and Water Assessment Tool (SWAT) and the small-scale watershed extended method (SWEM). The simulation results for four typical small catchments were extended to the entire watershed leading to estimates of the NPS load from 2001 to 2009. The results demonstrated that the NPS pollution load in the western area was the highest and that agricultural land was the primary pollutant source. The similar annual variation trends of runoff and sediment loads demonstrated that the sediment load was closely related to runoff. The loads of total nitrogen (TN) and total phosphorus (TP) were relatively stable from 2001 to 2007, except for high loads in 2006. The increase in pollution source strength was an important reason for the significant upward trend of TN and TP loads from 2008 to 2009. The rainfall from April to October contributed to the largest amount of runoff, sediment and nutrient loads for the year. The NPS load intensities in each sub-basin reveal large variations in the spatial distribution of different pollutants. It was shown that the temporal and spatial distributions of pollutant loads were positively correlated with the annual rainfall amounts and with human activities. Furthermore, this finding illustrates that conservation practices and nutrient management should be implemented in specific sites during special periods for the purpose of NPS pollution control in the TGRR. PMID:24946028

  8. Simulation of Sediment Transport Caused by Landslide at Nanhua Reservoir Watershed in Southern Taiwan

    NASA Astrophysics Data System (ADS)

    Lee, Ming-Hsi; Huang, Cong-Gi; Lin, Huan-Hsuan

    2016-04-01

    As a result of heavy rainfall, steep topography, young and weak geological formations, earthquakes, loose soils, slope land cultivation and other human disturbance, much area in Taiwan are prone to the occurrence of disastrous mass movements such as landslides and sediment disasters. During recent years, the extreme rainfall events brought huge amounts of rainfall and triggered severe changes in watershed environments. Typhoon Morakot in August 2009 caused severe landslides, debris flow, flooding and sediment disasters induced by record-break rainfall. The maximum rainfall of mountain area in Chiayi, Tainan, Kaohsiung and Pingtung County were over 2,900 mm. The study area is located at Nanhua reservoir watershed in southern Taiwan. The numerical model (HEC-RAS 4.1 and FLO-2D) will be used to simulate the sediment transport caused by landslide and the study will find out the separating location of erosion and deposition in the river, the danger area of riverbank, and the safety of the river terrace village under the return period of 50-year, 100-year and 200-year (such as Typhoon Morakot). The results of this study can provide for the disaster risk management of administrative decisions to lessen the impacts of natural hazards and may also be useful for time-space variation of sediment disasters caused by Climate Change.

  9. [Simulation as a method for solving medical technical problems].

    PubMed

    Wittmann, J; Obermayer, A

    1995-05-01

    The present paper first presents a comprehensive overview of the methodological aspects of simulation. It lists the various approaches to the modelling of real-world systems and emphasizes the respective constraints on the interpretation of model data. For this purpose, a classification of typical applications of simulation models is given, and is followed by a discussion of three structurally different ways of constructing a model. The correct matching of these modelling concepts to the particular application problem in hand is discussed, with special consideration being given to the drawing of conclusions from, and the interpretation of, the simulation results. To round of the classification, the simulation models are differentiated in accordance with the nature of their implementation. By defining application type, modelling concept and model implementation, an understanding of the validity and appropriate application of any industrial model is derived. PMID:7605886

  10. The Million-Body Problem: Particle Simulations in Astrophysics

    ScienceCinema

    Rasio, Fred [Northwestern University

    2009-09-01

    Computer simulations using particles play a key role in astrophysics. They are widely used to study problems across the entire range of astrophysical scales, from the dynamics of stars, gaseous nebulae, and galaxies, to the formation of the largest-scale structures in the universe. The 'particles' can be anything from elementary particles to macroscopic fluid elements, entire stars, or even entire galaxies. Using particle simulations as a common thread, this talk will present an overview of computational astrophysics research currently done in our theory group at Northwestern. Topics will include stellar collisions and the gravothermal catastrophe in dense star clusters.

  11. Numerical simulation of the environmental impact of hydraulic fracturing of tight/shale gas reservoirs on near-surface groundwater: Background, base cases, shallow reservoirs, short-term gas, and water transport

    PubMed Central

    Reagan, Matthew T; Moridis, George J; Keen, Noel D; Johnson, Jeffrey N

    2015-01-01

    Hydrocarbon production from unconventional resources and the use of reservoir stimulation techniques, such as hydraulic fracturing, has grown explosively over the last decade. However, concerns have arisen that reservoir stimulation creates significant environmental threats through the creation of permeable pathways connecting the stimulated reservoir with shallower freshwater aquifers, thus resulting in the contamination of potable groundwater by escaping hydrocarbons or other reservoir fluids. This study investigates, by numerical simulation, gas and water transport between a shallow tight-gas reservoir and a shallower overlying freshwater aquifer following hydraulic fracturing operations, if such a connecting pathway has been created. We focus on two general failure scenarios: (1) communication between the reservoir and aquifer via a connecting fracture or fault and (2) communication via a deteriorated, preexisting nearby well. We conclude that the key factors driving short-term transport of gas include high permeability for the connecting pathway and the overall volume of the connecting feature. Production from the reservoir is likely to mitigate release through reduction of available free gas and lowering of reservoir pressure, and not producing may increase the potential for release. We also find that hydrostatic tight-gas reservoirs are unlikely to act as a continuing source of migrating gas, as gas contained within the newly formed hydraulic fracture is the primary source for potential contamination. Such incidents of gas escape are likely to be limited in duration and scope for hydrostatic reservoirs. Reliable field and laboratory data must be acquired to constrain the factors and determine the likelihood of these outcomes. Key Points: Short-term leakage fractured reservoirs requires high-permeability pathways Production strategy affects the likelihood and magnitude of gas release Gas release is likely short-term, without additional driving forces PMID

  12. Application of advanced reservoir characterization, simulation, and production optimization strategies to maximize recovery in slope and basin clastic reservoirs, West Texas (Delaware Basin). Quarterly report, July 1 - September 30, 1996

    SciTech Connect

    Dutton, S.P.

    1996-10-01

    The objective of this project is to demonstrate that detailed reservoir characterization of slope and basin clastic reservoirs in sandstones of the Delaware Mountain Group in the Delaware Basin of West Texas and New Mexico is a cost effective way to recover a higher percentage of the original oil in place through strategic placement of infill wells and geologically based field development. Project objectives are divided into two major phases. The objectives of the reservoir characterization phase of the project are to provide a detailed understanding of the architecture and heterogeneity of two fields, the Ford Geraldine unit and Ford West field, which produce from the Bell Canyon and Cherry Canyon Formations, respectively, of the Delaware Mountain Group and to compare Bell Canyon and Cherry Canyon reservoirs. Reservoir characterization will utilize 3-D seismic data, high-resolution sequence stratigraphy, subsurface field studies, outcrop characterization, and other techniques. Once the reservoir- characterization study of both fields is completed, a pilot area of approximately 1 mi{sup 2} in one of the fields will be chosen for reservoir simulation. The objectives of the implementation phase of the project are to (1) apply the knowledge gained from reservoir characterization and simulation studies to increase recovery from the pilot area, (2) demonstrate that economically significant unrecovered oil remains in geologically resolvable untapped compartments, and (3) test the accuracy of reservoir characterization and flow simulation as predictive tools in resource preservation of mature fields. A geologically designed, enhanced-recovery program (CO{sup 2} flood, waterflood, or polymer flood) and well-completion program will be developed, and one to three infill wells will be drilled and cored. Accomplishments for this past quarter are discussed.

  13. Discrete Element Method Simulation of Nonlinear Viscoelastic Stress Wave Problems

    NASA Astrophysics Data System (ADS)

    Tang, Zhiping; Horie, Y.; Wang, Wenqiang

    2002-07-01

    A DEM(Discrete Element Method) simulation of nonlinear viscoelastic stress wave problems is carried out. The interaction forces among elements are described using a model in which neighbor elements are linked by a nonlinear spring and a certain number of Maxwell components in parallel. By making use of exponential relaxation moduli, it is shown that numerical computation of the convolution integral does not require storing and repeatedly calculating strain history, so that the computational cost is dramatically reduced. To validate the viscoelastic DM2 code1, stress wave propagation in a Maxwell rod with one end subjected to a constant stress loading is simulated. Results excellently fit those from the characteristics calculation. The code is then used to investigate the problem of meso-scale damage in a plastic-bonded explosive under shock loading. Results not only show "compression damage", but also reveal a complex damage evolution. They demonstrate a unique capability of DEM in modeling heterogeneous materials.

  14. Discrete Element Method Simulation of Nonlinear Viscoelastic Stress Wave Problems

    NASA Astrophysics Data System (ADS)

    Wang, Wenqiang; Tang, Zhiping; Horie, Y.

    2002-07-01

    A DEM(Discrete Element Method) simulation of nonlinear viscoelastic stress wave problems is carried out. The interaction forces among elements are described using a model in which neighbor elements are linked by a nonlinear spring and a certain number of Maxwell components in parallel. By making use of exponential relaxation moduli, it is shown that numerical computation of the convolution integral does not require storing and repeatedly calculating strain history, so that the computational cost is dramatically reduced. To validate the viscoelastic DM2 code[1], stress wave propagation in a Maxwell rod with one end subjected to a constant stress loading is simulated. Results excellently fit those from the characteristics calculation. The code is then used to investigate the problem of meso-scale damage in a plastic-bonded explosive under shock loading. Results not only show "compression damage", but also reveal a complex damage evolution. They demonstrate a unique capability of DEM in modeling heterogeneous materials.

  15. Bounded multi-scale plasma simulation: Application to sheath problems

    SciTech Connect

    Parker, S.E. ); Friedman, A.; Ray. S.L. ); Birdsall, C.K. )

    1993-08-01

    In our previous paper we introduced the multi-scale method, a self-consistent plasma simulation technique that allowed particles to have independent timesteps. Here we apply the method to one-dimensional electrostatic bounded plasma problems and demonstrate a significant reduction in computing time. We describe a technique to allow for variable grid spacing and develop consistent boundary conditions for the direct implicit method. Also discussed are criteria for specifying timestep size as a function of position in phase space. Next, an analytically solvable sheath problem is presented, and a comparison to simulation results in made. Finally, we show results for an ion acoustic shock front propagating toward a conducting wall. 20 refs., 16 figs., 2 tabs.

  16. Simulated monthly hydrologic data and estimated flood characteristics for Cherry Creek at a proposed reservoir site near Terry, Montana

    USGS Publications Warehouse

    Parrett, Charles; Johnson, D.R.

    1995-01-01

    A monthly hydrologic budget for water years 1937- 92 was developed for the proposed Cherry Creek Reservoir (maximum volume about 14,100 acre-feet). Monthly suspended-sediment loads and dissolved- solids concentrations in the reservoir and flood hydrographs and volumes having recurrence intervals of 25-, 50-, and 100-years were estimated. Monthly streamflow and precipitation were estimated using a mixed-station, record-extension procedure. Monthly suspended-sediment and dissolved-solids concentrations in the reservoir were estimated from regression relations between logarithms of concen- tration and streamflow. The simulation showed that flows that Cherry Creek generally were adequate to maintain the reservoir elevation above the minimum operating level for a seepage loss of 0 cubic ft per square. With a seepage loss rate of 3 cubic ft per square, diversions from the Yellowstone River were required for about on third of the months. Cumulative sediment deposition during the 56-year simulation period was about 138 acre-ft from Cherry Creek alone and was about 149 acre-ft when additional water was imported from the Yellowstone River. The concentration of dissolved solids in the reservoir reached a maximum value of about 2,540 mg/L for a seepage loss of 0 cubic ft per square. For a seepage loss of 3 cubic ft per square, water was imported from the Yellowstone River and the maximum concentration of dissolved solids was about 1,200 mg/L. Volumes for flood discharges were estimated from synthetic 24-hour duration storms that were used in a rainfall-runoff model (HEC-1).

  17. A reduced-order based CE-QUAL-W2 model for simulation of nitrate concentration in dam reservoirs

    NASA Astrophysics Data System (ADS)

    Noori, Roohollah; Yeh, Hund-Der; Ashrafi, Khosro; Rezazadeh, Najmeh; Bateni, Sayed M.; Karbassi, Abdulreza; Kachoosangi, Fatemeh Torabi; Moazami, Saber

    2015-11-01

    When the number of computational grids increases, water quality simulation complexities arise. Therefore, using a reduced order framework to express the variations of the objective parameter may facilitate the simulation task and also the interpretation of computational results. In this regard, a new reduced-order approach was proposed to link a water quality simulator model (CE-QUAL-W2, W2) with a data reduction technique (proper orthogonal decomposition, POD). The W2 model simulated spatio-temporal variations of nitrate in the Karkheh Reservoir, Iran. Thereafter, the POD model reduced the dimensions of simulated nitrate in the computational grids. The performance of the developed reduced-order model (ROM) results was evaluated through the comparison of the regenerated nitrate data by the model, and the simulated ones by W2. Findings indicated that the first four modes among 1825 computed ones by ROM conserved approximately over 91% of the nitrate variations. It means that the ROM was capable of showing the spatio-temporal variations of nitrate in the reservoir using the first few modes. Finally, confirmation of ROM indicated that the error of order of magnitude was less than 0.001, for nitrate ROM to regenerate nitrate data using 100 basis functions.

  18. Physical properties of rocks and aqueous fluids at conditions simulating near- and supercritical reservoirs

    NASA Astrophysics Data System (ADS)

    Kummerow, Juliane; Raab, Siegfried

    2016-04-01

    The growing interest in exploiting supercritical geothermal reservoirs calls for a thorough identification and understanding of physico-chemical processes occuring in geological settings with a high heat flow. In reservoir engineering, electrical sounding methods are common geophysical exploration and monitoring tools. However, a realistic interpretation of field measurements is based on the knowledge of both, the physical properties of the rock and those of the interacting fluid at defined temperature and pressure conditions. Thus, laboratory studies at simulated in-situ conditions provide a link between the field data and the material properties in the depth. The physico-chemical properties of fluids change dramatically above the critical point, which is for pure water 374.21 °C and 221.2 bar. In supercritical fluids mass transfer and diffusion-controlled chemical reactions are enhanced and cause mineral alterations. Also, ion mobility and ion concentration are affected by the change of physical state. All this cause changes in the electrical resistivity of supercritical fluids and may have considerable effects on the porosity and hydraulic properties of the rocks they are in contact with. While there are some datasets available for physical and chemical properties of water and single component salt solutions above their critical points, there exist nearly no data for electrical properties of mixed brines, representing the composition of natural geothermal fluids. Also, the impact of fluid-rock interactions on the electrical properties of multicomponent fluids in a supercritical region is scarcely investigated. For a better understanding of fluid-driven processes in a near- and supercritical geological environment, in the framework of the EU-funded FP7 program IMAGE we have measured (1) the electrical resistivity of geothermal fluids and (2) physical properties of fluid saturated rock samples at simulated in-situ conditions. The permeability and electrical

  19. A Condensed Disaggregation Model for Incorporating Parameter Uncertainty Into Monthly Reservoir Simulations

    NASA Astrophysics Data System (ADS)

    Stedinger, Jery R.; Pei, Daniel; Cohn, Timothy A.

    1985-05-01

    A condensed version of the Valencia-Schaake disaggregation model is developed which describes the distribution of monthly streamflow sequences using a set of coupled univariate regression models rather than a multivariate time series formulation. The condensed model has fewer parameters and is convenient for generating flow sequences which incorporate the intrinsic variability of streamflows and the uncertainty in the parameters of the annual and monthly streamflow models. The impact of parameter uncertainty on derived relationships between reservoir capacity and reservoir performance statistics is illustrated using required reservoir capacity (calculated with the sequent peak algorithm), system reliability, and the average total shortfall. Modeled sequences describe flows in the Rappahannock River in Virginia and the Boise River in Idaho. For high-reliability systems the results show that streamflow generation procedures which ignore model parameter uncertainty can grossly underestimate reservoir system failure rates and the severity of likely shortages, even if based on a 50-year record.

  20. Estimation of permeability of a sandstone reservoir by a fractal and Monte Carlo simulation approach: a case study

    NASA Astrophysics Data System (ADS)

    Vadapalli, U.; Srivastava, R. P.; Vedanti, N.; Dimri, V. P.

    2014-01-01

    Permeability of a hydrocarbon reservoir is usually estimated from core samples in the laboratory or from well test data provided by the industry. However, such data is very sparse and as such it takes longer to generate that. Thus, estimation of permeability directly from available porosity logs could be an alternative and far easier approach. In this paper, a method of permeability estimation is proposed for a sandstone reservoir, which considers fractal behavior of pore size distribution and tortuosity of capillary pathways to perform Monte Carlo simulations. In this method, we consider a reservoir to be a mono-dispersed medium to avoid effects of micro-porosity. The method is applied to porosity logs obtained from Ankleshwar oil field, situated in the Cambay basin, India, to calculate permeability distribution in a well. Computed permeability values are in good agreement with the observed permeability obtained from well test data. We also studied variation of permeability with different parameters such as tortuosity fractal dimension (Dt), grain size (r) and minimum particle size (d0), and found that permeability is highly dependent upon the grain size. This method will be extremely useful for permeability estimation, if the average grain size of the reservoir rock is known.

  1. Numerical simulation of the environmental impact of hydraulic fracturing of tight/shale gas reservoirs on near-surface groundwater: Background, base cases, shallow reservoirs, short-term gas, and water transport

    DOE PAGESBeta

    Reagan, Matthew T.; Moridis, George J.; Keen, Noel D.; Johnson, Jeffrey N.

    2015-04-18

    Hydrocarbon production from unconventional resources and the use of reservoir stimulation techniques, such as hydraulic fracturing, has grown explosively over the last decade. However, concerns have arisen that reservoir stimulation creates significant environmental threats through the creation of permeable pathways connecting the stimulated reservoir with shallower freshwater aquifers, thus resulting in the contamination of potable groundwater by escaping hydrocarbons or other reservoir fluids. This study investigates, by numerical simulation, gas and water transport between a shallow tight-gas reservoir and a shallower overlying freshwater aquifer following hydraulic fracturing operations, if such a connecting pathway has been created. We focus on twomore » general failure scenarios: (1) communication between the reservoir and aquifer via a connecting fracture or fault and (2) communication via a deteriorated, preexisting nearby well. We conclude that the key factors driving short-term transport of gas include high permeability for the connecting pathway and the overall volume of the connecting feature. Production from the reservoir is likely to mitigate release through reduction of available free gas and lowering of reservoir pressure, and not producing may increase the potential for release. We also find that hydrostatic tight-gas reservoirs are unlikely to act as a continuing source of migrating gas, as gas contained within the newly formed hydraulic fracture is the primary source for potential contamination. Such incidents of gas escape are likely to be limited in duration and scope for hydrostatic reservoirs. Reliable field and laboratory data must be acquired to constrain the factors and determine the likelihood of these outcomes.« less

  2. List-Based Simulated Annealing Algorithm for Traveling Salesman Problem

    PubMed Central

    Zhan, Shi-hua; Lin, Juan; Zhang, Ze-jun

    2016-01-01

    Simulated annealing (SA) algorithm is a popular intelligent optimization algorithm which has been successfully applied in many fields. Parameters' setting is a key factor for its performance, but it is also a tedious work. To simplify parameters setting, we present a list-based simulated annealing (LBSA) algorithm to solve traveling salesman problem (TSP). LBSA algorithm uses a novel list-based cooling schedule to control the decrease of temperature. Specifically, a list of temperatures is created first, and then the maximum temperature in list is used by Metropolis acceptance criterion to decide whether to accept a candidate solution. The temperature list is adapted iteratively according to the topology of the solution space of the problem. The effectiveness and the parameter sensitivity of the list-based cooling schedule are illustrated through benchmark TSP problems. The LBSA algorithm, whose performance is robust on a wide range of parameter values, shows competitive performance compared with some other state-of-the-art algorithms. PMID:27034650

  3. List-Based Simulated Annealing Algorithm for Traveling Salesman Problem.

    PubMed

    Zhan, Shi-hua; Lin, Juan; Zhang, Ze-jun; Zhong, Yi-wen

    2016-01-01

    Simulated annealing (SA) algorithm is a popular intelligent optimization algorithm which has been successfully applied in many fields. Parameters' setting is a key factor for its performance, but it is also a tedious work. To simplify parameters setting, we present a list-based simulated annealing (LBSA) algorithm to solve traveling salesman problem (TSP). LBSA algorithm uses a novel list-based cooling schedule to control the decrease of temperature. Specifically, a list of temperatures is created first, and then the maximum temperature in list is used by Metropolis acceptance criterion to decide whether to accept a candidate solution. The temperature list is adapted iteratively according to the topology of the solution space of the problem. The effectiveness and the parameter sensitivity of the list-based cooling schedule are illustrated through benchmark TSP problems. The LBSA algorithm, whose performance is robust on a wide range of parameter values, shows competitive performance compared with some other state-of-the-art algorithms. PMID:27034650

  4. CO2 Sequestration in Coalbed Methane Reservoirs: Experimental Studies and Computer Simulations

    SciTech Connect

    Muhammad Sahimi; Theodore T. Tsotsis

    2002-12-15

    One of the approaches suggested for sequestering CO{sub 2} is by injecting it in coalbed methane (CBM) reservoirs. Despite its potential importance for CO{sub 2} sequestration, to our knowledge, CO{sub 2} injection in CBM reservoirs for the purpose of sequestration has not been widely studied. Furthermore, a key element missing in most of the existing studies is the comprehensive characterization of the CBM reservoir structure. CBM reservoirs are complex porous media, since in addition to their primary pore structure, generated during coal formation, they also contain a variety of fractures, which may potentially play a key role in CO{sub 2} sequestration, as they generally provide high permeability flow paths for both CO{sub 2} and CH{sub 4}. In this report we present an overview of our ongoing experimental and modeling efforts, which aim to investigate the injection, adsorption and sequestration of CO{sub 2} in CBM reservoirs, the enhanced CH{sub 4} production that results, as well as the main factors that affect the overall operation. We describe the various experimental techniques that we utilize, and discuss their range of application and the value of the data generated. We conclude with a brief overview of our modeling efforts aiming to close the knowledge gap and fill the need in this area.

  5. Simulated annealing with probabilistic analysis for solving traveling salesman problems

    NASA Astrophysics Data System (ADS)

    Hong, Pei-Yee; Lim, Yai-Fung; Ramli, Razamin; Khalid, Ruzelan

    2013-09-01

    Simulated Annealing (SA) is a widely used meta-heuristic that was inspired from the annealing process of recrystallization of metals. Therefore, the efficiency of SA is highly affected by the annealing schedule. As a result, in this paper, we presented an empirical work to provide a comparable annealing schedule to solve symmetric traveling salesman problems (TSP). Randomized complete block design is also used in this study. The results show that different parameters do affect the efficiency of SA and thus, we propose the best found annealing schedule based on the Post Hoc test. SA was tested on seven selected benchmarked problems of symmetric TSP with the proposed annealing schedule. The performance of SA was evaluated empirically alongside with benchmark solutions and simple analysis to validate the quality of solutions. Computational results show that the proposed annealing schedule provides a good quality of solution.

  6. Solving the dynamic berth allocation problem by simulated annealing

    NASA Astrophysics Data System (ADS)

    Lin, Shih-Wei; Ting, Ching-Jung

    2014-03-01

    Berth allocation, the first operation when vessels arrive at a port, is one of the major container port optimization problems. From both the port operator's and the ocean carriers' perspective, minimization of the time a ship spends at berth is a key goal of port operations. This article focuses on two versions of the dynamic berth allocation problem (DBAP): discrete and continuous cases. The first case assigns ships to a given set of berth positions; the second one permits them to be moored anywhere along the berth. Simulated annealing (SA) approaches are proposed to solve the DBAP. The proposed SAs are tested with numerical instances for different sizes from the literature. Experimental results show that the proposed SA can obtain the optimal solutions in all instances of discrete cases and update 27 best known solutions in the continuous case.

  7. Numerical simulation of the electrical properties of shale gas reservoir rock based on digital core

    NASA Astrophysics Data System (ADS)

    Nie, Xin; Zou, Changchun; Li, Zhenhua; Meng, Xiaohong; Qi, Xinghua

    2016-08-01

    In this paper we study the electrical properties of shale gas reservoir rock by applying the finite element method to digital cores which are built based on an advanced Markov Chain Monte Carlo method and a combination workflow. Study shows that the shale gas reservoir rock has strong anisotropic electrical conductivity because the conductivity is significantly different in both horizontal and vertical directions. The Archie formula is not suitable for application in shale reservoirs. The formation resistivity decreases in two cases; namely (a) with the increase of clay mineral content and the cation exchange capacity of clay, and (b) with the increase of pyrite content. The formation resistivity is not sensitive to the solid organic matter but to the clay and gas in the pores.

  8. Dynamic simulation solves process control problem in Oman

    SciTech Connect

    1998-11-16

    A dynamic simulation study solved the process control problems for a Saih Rawl, Oman, gas compressor station operated by Petroleum Development of Oman (PDO). PDO encountered persistent compressor failure that caused frequent facility shutdowns, oil production deferment, and gas flaring. It commissioned MSE (Consultants) Ltd., U.K., to find a solution for the problem. Saih Rawl, about 40 km from Qarn Alam, produces oil and associated gas from a large number of low and high-pressure wells. Oil and gas are separated in three separators. The oil is pumped to Qarn Alam for treatment and export. Associated gas is compressed in two parallel trains. Train K-1115 is a 350,000 standard cu m/day, four-stage reciprocating compressor driven by a fixed-speed electric motor. Train K-1120 is a 1 million standard cu m/day, four-stage reciprocating compressor driven by a fixed-speed electric motor. Train K-1120 is a 1 million standard cu m/day, four-stage centrifugal compressor driven by a variable-speed motor. The paper describes tripping and surging problems with the gas compressor and the control simplifications that solved the problem.

  9. Simulation studies to evaluate the effect of fracture closure on the performance of naturally fractured reservoirs

    SciTech Connect

    Dauben, D.L.

    1991-07-15

    The study has two principal objectives: (1) To evaluate the effects of fracture closure on the recovery of oil and gas reserves from naturally fractured petroleum or natural gas reservoirs. (2) To evaluate procedures for improving the recovery of these reserves using innovative fluid injection techniques to maintain reservoir pressure and mitigate the impact of fracture closure. The total scope of the study has been subdivided into three main tasks: (1) Baseline studies (non-pressure sensitive fractures); (2)studies with pressure sensitive fractures; and (3) innovative approaches for improving oil recovery.

  10. Surrogate Reservoir Model

    NASA Astrophysics Data System (ADS)

    Mohaghegh, Shahab

    2010-05-01

    Surrogate Reservoir Model (SRM) is new solution for fast track, comprehensive reservoir analysis (solving both direct and inverse problems) using existing reservoir simulation models. SRM is defined as a replica of the full field reservoir simulation model that runs and provides accurate results in real-time (one simulation run takes only a fraction of a second). SRM mimics the capabilities of a full field model with high accuracy. Reservoir simulation is the industry standard for reservoir management. It is used in all phases of field development in the oil and gas industry. The routine of simulation studies calls for integration of static and dynamic measurements into the reservoir model. Full field reservoir simulation models have become the major source of information for analysis, prediction and decision making. Large prolific fields usually go through several versions (updates) of their model. Each new version usually is a major improvement over the previous version. The updated model includes the latest available information incorporated along with adjustments that usually are the result of single-well or multi-well history matching. As the number of reservoir layers (thickness of the formations) increases, the number of cells representing the model approaches several millions. As the reservoir models grow in size, so does the time that is required for each run. Schemes such as grid computing and parallel processing helps to a certain degree but do not provide the required speed for tasks such as: field development strategies using comprehensive reservoir analysis, solving the inverse problem for injection/production optimization, quantifying uncertainties associated with the geological model and real-time optimization and decision making. These types of analyses require hundreds or thousands of runs. Furthermore, with the new push for smart fields in the oil/gas industry that is a natural growth of smart completion and smart wells, the need for real time

  11. Variability of wet troposphere delays over inland reservoirs as simulated by a high-resolution regional climate model

    NASA Astrophysics Data System (ADS)

    Clark, E.; Lettenmaier, D. P.

    2014-12-01

    Satellite radar altimetry is widely used for measuring global sea level variations and, increasingly, water height variations of inland water bodies. Existing satellite radar altimeters measure water surfaces directly below the spacecraft (approximately at nadir). Over the ocean, most of these satellites use radiometry to measure the delay of radar signals caused by water vapor in the atmosphere (also known as the wet troposphere delay (WTD)). However, radiometry can only be used to estimate this delay over the largest inland water bodies, such as the Great Lakes, due to spatial resolution issues. As a result, atmospheric models are typically used to simulate and correct for the WTD at the time of observations. The resolutions of these models are quite coarse, at best about 5000 km2 at 30˚N. The upcoming NASA- and CNES-led Surface Water and Ocean Topography (SWOT) mission, on the other hand, will use interferometric synthetic aperture radar (InSAR) techniques to measure a 120-km-wide swath of the Earth's surface. SWOT is expected to make useful measurements of water surface elevation and extent (and storage change) for inland water bodies at spatial scales as small as 250 m, which is much smaller than current altimetry targets and several orders of magnitude smaller than the models used for wet troposphere corrections. Here, we calculate WTD from very high-resolution (4/3-km to 4-km) simulations of the Weather Research and Forecasting (WRF) regional climate model, and use the results to evaluate spatial variations in WTD. We focus on six U.S. reservoirs: Lake Elwell (MT), Lake Pend Oreille (ID), Upper Klamath Lake (OR), Elephant Butte (NM), Ray Hubbard (TX), and Sam Rayburn (TX). The reservoirs vary in climate, shape, use, and size. Because evaporation from open water impacts local water vapor content, we compare time series of WTD over land and water in the vicinity of each reservoir. To account for resolution effects, we examine the difference in WRF-simulated

  12. Coupling Hydraulic Fracturing Propagation and Gas Well Performance for Simulation of Production in Unconventional Shale Gas Reservoirs

    NASA Astrophysics Data System (ADS)

    Wang, C.; Winterfeld, P. H.; Wu, Y. S.; Wang, Y.; Chen, D.; Yin, C.; Pan, Z.

    2014-12-01

    Hydraulic fracturing combined with horizontal drilling has made it possible to economically produce natural gas from unconventional shale gas reservoirs. An efficient methodology for evaluating hydraulic fracturing operation parameters, such as fluid and proppant properties, injection rates, and wellhead pressure, is essential for the evaluation and efficient design of these processes. Traditional numerical evaluation and optimization approaches are usually based on simulated fracture properties such as the fracture area. In our opinion, a methodology based on simulated production data is better, because production is the goal of hydraulic fracturing and we can calibrate this approach with production data that is already known. This numerical methodology requires a fully-coupled hydraulic fracture propagation and multi-phase flow model. In this paper, we present a general fully-coupled numerical framework to simulate hydraulic fracturing and post-fracture gas well performance. This three-dimensional, multi-phase simulator focuses on: (1) fracture width increase and fracture propagation that occurs as slurry is injected into the fracture, (2) erosion caused by fracture fluids and leakoff, (3) proppant subsidence and flowback, and (4) multi-phase fluid flow through various-scaled anisotropic natural and man-made fractures. Mathematical and numerical details on how to fully couple the fracture propagation and fluid flow parts are discussed. Hydraulic fracturing and production operation parameters, and properties of the reservoir, fluids, and proppants, are taken into account. The well may be horizontal, vertical, or deviated, as well as open-hole or cemented. The simulator is verified based on benchmarks from the literature and we show its application by simulating fracture network (hydraulic and natural fractures) propagation and production data history matching of a field in China. We also conduct a series of real-data modeling studies with different combinations of

  13. Sparse approximation problem: how rapid simulated annealing succeeds and fails

    NASA Astrophysics Data System (ADS)

    Obuchi, Tomoyuki; Kabashima, Yoshiyuki

    2016-03-01

    Information processing techniques based on sparseness have been actively studied in several disciplines. Among them, a mathematical framework to approximately express a given dataset by a combination of a small number of basis vectors of an overcomplete basis is termed the sparse approximation. In this paper, we apply simulated annealing, a metaheuristic algorithm for general optimization problems, to sparse approximation in the situation where the given data have a planted sparse representation and noise is present. The result in the noiseless case shows that our simulated annealing works well in a reasonable parameter region: the planted solution is found fairly rapidly. This is true even in the case where a common relaxation of the sparse approximation problem, the G-relaxation, is ineffective. On the other hand, when the dimensionality of the data is close to the number of non-zero components, another metastable state emerges, and our algorithm fails to find the planted solution. This phenomenon is associated with a first-order phase transition. In the case of very strong noise, it is no longer meaningful to search for the planted solution. In this situation, our algorithm determines a solution with close-to-minimum distortion fairly quickly.

  14. Preliminary Three-Dimensional Simulation of Sediment and Cesium Transport in the Ogi Dam Reservoir using FLESCOT – Task 6, Subtask 2

    SciTech Connect

    Onishi, Yasuo; Kurikami, Hiroshi; Yokuda, Satoru T.

    2014-03-28

    After the accident at the Fukushima Daiichi Nuclear Power Plant in March 2011, the Japan Atomic Energy Agency and the Pacific Northwest National Laboratory initiated a collaborative project on environmental restoration. In October 2013, the collaborative team started a task of three-dimensional modeling of sediment and cesium transport in the Fukushima environment using the FLESCOT (Flow, Energy, Salinity, Sediment Contaminant Transport) code. As the first trial, we applied it to the Ogi Dam Reservoir that is one of the reservoirs in the Japan Atomic Energy Agency’s (JAEA’s) investigation project. Three simulation cases under the following different temperature conditions were studied: • incoming rivers and the Ogi Dam Reservoir have the same water temperature • incoming rivers have lower water temperature than that of the reservoir • incoming rivers have higher water temperature than that of the reservoir. The preliminary simulations suggest that seasonal temperature changes influence the sediment and cesium transport. The preliminary results showed the following: • Suspended sand, and cesium adsorbed by sand, coming into the reservoirs from upstream rivers is deposited near the reservoir entrance. • Suspended silt, and cesium adsorbed by silt, is deposited farther in the reservoir. • Suspended clay, and cesium adsorbed by clay, travels the farthest into the reservoir. With sufficient time, the dissolved cesium reaches the downstream end of the reservoir. This preliminary modeling also suggests the possibility of a suitable dam operation to control the cesium migration farther downstream from the dam. JAEA has been sampling in the Ogi Dam Reservoir, but these data were not yet available for the current model calibration and validation for this reservoir. Nonetheless these preliminary FLESCOT modeling results were qualitatively valid and confirmed the applicability of the FLESCOT code to the Ogi Dam Reservoir, and in general to other reservoirs in

  15. Numerical simulation of catastrophic flood: the case study of hypothetical failure of the Bielkowo hydro-power plant reservoir

    NASA Astrophysics Data System (ADS)

    Szydłowski, Michał; Szpakowski, Wojciech; Zima, Piotr

    2013-10-01

    The numerical modeling of flood wave propagation following the hypothetical breaks of the embankments of the Bielkowo hydro-power plant storage reservoir (Kolbudy II Reservoir) on the Radunia River in Poland has been presented. The results of computations were used to estimate the parameters of the flood waves, which are indispensable for the flood zone determination and mapping and then for the flood risk analysis. When estimating the reach and area of the inundation, related to the embankments failures, digital terrain model, and mathematical model of flood wave propagation are necessary. For the numerical simulations of flood, the mathematical model of free surface, two-dimensional unsteady water flow was applied. Four locations of potential breaks of the reservoir embankments were considered. The computed flood zones were presented on the flood hazard maps. The maps have been used by the local authorities and the dam owner to manage the flood risk related to hydro-power plants operations on the Radunia River. This type of research has been done for the first time for the water plant managed by the ENERGA Elektrownie Straszyn.

  16. A simulation research on evaluation of development in shale oil reservoirs by near-miscible CO2 flooding

    NASA Astrophysics Data System (ADS)

    Lai, Fengpeng; Li, Zhiping; Fu, Yingkun; Yang, Zhihao; Li, Hong

    2015-08-01

    Shale oil is a key resource that could mitigate the impending energy shortage in the future. Despite its abundance in China, studies on shale oil are still at the preliminary stage. Shale oil development through CO2 flooding has been successfully implemented in the United States. Therefore, the mechanics of CO2 flooding in shale oil reservoirs should be investigated. This study applies a simulation method to evaluate the development efficiency of CO2 flooding in shale oil reservoirs. Near-miscible CO2 flooding can effectively develop shale oil. After 20 years, recovery could improve by up to 9.56% as a result of depletion development under near-miscible CO2 flooding with 0.5% pore volume gas injection. Horizontal well injection is better than vertical well injection in terms of sweep efficiency and recovery. Cyclic gas injection is superior to continuous gas injection because the former reduces gas channelling. Thus, the use of horizontal wells with near-miscible cyclic gas injections has the potential to effectively develop shale oil reservoirs.

  17. Assessing Reservoir Depositional Environments to Develop and Quantify Improvements in CO2 Storage Efficiency. A Reservoir Simulation Approach

    SciTech Connect

    Okwen, Roland; Frailey, Scott; Leetaru, Hannes; Moulton, Sandy

    2014-09-30

    The storage potential and fluid movement within formations are dependent on the unique hydraulic characteristics of their respective depositional environments. Storage efficiency (E) quantifies the potential for storage in a geologic depositional environment and is used to assess basinal or regional CO2 storage resources. Current estimates of storage resources are calculated using common E ranges by lithology and not by depositional environment. The objectives of this project are to quantify E ranges and identify E enhancement strategies for different depositional environments via reservoir simulation studies. The depositional environments considered include deltaic, shelf clastic, shelf carbonate, fluvial deltaic, strandplain, reef, fluvial and alluvial, and turbidite. Strategies considered for enhancing E include CO2 injection via vertical, horizontal, and deviated wells, selective completions, water production, and multi-well injection. Conceptual geologic and geocellular models of the depositional environments were developed based on data from Illinois Basin oil fields and gas storage sites. The geologic and geocellular models were generalized for use in other US sedimentary basins. An important aspect of this work is the development of conceptual geologic and geocellular models that reflect the uniqueness of each depositional environment. Different injection well completions methods were simulated to investigate methods of enhancing E in the presence of geologic heterogeneity specific to a depositional environment. Modeling scenarios included horizontal wells (length, orientation, and inclination), selective and dynamic completions, water production, and multiwell injection. A Geologic Storage Efficiency Calculator (GSECalc) was developed to calculate E from reservoir simulation output. Estimated E values were normalized to diminish their dependency on fluid relative permeability. Classifying depositional environments according to

  18. Geological and petrophysical characterization of the Ferron Sandstone for 3-D simulation of a fluvial-deltaic reservoir. Quarterly technical progress report, January 1--March 31, 1997

    SciTech Connect

    Allison, M.L.

    1997-04-01

    The objective of this project is to develop a comprehensive, interdisciplinary, and quantitative characterization of a fluvial-deltaic reservoir which will allow realistic inter-well and reservoir-scale modeling to be constructed for improved oil-field development in similar reservoirs world-wide. The geological and petrophysical properties of the Cretaceous Ferron Sandstone in east-central Utah will be quantitatively determined. Both new and existing data will be integrated into a three-dimensional representation of spatial variations in porosity, storativity, and tensorial rock permeability at a scale appropriate for inter-well to regional-scale reservoir simulation. Results could improve reservoir management through proper infill and extension drilling strategies, reduction of economic risks, increased recovery from existing oil fields, and more reliable reserve calculations. Transfer of the project results to the petroleum industry is an integral component of the project.

  19. Geological and petrophysical characterization of the Ferron Sandstone for 3-D simulation of a fluvial-deltaic reservoir. Quarterly progress report, January 1, 1997--March 31, 1997

    SciTech Connect

    Allison, M.L.

    1997-05-01

    The objective of this project is to develop a comprehensive, interdisciplinary, and quantitative characterization of a fluvial-deltaic reservoir which will allow realistic inter-well and reservoir-scale modeling to be constructed for improved oil-field development in similar reservoirs world-wide. The geological and petrophysical properties of the Cretaceous Ferron Sandstone in east-central Utah will be quantitatively determined. Both new and existing data will be integrated into a three-dimensional representation of spatial variations in porosity, storativity, and tensorial rock permeability at a scale appropriate for inter-well to regional-scale reservoir simulation. Results could improve reservoir management through proper infill and extension drilling strategies, reduction of economic risks, increased recovery from existing oil fields, and more reliable reserve calculations. Transfer of the project results to the petroleum industry is an integral component of the project.

  20. Geological and petrophysical characterization of the Ferron Sandstone for 3-D simulation of a fluvial-deltaic reservoir. Technical progress report, January 1, 1995--March 31, 1995

    SciTech Connect

    Allison, M.L.

    1995-05-02

    The objective of this project is to develop a comprehensive, interdisciplinary, and quantitative characterization of a fluvial-deltaic reservoir which will allow realistic inter-well and reservoir-scale modeling to be developed for improved oil-field development in similar reservoirs world-wide. The geological and petrophysical properties of the Cretaceous Ferron Sandstone in east-central Utah will be quantitatively determined. Both new and existing data will be integrated into a three-dimensional representation of spatial variations in porosity, storativity, and tensorial rock permeability at a scale appropriate for inter-well to regional-scale reservoir simulation. Results could improve reservoir management through proper infill and extension drilling strategies, reduction of economic risks, increased recovery from existing oil fields, and more reliable reserve calculations. Transfer of the project results to the petroleum industry is an integral component of the project.

  1. Geological and petrophysical characterization of the Ferron Sandstone for 3-D simulation of a fluvial-deltaic reservoir. Quarterly report, July 1--September 30, 1994

    SciTech Connect

    Allison, M.L.

    1994-10-30

    The objective of this project is to develop a comprehensive, interdisciplinary, and quantitative characterization of a fluvial-deltaic reservoir which will allow realistic inter-well and reservoir-scale modeling to be developed for improved oil-field development in similar reservoirs world-wide. The geological and petrophysical properties of the Cretaceous Ferron Sandstone in east-central Utah will be quantitatively determined. Both new and existing data will be integrated into a 3-D representation of spatial variations in porosity, storativity, and tensorial rock permeability at a scale appropriate for inter-well to regional-scale reservoir simulation. Results could improve reservoir management through proper infill and extension drilling strategies, reduction of economic risks, increased recovery from existing oil fields, and more reliable reserve calculations. Transfer of the project results to the petroleum industry is an integral component of the project.

  2. Geological and petrophysical characterization of the Ferron Sandstone for 3-D simulation of a fluvial-deltaic reservoir. Technical progress report, July 1, 1996--September 30, 1996

    SciTech Connect

    Allison, M.L.

    1996-10-01

    The objective of this project is to develop a comprehensive, interdisciplinary, and quantitative characterization of a fluvial-deltaic reservoir which will allow realistic inter-well and reservoir-scale modeling to be constructed for improved oil-field development in similar reservoirs world-wide. The geological and petrophysical properties of the Cretaceous Ferron Sandstone in east-central Utah will be quantitatively determined. Both new and existing data will be integrated into a three-dimensional representation of spatial variations in porosity, storativity, and tensorial rock permeability at a scale appropriate for inter-well to regional-scale reservoir simulation. Results could improve reservoir management through proper infill and extension drilling strategies, reduction of economic risks, increased recovery from existing oil fields, and more reliable reserve calculations. Transfer of the project results to the petroleum industry is an integral component of the project.

  3. Mathematical simulation of temperatures in deep impoundments: verification tests of the Water Resources Engineers, Inc. model - Horsetooth and Flaming Gorge Reservoirs

    USGS Publications Warehouse

    King, D.L.; Sartoris, Jim J.

    1973-01-01

    Successful use of predictive mathematical models requires verification of the accuracy of the models by applying them to existing situations where the prediction can be compared with reality. A Corps of Engineers' modification of a deep reservoir thermal stratification model developed by Water Resources Engineers, Inc., was applied to two existing Bureau of Reclamation reservoirs for verification. Diffusion coefficients used for the Corps' Detroit Reservoir were found to apply to Horsetooth Reservoir in Colorado, for which very food computer input data were available. The Detroit diffusion coefficients gave a reasonable simulation of Flaming Gorge Reservoir in Wyoming and Utah, which has very complex and variable physical characteristics and for which only average-quality computer input data were available.

  4. Geological and petrophysical characterization of the Ferron Sandstone for 3-D simulation of a fluvial-deltaic reservoir. Quarterly progress report, October 1, 1997--December 31, 1997

    SciTech Connect

    Allison, M.L.

    1998-01-01

    The objective of this project is to develop a comprehensive, interdisciplinary, and quantitative characterization of a fluvial-deltaic reservoir which will allow realistic inter-well and reservoir-scale modeling to be constructed for improved oil-field development in similar reservoirs world-wide. The geological and petrophysical properties of the Cretaceous Ferron Sandstone in east-central Utah will be quantitatively determined. Both new and existing data will be integrated into a three-dimensional representation of spatial variations in porosity, storativity, and tensorial rock permeability at a scale appropriate for inter-well to regional-scale reservoir simulation. Results could improve reservoir management through proper infill and extension drilling strategies, reduction of economic risks, increased recovery from existing oil fields, and more reliable reserve calculations. Transfer of the project results to the petroleum industry is an integral component of the project.

  5. iTOUGH2-EOS1SC. Multiphase Reservoir Simulator for Water under Sub- and Supercritical Conditions. User's Guide

    SciTech Connect

    Magnusdottir, Lilja; Finsterle, Stefan

    2015-03-01

    Supercritical fluids exist near magmatic heat sources in geothermal reservoirs, and the high enthalpy fluid is becoming more desirable for energy production with advancing technology. In geothermal modeling, the roots of the geothermal systems are normally avoided but in order to accurately predict the thermal behavior when wells are drilled close to magmatic intrusions, it is necessary to incorporate the heat sources into the modeling scheme. Modeling supercritical conditions poses a variety of challenges due to the large gradients in fluid properties near the critical zone. This work focused on using the iTOUGH2 simulator to model the extreme temperature and pressure conditions in magmatic geothermal systems.

  6. Toward optical signal processing using photonic reservoir computing.

    PubMed

    Vandoorne, Kristof; Dierckx, Wouter; Schrauwen, Benjamin; Verstraeten, David; Baets, Roel; Bienstman, Peter; Van Campenhout, Jan

    2008-07-21

    We propose photonic reservoir computing as a new approach to optical signal processing in the context of large scale pattern recognition problems. Photonic reservoir computing is a photonic implementation of the recently proposed reservoir computing concept, where the dynamics of a network of nonlinear elements are exploited to perform general signal processing tasks. In our proposed photonic implementation, we employ a network of coupled Semiconductor Optical Amplifiers (SOA) as the basic building blocks for the reservoir. Although they differ in many key respects from traditional software-based hyperbolic tangent reservoirs, we show using simulations that such a photonic reservoir can outperform traditional reservoirs on a benchmark classification task. Moreover, a photonic implementation offers the promise of massively parallel information processing with low power and high speed. PMID:18648434

  7. Improved recovery from Gulf of Mexico reservoirs. Volume III (of 4): Characterization and simulation of representative resources. Final report, February 14, 1995--October 13, 1996

    SciTech Connect

    Kimbrell, W.C.; Bassiouni, Z.A.; Bourgoyne, A.T.

    1997-01-13

    Significant innovations have been made in seismic processing and reservoir simulation. In addition, significant advances have been made in deviated and horizontal drilling technologies. Effective application of these technologies along with improved integrated resource management methods offer opportunities to significantly increase Gulf of Mexico production, delay platform abandonments, and preserve access to a substantial remaining oil target for both exploratory drilling and advanced recovery processes. In an effort to illustrate the impact that these new technologies and sources of information can have upon the estimates of recoverable oil in the Gulf of Mexico, additional and detailed data was collected for two previously studied reservoirs: a South March Island reservoir operated by Taylor Energy and Gulf of Mexico reservoir operated by Mobil, whose exact location has been blind-coded at their request, and an additional third representative reservoir in the Gulf of Mexico, the KEKF-1 reservoir in West Delta Block 84 Field. The new data includes reprocessed 2-D seismic data, newly acquired 3-D data, fluid data, fluid samples, pressure data, well test data, well logs, and core data/samples. The new data was used to refine reservoir and geologic characterization of these reservoirs. Further laboratory investigation also provided additional simulation input data in the form of PVT properties, relative permeabilities, capillary pressures, and water compatibility. Geologic investigations were also conducted to refine the models of mud-rich submarine fan architectures used by seismic analysts and reservoir engineers. These results were also used, in part, to assist in the recharacterization of these reservoirs.

  8. Geology and petrophysical characterization of the Ferron Sandstone for 3-D simulation of a fluvial-deltaic reservoir. Annual report, October 1, 1996--September 30, 1997

    SciTech Connect

    Chidsey, T.C. Jr.; Anderson, P.B.; Morris, T.H.; Dewey, J.A. Jr.; Mattson, A.; Foster, C.B.; Snelgrove, S.H.; Ryer, T.A.

    1998-05-01

    The objective of the Ferron Sandstone (Utah) project is to develop a comprehensive, interdisciplinary, quantitative characterization of a fluvial-deltaic reservoir to allow realistic interwell and reservoir-scale models to be developed for improved oil-field development in similar reservoirs world-wide. Both new and existing data is being integrated into a 3-D model of spatial variations in porosity, storativity, and tensorial rock permeability at a scale appropriate for inter-well to regional-scale reservoir simulation. Simulation results could improve reservoir management through proper infill and extension drilling strategies, reduction of economic risks, increased recovery from existing oil fields, and more reliable reserve calculations. The project is divided into four tasks: (1) regional stratigraphic analysis, (2) case studies, (3) reservoirs models, and (4) field-scale evaluation of exploration strategies. The primary objective of the regional stratigraphic analysis is to provide a more detailed interpretation of the stratigraphy and gross reservoir characteristics of the Ferron Sandstone as exposed in outcrop. The primary objective of the case-studies work is to develop a detailed geological and petrophysical characterization, at well-sweep scale or smaller, of the primary reservoir lithofacies typically found in a fluvial-dominated deltaic reservoir. Work on tasks 3 and 4 consisted of developing two- and three-dimensional reservoir models at various scales. The bulk of the work on these tasks is being completed primarily during the last year of the project, and is incorporating the data and results of the regional stratigraphic analysis and case-studies tasks.

  9. Forward-Inverse Adaptive Techniques for Reservoir Characterization and Simulation: Theory and Applications

    SciTech Connect

    Doss, S D; Ezzedine, S; Gelinas, R; Chawathe, A

    2001-06-11

    A novel approach called Forward-Inverse Adaptive Techniques (FIAT) for reservoir characterization is developed and applied to three representative exploration cases. Inverse modeling refers to the determination of the entire reservoir permeability under steady state single-phase flow regime, given only field permeability, pressure and production well measurements. FIAT solves the forward and inverse partial differential equations (PDEs) simultaneously by adding a regularization term and filtering pressure gradients. An implicit adaptive-grid, Galerkin, numerical scheme is used to numerically solve the set of PDEs subject to pressure and permeability boundary conditions. Three examples are presented. Results from all three cases demonstrate attainable and reasonably accurate solutions and, more importantly, provide insights into the consequences of data undersampling.

  10. Geological and petrophysical characterization of the Ferron Sandstone for 3-D simulation of a fluvial-deltaic reservoir. Technical progress report, April 1--June 30, 1995

    SciTech Connect

    Allison, M.L.

    1995-07-28

    The objective of this project is to develop a comprehensive, interdisciplinary, and quantitative characterization of a fluvial-deltaic reservoir which will allow realistic inter-well and reservoir-scale modeling to be constructed for improved oil-field development in similar reservoirs world-wide. The geological and petrophysical properties of the Cretaceous Ferron Sandstone in east-central Utah will be quantitatively determined. Both new and existing data will be integrated into a three-dimensional representation of spatial variations in porosity, storativity, and tensorial rock permeability at a scale appropriate for inter-well to regional-scale reservoir simulation. Results could improve reservoir management through proper infill and extension drilling strategies, reduction of economic risks, increased recovery from existing oil fields, and more reliable reserve calculations. Technical progress this quarter is divided into regional stratigraphy, case studies, stochastic modeling and fluid-flow simulation, and technology transfer activities. The regional stratigraphy of the Ferron Sandstone outcrop belt from Last Chance Creek to Ferron Creek is being described and interpreted. Photomosaics and a database of existing surface and subsurface data are being used to determine the extent and depositional environment of each parasequence, and the nature of the contacts with adjacent rocks or flow units. For the second field season, detailed geological and petrophysical characterization of the primary reservoir lithofacies typically found in a fluvial-dominated deltaic reservoir, is continuing at selected case-study areas.