基于二维溃坝波数值模拟的生命损失估算 Life loss estimation caused by reservoir dam break based on 2D dam-break wave numerical simulation.pdfVIP

第4期 水利水运工程学报 N。．4 AND 2013年8月HYDRO-SCIENCEENGINEERING Aug．2013 基于二维溃坝波数值模拟的生命损失估算 沈照伟1，潘存鸿1，叶碎高1，王立辉2，鲁海燕1 摘要：通过界面数值通量以及底坡源项的计算求解尢碰撞二维Boltzmann方程，建立了二维溃坝波数值模型． 以杭州市余杭镇上游某水库为例，模拟了该水库溃坝后溃坝波的演进过程，模拟结果表明溃坝洪水淹没范围很 大，距坝约11km的下游，水深仍可达4．0 m以上．同时，应用通用的溃坝生命损失模型，在获取网格单元上的水 力学参数值(包括水深和最大流速)基础上，输入当地的地形、人口等基本数据，估算了水库溃坝对下游造成的 生命损失．该方法可用于其他水库的风险管理分析，为相关部门决策提供依据和参考． 关键词：溃坝波；数值模拟；生命损失 中图分类号：TVl22+．4；TV698．2+37 水库大坝溃决的风险是客观存在的，一旦溃决将对水库下游群众的生命和财产造成重大损失．因此，溃 坝生命损失评估对水库大坝安全运行及区域防洪减灾均具有重要意义． 些方法是在对众多已溃坝水库生命损失的统计和研究的基础上提出来的．BG法和DM法主要应用数理 统计回归方法，利用多个主要生命损失影响参数的经验公式来估算．由于溃坝资料有限，拟合效果可能也不 是很好．Assaf法是基于风险人口生还率的生命损失评价新模型的方法，由于目前对风险人口主观能动性、溃 坝洪水的互动性研究还不够深入，往往应用的条件还不够成熟．周克发等‘1o提出的生命损失模型更加符合 我国实际，该模型计算中必须考虑溃坝洪水演进过程的水深与流速．现有的一些洪灾损失模型主要按照一般 洪水的演进计算，而溃坝洪水与一般洪水计算相比主要有以下特点：(1)动边界处理时，在水量严格守恒基 础上，要求具有处理强间断情况下的动边界能力；(2)溃坝波水力条件变化梯度往往很大，计算时要求在不 规则地形下具有稳定性，同时具有捕捉激波的能力；(3)溃坝波计算区域很大，计算时要有局部加密功能，同 时要有较强的边界适应能力．构建边界适应性强，并具有模拟间断流动能力的二维溃坝波数值模型，是溃坝 波模拟计算的关键．早期的溃坝波数值模型大多为平底情况‘2。4J，近几年来，一些学者应用能模拟大梯度流 动的高性能格式建立了二维溃坝波数值模型，并在实际水库溃坝和河道溃堤中逐渐开始应用’5。7o．本文应用 FluxVector 有限体积一KFVS(Kinetic Riemann解方法计算动边界问题，以杭州市余杭区余杭镇为例，计算上游某水库溃坝洪水演进过程．结合现 有通用的溃坝生命损失模型，开展生命损失评估．在获取网格单元上的水力学参数值(包括水深和最大流 速)基础上，以区域地形、人口等数据为模型输人，分析了水库溃坝对下游余杭镇造成的生命损失．对水库大 坝安全运行及区域防洪减灾均具有重要的意义． 收稿日期：2013—03—05 基金项目：国家A然科学基金资助项目；浙江省重大科技计划资助项目(2007C13012) 作者简介：沈照伟(1975一)，男，浙江桐庐人，高级工程师，博士，主要从事防洪减灾及水库大坝安全评估 E-mail：shenzw@zjwater．gov．cn 万方数据 第4期 沈照伟，等：基于二维溃坝波数值模拟的生命损失估算 1研究方法 1．1溃坝波数值模型 3 无碰撞二维Bohzmann方程为‘8 dO-f。-+cx善4-Cy等+吼善+妒，丝Oy=o (，) 式中，为分子速度分布函数；乞为菇方向的分子速度；c，为在Y方向的分子速度；妒为外力作用、由非平底引 起的阻力项及重力项，织=g(S。。一％，妒，=g(s吖一勘)，g为重力加速度，％和勘为石和Y方向的阻力项， S。和S。，为石和Y方向的底坡项． 将Bohzmann方程乘以(1，c，，c，)7，并对分子速度空间积分，可以