海洋水文学

长江口宝钢码头溢油事故油膜漂移扩散数值模拟

  • 顾靖华 ,
  • 朱建荣 ,
  • 金智
展开
  • 河口海岸学国家重点实验室(华东师范大学), 上海 200241
朱建荣。e-mail:

顾靖华(1978—), 男, 上海市人, 工程师, 从事河口海岸学研究。email:

Copy editor: 林强

收稿日期: 2022-04-03

  修回日期: 2022-05-05

  网络出版日期: 2022-05-07

基金资助

上海市科委重点项目(21JC1402500)

Numerical simulation of oil film drift and diffusion after oil spill accident at the Baosteel wharf in the Changjiang Estuary

  • GU Jinghua ,
  • ZHU Jianrong ,
  • JIN Zhi
Expand
  • State Key Laboratory of Estuarine and Coastal Research (East China Normal University), Shanghai 200241, China
ZHU Jianrong. email:

Copy editor: LIN Qiang

Received date: 2022-04-03

  Revised date: 2022-05-05

  Online published: 2022-05-07

Supported by

Science and Technology Commission of Shanghai Municipality(21JC1402500)

摘要

溢油事故发生后, 油膜的漂移扩散会对周围水域和环境敏感目标造成污染。研究溢油事故后油膜漂移扩散, 可为溢油事故的处理提供理论指导。本文应用河口海岸三维水动力模式ECOM-si(semi-implicit estuarine, costal and ocean model), 耦合溢油漂移扩散模块, 模拟和分析长江口宝钢码头发生溢油事故后油膜的漂移扩散, 以及对环境敏感保护目标的影响。基于长江口崇明东滩气象站实测风速风向资料, 给出影响溢油漂移的主导风和不利风向。在冬季多年平均1月径流量11700m³·s-1、主导风、3个不利风和潮汐作用下, 数值模拟并详细给出了宝钢码头溢油事故发生后油膜的平面分布、油膜到达和离开4个水库取水口和饮用水水源保护区以及其他环境敏感保护区的时间、持续影响时间和最大油膜厚度。长江口宝钢码头溢油事故发生后, 油膜随涨落潮流作上下游来回振荡, 径流使油膜向海输运, 风使油膜朝风向方向漂移。在主导风北风5.6m³·s-1风速下, 油膜沿长江口南岸向下游漂移扩散, 小部分进入北槽南侧。在不利风向东南风4.0m³·s-1风速下, 油膜西北方向漂移, 聚集于南支北岸, 受径流作用沿南支和北港的中北侧向下游输运。在不利风向西北风4.8m³·s-1风速下, 西北风减轻了溢油事故点上游和北港、崇明东滩外侧敏感目标的影响, 加重了南港和南槽的环境敏感保护目标的影响。在不利风向西南风3.2m³·s-1风速下, 大部分油膜在北港中北侧向下游漂移扩散, 小部分油膜聚集在北槽的中北侧。西南风减轻了溢油事故点上游和下游南港、南槽环境敏感目标受到的影响, 加重了对北港和崇明东滩外侧敏感目标的影响。不同风向作用下油膜的分布和对环境敏感保护目标影响显著不同, 风在油膜漂移扩散中起着十分重要的作用。

本文引用格式

顾靖华 , 朱建荣 , 金智 . 长江口宝钢码头溢油事故油膜漂移扩散数值模拟[J]. 热带海洋学报, 2022 , 41(6) : 159 -170 . DOI: 10.11978/2022065

Abstract

The drift and diffusion of oil film after the oil spill accident will pollute the surrounding waters, and environmental sensitive protection targets. The study of oil film drift and diffusion after oil spill accident can provide theoretical guidance for the treatment of oil spill accident. Coupling the oil spill drift and diffusion module, the estuarine and coastal three-dimensional hydrodynamic model ECOM-si (semi-implicit estuarine, costal and ocean Model) was used to simulate and analyzed the drift and diffusion of oil film after the oil spill accident at the Baosteel wharf in the Changjiang Estuary, and its impact on the environmental sensitive protection targets in this paper. Based on the measured wind speed and direction data at weather station on the Chongming eastern shola in the Changjiang Estuary, the prevailing wind and unfavorable wind affecting the oil spill drift are given. Forced by the monthly mean river discharge of 11700 m³·s-1, prevailing wind, three unfavorable winds and tide in January, the horizontal distribution of oil film, time of the oil film reaching and leaving the water intakes of four reservoirs and drinking water source protection areas, and other environmental sensitive protection targets, duration and maximum oil film thickness of oil film after the oil spill accident at the Baosteel wharf was numerically simulated and analyzed in detail. After the oil spill accident, the oil film oscillates upstream and downstream with the flood and ebb current, meanwhile the runoff transports the oil film to the sea, and the wind makes the oil film drift towards the wind direction. Under the action of the prevailing wind, northerly wind of 5.6 m³·s-1, the oil film drifts and diffuses downstream along the south coast of the Changjiang Estuary, and a small part enters the south side of the North Passage. Under the unfavorable wind, southeasterly wind of 4.0 m³·s-1, the oil film drifts in the northwest and accumulates on the north coast of the South Branch. Forced by the runoff, it is transported downstream on the middle and north side of the South Branch and the North Channel. Under the unfavorable wind, northwesterly wind of 4.8 m³·s-1, it weakens the impact on the sensitive targets upstream of the oil spill accident site and the ones in the North Channel, outside Chongming eastern shoal, and aggravates the impact on the environmental sensitive protection targets in the North Channel and South Passage. Under the unfavorable wind, southwesterly wind of 3.2 m³·s-1, most of the oil film drifts and diffuses downstream on the middle and north side of the North Channel, and a small part of the oil film is on the middle and north side of the North Passage. It weakens the impact on environmental sensitive targets upstream of the oil spill accident site and the downstream ones in the South Channel and South Passage, and aggravates the impact on the sensitive targets in the North Channel and outside Chongming eastern shoal. The distribution of oil film and its impact on the environmental sensitive protection targets are significantly different under the force of different wind directions, indicating that wind plays a very important role in oil film drift and diffusion.

码头、港口、航道和桥墩等邻近水域是溢油事故多发之处。溢油事故发生后, 油膜的漂移扩散会对周围水域和敏感目标造成污染。因此, 研究溢油事故后的油膜漂移扩散, 可为溢油事故的处理提供油膜面积、到达敏感目标时间和油膜厚度变化等数据, 对于事故应对和环境保护具有实际意义。
研究溢油事故发生后油膜运动的主要方法有现场实验、水槽试验、理论分析和数值模拟。溢油事故后油膜漂浮于水面, 其运动主要是由水流的平流, 以及风引起的漂移和随机扩散决定, 还受蒸发、分散、乳化、沉降及生物降解等过程的影响。张和庆 等(2001)在珠江口通过现场溢油实验, 分析了油类入海后漂流和扩散过程。郭运武 等(2008)通过水槽实验, 分析了在不同风速和水流速度条件下油膜的扩展、漂移运动状况, 结果表明风在溢油油膜扩展、漂移中起着重要的作用。刘栋 等(2006)通过潮汐水槽模拟实验, 观察和记录油膜在潮汐水流中两个半日潮潮周期内的扩展、漂移情况, 分析了溢油油膜在潮汐水流中的扩展、漂移过程。娄厦 等(2008)理论分析了油在水中的行为和归宿, 概括了多个溢油模型的理论依据、基本假设和应用范围。Violeau等(2007)通过室内试验和应用平滑粒子水动力拉格朗日方法模拟溢油的漂移扩散。在溢油运动的研究和预测中, 数值模拟是一种应用最为广泛的研究方法。较早的一些模式应用垂向平均的二维模式较多, 例如娄安刚 等(2001)建立了胶州湾内的二维预报潮流场, 再考虑风场对海面溢油的影响, 模拟了溢油事故后溢油的油膜漂移轨迹。孙长青 等(2003)建立渤海湾二维溢油飘移数值模型, 计算在不同风场和潮流场的共同作用下, 溢油飘移的输移轨迹、可能抵岸的时间、地点及残留量等。李大鸣 等(2008)建立了在潮汐、风浪作用下油膜运动二维数学模型, 对大亚湾海区油码头瞬时溢油事故后的影响范围进行了模拟。Periáñez等(2008)建立了二维水动力模式, 模拟直布罗陀海峡潮流和余流, 并质点跟踪技术模拟潮流、混合和风对溢油运动的影响。垂向平均的二维数值模式不能真实反映风生流和密度流等在垂向上的变化, 从而影响模拟油膜漂移扩散的精度。因此, 将三维模式运用到油膜漂移扩散的模拟中, 能较为真实地模拟风生流和密度流等的垂向变化, 并更为准确地模拟表层流场对溢油漂移扩散的作用。例如, Wang等(2008)采用粒子方法在二维溢油模式的基础上发展三维溢油模式, 应用到渤海溢油事故的数值模拟。应用三维数值模式, 陈义中 等(2006)模拟了渤海原油码头溢油漂移扩散, 傅利辉 等(2008)模拟了长江口北支崇启大桥建成后流场变化及溢油漂移扩散, 陈士谦 等(2012)模拟了鸭绿江公路大桥溢油的漂移扩散。
长江口为三级分汊四口入海的大型潮汐河口, 环境敏感保护目标众多(表1)。环境敏感保护目标位于溢油事故发生点的不同方向, 不同的方向对应不同的敏感目标。数值预测油膜运动需要考虑主导风和不利风, 其对油膜的运动起着十分重要的作用。对溢油事故附近气象站实测风作统计分析, 可得出各个方向的频率和平均风速, 第一频率的风取为主导风。不利风向取决于溢油事故发生点和敏感目标的连线, 与指向敏感目标最接近的风向为不利风向。主导风和不利风可从实测风的统计结果给出。本文采用崇明东滩气象站实测的风资料(测站见图1), 给出主导风和不利风向。在此基础上, 应用长期改进和应用的河口海岸水动力三维数值模式ECOM-si(semi-implicit estuarine, costal and ocean model), 耦合溢油漂移扩散模块, 模拟并分析上海港宝钢原料区域码头船舶发生溢油事故后油膜的漂移扩散(溢油事故点位置见图1), 以及对环境敏感保护目标(图2)的影响。
表1 长江口环境敏感保护目标

Tab. 1 Environmental sensitive protection targets in the Changjiang Estuary

序号 保护目标名称 保护目标级别
1 长江刀鲚国家级水产种质资源保护区 国家级
2 青草沙饮用水水源保护区 一级保护区、二级保护区
3 陈行饮用水水源保护区 一级保护区、二级保护区
4 长江太仓浏河饮用水水源保护区 一级保护区、二级保护区
5 东风西沙饮用水水源保护区 一级保护区、二级保护区
6 上海吴淞炮台湾国家湿地公园 国家级
7 上海崇明东滩鸟类国家级自然保护区 国家级
8 上海九段沙湿地国家级自然保护区 国家级
9 上海市长江口中华鲟省级自然保护区 省级
图1 长江口形-势图及上海港宝钢原料区域码头溢油事故泄漏点

基于国家测绘地理信息局标准地图服务网站下载审图号为GS(2016)2556 标准地图制作。图中绿点为崇明东滩气象站; 红点为泄漏点; 黑点为水库取水口

Fig. 1 Sketch map of the Changjiang Estuary, and the site of oil leakage accident at the Baosteel raw material regional wharf of the Shanghai Port. The green dot marks the weather station at the Chongming eastern shoal; the red dot represents the oil leakage point, black dots stand for the water intakes of the reservoirs

图2 长江口环境敏感保护目标位置

基于国家测绘地理信息局标准地图服务网站下载审图号为GS(2016)2556 标准地图制作

Fig. 2 Positions of the environmental sensitive protection targets in the Changjiang Estuary

1 数值模式设置

采用作者课题组长期改进和应用的河口海岸水动力三维数值模式ECOM-si, 该模式在长江口水动力过程和盐水入侵等方面研究中经过了大量的率定和验证, 得到了广泛的应用, 取得了大量研究成果(Lyu et al, 2018; Chen et al, 2019; Zhu et al, 2020; Ma et al, 2022) 。在ECOM-si的基础上, 耦合油膜漂移扩散模块。该模块将油分成有限个粒子, 每个粒子的位置按以下公式计算:
$X={{X}_{0}}+\left( U+\alpha {{W}_{10}}\cos A+r\cos B \right)\Delta t$
$Y={{Y}_{0}}+\left( V+\alpha {{W}_{10}}\sin A+r\sin B \right)\Delta t$
式中X0Y0为某粒子的上一时刻坐标; XY为某质点的下一时刻坐标; UV分别为XY方向的流速分量, 由水动力数值模式计算得出; W10为海面上10m处的风速; A为风向; α为风漂移系数。在风的影响下, 油膜漂移速度的增加量约为风速的2%~3%, 漂移方向与风向成0°~40°夹角。由于油膜中心的漂移速度和方向是表面水流和风所引起的流速矢量之和, 在模拟中使用风漂移系数α为风与海流的角度修正系数, 一般在0.02~0.03之间, 本文取0.03; r为随机扩散项, r=RE, R为0~1之间的随机数, E为扩散系数, 计算公式为$E=\sqrt{6A_{\alpha}\Delta t}$, 其中Aα为水平扩散系数, 由水动力模型提供; B为随机扩散方向, 计算公式为B=2πR
在油膜的漂移扩散过程中, 还需考虑油膜的蒸发。油膜蒸发受溢油种类、气温、溢油面积、风速、太阳辐射和油膜厚度等因素的影响。本模块考虑溢油的蒸发公式为(Stiver et al,1984)
${{F}_{\text{v}}}=\ln \left[ 1+{{B}_{1}}\left( \frac{{{T}_{\text{G}}}}{T} \right)\theta \exp \left( {{A}_{1}}-\frac{{{B}_{1}}{{T}_{0}}}{T} \right) \right]\left[ \frac{T}{{{B}_{1}}{{T}_{\text{G}}}} \right]$
其中, Fv为溢油蒸发体积分数; T是环境温度(空气温度), 本文长江口冬季取283K; θ为无因次时间(蒸发暴露时间); A1B1T0TG为各种依赖于油品的常数, A1=6.3, B1=10.3; T0为油在Fv=0时的初始沸点温度(单位: K); TG为蒸馏曲线的坡度(单位: K); T0TG根据下列经验公式计算:
${{T}_{0}}=532.98-3.1295*\text{API}$
${{T}_{\text{G}}}=985.62-13.597*\text{API}$
其中API为比重指数, 对重油、重油和轻油可分别取值为16.9162、26.8458和36.1882。
蒸发公式中的θ由下式计算:
$\theta =K{{A}_{2}}t/V$
其中K是空气端传质系数, A2是溢油面积, t是时间, V是溢油泄漏的体积。针对公式实际使用中蒸发过快的现象, Fingas(1995)对K值进行了调整, 为了解决无风速情况下传质系数公式无法计算的问题, 孙宝楠 等(2010)假设在传质系数上加一个小量R, 从而传质系数可以表示为
$K=(0.0025-0.000021t)U_{10}^{0.78}+{{R}_{1}}$
其中U10为海表面10m风速。R1的范围在(0.1~10)×10-4m·s-1的量级(本文K值取10-4), t为溢油发生后的时间, 单位为h。每个油粒子初始时刻的体积(总溢油体积/划分总粒子数)减去t时刻的溢油蒸发量, 即为t时刻每个油粒子的体积。
溢油事故发生点为上海港宝钢原料区域码头石灰石及砂岩卸船码头, 坐标31°28′6″N, 121°25′33″E, 泄漏量1000t, 0#轻质柴油。在径流作用下, 油膜向海方向输运, 落憩时刻之后为涨潮, 相对于涨憩时刻, 油膜会有更长时间滞留在口门内, 为不利情况, 故假设溢油事故发生在大潮期间的落憩时刻。水动力模式先运行10d, 水位和流场准稳定后, 再释放油粒子。只要质点数足够多, 油膜漂移的模拟就越精确, 因此本文设置了一千万个油膜质点。每个质点代表一个恒定的体积大小, 每步根据网格的水平面积大小和所含质点数, 可以求出油膜在该网格的厚度值: 厚度=质点数×t时刻单个质点代表的体积/网格的水平面积。
模型计算范围包括整个长江河口、杭州湾和邻近海域(图3a), 东至125°E附近, 北至接近34°N, 南至27°30′N左右。模型水平方向上采用曲线非正交网格, 较好地拟合了长江河口的岸线, 并对南北支分汊口、深水航道等区域加密(图3b、c)。长江口内网格分辨率最高可至约100m, 长江口外网格分辨率最高可至10km左右。模型在垂直方向上采用σ坐标, 由表至底均匀分成10层。计算时间步长取为60s。由于长江河口潮滩分布较广, 模型采用干湿网格判断实现潮滩动边界, 最小临界水深阈值设置为0.2m。岸线和水深数据为2017年长江河口实测资料, 水深数据以黄海85高程为基面, 见图1
图3 模型计算范围和网格(a)以及局部放大的南北支分汊口附近网格(b)和南支口门附近网格(c)

红线为深水航道导堤和丁坝。基于国家测绘地理信息局标准地图服务网站下载审图号为GS(2016)2556标准地图制作

Fig. 3 Model calculation domain and grids (a), and the enlarged view of grids near the bifurcation between the South Branch and North Branch (b) and near the mouth of the South Branch. The red lines are the guide dikes and groins of the Deepwater Channel

本文模拟和分析冬季1月径流量和主导风、不利风作用下油膜的漂移扩散。选取大通作为上游边界, 模式的径流边界条件可直接采用大通水文站的月平均径流量资料。1950—2020年多年月平均长江入海径流量1月为11700m3·s-1。外海开边界潮位考虑由16个分潮(M2, S2, N2, K2, K1, O1, P1, Q1, MU2, NU2, T2, L2, 2N2, J1, M1和OO1) 的调和常数合成得到, 资料来自全球潮汐模式NAOTIDE(http://www.miz.nao.ac.jp)计算结果。
对主导风和不利风向, 本文采用长江口崇明东滩气象站(测站位置见图1)2007年至2021年14年的实测风速风向资料, 统计1月8个罗盘方向风的频率和平均风速。1月第一至第八风向频率分别为26.9%、26.1%、16.5%、8.3%、7.7%、6.3%、5.4%、2.8%, 风向分别是N、NW、NE、E、SE、W、S、SW, 风速分别是5.6、4.8、5.0、4.1、4.1、3.6、4.0、3.4m·s-1图4为1月风玫瑰图, 形象给出了8个方向风的频率和平均风速, 以及各个方向上不同风速(间隔2m·s-1)及其对应的频率。
图4 崇明东滩气象站2007—2021年统计得出的1月风玫瑰图

Fig. 4 Wind roses in January calculated from 2007 to 2021 at weather station on Chongming eastern shoal

主导风取为第一频率的风。1月第一频率(26.9%)风向为北风, 风速5.6m·s-1。不利风由溢油事故发生点和环境敏感保护目标位置确定。由于长江口环境敏感保护目标众多(图2), 需要给出多个不利风向。上海港宝钢原料区域码头船舶发生溢油事故发生点的上游敏感目标为陈行水库、太仓水库、东风西沙水库、刀鲚核心区和实验区, 不利风向为东南风, 为第5频率(7.7%), 风速4.0m·s-1。对下游东南向敏感目标为青草沙水库、刀鲚实验区、吴淞湿地和九段沙湿地, 不利风向为西北风, 为第2频率(26.1%), 风速4.8m·s-1。对下游东北向敏感目标为青草沙水库、刀鲚实验区、崇明东滩保护区和中华鲟保护区, 不利风向为西南风, 为第8频率(2.8%), 风速3.2m·s-1

2 油膜漂移扩散的数值模拟

本文模拟和分析上海港宝钢原料区域码头综合改造工程码头溢油事故发生后油膜的输运扩散及其对敏感目标的影响, 对模拟结果绘制油膜分布和水源地取水口油膜厚度随时间变化过程图进行分析。

2.1 主导风北风5.6m·s-1风速

油膜的运动主要受水流的平流、风的漂移和随机扩散作用。本文假设宝钢码头的溢油事故发生在大潮落憩时刻, 长江河口为半日潮, 事故发生后3h内为涨潮流, 在潮流的平流作用下向上游运动(图5)。油膜中心位于陈行水库水源地保护区, 厚度超过2μm。之后在径流的作用下, 油膜在随涨落潮流上下移动过程中向下游输运, 至12h油膜已进入南港, 中心位于外高桥附近(图6)。至72h, 油膜继续向下游漂移扩散, 大部分位于南港中下段和南槽中上段的南侧, 小部分进入北槽南侧(图7)。在油膜的运动过程中, 绝大部分位于长江口南岸, 在主导风北风的作用下向南漂移起着重要的作用。
图5 溢油事故后3小时油膜分布

基于国家测绘地理信息局标准地图服务网站下载审图号为GS(2016)2556 标准地图制作

Fig. 5 Distribution of oil film thickness after 3 hours of the oil spill accident

图6 溢油事故后24小时油膜分布

基于国家测绘地理信息局标准地图服务网站下载审图号为GS(2016)2556 标准地图制作

Fig. 6 Distribution of oil film thickness after 24 hours of the oil spill accident

图7 溢油事故后72小时油膜分布

基于国家测绘地理信息局标准地图服务网站下载审图号为GS(2016)2556 标准地图制作

Fig. 7 Distribution of oil film thickness after 72 hours of the oil spill accident

从长江口水库取水口油膜厚度随时间变化过程看(图8), 因受径流和北风的作用, 溢油事故对南支上端北侧的东风西沙水库取水口无影响。太仓水库位置南支中段南侧, 离溢油事故发生点较远, 同样受径流和北风的作用, 油膜最大厚度小于1μm, 影响微小。青草沙水库取水口位于溢油点东北侧不远处, 受北风作用, 油膜最大厚度小于1μm, 影响微小。陈行水库取水口位于溢油点上游不远处, 油膜2h随涨潮流到达取水口, 之后随着潮流周期性起伏, 最大厚度约为68μm。31h后, 油膜的影响消除, 期间最长持续影响10h, 总影响时间为29h。
图8 油膜厚度在东风西沙水库取水口(a)、太仓水库取水口(b)、陈行水库取水口(c)和青草沙水库取水口(d)随时间变化情况

黑线: 主导风(北风); 红线: 不利风向(东南风); 绿线: 不利风向(西北风); 紫线: 不利风向(西南风)

Fig. 8 Temporal variation of oil film thickness at the water intakes of the Dongfengxisha reservoir (a), Taicang reservoir (b), Chenhang reservoir (c) and Qingcaosha reservoir. Black line: prevailing wind (northerly wind); red line: unfavorable wind direction (southeasterly wind); green line: unfavorable wind direction (northwesterly wind); purple line: unfavorable wind direction (southwesterly wind)

对青草沙饮用水水源保护区, 溢油事故发生后11h油膜到达, 36h后离开, 持续影响25h, 期间最大厚度为1.42μm。对陈行饮用水水源保护区, 溢油事故发生后1.5h油膜到达, 53h后离开, 持续影响51.5h, 期间最大厚度为151.31μm。对长江太仓浏河饮用水水源保护区, 溢油事故发生后3.5h油膜到达, 5.5h后离开, 持续影响2h, 期间最大厚度为2.99μm。溢油事故对东风西沙饮用水水源保护区都无影响。
对长江刀鲚国家级水产种质资源核心保护区, 溢油事故发生后2.5h油膜到达, 29h后离开, 持续影响26.5h, 期间最大厚度为47.84μm。对长江刀鲚国家级水产种质资源实验保护区, 溢油事故发生后0.5h油膜到达, 影响时间超过溢油发生后的72h, 期间最大厚度为52.99μm。对上海吴淞炮台湾国家湿地公园, 溢油事故发生后8h油膜到达, 持续影响时间超过第72h, 期间最大厚度为32.70μm。对上海九段沙湿地国家级自然保护区, 溢油事故发生后44.5h油膜到达, 影响时间超过溢油发生后的72h, 期间最大厚度为1.47μm。溢油事故对上海崇明东滩鸟类国家级自然保护区、上海市长江口中华鲟省级自然保护区无影响, 缘于它们远离溢油点, 位于东北的北港口门北侧。

2.2 不利风向东南风4.0m·s-1风速

在不利风向东南风4.0m·s-1作用下, 油膜朝西北方向漂移, 聚集于南支北岸, 同时受径流作用向下游输运, 在溢油事故后72h油膜位于南支的上下段的中北侧和北港上段的中北侧(图9)。这与主导风北风作用下油膜的分布十分不同, 表明风在油膜漂移扩散中起着十分重要的作用, 迫使油膜克服径流朝风向的方向漂移。
图9 溢油事故后72h油膜分布

基于国家测绘地理信息局标准地图服务网站下载审图号为GS(2016)2556 标准地图制作

Fig. 9 Distribution of oil film thickness at 72 hours after the oil spill accident

从长江口水库取水口油膜厚度随时间变化过程看(图8), 对东风西沙水库取水口, 溢油事故发生后16h油膜到达, 影响时间超过溢油发生后的72h, 期间最大厚度为10μm。与主导风北风结果相比, 东南风使得油膜向西北漂移更远的距离。对太仓水库取水口, 溢油事故发生后4h油膜到达, 8h后离开, 持续影响4h, 期间最大厚度为1.8μm。对陈行水库取水口, 溢油事故发生后2h油膜到达, 12h后离开, 持续影响10h, 期间最大厚度为7μm。对青草沙水库取水口, 溢油事故发生后10h油膜到达, 39h后离开, 持续影响29h, 期间最大厚度为2.3μm。
对青草沙饮用水水源保护区, 溢油事故发生后9h油膜到达, 影响时间超过溢油发生后的72h, 期间最大厚度为10.45μm。对陈行饮用水水源保护区, 溢油事故发生后1h油膜到达, 影响时间超过溢油发生后的72h, 期间最大厚度为149.46μm。对长江太仓浏河饮用水水源保护区, 溢油事故发生后3h油膜到达, 55.5h后离开, 持续影响52.5h, 期间最大厚度为52.27μm。对东风西沙饮用水水源保护区, 溢油事故发生后17.5h油膜到达, 影响时间超过溢油发生后的72h, 期间最大厚度为37.89μm。
对长江刀鲚国家级水产种质资源核心保护区, 溢油事故发生后2h油膜到达, 影响时间超过溢油发生后的72h, 期间最大厚度为84.47μm。对长江刀鲚国家级水产种质资源实验保护区, 溢油事故发生后0.5h油膜到达, 影响时间超过溢油发生后的72h, 期间最大厚度为110.00μm。溢油事故对上海吴淞炮台湾国家湿地公园、上海崇明东滩鸟类国家级自然保护区、上海九段沙湿地国家级自然保护区、上海市长江口中华鲟省级自然保护区均无影响。这与油膜的平面分布图是一致的。

2.3 不利风向西北风4.8m·s-1风速

溢油事故后72h油膜分布与主导风北风作用下的分布相似, 这里不再给出。
从长江口水库取水口油膜厚度随时间变化过程看(图8), 西北风阻止了油膜向上游漂移, 溢油事故对东风西沙、太仓水库取水口无影响, 对陈行水库和青草沙水库影响微小, 最大油膜厚度小于1μm。与主导风北风结果相比, 西北风与北风作用下陈行水库取水口油膜厚度差别显著。
对青草沙饮用水水源保护区, 溢油事故发生后3h油膜到达, 23h后离开, 持续影响20h, 期间最大厚度为5.12μm。对陈行饮用水水源保护区, 溢油事故发生后10h油膜到达, 12.5h后离开, 持续影响2.5h, 期间最大厚度为4.90μm。对长江太仓浏河饮用水水源保护区和东风西沙饮用水水源保护区无影响。
对于其他环境敏感保护目标, 西北风减轻了溢油事故点上游和北港、崇明东滩外侧敏感目标的影响, 但加重了南港和南槽的环境敏感保护目标的影响。

2.4 不利风向西南风3.2m·s-1风速

溢油事故发生后, 油膜受涨潮流的作用向上游运动, 同时受西南风的作用向东北漂移, 大部分油膜绕过青草沙水库的西南角, 进入北港。在溢油事故后72h(图10), 油膜受向海的径流和西南风的作用, 大部分油膜在北港向下游漂移扩散, 位于崇明岛南岸的北港中段中北侧, 小部分油膜在北槽上端的中北侧。可见, 西南风使油膜漂移至北港, 并和径流一起作用向海输运。
图10 溢油事故后72h油膜分布

基于国家测绘地理信息局标准地图服务网站下载审图号为GS(2016)2556 标准地图制作

Fig. 10 Distribution of oil film thickness at 72 hours after the oil spill accident

从长江口水库取水口油膜厚度随时间变化过程看(图8), 溢油事故对东风西沙水库和太仓水库取水口无影响。对陈行水库取水口, 溢油事故发生后2h油膜到达, 8h后离开, 持续影响6h, 期间最大厚度为6.3μm。对青草沙水库取水口, 溢油事故发生后9h油膜到达, 54h后离开, 持续影响45h, 期间最大厚度为8.5μm。
对青草沙饮用水水源保护区, 溢油事故发生后8h油膜到达, 67h后离开, 持续影响59h, 期间最大厚度为17.79μm。对陈行饮用水水源保护区, 溢油事故发生后1h油膜到达, 41h后离开, 持续影响40h, 期间最大厚度为122.60μm。对长江太仓浏河饮用水水源保护区, 溢油事故发生后3.5h油膜到达, 6.5h后离开, 持续影响3h, 期间最大厚度为15.02μm。溢油事故对东风西沙饮用水水源保护区无影响。
对其他环境敏感保护目标的影响, 西南风减轻了溢油事故点上游和下游南港、南槽环境敏感目标的影响, 但加重了北港和崇明东滩外侧敏感目标的影响。

3 结论

本文应用河口海岸三维水动力模式ECOM-si, 耦合溢油漂移扩散模块, 模拟和分析长江口宝钢码头船舶发生溢油事故后油膜的漂移扩散, 以及对环境敏感保护目标的影响。基于长江口崇明东滩气象站2007—2021共14年的实测风速风向资料, 给出影响溢油漂移的主导风和不利风向。在冬季多年平均1月径流量11700m³·s-1和主导风、3个不利风作用下, 数值模拟和详细给出了上海港宝钢原料区域码头溢油事故发生后油膜的平面分布、油膜到达和离开4个水库取水口和饮用水水源保护区以及其他环境敏感保护区的时间、持续影响时间和最大油膜厚度。
长江口宝钢码头溢油事故发生后, 油膜的运动主要受水流的平流、风的漂移和随机扩散作用。在1月主导风北风5.6m·s-1风速作用下, 油膜沿长江口南岸向下游漂移扩散, 在溢油事故后第72h大部分位于南港中下段和南槽中上段的南侧, 小部分进入北槽南侧。在油膜的运动过程中, 绝大部分位于长江口南岸, 在主导风北风的作用下向南漂移起着重要的作用。在不利风向东南风4.0m·s-1风速作用下, 油膜西北方向漂移, 聚集于南支北岸, 同时受径流作用向下游输运, 在溢油事故后72h油膜位于南支的上下段的中北侧和北港上段的中北侧。在不利风向西北风4.8m·s-1风速作用下, 西北风减轻了溢油事故点上游和北港、崇明东滩外侧敏感目标受到的影响, 加重了对南港和南槽的环境敏感保护目标的影响。在不利风向西南风3.2m·s-1风速作用下, 在溢油事故后72h大部分油膜在北港向下游漂移扩散, 位于崇明岛南岸的北港中段中北侧, 小部分油膜在北槽上端的中北侧。西南风减轻了溢油事故点上游和下游南港、南槽环境敏感目标的影响, 加重了北港和崇明东滩外侧敏感目标的影响。油膜随涨落潮流作上下游来回振荡过程中, 径流使油膜向海输运, 风使油膜朝风向方向漂移。不同风向作用下油膜的分布和对环境敏感保护目标影响显著不同, 表明风在油膜漂移扩散中起着十分重要的作用。
[1]
陈士谦, 朱建荣, 2012. 鸭绿江公路大桥溢油漂移扩散三维数值模拟[J]. 华东师范大学学报(自然科学版), (6): 46-56.

CHEN SHIQIAN, ZHU JIANRONG, 2012. Three- dimensional numerical simulation of the drift and diffusion of oil film at the highway bridge in the Yalu River[J]. Journal of East China Normal University (Natural Science), (6): 46-56. (in Chinese with English abstract)

[2]
陈义中, 朱建荣, 陈昞睿, 等, 2006. 渤海原油码头溢油漂移扩散的数值模拟[J]. 华东师范大学学报(自然科学版), (6): 25-33.

CHEN YIZHONG, ZHU JIANRONG, CHEN BINGRUI, et al, 2006. Numerical simulation of the drift and diffusion of oil film at the crude oil wharf in the Bohai Sea[J]. Journal of East China Normal University (Natural Science), (6): 25-33. (in Chinese with English abstract)

[3]
傅利辉, 朱建荣, 陈昞睿, 等, 2008. 崇启大桥建成后流场变化及溢油的数值模拟. 华东师范大学学报(自然科学版), (4): 114-123.

FU LIHUI, ZHU JIANRONG, CHEN BINGRUI, et al, 2008. Numerical simulation of current change and diffusion of oil film after the Chongqi bridge construction[J]. Journal of East China Normal University (Natural Science), (4): 114-123. (in Chinese with English abstract)

[4]
郭运武, 刘栋, 钟宝昌, 等, 2008. 风对河道溢油扩展、漂移影响的实验研究[J]. 水动力学研究与进展, 23(4): 446-452.

GUO YUNWU, LIU DONG, ZHONG BAOCHANG, et al, 2008. Experimental investigations of wind effect on the oil spill spread and drift[J]. Chinese Journal of Hydrodynamics, 23(4): 446-452. (in Chinese with English abstract)

[5]
李大鸣, 陈海舟, 付庆军, 2008. 海上溢油数学模型的研究与应用[J]. 哈尔滨工程大学学报, 29(12): 1291-1297.

LI DAMING, CHEN HAIZHOU, FU QINJUN, 2008. Research and application of mathematical modeling to oil spills on the sea[J]. Journal of Harbin Engineering University, 29(12): 1291-1297. (in Chinese with English abstract)

[6]
刘栋, 林卫青, 钟宝昌, 等, 2006. 感潮河道溢油扩展、漂移特性实验[J]. 水动力学研究与进展, 21(6): 744-751.

LIU DONG, LIN WEIQING, ZHONG BAOCHANG, et al, 2006. Characteristic experiment of spread and transport of oil spill in tidal river[J]. Journal of Hydrodynamics, 21(6): 744-751. (in Chinese with English abstract)

[7]
娄安刚, 王学昌, 孙长青, 等, 2001. 胶州湾海面溢油轨迹的数值模拟[J]. 黄渤海海洋, 19(1): 1-8.

LOU ANGANG, WANG XUECHANG, SUN CHANGQING, et al, 2001. A simulation of oil spill trajectory on sea surface of Jiaozhou Bay[J]. Journal of Oceanography of Huanghai & Bohai Seas, 19(1): 1-8. (in Chinese with English abstract)

[8]
娄厦, 刘曙光, 2008. 溢油模型理论及研究综述[J]. 环境科学与管理, 33(10): 33-37, 61.

LOU SHA, LIU SHUGUANG, 2008. Review in theory and study of oil spill models[J]. Environmental Science and Management, 33(10): 33-37, 61. (in Chinese with English abstract)

[9]
孙宝楠, 杨永增, 于顺, 等, 2010. 基于分析法经典溢油蒸发模型的改进[J]. 海洋科学进展, 28(3): 364-371.

SUN BAONAN, YANG YONGZENG, YU SHUN, et al, 2010. Improvement for oil spill evaporation model based on classic analysis method[J]. Advances in Marine Science, 28(3): 364-371. (in Chinese with English abstract)

[10]
孙长青, 赵可胜, 郭耀同, 2003. 渤海湾海面溢油数值计算[J]. 海洋科学, 27(11): 63-67.

SUN CHANGQING, ZHAO KESHENG, GUO YAOTONG, 2003. Digital modeling for oil-spill in Bohai Bay[J]. Marine Sciences, 27(11): 63-67. (in Chinese with English abstract)

[11]
张和庆, 李福娇, 2001. 近海海面油类漂流扩散的研究和预测实践[J]. 热带气象学报, 17(1): 83-89.

ZHANG HEQING, LI FUJIAO, 2001. Study and practice of the oil drift and diffusion on the sea surface in the offshore[J]. Journal of Tropical Meteorology, 17(1): 83-89. (in Chinese with English abstract)

[12]
CHEN QING, ZHU JIANRONG, LYU HANGHANG, et al, 2019. Impacts of topography change on saltwater intrusion over the past decade in the Changjiang estuary[J]. Estuarine, Coastal and Shelf Science, 231: 106469.

DOI

[13]
FINGAS M F, 1995. A literature review of the physics and predictive modelling of oil spill evaporation[J]. Journal of Hazardous Materials, 42(2): 157-175.

DOI

[14]
LYU HANGHANG, ZHU JIANRONG, 2018. Impact of the bottom drag coefficient on saltwater intrusion in the extremely shallow estuary[J]. Journal of Hydrology, 557: 838-850.

DOI

[15]
MA RUI, ZHU JIANRONG, 2022. Water and Salt Transports in the Hengsha Channel of Changjiang Estuary[J]. Journal of Marine Science and Engineering, 10(1): 72.

DOI

[16]
PERIÁÑEZ R, PASCUAL-GRANGED A, 2008. Modelling surface radioactive, chemical and oil spills in the Strait of Gibraltar[J]. Computers & Geosciences, 34(2): 163-180.

DOI

[17]
STIVER W, MACKAY D, 1984. Evaporation rate of spills of hydrocarbons and petroleum mixtures[J]. Environmental Science & Technology, 18(11): 834-840.

DOI

[18]
WANG SHOUDONG, SHEN YONGMING, GUO YAKUN, et al, 2008. Three-dimensional numerical simulation for transport of oil spills in seas[J]. Ocean Engineering, 35(5-6): 503-510.

DOI

[19]
VIOLEAU D, BUVAT C, ABED-MERAIM K, et al, 2007. Numerical modelling of boom and oil spill with SPH[J]. Coastal Engineering, 54(12): 895-913.

DOI

[20]
ZHU JIANRONG, CHENG XINYUE, LI LINJIANG, et al, 2020. Dynamic mechanism of an extremely severe saltwater intrusion in the Changjiang estuary in February 2014[J]. Hydrology and Earth System Sciences, 24(10): 5043-5056.

DOI

文章导航

/