CN106875026B - 一种电力市场环境下中长期输电网拓展规划的混合性规划方法 - Google Patents

一种电力市场环境下中长期输电网拓展规划的混合性规划方法 Download PDF

Info

Publication number
CN106875026B
CN106875026B CN201510925496.8A CN201510925496A CN106875026B CN 106875026 B CN106875026 B CN 106875026B CN 201510925496 A CN201510925496 A CN 201510925496A CN 106875026 B CN106875026 B CN 106875026B
Authority
CN
China
Prior art keywords
line
channels
node
power
load
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201510925496.8A
Other languages
English (en)
Other versions
CN106875026A (zh
Inventor
王一哲
汤涌
董朝阳
顾卓远
韩奕
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
State Grid Gansu Electric Power Co Ltd
Economic and Technological Research Institute of State Grid Gansu Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
State Grid Gansu Electric Power Co Ltd
Economic and Technological Research Institute of State Grid Gansu Electric Power Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by State Grid Corp of China SGCC, China Electric Power Research Institute Co Ltd CEPRI, State Grid Gansu Electric Power Co Ltd, Economic and Technological Research Institute of State Grid Gansu Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201510925496.8A priority Critical patent/CN106875026B/zh
Publication of CN106875026A publication Critical patent/CN106875026A/zh
Application granted granted Critical
Publication of CN106875026B publication Critical patent/CN106875026B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Energy or water supply
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00—Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Landscapes

  • Business, Economics & Management (AREA)
  • Engineering & Computer Science (AREA)
  • Economics (AREA)
  • Human Resources & Organizations (AREA)
  • Strategic Management (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Marketing (AREA)
  • General Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Tourism & Hospitality (AREA)
  • Public Health (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Water Supply & Treatment (AREA)
  • Development Economics (AREA)
  • Game Theory and Decision Science (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明涉及一种电力市场环境下中长期输电网拓展规划的混合性规划方法，包括：步骤1：设定混合性规划模型的目标函数；步骤2：设定电网运行基本约束条件；步骤3：利用N‑1强制约束条件和N‑2柔性约束条件检测电网；步骤4：判断是否需要特定可靠性约束条件加强输电网以满足N‑1强制约束和N‑2柔性约束条件；步骤5：寻求最优规划方案进行仿真计算。本发明提供的技术方案可以为输电网中长期拓展规划提供评估新的评价指标、新的目标函数、更加合理全面的可靠性和经济性约束条件，该方法具有简单、实用、全面、操作性强的特点。

Description

一种电力市场环境下中长期输电网拓展规划的混合性规划 方法
技术领域
本发明涉及一种电力系统规划方法，具体涉及一种电力市场环境下中长期输电网拓展规划的混合性规划方法。
背景技术
随着我国电力改革的进一步深入和电力市场的快速发展，发电企业和电网公司都面临着新的挑战和机遇。电力市场的兴起也为输电网拓展规划带来了新特点。建立电力市场的主要目的：通过在电力工业中引入新的参与者、竞争者来打破电力系统的垄断，实现优化电力工业资源组合配置，提高电力系统的运行效率的目的。输电网的拓展规划不仅需要计及电网本身的不确定性因素，还需要在决策过程中充分考虑发电厂商、输电网公司、用户的利益需求，以及市场监管集团的管理职能等。因此，输电网拓展规划成为涉及金融、投资、电力工程技术等多个领域的交叉课题，变得更加复杂。
与此同时，国内外电网规划准则和方法在应对新一轮变化中已显不足。在电力工业集中管制的环境下，确定性规划准则是应用最为广泛的一种规划标准，主要的形式为将N-1、N-2作为可靠性约束条件。其优点主要有思路清晰，易于使用，对计算工具要求低；缺点有：1.只分析了故障的影响不考虑故障发生的概率，可能造成过度投资；2.不能很好地处理电网规划中的不确定因素。随着电力工业结构的重组、电力市场的不断发展与完善，概率性规划准则由于其在处理不确定因素的良好表现，也得到了迅速的发展和应用但是，概率性规划准则由于其采用了大量的统计数据，而统计数据本身的随机性会影响计算的准确性。
发明内容
为解决上述现有技术中的不足，本发明的目的是提供一种电力市场环境下中长期输电网拓展规划的混合性规划方法。
本发明的目的是采用下述技术方案实现的：
本发明提供一种电力市场环境下中长期输电网拓展规划的混合性规划方法，其改进之处在于，所述方法包括下述步骤：
步骤1：设定混合性规划模型的目标函数；(混合性规划模型是确定性规划模型和概率性规划模型的混合)；
步骤2：设定电网运行基本约束条件；
步骤3：利用N-1强制约束条件和N-2柔性约束条件检测电网；
步骤4：判断是否需要特定可靠性约束条件加强输电网以满足N-1强制约束和N-2柔性约束；
步骤5：寻求最优规划方案进行仿真计算。
进一步地，所述步骤1中，设定混合性规划模型的目标函数满足以下原则：
1)采用资本回收系数将整个设备(整个设备是指需要建设的设备的总称)的投资费用由一个时间点归算到设备寿命的每一年；
2)混合性规划模型的目标函数计算步长为一年，若有投资费用，则投资费用发生在年初，社会总效益和电能损失值在年终计算；
3)采用净现值计算规划周期内的社会总效益和投资费用；一个规划周期中社会总效益和投资费用均以规划周期开始时间为计算净现值的时间点；社会效益盈余等于社会总效益减去投资费用；
4)混合性规划模型通过不同的场景来模拟电力系统不同的负荷水平，同时赋予不同场景相应的比重系数用于衡量对应负荷在一年中的持续时间和负荷情况；混合性规划模型设置多个场景来模拟电力系统中各个典型的负荷水平及其持续时间，具体设置情况视具体的电网情况而定；
混合性规划模型的目标函数如下：
Figure BDA0000877580290000021
式中：P为特定故障发生概率；r为折现率；ωa为场景a的比重系数；Ωa为所有场景的集合；ΩD为所有负荷点的集合；Ωd为每个负荷所有区域的集合；a表示场景；Q表示发电机；d为每个负荷节点的负荷；h和b分别表示在场景a下负荷d的区域以及在场景a下发电机Q的区域；
Figure BDA0000877580290000022
为在场景a下负荷d在其h区域中的购电价格；
Figure BDA0000877580290000024
为在场景a下负荷d在其h区域中的负荷值；ΩG为所有发电机组的集合；ΩQ为发电机组所有区域的集合；
Figure BDA0000877580290000023
为在场景a下发电机Q在其b区域中的供电价格；
Figure BDA0000877580290000025
为在场景a下发电机Q在其b区域中的出力；
Figure BDA0000877580290000026
为 在场景a下负荷d在其h区域中损失单位电量对应的用户可靠性费用；
Figure BDA00008775802900000310
为在场景a下负荷d在其h区域中的损失电力；σ为资本回收系数；Cij为在通道i-j中建设一条线路的总费用；nij为在通道i-j中新建的线路数量；
为计算资本回收系数σ，设新建线路最少运行的寿命为t年，因此资本回收周期为t年，寿命周期t年的折现率为r,表达式如下式(2)：
Figure BDA0000877580290000031
Figure BDA00008775802900000311
为在场景a下负荷d在其h区域中的损失电力在一年中损失的电量时，采用期望缺供电量EENS进行计算，单位为MWh/year，其数值通过负荷持续曲线计算得到，如下式(3)：
Figure BDA0000877580290000032
式中：LPk-l为到达负荷节点k的最大潮流，单位为MW；LPk为节点k最大的用电负荷，单位为MW；φk为负荷节点k上有效负荷持续曲线。
进一步地，所述步骤2中，电网运行基本约束条件包括：
功率约束：
Figure BDA0000877580290000033
直流潮流计算：
Figure BDA0000877580290000034
线路潮流约束：
原有线路的潮流约束
Figure BDA0000877580290000035
新建线路的潮流约束
Figure BDA0000877580290000036
发电机机组出力约束：
Figure BDA0000877580290000037
节点上负荷损失值的约束：
0≤La≤d (9)
新建线路回数的约束：
Figure BDA0000877580290000038
其他：
Figure BDA0000877580290000039
式中：Ωa为所有场景的集合；S0为原网络分支节点的关联矩阵；
Figure BDA0000877580290000049
为在场景a下通道i-j中线路通过的潮流；f0a为在场景a下
Figure BDA00008775802900000415
的向量和；S为新建线路的网络分支节点的关联矩阵；fa为在场景a下通道i-j中新建线路通过的潮流，为场景a下
Figure BDA00008775802900000410
的向量和；
Figure BDA00008775802900000411
在场景a下通道i-j中新建每一条线路通过的潮流；ga为在场景a下节点上发电机的出力，即
Figure BDA00008775802900000412
的和，
Figure BDA00008775802900000420
为在场景a下发电机i在其b区域中的出力；La为在场景a节点上的电力损失，即
Figure BDA00008775802900000413
的和，
Figure BDA00008775802900000419
为在场景a下负荷d在其h区域中的损失电力；da为在场景a下节点的负荷值，即
Figure BDA00008775802900000414
的向量和，
Figure BDA00008775802900000418
为在场景a下负荷d在其h区域中的负荷值；γij为通道i-j中一条线路的导纳；
Figure BDA00008775802900000417
为通道i-j中原有的线路数量；
Figure BDA00008775802900000416
为在场景a下节点i的电压角度；
Figure BDA0000877580290000048
为原有网络线路的集合；
Figure BDA00008775802900000421
为通道i-j中一条线路的最大容量；
Figure BDA00008775802900000423
为发电机组的最大功率；d为每个负荷节点的负荷；
Figure BDA00008775802900000422
为通道i-j中线路的线路能够建设的最大条数；nij为非负整数；f0a和
Figure BDA00008775802900000424
为实数。
进一步地，所述步骤3中，N-1强制约束条件表述为：
S0mnf0mn+Smnfmn+gmn＝d (12)
Figure BDA0000877580290000041
Figure BDA0000877580290000042
Figure BDA0000877580290000043
Figure BDA0000877580290000044
Figure BDA0000877580290000045
Figure BDA0000877580290000046
Figure BDA0000877580290000047
(m,n)∈ψ (20)
Eq.(4)-(11) (21)
其中，S0mn为通道m-n中一条线路发生故障下原网络分支节点的关联矩阵；f0mn为通道m-n中一条线路发生故障下原有线路通过的潮流；f0mn为通道m-n中一条线路发生故障下
Figure BDA00008775802900000425
的向量和；Smn为新建线路的网络分支节点的关联矩阵；fmn为通道m-n中新建线路通过的潮流；gmn为通道m-n中一条线路发生故障下节点上发电机的出力；d为每个负荷节点的 负荷；
Figure BDA00008775802900000517
为在通道m-n中一条线路发生故障下剩余线路通过的潮流；γij为通道i-j中一条线路的导纳；
Figure BDA00008775802900000518
为通道i-j中原有的线路数量；
Figure BDA00008775802900000522
和
Figure BDA00008775802900000521
为在通道m-n中一条线路发生故障下节点i和节点j的电压角度；
Figure BDA00008775802900000519
为通道i-j中一条线路的最大容量；
Figure BDA00008775802900000520
在通道m-n中一条线路发生故障下通道i-j中新建每一条线路通过的潮流；
Figure BDA00008775802900000523
为发电机组的最大功率；
Figure BDA00008775802900000524
为通道i-j中线路的线路能够建设的最大条数；nij为非负整数；ψ为故障集合；
N-2柔性约束条件的表述为：
S0pqf0pq+Spqfpq+gpq+lpq＝d,ij＝pq≠st (22)
Figure BDA0000877580290000051
Figure BDA0000877580290000052
Figure BDA0000877580290000053
S0stf0st+Sstfst+gst+lst＝d,ij＝st≠pq (26)
Figure BDA0000877580290000054
Figure BDA0000877580290000055
Figure BDA0000877580290000056
S0pqstf0pqst+Spqstfpqst+gpqst+lpqst＝d,
(30)
ij＝pq＝st且
Figure BDA0000877580290000057
Figure BDA00008775802900000525
(31)
ij＝pq＝st且
Figure BDA0000877580290000058
Figure BDA0000877580290000059
且
Figure BDA00008775802900000510
Figure BDA00008775802900000511
且
Figure BDA00008775802900000512
Figure BDA00008775802900000513
Figure BDA00008775802900000514
Figure BDA00008775802900000515
Figure BDA00008775802900000516
0≤lpq≤d,0≤lst≤d,0≤lpqst≤d (38)
(p,q)∈ψ,(s,t)∈ψ (39)
Eq.(4)-(11) (40)
其中，S0pq为通道p-q中一条线路发生故障下原网络分支节点的关联矩阵；f0pq为通道p-q中一条线路发生故障下原有线路通过的潮流；f0pq为通道p-q中一条线路发生故障下
Figure BDA0000877580290000061
的向量和；Spq为新建线路的网络分支节点的关联矩阵；fpq为通道p-q中新建线路通过的潮流；gpq为通道p-q中一条线路发生故障下节点上发电机的出力；d为每个负荷节点的负荷；
Figure BDA0000877580290000062
为在通道p-q中一条线路发生故障下剩余线路通过的潮流；
Figure BDA0000877580290000064
和
Figure BDA0000877580290000063
为在通道p-q中一条线路发生故障下节点i和节点j的电压角度；S0st为通道s-t中一条线路发生故障下原网络分支节点的关联矩阵；fost为通道s-t中一条线路发生故障下原有线路通过的潮流；f0st为通道s-t中一条线路发生故障下
Figure BDA0000877580290000065
的向量和；Sst为新建线路的网络分支节点的关联矩阵；fst为通道s-t中新建线路通过的潮流；gst为通道s-t中一条线路发生故障下节点上发电机的出力；d为每个负荷节点的负荷；
Figure BDA0000877580290000066
为在通道s-t中一条线路发生故障下剩余线路通过的潮流；
Figure BDA0000877580290000067
和
Figure BDA0000877580290000068
为在通道s-t中一条线路发生故障下节点i和节点j的电压角度；S0pqst为通道i-j中两条线路发生故障下原网络分支节点的关联矩阵；f0pqst为通道i-j中两条线路发生故障下原有线路通过的潮流；f0pqst为通道i-j中两条线路发生故障下
Figure BDA0000877580290000069
的向量和；Spqst为新建线路的网络分支节点的关联矩阵；fpqst为通道i-j中新建线路通过的潮流；gpqst为通道i-j中两条线路发生故障下节点上发电机的出力；d为每个负荷节点的负荷；
Figure BDA00008775802900000610
通道i-j中两条线路发生故障下剩余线路通过的潮流；
Figure BDA00008775802900000612
和
Figure BDA00008775802900000611
为在通道i-j中两条线路发生故障下节点i和节点j的电压角度；γij为通道i-j中一条线路的导纳；
Figure BDA00008775802900000613
为通道i-j中原有的线路数量；
Figure BDA00008775802900000618
为通道i-j中一条线路的最大容量；
Figure BDA00008775802900000614
在通道m-n中一条线路发生故障下通道i-j中新建每一条线路通过的潮流；
Figure BDA00008775802900000619
为发电机组的最大功率；
Figure BDA00008775802900000615
为通道i-j中线路的线路能够建设的最大条数；上标pq,st表示在通道p-q和s-t中一回线路发生了故障，pqst表示在同一通道中两回线路同时发生故障，ψ为故障集合；lpq为通道p-q中一回线路发生故障下，电力系统中每个节点损失的电力；
Figure BDA00008775802900000616
为通道i-j中一回线路的最大容量；nij为通道i-j中新建的线路回数；
Figure BDA00008775802900000620
为通道i-j中新建线路中的功率；lst为通道s-t中一回线路发生故障下，电力系统中每个节点损失的电力；
Figure BDA00008775802900000617
为通道i-j中一回线路的最大容量；nij为通道i-j中新建的线路回数；
Figure BDA00008775802900000621
为通道i-j中新建线路中的功率；lpqst为通道i-j中两条线路发生故障下，电力系统中每个节点损失的电力；gpqst为通道i-j中两条线路发生故障下，节点上发电机的出力；
Figure BDA00008775802900000622
为通道i-j两回线路故障时，新建线路中的功率；
Figure BDA0000877580290000071
为通道i-j中一回线路的最大容量；nij为通道i-j中新建的线路回数；
Figure BDA0000877580290000072
为通道i-j中线路没有发生故障时通过的潮流；θi为通道i-j中线路没有发生故障时节点i的电压角度；θj为通道i-j中线路没有发生故障时节点j的电压角度；
Figure BDA0000877580290000074
为通道i-j中线路没有发生故障时通过的潮流；
Figure BDA0000877580290000076
为通道i-j中一条线路的最大容量；
Figure BDA0000877580290000075
为通道i-j中原有线路的回数；
Figure BDA0000877580290000073
为通道i-j中一条线路的最大容量；nij为通道i-j中新建的线路回数；fij为通道i-j中新建线路中的功率；
模型通过N-2柔性约束条件来计算在N-2故障情况下电力系统的负荷损失。其允许有损失负荷，但是允许损失负荷的大小由目标函数和其他约束条件确定，实现规划方案在N-2柔性约束条件下取得了经济性和可靠性的平衡；模型使用N-1强制性约束条件校验电网，电网在N-1故障情况下不允许出现负荷损失，保证得到的电网规划方案必然满足N-1准则。
进一步地，所述步骤4中，特定可靠性约束条件的表述如下式：
Lk＜T (41)
缺电时间期望的数值通过负荷持续曲线计算得到：
Figure BDA0000877580290000077
式中：Lk为节点k的缺电时间期望，单位为hours/year；T为缺电时间期望(LOLE)在特定可靠性约束中的阈值，单位为hours/year；LPk-l为到达节点k的最大潮流，单位为MW；LPk为节点k最大的用电负荷，单位为MW；φk为节点k上有效负荷持续曲线。
进一步地，所述步骤5中，混合性规划模型已对所有的约束条件进行了线性化处理，使用直流模型和最优潮流模型的方法，将输电网拓展规划问题转为混合整数线性规划问题，通过Matlab和YALMIP调用的商业软件CPLEX实现仿真计算。
本发明提供的技术方案具有的优异效果是：
本发明提供一种基于电力市场环境下应用于输电网拓展规划中的混合性规划方法，属于电力系统规划领域。基于确定性规划准则和概率性规划准则，提出了一种混合性规划模型，综合考虑N-1和N-2约束、经济性约束、可靠性约束、故障发生概率、可靠性成本、投资费用、期望缺供电量(EENS)以及缺电时间期望(LOLE)等，同时通过定义多个场景来模拟电力系统的负荷水平，并考虑电力市场的因素，将输电网拓展规划问题转为混合整数线性规划问题。所述方法可以为输电网中长期拓展规划提供评估新的评价指标、新的目标函数、更 加合理全面的可靠性和经济性约束条件，该方法具有简单、实用、全面、操作性强的特点。
附图说明
图1是本发明提供的电力市场环境下应用于输电网中长期拓展规划的混合性规划方法流程图；
图2是本发明提供的具体实施例的江西省南部500/220kV电网示意图；
图3是本发明提供的两种规划模型对应的结果示意图。
具体实施方式
下面结合附图对本发明的具体实施方式作进一步的详细说明。
以下描述和附图充分地示出本发明的具体实施方案，以使本领域的技术人员能够实践它们。其他实施方案可以包括结构的、逻辑的、电气的、过程的以及其他的改变。实施例仅代表可能的变化。除非明确要求，否则单独的组件和功能是可选的，并且操作的顺序可以变化。一些实施方案的部分和特征可以被包括在或替换其他实施方案的部分和特征。本发明的实施方案的范围包括权利要求书的整个范围，以及权利要求书的所有可获得的等同物。在本文中，本发明的这些实施方案可以被单独地或总地用术语“发明”来表示，这仅仅是为了方便，并且如果事实上公开了超过一个的发明，不是要自动地限制该应用的范围为任何单个发明或发明构思。
本发明提供一种基于电力市场环境下应用于输电网拓展规划中的混合性规划方法，流程图如图1所示，包括下述步骤：
步骤1：设定合理的目标函数，真实反映电网建设的要求：
在中长期输电网规划中，投资费用发生在某一个时间点，而社会效益和电量损失则发生在一个时间段上，两者之间不能直接比较。为了解决以上问题，本文采用以下解决方案：
1)采用资本回收系数(Capital Recovery Factor)将整个设备的投资费用由一个时间点归算到设备寿命的每一年。
2)混合性规划模型的目标函数计算步长为一年，且认为：如有投资费用，则投资费用发生在年初，社会总效益和电能损失值都在年终计算。
3)为了便于与传统的规划模型进行比较，本文中采用净现值的概念用于计算规划周期内的社会总效益和投资费用。一个规划周期中社会总效益和投资费用都以规划周期开始时间为计算净现值的时间点。社会效益盈余等于社会总效益减去投资费用。
4)模型通过不同的场景来模拟系统不同的负荷水平，同时赋予不同场景相应的比重系数
用于衡量对应负荷的在一年(8760小时)中的持续时间和负荷情况。
目标函数表达式如下：
Figure BDA0000877580290000091
式中：P为特定故障发生概率；r为折现率；ωa为场景a的比重系数；Ωa为所有场景的集合；ΩD为所有负荷点的集合；Ωd为每个负荷所有区域的集合；a表示场景；Q表示发电机；d为每个负荷节点的负荷；h和b分别表示在场景a下负荷d的区域以及在场景a下发电机Q的区域；
Figure BDA0000877580290000095
为在场景a下负荷d在其h区域中的购电价格；
Figure BDA0000877580290000097
为在场景a下负荷d在其h区域中的负荷值；ΩG为所有发电机组的集合；ΩQ为发电机组所有区域的集合；
Figure BDA0000877580290000096
为在场景a下发电机Q在其b区域中的供电价格；
Figure BDA0000877580290000098
为在场景a下发电机Q在其b区域中的出力；
Figure BDA00008775802900000910
为在场景a下负荷d在其h区域中损失单位电量对应的用户可靠性费用；
Figure BDA0000877580290000099
为在场景a下负荷d在其h区域中的损失电力；σ为资本回收系数；Cij为在通道i-j中建设一条线路的总费用；nij为在通道i-j中新建的线路数量；其中：
为计算资本回收系数σ，假设如下：新建线路最少运行的寿命为t年，因此资本回收周期为t年，该寿命周期中的折现率为r,表达式如下(2)：
Figure BDA0000877580290000092
Ek(EENS,MWh/year)为期望缺供电量，其数值可通过负荷持续曲线计算得到，见式(3)：
Figure BDA0000877580290000093
式中：LPk-l(MW)为到达节点k的最大潮流；LPk(MW)为节点k最大的用电负荷；φk为节点k上有效负荷持续曲线。
步骤2：设定电网运行基本约束条件
功率约束：
Figure BDA0000877580290000094
直流潮流计算：
Figure BDA0000877580290000101
线路潮流约束：
原有线路的潮流约束
Figure BDA0000877580290000102
新建线路的潮流约束
Figure BDA0000877580290000103
发电机机组出力约束：
Figure BDA0000877580290000104
节点上负荷损失值的约束：
0≤la≤d (9)
新建线路回数的约束：
Figure BDA0000877580290000105
其他：
Figure BDA0000877580290000106
式中：S0为原网络分支节点的关联矩阵；
Figure BDA00008775802900001013
为在场景a下通道i-j中线路通过的潮流；f0a为在场景a下
Figure BDA00008775802900001014
的向量和；S为新建线路的网络分支节点的关联矩阵；fa为在场景a下通道i-j中新建线路通过的潮流；ga为在场景a下节点上发电机的出力，即
Figure BDA00008775802900001020
的和；La为在场景a节点上的电力损失，即
Figure BDA00008775802900001015
的和；da为在场景a下节点的负荷值，即
Figure BDA00008775802900001021
的向量和；γij为通道i-j中一条线路的导纳；
Figure BDA00008775802900001016
为通道i-j中原有的线路数量；
Figure BDA00008775802900001022
为在场景a下节点i的电压角度；
Figure BDA00008775802900001012
为原有网络线路的集合；
Figure BDA00008775802900001017
为通道i-j中一条线路的最大容量；
Figure BDA00008775802900001023
为发电机组的最大功率；d为每个负荷节点的负荷；
Figure BDA00008775802900001018
为通道i-j中线路的线路能够建设的最大条数；nij为非负整数；f0a和
Figure BDA00008775802900001019
为实数。
步骤3：电网的N-1强制约束条件和N-2柔性约束条件的表述
N-1强制性约束条件的表述为：
S0mnf0mn+Smnfmn+gmn＝d (12)
Figure BDA0000877580290000107
Figure BDA0000877580290000108
Figure BDA0000877580290000109
Figure BDA00008775802900001010
Figure BDA00008775802900001011
Figure BDA0000877580290000111
Figure BDA0000877580290000112
(m,n)∈ψ (20)
Eq.(4)-(11) (21)
其中，S0mn为通道m-n中一条线路发生故障下原网络分支节点的关联矩阵；f0mn为通道m-n中一条线路发生故障下原有线路通过的潮流；f0mn为通道m-n中一条线路发生故障下
Figure BDA00008775802900001113
的向量和；Smn为新建线路的网络分支节点的关联矩阵；fmn为通道m-n中新建线路通过的潮流；gmn为通道m-n中一条线路发生故障下节点上发电机的出力；d为每个负荷节点的负荷；
Figure BDA00008775802900001114
为在通道m-n中一条线路发生故障下剩余线路通过的潮流；γij为通道i-j中一条线路的导纳；
Figure BDA00008775802900001115
为通道i-j中原有的线路数量；
Figure BDA00008775802900001121
和
Figure BDA00008775802900001120
为在通道m-n中一条线路发生故障下节点i和节点j的电压角度；
Figure BDA00008775802900001116
为通道i-j中一条线路的最大容量；
Figure BDA00008775802900001117
在通道m-n中一条线路发生故障下通道i-j中新建每一条线路通过的潮流；
Figure BDA00008775802900001119
为发电机组的最大功率；
Figure BDA00008775802900001118
为通道i-j中线路的线路能够建设的最大条数；nij为非负整数；ψ为故障集合；
N-2柔性约束条件的表述为：
S0pqf0pq+Spqfpq+gpq+lpq＝d,ij＝pq≠st (22)
Figure BDA0000877580290000113
Figure BDA0000877580290000114
Figure BDA0000877580290000115
S0stf0st+Sstfst+gst+lst＝d,ij＝st≠pq (26)
Figure BDA0000877580290000116
Figure BDA0000877580290000117
Figure BDA0000877580290000118
S0pqstf0pqst+Spqstfpqst+gpqst+lpqst＝d,
(30)
ij＝pq＝st且
Figure BDA0000877580290000119
Figure BDA00008775802900001122
(31)
ij＝pq＝st且
Figure BDA00008775802900001110
Figure BDA00008775802900001111
且
Figure BDA00008775802900001112
Figure BDA0000877580290000121
且
Figure BDA0000877580290000122
Figure BDA0000877580290000123
Figure BDA0000877580290000124
Figure BDA0000877580290000125
Figure BDA0000877580290000126
0≤lpq≤d,0≤lst≤d,0≤lpqst≤d (38)
(p,q)∈ψ,(s,t)∈ψ (39)
Eq.(4)-(11) (40)
其中，S0pq为通道p-q中一条线路发生故障下原网络分支节点的关联矩阵；f0pq为通道p-q中一条线路发生故障下原有线路通过的潮流；f0pq为通道p-q中一条线路发生故障下
Figure BDA0000877580290000127
的向量和；Spq为新建线路的网络分支节点的关联矩阵；fpq为通道p-q中新建线路通过的潮流；gpq为通道p-q中一条线路发生故障下节点上发电机的出力；d为每个负荷节点的负荷；
Figure BDA0000877580290000128
为在通道p-q中一条线路发生故障下剩余线路通过的潮流；θi pq和θj pq为在通道p-q中一条线路发生故障下节点i和节点j的电压角度；S0st为通道s-t中一条线路发生故障下原网络分支节点的关联矩阵；fost为通道s-t中一条线路发生故障下原有线路通过的潮流；f0st为通道s-t中一条线路发生故障下
Figure BDA00008775802900001211
的向量和；Sst为新建线路的网络分支节点的关联矩阵；fst为通道s-t中新建线路通过的潮流；gst为通道s-t中一条线路发生故障下节点上发电机的出力；d为每个负荷节点的负荷；
Figure BDA00008775802900001212
为在通道s-t中一条线路发生故障下剩余线路通过的潮流；
Figure BDA0000877580290000129
和
Figure BDA00008775802900001210
为在通道s-t中一条线路发生故障下节点i和节点j的电压角度；S0pqst为通道i-j中两条线路发生故障下原网络分支节点的关联矩阵；f0pqst为通道i-j中两条线路发生故障下原有线路通过的潮流；f0pqst为通道i-j中两条线路发生故障下
Figure BDA00008775802900001213
的向量和；Spqst为新建线路的网络分支节点的关联矩阵；fpqst为通道i-j中新建线路通过的潮流；gpqst为通道i-j中两条线路发生故障下节点上发电机的出力；d为每个负荷节点的负荷；
Figure BDA00008775802900001214
通道i-j中两条线路发生故障下剩余线路通过的潮流；
Figure BDA00008775802900001216
和
Figure BDA00008775802900001215
为在通道i-j中两条线路发生故障下节点i和节点j的电压角度；γij为通道i-j中一条线路的导纳；
Figure BDA00008775802900001217
为通道i-j中原有的线路数量；
Figure BDA00008775802900001219
为通道i-j中一条线路的最大容量；
Figure BDA00008775802900001218
在通道m-n中一条线路发生故障下通道i-j中新建每一条线路通过的潮流；
Figure BDA0000877580290000132
为发电机组的最大功率；
Figure BDA0000877580290000133
为通道i-j中线路的线路能够建设的最大条数；上标pq,st表示在通道p-q和s-t中一回线路发生了故障，pqst表示在同一通道中两回线路同时发生故障，ψ为故障集合；lpq为通道p-q中一回线路发生故障下，电力系统中每个节点损失的电力；
Figure BDA0000877580290000134
为通道i-j中一回线路的最大容量；nij为通道i-j中新建的线路回数；
Figure BDA0000877580290000135
为通道i-j中新建线路中的功率；lst为通道s-t中一回线路发生故障下，电力系统中每个节点损失的电力；
Figure BDA00008775802900001314
为通道i-j中一回线路的最大容量；nij为通道i-j中新建的线路回数；
Figure BDA0000877580290000136
为通道i-j中新建线路中的功率；lpqst为通道i-j中两条线路发生故障下，电力系统中每个节点损失的电力；gpqst为通道i-j中两条线路发生故障下，节点上发电机的出力；
Figure BDA00008775802900001310
为通道i-j两回线路故障时，新建线路中的功率；
Figure BDA00008775802900001311
为通道i-j中一回线路的最大容量；nij为通道i-j中新建的线路回数；
Figure BDA0000877580290000137
为通道i-j中线路没有发生故障时通过的潮流；θi为通道i-j中线路没有发生故障时节点i的电压角度；θj为通道i-j中线路没有发生故障时节点j的电压角度；
Figure BDA0000877580290000138
为通道i-j中线路没有发生故障时通过的潮流；
Figure BDA00008775802900001313
为通道i-j中一条线路的最大容量；
Figure BDA0000877580290000139
为通道i-j中原有线路的回数；
Figure BDA00008775802900001312
为通道i-j中一条线路的最大容量；nij为通道i-j中新建的线路回数；fij为通道i-j中新建线路中的功率；
注意：N-2柔性约束条件表述时，允许有损失负荷l的；而N-1强制性约束条件则不允许出现负荷损失。该约束条件保证了得到的电网规划方案必然满足N-1准则，同时在N-2柔性约束条件下取得了经济性和可靠性的平衡。
步骤4：特定可靠性约束条件的表述
约束条件如下(43):
Lk＜T (41)
式中：Lk(hours/year)为节点k的缺电时间期望；T(hours/year)为LOLE在该约束中的阈值。缺电时间期望的数值可通过负荷持续曲线计算得到，见式(44)：
Figure BDA0000877580290000131
该约束条件保证了电网的可靠性水平，此约束条件的阈值也需要依据国家的相关规定而制订。同时该约束条件也没有考虑故障的发生概率，只考虑故障的影响，因此，避免了某些小概率大影响事故的出现。
步骤5：仿真计算
模型所有的约束条件在寻优计算中都是平等性，寻优仿真计算的步长为一年，即得每一年 输电网拓展的最优规划方案，方案一定满足所有约束条件的。
实施例
本发明以江西省南部输电网络为研究对象，仿真计算其未来十年的中期规划方案，验证所提方法的有效性和经济性。
图2为江西省南部500/220kV电网，该测试系统有22个节点和4回500kV输电线路，每条通道可建最大线路数为3回。表1列出了系统中线路的相关参数。表2显示了每个发电机组的容量和负荷节点的初始负荷，注意：每个负荷都按负荷大小等分成了5个区域，同时赋予了5个购电价格；而每个发电机的出力则为一个整体区域，供电价格则为单一值。表3显示了不同场景的参数，例如：场景1的比重指数为0.4，它表示该场景代表的负荷水平的持续时间为一年(8760小时)中40％；场景1的负荷指数为0.35，其表示该场景代表的负荷水平为该负荷点在当年的最高负荷的35％。图3为系统的有效负荷持续曲线，由表3所给的数据换算得到。其中，节点16、17和21只需要满足N-1准则，不需要进行N-2校验。
表1江西省南部输电网络系统的线路参数
Figure BDA0000877580290000141
Figure BDA0000877580290000151
表2发电机组的容量和负荷节点的初始负荷与电价
Figure BDA0000877580290000152
Figure BDA0000877580290000161
表3不同场景的相关参数
Figure BDA0000877580290000162
Figure BDA0000877580290000171
表4江西省南部输电网络系统的相关参数
Figure BDA0000877580290000172
通过寻优计算得到规划方案如下：在负荷年增速为3％时，混合性规划模型在规划周期内没有线路需要建设，传统规划模型在第10年需要建设一回连接节点5和节点7的线路；在负荷年增长率为6％时，混合型规划模型在第8年需要建设一回连接节点5和节点7的线路，而传统规划模型则需要在第5年建设一回连接节点5至节点7的线路，在第7年建设一回连接节点3至节点4的线路；在负荷年增速为8％时，混合性规划模型要求在第6年和第9年分别建设连接节点5至节点7、节点3至节点4的两回线路，传统规划模型要求在第3年建设一回连接节点5和节点7的线路，第5年建设一回连接节点3至节点4的线路。
表5、表6和图3显示了两个规划模型对应3％，6％和8％负荷增速的社会总效益盈余和系统的严重指数。
表5两种模型对应三个负荷增速的社会总收益盈余(107￥)
Figure BDA0000877580290000173
表6两种模型对应三个负荷增速的系统的严重指数
Figure BDA0000877580290000174
通过以上数据的对比，我们可以发现：混合性规划模型在保证电网可靠性的前提下，合 理延后了连接节点5到7线路的建设时间，该回线路的建设可以使断面(3-4、5-6、5-7)满足N-2准则，同时从整个规划周期看，基于混合性规划模型下的电力系统的可靠性水平和传统规划模型相比没有明显的变化，将潜在的故障损失控制在一个合理的范围内。
以上实施例仅用以说明本发明的技术方案而非对其限制，尽管参照上述实施例对本发明进行了详细的说明，所属领域的普通技术人员依然可以对本发明的具体实施方式进行修改或者等同替换，这些未脱离本发明精神和范围的任何修改或者等同替换，均在申请待批的本发明的权利要求保护范围之内。

Claims (5)

1.一种电力市场环境下中长期输电网拓展规划的混合性规划方法，其特征在于，所述方法包括下述步骤：
步骤1：设定混合性规划模型的目标函数；
步骤2：设定电网运行基本约束条件；
步骤3：利用N-1强制约束条件和N-2柔性约束条件检测电网；
步骤4：判断是否需要特定可靠性约束条件加强输电网以满足N-1强制约束和N-2柔性约束；
步骤5：寻求最优规划方案进行仿真计算；
所述步骤1中，设定混合性规划模型的目标函数满足以下原则：
1)采用资本回收系数将新建线路的投资费用由一个时间点归算到线路寿命的每一年；
2)混合性规划模型的目标函数计算步长为一年，若有投资费用，则投资费用发生在年初，社会总效益和电能损失值在年终计算；
3)采用净现值计算规划周期内的社会总效益和投资费用；一个规划周期中社会总效益和投资费用均以规划周期开始时间为计算净现值的时间点；社会效益盈余等于社会总效益减去投资费用；
4)混合性规划模型通过不同的场景来模拟电力系统不同的负荷水平，同时赋予不同场景相应的比重系数用于衡量对应负荷在一年中的持续时间和负荷情况；
混合性规划模型的目标函数如下：
Figure FDA0002533757310000011
式中：P为特定故障发生概率；r为折现率；a表示场景，ωa为场景a的比重系数；Ωa为所有场景的集合；ΩD为所有负荷点的集合；Ωd为每个负荷所有区域的集合；Q表示发电机；d为每个负荷节点的负荷；h和b分别表示在场景a下负荷d的区域以及在场景a下发电机Q的区域；
Figure FDA0002533757310000012
为在场景a下负荷d在其h区域中的购电价格；
Figure FDA0002533757310000021
为在场景a下负荷d在其h区域中的负荷值；ΩG为所有发电机组的集合；ΩQ为发电机组所有区域的集合；
Figure FDA0002533757310000022
为在场景a下发电机Q在其b区域中的供电价格；
Figure FDA0002533757310000023
为在场景a下发电机Q在其b区域中的出力；
Figure FDA0002533757310000024
为在场景a下负荷d在其h区域中损失单位电量对应的用户可靠性费用；
Figure FDA0002533757310000025
为在场景a下负荷d在其h区域中的损失电力；σ为资本回收系数；Cij为在通道i-j中建设一条线路的总费用；nij为在通道i-j中新建的线路数量；
为计算资本回收系数σ，设新建线路最少运行的寿命为t年，因此资本回收周期为t年，寿命周期t年的折现率为r，表达式为：
Figure FDA0002533757310000026
Figure FDA0002533757310000027
为在场景a下负荷d在其h区域中的损失电力在一年中损失的电量时，采用期望缺供电量EENS进行计算，单位为MWh/year，其数值通过负荷持续曲线计算得到：
Figure FDA0002533757310000028
式中：LPk-l为到达负荷节点k的最大潮流，单位为MW；LPk为节点k最大的用电负荷，单位为MW；φk为负荷节点k上有效负荷持续曲线。
2.如权利要求1所述的混合性规划方法，其特征在于，所述步骤2中，电网运行基本约束条件包括：
功率约束：
Figure FDA0002533757310000029
直流潮流计算：
Figure FDA00025337573100000210
线路潮流约束：
原有线路的潮流约束
Figure FDA00025337573100000211
新建线路的潮流约束
Figure FDA00025337573100000212
发电机机组出力约束：
Figure FDA00025337573100000213
节点上负荷损失值的约束：
0≤La≤d (9)
新建线路回数的约束：
Figure FDA0002533757310000031
其他：
Figure FDA0002533757310000032
式中：Ωa为所有场景的集合；S0为原网络分支节点的关联矩阵；
Figure FDA0002533757310000033
为在场景a下通道i-j中线路通过的潮流；f0a为在场景a下
Figure FDA0002533757310000034
的向量和；S为新建线路的网络分支节点的关联矩阵；fa为在场景a下通道i-j中新建线路通过的潮流，为场景a下
Figure FDA0002533757310000035
的向量和；
Figure FDA0002533757310000036
在场景a下通道i-j中新建每一条线路通过的潮流；ga为在场景a下节点上发电机的出力，即
Figure FDA0002533757310000037
的和，
Figure FDA0002533757310000038
为在场景a下发电机i在其b区域中的出力；La为在场景a节点上的电力损失，即
Figure FDA0002533757310000039
的和，
Figure FDA00025337573100000310
为在场景a下负荷d在其h区域中的损失电力；da为在场景a下节点的负荷值，即
Figure FDA00025337573100000311
的向量和，
Figure FDA00025337573100000312
为在场景a下负荷d在其h区域中的负荷值；γij为通道i-j中一条线路的导纳；
Figure FDA00025337573100000313
为通道i-j中原有的线路数量；
Figure FDA00025337573100000314
为原有网络线路的集合；
Figure FDA00025337573100000315
为通道i-j中一条线路的最大容量；
Figure FDA00025337573100000316
为发电机组的最大功率；d为每个负荷节点的负荷；
Figure FDA00025337573100000317
为通道i-j中线路的线路能够建设的最大条数；nij为非负整数；f0a和
Figure FDA00025337573100000318
为实数。
3.如权利要求1所述的混合性规划方法，其特征在于，所述步骤3中，N-1强制约束条件表述为：
S0mnf0mn+Smnfmn+gmn＝d (12)
Figure FDA00025337573100000319
Figure FDA00025337573100000320
Figure FDA00025337573100000321
Figure FDA0002533757310000041
Figure FDA0002533757310000042
Figure FDA0002533757310000043
Figure FDA0002533757310000044
(m,n)∈ψ (20)
Eq.(4)-(11) (21)
其中，S0mn为通道m-n中一条线路发生故障下原网络分支节点的关联矩阵；f0mn为通道m-n中一条线路发生故障下原有线路通过的潮流；f0mn为通道m-n中一条线路发生故障下
Figure FDA0002533757310000045
的向量和；Smn为新建线路的网络分支节点的关联矩阵；fmn为通道m-n中新建线路通过的潮流；gmn为通道m-n中一条线路发生故障下节点上发电机的出力；d为每个负荷节点的负荷；fij 0mn为在通道m-n中一条线路发生故障下剩余线路通过的潮流；γij为通道i-j中一条线路的导纳；
Figure FDA0002533757310000046
为通道i-j中原有的线路数量；
Figure FDA0002533757310000047
和
Figure FDA0002533757310000048
为在通道m-n中一条线路发生故障下节点i和节点j的电压角度；
Figure FDA0002533757310000049
为通道i-j中一条线路的最大容量；
Figure FDA00025337573100000410
在通道m-n中一条线路发生故障下通道i-j中新建每一条线路通过的潮流；
Figure FDA00025337573100000411
为发电机组的最大功率；
Figure FDA00025337573100000412
为通道i-j中线路的线路能够建设的最大条数；nij为非负整数；ψ为故障集合；
N-2柔性约束条件的表述为：
S0pqf0pq+Spqfpq+gpq+lpq＝d,ij＝pq≠st (22)
Figure FDA00025337573100000413
Figure FDA00025337573100000414
Figure FDA00025337573100000415
S0stf0st+Sstfst+gst+lst＝d,ij＝st≠pq (26)
Figure FDA00025337573100000416
Figure FDA00025337573100000417
Figure FDA00025337573100000418
Figure FDA0002533757310000051
Figure FDA0002533757310000052
Figure FDA0002533757310000053
Figure FDA0002533757310000054
Figure FDA0002533757310000055
Figure FDA0002533757310000056
Figure FDA0002533757310000057
Figure FDA0002533757310000058
0≤lpq≤d,0≤lst≤d,0≤lpqst≤d (38)
(p,q)∈ψ,(s,t)∈ψ (39)
Eq.(4)-(11) (40)
其中，S0pq为通道p-q中一条线路发生故障下原网络分支节点的关联矩阵；
Figure FDA0002533757310000059
为通道p-q中一条线路发生故障下原有线路通过的潮流；f0pq为通道p-q中一条线路发生故障下
Figure FDA00025337573100000510
的向量和；Spq为新建线路的网络分支节点的关联矩阵；fpq为通道p-q中新建线路通过的潮流；gpq为通道p-q中一条线路发生故障下节点上发电机的出力；d为每个负荷节点的负荷；fij 0pq为在通道p-q中一条线路发生故障下剩余线路通过的潮流；
Figure FDA00025337573100000511
和
Figure FDA00025337573100000512
为在通道p-q中一条线路发生故障下节点i和节点j的电压角度；S0st为通道s-t中一条线路发生故障下原网络分支节点的关联矩阵；f0st为通道s-t中一条线路发生故障下
Figure FDA00025337573100000513
的向量和；Sst为新建线路的网络分支节点的关联矩阵；fst为通道s-t中新建线路通过的潮流；gst为通道s-t中一条线路发生故障下节点上发电机的出力；d为每个负荷节点的负荷；fij 0st为在通道s-t中一条线路发生故障下剩余线路通过的潮流；
Figure FDA00025337573100000514
和
Figure FDA00025337573100000515
为在通道s-t中一条线路发生故障下节点i和节点j的电压角度；S0pqst为通道i-j中两条线路发生故障下原网络分支节点的关联矩阵；f0pqst为通道i-j 中两条线路发生故障下原有线路通过的潮流；f0pqst为通道i-j中两条线路发生故障下
Figure FDA0002533757310000061
的向量和；Spqst为新建线路的网络分支节点的关联矩阵；fpqst为通道i-j中新建线路通过的潮流；gpqst为通道i-j中两条线路发生故障下节点上发电机的出力；d为每个负荷节点的负荷；fij 0pqst通道i-j中两条线路发生故障下剩余线路通过的潮流；
Figure FDA0002533757310000062
和
Figure FDA0002533757310000063
为在通道i-j中两条线路发生故障下节点i和节点j的电压角度；γij为通道i-j中一条线路的导纳；
Figure FDA0002533757310000064
为通道i-j中原有的线路数量；
Figure FDA0002533757310000065
为通道i-j中一条线路的最大容量；
Figure FDA0002533757310000066
在通道m-n中一条线路发生故障下通道i-j中新建每一条线路通过的潮流；
Figure FDA0002533757310000067
为发电机组的最大功率；
Figure FDA0002533757310000068
为通道i-j中线路的线路能够建设的最大条数；上标pq,st表示在通道p-q和s-t中一回线路发生了故障，pqst表示在同一通道中两回线路同时发生故障，ψ为故障集合；lpq为通道p-q中一回线路发生故障下，电力系统中每个节点损失的电力；
Figure FDA0002533757310000069
为通道i-j中一回线路的最大容量；nij为通道i-j中新建的线路回数；
Figure FDA00025337573100000610
为通道i-j中新建线路中的功率；lst为通道s-t中一回线路发生故障下，电力系统中每个节点损失的电力；
Figure FDA00025337573100000611
为通道i-j中一回线路的最大容量；nij为通道i-j中新建的线路回数；
Figure FDA00025337573100000612
为通道i-j中新建线路中的功率；lpqst为通道i-j中两条线路发生故障下，电力系统中每个节点损失的电力；gpqst为通道i-j中两条线路发生故障下，节点上发电机的出力；
Figure FDA00025337573100000613
为通道i-j两回线路故障时，新建线路中的功率；
Figure FDA00025337573100000614
为通道i-j中一回线路的最大容量；nij为通道i-j中新建的线路回数；
Figure FDA00025337573100000615
为通道i-j中线路没有发生故障时通过的潮流；θi为通道i-j中线路没有发生故障时节点i的电压角度；θj为通道i-j中线路没有发生故障时节点j的电压角度；
Figure FDA00025337573100000616
为通道i-j中线路没有发生故障时通过的潮流；
Figure FDA00025337573100000617
为通道i-j中一条线路的最大容量；
Figure FDA00025337573100000618
为通道i-j中原有线路的回数；
Figure FDA00025337573100000619
为通道i-j中一条线路的最大容量；nij为通道i-j中新建的线路回数；fij为通道i-j中新建线路中的功率；
模型通过N-2柔性约束条件来计算在N-2故障情况下电力系统的负荷损失， 其允许有损失负荷，但是允许损失负荷的大小由目标函数和其他约束条件确定，实现规划方案在N-2柔性约束条件下取得了经济性和可靠性的平衡；模型使用N-1强制性约束条件校验电网，电网在N-1故障情况下不允许出现负荷损失，保证得到的电网规划方案必然满足N-1准则。
4.如权利要求1所述的混合性规划方法，其特征在于，所述步骤4中，特定可靠性约束条件的表述如下式：
Lk＜T (41)
缺电时间期望的数值通过负荷持续曲线计算得到：
Figure FDA0002533757310000071
式中：Lk为节点k的缺电时间期望，单位为hours/year；T为缺电时间期望(LOLE)在特定可靠性约束中的阈值，单位为hours/year；LPk-l为到达节点k的最大潮流，单位为MW；LPk为节点k最大的用电负荷，单位为MW；φk为节点k上有效负荷持续曲线。
5.如权利要求1所述的混合性规划方法，其特征在于，所述步骤5中，混合性规划模型已对所有的约束条件进行了线性化处理，使用直流模型和最优潮流模型的方法，将输电网拓展规划问题转为混合整数线性规划问题，通过商业软件CPLEX实现仿真计算。
CN201510925496.8A 2015-12-14 2015-12-14 一种电力市场环境下中长期输电网拓展规划的混合性规划方法 Active CN106875026B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510925496.8A CN106875026B (zh) 2015-12-14 2015-12-14 一种电力市场环境下中长期输电网拓展规划的混合性规划方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510925496.8A CN106875026B (zh) 2015-12-14 2015-12-14 一种电力市场环境下中长期输电网拓展规划的混合性规划方法

Publications (2)

Publication Number Publication Date
CN106875026A CN106875026A (zh) 2017-06-20
CN106875026B true CN106875026B (zh) 2020-10-13

Family

ID=59178383

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510925496.8A Active CN106875026B (zh) 2015-12-14 2015-12-14 一种电力市场环境下中长期输电网拓展规划的混合性规划方法

Country Status (1)

Country Link
CN (1) CN106875026B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109523309B (zh) * 2018-10-26 2021-04-30 武汉华中思能科技有限公司 一种基于潮流分析计算的电力报价仿真方法及系统
CN109586280A (zh) * 2018-11-21 2019-04-05 国网冀北电力有限公司 一种机组交易极限电量确定方法及系统
CN113592673B (zh) * 2021-08-05 2024-03-22 国网冀北电力有限公司经济技术研究院 一种电力市场环境下的输电网规划成本效益评估方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103514570A (zh) * 2013-08-14 2014-01-15 国家电网公司 一种含分布式电源的配电网扩展规划综合优化方法
CN104143119A (zh) * 2014-07-29 2014-11-12 华北电力大学 一种多尺度分层蜂窝输电网及其规划方法
CN104732349A (zh) * 2015-03-30 2015-06-24 国家电网公司 一种电网规划方法
CN105141464A (zh) * 2015-09-24 2015-12-09 国网上海市电力公司 一种输配电协调的网架结构规划方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103514570A (zh) * 2013-08-14 2014-01-15 国家电网公司 一种含分布式电源的配电网扩展规划综合优化方法
CN104143119A (zh) * 2014-07-29 2014-11-12 华北电力大学 一种多尺度分层蜂窝输电网及其规划方法
CN104732349A (zh) * 2015-03-30 2015-06-24 国家电网公司 一种电网规划方法
CN105141464A (zh) * 2015-09-24 2015-12-09 国网上海市电力公司 一种输配电协调的网架结构规划方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
采用NSGA_混合智能算法的风电场多目标电网规划;王茜 等;《中国电机工程学报》;20110705;第31卷(第19期);17-24 *

Also Published As

Publication number Publication date
CN106875026A (zh) 2017-06-20

Similar Documents

Publication Publication Date Title
Bialek et al. Benchmarking and validation of cascading failure analysis tools
CN108183512B (zh) 一种接入新能源的电力系统的可靠性评估方法
Li et al. Component importance assessment of power systems for improving resilience under wind storms
CN110969347A (zh) 一种输电网结构形态评估方法
US20060190227A1 (en) Group-based BCU methods for on-line dynamical security assessments and energy margin calculations of practical power systems
CN103825272A (zh) 一种基于解析法的含分散风电的配电网可靠性确定方法
Wu et al. Transmission network expansion planning based on chronological evaluation considering wind power uncertainties
CN106875026B (zh) 一种电力市场环境下中长期输电网拓展规划的混合性规划方法
CN115456304A (zh) 考虑台风影响的海上风电场可靠性指标计算方法及装置
Xiyun et al. Wind power probability interval prediction based on bootstrap quantile regression method
Bellani et al. A reliability-centered methodology for identifying renovation actions for improving resilience against heat waves in power distribution grids
Cao et al. Reliability tracing of the integrated energy system using the improved shapley value
Gharebaghi et al. Inclusion of pre-existing undervoltage load shedding schemes in ac-qss cascading failure models
CN113837812A (zh) 节点电价联合概率预测方法及装置
CN106056305A (zh) 一种基于状态聚类的发电系统可靠性快速评估方法
CN103324795A (zh) 一种考虑站用电影响的直流系统可靠性评估方法
Kiel et al. Transmission line unavailability due to correlated threat exposure
CN106875077B (zh) 一种用于输电网中长期拓展规划的混合性规划方法
CN113177717B (zh) 一种基于影响增量灵敏度的输电系统韧性快速评估方法
CN115964951A (zh) 电力系统运行可靠性评估方法、装置及电子设备
Du et al. Hosting capacity assessment in distribution networks considering wind–photovoltaic–load temporal characteristics
Abbasghorbani Prioritization of transmission network components based on their failure impact on reliability of composite power systems
CN107546742B (zh) 一种日前计划潮流有功及电压计算误差分析方法
Papic et al. Practical experience in evaluating adequacy of generating capacity in the Western interconnection
Hamidovic et al. Load model identification and parameter estimation

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant