A fine-granularity scheduling algorithm for parallel XDraw viewshed analysis

Viewshed analysis is widely used in many terrain applications such as siting problem, path planning problem, and etc. But viewshed computation is very time-consuming, in particular for applications with large-scale terrain data. Parallel computing as a mainstream technique with the tremendous potential has been introduced to enhance the computation performance of viewshed analysis. This paper presents a revised parallel viewshed computation approach based on the existing serial XDraw algorithm in a distributed parallel computing environment. A layered data-dependent model for processing data dependency in the XDraw algorithm is built to explore scheduling strategy so that a fine-granularity scheduling strategy on the process-level and thread-level parallel computing model can be accepted to improve the efficiency of the viewshed computation. And a parallel computing algorithm, XDraw-L, is designed and implemented taken into account this scheduling strategy. The experimental results demonstrate a distinct improvement of computation performance of the XDraw-L algorithm in this paper compared with the coarse-partition algorithm, like XDraw-E which is presented by Song et al. (Earth Sci Inf 10(5):511–523, 2016), and XDraw-B that is the basic algorithm of serial XDraw. Our fine-granularity scheduling algorithm can greatly improve the scheduling performance of the grid cells between the layers within a triangle 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.     

Hydrology models approach to estimation of the groundwater recharge: case study in the Bulgarian Danube watershed

The groundwater (GW) makes an important part of a region runoff. GW bodies playing the role of accumulating reservoirs regulate the GW discharge enabling the river flow to have more uniform long-term distribution. Along with other important advantages, the GW offers the users stable water abstraction rate independent from the recharge rate. The GW recharge quantification belongs to the uneasy tasks in the water resource management. Applying the conventional methods needs multiyear observation records of the variation of the groundwater body (GWB) characteristics. The employment of hydrology models avoids that necessity but requires great amount of data related to the soil hydraulic properties, the land topography and cover of the GWB watershed and long-term records of the climatic effects. The paper presents an introduction of the mathematical model CLM3 into the GW recharge estimation problem. It is a complex and advanced model with adequate interpretation of the water-related processes in the soil and on the land surface under atmospheric effects. The input is available from NCEP/NCAR reanalysis atmosphere data and the International Geosphere-Biosphere Program (IGBP) data base. The model is applied to GW recharge assessment of the Bulgarian Danube district for the year 2013. The obtained monthly and yearly total district values and the areal distribution of the infiltration intensity are matched to the existing field observation-based estimates. The study shows that the CLM3 model approach leads to encouraging results. The method comes very useful with GWB lacking regime observation data as well as for GW recharge prognostic assessments under climatic scenarios.     

Accessing effects and signals of leakage from a CO<Subscript>2</Subscript> reservoir to a shallow freshwater aquifer by reactive transport modelling

CO₂ geological storage is a transitional technology for the mitigation of climate change. In the vicinity of potential CO₂ reservoirs in Hungary, protected freshwater aquifers used for drinking water supplies exist. Effects of disaster events of CO₂ escape and brine displacement to one of these aquifers have been studied by kinetic 1D reactive transport modelling in PHREEQC. Besides verifying that ion concentrations in the freshwater may increase up to drinking water limit values in both scenarios (CO₂ or brine leakage), total porosity of the rock is estimated. Pore volume is expected to increase at the entry point of CO₂ and to decrease at further distances, whereas it shows minor increase along the flow path for the effect of brine inflow. Additionally, electrical conductivity of water is estimated and suggested to be the best parameter to measure for cost-effective monitoring of both worst-case leakage scenarios.     

Study on the natural radioactivity level of stone coal-bearing strata in East China

Based on γ-radiation dose rate and radon concentration measurements and ²³⁸U, ²³²Th, ²²⁶Ra, and ⁴⁰K radionuclide testing, this study identifies the radioactive anomalies of stone coal-bearing strata in East China and evaluates the natural radioactivity levels in the air, solid, water and plant media in the typical area of the regional stone coal-bearing layers. The stone coal-bearing strata in East China occur in the lower Cambrian system along the margin of the Yangtze block; additionally, the radioactive anomaly area is sporadically distributed in the stone coal-bearing layers. The background values of ²³⁸U, ²³²Th, ²²⁶Ra, and ⁴⁰K are higher in the stone coal-bearing areas, and the spatial distribution of these natural radionuclides shows significant variability. ²³⁸U and ²²⁶Ra clearly accumulate in the coal, coal gangue and soil and are the main sources of the environmental radiation in coal mines. The γ-radiation shows a higher background value in the stone coal-bearing area, and this radioactive pollution cannot be ignored. Typically, the effective dose of γ-radiation exceeds the limit value of 5 mSv/a, and the total α and total β concentrations of the groundwater are 10–30 times the limit value at some points. The residents near the mining area are subjected to a higher radiation dose, and the groundwater, building materials, and plants have been contaminated by the radioactive pollution sporadically through time. It is necessary to strengthen the monitoring work of radioactive environments and to take appropriate control measures.     

Sealing efficiency analysis for shallow-layer caprocks in CO<Subscript>2</Subscript> geological storage

The CO₂ migrated from deeper to shallower layers may change its phase state from supercritical state to gaseous state (called phase transition). This phase transition makes both viscosity and density of CO₂ experience a sharp variation, which may induce the CO₂ further penetration into shallow layers. This is a critical and dangerous situation for the security of CO₂ geological storage. However, the assessment of caprock sealing efficiency with a fully coupled multi-physical model is still missing on this phase transition effect. This study extends our previous fully coupled multi-physical model to include this phase transition effect. The dramatic changes of CO₂ viscosity and density are incorporated into the model. The impacts of temperature and pressure on caprock sealing efficiency (expressed by CO₂ penetration depth) are then numerically investigated for a caprock layer at the depth of 800 m. The changes of CO₂ physical properties with gas partial pressure and formation temperature in the phase transition zone are explored. It is observed that phase transition revises the linear relationship of CO₂ penetration depth and time square root as well as penetration depth. The real physical properties of CO₂ in the phase transition zone are critical to the safety of CO₂ sequestration. Pressure and temperature have different impact mechanisms on the security of CO₂ geological storage.     

Socioenvironmental issues of river bed material extraction in the Himalayan piedmont (India)

The present study focuses on the Balason river running through the Himalayan piedmont zone (near Siliguri, India). The objective of the study is an assessment of the environmental effects of river bed material extraction by humans and the dependence of indigenous people on the river and its ecosystem services. The analysis is based on results of field work consisting of geodetic measurements of the river channel and interviews among the local community from the Nimtijot village. Historical hydrological data were also used for the study. The results of the investigation show that the Balason river is heavily affected by excessive exploitation of river bed material during dry season and the replenishment of extracted material in a monsoon season is not always sufficient. It leads to channel deepening. The local community working in the river heavily depends for its livelihood on continuing this activity. A decreasing amount of bed material to be extracted may lead to degradation of the strong relationship between the local community and their natural environment (river).     

Development of efficiently coupled thermo-hydro-mechanical model to predict hydraulic fracture morphology in heavy oil reservoirs

The aim of the study involves examining the effect of heavy oil viscosity on fracture geometry in detail by establishing a heavy oil fracturing model and conventional fracturing model based on thermal–hydraulic–mechanical (THM) coupled theory, Walther viscosity model, and K–D–R temperature model. We consider viscosity and density within the heavy oil fracturing model as functions of pressure and temperature while that as constants within the conventional fracturing model. A heavy oil production well is set as an example to analyze the differences between the two models to account for the thermo-poro-elastic effect. The results show that temperature exhibits the most significant influence on the heavy oil viscosity while the influence of pressure is the least. In addition, a cooling area with a width of 0–1 m and varied length is generated near the fracture. The heavy oil viscosity increases sharply in this area, thereby indicating an area of viscosity increment. The heavy oil viscosity increases faster and is closer to wellbore, and a high viscosity increment reduces the mobility of the heavy oil and prevents the fracturing fluid from entering into the reservoir. The special viscosity distribution results in significant differences in pore pressure, oil saturation, and changing trends between these two models. In the heavy oil reservoir fracturing model, the thermal effect completely exceeds the influence of pore elasticity, and the values of the fracture length, width, and static pressure exceed those calculated in the conventional fracturing model. Thus, a comparison of the measured values indicates that the results obtained by considering viscosity as a function of temperature and pressure are more accurate. Therefore, the results of this study are expected to provide good guidelines for the design of heavy oil fracturing.     

Synthesis and characterisation of a multifunctional oil-based drilling fluid additive

An oil-based drilling fluid additive H-DEA (or humic acid-cocamide diethanolamine) was synthesised using humic acid and cocamide diethanolamine as raw materials. The rheological behaviors of H-DEA showed that the synthesised product has the good properties in both decreasing the filtrate loss and improving rheology property of oil-based drilling fluids compared with other commercially available additives. Under the optimal additive amount of 3%, both API filtrate loss and yield point changed remarkably from 5.40 to 0.41 mL and 9.0 to 25.6 Pa, respectively. Furthermore, differential scanning calorimetry (DSC) showed that H-DEA has good thermal stability in a wide temperature range up to 170 °C. Infrared spectroscopy (IR) and rheological analysis revealed that the possible mechanism of the multifunctional effects may be attributed to the existing of high density of strong polar groups, hydrogen bonds, electrostatic forces, and intermolecular association on H-DEA molecular structure. The results of the study showed that the synthesised H-DEA can be potentially used as a multifunctional oil-based drilling fluid additive in oil-drilling excavation.