Abstract
The Xinchang site in northwest China has been selected as a potential deep geological repository for spent high-level radioactive waste (HLW) due to its low-permeability granite bedrock. Predicting the characteristics of long-term flow systems is of priority for improving the suitability assessment for a HLW disposal repository. In this study, a three-dimensional saturated–unsaturated groundwater flow numerical model, using the software TOUGH3/EOS9 with an improved evaporation module, was developed. It was subsequently used to reveal the regional groundwater flow patterns in the bedrock fracture network and to provide scientific support for the safety assessment of the potential HLW disposal repository at the Xinchang site. The regional groundwater flow system associated with the Xinchang site can be characterized by three hydrogeological units. The groundwater at Xinchang site cannot flow into the Hexi Corridor due to the high surface topography and low-permeability F11 fault, but it flows southward to Yuanyang Ditch and then discharges from the east and south boundaries. Changes in precipitation, fault permeability and fault connectivity only cause fluctuation of the water table in local areas, and do not impact the regional groundwater flow patterns. Topography is the main factor determining regional groundwater flow. The proposed model is expected to serve as a powerful tool for studying groundwater flow patterns in fractured bedrock regions and provides the foundation for deeper analysis of nuclides migration in terms of distance and time in the Xinchang area.
Résumé
Le site de Xinchang, dans le nord-ouest de la Chine, a été choisi comme site potentiel de stockage en couches géologiques profondes de déchets hautement radioactifs (DHR) dans un substratum granitique fracturé et de faible perméabilité. La prévision des caractéristiques des systèmes d'écoulement à long terme est une priorité pour améliorer l'évaluation de l'adéquation d'un dépôt de déchets DHR. Dans cette étude, un modèle numérique tridimensionnel d'écoulement des eaux souterraines en zones saturée et non saturée, utilisant le logiciel TOUGH3/EOS9 avec un module d'évaporation amélioré, a été développé. Il a ensuite été utilisé pour révéler les schémas régionaux d'écoulement des eaux souterraines et pour apporter un soutien scientifique à l'évaluation de la sûreté du dépôt potentiel de DHR sur le site de Xinchang. Le système régional d'écoulement des eaux souterraines associé au site de Xinchang peut être caractérisé par trois unités hydrogéologiques. Les eaux souterraines du site de Xinchang ne peuvent pas s'écouler dans le corridor de Hexi en raison de la topographie de surface élevée et de la faille F11 peu perméable ; mais elles s'écoulent vers le sud jusqu'au fossé de Yuanyang, puis se déversent à partir des limites est et sud. Les changements des précipitations, de la perméabilité des failles et leur connectivité n'entraînent que des fluctuations du niveau piézométrique dans des zones locales et n'ont pas d'impact sur les schémas d'écoulement des eaux souterraines à l'échelle régionale. La topographie est le principal facteur déterminant l'écoulement régional des eaux souterraines. Le modèle proposé devrait constituer un outil puissant pour l'étude des schémas d'écoulement des eaux souterraines dans les régions à substratum rocheux fracturé et fournir les bases d'une analyse plus approfondie de la migration des nucléides en termes de distance et de temps dans la région de Xinchang.
Resumen
El emplazamiento de Xinchang, en el noroeste de China, ha sido seleccionado como posible depósito geológico profundo de residuos radiactivos de alta actividad ("HLW", por sus siglas en inglés) almacenados en un basamento de granito fracturado de baja permeabilidad. La predicción de las características de los sistemas de flujo a largo plazo es prioritaria para mejorar la evaluación de la idoneidad de un depósito de HLW. En este estudio, se desarrolló un modelo numérico tridimensional de flujo de agua subterránea saturada-no saturada, utilizando el software TOUGH3/EOS9 con un módulo de evaporación mejorado. Posteriormente se utilizó para determinar los patrones de flujo de las aguas subterráneas regionales y proporcionar apoyo científico a la evaluación de la seguridad del posible depósito de HLW en el emplazamiento de Xinchang. El sistema regional de flujo de aguas subterráneas asociado al emplazamiento de Xinchang puede caracterizarse por tres unidades hidrogeológicas. Las aguas subterráneas del emplazamiento de Xinchang no pueden fluir hacia el corredor de Hexi debido a la elevada topografía de la superficie y a la falla F11 de baja permeabilidad, pero fluyen hacia el sur hasta la zona de la Fosa de Yuanyang y luego se descargan por los límites este y sur. Los cambios en las precipitaciones, la permeabilidad y la conectividad de las fallas sólo provocan la fluctuación del nivel freático en zonas locales y no afectan a los patrones regionales de flujo de las aguas subterráneas. La topografía es el principal factor que determina el flujo regional de aguas subterráneas. Se prevé que el modelo propuesto sirva de potente herramienta para estudiar los patrones de flujo de aguas subterráneas en regiones de basamento fracturado y proporcione la base para un análisis más profundo de la migración de nucleidos en términos de distancia y tiempo en la zona de Xinchang.
摘要
中国西北地区新场场址已被选为低渗透性断裂花岗岩基岩中高放射性废物(HLW)的潜在深地质库。预测长期水流系统的特征对于提高HLW处置库的适应性评估至关重要。本研究开发了一个三维饱和-非饱和地下水流数值模型,采用改进蒸发模块的TOUGH3/EOS9软件。随后利用该模型揭示了新场场址区域地下水流模式,并为新场场址潜在HLW处置库的安全评估提供科学支持。与新场场址相关的区域地下水流系统可以划分为三个水文地质单元。由于地表地形高差和低渗透性的F11断层,新场场址的地下水无法流入河西走廊,但会向南流向鸳鸯沟,然后从东部和南部边界排出。降水量、断层渗透性和断层连通性的变化只会导致局部地区水位的波动,不会影响区域地下水流模式。地形是决定区域地下水流的主要因素。该模型预计将成为研究断裂岩基区域地下水流模式的有力工具,并为在新场场址以距离和时间为基础进行更深入的核素迁移分析奠定基础。
Abstrato
O local de Xinchang, no noroeste da China, foi selecionado como um potencial repositório geológico profundo para resíduos radioativos de alto nível (RRAN) em leito rochoso de granito fraturado de baixa permeabilidade. Prever as características dos sistemas de fluxo de longo prazo é prioritário para melhorar a avaliação de adequação para um depósito de descarte de RRAN. Neste estudo, foi desenvolvido um modelo numérico tridimensional de fluxo de águas subterrâneas saturadas-insaturadas, usando o software TOUGH3/EOS9 com um módulo de evaporação aprimorado. Posteriormente, foi usado para revelar os padrões regionais de fluxo de águas subterrâneas e fornecer suporte científico para a avaliação de segurança do potencial repositório de descarte de HLW no local de Xinchang. O sistema regional de fluxo de águas subterrâneas associado ao local de Xinchang pode ser caracterizado por três unidades hidrogeológicas. A água subterrânea no local de Xinchang não pode fluir para o Corredor Hexi devido à alta topografia da superfície e à falha F11 de baixa permeabilidade, mas flui para o sul até a vala de Yuanyang e depois descarrega dos limites leste e sul. Mudanças na precipitação, permeabilidade de falhas e conectividade de falhas causam apenas flutuação do lençol freático em áreas locais e não afetam os padrões regionais de fluxo das águas subterrâneas. A topografia é o principal fator determinante do fluxo regional de águas subterrâneas. Espera-se que o modelo proposto sirva como uma ferramenta poderosa para estudar os padrões de fluxo de águas subterrâneas em regiões rochosas fraturadas e forneça a base para uma análise mais profunda da migração de nuclídeos em termos de distância e tempo na área de Xinchang.
Similar content being viewed by others
References
Albers BMC, Molson JW, Bense VF (2020) Parameter sensitivity analysis of a two-dimensional cryo-hydrogeological numerical model of degrading permafrost near Umiujaq (Nunavik, Canada). Hydrogeol J 28:905–919. https://doi.org/10.1007/s10040-020-02112-2
Barenblatt GI, Zheltov IP, Kochina IN (1960) Basic concepts in the theory of seepage of homogeneous liquids in fissured rocks. J Appl Math Mech 24(5):1286–1303
Belcher WR, Sweetkind DS (2010) Death Valley regional groundwater flow system, Nevada and California-Hydrogeologic framework and transient groundwater flow model. Professional Paper 1711. https://doi.org/10.3133/pp1711
Blessent D, Therrien R, Gable CW (2011) Large-scale numerical simulation of groundwater flow and solute transport in discretely-fractured crystalline bedrock. Adv Water Resour 34(12):1539–1552. https://doi.org/10.1016/j.advwatres.2011.09.008
Bodin J, Delay F, Marsily GD (2003a) Solute transport in a single fracture with negligible matrix permeability: 1. fundamental mechanisms. Hydrogeol J 11(4):418–433. https://doi.org/10.1007/s10040-003-0268-2
Bodin J, Delay F, Marsily GD (2003b) Solute transport in a single fracture with negligible matrix permeability: 2. mathematical formalism. Hydrogeol J 11(4):434–454. https://doi.org/10.1007/s10040-003-0269-1
Bodvarsson GS, Wu YS, Zhang K (2003) Development of discrete flow paths in unsaturated fractures at Yucca Mountain. J Contam Hydrol 62–63:23–42. https://doi.org/10.1016/S0169-7722(02)00177-8
Cao X, Hu L, Wang J, Wang J (2017a) Radionuclide transport model for risk evaluation of high-level radioactive waste in Northwestern China. Hum Ecol Risk Assess 23(8):2017–2032. https://doi.org/10.1080/10807039.2017.1361811
Cao X, Hu L, Wang J, Wang J (2017b) Regional groundwater flow assessment in a prospective high-level radioactive waste repository of China. Water 9(7):551. https://doi.org/10.3390/w9070551
Cao X, Hou D, Hu L (2020) Numerical simulation of regional groundwater flow in the Beishan area of Gansu. Hydrogeol Eng Geol 47(2):9–16. https://doi.org/10.16030/j.cnki.issn.1000-3665.201907022 (in Chinese)
Cvetkovic V, Painter S, Outters N, Selroos JO (2004) Stochastic simulation of radionuclide migration in discretely fractured rock near the Äspö Hard Rock Laboratory. Water Resour Res 40:W02404. https://doi.org/10.1029/2003WR002655
Dong Y, Li G, Li M (2009) Numerical modeling of the regional ground water flow in Beishan area Gansu Province. Chin Sci Bull 54(17):3112–3115. https://doi.org/10.1007/s11434-009-0344-7
Dong Y, Fu Y, Wang L, Wang J, Zhang Q, Zong Z, Zhou Z (2022) Regional groundwater flow pattern in Beishan, Hexi Corridor and Qilian Mountain. Bull Geol Sci Technol 41(1):79–89
Fu Y, Dong Y, Wang L, Bour O, Klepikova MV, Zong Z, Xu Z, Zhou Z (2022) Characteristics of hydraulic conductivity in mountain block systems and its effects on mountain block recharge: Insights from field investigation and numerical modeling. J Hydrol. 612:128184. https://doi.org/10.1016/j.jhydrol.2022.128184
Hartley L, Appleyard P, Baxter S, Hoek J, Joyce S, Mosley K, Williams T, Fox A, Cottrell M, Pointe PL, Gehör S, Darcel C, Goc RL, Aaltonen I, Vanhanarkaus O, Löfman J, Poteri A (2018) Discrete fracture network modelling (Version 3) in support of Olkiluoto site description 2018. Working report 2017-32. Posiva Oy, Eurajoki
Haukwa CB, Wu YS, Bodvarsson GS (2003) Modeling thermal-hydrological response of the unsaturated zone at Yucca Mountain, Nevada, to thermal load at a potential repository. J Contam Hydrol 62–63:529–552. https://doi.org/10.1016/S0169-7722(02)00188-2
He Y, Jiang X, Wang N, Zhang S, Ning T, Zhao Y, Hu Y (2019) Changes in mountainous runoff in three inland river basins in the arid Hexi Corridor, China, and its influencing factors. Sustain Cities Soc 50:101703. https://doi.org/10.1016/j.scs.2019.101703
Huang J, Ma J, Guan X, Li Y, He Y (2019) Progress in semi-arid climate change studies in China. Adv Atmos Sci 36(9):922–937. https://doi.org/10.1007/s00376-018-8200-9
James SC, Doherty JE, Eddebbarh AA (2009) Practical postcalibration uncertainty analysis: Yucca Mountain, Nevada. Ground Water 47(6):851–869. https://doi.org/10.1111/j.1745-6584.2009.00626.x
Joyce S, Hartley L, Applegate D, Hoek J, Jackson P (2014) Multi-scale groundwater flow modeling during temperate climate conditions for the safety assessment of the proposed high-level nuclear waste repository site at Forsmark, Sweden. Hydrogeol J 22(6):1233–1249. https://doi.org/10.1007/s10040-014-1165-6
Jung Y, Pau GSH, Finsterle S, Doughty CA (2018) TOUGH3 User’s Guide, Version 1.0. Lawrence Berkeley National Lab. (LBNL), Berkeley. https://doi.org/10.2172/1461175
Karimi FM, Gong B, Durlofsky LJ (2006) Generation of coarse-scale continuum flow models from detailed fracture characterizations. Water Resour Res 42(10):W10423. https://doi.org/10.1029/2006WR005015
Kelkar S, Ding M, Chu S, Robinson BA, Arnold B, Meijer A, Eddebbarh AA (2010) Modeling solute transport through saturated zone ground water at 10 km scale: example from the Yucca Mountain license application. J Contam Hydrol 117:7–25. https://doi.org/10.1016/j.jconhyd.2010.05.003
Kelkar S, Srinivasan G, Robinson BA, Roback R, Viswanathan H, Rehfeldt K, Tucci P (2013) Breakthrough of contaminant plumes in saturated volcanic rock: implications from the Yucca Mountain site. Geofluids 13:273–282. https://doi.org/10.1111/gfl.12035
Lan Y, Xiao H, Hu X, Ding H, Zou S, La C, Song J (2012) Study of temperature and precipitation change in upstream mountain area of the Hexi inland river basin since 1960s. Sci Cold Arid Reg 4(6):0522–0535. https://doi.org/10.3724/sp.J.1226.2012.00522
Li J, Su R, Zhou Z, Guo Y, Ji R, Zhang M, Wang H (2020) Estimation of shallow groundwater recharge using the environmental tracer chloride method in Beishan area. Arid Land Geol 43(1):135–143. https://doi.org/10.12118/j.issn.1000-6060.2020.01.16 (in Chinese)
Li L, Zhou Z, Shao J, Cui Y, Zhao J (2021) Advances in groundwater numerical simulation in deep geological disposal of high-level radioactive waste. Hydrogeol Eng Geol 48(6):13–23. https://doi.org/10.16030/j.cnki.issn.1000-3665.202010061 (in Chinese)
Li L, Zhang Q, Zhou Z, Cui Y, Shao J, Zhao Y (2022) Groundwater circulation patterns in bedrock aquifers from a pre-selected area of high-level radioactive waste repository based on two-dimensional numerical simulation. J Hydrol 610:127849. https://doi.org/10.1016/j.jhydrol.2022.127849
Liu HH, Haukwa CB, Ahlers CF, Bodvarsson GS, Flint AL, Guertal WB (2003) Modeling flow and transport in unsaturated fractured rock: an evaluation of the continuum approach. J Contam Hydrol 62–63:173–188. https://doi.org/10.1016/S0169-7722(02)00170-5
Liu Y, Wu C, Jia R, Huang J (2018) An overview of the influence of atmospheric circulation on the climate in arid and semi-arid region of central and east Asia. Sci China Earth Sci 61(9):1183–1194. https://doi.org/10.1007/s11430-017-9202-1
Löfman J, Mészáros F (2005) Simulation of hydraulic disturbances caused by the underground rock characterization facility in Olkiluoto. Posiva Report 2005-08. Posiva Oy, Olkiluoto
Long JCS, Gilmour P, Witherspoon PA (1985) A model for steady fluid flow in random three-dimensional networks of disc-shaped fractures. Water Resour Res 21(8):1105–1115
Mckinley IG, Alexander WR, Blaser PC (2007) Development of geological disposal concepts. Radioactivity in the Environment 9(6):41–76. https://doi.org/10.1016/S1569-4860(06)09003-6
Odén M, Niemi A, Tsang CF, Ohman J (2008) Regional channelized transport in fractured media with matrix diffusion and linear sorption. Water Resour Res 44(2):W02421. https://doi.org/10.1029/2006WR005632
Posiva Oy (2012a) Olkiluoto site description 2011. Posiva Report 2011-02. Posiva Oy, Eurajoki
Posiva Oy (2012b) Safety case for the disposal of spent nuclear fuel at Olkiluoto - assessment of radionuclide release scenarios for the repository system 2012. Posiva Report 2012-09. Posiva Oy, Eurajoki
Pulkkanen VM, Nordman H (2011) Effects of bedrock fractures on radionuclide transport near a vertical deposition hole for spent nuclear fuel. Posiva Report 2011-03. Posiva Oy, Eurajoki
Rechard RP, Arnold BW, Robinson BA, Houseworth JE (2014) Transport modeling in performance assessments for the Yucca Mountain disposal system for spent nuclear fuel and high-level radioactive waste. Reliab Eng Syst Saf 122:189–206. https://doi.org/10.1016/j.ress.2013.06.031
Richards LA (1931) Capillary conduction of liquids through porous mediums. Physics 1(5):318–333. https://doi.org/10.1063/1.1745010
Robinson BA, Chu S, Lu Z (2012) Simulation of radionuclide transport through unsaturated, fractured rock: application to Yucca Mountain, Nevada. Vadose Zone J 11(4). https://doi.org/10.2136/vzj2011.0142
Rui X (2004) Hydrology theory. China Water & Power Press, Beijing
Schwartz MO (2012) Modelling radionuclide transport in large fractured-media systems: the example of Forsmark, Sweden. Hydrogeol J 20(4):673–687. https://doi.org/10.1007/s10040-012-0837-3
Selroos J-O, Follin S (2014) Overview of hydrogeological safety assessment modeling conducted for the proposed high-level nuclear waste repository site at Forsmark, Sweden. Hydrogeol J 22(6):1229–1232. https://doi.org/10.1007/s10040-014-1163-8
SNL (Sandia National Laboratories) (2007) Saturated zone site-scale flow model. MDL-NBS-HS-000011 REV 03. SNL, Las Vegas
SNL (2008) Site-scale saturated zone transport. MDL-NBS-HS-000010 REV 03 AD 01. SNL, Las Vegas
Snow DT (1969) Anisotropic permeability of fractured media. Water Resour Res 5(6):1273–1289
Wang L (2015) Multi-scale groundwater numerical simulation study of regional-basin-site in Gansu Beishan area. University of Chinese Academy of Sciences, Beijing (in Chinese)
Wang J (2019) Progress of geological disposal of high-level radioactive waste in China in the 21st century. At Energy Sci Technol 53(10):2072–2082. https://doi.org/10.7538/yzk.2019.53.10.2072 (in Chinese)
Wang J, Su R, Chen W, Guo Y, Jin Y, Wen Z, Liu Y (2006) Deep geological disposal of high-level radioactive waste in China. Chin J Rock Mech Eng 25(4):649–658. https://doi.org/10.3321/j.issn:1000-6915.2006.04.001 (in Chinese)
Wang Y, Zhang Q, Wang S, Wang J, Yao Y (2017) Characteristics of agro-meteorological disasters and their risk in Gansu Province against the background of climate change. Nat Hazards 89(2):899–921. https://doi.org/10.1007/s11069-017-2999-8
Wang J, Chen L, Su R, Zhao X (2018) The Beishan underground research laboratory for geological disposal of high-level radioactive waste in China: planning, site selection, site characterization and in situ tests. J Rock Mech Geotech Eng 10(3):411–435
Ye M, Pan F, Wu YS, Hu BX, Shirley C, Yu Z (2007) Assessment of radionuclide transport uncertainty in the unsaturated zone of Yucca Mountain. Adv Water Resour 30:118–134. https://doi.org/10.1016/j.advwatres.2006.03.005
Zhang K, Wu YS, Houseworth JE (2006) Sensitivity analysis of hydrological parameters in modeling flow and transport in the unsaturated zone of Yucca Mountain, Nevada, USA. Hydrogeol J 14:1599–1619. https://doi.org/10.1007/s10040-006-0055-y
Zimmerman RW, Chen G, Hadgu T, Bodvarsson GS (1993) A numerical dual-porosity model with semianalytical treatment of fracture/matrix flow. Water Resour Res 29(7):2127–2137
Zyvoloski G, Kwicklis E, Eddebbarh AA, Arnold B, Faunt C, Robinson BA (2003) The site-scale saturated zone flow model for Yucca Mountain: calibration of different conceptual models and their impact on flow paths. J Contam Hydrol 62–63:731–750. https://doi.org/10.1016/S0169-7722(02)00190-0
Acknowledgements
Thanks to Beijing Research Institute of Uranium Geology for the in-situ hydraulic tests data. The constructive comments and suggestions from the anonymous reviewers are greatly appreciated and have helped improve the manuscript.
Funding
This research was sponsored by the special project for Decommissioning of Nuclear Facilities and Radioactive Waste Treatment by the State Administration of Science, Technology and Industry for National Defense under Grant Number: 2017-1405.
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflict of Interest
We declare that there are no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Below is the link to the electronic supplementary material.
Rights and permissions
Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
About this article
Cite this article
Li, L., Zhang, Q., Zhou, Z. et al. Regional groundwater flow pattern in fractured bedrock with a high-level radioactive waste repository: numerical simulation of the Xinchang site, northwestern China. Hydrogeol J 31, 2159–2176 (2023). https://doi.org/10.1007/s10040-023-02714-6
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10040-023-02714-6