一种模拟有机污染物时空迁移归趋分布的方法技术

本发明专利技术实施例公开了一种模拟有机污染物时空迁移归趋分布的方法，具体涉及模拟有机污染物迁移分布领域，包括以下步骤：S1、建立数据库；S2、对研究区域进行网格划分；S3、建立四级逸度模型：整合逸度计算公式、Z值计算公式、D值计算公式、平流计算公式，将这些公式通过Matlab编程平台进行程序化；S4、采用四级逸度模型并配合平流规则进行四级多介质逸度模型模拟。本发明专利技术通过对研究区域进行网格化，解决了空间异质性；为每个网格配置一套独立的参数文件，且对应模拟的年份，更符合实际情况；加入平流规则，污染物在网格间的迁移通过平流规则实现，模拟结果更接近实际情况；相对于现有的方法，更接近实际情况。

A method to simulate the temporal and spatial migration and fate distribution of organic pollutants

【技术实现步骤摘要】
一种模拟有机污染物时空迁移归趋分布的方法
本专利技术实施例涉及模拟有机污染物迁移分布领域，具体涉及模拟有机污染物时空迁移归趋分布的方法。
技术介绍
自20世纪以来，全球进入了快速农业化，工业化和城市化，人类对人工化学品的需求与日俱增，这些化学品通过各种途径释放到环境(大气、水体、土壤、沉积物和植物)中，其中有毒有害的化学品造成环境污染，对人类身体健康产生严重影响。持久性有机污染物(POPs)是一类全球普遍存在，具有环境持久性、生物积累性、生物毒性以及可远距离传输等特性的化学物质。传统对有机污染物把控的方法是通过对目标研究区域的环境监测，环境检测及风险评估。传统方法的优势是可以立刻得到有机污染物在环境中的一个瞬时浓度，但是有机污染物在环境中是个非稳态过程，很多自然因素会对它们在环境中的含量产生影响，传统方法无法对有机污染物进行持续性监测以及未来的预测。国内外研究：Mackay等[1]人发现有机污染物在环境多介质中的环境行为不是无规则的，自然界的定律决定他们的分配和迁移转化速率，从而建立了逸度模型。逸度模型基于质量平衡定律应用有毒有机化学品在多介质生态环境中的定量表达式建立数值模型，通过输入和输出，模拟预测污染物在现实环境中的行为归趋。田慧等人[2]使用逸度模型以广州为研究区域，模拟多溴联苯醚(PBDEs)在环境中迁移归趋情况。曹红英等人[3]构建四级逸度模型以北京为研究区域，模拟六六六(HCHs)四种异构体在环境中迁移随时间变化的情况；敖江婷等[4]人建立了IV及多介质动态逸度模型，模拟了1952-2010年期间黄河流域下游HCHs的4中同分异构体在该区域的归趋情况，模拟结果很好的揭示了HCHs在环境中随着农药的施用和金庸所形成的长期变化趋势。现有技术存在以下不足：以往的单箱式多介质逸度模型仅适合于研究空间特性均匀的小尺度区域，但是在研究大尺度区域时，这种单箱式模型的准确性和可靠性就比较低。因为在真实环境中，空间特性是并不均匀的：(1)地理环境因素包括不同的土地利用面积等具有空间异质性；(2)气候因素包括风向、风速、降雨量和温度等具有空间异质性；(3)污染物排放强度具有空间异质性，例如广东省的工业园区主要分布在珠三角洲，废旧电子拆解区分布在贵屿和清远。如果将整个区域统一在一起模拟计算，会将上述的具有空间异质性的因素的特点掩盖掉，对污染物的归趋模拟计算影响较大。与本专利技术相关之参考资料及文献如下：[1]MackayD.Multimediaenvironmentalmodels:thefugacityapproach[M].CRCpress，2001.[2]田慧，郭强，毛潇萱，etal.广州地区典型多溴联苯醚迁移和归趋行为模拟[J].中国环境科学，2014，34(3):758-765.[3]曹红英，梁涛，陶澍.北京地区50多年来有机氯农药迁移与残留的动态模拟与预测[J].中国科学D辑，2005，35(10).[4]敖江婷.Ⅳ级逸度模型对典型有机污染物环境行为的动态模拟[D].大连理工大学，2008.。
技术实现思路
为此，本专利技术实施例提供一种模拟有机污染物时空迁移归趋分布的方法，通过对研究区域进行网格化，解决了空间异质性；为每个网格配置一套独立的参数文件(污染物的物理化学参数，环境参数)，且对应模拟的年份，更符合实际情况；加入平流规则，污染物在网格间的迁移通过平流规则实现，模拟结果更接近实际情况；使用四级逸度模型，以月为一个模拟阶段，相对于现有的方法(以年为单位)，更接近实际情况，以解决现有技术中存在的问题。为了实现上述目的，本专利技术实施例提供如下技术方案：一种模拟有机污染物时空迁移归趋分布的方法，包括以下步骤：S1、建立数据库，包括：S1.1、首先，进行区域网格尺度的确定以及对各网格参数进行收集，然后通过收集的数据，得到污染物物理化学参数数据库、排放产量数据库、排放因子数据库、环境参数数据库、风速参数数据库、土地利用参数数据库、网格划分数据库、季节阶段数据库、温度和降雨速率数据库和不确定性分析数据库；S1.2、由上述步骤中得到的数据库组建整个模型数据库；S2、对研究区域进行网格划分；S3、建立四级逸度模型：整合逸度计算公式、Z值计算公式、D值计算公式、平流计算公式，将这些公式通过Matlab编程平台进行程序化；S4、采用四级逸度模型并配合平流规则进行四级多介质逸度模型模拟，包括：S4.1、调用已经编好的程序在Matlab上运行，得到不同网格在环境多介质中的浓度及迁移通量；S4.2、通过灵敏度、不确定性和拟合分析进行模型验证，分析参数有较大误差后，进行参数校正，优化模型；当分析的参数对比无较大误差，则进行污染物时空迁移、归趋特征、赋存浓度预测及暴露风险评估。进一步地，步骤S2中所述的网格划分，具体为利用arcgis将研究区域进行网格划分。进一步地，步骤S3中所述的整合逸度计算公式包括五个环境相的计算公式：大气：水体：土壤：V3Z3dF3/dt＝E3+D13F1+D53F3-(D31+D35+DR3)F3沉积物：V4Z4dF4/dt＝D24F2-(D42+DR4)F4植物：V5Z5dF5/dt＝D15F1+D35F3-(D51+DR5)F5其中1，2，3，4，5分别代表大气，水体，土壤，沉积物，植物；E代表污染物在环境中的排放速率，单位为mol/h；V为环境相的体积，单位为m3；Z为污染物在环境相内的逸度容量，单位为mol/(m3·Pa)；F是污染物在环境相内的逸度，单位为Pa；D为污染物在各环境间的迁移速率，单位为mol/(Pa·h)；DR代表污染物在环境相的降解速率，单位为mol/(Pa·h)；∑D(j,i)A和DA分别代表污染物为从网格j平流进网格i和从网格i平流出去的速率，单位为mol/(Pa·h)。进一步地，步骤S3中所述的Z值计算公式包括大气、水体、土壤、沉积物和植物等五个环境相的计算公式，且大气环境相包括气、颗粒相和总相共三个子环境相，水体环境相包括水、悬浮颗粒相和总相共三个子环境相，土壤环境相包括固体和总相两个子环境相，沉积物包括固体和总相两个子环境相，植物包括叶表皮相和总相两个子环境相；其中，每个子环境相对应一个Z值计算公式，具体包括：大气环境相中的气：Z1＝1/RT；大气环境相中的颗粒相：ZQ＝O13×Koa×ρq×Z1×0.001；大气环境相中的总相：ZA＝Z1×X11+ZQ×X13；水体环境相中的水：Z2＝1/H；水体环境相中的悬浮颗粒相：ZP＝O23×Koc×ρp×Z2×0.001；水体环境相中的总相：ZW＝Z2×X22+ZP×X23；土壤环境相中的固体：Z3＝O33×ρ3×Z2×Koc×0.001；土壤环境相中的总相：ZS＝Z1×X31+Z2×X32+Z3×X33；土壤环境相中的固体：Z4＝O43×ρ4×Z2×Koc×0.001；土壤本文档来自技高网...

【技术保护点】
1.一种模拟有机污染物时空迁移归趋分布的方法，其特征在于，包括以下步骤：/nS1、建立数据库，包括：/nS1.1、首先，进行区域网格尺度的确定以及对各网格参数进行收集，然后通过收集的数据，得到污染物物理化学参数数据库、排放产量数据库、排放因子数据库、环境参数数据库、风速参数数据库、土地利用参数数据库、网格划分数据库、季节阶段数据库、温度和降雨速率数据库和不确定性分析数据库；/nS1.2、由上述步骤中得到的数据库组建整个模型数据库；/nS2、对研究区域进行网格划分；/nS3、建立四级逸度模型：整合逸度计算公式、Z值计算公式、D值计算公式、平流计算公式，将这些公式通过Matlab编程平台进行程序化；/nS4、采用四级逸度模型并配合平流规则进行四级多介质逸度模型模拟，包括：/nS4.1、调用已经编好的程序在Matlab上运行，得到不同网格在环境多介质中的浓度及迁移通量；/nS4.2、通过灵敏度、不确定性和拟合分析进行模型验证，分析参数有较大误差后，进行参数校正，优化模型；当分析的参数对比无较大误差，则进行污染物时空迁移、归趋特征、赋存浓度预测及暴露风险评估。/n

【技术特征摘要】
1.一种模拟有机污染物时空迁移归趋分布的方法，其特征在于，包括以下步骤：
S1、建立数据库，包括：
S1.1、首先，进行区域网格尺度的确定以及对各网格参数进行收集，然后通过收集的数据，得到污染物物理化学参数数据库、排放产量数据库、排放因子数据库、环境参数数据库、风速参数数据库、土地利用参数数据库、网格划分数据库、季节阶段数据库、温度和降雨速率数据库和不确定性分析数据库；
S1.2、由上述步骤中得到的数据库组建整个模型数据库；
S2、对研究区域进行网格划分；
S3、建立四级逸度模型：整合逸度计算公式、Z值计算公式、D值计算公式、平流计算公式，将这些公式通过Matlab编程平台进行程序化；
S4、采用四级逸度模型并配合平流规则进行四级多介质逸度模型模拟，包括：
S4.1、调用已经编好的程序在Matlab上运行，得到不同网格在环境多介质中的浓度及迁移通量；
S4.2、通过灵敏度、不确定性和拟合分析进行模型验证，分析参数有较大误差后，进行参数校正，优化模型；当分析的参数对比无较大误差，则进行污染物时空迁移、归趋特征、赋存浓度预测及暴露风险评估。


2.根据权利要求1所述的一种模拟有机污染物时空迁移归趋分布的方法，其特征在于：步骤S2中所述的网格划分，具体为利用arcgis将研究区域进行网格划分。


3.根据权利要求1所述的一种模拟有机污染物时空迁移归趋分布的方法，其特征在于：步骤S3中所述的整合逸度计算公式包括五个环境相的计算公式：
大气：
水体：
土壤：V3Z3dF3/dt＝E3+D13F1+D53F3-(D31+D35+DR3)F3
沉积物：V4Z4dF4/dt＝D24F2-(D42+DR4)F4
植物：V5Z5dF5/dt＝D15F1+D35F3-(D51+DR5)F5
其中1，2，3，4，5分别代表大气，水体，土壤，沉积物，植物；E代表污染物在环境中的排放速率，单位为mol/h；V为环境相的体积，单位为m3；Z为污染物在环境相内的逸度容量，单位为mol/(m3·Pa)；F是污染物在环境相内的逸度，单位为Pa；D为污染物在各环境间的迁移速率，单位为mol/(Pa·h)；DR代表污染物在环境相的降解速率，单位为mol/(Pa·h)；∑D(j,i)A和DA分别代表污染物为从网格j平流进网格i和从网格i平流出去的速率，单位为mol/(Pa·h)。


4.根据权利要求1所述的一种模拟有机污染物时空迁移归趋分布的方法，其特征在于：步骤S3中所述的Z值计算公式包括大气、水体、土壤、沉积物和植物等五个环境相的计算公式，且大气环境相包括气、颗粒相和总相共三个子环境相，水体环境相包括水、悬浮颗粒相和总相共三个子环境相，土壤环境相包括固体和总相两个子环境相，沉积物包括固体和总相两个子环境相，植物包括叶表皮相和总相两个子环境相；其中，每个子环境相对应一个Z值计算公式，具体包括：
大气环境相中的气：Z1＝1/RT；
大气环境相中的颗粒相：ZQ＝O13×Koa×ρq×Z1×0.001；
大气环境相中的总相：ZA＝Z1×X11+ZQ×X13；
水体环境相中的水：Z2＝1/H；
水体环境相中的悬浮颗粒相：ZP＝O23×Koc×ρp×Z2×0.001；
水体环境相中的总相：ZW＝Z2×X22+ZP×X23；
土壤环境相中的固体：Z3＝O33×ρ3×Z2×Koc×0.001；
土壤环境相中的总相：ZS＝Z1×X31+Z2×X32+Z3×X33；
土壤环境...

【专利技术属性】
技术研发人员：亦如瀚，陈铭聪，高梓闻，杨东升，连逸轩，刘骝，
申请(专利权)人：暨南大学，
类型：发明
国别省市：广东;44