3.2地统计条件模拟的土壤Cu和Pb空间分布
仅通过1次条件模拟实现的结果在空间上存在较大不确定性，因此，利用多次模拟实现的期望值表示重金属的空间分布。基于地统计条件模拟和普通克里格插值的土壤Cu、Pb空间分布见图1。由图1可见，地统计条件模拟的Cu、Pb空间分布与普通克里格插值结果相似，普通克里格插值模拟结果相对平滑，地统计条件模拟因加入了随机因素，在空间上存在波动。从地统计条件模拟的方差分布(见图1(e)、图1(f))可知，在样点较稀疏和研究区边缘区域方差较大。总体来看，Cu含量在研究区域西北部相对较高，东南部较低。Pb在西部和中部区域有两个相对较高的区域，在东南部较低。
3.3土壤中Cu、Pb的污染概率
不同模拟结果显示了土壤重金属空间分布的随机性，单次模拟结果的污染区在空间上分布差异较大，因此仅基于一次模拟结果评估土壤重金属污染会存在较大的偏差。地统计条件模拟的优点是通过多次(≥1000次)的重复模拟，再现土壤重金属空间分布的波动性，基于多次模拟结果，评估土壤重金属空间分布的不确定性。进行土壤重金属污染评价时，需要采用污染概率的方法，即基于多次模拟获取的单点含量分布特征，估算超过污染阈值的概率P。
式中：n(x’)为模拟值大于污染阈值的个数，个;N为模拟次数，次。
基于地统计条件模拟和指示克里格插值的土壤Cu、Pb污染区域划分结果见图2。地统计条件模拟结果显示，研究区西部Cu污染概率较高，研究区西部和中部Pb污染概率较高。Cu、Pb的污染概率总体相对较低，这可能与研究区的污染区域较小，污染样点较少有关。地统计条件模拟结果与指示克里格插值的结果相似，但指示克里格插值的结果高估了污染样点周边的清洁样点。
3.4污染概率阈值的确定
