Successful application of multiscale methods in a real reservoir simulator environment

For the past 10 years or so, a number of so-called multiscale methods have been developed as an alternative approach to upscaling and to accelerate reservoir simulation. The key idea of all these methods is to construct a set of prolongation operators that map between unknowns associated with cells in a fine grid holding the petrophysical properties of the geological reservoir model and unknowns on a coarser grid used for dynamic simulation. The prolongation operators are computed numerically by solving localized flow problems, much in the same way as for flow-based upscaling methods, and can be used to construct a reduced coarse-scale system of flow equations that describe the macro-scale displacement driven by global forces. Unlike effective parameters, the multiscale basis functions have subscale resolution, which ensures that fine-scale heterogeneity is correctly accounted for in a systematic manner. Among all multiscale formulations discussed in the literature, the multiscale restriction-smoothed basis (MsRSB) method has proved to be particularly promising. This method has been implemented in a commercially available simulator and has three main advantages. First, the input grid and its coarse partition can have general polyhedral geometry and unstructured topology. Secondly, MsRSB is accurate and robust when used as an approximate solver and converges relatively fast when used as an iterative fine-scale solver. Finally, the method is formulated on top of a cell-centered, conservative, finite-volume method and is applicable to any flow model for which one can isolate a pressure equation. We discuss numerical challenges posed by contemporary geomodels and report a number of validation cases showing that the MsRSB method is an efficient, robust, and versatile method for simulating complex models of real reservoirs.     

Erosion-driven uplift of the modern Central Alps

We present a compilation of data of modern tectono-geomorphic processes in the Central European Alps which suggest that observed rock uplift is a response to climate-driven denudation. This interpretation is predominantly based on the recent quantification of basin-averaged Late Holocene denudation rates that are so similar to the pattern and rates of rock uplift rates as determined by geodetic leveling. Furthermore, a GPS data-based synthesis of Adriatic microplate kinematics suggests that the Central Alps are currently not in a state of active convergence. Finally, we illustrate that the Central Alps have acted as a closed system for Holocene redistribution of sediment in which the peri-Alpine lakes have operated as a sink for the erosional products of the inner Central Alps.While various hypotheses have been put forward to explain Central Alpine rock uplift (e.g. lithospheric forcing by convergence, mantle processes, or ice melting) we show with an elastic model of lithospheric deformation, that the correlation between erosion and rock uplift rates reflects a positive feedback between denudation and the associated isostatic response to unloading. Thus, erosion does not passively respond to advection of crustal material as might be the case in actively converging orogens. Rather, we suggest that the geomorphic response of the Alpine topography to glacial and fluvial erosion and the resulting disequilibrium for modern channelized and associated hillslope processes explains much of the pattern of modern denudation and hence rock uplift. Therefore, in a non-convergent orogen such as the Central European Alps, the observed vertical rock uplift is primarily a consequence of passive unloading due to erosion.     

Paleoclimatic changes in southern Africa during the intensification of Northern Hemisphere glaciation: evidence from ODP Leg 175 Site 1085

Pliocene age sediments from Ocean Drilling Program Leg 175, Site 1085-A and B in the Cape Basin were analyzed to investigate the impact of the intensification of Northern Hemisphere glaciation (INHG) on the South Atlantic Benguela Current system from 4 to 2 Ma. Proxies for productivity (concentrations and mass accumulation rates of total organic carbon, carbon to nitrogen ratios, percent calcium carbonate, and percent biogenic silica) as well as weight percent sand (a proxy for preservation or winnowing) peak at 3.2, 3.0, 2.4, and 2.25 Ma. Normative calculations of allied trace and major elemental determinations indicate synchronous increases in productivity peaks, as well as high concentrations and accumulations of terrigenous sediments. Coeval increases in hemipelagic sedimentation and productivity indicators could be the result of enhanced eolian sedimentation resulting from strengthened winds, leading to elevated rates of upwelling and enhanced productivity. However, rapid burial, as indicated by high sedimentation rates, could also enhance preservation. The very high concentrations (>30%) and accumulations (up to 60 g/cm²/kyr) limit the likelihood that eolian sedimentation was the only transport mechanism, invoking an additional fluvial source. Rapid burial by either eolian or fluvial transport links these intervals of enhanced preservation and productivity with continental climate changes resulting from (1) increased winds and/or dust availability due to higher aridity in the Namibia/northern South Africa region; (2) lowered sea-level related to increased ice volume; (3) increased sediment load due to wetter conditions in the continental interior; or (4) some combination. Peaks at 3.2, 2.4 and 2.25 Ma are coincident with maximum precession, suggesting a link between hemipelagic sedimentation and enhanced monsoonal circulation over southern Africa. The Site 1085 sedimentary record during the INHG seems to be controlled by low-latitude processes linked to precession rather than hig-latitude processes.     

The cosmic ray energy spectrum between 10 and 10 eV

The energy spectrum of cosmic rays with primary energies between 10¹⁴ eV and 10¹⁶ eV has been studied with the CASA-MIA air shower array. The measured differential energy spectrum is a power law (dj/dE ∝ E^−y) with spectral indices γ of 2.66±0.02 below approximately 10¹⁵ eV and 3.00±0.05 above. A new method is used for measuring primary energy derived from ground-based data in a compositionally insensitive way. In contrast with some previous reports, the “knee” of the energy spectrum does not appear sharp, but rather a smooth transition over energies from 10¹⁵ eV to 3.0 × 10¹⁵ eV.     

Numerical modeling of ocean hydrodynamics with variational assimilation of observational data

Models and methods of the numerical modeling of ocean hydrodynamics dating back to the pioneering works of A.S. Sarkisyan are considered, with emphasis on the formulation of problems and algorithms of mathematical modeling and the four-dimensional variational assimilation of observational data. An algorithm is proposed for studying the sensitivity of the optimal solution to observational data errors in a seasurface temperature assimilation problem in order to retrieve heat fluxes on the surface. An example of a solution of the optimal problem of the World Ocean hydrodynamics with the assimilation of climatic temperature and salinity observations is offered.     

Comparative assessment of baseline concentration of the heavy metals in the soils of Tehran (Iran) with the comprisable reference data

Industrial soil systems have become complex due to human activities, as they disturb the natural development of soil and add hazardous elements. Hence, there is a growing public concern over the accumulation of heavy metals in soil due to the rapid industrial development during the last decade in Iran. The objectives of the present study are to carry out comparative assessment of the heavy metals in soils of the Chitgar Industrial Area Tehran (Iran) with average world guidelines and to evaluate enrichment and normalized scatter coefficient. In this work the concentrations of seven heavy metals (Cr, Co, Cd, Cu, Pb, Ni and Zn) from the surface soils of the industrial area were quantified, based on 210 samples collected from 70 sampling stations in three consecutive seasons May 2007, November 2007 and May 2008 using 0.5 × l Km square mesh. Common comparison methods were employed to assess the quality of soils and guideline values were used to represent a desired level of elements. Comparison variables and enrichment factor showed that Pb has high-level value while normalized scatter coefficient demonstrates that Cd increases in soils more rapidly as compared to other elements. This study revealed that normalized scatter coefficient can be effectively used to evaluate soil pollution and is independent of the past.     

Climate Change Impacts for the Conterminous USA: An Integrated Assessment

Here we simulate dryland agriculture in the United States in order to assess potential future agricultural production under a set of general circulation model (GCM)-based climate change scenarios. The total national production of three major grain crops—corn, soybeans, and winter wheat—and two forage crops—alfalfa and clover hay—is calculated for the actual present day core production area (CPA) of each of these crops. In general, higher global mean temperature (GMT) reduces production and higher atmospheric carbon dioxide concentration ([CO₂]) increases production. Depending on the climatic change scenarios employed overall national production of the crops studied changes by up to plus or minus 25% from present-day levels. Impacts are more significant regionally, with crop production varying by greater than ±50% from baseline levels. Analysis of currently possible production areas (CPPAs) for each crop indicates that the regions most likely to be affected by climate change are those on the margins of the areas in which they are currently grown. Crop yield variability was found to be primarily influenced by local weather and geographic features rather than by large-scale changes in climate patterns and atmospheric composition. Future US agronomic potential will be significantly affected by the changes in climate projected here. The nature of the crop response will depend primarily on to what extent precipitation patterns change and also on the degree of warming experienced.     

Structural control of Au–Pd mineralization (Jacutinga): An example from the Cauê Mine, Quadrilátero Ferrífero, Brazil

Precambrian banded iron formations (BIFs) in the Quadrilátero Ferrífero host a special kind of Au–Pd mineralization known as Jacutinga. The main orebodies are hosted within the Cauê Syncline, a SW-verging fold that involves Paleoproterozoic metasedimentary rocks in the Itabira District, a regional synclinorium with BIFs in the core of synclinal folds in the northeastern part of Quadrilatéro Ferrífero, Minas Gerais. Structural analysis reveals two important features of the district: the polydeformed character of the rocks and the importance of brittle structures in the control of the orebodies. Two deformational events are recognized in this area. The first event developed the main foliation, S1, that is the enveloping surface of the Cauê Syncline. The second event is better defined in the northern boundary of the structure where it is represented by a right-lateral wrench fault zone that has developed a foliation, S2, that truncates S1. This wrench fault was also responsible for the development of a system of fractures (Frm) that host the Au–Pd mineralization. The auriferous bodies of Cauê Syncline (Y, X, Área Central, Aba Norte, Noroeste and Aba Leste/Aba Leste Inferior) were generated during this second event. Shear fractures (R, R′ and P) and tension fractures (T) developed in response to the wrench fault system under brittle conditions. The best-developed, and most commonly mineralized fractures are R and T in all auriferous bodies. Elsewhere, the best mineralization occurs in the contacts of hematite bodies (soft/hard) and intrusive rocks with fractured itabirites. Other mineralization (Aba Norte, Área Central and X) is hosted on the contacts of other units.A system of fractures, as well as their intersections, thus represents the structural control on Jacutinga bodies and is responsible for the geometry of the orebodies. Of importance, there is no control by mineral/stretching lineations, fold axes and other ductile structure on the geometry and plunge of the orebodies.     

Iterative implicit integration procedure for hybrid simulation of large nonlinear structures

A fully implicit iterative integration procedure is presented for local and geographically distributed hybrid simulation of the seismic response of complex structural systems with distributed nonlinear behavior. The purpose of this procedure is to seamlessly incorporate experimental elements in simulations using existing fully implicit integration algorithms designed for pure numerical simulations. The difficulties of implementing implicit integrators in a hybrid simulation are addressed at the element level by introducing a safe iteration strategy and using an efficient procedure for online estimation of the experimental tangent stiffness matrix. In order to avoid physical application of iterative displacements, the required experimental restoring force at each iteration is estimated from polynomial curve fitting of recent experimental measurements. The experimental tangent stiffness matrix is estimated by using readily available experimental measurements and by a classical diagonalization approach that reduces the number of unknowns in the matrix. Numerical and hybrid simulations are used to demonstrate that the proposed procedure provides an efficient method for implementation of fully implicit numerical integration in hybrid simulations of complex nonlinear structures. The hybrid simulations presented include distributed nonlinear behavior in both the numerical and experimental substructures. Copyright © 2010 John Wiley & Sons, Ltd.