CN106339586A - 芥子气污染溯源方法及预报方法 - Google Patents

Abstract

本发明提供一种化学污染物溯源及实时预报方法，芥子气等化学毒剂弹在河流潜水打捞、挖掘等清理工作时，如果因泄漏、爆炸造成芥子气进入河流中，则会产生较大环境风险。因此，作业时首先要了解作业区的实时水文情况，在不利作业的水文情况发生前，采取停止作业、加强防范等措施；另外万一发生芥子气进入河流中时，必须精确了解污染态势，采用合适的措施。

Description

芥子气污染溯源方法及预报方法

一、技术领域

本发明涉及一种遗弃化学武器(主要涉及芥子气)河流化学风险精准预测模型及潜水作业环境风险实时预报技术，不仅可以对河流中化学武器的扩散风险进行精准预测，也可以对潜水作业风险进行实时预报，以保障环境和居民、操作人员的人身安全。

二、背景技术

侵华日军在中国东北地区遗留下大量的芥子气毒剂和其他类型的化学武器，曾引发了多起中毒事件,至今仍威胁中国的生态环境和居民生命安全。芥子气是一种危害性巨大的糜烂性化学毒剂，在正常气候下，仅0.2mg/l的浓度就可以使人受到毒害，且性质稳定，进入环境中不易自然分解，必须使用特殊手段进行处理。

目前，国内外学者对芥子气的研究主要集中在芥子气的识别、毒理机制、检测分析方法上，对芥子气在自然水体中的扩散模拟方面的研究尚未见报道。

本发明采用多维精细网格化数值模拟方法结合GIS对化学武器潜水挖掘作业河道流场进行模拟，针对多种风险情景进行了分析，计算出芥子气在爆炸、泄漏等不同情景时在水体中扩散的影响，界定了其污染在整个作业河流内随时间的变化在空间分布，评估作业河流的化学风险。

芥子气毒剂弹在河流潜水打捞、挖掘等清理工作时，如果因泄漏、爆炸造成芥子气进入河流中，则会产生较大环境风险。因此，作业时首先要了解作业区的实时水文情况，在不利作业的水文情况发生前，采取停止作业、加强防范等措施；另外万一发生芥子气进入河流中时，必须精确了解污染态势，采用合适的措施。

三、发明内容

基于上述技术问题，本发明根据其中一个发明目的提供一种芥子气污染溯源方法，包括：与GIS地理信息系统联网的主机系统，通过无线网络与主机系统通信的若干无线监测设备，无线监测设备不断获取监测区域的个中数据并及时对主机系统的数据进行实时的更改和修正，主机系统根据所获取的GIS数据建立实时预测模型库系统，所述实时预测模型库系统也不断根据所获取的无线监测设备对其中的预测模型进行实时的调整和分析，所述分析方法步骤如下：

步骤1：获取GIS地理信息系统中作业区整体水环境数据并进行基本属性情况调查；

步骤2：根据所获取的上述数据建立实时预测模型库系统，包括二、三维水流及水质模型；

步骤3：选取目标河段设置合适的控制点，置放可与实时预测模型库系统建立实时通信的无线监测设备，建立实时水文站点，实时发送所监测得到的数据，其中该无线监测设备可用于监测化学污染物芥子气的浓度；

步骤4：利用实时水文站点数据，结合上游远端水文站的信息，自动生成动态欧拉流场，并预测作业区的水文变化情况，包括流量、流速、浊度、水温等指标。

步骤5：根据水动力情况，模拟发生事故的污染态势分析；

步骤6：根据步骤4及5结果，形成水文及风险预报报告，提供指挥部指导现场作业。

基于本发明进一步的发明目的，其中还包括：步骤7：一旦发现存在芥子气泄漏，立刻对发现芥子气的无线监测设备进行准确定位，并对其下游水体实施截流，再结合水流参数和不同无线监测设备的数据进行综合分析定位，确认污染源的大致位置进行搜寻，对污染水体进行特殊处理。

根本本发明的另一发明目的，其中提供一种芥子气污染预报方法，包括：与GIS地理信息系统联网的主机系统，通过无线网络与主机系统通信的若干无线监测设备，无线监测设备不断获取监测区域的个中数据并及时对主机系统的数据进行实时的更改和修正，主机系统根据所获取的GIS数据建立实时预测模型库系统，所述实时预测模型库系统也不断根据所获取的无线监测设备对其中的预测模型进行实时的调整和分析，所述预报方法步骤如下：

步骤1：根据事故特征，设置风险情景模式：

(1)最大风险情景：发生爆炸情况下，芥子气以细微粒子的形式进入水体扩散；

(2)一般风险情景：发生泄漏等意外情况造成芥子气进入河流水体。

步骤2：利用GIS地理信息系统所获得的地理和水文数据，基于所述数据和水环境条件，建立三维污染扩散预测模型，所建立模型涉及以下算法：

(1)动量方程：

$\begin{array}{c}{\∂}_t\left(\mathrm{mHu}\right)+{\∂}_x\left(m_y\mathrm{Huu}\right)+{\∂}_y\left(m_x\mathrm{Hvu}\right)+{\∂}_z\left(\mathrm{mwu}\right)-(\mathrm{mf}+v{\∂}_x{m}_y-u{\∂}_y{m}_x)\mathrm{Hv}\\ =-m_{y}H{\∂}_x(\mathrm{g\δ}+p)-m_y({\∂}_{x}h-z{\∂}_{x}H){\∂}_{z}p+{\∂}_z\left(m{H}^{-1}{A}_v{\∂}_{z}u\right)+Q_u\end{array}---(8-2)$

${\∂}_t\left(\mathrm{mHv}\right)+{\∂}_x\left(m_y\mathrm{Huv}\right)+{\∂}_y\left(m_x\mathrm{Hvv}\right)+{\∂}_z\left(\mathrm{mwv}\right)+(\mathrm{mf}+v{\∂}_x{m}_y-u{\∂}_y{m}_x)\mathrm{Hu}$

$=-m_{x}H{\∂}_y(\mathrm{g\δ}+p)-m_x({\∂}_{y}h-z{\∂}_{y}H){\∂}_{z}p+{\∂}_z\left({\mathrm{mH}}^{-1}{A}_v{\∂}_{z}v\right)+Q_v---(8-3)$

(2)连续方程：

${\∂}_t\left(\mathrm{m\δ}\right)+{\∂}_x\left(m_{y}H{\∫}_0^1\mathrm{udz}\right)+{\∂}_y\left(m_{x}H{\∫}_0^1\mathrm{vdz}\right)=0---(8-4)$

(3)压强、密度、温盐状态方程：

${\∂}_{z}p=-\mathrm{gh}(\mathrm{\ρ}-{\mathrm{\ρ}}_0){{\mathrm{\ρ}}_0}^{-1}=-\mathrm{gHb}---(8-5)$

ρ＝ρ(p，S，T) (8-6)

$\begin{array}{c}{\∂}_t\left(\mathrm{mHS}\right)+{\∂}_x\left(m_y\mathrm{HuS}\right)+{\∂}_y\left(m_x\mathrm{HvS}\right)+{\∂}_z\left(\mathrm{mwS}\right)\\ ={\∂}_z\left(m{H}^{-1}{A}_b{\∂}_{z}S\right)+Q_S\end{array}---(8-7)$

$\begin{array}{c}{\∂}_t\left(\mathrm{mHT}\right)+{\∂}_x\left(m_y\mathrm{HuT}\right)+{\∂}_y\left(m_x\mathrm{HvT}\right)+{\∂}_z\left(\mathrm{mwT}\right)\\ ={\∂}_z\left(m{H}^{-1}{A}_b{\∂}_{z}T\right)+Q_T\end{array}---(8-8)$

(4)物质输运方程：

$\begin{array}{c}{\∂}_{t}C+{\∂}_x\left(\mathrm{uC}\right)+{\∂}_y\left(\mathrm{vC}\right)+{\∂}_z\left(\mathrm{wC}\right)\\ ={\∂}_x\left(K_x{\∂}_{x}C\right)+{\∂}_y\left(K_y{\∂}_{y}C\right)+{\∂}_z\left(K_z{\∂}_{z}C\right)+S_C\end{array}---(8-9)$

在上述式(8-2)—式(8-9)中：是时间；u、v、w分别是x、y、z三个方向的流速；mx、my是坐标变化张量(分别是度量张量对角元素的平方根)；f为科氏系数；δ是水位，h是河床高程，H为水深，且H＝h+ε；ρ是密度，S是盐度，T是水温，p为压力；A_v和A_b分别为垂向紊动粘性系数和扩散系数；Q_u和Q_v为动量方程的源汇项(包括水平扩散项等)；Q_S和Q_T项分别是盐度和水温的外部源汇项；b代表参考密度的偏移值；C为污染物浓度，K_v和K_H分别为垂向和水平的紊动扩散系数，S_C为源汇项。

步骤3：利用作业场水文实时监测数据，实时调整参数，并计算污染物的扩散情况；

步骤4：把步骤3结果反馈空间决策模型，并通过GIS可视化表达，通过空间分析，找到影响的河流范围及相关行政区域。

步骤5：根据受影响区域受到的影响程度、时间等要素，估计其受到的损失，形成风险图，并提出防范措施。

根据本发明另一发明目的，其中提供一种芥子气污染溯源后再进行芥子气污染预报的方法，先执行如前述方法的污染源溯源，再进行如前述方法的预测预报分析，对污染源附近或扩散方向的人群进行预报预警。

四、附图说明

图1是本发明的模型库基本结构；

图2是本发明实时预测模型库系统的主要功能；

图3是本发明预报分析方法的主要流程。

五、具体实施方式

实施例1：

下面结合附图1-2具体描述本发明预测化学武器(主要是芥子气)泄漏/污染并溯源的具体技术方案：

化学武器(主要是芥子气)污染溯源系统包括：与GIS地理信息系统联网的主机系统，通过无线网络与主机系统通信的若干无线监测设备，无线监测设备不断获取监测区域的个中数据并及时对主机系统的数据进行实时的更改和修正，主机系统根据所获取的GIS数据建立实时预测模型库系统，所述实时预测模型库系统也不断根据所获取的无线监测设备对其中的预测模型进行实时的调整和分析，所述分析方法步骤如下：

步骤1：获取GIS地理信息系统中作业区整体水环境数据并进行基本属性情况调查；

步骤2：根据所获取的上述数据建立实时预测模型库系统，包括二、三维水流及水质模型；

步骤3：选取目标河段设置合适的控制点，置放可与实时预测模型库系统建立实时通信的无线监测设备，建立实时水文站点，实时发送所监测得到的数据，其中该无线监测设备可用于监测化学污染物芥子气的浓度；

步骤4：利用实时水文站点数据，结合上游远端水文站的信息，自动生成动态欧拉流场，并预测作业区的水文变化情况，包括流量、流速、浊度、水温等指标。

步骤5：根据水动力情况，模拟发生事故的污染态势分析；

步骤6：根据步骤4及5结果，形成水文及风险预报报告，提供指挥部指导现场作业。

步骤7：一旦发现存在芥子气泄漏，立刻对发现芥子气的无线监测设备进行准确定位，并对其下游水体实施截流，再结合水流参数和不同无线监测设备的数据进行综合分析定位，确认污染源的大致位置进行搜寻，对污染水体进行特殊处理。

上述技术方案能达成以下技术效果：

1)研究充分国家各级水文、水质监测网络，结合控制段面，通过多模型结合的实时预测模型库系统，实现了作业区水文情况的精确控制；

2)实现了对作业场水文实时预报，实时对不利情况(水位、温度及浊度等指标不利于作业)进行预警，有效地控制了现场作业的风险，保障了作业人员的安全；

3)水文及风险预报报告可以为区域风险预警、风险防控提供明准确的信息，对区域预警方案制定、应急资源的有效调配具有明确的指导作用。

实施例2：

下面结合附图具体描述对河流中化学武器(主要是芥子气)例如芥子气的风险进行实时预报的技术方案：

化学武器(主要是芥子气)污染溯源系统包括：与GIS地理信息系统联网的主机系统，通过无线网络与主机系统通信的若干无线监测设备，无线监测设备不断获取监测区域的个中数据并及时对主机系统的数据进行实时的更改和修正，主机系统根据所获取的GIS数据建立实时预测模型库系统，所述实时预测模型库系统也不断根据所获取的无线监测设备对其中的预测模型进行实时的调整和分析，所述分析预报方法步骤如下

步骤1：根据事故特征，设置风险情景模式：

(1)最大风险情景：发生爆炸情况下，芥子气以细微粒子的形式进入水体扩散；

(2)一般风险情景：发生泄漏等意外情况造成芥子气进入河流水体。

步骤2：利用GIS地理信息系统所获得的地理和水文数据，基于所述数据和水环境条件，建立三维污染扩散预测模型，所建立模型涉及以下算法：

(1)动量方程：

$\begin{array}{c}{\∂}_t\left(\mathrm{mHu}\right)+{\∂}_x\left(m_y\mathrm{Huu}\right)+{\∂}_y\left(m_x\mathrm{Hvu}\right)+{\∂}_z\left(\mathrm{mwu}\right)-(\mathrm{mf}+v{\∂}_x{m}_y-u{\∂}_y{m}_x)\mathrm{Hv}\\ =-m_{y}H{\∂}_x(\mathrm{g\δ}+p)-m_y({\∂}_{x}h-z{\∂}_{x}H){\∂}_{z}p+{\∂}_z\left(m{H}^{-1}{A}_v{\∂}_{z}u\right)+Q_u\end{array}---(8-2)$

$\begin{array}{c}{\∂}_t\left(\mathrm{mHv}\right)+{\∂}_x\left(m_y\mathrm{Huv}\right)+{\∂}_y\left(m_x\mathrm{Hvv}\right)+{\∂}_z\left(\mathrm{mwv}\right)-(\mathrm{mf}+v{\∂}_x{m}_y-u{\∂}_y{m}_x)\mathrm{Hu}\\ =-m_{x}H{\∂}_y(\mathrm{g\δ}+p)-m_x({\∂}_{y}h-z{\∂}_{y}H){\∂}_{z}p+{\∂}_z\left(m{H}^{-1}{A}_v{\∂}_{z}v\right)+Q_v\end{array}---(8-3)$

(2)连续方程：

${\∂}_t\left(\mathrm{m\δ}\right)+{\∂}_x\left(m_{y}H{\∫}_0^1\mathrm{udz}\right)+{\∂}_y\left(m_{x}H{\∫}_0^1\mathrm{vdz}\right)=0---(8-4)$

(3)压强、密度、温盐状态方程：

${\∂}_{z}p=-\mathrm{gh}(\mathrm{\ρ}-{\mathrm{\ρ}}_0){{\mathrm{\ρ}}_0}^{-1}=-\mathrm{gHb}---(8-5)$

ρ＝ρ(p，S，T) (8-6)

$\begin{array}{c}{\∂}_t\left(\mathrm{mHS}\right)+{\∂}_x\left(m_y\mathrm{HuS}\right)+{\∂}_y\left(m_x\mathrm{HvS}\right)+{\∂}_z\left(\mathrm{mwS}\right)\\ ={\∂}_z\left(m{H}^{-1}{A}_b{\∂}_{z}S\right)+Q_S\end{array}---(8-7)$

$\begin{array}{c}{\∂}_t\left(\mathrm{mHT}\right)+{\∂}_x\left(m_y\mathrm{HuT}\right)+{\∂}_y\left(m_x\mathrm{HvT}\right)+{\∂}_z\left(\mathrm{mwT}\right)\\ ={\∂}_z\left(m{H}^{-1}{A}_b{\∂}_{z}T\right)+Q_T\end{array}---(8-8)$

(4)物质输运方程：

$\begin{array}{c}{\∂}_{t}C+{\∂}_x\left(\mathrm{uC}\right)+{\∂}_y\left(\mathrm{vC}\right)+{\∂}_z\left(\mathrm{wC}\right)\\ ={\∂}_x\left(K_x{\∂}_{x}C\right)+{\∂}_y\left(K_y{\∂}_{y}C\right)+{\∂}_z\left(K_z{\∂}_{z}C\right)+S_C\end{array}---(8-9)$

在上述式(8-2)—式(8-9)中：是时间；u、v、w分别是x、y、z三个方向的流速；mx、my是坐标变化张量(分别是度量张量对角元素的平方根)；f为科氏系数；δ是水位，h是河床高程，H为水深，且H＝h+ε；ρ是密度，S是盐度，T是水温，p为压力；A_v和A_b分别为垂向紊动粘性系数和扩散系数；Q_u和Q_v为动量方程的源汇项(包括水平扩散项等)；Q_S和Q_T项分别是盐度和水温的外部源汇项；b代表参考密度的偏移值；C为污染物浓度，K_v和K_H分别为垂向和水平的紊动扩散系数，S_C为源汇项。

步骤3：利用作业场水文实时监测数据，实时调整参数，并计算污染物的扩散情况；

步骤4：把步骤3结果反馈空间决策模型，并通过GIS可视化表达，通过空间分析，找到影响的河流范围及相关行政区域。

步骤5：根据受影响区域受到的影响程度、时间等要素，估计其受到的损失，形成风险图，并提出防范措施。

上述技术方案能达成以下技术效果：

1)通过应用三维水动力水质模型，对芥子气在河流中的扩散情况进行了预测分析，填补了国内外对芥子气在自然水体中的扩散模拟方面的研究空白。

2)利用三维模型分水层面的计算分析芥子气污染的扩散情况，可以为有效地控制防范芥子气风险提供精准的信息，支持现场作业的风险控制，大大提高现场工作能力。

3)通过模型与GIS相结合，研究制定了区域风险分布图，提高了作业、管理等相关人员对此风险的认识，对区域预警方案制定、应急资源的配置、应急决策有很重要的意义。

也可以将上述实施例1和实施例2结合在一起，形成实时监测和预防一体的技术方案。

虽然先前的描述和附图描述了本发明的优选实施例，但是可以理解：在不脱离本发明的精神的情况下，在此可以产生各种附加、修改和替换。本领域普通技术人员很清楚：在不脱离本发明的精神或本质特性的情况下，可以以其他特殊形式、结构、布置、比例、以及利用其他元件、材料和部件来实现本发明。本领域的技术人员将意识到：本发明可以使用发明实际中使用的结构、布置、比例、材料以及部件和其他的许多修改，这些修改在不脱离本发明的原理的情况下而特别适应于特殊环境和操作需求。因此，当前公开的实施例在所有方面应被理解为说明性的而非对其请求保护的范围的限制。

Claims (4)

1.一种芥子气污染溯源方法，包括：与GIS地理信息系统联网的主机系统，通过无线网络与主机系统通信的若干无线监测设备，无线监测设备不断获取监测区域的个中数据并及时对主机系统的数据进行实时的更改和修正，主机系统根据所获取的GIS数据建立实时预测模型库系统，所述实时预测模型库系统也不断根据所获取的无线监测设备对其中的预测模型进行实时的调整和分析，所述分析方法步骤如下：

步骤1：获取GIS地理信息系统中作业区整体水环境数据并进行基本属性情况调查；

步骤2：根据所获取的上述数据建立实时预测模型库系统，包括二、三维水流及水质模型；

步骤3：选取目标河段设置合适的控制点，置放可与实时预测模型库系统建立实时通信的无线监测设备，建立实时水文站点，实时发送所监测得到的数据，其中该无线监测设备可用于监测化学污染物芥子气的浓度；

步骤4：利用实时水文站点数据，结合上游远端水文站的信息，自动生成动态欧拉流场，并预测作业区的水文变化情况，包括流量、流速、浊度、水温等指标。

步骤5：根据水动力情况，模拟发生事故的污染态势分析；

步骤6：根据步骤4及5结果，形成水文及风险预报报告，提供指挥部指导现场作业。

2.如权利要求1的芥子气污染溯源方法，其中还包括：步骤7：一旦发现存在芥子气泄漏，立刻对发现芥子气的无线监测设备进行准确定位，并对其下游水体实施截流，再结合水流参数和不同无线监测设备的数据进行综合分析定位，确认污染源的大致位置进行搜寻，对污染水体进行特殊处理。

3.一种芥子气污染预报方法，包括：与GIS地理信息系统联网的主机系统，通过无线网络与主机系统通信的若干无线监测设备，无线监测设备不断获取监测区域的个中数据并及时对主机系统的数据进行实时的更改和修正，主机系统根据所获取的GIS数据建立实时预测模型库系统，所述实时预测模型库系统也不断根据所获取的无线监测设备对其中的预测模型进行实时的调整和分析，所述预报方法步骤如下：

步骤1：根据事故特征，设置风险情景模式：

(1)最大风险情景：发生爆炸情况下，芥子气以细微粒子的形式进入水体扩散；

(2)一般风险情景：发生泄漏等意外情况造成芥子气进入河流水体。

步骤2：利用GIS地理信息系统所获得的地理和水文数据，基于所述数据和水环境条件，建立三维污染扩散预测模型，所建立模型涉及以下算法：

(1)动量方程：

$\begin{array}{c}{\∂}_t\left(mHu\right)+{\∂}_x\left(m_{y}Huu\right)+{\∂}_y\left(m_{x}Hvu\right)+{\∂}_z\left(mwu\right)-(mf+v{\∂}_x{m}_y-u{\∂}_y{m}_x)Hv\\ =-m_{y}H{\∂}_x(g\mathrm{\δ}+p)-m_y({\∂}_{x}h-z{\∂}_{x}H){\∂}_{z}p+{\∂}_z\left({\mathrm{mH}}^{-1}{A}_v{\∂}_{z}u\right)+Q_u\end{array}---(8-2)$

$\begin{array}{c}{\∂}_t\left(mHv\right)+{\∂}_x\left(m_{y}Huv\right)+{\∂}_y\left(m_{x}Hvv\right)+{\∂}_z\left(mwv\right)-(mf+v{\∂}_x{m}_y-u{\∂}_y{m}_x)Hu\\ =-m_{x}H{\∂}_y(g\mathrm{\δ}+p)-m_x({\∂}_{y}h-z{\∂}_{y}H){\∂}_{z}p+{\∂}_z\left({\mathrm{mH}}^{-1}{A}_v{\∂}_{z}v\right)+Q_v\end{array}---(8-3)$

(2)连续方程：

$a_t\left(m\mathrm{\δ}\right)+{\∂}_x\left(m_{y}H{\∫}_0^{1}udz\right)+{\∂}_y\left(m_{x}H{\∫}_0^{1}vdz\right)=0---(8-4)$

(3)压强、密度、温盐状态方程：

${\∂}_{z}p=-gH(\mathrm{\ρ}-{\mathrm{\ρ}}_0){{\mathrm{\ρ}}_0}^{-1}=-gHb---(8-5)$

ρ＝ρ(p，S，T) (8-)

$\begin{array}{c}{\∂}_t\left(mHS\right)+{\∂}_x\left(m_{y}HuS\right)+{\∂}_y\left(m_{x}HvS\right)+{\∂}_z\left(mwS\right)\\ ={\∂}_z\left({\mathrm{mH}}^{-1}{A}_b{\∂}_{z}S\right)+Q_S\end{array}---(8-7)$

$\begin{array}{c}{\∂}_t\left(mHT\right)+{\∂}_x\left(m_{y}HuT\right)+{\∂}_y\left(m_{x}HvT\right)+{\∂}_z\left(mwT\right)\\ ={\∂}_z\left({\mathrm{mH}}^{-1}{A}_b{\∂}_{z}T\right)+Q_T\end{array}---(8-8)$

(4)物质输运方程：

$\begin{array}{c}{\∂}_{t}C+{\∂}_x\left(uC\right)+{\∂}_y\left(vC\right)+{\∂}_z\left(wC\right)\\ ={\∂}_x\left(K_x{\∂}_{x}C\right)+{\∂}_y\left(K_y{\∂}_{y}C\right)+{\∂}_z\left(K_z{\∂}_{z}C\right)+S_C\end{array}---(8-9)$

在上述式(8-2)—式(8-9)中：是时间；u、v、w分别是x、y、z三个方向的流速；mx、my是坐标变化张量(分别是度量张量对角元素的平方根)；f为科氏系数；δ是水位，h是河床高程，H为水深，且H＝h+ε；ρ是密度，S是盐度，T是水温，p为压力；A_v和A_b分别为垂向紊动粘性系数和扩散系数；Q_u和Q_v为动量方程的源汇项(包括水平扩散项等)；Q_S和Q_T项分别是盐度和水温的外部源汇项；b代表参考密度的偏移值；C为污染物浓度，K_v和K_H分别为垂向和水平的紊动扩散系数，S_C为源汇项。

步骤3：利用作业场水文实时监测数据，实时调整参数，并计算污染物的扩散情况；

步骤4：把步骤3结果反馈空间决策模型，并通过GIS可视化表达，通过空间分析，找到影响的河流范围及相关行政区域。

步骤5：根据受影响区域受到的影响程度、时间等要素，估计其受到的损失，形成风险图，并提出防范措施。

4.一种芥子气污染溯源后再进行芥子气污染预报的方法，先执行如权利要求1所述的污染源溯源，再进行如权利要求3所述的预测预报分析，对污染源附近或扩散方向的人群进行预报预警。