一种水电外送交易优化模型的构建方法
技术领域
[0001] 本
发明涉及电
力自动化领域,具体涉及一种水电外送交易优化模型的构建方法。
背景技术
[0002] 近年来,水力发电由于其清洁可再生等优势在世界范围内发展迅速,我国水电装机容量更是发展迅猛,已建设或者正在建设多个大规模水电站。但一般水资源富集区多集中在负荷水平不高的西部地区,往往需要外送到
电力负荷需求高的中东部地区进行消纳。
[0003] 水电加入到电力市场中进行交易,需要考虑水电站运行特性对交易的影响。同时,由于电量外送交易会对交易途经
电网产生影响,比如需要占用一定的输电通道,对途经地区自身区内的输电带来的影响,并带来新增的网络损耗。由于水利资源存在季节特性,在水电参与交易的情况下,在丰水期经常会存在大量的弃水现象,严重浪费水资源,而水电丰水期外送也容易造成输电通道阻塞等问题,故需要对水电外送交易的情况分时期进行分析。
发明内容
[0004] 为了使水电得到充分利用,并保证外送通道的建设运行成本,本发明提供一种水电外送交易优化模型的构建方法,充分考虑水利资源特性,综合分析了水利资源自身特性及季节性因素对水电站发电的影响,并考虑在典型运行时期下潮流约束对外送的影响,针对水电外送交易进行优化,促进水电跨区外送消纳。
[0005] 为了实现上述发明目的,本发明采取如下技术方案:
[0006] 本发明提供一种水电外送交易优化模型的构建方法,所述方法包括:
[0007] 计算水电站库容和发电量;
[0008] 建立基于电量的电力网络模型,进而得到线路剩余传输极限;
[0009] 建立水电外送交易优化模型;
[0010] 将建立的水电外送交易优化模型转化为混合整数规划进行求解。
[0011] 在计算水电站库容和发电量之前,获取水电站水利参数、水电外送交易信息和网架信息;
[0012] 所述水电站水利参数包括水电站的水位、库容、
水头和流量;
[0013] 所述水电外送交易信息包括购电
节点信息、水电站信息、购电方报价、购电量极限、售电量极限、网损电价信息和输电价格信息;
[0014] 所述网架信息包括拓扑结构、线路阻抗、发
电机节点
位置、负荷节点位置、
变压器线路的非标准变比。
[0015] 水电站库容的计算包括:
[0016] 设水电站m在t时刻的库容为Vm(t),有:
[0017] Vm(t)=Vm(t-1)+(Xm,in(t)-Xm,out(t))*Δt
[0018] 其中,Vm(t-1)表示水电站m在t-1时刻的库容,Xm,in(t)为水电站m在t时刻的入库流量,Xm,out(t)水电站m在t时刻的出库流量,Δt为时段时长。
[0019] 水电站发电量的计算包括:
[0020] 设 为水电站m在t时刻的发电功率,有:
[0021]
[0022] 其中, 表示水电站m在t时刻的流量; 为水电站m在t时刻的水电转换系数,其表示为 ηm表示水电站m的机组效率,Hm(t)表示水电站m在t
时段的水头。
[0023] 所述建立基于电量的电力网络模型,进而得到线路剩余传输极限之前包括:
[0024] 根据水利资源的季节特性,将水力资源发电情况分为丰水期和枯水期。
[0025] 所述建立基于电量的电力网络模型包括:
[0026] 根据电力网络拓扑结构、电力网络参数建立电力网络模型;
[0027] 所述电力网络参数包括线路阻抗、变压器线路的非标准变比、发电机节点、负荷节点位置信息以及线路传输极限。
[0028] 基于建立的基于电量的电力网络模型得到线路剩余传输极限包括:
[0029] 计算线路单位小时检修缺损电量和线路保留的安全备用裕量,再从线路传输极限中扣除线路单位小时检修缺损电量和线路保留的安全备用裕量,得到线路剩余传输极限P′max;有:
[0030] P′max=Pmax-Rbranch-Ubranch
[0031] 其中,Pmax为线路传输极限,Ubranch为线路单位小时检修缺损电量,Rbranch为线路保留的安全备用裕量。
[0032] 所述线路单位小时检修缺损电量的计算方法为线路检修时长除以全年小时数,再乘以线路传输极限;有:
[0033]
[0034] 其中,Ubranch为线路单位小时检修缺损电量;Tbranch为线路检修时长,单位为小时;Pmax为线路传输极限。
[0035] 所述线路保留的安全备用裕量表示为:
[0036] Rbranch=δbranch*Pmax
[0037] 其中,Rbranch为线路保留的安全备用裕量,δbranch为线路的安全备用系数,Pmax为线路传输极限。
[0038] 所述基于建立的基于电量的电力网络模型得到线路剩余传输极限之后包括,确定断面剩余传输极限,有:
[0039] P′ITR(r),max=PITR(r),max-∑Rbranch-∑Ubranch
[0040] 其中,P′ITR(r),max为断面r的剩余传输极限,PITR(r),max为断面r的传输极限,∑Rbranch为断面r上所有线路保留的安全备用裕量之和,∑Ubranch为断面r上所有线路单位小时检修缺损电量之和,Rbranch为线路保留的安全备用裕量,Ubranch为线路单位小时检修缺损电量。
[0041] 所述建立水电外送交易优化模型包括:
[0042] 确定水电外送交易优化模型的目标函数和约束条件。
[0043] 所述确定水电外送交易优化模型的目标函数包括:
[0044] 以水电外送交易的综合收益最大为优化目标,所述水电外送交易的综合收益指购电方购电
费用减去水电站成本以及
附加费用;所述附加费用包括电量传输过程中利用线路的输电费用和在线路上产生的网络损耗费用;
[0045] 用F表示水电外送交易的综合收益,有:
[0046] MAX F=∑bn(Qbuy,n)*Qbuy,n-∑sm(Qsell,m)*Qsell,m-Floss-Ftran[0047] 其中,∑bn(Qbuy,n)*Qbuy,n为购电方购电费用,bn(Qbuy,n)为购电节点n的关于电量的购买价格函数,Qbuy,n为购电节点n的购电量,且Qbuy,n=Pbuy,n*T,其中T为水电外送交易的时间跨度,Pbuy,n为购电节点n单位小时的购电量;
[0048] ∑sm(Qsell,m)*Qsell,m为水电站成本,sm(Qsell,m)为水电站m的关于电量的成本函数,Qsell,m为水电站m的售电量,且Qsell,m=Psell,m*T,Psell,m为水电站m单位小时的售电量;
[0049] Floss为电量传输过程中在线路上产生的网络损耗费用,Ftran为电量传输过程中利用线路的输电费用。
[0050] 电量传输过程中在线路上产生的网络损耗费用Floss表示为:
[0051] Floss=p*∑αi-j*T*ΔP
[0052] 其中,T为水电外送交易的时间跨度,ΔP为单位小时内线路流过的水电外送交易电量,T*ΔP为线路流过的水电外送交易电量,α为线路的线损率,p为单位电量的网络损耗费用。
[0053] 电量传输过程中利用线路的输电费用Ftran表示为:
[0054] Ftran=∑β*T*ΔP
[0055] 其中,T为水电外送交易的时间跨度,ΔP为单位小时内线路流过的水电外送交易电量,T*ΔP为线路流过的水电外送交易电量,β为线路流过单位电量所需缴纳的费用。
[0056] 水电外送交易优化模型的约束条件包括水电站库容约束、水电站流量约束、节点电量平衡约束、线路电量剩余传输极限约束、断面电量剩余传输极限约束、购电量约束和售电量约束。
[0057] 所述水电外送交易优化模型的约束条件包括:
[0058] 所述水电站库容约束表示为:
[0059]
[0060] 其中, 为水电站m在y时段时的平均库容, 为水电站m最小库容,为水电站m的最大库容;
[0061] 所述水电站流量约束表示为:
[0062]
[0063] 其中, 为水电站m在y时段时的平均发
电流量, 为水电站m最小发电流量, 为水电站m的最大发电流量;
[0064] 所述节点电量平衡约束包括水电站电量平衡约束、购电节点电量平衡约束和转运节点电量平衡约束,分别表示为:
[0065]
[0066]
[0067]
[0068] 其中,Psell,m为水电站m单位小时的售电量,i∈m表示转运节点i与水电站m相连,i≠m表示转运节点i与水电站m不重合,ΔPm-i表示单位小时内水电站m流向转运节点i的电量;
[0069] Pbuy,n表示购电节点n单位小时的购电量,i∈n表示转运节点i与购电节点n相连,i≠n表示转运节点i与购电节点n不重合,ΔPn-i表示单位小时内购电节点n流向转运节点i的电量;
[0070] ΔPk-j表示单位小时内转运节点k流向转运节点i的电量,i∈k表示转运节点i与转运节点k相连,i≠k表示转运节点i与转运节点k不重合;
[0071] 所述线路剩余传输极限约束表示为:
[0072] |P0+ΔP|≤P′max
[0073] 其中,P0为线路的初始潮流,ΔP为单位小时内线路流过的水电外送交易电量,P′max为线路剩余传输极限;
[0074] 所述断面剩余传输极限约束表示为:
[0075] ∑|P0+ΔP|≤P′ITR(r).max
[0076] 其中,P0为线路的初始潮流,ΔP为单位小时内线路流过的水电外送交易电量,P′ITR(r),max为断面r的剩余传输极限;
[0077] 所述购电量约束表示为:
[0078] 0≤Pbuy,n≤Pbuy,n,max
[0079] 其中,Pbuy,n为购电节点n单位小时的购电量,Pbuy,n,max为购电节点n的最大购电量;
[0080] 所述售电量约束表示为:
[0081] 0≤Psell,m≤Psell,m,max
[0082] 其中,Psell,m为水电站m单位小时的售电量,Psell,m,max为水电站m的最大售电量。
[0083] 所述将建立的水电外送交易优化模型转化为混合整数规划进行求解包括:
[0084] 采用混合整数规划通过计算机求解水电外送交易优化模型,求解结果包括水电站流量、售电方和购电方的成交电量、成交价格、成交单价、水电站收益、网络损耗费用和输电费用。
[0085] 与最接近的
现有技术相比,本发明提供的技术方案具有以下有益效果:
[0086] 1.本发明充分体现出以水电站作为外送交易的卖方时,水电站运行特性与发电量之间的关系,建立了考虑直流潮流约束和传输成本的水电外送交易优化模型,利于最大限度的利用水电资源在跨区范围进行消纳,减少弃水,提高水电站的收益;
[0087] 2.本发明具体分析了水利资源的季节特性,将年电量划分为丰水期和枯水期两个时间段,符合实际情况,更能体现出水电资源夏季水电大发对通道的阻塞情况。同时,以两个时间段组成的年电量形式完成交易,利于双边交易中电量合同的签订。另外,将丰水期和枯水期进行统一优化,规定全年总发电量,实现丰水期和枯水期的电量耦合关系;
[0088] 3.本发明在考虑传统的电量交易约束外,还考虑了直流潮流约束对外送交易的影响,相比于网络流的方法,更加符合电力系统实际运行情况;从电量长期交易的
角度,在对丰水期和枯水期进行统一优化时,分别考虑两个时期直流潮流约束的影响;
[0089] 4.本发明明确提出了计算网络损耗费用方法,并根据各线路的
电阻值计算折算得到各线路线损率,而不是全网都采用同一个网损率,有效地反应了各线路的特征,有利于在交易过程中更精确体现网络损耗对交易的影响;同时,在已有的交易模型中常常没有考虑输电费用,本模型同时将网络损耗费用和输电费用纳入优化目标中,以达到更加合理优化交易的目的;
[0090] 5.本发明从实际情况考虑,在交易的传输网络模型中考虑了线路检修计划和安全裕度备用,基于此提出线路剩余传输极限的定义,更贴近实际电网运行,使交易通道利用更加合理。
附图说明
[0091] 图1是本发明
实施例中水电外送交易优化模型的构建方法
流程图;
[0092] 图2是本发明实施例中IEEE 30算例系统拓扑结构图。
具体实施方式
[0093] 下面结合附图对本发明作进一步详细说明。
[0094] 本发明提供一种水电外送交易优化模型的构建方法,所述方法包括:
[0095] 计算水电站库容和发电量;
[0096] 建立基于电量的电力网络模型,进而得到线路剩余传输极限;
[0097] 建立水电外送交易优化模型;
[0098] 将建立的水电外送交易优化模型转化为混合整数规划进行求解。
[0099] 在计算水电站库容和发电量之前,获取水电站水利参数、水电外送交易信息和网架信息;
[0100] 所述水电站水利参数包括水电站的水位、库容、水头和流量;
[0101] 所述水电外送交易信息包括购电节点信息、水电站信息、购电方报价、购电量极限、售电量极限、网损电价信息和输电价格信息;
[0102] 所述网架信息包括拓扑结构、线路阻抗、发电机节点位置、负荷节点位置、变压器线路的非标准变比。
[0103] 水电站库容的计算包括:
[0104] 设水电站m在t时刻的库容为Vm(t),有:
[0105] Vm(t)=Vm(t-1)+(Xm,in(t)-Xm,out(t))*Δt
[0106] 其中,Vm(t-1)表示水电站m在t-1时刻的库容,Xm,in(t)为水电站m在t时刻的入库流量,Xm,out(t)水电站m在t时刻的出库流量,Δt为时段时长。
[0107] 水电站发电量的计算包括:
[0108] 设 为水电站m在t时刻的发电功率,有:
[0109]
[0110] 其中, 表示水电站m在t时刻的流量; 为水电站m在t时刻的水电转换系数,其表示为 ηm表示水电站m的机组效率,Hm(t)表示水电站m在t
时段的水头。
[0111] 所述建立基于电量的电力网络模型,进而得到线路剩余传输极限之前包括:
[0112] 根据水利资源的季节特性,将水力资源发电情况分为丰水期和枯水期。
[0113] 所述建立基于电量的电力网络模型包括:
[0114] 根据电力网络拓扑结构、电力网络参数建立电力网络模型;
[0115] 所述电力网络参数包括线路阻抗、变压器线路的非标准变比、发电机节点、负荷节点位置信息以及线路传输极限。
[0116] 基于建立的基于电量的电力网络模型得到线路剩余传输极限包括:
[0117] 计算线路单位小时检修缺损电量和线路保留的安全备用裕量,再从线路传输极限中扣除线路单位小时检修缺损电量和线路保留的安全备用裕量,得到线路剩余传输极限P′max;有:
[0118] P′max=Pmax-Rbranch-Ubranch
[0119] 其中,Pmax为线路传输极限,Ubranch为线路单位小时检修缺损电量,Rbranch为线路保留的安全备用裕量。
[0120] 所述线路单位小时检修缺损电量的计算方法为线路检修时长除以全年小时数,再乘以线路传输极限;有:
[0121]
[0122] 其中,Ubranch为线路单位小时检修缺损电量;Tbranch为线路检修时长,单位为小时;Pmax为线路传输极限。
[0123] 所述线路保留的安全备用裕量表示为:
[0124] Rbranch=δbranch*Pmax
[0125] 其中,Rbranch为线路保留的安全备用裕量,δbranch为线路的安全备用系数,Pmax为线路传输极限。
[0126] 所述基于建立的基于电量的电力网络模型得到线路剩余传输极限之后包括,确定断面剩余传输极限,有:
[0127] P′ITR(r),max=PITR(r),max-∑Rbranch-∑Ubranch
[0128] 其中,P′ITR(r),max为断面r的剩余传输极限,PITR(r),max为断面r的传输极限,∑Rbranch为断面r上所有线路保留的安全备用裕量之和,∑Ubranch为断面r上所有线路单位小时检修缺损电量之和,Rbranch为线路保留的安全备用裕量,Ubranch为线路单位小时检修缺损电量。
[0129] 所述建立水电外送交易优化模型包括:
[0130] 确定水电外送交易优化模型的目标函数和约束条件。
[0131] 所述确定水电外送交易优化模型的目标函数包括:
[0132] 以水电外送交易的综合收益最大为优化目标,所述水电外送交易的综合收益指购电方购电费用减去水电站成本以及附加费用;所述附加费用包括电量传输过程中利用线路的输电费用和在线路上产生的网络损耗费用;
[0133] 用F表示水电外送交易的综合收益,有:
[0134] MAX F=∑bn(Qbuy,n)*Qbuy,n-∑sm(Qsell,m)*Qsell,m-Floss-Gtran[0135] 其中,∑bn(Qbuy,n)*Qbuy,n为购电方购电费用,bn(Qbuy,n)为购电节点n的关于电量的购买价格函数,Qbuy,n为购电节点n的购电量,且Qbuy,n=Pbuy,n*T,其中T为水电外送交易的时间跨度,Pbuy,n为购电节点n单位小时的购电量;
[0136] ∑sm(Qsell,m)*Qsell,m为水电站成本,sm(Qsell,m)为水电站m的关于电量的成本函数,Qsell,m为水电站m的售电量,且Qsell,m=Psell,m*T,Psell,m为水电站m单位小时的售电量;
[0137] Floss为电量传输过程中在线路上产生的网络损耗费用,Ftran为电量传输过程中利用线路的输电费用。
[0138] 电量传输过程中在线路上产生的网络损耗费用Floss表示为:
[0139] Floss=p*∑αi-j*T*ΔP
[0140] 其中,T为水电外送交易的时间跨度,ΔP为单位小时内线路流过的水电外送交易电量,T*ΔP为线路流过的水电外送交易电量,α为线路的线损率,p为单位电量的网络损耗费用。
[0141] 电量传输过程中利用线路的输电费用Ftran表示为:
[0142] Ftran=∑β*T*ΔP
[0143] 其中,T为水电外送交易的时间跨度,ΔP为单位小时内线路流过的水电外送交易电量,T*ΔP为线路流过的水电外送交易电量,β为线路流过单位电量所需缴纳的费用。
[0144] 水电外送交易优化模型的约束条件包括水电站库容约束、水电站流量约束、节点电量平衡约束、线路电量剩余传输极限约束、断面电量剩余传输极限约束、购电量约束和售电量约束。
[0145] 所述水电外送交易优化模型的约束条件包括:
[0146] 所述水电站库容约束表示为:
[0147]
[0148] 其中, 为水电站m在y时段时的平均库容, 为水电站m最小库容,为水电站m的最大库容;
[0149] 所述水电站流量约束表示为:
[0150]
[0151] 其中, 为水电站m在y时段时的平均发电流量, 为水电站m最小发电流量, 为水电站m的最大发电流量;
[0152] 所述节点电量平衡约束包括水电站电量平衡约束、购电节点电量平衡约束和转运节点电量平衡约束,分别表示为:
[0153]
[0154]
[0155]
[0156] 其中,Psell,m为水电站m单位小时的售电量,i∈m表示转运节点i与水电站m相连,i≠m表示转运节点i与水电站m不重合,ΔPm-i表示单位小时内水电站m流向转运节点i的电量;
[0157] Pbuy,n表示购电节点n单位小时的购电量,i∈n表示转运节点i与购电节点n相连,i≠n表示转运节点i与购电节点n不重合,ΔPn-i表示单位小时内购电节点n流向转运节点i的电量;
[0158] ΔPk-j表示单位小时内转运节点k流向转运节点i的电量,i∈k表示转运节点i与转运节点k相连,i≠k表示转运节点i与转运节点k不重合;
[0159] 所述线路剩余传输极限约束表示为:
[0160] |P0+ΔP|≤P′max
[0161] 其中,P0为线路的初始潮流,ΔP为单位小时内线路流过的水电外送交易电量,P′max为线路剩余传输极限;
[0162] 所述断面剩余传输极限约束表示为:
[0163] ∑|P0+ΔP|≤P′ITR(r).max
[0164] 其中,P0为线路的初始潮流,ΔP为单位小时内线路流过的水电外送交易电量,P′ITR(r),max为断面r的剩余传输极限;
[0165] 所述购电量约束表示为:
[0166] 0≤Pbuy,n≤Pbuy,n,max
[0167] 其中,Pbuy,n为购电节点n单位小时的购电量,Pbuy,n,max为购电节点n的最大购电量;
[0168] 所述售电量约束表示为:
[0169] 0≤Psell,m≤Psell,m,max
[0170] 其中,Psell,m为水电站m单位小时的售电量,Psell,m,max为水电站m的最大售电量。
[0171] 所述将建立的水电外送交易优化模型转化为混合整数规划进行求解包括:
[0172] 采用混合整数规划通过计算机求解水电外送交易优化模型,求解结果包括水电站流量、售电方和购电方的成交电量、成交价格、成交单价、水电站收益、网络损耗费用和输电费用。
[0173] 本例采用如图2所示的IEEE 30系统进行说明。水电站信息见表1、网络线路数据、节点数据分别见表2和表3,数据以标幺值的方式表示,其中基准容量SB=100MVA,基准
电压UB=Uav,基准交易时间长度设为一年,其中包括丰水期和枯水期,丰水期:6月~10月共5个月;其余7个月为枯水期。
[0174] 表1
[0175]
[0176] 表1给出了水电站的节点编号信息,以及安全运行约束等信息。另外,水电站的库容变化不能超出允许的范围。
[0177] 表2
[0178]
[0179]
[0180] 表3
[0181]
[0182]
[0183] 水电站1和2的全年最大售电量总共为262800MWh,购电方节点26、29的最大买进电量皆为262800MWh,假设丰水期外送发电量不低于全年发电量的70%,水电站1,水电站2的报价函数分别为: f2(S2)=1.2*S2+240(元/MWh);购电节点26,购电节点29的报价函数分别为:g26(B26)=-1.3*B26+480(元/MWh),g29(B29)=-1.2*B29+
470(元/MWh)。每条线路输电费用设为10元/MWh;单位电量的网损费用设为400元/MWh,全网网损率取1.0%;所有线路的检修计划时长均为64h;系统安全备用裕量系数5%。
[0184] 经过数据准备,建立水电外送交易优化模型后,求解得到如下丰水期交易结果和枯水期交易结果,分别如表4和表5:
[0185] 表4
[0186]
[0187] 表5
[0188]
[0189] 注:购售电费用差为购电方购电费用减去水电站成本;交易净收益为购电方购电费用减去水电站成本后再减去附加费用。
[0190] 最后应当说明的是:以上实施例仅用以说明本发明的技术方案而非对其限制,所属领域的普通技术人员参照上述实施例依然可以对本发明的具体实施方式进行
修改或者等同替换,这些未脱离本发明精神和范围的任何修改或者等同替换,均在
申请待批的本发明的
权利要求保护范围之内。