[发明专利]一种模拟有机污染物时空迁移归趋分布的方法

权利要求书

1.一种模拟有机污染物时空迁移归趋分布的方法，其特征在于，包括以下步骤：

S1、建立数据库，包括：

S1.1、首先，进行区域网格尺度的确定以及对各网格参数进行收集，然后通过收集的数据，得到污染物物理化学参数数据库、排放产量数据库、排放因子数据库、环境参数数据库、风速参数数据库、土地利用参数数据库、网格划分数据库、季节阶段数据库、温度和降雨速率数据库和不确定性分析数据库；

S1.2、由上述步骤中得到的数据库组建整个模型数据库；

S2、对研究区域进行网格划分；

S3、建立四级逸度模型：整合逸度计算公式、Z值计算公式、D值计算公式、平流计算公式，将这些公式通过Matlab编程平台进行程序化；

S4、采用四级逸度模型并配合平流规则进行四级多介质逸度模型模拟，包括：

S4.1、调用已经编好的程序在Matlab上运行，得到不同网格在环境多介质中的浓度及迁移通量；

S4.2、通过灵敏度、不确定性和拟合分析进行模型验证，分析参数有较大误差后，进行参数校正，优化模型；当分析的参数对比无较大误差，则进行污染物时空迁移、归趋特征、赋存浓度预测及暴露风险评估；

其中，步骤S3中所述的整合逸度计算公式包括五个环境相的计算公式：

大气：

水体：

土壤：V₃Z₃dF₃/dt＝E₃+D₁₃F₁+D₅₃F₃-(D₃₁+D₃₅+D_R3)F₃

沉积物：V₄Z₄dF₄/dt＝D₂₄F₂-(D₄₂+D_R4)F₄

植物：V₅Z₅dF₅/dt＝D₁₅F₁+D₃₅F₃-(D₅₁+D_R5)F₅

其中1，2，3，4，5分别代表大气，水体，土壤，沉积物，植物；E代表污染物在环境中的排放速率，单位为mol/h；V为环境相的体积，单位为m³；Z为污染物在环境相内的逸度容量，单位为mol/(m³·Pa)；F是污染物在环境相内的逸度，单位为Pa；D为污染物在各环境间的迁移速率，单位为mol/(Pa·h)；D_R代表污染物在环境相的降解速率，单位为mol/(Pa·h)；∑D(j,i)_A和D_A分别代表污染物为从网格j平流进网格i和从网格i平流出去的速率，单位为mol/(Pa·h)；

步骤S3中所述的Z值计算公式包括大气、水体、土壤、沉积物和植物五个环境相的计算公式，且大气环境相包括气、颗粒相和总相共三个子环境相，水体环境相包括水、悬浮颗粒相和总相共三个子环境相，土壤环境相包括固体和总相两个子环境相，沉积物包括固体和总相两个子环境相，植物包括叶表皮相和总相两个子环境相。

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

3.根据权利要求1所述的一种模拟有机污染物时空迁移归趋分布的方法，其特征在于：步骤S4中所述的四级多介质逸度模型模拟，以月为一个模拟阶段。