Excel VBA在分级控制水库调洪计算中的运用.pdfVIP

第28卷第2期 水电能源科学 V01．28No．2 201 Water 0年2月 ResourcesandPower Feb．2010 文章编号：1000．7709(2010)02—0038—03 Excel VBA在分级控制水库调洪计算中的运用 谭毅源 邱 钢 (楚雄欣源水利电力勘察设计有限责任公司，云南楚雄675000) 摘要：针对分级控制水库的泄洪特点，使用ExcelVBA编制了调洪计算数值解程序，计算速度快、精度高、结 果直观，应用于下口坝水库工程中，取得了较好的工程效益。 关键词：ExcelVBA；水库；分级控制；调洪计算；程序 中图分类号：TVl22；TV697．1 文献标志码：B 传统的水库调洪计算方法有试算法、图解法 式中，Z(t)为水位；Q(f)为人库流量；口(Z)为出库 或相关图法，试算前需绘制水位一下泄流量和水 流量；，(名)为水库水面面积，是水位的函数；Zo为 位一库容关系曲线，试算过程中需反复查读曲线 和计算，工作量大，且易引起误差和错误。分级控 制调洪计算的试算工作量则又比常规调洪计算大 的初值问题，采用定步长的龙格一库塔方法啪，其 数倍，使用ExcelVBA编制数值方法求解的程 序，可大幅缩小计算工作量、提高计算精度。利用 Excel的界面输入输出数据，既直观又方便，且可 将计算结果链接于曲线图，实现动态观察水位流 量的变化过程；Excel自身强大的计算功能又可 K·=丁啦铲 省去许多程序代码的编写工作，易于被工程人员 所掌握。 。，个Q(￡，，+虿T)一q(乙，+争) 心钉—‘=疆■一，(Z，l+睾) 1 数值方法 1．1调洪计算原理和方法 K)譬+--．Z(q一)T虿+t，∞， 3一丁—————兰————F——三 洪水入库后的流态属于明渠非恒定流，水库 厂(乙-+未) 沿程的水位、流速和过水断面等均随时变化，明渠 非恒定流的基本方程可用水力学的圣维南方程组 Ⅳ一个Q(“)一q(Z，广l+K3) n4—1—]冠=了i厂 表示。该偏微分方程组一般难以求得精确的解析 式中，丁为洪水流量时段间隔；，l=1，2，…，_『。 解，实用中多采用瞬态法、直接差分法和特征线法 1．3泄水建筑物泄流量公式 等近似解法。其中，瞬态法最为简单实用，用有限 当泄流为孔口出流时删： 差分代替微分值并加以简化，忽略了洪水入库至 口=矽P、／29(Z—Cl—P)(3) 泄洪建筑物间的行进时间、沿程流速变化及动库 式中，岸为流量系数；b为泄流孔口宽；P为泄流孔 容等的影响，得到水量平衡方程详见文献[1]。 口高；Z为水位；c1为泄流孔口底槛高程。 1．2水量平衡分方程的数值解 当泄水建筑物为堰流时‘4]： 求数值解采用水库水量平衡的微分方程[2]： (4)