Hydrogeological assessment of the Tabarteh fault zone by physicochemical analysis and multivariate statistical methods

Hydrogeologically, faults may impede, conduit, exert no influence, or may play a combination of these roles on groundwater flow. The object of this paper is to study the hydrogeological role of the Tabarteh fault, which is located on the border of Zagros and Central Iran tectonic zones in an alluvial aquifer. The recorded data of water table levels, chemical parameters, and discharge rate of wells, in addition to geological maps and geophysical results, were collected and evaluated. The outcrop of travertine in limited areas and the emergence of a few small springs within the alluvium show a barrier role of the fault in the groundwater flow. The spatial analysis of chemical components, head time series, and groundwater flow direction assessment demonstrated that the fault acts as both a barrier and a non-barrier in different sections. The multivariate statistical methods of cluster and discriminant analyses also confirm the dual role of the fault.     

New occurrence of potential phosphate resource in northeast Jordan

An exploration program aided by field investigation, exploration drilling, detailed sampling, lithological and petrological studies, geochemical investigation, and resource calculation leads to the discovery of a potential phosphate resource in northeast Jordan close to the international border with Saudi Arabia and Iraq. The studied phosphate is of the Middle Eocene age that belongs to the Wadi Shallala Formation. It is equivalent to the phosphate deposits recorded in the lower part of the Umm Wual Formation in the Turayf region of Saudi Arabia and the Eocene Ratga Formation in the Ethna phosphate deposit west of Iraq. The phosphorites in the region are broadly similar in mineralogical composition and geochemical affinities. X-ray diffraction indicates the presence of francolite with variable amounts of calcite and quartz. Most samples consist of phosphate clasts embedded in carbonate and silica matrix and cement. P₂O₅ content is up to 32.3 % with an average equal to 18.6 %. The impurity is caused by the presence of variable amounts of SiO₂ and CaO. The F% and F/P₂O₅ ratio in studied phosphates is lower compared with that in phosphates from Jordan and Saudi Arabia. The geological and geochemical results were integrated for resource estimation. Three high-grade phosphate layers with ≥23 % P₂O₅ were considered in the calculations. The phosphate resource is classified as an inferred resource. The total volume of the resource is about 649 million tons. The average P₂O₅ content is 24.57, and the stripping ratio is 1:5.8.     

Textural evidences for dissolution of silica-rich rocks of the Ypresian phosphatic series,Gafsa-Metlaoui basin,southwestern Tunisia: implication of biogenic silica supply on genesis of fibrous clays

The present work aimed to determine the mineralogical composition of Ypresian series and to clarify the influence of the dissolution of siliceous frustules on the genesis of fibrous clay minerals. The specimens sampled from CPG trench are mainly constituted of silica-rich rocks at Mides area located at the western part of Gafsa-Metlaoui basin. The samples were investigated using X-ray diffraction (XRD) and scanning electron microscopy (SEM) in order to determine texture of constituents. The data obtained indicate that the bulk rock samples are mainly made up of opal CT and clay minerals. The latter consists of palygorskite-sepiolite minerals associated with smectite and few amount of illite. The trend of mineralogical composition indicates that fibrous clay minerals are more concentrated at the upper part. SEM observations indicated that palygorskite mineral appears as thread-like facies, which surround foliated texture of smectite in the lower part of the Mides section, although with the low Mg activity confirmed by the absence of dolomite. But, at the upper part of the Mides section, SEM observations revealed the occurrence of siliceous frustules, which have numerous dissolved areas and replacement of carbonate tests by silica. The dissolution saturated the depositional environment with silica which is required for the formation of palygorskite and sepiolite minerals, in addition to high Mg activity confirmed by the presence of dolomite in the bulk rock, which is required basically for the formation of sepiolite. Although the genesis mode of palygorskite and sepiolite is similar with very little difference, the genesis of sepiolite needs a high alkalinity than the formation of palygorskite.     

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.     

Prediction of water inflow to mechanized tunnels during tunnel-boring-machine advance using numerical simulation

An accurate estimate of the groundwater inflow to a tunnel is one of the most challenging but essential tasks in tunnel design and construction. Most of the numerical or analytical solutions that have been developed ignore tunnel seepage conditions, material properties and hydraulic-head changes along the tunnel route during the excavation process, leading to inaccurate prediction of inflow rates. A method is introduced that uses MODFLOW code of GMS software to predict inflow rate as the tunnel boring machine (TBM) gradually advances. In this method, the tunnel boundary condition is conceptualized and defined using Drain package, which is simulated by dividing the drilling process into a series of successive intervals based on the tunnel excavation rates. In addition, the drain elevations are specified as the respective tunnel elevations, and the conductance parameters are assigned to intervals, depending on the TBM type and the tunnel seepage condition. The Qomroud water conveyance tunnel, located in Lorestan province of Iran, is 36 km in length. Since the Qomroud tunnel involved groundwater inrush during excavating, it is considered as a good case study to evaluate the presented method. The groundwater inflow to this tunnel during the TBM advance is simulated using the proposed method and the predicted rates are compared with observed rates. The results show that the presented method can satisfactorily predict the inflow rates as the TBM advances.     

Subduction-related mafic to felsic magmatism in the Malayer–Boroujerd plutonic complex,western Iran

The Malayer–Boroujerd plutonic complex (MBPC) in western Iran, consists of a portion of a magmatic arc built by the northeast verging subduction of the Neo-Tethys plate beneath the Central Iranian Microcontinent (CIMC). Middle Jurassic-aged felsic magmatic activity in MBPC is manifested by I-type and S-type granites. The mafic rocks include gabbroic intrusions and dykes and intermediate rocks are dioritic dykes and minor intrusions, as well as mafic microgranular enclaves (MMEs). MBPC Jurassic-aged rocks exhibit arc-like geochemical signatures, as they are LILE- and LREE-enriched and HFSE- and HREE-depleted and display negative Nb–Ta anomalies. The gabbro dykes and intrusions originated from metasomatically enriched garnet-spinel lherzolite [Degree of melting (f_mel) ~ 15%] and exhibit negative Nd and positive to slightly negative ε_Hf(T) (+ 3.0 to ? 1.6). The data reveal that evolution of Middle Jurassic magmatism occurred in two stages: (1) deep mantle-crust interplay zone and (2) the shallow level upper crustal magma chamber. The geochemical and isotopic data, as well as trace element modeling, indicate the parent magma for the MBPC S-type granites are products of upper crustal greywacke (f_mel: 0.2), while I-type granites formed by partial melting of amphibolitic lower crust (f_mel: 0.25) and mixing with upper crustal greywacke melt in a shallow level magma chamber [Degree of mixing (f_mix): 0.3]. Mixing between andesitic melt leaving behind a refractory dense cumulates during partial crystallization of mantle-derived magma and lower crustal partial melt most likely produced MMEs (f_mix: 0.2). However, enriched and moderately variable ε_Nd(T) (? 3.21 to ? 4.33) and high (⁸⁷Sr/⁸⁶Sr)_i (0.7085–0.7092) in dioritic intrusions indicate that these magmas are likely experienced assimilation of upper crustal materials. The interpretations of magmatic activity in the MBPC is consistent with the role considered for mantle-derived magma as heat and mass supplier for initiation and evolution of magmatism in continental arc setting, elsewhere.