空间插值方法在土壤污染物研究方面的应用.pdf

油气田环境保护2023年8 月ENVIRONMENTALPROTECTION OFOIL&GASFIELDSVol.33No.425空间插值方法在土壤污染物研究方面的应用张慧斌12吴倩1杨磊”李薇刘晓丹王占生（1.中国石油集团安全环保技术研究院有限公司；2.华北电力大学环境科学与工程学院；3.中国石油天然气股份有限公司辽河石化分公司）摘要空间插值方法是土壤污染物信息分析的常用方法，文章分析了4种常用空间插值方法及其特点，梳理了空间插值方法在土壤污染物空间分布、采样方案、制图、源解析和风险评价方面的应用，为土壤污染物研究提供思路，最后总结了空间插值方法在土壤污染物研究中应用的优势，提出了空间插值方法与信息技术结合、加强空间插值方法之间组合应用等未来发展方向。关键词空间插值方法；土壤污染物；信息技术DOl:10.3969/j.issn.1005-3158.2023.04.006Application of Spatial Interpolation Method in the Investigation on Soil PollutantsZhang Huibinl:22.North China Electric Pozver University;3.PetroChina Liaohe Petrochemical Company)ABSTRACT Spatial interpolation methods were commonly used to analyze soil pollution information.This paper analyzed four common spatial interpolation methods and their characteristics;summarized theapplication of spatial interpolation methods in soil pollutant spatial distribution,sampling schemes,mapping,source analysis,and risk assessment;provided new ideas for the research on soil pollution.Meanwhile,the advantages of spatial interpolation methods in soil pollution research were summarized,andfuture development directions such as combining spatial interpolation methods with information technologyand strengthening the combined application of different spatial interpolation methods were proposed.KEYWORDS spatial interpolation method;soil pollutants;information technology0引言近年来，我国土壤污染日益严峻，工业蓬勃发展产生的废物、废水，化肥、杀虫剂的大量使用等是造成土壤污染的主要因素。对土壤污染状况进行调查研究，掌握土壤污染物超标情况，了解土壤污染物的空间分布特征、污染来源等信息，对土壤污染的防治工作有重要意义1-2 。空间插值方法作为空间分析技术重要的分支，可以有效解决土壤污染物采样数据不完整且不能连续表达的问题，依据已知的数据，通过某种数学关系推出未知数据3 ，对土壤文章编号：10 0 5-3 15 8（2 0 2 3)0 4-0 0 2 5-0 5Wu QianYang Lei?Li Wei?Liu XiaodanWang Zhansheng(1.CNPC Research Institute of Sa fety&Environment Technology;污染物进行空间插值研究，能够实现研究区域的空间连续预测，可以降低成本，有助于掌握土壤污染物信息。空间插值方法基于“地理学第一定律的基本原理”，通过已知要素的空间信息及关联特征，求未知要素的特征-5 。主要应用在土壤污染物空间分布、土壤污染物采样方案、土壤污染物制图、土壤污染物源解析和土壤污染物风险评价5 个方面。本文基于对空间插值方法及其特点的总结，论述其在土壤污染信息分析中的作用，以期为土壤污染的基金项目：国家重点研发计划，污染场地挥发类有机污染物传输机制与扩散通量项目“复杂场地挥发类有机污染物多相传输扩散的环境风险评估及量化表征（2 0 2 0 YFC1807103）。第一作者：张慧斌，华北电力大学2 0 2 0 级在读硕士，研究方向：地统计学及空间分析。通信地址：北京市昌平区北农路2 号华北电力大学，10 2 2 0 6。E-。通讯作者：吴倩，2 0 0 4年毕业于北京大学环境科学专业，高级工程师，现在中国石油集团安全环保技术研究院有限公司从事土壤与固废处理处置研究工作。通信地址：北京市昌平区黄河北街1号院1号楼，10 2 2 0 6。E-mail：q i a n w u c n p c.c o m.c n。26防治工作提供指导。1常用空间插值方法及其特点目前常用的空间插值方法有反距离加权法、径向基函数法、普通克里格法、局部多项式法。反距离加权法是一种加权平均内插法，插值精度依赖于反距离的幂值；径向基函数法是一类多变量插值函数的插值方法，可以反映研究目标的整体和局部变化趋势；普插值方法优点距离越近影响程度越大，根反距离加权法据点的远近程度赋予相关方法简单，考虑了空间自相关性权重进行加权平均径向基函数所生成的曲面不径向基函数法仅能够反映整体变化趋势而且可以反映局部变化趋势对空间分布的数据求线性普通克里格法最优、无偏内插估计的一种插值方法使用既定邻域内的点对指局部多项式法定阶数的多项式进行拟合2空间插值方法在土壤污染物研究方面的应用2.1在土壤污染物空间分布研究中的应用土壤污染物数据具有空间维度特性。传统的统计学在研究土壤污染物变化趋势时，由于忽视了土壤污染物的空间属性，从而很难探究其空间变化规律，而空间插值方法充分考虑了变量的空间属性，在研究土壤污染物空间分布特征方面有很好的效果6 1。空间插值方法在小尺度、中尺度以及大尺度区域内土壤污染物空间分布特征研究中都有应用。陈晓晨等7 运用ArGIS软件地统计分析模块中普通克里格插值对小尺度下首钢厂区附近土壤表层重金属进行了插值分析；郑一等8 以天津地区表层土壤中PAHs为例,对其进行中尺度普通克里格插值研究，发现研究区存在一定程度的PAHs污染；浦江等9 对杞麓湖径流区7 种重金属进行大尺度反距离加权插值分析，结果显示，杞麓湖径流区Cd和Hg污染比较严重。近年来，随着研究区内监测数据的累积，许多学者开始研究土壤污染物时空分布问题。土壤污染物时空分析的基础是获得多时期的土壤样品数据，而数油气田环境保护技术研究与应用通克里格法利用原始数据和半方差函数的结构性，对未采样点的区域化变量进行无偏最优的一种插值方法；局部多项式法是一种局部加权最小二乘拟合方法，使用既定邻域内的点对指定阶数的多项式进行拟合。表1分析了这4种常用空间插值方法的特点、优点、缺点和适用性，学者可以根据研究方向需要，选择不同插值方法进行研究。表1空间插值主要方法及其特点特点计算格式简单、计算速度快、逼近效果理想，对数据误差不敏感，有良好的平滑性既考虑了插值点和采样点之间位置关系，也考虑了变量的空间相关性用趋势面考虑全部数据点反应趋势性变化；用距离法反映局部特征，可以准确地体现平滑表面据的获取需要大量的人力和物力，空间插值方法在不增加测量成本的前提下，可以利用已知数据很好地了解土壤污染物时空分布特征。夏敏等10 1运用GIS数据统计分析方法和地统计学空间插值方法相结合对河南省封丘县土壤重金属污染地块时空分布研究，明确了土壤重金属历史变化趋势及空间分布特征；马妍等11运用普通克里格插值法对2 0 0 0 一2 0 2 0 年2 0 年间我国表层土壤PAHs污染数据进行插值分析，很好地描述了表层土壤PAHs在时空上的变异特征，反映了土壤PAHs多个时期的分布变化规律。2.2 在土壤污染物采样方案中的应用选择合理的采样方案在土壤污染物研究中至关重要12-141，对研究区域内污染物的空间分布影响很大15 。合理的采样方案应该考虑采样数目、位置和采样单元等因素16-17 ，使采样方案具有代表性、合理性和低成本性。传统采样方案多数不考虑土壤空间变异属性，在采样过程中主观性较强，一般需要通过增加采样点数量来提高预测的精度18 ，采样效率低，难以满足实际工作需要。空间插值方法可以弥补传统采样存在的缺陷，它充分考虑了研究区域土壤的空Vol.33No.4缺点适用性主要依赖反距离的幂值，常规采样或不规幂值越高，插值越精确；易则间隔采样的稀受极值影响疏数据快速插值有高程数据和相易产生过于理想化的偏差关属性的快速插值保证了数据的局部最优，数据量大，适合不能保证数据的总体最用普通克里格插优；一些有意义的异常值值法，精度高可能被“光滑效应”光滑掉常规采样或不规数值很少过原监测点，易则间隔采样的稀出现条带现象疏数据快速插值2023年6 月间变异属性，能通过预先采集少数样点计算空间变异模型相关参数和插值预测误差，确定合理的采样单元和采样空间位置19 ，这一方法能够在满足一定精度的前提下，减少采样数量，节约成本。在国外，研究人员很早开始研究土壤污染采样方案问题。McBratney和Websrer181在19 8 3年运用克里格最小估计方差所得出的土壤采样效率比传统采样方法的采样效率高出39 倍；Trangmar201、Chang等2 1 发现在相同采样精度下，普通克里格法设计的采样方案所需的样品要比传统采样方法少。在国内，阎波杰等2 2 通过设置不同的采样密度，基于不同空间插值方法优化采样数目使插值精度最大化；胡清等2 3 用不同空间插值方法分析对比了快速检测点位与采样策略优化后实验室检测点位两种采样方案，发现两种采样方案的插值结果相似，都可以反映污染物分布趋势，但是第二种采样方案优化了采样点数量和位置，降低了采样成本；李俊颖2 4、李润林2 5、李楠等2 6 研究表明，空间插值方法与辅助变量相结合，在保证插值精度的前提下，能够有效优化采样数；王宇航等2 7 通过交叉验证法研究了采样点数目与反距离加权插值精度的关系，研究发现，随着采样点数目减少，插值误差逐渐增加，采样点减少到总采样点数目的7 0%左右，插值精度会发生突变。2.3在土壤污染物制图中的应用土壤污染物制图在环境领域应用时间较长，能够快速的、直观的表达污染物分布，通常用来描述土壤污染物的空间分布特征，便于分析污染物在土壤内的污染来源、迁移转化及变化趋势2 8 。空间插值方法是目前用于土壤污染物制图的一种重要方法，绘制空间分布图时能够保留土壤污染物空间分布信息，还可以减弱最大值和最小值对土壤污染物空间分布影响，准确描述土壤污染物空间变异信息，使制图结果更符合实际情况。在国外，Imperato等2 9 利用普通克里格法绘制了意大利那不勒斯市区表层土壤中重金属铜、锌、镉、铅等值线图，直观刻画了土壤污染物空间分布信息；White等30 1用地统计学分析美国周边土壤中Zn的采样数据,得到了合理的的指数半变异函数，在此基础上用普通克里格插值法绘制了土壤重金属Zn的空间分布图，减弱了最大值和最小值的影响；Alimh等311在GIS软件中利用反距离权重插值模块，绘制了研究区内农业地表土壤样品中重金属含量的单元素分布图，准确描述了土壤重金属空间变异信息。国内也有很多学者使用不同的空间插值方法开张慧斌等：空间插值方法在土壤污染物研究方面的应用27展土壤污染物制图的相关研究工作。王昌全等32 等通过构建遗传径向基函数神经网络绘制土壤有效锌空间插值图，发现该方法扩大了对土壤相关空间信息的搜索范围,在一定程度上避免了“平滑效应”；张路伟等33 以空间差异度为评价标准，分析比较了3种不同空间插值方法的土壤属性制图效果，发现普通克里格法精度最高，减弱了极值的影响。2.4 在土壤污染物源解析中的应用土壤污染物具有降解慢、多种形态、潜伏时间长、毒性强和隐藏性等特点34。污染物一旦进人土壤后会发生一系列物理、化学和生物过程，严重危害生物体及生态系统，很难进行控制和治理，因此，准确解析土壤污染物来源是控制污染最有效的措施之一。常用的污染物源解析法有主成分分析法、正定矩阵因子分解法和同位素比值法等，这些方法忽略了土壤污染物的空间属性和空间分布信息，无法满足实际工作的需要35 。源解析法常与空间插值法联合使用，可以获得不同污染源空间分布情况。土壤中污染物的自然来源以土壤母质为主，人为来源包括工业源、生活源和农业源等35 。在国外，Hoehun等36 以美国阿拉巴马州安尼斯顿地区的工业区和住宅区土壤中11种重金属为研究对象，采用主成分分析法和克里格插值法，实现了对土壤重金属的来源解析以及空间分布表征。在国内，陈丹青等37 利用普通克里格插值法描述了Cd和Hg来源贡献率在研究区域内的空间分布特征。2.5在土壤污染物风险评价中的应用土壤是一个时空连续的变异体，具有很强的空间变异特征，传统的土壤污染风险评价方法操作简单，但权重确定主观性强，只能给出定性分析结果，不能准确反映土壤污染风险评价的不确定性本质。空间插值方法可以对风险评价结果进行估值预测，准确找到存在风险的区域，避免出现决策失误的风险。王如松等38 提出了基于风险评价的土地治理措施，应用反距离加权对不同人体健康风险概率进行空间插值，发现污染物对儿童的健康风险影响较大；路一帆等39 以苏州某地区铅蓄电池厂遗留的污染地块作为研究对象，用普通克里格插值法对场地致癌风险和非致癌风险进行插值模拟，发现该区域非致癌风险基本忽略，但存在一定的致癌风险；毕如田等40 1以我国某焦化场地中特征性污染物苯并（a)芘为研究对象，参照北京市污染场地土壤筛选值标准，运用空间插值组合预测模型界定了土壤污染范围。283结论与展望1)空间插值方法可以将离散的、有限的采样数据插值成连续的空间信息，绘制出土壤污染物分布图，分析出土壤污染物的空间分布特征。2)应用空间插值方法能充分考虑到土壤污染物的空间相关性，提高采样方案的合理性和科学性，降低采样成本。3)从地理学理论的角度分析土壤污染物与环境因素之间的相关性，结合场地历史信息，揭示导致土壤污染的污染源，为土壤修复提供依据。空间插值方法在土壤污染物研究领域已经应用了5 0 多年，取得了飞速发展。未来在该领域的可能发展方向如下：1)随着空间插值方法与大数据技术、遥感技术、云服务技术和全球定位技术等结合，快速挖掘出土壤污染物的空间数据信息，开发预测精度更高的空间插值模型。2)空间插值方法组合使用不断增加。单一空间插值方法不能适应复杂数据和特定场景的缺点，根据数据的特征选择合适的空间插值方法进行组合能够有效提高插值精度。3)随着土壤采样技术和信息技术的发展，使用三维空间插值方法准确刻画土壤污染物空间迁移转化和时空变异特征。参考文献1许学慧，兰亭,刘海云，等.活化剂添加对重金属污染土壤性状的影响J.干旱区资源与环境，2 0 2 2，36（2）：105-109.2彭雨璇，谢恩泽，王美艳，等.场地土壤污染信息三维空间插值研究进展及展望J.土壤通报，2 0 2 1，5 2（5）：1244-1250.3李海涛，邵泽东.空间插值分析算法综述J.计算机系统应用,2 0 19,2 8(7):1-8.4崔清松.空间插值算法在地质建模中的应用D.成都：西南石油大学，2 0 10：5 1-5 2.5王士博，王勇.空间插值方法在恶性肿瘤研究中的应用J.中华肿瘤防治杂志，2 0 2 1,2 8（5）：40 1-40 6.6冯经昆.西南山区某垃圾焚烧厂周边土壤重金属空间分布及其污染评价D.桂林：广西师范大学，2 0 14：2 0-2 2.7陈晓晨，崔岩山.城市表层土壤中重金属的小尺度空间分布一以首钢厂区附近小区域为例门.中国科学院研究生院学报，2 0 10,2 7（2）：17 6-18 3.8郑一，王学军，刘瑞民，等.天津地区土壤多环芳烃的克里格插值与污染评价J.中国环境科学，2 0 0 3，2 3（2）：油气田环境保护技术研究与应用113-116.9浦江，张翠萍，刘淑娟，等.杞麓湖径流区不同湿地沉积物重金属污染特征及潜在生态风险评价.农业资源与环境学报，2 0 2 1,38（5）：7 5 5-7 6 3.10夏敏,赵炳梓,张佳宝,基于GIS的黄淮海平原典型潮土区土壤重金属积累研究.土壤学报，2 0 13,5 0（4)：6 8 4-6 9 2.11马妍，程芦，阮子渊，等.近2 0 年中国表层土壤中多环芳烃时空分布特征及源解析.环境科学，2 0 2 1，42（3）：1065-1072.12 ROCCO C,DURO I,ROSA DI,et al.Composite vs.dis-crete soil sampling in assessing soil pollution of agricul-tural sites affected by solid waste disposalJ.Journal ofgeochemical exploration,2016,170:30-38.13 ANDRONIKOV S V,DAVIDSON D A,SPIERS R B.Variability in contamination by heavy metals:samplingimplicationsJJ.Water air and soil pollution,2000,120(1-2):29-45.14顾振飞，曹见飞，王翠秀，等.考虑人为因素影响的重金属含量分析用土壤样品采集方法研究一一以龙口为例土壤通报,2 0 2 1,5 2（3：6 7 0-6 7 8.15杨琳,朱阿兴，秦承志，等.一种基于样点代表性等级的土壤采样设计方法.土壤学报，2 0 11，48（5）：9 38-9 46.16张宏凯，左锐，王金生，等.加油站泄漏污染物的迁移分布规律.中国环境科学，2 0 18，38（4）：15 32-15 39.17张雷，张帝，谯兴国，等.石油类污染物在季节性冻土中迁移转化规律研究J.环境保护科学，2 0 19，45（3）：106-109.18 MCBRATNEY AB,WEBSTER R.How many observa-tions are needed for regional estimation of soil propertiesJJ.Soil science,1983,135(3):177-183.19胡碧峰，王佳昱，傅婷婷，等.空间分析在土壤重金属污染研究中的应用J.土壤通报，2 0 17，48（4）：10 14-10 2 4.20VANGROENIGENJW,STEIN A,ZUURBIER R.Optimization of environmental sampling using interactiveGISJJ.Soil technology,1997,10(2):83-97.21J CHANG Y H,SCRIMSHAW M D,EMMERSON R HC,et al.Geostatistical analysis of sampling uncertaintyat the tollesbury managed retreat site in blackwater estu-ary,essex,UK:kriging and cokriging approach to mini-mise sampling densityJ.Science of the total environ-ment,1998,221(1):43-57.22阎波杰，潘瑜春，赵春江.区域土壤重金属空间变异及合理采样数确定J.农业工程学报，2 0 0 8，2 4（增刊）：260-264.23胡清，童立志，王宏，等.基于现场快速检测技术的土壤采样优化策略J.环境工程，2 0 2 0，38（12）：16 3-16 7.2 4李俊颖，吴克宁，宋文.农用地土壤重金属空间变异及合理采样数确定J.中国农业资源与区划，2 0 2 0,41（9）：64-71.Vol.33No.42023年6 月25李润林，姚艳敏，唐鹏钦，等.县域耕地土壤锌含量的协同克里格插值及采样数量优化.土壤通报，2 0 13，44(4):830-838.26李楠，徐东瑞，吴杨洁.土壤养分含量的协同克里格法插值研究.浙江农业学报，2 0 11,2 3（5）：10 0 1-10 0 6.27王宇航，缪亚敏，杨昕.采样点数目对反距离加权插值结果的敏感性分析.地理信息世界，2 0 12，10（4）：31-35.2 8 邓峰.KRIGING方法及其在环境科学中的应用J.重庆环境科学，19 9 0（3）：33-37.29 IMPERATO M,ADAMO P,NAIMO D,et al.Spatialdistribution of heavy metals in urban soils of Naples city(Italy）JJ.E n v ir o n m e n ta l p o llu tio n,2 0 0 3,12 4（2):247-256.30 WHITE J G,WELCH R M,NORVELL W A.Soil zincmap of the USA using geostatistics and geographic infor-mation systemsJJ:Soil science society of America jour-nal,1997,61(1):185-194.31 ALI M H,MUSTAFA A A,EL-SHEIKH A A.Geo-chemistry and spatial distribution of selected heavy met-als in surface soil of Sohag,Egypt:a multivariate statisti-cal and GIS approachJ.Environmental earth sciences,2016,75(18):247-256.32董敏，王昌全，李冰，等.基于GARBF神经网络的土壤有效锌空间插值方法研究J.土壤学报，2 0 10，47（1）：42-50.33张路伟，赵彦锋，程道全，等.基于空间差异度的不同插张慧斌等：空间插值方法在土壤污染物研究方面的应用2021,16(5):259-270.35陈雅丽，翁莉萍，马杰，等.近十年中国土壤重金属污染源解析研究进展J.农业环境科学学报，2 0 19,38（10）：2219-2238.36 HOEHUN H,OLSON J R,LING B,et al.Analysis ofheavy metal sources in soil using kriging interpolation onprincipal componentsJ.Environmental science&tech-nology,2014,48(9):4999-5007.37陈丹青，谢志宜，张雅静，等.基于PCA/APCS和地统计学的广州市土壤重金属来源解析.生态环境学报，2016,25(6):1014-1022.38阳文锐，王如松，黄锦楼，等.反距离加权插值法在污染场地评价中的应用.应用生态学报，2 0 0 7（9）：2 0 13-2 0 18.39路一帆，陆胤,蔡慧，等.铅蓄电池厂遗留场地重金属污染分析及健康风险评价.环境工程，2 0 2 2，40(1)：135-140.40刘庚,毕如田，张朝，等.某焦化场地苯并（a)芘污染空间分布范围预测的不确定性分析.环境科学学报，2013,33(2):587-593.（修回日期2 0 2 2-0 9-30）（录用日期月2 0 2 2-10-10)（编辑王蕊）29值方法土壤属性制图效果及影响因素分析一以新乡县为例刀.河南农业科学，2 0 13，42（11：7 6-8 0.34马喆，王美娥，霍彦慧，等.场地土壤重金属镉和铅复合污染毒性阈值的推导及其应用.生态毒理学报，（上接第2 4页）2)相比于单一微生物修复，硫酸钠/过氧化钙预氧化-微生物联合修复实验两者联合修复方法对土壤石油烃的降解效果更好，污染土壤中石油烃降解率可达8 5.13%，较直接进行生物降解提高了39.6 6%。3)在微生物修复过成中，第二周、第三周时间内土壤中石油烃降解速率最大，微生物数量最多,生物修复六周后，土壤中石油烃含量的变化不再明显。修复后土壤中残留石油烃仍以C26C4o 大分子石油烃为主，说明大分子石油烃生化效果较差。联合修复后土壤的pH值由9.31降低至7.2 1，土壤性质得到了改良。4)过硫酸钠/过氧化钙预氧化-微生物联合修复方式在石油烃污染土壤修复中具有良好的应用前景。参考文献1张金永，叶倩，王明新，等.机械化学法修复柴油污染土壤的效率、产物及影响.环境科学学报,2 0 2 1,41(2）：1-8.2吴昊，孙丽娜，王辉，等.活化过硫酸钠原位修复石油类污染土壤研究进展.环境化学，2 0 15，34（11）：2 0 8 5-2 0 9 5.3骆永明，滕应.中国土壤污染与修复科技研究进展和展望.土壤学报，2 0 2 0,5 7（5）：1137-1142.4KULIK N,GOI A,TRAPIDO M,et al.Degradation of pol-ycyclic aromatic hydrocarbons by combined chemical pre-oxi-dation and bioremediation in creosote contaminated soil J.Journal of environmental management,2006,78(4):382-391.5高春阳，韩占涛，陈昌照，等.过氧化氢、过碳酸钠活化过硫酸钠氧化一一微生物降解联用技术修复原油污染土壤.化工学报，2 0 19,39（1)：7 1-7 6.6中华人民共和国生态环境部，国家市场监督管理总局.土壤环境质量建设用地土壤污染风险管控标准（试行)：GB366002018S.北京：中国环境科学出版社，2 0 18.7吴昊，孙丽娜，王辉，等.CaO2/Fe2+活化过硫酸钠对石油类污染土壤的修复效果.环境化学，2 0 16,35（4)：6 2 3-6 2 8.（修回日期2 0 2 2-10-0 2）（录用日期月2 0 2 2-11-0 3)（编辑王蕊）