Large‐scale poroelastic fractured reservoirs modeling using the fast multipole displacement discontinuity method

An effective approach to modeling the geomechanical behavior of the network and its permeability variation is to use a poroelastic displacement discontinuity method (DDM). However, the approach becomes rather computationally intensive for an extensive system of cracks, particularly when considering coupled diffusion/deformation processes. This is because of additional unknowns and the need for time‐marching schemes for the numerical integration. The Fast Multipole Method (FMM) is a technique that can accelerate the solution of large fracture problems with linear complexity with the number of unknowns both in memory and CPU time. Previous works combining DDM and FMM for large‐scale problems have accounted only for elastic rocks, neglecting the fluid leak‐off from the fractures into the matrix and its influence on pore pressure and stress field. In this work we develop an efficient geomechanical model for large‐scale natural fracture networks in poroelastic reservoirs with fracture flow in response to injection and production operations. Accuracy and computational performance of the proposed method with those of conventional poroelastic DDM are compared through several case studies involving up to several tens of thousands of boundary elements. The results show the effectiveness of the FMM approach to successfully evaluate field‐scale problems for the design of exploitation strategies in unconventional geothermal and petroleum reservoirs. An example considering faults reveals the impact of reservoir compartmentalization because of sealing faults for both geomechanical and flow variables under elastic and poroelastic rocks. Copyright © 2015 John Wiley & Sons, Ltd.     

Watershed delineation using hydrographic features and a DEM in plain river network region

Watershed delineation is a required step when conducting any spatially distributed hydrological modelling. Automated approaches are often proposed to delineate a watershed based on a river network extracted from the digital elevation model (DEM) using the deterministic eight‐neighbour (D8) method. However, a realistic river network cannot be derived from conventional DEM processing methods for a large flat area with a complex network of rivers, lakes, reservoirs, and polders, referred to as a plain river network region (PRNR). In this study, a new approach, which uses both hydrographic features and DEM, has been developed to address the problems of watershed delineation in PRNR. It extracts the river nodes and determines the flow directions of the river network based on a vector‐based hydrographic feature data model. The river network, lakes, reservoirs, and polders are then used to modify the flow directions of grid cells determined by D8 approach. The watershed is eventually delineated into four types of catchments including lakes, reservoirs, polders, and overland catchments based on the flow direction matrix and the location of river nodes. Multiple flow directions of grid cells are represented using a multi‐direction encoding method, and multiple outflows of catchments are also reflected in the topology of catchments. The proposed approach is applied to the western Taihu watershed in China. Comparisons between the results obtained from the D8 approach, the ‘stream burning’ approach, and those from the proposed approach clearly demonstrate an improvement of the new approach over the conventional approaches. This approach will benefit the development of distributed hydrological models in PRNR for the consideration of different types and multiple inlets and outlets of catchments. Copyright © 2015 John Wiley & Sons, Ltd.     

Variability in the vertical hyporheic water exchange affected by hydraulic conductivity and river morphology at a natural confluent meander bend

River confluences and their associated tributaries are key morphodynamic nodes that play important roles in controlling hydraulic geometry and hyporheic water exchange in fluvial networks. However, the existing knowledge regarding hyporheic water exchange associated with river confluence morphology is relatively scarce. On January 14 and 15, 2016, the general hydraulic and morphological characteristics of the confluent meander bend (CMB) between the Juehe River and the Haohe River in the southern region of Xi'an City, Shaanxi Province, China, were investigated. The patterns and magnitudes of vertical hyporheic water exchange (VHWE) were estimated based on a one‐dimensional heat steady‐state model, whereas the sediment vertical hydraulic conductivity (K_v) was calculated via in situ permeameter tests. The results demonstrated that 6 hydrodynamic zones and their extensions were observed at the CMB during the test period. These zones were likely controlled by the obtuse junction angle and low momentum flux ratio, influencing the sediment grain size distribution of the CMB. The VHWE patterns at the test site during the test period mostly showed upwelling flow dominated by regional groundwater discharging into the river. The occurrence of longitudinal downwelling and upwelling patterns along the meander bend at the CMB was likely subjected to the comprehensive influences of the local sinuosity of the meander bend and regional groundwater discharge and finally formed regional and local flow paths. Additionally, in dominated upwelling areas, the change in VHWE magnitudes was nearly consistent with that in K_v values, and higher values of both variables generally occurred in erosional zones near the thalweg paths of the CMB, which were mostly made up of sand and gravel. This was potentially caused by the erosional and depositional processes subjected to confluence morphology. Furthermore, lower K_v values observed in downwelling areas at the CMB were attributed to sediment clogging caused by local downwelling flow. The confluence morphology and sediment K_v are thus likely the driving factors that cause local variations in the VHWE of fluvial systems.     

基于网络分析法的武汉市景观生态网络格局特征分析

从网络分析法的角度对武汉市水系网络、道路网络、绿地网络及其复合网络构成的景观生态网络格局进行分析。根据总体景观生态网络格局特征、分区统计景观生态网络格局特征、分等级统计景观生态网络格局特征分析结果显示,道路网络和绿地网络的α指数、β指数、γ指数都明显好于水系网络,相比较而言,三网复合网络都有比独立网络更好的连接度,分区统计结果显示,主城区复合网络的三项指标均处于最高,网络连通度最好,分等级统计结果显示,等级道路网络的α指数、β指数、γ指数均高于全部道路网络,等级水系和等级绿地网络的α指数、β指数、γ指数均低于全部水系和全部绿地网络。     

基于网络GPS-RTK技术的地籍测量研究

地籍测量工作是地籍管理和地籍信息系统建设的基础,在地籍测量土地状况利用调查中,查清各类土地的分布和利用状况需要有一个准确的定位。随着GPS测量技术的飞速发展和不断成熟,网络GPS-RTK技术以其高精度和高自动化优势在地籍测量中的作用越来越重要。在手机网络流行的今天,网络GPS的应用也越来越便利。本文结合河南理工大学新校区地籍测量工作,首先介绍网络RTK组成及定位原理,分析了网络RTK测图的关键技术,重点说明了网络GPS-RTK的数据采集作业方法以及数据处理和在地籍测量中的使用步骤,得出了网络GPS是一个快速和准确的测量方法,具有较强的实用性,可以满足地籍测量的需要。     

Matlab在矿区GPS网间接平差中的应用

针对矿区GPS网间接平差,首先,详细地分析了GPS网间接平差原理,按一定格式将GPS网中的观测数据编制成观测文件,其次,按观测文件的格式编制算法将数据读入并存储,然后用间接平差原理处理得到平差结果,最后绘制网形和误差椭圆并生成平差报告,为矿区控制网的建立和矿区生产管理提供了便利。实验验证了本文程序的实用性。     

含AgI焰剂成冰性能与物化特征的实验研究

利用CAMS的1 m~3等温云室系统筛选出新型高效AgI焰剂WMC-IN-001和WMG-IN-002。检测结果表明,它们具有较高的成核率,在-15℃时达到10~(15)g~(-1)AgI量级,尤其在-7℃时WMC-IN-001的成核率仍可达到10~(14)g~(-1)AgI量级。同时给出对节银剂配方和2011年市场上主要的几种催化剂的检测结果进行对比。WMC-IN-001和WMC-IN-002的成冰速率较慢,在各检测温度的成冰速率差异较小,均在40~55 min。利用冷场发射扫描电镜和能谱仪对WMC-IN-001燃烧产生的气溶胶粒子作了物化特征分析,粒子分布在0.02~0.60μm,具有两个典型的模态:0.02~0.10μm的较小的粒子和0.20~0.55μm的较大的粒子,均立方直径为0.2472μm。WMC-IN-001气溶胶粒子明显偏大,小粒子相对较少,这可能是其成冰速率偏慢的原因之一。     

从平衡-封闭-静态转向非平衡-开放-动力学研究-热液矿床成因的若干新概念概述

热液金属矿床成因研究过程中,观测与实验始终是密切结合的。上世纪70年代,平衡热力学的实验数据的快速积累,使人们用热力学理论计算可以预测和反演矿石和岩石的成因。但是,没有矿物-水溶液的反应速率数据,又没有与流体力学的结合,搞清楚矿石成因是困难的。七、八十年代,开始研究矿物与水溶液的反应动力学实验。科学家们开始瞄准了从平衡-封闭-静态转向非平衡-开放-动力学研究的这个大方向。1992年我们建立地球化学动力学开放研究实验室。研究高温高压矿物与水反应速率,发现固液的开放体系的自组织现象。实验发现温度影响矿物的各个元素反应速率改变,发现在跨越水临界态时矿物与水反应速率涨落、在近临界的气-液两相不混溶区一些金属进入气相、超临界流体的氧化作用及特别的溶剂性能影响矿物溶解性质。实验证实:临界态区流体与矿石成因有关。水岩相互作用的反应动力学实验温度从低温到550度,揭示矿石的金属来源、迁移、金属与蚀变分带机制。一大批大于300度的矿物与水反应动力学实验在国际界是少有报道的。九十年代,超高压的科学发展,与同步辐射光源的技术进步的结合,使固体地球科学又迈向了地球深内部。我们发展了高温高压流体性质的原位直测(测量850℃水溶液)红外谱,发现深部流体的新性质:气液两相流体的新结构,在临界温度区(300~400℃),水分子氢键网络的破坏受压力影响不大(23MPa~3GPa),同时,出现水的高电导率。研发新仪器为开放-流动-非平衡的反应动力学实验与极端条件下物质性质的直接观测结合,在科学前沿领域开辟了创新道路。     

Determination of the factors governing soil erodibility using hyperspectral visible and near-infrared reflectance spectroscopy

Soil erodibility, which is difficult to estimate and upscaling, was determined in this study using multiple spectral models of soil properties (soil organic matter (SOM), water-stable aggregates (WSA) > 0.25 mm, the geometric mean radius (D_g)). Herein, the soil erodibility indicators were calculated, and soil properties were quantitatively analyzed based on laboratory simulation experiments involving two selected contrasting soils. In addition, continuous wavelet transformation was applied to the reflectance spectra (350–2500 nm) of 65 soil samples from the study area. To build the relationship, the soil properties that control erodibility were identified prior to the spectral analysis. In this study, the SOM, D_g and WSA >0.25 mm were selected to represent the most significant soil properties controlling erodibility and describe the erodibility indicator based on a logarithmic regression model as a function of SOM or WSA > 0.25 mm. Five, six and three wavelet features were observed to calibrate the estimated soil properties model, and the best performance was obtained with a combination feature regression model for SOM (R² = 0.86, p < 0.01), D_g (R² = 0.79, p < 0.01) and WSA >0.25 mm (R² = 0.61, p < 0.01), respectively. One part of the wavelet features captured amplitude variations in the broad shape of the reflectance spectra, and another part captured variations in the shape and depth of the soil dry substances. The wavelet features for the validated dataset used to predict the SOM, WSA >0.25 mm and D_g were not significantly different compared with the calibrated dataset. The synthesized spectral models of soil properties, and the formation of a new equation for soil erodibility transformed from the spectral models of soil properties are presented in this study. These results show that a spectral analytical approach can be applied to complex datasets and provide new insights into emerging dynamic variation with erodibility estimation.