温度对75Se（Ⅳ）在我国北山花岗岩中扩散影响的再思考.pdf

本文对S e在我国高放废物地质处置库预选围岩北山花岗岩中的扩散研究进行了总结,包括研究方法、获取的有关数据以及尚待开展的工作等,重点分析了温度对 S e()在北山花岗岩中扩散的影响及可能的影响机制,旨在为国内同行提供有益参考.关键词:扩散;S e;北山花岗岩;温度中图分类号:T L 文献标志码:A文章编号:()d o i:/h h x Y X R e t h i n k i n go nI n f l u e n c eo fT e m p e r a t u r eo nD i f f u s i o no f S e()i nB e i s h a nG r a n i t eWANGC h u n l i,YANGX i a o y u,X I AND o n g f a n,Z HUZ h a o w e n,N u c l e a ra n dR a d i a t i o nS a f e t yC e n t e r,B e i j i n g ,C h i n a;B e i j i n gN a t i o n a lR e s e a r c hC e n t e r f o rM o l e c u l a rS c i e n c e,R a d i o c h e m i s t r y&R a d i a t i o nC h e m i s t r yK e yL a b o r a t o r yf o rF u n d a m e n t a lS c i e n c e,C o l l e g eo fC h e m i s t r ya n dM o l e c u l a rE n g i n e e r i n g,P e k i n gU n i v e r s i t y,B e i j i n g ,C h i n aA b s t r a c t:T h er e s e a r c h e so nt h ed i f f u s i o no fk e yn u c l i d eS e i nB e i s h a ng r a n i t e,w h i c hi s t h ep r e s e l e c t e ds u r r o u n d i n gr o c ko fC h i n ash i g h l e v e l r a d i o a c t i v ew a s t eg e o l o g i c a l r e p o s i t o r y,w e r es u mm a r i z e d,i n c l u d i n gr e s e a r c hm e t h o d s,r e l e v a n t d a t ao b t a i n e da n dw o r kt ob e c a r r i e do u t M o r e o v e r,t h e i n f l u e n c eo f t e m p e r a t u r eo nt h ed i f f u s i o no f S e()i nB e i s h a ng r a n i t ea n dt h ep o s s i b l e i n f l u e n c em e c h a n i s mw e r ea n a l y z e de m p h a t i c a l l y,i no r d e r t op r o v i d eu s e f u lr e f e r e n c e s f o rd o m e s t i cp e e r s K e yw o r d s:d i f f u s i o n;S e;B e i s h a ng r a n i t e;t e m p e r a t u r e当前,基于多重屏障理念的深地质处置是国际公认的高水平放射性废物(以下简称高放废物)安全处置的最有效方式.多重屏障系统通常由工程屏障和天然地质屏障组成,以将废物与生物圈进行有效隔离.为确保高放废物能在其放射性衰减到自然水平所需的长时间范围内得到有效安全隔离,处置库建造前的安全预评价至关重要.而作为高放废物处置库系统中的最后一道屏障,围岩(天然屏障)可有效阻滞地下水流中的放射性核素.因此,研究关键放射性核素在处置库围岩中的迁移行为是处置库安全预评价的一个重要方面.国际上许多国家将花岗岩等结晶岩作为深地质处置库的潜在围岩.位于甘肃省北山地区的北山花岗岩被认为是中国潜在深地质处置库的候选围岩.在高放废物地质处置的安全评价中,S e(半衰期为 a)的阴离子被认为是为数不多的对未来数十万年内释放到生物圈的最终剂量有显著贡献的长半衰期放射性核素之一,因此 S e的迁移行为引起了研究者的广泛关注.作为一种氧化还原敏感元素,S e在环境中的溶解度和流动性在很大程度上取决于其氧化状态和化学形态,而其氧化状态和化学形态受周围氧化还原电位(E h)p H条 件 控 制.以 阴 离 子 形 式 存 在 的S e()和S e()更易溶解和迁移,而低氧化态的S e(,及)通常以较难溶解的固体形式存在.为获取可靠的参数并建立有效模型,以对中国的潜在深地质处置库进行安全预评价,研究溶液状态下 S e的优势种态在北山花岗岩中的迁移行为是非常必要的.由于大量高释热率放射性核素的存在,高放废物固化体将产生大量热量并引起温度升高,而处置库是一个相对密闭的体系,散热困难,因此在处置的前一百年期间,废物罐甚至处置库近场的温度会达到 甚至更高.增加的热负荷不仅会影响容器的化学热水力学行为,还会影响附近围岩的物理和化学性质 .国际上关于温度对核素在处置库预选围岩中迁移影响的研究相对较少,仅有少量研究集中在部分核素在膨润土及黏土岩中的扩散 ,而在结晶类岩石如花岗岩中的扩散研究则少之又少.北京大学的陈涛 及黎春 等分别研究了温度()对I和T c O在 花 岗 岩 中 扩 散 的 影 响,并 用 修 正 的N e r n s t方程对实验结果进行了一定的解释.为研究 更 接 近 处 置 库 实 际 条 件 的 高 温(如 接 近)对放射性核素迁移扩散的影响,北京大学刘春立教授课题组 专门自主研发了高温扩散装置,通过研究 S e()在我国北山花岗岩中的扩散行为,证实了自主研发高温扩散装置的可靠性,并获取了 下S e()在我国北山花岗岩中的扩散数据.本文通过系统回顾不同温度下S e在我国北山花岗岩中的扩散行为,结合有关领域的最新研究进展,探 讨 了 不 同 温 度 尤 其 是 高 温 条 件 下 S e()在北山花岗岩中的扩散和迁移机制,并对后续工作提出了建议.高放废物地质处置库安全评价基础数据 扩散系数的测定 关键放射性核素(R N)在我国预选处置库围岩 北山花岗岩中的扩散系数是对潜在高放废物地质处置库开展安全评价所需的关键基础数据.当前,国内外可测定核素在介质中的扩散系数的方法主要有:通透扩散法(T h r o u g h d i f f u s i o n)、内扩散法(I n d i f f u s i o n)、外扩散法(O u t d i f f u s i o n)、柱法实验(C o l u m nt e s t)等,其中,利用小而薄的样品岩片进行的T h r o u g h d i f f u s i o n和I n d i f f u s i o n法是研究基质扩散的经典实验室方法 .T h r o u g h d i f f u s i o n(通透扩散法)T h r o u g h d i f f u s i o n法是目标核素在浓差作用下,在扩散装置中从高浓度一侧穿透扩散介质进入到低浓度一侧的实验方法,是进行核素在环图P S I的T h r o u g h d i f f u s i o n法实验装置图 F i g E x p e r i m e n t a l s e t u pf o rT h r o u g h d i f f u s i o na p p a r a t u so fP S I 境介质中扩散行为研究的经典方法之一.图所示为一种典型的T h r o u g h d i f f u s i o n实验装置,为瑞士保罗谢勒研究所(P a u lS c h e r r e rI n s t i t u t e,简称P S I)设计.此扩散装置由两个储液池和一个扩散池组成,目标扩散介质固定于扩散池中间,两个储液池分别为高核素浓度的源液池和低核素浓度的取样池,借助蠕动泵的作用使储液池中的溶液进入扩散池介质两侧并实现整个系统的持续循环,如此,目标核素即从高浓度一侧向低浓度一侧进行扩散.借助T h r o u g h d i f f u s i o n法,可计算出核素在目标介质中的有效扩散系数(De)及吸附分配系数(Kd)等关键参数,而且可以通过改变溶液/介质体系的实验环境来探讨第期王春丽等:温度对 S e()在我国北山花岗岩中扩散影响的再思考相关因素如温度、p H、离子强度等对核素在目标介质中扩散和吸附行为的影响.I n/O u t d i f f u s i o n(内/外扩散法)I n d i f f u s i o n法一般将扩散介质制成柱状或饼状,使核素溶液接触扩散介质的一侧,由溶液向介质内扩散,并通过测量目标核素在扩散介质中扩散但未穿透时在介质中的浓度分布来对核素在介质中的扩散行为进行研究.O u t d i f f u s i o n与I n d i f f u s i o n法相对应,是将放射源埋在扩散介质中间,使目标核素在扩散介质中从内向外扩散,进而测定核素在介质中扩散但未穿透时其在介质中的浓度分布.相对而言,I n d i f f u s i o n法是更为常用的过渡态实验方法.毛细管法是一种典型的I n d i f f u s i o n法,其利用小巧的毛细管作为实验工具从而实现了传统I n d i f f u s i o n法中扩散装置的小型化及便捷化.王祥科等 首先将毛细管法应用到核素迁移的研究中,其方法如图所示:将岩石介质粉末填入毛细管,进行压实处理后,将毛细管浸泡在装有核素的溶液中,一般一端封口,另一端敞开,使核素在介质中进行扩散,扩散完成后将毛细管进行定长切割并测量放射性计数分布,从而推算出放射性核素在对应介质中的扩散系数及吸附分配系数等参数.I n d i f f u s i o n法只能得到表观扩散系数Da,无法直接获取有效扩散系数De,而Da是一种全局参数,在反映扩散过程的内在机理方面不如De.图毛细管法研究核素迁移的示意图 F i g D i a g r a mf o r t h es t u d yo f r a d i o n u c l i d em i g r a t i o nb yc a p i l l a r ym e t h o d C o l u m nt e s t(柱法实验)柱法实验是通过将目标扩散介质填充成柱后,从扩散柱的一侧以恒定流速注入核素溶液,而在扩散柱的另一侧检测流出液中核素的浓度,得到目标核素的流出曲线,进而拟合出核素在目标扩散介质中的弥散系数Dd及吸附分配系数Kd等关键参数.柱法实验可通过改变注入溶液的流速及其它条件来考察水动力学及相关因素的影响.如图所示,为法国B a z e r B a c h i等 设计的一种柱法实验装置.高温条件下扩散系数的测定 高温扩散实验装置 利用上述实验方法制作的常规扩散实验装置,一般通过调整实验装置所处环境的温度来实现对实验温度的控制.然而,受实验装置材料及工艺的限制,实验温度仅能控制在一定范围内(一般低于),无法满足开展高温扩散实验的需求.对于更高温度条件下的扩散实验,北京大学刘春立教授课题组 专门自主研发了高温扩散装置(发明专利号C N B)且首次成功在国内使用.该高温扩散装置依据通透扩散法理念制作,实验运行体系借鉴了国际上通用的瑞士P S I的扩散实验 设计,其扩散 实验原理图 参见图柱法实验研究核素迁移行为的示意图 F i g E x p e r i m e n t a l s e t u pf o r c o l u m nt e s t 图,实物运行图如图所示.核化学与放射化学第 卷图高温扩散装置实验运行实物图 F i g P i c t u r e so f t h ep r e l i m i n a r yh i g h t e m p e r a t u r ed i f f u s i o ns e t u p 该高温扩散装置的制作及组装工艺为:)采用具有良好化学稳定性的聚四氟乙烯(P T F E)制成内胆,外加具有固定与压紧双重功能的金属外框,扩散内胆安置于所述金属外框开口处,并由金属外框的顶板顶紧以固定在外框内,实物如图所示;)其扩散池内胆如图所示,实际使用中将中心扩散岩片两端分别用O型垫片夹住后放入扩散池体和相应卡槽内,并施力将扩散池体与合为一体并压紧,使得岩片两侧形成两个扩散池.扩散接口,扩散池体,O型垫圈,中心扩散岩片,扩散池体图扩散池内胆装配图及零件分解图 F i g G e n e r a l a n dp a r t i a l a s s e m b l yd r a w i n g so f t h e i n n e rd i f f u s i o nc o n t a i n e r 扩散池体组装过程中以岩片为分界形成两侧扩散池体,分别通过扩散池内胆上的模块化接口与设定体积的储液池进行连接,从而形成外接源液池和取样池,并借助蠕动泵的循环使两个扩散池腔体充满溶液,形成完整的源液池及取样池循环系统.将高温扩散装置的主体扩散池放入烘箱中,通过烘箱进行温度调控从而实现对设定高温的控制.借助高温扩散装置,可获取核素在目标扩散介质中的有效扩散系数De及吸附分配系数Kd等关键参数.不同温度条件下S e在我国北山花岗岩中的扩散核工业北京地质研究院王驹在我国高放废物地质处置库选址及场址评价等方面做了大量工作,为国内同行提供了宝贵的花岗岩及地下水样品.在此基础上,国内许多研究院所和高校 开 展 了 核 素 在 北 山 花 岗 岩 上 的 迁 移 研究.不同温度条件下 S e在我国北山花岗岩中的扩散数据结果列于表.实际研究工作中,考虑 S e半衰期太长、不易获得,且后期废物处理处置问题突出,而对于相对较重且质量数差异较小的同位素而言,其化学性质基本一致,因此一般采用半衰期较短的 S e(T/d)作为替代核素进行实验.由表可见,当前国内对 S e在北山花岗岩中扩散研究所采用的实验方法主要为通透扩散法及毛细管法,而通透扩散法作为核素在环境介质中扩散行为研究的经典方法,且实验中采用更接近处置库条件下的北山花岗岩岩片进行实验,因而应用更为广泛.从所采用核素 S e的种态来看,不同实验中普遍选择S e()作为研究对象,主要是因为在S e的所有种态中,S e()和S e()均具有高溶解性且更易迁移,而较之于S e(),S e()是高放废物处置深地质条件下S e的主要存在种态,因此选用S e()进行实验得到的De可为后续安全评价提供更为准确和安全的数据输入.表汇总的扩散系数结果表明:当采用通透扩散法时,不同实验条件下 S e()在北山花岗岩中的有效扩散系数De的数值范围为 m/s,大多在 m/s量级范围,受温度及p H影响较大;当采用毛细管法时,不同p H条件下 S e()在北山花岗岩中的表观扩散系数Da的数值范围为 m/s,大多在 m/s量级范围,受p H影响较大.同时,毛细管法测定的扩散系数普遍比近似条件下通透扩散法得到的结果第期王春丽等:温度对 S e()在我国北山花岗岩中扩散影响的再思考表不同温度条件下S e在北山花岗岩中的扩散系数T a b l eD i f f u s i o nc o e f f i c i e n t so fS eo nB e i s h a ng r a n i t eu n d e rd i f f e r e n t t e m p e r a t u r e sS e的种态/价态S e的浓度实验方法实验条件温度/实验体系(背景溶液)离子强度I/(m o lL)p H扩散实验总时间/d扩散系数/(ms)文献 S e()示踪量 S e()(源液活度浓度:B q/m L)T D 北山花岗岩岩片(B S 钻孔,m)/北山地下水 De:De:De:De:De:S e()示踪量 S e():m o l/LT D(常规扩散装置)T D(高温扩散装置)北山花岗岩岩片(B S 钻孔,m)/N a C l O水溶液 De:De:De:De:De:De:De:S e()示踪量 S e():m o l/LT D 北山花岗岩岩片(B S 钻孔,m)/N a C l O水溶液 De:De:De:De:De:S e()N aS e O(S e示踪):m o l/L毛细管法(I D)破碎北山花岗岩(B S 钻孔,m)/N a C l O水溶液 Da:Da:Da:Da:Da:Da:Da:Da:S e()示踪量 S e():m o l/LT D 北山花岗岩岩片(B S /钻孔,m)/N a C l O水溶液 De:()S e O示踪量 S e O(源液活度计数:/(sm L)T D 花岗岩岩片(取自北京郊区)/花岗岩的平衡水溶液 Di:S e O示踪量 S e O(源液活度计数:/(sm L)T D 花岗岩岩片(取自北京郊区)/花岗岩的平衡水溶液 Di:S e O示踪量 S e O(源液活度浓度:B q/m L)T D 花岗岩岩片/花岗岩的平衡水溶液 Di:注:)T D:通透扩散法(T h r o u g h d i f f u s i o n);I D:内扩散法(I n d i f f u s i o n);Da:表观扩散系数(a p p a r e n td i f f u s i o nc o e f f i c i e n t,m/s);De:有效扩散系数(e f f e c t i v ed i f f u s i o nc o e f f i c i e n t,m/s);Di:固有扩散系数(i n t r i n s i cd i f f u s i o nc o e f f i c i e n t,在数值上与De相等);)通透扩散法中,为采用常规扩散装置()获取的数值结果,为采用高温扩散装置()获取的数值结果核化学与放射化学第 卷更大(高出近一个数量级),可能是因为毛细管法采用破碎的北山花岗岩作为扩散介质、而通透扩散法采用的是完整北山花岗岩岩片、毛细管内破碎花岗岩的孔隙率要远大于花岗岩片的孔隙率所致.整体而言,较之于毛细管法,通透扩散法获取的扩散系数的直接数据结果更能反映处置库实际条件下的核素迁移特性.需要说明的是,在文献 的工作中,作者分别采用常规扩散实验装置及高温扩散实验装置同步开展了 条件下的扩散实验,且两套实验装置各设置了平行样.实验结果表明,两套实验装置在运行状态均良好的前提下,获取的扩散数据结果高度一致,数据相互印证,证实了北京大学刘春立教授课题组在国内首次使用的高温扩散实验装置获取的实验数据是可靠的.基于温度对 S e()在北山花岗岩中扩散影响的数据结果,总体来看,S e()在北山花岗岩中的扩散随温度的升高而增强,在 范围内,整体受温度影响不显著,De在数值上基本维持在同一数量级;而当温度增大到 时,De有明显的增大,甚至高出一个数量级.因温度高于 的扩散系数仅有文献 报道,在国内外其他文献中未查询到相关数据,考虑到实验误差或实验方法的严谨性等问题,此实验结果是否广泛适用于不同类型的花岗岩尚需更多实验数据证实.温度对 S e()在北山花岗岩中扩散的影响由表可见,有关温度对 S e()在我国北山花岗岩中扩散的影响,实验温度主要在 ,除文献,的工作系统研究了温度对 S e()在北山花岗岩中扩散的影响外,其他工作的实验温度均集中在室温条件,且实验条件各异.为深入探讨不同温度下 S e()在北山花岗岩中的扩散机制,此处以文献 中的数据为基础,系统分析温度对 S e()在北山花岗岩中扩散的可能影响机理,评估后续开展进一步工作的可能方向,以期为深地质处置库的安全评价提供更有力支撑.温度对 S e()在北山花岗岩中扩散的影响将文献 中不同温度条件下 S e()在北山花岗岩中的有效扩散系数De作图,如图所示.整体来看,S e()在北山花岗岩中的De随温度的升高而增大.值得注意的是,在 范围内,De基本维持在同一量级,变化不显著;而当温度增大到 时,De明显增大(高出近一个数量级).根据伍涛等 的报道,温度对分子扩散的影响可用A r r h e n i u s经验公式进行描述:DeAeEa/R T式中:A为指前因子,m/s;R为摩尔气体常数,J/(m o lK);T为热力学温度,K;Ea为表观活化能,J/m o l.Ea的数值大小与物质在介质中的扩散机理有关.分别用A r r h e n i u s公式对不同温度范围内 S e()在北山花岗岩中的De数据点进行拟合,如图所示.可见,尽管不同温度范围内的De数值均可用A r r h e n i u s公式得到较好的拟合结果,但 的点相对于 范围内的拟合曲线有明显的偏移,而且用A r r h e n i u s公式拟合出的不同温度范围内活化能的数值结果也存在较大差别,如表所示,进一步证实了 得到的数据点的特殊性.图 S e()在北山花岗岩中的De随温度的变化曲线 F i g Dev a l u e so f S e()i nB e i s h a ng r a n i t ea sa f u n c t i o no f t e m p e r a t u r e 第期王春丽等:温度对 S e()在我国北山花岗岩中扩散影响的再思考 拟合,拟合图不同温度下De的A r r h e n i u s公式拟合结果 F i g R e s u l t so fDea td i f f e r e n t t e m p e r a t u r e s f i t t e db yt h eA r r h e n i u se q u a t i o n 表A r r h e n i u s公式拟合DeT曲线获取的数据结果 T a b l eP a r a m e t e r so fDea td i f f e r e n t t e m p e r a t u r e sf i t t e db yA r r h e n i u se q u a t i o n 温度/Ea/(k Jm o l)相关系数r 温度对S e()在北山花岗岩中扩散影响的机理探讨鉴于文献 的工作中扩散体系所处溶液的化学环境保持不变,为探究温度对 S e()在北山花岗岩中扩散的影响机制,重点对实验条件下的溶液及分子状态、S e的种态及北山花岗岩本身的结构情况开展进一步分析.)温度变化对溶液及分子状态的影响一方面,热力学温度的升高,会降低水的黏度及密度,增强S e()在水中的溶解性能;另一方面,温度升高将导致整个体系中S e分子的热运动加快,促进S e在花岗岩中的扩散.因此,温度的升高将对S e在北山花岗岩中的扩散有一个整体的促进作用.)实验温度范围内S e的种态分析鉴于实验过程中监测到的溶液体系的p H和E h值无显著变化,基于p H及E h监测结果,利用CHEM S P E C软件 计算不同实验温度下溶液体系中S e的种态分布情况(图).图结果显示,不同温度下溶液中S e的优势种态始终为H S e O和S e O,温度的改变并未对S e的存在种态造成显著影响.故温度对S e()在北山花岗岩中扩散的影响并非由S e的种态变化引起.图不同温度下S e在扩散溶液中的种态分布情况(p H)F i g S p e c i a t i o nd i s t r i b u t i o no fS e i nd i f f u s i o ns o l u t i o na td i f f e r e n t t e m p e r a t u r e s(p H)实验温度范围内北山花岗岩的内部结构分析文献 采用热重分析仪在 下对花岗岩进行热重分析(T GA)时未观测到可识别的质量损失,故此温度范围内,花岗岩应该不存在矿物分解作用.赵星光等 的研究结果也表明,在 的升温范围内,花岗岩的矿物组成未发生改变.通常,溶质在多孔介质中的扩散受介质的孔隙率()、孔隙的曲折度()和收缩度()以及溶质在水溶液中的分子扩散系数(Dm)影响.溶质在多孔介质中扩散的有效扩散系数De及其在水溶液中的分子扩散系数Dm间存在如下关系(Dp为孔隙扩散系数):De DpDmFfDm()式中,/作为一个整体被称为形成因子(f o r m a t i o nf a c t o r,Ff).形成因子Ff 综合考虑了多孔介质的曲折度、收缩度、孔隙率等参数,在一定程度上反映了多孔介质的整体、特定结构特征.因此,形成因子只依赖于孔隙介质的内部结构,可在整体或宏观上概括性地反映岩石本身的性能,一般同一种岩石的形成因子应该近似相同.鉴于北山花岗岩的形成因子(Ff)可用来验证花岗岩的内部结构是否发生改变,文献 中也对实验温度范围内花岗岩的形成因子进行了分析.Ff的数值结果如表所示,Ff随温度变化的示意图如图所示.从北山花岗岩在不同温度下的Ff的数值结果来看,当实验考察温度为 及以下时,Ff虽核化学与放射化学第 卷有波动但变化不显著;而当温度达 时,Ff有较显著的变化.表不同温度下北山花岗岩的形成因子 T a b l eF o r m a t i o nf a c t o ro fB e i s h a ng r a n i t ea td i f f e r e n t t e m p e r a t u r e s 温度/Ffl gFf()()()()()()()注:为采用常规扩散装置()获取的数值结果,为采用高温扩散装置()获取的数值结果图不同温度下北山花岗岩的形成因子(Ff)F i g F o r m a t i o nf a c t o ro fB e i s h a ng r a n i t ea td i f f e r e n t t e m p e r a t u r e s 北山花岗岩的Ff值与文献结果的比较列于表.条件下文献 获取的北山花岗岩的l gFf与瑞典S K B等 报道的结果基本一致,并且与之前同课题组陈涛等 使用的相同花岗岩相比,l gFf值几乎无差异,这也验证了实验数据的可靠性.而当温度达 时,l gFf则偏差较大,此前的文献报道中也未见与此结果接近的数值.推测在温度相对较高的 时,花岗岩在扩散实验的研究体系下其内部结构如曲折度、收缩度、有效孔隙率等发生了部分改变.花岗岩是由不同矿物(长石、石英、云母等)颗粒组成的非均质岩体,不同矿物颗粒具有不同的热膨胀性质和热膨胀各向异性.贺琦 的研究表明,随温度升高,花岗岩内部不同矿物的不均匀热膨胀导致颗粒之间的应力集中,造成原始裂隙压缩、闭合,孔隙度降低;随温度进一步升高,不同矿物间热膨胀差异性进一步体现而发生热破裂,导致数量更多、尺度更大的热致裂隙形成.表文献 获取的北山花岗岩的Ff与文献值的比较T a b l eFfo f S e Oi nB e i s h a ng r a n i t ei nc o m p a r i s o nw i t hl i t e r a t u r ed a t a文献花岗岩类型l gFfS K B F o r s m a r k花岗闪长岩 结晶岩 L a x e m a r花岗岩 瑞典花岗岩 瑞典结晶岩 S i m p e r v a r p花岗岩 黎春等 北山花岗岩 陈涛等 北山花岗岩 文献 ()北山花岗岩 文献 ()北山花岗岩 对于非均质岩体,只要温度改变就会在岩体内部产生热应力,当热应力超过岩体内部晶粒/矿物与胶结物的承载力/黏结强度时,岩石内部就会产生热破裂 .热破裂的实质是裂纹扩展过程中热能的积累和释放过程.岩石热破裂使得岩体裂隙进一步发育,形成更好的孔隙裂隙通道.赵阳升等 测得采自山东平邑的鲁灰花岗岩的热破裂临界温度约为 .W a n g等 的研究表明,美国W e s t e r l y花岗岩的热破裂临界温度为;陈颙等 发现,W e s t e r l y花岗岩的热破裂临界温度为 .W a n g等 的研究表明,花岗岩的热致裂隙形成于 左右.贺琦等 证实,中等温度()荷载范围内北山花岗岩产生的裂隙数量和尺度随温度升高而增大,温度改变了花岗岩内部结构,使其更易发生破裂或发育大尺度裂纹.左建平等 针对性地研究了我国北山花岗岩(取自B S 号钻孔,埋深 m)的热破裂情况:通过实时观察 范围内北山花岗岩的热开裂过程证实,北山花岗岩热破裂临界温度为 (可能存在 的波动范围).当温度超过 时,北山花岗岩的原生裂隙会进一步发育和扩展,并诱发一些新的热裂纹;且随温度升高,热裂纹不仅数量不断增加,而且其长度和宽度均相应有所增加.在 时热破裂裂纹的数量、长度和宽度均有较大程度的增长,是热破裂裂纹快速扩展的一个阶段.需要说明的是,由于各矿物颗粒的热膨胀系第期王春丽等:温度对 S e()在我国北山花岗岩中扩散影响的再思考数差异及各向异性,以及不同矿物的沉积条件均有所不同,因此每个实验样品的热破裂临界温度基本均有差异,热破裂临界温度只具有统计意义,通常而言不是定值.基于已报道工作 ,花岗岩的热破裂临界温度范围为 .结合文献报道及实验结果,推测文献 所用北山花岗岩的热破裂临界温度在 范围内.赵阳升等 对鲁灰花岗岩热破裂与渗透性相关规律的研究表明,在花岗岩热破裂临界温度附近,其渗透率将随温度升高(可达 的温升区间)而增大直至达到峰值,而后随温度继续升高而减小至峰谷.达热破裂临界值后的温升过程中,花岗岩的热破裂呈现间断性及多期性变化.鲁灰花岗岩在 之前有个热破裂剧烈期:、温度段,在热破裂剧烈期的温度段内,其渗透率均呈现一个峰值区间,在热破裂平静期的前期,渗透率依然维持较高水平,而后缓慢降至峰谷区间;其机制是随温度升高,花岗岩矿物颗粒热膨胀使先前热破裂产生的裂纹宽度减小,部分裂隙甚至闭合,裂隙整体连通性降低,使花岗岩的渗透率较峰值降低,但仍维持较高的渗透率水平,当另一个热破裂高峰出现时,渗透率又增大,如此经历数次热破裂积累后,使花岗岩的渗透率愈来愈大,也同步伴随着花岗岩的更加破裂.因此,随温度升高,花岗岩的渗透率并非单一线性变化,而存在一定的峰谷波动,但整体随温度升高呈上升趋势.)温度对S e()在北山花岗岩扩散影响的可能机理综上,基于文献 的数据结果,考虑温度对北山花岗岩结构的可能影响,对 S e()在北山花岗岩中的扩散随温度变化的机理推测如下:()范围内.随温度升高:一方面,溶液中的分子热运动增强,S e()扩散增强;另一方面,花岗岩的不同组成矿物发生持续热膨胀,导致花岗岩内部存在的原始裂隙被压缩甚至部分闭合,孔隙度降低,导致S e()的扩散被部分抑制.因此,受上述两方面因素协同影响,S e()在北山花岗岩中的有效扩散系数与温度并不呈单调线性关系,虽随温度升高而扩散增强,但整体受温度影响不显著.()范围内.此温度区间内达到实验用北山花岗岩的热破裂临界温度,北山花岗岩的原生裂隙进一步发育和扩展,并诱发新的热致裂隙,时裂隙进一步发育,已形成数量更多、尺度更大的热致裂隙,进而形成更多的孔隙裂隙通道,渗透率持续增大,从而使扩散溶液可通过的孔隙变多即扩散通道增多,进而促进S e()在北山花岗岩中的扩散,使 时S e()在北山花岗岩中的扩散较之 显著增强.同时,叠加随温度升高溶液分子热运动增强引起的S e()扩散增强,整体表现为 时S e()在北山花岗岩中的有效扩散系数显著增大,甚至比 高出一个数量级.()温度高于 .当北山花岗岩的温度进一步升高到 时,花岗岩的热破裂更为剧烈,该热破裂高峰会使北山花岗岩裂隙愈多、渗透率进一步增大,从而导致S e()在北山花岗岩中的扩散进一步增强.综合考虑热破裂剧烈程度与渗透率大小间的峰谷关系,推测随温度升高(),虽然S e()在北山花岗岩中的有效扩散系数在部分温度区间内可能会有波动,但整体会进一步持续增大.总结与展望本文对国内开展的S e在北山花岗岩中的扩散研究方法及获取的有关数据进行了总结和评述,重点分析了不同温度条件下 S e()在北山花岗岩中的可能扩散机理,旨在为国内同行提供有益参考.现有统计数据表明,在 ,S e()在北山花 岗 岩 中 的 扩 散 整 体 随 温 度 升 高 而 增强,特别是当温度达到 时,其有效扩散系数De比 时高出近一个数量级,扩散显著增强.结合不同温度下北山花岗岩形成因子Ff的数值变化及其热破裂临界温度(),推测当温度达 时,北山花岗岩已发生热破裂,其内部结构发生改变,渗透率明显变大,进而扩散显著增强.综合考虑目标元素S e与扩散介质花岗岩的属性随温度的变化情况,推测 S e()在北山花岗岩中的扩散随温度升高而加快,主要受S e分子热运动及北山花岗岩内部孔隙结构变化的综合影响:()分子热运动:温度升高加速了S e分子在花岗岩孔隙溶液中的热运动,扩散加快;()北山花 岗 岩 孔 隙 结 构 变 化:随 温 度 升 高,北山花岗岩不同组成矿物发生持续热膨胀,花岗岩内部的原始裂隙被压缩甚至部分闭合,核化学与放射化学第 卷S e()的扩散被部分抑制;,北山花岗岩发生热破裂,原生裂隙进一步发育和扩展,并诱发新的、数量更多、尺度更大的热致裂隙,形成更多孔隙裂隙通道,渗透率持续增大,S e()的扩散显著增强;()预测随温度继续升高,S e()在北山花岗岩中的扩散整体将进一步增强(部分温度区间可能会有波动),De甚至将继续产生数量级上的差异.为了在上万年的时间尺度内对R N在高放废物地质处置库中的迁移行为进行预测,必须建立相关模型,而参数获取是模型建立的根基.有效扩散系数作为处置库安全评价的最关键参数之一,在高放废物地质处置上万年的评价周期内,其数值变化对评价结果的灵敏度将产生重要影响,获得尽可能接近处置库实际条件下的数值结果对保证评价的有效性至关重要.采用 温度条件下的De值来评估 S e在北山花岗岩中的迁移行为,可以为潜在处置库的安全预评价提供更安全/保守的结果.特别是在接近废物桶的近场环境下,高放废物释热导致的花岗岩温度升高可能会达到花岗岩的热破裂临界温度,而热破裂形成的裂隙是核素迁移的关键通道,可显著促进核素的扩散迁移,因此处置库环境下花岗岩的热破裂对核素在近场环境下的迁移影响应予以特别关注.对于该领域的后续工作,提出以下建议:()鉴于当前S e在北山花岗岩中扩散的考察温度最高达,与文献报道的可能近场温度 还存在一定差距,建议后续拓展高温扩散实验的高温范围,进一步评估高温对 S e在北山花岗岩中扩散的实际影响,为处置库安全评价提供更准确的数据输入;()对于温度对S e在北山花岗岩中扩散机理的探索,可结合可用的先进微观分析技术,对实验温度条件下花岗岩的内部结构进行追踪分析,验证上述推理的合理性,并对更广阔温度及实验条件下的结果给出更科学合理的预测,建立相关理化模型;()考虑拟进行深地质处置废物中的核素类型及S e在北山花岗岩中扩散的温度敏感性,可拓宽需要考察温度对其在北山花岗岩中扩散影响的R N种类(如T c、超铀核素等),评估温度影响的共性及差异性,系统总