钾是农作物生长必需的三大营养元素之一^[1].我国呈现出不同程度的土壤缺钾现状, 且南方较北方的缺钾情况严重, 如广东的水稻土壤、旱地土壤缺钾面积均达到了90%以上^[2].同时, 可利用的水溶性钾矿资源短缺, 供给量仅占需求量的30% ~ 40%, 进口已成为钾肥的主要来源^[3].我国的非水溶性钾矿资源却非常丰富, 达到了100亿t^[4].充分开发利用这部分资源, 将极大缓解我国的缺钾现状.我国对非水溶性钾矿已做了较多的研究^[5-6], 然而大部分的提钾工艺复杂, 成本高, 难以推广.近年来, 基于活化概念的理化促释技术提供了一条新型研发思路^[7].已有的研究结果表明, 活化后钾长石的钾释放量显著提高, 可部分替代氯化钾而肥效不减^[8-9].开展活化钾矿的钾素释放机理与规律的理论研究, 尤其是对动力学规律的研究, 对其肥力评价指标的建立具有重要意义.

钾的释放既受含钾矿物类型的影响, 也受钾释放的溶液环境(如各种离子种类和浓度)的影响^[10].范钦桢^[11]发现铵态氮肥中的NH₄⁺会抑制土壤中非交换态钾和结构钾的释放, 土壤中常见的阳离子有NH₄⁺、Na⁺、Ca²⁺等, 而这些阳离子同样会影响非水溶性钾矿的释放.王瑾^[12]研究了不同阳离子盐溶液对黑云母、白云母、正长石等非水溶性钾矿钾释放的影响, 发现随着离子种类、钾矿类型的不同, 钾的释放也呈现不同的释放特征.本研究选取2种钾矿及其QN活化钾矿为研究材料, 采用NH₄Cl和NaCl这2种阳离子盐溶液做浸提剂对其连续振荡提取, 建立活化钾矿的钾释放动力学模型, 旨在从动力学角度研究活化钾矿的高效释放特征, 为评价活化钾矿的植物有效性提供理论依据.

1.   材料与方法

1.1   材料

供试的非水溶性钾矿为钾长石和富钾页岩, 其中, 钾长石w (K₂O)为8.57%, 取自广东五华; 富钾页岩w(K₂O)为11.42%, 取自河北张家口.2种钾矿经风干、磨细后过100目筛备用.

活化钾矿的制备方法:分别称取上述钾矿20 g, 加入质量分数为5%的QN活化剂, 再加入2 mL蒸馏水, 混合研磨5 min, 风干、磨细, 过100目筛备用.其中, QN活化剂为含Na⁺、不含K⁺的无机活化剂.

浸提剂分别为10 mmol·L^-1的NH₄Cl和NaCl.

1.2   方法

准确称取钾矿及活化钾矿0.500 0 g于离心管中, 分别加入50 mL不同的浸提剂溶液, 对照加入去离子水(H₂O), 摇匀, 在振荡机上振荡15 min, 取出后5 000 r·min^-1离心.倒出全部上清液, 用火焰光度法测定溶液钾的含量.残渣中分别加入50 mL上述溶液, 重复浸提步骤, 钾长石、活化钾长石连续提取10次, 页岩、活化页岩连续提取15次.每个处理设3个重复.上述提取次数均根据实际浸提过程中到达平衡附近的时间确定.

1.3   释放动力学模型

一级动力学模型:y = a-ae^-bx,

双常数模型:y = ax^b,

扩散模型:y = a+ bx^0.5,

Elovich模型:y = a+ blnx.

上述模型中, x为浸提时间, y为钾矿的累积释钾量, a、b为模型常数^[13-16].

1.4   数据处理方法

数据的处理、分析和制图分别采用Excel、Spass13.0、Matlab7.1等软件.

2.   结果与分析

2.1   连续振荡条件下活化钾矿的钾释放

如图 1所示, 各浸提剂浸提钾矿的释放均表现为前期快速, 之后缓慢释放的变化趋势, 其中NH₄Cl浸提时, 钾矿及活化钾矿在30 min左右即完成了快速释放, 进入了缓慢释放阶段.NH₄Cl和NaCl在浸提钾长石、活化钾长石时, 均在45 min左右达到了缓慢释放阶段; H₂O在浸提富钾页岩、活化页岩时, 快速释放阶段为0 ~ 90 min, 90 min之后为缓慢释放阶段, 而NaCl浸提到45 min左右时, 富钾页岩、活化页岩即进入了缓慢释放阶段.

图  1  不同浸提剂下活化钾矿的钾释放变化

Figure  1.  Variations of K releases of activated insoluble potassium ores with different cation solutions

下载: 全尺寸图片 幻灯片

连续浸提的过程中, 活化钾矿与钾矿表现出一致的浸提规律.初始阶段, 钾长石、活化钾长石的钾释放量表现为:NH₄Cl＞NaCl＞H₂O, 随着浸提时间的延长, 释钾量逐渐减少, 到达释钾平衡附近时, 3种浸提剂的累积释钾量表现为:NaCl＞NH₄Cl≈H₂O; 富钾页岩、活化页岩的整个动态释钾过程均表现为: H₂O＞NaCl＞NH₄Cl.

对钾矿及其活化钾矿的累积释钾量分析可知, 在H₂O、NH₄Cl、NaCl浸提下, 活化钾长石的累积释钾量分别是钾长石的2.3、2.0和1.7倍, 活化页岩的累积释钾量分别是富钾页岩的2.5、3.3和2.5倍, 所以活化钾矿的累积释钾能力大于钾矿.由图 1还可以看出, 富钾页岩的累积释钾能力大于钾长石、活化页岩的累积释钾能力大于活化钾长石.

2.2   钾矿的释钾动力学模型

由图 1中钾矿释钾的动态数据, 建立不同浸提剂钾矿的钾释放动力学模型, 拟合结果如表 1所示.其中, 模型拟合的优劣取决于拟合性, 即计算值与实测值的符合程度, 常用相关系数(R)和标准差(S)来评定, R越大、S越小拟合性越好.由表 1可以看出, 除钾长石的一级动力学拟合方程R达显著水平外, 其余拟合方程的R均达到极显著水平, 相关系数在0.698 6 ~ 0.997 3之间.累积释钾量的计算值与实测值之间的S在11.35 ~ 295.90之间, 表明4个模型均能很好的拟合钾矿及活化钾矿的动态释钾过程.

表  1  连续振荡条件下活化钾矿释钾的动力学模型¹⁾

Table  1.  The kinetics model of K release of activated potassium ores with successive extraction

下载: 导出CSV 

| 显示表格

双常数模型、一级动力学模型、扩散模型和Elovich模型拟合钾长石的累积释钾量, 拟合R的平均值分别为0.971 3、0.832 2、0.953 3和0.978 2;拟合S的平均值分别为15.48、22.99、19.49和13.30, 双常数模型和Elovich模型两者间的R和S基本没有差别, 所以, 钾长石的最优释放动力学模型是Elovich模型或双常数模型.

对于活化钾长石, 4种模型拟合R的平均值分别为0.948 4、0.920 1、0.900 3和0.958 1;拟合S的平均值分别为22.35、26.28、31.07和20.50, 与钾长石的结果类似, 活化钾长石的最优释放动力学模型亦为Elovich模型或双常数模型.

对页岩, 4种模型拟合R的平均值分别为0.979 4、0.906 4、0.974 6和0.981 7;S的平均值分别为:64.68、88.84、72.36和51.36, 钾的释放动力学模型拟合性表现为:Elovich模型＞双常数模型＞扩散模型＞一级动力学模型.所以, 页岩的最优释放动力学模型为Elovich模型.

对活化页岩, 4种模型拟合R的平均值分别为0.974 3、0.849 7、0.956 2和0.981 6;S的平均值分别为126.70、194.30、169.11和95.79, 动力学方程拟合性:Elovich模型＞双常数模型＞扩散模型＞一级动力学模型.所以, 活化页岩的最优释放动力学模型是Elovich模型.

4种模型拟合时, H₂O、NH₄Cl、NaCl浸提钾长石的方程S均表现为:NaCl＞H₂O＞NH₄Cl, 活化钾长石亦表现出相同规律.一级动力学模型拟合时, 3种浸提剂之间的S相差不大, 说明除了一级动力学模型, 其他3种动力学模型拟合钾长石、活化钾长石时, NH₄Cl浸提下的模型拟合性均优于NaCl.4种模型拟合页岩时浸提剂间的S均表现出:H₂O＞NaCl＞NH₄Cl, 浸提活化页岩时亦表现出相同规律, 说明4种模型拟合下, NH₄Cl浸提页岩、活化页岩的拟合性均优于NaCl.

2.3   动力学模型参数与钾矿释钾关系

由上述分析可知, 双常数模型、一级动力学模型、扩散模型和Elovich模型均具有较好的拟合性, 模型中的参数对于活化钾矿中钾素在盐溶液持续作用下的释放特征具有重要的意义.

拟合方程的参数见表 2.双常数方程的参数a表示释放过程的初始瞬时速率^[17-18], a值越大, 钾矿释钾的初始瞬时速率越大.通过比较双常数方程的a值可以看出, 活化钾矿的钾初始释放速率显著高于未活化钾矿, 在H₂O、NH₄Cl和NaCl的浸提条件下, 活化钾长石的钾初始释放速率分别比钾长石增加了3.7、1.1和1.6倍; 活化页岩的钾初始释放速率分别比富钾页岩增加了4.7、6.0和4.0倍.

表  2  活化钾矿的动力学模型拟合参数

Table  2.  The kinetics model fitting parameters of activated potassium ores

下载: 导出CSV 

| 显示表格

对Elovich方程求导, 可得到Elovich速率方程: y = b/x, 由b值可以求出任一时间的释放速率^{[17, 19-20]}.b值越大, 钾矿释钾的速率越大.表 2可以看出, NH₄Cl浸提钾矿及活化钾矿的b显著小于其他浸提剂, 表明NH₄Cl浸提的钾释放速率远小于其他浸提剂.比较钾长石和活化钾长石的b发现, 除了NaCl浸提活化钾长石的b小于钾长石, H₂O和NH₄Cl浸提活化钾长石的b分别比钾长石增加了10.3%和71.2%.H₂O、NH₄Cl和NaCl浸提活化页岩的b分别比页岩增加了61.7%、102.4%和86.3%.

一级动力学方程的a值表示释放过程的最大平衡释放量^{[16, 22-23]}, 由表 2可以看出, H₂O、NH₄Cl和NaCl连续浸提活化钾长石钾的一级动力学方程的a分别比钾长石增加了139.6%、99.6%和73.3%;活化页岩钾的a分别比富钾页岩增加了1.3、2.4和1.6倍.

动力学参数的分析可以看出, 活化后钾矿中钾的初始释放速率、平均释放速率以及最大平衡释放量均有所提高.其中, 初始释放速率增大是活化钾矿中钾的最大平衡释放量增多的主要原因.

3.   讨论与结论

各浸提剂下, 钾矿的累积释钾能力及动力学参数均表现出了富钾页岩大于钾长石、活化页岩大于活化钾长石, 说明矿物类型不同, 其释钾能力差异较大.这可能与矿物结构有关, 钾长石矿物结构为无水架状结构铝硅酸盐矿物, 钾原子的位置位于晶格内部, 钾原子落在10个氧原子所组成的穴中, 与6个氧原子相距0.285 nm, 因此阻碍了钾的释放^[10].页岩结构较为复杂, 如辽宁省朝阳地区的页岩包含了多种以独立矿物形式出现的含钾矿物, 主要成分是钾长石, 约80%的钾赋存于钾长石中, 其次为白云母、伊利石等, 约20%的钾赋存于云母类矿物中^[24], 而长石类与云母类相比, 由于钾离子处在相邻四面体的空隙中, 释钾较为困难^[15].

本研究结果表明, 钾长石、活化钾长石的钾初始释放量表现为:NH₄Cl＞NaCl＞H₂O, 累积释钾量表现为:NaCl＞NH₄Cl≈H₂O.一开始, NH₄⁺的提取能力强于Na⁺, 之后提取能力减弱, 这可能是因为与Na⁺、Ca²⁺等水化半径较大的离子相比, NH₄⁺与K⁺有几乎相同的离子半径和水化能^[25-26], 甚至在电性、化合价、释放与固定机制等多方面都具有相似之处^[27], 故NH₄⁺更容易置换矿物表面及边缘、楔形位点吸附的钾, 当矿物的速效钾含量较高时, 交换能力最强^[28].随着浸提时间的延长, 矿物表面和边缘处吸附的钾较少, 水化半径较大的Na⁺提钾能力大于NH₄⁺, 这是因为Na⁺虽然因其水化半径较大不易置换楔形位点上吸附的K⁺, 但能置换一部分矿物晶层表面吸附的K^{+ [29-30]}; Na⁺、Ca²⁺等水化半径较大的盐离子可以撬开矿物晶层, 使得一部分易释放的非交换性钾释放出来^[28].本试验中, H₂O浸提下页岩、活化页岩的累积释钾量大于2种盐溶液, 这可能是因为页岩复杂的结构中存在大量的层状结构, K⁺通过直接的扩散而不需要离子交换就可以释放出来.

通过动力学模型来拟合钾矿释钾的过程并描述其释钾规律鲜见报道.王瑾等^[15]研究发现用不同有机酸连续浸提黑云母、正长石等钾矿, 其释放的最优动力学模型为双常数模型或Elovich模型.为了准确、定量地描述某种元素的动态释放过程, 有必要建立或引用各种数学模型, 吕晓男等^[30]通过多个动力学模型拟合电超滤方法下土壤钾释放的动态过程发现, Elovich方程的参数b和双常数方程的参数a与土壤速效钾和大麦相对产量之间存在显著或极显著相关.本试验建立了活化钾矿的释放动力学模型, 同时对动力学模型的参数进行了比较, 发现活化后钾矿的钾初始释放速率、平均释放速率、最大平衡释放量均显著增大, 并从动力学角度定量评价了活化后钾矿的释钾效果.至于参数能否作为评价其生物有效性的指标还需进一步的研究.

本试验得到如下结论:1)浸提剂浸提活化钾长石的钾初始释放量表现为:NH₄Cl＞NaCl＞H₂O; 累积释放量表现为:NaCl＞NH₄Cl≈H₂O; 浸提活化页岩的整个动态释钾过程均表现为:H₂O＞NaCl＞NH₄Cl.在盐溶液的连续振荡浸提下, 活化钾矿的累积释钾能力大于钾矿, 富钾页岩及活化页岩的累积释钾能力大于相同处理下的钾长石.2)钾长石、活化钾长石钾的最优释放动力学模型是Elovich模型或双常数模型; 页岩、活化页岩的最优动力学模型是Elovich模型.除了一级动力学模型拟合钾长石、活化钾长石时, NH₄Cl、NaCl浸提剂间拟合性无差异外, 其他情况时, NH₄Cl浸提下的模型拟合性均优于NaCl.3)通过动力学模型参数比较, 从动力学角度定量评价了钾矿的QN活化效果, 结果表明, 活化后钾矿的钾初始释放速率、平均释放速率以及最大平衡释放量均有显著提高.