CN111027262A - 一种河流突发性水污染分析与风险图编制方法 - Google Patents
一种河流突发性水污染分析与风险图编制方法 Download PDFInfo
- Publication number
- CN111027262A CN111027262A CN201911125035.7A CN201911125035A CN111027262A CN 111027262 A CN111027262 A CN 111027262A CN 201911125035 A CN201911125035 A CN 201911125035A CN 111027262 A CN111027262 A CN 111027262A
- Authority
- CN
- China
- Prior art keywords
- water
- concentration
- pollutant
- time
- pollution
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000003911 water pollution Methods 0.000 title claims abstract description 54
- 238000000034 method Methods 0.000 title claims abstract description 41
- 238000004458 analytical method Methods 0.000 title claims abstract description 21
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 137
- 238000004088 simulation Methods 0.000 claims abstract description 36
- 230000008878 coupling Effects 0.000 claims abstract description 10
- 238000010168 coupling process Methods 0.000 claims abstract description 10
- 238000005859 coupling reaction Methods 0.000 claims abstract description 10
- 238000003903 river water pollution Methods 0.000 claims abstract description 9
- 239000003344 environmental pollutant Substances 0.000 claims description 84
- 231100000719 pollutant Toxicity 0.000 claims description 84
- 238000009792 diffusion process Methods 0.000 claims description 27
- 230000008859 change Effects 0.000 claims description 19
- 239000003086 colorant Substances 0.000 claims description 16
- 230000008569 process Effects 0.000 claims description 16
- 238000012545 processing Methods 0.000 claims description 9
- 238000010586 diagram Methods 0.000 claims description 8
- 239000000126 substance Substances 0.000 claims description 8
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 7
- 239000001301 oxygen Substances 0.000 claims description 7
- 229910052760 oxygen Inorganic materials 0.000 claims description 7
- 238000003672 processing method Methods 0.000 claims description 6
- 230000002262 irrigation Effects 0.000 claims description 5
- 238000003973 irrigation Methods 0.000 claims description 5
- 239000003621 irrigation water Substances 0.000 claims description 5
- 238000011160 research Methods 0.000 claims description 5
- 239000010865 sewage Substances 0.000 claims description 5
- 239000002352 surface water Substances 0.000 claims description 5
- 206010039203 Road traffic accident Diseases 0.000 claims description 4
- 230000005540 biological transmission Effects 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 239000000356 contaminant Substances 0.000 claims description 3
- 238000010790 dilution Methods 0.000 claims description 3
- 239000012895 dilution Substances 0.000 claims description 3
- 239000006185 dispersion Substances 0.000 claims description 3
- 230000007613 environmental effect Effects 0.000 claims description 3
- 238000013508 migration Methods 0.000 claims description 3
- 230000005012 migration Effects 0.000 claims description 3
- 238000012544 monitoring process Methods 0.000 claims description 3
- 238000012502 risk assessment Methods 0.000 claims description 3
- 238000011144 upstream manufacturing Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 2
- 230000015556 catabolic process Effects 0.000 claims description 2
- 238000006731 degradation reaction Methods 0.000 claims description 2
- 238000009877 rendering Methods 0.000 claims description 2
- 230000003068 static effect Effects 0.000 claims description 2
- 238000010276 construction Methods 0.000 description 3
- 230000006378 damage Effects 0.000 description 3
- 238000007726 management method Methods 0.000 description 3
- 238000012800 visualization Methods 0.000 description 3
- XKMRRTOUMJRJIA-UHFFFAOYSA-N ammonia nh3 Chemical compound N.N XKMRRTOUMJRJIA-UHFFFAOYSA-N 0.000 description 2
- 230000001788 irregular Effects 0.000 description 2
- 230000002265 prevention Effects 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000011109 contamination Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000003651 drinking water Substances 0.000 description 1
- 235000020188 drinking water Nutrition 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/18—Water
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/18—Water
- G01N33/1806—Biological oxygen demand [BOD] or chemical oxygen demand [COD]
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Engineering & Computer Science (AREA)
- Food Science & Technology (AREA)
- Medicinal Chemistry (AREA)
- Physics & Mathematics (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Biodiversity & Conservation Biology (AREA)
- Biomedical Technology (AREA)
- Emergency Medicine (AREA)
- Molecular Biology (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种河流突发性水污染分析与风险图编制方法，其特征在于：所述分析方法的步骤为：1)突发性水污染情景设定；2)建立水动力水质耦合模型；3)模拟结果分析；4)水污染风险图绘制。本发明方法集风险源识别，数值模拟及事故预警为一体，精确地模拟了水污染指标如区域受影响的时间及污染团到达特定断面的时间，为水环境管理者制定应急预案及事故后处理提供技术支持。
Description
技术领域
本发明属于水环境的应急防灾领域，具体涉及一种河流突发性水污染分析与风险图编制方法。
背景技术
随着社会经济的不断发展，城市规划建设的进程不断加快，我国的水环境问题日益恶化，尤其是突发性水污染事件频现，对水环境和供水系统造成了严重、恶劣的影 响。由于突发性水污染事故具有不确定性，事故预测和风险规避有一定的难度，且事 故一旦发生，影响范围极广，危害极大。因此，建立完善的突发性水污染状况的应急 体系，完成水污染风险图的绘制成为水污染应急防治亟待解决的重点问题。
突发性水污染问题具有事故发生的不确定性、危害的紧急性、影响的长期性等特点，现有研究主要通过风险源识别与评价、建立水质模型、水污染预警等方法来模拟 事故发生的过程及危害，取得了一定成果，但也存在一定的局限性：①传统的数值模 拟研究重点关注浓度的变化，对其他水污染指标如区域受影响的时间及污染团到达特 定断面的时间关注较少，应急预案的实施效率较低；②一维水质模型不能反映污染物 横向的浓度分布，且模拟的精度不足；③以往绘制的污染风险图中包含信息不全面， 不规范，且可视化效果较差。
本发明将针对以上问题，提供一种基于一、二维水动力水质耦合模型的丰、枯水时期的河流突发性水污染的模拟分析与风险图编制方法。
通过对公开专利文献的检索，并未发现与本专利申请相似的公开专利文献。
发明内容
本发明的目的在于克服现有技术的不足，，提供一种河流突发性水污染分析与风险图编制方法，该方法集风险源识别，数值模拟及事故预警为一体，精确地模拟了水 污染指标如区域受影响的时间及污染团到达特定断面的时间，为水环境管理者制定应 急预案及事故后处理提供技术支持；绘制的水污染风险图能够实现水污染信息的可视 化，提升我国水污染风险管理水平，促进我国水环境治理和生态文明建设。
本发明解决其技术问题是通过以下技术方案实现的：
一种河流突发性水污染分析与风险图编制方法，其特征在于：所述方法的步骤为：
1)突发性水污染情景设定：根据对突发性水污染特点及影响分析，设定两种突 发性水污染事件，分别为企业污水泄露和交通运输事故，研究水域为河流干流；
2)建立水动力水质耦合模型：水质模型用于模拟污染物在水体中迁移、扩散及 衰减转化过程，运用一、二维模型对步骤1)中的情景进行模拟，分析河道水污染风 险；
一维水动力水质模型
水动力基本方程为：
水质模型基本方程为：
其中：C为整个断面的平均浓度，单位为mg/m3；
Q为流量，单位为m3/s；
A为过水断面面积，单位为m2；
u为断面平均流速，单位为m/s；
D为纵向弥散系数，单位为m2/s；
S为各种源和漏的代数和；
b)二维水动力水质模型
水流连续方程为：
水流动量方程为：
其中：t为时间；
x、y为空间坐标系；
η为水位；
h为静止水深；
u为X方向的流速；
v为Y方向的流速；
pa为当地气压；
ρ为水密度；
ρ0为参考水密度；
f为科里奥利参量；
Txx、Txy、Tyx、Tyy为水平粘滞应力；
Isx、Ibx、Ixy、Iby为水流切应力分量；
us、vs为源汇项水流流速；
水质模型基本方程为：
其中：C为污染物的浓度，单位为mg/l；
t为时间，单位为d；
ux、uy分别为x、y方向的水流速度，单位为m/s；
Dx、Dy分别为x、y方向的扩散(离散)系数，m2/s；
K为污染物的降解速率常数，单位为s-1；
3)模拟结果分析：
a)一维模拟结果分析：
①取水口污染物分析：河流干流沿岸一般分布着大量的取水口，在假定污染源位置的前提下，根据一维水动力水质耦合模型模拟各取水口的浓度变化过程线；
根据《农田灌溉水质标准》GB5084-2005规范中对农田灌溉用水水质基本控制项目浓度的规定，化学需氧量对于水作作物应小于150mg/l，旱作作物要控制在200mg/l 以内。本发明约定当取水口处COD浓度首次达到150mg/l的时间为污染物到达时间， 当取水口处COD浓度下降到小于150mg/l的时刻与到达时刻的时间差，即为影响时 间；
②河道污染风险分析：根据《地表水环境质量标准》中对水质类别的划分，Ⅰ类、 Ⅱ类水中COD的浓度小于15mg/l、Ⅲ类水COD浓度介于15mg/l-20mg/l之间，Ⅳ 类水COD浓度介于20mg/l-30mg/l之间，Ⅴ类水中COD浓度介于30mg/l-40mg/l， 大于40mg/l为劣Ⅴ类水；在设定突发性水污染事故的情景下，可以绘制不同位置不 同时刻的COD浓度变化曲线；
b)二维模拟结果分析：在设定突发性水污染事故的情景下，二维模型考虑了横 向扩散，在水动力条件下，污染团沿着主河槽向下游迁移扩散，模拟污染物平面上的 浓度分布，绘制不同位置不同时刻的COD浓度变化过程平面图；
4)水污染风险图绘制：参照洪水风险图的种类，将河流水污染风险图分为污染 物浓度分布图、污染物到达时间图、影响时间图等三类，并在其上叠加监测站网、污 染源、取水工程信息及水源地分布等信息；将计算得到的污染物浓度分布、污染物到 达时间及污染物影响时间矢量化，与其他的矢量图层叠加，最后还要对图中的符号、 文字及颜色等进行设计、渲染；
(1)一维结果的矢量化：模拟计算的结果以沿程断面浓度变化的形式表征污染 物的扩散及传输规律，提取不同断面的污染物最大浓度及不同时刻的浓度，并认为浓 度沿河道中心线分布；经过差值加密的断面与河道中心线矢量图层叠加，并给各个交 点赋予浓度值。结合浓度分级标准，以浓度值字段为对象，将浓度值进行分级，并与 河道面矢量图层叠加，分级赋予不同的颜色，形成最大浓度的空间分布图及不同时刻 的污染物浓度分布图；
根据模拟的结果分析各个河道断面受到污染的起始时刻，并按照一定的间隔将时间分级，在GIS中处理河道图层shp，按照分级点的控制距离将不同的时间类别赋予 不同的颜色，即得到到达时间分布图层。影响时间分布层的处理方法与到达时间处理 方法相同；
(2)二维结果的矢量化：二维模拟结果存储在沿河道分布的网格中，每个网格 点按照时间序列存储了该点的空间位置信息和污染物的浓度值；根据网格的坐标信 息，创建空间矢量图层；污染物浓度分布是空间中的连续场，为了实现污染物浓度扩 散的连续二维场，在GIS中对浓度矢量图层进行IDW差值处理，按照分级标准将不 同的分级赋予不同的颜色，即得到污染物浓度分布矢量图；通过C++编程，提取分析 污染物的到达时间及影响历时，以网格编号作为不同文件之间转换的连接的唯一属 性，将分析提取的文件处理成与河道网格对应的矢量图层，最后按照分级标准，赋予 不同的颜色，表征污染物影响的变化；
(3)根据步骤3)中的一、二维模拟结果，分别进行矢量化，绘制最大浓度分 布图、不同时刻浓度分布图、污染物到达时间及影响时间分布图等不同种类的突发性 水污染风险图。
而且，所述情景设置的参数为：
(a)水文条件：考虑到不同时期上游来流量大小不同，污染物扩散和稀释的能 力不同，对下游各水取水口的应急时间影响也较大；设定的两种水文条件为枯水期和 丰水期的流量过程；
(b)水质条件：根据当地水环境现状的分析，确定当地水质状况，设定浓度本 底值。
本发明的优点和有益效果为：
1、本发明一种河流突发性水污染分析与风险图编制方法，集风险源识别，数值 模拟及事故预警为一体，精确地模拟了水污染指标如区域受影响的时间及污染团到达 特定断面的时间，为水环境管理者制定应急预案及事故后处理提供技术支持。
2、本发明一种河流突发性水污染分析与风险图编制方法，在一维水动力水质耦合模型的基础上，建立了二维水动力水质耦合模型，模拟了污染物横、纵向上的浓度 分布，反映了污染物扩散的客观规律。
3、本发明一种河流突发性水污染分析与风险图编制方法，同时基于GIS平台绘 制的水污染风险图能够实现水污染信息的可视化，提升我国水污染风险管理水平，促 进我国水环境治理和生态文明建设。
附图说明
图1为本发明的流程框图。
具体实施方式
下面通过具体实施例对本发明作进一步详述，以下实施例只是描述性的，不是限定性的，不能以此限定本发明的保护范围。
如图1所所示，一种河流突发性水污染分析与风险图编制方法，其特征在于：所 述方法的步骤为：
1)突发性水污染情景设定：根据对突发性水污染特点及影响分析，突发性水污 染一般是由企业污水泄露、交通运输事故、输油管道破裂等因素引起。本发明设定两 种突发性水污染事件，分别为企业污水泄露和交通运输事故，研究水域为河流干流。 通过对污染物原因的分析，河道主要污染物为COD及氨氮，考虑到农业灌溉用水对 氨氮的浓度没有限制，所以本发明以COD作为代表污染物。设定上述两种方案污染 源中含有大量的COD。其中企业排放的污水中COD浓度为5000mg/l，排放持续时间 2h；而运输化学物品的车辆泄露的污染物中COD浓度设定为8000mg/l，排放持续时 间为0.5h。
情景设置中：
(1)水文条件：考虑到不同时期上游来流量大小不同，污染物扩散和稀释的能 力不同，对下游各水取水口的应急时间影响也较大。本发明设定两种水文条件：枯水 期和丰水期的流量过程。
(2)水质条件：据黄河水环境现状的分析，黄河整体水质保持在Ⅲ类水以上， 设定浓度本底值COD为20mg/l。
2)建立水动力水质耦合模型：水质模型用于模拟污染物在水体中迁移、扩散及 衰减转化过程。本发明主要运用一、二维模型对上述不同情景进行模拟，分析河道水 污染风险。
(1)一维水动力水质模型
①本发明利用有限差分法中的显格式求解一维水动力基本方程。具体过程如下：
其中：
为已知层j上相邻3点的带权平均数。设权重分配系数为∝， 则水动力方程可离散为：
水质模型基本方程
针对较短河段，一般假定水流处于稳定状态，断面沿程变化不大。所以，水质模 型基本方程可简化为如下方程：
式中：K1为污染物衰减率系数，单位为1/d。
应用隐式差分方法对方程进行求解，f表示水质要素变化函数
将上面三式带入水质方程得到差分方程：
通过边界条件、初始条件对上式求解。
(2)二维水动力方程
①实际模拟河流的边界不规则，有限体积法采用非结构化网格，在解决不规则边界时，能够真实的模拟各种水力要素。本发明采用有限体积法对方程组进行离散。
利用高斯定理，将体积分转化成面积分，可得到有限体积法的基本方程
式中，V为每一个三角形网格，n为单位边界A的外法向单位向量，dA为面积分 微元，dV为体积分微元。
对控制体单元取平均，上式可化为
式中，△V为每个单元的体积，m为总单元总数，Aj为j的面积，Fn.j为法向通 量。
②水质方程模型参数
水质参数是影响模型精度决定水质的重要因素。经过概化的模型需要确定的横向扩散系数和纵向扩散系数。其中横向扩散系数是二维模型区别于一维模型的重要参 数，可采用天然河流横向扩散系数常用的经验公式：
Dy＝∝Hu*
式中：H为河流的平均水深，m；
u*为摩阻流速，m/l；
i为河流的平均水力坡度；
g为重力加速度，m2/s；
∝—量纲为1的横向扩散系数。
3)模拟结果分析：
(1)一维模拟结果分析
①取水口污染物分析
河流干流沿岸一般分布着大量的取水口，在假定污染源位置的前提下，一维水动力水质耦合模型可以模拟出各取水口的浓度变化过程线。
根据《农田灌溉水质标准》GB5084-2005规范中对农田灌溉用水水质基本控制项目浓度的规定，化学需氧量对于水作作物应小于150mg/l，旱作作物要控制在200mg/l 以内。本发明约定当取水口处COD浓度首次达到150mg/l的时间为污染物到达时间， 当取水口处COD浓度下降到小于150mg/l的时刻与到达时刻的时间差，即为影响时 间。
②河道污染风险分析
本发明设置的本底浓度值为Ⅲ类水中COD的浓度20mg/l，即认为河道某处的 COD浓度首次大于20mg/l的时间为到达时间。根据污染物扩散的规律，COD的浓度 会先升高，到达峰值浓度，随后浓度会下降，当COD浓度下降到小于20mg/l，即认 为此时刻为污染物离开的时间，该时间与到达时间的差值，即为影响时间。
由于黄河沿岸的取水口都用于农业灌溉，农业灌溉对水质的要求远远低于饮用水的标准，低于地表水水质标准要求。根据《地表水环境质量标准》中对水质类别的划 分，Ⅰ类、Ⅱ类水中COD的浓度小于15mg/l、Ⅲ类水COD浓度介于15mg/l-20mg/l 之间，Ⅳ类水COD浓度介于20mg/l-30mg/l之间，Ⅴ类水中COD浓度介于30mg/l- 40mg/l，大于40mg/l为劣Ⅴ类水。在设定突发性水污染事故的情景下，可以绘制不 同位置不同时刻的COD浓度变化曲线。
(2)二维模拟结果分析
在设定突发性水污染事故的情景下，除了绘制不同位置不同时刻的COD浓度变 化曲线外，二维模型考虑了横向扩散，在水动力条件下，污染团沿着主河槽向下游迁 移扩散，模拟污染物平面上的浓度分布，绘制不同位置不同时刻的COD浓度变化过 程平面图。
4)水污染风险图绘制：水污染风险图是反应河流突发性水污染风险专题图；参 照洪水风险图的种类，将河流水污染风险图分为污染物浓度分布图、污染物到达时间 图、影响时间图等三类，并在其上叠加监测站网、污染源、取水工程信息及水源地分 布等信息。
水污染风险图的绘制主要是在污染源识别、污染物分析、情景设定、水污染事件数值模拟基础上进行研究。一、二维水质模型模拟结果的表现形式不同，绘制水污染 风险图的技术方法也存在一定的差异。但主要的思路相同，将计算得到的污染物浓度 分布、污染物到达时间及污染物影响时间矢量化，与其他的矢量图层叠加。最后还要 对图中的符号、文字及颜色等进行设计、渲染等。
(1)一维结果的矢量化
模拟计算的结果以沿程断面浓度变化的形式表征污染物的扩散及传输规律，提取不同断面的污染物最大浓度及不同时刻的浓度，并认为浓度沿河道中心线分布。经过 差值加密的断面与河道中心线矢量图层叠加，并给各个交点赋予浓度值。结合浓度分 级标准，以浓度值字段为对象，将浓度值进行分级，并与河道面矢量图层叠加，分级 赋予不同的颜色，形成最大浓度的空间分布图及不同时刻的污染物浓度分布图。
根据模拟的结果分析各个河道断面受到污染的起始时刻，并按照一定的间隔将时间分级，在GIS中处理河道图层shp，按照分级点的控制距离将不同的时间类别赋予 不同的颜色，即得到到达时间分布图层。影响时间分布层的处理方法与到达时间处理 方法相同。
(2)二维结果的矢量化
二维模拟结果存储在沿河道分布的网格中，每个网格点按照时间序列存储了该点的空间位置信息和污染物的浓度值。根据网格的坐标信息，创建空间矢量图层。污染 物浓度分布是空间中的连续场，为了实现污染物浓度扩散的连续二维场，在GIS中 对浓度矢量图层进行IDW差值处理，按照分级标准将不同的分级赋予不同的颜色， 即得到污染物浓度分布矢量图。通过C++编程，提取分析污染物的到达时间及影响历 时，以网格编号作为不同文件之间转换的连接的唯一属性，将分析提取的文件处理成 与河道网格对应的矢量图层，最后按照分级标准，赋予不同的颜色，表征污染物影响 的变化。
(3)根据步骤3)中的一、二维模拟结果，分别绘制不同种类的突发性水污染 风险图。
①最大浓度分布图
提取设定情景模拟的整个过程各个断面的COD浓度最大值，综合考虑《地表水 环境质量标准》和《农田灌溉水质标准》GB5084-2005及模拟计算的浓度，将浓度按 照标准划分类别，分类赋予不同的颜色，表征不同的污染程度。
②不同时刻浓度分布图
绘制不同时刻的污染物浓度分布情况，便于掌握河道水质的污染变化情况。此情景下的突发性水污染事件在模拟开始后48h，河道中COD的浓度基本恢复到Ⅲ类水 标准。绘制2h、6h、12h、24h及48h的COD浓度分布图。
③污染物到达时间及影响时间分布图
根据河道污染分析结果，将到达时间及影响时间综合分级，确定污染物到达时间及影响时间分布图。
尽管为说明目的公开了本发明的实施例和附图，但是本领域的技术人员可以理解：在不脱离本发明及所附权利要求的精神和范围内，各种替换、变化和修改都是可 能的，因此，本发明的范围不局限于实施例和附图所公开的内容。
Claims (2)
1.一种河流突发性水污染分析与风险图编制方法，其特征在于：所述方法的步骤为：
1)突发性水污染情景设定：根据对突发性水污染特点及影响分析，设定两种突发性水污染事件，分别为企业污水泄露和交通运输事故，研究水域为河流干流；
2)建立水动力水质耦合模型：水质模型用于模拟污染物在水体中迁移、扩散及衰减转化过程，运用一、二维模型对步骤1)中的情景进行模拟，分析河道水污染风险；
a)一维水动力水质模型
水动力基本方程为：
水质模型基本方程为：
其中：C为整个断面的平均浓度，单位为mg/m3；
Q为流量，单位为m3/s；
A为过水断面面积，单位为m2；
u为断面平均流速，单位为m/s；
D为纵向弥散系数，单位为m2/s；
S为各种源和漏的代数和；
b)二维水动力水质模型
水流连续方程为：
水流动量方程为：
其中：t为时间；
x、y为空间坐标系；
η为水位；
h为静止水深；
u为X方向的流速；
v为Y方向的流速；
pa为当地气压；
ρ为水密度；
ρ0为参考水密度；
f为科里奥利参量；
Txx、Txy、Tyx、Tyy为水平粘滞应力；
Isx、Ibx、Ixy、Iby为水流切应力分量；
us、vs为源汇项水流流速；
水质模型基本方程为：
其中：C为污染物的浓度，单位为mg/l；
t为时间，单位为d；
ux、uy分别为x、y方向的水流速度，单位为m/s；
Dx、Dy分别为x、y方向的扩散(离散)系数，m2/s；
K为污染物的降解速率常数，单位为s-1；
3)模拟结果分析：
a)一维模拟结果分析：
①取水口污染物分析：河流干流沿岸一般分布着大量的取水口，在假定污染源位置的前提下，根据一维水动力水质耦合模型模拟各取水口的浓度变化过程线；
根据《农田灌溉水质标准》GB5084-2005规范中对农田灌溉用水水质基本控制项目浓度的规定，化学需氧量对于水作作物应小于150mg/l，旱作作物要控制在200mg/l以内。本发明约定当取水口处COD浓度首次达到150mg/l的时间为污染物到达时间，当取水口处COD浓度下降到小于150mg/l的时刻与到达时刻的时间差，即为影响时间；
②河道污染风险分析：根据《地表水环境质量标准》中对水质类别的划分，Ⅰ类、Ⅱ类水中COD的浓度小于15mg/l、Ⅲ类水COD浓度介于15mg/l-20mg/l之间，Ⅳ类水COD浓度介于20mg/l-30mg/l之间，Ⅴ类水中COD浓度介于30mg/l-40mg/l，大于40mg/l为劣Ⅴ类水；在设定突发性水污染事故的情景下，可以绘制不同位置不同时刻的COD浓度变化曲线；
b)二维模拟结果分析：在设定突发性水污染事故的情景下，二维模型考虑了横向扩散，在水动力条件下，污染团沿着主河槽向下游迁移扩散，模拟污染物平面上的浓度分布，绘制不同位置不同时刻的COD浓度变化过程平面图；
4)水污染风险图绘制：参照洪水风险图的种类，将河流水污染风险图分为污染物浓度分布图、污染物到达时间图、影响时间图等三类，并在其上叠加监测站网、污染源、取水工程信息及水源地分布等信息；将计算得到的污染物浓度分布、污染物到达时间及污染物影响时间矢量化，与其他的矢量图层叠加，最后还要对图中的符号、文字及颜色等进行设计、渲染；
(1)一维结果的矢量化：模拟计算的结果以沿程断面浓度变化的形式表征污染物的扩散及传输规律，提取不同断面的污染物最大浓度及不同时刻的浓度，并认为浓度沿河道中心线分布；经过差值加密的断面与河道中心线矢量图层叠加，并给各个交点赋予浓度值。结合浓度分级标准，以浓度值字段为对象，将浓度值进行分级，并与河道面矢量图层叠加，分级赋予不同的颜色，形成最大浓度的空间分布图及不同时刻的污染物浓度分布图；
根据模拟的结果分析各个河道断面受到污染的起始时刻，并按照一定的间隔将时间分级，在GIS中处理河道图层shp，按照分级点的控制距离将不同的时间类别赋予不同的颜色，即得到到达时间分布图层。影响时间分布层的处理方法与到达时间处理方法相同；
(2)二维结果的矢量化：二维模拟结果存储在沿河道分布的网格中，每个网格点按照时间序列存储了该点的空间位置信息和污染物的浓度值；根据网格的坐标信息，创建空间矢量图层；污染物浓度分布是空间中的连续场，为了实现污染物浓度扩散的连续二维场，在GIS中对浓度矢量图层进行IDW差值处理，按照分级标准将不同的分级赋予不同的颜色，即得到污染物浓度分布矢量图；通过C++编程，提取分析污染物的到达时间及影响历时，以网格编号作为不同文件之间转换的连接的唯一属性，将分析提取的文件处理成与河道网格对应的矢量图层，最后按照分级标准，赋予不同的颜色，表征污染物影响的变化；
(3)根据步骤3)中的一、二维模拟结果，分别进行矢量化，绘制最大浓度分布图、不同时刻浓度分布图、污染物到达时间及影响时间分布图等不同种类的突发性水污染风险图。
2.根据权利要求1所述的河流突发性水污染分析与风险图编制方法，其特征在于：所述情景设置的参数为：
(a)水文条件：考虑到不同时期上游来流量大小不同，污染物扩散和稀释的能力不同，对下游各水取水口的应急时间影响也较大；设定的两种水文条件为枯水期和丰水期的流量过程；
(b)水质条件：根据当地水环境现状的分析，确定当地水质状况，设定浓度本底值。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911125035.7A CN111027262A (zh) | 2019-11-18 | 2019-11-18 | 一种河流突发性水污染分析与风险图编制方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911125035.7A CN111027262A (zh) | 2019-11-18 | 2019-11-18 | 一种河流突发性水污染分析与风险图编制方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN111027262A true CN111027262A (zh) | 2020-04-17 |
Family
ID=70200263
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911125035.7A Pending CN111027262A (zh) | 2019-11-18 | 2019-11-18 | 一种河流突发性水污染分析与风险图编制方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111027262A (zh) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111724058A (zh) * | 2020-06-17 | 2020-09-29 | 北京林业大学 | 基于取水口水质安全的河流型饮用水源地保护系统 |
| CN112541307A (zh) * | 2020-12-22 | 2021-03-23 | 中国水利水电科学研究院 | 一种计算有毒物质泄漏位置风险程度的方法 |
| CN114648167A (zh) * | 2022-04-02 | 2022-06-21 | 重庆市生态环境科学研究院 | 一种基于山地流域污染负荷的水质预警系统及预警方法 |
| CN117094524A (zh) * | 2023-09-13 | 2023-11-21 | 北京化工大学 | 一种基于云平台的碳收集及服务系统 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105550799A (zh) * | 2015-12-07 | 2016-05-04 | 南通大学 | 流域跨界区水环境风险减免技术体系 |
-
2019
- 2019-11-18 CN CN201911125035.7A patent/CN111027262A/zh active Pending
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105550799A (zh) * | 2015-12-07 | 2016-05-04 | 南通大学 | 流域跨界区水环境风险减免技术体系 |
Non-Patent Citations (1)
| Title |
|---|
| 封桂敏: "《黄河宁夏段突发性水污染风险研究》", 《中国优秀硕士学位论文全文数据库工程科技Ⅰ辑》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111724058A (zh) * | 2020-06-17 | 2020-09-29 | 北京林业大学 | 基于取水口水质安全的河流型饮用水源地保护系统 |
| CN111724058B (zh) * | 2020-06-17 | 2023-06-16 | 北京林业大学 | 基于取水口水质安全的河流型饮用水源地保护系统 |
| CN112541307A (zh) * | 2020-12-22 | 2021-03-23 | 中国水利水电科学研究院 | 一种计算有毒物质泄漏位置风险程度的方法 |
| CN114648167A (zh) * | 2022-04-02 | 2022-06-21 | 重庆市生态环境科学研究院 | 一种基于山地流域污染负荷的水质预警系统及预警方法 |
| CN114648167B (zh) * | 2022-04-02 | 2024-01-23 | 重庆市生态环境科学研究院 | 一种基于山地流域污染负荷的水质预警系统及预警方法 |
| CN117094524A (zh) * | 2023-09-13 | 2023-11-21 | 北京化工大学 | 一种基于云平台的碳收集及服务系统 |
| CN117094524B (zh) * | 2023-09-13 | 2024-03-19 | 北京化工大学 | 一种基于云平台的碳收集及服务系统 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111027262A (zh) | 一种河流突发性水污染分析与风险图编制方法 | |
| Rui et al. | GIS-based emergency response system for sudden water pollution accidents | |
| He et al. | Development and application of a water pollution emergency response system for the Three Gorges Reservoir in the Yangtze River, China | |
| CN107798467B (zh) | 基于深度学习的水污染突发事故快速应急评估和决策方法 | |
| CN104392100A (zh) | 基于水质在线监测系统的污染源扩散预警方法 | |
| CN109886568A (zh) | 一种危化品道路运输流域水环境风险评估方法及系统 | |
| CN115204060B (zh) | 一种饮用水源地水质多维模型快速预警方法及系统 | |
| CN103489036B (zh) | 基于混合空间关联分析的应急预警网格划分系统及方法 | |
| CN109740971A (zh) | 一种危化品生产企业突发流域水环境的风险评估方法及系统 | |
| Yerramilli et al. | Air pollution, modeling and GIS based decision support systems for air quality risk assessment | |
| Liu et al. | Research on transport and weathering of oil spills in Jiaozhou Bight, China | |
| Cheng et al. | Unloading characteristics of sand-drift in wind-shallow areas along railway and the effect of sand removal by force of wind | |
| CN109816264A (zh) | 一种区域突发流域水环境的风险评估方法及系统 | |
| Zhai et al. | Hazardous chemical accident prediction for drinking water sources in Three Gorges Reservoir | |
| Liu et al. | Prediction and countermeasures of heavy metal copper pollution accident in the Three Gorges Reservoir Area | |
| Sivanandan et al. | A study on plume dispersion characteristics of two discrete plume stacks for negative temperature gradient conditions | |
| Xu et al. | CFD simulation of chemical gas dispersion under atmospheric boundary conditions | |
| Issakhov et al. | Numerical study of the dynamics of particles motion with different sizes from coal-based thermal power plant | |
| Oaks et al. | Lagrangian dynamics of contaminant particles released from a point source in New York City | |
| Zeng et al. | A two-dimensional water-quality model for a winding and topographically complicated river | |
| Li et al. | A real-time GIS platform for high sour gas leakage simulation, evaluation and visualization | |
| Mansor | Computational Modelling of Open Channel Flow Around Bridge Piers | |
| Cui et al. | System development for storm surge hazard assessment based on WebGIS for Tianjin Binhai New Area | |
| Chin | Outfall dilution: the role of a far-field model | |
| CN114037317A (zh) | 一种基于决策树模型的地下管网运行能力评估方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200417 |
|
| RJ01 | Rejection of invention patent application after publication |