[发明专利]一种土壤重金属浓度随时间变化的估算方法

权利要求书

1.一种土壤重金属浓度随时间变化的估算方法，其特征在于：该浓度随时间变化的估算方法的形式为：

C_t＝C₀-E*t (1)

E＝λ*C₀*R^a1*K^a2*LS^a3*G^a4*P^a5/z (2)

式(1)中：C_t为土壤中某重金属在t年后的浓度，mg/kg；C₀为其初始值，mg/kg；E为重金属年流失量，mg/kg/年；

式(2)中：λ为坡地养分传输系数；R为降雨侵蚀力因子，MJ*mm/(hm²*h*a)；K为土壤可蚀性因子，(t*hm²*h)/(hm²*MJ*mm)；LS为地形因子；G为植被覆盖与管理因子；P为水土保持措施因子；z为转换系数，t*m^-2，取z＝8.1，LS、G、P为无量纲因子；a1、a2、a3、a4、a5为经验指数；

所述LS的计算方式为：

LS＝1.5*(H/3.84)^0.8

其中，H为区域海拔；

所述K值的计算方式为：

K＝{0.2+0.3exp[0.0256SAN(1-SIL/100)]}*[SIL/(CLA+SIL)]^0.3*[1-0.25TOC/(TOC+2.16)]*{1-0.7SN1/[SN1+exp(-5.51+22.9*SN1)]}

其中，SAN为砂粒含量，％，SIL为粉粒含量，％，CLA为黏粒含量，％，TOC为有机碳含量；SN1的计算公式为SN1＝1–SAN/100；

所述λ、R、G、P及a1、a2、a3、a4、a5取经验数值：

λ＝0.06，R＝200，G＝0.18，P＝0.6；

a1＝0.85，a2＝1.1，a3＝0.9，a4＝1.1，a5＝1.25。

2.采用权利要求1所述的土壤重金属浓度随时间变化的估算方法的土壤重金属预测方法，其特征在于：具体预测步骤如下：

步骤一：通过实验和查询得到土壤重金属初始值、区域海拔、土壤粒度和有机碳TOC含量的参数；

步骤二：将这些参数依次整理到excel表格内；

步骤三：在Matlab软件上编程，将浓度随时间变化的估算方法公式程序化；

步骤四：通过步骤三中已编好程序的Matlab软件自动读取并计算步骤二中输入至excel表格内的参数数据；

步骤五：最后输入预测年限，运行，得出结果；

步骤六：将步骤五中所得结果与实测值比较，验证拟合程度。

3.根据权利要求2所述的土壤重金属预测方法，其特征在于：所述土壤粒度、TOC含量通过实验测得，区域海拔通过GoogleEarth查询得到。