横观各向同性饱和土体三维粘弹性动力分析

采用针对横观各向同性饱和土体u-w形式三维粘弹性动力方程,考虑土骨架的粘弹性性质且基于粘弹性理论,通过运用Fourier 展开、Laplace 和Hankel 积分变换方法和引入中间变量,将含有粘弹性参数的六元二阶偏微分运动控制方程组,化为2组各含4个未知变量的常微分方程组,从而给出了柱坐标系下粘弹性横观各向同性饱和土体在非轴对称动力荷载作用下的瞬态反应的土骨架位移分量、孔隙流体相对于土骨架的位移分量瞬态反应一般解。在此基础上,引入初始条件和边界条件,对垂直向和水平向动力荷载作用下半空间边值问题进行了求解。根据动力时域解答的一般解,利用Laplace和Hankel 数值逆变换技术,编制了相应的数值计算程序。并进行了实例验证和弹性、粘弹性解的对比分析。结果表明,在进行横观各向同性饱和土体动力分析时,考虑土骨架的粘弹性是必要的。     

岩盐中不溶物含量对储气库腔体有效空间大小的影响

主要研究了岩盐不溶物含量对岩盐储气库水溶建腔有效空间大小的影响。根据水溶建腔机理,建立岩盐储气库水溶建腔数学模型,通过岩盐溶腔数学模型,分别对岩盐中不溶物含量不同的地层条件下,岩盐储气库水溶建腔进行动态模拟数值计算,综合分析岩盐中不溶物含量对岩盐储气库水溶建腔的影响规律。     

Analytical solutions for the consolidation of unsaturated foundation with prefabricated vertical drain

Based on the consolidation theory raised by Fredlund, the solutions for the equal-strain consolidation of unsaturated foundation with the prefabricated vertical drain considering smear effect and drain resistance are analytically formulated in this paper. Firstly, governing equations for excess pore pressures (i.e., excess pore-air and pore-water pressures) under the equal-strain hypothesis are derived with the introduction of radial boundary conditions. Afterwards, the obtained coupled equations are solved by applying general integration, decoupling process, and Fourier sine series expansion. The smear coefficients and factors of drain resistance corresponding to air and water phases are both captured explicitly in the final solutions. Furthermore, the degenerated solutions are employed to verify the reliability of the current solutions. Finally, a parametric study is conducted to study the consolidation characteristics of the proposed foundation model against modeling sizes (S and N), smear coefficients (α_a and α_w), and drain resistance factors (G_a and G_w).     

Similarity search and pattern discovery in hydrological time series data mining

The rapid development of data mining provides a new method for water resource management, hydrology and hydroinformatics research. In the paper, based on data mining theory and technology, we analyse hydrological daily discharge time series of the Shaligunlanke Station in the Tarim River Basin in China from the year 1961 to 2000. Firstly, according to the four monthly statistics, namely mean monthly discharge, monthly maximum discharge, monthly amplitude and monthly standard deviation, K‐mean clustering was used to segment the annual process of the daily discharge. The clustering result showed that the annual process of the daily discharge can be divided into five segments: snowmelt period I (April), snowmelt period II (May), rainfall period I (June–August), rainfall period II (September) and dry period (October–December and January–March). Secondly, dynamic time warping (DTW), which is a different distance metric method from the traditional Euclidian distance metric, was used to look for similarities in the discharge process. On the basis of the similarity matrix, the similar discharge processes can be mined in each period. Thirdly, agglomerative hierarchical clustering was used to cluster and discover the discharge patterns in terms of the autoregressive model. It was found that the discharge had a close relationship with the temperature and the precipitation, and the discharge processes were more similar under the same climatic condition. Our study shows that data mining is a feasible and efficient approach to discover the hidden information in the historical hydrological data and mining the implicative laws under the hydrological process. Copyright © 2010 John Wiley & Sons, Ltd.     

A simple expression for the bulk field capacity of a sloping soil horizon

Field capacity is a commonly used soil parameter in surface water hydrological models, loosely defined as the moisture content of a soil after drainage. The most commonly applied expression for field capacity is defined as the remaining water in a vertical soil column subject to 1/3 atm. of pressure head. While this quantification is sufficient in some cases, the definition is not consistent with the use of bulk field capacity in calculations of lateral drainage from hillslopes, as required by some surface soil parameterizations, nor does it address additional complications arising from differences in soil texture or sample size. Here, a simple alternative expression for bulk field capacity in a sloping or vertical soil is derived directly from Richards equation with the use of the Brooks‐Corey characteristics. It is demonstrated that this expression is consistent with data acquired from vertical soil columns, but may be extended to additional situations commonly found in surface water models and land surface schemes. The calculation of bulk field capacity requires only the Brooks‐Corey pore size distribution index, soil air‐entry pressure, and hillslope length and slope, and may be considered a physically based alternative to pedotransfer function or lookup table approaches. Copyright © 2010 John Wiley & Sons Ltd and Crown in the right of Canada.     

A high‐resolution local RBF collocation method for steady‐state poroelasticity and hydromechanical damage analysis

In this work, we describe a meshless numerical method based on local collocation with RBFs for the solution of the poroelasticity equation. The RBF finite collocation approach forms a series of overlapping nodal stencils, over which an RBF collocation is performed. These local collocation systems enforce the governing PDE operator throughout their interior, with the intersystem communication occurring via the collocation of field variables at the stencil periphery. The method does not rely on a generalised finite differencing approach, whereby the governing partial differential operator is reconstructed at the global level to drive the solution of the PDE. Instead, the PDE governing and boundary operators are enforced directly within the local RBF collocation systems, and the sparse global assembly is formed by reconstructing the value of the field variables at the centrepoint of the local stencils. In this way, the solution of the PDE is driven entirely by the local RBF collocation, and the method more closely resembles the approach of the full‐domain RBF collocation method. By formulating the problem in this fashion, high rates of convergence may be attained without the computational cost and numerical ill‐conditioning issues that are associated with the full‐domain RBF collocation approach. An analytical solution is formulated for a 2D poroelastic fluid injection scenario and is used to verify the proposed implementation of the method. Highly accurate solutions are produced, and convergence rates in excess of sixth order are observed for each field variable (i.e. pressure and displacement) and field‐variable derivative (i.e. pressure gradients and stresses). The stress and displacement fields resulting from the solution of the poroelasticity equation are then used to describe the formation and propagation of microfractures and microfissures, which may form in the presence of large shear strain, in terms of a continuous damage variable which modifies the mechanical and hydraulic properties of the porous medium. The formation of such hydromechanical damage, and the resulting increase in hydraulic conductivity, is investigated for a pressurised injection into sandstone. Copyright © 2014 John Wiley & Sons, Ltd.     

A soft hydrological monitoring approach for comparing runoff on a network of small poorly gauged catchments

Catchments in many parts of the world are either ungauged or poorly gauged, and the dominant processes governing their streamflow response are still poorly understood. The analysis of runoff coefficients provides essential insight into catchment response, particularly if both range of catchments and a range of events are compared. This paper investigates how well the hydrological runoff of 11 small, poorly gauged catchments with ephemeral streams (0·1‐0·6 km²) can be compared using estimated runoff with the associated uncertainty. Data of rainfall and water depth at a catchment's outlet were recorded using automatic logging equipment during 2008‐2009. The hydrological regime is intermittent and the annual precipitation ranged between 569 and 727 mm. Discharge was estimated using Manning's equation and channel cross‐section measurements. Innovative work has been performed under controlled experimental conditions to estimate Manning's coefficient values for the different cover types observed in studied streams: non‐aquatic vegetations (giant reed, bramble and thistle), grass and coarse granular deposits. The results show that estimates derived using roughness coefficients differ from those previously established for larger streams with aquatic vegetation. Catchment runoff was compared at both the event and the annual scale. The results indicate significant variability between the catchment's responses. This variability allows for classification in spite of all the uncertainty associated with runoff estimation. This study highlights the potential of using a network of poorly gauged catch ments. From almost no catchment understanding the proposed methodology allows to compare poorly gauged catchments and highlights similarity/dissimilarity between catchment responses. Copyright © 2011 John Wiley & Sons, Ltd.     

The effect of temperature‐humidity similarity on Bowen ratios,dimensionless standard deviations,and mass transfer coefficients over a lake

Similarity between heat and water vapor turbulent transport in the Atmospheric Surface Layer has been the basis of many engineering models to calculate surface fluxes, including the widely applied Bowen ratio equation, for a long time. Modernly, it is best understood within the context of Monin‐Obkhov Similarity Theory (MOST). In this work we study similarity between temperature and humidity, the Bowen ratio, and turbulent mass and heat transfer coefficients over a tropical lake in Brazil (Furnas Lake). The analysis was partly based on the concept of ‘Surface flux numbers’ recently proposed to diagnose scalar similarity, and considered wind directions and flux footprints. A period of 50 days of 30‐min. micrometeorological runs was used. Several cases of dissimilar temperature‐humidity behavior were found in the data. Both footprint extent and an aggregate temperature‐humidity Surface flux number turned out to be insufficient to diagnose these situations, but separate flux numbers for each scalar were able to diagnose their individual conformity to MOST. Overall, temperature displayed consistently larger relative variances and fluxes in comparison with humidity. The results highlight the need of careful analysis when measurements are made at sites close to land, when flux footprints may extend over there, indicating the possibility of advection effects. Copyright © 2016 John Wiley & Sons, Ltd.     

Mechanical and transport constitutive models for fractures subject to dissolution and precipitation

Transient changes in the permeability of fractures in systems driven far‐from‐equilibrium are described in terms of proxy roles of stress, temperature and chemistry. The combined effects of stress and temperature are accommodated in the response of asperity bridges where mineral mass is mobilized from the bridge to the surrounding fluid. Mass balance within the fluid accommodates mineral mass either removed from the flow system by precipitation or advection, or augmented by either dissolution or advection. Where the system is hydraulically closed and initially at equilibrium, reduction in aperture driven by the effects of applied stresses and temperatures will be augmented by precipitation on the fracture walls. Where the system is open, the initial drop in aperture may continue, and accelerate, where the influent fluid is oversaturated with respect to the equilibrium mineral concentration within the fluid, or may reverse, if undersaturated. This simple zero‐dimensional model is capable of representing the intricate behavior observed in experiments where the feasibility of fracture sealing concurrent with net dissolution is observed. This zero‐order model is developed as a constitutive model capable of representing key aspects of changes in the transport parameters of the continuum response of fractured media to changes in stress, temperature and chemistry. Copyright © 2009 John Wiley & Sons, Ltd.     

Micromechanical modeling for rate‐dependent behavior of salt rock under cyclic loading

The dependence of rock behavior on the deformation rate is still not well understood. In salt rock, the fundamental mechanisms that drive the accumulation of irreversible deformation, the reduction of stiffness, and the development of hysteresis during cyclic loading are usually attributed to intracrystalline plasticity and diffusion. We hypothesize that at low pressure and low temperature, the rate‐dependent behavior of salt rock is governed by water‐assisted diffusion along grain boundaries. Accordingly, a chemo‐mechanical homogenization framework is proposed in which the representative elementary volume (REV) is viewed as a homogeneous polycrystalline matrix that contains sliding grain‐boundary cracks. The slip is related to the mass of salt ions that diffuse along the crack surface. The relationship between fluid inclusion‐scale and REV‐scale stresses and strains is established by using the Mori–Tanaka homogenization scheme. It is noted from the model that a lower strain rate and a larger number of sliding cracks enhance stiffness reduction and hysteresis. Thinner sliding cracks (i.e., thinner brine films) promote stiffness reduction and accelerate stress redistributions. The larger the volume fraction of the crack inclusions, the larger the REV deformation and the larger the hysteresis. Results presented in this study shed light on the mechanical behavior of salt rock that is pertinent to the design of geological storage facilities that undergo cyclic unloading, which could help optimize the energy production cycle with low carbon emissions.