简介《工程水文分析计算集成应用软件PHAC v2.28使用说明书》由贵州省水利水电勘测设计研究院编写面向水利水电工程技术人员、水文分析计算人员及相关专业学生系统讲解该集成软件在水文频率分析、暴雨洪水计算、水库库容与调洪演算等场景中的使用方法。压缩包内仅含1个PDF文档大小约4.43MB全文覆盖软件安装、系统启动、Windows窗口操作、数据库通用操作及15个计算单元的具体操作流程并对P—Ⅲ型频率曲线分析、河道加权平均比降、水库泥沙淤积、河道水面线推算等模块作了详细说明。目前已有791人学习适合需要借助PHAC完成工程水文分析计算、或希望掌握该软件操作细节的用户。借助这份说明书读者可快速了解软件模块化结构、软硬件环境要求并参照计算过程简述与图形绘制环境说明完成常用水文计算任务。1. 工程水文分析计算集成应用软件 v2.28 到底解决什么问题工程水文分析计算集成应用软件 v2.28 对应的是一类把设计洪水全流程搬进同一工程文件的工具从暴雨统计参数、产汇流计算到调洪演算和成果整编。真正用过的人会告诉你价值不在算法多先进而在中间过程可追溯——老办法用 Excel 排 P-III 频率曲线算完调洪很难说清楚哪个参数改过之后影响了最终水位。v2.28 把工程文件、成果表和过程日志串成一条线适合水利设计院、防汛部门、咨询单位里做防洪评价、水库安全鉴定和排涝规划的人。新版大多在成果输出和数据接口上做了调整升级上来的老用户更容易踩环境兼容的坑。2. 装 v2.28 之前先核对环境部署完做一次命令行验证2.1 装 v2.28 之前先核对三样东西数据库位数、授权组件、运行库这类软件在工程单位落地时最容易出的问题不是安装包损坏而是环境不干净。我一般先确认三件事数据库实例是 32 位还是 64 位授权组件加密锁驱动或授权服务是否安装以及系统里有没有老版本残留的运行库。v2.28 如果带空间分析模块还牵扯到 GIS 组件的注册权限普通用户登录时往往注册失败。检查操作系统位数的命令如下systeminfo | findstr /C:OS Name /C:System Type/C:是 findstr 指定搜索字符串的参数System Type一行会显示 x64-based PC 或 x86-based PC。数据库方面如果是 SQL Server用 sqlcmd 验证实例连通性sqlcmd -S .\SQLEXPRESS -E -Q SELECT VERSION-S指定实例名-E表示用 Windows 身份认证-Q执行完查询后立即退出。连不上时先看服务是否启动net start | findstr /i sql。这里要说清楚命令只做连通性验证v2.28 建库脚本一般在第一次启动时执行不要在装软件前手工建库库名和字符集不对反而起不来。2.2 用命令验证安装结果而不是只双击图标装上之后很多人习惯双击图标转圈几秒后报“许可错误”或“无法连接数据库”才回头排查。我习惯装完先不启动界面按下面顺序过一遍sc query type service state all | findstr /i sentinel hasp safenetsc query不加斜杠是 Windows 服务控制命令的标准写法/i是 findstr 的忽略大小写参数type service state all里的等号后面必须带空格。厂商不同授权服务名也不一样常见关键字包括 Sentinel、HASP 和 Safenet如果一条都查不到再打开“服务”管理界面找找厂商缩写开头的服务名。除此之外建议执行一次程序完整性校验cd /d C:\Program Files\HydroCalc\v2.28 dir *.dll *.exe/d切换盘符和目录这里只是举例实际安装路径通常在安装时可选。你看到 exe 和 dll 同目录不代表组件都注册成功真正的验证是打开软件进任意一个工程能正常弹出频率适线图窗口才算过。我遇到过多次目录文件齐全、界面可以打开、但一点“计算”就闪退的情况最后定位到是 VC 运行库缺失装好对应版本后恢复正常。2.3 从目录结构和日志判断软件是否进入可用状态v2.28 这类集成软件安装目录下通常会有 bin、config、log、template 几个子目录不同模块的数据文件不会堆在根目录。下面是一份常见目录安排不同厂家略有差异但思路一致目录一般存放内容出错时看什么bin主程序、算法动态库exe 和关键 dll 是否在同一层级是否存在位数混装config数据库连接串、模块开关连接串里的实例名、账号密码、字符集是否与现场一致log运行日志、算法中间日志计算报错时找当天日志看异常堆栈里的类名templateWord/Excel 成果模板模板里的书签或域代码是否被改动过project用户工程文件备份时整个目录打包不要只拷 .hyd 单文件安装完成后打开 log 目录看最新的日志文件能搜到类似init engine ok、license verified的字样基本可以判断授权和算法引擎都已加载。日志目录没有写入权限时软件会静默跳过日志输出这也是个坑——排查问题时会发现什么都没留下所以最好把安装目录的写入权限放给实际运行账号。3. 从雨量到调洪成果v2.28 的完整计算链路3.1 录入流域特征与暴雨参数时的口径问题v2.28 的工程文件一般分成基础资料、计算公式、成果三块。基础资料里的流域面积 F、主河长 L、河道比降 J这三个参数直接进推理公式和单位线计算取值口径必须和所依据的水文手册一致。常见误用是设计暴雨参数 A、C、b、n 直接从手册抄却不注意手册适用的历时区间比如某地区手册公式只对 10 分钟到 24 小时区间有效往里录入 3 天暴雨就明显外延算出的设计流量会异常偏大。录入界面里通常还要求填汇流参数 m这个值在小流域推理公式里很敏感m 从 0.7 调到 1.0洪峰流量可能差 30% 以上。v2.28 里这个参数一般放在“产汇流参数”页签旁边有手册参考值但每省分区图不同参考值只对应典型分区不能直接替代综合成果。我一般把参数来源写进“备注”字段导出成果时带着来源后续审查会用到。设计暴雨时程分配也在这里选是采用当地暴雨图集推荐的典型过程还是按均匀概化。这两者算出来的洪峰流量数值差异很大选错后调洪结果也会错。3.2 设计洪水推算从 P-III 频率适线到产汇流计算参数录完后进入频率分析模块。v2.28 里常用做法是矩法估计初值再人工适线。界面允许选目标函数最小绝对离差和、最小二乘、相对离差最小等。同一个序列不同目标函数适出来的 Cv、Cs/Cv 有差异工程上一般按规范要求以目估适线为主、数值适线为辅。操作时注意一个细节频率曲线图上的经验点据可以手工拖动刷新重算这是这类软件的通用交互目的是让工程师判断“尾部点据”是否脱离群点。频率计算完成后生成设计雨力或设计流量接着进入产汇流计算。以推理公式法为例核心是联解洪峰流量 Qm 与汇流历时 τQm 0.278 * (S / τ^n) * F * μ 修正软件界面里这几个字母都有对应输入框S 是设计雨力n 是暴雨递减指数μ 是损失参数。需要特别留意的是公式中径流系数或损失参数的取值v2.28 里有“初损后损法”和“蓄满产流”两种产流模式选择前要确认流域下垫面类型。计算中间结果可以在“洪水过程线”页面看到逐时段的流量序列导出前先目测一遍洪峰数值有没有量级错误比最后统一检查高效得多。查中间结果表用 SQL 直接看SELECT period_no, rainfall_mm, runoff_mm, flow_m3s FROM calc_flow_process WHERE project_id P2024001 ORDER BY period_no;这段查询筛选出指定工程的逐时段雨量、产流量和流量值ORDER BY period_no保证按时段顺序输出用于核对峰值所在时段是否与降雨峰值时段匹配。实际中如果洪峰出现时间比降雨峰值滞后太多先检查汇流参数 m 是否偏大再看时段步长是否统一。3.3 调洪演算的关键边界条件与结果输出调洪演算环节v2.28 按库水位-库容、泄流能力曲线和入库洪水过程求解水量平衡方程。开始之前把三条曲线输对水位库容曲线、泄洪洞/溢洪道泄流曲线、起调水位。起调水位这个值经常被填错溢洪道无闸门时取堰顶高程有闸门时取汛限水位混淆之后調洪结果会明显偏离实际。计算参数里有时间步长和迭代精度。时间步长一般取 1 小时对涨水段很陡的小流域可以改到 0.5 小时代价是计算耗时增加但 v2.28 的调洪模块对这类规模的计算基本是秒级完成放开用。迭代精度控制的是试算收敛条件默认 0.01 米即可不要为了“精确”设到 0.0001 米一方面意义不大另一方面可能因数值震荡收敛不了。调洪算完主界面给出最高库水位、最大下泄流量和峰现时间。这三个数字和手算是可以对上的我一般先用简化三角形概化过程粗算一遍如果和 v2.28 的结果差距在 5% 以内说明边界输入基本没问题差距过大就先查起调水位再查泄流曲线的量纲。成果输出时选“调洪计算书”软件自动生成包含原始数据、过程和结论的文档这时检查一下单位显示是否为米、立方米每秒、小时导出 PDF 前把精度设成小数点后两位避免报告里出现一长串小数。4. 看懂 v2.28 的成果数据表导入导出才不吃亏4.1 成果数据表的分层结构与关联关系v2.28 一般用数据库表存工程数据和中间结果按“工程 → 方案 → 计算工况”三层组织。一个工程文件里可以建多个方案每个方案对应一组独立的参数集。理解这个分层结构比记住每个字段都重要审查报告时需要对比“原方案”和“加高方案”直接在软件里切方案比导出 Excel 手工对比可靠得多。常用数据表可以按功能分成以下几类表类型表名示例关键字段说明工程信息t_projectproject_id, project_name, basin_area工程级基本信息一个工程一行方案参数t_scheme_paramscheme_id, param_type, param_value不同方案的参数取值param_type 区分雨力、汇流参数等频率成果t_freq_resultscheme_id, return_period, design_value各重现期设计值单位需查字段单位表洪水过程t_flood_processperiod_no, flow, time_tag逐时段流量用于画过程线调洪成果t_routing_resultmax_level, max_discharge, peak_time调洪计算的特征值结果字段命名各版本有差异但结构上逃不出这套框架。拿到一个新工程文件我建议先打开表清单找到t_project开头的主表查project_id再顺着关联字段找方案表和成果表比在界面上点击深入要快。4.2 导入导出最容易踩的三个坑单位、编码、坐标v2.28 支持从 Excel 导入水库特征曲线也支持把成果导出成 Excel、CSV 和 Word。几个坑常年有人踩。第一Excel 模板的表头说明行不能删。导入模板通常前两行是字段说明和单位注释从第三行开始才是数据很多人在整理数据时把前两行删掉软件读不到表头就报“列数不匹配”。正确做法是保留模板原文件只在下方的数据区域做覆盖。第二CSV 编码问题。v2.28 在中文 Windows 下默认按 GBK 读写 CSV但很多自动生成的 CSV 是 UTF-8 格式直接导入会中文乱码或数字错位。用文本编辑器转换编码后再导入with open(input_utf8.csv, r, encodingutf-8-sig) as f: text f.read() with open(input_gbk.csv, w, encodinggbk) as f: f.write(text)utf-8-sig会去掉文件头部的 BOM 标记这是 UTF-8 带签名文件与无签名文件的主要差异。转换完打开文件目测一遍重点看中文列是否正常再导入软件。第三坐标和椭球参数。带空间分析能力的 v2.28 模块导入流域边界或 DEM 时投影坐标系不统一会导致面积计算差出好几倍。我处理过一份数据流域边界用的是 CGCS2000 3 度带DEM 是西安 80 坐标系软件识别不出差异直接按投影坐标算面积结果明显偏大。之后我坚持在导入前先统一成同一种坐标系并用软件的“属性查询”查看边界范围值。4.3 把计算成果送进报告模板的常见做法成果导出到 Word 报告v2.28 的模板机制一般是书签替换或域代码引用。用书签替换时先在 Word 模板里插入书签软件导出时把计算值写入书签位置。需要注意书签名称不要在模板里乱改v2.28 配置文件中定义了书签与数据字段的映射关系改了书签名就导不出数值。模板里有汇总表、过程表和附图三个部分。汇总是特征值过程是逐时段数据附图是频率曲线和调洪过程线。导出之前先检查模板中表格的行列数是否与软件约定一致比如过程表模板预设 100 行但你的计算时段超过 100导出后超出的时段会被丢弃。比较稳妥的方式是先在软件里把过程数据导成 Excel然后用 Word 的邮件合并功能把 Excel 数据带进报告模板绕开模板行数限制。v2.28 自带的模板导出适合快速出稿修改多的报告建议走外部数据链路。5. 高频报错与调参技巧拿脚本复核 v2.28 的成果5.1 频率曲线不收敛时先查这三个地方频率适线“算不出来”或“曲线跑到图外”这是 v2.28 使用者最常反馈的问题。我的排查顺序是固定的。先看输入序列里有没有 0 值或负值年最大降雨序列里出现 0 会导致矩法估计均值严重偏移处理方式是检查原始资料缺失年份要明确标注不要用 0 填充。再看 Cv 初值是否超出合理范围把矩法初值改成人工输入的 Cv0.5、Cs/Cv3.5多数能起算之后再细微调整。第三看样本容量序列短于 15 年时 P-III 适线的外延部分没有统计意义软件能算但结果不稳报告里要说明。5.2 用一段 Python 脚本对调洪成果做水量平衡复核调洪计算的本质是水量平衡所以最可靠的复核办法就是把入库、出库和库容变化重新对一遍。v2.28 导出的调洪过程表通常包含时段号、入库流量、出库流量和库水位下面这段脚本用相邻时段平均流量法检查水量闭合import csv inflow [] outflow [] levels [] with open(routing_result.csv, encodinggbk) as f: reader csv.DictReader(f) for row in reader: inflow.append(float(row[inflow])) outflow.append(float(row[outflow])) levels.append(float(row[level])) dt 3600 # 时间步长单位秒 total_in sum((inflow[i] inflow[i-1]) / 2 * dt for i in range(1, len(inflow))) total_out sum((outflow[i] outflow[i-1]) / 2 * dt for i in range(1, len(outflow))) diff total_in - total_out print(f入库总量 {total_in:.1f} m3, 出库总量 {total_out:.1f} m3, f差值 {diff:.1f} m3)起步逻辑是每个相邻时段取流量平均值乘以时间步长得该时段水量累加得到总量。encodinggbk对应 v2.28 导出的中文 CSV 默认编码文件是 UTF-8 时要改成 utf-8。运行后差值大于入库总量的 2%先查时间步长是不是 1 小时但解析成了别的单位再查库容曲线插值是否有明显台阶。闭合差正常才有把握说调洪成果可信。5.3 复核通过后再看峰现时间和单位换算水量平衡通过后还有三个点值得花半分钟确认。峰现时间与入库洪峰出现时间的先后关系调洪后出库峰现应该落后于入库峰现这是在所有水库调洪里都成立的物理约束如果软件算出来出库洪峰提前了基本是泄流曲线方向写反。最高库水位对应的库容增幅应该等于该时刻累计入库减去累计出库这个值和脚本里的diff是同源逻辑只是截取时段不同。单位换算上CSV 导出表头里库容是万立方米还是立方米流量是立方米每秒还是升每秒字段名不会写清楚靠数值量级判断一座中型水库的调洪库容通常千万立方米量级如果导出的数字只有四位先怀疑单位是万立方米。把这三项检查固化成习惯v2.28 的成果拿去评审才经得起追问。本文还有配套的精品资源点击获取