相关和频谱分析法在岩溶系统中的应用研究综述

岩溶含水层由多重含水介质组成,具有高度的非均质性,从系统角度出发,将岩溶系统看作是一个过滤器,将降雨作为输入信号,泉流量、水位、电导率、浊度等物理化学参数作为输出量进行相关和频谱分析是认识岩溶系统内部结构特征的有效方法。文章综述了相关和频谱分析法在岩溶系统研究中的应用,首先介绍了相关和频谱分析法基本理论的形成、发展过程和研究意义,然后总结了如何根据不同的研究目的选择相关和频谱分析法使用的时间序列的长度、周期和类型,最后阐述了该方法应用于岩溶系统分析的局限性。     

基于数值模拟的汾河二库对晋祠泉岩溶水系统渗漏补给作用研究

汾河二库的下闸蓄水,改变了晋祠泉域汾河地表水对岩溶水的渗漏补给条件,对区域岩溶水流场变化有着长远的影响。本研究通过系统梳理库区水文地质条件,勾画地表水渗漏补给岩溶水范围,考虑地层岩性、蓄水高度等因素,确定二库不同蓄水水位下晋祠泉域获渗漏量占比;构建晋祠泉三维岩溶地下水流模型,研究二库渗漏对泉域岩溶地下水补给作用,并预测不同蓄水情景下晋祠泉口水位变化趋势。结果表明：晋祠泉域获渗漏量占比随二库蓄水位非线性变化,最低为92.8%,在二库蓄水位达902 m后稳定在93.7%,而渗漏量与二库蓄水高度呈正相关;泉域岩溶水径流区受渗漏补给作用最显著,其次为排泄区、北部补给区;渗漏补给作用下,地下水位回升值受二库水位和补水距离影响;维持设计水位（905.7 m）和2017年水位（895.9~902.4 m）时,预计晋祠泉口水位分别在2021年7月和2023年1月达到泉口高程。     

关于发展我国地下流体前兆流动观测问题的建议

从主机故障、传感器故障、避雷系统故障、供电系统故障4个方面,对天津流体台网水温观测系统故障事件及其跟踪分析进行讨论,提出规范化的排查思路及记录处理方法,可供兄弟台站作为同类数据跟踪分析工作的参考。     

Precessional forcing of tropical vegetation carbon storage

Since the Mid Pleistocene Revolution, which occurred about one million years ago, global temperatures have fluctuated with a quasi‐periodicity of ca. 100 ka. The pattern of past change in the extent of woodlands, and therefore by inference vegetation carbon storage, has been demonstrated to have a strong positive link with this global temperature change at high and mid latitudes. However, understanding of climate systems and ecosystem function indicates that the pattern of woodland change at low latitudes may follow a fundamentally different pattern. We present output from the intermediate complexity model GENIE‐1, comprising a single transient simulation over the last 800 ka and a 174‐member ensemble of 130 ka transient simulations over the last glacial cycle. These simulations suggest that while vegetation carbon storage in mid–high northern latitudes robustly follows the characteristic ca. 100 ka cycle, this signal is not a robust feature of tropical vegetation, which is subject to stronger direct forcing by the precessional (21 ka) orbital cycle (albeit with a highly uncertain response). We conclude that the correlation of palaeoenvironmental records from low latitudes with global temperature change must be done with caution. Copyright © 2011 John Wiley & Sons, Ltd.     

The influence of coupled physical swelling and chemical reactions on deformable geomaterials

Coupled thermo‐hydro‐mechanical‐chemical modelling has attracted attention in past decades due to many contemporary geotechnical engineering applications (e.g., waste disposal, carbon capture and storage). However, molecular‐scale interactions within geomaterials (e.g., swelling and dissolution/precipitation) have a significant influence on the mechanical behaviour, yet are rarely incorporated into existing Thermal‐Hydro‐Mechanical‐Chemical (THMC) frameworks. This paper presents a new coupled hydro‐mechanical‐chemical constitutive model to bridge molecular‐scale interactions with macro‐physical deformation by combining the swelling and dissolution/precipitation through an extension of the new mixture‐coupling theory. Entropy analysis of the geomaterial system provides dissipation energy, and Helmholtz free energy gives the relationship between solids and fluids. Numerical simulation is used to compare with the selected recognized models, which demonstrates that the swelling and dissolution/precipitation processes may have a significant influence on the mechanical deformation of the geomaterials.     

Elastic behavior of weakly cemented contact

We investigate the elastic behavior of weakly cemented contact. We show that the radial distribution of stresses and the stiffnesses of a cemented contact are governed by the ratio a/RΛ, where R, a, and Λ are, respectively, the grain radius, the contact size, and the ratio of the elastic moduli of cement and grains. Moreover, we show that a cemented contact is always less stiff than a Hertzian contact having a similar size. Finally, we propose accurate approximate expressions of the contact stiffnesses. Copyright © 2014 John Wiley & Sons, Ltd.     

Why asymptotic limit of the acid dissolution capacity can lead to a sharp dissolution front in chemical dissolution of porous rocks?

The use of the asymptotic limit can greatly simplify the theoretical analysis of chemical dissolution front instabilities in fluid‐saturated rocks and therefore make it possible to obtain mathematical solutions, which often play a crucial role in understanding the propagation behavior of chemical dissolution fronts in chemical dissolution systems. However, there has been a debate in recent years that the asymptotic limit of the acid dissolution capacity (i.e., the acid dissolution capacity number approaching zero) alone cannot lead to a sharp dissolution front of the Stefan type in the acidization dissolution system, in which the dissolvable minerals of carbonate rocks are chemically dissolved by the injected acid flow. The acid dissolution capacity number is commonly defined as the ratio of the volume of the carbonate rock dissolved by an acid to that of the acid. In this paper, we use four different proof methods, including (i) direct use of the fundamental concepts; (ii) use of the mathematical governing equations of an acidization dissolution system; (iii) use of the different time scaling approach; and (iv) use of a moving coordinate system approach, to demonstrate that the asymptotic limit of the acid dissolution capacity can indeed lead to sharp dissolution fronts of the Stefan type in acidization dissolution systems on a much larger time scale (than the dissolution time scale). Our new finding is that on the reaction time scale, the condition of the conventional time derivative of porosity approaching zero alone can ensure that the acidization dissolution front has a sharp shape of the Stefan type. Copyright © 2017 John Wiley & Sons, Ltd.     

A network theory approach for a better understanding of overland flow connectivity

Hydrological connectivity describes the physical coupling (linkages) of different elements within a landscape regarding (sub‐) surface flows. A firm understanding of hydrological connectivity is important for catchment management applications, for example, habitat and species protection, and for flood resistance and resilience improvement. Thinking about (geomorphological) systems as networks can lead to new insights, which has also been recognized within the scientific community, seeing the recent increase in the use of network (graph) theory within the geosciences. Network theory supports the analysis and understanding of complex systems by providing data structures for modelling objects and their linkages, and a versatile toolbox to quantitatively appraise network structure and properties. The objective of this study was to characterize and quantify overland flow connectivity dynamics on hillslopes in a humid sub‐Mediterranean environment by using a combination of high‐resolution digital‐terrain models, overland flow sensors and a network approach. Results showed that there are significant differences between overland flow connectivity on agricultural areas and semi‐natural shrubs areas. Significant positive correlations between connectivity and precipitation characteristics were found. Significant negative correlations between connectivity and soil moisture were found, most likely because of soil water repellency and/or soil surface crusting. The combination of structural networks and dynamic networks for determining potential connectivity and actual connectivity proved a powerful tool for analysing overland flow connectivity. Copyright © 2016 John Wiley & Sons, Ltd.     

An analytical study on three‐dimensional versus two‐dimensional water table‐induced flow patterns in a Tóthian basin

Although it has been increasingly acknowledged that groundwater flow pattern is complicated in the three‐dimensional (3‐D) domain, two‐dimensional (2‐D) water table‐induced flow models are still widely used to delineate basin‐scale groundwater circulation. However, the validity of 2‐D cross‐sectional flow field induced by water table has been seldom examined. Here, we derive the analytical solution of 3‐D water table‐induced hydraulic head in a Tóthian basin and then examine the validity of 2‐D cross‐sectional models by comparing the flow fields of selected cross sections calculated by the 2‐D cross‐sectional model with those by the 3‐D model, which represents the “true” cases. For cross sections in the recharge or discharge area of the 3‐D basin, even if head difference is not significant, the 2‐D cross‐sectional models result in flow patterns absolutely different from the true ones. For the cross section following the principal direction of groundwater flow, although 2‐D cross‐sectional models would overestimate the penetrating depth of local flow systems and underestimate the recharge/discharge flux, the flow pattern from the cross‐sectional model is similar to the true one and could be close enough to the true one by adjusting the decay exponent and anisotropy ratio of permeability. Consequently, to determine whether a 2‐D cross‐sectional model is applicable, a comparison of hydraulic head difference between 2‐D and 3‐D solutions is not enough. Instead, the similarity of flow pattern should be considered to determine whether a cross‐sectional model is applicable. This study improves understanding of groundwater flow induced by more natural water table undulations in the 3‐D domain and the limitations of 2‐D models accounting for cross‐sectional water table undulation only.