A system-theoretical approach to selective grid coarsening of reservoir models

From a system-theoretical point of view and for a given configuration of wells, there are only a limited number of degrees of freedom in the input–output dynamics of a reservoir system. This means that a large number of combinations of the state variables (pressure and saturation values) are not actually controllable and observable from the wells, and accordingly, they are not affecting the input–output behavior of the system. In an earlier publication, we therefore proposed a control-relevant upscaling methodology that uniformly coarsens the reservoir. Here, we present a control-relevant selective (i.e. non-uniform) coarsening (CRSC) method, in which the criterion for grid size adaptation is based on ranking the grid block contributions to the controllability and observability of the reservoir system. This multi-level CRSC method is attractive for use in iterative procedures such as computer-assisted flooding optimization for a given configuration of wells. In contrast to conventional flow-based coarsening techniques our method is independent of the specific flow rates or pressures imposed at the wells. Moreover the system-theoretical norms employed in our method provide tight upper bounds to the ‘input–output energy’ of the fine and coarse systems. These can be used as an a priori error-estimate of the performance of the coarse model. We applied our algorithm to two numerical examples and found that it can accurately reproduce results from the corresponding fine-scale simulations, while significantly speeding up the simulation.     

Glide systems of hematite single crystals in deformation experiments

The critical resolved shear stresses (CRSSs) of hematite crystals were determined in compression tests for r-twinning, c-twinning and {a}<m>-slip in the temperature range 25 °C to 400 °C, at 400 MPa confining pressure, and a strain rate of 10^− 5 s^− 1 by Hennig-Michaeli, Ch., Siemes, H., 1982. Experimental deformation of hematile crstals betwen 25 °C and 400 °C at 400 MPa confining pressure. In: Schreyer, W. (Ed.) High Pressure Research in Geoscience, Schweizerbart'sche Verlagsbuchhandlung, Stuttgart, p. 133–150. In the present contribution newly performed experiments on hematite single crystals at temperatures up to 800 °C at strain rates of 10^− 5 s^− 1 and 300 MPa confining pressure extends the knowledge about the CRSS of twin and slip modes. Optical observations, neutron diffraction goniometry, SEM forescatter electron images and electron backscatter diffraction are applied in order to identify the glide modes. Both twinning systems and {a}<m>-slip were confirmed by these methods. Besides the known glide systems the existence of the (c)<a>-slip system could be stated. Mechanical data establish that the CRSS of r-twinning decreases from 140 MPa at 25 °C to  5 MPa at 800 °C and for {a}<m>-slip from > 560 MPa at 25 °C to  40 MPa at 700 °C. At room temperature the CRSS for c-twinning is around 90 MPa and at 600 °C  60 MPa. The data indicate that the CRSSs above 200 °C seem to be between the values for r-twinning and {a}<m>-slip. For (c)<a>-slip only the CRSS at 600 °C could be evaluated to  60 MPa. Exact values are difficult to determine because other glide systems are always simultaneously activated.     

Preliminary diagnoses of four new littoral species of Marine Crustacea isopoda (family sphaeromatidae)

Preliminary diagnoses are given for four new species belonging to Exosphaeroma, Dynamenella, Dynamenopsis and Amphoroidea.     

The Nearby Field Galaxy Survey

We report on our observing programme to obtain integrated spectrophotometry, intermediate- and high-resolution major-axis spectra, and UBR surface photometry of a representative sample of ∼ 200galaxies in the nearby field. The main goal of this programme is to provide a comparison sample for high-redshift studies and to study the variation in star formation rates (SFRs), star formation history (SFH),excitation, metallicity, and internal kinematics over a wide range of galaxy luminosities and morphological types. In particular, we extend the work of Kennicutt (1992) to lower-luminosity systems. We present the main results of our analysis so far. In these proceedings, we condense the two atlases presented of 1) images and radial surface brightness profiles and colour profiles, and 2) images and integrated spectra into several example images, profiles and spectra, showing the general trends observed. For the original atlases we refer the reader to the electronic distribution on CDROM or as available on the Web at http://www.astro.rug.nl/∼nfgs/. This revised version was published online in July 2006 with corrections to the Cover Date.     

Preliminary Results on the 3-dimensional Seismic Structure of the Lithosphere under the USGS Central California Seismic Array

About 1500 readings of teleseismic P -time residuals obtained from the US Geological Survey seismograph network in central California have been used to obtain a three-dimensional image of seismic velocity anomalies for this area by the method of Aki, Christoffersson & Husebye We found that the California network is less suitable than the LASA and NORSAR arrays for this kind of studies because of its greater proportion of peripheral blocks in which the resolution is very poor for the stochastic inverse solution and the random error effect is severe for the generalized inverse solution. Nevertheless, the resultant velocity anomalies show a remarkable correlation with the San Andreas fault zone to a depth of 75 km. The anomaly pattern changes drastically as the depth exceeds 75 km, suggesting that the asthenosphere has been reached.     

Use of Spherical Trigonometry to Compute Near-Well Flow through Irregular Grid Block Boundaries

This paper concerns the computation of near-well flow in numerical reservoir simulation with unstructured grids. In particular, it uses spherical trigonometry to derive analytical expressions for the flow towards a well modeled as either a number of point sources or a constant-flux line source. The expression for the point source representation is based on projections of the grid block boundaries on spheres with unit radius around the sources. The expression for the line source is based on projection on a prolate spheroid. The computation of the surface area is done through transformation to prolate spheroidal coordinates and subsequent projection on a sphere at infinity. The point source expression for a single source is exact for grid block boundaries with straight edges; the line source expression is an approximation. Both representations are fully volume conserving, such that the sum of the fluxes through the grid block boundaries surrounding a source adds up exactly to the total source flow rate. Both representations can be used to accurately model complicated wells in the form of segments. The point source representation is simpler to implement and not necessarily less accurate than the line source representation.     

Sensitivity of last glacial maximum climate to sea ice conditions in the Nordic Seas

Published reconstructions of last glacial maximum (LGM) sea surface temperatures and sea ice extent differ significantly. We here test the sensitivity of simulated North Atlantic climates to two different reconstructions by using these reconstructions as boundary conditions for model experiments. An atmospheric general circulation model has been used to perform two simulations of the (LGM) and a modern-day control simulation. Standard (CLIMAP) reconstructions of sea ice and sea surface temperatures have been used for the first simulation, and a set of new reconstructions in the Nordic Seas/Northern Atlantic have been used for the second experiment. The new reconstruction is based on 158 core samples, and represents ice-free conditions during summer in the Nordic Seas, with accordingly warmer sea surface temperatures and less extensive sea ice during winter as well. The simulated glacial climate is globally 5.7 K colder than modern day, with the largest changes at mid and high latitudes. Due to more intense Hadley circulation, the precipitation at lower latitudes has increased in the simulations of the LGM. Relative to the simulation with the standard CLIMAP reconstructions, reduction of the sea ice in the North Atlantic gives positive local responses in temperature, precipitation and reduction of the sea level pressure. Only very weak signatures of the wintertime Icelandic Low occur when the standard CLIMAP sea surface temperature reconstruction is used as the lower boundary condition in LGM. With reduced sea ice conditions in the Nordic Seas, the Icelandic Low becomes more intense and closer to its present structure. This indicates that thermal forcing is an important factor in determining the strength and position of the Icelandic Low. The Arctic Oscillation is the most dominant large scale variability feature on the Northern Hemisphere in modern day winter climate. In the simulation of the LGM with extensive sea ice this pattern is significantly changed and represents no systematic large scale variability over the North Atlantic. Reduction of the North Atlantic sea ice extent leads to stronger variability in monthly mean sea level pressure in winter. The synoptic variability appears at a lower level in the simulation when standard reconstructions of the sea surface in the LGM are used. A closer inspection of storm tracks in this model experiment shows that that the synoptic lows follow a narrow band along the ice edge during winter. The trajectories of synoptic lows are not constrained to the sea ice edge to the same degree when the sea ice extent is reduced. Seasonally open waters in the Nordic Seas in the new reconstruction apparently act as a moisture source, consistent with the current understanding of the rapid growth of the Fennoscandian and Barents Ice Sheets, during the LGM. The signal from the intensified thermal forcing in the North Atlantic in Boreal winter is carried zonally by upper tropospheric waves, and thus generates non-local responses to the changed sea ice cover.