A spatio-temporal three-dimensional conceptualization and simulation of Dera Ismail Khan alluvial aquifer in visual MODFLOW: a case study from Pakistan

Dera Ismail Khan (DIK) is situated in the Lower Indus Basin of Pakistan. The land use has been changed in the canal command area due to irrigation activities near the Indus River. To check the current status and predict the groundwater levels in the area, the unconfined aquifer has been simulated in Visual MODFLOW for a period of 35 years, i.e., from 1985 to 2020. The 2900-km² area has been modeled with a grid of 500 by 500 m and the depth set to 100 m. The aquifer in the study area has been divided vertically and laterally into three and ten zones, respectively, for the characterization. Water wells and streams were used as the sinks and hydrologic boundaries, respectively. The model was successfully calibrated in steady and the non-steady state. The simulation revealed that the whole simulation can be divided into two phases, i.e., before and after the construction of the Chashma Right Bank Canal (CRBC), whereas the results were summarized in the form of water table depth maps and groundwater budget calculations. To determine the groundwater sustainability, a conjunctive use scenario has been employed to simulate the aquifer dynamics till 2020. The simulation revealed incremental drawdowns till the end.     

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.     

A method of alternating characteristics with application to advection-dominated environmental systems

We present a numerical method for solving a class of systems of partial differential equations (PDEs) that arises in modeling environmental processes undergoing advection and biogeochemical reactions. The salient feature of these PDEs is that all partial derivatives appear in linear expressions. As a result, the system can be viewed as a set of ordinary differential equations (ODEs), albeit each one along a different characteristic. The method then consists of alternating between equations and integrating each one step-wise along its own characteristic, thus creating a customized grid on which solutions are computed. Since the solutions of such PDEs are generally smoother along their characteristics, the method offers the potential of using larger time steps while maintaining accuracy and reducing numerical dispersion. The advantages in efficiency and accuracy of the proposed method are demonstrated in two illustrative examples that simulate depth-resolved reactive transport and soil carbon cycling.     

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.     

Health risk of heavy metals from vegetables irrigated with sewage water in peri-urban of Dera Ismail Khan,Pakistan

Health risks of heavy metals in vegetables irrigated with sewage water were investigated in the present study. The findings indicated a massive accumulation of heavy metals in soil and vegetables collected from Dera Ismail khan, Pakistan. The concentration levels of heavy metal in vegetables grown on soil irrigated with untreated sewage water were significantly higher at (P ≤ 0.001) than in vegetables grown on fresh-water-irrigated soil and proceeded the recommended limits of World health organization. Moreover, the findings also indicated that the adults and children consuming such vegetables ingested a large proportion of the selected toxic metals. Health risk index was greater than one for Pb and Cd in all the selected vegetables and was greater than one for Ni in three vegetables like Spinacia oleracea, Benincasa fistulosa and Lactuca sativa. Health risk assessment would be a useful tool for information regarding any threats of heavy metals contamination in vegetables.     

High-performing antifouling bacterial consortium for submerged membrane bioreactor treating synthetic wastewater

Biofouling, associated with membranes, is considered as a major operational challenge in membrane bioreactor (MBR) technology. Interrupting the process for the formation of biofilm by the action of interspecies quorum quenching (QQ) has received a significant attention since recent years. An antifouling bacterial consortium was identified to improve biofouling inhibition performance during MBR operation. For this purpose, various QQ bacteria were isolated from laboratory-scale MBR using enrichment culture method and identified via 16S rRNA. Potential quenching strains including Enterobacter cloaca, Delftia sp., and Pseudomonas sp. were utilized to control biofouling in the MBR operated in the continuous mode for 38 days. Three laboratory-scale MBRs, including two MBRs with different anti fouling consortium and a control, were operated in parallel under similar operating conditions. Biofouling control by QQ bacteria was compared based on the membrane permeability and EPS secretion from biofilm on the membrane. Both MBRs with antifouling consortium (AC-MBRs) experienced around three times less biofouling as compared to conventional MBR leading to significant decrease in acyl homoserine lactones (AHLs) concentration in the biocake. More than 90, 45, and 49% of COD, NH₄–N, and PO₄ ^3?–P removal efficiencies elucidate that QQ bacterial consortium could effectively reduce membrane biofouling without compromising the MBR efficiency. Comparatively lower concentration of bound EPS in AC-MBRs restricted the bacterial adhesion to membrane resulting in enhanced membrane permeability depicting that a broader range of signal molecules could be hydrolyzed using antifouling consortium than single or no QQ strain in the submerged MBR.