Integration of local–global upscaling and grid adaptivity for simulation of subsurface flow in heterogeneous formations

We propose a methodology, called multilevel local–global (MLLG) upscaling, for generating accurate upscaled models of permeabilities or transmissibilities for flow simulation on adapted grids in heterogeneous subsurface formations. The method generates an initial adapted grid based on the given fine-scale reservoir heterogeneity and potential flow paths. It then applies local–global (LG) upscaling for permeability or transmissibility [7], along with adaptivity, in an iterative manner. In each iteration of MLLG, the grid can be adapted where needed to reduce flow solver and upscaling errors. The adaptivity is controlled with a flow-based indicator. The iterative process is continued until consistency between the global solve on the adapted grid and the local solves is obtained. While each application of LG upscaling is also an iterative process, this inner iteration generally takes only one or two iterations to converge. Furthermore, the number of outer iterations is bounded above, and hence, the computational costs of this approach are low. We design a new flow-based weighting of transmissibility values in LG upscaling that significantly improves the accuracy of LG and MLLG over traditional local transmissibility calculations. For highly heterogeneous (e.g., channelized) systems, the integration of grid adaptivity and LG upscaling is shown to consistently provide more accurate coarse-scale models for global flow, relative to reference fine-scale results, than do existing upscaling techniques applied to uniform grids of similar densities. Another attractive property of the integration of upscaling and adaptivity is that process dependency is strongly reduced, that is, the approach computes accurate global flow results also for flows driven by boundary conditions different from the generic boundary conditions used to compute the upscaled parameters. The method is demonstrated on Cartesian cell-based anisotropic refinement (CCAR) grids, but it can be applied to other adaptation strategies for structured grids and extended to unstructured grids.     

Accurate local upscaling with variable compact multipoint transmissibility calculations

We propose a new single-phase local upscaling method that uses spatially varying multipoint transmissibility calculations. The method is demonstrated on two-dimensional Cartesian and adaptive Cartesian grids. For each cell face in the coarse upscaled grid, we create a local fine grid region surrounding the face on which we solve two generic local flow problems. The multipoint stencils used to calculate the fluxes across coarse grid cell faces involve the six neighboring pressure values. They are required to honor the two generic flow problems. The remaining degrees of freedom are used to maximize compactness and to ensure that the flux approximation is as close as possible to being two-point. The resulting multipoint flux approximations are spatially varying (a subset of the six neighbors is adaptively chosen) and reduce to two-point expressions in cases without full-tensor anisotropy. Numerical tests show that the method significantly improves upscaling accuracy as compared to commonly used local methods and also compares favorably with a local–global upscaling method.     

一种测求水井含水层导水系数的新方法

本文给出了一种测求水井含水层导水系数的新方法。利用Cooper理论和振动理论,通过简单的试验,可以测求出水井含水层的导水系数。用这种方法计算出珍珠泉井含水层的导水系数为2439m~2/d,用抽水试验法测得该系数为2618m~2/d,两者符合得较好。     

井水位振荡试验及其结果

通过鲁04井的频率特性试验,求出了该井含水层系统的频率特性参数,固有振动周期为40.59s,阻尼系数为0.0362s-1。并且计算出鲁04井含水层的导水系数为1.8×10-3m2/s。利用这些结果解释了鲁04井对地震波响应比珍珠泉井差的原因     

潜水含水层导水系数空间分布特征的初步分析

以新疆三工河流域浅层含水层为研究对象,采用地质统计学中的原理与方法,初步分析了该含水层导水系数的空间分布特征.结果表明,1)该含水层导水系数在空间上的变化服从正态概率分布;2)含水层的导水性能在空间上属于中等变异性;3)导水系数在空间位置的一定范围内具有相关性,东西方向的相关尺度a=5.0km,南北方向的相关尺度a=7.8km;4)该含水层的导水性能为各向异性的.文中的研究方法与步骤可以用来分析其它冲洪积浅层含水层导水系数的空间分布特征.     

中国食管癌与土壤环境中元素的相关性

利用土壤元素资料 2 0 6 94 4 8个数据 ,食管癌死亡调查资料 1750 80例 ,研究了食管癌死亡率与人群生存的土壤环境中元素的相关性。结果表明 ,铅、硒、铷、铯、镁、钙、镓、铊、钍、铀、锡、铪、铋、钼、溴、碘元素与食管癌死亡率有相关性     

Nonlinear two-point flux approximation for modeling full-tensor effects in subsurface flow simulations

Subsurface flow models can exhibit strong full-tensor anisotropy due to either permeability or grid nonorthogonality effects. Upscaling procedures, for example, generate full-tensor effects on the coarse scale even for cases in which the underlying fine-scale permeability is isotropic. A multipoint flux approximation (MPFA) is often needed to accurately simulate flow for such systems. In this paper, we present and apply a different approach, nonlinear two-point flux approximation (NTPFA), for modeling systems with full-tensor effects. In NTPFA, transmissibility (which provides interblock connections) is determined from reference global flux and pressure fields for a specific flow problem. These fields can be generated using either fully resolved or approximate global simulations. The use of fully resolved simulations leads to an NTPFA method that corresponds to global upscaling procedures, while the use of approximate simulations gives a method corresponding to recently developed local–global techniques. For both approaches, NTPFA algorithms applicable to both single-scale full-tensor permeability systems and two-scale systems are described. A unified framework is introduced, which enables single-scale and two-scale problems to be viewed in a consistent manner. Extensive numerical results demonstrate that the global and local–global NTPFA techniques provide accurate flow predictions over wide parameter ranges for both single-scale and two-scale systems, though the global procedure is more accurate overall. The applicability of NTPFA to the simulation of two-phase flow in upscaled models is also demonstrated.     

水井含水层导水系数及其对地震波的响应

利用弹性理论和渗流理论研究了井水位振荡试验,给出了一种计算水井含水层渗流特性参数的新方法。该方法不是用井水位的恢复曲线,而是用井水位的振荡曲线计算出水井含水层的导水系数。对５ 口水井的导水系数及其对地震波的响应进行了统计,结果表明：水井含水层的导水系数越大,井水位对地震波响应的幅度越大,响应的次数越多。     

测试水井导水性能的简便试验

这种测试水井含水层导水系数的新方法是利用井口水面空气压力阶变试验,计算出聊城井含水层的导水系数为14.4m2/d,而通过抽水试验得出的导水系数为43.6m2/d,二者有相同的数量级,但新方法比抽水试验简单易行和经济快速,用它来解释地下水微动态更合理,容易推广应用