长江干流铜陵段突发水污染情景模拟分析

齐于顺, 刘仁志, 张启月, 崔保山, 郭忠, 王飞, 方丽, 方陵

齐于顺,刘仁志,张启月,等.长江干流铜陵段突发水污染情景模拟分析[J].环境工程技术学报,2022,12(2):607-614. DOI: 10.12153/j.issn.1674-991X.20210662
引用本文: 齐于顺,刘仁志,张启月,等.长江干流铜陵段突发水污染情景模拟分析[J].环境工程技术学报,2022,12(2):607-614. DOI: 10.12153/j.issn.1674-991X.20210662
QI Y S,LIU R Z,ZHANG Q Y,et al.Scenario simulation analysis of abrupt water pollution in Tongling section of the mainstream of the Yangtze River[J].Journal of Environmental Engineering Technology,2022,12(2):607-614. DOI: 10.12153/j.issn.1674-991X.20210662
Citation: QI Y S,LIU R Z,ZHANG Q Y,et al.Scenario simulation analysis of abrupt water pollution in Tongling section of the mainstream of the Yangtze River[J].Journal of Environmental Engineering Technology,2022,12(2):607-614. DOI: 10.12153/j.issn.1674-991X.20210662

长江干流铜陵段突发水污染情景模拟分析

基金项目: 铜陵市长江生态环境保护修复驻点跟踪研究项目(2020CGSF134);国家自然科学基金项目(52170186)
详细信息
    作者简介:

    齐于顺(1998—),男,硕士研究生,主要从事环境风险管理研究,202021180054@mail.bnu.edu.cn

    通讯作者:

    刘仁志(1974—),男,副教授,主要从事区域环境承载力、环境风险、环境评价规划与管理研究, liurenzhi@bnu.edu.cn

  • 中图分类号: X522

Scenario simulation analysis of abrupt water pollution in Tongling section of the mainstream of the Yangtze River

  • 摘要: 突发水环境事件情景模拟分析对风险防控与应急处置至关重要。在对铜陵市沿江化工企业风险物质泄漏情景分析的基础上,采用MIKE21二维水动力-水质耦合模型对长江干流铜陵段突发水污染事故进行了模拟。模型验证结果表明,构建的水动力模型满足精度要求。对某化工企业储罐3种泄漏情景下的模拟结果表明:同一水文期,风险物质对下游敏感受体的污染程度和污染时间主要与泄漏总量有关,储罐完全泄漏情景对下游敏感受体影响最大,其最大污染峰团浓度高于罐体20%管径破裂和罐体100%管径破裂泄漏情景2~3个数量级;不同水文期,泄漏物质在丰水期到达下游敏感受体的时间最短,储罐完全泄漏时风险物质到达下游三水厂取水口用时75 min,而在平水期和枯水期分别为103和111 min,同时影响时间更长、浓度更高,下游五水厂取水口和铜陵市出境断面在平水期和枯水期先后有2次污染峰团到达。
    Abstract: Scenario simulation analysis of abrupt water pollution is important for risk prevention and control and emergency response. Based on the analysis of risky substance leakage scenarios of chemical enterprises, MIKE21, a two-dimensional hydrodynamic-water quality coupling model, was used to simulate the abrupt water pollution accident in the Tongling section of the mainstream of the Yangtze River. The model verification showed the constructed hydrodynamic model could meet the accuracy requirements. The simulation results of three different leakage scenarios for one chemical company's storage tank showed that: in the same hydrological period, the pollution degree and time of risky substances to downstream sensitive receptors were mainly related to the total amount of leakage, and it had the most serious impact on downstream sensitive receptors when the storage tank was leaked entirely, and the maximum pollution peak was 2~3 orders of magnitude higher than 20% and 100% pipe diameter leakage. In different hydrological periods, risk substances took the shortest time to get to downstream sensitive receptors in the wet period, and it took 75 min for the risky material to reach the water intake of the Third Water Plant when the storage tank was leaked completely; but it took 103 min and 111 min during the normal and dry periods, respectively, and the impact time was longer and the concentration was higher. And there were two pollution peaks successively at the water intake of the downstream plant and the outbound section of the city in the normal and dry periods.
  • 突发环境事件是指由于污染物排放或自然灾害、生产安全事故等因素,导致有毒有害物质大量进入环境介质,突然造成或可能造成环境质量下降,危及公众身体健康和财产安全,需要采取紧急措施予以应对的事件[1]。与其他普通污染事件不同的是,突发环境事件具有随机性、复杂性、高强度、高风险的特点[2],会在短时间内造成周边环境污染物超标,甚至会造成长期影响[3]。其中,突发水环境事件最易发生,李旭等[4]对2011—2017年突发环境污染事件的受体类型进行统计发现,水是主要污染受体,占比高达85.59%。长江沿岸城市“化工围江”、上下游排污口与取水口交错分布等布局性风险突出,导致长江水环境风险形势不容乐观[5-6]。此外,由于长江水体流动性强的特点,突发水环境事件的影响范围更大,响应时间更短,如果不能及时妥善处置,短时间内可能导致饮用水水源地污染、水域生态环境破坏,甚至造成跨界污染,产生极大的不良社会影响[7-8]。目前突发水环境事件仍是影响供水和人民安全的重要威胁因素[9-10]。因此,有必要通过多种方式提前做好突发水环境事件的预测预警与应急准备,有效防范突发水污染事故,降低突发水环境风险[11]

    水环境模型可对水环境系统及其内部发生的复杂过程进行定量化描述[12],模拟污染物在一定时间和空间范围内的迁移转化过程,计算污染物到达下游环境敏感目标的时间、污染团的历经时间和污染物的浓度分布等[13],有助于帮助应急管理人员了解事故可能的发展趋势并做出合理的处置。常用的河流水质模型主要有EFDC、WASP、QUAL2E、QUAL2K、MIKE、BASINS、SWAT等[13-17]。其中,丹麦水力研究所开发的平面二维模型MIKE21,因用户界面友好,具有多种计算网格、模块供用户选择,且有强大的前、后处理功能,被广泛用于河流[18-19]、湖泊[20-21]、水库[22]、海湾[23]等水质模拟。如袁玥等[18]应用MIKE21模型对长江蕲春段非正常排污情况进行了模拟,舒长莉等[24]应用MIKE21模型对赣江南昌段突发污染事故进行了模拟,均实现了对污染影响后果的量化表达。为贯彻落实习近平总书记“不搞大开发,共抓大保护”的重要指示精神,着力解决长江水环境、水生态、水风险突出问题[25],为长江经济带高质量发展提供保证[26],笔者以长江生态环境保护修复驻点跟踪研究工作为基础,利用MIKE21构建长江干流铜陵段突发水污染事故水动力-水质耦合模型,考虑了多种可能的泄漏情景,并对不同水文期进行模拟分析,实现对突发水污染事故的动态模拟,以期为长江干流铜陵段水域突发水污染事故风险防范、预测预警和应急处置提供参考。

    铜陵市位于安徽省中南部、长江下游(117°04’E~118°09’E,30°38’N~31°09’N),是长江经济带重要节点城市。长江干流铜陵段全长142.6 km,过境的长江水资源丰富,是铜陵市主要供水来源,最小日平均流量为10 700 m3/s,最大日平均流量为68 400 m3/s,年平均流量为29 250 m3/s(2019年),水量年内变化显著。研究河段内环境敏感目标包括铜陵市淡水豚国家级自然保护区和3处集中式生活饮用水水源地(铜陵市水厂、新三水厂和第五水厂)。淡水豚国家级自然保护区范围上始铜陵市枞阳县老洲镇,下至铜陵市义安区金牛渡,全长58 km,其中核心区、缓冲区和实验区面积分别为9 534、6 360、15 624 hm2[27]。横港扫把沟地区紧临长江,是铜陵市最早的工业聚集区之一,由于历史原因,该区域功能布局不合理,工业用地与居住用地混杂,多家企业距长江直线距离不足1 km。目前该区域仍存留多家风险企业,也没有成立工业园区管委会,这给下游集中式饮用水水源地、淡水豚国家级自然保护区以及长江水体水质安全带来重大威胁。研究区概况如图1所示。

    图  1  研究区概况
    Figure  1.  Schematic map of the study area

    研究区水文数据采用2019年大通水文站实测逐日流量数据和横港、荻港2个水位站实测逐日水位数据;地形文件为1∶25 000实测水下地形数据;风险企业、饮用水水源保护区和淡水豚保护区等相关资料由铜陵市生态环境局提供。

    长江铜陵段江面宽度远大于深度,且不存在明显分层现象,因此采用二维水动力模型可以满足研究需要[28]。选取MIKE21 Flow Model FM模型进行模拟,二维水动力模型基于三向不可压缩和Reynolds值均布的Navier-Stokes方程,并服从于Boussinesq假定和静水压力的假定,水动力控制方程组(浅水方程)为[24,29-31]

    $$ \frac{\partial h}{\partial t}+\frac{\partial h{\stackrel-{u}}}{\partial x}+\frac{\partial h{\stackrel-{v}}}{\partial y}=hS $$ (1)
    $$ \begin{aligned} &\frac{\partial h{\stackrel-{u}}}{\partial t}+\frac{\partial h{\stackrel-{u}}^{2}}{\partial x}+\frac{\partial h{\stackrel-{uv}}}{\partial y}=\\ &f{\stackrel-{v}}h-gh\frac{\partial \eta }{\partial x}-\frac{h}{{\rho }_{0}}\frac{\partial {P}_{{\rm{a}}}}{\partial x}-\frac{g{h}^{2}}{2{\rho }_{0}}\frac{\partial \rho }{\partial x}+\frac{{\tau }_{sx}}{{\rho }_{0}}-\frac{{\tau }_{bx}}{{\rho }_{0}}-\\ &\frac{1}{{\rho }_{0}}\left(\frac{{\partial s}_{xx}}{\partial x}+\frac{{\partial s}_{xy}}{\partial y}\right)+\frac{\partial }{\partial x}\left(h{T}_{xx}\right)+\frac{\partial }{\partial y}\left(h{T}_{xy}\right)+h{u}_{s}S \end{aligned} $$ (2)
    $$ \begin{aligned} &\\ &\frac{\partial h\stackrel-{v}}{\partial t}+\frac{\partial h{\stackrel-{v}}^{2}}{\partial y}+\frac{\partial h\stackrel-{uv}}{\partial x}=\\ &-f\stackrel-{u}h-gh\frac{\partial \eta }{\partial y}-\frac{h}{{\rho }_{0}}\frac{\partial {P}_{{\rm{a}}}}{\partial y}-\frac{g{h}^{2}}{2{\rho }_{0}}\frac{\partial \rho }{\partial y}+\frac{{\tau }_{sy}}{{\rho }_{0}}-\frac{{\tau }_{by}}{{\rho }_{0}}-\\ &\frac{1}{{\rho }_{0}}\left(\frac{{\partial s}_{yx}}{\partial x}+\frac{{\partial s}_{yy}}{\partial y}\right)+\frac{\partial }{\partial x}\left(h{T}_{yx}\right)+\frac{\partial }{\partial y}\left(h{T}_{yy}\right)+h{v}_{s}S \end{aligned} $$ (3)

    式中:t为时间;xy为笛卡尔坐标;$ \eta $为水位,m;d为静止水深,m;$ h $为总水深$ ,h=d+\eta $,m;uv分别为xy方向上的速度分量;f为哥氏力系数,$f= $$ 2\omega {\rm{sin}}\varphi$,其中$ \omega $为地球自转角速度,$ \varphi $为当地纬度;g为重力加速度,取9.81 ${\rm{m}}/{{\rm{s}}}^{2}$$ \rho $$ {\rho }_{0} $分别为水体和参考水的密度,${\rm{kg}}/{{\rm{m}}}^{3}$Pa为大气压力,Pa;$ {s}_{xx} $$ {s}_{xy} $$ {s}_{yx} $$ {s}_{yy} $分别为x轴的法线方向、y轴的法线方向、$ xy $表面上的切线方向和$ yx $表面上的切线方向的辐射应力分量;$ {A}_{x} $$ {A}_{y} $为应力分项;$ {T}_{xx} $$ {T}_{xy} $$ {T}_{yx} $$ {T}_{yy} $为黏滞切应力分量;S为源汇项;$ {u}_{s} $$ {v}_{s} $为源项水流流速,m/s;$ \bar{u} $$ \bar{v} $为沿水深平均的流速;τsxτsy为水面风应力在xy 方向上的分量; τbxτby为河床( 海底) 底部应力的分量。

    $ \bar{u} $$ \bar{v} $计算公式如下:

    $$ h\stackrel-{u}={\int }_{-d}^{n}udz,h\stackrel-{v}={\int }_{-d}^{n}vdz $$ (4)

    水质模块建立于二维空间下的输移扩散方程为:

    $$ \frac{\partial C}{\partial t}+u\frac{\partial C}{\partial x}+v\frac{\partial C}{\partial y}={D}_{x}\frac{{\partial }^{2}C}{\partial {x}^{2}}+{D}_{y}\frac{{\partial }^{2}C}{\partial {y}^{2}}-P+S $$ (5)

    式中:C为各典型污染物浓度;$ {D}_{x} $$ {D}_{y} $为各典型污染物在xy上的扩散系数;P为各典型污染物降解项。

    MIKE21模型采用有限体积法对计算区域进行空间离散,将该连续统一体细分为若干个不重叠的三角形或四边形单元。模型计算的时间和精度取决于计算数值方法所使用的求解格式精度,浅水方程的时间积分和输移扩散方程基于半隐格式求解,相应平流项采用显式格式求解,垂直对流项采用全隐格式求解。受显式格式稳定性的限制,为保持模型计算的稳定性,模型中时间步长的设定必须保证CFL(Courant-Friedrich Levy)数小于1。浅水方程和输移扩散方程在笛卡尔坐标上的CFL分别定义为:

    $$ {\mathrm{C}\mathrm{F}\mathrm{L}}_{{\rm{HD}}}=\left(\sqrt{gh}+\left|u\right|\right)\frac{\Delta t}{\Delta x}+\left(\sqrt{gh}+\left|v\right|\right)\frac{\Delta t}{\Delta y} $$ (6)
    $$ {\mathrm{C}\mathrm{F}\mathrm{L}}_{\mathrm{A}\mathrm{D}}=\left|u\right|\frac{\Delta t}{\Delta x}+\left|v\right|\frac{\Delta t}{\Delta y} $$ (7)

    式中:$ \Delta x $$ \Delta y $xy方向的特征长度;$ \Delta t $为时间间距。$ \Delta x $$ \Delta y $近似于三角形网格的最小边长,水深和流速值为三角形网格中心的取值。

    利用MIKE21建立水动力-水质耦合模型,模型构建流程如图2所示。基于铜陵市水利工程图和实测水下地形数据,利用ArcMap与AutoCAD软件提取生成长江铜陵段水陆边界线(Land. xyz文件),并读取相应水深散点数据(Water. xyz文件)。在Mesh Generator 中导入水陆边界线(Land. xyz文件)并进行光滑处理,采用非结构三角网格对河道地形进行处理,然后生成网格;插入地形水深散点数据(Water. xyz文件)并进行地形插值后导出mesh文件(图3),计算区域共有10 896个三角网格,6 886个网格节点。模型上边界采用2019年大通水文站逐日流量时间序列,下边界采用2019年荻港水位站逐日水位时间序列。选择Flow Model FM建立水动力-水质耦合模型,经过多次调试后选择合适参数;将模拟时长设为24~30 h,时间步长选择60 s;为保证模型稳定运行,设置CFL为0.8,最大时间步长为60 s;干水深$({h}_{{\rm{dry}}})为 0.005\;{\rm{m}}$,淹没水深$({h}_{{\rm{flood}}})为0.05\;{\rm{m}}$,湿水深${(h}_{{\rm{wet}}})为 $$ 0.1\;{\rm{m}}$;底床摩擦力采用曼宁系数,取值32 m1/3/s;初始水位根据不同水文期选择模拟对应开始日期水位值;其他参数采用推荐或默认值。在对流扩散模块,添加需要模拟的风险物质硫酸,从更安全的角度考虑预测结果,将硫酸降解系数设为0[32];并在源汇项中添加风险物质初始浓度和泄漏流量等相关数据。

    图  2  模型构建流程
    Figure  2.  Model construction flowchart
    图  3  水动力模型网格地形高程
    Figure  3.  Grid terrain elevation of hydrodynamic model

    因长江铜陵段水文情势变化显著,有明显的丰、平、枯水文情势节律性变化,且枯水期部分时间段水位站处于不工作状态。故对丰、平、枯水文时期部分时间段进行模拟,并选取2019年横港水位站逐日水位数据进行验证,将实测水位数据与模拟水位数据进行对比。经对照发现(图4),模拟水位与实际水位数值相差不大,平均相对误差为2.4%,曲线趋势基本相似,波峰与波峰相对,波谷与波谷相对。模拟值与实测值吻合性良好,可作为水质模块的水动力基础。

    图  4  横港水位站2019年模拟水位与实测水位
    Figure  4.  Simulated and measured water level of Henggang water-level station in 2019

    横港扫把沟地区多家风险企业的储罐内存在大量风险物质,给下游集中式饮用水水源地、淡水豚国家级自然保护区以及长江水体水质带来巨大安全隐患。针对某化工企业硫酸储罐,设定3种泄漏情景,分别为罐体20%管径破裂、罐体100%管径破裂和储罐完全泄漏,泄漏时间均设为10 min。在丰水期、平水期和枯水期对3种泄漏情景分别进行模拟,具体泄漏时间分别设定为7月10日、4月16日和10月12日的06:00。所设定的泄漏情景如表1所示。

    表  1  3种设定泄漏情景
    Table  1.  Three supposed leakage situations
    水文期泄漏情景泄漏时间/min
    丰水期/平水期/枯水期情景A:罐体20%管径破裂10
    情景B:罐体100%管径破裂10
    情景C:储罐完全泄漏10
    下载: 导出CSV 
    | 显示表格

    硫酸罐体20%管径和100%管径破裂时泄漏速率用伯努利方程[33]计算(液体在泄漏口没有急骤蒸发),公式如下:

    $$ {Q}_{\mathrm{L}}={C}_{{\rm{d}}}A\rho \sqrt{\frac{2(P_{1}-{P}_{0})}{\rho }+2gh} $$ (8)

    式中:$ {Q}_{\mathrm{L}} $为液体泄漏速率,$ \mathrm{k}\mathrm{g}/\mathrm{s} $$P_{1}$为容器内介质压力,Pa;$ {P}_{0} $为环境压力,Pa;$ \rho $为泄漏液体密度,$ \mathrm{k}\mathrm{g}/{\mathrm{m}}^{3} $$ h $为裂口之上液位高度,$ \mathrm{m} $$ {C}_{\mathrm{d}} $为液体泄漏系数,取0.62;$ A $为裂口面积,$ {\mathrm{m}}^{2} $

    根据实地调研,判断突发泄漏情况下的可能入江位置如图1所示,将不同泄漏情景下计算得到的泄漏流量-时间序列等相关数据输入构建的水动力-水质耦合模型,可得到不同水文期的风险物质迁移扩散分布结果。

    采用MIKE21水动力-水质耦合模型对3种泄漏情景在丰水期、平水期和枯水期的风险物质迁移扩散情况分别进行了模拟,定量化表达出风险物质到达下游环境敏感目标的时间、最大污染峰团浓度和污染持续时间等。泄漏硫酸对长江水质的影响可通过对水体pH的影响来反映[32],GB 3838—2002《地表水环境质量标准》规定地表水Ⅰ类~Ⅴ类水域水体pH为6~9,当pH为6时,通过计算得到对应硫酸的浓度为0.049 mg/L,即硫酸浓度大于0.049 mg/L时会对下游各风险受体和长江水质造成影响。通过水质模型模拟结果可得到储罐泄漏后各时间节点的浓度,进而预测和判断影响程度。下面针对同一水文期不同泄漏情景和同一泄漏情景在不同水文期的水质模拟结果进行说明。

    以3种不同泄漏情景在丰水期模拟结果为例进行说明,模拟结果如图5~图7所示。3种泄漏情景下最大污染峰团到达下游敏感受体所需时间相同,到达三水厂取水口、市水厂取水口、淡水豚核心区、五水厂取水口和市出境断面的时间分别为89、93、332、677、751 min,这是由于同一时间段的水动力结果相同。但在储罐完全泄漏情景下,风险物质到达下游敏感受体时间稍早于其他2种泄漏情景,风险物质完全流过敏感受体时间稍晚于其他2种泄漏情景,即影响时间更长。此外,该情景下最大污染峰团浓度也明显高于其他2种泄漏情景下2~3个数量级。在三水厂取水口处,20%管径破裂时最大污染峰团浓度为1.891 mg/L,100%管径破裂时最大污染峰团浓度为47.273 mg/L,储罐完全泄漏时最大污染峰团浓度为1 381.39 mg/L,这是由于不同泄漏情景下风险物质泄漏总量不同。随着时间的延长,风险物质被稀释,到达下游各敏感受体时的最大污染峰团浓度会逐渐降低。在3种不同泄漏情景下,风险物质到达市出境断面时的最大污染峰团浓度分别为0.090 5、2.262 4和66.123 4 mg/L,均超过0.049 mg/L的标准要求,会造成不同程度的跨界污染。淡水豚核心保护区范围最大,其受影响持续时间也最长,分别为251、313、343 min。

    图  5  丰水期罐体20%管径破裂时水质模拟结果
    Figure  5.  Simulation results of water quality when 20% pipe diameter of tank breaks in wet period
    图  6  丰水期罐体100%管径破裂时水质模拟结果
    Figure  6.  Simulation results of water quality when 100% pipe diameter of tank breaks in wet period
    图  7  丰水期储罐完全泄露时水质模拟结果
    Figure  7.  Water quality simulation results when the tank breaks completely in wet period
    图  8  平水期储罐完全泄露时水质模拟结果
    Figure  8.  Water quality simulation results when the tank breaks completely in normal period
    图  9  枯水期储罐完全泄露时水质模拟结果
    Figure  9.  Water quality simulation results when the tank breaks completely in dry period

    以不同水文期储罐完全泄漏情景的模拟结果为例进行说明,模拟结果如图7~图9所示。根据模拟结果可知,丰水期风险物质到达下游敏感受体时间最短,其中风险物质到达三水厂取水口仅用75 min,影响持续时间为43 min。相较于丰水期,平水期和枯水期的流速更小,风险物质需更长时间到达下游敏感受体,如平水期、枯水期风险物质到达三水厂取水口的时间分别为103和111 min。但在平水期和枯水期,风险物质对敏感受体的影响持续时间会更长,其对三水厂取水口的影响持续时间分别为65和72 min。同时因流量更小,风险物质在市水厂取水口、三水厂取水口和淡水豚核心区的最大污染峰团浓度会更高,枯水期的浓度分别为1 883、1 736和426 mg/L。而到达下游更远的五水厂取水口时,丰水期、平水期、枯水期浓度依次降低,其主要原因可能是停留时间较长,风险物质被稀释[24]。平水期和枯水期,在五水厂取水口、市出境断面处会先后出现2次污染峰团,这是由于江心洲两侧水体流速相差较大,部分风险物质随右侧支流先行到达敏感受体,随后左侧支流中的风险物质再次到达敏感受体。

    综上,风险物质对下游风险受体的影响程度与泄漏总量和水文期密切相关,所得结论与长江其他段类似研究[29,31]相同。同一水文期,风险物质泄漏总量越大,对下游敏感受体的影响时间越长,污染团浓度也越高。泄漏风险物质在平水期和枯水期时,对受体的影响时间更长;在丰水期时,到达受体的时间更短,即应急反应时间更少。

    (1)利用MIKE21水动力-水质耦合模型,对长江铜陵段某化工企业储罐内风险物质在不同水文期、3种不同泄漏情景下的突发水污染事故进行了模拟,所建立的水动力模型模拟水位结果与实测结果吻合性良好,曲线趋势基本相似,表明模型选取参数较为合理,可作为水质模块的水动力基础。

    (2)在设定的3种不同泄漏情景下,储罐完全泄漏时因其泄漏总量大,导致风险物质浓度最高,影响时间最长。在不同水文期的模拟结果显示,丰水期风险物质到达下游敏感受体的时间最短,留给应急人员处置时间也最少;平水期和枯水期风险物质的最大污染峰团浓度会更高,影响时间会更长,且在五水厂取水口和市出境断面处先后有2次污染峰团到达。

    (3)本研究基于MIKE21构建了适用于长江干流铜陵段的水动力-水质模型,通过对风险企业内多种风险物质在不同水文情势和泄漏情景下的模拟分析,可建立当地的模拟预测预警资料库,为突发水污染事故的预测预警、应急处置提供决策参考。结合本次模拟结果,因风险物质在丰水期到达受体时间更短,平水期和枯水期对受体影响时间更长,当接到突发水环境事件报告时,在丰水期应立即采取措施保护风险受体,在平水期和枯水期则首先对泄漏风险源进行控制,从而减少影响时长。此外,有必要通过加强日常环境安全检查等措施,及时发现并整改问题,从而有效降低对水环境的影响程度。后续研究还需采集更为丰富的水文数据、水质监测数据等,不断进行完善以进一步提高模型模拟的准确性。

  • 图  1   研究区概况

    Figure  1.   Schematic map of the study area

    图  2   模型构建流程

    Figure  2.   Model construction flowchart

    图  3   水动力模型网格地形高程

    Figure  3.   Grid terrain elevation of hydrodynamic model

    图  4   横港水位站2019年模拟水位与实测水位

    Figure  4.   Simulated and measured water level of Henggang water-level station in 2019

    图  5   丰水期罐体20%管径破裂时水质模拟结果

    Figure  5.   Simulation results of water quality when 20% pipe diameter of tank breaks in wet period

    图  6   丰水期罐体100%管径破裂时水质模拟结果

    Figure  6.   Simulation results of water quality when 100% pipe diameter of tank breaks in wet period

    图  7   丰水期储罐完全泄露时水质模拟结果

    Figure  7.   Water quality simulation results when the tank breaks completely in wet period

    图  8   平水期储罐完全泄露时水质模拟结果

    Figure  8.   Water quality simulation results when the tank breaks completely in normal period

    图  9   枯水期储罐完全泄露时水质模拟结果

    Figure  9.   Water quality simulation results when the tank breaks completely in dry period

    表  1   3种设定泄漏情景

    Table  1   Three supposed leakage situations

    水文期泄漏情景泄漏时间/min
    丰水期/平水期/枯水期情景A:罐体20%管径破裂10
    情景B:罐体100%管径破裂10
    情景C:储罐完全泄漏10
    下载: 导出CSV
  • [1] 国家突发环境事件应急预案[EB/OL]. (2019-12-27)[2021-10-17]. https://www.mee.gov.cn/ywgz/hjyj/yjzb/201912/t20191227_751708.shtml.
    [2] 宋永会, 袁鹏, 彭剑峰, 等.突发环境事件风险源识别与监控技术创新进展: (Ⅰ)环境风险源识别技术与应用[J]. 环境工程技术学报,2015,5(5):347-352. DOI: 10.3969/j.issn.1674-991X.2015.05.055

    SONG Y H, YUAN P, PENG J F, et al. Technological innovation progress of risk sources identification and monitoring of sudden environmental pollution accidents: (Ⅰ) risk sources identification technologies and applications[J]. Journal of Environmental Engineering Technology,2015,5(5):347-352. DOI: 10.3969/j.issn.1674-991X.2015.05.055

    [3] 陈成健.突发环境事件应急与管理研究[J]. 资源节约与环保,2018(7):122. DOI: 10.3969/j.issn.1673-2251.2018.07.098
    [4] 李旭, 吕佳佩, 裴莹莹, 等.国内突发环境事件特征分析[J]. 环境工程技术学报,2021,11(2):401-408.

    LI X, LÜ J P, PEI Y Y, et al. Analysis of the characteristics of environmental emergencies in China[J]. Journal of Environmental Engineering Technology,2021,11(2):401-408.

    [5] 李海生, 王丽婧, 张泽乾, 等.长江生态环境协同治理的理论思考与实践[J]. 环境工程技术学报,2021,11(3):409-417. DOI: 10.12153/j.issn.1674-991X.20210071

    LI H S, WANG L J, ZHANG Z Q, et al. Theoretical thought and practice of eco-environment synergistic management in the Yangtze River[J]. Journal of Environmental Engineering Technology,2021,11(3):409-417. DOI: 10.12153/j.issn.1674-991X.20210071

    [6] 刘录三, 黄国鲜, 王璠, 等.长江流域水生态环境安全主要问题、形势与对策[J]. 环境科学研究,2020,33(5):1081-1090.

    LIU L S, HUANG G X, WANG F, et al. Main problems, situation and countermeasures of water eco-environment security in the Yangtze River Basin[J]. Research of Environmental Sciences,2020,33(5):1081-1090.

    [7] 陈思莉, 黄大伟, 张政科, 等.流域突发水污染事件应急处置工程削污技术[J]. 环境工程学报,2021,15(7):2233-2238. DOI: 10.12030/j.cjee.202008108

    CHEN S L, HUANG D W, ZHANG Z K, et al. Pollution reduction technologies for emergent water pollution disposal in river basins[J]. Chinese Journal of Environmental Engineering,2021,15(7):2233-2238. DOI: 10.12030/j.cjee.202008108

    [8] 李春晖, 田雨桐, 赵彦伟, 等.突发水污染风险评价与应急对策研究进展[J]. 农业环境科学学报,2020,39(6):1161-1167. DOI: 10.11654/jaes.2020-0411

    LI C H, TIAN Y T, ZHAO Y W, et al. Research progress on risk assessment and emergency countermeasures of sudden water pollution[J]. Journal of Agro-Environment Science,2020,39(6):1161-1167. DOI: 10.11654/jaes.2020-0411

    [9]

    PAN F H, WANG Y B, ZHANG X X. Emergency measure of soft isolation controlling pollution diffusion response to sudden water pollution accidents[J]. Water Science and Technology:A Journal of the International Association on Water Pollution Research,2019,80(7):1238-1248. DOI: 10.2166/wst.2019.368

    [10]

    CHENG T J, WANG P P, LU Q Y. Risk scenario prediction for sudden water pollution accidents based on Bayesian networks[J]. International Journal of System Assurance Engineering and Management,2018,9(5):1165-1177.

    [11] 张珂, 刘仁志, 张志娇, 等.流域突发性水污染事故风险评价方法及其应用[J]. 应用基础与工程科学学报,2014,22(4):675-684.

    ZHANG K, LIU R Z, ZHANG Z J, et al. A method of environmental risk assessment for abrupt water pollution accidents in river basin[J]. Journal of Basic Science and Engineering,2014,22(4):675-684.

    [12] 朱瑶, 梁志伟, 李伟, 等.流域水环境污染模型及其应用研究综述[J]. 应用生态学报,2013,24(10):3012-3018.

    ZHU Y, LIANG Z W, LI W, et al. Watershed water environment pollution models and their applications: a review[J]. Chinese Journal of Applied Ecology,2013,24(10):3012-3018.

    [13] 李林子, 钱瑜, 张玉超.基于EFDC和WASP模型的突发水污染事故影响的预测预警[J]. 长江流域资源与环境,2011,20(8):1010-1016.

    LI L Z, QIAN Y, ZHANG Y C. Forecasting and warning the accidental water pollution effect based on the EFDC and WASP[J]. Resources and Environment in the Yangtze Basin,2011,20(8):1010-1016.

    [14]

    ZHENG H Z, LEI X H, SHANG Y Z, et al. Sudden water pollution accidents and reservoir emergency operations: impact analysis at Danjiangkou Reservoir[J]. Environmental Technology,2018,39(6):787-803. DOI: 10.1080/09593330.2017.1311945

    [15] 周华.河流综合水质模型QUAL2K应用研究[J]. 中国水利水电科学研究院学报,2010,8(1):71-75. DOI: 10.3969/j.issn.1672-3031.2010.01.012

    ZHOU H. Application research of a comprehensive river water quality model QUAL2K[J]. Journal of China Institute of Water Resources and Hydropower Research,2010,8(1):71-75. DOI: 10.3969/j.issn.1672-3031.2010.01.012

    [16] 郭晓明, 乔明.MIKE21模型在排污口地表水影响预测方面的应用[J]. 四川环境,2021,40(3):117-123.

    GUO X M, QIAO M. Application of MIKE21 model in prediction of environmental impact of sewage outlet on surface water[J]. Sichuan Environment,2021,40(3):117-123.

    [17] 王晓青, 李哲.SWAT与MIKE21耦合模型及其在澎溪河流域的应用[J]. 长江流域资源与环境,2015,24(3):426-432. DOI: 10.11870/cjlyzyyhj201503011

    WANG X Q, LI Z. SWAT and MIKE21 coupled models and their application in the Pengxi Watershed[J]. Resources and Environment in the Yangtze Basin,2015,24(3):426-432. DOI: 10.11870/cjlyzyyhj201503011

    [18] 袁玥, 师懿, 程胜高.基于二维动态水质模型的长江蕲春段水污染扩散分析[J]. 安全与环境学报,2016,16(1):212-217.

    YUAN Y, SHI Y, CHENG S G. Analysis of water pollution dispersion at Qichun section of the Yangtze River based on the 2-D dynamic water quality model[J]. Journal of Safety and Environment,2016,16(1):212-217.

    [19]

    WU D A, SHAO Y C, PAN J X. Study on activities and concentration of saline groupin the south branch in Yangtze River Estuary[J]. Procedia Engineering,2015,116:1085-1094. DOI: 10.1016/j.proeng.2015.08.349

    [20] 宫雪亮, 孙蓉, 芦昌兴, 等.基于MIKE21的南四湖上级湖水量水质响应模拟研究[J]. 中国农村水利水电,2019(1):70-76. DOI: 10.3969/j.issn.1007-2284.2019.01.014

    GONG X L, SUN R, LU C X, et al. Research on the water quality and quantity response based on MIKE21 in the upper-reach lakes of Nansi Lake[J]. China Rural Water and Hydropower,2019(1):70-76. DOI: 10.3969/j.issn.1007-2284.2019.01.014

    [21]

    ZHU C J, LIANG Q, YAN F, et al. Reduction of waste water in Erhai Lake based on MIKE21 hydrodynamic and water quality model[J]. The Scientific World Journal,2013,2013:958506.

    [22]

    XU M J, YU L, ZHAO Y W, et al. The simulation of shallow reservoir eutrophication based on MIKE21: a case study of Douhe Reservoir in North China[J]. Procedia Environmental Sciences,2012,13:1975-1988. DOI: 10.1016/j.proenv.2012.01.191

    [23] 常赜, 张琼海, 姜宇, 等.基于MIKE21模型的入海排污口对澳门近岸海域水质影响分析[J]. 人民珠江,2021,42(3):70-80. DOI: 10.3969/j.issn.1001-9235.2021.03.010

    CHANG Z, ZHANG Q H, JIANG Y, et al. Analysis of the impact of sewage outfall on the water quality of coastal waters of Macao based on MIKE21 model[J]. Pearl River,2021,42(3):70-80. DOI: 10.3969/j.issn.1001-9235.2021.03.010

    [24] 舒长莉, 李林, 冯韬.基于MIKE21的河道饮用水源地突发污染事故模拟: 以赣江南昌段为例[J]. 人民长江,2019,50(3):73-77.

    SHU C L, LI L, FENG T. Simulation of emergent water pollution accident in river-type drinking water sources based on MIKE21: case of Nanchang reach of Ganjiang River[J]. Yangtze River,2019,50(3):73-77.

    [25] 张铃松, 丁文慧, 孟凡生, 等.长江生态环境保护修复咸宁市驻点跟踪研究分析[J]. 环境与可持续发展,2019,44(5):48-50.

    ZHANG L S, DING W H, MENG F S, et al. Research and analysis of stationary tracking in Xianning City of Yangtze River ecological environment protection and remediation[J]. Environment and Sustainable Development,2019,44(5):48-50.

    [26] 王金南, 孙宏亮, 续衍雪, 等.关于“十四五”长江流域水生态环境保护的思考[J]. 环境科学研究,2020,33(5):1075-1080.

    WANG J N, SUN H L, XU Y X, et al. Water eco-environment protection framework in the Yangtze River Basin during the 14th Five-Year Plan Period[J]. Research of Environmental Sciences,2020,33(5):1075-1080.

    [27] 张西斌.铜陵淡水豚国家级自然保护区长江段物种资源现状调查研究[J]. 现代农业科技,2017(11):227-229. DOI: 10.3969/j.issn.1007-5739.2017.11.141
    [28]

    WARREN I R, BACH H K. MIKE 21: a modelling system for estuaries, coastal waters and seas[J]. Environmental Software,1992,7(4):229-240. DOI: 10.1016/0266-9838(92)90006-P

    [29] 王青青, 吴昊, 王腊春, 等.突发污染事故对水源地的影响研究: 以长江南京段为例[J]. 生态与农村环境学报,2019,35(12):1557-1563.

    WANG Q Q, WU H, WANG L C, et al. Study on the impact of sudden pollution on water resources: a case study of Nanjing section of Yangtze River[J]. Journal of Ecology and Rural Environment,2019,35(12):1557-1563.

    [30] 刘晨辉, 刘思飔, 李丹, 等.基于MIKE21模型的长江中环排污口水质影响分析[J]. 中国农村水利水电,2020(1):72-76. DOI: 10.3969/j.issn.1007-2284.2020.01.016

    LIU C H, LIU S S, LI D, et al. An analysis of the impact of Zhonghuan sewage outlet on water quality of the Yangtze River based on MIKE21 model[J]. China Rural Water and Hydropower,2020(1):72-76. DOI: 10.3969/j.issn.1007-2284.2020.01.016

    [31] 龙天渝, 屈瑶.长江重庆主城区段突发水污染事故模拟研究[J]. 重庆大学学报,2020,43(10):90-103.

    LONG T Y, QU Y. Simulation of an emergent water pollution accident in the portion of the Yangtze River running through the main urban areas of Chongqing[J]. Journal of Chongqing University,2020,43(10):90-103.

    [32] 刘晓东.长江下游水域码头硫酸泄漏事故风险评价[J]. 水资源保护,2009,25(3):76-79. DOI: 10.3969/j.issn.1004-6933.2009.03.020

    LIU X D. Risk assessment of vitriol leakage accident at wharf in downriver reach of Yangtze River[J]. Water Resources Protection,2009,25(3):76-79. DOI: 10.3969/j.issn.1004-6933.2009.03.020

    [33] 生态环境部. 建设项目环境风险评价技术导则: HJ 169—2018[S]. 北京: 中国环境出版社, 2019.
  • 期刊类型引用(10)

    1. 耿玉莹,钱秀红,金炳琪,潘维艳,徐征和,窦祥洲. 基于MIKE21的小清河感潮河段水质数值模拟研究. 人民黄河. 2024(10): 96-102 . 百度学术
    2. 文生垚. 水利水电工程施工期污染水体泥水同步净化技术研究. 环境科学与管理. 2024(10): 102-105 . 百度学术
    3. 任昶宇,刘强,张璐涛,史箴,曹阳,俸强,席英伟. 沱江三氯甲烷浓度异常事件中污染物趋势分析. 中国沼气. 2023(02): 59-64 . 百度学术
    4. 洪慧,张强,裴云霞,包美玲,邓俊,陈岷轩. 生态环境损害鉴定证据采集与认定研究. 环境污染与防治. 2023(05): 724-729+735 . 百度学术
    5. 马铁梅,付俊峰,孙英. 基于MIKE21的异龙湖水动力优化与水质改善研究. 水电能源科学. 2023(07): 54-58 . 百度学术
    6. 杨啸宇,王欣,丁磊,王逸飞,周玲霞,缴健,臧英平. 长江南京段不同径流—潮汐动力作用下突发污染事故对饮用水水源地的影响研究. 环境监控与预警. 2023(04): 45-51 . 百度学术
    7. 李建勋,仇池,王超,李庆. 代表性突发水污染情景的层级化涌现机制. 环境工程. 2023(S2): 1167-1172+1177 . 百度学术
    8. 吴锡峰. 九龙江雁石溪“一河一策一图”环境应急响应方案技术研究. 海峡科学. 2023(12): 89-93 . 百度学术
    9. 王盼新,宋玉栋,吴昌永,袁玥,黄琪,付丽亚,毛率先. 关于突发水污染事件应急处置技术需求的几点思考. 环境工程技术学报. 2022(06): 1972-1977 . 本站查看
    10. 陈铭瑞,靳燕国,刘爽,张召,雷晓辉,朱杰. 明渠突发水污染事故段及下游应急调控. 南水北调与水利科技(中英文). 2022(06): 1188-1196 . 百度学术

    其他类型引用(7)

图(9)  /  表(1)
计量
  • 文章访问数:  521
  • HTML全文浏览量:  344
  • PDF下载量:  37
  • 被引次数: 17
出版历程
  • 收稿日期:  2021-11-10
  • 网络出版日期:  2022-04-01
  • 刊出日期:  2022-03-19

目录

/

返回文章
返回