环境优先控制污染物风险评估推荐值的确定方法与流程

1.本发明涉及环境与健康管理技术领域，尤其涉及一种环境优先控制污染物风险评估推荐值的确定方法。


背景技术：

2.环境优先控制污染物是环境管理和环境保护中的有效技术手段，对于环境优先控制污染物，尤其是已经有明确名单的环境介质，基于人体健康确定其在相应环境介质中的风险评估推荐值，有助于快速判断环境介质对于人体是否具有健康危害，同时对于环境健康管理工作及人群健康防护具有重要的意义。
3.我国目前对于环境优先控制污染物的筛选及确定具有大量的相关研究，但针对于其基于人体健康的风险控制值还少有报道。就目前生态环境部公开的技术规范和资料而言，仅公开发布了针对于地表水环境的《人体健康水质基准制定技术指南》(hj837-2017)和针对于土壤环境介质的《土壤环境质量建设用地土壤污染风险管控标准(试行)》(gb36600-2018)及《土壤环境质量农用地土壤污染风险管控标准(试行)》(gb 15618-2018)，对于水环境、大气环境等目前均没有成熟的风险筛选值和风险管控值供环境管理使用，且不同环境介质的基于保护对象的不同，制订方法也不同，其中hj837-2017和gb 36600-2018基于人体健康制订，而gb 15618-2018基于农产品安全制订。
4.通过资料调研发现，目前对于优先控制污染物基于人体健康的风险评估推荐值的确定存在以下几点问题：一是研究较少，大多针对于污染场地的修复，对于其他环境介质的研究较少，且没有相关完整的方法可在实际应用中参考使用；二是参数的使用没有较合理的原则，大多参数使用采用国外参数或国家统一参数，没有考虑实际的暴露接触情况和当地的实际环境情况，使得计算结果可能存在偏差；三是对于优控污染物来源的迁移转化没有纳入研究中，目前仅有的针对于污染场地的推荐值研究中，仅考虑到污染场地本身，并未考虑污染场地污染物的来源，未考虑到其来源的风险控制值。


技术实现要素：

5.鉴于上述的分析，本发明实施例旨在提供一种环境优先控制污染物风险评估推荐值的确定，用以解决现有基于人体健康的风险评估推荐值的确定中涉及的环境介质系统框架不全面、未考虑污染物来源，以及参数选取存在偏差的问题。
6.本发明实施例提供了一种环境优先控制污染物风险评估推荐值的确定方法，包括如下步骤：根据选择的待确定风险评估推荐值的环境优先控制污染物，以及对应筛选的环境介质，提取暴露特征；
7.根据暴露特征和选择的健康效应，获取健康效应参数及参数值；
8.根据环境介质和暴露特征中的暴露途径，获取环境介质在每个暴露途径下的暴露计算模型，并根据暴露计算模型和选择的健康效应，获取风险控制值计算模型；
9.根据暴露计算模型和风险控制值计算模型，提取出除健康效应参数之外的暴露参
数并获取暴露参数的参数值；
10.根据健康效应参数和暴露参数的参数值，计算出环境介质的风险控制值，作为环境优先控制污染物的风险评估推荐值。
11.基于上述方法的进一步改进，方法还包括：
12.根据暴露特征中的排放来源，获取环境介质的每个排放来源的迁移转化修正系数；
13.将环境介质的风险控制值与对应的每个排放来源的迁移转化修正系数相乘，得到环境介质的污染源风险控制值，与环境介质的风险控制值一起作为环境优先控制污染物的风险评估推荐值。
14.基于上述方法的进一步改进，环境介质包括空气环境、土壤环境和水体环境；其中，空气环境包括环境空气和室内空气，土壤环境包括建设用地和农用地，水体环境包括地表水和地下水；
15.暴露特征的提取范围包括：暴露途径、暴露人群、排放来源、用地类型和迁移途径；其中，暴露途径包括经消化道、经呼吸道和经皮肤接触，暴露人群包括一般人群、儿童和特殊人群，排放来源包括经空气排放、经水体排放和经固废排放，用地类型包括建设用地和农用地，迁移途径包括废气、废水和废渣。
16.基于上述方法的进一步改进，健康效应包括致癌效应和非致癌效应；
17.根据暴露特征和选择的健康效应，获取健康效应参数及参数值，包括：
18.如果选择的健康效应是致癌效应，则获取可接受致癌风险，以及，根据暴露特征中的暴露途径，获取在环境介质中各暴露途径致癌因子，将可接受致癌风险和各暴露途径致癌因子作为健康效应参数，将健康效应参数及参数值，用于致癌效应风险控制值计算模型；
19.如果选择的健康效应是非致癌效应，则获取可接受非致癌风险，以及，根据暴露特征中的暴露途径，获取在环境介质中各暴露途径参考剂量，作为健康效应参数，将健康效应参数及参数值，用于非致癌效应风险控制值计算模型；
20.如果选择的健康效应是致癌效应和非致癌效应，则同时获取上述两种情况的健康效应参数及参数值，分别用于致癌效应风险控制值和非致癌效应风险控制值计算模型。
21.基于上述方法的进一步改进，根据环境介质和暴露特征中的暴露途径，获取环境介质在每个暴露途径下的暴露计算模型，包括：
22.如果环境介质是空气环境，暴露途径是经呼吸道，得到环境空气中污染物经呼吸道暴露量和室内空气中污染物经呼吸道暴露量的暴露计算模型；
23.如果环境介质是水体环境，暴露途径是经消化道和经皮肤接触，得到饮用地表水时污染物经消化道暴露量、摄入地表水中水产品时污染物经消化道暴露量、使用地表水时污染物经皮肤接触暴露量、饮用地下水时污染物经消化道暴露量、使用地下水时污染物经皮肤接触暴露量的暴露计算模型；
24.如果环境介质是土壤环境，暴露途径是经消化道、经皮肤接触和经呼吸道，得到人群暴露于建设用地时污染物经消化道暴露量、人群暴露于农用地时污染物经消化道暴露量、建设用地/农用地人群非气态污染物经呼吸道暴露量、建设用地/农用地人群表层土壤气态污染物室外经呼吸道暴露量、建设用地/农用地人群下层土壤气态污染物室外经呼吸道暴露量、建设用地人群地下层土壤气态污染物室内经呼吸道暴露量和人群建设用地/农
用地污染物经皮肤接触暴露量的暴露计算模型。
25.基于上述方法的进一步改进，根据暴露计算模型和选择的健康效应，获取风险控制值计算模型包括：
26.基于空气环境中的多种暴露情景，根据对应的多个暴露计算模型，得到选择的健康效应的风险控制值计算模型；
27.基于水体环境中每种水体类型的多种暴露途径，根据对应的多个暴露计算模型，以及各暴露途径下暴露量的叠加，得到选择的健康效应的综合风险控制值计算模型；
28.基于土壤环境中每种用地类型的多种暴露途径，根据对应的多个暴露计算模型，以及各暴露途径下暴露量的叠加，得到选择的健康效应的综合风险控制值计算模型。
29.基于上述方法的进一步改进，风险控制值计算模型中包括人群在每种暴露途径和每种暴露途径的每种方式下的暴露贡献比。
30.基于上述方法的进一步改进，获取暴露参数的参数值，包括：根据数据来源和数据范围识别暴露参数的参数值是否可使用，如果不可使用，则通过蒙特卡洛计算方法得到暴露参数的参数值；
31.数据来源的可使用条件包括：来源于国家或地区大型调查或具有严格质量控制措施的结果，并且，来源于近10年内的研究结果，并且，来源于权威性的公开发布的数据；
32.数据范围的可使用条件包括：数据所属区域范围在所选择的计算区域范围或相邻区域范围内。
33.基于上述方法的进一步改进，识别暴露参数的参数值是否可使用，包括：
34.如果暴露参数的参数值具有国家或地区推荐值，且同时满足数据来源和数据范围的可使用条件，则暴露参数的参数值可使用；
35.如果暴露参数的参数值具有国家或地区推荐值，但不能同时满足数据来源和数据范围的可使用条件，则属于不可使用，将暴露参数加入待计算集合；
36.如果暴露参数的参数值不具有国家或地区推荐值，但同时满足数据来源和数据范围的可使用条件，则属于不可使用，将暴露参数加入待计算集合；
37.如果暴露参数无参数值，或者参数值不具有国家或地区推荐值，且不满足数据来源或数据范围的可使用条件，则属于异常数据。
38.基于上述方法的进一步改进，通过蒙特卡洛计算方法得到暴露参数的参数值，包括：
39.获取待计算集合中暴露参数的原始数据，并识别原始数据是否具有数据分布特征，如果不具有数据分布特征，根据预置的分布类型，依次校验原始数据是否符合该分布类型，最终确定出数据分布特征；
40.根据直接得到或确定出来的数据分布特征，进行预置次数的运算模拟，将模拟运算后的结果，作为暴露参数的最终参数值。
41.与现有技术相比，本发明至少可实现如下有益效果之一：
42.1、基于对人群健康的保护，针对环境优先控制污染物提出了各环境介质中的风险评估推荐值的具体确定方法，在方法框架中考虑到计算参数的使用，嵌入了蒙特卡洛计算方法，获取计算所需参数值的最优值，提高了风险评估推荐值的准确度与可靠性；
43.2、在计算模块中考虑到计算优先污染物来源的控制值的情况，提出了采用污染物
迁移转化系数修正的方法来获取污染源的风险控制值，可同时获取环境介质和污染源中环境优先控制污染物的风险控制值；
44.3、同时给出了几种不同环境介质的风险控制值的计算方法，对于人群环境优先控制污染物的健康风险控制值的计算提供了便捷，可在短时间内获取不同环境介质中的控制值，对于人群环境健康管理非常有用，在提供参考和支撑的同时，为环境优先控制污染物的控制提供借鉴，且可为保障人群健康给出科学指导。
45.本发明中，上述各技术方案之间还可以相互组合，以实现更多的优选组合方案。本发明的其他特征和优点将在随后的说明书中阐述，并且，部分优点可从说明书中变得显而易见，或者通过实施本发明而了解。本发明的目的和其他优点可通过说明书以及附图中所特别指出的内容中来实现和获得。
附图说明
46.附图仅用于示出具体实施例的目的，而并不认为是对本发明的限制，在整个附图中，相同的参考符号表示相同的部件。
47.图1为本发明实施例中环境优先控制污染物风险评估推荐值的确定方法的流程图；
48.图2为本发明实施例中数据库的构建流程图；
49.图3为本发明实施例中环境介质筛选与暴露特征提取流程图；
50.图4为本发明实施例中空气环境与水体环境暴露计算模型及风险控制值计算模型对应选择的流程图；
51.图5为本发明实施例中土壤环境暴露计算模型及风险控制值计算模型对应选择的流程图。
具体实施方式
52.下面结合附图来具体描述本发明的优选实施例，其中，附图构成本技术一部分，并与本发明的实施例一起用于阐释本发明的原理，并非用于限定本发明的范围。
53.本发明的一个具体实施例，公开了一种环境优先控制污染物风险评估推荐值的确定方法，如图1所示，包括如下步骤：
54.s10：收集数据，构建数据库；
55.需要说明的是，在确定风险评估推荐值的过程中需要基于相关数据进行选择和运算，首先需要针对可能用到的参数数据进行资料收集，收集内容如图2所示，包括优控污染物相关、企业及排放数据、污染物迁移转化数据、污染物性质数据、人群暴露数据。其中，优控污染物相关资料收集包括：针对于空气、水体、土壤等环境介质或化学品中物质的相关优先控制污染物名录，以及相关环境背景值或环境介质质量浓度限值；企业及排放数据资料收集包括：环境优先控制污染物可能的排放企业类型、企业位置分布、污染物排放类型及污染物的排放量等数据资料；污染物迁移转化数据资料收集包括：污染物的迁移途径及其相关迁移参数等；污染物性质数据资料收集包括：污染物毒性参数及其理化性质；人群暴露数据资料收集包括：基础的暴露参数，比如身体特征参数、摄入量参数、行为模式参数等常规暴露参数，以及包含各介质各途径暴露接触贡献比的总暴露数据。资料收集过程中以收集
本国相关结果为主，如本国数据缺失，则链接其他国家数据，示例性地，污染物毒性数据可链接美国epa相关数据库。
56.数据资料收集完成后，将数据导入分类存储单元，进行分类存储，最终构建优控污染物数据库、污染物来源及排放数据库、污染物迁移数据库、污染物性质数据库和暴露参数数据库，存储于数据中心，用于后续选择和计算的调用。
57.s11：根据选择的待确定风险评估推荐值的环境优先控制污染物，以及对应筛选的环境介质，提取暴露特征；
58.需要说明的是，基于已建立的优控污染物数据库，首先选择待确定风险评估推荐值的环境优先控制污染物，然后对选择的优先控制污染物风险控制值筛选环境介质，包括空气环境、土壤环境、水体环境。其中，空气环境包括环境空气和室内空气，土壤环境包括建设用地和农用地，水体环境包括地表水和地下水。
59.如果选择的环境优先控制污染物是化学品中的，则可以根据化学品中优先控制污染物的赋存介质选择空气环境、土壤环境和水体环境中的1种或多种介质开展计算。
60.基于污染物迁移数据库，暴露特征的提取范围包括：暴露途径、暴露人群、排放来源、用地类型和迁移途径；其中，暴露途径包括经消化道、经呼吸道和经皮肤接触，暴露人群包括一般人群、儿童和特殊人群，排放来源包括经空气排放、经水体排放和经固废排放，用地类型包括建设用地和农用地，迁移途径包括废气、废水和废渣。
61.需要说明的是，暴露人群指需要进行计算的人群，可以根据需求从中进行选择，暴露途径和环境介质用于确定暴露计算模型，用地类型仅在选择土壤环境时使用，排放来源和迁移途径在计算污染源风险控制值时使用。
62.提取暴露特征时，基于环境介质与暴露特征的关联关系，预置为筛选规则，根据选择的环境介质，自动输出暴露特征。
63.通过图3说明该筛选规则，在图中，空气环境标记为
①
，土壤环境标记为
②
，水体环境标记为
③
，根据图中连线上
①②③
的标记，对应得到每种选择的环境介质的暴露特征。需要说明的是，对污染物是化学品的，如选择的是单一的环境介质，则根据单一环境介质对应的标记得到暴露特征；如选择的是多种环境介质，则根据选择的每一个单一的环境介质对应的标记，得到暴露特征后合并作为最终提取的暴露特征。
64.具体来说，对空气环境提取的暴露特征包括：经呼吸道(暴露途径)，一般人群、儿童和特殊人群(暴露人群)，经空气排放(排放来源)，废气(迁移途径)；对土壤环境提取的暴露特征包括：经消化道、经呼吸道和经皮肤接触(暴露途径)，一般人群、儿童和特殊人群(暴露人群)，经水体排放和经固废排放(排放来源)，建设用地和农用地(用地类型)，废气、废水和废渣(迁移途径)；对水体环境提取的暴露特征包括：经消化道和经皮肤接触(暴露途径)，一般人群、儿童和特殊人群(暴露人群)，经水体排放和经固废排放(排放来源)，废水和废渣(迁移途径)。
65.s12：根据暴露特征和选择的健康效应，获取健康效应参数及参数值；
66.需要说明的是，健康效应包括致癌效应和非致癌效应；根据步骤s11中选择的待确定风险评估推荐值的环境优先控制污染物，调用污染物性质数据库，获取对应的健康效应参数及参数值。
67.具体来说，根据暴露特征和选择的健康效应，获取健康效应参数及参数值，包括：
68.如果选择的健康效应是致癌效应，则获取可接受致癌风险，以及，根据暴露特征中的暴露途径，获取在环境介质中各暴露途径致癌因子，将可接受致癌风险和各暴露途径致癌因子作为健康效应参数，将健康效应参数及参数值，用于致癌效应风险控制值计算模型；
69.示例性地，环境介质是空气环境，暴露途径是经呼吸道，则获取呼吸吸入致癌斜率因子。
70.如果选择的健康效应是非致癌效应，则获取可接受非致癌风险，以及，根据暴露特征中的暴露途径，获取在环境介质中各暴露途径参考剂量，作为健康效应参数，将健康效应参数及参数值，用于非致癌效应风险控制值计算模型；
71.示例性地，环境介质是空气环境，暴露途径是经呼吸道，则获取呼吸吸入参考剂量。
72.如果选择的健康效应是致癌效应和非致癌效应，则同时获取上述两种情况的健康效应参数及参数值，分别用于致癌效应风险控制值和非致癌效应风险控制值计算模型。
73.s13：根据环境介质和暴露特征中的暴露途径，获取环境介质在每个暴露途径下的暴露计算模型，并根据暴露计算模型和选择的健康效应，获取风险控制值计算模型；
74.需要说明的是，暴露计算模型按照环境介质分类，以不同暴露途径为计算方式，实际运算过程中针对不同暴露人群，对各类具有人群特征差异的参数，采用各类人群的参数分别进行计算，模型公式中不再体现不同的暴露人群。
75.具体来说，每种环境介质在每个暴露途径下的暴露计算模型，包括：
76.s131：如果环境介质是空气环境，暴露途径是经呼吸道，得到环境空气中污染物经呼吸道暴露量和室内空气中污染物经呼吸道暴露量的暴露计算模型；
77.需要说明的是，环境空气(即室外空气)和室内空气是两种暴露情景。
78.a)环境空气中污染物经呼吸道暴露量的暴露计算模型，表示人群在室外暴露于环境空气中污染物的情况，计算公式是：
[0079][0080]
其中，ajer
out
为环境空气中污染物经呼吸道暴露量，m3/(kg
·
d)；ir为呼吸量，m3/h；et
out
为室外活动时间，h/d；ef
out
为室外环境暴露频率，d/a；ed为暴露持续时间，a；bw为体重，kg；at为平均暴露时间，d。
[0081]
需要说明的是，h表示小时，d表示日，a表示年。
[0082]
b)室内空气中污染物经呼吸道暴露量的暴露计算模型，表示人群在室内环境中暴露于室内空气中污染物的情况，计算公式是：
[0083][0084]
其中，aier
in
为室内空气中污染物经呼吸道暴露量，m3/(kg
·
d)；ir为呼吸量，m3/h；et
in
为室内活动时间，h/d；ef
in
为室内环境暴露频率，d/a；ed为暴露持续时间，a；bw为体重，kg；at为平均暴露时间，d。
[0085]
s132：如果环境介质是水体环境，暴露途径是经消化道和经皮肤接触，得到饮用地表水时污染物经消化道暴露量、摄入地表水中水产品时污染物经消化道暴露量、使用地表水时污染物经皮肤接触暴露量、饮用地下水时污染物经消化道暴露量、使用地下水时污染
物经皮肤接触暴露量的暴露计算模型；
[0086]
需要说明的是，水体环境的水体类型包括地表水和地下水，其中，地表水涉及到的主要为人群直接暴露为饮水、皮肤接触的用水暴露以及摄入水产品的间接暴露，地下水主要涉及的是人群饮水暴露及皮肤接触的用水暴露，下面分别进行说明。
[0087]
1)经消化道暴露
[0088]
a)饮用地表水时污染物经消化道暴露量的暴露计算模型，表示人群通过饮用和摄入地表水水产品暴露污染物的情况，计算公式是：
[0089][0090]
其中，swier
oral_d
为饮用地表水时污染物经消化道暴露量，l/(kg
·
d)；irw为饮水摄入量，l/d；ef为暴露频率，d/a；ed为暴露持续时间，a；bw为体重，kg；at为平均暴露时间，d。
[0091]
b)摄入地表水中水产品时污染物经消化道暴露量的暴露计算模型，表示人群通过摄入地表水水产品暴露污染物的情况，计算公式是：
[0092][0093]
其中，swier
oral_i
为摄入地表水中水产品时污染物经消化道暴露量，l/(kg
·
d)；ir
t
为第t营养级水产品摄入量(t为2、3和4级)，kg/d；baf
t
为第t营养级水产品的生物累积系数；ef为暴露频率，d/a；fi为吸收因子，无量纲；ed为暴露持续时间，a；bw为体重，kg；at为平均暴露时间，d。
[0094]
2)经皮肤接触暴露
[0095]
a)使用地表水时污染物经皮肤接触暴露量的暴露计算模型，表示人群通过用水皮肤接触地表水暴露污染物的情况，计算公式是：
[0096][0097]
其中，swier
dermal
为使用地表水时污染物经皮肤接触暴露量，l/(kg
·
d)；sa
w-dermal
为接触水的皮肤表面积，cm2；cf为体积转换因子，1l/1000cm3；et为暴露时间，h/d；ef为暴露频率，d/a；ed为暴露持续时间，a；bw为体重，kg；at为平均暴露时间，d。
[0098]
b)饮用地下水时污染物经消化道暴露量的暴露计算模型，表示人群通过饮用地下水暴露污染物的情况，计算公式是：
[0099][0100]
其中，gwier
oral_d
为饮用地下水时污染物经消化道暴露量，l/(kg
·
d)；irw为饮水摄入量，l/d；ef为暴露频率，d/a；ed为暴露持续时间，a；bw为体重，kg；at为平均暴露时间，d。
[0101]
c)使用地下水时污染物经皮肤接触暴露量的暴露计算模型，表示人群通过用水皮肤接触地下水暴露污染物的情况，计算公式是：
[0102][0103]
其中，gwier
dermal
为使用地下水时污染物经皮肤接触暴露量，l/(kg
·
d)；sa
w-dermal
为接触水的皮肤表面积，cm2；pc为皮肤渗透常数，cm/h；cf为体积转换因子，1l/1000cm3；et为暴露时间，h/d；ef为暴露频率，d/a；ed为暴露持续时间，a；bw为体重，kg；at为平均暴露时间，d。
[0104]
s133：如果环境介质是土壤环境，暴露途径是经消化道、经皮肤接触和经呼吸道，得到人群暴露于建设用地时污染物经消化道暴露量、人群暴露于农用地时污染物经消化道暴露量、建设用地/农用地人群非气态污染物经呼吸道暴露量、建设用地/农用地人群表层土壤气态污染物室外经呼吸道暴露量、建设用地/农用地人群下层土壤气态污染物室外经呼吸道暴露量、建设用地人群地下层土壤气态污染物室内经呼吸道暴露量和人群建设用地/农用地污染物经皮肤接触暴露量的暴露计算模型；
[0105]
需要说明的是，土壤环境涉及到的土壤主要分为两大类，一类是建设用地，按照《土地利用现状分类》(gb/t 21010-2017)分类及《土壤环境质量建设用地土壤污染风险管控标准(试行)》(gb 36600-2018)将建设用地分为住宅/公共用地、商业服务用地以及工业用地三大类；另一类是农用地。不同用地类型下，人群暴露途径不同。
[0106]
1)经消化道暴露
[0107]
a)人群暴露于建设用地时污染物经消化道暴露量的暴露计算模型，表示建设用地人群经消化道暴露的情况，计算公式是：
[0108][0109]
其中，soierb为人群暴露于建设用地时污染物经消化道暴露量，kg/(kg
·
d)；i为人群第i类时期，包括儿童期和成人期；ir
s_i
为人群第i类时期土壤摄入量，mg/d；cfs为转换因子，10-6
kg/mg；fi为吸收因子，无量纲；efi为人群第i类时期暴露频率，d/a；edi为人群第i类时期暴露持续时间，a；bwi为人群第i类时期体重，kg；ati为人群第i类时期平均暴露时间，d。
[0110]
需要说明的是，人群第i类时期与实际运算时用地类型和针对的暴露人群有关，根据人群与不同用地类型的暴露接触方式或不同的活动模式来选择需要计算的暴露期。示例性地，如果是建设用地中的工业用地，由于成人在该类型用地上的暴露周期长、暴露频率高，则该用地类型下的暴露人群主要针对成人，则只需要计算成人期的暴露量；如果是建设用地中的住宅/公共用地，则同时计算儿童期和成人期的暴露量。
[0111]
b)人群暴露于农用地时污染物经消化道暴露量的暴露计算模型，表示农用地人群经消化道暴露的情况，计算公式是：
[0112][0113]
其中，soiera为人群暴露于农用地时污染物经消化道暴露量，kg/(kg
·
d)；i为人群第i类时期，包括儿童期和成人期；ir
s_i
为人群第i类时期土壤摄入量，mg/d；ir
sf_i
为人群第i类时期农产品摄入量，mg/d；b
vi
为人群第i类时期摄入农产品中污染物由土壤至农作物的转换因子，kg/kg；fi为吸收因子，无量纲；cfs为转换因子，10-6
kg/mg；efi为人群第i类时期暴露频率，d/a；edi为人群第i类时期暴露持续时间，a；bwi为人群第i类时期体重，kg；ati为人群第i类时期平均暴露时间，d。
[0114]
2)经呼吸道接触
[0115]
a)建设用地/农用地人群非气态污染物经呼吸道暴露量的暴露计算模型，表示建设用地/农用地人群暴露于非气态环境优先控制污染物经呼吸道暴露的情况，计算公式是：
[0116][0117]
其中，saier
pb/pa
表示建设用地人群非气态污染物经呼吸道暴露量saier
pb
和农用地人群非气态污染物经呼吸道暴露量saier
pa
，kg/(kg
·
d)；i为人群第i类时期，包括儿童期和成人期；iri为人群第i类时期呼吸量，m3/d；edi为人群第i类时期暴露持续时间，a；ef
out_i
为人群第i类时期室外暴露频率，d/a；ef
in_i
为人群第i类时期室内暴露频率，d/a；f
out
为室外空气中来自土壤的颗粒物所占比例，无量纲；f
in
为室内空气中来自土壤的颗粒物所占比例，无量纲；bwi为人群第i类时期体重，kg；ati为人群第i类时期平均暴露时间，d；pef为起尘因子，m3/kg。
[0118]
b)建设用地/农用地人群室外表层土壤气态污染物经呼吸道暴露量的暴露计算模型，表示建设用地/农用地人群暴露于室外表层土壤气态环境优先控制污染物经呼吸道暴露的情况，计算公式是：
[0119][0120]
其中，saier
gosb/gosa
为建设用地人群室外表层土壤气态污染物经呼吸道暴露量saier
gosb
或农用地人群室外表层土壤气态污染物经呼吸道暴露量saier
gosa
，kg/(kg
·
d)；i为人群第i类时期，包括儿童期和成人期；iri为人群第i类时期呼吸量，m3/d；edi为人群第i类时期暴露持续时间，a；ef
out_i
为人群第i类时期室外暴露频率，d/a；vf
out_s
为表层土壤中污染物扩散进入室外空气的挥发因子，kg/m3；bwi为人群第i类时期体重，kg；ati为人群第i类时期平均暴露时间，d。
[0121]
c)建设用地/农用地人群室外下层土壤气态污染物经呼吸道暴露量的暴露计算模型，表示建设用地/农用地人群暴露于室外下层土壤气态环境优先控制污染物经呼吸道暴露的情况，计算公式是：
[0122][0123]
其中，saier
gogb/goga
为建设用地人群室外下层土壤气态污染物经呼吸道暴露量saier
gogb
或农用地人群室外下层土壤气态污染物经呼吸道暴露量saier
goga
，kg/(kg
·
d)；i为人群第i类时期，包括儿童期和成人期；iri为人群第i类时期呼吸量，m3/d；edi为人群第i类时期暴露持续时间，a；ef
out_i
为人群第i类时期室外暴露频率，d/a；vf
out_g
为下层土壤中污染物扩散进入室外空气的挥发因子，kg/m3；bwi为人群第i类时期体重，kg；ati为人群第i类时期平均暴露时间，d。
[0124]
d)建设用地人群室内下层土壤气态污染物经呼吸道暴露量的暴露计算模型，表示建设用地人群暴露于室内下层土壤气态环境优先控制污染物经呼吸道暴露的情况，计算公式是：
[0125][0126]
其中，saier
gigb
为建设用地人群室内下层土壤气态污染物经呼吸道暴露量，kg/
(kg
·
d)；i为人群第i类时期，包括儿童期和成人期；iri为人群第i类时期呼吸量，m3/d；edi为人群第i类时期暴露持续时间，a；ef
in_i
为人群第i类时期室内暴露频率，d/a；vf
in_g
为下层土壤中污染物扩散进入室外空气的挥发因子，kg/m3；bwi为人群第i类时期体重，kg；ati为人群第i类时期平均暴露时间，d。
[0127]
3)经皮肤接触
[0128]
人群建设用地/农用地污染物经皮肤接触暴露量的暴露计算模型，表示建设用地/农用地人群暴露于土壤环境优先控制污染物经皮肤接触暴露的情况，计算公式是：
[0129][0130]
其中，sdcer
b/a
为人群建设用地污染物经皮肤接触暴露量sdcerb或农用地污染物经皮肤接触暴露量sdcera，kg/(kg
·
d)；i为人群第i类时期，包括儿童期和成人期；sai为人群第i类时期接触土壤皮肤表面积，cm2/event；afi为人群第i类时期皮肤对土壤的黏附因子，mg/cm2；edi为人群第i类时期暴露持续时间，a；efi为人群第i类时期暴露频率，d/a；evi为人群第i类时期每日皮肤接触土壤频率，event/d；absd为皮肤吸收系数，无量纲；bwi为人群第i类时期体重，kg；cfd为转换因子，10-6
kg/mg；ati为人群第i类时期平均暴露时间，d。
[0131]
s134：如果选择的环境优先控制污染物是化学品中的，则可以根据化学品中优先控制污染物的赋存介质选择空气环境、土壤环境和水体环境中的1种或多种介质开展相关计算，包括上述对应的暴露计算模型。
[0132]
与现有技术相比，明确了不同环境介质中优先控制污染物风险控制值计算过程中的暴露量计算内容，并给出了完整的计算方法，丰富完善了计算系统框架。
[0133]
通过步骤s131-s134，得到了暴露计算模型，再根据选择的健康效应，获取风险控制值计算模型，包括两类：致癌健康效应的风险控制值计算模型和非致癌健康效应的风险控制值计算模型。
[0134]
根据每种环境介质的每个暴露计算模型，得到选择的健康效应的风险控制值计算模型，风险控制值计算模型包括以下三个维度：
[0135]
基于空气环境中的多种暴露情景，根据对应的多个暴露计算模型，得到选择的健康效应的风险控制值计算模型；
[0136]
基于水体环境中每种水体类型的多种暴露途径，根据对应的多个暴露计算模型，以及各暴露途径下暴露量的叠加，得到选择的健康效应的综合风险控制值计算模型；
[0137]
基于土壤环境中每种用地类型的多种暴露途径，根据对应的多个暴露计算模型，以及各暴露途径下暴露量的叠加，得到选择的健康效应的综合风险控制值计算模型。
[0138]
具体地，为了清晰展示出不同健康效应下，不同环境介质的不同暴露计算模型对应的风险控制值计算模型，通过图4和图5，结合下面的计算公式进行说明。
[0139]
1)致癌健康效应的风险控制值计算模型
[0140]
1.1)空气环境，包括图4中由a1-1、a1-2和a2-1、a2-2路径选择对应的风险控制值计算模型。
[0141]
·
a1-1：基于经呼吸道途径致癌效应的室外空气人群健康风险控制值计算模型，计算公式是：
[0142][0143]
其中，cal
out
为基于经呼吸道途径致癌效应的室外空气人群健康风险控制值，μg/m3；acr为可接受致癌风险，无量纲；aier
out
为环境空气中污染物经呼吸道暴露量，m3/(kg
·
d)；rsc
aout
为人群室外空气暴露贡献比；sf
inh
为呼吸吸入致癌斜率因子，(kg
·
d)/mg，ir为呼吸量，m3/d；iur为呼吸吸入单位致癌因子，m3/mg；bw为体重，kg。
[0144]
·
a2-1：基于经呼吸道途径致癌效应的室内空气人群健康风险控制值计算模型，计算公式是：
[0145][0146]
其中，cal
in
为基于经呼吸道途径致癌效应的室内空气人群健康风险控制值，μg/m3；acr为可接受致癌风险，无量纲；aier
in
为室内空气中污染物经呼吸道暴露量，m3/(kg
·
d)；rsc
ain
为人群室内空气暴露贡献比；sf
inh
为呼吸吸入致癌斜率因子，(kg
·
d)/mg，ir为呼吸量，m3/d；iur为呼吸吸入单位致癌因子，m3/mg；bw为体重，kg。
[0147]
1.2)水体环境，包括图4中由b1-11、b1-21、b1-31、b2-11、b2-21、b2-31、b3-11和b3-21路径选择所对应的地表水致癌效应风险控制值计算模型，以及，由c1-11、c1-21、c2-11和c2-21路径选择所对应的地下水致癌效应风险控制值计算模型。
[0148]
·
b1-11：仅考虑地表水饮用时，基于饮用经消化道途径致癌效应的地表水人群健康风险控制值计算模型，计算公式是：
[0149][0150]
其中，cswl
oral_d
为基于饮用经消化道途径致癌效应的地表水人群健康风险控制值，μg/l；acr为可接受致癌风险，无量纲；swier
oral_d
为饮用地表水时污染物经消化道暴露量，l/(kg
·
d)；rsc
sworal_d
为人群饮用地表水暴露贡献比；sf
oral
为经口摄入致癌斜率因子，(kg
·
d)/mg。
[0151]
·
b2-11：仅考虑摄入地表水中水产品时，基于摄入水产品经消化道途径致癌效应的人群健康风险控制值计算模型，计算公式是：
[0152][0153]
其中，cswl
oral_i
为基于摄入水产品经消化道途径致癌效应的人群健康风险控制值，μg/l；acr为可接受致癌风险，无量纲；swier
oral_i
为摄入地表水中水产品时污染物经消化道暴露量，l/(kg
·
d)；rsc
sworal_i
为人群摄入地表水中水产品时的暴露贡献比；sf
oral
为经口摄入致癌斜率因子，(kg
·
d)/mg。
[0154]
·
b1-21和b2-21：同时考虑饮用和摄入地表水水产品时，经消化道途径致癌效应的地表水人群健康风险控制值计算模型，计算公式是：
[0155][0156]
其中，cswl
oral
为基于饮用和摄入水产品经消化道途径致癌效应的地表水人群健
康风险控制值，μg/l；其他参数见式2-4和式2-5。
[0157]
·
b3-11：基于经皮肤接触途径致癌效应的地表水人群健康风险控制值计算模型，计算公式是：
[0158][0159]
其中，cswl
dermal
为基于经皮肤接触途径致癌效应的地表水人群健康风险控制值，μg/l；acr为可接受致癌风险，无量纲；swier
dermal
为使用地表水时污染物经皮肤接触暴露量，l/(kg
·
d)；rsc
swsdermal
为人群经皮肤接触地表水时的暴露贡献比；sf
dermal
为皮肤接触致癌斜率因子，(kg
·
d)/mg。
[0160]
·
b1-31、b2-31和b3-21：同时考虑饮用、摄入水产品及用水皮肤接触地表水时，基于2种地表水暴露途径综合致癌效应的地表水人群健康风险控制值计算模型，计算公式是：
[0161][0162]
其中，cswl为基于2种地表水暴露途径综合致癌效应的地表水人群健康风险控制值，μg/l；其他参数见式2-3至式2-6。
[0163]
·
c1-11：基于饮用经消化道途径致癌效应的地下水人群健康风险控制值计算模型，计算公式是：
[0164][0165]
其中，cgwl
oral_d
为基于饮用经消化道途径致癌效应的地下水人群健康风险控制值，μg/l；acr为可接受致癌风险，无量纲；gwier
oral_d
为饮用地下水时污染物经消化道暴露量，l/(kg
·
d)；rsc
gworal_d
为人群饮用地下水暴露贡献比；sf
oral
为经口摄入致癌斜率因子，(kg
·
d)/mg。
[0166]
·
c2-11：基于经皮肤接触途径致癌效应的地下水人群健康风险控制值计算模型，计算公式是：
[0167][0168]
其中，cgwl
dermal
为基于经皮肤接触途径致癌效应的地下水人群健康风险控制值，μg/l；acr为可接受致癌风险，无量纲；gwier
dermal
为使用地下水时污染物经皮肤接触暴露量，l/(kg
·
d)；rsc
gwdermal
为人群经皮肤接触地下水时的暴露贡献比；sf
dermal
为皮肤接触致癌斜率因子，(kg
·
d)/mg。
[0169]
·
c1-21和c2-21：同时考虑饮用及用水皮肤接触地下水时，基于2种地下水暴露途径综合致癌效应的地下水人群健康风险控制值计算模型，计算公式是：
[0170][0171]
其中，cgwl为基于2种地下水暴露途径综合致癌效应的地下水人群健康风险控制值，μg/l；其他参数见式2-8和式2-9。
[0172]
1.3)土壤环境，对于建设用地土壤，包括图5中由d1-11、d2-11、d2-21、d3-11、d3-21、d4-11、d4-21、d5-11、d5-21、d6-11、d6-21路径选择所对应的建设用地土壤致癌效应风
险控制值计算模型；对于农用地土壤，包括图5中由e1-11、e2-11、e2-21、e3-11、e3-21、e4-11、e4-21、e6-11、e6-21路径选择所对应的农用地土壤致癌效应风险控制值计算模型。
[0173]
·
d1-11：基于经消化道途径致癌效应的建设用地人群健康风险控制值计算模型，计算公式是：
[0174][0175]
其中，csl
oral_d
为基于经消化道途径致癌效应的建设用地人群健康风险控制值，mg/kg；acr为可接受致癌风险，无量纲；soierb为暴露于建设用地时污染物经消化道暴露量，kg/(kg
·
d)；rsc
soral_b
为人群建设用地经消化道暴露贡献比；sf
oral
为经口摄入致癌斜率因子，(kg
·
d)/mg。
[0176]
·
e1-11：基于经消化道途径致癌效应的农用地人群健康风险控制值计算模型，计算公式是：
[0177][0178]
其中，csl
oral_d
为基于经消化道途径致癌效应的农用地人群健康风险控制值，mg/kg；acr为可接受致癌风险，无量纲；soiera为暴露于农用地时污染物经消化道暴露量，kg/(kg
·
d)；rsc
soral_a
为人群农用地经消化道暴露贡献比；sf
oral
为经口摄入致癌斜率因子，(kg
·
d)/mg。
[0179]
·
d2-11：非气态污染物基于经呼吸道途径致癌效应的建设用地土壤人群健康风险控制值计算模型，计算公式是：
[0180][0181]
其中，csl
inh-pb
为非气态污染物基于经呼吸道途径致癌效应的建设用地土壤人群健康风险控制值，mg/kg；acr为可接受致癌风险，无量纲；saier
pb
为建设用地人群非气态污染物经呼吸道暴露量，kg/(kg
·
d)；rsc
sinh_b
为人群建设用地非气态污染物经呼吸道暴露贡献比；sf
inh
为呼吸吸入致癌斜率因子，(kg
·
d)/mg，ir为呼吸量，m3/d；iur为呼吸吸入单位致癌因子，m3/mg；bw为体重，kg。
[0182]
·
e2-11：非气态污染物基于经呼吸道途径致癌效应的农用地土壤人群健康风险控制值计算模型，计算公式是：
[0183][0184]
其中，csl
inh-pa
为非气态污染物基于经呼吸道途径致癌效应的农用地土壤人群健康风险控制值，mg/kg；acr为可接受致癌风险，无量纲；saier
pa
为农用地人群非气态污染物经呼吸道暴露量，kg/(kg
·
d)；rsc
sinh_a
为人群农用地非气态污染物经呼吸道暴露贡献比；sf
inh
为呼吸吸入致癌斜率因子，(kg
·
d)/mg，ir为呼吸量，m3/d；iur为呼吸吸入单位致癌因子，m3/mg；bw为体重，kg。
[0185]
·
d3-11、d4-11和d5-11：气态污染物基于经呼吸道途径致癌效应的建设用地土壤人群健康风险控制值计算模型，计算公式是：
[0186][0187]
其中，csl
inh-ab
为气态污染物基于经呼吸道途径致癌效应的建设用地土壤人群健康风险控制值，mg/kg；acr为可接受致癌风险，无量纲；saier
gosb
为建设用地或农用地人群表层土壤气态污染物室外经呼吸道暴露量，kg/(kg
·
d)；saier
gogb
为建设用地人群下层土壤气态污染物室外经呼吸道暴露量，kg/(kg
·
d)；saier
gigb
为建设用地人群地下层土壤气态污染物室内经呼吸道暴露量，kg/(kg
·
d)；rsc
sinh_gosb
为人群建设用地室外表层土壤气态污染物经呼吸道暴露贡献比；rsc
sinh_gogb
为人群建设用地室外下层土壤气态污染物经呼吸道暴露贡献比；rsc
sinh_gigb
为人群建设用地室内下层土壤气态污染物经呼吸道暴露贡献比；sf
inh
为呼吸吸入致癌斜率因子，(kg
·
d)/mg，
[0188]
·
e3-11和e4-11：气态污染物基于经呼吸道途径致癌效应的农用地土壤人群健康风险控制值计算模型，计算公式是：
[0189][0190]
其中，csl
inh-aa
为气态污染物基于经呼吸道途径致癌效应的农用地土壤人群健康风险控制值，mg/kg；acr为可接受致癌风险，无量纲；saier
gosa
为农用地人群表层土壤气态污染物室外经呼吸道暴露量，kg/(kg
·
d)；saier
goga
为农用地人群下层土壤气态污染物室外经呼吸道暴露量，kg/(kg
·
d)；rsc
sinh_gosa
为人群农用地室外表层土壤气态污染物经呼吸道暴露贡献比；rsc
sinh_goga
为人群农用地室外下层土壤气态污染物经呼吸道暴露贡献比；sf
inh
为呼吸吸入致癌斜率因子，(kg
·
d)/mg，d)/mg，
[0191]
·
d2-21、d3-21、d4-21和d5-21：优选地，对于建设用地环境优先控制污染物，考虑在非气态和气态中均存在的情况，其经呼吸道的致癌效应风险控制值计算模型，计算公式是：
[0192][0193]
其中，csl
inh-agb
为存在于非气态和气态中的污染物基于经呼吸道途径致癌效应的建设用地土壤人群健康风险控制值，mg/kg；其他参数见式2-13和式2-15。
[0194]
·
e2-21、e3-21和e4-21：优选地，对于农用地环境优先控制污染物，考虑在非气态和气态中均存在的情况，其经呼吸道的致癌效应风险控制值计算模型，计算公式是：
[0195][0196]
其中，csl
inh-aga
为为存在于非气态和气态中的污染物基于经呼吸道途径致癌效应的农用地土壤人群健康风险控制值，mg/kg；其他参数见式2-14和式2-16。
[0197]
·
d6-11：基于皮肤接触途径致癌效应的建设用地土壤人群健康风险控制值计算模型，计算公式是：
[0198]
[0199]
其中，csl
dermal-b
为基于皮肤接触途径致癌效应的建设用地土壤人群健康风险控制值，mg/kg；acr为可接受致癌风险，无量纲；sdcerb为人群建设用地污染物经皮肤接触暴露量，kg/(kg
·
d)；rsc
sder_b
为人群经皮肤接触建设用地土壤时的暴露贡献比；sf
dermal
为皮肤接触致癌斜率因子，(kg
·
d)/mg。
[0200]
·
e6-11：基于皮肤接触途径致癌效应的农用地土壤人群健康风险控制值计算模型，计算公式是：
[0201][0202]
其中，csl
dermal-a
为基于皮肤接触途径致癌效应的农用地土壤人群健康风险控制值，mg/kg；acr为可接受致癌风险，无量纲；sdcera为人群农用地污染物经皮肤接触暴露量，kg/(kg
·
d)；rsc
sder _a
为人群经皮肤接触农用地土壤时的暴露贡献比；sf
dermal
为皮肤接触致癌斜率因子，(kg
·
d)/mg。
[0203]
·
d6-21：优选地，对于建设用地基于全部暴露途径的综合致癌效应的人群健康风险控制值计算模型，计算公式是：
[0204][0205]
其中，cslb为基于全部暴露途径综合致癌效应的建设用地土壤人群健康风险控制值，mg/kg；其他参数见式2-11、式2-17和式2-19。
[0206]
·
e6-21：优选地，对于农用地基于全部暴露途径的综合致癌效应的人群健康风险控制值计算模型，计算公式是：
[0207][0208]
其中，csla为基于全部暴露途径综合致癌效应的农用地土壤人群健康风险控制值，mg/kg；其他参数见式2-12、式2-18和式2-20。
[0209]
2)非致癌健康效应的风险控制值计算模型
[0210]
2.1)空气环境，包括图4中由a1-2和a2-2路径选择对应的风险控制值的计算模型。
[0211]
·
a1-2：基于经呼吸道途径非致癌效应的室外空气人群健康风险控制值计算模型，计算公式是：
[0212][0213]
其中，ncal
out
为基于经呼吸道途径非致癌效应的室外空气人群健康风险控制值，μg/m3；ahq为可接受非致癌风险，无量纲；aier
out
为环境空气中污染物经呼吸道暴露量，m3/(kg
·
d)；rsc
aout
为人群室外空气暴露贡献比；rfd
inh
为呼吸吸入参考剂量，mg/(kg
·
d)，rfc为呼吸吸入参考浓度，mg/m3，ir为呼吸量，m3/d；bw为体重，kg。
[0214]
·
a2-2：基于经呼吸道途径非致癌效应的室内空气人群健康风险控制值计算模型，计算公式是：
[0215]
[0216]
其中，ncal
in
为基于经呼吸道途径非致癌效应的室内空气人群健康风险控制值，μg/m3；ahq为可接受非致癌风险，无量纲；aier
in
为室内空气中污染物经呼吸道暴露量，m3/(kg
·
d)；rsc
ain
为人群室内空气暴露贡献比；rfd
inh
为呼吸吸入参考剂量，mg/(kg
·
d)，rfc为呼吸吸入参考浓度，mg/m3，ir为呼吸量，m3/d；bw为体重，kg。
[0217]
2.2)水体环境，包括图4中由b1-12、b1-22、b1-32、b2-12、b2-22、b2-32、b3-12和b3-22路径选择所对应的地表水致癌效应风险控制值计算模型，以及，由c1-12、c1-22、c2-12和c2-22路径选择所对应的地下水致癌效应风险控制值计算模型。
[0218]
·
b1-12：仅考虑地表水饮用时，基于饮用经消化道途径非致癌效应的地表水人群健康风险控制值计算模型，计算公式是：
[0219][0220]
其中，ncswl
oral_d
为基于饮用经消化道途径非致癌效应的地表水人群健康风险控制值，μg/l；ahq为可接受非致癌风险，无量纲；swier
oral_d
为饮用地表水时污染物经消化道暴露量，l/(kg
·
d)；rsc
sworal_d
为人群饮用地表水暴露贡献比；rfd
oral
为经口摄入参考剂量，mg/(kg
·
d)。
[0221]
·
b2-12：仅考虑摄入地表水中水产品时，基于摄入水产品经消化道途径非致癌效应的人群健康风险控制值计算模型，计算公式是：
[0222][0223]
其中，ncswl
oral_i
为基于摄入水产品经消化道途径非致癌效应的人群健康风险控制值，μg/l；ahq为可接受非致癌风险，无量纲；swier
oral_i
为摄入地表水中水产品时污染物经消化道暴露量，l/(kg
·
d)；rsc
sworal_i
为人群摄入地表水中水产品时的暴露贡献比；rfd
oral
为经口摄入参考剂量，mg/(kg
·
d)。
[0224]
·
b1-22和b2-22：同时考虑饮用和摄入地表水水产品时，经消化道途径非致癌效应的地表水人群健康风险控制值计算模型，计算公式是：
[0225][0226]
其中，ncswl
oral
为基于饮用和摄入水产品经消化道途径非致癌效应的地表水人群健康风险控制值，μg/l；其他参数见式3-3和式3-4。
[0227]
·
b3-12：基于经皮肤接触途径非致癌效应的地表水人群健康风险控制值计算模型，计算公式是：
[0228][0229]
其中，ncswl
dermal
为基于经皮肤接触途径非致癌效应的地表水人群健康风险控制值，μg/l；ahq为可接受非致癌风险，无量纲；swier
dermal
为使用地表水时污染物经皮肤接触暴露量，l/(kg
·
d)；rsc
swsdermal
为人群经皮肤接触地表水时的暴露贡献比；rfd
dermal
为皮肤接触参考剂量，mg/(kg
·
d)。
[0230]
·
b1-32、b2-32和b3-22：同时考虑饮用、摄入水产品及用水皮肤接触地表水时，基
于2种地表水暴露途径综合非致癌效应的地表水人群健康风险控制值计算模型，计算公式是：
[0231][0232]
其中，ncswl为基于2种地表水暴露途径综合非致癌效应的地表水人群健康风险控制值，μg/l；其他参数式3-3见至式3-6。
[0233]
·
c1-12：基于饮用经消化道途径非致癌效应的地下水人群健康风险控制值计算模型，计算公式是：
[0234][0235]
其中，ncgwl
oral_d
为基于饮用经消化道途径非致癌效应的地下水人群健康风险控制值，μg/l；ahq为可接受非致癌风险，无量纲；gwier
oral_d
为饮用地下水时污染物经消化道暴露量，l/(kg
·
d)；rsc
gworal_d
为人群饮用地下水暴露贡献比；rfd
oral
为经口摄入参考剂量，mg/(kg
·
d)。
[0236]
·
c2-12：基于经皮肤接触途径非致癌效应的地下水人群健康风险控制值计算模型，计算公式是：
[0237][0238]
其中，ncgwl
dermal
为基于经皮肤接触途径非致癌效应的地下水人群健康风险控制值，μg/l；ahq为可接受非致癌风险，无量纲；gwier
dermal
为使用地下水时污染物经皮肤接触暴露量，l/(kg
·
d)；rsc
gwdermal
为人群经皮肤接触地下水时的暴露贡献比；rfd
dermal
为皮肤接触参考剂量，mg/(kg
·
d)。
[0239]
·
c1-22和c2-22：同时考虑饮用及用水皮肤接触地下水时，基于2种地下水暴露途径综合非致癌效应的地下水人群健康风险控制值计算模型，计算公式是：
[0240][0241]
其中，ncgwl为基于2种地下水暴露途径综合非致癌效应的地下水人群健康风险控制值，μg/l；其他参数见式3-8和式3-9。
[0242]
2.3)土壤环境，对于建设用地土壤，包括图5中由d1-12、d2-12、d2-22、d3-12、d3-22、d4-12、d4-22、d5-12、d5-22、d6-12、d6-22路径选择所对应的建设用地土壤非致癌效应风险控制值计算模型；对于农用地土壤，包括图5中由e1-12、e2-12、e2-22、e3-12、e3-22、e4-12、e4-22、e6-12、e6-22路径选择所对应的农用地土壤非致癌效应风险控制值计算模型。
[0243]
·
d1-12：基于经消化道途径非致癌效应的建设用地人群健康风险控制值计算模型，计算公式是：
[0244][0245]
其中，ncsl
orai_d
为基于经消化道途径非致癌效应的建设用地人群健康风险控制
值，mg/kg；ahq为可接受非致癌风险，无量纲；soierb为暴露于建设用地时污染物经消化道暴露量，kg/(kg
·
d)；rsc
soral_b
为人群建设用地经消化道暴露贡献比；rfd
oral
为经口摄入参考剂量，mg/(kg
·
d)。
[0246]
·
e1-12：基于经消化道途径非致癌效应的农用地人群健康风险控制值计算模型，计算公式是：
[0247][0248]
其中，ncsl
oral_a
为基于经消化道途径非致癌效应的农用地人群健康风险控制值，mg/kg；ahq为可接受非致癌风险，无量纲；soiera为暴露于农用地时污染物经消化道暴露量，kg/(kg
·
d)；rsc
soral_a
为人群农用地经消化道暴露贡献比；rfd
oral
为经口摄入参考剂量，mg/(kg
·
d)。
[0249]
·
d2-12：非气态污染物基于经呼吸道途径非致癌效应的建设用地土壤人群健康风险控制值计算模型，计算公式是：
[0250][0251]
其中，ncsl
inh-pb
为非气态污染物基于经呼吸道途径非致癌效应的建设用地土壤人群健康风险控制值，mg/kg；ahq为可接受非致癌风险，无量纲；saier
pb
为建设用地人群非气态污染物经呼吸道暴露量，kg/(kg
·
d)；rsc
sinh_b
为人群建设用地非气态污染物经呼吸道暴露贡献比；rfd
inh
为呼吸吸入参考剂量，mg/(kg
·
d)，rfc为呼吸吸入参考浓度，mg/m3，ir为呼吸量，m3/d；bw为体重，kg。
[0252]
·
e2-12：非气态污染物基于经呼吸道途径非致癌效应的农用地土壤人群健康风险控制值计算模型，计算公式是：
[0253][0254]
其中，ncsl
inh-pa
为非气态污染物基于经呼吸道途径非致癌效应的农用地土壤人群健康风险控制值，mg/kg；ahq为可接受非致癌风险，无量纲；saier
pa
为农用地人群非气态污染物经呼吸道暴露量，kg/(kg
·
d)；rsc
sinh_a
为人群农用地非气态污染物经呼吸道暴露贡献比；rfd
inh
为呼吸吸入参考剂量，mg/(kg
·
d)，rfc为呼吸吸入参考浓度，mg/m3，ir为呼吸量，m3/d；bw为体重，kg。
[0255]
·
d3-12、d4-12和d5-12：气态污染物基于经呼吸道途径非致癌效应的建设用地土壤人群健康风险控制值计算模型，计算公式是：
[0256][0257]
其中，ncsl
inh-ab
为气态污染物基于经呼吸道途径非致癌效应的建设用地土壤人群健康风险控制值，mg/kg；ahq为可接受非致癌风险，无量纲；saier
gosb
为建设用地或农用地人群表层土壤气态污染物室外经呼吸道暴露量，kg/(kg
·
d)；saier
gogb
为建设用地人群下层土壤气态污染物室外经呼吸道暴露量，kg/(kg
·
d)；saier
gigb
为建设用地人群地下层土
壤气态污染物室内经呼吸道暴露量，kg/(kg
·
d)；rsc
sinh_gosb
为人群建设用地室外表层土壤气态污染物经呼吸道暴露贡献比；rsc
sinh_gogb
为人群建设用地室外下层土壤气态污染物经呼吸道暴露贡献比；rsc
sinh_gigb
为人群建设用地室内下层土壤气态污染物经呼吸道暴露贡献比；rfd
inh
为呼吸吸入参考剂量，mg/(kg
·
d)。
[0258]
·
e3-12和e4-12：气态污染物基于经呼吸道途径非致癌效应的农用地土壤人群健康风险控制值计算模型，计算公式是：
[0259][0260]
其中，ncsl
inh-aa
为气态污染物基于经呼吸道途径非致癌效应的农用地土壤人群健康风险控制值，mg/kg；ahq为可接受非致癌风险，无量纲；saier
gosa
为农用地人群表层土壤气态污染物室外经呼吸道暴露量，kg/(kg
·
d)；saier
goga
为农用地人群下层土壤气态污染物室外经呼吸道暴露量，kg/(kg
·
d)；rsc
sinh_gosa
为人群农用地室外表层土壤气态污染物经呼吸道暴露贡献比；rsc
sinh_goga
为人群农用地室外下层土壤气态污染物经呼吸道暴露贡献比；rfd
inh
为呼吸吸入参考剂量，mg/(kg
·
d)。
[0261]
·
d2-22、d3-22、d4-22和d5-22：优选地，对于建设用地环境优先控制污染物，考虑在非气态和气态中均存在的情况，其经呼吸道的非致癌效应风险控制值计算模型，计算公式是：
[0262][0263]
其中，ncsl
inh-agb
为存在于非气态和气态中的污染物基于经呼吸道途径非致癌效应的建设用地土壤人群健康风险控制值，mg/kg；其他参数见式3-13和式3-15。
[0264]
·
e2-22、e3-22和e4-22：优选地，对于农用地环境优先控制污染物，考虑在非气态和气态中均存在的情况，其经呼吸道的非致癌效应风险控制值计算模型，计算公式是：
[0265][0266]
其中，ncsl
inh-aga
为为存在于非气态和气态中的污染物基于经呼吸道途径非致癌效应的农用地土壤人群健康风险控制值，mg/kg；其他参数见式3-14和式3-16。
[0267]
·
d6-12：基于皮肤接触途径非致癌效应的建设用地土壤人群健康风险控制值计算模型，计算公式是：
[0268][0269]
其中，ncsl
dermal-b
为基于皮肤接触途径非致癌效应的建设用地土壤人群健康风险控制值，mg/kg；ahq为可接受非致癌风险，无量纲；sdcerb为人群建设用地污染物经皮肤接触暴露量，kg/(kg
·
d)；rsc
sder_b
为人群经皮肤接触建设用地土壤时的暴露贡献比；rfd
dermal
为皮肤接触参考剂量，mg/(kg
·
d)。
[0270]
·
e6-12：基于皮肤接触途径非致癌效应的农用地土壤人群健康风险控制值计算模型，计算公式是：
[0271]
[0272]
其中，ncsl
dermal-a
为基于皮肤接触途径非致癌效应的农用地土壤人群健康风险控制值，mg/kg；ahq为可接受非致癌风险，无量纲；sdcera为人群农用地污染物经皮肤接触暴露量，kg/(kg
·
d)；rsc
sder_a
为人群经皮肤接触农用地土壤时的暴露贡献比；rfd
dermal
为皮肤接触参考剂量，mg/(kg
·
d)。
[0273]
·
d6-22：优选地，对于建设用地基于全部暴露途径的综合非致癌效应的人群健康风险控制值计算模型，计算公式是：
[0274][0275]
其中，ncslb为基于全部暴露途径综合非致癌效应的建设用地土壤人群健康风险控制值，mg/kg；其他参数见式3-11、式3-17和式3-19。
[0276]
·
e6-22：优选地，对于农用地基于全部暴露途径的综合非致癌效应的人群健康风险控制值计算模型，计算公式是：
[0277][0278]
其中，ncsla为基于全部暴露途径综合非致癌效应的农用地土壤人群健康风险控制值，mg/kg；其他参数见式3-12、式3-18和式3-20。
[0279]
与现有技术相比，风险控制值计算模型中包括人群在每种暴露途径和每种暴露途径的每种方式下的暴露贡献比，增加了实际的暴露接触情况，使风险评估推荐值更准确，有助于快速判断环境介质对于人体是否具有健康危害，同时降低了在计算过程中过度高估或者低估各环境介质或暴露途径带来的计算不确定性。
[0280]
s14：根据暴露计算模型和风险控制值计算模型，提取出除健康效应参数之外的暴露参数并获取暴露参数的参数值；
[0281]
需要说明的是，根据暴露计算模型和风险控制计算模型，可以得到模型中的参数，其中，健康效应参数，比如可接受致癌风险(acr)、可接受非致癌风险(ahq)、呼吸吸入参考剂量(rfd
inh
)、经口摄入参考剂量(rfd
oral
)、皮肤接触参考剂量(rfd
dermal
)、呼吸吸入致癌斜率因子(sf
inh
)、经口摄入致癌斜率因子(sf
oral
)、皮肤接触致癌斜率因子(sf
dermal
)、呼吸吸入单位致癌因子(iur)、呼吸吸入参考浓度(rfc)参数的参数值在步骤s12中通过污染物性质数据库可以获取到。本步骤针对除健康效应参数之外的暴露参数，通过调用暴露参数数据库，获取暴露参数的参数值。
[0282]
需要说明的是，暴露参数数据库中包含采集的各暴露参数的原始数据和参数值，在获取参数值时，根据当前运算时选择的某类暴露人群及所需计算区域范围进行查询，使参数值与所选择计算的人群相关。
[0283]
考虑到有些参数值虽然具有国家或地区推荐值，但是距离现在时间很久了，缺少参考价值，或者，有些参数值不具有权威性，或者并不是当地实际环境情况，这些都使计算结果存在偏差，因此基于数据库获取到参数值后，根据数据来源和数据范围识别暴露参数
的参数值是否可使用，如果不可使用，则通过蒙特卡洛计算方法得到暴露参数的参数值。具体来说，数据来源的可使用条件包括：来源于国家或地区大型调查或具有严格质量控制措施的结果，并且，来源于近10年内的研究结果，并且，来源于权威性的公开发布的数据；数据范围的可使用条件包括：数据所属区域范围在所选择的计算区域范围或相邻区域范围内。
[0284]
识别方法包括：
[0285]
如果暴露参数的参数值具有国家或地区推荐值，且同时满足数据来源和数据范围的可使用条件，则暴露参数的参数值可使用；
[0286]
如果暴露参数的参数值具有国家或地区推荐值，但不能同时满足数据来源和数据范围的可使用条件，则属于不可使用，将暴露参数加入待计算集合；
[0287]
如果暴露参数的参数值不具有国家或地区推荐值，但同时满足数据来源和数据范围的可使用条件，则属于不可使用，将暴露参数加入待计算集合；
[0288]
如果暴露参数无参数值，或者参数值不具有国家或地区推荐值，且不满足数据来源或数据范围的可使用条件，则属于异常数据。
[0289]
对异常数据可开展实际调查后获取相关计算值。
[0290]
对获取的不可使用的参数值，通过蒙特卡洛计算方法得到暴露参数的参数值，包括：
[0291]
获取待计算集合中暴露参数的原始数据，并识别原始数据是否具有数据分布特征，如果不具有数据分布特征，根据预置的分布类型，依次校验原始数据是否符合该分布类型，最终确定出数据分布特征；
[0292]
根据直接得到或确定出来的数据分布特征，进行预置次数的运算模拟，将模拟运算后的结果，作为暴露参数的最终参数值。
[0293]
示例性地，根据平均数、中位数和方差等统计数据判断数据分布类型，或者利用第三方校验函数判断数据分布类型。
[0294]
s15：根据健康效应参数和暴露参数的参数值，计算出环境介质的风险控制值，作为环境优先控制污染物的风险评估推荐值。
[0295]
具体来说，将健康效应参数的参数值、可使用的暴露参数的参数值和计算得到的待计算集合中暴露参数的参数值代入暴露计算模型和风险控制值计算模型中，计算得到致癌效应风险控制值和/或非致癌效应风险控制值，作为环境优先控制污染物在对应的环境介质的风险评估推荐值。
[0296]
优选地，考虑是否计算污染源的风险控制值，对于需要计算污染源风险控制值的优先控制污染物，引入迁移转化系数进行修正。
[0297]
具体来说，在计算前，提供“是否计算污染源控制值”的选择，如果选择否，则通过步骤s15仅计算优先控制污染物在环境介质中的人群健康风险控制值，运算结束后输出环境介质风险控制值；如果选择是，则先通过步骤s15计算优先控制污染物的环境介质人群健康风险控制值，再引入迁移转化系数的修正，根据所计算的环境介质的不同，采用不同环境介质中污染物排放源的迁移转化系数进行修正，修正步骤s15的计算结果，得到污染源风险控制值，最后输出环境介质人群健康风险控制值和污染源风险控制值，至此计算结束。
[0298]
需要说明的是，在步骤s11中对选择的每种环境介质提取到暴露特征，根据其中的排放来源，获取对应的迁移转化修正系数。
[0299]
对于污染源的每种排放来源的迁移转化修正系数根据已有资料或污染物排放后的迁移给出。
[0300]
具体来说，空气中的污染物的迁移转化可以根据污染物的扩散模型，根据所计算区域内的污染物浓度含量及污染物扩散至该区域后的落点浓度计算得到该污染物排放源对该区域内污染物的贡献来获取污染物迁移转化系数。
[0301]
土壤中的污染物可能来自于水体排放或固废排放的迁移转化，其中对于水体排放污染物迁移至土壤可根据水排放点的水处理系数、排放后水体中污染物的下渗扩散系数等参数根据水体污染物在土壤中的淋溶渗透模型等，反推计算得到水体排放点对所计算区域土壤中污染物浓度的贡献比例，作为水体排放的迁移转化系数；对于固体废物排放迁移至土壤中污染物的迁移转化系数，与水体排放类似，结合风化、雨水淋溶下渗对固废排放后对土壤中污染物的影响，结合土壤中污染物的扩散模型，得到固废排放后对土壤中污染物的迁移转化系数。
[0302]
水体中污染物的来源于水体排放和固废排放的迁移转化系数的获取与土壤类似，对于经水体排放，利用排放点的水处理系数和污染物在水体中的水质扩散模型获取得到排放点至计算区域的浓度，根据得到的计算区域内污染物的浓度和排放点浓度的比值获取迁移转化系数，同样的对于固体废物的排放，涉及到雨水淋溶下渗对水体的影响，根据淋溶下渗模型计算获取固废排放对于水体污染物的迁移转化系数。
[0303]
此外，对于不同类型的污染物而言，在获取迁移转化系数的过程中可能涉及多个迁移方式，如挥发性污染物还包括污染物的挥发过程。
[0304]
将环境介质的风险控制值与对应的每个排放来源的迁移转化修正系数相乘，得到环境介质的污染源风险控制值。
[0305]
具体来说，致癌健康效应的污染源风险控制值计算公式是：
[0306]
ccsls
p_j_f
＝致癌效应风险控制值
p_j
×
cfs
p_f
ꢀꢀꢀꢀꢀꢀ
式3-23
[0307]
其中，ccsls
p_j_f
为第p类环境介质下第j个途径的第f个来源的人体健康致癌效应的污染源风险控制值；cfs
p_f
为第p类环境介质的第f个来源的迁移转化修正系数；致癌效应风险控制值
p_j
为第p类环境介质下第j个途径的致癌效应风险控制值。
[0308]
非致癌健康效应的污染源风险控制值计算公式是：
[0309]
hcsls
p_j_f
＝非致癌效应风险控制值
p_j
×
cfs
p_f
ꢀꢀꢀꢀ
式3-24
[0310]
其中，hcsls
p_j_f
为第p类环境介质下第j个途径的第f个来源的人体健康非致癌效应的污染源风险控制值；cfs
p_f
为第p类环境介质的第f个来源的迁移转化修正系数；非致癌效应风险控制值
p_j
为第p类环境介质下第j个途径的非致癌效应风险控制值。
[0311]
采用污染物迁移转化系数修正的方法来获取污染源的风险控制值，同时获取环境介质和污染源中环境优先控制污染物的风险控制值，不仅可对环境中优先控制污染物对人体的健康影响获取直观快速的判断，同时也可获取周边区域污染物排放源对所关注区域内人群健康是否产生影响，以便及时的采取管理措施，降低环境及污染来源的排放水平，降低人群健康风险。
[0312]
与现有技术相比，本实施例提供的一种环境优先控制污染物风险评估推荐值的确定方法，基于对人群健康的保护，针对环境优先控制污染物提出了各环境介质中的风险评估推荐值的具体确定方法，在方法框架中考虑到计算参数的使用，嵌入了蒙特卡洛计算方
法，获取计算所需参数值的最优值，提高了风险评估推荐值的准确度与可靠性；在计算模块中考虑到计算优先污染物来源的控制值的情况，提出了采用污染物迁移转化系数修正的方法来获取污染源的风险控制值，可同时获取环境介质和污染源中环境优先控制污染物的风险控制值；同时给出了几种不同环境介质的风险控制值的计算方法，对于人群环境优先控制污染物的健康风险控制值的计算提供了便捷，可在短时间内获取不同环境介质中的控制值，对于人群环境健康管理非常有用，在提供参考和支撑的同时，为环境优先控制污染物的控制提供借鉴，且可为保障人群健康给出科学指导。
[0313]
本领域技术人员可以理解，实现上述实施例方法的全部或部分流程，可以通过计算机程序来指令相关的硬件来完成，所述的程序可存储于计算机可读存储介质中。其中，所述计算机可读存储介质为磁盘、光盘、只读存储记忆体或随机存储记忆体等。
[0314]
以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。

技术特征：
1.一种环境优先控制污染物风险评估推荐值的确定方法，其特征在于，包括如下步骤：根据选择的待确定风险评估推荐值的环境优先控制污染物，以及对应筛选的环境介质，提取暴露特征；根据所述暴露特征和选择的健康效应，获取健康效应参数及参数值；根据所述环境介质和暴露特征中的暴露途径，获取所述环境介质在每个暴露途径下的暴露计算模型，并根据暴露计算模型和所述选择的健康效应，获取风险控制值计算模型；根据所述暴露计算模型和风险控制值计算模型，提取出除所述健康效应参数之外的暴露参数并获取所述暴露参数的参数值；根据健康效应参数和暴露参数的参数值，计算出所述环境介质的风险控制值，作为所述环境优先控制污染物的风险评估推荐值。2.根据权利要求1所述的环境优先控制污染物风险评估推荐值的确定方法，其特征在于，所述方法还包括：根据所述暴露特征中的排放来源，获取所述环境介质的每个排放来源的迁移转化修正系数；将所述环境介质的风险控制值与对应的每个排放来源的迁移转化修正系数相乘，得到所述环境介质的污染源风险控制值，与环境介质的风险控制值一起作为所述环境优先控制污染物的风险评估推荐值。3.根据权利要求1或2所述的环境优先控制污染物风险评估推荐值的确定方法，其特征在于，所述环境介质包括空气环境、土壤环境和水体环境；其中，空气环境包括环境空气和室内空气，土壤环境包括建设用地和农用地，水体环境包括地表水和地下水；所述暴露特征的提取范围包括：暴露途径、暴露人群、排放来源、用地类型和迁移途径；其中，暴露途径包括经消化道、经呼吸道和经皮肤接触，暴露人群包括一般人群、儿童和特殊人群，排放来源包括经空气排放、经水体排放和经固废排放，用地类型包括建设用地和农用地，迁移途径包括废气、废水和废渣。4.根据权利要求3所述的环境优先控制污染物风险评估推荐值的确定方法，其特征在于，所述健康效应包括致癌效应和非致癌效应；所述根据所述暴露特征和选择的健康效应，获取健康效应参数及参数值，包括：如果选择的健康效应是致癌效应，则获取可接受致癌风险，以及，根据暴露特征中的暴露途径，获取在所述环境介质中各暴露途径致癌因子，将所述可接受致癌风险和各暴露途径致癌因子作为健康效应参数，将所述健康效应参数及参数值，用于致癌效应风险控制值计算模型；如果选择的健康效应是非致癌效应，则获取可接受非致癌风险，以及，根据暴露特征中的暴露途径，获取在所述环境介质中各暴露途径参考剂量，作为健康效应参数，将所述健康效应参数及参数值，用于非致癌效应风险控制值计算模型；如果选择的健康效应是致癌效应和非致癌效应，则同时获取上述两种情况的健康效应参数及参数值，分别用于致癌效应风险控制值和非致癌效应风险控制值计算模型。5.根据权利要求4所述的环境优先控制污染物风险评估推荐值的确定方法，其特征在于，根据所述环境介质和暴露特征中的暴露途径，获取所述环境介质在每个暴露途径下的暴露计算模型，包括：
如果环境介质是空气环境，暴露途径是经呼吸道，得到环境空气中污染物经呼吸道暴露量和室内空气中污染物经呼吸道暴露量的暴露计算模型；如果环境介质是水体环境，暴露途径是经消化道和经皮肤接触，得到饮用地表水时污染物经消化道暴露量、摄入地表水中水产品时污染物经消化道暴露量、使用地表水时污染物经皮肤接触暴露量、饮用地下水时污染物经消化道暴露量、使用地下水时污染物经皮肤接触暴露量的暴露计算模型；如果环境介质是土壤环境，暴露途径是经消化道、经皮肤接触和经呼吸道，得到人群暴露于建设用地时污染物经消化道暴露量、人群暴露于农用地时污染物经消化道暴露量、建设用地/农用地人群非气态污染物经呼吸道暴露量、建设用地/农用地人群表层土壤气态污染物室外经呼吸道暴露量、建设用地/农用地人群下层土壤气态污染物室外经呼吸道暴露量、建设用地人群地下层土壤气态污染物室内经呼吸道暴露量和人群建设用地/农用地污染物经皮肤接触暴露量的暴露计算模型。6.根据权利要求5所述的环境优先控制污染物风险评估推荐值的确定方法，其特征在于，所述根据暴露计算模型和所述选择的健康效应，获取风险控制值计算模型包括：基于空气环境中的多种暴露情景，根据对应的多个暴露计算模型，得到选择的健康效应的风险控制值计算模型；基于水体环境中每种水体类型的多种暴露途径，根据对应的多个暴露计算模型，以及各暴露途径下暴露量的叠加，得到选择的健康效应的综合风险控制值计算模型；基于土壤环境中每种用地类型的多种暴露途径，根据对应的多个暴露计算模型，以及各暴露途径下暴露量的叠加，得到选择的健康效应的综合风险控制值计算模型。7.根据权利要求6所述的环境优先控制污染物风险评估推荐值的确定方法，其特征在于，所述风险控制值计算模型中包括人群在每种暴露途径和每种暴露途径的每种方式下的暴露贡献比。8.根据权利要求1或2或7所述的环境优先控制污染物风险评估推荐值的确定方法，其特征在于，所述获取暴露参数的参数值，包括：根据数据来源和数据范围识别所述暴露参数的参数值是否可使用，如果不可使用，则通过蒙特卡洛计算方法得到暴露参数的参数值；所述数据来源的可使用条件包括：来源于国家或地区大型调查或具有严格质量控制措施的结果，并且，来源于近10年内的研究结果，并且，来源于权威性的公开发布的数据；所述数据范围的可使用条件包括：数据所属区域范围在所选择的计算区域范围或相邻区域范围内。9.根据权利要求8所述的环境优先控制污染物风险评估推荐值的确定方法，其特征在于，所述识别所述暴露参数的参数值是否可使用，包括：如果所述暴露参数的参数值具有国家或地区推荐值，且同时满足数据来源和数据范围的可使用条件，则所述暴露参数的参数值可使用；如果所述暴露参数的参数值具有国家或地区推荐值，但不能同时满足数据来源和数据范围的可使用条件，则属于不可使用，将所述暴露参数加入待计算集合；如果所述暴露参数的参数值不具有国家或地区推荐值，但同时满足数据来源和数据范围的可使用条件，则属于不可使用，将所述暴露参数加入待计算集合；如果所述暴露参数无参数值，或者参数值不具有国家或地区推荐值，且不满足数据来
源或数据范围的可使用条件，则属于异常数据。10.根据权利要求9所述的环境优先控制污染物风险评估推荐值的确定方法，其特征在于，所述通过蒙特卡洛计算方法得到暴露参数的参数值，包括：获取所述待计算集合中暴露参数的原始数据，并识别所述原始数据是否具有数据分布特征，如果不具有数据分布特征，根据预置的分布类型，依次校验原始数据是否符合该分布类型，最终确定出数据分布特征；根据直接得到或确定出来的数据分布特征，进行预置次数的运算模拟，将模拟运算后的结果，作为所述暴露参数的最终参数值。

技术总结
本发明涉及一种环境优先控制污染物风险评估推荐值的确定方法，属于环境与健康管理技术领域，解决了现有环境介质系统框架不全面、未考虑污染物来源和参数选取存在偏差的问题。包括根据选择的环境优先控制污染物和环境介质提取暴露特征；获取健康效应参数及参数值；获取环境介质在每个暴露途径下的暴露计算模型，根据暴露计算模型和选择的健康效应获取风险控制值计算模型；根据暴露计算模型和风险控制值计算模型，提取出除健康效应参数之外的暴露参数并获取暴露参数值；根据健康效应参数和暴露参数的参数值，计算出环境介质的风险控制值，作为环境优先控制污染物的风险评估推荐值。实现了准确获取多种环境介质中环境优先控制污染物风险评估推荐值。制污染物风险评估推荐值。制污染物风险评估推荐值。


技术研发人员：赵秀阁 王丹璐 李政蕾 王海燕 徐建 李琴 张婷
受保护的技术使用者：中国环境科学研究院
技术研发日：2022.03.22
技术公布日：2023/10/6