FASTREACT – An efficient numerical framework for the solution of reactive transport problems

In the framework of safety assessment studies for geological disposal, large scale reactive transport models are powerful inter-disciplinary tools aiming at supporting regulatory decision making as well as providing input to repository engineering activities. Important aspects of these kinds of models are their often very large temporal and spatial modelling scales and the need to integrate different non-linear processes (e.g., mineral dissolution and precipitation, adsorption and desorption, microbial reactions and redox transformations). It turns out that these types of models may be computationally highly demanding. In this work, we present a Lagrangian-based framework, denoted as FASTREACT, that aims at solving multi-component-reactive transport problems with a computationally efficient approach allowing complex modelling problems to be solved in large spatial and temporal scales. The tool has been applied to simulate radionuclide migration in a synthetic heterogeneous transmissivity field and the results have been successfully compared with those obtained using a standard Eulerian approach. Finally, the same geochemical model has been coupled to an ensemble of realistic three-dimensional transport pathways to simulate the migration of a set of radionuclides from a hypothetical repository for spent nuclear fuel to the surface. The results of this modelling exercise, which includes key processes such as the exchange of mass between the conductive fractures and the matrix, show that FASTREACT can efficiently solve large-scale reactive transport models.     

云南省地质灾害的发育特征及防治对策与建议

由于地质环境恶化，近年来崩塌、滑坡、泥石流、地面塌陷等灾害在云南频繁发生，严重危及灾区人民的生命财产安全并严重影响经济建设的顺利进行。从云南的环境现状，地震及新构造运动，地层岩性及其岩土体工程地质特性，气象、水文等自然环境格局的背景因素(内因及诱因)入手，举例分析影响我省环境恶化的外因，在此基础上，提出防灾减灾的对策建议。     

An intermediate-scale model for thermal hydrology in low-relief permafrost-affected landscapes

Integrated surface/subsurface models for simulating the thermal hydrology of permafrost-affected regions in a warming climate have recently become available, but computational demands of those new process-rich simu- lation tools have thus far limited their applications to one-dimensional or small two-dimensional simulations. We present a mixed-dimensional model structure for efficiently simulating surface/subsurface thermal hydrology in low-relief permafrost regions at watershed scales. The approach replaces a full three-dimensional system with a two-dimensional overland thermal hydrology system and a family of one-dimensional vertical columns, where each column represents a fully coupled surface/subsurface thermal hydrology system without lateral flow. The system is then operator split, sequentially updating the overland flow system without sources and the one-dimensional columns without lateral flows. We show that the app- roach is highly scalable, supports subcycling of different processes, and compares well with the corresponding fully three-dimensional representation at significantly less computational cost. Those advances enable recently developed representations of freezing soil physics to be coupled with thermal overland flow and surface energy balance at scales of 100s of meters. Although developed and demonstrated for permafrost thermal hydrology, the mixed-dimensional model structure is applicable to integrated surface/subsurface thermal hydrology in general.     

Evaluation of retention properties of a semi-synthetic fractured block from modelling at performance assessment time scales (Äspö Hard Rock Laboratory, Sweden)

Modelling of radionuclide transport in fractured media is a primary task for safety evaluation of a deep nuclear waste repository. A performance assessment (PA) model has been derived from site characterization data with the aim of improving confidence for quantifying transport of sorbing radionuclides at a safety time scale of 10⁶ y. The study was conducted on a 200?×?200?×?200 m semi-synthetic fractured block, providing a realistic system derived from the Äspö Hard Rock Laboratory (Sweden) dataset. The block includes 5,632 fractures ranging from 0.5 to 100 m in length and a heterogeneous matrix structure (fracture coating, gouge, mylonite, altered and non-altered diorite). The PA model integrates steady-state flow conditions and transport of released radionuclides during the safety time scale. An original simulation method was developed involving Eulerian flow and transport within fracture planes with a mixed hybrid finite element scheme and a semi-analytical source term to account for heterogeneous matrix diffusion. Total mass flux of radionuclides (conservative to strongly sorbing) was computed. A method to simplify the system was demonstrated, leading to a major path of 12 fractures. This is mainly due to the low connectivity of the fracture network. Matrix diffusion and sorption proved to have major impact on block retention properties for PA conditions.     

Mathematical modelling of groundwater flow at Sellafield, UK

Sellafield in West Cumbria was a potential site for the location of the UK's first underground repository for radioactive, intermediate level waste (ILW). The repository was to lie around 650 m beneath the ground surface within rocks of the Borrowdale volcanic group (BVG), a thick suite of SW dipping, fractured, folded and metamorphosed Ordovician meta-andesites and ignimbrites. These are overlain by an onlapping sequence of Carboniferous and Permo-Triassic sediments. In situ borehole measurements showed that upward trending fluid pressure gradients exist in the area of the potential repository site, and that there are three distinct fluid types in the subsurface; fresh, saline and brine (at depth, to the west of the site). Simulations of fluid flow in the Sellafield region were undertaken with a 2D, steady-state, coupled fluid and heat flow simulation code (OILGEN). In both simplified and geologically complex models, topographically driven flow dominated the regional hydrogeology. Fluids trended persistently upwards through the potential repository site. The dense brine to the west of the site promoted upward deflection of topographically driven groundwaters. The inclusion in hydrogeological models of faults and variably saline sub-surface fluids was essential to the accurate reproduction of regional hydraulic head variations. Sensitivity analyses of geological variables showed that the rate of groundwater flow through the potential repository site was dependent upon the hydraulic conductivity of the BVG, and was unaffected by the hydraulic conductivity of other hydrostratigraphic units. Calibration of the model was achieved by matching simulated subsurface pressures to those measured in situ. Simulations performed with BVG hydraulic conductivity 100 times the base case median value provided the “best-fit” comparison between the calculated equivalent freshwater head and that measured in situ, regardless of the hydraulic conductivity of other hydrostratigraphic units. Transient mass transport simulations utilising the hydraulic conductivities of this “best fit” simulation showed that fluids passing through the potential repository site could reach the surface in 15 000 years. Simple safety case implications drawn from the results of the study showed that the measured BVG hydraulic conductivity must be less than 0.03 m year⁻¹ to be simply declared safe. Recent BVG hydraulic conductivity measurements showed that the maximum BVG hydraulic conductivity is around 1000 times this safety limit.     

Groundwater flow modeling for effective implementation of landslide stability enhancement measures

This paper deals with groundwater hydrology at a prominent fracture zone landslide slope (Nuta–Yone landslides) in Japan with an objective to explore an efficient method for the application of landslide stability enhancement measures. The correlation analyses between the hydrological parameters and ground surface movement data at this landslide resulted in low correlation values indicating that the geological formation of the area is extremely complex. For the purpose of understanding the groundwater flow behavior in the landslide area, a three-dimensional transient groundwater flow model was prepared for a part of the landslide slope, where the levels of effectiveness of applied landslide stability enhancement measures (in the form of multilayered deep horizontal drains) are different, and was calibrated against the measured water surface elevations at different piezometer locations. The parameter distributions in the calibrated model and the general directions of the groundwater flow in terms of flow vectors and the results of particle tracking at the model site were interpreted to understand the reasons for variations in effectiveness of existing landslide stability enhancement measures and to find potentially better locations for the implementation of future landslide stability enhancement measures. From the modeling results, it was also understood that groundwater flow model can be effectively used in better planning and locating the landslide stability enhancement measures.     

Assessing impact of urban impervious surface on watershed hydrology using distributed object-oriented simulation and spatial regression

In this study, we investigated the relationship between watershed characteristics and hydrology using high spatial resolution impervious surface area (ISA), hydrologic simulations and spatial regression. We selected 20 watersheds at HUC 12 level with different degrees of urbanization and performed hydrologic simulation using a distributed object-oriented rainfall and runoff simulation model. We extracted the discharge per area and ratio of runoff to base flow from simulation results and used them as indicators of hydrology pattern. We derived percentage of ISA, distance from ISA to streams, and stream density as the watershed characteristics to evaluate the relationship with hydrology pattern in watersheds using ordinary least square, spatial error and spatial lag regression models. The comparison indicates that spatial lag regression model can achieve better performance for the evaluation of relationship between ratio of runoff to base flow and watershed characteristics, and that three models provide similar performance for the evaluation of relationship between discharge per area and watershed characteristics. The results from regression analyses demonstrate that ISA plays an important role in watershed hydrology. Ignorance of spatial dependence in analyses will likely cause inaccurate evaluation for relationship between ISA and watershed hydrology. The hydrologic model, regression methods and relationships between watershed characteristics and hydrology pattern provide important tools and information for decision makers to evaluate the effect of different scenarios in land management.     

An initial screening tool for water resource contamination due to development in the Olympic Park 2012 site,London

Groundwater is the primary source of potable water in southeast England. Its protection in urban environments is of paramount importance. Following a scoping study the British Geological Survey (BGS) established a project to develop an initial screening tool (IST) to assist the planning community in the assessment of the potential risk to ground and surface waters from contaminants mobilised as a consequence of redevelopment. The tool has been designed in the context of the source-pathway-receptor paradigm that informs Part IIa of the UK Environmental Protection Act (1990). Building on the work of previous screening tools and in particular ConSEPT, a BGS contaminated site evaluation and prioritisation tool, the IST incorporates significant refinements to scoring methodologies and takes the prioritisation approach into the 3-D environment. Implemented as a customised GIS application and utilising surfaces extracted from 3-D geological modelling, the tool collates and interrogates a range of geoscientific information, including contaminant scale, geological, historic land use, groundwater level and hydrogeological domain data. The IST facilitates the ranking of various proposed development scenarios through a semi-quantitative assessment of contamination potential, via a number of pollutant linkages, providing planners with reports on the type, spatial distribution and hazards associated with potential contaminant sources within their area. To achieve this, a range of evaluation factors applied to the sources, pathways and receptors are scored through a combination of spatial and attribute queries, then assessed on the basis of potential linkages. The initial research area selected for the application of the IST was the Olympic Park site, London.     

Quantitative assessment of radionuclide retention in the Quaternary sediments/granite interface of the Fennoscandian shield (Sweden)

The Quaternary sediments representing the interface between the granite host rock and the Earth surface are of paramount importance when determining the potential cycling of anthropogenic and natural radionuclides in near-surface systems. This is particularly true in the case of high-level nuclear waste (HLNW) repositories placed in granite. In this work a modelling procedure is presented to quantitatively determine the retention capacity of a Quaternary till in the Forsmark area, which has been recently selected to host the deep geologic storage of HLNW in Sweden. Reactive transport numerical models have been used to simulate the intrusion of a deep groundwater carrying radionuclides potentially released from a repository into a Quaternary till. Four radionuclides (²³⁵U, ¹³⁵Cs, ²²⁶Ra and ⁹⁰Sr) have been selected according to their different geochemical behaviour and potential dose relevance to surface ecosystems. Numerical results indicate that repository-derived: (i) U will have a minor impact in the till, mainly due to the high natural concentration of U and its adsorption on ferrihydrite; (ii) Cs will be efficiently retained by cation exchange on illite; (iii) Ra will be retained via co-precipitation with barite; and although (iv) Sr will be retained via co-precipitation with calcite and cation exchange on illite, the retention capacity of the Quaternary till for Sr is limited.     

耦合湿地模块的流域水文模型模拟效率评价

耦合湿地模块的流域水文模型是开展流域湿地变化的水文效应及其水文功能评估的有效工具,其模拟效率直接关系到模拟精度和应用价值。选取中国湿地主要分布区之一的嫩江流域,利用PHYSITEL/HYDROTEL模型平台,构建孤立湿地和河滨湿地模块并与流域水文模型耦合,并从拟合优度指数和模拟效率角度评价了耦合湿地模块的流域水文模型模拟效率。结果表明,耦合湿地模块的流域水文模型,拟合和验证期间模型的拟合度指数均有所提高,嫩江富拉尔基和大赉水文站控制流域的Nash-Sutcliffe系数和克林效率系数平均提高了3.08%和4.64%,均方根误差和相对偏差平均降低了12.72%和55.93%,且整个研究时段模拟结果的流量指数总体更接近于观测的径流数据。可见,耦合湿地模块的流域水文模型可提高水文过程模拟精度,更好地定量评估湿地水文功能,为流域水资源精准调控与湿地修复保护提供重要支撑。     

A model approach to radioactive waste disposal at Sellafield

One of the great environmental problems of our age is the safe disposal of radioactive waste for geological time periods. Britain is currently investigating a potential site for underground burial of waste, near the Sellafield nuclear plant. Future leakage of radionuclides depends greatly on subsurface water flows; these must be understood from the past, to predict hydrogeology 10⁴–10⁵ years into the future. We have taken information from the present-day, published by the government company Nirex, and used a finite-element steady-state fluid flow computer code to examine water flows in the subsurface. We find that flow directions at the planned Repository are persistently upwards, and that geologically significant flow rates could occur. How rates are particularly sensitive to uncertainties of rock permeability (conductivity) measurements made from site investigation boreholes. The hydrogeology at this site needs longer term investigation before a confident and credible prediction can be made.     

地下水饮用水水源地保护区划分关键技术研究

依据《中华人民共和国水污染防治法》和《饮用水水源保护区划分技术规范》(HJ/T338-2007),对饮用水源保护区进行科学合理划分是实现水资源综合开发、合理利用、积极保护、科学管理的基础工作,也是运用法律、行政、经济手段强化水资源管理的根本保证,是防治水污染、保护水资源的重要措施。饮用水源保护区划分涉及经济社会发展、人类活动影响,地表、地下水资源数量,地表、地下水资源水质,陆地水文、工程地质水文,以及水利、环保等学科内容,是一项技术性较强,涉及学科较多的技术工作。通过饮用水源保护区划分涉及关键技术研究和保护区划分工作的开展,合理确定水源的保护范围和保护目标,为地表、地下水资源的合理开发、保护、利用提供科学依据,以保障供水安全、生态与环境安全。     

Development and testing of a contamination potential mapping system for a portion of the General Separations Area, Savannah River Site, South Carolina

 A methodology was developed to evaluate and map the contamination potential or aquifer sensitivity of the upper groundwater flow system of a portion of the General Separations Area (GSA) at the Department of Energy's Savannah River Site (SRS) in South Carolina. A Geographic Information System (GIS) was used to integrate diverse subsurface geologic data, soils data, and hydrology utilizing a stack-unit mapping approach to construct mapping layers. This is the first time that such an approach has been used to delineate the hydrogeology of a coastal plain environment. Unit surface elevation maps were constructed for the tops of six Tertiary units derived from over 200 boring logs. Thickness or isopach maps were created for five hydrogeologic units by differencing top and basal surface elevations. The geologic stack-unit map was created by stacking the five isopach maps and adding codes for each stack-unit polygon. Stacked-units were rated according to their hydrogeologic properties and ranked using a logarithmic approach (utility theory) to establish a contamination potential index. Colors were assigned to help display relative importance of stacked-units in preventing or promoting transport of contaminants. The sensitivity assessment included the effects of surface soils on contaminants which are particularly important for evaluating potential effects from surface spills. Hydrogeologic/hydrologic factors did not exhibit sufficient spatial variation to warrant incorporation into contamination potential assessment. Development of this contamination potential mapping system provides a useful tool for site planners, environmental scientists, and regulatory agencies. Received: 1 April 1997 · Accepted: 4 November 1997     

GPS measurement of land subsidence in Gandhinagar,Gujarat (Western India), due to groundwater depletion

A continuous GPS measurement site, ISRR, at Gandhinagar (Western India), documented ~?5 mm/year of surface subsidence rate during 2009–2016. Preliminary modelling using an analytical solution indicates that the observed surface subsidence rate at the ISRR site is consistent with the ground water depletion reported from Gandhinagar. An assessment of data from GPS sites at Lucknow and Varanasi in the Indo-Gangetic plains in Northern India does not indicate any significant subsidence at these sites which is also consistent with the in situ observations of insignificant depletion of ground water in the region.     

Investigation of groundwater flow in karst areas using component separation of natural potential measurements

 The natural (electrical) potential (NP) method – also known as self-potential, spontaneous potential and streaming potential (SP) – has been used to locate areas of groundwater flow in karst terrane. NP is the naturally occurring voltage at the ground surface resulting from ambient electrical currents within the earth. The measurement of NP can be used to characterize groundwater flow in karst terrane because electrical potential gradients are generated by the horizontal flow of water along fractures or conduits and the vertical infiltration of water into fractures or shafts. NP data from a site on the Mitchell Plain of southern Indiana, USA, revealed that NP data can be decomposed into three components: topographic effect, residual NP and noise. At this site, NP was inversely proportional to elevation, but the correlation varied with time. The topographic correction factor varied from –2.5 to –1.2 mV/m (NP change per unit elevation increase), with an average linear correlation coefficient (R) of 0.95. Because the site slopes toward an adjacent creek that is the local groundwater discharge zone, one possible explanation for this effect is a streaming-potential mechanism generated by groundwater movement toward the creek. The residual NP data revealed three negative anomalies at the survey area. Two of them coincide with sinkholes. A part of the third anomaly is coincident with a small valley, and concentrated infiltration does occur at this elevation in other valleys at the site, as evidenced by the existence of sinkholes. However, the dispersed, low-magnitude nature of the third anomaly does not prove the existence of concentrated groundwater recharge activity. Received: 18 March 1998 · Accepted: 27 April 1998     

Spatial management of geotechnical data for site selection

A ground characterisation exercise was conducted in order to identify sites suitable for a blast furnace at the Redcar steelworks in northeast England. The purpose of this exercise, conducted some 20 years after the investigations were carried out and a decade after the steelworks had been constructed, was to develop techniques for managing spatial data using Geographical Information Systems (GIS) and geostatistics.

Borehole details, from the factual site investigation reports, were entered into a database using a relational structure with linked files. The borehole location and lithology files were keyed in using codes for geological material, weathering grade, strength, borehole type and verification status. Laboratory and field test data were input by scanning tables in the ground investigation reports and processing the scanned images using optical character recognition software.

The XBase programming language was used to execute graphing, GIS and geostatistics programs; thereby enabling repetitive tasks to be carried out by the computer. This both saved time and reduced errors and inconsistencies that could result from repeated manual execution of the same program. Some simple spatial operations could also be carried out entirely from within dBase, an example being the use of moving-window statistics employing rectangular windows aligned parallel to the axes of the coordinate system to study local stationarity.

GIS have been found to offer an effective tool for analysing, querying and presenting spatial site investigation information. Geostatistical kriging and simulation techniques offer powerful spatial modelling tools for creating visualisations of the spatial variability of parameters measured at points during a ground investigation.

The Redcar case history demonstrates how a database management system (DBMS) can be used to manipulate site investigation data as part of a modelling and interpretation exercise. DBase source code allows others to follow the same path or build upon the foundations herein.     

Certain features in zonation of dispersion halos of Akhtal barite-polymetallic ore deposit,Northern Armenia

The most higly contaminated surface areas from cesium-137 fallout from the April 1986 accident at the Chernobyl' nuclear power station in Ukraine occur within the 30-km radius evacuation zone set up around the station, and an 80-km lobe extending to the west-southwest. Lower levels of contamination extend 300 km to the west of the power station. The deposition of this radioactive dust on the surface and the subsequent entombment of the damaged reactor effectively result in the de facto establishment of an above-ground nuclear waste storage site. This site is located on a thick sedimentary sequence of loose, mostly clastic deposits, with a shallow (generally 3-5 m) water table. The geology, the presence of surface water, a shallow water table, and leaky aquifers at depth make this an unfavorable environment for the long-term containment and storage of the radioactive debris. An understanding of the general geology and hydrology of the area is important to assess the environmental impact of this unintended waste storage site, and to evaluate the potential for radionuclide migration through the soil and rock and into subsurface aquifers and nearby rivers.     

Easy to Use Empirical Correlations for Liquefaction and No Liquefaction Conditions

Considerable damages during an earthquake (EQ) are the consequence of in situ soil losing its shear strength which is popularly known as liquefaction. A number of methodologies are available to quantify the safety of a site against liquefaction occurrence. Widely accepted recent methodologies follow iterative process making it cumbersome for the field engineer. In the present work, empirical correlations are proposed in accordance with widely accepted methodology, analysing the effect of various parameters such as overburden pressure, fines content (FC), factor of safety (FOS) etc. These proposed correlations are easy to use for the designers and the field engineers to determine the liquefaction potential of a site. Considering data from 207 global sites, proposed correlations are validated by comparing with standard methodology. Three different graphical validations are presented supporting that the results based on the proposed correlations are closely matching with the standard methodology. In case a site is found susceptible to liquefaction, so far no correlations are available to determine the shear strength required to be achieved after ground improvement which will ensure no liquefaction during future EQ. Proposed correlations in this work can also be used easily to determine improved shear strength required for a known FOS, FC and EQ magnitude (M) from ground improvement. Two flowcharts explaining the use of proposed correlations to determine FOS of a site and improvement shear strength required for a liquefied site from ground improvement respectively are developed in this work. Based on the second flowchart, determination of shear strength required from ground improvement are done for 45 random sites out of 207 liquefied sites during worldwide EQ in this work.     

Numerical modelling of compensation grouting above shallow tunnels

This paper describes the development of a numerical model for compensation grouting which is a useful technique for the protection of surface structures from the potentially damaging movements arising from tunnel construction. Pipes are inserted into the ground between the tunnel and the overlaying structure from an access shaft. Buildings on the surface are instrumented and movements are carefully monitored. Once the deformations exceed a certain Trigger Level, grout is injected into the ground to prevent damage. In the finite element model described here, compensation grouting is modelled by applying an internal pressure to zero‐thickness interface elements embedded in the mesh. An ‘observational algorithm’ is used, where the deformations of the surface are monitored and used to control the injection process. Example analyses of compensation grouting are given for three‐dimensional tunnel construction underneath a greenfield site. Different strategies are used to control the injection process and their effectiveness in preventing surface movement is assessed. The numerical model is shown to replicate general behaviour expected in the field and is capable of modelling the control of ground surface movements at a greenfield site. Copyright © 2005 John Wiley & Sons, Ltd.