Multidimensional upstream weighting for multiphase transport on general grids

The governing equations for multiphase flow in porous media have a mixed character, with both nearly elliptic and nearly hyperbolic variables. The flux for each phase can be decomposed into two parts: (1) a geometry- and rock-dependent term that resembles a single-phase flux; and (2) a mobility term representing fluid properties and rock–fluid interactions. The first term is commonly discretized by two- or multipoint flux approximations (TPFA and MPFA, respectively). The mobility is usually treated with single-point upstream weighting (SPU), also known as dimensional or donor cell upstream weighting. It is well known that when simulating processes with adverse mobility ratios, SPU suffers from grid orientation effects. An important example of this, which will be considered in this work, is the displacement of a heavy oil by water. For these adverse mobility ratio flows, the governing equations are unstable at the modeling scale, rendering a challenging numerical problem. These challenges must be addressed in order to avoid systematic biasing of simulation results. In this work, we present a framework for multidimensional upstream weighting for multiphase flow with buoyancy on general two-dimensional grids. The methodology is based on a dual grid, and the resulting transport methods are provably monotone. The multidimensional transport methods are coupled with MPFA methods to solve the pressure equation. Both explicit and fully implicit approaches are considered for time integration of the transport equations. The results show considerable reduction of grid orientation effects compared to SPU, and the explicit multidimensional approach allows larger time steps. For the implicit method, the total number of non-linear iterations is also reduced when multidimensional upstream weighting is used.     

Errors in the upstream mobility scheme for countercurrent two-phase flow in heterogeneous porous media

In reservoir simulation, the upstream mobility scheme is widely used for calculating fluid flow in porous media and has been shown feasible for flow when the porous medium is homogeneous. In the case of flow in heterogeneous porous media, the scheme has earlier been shown to give erroneous solutions in approximating pure gravity segregation. Here, we show that the scheme may exhibit larger errors when approximating flow in heterogeneous media for flux functions involving both advection and gravity segregation components. Errors have only been found in the case of countercurrent flow. The physically correct solution is approximated by an extension of the Godunov and Engquist–Osher flux. We also present a new finite volume scheme based on the local Lax–Friedrichs flux and test the performance of this scheme in the numerical experiments.     

On several numerical strategies to solve Richards’ equation in heterogeneous media with finite volumes

Computational Geosciences - We propose several numerical approaches building on upstream mobility two-point flux approximation finite volumes to solve Richards’ equation in domains made of...     

Linear stability analysis of the explicit treatment of mobilities in non-Newtonian and non-Darcy porous media flow simulation

A von Neumann stability analysis of the discretized conservation equation for single-phase porous media flows is performed, where non-Newtonian and non-Darcy effects are accounted for using a velocity (or mass flux)-dependent mobility factor. Comprehensive results in three dimensions for two low-order finite-volume discretizations typically encountered in reservoir simulation are provided, based on edge-centered and upstream cell-centered mobility calculations. It is found that common semi-implicit schemes, where the pressure gradient driving the flow is taken implicitly while the velocity-dependent mobility is evaluated explicitly, are subject to restrictions on the logarithmic derivative of mobility with respect to velocity. A remarkable new result is nevertheless obtained: for any physically acceptable strength of non-Newtonian and non-Darcy effects, there exists a stable and explicit method to evaluate the mobility, rendering the need to implement costly fully implicit schemes more difficult to justify.     

Circulation and mixing due to tidal forcing in the Bertioga Channel, São Paulo, Brazil

Bertioga Channel is a partially mixed (type 2) tidal estuary on the coastal plain of São Paulo, Brazil. Hourly current and salinity measurements during neap and spring tides in July 1991 yielded information about the physical structure of the system. Peak along-channel velocities varied from 40 cm s^?1 to 60 cm s^?1 during flood tides and from 70 cm s^?1 to 100 cm s^?1 during ebb tides. Net vertical velocity profiles indicate that the net current reverses directions at a depth of 2.5–3.0 m in the halocline. Due to appreciable fortnightly tidal modulation, the estuary alternates from being highly stratified (type 2b) during neap tides, with advection and diffusion contributing equally to the net upstream salt flux, to being moderately stratified (type 2a) during spring tides, when 90% of the net upstream salt transport is the result of effective tidal diffusion. Decomposition of the salt flux indicates that the relative contribution to the upstream salt transport by gravitational circulation shear is greater than the oscillatory tidal flux by a factor of 2.6 during neap tides. The oscillatory tidal flux is generated by the correlation of the tidal components of the u-velocity and salinity and is responsible for approximately the same amount of upstream salt transport, during neap and spring tides. However, during spring tides, this oscillatory term is greater than the other salt flux terms by a factor of 1.4. The total salt transport, through a unit width of the section perpendicular to the flow, was within 2% of the sum of the seven major decomposed, advective and dispersive terms. On the assumption that the Bertioga Channel is laterally homogeneous, the results also indicate that the estuary is not in steady state with respect to salt flux.     

秦岭泥盆系硅岩岩石类型对比研究

硅岩可以通过化学沉积（真溶液或胶体沉淀）和内碎屑堆积而成。由于热水有较强的化学活动性，因而在水岩界面附近也会形成交代成因的硅岩。在许多情况下，硅岩既具有沉积岩的岩石学特征，又具有交代岩的岩石学特征，其成岩环境十分特殊。因此，虽然有许多地球化学方法可以区分不同类型的硅岩，但岩石学特征是最基本的和最直接的标志。     

渤海湾盆地古近纪沙三、四段沉积期古盐度不均衡性研究

渤海湾盆地不同凹陷区钻井岩芯资料显示,古近纪沙三、沙四段沉积的泥岩、灰岩、云岩、膏岩多呈现相互间隔的韵律层分布。但在不同凹陷及各洼陷内,其发育的厚度、缺失程度以及韵律的完整性却是不尽一致的。通过对孢粉资料的分析发现,当时的古气候在全盆地基本趋于稳定,干、湿程度的微弱变化不足以引发沉积剖面上如此大的岩性差异。对比下辽河、东营、沾化和东濮凹陷的沉积记录和古盐度数据,发现各地岩性组合的不一致性主要是由于各地湖盆古盐度差异所致。综合各项数据分析推断,沉积时期盆地东北部隆-堑夹持区下辽河凹陷(上游)的动态淡水不断向盆内注入,最终在西南部低洼区东濮凹陷(下游)汇集,造成了盆内咸化环境产生了较大的不均衡性。上游区咸化环境被淡水稀释,盐岩类型单一、厚度小;下游区则由于较稳定的高盐度环境,在遇干旱气候时,盐、膏层大量析出,而位处中游的东营、沾化凹陷则呈现由淡趋咸的过渡状态。高盐度咸化环境有利于有机质的保存,能够促成烃类的早生早排,对规模油气资源量的形成极为有利。     

采场结构体回采可移动性分析

采场范围内存在的结构体严重影响地下矿山安全生产,分析不同回采顺序下采场结构体可移动性具有重要意义。以块体理论为基础,结合空场法开采的采场空间属性,分析了结构体在采场内的一般嵌套模式。运用GeneralBlock程序建立了采场结构体模型,对采场内结构体的类型、体积和位置等元素进行了解构。利用矢量方法分析了不同回采顺序下采场结构体可移动性,并利用力学方法分析了考虑岩桥破坏时的矢量不可动结构体的回采移动性,推导了结构体破坏岩桥移动的条件公式。以铜坑矿锌多金属矿体1#矿块采场为样本,进行了回采可移动性实例分析,结果表明：一端至另一端后退回采时,4个结构体滑动;中间至两端后退回采时,7个结构体滑动,1个结构体破坏岩桥后滑动;上向水平分层回采时,4个结构体滑动,2个结构体掉落。研究成果可为地下矿山裂隙矿岩空场法回采顺序的合理选择及采场安全管理提供技术指导。     

Landslide at Su-Hua Highway 115.9k triggered by Typhoon Megi in Taiwan

This study focused on the landslide case at Su-Hua Highway 115.9k, Taiwan. A preliminary investigation was conducted on geomorphologic features change and landslide mechanisms using digital elevation models, geographical maps, and remote sensing images at different times in conjunction with geological surveys and analysis results. Using the results of geological surveys and physical model experiments, we constructed a discrete element method to simulate the process of landslide movement. The results revealed deformation in the metamorphic rock slopes upstream of 115.9k. The slopes around the erosion gully upstream presented visible slope toes cutting and tension cracks at the crest as well as unstable rock masses. According to the results of numerical simulation for typhoon Megi event, intense rains could induce slippage in the rock debris/masses in the source area, initially at a speed of 5–20 m/s. Subsequently, steeper terrain could cause the rock debris/masses to accelerate to form a high-speed (>30 m/s) debris slide quickly moving downstream to form an alluvial fan downstream by the sea.     

一种高边坡岩体卸荷分带方法的探讨

岩体卸荷是水利水电工程建设中常见的工程地质问题之一。本文介绍了高边坡岩体卸荷的背景条件、作用机理、卸荷岩体特征以及常见的卸荷带划分方法。针对河谷斜坡岩体从表部到深部 ,岩体裂隙发育程度、张开程度逐渐降低的特征 ,以黄河上游某水电站高边坡印支期花岗岩体的现场调查资料为依托 ,应用岩体裂隙发育程度、张开程度量化指标对斜坡岩体卸荷进行分带 ,较好地将岩体卸荷分带与岩体结构、岩体工程特性紧密结合起来 ,为以后其它边坡岩体的工程实践提供了借鉴.     

Cosmogenic nuclide methods for measuring long-term rates of physical erosion and chemical weathering

Understanding the evolution of geochemical and geomorphic systems requires measurements of long-term rates of physical erosion and chemical weathering. Erosion and weathering rates have traditionally been estimated from measurements of sediment and solute fluxes in streams. However, modern sediment and solute fluxes are often decoupled from long-term rates of erosion and weathering, due to storage or re-mobilization of sediment and solutes upstream from the sampling point. Recently, cosmogenic nuclides such as ¹⁰Be and ²⁶Al have become important new tools for measuring long-term rates of physical erosion and chemical weathering. Cosmogenic nuclides can be used to infer the total denudation flux (the sum of the rates of physical erosion and chemical weathering) in actively eroding terrain. Here we review recent work showing how this total denudation flux can be partitioned into its physical and chemical components, using the enrichment of insoluble tracers (such as Zr) in regolith relative to parent rock. By combining cosmogenic nuclide measurements with the bulk elemental composition of rock and soil, geochemists can measure rates of physical erosion and chemical weathering over 1000- to 10,000-year time scales.     

近距离上保护层开采瓦斯运移规律数值分析

采动裂隙是瓦斯运移的通道,搞清瓦斯运移规律是瓦斯治理的前提。在考虑岩石动态破坏过程和含瓦斯煤岩渗流-应力-损伤耦合的基础上,结合平煤五矿实际地质条件和开采工艺,建立了数值计算模型,应用RFPA-Gas程序模拟了近距离上保护层采动顶底板岩层变形破坏、裂隙演化规律与瓦斯运移规律。模拟结果较好地再现了保护层开采过程中煤岩层应力变化、顶底板损伤及裂隙演化过程,得到了上覆岩层移动的“上三带”（冒落带、裂隙带和弯曲下沉带）和底板变形的“下两带”（底板变形破坏带和弹塑性变形带）。得到了被保护层瓦斯流量分布、瓦斯压力分布和透气系数的变化规律,卸压煤层瓦斯透气性增大了2 500倍,得到了煤壁下方压缩区和膨胀区之间的张剪瓦斯渗流通道,并将保护层底板压缩区和膨胀区的瓦斯渗流特征提炼出来：压缩区对应的是渗流减速减量区、膨胀区由卸压膨胀陡变区和卸压膨胀平稳区组成,分别对应着渗流急剧增速增量区和渗流平稳增量区。指出卸压膨胀陡变区是瓦斯突出危险区,为近距离保护层开采瓦斯治理指明了方向。实践表明,瓦斯治理效果显著。     

A phase-by-phase upstream scheme that converges to the vanishing capillarity solution for countercurrent two-phase flow in two-rock media

We discuss the convergence of the upstream phase-by-phase scheme (or upstream mobility scheme) towards the vanishing capillarity solution for immiscible incompressible two-phase flows in porous media made of several rock types. Troubles in the convergence were recently pointed out by Mishra and Jaffré (Comput. Geosci. 14, 105–124, 2010) and Tveit and Aavatsmark (Comput. Geosci. 16, 809–825, 2012). In this paper, we clarify the notion of vanishing capillarity solution, stressing the fact that the physically relevant notion of solution differs from the one inferred from the results of Kaasschieter (Comput. Geosci. 3, 23–48, 1999). In particular, we point out that the vanishing capillarity solution depends on the formally neglected capillary pressure curves, as it was recently proven in by Andreianov and Cancès (Comput. Geosci. 17, 551–572, 2013). Then, we propose a numerical procedure based on the hybridization of the interfaces that converges towards the vanishing capillarity solution. Numerical illustrations are provided.     

湘东锡田钨锡矿区成岩成矿时代研究

湘东锡田钨锡矿是近年来在找矿上有重大突破的地区。笔者在分析前人资料的基础上,应用Rb-Sr全岩法和Re-Os辉钼矿法对锡田花岗岩及产于其中的云英岩石英脉型钨锡矿进行了精确定年。结果表明:锡田岩体主体、补体和晚期花岗岩形成年龄分别为:165±16Ma(2σ)、151±24Ma(2σ)、114±14Ma(2σ);早期云英岩石英脉型钨锡矿辉钼矿Re-Os等时线年龄为150±2.7Ma。研究显示该区可能还残留有印支期花岗岩,钨锡成矿与燕山期花岗岩侵入活动有关。不同类型钨锡矿生成顺序为:矽卡岩型→早期云英岩石英脉型→晚期云英岩石英脉型→破碎带蚀变岩型。     

韧性剪切带内流体作用及成分变化与体积变化

韧性剪切带中的大量流体活动,影响了变形机制和矿物成分的变化。由于流体的渗透,造成某些元素的迁移,形成剪切带的物质成分变化和体积变化。流体与岩石的相互作用则是剪切带发展过程中的重要作用。应用成分变化计算体积变化已成为研究的热点问题     

铲子坪金矿床金的赋存状态及分布规律

铲子坪金矿是产于震旦系中的破碎带蚀变岩型金矿床,金主地以独立的金矿物形式存在,金的赋存状态及分布与硫化物的发育程度密切相关。金优先赋存于硫化矿物中,由于各矿带矿物组合及含量的变化,使金的赋存状态及分布也发生相应的改变。     

岩溶试片的岩性差异对估算岩溶速率和碳通量的影响

研究岩性差异对溶蚀速率的影响有助于提高溶蚀试片法估算岩溶碳汇强度的精确度。本文以贵州省普定县为研究区,将埋放地的主要基岩类型（石灰岩与白云岩）制成标准尺寸的试片,并将其埋设于不同土地利用类型和土壤深度下,经过4个水文年的监测后将估算的结果与前人在同一区域使用标准溶蚀试片的研究结果进行对比分析,结果表明:（1）在相同气候和土壤环境条件下,岩性对溶蚀试片的溶蚀速率有显著影响,且溶蚀速率与岩石中CaO含量呈正相关关系,与MgO含量呈负相关关系;（2）石灰岩与白云岩试片溶蚀速率的差异程度受土地利用及埋放深度的调控,整体上石灰岩溶蚀速率比白云岩溶蚀速率大14%;（3）不同岩性试片估算的岩溶碳汇强度相差较大,标准溶蚀试片估算的结果比埋放地基岩试片估算的结果高。故使用溶蚀试片法估算区域岩溶碳通量时应考虑埋放地基岩的岩性,或者对基于标准溶蚀试片的估算结果进行校正,才能准确反映区域尺度真实的岩溶碳通量大小。      

A new device for investigating the hydro‐mechanical properties of rock joints

A new experimental system composed of two devices is presented for studying the hydraulic conductivity of rock joints and other interfaces in relation with their mechanical loading. This system allows for investigations on the anisotropy of the hydraulic conductivity, as well as for its heterogeneity. The first device controls the injection of a fluid, regulated in pressure or discharge near the centre of the sample of rock joint, which produces a quasi‐radial flow within the joint. The output flow is measured in real time using a sectorized peripheral membrane for acquisition of data about the directionality or anisotropy of the flow. The second device allows for the measurement of the evolution of the morphology and of the void map of the rock joints during their mechanical loading (compression or shearing). It is a laser beam moved by an (x, y) frame. This frame is installed within the shear box at several stages of the loading, after temporarily moving one rock wall away from the other one, but without unmounting them from the shear box. After the presentation of the experimental system, and of its performances, we give some examples of test results of replicas of rock joints, loaded in compression and shearing. First, these results show that the change in global transmissivity of the samples is strongly correlated with the normal relative displacement of the rock walls. Several statistical analyses of the change in the void map, taking into account the damage of the asperities of the joint, are presented and compared with the results of the directional output flows. Today further studies involving a FEM as well as a finite volume modelling of the flow within the joints, based on the measured void map are in progress for a more precise comparison to the experimental results. Copyright © 2003 John Wiley & Sons, Ltd.     

江西相山矿田主要铀矿化类型及其地球化学特征对比研究

江西相山矿田发现近30个铀矿床,矿化与浅成超浅成侵入岩(次火山岩)关系密切,不管是划归火山岩型还是斑岩型矿床类型,均属于热液型矿床。从热液矿床的成矿构造特点分析,相山矿田铀矿化主要有隐爆碎屑岩型和蚀变岩型两大类,后者主要有水云母化蚀变岩型和钠长石化蚀变岩型。通过收集资料,结合笔者完成的有关课题成果数据,对各铀矿化类型进行了地球化学的对比研究,结果表明：隐爆碎屑岩型矿石品位较高,U含量大多大于1%,SiO₂、Na₂O含量较低,P₂O₅含量较高,K₂O/Na₂O比值平均11.51,主要伴生元素有Hf、Sb、Cu、Pb、Zn、Zr等;水云母化蚀变岩型矿化的K₂O/Na₂O比值平均4.83,伴生元素有W、Pb、Mo、Th、Sb等;钠长石化蚀变岩型矿化的K₂O/Na₂O比值平均0.19,U与CaO、P₂O₅关系密切,伴生元素有Sr、Zr、Hf、Sc、W等。从钠长石化蚀变岩型,到水云母化蚀变岩型,到隐爆碎屑岩型,即随矿化岩石U品位的增高,稀土配分曲线呈现出从右倾向左倾有规律的变化。