CN113297531B - 一种完美相位估计下直流潮流方程的求解方法 - Google Patents
一种完美相位估计下直流潮流方程的求解方法 Download PDFInfo
- Publication number
- CN113297531B CN113297531B CN202110430057.5A CN202110430057A CN113297531B CN 113297531 B CN113297531 B CN 113297531B CN 202110430057 A CN202110430057 A CN 202110430057A CN 113297531 B CN113297531 B CN 113297531B
- Authority
- CN
- China
- Prior art keywords
- direct current
- flow equation
- power flow
- phase estimation
- current 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
Images
Classifications
-
- 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/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- 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
- H02J1/00—Circuit arrangements for dc mains or dc 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Power Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Operations Research (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Optical Modulation, Optical Deflection, Nonlinear Optics, Optical Demodulation, Optical Logic Elements (AREA)
- Complex Calculations (AREA)
Abstract
本发明提出一种完美相位估计下直流潮流方程的求解方法,该方法根据前轮迭代的量子线路的测量结果设计下一轮迭代的量子线路,在迭代结束后对底部量子寄存器执行测量,最后从经典寄存器中提取信息求解直流潮流方程。本发明通过迭代的方式节省了求解直流潮流方程过程中量子比特资源的消耗,可在完美相位估计情况下基于小规模量子系统求解大规模直流潮流方程,避免因量子比特资源不足而无法利用量子程序求解直流潮流方程的问题,提供了一种基于量子比特资源受限的现有硬件平台求解直流潮流方程的方法。
Description
技术领域
本发明属于电力系统领域,涉及一种完美相位估计下直流潮流方程的求解方法。
背景技术
电力系统潮流计算是电力系统最基本的计算,也是最重要的计算。所谓潮流计算,就是已知电网的接线方式、参数及运行条件,计算电力系统稳态运行时各母线电压、各支路电流与功率及网损。直流潮流模型是把潮流计算中非线性模型经过简化形成线性模型,从而使计算变得简单,在一些诸如过负荷校验计算等对求解精度要求不高的场景下具有巨大优势。
美国Denver大学的Eskandarpour等在《Quantum Computing Solution of DCPower Flow》一文中提出把潮流方程线性化后,使用Harrow,Hassidim以及Lloyd于2009年提出的HHL算法实现计算加速并进行了简单的算例演示,该算法可将求解线性方程组的时间复杂度降低至O(log(N)s2κ2/ε),与当前最好的经典算法相比达到了指数级别的加速。
由于HHL算法将求得的直流潮流方程的解存储在量子寄存器中,通过测量提取出的实际上是归一化的解的信息,若要得到实际解,还需要通过回代的方式求出归一化常数,因此发展能直接得到非归一化的解的算法,具有显著的意义。
同时HHL算法中所需量子比特数目随着直流潮流方程规模增大而增多,其大规模实际应用要建立在已有通用量子计算机的前提上,而当前处于有噪声中程量子时代,量子比特资源是受限的(数目在50~100个,同时存在量子噪音),故在此前提下发展所需量子比特数较少的直流潮流方程求解方法具有重要的现实意义和应用前景。
发明内容
本发明所要解决的技术问题是现有求解直流潮流方程的量子算法需要占用大量的量子比特资源,因此提供一种基于完美迭代相位估计求解直流潮流方程的方法。直流潮流方程为P=Bθ,其中P为有功功率,B为电纳矩阵,θ为要求解的相角。当B的特征值可以被有限位二进制数完整编码即相位可以被完美估计时,通过对迭代过程中量子系统的状态执行测量并进行信息提取及处理,得到直流潮流方程的解θ。
本发明为解决上述技术问题采用以下技术方案:
1.一种完美相位估计下直流潮流方程的求解方法,其特征在于,该方法为基于完美迭代相位估计的经典-量子混合算法,通过经典计算机处理完美迭代相位估计的量子线路的测量结果求得直流潮流方程的解,主要步骤为:
(1)提出根据量子系统规模设计迭代相位估计路线,分为单辅助量子比特路线,N辅助量子比特路线以及任意辅助量子比特路线;
(2)提取直流潮流方程中电纳矩阵B的特征值λj及相角θ在电纳矩阵B各特征向量上的投影|βj|;
(3)提取电纳矩阵B特征向量各元素取绝对值后的|uj|;
(4)进行βjuj正负号校准;
3.完美相位估计下直流潮流方程的求解方法,其特征在于,所述步骤(3)中|uj|的提取,通过对系统状态进行后选择操作,当m轮迭代的测量结果为λj时,对底部寄存器执行测量,对经典寄存器存储的信息进行统计,提取出|uj|。
附图说明
图1是基于单辅助量子比特混合算法的量子线路图。
图2是基于N辅助量子比特混合算法的量子线路图。
图3是基于任意辅助量子比特混合算法的量子线路图。
具体实施方式
本发明根据量子系统规模设计三种迭代相位估计路线。当量子比特资源匮乏时,如实施例一所示,设计基于单辅助量子比特的迭代相位估计路线;当顶部寄存器可分配量子比特数大于N时,如实施例二所示,设计基于N辅助量子比特的迭代相位估计路线;而实施例三给出基于任意辅助量子比特的迭代相位估计路线,具有更高灵活性。
1.实施例一。参见图1,图1是基于单辅助量子比特混合算法的量子线路图,求解直流潮流方程的流程如下:
(1)如图1所示,执行第一轮迭代,此后将量子寄存器初始化,其中其中B即直流潮流方程P=Bθ中的电纳矩阵B,下同。在此轮迭代中,Rz门的旋转参数ωm=-2π(0.0)2,其中(0.0)2表示二进制小数0.0,下同。将第一轮迭代中测量结果为0的情况记为Cm_0,测量结果为1的情况记为Cm_1,其中m表示电纳矩阵特征值从左往右第m位。将Cm_0发生的概率记为Pm_0,Cm_1发生的概率记为Pm_1;
(2)执行第二轮迭代,迭代后将量子寄存器初始化。在该轮迭代中,根据以下三种情形来设计本轮迭代中Rz门的旋转参数:情况一,Pm_0=1,Pm_1=0,将ω(m-1)=-2π(0.00)2设计为下一轮迭代的旋转参数;情况二,Pm_0=0,Pm_1=1,将ω(m-1)=-2π(0.01)2设计为下一轮迭代的旋转参数;情况三,Pm_1≠0,Pm_1≠0,在此处进行标记,并将当前实验记为实验1。为实验1中的下一轮迭代设计旋转参数ω(m-1)=-2π(0.00)2,并设计实验2第二轮迭代中旋转参数ω(m-1)=-2π(0.01)2。记第一次与第二次迭代的测量结果为00的情况为C(m-1)m_00,并将其对应的概率表示为P(m-1)m_00。记N维电纳矩阵的N个特征值λj=φj=0.φj1φj2…φjm,经过第二轮迭代,得到和
(4)在m次迭代的测量结果为λj=φj=(0.φj1φj2…φj(m-1)φjm)2时,对底部寄存器执行测量,提取出电纳矩阵B的特征向量uj各元素取绝对值后的|uj|;
由于|βjuj|=|βj|×|uj|,具体数值已知,通过对其系数正负号执行遍历进行验证,得到βjuj的正负号信息。
2.实施例二。参见图2,图2是基于N辅助量子比特混合算法的量子线路图,其中顶部寄存器中含有anc_q[0]至anc_q[N-1]的N个量子比特,求解直流潮流方程的流程如下:
(1)如图2所示,执行第一轮迭代。在该轮迭代中,仅使用顶部寄存器中的一个辅助量子比特,旋转参数设置为ωm_0=-2π(0.0)2,提取信息Cm_0,Cm_1,Pm_0和Pm_1;
(2)根据以下三种情况设计第二轮迭代的量子线路图以及旋转参数:情况一,Pm_0=1,Pm_1=0,第二轮迭代中,顶部寄存器中仅使用anc_q[0]一个量子比特,旋转参数设计为ω(m-1)_0=-2π(0.00)2,以其作为anc_q[0]上Rz门的旋转参数;情况二,Pm_0=0,Pm_1=1,第二轮迭代中,顶部寄存器中仅使用anc_q[0]一个量子比特,旋转参数设计为ω(m-1)_0=-2π(0.01)2,以其作为anc_q[0]上Rz门的旋转参数;情况三,Pm_1≠0,Pm_1≠0,第二轮迭代中,顶部寄存器中使用anc_q[0]和anc_q[1]两个量子比特,将旋转参数设计为ω(m-1)_0=-2π(0.00)2和ω(m-1)_1=-2π(0.01)2,将ω(m-1)_0作为anc_q[0]上Rz门的旋转参数,将ω(m-1)_1作为第二轮迭代中anc_q[1]上Rz门的旋转参数。通过此轮迭代,提取出信息P(m-1)m_00,P(m-1)m_10,P(m-1)m_01和P(m-1)m_11。
(3)以此类推,如图2所示,假设第一个测量结果的分歧产生在第x1轮迭代,则在后一轮迭代将anc_q[1]投入使用,假设第N-1个测量结果的分歧产生在第xN-1轮迭代,则第xN-1轮迭代后,顶部寄存器中的N个量子比特投入使用。经过m轮迭代后,提取信息λj=φj=0.φj1φj2…φj(m-1)φjm以及其对应的概率随之得到|βj|。
(4)采用实施例一中步骤(4)给出的|uj|提取方案,提取|uj|;
(5)采用实施例一中步骤(5)给出的βjuj正负号校准方案,提取βjuj整体信息;
(6)采用实施例一中步骤(6)给出的求解方案,求直流潮流方程的解。
3.实施例三。参见图3,图3是基于任意辅助量子比特混合算法的量子线路图,其中顶部寄存器中含有anc_q[0]至anc_q[n-1]的n个量子比特,求解直流潮流方程的流程如下:
(1)第一轮迭代中,对旋转参数可能的方案执行遍历,每个实验旋转参数的设置如下表所示;
(2)第二轮迭代中,记ω=-2π(nbin)2,其中nbin为二进制小数,以不同的φj(m-n+1)φj(m-n+2)…φjm为nbin从右往左的n项,从右往左第n+1项起按照第一轮迭代中遍历的思想设置旋转参数;当m可被n整除时,经过m/n轮迭代,提取所有特征值信息λj=φj=0.φj1φj2…φj(m-1)φjm以及开方后得到|βj|。
(3)采用实施例一中步骤(4)给出的|uj|提取方案,提取|uj|;
(4)采用实施例一中步骤(5)给出的βjuj正负号校准方案,提取βjuj整体信息;
(5)采用实施例一中步骤(6)给出的求解方案,求直流潮流方程的解。
Claims (1)
1.一种完美相位估计下直流潮流方程的求解方法,其特征在于,该方法为基于完美迭代相位估计的经典-量子混合算法,通过经典计算机处理完美迭代相位估计的量子线路的测量结果求得直流潮流方程的解;直流潮流方程为P=Bθ;其中P为有功功率;B为电纳矩阵;θ为要求解的相角;提出根据量子系统规模设计迭代相位估计线路;量子比特资源匮乏时,设计基于单辅助量子比特的迭代相位估计线路;顶部寄存器分配量子比特数大于N时,设计基于N辅助量子比特的迭代相位估计线路;具有更高灵活性的基于任意辅助量子比特的迭代相位估计线路,求解直流潮流方程的流程如下:
(2)提取电纳矩阵B特征向量各元素取绝对值后的|uj|:在m次迭代的测量结果为λj=φj=(0.φj1φj2…φj(m-1)φjm)2时,对底部寄存器执行测量,提取出电纳矩阵B的特征向量uj各元素取绝对值后的|uj|;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110430057.5A CN113297531B (zh) | 2021-04-21 | 2021-04-21 | 一种完美相位估计下直流潮流方程的求解方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110430057.5A CN113297531B (zh) | 2021-04-21 | 2021-04-21 | 一种完美相位估计下直流潮流方程的求解方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113297531A CN113297531A (zh) | 2021-08-24 |
CN113297531B true CN113297531B (zh) | 2022-09-02 |
Family
ID=77320037
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110430057.5A Active CN113297531B (zh) | 2021-04-21 | 2021-04-21 | 一种完美相位估计下直流潮流方程的求解方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113297531B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117977600B (zh) * | 2024-04-02 | 2024-06-18 | 合肥工业大学 | 电力系统多潮流集成并行量子计算方法、系统和存储介质 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105305442A (zh) * | 2015-11-30 | 2016-02-03 | 河海大学常州校区 | 基于量子遗传算法的多目标智能配电网自愈恢复方法 |
CN107317338A (zh) * | 2017-08-30 | 2017-11-03 | 广东工业大学 | 一种电力系统的最优潮流计算方法及装置 |
CN107437811A (zh) * | 2017-09-13 | 2017-12-05 | 广西大学 | 电力系统暂态稳定约束最优潮流并行计算方法 |
CN108400592A (zh) * | 2018-03-19 | 2018-08-14 | 国网江西省电力有限公司电力科学研究院 | 一种计及潮流约束的配电网区间状态估计算法 |
CN109767007A (zh) * | 2018-12-10 | 2019-05-17 | 东南大学 | 一种基于量子计算的最小均方误差检测方法 |
CN111095307A (zh) * | 2017-09-22 | 2020-05-01 | 国际商业机器公司 | 用于量子计算机器的硬件高效的变分量子本征值求解器 |
CN112232512A (zh) * | 2020-09-08 | 2021-01-15 | 中国人民解放军战略支援部队信息工程大学 | 量子计算模拟平台及线性方程组量子求解模拟方法、系统 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20200342330A1 (en) * | 2019-04-26 | 2020-10-29 | International Business Machines Corporation | Mixed-binary constrained optimization on quantum computers |
-
2021
- 2021-04-21 CN CN202110430057.5A patent/CN113297531B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105305442A (zh) * | 2015-11-30 | 2016-02-03 | 河海大学常州校区 | 基于量子遗传算法的多目标智能配电网自愈恢复方法 |
CN107317338A (zh) * | 2017-08-30 | 2017-11-03 | 广东工业大学 | 一种电力系统的最优潮流计算方法及装置 |
CN107437811A (zh) * | 2017-09-13 | 2017-12-05 | 广西大学 | 电力系统暂态稳定约束最优潮流并行计算方法 |
CN111095307A (zh) * | 2017-09-22 | 2020-05-01 | 国际商业机器公司 | 用于量子计算机器的硬件高效的变分量子本征值求解器 |
CN108400592A (zh) * | 2018-03-19 | 2018-08-14 | 国网江西省电力有限公司电力科学研究院 | 一种计及潮流约束的配电网区间状态估计算法 |
CN109767007A (zh) * | 2018-12-10 | 2019-05-17 | 东南大学 | 一种基于量子计算的最小均方误差检测方法 |
CN112232512A (zh) * | 2020-09-08 | 2021-01-15 | 中国人民解放军战略支援部队信息工程大学 | 量子计算模拟平台及线性方程组量子求解模拟方法、系统 |
Non-Patent Citations (7)
Title |
---|
Arbitrary accuracy iterative phase estimation algorithm as a two qubit benchmark;M. Dobsicek等;《arXiv:quant-ph/0610214v3》;20070712;1-4 * |
Blaga N. Todorova ; René Steijl.Quantum algorithm for the collisionless Boltzmann equation.《Journal of Computational Physics》.2020,1-10. * |
Hybrid algorithms to solve linear systems of equations with limited qubit resources;Fang Gao等;《arXiv》;20210629;1-22 * |
Quantum Computing Solution of DC Power Flow;Rozhin Eskandarpour;《arXiv》;20201006;1-19 * |
Quantum Power Flow;Fei Feng等;《arXiv:2104.04888v1》;20210411;1-4 * |
Solving DC Power Flow Problems Using Quantum and Hybrid algorithms;Fang Gao等;《arXiv》;20220113;1-17 * |
基于核自旋的量子计算实验研究;李兆凯;《中国博士学位论文全文数据库 (基础科学辑)》;20141015(第10期);A005-7 * |
Also Published As
Publication number | Publication date |
---|---|
CN113297531A (zh) | 2021-08-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Kerenidis et al. | q-means: A quantum algorithm for unsupervised machine learning | |
WO2022257316A1 (zh) | 量子体系基态能量估计方法及系统 | |
Cortese et al. | Loading classical data into a quantum computer | |
Cabral et al. | Multivariate mixture modeling using skew-normal independent distributions | |
Liang et al. | An equivalent measure of partial correlation coefficients for high-dimensional gaussian graphical models | |
Escanciano et al. | Uniform convergence of weighted sums of non and semiparametric residuals for estimation and testing | |
CN107480694B (zh) | 基于Spark平台采用两次评价的加权选择集成三支聚类方法 | |
US20170061279A1 (en) | Updating an artificial neural network using flexible fixed point representation | |
CN113592093B (zh) | 量子态制备电路生成方法、装置、量子操作芯片及设备 | |
Butka et al. | On equivalence of conceptual scaling and generalized one-sided concept lattices | |
Lavielle et al. | An improved SAEM algorithm for maximum likelihood estimation in mixtures of non linear mixed effects models | |
Liang et al. | Alternating iterative methods for solving tensor equations with applications | |
Li et al. | Application of distributed semi-quantum computing model in phase estimation | |
CN115577791B (zh) | 基于量子系统的信息处理方法及装置 | |
CN113297531B (zh) | 一种完美相位估计下直流潮流方程的求解方法 | |
Huang et al. | Parameter estimation of fractional-order chaotic systems by using quantum parallel particle swarm optimization algorithm | |
Zhou et al. | Quantum image scaling based on bilinear interpolation with decimals scaling ratio | |
Krämer et al. | Splitting methods for nonlinear Dirac equations with Thirring type interaction in the nonrelativistic limit regime | |
Moerkotte et al. | Proximal operator of quotient functions with application to a feasibility problem in query optimization | |
Wang et al. | Hole–particle correspondence in CI calculations | |
Muñoz-Coreas et al. | T-count optimized quantum circuits for bilinear interpolation | |
CN115577792A (zh) | 基于量子系统的信息处理方法及装置 | |
EP3735659A1 (en) | Optimization of learning network in equivalent class space | |
JP2012079225A (ja) | 協調フィルタリング処理方法およびプログラム | |
Wen et al. | Sparse solution of nonnegative least squares problems with applications in the construction of probabilistic Boolean networks |
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 |