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.     

The measurement of soil gases and shallow temperature for determination of active faults in a geothermal area: a case study from Ömer–Gecek,Afyonkarahisar (West Anatolia)

Afyonkarahisar is a very important geothermal province of western Anatolia and has low and medium enthalpy geothermal areas. This study has been carried out for the preparation of distribution maps of soil gases (radon and carbon dioxide) and shallow soil temperature and the exploration of permeable tectonic regions associated with geothermal systems and reveal the origins of radon and carbon dioxide gases. The western district of the study area is characterized by the high radon concentration (168.30 kBq/m³), carbon dioxide ratio (0.30%), and soil temperature (21.0 °C) values. Fethibey and Demirçevre faults, which allow the circulation of geothermal fluids, have been detected in the distribution maps of radon, carbon dioxide, and shallow depth temperature and the directions of the curves in these maps correspond to the strikes of Demirçevre faults. The effect of the fault plays an important role in the change of carbon dioxide concentration along the W-E directional geological section prepared to determine the change of soil gas and shallow depth temperature values depending on lithological differences, fault existence, and geothermal reservoir depth. On the other hand, it was determined that Rn²²² concentration and soil temperature changed as a function of geothermal reservoir depth or lithological difference. Tuffs in Köprülü volcano-sedimentary units are the main source of radon due to their higher uranium contents. Besides, the carbon dioxide in Ömer–Gecek soils has geothermal origin because of the highest carbon dioxide content (99.3%) in non-condense gas. The similarities in patterns of soil temperature, radon, and carbon dioxide indicate that the variation in soil temperatures is related to radon and carbon dioxide emissions. It is concluded that soil gas and temperature measurements can be used to determine the active faults in the initial stage of geothermal exploration successfully.     

Uncertainties in Grout Penetrability Measurements; Evaluation and Comparison of Filter pump,Penetrability Meter and Short Slot

To measure grout penetrability in fractured hard rock, various measuring instruments have been developed over the years. Penetrability meter and Filter pump have been designed to use in both the lab and the field. Short slot has been applicable mainly in the lab due to its complexity. The fact, that these instruments have been built based on different assumptions, limitations, and test conditions, makes their results occasionally in contradict. Deficiency in design of the instruments as well as the methods of evaluating grout penetrability is additionally a basis for uncertainty in results. This study is an experimental effort to determine and thoroughly perceive the nature of the most governing uncertainties in grout penetrability measurements. The test apparatus, procedure, and method used to evaluate the grout penetrability in both Penetrability meter and Filter pump were thus modified. The aim was to control the corresponding uncertainties and make their limitations and test conditions as similar as possible with the ones in Short slot. The results suggested that to obtain a more realistic evaluation of the grout penetrability, measurement should be accomplished at both the high and the low pressures with sufficient grout volume using Short slot. Moreover, application of both Filter pump and Penetrability meter is no longer recommended due to the revealed uncertainties.     

2-D delineation of individual citrus trees from UAV-based dense photogrammetric surface models

One of the challenges of remote sensing and computer vision lies in the three-dimensional (3-D) reconstruction of individual trees by using automated methods through very high-resolution (VHR) data sets. However, a successful and complete 3-D reconstruction relies on precise delineation of the trees in two dimensions. In this paper, we present an original approach to detect and delineate citrus trees using unmanned aerial vehicles based on photogrammetric digital surface models (DSMs). The symmetry of the citrus trees in a DSM is handled by an orientation-based radial symmetry transform which is computed in a unique way. Next, we propose an efficient strategy to accurately build influence regions of each tree, and then we delineate individual citrus trees through active contours by taking into account the influence region of each canopy. We also present two efficient strategies to filter out erroneously detected canopy regions without having any height thresholds. Experiments are carried out on eight test DSMs composed of different types of citrus orchards with varying densities and canopy sizes. Extensive comparisons to the state-of-the-art approaches reveal that our proposed approach provides superior detection and delineation performances through supporting a nice balance between precision and recall measures.     

Empirical assessment of rockfall and debris flow risk along the Karakoram Highway,Pakistan

Natural Hazards - The Karakoram Highway links north Pakistan with southwest China. It passes through unique geomorphological, geological and tectonic setting. This study focused 200-km-long section...