Three dimensional simulation of fluid flow in X‐ray CT images of porous media

A numerical scheme is developed in order to simulate fluid flow in three dimensional (3‐D) microstructures. The governing equations for steady incompressible flow are solved using the semi‐implicit method for pressure‐linked equations (SIMPLE) finite difference scheme within a non‐staggered grid system that represents the 3‐D microstructure. This system allows solving the governing equations using only one computational cell. The numerical scheme is verified through simulating fluid flow in idealized 3‐D microstructures with known closed form solutions for permeability. The numerical factors affecting the solution in terms of convergence and accuracy are also discussed. These factors include the resolution of the analysed microstructure and the truncation criterion. Fluid flow in 2‐D X‐ray computed tomography (CT) images of real porous media microstructure is also simulated using this numerical model. These real microstructures include field cores of asphalt mixes, laboratory linear kneading compactor (LKC) specimens, and laboratory Superpave gyratory compactor (SGC) specimens. The numerical results for the permeability of the real microstructures are compared with the results from closed form solutions. Copyright © 2004 John Wiley & Sons, Ltd.     

Isotope hydrological study of mean transit time in the granitic Strengbach catchment (Vosges massif,France): application of the FlowPC model with modified input function

Measurements of ¹⁸O concentrations in precipitation, soil solution, spring and runoff are used to determine water transit time in the small granitic Strengbach catchment (0·8 km²; 883–1146 m above sea level) located in the Vosges Mountains of northeastern France. Water transit times were calculated by applying the exponential, exponential piston and dispersion models of the FlowPC program to isotopic input (rainfall) and output (spring and stream water) data sets during the period 1989–95. The input function of the model was modified compared with the former version of the model and estimated by a deterministic approach based on a simplified hydrological balance. The fit between observed and calculated output data showed marked improvements compared with results obtained using the initial version of the model. An exponential piston version of the model applied to spring water indicates a 38·5 month mean transit time, which suggests that the volume in the aquifer, expressed in water depth, is 2·4 m. A considerable thickness (>45 m) of fractured bedrock may be involved for such a volume of water to be stored in the aquifer. Copyright © 2005 John Wiley & Sons, Ltd.     

Reef exploration in the East Sengkang Basin,Sulawesi, Indonesia

Reconnaissance seismic shot in 1971/72 showed a number of well defined seismic anomalies within the East Sengkang Basin which were interpreted as buried reefs. Subsequent fieldwork revealed that Upper Miocene reefs outcropped along the southern margin of the basin. A drilling programme in 1975 and 1976 proved the presence of shallow, gas-bearing, Upper Miocene reefs in the northern part of the basin. Seismic acquisition and drilling during 1981 confirmed the economic significance of these discoveries, with four separate accumulations containing about 750 × 10⁹ cubic feet of dry gas in place at an average depth of 700 m. Kampung Baru is the largest field and contains over half the total, both reservoir quality and gas deliverability are excellent. Deposition in the East Sengkang Basin probably started during the Early Miocene. A sequence of Lower Miocene mudstones and limestones unconformably overlies acoustic basement which consists of Eocene volcanics. During the tectonically active Middle Miocene, deposition was interrupted by two periods of deformation and erosion. Carbonate deposition became established in the Late Miocene with widespread development of platform limestones throughout the East Sengkang Basin. Thick pinnacle reef complexes developed in the areas where reef growth could keep pace with the relative rise in sea level. Most reef growth ceased at the end of the Miocene and subsequent renewed clastic sedimentation covered the irregular limestone surface. Late Pliocene regression culminated in the Holocene with erosion. The Walanae fault zone, part of a major regional sinistral strike-slip system, separates the East and West Sengkang Basins. Both normal and reverse faulting are inferred from seismic data and post Late Pliocene reverse faulting is seen in outcrop.     

Petroleum exploration potential of Tamtsag Basin

The Tamtsag Basin is located in the extreme eastern portion of the Mongolia. The Basin and its counterpart in China (the Hailar Basin) are united a whole basin on the structural setting. In recent years, the Tamtsag Basin attracts more and more attention with the important exploration discovered in the 19th block by SOCO and in Hailar Basin of China. This paper discusses the exploration potential of Tamtsag Basin from the viewpoint of petroleum geology.     

Tephrochronology and tephrostratigraphy of two Pleistocene continental fossiliferous successions from central Italy

The stratigraphical context of two Middle Pleistocene fossiliferous palaeosols from Central Italy (Abruzzo and Tuscany) have been studied. Small mammals and molluscs occur in both palaeosols, which are covered by tephra layers that were analysed using an interdisciplinary approach. Application of fission‐track dating to apatites separated from the Case Picconetto tephra (Pescara, Abruzzo), yielded an age of 0.48 ± 0.04 Ma, indistinguishable from those previously determined for the Campani Quarry (Lower Valdarno, Tuscany) (0.46 ± 0.05 Ma and 0.48 ± 0.05 Ma). Geochemical and petrographic investigations indicate that these tephra originated from different volcanoes, the Alban Hills Volcanic Complex and the Vico Volcano (Latium) respectively. Small mammal and mollusc assemblages indicate different palaeoclimatic and palaeoenvironmental conditions for the Case Picconetto and Campani Quarry palaeosols. Warm and humid conditions can be inferred for the Campani Quarry site, whereas open and cold conditions can be inferred for Case Picconetto. On the basis of faunal data, fission‐track dates and attribution of tephra to specific volcanic eruptions, we suggest a correlation of these faunas with marine oxygen isotope stage 14 (Case Picconetto) and with marine oxygen isotope stage 11 (Campani Quarry), respectively. Copyright © 2003 John Wiley & Sons, Ltd.     

The response of an elastic half‐space under a momentary shear line impulse

The response of an ideal elastic half‐space to a line‐concentrated impulsive vector shear force applied momentarily is obtained by an analytical–numerical computational method based on the theory of characteristics in conjunction with kinematical relations derived across surfaces of strong discontinuities. The shear force is concentrated along an infinite line, drawn on the surface of the half‐space, while being normal to that line as well as to the axis of symmetry of the half‐space. An exact loading model is introduced and built into the computational method for this shear force. With this model, a compatibility exists among the prescribed applied force, the geometric decay of the shear stress component at the precursor shear wave, and the boundary conditions of the half‐space; in this sense, the source configuration is exact. For the transient boundary‐value problem described above, a wave characteristics formulation is presented, where its differential equations are extended to allow for strong discontinuities which occur in the material motion of the half‐space. A numerical integration of these extended differential equations is then carried out in a three‐dimensional spatiotemporal wavegrid formed by the Cartesian bicharacteristic curves of the wave characteristics formulation. This work is devoted to the construction of the computational method and to the concepts involved therein, whereas the interpretation of the resultant transient deformation of the half‐space is presented in a subsequent paper. Copyright © 2003 John Wiley & Sons, Ltd.     

Modelling of deformation response and chemo‐mechanical coupling in chalk

A constitutive relation is derived for describing the mechanical response of chalk. The approach is based on a phenomenological framework which employs chemo‐plasticity. The properties of the material are assumed to be affected by the physico‐chemical processes that occur through the interaction between the skeleton and the pore fluid. The underlying mechanism is discussed by invoking a micromechanical analysis. The performance of the framework is illustrated by examining the evolution of mechanical characteristics in the presence of different pore fluids. Copyright © 2006 John Wiley & Sons, Ltd.     

Use of multi‐platform,multi‐temporal remote‐sensing data for calibration of a distributed hydrological model: an application in the Arno basin,Italy

Images from satellite platforms are a valid aid in order to obtain distributed information about hydrological surface states and parameters needed in calibration and validation of the water balance and flood forecasting. Remotely sensed data are easily available on large areas and with a frequency compatible with land cover changes. In this paper, remotely sensed images from different types of sensor have been utilized as a support to the calibration of the distributed hydrological model MOBIDIC, currently used in the experimental system of flood forecasting of the Arno River Basin Authority. Six radar images from ERS‐2 synthetic aperture radar (SAR) sensors (three for summer 2002 and three for spring–summer 2003) have been utilized and a relationship between soil saturation indexes and backscatter coefficient from SAR images has been investigated. Analysis has been performed only on pixels with meagre or no vegetation cover, in order to legitimize the assumption that water content of the soil is the main variable that influences the backscatter coefficient. Such pixels have been obtained by considering vegetation indexes (NDVI) and land cover maps produced by optical sensors (Landsat‐ETM). In order to calibrate the soil moisture model based on information provided by SAR images, an optimization algorithm has been utilized to minimize the regression error between saturation indexes from model and SAR data and error between measured and modelled discharge flows. Utilizing this procedure, model parameters that rule soil moisture fluxes have been calibrated, obtaining not only a good match with remotely sensed data, but also an enhancement of model performance in flow prediction with respect to a previous calibration with river discharge data only. Copyright © 2006 John Wiley & Sons, Ltd.     

Chemical alteration of tephra in the depositional environment: theoretical stability modelling

The study of the chemical stability of vitreous material in aqueous media is well‐established. There has to date been little consideration of the implications of variations in the chemical durability of tephra in Quaternary tephrochronology. Chemical alteration can take the form of cationic leaching from the matrix, or complete destruction of the silica network, either of which could constrain the ability to chemically identify distal tephra. Here we apply established models of vitreous durability to the published chemical analyses of a large number of Icelandic tephras in order to predict their relative durabilities under equivalent conditions. This suggests that some important tephras have relatively poor chemical stability, and that rhyolitic tephras are, in general, more stable than basaltic. We conclude that tephras should be expected to show predictable differential chemical stability in the post‐depositional environment. Copyright © 2003 John Wiley & Sons, Ltd.     

Mezen Bay—a key area for understanding Weichselian glaciations in northern Russia

Sediment successions in coastal cliffs around Mezen Bay, southeastern White Sea, record an unusually detailed history of former glaciations, interstadial marine and fluvial events from the Weichselian. A regional glaciation model for the Weichselian is based on new data from the Mezen Bay area and previously published data from adjacent areas. Following the Mikulinian (Eemian) interglacial a shelf‐centred glaciation in the Kara Sea is reflected in proglacial conditions at 100–90 ka. A local ice‐cap over the Timan ridge existed between 75 and 65 ka. Renewed glaciation in the Kara Sea spread southwestwards around 60 ka only, interrupted by a marine inundation, before it advanced to its maximum position at about 55–50 ka. After a prolonged ice‐free period, the Scandinavian ice‐sheet invaded the area from the west and terminated east of Mezen Bay about 17 ka. The previously published evidence of a large ice‐dammed lake in the central Arkhangelsk region, Lake Komi, finds no support in this study. Copyright © 2003 John Wiley & Sons, Ltd.