A港区海上溢油事故数值模拟预测.docxVIP

A港区海上溢油事故数值模拟预测 摘要：为研究A港区溢油事故的影响，通过MIKE21软件建模、并综合考虑风向、风速和潮流等因素设定6种不同情景进行模拟预测。结果显示，在不利风条件下，72h油膜扩散距离为0.7～39.50km，扫海面积可达0.03～142.97km 0 引言随着海洋石油开采和海洋运输业的发展，海上溢油事故的发生频率也相应增加。石油中含有烃类等有毒有害物质，一旦泄漏，会对海洋环境造成不良影响。因此，开展海上溢油数值模拟研究，科学预测海洋溢油扩散特征，对海上溢油防范、扩散控制、应急处理及减少溢油对海洋生态环境的影响具有重大意义A港位于某开发区内，主要为开发区及周边地区的物资运输服务，是一座综合性港口和现代物流基地。一旦发生溢油事故，将对该海域生态环境乃至沿岸生产活动造成不良影响1 研究区环境概况A港水域由B湾和C湾组成，周边的海洋功能区包括海洋保护区、保留区、农渔业区、休闲娱乐区等。主要环境保护目标包括Ⅰ红树林保护区、Ⅱ自然保护区、海洋生态红线区和渔业养殖区等。其中，Ⅰ红树林保护区主要保护对象为红树林生态系统。植物群落类型以红海榄群落为主。植物种类有红海榄、木榄、秋茄、角果木、桐花树、海莲、海漆、榄李、白骨壤、桐花树、银叶树及黄槿等；Ⅱ自然保护区主要保护对象为保护区、周边近海珊瑚礁及其生态环境。2 模型的建立与验证以MIKE21软件建立研究海域的二维水动力模型，然后耦合油粒子模型进行溢油模拟计算。丹麦水环境模型集MIKE21是二维自由水面流动模拟系统工程软件包，适用于湖泊、河口、海湾和海岸地区的水力及水质的平面二维仿真模拟。2.1 计算域的确定及网格划分本文建立的海域数学模型采用非结构三角网格进行模拟，整个模拟区域由6 629个节点和12 029个三角单元组成，在工程附近进行加密，加密后网格空间步长约20～50m。计算域模型网格及水深地形见图1。2.2 参数设置模型计算时间步长根据CFL条件进行动态调整，确保模型计算稳定进行，最小时间步长0.05s。底床糙率通过曼宁系数进行控制，曼宁系数M取20～60m2.3 边界条件与初始条件1）开边界条件：即水域边界条件。在此边界上，或者给定流速，或者给定潮位。本文开边界给定潮位，即：η=η（x,y,t)2）闭边界条件：即水陆交界条件。在该边界上，水质点的法向流速为0，即：V3）初始条件本模型初始条件采用稳态启动，初值的误差将随着计算的进行而逐步消失，文中模型验证取计算稳定后的模拟结果开边界的选择采取潮位驱动，潮位数据由TPX-OChinaSeas潮波预报模型获得，该模型是2010年美国俄勒冈州立大学在TPXO全球潮波预报模型的基础上建立的、针对中国海域的高精度潮波预报模型，共包含8个分潮，其网格精度为1/30°×1/30°2.4 模型验证潮位验证数据采用L3站位和B1临时验潮站两个站位的实测数据，监测时间分别为2014年10月14日—15日和2018年7月27—28日。潮流验证站点包括C1～C6共6个站位，此6个站位的潮流验证数据为2018年7月27—28日的连续观测数据。潮位验证结果见图2，潮流流向和流速验证见图3。验证结果表明各站位对应观测点上潮位和潮流模拟结果与实测潮位和潮流资料基本吻合，时间偏差和平均流速偏差符合技术规程。预测流场能够较好地反映周边海域潮流状况，可以为后续溢油模型提供较可靠的动力场。2.5“油粒子”模型介绍溢油进入海洋水体后，在自身重力及海洋水体物理化学的作用下，同时发生扩展、漂移、扩散、蒸发、乳化、溶解等风化过程3 模拟预测结果与分析3.1 预测模式中有关参数的选取船舶溢油事故多发地点主要集中在码头、航道等近岸海域，因此溢油位置选择在A港码头前沿。A港多靠泊集装箱船，来往船只较多，难免发生碰撞溢油事故，本次考虑典型船舶的储油舱单舱全部破裂溢油，溢油量为约380m3.2 模拟情景模拟情景的选取综合考虑风向、风速和潮流因素。港区附近敏感目标主要为Ⅰ红树林保护区、Ⅱ自然保护区以及海洋生态红线区，因此不利风向是W向及NE向，风速考虑年平均风速和最大风速。模拟情景组合见表1。3.3 模拟预测结果分析不同条件下，溢油预测结果见图4，溢油影响数值分析结果见表2。1）不同的风向风速和潮汐情况下，漂移距离、扫海面积与残存油量不同。均风条件下，油膜扩散距离9.85～27.05km，扫海面积10.05～130.57km2）不利风向情景下，溢油很快（0.5h）抵岸。溢油抵达海洋生态红线区的最短时间为5.5h，抵达Ⅱ自然保护区的最短时间为7h，抵达Ⅰ红树林保护区的最短时间为22.5h。溢油发生72h后，残油量仍高达60%以上。若无有效应急措施，溢油将会对