管道水力摩阻系数的敏感性分析

1 研究背景

管道输水具有效率高、环境影响小、拆迁占地少、运行调度安全可靠、维护管理简单等特点，广泛应用于城市给排水、农业灌溉、跨流（区）域调水等水利工程中。管道水力计算公式及其摩阻系数的合理选用是流体输送工程设计的关键，其选取合理与否直接影响水力计算成果的精度，进而影响工程的总体布局、设计规模、投资乃至运行费用。现有规范和工程设计常用的水力计算公式及相应的摩阻系数有［1-2］：达西-魏斯巴哈（Darcy-Weisbach）公式和达西摩阻系数λ，海森-威廉（Hazen-Williams）公式和海森-威廉系数Ch，谢才-曼宁（Chézy-Manning）公式和曼宁糙率系数n。

大量研究和工程实践表明水力摩阻系数不仅与管道内壁的粗糙度有关，而且受管径D、流速V、水流黏滞系数ν等的影响［3-4］。然而，国内现有规范给出的摩阻系数取值范围较宽，对摩阻系数影响因素的考虑不全面，导致工程设计，尤其是大口径输水管道的摩阻系数精确取值困难，部分工程的水力计算误差较大。例如，新疆北疆供水工程小洼槽倒虹吸，采用DN3100 玻璃钢夹砂管，曼宁糙率系数设计值n=0.0090，工程运行后实测n=0.0106，比设计值大18%，导致上下游进出口水位差比设计值增加0.6 m，不能满足设计输水流量要求［5］；三个泉倒虹吸采用钢管和PCCP 管的组合方案，其中DN2700 环氧内衬钢管的曼宁糙率系数设计值n=0.0120，实测n=0.0098，比设计值小23%，DN2800PCCP 管道的曼宁糙率系数设计值n=0.0135，实测n=0.0108，比设计值小25%，实际过流能力大于设计值，原设计方案偏于保守［6］；南水北调中线北京段DN4000PCCP 管道，曼宁糙率系数设计值n=0.0120，实测n=0.0101，比设计值小15%，实际运行水头损失（5.30 m）仅为设计值（7.64 m）的69%［7］。

土洞开挖采用全断面开挖的方式，“短进尺、勤支护”，利用挖掘机直接开挖，必要时在周边预留10～20 cm的土层人工开挖，洞口10 m范围由于埋深浅（3～10 m），每轮循环进尺控制在1.5 m以内；随后在施工中根据现场实践经验结合变形监测结果，加大循环进尺长度，到埋深较大段，最大单次循环进尺达到3 m，顶拱出现局部坍塌和掉块，说明土层不稳定，循环进尺达到最大长度。

C、D泊位主要装卸货种为柴油，其他危险品均为少量采用瓶装运输（如乙炔、油漆等），按装卸丙类危险品的二级码头进行消防设计。

因此，有必要通过建立不同水力摩阻系数间的换算关系，系统分析各摩阻系数对不同因素的敏感程度和变化规律，进而提出工程中水力计算公式选用和摩阻系数取值应注意的问题。

2 管道水力计算公式

2.1 达西-魏斯巴哈公式

式中：hf为管道沿程水头损失，m；λ为达西摩阻系数；D 为管道内径，m；L 为管段长度，m；V 为管道平均流速，m/s；g 为重力加速度，m/s2。达西-魏斯巴哈公式适用于不同流体（气体或液体）的不同流态（包括层流和紊流），达西摩阻系数λ通常采用柯尔勃洛克-怀特（Colebrook - White）公式计算：

青岛市作为全国第一批医养结合养老试点城市，其医养结合的发展进程一直在加快。根据研究结果以及对青岛市医养结合养老服务的调研，本文提出以下对策建议。

式中：k 为当量粗糙度，m；Re 为雷诺数。柯尔勃洛克-怀特公式是计算紊流区达西摩阻系数λ的半理论半经验公式，其适用范围涵盖了光滑管区、紊流过渡区和粗糙管区，大量的试验结果表明该式与实际商用管道的阻力试验结果吻合良好，被公认为是管道水力计算的标准方程［8］。

当水流处于紊流粗糙区时，雷诺数Re足够大，式（2）可简化为尼古拉兹（Nikurasde）粗糙管区公式：

2.2 海森-威廉公式

那本书，方玫一直没有丢，走哪带哪。许多年后她在车站包被偷了，里面的东西她全然不在乎，唯独那本书的丢失让她难过了很久。

式中：Ch为海森-威廉系数；Q 为管道流量，m3/s。海森-威廉公式是在D＜1.8m 工业管道的大量测试数据基础上建立的经验公式，适用范围为部分紊流过渡区，Ch 的取值范围介于100 ～157 之间［9-10］，该式在我国市政给排水管网设计中应用较广。

2.3 谢才公式

式中：C 为谢才系数，m1/2/s；R 为水力半径，m。谢才系数C 通常采用曼宁（Manning）公式计算：

式中：n 为曼宁糙率系数。曼宁经验公式最初来源于明渠均匀流的水力计算，用于管道时仅适用于紊流粗糙区或非满管重力流输水。然而，由于一直以来的工程设计习惯，以及公式简单的指数关系，我国水利工程中普遍使用该式进行水力计算，即使管内水流处于紊流过渡区流态［11］。

上述各摩阻系数之间可相互进行换算，其关系式为：

3 水力摩阻系数的敏感性分析

以达西-魏斯巴哈公式和柯尔勃洛克-怀特公式为水力计算标准，依据各摩阻系数间换算关系式（7）和式（8），分析管道内壁当量粗糙度k、管径D、流速V、水温（水流黏滞系数ν）等因素对摩阻系数取值的敏感性影响。

3.1 摩阻系数的相对变化率关系式 上述各水力计算公式可写为水力坡降J 与摩阻系数、管径D 和流速V 的关系式，对水力坡降J 求全微分可得：

式中：J 为水力坡降，J= hf/L；由微分运算性质可知式（9）—式（11）表示水力坡降J 与各水力摩阻系数、流速V 和管径D 之间的相对变化率关系。特别地，当流速V 和管径D 一定时各摩阻系数的相对变化率近似满足：

3.双腔吸虫病。(1)定期驱虫；(2)粪便堆积发酵处理；(3)加强牛群管理,避免到利于中间宿主陆地螺及蚂蚁生存条件的地区放牧。

上式表明：流速V 和管径D 一定时，管道内壁越粗糙，水力坡降J、达西摩阻系数λ和曼宁糙率系数n 越大，而海森-威廉系数Ch越小；曼宁糙率系数n 和海森-威廉系数Ch取值的1%相对误差将引起水力坡降J 或沿程水头损失约2%的相对误差，计算误差被成倍放大，这也进一步说明了摩阻系数合理取值对输水工程精确水力计算的重要性。

3.2 当量粗糙度k 图1为DN1000mm，水温15 ℃，流速范围0.3 ～3.0 m/s（考虑工程运行的技术经济性，输水管道流速通常不超过3.0m/s［1］），当量粗糙度k 分别为0.00、0.10、0.20 和0.30 mm 时各摩阻系数随流速的变化，表1为流速V=1.0 m/s 和V=3.0 m/s 的摩阻系数值及其相对变率。由图表可知：流速V 和管径D 一定时，管道的当量粗糙度k 增大，相对粗糙度（k/D）增大，水力坡降J、达西摩阻系数λ和曼宁糙率n 值增大，海森-威廉系数Ch减小，摩阻系数间的相对变化率关系近似遵循式（12）；各摩阻系数的相对变化率随流速的增大而增大，随当量粗糙度k 的增大而减小。

在教学过程中，教师往往是灌输式教学，在课堂中滔滔不绝，从未考虑学生对物理是否理解.最终教师讲了很多，学生听了很多却没有实际的意义.探究性一部分内容要求的也是学生的自主性，在新的课堂教学中学会转换角色，让学生成为教学的重点，给学生空间去探究和解决物理难题，教师则担任从旁指导的角色，在课堂中引导学生去探究物理，去发现物理，去理解物理.给学生一定的思考空间，让学生通过自己的双手真正意义上的理解物理，吸引学生对物理学科的兴趣.

表1 当量粗糙度k 对摩阻系数的敏感性影响

V/(m/s)ΔJ/J/%λ Δλ/λ/%n Δn/n/%Ch ΔCh/Ch/%1.0 14.06 8.78 6.45 14.06 8.78 6.45 6.80 4.30 3.17-6.86-4.44-3.32 3.0 k/mm 0.00 0.10 0.20 0.30 0.00 0.10 0.20 0.30 J/(m/km)0.608 0.693 0.754 0.803 4.559 5.792 6.478 6.984 27.06 11.83 7.81 0.0119 0.0136 0.0148 0.0158 0.0099 0.0126 0.0141 0.0152 27.06 11.83 7.81 0.0098 0.0105 0.0109 0.0112 0.0089 0.0101 0.0106 0.0111 12.72 5.75 3.83 153.8 143.2 136.8 132.3 155.4 136.6 128.6 123.4-12.13-5.86-3.98

3.3 管径D 图2为当量粗糙度k=0.10 mm，流速范围0.3 ～3.0 m/s，DN100 ～DN2600 不同管径的摩阻系数计算结果。大量试验研究和现场检测结果表明，水泥砂浆内衬球墨铸铁管［12-13］和PCCP 管道［14-15］的当量粗糙度近似为k=0.10 mm。由图可知，当量粗糙度k 一定时，在相同流速V 和水温下，管径D增加，相对粗糙度（k/D）减小，水力坡降J 和达西摩阻系数λ减小，而曼宁糙率系数n 和海森-威廉系数Ch均增加。

曼宁糙率n 随管径D 的变化规律与《SL 702-2015 预应力钢筒混凝土管道技术规范》［16］中“同种管材，管径D 越大，粗糙系数n 越小”的论述相矛盾。由微分关系式（11）可知，水力坡降J 随着曼宁糙率n 和流速V 的增大而增大，但随着管径D 的增大而减小。由于管径D 变化引起的曼宁糙率n 的相对变化率要远小于D 的相对变化率（参见表2），因此，虽然曼宁糙率n 随着管径D 增大而略有增大，但是水力坡降J 总是随着管径D 的增大而减小。只有当管径D 一定时，水力坡降J 才随着曼宁糙率n 的增大而增大。显然，规范［16］对曼宁糙率n 与管径D 之间的变化规律考虑不足。

再次，还深化了对“评价”和“编辑”这两个阶段联系机制的认识。本研究认为，失地农民对“项目质量”“机构实力”“师资力量”的特征感知越积极，对“就业素质提升”这一选项结果的前景价值评价也会越高。鉴于“项目质量”“机构实力”“师资力量”隶属于“编辑”阶段，而对“就业素质提升”的前景价值评价隶属于“评价”阶段，因此，失地农民对“项目质量”“机构实力”和“师资力量”的特征感知，同对“就业素质提升”这一选项结果的前景价值评价之间的关系，所反映的其实就是“评价”和“编辑”这两个阶段的联系机制。

图3为当量粗糙度k=0.10 mm，不同流速、不同管径的曼宁糙率系数n 和海森-威廉系数Ch。由图可知：流速V=1.0 m/s 时，管径由DN100 变化到DN2600，曼宁糙率系数n 从0.0092 增大到0.0112，增幅为22%；海森-威廉系数Ch从134 增大到145，增幅为8%；管径D 对曼宁糙率n 的敏感性影响大于对海森-威廉系数Ch的影响，尤其是当D＞1.0 m 时，n 的增幅为7%，而Ch的增幅仅为1%。工程设计阶段，若不考虑管径D 的影响，对所有管径的曼宁糙率系数n 和海森-威廉系数Ch取固定值，则根据摩阻系数间的相对变化率关系式（12），沿程水头损失计算误差最大可分别达到44%（谢才-曼宁公式）和16%（海森-威廉公式）。

图1 当量粗糙度k 对摩阻系数的敏感性影响

表2 管径D 对水力坡降J 和曼宁糙率n 的敏感性影响

n DN900 DN1000 DN1100 ΔD/D-10.00%0.01037 0.01045 0.01052 Δn/n-0.77%ΔJ/J+13.26%备注k=0.10mm，V=1.0m/s+10.00%+0.67%J/(m/km)0.786 0.694 0.619-10.81%

3.4 流速V 如图2所示，管内水流处于紊流过渡区时，管道流速V（或流量Q）增大，相同管径D 和水温下的雷诺数Re 增大，相对粗糙度（k/D）一定时，水力坡降J 增大，达西摩阻系数λ和曼宁糙率n 减小，海森-威廉系数Ch则先增大后减小。如图3，流速从1.0 m/s 增大到2.0 m/s，曼宁糙率系数n 减小约2.71%，海森-威廉系数Ch减小约2.56%。当水流进入紊流粗糙区后，达西摩阻系数λ仅与管道的相对粗糙度（k/D）有关，此时达西摩阻系数λ和曼宁糙率n 为一常数。

3.5 水温T（运动黏滞系数ν） 表3为DN1000 mm，当量粗糙度k=0.10 mm，水温分别为10 ℃（ν=1.310×10-6m2/s）、15℃（ν=1.145×10-6m2/s）、20℃（ν=1.009×10-6m2/s），流速V=1.0 m/s 和V=3.0 m/s 时各摩阻系数值及其相对变率。由表可知，管道相对粗糙度（k/D）一定时，相同流速V 和管径D，温度T 升高，运动黏滞系数ν减小，水流雷诺数Re 增大，水力坡降J、达西摩阻系数λ和曼宁糙率n 均减小，海森-威廉系数Ch增大，各摩阻系数的相对变化率关系近似满足式（12），其相对变化幅值随流速V 增大而减小。

图2 管径D 对摩阻系数的敏感性影响

图3 管径D 和流速V 对摩阻系数的敏感性影响

表3 水温T 对摩阻系数的敏感性影响

V/(m/s)ΔJ/J/%λ Δλ/λ/%n Δn/n/%Ch ΔCh/Ch/%1.0-1.14-1.15-1.23-1.18-0.67-0.57 0.70 0.63 3.0 T/℃10 15 20 10 15 20 J/(m/km)0.702 0.694 0.686 5.830 5.792 5.760-0.65-0.55 0.01378 0.01361 0.01345 0.01271 0.01263 0.01256-0.63-0.55 0.01052 0.01045 0.01039 0.01010 0.01007 0.01004-0.30-0.30 142.2 143.2 144.1 136.1 136.6 137.0 0.37 0.29

3.6 摩阻系数的敏感性变化规律 基于柯尔勃洛克-怀特公式和各摩阻系数间换算关系，得出管道摩阻系数随当量粗糙度k、管径D、流速V 和水温T 的敏感性变化规律，如表4。对于同种管材（假定其当量粗糙度k 为一定值），管径D 的取值变化范围最大，其对摩阻系数的敏感性影响也最大，其次为流速V 的影响，水温T 的影响最小。

表4 摩阻系数随各敏感性因素的变化规律

影响因素k↗D↗V↗T↗J↗↘↗↘λ↗↘↘↘n↗↗↘↘Ch↘↗↘↗

4 实例分析与工程应用

4.1 实例分析 表5为典型工程曼宁糙率系数n 和海森-威廉系数Ch的取值比较，表中实测值为工程现场测量值，设计值为工程设计取值，理论值为采用柯尔勃洛克-怀特公式求得达西摩阻系数λ后依据摩阻系数间关系式（7）和式（8）换算得到。由表可知：现有工程设计中曼宁糙率n 取值的随意性较大，考虑管径D 和流速V 影响的理论计算n 值与实测值接近，计算精度有较大地提高；海森-威廉系数Ch的经验取值较为成熟，海森-威廉公式比谢才-曼宁公式更适宜于管道输水工程的水力计算。

表5 典型工程案例分析［5-7］

注：（1）不同管材的当量粗糙度k 采用文献［6］的值；（2）玻璃钢管的海森-威廉系数采用《M45 玻璃钢管设计》给出的Ch=150；（3）PCCP 管的海森-威廉系数采用《M9 混凝土压力管》给出的1.37m＜D＜2.89 时Ch=145，D＞2.89 时Ch =150。

工程名管材D/mm V/(m/s)k/mm n Ch实测值设计值设计误差理论值理论误差实测值规范值规范误差理论值理论误差小洼槽倒虹吸玻璃钢管3100 2.252 0.01 0.0106 0.0090-15.09%0.0102-3.77%146 150 2.74%152 4.11%三个泉倒虹吸环氧内衬钢管2700 2.663 0.01 0.0098 0.0120 22.45%0.0100 2.04%155 150-3.23%152-1.94%三个泉倒虹吸PCCP 管2800 2.477 0.10 0.0108 0.0135 25.00%0.0109 0.93%141 145 2.84%139-1.42%南水北调北京段PCCP 管4000 1.592 0.10 0.0101 0.0120 18.81%0.0114 12.87%162 150-7.41%142-12.35%

4.2 水力计算公式的选用 对于不同管材、不同管径的水力计算应优先采用达西-魏斯巴哈公式和柯尔勃洛克-怀特公式，这也是国内外规范推荐的主要公式，其计算精度较高，适用于紊流的所有流态，并在科学研究和工程实践中得到广泛验证。

海森-威廉公式适用于常温下清水输送管道的紊流过渡区，计算较为简便快捷，计算精度可满足工程要求。海森-威廉系数Ch的取值需考虑管径D 的影响，美国水行业协会（AWWA）技术手册《M9 混凝土压力管》［17］给出Ch=139.3+6.65D，《M11 钢管》［18］给出Ch=130.0+ 6.30D。图3（b）表明管径D 大于1.0 m 后，海森-威廉系数Ch随管径D 的变化趋于平缓，若采用上述规范给出的线性拟合公式，其误差将随着管径D 的增大而增大。对于k=0.10 mm、V=1.0 m/s 的水泥砂浆内衬管道其经验拟合关系式为Ch=147.25-4.13/D0.5。

随着制管工艺和水平的发展，现有大型输水工程的管道内壁均匀、光滑，水流多处于紊流过渡区，流态超出谢才-曼宁公式的适用范围，并且现有规范规定的曼宁糙率n 取值范围较宽，如规范［16］指出立式震捣法制造的PCCP 管的曼宁糙率系数n 取0.010 ～0.0125，变幅为25%，相应水头损失的变幅将达到50%。由于未考虑管径D 和流速V 的影响，造成实际工程中n 的精确取值困难，甚至较大的水头损失计算误差。

特别地，当α=0，即每个节点的权值为1，和f=0，即网络为无向时，节点上的负载就等于文献[6,8,14]中被广泛使用的介数.

4.3 摩阻系数取值的裕度 管道摩阻系数的取值不仅要考虑内壁绝对粗糙度的影响，而且需要考虑弯管、接头、阀门等的局部水头损失影响，以及长期运行后由于化学和生物作用产生的腐蚀、沉淀、结垢等对管壁粗糙度的影响［19］。实践中需综合考虑上述因素，结合工程经验，对摩阻系数的选取预留适当的裕度，以保证输水工程在长期运行过程能够达到设计流量［20］。考虑裕度的摩阻系数取值可参考关系式（12），假设工程设计中管线的水头损失预留10%的裕度，则达西摩阻系数λ需增大约10%，曼宁糙率n 需增大约5%，海森-威廉系数Ch需减小约5.4%。

5 结论

随着经济和社会发展，大口径、新材料管道在水利工程中得到广泛应用，本文对影响管道输水能力的关键参数——水力摩阻系数进行了系统的敏感性分析，得出如下结论：

（1）管道水力计算应优先采用达西-魏斯巴哈公式和柯尔勃洛克-怀特公式，尽管后者为隐式方程，需迭代求解，但在当前计算机技术发展的条件下，这不应成为该公式在工程中普遍使用的障碍。

所有患者均仰卧于体部的真空垫进行CT定位扫描及接受射波刀(Accuray Incorporated, Sunnyvale, CA))放射治疗。肿瘤靶区(gross tumor volume, GTV)定义为影像学所示病灶靶区，计划靶区(planning target volume, PTV)定义为肿瘤靶区在x、y、z轴方向各外放2～3 mm。总剂量30～36 Gy(5～8次分割)。至少90% PTV被处方剂量线包绕。

（2）海森-威廉公式和曼宁公式均为经验公式，由于误差传递，海森-威廉系数Ch和曼宁糙率系数n 的取值误差将引起沿程水头损失计算误差的成倍放大，因此，二者的合理选取非常重要。海森-威廉系数Ch和曼宁糙率系数n 的经验取值不仅要考虑管道内壁粗糙度的影响，流速V 和管径D 的影响也不能忽略，对于同种管材，二者均随管径D 增大而增大，随流速V 增大而减小。相比较而言，海森-威廉系数Ch对不同管径D 的适应性优于曼宁糙率系数n，例如当管径D＞1.0m 时，n 的变幅为7%，而Ch的变幅仅为1%。

（3）工程设计应慎用谢才-曼宁公式。现有规范给出的曼宁糙率系数n 取值通常不考虑管径D 和流速V 的影响，且取值范围较宽，n 的精确取值困难，甚至可能影响工程效率和效益的发挥。

参 考 文 献：

［1］GB 50013-2006 .室外给水设计规范［S］.

［2］王雪原，黄慎勇，付忠志.长距离输水管道水力计算公式的选用［J］.给水排水，2006，32（10）：32-35 .

［3］KAMAND F Z .Hydraulic friction factors for pipe flow［J］.Journal of Irrigation & Drainage Engineering，1988，114（2）：311-323 .

［4］BOMBARDELLI F A，GARCIA M H .Hydraulic design of large diameter pips［J］.Journal of Hydraulic Engineering，2003，129（11）：839-846 .

［5］蒲振旗，徐元禄，周骞.玻璃钢管道糙率值实证分析［C］//中国水利学会第五届青年科技论坛论文集.2012：288-295 .

［6］李新，谢晓勇.大型倒虹吸工程水头损失及水力计算［J］.人民长江，2017，48（20）：71-75 .

［7］贾天云.北京市南水北调工程PCCP 输水管道糙率分析［C］//中国水利学会2015 学术年会，2015：577-581 .

［8］COLEBROOK C F，WHITE C M .Experiments with fluid friction in roughened pipes［J］.Proceedings of the Royal Society A-Mathematical，Physical and Engineering Sciences，1937，161（904）：367-381 .

［9］DISKIN M H .The limits of applicability of the Hazen-Williams formula［C］//Houille Blanche，1960（6）：720-723 .

［10］TRAVIS Q B，MAYS L W .Relationship between Hazen-William and Colebrook-White Roughness Values［J］.Journal of Hydraulic Engineering，ASCE，2007，133（11）：1270-1273 .

［11］刘进，王东黎，杨进新.南水北调中线PCCP 管道的摩阻损失计算分析［J］.南水北调与水利科技，2008，（1）：231-234 .

［12］久保田铁工株式会社.球墨铸铁管手册［M］.北京：中国金属学会铸铁管委员会，1994 .

［13］杨开林，郭永鑫，付辉，等.管道当量粗糙度的率定及不确定度［J］.水利学报，2012，43（12）：1397-1404 .

［14］常胜，牧振伟，李新.大口径预应力钢筒混凝土管当量粗糙度取值探究［J］.水电能源科学，2016，34（5）：124-127，123 .

［15］郑双凌，马吉明，南春子，等.预应力钢筒混凝土管（PCCP）的阻力系数与粗糙度研究［J］.水力发电学报，2012（3）：126-130 .

［16］SL 702-2015 .预应力钢筒混凝土管道技术规范［S］.

［17］AWWA Manual M9，Concrete Pressure Pipe［S］.American Water Works Association，2008 .

［18］AWWA Manual M11，Steel Pipe［S］.American Water Works Association Manual，2004 .

［19］李莎.供水管道管壁腐蚀及摩阻系数研究［D］.重庆：重庆大学，2013 .

［20］GB50316-2008 .工业金属管道设计规范［S］.