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.     

鄂尔多斯盆地靖边气田马五_(1+2)气藏储集单元研究

综合地震、地质、测井等信息,考虑沉积、成岩、地貌、生产动态等多因素,在沟槽识别和致密带研究基础上,通过井间干扰、油层压力变化和试井边界响应分析,按照局部、就近原则对靖边气田马五1+2气藏进行了储集单元划分。共划分出111个储集单元,其中独立型69个,复合型42个。这些储集单元的面积多介于10～50km2之间,单个储集单元的井数一般少于6口。在对气藏储层特征综合认识的基础上,确定流动单元划分思路与流程,并认为关键环节包括古沟槽识别、致密带研究、动静结合确定储集单元等3个方面。储集单元一方面受控于被致密岩相包围的有利成岩相带,宏观上与弱—中等充填泥粉晶白云岩相对应性好;另一方面又受各级古沟槽夹持,在两个一级古沟槽之间往往共生多个储集单元,而这些储集单元边界又严格地受二级或三级沟槽切割控制。     

Laboratory simulation of geogas microbubble flow

Preliminary laboratory tests provided first data on the behavior of gas microbubbles through porous media in the framework of the geogas theory. Under experimented conditions with laboratory equipment arranged for pressure controlled gas-tracer injection and sampling, gas microbubbles moved up to ten times faster than singlephase flow in dry media under the same injection pressure. Microbubbles were determined to be very sensitive to changes in injection pressure and their terminal velocity seems to be described with good approximation by the Stokes formula. The capability of microbubbles to pick up and transport upward for short distances solid ultra-small particles (metallic and radionuclide compounds) has been proved. Results are consistent with a time-dependent process linked to the transport properties of microbubbles (e.g., flotation), such as that reported by some authors.     

Phase transitions in the ISM a source of dissipative behaviour

By considering a simple fluid model, we investigate the role of phase transitions in the ISM on the galaxy- scale gas dynamics. Cooling and heating timescales in the ISM are typically shorter than typical galactic rotation timescales, so the individual phases in the ISM can be assumed to be in temperature equilibrium with the radiation field. Using this approximation we can construct an equation of state which depends upon the average density and mass fractions in the individual phases. Previous studies suggest that there is an equilibrium phase fraction as a function of pressure. We incorporate evolution towards this equilibrium state as a relaxation term with a time to obtain equilibrium . We derive a condition in terms of a critical Mach number when one dimensional shocks should be continuous. For small values of the relaxation time we show that the relaxation term acts like a viscosity. We show with one dimensional simulations that increasing causes shocks to become smoother. This work suggests that phase changes can strongly effect the gas dynamics of the ISM across spiral arms and bars.     

High-resolution shallow seismic identification of gas escape features in the sediments of Loch Tay, Scotland: tectonic and microbiological associations

High-resolution seismic reflection profiles (3·5 kHz) have revealed the presence of extensive interstitial gas accumulation within the sedimentary sequences of Loch Tay, Scotland, as identified by acoustic turbidity masking the seismic stratigraphy. Within the central section of the loch, in the deepest water area directly above the zone of the seismically active Loch Tay Fault, focused flows of gas through the sediment pile to the loch bed via chimneys and pockmarks, together with gas seeps within the water column, have been identified. Microbiological observations indicate that the gas is biogenic CH₄, produced by both chemoautotrophic (which use CO₂ as a source of carbon and H₂ as a source of energy) and aceticlastic species (which use acetate as a source of carbon and energy) of methanogens in the fine-grained, organic rich deposits that have been focused into the zone of accumulation in the deep central part of the loch. The spatial distribution of the gas escape features suggests that earthquake movements along the Loch Tay Fault are responsible for facilitating focused gas escape in this part of the loch, by the creation of new pathways and conduits through the sediment pile, along which gas can migrate upwards and exit into the water column. Relict pockmarks and associated chimneys identified in the seismic records indicate that gas escape has been taking place since Pleistocene times though the precise timings cannot be ascertained. This is the first time that such features have been reported from a lake in the UK.     

Rock fracturing by explosive energy: review of state-of-the-art

A study of the dynamic rock fracture initiation and propagation due to explosive energy is presented through a detailed state-of-the-art review. Explosive energy dissipation in crushing and fracturing is examined and the various means to enhance the explosive energy utilization for dynamic rock fracturing are reviewed. The study highlights the need for a better understanding of the dynamic fracturing process particularly in the presence of in situ stresses in the rock mass.     

冲绳海槽玄武岩中中酸性残余熔体研究及其岩石学意义

与洋壳有关的酸性岩由于对了解幔源岩浆的演化以及判别古老蛇绿岩套及其构造位置的重要意义而倍受岩石学家和构造地质学家的关注。本对发现于冲绳海槽玄武岩基质中的中酸性到酸性残余熔体进行了详细研究，它们提供了幔源玄武岩浆结晶分异形成酸性岩浆的直接证据。在细小的基质矿物间分布有一种玻璃质的残余熔体，其成分随距冷凝边距离(L)的增加而越来越酸性。在SiO2对Na2O K2O图解上，残余熔体的投影点从玄武岩到英安岩均有分布，反映了一个连续的演化系列。在AFM图解上，残余熔体表现出与Thingmuli火山岩系列类似的拉斑玄武岩系列的演化趋势。我们的研究表明：残余熔体的演化受结晶分异作用控制。在早期结晶阶段，辉石的结晶起主导作用，结果造成残余熔体中SiO2、Al2O3，Na2O的迅速增加，FeO、MgO、CaO迅速降低。在晚期结晶阶段，斜长石成为主导结晶相，导致残余熔体中Al2O3，Na2O的迅速消减。Al2O3、Na2O从增加到降低的转变出现在SiO2=62％左右。在L=27．5mm处，85～90％的基质岩浆已发生了结晶作用，导致残余熔体中SiO2含量达到69～70％，而且此处还新出现了一种富FeTi的氧化物。该玄武岩中残余熔体和基质矿物的成分及演化特征分别与Thingmuli火山岩系列中酸性端元的组成相似，在也佐证了Thingmuli火山岩系列是幔源岩浆结晶分异的产物.     

金刚石工具富铁胎体掺杂稀土的研究

稀土的加入量、加入形态和在混料中的均匀弥散性直接影响热压富铁金刚石复合材料的性能。改进的稀土掺杂工艺，保证了稀土在胎体中的均匀弥散性；通过试验研究了稀土的加入量与富铁胎体的抗弯强度、抗冲击韧性和孔隙率的关系，从而确定了稀土的最优加入量。通过差热分析试验，认为稀土可以改变富铁胎体的热物理特性。     

滇西北地区富碱岩体(脉)地质学及岩石化学特征

喜山期富碱岩体广布于滇西北地区 ,通过对其地质学和岩石化学特征的研究 ,认为富碱岩体 (脉 )由正长 (斑 )岩、二长斑岩、花岗斑岩及相关的火山岩组成 ,具有多次活动的特征 ,并形成区内富有特色的岩浆岩系列 ;该系列岩浆岩以富碱高钾为突出特点 ,具有正常的岩浆演化趋势 ,不同岩带岩石地质学和岩石化学方面有区域性变化规律。     

黄铁矿中富铜闪锌矿的斑点状结构

内蒙古大井矿床中产出一种奇特的闪锌矿斑点结构，闪锌矿斑点包裹于黄铁矿中，通过显微镜研究和电子探针分析，认为该闪锌矿斑点属罕见的富铜闪锌矿，ω(Cu)为3．87％—6．99％，平均5．37％，其成因与在黄铜矿中的闪锌矿斑点一样，是高温下出溶作用的产物。