一种非完美相位估计下直流潮流方程的求解方法与流程

1.本发明属于电力系统领域，涉及一种非完美相位估计下直流潮流方程的求解方法。


背景技术：

2.电力系统潮流的计算以及分析是电力系统运行和规划工作的基础。而所谓潮流计算，即在已知电网的接线方式、参数运行条件的条件下检验所提出的网络规划方案能否满足各种运行方式的要求。而直流潮流法就是把非线性潮流问题简化为线性直流电路问题，具有精度稍差，但是方便快捷的特点，非常适合处理断线分析和线性规划处理的优化问题，具有广泛的应用。随着电力系统规模的不断扩大，探索更高效的直流潮流方程求解方法变得日益重要。
3.eskandarpour等在《quantum computing solution of dc power flow》提出使用hhl量子算法求解直流潮流方程，相较于当前最好的经典算法可以达到指数级别的加速。而hhl算法仅能得到归一化的解的信息，且当电纳矩阵的特征值无法被有限位二进制数完整编码即相位不可以被完美估计时，如nielsen在《quantum computation and quantum information》一书中所述，相位估计存在失败率，而相位估计的失败会导致使用hhl算法求得的直流潮流的解与实际解产生偏差，故在求解直流潮流方程的过程中，对量子比特资源的较高占用使该算法在量子比特资源受限的条件下无法得到有效应用。
4.进而提出的基于完美迭代相位估计求解直流潮流的方法，可以求得未归一化的解，且可以有效节约量子比特资源，但其目前仅可基于完美相位估计对直流潮流方程进行求解，具有较大的局限性。故发展出非完美相位估计下直流潮流方程得求解方法具有重要的现实意义和应用前景。


技术实现要素：

5.本发明所要解决的技术问题是现有求解直流潮流方程的量子方法需要占用大量的量子比特资源，而基于完美迭代相位估计算法求解直流潮流方程的方法目前仅可用于电纳矩阵b的特征值可被二进制数完美表示的情形，因此提供一种基于非完美相位估计求解直流潮流方程的方法。直流潮流方程为p＝bθ，其中p为节点注入有功功率，b为n维电纳矩阵，θ为要求解的相角，记λ为电纳矩阵b的特征值，λ
j
为电纳矩阵b的第j个特征值，u
j
为电纳矩阵b的特征值λ
j
对应的单位特征向量，u
jp
为单位特征向量u
j
的第p个元素，p
p
为节点注入有功功率p的第p个元素，β
j
为相角θ在电纳矩阵b各单位特征向量上的投影。该方法的使用需要顶部寄存器、底部寄存器两个量子寄存器以及一个经典寄存器，其中顶部寄存器中量子比特分为决定单个相位估计模块精度的精度比特和用以提高相位估计成功率的冗余比特。在单个相位估计模块中分配用于提高相位估计成功率的冗余比特，再通过对多个相位估计模块中顶部寄存器中精度比特的状态执行测量并进行信息提取及处理，得到近似解再根据此近似解更新待输出的直流潮流方程的解θ
*
，使其趋近于实际解。
6.本发明为解决上述技术问题采用以下技术方案：
7.1.一种非完美相位估计下直流潮流方程的求解方法，其特征在于，该方法为基于多个相位估计模块求解方程组，再通过外循环减小误差的经典
‑
量子混合方法，通过经典计算机处理多个非完美相位估计模块的量子线路的测量结果求得直流潮流方程的解，主要步骤为：
8.(1)构造直流潮流方程p＝bθ，初始化待输出的解θ
*
＝0；
9.(2)分配顶部量子寄存器、底部量子寄存器以及经典寄存器并设计相位估计模块的数量，根据直流潮流方程中电纳矩阵特征值的范围以及量子系统的规模，设计相位估计模块数以及单个相位估计模块顶部量子寄存器中精度比特和冗余比特的数量；
10.(3)映射有功功率p至底部量子寄存器，记为|c
·
p>；
11.(4)基于多个非完美相位估计模块求解方程组，得到近似解
12.(5)更新待输出的解
13.(6)当误差ξ小于误差阈值ξ
max
时，输出直流潮流方程的解θ
*
；反之则更新并重复步骤(3)至步骤(6)。
14.2.非完美相位估计下直流潮流方程的求解方法，其特征在于，所述步骤(4)中基于多个非完美相位估计模块求解方程组，主要步骤为：
15.(1)提取直流潮流方程中电纳矩阵b特征值λ
j
的近似值及|c
·
p>在电纳矩阵b各单位特征向量上的投影|β
j
|的近似值
16.(2)提取电纳矩阵b单位特征向量u
j
各元素取绝对值后的|u
jp
|的近似值
17.(3)通过进行正负号校准；
18.(4)根据步骤(1)和(3)提取的和通过计算得到近似解
附图说明
19.图1是非完美相位估计下直流潮流方程的求解流程图；
20.图2是基于多个非完美相位估计模块求解方程组方法的量子线路图。
具体实施方式
21.本发明根据量子系统的规模、解的精度需求以及相位估计失败率上限来设计相位估计模块的数量以及单个相位估计模块内精度比特与冗余比特的数量。
22.1.实施例一。参见图1，图1是非完美相位估计下直流潮流方程的求解流程图，在电纳矩阵b的特征值无法被有限位二进制数完整表示时，求解直流潮流方程的流程如下：
23.(1)如图1所示，构造直流潮流方程p＝bθ，同时初始化θ
*
＝0，其中p为有功功率，b为电纳矩阵，θ为待求的相角，θ
*
为待输出的解，下同；
24.(2)如图1所示，分配顶部、底部量子寄存器及经典寄存器。当要估计的特征值λ
j
的精度需求为其二进制形式最高位开始的后m位时,带截断误差的近似特征值为使用m/n
accur
个相位估计模块，其中单个相位估计模块中顶部寄存器共分配n
top_q
个量子比特，精度比特的数量为n
accur
，冗余比特的数量为n
redun
；
25.(3)如图1所示，将有功功率p进行归一化处理为c
·
p，进而映射至底部量子寄存器，构造底部量子寄存器初态|c
·
p>，其中c为归一化常数；
26.(4)如图1所示，基于多个非完美相位估计模块求解方程组，得到近似解
27.(5)如图1所示，更新待输出的解
28.(6)如图1所示，判断误差是否超出上限。记误差若误差ξ小于误差上限ξ
max
，输出直流潮流方程的解θ
*
；反之则更新并重复步骤(3)至步骤(6)。
29.2.实施例一步骤(4)中基于多个非完美相位估计模块求解方程组的方法。参见图2，图2是基于多个非完美相位估计模块求解方程组方法的量子线路图，qpe
‑
m表示第m个相位估计模块。其求解方程组的步骤如下：
30.(1)如图2所示，将系统状态初始化为执行qpe
‑
1，其中b为直流潮流方程p＝bθ中的电纳矩阵b，下同。对顶部寄存器中n
accur
个精度比特执行测量并将测量结果存储至经典寄存器，此时经典寄存器以的成功率存储电纳矩阵b特征值λ的高n
accur
位；
31.(2)如图2所示，初始化系统状态，执行qpe
‑
2，对顶部寄存器中n
accur
个精度比特执行测量并存储至经典寄存器后，经典寄存器以的成功率存储电纳矩阵b特征值λ的高2n
accur
位；
32.(3)如图2所示，每次执行相位估计前初始化系统状态，依次执行m/n
accur
个相位估计模块，并完成每个相位估计模块中对精度比特的测量以及测量结果的存储，经典寄存器以的成功率存储电纳矩阵b特征值λ的高m位，即存储电纳矩阵b带截断误差的近似特征值对测量结果进行统计，测量结果为的概率为提取出信息和
33.(4)如图2所示，在m/n
accur
个相位估计模块精度比特的测量结果为时，对底部寄存器执行测量，提取电纳矩阵b单位特征向量u
j
各元素取绝对值后的|u
jp
|的近似值
34.(5)根据联立如下方程组：
[0035][0036]
其中p
p
表示p的第p个元素，由于具体数值已知，通过遍历系数的正负号，得到的正负号信息；
[0037]
(6)通过等式求得解的近似值