土壤调查加密采样布点优化方法

1.本发明涉及场地土壤环境调查和场地土壤生态环境损害调查技术领域，尤其涉及一种土 壤调查加密采样布点优化方法。

背景技术：

2.土壤污染已经形成点、面并发，并逐步向立体与综合演化的时空格局，呈现出流域性和 区域化发展的态势，以及与食物链复合交叉的特征。

3.土壤污染不仅造成了土壤资源的严重浪费，妨碍城市化建设和社会经济发展，还严重威 胁当地居民和周边人群的身体健康，为当地生态环境建设与恢复造成了严重的经济损失。因 此，污染场地治理修复刻不容缓。而污染场地环境调查工作是污染场地治理修复和生态环境 损害司法鉴定等工作的第一步，也是关键一步，是了解场地土壤污染现状的重要手段，调查 结果的准确程度直接关系到后续场地管理的经济、社会和技术效益，直接影响环生态境损害 司法鉴定结果和赔偿金额。制定科学高效的土壤调查采样布点方案不仅会节省大量的人力物 力，还会大大提高土壤调查结果的准确性和科学性。而现有的国家技术导则中，缺乏科学精 准的土壤污染调查采样点布设方法和采样点布设数量的确定方法，造成实际采样工作中无法 保证采集样品的代表性，经常需要通过密集的采样密度，来保证场地污染情况的真实性，从 而保证土壤污染调查的质量。由于不同的土壤调查采样方案会得到不同的土壤采样调查结果， 相对于实验室分析误差，不科学的土壤采样方案带来的误差影响更为显著。如今土壤污染调 查采样工作中的主要问题在于，野外采样工序复杂且经常受环境及人为因素的限制；采样点 位的布设密度缺乏统一标准，密集的采样密度虽然可以详细的反映场地污染情况，但需要耗 费较大的人力物力，也会参杂很多不具备代表性的样点，增加后续实验室化验的工作量，稀 疏的采样密度尽管可以节省费用，但容易丢失重要的场地污染信息。因此，提出科学高效合 理的土壤污染调查采样布点方法在土壤环境质量管理的技术和标准层面均具有重要意义。科 学高效的土壤污染调查采样布点方法应满足采样费用少、代表性强、采样精度高的要求，将 采样布点调查成本与调查不确定性导致的经济损失降到最低，在所需精度范围内用尽可能少 的采样点快速准确掌握场地土壤污染信息。

技术实现要素：

4.针对上述现有技术的不足，本发明提出一种土壤调查加密采样布点优化方法。

5.一种污染场地土壤调查加密采样布点优化方法，包括以下步骤：

6.步骤1：收集疑似污染场地或待查场地信息及历史资料，包括地理位置、地形地貌、土 壤类型、水文地质资料、生物资源分布、污染物资料、用地状况、生产工艺、污染物排放情 况、产品、原辅材料及中间产物；布设场地土壤初步调查采样点，采集场地土壤初步调查采 样点样品，检测样品中污染物量，得到样品的检测数据；根据样品的检测数据，初步判断场 地土壤中潜在污染物、污染范围及污染程度；

7.步骤2：对初步调查样品检测数据进行预处理：采用统计学方法spss分析初步调查采样 点样品检测数据分布情况，检验检测数据分布是否符合正态分布，若不符合正态分布需对检 测数据进行正态转换；若检测数据满足正态分布则执行步骤3，若不满足对检测数据进行正 态转换后执行步骤3；

8.步骤3：根据场地环境调查和生态环境损害调查确定相应的污染物污染阈值；

9.步骤3.1：场地环境调查：采用当地背景值、国家相关标准或土壤污染风险筛选值作为污 染物污染阈值；

10.步骤3.2：生态环境损害调查：首先根据初步调查采样点样品检测数据，利用原场点位法 计算土壤环境基线值，将计算出的土壤基线值确定为污染阈值；如无法采用原场点位法计算 土壤环境基线值，分析对照点值、当地背景值、国家相关标准和污染历史，确定污染阈值；

11.其中，利用原场点位法计算土壤环境基线值具体步骤为：

12.步骤3.2.1：收集场地相关资料，如包括地理位置、地形地貌、土壤类型、水文地质资料、 生物资源分布、污染物资料、用地状况、生产工艺、污染物排放情况、产品、原辅材料及中 间产物，全面分析资料后，进行布点采样；

13.步骤3.2.1.1：当场地条件为污染分布或污染范围明确时，利用收集场地的相关资料，确 定场地基线计算区域，即未污染或微污染区域；在该区域进行布点采样，对采样点样品进行 检测；

14.其中，布点采样为正常采样或低密度采样；

15.步骤3.2.1.2：当场地污染分布不明确或者污染范围未知时，收集污染场地相关资料，初 步确定场地基线计算区域，即未污染或微污染区域；在该区域进行布点采样，获取区域检测 数据；接着计算区域内污染物整体水平，确保计算区域污染物整体水平处于整个场地的最低 水平；

16.其中，布点采样为高密度采样；

17.步骤3.2.2：对区域样品检测数据进行分析，确定基线值，具体步骤如下；

18.步骤3.2.2.1：如果未检出物质，那么以检出限的1/2作为基线值；

19.步骤3.2.2.2：如果检出物质，使用累积频率分布法对到区域样品检测数据进行分析，得 到到检测数据分布图，分布图存在拐点的物质，采用低拐点作为基线值；分布图类似于直线 的物质，采用平均值作为基线值；分布图既不存在低拐点，也不呈现直线分布的物质，进一 步利用正态分布检验法确定基线；

20.其中，检测数据分布图呈正态分布时，采用平均值作为基线值；样点数据呈非正态分布 时，采用中位值作为基线值；

21.步骤4：采用指示克里格插值法预测场地土壤污染物含量超过污染阈值的概率，绘制污 染概率预测图，并选择污染不确定区域为重点加密布点区域；

22.将符合正态分布的初步调查采样点位数据添加到arcgis软件的地统计模块中，通过指示 克里格插值法预测土壤污染物含量超过污染阈值的概率，污染概率取值范围为0～1；概率值 越高则代表污染超标的可能性越大，反之则为清洁土壤；概率值介于高值和低值之间即为污 染不确定区域，并选择污染不确定区域为重点加密布点区域；

23.步骤5：利用差变系数法对重点加密布点区域的最优加密样本量进行计算，即计算

采样点样品，检测样品中污染物量，得到样品的检测数据；根据样品的检测数据，初步判断 场地土壤中潜在污染物、污染范围及污染程度；

47.本化工厂占地面积约96万m2，主要生产铬盐、农药、化学试剂等产品，化工厂所在地 地层结构主要为杂填土(主要由黏性土、碎石、混粒砂等组成)和粉质黏土。厂区内西北角 长期堆存铬渣，对周边环境及土壤存在一定的铬污染。已知场地内可能存在的污染物为六价 铬、铅、砷、汞、氯仿、苯、1,1,2-三氯乙烷、氯乙烯、α-六六六、β-六六六和六氯苯。

48.结合场地实际情况，运用系统网格法在研究区叠加200m

×

200m的网格布点，共选取初步 调查点位15个，采样深度为50cm，每个样点均采用五点法采样。初步调查采样样品经检验 发现六价铬全部超标，因此本实施例以六价铬为特征污染物进行研究。

49.步骤2：对初步调查样品检测数据进行预处理：采用统计学方法spss分析初步调查采样 点样品检测数据分布情况，检验检测数据分布是否符合正态分布，若不符合正态分布需对检 测数据进行正态转换；若检测数据满足正态分布则执行步骤3，若不满足对检测数据进行正 态转换后执行步骤3；图2为本发明数据预处理流程示意图。

50.本实施例中，由于重金属污染在土壤中的扩散和迁移具不确定性和不可逆性等特点，为 了降低重金属污染对土壤污染调查结果的影响，需要在加密采样布点优化方法研究前进行数 据预处理。图3为未经过数据预处理与数据预处理后的六价铬污染空间分布预测误差对比图。 由图3可知，数据预处理可以降低重金属污染分布预测的误差，极大地提高采样布点结果的 准确性，是土壤污染调查工作中必不可少的一项。借助spss 13.0软件对初步采样点位样品 检测数据进行统计分析和正态检验，得出六价铬含量数据的最大值、最小值、均值、标准差、 变异系数及数据分布形态等分析指标，若初步采样点位数据呈非正态分布，需将数据进行正 态变换使数据的分布相对连续且接近正态分布，才能进行下一步的污染概率预测。

51.步骤3：根据场地环境调查和生态环境损害调查确定相应的污染物污染阈值；

52.步骤3.1：场地环境调查：采用当地背景值、国家相关标准或土壤污染风险筛选值作为污 染物污染阈值；

53.步骤3.2：生态环境损害调查：首先根据初步调查采样点样品检测数据，利用原场点位法 计算土壤环境基线值，将计算出的土壤基线值确定为污染阈值；如无法采用原场点位法计算 土壤环境基线值，分析对照点值、当地背景值、国家相关标准和污染历史，确定污染阈值；

54.其中，利用原场点位法计算土壤环境基线值具体步骤为：

55.步骤3.2.1：收集场地相关资料，如包括地理位置、地形地貌、土壤类型、水文地质资料、 生物资源分布、污染物资料、用地状况、生产工艺、污染物排放情况、产品、原辅材料及中 间产物，全面分析资料后，进行布点采样；

56.步骤3.2.1.1：当场地条件为污染分布或污染范围明确时，利用收集场地的相关资料，确 定场地基线计算区域，即未污染或微污染区域；在该区域进行布点采样，对采样点样品进行 检测；

57.其中，布点采样为正常采样或低密度采样；

58.步骤3.2.1.2：当场地污染分布不明确或者污染范围未知时，收集污染场地相关资料，初 步确定场地基线计算区域，即未污染或微污染区域；在该区域进行布点采样，获取区

域检测 数据；接着计算区域内污染物整体水平，确保计算区域污染物整体水平处于整个场地的最低 水平；

59.其中，布点采样为高密度采样；

60.步骤3.2.2：对区域样品检测数据进行分析，确定基线值，具体步骤如下；

61.步骤3.2.2.1：如果未检出物质，那么以检出限的1/2作为基线值；

62.步骤3.2.2.2：如果检出物质，使用累积频率分布法对到区域样品检测数据进行分析，得 到到检测数据分布图，分布图存在拐点的物质，采用低拐点作为基线值；分布图类似于直线 的物质，采用平均值作为基线值；分布图既不存在低拐点，也不呈现直线分布的物质，进一 步利用正态分布检验法确定基线；

63.其中，检测数据分布图呈正态分布时，采用平均值作为基线值；样点数据呈非正态分布 时，采用中位值作为基线值；

64.本实施例利用原场点位法计算土壤环境基线值具体步骤如表1所示；

65.表1原场点位法基线的确定；

[0066][0067]

不同场地条件的基线确定步骤不同，主要分为场地功能分区明显，污染分布或污染范围 明确和污染源位置明确，污染分布不明确或者污染范围未知两种情况，第二种情况相对于第 一种情况，主要难点是基线计算区域的确定和计算区域合理性核实。数据分析方法优先采用 累积频率分布法，能较好反映数据的数据变化情况，其次采用正态分布检验法，方法简便， 但准确性较低。

[0068]

由于六价铬所有区域都存在超标情况，虽然总铬不是所有点位都超标，但数据大

部分远 超背景值。因此重点对as和pb计算基线，尽量选择深层土壤的数据。采用表1基线确定步 骤进行基线确定，图4为10m处数据的累积频率分布图。as元素累积频率分布曲线存在低拐 点，低拐点值为6.8mg/kg，小于低拐点的平均值为5.93mg/kg。pb元素累积频率分布近似直 线，取平均值13.38mg/kg作为基线。总铬和六价铬采用表1计算出基线，结果如表2所示。

[0069]

表2基线计算结果

[0070][0071]

步骤4：采用指示克里格插值法预测场地土壤污染物含量超过污染阈值的概率，绘制污 染概率预测图，并选择污染不确定区域为重点加密布点区域；

[0072]

将符合正态分布的初步调查采样点位数据添加到arcgis软件的地统计模块中，通过指示 克里格插值法预测土壤污染物含量超过污染阈值的概率，污染概率取值范围为0～1；概率值 越高则代表污染超标的可能性越大，反之则为清洁土壤；概率值介于高值和低值之间即为污 染不确定区域，并选择污染不确定区域为重点加密布点区域；由于采样调查结果的不确定性 主要出现在污染不确定区域即污染物含量过渡区域，所以本发明中提取污染不确定区域作为 重点加密布点区域。

[0073]

本实施例绘制污染概率预测图并确定重点加密布点区域；

[0074]

以初查点位数据为基础，土壤六价铬背景值8.4mg/kg作为污染阈值，假设初查点位数据 超过污染阈值则认定为受到污染，初查点位数据小于污染阈值则认定为未受到污染。选用指 示克里格插值法绘制污染概率预测图，模拟预测研究区内土壤中六价铬污染超过污染阈值的 概率，为了详细划分污染区域、未污染区域与污染不确定区域，将研究区污染概率分为五个 区段，分别为(0，0.3)、(0.3，0.5)、(0.5，0.75)、(0.75，0.9)和(0.9，1)。图5 为污染概率预测图。由图5可知，(0.9，1)区间的概率值较大，污染超标的可能性较大， 假设为污染土壤区域，(0，0.3)区间的概率值较小，污染超标的可能性较小，假设为清洁 土壤区域，(0.3，0.9)区间为污染过渡区域，是污染预测误差较大区域，因此提取污染概率 区间(0.3，0.9)所在区域为重点加密布点研究区域。

[0075]

步骤5：利用差变系数法对重点加密布点区域的最优加密样本量进行计算，即计算场地 地块污染物平均浓度与场地土壤污染阈值的差异与“估计标准差”的比值；根据这一比值， 确定各污染概率值区间合理的加密采样布点数量；

[0076]

比值用τ表示，计算公式如下：

[0077][0078]

式中，μ1表示地块污染物平均浓度；s0为步骤3所述污染阈值，σ为地块污染物估计 标准差；在τ值较大的区域分配较少的加密点位，在τ值较小的区域分配较多的加密点位。

[0079]

本发明可通过计算各污染概率区间差变系数的大小，给出合理分配各污染概率区

间加密 布点数量的建议。在τ值较大的区域分配较少的加密点位，在τ值较小的区域分配较多的加 密点位。

[0080]

本实施例比较各污染概率区间布点效率；

[0081]

在五个污染概率区间(0，0.3)、(0.3，0.5)、(0.5，0.75)、(0.75，0.9)和(0.9， 1)所在区域内分别选取10个点位作为基础布设点位，并各加密布设10个点位，比较分析各 污染概率区间加密布点效率。结果表明：各污染概率区间布点效率大不相同，清洁区域和污 染区域加密布点效率很低，污染概率在(0.5，0.9)区间时，布点效率较高，说明在污染过渡 区域内污染概率大于0.5的区间内加密布点，得到的结果有较高准确性。图6为各污染概率 区间内平均布设10个点的点位图，图7为各污染概率区间内平均布设20个点的点位图，表 3为各污染概率区间布点效率。

[0082]

其中，布点精度计算公式如下：

[0083][0084]

布点效率计算公式如下：

[0085][0086]

a0—场地真实污染面积；a1—初查点位预测出的污染面积；ai—加密点位预测出的污染 面积。p1—初查布点的布点精度；pi—加密布点的布点精度；ei—加密布点的平均效率；n— 加密点位数量。表3为各污染概率区间布点效率。

[0087]

表3各污染概率区间布点效率；

[0088][0089][0090]

本实施例利用差变系数法给出各污染概率区间合理的加密采样布点数量建议；

[0091]

差变系数法的计算公式如下：

[0092][0093]

其中，μ1为地块污染物平均浓度；s0为土壤污染物基线值或背景值，σ为估计标准差。 差变系数越大，所需样本量越小；差变系数越小，所需样本量越大。

[0094]

利用差变系数法对各污染概率区间合理加密样本量大小估计计算，以上述各污染概率区 间分别布设20个点位为样本数据，带入差变系数公式，给出各污染概率区间合理的加密采样 布点数量建议。结果表明：为了得到准确可靠的污染范围预测结果，各污染概率区间最佳布 点数量从大到小排序为：(0.5，0.75)，(0.75，0.9)，(0.3，0.5)，(0，0.3)，(0.9， 1)。因此在使用加密采样布点优化方法时，可在(0.5，0.75)区间内分配较多的加密点位， (0.75，0.9)和(0.3，0.5)区间内分配较少的加密点位，(0，0.3)和(0.9，1)区间内可 在明显

清洁区域或污染源所在区域选择性加密布点。表4为各污染概率区间最优加密样本量 排序。

[0095]

表4各污染概率区间最优加密样本量排序

[0096][0097]

本方法布点优化方法与其他方法对比分析：

[0098]

为验证加密采样布点优化方法的高效快捷，以初查15个点位为基础，选取除初查点位外 的50个点位制定对照布点方案进行对比分析：

[0099]

方案一：普通均匀加密布点法，在整个研究区域内均匀布设采样点位，如图6和7所示。

[0100]

方案二：污染概率加密布点法，在(0.9，1)和(0，0.3)概率区间所在区域内分别布设 10个点，在污染不确定区域(0.3，0.9)内均匀布设30个点位。点位布设图如图8所示。

[0101]

方案三：加密采样布点优化方法，根据差变系数大小分配点位，即在差变系数小的区域 多分配点位。根据表2计算结果，点位分配方式为(0.5，0.75)内布设20个点位，(0.75， 0.9)内布设15个点位，(0.3，0.5)、(0.9，1)和(0，0.3)内分别布设5个点位，点位 布设图如图9所示。方案三为本发明方法。

[0102]

以布点精度和布点效率作为采样精度评价指标，分析各布点方案布点结果，表5为各布 点方案结果统计表。结果表明：本发明所述的一种污染场地土壤调查加密采样布点优化方法 在污染范围预测、布点精度和布点效率上都要好于普通均匀加密布点法和污染概率加密布点 法。

[0103]

表5各加密布点方案结果统计表

[0104][0105]

案例2：以某制药厂为例，对本发明实施例进行详细说明；

[0106]

步骤1：根据污染场地信息及历史资料，设置污染场地土壤初步调查采样布点，采集污 染场地土壤初步调查采样点位数据，初步判断污染场地土壤中潜在污染物、污染范围及污染 程度；

[0107]

本制药厂占地面积约26万m2，所在地土壤类型包括杂填土、粉土、粉质黏土、粉砂、 细砂、中砂、砾砂和圆砾等。

[0108]

由于地块内原土被大量挖掘出场，部分区域的挖掘深度超过15m，使该地块的土壤受到 了很大的扰动。各个地块的表层均有大量外来土壤，这些外来土壤中也可能含有各类污染物， 且本制药厂生产历史较长，在生产和污染排放过程中污染物容易在土壤中累积，可能存在的 污染物为高环多环芳烃(苯并(a)蒽、苯并(b)荧蒽、苯并(k)荧蒽、苯并(a)芘、茚并(1,2,3-cd) 芘、二苯并(a,h)蒽)、氯仿和铜、铅、镉、砷、汞和镍等。

[0109]

结合场地实际情况，在研究区上叠加150m

×

150m的网格布点，共选取初查点位15个点， 采样深度为50cm，初步调查采样点位分布图见如图10。以15个初步采样点位数据为基础， 针对主要污染物镉进行研究。

[0110]

步骤2：对初步调查采样点位数据进行预处理：采用统计学方法利用spss软件分析初步 调查采样点位数据分布情况，检验数据是否符合正态分布，若不符合正态分布需对数据进行 正态转换以降低采样布点结果误差；表6为土壤重金属初始含量数据和转换后数据的偏度和 峰度

[0111]

表6土壤重金属镉初始含量数据和转换后数据的偏度和峰度

[0112]

[0113]

在spss 13.0软件中对比分析了各种转换方法对六种重金属元素的适用性，发现对数转 换(log10)方法较适合本研究区土壤重金属含量数据的正态转换，且转换后数据的偏度和峰 度值均接近正态分布数据的偏度和峰度值。

[0114]

以初查点位数据为基础，选用指示克里格法预测研究区镉含量超过污染阈值的概率，以 该地区镉背景值0.16mg/kg为污染阈值，假设采样点位中镉含量超过0.16mg/kg则认定为受到 污染，镉含量小于0.16mg/kg则认定为未受到污染，以此为划分污染区域与未污染区域的研 究依据。以初步采样点位数据为基础的污染概率预测图如图11所示。

[0115]

将图11中的污染概率详细均分为10个区间，分别为：(0，0.1)、(0.1，0.2)、(0.2， 0.3)、(0.3，0.4)、(0.4，0.5)、(0.5，0.6)、(0.6，0.7)、(0.7，0.8)、(0.8，0.9) 和(0.9，1)，由于(0，0.2)区间的概率值较小，此区域假设为清洁土壤区域，(0.9，1) 区间概率值大，此区域假设为污染土壤区域。(0.2，0.9)区间为污染过渡区域即污染不确定 区域，由于此区域为污染预测误差较大区域，因此提取污染概率区间(0.2，0.9)所在区域为 重点加密布点研究区域。

[0116]

将布点优化方法与其他方法对比分析：

[0117]

为验证加密采样布点优化方法的高效快捷，以初查15个点位为基础，选取除初查点位外 的50个点位制定对照布点方案进行对比分析：

[0118]

方案一：系统网格法布点方案，在整个研究区域内均匀布设60个采样点位。点位布设图 如图12所示。

[0119]

方案二：污染概率加密布点法，在(0.9，1)和(0，0.2)概率区间所在区域内分别布设 10个点，在污染不确定区域(0.2，0.9)内均匀布设40个点位。点位布设图如图13所示。

[0120]

方案三：加密采样布点优化方法，在(0，0.1)、(0.1，0.2)、(0.4，0.5)和(0.5， 0.6)内各布设3个点位，(0.2，0.3)和(0.3，0.4)内各布设2个点位，在(0.6，0.7)、 (0.7，0.8)、(0.8，0.9)和(0.9，1)内分别布设8个点位，点位布设图如图14所示。方 案三为发明方法。

[0121]

以布点精度和布点效率作为采样精度评价指标，分析各布点方案布点结果，表7为各布 点方案结果统计表。结果表明：本发明所述的一种污染场地土壤调查加密采样布点优化方法 在污染范围预测、布点精度和布点效率上都要好于普通均匀加密布点法和污染概率加密布点 法。

[0122]

表7各加密采样布点方案结果统计

[0123][0124]

本发明在污染概率加密布点法的基础上加入了差变系数法，对土壤污染调查加密布点方 法进行了优化，以初步调查点位数据为基础，通过指示克里格插值法模拟预测污染物含量超 过污染阈值的概率，绘制污染概率预测图并确定加密布点研究区域，在污染概率加密布点法 的基础上引入差变系数法，得出各污染概率区间合理加密布点数量大小排序，对各污染概率 区间的点位分配方式给出指导性建议。本发明在本实施例中能有针对性的以较少采样点位快 速准确掌握场地污染情况，大幅度提高了污染场地土壤调查的布点精度和布点效率。技术特征：

1.一种土壤调查加密采样布点优化方法，其特征在于，包括以下步骤：步骤1：收集疑似污染场地或待查场地信息及历史资料，包括地理位置、地形地貌、土壤类型、水文地质资料、生物资源分布、污染物资料、用地状况、生产工艺、污染物排放情况、产品、原辅材料及中间产物；布设场地土壤初步调查采样点，采集场地土壤初步调查采样点样品，检测样品中污染物量，得到样品的检测数据；根据样品的检测数据，初步判断场地土壤中潜在污染物、污染范围及污染程度；步骤2：对初步调查样品检测数据进行预处理：采用统计学方法spss分析初步调查采样点样品检测数据分布情况，检验检测数据分布是否符合正态分布，若不符合正态分布需对检测数据进行正态转换；若检测数据满足正态分布则执行步骤3，若不满足对检测数据进行正态转换后执行步骤3；步骤3：根据场地环境调查和生态环境损害调查确定相应的污染物污染阈值；步骤4：采用指示克里格插值法预测场地土壤污染物含量超过污染阈值的概率，绘制污染概率预测图，并选择污染不确定区域为重点加密布点区域；将符合正态分布的初步调查采样点位数据添加到arcgis软件的地统计模块中，通过指示克里格插值法预测土壤污染物含量超过污染阈值的概率，污染概率取值范围为0～1；概率值越高则代表污染超标的可能性越大，反之则为清洁土壤；概率值介于高值和低值之间即为污染不确定区域，并选择污染不确定区域为重点加密布点区域；步骤5：利用差变系数法对重点加密布点区域的最优加密样本量进行计算，即计算场地地块污染物平均浓度与场地土壤污染阈值的差异与“估计标准差”的比值；根据这一比值，确定各污染概率值区间合理的加密采样布点数量。2.根据权利要求1所述的一种土壤调查加密采样布点优化方法，其特征在于，所述步骤3具体为：步骤3.1：场地环境调查：采用当地背景值、国家相关标准或土壤污染风险筛选值作为污染物污染阈值；步骤3.2：生态环境损害调查：首先根据初步调查采样点样品检测数据，利用原场点位法计算土壤环境基线值，将计算出的土壤基线值确定为污染阈值；如无法采用原场点位法计算土壤环境基线值，分析对照点值、当地背景值、国家相关标准和污染历史，确定污染阈值。3.根据权利要求1所述的一种土壤调查加密采样布点优化方法，其特征在于，所述步骤3.2具体为：利用原场点位法计算土壤环境基线值具体步骤为：步骤3.2.1：收集场地相关资料，如包括地理位置、地形地貌、土壤类型、水文地质资料、生物资源分布、污染物资料、用地状况、生产工艺、污染物排放情况、产品、原辅材料及中间产物，全面分析资料后，进行布点采样；步骤3.2.1.1：当场地条件为污染分布或污染范围明确时，利用收集场地的相关资料，确定场地基线计算区域，即未污染或微污染区域；在该区域进行布点采样，对采样点样品进行检测；其中，布点采样为正常采样或低密度采样；步骤3.2.1.2：当场地污染分布不明确或者污染范围未知时，收集污染场地相关资料，

初步确定场地基线计算区域，即未污染或微污染区域；在该区域进行布点采样，获取区域检测数据；接着计算区域内污染物整体水平，确保计算区域污染物整体水平处于整个场地的最低水平；其中，布点采样为高密度采样；步骤3.2.2：对区域样品检测数据进行分析，确定基线值，具体步骤如下；步骤3.2.2.1：如果未检出物质，那么以检出限的1/2作为基线值；步骤3.2.2.2：如果检出物质，使用累积频率分布法对到区域样品检测数据进行分析，得到到检测数据分布图，分布图存在拐点的物质，采用低拐点作为基线值；分布图类似于直线的物质，采用平均值作为基线值；分布图既不存在低拐点，也不呈现直线分布的物质，进一步利用正态分布检验法确定基线；其中，检测数据分布图呈正态分布时，采用平均值作为基线值；样点数据呈非正态分布时，采用中位值作为基线值。4.根据权利要求1所述的一种土壤调查加密采样布点优化方法，其特征在于，所述步骤5比值用τ表示，计算公式如下：式中，μ1表示地块污染物平均浓度；s0为步骤3所述污染阈值，σ为地块污染物估计标准差；在τ值较大的区域分配较少的加密点位，在τ值较小的区域分配较多的加密点位。

技术总结

本发明提供一种土壤调查加密采样布点优化方法，涉及场地土壤环境调查和场地土壤生态环境损害调查技术领域；收集场地信息及历史资料；布设初步调查采样点，采样检测；检测数据预处理；根据初步采样数据，利用原场点位法计算土壤环境基线值；分析污染物土壤基线值或当地背景值、国家相关标准、污染历史，结合调查目的，确定污染阈值；采用指示克里格插值法预测土壤污染概率，绘制污染概率预测图，确定重点加密布点区域；利用差变系数法对加密布点区域的最优加密样本量进行估计计算，确定场地土壤调查加密采样布点优化方案。本发明引入差变系数法和污染概率加密布点法，对场地土壤调查采样布点方法优化，实现了高效、准确的少采样目标。标。

技术研发人员：唐玉兰 刘强 郭祁萱 李婷 傅金祥 刘钰 王洪达 李汶龙 李继伟

受保护的技术使用者：沈阳建筑大学

技术研发日：2022.02.23

技术公布日：2022/5/30