基于地理信息系统的农业地质环境监测研究
——以江西省余干县为例
2016-09-08刘德元杨嘉琳孙新华岑世宏张卓然
刘德元，杨嘉琳，孙新华，岑世宏，张卓然
(1.河南省煤田地质局资源环境调查中心，河南郑州 450003；2.华北水利水电大学外国语学院，河南郑州 450046；3.郑州工业应用技术学院，河南郑州 451150；4.河南省煤田地质局物探测量队，河南郑州 450000)
基于地理信息系统的农业地质环境监测研究
——以江西省余干县为例
刘德元1，杨嘉琳2，孙新华1，岑世宏3，张卓然4
(1.河南省煤田地质局资源环境调查中心，河南郑州 450003；2.华北水利水电大学外国语学院，河南郑州 450046；3.郑州工业应用技术学院，河南郑州 451150；4.河南省煤田地质局物探测量队，河南郑州 450000)
[目的]运用地理信息系统监测与评价农业土壤环境状况。[方法]以余干县为例，将数据空间统计方法与模糊综合评价法融入农业地质环境监测方法中，运用ArcGIS 10.2软件分析了土壤养分分布情况。[结果]余干县表层95.32%的面积土壤为一类土壤(未被污染土壤)，4.68%的面积为二类土壤(轻度污染土壤)，土壤整体上缺乏速效钾。[结论]地理信息技术为农业地质环境监测与评价提供了新方法。
地理信息系统；空间统计；土壤有机质；农业地质环境
地理信息系统(GIS)是一门综合性学科，涉及内容广泛，是在计算机软件、硬件支持下对空间地理分布数据进行采集、存储、分析、运算、显示描述的技术系统[1]。GIS技术可借助于计算机工具将地理分析功能与地图这种视觉化效果集成在一起，得到数据规律的直观显示[2]。地理信息作为一种特殊信息，源于地理数据，地理数据表征环境中的重要要素数量、质量及分布规律。地理数据包括属性特征数据、空间位置特征数据以及时域特征数据，这3种特征数据构成了地理空间分析三大要素[3]。近年来，处理地理信息的技术日益完善，主要包括数字化技术、存储技术、空间分析技术、环境预测模拟计算、可视化技术[4]。现代农业发展依赖于农业地质背景与土壤营养及生物工程技术。农业地质可通过优化土壤常量组分与微量元素比例[5]，使其向绿色农业、精细农业、有机农业发展。笔者结合GIS技术中的空间分析技术及可视化技术，以江西省余干县土壤元素为调查对象，以区域地球化学勘查技术为主要手段，对地质条件进行监测，旨在为农业地质环境监测技术提供科学依据。
1 材料与方法
1.1样品采集对余干县境内土壤进行样本采集化验。表层土壤样品采集密度为每4 km2(0～20 cm)采集1个点，深层土壤采集密度为每16 km2(150～200 cm)采集1个点，共采集得到表层和深层土壤样品752件。土壤化学样品分析指标包括As、Cd、Cr 、Cu、Hg、N、Ni、P、Pb、S、Zn、MgO、CaO、K2O与pH，通过化学方法测得各样本土壤中的元素含量。表1 为各元素化学特性参照标准。
表1 各元素参照标准值
1.2数据的空间统计方法空间统计方法适合处理地理区域及离散化的统计数据，空间统计核心技术是空间数据分析，主要包括空间权重矩阵构建、空间相关性分析及数据检验[6]。通过空间数据统计分析方法，研究余干县土壤金属含量情况，便于全面了解余干县农业地质情况。空间自相关性反映区域单元与邻近区域单元同一属性值的相关程度，也反映2个观测数据之间的潜在依赖性。在地理空间中，正空间自相关是指位置紧密连接的观测单元倾向于相似(图1a)，负空间自相关是指位置相近的观测单元与位置较远观测单元进行比较(图1b)。
注：a为正空间自相关；b为负空间自相关。Note: a was the positive spatial autocorrelation; and b was negative spatial autocorrelation. 图1 空间自相关示意Fig.1 Spatial autocorrelation
GIS系统生成的拓扑关系提供了空间临近关系度量，权重矩阵建立时采取邻接标准或距离标准，如式(1)所示，n个位置空间邻接关系通常用1个二元对称空间矩阵Wn×n:
(1)
定义二进制空间权重矩阵如下：
(2)
空间权重矩阵选取二进制空间权重矩阵，以实际空间邻接关系为例说明空间矩阵定义，图2中空间排列情况为9个区域。
图2 空间排列邻接关系Fig.2 Diagram of spatial arrangement adjacency relation
1.3模糊综合评价对于空间数据分析，采用模糊综合评价法。模糊数学理论更具备优势，模糊集、隶属度及隶属函数是模糊数学三大概念[7]。模糊集是指边界模糊的集合，其各元素关联关系用隶属度表示，隶属度∈[0,1]，表示属于该模糊集，隶属度为1，表示完全不属于该模糊集。隶属函数表示元素与隶属度之间的解析函数，针对不同模糊集建立不同隶属函数。针对该研究的重金属元素含量建立的隶属函数如下式所示:
(3)
式中，yi为第i个因素隶属度；ui为样品i金属含量观测值；ci为标准值(表1)；ut为标准下限值。在模糊综合评价方法中，一般包涵多个评价指标，笔者选取4个评价指标，分别为土壤有机质、碱解氮、有效磷、速效钾。模糊评价模型有3种：加法模型、乘法模型以及综合模型，笔者选取加法模型，计算公式如下：
(4)
式中，W为综合评价值；ai为评价指标权重；pi为评价指标隶属度。根据GIS空间统计方法，对不同评价指标进行图层叠加，进而形成一个综合评价指标。
2 结果与分析
2.1土壤化学特性
2.1.1土壤pH。从图3、4可以看出，深层土壤与表层土壤pH具有相似的分布规律，深层土壤呈弱酸性，表层土壤呈酸性，由表及深酸性逐步减弱，并没有突变等异常现象，符合自然界规律，说明余干县土壤pH整体稳定，并未受人类活动影响。在区域分布上，余干县北部深层土壤呈中性(pH6.4～7.5)，社庚—九龙—大溪—白马桥—古掉—东塘—区域呈酸性(pH4.4～5.5)，其余地方呈弱酸性(pH5.4～6.5)。
图3 表层土壤pH分布情况Fig.3 Distribution of pH in surface soil
图4 深层土壤pH分布情况Fig.4 Distribution of pH in deep soil
2.1.2土壤全氮。查阅资料可知，余干县深层土壤氮含量背景值为586.0 mg/kg，标准离差为282，表层土壤背景含量为1 369.6 mg/kg，标准离差为497.2。从图5、6可以看出，余干县表层至深层土壤全氮出现了强烈的富集现象，这是由于人类补给造成的结果。在空间分布上，表层与深层土壤全氮含量具备不同区域化学分布特性，这是由于深层土壤全氮含量受土母岩成岩环境制约，表层受到人为等外部干扰较大。从地域来看，余干县深层土壤全氮含量最高值出现在康山垦殖场与黄金埠，表层土壤全氮含量最高值出现在古埠镇，向外全氮含量逐步降低。
图5 表层土壤全氮分布情况Fig.5 Distribution of total nitrogen in surface soil
图6 深层土壤全氮分布情况Fig.6 Distribution of total nitrogen in deep soil
2.2土壤环境质量评价该研究根据相关规范《土壤环境质量标准》选取表层土壤中的Cd、Re、As、Cu、Pb、Cr、Zn、Ni微量元素作为评价指标，对表层土壤环境进行综合评价，运用内梅罗综合污染指数评价法[8]，计算土壤污染指数：
(5)
式中，P为土壤污染综合指数；Pi为土壤i元素污染指数；n为土壤参与评价的污染元素种类；(Pi,max)2为元素污染指数最大值的平方。
《土壤环境质量标准》中规定P<1，为未受污染；1
图7 余干县表层土壤环境质量Fig.7 Comprehensive evaluation of surface soil environmental quality in Yugan County
根据表层土壤有机质、N、P、K数据分析，应用ArcGIS中的Summrize功能统计分析土壤养分丰缺情况。
由表2和图8～11可知土壤有机质、碱解氮、有效磷、速效钾含量及区域分布情况。①有机质。余干县区内土壤有机质含量总体较高，88.18%的土壤为适度-富足有机质土壤，主要分布在古埠镇、玉亭镇、瑞洪镇、三塘乡等地区，鸳鸯港乡土壤有机质相对缺乏。②碱解氮。区域总面积的 97.28%为适度-富足碱解氮地区，主要分布在玉亭镇、瑞洪镇地区，碱解氮缺乏地区为鸳鸯港乡。③有效磷。区域总面积的 87.84%为适度-富足有效磷地区，主要分布在康山乡、瑞洪镇、信丰垦殖场，洪家嘴乡土壤有效磷含量偏低。④速效钾。仅区域总面积的 25.87%为适度-富足速效钾地区，主要分布在瑞洪镇地区，大多数地区表现为缺乏现状，杨掉乡土壤速效钾含量最低，占区域总面积的2.85%。
表2 余干县土壤养分丰缺状况
图8 土壤有机质含量分布情况Fig.8 Distribution of organic matter content in soil
图9 土壤碱解氮含量分布情况Fig.9 Distribution of the content of alkali-hydrolyzable nitrogen in soil
图10 土壤有效磷含量分布情况Fig.10 Distribution of available phosphorus content in soil
图11 土壤速效钾含量分布情况Fig.11 Distribution of the content of rapidly available potassium in soil
3 结论
(1)将数据空间统计方法与模糊综合评价引入农业地质环境监测方法中，结合ArcGIS10.2软件得出了各土壤养分分布情况，为农业地质环境监测提供了新方法。
(2)监测结果表明，余干县表层土壤总体环境质量优良，一类土壤(未污染土壤)分布面积为2 227.82 km2，占总面积的95.32%，二类土壤(轻度污染土壤)分布面积为109.28 km2，占总面积的4.68%。对土壤有机质、碱解氮、有效磷、速效钾含量以及区域分布情况进行了研究，结果表明余干县农业地质环境整体上缺乏速效钾。
[1] 甘洪川. 基于ArcGIS的数字校园地理信息系统的设计与实现[D].成都：电子科技大学,2010.
[2] 周成虎. 全空间地理信息系统展望[J].地理科学进展,2015(2):129-131.
[3] 杰克·丹哲芒,马劲武.地理信息系统:设计未来[J].中国园林,2010(4):19-26.
[4] 胡祎.地理信息系统(GIS)发展史及前景展望[D].北京：中国地质大学(北京),2011.
[5] 牟琳.地理信息系统在高中地理合作学习中的实践与探索[D].成都：四川师范大学,2012.
[6] 任忠慧.基于地理信息系统(GIS)的我国都市农业环境问题研究[D].武汉：华中师范大学,2014.
[7] 闫丽娟.校园三维地理信息系统的建模与可视化研究[D].大连：大连理工大学,2013.
[8] MAKANGA P T,SCHUURMAN N,VON DADELSZEN P,et al.A scoping review of geographic information systems in maternal health[J].International journal of gynecology and obstetrics,2016,134(1)：13-17.
Agricultural Geological Environment Monitoring Based on Geographic Information System—A Case of Yugan County in Jiangxi Province
LIU De-yuan1,YANG Jia-lin2, SUN Xin-hua1et al
(1.Resources and Environment Investigation center of Henan Provincial Coalfield Geology Bureau，Zhengzhou，Henan 450003; 2.School of Foreign Studies,North China University of Water Resources and Electric Power, Zhengzhou, Henan 450046)
[Objective] To monitor and evaluate the agricultural soil environment status by Geographic Information System. [Method] With Yugan County as a case, spatial statistic method and fuzzy comprehensive evaluation method were integrated into the agricultural geological environment monitoring method. ArcGIS10.2 software was used to analyze the nutrient distribution in soil. [Result] 95.32% surface area of soil in Yugan County belonged to first-grade soil (not polluted); 4.68% area was second-grade soil (slightly polluted). Soil was lack of rapidly available potassium in Yugan County. [Conclusion] Geographic Information System provides new method for the monitoring and evaluation of agricultural geological environment.
Geographic Information System; Spatial statistics; Soil organic matter; Agricultural geological environment
刘德元(1964- )，男，河北魏县人，高级工程师，从事地质工程研究。
2016-06-08
S 181.3
A
0517-6611(2016)20-056-04
