[发明专利]一种生态类自然保护地的保护成效评估方法及系统

权利要求书

1.一种生态类自然保护地的保护成效评估方法，其特征在于，包括：

获取所述自然保护地的第一时刻的状态和第二时刻的状态；所述第一时刻为监测评估周期的初期，所述第二时刻为监测评估周期的末期；所述监测评估周期为对自然保护地进行保护成效确定过程中监测自然保护地的状态的周期；

获取表征所述第一时刻的状态和所述第二时刻的状态的指标；

根据所述指标的变量，计算所述自然保护地保护成效评估的多种参数的变化率；

将每种参数的变化率代入对应的预先构建的指数模型中，得到各种参数指标的保护成效指数；

将各保护成效指数求和得到综合保护成效值；

所述根据所述指标的变量，计算所述自然保护地保护成效评估的多种参数的变化率，具体包括：

利用公式计算景观格局有效性指数，其中x₁为景观格局有效性指数；EP为景观格局指数，t为第一时刻；t+1为第二时刻；m₁为保护性景观镶嵌体总数；i₁，j₁均为保护性景观镶嵌体的顺序号；S_i1为第i₁个保护性景观镶嵌体的面积，S_j1为第j₁个保护性景观镶嵌体的面积；TS为自然保护地的总面积；FI为景观单元破碎化指数，介于0-1之间，FI的值越大，保护性景观镶嵌体破碎化程度越大，反之越小；

利用公式计算生态系统转化有效性指数，其中x₂为生态系统转化有效性指数；为评估初期第i₂类生态系统转化为评估末期第j₂类生态系统的面积；m₂、n₂分别为第一时刻和第二时刻生态系统类型的数量；Q为不同生态系统类型的重要性赋值；

当生态系统为非湿地生态系统时，利用公式计算生态系统质量有效性指数，其中x₃₁为生态系统为非湿地生态系统时的生态系统质量有效性指数，n₃为像元数量，为像元i₃的用遥感手段获得的归一化植被指数变化率，NDVI_st为第一时刻像元值被生长季的归一化植被指数值，NDVI_s(t+1)为第二时刻像元植被生长季的归一化植被指数值；当生态系统为湿地生态系统时，利用公式计算生态系统质量有效性指数，其中x₃₂为生态系统为湿地生态系统时的生态系统质量有效性指数，为第一时刻湿地丰水期的面积或水质指标的第i′₃个变量，为第二时刻湿地丰水期的面积或水质指标的第i′₃个变量；i′₃取值为1、2、……、m′₃，m′₃为湿地水面、水质的变量个数；

利用公式计算森林生长有效性指数，其中x₄为森林生长有效性指数；为森林生长有效性指数第i₄项次级指标变化率，次级指标包括郁闭度E_y和蓄积量E_x，其中蓄积量根据胸径和树高，查阅二元材积表确定；为次级指标权重，w₁取0.4，w₂取0.6；

利用公式计算植被稳定有效性指数，其中x₅为森林植被稳定有效性指数，H为物种多样性的测算指标，为物种i₅的个体数占总个体数的比例，i₅＝1,2,…,n₅；n₅为物种总数；

利用公式计算保护对象物种数量有效性指数，其中x₆为保护对象物种数量有效性指数；为自然保护地第i₆种保护物种的权重，根据保护物种的相对重要性程度确定；n₆为自然保护地保护物种的种类数；为第一时刻第i₆种保护物种的种群数量或分布面积；为第二时刻第i₆种保护物种的种群数量或分布面积；

利用公式计算干扰变化有效性指数，其中，x₇为干扰变化有效性指数；I为人类活动干扰指数，为功能区系数，即人类干扰土地类型处于自然保护地第j₇功能区的干扰赋值；为第i₇种人类干扰土地类型的干扰权重，自然保护地人类干扰土地类型包括采矿用地、工业用地、交通运输用地、旅游用地、住宅用地和农业生产用地，权重值分别取0.42、0.22、0.16、0.12、0.05和0.03；为第i₇种人类干扰土地类型在功能区j₇中的面积；

所述将每种参数的变化率代入对应的预先构建的指数模型中，得到各种参数指标的保护成效指数，具体包括：

将所述景观格局有效性指数代入景观格局有效性指数模型F(x₁)＝15/(1+EXP(-1.635-0.757*x₁))中，得到景观格局保护成效指数；

将所述生态系统转化有效性指数代入生态系统转化有效性指数模型F(x₂)＝15/(1+EXP(-1.635-1.557*x₂))中，得到生态系统转化保护成效指数；

当生态系统为非湿地生态系统时，将所述生态系统质量有效性指数代入非湿地生态系统质量有效性指数模型F(x₃₁)＝18/(1+EXP(-1.21-0.69*x₃₁))中，得到非湿地生态系统质量保护成效指数；当生态系统为湿地生态系统时，将所述生态系统质量有效性指数代入湿地面积变化指数模型F(x₃₂)＝10/(1+EXP(-1.35-0.65*x₃₂))中，得到湿地面积保护成效指数；

将所述森林生长有效性指数代入森林生长有效性指数模型F(x₄)＝10/(1+EXP(-0.73-0.54*x₄))中，得到森林生长保护成效指数；

将所述植被稳定有效性指数代入植被稳定有效性指数模型F(x₅)＝10/(1+EXP(-1.83-0.81*x₅))中，得到植被稳定保护成效指数；

将所述保护对象物种数量有效性指数代入保护对象物种数量有效性指数模型F(x₆)＝18/(1+EXP(-1.23-0.69*x₆))中，得到保护对象物种数量保护成效指数；

将所述干扰变化有效性指数代入干扰变化有效性指数模型F(x₇)＝14/(1+EXP(-1.135+0.71*x₇))中，得到人类活动干扰指数。