不确定性分析技术在某环保配套工程安全评价中的应用研究.pdf
《不确定性分析技术在某环保配套工程安全评价中的应用研究.pdf》由会员分享,可在线阅读,更多相关《不确定性分析技术在某环保配套工程安全评价中的应用研究.pdf(11页珍藏版)》请在咨信网上搜索。
1、第 43 卷第 4 期2023 年7 月辐射防护Radiation ProtectionVol.43No.4July2023辐射防护评价不确定性分析技术在某环保配套工程安全评价中的应用研究刘兴伟,王旭宏,吕涛,李星宇,李昶,康宝伟,王馨,夏加国(中国核电工程有限公司,北京 100840)摘要:为解决放射性废物处置安全评价结果受认知水平和长时间尺度等因素引起的不确定性问题,基于某环保配套工程项目开展确定性模型安全评价,然后基于概率论方法、拉丁超立方抽样技术生成 1 000 个样本,开展不确定性分析和灵敏度分析工作。结果表明确定性模型相对于不确定性模型,对近场释放率计算结果偏高,对地质圈释放率计算
2、结果良好,对 Mo-93、Ni-59 产生的照射剂量估算较好,对 I-129 产生的照射剂量代表性中等偏下。灵敏度分析结果表明工程屏障中混凝土相关参数(如:混凝土有效扩散系数、混凝土中核素分配系数、混凝土屏障厚度)、废物浸出率、初始活度和核素在地质圈的分配系数为模型释放率峰值主要影响参数,可为后续工程设计与野外调查遴选出重点参数。关键词:放射性废物处置;不确定性分析;灵敏度分析;岩洞型处置;库室模型中图分类号:TL942+.1文献标识码:A收稿日期:2022-04-26作者简介:刘兴伟(1990),男,2013 年本科毕业于中国地质大学(北京)水文与水资源工程专业,2018 年博士毕业于中国地
3、质大学(北京)水文地质学专业(硕博连读),高级工程师。E-mail:liuxingwei0912 放射性废物处置是完整核工业产业不可缺失的一环,相应的放射性废物处置设施长期安全日益受到关注,通常采用安全评价方法论证放射性废物处置设施的长期安全。安全评价的时间尺度应涵盖处置设施对放射性废物的包容、隔离期。我国相关法规规定:对于低、中水平放射性废物包容,隔离期不少于 300 年1,对于高水平放射性废物包容,隔离期不少于 1 万年2。长时间尺度安全评价受自然系统演变、人类现有知识水平、研究程度和参数尺度效应等多方面影响,不可避免的存在不确定性3,不确定性分析的好坏将直接影响安全评价结果的可信程度。国
4、际原子能机构(IAEA)相关安全标准 SSG-234、通用安全要求 GSR Part4 5及我国低、中水平放射性固体废物近地表处置安全规定(GB 91322018)1、放射性废物处置安全全过程系统分析(NNSA-HAJ-00012020)6等法规标准均已将不确定性分析列为安全评价及安全全过程系统分析的重要组成部分。在放射性废物处置领域,通常将不确定性分为:景象不确定性、模型不确定性和参数不确定性6。目前在我国放射性废物处置领域,不确定分析技术虽在理论研究方面引起了部分学者的注意3,7-10,但我国已完成的放射性废物处置设施安全评价中,多采用相关参数的保守值处理不确定性,基于概率论方法的不确定性
5、分析技术在实际工程中完整迁移评价的应用研究相对较少11,对于不确定性安全评价流程及不确定性的展示方式尚未有统一认识。本研究以某环保配套工程为例,采取概率论的方法考虑参数不确定性,对放射性废物处置工程分别进行确定论及概率论方法的安全评价,并进行参数灵敏度及相关性分析,对不确定性安全评价及不确定性展示方式进行初步探索。1模型概述本次研究以某环保配套项目为例开展不确定性安全评价。该环保配套项目主要用于低中放废物最终处置,为地下岩洞型,共有 15 条处置巷道,位于海岛下方-180 m 岩体中。根据勘察资料显示,场址区主要为花岗斑岩及长石斑岩。地下水类型以基岩裂隙水为主,地下水流速较为缓慢,呈523 辐
6、射防护第 43 卷第 4 期“孤岛”模型式地下水系统,处置区域主要接受大气降水及侧向径流补给,最终排至海洋(如图 1所示)。图 1某处置设施关闭后核素迁移途径示意图Fig.1Schematic diagram of nuclide transport path post closure处置设施关闭后,处置巷道使用回填材料(水泥砂浆)充满。在近场放射性核素释放主要以扩散为主,从水泥固化体中浸出进入水泥砂浆及混凝土。考虑处置库内废物桶直接放置于混凝土屏障上,因此近场有两条释出途径:一条为废物桶水泥砂浆周边及顶部混凝土屏障工程扰动区;另一条为废物桶底部混凝土屏障工程扰动区。在地质圈,放射性核素随地下
7、水迁移,主要发生对流和弥散作用。最终进入海洋生态系统,经海水稀释后,对近海区域鱼和甲壳类生物产生污染,从而对食用受污染的海洋生物人体产生内照射。因放射性废物处置工程较为复杂,涉及屏障多、时间跨度大,不可避免的存在不确定性。对于岩洞型放射性废物处置库评价过程中的不确定性同样可总体上分为景象不确定性、模型不确定性和参数不确定性。景象不确定性:工程屏障(如混凝土)可能处于不同的水化学组分演变过程中、不同的地下水流速景象等。可通过特征(features)、事件(events)和过程(processes)进行梳理,从而构建不同的景象,降低不确定性的影响。模型不确定性:核素在裂隙介质中迁移可概化为仅考虑对
8、流、考虑对流和扩散、考虑对流扩散和吸附三种形式,可根据研究精度的需要采用多种模型相互验证降低不确定性影响。参数不确定性:可能来自随机误差、实验尺度效应、测量设备及人员等方面,可通过增加数据量减少随机误差的影响,对于部分参数可采取保守值或概率论的处理方法。本次研究将重点以概率论的方法了解参数不确定性对评价结果的影响。相 关 测 试 使 用Ecolego软 件,是 基 于MATLAB/Simulink 模块开发的用于放射性废物管理安全评价和动态模拟计算软件,软件中内嵌多种解法器,并包含不确定性分析与灵敏度分析工具包,目前已广泛应用于国内外放射性废物处置安全评价中3,12。2模型公式与参数2.1输入
9、源项因废物源项存在较大的不确定性,根据预估关闭后总量均分到 15 条巷道中,并对单条巷道活度进行 50%和 200%的变化。单条巷道关闭时预估源项如表 1 所示。相关初始数据源于项目输入文件,概率密度函数的选取主要依据各类放射性废物处置相关文献。2.2近场迁移公式与参数本次研究将近场定义为废物体、水泥砂浆、混凝土屏障和工程开挖后的工程扰动区(EDZ)。因623刘兴伟等:不确定性分析技术在某环保配套工程安全评价中的应用研究 表 1废物源项核素与活度Tab.1Radionuclide and activity of waste source term核素半衰期(a)单条巷道参考活度(Bq)概率密度
10、函数取值范围(Bq)Co-605.36.921015均匀分布3.4610151.381016Cs-1373.0103.491014均匀分布1.7410146.981014I-1291.61072.49106均匀分布1.251064.98106Sr-902.9107.93109均匀分布3.971091.591010C-145.71031.51109均匀分布7.551083.02109Ni-597.61042.291013均匀分布1.1510134.581013Ni-631.01022.071015均匀分布1.0410154.141015Mo-934.11031.721012均匀分布8.60101
11、13.441012Tc-992.11058.071010均匀分布4.0410101.611011区域地下水流速相对较低,在近场仅考虑核素在各库室间的扩散作用,相关公式见参考文献3,12。工程扰动区(EDZ)的孔隙度看作 1,相关阻滞因子可忽略。由工程扰动区进入地质圈的迁移系数可表示为:EDZGeo=QEDZVEDZ(1)式中,QEDZ为工程扰动区的流量,m3/a;VEDZ为工程扰动区体积,m3。放射性核素由废物桶向工程屏障浸出率取值参考低、中水平放射性废物固化体性能要求 水泥固化体(GB 14569.12011)13中的水泥固化体 42 天浸出率限值(如表 2 所示),近场迁移参数取值见表 3
12、4。表 2核素浸出率取值范围Tab.2Value range of nuclide leaching rate核素参考值(cm/d)概率密度函数取值范围(cm/d)Co-602.010-3对数分布148.610-72.010-3 15C-1374.010-3对数分布141.810-44.010-3 15Sr-901.010-3对数分布141.010-51.010-3 16,17其他核素1.010-5对数分布141.010-61.010-5表 3近场迁移参数 18-19Tab.3Near-field transport parameters 18-19参数名称1)参考值单位概率密度函数取值范围L
13、_G0.02m均匀分布0.020.64L_C1m均匀分布0.51.5n_G0.3均匀分布0.30.5n_C0.11均匀分布0.110.5De_G9.4610-3m2/a对数均匀1.110-23.1510-218De_C9.4610-5m2/a对数均匀203.1510-53.1510-420L_EDZ0.5m均匀分布0.32.7 1)L 表示核素迁移方向上的距离,n 表示孔隙度,De 表示有效扩散系数,G 表示水泥砂浆,C 表示混凝土,EDZ 表示工程扰动区。2.3地质圈迁移公式与参数目前国际上对于裂隙岩体有两种概化方案:一种用等效多孔介质,另一种作为随机裂隙处理。结合现有数据情况,本次研究将岩
14、体作为等效多孔介质处理。核素在地质圈主要发生对流和弥散两种迁移形式。相关迁移公式参见文献3,12,参数取值如表 5 所示。2.4生物圈计算公式与参数放射性核素由地质圈释放进入海洋,首先被海水稀释,后被生活在该海域的鱼类和甲壳类海洋生物富集,该海域周边公众食入受污染的鱼类和甲壳类生物而受到照射。海水中放射性核素浓度Cwki=QiCki/(Adisvrock)(2)式中,Cwki为 k 海域海水中放射性核素 i 的浓度,Bq/m3;Qi为放射性核素 i 经地下水进入海洋中的年释放率,Bq/a;Adis为含核素地下水进入海洋时的排泄面积,m2;Cki 为 k 海域海水中放射性核素 i的稀释因子,无量
15、纲,每年流入到海洋的地下水水量由地下水迁移速度与含水层截面积进行计算,根据以往工程的经验,本次研究稀释因子保守取1.010-5;vrock为地下水渗流速度,m/a。海产品中放射性核素浓度Cpki=CwkiBpi(3)式中,Cpki为 k 海域内的海产品中放射性核素 i 浓度,Bq/kg;Bpi为海产品 p 中放射性核素 i 浓集因子,m3/kg。食入海产品所致的个人内照射剂量Dep=CpkiUpDFei(4)式中,Dep为公众个人食入 k 海域内海产品 p 所致723 辐射防护第 43 卷第 4 期表 4核素分配系数 KdTab.4The Kd of nuclides材质水泥砂浆混凝土地质圈参
16、数名称参考值(m3/kg)取值范围(m3/kg)参考值(m3/kg)取值范围(m3/kg)概率密度函数1)参考值(m3/kg)取值范围(m3/kg)概率密度函数C-146.9410-1min=0.243,max=6.94,mode=0.942.1910-1min=0.153,max=4.38,mode=0.438对数三角05.010-45.010-3 21对数均匀Co-602.0810-3min=0.00139,max=0.139,mode=0.01392.010-2min=0.0035,max=0.088,mode=0.0088对数三角209.310-32.110-52.710-2对数均匀2
17、0Ni-591.0410-2min=6.9410-3,max=0.694,mode=0.032.010-2min=0.018,max=0.44,mode=0.044对数三角205.510-32.410-47.010-2对数均匀20Ni-631.0410-2min=6.9410-3,max=0.694,mode=0.032.010-2min=0.018,max=0.44,mode=0.044对数三角201.210-21.210-22.610-2对数均匀20Sr-909.810-4min=3.4710-5,max=6.94,mode=3.4710-45.010-3min=0.0011,max=0.
18、022,mode=0.0066对数三角201.210-41.210-45.510-3对数均匀20Mo-931.0410-3min=1.0410-4,max=1.1510-2,mode=1.0410-34.610-4min=6.5710-5,max=0.00723,mode=6.5710-4对数三角3.010-49.410-99.210-4对数均匀20Tc-993.4710-4min=3.4710-5,max=3.4710-3,mode=3.4710-45.010-4min=2.210-5,max=4.4,mode=2.210-4对数三角204.010-12.010-28.0对数均匀20I-12
19、93.4710-4min=0.0,max=3.4710-3,mode=3.4710-40.0min=0.0,max=0.0022,mode=2.210-4三角分布2001.010-103.010-6对数均匀20Cs-1376.9410-4min=3.4710-5,max=0.104,mode=6.9410-41.010-4min=2.210-5,max=0.011,mode=4.410-4对数三角205.410-22.310-22.110-1对数均匀201)水泥砂浆和混凝土采用相同形态的概率密度函数,未特别注释相关参数范围引自项目勘察报告。表 5地质圈迁移参数Tab.5Geosphere tr
20、ansport parameters参数名称参考值单位概率密度函数取值范围流速 v_Rock7.310-2m/a均匀分布3.410-21.610-1纵向弥散度 L1.110-2m均匀分布1.110-31.110-2岩体孔隙度 e0.02正态分布均值=0.02,标准差=0.005地质圈迁移距离 L地质500m均匀分布320954的有效剂量,Sv/a;Cpki为 k 海域内的海产品中放射性核素 i 浓度,Bq/kg;Up为公众个人的海产品 p消费量,kg/a;DFei为因食入海产品 p 放射性核素i 对公众个人的有效剂量转换因子,Sv/Bq。含核素地下水进入海洋时的排泄面积取均匀分布,参考值 4
21、000 m2,最小值 3 000 m2,最大值5 000 m2;鱼 类 与 甲 壳 类 生 物 浓 集 因 子 引 用Resrad 6.0 参 数22,食 入 海 产 品 质 量 分 别 取:5.4 kg/a(鱼类),0.9 kg/a(甲壳类)23。有效剂量转换因子采用电离辐射防护与辐射源安全基本标准GB 188712002 中成人组食入转换因子数据。2.5运行参数设置首先进 行确定性模 型运算,相关 参数采用表 15 中的参考值,后续分析称为确定性模型结果。随后开展不确定性及灵敏度分析试验,相关参数采用随 机种子方法 生成,使用拉 丁超立方(Latin Hypercube)方法进行采样,共生
22、成 1 000 组参数,使用 TR-BDF2 解法器进行求解24。823刘兴伟等:不确定性分析技术在某环保配套工程安全评价中的应用研究 3结果与讨论3.1确定性结果分析3.1.1近场释放率确定性分析图 2 为各核素在确定性模型的近场释放率。可以看出,处置设施关闭后 700 年近场主要释放核素为 Cs-137,700 年至 4 100 年主要释放核素为Mo-93 和 Ni-59,其中以 Mo-93 为主。4 100 年后Ni-59 近场释放率继续上升,约在 80 000 年达到峰值,为近场区域释放的最大峰值 1.06107 Bq/a。近场中 Co-60 在处置设施关闭 20 年前后首先达到释放率
23、峰值。图 2确定性模型近场释放率Fig.2The release rate of near-field ofthe deterministic model3.1.2地质圈释放率确定性分析图 3 为各核素在地质圈释放率的确定性分析结果。相对于近场释放率来说,地质圈释放率数值明显降低,约下降 2 个数量级。同时,释出的核素种类也有所减少。在 300 年至 8 000 年,地质圈主要释出核素为 I-129,8 000 至 80 000 年间地质圈主要释出核素为 Mo-93,随后为 Ni-59。地质圈释放峰值发生在 32 万年前后,滞后于近场释放率峰值约 24 万年。地质圈最早达到释放率峰值的核素为
24、C-14,释放率峰值约发生在 1.5 万年,约滞后于 C-14 近场释放率峰值 400 年,说明地质圈对核素有较好的阻滞作用。I-129 最先释出,主要由于地质圈对 I-129 的吸附能力较弱,如表 4 所示核素分配系数为 0。由图 3 可以看出地质圈释出核素均为长衰变周期核素,主要原因一方面为地质圈吸附阻滞作用,另一方面由于核素自身衰变原因,短寿命核素随着时间推移衰变至很低的活度,因此地质圈释出均为长寿命核素。图 3地质圈释放率确定性分析结果Fig.3The release rate of geosphere ofthe deterministic model3.1.3生物圈照射剂量确定性分
25、析图 4 为各核素确定性模型的有效剂量。可以看出,照射剂量出现两次峰值,第一次峰值出现在处置设施关闭后 2 万年左右,第二次峰值出现在处置设施关闭后 32 万年左右,且高于第一次峰值。两次峰值的主要贡献核素不同,第一次峰值主要贡献核素为 Mo-93,第二次峰值主要贡献核素为 Ni-59。两次峰值均远小于低、中水平放射性固体废物处置安全规定(GB 91322018)中对公众有效照射剂量 0.25 mSv 的限值标准,总照射剂量峰值与剂量限值相差约 8 个数量级。同时可以看出,关键居民组食入鱼类产生的受照剂量均高于食入甲壳类食物产生的照射剂量。其中,食入含 I-129 的鱼类产生的照射剂量比食入含
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 不确 定性分析 技术 环保 配套 工程 安全评价 中的 应用 研究
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。