A Dual Delineation Approach to Untangle Interferences in Well Doublets Systems

Traditional numerical methods for the delineation of wellhead protection areas span deterministic and probabilistic approaches. They provide time-related capture zones. However, none of the existing approaches identifies the groundwater contribution areas related to each source or sink. In this work, the worthiness of the so-called double delineation approach was extended. This task was achieved by simple postprocessing of its dual outputs leading to a highly efficient screening tool. In the particular context of geothermal resources management through the well doublets of the Dogger aquifer in the Paris Basin (France), the approach was extended to forecast the compositional heat breakthrough at production wells. Hence, cold-water breakthrough and temperature decline in production wells are timely assessed in low-enthalpy geothermal reservoirs. The method quantifies how groundwater volumes are moving through space and time between any couple of source and sink. It provides unprecedented tools advancing the enhanced understanding of water resources systems functioning. It is highly recommended to implement the presented concepts in the current and future generations of community groundwater models.     

Lithology-dependent seismic anisotropic amplitude variation with offset correction in transversely isotropic media

Analysis of amplitude variation with offset is an essential step for reservoir characterization. For an accurate reservoir characterization, the amplitude obtained with an isotropic assumption of the reservoir must be corrected for the anisotropic effects. The objective is seismic anisotropic amplitude correction in an effective medium, and, to this end, values and signs of anisotropic parameter differences (Δδ and Δε) across the reflection interfaces are needed. These parameters can be identified by seismic and well log data. A new technique for anisotropic amplitude correction was developed to modify amplitude changes in seismic data in transversely isotropic media with a vertical axis of symmetry. The results show that characteristics of pre-stack seismic data, that is, amplitude variation with offset gradient, can be potentially related to the sign of anisotropic parameter differences (Δδ and Δε) between two layers of the reflection boundary. The proposed methodology is designed to attain a proper fit between modelled and observed amplitude variation with offset responses, after anisotropic correction, for all possible lithofacies at the reservoir boundary. We first estimate anisotropic parameters, that is, δ and ε, away from the wells through Backus averaging of elastic properties resulted from the first pass of isotropic pre-stack seismic inversion, on input data with no amplitude correction. Next, we estimate the anisotropic parameter differences at reflection interfaces (values and signs of Δδ and Δε). We then generate seismic angle gather data after anisotropic amplitude correction using Rüger's equation for the P-P reflection coefficient. The second pass of isotropic pre-stack seismic inversion is then performed on the amplitude-corrected data, and elastic properties are estimated. Final outcome demonstrates how introduced methodology helps to reduce the uncertainty of elastic property prediction. Pre-stack seismic inversion on amplitude-corrected seismic data results in more accurate elastic property prediction than what can be obtained from non-corrected data. Moreover, a new anisotropy attribute (ν) is presented for improvement of lithology identification.     

Transient and Transition Factors in Modeling Permafrost Thaw and Groundwater Flow

Permafrost covers approximately 24% of the Northern Hemisphere, and much of it is degrading, which causes infrastructure failures and ecosystem transitions. Understanding groundwater and heat flow processes in permafrost environments is challenging due to spatially and temporarily varying hydraulic connections between water above and below the near-surface discontinuous frozen zone. To characterize the transitional period of permafrost degradation, a three-dimensional model of a permafrost plateau that includes the supra-permafrost zone and surrounding wetlands was developed. The model is based on the Scotty Creek basin in the Northwest Territories, Canada. FEFLOW groundwater flow and heat transport modeling software is used in conjunction with the piFreeze plug-in, to account for phase changes between ice and water. The Simultaneous Heat and Water (SHAW) flow model is used to calculate ground temperatures and surface water balance, which are then used as FEFLOW boundary conditions. As simulating actual permafrost evolution would require hundreds of years of climate variations over an evolving landscape, whose geomorphic features are unknown, methodologies for developing permafrost initial conditions for transient simulations were investigated. It was found that a model initialized with a transient spin-up methodology, that includes an unfrozen layer between the permafrost table and ground surface, yields better results than with steady-state permafrost initial conditions. This study also demonstrates the critical role that variations in land surface and permafrost table microtopography, along with talik development, play in permafrost degradation. Modeling permafrost dynamics will allow for the testing of remedial measures to stabilize permafrost in high value infrastructure environments.     

Research note: deblended-data reconstruction using generalized blending and deblending models

We introduce a concept of generalized blending and deblending, develop its models and accordingly establish a method of deblended-data reconstruction using these models. The generalized models can handle real situations by including random encoding into the generalized operators both in the space and time domain, and both at the source and receiver side. We consider an iterative optimization scheme using a closed-loop approach with the generalized blending and deblending models, in which the former works for the forward modelling and the latter for the inverse modelling in the closed loop. We applied our method to existing real data acquired in Abu Dhabi. The results show that our method succeeded to fully reconstruct deblended data even from the fully generalized, thus quite complicated blended data. We discuss the complexity of blending properties on the deblending performance. In addition, we discuss the applicability to time-lapse seismic monitoring as it ensures high repeatability of the surveys. Conclusively, we should acquire blended data and reconstruct deblended data without serious problems but with the benefit of blended acquisition.     

Stress distribution from railway track over geogrid reinforced ballast underlain by clay

Railway ballast forms a major component of a conventional rail track and is used to distribute the load to the subgrade, providing a smooth running surface for trains. It plays a significant role in providing support for the rail track base and distributing the load to the weaker layer underneath. Ballast also helps with drainage, which is an important factor for any type of transportation structure, including railroads. Over time, ballast progressively deforms and degrades under dynamic loading and loses its strength. In this study, extensive laboratory tests were conducted to investigate the effect of load amplitude, geogrid position, and number of geogrid layers, thickness of ballast layer and clay stiffness on the behavior of the reinforced ballast layer and induced strains in a geogrid. A half full-scale railway was constructed for carrying out the tests, which consisted of two rails 800 mm in length with three wooden sleepers(900 mm × 10 mm × 10 mm). Three ballast thicknesses of 200, 300 and 400 mm were used in the tests. The ballast was overlying 500 mm thickness clay in two states, soft and stiff. The tests were carried out with and without geogrid reinforcement; the tests were performed in a well-tied steel box of 1.5 m length ×1 m width ×1 m height. Laboratory tests were conducted to investigate the response of the ballast and the clay layers where the ballast was reinforced by a geogrid. Settlement in ballast and clay, soil pressure and pore water pressure induced in the clay were measured in reinforced and unreinforced ballast cases. It was concluded that the amount of settlement increased as the simulated train load amplitude increased, and there was a sharp increase in settlement up to cycle 500. After that, there was a gradual increase that leveled out between, 2500 to 4500 cycles depending on the frequency used. There was a slight increase in the induced settlement when the load amplitude increased from 0.5 to 1 ton but it was higher when the load amplitude increased to 2 tons. The increased amount in settlement depended on the existence of the geogrid and other parameters studied. The transmitted average vertical stress for ballast thicknesses of 30 cm and 40 cm increased as the load amplitude increased, regardless of the ballast reinforcement for both soft and stiff clay. The position of the geogrid had no significant effect on the transmitted stresses. The value of the soil pressure and pore water pressure on ballast thicknesses of 20 cm was higher than for 30 cm and 40 cm thicknesses. This meant that the ballast attenuated the induced waves. The soil pressure and pore water pressure for reinforced and unreinforced ballast was higher in stiff clay than in soft clay.     

Possibility of the Application of Danube,Sava and Morava River Area Shallow Aquifer Groundwater (Serbia) in the Irrigation Purposes

Water Resources - This paper presents the results of the groundwater quality assessment, on the Serbian Danube, Sava and Morava River water area, for irrigation purposes. The analysis was done at...     

Vorticity analysis in the Zagros orogen,Shiraz area,Iran

International Journal of Earth Sciences - Quantitative vorticity analyses in orogenic belts are essential for studying the kinematics of deformation and can be performed using a range of methods....     

Types and nature of fracture associated with Late Ordovician paleochannels of glaciofluvial Sarah Formation,Qasim region,Central Saudi Arabia

This study evaluates the Late Ordovician glaciofluvial deposits of the Sarah Formation and equivalent outcroppings in north, central, and southwestern Saudi Arabia. The Sarah Formation also covers a wide area in the subsurface and is considered as an important target for unconventional tight gas reservoir. Defining the fracture types, nature, and distribution in outcrop scale might help to establish a successful fracture simulation model and behavior for the Sarah tight gas reservoir in the subsurface. This study investigates fracture characteristics for the Sarah Formation at Sarah paleochannel outcrops. The study revealed three sets of fractures, which have EW, NS, and SE-NW directions, and these fractures vary from open, resistive, and filled to resistive fractures. The closed fractures are filled with ferruginated iron oxides and gypsum. The filled fractures (the thrust boundary) are found in the study area at the SE-NW strike fracture set, while open and resistive fractures are found mainly at S-N and E-W fracture sets, respectively. The syn-depositional filled fractures (iron oxides) are considered as the younger fracture sets while the open and resistive fractures are post-depositional fractures which may have resulted from uplift or tectonic movement. A general model representing the fracture pattern and the thrusting boundaries due to glacial movement was constructed. It has been noticed that the systematic occurrence of filled fractures (thrust boundaries) described the boundaries between different glacial events, which act as a fluid barrier (filled fractures) and decrease the reservoir quality. The finding of this study might be utilized as a guide and lead for exploration in the subsurface Sarah glacial deposits. It will also help to understand and speculate the nature pattern and distribution of fractures with the Sarah Formation.     

Application of a 3D photorealistic model for the geological analysis of the Permian carbonates (Khuff Formation) in Saudi Arabia

Three dimensional (3D) photorealistic models of geological outcrops have the potential to enhance the teaching of earth sciences by providing scale models in a virtual reality environment. These models can be run on low-cost desktop computers. Photorealistic models for geological outcrops are a digital illustration of outcrop photographs with either a point cloud representation or Triangular Irregular Network (TIN) mesh of the outcrop surface. The level of detail for these models is dependent on the target resolutions (physical and optical) that were used during data acquisition. In addition, the technique in which the data is rendered as a digital model affects the level of detail that can be observed by the geologists. A colored point cloud representation is suitable for large-scale features, but fine details are lost when the geologist zooms in to view the model close up. In contrast, a photorealistic model that is constructed from photographs draped onto a triangle mesh surface derived from Light Detection and Ranging (LiDAR) point clouds provides a level of detail that is restricted only by the resolution of the photographs.     

Discussion of “Interfacing the geographic information system,remote sensing,and the soil conservation service-curve number method to estimate curve number and runoff volume in the ASIR region of Saudi Arabia” by Fawzi S. Mohammad,Jan Adamowski

The soil conservation service curve number (SCS-CN) method is one of the most commonly used methods to compute the direct runoff from a rainfall event. Since the method was established, numerous researches were undertaken to improve the method through accurate estimation of its parameter and especially the curve number (CN). However, the essence of the SCS method, as an event-based Hortonian mechanism method, remained unchanged. The main assumption of the method related to the rainfall input is that the rainfall is continuous in time and uniform over the watershed. Mohammad and Adamowski (2015) paper apparently used the SCS method to estimate the annual runoff using the annual rainfall as one cumulative rainfall input value, which is a violation of the event-based principle of the method and of the assumption of the continuity of the rainfall event.To re-estimate the average annual runoff more realistically for the Asir region, Saudi Arabia, daily rainfall data from 14 rainfall stations are used for calculating the resulting runoff depths, on a daily event-by-event rainfall basis, throughout the whole simulation period. The resulting runoff depths are added for each year, and the total cumulative annual runoff values for each year are averaged to get the average annual runoff. The runoff values based on the previously mentioned procedure are an upper limit of the actual average annual runoff as the underlying SCS equations discard evaporation and similar long-term losses. Nevertheless, the average runoff values obtained in the discussion paper are an order of magnitude (at least five to tenfold) lower than the ones of the original paper. An equation is proposed to obtain a more realistic estimate of the average annual runoff, to be used with the average annual rainfall as an input, if the annual value is the only available rainfall information.