岸堤水库洪水预报及调洪演算软件使用说明书

1、岸堤水库雨洪资源解析 使用说明书二0一五年六月一日作者:文华:*:fblwh150163.目录第一章概述（3第二章功能简介（5第一节功能特点（5第二节软件画面（6第三节运算功能（7第四节气象云图及气象雷达（13第三章数学模型（14第一节洪水模型（141、瞬时单位线（142、CAMMADIST 函数语法（153、CAMMADIST 函数应用（164、流域洪水错时叠加（17第二节洪水传播（18第三节泄量模型（191、闸门出流（192、推求水面线（213、闸门泄量（22第四节调洪演算（22第五节控运案（23第四章扩展性设计（23第五章调洪实例（29第六章课目攻关概况（30第七章使用说明书（31第一节

2、洪水预报（31第二节调洪演算 （33第三节其他计算（33附件课题研发小组成员错误!未定义书签。第一章概述控制和预见洪水，让洪水变为一种资源,实现科学预见、动态管理、合理利用，是 本课题的研究对象。科学控制洪水，真正能够对洪水运用自如，其首要问题是准确解析、及时预报，掌 握洪水动态。但目前实际应用中，对水库防洪兴利控制运用，还仅限于依靠库水位的 变化，结合下游河道的承受能力，试探性的调节洪水，这种洪水调整模式，具有较大的盲 目性，理论面的支撑相对不足。当前,各水库防汛主体单位，均制定了相应的水库控制运用案。如岸堤水库 防洪调度图（图1,但这些案的编制和批复仅表现为粗线条和原则性的界定，是在进行

3、大量假定的基础上进行编制的,应用中的可操作性相对欠缺，在实践中仅具有指导意 义。（图1洪水调度控制案的编制,偏离实际应用，存在的突出问题，主要表现在以下几个面:1、假定了降雨的空间分配是均匀的，即整个流域降雨分布是均等的。但实际降 雨,特别是流域面积稍大的水库，降雨的空间分布几乎不可能是均等。2、事先拟定了24小时降雨在1日各时段上的雨量分配。但实际降雨在时段 上的分配，是个随机的不确定因素。3、控制运用案的编制起调水位为汛中限制水位，但实际降雨前的库水位，却几 乎不可能恰巧是汛中限制水位。4、所有闸门同开度启用，与实际控制运用也不相符。5、控运案的编制基础是净雨量，但实际操作中我们所能掌握的

4、是降雨量。把降 雨量合理的转换为净雨量，也成了水库实际控制运用要解决的问题。岸堤水库雨洪资源动态解析预报平台,较为圆满的解决了以上五面的问题。第二章功能简介第一节功能特点1、采用.dll 数据调用模式，获取水库基本数据。软件具有非常理想的可扩展性, 支持任某水库的洪水资源解析。2、流域洪水模型参数,支持外部修订（需慎重修订;上游中小型水库、塘坝等的 拦蓄作用，可通过修订流域权重参数，简化上游拦蓄效果的修正计算。3、数据支持电子表格输出和图形输出。可导出对应某次降雨的洪水过程线文 件、对应某泄洪案的库水位动态变化过程线文件、调洪演算成果文件及其他数据的 导出文件等。4、支持具有指导意义的控制运用

5、案的编制计算。v2.0及以后版本，还支持演 算过程数据文件（*.xls的导出。5、支持库水位库容水面积对应关系转换（附:溢流堰泄量;支持库水位闸 门开度溢洪道泄量对应关系转换;支持降雨信息、降雨重现期等其他数据转换，支 持逆推进库流量等计算功能。6、界面采用GDI +函数绘制;强化鼠标移动算法,瞬间获得流域、计算钮等虚拟 按键操作指令;频繁调用的关键代码，特别是v2.20及以后版本,大量的嵌入了汇编语言代码，直接操作 计算机CPU和存等硬件设备,数据处理快速、高效。7、全面支持包括调洪规则等大量的信息自定义。第二节软件画面1、启动画面系透明.png文件（图2,其背景系洪水创意手绘图。（图22、

6、主画面背景系岸堤水库流域布局图（图3（图3第三节运算功能1、前期降雨情况数据输入功能（图4前期降雨影响输入模式下，还可通过右键菜单，调用软件置的，前期降雨影响量 计算器，进行前期降雨影响计算;支持快速自动录入Pa值。（图42、流域降雨数据，按真实时程分配输入（图5（图53、洪水过程线屏幕输出（图6,支持数据文件导出（*.xls和*.png。（图64、调洪库水位过程线屏幕输出（图7,支持数据文件导出;支持任意起调水位和 任意闸门操作形式下的洪水控制运用计算。（图74、控制运用案编制计算模块（图8。支持任意起调水位（默认为限制水位设置； 支持外置调洪规则设置。软件作者配置的调洪规则意义为:起调水位

7、以下不泄洪，防 洪高水位以下，按下游河道安全泄量调洪（保下游，库水位超过防洪高水位后，敞闸泄洪 （保工程;并依据综合利用水库调度通则（水利部水管199361号对防洪调度图的 编制要求,对降雨做如下假定：假定1:流域均匀净雨量;假定2:流域净雨的时段分配，恒定匹配于水库所在地域；假定3:所有闸门同高度启用。软件V2.0及以后版本的调洪演算模块，可导入自定义洪水;可导出调洪数据文件, 演算过程清晰直观。（图85、库水位库容水面积，对应关系表（图9。（图9 6、闸门泄量计算功能区（图10（图107、进库流量逆推功能（图11（图11 8、降雨量重现期（图12（图129、软件的v2.0及以后版本,还设计

8、了，脚本调试器。调试器直接运行脚本文件 VbThguize.txt，并输出计算结果，软件使用者可验证自己编写的脚本文件是否正确。10、软件v2.0及以后版本，前期降雨影响的计算是以陆地蒸发系数的模式予以 外部设置，文件名是在配置文件中予以指定YueyinxiangPa.txt的形式外置，允用户根 据当地的气候条件，参照已有资料随时予以调整和更新。11、外置调洪规则脚本（VbThguize.txt，功能强大，对于调洪演算和洪水预报，可 轻松的实现闸门的动态操作（详见:第四章扩展性设计。12、v2.20及以后版本,支持洪水时段长自定义。13、v2.24版本还支持小流域洪水数据导出。14、v3.00

9、版本的洪水河道传播，采用了被国外广泛应用的马斯京根算法。第四节气象云图及气象雷达1、，风云二号气象卫星实况数据调用（图132、，当地、当前气象雷达数据调用（图14（图13（图14第三章数学模型第一节洪水模型1、瞬时单位线为确保软件解析功能在各水库之间的通用性和可拓展性，洪水单位线采用传统 的瞬时单位线。瞬时单位线的理论依据就是把洪水过程看作是流域n个相同的线性水库串联 调节后的汇流结果,是个已经过理论论证和推导，并被大家所接受的成熟结论。其数学表达式：11u(t=(n t kt ek n k -/)I其中:n二线性水库数量,反应流域调节能力；(n r= n的伽玛函数；K二线性水库的调节系数。e

10、二常数(自然对数的底数瞬时单位线S曲线数学表达式：101S(t=(n t ktkt ed t k n k -( lr打常规利于瞬时单位线推求洪水过程线时,我们就是通过S曲线，计算无因次时段 单位线，最后转换为有因次时段单位线。这种汇流计算,不仅是工作量大，而且通过 S(t/k ,n推求参数的时段单位线时,一面因使用中间参数(t/k增加了舍入误差,另一面 (t/k和n的双向插，增加了法误差。计算结果的精度受到了很大影响。微软Excel办公组件为我们提供了 CAMMADIST函数是洪水计算较为便捷理 想的计算案。2、CAMMADIST函数语法CAMMADIST(x,alpha,beta,cumul

11、ative式中:x是自变量，用来计算伽玛分布函数 的值；Alpha和beta均为相关参数。Cumulative是个逻辑参数：Cumulative =false时,返回概率密度函数：xf x x e apaappa-=rCumulative =true时,返回累积分布函数:11F(,(xxx x ed x apaappa-= rj2C 1212M M FJ RT 0.196F km J m m R mm T h M 3M M 3aa- - - - 一K省瞬时单位线参数经验公式：二山丘区二流域面积(河道干流平均坡度(/净雨量(净雨历时(若:则 3 若:则 312M n k 1M n kM n单位线

12、参数与、的关系：3、CAMMADIST函数应用洪预报解析平台,就是调用微软为我们提供的CAMMADIST函数，把,x a参数， 分别用瞬时单位线的,t n k 代替，当 Cumulative =false 时得到 u(t 函数；当Cumulative =true时得到S(t函数。因此洪预报解析平台要求用户必须安装 有微软的Excel办公组件,这是个大多数计算机用户都可满足的一般条件。(1(,CAMMADIST(,1CAMMADIST(,110q(,(,3.610Q(CAMMADIST(1,1CAMMADIST(, ,13.6mi i u t t t n k t t n k Ft t u t t

13、 tF t h t i n k t i t n k t -= - 二 二.- + .-.学无因次时段单位线:有因次时段单位线:流域汇流计算公式:4、流域洪水错时叠加CAMMADIST函数本身并不能解决降雨在空间上的不均匀分布问题。透过微分 观点,我们可以把全流域进行细分，分割成若干个微小流域，我们完全可以认为微小流 域的降雨是均布的。鉴于目前技术条件和设备配置，我们还无法得到无限微小流域 的降雨量，因此，流域的细分，不能真正的进行无限分割，必须结合雨量观测设备在全流 域的真实布局情况予以分割。利于微软的CAMMADIST函数计算洪水过程，流域划分越细计算结果越合理。 按照水库全流域雨量检测设备

14、的分布，进行了二次分割，把整个水库流域,细分成了多 个小流域，依次计算洪水并经过河道传播和叠加计算,最终得到流域坝前洪水过程 线。岸堤水库上游在水库整个流域,分布着4个水文站和5个雨量观测点，按照这些 雨量观测设备的布局，把整个流域又分割成了 9个小流域。这些小流域的洪水，通过 调用微软的CAMMADIST函数得到相应过程线，各个小流域洪水过程线经主河道传 播后,错时段叠加得到洪水演算所需要的坝前洪水动态数据。第二节洪水传播流域洪水在河道中的传播,采用马斯京根算法。马斯京根法假定,对某特定河道 段,其河道槽蓄量S与流量Q之间，存在着一种固有关系,即槽蓄曲线。槽蓄曲线反映河段的水力学特性。（S

15、Q K 1Q 1K x K Q Q x Q x Q x Q x f Q - - +- +- fff下上下上下上二，二（式 中:=（示储流量二流量比重因子河段汇流时间（物理意义建立如下的河道水量平衡程式：(2111S S 2211S 22Q Q t Q Q t Q Q t Q Q t *+ - + - + - + 下 1 下 2 上 1 上 2 下 1下2上1上2二二联立槽蓄曲线求解：1111222111222t kx t kx k kx tQ Q Q Q k kx t k kx t k kx t* - +- + + - + - + - + 下 2 下上 2 上 1二0101 1.0t kx C

16、 k kx t t kx C C C C k kx t k kx t C k kx t A- (I - + Al A+ I+ + - + A I- A II-+.U-令-注=0121Q C Q C Q C Q *. + +. 下2下上2上1则-即马斯京根法河道洪水演算程。第三节泄量模型泄量与库水位、开度关系，在合理确定相关参数后，其推求过程的计算工作量异 常庞大。泄量计算首先是假定已知闸前水头，然后计算闸门泄量，并利于推求水面线 法，逐断面试算得出对应库水位值。对应于某个库水位的泄量值，更是需要经大量的试算过程求得。1、闸门出流已知闸前水头，计算闸门泄量，关键是确定对应开度下的闸门水流形态（流

17、或堰 流。水力学对该临界点确定的比较笼统弧形闸门00.75e H -、平板闸门0.65e H -确定为临界点。准确的说，临界点是对应于某特定库水位，闸门堰流泄量等于流泄量时0e H值。 经过对岸堤水库弧形闸门的计算，不同库水位下的临界点不是一个固定数据,而是一 个围值:00.710.76eH -洪预报动态解析平台，计算闸门泄量时，首先就是确定对应该特定库水位下的临 界点。堰流泄量：32Q mB H 堰其中:.2(1k Hn nb生=-+ -（侧向收缩系数B二溢洪道净宽（m H二堰顶水头（m m二垂直收缩综合系数（驼峰堰二0.42 n二开启数b二单净宽 k 6二边墩形状系数0g二中墩形状系数流泄

18、量:0s Q o四二岸堤水库溢洪闸流量系数H四二-（适应围H软件系按照堰流等于流的临界条件（绝对条件,通过试算推求对应闸门开度和闸前水头,最终得到对应某库水位真实的临界 点0e H值。计算某闸前水头的泄量时，首先依据该临界点，确定相应泄量公式(流、堰流计算 特定闸前下的泄量值，再利用推求水面线的法，计算对应于该特定闸前水头下的库水 位值。这也是个一系列的试算过程。2、推求水面线即:利于闸前水头和求得的泄量值，通过推求闸前明渠水面线，计算对应库水位。溢洪道闸前明渠水面线需逐个断面依次推求。计算公式采用水文调查规 (SL196-97附录C天然河道水面曲线计算公式：IU+ +=g V g V L K

19、 Q K Q Z Z L u L u L u 221(21222222 a式中:L u Z Z、上、下断面水位,m ;Q 流量,m 3/s ;L u V V、上、下断面平均流速,m/s ;L u K K、上、下断面流量模数,m 3/s ;3/21AR nK =; n 糙率；A 过水面积,m 2;R水力半径,m ;a 断面扩散系数，若L u V V ,5.0=a; L上、下断面间 距,m ;g 重力加速度,m/s 2;g =9.81m/s 2。3、闸门泄量泄量计算采用递归算法。设定闸前水头初始值等于库水位，计算闸门出流和对 应的明渠水面线,得到一个新库水位值，利用与其初始值差反复校正库水位值,当

20、误差 小于0.005m时，递归算计算终止,得到该库水位下的泄量值。第四节调洪演算调洪演算就是个水量平衡计算过程。即某一时段，进入水库的水量与水库下泄 水量之差,应等于该时段水库蓄水量的变化量。用公式表示：12212122V V t q q t Q Q -二+-+式中:1Q、2Q 时段初、末入库流量；1q、2q 时段初、末出库流量；1V、2V 时段初、末水库蓄水量；t A一时段长，取t A=1.0h ;洪预报解析平台的演算过程就是根据入库洪水过程及水库的泄量关系（Vq，按照水量平衡关系，采f用递归试算法完成。第五节控运案水库控制运用案的制定，是在调洪演算的基础上进行的。为了满足用户需求，洪 预报

21、解析平台，开发了调洪演算专区模块。软件V2.0版本的调洪演算模块，不仅可 导入自定义洪水，而且还支持演算过程数据文件的导出，可直观的校对演算数据，它也 是顺利通过控制运用案审批的重要依据。第四章扩展性设计1、流域参数支持外部修订（图15确保水库工程安全，是制定水库控制运用计划的基本原则,但这一原则与水库所 在地域的历史水文资料有较大差距。为确保工程安全，岸堤水库在制定控制运用案时，就是通过省水文图集查得 H24=120mm,Cv=0.8,Cs =3.5xCv;但分析岸堤水库所在地域的历史水文资料:H24= 105.4mm,Cv=0.45,Cs=2.5xCvo考虑到各水库均可能存在以上问题，软件

22、特使用者保留了流域参数外部修订功能。使用者可根据你的需要随时修改这些参 数。（图152、模型可拓展到任意水库的洪预报解析岸堤水库基本数据,均保存于软件运行目录下的,xx库dll水库.dll文件中（文件 保存目录可在软件配置文件中设置，这就保障了洪预报解析平台的可拓展性。对于 其他任意水库，只要编辑制作相应的dll基本数据文件，即可完成平台在水库间的拓 展。3、调洪规则外置（1编制调洪规则Jscript或VBScript脚本（附例置“重要目录”（xx库dll下。注:重要目录名在配置文件中定义。（2 配置“default.ini”文件（例目录临时目录二data重要目录二xx库dll短文件调洪规则二

23、VbThguize.txt库数据二水库.dll（3脚本（调洪规则参数：库水位、起调水位、闸门拟定泄量、固定引用量、溢流堰（坝 泄量、调洪第i时段、进库流量、防洪高水位、安全泄量、调 洪阶段进程标志（0、1、2。阶段进程标志说明:0二调洪初1二蓄洪期2二泄洪期（4调洪规则功能强大。允自行编制闸门泄量函数（闸门变开度、变开启数支持外部修订防汛指标等泄量校正固定引用二0标志二1说明：0.拟定泄量1.拟定泄量+脚本校正2.脚本泄量（例（5陆地月蒸发系数设置。4、洪水河道传播槽蓄参数自定义（v3.0版本小流域河道槽蓄参数依次为:洪水传播时间、槽蓄权重因子和洪水拓展成员数, 并保存为,.txt格式文件置于

24、,xx库dll目录下（保存目录可在软件配置文件中设置，相 关短文件名在流域参数修订窗口中定义。5、调洪控运案Pa值、P+PaR线号预定义（控运案Pa预定义净雨n=12时段长二1陆地蒸发文件二YueyinxiangPa.txt说明:Pa=前期影响（mmn=（P+paR线号Tc二单位时段长度（小时陆地月蒸发k值表二前期降雨影响系数表的短文件名，文件结构（可略：（系数1；系数126、洪水基流数据预定义河道k=1 （洪水折算调整系数，建议暂取1.07、固定引用泄量Q外置预定义(如:输水洞泄量+城市供水引用等泄量校正固定引用二0 (例8、版本v2.1还支持防汛指标外置自定义(优先防汛指标；限制水位二.0

25、；防洪高水位二177.8；安全泄量=2000vbScript:XieLiangXiudingfunction XieLiangXiuding ( kushuiwei , qitiaoshuiwei , nidingxieliang , gudingyinyong , yiliuyanxieliang , i, jinkuliuliang , fanghonggaoshuiwei , anquanxieliang , jieduanif(kushuiwei fanghonggaoshuiwei +0.005 ThenXieLiangXiuding = nidingxieliangExit fun

26、ctionend ifif(jieduan=1 Thenif(nidingxieliang jinkuliuliang ThenXieLiangXiuding = jinkuliuliangelseXieLiangXiuding = nidingxieliangend ifif(XieLiangXiuding anquanxieliang ThenXieLiangXiuding = anquanxieliang end ifExit functionend ifif(jieduan=2 Thenif(nidingxieliang anquanxieliang ThenXieLiangXiudi

27、ng = anquanxieliang elseXieLiangXiuding = nidingxieliangend ifend ifJScript:XieLiangXiudingfunction XieLiangXiuding( kushuiwei , qitiaoshuiwei , nidingxieliang , gudingyinyong , yiliuyanxieliang , i, jinkuliuliang , fanghonggaoshuiwei , anquanxieliang , jieduan if(kushuiwei jinkuliuliang & jieduan =

28、 0 XieLiangXiuding = jinkuliuliang;if(nidingxieliang anquanxieliang & jieduan = 1 XieLiangXiuding = anquanxieliangreturn XieLiangXiuding;第五章调洪实例指时17&. 8&52382. 342QOO. 0014时17&. 7551665. 172000.00均时17&.50945. 882000,0012.14& 5S0. 37调洪累计泄水量（万80366.00调洪累计来水量（万97996.00第六章课目攻关概况2011.102012.02研发计算机存数据库。2

29、012.092014.08研究GDI +函数，攻克洪预报解析平台的界面绘制问题。2012.102012.12 研究通过 WinHttpRequest 对象访问 Internet。2013.072013.12攻克非线性相关函数用数学表达式，对已知多点用平滑曲线连 接问题；2014.012014.06用汇编代码实现非线性函数数学表达式的计算问题。2014.052014.07研究Excel的伽玛函数、概率密度函数和累积函数，解决洪水 计算问题。2014.102015.05构建岸堤水库，洪预报动态解析平台。2015.12月份软件升级到v2.0v2.25。2016.1月份软件升级到v3.00，洪水传播采用马斯京根河道洪水演算程式计算。第七章使用说明书第一节洪水预报1、修订调洪（泄洪规则：软件外置调洪规 则保存于.xx库dllVbThguize.txt目录（通过修改配置文件default.ini的相关信 息，可修改该目录；规则支持VBScript 和Jscript 语言，编写语法必须谨，并在英 文状态下输入。软件自带脚本意义如下：a.库水位低于起调水位时，不泄洪。b. 调洪初期，来水量小于泄水量时，按来水量泄洪。c.库水位小于防洪高水位时，按 下游安全泄量泄洪。d.库水位高于防洪高水位时，按照闸门拟定操作案泄洪。2、 修订流域参数：鼠标在对应流域，右键单击弹出流域参数修订窗口，在