Geostatistics as a groundwater exploration planning tool: case of a brackish-saline aquifer

This paper discusses a geostatistical approach to model a groundwater aquifer in 3-D. The study aims at utilizing geostatistics as a tool for characterizing zones of better-water quality in a brackish-saline aquifer. In particular, the geostatistical model was constructed to characterize the aquifer’s salinity, represented by total dissolved solids (TDS), using logs of porosity and resistivity. Quality-checked estimated TDS vertical profiles were employed to construct and model horizontal and vertical semivariograms. Parameters of semivariogram models were used to develop both the kriging plan and the generated model. Results of this modeling process are shown in the form of horizontal salinity distribution maps. The aquifer was sliced into 20 layers, each 20 m thick, to represent its overall thickness. Salinity layers maps reflect vertical stratification of TDS concentrations in the aquifer and show that water quality deteriorates with depth and toward the northern part of the aquifer. Relatively better-quality water (TDS ≤10,000 mg/l) can be found at depths between 100 and 250 m below the aquifer’s top in both eastern and southeastern parts. Water in the same interval to the western and southwestern parts reflects the presence of higher TDS concentration. From a planning point of view, it is more feasible to target the eastern part of the aquifer for pumping and desalination purposes. In addition, the generated model could be utilized as an initial condition for flow simulation.