Study on electric variations of media in epicentral area by geomagnetic transfer functions

ＳｔｕｄｙｏｎｅｌｅｃｔｒｉｃｖａｒｉａｔｉｏｎｓｏｆｍｅｄｉａｉｎｅｐｉｃｅｎｔｒａｌａｒｅａｂｙｇｅｏｍａｇｎｅｔｉｃｒａｎｓｆｅｒｆｕｎｃｔｉｏｎｓＸｉａｏ－ＰｉｎｇＺＥＮＧ；（曾小苹）Ｙｕｎ－ＦａｎｇＬＩＮＩ；（林云芳）Ｚｈｏｎｇ－ＪｉｅＺＨ...     

横向非均匀介质中弹性地震波传播的数值模拟

应用混合变量弹性动力学方程和线性常微分方程组的矩阵指数解法，将层状介质中广泛应用的弹性波传播矩阵解法推广至横向非均匀介质，给出了一种可计算复杂地质体中弹性波传播的广义传播矩阵数值解法。该方法可模拟任意震源及所产生的各种体波、面波，数值结果表明具有很高的计算精度。     

Micro-scale modelling in the study of plume evolution in heterogeneous media

The migration of contaminants in heterogeneous aquifers involves dispersive processes that act at different scales. The interaction of these processes as a plume evolves can be studied by micro-scale modelling whereby two scales, a local- or micro-scale and an aquifer- or macro-scale, are covered simultaneously. Local-scale dispersive processes are represented through the local dispersion coefficient in the transport equation, while large-scale dispersion due to heterogeneities is represented through the resolution of the flow field and the diffusive exchange between streamtubes. The micro-scale model provides both the high degree of resolution compatible with local-scale processes, and the extent required for the approach to asymptotic conditions, using grids of up to a million nodal points. The model is based on the dual potential-streamfunction formulation for flow, and the transport problem is formulated in a natural coordinate system provided by the flownet. Simulations can be used to verify stochastic theories of dispersion, without the restrictive assumptions inherent in the theory. For the two-dimensional case, results indicate convergence of the effective dispersivity to the theoretical macrodispersivity value. Convergence takes place within a travel distance of about 50 correlation lengths of the hydraulic conductivity field. However, the approach taken to asymptotic conditions, as well as the macrodispersivity value, may differ for different realizations of the same medium. The influence of early-time events such as plume splitting on the asymptotic convergence remains to be investigated.     

Risk Management for Petroleum Reservoir Production: A Simulation-Based Study of Prediction

We consider numerical solutions of the Darcy and Buckley–Leverett equations for flow in porous media. These solutions depend on a realization of a random field that describes the reservoir permeability. The main content of this paper is to formulate and analyze a probability model for the numerical coarse grid solution error. We explore the extent to which the coarse grid oil production rate is sufficient to predict future oil production rates. We find that very early oil production data is sufficient to reduce the prediction error in oil production by about 30%, relative to the prior probability prediction.     

考虑土体固液二相性质的地下管线地震反应研究

本文将管线周围的土体视为固-液二相介质，采用考虑地震过程中土体内孔隙水压力的增长和消散的有效应力分析方法，并以非线性的土体本构模型模拟土体动力特性的改变以及在孔隙水压作用下的变形，以求接近土体的真实地震反应。研究表明随着土体软化程度的不同，受到土体约束的管线应力亦发生相应变化。文中详细研究了地震动峰值、管线直径、行波波速和入射角度等对管线应力的影响。     

Density-dependent dispersion in heterogeneous porous media Part I: A numerical study

In this paper, we describe carefully conducted numerical experiments, in which a dense salt solution vertically displaces fresh water in a stable manner. The two-dimensional porous media are weakly heterogeneous at a small scale. The purpose of these simulations, conducted for a range of density differences, is to obtain accurate concentration profiles that can be used to validate nonlinear models for high-concentration-gradient dispersion. In this part we focus on convergence of the computations, in numerical and statistical sense, to ensure that the uncertainty in the results is small enough.Concentration variances are computed, which give estimates of the uncertainty in local concentration values. These local variations decrease with increasing density contrast. For tracer transport, obtained longitudinal dispersivities are in accordance with analytical findings. In the case of high-density contrasts, stabilizing gravity forces counteract the growth of dispersive fingers, decreasing the effective width of the transition zone. For small log-permeability variances, the decrease of the apparent dispersivity that is found is in agreement with laboratory results for homogeneous columns.     

3D image-based characterization of fluid displacement in a Berea core

Improved network flow models require the incorporation of increasingly accurate geometrical characterization of the microscale pore structure as well as greater information on fluid–fluid interaction (interfaces) at pore scales. We report on three dimensional (3D) pore scale medium characterization, absolute permeability computations for throat structures, and pore scale residual fluid distribution in a Berea core. X-ray computed microtomography combined with X-ray attenuating dopants is used to obtain 3D images of the pore network and to resolve phase distributions in the pore space.     

Inversion of interference hydraulic pumping tests in both homogeneous and fractal dual media

This work proposes a complete method for automatic inversion of data from hydraulic interference pumping tests based on both homogeneous and fractal dual-medium approaches. The aim is to seek a new alternative concept able to interpret field data, identify macroscopic hydraulic parameters and therefore enhance the understanding of flow in porous fractured reservoirs. Because of its much contrasted sensitivities to parameters, the dual-medium approach yields an ill-posed inverse problem that requires a specific optimization procedure including the calculation of analytical sensitivities and their possible re-scaling. Once these constraints are fulfilled, the inversion proves accurate, provides unambiguous and reliable results. In the fractal context inverting several drawdown curves from different locations at the same time reveals more accurate. Finally, hydraulic parameters drawn from inversion should be taken into account to improve in various situations the conditioning of up-scaled flow in fractured rocks.     

Coarsening of three-dimensional structured and unstructured grids for subsurface flow

We present a generic, semi-automated algorithm for generating non-uniform coarse grids for modeling subsurface flow. The method is applicable to arbitrary grids and does not impose smoothness constraints on the coarse grid. One therefore avoids conventional smoothing procedures that are commonly used to ensure that the grids obtained with standard coarsening procedures are not too rough. The coarsening algorithm is very simple and essentially involves only two parameters that specify the level of coarsening. Consequently the algorithm allows the user to specify the simulation grid dynamically to fit available computer resources, and, e.g., use the original geomodel as input for flow simulations. This is of great importance since coarse grid-generation is normally the most time-consuming part of an upscaling phase, and therefore the main obstacle that has prevented simulation workflows with user-defined resolution. We apply the coarsening algorithm to a series of two-phase flow problems on both structured (Cartesian) and unstructured grids. The numerical results demonstrate that one consistently obtains significantly more accurate results using the proposed non-uniform coarsening strategy than with corresponding uniform coarse grids with roughly the same number of cells.