2. 国家海洋局东海分局, 上海 200137;
3. 国家海洋局东海环境监测中心, 上海 201206
近年来随着我国沿海经济的迅猛发展,国内对危险化学品的需求及产量每年大约以10%的幅度增加[1],沿海石化区等以化学工业为主导的工业园区建设正在加快。但是同时,近岸临海化工产业逐渐成为危险化学品事故高发区,海上化学品泄漏的环境风险问题已经成为威胁我国海洋生态安全的重要因素之一。危险化学品泄漏等事故的发生往往会导致重大的人身伤亡和财产损失,对周围环境造成严重的破坏,如造成大气污染、水质污染等,不仅直接威胁当地生态安全,也会对周围群众的生活生产带来严重威胁[2]。
宁波-舟山沿海建设有国家石油战略储备基地以及一系列下游化工企业集群,重点发展石油化工和能源等工业,环境风险形势极其严峻。在我国已有学者模拟了宁波-舟山内海域的潮流,预测分析了含油污染物、COD、无机氮、活性磷酸盐等随海流的输运特征[3-5],但对宁波-舟山海域突发性危化品泄漏扩散模拟的研究尚不多见。本研究建立了该海域的溶解型危化品泄漏扩散模型、分析了危化品泄漏后24 h内的扩散行为,为危化品泄漏后的应急措施提供科学依据,具有较大的学术价值和现实意义。
1 宁波-舟山海域水动力数值模拟 1.1 模型设置FVCOM(Finite Volume Community Ocean Model)模型由CHEN等[6-7]发展起来,是三角网格、自由表面、原始方程、有限体积的海岸大洋数值模型,主要包括水质模块、生态模块、泥沙输运模块、流场-波浪-泥沙耦合模块等等。该模型结合了有限元法和有限差分法的优点,适合模拟浅海复杂边界[8]。模型水平方向上采用非结构化的三角形网格,可对地形复杂的区域局部加密,以更好的模拟复杂的岸线;垂向采用σ坐标,有助于处理不规则的海底地形。在潮间带,模型采用干湿网格技术考虑了潮滩对潮流的影响。因此,利用FVCOM模型非常有利于模拟具有复杂岸线和地形的宁波-舟山海域[9]。FVCOM模式在模拟近岸浅海复杂边界以及数值处理方法上的优势,已在国外多个河口、港湾以及中国渤海、长江口等多处海域得到验证[10-14]。此次研究区域为宁波-舟山海域,计算范围包括了长江口、杭州湾以及邻近海域。地形资料采用中国人民解放军海军司令部航海保证部2013年东海海图数字化结果,计算基准面为1985国家高程基准,深度为理论最低潮面。模型采用SMS软件生成三角网格(图 1),并对宁波-舟山海域进行了加密,共有50 098个网格点,94 345个三角单元,最小分辨率为200 m。开边界共有104个网格点,由潮位驱动,本次计算考虑M2、S2、K1、O1、N2、K2、P1以及Q1八个分潮。上游径流数据采用长江水文网(http://www.cjh.com.cn/)发布的大通水文站资料。风场数据从欧洲中期天气预报中心(ECMWF)官网(https://www.ecmwf.int/)下载获得,周期为6 h,该风场数据可靠性已经被众多学者所验证。
模型采用sigma坐标,垂向分为20个σ层,计算所用的初始场由东海海图图集上的月平均数据进行数字化并插值得到,为了更好得模拟近岸垂向混合,采用Mellor-Yamada2.5阶湍流闭合模型,计算时间为2016年8月1日到2016年8月31日,计算结果每10 min输出1次,计算外模态时间步长为1 s,内模态时间步长10 s,模型采用冷启动,设定初始为0的水位和流速。
1.2 模型率定和验证潮位数据采用北仑港站(29°54′N,122°07′E)、镇海站(29°59′N,121°44′E)2016年8月26日至8月28日(48 h)数据,验证结果如图 2所示,其中实线为计算值,黑点为观测值。从结果来看,潮位模拟结果吻合良好,趋势一致,误差均小于10%,可以认为模拟比较准确。
海流数据采用宁波中心站2016年在A~D 4个站点(A:30°01′24″N,121°46′13″E;B:29°59′00″N,121°46′38″E;C:29°58′30″N,122°11′5″E;D:29°58′17″N,122°12′56″E)测得的数据。验证时间为2016年8月26日至8月28日,将实测流速和流向的资料与计算值作对比,对比结果如图 4所示,通过验证,计算结果与实测结果吻合较好,流速大小、趋势符合实际情况,流向误差基本在15°以内,可以认为计算结果精确可靠,可用于该区域流场的模拟,满足危化品污染扩散计算的需求。
在对潮位、海流流速、流向充分验证的基础上,对宁波-舟山海域的潮流进行了模拟计算,其特征如下:宁波-舟山海域水流基本为往复流,涨潮时,水流从穿山半岛东侧,舟山岛南侧进入,沿着册子水道到达金塘岛,被金塘岛和册子岛分为3股具有较大流速的水流,位于舟山岛、册子岛、金塘岛以及宁波北仑区之间,大潮涨急时刻其流速可达2.5~3.0 m/s(图 5)。落潮时,水流由宁波-舟山西北侧灰鱉洋海域进入,到达金塘岛和册子岛同样被分为3股水流,并在大榭岛北面册子水道交汇,速度可达2.0~2.5 m/s,方向与涨潮时刻相反(图 6)。从整体上看,两个时刻都能明显看出册子水道上形成了一股较强水流,其流速比镇海码头大0.5~1.5 m/s。大潮期涨憩时刻和落憩时刻流速均在0.5 m/s左右。涨憩时刻流速略小于落憩时刻流速。
小潮期间流速明显小于大潮期,最大流速只有2 m/s左右。篇幅限制,本文只给出了大潮期间涨急、落急的流场图。
2 宁波-舟山海域溶解型危化品扩散模拟 2.1 研究区域概况宁波主要的化工品和油品生产企业、企业仓储、码头和输送管道主要集中在镇海、北仑和大榭3个区域,大型港口航道等运输路线也集中在宁波舟山港,因此,风险源调查主要选取镇海、北仑和大榭3个区域。根据统计镇海区共调查到主要的化工品码头企业4家,4个码头,11个泊位,码头的化工品和油品的年接卸量有4 991 548.43万吨,化工品和油品的品种主要有64种。北仑区共调查到主要的化工品和油品码头企业7家,7个码头,20个泊位,码头化工品和油品的年接卸量有4 544.87万吨,化工品和油品的品种主要有29种。大榭开发区共调查到主要的化工品和油品码头企业11家,14个码头,16个泊位,码头化工品和油品的年接卸量有126 723.92万吨,化工品和油品的品种主要有30种[15]。
由于条件限制,无法对所有溶解型危化品进行模拟,因此此次研究选取了烷烃中的环氧氯丙烷作为溶解型化学品。环氧氯丙烷为无色液体,其在25℃环境下的理化性质如表 1所示,密度略大于水,中等毒性,能溶于水,其液态在常温下蒸发性较强[16],环氧氯丙烷辛醇-水分配系数小,在水和悬浮物间分配过程不明显,不易沉积在底质上。《海水水质标准》中无环氧氯丙烷标准,在地表水环境质量标准中环氧氯丙烷限值为0.02 mg/L,取该限值作为环氧氯丙烷阈值,作为安全浓度。
此次溶解型危化品扩散模拟采用FVCOM的DYE(染色追踪模块)模块计算,该模块可用于溶于水的保守型危化品扩散模拟。危化品泄漏地点选取在较具代表性的镇海码头以及册子水道,泄漏层为表层,泄漏持续时间10 min,具体位置见图 7。
由于每个时刻潮流情况各不相同,因此本次研究选取影响范围最大情况,将泄漏时间设置在大潮期涨潮时刻以及大潮期落潮时刻,其他情况影响范围均小于此。源强采用《上海港船舶污染事故风险评估研究》给出的界定,把海上船舶污染事故分为4级,即一般(Ⅳ级)、较大(Ⅲ级)、重大(Ⅱ级)、特别重大(Ⅰ级),具体污染物溢出量的界定如下:一般事故(Ⅳ级):污染物溢出量在100 t以下;较大事故(Ⅲ级):污染物溢出量在100 t以上500 t以下;重大事故(Ⅱ级):污染物溢出量在500 t以上1 000 t以下;特别重大事故(Ⅰ级):污染物溢出量在1 000 t以上[17]。
从分类来看,100 t和500 t正好是一般事故、较大事故和重大事故溢出量的分界线,可见该源强大小的选择具有合理性;此外1 000 t为特别重大事故分界线,因此也选择1 000 t严重情况作为本研究中泄漏事故源强之一。
根据上述的泄漏时间、泄漏地点,泄漏源3个因素,制定具体计算方案如下表 2所示。
按照上述方案组进行模拟计算,由于表层以下的环氧氯丙烷浓度要低于表层,因此本文只给出表层扩散情况并加以分析。环氧氯丙烷泄漏24 h后得到以下结果,如表 3所示。由表可知,100 t环氧氯丙烷泄漏24 h后,海域中已经不存在超安全浓度(0.02 ml/L)的区域;500 t环氧氯丙烷泄漏24 h后,除涨潮时刻泄漏的情况依旧存在10 km2左右的超安全浓度水域,其他情况环氧氯丙烷也已扩散至0.02 mg/L以下;1 000 t泄漏24 h后,所有泄漏方案海域中均存在较大面积的超安全浓度水域;为研究影响相对较大的情况,选取泄漏量为1 000 t作为典型,其余计算方案影响范围均小于此。1 000 t泄漏后海域环氧异丙烷最高浓度下降趋势如图 8所示,扩散情况如图 9至12所示,图中包裹着扩散区域的实线为0.02 mg/L包络线。
当泄漏事故发生在镇海码头,泄漏时间为涨潮时刻,前6 h,环氧氯丙烷以镇海码头为中心均匀向四周扩散,且0.02 mg/L包络线包裹区域内浓度呈现较高水平,随着扩散稀释,包络线逐渐扩大,但在大榭岛北仑区之间依旧存在污染浓度较高的水域,如图中的12 h和18 h时刻。泄漏24 h后,环氧氯丙烷已遍布北仑区至舟山岛之间,0.02 mg/L包络线出现分离,分别位于北仑区北侧以及穿山半岛东侧,总面积达到190 km2左右(图 9)。落潮时刻泄漏扩散较涨潮时刻快2~6 h(图 8),泄漏6 h内,环氧氯丙烷沿着镇海沿岸迅速移动到大榭岛,在大榭岛与北仑区之间形成较高浓度水域,随后随着涨落潮水流来回移动,0.02 mg/L包络线向四周展开,到18 h左右,包络线出现分离的情况,分为3块,分别位于北仑沿岸一带,大猫岛东侧以及北仑区南侧,24 h后,0.02 mg/L包络线包裹着大榭岛周围海域,面积为140 km2左右(图 10)。
图 11和图 12分别是册子水道上涨落潮时刻发生泄漏的扩散情况。涨潮时刻,泄漏的环氧氯丙烷立即沿着金塘岛和册子岛之间的大流速水流扩散到金塘岛西北侧灰鱉样水域,并随着涨落潮的水流,沿着册子水道来回移动,最远到达30.03°N左右,距泄漏点20 km左右处,在泄漏24 h后,0.02 mg/L包络线位于金塘岛、册子岛以及舟山岛之间,覆盖面积达100 km2左右。与镇海码头泄漏相同,落潮时刻泄漏明显快于涨潮时刻(图 8),当泄漏发生时,环氧氯丙烷沿着册子水道的水流,迅速向南迁移至大榭岛东侧,随着污染继续扩散,0.02 mg/L包络线发生变形、分离,在18 h左右,包络线出现在两处,分别位于金塘岛与镇海区之间,以及大榭岛东侧,扩散24 h后,除了在金塘岛西北侧存在一块约30 km2的超安全浓度水域,其他水域内环氧氯丙烷浓度基本在0.02 mg/L以下。
此外,对比两个不同的泄漏地点,册子水道由于具有比镇海码头大0.5~1.5 m/s的流速,环氧氯丙烷泄漏扩散速度快于镇海码头,如图 8所示,扩散12 h内,表示册子水道泄漏的两条实线下降速度明显快于镇海码头,12 h后,4种泄漏情况最高污染浓度的下降趋势均趋于稳定。
3 结论对宁波-舟山海域的水动力情况进行了模拟。结果与实测资料对比,误差较小,模拟比较准确。在充分率定检验的基础上,选取环氧氯丙烷为典型,利用FVCOM数值模型对宁波-舟山海域溶解型危化品泄漏行为进行了模拟。
模拟结果表明:从不同事故等级来看,一般事故(100 t)和较大事故(500 t)危化品泄漏24 h后,环氧氯丙烷浓度基本扩散至安全浓度(0.02 mg/L)以下,重大事故泄漏24 h后仍存在较大(30~190 km2)超安全浓度污染水域;从不同扩散潮时来看,无论在哪个地点泄漏,落潮时刻危化品扩散均快于涨潮时刻,落潮时刻扩散至安全浓度以下较涨潮时刻快2~6 h;从不同泄漏点来看,册子水道由于具有比镇海码头大0.5~1.5 m/s的流速,扩散明显快于镇海码头,其中在册子水道泄漏时,危化品基本沿水道扩散,在镇海码头泄漏时,危化品沿着宁波沿岸扩散,且易在北仑区沿岸、北仑和大榭岛之间滞留。
以上,对于宁波-舟山海域溶解型危化品泄漏事故,本文建立的该海域溶解型危化品泄漏扩散模拟能为危化品泄漏后的应急措施提供科学依据, 应急监测的单位可根据不同情况的预测模拟迅速作出响应,在泄露最初的一个潮周期内采取相应措施。
[1] |
虞志娟, 沈家. 基于GIS技术的危险化学品物流系统的风险管理[J]. 物流技术, 2008, 27(4): 55-58. YU Z J, SHEN J. Risk management of the dangerous chemicals logistics system based in GIS technology[J]. Logistics Technology, 2008, 27(4): 55-58. |
[2] |
胡艳菊, 李彦海, 胡忆沩. 危险化学品突发泄漏事故应急决策系统的开发与应用[J]. 中国安全生产科学技术, 2010, 6(4): 109-113. HU Y J, LI Y H, HU Y W. Development and application of the emergency decision system for dangerous chemicals burst leaking accident[J]. Journal of Safety Science and Technology, 2010, 6(04): 109-113. |
[3] |
耿兆铨, 倪勇强, 程杭平. 溢油和含油污水对舟山海域水质影响的数值分析[J]. 水动力学研究与进展(A辑), 1991, 6: 46-54. GENG Z Q, NI Y Q, CHENG H P. Numerical Study of Oil Pollution on Sea Region[J]. Journal of Hydrodynamics, 1991, 6: 46-54. |
[4] |
曹欣中, 唐龙妹. 宁波、舟山内海域实测海流分析及潮流场的数值模拟[J]. 东海海洋, 1996, 41(2): 1-9. CAO X Z, TANG L M. The analysis of the observed sea current and the numerical model of the tidal current field for the inner sea area if Ningbo and Zhoushan[J]. Donghai Marine Science, 1996, 41(2): 1-9. |
[5] |
蔡惠文, 孙英兰, 张越美, 等. 宁波-舟山海域污染物扩散的数值模拟[J]. 中国海洋大学学报(自然科学版), 2006(6): 975-980. CAI H W, SUN Y L, ZHANG Y M, et al. Numerical Simulation of Pollutant Diffusion in the Ningbo-Zhoushan Sea Area[J]. Periodical of Ocean University of China, 2006(6): 975-980. |
[6] |
CHEN C, LIU H, BEARDSLEY R C. An unstructured grid, finite-volume, three-dimensional, primitive equationsocean model:Application to coastal ocean and estuaries[J]. Journal of Atmospheric and Oceanic Technology, 2003, 20(1): 159-186. DOI:10.1175/1520-0426(2003)020<0159:AUGFVT>2.0.CO;2 |
[7] |
CHEN C, BEARDSLEY R C, COWLES G. An unstructured grid, finite-volume coastal ocean model (FVCOM)system[J]. Oceanography, 2006, 19(1): 78-89. DOI:10.5670/oceanog |
[8] |
XIU G W, HAN S T. Coupling of CFD model and FVCOM to predict small-scale coastal flows[J]. Journal of Hydrodynamics, 2010, 22(5): 284-289. DOI:10.1016/S1001-6058(09)60208-0 |
[9] |
韩松林, 梁书秀, 孙昭晨. 基于FVCOM的象山港海域潮汐潮流与温盐结构特征数值模拟[J]. 水道港口, 2014, 35(5): 481-488. HAN S L, LIANG S X, SUN Z C. Numerical simulation of tides, tidal currents and temperature-salinity structures in Xiangshan Bay based on FVCOM[J]. Journal of Waterway and Harbor, 2014, 35(5): 481-488. |
[10] |
YU Z, CHEN C S, ROBERT C, et al. Studies of the Canadian Arctic Archipelago water transport and its relationship to basin-local forcings:Results from AO-FVCOM[J]. Journal of Geophysical Research:Oceans, 2016, 121(6): 4392-4415. DOI:10.1002/2016JC011634 |
[11] |
MILES C, JENNIFER A, YERUBANDI R. Sensitivity of the Simulated Kingston Basin-Lake Ontario Summer Temperature Profile using FVCOM[J]. Atmosphere-Ocean, 2013, 51(3): 319-331. DOI:10.1080/07055900.2013.800017 |
[12] |
XUE P F, CHEN C S, QI J H, et al. Mechanism studies of seasonal variability of dissolved oxygen in Mass Bay:A multi-scaleFVCOM/UG-RCA application[J]. Journal of Marine Systems, 2014, 131: 102-119. DOI:10.1016/j.jmarsys.2013.12.002 |
[13] |
LI B, KISEI R, CHEN Y, et al. Assessing the quality of bottom water temperatures from the Finite-Volume Community Ocean Model (FVCOM) in the Northwest Atlantic Shelf region[J]. Journal of Marine Systems, 2017, 173: 21-30. DOI:10.1016/j.jmarsys.2017.04.001 |
[14] |
WAN Z H, REN Z F, WANG X C, et al. Application of unstructured grid finite volume coastal ocean model(FVCOM)to the Yangtze River hypoxic zone[J]. International Journal of Numerical Methods for Heat & Fluid Flow, 2016, 26(8): 2410-2418. |
[15] |
杨耀芳, 叶然, 周巳颖, 等. 宁波海域环境风险源识别与分级研究[J]. 海洋开发与管理, 2014, 31(7): 119-124. YANG Y F, YE R, ZHOU S Y, et al. Ningbo Marine environment risk source recognition and classification research[J]. Ocean Development and Management, 2014, 31(7): 119-124. |
[16] |
北京化工研究院. 化学品毒性、法规、环境数据手册[M]. 北京: 中国环境科学出版社, 1992: 59-62. BEIJING RESEARCH INSTITUTE OF CHEMICAL INDUSTRY. Handbook of Chemical Toxicity, Legislation, Environmental Data[M]. BEIJING: China Environmental Science Press, 1992: 59-62. |
[17] |
陈维. 上海港船舶溢漏油污染事故风险评估研究[J]. 航海技术, 2010(6): 75-77. CHEN W. Study on Risk Assessment of Ship Spill Oil Pollution Accident in Shanghai Port[J]. Navigation Technology, 2010(6): 75-77. |
2. East China Sea Branch Of State Oceanic Administration, Shanghai 200137, China;
3. East China Sea Environmental Monitoring Center, S. O. A., Shanghai 201206, China