A multidirectional semi‐analytical method for analysis of laterally loaded pile groups in multi‐layered elastic strata

A semi‐analytical method for calculating the response of single piles and pile groups subjected to lateral loading is developed in this paper. Displacements anywhere in the soil domain are tied to the displacements of the piles through decay functions. The principle of virtual work and the calculus of variations are used to derive the governing differential equations that describe the response of the piles and soil. The eigenvalue method and the finite difference technique are used to solve the system of coupled differential equations for the piles and soil, respectively. The proposed method takes into account the soil surface displacement along and perpendicular to the loading direction and produces displacement fields that are very close to those produced by the finite element method but at lower computational effort. Compared with the previous method that considered only the soil displacement along the loading direction, accounting for the multi‐directional soil displacement field produces responses for the piles and soil that are closer to those approximated by the finite element method. Copyright © 2016 John Wiley & Sons, Ltd.     

Groundwater mixing between different aquifer types in a complex structural setting discerned by elemental and stable isotope geochemistry

Karst aquifers are well known for their intricate stratigraphy and geologic structures, which make groundwater characterization challenging because flowpaths and recharge sources are complex and difficult to evaluate. Geochemical data, collected from ten closely spaced production wells constructed in two karst aquifers (Bangor Limestone (Mb) and Tuscumbia Limestone/Fort Payne Chert (Mftp)) in Trussville, north‐central Alabama, illustrate two distinctive groundwater end‐members: (1) higher major ion, dissolved inorganic carbon, conductivity, alkalinity concentrations, heavier δ¹³C ratios (max: −10.2 ± 0.2‰ Vienna Pee Dee Belemnite (PDB)) and lower residence times (mean: 19.5 ± 2 years, n = 2) of groundwater in the Mb aquifer and (2) lower constituent concentrations, lighter δ¹³C ratios (min: −13.4 ± 0.2‰ PDB) and longer residence times of groundwater (mean: 23.6 ± 2 years, n = 4) in the Mftp aquifer. Summer and fall data and the binary mixing model show aquifer inter‐flow mixing along solution fractures and confirms the distinctive groundwater geochemistry of the two aquifers. Lowering of static water levels over the summer (drawdown from 2 to 5.2 m) leads to more reducing groundwater conditions (lower Eh values) and slightly enriched δ¹⁸O and δD ratios during the fall [δ¹⁸O: −4.8 ± 0.1 to −5.4 ± 0.1‰ Vienna Standard Mean Oceanic Water (VSMOW), n = 9; δD: −25.4 ± 1 to −27.4 ± 1‰ VSMOW, n = 9] when compared with summer season samples (δ¹⁸O: −5.1 ± 0.1 to −5.7 ± 0.1‰ VSMOW, n = 11; δD: −25.0 ± 1 to −30.6 ± 1‰ VSMOW, n = 11). GIS analyses confirm the localized origin of recharge to the investigated aquifers. The combination of GIS, field parameters and geochemistry analyses can be successfully used to identify recharge sources, evaluate groundwater flow and transport pathways and to improve understanding of how groundwater withdrawals impact the sustainability and susceptibility to contamination of karst aquifers. Copyright © 2015 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.     

Steady‐state analytical solutions of flow and deformation coupling due to a point sink within a finite fluid‐saturated poroelastic layer

An exact steady‐state closed‐form solution is presented for coupled flow and deformation of an axisymmetric isotropic homogeneous fluid‐saturated poroelastic layer with a finite radius due to a point sink. The hydromechanical behavior of the poroelastic layer is governed by Biot's consolidation theory. Boundary conditions on the lateral surface are specifically chosen to match the appropriate finite Hankel transforms and simplify the transforms of the governing equations. Ordinary differential equations in the transformed domain are solved, and then the analytical solutions in the physical space for the pore pressure and the displacements are finally obtained by using finite Hankel inversions. The analytical solutions at some special locations such as the top and bottom surfaces, lateral surface, and the symmetrical axis are given and analyzed. And a case study for the consolidation of a water‐saturated soft clay layer due to pumping is conducted. The analytical solution is verified against the finite element solution. Meanwhile, an analysis of coupled hydromechanical behavior is carried out herein. The presented analytical solution is an exact solution to the practical poroelastic problem within an axisymmetric finite layer. It can provide us a better understanding of the poroelastic behavior of the finite layer due to fluid extraction. Besides, it can be applied to calibrate numerical schemes of axisymmetric poroelasticity within finite domains. Copyright © 2017 John Wiley & Sons, Ltd.     

Optimal routing for waste collection: a case study in Singapore

ABSTRACT

To collect and dispose growing amounts of municipal solid waste (MSW) changed to be a hot topic along with the rapid urbanization in past decades. Cities are more and more dependent on the incineration instead of landfilling due to the cost-efficiency and environmental concerns. Considering the limited number of incineration plants and complicated situation of transportation in both spatial and temporal dimensions in different cities, the optimal routing for waste collection turns to be meaningful research topic. In this research, the ant colony optimization (ACO)-based multi-objective routing model coupled with min-max model and Dijkstra’s algorithm is proposed to address the question of which route to take from these waste-generating points to the target incineration plant(s) considering travel time, accident probability (black spots), and population exposure, so as to support the routing decision-making. The model is successfully implemented in Singapore and the effectiveness of the model has also been justified. Besides, few limitations of this research have also been discussed, some of which would also be the future directions of our research, especially the design and integration of a web-based routing decision-making support system.     

An Analytical Solution for Dynamic Response of the Plate with Different Impedances Subjected to Underwater Explosion

In this paper, the transmitted part of the incident wave is considered to revise Taylor’s solution, which is used to extend its application for analytical models to predict the response of the plate with different material properties. The influence of the material properties and the boundary condition of the plate on fluid and structural dynamics is systematically investigated. The analytical results are compared with those of detailed dynamic FE simulations and the two are in good agreement. The results indicate that the analytical method is valid and suitable for the plates with different material properties subjected to underwater explosion. It is found that Taylor’s results of the plate with small impedance are invalid, which indicates a potential application field for the analytical method.     

Effect of Groundwater Fluctuations on Pore Pressures and Earth Pressures on Coastal Excavation Retaining Walls

Pore water and earth pressures acting on retaining structures are investigated using an efficient coastal double-layered excavation model to determine offshore excavation responses to groundwater fluctuations outside foundation pits. Total pore water pressure includes excess pore water pressure (due to groundwater fluctuations) and steady pore water pressure (due to steady seepage) determined using one-dimensional consolidation theory of double-layered soil and one-dimensional steady-state flow theory, respectively. Rankine's active and passive earth pressures are obtained from pore water pressure. This method is applicable to arbitrary groundwater fluctuation conditions. How physical parameters affect pore water pressure is numerically investigated using examples, demonstrating the method's practicality for calculating pore water and earth pressures.     

Is Ziarat a Potential Site for Conventional or Unconventional Energy Resources?

A structural interpretation of the Ziarat block in the Balochistan region (a part of the Suleiman Fold and Thrust Belt) has been carried out using seismic and seismological data. Seismic data consists of nine 2.5D pre‐stack migrated seismic lines, whereas the seismological data covers the Fault Plane Solution and source parameters. Structural interpretation describes two broad fault sets of fore and back thrusts in the study area that have resulted in the development of pop‐up structures, accountable for the structural traps and seismicity pattern in terms of seismic hazard. Seismic interpretation includes time and depth contour maps of the Dungan Formation and Ranikot group, while seismological interpretation includes Fault Plane Solution, that is correlated with a geological and structural map of the area for the interpretation of the nature of the subsurface faults. Principal stresses are also estimated for the Ranikot group and Dungan Formation. In order to calculate anisotropic elastic properties, the parameters of the rock strength of the formations are first determined from seismic data, along with the dominant stresses (vertical, minimum horizontal, and maximum horizontal). The differential ratio of the maximum and minimum horizontal stresses is obtained to indicate optimal zones for hydraulic fracturing, and to assess the potential for geothermal energy reservoir prospect generation. The stress maps indicate high values towards the deeper part of the horizon, and low towards the shallower part, attributed to the lithological and structural variation in the area. Outcomes of structural interpretation indicate a good correlation of structure and tectonics from both seismological and seismic methods.     

成层软土地基预固结处理后桩基水平承载力估算方法

目前通过对软土地基预加固处理来提高桩基水平承载力已被工程界认可,但如何在工程前期设计过程中估算软土地基预处理后桩基水平承载力提高值仍是技术难点。基于此,参考Bowles[1]的地基土水平抗力计算式,同时考虑成层软土地基预排水固结处理影响,通过数学推导,推求出根据原状软土室内土工试验抗剪强度指标及预加固处理时间,估算软土地基预处理后桩基水平承载力提高值的实用计算方法。考虑桩侧土弹塑性屈服影响,推导出成层软土中水平受荷桩弹塑性解析解及塑性区深度的计算式,给出了桩顶水平位移、桩身最大弯矩的无量纲计算式及相关计算源代码。依托于浙江省某水闸桩基工程案例,根据提出的计算方法对桩基水平承载力、桩顶水平位移及桩身最大弯矩等性状进行预估计算,并与地基预处理前、后现场试桩检测值进行验证对比,认为桩基水平承载力、桩顶水平位移及桩身最大弯矩等预估计算成果与工程现场试桩的检测值较接近,对类似工程设计具有较好的参考价值。