CN103985066B - 一种基于混合潮流的电力***静态风险评估方法 - Google Patents
一种基于混合潮流的电力***静态风险评估方法 Download PDFInfo
- Publication number
- CN103985066B CN103985066B CN201410213868.XA CN201410213868A CN103985066B CN 103985066 B CN103985066 B CN 103985066B CN 201410213868 A CN201410213868 A CN 201410213868A CN 103985066 B CN103985066 B CN 103985066B
- Authority
- CN
- China
- Prior art keywords
- load
- power
- formula
- flow
- node
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000012502 risk assessment Methods 0.000 title claims abstract description 16
- 230000003068 static effect Effects 0.000 title claims abstract description 16
- 239000011159 matrix material Substances 0.000 claims description 33
- 230000007257 malfunction Effects 0.000 claims description 30
- 230000035945 sensitivity Effects 0.000 claims description 18
- 238000004364 calculation method Methods 0.000 claims description 13
- 238000004458 analytical method Methods 0.000 claims description 12
- 238000011156 evaluation Methods 0.000 claims description 12
- 238000002347 injection Methods 0.000 claims description 12
- 239000007924 injection Substances 0.000 claims description 12
- 230000009467 reduction Effects 0.000 claims description 9
- 239000013598 vector Substances 0.000 claims description 9
- 230000008859 change Effects 0.000 claims description 8
- 238000005516 engineering process Methods 0.000 claims description 7
- 239000012141 concentrate Substances 0.000 claims description 6
- 239000000203 mixture Substances 0.000 claims description 5
- 239000000243 solution Substances 0.000 claims description 4
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 238000013178 mathematical model Methods 0.000 claims description 3
- 235000013399 edible fruits Nutrition 0.000 claims description 2
- 230000008569 process Effects 0.000 abstract description 3
- 230000008901 benefit Effects 0.000 description 3
- 230000005611 electricity Effects 0.000 description 3
- 230000015556 catabolic process Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000000205 computational method Methods 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种基于混合潮流的电力***静态风险评估方法,包括:(1)根据当前电力***结构和规划方案,建立相应的规划方案模型;(2)根据设备数据及电网结构分析电网预计发生的故障及概率,建立预想故障集;(3)利用基于混合潮流的风险状态计算预想故障集中所有故障对应的潮流,对于存在潮流越界的故障,将负荷损失作为故障后果;(4)根据故障状态的负荷损失与概率计算***风险指标;(5)根据***风险指标寻找***的薄弱环节。本发明利用混合潮流克服了在风险评估过程中,复杂大电网以及高阶故障中单独应用交流潮流法遇到的潮流不收敛问题。
Description
技术领域
本发明属于电网规划领域,并涉及电力***风险评估领域
背景技术
传统电网规划设计研究大多数集中于在满足未来负荷需求的基础上,确定待建输电线路的数量以及位置以使总投资成本最低这样一个问题。其中,在对规划方案进行比较和安全稳定校核时,一般采用确定性的方法对规划方案进行可靠性评估。确定性方法是在预期故障发生的情况下,研究电网的可靠性水平。最常用的确定性方法是***N-1或N-K方法,即在任意一个或K个元件(如发电机、线路、变压器等)故障退出运行时,计算***中线路的潮流和节点电压,校验其是否发生过载或电压越界现象。确定性方法只是对***的安全性作粗略估计,提出改进薄弱环节的措施,但是该方法只是考虑了预想事故发生的后果,忽视了事故发生的概率,从而容易导致对发生后果严重而发生概率低的事件过于悲观,而对发生概率高影响小的事件过于乐观。
将风险评估的理论引入电网规划中,同时考虑故障发生的概率及其产生的后果,来综合地评定电力***的风险,并将风险转化为费用,将电力***的风险水平和风险费用纳入到电网规划设计范围内。
在静态风险评估中,判断***失效状态和估计故障损失后果都依赖于潮流计算程序进行的计算分析。然而在工程实际应用时,由于***规模的庞大和复杂性,以及大量直流输电线路接入交流电网运行,使得现有交流潮流计算方法的收敛性很难得到保障。尤其是在基于负荷水平预测的电网发展规划、或者是针对高阶故障进行分析时,往往出现潮流计算不收敛导致风险后果无法评估的问题。
复杂大电网的潮流计算收敛性问题一直是世界性的难题,至今也并无有效的方法能够保证任意运行方式下的电网潮流计算收敛。然而从电力***风险评估的要求考虑,实际并不需要十分精确的潮流计算结果,因此可以将交流潮流法和直流潮流法相结合,用于计算电力***的潮流。
直流潮流法虽然不能计算节点电压的幅值,有功潮流也存在部分误差,但与采用牛顿-拉夫逊法或P-Q解耦法的交流潮流比较,其有功功率的平均误差不大(对于高压电网,误差一般在3-5%左右)。直流潮流求解时无需迭代,占用内存少,计算时间大大降低,潮流的收敛性也优于交流潮流。在计算故障后***的潮流时,对于交流潮流计算发散的预想事故,采用直流潮流方法,近似计算得到电网中各线路的有功潮流;假定功率因数不变,则可进一步基于有功潮流折算得到各线路的视在功率潮流,再结合直流最优潮流便可计算出消除越限状态所需的最小切负荷量,并将其作为故障损失的后果。
发明内容
针对传统电网规划评估方法的不足,本发明提供一种基于混合潮流的电力***静态风险评估方法,利用风险评估技术,结合故障发生的概率与其产生的后果严重性,综合地对电力***安全水平做出评价,并且在潮流计算过程中,结合交流潮流法和直流潮流法的优点,利用混合潮流克服在复杂大电网以及高阶故障中单独应用交流潮流法遇到的潮流不收敛问题。
为了解决上述技术问题,本发明一种基于混合潮流的电力***静态风险评估方法,包括以下步骤:
步骤一、根据当前电力***结构和未来需求制定规划方案,并建立相应的规划方案模型;
步骤二、根据设备数据及电网结构分析电网预计发生的故障及其故障概率,建立预想故障集;
步骤三、利用基于混合潮流的风险状态计算预想故障集中所有故障对应的潮流,对于存在潮流越界的故障,利用基于灵敏度的最优安全控制算法计算用于消除越限状态所需的切机切负荷量,该负荷损失即为故障造成的后果;具体内容包括:
步骤1)从预扫描故障集中提取故障状态fk,利用交流潮流程序计算该故障状态fk下的***潮流;
步骤2)判断步骤1)的交流潮流计算结果是否收敛,若收敛,则执行步骤4),否则执行步骤3);
步骤3)利用直流潮流计算故障状态fk下的***直流潮流,假设:线路电阻比线路电抗小10倍以上,则线路i-j的电纳为:
式(1)中,xij为线路i-j的电抗,
线路节点i与线路节点j之间的电压相角差δij小于10°,则:
sinδij≈δi-δj (2)
cosδij≈1 (3)
式(2)中,δi为节点i的电压相角,δj为节点j的电压相角,
线路节点i与线路节点j的对地电纳bi0与bj0忽略不计,即:
bi0=bj0≈0 (4)
所有节点的电压幅值的标幺值假设为1,则一条线路的线路潮流Pij是:
节点注入的有功功率表达式为:
P=B′δ (6)
式(6)中,矩阵B′为n-1阶方矩阵,n为***节点数,矩阵B′中的元素 P为n-1阶列向量,其元素为除平衡节点外各节点的有功功率注入量;δ为n-1阶列向量,其元素为除平衡节点外各节点的电压相角;
合并式(5)和式(6)得到节点注入有功功率和线路潮流间的线性关系:
Tp=AP (7)
式(7)中,Tp为线路潮流向量,其元素是线路潮流Pij;矩阵A是节点注入有功功率与线路潮流之间的关系矩阵,矩阵A的维数为L×(n-1),其中,L表示故障后***的线路数;由矩阵B′直接计算矩阵A,假设线路l的两个节点编号为i和j,当l=1,…,L时,矩阵A中的第l行元素由以下方程组解出:
B′Al=C (8)
式(8)中,1/xij为第i个元素,-1/xij为第j个元素;
步骤4)判断故障状态fk下的***潮流是否存在越限,若存在,则该故障状态为失效状态,计算为消除潮流越限所需的切机切负荷量,并将负荷损失作为故障造成的后果;
若失效状态fk的***潮流是通过交流潮流程序计算得到的,则需要通过静态安全最优控制模型来计算预想故障所造成的负荷损失,将发电机和负荷的控制费用最小值作为目标函数,利用发电机有功功率和负荷有功功率作为控制手段,引入静态灵敏度技术将非线性静态安全最优控制模型转化成线性静态安全最优控制模型,此线性模型如下所示:
式(9)表示对某一故障的最小控制成本;式(10)表示越限线路ij的潮流约束;式(11)表示平衡节点有功功率平衡约束;式(12)表示发电机有功功率变化上下限;式(13)表示负荷有功功率变化上下限;
式(9)至式(13)中:
y、z分别为发电机节点个数和负荷节点个数;
ΔPga、ΔPLa分别为发电机和负荷有功功率变化量;
Cga、CLa分别为发电机和负荷的控制费用;
i、j为潮流越限线路ij首末端节点编号;
分别为线路ij上电流对发电机a有功功率和负荷a有功功率的灵敏度系数,对于变压器线路可以采用视在功率的灵敏度系数;
分别为平衡节点s有功功率对发电机a有功功率和负荷a有功功率的灵敏度系数;
和各灵敏度系数均通过对潮流方程求偏导数得到;
该最优控制模型是典型的线性规划问题,采用单纯型法进行求解,将Cga设置为0,即对于某一失效状态,优先进行发电机调整,若只通过发电调整无法消除越限,再进行负荷削减操作;利用线性静态安全最优控制模型计算所得最优解ΔPga与ΔPLa即为消除潮流越限所需的切机切负荷量;
若失效状态fk的***潮流是通过直流潮流计算得到的,利用基于直流潮流的最优潮流模型近似估计故障造成的负荷损失量,其数学模型如下所示:
s.t.Tk=Ak(Pg-PL+△PLa) (15)
|Tk|≤Tmax (19)
式(17)表示发电机有功功率的上下限,式(18)表示负荷有功功率变化上下限,式(19)表示***潮流的限制;
式(14)至式(19)中:
y、z分别为发电机节点个数和负荷节点个数;
ΔPLa为负荷有功功率变化量;
CLa为负荷损失对应的费用,
Tk是故障后***的有功潮流矩阵;
Pg和PL分别发电机输出和负荷功率矩阵;
Ak为故障状态fk所对应的有功潮流和注入功率间的关系矩阵;
步骤5)判断是否完成预想故障集中所有故障状态的混合潮流计算及失效状态分析,如果完成,则执行步骤6),否则令k=k+1,并返回步骤1);
步骤6)存储分析结果,以便后续操作。
步骤四、利用步骤三的计算的切机切负荷量和相应的故障概率计算***风险指标;
步骤五、根据***风险指标对电力***进行综合的评价,从而寻找***的薄弱节点与薄弱线路。
本发明中,若在步骤一中建立了多个规划方案,将利用***风险指标进行规划方案间的横向比较,最终确定出最优的规划方案。
与现有技术相比,本发明的有益效果是:
本发明在电网安全性评估中引入了风险理论,将故障发生的概率与故障产生的后果相结合,综合的对***风险水平做出评价,避免了对后果严重而发生概率低的事件做出过于悲观的评价,而对发生概率高影响小的事件过于乐观。本风险评估方法基于混合潮流理论,结合了交流潮流与直流潮流的优点,有效的避免了交流潮流计算过程中结果不收敛的问题。
附图说明
图1是本发明提供的基于混合潮流的电力***静态风险评估方法整体实施流程图;
图2是本发明提供的基于混合潮流的故障扫描方法流程图;
图3是实施算例***的结构示意图,该实施算例中所用节点均已圈出,并用英文字母进行编号。
具体实施方式
下面结合附图和具体实施例对本发明技术方案作进一步详细描述。
本发明提出了一种基于混合潮流的电力***静态风险评估技术,其实施流程图如图1和图2所示,现以图3所示的算例网络为例,详细说明如下:
步骤一、根据规划区当前电力***结构和未来需求制定规划方案,并建立相应的规划方案模型;实施算例电网预计未来用电需求及装机容量如表1所示,按照预测结果制定规划方案如图3所示。
表1电网用电需求及装机容量预测结果表
步骤二、根据设备数据及电网结构分析电网预计发生的故障及其故障概率,建立预想故障集;根据规划方案模型中的设备数据及电网结构,利用状态枚举法建立相应的预想故障集F={(fk,pk)|k=1,2…N},并另k=1;
根据图3所示规划方案,考虑电网所有500kv线路故障,并扫描到3阶故障,利用状态枚举法建立相应的预想故障集,该故障集包括很多个故障状态,部分故障状态如表2所示。由于本算例***网络较大,仅将算例中用到的节点在图3中用英文简写进行标注。
表2部分故障状态及其故障概率
步骤三、利用基于混合潮流的风险状态计算预想故障集中所有故障对应的潮流,对于存在潮流越界的故障,利用基于灵敏度的最优安全控制算法计算用于消除越限状态所需的切机切负荷量,该负荷损失即为故障造成的后果;
该步骤三的具体内容包括:
步骤1)从预扫描故障集中提取故障状态fk,利用交流潮流程序计算该故障状态fk下的***潮流;
步骤2)判断步骤1)的交流潮流计算结果是否收敛,若收敛,则执行步骤4),否则执行步骤3);
针对图3所示的规划方案,由交流潮流计算结果的收敛情况决定故障状态的潮流计算类型,部分故障状态的潮流计算情况如表3所示:
表3规划方案故障状态对应潮流计算情况
步骤3)利用直流潮流计算故障状态fk下的***直流潮流,假设:线路电阻比线路电抗小得多,小10倍以上,则线路i-j的电纳为:
式(1)中,xij为线路i-j的电抗,
线路节点i与线路节点j之间的电压相角差δij小于10°,则:
sinδij≈δi-δj (2)
cosδij≈1 (3)
式(2)中,δi为节点i的电压相角,δj为节点j的电压相角,
线路节点i与线路节点j的对地电纳bi0与bj0忽略不计,即:
bi0=bj0≈0 (4)
所有节点的电压幅值的标幺值假设为1,则一条线路的线路潮流Pij是:
节点注入的有功功率表达式为:
P=B′δ (6)
式(6)中,矩阵B′为n-1阶方矩阵,n为***节点数,矩阵B′中的元素 P为n-1阶列向量,其元素为除平衡节点外各节点的有功功率注入量;δ为n-1阶列向量,其元素为除平衡节点外各节点的电压相角;
合并式(5)和式(6)得到节点注入有功功率和线路潮流间的线性关系:
Tp=AP (7)
式(7)中,Tp为线路潮流向量,其元素是线路潮流Pij;矩阵A是节点注入有功功率与线路潮流之间的关系矩阵,矩阵A的维数为L×(n-1),其中,L表示故障后***的线路数;由矩阵B′直接计算矩阵A,假设线路l的两个节点编号为i和j,当l=1,…,L时,矩阵A中的第l行元素由以下方程组解出:
B′Al=C (8)
式(8)中,1/xij为第i个元素,-1/xij为第j个元素;
步骤4)判断故障状态fk下的***潮流是否存在越限,若存在,则该故障状态为失效状态,计算为消除潮流越限所需的切机切负荷量,并将负荷损失作为故障造成的后果;
若失效状态fk的***潮流是通过交流潮流程序计算得到的,则需要通过静态安全最优控制模型来计算预想故障所造成的负荷损失,将发电机和负荷的控制费用最小值作为目标函数,利用发电机有功功率和负荷有功功率作为控制手段,引入静态灵敏度技术将非线性静态安全最优控制模型转化成线性静态安全最优控制模型,此线性模型如下所示:
式(9)表示对某一故障的最小控制成本;式(10)表示越限线路ij的潮流约束;式(11)表示平衡节点有功功率平衡约束;式(12)表示发电机有功功率变化上下限;式(13)表示负荷有功功率变化上下限;
式(9)至式(13)中:
y、z分别为发电机节点个数和负荷节点个数;
ΔPga、ΔPLa分别为发电机和负荷有功功率变化量;
Cga、CLa分别为发电机和负荷的控制费用;
i、j为潮流越限线路ij首末端节点编号;
分别为线路ij上电流对发电机a有功功率和负荷a有功率的灵敏度系数,对于变压器线路可以采用视在功率的灵敏度系数;
分别为平衡节点s有功功率对发电机a有功功率和负荷a有功功率的灵敏度系数;
和各灵敏度系数均通过对潮流方程求偏导数得到;
该最优控制模型是典型的线性规划问题,采用单纯型法进行求解,将Cga设置为0,即对于某一失效状态,优先进行发电机调整,若只通过发电调整无法消除越限,再进行负荷削减操作;利用线性静态安全最优控制模型计算所得最优解ΔPga与ΔPLa即为消除潮流越限所需的切机切负荷量;
若失效状态fk的***潮流是通过直流潮流计算得到的,利用基于直流潮流的最优潮流模型近似估计故障造成的负荷损失量,其数学模型如下所示:
s.t.Tk=Ak(Pg-PL+△PLa) (15)
|Tk|≤Tmax (19)
式(17)表示发电机有功功率的上下限,式(18)表示负荷有功功率变化上下限,式(19)表示***潮流的限制;
式(14)至式(19)中:
y、z分别为发电机节点个数和负荷节点个数;
ΔPLa为负荷有功功率变化量;
CLa为负荷损失对应的费用,
Tk是故障后***的有功潮流矩阵;
Pg和PL分别发电机输出和负荷功率矩阵;
Ak为故障状态fk所对应的有功潮流和注入功率间的关系矩阵;
针对图3所示的规划方案,对于存在越限情况的故障状态,在进行切机切负荷操作时,优先切除发电机,通过发电调整来消除越限,若仅仅通过切机操作无法消除越限,那么再进行切负荷操作。***负荷损失状况如表4所示。
表4***负荷损失情况
步骤5)判断是否完成预想故障集中所有故障状态的混合潮流计算及失效状态分析,如果完成,则执行步骤6),否则令k=k+1,并返回步骤1);
步骤6)存储分析结果,以便后续操作。
步骤四、利用步骤三的计算的切机切负荷量和相应的故障概率计算***风险指标;
本发明方法中用到的风险指标如下:
1)负荷削减概率PLC(probability of load curtailments)
负荷削减概率是指在所研究的时期内出现负荷削减的时间概率总和。电力***会因多种原因导致不得不切除部分负荷,PLC指标则表征了***出现负荷削减状态的总概率,其计算公式如下:
其中,S为有负荷削减的***状态集合。
2)负荷削减频率EFLC(expected frequency of load curtailments)
负荷削减频率则从频率角度来描述***的负荷削减状态,它表征了电力***在一年的时间内发生负荷削减的次数,该指标可由下式计算:
其中,m为元件总数,λi为第i个元件离开状态k的转移率。
3)负荷削减期望持续时间EDLC(expected duration of load curtailments)
将负荷削减概率指标乘以周期小时数,可得到负荷削减的周期时长。若周期为一年,则EDLC称为负荷削减年度小时数,它表征***在一年内发生负荷削减的总期望时长。
EDLC=PLC×8760(小时/年) (22)
4)负荷削减平均持续时间ADLC(average duration of load curtailments)
将负荷削减的周期持续总时间除以周期负荷削减次数,得到的值即为每次负荷削减状态的期望持续小时数,如下式所示。
5)期望缺供电量EENS(expected energy not supplied)
期望缺供电量EENS是计算负荷损失的另一个重要指标,它表征在一年内,由于停运所造成电量损失的平均值,其表达式为:
其中,ΔPLk为故障状态fk对应的负荷削减量。由于EENS是能量指标,对可靠性经济评估、最优可靠性、***规划等均具有重要意义,因此EENS是风险评估中非常重要的指标。
7)严重程度指标SI(severity index)
SI=EENS×60/L(***分) (25)
其中,1个***分相当于在最大负荷时全***停电1分钟,是对***故障的严重程度的一种度量。
图3所示***的风险指标如表5所示。
表5***风险指标
PLC | EENS | EFLC | EDLC | SI |
6.82E-08 | 1.81 | 1.13E-04 | 5.30 | 3.17E-04 |
步骤五、根据***风险指标对电力***进行综合的评价,从而寻找***的薄弱节点与薄弱线路。
由风险指标可以看出,***的SI指标很小,***安全水平很高。虽然该规划方案的充裕度很高,但仍存在薄弱环节。通过对***各线路、节点风险指标进行分析,可以看出,由于MM区域通过双回线与CH和BJ相联系,当BJ-MM双回线故障、CH-MM双回线一回故障,或者CH-MM双回线故障、BJ-MM双回线一回故障时,均会造成MM区域仅剩一条线路与外部联系。MM区负荷较重,仅靠一条线路不足以满足其负荷要求,因此必然造成该区域负荷损失,这一点在规划设计过程中应当予以重视。
通过实施算例,可以得出,本发明利用风险评估技术,结合故障发生的概率与其产生的后果严重性,能够综合地对电力***安全水平做出评价,并且在潮流计算过程中,结合交流潮流法和直流潮流法的优点,利用混合潮流克服了在复杂大电网以及高阶故障中单独应用交流潮流法遇到的潮流不收敛问题。
尽管上面结合附图对本发明进行了描述,但是本发明并不局限于上述的具体实施方式,上述的具体实施方式仅仅是示意性的,而不是限制性的,本领域的普通技术人员在本发明的启示下,在不脱离本发明宗旨的情况下,还可以做出很多变形,这些均属于本发明的保护之内。
Claims (2)
1.一种基于混合潮流的电力***静态风险评估方法,其特征在于,该方法包括以下步骤:
步骤一、根据当前电力***结构和未来需求制定规划方案,并建立相应的规划方案模型;
步骤二、根据设备数据及电网结构分析电网预计发生的故障及其故障概率,建立预想故障集;
步骤三、利用基于混合潮流的风险状态计算预想故障集中所有故障对应的潮流,对于存在潮流越界的故障,利用基于灵敏度的最优安全控制算法计算用于消除越限状态所需的切机切负荷量,该负荷损失即为故障造成的后果;具体内容包括:
步骤1)从预扫描故障集中提取故障状态fk,利用交流潮流程序计算该故障状态fk下的***潮流;
步骤2)判断步骤1)的交流潮流计算结果是否收敛,若收敛,则执行步骤4),否则执行步骤3);
步骤3)利用直流潮流计算故障状态fk下的***直流潮流,假设:线路电阻比线路电抗小10倍以上,则线路i-j的电纳为:
式(1)中,xij为线路i-j的电抗,
线路节点i与线路节点j之间的电压相角差δij小于10°,则:
sinδij≈δi-δj (2)
cosδij≈1 (3)
式(2)中,δi为节点i的电压相角,δj为节点j的电压相角,
线路节点i与线路节点j的对地电纳bi0与bj0忽略不计,即:
bi0=bj0≈0 (4)
所有节点的电压幅值的标幺值假设为1,则一条线路的线路潮流Pij是:
节点注入的有功功率表达式为:
P=B′δ (6)
式(6)中,矩阵B′为n-1阶方矩阵,n为***节点数,矩阵B′中的元素 P为n-1阶列向量,其元素为除平衡节点外各节点的有功功率注入量;δ为n-1阶列向量,其元素为除平衡节点外各节点的电压相角;
合并式(5)和式(6)得到节点注入有功功率和线路潮流间的线性关系:
Tp=AP (7)
式(7)中,Tp为线路潮流向量,其元素是线路潮流Pij;矩阵A是节点注入有功功率与线路潮流之间的关系矩阵,矩阵A的维数为L×(n-1),其中,L表示故障后***的线路数;由矩阵B′直接计算矩阵A,假设线路l的两个节点编号为i和j,当l=1,…,L时,矩阵A中的第l行元素由以下方程组解出:
B′Al=C (8)
式(8)中,1/xij为第i个元素,-1/xij为第j个元素;
步骤4)判断故障状态fk下的***潮流是否存在越限,若存在,则该故障状态为失效状态,计算为消除潮流越限所需的切机切负荷量,并将负荷损失作为故障造成的后果;
若失效状态fk的***潮流是通过交流潮流程序计算得到的,则需要通过静态安全最优控制模型来计算预想故障所造成的负荷损失,将发电机和负荷的控制费用最小值作为目标函数,利用发电机有功功率和负荷有功功率作为控制手段,引入静态灵敏度技术将非线性静态安全最优控制模型转化成线性静态安全最优控制模型,此线性模型如下所示:
式(9)表示对某一故障的最小控制成本;式(10)表示越限线路ij的潮流约束;式(11)表示平衡节点有功功率平衡约束;式(12)表示发电机有功功率变化上下限;式(13)表示负荷有功功率变化上下限;
式(9)至式(13)中:
y、z分别为发电机节点个数和负荷节点个数;
ΔPga、ΔPLa分别为发电机和负荷有功功率变化量;
Cga、CLa分别为发电机和负荷的控制费用;
i、j为潮流越限线路ij首末端节点编号;
分别为线路ij上电流对发电机a有功功率和负荷a有功功率的灵敏度系数,对于变压器线路可以采用视在功率的灵敏度系数;
分别为平衡节点s有功功率对发电机a有功功率和负荷a有功功率的灵敏度系数;
和各灵敏度系数均通过对潮流方程求偏导数得到;
该最优控制模型是典型的线性规划问题,采用单纯型法进行求解,将Cga设置为0,即对于某一失效状态,优先进行发电机调整,若只通过发电调整无法消除越限,再进行负荷削减操作;利用线性静态安全最优控制模型计算所得最优解ΔPga与ΔPLa即为消除潮流越限所需的切机切负荷量;
若失效状态fk的***潮流是通过直流潮流计算得到的,利用基于直流潮流的最优潮流模型近似估计故障造成的负荷损失量,其数学模型如下所示:
s.t. Tk=Ak(Pg-PL+△PLa) (15)
|Tk|≤Tmax (19)
式(17)表示发电机有功功率的上下限,式(18)表示负荷有功功率变化上下限,式(19)表示***潮流的限制;
式(14)至式(19)中:
y、z分别为发电机节点个数和负荷节点个数;
ΔPLa为负荷有功功率变化量;
CLa为负荷损失对应的费用,
Tk是故障后***的有功潮流矩阵;
Pg和PL分别发电机输出和负荷功率矩阵;
Ak为故障状态fk所对应的有功潮流和注入功率间的关系矩阵;
步骤5)判断是否完成预想故障集中所有故障状态的混合潮流计算及失效状态分析,如果完成,则执行步骤6),否则令k=k+1,并返回步骤1);
步骤6)存储分析结果,以便后续操作;
步骤四、利用步骤三的计算的切机切负荷量和相应的故障概率计算***风险指标;
步骤五、根据***风险指标对电力***进行综合的评价,从而寻找***的薄弱节点与薄弱线路。
2.根据权利要求1所述基于混合潮流的电力***静态风险评估方法,其中,步骤一中建立多个规划方案,利用***风险指标进行规划方案间的横向比较,最终确定出最优的规划方案。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410213868.XA CN103985066B (zh) | 2014-05-20 | 2014-05-20 | 一种基于混合潮流的电力***静态风险评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410213868.XA CN103985066B (zh) | 2014-05-20 | 2014-05-20 | 一种基于混合潮流的电力***静态风险评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103985066A CN103985066A (zh) | 2014-08-13 |
CN103985066B true CN103985066B (zh) | 2017-05-03 |
Family
ID=51277024
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410213868.XA Active CN103985066B (zh) | 2014-05-20 | 2014-05-20 | 一种基于混合潮流的电力***静态风险评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103985066B (zh) |
Families Citing this family (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104376506B (zh) * | 2014-11-17 | 2017-06-06 | 三峡大学 | 一种基于确定性分析的电力***风险评估方法 |
CN104504246B (zh) * | 2014-12-05 | 2017-11-21 | 清华大学 | 基于环‑辐网解耦的可靠性快速评估算法 |
CN105022930A (zh) * | 2015-08-10 | 2015-11-04 | 国网上海市电力公司 | 基于状态枚举法的输电***风险评估方法 |
CN105046026B (zh) * | 2015-08-28 | 2019-02-12 | 南方电网科学研究院有限责任公司 | 一种计及节点功率注入和电网拓扑的输电能力评估方法 |
CN105186513B (zh) * | 2015-10-16 | 2017-05-17 | 山东大学 | 一种基于时序最优潮流法的长期电压稳定预防控制方法 |
CN106651626A (zh) * | 2015-10-29 | 2017-05-10 | 中国电力科学研究院 | 一种判定电网可靠性水平的动态可靠性概率指标确定方法 |
CN105470951B (zh) * | 2015-12-10 | 2017-10-31 | 浙江大学 | 一种基于大数据挖掘技术的配电网可靠供电能力评估方法 |
CN105654245B (zh) * | 2015-12-31 | 2019-11-29 | 太原理工大学 | 基于负荷不确定性建模的电网静态电压稳定风险评估方法 |
CN105576655B (zh) * | 2016-01-13 | 2018-08-24 | 国网山西省电力公司大同供电公司 | 一种基于可靠性价值增量的关键线路识别方法和*** |
CN105741016B (zh) * | 2016-01-22 | 2020-10-27 | 中国电力科学研究院 | 一种用于中期电网规划的静态可靠性概率指标获取方法 |
CN107221945B (zh) * | 2017-05-15 | 2019-09-27 | 中国电力科学研究院 | 一种特高压直流线路预想故障辅助决策方法和装置 |
CN107196300B (zh) * | 2017-07-18 | 2020-08-18 | 西安交通大学 | 一种考虑事故负荷削减量的变压器重要程度评估方法 |
CN109038550B (zh) * | 2018-07-16 | 2022-01-18 | 贵州电网有限责任公司 | 基于电压静态稳定性的电力***自复性指标计算方法 |
CN109447512B (zh) * | 2018-11-13 | 2021-12-03 | 重庆大学 | 基于均匀设计的大电网可靠性评估方法 |
CN112769124B (zh) * | 2019-11-05 | 2023-09-05 | 华北电力大学(保定) | 一种基于潮流转移和追踪的电力***快速运行风险评估方法 |
CN110969355B (zh) * | 2019-12-03 | 2023-06-09 | 重庆大学 | 增量风险事件的筛选方法、装置及计算机可读介质 |
CN111652469B (zh) * | 2020-04-29 | 2022-11-08 | 国网河南省电力公司电力科学研究院 | 一种交直流混联电网的薄弱环节辨识方法和*** |
CN112734209B (zh) * | 2020-12-31 | 2022-07-29 | 国网山东省电力公司青岛供电公司 | 负荷密集型城市电网安全稳定控制风险预警方法及*** |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012109548A2 (en) * | 2011-02-10 | 2012-08-16 | Eaton Corporation | Equipment-related risk assessment from electrical event analysis |
CN102737286A (zh) * | 2012-04-23 | 2012-10-17 | 中国电力科学研究院 | 一种地区电网在线风险分析***和方法 |
CN102738791A (zh) * | 2012-06-12 | 2012-10-17 | 中国电力科学研究院 | 一种基于拓扑比较的在线静态安全分析方法 |
CN103218754A (zh) * | 2013-04-23 | 2013-07-24 | 广东电网公司电力调度控制中心 | 一种电网调度操作的风险检测方法及装置 |
CN103279807A (zh) * | 2013-05-06 | 2013-09-04 | 国家电网公司 | 一种恶劣天气下电网静态风险评估方法 |
-
2014
- 2014-05-20 CN CN201410213868.XA patent/CN103985066B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012109548A2 (en) * | 2011-02-10 | 2012-08-16 | Eaton Corporation | Equipment-related risk assessment from electrical event analysis |
CN102737286A (zh) * | 2012-04-23 | 2012-10-17 | 中国电力科学研究院 | 一种地区电网在线风险分析***和方法 |
CN102738791A (zh) * | 2012-06-12 | 2012-10-17 | 中国电力科学研究院 | 一种基于拓扑比较的在线静态安全分析方法 |
CN103218754A (zh) * | 2013-04-23 | 2013-07-24 | 广东电网公司电力调度控制中心 | 一种电网调度操作的风险检测方法及装置 |
CN103279807A (zh) * | 2013-05-06 | 2013-09-04 | 国家电网公司 | 一种恶劣天气下电网静态风险评估方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103985066A (zh) | 2014-08-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103985066B (zh) | 一种基于混合潮流的电力***静态风险评估方法 | |
CN103985065B (zh) | 一种基于故障预扫描的电力***风险评估方法 | |
Shao et al. | Fast evaluation of commutation failure risk in multi-infeed HVDC systems | |
Williams et al. | Probabilistic load flow modeling comparing maximum entropy and Gram-Charlier probability density function reconstructions | |
CN105406470B (zh) | 一种基于开关边界分区的有源配电网可靠性评估方法 | |
Cao et al. | Chance-constrained optimization-based unbalanced optimal power flow for radial distribution networks | |
CN102637289A (zh) | 含大规模风电的电力***规划方案的安全性价值评估方法 | |
CN101592700A (zh) | 基于事故链的大电网连锁故障分析方法 | |
CN102222170B (zh) | 一种换流变***可靠性评估的马尔可夫状态空间图方法 | |
CN103366220A (zh) | 电力***的运行风险评估方法 | |
CN104599189A (zh) | 一种计入电力***运行方式的电网规划方案风险评估方法 | |
CN107181253B (zh) | 一种基于电网动态可靠性概率指标的电网规划方法 | |
CN103995982A (zh) | 一种考虑机组随机故障的概率潮流计算方法 | |
CN104899798B (zh) | 一种考虑旋转备用的风电接入***的暂态风险控制方法 | |
CN102570450A (zh) | 一种适用于复杂电网的静态可靠性评估方法 | |
CN104037760B (zh) | 一种不确定注入功率的电力***预想事故选择方法 | |
CN104617575B (zh) | 一种用于静态安全分析预想故障下潮流不可解评估方法 | |
Su et al. | Novel voltage-to-power sensitivity estimation for phasor measurement unit-unobservable distribution networks based on network equivalent | |
CN104868478A (zh) | 一种电网紧急状态下的启动动态分区方案的方法 | |
CN104578050B (zh) | 一种电网暂态稳定强相关输电断面的识别方法 | |
CN109767113A (zh) | 一种基于决策树的动态安全评估及预防控制综合方法 | |
CN101841154A (zh) | 输电网严重故障后的电压稳定裕度实时评估与最优控制方法 | |
CN104283208B (zh) | 大型电网概率性可用输电能力的分解协调计算方法 | |
CN110994614A (zh) | 一种直流输电***的可靠性评估方法、***及设备 | |
Abdelghany et al. | Reliability analysis of smart grids using formal methods |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |