一种流域生态输水成效评估方法

本发明涉及生态环境治理，具体涉及一种流域生态输水成效评估方法。

背景技术：

1、干旱区内陆湖泊研究表明，生态输水对维持干旱区湖泊面积，以及恢复当地生态有着直接联系。为保护干旱区生态安全，可以在塔里木河、石羊河和黑河流域等地实施生态输水。生态输水工程特别是在塔里木河流域等干旱地区，在恢复和维持世界最干旱地区之一的生态平衡方面发挥着关键的作用。现有方法通过研究塔里木河生态输水影响下的湖泊面积、地下水位和植被的变化，发现生态输水工程使得水体和植被面积增加，植被变化明显，沙地面积减少。而对于塔里木河生态输水后对湖泊生态环境质量变化的研究主要有以下三个方面：(1)基于植被响应，从天然植被的净初级生产力(npp)、总初级生产力(gpp)、植被水分利用效率(wue)、归一化植被指数(ndvi)等生态系统功能指标的变化来评价输水后植被恢复情况；(2)基于水文响应，从地下水位埋深变化、地下水质变化、地表水体面积变化等方面来评估输水后流域水文状况的改变。(3)从土壤养分、土壤湿度来评估输水后土壤质量的状况。这些研究发现，生态输水后流域局部地区的天然植被物种数量、长势、植被覆盖度逐渐恢复，生物多样性不断增加；多数输水区水文条件明显改善，水质状况很大提升，地下水位明显抬升；土壤均呈湿润化发展，距河道越近的地方土壤养分含量越高。

2、虽然这些研究共同强调了生态输水工程在恢复和维持干旱地区生态平衡中发挥的关键作用，但是对于台特玛湖生态环境的变化综合评估并没有一个完整体系，无法对当地生态环境的变化进行量化比较，同时没有结合台特玛湖面积改变的驱动因素进行研究，难以建立生态输水与台特玛湖之间的相互关系，从而导致生态输水成效评估不准确。

技术实现思路

1、针对现有技术中的上述不足，本发明提供了一种流域生态输水成效评估方法，通过引入湿度、绿度、干度以及热度指标，构建用于评估流域生态输水成效的遥感生态指数模型，用以解决难以建立生态输水与研究区域即台特玛湖之间的相互关系以及生态输水成效评估不准确的问题。

2、为了达到上述发明目的，本发明采用的技术方案为：

3、一种流域生态输水成效评估方法，包括以下步骤：

4、s1、获取研究区域的流域遥感影像与历史实测水文数据并进行数据校正，计算流域的湿度、绿度、干度以及热度指标；

5、s2、根据步骤s1中计算的流域的湿度、绿度、干度以及热度指标，采用主成分分析法对各指标赋予权重，构建用于评估流域生态输水成效的遥感生态指数模型；

6、s3、根据步骤s2中构建的用于评估流域生态输水成效的遥感生态指数模型，获取不同年份的用于评估流域生态输水成效的遥感生态指数，并采用非参数化趋势度分析方法进行拟合效果分析，得到生态环境质量变化趋势，用于评估流域生态输水成效。

7、进一步地，步骤s1具体包括：

8、s11、获取研究区域的流域遥感影像与历史实测水文数据，利用历史实测水文数据对流域遥感影像进行数据校正，得到校正的流域遥感影像数据；

9、s12、根据步骤s11中得到的校正的流域遥感影像数据，计算绿度指标，即：

10、ndvi＝(nir-red)/(nir+red)

11、其中，ndvi表示绿度指标，nir表示流域遥感影像数据的近红外光，red表示流域遥感影像数据的红光波段；

12、s13、根据步骤s11中得到的校正的流域遥感影像数据，计算热度指标，即：

13、lst＝α×dn-β

14、其中，lst表示热度指标，dn表示流域遥感影像数据的陆地地表温度的温度值，α、β表示常量系数；

15、s14、根据步骤s11中得到的校正的流域遥感影像数据，计算湿度指标，即：

16、wetness＝α1ρ1+α2ρ2+α3ρ3+α4ρ4+α5ρ5+α6ρ6+α7ρ7

17、其中，ρ1、ρ2、ρ3、ρ4、ρ5、ρ6、ρ7分别表示流域遥感影像数据的各地表反射率波段，α1、α2、α3、α4、α5、α6、α7分别表示各地表反射率波段对应的常量系数；

18、s15、根据步骤s11中得到的校正的流域遥感影像数据，计算干度指标，即：

19、

20、其中，ndbsi表示干度指标，bsi表示裸土指数，ibi表示建筑指数，s1表示流域遥感影像数据的短波红外光，red表示流域遥感影像数据的红光，blue表示流域遥感影像数据的蓝光波段，green表示流域遥感影像数据的蓝光。

21、进一步地，步骤s2具体包括：

22、s21、将湿度、绿度、干度以及热度指标进行归一化处理，得到归一化处理的湿度、绿度、干度以及热度指标；

23、s22、对步骤s21中归一化处理的湿度、绿度、干度以及热度指标进行掩膜处理，且进行主成分分析转换，得到第一主成分结果，利用1减去第一主成分结果，构建用于评估流域生态输水成效的最初的遥感生态指数模型，即：

24、rsei0＝1-{pca1[ρ(ndvi′,lst′,wetness′,ndbsi′)]}

25、其中，rsei0表示用于评估流域生态输水成效的最初的遥感生态指数，pca1表示第一主成分结果，ρ表示掩膜处理，ndvi′、lst′、wetness′、ndbsi′分别表示归一化的绿度、热度、湿度以及干度指标；

26、s23、对步骤s22中构建的用于评估流域生态输水成效的最初的遥感生态指数进行归一化处理，得到用于评估流域生态输水成效的遥感生态指数模型，即：

27、rsei＝(rsei0-rseimin)/(rseimax-rseimin)

28、其中，rsei表示用于评估流域生态输水成效的遥感生态指数，rseimin表示遥感生态指数的最小值，rseimax表示遥感生态指数的最大值。

29、进一步地，步骤s3具体包括：

30、s31、根据步骤s2中构建的用于评估流域生态输水成效的遥感生态指数模型，分别获取不同年份的用于评估流域生态输水成效的遥感生态指数；

31、s32、根据步骤s32中获取的不同年份的用于评估流域生态输水成效的遥感生态指数，采用非参数化非参数化趋势度分析方法对不同年份的用于评估流域生态输水成效的遥感生态指数进行非正态分布与正态分析的序列数据进行拟合效果分析，得到生态环境质量变化趋势，即：

32、

33、其中，β表示生态环境质量变化趋势，median(·)表示统计函数，xj、xi分别表示第j、第i年的遥感生态指数，表示任意符；

34、s33、判断步骤s32中得到的生态环境质量变化趋势是否大于0，若是，则表明生态环境质量增加，即流域生态输水成效好，否则，表明生态环境质量减低，流域生态输水成效差，生态退化。

35、本发明具有以下有益效果：

36、本发明所提出的一种流域生态输水成效评估方法，利用历史实测水文数据对流域遥感影像数据进行数据校正，提高了数据精度；其次，利用校正的流域遥感影像数据，通过计算流域的湿度、绿度、干度以及热度指标，利用空间主成分分析方法对指标赋予权重，从而构建用于评估流域生态输水成效的遥感生态指数模型，能够较好地反应出研究区域的流域生态环境质量变化；同时，采用非参数化趋势度分析方法对不同年份的遥感生态指数进行拟合效果分析，提高了遥感生态指数的精度，能够综合与客观的评估流域生态输水成效，提高了评估效率和准确度。