Coupled finite‐element simulation of injection well testing in unconsolidated oil sands reservoir

This paper presents a finite‐element (FE) model for simulating injection well testing in unconsolidated oil sands reservoir. In injection well testing, the bottom‐hole pressure (BHP) is monitored during the injection and shut‐in period. The flow characteristics of a reservoir can be determined from transient BHP data using conventional reservoir or well‐testing analysis. However, conventional reservoir or well‐testing analysis does not consider geomechanics coupling effects. This simplified assumption has limitations when applied to unconsolidated (uncemented) oil sands reservoirs because oil sands deform and dilate subjected to pressure variation. In addition, hydraulic fracturing may occur in unconsolidated oil sands when high water injection rate is used. This research is motivated in numerical modeling of injection well testing in unconsolidated oil sands reservoir considering the geomechanics coupling effects including hydraulic fracturing. To simulate the strong anisotropy in mechanical and hydraulic behaviour of unconsolidated oil sands induced by fluid injection in injection well testing, a nonlinear stress‐dependent poro‐elasto‐plastic constitutive model together with a strain‐induced anisotropic permeability model are formulated and implemented into a 3D FE simulator. The 3D FE model is used to history match the BHP response measured from an injection well in an oil sands reservoir. Copyright © 2013 John Wiley & Sons, Ltd.     

Hydraulic diffusivity and macrodispersivity calculations embedded in a geographic information system

Abstract

A numerical technique is presented whereby aquifer hydraulic diffusivities (D) and macrodispersivities (α) are calculated by linear equations rewritten from flow and solute transport differential equations. The approach requires a GIS to calculate spatial and temporal hydraulic head (h) and solute concentration gradients. The model is tested in Portugal, in a semi-confined aquifer periodically monitored for h and chloride/sulphate concentrations. Average D (0.46 m²/s) and α (1975 m) compare favourably with literature results. The relationship between α and scale (L) is also investigated. In this context, two aquifer groups could be identified: the first group is heterogeneous at the “macroscopic” scale (solute travelled distances <1 km), but homogeneous at the “megascopic” scale. The overall scale dependency in this case is given by an equation of logarithmic type. The second group is heterogeneous at the macroscopic and megascopic scales, with a scale dependency of linear type.

Citation Pacheco, F.A.L., 2013. Hydraulic diffusivity and macrodispersivity calculations embedded in a geographic information system. Hydrological Sciences Journal, 58 (4), 930–944.     

Modelling of groundwater flow in heterogeneous porous media by finite element method

The HySuf‐FEM code (Hydrodynamic of Subsurface Flow by Finite Element Method) is proposed in this article in order to estimate the spatial variability of the transmissivity values of the Berrechid aquifer (Morocco). The calibration of the model is based on the hydraulic head, hydraulic conductivity and recharge. Three numerical tests are used to validate the model and verify its convergence. The first test case consists in using the steady analytical solution of the Poisson equation. In the second, the model has been compared with the hydrogeological system which is characterized by an unconfined monolayer (isotropic layer) and computed by using PMWIN‐MODFLOW software. The third test case is based on the comparison between the results of HySuf‐FEM and the multiple cell balance method in the aquifer system with natural boundaries case. Good agreement between the Hydrodynamic of Subsurface Flow, the numerical tests and the spatial distribution of the thickening of the hydrogeological system is deduced from the analysis and the interpretations of hydrogeological wells. Copyright © 2011 John Wiley & Sons, Ltd.     

Geotechnical assessment of the bearing capacity of coal mine floors

Summary This investigation concerns 3 longwall faces having variable strata and mining conditions, but supported by 202 tonnes hydraulic 6-leg chocks. Bearing capacity tests along the faces were conducted to evaluate the factors influencing floor deformation and failure. The effect of size, shape and perimeter of the base plates, thickness of the floor layer, time and moisture on the ultimate bearing capacity of the floor was measured and discussed. Application of this work is the prediction of stability and support performance of face ends, as well as the design of support systems and ground control on production faces.     

北京市北山泥石流分布规律及形成过程

北京北山泥石流的分布受构造控制，呈北东方向展布。其分布在时间上和空间上均有一定的规律可寻。一条泥石流沟谷的形成从开始到结束，概括分为形成、位移和堆积。区内泥石流的形成主要有三种形式：以水力侵蚀为主，以坡面侵蚀为主和以沟源崩、滑塌形式触发沟床物质活动而形成的泥石流。     

南水北调中线工程交叉建筑物水毁风险分析

南水北调中线工程输水总干渠上需建数以百计的河渠交叉建筑物。总干渠经过地区为暴雨多发区，只要其中任一座遭遇特大洪水而失事，干渠供水将会受到影响。对这些交叉建筑物进行水毁风险分析，是此项工程规划、设计和运行管理的供输水风险评价和可靠度分析必不可少的组成部分。主要从水文风险要素的辨识入手，采用概率组合法，建立了整个总干渠洪水水毁风险计算的框架，并提出了二维复合事件的风险计算模型。所建模型巧妙地解决了各交     

过山气流与大兴安岭焚风

本文利用二维大气中尺度数值模式,模拟了冬季1月份大兴安岭东坡的焚风现象,并讨论了该焚风的成因。     

The hydrological response of soil surfaces to rainfall as affected by cover and position of rock fragments in the top layer

Rainfall experiments have been conducted in the laboratory in order to assess the hydrological response of top soils very susceptible to surface sealing and containing rock fragments in different positions with respect to the soil surface. For a given cover level, rock fragment position in the top soil has an ambivalent effect on water intake and runoff generation. Compared to a bare soil surface rock fragments increase water intake rates as well as time of runoff concentration and decrease runoff volume if they rest on the soil surface. For the same cover level, rock fragments reduce infiltration rate and enhance runoff generation if they are well embedded in the top layer. The effects of rock fragment position on infiltration rate and runoff generation are proportional to cover percentage. Micromorphological analysis and measurements of the saturated hydraulic conductivity of bare top soils and of the top layer underneath rock fragments resting on the soil surface reveal significant differences supporting the mechanism proposed by Poesen (1986): i.e. runoff generated as rock flow or as Horton overland flow can (partly) infiltrate into the unsealed soil surface under the rock fragments, provided that they are not completely embedded in the top layer. Hence, rock fragment position, beside other rock fragment properties, should be taken into account when assessing the hydrological response of soils susceptible to surface sealing and containing rock fragments in their surface layers. A simple model, based on the proportions of bare soil surface, soil surface occupied by embedded rock fragments, and soil surface covered with rock fragments resting on the soil surface, describes the runoff coefficient data relatively well.     

A swallow-tail type catastrophic model of earthquake process

（殷有泉）（杜静）Ａｓｗａｌｌｏｗ－ｔａｉｌｔｙｐｅｃａｔａｓｔｒｏｐｈｉｃｍｏｄｅｌｏｆｅａｒｔｈｑｕａｋｅｐｒｏｃｅｓｓ￥Ｙｏｕ－ＱｕａｎＹＩＮａｎｄＪｉｎｇＤＵ（ＰｅｋｉｎｇＵｎｉｖｅｒｓｄｙ，Ｂｅｉｊｉｎｇ１００８７１，Ｃｈｉｎａ）（Ｎｏｒｔ...     

水压致裂煤层裂缝发育特点的研究

本文对煤岩基本力学性质、煤层水压致裂缝形成条件、裂缝形态以及裂缝开裂角方位等基本理论进行了研究与探讨，通过大量煤岩力学性质测试，证实了试验区目的层煤岩弹性模量低，泊松比较高，脆性大，易破碎，易压缩。文章还得出了目的层煤岩Ｍｏｈｒ断裂准则二次抛物线型包络线，煤层水压致裂裂缝形式判断，裂缝开裂角方位的计算公式以及有关结论。