发明名称 一种梯级泵站输水系统日优化调控方法
摘要 本发明涉及一种梯级泵站输水系统日优化调控方法,包括:建立水力学仿真模型和计算水头损失;对各泵站内机组流量进行优化分配;对梯级泵站的各泵站扬程进行优化分配;对一天内的各时段调水流量进行优化分配。本发明多约束、多目标的梯级泵站输水系统优化调控问题为研究对象,对调蓄池、渠道和梯级泵站进行综合考虑,合成水力学仿真模型和大系统分解-协调法,建立此类调水工程日优化调控模式,为实际调控提供科学依据和技术支撑。与现有技术相比,本发明包含调蓄池蓄量约束、渠段水头损失和梯级泵站总扬程等约束,使工程调控约束条件更符合此类调水工程输水系统,并以弃水和梯级泵站总电费之和最小为目标函数,使输水系统经济效益最大化。
申请公布号 CN105574605A 申请公布日期 2016.05.11
申请号 CN201510897648.8 申请日期 2015.12.08
申请人 中国水利水电科学研究院;北京市南水北调团城湖管理处 发明人 郑和震;房彦梅;吴辉明;雷晓辉;田雨;王旭;廖卫红;张云辉
分类号 G06Q10/04(2012.01)I;G06Q50/06(2012.01)I 主分类号 G06Q10/04(2012.01)I
代理机构 北京国林贸知识产权代理有限公司 11001 代理人 李桂玲;王敬智
主权项 一种梯级泵站输水系统日优化调度方法,其特征在于,所述方法的步骤如下:建立水力学仿真模型和计算水头损失的步骤:构建一维非恒定流水力学仿真模型,采用Preissmann四点时空偏心格式对方程组进行离散,用双扫描法求解,得到各渠段的水头损失;对各泵站内机组流量进行优化分配的步骤:本步骤包含m个泵站的优化,单个泵站的机组流量优化分配方法如下:①阶段变量:阶段变量采用序列数字i=1,2,……,n来表示;②状态变量:选取第i阶段至最末阶段n的累计流量作为状态变量:<maths num="0001"><math><![CDATA[<mrow><msub><mi>S</mi><mi>i</mi></msub><mo>=</mo><msubsup><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mi>i</mi></mrow><mi>n</mi></msubsup><msub><mi>q</mi><mi>j</mi></msub></mrow>]]></math><img file="FDA0000870797680000011.GIF" wi="254" he="94" /></maths>式中,j为1~n中的一个泵站机组;③决策变量:采用每台机组的流量q<sub>i</sub>作为决策变量;④状态转移方程:表示单个泵站中第i+1阶段的状态变量S<sub>i+1</sub>与第i阶段的状态变量S<sub>i</sub>和决策变量q<sub>i</sub>之间的关系:S<sub>i+1</sub>=S<sub>i</sub>‑q<sub>i</sub>式中,S<sub>1</sub>=Q<sub>k</sub>,S<sub>n+1</sub>=0;Q<sub>k</sub>为调水量;⑤目标函数:对于单个泵站,只考虑机组的电费,其目标函数的表达式为:<maths num="0002"><math><![CDATA[<mrow><mi>min</mi><mi>F</mi><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mfrac><mrow><msub><mi>&rho;gq</mi><mi>i</mi></msub><msub><mi>H</mi><mi>i</mi></msub></mrow><msub><mi>&eta;</mi><mrow><mi>s</mi><mi>t</mi><mi>i</mi></mrow></msub></mfrac><mo>&times;</mo><mi>&Delta;</mi><mi>t</mi><mo>&times;</mo><mi>c</mi></mrow>]]></math><img file="FDA0000870797680000012.GIF" wi="558" he="134" /></maths>ρ为水的密度,q<sub>i</sub>为第i台机组的出水量,H<sub>i</sub>为第i台机组的扬程,η<sub>sti</sub>为第i台机组的效率,c为电价,minF为在Δt时间范围内的最小费用函数,Δt为一个时间段,g为重力加速度;⑥约束条件:总流量约束:<maths num="0003"><math><![CDATA[<mrow><msub><mi>Q</mi><mi>k</mi></msub><mo>=</mo><msubsup><mi>&Sigma;</mi><mn>1</mn><mi>n</mi></msubsup><msub><mi>q</mi><mi>i</mi></msub></mrow>]]></math><img file="FDA0000870797680000021.GIF" wi="237" he="87" /></maths>机组过流能力约束:q<sub>imin</sub>≤q<sub>i</sub>≤q<sub>imax</sub>式中,Q<sub>k</sub>为调水总流量,q<sub>imin</sub>和q<sub>imax</sub>为第i台机组的最小和最大抽水流量;⑦求解方法:用逆序解法解此问题:F<sub>n+1</sub>(S<sub>n+1</sub>)=0<maths num="0004"><math><![CDATA[<mrow><msubsup><mi>F</mi><mi>i</mi><mo>*</mo></msubsup><mrow><mo>(</mo><msub><mi>S</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>=</mo><mi>m</mi><mi>i</mi><mi>n</mi><mo>{</mo><msub><mi>L</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>S</mi><mi>i</mi></msub><mo>,</mo><msub><mi>q</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>+</mo><msubsup><mi>F</mi><mrow><mi>i</mi><mi>+1</mi></mrow><mo>*</mo></msubsup><mrow><mo>(</mo><msub><mi>S</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mo>}</mo></mrow>]]></math><img file="FDA0000870797680000022.GIF" wi="750" he="87" /></maths>式中,L<sub>i</sub>(S<sub>i</sub>,q<sub>i</sub>)为某阶段的费用函数,<img file="FDA0000870797680000023.GIF" wi="174" he="71" />为某阶段的最小费用函数。根据上述递推方程式,按逆序逐阶段递推;对梯级泵站的各泵站扬程进行优化分配的步骤:对梯级泵站的各泵站扬程进行优化分配:①阶段变量:阶段变量采用序列数字j=1,2,……,m来表示,用泵站编号作为阶段变量,m为投入运行的泵站数;②状态变量:选取第j阶段至最末阶段m的累计扬程作为状态变量:<maths num="0005"><math><![CDATA[<mrow><msub><mi>S</mi><mi>j</mi></msub><mo>=</mo><msubsup><mi>&Sigma;</mi><mrow><mi>x</mi><mo>=</mo><mi>j</mi></mrow><mi>m</mi></msubsup><msub><mi>H</mi><mi>x</mi></msub></mrow>]]></math><img file="FDA0000870797680000024.GIF" wi="284" he="90" /></maths>式中,x为m~j中的变化量;③决策变量:采用每台机组的流量H<sub>j</sub>作为决策变量;④状态转移方程:表示梯级泵站中第j+1阶段的状态变量S<sub>j+1</sub>与第j阶段的状态变量S<sub>j</sub>和决策变量H<sub>j</sub>之间的关系:S<sub>j+1</sub>=S<sub>j</sub>‑H<sub>j</sub>式中,S<sub>1</sub>=H<sup>*</sup>,S<sub>m+1</sub>=0;H<sup>*</sup>表示梯级泵站总扬程;⑤目标函数:对于梯级泵站,考虑各泵站各机组的电费,其目标函数的表达式为:<maths num="0006"><math><![CDATA[<mrow><mi>min</mi><mi>F</mi><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mfrac><mrow><msub><mi>&rho;gq</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><msub><mi>H</mi><mi>j</mi></msub></mrow><mrow><msub><mi>&eta;</mi><mrow><mi>s</mi><mi>t</mi></mrow></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mrow></mfrac><mo>&times;</mo><mi>&Delta;</mi><mi>t</mi><mo>&times;</mo><mi>c</mi></mrow>]]></math><img file="FDA0000870797680000031.GIF" wi="645" he="151" /></maths>式中,q<sub>i,j</sub>是第j个泵站内第i个机组的出水量,η<sub>st</sub>(i,j)是第j个泵站内第i个机组的效率;⑥约束条件:单个泵站的扬程之和等于总扬程:<maths num="0007"><math><![CDATA[<mrow><msubsup><mi>&Sigma;</mi><mn>1</mn><mi>m</mi></msubsup><msub><mi>H</mi><mi>j</mi></msub><mo>=</mo><msup><mi>H</mi><mo>*</mo></msup><mo>=</mo><msub><mi>Z</mi><mi>m</mi></msub><mo>-</mo><msub><mi>Z</mi><mn>0</mn></msub><mo>+</mo><msubsup><mi>&Sigma;</mi><mrow><mn>1</mn><mo>,</mo><mn>2</mn></mrow><mrow><mi>m</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>m</mi></mrow></msubsup><msub><mi>h</mi><mrow><mi>j</mi><mo>,</mo><mi>j</mi><mo>+</mo><mn>1</mn></mrow></msub></mrow>]]></math><img file="FDA0000870797680000032.GIF" wi="741" he="102" /></maths>式中,Z<sub>m</sub>为最后一级泵站站后水位,Z<sub>0</sub>为第一级泵站站前水位,h<sub>j,j+1</sub>为第j级泵站和第j+1泵站之间渠段的水头损失;⑦求解方法:用逆序解法解此问题:F<sub>m+1</sub>(S<sub>m+1</sub>)=0<maths num="0008"><math><![CDATA[<mrow><msubsup><mi>F</mi><mi>j</mi><mo>*</mo></msubsup><mrow><mo>(</mo><msub><mi>S</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>=</mo><mi>m</mi><mi>i</mi><mi>n</mi><mo>{</mo><msub><mi>L</mi><mi>j</mi></msub><mrow><mo>(</mo><msub><mi>S</mi><mi>j</mi></msub><mo>,</mo><msub><mi>q</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>+</mo><msubsup><mi>F</mi><mrow><mi>j</mi><mo>+</mo><mn>1</mn></mrow><mo>*</mo></msubsup><mrow><mo>(</mo><msub><mi>S</mi><mrow><mi>j</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mo>}</mo></mrow>]]></math><img file="FDA0000870797680000033.GIF" wi="789" he="95" /></maths>式中,L<sub>j</sub>(S<sub>j</sub>,q<sub>j</sub>)为某阶段的费用函数,<img file="FDA0000870797680000034.GIF" wi="175" he="93" />为某阶段的最小费用函数;对一天内的各时段调水流量进行优化分配的步骤:针对不给定日调水总量或给定调水总量,对一天内的各时段调水流量进行优化分配。
地址 100038 北京市海淀区玉渊潭南路1号A座9层966