Review of groundwater numerical simulation for deep geological disposal of high-level radioactive waste
-
摘要: 开展地下水数值模拟研究是高放废物处置场地安全评价的重要组成部分，然而深地质处置介质类型的复杂性、基岩深部资料的相对匮乏性导致模拟结果存在不确定性，如何刻画深部地下水动力场并评估可能引起的风险已成为高放废物处置安全评价中重点关注的问题。在大量文献调研的基础上，综述了世界典型国家高放废物深地质处置场地的地下水数值模拟与不确定性分析应用，并归纳总结该领域研究经验，得到以下认识：（1）深地质处置场深部构造、裂隙的发育与展布决定了地下水循环条件，探究适用于基岩裂隙地区新的水文地质试验方法是提高地下水数值模型仿真性的基础；（2）不同尺度模型融合是解决深地质处置地下水模拟的有效技术方法，区域尺度多采用等效连续介质法，场地尺度使用等效连续多孔介质和离散裂隙网络耦合模型，处置库尺度使用离散裂隙网络方法，其次需重点关注未来大时间尺度下放射性核素在地质体中的迁移转化规律，模拟预测场址区域地下水环境长期循环演变对核素迁移的潜在影响；（3）考虑到不同的处置层主岩岩性以及在多介质中发生的THMC（温度场—渗流场—应力场—化学场）过程，目前国内外常用的地下水模拟软件有：Porflow、Modflow、GMS及MT3DMS等用于模拟孔隙或等效连续介质，Connectflow、Feflow及FracMan等用于模拟地下水和核素在结晶岩、花岗岩等裂隙中的迁移，TOUGH系列软件主要应用于双重介质的水流、溶质及热运移模拟；（4）指导开展有针对性的模型和参数的不确定性分析工作，减少投入工作量，提高模型精度，并可针对处置库长期演变、废物罐失效、极端降雨等多情景预测模拟，为处置库安全评价及设计提供基础数据支撑；（5）针对我国深地质处置地下水数值模拟研究现状，下一步应加强区域地质、水文地质、裂隙测量以及现场试验等相关的调查及监测工作，多介质耦合、多场耦合模拟及不确定性分析研究将会是未来的研究重点。Abstract: Groundwater numerical simulation is an important part of the safety assessment for high-level radioactive waste (HLW) disposal sites. The complexity of deep geological disposal media and the relative lack of deep bedrock data lead to uncertainties in the simulation results. How to characterize the deep groundwater dynamic field and evaluate the possible risks has become a key issue in the safety assessment of HLW disposal. Based on a lot of literature investigation, this paper reviews the application cases of groundwater numerical simulation and uncertainty analysis of HLW deep geological disposal sites in typical countries, and summarizes the research experience in this field. The results show that (1) the structure and fracture development and distribution of deep geologic repositories determine the groundwater circulation conditions, and exploration of the new hydrogeological test method suitable for fractured bedrock areas is the basis of improving the accuracy of groundwater numerical simulation. (2) The integration of different scale models is an effective technical method to solve the groundwater simulation for deep geological disposal. The equivalent continuum method is mostly used at a regional scale, the coupling model of equivalent continuous porous medium and discrete fracture network is used at a site scale, and the discrete fracture network method is used at a repository scale. Furthermore, it is necessary to pay more attention to the migration and transformation of radionuclides in geological formations, to simulate and predict the potential impact of the long-term groundwater environment evolution on radionuclide migration in the site area. (3) Considering different host rocks in the disposal layers and the thermal-hydrogeological-mechanical-chemical process occurring in multiple media, the commonly used software for groundwater numerical simulation in HLW deep geological disposal includes Porflow, Modflow, GMS and MT3DMS, which are used to deduct pores or fractures into equivalent continuum, and groundwater flow and nuclides migration in fractures of crystalline rock and granite can be generalized by using Connectflow, Feflow and Fracman, whereas TOUGH series is mainly used to simulate flow, solute and heat transport in dual media. (4) Targeted uncertainty analysis of the models and parameters should be carried out to reasonably reduce workload and improve the model accuracy. Moreover, the prediction simulation of repository long-term evolution, waste canister failure, extreme rainfall and other scenarios can provide basic data support for the safety assessment and design of the repository. (5) In view of the current researches on groundwater numerical simulation of HLW deep geological disposal in China, it is necessary to strengthen the investigation and monitoring of geology, hydrogeology, fissure measurement and field test in the next step, and the multi-medium coupling model, multi-field coupling simulation and uncertainty analysis will be the focus of future researches.
-
地质遗迹是由各种内外地质作用形成、发展并遗留下来的珍贵的、不可再生的地质现象，不仅为人类探索地质历史时期地球古气候、古环境、生命演化等过程提供可靠证据，同时也造就了绚丽多彩的秀美山川、大江大河、湿地平原等自然景观[1-2]。由成因上相同或相关，空间上相连的地质遗迹构成的区域或范围即为地质遗迹点[1]。
人类活动已经成为重要地质营力，人-地相互作用影响也越来越明显[3-4]，尤其是对于区域地质地貌（包括地质遗迹）的改造亦成为城镇化发展过程的表征之一。作为中国城镇化高水平发展代表的京津冀城市群之一的天津，其地质遗迹资源大多位于城区周边人类生产、生活较频繁的地区，受人为干扰较大。对地质遗迹分布边界和演化规律的调查和研究可以为科学划定地质遗迹分区和保护范围、妥善处理地质遗迹保护与开发的关系、合理规划城市产业布局、有效缓解人类活动和地质遗迹保护之间的冲突和矛盾等提供基础支撑。
天津市地质遗迹调查工作尚处于早期，还未出台正式的规范性文件。由于对地质遗迹点概念认识上存在差异，在进行地质遗迹调查时，常常局限于蓟州北部基岩出露区的地质公园或保护区内，调查和评价对象一般为某一地质遗迹点的特征点，导致存在地质遗迹分类方案不够规范、地质遗迹特征点数量多但却难以反映天津市地质遗迹整体状况以及缺乏对地质遗迹的整体评价的问题。本次地质遗迹调查工作不同于以往，一是作为行业标准《地质遗迹调查规范》（DZ/T 0303—2017）颁布实施以后的第一次规范性调查，对天津市地质遗迹资源进行了全面调查和梳理，查清了重要地质遗迹特征、规模、分布范围、演化历史等[5]，采用了规范的地质遗迹分类方案；二是调查研究的空间尺度由地质公园和保护区，上升为省、市乃至国家层面，对观察点和特征点进行筛选、梳理等，完整、科学、规范、全面地对地质遗迹点进行认定、描述、分类；三是工作过程中采用了定性评价和定量评价相结合的综合分级评价方法，使地质遗迹评价结果更可靠。
1. 研究区概况
天津市总面积1.37万 km2，陆域面积约为1.19万 km2，海域面积约为0.18万km2。北部山区属燕山山地，分布在蓟州区北部，面积约占陆域面积的7%，出露新太古界、中新元古界及下古生界地层，其中，以中新元古界为主，新太古界及下古生界仅零星分布；南部平原属华北平原的一部分，为堆积平原，分布在山区以南广大地区，约占陆域面积的93%，被第四系覆盖，第四系之下发育中新元古界至新近系地层；东南部濒临渤海湾，属海积低平原,海岸线长度约360 km，主要由滨海盐滩、潟湖洼地、沼泽和潮滩组成，以淤泥质黏土为主。
2. 重要地质遗迹分类及其分布
在充分搜集天津市地质遗迹有关资料基础上，初步筛选出具有重要科学、科普和观赏价值的重要地质遗迹作为研究对象。通过野外调查和综合研究分析，对初选的地质遗迹进行分类和综合分级评价。由于地质遗迹在漫长的地质历史时期，经过多期复杂的地质变化，因此地质遗迹具有改造成因多样、自身具备多重科学价值属性的特点，这都将决定地质遗迹的最终分类和综合分级评价结果。
经调查研究，天津市重要地质遗迹点为15处，根据《地质遗迹调查规范》（DZ/T 0303—2017）中地质遗迹分类方案和最能反映地质遗迹科学价值的成因等，分为3大类，7类和10亚类（表1）。
表 1 天津市重要地质遗迹点一览表Table 1. List of important geological sites in Tianjin大类 类 亚类 遗迹名称 出露规模/km2 基础地
质大类地层剖面 层型 蓟县剖面（即中新元
古界标准剖面）8.9 重要岩矿石产地 典型矿床类露头 锰方硼石矿床 − 地貌景
观大类水体地貌 湿地 七里海湿地 68 北大港湿地 349 团泊洼湿地 60 大黄堡湿地 112 海岸地貌 海积地貌 渤海湾贝壳堤 − 渤海湾牡蛎礁 − 天津沿海滩涂 370 岩土体
地貌侵入岩地貌 盘山花岗岩地貌 12 碎屑岩地貌 蓟州区石英砂岩峰
林峡谷地貌24 碳酸盐岩地貌 九龙山碳酸盐岩地貌 5 构造地貌 飞来峰 府君山飞来峰地貌 3 断层崖 黄崖关断崖地貌 12 地质灾
害大类地质灾害
遗迹地面沉降 天津市地面沉降 5150 天津市地质遗迹涵盖地层剖面、重要岩矿石产地、岩土体地貌、构造地貌、水体地貌、海岸地貌、地质灾害遗迹，其中前四者分布在天津市北部基岩出露区，后三者分布在天津市南部第四系深覆盖区。
3. 重要地质遗迹
3.1 地层剖面
蓟县剖面（即蓟县中、新元古界标准地层剖面简称）位于蓟州区北部山区，起点从常州沟村到蓟州城区附近，绵延十数公里，展布在燕山南麓，出露面积约9 km2。蓟县剖面地层厚度约9200 m，底界与太古界分界，见图1（a），距今16.5亿 a[6-7]，顶界与寒武系昌平组分界，见图1（b），距今7.8亿 a[8]。地层划分为4个系：长城系、蓟县系、待建系、青白口系；12个组：常州沟组、串岭沟组、团山子组、大红峪组、高于庄组、杨庄组、雾迷山组、洪水庄组、铁岭组、下马岭组、龙山组、景儿峪组。除下马岭组和龙山组之外，其他各组均为层型所在地。地层岩性主要为碳酸盐岩，占岩石种类的73.53%，其次是碎屑岩，占17.6%，其余为黏土质岩石、火山岩和火山碎屑岩，因此，蓟县剖面属于碳酸盐岩相。蓟县剖面的古生物化石，以页岩相微古植物化石、燧石相微化石、炭质宏观藻类化石和隐藻化石中的叠层石资源最为丰富，尤其是宏观藻类化石群的发现，对研究生命起源具有重大意义。现已建立国家级自然保护区。
3.2 重要岩矿石产地
锰方硼石是一种很少见的无水氯硼酸盐矿物，分子式为Mn3B7O13Cl，于1962年在美国德克萨斯州钱伯斯镇发现，并命名为钱伯斯矿[9]。1971年，在天津蓟县发现金伯斯矿床，并翻译为锰方硼石矿。锰方硼石矿床位于蓟州区穿芳峪镇东水厂附近，矿层赋存于蓟县系高于庄组二段泥晶白云岩中，矿体成透镜状断续成层分布，东西延长达30 km。矿层中致密块状富矿体呈扁豆状、饼状、串珠状、筒状、散点状等断续排列，矿体走向280°~310°，倾向SW，倾角40°~60°，由于该矿极易风化，鲜见出露。目前，该矿床处于天津蓟县国家地质公园内，是天津市政府明令禁止开采的重点保护矿种。
3.3 水体地貌
天津市地处北运河、永定河、大清河、子牙河、南运河5条河流下游，河道纵横交错、支流众多，地表水资源丰富。全新世时期，天津市经历了多期海侵、海退事件，海侵最西处直达白洋淀；随着全球气候逐渐变冷，渤海逐渐退去，海水入侵逐渐减少，在渤海湾北岸形成大量的潟湖、沼泽等，在大气降水与河水注入的混合作用下，潟湖的盐度逐渐降低，逐渐成为了湿地。在没有人类干扰的情况下，湿地演化比较稳定，自然更替相对比较缓慢，但是在自然环境、人类活动和政策因素驱动下，天津市湿地类型有加速由自然湿地向人工湿地转化、空间上有向东南部转移的趋势[10-11]。
天津市湿地主要包括近海及海岸湿地、河流湿地、湖泊湿地、沼泽和沼泽化草甸湿地、人工湿地5大类型，面积约占天津总面积的20%，居北方省市第一，较重要的为七里海、北大港、大黄堡、团泊洼湿地等，其中七里海湿地最具地学研究价值和观赏价值。七里海湿地位于宁河区境内，面积约68 km2，为典型的古潟湖遗迹，形状呈凸字形，是由堤岸、贝壳、碎屑沙等组成的退海地，潮白新河纵穿七里海，将七其分隔为东西两海，东海为水库和苇地、西海为苇地。七里海湿地处于华北平原强烈沉降带东北部，新生代以来持续下沉，地貌上介于黄河海河三角洲与滦河三角洲之间相对低洼区。全新世时期，随着最后一次海退，七里海地区变成只能高潮水到达的潟湖，约隋唐时期，逐渐脱离海水影响，最终演化成淡水湖泊湿地。
3.4 海岸地貌
3.4.1 贝壳堤
渤海湾贝壳堤（图2）是由波浪或高潮水将贝壳和贝壳碎片从浅海泥砂中簸选出来，在海边的泥滩上逐渐堆积加高，形成上凸的堤状堆积体，一般埋藏在地下0~10 m，堤宽数十米，高数米，长数十公里，极少露头，目前研究发现共有4~6道贝壳堤沿渤海湾西北侧[12]，与海岸线平行分布，从滨海新区一直延伸到河北省黄骅市。目前，贝壳堤已处于国家级古海岸与湿地自然保护区内。渤海湾贝壳堤形成于7000 a 前，是由海生贝类动物在海潮推动下，逐渐堆积而成的古渤海岸线的标志。物质组成主要为海相软体动物双壳类的贝壳和碎片、少量的腹足类壳，以及他们的碎屑。
![]()
图 2 古林古海岸遗迹博物馆原地保存的贝壳堤Figure 2. Shell ridges preserved in situ in Gulin Ancient Coastal Heritage Museum3.4.2 牡蛎礁滩
渤海湾的西北岸，潮白河两岸及蓟运河以西，自早全新世最大海侵以来, 发育了规模巨大、数量众多的埋藏牡蛎礁体，往往形成于略高于自然地坪的垄岗状高地——俗称牡蛎（礁）岭, 是该地区重要地质现象[13-14]。渤海湾牡蛎礁在宁河区境内多处分布，是深埋于地下3~7 m不等的牡蛎遗骸的堆积体，形成于距今2000~7000 a，由长重蛎、近江重蛎组成，厚达5 m，据《天津古海岸与湿地国家级自然保护区保护对象调查与研究》[15]测算天津市牡蛎礁储量约为2.4×107 m3。牡蛎礁滩内部形态一般以直立与倒伏为主，气候温暖，牡蛎正常生长，则以直立矗立生长；气候变冷，海水退却，牡蛎大面积死亡、倒伏，则形成倒伏状，这也揭示了气候变化的规律，为研究古气候、古环境等提供证据。目前，牡蛎礁滩处于国家级古海岸与湿地自然保护区内。
3.4.3 天津沿海滩涂
海岸滩涂（简称海涂）位于陆地与海洋之间狭长的潮间地带，通常指海岸线至理论深度基准面——零米线间低潮时出露的滩地。天津市渤海湾海岸带位于渤海湾西—西北岸，海岸带宽广低平，形态单一，呈弧形带状分布，为典型的粉砂淤泥质平原海岸，由于人工填海、防潮堤的普遍修建，自然发展的海岸线已面目全非，补给泥沙量不足，导致淤泥质海岸处于不断侵蚀后退状态，并逐渐演化为砂质海岸[16]。目前，沿海滩涂尚未采取有效保护措施。
3.5 岩土体地貌和构造地貌
天津市北部地区主要出露地层为中新元古界地层，其中，北部的常州沟组石英砂岩在地质构造及流水侵蚀下形成石英砂岩峰林峡谷地貌，面积约24 km2，如刀砍斧凿、鬼斧神工般秀美；中部发育的高于庄组灰质白云岩构成的九龙山峰丛地貌景观，面积约5 km2，形成典型的低矮山峰、峰丛、石芽等地貌单元；蓟州北部与兴隆交界处在黄崖关断裂带的影响下形成了黄崖关断崖地貌，面积约12 km2，与盘桓而上的长城相映成辉；西部出露印支期火山活动形成的盘山侵入岩地貌，面积约12 km2，岩体节理裂隙发育，球形风化强烈，形成怪石嶙峋的石蛋山地貌；南部蓟州城区附近推覆构造形成的府君山飞来峰地貌，面积约3 km2，造成此地貌的地质活动被我国著名地质学家孙云铸命名为“蓟县运动”。目前，以上地质遗迹资源处于天津蓟县国家地质公园范围内。
3.6 地质灾害遗迹
天津市南部平原区松散沉积层巨厚，地下水开采集中，地面沉降地质灾害发育，是天津市最主要的地质灾害，也是全国地面沉降地质灾害主要代表区域[17]。目前，沉降区基本覆盖天津市南部平原区，呈不规则椭圆形，最大累计沉降量为3 m多，主要形成汉沽、杨村、西青—静海—津南区3个沉降中心。由于自然因素与人为因素，导致新近系、古近系、第四系地层产生沉降，继而造成地面高程损失（图3）、降低河道行洪和排泄能力、加重风暴潮危害、威胁轨道交通安全等一系列问题，损失大，影响持久，但随着控沉工作力度加强，这一趋势将放缓。
![]()
图 3 地面沉降导致的管道一侧抬升Figure 3. Lifting of one side of the pipeline caused by land subsidence4. 地质遗迹综合分级评价
本次地质遗迹综合分级评价的目的是为了建立全国地质遗迹保护名录，为地质遗迹保护区划、分级保护和开发提供政策依据。重要地质遗迹综合分级评价结果可根据地质遗迹的科学价值、观赏价值和保护价值分为世界级、国家级、省级3个等级。
地质遗迹评价方法众多，评价结果也有所不同，总体来说可以分为定性评价和定量评价。定性评价是大多数习惯采用的地质遗迹评价方法，但评价结果往往带有一定的主观性。定量评价一定程度可弥补定性评价的主观性缺点，但各评价指标的模糊得分仍然离不开定性分析。因此，本文采用定性评价和定量评价相结合的方式进行综合分级评价。
4.1 地质遗迹定性评价
4.1.1 科研价值
天津市地质遗迹资源科研价值主要体现在地层学、矿床学、生态地质学、第四纪地质学等学科上。如蓟县剖面对区域地层对比、生命的形成及进化、地球的演变以及古地理环境恢复等方面具有重大的科学价值，尤其是随着锆石U-Pb同位素测年技术的进步，蓟县剖面相关地层首次获得直接定年，为重新厘定年代格架、对比等，提供了新的、精确的年代学锚点[18-19]。锰方硼石是一种罕见的锰氯硼酸盐矿物，对该矿物的矿物学特征和成因研究填补了我国乃至世界锰方硼石矿物学和矿床学研究的空白[20]。渤海西北岸的贝壳堤、牡蛎礁、海岸滩涂及湿地等，记录了过去7000 a 间“沧海变桑田”过程中的数道古海岸线，是海陆变迁的重要产物和证据，众多国内外相关专业学者在此开展研究、调查工作，相继在国内外学术期刊上发表了众多关于海洋、地质、地理等系统海岸演变、古气候、湿地生态等学科的研究成果和文章。
4.1.2 观赏价值
天津市南部平原区湿地是全新世以来海侵、海退地质事件形成的独特的景观资源，这里大面积生长的芦苇，在天津乃至华北地区均属罕见，动植物资源丰富，拥有良好的湿地生态系统。每年大量候鸟在此汇集，其中包括多种珍稀濒危鸟种，是鸟类南北迁徙的中转站。北部地质遗迹出露区，位于燕山山脉南麓，自然风光优美，人文历史悠久，地质构造形成的峰林峡谷地貌、断层崖和花岗岩山峰与皇家文化、宗教文化和长城等历史文化遗迹交相呼应，形成了人文古迹与青山秀水相融合的景观特色，极具观赏价值。
4.1.3 保护价值
蓟县剖面从典型性、自然性、系统性和完整性上，均属世界一流，真实地记录着地球演化距今约16.5～7.8亿 a 间的地质演化史，赋存着能反映当时古地理、古气候、古生物、古构造、古地磁等的大量自然信息以及多种金属、非金属矿产资源。锰方硼石矿，在世界上分布稀少，只在天津市蓟县形成了世界上唯一的具有开发利用价值的锰方硼石矿床[21]，是一种重要的化工原料，在电子工业和原子能工业上都有重要应用价值，开发应用前景巨大，政府已明文禁止开采。渤海湾西北岸的贝壳堤、牡蛎礁和湿地是古海岸线的标志，对研究古海洋、古气候、古环境等科研领域和调节区域气候、净化水质、保护生物多样性具有重要作用，在国际上占有重要位置。
4.2 地质遗迹定量评定
4.2.1 评价指标与方法
地质遗迹评价指标及其权重的合理确定关系着评价结果的准确性。目前，地质遗迹定量评价不外乎从自然属性、功能价值、开发条件等几个方面，但具体评价时选取的评价指标各有不同。评价指标过少，使得结果不够全面，评价指标过多，指标间存在包含或交叉关系，使得评价属性冗余[22]。本文根据天津地质遗迹的特点，通过反复权衡，选取能最大程度体现遗迹价值的因子，去除不相关和冗余因子，得到了如表2所示的评价指标体系。
表 2 天津市地质遗迹评价指标、权重及因子含义Table 2. The evaluation indexes, weights and factor meanings of geological heritage in Tianjing评价综合层 权重 评价项目层 权重 评价因子层 权重 因子含义 地质遗迹价值 0.581 观赏价值 0.154 稀有性 0.062 在国内外出现的几率 奇特性 0.037 形态特征特殊性 完整性 0.027 自然状态、保存状况 愉悦度 0.028 艺术、造型的美观程度 科学价值 0.289 科学研究 0.171 科学研究程度 科普教育 0.118 科普教育程度 文化价值 0.138 历史文化 0.054 历史文化价值 宗教传说 0.084 宗教民俗价值 地质遗迹规模与组合 0.257 遗迹规模 0.130 遗迹面积 0.058 遗迹面积大小 宏伟程度 0.072 遗迹的数量与密度 地域组合 0.127 多样性 0.082 地质遗迹特征点的丰富程度 协调性 0.045 不同遗迹之间的配合程度 地质遗迹外部因素 0.162 资源影响力 0.107 社会认知度 0.056 社会对地质遗迹的认知程度 社会影响力 0.051 地质遗迹对人们的影响程度 环境状况 0.055 环境质量与安全性 0.033 生态环境质量、自然或人为威胁因素 生态旅游环境承载力 0.022 能承受的最大生态旅游活动强度 本文采用层次分析法（AHP）进行评价指标体系中各评价因子权重的确定。层次分析法是通过将评价体系分解成若干层次，通过专家两两对比下一层因子相对于上一层次因子的重要性而逐层判断评分，给出相对重要性定量指标，建立判断矩阵，进而计算特征向量并做一致性检验[23]。运用层次分析法，得出了天津地质遗迹定量评价的各因子权重如表2所示。
4.2.2 评价模型
在确定指标体系中各因子权重后，运用模糊数学法确定综合评分。首先，邀请5位地质遗迹研究相关专家对各处地质遗迹各项指标打分，打分标准采用模糊数学百分制计分法。打分之前，提供每处地质遗迹的详细调查资料，对每处地质遗迹的各要素指标有图文和数据描述。评分标注划分为Ⅰ极好[85,100]、Ⅱ很好[70,85)、Ⅲ较好[55,70)、Ⅳ一般[0,55]等4档。
每位专家依据评价因子打分，选取一个集中度较高的分值来确定每个因子的模糊得分，将各因子得分乘以各自的权重得出该因子的最终得分，把各因子的最终得分相加，计算出地质遗迹的总得分，总分越高，其价值越大，公式如下：
式中：F ——地质遗迹综合评分；
Si——第 i 个评价因子的模糊得分值;
Wi——第 i 个评价因子的权重值；
n——评价因子的数目。
4.2.3 定量评价结果
通过上述方法和模型，对主要地质遗迹的综合分级评价结果值进行等级划分，其依据标准为：世界级（极高的遗迹价值，得分为[85,100]）、国家级（很高的遗迹价值，得分为[70,85)）、省级（较好的遗迹价值，得分为[55,70)）、地方级（一般的遗迹价值，得分[0,55]）。最终结果为世界级4处，国家级3处，省级8处（表3）。
表 3 天津市重要地质遗迹评价等级Table 3. Evaluation grade of important geological heritages in Tianjin评价等级 地质遗迹名称 世界级 蓟县剖面（即中新元古界标准剖面）、锰方硼石矿床、渤海湾贝壳堤、渤海湾牡蛎礁 国家级 七里海湿地、天津市沿海滩涂、天津市地面沉降 省级 盘山花岗岩地貌、蓟州区石英砂岩峰林峡谷地貌、九龙山碳酸盐岩地貌、府君山飞来峰地貌、黄崖关断崖地貌、北大港湿地、团泊洼湿地、大黄堡湿地 5. 结论及建议
（1）天津市地域面积较小，地形地貌、地质背景简单，在漫长的地质历史过程中，共发育了15处重要地质遗迹，按照《地质遗迹调查规范》（DZ/T 0303—2017）分为3大类、7类和10亚类。3大类为：基础地质大类、地貌景观大类和地质灾害大类；7类为：地层剖面类、重要岩矿石产地类、岩土体地貌类、构造地貌类、水体地貌类、海岸地貌类以及地质灾害遗迹类；10亚类为：层型（典型剖面）、典型矿床类露头、侵入岩地貌、碎屑岩地貌、碳酸盐岩地貌、湿地、飞来峰、断层崖、海积地貌、地面沉降。
（2）天津市重要地质遗迹分布范围基本覆盖天津市全域，其中地层剖面类、重要岩矿石产地类、岩土体地貌类、构造地貌类地质遗迹分布在北部基岩出露区，水体地貌类、海岸地貌类以及地质灾害遗迹类分布在南部平原区。
（3）通过定性和定量相结合的方式综合分级评价地质遗迹等级为：世界级4处，蓟州区蓟县剖面、锰方硼石矿床，渤海湾西北岸贝壳堤、牡蛎礁；国家级3处，宁河区七里海湿地、滨海新区沿海滩涂和全市范围内的地面沉降；省级8处，蓟州区盘山花岗岩地貌、石英砂岩峰林峡谷地貌、九龙山碳酸盐岩地貌、府君山飞来峰地貌、黄崖关断崖地貌，滨海新区北大港湿地、静海区团泊洼湿地以及武清区大黄堡湿地。
综上，建议科学合理划定地质遗迹保护边界，适度开展地学旅游和地学科普活动；从多学科多角度、以通俗易懂的方式解读地质遗迹，增强群众保护地质遗迹意识；颁布地质遗迹保护规定，加强地质遗迹巡查监测，及时掌握地质遗迹及保护设施保存现状；结合村镇发展规划，积极申报地质文化村，打造地学科普基地等。
-
表 1 高放废物深地质处置选址及软件使用情况
Table 1 Site selection and simulation software for HLW deep geological disposal
国家 场址 地下实验室 围岩类型 主要用到的模拟软件 中国 新场 北山 花岗岩 Modflow，PMWIN，GoldSim，Phreeqc，EQ3/6，Rockflow，TOUGH系列软件 美国 Yucca Mountain ESF 凝灰岩 Modflow，FEHM，PEST，TOUGH系列软件 瑞典 Forsmark Äspö 花岗岩 Modflow，Darcy Tools，HydroGeoSphere，Connectflow，MIKE 11，MIKE SHE，TOUGH系列软件 法国 Meuse/Haute Marne Bure 黏土岩 Porflow，TOUGH2-MP，Phreeqc，Code-Bright，Comsol Multiphysics 芬兰 Olkiluoto ONKALO 花岗岩 Feflow，Phreeqc，Connectflow，GoldSim，FracMan，HydroGeoSphere 
下载: 导出CSV
-
[1] 杜祥琬. 对我国核能发展战略的几点思考[J]. 中国核电,2012,5(3):194 − 196. [DU Xiangwan. Thought about China’s nuclear power development strategy[J]. China Nuclear Power,2012,5(3):194 − 196. (in Chinese) [2] LIU H H, DOUGHTY C, BODVARSSON G S. An active fracture model for unsaturated flow and transport in fractured rocks[J]. Water Resources Research,1998,34(10):2633 − 2646. DOI: 10.1029/98WR02040
[3] BODVARSSON G S, WU Y S, ZHANG K N. Development of discrete flow paths in unsaturated fractures at Yucca Mountain[J]. Journal of Contaminant Hydrology,2003,62/63:23 − 42. DOI: 10.1016/S0169-7722(02)00177-8
[4] HAUKWA C B, WU Y-S, BODVARSSON G S. Modeling thermal-hydrological response of the unsaturated zone at Yucca Mountain, Nevada, to thermal load at a potential repository[J]. Journal of Contaminant Hydrology,2003,62−63:529 − 552. DOI: 10.1016/S0169-7722(02)00188-2
[5] ZHANG K, WU Y-S, BODVARSSON G S. Parallel computing simulation of fluid flow in the unsaturated zone of Yucca Mountain, Nevada[J]. Journal of Contaminant Hydrology,2003,62-63:381 − 399. DOI: 10.1016/S0169-7722(02)00159-6
[6] YE M, PAN F, WU Y-S, et al. Assessment of radionuclide transport uncertainty in the unsaturated zone of Yucca Mountain[J]. Advances in Water Resources,2007,30(1):118 − 134. DOI: 10.1016/j.advwatres.2006.03.005
[7] RECHARD R P, COTTON T A, VOEGELE M D. Site selection and regulatory basis for the Yucca Mountain disposal system for spent nuclear fuel and high-level radioactive waste[J]. Reliability Engineering & System Safety,2014,122:7 − 31.
[8] RECHARD R P. Results from past performance assessments for the Yucca Mountain disposal system for spent nuclear fuel and high-level radioactive waste[J]. Reliability Engineering & System Safety,2014,122(13):207 − 222.
[9] POHLMANN K, HASSAN A, CHAPMAN J. Description of hydrogeologic heterogeneity and evaluation of radionuclide transport at an underground nuclear test[J]. Journal of Contaminant Hydrology,2000,44(3/4):353 − 386.
[10] ARNOLD B W, KUZIO S P, ROBINSON B A. Radionuclide transport simulation and uncertainty analyses with the saturated-zone site-scale model at Yucca Mountain, Nevada[J]. Journal of Contaminant Hydrology,2003,62/63:401 − 419. DOI: 10.1016/S0169-7722(02)00158-4
[11] ZYVOLOSKI G, KWICKLIS E, EDDEBBARH A A, et al. The site-scale saturated zone flow model for Yucca Mountain: calibration of different conceptual models and their impact on flow paths[J]. Journal of Contaminant Hydrology,2003,62/63(2):731 − 750.
[12] KELKAR S, DING M, CHU S, et al. Modeling solute transport through saturated zone ground water at 10 km scale: Example from the Yucca Mountain license application[J]. Journal of Contaminant Hydrology,2010,117:7 − 25. DOI: 10.1016/j.jconhyd.2010.05.003
[13] KELKAR S, SRINIVASAN G, ROBINSON B A, et al. Breakthrough of contaminant plumes in saturated volcanic rock: implications from the Yucca Mountain site[J]. Geofluids,2013,13(3):273 − 282. DOI: 10.1111/gfl.12035
[14] JOYCE S, HARTLEY L, APPLEGATE D, et al. 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[J]. Hydrogeology Journal,2014,22(6):1233 − 1249. DOI: 10.1007/s10040-014-1165-6
[15] SELROOS J, FOLLIN S. Overview of hydrogeological safety assessment modeling conducted for the proposed high-level nuclear waste repository site at Forsmark, Sweden[J]. Hydrogeology Journal,2014,22(6):1229 − 1232. DOI: 10.1007/s10040-014-1163-8
[16] FOLLIN S, HARTLEY L, RHÉN I, et al. A methodology to constrain the parameters of a hydrogeological discrete fracture network model for sparsely fractured crystalline rock, exemplified by data from the proposed high-level nuclear waste repository site at Forsmark, Sweden[J]. Hydrogeology Journal,2014,22(2):313 − 331. DOI: 10.1007/s10040-013-1080-2
[17] CVETKOVIC V, PAINTER S, OUTTERS N, et al. Stochastic simulation of radionuclide migration in discretely fractured rock near the Äspö Hard Rock Laboratory[J]. Water Resources Research,2004,40(2):1 − 16.
[18] WERNER K, BOSSON E, BERGLUND S. Flow and radionuclide transport from rock to surface systems: characterization and modelling of potential repository sites in Sweden[C]//Proceedings of the 11th international conference on environmental remediation and radioactive waste management, September 2−6, 2007. Bruges: Belgium, 2009: 867 − 872.
[19] SCHWARTZ M O. Modelling radionuclide transport in large fractured-media systems: the example of Forsmark, Sweden[J]. Hydrogeology Journal,2012,20(4):673 − 687. DOI: 10.1007/s10040-012-0837-3
[20] BOSSON E, SELROOS J O, M STIGSSON, et al. Exchange and pathways of deep and shallow groundwater in different climate and permafrost conditions using the Forsmark site, Sweden, as an example catchment[J]. Hydrogeology Journal,2013,21(1):225 − 237. DOI: 10.1007/s10040-012-0906-7
[21] VIDSTRAND P, FOLLIN S, SELROOS J O, et al. Groundwater flow modeling of periods with periglacial and glacial climate conditions for the safety assessment of the proposed high-level nuclear waste repository site at Forsmark, Sweden[J]. Hydrogeology Journal,2014,22(6):1251 − 1267. DOI: 10.1007/s10040-014-1164-7
[22] MONTES-H G, FRITZ B, CLEMENT A, et al. Modeling of transport and reaction in an engineered barrier for radioactive waste confinement[J]. Applied Clay Science,2005,29(3/4):155 − 171.
[23] 李露露, 周志超, 邵景力, 等. 高放废物深地质处置地下水流数值模拟方法研究进展[J]. 水文地质工程地质,2021,48(6):13 − 23. [LI Lulu, ZHOU Zhichao, SHAO Jingli, et al. Advances in groundwater numerical simulation in deep geological disposal of high-level radioactive waste[J]. Hydrogeology & Engineering Geology,2021,48(6):13 − 23. (in Chinese with English abstract) [24] LÖFMAN J. Site scale groundwater flow in Olkiluoto[R]. Helsinki: Posiva Oy, 1999.
[25] KATTILAKOSKI E, KOSKINEN L. Regional-to-site scale groundwater flow in Romuvaara[R]. Helsinki: Posiva Oy, 1999.
[26] LÖFMAN J. Site scale groundwater flow in Olkiluoto- Complementary simulations[R]. Helsinki: Posiva Oy, 2000.
[27] POSIVA O Y. Safety case for the disposal of spent nuclear fuel at Olkiluoto-performance assessment[R]. Eurajoki: Posiva Oy, 2012.
[28] NIEMI A, KOUTIO K, VAITTINEN T, et al. Estimation of block conductivities from hydrologically calibrated fracture networks - description of methodology and application to Romuvaara investigation area[R]. Helsinki: Posiva Oy, 1999.
[29] LANYON G W, MARSCHALL P. Discrete fracture network modelling of a KBS-3H repository at Olkiluoto[R]. Olkiluoto: Posiva Oy, 2006.
[30] POSIVA OY. Olkiluoto site description 2011[R]. Eurajoki: Posiva Oy, 2012.
[31] PULKKANEN V-M, NORDMAN H. Effects of bedrock fractures on radionuclide transport near a vertical deposition hole for spent nuclear fuel[R]. Eurajoki: Posiva Oy, 2011.
[32] POSIVA O Y. Safety case for the disposal of spent nuclear fuel at Olkiluoto - assessment of radionuclide release scenarios for the repository system 2012[R]. Eurajoki: Posiva Oy, 2012.
[33] POTERI A, LAITINEN M. Site-to-canister scale flow and transport in Hästholmen, Kivetty, Olkiluoto and Romuvaara[R]. Helsinki: Posiva Oy, 1999.
[34] BLESSENT D, THERRIEN R, GABLE C W. Large-scale numerical simulation of groundwater flow and solute transport in discretely-fractured crystalline bedrock[J]. Advances in Water Resources,2011,34(12):1539 − 1552. DOI: 10.1016/j.advwatres.2011.09.008
[35] OPHORI D U. A simulation of large-scale groundwater flow and travel time in a fractured rock environment for waste disposal purposes[J]. Hydrological Processes,2004,18(9):1579 − 1593. DOI: 10.1002/hyp.1407
[36] 王驹, 徐国庆. 中国高放废物深地质处置研究[J]. 水文地质工程地质,1998,25(5):7 − 10. [WANG Ju, XU Guoqing. Study on deep geological disposal of high-level radioactive waste in China[J]. Hydrogeology & Engineering Geology,1998,25(5):7 − 10. (in Chinese) [37] 苏锐. 低渗透裂隙介质渗透特征评价技术及其应用研究[D]. 北京: 核工业北京地质研究院, 2007 SU Rui. Studies of hydraulic characterization of low-permeable fractured medium and its application[D]. Beijing: Beijing Institute of Geology for Nuclear Industry, 2007. (in Chinese with English abstract)
[38] 董艳辉, 李国敏. 甘肃北山区域地下水随机模拟研究[C]//废物地下处置学术研讨会. 第三届废物地下处置学术研讨会论文集. 北京: 中国岩石力学与工程学会, 2010: 53-58 DONG Yanhui, LI Guomin. Stochastic simulation of regional groundwater flow in Beishan area[C]//Proceedings of the third symposium on underground waste disposal. Beijing: Chinese Society for Rock Mechanics & Engineering, 2010: 53-58. ( in Chinese with English abstract
[39] 王礼恒. 甘肃北山区域-盆地-岩体多尺度地下水数值模拟研究[D]. 北京: 中国科学院大学, 2015 WANG Liheng. Multi-scale groundwater numerical simulation study of regional-basin-site in Gansu Beishan area[D]. Beijing: University of Chinese Academy of Sciences, 2015. (in Chinese with English abstract)
[40] 魏亚强, 董艳辉, 周鹏鹏, 等. 基于离散裂隙网络模型的核素粒子迁移数值模拟研究[J]. 水文地质工程地质,2017,44(1):123 − 130. [WEI Yaqiang, DONG Yanhui, ZHOU Pengpeng, et al. Numerical simulation of radionuclide particle tracking based on discrete fracture network[J]. Hydrogeology & Engineering Geology,2017,44(1):123 − 130. (in Chinese with English abstract) [41] 尹文杰, 胡立堂, 王景瑞. 基于GRACE重力卫星的甘肃北山地区地下水储量变化规律研究[J]. 水文地质工程地质,2015,42(4):29 − 34. [YIN Wenjie, HU Litang, WANG Jingrui. Changes of groundwater storage variation based on GRACE data at the Beishan area, Gansu Province[J]. Hydrogeology & Engineering Geology,2015,42(4):29 − 34. (in Chinese with English abstract) [42] CAO X, HU L, WANG J, et al. Radionuclide transport model for risk evaluation of high-level radioactive waste in Northwestern China[J]. Human and Ecological Risk Assessment,2017,23(8):2017 − 2032. DOI: 10.1080/10807039.2017.1361811
[43] 曹潇元, 侯德义, 胡立堂. 甘肃北山区域地下水流数值模拟研究[J]. 水文地质工程地质,2020,47(2):9 − 16. [CAO Xiaoyuan, HOU Deyi, HU Litang. Numerical simulation of regional groundwater flow in the Beishan area of Gansu[J]. Hydrogeology & Engineering Geology,2020,47(2):9 − 16. (in Chinese with English abstract) [44] 肖丰, 王驹, 郭永海, 等. 中国高放废物处置库甘肃北山预选区水文地质研究进展[J]. 铀矿地质,2011,27(3):185 − 192. [XIAO Feng, WANG Ju, GUO Yonghai, et al. Progress of hydrogeological studies in Beishan preselected area of disposal repository for high level radioactive waste in China[J]. Uranium Geology,2011,27(3):185 − 192. (in Chinese with English abstract) DOI: 10.3969/j.issn.1000-0658.2011.03.009 [45] 苏锐, 王驹, 陈伟明, 等. 放射性核素在CRP-GEORC地质处置库远场中迁移的灵敏性与不确定性分析[C]//废物地下处置学术研讨会. 第二届废物地下处置学术研讨会论文集. 北京: 中国岩石力学与工程学会, 2008: 327 − 336 SU Rui, WANG Ju, CHEN Weiming, et al. Sensitivity analysis and uncertainty simulation of the migration of radionuclide in the system of geological disposal-CRP-GEORC model[C]//Symposium on underground waste disposal. proceedings of the second conference on underground waste disposal. Beijing: Chinese Society for Rock Mechanics & Engineering, 2008: 327 − 336. (in Chinese with English abstract)
[46] 吴晓东. 高放废物在地质处置中核素迁移数值模拟研究[D]. 抚州: 东华理工大学, 2012 WU Xiaodong. Nuclide transport numerical simulation of high level radioactive waste geological disposal[D]. Fuzhou: East China Institute of Technology, 2012. (in Chinese with English abstract)
[47] 凌辉, 王驹, 陈伟明. 花岗岩天然屏障对放射性核素的阻滞性能模拟研究[J]. 中国科技论文,2017,12(21):2507 − 2511. [LING Hui, WANG Ju, CHEN Weiming. Simulation study on retardation of nuclides by granite geological barriers[J]. China Sciencepaper,2017,12(21):2507 − 2511. [48] 凌辉, 王驹, 陈伟明. 高放废物地质处置算井子候选场址核素迁移模拟研究[J]. 铀矿地质,2018,34(2):118 − 123. [LING Hui, WANG Ju, CHEN Weiming. Preliminary evaluation of radioactive effect of candidate sites in Suanjingzi sub-area for geological disposal of high-level radioactive waste[J]. Uranium Geology,2018,34(2):118 − 123. [49] 彭志娟. 高放废物处置库EDZ中核素129I、135Cs的迁移模拟研究[D]. 抚州: 东华理工大学, 2019 PENG Zhijuan. Numerical simulation of 129I & 135Cs nuclide migration in EDZ of the high-level radioactive waste repository[D]. Fuzhou: East China Institute of Technology, 2019. (in Chinese with English abstract)
[50] 李楚. 镅在北山预选处置场研究区的赋存形态及迁移行为探究[D]. 抚州: 东华理工大学, 2019 LI Chu. Study on the occurrence and migration behavior of Am in the research area of Beishan pre-selected disposal site[D]. Fuzhou: East China Institute of Technology, 2019. (in Chinese with English abstract)
-
期刊类型引用(8)
1. 王凯. 物探技术下废物处置场地下水污染地质勘查方法研究. 环境科学与管理. 2025(06): 139-144 . 
百度学术
2. 刘国东，杜成鸿，侯杰，杨梦溪，陈宇，谢杨. 基于用高斯分布生成管道的岩溶地下水流数值模拟. 吉林大学学报(地球科学版). 2024(02): 592-603 . 百度学术
3. 王瑞青，陈超，谢添，陈文杰，朱君. ~(90)Sr在微风化二长花岗岩中迁移实验及模拟研究. 核化学与放射化学. 2024(03): 258-263 . 百度学术
4. 刘栋永. 人工智能在水文地质调查过程中的战略意义探索. 世界有色金属. 2024(10): 166-168 . 百度学术
5. 王志秀，李亚松，郝奇琛，张媛静，肖勇，刘春雷. 基于盐度动态模拟估算潮间带地下淡水排泄量. 水文地质工程地质. 2024(05): 56-67 . 本站查看
6. 闫成宏. 某预选低中放废物处置场岩土层渗透性能综合研究. 城市建设理论研究(电子版). 2024(34): 79-81 . 百度学术
7. 王顿，裴丽欣，张礼中，樊连杰，卢丽，李习文，刘潇桐，李俊楠，梁林德，白雪冬. 基于文献计量学近20年国内外地下水数值模拟研究进展及展望. 环境工程. 2023(S1): 240-247 . 百度学术
8. 李婕，赵键，刘海峰，朱永昌，王华，张行泉，郑奎，霍冀川. 锆基钙钛矿结构氧化物对Cs的玻璃陶瓷固化. 非金属矿. 2023(04): 88-91 . 百度学术
其他类型引用(9)
计量
- 文章访问数: 598
- HTML全文浏览量: 253
- PDF下载量: 334
- 被引次数: 17
邮件订阅
RSS