Radial basis function neural network for hydrologic inversion: an appraisal with classical and spatio-temporal geostatistical techniques in the context of site characterization

This paper investigates three techniques for spatial mapping and the consequential hydrologic inversion, using hydraulic conductivity (or transmissivity) and hydraulic head as the geophysical parameters of concern. The data for the study were obtained from the Waste Isolation and Pilot Plant (WIPP) site and surrounding area in the remote Chihuahuan Desert of southeastern New Mexico. The central technique was the Radial Basis Function algorithm for an Artificial Neural Network (RBF-ANN). An appraisal of its performance in light of classical and temporal geostatistical techniques is presented. Our classical geostatistical technique of concern was Ordinary Kriging (OK), while the method of Bayesian Maximum Entropy (BME) constituted an advanced, spatio-temporal mapping technique. A fusion technique for soft or inter-dependent data was developed in this study for use with the neural network. It was observed that the RBF-ANN is capable of hydrologic inversion for transmissivity estimation with features remaining essentially similar to that obtained from kriging. The BME technique, on the other hand, was found to reveal an ability to map localized lows and highs that were otherwise not as apparent in OK or RBF-ANN techniques.     

Groundwater management scenarios for the Biyadh-Wasia aquifer systems in the eastern part of Riyadh region,Saudi Arabia

Groundwater stored in the deep seated sedimentary aquifers is the most important source of water supply. The lack of sufficient groundwater recharge and the overdependence on groundwater might lead to unavailability of this precious natural resource if proper management practices are not adopted. Finite difference modelling using the MODFLOW program was carried out in the east of Riyadh city to simulate the groundwater level conditions under different abstraction scenarios. The simulated aquifer system combines the Wasia and Biyadh aquifer (composed of sandstone) and Aruma aquifer (limestone) which lies between 24°30'00"- 25°30'00" N and longitudes 47°00'00"- 48°00'00"E. The transmissivity and storage coefficient values of Biyadh aquifer are 7.0x10^-3 to 7.0x10^-2 m²/day and 3.7x10^-4 to 9.4x10^-4 respectively. The transmissivity and storage coefficient values of Wasia aquifer ranges from 6.7x10^-3 to 8.5x10^-2 m²/day and 2x10^-4 to 2.3x10^-4 respectively. The model calibration involved altering the values of model input parameters to match field conditions within certain acceptable limits to forecast the aquifer response over a period of 35 years (2015-2050). The modelling grid consisted of 20 and 24 columns with the grids spacing of 4 km for the small grids and 6 km for large grids. The results showed that though the Wasia aquifer was productive, it showed a large decline in water levels if water abstraction continued at the present rate. If the existing trends of groundwater withdraw continues; the piezometric heads in Wasia and Biyadh aquifers will decline by the year 2050. A reduction in 25% of the existing groundwater pumping rate in the well field will minimize the rate of groundwater decline in the aquifer to a considerable extent.     

Environmental assessment of soil,groundwater, and surface water quality in the south of the Riyadh,Saudi Arabia

To understand the impact of the influence of treated wastewater, a study was undertaken in the downstream side of Wadi Hanifa in the southern part of Riyadh City in Saudi Arabia. Chemical elements from 17 groundwater samples, 9 surface water samples (treated wastewater), and 14 soil samples were analyzed. Water facies analyses showed that both groundwater and surface water belong to the SO₄-Cl class. However, the groundwater is characterized by high salinity (average value of 3547 mg/l), which may be result of the greater rock-water interaction and limited rainfall recharge. The NO₃ values are also high in the groundwater samples (average value of 40 mg/l) and are mainly attributed to the agricultural practices in the study area. The surface water samples (treated wastewater from the plant) shows an average salinity value of 1232 mg/l and is better suited for irrigation purposes. Heavy elements analyzed in the soil samples show high concentrations of all the elements except Mn and Ni as compared to their background concentration. Enrichment factor and index of geoaccumulation calculated from the soil samples indicate that they are severely enriched with respect to Cd and Se. The spatial distribution maps were prepared based on kriging interpolation technique to estimate the concentrations of the analyzed elements at unknown locations. The treated wastewater in the study area is better suited for agricultural and domestic activities as compared to groundwater.     

Steered migration in hard rock environments

Hard rock seismic exploration normally has to deal with rather complex geological environments. These types of environments are usually characterized by a large number of local heterogeneity (e.g., faults, fracture zones, and steeply dipping interfaces). The seismic data from such environments often have a poor signal‐to‐noise ratio because of the complexity of hard rock geology. To be able to obtain reliable images of subsurface structures in such geological conditions, processing algorithms that are capable of handling seismic data with a low signal‐to‐noise ratio are required for a reflection seismic exploration. In this paper, we describe a modification of the 3D Kirchhoff post‐stack migration algorithm that utilizes coherency attributes obtained by the diffraction imaging algorithm in 3D to steer the main Kirchhoff summation. The application to a 3D synthetic model shows the stability of the presented steered migration to the presence of high level of the random noise. A test on the 3D seismic volume, acquired on a mine site located in Western Australia, reveals the capability of the approach to image steep and sharp objects such as fracture and fault zones and lateral heterogeneity.     

Group velocity dispersion analysis in northern Peninsular Malaysia

The crustal structure beneath three seismic stations over Malaysia has been investigated with the application of the group velocity dispersion analysis of the northern Sumatra earthquake data which occurred on 06 April 2010. Eighteen crustal layer models are constructed to assess the structure. Group velocity dispersions have been computed for the recorded earthquake data using a graphical method and modified Haskell matrix method for the models. Both dispersions have been presented for the interpretation of crustal layers. Findings have shown four major crustal layers having thicknesses of 2.5–4.0, 2.0–5.5, 5.0–8.0, and 8.5–9.0 km, while in Terengganu, it has shown three layers. Density, shear, and compressional wave velocities used in models have suggested that the crustal structure of the northern part of Peninsular Malaysia is crystalline. Major crustal minerals are of quartz, plagioclase, and mica. Most layers seem to have upward directions toward Perak from Kedah and Terengganu.     

Is correlation dimension a reliable proxy for the number of dominant influencing variables for modeling risk of arsenic contamination in groundwater?

The correlation dimension (CD) of a time series provides information on the number of dominant variables present in the evolution of the underlying system dynamics. In this study, we explore, using logistic regression (LR), possible physical connections between the CD and the mathematical modeling of risk of arsenic contamination in groundwater. Our database comprises a large-scale arsenic survey conducted in Bangladesh. Following the recommendation by Hossain and Sivakumar (Stoch Environ Res Risk Assess 20(1–2):66–76, 2006), who reported CD values ranging from 8 to 11 for this database, 11 variables are considered herein as indicators of the aquifer’s geochemical regime with potential influence on the arsenic concentration in groundwater. A total of 2,048 possible combinations of influencing variables are considered as candidate LR risk models to delineate the impact of the number of variables on the prediction accuracy of the model. We find that the uncertainty associated with prediction of wells as safe and unsafe by LR risk model declines systematically as the total number of influencing variables increases from 7 to 11. The sensitivity of the mean predictive performance also increases noticeably for this range. The consistent reduction in predictive uncertainty coupled with the increased sensitivity of the mean predictive behavior within the universal sample space exemplify the ability of CD to function as a proxy for the number of dominant influencing variables. Such a rapid proxy, based on non-linear dynamic concepts, appears to have considerable merit for application in current management strategies on arsenic contamination in developing countries, where both time and resources are very limited.     

Numerical simulation of fabric anisotropy and strain localization of sand under simple shear

This paper studies the effects of initial fabric anisotropy of dry sand in simple shear deformation. The effects of anisotropy are taken into consideration through the modification of the mobilized friction in the Mohr–Coulomb‐type yield surface as a function of a fabric parameter. In addition, the constitutive model uses a gradient term that directly incorporates the effects of material length scale. The constitutive formulation is implemented into ABAQUS finite element code and used to simulate shearing of the dry sand under various conditions of simple shear. The numerical simulations show that while the shear stress response is affected by fabric anisotropy, its effects on strain localization in simple shear are minimal. This is in contrast to other devices such as the biaxial shear. The strain localization in simple shear is controlled more by the imposed boundary conditions. The use of material length scale is shown to remove the effects of strain localization in the shearing response. Copyright © 2009 John Wiley & Sons, Ltd.     

Origin of Quartz Cement in the Paleogene Sandstone Reservoir of Shahejie Formation, Bohai Bay Basin, China

The precipitation of authigenic quartz plays a significant role to reduce the reservoir characteristics and enhance the stiffness of the rock.The Es₁ sandstone of Shahejie Formation is acting as a significant hydrocarbon producing rock in the Nanpu Sag.Various methods like thin section petrography,cathodoluminescence(CL),scanning electron microscope(SEM,with EDS),and electron microprobe analysis has been used to reveal the origin of quartz cement as well as to evaluate the effect of quartz cement on reservoir quality.The studied sandstone is classified as immature to mature feldspathic litharenite and lithic arkose and consists of quartz,feldspar,rock fragments and micas.Petrographic studies and SEM analysis shows that the authigenic quartz is acting a significant cement that reduces the reservoir quality.Whereas clay minerals(kaolinite and mixed layer illite to smectite)are dominant in the Es₁ sandstone,that can reduce the reservoir quality.SEM,CL and thin section analysis reveal that there are two stages of quartz cement in the studied samples;that are pore filling authigenic cement and quartz overgrowth cement.Fluid inclusion homogenization temperatures shows that stages of quartz cement were developed with continuous process from 70℃ to 130℃.Quartz cements were generally originated from I/S reaction,feldspar dissolution,conversion of rock fragments and pressure solution.Feldspar dissolution(K-feldspar)and kaolinite to illite reaction is an insignificant silica source for the silica cement which is internally precipitated in a close system with diffusion transporting mechanism.Overall,quartz cement significantly enhance the rock strengthen and brittleness effectively as well as it reduce the overall reservoir quality.     

Transport phenomena in kaolinite clay: Molecular simulation,homogenization analysis and similitude law

Kaolinite is a common clay mineral. It is a nanomaterial with a platelet crystalline structure. In order to analyze the behavior of kaolinite, its microscopic structure and material properties must be specified correctly. A molecular dynamics (MD) simulation is used for determining the microscale properties of hydrated kaolinite, and these properties are introduced into a multiscale homogenization analysis (HA). We previously developed such an MD/HA technique to investigate seepage, diffusion, sorption and consolidation in bentonite clay (Proceedings of the Science Basis for Nuclear Waste Management, Davos, Switzerland, vol. XXI. Material Research Society: Warrendale, PA, 1997; 359–366; Eng. Geol. 1999; 54 :21–31; Eng. Geol. 2001; 60 :127–138; Coupled Thermo‐Hydro‐Mechanical‐Chemical Processes in Geo‐systems. Elsevier: Amsterdam, 2005; 457–464). We here apply the method to kaolinite clay to investigate the permeability, diffusion and related similitude law. The obtained results are supported by existing experimental data. Copyright © 2008 John Wiley & Sons, Ltd.