CN114707337B - 一种放射性废物处置过程中的核素迁移模拟方法及系统 - Google Patents
一种放射性废物处置过程中的核素迁移模拟方法及系统 Download PDFInfo
- Publication number
- CN114707337B CN114707337B CN202210377800.XA CN202210377800A CN114707337B CN 114707337 B CN114707337 B CN 114707337B CN 202210377800 A CN202210377800 A CN 202210377800A CN 114707337 B CN114707337 B CN 114707337B
- Authority
- CN
- China
- Prior art keywords
- nuclide
- migration
- migration process
- component
- waste disposal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 230000005012 migration Effects 0.000 title claims abstract description 419
- 238000013508 migration Methods 0.000 title claims abstract description 419
- 238000000034 method Methods 0.000 title claims abstract description 358
- 230000008569 process Effects 0.000 title claims abstract description 292
- 239000002901 radioactive waste Substances 0.000 title claims abstract description 111
- 238000004088 simulation Methods 0.000 title claims abstract description 79
- 230000004888 barrier function Effects 0.000 claims abstract description 123
- 238000013178 mathematical model Methods 0.000 claims abstract description 35
- 238000011156 evaluation Methods 0.000 claims abstract description 30
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 109
- 239000011148 porous material Substances 0.000 claims description 89
- 230000004907 flux Effects 0.000 claims description 57
- 239000011435 rock Substances 0.000 claims description 41
- 238000009792 diffusion process Methods 0.000 claims description 39
- 239000006185 dispersion Substances 0.000 claims description 24
- 230000008859 change Effects 0.000 claims description 23
- 230000009969 flowable effect Effects 0.000 claims description 21
- 230000005258 radioactive decay Effects 0.000 claims description 13
- 229920006395 saturated elastomer Polymers 0.000 claims description 8
- 230000004941 influx Effects 0.000 claims description 7
- 230000007774 longterm Effects 0.000 abstract description 6
- 239000011159 matrix material Substances 0.000 description 21
- 238000001179 sorption measurement Methods 0.000 description 16
- 241000894007 species Species 0.000 description 12
- 239000011521 glass Substances 0.000 description 10
- 239000002927 high level radioactive waste Substances 0.000 description 9
- 239000000463 material Substances 0.000 description 9
- 210000004027 cell Anatomy 0.000 description 7
- 238000010586 diagram Methods 0.000 description 7
- 239000010438 granite Substances 0.000 description 7
- 238000009375 geological disposal Methods 0.000 description 6
- 239000003673 groundwater Substances 0.000 description 6
- 229910000975 Carbon steel Inorganic materials 0.000 description 5
- 239000010962 carbon steel Substances 0.000 description 5
- 238000004090 dissolution Methods 0.000 description 5
- 238000009412 basement excavation Methods 0.000 description 4
- 239000000440 bentonite Substances 0.000 description 4
- 229910000278 bentonite Inorganic materials 0.000 description 4
- SVPXDRXYRYOSEX-UHFFFAOYSA-N bentoquatam Chemical compound O.O=[Si]=O.O=[Al]O[Al]=O SVPXDRXYRYOSEX-UHFFFAOYSA-N 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 239000000126 substance Substances 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 238000010276 construction Methods 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- 230000012010 growth Effects 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 3
- 239000007787 solid Substances 0.000 description 3
- 238000012546 transfer Methods 0.000 description 3
- 239000002699 waste material Substances 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 230000007797 corrosion Effects 0.000 description 2
- 238000005260 corrosion Methods 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000012854 evaluation process Methods 0.000 description 2
- 239000000945 filler Substances 0.000 description 2
- -1 force Substances 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000017111 nuclear migration Effects 0.000 description 2
- 210000003429 pore cell Anatomy 0.000 description 2
- 238000001556 precipitation Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000002285 radioactive effect Effects 0.000 description 2
- 241000282412 Homo Species 0.000 description 1
- 241000282414 Homo sapiens Species 0.000 description 1
- 102100021010 Nucleolin Human genes 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 229910052768 actinide Inorganic materials 0.000 description 1
- 150000001255 actinides Chemical class 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000004913 activation Effects 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000031018 biological processes and functions Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000003139 buffering effect Effects 0.000 description 1
- 230000012292 cell migration Effects 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 238000012412 chemical coupling Methods 0.000 description 1
- 239000000084 colloidal system Substances 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000013210 evaluation model Methods 0.000 description 1
- 230000004992 fission Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 230000033001 locomotion Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 108010044762 nucleolin Proteins 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 230000000191 radiation effect Effects 0.000 description 1
- 238000004055 radioactive waste management Methods 0.000 description 1
- 230000000979 retarding effect Effects 0.000 description 1
- 238000007711 solidification Methods 0.000 description 1
- 230000008023 solidification Effects 0.000 description 1
- 239000002915 spent fuel radioactive waste Substances 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 238000011144 upstream manufacturing Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Tourism & Hospitality (AREA)
- General Physics & Mathematics (AREA)
- Economics (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- General Health & Medical Sciences (AREA)
- General Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Educational Administration (AREA)
- Development Economics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Processing Of Solid Wastes (AREA)
Abstract
本发明涉及一种放射性废物处置过程中的核素迁移模拟方法及系统。所述方法包括：对放射性废物处置设施多重屏障系统中核素迁移过程进行描述，确定核素在多重屏障系统不同单元中的迁移过程；结合万年以上的安全评价要求，以及核素在多重屏障系统不同单元中迁移过程的差异性，对核素迁移过程进行概化；在概化后的核素迁移过程基础上，建立核素迁移过程的数学模型；采用核素迁移过程的数学模型对核素在所述放射性废物处置设施多重屏障系统不同单元中的迁移过程进行模拟。采用本发明方法能够在万年尺度上实现对核素在放射性废物处置设施多重屏障系统中迁移的模拟，为放射性废物处置设施长期安全评价提供可行的解决方案。
Description
技术领域
本发明涉及放射性废物处置技术领域，特别是涉及一种放射性废物处置过程中的核素迁移模拟方法及系统。
背景技术
放射性废物安全处置是放射性废物管理的关键一环，是确保核能可持续发展的重要保障。在复杂地下环境长期作用下，放射性核素可能会出现脱稳，从放射性废物处置设施中释放并迁移到地质环境中，造成局部地下水的污染，甚至随地下水迁移到达生物圈引起生态环境安全风险。因此，必须通过安全评价，对放射性核素的释放和迁移过程进行模拟预测，以评价放射性废物处置设施的长期安全性，为其选址、设计优化和安全防护提供支撑。
相较于常规核电站和污染治理工程数十年的评价周期，放射性废物处置设施核素迁移模拟预测的时间尺度在万年以上。万年尺度上的安全评价具有极大的不确定性，一般无法对所有的核素迁移过程进行精确的模拟，需要对多重屏障系统中的核素迁移过程进行概化和近似，或以保守估计的方式实现对核心迁移过程的模拟预测。现有的以(精确)现实模拟为核心的核素迁移模拟解决方案，无法直接应用到放射性废物处置设施安全评价中。与此同时，已有核素迁移模拟方法往往只适用于单一的固体介质(孔隙介质或裂隙介质)，而放射性废物处置多重屏障系统包含源项、工程屏障、天然屏障和生物圈等子系统，其核素迁移路径上可能同时包含孔隙介质和裂隙介质。因此，有必要建立一种适应万年尺度安全评价需要的核素迁移模拟方法，为放射性废物处置安全评价提供可行的解决方案。
发明内容
本发明的目的是提供一种放射性废物处置过程中的核素迁移模拟方法及系统，以在万年尺度上实现对核素在放射性废物处置设施多重屏障系统中迁移的全过程可靠模拟。
为实现上述目的，本发明提供了如下方案：
一种放射性废物处置过程中的核素迁移模拟方法，包括：
对放射性废物处置设施多重屏障系统中核素迁移过程进行描述，确定核素在所述放射性废物处置设施多重屏障系统不同单元中的迁移过程；所述放射性废物处置设施多重屏障系统不同单元包括源项、工程屏障、天然屏障和生物圈；
结合万年以上的安全评价要求，根据所述核素在放射性废物处置设施多重屏障系统不同单元中的迁移过程的差异性，对核素迁移过程进行概化；
在所述概化后核素迁移过程的基础上，建立核素迁移过程的数学模型；所述核素迁移过程的数学模型包括孔隙单元迁移模型、含水通道迁移模型和单裂隙迁移模型；
采用所述核素迁移过程的数学模型对所述核素在所述放射性废物处置设施多重屏障系统不同单元中的迁移过程进行模拟。
可选地，所述对放射性废物处置设施多重屏障系统中核素迁移过程进行描述，确定核素在所述放射性废物处置设施多重屏障系统不同单元中的迁移过程，具体包括：
对放射性废物处置设施多重屏障系统中核素迁移过程进行描述，确定核素在所述放射性废物处置设施多重屏障系统不同单元中的迁移过程；所述迁移过程包括放射性衰变的物理过程以及核素在所述放射性废物处置设施多重屏障系统不同单元中迁移的近场过程和远场过程。
可选地，所述结合万年以上的安全评价要求，根据所述核素在放射性废物处置设施多重屏障系统不同单元中的迁移过程的差异性，对核素迁移过程进行概化，具体包括：
结合万年以上的安全评价要求，根据所述核素在放射性废物处置设施多重屏障系统不同单元中的迁移过程的差异性，将核素在孔隙介质中的迁移过程概化为核素在均匀连续的孔隙单元间的迁移过程，将核素在弱导水性裂隙介质中以对流-扩散形式为主的迁移过程概化为核素随地下水在一维连续通道中的迁移过程；将核素在导水性裂隙介质中以对流-弥散形式为主的迁移过程概化为核素随地下水在一维裂隙通道中的迁移过程。
可选地，所述在所述概化后核素迁移过程的基础上，建立核素迁移过程的数学模型，具体包括：
在所述概化后核素迁移过程的基础上，建立所述孔隙单元迁移模型
其中M′mi为网格m中核素i的质量增加率；Mmi为网格m中核素i的质量；Dri为核素i的衰变速率；NJi表示核素i的直系母核素的数量；Mmj为网格m中母核素j的质量；Drj为母核素j的衰变速率；Pj-i表示衰变为i的母核素j的比例；Sj-i表示每摩尔j衰变产生i的摩尔化学计量比；Wi表示核素i的摩尔质量；Wj表示母核素j的摩尔质量；NEm为与单元m相连接的质量通量连接数量；fci为通过质量通量连接l的核素i的流入通量；Wmi为外部源汇项对网格m中核素i的直接补给量；
在所述概化后核素迁移过程的基础上，建立所述含水通道迁移模型
其中m′is为单元i中组分s的质量增加率；mis为单元i中组分s的质量；Drs为组分s的衰变变化速率；NPs表示组分s的直系父组分p的数目；mip为单元i中父组分p的质量；Drp为父组分p的衰变变化速率；fp-s表示衰变为组分s的父组分p的百分比；Sp-s表示每摩尔组分p衰变产生组分s的摩尔化学计量比；Ws表示组分s的摩尔质量；Wp表示父组分p的摩尔质量；NEi为与单元i相连接的质量通量连接c的数量；fcs为通过质量通量连接c的组分s的流入通量；
在所述概化后核素迁移过程的基础上，建立所述单裂隙迁移模型
其中Fk表示组分k在裂隙通道出口端的通量；q表示裂隙通道中的体积流速；ck表示裂隙通道中可流动区域内单位有效饱和孔隙空间中的组分k在活动区多孔介质的有效孔隙中的平均浓度；Sp表示可流动区域的横截面积；Dk表示组分k可流动区内的有效扩散系数；α表示通道的弥散度；L表示通道的长度；x表示通道内的位置。
可选地，所述采用所述核素迁移过程的数学模型对所述核素在所述放射性废物处置设施多重屏障系统不同单元中的迁移过程进行模拟，具体包括：
采用所述孔隙单元迁移模型模拟核素在孔隙介质中的迁移过程；所述孔隙介质包括源项、工程屏障、天然屏障、生物圈；
采用所述含水通道迁移模型模拟核素在弱导水性裂隙介质中的迁移过程；所述弱导水性裂隙介质包括导水性极低的围岩；
采用所述单裂隙迁移模型模拟核素在导水性裂隙介质中的迁移过程；所述导水性裂隙介质包括具有一定导水性的围岩裂隙。
一种放射性废物处置过程中的核素迁移模拟系统，包括：
核素迁移过程描述模块，用于对放射性废物处置设施多重屏障系统中核素迁移过程进行描述，确定核素在所述放射性废物处置设施多重屏障系统不同单元中的迁移过程；所述放射性废物处置设施多重屏障系统不同单元包括源项、工程屏障、天然屏障和生物圈；
核素迁移过程概化模块，用于结合万年以上的安全评价要求，根据所述核素在放射性废物处置设施多重屏障系统不同单元中的迁移过程的差异性，对核素迁移过程进行概化；
核素迁移过程数学模型建立模块，用于在所述概化后核素迁移过程的基础上，建立核素迁移过程的数学模型；所述核素迁移过程的数学模型包括孔隙单元迁移模型、含水通道迁移模型和单裂隙迁移模型；
核素迁移过程模拟模块，用于采用所述核素迁移过程的数学模型对所述核素在所述放射性废物处置设施多重屏障系统不同单元中的迁移过程进行模拟。
可选地，所述核素迁移过程描述模块具体包括：
核素迁移过程描述单元，用于对放射性废物处置设施多重屏障系统中核素迁移过程进行描述，确定核素在所述放射性废物处置设施多重屏障系统不同单元中的迁移过程；所述迁移过程包括放射性衰变的物理过程以及核素在所述放射性废物处置设施多重屏障系统不同单元中迁移的近场过程和远场过程。
可选地，所述核素迁移过程概化模块具体包括：
核素迁移过程概化单元，用于结合万年以上的安全评价要求，根据所述核素在放射性废物处置设施多重屏障系统不同单元中的迁移过程的差异性，将核素在孔隙介质中的迁移过程概化为核素在均匀连续的孔隙单元间的迁移过程，将核素在弱导水性裂隙介质中以对流-扩散形式为主的迁移过程概化为核素随地下水在一维连续通道中的迁移过程；将核素在导水性裂隙介质中以对流-弥散形式为主的迁移过程概化为核素随地下水在一维裂隙通道中的迁移过程。
可选地，所述核素迁移过程数学模型建立模块具体包括：
孔隙单元迁移模型建立单元，用于在所述概化后核素迁移过程的基础上，建立所述孔隙单元迁移模型
其中M′mi为网格m中核素i的质量增加率；Mmi为网格m中核素i的质量；Dri为核素i的衰变速率；NJi表示核素i的直系母核素的数量；Mmj为网格m中母核素j的质量；Drj为母核素j的衰变速率；Pj-i表示衰变为i的母核素j的比例；Sj-i表示每摩尔j衰变产生i的摩尔化学计量比；Wi表示核素i的摩尔质量；Wj表示母核素j的摩尔质量；NEm为与单元m相连接的质量通量连接数量；fci为通过质量通量连接l的核素i的流入通量；Wmi为外部源汇项对网格m中核素i的直接补给量；
含水通道迁移模型建立单元，用于在所述概化后核素迁移过程的基础上，建立所述含水通道迁移模型
其中m′is为单元i中组分s的质量增加率；mis为单元i中组分s的质量；Drs为组分s的衰变变化速率；NPs表示组分s的直系父组分p的数目；mip为单元i中父组分p的质量；Drp为父组分p的衰变变化速率；fp-s表示衰变为组分s的父组分p的百分比；Sp-s表示每摩尔组分p衰变产生组分s的摩尔化学计量比；Ws表示组分s的摩尔质量；Wp表示父组分p的摩尔质量；NEi为与单元i相连接的质量通量连接c的数量；fcs为通过质量通量连接c的组分s的流入通量；
单裂隙迁移模型建立单元，用于在所述概化后核素迁移过程的基础上，建立所述单裂隙迁移模型
其中Fk表示组分k在裂隙通道出口端的通量；q表示裂隙通道中的体积流速；ck表示裂隙通道中可流动区域内单位有效饱和孔隙空间中的组分k在活动区多孔介质的有效孔隙中的平均浓度；Sp表示可流动区域的横截面积；Dk表示组分k可流动区内的有效扩散系数；α表示通道的弥散度；L表示通道的长度；x表示通道内的位置。
可选地，所述核素迁移过程模拟模块具体包括：
孔隙介质迁移过程模拟单元，用于采用所述孔隙单元迁移模型模拟核素在孔隙介质中的迁移过程；所述孔隙介质包括源项、工程屏障、天然屏障、生物圈；
弱导水性裂隙介质迁移过程模拟单元，用于采用所述含水通道迁移模型模拟核素在弱导水性裂隙介质中的迁移过程；所述弱导水性裂隙介质包括导水性极低的围岩；
导水性裂隙介质迁移过程模拟单元，用于采用所述单裂隙迁移模型模拟核素在导水性裂隙介质中的迁移过程；所述导水性裂隙介质包括具有一定导水性的围岩裂隙。
根据本发明提供的具体实施例，本发明公开了以下技术效果：
本发明提供了一种放射性废物处置过程中的核素迁移模拟方法及系统，所述方法包括：对放射性废物处置设施多重屏障系统中核素迁移过程进行描述，确定核素在所述放射性废物处置设施多重屏障系统不同单元中的迁移过程；所述放射性废物处置设施多重屏障系统不同单元包括源项、工程屏障、天然屏障和生物圈；结合万年以上的安全评价要求，根据所述核素在放射性废物处置设施多重屏障系统不同单元中的迁移过程的差异性，对核素迁移过程进行概化；在所述概化后核素迁移过程的基础上，建立核素迁移过程的数学模型；所述核素迁移过程的数学模型包括孔隙单元迁移模型、含水通道迁移模型和单裂隙迁移模型；采用所述核素迁移过程的数学模型对所述核素在所述放射性废物处置设施多重屏障系统不同单元中的迁移过程进行模拟。采用本发明方法可以有效克服对多重屏障系统中核素迁移过程进行可靠性模拟的技术难题，能够在万年尺度上实现对核素在放射性废物处置设施多重屏障系统中迁移的全过程可靠模拟，为放射性废物处置设施长期安全评价提供可行的解决方案。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明一种放射性废物处置过程中的核素迁移模拟方法的流程图；
图2为本发明实施例提供的高放废物深地质处置概念示意图；
图3为本发明实施例提供的高放废物处置设施核素迁移涉及的主要过程的示意图；
图4为本发明实施例提供的地下水运移与核素迁移路径的示意图；
图5为本发明一种放射性废物处置过程中的核素迁移模拟系统的结构图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明的目的是提供一种放射性废物处置过程中的核素迁移模拟方法及系统，以在万年尺度上实现对核素在放射性废物处置设施多重屏障系统中迁移的全过程可靠模拟。
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1为本发明一种放射性废物处置过程中的核素迁移模拟方法的流程图。参见图1，本发明一种放射性废物处置过程中的核素迁移模拟方法包括：
步骤101：对放射性废物处置设施多重屏障系统中核素迁移过程进行描述，确定核素在所述放射性废物处置设施多重屏障系统不同单元中的迁移过程。
放射性废物处置设施多重屏障系统(以下简称多重屏障系统)由内而外依次是源项、工程屏障、天然屏障(天然围岩)、生物圈等。对核素迁移过程进行概化前，需要对与安全相关的核素迁移过程(包括核素衰变、物理、化学、生物等过程)进行全面描述。
图2为本发明实施例提供的高放废物深地质处置概念示意图。参见图2，下面以高放废物深地质处置设施为例，并考虑以花岗岩场址为参考场址进行高放废物处置设施多重屏障系统核素迁移模型构建和核素迁移模拟。
首先从处置概念出发，对多重屏障系统中核素迁移过程进行探讨，为过程概化和模型建立提供参考。
高放废物地质处置的长期安全建立在多重屏障系统基础之上。现阶段考虑位于地下500-600米饱水层花岗岩围岩中的矿山式处置设施，其中源项为玻璃固化体，由内而外依次被废物处置容器(碳钢)、缓冲回填材料(膨润土)、围岩(花岗岩)所包覆。选址过程保证了处置设施处置坑附近围岩的适宜性，处置设施及其处置坑安全距离外的围岩中存在导水节理带、导水断裂带。受开挖的影响，缓冲回填材料与完整围岩之间存在一个开挖扰动带。
玻璃固化体中除了玻璃外，主要成分是经过后处理的轻水反应堆乏燃料，包含锕系元素、裂变产物和活化产物，其具体成分由反应堆设计、冷却时间、后处理技术、处置安全要求等因素共同决定，大致包含四十余种核素。总体上，可以按照时间顺序，以核素随地下水迁移路径为主线，把与安全评价相关的物理、化学、生物过程(放射性衰变等为贯穿始终的过程)，大致归类为近场过程和远场过程进行梳理。
图3为本发明实施例提供的高放废物处置设施核素迁移涉及的主要过程的示意图。参见图3，核素在放射性废物处置设施多重屏障系统不同单元中的迁移过程包括放射性衰变的物理过程以及核素在放射性废物处置设施多重屏障系统不同单元中迁移的近场过程和远场过程。
(1)放射性衰变
放射性衰变是安全评价全过程都要考虑的物理过程，其影响有二，其一是废物体释热功率，其二是各核素的浓度和盘存量。
(2)近场过程
当前处置概念下，近场是指工程屏障(玻璃固化体、碳钢处置容器、膨润土缓冲回填材料)和近场花岗岩围岩(开挖扰动带及其受热影响的围岩)。近场过程主要涉及如下事件链：地下水从近场围岩流过缓冲材料，与碳钢接触；地下水腐蚀碳钢处置容器，并最终导致处置容器失效；地下水流入处置容器，与玻璃固化体接触；放射性核素随玻璃固化体的溶解而释出；释出量大于溶解度限制的核素形成二次沉淀；溶解的核素通过碳钢处置容器进入膨润土缓冲回填材料，腐蚀产物可能对核素有一定的吸附作用；溶解的核素通过缓冲回填材料进入近场围岩，缓冲回填材料对核素产生吸附并过滤掉玻璃固化体溶解产生的胶体；溶解的核素随地下水进入远场围岩。
处置设施区域范围内的水文地质条件属于场址调查的范畴。处置设施在建设、运行、关闭后阶段都会影响地下水流特性，进而影响近场中的核素迁移，主要体现在三个方面：第一，处置巷道、硐室会改变围岩中的导水通道；第二、处置设施整体释热会以热水耦合的方式影响处置库附近地下水流场；第三，单个处置坑附近由于膨润土吸水膨胀等原因，以热、水、力、化耦合的方式影响局部的地下水流动。
(3)远场过程
远场是指处置设施产生的热、水、力、化等影响小到可以忽略不计的岩石圈区域。核素随地下水由近场进入远场发生数千年以后，因而远场过程相较近场过程存在更大的不确定性，一定程度上需要保守估计。
图4为本发明实施例提供的地下水运移与核素迁移路径的示意图。参考图4中的核素迁移路径，可以把远场过程依次分为：核素在完整围岩、导水节理带、导水断裂带中的迁移这三个阶段，其中导水节理带和导水断裂带主要是依据裂隙密集度和裂隙尺度的不同进行划分。其中第一个阶段核素运移以扩散为主。而后两个阶段核素运移的主导机制则是花岗岩不同结构面(裂隙网络)中的对流弥散过程。裂隙表面和裂隙填充物对核素的吸附一定程度上会阻滞核素的迁移。核素会沿垂直于裂隙的方向向岩体内扩散，这种基质扩散也会对核素迁移有阻滞作用。低孔隙度花岗岩岩体内的孔隙水占总含水量的很大部分，因此在裂隙水流速较慢的条件下，基质扩散的阻滞作用十分重要。流动的裂隙水与静止的孔隙水在裂隙网络中的相互作用可以用双孔隙模型来描述。完整花岗岩中孔隙的连通性是一个十分重要的问题，只有与裂隙连通的孔隙才能发生基质扩散作用。
核素通过远场释放到达生物圈，从而可能对人类造成辐射影响。对处置系统安全的评价主要是预测放射性核素释放率和个人有效剂量率，因此安全评价模型就必须模拟并预测核素从处置设施到生物圈的迁移。核素迁移模拟的过程有玻璃固化体溶解，核素释放，核素溶解/沉淀，核素放射性衰变/内增，核素在废物处置容器腐蚀产物和缓冲层中的扩散和吸附，核素在近场和远场岩石中的平流、扩散/弥散、吸附、基质扩散和吸附等，如表1所示。
表1高放废物地质处置系统组成部分和作用过程
步骤102：结合万年以上的安全评价要求，根据所述核素在放射性废物处置设施多重屏障系统不同单元中的迁移过程的差异性，对核素迁移过程进行概化。
放射性废物处置设施多重屏障系统包含源项、工程屏障、天然屏障和生物圈等子系统，核素随地下水在多重屏障系统中的迁移，涉及物理、化学、生物等多场多尺度的耦合过程，迁移机理复杂，且放射性衰变贯穿始终。放射性废物处置设施安全评价开展过程中，一般无法对所有的核素迁移过程进行精确的模拟，这就需要在评价的过程中对多重屏障系统中的核素迁移过程进行概化和近似，或以保守估计的方式实现对核心迁移过程的模拟预测，并以此为基础对核素迁移通道进行模块化的拆分，形成数学上、数值方法和算法上均可实现的计算功能要求，以便实现对多重屏障系统全过程中核素迁移过程的可靠模拟。
核素迁移模型的代表性及其可靠性需要满足放射性废物处置安全评价时间尺度要求。比如，对于高放废物地质处置，核素迁移模拟的时间一般在万年以上。因此，需要对万年时间尺度上的核素迁移过程进行保守估计(概化)。
与安全评价紧密相关的过程，时间空间尺度跨度大、内容丰富、综合性极强。核素迁移过程的概化、近似、模块拆分、保守估计等，是对放射性废物处置设施多重屏障系统的物理、化学、生物等过程的简化描述，是以相对简单的方式确保安全评价得以有效开展的关键步骤。
放射性废物处置设施多重屏障系统由内而外依次是源项、工程屏障、天然围岩、生物圈等。源项、工程屏障和生物圈一般为孔隙介质，天然围岩一般为孔隙介质或裂隙介质(其中，裂隙介质又细分为弱导水性裂隙和导水性裂隙)。核素随地下水由近场进入远场，从远场进入生物圈一般发生在数千年以后，因而远场过程和生物圈相较近场过程存在更大的不确定性，一定程度上需要保守估计，为各个子系统相应的迁移模型的模块化拆分提供依据。
在全面描述多重屏障系统核素迁移过程的基础上，对万年时间尺度上的核素迁移过程进行概化。多重屏障系统由内而外依次是源项、工程屏障、天然围岩、生物圈等。源项、工程屏障和生物圈一般为孔隙介质，天然围岩一般为孔隙介质或裂隙介质(其中，裂隙介质又细分为导水性强的裂隙和导水性弱的裂隙)。核素在孔隙介质中的迁移过程，可概化为核素在均匀连续的孔隙单元间的迁移过程；核素在弱导水性裂隙介质中以对流-扩散形式为主的迁移过程，可概化为核素随地下水在一维连续通道中的迁移过程；核素在导水性裂隙介质中以对流-弥散形式为主的迁移过程，可概化为核素随地下水在一维裂隙通道中的迁移过程。
步骤103：在所述概化后核素迁移过程的基础上，建立核素迁移过程的数学模型。
此步骤涉及核素在包含源项、工程屏障、天然屏障、生物圈等的多重屏障系统中迁移模型的构建，可为放射性废物处置设施的安全评价提供计算支撑。在概化后核素迁移过程的基础上，对核素迁移模型进行模块化划分。根据孔隙介质和不同导水性裂隙介质中核素迁移模拟过程的差异性，对核素迁移模型进行模块化划分，用孔隙单元迁移模型、含水通道迁移模型、单裂隙迁移模型分别模拟核素在孔隙介质、导水性弱的裂隙介质、导水性强的裂隙介质中的迁移过程。
多重屏障系统核素迁移路径上伴随着核素衰变、吸附、核素迁移等，但各个屏障单元中核素迁移的形式存在一定的差异性。对于源项单元，侧重对核素溶解的模拟；对于工程屏障单元，侧重对极低水力条件下核素扩散的模拟；对于低水力条件下的天然屏障，侧重对对流-扩散的模拟；对于极低水力条件下的天然屏障，侧重对对流、弥散、围岩基质扩散的模拟；对于生物圈单元，侧重对核素迁移及其剂量吸附的模拟。因此，需要针对核素迁移路径上不同的迁移形式，进行核素迁移过程数学模型的设计，以便构建适合不同单元的核素迁移模型。
针对放射性核素随地下水在介质中混合传输过程的模拟需要，构建孔隙单元迁移模型以实现对核素衰变、扩散、吸附等过程的模拟，应用对象包括源项、工程屏障、天然屏障、生物圈等；针对核素随地下水在一维连续通道的迁移，构建含水通道迁移模型以实现对核素衰变、对流、弥散等过程的模拟，应用对象主要是导水性极低的围岩；针对核素随地下水在一维通道的迁移，构建单裂隙迁移模型以实现对核素衰变、对流、弥散、基质扩散、吸附等过程的模拟，应用对象主要是具有一定导水性的围岩裂隙。实际使用中，通过对各类型核素迁移模型的组合表达，能够对现实环境进行综合模拟，可实现对核素在源项、工程屏障、天然屏障和生物圈中的迁移过程的模拟。
上述孔隙单元迁移模型、含水通道迁移模型、单裂隙迁移模型，模拟的对象不同，侧重的核素迁移形式一般也不同，其对应的主控方程往往差异明显，接下来针对上述模型，构建与之对应的数学表达模型及其求解方法。
1、孔隙单元迁移模型
单元网格m中核素i的质量变化遵从如下的质量守恒控制方程：
其中，M′mi为网格m中核素i的质量增加率[g/s]；Mmi为网格m中核素i的质量[g]；Dri为核素i的衰变速率[1/s]；NJi表示核素i的直系母核素的数量；Mmj为网格m中母核素j的质量[g]；Drj为母核素j的衰变速率[1/s]；Pj-i表示衰变为i的母核素j的比例；Sj-i表示每摩尔j衰变产生i的摩尔化学计量比；Wi表示核素i的摩尔质量(数值上等于相对分子质量)[g/mol]；Wj表示母核素j的摩尔质量[g/mol]；NEm为与单元m相连接的质量通量连接l的数量；fci为通过质量通量连接l的核素i的流入通量[g/s]；Wmi为外部源汇项对网格m中核素i的直接补给量[g/s]。
由质量平衡控制方程(1)可以看出，网格m中核素的质量变化主要由自身衰变项、内生长项、质量通量项以及源汇项四个部分组成。对于孔隙单元连接各核素所构成的方程组，可采用向后的时间隐式差分对单元网格进行求解：
M(t+Δt)＝[I-ΔtDI-ΔtFL]-1M(t) (2)
M(t+Δt)为包含衰变族中所有组分的质量变化量矩阵(g/s)；I表示单元网络排出的质量在其离开的那一刻就不再衰减；Δt为时间变化率；M(t)为包含衰变族中所有组分的质量矩阵(g)；DI为包含衰变组中衰变与增长率的矩阵(1/a)；FL为对流引起的通过质量通量连接的质量变化率矩阵(1/a)。
每个网络求解的方程式是根据节点定义的，其中一个节点代表单元网格中的一个单元(主节点)和该单元网格向其他模块排出物质的临时监测点(辅助节点)。为需要计算通量的所有质量通量链接定义了一个辅助节点。增加观测点M2，观测该孔隙单元的质量输出量。以下为是GoldSim为每个衰变链求解的基本方程式
当质量进入一个下游模块中时，衰减过程恢复。假定从上游模块排出的组分在应用于下游通道时在时间步长上均匀分布。
公式(3)的求解方式如下：
求解方法采用主元高斯-约当(G-J)消去法作为参考。
2、含水通道迁移模型
含水通道迁移模型实质上是孔隙单元网格以适当的方式链接而成的含水层连续流通道，构成含水层连续流通道的某一单元网格的控制方程如公式(5)所示：
其中，m′is为单元i中组分s的质量增加率(g/s)；mis为单元i中组分s的质量(g)；Drs为组分s的衰变变化速率；NPs表示组分s的直系父组分p的数目；mip为单元i中父组分p的质量(g)；Drp为组分p的衰变变化速率(1/a)；fp-s表示衰变为组分s的父组分p的百分比；Sp-s表示每摩尔组分p衰变产生组分s的摩尔化学计量比；Ws表示组分s的分子量/原子量(又叫摩尔质量，数值上等于相对分子质量)(g/mol)；Wp表示父组分p的摩尔质量(g/mol)；NEi为与单元i相连接的质量通量连接数量；fcs为通过质量通量连接c的组分s的流入通量。
含水通道迁移模型实质上是孔隙单元网格以适当的方式链接而成的含水通道连续流通道。因此，含水通道模型的求解算法与孔隙单元迁移模型相同，采用向后差分的时间隐式差分对含水通道模型进行求解。
3、单裂隙迁移模型
单裂隙迁移模型是对自然条件下某裂隙中水流运动与核素运移的简化描述，主要假设核素在平行板裂隙中流动，并遵循达西定律和质量守恒。在模拟过程中考虑以下过程：在裂隙系统中，核素运移包含对流迁移、分子扩散、纵向机械弥散等机制；核素从裂隙区向基质区扩散；裂隙区和基质区填充物对核素的吸附作用；核素具有放射性衰变等。其中对流过程通过质量守恒建立，分子扩散弥散遵循Fick扩散定律，即单位时间单位面积上的浓度通量与浓度梯度呈正比。吸附过程的简化为等温平衡吸附，并且在公式中用阻滞系数(延迟因子)来代表吸附过程。放射性衰变过程均考虑母体衰变量和自身衰变量，通过求解方程组的形式来简化获取母体衰变量的过程。
根据质量守恒定律，通道末端流动流体中组分浓度数学通式为：
式中，FK表示组分k在裂隙通道出口端的通量(kg/s)；q表示裂隙通道中的体积流速(m3/s)；ck表示裂隙通道中可流动区域内单位有效饱和孔隙空间中的组分k在活动区多孔介质的有效孔隙中的平均浓度(kg/m3)；Sp表示可流动区域的横截面积(m2)；Dk表示组分k可流动区内的有效扩散系数(m2/s)；α表示通道的弥散度(m)；L表示通道的长度(m)；x表示通道内的位置(m，x＝L表示通道出口端位置)。
使用拉普拉斯变换求解单裂隙迁移模型(6)。拉普拉斯变换解决方案技术的中心概念是对于对流占优的系统，排放过程是输入过程的线性函数。
针对具有带有任意数量和类型存储区域的移动区域，开发了管状通道传递矩阵的一般解决方案。这是通过使用传递矩阵Fz表示从裂隙区域到基质扩散区域的通量，通过参考一般裂隙区域与基质扩散区域之间的交换通量来实现的。
考虑无穷小单位长度裂隙区域的体积为S。假定溶质在每个横截面上均匀混合，则浓度可以表示为沿路径x的距离的函数。可移动(溶解的物质)的总量为SpSAmpc，并且总质量为SpSAmpARc。
用于移动区域质量守恒的一般微分方程定义为(总质量变化率)＝(由于衰减引起的变化率)+(由于水流通量引起的变化率)+(由于流出到存储区引起的变化率)：
其中Amp为可移动区域中多孔介质可用孔隙度的对角矩阵；AR为活动区域中延迟因子的对角矩阵，即活动区域中总质量与活动质量之比；RD为衰变/生长速率因子矩阵(T-1)；vq为达西流速(L/T)；Ad为有效扩散系数对角矩阵(L2T-1)；d为弥散系数，等于αq(L2/T)；I为悬浮固体速度放大因子的对角矩阵；z为通道的存储区数；Rz为从流动区到不移动区编号z的质量转移速率(M/T)。
因此，所述步骤103在所述概化后核素迁移过程的基础上，建立核素迁移过程的数学模型，具体包括：
在所述概化后核素迁移过程的基础上，建立所述孔隙单元迁移模型
其中M′mi为网格m中核素i的质量增加率；Mmi为网格m中核素i的质量；Dri为核素i的衰变速率；NJi表示核素i的直系母核素的数量；Mmj为网格m中母核素j的质量；Drj为母核素j的衰变速率；Pj-i表示衰变为i的母核素j的比例；Sj-i表示每摩尔j衰变产生i的摩尔化学计量比；Wi表示核素i的摩尔质量；Wj表示母核素j的摩尔质量；NEm为与单元m相连接的质量通量连接数量；fci为通过质量通量连接l的核素i的流入通量；Wmi为外部源汇项对网格m中核素i的直接补给量；
在所述概化后核素迁移过程的基础上，建立所述含水通道迁移模型
其中m′is为单元i中组分s的质量增加率；mis为单元i中组分s的质量；Drs为组分s的衰变变化速率；NPs表示组分s的直系父组分p的数目；mip为单元i中父组分p的质量；Drp为父组分p的衰变变化速率；fp-s表示衰变为组分s的父组分p的百分比；Sp-s表示每摩尔组分p衰变产生组分s的摩尔化学计量比；Ws表示组分s的摩尔质量；Wp表示父组分p的摩尔质量；NEi为与单元i相连接的质量通量连接c的数量；fcs为通过质量通量连接c的组分s的流入通量；
在所述概化后核素迁移过程的基础上，建立所述单裂隙迁移模型
其中Fk表示组分k在裂隙通道出口端的通量；q表示裂隙通道中的体积流速；ck表示裂隙通道中可流动区域内单位有效饱和孔隙空间中的组分k在活动区多孔介质的有效孔隙中的平均浓度；Sp表示可流动区域的横截面积；Dk表示组分k可流动区内的有效扩散系数；α表示通道的弥散度；L表示通道的长度；x表示通道内的位置。
其中孔隙单元迁移模型可以实现对核素衰变、扩散、吸附等过程的模拟，应用对象包括源项、工程屏障、天然屏障、生物圈等；含水通道迁移模型可以实现对核素衰变、对流、扩散、弥散等过程的模拟，应用对象主要是导水性极低的围岩；单裂隙迁移模型可以实现对核素衰变、对流、弥散、基质扩散、吸附等过程的模拟，应用对象主要是具有一定导水性的围岩裂隙。通过对各类型核素迁移模型的组合表达，可实现对核素在源项、工程屏障、天然围岩和生物圈中的迁移过程的模拟。
所述步骤103分别构建了与孔隙单元、含水通道、单裂隙等迁移模型相对应的主控方程和求解算法，实现不同类型迁移模型的数值表达，可以对核素在不同固体介质单元中的核素迁移的结果进行定量化表达。对于放射性废物处置设施，核素从源项中释放出来后，扩散迁移经过工程屏障到达天然屏障，并迁移通过天然屏障到达生物圈，对上述迁移过程及其结果进行计算与分析，是定量表达处置设施对人类和自然环境的潜在危害的核心。
步骤104：采用所述核素迁移过程的数学模型对所述核素在所述放射性废物处置设施多重屏障系统不同单元中的迁移过程进行模拟。
根据孔隙介质和不同导水性裂隙介质中核素迁移模拟过程的差异性，对核素迁移模型进行模块化划分，用孔隙单元迁移模型对核素在孔隙介质中的迁移过程进行模拟，实现对核素在孔隙单元中衰变、扩散、吸附等过程的模拟，应用对象包括源项、工程屏障、天然屏障、生物圈等；用含水通道迁移模型对核素在弱导水性裂隙介质中的迁移过程进行模拟，实现对核素衰变、对流、扩散、弥散等过程的模拟，应用对象主要是弱导水性的围岩；用单裂隙迁移模型对核素在导水性裂隙中迁移过程进行模拟，实现对核素衰变、对流、弥散、基质扩散、吸附等过程的模拟，应用对象主要是具有一定导水性的裂隙围岩。
因此，所述步骤104采用所述核素迁移过程的数学模型对所述核素在所述放射性废物处置设施多重屏障系统不同单元中的迁移过程进行模拟，具体包括：
采用所述孔隙单元迁移模型模拟核素在孔隙介质中的迁移过程；所述孔隙介质包括源项、工程屏障、天然屏障、生物圈；
采用所述含水通道迁移模型模拟核素在弱导水性裂隙介质中的迁移过程；所述弱导水性裂隙介质包括导水性极低的围岩；
采用所述单裂隙迁移模型模拟核素在导水性裂隙介质中的迁移过程；所述导水性裂隙介质包括具有一定导水性的围岩裂隙。
本发明涉及的一种放射性废物处置过程中的核素迁移模拟方法，可以用孔隙单元模型对核素在工程屏障中的迁移过程进行模拟，得到的计算结果作为含水通道模型的输入；接着用含水通道模型对核素在工程屏障外围弱导水的裂隙介质中的迁移过程进行模拟，得到的计算结果作为单裂隙模型的输入；然后用单裂隙模型对核素在导水性裂隙介质中的迁移过程进行模拟，得到的计算结果作为生物圈模拟的输入条件。最终实现对核素在放射性废物处置设施多重屏障系统中全过程的模拟，为放射性废物处置设施长期安全评价提供可行的解决方案。本发明提供的方法能够有效的实现对万年时间尺度上核素迁移过程的模拟，对推进放射性废物地质处置安全评价具有重要意义。
基于本发明提供的一种放射性废物处置过程中的核素迁移模拟方法，本发明还提供一种放射性废物处置过程中的核素迁移模拟系统，参见图5，所述系统包括：
核素迁移过程描述模块501，用于对放射性废物处置设施多重屏障系统中的核素迁移过程进行描述，确定核素在所述放射性废物处置设施多重屏障系统不同单元中的迁移过程；所述放射性废物处置设施多重屏障系统不同单元包括源项、工程屏障、天然屏障和生物圈；
核素迁移过程概化模块502，用于结合万年以上的安全评价要求，根据所述核素在放射性废物处置设施多重屏障系统不同单元中迁移过程的差异性，对核素迁移过程进行概化；
核素迁移过程数学模型建立模块503，用于在所述概化后核素迁移过程的基础上，建立核素迁移过程的数学模型；所述核素迁移过程的数学模型包括孔隙单元迁移模型、含水通道迁移模型和单裂隙迁移模型；
核素迁移过程模拟模块504，用于采用所述核素迁移过程的数学模型对所述核素在所述放射性废物处置设施多重屏障系统不同单元中的迁移过程进行模拟。
其中，所述核素迁移过程描述模块501具体包括：
核素迁移过程描述单元，用于对放射性废物处置设施多重屏障系统中核素迁移过程进行描述，确定核素在所述放射性废物处置设施多重屏障系统不同单元中的迁移过程；所述迁移过程包括放射性衰变的物理过程以及核素在所述放射性废物处置设施多重屏障系统不同单元中迁移的近场过程和远场过程。
所述核素迁移过程概化模块502具体包括：
核素迁移过程概化单元，用于结合万年以上的安全评价要求，根据所述核素在放射性废物处置设施多重屏障系统不同单元中迁移过程的差异性，将核素在孔隙介质中以扩散形式迁移为主的迁移过程概化为核素随地下水在均匀连续的孔隙介质中传输；将核素在弱导水性裂隙介质中以对流-扩散形式为主的迁移过程概化为核素随地下水在裂隙介质中传输的一维连续通道；将核素在导水性裂隙介质中以对流-弥散形式为主的迁移过程概化为核素随地下水在裂隙介质中传输的一维通道。
所述核素迁移过程数学模型建立模块503具体包括：
孔隙单元迁移模型建立单元，用于在所述概化后核素迁移过程的基础上，建立所述孔隙单元迁移模型
其中M′mi为网格m中核素i的质量增加率；Mmi为网格m中核素i的质量；Dri为核素i的衰变速率；NJi表示核素i的直系母核素的数量；Mmj为网格m中母核素j的质量；Drj为母核素j的衰变速率；Pj-i表示衰变为i的母核素j的比例；Sj-i表示每摩尔j衰变产生i的摩尔化学计量比；Wi表示核素i的摩尔质量；Wj表示母核素j的摩尔质量；NEm为与单元m相连接的质量通量连接数量；fci为通过质量通量连接l的核素i的流入通量；Wmi为外部源汇项对网格m中核素i的直接补给量；
含水通道迁移模型建立单元，用于在所述概化后核素迁移过程的基础上，建立所述含水通道迁移模型
其中m′is为单元i中组分s的质量增加率；mis为单元i中组分s的质量；Drs为组分s的衰变变化速率；NPs表示组分s的直系父组分p的数目；mip为单元i中父组分p的质量；Drp为父组分p的衰变变化速率；fp-s表示衰变为组分s的父组分p的百分比；Sp-s表示每摩尔组分p衰变产生组分s的摩尔化学计量比；Ws表示组分s的摩尔质量；Wp表示父组分p的摩尔质量；NEi为与单元i相连接的质量通量连接c的数量；fcs为通过质量通量连接c的组分s的流入通量；
单裂隙迁移模型建立单元，用于在所述概化后核素迁移过程的基础上，建立所述单裂隙迁移模型
其中Fk表示组分k在裂隙通道出口端的通量；q表示裂隙通道中的体积流速；ck表示裂隙通道中可流动区域内单位有效饱和孔隙空间中的组分k在活动区多孔介质的有效孔隙中的平均浓度；Sp表示可流动区域的横截面积；Dk表示组分k可流动区内的有效扩散系数；α表示通道的弥散度；L表示通道的长度；x表示通道内的位置。
所述核素迁移过程模拟模块504具体包括：
孔隙介质迁移过程模拟单元，用于采用所述孔隙单元迁移模型模拟核素在孔隙介质中的迁移过程；所述孔隙介质包括源项、工程屏障、天然屏障、生物圈；
弱导水性裂隙介质迁移过程模拟单元，用于采用所述含水通道迁移模型模拟核素在弱导水性裂隙介质中的迁移过程；所述弱导水性裂隙介质包括导水性极低的围岩；
导水性裂隙介质迁移过程模拟单元，用于采用所述单裂隙迁移模型模拟核素在导水性裂隙介质中的迁移过程；所述导水性裂隙介质包括具有一定导水性的围岩裂隙。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。
Claims (8)
1.一种放射性废物处置过程中的核素迁移模拟方法，其特征在于，包括：
对放射性废物处置设施多重屏障系统中的核素迁移过程进行描述，确定核素在所述放射性废物处置设施多重屏障系统不同单元中的迁移过程；所述放射性废物处置设施多重屏障系统不同单元包括源项、工程屏障、天然屏障和生物圈；
结合万年以上的安全评价要求，根据所述核素在放射性废物处置设施多重屏障系统不同单元中迁移过程的差异性，对核素迁移过程进行概化；
所述结合万年以上的安全评价要求，根据所述核素在放射性废物处置设施多重屏障系统不同单元中迁移过程的差异性，对核素迁移过程进行概化，具体包括：
结合万年以上的安全评价要求，根据所述核素在放射性废物处置设施多重屏障系统不同单元中迁移过程的差异性，将核素在孔隙介质中的迁移过程概化为核素在均匀连续的孔隙单元间的迁移过程，将核素在弱导水性裂隙介质中以对流-扩散形式为主的迁移过程概化为核素随地下水在一维连续通道中的迁移过程；将核素在导水性裂隙介质中以对流-弥散形式为主的迁移过程概化为核素随地下水在一维裂隙通道中的迁移过程；
在所述概化后核素迁移过程的基础上，建立核素迁移过程的数学模型；所述核素迁移过程的数学模型包括孔隙单元迁移模型、含水通道迁移模型和单裂隙迁移模型；
采用所述核素迁移过程的数学模型对所述核素在所述放射性废物处置设施多重屏障系统不同单元中的迁移过程进行模拟。
2.根据权利要求1所述的方法，其特征在于，所述对放射性废物处置设施多重屏障系统中核素迁移过程进行描述，确定核素在所述放射性废物处置设施多重屏障系统不同单元中的迁移过程，具体包括：
对放射性废物处置设施多重屏障系统中核素迁移过程进行描述，确定核素在所述放射性废物处置设施多重屏障系统不同单元中的迁移过程；所述迁移过程包括放射性衰变的物理过程以及核素在所述放射性废物处置设施多重屏障系统不同单元中迁移的近场过程和远场过程。
3.根据权利要求2所述的方法，其特征在于，所述在所述概化后核素迁移过程的基础上，建立核素迁移过程的数学模型，具体包括：
在所述概化后核素迁移过程的基础上，建立所述孔隙单元迁移模型
其中M′mi为网格m中核素i的质量增加率；Mmi为网格m中核素i的质量；Dri为核素i的衰变速率；NJi表示核素i的直系母核素的数量；Mmj为网格m中母核素j的质量；Drj为母核素j的衰变速率；Pj-i表示衰变为i的母核素j的比例；Sj-i表示每摩尔j衰变产生i的摩尔化学计量比；Wi表示核素i的摩尔质量；Wj表示母核素j的摩尔质量；NEm为与单元m相连接的质量通量连接数量；fci为通过质量通量连接l的核素i的流入通量；Wmi为外部源汇项对网格m中核素i的直接补给量；
在所述概化后核素迁移过程的基础上，建立所述含水通道迁移模型
其中m′is为单元i中组分s的质量增加率；mis为单元i中组分s的质量；Drs为组分s的衰变变化速率；NPs表示组分s的直系父组分p的数目；mip为单元i中父组分p的质量；Drp为父组分p的衰变变化速率；fp-s表示衰变为组分s的父组分p的百分比；Sp-s表示每摩尔组分p衰变产生组分s的摩尔化学计量比；Ws表示组分s的摩尔质量；Wp表示父组分p的摩尔质量；NEi为与单元i相连接的质量通量连接c的数量；fcs为通过质量通量连接c的组分s的流入通量；
在所述概化后核素迁移过程的基础上，建立所述单裂隙迁移模型
其中Fk表示组分k在裂隙通道出口端的通量；q表示裂隙通道中的体积流速；ck表示裂隙通道中可流动区域内单位有效饱和孔隙空间中的组分k在活动区多孔介质的有效孔隙中的平均浓度；Sp表示可流动区域的横截面积；Dk表示组分k可流动区内的有效扩散系数；α表示通道的弥散度；L表示通道的长度；x表示通道内的位置。
4.根据权利要求3所述的方法，其特征在于，所述采用所述核素迁移过程的数学模型对所述核素在所述放射性废物处置设施多重屏障系统不同单元中的迁移过程进行模拟，具体包括：
采用所述孔隙单元迁移模型模拟核素在孔隙介质中的迁移过程；所述孔隙介质包括源项、工程屏障、天然屏障、生物圈；
采用所述含水通道迁移模型模拟核素在弱导水性裂隙介质中的迁移过程；所述弱导水性裂隙介质包括导水性极低的围岩；
采用所述单裂隙迁移模型模拟核素在导水性裂隙介质中的迁移过程；所述导水性裂隙介质包括具有一定导水性的围岩裂隙。
5.一种放射性废物处置过程中的核素迁移模拟系统，其特征在于，包括：
核素迁移过程描述模块，用于对放射性废物处置设施多重屏障系统中核素迁移过程进行描述，确定核素在所述放射性废物处置设施多重屏障系统不同单元中的迁移过程；所述放射性废物处置设施多重屏障系统不同单元包括源项、工程屏障、天然屏障和生物圈；
核素迁移过程概化模块，用于结合万年以上的安全评价要求，根据所述核素在放射性废物处置设施多重屏障系统不同单元中迁移过程的差异性，对核素迁移过程进行概化；
所述核素迁移过程概化模块具体包括：
核素迁移过程概化单元，用于结合万年以上的安全评价要求，根据所述核素在放射性废物处置设施多重屏障系统不同单元中迁移过程的差异性，将核素在孔隙介质中的迁移过程概化为核素在均匀连续的孔隙单元间的迁移过程，将核素在弱导水性裂隙介质中以对流-扩散形式为主的迁移过程概化为核素随地下水在一维连续通道中的迁移过程；将核素在导水性裂隙介质中以对流-弥散形式为主的迁移过程概化为核素随地下水在一维裂隙通道中的迁移过程；
核素迁移过程数学模型建立模块，用于在所述概化后核素迁移过程的基础上，建立核素迁移过程的数学模型；所述核素迁移过程的数学模型包括孔隙单元迁移模型、含水通道迁移模型和单裂隙迁移模型；
核素迁移过程模拟模块，用于采用所述核素迁移过程的数学模型对所述核素在所述放射性废物处置设施多重屏障系统不同单元中的迁移过程进行模拟。
6.根据权利要求5所述的系统，其特征在于，所述核素迁移过程描述模块具体包括：
核素迁移过程描述单元，用于对放射性废物处置设施多重屏障系统中核素迁移过程进行描述，确定核素在所述放射性废物处置设施多重屏障系统不同单元中的迁移过程；所述迁移过程包括放射性衰变的物理过程以及核素在所述放射性废物处置设施多重屏障系统不同单元中迁移的近场过程和远场过程。
7.根据权利要求6所述的系统，其特征在于，所述核素迁移过程数学模型建立模块具体包括：
孔隙单元迁移模型建立单元，用于在所述概化后核素迁移过程的基础上，建立所述孔隙单元迁移模型
其中M′mi为网格m中核素i的质量增加率；Mmi为网格m中核素i的质量；Dri为核素i的衰变速率；NJi表示核素i的直系母核素的数量；Mmj为网格m中母核素j的质量；Drj为母核素j的衰变速率；Pj-i表示衰变为i的母核素j的比例；Sj-i表示每摩尔j衰变产生i的摩尔化学计量比；Wi表示核素i的摩尔质量；Wj表示母核素j的摩尔质量；NEm为与单元m相连接的质量通量连接数量；fci为通过质量通量连接l的核素i的流入通量；Wmi为外部源汇项对网格m中核素i的直接补给量；
含水通道迁移模型建立单元，用于在所述概化后核素迁移过程的基础上，建立所述含水通道迁移模型
其中m′is为单元i中组分s的质量增加率；mis为单元i中组分s的质量；Drs为组分s的衰变变化速率；NPs表示组分s的直系父组分p的数目；mip为单元i中父组分p的质量；Drp为父组分p的衰变变化速率；fp-s表示衰变为组分s的父组分p的百分比；Sp-s表示每摩尔组分p衰变产生组分s的摩尔化学计量比；Ws表示组分s的摩尔质量；Wp表示父组分p的摩尔质量；NEi为与单元i相连接的质量通量连接c的数量；fcs为通过质量通量连接c的组分s的流入通量；
单裂隙迁移模型建立单元，用于在所述概化后核素迁移过程的基础上，建立所述单裂隙迁移模型
其中Fk表示组分k在裂隙通道出口端的通量；q表示裂隙通道中的体积流速；ck表示裂隙通道中可流动区域内单位有效饱和孔隙空间中的组分k在活动区多孔介质的有效孔隙中的平均浓度；Sp表示可流动区域的横截面积；Dk表示组分k可流动区内的有效扩散系数；α表示通道的弥散度；L表示通道的长度；x表示通道内的位置。
8.根据权利要求7所述的系统，其特征在于，所述核素迁移过程模拟模块具体包括：
孔隙介质迁移过程模拟单元，用于采用所述孔隙单元迁移模型模拟核素在孔隙介质中的迁移过程；所述孔隙介质包括源项、工程屏障、天然屏障、生物圈；
弱导水性裂隙介质迁移过程模拟单元，用于采用所述含水通道迁移模型模拟核素在弱导水性裂隙介质中的迁移过程；所述弱导水性裂隙介质包括导水性极低的围岩；
导水性裂隙介质迁移过程模拟单元，用于采用所述单裂隙迁移模型模拟核素在导水性裂隙介质中的迁移过程；所述导水性裂隙介质包括具有一定导水性的围岩裂隙。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210377800.XA CN114707337B (zh) | 2022-04-12 | 2022-04-12 | 一种放射性废物处置过程中的核素迁移模拟方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210377800.XA CN114707337B (zh) | 2022-04-12 | 2022-04-12 | 一种放射性废物处置过程中的核素迁移模拟方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114707337A CN114707337A (zh) | 2022-07-05 |
| CN114707337B true CN114707337B (zh) | 2022-11-29 |
Family
ID=82173390
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210377800.XA Active CN114707337B (zh) | 2022-04-12 | 2022-04-12 | 一种放射性废物处置过程中的核素迁移模拟方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114707337B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116384045B (zh) * | 2023-01-16 | 2023-12-05 | 中核第四研究设计工程有限公司 | 采用三维数值模拟确定土壤Ra-226源项调查范围的方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106815390A (zh) * | 2016-05-16 | 2017-06-09 | 中国辐射防护研究院 | 高放废物玻璃固化体核素迁移水解反应动力学模型的构建方法 |
| CN109325017A (zh) * | 2018-10-25 | 2019-02-12 | 中国辐射防护研究院 | 一种放射性核素迁移扩散数值模拟系统的数据库设计方法 |
| CN109614634A (zh) * | 2018-10-25 | 2019-04-12 | 中国辐射防护研究院 | 一种水环境中放射性核素迁移扩散数值模拟方法及系统 |
| CN110555215A (zh) * | 2018-05-30 | 2019-12-10 | 中国辐射防护研究院 | 一种近场核素迁移模型中edz区域的等效水流量模型 |
| CN113155679A (zh) * | 2021-05-07 | 2021-07-23 | 吉林大学 | 岩体裂隙中放射性核素吸附运移模拟的实验装置及实验方法 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8298513B2 (en) * | 2007-02-05 | 2012-10-30 | Osaka University | Hexatriene-β-carbonyl compound |
| CN106815460A (zh) * | 2016-10-11 | 2017-06-09 | 中国辐射防护研究院 | 一种离散裂隙网络评价放射性核素在岩石裂隙中迁移的方法 |
| CN107704682B (zh) * | 2017-09-30 | 2021-08-10 | 西南科技大学 | 基于概率用于核素近场、远场迁移评估的空间域描述方法 |
| CN113866355B (zh) * | 2021-09-13 | 2023-10-10 | 中国人民解放军63653部队 | 处置库多重屏障中水岩作用与核素运移模拟实验方法 |
-
2022
- 2022-04-12 CN CN202210377800.XA patent/CN114707337B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106815390A (zh) * | 2016-05-16 | 2017-06-09 | 中国辐射防护研究院 | 高放废物玻璃固化体核素迁移水解反应动力学模型的构建方法 |
| CN110555215A (zh) * | 2018-05-30 | 2019-12-10 | 中国辐射防护研究院 | 一种近场核素迁移模型中edz区域的等效水流量模型 |
| CN109325017A (zh) * | 2018-10-25 | 2019-02-12 | 中国辐射防护研究院 | 一种放射性核素迁移扩散数值模拟系统的数据库设计方法 |
| CN109614634A (zh) * | 2018-10-25 | 2019-04-12 | 中国辐射防护研究院 | 一种水环境中放射性核素迁移扩散数值模拟方法及系统 |
| CN113155679A (zh) * | 2021-05-07 | 2021-07-23 | 吉林大学 | 岩体裂隙中放射性核素吸附运移模拟的实验装置及实验方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114707337A (zh) | 2022-07-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Černý | Modeling of radionuclide transport in porous media: A review of recent studies | |
| CN114707337B (zh) | 一种放射性废物处置过程中的核素迁移模拟方法及系统 | |
| Adinarayana et al. | Modelling of calcium leaching and its influence on radionuclide migration across the concrete engineered barrier in a NSDF | |
| Du et al. | Uncertainty and sensitivity analysis of radionuclide migration through fractured granite aquifer | |
| Samper et al. | Inverse modeling of tracer experiments in FEBEX compacted Ca-bentonite | |
| Kim et al. | A modularized numerical framework for the process-based total system performance assessment of geological disposal systems | |
| Cvetkovic | Statistical formulation of generalized tracer retention in fractured rock | |
| Kelkar et al. | Modeling solute transport through saturated zone ground water at 10 km scale: Example from the Yucca Mountain license application | |
| Zuber | Review of existing mathematical models for interpretation of tracer data in hydrology | |
| Malkovsky et al. | Numerical analysis of safety of a borehole repository for vitrified high-level nuclear waste | |
| Uchida et al. | An empirical probabilistic approach for constraining the uncertainty of long-term solute transport predictions in fractured rock using in situ tracer experiments | |
| Cadini et al. | Monte Carlo simulation of radionuclide migration in fractured rock for the performance assessment of radioactive waste repositories | |
| Cvetkovic et al. | Comparative Measures of Radionuclide Containment in the Crystalline Geophere | |
| Russo | Prediction of the migration of several radionuclides in ocean sediment with the computer code IONMIG: a preliminary report | |
| CN114722608B (zh) | 一种核素在裂隙介质中衰变迁移的模拟方法及系统 | |
| Elert et al. | Assessment model validity document FARF31 | |
| Onishi et al. | Preliminary Three-Dimensional Simulation of Sediment and Cesium Transport in the Ogi Dam Reservoir using FLESCOT–Task 6, Subtask 2 | |
| LaForce et al. | DECOVALEX-2023 Task F Specification Revision 9. | |
| Rasmussen | Solute transport in saturated fractured media | |
| Robinson | Particle tracking model and abstraction of transport processes | |
| CN114674716A (zh) | 一种核素在孔隙单元中衰变迁移的模拟方法及系统 | |
| Wu et al. | A 3-D hydrodynamic dispersion model for modeling tracer transport in Geothermal Reservoirs | |
| Migranova et al. | MATHEMATICAL MODELING OF RADIONUCLIDES SPREADING IN SOIL WITH DECAY | |
| Lee et al. | Current status of KURT and its long-term experimental research programme | |
| Shahkarami et al. | Simulation of Helium transport in fractured rocks: Implementation of a dual continuum model in DarcyTools |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |