Using spatial statistics to identify the uranium hotspot in groundwater in the mid-eastern Gangetic plain,India

The occurrence of uranium in groundwater is of particular interest due to its toxicological and radiological properties. It has been considered as a relevant contaminant for drinking water even at a low concentration. Uranium is a ubiquitously occurring radionuclide in the environment. Four hundred and fifty-six (456) groundwater samples from different locations of five districts of South Bihar (SB) were collected and concentrations of uranium (U) were analyzed using a light-emitting diode (LED) fluorimetric technique. Uranium concentrations in groundwater samples varied from 0.1 µg l^?1 to 238.2 µg l^?1 with an average value of 12.3 µg l^?1 in five districts of Bihar in the mid-eastern Gangetic plain. This study used hot spot spatial statistics to identify the distribution of elevated uranium concentration in groundwater. The hypothesis whether spatial distribution of high value and low value of U is more likely spatially clustered due to random process near a uranium hotspot in groundwater was tested based on z score and Getis-Ord Gi* statistics. The method implemented in this study, can be utilized in the field of risk assessment and decision making to locate potential areas of contamination.     

Laboratory investigations of inert gas flow behaviors in compact sandstone

The seepage evolution behavior of compact rock is significant for the stability and safety of many engineering applications. In this research, both hydrostatic and triaxial compression tests were conducted on compact sandstone using an inert gas, namely argon. A triaxial compression test with a water permeability measurement was carried out to study the difference between the gas permeability and water permeability evolutions during the complete stress–strain process. Based on the experimental data, the hydrostatic stress-dependent gas permeability was discussed firstly. A second-order function was proposed to predict and explain the gas slippage effect. The mechanical properties and crack development of the sandstone samples were discussed to better understand the permeability evolution with crack growth during the complete stress–strain process. The results show that the gas permeability evolution can be divided into five stages according to the different crack growth stages. Then, the permeability changes in the crack closure stress \( \sigma_{\text{cc}} \), crack initiation stress \( \sigma_{\text{ci}} \), crack damage stress \( \sigma_{\text{cd}} \) and peak stress \( \sigma_{\text{p}} \) with confining pressures were analyzed. Finally, we found that the difference between the corrected gas permeability and water permeability can be attributed to the interaction between the water and sandstone grains.     

Interior boundary-aligned unstructured grid generation and cell-centered versus vertex-centered CVD-MPFA performance

Grid generation for reservoir simulation must honor classical key constraints and be boundary aligned such that control-volume boundaries are aligned with geological features such as layers, shale barriers, fractures, faults, pinch-outs, and multilateral wells. An unstructured grid generation procedure is proposed that automates control-volume and/or control point boundary alignment and yields a PEBI-mesh both with respect to primal and dual (essentially PEBI) cells. In order to honor geological features in the primal configuration, we introduce the idea of protection circles, and to generate a dual-cell feature based grid, we construct halos around key geological features. The grids generated are employed to study comparative performance of cell-centred versus cell-vertex control-volume distributed multi-point flux approximation (CVD-MPFA) finite-volume formulations using equivalent degrees of freedom. The formulation of CVD-MPFA schemes in cell-centred and cell-vertex modes is analogous and requires switching control volume from primal to dual or vice versa together with appropriate data structures and boundary conditions. The relative benefits of both types of approximation, i.e., cell-centred versus vertex-centred, are made clear in terms of flow resolution and degrees of freedom required.     

Constraint energy minimizing generalized multiscale finite element method in the mixed formulation

This paper presents a novel mass-conservative mixed multiscale method for solving flow equations in heterogeneous porous media. The media properties (the permeability) contain multiple scales and high contrast. The proposed method solves the flow equation in a mixed formulation on a coarse grid by constructing multiscale basis functions. The resulting velocity field is mass-conservative on the fine grid. Our main goal is to obtain first-order convergence in terms of the mesh size which is independent of local contrast. This is achieved, first, by constructing some auxiliary spaces, which contain global information that cannot be localized, in general. This is built on our previous work on the generalized multiscale finite element method (GMsFEM). In the auxiliary space, multiscale basis functions corresponding to small (contrast-dependent) eigenvalues are selected. These basis functions represent the high-conductivity channels (which connect the boundaries of a coarse block). Next, we solve local problems to construct multiscale basis functions for the velocity field. These local problems are formulated in the oversampled domain, taking into account some constraints with respect to auxiliary spaces. The latter allows fast spatial decay of local solutions and, thus, allows taking smaller oversampled regions. The number of basis functions depends on small eigenvalues of the local spectral problems. Moreover, multiscale pressure basis functions are needed in constructing the velocity space. Our multiscale spaces have a minimal dimension, which is needed to avoid contrast dependence in the convergence. The method’s convergence requires an oversampling of several layers. We present an analysis of our approach. Our numerical results confirm that the convergence rate is first order with respect to the mesh size and independent of the contrast.     

Residual-state creep of clastic soil in a reactivated slow-moving landslide in the Three Gorges Reservoir Region,China

We study the creep properties of clastic soil in residual state. The intact samples are taken from a reactivated slow-moving landslide in the Three Gorges Reservoir Region, China. Firstly, the patterns of the landslide movement are analysed based on recent monitoring data, which indicate that the soil within the shear zone is undergoing two deformation processes: a creep phase, characterised by different creep rates, and a dormant phase. We then study the creep behaviour of the soil samples through a series of ring shear creep tests under various shear stress conditions. The creep response depends strongly on the ratio of the shear stress to the residual strength, and the normal effective stress, whereas the creep rate decreases due to strength regain. The long-term strength of the clastic soil is close to the residual strength. Therefore, the residual strength obtained from conventional shear test, which is less time consuming than creep test, can be used in long-term stability analyses of creeping landslides.     

Research on the loss of group residential buildings under fierce winds

The prediction and assessment of the loss of group residential buildings under fierce winds are not only an important but also a very basic work. Many rural buildings, especially situated in southeastern of China, have suffered great damage and losses during such wind events in past 10 years. One new methodology is proposed in this study in order to estimate and evaluate the loss reasonably based on group building types and distribution, wind intensities and directions. This methodology comprehensively considers the mechanisms of interaction between the roof tiles, roof panels, doors and windows, which are more likely damaged seriously during fierce winds, the physical properties of those components, the wind pressure coefficients on the surface of group buildings in different wind directions, etc. Then Monte Carlo simulation is used to estimate the damage of group residential buildings under simulated fluctuating winds with different speeds and measured typhoon records. The simulation results indicate that our proposed quantified assessment method can be utilized for official, developers, architects, designers and homeowners to estimate and mitigate the losses that might be experienced during wind-related disasters either pre- or post-construction.     

Local site effect of a clay site in Shiraz based on seismic hazard of Shiraz Plain

Propagation of seismic waves through soil layers would drastically change the frequency content and amplitude-based features of ground motions at the surface. These alterations are known as seismic site effects. Computation of site effects of high-populated areas such as large cities is of great importance (e.g., it is used in development of seismic microzonation of a region). Shiraz is one of the most populous cities of Iran and is located in a high seismic hazardous region. A representative clay site in this city is selected to assess local site effects. The time series and random vibration theory procedure in the frequency domain are implemented to analyze the aforementioned site. Furthermore, the nonlinear dynamic soil behavior is simulated by the equivalent linear method and the nonlinear method via DEEPSOIL program. Three types of soil column uncertainties such as shear wave velocity, modulus reduction, and damping ratio of soil layers as well as depth of underlying rock half-space (D _bed) are considered herein. The mean amplification and standard deviation of natural logarithm of amplification factors are computed for a variety of analysis types. The results of the current study show that the computed mean and standard deviation of amplification factor in ln units by considering only V _S uncertainty are in good agreement with the corresponding ones by considering V _S and modulus reduction and damping ratio variabilities simultaneously for the studied site. Furthermore, it seems that the effect of bedrock depth in definition of spectral shapes of the Iranian seismic building code should be taken into account.