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.     

Shear waves at a corrugated interface between anisotropic elastic and visco-elastic solid half-spaces

Rayleigh’s method of approximation is employed to find out the reflection and transmission coefficients due to an incident plane SH wave at a corrugated interface between a laterally and vertically inhomogeneous anisotropic elastic solid half-space and a laterally and vertically inhomogeneous isotropic visco-elastic solid half-space. The lateral and vertical inhomogeneities are described by the exponential variations of elastic parameters. The formulae of reflection and transmission coefficients are derived in closed form for the first-order approximation of the corrugation. The effects of the corrugation of the interface, the inhomogeneity, the anisotropy, the visco-elasticity and the frequency of the incident wave on these coefficients are studied analytically and numerically for a specific model containing a periodic interface. The results of earlier workers have been reduced as particular cases from the present formulation.     

裂隙岩体地下水溶质运移的尺度问题

核废料的安全贮存、地下污染的预测控制等目前面临的许多重大问题要求对大尺度裂隙岩体地下水溶质运移问题进行研究,然而裂隙岩体地下水溶质运移问题的尺度效应给这项研究带来了极大的挑战.本文介绍了裂隙岩体地下水溶质运移的多尺度概念模型及其物理和数学模拟方法,并探讨了模拟尺度与观测尺度、预测尺度、裂隙介质尺度之间的关系及尺度效应的分形特征.     

横向各向同性介质中的弹性波方程逆时偏移及其成像条件

作者在本文中导出了在横向各向同性介质情况下，二维弹性波逆时传播的高阶差分格式，实现了弹性波在数值空间中的逆时延拓，从snell定量出发，采用插值法求取网格空间中各点的直达波旅行时，以此作为弹性波逆时偏移的成像条件，实现多波多分量资料的逆时偏移，数值试验表明，逆时偏移能够使地表接收到的波场准确归位，提高资料处理的精度。     

Scattering of seismic waves by cracks in multi-layered geological regions: II. Numerical results

The mechanical model for plane strain, time-harmonic seismic wave propagation problems in cracked, multi-layered geological regions with surface topography and non-parallel interfaces was described in the first part of this work. Here, this model is used to investigate the response of such a region to the presence of traveling elastic waves generated by a seismic source. The computational methodology that was developed in the first part is based on a combination of both the regular (displacement-based) and the hypersingular (traction-based) Boundary Integral Equation Method (BIEM). First, the accuracy and convergence characteristics of this hybrid BIEM are studied. Then, a series of problems involving four different configurations of a reference geological deposit with both interface and internal cracks are solved, for a loading that is due to a seismically-induced pressure wave propagating upwards from the underlying rigid half-plane. The purpose of the numerical study is to investigate the influence of various key parameters of the problem, such as frequency and incidence angle of the incoming wave, size of the surface relief, location and size of the buried cracks, interaction effects between cracks and finally the presence of layers, on both the scattered displacement field and the stress concentration field.     

Rocking response of flexible circular foundations on layered media

The analysis of the response of a flexible circular foundation on layered media due to an arbitrarily distributed vertical loading is presented. The analysis is based on the ‘ring method’ approach, i.e. discretization of the foundation in a set of concentric rings. The arbitrarily distributed loading is expanded in the circumferential direction in a Fourier series. The influence coefficient matrix of soil for each element of the series is evaluated utilizing the stiffness matrix approach. The stiffness matrix of the foundation is obtained from the finite difference energy method approach. Numerical examples illustrate the influence of several soil-foundation parameters on the rocking response of a foundation. Results are presented in terms of displacement and soil reaction distributions and impedance functions point to significantly different responses of flexible and rigid foundations.     

对称化测量高频电磁波测井的优化组合研究

采用对称化测量可克服现有单发双收高频电磁波测井仪器（ＤＰＴ）由于收－发不对称结构导致的幅度比测井曲线严重失真的缺陷，提高ＤＰＴ的纵向分辨率．通过对二线非均质测井环境下ＤＰＴ响应的数值模拟，指出选用原频率较低、源距较长仪器的相位差测量与源频率较高、源距较短仪器的幅度比测量进行组合，可获得高品质的视介电常数和视电阻率曲线，显著提高ＤＰＴ的资料解释能力．这里源频率的恰当选择起着关键作用．