CN115800269B - 配电网量测数据驱动的电压功率灵敏度拓扑约束估计方法 - Google Patents
配电网量测数据驱动的电压功率灵敏度拓扑约束估计方法 Download PDFInfo
- Publication number
- CN115800269B CN115800269B CN202211695968.1A CN202211695968A CN115800269B CN 115800269 B CN115800269 B CN 115800269B CN 202211695968 A CN202211695968 A CN 202211695968A CN 115800269 B CN115800269 B CN 115800269B
- Authority
- CN
- China
- Prior art keywords
- voltage
- power sensitivity
- sensitivity
- node
- voltage power
- 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
- 230000035945 sensitivity Effects 0.000 title claims abstract description 73
- 238000000034 method Methods 0.000 title claims abstract description 57
- 238000005259 measurement Methods 0.000 title claims abstract description 21
- 238000013178 mathematical model Methods 0.000 claims abstract description 10
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 3
- 239000011159 matrix material Substances 0.000 claims description 35
- 238000002347 injection Methods 0.000 claims description 10
- 239000007924 injection Substances 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 9
- 238000005457 optimization Methods 0.000 claims description 5
- 238000005070 sampling Methods 0.000 claims description 3
- 230000009467 reduction Effects 0.000 abstract description 3
- 239000000243 solution Substances 0.000 description 7
- 230000000694 effects Effects 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 238000007792 addition Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000007429 general method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
Images
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明涉及电压功率灵敏度估计方法领域,具体涉及一种配电网量测数据驱动的电压功率灵敏度拓扑约束估计方法,包括如下步骤:对配电网电压功率灵敏度进行成分分解,将电压功率灵敏度分解为主成分和次成分;采集配电网节点量测数据和网络拓扑信息;建立电压功率灵敏度主成分估计的数据驱动的回归数学模型,并进行求解;建立电压功率灵敏度次成分估计的数据驱动的回归数学模型,并进行求解;将回归数学模型所得的主成分与次成分叠加得到最终的电压功率灵敏度。本发明通过主成分的变量降维和次成分的变量大小约束,显著降低了多重共线性问题对数据驱动方法的负面影响,相比现有方法精度大幅提升。
Description
技术领域
本发明涉及电压功率灵敏度估计方法技术领域,具体涉及配电网量测数据驱动的电压功率灵敏度拓扑约束估计方法。
背景技术
电压功率灵敏度的获得分为基于物理模型和基于数据驱动的两种方法:物理模型法需已知配电网络的详细参数并构建潮流模型,一般基于牛顿法潮流的雅克比矩阵求逆得到;数据驱动法无需配电网络的物理模型和具体参数,由节点电压、节点注入功率的量测数据通过求解回归数学模型得到。在实际工程中,配电网各元件的准确参数难以获得,且会随运行工况的变化而改变,因此,物理模型法存在一定的局限。随着智能电表、微型同步相量测量单元的快速发展,可采集配电网各节点电压和注入功率的高密度时序数据。因此,基于数据驱动的电压功率灵敏度估计方法近年来得到国内外学者的广泛关注。由于实际配电网采集的数据存在突出的多重共线性问题,即不同节点的注入功率具有显著的相关性,将导致电压功率灵敏度估计的精度急剧下降,是数据驱动方法工程应用的主要障碍。
现有数据驱动的配电网电压功率灵敏度估计方法有最小二乘回归、岭回归、Lasso回归等。岭回归、Lasso回归等通过在最小二乘回归模型中添加正则项的方法为现有的主流方法,岭回归、Lasso回归通过在最小二乘回归模型的基础上添加正则项,相比最小二乘回归能够有效提升多重共线性问题下的精度,是一种处理多重共线性问题的通用方法,但其并未考虑配电网的独特特征。由于配电网电压功率灵敏度估计的回归变量维度很高、多重共线性问题非常突出,致使添加正则项的方法获得的结果仍达不到高精度要求,目前仍难以应用于实际工程。
发明内容
为了克服上述的技术问题,本发明的目的在于提供配电网量测数据驱动的电压功率灵敏度拓扑约束估计方法,充分利用可获得的配电网信息,为有效克服多重共线性问题的负面影响,有针对的提出嵌入配电网拓扑约束的估计方法,具体包括了主、次成分分解,建立主、次成分估计的回归数学模型并依次求解等步骤,从而实现配电网电压功率灵敏度估计精度的大幅度提升。
本发明的目的可以通过以下技术方案实现:
配电网量测数据驱动的电压功率灵敏度拓扑约束估计方法,包括如下步骤:
步骤1:对配电网电压功率灵敏度进行成分分解,基于线性DistFlow模型将电压功率灵敏度分解为主成分和次成分;
步骤2:采集配电网节点量测数据和网络拓扑信息;
步骤3:建立电压功率灵敏度主成分估计的数据驱动的回归数学模型,并进行求解;
步骤4:建立电压功率灵敏度次成分估计的数据驱动的回归数学模型,并进行求解;
步骤5:将所述步骤3和步骤4所得的主成分与次成分叠加得到最终的电压功率灵敏度。
所述步骤1具体为:基于配电网线性DistFlow模型得到节点电压幅值平方与节点注入功率之间的线性关系,将电压有功功率灵敏度、电压无功功率灵敏度近似为下式:
其中,F为降阶支路-节点关联矩阵的逆矩阵,Dr、Dx分别为支路电阻、支路电抗的对角矩阵,FDrFT、FDxFT定义为电压有功功率灵敏度、电压无功功率灵敏度的主成分,与精确值之间的误差定义为次成分;
进一步的,所述步骤2具体为:由配电网量测装置采集获得节点电压Vi,t、节点有功功率Pk,t、节点无功功率Qk,t数据,经差分计算得到各时刻的增量ΔVi,t、ΔPk,t、ΔQk,t,其中i、k表示节点编号,t表示时刻;由配电管理系统获得网络拓扑信息,形成支路-节点关联矩阵M,将配电网根节点所在的列删除形成降阶支路-节点关联矩阵Mr,对Mr求逆得到矩阵F。
式中,Rj、Xj分别为对角矩阵Dr、Dx的第j个元素,为决策变量;分别为矩阵中的元素;Fij、Fkj为矩阵F中的元素,由降阶支路-节点关联矩阵求逆得到,即所嵌入的网络拓扑信息;ΔVi,t、ΔPk,t、ΔQk,t分别为t时刻节点i的电压幅值、节点k的注入有功、无功功率相对t-1时刻的增量,由量测数据经差分计算得到;m为数据采样的时刻数量;n为非根节点的数量;λ为线路电抗与电阻的比值上限;式(3a)为目标函数,式(3b)-(3c)为主成分中各元素基于拓扑信息的约束条件;式(3d)-(3f)为配电线路电抗、电阻以及其比值的约束条件,以描述配电线路电抗与电阻之比在一定范围内;采用数学优化求解器求解式(3a)-(3f)所示的二次规划问题得到电压功率灵敏度主成分
式中,分别为矩阵中的元素,为决策变量;由步骤3求解得到,在本问题中为常量;ξ为小于1且大于0的系数;其他符号含义同式(3a)-(3f);式(4a)为目标函数,式(4b)-(4c)为次成分相比主成分为微小量的约束条件;采用数学优化求解器求解式(4a)-(4c)所示的二次规划问题得到电压功率灵敏度次成分
本发明的有益效果:(1)通过主成分的变量降维和次成分的变量大小约束,显著降低了多重共线性问题对数据驱动方法的负面影响,相比现有方法精度大幅提升。(2)参数的物理意义明确、敏感性低,易于设定、无需反复调参,且计算耗时与现有方法相近,具有良好的工程应用前景。
附图说明
图1是本发明方法流程示意图。
具体实施方式
下面将结合本发明实施例,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
请参阅图1,本发明配电网量测数据驱动的电压功率灵敏度拓扑约束估计方法,包括如下步骤:
步骤1:配电网电压功率灵敏度主、次成分分解,详细过程如下:
基于配电网线性DistFlow模型可得到节点电压幅值平方与节点注入功率之间的线性关系,该模型忽略了网络损耗。进一步,考虑电压幅值V在1.0pu附近,电压有功功率灵敏度、电压无功功率灵敏度可近似为式(1)。
其中,F为降阶支路-节点关联矩阵的逆矩阵,Dr、Dx分别为支路电阻、支路电抗的对角矩阵。FDrFT、FDxFT与精确的电压有功功率灵敏度、电压无功功率灵敏度之间存在误差,但相比精确值为微小量。因此,FDrFT、FDxFT定义为电压有功功率灵敏度、电压无功功率灵敏度的主成分,与精确值之间的误差定义为次成分。
步骤2:采集配电网节点量测数据和网络拓扑信息,详细过程如下:由配电网量测装置采集获得节点电压Vi,t、节点有功功率Pk,t、节点无功功率Qk,t数据,经差分计算得到各时刻的增量ΔVi,t、ΔPk,t、ΔQk,t,其中i、k表示节点编号,t表示时刻;由配电管理系统获得网络拓扑信息,形成支路-节点关联矩阵M,将配电网根节点所在的列删除形成降阶支路-节点关联矩阵Mr,对Mr求逆得到矩阵F;
步骤3:建立电压功率灵敏度主成分估计的数据驱动的回归数学模型,并进行求解,详细过程如下:
式中,Rj、Xj分别为对角矩阵Dr、Dx的第j个元素,为决策变量;分别为矩阵中的元素;Fij、Fkj为矩阵F中的元素,由降阶支路-节点关联矩阵求逆得到,即所嵌入的网络拓扑信息;ΔVi,t、ΔPk,t、ΔQk,t分别为t时刻节点i的电压幅值、节点k的注入有功、无功功率相对t-1时刻的增量,由量测数据经差分计算得到;m为数据采样的时刻数量;n为非根节点的数量;λ为线路电抗与电阻的比值上限,一般低于4,可由运行经验进一步缩小范围。式(3a)为目标函数,式(3b)-(3c)为主成分中各元素基于拓扑信息的约束条件;式(3d)-(3f)为配电线路电抗、电阻以及其比值的约束条件,以描述配电线路电抗与电阻之比在一定范围内。
步骤4:建立电压功率灵敏度次成分估计的数据驱动的回归数学模型,并进行求解,详细过程如下:
式中,分别为矩阵中的元素,为决策变量;由步骤3求解得到,在本问题中为常量;ξ为小于1且大于0的系数;其他符号含义同式(3a)-(3f)。式(4a)为目标函数,式(4b)-(4c)为次成分相比主成分为微小量的约束条件。
作为本发明进一步的实施例,采用IEEE 33节点系统对本发明方案进行验证,节点1为连接上级电网的根节点,其余32个节点的注入功率由实际系统的微型同步相量测量单元采集获得,共3600组数据,并按均值与IEEE 33节点原始功率数据相等为原则将采集的实际数据按比例缩放,通过潮流计算获得对应的节点电压数据,即模拟IEEE 33节点系统实际运行并采集量测数据的过程。电压、功率基准值分别为12.66kV、1MVA。所采集的注入功率及其时序增量数据皆具有显著的多重共线性问题。根据所述的节点电压、节点注入功率量测数据测试本发明方案的有效性,以基于精确系统参数的雅克比矩阵求逆的实际值作为基准,并与常规最小二乘回归、岭回归、Lasso回归比较。
1、估计误差对比
将本发明方案的参数λ、ξ分别设置为4、0.1;岭回归的正则项(目标函数加权的二范数正则项)参数β经调优后设置为8×10-11;Lasso回归的正则项(一范数不等式约束正则项)参数γ经调优后设置为60。
如表1所示为电压功率灵敏度矩阵估计误差的对比。岭回归和Lasso回归精度接近,Lasso回归的归一化均方误差、最大误差低于岭回归,岭回归和Lasso回归的归一化均方误差分别为常规最小二乘回归的18%、14%。本发明方案的精度相比岭回归和Lasso回归进一步提高了一个数量级,归一化均方误差、最大误差、误差中位数分别仅为常规最小二乘回归的0.5%、8%、5%。综上分析,本发明方案的精度相比现有技术方案具有数量级的提升。
表1电压功率灵敏度矩阵的估计误差对比
2、参数敏感性对比
进一步测试不同参数对本发明方案估计精度的影响,并与岭回归和Lasso回归方法对比,分析参数敏感性。如表2所示是本发明方案在不同参数组合下的误差,当λ取极端值100时,误差较最优参数组合仅增加17%;ξ取极端值0、1时,误差较最优参数组合仅分别增加13%、42%。表2显示ξ在0.1-0.4范围对精度几乎没有影响,根据ξ的物理意义,ξ不应设置过大,否则与次成分的定义不符。表2显示λ取1时误差最小,原因是IEEE 33节点系统中多数线路的电抗与电阻的比值接近1,该参数可由配电线路参数经验值获得。表3所示为岭回归和Lasso回归在不同参数设定下的误差,由于这两种方法的参数没有明确物理意义,其合适范围难以预先确定,需要反复调参。本发明方案中,参数λ是线路电抗与电阻的比值约束,参数ξ是次成分与主成分的比值约束,参数的合适范围易确定,即使设置为极端值时也不会显著降低精度,误差仍远低于最优参数下的岭回归和Lasso回归。
表2本发明方案在不同参数组合下的误差
表3岭回归和Lasso回归在不同参数下的误差
3、耗时对比
表4所示为不同方法的耗时对比,由于本发明方案将估计过程分为主成分和次成分估计两个阶段,第一阶段耗时(0.3s)与常规最小二乘回归一致,第二阶段耗时(0.24s)略低于常规最小二乘回归,总体耗时比Lasso回归略高,耗时数量级与现有方法相同。此外,第一阶段估计所得主成分的精度已显著高于其他方法,即表3中ξ=0时的精度。因此,验证了本发明方案的效率。
表4不同方法的耗时对比
综上,本发明方案具有如下优点:
(1)通过主成分的变量降维和次成分的变量大小约束,显著降低了多重共线性问题对数据驱动方法的负面影响,相比现有方法精度大幅提升。
(2)参数的物理意义明确、敏感性低,易于设定、无需反复调参,且计算耗时与现有方法相近,具有良好的工程应用前景。
以上内容仅仅是对本发明所作的举例和说明,所属本技术领域的技术人员对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,只要不偏离发明或者超越本权利要求书所定义的范围,均应属于本发明的保护范围。
Claims (1)
1.配电网量测数据驱动的电压功率灵敏度拓扑约束估计方法,其特征在于,包括如下步骤:
步骤1:对配电网电压功率灵敏度进行成分分解,基于线性DistFlow模型将电压功率灵敏度分解为主成分和次成分;
步骤2:采集配电网节点量测数据和网络拓扑信息;
步骤3:建立电压功率灵敏度主成分估计的数据驱动的回归数学模型,并进行求解;
步骤4:建立电压功率灵敏度次成分估计的数据驱动的回归数学模型,并进行求解;
步骤5:将所述步骤3和步骤4所得的主成分与次成分叠加得到最终的电压功率灵敏度;
所述步骤1具体为:基于配电网线性DistFlow模型得到节点电压幅值平方与节点注入功率之间的线性关系,将电压有功功率灵敏度、电压无功功率灵敏度近似为下式:
其中,F为降阶支路-节点关联矩阵的逆矩阵,Dr、Dx分别为支路电阻、支路电抗的对角矩阵,FDrFT、FDxFT定义为电压有功功率灵敏度、电压无功功率灵敏度的主成分,与精确值之间的误差定义为次成分;
所述步骤2具体为:
由配电网量测装置采集获得节点电压Vi,t、节点有功功率Pk,t、节点无功功率Qk,t数据,经差分计算得到各时刻的增量ΔVi,t、ΔPk,t、ΔQk,t,其中i、k表示节点编号,t表示时刻;由配电管理系统获得网络拓扑信息,形成支路-节点关联矩阵M,将配电网根节点所在的列删除形成降阶支路-节点关联矩阵Mr,对Mr求逆得到矩阵F;
所述步骤3具体为:
式中,Rj、Xj分别为对角矩阵Dr、Dx的第j个元素,为决策变量;分别为矩阵中的元素;Fij、Fkj为矩阵F中的元素,由降阶支路-节点关联矩阵求逆得到,即所嵌入的网络拓扑信息;ΔVi,t、ΔPk,t、ΔQk,t分别为t时刻节点i的电压幅值、节点k的注入有功、无功功率相对t-1时刻的增量,由量测数据经差分计算得到;m为数据采样的时刻数量;n为非根节点的数量;λ为线路电抗与电阻的比值上限;式(3a)为目标函数,式(3b)-(3c)为主成分中各元素基于拓扑信息的约束条件;式(3d)-(3f)为配电线路电抗、电阻以及其比值的约束条件;采用数学优化求解器求解式(3a)-(3f)所示的二次规划问题得到电压功率灵敏度主成分
所述步骤4具体为:
式中,分别为矩阵中的元素,为决策变量;由步骤3求解得到,为常量;ξ为小于1且大于0的系数;其他符号含义同式(3a)-(3f);式(4a)为目标函数,式(4b)-(4c)为次成分相比主成分为微小量的约束条件;采用数学优化求解器求解式(4a)-(4c)所示的二次规划问题得到电压功率灵敏度次成分
所述步骤5具体为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211695968.1A CN115800269B (zh) | 2022-12-28 | 2022-12-28 | 配电网量测数据驱动的电压功率灵敏度拓扑约束估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211695968.1A CN115800269B (zh) | 2022-12-28 | 2022-12-28 | 配电网量测数据驱动的电压功率灵敏度拓扑约束估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115800269A CN115800269A (zh) | 2023-03-14 |
CN115800269B true CN115800269B (zh) | 2023-05-05 |
Family
ID=85428025
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211695968.1A Active CN115800269B (zh) | 2022-12-28 | 2022-12-28 | 配电网量测数据驱动的电压功率灵敏度拓扑约束估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115800269B (zh) |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107039984B (zh) * | 2017-06-19 | 2020-04-14 | 国网重庆市电力公司电力科学研究院 | 一种无功电压分区方法 |
CN107577870B (zh) * | 2017-09-04 | 2020-10-27 | 天津大学 | 基于同步相量量测的配电网电压功率灵敏度鲁棒估计方法 |
CN114744631A (zh) * | 2022-03-30 | 2022-07-12 | 河北工业大学 | 一种基于非pmu配电网的数据驱动电压估计方法 |
-
2022
- 2022-12-28 CN CN202211695968.1A patent/CN115800269B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN115800269A (zh) | 2023-03-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107104442B (zh) | 计及参数模糊性的含风电场电力系统概率潮流计算方法 | |
CN109063276B (zh) | 适用于长时域风速波动的风电场动态等值建模方法 | |
CN106372339A (zh) | 电力电子化电力系统的多速率仿真方法及装置 | |
CN108054757B (zh) | 一种内嵌无功和电压的n-1闭环安全校核方法 | |
CN104239982A (zh) | 一种基于时间序列和神经网络的建筑节假日能耗预测方法 | |
CN110428185B (zh) | 基于伪量测模型的电-热互联综合能源系统抗差状态估计方法 | |
CN110765703B (zh) | 一种风电场聚合特性建模方法 | |
CN107257130B (zh) | 基于区域量测解耦的低压配电网损耗计算方法 | |
CN115146538A (zh) | 基于消息传递图神经网络的电力系统状态估计方法 | |
CN112001112A (zh) | 基于bp神经网络的线性约束黑箱优化问题的优化方法 | |
CN115800269B (zh) | 配电网量测数据驱动的电压功率灵敏度拓扑约束估计方法 | |
CN107204616B (zh) | 基于自适应稀疏伪谱法的电力系统随机状态估计方法 | |
Gokpinar et al. | Generalization of inclusion probabilities in ranked set sampling | |
CN110707703B (zh) | 基于改进Nataf变换的含高维相关不确定源的概率潮流高效计算方法 | |
CN111597494A (zh) | 一种基于非平稳时间序列分解的统计降尺度方法 | |
CN107918920B (zh) | 多光伏电站的出力关联性分析方法 | |
CN114744631A (zh) | 一种基于非pmu配电网的数据驱动电压估计方法 | |
CN114221334B (zh) | 一种基于图神经网络的快速状态估计方法 | |
CN113783186B (zh) | 一种考虑配电网拓扑结构变化的电压预测方法 | |
CN106251238B (zh) | 风电场建模序列离散化步长的选取及模型误差分析方法 | |
CN114236236A (zh) | 一种基于区间动态状态估计的谐波源定位方法 | |
CN109494747A (zh) | 一种基于交替梯度算法的电网概率潮流计算方法 | |
CN106887844B (zh) | 一种基于稀疏恢复的快速分解配电网三相状态估计方法 | |
CN111797564A (zh) | 一种高维分布式光伏出力的相关性样本获取方法及系统 | |
CN110751290A (zh) | 基于人工智能数据驱动的主要断面预测方法 |
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 |