棱化学与放射化学Ma.pdfVIP

笛 ，7卷 革 ，期 棱 化 学 与 放 射 化 学 2000 年 5 月 lournal of Nudear and Radiechemistry May，2009 文章编号 ：0253．995o(2o0o)o2—0080—07 花 岗岩体单裂 隙中核素迁移数学模型 Ⅲ ．扩散模型及其有 限单元法解 P6州 7 -v ， 叮 0f 苏 锐 ，李春江。，王 驹 ，高宏成 (1．校工业北京地质研究院 环保 中心，北京 100029；2．北京大学 技术物理系，北京 100871) 摘要 ：以北京市阳坊地区地表新鲜花 岗岩为实验介 质，在含有天然裂 隙面切 割的花岗岩柱 中进行 了核紊扩散实验。考虑 了花岗岩基质扩散、裂睬表 面吸附和基质 吸附等因素对核素扩敬 的影响， 着重探讨 了花岗岩基质扩敬对核素在裂隙中迁移的影响。建立 了相应 的二维扩敬数学模型，运甩 有限单元法求得其数值解。将数值解 、解析解和实验值进行 了对 比。三者较吻台 ，说明所建立的数 学模型能很好的反映实际实验条件 ，模型是正确的，所进行 的数值离散是合理的。由数学模 型估 算的迁移参数为 ：基质域中，扩散系数 D =5．915×10 ，s，孔隙度 e=0．019，阻滞系数 Rd= ， b 1．46×1o3；裂隙中，扩敬系数 Df=118×10 rn2／s，等效孔酵度 E=0舟7，表面 阻滞系数 R = 。 l3o．6。 关 建 中国分类词号：天 苎 ；，塑 j皇璺兰量 ， 两 ． ：P641 ；TL~)42 ．1 叟献标识鸡 ：A 文献 i1．2]对锶 、铯 、碘 、硒等在花 岗岩 中的扩散 、吸附、弥散 ．以及相应数学模型进行 了研 究。由于实验所用花 岗岩均被严重扰动 (破碎、研磨、劈开等)，增加了大量的人为因素 ，而且把 核素复杂的迁移过程简化为一维模型，很难模拟现场核素迁移过程 ，所获得的迁移参数不宜用 于现场场址评价研究。文献 [3，4]分别论述 了核素在天然单裂隙中迁移的脉冲模型和连续注 入模型。本文用等效多孔介质模型代替平行板模型模拟天然单裂踪系统，建立核素二维扩散 模型，并用有限元法求 出其数值解，探讨花岗岩天然单裂隙中核素扩散特征 ，了解花岗岩介质 中校素在不同水动力条件下的迁移特征。 1 物理模型 花 岗岩柱及简易实验装置纵剖面图示于图1。岩柱直径D=50rnrn、长 L=100lllr12，水平 天然单裂隙通过岩柱中心轴 ，裂隙隙径 26远小于直径D 和长度 L，岩柱经蒸馏水充分饱和。 收稿 日期 ：1999-09．14； 修订 日期：1999—12-24 基金项 目；援工业科学基金资助项 目(J94P61056) 作者筒介；苏 镜 (1974一)，男 ，四川成都人 ，研士研究生，高放废物地质处置专业 第 2期 苏 锐等 ：花岗岩体单裂隙中核素迁移数学模型Ⅲ 8l 一 一 藏 量 槽 0 。一槽 图 1 岩柱及实验装置 Fig1 Reekcylinderandexperimentalsetup 图2 几何模型 F 2 Geometrica[r】10[jd 假设 ： (1)把天然裂隙系统和岩块看作均匀多孔介质系统 ，如 图2所示 。 (2)裂 隙和岩块之 间只有溶质交换 ，无水量交换 。溶质交换量 由Fick第一定律给 出，即： r = 一 其