[发明专利]一种不同土地功能区域的生态环境受损动态评估方法

权利要求书

1.一种不同土地功能区域的生态环境受损动态评估方法，其特征在于，包括以下步骤：

S1、对研究区域进行单元栅格划分，对划分后的不同土地功能区域的生态环境系统，建立矿区生态环境受损评估动态模型；

S2、根据研究区域的土地类型，采用AHP法建立不同土地类型对地表形变指数的敏感性评价体系，得到各个指标权重；

S3、通过各个指标权重对矿区生态环境受损评估动态模型进行赋值，采用ARCGIS对赋值后的矿区生态环境受损评估动态模型进行赋值求解，得到不同土地功能区域的生态环境受损值。

2.根据权利要求1所述的不同土地功能区域的生态环境受损动态评估方法，其特征在于，所述步骤S1包括以下分步骤：

S11、对研究区域进行单元栅格划分，得到多个栅格单元；

S12、采用差分干涉测量方法D-InSAR采集每个栅格单元的地表形变值，并对地表形变值进行无量纲化处理，得到矿区地表形变指数j；

S13、通过ArcGIS处理遥感数据，得到植被覆盖指数值，并基于植被覆盖指数值，得到矿区地表植被受损指数k；

S14、根据待建立的矿区生态环境受损评估动态模型的最大值，得到常系数c；

S15、根据矿区地表形变指数j、矿区地表植被受损指数k和常系数c，对划分后的不同土地功能区域的生态环境系统，建立矿区生态环境受损评估动态模型。

3.根据权利要求2所述的不同土地功能区域的生态环境受损动态评估方法，其特征在于，所述步骤S12中计算矿区地表形变指数的公式为：

其中，j为矿区地表形变指数，D为每个栅格单元的地表形变值，D_max为整体最大沉降值。

4.根据权利要求2所述的不同土地功能区域的生态环境受损动态评估方法，其特征在于，所述步骤S13中计算矿区地表植被受损指数的公式为：

k＝1-NDVI

其中，k为矿区地表植被受损指数，NDVI为不同栅格单元的植被覆盖指数值，NDVI的值域为[0，1]。

5.根据权利要求2所述的不同土地功能区域的生态环境受损动态评估方法，其特征在于，所述步骤S14中常系数c的计算公式为：

c＝M+N

其中，M为矿区地表形变指数的系数，N为矿区地表植被受损指数的系数。

6.根据权利要求2所述的不同土地功能区域的生态环境受损动态评估方法，其特征在于，所述步骤S15中矿区生态环境受损评估动态模型为：

其中，y为矿区生态环境受损评估动态模型，M为矿区地表形变指数的系数，N为矿区地表植被受损指数的系数，D为每个栅格单元的地表形变值，D_max为整体最大沉降值，NDVI为不同栅格单元的植被覆盖指数值。

7.根据权利要求1所述的不同土地功能区域的生态环境受损动态评估方法，其特征在于，所述步骤S2中研究区域的土地类型包括：林草地、耕地、水域、裸地、建设用地和采煤用地。

8.根据权利要求1所述的不同土地功能区域的生态环境受损动态评估方法，其特征在于，所述步骤S3中赋值后的矿区生态环境受损评估动态模型为：

其中，y₁为林草地的生态环境受损值，y₂为耕地的生态环境受损值，y₃为水域的生态环境受损值，y₄为裸地的生态环境受损值，y₅为建设用地的生态环境受损值，y₆为采煤用地的生态环境受损值，{w₁，w₂，w₃，w₄，w₅，w₆}为各个指标，{W₁，W₂，W₃，W₄，W₅，W₆}为各个指标权重，D为每个栅格单元的地表形变值，D_max为整体最大沉降值，NDVI为不同栅格单元的植被覆盖指数值。