基于"3S"技术的于田绿洲湿地动态变化研究

以极端干旱区典型绿洲———于田绿洲为研究区,采用遥感(RS)、地理信息系统(GIS)和全球定位系统(GPS)相结合的方法,进行了湿地资源调查,查明了于田绿洲湿地类型、区域分布及面积,应用模型和分形理论,探讨了于田绿洲湿地的动态变化特征,并分析了其发生动态变化的原因。研究结果对于干旱区绿洲湿地的合理开发和可持续发展有着重要现实意义。     

Impacts of groundwater flow on the evolution of a thermokarst lake in the permafrost–dominated region on the Qinghai–Tibet plateau

As a result of global warming induced permafrost degradation in recent decades, thermokarst lakes in the Qinghai–Tibet plateau (QTP) have been regulating local hydrological and ecological processes. Simulations with coupled moisture–heat numerical models in the Beiluhe basin (located in the hinterland of permafrost regions on the QTP) have provided insights into the interaction between groundwater flow and the freeze–thaw process. A total of 30 modified SUTRA scenarios were established to examine the effects of hydrodynamic forces, permeability, and climate on thermokarst lakes. The results indicate that the hydrodynamic condition variables regulate the permafrost degradation around the lakes. In case groundwater recharges to the lake, a low–temperature groundwater flow stimulates the expansion of the surrounding thawing regions through thermal convection. The thawing rate of the permafrost underlying the lake intensifies when groundwater is discharged from the lake. Under different permeability conditions, spatiotemporal variations in the active layer thickness significantly influence the occurrence of an open talik at the lake bottom. A warmer and wetter climate will inevitably lead to a sharp decrease in the upper limit of the surrounding permafrost, with a continual decrease in the duration of open talik events. Overall, our results underscore that comprehensive consideration of the relevant hydrologic processes is critical for improving the understanding of environmental and ecological changes in cold environments.     

Numerical modeling of toxic nonaqueous phase liquid removal from contaminated groundwater systems: mesh effect and discretization error estimation

Numerical modeling has now become an indispensable tool for investigating the fundamental mechanisms of toxic nonaqueous phase liquid (NAPL) removal from contaminated groundwater systems. Because the domain of a contaminated groundwater system may involve irregular shapes in geometry, it is necessary to use general quadrilateral elements, in which two neighbor sides are no longer perpendicular to each other. This can cause numerical errors on the computational simulation results due to mesh discretization effect. After the dimensionless governing equations of NAPL dissolution problems are briefly described, the propagation theory of the mesh discretization error associated with a NAPL dissolution system is first presented for a rectangular domain and then extended to a trapezoidal domain. This leads to the establishment of the finger‐amplitude growing theory that is associated with both the corner effect that takes place just at the entrance of the flow in a trapezoidal domain and the mesh discretization effect that occurs in the whole NAPL dissolution system of the trapezoidal domain. This theory can be used to make the approximate error estimation of the corresponding computational simulation results. The related theoretical analysis and numerical results have demonstrated the following: (1) both the corner effect and the mesh discretization effect can be quantitatively viewed as a kind of small perturbation, which can grow in unstable NAPL dissolution systems, so that they can have some considerable effects on the computational results of such systems; (2) the proposed finger‐amplitude growing theory associated with the corner effect at the entrance of a trapezoidal domain is useful for correctly explaining why the finger at either the top or bottom boundary grows much faster than that within the interior of the trapezoidal domain; (3) the proposed finger‐amplitude growing theory associated with the mesh discretization error in the NAPL dissolution system of a trapezoidal domain can be used for quantitatively assessing the correctness of computational simulations of NAPL dissolution front instability problems in trapezoidal domains, so that we can ensure that the computational simulation results are controlled by the physics of the NAPL dissolution system, rather than by the numerical artifacts. Copyright © 2014 John Wiley & Sons, Ltd.     

Dynamic storage: a potential metric of inter‐basin differences in storage properties

The potential for dynamic storage to serve as a metric of basin behaviour was assessed using data from five drainage basins with headwaters on the thick sand and gravel deposits of the Oak Ridges Moraine in southern Ontario, Canada. Dynamic storage was directly correlated with the ratio of variability of δ²H in streamflow relative to that in precipitation. This ratio has previously been shown to be inversely related to basin mean transit time (MTT), suggesting an inverse relationship between dynamic storage and MTT for the study basins. Dynamic storage was also directly correlated with interannual variability in stream runoff, baseflow and baseflow:runoff ratio, implying that basins with smaller dynamic storage have less interannual variability in their streamflow regimes. These preliminary results suggest that dynamic storage may serve as a readily derived and useful metric of basin behaviour for inter‐basin comparisons. Copyright © 2016 John Wiley & Sons, Ltd.     

Elastoplastic model with loading memory surfaces (LMS) for monotonic and cyclic behaviour of geomaterials

A new elastoplastic model called loading memory surface based on the critical state concept and the multi‐surface framework is proposed for geomaterials. The model uses a hypoelastic formulation and two plastic mechanisms. The formulations of the model are made in three‐dimensional stress–strain space and work under both monotonic and cyclic loadings. A newly introduced formalism makes it possible to obtain the cyclic response directly from the monotonic loading one. This formalism gives a three‐dimensional generalization of the well‐known Masing rule. The model has been validated against test results of Hostun sand under several conditions: monotonic and cyclic, drained and undrained, tests in compression and in extension, and at different confining pressures and different densities. Copyright © 2014 John Wiley & Sons, Ltd.     

Elastic lateral dynamic impedance functions for a rigid cylindrical shell type foundation

Elastic lateral dynamic impedance functions are defined as the ratio of the lateral dynamic force/moment to the corresponding lateral displacement/rotation at the top ending of a foundation at very small strains. Elastic lateral dynamic impedance functions have a defining influence on the natural frequencies of offshore wind turbines supported on cylindrical shell type foundations, such as suction caissons, bucket foundations, and monopiles. This paper considers the coupled horizontal and rocking vibration of a cylindrical shell type foundation embedded in a fully saturated poroelastic seabed in contact with a seawater half‐space. The formulation of the coupled seawater–shell–seabed vibration problem is simplified by treating the shell as a rigid one. The rigid shell vibration problem is approached by the integral equation method using ring‐load Green's functions for a layered seawater‐seabed half‐space. By considering the boundary conditions at the shell–soil interface, the shell vibration problem is reduced to Fredholm integral equations. Through an analysis of the corresponding Cauchy singular equations, the intrinsic singular characteristics of the problem are rendered explicit. With the singularities incorporated into the solution representation, an effective numerical method involving Gauss–Chebyshev method is developed for the governing Fredholm equations. Selected numerical results for the dynamic contact load distributions, displacements of the shell, and lateral dynamic impedance functions are examined for different shell length–radius ratio, poroelastic materials, and frequencies of excitation. Copyright © 2016 John Wiley & Sons, Ltd.     

Dynamic analysis of strain localization in water‐saturated clay using a cyclic elasto‐viscoplastic model

A computational framework is presented for dynamic strain localization and deformation analyses of water‐saturated clay by using a cyclic elasto‐viscoplastic constitutive model. In the model, the nonlinear kinematic hardening rule and softening due to the structural degradation of soil particles are considered. In order to appropriately simulate the large deformation phenomenon in strain localization analysis, the dynamic finite element formulation for a two‐phase mixture is derived in the updated Lagrangian framework. The shear band development is shown through the distributions of viscoplastic shear strain, the axial strain, the mean effective stress, and the pore water pressure in a normally consolidated clay specimen. From the local stress–strain relations, more brittleness is found inside the shear bands than outside of them. The effects of partially drained conditions and mesh‐size dependency on the shear banding are also investigated. The effect of a partially drained boundary is found to be insignificant on the dynamic shear band propagation because of the rapid rate of applied loading and low permeability of the clay. Using the finer mesh results in slightly narrower shear bands; nonetheless, the results manifest convergency through the mesh refinement in terms of the overall shape of shear banding and stress–strain relations. Copyright © 2013 John Wiley & Sons, Ltd.     

Vertical vibration of a rigid strip footing on viscoelastic half-space

This paper presents an analytical method for modeling the dynamic response of a rigid strip footing subjected to vertical-only loads. The footing is assumed to rest on the surface of a viscoelastic half-space; therefore, effects of hysteretic soil damping on the impedance of the foundation and the generated ground vibrations are considered in the solution. To solve the mixed boundary value problem, we use the Fourier transform to cast a pair of dual integral equations providing contact stresses, which are solved by means of Jacobi orthogonal polynomials. The resulting soil and footing displacements and stresses are obtained by means of the Fourier inverse transform. The solution provides more realistic estimates of footing impedance, compared to existing solutions for elastic soil, as well as of the attenuation of ground vibrations with distance of the footing. The latter is important for the estimation of machine vibration effects on nearby structures and installations.     

Influence of a coarse interlayer on seawater intrusion and contaminant migration in coastal aquifers

Vertical 2D slice laboratory experiments were carried out in homogenous and layered sand tanks to elucidate the effects of a highly permeable (coarse‐grained sand) interlayer on seawater intrusion and transport of contaminants to a coastal sea. Tidal fluctuations produced oscillations in the seawater–freshwater transition zone, fluctuations of the contaminant infiltration rate and a zigzag contaminant plume outline. The seawater wedge became discontinuous at the (vertical) edges of the interlayer because of increased lateral movement of the seawater–freshwater interface within the interlayer. The contaminant plume formed a tail within the interlayer depending on the tidal stage, and similar to the wedge, its movement was accentuated. A simple analytical model that neglected vertical flow reliably predicted steady‐state seawater intrusion into the coastal aquifer. Numerical modeling was used to gain insight into the groundwater hydrodynamics and contaminant migration. The numerical results confirmed the experimental findings, i.e. that a highly permeable interlayer can provide a rapid transit path for contaminants to reach the seaward boundary and that the interlayer amplifies the effects of tidal fluctuations, resulting in wider transition zones for the seawater wedge and contaminant plume. Numerical simulations further showed that, with increasing interlayer hydraulic conductivity, the maximum seawater intrusion distance inside the interlayer increases approximately linearly. For the fixed‐head contaminant injection condition used, the model showed that contaminant infiltration increases approximately logarithmically with increasing interlayer hydraulic conductivity (other factors held fixed). Copyright © 2013 John Wiley & Sons, Ltd.     

Dynamic response of a pile considering the interaction of pile variable cross section with the surrounding layered soil

Assuming that the pile variable cross section interacts with the surrounding soil in the same way as the pile toe does with the bearing stratus, the interaction of pile variable cross section with the surrounding soil is represented by a Voigt model, which consists of a spring and a damper connected in parallel, and the spring constant and damper coefficient are derived. Thus, a more rigid pile–soil interaction model is proposed. The surrounding soil layers are modeled as axisymmetric continuum in which its vertical displacements are taken into account and the pile is assumed to be a Rayleigh–Love rod with material damping. Allowing for soil properties and pile defects, the pile–soil system is divided into several layers. By means of Laplace transform, the governing equations of soil layers are solved in frequency domain, and a new relationship linking the impedance functions at the variable‐section interface between the adjacent pile segments is derived using a Heaviside step function, which is called amended impedance function transfer method. On this basis, the impedance function at pile top is derived by amended impedance function transfer method proposed in this paper. Then, the velocity response at pile top can be obtained by means of inverse Fourier transform and convolution theorem. The effects of pile–soil system parameters are studied, and some conclusions are proposed. Then, an engineering example is given to confirm the rationality of the solution proposed in this paper. Copyright © 2017 John Wiley & Sons, Ltd.