Spatial multi-criteria analysis for selecting potential sites for aquifer recharge via harvesting and infiltration of surface runoff in north Jordan

In (semi-)arid regions, available water resources are scarce and groundwater resources are often overused. Therefore, the option to increase available water resources by managed aquifer recharge (MAR) via infiltration of captured surface runoff was investigated for two basins in northern Jordan. This study evaluated the general suitability of catchments to generate sufficient runoff and tried to identify promising sites to harvest and infiltrate the runoff into the aquifer for later recovery. Large sets of available data were used to create regional thematic maps, which were then combined to constraint maps using Boolean logic and to create suitability maps using weighted linear combination. This approach might serve as a blueprint which could be adapted and applied to similar regions. The evaluation showed that non-committed source water availability is the most restricting factor for successful water harvesting in regions with <200 mm/a rainfall. Experiences with existing structures showed that sediment loads of runoff are high. Therefore, the effectiveness of any existing MAR scheme will decrease rapidly to the point where it results in an overall negative impact due to increased evaporation if maintenance is not undertaken. It is recommended to improve system operation and maintenance, as well as monitoring, in order to achieve a better and constant effectiveness of the infiltration activities.     

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.     

Seismic stratigraphy of the Mianwali and Bannu depressions,north-western Indus foreland basin

International Journal of Earth Sciences - Regional seismic reflection profiles, deep exploratory wells, and outcrop data have been used to study the structure and stratigraphic architecture of the...     

Numerical solution and convergence analysis of steam injection in heavy oil reservoirs

In this paper, the numerical methods for solving the problem of steam injection in the heavy oil reservoirs are presented. We consider a 3-dimensional model of 3-phase flow, oil, water, and steam, with the effect of 3-phase relative permeability. Interphase mass transfer of water and steam is considered; oil is assumed nonvolatile. We apply the simultaneous solution approach to solve the corresponding nonlinear discretized partial differential equation in the fully implicit form. The convergence of finite difference scheme is proved by the Rosinger theorem. The heuristic Jacobian-Free-Newton-Krylov (HJFNK) method is proposed for solving the system of algebraic equations. The result of this proposed numerical method is well compared with some experimental results. Our numerical results show that the first iteration of the full approximation scheme (FAS) provides a good initial guess for the Newton method. Therefore, we propose a new hybrid-FAS-HJFNK method while there is no steam in the reservoir. The numerical results show that the hybrid-FAS-HJFNK method converges faster than the HJFNK method.     

Development of an efficient surrogate model based on aquifer dimensions to prevent seawater intrusion in anisotropic coastal aquifers,case study: the Qom aquifer in Iran

Coastal aquifers are usually exposed to saltwater intrusion. Therefore, groundwater extracted from these aquifers should be regulated considering their dimensions and effective parameters. In this paper, optimum discharge from a large number of exploitation wells is evaluated according to variations of width, length, and anisotropy coefficient in the Qom aquifer near the salt lake in central Iran. First, the wells are divided into clusters to decrease the number of decision variables. Then, the location and discharge from each cluster is obtained using SEAWAT and charged system search (CSS) simulation–optimization model with the assumption of three-dimensional variable density flow. The maximum discharge considering various anisotropy rates is computed based on different values of lengths and widths of the aquifer. Finally, an M5-tree model is trained using the obtained samples to derive a linear relationship between input and output data. Based on the results, for various ranges of width and length of an aquifer with impermeable boundaries, different linear equations for optimum discharge are obtained. Also, it was found that for an aquifer with a small width, the critical discharge is a function of the length while the effect of the boundaries is negligible. Sensitivity analysis of the anisotropy coefficient reveals that with increasing the anisotropy rate, thickness and slope of the transition zone decrease and as the maximum discharge increases consequently. However, the sensitivity of the discharge to anisotropy rate is not remarkable. A comparison between the results of this study with those of the analytical method based on sharp interface assumption is carried out. For the critical condition, the best agreement between analytical equation (\(\overline {L} =0.87\overline {W} +0.62\)) and proposed method (\(\overline {L} =0.83\overline {W} - 1.41\)) is achieved for the anisotropic aquifer when the 50% isochlor is assumed as the measure of salt water intrusion.     

Structural properties of Urfa stone doped with erbium oxide

In this study, we investigated the structural properties of Urfa stone (US) doped with erbium oxide (Er₂O₃). Solid US was powdered by using an agate mortar, and its elemental composition was determined using inductive coupling plasma (ICP) methods. Varying amounts of Er₂O₃ (5, 10, 20, 30, and 40%) were added as a dopant to the US powder using mechanical alloying methods. The resultant samples were sintered at 1000 °C for 1 h. The structural properties of the Er₂O₃-doped US samples were subsequently investigated using X-ray diffraction (XRD), Fourier-transform infrared spectrometry (FTIR), and photoluminescence methods. Results from the XRD analysis of the Er₂O₃-doped US powder indicated two crystalline phases: (1) calcium oxide (CaO) or lime and (2) Er₂O₃. After the samples were sintered at 1000 °C, CaO, Er₂O₃, calcium carbonate (CaCO₃), and mixed crystalline phases were observed. Results from the FTIR analysis of the Er₂O₃-doped US samples indicated absorption bands at 711.91, 872.08, and 1396.87 cm^?1 in the spectra. Finally, photoluminescence analysis results indicated a shift in the emission and excitation bands to longer and shorter wavelengths, respectively, in the solid state (non-aqueous media) US-Er complex.