基于ArcGIS与分形理论确定土壤重金属异常下限值——以喀斯特小流域Cd为例

1、.基于ArcGIS与分形理论确定土壤重金属异常下限值以喀斯特小流域Cd为例-农学论文基于ArcGIS与分形理论确定土壤重金属异常下限值以喀斯特小流域Cd为例 宋书巧，胡 伟（广西师范学院地理科学与规划学院，南宁 530001）摘要：土壤重金属元素异常下限值的确定对区域土壤环境评价具有重要意义，本研究基于ArcGIS技术与分形理论对喀斯特小流域土壤重金属Cd元素的异常下限值进行研究，通过含量面积方法确定了喀斯特小流域土壤重金属Cd元素异常下限值为107 mg/kg。通过与传统方法（平均值加两倍标准离差）和85累计频率法进行对比，表明分形方法确定的异常区域是有效的、合理的。关键词：分形；Cd；异常

2、下限；土壤；喀斯特；流域中图分类号：X825 文献标识码：A 文章编号：04398114（2015）061361-04DOI:10.14088/ki.issn0439-8114.2015.06.020Determing Minimum Thresholds of Heavy Metals in Soil Based on ArcGIS and Fractal MethodTaking Cd in Karst Basin as an ExampleSONG Shuqiao， HU Wei（School of Geography and Planning， Guangxi Teachers Edu

3、cation University， Nanning 530001， China）Abstract： It is significant to determine the minimum thresholds of heavy metals in soil for regional environmental quality The minimum threshold of Cd in soil of karst basin was studied based on ArcGIS and fractal method The concentrationarea fractal method w

4、as used to determine the minimum threshold of Cd in soil of karst basin， which was 107 mg/kg Compared with the results obtained by traditional methods（x2S） and 85 of the cumulative frequency methods， the anomaly extension determined by fractal method was effective and reasonableKey Words： fracal； Cd

5、 element； minimum threshold； soil； Karst； basin收稿日期：20141209基金项目：2014广西研究生教育创新计划项目（YCSW2014184）；广西软科学研究计划课题（1121700231）作者简介：宋书巧（1964），女，山东烟台人，博士，硕士生导师，主要从事土壤环境质量评价的研究，（电话电子信?a class=“_cf_email_” href=“/cdn-cgi/l/email-protection” data-cfemail=“46a2e5ef35352e333706777075”>emailprotect

6、edcom；胡伟（1988），男，湖北仙桃人，在读硕士，研究方向为区域环境科学，（电话电子信?a class=“_cf_email_” href=“/cdn-cgi/l/email-protection” data-cfemail=“18fcbbb1707a706d6f7d71586b717679”>emailprotectedcom。分形理论是当今世界十分风靡和活跃的新理论、新学科，分形的概念是美籍数学家曼德布罗特（BBMandelbort）首先提出的1，分形理论目前已经成为定量描述地理和环境现象的有力工具。用分形理论刻画自然界中一些不规则、不稳定或具有高度复

7、杂结构的现象，能够以少量数据描述总体情况24。土壤具有不同程度的空间变异性，导致其分布的重金属元素也存在空间变异性。确定土壤重金属异常下限值是研究土壤重金属污染物空间分布特征的关键步骤之一。以往通过计算研究区域元素含量平均值与2倍标准离差之和作为元素的地球化学异常下限值，这种方法对土壤重金属元素异常下限的判断，往往不能准确地反映土壤中重金属元素异常状况57。地球化学数据是一个具有低维吸引子的混沌系统，背景和异常的形成是两个独立的过程8。利用分形方法研究区域土壤重金属含量的异常下限值是可行的。本研究以喀斯特流域为研究区域，计算流域内土壤重金属Cd含量异常下限值，并分析其合理性与有效性。1 材料与

8、方法11 研究区域概况流域地处亚热带，属南亚热带季风型气候，日照充足，气候温和，雨量充沛，年平均气温在213 ，极端最高温度389 ，最低温度07 。流域面积106 km2，最大流量为69 m3/s，正常流量为1 m3/s。以县城为界，河流分为上、中、下游，县城位于中游。流域内分布着3种地貌，即东部喀斯特峰丛洼地、中部喀斯特岩溶谷地、西南部中低山地区，其中喀斯特地区占整个流域面积的75以上。主要土壤有石灰土、水稻土、红壤，石灰土以棕色石灰土为主，石灰土分布占流域面积的75以上。12 样品采集与测定采用随机网格法取较为稳定地块的表层土壤（图1），每个采样单元内利用梅花布点法，采集020 cm表层

9、土的样品混合为一个样，共采集52个土样。采集后的土壤样品在室温下自然风干后研磨，过200目筛。处理好的土壤样品采用反王水和氢氟酸微波消解，利用ICPMS（Vista MPX， Varian， USA）测试样品Cd元素总量。所有样品均由广西地质矿产测试研究中心测试，土壤样品中Cd元素总量范围在0103300 mg/kg之间。 13 研究方法CA分形方法（浓度面积法）是分形理论中用于计算土壤重金属元素异常下限值的方法之一。从分形理论来看，区域土壤重金属异常时，异常值可能遵循分形分布而不是正态或对数正态分布。符合分形分布的特点要求，大于等于元素的某一含量值与该含量值的分布范围（面积）存在幂函数关系，

10、即9，10：N(a)a-D式中，a表示元素含量；D0，为分维数；N(a)N（a）表示含量大于等于a的分布面积。区域土壤重金属元素含量异常，在双对数坐标下绘制含量与面积的散点图时，其散点大致分布在两条直线上，采用分段拟和的方法，在两个区间用最小二乘法进行回归分析，并用最优化方法确定分界点。这里所获得的分界点的地球化学意义可以看成是元素含量在空间上至少存在两个层次的分布，即小于分界点的含量为元素含量的背景分布，大于分界点含量为元素含量的异常分布，该点对应的值为元素含量分布的异常下限值1116。CA分形法难点在于计算不同元素含量对应的面积，这点可以通过在ArcGIS平台进行等值线的描绘、面积与含量的

11、统计来解决，其具体过程如下：1）在ArcGIS软件平台中，利用普通克里金插值法将已经实测重金属元素含量数据进行插值处理并检验插值的正确性，然后生成重金属元素含量分布的栅格文件。 2）利用ArcGIS等值线功能，使用重金属元素含量栅格文件生成线文件，再将线文件与研究区范围叠加生成等值面。 3）对等值面进行条件检索，检索不同的a值（含量值），对大于等于所检索 C 值的等值线区域进行面积统计。 4）不断重复上述步骤，选取不同的a的值，统计其等值线区域面积。 5）对不同的a的值及其对应的面积数据取对数，进行最小二乘法拟和，其分界点对应的a的值即为元素的异常下限值。2 结果与分析21 CA分形方法对异常

12、下限值的确定运用浓度面积法对流域表层土壤中Cd元素的多维分形进行探讨，ArcGIS处理过的元素含量与面积在双对数坐标系统下的分布符合两段分布模型，用最小二乘法拟合为两段直线（图2），两段直线的方程分别是： y0128 5x12855（01x107）（1）y1252 4x13545（107x3300）（2）利用最优法求解两条直线在相应区间的剩余平方和（EE1E2）是107，以上两条直线的方程均通过一致性检验，方程1中R2为0903；方程2中R2为0908。两条直线拟合度高。分界点所对应Cd元素含量为107 mg/kg，由此得到该流域土壤中Cd元素的异常下限值为107 mg/kg。 22 传统法计

13、算下限值传统计算方法是建立在数据符合正态或对数正态分布基础上，但在实际工作中没有任何数据完全服从正态分布。所以对于测试数据一般要进行离群点数据（最高值、最低值）的迭代处理，一般认为不在X±3S区间内的值为异常值，应该剔除。处理后的数据平均值与两倍标准差值和为元素的异常下限。实测数据描述统计量见表1。传统方法计算的Cd元素异常下限值为1468 mg/kg。23 85累计频率法计算下限值累计频率就是一个数值的频率和比它的频率高的数值频率的总和。设x1x2xm不重复的样本值，mn。把样本值小于或等于样本的数据xi的频率累加，得到小于或等于xi的累积频率。再把累积频率85时xi的值作为异常下

14、限，该法计算的Cd元素异常下限值为702 mg/kg。3 讨论与结论31 讨论三种方法对Cd元素的异常下限值计算的对比，见表2。采用CA分形法，对流域土壤Cd元素异常的下限值进行计算，其异常值的下限为107 mg/kg ，其异常区域占整个流域面积的6875；采用传统方法计算流域土壤Cd元素异常的下限值为1468 mg/kg，其异常区域占整个流域面积的061；采用85累计频率法对流域土壤Cd元素异常的下限值进行计算，其异常值的下限为702 mg/kg，其异常区域占整个流域面积的1294。其分布范围如图3所示。采用传统方法计算的异常值，对数据的要求较高，需要数据大致服从正态分布，而流域内实测土壤中

15、Cd元素的含量基本不符合正态分布，从而导致该法的计算值远高于我国土壤环境质量标准（B 156181995）中规定土壤中Cd元素的三级标准10 mg/kg。采用85累计频率法，同样受到实测数据的影响，从而导致其计算结果也超过土壤中Cd元素的的三级标准，以上两种方法在未受到污染的区域使用，其结果明显偏大，而采用CA分形的方法计算流域内土壤中Cd元素异常的下限，考虑了土壤中地球化学元素的空间相关性，将元素背景值与异常看作两个分开的过程。并且其异常区域分布范围主要位于流域的中游（县城），说明流域中游的人为干扰是土壤Cd元素含量偏高原因之一，此外喀斯特地区中Cd元素本底值偏高17，是该流域土壤Cd元素异

16、常下限值略超过土壤环境三级标准的原因之一。以上都能够说明利用CA法计算流域土壤Cd元素异常下限值是可行的，是符合实际的。32 结论流域土壤重金属Cd元素含量与面积的对数关系符合两段分布模型，流域内土壤重金属Cd元素的分布具有分形特征，基于ArcGIS平台，利用CA分形方法计算流域内土壤重金属Cd元素异常下限值可行，其异常值为107 mg/kg。通过对比分析，基于ArcGIS和分形方法确定流域内土壤重金属Cd元素异常下限值是有效的、合理的；相较传统的方法和85累计频率法，其计算的结果具有更高的可信度。 参考文献： 1 MANDELBORT B B How long is the coast of

17、 Britain statistical selfsimilarity and fractional dimensionJScience，1967，156（3775）：6366382 朱晓华，王 建分形理论在地理学中的应用现状和前景展望J大自然探索，1999（3）：42463 周广柱，肖唐付，刘意章，等多重分形理论在环境科学领域的研究进展J矿物岩石地球化学通报，2013，32（1）：1071134 柯贤忠，程绪江，谢淑云，等基于多重分形的表层土壤中重金属元素来源分析以海南省琼海市为例J地质科技情报，2010，29（1）：971025 施俊法，向运川，王春宁区域地球化学异常空间分形结构及其意义以

18、浙江省诸暨地区区域地球化学数据为例J矿物学报，2000，20（1）：68726 李建东，龚庆杰，窦金龙，等分形理论在勘查地球化学中的应用J中国矿业，2006，15（11）：91947 熊 超， 葛良全， 罗耀耀， 等 多种确定地球化学异常下限方法的比较J 四川有色金属，2012（1）：5255，618 王瑞廷， 毛景文， 任小华， 等区域地球化学异常评价的现状及其存在的问题J中国地质，2005，32（1）：1681759 谢淑云，鲍征宇，秦 兵，等成都盆地浅层土壤中元素来源研究J土壤通报，2007，38（5）：89690010 袁 峰，张颖慧，周涛发，等基于分形方法的合肥大兴地区土壤中污染元素Hg异常下限J生态环境，2008，17（6）：2259226311 张 焱，成秋明，周永章，等分形插值在地球化学数据中的应用J中山大学学报（自然科学版），2011，50（1）：13313712 张法升，刘作新分形理论及其在土壤空间变异研究中的应用J应用生态学报，2011，22（5）：1351135813 疏志明，王雄军，赖健清，等分形理论在太原盆地土壤重金属元素分析中的应用J物探与化探，2009，（2）：1