CN110601185B - 一种基于统一潮流模型的综合能源***薄弱点辨识方法 - Google Patents
一种基于统一潮流模型的综合能源***薄弱点辨识方法 Download PDFInfo
- Publication number
- CN110601185B CN110601185B CN201910876207.8A CN201910876207A CN110601185B CN 110601185 B CN110601185 B CN 110601185B CN 201910876207 A CN201910876207 A CN 201910876207A CN 110601185 B CN110601185 B CN 110601185B
- Authority
- CN
- China
- Prior art keywords
- matrix
- comprehensive energy
- energy system
- rate
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 52
- 239000011159 matrix material Substances 0.000 claims abstract description 122
- 230000002159 abnormal effect Effects 0.000 claims abstract description 24
- 230000008878 coupling Effects 0.000 claims description 14
- 238000010168 coupling process Methods 0.000 claims description 14
- 238000005859 coupling reaction Methods 0.000 claims description 14
- 238000001228 spectrum Methods 0.000 claims description 11
- 238000007781 pre-processing Methods 0.000 claims description 9
- 230000003416 augmentation Effects 0.000 claims description 8
- 238000010606 normalization Methods 0.000 claims description 8
- 230000003595 spectral effect Effects 0.000 claims description 8
- OIGNJSKKLXVSLS-VWUMJDOOSA-N prednisolone Chemical compound O=C1C=C[C@]2(C)[C@H]3[C@@H](O)C[C@](C)([C@@](CC4)(O)C(=O)CO)[C@@H]4[C@@H]3CCC2=C1 OIGNJSKKLXVSLS-VWUMJDOOSA-N 0.000 claims description 7
- 238000005516 engineering process Methods 0.000 claims description 6
- 238000000926 separation method Methods 0.000 claims description 6
- 230000003190 augmentative effect Effects 0.000 claims description 5
- 238000011156 evaluation Methods 0.000 claims description 4
- 230000011218 segmentation Effects 0.000 claims description 4
- 238000002940 Newton-Raphson method Methods 0.000 claims description 3
- 238000010438 heat treatment Methods 0.000 claims description 3
- 238000010223 real-time analysis Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 abstract description 10
- 239000007789 gas Substances 0.000 description 20
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 8
- 230000003321 amplification Effects 0.000 description 6
- 238000010276 construction Methods 0.000 description 6
- 238000003199 nucleic acid amplification method Methods 0.000 description 6
- 230000004044 response Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000012986 modification Methods 0.000 description 5
- 230000004048 modification Effects 0.000 description 5
- 239000003345 natural gas Substances 0.000 description 4
- 238000010248 power generation Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 238000007405 data analysis Methods 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 238000012163 sequencing technique Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000008929 regeneration Effects 0.000 description 1
- 238000011069 regeneration method Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000010206 sensitivity analysis Methods 0.000 description 1
- 230000000087 stabilizing effect Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于统一潮流模型和随机矩阵的综合能源***薄弱点辨识方法,包括以下步骤:(1)建立综合能源***统一潮流模型;(2)选取潮流数据并建立随机矩阵;(3)标准化预处理随机矩阵;(4)识别***异常状态;(5)辨识***薄弱点。本发明的方法可以。本发明将随机矩阵理论引入综合能源***薄弱辨识中,以数据驱动的方式分析***运行状态,不需要对综合能源***物理结构进行辨识,避免了对综合能源***的建模过程,克服了综合能源***物理结构和建模过程的复杂性,提高***状态识别和薄弱点辨识的准确性和快速性。
Description
技术领域
本发明涉及综合能源***薄弱环节辨识技术领域,具体涉及一种基于统一潮流模型的综合能源***薄弱点辨识方法。
背景技术
随着化石能源日益稀缺和环境恶化,运用综合能源***实现对电、热、气等多种能源的一体化规划设计和运行优化,提高能源的综合利用效率,已经成为人类社会能源领域变革的必由之路。然而多种能源***的紧密耦合,使得IES(综合能源***)需考虑不同能源***的线路或管道传输速率、不同能源的相互转换及供能效益等方面。这些使综合能源***较于单一***能源,物理结构更复杂,模型构建更困难,从而使综合能源***薄弱分析更困难。因此,对IES薄弱点能快速有效地辨识具有重要的理论和应用研究价值。
现有技术中,对于单一电力***的薄弱辨识主要有潮流熵、复杂网络理论、信息熵等方法。潮流熵结合熵理论,定量分析了***功率分布的不均衡性;复杂理论基于***结构,可构造网络连锁故障模型。这些方法都是基于***物理结构来对***运行状态进行分析的。对于单一热力***和单一天然气***的薄弱辨识主要有敏感性分析、重要度、可靠性跟踪等方法。敏感性分析计算敏感性系数,作为***薄弱指标;可靠性跟踪基于最小割集对***元件进行不可靠性分配,判断元件薄弱性。这些方法都要对管网***各层级进行可靠性计算分析,才能评价***和元件的薄弱性。
本申请发明人在实施本发明的过程中,发现现有技术的方法,至少存在如下技术问题:
对于综合能源***的薄弱辨识,由于***结构复杂且模型构建困难,现有技术中对IES薄弱性分析的研究只处于初步探索阶段。目前可利用潮流数据,计算节点薄弱指标,对含电-气耦合的能源***进行初步研究。该方法表明从数据分析角度对***进行薄弱辨识是可行的,但对于大型或含更多能源耦合的***,该方法潮流模型复杂,仍难以实现快速响应。
由此可知,现有技术中的方法存在模型复杂和响应速度较慢的技术问题。
发明内容
有鉴于此,本发明提供了一种基于统一潮流模型的综合能源***薄弱点辨识方法,用以解决或者至少部分解决现有技术中的方法存在的模型复杂和响应速度较慢的技术问题。
为了解决上述技术问题,本发明提供了一种基于统一潮流模型的综合能源***薄弱点辨识方法,包括:
步骤S1:分别构建电***模型、热***模型以及气***模型,加入耦合元件模型后建立综合能源***统一潮流模型,其中,综合能源***统一潮流模型包括不平衡量的初值,用以表征综合能源***运行状态;
步骤S2:采集综合能源***统一潮流模型中不平衡量初值的历史数据和实时数据,并利用实时分离窗技术构建随机矩阵;
步骤S3:对随机矩阵中的所有数据进行标准化预处理;
步骤S4:计算随机矩阵的M-P率和圆环率,基于M-P率和圆环率对综合能源***进行异常状态识别;
步骤S5:对随机矩阵中与待分析节点对应的数据进行复制,构建增广矩阵,并计算增广矩阵的平均谱半径MSR,MSR用以表征特征值的统计特性,根据平均谱半径的分布情况以及步骤S4中识别的异常状态对综合能源***薄弱点进行辨识。
在一种实施方式中,步骤S1具体包括:
步骤S1.1:分别构建电***模型、热***模型以及气***模型;
步骤S1.2:采用耦合元件将电***模型、热***模型以及气***模型建综合能源***统一潮流模型,综合能源***统一潮流模型的不平衡量ΔF(χ)表示为:
其中,ΔF(χ)为不平衡量;ΔP和ΔQ为电***的有功和无功偏差;ΔΦ、Δph、ΔTs、ΔTr分别为热***的热功率、回路压降、供热温度、回热温度偏差;Δνq和Δpg为气***的管道流量和回路压降偏差;上标sp代表已知的给定值;下标load代表负荷节点;Cs、Cr为与供、回热网络结构和管道流量有关的矩阵,bs、br为与供、回热温度有关的列向量,
综合能源***的状态变量χ为:
利用牛顿拉夫逊方法,求解综合能源***的不平衡量和状态变量:
χ(i+1)=χ(i)-J-1ΔF(χ)
式中:i是迭代数;J是雅可比矩阵。
在一种实施方式中,步骤S2具体包括:
步骤S2.1:将每一个时间点的不平衡量初值构成一个列向量,x(ti)=[x1(ti),x2(ti),...xN(ti)]T;
步骤S2.2:采集不平衡量初值列向量的历史数据和实时数据;
步骤S2.3:在历史数据的基础上进行实时分析的数据分割,通过采集第i时刻和第i时刻之前长度为T-1的数据后构建随机矩阵XN×T:
XN×T=[x(ti-T+1)...x(ti)]
其中,T表示窗口宽度。
在一种实施方式中,步骤S3具体为采用下述方式对数据进行标准化预处理:
在一种实施方式中,步骤S4具体包括:
步骤S4.3:对于M-P率,根据高斯核函数与M-P率对比情况,判定***是否存在异常状态,对于圆环率,根据特征值的分布情况判定***是否存在异常状态。
在一种实施方式中,步骤S5具体包括:
步骤S5.2:计算增广矩阵的平均谱半径MSR,MSR用以表征特征值的统计特性,具体为复平面上所有特征值距离中心点的平均值rMSR:
其中,λi(i=1,2,...N)为矩阵特征值;
步骤S5.3:采用熵理论,将MSR的分布情况数值化,熵用于度量一个***的无序程度,定义为:
其中,C为常数;l为状态数;p(ωi)(i=1,...,l)为状态出现概率;
结合MSR的分布情况和步骤S4中识别的异常状态,计算每个节点的评价值S,获得薄弱点辨识结果;
其中,M为权重因子;C为常数,S值表征了节点对***的影响程度,S值越大,则节点薄弱程度越高。
本申请实施例中的上述一个或多个技术方案,至少具有如下一种或多种技术效果:
本发明提供了一种基于统一潮流模型的综合能源***薄弱点辨识方法,首先,分别构建电***模型、热***模型以及气***模型,加入耦合元件模型后建立综合能源***统一潮流模型;然后,采集综合能源***统一潮流模型中不平衡量初值的历史数据和实时数据,并利用实时分离窗技术构建随机矩阵;接着,对随机矩阵中的所有数据进行标准化预处理;接下来计算随机矩阵的M-P率和圆环率,基于M-P率和圆环率对综合能源***进行异常状态识别;最后对随机矩阵中与待分析节点对应的数据进行复制,构建增广矩阵,并计算增广矩阵的平均谱半径MSR,MSR用以表征特征值的统计特性,根据平均谱半径的分布情况以及步骤S4中识别的异常状态对综合能源***薄弱点进行辨识。
由于本发明将随机矩阵理论引入综合能源***薄弱辨识中,以数据驱动的方式分析***运行状态,不需要对综合能源***物理结构进行辨识,避免了对综合能源***的建模过程,无需考虑***物理结构,模型构建简单,可实现对复杂***薄弱点辨识的快速响应。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明提供的一种基于统一潮流模型的综合能源***薄弱点辨识方法的整体流程示意图;
图2为本发明实施例中综合能源***示意图;
图3为本发明实施例中综合能源***算例模型图;
图4为本发明实施例中***初始状态的M-P率;
图5为本发明实施例中***初始状态的圆环率;
图6为本发明实施例中***负荷开始逐渐增加时的M-P率;
图7为本发明实施例中***负荷开始逐渐增加时的圆环率;
图8为本发明实施例中***接近崩溃的M-P率;
图9为本发明实施例中***接近崩溃的圆环率;
图10为本发明实施例中热节点1的MSR;
图11为本发明实施例中热节点9的MSR;
图12为本发明实施例中电节点10的圆环率;
图13为本发明实施例中电节点4的圆环率;
图14为本发明实施例中气节点5的圆环率;
图15为本发明实施例中热节点9的圆环率;
具体实施方式
本发明的目的是提出一种基于统一潮流模型和随机矩阵的电-热-气综合能源***薄弱点辨识方法,以数据驱动角度,对***无需进行物理结构辨识,同时也避免了***建模过程,为综合能源***薄弱辨识提供了新思路,可以简化建模的过程,并提高***状态识别和薄弱点辨识的准确性和快速性。
为达到上述目的,本发明的主要构思如下:
随机矩阵理论是一种具有高度普适性的大数据分析方法。本发明从数据驱动的角度,将随机矩阵理论引入到综合能源***中,提出一种基于统一潮流模型的综合能源***薄弱点辨识方法,以改进的综合能源***算例仿真,验证了所提方法的有效性。该方法无需考虑***物理结构,模型构建简单,可实现对复杂***薄弱点辨识的快速响应。
本发明的方法具体包括:(1)建立综合能源***统一潮流模型;(2)选取潮流数据并建立随机矩阵;(3)标准化预处理随机矩阵;(4)识别***异常状态;(5)辨识***薄弱点;本发明将随机矩阵理论引入综合能源***薄弱辨识中,以数据驱动的方式分析***运行状态,不需要对综合能源***物理结构进行辨识,避免了对综合能源***的建模过程,为综合能源***薄弱辨识提供了新思路。
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本实施例提供了一种基于统一潮流模型的综合能源***薄弱点辨识方法,请参见图1,该方法包括:
步骤S1:分别构建电***模型、热***模型以及气***模型,加入耦合元件模型后建立综合能源***统一潮流模型,其中,综合能源***统一潮流模型包括不平衡量的初值,用以表征综合能源***运行状态。
在一种实施方式中,步骤S1具体包括:
步骤S1.1:分别构建电***模型、热***模型以及气***模型;
步骤S1.2:采用耦合元件将电***模型、热***模型以及气***模型建综合能源***统一潮流模型,综合能源***统一潮流模型的不平衡量ΔF(χ)表示为:
其中,ΔF(χ)为不平衡量;ΔP和ΔQ为电***的有功和无功偏差;ΔΦ、Δph、ΔTs、ΔTr分别为热***的热功率、回路压降、供热温度、回热温度偏差;Δνq和Δpg为气***的管道流量和回路压降偏差;上标sp代表已知的给定值;下标load代表负荷节点;Cs、Cr为与供、回热网络结构和管道流量有关的矩阵,bs、br为与供、回热温度有关的列向量,
综合能源***的状态变量χ为:
利用牛顿拉夫逊方法,求解综合能源***的不平衡量和状态变量:
χ(i+1)=χ(i)-J-1ΔF(χ)
式中:i是迭代数;J是雅可比矩阵。
在具体的实施过程中,各个模型的构建方式如下:
1、电力***模型的构建
传统电力***潮流模型,节点的功率表达为:
式中:Ne为节点数;V为节点电压;I为节点电流;Y为节点导纳矩阵;S为复功率;P为有功功率;Q为无功功率。
2、热力***模型的构建
1)水力模型
节点满足流量连续性方程:
各回路满足回路压力方程:
Bhhf=0
式中:Bh为回路-支路关联矩阵;hf为压头损失向量。
各管道流量和压头损失之间的关系为:
式中:Kh为管道阻力系数矩阵。
2)热力模型
热力模型用于确定每个节点的温度。每个热负荷节点有:供热温度Ts,输出温度To,回热温度Tr。各节点的热功率Φ为:
式中:Cp为水的比热。采用温降方程表示管道始末端温度的关系:
式中:Lh为管道长度;Tstart和Tend为管道起始和末端节点的温度;Ta为环境温度;λh为管道热传导系数。热水在节点处混合前后的关系为:
3、天然气***模型
天然气模型与热力网络水力模型相似,满足流量连续性方程:
Agνg=νq
式中:Ag为节点-支路关联矩阵;νg为各管道气体流量;νq为各节点流出气体流量。回路压力方程表示为:
BgΔp=0
式中:Bg为回路-支路关联矩阵;Δp为管道压降向量。气体流量和管道压降之间的关系为:
式中:Kg为管道阻力系数,指数k与气网压力水平有关。
具体来说,本发明构建的模型包括电力***潮流模型、热***中的水力和热力潮流模型、天然气***潮流模型。
耦合元件是实现不同能源之间的转换的设备,包括热电联产、燃气锅炉、电锅炉等,其中热电联产(combined heat andpower,CHP)可利用燃气的高品位能量发电,低品位能量供热供冷,用能效率可达80%以上,是现今综合能源***中最具商业前景的一种运营模式,本发明建立含CHP机组的耦合模型;
式中:PCHP和ΦCHP为CHP机组输出的电功率和热功率;Fin为耗气量;cm为CHP机组的热电比例系数;ηe为CHP机组的发电效率。
步骤S2:采集综合能源***统一潮流模型中不平衡量初值的历史数据和实时数据,并利用实时分离窗技术构建随机矩阵。
在一种实施方式中,步骤S2具体包括:
步骤S2.1:将每一个时间点的不平衡量初值构成一个列向量,x(ti)=[x1(ti),x2(ti),...xN(ti)]T;
步骤S2.2:采集不平衡量初值列向量的历史数据和实时数据;
步骤S2.3:在历史数据的基础上进行实时分析的数据分割,通过采集第i时刻和第i时刻之前长度为T-1的数据后构建随机矩阵XN×T:
XN×T=[x(ti-T+1)...x(ti)]
其中,T表示窗口宽度。
具体来说,把获取的历史数据和实时数据保留下来,添加高斯白噪声,用来模拟实际运行获得的数据构建矩阵X:
X=[x(t1),x(t2),...x(ti)...]
然后采用实时分离窗技术进行数据分割,采集第i时刻和第i时刻之前长度为T-1的数据,共同构建随机矩阵。
步骤S3:对随机矩阵中的所有数据进行标准化预处理。
具体来说,通过对随机矩阵标准化预处理包括对矩阵中所有数据进行去量纲化和数值归一化,可以使得各个指标具有可比性。
在一种实施方式中,步骤S3具体为采用下述方式对数据进行标准化预处理:
步骤S4:计算随机矩阵的M-P率和圆环率,基于M-P率和圆环率对综合能源***进行异常状态识别。
在一种实施方式中,步骤S4具体包括:
步骤S4.3:对于M-P率,根据高斯核函数与M-P率对比情况,判定***是否存在异常状态,对于圆环率,根据特征值的分布情况判定***是否存在异常状态。
具体来说,步骤S4.1中,对于N×T阶非Hermitian矩阵X,当N,T→∞且比例系数c=N/T∈(0,1]时,该矩阵的样本协方差矩阵当矩阵X中的元素是独立同分布且为随机变量,并满足均值μ=0、方差σ2<∞时,相应协方差矩阵的经验谱分布收敛于M-P率。
由于不平衡量的初值ΔF(0)(χ)数据为实数,需要将实数域的特征值映射到复数域,即进行奇异化处理。本发明采用酉矩阵U求取的奇异值等价矩阵U为满足Haar分布的酉阵,且可知对于L个该矩阵,考虑其乘积并进行标准化处理:
在具体实施过程中,步骤S4.3中,对于M-P率,本实施例采用高斯核函数与M-P率对比,其中,高斯核函数是对特征值分布进行非参数估计:
式中:fn(x)为非参数函数;h为宽带参数;K(.)为核函数。
如果二者曲线吻合较好,则***稳定,若二者相差较大,则***存在异常情况。高斯核函数可以将数据映射到无穷维,是某种沿径向对称的标量函数。
对于圆环率,将特征值画在复平面上,若所有特征值都在圆环内,则***稳定;若仅有极少特征值(小于2.5%)落在圆环外,则认为是统计误差。相反,当出现预设数量的特征值落在内环时,可判断***出现异常。
步骤S5:对随机矩阵中与待分析节点对应的数据进行复制,构建增广矩阵,并计算增广矩阵的平均谱半径MSR,MSR用以表征特征值的统计特性,根据平均谱半径的分布情况以及步骤S4中识别的异常状态对综合能源***薄弱点进行辨识。
具体来说,步骤S5是利用随机矩阵中的MSR、增广矩阵,并结合熵理论对综合能源***薄弱点有效辨识。
在一种实施方式中,步骤S5具体包括:
步骤S5.2:计算增广矩阵的平均谱半径MSR,MSR用以表征特征值的统计特性,具体为复平面上所有特征值距离中心点的平均值rMSR:
其中,λi(i=1,2,...N)为矩阵特征值;
步骤S5.3:采用熵理论,将MSR的分布情况数值化,熵用于度量一个***的无序程度,定义为:
其中,C为常数;l为状态数;p(ωi)(i=1,...,l)为状态出现概率;
结合MSR的分布情况和步骤S4中识别的异常状态,计算每个节点的评价值S,获得薄弱点辨识结果;
其中,M为权重因子;C为常数,S值表征了节点对***的影响程度,S值越大,则节点薄弱程度越高。
具体来说,增广矩阵使MSR结果偏向于该节点的数据特征,本发明采用熵理论,将MSR的分布情况数值化。熵用于度量一个***的无序程度。然后选取节点的MSR波动值作为状态量,和分别为正常和异常状态下第i个增广矩阵的MSR值。对ωi进行归一化处理:
其中,μi为第i个节点的MSR概率化数值;nh为矩阵行数。
接下来则计算每个节点的评价值S,由于每个节点会对***产生不同程度的影响,其薄弱性可利用MSR计算权重因子:
S值表征了节点对***的影响程度,S值越大则MSR值分布越无序,节点状态变化对***状态影响越大,节点薄弱程度越高。
下面结合具体实例和附图,对本发明的技术方案作进一步具体的说明。
如附图1所示,一种基于统一潮流模型和随机矩阵的电-热-气综合能源***薄弱点辨识方法,包括以下步骤:
步骤1,单独构建电***13节点潮流模型、热***13节点潮流模型、气***7节点潮流模型,加入2个CHP机组作为耦合元件的耦合元件模型,共同建立33节点综合能源***的统一潮流模型。附图2为综合能源***示意图,附图3为综合能源***算例模型图;
步骤2,采集综合能源***统一潮流模型中不平衡量初值的数据。通过不断改变***负荷大小来改变***状态,使得不平衡量初值不断被刷新,累计历史数据和实时数据。这些数据以时间为顺序排序,加上综合能源***中节点特征具有空间特性,构成具有时空特性的数据源。利用实时分离窗技术将数据源分解为一个个随机矩阵;
步骤3,对随机矩阵中的所有数据进行标准化预处理,即去量纲化和数值归一化;
步骤4,通过改变***负荷使***状态发生变化,计算M-P率和圆环率,分析M-P率曲线和圆环率特征值分布随***状态变化而变化的过程。
附图4是***初始状态的M-P率,附图5是***初始状态的圆环率,附图6是***负荷开始逐渐增加时的M-P率,附图7是***负荷开始逐渐增加时的圆环率。附图8是***接近崩溃的M-P率,附图9是***接近崩溃的圆环率。
对于M-P率,高斯核函数与其对比。若二者吻合较好,则***稳定,若二者相差较大,则***存在异常情况。对于圆环率,若所有特征值都在圆环内,则***稳定;当出现数量可观的特征值落在内环时,可判断***出现异常;
步骤5,改变某一节点的负荷大小,分别构建以每一节点为重点观察对象的增广矩阵,计算增广矩阵的MSR,分析某一节点变化对***的影响。附图10是热节点1的MSR,附图11是热节点9的MSR;
结合熵理论计算综合能源***中所有节点的薄弱值并进行排序。本发明通过圆环率对该薄弱点排序进行了验证,即改变某一节点的负荷大小,观察***同一时刻下的圆环率,验证表明了该方法的有效性。附图12是电节点10的圆环率,附图13是电节点4的圆环率,附图14是气节点5的圆环率,附图15是热节点9的圆环率;
本发明从数据驱动的角度,将随机矩阵理论引入到综合能源***中,提出一种基于统一潮流模型的综合能源***薄弱点辨识方法,以改进的综合能源***算例仿真,验证了所提方法的有效性。该方法无需考虑***物理结构,模型构建简单,可实现对复杂***薄弱点辨识的快速响应。
尽管已描述了本发明的优选实施例,但本领域内的技术人员一旦得知了基本创造性概念,则可对这些实施例做出另外的变更和修改。所以,所附权利要求意欲解释为包括优选实施例以及落入本发明范围的所有变更和修改。
显然,本领域的技术人员可以对本发明实施例进行各种改动和变型而不脱离本发明实施例的精神和范围。这样,倘若本发明实施例的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。
Claims (4)
1.一种基于统一潮流模型的综合能源***薄弱点辨识方法,其特征在于,包括:
步骤S1:分别构建电***模型、热***模型以及气***模型,加入耦合元件模型后建立综合能源***统一潮流模型,其中,综合能源***统一潮流模型包括不平衡量的初值,用以表征综合能源***运行状态;
步骤S2:采集综合能源***统一潮流模型中不平衡量初值的历史数据和实时数据,并利用实时分离窗技术构建随机矩阵;
步骤S3:对随机矩阵中的所有数据进行标准化预处理;
步骤S4:计算随机矩阵的M-P率和圆环率,基于M-P率和圆环率对综合能源***进行异常状态识别;
步骤S5:对随机矩阵中与待分析节点对应的数据进行复制,构建增广矩阵,并计算增广矩阵的平均谱半径MSR,MSR用以表征特征值的统计特性,根据平均谱半径的分布情况以及步骤S4中识别的异常状态对综合能源***薄弱点进行辨识;
其中,步骤S1具体包括:
步骤S1.1:分别构建电***模型、热***模型以及气***模型;
步骤S1.2:采用耦合元件将电***模型、热***模型以及气***模型建综合能源***统一潮流模型,综合能源***统一潮流模型的不平衡量ΔF(χ)表示为:
其中,ΔF(χ)为不平衡量;ΔP和ΔQ为电***的有功和无功偏差;ΔΦ、Δph、ΔTs、ΔTr分别为热***的热功率、回路压降、供热温度、回热温度偏差;Δνq和Δpg为气***的管道流量和回路压降偏差;上标sp代表已知的给定值;下标load代表负荷节点;Cs、Cr为与供、回热网络结构和管道流量有关的矩阵,bs、br为与供、回热温度有关的列向量,
综合能源***的状态变量χ为:
利用牛顿拉夫逊方法,求解综合能源***的不平衡量和状态变量:
χ(i+1)=χ(i)-J-1ΔF(χ)
式中:i是迭代数;J是雅可比矩阵;
步骤S4具体包括:
步骤S4.3:对于M-P率,根据高斯核函数与M-P率对比情况,判定***是否存在异常状态,对于圆环率,根据特征值的分布情况判定是否存在异常状态。
2.如权利要求1所述的方法,其特征在于,步骤S2具体包括:
步骤S2.1:将每一个时间点的不平衡量初值构成一个列向量,x(ti)=[x1(ti),x2(ti),…xN(ti)]T;
步骤S2.2:采集不平衡量初值列向量的历史数据和实时数据;
步骤S2.3:在历史数据的基础上进行实时分析的数据分割,通过采集第i时刻和第i时刻之前长度为T-1的数据后构建随机矩阵XN×T:
XN×T=[x(ti-T+1)...x(ti)]
其中,T表示窗口宽度。
4.如权利要求1所述的方法,其特征在于,步骤S5具体包括:
步骤S5.2:计算增广矩阵的平均谱半径MSR,MSR用以表征特征值的统计特性,具体为复平面上所有特征值距离中心点的平均值rMSR:
其中,λi(i=1,2,...N)为矩阵特征值;
步骤S5.3:采用熵理论,将MSR的分布情况数值化,熵用于度量一个***的无序程度,定义为:
其中,C为常数;l为状态数;p(ωi)(i=1,...,l)为状态出现概率;
结合MSR的分布情况和步骤S4中识别的异常状态,计算每个节点的评价值S,获得薄弱点辨识结果;
其中,M为权重因子;C为常数,S值表征了节点对***的影响程度,S值越大,则节点薄弱程度越高。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910876207.8A CN110601185B (zh) | 2019-09-17 | 2019-09-17 | 一种基于统一潮流模型的综合能源***薄弱点辨识方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910876207.8A CN110601185B (zh) | 2019-09-17 | 2019-09-17 | 一种基于统一潮流模型的综合能源***薄弱点辨识方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110601185A CN110601185A (zh) | 2019-12-20 |
CN110601185B true CN110601185B (zh) | 2021-07-06 |
Family
ID=68860268
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910876207.8A Expired - Fee Related CN110601185B (zh) | 2019-09-17 | 2019-09-17 | 一种基于统一潮流模型的综合能源***薄弱点辨识方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110601185B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111509707A (zh) * | 2020-04-24 | 2020-08-07 | 贵州电网有限责任公司 | 一种能源互联网热力网络潮流标幺值模型建模方法 |
CN111400937B (zh) * | 2020-04-29 | 2023-12-01 | 上海电力大学 | 一种综合能源***潮流计算方法 |
CN113255105B (zh) * | 2021-04-26 | 2022-10-18 | 上海电力大学 | 双向耦合网络结构的电气热综合能源***潮流计算方法 |
CN113609651B (zh) * | 2021-07-13 | 2024-07-19 | 华南理工大学 | 基于统一潮流模型的电-气互联***作用机理分析方法 |
CN113537844B (zh) * | 2021-09-15 | 2021-12-17 | 山东大学 | 基于随机矩阵的区域能源互联网负荷行为分析方法及*** |
CN114154855A (zh) * | 2021-12-02 | 2022-03-08 | 深圳供电局有限公司 | 一种区域配电网状态分析方法 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9639897B2 (en) * | 2013-11-12 | 2017-05-02 | Bank Of America Corporation | Predicting economic conditions |
CN106992551B (zh) * | 2017-06-05 | 2019-06-18 | 合肥工业大学 | 基于模糊c均值和差分进化混合算法的光伏逆变控制器参数辨识方法 |
CN107947151B (zh) * | 2017-11-02 | 2021-01-08 | 广西电网有限责任公司电力科学研究院 | 一种电力***关键节点的识别方法 |
CN108039726B (zh) * | 2017-12-12 | 2021-02-09 | 国网山东省电力公司德州供电公司 | 一种基于多智能体***的能源局域网分布式协同控制方法 |
CN108416507B (zh) * | 2018-02-07 | 2021-08-20 | 天津大学 | 面向电-气耦合综合能源***的静态灵敏度分析方法 |
CN108764603A (zh) * | 2018-03-30 | 2018-11-06 | 广东电网有限责任公司 | 一种基于高维随机矩阵的用户窃电行为辨识方法 |
CN109347095A (zh) * | 2018-10-29 | 2019-02-15 | 国网山东省电力公司枣庄供电公司 | 考虑多能互补特性的主动配电网分层分布式协调控制器 |
CN109274097A (zh) * | 2018-11-16 | 2019-01-25 | 四川大学 | 一种基于随机矩阵理论的电力***暂态稳定态势评估方法 |
CN109800906A (zh) * | 2018-12-25 | 2019-05-24 | 天津大学 | 面向新能源消纳的配水网与配电网联合经济调度方法 |
-
2019
- 2019-09-17 CN CN201910876207.8A patent/CN110601185B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN110601185A (zh) | 2019-12-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110601185B (zh) | 一种基于统一潮流模型的综合能源***薄弱点辨识方法 | |
CN110441065B (zh) | 基于lstm的燃气轮机在线检测方法与装置 | |
Wang et al. | Wind power curve modeling and wind power forecasting with inconsistent data | |
Xue et al. | Multi-step ahead forecasting of heat load in district heating systems using machine learning algorithms | |
CN109992921B (zh) | 一种燃煤电厂锅炉热效率的在线软测量方法及*** | |
Jung et al. | Prediction of building energy consumption using an improved real coded genetic algorithm based least squares support vector machine approach | |
CN106529719B (zh) | 基于粒子群优化算法风速融合的风电功率预测的方法 | |
Huang et al. | Feature selection and hyper parameters optimization for short-term wind power forecast | |
Sedano et al. | A soft computing method for detecting lifetime building thermal insulation failures | |
CN112101480A (zh) | 一种多变量聚类与融合的时间序列组合预测方法 | |
Elfaki et al. | Prediction of electrical output power of combined cycle power plant using regression ANN model | |
Wang et al. | Wind power curve modeling with asymmetric error distribution | |
Gao et al. | Comprehensive study on sensitive parameters for chiller fault diagnosis | |
CN109272258B (zh) | 基于K-means聚类的区域风光发电资源评估方法 | |
Zhang et al. | Energy modeling using an effective latent variable based functional link learning machine | |
Effendy et al. | The prediction of the oxygen content of the flue gas in a gas-fired boiler system using neural networks and random forest | |
CN113762591A (zh) | 一种基于gru和多核svm对抗学习的短期电量预测方法及*** | |
Peng et al. | An improved 2DCNN with focal loss function for blade icing detection of wind turbines under imbalanced SCADA data | |
CN110414734B (zh) | 一种计及风资源利用率预测评估的方法 | |
Tian et al. | Causal network construction based on convergent cross mapping (CCM) for alarm system root cause tracing of nonlinear industrial process | |
Lin et al. | The study on classification and prediction for data mining | |
CN116595845A (zh) | 一种螺旋翅片管式换热器的结构优化方法 | |
CN106816871B (zh) | 一种电力***状态相似性分析方法 | |
Gu et al. | A modified clustering procedure for energy consumption monitoring in the steam turbine considering volume effect | |
Hua et al. | A novel sampled-data asynchronous CatBoost model with knowledge-based guidance for product quality estimation in the industrial processes |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210706 |
|
CF01 | Termination of patent right due to non-payment of annual fee |