本发明涉及污染场地风险评估、信息处理、计算机程序及环境治理等技术领域，具体涉及一种用于污染场地逐年化风险评估系统的架构方法。
背景技术：
污染场地风险评估方法可以被定义为系统的程序，包括几个步骤：危害识别、暴露评估、毒性评估、风险表征及风险控制值的计算。由于风险评估涉及许多变量，为避免计算过程中的操作性失误，有助于减少错误和加快计算速度的计算机软件应运而生。虽然许多软件都被广泛用于风险评估研究，但它们都有一些特定的缺点，例如有些软件只考虑有机污染物，而另一些则需要了解污染物随时间所经历的变化。然而，所有软件的一个共同缺点是：风险评估考虑的是固定的暴露时间(即评估的是人体一生的健康风险水平)，而不是可变的时间间隔(可用于评估实时的人体健康风险水平)。没有考虑污染物浓度随时间变化及生物有效性等因素会影响污染物实际对人体产生的健康风险，最终导致评估结果过于保守，修复成本昂贵、经济效益低。
技术实现要素：
本发明的目的在于提供一种用于污染场地逐年化风险评估系统的架构方法，对污染场地进行临时分析，考虑在同一地点不同时期测量的灵活暴露时间和污染物浓度值。
本发明解决上述技术问题所采用的方案是：一种用于污染场地逐年化风险评估系统的架构方法，包括以下步骤：
(a)、建立污染场地风险识别领域的污染场地数据库，包括污染场地的土工参数、暴露情景、各采样点不同时期测量的灵活暴露时间和污染物浓度值；
(b)、建立污染场地风险识别领域的风险评估筛选值标准库，包括各地区风险评估筛选值的标准、《土壤环境质量建设用地土壤污染风险管控标准(试行)》(gb36600—2018)和《地下水质量标准》(gb/t14848—2017)，并收集各地区的土壤环境背景值；
(c)、建立污染场地风险识别领域的风险评估模型库，包括rbca模型、clea模型、csoil模型及《污染场地风险评估技术导则》(hj25.3—2014)模型；
(d)、建立污染场地风险识别领域的公共参数库，包括污染物毒理性参数和风险评估模型参数，并收集不同国家的推荐值；
(e)、采用计算机逻辑处理工具对步骤(a)、(b)、(c)、(d)的数据库之间的关系进行处理；
(f)、基于逻辑处理工具，建立污染场地风险识别过程的计算过程模型，定义计算函数、模型参数，即根据使用者要求来选择具体的风险评估模型，按照逐年化风险评估计算引擎可执行的规则建立计算控制指令集；
(g)、建立污染场地风险识别的基于生物有效性的修正模型，对评估模型进行修正；
(h)、通过逐年化风险评估计算引擎，进行不同数据量的计算；
(i)、以图文和三维形式分别输出单一污染物经单一暴露途径的逐年累积风险水平、单一污染物经所有暴露途径的逐年累积风险水平及场地所有污染物经所有暴露途径的逐年累积风险水平结果；
(j)、输出各污染物风险控制值；
其中，上述(a)、(b)、(c)、(d)不分先后。
优选地，所述步骤(a)中，将暴露情景、各采样点不同时期测量的灵活暴露时间和污染物浓度值分别设置为两个主要的输入文件：暴露情景文件和污染物浓度文件；其中暴露情景文件包括所有与暴露情景有关的信息，包括致癌和非致癌污染物的数量、采样点数量、进行风险评估的灵活暴露时间、与经口摄入和皮肤接触等暴露途径有关的暴露参数和不确定性参数；在污染物浓度文件中，数据以矩阵的形式放置，列数代表采样点个数，行中放置不同时间段测量的污染物浓度值；土工参数是指对污染场地进行土工试验得到的土壤有机质含量fom、土壤中水的渗透系数i、土壤含水率pws、土壤容重ρb及土壤颗粒密度ρs等。各污染点不同时间的污染物浓度值通过现场采样、浓度分层检测得到，该数据库所涉及的数据均根据具体项目采集。
优选地，所述步骤(b)中，风险评估筛选值标准库包括污染场地风险评估筛选值地方标准：北京的《场地土壤环境风险评价筛选值》(db11/t811—2011)、浙江的《污染场地风险评估技术导则》(db33/t892—2013，附录-土壤风险评估筛选值)、广东的《土壤重金属风险评价筛选值珠江三角洲》(db44/t1415—2014)、上海的《上海市场土壤环境健康风险评估筛选值》(试行，2015)、湖南的《重金属污染场地土壤修复标准》(db43/t1125—2016)、重庆的《场地土壤环境风险评估筛选值》(db50/t725—2016)、《土壤环境质量建设用地土壤污染风险管控标准(试行)》(gb36600—2018)、《地下水质量标准》(gb/t14848—2017)以及各地区土壤环境背景值；将各标准及土壤环境背景值以excel表格形式存储于数据库中，以便匹配使用来筛选污染物。
优选地，所述步骤(c)中，rbca模型、clea模型、csoil模型及《污染场地风险评估技术导则》(hj25.3—2014)模型中涉及的暴露时间ed为灵活的暴露时间δt，以根据使用者测量污染场地污染物浓度的时间间隔为准，当每年测一次场地污染物的浓度时，则δt＝1a，其中a表示单位年。
在污染场地风险评估的计算中，特别是场地污染物浓度随时间发生变化的情况下，将暴露时间分为多段可变的时间间隔，分别测量不同时间对应的污染物浓度再进行评估。可根据结果直观判断出场地现阶段的污染情况，对于当前污染较轻的场地，采用基于风险的场地管理办法，如通过工程或自然的手段降低污染物的浓度或切断主要暴露途径(如阻隔技术)，把风险降低到可以保护人体健康和生态安全的限度之内，不一定非要进行大规模的修复。
以hj25.3—2014模型中非敏感用地方式的情况为例，当δta＝1a时：
在步骤(c)中，暴露评估模型为：
对于单一污染物的致癌效应，经口摄入土壤途径对应的土壤暴露量为：
式中：oiserca(tj)—j年内经口摄入的土壤暴露量(致癌效应)，kg土壤·kg-1体重·d-1；osira—成人每日摄入土壤量，mg·d-1；efa—成人暴露频率，d·a-1；abso—经口摄入吸收效率因子，无量纲；δta—灵活的成人暴露期(此例中数值为1，表示每年测一次场地污染物的浓度)，a；bwa—成人体重，kg；atca—致癌效应平均时间，d。
对于单一污染物的非致癌效应，经口摄入土壤途径对应的土壤暴露量为：
式中：oisernc(tj)—j年内经口摄入的土壤暴露量(非致癌效应)，kg土壤·kg-1体重·d-1；atnc—非致癌效应平均时间，d；其余参数含义同上。
对于单一污染物的致癌效应，皮肤接触土壤途径对应的土壤暴露量为：
式中：dcserca(tj)—j年内皮肤接触途径对应的土壤暴露量(致癌效应)，kg土壤·kg-1体重·d-1；saea—成人暴露皮肤表面积，cm2；ssara—成人皮肤表面土壤粘附系数，mg·cm-2；ev—每日皮肤接触事件频率，次·d-1；其余参数含义同上。
对于单一污染物的非致癌效应，皮肤接触土壤途径对应的土壤暴露量为：
式中：dcsernc(tj)—j年内皮肤接触途径对应的土壤暴露量(非致癌效应)，kg土壤·kg-1体重·d-1；其余参数含义同上。
对于单一污染物的致癌效应，吸入土壤颗粒物途径对应的土壤暴露量为：
式中：piserca(tj)—j年内吸入土壤颗粒物对应的土壤暴露量(致癌效应)，kg土壤·kg-1体重·d-1；pm10—空气中可吸入浮颗粒物含量，mg·m-3；daira—成人每日空气呼吸量，m3·d-1；piaf—吸入土壤颗粒物在体内滞留比例，无量纲；fspi—室内空气中来自土壤的颗粒物所占比例，无量纲；fspo—室外空气中来自土壤的颗粒物所占比例，无量纲；efoa—成人的室外暴露频率，d·a-1；efia—成人的室内暴露频率，d·a-1；其余参数含义同上。
对于单一污染物的非致癌效应，吸入土壤颗粒物途径对应的土壤暴露量为：
式中：pisernc(tj)—j年内吸入土壤颗粒物对应的土壤暴露量(非致癌效应)，kg土壤·kg-1体重·d-1；其余参数含义同上。
对于单一污染物的致癌效应，吸入室外空气中来自表层土壤的气态污染物对应的土壤暴露量为：
式中：ioverca1(tj)—j年内吸入室外空气中来自表层土壤的气态污染物对应的土壤暴露量(致癌效应)，kg土壤·kg-1体重·d-1；vfsuroa—表层土壤中污染物扩散进入室外空气的挥发因子，kg·m-3；其余参数含义同上。
对于单一污染物的非致癌效应，吸入室外空气中来自表层土壤的气态污染物对应的土壤暴露量为：
式中：iovernc1(tj)—j年内吸入室外空气中来自表层土壤的气态污染物对应的土壤暴露量(非致癌效应)，kg土壤·kg-1体重·d-1；其余参数含义同上。
对于单一污染物的致癌效应，吸入室外空气中来自下层土壤的气态污染物对应的土壤暴露量为：
式中：ioverca2(tj)—j年内吸入室外空气中来自下层土壤的气态污染物对应的土壤暴露量(致癌效应)，kg土壤·kg-1体重·d-1；vfsuboa—下层土壤中污染物扩散进入室外空气的挥发因子，kg·m-3；其余参数含义同上。
对于单一污染物的非致癌效应，吸入室外空气中来自下层土壤的气态污染物对应的土壤暴露量为：
式中：iovernc2(tj)—j年内吸入室外空气中来自下层土壤的气态污染物对应的土壤暴露量(非致癌效应)，kg土壤·kg-1体重·d-1；其余参数含义同上。
对于单一污染物的致癌效应，吸入室外空气中来自地下水的气态污染物对应的地下水暴露量为：
式中：ioverca3(tj)—j年内吸入室外空气中来自地下水的气态污染物对应的地下水暴露量(致癌效应)，l地下水·kg-1体重·d-1；vfgwoa—地下水中污染物扩散进入室外空气的挥发因子，l·m-3；其余参数含义同上。
对于单一污染物的非致癌效应，吸入室外空气中来自地下水的气态污染物对应的地下水暴露量为：
式中：iovernc3(tj)—j年内吸入室外空气中来自地下水的气态污染物对应的地下水暴露量(非致癌效应)，l地下水·kg-1体重·d-1；其余参数含义同上。
对于单一污染物的致癌效应，吸入室内空气中来自下层土壤的气态污染物对应的土壤暴露量为：
式中：iiverca1(tj)—j年内吸入室内空气中来自下层土壤的气态污染物对应的土壤暴露量(致癌效应)，kg土壤·kg-1体重·d-1；vfsubia—下层土壤中污染物扩散进入室内空气的挥发因子，kg·m-3；其余参数含义同上。
对于单一污染物的非致癌效应，吸入室内空气中来自下层土壤的气态污染物对应的土壤暴露量为：
式中：iivernc1(tj)—j年内吸入室内空气中来自下层土壤的气态污染物对应的土壤暴露量(非致癌效应)，kg土壤·kg-1体重·d-1；其余参数含义同上。
对于单一污染物的致癌效应，吸入室内空气中来自地下水的气态污染物对应的地下水暴露量为：
式中：iiverca2(tj)—j年内吸入室内空气中来自地下水的气态污染物对应的地下水暴露量(致癌效应)，l地下水·kg-1体重·d-1；vfgwia—地下水中污染物扩散进入室外空气的挥发因子，l·m-3；其余参数含义同上。
对于单一污染物的非致癌效应，吸入室内空气中来自地下水的气态污染物对应的地下水暴露量为：
式中：iivernc2(tj)—j年内吸入室内空气中来自地下水的气态污染物对应的地下水暴露量(非致癌效应)，l地下水·kg-1体重·d-1；其余参数含义同上。
对于单一污染物的致癌效应，饮用地下水途径对应的地下水暴露量为：
式中：cgwerca(tj)—j年内饮用受影响地下水对应的地下水暴露量(致癌效应)，l地下水·kg-1体重·d-1；gwcra—成人每日饮水量，l地下水·d-1；其余参数含义同上。
对于单一污染物的非致癌效应，饮用地下水途径对应的地下水暴露量为：
式中：cgwernc(tj)—j年内饮用受影响地下水对应的地下水暴露量(非致癌效应)，l地下水·kg-1体重·d-1；其余参数含义同上。
在步骤(c)中，计算致癌风险的模型则为：
经口摄入土壤途径的致癌风险为：
式中：crois—j年内经口摄入土壤途径的累积致癌风险，无量纲；csur(δta)—每段灵活的成人暴露期(此例中指1a)所对应的场地表层土壤中污染物浓度，mg·kg-1；sfo—经口摄入致癌斜率因子，mg污染物·kg-1体重·d-1；其余参数含义同上。
皮肤接触土壤途径的致癌风险为：
式中：crdcs—j年内皮肤接触土壤途径的累积致癌风险，无量纲；sfd—皮肤接触致癌斜率因子，mg污染物·kg-1体重·d-1；其余参数含义同上。
吸入土壤颗粒物途径的致癌风险为：
式中：crpis-j年内吸入土壤颗粒物途径的累积致癌风险，无量纲；sfi-呼吸吸入致癌斜率因子，mg污染物·kg-1体重·d-1；其余参数含义同上。
吸入室外空气中来自表层土壤的气态污染物途径的致癌风险为：
式中：criov1-j年内吸入室外空气中来自表层土壤的气态污染物途径的累积致癌风险，无量纲；其余参数含义同上。
吸入室外空气中来自下层土壤的气态污染物途径的致癌风险为：
式中：criov2-j年内吸入室外空气中来自下层土壤的气态污染物途径的累积致癌风险，无量纲；csub(δta)-每段灵活的成人暴露期(此例中指1a)所对应的场地下层土壤中污染物浓度，mg·kg-1；其余参数含义同上。
吸入室内空气中来自下层土壤的气态污染物途径的致癌风险为：
式中：criiv1-j年内吸入室内空气中来自下层土壤的气态污染物途径的累积致癌风险，无量纲；其余参数含义同上。
土壤中单一污染物经所有暴露途径的总致癌风险为：
cragg(tj)＝crois(tj)+crdcs(tj)+crpis(tj)+criov1(tj)+criov2(tj)+criiv1(tj)
式中：cragg(tj)-j年内土壤中单一污染物经所有暴露途径的累积总致癌风险，无量纲；其余参数含义同上。
土壤中所有污染物经所有暴露途径的总致癌风险为：
式中：crcuml(tj)-j年内土壤中所有污染物经所有暴露途径的累积总致癌风险，无量纲。
吸入室外空气中来自地下水的气态污染物途径的致癌风险为：
式中：criov3—j年内吸入室外空气中来自地下水的气态污染物途径的累积致癌风险，无量纲；cgw(δta)—每段灵活的成人暴露期(此例中指1a)所对应的场地地下水中污染物浓度，mg·kg-1；其余参数含义同上。
吸入室内空气中来自地下水的气态污染物途径的致癌风险为：
式中：criiv2—j年内吸入室内空气中来自地下水的气态污染物途径的累积致癌风险，无量纲；其余参数含义同上。
饮用地下水途径的致癌风险为：
式中：crcgw—j年内饮用地下水途径的累积致癌风险，无量纲。
地下水中单一污染物经所有暴露途径的总致癌风险为：
cragg(tj)＝criov3(tj)+criiv2(tj)+crcgw(tj)
式中：cragg(tj)—j年内地下水中单一污染物经所有暴露途径的累积总致癌风险，无量纲；其余参数含义同上。
地下水中所有污染物经所有暴露途径的总致癌风险为：
式中：crcuml(tj)—j年内地下水中所有污染物经所有暴露途径的累积总致癌风险，无量纲。
在步骤(c)中，计算危害商的模型则为：
经口摄入土壤途径的危害商为：
式中：hqois(tj)—j年内经口摄入土壤途径的累积危害商，无量纲；rfdo—经口摄入参考剂量，mg污染物·kg-1体重·d-1；saf—暴露于土壤的参考剂量分配系数，无量纲；其余参数含义同上。
皮肤接触土壤途径的危害商为：
式中：hqdcs(tj)-j年内皮肤接触土壤途径的累积危害商，无量纲；rfdd-皮肤接触参考剂量，mg污染物·kg-1体重·d-1；其余参数含义同上。
吸入土壤颗粒物途径的危害商为：
式中：hqpis(tj)-j年内吸入土壤颗粒物途径的累积危害商，无量纲；rfdi-呼吸吸入参考剂量，mg污染物·kg-1体重·d-1；其余参数含义同上。
吸入室外空气中来自表层土壤的气态污染物途径的危害商为：
式中：hqiov1(tj)-j年内吸入室外空气中来自表层土壤的气态污染物途径的累积危害商，无量纲；其余参数含义同上。
吸入室外空气中来自下层土壤的气态污染物途径的危害商为：
式中：hqiov2(tj)-j年内吸入室外空气中来自下层土壤的气态污染物途径的累积危害商，无量纲；其余参数含义同上。
吸入室内空气中来自下层土壤的气态污染物途径的危害商为：
式中：hqiiv1(tj)-j年内吸入室内空气中来自下层土壤的气态污染物途径的累积危害商，无量纲；其余参数含义同上。
土壤中单一污染物经所有暴露途径的危害指数为：
hiagg(tj)＝hqois(tj)+hqdcs(tj)+hqpis(tj)+hqiov1(tj)+hqiov2(tj)+hqiiv1(tj)
式中：hiagg(tj)-j年内土壤中单一污染物经所有暴露途径的累积危害指数，无量纲；其余参数含义同上。
土壤中所有污染物经所有暴露途径的危害指数为：
式中：hicuml(tj)—j年内土壤中所有污染物经所有暴露途径的累积危害指数，无量纲。
吸入室外空气中来自地下水的气态污染物途径的危害商为：
式中：hqiov3(tj)—j年内吸入室外空气中来自地下水的气态污染物途径的累积危害商，无量纲；其余参数含义同上。
吸入室内空气中来自地下水的气态污染物途径的危害商为：
式中：hqiiv2(tj)—j年内吸入室内空气中来自地下水的气态污染物途径的累积危害商，无量纲；其余参数含义同上。
饮用地下水途径的危害商为：
式中：hqcgw(tj)—j年内饮用地下水途径的累积危害商，无量纲；waf—暴露于地下水的参考剂量分配系数，无量纲；其余参数含义同上。
地下水中单一污染物经所有暴露途径的危害指数为：
hiagg(tj)二hqiov3(tj)+hqiiv2(tj)+hqcgw(tj)
式中：hiagg(tj)—j年内地下水中单一污染物经所有暴露途径的累积危害指数，无量纲；其余参数含义同上。
地下水中所有污染物经所有暴露途径的危害指数为：
式中：hicuml(tj)—j年内地下水中所有污染物经所有暴露途径的累积危害指数，无量纲。
在步骤(c)中，计算土壤和地下水风险控制值的模型则为：
基于上述6种土壤暴露途径综合致癌效应的土壤风险控制值为：
式中：rcvs—单一污染物基于6种土壤暴露途径综合致癌效应的土壤风险控制值，mg·kg-1；acr—可接受致癌风险，无量纲；其余参数含义同上。
基于上述6种土壤暴露途径综合非致癌效应的土壤风险控制值为：
式中：hcvs—单一污染物基于6种土壤暴露途径综合非致癌效应的土壤风险控制值，mg·kg-1；ahq—可接受危害商，无量纲；其余参数含义同上。
基于上述3种地下水暴露途径综合致癌效应的地下水风险控制值为：
式中：rcvg—单一污染物基于3种地下水暴露途径综合致癌效应的地下水风险控制值，mg·kg-1；其余参数含义同上。
基于上述3种地下水暴露途径综合非致癌效应的地下水风险控制值为：
式中：hcvg—单一污染物基于3种地下水暴露途径综合非致癌效应的地下水风险控制值，mg·kg-1；其余参数含义同上。
优选地，所述步骤(d)中，污染物毒理性参数包括：经口摄入致癌斜率因子sfo、经口摄入参考剂量rfdo、呼吸吸入参考浓度rfc、消化道吸收因子absgi、皮肤接触吸收效率因子absd、无量纲亨利常数h′、空气中扩散系数da、水中扩散系数dw、水溶解度s等；风险评估模型参数包括：暴露频率ef、体重bw、身高h、每日空气呼吸量dair、每日饮用水量gwcr、每日摄入土壤量osir、每日皮肤接触事件频率ev、经口摄入吸收因子abso、单一污染物可接受致癌风险acr、可接受危害商ahq、平均时间at等。
优选地，所述步骤(h)中，所述的计算引擎在接受计算控制指令集后，引擎计算单元，根据计算资源、计算方式和依赖关系建立计算调度模型，然后进行计算处理区调度计算。
优选地，所述步骤(g)中，基于生物有效性对评估模型进行修正是指根据生物有效性研究方法得到的修正系数用于风险评估模型的修正；
修正系数：
式中：cbio为进入人体循环系统的污染物浓度，mg·kg-1；csoil为土壤中的污染物浓度，mg·kg-1。
则修正后的土壤中单一污染物经所有暴露途径的总致癌风险为：
cragg(tj)＝[crois(tj)+crdcs(tj)+crpis(tj)+criov1(tj)+criov2(tj)+criiv1(tj)]×cf
对其它用于风险表征的计算模型的修正依此类推。
本发明的有益效果为：本发明提出了一种逐年化计算污染场地风险评估系统的架构方法,对污染场地进行临时分析，考虑在同一地点不同时期测量的灵活暴露时间和污染物浓度值。在污染物浓度随时间变化的地区发生事故或进行风险评估时，这一特性非常有用。这些信息将使使用者不仅知道人体在什么时期会处于有效的危险之中，而且能够在污染尚未达到临界水平的时期采取更一致的风险管控措施。
与现有技术相比，本发明的方法基于计算机技术实现对污染场地的逐年化风险评估，并考虑到生物有效性对实际风险水平产生的影响，通过图文和三维形式展现污染场地在灵活的暴露时间和污染物浓度随时间发生变化的情况下的风险水平结果。本发明中的系统架构方法允许更改暴露时间，这对于在更短的时间间隔内评估风险非常有用，且该方法进行的时空分析使得对人体健康风险的评估更加现实和详细。
附图说明
图1为根据本发明方法所形成的污染场地逐年化风险评估系统架构；
图2为本发明逐年化风险评估计算流程图；
图3为本发明污染物浓度输入文件示例。
具体实施方式
为更好的理解本发明，下面的实施例是对本发明的进一步说明，但本发明的内容不仅仅局限于下面的实施例。
如图1-3所示，本发明污染场地逐年化风险评估的实现方法主要包括建立污染场地数据库、风险评估筛选值标准库、风险评估模型库、公共参数库；采用计算机逻辑处理工具对上述数据库之间的关系进行处理；基于逻辑处理工具，建立污染场地风险识别过程的计算过程模型；建立基于生物有效性的修正模型；采用计算引擎进行逐年化风险评估计算等。下面对各个数据库及模型的建立进行说明：
所述的污染场地数据库包括污染场地土工参数、暴露情景及各污染点不同时间的污染物浓度值。由于污染场地数据库涉及的所有数据均与具体污染场地相关，为避免数据混淆，将暴露情景、各采样点不同时期测量的灵活暴露时间和污染物浓度值分别设置为两个主要的输入文件：暴露情景文件和污染物浓度文件。其中暴露情景文件包括所有与暴露情景有关的信息，如变量：致癌和非致癌污染物的数量、采样点数量、进行风险评估的灵活暴露时间、与经口摄入和皮肤接触等暴露途径有关的暴露参数和不确定性参数；在污染物浓度文件中，数据以矩阵的形式放置，列数代表采样点个数，行中放置不同时间段测量的污染物浓度值。各污染点不同时间的污染物浓度值通过现场采样、浓度分层检测得到，该数据库所涉及的数据均根据具体项目采集。
所述的风险评估筛选值标准库包括污染场地风险评估筛选值地方标准：北京的《场地土壤环境风险评价筛选值》(db11/t811—2011)、浙江的《污染场地风险评估技术导则》(db33/t892—2013，附录—土壤风险评估筛选值)、广东的《土壤重金属风险评价筛选值珠江三角洲》(db44/t1415—2014)、上海的《上海市场土壤环境健康风险评估筛选值》(试行，2015)、湖南的《重金属污染场地土壤修复标准》(db43/t1125—2016)、重庆的《场地土壤环境风险评估筛选值》(db50/t725—2016)等；《土壤环境质量建设用地土壤污染风险管控标准(试行)》(gb36600—2018)；《地下水质量标准》(gb/t14848—2017)以及各地区土壤环境背景值。
所述的风险评估模型库包括美国的rbca模型、英国的clea模型、荷兰的csoil模型及我国的《污染场地风险评估技术导则》(hj25.3—2014)模型
所属的公共参数库包括污染物毒理性参数和风险评估模型参数。其中污染物毒理性参数有经口摄入致癌斜率因子sfo、经口摄入参考剂量rfdo、呼吸吸入参考浓度rfc、消化道吸收因子absgi、皮肤接触吸收效率因子absd、无量纲亨利常数h′、空气中扩散系数da、水中扩散系数dw、水溶解度s等；风险评估模型参数有暴露频率ef、体重bw、身高h、每日空气呼吸量dair、每日饮用水量gwcr、每日摄入土壤量osir、每日皮肤接触事件频率ev、经口摄入吸收因子abso、单一污染物可接受致癌风险acr、可接受危害商ahq、平均时间at等。
采用计算机逻辑处理工具对上述数据库之间的关系进行处理；
基于逻辑处理工具，建立污染场地风险识别过程的计算过程模型，定义计算函数、模型参数，即根据使用者要求来选择具体的风险评估模型，按照逐年化风险评估计算引擎可执行的规则建立计算控制指令集；
所述的基于生物有效性的修正模型通过建模工具，将根据生物有效性研究方法得到的修正系数用于风险评估模型的修正。
所述的计算引擎在接受计算控制指令集后，引擎计算单元，根据计算资源、计算方式和依赖关系建立计算调度模型，然后进行计算处理区调度计算，并以图文、三维及表格形式输出计算结果。
具体计算流程如图2所示。由于执行风险评估计算需要许多变量，因此必须以可靠的方式组织所有参数，以便能够快速有效地访问系统程序的例程。因此，在污染场地数据库中使用两个主要输入文件。为了增加系统程序的灵活性和实用性，通过输入文件选择暴露途径和涉及的变量，以避免源代码中的更改。
两个主要的输入文件是指：暴露情景文件和污染物浓度文件。其中在暴露情景文件中，输入所有与暴露情景有关的信息，如变量：致癌和非致癌污染物的数量、采样点数量、进行逐年化风险评估的灵活暴露时间、与经口摄入和皮肤接触等暴露途径有关的暴露参数和不确定性参数；在污染物浓度文件中，数据以矩阵的形式放置，列数代表采样点个数，行中放置不同时间段测量的污染物浓度值。
第一步：输入暴露情景文件，子程序读取暴露情景文件以获得该污染场地的与暴露情景有关的信息并进行处理和存储，同时获取污染场地数据库中的土工参数；
第二步：输入污染物浓度文件，子程序读取污染场地中不同时间段测量的污染物浓度值，并自行执行与对应的风险评估筛选值进行比较的命令，对污染物浓度值进行处理和存储。
第三步：子程序从公共参数库中读取污染场地逐年化风险评估计算所需要的污染物毒理性参数和风险评估模型参数信息并存储。
第四步：主要程序调度处理各子程序的信息，从风险评估模型库选择匹配的计算模型并根据污染场地信息进行模型修正。
第五步：计算处理器根据计算资源、计算方式和依赖关系自动进行调度计算，将计算结果存储到结果数据中。
第六步：最后根据使用者要求，读取计算结果，以图文、三维及表格形式输出计算结果，结束程序。
上述对本发明计算机程序的实现流程进行了说明，通过计算机技术和逻辑处理工具等实现了污染场地的逐年化风险评估。风险评估结果将清晰明了地展示各段时间间隔内及累积的风险水平，以便于使用者在污染尚未达到临界水平的时期采取更一致的减轻行动。
污染物浓度文件中，每个特定采样点不同采样时间的浓度值以矩阵形式排列，列数代表采样点个数，行中放置不同时间段测量的污染物浓度值。如图3所示，该图仅作为示例，有关数据均无实际意义。
以上所述是本发明的优选实施方式而已，当然不能以此来限定本发明之权利范围，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和变动，这些改进和变动也视为本发明的保护范围。