乌江渡水电站产汇流分析及应用

ＨｏｎｇＳｈｕｉＲｉｖｅｒＶｏｌ．３８Ｎｏ．５Ｏｃｔ．２０１９

乌江渡水电站产汇流分析及应用

何　凡

（贵州乌江水电开发有限公司乌江渡发电厂，贵州　遵义　５６３１０４）

摘　要：乌江渡水电站位于乌江中游，为了进一步了解乌江流域产流机制，加强应对突发意外情况时对流量变化的掌控，做好乌江渡水电站防洪预报工作，笔者从产流机制、产流方法、汇流机理和汇流方法对乌江流域进行深入的分析，将蓄满产流模型应用在乌江渡水库上，借用ＶＢ编程工具计算入库流量。举例计算了２０１８年一场降雨产生的径流，结果与实测数据对比，误差在可接受范围内，表明此编程模型适用于乌江渡水电站。关键词：乌江渡水电站；产流；汇流；蓄满产流；超渗产流中图分类号：ＴＶ１２１．７

文献标识码：Ａ

文章编号：１００１－４０８Ｘ（２０１９）０５－００１０－０５

１　概述

１．１　流域概况

流域地处湿润的亚热带季风气候区，冬季主要受北方西伯利亚冷空气影响，夏季受印度洋孟加拉湾的西南暖湿气流和西太平洋的海洋性气候影响，温润多雨，同时也因地处云贵高原，造成地势较高的高原面与地势较低洼的河谷地带的气候有明显差异。

降雨年际间变化大，年内分配不均，多年平均降雨量为１４００～１６００ｍ，流域大暴雨多发生在６

Ｉｍ（流域最大初损）ＵＭ（不透水面积比例）

０．０１４０ｍｍ

ＦＣ（稳定下渗率）ＷＵＭ（上层最大缺水量）

２０ｍｍ０．５ｍｍ／ｈ

月—７月，持续时间为１～３天，昼夜温差大。７月降雨量占年降雨量的２０．５％左右。降雨在时空分布很不均匀，西部为少雨地区，雨量不足１０００ｍｍ，东部雨量相对较多，降雨在地区上分布呈自下游向上游递减趋势，乌江渡平均年降雨量９８２．１ｍｍ。流域控制面积为２７７９０ｋｍ２，共有１５个雨量站且间大于白天。１．２　流域基本资料所示。

流域面积为２７７９０ｋｍ２。流域其他资料如表１

Ｋ（折减系数）ＷＬＭ（下层最大缺水量）

４０ｍｍ０．８

ＥＫ（蒸发系数）Ｂ（蓄水容量曲线系数）

０．２５０．６

从索风营到乌江渡区间集水面积为５９２８ｋｍ２，

分布较均匀。由于受到地形影响，暴雨强度一般夜

表１　流域基础资料表

２　乌江流域产流机制

从降水（暴雨）到产生水流，接着水流汇至流域出口断面的整个过程被称为产汇流过程。这一动态物理过程错综复杂：一部分雨水以径流的形式汇集到流域出口，另一部分经由蒸散发返回大气，通常把不转化为径流的那一部分雨水称为损失。一次降雨过程的总雨量扣去损失得到径流量，也称为净雨量，这一过程就是产流［１］，如图１所示。

　　本文基于乌江流域现有资料，对比分析蓄满和超渗条件，进而确定乌江渡水电站所在位置的产流机制。

图１　产流过程图

２．１　超渗产流和蓄满产流２．１．１　概念

径流是在降雨、植物截留、填洼、蒸发及下渗等

　　收稿日期：２０１９－０４－２５；修回日期：２０１９－０６－０６

　　作者简介：何　凡（１９９４），女，四川南充人，工学学士，研究方向：水文水资源工程专业，Ｅ－ｍａｉｌ：９０８０６２０８３＠ｑｑ．ｃｏｍ。

１０

何　凡：乌江渡水电站产汇流分析及应用　

几个过程组合下产生的。

１）超渗产流：指地面径流产生的原因是同期降水量大于同期植物截留量、填洼量、雨期蒸发量及下渗量等的总和，多余出来的水量产生了地面径流。植物截留量、雨期蒸发量、填洼量一般较小，而下渗量一般较大，下渗量的时空变化一般表现为：同一种土壤情况下，土壤干燥，下渗能力强；土壤湿润，下渗能力小。产生下渗的条件是产流界面为地面（包气带的上界面）、有水源（降水）、降雨强度大于下渗能力。

２）蓄满产流：又称超蓄产流。因降水使土壤包气带和饱水带基本饱和而产生径流的方式，是降雨径流的产流方式之一。在降雨量较充沛的湿润、半湿润地区，地下潜水位较高，土壤前期含水量大，当一次降雨量大，历时长，降水满足植物截留、入渗、填洼损失后，损失不再随降雨延续而显著增加，土壤基本饱和，从而广泛产生地表径流［２］。

具体定性参数如表２所示。

表２　产流方式对比表

编号１２３４５６７８９

对比分析内容多年平均降雨量多年平均径流系数降雨强度

蓄满产流＞１０００ｍｍ＞０．４大小

超渗产流＜４００ｍｍ＜０．２小大

　　４）地下径流：流域内富含大量地下水，地下水径流占整体径流比例较大。２．２．２　定量分析

１）径流系数：同一流域面积、同一时段内径流深度（Ｒ）与降水量（Ｐ）的比值称为径流系数，以小数或百分数计，表示降水量中形成径流的比例，它综合反映了流域内自然地理要素对径流的影响。计算式为：α＝Ｒ／Ｐ，式中α为径流系数，Ｒ为径流深度，Ｐ为降水深度。α值在０～１之间，湿润地区α值大，干旱地区α值小。

乌江渡年径流系数＝乌江渡坝上断面年径流深／乌江渡坝上断面年降雨量。运用集控系统查取３年数据如表３所示。

表３　乌江渡实测资料表

２０１５年２４６１０６９４４０４１５．１６９５７．８０．４３３４５１６６

２０１６年１７８０８３７９２０３００．４１９０９．５０．３３０３０２３６４

参数年径流量／ｍ３年径流深／ｍｍ年降雨量／ｍｍ年径流系数

２０１７年３７２．４４１８６２３９１６．６０．４０６３２９７６５２２０７８３５３６０

２．１．２　定性参数

　　表２中第２行说明年径流系数大于０．４则为蓄满产流。通过表３计算结果可以看出，索风营到乌江渡区间年径流系数大于０．４，满足蓄满产流定义。

２）降雨量与流量关系：查取区间降雨以及区间断面流量做出２０１５年的降雨量与流量关系图，如图２所示。

流量过程线不对称系数影响产流因素表层土质结构缺水量地下径流

初始土湿度和降雨量疏松、不易超渗小、易蓄满比例大

初始土湿度和降雨强度密实、易超渗大、不易蓄满比例小与降雨强度关系密切

产流与降雨特征关系

与降雨量关系密切

２．２　乌江渡区域流域产流方式论证２．２．１　定性分析

１）多年平均降雨量：查阅资料表明流域多年２）气候：乌江流域地处湿润地区，气候暖３）表层土质结构：流域内土质疏松，透水性根据表２进行定性分析：

图２　２０１５年降雨量－流量相关图

平均降雨量为１４００～１６００ｍｍ，大于１０００ｍｍ，满足蓄满产流条件。

　　由图２中可看出，流量与降雨量关系密切，降雨多则流量大。且一次洪水的流量过程陡涨缓落，持续时间较长，满足蓄满产流条件。

３）统计１９８８年—１９９１年实测流域降雨量和径流深，算出径流系数，并对其进行相关性分析，结果如表４所示。图３表示降雨量与径流深两者的相关系数Ｒ２＝０．６９６７。由表２第９行可看出，产生的径流量与降雨量关系越密切则表示其是蓄满产流，计算结果相关系数大于０．５，则表示径流量与降雨量为强相关，满足蓄满产流定义。

６月１７１．４６６．０３０．３９

７月８３．９２８．１８０．３４

８月１６４．１６２．２４０．３８

９月８６．６２６．６００．３１

１０月７０４３．２００．６２

１１月１５．１１７．４４１．１５

１２月３０．５０．９１０．０３

一年９５７．８３８５．８０．４０

热，雨量充沛，汛期为５月—８月。

好，植被种类多且良好，不易发生超渗。

分类降雨量／ｍｍ

２０１５年

径流深／ｍｍ径流系数

１月３３．９１３．４００．４０

２月１９．２８．８４０．４６

表４　各年各月实测降雨径流表

３月２８．８８．０１０．２８

４月７５．６１４．５３０．１９

５月１７８．７９６．３６０．５４

１１

　红水河２０１９年第５期

表４（续）

分类降雨量／ｍｍ

２０１６年

径流深／ｍｍ径流系数降雨量／ｍｍ

２０１７年

径流深／ｍｍ径流系数

１月３１．８１１．９６０．３８１９．５４．４４０．２３

２月１２．２８．７１０．７１１７．７５．８９０．３３

３月７２．２２０．３３０．２８３８１０．９８０．２９

４月１３１．５５２．１３０．４０４６．６８．６６０．１９

５月６１．３２７．１９０．４４６４．３１９．４２０．３０

６月２３０．６５１．４１０．２２２８３．２１４５．７２０．５１

７月９３９７．２２１．０５８７．２８０．６１０．９２

８月１３１．４３０．０４０．２３１０９．３３２．７５０．３０

９月１９．１８．４５０．４４１３３．５５２．２５０．３９

１０月６６．２７．７８０．１２８３．９２９．１２０．３５

１１月４５．２９．６２０．２１１５．８１０．０００．６３

１２月１５．１６．２７０．４２１７．５１１．０９０．６３

一年９０９．６３３１．１２０．３６９１６．５４１０．９４０．４５

图３　降雨径流相关性分析图

　的气象条件　总结：不论是定量分析还是定性分析、下垫面条件和流量过程来看，，从流域该流域降雨径流关系具有蓄满产流的特点，按蓄满产流建立产流量预报方案。２．３　产流计算公式２．３．１　降雨量计算

面雨量的计算方法为加权平均法：在对流域基本情况（如面积、地类、坡度、坡向、海拔等）进行勘察的基础上，在每个地类上选择有代表性的地点作为降水观测点，把每个测点控制的地类面积作为各测点降水量的权重［３］。

１５由现有资料得知：该地区雨量站分布均匀，有

泰森加权平均法个雨量站，乌江渡电站在修建时面雨量计算选择。１５个雨量站权重系数如表５所示。

表５　１５个雨量站的权重系数表

站名赫章对江洪家渡铁索桥鸭池河狮石子火石坝息烽权重０．１４６０．０９６０．０６４０．０９８０．０２８０．０１７０．０４９０．０３０站名乌江渡徐家渡

织金金沙毕节水城纳雍权重０．０６９

０．０３８

０．０６０

０．０７０

０．０８０

０．０４５

０．１０４

２．３．２　蒸发量计算

一般来讲，流域蒸发量主要取决于土壤蒸发量。在一定的气象条件下，裸土含水量从饱和到干燥的蒸发过程，大致呈现为图４所示的３个阶段特征。从图４可看出，蒸发随土壤含水率θ变化关系可概化为：　　当θ≥θａ时，

ＥＥＳ

＝Ｐ１．０；　　当θＥＥＳ

ａ≥θ≥θｂ时，

＝Ｐｆ（θ）＝αθ；

１２

　　当θ＜θＥｂ时，

ＥＳ

＝Ｐ

Ｃ；

式中：θａ、θｂ为第一、第二临界点土壤含水率，ＥＳ为土壤蒸发量，ＥＰ为流域蒸发能力，Ｃ为蒸发扩散

系数，α为系数。

图４　蒸发过程图

　壤含水率不宜用于水量平衡式的产流量计算　上式表示了蒸发与土壤含水率的关系。由于土，把蒸发与土壤含水率的关系转化为与土壤含水量的关系。

三层蒸发模式按照先上层后下层的次序，具体计算为　　１）：

当ＷＵ＋Ｐ＞ＥＰ时：

ＥＵ＝ＥＰ，ＥＬ＝Ｅ２）当ＷＵ＋Ｐ＜Ｅ０，ＥＤ＝０

Ｐ，ＷＬ≥Ｃ∙ＷＬＭ时Ｕ＝ＷＵ＋Ｐ，ＥＬ＝（ＥＰ－ＥＵ）ＷＬ

，：

ＷＬＭ３）当ＷＵ＋Ｐ＜ＥＷＬＭ

ＥＤ＝０

Ｐ，Ｃ（ＥＰ－ＥＥ时：

Ｕ）≤ＷＬ＜Ｃ∙Ｕ＝ＷＵ＋Ｐ，ＥＬ＝ＣＥ４）ＷＵ当＋ＷＵＰ，＋ＥＰ＜Ｅ（ＷＬＥＰ－＜ＥＣＵ（），ＥＥＤ＝０Ｐ，Ｐ－ＥＵ）时：Ｕ＝Ｌ＝ＷＬ，ＥＤ＝Ｃ（ＥＰ－ＥＵＰ为降雨量，Ｅ）－ＥＬ式中：ＫＰ为流域蒸发能力（ＥＰ＝Ｅ×

ＥＣ，Ｅ为水面蒸发，ＫＣ为蒸发折算系数），ＥＵ深层蒸发量，ＷＵ、ＷＬ、、ＥＬＷＤ、Ｄ分别为上层、下层、分别为上层、下层、深层土壤含水量，ＷＬＭ为下层土壤含水容量，Ｃ为蒸发扩散系数［４］２．３．３　降雨产流量计算

。

蓄满产流，即任一地点上，土壤含水量达到田间持水量前，降雨量全部补充土壤含水量，不发生产流；当土壤蓄满后，其后续降雨量全部产生径流。具体流程如图５所示。

ａ＝ＷＭＭæçＷè１－（１－ＷＭ）ｂ１＋１ö÷ø

何　凡：乌江渡水电站产汇流分析及应用　

时，常常将壤中流进一步划分为两部分，快速部分和地面径流合在一起，成为直接径流；而慢速部分则和地下径流合并，统称为地下径流。３．１．２　计算公式

通过稳定下渗率ＦＣ可划分产流中的直接径流和地下径流。

图５　产流过程图

发生产流时：

１）当ＰＥ＋ａ＜ＷＭＭ时，Ｒ＝ＰＥ＋Ｗ－ＷＭ＋

（ＰＥ＋ａ）（ｂ＋１）

ＷＭ（１－）；

ＷＭＭ

２）当ＰＥ＋ａ≥ＷＭＭ时，Ｒ＝ＰＥ＋Ｗ－ＷＭ；其中：ａ为初始土壤含水量，ＷＭ为流域平均需水容量，ＷＭＭ为流域最大蓄水容量，Ｗ为土壤含水量，ｂ为流域土壤含水量分布指数，ＰＥ为净雨，Ｒ为产流量。

　　当ＰＥ＞０时，发生产流；否则不发生产流。

　　１）当ＰＥ≤ＦＣ时，ＲＤ＝０，ＲＧ＝Ｒ；

２）当ＰＥ＞ＦＣ时，ＲＧ＝ＦＣ

Ｒ

，ＲＤ＝Ｒ－ＲＧ；ＰＥ

式中：ＦＣ为稳定下渗率，ＲＤ为直接径流，ＲＧ为地下径流。

次洪的各水源分量为：　　总地下径流：ＲＧ＝　　总直接径流：ＲＳ＝３．２　卷积差分法

差分形式算出流域出口的流量过程［５］。计算公式为：　　直接径流：ＱＳ（Ｉ）＝ＲＳ（Ｉ）∙Ｕ；

根据流域净雨和流域径流单位线，采用卷积的

ｆｃ∑ＰＥｉ

ｒｉ

ｉ

＋

ＰＥｉ＞ＦＣ

ＰＥｉ＞ＦＣ

（ＰＥｉ∑ｉ

ＰＥｉ≤ＦＣ

ｒｉ；∑ｉ

ＰＥｉｒｉ

。

－ｆｃ）

３　汇流

流域内各部分径流向流域出口断面的汇集过程，称为流域汇流（见图６）。通常可以把流域分成坡地和河网两个基本部分，因此流域汇流也可以分为坡地汇流与河网汇流两部分。

地下径流：ＱＧ（Ｉ）＝ＣＧ∙ＲＧ（Ｉ）∙Ｕ＋（１－ＣＧ）∙

ＱＧ（Ｉ－１）；总流量：Ｑ（Ｉ）＝ＱＳ（Ｉ）＋ＱＧ（Ｉ）；

式中：ＱＳ为地表径流流量，第一个Ｕ为单位线纵坐标，第二个Ｕ为折算系数，ＱＧ为地下径流流量，ＣＧ为Δｔ／（ｋｇ＋０．５Δｔ）（其中，Δｔ为计算时段，ｋｇ为蓄泄系数），Ｑ为总流量。

４　编程计算

图６　汇流过程图

３．１　划分水源３．１．１　原理

地面径流由坡面直接汇入河网，运动速度快，形成的流量过程呈陡涨陡落，是涨洪和洪峰附近流量过程的主体部分；地下径流是由渗透到浅水面的水流缓慢流出，运动速度慢，汇流时间长，变化平缓，是洪水退水尾部段流量的主要成分，且常延续至后继洪水过程中；而壤中流介于两者之间。在计算

根据上述原理，降雨到径流过程复杂，数据计算量大，选择用ＶＢ建立模型编程计算。具体程序见有关材料。

５　应用与结论

５．１　应用

利用所编程序计算２０１８年５月的一场降雨，原始数据如表６所示，单位线如表７所示，计算结果如表８所示。

表６　２０１８年５月一场降雨数据表

项目降雨量／ｍｍ索风营出库流量／（ｍ３／ｓ）

５月１９日５月１９日５月１９日５月２０日５月２０日５月２０日５月２０日５月２０日５月２０日５月２０日５月２０日５月２０日２１：００２２：００２３：００００：０００１：０００２：０００３：０００４：０００５：０００６：０００７：０００８：００１．３６７６６５．６４３

４．０６５７８９．７１３

４．９５７８２２．５７

２．７８２３４５．６３２

３．２９０．３８１

２．２６９０．２０４

０．９０６０

０．０２７０

００

０．１１４０

００

０．３５８１３０．８９９

１３

　红水河２０１９年第５期

表７　单位线表

时段（６ｈ）

００

１４０

２５３５

３６７３

４５４８

５４１３

６３１３

７２４０

８１８１

９１３７

１０１０２

１１７２

１２４７

１３２８

表８　计算结果数据表

项目降雨量／ｍｍ索风营出库流量／（ｍ３／ｓ）计算入库流量／（ｍ３／ｓ）

５月１９日５月１９日５月１９日５月２０日５月２０日５月２０日５月２０日５月２０日５月２０日５月２０日５月２０日５月２０日２１：００２２：００２３：００００：０００１：０００２：０００３：０００４：０００５：０００６：０００７：０００８：００１．３７６６５．６４６６５．６４

４．０７７８９．７１７８９．７１７８７．９４

４．９６８２２．５７９３２．１０７６３．２８

２．７８３４５．６３７５８．３８８２９．５７

３．２９０．３８７９７．９０８８４．１６

２．２７０．２０６７８．４３７４６．０３

０．９１０．００４５１．６０３６１．５６

０．０３０．００５７．２０１．９３

０．０００．００４０．３２３６．４３

０．１１０．００２１．５０１６９．８９

０．０００．００１．４４５０．２３

０．３６１３０．９０１３０．９０１２７．８５

系统显示入

７１８．６４

库流量／（ｍ３／ｓ）

５．２　误差分析５．２．１　模型结构误差

预报方案采用的预报方法或水文预报模型都是对实际流域进行概化后建立的，例如将空间分布的水文气象要素表示成流域上平均值的处理；将水文循环过程进行人为分离，并用不同程度简化的模型予以描述等，这样不完善性处理引起的误差称为模型结构误差。５．２．２　模型参数误差

由于本文所用参数全来源于大坝年鉴资料，时间久远，甚至很多用的是大坝修建前测量的数据，不具有代表性，必然会存在误差，如若想要得到精确的预报，还需重新定参数。５．２．３　下垫面条件改变造成误差

由于年代久远和上游修建水库的影响，大坝上游河道下垫面条件早就发生改变，造成误差。５．２．４　模型输入资料误差

进行水文预报所输入的资料通常是降雨、流量和蒸发量，这些资料或由实测获得，或根据天气预

报估算得到。前者存在着测验和时段统计平均误差，后者则存在着相当大的预报误差。５．３　结论

通过实际计算，误差在要求之内，这表明蓄满产流、卷积法适用于乌江渡水电站的实际计算，遇到特殊情况（例如：入库流量数据异常、停电、断网等不可预测情况）时，可以通过本文介绍的方法计算入库流量，及时掌握最新信息，第一时间发现险情，预报水位涨势。

参考文献：

［１］　赵人俊．流域水文模拟———新安江模型与陕北模型［２］　崔泰昌，陆建华．试论蓄满产流模型与超渗产流模型

［Ｊ］．山西水利科技，２０００（３）：１３－１５．［３］　王安琪，瞿宁玲，陈元芳，等．上海市杨浦区设计面暴雨［４］　包为民．水文预报［Ｍ］．３版．北京：中国水利水电出版

社，２００６：１５４－１５５．

量与设计暴雨雨型研究［Ｊ］．水电能源科学，２０１８，３６（８）：１－４．

［Ｍ］．北京：水利电力出版社，１９８４．

［５］　董向东．流域汇流单位线洪水预报应用分析［Ｊ］．科技创

新与应用，２０１９（４）：１６６－１６７．

ＡｎａｌｙｓｉｓａｎｄＡｐｐｌｉｃａｔｉｏｎｏｆＲｕｎｏｆｆＧｅｎｅｒａｔｉｏｎａｎｄＦｌｏｗＣｏｎｃｅｎｔｒａｔｉｏｎ

ｉｎＷｕｊｉａｎｇｄｕＨｙｄｒｏｐｏｗｅｒＳｔａｔｉｏｎ

ＨＥＦａｎ

ＷｕｊｉａｎｇｄｕＰｏｗｅｒＰｌａｎｔｏｆＧｕｉｚｈｏｕＷｕｊｉａｎｇＨｙｄｒｏｐｏｗｅｒＤｅｖｅｌｏｐｍｅｎｔＣｏ． Ｌｔｄ． Ｚｕｎｙｉ Ｇｕｉｚｈｏｕ ５６３１０４ Ａｂｓｔｒａｃｔ ＷｕｊｉａｎｇｄｕＨｙｄｒｏｐｏｗｅｒＳｔａｔｉｏｎｉｓｌｏｃａｔｅｄｉｎｔｈｅｍｉｄｄｌｅｒｅａｃｈｅｓｏｆＷｕｊｉａｎｇＲｉｖｅｒ．ＩｎｏｒｄｅｒｔｏｆｕｒｔｈｅｒｕｎｄｅｒｓｔａｎｄｔｈｅｍｅｃｈａｎｉｓｍｏｆｒｕｎｏｆｆｇｅｎｅｒａｔｉｏｎｉｎＷｕｊｉａｎｇＲｉｖｅｒＢａｓｉｎ ｓｔｒｅｎｇｔｈｅｎｔｈｅｃｏｎｔｒｏｌｏｆｆｌｏｗｃｈａｎｇｅｉｎｒｅｓｐｏｎｓｅｔｏｕｎｅｘｐｅｃｔｅｄｓｉｔｕａｔｉｏｎｓ ａｎｄｄｏａｇｏｏｄｊｏｂｉｎｆｌｏｏｄｃｏｎｔｒｏｌａｎｄｆｏｒｅｃａｓｔｉｎｇｏｆＷｕｊｉａｎｇｄｕＨｙｄｒｏｐｏｗｅｒＳｔａｔｉｏｎ ｔｈｅａｕｔｈｏｒｍａｋｅｓａｔｈｏｒｏｕｇｈａｎａｌｙｓｉｓｏｆｔｈｅｒｕｎｏｆｆｇｅｎｅｒａｔｉｏｎｍｅｃｈａｎｉｓｍ ｒｕｎｏｆｆｇｅｎｅｒａｔｉｏｎｍｅｔｈｏｄ ｆｌｏｗｃｏｎｃｅｎｔｒａｔｉｏｎｍｅｃｈａｎｉｓｍａｎｄｆｌｏｗｃｏｎｃｅｎｔｒａｔｉｏｎｍｅｔｈｏｄｉｎＷｕｊｉａｎｇＲｉｖｅｒＢａｓｉｎ ａｎｄａｐｐｌｉｅｓｔｈｅｒｕｎｏｆｆｇｅｎｅｒａｔｉｏｎｕｎｄｅｒｓａｔｕｒａｔｅｄｃｏｎｄｉｔｉｏｎｍｏｄｅｌｔｏＷｕｊｉａｎｇｄｕＨｙｄｒｏｐｏｗｅｒＳｔａｔｉｏｎ ａｎｄｔｈｅｒｅｓｅｒｖｏｉｒｉｎｆｌｏｗｉｓｃａｌｃｕｌａｔｅｄｂｙｕｓｉｎｇＶＢｐｒｏｇｒａｍｍｉｎｇｔｏｏｌ．Ａｎｅｘａｍｐｌｅｉｓｇｉｖｅｎｔｏｃａｌｃｕｌａｔｅｔｈｅｒｕｎｏｆｆｇｅｎｅｒａｔｅｄｂｙａｒａｉｎｆａｌｌｉｎ２０１８ ｔｈｅｃｏｍｐａｒｉｓｏｎｂｅｔｗｅｅｎｔｈｅｃａｌｃｕｌａｔｅｄｒｅｓｕｌｔｓａｎｄｔｈｅｍｅａｓｕｒｅｄｄａｔａｓｈｏｗｓｔｈａｔｔｈｅｅｒｒｏｒｉｓａｃｃｅｐｔａｂｌｅ ａｎｄｔｈｅｐｒｏｇｒａｍｍｉｎｇｍｏｄｅｌｉｓｓｕｉｔａｂｌｅｆｏｒＷｕｊｉａｎｇｄｕＨｙｄｒｏｐｏｗｅｒＳｔａｔｉｏｎ．

Ｋｅｙｗｏｒｄｓ ＷｕｊｉａｎｇｄｕＨｙｄｒｏｐｏｗｅｒＳｔａｔｉｏｎ ｒｕｎｏｆｆｇｅｎｅｒａｔｉｏｎ ｆｌｏｗｃｏｎｃｅｎｔｒａｔｉｏｎ ｒｕｎｏｆｆｇｅｎｅｒａｔｉｏｎｕｎｄｅｒｓａｔｕｒａｔｅｄｃｏｎｄｉｔｉｏｎ ｒｕｎｏｆｆｇｅｎｅｒａｔｉｏｎｆｒｏｍｅｘｃｅｓｓｒａｉｎ１４