CN115347553B - 基于线性回归初值估计的配电网线路参数辨识方法 - Google Patents
基于线性回归初值估计的配电网线路参数辨识方法 Download PDFInfo
- Publication number
- CN115347553B CN115347553B CN202210862077.4A CN202210862077A CN115347553B CN 115347553 B CN115347553 B CN 115347553B CN 202210862077 A CN202210862077 A CN 202210862077A CN 115347553 B CN115347553 B CN 115347553B
- Authority
- CN
- China
- Prior art keywords
- branch
- mid
- power
- end node
- 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 39
- 238000012417 linear regression Methods 0.000 title claims abstract description 36
- 238000005457 optimization Methods 0.000 claims abstract description 33
- 238000005259 measurement Methods 0.000 claims description 22
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000012897 Levenberg–Marquardt algorithm Methods 0.000 claims description 4
- 125000004122 cyclic group Chemical group 0.000 claims description 4
- 239000013598 vector 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
- 230000009466 transformation Effects 0.000 claims description 3
- 230000001131 transforming effect Effects 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract description 4
- 238000006467 substitution reaction Methods 0.000 abstract description 4
- 230000002411 adverse Effects 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 5
- 238000010276 construction Methods 0.000 description 4
- 238000012360 testing method Methods 0.000 description 3
- 238000012795 verification Methods 0.000 description 3
- 230000004075 alteration Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000004146 energy storage Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
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
-
- 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
- H02J3/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
- H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
-
- 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
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/10—Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
-
- 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
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Measurement Of Resistance Or Impedance (AREA)
Abstract
本发明公开基于线性回归初值估计的配电网线路参数辨识方法,包括以下步骤:1)获取电网的基础数据;2)建立配变局部变量多元非线性模型;3)建立含初值的最小二乘优化模型,并根据含初值的最小二乘优化模型计算出支路末端节点的阻抗;4)递推计算得到配电网线路参数。本发明考虑了变量替代的线性回归策略,可以避免非线性最小二乘估计初值的不良影响,使参数辨识不仅精度高而且更加稳定。
Description
技术领域
本发明涉及配电网参数辨识领域,具体是基于线性回归初值估计的配电网线路参数辨识方法。
背景技术
分布式电源、电动汽车以及用户侧储能正大规模接入配电网,对配电网规划和运行的安全性和可靠性提出新的挑战。为了应对分布式电源以及电动汽车等接入带来的随机性和波动性,需要对配电网络进行实时的拓扑重构和智能调度,从而实现配电网安全运行,需要进行更加精细化的无功控制,保证电能质量,这需要建立在对配电网拓扑结构和线路参数有较为准确判断的基础上。配电网的运行状态估计和优化,同样对拓扑和参数信息的准确判断提出了要求。因此,精细化的拓扑和线路参数信息是配电网高阶优化的前提条件。
目前许多研究无法辨识线路参数,且需要电压相角信息,但对于我国实际的配电网而言,同步相角测量单元PMU安装有限,对方法的数据量支撑不足。
传统的配电网参数辨识方法主要存在以下缺点:传统的配电网参数辨识方法常通过非线性电压降落方程来辨识两两节点拓扑关系和参数,然后求解采用最小二乘方法进行求解。但是,非线性方程采用最小二乘方法求解,其初值对结果精度影响较大。若初值设置不当,则难以取得良好效果。
发明内容
本发明的目的是提供基于线性回归初值估计的配电网线路参数辨识方法,包括以下步骤:
1)获取电网的基础数据。
所述电网基础数据包括电网拓扑结构、可观测节点数量Nobs、节点的有功功率量测量Pi、节点的无功功率量测量Qi、节点的电压幅值量测Vi。
2)建立配变局部变量多元非线性模型。
建立配变局部变量多元非线性模型的步骤包括:
2.1)建立始端电压V0和待求末端电压V1、待求末端电压V2之间的关系式,即:
式中,P1、P2分别为待求末端节点1和待求末端节点2的有功功率。Q1、Q2分别为待求末端节点1和待求末端节点2的无功功率。R1、R2分别为待求末端节点1和待求末端节点2的等效电阻。X1、X2分别为待求末端节点1和待求末端节点2的等效电抗。
2.2)根据公式(1)建立并联节点对待求末端节点1和待求末端节点2的电气参量准确关系式,即:
式中,参数f1(R1,X1)、参数f1(R1,X1)分别如下所示:
2.3)建立N个时刻下,任意支路末端节点i和支路末端节点j的连接判据,即:
式中,Ri、Rj分别为支路末端节点i和支路末端节点j的等效电阻。Xi、Xj分别为支路末端节点i和支路末端节点j的等效电抗。
3)建立含初值的最小二乘优化模型,并根据含初值的最小二乘优化模型计算出支路末端节点的阻抗。
建立含初值的最小二乘优化模型的步骤包括:
3.1)设参数状态变量x包括电压、功率,以常数量k替换电阻、电抗,即:
式中,x1、x2、x3、x4、x5、x6为状态变量。k1、k2、k3、k4、k5、k6为常数量。
3.2)更新参数fi(Ri,Xi)、参数fi(Rj,Xj),得到:
3.3)将公式(5)变换成线性回归方程,得到:
3.4)建立配变局部变量的多元线性回归模型,即:
式中,y和xi,i=1,2,…,6为T维列向量。T为历史正常断面数。
3.5)根据公式(10),建立含初值的最小二乘优化模型,即:
Constraints=Rj,Xj,Ri,Xi>0
3.6)根据含初值的最小二乘优化模型计算出末端节点的阻抗,步骤包括:
3.6.1)利用配变局部变量多元非线性模型计算出常数量k1、k2、k3、k4、k5、k6的初值。
3.6.2)将常数量k1、k2、k3、k4、k5、k6的初值输入到含初值的最小二乘优化模型中,求解得到支路末端节点i的电阻Ri和电抗Xi、支路末端节点j的电阻Ri和电抗Xi。
求解含初值的最小二乘优化模型的方法包括Levenberg-Marquardt算法(列文伯格-马夸尔特法)。
4)递推计算得到配电网线路参数。
所述配电网线路参数包括馈线末端支路Zend的状态变量Umid+1,Pmid+1,Qmid+1,Uend,Pend,Qend、馈线中段支路Zmid的状态变量Umid,Umid+1,Pmid+1,Qmid+1。为馈线末端支路Yend首端的有功功率和无功功率;Pend,Qend为馈线末端支路Yend末端的有功功率和无功功率。
递推计算得到配电网线路参数的步骤包括:
4.1)量测馈线末端支路Zend的功率量测量Umid+1,Pmid+1,Qmid+1,Uend,Pend,Qend,并根据配变局部变量多元非线性模型和含初值的最小二乘优化模型计算得到馈线末端支路Zend的电阻Rend和电抗Xend。其中,Umid+1、Uend为电压。Pmid+1、Pend为有功功率。Qmid+1、Qend为无功功率。
4.2)计算馈线末端支路Zend首端的有功功率和无功功率/>即:
4.3)从馈线末端支路Zend开始,自下而上逐级进行参数估计、循环变量递推计算得到全网的拓扑及支路参数信息。
本发明的技术效果是毋庸置疑的,本发明方考虑了变量替代的线性回归策略,可以避免非线性最小二乘估计初值的不良影响,使参数辨识不仅精度高而且更加稳定。
附图说明
图1为线路电压降相量图;
图2为配变局部变量多元线性回归模型的递推计算图;
图3为方法流程图;
图4为IEEE-33节点系统图;
图5为33rds节点系统图;
图6为69rd节点系统图;
图中1-69均表示节点。
具体实施方式
下面结合实施例对本发明作进一步说明,但不应该理解为本发明上述主题范围仅限于下述实施例。在不脱离本发明上述技术思想的情况下,根据本领域普通技术知识和惯用手段,做出各种替换和变更,均应包括在本发明的保护范围内。
实施例1:
参见图1至图3,基于线性回归初值估计的配电网线路参数辨识方法,包括以下步骤:
1)获取电网的基础数据。
所述电网基础数据包括电网拓扑结构、可观测节点数量Nobs、节点的有功功率量测量Pi、节点的无功功率量测量Qi、节点的电压幅值量测Vi。
2)建立配变局部变量多元非线性模型。
建立配变局部变量多元非线性模型的步骤包括:
2.1)建立始端电压V0和待求末端电压V1、待求末端电压V2之间的关系式,即:
式中,P1、P2分别为待求末端节点1和待求末端节点2的有功功率。Q1、Q2分别为待求末端节点1和待求末端节点2的无功功率。R1、R2分别为待求末端节点1和待求末端节点2的等效电阻。X1、X2分别为待求末端节点1和待求末端节点2的等效电抗。
2.2)根据公式(1)建立并联节点对待求末端节点1和待求末端节点2的电气参量准确关系式,即:
式中,参数f1(R1,X1)、参数f1(R1,X1)分别如下所示:
2.3)建立N个时刻下,任意支路末端节点i和支路末端节点j的连接判据,即:
式中,Ri、Rj分别为支路末端节点i和支路末端节点j的等效电阻。Xi、Xj分别为支路末端节点i和支路末端节点j的等效电抗。
3)建立含初值的最小二乘优化模型,并根据含初值的最小二乘优化模型计算出支路末端节点的阻抗。
建立含初值的最小二乘优化模型的步骤包括:
3.1)设参数状态变量x包括电压、功率,以常数量k替换电阻、电抗,即:
式中,x1、x2、x3、x4、x5、x6为状态变量。k1、k2、k3、k4、k5、k6为常数量。
3.2)更新参数fi(Ri,Xi)、参数fi(Rj,Xj),得到:
3.3)将公式(5)变换成线性回归方程,得到:
3.4)建立配变局部变量的多元线性回归模型,即:
式中,y和xi,i=1,2,…,6为T维列向量。T为历史正常断面数。
3.5)根据公式(10),建立含初值的最小二乘优化模型,即:
Constraints=Rj,Xj,Ri,Xi>0
式中,obj表示模型目标函数,Constraints表示约束条件。
3.6)根据含初值的最小二乘优化模型计算出末端节点的阻抗,步骤包括:
3.6.1)利用配变局部变量多元非线性模型计算出常数量k1、k2、k3、k4、k5、k6的初值。
3.6.2)将常数量k1、k2、k3、k4、k5、k6的初值输入到含初值的最小二乘优化模型中,求解得到支路末端节点i的电阻Ri和电抗Xi、支路末端节点j的电阻Ri和电抗Xi。
求解含初值的最小二乘优化模型的方法包括Levenberg-Marquardt算法(列文伯格-马夸尔特法)。
4)递推计算得到配电网线路参数。
所述配电网线路参数包括馈线末端支路Zend的状态变量Umid+1,Pmid+1,Qmid+1,Uend,Pend,Qend、馈线中段支路Zmid的状态变量Umid,Umid+1,Pmid+1,Qmid+1。为馈线末端支路Yend首端的有功功率和无功功率;Pend,Qend为馈线末端支路Yend末端的有功功率和无功功率。
递推计算得到配电网线路参数的步骤包括:
4.1)量测馈线末端支路Zend的功率量测量Umid+1,Pmid+1,Qmid+1,Uend,Pend,Qend,并根据配变局部变量多元非线性模型和含初值的最小二乘优化模型计算得到馈线末端支路Zend的电阻Rend和电抗Xend。其中,Umid+1、Uend为电压。Pmid+1、Pend为有功功率。Qmid+1、Qend为无功功率。
4.2)计算馈线末端支路Zend首端的有功功率和无功功率/>即:
4.3)从馈线末端支路Zend开始,自下而上逐级进行参数估计、循环变量递推计算得到全网的拓扑及支路参数信息。
实施例2:
基于线性回归初值估计的配电网线路参数辨识方法,包括以下步骤:
1)获取电网的基础数据
所述电网基础数据包括电网拓扑结构、所有可观测节点数量Nobs、所有节点的有功功率量测量Pi、节点的无功功率量测量Qi、节点的电压幅值量测Vi。
2)构建配变局部变量多元非线性模型
V0为始端电压,V1、V2为待求末端电压。根据图1所示线路电压降相量图,可得节点0、1的电压幅值关系为:
上面二式减可得并联节点对(1,2)的电气参量准确关系式:
根据上述关系可得N个时刻下,任意支路末端节点i,j的连接判据为:
3)构建含初值的最小二乘优化模型
将用y替换;电压、功率用状态变量x替换;电阻、电抗用常数变量k替换。可得:
即(5)可以变换成线性回归方程:
于是可构建配变局部变量的多元线性回归模型为:
式中,y和xi,i=1,2,…,6为T维列向量,由配变局部状态变量的历史正常断面量测Ui,Pi,Qi,Uj,Pj,Qj按式(10)。求得,其中T为历史正常断面数。求解上述多元线性回归模型即可估计出回归系数进而通过求解式(7)即可对等值模型中两支路的阻抗Ri,Xi,Rj,Xj进行精细化估计。
由式(10)可以构建含初值的最小二乘优化模型:
Constraints=Rj,Xj,Ri,Xi>0
为步骤2计算得到的初值,带入式(11)用Levenberg-Marquardt算法求解上述含初值的优化模型,即可求得精细化辨识的Rj,Xj,Ri,Xi。
4)配变局部变量多元线性回归模型的递推计算
如图2所示,配变局部变量多元线性回归模型参数估计的过程中,对于馈线末端支路Zend,变量Umid+1,Pmid+1,Qmid+1,Uend,Pend,Qend可分别由配变mid+1和配变end的量测已知,将其代入式(10),然后求解式(10)即可估计末端支路Zend的参数。然而对于处于馈线中段支路Zmid的Umid,Umid+1可由配变mid和配变mid+1的量测得知,而馈线中段支路Zmid末端量测Pmid+1,Qmid+1需要叠加馈线末端支路Zend支路的功率损耗,无法直接对Zmid的参数进行估计。为此需要首先对馈线末端支路Zend的功率损耗进行计算。
由于馈线末端两配变的电压,功率量测Umid+1,Pmid+1,Qmid+1,Uend,Pend,Qend是已知的。首先执行馈线末端支路Zend的参数估计,得到Rend,Xend。再计算馈线末端支路Zend首端的功率 如下式所示,替换原功率Pmid+1,Qmid+1。之后再执行步骤1所述的两两支路的参数估计。
因此,只需从馈线末端支路Zend开始,自下而上逐级进行参数估计、循环变量递推计算即可得到全网的拓扑及支路参数信息,总体方法流程见图3所示。
实施例3:
参见图4,基于线性回归初值估计的配电网线路参数辨识方法的验证试验,内容包括:
为验证实施例1-2所提计及对地电容并考虑线路长度估计的配电网线路参数辨识正确性及有效性,故针对IEEE 33节点系统进行仿真分析。该系统额定电压为12.66kV,线路均为电缆线路,共有31个负荷端点,0号节点为电源点,5号节点为无量测的汇集节点。其中,负荷数据采用某配网实际用户数据,将其规整到33节点系统标准算例的负荷水平,保留其负荷波动特性。然后,利用计及对地电容并考虑线路长度估计的配电网线路参数辨识方法辨识IEEE 33节点系统各线路参数。
1)构建配变局部变量多元非线性模型
上面二式减可得并联节点对(1,2)的电气参量准确关系式:
根据上述关系可得N个时刻下,任意支路末端节点i,j的连接判据为:
2)构建含初值的最小二乘回归模型
将用y替换;电压、功率用状态变量x替换;电阻、电抗用常数变量k替换。可得:
即(5)可以变换成线性回归方程:
于是可构建配变局部变量的多元线性回归模型为:
由式(10)可以构建含初值的最小二乘优化模型:
Constraints=Rj,Xj,Ri,Xi>0
3)多元非线性电压降落方程的递推计算
/>
表1 33rd配电网线路参数辨识结果
/>
可以看出,通过线性回归获得初值的最小二乘算法辨识精度有较好的提升。
实施例4:
参见图5,基于线性回归初值估计的配电网线路参数辨识方法的验证试验,内容包括:
33rd网络与33rds网络区别在于:26节点是无负荷的串联节点,因此并不会出现在结果中,并且各支路的阻抗值与33rds网络是不同。
1)构建配变局部变量多元非线性模型
上面二式减可得并联节点对(1,2)的电气参量准确关系式:
根据上述关系可得N个时刻下,任意支路末端节点i,j的连接判据为:
2)构建含初值的最小二乘回归模型
将用y替换;电压、功率用状态变量x替换;电阻、电抗用常数变量k替换。可得:
即(5)可以变换成线性回归方程:
于是可构建配变局部变量的多元线性回归模型为:
由式(0.1)可以构建含初值的最小二乘优化模型:
Constraints=Rj,Xj,Ri,Xi>0
3)多元非线性电压降落方程的递推计算
表2 33rds配电网线路参数辨识结果
/>
可以看出33节点rds网络,更换了线路阻抗数值及包含串联无负荷节点,误差在允许范围之内。
实施例5:
参见图6,基于线性回归初值估计的配电网线路参数辨识方法的验证试验,内容包括:
与33rd,33rds网络区别在于:69rd系统的节点数更多;且含24、22、18、14等无负荷的串联节点和2、3等无负荷并联节点;辨识难度最大。因为无负荷节点,即为无量测点,这些节点的信息需要从其他节点推断得出。
1)构建配变局部变量多元非线性模型
上面二式减可得并联节点对(1,2)的电气参量准确关系式:
根据上述关系可得N个时刻下,任意支路末端节点i,j的连接判据为:
2)构建含初值的最小二乘回归模型
将用y替换;电压、功率用状态变量x替换;电阻、电抗用常数变量k替换。可得:
/>
即(5)可以变换成线性回归方程:
于是可构建配变局部变量的多元线性回归模型为:
由式(10)可以构建含初值的最小二乘优化模型:
Constraints=Rj,Xj,Ri,Xi>0
3)多元非线性电压降落方程的递推计算
表3 69rd配电网线路参数辨识结果
/>
/>
/>
可以看出69节点rd网络,增加了节点数量以及更换了线路阻抗数值及包含串联、并联无负荷节点。部分线路辨识误差较33节点系统增加,但误差仍然在允许范围之内。
Claims (5)
1.基于线性回归初值估计的配电网线路参数辨识方法,其特征在于,包括以下步骤:
1)获取电网的基础数据;
2)建立配变局部变量多元非线性模型;
3)建立含初值的最小二乘优化模型,并根据含初值的最小二乘优化模型计算出支路末端节点的阻抗;
4)递推计算得到配电网线路参数;
建立配变局部变量多元非线性模型的步骤包括:
2.1)建立始端电压V0和待求末端电压V1、待求末端电压V2之间的关系式,即:
式中,P1、P2分别为待求末端节点1和待求末端节点2的有功功率;Q1、Q2分别为待求末端节点1和待求末端节点2的无功功率;R1、R2分别为待求末端节点1和待求末端节点2的等效电阻;X1、X2分别为待求末端节点1和待求末端节点2的等效电抗;
2.2)根据公式(1)建立并联节点对待求末端节点1和待求末端节点2的电气参量准确关系式,即:
V1 2-V2 2=-f1(R1,X1)+f2(R2,X2) (2)
式中,参数f1(R1,X1)、参数f1(R1,X1)分别如下所示:
2.3)建立N个时刻下,任意支路末端节点i和支路末端节点j的连接判据,即:
式中,Ri、Rj分别为支路末端节点i和支路末端节点j的等效电阻;Xi、Xj分别为支路末端节点i和支路末端节点j的等效电抗;
建立含初值的最小二乘优化模型的步骤包括:
3.1)设参数状态变量x包括电压、功率,以常数量k替换电阻、电抗,即:
式中,x1、x2、x3、x4、x5、x6为状态变量;k1、k2、k3、k4、k5、k6为常数量;
3.2)更新参数fi(Ri,Xi)、参数fj(Rj,Xj),得到:
3.3)将公式(5)变换成线性回归方程,得到:
3.4)建立配变局部变量的多元线性回归模型,即:
式中,y和xi,i=1,2,…,6为T维列向量;T为历史正常断面数;
3.5)根据公式(10),建立含初值的最小二乘优化模型,即:
Constraints=Rj,Xj,Ri,Xi>0
根据含初值的最小二乘优化模型计算出末端节点的阻抗的步骤包括:
s1)利用配变局部变量多元非线性模型计算出常数量k1、k2、k3、k4、k5、k6的初值;
s2)将常数量k1、k2、k3、k4、k5、k6的初值输入到含初值的最小二乘优化模型中,求解得到支路末端节点i的电阻Ri和电抗Xi、支路末端节点j的电阻Ri和电抗Xi。
2.根据权利要求1所述的基于线性回归初值估计的配电网线路参数辨识方法,其特征在于,所述电网基础数据包括电网拓扑结构、可观测节点数量Nobs、节点的有功功率量测量Pi、节点的无功功率量测量Qi、节点的电压幅值量测Vi。
3.根据权利要求1所述的基于线性回归初值估计的配电网线路参数辨识方法,其特征在于,求解含初值的最小二乘优化模型的方法包括Levenberg-Marquardt算法。
4.根据权利要求3所述的基于线性回归初值估计的配电网线路参数辨识方法,其特征在于,所述配电网线路参数包括馈线末端支路Yend的状态变量Umid+1,Pmid+1,Qmid+1,Uend,Pend,Qend、馈线中段支路Ymid的状态变量Umid,Umid+1,Pmid+1,Qmid+1;为馈线末端支路Yend首端的有功功率和无功功率;Pend,Qend为馈线末端支路Yend末端的有功功率和无功功率。
5.根据权利要求4所述的基于线性回归初值估计的配电网线路参数辨识方法,其特征在于,递推计算得到配电网线路参数的步骤包括:
1)量测馈线末端支路Zend的功率量测量Umid+1,Pmid+1,Qmid+1,Uend,Pend,Qend,并根据配变局部变量多元非线性模型和含初值的最小二乘优化模型计算得到馈线末端支路Zend的电阻Rend和电抗Xend;其中,Umid+1、Uend为电压;Pmid+1、Pend为有功功率;Qmid+1、Qend为无功功率;
2)计算馈线末端支路Yend首端的有功功率和无功功率/>即:
3)从馈线末端支路Zend开始,自下而上逐级进行参数估计、循环变量递推计算得到全网的拓扑及支路参数信息。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210862077.4A CN115347553B (zh) | 2022-07-20 | 2022-07-20 | 基于线性回归初值估计的配电网线路参数辨识方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210862077.4A CN115347553B (zh) | 2022-07-20 | 2022-07-20 | 基于线性回归初值估计的配电网线路参数辨识方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115347553A CN115347553A (zh) | 2022-11-15 |
CN115347553B true CN115347553B (zh) | 2024-05-24 |
Family
ID=83950494
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210862077.4A Active CN115347553B (zh) | 2022-07-20 | 2022-07-20 | 基于线性回归初值估计的配电网线路参数辨识方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115347553B (zh) |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112600191A (zh) * | 2020-10-29 | 2021-04-02 | 河北工业大学 | 基于非线性回归的低压配电网参数和拓扑联合辨识方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP7021185B2 (ja) * | 2016-04-22 | 2022-02-16 | デプシス ソシエテ アノニム | 電力ネットワークの複数の測定ノード間の相互電圧感度係数を求める方法 |
-
2022
- 2022-07-20 CN CN202210862077.4A patent/CN115347553B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112600191A (zh) * | 2020-10-29 | 2021-04-02 | 河北工业大学 | 基于非线性回归的低压配电网参数和拓扑联合辨识方法 |
Non-Patent Citations (1)
Title |
---|
一种基于最小二乘法的配网设备参数辨识方法;吴凯;袁慧宏;杨斌;王瑶;张怡心;王龙;郭松;;能源与节能;20190520(第05期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN115347553A (zh) | 2022-11-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107453357B (zh) | 一种基于分层求解的配电网状态估计方法 | |
KR101219545B1 (ko) | 전력계통에서의 최적화 기법을 적용한 파라미터 추정 방법 | |
CN110705107B (zh) | 一种配电网电压评估方法、系统、设备和存储介质 | |
CN114184876B (zh) | 一种直流偏磁监测、评价及大地模型校正平台 | |
CN109752629B (zh) | 一种电网量测问题智能诊断方法及系统 | |
Dzafic et al. | High performance state estimation for smart grid distribution network operation | |
CN112288303A (zh) | 确定线损率的方式、装置 | |
CN111969604A (zh) | 基于实测数据的台区线损动态计算方法及装置 | |
CN110350515A (zh) | 一种适用于故障电流分析的柔性直流电网建模方法 | |
CN115685046A (zh) | 互感器计量异常识别方法、装置、设备及存储介质 | |
CN111553080A (zh) | 一种配电台区负荷动态等值非机理模型参数闭环辨识方法 | |
CN115000947A (zh) | 基于智能电表量测的配电网拓扑结构与线路参数辨识方法 | |
CN115347553B (zh) | 基于线性回归初值估计的配电网线路参数辨识方法 | |
CN110867861B (zh) | 一种配电网理论线损精确定界方法 | |
CN110907720B (zh) | 基于pmu量测的短线路同杆并架双回线完整参数辨识方法 | |
CN105977963B (zh) | 一种基于两阶段法的主配网一体化分布式状态估计方法 | |
CN110865328B (zh) | 基于ami的智能电表相位识别、拓扑辨识及阻抗估计方法 | |
CN109586289B (zh) | 一种配电网多时间尺度递归动态状态估计方法及系统 | |
CN116805029A (zh) | 一种基于电网直流偏磁电流实测的大地电位分布计算方法、装置、系统及存储介质 | |
CN113435065B (zh) | 一种基于“黑箱”dc/dc变换器的等效阻抗建模方法 | |
CN107767060B (zh) | 配网线路理论线损计算系统及方法 | |
CN115333083A (zh) | 计及对地电容并考虑线路长度估计的配电网线路参数辨识方法 | |
CN114626959A (zh) | 考虑多源量测的配电网二阶锥规划鲁棒状态估计方法 | |
CN109270347B (zh) | 一种针对同步相量量测的电网可观测性分析方法 | |
Cataliotti et al. | Uncertainty evaluation of a backward/forward load flow algorithm for a MV smart grid |
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 |