[发明专利]海水入侵与土壤盐渍化的监测评价方法

权利要求书

1.一种海水入侵与土壤盐渍化的监测评价方法，其特征在于，包括：

获取被测区域的土壤电导率ECe；

获取被测区域的地下水水位L和地下水电导率ECw；

根据地下水水位L和地下水电导率ECw，利用麦夸特法拟合地下水水位L与地下水电导率ECw关系，拟合生成土壤电导率拟合值

在被测区域选取N个监测点，分别获取每一个监测点的土壤电导率ECe_i、土壤电导率拟合值以及N个监测点的土壤电导率拟合结果平均值利用公式：

计算出海水入侵与土壤盐渍化灾害链效应的评价结果，其中，为评价因子。

2.根据权利要求1所述的海水入侵与土壤盐渍化的监测评价方法，其特征在于，根据所述评价因子的数值大小，判断海水入侵对土壤盐渍化的影响程度，且评价因子的数值越大，影响程度越高。

3.根据权利要求1所述的海水入侵与土壤盐渍化的监测评价方法，其特征在于，所述土壤电导率拟合值的拟合生成过程为：

首先，确定土壤电导率拟合值的拟合公式：

其中，P₁,P₂,…,P₁₁为拟合参数；

其次，对拟合参数P₁,P₂,…,P₁₁进行求解，步骤为：

(1)建立残差平方和函数：

其中，为第i个监测点的土壤电导率拟合值；ECe_i为第i个监测点的土壤电导率；L_i为第i个监测点的地下水水位；ECw_i为第i个监测点的地下水电导率；

(2)对所述拟合参数P₁,P₂,…,P₁₁赋初始值，记为P₁⁰,P₂⁰,…,P₁₁⁰；

(3)令P_j-P_j⁰＝Δ_j，j＝1,2,…,m，m＝11；

其中，P_j是第j个待定的拟合参数，P_j⁰是P_j的初始值；

(4)计算a_jk和a_jy；其中，

(5)根据最小二乘法原则，Δ_j,j＝1,2,…,m满足线性方程组：

其中，d为阻尼因子，取值范围为0～1，求解所述线性方程组，计算出一组拟合参数P₁',P'₂,…,P₁₁'；

(6)将计算出的所述一组拟合参数P₁',P'₂,…,P₁₁'作为初始值P₁⁰,P₂⁰,…,P₁₁⁰，返回步骤(3)进行迭代运算，并在每一次迭代运算过程中，根据计算出的拟合参数P₁',P'₂,…,P₁₁'计算残差平方和Q；当(Q_M-Q_M-1)/Q_M＜10^-3时，停止计算，将第M次迭代运算过程中计算出的拟合参数作为最终的拟合参数P₁,P₂,…,P₁₁；其中，Q_M为第M次迭代运算的残差平方和。

4.根据权利要求3所述的海水入侵与土壤盐渍化的监测评价方法，其特征在于，所述N>10。