Geothermal reservoir engineering research

The Stanford University research program on the study of stimulation and reservoir engineering of geothermal resources commenced as an interdisciplinary program in September, 1972. The broad objectives of this program have been: (1) the development of experimental and computational data to evaluate the optimum performance of fracture-stimulated geothermal reservoirs; (2) the development of a geothermal reservoir model to evaluate important thermophysical, hydrodynamic, and chemical parameters based on fluid-energy-volume balances as part of standard reservoir engineering practice; and (3) the construction of a laboratory model of an explosion-produced chimney to obtain experimental data on the processes of in-place boiling, moving flash fronts, and two-phase flow in porous and fractured hydrothermal reservoirs

Optimized energy calculation in lattice systems with long-range interactions

We discuss an efficient approach to the calculation of the internal energy in numerical simulations of spin systems with long-range interactions. Although, since the introduction of the Luijten-Bl\"ote algorithm, Monte Carlo simulations of these systems no longer pose a fundamental problem, the energy calculation is still an O(N^2) problem for systems of size N. We show how this can be reduced to an O(N logN) problem, with a break-even point that is already reached for very small systems. This allows the study of a variety of, until now hardly accessible, physical aspects of these systems. In particular, we combine the optimized energy calculation with histogram interpolation methods to investigate the specific heat of the Ising model and the first-order regime of the three-state Potts model with long-range interactions.Comment: 10 pages, including 8 EPS figures. To appear in Phys. Rev. E. Also available as PDF file at http://www.cond-mat.physik.uni-mainz.de/~luijten/erikpubs.htm

A study of Dykstra-Parsons curves

The Dykstra-Parsons method for prediction of oil recovery by water flooding is a well known technique which has been used by the petroleum industry since 1945. The present work carries their study further, solving the same problem of calculating coverage for certain values of permeability variation having water-oil-ratio and mobility ratio as fixed parameters. The work herein, instead of using 50 layers, uses 200. Also a more precise theoretical approach to the problem is given. Because of these differences the resulting curves are slightly modified. In a second part, the authors deal with empirical simplifications with considerable success. The idea was to collapse the data and curves obtained in the first part into a single curve which covers most of the range of variables commonly seen in reservoir displacements

Universal and Non-Universal First-Passage Properties of Planar Multipole Flows

The dynamics of passive Brownian tracer particles in steady two-dimensional potential flows between sources and sinks is investigated. The first-passage probability, $p(t)$, exhibits power-law decay with a velocity-dependent exponent in radial flow and an order-dependent exponent in multipolar flows. For the latter, there also occur diffusive ``echo'' shoulders and exponential decays associated with stagnation points in the flow. For spatially extended dipole sinks, the spatial distribution of the collected tracer is independent of the overall magnitude of the flow field.Comment: 7 pages, LaTe

Techno-economic and risk evaluation of a thermal recovery project

Field production data were studied, to derive an overall energy balance for the steamflood, to calculate the steamflood capture efficiency and predict future steamflood performance. Heat-losses due to produced fluids were also calculated. Predicted production schedules from the model were history-matched with field production data The reservoir parameters (porosity, {phi}, net thickness, h{sub n}, initial oil saturation, S{sub oi}, and residual oil saturation, S{sub or}) were evaluated statistically using both Gaussian and triangular distributions. These resulted in distributed recovery predictions. The Gaussian distributions behaved as predicted; but of great importance, the skewed triangular distributions also behaved in much the same manner. The results fit closely with predictions using logical formulas to predict expected values, peak values and standard variations of recoveries. This result is important, for it indicates that complete Monte-Carlo simulations may not be necessary. All steamflood calculations were carried out using a PC-based spreadsheet program. The major results were as follows: The capture efficiency of the Wilmington steamflood was calculated at 60%. This is an acceptable value, taking into account the reservoir geometry and history. The calculated heat balance showed high heat-loss to adjacent formations and through produced fluids. Of the cumulative heat injected at the time of the study, 21% had been lost to vertical conduction and 21% through produced fluids. Predicted production schedules indicated that up to 43% of the oil in place (at steamflood initiation) could be recovered by the steamflood

CT imaging techniques for two-phase and three-phase in-situ saturation measurements

The aim of this research is to use the SUPRI 3D steam injection laboratory model to establish a reliable method for 3-phase in-situ saturation measurements, and thereafter investigate the mechanism of steamflood at residual oil saturation. Demiral et al. designed and constructed a three dimensional laboratory model that can be used to measure temperature, pressure and heat loss data. The model is also designed so that its construction materials are not a limiting factor for CT scanning. We have used this model for our study. In this study, we saturated the model with mineral oil, and carried out waterflood until residual oil saturation. Steamflood was then carried out. A leak appeared at the bottom of the model. Despite this problem, the saturation results, obtained by using 2-phase and 3-phase saturation equations and obtained from the Cat scanner, were compared with the saturations obtained from material balance. The errors thus obtained were compared with those obtained by an error analysis carried out on the saturation equations. This report gives details of the experimental procedures, the data acquisition and data processing computer programs, and the analysis of a steamflood experiment carried out at residual oil saturation

Approximate Quantum Fourier Transform and Decoherence

We discuss the advantages of using the approximate quantum Fourier transform (AQFT) in algorithms which involve periodicity estimations. We analyse quantum networks performing AQFT in the presence of decoherence and show that extensive approximations can be made before the accuracy of AQFT (as compared with regular quantum Fourier transform) is compromised. We show that for some computations an approximation may imply a better performance.Comment: 14 pages, 10 fig. (8 *.eps files). More information on http://eve.physics.ox.ac.uk/QChome.html http://www.physics.helsinki.fi/~kasuomin http://www.physics.helsinki.fi/~kira/group.htm

CT measurements of two-phase flow in fractured porous media

The simulation of flow in naturally fractured reservoirs commonly divides the reservoir into two continua - the matrix system and the fracture system. Flow equations are written presuming that the primary flow between grid blocks occurs through the fracture system and that the primary fluid storage is in the matrix system. The dual porosity formulation of the equations assumes that there is no flow between matrix blocks while the dual permeability formulation allows fluid movement between matrix blocks. Since most of the fluid storage is contained in the matrix, recovery is dominated by the transfer of fluid from the matrix to the high conductivity fractures. The physical mechanisms influencing this transfer have been evaluated primarily through numerical studies. Relatively few experimental studies have investigated the transfer mechanisms. Early studies focused on the prediction of reservoir recoveries from the results of scaled experiments on single reservoir blocks. Recent experiments have investigated some of the mechanisms that are dominant in gravity drainage situations and in small block imbibition displacements. The mechanisms active in multiphase flow in fractured media need to be further illuminated, since some of the experimental results appear to be contradictory. This report describes the design, construction, and preliminary results of an experiment that studies imbibition displacement in two fracture blocks. Multiphase (oil/water) displacements will be conducted at the same rate on three core configurations. The configurations are a compact core, a two-block system with a 1 mm spacer between the blocks, and a two-block system with no spacer. The blocks are sealed in epoxy so that saturation measurements can be made throughout the displacement experiments using a Computed Tomography (CT) scanner

Elucidation of Beta-Oxidation Pathways in Ralstonia Eutropha H16 by Examination of Global Gene Expression

Ralstonia eutropha H16 is capable of growth and polyhydroxyalkanoate production on plant oils and fatty acids. However, little is known about the triacylglycerol and fatty acid degradation pathways of this bacterium. We compare whole-cell gene expression levels of R. eutropha H16 during growth and polyhydroxyalkanoate production on trioleate and fructose. Trioleate is a triacylglycerol that serves as a model for plant oils. Among the genes of note, two potential fatty acid β-oxidation operons and two putative lipase genes were shown to be upregulated in trioleate cultures. The genes of the glyoxylate bypass also exhibit increased expression during growth on trioleate. We observed that single β-oxidation operon deletion mutants of R. eutropha could grow using palm oil or crude palm kernel oil as the sole carbon source, regardless of which operon was present in the genome, but a double mutant was unable to grow under these conditions. A lipase deletion mutant did not exhibit a growth defect in emulsified oil cultures but did exhibit a phenotype in cultures containing nonemulsified oil. Mutants of the glyoxylate shunt gene for isocitrate lyase were able to grow in the presence of oils, while a malate synthase (aceB) deletion mutant grew more slowly than wild type. Gene expression under polyhydroxyalkanoate storage conditions was also examined. Many findings of this analysis confirm results from previous studies by our group and others. This work represents the first examination of global gene expression involving triacylglycerol and fatty acid catabolism genes in R. eutropha.Malaysia-MIT Biotechnology Partnership Programm

Overview: Will Great Power Politics Threaten Arctic Sustainability?

The eight Arctic states (Canada, Denmark/Greenland, Finland, Iceland, Norway, Sweden, the Russian Federation, and the United States) wish to maintain a position of preeminence when it comes to dealing with matters of Arctic Ocean governance. However, major non-Arctic states, while recognizing the sovereign rights of the coastal states in their economic zones and on their continental shelves, have growing interests in the maritime Arctic relating to activities such as commercial shipping, oil and gas development, fishing and ship-based tourism. They are increasingly claiming to have a legitimate interest in being consulted when it comes to addressing matters relating to the governance of such activities. Many questions arise from these new formulations: What is the appropriate mechanism for introducing the concerns of the non-Arctic North Pacific countries in forums dealing with Arctic issues? Can they achieve a significant voice in the deliberations of the Arctic Council? Can they address these issues through other mechanisms