Optimal and sustainable extraction of groundwater in coastal aquifers

Four examples are investigated for the optimal and sustainable extraction of groundwater from a coastal aquifer under the threat of seawater intrusion. The objectives and constraints of these management scenarios include maximizing the total volume of water pumped, maximizing the profit of selling water, minimizing the operational and water treatment costs, minimizing the salt concentration of the pumped water, and controlling the drawdown limits. The physical model is based on the density-dependent advective-dispersive solute transport model. Genetic algorithm is used as the optimization tool. The models are tested on a hypothetical confined aquifer with four pumping wells located at various depths. These solutions establish the feasibility of simulating various management scenarios under complex three-dimensional flow and transport processes in coastal aquifers for the optimal and sustainable use of groundwater.     

Collapsing loess-like soil in western Saudi Arabia

Earth fissures have developed at Wadi Al-Yutamah, western Saudi Arabia. The fissues are associated with land subsidence which is considered to be due to both rapid draw down of the ground-water level and hydrocompaction of the wadi soil after flooding. This phenomenon is relatively recent in the area.The wadi soil was investigated and classified in the field, and disturbed and undisturbed samples were collected for laboratory testing and analysis. The engineering properties of the wadi soils were determined, including in situ field density, specific gravity, liquid limit, plastic limit, plasticity index, shrinkage limit, and consolidation characteristics such as total settlement and coefficient of subsidence (collapse).Four types of silty soil of different colours were identified, in addition to yellowish-brown sand (representing the dunes) which accumulated on the surface of the study area. The silty soils include yellowish-brown silty sand, yellowish-gray silt, pale brown silt and yellowish-brown clayey silt. The results of X-ray diffraction analysis indicated that the clay is mainly kaolinite and illite with minor smectite. The dominant soil type in the study area is silt of low plasticity, high void ratio and low density which decreased with depth. This soil was classified as loess-like materials.The studied soils are of a collapsing type, and settlement is greatly increased by excessive wetting under constant pressure. The calculated coefficients of subsidence (collapse) of the soils at different depths generally increased with decreasing soil density and ranged between 3·1% and 10·8%. The wadi soils are considered to pose a moderate problem when wetted.     

DHI evaluation by combining rock physics simulation and statistical techniques for fluid identification of Cambrian-to-Cretaceous clastic reservoirs in Pakistan

The use of seismic direct hydrocarbon indicators is very common in exploration and reservoir development to minimise exploration risk and to optimise the location of production wells. DHIs can be enhanced using AVO methods to calculate seismic attributes that approximate relative elastic properties. In this study, we analyse the sensitivity to pore fluid changes of a range of elastic properties by combining rock physics studies and statistical techniques and determine which provide the best basis for DHIs. Gassmann fluid substitution is applied to the well log data and various elastic properties are evaluated by measuring the degree of separation that they achieve between gas sands and wet sands. The method has been applied successfully to well log data from proven reservoirs in three different siliciclastic environments of Cambrian, Jurassic, and Cretaceous ages. We have quantified the sensitivity of various elastic properties such as acoustic and extended elastic (EEI) impedances, elastic moduli (K _sat and K _sat–μ), lambda–mu–rho method (λρ and μρ), P-to-S-wave velocity ratio (V _P/V _S), and Poisson’s ratio (σ) at fully gas/water saturation scenarios. The results are strongly dependent on the local geological settings and our modeling demonstrates that for Cambrian and Cretaceous reservoirs, K _sat–μ, EEI, V _P/V _S, and σ are more sensitive to pore fluids (gas/water). For the Jurassic reservoir, the sensitivity of all elastic and seismic properties to pore fluid reduces due to high overburden pressure and the resultant low porosity. Fluid indicators are evaluated using two metrics: a fluid indicator coefficient based on a Gaussian model and an overlap coefficient which makes no assumptions about a distribution model. This study will provide a potential way to identify gas sand zones in future exploration.     

Stress distribution from railway track over geogrid reinforced ballast underlain by clay

Railway ballast forms a major component of a conventional rail track and is used to distribute the load to the subgrade, providing a smooth running surface for trains. It plays a significant role in providing support for the rail track base and distributing the load to the weaker layer underneath. Ballast also helps with drainage, which is an important factor for any type of transportation structure, including railroads. Over time, ballast progressively deforms and degrades under dynamic loading and loses its strength. In this study, extensive laboratory tests were conducted to investigate the effect of load amplitude, geogrid position, and number of geogrid layers, thickness of ballast layer and clay stiffness on the behavior of the reinforced ballast layer and induced strains in a geogrid. A half full-scale railway was constructed for carrying out the tests, which consisted of two rails 800 mm in length with three wooden sleepers(900 mm × 10 mm × 10 mm). Three ballast thicknesses of 200, 300 and 400 mm were used in the tests. The ballast was overlying 500 mm thickness clay in two states, soft and stiff. The tests were carried out with and without geogrid reinforcement; the tests were performed in a well-tied steel box of 1.5 m length ×1 m width ×1 m height. Laboratory tests were conducted to investigate the response of the ballast and the clay layers where the ballast was reinforced by a geogrid. Settlement in ballast and clay, soil pressure and pore water pressure induced in the clay were measured in reinforced and unreinforced ballast cases. It was concluded that the amount of settlement increased as the simulated train load amplitude increased, and there was a sharp increase in settlement up to cycle 500. After that, there was a gradual increase that leveled out between, 2500 to 4500 cycles depending on the frequency used. There was a slight increase in the induced settlement when the load amplitude increased from 0.5 to 1 ton but it was higher when the load amplitude increased to 2 tons. The increased amount in settlement depended on the existence of the geogrid and other parameters studied. The transmitted average vertical stress for ballast thicknesses of 30 cm and 40 cm increased as the load amplitude increased, regardless of the ballast reinforcement for both soft and stiff clay. The position of the geogrid had no significant effect on the transmitted stresses. The value of the soil pressure and pore water pressure on ballast thicknesses of 20 cm was higher than for 30 cm and 40 cm thicknesses. This meant that the ballast attenuated the induced waves. The soil pressure and pore water pressure for reinforced and unreinforced ballast was higher in stiff clay than in soft clay.     

Effect of near-fault earthquakes with forward directivity on telecommunication towers

In this paper, the effect of pulse-type motions caused by forward directivity that can release huge amounts of energy in a short time period is studied on a telecommunication tower. Since telecommunication towers have longer periods, they are not as affected by seismic forces. Nevertheless, near source earthquakes characterized by high velocity and velocity pulses can change the behavior of these structures. For this reason, a telecommunication tower located near active faults was selected in this study. Considering the probable earthquake magnitude at the site and the distance of the tower from adjacent faults, nine simulated pulses and three near-fault earthquake records with forward directivity are selected and applied to a 3D fi nite element model of the tower. The results of nonlinear dynamic analysis, i.e., displacements and damage in the tower, indicate that the maximum displacement and drift ratio of the tower under the pulses are obviously affected by the ratio of the structure period to pulse period. When this ratio is decreased and close to 1.0, the maximum displacement and drift ratio are sharply increased and cause large displacements in the tower.     

Hydrological responses to climatic variability in a cold agricultural region

Extended severe dry and wet periods are frequently observed in the northern continental climate of the Canadian Prairies. Prairie streamflow is mainly driven by spring snowmelt of the winter snowpack, whilst summer rainfall is an important control on evapotranspiration and thus seasonality affects the hydrological response to drought and wet periods in complex ways. A field‐tested physically based model was used to investigate the influences of climatic variability on hydrological processes in this region. The model was set up to resolve agricultural fields and to include key cold regions processes. It was parameterized from local and regional measurements without calibration and run for the South Tobacco Creek basin in southern Manitoba, Canada. The model was tested against snow depth and streamflow observations at multiple scales and performed well enough to explore the impacts of wet and dry periods on hydrological processes governing the basin scale hydrological response. Four hydro‐climatic patterns with distinctive climatic seasonality and runoff responses were identified from differing combinations of wet/dry winter and summer seasons. Water balance analyses of these patterns identified substantive multiyear subsurface soil moisture storage depletion during drought (2001–2005) and recharge during a subsequent wet period (2009–2011). The fractional percentage of heavy rainfall days was a useful metric to explain the contrasting runoff volumes between dry and wet summers. Finally, a comparison of modeling approaches highlights the importance of antecedent fall soil moisture, ice lens formation during the snowmelt period, and peak snow water equivalent in simulating snowmelt runoff.     

Optimization of nonconventional wells under uncertainty using statistical proxies

The determination of the optimal type and placement of a nonconventional well in a heterogeneous reservoir represents a challenging optimization problem. This determination is significantly more complicated if uncertainty in the reservoir geology is included in the optimization. In this study, a genetic algorithm is applied to optimize the deployment of nonconventional wells. Geological uncertainty is accounted for by optimizing over multiple reservoir models (realizations) subject to a prescribed risk attitude. To reduce the excessive computational requirements of the base method, a new statistical proxy (which provides fast estimates of the objective function) based on cluster analysis is introduced into the optimization process. This proxy provides an estimate of the cumulative distribution function (CDF) of the scenario performance, which enables the quantification of proxy uncertainty. Knowledge of the proxy-based performance estimate in conjunction with the proxy CDF enables the systematic selection of the most appropriate scenarios for full simulation. Application of the overall method for the optimization of monobore and dual-lateral well placement demonstrates the performance of the hybrid optimization procedure. Specifically, it is shown that by simulating only 10% or 20% of the scenarios (as determined by application of the proxy), optimization results very close to those achieved by simulating all cases are obtained.     

A middle–late Ediacaran volcano‐sedimentary record from the eastern Arabian‐Nubian shield

The Ediacaran Jibalah Group comprises volcano‐sedimentary successions that filled small fault‐bound basins along the NW–SE‐trending Najd fault system in the eastern Arabian‐Nubian Shield. Like several other Jibalah basins, the Antaq basin contains exquisitely preserved sedimentary structures and felsic tuffs, and hence is an excellent candidate for calibrating late Ediacaran Earth history. Shallow‐marine strata from the upper Jibalah Group (Muraykhah Formation) contain a diversity of load structures and intimately related textured organic (microbial) surfaces, along with a fragment of a structure closely resembling an Ediacaran frond fossil and a possible specimen of Aspidella. Interspersed carbonate beds through the Muraykhah Formation record a positive δ¹³C shift from ?6 to 0‰. U‐Pb zircon geochronology indicates a maximum depositional age of ~570 Ma for the upper Jibalah Group, consistent with previous age estimates. Although this age overlaps with that of the upper Huqf Supergroup in nearby Oman, these sequences were deposited in contrasting tectonic settings on opposite sides of the final suture of the East African Orogen.     

System identification of instrumented bridge systems

Several recorded motions for seven bridge systems in California during recent earthquakes were analysed using parametric and non‐parametric system identification (SI) methods. The bridges were selected considering the availability of an adequate array of accelerometers and accounting for different structural systems, materials, geometry and soil types. The results of the application of SI methods included identification of modal frequencies and damping ratios. Excellent fits of the recorded motion in the time domain were obtained using parametric methods. The multi‐input/single‐output SI method was a suitable approach considering the instrumentation layout for these bridges. Use of the constructed linear filters for prediction purposes was also demonstrated for three bridge systems. Reasonable prediction results were obtained considering the various limitations of the procedure. Finally, the study was concluded by identifying the change of the modal frequencies and damping of a particular bridge system in time using recursive filters. Copyright © 2003 John Wiley & Sons, Ltd.