A new computational strategy for solving two‐phase flow in strongly heterogeneous poroelastic media of evolving scales

We develop a new computational methodology for solving two‐phase flow in highly heterogeneous porous media incorporating geomechanical coupling subject to uncertainty in the poromechanical parameters. Within the framework of a staggered‐in‐time coupling algorithm, the numerical method proposed herein relies on a Petrov–Galerkin postprocessing approach projected on the Raviart–Thomas space to compute the Darcy velocity of the mixture in conjunction with a locally conservative higher order finite volume discretization of the nonlinear transport equation for the saturation and an operator splitting procedure based on the difference in the time‐scales of transport and geomechanics to compute the effects of transient porosity upon saturation. Notable features of the numerical modeling proposed herein are the local conservation properties inherited by the discrete fluxes that are crucial to correctly capture the fingering patterns arising from the interaction between heterogeneity and nonlinear viscous coupling. Water flooding in a poroelastic formation subject to an overburden is simulated with the geology characterized by multiscale self‐similar permeability and Young modulus random fields with power‐law covariance structure. Statistical moments of the poromechanical unknowns are computed within the framework of a high‐resolution Monte Carlo method. Numerical results illustrate the necessity of adopting locally conservative schemes to obtain reliable predictions of secondary recovery and finger growth in strongly heterogeneous deformable reservoirs. Copyright © 2011 John Wiley & Sons, Ltd.     

A Markov Chain Model for Subsurface Characterization: Theory and Applications

This paper proposes an extension of a single coupled Markov chain model to characterize heterogeneity of geological formations, and to make conditioning on any number of well data possible. The methodology is based on the concept of conditioning a Markov chain on the future states. Because the conditioning is performed in an explicit way, the methodology is efficient in terms of computer time and storage. Applications to synthetic and field data show good results.     

Evaluation of heterogeneity and field-scale groundwater flow regime in a leaky till aquitard

Recent work in southern Ontario, Canada, demonstrates anomalously high vertical groundwater flow velocities (>1 m/year) through a thick (as much as 60 m), sandy silt till aquitard (Northern till), previously assumed to be of very low permeability (hydraulic conductivity <10^–10 m/s). Rapid recharge is attributed to the presence of fractures and sedimentary heterogeneities within the till, but the field-scale flow regime is poorly understood. This study identifies the nature of physical groundwater pathways through the till and provides estimates of the associated groundwater fluxes. The aquitard groundwater flow system is characterized by integrating details of the outcrop and subsurface sedimentary characteristics of the till with field-based hydrogeologic investigation and numerical modeling. Outcrop and subsurface data identify a composite internal aquitard stratigraphy consisting of tabular till beds (till elements) separated by laterally continuous sheet-like sands and gravels (interbeds) and boulder pavements. Individual till elements contain sedimentary heterogeneities, including discontinuous sand and gravel lenses, vertical sand dikes, and zones of horizontal and vertical fractures. Hydrogeologic field investigations indicate a three-layer aquitard flow system, consisting of upper and lower zones of more hydraulically active and heterogeneous till separated by a middle unit of relatively lower hydraulic conductivity. Groundwater pathways and fluxes in the till were evaluated using a two-dimensional aquitard/aquifer flow model which indicates a step-wise flow mechanism whereby groundwater moves alternately downward along vertical pathways (fractures, sedimentary dikes) and laterally along horizontal sand interbeds within the till. This model is consistent with observed hydraulic-head and isotope profiles, and the presence of tritiated pore waters at various depths throughout the till. Simulations suggest that a bulk aquitard vertical hydraulic conductivity on the order of 1×10^–9 m/s is required to reproduce observed hydraulic-head and tritium profiles. Electronic Publication     

横向不均匀介质中的面波及地球内部速度结构研究

近十年中地震面波的理论和应用研究有很大进展。由于资料增加、反演方法的改进及利用了振幅和波形等更多的面波信息，对于全球及区域上地幔速度结构及各向异性有了进一步的认识。同时对于显著不均匀三维介质中面波的传播特征，在理论探讨和实际观测方面取得了较多成果。本文对于面波研究的最新成果作了介绍。     

秦岭─大别造山带及其南北缘地震层析成像

利用秦岭─大别造山带及其毗邻地区３１０个地震台站记录到的区域地震２３６００条Ｐ波到时数据，重建了该区地壳和上地幔三维速度图像。结果表明：１．秦岭─大别造山带及其毗邻地区地壳和上地幔存在显著的横向不均匀性，直至１１０ｋｍ深度处依然明显。２．地壳上部的速度图像与地表地质构造密切相关：造山带隆起区显著高速；盆地及坳陷区明显低速。由速度鲜明对比勾勒出的秦岭─大别造山带南界基本上位于扬子北缘主边断裂带上。３．中地壳的速度图像表明，造山带内部的一些低速区对应于一些大型推覆构造。４．４０＋０ｋｍ深度处的速度图像反映了该区莫霍界面深度的起伏。大致以１０７°Ｅ为界，以东地区地壳厚度小于４０ｋｍ，以西地区大于４０ｋｍ，且呈现出往西地壳逐渐加厚的趋势。５．位于滦川、商县、丹凤的北秦岭构造带，上地幔顶部出现低速异常，异常速度值约为７．３９-７．５５ｋｍ／ｓ。结合地球物理测深的结果，可能是由下地壳、上地幔顶部的热过程所致。     

Free oscillations of a laterally heterogenous and anelastic earth

We calculate normal modes of the laterally heterogeneous and anelastic earth model by using the variational method to include the coupling of the modes due to the asphericity of the earth. If the aspherical anelasticity correlates with the heterogeneity of elastic velocity structure, the quality factorQ of the split singlets has a correlation with the eigenfrequency. This can cause a center frequency shift of the spectral peak with time. We perform a synthetic experiment to examine whether the magnitude of the shift can become an observable for the realistic lateral heterogeneity model of anelasticity. The result of the experiment reveals that the shift of the center frequency is consistent with the initial estimate for the fundamental spheroidal modes used in the experiment. We then examine the actual seismograms of the June 9, 1994, Bolivian earthquake to determine if this shift of center frequency can be observed. Although the amount of the center frequency shift of each multiplet is large, there is no consistent shift of the center frequency that is predicted in the synthetic experiment.     

A resumé of nonlocal transport theories

Analysis from a number of different perspectives has shown diffusion and dispersion in natural porous formations to generally be nonlocal in character, i.e., the mass balance involves integro-partial differential equations. Only in certain asymptotic limits do these laws localize to classical partial differential equations. Compiled within is a resume of nonlocal laws that our group has developed over the last few years for systems with physical, chemical and biological heterogeneity. Analytical tools used to obtain these laws are nonequilibrium and equilibrium statistical mechanics, and first-order spectral-perturbation methods. This paper is an expansion of the material presented at the Waterloo conference held in the memory of Dr. Unny.     

Deglacial land emergence and lateral upper-mantle heterogeneity in the Svalbard Archipelago—II. Extended results for high-resolution load models

In Paper I (Breuer & Wolf 1995), a preliminary interpretation of the postglacial land emergence observed at a restricted set of six locations in the Svalbard Archipelago was given. The study was based on a simple model of the Barents Sea ice sheet and suggested increases in lithosphere thickness and asthenosphere viscosity with increasing distance from the continental margin.
In the present paper, the newly developed high-resolution load model. BARENTS-2, and land-uplift observations from an extended set of 25 locations are used to study further the possibility of resolving lateral heterogeneity in the upper mantle below the northern Barents Sea. A comparison of the calculated and observed uplift values shows that the lithosphere thickness is not well resolved by the observations, although values above 110 km are most common for this parameter. In contrast to this, there are indications of a lateral variation of asthenosphere viscosity. Whereas values in the range 10¹⁸-10²⁰Pas are inferred for locations close to the continental margin, 10²⁰-10²¹ Pa s are suggested further away from the margin.
A study of the sensitivity of the values found for lithosphere thickness and asthenosphere viscosity to modifications of load model BARENTS-2 shows that such modifications can be largely accommodated by appropriate changes in lithosphere thickness, whereas the suggested lateral variation of asthenosphere viscosity is essentially unaffected. An estimate of the influence of the Fennoscandian. ice sheet leads to the conclusion that its neglect results in an underestimation of the thickness of the Barents Sea ice sheet by about 10 per cent.     

Horizontal layered structure with heterogeneity beneath continents and island arcs from particle orbits of long-period P waves

We investigate the particle orbits of long-period (about 20 s) P waves observed with the global seismic network. By analysing 84 three-component seismograms recorded at 25 stations from 60 earthquakes occurring beneath 300 km, we quantitatively evaluate the orbits by three sets of eigenvalues and eigenvectors, using a covariance matrix method. The eigenvalues for P waves recorded at stations located on continents are explained by the standard horizontal layered structure model (iasp91). On the other hand, the orbits observed at stations close to island arcs are affected not only by the horizontal layered structure but also by heterogeneity due to subducting plates, mantle diapirs and so on. On the basis of a single-scattering model for a plane P wave, we quantify the heterogeneities by an isotropic scattering coefficient g0. Fitting the theoretical eigenvalues to the observed ones, we estimate g0 for the crust and upper mantle beneath continents to be less than 0.0005 km^-1, and the mean g₀ for the structure beneath island arcs to be about 0.0015 to 0.003 km^-1.     

Towards knowledge-based geovisualisation using Semantic Web technologies: a knowledge representation approach coupling ontologies and rules

ABSTRACT

Geovisualisation is a knowledge-intensive art in which both providers and users need to possess a wide range of knowledge. Current syntactic approaches to presenting visualisation information lack semantics on the one hand, and on the other hand are too bespoke. Such limitations impede the transfer, interpretation, and reuse of the geovisualisation knowledge. In this paper, we propose a knowledge-based approach to formally represent geovisualisation knowledge in a semantically-enriched and machine-readable manner using Semantic Web technologies. Specifically, we represent knowledge regarding cartographic scale, data portrayal and geometry source, which are three key aspects of geovisualisation in the contemporary web mapping era, coupling ontologies and semantic rules. The knowledge base enables inference for deriving the corresponding geometries and portrayals for visualisation under different conditions. A prototype system is developed in which geospatial linked data are used as underlying data, and some geovisualisation knowledge is formalised into a knowledge base to visualise the data and provide rich semantics to users. The proposed approach can partially form the foundation for the vision of web of knowledge for geovisualisation.