基于Cholesky分解法的LHS放射性废物处置场安全不确定性分析.pdf

第 卷第期原子能科学技术V o l ,N o 年月A t o m i cE n e r g yS c i e n c ea n dT e c h n o l o g yA p r 基于C h o l e s k y分解法的L H S放射性废物处置场安全不确定性分析赵润才,玉宇,陈涛(华北电力大学 核科学与工程学院,北京 )摘要:放射性废物处置是一项与国土环境、公众安全、核工业健康以及可持续发展有关的重大问题.安全全过程系统分析是保障放射性废物处置设施从选址、建设、运行到关闭后安全性的重要手段,不确定性分析是其中重要一环.环境变化、人员行为等事前无法控制的外部因素都将对放射性废物处置设施的安全产生重大影响,需要对其进行不确定性评估.在放射性废物处置库的不确定性分析中,参数不确定性分析的计算过程相较于常见的蒙特卡罗模拟的运用场景,其涉及输入的随机参数多、运用模型庞杂,势必需求更小的抽样样本以减少运算时间、提高抽样效率.拉丁超立方抽样(LH S)是不确定性分析中常用的方法,但该方法应用于多维抽样时由于排序质量较低,使得小样本条件下的相关性要求不能得到满足.本文采用C h o l e s k y分解法对拉丁超立方抽样过程中的排列构造过程进行了改进,通过对排列矩阵各行向量进行解耦,以最小化其各维度间的相关性.此改进方案显著降低了拉丁超立方抽样对样本相关性的影响,加速了计算结果的收敛速度.在本文的使用场景下,改进后的抽样方法只需要使用改进前所需样本规模的/,提高了计算效率.关键词:不确定性分析;C h o l e s k y分解法;拉丁超立方抽样;放射性废物处置中图分类号:T L 文献标志码:A文章编号:()收稿日期:;修回日期:通信作者:玉宇d o i:/y z k y o u x i a n C h o l e s k y f a c t o r i z a t i o n b a s e dL a t i nH y p e r c u b eS a m p l i n gf o rS a f e t yU n c e r t a i n t yA n a l y s i so fR a d i o a c t i v eW a s t eD i s p o s a l S i t e sZ HAOR u n c a i,YUY u,CHE NT a o(S c h o o l o fN u c l e a rS c i e n c ea n dE n g i n e e r i n g,N o r t hC h i n aE l e c t r i cP o w e rU n i v e r s i t y,B e i j i n g ,C h i n a)A b s t r a c t:S a f e t y f o c u s e dp r o c e s ss y s t e ma n a l y s i s i s a n i m p o r t a n t a p p r o a c ht oe n s u r e t h es a f e t yo f r a d i o a c t i v ew a s t ed i s p o s a l f a c i l i t i e s t h r o u g h o u t t h e i r l i f e c y c l e,i n c l u d i n gs i t i n g,c o n s t r u c t i o n,o p e r a t i o n,a n dc l o s u r e I nr e c e n ty e a r s,t h eu n c e r t a i n t yi nt h ea n a l y s i sp r o c e s s r e c e i v e s i n c r e a s i n ga t t e n t i o n,a n du n c e r t a i n t ya n a l y s i sa n dm a n a g e m e n tb e c o m ei n t e g r a l p a r t so fs a f e t y f o c u s e dp r o c e s ss y s t e ma n a l y s i sa n ds a f e t ye v a l u a t i o n I n t e r n a t i o n a ld o c u m e n t ss u c ha st h eS a f e t y F o c u s e d P r o c e s sS y s t e m A n a l y s i sa n d S a f e t yE v a l u a t i o nf o rR a d i o a c t i v eW a s t eD i s p o s a lb yI A E Aa n dd o m e s t i cr e g u l a t i o n s l i k et h eS a f e t yR e g u l a t i o n sf o rN e a r S u r f a c eD i s p o s a l o fL o wa n dI n t e r m e d i a t e L e v e lR a d i o a c t i v eS o l i dW a s t ee m p h a s i z e t h e i m p o r t a n c eo fu n c e r t a i n t ya n a l y s i s i nt h e s ep r o c e s s e s P r o p e r l ya d d r e s s i n gu n c e r t a i n t yi ns a f e t y f o c u s e dp r o c e s ss y s t e m a n a l y s i sa n ds a f e t ye v a l u a t i o ni sc r u c i a lf o re n s u r i n gt h ec r e d i b i l i t yo ft h ee v a l u a t i o nr e s u l t s T h e r e f o r e,c o n d u c t i n gu n c e r t a i n t ya n a l y s i s o nc a s e s t u d i e s f o r l o n g t e r ms a f e t ye v a l u a t i o n s i sp a r t i c u l a r l y i m p o r t a n t E x t e r n a l f a c t o r s,s u c ha s e n v i r o n m e n t a l c h a n g e s a n dh u m a nb e h a v i o r,w h i c hc a n n o tb ec o n t r o l l e d i na d v a n c e,s i g n i f i c a n t l y i m p a c t t h es a f e t yo f v a r i o u sn u c l e a rf a c i l i t i e s,i n c l u d i n gr a d i o a c t i v ew a s t ed i s p o s a lf a c i l i t i e s,a n dt h u sr e q u i r eu n c e r t a i n t ya s s e s s m e n t M a n yc o u n t r i e sa n d i n t e r n a t i o n a l o r g a n i z a t i o n s i nt h en u c l e a r i n d u s t r yc o n d u c t e dr e l a t e dr e s e a r c h,f o c u s i n go np r o b a b i l i t y,r i s kg u i d a n c e,a n du n c e r t a i n t ya n a l y s i s T h e yd e v e l o p e ds o m eg e n e r a l e v a l u a t i o nm e t h o d s,s u c ha sM o n t eC a r l os i m u l a t i o na n dL a t i nh y p e r c u b es a m p l i n g,f o ru n c e r t a i n t ya n a l y s i sa c c o r d i n gt os p e c i f i cc i r c u m s t a n c e s I nt h eu n c e r t a i n t ya n a l y s i so f r a d i o a c t i v ew a s t ed i s p o s a l r e p o s i t o r i e s,t h ec o m p u t a t i o n a l p r o c e s so fp a r a m e t e ru n c e r t a i n t ya n a l y s i si sm o r ec o m p l e xc o m p a r e dt ot h ec o mm o na p p l i c a t i o ns c e n a r i o so fM o n t eC a r l os i m u l a t i o n I t i n v o l v e sal a r g en u m b e ro fr a n d o mi n p u tp a r a m e t e r sa n dc o m p l e xm o d e l s,w h i c hn e c e s s i t a t e ss m a l l e rs a m p l es i z e st or e d u c ec o m p u t a t i o nt i m ea n d i m p r o v es a m p l i n ge f f i c i e n c y I nt h i sc a s e,t h es a m p l i n gp r o c e s sw i t hs m a l ls a m p l es i z e sa n dm u l t i p l er a n d o mp a r a m e t e r sw i l l i n e v i t a b l yf a c eac h a l l e n g e:Wh e t h e r t h eu n i f o r m i t ya n dc o r r e l a t i o nb e t w e e nt h ed i m e n s i o n so f t h es a m p l e ds a m p l e sm e e t t h e c o m p u t a t i o n a l r e q u i r e m e n t s T oa d d r e s s t h e i s s u eo f u n i f o r m i t y i nt h es a m p l i n gp r o c e s s,a d o p t i n gL a t i nh y p e r c u b es a m p l i n g(LH S)c o u p l e dw i t h M o n t eC a r l os i m u l a t i o n i sac o mm o nc h o i c e I nt h i s s t u d y,LH Sw a s i m p r o v e db ya p p l y i n gt h eC h o l e s k yd e c o m p o s i t i o nm e t h o di nt h ep a r a m e t e ru n c e r t a i n t ya n a l y s i so fs a f e t y f o c u s e dp r o c e s ss y s t e ma n a l y s i s f o rr a d i o a c t i v ew a s t ed i s p o s a l r e p o s i t o r i e s T h e i m p r o v e dLH Sm e t h o dt h r o u g ht h i sa p p r o a c hc a nu s es m a l l e rs a m p l e st om e e t t h ec o r r e l a t i o nr e q u i r e m e n t so f t h e c o m p u t a t i o n,t h e r e b ye n h a n c i n gs a m p l i n ge f f i c i e n c y I n t h e s c e n a r i oo f t h i ss t u d y,t h e i m p r o v e ds a m p l i n gm e t h o do n l yr e q u i r e so n e t e n t ho f t h es a m p l es i z en e e d e di nt h eo r i g i n a lm e t h o d K e yw o r d s:u n c e r t a i n t ya n a l y s i s;C h o l e s k yf a c t o r i z a t i o nm e t h o d;L a t i nh y p e r c u b es a m p l i n g;r a d i o a c t i v ew a s t ed i s p o s a l安全全过程系统分析是保障放射性废物处置设施从选址到建设、运行乃至关闭后安全性的重要手段,近年来,分析过程中的不确定性逐渐受到关注,不确定性分析与管理已成为安全全过程系统分析和安全评价的重要组成部分,国际原子能机构 放射性废物处置安全全过程系统分析和安全评价、国内 低、中水平放射性固体废物近地表处置安全规定 等文件均对此提出了要求.如何妥善处理安全全过程系统分析和安全评价中的不确定性,将会影响评价结果的可信程度,因此针对长时间尺度的安全评价,开展对案例的不确定性分析研究就显得尤为重要 .环境变化、人员行为等事前无法控制的外部因素都将对包含放射性废物处置设施在内的各类核设施的安全产生重大影响,需要对其进行不确定性评估.世界各国以及核行业相关的国际组织均开展了相关的研究,重点关注概率、风险指引、不确定性分析等内容,并根据实际情况开发了一些通用的评价方法,用于不确定性分析,如蒙特卡罗模拟、拉丁超立方抽样等 .在放射性废物处置库的不确定性分析中,参数不确定性分析的计算过程相较于常见的蒙特卡罗模拟的运用场景,其涉及输入的随机参数多、运用模型庞杂,势必需求更小的抽样样本以减少运算时间、提高抽样效率.在此情况下,小样本多随机参数的抽样过程必然会涉及一个原子能科学技术第 卷问题,即抽样样本各维度间的均匀性与相关性能否满足计算所需.为解决抽样的均匀性问题,采用拉丁超立方抽样(LH S)的蒙特卡罗模拟是常见的选择,本文针对放射性废物处置库安全全过程分析中参数不确定性分析过程,尝试采用C h o l e s k y分解法改进拉丁超立方抽样.经此方法改进的拉丁超立方抽样可以使用更小的样本以满足计算所需的相关性要求,进而提高抽样效率.不确定性分析方法 拉丁超立方抽样与其存在的问题为解决抽样的均匀性问题,采用拉丁超立方抽样的蒙特卡罗模拟是常见的选择.该方法已广泛应用于电网计算、飞行器设计、建筑材料设计等领域,是 世纪 年代末,由M c K a y、B e c k m a n和C o n o v e r提出的一种基于分层抽样的多维抽样方法 .其原理是将样本空间划分成数个独立的子空间,并在每个子空间中随机抽取个样本点.这样,每个子空间都有个样本点,而且样本点之间是独立的.其具有每个维度的单个子空间只抽样次的特点,可减少抽样样本在样本空间的局部聚集,大幅提高抽样的均匀性,进而提高抽样效率.目前存在的主要问题是拉丁超立方抽样没有解决抽样样本的相关性要求.理想的抽样样本各随机变量间应是相互独立的.这是因为研究中涉及的随机参数之间的相关性可以分为两种:相关性已知的参数和相关性未知的参数.而相关性已知的参数可以进行协变量正交化等方式进行处理以达到相互独立,进而被囊括在同一物理模型中,如处置库中核素的库存量;相关性未知的参数则应假设其为相互独立的.对于拉丁超立方抽样,其在多维的情况下,会对各维度分别生成一个随机排列,需要生成随机排列矩阵PK N,其每行的元素值代表采样矩阵SKN对应行元素的位置,按照序列对每层进行抽样获得原始样本RK N,再对各独立维度进行整合(图),即每行的样本按照其在相应行的排列进行映射.图多维情况下拉丁超立方抽样示意图F i g S c h e m a t i cd i a g r a mo fL a t i nh y p e r c u b es a m p l i n gu n d e rm u l t i d i m e n s i o n a l s i t u a t i o n在此过程会产生一些相关性结构,如排列矩阵中两行元素出现相似的排序,使得部分样本空间没有采样等,如图所示,此种情况即可称为排序质量低.这种情况会影响抽样样本各维度间的相关性,使得抽样样本的相关性特征偏离计算要求,导致实际计算结果与理论结果的偏差大于预期.这种排序质量低的情形会随样本数量的增加而减轻,为了减少排序质量低的抽样结果对计算结果的影响,需要更多的抽样次数,即抽样效率会降低.改进的方案本文使用C h o l e s k y分解法协助拉丁超立图多维情况下低排序质量示例F i g E x a m p l eo f l o w o r d e rq u a l i t yi nm u l t i d i m e n s i o n a l s i t u a t i o n第期赵润才等:基于C h o l e s k y分解法的LH S放射性废物处置场安全不确定性分析方法以解决小样本条件下如何满足相关性要求的问题.从本质上,C h o l e s k y分解是一种分解矩阵的方法.C h o l e s k y分解可以把矩阵分解为一个下三角矩阵以及它的共轭转置矩阵的乘积.C h o l e s k y分解要求被分解的矩阵为半正定矩阵.如若被分解的矩阵是正定的,则分解结果是唯一的.C h o l e s k y分解已广泛用于金融工程、人工智能等领域,主要用于生成符合要求的独立随机变量,也用于相关的蒙特卡罗模拟与广义最小二乘回归.本文利用其作为一种降低矩阵各行之间相关性的方法,C h o l e s k y分解的主要思想是:修正排列矩阵PKN,利用C h o l e s k y分解构造近似正交的新排列矩阵GKN,使得各行元素的大小不变的同时,对采样矩阵SKN中各元素的位置进行重新排列,从而降低其行相关性.其本质是将修正前排列矩阵的各行视作具有相关性的另一组随机变量.利用C h o l e s k y分解对PK N进行线性变换,以对其各行向量进行解耦.排列矩阵PKN的相关系数矩阵P是对称且半正定的,可采用C h o l e s k y分解法构造排列矩阵.半正定性的定义为:任意非零向量x,都有xTA x,则矩阵A为半正定矩阵.对任意随机变量矩阵XX,X,Xn,设其均值矩阵为,n,协方差矩阵为Dd,d,dn.如 果 构 造 一 个 矩 阵:XXd,Xd,Xnndn,则 相 关 系 数 矩 阵XTX.故 对 于 任 意 非 零 向 量有:TTXTXXD.C h o l e s k y分解法构造新排列矩阵GK N的步骤如下:)初始化排列矩阵PK N,该矩阵的每行由整数l,N的随机排列组成,如果原先存在一个排列矩阵PK N,此步骤可跳过.)随机排列矩阵PK N两行之间的相关性可用一个K阶的相关系数矩阵 来表示:Pi j,i,K;j,K()i j为对应两行之间的相关系数:i jKk(xi,kxi)(xj,kxj)Kk(xi,kxi)Kk(xj,kxj)()其中:xi,k和xj,k分别表示矩阵第i行和第j行的第k个元素;xi和xj分别表示第i行和第j行的平均值.随机变量的相关系数矩阵必然是一个半正定对称矩阵,因此可以对其进行C h o l e s k y分解,得到一个实数的非奇异下三角方阵D,并满足:PD DT()由于D是非奇异的,非奇异矩阵必是可逆的,可以利用D构造一个列相关性较小的矩阵:GK ND PKN()GKN的行相关性小于PKN,由于GK N的元素并非正整数,不能直接作为排序,可以用其中各元素的大小关系作为新的排列顺序,用其大小关系构成的新排列相关性虽与GKN不能等同,但相较于PK N相关性有显著下降.算例分析图案例景象示意图F i g L l u s t r a t i v ed i a g r a mo f c a s es c e n a r i o 算例的景象设计本文设计的景象示于图.该景象是一个低放处置库的事故景象,假定从处置设施关闭后开始评估.基于 核安全法规技术文件放射性废 物 处 置 安 全 全 过 程 系 统 分 析(NN S A HA J )中对于模型的分类,本模型属于保守模型.具体景象设定如下.)处置库遭遇外部突发事件,如地震、地质运动、人类钻探活动等 ,导致废物包装破原子能科学技术第 卷损、工程屏障失效、地下水侵入处置库并将破损废物包装中的放射性核素直接浸出.为便于研究计算方法,相较实际情况进行了简化,其包含的不合理性不影响计算结果对抽样方法的验证.首先,考虑到 年的核素迁移时长远大于核素从破损废物包装中浸出所需的时间,加之包装破损的形态复杂、破损原因多样,故本文在近场过程中进行了简化,所用模型的浸出过程为直接浸出.其次,考虑到美国的评估经验,放射性核素在包气带中的流动对处置库的最终性能影响较小,因此在本案例中也忽略核素在包气带中的迁移,而是假设整个处置库都位于地下水的饱和带中,放射性核素直接在饱和带中迁移扩散,且在整个迁移过程中,地下水的流速保持不变.)景象所经历的年限设置为处置库正式关闭后的 年.此时处置库处于无人运行的被动状态,且由于此时包装与屏障已破损,之后的人类入侵事件无法进一步影响后果,故不考虑在此时间点后的人类入侵事件.关于年限的确切起点,本研究参照I A E A给出的处置库案例设定时间框架,认为“处置库彻底封存、正式关闭”即为“取消主动控制”.并依照本研究相关项目任务书的规划要求,将经历年限的考察范围增加至 年.I A E A给出的处置库案例的时间框架如表所列.表 I A E A给出的处置库案例的时间框架T a b l eT i m e f r a m e f o rd i s p o s a l r e p o s i t o r yc a s e sp r o v i d e db yI A E A阶段时间/a备注废物处理阶段 厂址内的活动要保证关闭后的风险机构或政府主动控制期 有围栏和巡逻,以防止动物和未经许可的人员进入取消主动控制 被动控制期(地方或国家)覆盖层顶部局部损坏,可能无法修复不受控制(假定资料记录全部丢失)拱顶暴露,近场中可能形成多孔流质,废物包装逐渐暴露并降解)景象考虑了种放射性核素,分别为 I、C s、S r,均为常见的核电厂裂变产物.其中,I是长寿命核素的代表,也是在环境中迁移较快的核素,在低放废物中其活度浓度不超过 MB q/k g.C s和 S r的半衰期都在 年左右,属于中等寿命核素,在环境中迁移较慢,土壤等介质对这两个核素的吸附能力较强.它们在低放废物中的活度浓度不超过G B q/k g.)放射性废物处置场的库存量不定,因此放射性废物的库存量设置为输入参数.)关键人群摄入途径为饮用受污染地下水、接触污染土壤和食用污染农作物种途径.关键人群距处置库 k m.算例公式选取)近场浸出处置库关闭时,存放的放射性废物中有部分包装有破损.由于处置库在关闭时,工程屏障系统失效,地下水进入处置库,这些地下水接触到包装破损的废物,将其中的放射性核素浸出,放射核素在近场浸出的浓度为:CN,iQiRiFW()其中:CN,i为近场被地下水浸出后核素i在地下水中的浓度,B q/m;Qi为核素i在处置库关闭时的存放总量,B q;FW为进入处置库的地下水流量,m/a;Ri为核素i的浸出速率,即单位时间内浸出的放射性核素 占总量的比 例,a;为破损包装的比例.)远场迁移放射性核素在地质介质中的迁移扩散称为远场迁移.远场迁移的机理较明确,包括地下水的对流输运,污染物的机械弥散、分子扩散,地质介质对放射性核素的吸附、阻滞作用,放射性衰变和其他衰减过程,以及地下水的源、汇项对迁移的影响.本文主要讨论地下水的对流输运、核素的机械弥散、介质对核素的吸附,以及放射性衰变等过程.第期赵润才等:基于C h o l e s k y分解法的LH S放射性废物处置场安全不确定性分析一般用对流扩散方程来描述污染物在地下水中的迁移行为,该方程的形式如下:Rd,iCF,itDL,iCF,ixuCF,ixiRd,iCF,iI()Rd,iKd,i()其中:Rd,i为核素i的阻滞系数,表征核素在固相上的吸附对其迁移行为的影响;Kd,i为核素i在固相介质上的分配系数,m/k g;为固相介质的密度,k g/m;为固相介质的孔隙率;CF,i为核素i在远场地下水中的浓度,B q/m;t为核素迁移时间,a;x为在地下水的流动方向下游距处置库的直线距离,m;u为地下水的流速,m/a,其大小与流入处置库的地下水量有关,计算公式为uFW/A,A为处置库处置单元的横截面积(m);i为核素i的衰变常量,a;I为单位时间单位体积多孔介质中所得到的污染物质量,B q/(ma),反映的是各种源、汇项的作用,若不考虑各种源、汇项的作用,则I;DL,i为核素i在地质介质中的弥散系数,m/a,污染物的水动力弥散包括机械弥散和分子扩散.机械弥散是由于多孔介质空隙和固体骨架的存在而造成流体的微观速度在空隙中的分布无论是大小还是方向都不均匀的现象.污染物在介质中的机械弥散能力不仅与介质的性质有关,而且与地下水的流动速度有关.当地下水的流速为m/a时,不存在机械弥散,则DL,i为核素i在地质介质中的扩散系数.分子扩散是指在浓度梯度作用下,污染物从高浓度向低浓度位置的扩散.介质中原先不存在 I、C s、S r,从近场浸出的核素浓度CN恒定,地下水的流速恒定,介质是均匀和各向同性的,迁移发生在地下水流的下游,则上述对流扩散方程存在如下解析解:CF,iCN,ie x p(uw)xDL,ie r f c(Rd,ixwtRd,iDL,it)e x p(uw)xDL,ie r f c(Rd,ixwtRd,iDL,it)()wuiRd,iDL,i()剂量转换由于处置库中放射性核素的迁移,作为饮用水的井水可能会受到污染,直接饮用被污染的地下水会导致人体受到内照射.受地下水污染的影响,周边的土壤中也会出现从处置库中释放出的核素,这些土壤中的核素会对人体产生外照射.农作物从土壤中吸收并累积放射性核素,食用这些农作物人体会受到内照射.根据本文的景象设计,从处置库中释放出的放射性核素对于关键人群的照射途径有个:一是受污染土壤的直接外照射;二是直接饮用受污染的地下水的内照射;三是食用污染土地上生长的作物所导致的内照射.这个途径的剂量计算方法如下.外照射剂量DE,i(S v/a)的计算公式如下:DE,iCg,iL DF s,iOf()Cg,iKd,iCF,i()式中:Cg,i为核素i在土壤中的浓度,B q/k g,主要从污染地下水中通过吸附作用转移而来;L为对外照射剂量有贡献的表层土壤的厚度,m;DF s,i为核素i的外照射剂量的转换系数,(S v/a)/(B q/m);Of为关键人群受土壤外照射的时间份额.直接饮用受污染地下水所导致的内照射剂量DW,i的计算公式如下:DW,iCF,iHWDFW,i()式中:HW为关键人 群每年的饮 水量,m/a;DF W,i为核素i的内照射剂量的转换系数,S v/B q.食用污染土地上生长的作物导致的内照射剂量DV,i的计算公式如下:DV,iCV,iHVDFW,i()CV,iFV,iCg,i()式中:CV,i为所食用作物中放射性核素i的浓度,B q/k g(鲜重),其主要来源于作物从土壤中吸收、富集的放射性核素;FV,i为放射性核素i从土壤中转移到作物中的转移系数;HV为关键人群每年食用作物的量,k g/a.总剂量Dt o t(S v/a)的计算公式为:Dt o ti(DE,iDW,iDV,i)()算例参数选取)输入参数根据本案例的景象设置,算例模型分为个模块,分别是近场浸出模块、远场迁移模块和原子能科学技术第 卷剂量转换模块.根据模型需求,本文输入参数的取值范围和概率密度分布如表 所列.)输出参数本文算例共 个输出参数,其中远场地下水核素总浓度(CF)和总有效剂量(Dt o t)是两个最重要的输出参数,其余输出参数分别为各核素在远场地下水中的浓度与按照核素与途径进行区分的有效吸收剂量.表输入参数的取值范围和概率密度分布T a b l eR a n g eo fv a l u ea n dp r o b a b i l i t yd e n s i t yd i s t r i b u t i o nf o r i n p u tp a r a m e t e r参数核素取值范围概率分布Qi,B q I 均匀分布 C s S r FW,m/a 对数正态分布 均匀分布Ri,a I 对数均匀分布 C s S r DL,i,m/a I 对数均匀分布 C s S r ,k g/m 常数 对数正态分布Kd,i,m/k g I 对数均匀分布 C s S r i,a I 常数 C s S r A,m 常数L,m常数DF s,i,(S v/a)/(B q/m)I 常数 C s S r Of 常数HW,m/a 常数DFW,i,S v/B q I 常数 C s S r FV,i I 对数均匀分布 C s S r HV,k g/a 常数第期赵润才等:基于C h o l e s k y分解法的LH S放射性废物处置场安全不确定性分析按 低中水平放射性固体废物的浅地层处置规定(G B )与 生活饮用水卫生标准(G B )的要求,输出参数CF和Dt o t可以作为安全指标.结果分析 案例结果为验证本文所提出的方法,以拉丁超立方 次抽样结果作为标准,依照算例进行了改进前后 次与 次抽样计算.使用拉丁超立方法计算所得各主要参数的均值与标准差如表、所列.由表、可见,关键人群的总吸收剂量低于 低中水平放射性固体废物的浅 地 层 处 置 规 定(G B )中 的 m S v/a剂量限值,地下水核素浓度也低于 生活饮用水卫生标准(G B )中 B q/L的总放射性指标指导值.表、的绝大多数模拟结果中,CF和Dt o t的数值趋近于是正常的,处置场在设计上应具有较高的保守性.即使景象设定较为严苛,绝大多数情况下也应是安全的,一旦出现极小概率事件才会有一定的危害性,在数据上的表现就是计算结果的变异系数极高,均值较标准差小个数量级左右.依照表、中的计算结果,并不能判断在本文设定的情境下,处置库的安全水平是否符合标准,引入C h o l e s k y分解法的改进效果也不够直观明显,因此对计算所得结果进一步处理,生成互补累积分布函数(C C D F),以真实判断安全水平与改进结果.次抽样条件所生成的C C D F曲线如图、所示.表 次抽样下改进前后的计算结果T a b l eC a l c u l a t i o nr e s u l t sb e f o r ea n da f t e r i m p r o v e m e n tb a s e do n s a m p l i n g统计量LH S法C h o l e s k y法均值标准差均值标准差Dt o t,m S v/a DI,m S v/a DC s,m S v/a DS r,m S v/a CF t o t,B q/L CF I,B q/L CF C s,B q/L CF S r,B q/L 表 次抽样下改进前后的计算结果T a b l eC a l c u l a t i o nr e s u l t sb e f o r ea n da f t e r i m p r o v e m e n tb a s e do n s a m p l i n g统计量LH S法C h o l e s k y法均值标准差均值标准差Dt o t,m S v/a DI,m S v/a DC s,m S v/a DS r,m S v/a CF,B q/L CF,I,B q/L CF,C s,B q/L CF,S r,B q/L 对于远场地下水中核素浓度CF,以保守性作为选取标准.以 生活饮用水卫生标准 中对饮用水的要求来约束应足够保守,其对饮用水图远场地下水中核素浓度的C C D F曲线F i g C C D Fc u r v e f o r r a d i o n u c l i d ec o n c e n t r a t i o n i ng r o u n d w a t e r i nf a r f i e l d原子能科学技术第 卷图总吸收剂量的C C D F曲线F i g C C D Fc u r v e f o r t o t a l a b s o r b e dd o s e的要求为总放射性指标指导值为 B q/L,总放射性指标指导值为B q/L,换算后约为 B q/m和 B q/m.在以保守性为前提下,以 B q/m为阈值,从图 b可见,大于该阈值的概率约为 ,即P(CF B q/m).同样从图 b可见,在 置信水平下总浓度大致小于 B q/m.对于吸收剂量Dt o t,根据 低中水平放射性固体废 物 的 浅 地 层 处 置 规 定,其 阈 值 设 为 m S v/a.从图可看出,大于该阈值的概率为;在 置信水平下,换算后总剂量大致小于 m S v/a.改进结果r m s可以改进过程对随机变量的相关性影响,较小的相关性代表着更少的相关性结构.故本研究将其引入以评价改进结果.整体的相关性可用相关系数均方根r m s表示:r m sKiKji jKK(K)()为了对比改进前后的分析结果,除相关性均方根r m s外,再引入以拉丁超立方抽样 次为基准的误差度指标:a c ca c c()其中:为用于对比的计算结果均值;a c c为作为基准的拉丁超立方抽样 次的计算结果.改进前后每个输出参数均可得出一个误差度指标.改进前后的相关系数均方根r m s及误差度指标列于表.表改进前后的衡量指标T a b l e I m p r o v e m e n t e v a l u a t i o nm e t r i c sb e f o r ea n da f t e r i m p r o v e m e n t衡量指标L H SC h o l e s k y法 次 次 次 次r m s Dt o t D