摘要：鉴于传统调查布点方法对土壤污染物平均含量的估计精度较高，但对污染区范围的估计精度不能满足修复治理需求的问题，同时，为了准确估计土壤污染区面积及其空间位置，提出了土壤污染详查样点优化方法.在初步调查的基础上，首先利用地统计条件模拟方法预测土壤污染概率，基于污染概率和土壤污染物含量局部空间变异确定加密布点的优先区域，并根据污染物含量的空间变化趋势布设样点，然后根据优化后的土壤污染调查布点方案估计污染区面积和空间位置.最后，以某Cd污染场地为例，验证了布点优化方法的效果.结果表明，土壤污染调查加密布点方法显著提高了污染区面积及其空间位置的估计精度.案例场地土壤Cd污染区面积的预测误差为4.10%，污染区空间位置的精度为86.10%.土壤污染调查布点优化方法在保证土壤污染调查精度的同时，相比于传统调查方法可显著降低土壤污染调查的样本量.
关键词：土壤污染调查；地统计条件模拟；污染概率；局部空间变异；污染区范围；布点优化；
作者简介：谢云峰(1981—)，男，副研究员(博士)；E-mail：xieyf@craes.org.cn；
1引言(Introduction)
土壤采样调查是获取土壤污染物空间分布信息最重要的手段，采样调查结果的精度直接影响污染风险评价结果的准确性和风险管理决策的合理性.土壤污染调查包括土壤样点布设、样品采集、污染物含量分析等环节.实际工作中，通常认为污染物分析方法的准确性是影响污染物调查准确性的最主要因素(Crumblingetal.，2001)，而忽略了土壤采样布点方案的重要性.大量研究表明，污染物在土壤中的空间分布表现出明显的空间变异性，人类活动影响越大的区域，局部变异程度越大(Thompson，1996；丛鑫等，2009；杜平等，2006；张娟等，2014；郑一等，2003).针对空间变异性较大的环境要素，样点布设方案是影响调查结果准确性最主要的因素之一.Jenkins等(1997)对土壤中三硝基甲苯污染的调查结果表明，至少95%的变异度(统计方差)是由采样位置导致，而含量分析(室内分析和现场分析)手段对变异度的贡献不超过5%.其他类似研究也表明，土壤采样导致的不确定对污染物含量测定不确定性的贡献超过50%(Argyrakietal.，1997；Theocharopoulosetal.，2001；Jenkinsetal.，1999).因此，科学合理的土壤采样布点方案对保障污染调查结果的精度非常重要.现有的土壤污染调查布点方法主要包括判断性采样和非判断性采样(姜成晟等，2009)，其中，判断性采样主要根据已有先验知识设计采样布点方案，并在潜在的高污染风险区域加大采样密度(UKEnvironmentAgency，2000)；当缺乏场地污染物分布的背景信息时，就只能采取非判断性采样方法，如随机采样、均匀网格布点采样等(Thompsonetal.，1995；USEPA，1989).传统的土壤污染调查布点方法主要用于对污染物总体(平均含量)的最佳估计(Brusetal.，1999)，样本量主要取决于污染物含量的空间变异程度.土壤污染治理过程中，污染调查主要关注目标污染物的超标程度及污染区范围.因此，以总体估计为目标的传统土壤污染调查布点方法对土壤污染范围的估计精度通常不能满足修复决策的需求(刘庚等，2013；谢云峰等，2010).近年来，应用地统计学方法来提高土壤污染调查精度已成为研究热点之一(D'Or，2005；Demougeot-Renardetal.，2004；Juangetal.，2005；VanGroenigenetal.，1999；VanToorenetal.，1997)，该方法基于土壤污染物空间分布的自相关性，优化土壤调查布点空间布局，可提高土壤污染调查效率(Burgessetal.，1981；Demougeot-Renardetal.，2004；Englundetal.，1993；阎波杰等，2008；赵倩倩等，2012).虽然基于地统计学和条件模拟方法的样点布设方法效率最高(Jonesetal.，2003)，但在土壤污染调查过程中却很少用于土壤污染调查布点优化(Verstraeteetal.，2008).
为了获得准确的土壤中污染物空间分布信息，土壤污染调查通常包括污染初步调查、污染详查等多个阶段.初步调查的主要目的是识别土壤主要污染物及潜在污染区域，通常样本量较少.污染详查是在初步调查基础上，在潜在的污染区域增加样点，确定污染区的范围及其污染程度.土壤污染调查方案的误差主要包括污染区被低估和清洁区被高估(Marchantetal.，2013；Ramseyetal.，2002)，其中，前者会导致污染区面临的污染风险不能得到有效控制，后者会导致不必要的修复投入.为了获取准确的污染区信息，通常需要增加样本量，但这会导致采样分析成本的增加.高效的采样方案是将采样调查成本与调查不确定性导致的经济损失的总成本降到最低(Ramseyetal.，2002).采样方案优化的目的就是要寻求降低污染修复不确定性的最佳样本量(Demougeot-Renardetal.，2004).土壤污染物的空间分布受污染来源、环境条件、污染物性质等因素的综合影响，其在空间上表现出不同程度的空间相关性和变异性，对土壤污染物空间变异性的描述准确与否是影响调查结果的关键.本研究结合土壤污染调查的特定需求，提出基于污染概率和污染物局部空间变异特征的土壤污染调查加密布点方法，以提高土壤污染调查方法对污染区范围和污染程度的估计精度，并为土壤污染调查提供方法学支持.
2土壤污染调查加密布点方法(Samplingdesignoptimizationprocedurefordetailedsoilpollutioninvestigation)
土壤污染调查结果的不确定性主要出现在污染物含量过渡区域(刘庚等，2013；谢云峰等，2010；Xieetal.，2011)，为此，该研究针对污染调查结果的不确定性，提出土壤污染调查加密布点的工作流程和方法(见图1).土壤污染调查加密布点的2个核心问题分别为确定需要加密布点区域和样点布设方法.
2.1加密布点区域的确定方法
由于土壤污染治理仅关注污染物含量超过相关环境标准或修复目标值的区域，因此，提高污染区范围的估计精度就显得尤为重要，加密布点法正是基于这一需求而提出.由于初步调查阶段已经获得了一定的污染物分布信息，所以在加密详查阶段只需要针对污染分布信息不确定性较大的区域进行补充调查即可，其中，不确定性区域是指污染物空间分布精度低于修复决策需求精度的区域.
为了定量评估土壤污染调查的不确定性，该研究引入土壤污染概率方法.基于初步调查数据，利用概率制图方法预测土壤污染物超过环境标准或修复目标值的概率，常用的概率制图方法有地统计条件模拟方法、指示克里格方法等.其中，地统计条件模拟方法包括多种模拟算法，如序贯高斯模拟、序贯指示模拟等.污染概率的取值范围为0~1，概率值越高，可优先判定为污染土壤；相反，污染概率值越低，可优先判定为清洁土壤.概率制图结果中，概率值介于高值和低值之间者即为不确定性区域，需要进一步补充调查确认.假定某污染土壤地统计条件模拟的污染概率阈值范围为0.1~0.8，设定污染概率阈值和清洁概率阈值分别为0.5和0.3，则污染概率为0.5~0.8者为污染区域，0.1~0.3者为清洁区域，0.3~0.5者即是需要加密调查的区域.
不确定性区域污染概率值较低的可能原因为：①区域内污染物含量较低；②区域属污染区域，并且样本量较少.为了进一步探究其具体原因，该研究引入局部变异特征方法.基于初步调查数据，分析土壤污染物含量的局部变异特征(包括变异系数、方差、自相关性等)，如果局部变异性较大，表明土壤中污染物含量空间分布差异较大；反之，则表明污染物含量空间分布差异较小.对于局部变异性较大者，通常是污染物含量高值区向低值区的过渡区域，也是调查结果不确定性较大的区域；对于变异性较小者，通常是高值集中或低值集中的区域，调查结果的可靠性较高.因此，根据土壤污染物的局部变异系数，将土壤污染调查结果划分为不确定性区域和确定性区域.假定某污染土壤局部变异系数为20%~200%，设定变异系数阈值为100%，则变异系数为100%~200%者为不确定性区域；低于100%者为确定性区域.
综合污染概率和局部空间变异系数确定的污染调查不确定性区域，即为污染调查加密布点的目标区域.
2.2不确定性区域样点布设方法
不确定性区域样点布设包括加密样点的数量和样点的空间位置.其中，加密样点数量主要与不确定性程度相关，不确定性较大的区域，加密布设的样本量也较大；样点的空间位置主要与污染物含量空间变化趋势相关，主要利用趋势分析方法分析土壤污染物空间变化规律，沿着土壤污染物含量变化的方向布设加密样点.
本研究提出的污染调查加密布点方法的主要目的是为提高污染区范围的估计精度.在初步调查结果的基础上，结合污染概率和局部变异系数方法确定加密布点的目标区域，再根据土壤污染物含量分布的空间变异性及其变化趋势，确定加密样点的布设方案.该方法可优化加密布点的位置，降低加密布点的数量，提高加密布点的效果，从而在保证调查精度的前提下，降低调查成本.
3加密布点方法案例验证(Validationofthesamplingdesignoptimizationprocedurefordetailedsoilpollutioninvestigation)
3.1案例区概况
案例数据来源于某重金属污染场地，场地面积约14.50km2.按照200m间隔进行均匀采样，在部分高污染区域适当增加样本量，共采集359个土壤样品.土壤污染调查结果表明，土壤重金属Cu、Pb、As、Cd等污染物都存在不同程度的污染.以该场地土壤Cd污染为例，开展土壤污染调查详查加密布点优化方法研究.
3.2样点加密布点方案
案例验证研究过程中并不实际开展土壤污染初步调查布点取样，以及初步调查结果分析和详查加密布点工作.而是利用案例场地已有的359个调查数据，采用空间抽样的思路，模拟开展土壤污染初步调查和加密详查布点过程.具体操作步骤为：首先基于案例数据的359个样点数据，进行模拟的土壤污染初步调查.根据图1的工作流程可知，土壤污染物空间变异特征研究和土壤污染不确定性区域确定是土壤详查加密布点的2个最重要的环节.地统计学的半方差分析方法是最常用的空间变异特征研究手段之一，为了获取比较准确的土壤污染物的空间分布规律，需要有足够的样本量.因此，在初步调查阶段，将研究区域划分为10×10的网格，落在网格内的土壤样点作为初步调查样点，当网格内有多个土壤样点时，随机选取其中一个，由此共获得土壤初步调查样点97个，样点间平均距离约为386m.在初步调查的97个样点数据的基础上，利用本研究提出的加密布点方法进行加密布点.具体步骤为：基于初步调查数据，利用地统计学方法分析场地土壤Cd含量(w(Cd))的空间分布规律.利用条件模拟方法预测该场地土壤Cd污染概率(图2a).基于污染概率预测结果，设定污染概率阈值(Pt)和清洁概率阈值(Ct)，污染区域确定方法如式(1)所示.土壤Cd污染概率阈值和清洁概率阈值分别设定为0.8和0.2，基于污染概率划定的不确定性区域见图2b；在此基础上，结合土壤污染局部变异特征(图2c)，将局部变异性大于变异系数阈值(CVt)的区域划定为不确定性区域(图2d)，变异系数阈值设定为局部变异系数最大值的75%(式(2)).综合污染概率和局部变异系数的结果，即为土壤污染详查布点的优先区域，根据土壤污染物空间结构分析结果，沿着污染物含量变化的方向确定加密样点的位置(图3a).由于该研究是模拟研究，如果在最佳的采样位置没有样点数据，就选择邻近样点作补充，土壤详查加密样点为57个，布点方案见图3b.将加密布点后的污染调查结果与案例场地359个数据获得的结果进行对比，评价加密布点的效果.
式中，Rp为污染概率分区，Z(x)为条件模拟预测的土壤污染物含量，Zc为土壤污染评价标准，Pt为污染概率阈值，Ct为清洁概率阈值，Rcv为污染变异系数分区，CVx为局部变异系数，CVt为变异系数阈值.
3.3数据处理方法
利用GS+7.0软件进行土壤污染物含量的空间结构特征分析.样点污染物含量局部变异特征是在样点VORONOI图的基础上，借助ArcGIS10.1的GeostatisticalAnalyst工具，计算每个样点及其邻近样点的变异系数.采样网格、初步调查样点设计及所有空间制图均在ArcGIS10.1软件中实现.土壤污染物含量条件模拟及污染概率计算在GSLIB(GeostatisticalSoftwareLibrary)中实现(Journeletal.，1998).地统计学条件模拟方法较多，该研究采用最常用的算法之一序贯高斯模拟方法(SequentialGaussianSimulation，SGS)(谢云峰等，2015)，该方法算法简单、灵活、计算方便，其基本思路为：根据现有样点数据计算待模拟点污染物浓度的条件概率分布，从该分布中随机取值作为模拟实现；将得到的每一个模拟值，连原始样点数据一起作为条件数据，进入下一个点的模拟.
3.4结果与讨论
3.4.1土壤Cd统计特征的估计精度
由表1可见，土壤Cd污染初步调查样点(97个)与污染详查样点加密后(154个)的统计特征很相似，平均值差异仅为0.01mg·kg-1.加密详查后样本的变异系数降低.与总体样本相比，初步调查和加密详查这2个阶段采样的Cd平均值都偏高，误差为5.40%.变异系数较总体分别降低2.79%和6.71%.初步调查平均值的估计精度较高，而加密详查并没有进一步提高平均值的估计精度.在污染详查阶段，由于在土壤污染空间变异较大的区域增加了样点，因此，其变异系数降低.
3.4.2土壤污染区面积的估计精度
土壤污染调查重点关注的是污染信息的识别精度.初步调查和加密详查阶段，根据样点w(Cd)超标率(表2)估算的污染区面积所占比例分别为68.04%和70.13%，比所有样本的估算结果分别高3.14%和5.23%.污染概率预测结果表明，当污染概率阈值为0.8时，污染概率预测的污染区面积所占比例在53.58%~57.84%之间，比样点超标率估计结果低7.06%~16.39%.基于超标率估算污染区面积，意味着当某个采样网格内的土壤样点污染物含量超标时，则判定该网格超标.样点加密详查后，增加的样点都位于污染概率较高的区域，因此，总体样点中污染区域样点的比例增加，导致污染面积估计结果增大.
初步调查和加密详查这2个阶段估计的污染区面积非常接近，样点加密后污染概率预测的污染区面积仅增加0.16%，初步调查与加密详查估算的面积均小于总体样本的估计结果，污染面积低估4.10%.为了评价污染区范围空间位置的预测精度，将不同采样阶段预测的污染区范围与总体样本预测的结果进行空间差值运算，并根据差值结果将污染区空间位置预测精度分为相同、低估和高估3种情况.相同表示污染程度预测结果一致，低估表示污染区被预测为清洁区，高估表示清洁区被预测为污染区(图4).从污染区的空间位置精度来看，初步调查污染区面积预测的准确度为79.35%，分别有12.45%的区域污染程度被低估，8.20%的区域污染程度被高估.加密详查后，污染区面积预测的准确度提高到86.10%，污染程度被低估和被高估的面积分别降至9.00%和4.90%.
土壤Cd平均值估计结果表明，在初步调查阶段，其估计精度就已达到94.00%以上，而污染区的估计精度仅为79.35%.表明在土壤污染调查过程中，平均值或土壤污染统计特征的估计精度，并不能反映污染区范围的估计精度.土壤污染治理过程中，污染区空间分布信息比平均值更重要，直接影响到修复成本的估计.本研究提出的土壤污染详查加密布点方法，在保证土壤污染总体平均含量估计精度的前提下，显著提高了污染区面积的估计精度；加密详查后，污染区面积的估计误差为4.10%，空间位置精度为86.10%，比初步调查精度提高了6.75%；土壤污染调查的样本量显著降低，初步调查和加密详查的样本量仅为总体的42.90%.
本研究的样点优化思路是在不确定性较大的区域内增加样点，不确定性区域的界定标准为条件模拟的污染概率和局部变异系数.从图2可知，不确定性区域主要分布在污染区边缘，在这些区域增加样点密度，能显著提高污染区空间位置精度.初步调查过程中，污染程度被低估时，污染区域被误判为清洁区域(见图4左下角和左上角的绿色区域)；样点优化过程中，清洁区域不会补充调查样点，因此，优化后的结果仍然是被低估.污染程度被低估与初步调查布点、污染概率阈值选择有关.由于没有污染物分布相关的背景信息，网格随机采样布点法对总体平均含量和变异程度的预测精度较高，对局部污染信息的预测精度较低.在初步调查前，收集场地污染源排放、土地利用方式、土壤理化性质、水文地质条件等影响污染物空间分布的相关信息，辅助调查样点设计，可以提高对污染区识别的精度(Falketal.，2011).污染概率阈值选择对加密点的空间分布有较大影响，如果选择的污染概率阈值过低，就会导致被高估的区域不能被识别；概率阈值过高，则会导致不确定性区域增大，需加密的样本过多，从而降低加密效率.本研究为了获取较大的不确定性区域，选择了较高的污染概率阈值和较低的清洁概率阈值，用于检验样点优化方案的效率.在具体应用中，应结合研究区的特点和调查目标，选择适宜的污染概率阈值，进一步提高样点优化方案的效率.加密详查样点优化过程中，基于污染概率和局部变异系数筛选出不确定性较大的区域，该研究并没有在这些区域增加样点，而是根据已有的样点数据，基于距离邻近原则，用邻近样点替代最佳位置的样点.增加的样点在空间位置上并不是最优化的，这可能会降低样点优化的效率.实际应用中在最佳的空间位置补充样点，应该会取得更好的调查效果.
本研究提出的加密布点方法的核心是在污染预测结果不确定性的区域，根据污染物空间分布规律补充调查样点.如图1所示，在污染物空间分布、污染概率预测、预测结果不确定性评价等阶段都应用了地统计学方法.根据地统计学方法的基本假设，应用该方法时要求污染物空间分布具有显著的空间自相关性.大量的研究结果表明，重金属、多环芳烃等污染物在土壤中的空间分布都表现出明显的空间相关性(胡克林，2004；郑一等，2003).因此，地统计学方法是适用的.对某些污染物，如化工场地的氯代烃污染等，这类污染物主要是通过泄漏释放到土壤中，然后通过土壤孔隙进一步向下迁移.在水平空间上，存在泄漏的区域就会检出污染物，没有泄漏的区域就不存在污染(韩春梅等，2009)，因此，这类污染物在空间上自相关性较差，本研究提出的加密布点方法就不适用.土壤中污染物空间分布受污染源分布及释放特征、区域环境条件、污染物性质及环境行为特点等多种因素的综合影响，在不同尺度上会表现出不同的空间分布规律.针对具体区域开展污染调查时，需综合考虑污染物空间分布的影响因素，同时可借鉴前期研究和其它类似研究的成果，初步分析土壤中污染物的空间分布特征，在此基础上，进行初步调查布点.基于初步调查结果，应用地统计学方法研究污染物空间分布规律，如果污染物具有较好的空间自相关性，就可以采用本研究的方法进行加密布点优化，否则，本研究的方法就不适用.加密布点是在初步调查结果的基础上，通过辨识污染物的空间分布规律，结合污染调查的要求，开展详查布点优化.因此，初步调查的可靠性会直接影响加密布点的效果.地统计学应用半方差分析研究污染物的空间自相关性.相关研究表明，样点数量和空间分布会直接影响半方差分析结果的准确性(Goovaerts，1999).从样点数量来看，由于污染物类型、研究区域条件的差异，不同研究的结论不太一致，通常认为样点数小于60时，难以获得较准确的半方差(秦耀东，1998).在具体应用时，可根据半方差函数的拟合效果，评估样点数是否足够.从样点空间分布来看，为评估污染物在不同距离和不同方向上的空间分异规律，初步调查样点应尽可能在研究区域内均匀分布，在不同距离和方向上都有足够的样点数用于分析污染物的空间分布规律，可帮助提高加密布点优化的效率.
4结论(Conclusions)
1)土壤污染调查布点方法对土壤污染物含量的估计精度较高，案例场地土壤中Cd平均值的预测误差为5.40%，变异系数的预测误差为6.71%.
2)土壤污染调查布点方法显著提高了污染区面积和污染区空间位置的估计精度，案例场地土壤Cd污染区面积的预测误差为4.10%，污染区空间位置的精度为86.10%.
3)在初步调查的基础上，基于污染概率和局部变异系数优化土壤污染调查布点，在保证土壤污染调查精度的同时可显著降低土壤污染调查的样本量.
