近年来，随着工业化、城市化和农业集约化的快速发展，我国土壤环境所面临的情况不容乐观，其中重金属污染问题尤为突出[1-5]。土壤重金属作为土地生态系统中的慢变量，具有隐蔽性、积累性、地域性且不易降解、易富集等特点[6-8]。工厂废水、废气，机动车尾气的排放以及农业活动化肥、农药的施用，更加剧了土壤重金属污染的局面。据2014年4月环保部和国土资源部联合发布的《全国土壤污染状况调查公报》报道，我国土壤无机污染物（重金属）超标点位数占全部超标点位的82.8%，其中Cr的点位超标率为1.1%[9]。土壤中的Cr主要来源于大气中重金属Cr的沉降，农药、化肥和塑料薄膜的使用以及污水灌溉等[10]。Cr作为五毒元素之一，它不仅能抑制土壤微生物的活性，降低土壤中的生物多样性，还会通过饮用水、食物链等方式对生物体以及生态环境造成严重影响[11-14]。因此，有必要对区域尺度上土壤重金属Cr的污染分布特征及其影响因素进行系统研究。

整个夜晚，甲洛洛都同情着丁主任，感动着丁主任对自己说的话，觉得自己的份量一下加重了，心口也热烘烘的，他想着老婆本玛对自己敬重的眼神，不由得呵呵地笑着：甲洛洛，你可真是个老好人啊！

围绕农田土壤重金属分布特征，很多学者已开展大量研究[15-16]，张小敏等[17]探讨了我国农田土壤重金属 Pb、Cd、Cu、Zn和 Cr的空间富集现状，表明我国重金属空间分布具有明显的区域特征；王彬武等[18]对耕地土壤重金属的时空变化进行研究，运用潜在生态风险和指示克里格法对比分析不同年份耕地表层土壤重金属含量；韩平等[19]对北京顺义区土壤重金属污染生态风险进行评估研究；索琳娜等[20]研究表明，北京市菜地重金属 Cd、Cr、Cu、Zn 存在累积污染风险，但总体土壤环境安全；王学军等[21]应用克里格方法对北京东郊污灌土壤表层重金属含量进行了分析，采用单因子和多因子评价，结果表明研究区土壤属于轻度污染或无污染；赵倩等[22]运用多元统计分析和潜在风险指数方法，研究土层中重金属的分布特征与潜在生态风险。鉴于目前对于北京市土壤重金属污染研究大多集中于水平空间、不同土壤类型以及不同功能区方面，本文在已有研究的基础上，以北京平原区农田土壤为研究对象，基于所监测数据及土壤基本理化性质和相关环境因子，研究不同土壤深度Cr水平和垂直分布特征，并从不同功能区、不同种植模式，以及不同母质类型等方面探究该区域土壤重金属Cr的影响因素，旨在为今后该区域人为干预以及重金属Cr预警模型建立和污染治理提供科学参考。

1 材料与方法

1.1 研究区概况

北京市位于 115.7°～117.4°E、39.4°～41.6°N，地处华北平原北部，背靠燕山，面积达1.641万km2，其中山区面积约为10 068 km2，约占全市总面积的61.2%。地势西北高、东南低，平均海拔为43.5 m。北京为典型的北温带大陆性半湿润季风气候，夏季高温多雨，冬季寒冷干燥，春、秋短促。年平均温度8～12℃，降水量600～700 mm。全年无霜期180～200 d，西部山区较短。

首先，在设计建筑的外观时，相关的设计人员可以利用色彩的组合营造出新的建筑风格，使整个建筑展现出不同的风格，利用色彩的搭配可以有更好的艺术效果，转变传统的陈旧风格和束缚，来使整个建筑变的主次分明，增加整个建筑的艺术魅力。除此之外，相关的设计人员可以采用两种不同颜色叠加搭配的方式来增强整个建筑的视觉效果，因为色彩可以带来视觉上的享受和冲击，让人们可以直观的感受到建筑物的美观性。

由表3可见，0～25 cm及25～50 cm土壤中不同种植模式下Cr含量均值较为接近，且3种不同种植模式下重金属Cr的含量均已超过北京市土壤重金属Cr环境背景值。果园的土壤Cr含量在0～25 cm土壤深度中高于耕地，25～50 cm土壤深度中含量为最高。原因在于一些农业污染源如肥料、农药的施用，鸡粪是主要有机肥，鸡粪中Cr的含量达到20.5 mg·kg-1[33]，磷肥中含有一定量的Cr元素[34]，相比耕地和菜地，果园需肥量较大，因此造成该区域果园重金属Cr含量偏高。此外，我国土壤有机质普遍较低，土壤磷素也普遍偏少，因此肥料的投入增加是导致研究区Cr累积的重要原因。河水和灌溉用水中Cr的沉淀易被土壤吸附，污水灌溉与该区域Cr含量富集程度较高也具有一定的关联性。

1.2 数据的来源及处理

本文以北京市农业用地为研究对象，采样时间为2015年6月，针对不同种植模式、不同土壤深度（0～25、25～50 cm）选取94个采样点。每个样点选取中心点及相邻两角，并将同一深度土样混合作为代表该点的样品，样品均经野外剔除与金属采样器械接触部分。94个采样点覆盖北京平原区各区县，共分布于九区一县，自北向南依次包括水源保护区的密云、平谷、怀柔、昌平，农业保障区的顺义、通州、大兴和房山，以及都市生活区的海淀、丰台。采样点分布位置见图1。利用GPS定位采样点，并记录采样点的功能区属性、种植模式、土壤类型和土质等基本属性，检测了土壤样品中重金属Cr的含量。土壤中重金属Cr采用火焰原子吸收分光光度法进行测定，分析过程采用平行双样、加标进行质量控制。

  

图1 研究区位置及采样点位置分布图Figure 1 Location of the study area and distribution of sampling points

采用SPSS 21.0进行描述性统计分析、方差分析，计算变异系数，并结合ArcGIS 10.2软件对样点数据进行插值处理分析，预测其空间分布特征。

1.3 研究方法

1.3.1 单因子指数法

采用单因子指数法对土壤中不同深度的重金属Cr污染强度进行评价，评价标准见表1，计算公式为：

不同土壤深度重金属Cr的单因子污染指数Pi在不同区域的数值存在差异。0～25 cm土壤深度，平谷区 Cr污染指数 Pi（ave）最低为 0.66，说明水源保护区的平谷土壤处于清洁水平；房山、昌平两区Pi(ave)分别为0.984和0.848，属尚清洁水平，但由于其处于警戒线位置，应予以重视；其余各区县均超过警戒线。25～50 cm土壤中，除昌平和平谷处于尚清洁状态，其余各区县均达到轻度污染水平，其中以海淀区重金属Cr的 Pi(ave)最高，为 1.665。在 0～25 cm 及 25～50 cm 土壤中，受污染的样点数分别占总样点数的73.4%和77.7%，整体比例较高。

 

该研究区主要种植模式为果园、耕地、菜地3种类型，3种种植模式下不同土壤深度（0～25、25～50 cm）重金属Cr含量数值统计见表3。通过方差齐性检验及多重比较分析，结果显示不同种植模式组间无显著性差异。

 

表1 土壤评价分级标准[23]Table 1 Evaluation of soil classification standard[23]

  

等级 环境质量指数 污染程度 污染水平1 Pi≤0.7 安全 清洁2 0.7＜Pi≤1 警戒线 尚清洁3 1＜Pi≤2 轻度污染 土壤轻度污染4 2＜Pi≤3 中度污染 土壤受到污染5 Pi＞3 重度污染 土壤受污染已较严重

1.3.2 经验贝叶斯克里格插值

经验贝叶斯克里格法（EBK）是一种常用的地统计插值方法，目前常用的克里格方法是通过已知的数据位置计算半变异函数，并使用此单一变异函数在未知位置进行预测，具有一定的平滑效应，不能够完全反映出土壤本身的性质特点[24-26]。本文所采用的经验贝叶斯克里格法与其他克里格方法有所不同，它通过估计基础半变异函数来说明所引入的误差，可准确预测一般程度上不稳定的数据，预测标准误差比其他克里格方法准确，尤其对于小型数据集更准确。因此，针对研究区的实际情况，该方法更适用于土壤重金属Cr污染空间分布的研究。

受动载影响，左帮片冒严重，因此左帮锚杆在安装后，1、3号受力逐渐减小，甚至趋于0，基本失效(图10)。2号未失效，受力约15 kN。这与数值分析中左帮锚杆受力较小相吻合。

为了了解参保人员的基本情况，笔者随机走访了几个村进行实地调研，发放调查问卷120份，调查对象为男性70人，女性50人，年龄段主要集中在31-59岁，占调查总人数的91.7%；60岁以上的10人。受教育情况为，小学及以下、初中和高中的人数分别为61人、47人和12人。经济来源情况：纯务农的63人，这部分人员月均收入均在1000元以下；手工艺者27人，月均收入约在2000—4000元，经商者12人，依靠子女供养的18人。

2 结果与讨论

2.1 土壤重金属Cr描述性统计分析

利用 SPSS 21.0对各区县不同深度（0～25、25～50 cm）土壤重金属Cr含量特征进行描述性统计分析，结果见表2。

 

表2 土壤重金属Cr含量描述性统计Table 2 Descriptive statistics of soil heavy metal Cr

  

注：北京市土壤重金属Cr环境背景值为29.8 mg·kg-1[27]。超标率为各区县超标样点数与各区县总样点数的比值。Note：The environmental background value of soil heavy metal Cr in Beijing is 29.8 mg·kg-1[27].The rate of exceeding standard is the ratio of the standard sample points of each county and the total sample number of each county.

 

地区 计算样本数 Cr含量（0～25 cm） Cr含量（25～50 cm）均值/mg·kg-1 变异系数 超标率/% Pi（ave） 均值/mg·kg-1 变异系数 超标率/% Pi（ave）顺义 17 34.159±10.460 0.306 58.82 1.146 34.424±11.632 0.338 52.94 1.155房山 4 29.325±10.138 0.346 25.00 0.984 32.525±10.748 0.330 50.00 1.091昌平 9 25.278±13.667 0.541 22.22 0.848 27.800±15.061 0.542 33.33 0.933大兴 9 37.667±17.862 0.474 66.67 1.264 32.029±12.200 0.381 55.56 1.075怀柔 10 34.800±10.850 0.312 60.00 1.168 32.230±8.174 0.254 60.00 1.082海淀 4 57.725±18.055 0.313 100.00 1.937 49.625±10.471 0.211 100.00 1.665密云 13 41.508±23.154 0.558 76.92 1.393 45.046±17.604 0.391 92.31 1.512平谷 12 19.679±12.351 0.628 25.00 0.660 23.737±28.027 1.181 16.67 0.797丰台 6 47.217±6.298 0.133 100.00 1.585 44.500±7.865 0.177 100.00 1.493通州 10 36.810±9.675 0.263 70.00 1.235 35.980±10.327 0.287 60 1.207

本研究将北京市平原区农用地划分为3种不同的功能分区：包括水源保护区的密云、平谷、怀柔、昌平，农业保障区的顺义、通州、大兴和房山，以及都市生活区的海淀、丰台，对0～25、25～50 cm不同土壤深度的重金属Cr含量的分布特点进行比较分析。

激荡40年，每一位被采访者，讲述过往，满满正能量。看的出来，距离与对比产生幸福感。当我们求教人生经验与感悟时，谈的最多的是：踏实态度。也许全球化的当下，太多的选择性让我们显得不踏实。其次的告诫：学习能力。适应性强主要在于学习性强，而自学力尤为重要。反复地强调：合作精神。今天单打独斗不可能成功，靠合作、靠团队，才可撑起一片天。

2.2 土壤重金属Cr含量空间分布特征

利用ArcGIS软件中的地统计分析工具，采用EBK法得到不同土壤深度重金属Cr的空间分布图（图2），整体了解不同土壤深度Cr随区域变化的空间分布趋势。EBK子集大小为100，重叠因子为1，模拟次数为100次。邻域类型为半滑圆形，平滑系数为0.2。

图2a、图2b分别为研究区0～25 cm和25～50 cm深度土壤重金属Cr含量分布图，其含量范围分别为16.525～59.551 mg·kg-1，12.055～84.849 mg·kg-1。整体上，0～25 cm和25～50 cm深度土壤Cr水平分布格局基本一致，均为都市生活区的海淀、丰台含量居高，并以城区周围为中心向外逐渐降低，可见人类活动导致了部分土壤中Cr的积累，对土壤Cr含量造成一定影响[29]。密云北部在不同土壤深度下重金属Cr含量分布显著高于周边地区，这主要与该地区有一个Cr矿有关。25～50 cm深度下，密云南部、顺义东部及平谷西部位置有明显的Cr富集现象，由于该区域土壤环境质量与北京市生活用水的水质有较大关系，因此Cr含量过高是否会对水质产生影响则值得关注。插值分析与表2所述结果基本一致，预测结果与实测的拟合程度较高。

834 Research status and application prospect of artificial intelligence technology in lung tumors

为检验垂直方向重金属Cr的含量间相似性和来源差异，对0～25 cm及25～50 cm土壤深度的重金属Cr含量进行相关性分析（相关系数越接近1或-1，相关性越强；相关系数越接近0，相关程度则越弱[30]）。结果显示0～25 cm和25～50 cm土壤深度重金属Cr含量呈显著正相关，相关系数为0.579。因此，北京市平原区农田土壤垂直方向Cr的含量存在相关性。

  

图2 不同土壤深度Cr空间分布图Figure 2 Different soil depths spatial distribution pattern of Cr

2.3 土壤重金属Cr污染程度评价

采用ArcGIS中常用的地统计分析工具EBK对研究区不同土壤深度重金属Cr的环境质量指数Pi插值，并进行等级划分（图3），统计出不同等级重金属Cr所占面积，从而判断研究区不同土壤深度整体Cr的富集程度。EBK子集大小为100，重叠系数为1，模拟次数为300，领域类型为平滑圆形，平滑系数为0.2。

  

图3 不同土壤深度Cr污染等级分布图Figure 3 The distribution of Cr pollution grade in different soil depths

图3为北京平原区不同土壤深度Cr污染等级分布图。结合表1，共分为4个等级，从研究区整体看，不同土壤深度重金属Cr污染指数分布格局大体一致，其中处于警戒线以内的区域占研究区总面积的19.16%（0～25 cm）和 14.81%（25～50 cm），轻度污染的地区整体比例明显高于其他等级，分别占研究区总面积的 80.35%（0～25 cm）和 84.84%（25～50 cm）；处于中度污染及以上等级的区域面积仅占总体的0.03%（0～25 cm）和 0.36%（25～50 cm）；保持在清洁水平的区域基本位于水源保护区的平谷区内。分析垂直方向不同土壤深度Cr污染等级分布可见，除清洁水平外，其余各等级25～50 cm土壤深度重金属Cr的富集程度均高于0～25 cm，其中0～25 cm土壤深度中都市生活区的海淀存在2.52 km2区域处于中度污染，25～50 cm深度密云北部及其与顺义、平谷交界处均出现了Cr的富集，达到中度污染等级。研究认为该区域土壤Cr含量受母质影响较大[31]，与郑袁明等[28]的研究结果一致。

以2017年7月为分水岭实施风险管理，随机选取2017年1月～7月在我院接受治疗的骨科患者136例作为对照组，另选取2017年8月～2018年9月我院骨科收治的患者136例作为对策组。其中，对照组男81例，女55例，年龄18～54岁，平均年龄（32.1±2.4）岁；对策组男79例，女57例，年龄22～46岁，平均年龄（31.6±1.7）岁。两组患者性别、年龄等一般资料比较，差异均无统计学意义（P＞0.05）。参与护理人员前后一致，4名护士，1名护师，2名护理主管，具有较好的沟通能力以及自愿参与试验。

2.4 土壤重金属Cr含量影响因素分析

2.4.1 不同功能区土壤重金属Cr含量分析

由表2可知，大兴、怀柔、海淀、丰台、通州5个区域0～25 cm深度的土壤Cr含量均高于25～50 cm深度，顺义、房山、昌平、密云、平谷5个区Cr含量随土壤深度的增加而增加。除平谷和昌平，其他各地区土壤Cr含量均已超过北京市土壤重金属Cr环境背景值（29.8 mg·kg-1）。0～25 cm 深度，海淀区和丰台区 Cr含量最高，分别为（57.725±18.055）mg·kg-1和（47.217±6.298）mg·kg-1；25～50 cm 土壤中，海淀区 Cr含量为（49.625±10.471）mg·kg-1，位居首位，海淀和丰台两区的点位超标率均达100%。平谷25～50 cm土壤Cr变异系数为1.181，属强变异程度。为进一步分析不同区县土壤Cr含量是否存在显著差异，分别对不同区县0～25 cm及25～50 cm土壤深度的数据进行方差齐性检验及多重比较，检验结果表明，0～25 cm土壤中丰台与昌平、平谷两地，以及平谷和通州地区Cr含量存在显著性差异（P＜0.05），丰台位于市区周边，在区域农业生产及人类活动频繁的干扰下，该地区土壤Cr的聚集趋势较为明显[22]，通州处在研究区东南部，该地区主要是冲积平原，其Cr含量偏高与冲积过程有一定关系[28]，水源保护区的昌平与平谷受人类活动扰动相对较少，因此上述不同区县间土壤Cr含量存在差异；25～50 cm土壤深度各区县Cr含量无显著差异。垂直方向Cr分布不同于以往随深度增加而降低的规律，表现出人为和自然地质的双重叠加影响，与赵倩等[22]的研究结果一致。

  

图4 不同功能区土壤重金属Cr含量平均值Figure 4 The average of soil heavy metal Cr content in different functional zones

 

同一土壤深度下不同字母表示在0.05水平具有显著性差异Different letters at the same soil depth indicate significant differences at the 0.05 level

图4为北京市平原区农用地不同功能区土壤Cr含量平均值。如图4所示，都市生活区和农业保障区中0～25 cm土壤重金属Cr含量均高于25～50 cm。Cr含量在不同土壤深度中排序均为都市生活区＞农业保障区＞水源保护区，且0～25 cm土壤中数值变化更为明显。影响土壤重金属含量与空间分布趋势的因素主要有两大类：一是母质及土壤自身的性质；二是人为活动的影响[32]。水源保护区受人类活动扰动较少，因此土壤Cr含量极可能受其母质影响，pH值、黏粒含量等会影响土壤对Cr吸附以及重金属本身的迁移，从而造成土壤中Cr含量发生变化。相比较水源保护区与农业保障区，都市生活区受人类活动影响较大，进而导致了部分土壤中Cr的积累。为进一步分析不同功能区土壤Cr含量是否有显著差异，分别对0～25 cm及25～50cm土壤深度的数据进行方差齐性检验，检验结果表明，0～25 cm数据满足方差齐性，可采用LSD法进行多重比较，而25～50 cm数据不满足方差齐性，因此采用Games-Howell（A）法进行多重比较，检验结果显示组间存在差异（图4）。通过采用单因素方差分析，结果显示不同功能区对不同土壤深度重金属Cr的含量均存在显著性影响（P＜0.05）。

2.4.2 不同种植模式下土壤重金属Cr含量分析

式中：Pi为土壤中污染物i的环境质量指数；Ci为污染物 i的实测浓度，mg·kg-1；Si为污染物 i的评价标准，mg·kg-1。

 

表3 不同种植模式下土壤重金属Cr含量平均值（mg·kg-1）Table 3 The average of soil heavy metal Cr content under different planting patterns（mg·kg-1）

  

土壤深度 果园（n=27） 耕地（n=36） 菜地（n=31）0～25 cm 35.792±22.398 33.277±11.539 36.448±15.332 25～50 cm 35.494±21.836 34.903±14.192 34.295±13.680

根据北京市数字高程，结合农用地粮田、菜地（露天菜地、设施菜地）和果园的空间分布情况，以100 m等高线为界线，将高程≤100 m的区域作为本文的研究区。研究区位于北京东南部，面积为7 779.12 km2。研究区的土壤质地分为轻壤、砂壤、砂土、中壤、重壤5种，其中以轻壤质为主，占到研究区总面积将近50%，其次是砂壤质和中壤质。根据北京市数字土壤系统，研究区可分为农业保障区、都市生活区、水源保护区3类功能分区，不同功能区土壤类型均以潮土、褐土为主，还存在有少量水稻土。主要成土母质包括冲积物、洪冲积物、钙质岩类风化物、黄土物质、红黄土性物质、页片岩类风化物等。

2.4.3 不同母质类型下土壤重金属Cr含量分析

母质是形成土壤的物质基础，是土壤的骨架和矿物质的来源。其机械组成影响土壤的机械组成，化学成分对土壤形成、性质和肥力均有显著影响[35]。研究区不同区县其土壤母质类型有所差异，大兴区母质类型主要为冲积物，顺义、昌平、通州、海淀和丰台地区均以冲积物和洪冲积物为主，其中海淀区存在少量的红黄土性物质，丰台地区存在部分钙质盐类风化物，密云和平谷地区母质类型较为复杂，密云地区为红黄土性物质、冲积物、洪冲积物、水库、长石岩类风化物，平谷地区则以钙质岩类风化物、洪冲积物、冲积物为主要母质成分。结合北京土壤环境背景特点分析不同母质类型下（主要分为冲积物、钙质岩类风化物、红黄土性物质、洪冲积物、黄土物质和其他6种类型）土壤重金属Cr的含量变化情况，见图5。

  

图5 不同母质类型下土壤重金属Cr含量平均值Figure 5 The average of soil heavy metal Cr content in different parent material types

由图5可知，冲积物、红黄土性物质、洪冲积物、黄土物质和其他这5种母质类型下在不同土壤深度重金属Cr含量均超过北京市土壤重金属Cr环境背景值（29.8 mg·kg-1），其中在 25～50 cm 土层中红黄土性物质所形成的土壤重金属Cr含量最高，为40.975 mg·kg-1，是北京土壤重金属Cr环境背景值的1.375倍。钙质岩类风化物和红黄土性物质中，25～50 cm土层深度Cr的含量高于0～25 cm，而平谷和密云地区的母质分别存在大量钙质岩类风化物和红黄土性物质，因而该区域部分样点25～50 cm中Cr含量高于0～25 cm，这与图2通过EBK插值得到不同土壤深度重金属Cr的空间分布结果一致。多重比较分析表明，不同母质类型对各区县不同土层中土壤重金属Cr含量无显著影响，单因素方差分析显示不同土壤深度组间无显著性差异。

3 结论

（1）从研究区整体看，不同土壤深度重金属Cr含量存在差异，海淀和丰台Cr平均含量最高。除昌平和平谷外，其余各区县Cr污染指数均已超过警戒线，应予以重视。空间分布上土壤Cr以城区周围为中心向外逐渐降低，不同土壤深度Cr分布格局基本一致，垂直方向上Cr含量存在显著相关性，其分布不同于以往随深度增加而降低的规律，表现出人为和自然地质的双重叠加影响。从Cr的污染等级分布看，处于轻度污染等级的区域占研究区总体的80.35%（0～25 cm深度）及84.84%（25～50 cm深度）。整个研究区土壤Cr存在富集现象，有一定污染风险。密云与平谷交界处土壤环境质量对北京市生活用水的水质影响较大，因此Cr含量过高是否会对水质产生影响则值得关注。

如果能够直观地读出细胞的膜电位，那么科学家们就可以对抑制神经元活化的电信号进行研究，而不是仅仅将研究局限于激活信号上。法国马赛港地中海神经生物学研究所的神经生物学家罗莎·柯萨特（Rosa Cossart）指出：由于钙离子成像等方法无法记录抑制性信号，我们并不清楚这类型号是如何影响脑活动的。

（2）针对不同功能区，土壤Cr含量在不同土壤深度中排序均为都市生活区＞农业保障区＞水源保护区，不同功能区对不同土壤深度Cr含量存在显著影响（P＜0.05）。不同种植模式下，不同土壤深度Cr含量无显著性差异，由于其含量均已超过北京市土壤重金属Cr环境背景值，而农药、肥料的投入增加是导致研究区Cr累积的重要原因，因此今后化肥农药施用及污水灌溉等方面应严格要求。

（3）本文主要针对土壤Cr的空间分布规律以及Cr污染评价方面进行研究，从定性角度分析不同功能区、不同种植模式及不同母质类型对不同土壤深度重金属Cr的影响，揭示了土壤Cr的分布特征，但并没有针对平原区农田土壤特点定量化地深入研究土壤Cr的来源，在后续研究中将继续探讨Cr分布特点并进行有关定量的析源分析，以期为农田土壤重金属污染防治提供科学依据。

参考文献：

[1]李 冰,王昌全,谭 婷,等.成都平原土壤重金属区域分布特征及其污染评价[J].核农学报,2009,23（2）：308-315.LI Bing,WANG Chang-quan,TAN Ting,et al.Regional distribution and pollution evaluation of heavy metal pollution in topsoils of the Chengdu Plain[J].Journal of Nuclear Agricultural Sciences,2009,23（2）：308-315.（in Chinese）

[2]陈红亮,谭 红,谢 锋,等.遵义东南部地区农业土壤重金属分布特征及风险评价[J].核农学报,2008,22（1）：105-110.CHEN Hong-liang,TAN Hong,XIE Feng,et al.Distribution and risk evaluation of heavy metals in agricultural soils of southeastern Zunyi City,China[J].Journal of Nuclear Agricultural Sciences,2008,22（1）：105-110.（in Chinese）

[3]Macfarlane G R,Koller C E,Blomberg S P.Accumulation and partitioning of heavy metals in mangroves：A synthesis of field-based studies[J].Chemosphere,2007,69（9）：1454-1464.

[4]田 媛,郭希娟,刘效兰.北京市不同功能区土壤重金属污染探究[J].环境科学与技术,2010,33（12F）：83-86.TIAN Yuan,GUO Xi-juan,LIU Xiao-lan.Exploration on the pollution of topsoil heavy metal in the different areas of Beijing[J].Environmental Science&Technology,2010,33（12F）：83-86.（in Chinese）

[5]江 宏,马友华,尹国庆,等.安徽省某县农田土壤重金属污染及潜在生态风险评价[J].农业资源与环境学报,2017,34（4）：352-359.JIANG Hong,MA You-hua,YIN Guo-qing,et al.Evaluation on heavy metal pollution and potential ecological risk of farmland soils in a county of Anhui Province,China[J].Journal of Agricultural Resources and Environment,2017,34（4）：352-359.（in Chinese）

[6]Lin Y P.Multivariate geostatistical methods to identify and map spatial variations of soil heavy metals[J].Environmental Geology,2002,42（1）：1-10.

[7]翟 航.长春市土壤重金属分布规律及土壤环境质量评价研究[D].长春：吉林大学,2007.ZHAI Hang.Research on distribution characteristic of the soil heavy metals and assessment for the soil environmental quality in Changchun[D].Changchun：Jilin University,2007.（in Chinese）

[8]Li Z,Ma Z,van der Kuijp T J,et al.A review of soil heavy metal pollution from mines in China：Pollution and health risk assessment[J].Science of the Total Environment,2014,468/469：843-853.

[9]陈同斌,雷 梅,杨 军,等.关于重金属污染土壤风险控制区划的研究与建议[J].中国科学院院刊,2014,29（3）：321-326.CHEN Tong-bin,LEI Mei,YANG Jun,et al.Discussion on zoning of soil environmental risk control and remediation contaminated by heavy metals on regional scale[J].Bulletin of Chinese Acadamy of Sciences,2014,29（3）：321-326.（in Chinese）

[10]Mapanda F,Mangwayana E N.The effects of long-term irrigation using water on heavy metal contents of soils under vegetables[J].Agriculture Ecosystems&Environment,2005,107（2）：151-165.

[11]崔德杰,张玉龙.土壤重金属污染现状与修复技术研究进展[J].土壤通报,2004,35（3）：366-370.CUI De-jie,ZHANG Yu-long.Current situation of soil contamination by heavy metals and research advances on the remediation techniques[J].Chinese Joumal of Soil Science,2004,35（3）：366-370.（in Chinese）

[12]吴 昊.重金属铬污染土壤治理方法研究进展[J].农业科技与装备,2015（8）：8-10.WU Hao.Research progress of soil controlling methods for heavy metal chromium pollution[J].Agricultural Science＆ Technology and Equipment,2015（8）：8-10.（in Chinese）

[13]He Z L,Yang X E,Stoffella P J.Trace elements in agroecosystems and impacts on the environment[J].Journal of Trace Elements in Medicine&Biology,2005,19（2/3）：125-140.

[14]Kumpiene J,Lagerkvist A,Maurice C.Stabilization of As,Cr,Cu,Pb and Zn in soil using amendments：A review[J].Waste Management,2008,28（1）：215-225.

[15]Jiao W,Chen W,Chang A C,et al.Environmental risks of trace elements associated with long-term phosphate fertilizers applications：A review[J].Environmental Pollution,2012,168（1）：44-53.

[16]Wang C,Yang Z,Zhong C,et al.Temporal-spatial variation and source apportionmentofsoilheavymetalsintherepresentativeriver-alluviation depositional system[J].Environmental Pollution,2016,216：18-26.

[17]张小敏,张秀英,钟太洋,等.中国农田土壤重金属富集状况及其空间分布研究[J].环境科学,2014,35（2）：692-703.ZHANG Xiao-min,ZHANG Xiu-ying,ZHONG Tai-yang,et al.Spatial distribution and accumulation of heavy metal in arable land soil of China[J].Environmental Science,2014,35（2）：692-703.（in Chinese）

[18]王彬武,李 红,蒋红群,等.北京市耕地土壤重金属时空变化特征初步研究[J].农业环境科学学报,2014,33（7）：1335-1344.WANG Bin-wu,LI Hong,JIANG Hong-qun,et al.Spatio-temporal variation of soil heavy metals in agricultural land in Beijing,China[J].Journal of Agro-Environment Science,2014,33（7）：1335-1344.（in Chinese）

[19]韩 平,王纪华,冯晓元,等.北京顺义区土壤重金属污染生态风险评估研究[J].农业环境科学学报,2015,34（1）：103-109.HAN Ping,WANG Ji-hua,FENG Xiao-yuan,et al.Ecological risk assessment of heavy metals in soils in Shunyi,Beijing[J].Journal of Agro-Environment Science,2015,34（1）：103-109.（in Chinese）

[20]索琳娜,刘宝存,赵同科,等.北京市菜地土壤重金属现状分析与评价[J].农业工程学报,2016,32（9）：179-186.SUO Lin-na,LIU Bao-cun,ZHAO Tong-ke,et al.Evaluation and analysis of heavy metals in vegetable field of Beijing[J].Transactions of the Chinese Society of Agricultural Engineering,2016,32（9）：179-186.（in Chinese）

[21]王学军,席 爽.北京东郊污灌土壤重金属含量的克立格插值及重金属污染评价[J].中国环境科学,1997,17（3）：225-228.WANG Xue-jun,XI Shuang.Kriging analysis and heavy metal pollution assessment for soil from eastern suburb of Beijing City[J].China Environmental Science,1997,17（3）：225-228.（in Chinese）

[22]赵 倩,马 琳,刘翼飞,等.北京东南郊典型地层重金属分布特征与潜在生态风险[J].环境科学,2016,37（5）：1931-1937.ZHAO Qian,MA Lin,LIU Yi-fei,et al.Distribution characteristics and potential ecological hazards assessment of soil heavy metals in typical soil profiles in southeast suburb of Beijing[J].Environmental Science,2016,37（5）：1931-1937.（in Chinese）

[23]李良忠,杨 彦,蔡慧敏,等.太湖流域某农业活动区农田土壤重金属污染的风险评价[J].中国环境科学,2013,33（s1）：60-65.LI Liang-zhong,YANG Yan,CAI Hui-min,et al.Health risk assessment of heavy metal pollution in farmland soils in Taihu Lake Basin agricultural area[J].China Environmental Science,2013,33（s1）：60-65.（in Chinese）

[24]刘琼峰,李明德,段建南,等.农田土壤铅、镉含量影响因素地理加权回归模型分析[J].农业工程学报,2013,29（3）：225-234.LIU Qiong-feng,LI Ming-de,DUAN Jian-nan,et al.Analysis on influence factors of soil Pb and Cd in agricultural soil of Changsha suburb based on geographically weighted regression model[J].Transactions of the Chinese Society of Agricultural Engineering,2013,29（3）：225-234.（in Chinese）

[25]Yang P G,Mao R Z,Shao H B,et al.The spatial variability of heavy metal distribution in the suburban farmland of Taihang Piedmont Plain,China.[J].Comptes Rendus Biologies,2009,332（6）：558-566.

[26]刘霈珈,吴克宁,罗 明,等.农用地土壤重金属超标评价与安全利用分区[J].农业工程学报,2016,32（23）：254-262.LIU Pei-jia,WU Ke-ning,LUO Ming,et al.Evaluation of agricultural land soil heavy metal elements exceed standards and safe utilization zones[J].Transactions of the Chinese Society of Agricultural Engineering,2016,32（23）：254-262.（in Chinese）

[27]陈同斌,郑袁明,陈 煌,等.北京市土壤重金属含量背景值的系统研究[J].环境科学,2004,25（1）：117-122.CHEN Tong-bin,ZHENG Yuan-ming,CHEN Huang,et al.Background concentrations of soil heavy metals in Beijing[J].Environmental Science,2004,25（1）：117-122.（in Chinese）

[28]郑袁明,陈 煌,陈同斌,等.北京市土壤中Cr,Ni含量的空间结构与分布特征[J].第四纪研究,2003,23（4）：436-445.ZHENG Yuan-ming,CHEN Huang,CHEN Tong-bin,et al.Spatial distribution patterns of Cr and Ni soils of Beijing[J].Quaternary Sciences,2003,23（4）：436-445.（in Chinese）

[29]郑袁明,陈同斌,陈 煌,等.北京市近郊区土壤镍的空间结构及分布特征[J].地理学报,2003,58（3）：470-476.ZHENG Yuan-ming,CHEN Tong-bin,CHEN Huang,et al.The spatial structure and distribution of Ni contents in soils of suburbs of Beijing[J].Acta Geographica Sinica,2003,58（3）：470-476.（in Chinese）

[30]钟 巧,王勇辉,焦 黎.夏尔希里地区土壤重金属含量特征及空间变异分析[J].水土保持研究,2016,23（3）：360-365.ZHONG Qiao,WANG Yong-hui,JIAO Li,et al.Characteristics and spatial variability of the soil heavy metals contents in Xiaerxili area[J].Research of Soil and Water Conservation,2016,23（3）：360-365.（in Chinese）

[31]Pierce F J,Dowdy R H,Grigal D F.Concentrations of six trace metals in some major minnesota soil series[J].Journal of Environmental Quality,1982,11（3）：416-422.

[32]赵 伟.重庆耕地土壤重金属Hg、Cr含量的空间变异性研究[J].河南农业科学,2009,38（11）：72-76.ZHAO Wei.Study of spatial variability of soil heavy metals Hg,Cr content in Chongqing[J].Journal of Henan Agricultural Sciences,2009,38（11）：72-76.（in Chinese）

[33]汤 民,张进忠,张 丹,等.果园土壤重金属污染调查与评价——以重庆市金果园为例[J].中国农学通报,2011,27（14）：244-249.TANG Min,ZHANG Jin-zhong,ZHANG Dan,et al.Pollution investigation and assessment of heavy metals in orchard soil：A case study in golden orchard of Chongqing[J].Chinese Agricultural Science Bulletin,2011,27（14）：244-249.（in Chinese）

[34]杨世琦,刘国强,张爱平,等.典型区域果园表层土壤5种重金属累积特征[J].生态学报,2010,30（22）：6201-6207.YANG Shi-qi,LIU Guo-qiang,ZHANG Ai-ping,et al.Analysis on accumulation characters of 5 kinds of heavy metals in orchard topsoil in typical production zone in China[J].Acta Ecologica Sinica,2010,30（22）：6201-6207.（in Chinese）

[35]邵学新,黄 标,赵永存,等.长江三角洲典型地区土壤中重金属的污染评价[J].环境化学,2008,27（2）：218-221.SHAO Xue-xin,HUANG Biao,ZHAO Yong-cun,et al.Pollution assessment of soil heavy metals in a representative area of the Yangtze River Delta region[J].Environmental Chemisiry,2008,27（2）：218-221.（in Chinese）