CN102622670A - 一种突发水环境污染事故涉嫌风险源的溯源方法 - Google Patents

Abstract

针对突发水环境污染事故涉嫌风险源难以快速排查问题，本发明建立了一种突发水环境污染事故涉嫌风险源溯源方法，其步骤包括：(1)企业风险信息数据库构建；(2)污染事故应急监测及信息获取；(3)事故涉嫌风险源初步筛选；(4)突发事故涉嫌风险源定位；(5)事故涉嫌风险源二次筛选；(6)涉嫌风险源事故场景符合性分析；(7)涉嫌风险源现场取证及验证。本发明建立的标准化排查方法解决了以往水环境污染事故涉嫌风险源排查过程中存在的主观、滞后及工作量大等问题，便于环保主管部门实现水环境污染事故应急过程的规范化。

Description

一种突发水环境污染事故涉嫌风险源的溯源方法

技术领域

本发明涉及一种专用于突发水环境污染事故涉嫌风险源的溯源方法，属于环境风险影响评估技术领域。

背景技术

近年来我国逐步进入突发性水环境污染事件密集高发期，统计显示自2006年我国平均每两天发生一起环境突发事故，其中70％是水污染事故。突发水环境污染事故已成为危害人身健康、破坏生态环境的重要因素，严重威胁了我国环境、经济及社会的和谐发展。如2005年松花江硝基苯污染事故导致300多万人饮水困难；2008年云南澄江阳宗海含砷废水泄漏污染事件导致下游农业和渔业严重损失；2009年山东含砷废水污染事故严重影响了下游50余万人的健康和生产。相关污染事故均加剧了流域水质恶化，导致了严重的经济损失和生态破坏，还造成了区域社会恐慌。随着水环境污染事故的频繁发生，如何在事故应急过程中及时、客观、准确识别涉嫌风险源、明确企业责任、避免事故的进一步恶化成为事故应急工作的重要内容，也成为降低事故再次发生的关键。

然而，目前国内外缺乏规范、合理、有效的水环境污染事故涉嫌风险源排查办法。大部分环境污染事故涉嫌风险源的排查主要基于企业自行上报、群众举报，以及事故调查者个人经验，通常需要完成大量现场调查及取证工作。由于现有排查方法具有不规范性及滞后性，缺乏科学性，导致排查结果的权威性及时效性往往受到广泛质疑，这成为环境污染事故涉嫌风险源排查中遇到的最棘手问题。

环境污染事故涉嫌风险源的溯源技术完全不同于现有的环境风险源识别技术，也不同于环境事故预警技术。环境污染事故风险源识别及预警技术主要针对某一固定的风险单元，通过最大可能事故场景模拟，评估事故未来某一个(段)时间节点可能产生的危害性后果。而涉嫌风险源的溯源技术主要基于事故发生过程中某几个时间节点污染带的变化，采用数学模型反推事故发生的初始位置、事故类型及爆发时间。由于事故爆发时间、污染带位置及长度等具有不确定性及关联性，导致事故源头的准确定位难以实现。

目前已有的溯源技术多集中在通信及互联网领域(如CN100574521C；CN101894345A和CN101252592A)，以及采用基因技术对农产品或肉类制品来源进行追踪(CN102142116A；CN102063692；CN101717820A)，用于突发环境污染事故风险源溯源的技术比较少(CN102116734A；CN101789162A)。仅有个别研究提出采用同位素分析或指纹图谱的方法，对比事故化学品与风险源化学品组成差异，来实现涉嫌风险源的排查。该方法能显著提升经验判别法的准确性，但该方法仅适用于某一类特征污染物，如重金属、石油类等，且技术监测方法复杂、排查工作量大、监测工作量大、缺乏标准程序，目前难以直接用于事故涉嫌风险源的逐一排查。

由于缺乏切实可行、规范标准的水环境污染事故涉嫌风险源排查方法，当前国内外环保管理部门仍主要依靠经验法逐一检查来排查事故涉嫌风险源，导致排查结果准确性、时效性均受到质疑，如何规划、科学、合理、有效的排查事故涉嫌风险源成为国内外事故应急管理部门迫切需要解决的问题。

发明内容

本发明规范了水环境污染事故涉嫌风险源的排查方法，建立了基于污染事故复杂性逐级解析的涉嫌风险源信息获取、多级筛选、符合性分析及取证验证的标准化程序，解决了以往污染事故涉嫌风险源排查过程中方法主观、滞后、工作量大的问题，减小了事故各方对排查结果的质疑。

本发明的技术方案如下：

1.一种突发水环境污染事故涉嫌风险源的溯源方法，其特征在于：其具体操作步骤如下，步骤一、风险源风险信息数据库构建

第一步：建立基于WebGis的流域风险源分布软件***；

***包括风险源名称、经济产值、所属行业、风险源坐标、污水排放口坐标、受纳水体等风险源基础信息；

第二步：建立定时更新的风险源风险物质数据库；

通过企业定期申报、不定期现场检查等更新手段，建立流域排污风险源风险物质数据库，主要包括：有毒有害原料、辅料、产品和废弃物，并定期更新风险源风险物质名称、储量、存储方式、监控手段、存储场所等信息；

第三步：建立流域高风险化学品组成信息库；

采用环境风险源辨识法，筛选流域内风险源具有的高毒有机物、重金属和油类等高风险化学品，采用内标法测定此类风险化学品的典型标志物，建立流域高风险化学品组成信息库；

步骤二、污染事故应急监测及信息获取

第一步：获取事故污染物组成；

接到举报或发现水质异常后，开展事故应急监测，测定事故污染物主要组成；

第二步：获取污染带变化信息；

采用流动监测(车、船等)、简易监测、低空航测、遥感等监测手段，在事故发现地点上下游获取污染带扩散速率、长度、宽度，以及浓度分布；

第三步：计算污染物泄漏总量；

基于第二步获取的污染带变化信息，利用公式(1)计算事故泄漏的污染物总量；

$M={\mathrm{\Σ}}_1^n(C_i\·L_i)---\left(1\right)$

其中，M-污染物泄漏总量，t；C_i-监测点位浓度，mg/L；L_i-监测点位间隔距离，km；n-监测点位个数，通常大于10；

步骤三、事故涉嫌风险源初步筛选

依据从步骤二中获取的污染物种类及泄漏量等事故信息，利用排除法从流域风险源风险物质数据库和流域高风险化学品组成信息库中筛选事故涉嫌风险源，初步确定嫌疑风险源范围；

步骤四、突发事故涉嫌风险源定位

第一步：事故最大可能类型分析；

根据污染物的理化性质、污染带物质分布特点，判定事故的最大可能类型，分析事故是否属于废水非正常排放事故、液态化学品直接泄漏事故或液态化学次生泄漏事故等；依据案例分析结果，预测事故的初始污染带长度范围；

第二步：事故最早扩散时间分析；

依据河流基本参数及步骤二中获得的污染带变化数据，采用式(2)计算污染带长度与扩散时间的关系，并确定污染事故的最早扩散时间；

$\frac{L}{t}=\∝\·\ln \frac{M}{C_0\·t^{0.5}}-\mathrm{\β}\·u_0-k---\left(2\right)$

其中，L-不同时间节点污染带长度，m；t-污染带扩散时间，h；C_o-污染带监测的最高浓度，mg/L；u_o-污染带长度扩展速率，m/h；α、β和k-河流常数；

第三步：事故最大扩散范围分析；

依据事故最早扩散时间反推事故最大可能发生地点，按照式(3)计算事故的最大扩散距离；

$X=u_{x}t+0.63\sqrt{\frac{t}{H}\ln \left(\frac{M}{1.12\sqrt{H\·t}\·B\·C_0}\right)}---\left(3\right)$

其中，X-污染物泄漏最大长度，m；u_x-河流流速，m/s；t-污染带扩散时间，h；M-污染物泄露总量；C₀-污染带监测的最高浓度，mg/L；B-河流平均宽度，m；H-河流平均深度，m；

步骤五、事故涉嫌风险源二次筛选

基于步骤四中获得的事故最大扩散范围，依据风险源分布的WebGis***，对第一次筛选后涉嫌风险源进行第二次筛选，进一步缩小事故涉嫌风险源范围；

步骤六、涉嫌风险源事故场景符合性分析

第一步：风险源事故场景模拟；

获取步骤五中二次筛选后涉嫌风险源的环境风险源信息，根据风险源最大可能事故场景分析，模拟风险源同类事故中扩散时间、泄漏量和扩散距离的关系；

第二步：事故场景符合性分析；

对比实际事故监测结果与涉嫌风险源事故场景模拟结果，按照公式(4)计算涉嫌风险源事故场景的符合性：

$f(X,M,t)=\mathrm{\γ}\frac{X}{X_0}+\mathrm{\δ}\frac{M}{M_0}+\mathrm{\ϵ}\frac{t}{t_0}---\left(4\right)$

其中，f(X，M，t)-事故场景符合性因子；X，M，t(X₀，M₀，t₀)-分别为涉嫌风险源事故模拟场景(及实际监测)的事故扩散距离、泄漏量和扩散时间；γ、δ和ε-权重因子，且γ+δ+ε＝1；

第三步：涉嫌风险源事故场景符合性排序；

基于第二步计算出的事故场景符合因子大小，建立事故涉嫌风险源的符合性排序清单；

步骤七、涉嫌风险源现场取证及验证

依据建立的事故涉嫌风险源的符合性排序清单，依次开展涉嫌风险源的现场调查、取证及后果验证；获取涉嫌风险源中与事故化学品相同的材料，采取稳定同位素分析法，分析化学品的组成或特征指纹图库；通过样品成分比对分析，确定事故风险源，明确企业责任。

2.如权利要求1所述的一种突发水环境污染事故涉嫌风险源的溯源方法，其特征在于：步骤三中事故涉嫌污染风险源的初步筛选应先排除不存储、不使用及不生产事故化学品的风险源，再排除事故化学品存储量、使用及生产量超过事故泄漏量的风险源。

3.如权利要求1所述的一种突发水环境污染事故涉嫌风险源的溯源方法，其特征在于：当经步骤三(或步骤五)排除后涉嫌风险源数小于5时，可直接超越到步骤七开展涉嫌风险源的现场调查、取证及验证。

4.如权利要求1所述的一种突发水环境污染事故涉嫌风险源的溯源方法，其特征在于：风险源事故场景模拟由公式(5)和公式(6)依次确定：

确定最大可信事故源强，依据单位时间内液态有毒有害危险化学品进入地表水体中的含量，采用下式进行计算：

Q_源强＝λ×Q_max×P_泄漏比例/t    (5)

其中，Q_max-有毒有害危险化学品最大储量，t；P_泄漏比例-有毒有害危险化学品的泄漏比例；λ-液态泄露污染物自流入敏感地表水体的比例；t-泄漏时间，min；

对于水污染事故风险源导致的地表水体污染，运用二维扩散模型计算事故危害范围，表达式如下：

其中，ρ(x，y，t)-泄漏点下游距离x处、t时溶解态污染物浓度，mg/L；Q_源强-污染源源强，g；D_x/D_y-纵向及横向离散系数，m²/s；u_x/u_y-横向及纵向流速，m/s；h-平均水深，m；K-降解速率，g/g·d。

附图说明

图1为突发水环境污染事故涉嫌风险源的溯源方法工作流程图

具体实施方式

下面将结合附图和实施例对本发明做进一步的详细说明：

本方法涉及的一种突发水环境污染事故涉嫌风险源的溯源方法，适用于河流类水体中突发水污染事故中涉嫌风险源的逐级排查，包括：向环境排放污染物的单位，生产、贮存、经营、使用危险物质的单位或产生、收集、利用、处置危险废物等可能发生环境污染事件的单位。本方法不适用于长期累积性事故风险源的识别，也不适合车、船等移动源的事故溯源。

根据各流域涉嫌风险源的复杂程度，事故涉嫌风险源的排查也分类、分级逐次开展。对于事故涉嫌风险源少、污染物质专属性强的流域，简单的一级或二级筛选即可实现涉嫌风险源有效排查，为简化步骤可直接从步骤三或步骤五直接超越到步骤七进行风险源的验证；对于涉嫌风险源数量多、污染物比较常见的流域，需要按步骤依次排查。

本发明提出的突发水环境污染事故涉嫌风险源的溯源方法，包括以下几个步骤：

步骤一、风险源风险信息数据库构建

1.建立基于WebGis的流域风险源分布软件***。风险源的WebGid分布体系是进行涉嫌风险源溯源的基础。只有建立尽可能详尽的风险源分布数据库，才能避免涉嫌风险源排查中产生涉嫌企业漏报、错报等问题。该***中应包含所有可能用于环境污染事故风险源管理的要素，如风险风险源名称、经济产值、所属行业、风险源坐标、污水排放口坐标、受纳水体等。

2.建立定时更新的风险源风险物质数据库。通过企业定期申报、现场排查等多种方式，建立包含河流周边风险源主要有毒有害原辅料、产品、废弃物的基本信息数据库；为保障风险物质信息的有效性，应定时更新风险源环境风险污染名称、数量、所处场所等基本信息；

3.建立流域高风险化学品组成信息库。通过临界值法，筛选流域风险源内具有高毒的有机物、重金属和油类等的高风险单元，采用内标法测定风险化学品中典型标志物，建立流域高风险化学品组成库；在高风险污染事故发生时便于缩小重点排查单元。

步骤二、污染事故应急监测及信息获取

1.获取事故污染物组成。开展事故应急监测，确定事故污染物主要组成；具体分析方法可以根据流域主要突发环境污染事故特点及监测水平来确定，也可以采用GC-MS、稳定同位素法，或采用简单的化学定量方法；

2.获取污染带变化信息。在事故发现地上下游，采用流动监测(车、船等)、简易监测、低空航测、遥感等监测手段，获取污染带扩散速率、长度、宽度，以及浓度分布；

3.计算污染物泄漏总量。针对第二步获取的污染带基本信息，通过公式(1)计算事故泄漏污染物的总量；

$M={\mathrm{\Σ}}_1^n(C_i\·L_i)---\left(1\right)$

M-污染物泄漏总量，t；C_i-监测点位浓度，mg/L；L_i-监测点位间隔距离，km；n-监测点位个数，通常大于10。

在实际事故涉嫌风险源排查中，考虑到污染物的降解及吸附等物理化学作用，其在顺河流迁移过程中总量会有一定的损失，因而对计算出的污染物泄漏总量进行修正。修正方法可以根据式(2)计算：

M’＝a·b·M    (2)

a为污染物降解系数，b为污染物吸附系数

步骤三、事故涉嫌风险源初步筛选

实现涉嫌风险源依据获取的污染物类型及泄漏量等事故信息，利用排除法从步骤一中建立的流域风险源风险污染基础数据库和流域高风险化学品组成信息库中筛选流域内涉嫌污染风险源，初步确定嫌疑风险源范围。事故涉嫌污染风险源的初步筛选应先筛选存储、使用及生产该化学物质的风险源，再筛选事故化学品存储量、使用及生产量超过事故泄漏量的风险源。如排查后涉嫌风险源低于5家，可不进行涉嫌风险源的定位及二次筛选，直接超越到步骤七开展涉嫌风险源的验证工作。

步骤四、突发事故涉嫌风险源定位

1.事故最大可能类型分析。事故最大可能类型的选取将直接决定模型中各参数的选取，进而影响污染带的分布、事故符合性分析结果，并进而对涉嫌风险源排序具有一定的影响。因而事故最大可能类型的选取也是涉嫌风险源排查中重要的一环，其具体操作方法按照水环境风险源识别方法开展。根据污染带物质分布形式，以及污染物的理化性质，判定事故的最大可能类型；水污染事故风险源出现的事故类型主要有：废水非正常排放事故、液态化学品直接泄漏事故、液态化学次生泄漏事故等。

2.事故最早扩散时间分析。根据河流参数及最大可能事故类型分析，依据在步骤二中获得的污染带变化数据基础上，采用式(3)确定扩散带长度与扩散时间的关系。并基于最大可能事故类型，计算最初的污染带长度，估算污染事故的最早扩散时间；

$\frac{L}{t}=\∝\·\ln \frac{M}{C_0\·t^{0.5}}-\mathrm{\β}\·u_0-k---\left(3\right)$

其中，L-不同时间节点污染带长度，m；t-污染带扩散时间，h；C_o-污染带监测的最高浓度，mg/L；u_o-污染带长度扩展速率，m/h；α、β和k-河流常数；

3.事故最大扩散范围分析。事故最大扩散范围是基于事故发现地点，涉嫌风险源在河流上游可能存在的最远位置。事故最大扩散范围是筛选涉嫌风险源的重要依据，只有纳入到事故危害范围内的风险源才会参与到后续风险源的排查中，这有利于避免排查范围的扩大化。依据步骤四第二步确定的事故最早扩散时间，按照式(4)反推事故的最大扩散距离；

$X=u_{x}t+0.63\sqrt{\frac{t}{H}\ln \left(\frac{M}{1.12\sqrt{H\·t}\·B\·C_0}\right)}---\left(4\right)$

其中，X-污染物泄漏最大长度，m；u_x-河流流速，m/s；t-污染带扩散时间，h；M-污染物泄露总量；C₀-污染带监测的最高浓度，mg/L；B-河流平均宽度，m；H-河流平均深度，m；

在实际工作中，事故最大扩散范围的排查还应与已构建的河流监测体系、河流感性指标变化相结合，并对事故最大扩散范围进行修正。此外，根据污染带迁移过程中在线监测或流域现场调查获取的资料，同样可根据式(4)计算事故最近扩散距离，即确定从事故发现地到涉嫌风险源间的最近距离。通过将事故最大扩散范围中扣除事故最近扩散距离即可得到涉嫌风险源的最大可能事故发生范围。

步骤五、事故涉嫌风险源二次筛选

基于步骤四中获取的最大可能事故发生范围，依据流域风险源的WebGis分布***，通过风险源及污水排放口坐标分布，完成事故发现地点上游涉嫌风险源的二次筛选，进一步缩小事故涉嫌风险源范围。当筛选后流域涉嫌风险源数小于5家以后，可直接超越到步骤七进行涉嫌风险源的现场取证及验证。

步骤六、涉嫌风险源事故场景符合性分析

1.风险源事故场景模拟。基于二次筛选后涉嫌风险源的风险物质储量、存贮形式、状态等信息，模拟该风险源的最大可信事故场景，模拟分析事故发生后涉嫌风险源内污染物泄漏时间、泄漏量和污染带分布的关系。最大可信事故源强可依据单位时间内液态有毒有害危险化学品进入地表水体中的含量，采用下式进行计算：

Q_源强＝λ×Q_max×P_泄漏比例/t         (5)

其中，Q_max为有毒有害危险化学品最大储量，t；P_泄漏比例为有毒有害危险化学品的泄漏比例；λ为液态泄露污染物自流入敏感地表水体的比例；t为泄漏时间，min；

对于水污染事故风险源导致的地表水体污染，运用二维扩散模型计算事故危害范围，表达式如下：

其中，ρ(x，y，t)为泄漏点下游距离x处、t时溶解态污染物浓度，mg/L；Q_源强为污染源源强，g；D_x/D_y为纵向及横向离散系数，m²/s；u_x/u_y为横向及纵向流速，m/s；h为平均水深，m；K为降解速率，g/g·d；

2.事故场景符合性分析。对比实际事故监测结果与涉嫌风险源事故场景模拟结果，按照公式(7)计算涉嫌风险源事故场景的符合性：

$f(X,M,t)=\mathrm{\γ}\frac{X}{X_0}+\mathrm{\δ}\frac{M}{M_0}+\mathrm{\ϵ}\frac{t}{t_0}---\left(7\right)$

其中，f为事故场景符合因子；X，M，t(X₀，M₀，t₀)分别为模拟事故场景(或实际事故)的扩散距离、泄漏量和扩散时间；γ、δ和ε为权重因子，且γ+δ+ε＝1；

3.涉嫌风险源事故场景符合性排序：基于步骤六第二步计算的事故场景符合因子大小，建立事故涉嫌风险源的符合性排序清单。

步骤七、涉嫌风险源现场取证及验证

依据建立的事故涉嫌风险源的符合性排序清单，依次开展涉嫌风险源的现场调查、取证及后果验证；获取涉嫌风险源中与事故化学品相同的材料，采取稳定同位素分析法，分析化学品的组成或特征指纹图库；通过样品成分比对分析，确定事故风险源，明确企业责任。

Claims (4)

1.一种突发水环境污染事故涉嫌风险源的溯源方法，其特征在于：其具体操作步骤如下，

步骤一、风险源信息数据库构建

第一步：建立基于WebGis的流域风险源分布软件***；

***包括风险源名称、经济产值、所属行业、风险源坐标、污水排放口坐标、受纳水体等风险源基础信息；

第二步：建立定时更新的风险源风险物质数据库；

通过企业定期申报、不定期现场检查等更新手段，建立流域风险源风险物质数据库，主要包括：有毒有害原料、辅料、产品和废弃物，并定期更新风险源风险物质名称、储量、存储方式、监控手段、存储场所等信息；

第三步：建立流域高风险化学品组成信息库；

采用环境风险源辨识法，筛选流域内风险源具有的高毒有机物、重金属和油类等高风险化学品，采用内标法测定此类风险化学品的典型标志物，建立流域高风险化学品组成信息库；

步骤二、污染事故应急监测及信息获取

第一步：获取事故污染物组成；

接到举报或发现水质异常后，开展事故应急监测，测定事故污染物主要组成；

第二步：获取污染带变化信息；

采用流动监测(车、船等)、简易监测、低空航测、遥感等监测手段，在事故发现地点上下游获取污染带扩散速率、长度、宽度，以及浓度分布；

第三步：计算污染物泄漏总量；

基于第二步获取的污染带变化信息，利用公式(1)计算事故泄漏的污染物总量；

$M={\mathrm{\Σ}}_1^n(C_i\·L_i)---\left(1\right)$

其中，M-污染物泄漏总量，t；C_i-监测点位浓度，mg/L；L_i-监测点位间隔距离，km；n-监测点位个数，通常大于10；

步骤三、事故涉嫌风险源初步筛选

依据从步骤二中获取的污染物种类及泄漏量等事故信息，利用排除法从流域风险源风险物质数据库和流域高风险化学品组成信息库中筛选事故涉嫌风险源，初步确定嫌疑风险源范围；

步骤四、突发事故涉嫌风险源定位

第一步：事故最大可能类型分析；

根据污染物的理化性质、污染带物质分布特点，判定事故的最大可能类型，分析事故是否属于废水非正常排放事故、液态化学品直接泄漏事故或液态化学次生泄漏事故等；依据案例分析结果，预测事故的初始污染带长度范围；

第二步：事故最早扩散时间分析；

依据河流基本参数及步骤二中获得的污染带变化数据，采用式(2)计算污染带长度与扩散时间的关系，并确定污染事故的最早扩散时间；

$\frac{L}{t}=\∝\·\ln \frac{M}{C_0\·t^{0.5}}-\mathrm{\β}\·u_0-k---\left(2\right)$

其中，L-不同时间节点污染带长度，m；t-污染带扩散时间，h；C_o-污染带监测的最高浓度，mg/L；u_o-污染带长度扩展速率，m/h；α、β和k-河流常数；

第三步：事故最大扩散范围分析；

依据事故最早扩散时间反推事故最大可能发生地点，按照式(3)计算事故的最大扩散距离；

$X=u_{x}t+0.63\sqrt{\frac{t}{H}\ln \left(\frac{M}{1.12\sqrt{H\·t}\·B\·C_0}\right)}---\left(3\right)$

其中，X-污染物泄漏最大长度，m；u_x-河流流速，m/s；t-污染带扩散时间，h；M-污染物泄露总量；C₀-污染带监测的最高浓度，mg/L；B-河流平均宽度，m；H-河流平均深度，m；

步骤五、事故涉嫌风险源二次筛选

基于步骤四中获得的事故最大扩散范围，依据风险源分布的WebGis***，对第一次筛选后涉嫌风险源进行第二次筛选，进一步缩小事故涉嫌风险源范围；

步骤六、涉嫌风险源事故场景符合性分析

第一步：风险源事故场景模拟；

获取步骤五中二次筛选后涉嫌风险源的环境风险源信息，根据风险源最大可能事故场景分析，模拟风险源同类事故中扩散时间、泄漏量和扩散距离的关系；

第二步：事故场景符合性分析；

对比实际事故监测结果与涉嫌风险源事故场景模拟结果，按照公式(4)计算涉嫌风险源事故场景的符合性：

$f(X,M,t)=\mathrm{\γ}\frac{X}{X_0}+\mathrm{\δ}\frac{M}{M_0}+\mathrm{\ϵ}\frac{t}{t_0}---\left(4\right)$

其中，f(X，M，t)-事故场景符合性因子；X，M，t(X₀，M₀，t₀)-分别为涉嫌风险源事故模拟场景(及实际监测)的事故扩散距离、泄漏量和扩散时间；γ、δ和ε-权重因子，且γ+δ+ε＝1；

第三步：涉嫌风险源事故场景符合性排序；

基于第二步计算出的事故场景符合因子大小，建立事故涉嫌风险源的符合性排序清单；步骤七、涉嫌风险源现场取证及验证

依据建立的事故涉嫌风险源的符合性排序清单，依次开展涉嫌风险源的现场调查、取证及后果验证；获取涉嫌风险源中与事故化学品相同的材料，采取稳定同位素分析法，分析化学品的组成或特征指纹图库；通过样品成分比对分析，确定事故风险源，明确企业责任。

2.如权利要求1所述的一种突发水环境污染事故涉嫌风险源的溯源方法，其特征在于：步骤三中事故涉嫌风险源的初步筛选应先排除不存储、不使用及不生产事故化学品的风险源，再排除事故化学品存储量、使用及生产量超过事故泄漏量的风险源。

3.如权利要求1所述的一种突发水环境污染事故涉嫌风险源的溯源方法，其特征在于：当经步骤三(或步骤五)排除后涉嫌风险源数小于5时，可直接超越到步骤七开展涉嫌风险源的现场调查、取证及验证。

4.如权利要求1所述的一种突发水环境污染事故涉嫌风险源的溯源方法，其特征在于：风险源事故场景模拟由公式(5)和公式(6)依次确定：

确定最大可信事故源强，依据单位时间内液态有毒有害危险化学品进入地表水体中的含量，采用下式进行计算：

Q_源强＝λ×Q_max×P_泄漏比例/t         (5)

其中，Q_max-有毒有害危险化学品最大储量，t；P_泄漏比例-有毒有害危险化学品的泄漏比例；λ-液态泄露污染物自流入敏感地表水体的比例；t-泄漏时间，min；

对于水污染事故风险源导致的地表水体污染，运用二维扩散模型计算事故危害范围，表达式如下：

其中，ρ(x，y，t)-泄漏点下游距离x处、t时溶解态污染物浓度，mg/L；Q_源强-污染源源强，g；D_x/D_y-纵向及横向离散系数，m²/s；u_x/u_y-横向及纵向流速，m/s；h-平均水深，m；K-降解速率，g/g·d。

Images