基于MIKE21FM北京市南水北调配套工程.docVIP

基于MIKE21FM北京市南水北调配套工程 　　摘要：大宁水库库尾位于交通节点附近，存在交通事故引发突发性水污染事故的风险。从风险防范与管理角度，利用MIKE21FM模型，以5 t硝基苯泄露为例，考虑风、水库调度因素，设定5个情景开展突发事故污染物扩散规律研究。结果表明：（1）不同情景污染物扩散速度和超标水域面积相差较大。无风、不调度情景下，库区水体接近静止，污染物扩散速度最慢；水库向大宁调压池输水20 m3/s，西北风12.3 m/s情景下，污染物扩散速度最快。（2）与常规调度相比，风速和风向对污染物的扩散影响更大。同样风速情况下，西北风较东北风更利于污染物的扩散。（3）正常情况下（即风速不超过2.56 m/s，调水不超过20 m3/s），5 t硝基苯泄露应急响应时间在15 h以上，最不利情况下（风速为20年最大平均风速12.3 m/s且调水20 m3/s），应急响应时间为4～9 h。根据研究结果，一方面建议加强库尾交通疏导，加固桥梁护栏，降低交通事故发生概率；另一方面，建议在现有监测站点规划基础上，新增库尾自动监测站一处。 　　关键词：北京市南水北调配套工程；大宁水库；突发水污染事故模拟；MIKE21FM 　　中图分类号：X52；TV68 文献标识码：A 文章编号：1672-1683（2013）04-0067-05 　　大宁调蓄水库工程作为北京市南水北调配套工程的重要组成部分，在调蓄来水、工程检修、事故应急供水方面发挥着至关重要的作用[1]。由于水库库尾（库区北端）位于交通节点附近，车流量大，交通事故导致化学品泄漏引发水污染事故的风险相对较大。在当前突发性水污染事故频发的背景下[2]，考虑到南水北调水质安全的重要性及研究区地理位置的特殊性，开展突发性水污染事故应急等相关方面的研究是十分必要的。 　　在水污染事故应急过程中发生时，管理部门迫切需要了解污染带的迁移状况以及污染物浓度在时间、空间两个维度上的变化，以便采取有效的应急措施[3]。当前对于此问题的研究，受到区域特点、模型简化、数据限制等因素影响，多集中于河流水污染事故的模拟，大多采用一维模型[45]。但一维模型只适用于河宽较窄的中小型河流，当研究区为较大河流或湖库时，一维模型已不能满足研究需要。MIKE21FM是目前应用较为广泛的二维模型，在平面二维自由表面流数值模拟方面具有强大的功能。本文基于MIKE21FM开展大宁水库水污染事故的模拟分析，旨在探明污染物的时空变化规律，明确污染影响范围和程度，为后续水库风险应急管理提供依据。 　　1 水库概况 　　大宁水库位于北京市房山区大宁村北，京港澳高速公路（京石高速）杜家坎至赵辛店路段东侧，与南部的稻田水库、马厂水库组成永定河滞洪水库，见图1。上游为小哑叭河，下游为小清河。南水北调总干渠永定河倒虹吸由大宁水库库区穿过，大宁调压池位于大宁水库副坝南侧，通过退水暗涵与之相连。改建后库区总库容为4 611万m3，防洪库容3 611万m3，正常蓄水位为56.40 m。当水位低于52.46 m，为向总干渠供水，需启用位于水库副坝北侧130 m处的抽水泵站。 　　大宁水库库尾正是京石高速、G107等公路线及老京广正线、二七厂等铁路线交汇处，距杜家坎交通节点直线距离不足1 km。库尾处有公路桥梁1座，铁路桥梁2座，交通繁忙，车流量较大；水库坝前为大宁调压池，无正式公路，基本无外来车辆。因此在库尾发生交通事故并引发突发污染事故的可能性远大于坝前。 　　2 水质模型构建 　　2.1 模型选择与概述 　　大宁水库在56 m正常蓄水位时平均水深不足10 m，不存在明显分层现象，采用平面二维模型可以满足研究需要。研究区边界不规则，采用传统矩形网格模拟可能会影响模拟结果。经模型比较，选择使用MIKE21软件包中的MIKE21FM模型。MIKE21由丹麦水力研究所开发，在平面二维自由表面流数值模拟方面具有强大的功能，具有用户界面友好、前后处理功能强大、数值计算方法科学等诸多优点[6]，近年来在国内得以广泛应用。其研究对象涵盖海洋、河口、河流、湖库等水体，研究类型包括单纯的数值模拟[78]、洪水淹没分析[910]、温排水影响分析[1112]、水质模拟与预测[1314]等几类。但在突发事故模拟方面，相关应用并不多见。 　　2.2 模型地形文件构建 　　以分辨率为30 m的大宁水库DEM为数据源，在ArcGIS平台的支持下进行剪裁后导出成xyz格式；参照《北京市南水北调配套工程大宁调蓄水库工程初步设计报告》（以下简称“设计报告”），结合现场勘察结果，对地形数据进行初步修正；利用MIKE Zero中的Mesh Generator工具建立水库三角网格，导入修正后的xyz数据，形成模拟地形，见图2。模型建立后，对不同水位下库容进行模拟，