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

岸堤水库雨洪资源解析书二O一五年六月一日作者:李文华电话:135“邮箱:fblwh150 目录第一章概述（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第一节洪水预报（31第二节调洪演算 （33第三节其他计算（33附件课题研发小组成员名单错误!未定义书签。第一章概述控制和预见洪水，让洪水变为一种资源,实现科学预见、动态管理、合理利用，是 本课题的研究对象。科学控制洪水，真正能够对洪水运用自如，其首要问题是准确解析、及时预报，掌 握洪水动态。但目前实际应用中，对水库防洪兴利控制运用，还仅限于依靠库水位的 变化，结合下游河道的承受能力，试探性的调节洪水，这种洪水调整模式，具有较大的盲 目性,理论方面的支撑相对不足。当前,各水库防汛主体单位，均制定了相应的水库控制运用方案。如岸堤水 库防洪调度图（图1，但这些方案的编制和批复仅表现为粗线条和原则性的界定，是在 进行大量假定的基础上进行编制的，应用中的可操作性相对欠缺，在实践中仅具有指 导意义。（图1洪水调度控制方案的编制,偏离实际应用，存在的突出问题,主要表现在以下几个 方面：1、假定了降雨的空间分配是均匀的，即整个流域降雨分布是均等的。但实际降 雨,特别是流域面积稍大的水库，降雨的空间分布几乎不可能是均等。2、事先拟定了 24小时降雨在1日内各时段上的雨量分配。但实际降雨在时段 上的分配，是个随机的不确定因素。3、控制运用方案的编制，起调水位为汛中限制水位，但实际降雨前的库水位，却 几乎不可能恰巧是汛中限制水位。4、所有闸门同开度启用，与实际控制运用也不相符。5、控运方案的编制基础是净雨量，但实际操作中我们所能掌握的是降雨量。把 降雨量合理的转换为净雨量，也成了水库实际控制运用要解决的问题。岸堤水库雨洪资源动态解析预报平台，较为圆满的解决了以上五方面的问题。第二章功能简介第一节功能特点1、采用.dl 1数据调用模式，获取水库基本数据。软件具有非常理想的可扩展性, 支持任何某水库的洪水资源解析。2、流域洪水模型参数,支持外部修订（需慎重修订;上游中小型水库、塘坝等的 拦蓄作用，可通过修订流域权重参数，简化上游拦蓄效果的修正计算。3、数据支持电子表格输出和图形输出。可导出对应某次降雨的洪水过程线文 件、对应某泄洪方案的库水位动态变化过程线文件、调洪演算成果文件及其他数据 的导出文件等。4、支持具有指导意义的控制运用方案的编制计算。v2.0及以后版本,还支持 演算过程数据文件（*.xls的导出。5、支持库水位库容水面积对应关系转换（附:溢流堰泄量;支持库水位闸 门开度溢洪道泄量对应关系转换;支持降雨信息、降雨重现期等其他数据转换，支 持逆推进库流量等计算功能。6、界面采用GDI+函数绘制;强化鼠标移动算法,瞬间获得流域、计算钮等虚拟 按键操作指令;频繁调用的关键代码，特别是v2.20及以后版本,大量的嵌入了汇编语言代码，直接操作 计算机CPU和内存等硬件设备,数据处理快速、高效。7、全面支持包括调洪规则等大量的信息自定义。第二节软件画面1、启动画面系透明.png文件（图2,其背景系洪水创意手绘图。（图22、主画面背景系岸堤水库流域布局图（图3（图3第三节运算功能1、前期降雨情况数据输入功能（图4前期降雨影响输入模式下，还可通过右键菜单，调用软件内置的，前期降雨影响量 *计算器，进行前期降雨影响计算;支持快速自动录入Pa值。（图42、流域降雨数据，按真实时程分配输入（图5（图53、洪水过程线屏幕输出（图6,支持数据文件导出（*.xls和*.png。（图64、调洪库水位过程线屏幕输出（图7,支持数据文件导出;支持任意起调水位和 任意闸门操作形式下的洪水控制运用计算。（图74、控制运用方案编制计算模块（图8。支持任意起调水位（默认为限制水位设置; 支持外置调洪规则设置。软件作者配置的调洪规则意义为:起调水位以下不泄洪，防 洪高水位以下，按下游河道安全泄量调洪（保下游,库水位超过防洪高水位后,敞闸泄 洪（保工程;并依据综合利用水库调度通则（水利部水管199361号对防洪调度图 的编制要求,对降雨做如下假定：假定1:流域均匀净雨量;假定2:流域净雨的时段分配，恒定匹配于水库所在地域；假定3:所有闸门同高度启用。软件V2.0及以后版本的调洪演算模块，可导入自定义洪水;可导出调洪数据文件, 演算过程清晰直观。.（图85、库水位库容水面积，对应关系表（图9。（图9 6、闸门泄量计算功能区（图10（图107、进库流量逆推功能（图11（图11 8、降雨量重现期（图12（图129、软件的v2.0及以后版本，还设计了，脚本调试器'。调试器直接运行脚本文件 VbThguize.txt,并输出计算结果，软件使用者可验证自己编写的脚本文件是否正确。10、软件v2.0及以后版本,前期降雨影响的计算，是以陆地蒸发系数的模式予以 外部设置，文件名是在配置文件中予以指定YueyinxiangPa.txt的形式外置，允许用户 根据当地的气候条件，参照已有资料随时予以调整和更新。11、外置调洪规则脚本（VbThguize.txt，功能强大，对于调洪演算和洪水预报，可 轻松的实现闸门的动态操作（详见:第四章扩展性设计。12、v2.20及以后版本，支持洪水时段长自定义。13、v2.24版本还支持小流域洪水数据导出。14、v3.00版本的洪水河道传播，采用了被国内外广泛应用的马斯京根算法。第四节气象云图及气象雷达1、，风云二号'气象卫星实况数据调用（图132、，当地、当前'气象雷达数据调用（图14（图13（图14第三章数学模型第一节洪水模型1、瞬时单位线为确保软件解析功能在各水库之间的通用性和可拓展性，洪水单位线采用传统 的瞬时单位线。瞬时单位线的理论依据就是把洪水过程看作是流域n个相同的线性水库串联调 节后的汇流结果,是个已经过理论论证和推导，并被大家所接受的成熟结论。其数学表达式：11u(t=(n t kt ek n k -/)IA其中:n =线性水库数量,反应流域调节能力；(n r= n的伽玛函数；K =线性水库的调节系数。e =常数(自然对数的底数瞬时单位线S曲线数学表达式：101S(t=n t ktkt ed t k n k -(. iA打常规利于瞬时单位线推求洪水过程线时,我们就是通过S曲线，计算无因次时段 单位线，最后转换为有因次时段单位线。这种汇流计算，不仅是工作量大，而且通过 S(t/k ,n推求参数的时段单位线时，一方面因使用中间参数(t/k增加了舍入误差，另一 方面(t/k和n的双向内插，增加了方法误差。计算结果的精度受到了很大影响。微软Excel办公组件，为我们提供了 CAMMADIST函数，是洪水计算较为便捷理 想的计算方案。2、CAMMADIST函数语法CAMMADIST(x，alpha，beta，cumulative式中:x是自变量，用来计算伽玛分布函数 的值；Alpha和beta均为相关参数。Cumulative是个逻辑参数：Cumulative =false时,返回概率密度函数：xf x x e apaappa-=rCumulative =true时,返回累积分布函数:11F(,(xxx x ed x apaappa-=rj0.330.270.20.171C2C 1212M M FJ RT 0.196F km J m m R mm T h M 3M M 3aa-，一M=>n山东省瞬时单位线参数经验公式：=山丘区=流域面积(河道干流平均坡度(/净雨量(净雨历时(若:则 3 若:则 312M n k 1M n kM n = =单位线参数与、的关系：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 F t t u t t tF t h t i n k t i t n k t -= - A= =.-+.-. ¥无因次时段单位线:有因次时段单位线:流域汇流计算公式：4、流域洪水错时叠加CAMMADIST函数本身并不能解决降雨在空间上的不均匀分布问题。透过微 分观点，我们可以把全流域进行细分，分割成若十个微小流域，我们完全可以认为微小 流域内的降雨是均布的。鉴于目前技术条件和设备配置，我们还无法得到无限微小 流域的降雨量，因此，流域的细分，不能真正的进行无限分割，必须结合雨量观测设备在 全流域内的真实布局情况予以分割。利于微软的CAMMADIST函数计算洪水过程，流域划分越细,计算结果越合 理。按照水库全流域雨量检测设备的分布,进行了二次分割，把整个水库流域,细分成 了多个小流域，依次计算洪水并经过河道传播和叠加计算，最终得到流域坝前洪水过 程线。岸堤水库上游在水库整个流域内,分布着4个水文站和5个雨量观测点，按照这 些雨量观测设备的布局,把整个流域又分割成了 9个小流域。这些小流域的洪水，通 过调用微软的CAMMADIST函数得到相应过程线，各个小流域洪水过程线经主河道 传播后，错时段叠加得到洪水演算所需要的坝前洪水动态数据。第二节洪水传播流域洪水在河道中的传播,采用马斯京根算法。马斯京根法假定,对某特定河道 段,其河道槽蓄量S与流量Q之间，存在着一种固有关系,即槽蓄曲线。槽蓄曲线反映河段的水力学特性。(''S Q K 1Q 1K x K Q Q x Q x Q x Q x f Q - - +- - +-下上下上下上=,=(式 中:=(示储流量=流量比重因子河段汇流时间(物理意义建立如下的河道水量平衡方程式：(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*. -. +-_. + + -+ -+ -+ AT 2 下上 2 上 1=0101220.50.50.5 1.00.50.50.5t kx C k kx t t kx C C C C k kx t k kx t C k kx t A-fI -+ A I A+ I+ I - , A I- A II -+ A I=令:=注：=0121Q C Q C Q C Q * + +下 2 下上 2 上 1 则:=即:马斯京根法河道洪水演算方程。第三节泄量模型泄量与库水位、开度关系，在合理确定相关参数后,其推求过程的计算工作量异 常庞大。泄量计算首先是假定已知闸前水头，然后计算闸门泄量，并利于推求水面线 方法，逐断面试算得出对应库水位值。对应于某个库水位的泄量值，更是需要经大量的试算过程求得。1、闸门出流已知闸前水头，计算闸门泄量，关键是确定对应开度下的闸门水流形态（孔流或堰 流。水力学对该临界点确定的比较笼统:弧形闸门00.75e H =、平板闸门0.65e H =确定为临界点。准确的说，临界点是对应于某特定库水位，闸门堰流泄量等于孔流泄量时0e H 值。经过对岸堤水库弧形闸门的计算，不同库水位下的临界点不是一个固定数据，而 是一个范围值:00.710.76eH -洪预报动态解析平台，计算闸门泄量时，首先就是确定对应该特定库水位下的临 界点。堰流泄量:32Q mB H £ 堰%其中:010.2（1k Hn nb£生=-+-（侧向收缩系数B =溢洪道净宽（m H 二堰顶水头（mm =垂直收缩综合系数（驼峰堰=0.42 n =开启孔数b 二单孔净宽k孑边墩形状系数0孑中墩形状系数孔流泄量:0s Q o四孔=%岸堤水库溢洪闸流量系数00.6850.19eH四=-（适应范围0.10.75eH软件系按照堰流等于孔流的临界条件（绝对条件，通过试算推求对应闸门开度和闸前水头，最终得到对应某库水位真实的临界 点0e H值。计算某闸前水头的泄量时，首先依据该临界点，确定相应泄量公式（孔流、堰流计 算特定闸前下的泄量值，再利用推求水面线的方法，计算对应于该特定闸前水头下的 库水位值。这也是个一系列的试算过程。2、推求水面线即:利于闸前水头和求得的泄量值，通过推求闸前明渠水面线，计算对应库水位。溢洪道闸前明渠水面线需逐个断面依次推求。计算公式采用水文调查规范（SL196-97附录C天然河道水面曲线计算公式：IU)"/-IU1 "/+=g V g V L K Q K Q Z Z L u L u L u 221(21222222a式中: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 <,0=a;若 L u V V >,5.0=a; L 上、下断面间 距,m ;g 重力加速度,m/s 2;g =9.81m/s 2。3、闸门泄量泄量计算采用递归算法。设定闸前水头初始值等于库水位，计算闸门出流和对 应的明渠水面线,得到一个新库水位值，利用与其初始值差反复校正库水位值,当误差 小于0.005m时，递归算计算终止,得到该库水位下的泄量值。第四节调洪演算调洪演算就是个水量平衡计算过程。即某一时段，进入水库的水量与水库下泄 水量之差，应等于该时段内水库蓄水量的变化量。用公式表示：12212122V V t q q t Q Q -=+-+式中:1Q、2Q 时段初、末入库流量；1q、2q时段初、末出库流量；1V、2V时段初、末水库蓄水量；t A时段长，取t A=1.0h ;洪预报解析平台的演算过程就是根据入库洪水过程及水库的泄量关系（Vq，按照水量平衡关系，采用递归试算法完成。第五节控运方案水库控制运用方案的制定，是在调洪演算的基础上进行的。为了满足用户需求, 洪预报解析平台,开发了调洪演算专区模块。软件V2.0版本的调洪演算模块，不仅 许可导入自定义洪水，而且还支持演算过程数据文件的导出，可直观的校对演算数据, 它也是顺利通过控制运用方案审批的重要依据。第四章扩展性设计1、流域参数支持外部修订（图15确保水库工程安全，是制定水库控制运用计划的基本原则,但这一原则与水库所 在地域的历史水文资料有较大差距。为确保工程安全，岸堤水库在制定控制运用方案时，就是通过山东省水文图 集查得 H24=120mm,Cv=0.8,Cs =3.5xCv;但分析岸堤水库所在地域的历史水文资料:H24= 105.4mm,Cv=0.45,Cs=2.5xCv。考虑到各水库均可能存在以上问题，软件特会合使用者保留了流域参数外部修订功能。使用者可根据你的需要随时修改这些参 数。（图152、模型可拓展到任意水库的洪预报解析岸堤水库基本数据，均保存于软件运行目录下的,xx库dii水库.dir文件中（文件 保存目录可在软件配置文件中设置，这就保障了洪预报解析平台的可拓展性。对于其他任意水库，只要编辑制作相应的dll基本数据文件，即可完成平台在水库间的拓 展。3、调洪规则外置（1编制调洪规则Jscript或VBScript脚本（附例置“重要目录”（xx库dll下。注:重要目录名在配置文件中定义。（2 配置 “default.ini” 文件（例目录临时目录=data重要目录=xx库dll短文件调洪规则=VbThguize.txt库数据=水库.dll（3脚本（调洪规则参数：库水位、起调水位、闸门拟定泄量、固定引用量、溢流堰（坝泄量、调洪第i时段、进库流量、防洪高水位、安全泄量、调洪阶段进程标志（0、1、2。阶段进程标志说明:0=调洪初1=蓄洪期2=泄洪期（4调洪规则功能强大。允许自行编制闸门泄量函数（闸门变开度、变开启孔数支持外部修订防汛指标等泄量校正固定引用=0标志=1说明：0.拟定泄量1.拟定泄量+脚本校正2.脚本泄量（例（5陆地月蒸发系数设置。4、洪水河道传播槽蓄参数自定义（v3.0版本小流域河道槽蓄参数依次为:洪水传播时间、槽蓄权重因子和洪水拓展成员数, 并保存为,.txf格式文件，置于夕库dll'目录下（保存目录可在软件配置文件中设置,相 关短文件名在流域参数修订窗口中定义。5、调洪控运方案Pa值、P+PaR线号预定义（控运方案Pa预定义净雨n=12时段长=1陆地蒸发文件=YueyinxiangPa.txt说明:Pa=前期影响（mmn=（P+paR线号Tc=单位时段长度（小时陆地月蒸发k值表=前期降雨影响系数表的短文件名，文件结构（可略：（系数1系数126、洪水基流数据预定义河道k=1 （洪水折算调整系数，建议暂取1.0基流=27、固定引用泄量Q外置预定义（如:输水洞泄量+城市供水引用等泄量校正固定引用=0（例8、版本v2.1还支持防汛指标外置自定义（优先防汛指标；限制水位=173.0;防洪高水位=177.8；安全泄量=2000vbScript:XieLiangXiudingfunction XieLiangXiuding （ kushuiwei , qitiaoshuiwei , nidingxieliang , gudingyinyong , yiliuyanxieliang , i, jinkuliuliang , fanghonggaoshuiwei , anquanxieliang , jieduanif（kushuiwei < qitiaoshuiwei -0.005 ThenXieLiangXiuding = 0Exit functionend ifif(kushuiwei > fanghonggaoshuiwei +0.005 ThenXieLiangXiuding = nidingxieliangExit functionend ifif(jieduan=1 Thenif(nidingxieliang > jinkuliuliang ThenXieLiangXiuding = jinkuliuliangelseXieLiangXiuding = nidingxieliangend ifif(XieLiangXiuding > anquanxieliang ThenXieLiangXiuding = anquanxieliangend ifExit functionend ifif(jieduan=2 Thenif(nidingxieliang > anquanxieliang ThenXieLiangXiuding = anquanxieliangelseXieLiangXiuding = nidingxieliangend ifend ifJScript:XieLiangXiudingfunction XieLiangXiuding( kushuiwei , qitiaoshuiwei , nidingxieliang , gudingyinyong , yiliuyanxieliang , i, jinkuliuliang , fanghonggaoshuiwei , anquanxieliang , jieduan if(kushuiwei <= qitiaoshuiwei + 0.005 XieLiangXiuding = 0;else if(nidingxieliang > jinkuliuliang && jieduan = 0 XieLiangXiuding = jinkuliuliang;if(nidingxieliang > anquanxieliang && jieduan = 1 XieLiangXiuding = anquanxieliangreturn XieLiangXiuding;第五章调洪实例If*商加II 1 lilL6 3CaCJ16 Bt. E"m帝.姓175. W5WEvT77M.5122E. MifairIWJ3 用厂0_邪m脸1涪wi沱ao 】梏.画&1.L?1FS.M rFsijinn5史M!Um wJ17 17. E53 h EM E'«K243.3WL.CU口用irr.ia临殖55的队54AJ7k15S2KCL 匹MUL.'H.a si瑚皿脸tS5<iT.±73物：.笫缶用i&i：ran,如汕虻布M33li5a!M!IKl)E-4 II 1 LIDHmssrriG_i3L79tl.1H1KL9J-时i-1.- =：藉3-179£±.!»弓跖政Wraw舟n醐出腭ITS. E(ECU.飞VQM.JIJl.LE,时E«U.t-1汹二:H5H5151. «JKlS.TiE-IRK . L3昭电匹LK«9.D72S7ZK商St 7-L1皿齐«：I1.5IMKL. I.-513时176.8652382.342000. DO12.1414时 '17&.7551665.172000.003. 5815时*176.6350945.882000, DO0.37技调洪累计泄水量（万方80366.00调洪累计来水量（万方97996.00第六章课目攻关概况2011.102012.02研发计算机内存数据库。2012.092014.08研究GDI+函数，攻克洪预报解析平台的界面绘制问题。2012.102012.12 研究通过 WinHttpRequest 对象访问 Interneto2013.072013.12攻克非线性相关函数用数学表达式，对已知多点用平滑曲线连 接问题；2014.012014.06用汇编代码实现非线性函数数学表达式的计算问题。2014.052014.07研究Excel的伽玛函数、概率密度函数和累积函数，解决洪水 计算问题。2014.102015.05构建岸堤水库，洪预报'动态解析平台。2015.12月份软件升级到v2.0v2.25o2016.1月份软件升级到v3.00，洪水传播采用马斯京根河道洪水演算方程式计 算。第七章使用说明书第一节洪水预报1、修订调洪（泄洪规则：软件外置调洪规 则保存于.xx库dll VbThguize.txt目录（通过修改配置文件default.ini的相关信 息，可修改该目录；规则支持VBScript和Jscript语言，编写语法必须严谨，并在 英文状态下输入。软件自带脚本意义如下：a.库水位低于起调水位时，不泄洪。b. 调洪初期，来水量小于泄水量时，按来水量泄洪。c.库水位小于防洪高水位时，按 下游安全泄量泄洪。d.库水位高于防洪高水位时，按照闸门拟定操作方案泄洪。2、修订流域参数：鼠标在对应流域，右键单击弹出流域参数修订窗口，在对应栏 目中修改即可。温馨提示：建议慎重修改，最好是由水文专家指导修改。3、输入 各流域的降雨信息：鼠标在对应流域，左键单击弹出降雨信息输入窗口，按栏目 意义输入即可。最终可通过点击降雨、面雨、净雨按钮，检查校对降雨输入信 息的正确性。4、输入前期影响量。可输入经验数据，也可在降雨前期影响输入模 式下，右键选择降雨前期影响计算器，输入相应的前期降雨情况（最好能输入 前期34个月的所有降雨信息）通过计算确定前期影响量。支持快速赋值。温馨提示：前期降雨影响计算器，系根据外置的陆地月蒸发k值表 =YueyinxiangPa.txt，依次叠加计算得来。因此计算前期前，要依据当地的气候条 件，事先合理设置各月陆地蒸发系数（山东省境内可直接采用软件自带数据文 件）。a.土壤饱和含水量自动调用流域参数中的对应数据，但支持即时修订。b. 计算终止日系调用的当前时间，需注意是否和计算意图想匹配。5、查看或导出洪 水过程线（可略）：点击根据流域实际降雨推求坝前洪水过程线查看洪水过 程。温馨提示：右键单击，弹出对应的保存数据菜单。6、设定闸门操作方案：可 选定固定闸门泄量或固定闸门开度方案。a.闸门孔号组系闸门开启孔数组合， 孔号间用分号分割（支持其他字符的分割）b.固定闸门泄量方案时，如果闸门开 启孔数不能满足预定的固定泄量，则按实际闸门泄量数据演算。c.固定闸门开店 方案时，固定开度系组合数据。数据需要和闸门开启孔数组相对应，数据集合用 分号分割（支持其他字符的分割）。允许简化输入（略。7、点击推求最高雍水 位按钮，查看对应泄洪方案下的库水位动态变动过程曲线及相应的调洪演算成果数据。温馨提示：右键单 击，弹出对应的保存数据菜单。第二节调洪演算1、修订调洪（泄洪规则：同洪水 预报。2、修订流域参数：同洪水预报。温馨提示：调洪演算时，一般是按照水库 全流域进行洪水计算，但软件把水库上游流域，按照上游测雨设备的布局，进行 了二次细分。故此可以用某一流域模拟替代全流域，即：定义该流域的加权系 数是1.0，其他流域相应的加权系数为0.0即可。3、点击调洪演算按钮进行 调洪演算（分别支持降雨和净雨输入）输出调洪演算综合成果。4、v2.0及其以后 版本，支持外部导入已有洪水过程线；调洪演算模式下，右键可保存调洪演算过 程数据。第三节其他计算包括库水位、库容、水面积转换；闸门泄量计算；逆推 进库流量；降雨量、频率转换等。（略）