CN114792048B - 一种钚在矿物表面动态解吸模型的建立方法 - Google Patents

一种钚在矿物表面动态解吸模型的建立方法 Download PDF

Info

Publication number
CN114792048B
CN114792048B CN202210390800.3A CN202210390800A CN114792048B CN 114792048 B CN114792048 B CN 114792048B CN 202210390800 A CN202210390800 A CN 202210390800A CN 114792048 B CN114792048 B CN 114792048B
Authority
CN
China
Prior art keywords
desorption
plutonium
concentration
function
dynamic
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
Application number
CN202210390800.3A
Other languages
English (en)
Other versions
CN114792048A (zh
Inventor
柴娜娜
梁威
王卫宪
徐辉
汤海松
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
63653 Troops of PLA
Original Assignee
63653 Troops of PLA
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 63653 Troops of PLA filed Critical 63653 Troops of PLA
Priority to CN202210390800.3A priority Critical patent/CN114792048B/zh
Publication of CN114792048A publication Critical patent/CN114792048A/zh
Application granted granted Critical
Publication of CN114792048B publication Critical patent/CN114792048B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/02Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种钚在矿物表面动态解吸模型建立的方法,通过开展不同流速下钚的动态解吸实验,获取钚解吸实验的关键参数,推导钚在矿物表面解吸动力学反应方程,建立钚的动态解吸模型,本方法原理简单,可适用于多种场景下核素钚在各岩矿表面的解吸行为研究,并根据该解吸模型获得核素在不同流速下随时间变化的解吸量,为实现大尺度时间范围内钚的解吸行为预测提供支撑。

Description

一种钚在矿物表面动态解吸模型的建立方法
技术领域
本发明属于核环境科学领域,涉及一种钚在矿物表面解吸模型的建立方法,该方法可应用于地下水环境体系、放射性处置库等涉及核素吸附/解吸行为领域的研究和预测。
背景技术
放射性废物是人们利用核能过程中产生的主要废物,它具有放射性活度大(占乏燃料后处理废物中98%以上的放射性),比活度大(≥1010Bq/kg),半衰期长(239Pu-T1/2:2.44×104年;237Np-T1/2:2.14×106年),毒性大等特性;若处置不当将严重危害人类生存环境,该问题尤其在使用核能的国家引起了广泛的关注,并投入了大量的人力、物力和财力。目前普遍认可的做法就是深地质处置,即通过人为设置的工程屏障和天然屏障阻滞放射性废物迁移到环境中去。但是,由于地壳运动等不可控因素有可能导致地下水对处置场进行浸泡,从而导致核素-水-矿物三者之间发生相互作用,引发核素的泄露。这种作用包括氧化还原、吸附/解吸、溶解/沉淀等,大量学者研究核素在各种矿物、黏土表面的吸附行为,目的就是为了了解各矿物体系与放射性核素的相互作用,估计核素在地质材料中的迁移速度,为各种以地质材料为基础的处置库安全评价提供必要的参考。然而,光有吸附行为不足以预测核素的迁移,决定核素迁移量的是解吸行为,但是解吸反应机理和动态解吸模型在国内研究中很少,尤其是针对于核素钚的解吸,几乎没有,钚是高放废物中的关键核素之一,其特殊的化学性质和重大国防意义是处置库安全评价过程中及其重要且不可缺少的工作,因此,建立一种钚在矿物表面动态解吸模型,不仅有益于开展钚-岩矿-地下水相互作用研究,更是我国废物处理处置技术发展的需求,具有深远的实用价值。
钚的解吸模型建立存在以下几个问题:
①缺少可借鉴的解吸模型。国内针对核素在矿物表面的吸附行为研究很多,但是解吸行为研究很少且几乎没有涉及解吸机理和动态模型推导,这对于预测核素在地下水环境中的解吸量和迁移行为都非常不利。
②缺少钚解吸行为研究数据。国内针对钚的解吸行为研究非常少,大多数吸附/解吸研究都是用其它核素来间接的替代钚,例如Np、Eu等,但是众所周知,钚的化学行为是非常复杂的,这种用相似核素代替的做法并不能完全代表钚的解吸行为。
③尚未建立多种因素耦合影响下钚的解吸模型。目前,国内现有研究涉及到钚的解吸,大多是单变量下(例如:溶液pH值、Fe2+浓度等)钚的解吸量变化分析,脱离了场地的实际情况,没有整体考虑环境中氧化还原作用、流速等多种变量耦合情况下钚的解吸行为变化,也并未建立相应模型,更是无法预测钚的长期解吸行为。
发明内容
本发明针对现有技术不能满足钚在多种场景下的解吸行为研究和预测需求,提供一种钚动态解吸模型的建立方法,并在一系列解吸实验、解吸机理和动力学表达式的基础上,建立钚动态解吸的数学模型,为后续钚在多种场景环境下的解吸行为的预测提供有力支撑。
一种钚在矿物表面动态解吸模型的建立方法,包括以下步骤:
(1)实验样品制备:将吸附了钚的矿石解吸样品进行粉碎、研磨、制成粒径为200目的粉末,为保证钚的解吸行为涉及氧化还原,而非单一的物理吸附,因此选用的解吸样品保证预先吸附钚的时间至少达到半年以上。
(2)开展解吸实验:将上述样品称量一定质量并填装在解吸柱内通过解吸液进行解吸,解吸柱两端装有一定孔隙尺寸的聚四氟乙烯滤水板,再由螺纹堵盖和密封圈进行密封,螺纹堵盖两端接有硅胶管,分别作为解吸液注入口和接液口。解吸液在进入解吸柱之前要经过一个蠕动泵来改变流速,可以模拟地下水水流速度。
解吸液为NaHCO3-NaCl溶液,其中NaHCO3浓度0.5~0.7mM,NaCl浓度0.01~0.1M。控制解吸液流速为0.04~0.4mL/min,将解吸液不断注入解吸柱内进行钚动态解吸实验,接液口连接自动集液器收集解吸液,并通过质谱进行测量,获得一系列不同流速及解吸时间下的钚浓度,再通过放化分析流程获得解吸液中钚的价态分布,得到钚的主要存在形态为Pu(V+VI)。
(3)解吸数据整理:将获得的一系列钚浓度点与解吸时间作图,获得不同流速下的钚动力学解吸曲线。根据钚浓度和钚价态,推导出钚解吸机理。
(4)解吸过程分析及控制方程推导:根据质量守恒定律,动态解吸实验中Pu的总浓度等于流出浓度和固相吸附浓度之和,即:
C in C out —解吸实验中Pu流入和流出的浓度(mol/L);
j-解吸液流动速率(L/min);
V-流动解吸装置的体积(L);
M-固相解吸样品质量(g);
S-固相解吸样品表面吸附的Pu浓度(mol/g);
其中Pu的流入浓度为0,将上式写成浓度随时间变化的微分表达式:
同时,Pu在矿石表面的吸附行为符合Freundlich方程式,有:S=K×Cn,因此可以得到S和C之间的关系:
k 1 —正向动力学反应常数(/min);
k -1 —逆向动力学反应常数(/min);
根据动态柱解吸实验:
两边求导,得二阶微分方程式:
最终可以得到如下钚的动态解吸二阶微分等式。
(5)拟合数据获得试函数:将整个解吸曲线根据流速的不同划分为多段,文献报道中,解吸规律符合下式变化规律,因此选择此形式函数利用origin软件对解吸数据进行拟合,获得函数中的a、b、m值,将每段获得的函数表达式作为对应曲线的试函数;
本研究最终目的是获得微分方程中的动力学参数(k 1 k -1 n),因此将已知的试函数带入上述二阶微分方程后,得到误差函数f(t):
对得到的上述函数在整个解吸时间范围内积分,得到一个以k 1 、k -1 n为自变量的函数g(k 1 k -1 n):
(6)利用优化算法获得解吸参数:将g(k 1 k -1 n)函数作为目标函数,根据实验参数理论变化范围给出如下约束条件,建立如下优化模型M,并使用python skit-opt库中的遗传算法进行参数优化,获得了多组参数(k 1 k -1 n),此为多组局部最优解;将得到的参数代入二阶微分方程中进行求解验证,再通过Runge-Kutta算法求取该方程的数值解,用matlab中的ode函数实现这一算法,最终获得全局最优的二阶微分方程中的解吸参数,建立钚的动态解吸数学模型;
/>,其中,C-解吸液中的钚浓度(mol/L),j-解吸液流动速率(L/min),M-固相解吸样品质量(g)。
综上所述,本发明通过开展不同流速下钚的动态解吸实验,获取钚解吸实验的关键参数,推导钚在矿物表面解吸动力学反应方程,建立钚的动态解吸模型,本方法原理简单,可适用于多种场景下核素钚在各岩矿表面的解吸行为研究,并根据该解吸模型获得核素在不同流速下随时间变化的解吸量,为实现大尺度时间范围内钚的解吸行为预测提供支撑。
附图说明
图1为本发明连续条件下钚的动态解吸曲线;
图2为本发明连续条件下钚的解吸实验数据与模拟计算值。
具体实施方式
下面通过具体实施例对本发明钚在矿物表面动态解吸模型的建立方法进行详细说明。
一种钚在矿物表面动态解吸模型的建立方法,包括以下步骤:
(1)将吸附了钚的花岗岩矿石样品进行干燥、粉碎、研磨、制成粒径为200目的粉末,为确保钚的解吸行为涉及氧化还原反应而非单一的物理吸附,因此,选用的解吸样品保证预先吸附钚的时间至少达到半年以上。
(2)将上述样品称量约1g填装在解吸柱内,检验解吸柱密闭性。解吸柱两端装有一定孔隙尺寸的聚四氟乙烯滤水板,再由螺纹堵盖和密封圈进行密封,螺纹堵盖两端接有硅胶管,分别作为解吸液注入口和接液口。解吸液在进入解吸柱之前要经过一个蠕动泵来改变流速,可以模拟地下水水流速度。解吸液为NaHCO3-NaCl溶液,其中NaHCO3的浓度为0.5mM,NaCl的浓度为0.1M。
(3)通过外接蠕动泵控制流速为0.4、0.1和0.04ml/min,将解吸液不断注入解吸柱内进行速率控制下的钚动态解吸,接液端口连接自动集液器以便定时定量收集解吸液,每个流速下的解吸总时间依次设为1250min、5700min和15625min,通过质谱对解吸液中钚浓度进行测量,再通过放化分析流程获得解吸液中钚的价态分布,得到钚的主要存在形态为Pu(V+VI)。
(4)将获得的一系列钚浓度点与解吸时间作图C~t(图1),获得三个不同流速下的钚动态解吸曲线。根据钚浓度和钚价态,推导出钚解吸机理如下,反应1和3分别代表不同价态钚的动态吸附/解吸平衡,反应2代表不同价态钚之间的氧化还原转化,反应4代表钚在矿物表面的永久吸附。
(5)根据质量守恒定律,动态解吸实验中Pu的总浓度等于流出浓度和固相吸附浓度之和,即:
C in C out —解吸实验中Pu流入和流出的浓度(mol/L);
j-解吸液流动速率(L/min);
V-流动解吸装置的体积(L);
M-固相解吸样品质量(g);
S-固相解吸样品表面吸附的Pu浓度(mol/g);
其中钚的流入浓度为0,即:C in =0,将上式写成浓度随时间变化的微分表达式,:
同时,根据Pu在花岗岩表面的吸附行为符合Freundlich方程式,有:S=K Cn,因此可以得到S和C之间的关系:
k 1 —正向动力学反应常数(/min);
k -1 —逆向动力学反应常数(/min);
根据动态柱解吸实验:
两边求导,得二阶微分方程式:
最终得到如下钚的动态解吸二阶微分等式:
(6)将整个解吸曲线根据流速的不同划分为多段,文献报道中,解吸规律符合下式变化规律,因此选择此形式函数利用origin软件对解吸数据进行拟合,获得函数中的a、b、m值,将每段获得的函数表达式作为对应曲线的试函数;
本研究最终目的是获得微分方程中的动力学参数(k 1 k -1 n),因此将已知的试函数带入上述二阶微分方程后,得到误差函数f(t),且得到的f(t)并不满足目标函数的条件,即f(t)≠0
对得到的上述函数在整个解吸时间范围内积分,得到一个以k1、k-1和n为自变量的函数g(k 1 k -1 n):
(7)将g(k 1 k -1 n)函数作为目标函数,根据实验参数理论变化范围给出如下约束条件,建立如下优化模型M,并使用python skit-opt库中的遗传算法进行参数优化,获得了多组参数(k 1 k -1 n),此为多组局部最优解;将得到的参数代入二阶微分方程中进行求解验证,再通过Runge-Kutta算法求取该方程的数值解,用matlab中的ode函数实现这一算法,最终获得全局最优的二阶微分方程中的解吸参数,建立钚的动态解吸数学模型。根据解吸模型计算获得的钚解吸浓度能够很好的与实验值拟合(图2);
钚的动态解吸数学模型:
,其中,C-解吸液中的钚浓度(mol/L),j-解吸液流动速率(L/min),M-固相解吸样品质量(g)。

Claims (2)

1.一种钚在矿物表面动态解吸模型的建立方法,其特征在于,包括以下步骤:
(1)实验样品制备:将吸附了钚的矿石解吸样品进行粉碎、研磨、制成粒径为200目的粉末,解吸样品预先吸附钚的时间达到半年以上;
(2)开展解吸实验:将上述解吸样品填装在解吸柱内,控制解吸液流速并将其不断注入解吸柱内进行钚动态解吸,收集不同流速及解吸时间下的解吸液,并通过质谱进行测量,获得不同流速及解吸时间下解吸液中的钚浓度;
(3)解吸数据整理:将上述获得钚浓度与解吸时间作图,获得不同流速下的钚动力学解吸曲线;
(4)解吸过程分析及控制方程推导:根据质量守恒定律,钚动态解吸过程中Pu的总浓度等于流出浓度和固相吸附浓度之和,即:
C in C out —解吸过程中Pu流入和流出的浓度,单位mol/L;
j-解吸液流动速率,单位L/min;
V-流动解吸装置的体积,单位L;
M-固相解吸样品质量,单位g;
S-固相解吸样品表面吸附的Pu浓度,单位mol/g;
其中Pu的流入浓度为0,将上式写为浓度随时间变化的微分表达式:
Pu在矿石表面的吸附行为符合Freundlich方程式:S=K×Cn,得到S和C之间的关系:
k 1 —正向动力学反应常数(/min);
k -1 —逆向动力学反应常数(/min);
代入浓度随时间变化的微分表达式得:
根据动态柱解吸实验过程:
两边求导,得二阶微分方程式:
最终得到如下钚的动态解吸二阶微分等式:
(5)拟合数据获得试函数:将整个解吸曲线根据流速的不同划分为多段,解吸规律符合下式变化规律,因此选择此形式函数利用origin软件对解吸数据进行拟合,获得函数中的a、b、m值,将每段获得的函数表达式作为对应曲线的试函数;
将已知的试函数带入上述二阶微分方程后,得到误差函数f(t):
对得到的上述函数在整个解吸时间范围内积分,得到以k1、k-1和n为自变量的函数g(k 1
k -1 n):
(6)利用优化算法获得解吸参数:将g(k 1 k -1 n)函数作为目标函数,根据实验参数理论变化范围给出如下约束条件,建立如下优化模型M,并使用python skit-opt库中的遗传算法进行参数优化,获得了多组参数(k 1 k -1 n),此为多组局部最优解;将得到的参数代入二阶微分方程中进行求解验证,再通过Runge-Kutta算法求取该方程的数值解,用matlab中的ode函数实现这一算法,最终获得全局最优的二阶微分方程中的解吸参数,建立钚的动态解吸数学模型;
钚的动态解吸数学模型:
,其中,
C-解吸液中的钚浓度,单位mol/L;j-解吸液流动速率,单位L/min;M-固相解吸样品质量,单位g。
2.根据权利要求1所述一种钚在矿物表面动态解吸模型的建立方法,其特征在于:步骤(2)中,所述解吸液为NaHCO3-NaCl溶液,其中NaHCO3浓度0.5~0.7mM,NaCl浓度0.01~0.1M;控制解吸液流速为0.04~0.4mL/min。
CN202210390800.3A 2022-04-14 2022-04-14 一种钚在矿物表面动态解吸模型的建立方法 Active CN114792048B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210390800.3A CN114792048B (zh) 2022-04-14 2022-04-14 一种钚在矿物表面动态解吸模型的建立方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210390800.3A CN114792048B (zh) 2022-04-14 2022-04-14 一种钚在矿物表面动态解吸模型的建立方法

Publications (2)

Publication Number Publication Date
CN114792048A CN114792048A (zh) 2022-07-26
CN114792048B true CN114792048B (zh) 2024-04-02

Family

ID=82462411

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210390800.3A Active CN114792048B (zh) 2022-04-14 2022-04-14 一种钚在矿物表面动态解吸模型的建立方法

Country Status (1)

Country Link
CN (1) CN114792048B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2828361A1 (en) * 2011-02-25 2012-12-13 William Marsh Rice University Sorption and separation of various materials by graphene oxides
CN105023625A (zh) * 2015-06-10 2015-11-04 北京大学 放射性有机废液中微量铀和/或钚的回收方法
CN112285226A (zh) * 2020-10-16 2021-01-29 中国人民解放军63653部队 废液中Pu-239、Sr-90、Cs-137快速联合分析方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2828361A1 (en) * 2011-02-25 2012-12-13 William Marsh Rice University Sorption and separation of various materials by graphene oxides
CN105023625A (zh) * 2015-06-10 2015-11-04 北京大学 放射性有机废液中微量铀和/或钚的回收方法
CN112285226A (zh) * 2020-10-16 2021-01-29 中国人民解放军63653部队 废液中Pu-239、Sr-90、Cs-137快速联合分析方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
放射性废液中钚的耐辐射真菌吸附分离研究;徐辉;王耀芹;王卫宪;韩小元;梁威;翟秀芳;李伟平;王玮;孙应龙;李磊;袁超;王邦达;王毅;;原子能科学技术;20170520(第05期);第775-782页 *
环境土壤样品微波消解及钚的ICP-MS测量;沈忠等;《原子能科学技术》;20190130;第53卷(第04期);第603-610页 *

Also Published As

Publication number Publication date
CN114792048A (zh) 2022-07-26

Similar Documents

Publication Publication Date Title
McKinley et al. Radionuclide sorption database for Swiss safety assessment
Li et al. Study on dynamic adsorption characteristics of broken coal gangue to heavy metal ions under leaching condition and its cleaner mechanism to mine water
Jia et al. Investigation on gas drainage effect under different borehole layout via 3D monitoring of gas pressure
Jedináková-Křižová Migration of radionuclides in the enviroment
CN114792048B (zh) 一种钚在矿物表面动态解吸模型的建立方法
CN115165680A (zh) 一种铀尾矿渣-含水层耦合淋溶实验装置及方法
Windt et al. Modeling spent nuclear fuel alteration and radionuclide migration in disposal conditions
John et al. Development of New Processes for the Decommissioning Decontamination and for Treatment and Disposal of the Secondary Low-and Intermediate-Level Radioactive Waste
Miller et al. Simplified behaviors from increased heterogeneity: I. 2-D uranium transport experiments at the decimeter scale
Ferroni et al. Radionuclides migration in the far field of geological repositories: a numerical example
Kim et al. Review for mechanisms of gas generation and properties of gas migration in SNF (spent nuclear fuel) repository site
Bayliss et al. Radioelement behaviour in a cementitious environment
Langmuir Overview of coupled processes with emphasis in geochemistry
Patel et al. United States launches new direction to manage nuclear waste
Wang et al. Used Fuel Disposition in Crystalline Rocks: FY16 Progress Report
Fan et al. China: Experience of radioactive waste (RAW) management
Neretnieks The impact of alpha-radiolysis on the release of radionuclides from spent fuel in a geologic repository
Zheng et al. Coupling flow and electric fields to simulate migration and remediation of uranium in groundwater remediated by electroosmosis and a permeable reactive bio-barrier
Zhang et al. Adsorption-desorption of 241Am (Ⅲ) on montmorillonite colloids and quartz sand: Effects of pH, ionic strength, colloid concentration and grain size
Clayton et al. BNFL Lysimeter Proaramme to Investigate the Leachina of Radlonuclides from Low-Level Radioactive Waste
Longcheng Oxidative dissolution of spent fuel and release of nuclides from a copper/iron canister. Model developments and applications
Lucille et al. Modeling of geochemical compatibility of near field materials in terms of radionuclide retention properties
Brewitz et al. Long-term performance of spent fuel in geological repositories
Liu Oxidative dissolution of spent fuel and release of nuclides from a copper/iron canister. Model developments and applications
Tewhey et al. Studies of radionuclide availability and migration at the Nevada Test Site relevant to radioactive waste disposal

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