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.     

Isolation of <Emphasis Type="Italic">Bacillus cereus</Emphasis> from botanical soil and subsequent biodegradation of waste engine oil

Waste engine oil causes a vital environmental pollution when it spill during change and transportation and products of waste engine oil causes lethal effects to the living systems. Thus, abiotic and biotic approaches are being extensively used for removal of waste engine oil pollution. Therefore in present study, waste engine oil degradation was accomplished by a new bacterial culture, isolated from the soil by an enrichment technique. Morphological, biochemical and gene sequence analysis revealed that isolate was Bacillus cereus. Subsequently, biodegradation potential of B. cereus for waste engine oil was studied. Experimental variables, such as pH, substrate concentration, inoculum size, temperature and time on the biodegradation, were checked in mineral salt medium. The biodegradation efficiency of B. cereus was determined by gravimetry, UV–visible spectrophotometry and gas chromatography. In addition, waste engine oil was also characterized by GC–MS and FTIR for its major constituents, which showed total 38 components in waste engine oil, including hopanes, benzopyrene, long-chain aliphatic hydrocarbons, dibenzothiophenes, biphenyl and their derivatives. Results of successive biodegradation indicated that B. cereus was capable to degrade 1% of waste engine oil with 98.6% degradation potential at pH 7 within 20 days. Hence, B. cereus presents an innovative tool for removing the engine oil from the contaminated area.     

Impact of drilling mud on chemistry and microbiology of an Upper Triassic groundwater after drilling and testing an exploration well for aquifer thermal energy storage in Berlin (Germany)

After completion of an exploration well, sandstones of the Exter Formation were hydraulically tested to determine the hydraulic properties and to evaluate chemical and microbial processes caused by drilling and water production. The aim was to determine the suitability of the formation as a reservoir for aquifer thermal energy storage. The tests revealed a hydraulic conductivity of 1–2 E-5 m/s of the reservoir, resulting in a productivity index of 0.6–1 m³/h/bar. A hydraulic connection of the Exter Formation to the overlaying, artesian “Rupelbasissand” cannot be excluded. Water samples were collected for chemical and microbiological analyses. The water was similarly composed as sea water with a maximum salinity of 24.9 g/L, dominated by NaCl (15.6 g/L Cl and 7.8 g/L Na). Until the end of the tests, the water was affected by drilling mud as indicated by the high pH (8.9) and high bicarbonate concentration (359 mg/L) that both resulted from the impact of sodium carbonate (Na₂CO₃) additives. The high amount of dissolved organic matter (>?58 mg/L) and its molecular-weight distribution pattern indicated that residues of cellulose, an ingredient of the drilling mud, were still present at the end of the tests. Clear evidence of this contamination gave the measured uranine that was added as a tracer into the drilling mud. During fluid production, the microbial community structure and abundance changed and correlated with the content of drilling mud. Eight taxa of sulfate-reducing bacteria, key organisms in processes like bio-corrosion and bio-clogging, were identified. It can be assumed that their activity will be affected during usage of the reservoir.     

Geospatial assessment of soil erosion intensity and sediment yield: a case study of Potohar Region,Pakistan

Estimation of spatial extent of soil erosion, one of the most serious forms of land degradation, is critical because soil erosion has serious implications on soil fertility, water ecosystem, crop productivity and landscape beauty. The primary objective of the current study was to assess and map the soil erosion intensity and sedimentation yield of Potohar region of Pakistan. Potohar is the rainfed region with truncated and complex topography lying at the top of the Indus Basin, the world’s largest irrigation networks of canals and barrages. Spatially explicit Revised Universal Soil Loss Equation (RUSLE) Model integrated with Remote Sensing-GIS techniques was used for detecting/mapping of erosion prone areas and quantification of soil losses. The results show that the Potohar region is highly susceptible to soil erosion with an average annual soil loss of 19 tons ha^?1 year^?1 of which the maximum erosion (70–208 tons ha^?1 year^?1) was near the river channels and hilly areas. The sediment yield due to the erosion is as high as 148 tons ha^?1 year^?1 with an average of 4.3 tons ha^?1 year^?1. It was found that 2.06% of the total area falls under severe soil erosion, 13.34% under high erosion, 15.35% under moderate soil erosion while 69.25% of the area lies in the low (tolerable) soil erosion. Chakwal and Jhelum districts of the region are seriously affected by erosion owing to their topography and soil properties. The information generated in this study is a step forward towards proper planning and implementation of strategies to control the erosion and for protection of natural resources. It is, hence, necessary that suitable water harvesting structures be made to control water to prevent soil erosion and provision of water in the lean season in this region. Tree plantation and other erosion control practices such as strip cropping can also minimize soil erosion in this region.     

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.     

Hydrothermal alteration and evolution of Zr-Th-U-REE mineralization in the microgranite of Wadi Ras Abda,North Eastern Desert,Egypt

Ras Abda plutonic suite, North Eastern Desert of Egypt, consists predominantly of Neoproterozoic calc-alkaline older granites. Minor exposures of pink microgranite are occurring along Wadi Ras Abda within the older granites. Previous studies on this area demonstrated that the microgranite is altered in some parts and contains anomalous concentrations of rare metal elements (Zr, Th, and U). These altered and mineralized zones are re-assessed using field observations, chemical analysis, and by the application of various transmitted light and electron microscopic techniques. The rare metals exist as mineral segregation grew freely into open cavities of the microgranite and concordant with the NNE strike-slip fault movement. The mineralized zones contain an assemblage of secondary magnetite, zircon, uranothorite, columbite-(Mn), fergusonite-(Y), and allanite-(Ce). The extreme abundance of zircon in the mineralized zone, along with other evidence, indicates a hydrothermal origin of this zircon together with associated rare metals. The geochemical investigation and mass balance calculations revealed extreme enrichment of Zr, Th, U, Y, Nb, Ta, and REE. Post-magmatic hydrothermal alterations resulted in such pronounced chemical and mineralogical heterogeneity. The hydrothermal fluids are thought to be oxidizing, alkaline and of medium temperature (>?250 °C). The average contents of the elements Zr (1606 ppm), Th (1639 ppm), U (306 ppm), Nb (955 ppm), and REE (1710 ppm) in the mineralized microgranite reach sub-economic levels and could be a potential source of these elements.     

2D seismic reflection contribution to structural and geometric study of Cenozoic aquifer systems in the Sisseb El-Alem basin,central-eastern Tunisia

The Sisseb-ElAlem basin in central-eastern Tunisia contains several aquifers intensely used in agriculture. Highly exploited shallow aquifers are vulnerable to pollution and drought and reveal high salinity. Consequently, the new hydrogeological surveys were oriented to deep aquifers. The present work corresponds to study of 2D seismic sections, calibrated by water and petroleum wells, and geological outcrops. It highlights geometry and faulting of buried and preserved structures in depth and leads to identify the spatial distribution of Cenozoic potential reservoir layers in the different paleogeographic domains. The Sisseb-ElAlem basin is affected mainly by near E- and N-trending faults. Iso-depth and iso-thickness mapping were used to follow lateral and vertical variations of sedimentary series related to tectonics. The study area consists of Sisseb, Sbikha and El Ktifa uplifts, and El Alem, and Bled Saadia subsiding synclines. The Oligocene-Miocene reservoir series are exploited locally at the high borders. They reveal thickening and deepening southward in the El Alem-Bled Saadia domain, where they can be exploited at about 400 to 2000 m depths.This study should contribute to the orientation of water research to deeper aquifers, and therefore, ensure sufficient drinking and irrigation water supply mainly in areas called for sustainable development.     

Investigation of monotonic and cyclic behavior of sand using a bounding surface plasticity model

Developing the pore water pressures in loose to medium sands below the water table may lead to liquefaction during earthquakes. The simulation of liquefaction (cyclic mobility and flow liquefaction) in sandy soils is one of the major challenges in constitutive modeling of soils. This paper presents the simulation of sand behavior using a critical state bounding surface plasticity model (Dafalias and Manzari’s model, 2004) during monotonic and cyclic loading. The drained, undrained, and cyclic triaxial tests were simulated using Dafalias and Manzari’s model. The simulation results showed that the model predicts behavior of sand, reasonably well. Also, for CSR?<?0.2, number of cycles for liquefaction is significantly increased. The residual strength of Babolsar sand is produced when it is deformed to an axial strain of 20 to 25%.     

Characterization of aerosol type based on aerosol optical properties over Baghdad,Iraq

Aerosol optical depth (AOD), Angstrom exponent (AE), and ozone monitoring instrument aerosols index (OMI-AI) data, derived from MODerate Resolution Imaging Spectroradiometer (MODIS) and OMI sensor on board NASA’s Aqua satellite and NASA-Aura satellite platforms, have been analyzed and classified over Baghdad, Iraq, for an 8-year period (2008–2015). In order to give an obvious understanding of temporal inconsistency in the characteristics and classification of aerosols during each season separately, PREDE POM-02 sky radiometer measurements of AOD, carried out during a 2-year period (2014–2015), were compared with MODIS–Aqua AODs. On seasonal bases, MODIS–Aqua AODs corroborate well with ground-based measurements, with correlation coefficients ranging between 0.74 and 0.8 and RMSE ranging from 0.097 to 0.062 during spring and autumn seasons respectively. The overall satellite- and ground-based measurement comparisons showed a good agreement with correlation coefficients of 0.78 and RMSE of 0.066. These results suggest that MODIS–Aqua gives a good estimate of AOD. Analysis of MODIS–Aqua data for the 8-year period showed that the overall mean AOD, AE, and OMI-AI over Baghdad were 0.44?±?0.16, 0.77?±?0.29, and 1.34?±?0.33 respectively. AOD records presented a unique peak which was extended from mid-spring (April) to mid-summer (July) while the AE annual variability indicated a more complicated behavior with minimum values during the period from late spring (May) to early autumn (September). The maximum AOD and OMI-AI values occurred during summer while their minimum values occurred during winter. The AE showed an opposite behavior to that of AOD such that the highest AE values occurred during autumn and winter and the lowest values happened during spring and summer. This behavior may be attributed to the domination of coarse aerosol particles during autumn and winter seasons and fine aerosol particles during spring and summer seasons. A Hybrid Single-Particle Lagrangian Integrated Trajectory model was utilized to determine the source of air mass transport and to recognize the variability of aerosol origin regions. Finally, AOD, AE, and OMI-AI values have been employed to identify several aerosol types and to present seasonal heterogeneity in their contribution based on their origins.