CN107749627A - 基于改进匹配追踪的智能配电网潮流雅可比矩阵估计方法 - Google Patents
基于改进匹配追踪的智能配电网潮流雅可比矩阵估计方法 Download PDFInfo
- Publication number
- CN107749627A CN107749627A CN201711118202.6A CN201711118202A CN107749627A CN 107749627 A CN107749627 A CN 107749627A CN 201711118202 A CN201711118202 A CN 201711118202A CN 107749627 A CN107749627 A CN 107749627A
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- vector
- matrix
- mover
- 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.)
- Granted
Links
- 239000011159 matrix material Substances 0.000 title claims abstract description 144
- 238000000034 method Methods 0.000 title claims abstract description 34
- 239000013598 vector Substances 0.000 claims abstract description 156
- 238000005259 measurement Methods 0.000 claims description 83
- 238000004364 calculation method Methods 0.000 claims description 15
- 238000011084 recovery Methods 0.000 claims description 13
- 230000001360 synchronised effect Effects 0.000 claims description 8
- 230000035945 sensitivity Effects 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005457 optimization Methods 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
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- 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)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
一种基于改进匹配追踪的智能配电网潮流雅可比矩阵估计方法：获取配电网的节点数，输入网络最大度的保守估计值，设定残差阈值及最大迭代次数；获取系统各节点的当前时刻和历史量测数据；将历史量测数据分别与当前时刻测数据做差，得到多组变化向量来构造传感矩阵；得到电压相角和电压幅值的相关系数向量；更新传感矩阵列号索引集合；利用传感矩阵列号索引集合求最小二乘解，更新残差向量；判断更新后的残差向量的2范数是否小于残差阈值；得到残差相关系数向量，更新传感矩阵列号索引集合；求最小二乘解，重新构造列号索引集合，再次更新残差向量；输出最小二乘估计结果，根据传感矩阵列号索引集合恢复出2N维向量作为雅可比矩阵一个元素的估计结果并输出。
Description
技术领域
本发明涉及一种雅可比矩阵估计方法。特别是涉及一种基于改进匹配追踪的智能配电网潮流雅可比矩阵估计方法。
背景技术
大量分布式电源、电动汽车、电力电子装置的接入，给配电网的运行带来更大灵活性的同时，也对配电网的可观性和可控性提出了更高的要求。同步相量量测技术在配电网层面的发展，能够有效提高智能配电网的可观性，也为解决智能配电网的运行控制问题提供了新的思路。通过潮流雅可比矩阵求逆可以得到电压-功率灵敏度关系，而电压-功率灵敏度能够准确反映系统中电压与功率的变化关系，在配电系统的运行控制中发挥着重要的作用。采用离线潮流计算的电力系统灵敏度参数，由于存在元件参数不准确、信息更新不及时、难以追踪系统运行点及相关拓扑变化等缺陷，往往会造成较大的计算误差，极大地影响了电压控制的效果。同步相量量测技术的发展使得通过量测数据实现系统灵敏度参数的在线估计成为可能。
同步相量量测单元能够实现有功功率、无功功率、电压相角、电压幅值、系统频率等电气量的高精度同步量测，已经广泛应用于电力系统的参数辨识与运行控制等方面。通过同步相量量测装置采集的多时间段面的量测信息，采用最小二乘估计方法，能够实现潮流雅可比矩阵的精确估计。同时根据潮流方程的特点，雅可比矩阵本身具有较大的稀疏性，在参数估计问题中考虑稀疏性能够在减少量测组数需求的同时提高估计的精度。利用稀疏性估计系统参数，需要将原最小二乘估计问题转换为稀疏恢复问题，正交匹配追踪类的贪婪算法是解决稀疏恢复问题的有效手段。但是现有的正交匹配追踪类算法在求解稀疏恢复问题时，要求传感矩阵满足有限等距条件，恢复的成功率依赖于传感矩阵中各列间具有较小的相关性，但是雅可比矩阵估计问题的传感矩阵中各个节点电压变化向量间具有较强的相关性，从而造成采用一般稀疏恢复方法在求解恢复问题时，在量测组数较少的情况下无法保证恢复的成功率。因此结合雅可比矩阵的特殊性，需要提出更有针对性的改进方法，以进一步提高求解效率，降低问题对量测组数的依赖性。
发明内容
本发明所要解决的技术问题是，提供一种在量测数据较少的情况下能够实现潮流雅可比矩阵精确估计的基于改进匹配追踪的智能配电网潮流雅可比矩阵估计方法。
本发明所采用的技术方案是：一种基于改进匹配追踪的智能配电网潮流雅可比矩阵估计方法，包括如下步骤：
1)获取配电网的节点数，将源节点编号为0，其他节点依次编号为1,…,i,…,N，输入网络最大度的保守估计值dmax，设定残差阈值ε及最大迭代次数M；
2)获取系统各节点同步相量量测装置的有功功率、无功功率、电压幅值和电压相角当前时刻的量测数据和C组历史量测数据，其中，C是大于1的整数；
3)将C组历史量测数据分别与当前量测值做差，得到C组有功功率、无功功率、电压幅值和电压相角的变化向量，利用节点1～节点N电压量测的变化向量构造传感矩阵，初始化潮流雅可比矩阵的行号m＝1；
4)若m不大于N，选取节点i的有功功率量测变化向量初始化残差向量，分别计算传感矩阵的第q列和第q+N列与传感矩阵所有列的相关性，得到电压相角的相关系数向量uθ和电压幅值的相关系数向量uU，此时i＝m，q＝m；若m大于N则用节点i的无功功率量测变化向量初始化残差向量，分别计算传感矩阵的第q-N列和第q列与传感矩阵所有列的相关性，得到电压相角的相关系数向量uθ和电压幅值的相关系数向量uU，此时i＝m-N，q＝m；初始化传感矩阵列号索引集合为空集，初始化迭代次数n＝1；
5)分别选取电压相角的相关系数向量uθ和电压幅值的相关系数向量uU中最大的z个数值，其中z＝dmax+1，将z个数值对应的传感矩阵中的列号索引构成集合Ωn，更新传感矩阵列号索引集合；
6)利用传感矩阵列号索引集合求最小二乘解，更新残差向量；
7)若更新后的残差向量的2范数R小于残差阈值ε，则进入步骤11)；否则：如果迭代次数n＝1，分别选取电压相角的相关系数向量uθ和电压幅值的相关系数向量uU中最大的2个值，将2个值对应的传感矩阵中的索引构成集合Ωn，更新列号索引集合，迭代次数n＝n+1，返回步骤6)；如果迭代次数n＝2，则迭代次数n＝n+1，进入步骤8)；
8)分别计算更新后的残差向量与传感矩阵各列的相关性，得到残差相关系数向量u，选取残差相关系数向量u中最大的2列，将2列对应的传感矩阵中的索引构成集合Ωn，更新传感矩阵列号索引集合；
9)利用传感矩阵列号索引集合求最小二乘解，选取最小二乘解中绝对值最大的4项，将4项对应的传感矩阵中的索引更新集合Ωn，重新构造列号索引集合，再次更新残差向量；
10)若再次更新的残差向量的2范数R小于残差阈值ε或者迭代次数超过设定的最大迭代次数M，则进入步骤11)；否则n＝n+1，返回步骤8)；
11)输出最小二乘估计结果，根据传感矩阵列号索引集合恢复出2N维向量作为雅可比矩阵第m行的估计结果，m＝m+1；若m大于2N，停止迭代，输出雅可比矩阵估计结果，否则，返回步骤4)。
步骤3)中
(1)所述的C组有功功率、无功功率、电压幅值和电压相角的变化向量表示为：
ΔPi[k]＝Pi(k)-Pi(0)、ΔQi[k]＝Qi(k)-Qi(0)、ΔVi[k]＝Vi(k)-Vi(0)和Δθi[k]＝θi(k)-θi(0)，k＝1,2,…,C，Pi(0)、Qi(0)、θi(0)、Vi(0)分别表示节点i当前时刻的有功功率、无功功率、电压相角和电压幅值的量测值；Pi(k)、Qi(k)、θi(k)、Vi(k)分别表示节点i第k组的历史量测值；
(2)所述构造传感矩阵A如下：
式中，表示由电压相角和电压幅值量测变化向量构成的矩阵，Δθi＝[Δθi[1],…,Δθi[C]]T表示节点i的C组电压相角量测变化量组成的列向量，ΔVi＝[ΔVi[1],…,ΔVi[C]]T表示节点i的C组电压幅值量测变化量组成的列向量；表示中的元素，Ap,q表示传感矩阵A中第p行第q列的元素。
步骤4)中
(1)所述的选取节点i的有功功率量测变化向量初始化残差向量表示为：
r0＝ΔPi
式中，r0表示初始的残差向量，ΔPi＝[ΔPi[1],…,ΔPi[C]]T表示节点i的C组有功功率量测变化量组成的列向量；
(2)用节点i的无功功率量测变化向量初始化残差向量表示为：
r0＝ΔQi
式中，ΔQi＝[ΔQi[1],…,ΔQi[C]]T表示节点i的C组无功功率量测变化量组成的列向量。
(3)电压相角的相关系数向量uθ的计算方法为：
若1≤m≤N：
uθ＝abs[ATAq]
若N＜m≤2N：
uθ＝abs[ATAq-N]
式中，abs[·]表示取绝对值运算，A为传感矩阵，Aq和Aq-N分别表示矩阵A的第q和q-N列；
(4)电压幅值的相关系数向量uU的计算方法为：
若1≤m≤N：
uU＝abs[ATAq+N]
若N＜m≤2N：
uU＝abs[ATAq]
步骤5)所述的更新传感矩阵列号索引集合表示为：
Λn＝Λn-1∪Ωn
式中，Λn为第n次迭代的列号索引集合，当n＝1时，Λn-1表示初始的列号索引集合，Ωn表示新选入的索引构成的集合。
步骤6)包括：
(1)最小二乘解表示为：
式中，表示第n次迭代时的最小二乘解，r0为初始的残差向量，表示第n次迭代的列号索引集合Λn中元素对应的传感矩阵各列构成的矩阵；
(2)更新后的残差向量表示为：
式中，rn为第n次迭代的残差向量。
步骤8)所述的残差相关系数向量u计算方式为：
u＝abs[ATrn-1]
式中，u为相关系数向量，abs[·]表示取绝对值运算，rn-1表示第n-1次迭代时的残差向量，当n＝1时，rn-1表示初始残差向量，A为传感矩阵。
步骤9)中
(1)所述的重新构造列号索引集合为
Λn＝Ωn
(2)所述的再次更新残差向量表示为：
式中，Λn为第n次迭代的列号索引集合，Ωn表示新选入的索引构成的集合，表示集合Λn中元素对应的传感矩阵各列构成的矩阵，rn为第n次迭代的残差向量，r0为初始的残差向量。
步骤11)所述的根据传感矩阵列号索引集合恢复出2N维向量作为雅可比矩阵第m行的估计结果表示为：
式中，表示第n次迭代时的最小二乘解，表示2N维恢复向量，表示的第g个元素，表示列号索引集合Λn中各元素所对应的向量的元素组成的向量，对于其他不在列号索引集合中的元素取值为0；矩阵第q列的模值， 表示雅可比矩阵Y的第m行第q个元素的估计解，其中，g＝q。
本发明的基于改进匹配追踪的智能配电网潮流雅可比矩阵估计方法，利用了矩阵的稀疏性质，采用压缩感知技术将估计问题转化成稀疏恢复问题，通过求解稀疏恢复问题，以少量的量测组数实现了雅可比矩阵的精确估计。并且针对现有稀疏恢复算法在求解雅可比矩阵估计问题的局限性，考虑电压幅值和相角量测变化向量之间的相关关系对算法进行了改进，进一步减少了对量测组数的要求，提高了估计的成功率和计算效率。
附图说明
图1是本发明基于改进匹配追踪的智能配电网潮流雅可比矩阵估计方法的流程图；
图2是IEEE 33节点算例图。
具体实施方式
下面结合实施例和附图对本发明的基于改进匹配追踪的智能配电网潮流雅可比矩阵估计方法做出详细说明。
如图1所示,本发明的基于改进匹配追踪的智能配电网潮流雅可比矩阵估计方法，包括如下步骤：
1)，获取配电网的节点数，将源节点编号为0，其他节点依次编号为1,…,i,…,N，输入网络最大度的保守估计值dmax，设定残差阈值ε及最大迭代次数M；
2)获取系统各节点同步相量量测装置的有功功率、无功功率、电压幅值和电压相角当前时刻的量测数据和C组历史量测数据，其中，C是大于1的整数；
所述的C组历史量测数据生成方法如下：
(1)采用如下的公式生成节点i的第k组有功功率量测数据,
式中，Pi(k)表示节点i的第k组有功功率，Pi(0)表示当前节点i的有功功率量测，是服从均值为0正态分布的随机数，分别用来模拟不同量测时刻相对于当前时刻的功率变化和量测误差。
(2)采用如下公式生成节点i的第k组无功量测数据,
式中，Qi(k)表示节点i的第k组无功功率，Qi(0)表示当前节点i的无功功率量测，是服从均值为0正态分布的随机数。
(3)在得到节点i的第k组有功功率和无功功率量测数据后，通过潮流计算求得对应的电压相角θi(k)和幅值Vi(k)作为节点i的第k组电压相角和幅值量测值。
3)将C组历史量测数据分别与当前量测值做差，得到C组有功功率、无功功率、电压幅值和电压相角的变化向量，利用节点1～节点N电压量测的变化向量构造传感矩阵，初始化潮流雅可比矩阵的行号m＝1；其中,
所述的C组有功功率、无功功率、电压幅值和电压相角的变化向量表示为：
ΔPi[k]＝Pi(k)-Pi(0)、ΔQi[k]＝Qi(k)-Qi(0)、ΔVi[k]＝Vi(k)-Vi(0)和Δθi[k]＝θi(k)-θi(0)，k＝1,2,…,C，Pi(0)、Qi(0)、θi(0)、Vi(0)分别表示节点i当前时刻的有功功率、无功功率、电压相角和电压幅值的量测值；Pi(k)、Qi(k)、θi(k)、Vi(k)分别表示节点i第k组的历史量测值；
(2)所述构造传感矩阵A如下：
式中，表示由电压相角和电压幅值量测变化向量构成的矩阵，Δθi＝[Δθi[1],…,Δθi[C]]T表示节点i的C组电压相角量测变化量组成的列向量，ΔVi＝[ΔVi[1],…,ΔVi[C]]T表示节点i的C组电压幅值量测变化量组成的列向量；表示中的元素，Ap,q表示传感矩阵A中第p行第q列的元素。
4)若m不大于N，选取节点i的有功功率量测变化向量初始化残差向量，分别计算传感矩阵的第q列和第q+N列与传感矩阵所有列的相关性，得到电压相角的相关系数向量uθ和电压幅值的相关系数向量uU，此时i＝m，q＝m；若m大于N则用节点i的无功功率量测变化向量初始化残差向量，分别计算传感矩阵的第q-N列和第q列与传感矩阵所有列的相关性，得到电压相角的相关系数向量uθ和电压幅值的相关系数向量uU，此时i＝m-N，q＝m；初始化传感矩阵列号索引集合为空集，初始化迭代次数n＝1；其中，
(1)所述的选取节点i的有功功率量测变化向量初始化残差向量表示为：
r0＝ΔPi
式中，r0表示初始的残差向量，ΔPi＝[ΔPi[1],…,ΔPi[C]]T表示节点i的C组有功功率量测变化量组成的列向量；
(2)用节点i的无功功率量测变化向量初始化残差向量表示为：
r0＝ΔQi
式中，ΔQi＝[ΔQi[1],…,ΔQi[C]]T表示节点i的C组无功功率量测变化量组成的列向量。
(3)电压相角的相关系数向量uθ的计算方法为：
若1≤m≤N：
uθ＝abs[ATAq]
若N＜m≤2N：
uθ＝abs[ATAq-N]
式中，abs[·]表示取绝对值运算，A为传感矩阵，Aq和Aq-N分别表示矩阵A的第q和q-N列；
(4)电压幅值的相关系数向量uU的计算方法为：
若1≤m≤N：
uU＝abs[ATAq+N]
若N＜m≤2N：
uU＝abs[ATAq]
5)分别选取电压相角的相关系数向量uθ和电压幅值的相关系数向量uU中最大的z个数值，其中z＝dmax+1，将z个数值对应的传感矩阵中的列号索引构成集合Ωn，更新传感矩阵列号索引集合；
所述的更新传感矩阵列号索引集合表示为：
Λn＝Λn-1∪Ωn
式中，Λn为第n次迭代的列号索引集合，当n＝1时，Λn-1表示初始的列号索引集合，Ωn表示新选入的索引构成的集合。
6)利用传感矩阵列号索引集合求最小二乘解，更新残差向量；包括：
(1)最小二乘解表示为：
式中，表示第n次迭代时的最小二乘解，r0为初始的残差向量，表示第n次迭代的列号索引集合Λn中元素对应的传感矩阵各列构成的矩阵；
(2)更新后的残差向量表示为：
式中，rn为第n次迭代的残差向量。
7)若更新后的残差向量的2范数R小于残差阈值ε，则进入步骤11)；否则：如果迭代次数n＝1，分别选取电压相角的相关系数向量uθ和电压幅值的相关系数向量uU中最大的2个值，将2个值对应的传感矩阵中的索引构成集合Ωn，更新列号索引集合，迭代次数n＝n+1，返回步骤6)；如果迭代次数n＝2，则迭代次数n＝n+1，进入步骤8)；
8)分别计算更新后的残差向量与传感矩阵各列的相关性，得到残差相关系数向量u，选取残差相关系数向量u中最大的2列，将2列对应的传感矩阵中的索引构成集合Ωn，更新传感矩阵列号索引集合；
所述的残差相关系数向量u计算方式为：
u＝abs[ATrn-1]
式中，u为相关系数向量，abs[·]表示取绝对值运算，rn-1表示第n-1次迭代时的残差向量，当n＝1时，rn-1表示初始残差向量，A为传感矩阵。
9)利用传感矩阵列号索引集合求最小二乘解，选取最小二乘解中绝对值最大的4项，将4项对应的传感矩阵中的索引更新集合Ωn，重新构造列号索引集合，再次更新残差向量；其中,
所述的重新构造列号索引集合为
Λn＝Ωn
所述的再次更新残差向量表示为：
式中，Λn为第n次迭代的列号索引集合，Ωn表示新选入的索引构成的集合，表示集合Λn中元素对应的传感矩阵各列构成的矩阵，rn为第n次迭代的残差向量，r0为初始的残差向量。
10)若再次更新的残差向量的2范数R小于残差阈值ε或者迭代次数超过设定的最大迭代次数M，则进入步骤11)；否则n＝n+1，返回步骤8)；
11)输出最小二乘估计结果，根据传感矩阵列号索引集合恢复出2N维向量作为雅可比矩阵第m行的估计结果，m＝m+1；若m大于2N，停止迭代，输出雅可比矩阵估计结果，否则，返回步骤4)。
所述的根据传感矩阵列号索引集合恢复出2N维向量作为雅可比矩阵第m行的估计结果表示为：
式中，表示第n次迭代时的最小二乘解，表示2N维恢复向量，表示的第g个元素，表示列号索引集合Λn中各元素所对应的向量的元素组成的向量，对于其他不在列号索引集合中的元素取值为0；矩阵第q列的模值， 表示雅可比矩阵Y的第m行第q个元素的估计解，其中，g＝q。
下面给出具体实例：
首先输入IEEE 33节点算例网络拓扑连接关系如图2所示，其中节点0为平衡节点，其他节点1～32为PQ节点，系统的基准容量为1MVA，基准电压为12.66kV，各个PQ节点的当前功率量测如表1所示。输入网络的最大度的估计值为4，模拟量测功率变化和误差随机数的标准差分别设为0.01和0.025％。采用下式计算雅可比矩阵的误差。
式中，分别表示雅可比矩阵参数的估计值，Ji,j为采用精确线路参数的计算值。
为验证本发明方法的先进性，采取如下两种场景进行分析：
场景1，设置量测组数为30、35、40、45、50、55、60，分别采用最小二乘方法、正交匹配追踪算法和本专利所提出算法进行雅可比矩阵估计；
场景2，设置量测组数为40、45、50、55、60，采用本发明的方法对雅可比矩阵进行估计；设置量测组数为100、500、1000、1500、2000，直接采用最小二乘方法进行估计。
场景1估计结果如表2所示，场景2估计结果如表3和表4所示。
执行优化计算的计算机硬件环境为Intel(R)Xeon(R)CPU E5-1620，主频为3.70GHz，内存为32GB；软件环境为Windows 7操作系统，采用MATLAB的MATPOWER工具包计算潮流。
从表2中可以看出，当量测组数少于60组时最小二乘计算方法无法实现雅可比矩阵任意一行的估计，正交匹配追踪算法能够利用雅可比矩阵的稀疏性，实现雅可比矩阵中一部分行的估计，而且随着量测组数的增加，能够成功估计的行数增加，但是即使量测组数增加为60组时，仍然不能实现所有行的估计；本发明的方法在量测组数增加到40时，即可实现雅可比矩阵所有行的估计，减少了估计问题对量测组数的依赖，提高了估计的成功率。通过表3和表4的对比可以看出，本发明所提出的估计方法在减少对量测组数依赖的同时也能够保证估计的精确性，当量测组数为60组时的估计精度与采用最小二乘方法2000组量测的估计精度基本相同。
表1IEEE 33节点算例PQ节点当前功率量测值
表2场景1的估计结果
表3场景2的估计结果
表4场景2的估计结果
Claims (8)
1.一种基于改进匹配追踪的智能配电网潮流雅可比矩阵估计方法，其特征在于，包括如下步骤：
1)获取配电网的节点数，将源节点编号为0，其他节点依次编号为1,…,i,…,N，输入网络最大度的保守估计值dmax，设定残差阈值ε及最大迭代次数M；
2)获取系统各节点同步相量量测装置的有功功率、无功功率、电压幅值和电压相角当前时刻的量测数据和C组历史量测数据，其中，C是大于1的整数；
3)将C组历史量测数据分别与当前量测值做差，得到C组有功功率、无功功率、电压幅值和电压相角的变化向量，利用节点1～节点N电压量测的变化向量构造传感矩阵，初始化潮流雅可比矩阵的行号m＝1；
4)若m不大于N，选取节点i的有功功率量测变化向量初始化残差向量，分别计算传感矩阵的第q列和第q+N列与传感矩阵所有列的相关性，得到电压相角的相关系数向量uθ和电压幅值的相关系数向量uU，此时i＝m，q＝m；若m大于N则用节点i的无功功率量测变化向量初始化残差向量，分别计算传感矩阵的第q-N列和第q列与传感矩阵所有列的相关性，得到电压相角的相关系数向量uθ和电压幅值的相关系数向量uU，此时i＝m-N，q＝m；初始化传感矩阵列号索引集合为空集，初始化迭代次数n＝1；
5)分别选取电压相角的相关系数向量uθ和电压幅值的相关系数向量uU中最大的z个数值，其中z＝dmax+1，将z个数值对应的传感矩阵中的列号索引构成集合Ωn，更新传感矩阵列号索引集合；
6)利用传感矩阵列号索引集合求最小二乘解，更新残差向量；
7)若更新后的残差向量的2范数R小于残差阈值ε，则进入步骤11)；否则：如果迭代次数n＝1，分别选取电压相角的相关系数向量uθ和电压幅值的相关系数向量uU中最大的2z个值，将2z个值对应的传感矩阵中的索引构成集合Ωn，更新列号索引集合，迭代次数n＝n+1，返回步骤6)；如果迭代次数n＝2，则迭代次数n＝n+1，进入步骤8)；
8)分别计算更新后的残差向量与传感矩阵各列的相关性，得到残差相关系数向量u，选取残差相关系数向量u中最大的2z列，将2z列对应的传感矩阵中的索引构成集合Ωn，更新传感矩阵列号索引集合；
9)利用传感矩阵列号索引集合求最小二乘解，选取最小二乘解中绝对值最大的4z项，将4z项对应的传感矩阵中的索引更新集合Ωn，重新构造列号索引集合，再次更新残差向量；
10)若再次更新的残差向量的2范数R小于残差阈值ε或者迭代次数超过设定的最大迭代次数M，则进入步骤11)；否则n＝n+1，返回步骤8)；
11)输出最小二乘估计结果，根据传感矩阵列号索引集合恢复出2N维向量作为雅可比矩阵第m行的估计结果，m＝m+1；若m大于2N，停止迭代，输出雅可比矩阵估计结果，否则，返回步骤4)。
2.根据权利要求1所述的基于改进匹配追踪的智能配电网潮流雅可比矩阵估计方法，其特征在于，步骤3)中
(1)所述的C组有功功率、无功功率、电压幅值和电压相角的变化向量表示为：
ΔPi[k]＝Pi(k)-Pi(0)、ΔQi[k]＝Qi(k)-Qi(0)、ΔVi[k]＝Vi(k)-Vi(0)和Δθi[k]＝θi(k)-θi(0)，k＝1,2,…,C，Pi(0)、Qi(0)、θi(0)、Vi(0)分别表示节点i当前时刻的有功功率、无功功率、电压相角和电压幅值的量测值；Pi(k)、Qi(k)、θi(k)、Vi(k)分别表示节点i第k组的历史量测值；
(2)所述构造传感矩阵A如下：
<mrow>
<mover>
<mi>A</mi>
<mo>~</mo>
</mover>
<mo>=</mo>
<mo>&lsqb;</mo>
<msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&Delta;&theta;</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>...</mo>
<mi>N</mi>
</mrow>
</msub>
<msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&Delta;V</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mi>N</mi>
</mrow>
</msub>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mo>|</mo>
<msub>
<mover>
<mi>A</mi>
<mo>~</mo>
</mover>
<mi>q</mi>
</msub>
<mo>|</mo>
<mo>=</mo>
<msqrt>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>p</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>C</mi>
</munderover>
<msubsup>
<mover>
<mi>A</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>p</mi>
<mo>,</mo>
<mi>q</mi>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</msqrt>
</mrow>
<mrow>
<msub>
<mi>A</mi>
<mrow>
<mi>p</mi>
<mo>,</mo>
<mi>q</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mover>
<mi>A</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>p</mi>
<mo>,</mo>
<mi>q</mi>
</mrow>
</msub>
<mrow>
<mo>|</mo>
<msub>
<mover>
<mi>A</mi>
<mo>~</mo>
</mover>
<mi>q</mi>
</msub>
<mo>|</mo>
</mrow>
</mfrac>
</mrow>
式中，表示由电压相角和电压幅值量测变化向量构成的矩阵，Δθi＝[Δθi[1],…,Δθi[C]]T表示节点i的C组电压相角量测变化量组成的列向量，ΔVi＝[ΔVi[1],…,ΔVi[C]]T表示节点i的C组电压幅值量测变化量组成的列向量；表示中的元素，Ap,q表示传感矩阵A中第p行第q列的元素。
3.根据权利要求1所述的基于改进匹配追踪的智能配电网潮流雅可比矩阵估计方法，其特征在于，步骤4)中
(1)所述的选取节点i的有功功率量测变化向量初始化残差向量表示为：
r0＝ΔPi
式中，r0表示初始的残差向量，ΔPi＝[ΔPi[1],…,ΔPi[C]]T表示节点i的C组有功功率量测变化量组成的列向量；
(2)用节点i的无功功率量测变化向量初始化残差向量表示为：
r0＝ΔQi
式中，ΔQi＝[ΔQi[1],…,ΔQi[C]]T表示节点i的C组无功功率量测变化量组成的列向量。
(3)电压相角的相关系数向量uθ的计算方法为：
若1≤m≤N：
uθ＝abs[ATAq]
若N＜m≤2N：
uθ＝abs[ATAq-N]
式中，abs[·]表示取绝对值运算，A为传感矩阵，Aq和Aq-N分别表示矩阵A的第q和q-N列；
(4)电压幅值的相关系数向量uU的计算方法为：
若1≤m≤N：
uU＝abs[ATAq+N]
若N＜m≤2N：
uU＝abs[ATAq]
4.根据权利要求1所述的基于改进匹配追踪的智能配电网潮流雅可比矩阵估计方法，其特征在于，步骤5)所述的更新传感矩阵列号索引集合表示为：
Λn＝Λn-1∪Ωn
式中，Λn为第n次迭代的列号索引集合，当n＝1时，Λn-1表示初始的列号索引集合，Ωn表示新选入的索引构成的集合。
5.根据权利要求1所述的基于改进匹配追踪的智能配电网潮流雅可比矩阵估计方法，其特征在于，步骤6)包括：
(1)最小二乘解表示为：
<mrow>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mi>n</mi>
</msub>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>A</mi>
<msub>
<mi>&Lambda;</mi>
<mi>n</mi>
</msub>
<mi>T</mi>
</msubsup>
<msub>
<mi>A</mi>
<msub>
<mi>&Lambda;</mi>
<mi>n</mi>
</msub>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msubsup>
<mi>A</mi>
<msub>
<mi>&Lambda;</mi>
<mi>n</mi>
</msub>
<mi>T</mi>
</msubsup>
<msub>
<mi>r</mi>
<mn>0</mn>
</msub>
</mrow>
式中，表示第n次迭代时的最小二乘解，r0为初始的残差向量，表示第n次迭代的列号索引集合Λn中元素对应的传感矩阵各列构成的矩阵；
(2)更新后的残差向量表示为：
<mrow>
<msub>
<mi>r</mi>
<mi>n</mi>
</msub>
<mo>=</mo>
<msub>
<mi>r</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<msub>
<mi>A</mi>
<msub>
<mi>&Lambda;</mi>
<mi>n</mi>
</msub>
</msub>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mi>n</mi>
</msub>
</mrow>
式中，rn为第n次迭代的残差向量。
6.根据权利要求1所述的基于改进匹配追踪的智能配电网潮流雅可比矩阵估计方法，其特征在于，步骤8)所述的残差相关系数向量u计算方式为：
u＝abs[ATrn-1]
式中，u为相关系数向量，abs[·]表示取绝对值运算，rn-1表示第n-1次迭代时的残差向量，当n＝1时，rn-1表示初始残差向量，A为传感矩阵。
7.根据权利要求1所述的基于改进匹配追踪的智能配电网潮流雅可比矩阵估计方法，其特征在于，步骤9)中
(1)所述的重新构造列号索引集合为
Λn＝Ωn
(2)所述的再次更新残差向量表示为：
<mrow>
<msub>
<mi>r</mi>
<mi>n</mi>
</msub>
<mo>=</mo>
<msub>
<mi>r</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<msub>
<mi>A</mi>
<msub>
<mi>&Lambda;</mi>
<mi>n</mi>
</msub>
</msub>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>A</mi>
<msub>
<mi>&Lambda;</mi>
<mi>n</mi>
</msub>
<mi>T</mi>
</msubsup>
<msub>
<mi>A</mi>
<msub>
<mi>&Lambda;</mi>
<mi>n</mi>
</msub>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msubsup>
<mi>A</mi>
<msub>
<mi>&Lambda;</mi>
<mi>n</mi>
</msub>
<mi>T</mi>
</msubsup>
<msub>
<mi>r</mi>
<mn>0</mn>
</msub>
</mrow>
式中，Λn为第n次迭代的列号索引集合，Ωn表示新选入的索引构成的集合，表示集合Λn中元素对应的传感矩阵各列构成的矩阵，rn为第n次迭代的残差向量，r0为初始的残差向量。
8.根据权利要求1所述的基于改进匹配追踪的智能配电网潮流雅可比矩阵估计方法，其特征在于，步骤11)所述的根据传感矩阵列号索引集合恢复出2N维向量作为雅可比矩阵第m行的估计结果表示为：
<mrow>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<msub>
<mi>&Lambda;</mi>
<mi>n</mi>
</msub>
</msub>
<mo>=</mo>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mi>n</mi>
</msub>
</mrow>
<mrow>
<msub>
<mover>
<mi>Y</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>m</mi>
<mo>,</mo>
<mi>q</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mi>g</mi>
</msub>
<mrow>
<mo>|</mo>
<msub>
<mover>
<mi>A</mi>
<mo>~</mo>
</mover>
<mi>q</mi>
</msub>
<mo>|</mo>
</mrow>
</mfrac>
</mrow>
式中，表示第n次迭代时的最小二乘解，表示2N维恢复向量，表示的第g个元素，表示列号索引集合Λn中各元素所对应的向量的元素组成的向量，对于其他不在列号索引集合中的元素取值为0；矩阵第q列的模值， 表示雅可比矩阵Y的第m行第q个元素的估计解，其中，g＝q。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711118202.6A CN107749627B (zh) | 2017-11-13 | 2017-11-13 | 基于改进匹配追踪的智能配电网潮流雅可比矩阵估计方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711118202.6A CN107749627B (zh) | 2017-11-13 | 2017-11-13 | 基于改进匹配追踪的智能配电网潮流雅可比矩阵估计方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107749627A true CN107749627A (zh) | 2018-03-02 |
| CN107749627B CN107749627B (zh) | 2020-09-11 |
Family
ID=61250724
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711118202.6A Active CN107749627B (zh) | 2017-11-13 | 2017-11-13 | 基于改进匹配追踪的智能配电网潮流雅可比矩阵估计方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107749627B (zh) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110221168A (zh) * | 2019-05-15 | 2019-09-10 | 广西电网有限责任公司电力科学研究院 | 一种主导谐波源定位和谐波污染传播路径追踪方法 |
| CN111766445A (zh) * | 2020-07-11 | 2020-10-13 | 天津大学 | 一种配电网同步相量相位量测的连续化转换和恢复方法 |
| CN111865325A (zh) * | 2020-07-10 | 2020-10-30 | 山东云海国创云计算装备产业创新中心有限公司 | 一种压缩感知信号重构方法、装置及相关设备 |
| CN113221307A (zh) * | 2021-06-11 | 2021-08-06 | 南方电网科学研究院有限责任公司 | 基于压缩感知算法的配电网暂态运行拓扑识别方法及装置 |
| CN116360352A (zh) * | 2022-12-02 | 2023-06-30 | 山东和信智能科技有限公司 | 一种发电厂用智能控制方法及系统 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003134666A (ja) * | 2001-10-16 | 2003-05-09 | Takeo Kawamura | 非線形計画法による最適潮流計算システム |
| CN106340875A (zh) * | 2016-10-26 | 2017-01-18 | 珠海许继芝电网自动化有限公司 | 一种配电网多相状态估计方法 |
| CN107194131A (zh) * | 2017-07-10 | 2017-09-22 | 大连海事大学 | 基于Matlab稀疏矩阵的极坐标牛顿法潮流计算方法 |
-
2017
- 2017-11-13 CN CN201711118202.6A patent/CN107749627B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003134666A (ja) * | 2001-10-16 | 2003-05-09 | Takeo Kawamura | 非線形計画法による最適潮流計算システム |
| CN106340875A (zh) * | 2016-10-26 | 2017-01-18 | 珠海许继芝电网自动化有限公司 | 一种配电网多相状态估计方法 |
| CN107194131A (zh) * | 2017-07-10 | 2017-09-22 | 大连海事大学 | 基于Matlab稀疏矩阵的极坐标牛顿法潮流计算方法 |
Non-Patent Citations (1)
| Title |
|---|
| 赵冠琨: "直流配电网潮流计算模型及算法", 《电力建设》 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110221168A (zh) * | 2019-05-15 | 2019-09-10 | 广西电网有限责任公司电力科学研究院 | 一种主导谐波源定位和谐波污染传播路径追踪方法 |
| CN110221168B (zh) * | 2019-05-15 | 2021-03-09 | 广西电网有限责任公司电力科学研究院 | 一种主导谐波源定位和谐波污染传播路径追踪方法 |
| CN111865325A (zh) * | 2020-07-10 | 2020-10-30 | 山东云海国创云计算装备产业创新中心有限公司 | 一种压缩感知信号重构方法、装置及相关设备 |
| CN111766445A (zh) * | 2020-07-11 | 2020-10-13 | 天津大学 | 一种配电网同步相量相位量测的连续化转换和恢复方法 |
| CN111766445B (zh) * | 2020-07-11 | 2022-09-23 | 天津大学 | 一种配电网同步相量相位量测的连续化转换和恢复方法 |
| CN113221307A (zh) * | 2021-06-11 | 2021-08-06 | 南方电网科学研究院有限责任公司 | 基于压缩感知算法的配电网暂态运行拓扑识别方法及装置 |
| CN116360352A (zh) * | 2022-12-02 | 2023-06-30 | 山东和信智能科技有限公司 | 一种发电厂用智能控制方法及系统 |
| CN116360352B (zh) * | 2022-12-02 | 2024-04-02 | 山东和信智能科技有限公司 | 一种发电厂用智能控制方法及系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107749627B (zh) | 2020-09-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107749627B (zh) | 基于改进匹配追踪的智能配电网潮流雅可比矩阵估计方法 | |
| CN108199375B (zh) | 基于同步相量量测的智能配电网拓扑辨识方法 | |
| CN107742885B (zh) | 基于正则匹配追踪的配电网电压功率灵敏度估计方法 | |
| CN107133406B (zh) | 一种电力系统静态电压稳定域边界的快速搜索方法 | |
| CN108462181B (zh) | 考虑稀疏性的智能配电网潮流雅可比矩阵鲁棒估计方法 | |
| CN107658881A (zh) | 基于戴维南等值方法的电压稳定临界点判断方法 | |
| CN107204617B (zh) | 基于线性规划的直角坐标形式的区间潮流计算方法 | |
| CN115932702B (zh) | 基于虚拟标准器的电压互感器在线运行校准方法及装置 | |
| CN109088407B (zh) | 基于深度信念网络伪量测建模的配电网状态估计方法 | |
| CN109754013B (zh) | 一种基于无迹卡尔曼滤波的电力系统混合量测融合方法 | |
| CN107991647A (zh) | 基于WiFi和UWB的室内定位方法及系统、存储介质 | |
| CN108448585B (zh) | 一种基于数据驱动的电网潮流方程线性化求解方法 | |
| CN105548718A (zh) | 一种基于混合整体最小二乘法的系统谐波阻抗计算方法 | |
| CN106443246A (zh) | 基于pmu量测数据的小干扰稳定参数的在线辨识方法 | |
| CN104102836A (zh) | 一种电力系统快速抗差状态估计方法 | |
| CN107634516A (zh) | 一种基于灰色‑马尔可夫链的配网状态估计方法 | |
| CN102280877B (zh) | 一种多量测断面的电力系统不良支路参数辨识方法 | |
| CN114168574A (zh) | 一种面向工业负荷的数据缺失处理方法及装置 | |
| CN116706921A (zh) | 基于hhl算法的量子牛顿-拉夫逊法潮流计算方法和系统 | |
| CN103312297B (zh) | 一种迭代扩展增量卡尔曼滤波方法 | |
| CN108649585B (zh) | 一种电力系统静态电压稳定域边界快速搜索的直接法 | |
| CN113158135A (zh) | 一种含噪声的暂降源定位数据缺失值估计方法 | |
| CN110703038B (zh) | 一种适用于风机接入配电网的谐波阻抗估算方法 | |
| CN111639463A (zh) | 一种基于XGBoost算法的电力系统扰动后频率特征预测方法 | |
| CN113553538B (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 | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20210831 Address after: No.69 Feitian Avenue, Airport Economic Development Zone, Jiangning District, Nanjing City, Jiangsu Province Patentee after: Nanjing Shoufeng Smart Power Research Institute Co.,Ltd. Address before: 300192 92 Weijin Road, Nankai District, Tianjin Patentee before: Tianjin University |