CN113962108A - 一种电力系统区间抗差状态估计方法、系统及装置 - Google Patents

一种电力系统区间抗差状态估计方法、系统及装置 Download PDF

Info

Publication number
CN113962108A
CN113962108A CN202111319373.1A CN202111319373A CN113962108A CN 113962108 A CN113962108 A CN 113962108A CN 202111319373 A CN202111319373 A CN 202111319373A CN 113962108 A CN113962108 A CN 113962108A
Authority
CN
China
Prior art keywords
interval
state estimation
solving
state variable
establishing
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.)
Pending
Application number
CN202111319373.1A
Other languages
English (en)
Inventor
曾顺奇
李欣
王斐
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Guangzhou Power Supply Bureau of Guangdong Power Grid Co Ltd
Original Assignee
Guangzhou Power Supply Bureau of Guangdong Power Grid Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Guangzhou Power Supply Bureau of Guangdong Power Grid Co Ltd filed Critical Guangzhou Power Supply Bureau of Guangdong Power Grid Co Ltd
Priority to CN202111319373.1A priority Critical patent/CN113962108A/zh
Publication of CN113962108A publication Critical patent/CN113962108A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Electricity, gas or water supply
    • 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/04—Constraint-based CAD
    • 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]

Abstract

本发明涉及电力系统状态评估技术领域，公开了一种电力系统区间抗差状态估计方法、系统及装置。本发明在对建立的区间约束模型进行求解时，根据得到的状态变量区间值近似解确定初值，并以求解过程结束后输出的状态变量可行域为约束，以量测点在一定置信水平上的相对偏差最小化作为目标函数，建立对应的状态估计模型，进而对状态估计模型进行优化求解，得到状态估计结果；本发明通过取状态变量的可行域为约束，以状态变量可行域设置初值，无需将潮流计算结果作为初值，可有效缩小状态变量的求解范围和区间范围，减少计算量，从而提高计算效率和收敛速度。

Description

一种电力系统区间抗差状态估计方法、系统及装置
技术领域
本发明涉及电力系统状态评估技术领域，尤其涉及一种电力系统区间抗差状态估计方法、系统及装置。
背景技术
状态估计技术可以通过可观测的信息评估电力系统中节点的运行状况，以及电压和功率的状态值，从而为最优潮流、无功优化等应用提供可靠、真实的数据信息。目前，基于加权最小二乘法的状态估计方法被广泛使用，但是这种状态估计方法容易受到不良数据的影响，导致鲁棒性较差。而基于测量不确定性理论，在确定误差分布的前提下，不确定性测度可以通过定量评价的方法给出置信区间，并以区间的形式研究状态估计算法的鲁棒性。当量测点的状态估计值在范围内时，则表示该量测点为有效量测点。
配电网的不确定性一直是状态估计方法中的难点。现有技术通过区间分析方法对状态估计值进行求解时，是以潮流计算值为初值或参考潮流计算结果给出初值，在实际应用中存在计算效率低、收敛速度慢的问题。
发明内容
本发明提供了一种电力系统区间抗差状态估计方法、系统及装置，解决了现有区间分析方法计算效率低、收敛速度慢的技术问题。
本发明第一方面提供一种电力系统区间抗差状态估计方法，包括：
建立区间约束模型；
对所述区间约束模型进行求解，其中根据得到的状态变量区间值近似解确定初值，当满足收敛准则时结束求解过程，输出状态变量可行域；
根据状态变量可行域及潮流约束建立约束关系式，以量测点在一定置信水平上的相对偏差最小化作为目标函数，根据所述约束关系式和所述目标函数建立状态估计模型；
求解所述状态估计模型，得到状态估计结果。
根据本发明第一方面的一种能够实现的方式，所述建立区间约束模型，包括：
根据量测点在一定置信水平上的相对偏差建立量测点评价规则；
根据所述量测点评价规则建立用于判断量测点是否有效的评价函数。
根据本发明第一方面的一种能够实现的方式，所述根据量测点在一定置信水平上的相对偏差建立量测点评价规则，包括：
设置所述相对偏差的计算公式为：
Figure BDA0003344670170000021
式中，di表示量测点i在一定置信水平上的相对偏差，hi(x)为量测点i与状态变量x之间的关系函数，zi为量测点i的量测值，Ui为所述一定置信水平对应的量测不确定度；
设置量测点评价规则为：若|di|≤1，则量测点i为有效量测点；若|di|>1，则量测点i为无效量测点。
根据本发明第一方面的一种能够实现的方式，所述评价函数的表达式为：
Figure BDA0003344670170000022
式中，参数λ、k的取值根据工程经验设置，λ的取值范围为1～5，k的取值范围为2～4。
根据本发明第一方面的一种能够实现的方式，所述对所述区间约束模型进行求解，其中根据得到的状态变量区间值近似解确定初值，当满足收敛准则时结束求解过程，包括：
将所述评价函数的表达式转化为基于雅可比矩阵的线性规划方程；
采用基于迭代运算的线性规划方法求解所述线性规划方程，得到状态变量区间值近似解；
以状态变量区间值近似解的上限值和下限值的平均值作为初值，不断迭代求解修正所述线性规划方程，直到达到收敛准则时结束求解过程。
根据本发明第一方面的一种能够实现的方式，所述状态估计模型为：
Figure BDA0003344670170000031
s.t.g(x)＝0
Figure BDA0003344670170000032
式中，di表示量测点在一定置信水平上的相对偏差，N为量测点数量，x表示状态变量，
Figure BDA0003344670170000033
为状态变量可行域，g(x)＝0为潮流约束。
本发明第二方面提供一种电力系统区间抗差状态估计系统，包括：
区间约束模型建立模块，用于建立区间约束模型；
区间约束模型求解模块，用于对所述区间约束模型进行求解，其中根据得到的状态变量区间值近似解确定初值，当满足收敛准则时结束求解过程，输出状态变量可行域；
状态估计模型建立模块，用于根据状态变量可行域及潮流约束建立约束关系式，以量测点在一定置信水平上的相对偏差最小化作为目标函数，根据所述约束关系式和所述目标函数建立状态估计模型；
状态估计优化求解模块，用于求解所述状态估计模型，得到状态估计结果。
根据本发明第二方面的一种能够实现的方式，所述区间约束模型建立模块包括：
量测点评价规则建立单元，用于根据量测点在一定置信水平上的相对偏差建立量测点评价规则；
评价函数建立单元，用于根据所述量测点评价规则建立用于判断量测点是否有效的评价函数。
根据本发明第二方面的一种能够实现的方式，所述量测点评价规则建立单元包括：
相对偏差定义子单元，用于设置所述相对偏差的计算公式为：
Figure BDA0003344670170000034
式中，di表示量测点i在一定置信水平上的相对偏差，hi(x)为量测点i与状态变量x之间的关系函数，zi为量测点i的量测值，Ui为所述一定置信水平对应的量测不确定度；
量测点评价规则设置子单元，用于设置量测点评价规则为：若|di|≤1，则量测点i为有效量测点；若|di|>1，则量测点i为无效量测点。
根据本发明第二方面的一种能够实现的方式，所述评价函数的表达式为：
Figure BDA0003344670170000041
式中，参数λ、k的取值根据工程经验设置，λ的取值范围为1～5，k的取值范围为2～4。
根据本发明第二方面的一种能够实现的方式，所述区间约束模型求解模块包括：
转化单元，用于将所述评价函数的表达式转化为基于雅可比矩阵的线性规划方程；
求解单元，用于采用基于迭代运算的线性规划系统求解所述线性规划方程，得到状态变量区间值近似解；
迭代单元，用于以状态变量区间值近似解的上限值和下限值的平均值作为初值，不断迭代求解修正所述线性规划方程，直到达到收敛准则时结束求解过程。
根据本发明第二方面的一种能够实现的方式，所述状态估计模型为：
Figure BDA0003344670170000042
s.t.g(x)＝0
Figure BDA0003344670170000043
式中，di表示量测点在一定置信水平上的相对偏差，N为量测点数量，x表示状态变量，
Figure BDA0003344670170000044
为状态变量可行域，g(x)＝0为潮流约束。
本发明第三方面提供了一种电力系统区间抗差状态估计装置，包括：
存储器，用于存储指令；其中，所述指令为可实现如上任意一项能够实现的方式所述的电力系统区间抗差状态估计方法的指令；
处理器，用于执行所述存储器中的指令。
本发明第四方面一种计算机可读存储介质，所述计算机可读存储介质上存储有计算机程序，所述计算机程序被处理器执行时实现如上任意一项能够实现的方式所述的电力系统区间抗差状态估计方法。
从以上技术方案可以看出，本发明具有以下优点：
本发明上述实施例，在对建立的区间约束模型进行求解时，根据得到的状态变量区间值近似解确定初值，并以求解过程结束后输出的状态变量可行域为约束，以量测点在一定置信水平上的相对偏差最小化作为目标函数，建立对应的状态估计模型，进而对状态估计模型进行优化求解，得到状态估计结果；本发明实施例通过取状态变量的可行域为约束，以状态变量可行域设置初值，无需将潮流计算结果作为初值，可有效缩小状态变量的求解范围和区间范围，减少计算量，从而提高计算效率和收敛速度；
本发明实施例相比于传统的加权最小二乘状态估计，能够解决不良数据影响下区间估计精度下降的问题，有效抵抗量测坏数据对计算结果的影响，具有更高的鲁棒性和工程实用性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其它的附图。
图1为本发明一个可选实施例提供的一种电力系统区间抗差状态估计方法的流程图；
图2为本发明一个可选实施例提供的一种电力系统区间抗差状态估计系统的结构连接示意图。
附图说明：
1-区间约束模型建立模块；2-区间约束模型求解模块；3-状态估计模型建立模块；4-状态估计优化求解模块。
具体实施方式
本发明实施例提供了一种电力系统区间抗差状态估计方法、系统及装置，用于解决现有区间分析方法计算效率低、收敛速度慢的技术问题。
为使得本发明的发明目的、特征、优点能够更加的明显和易懂，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，下面所描述的实施例仅仅是本发明一部分实施例，而非全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例，都属于本发明保护的范围。
请参阅图1，图1示出了本发明实施例提供的一种电力系统区间抗差状态估计方法的流程图。
本发明提供的一种电力系统区间抗差状态估计方法，包括步骤S1-S4。
步骤S1，建立区间约束模型。
在一种能够实现的方式中，所述建立区间约束模型，包括：
根据量测点在一定置信水平上的相对偏差建立量测点评价规则；
根据所述量测点评价规则建立用于判断量测点是否有效的评价函数。
在一种能够实现的方式中，所述根据量测点在一定置信水平上的相对偏差建立量测点评价规则，包括：
设置所述相对偏差的计算公式为：
Figure BDA0003344670170000061
式中，di表示量测点i在一定置信水平上的相对偏差，hi(x)为量测点i与状态变量x之间的关系函数，zi为量测点i的量测值，Ui为所述一定置信水平对应的量测不确定度，其与测量装置的精度有关；
设置量测点评价规则为：若|di|≤1，则量测点i为有效量测点；若|di|>1，则量测点i为无效量测点。
本发明实施例的设置量测点评价规则简单便捷，易于实施。
在一种能够实现的方式中，所述评价函数的表达式为：
Figure BDA0003344670170000071
式中，参数λ、k的取值根据工程经验设置，λ的取值范围为1～5，k的取值范围为2～4。
步骤S2，对所述区间约束模型进行求解，其中根据得到的状态变量区间值近似解确定初值，当满足收敛准则时结束求解过程，输出状态变量可行域。
在一种能够实现的方式中，所述对所述区间约束模型进行求解，其中根据得到的状态变量区间值近似解确定初值，当满足收敛准则时结束求解过程，包括：
将所述评价函数的表达式转化为基于雅可比矩阵的线性规划方程；
采用基于迭代运算的线性规划方法求解所述线性规划方程，得到状态变量区间值近似解；
以状态变量区间值近似解的上限值和下限值的平均值作为初值，不断迭代求解修正所述线性规划方程，直到达到收敛准则时结束求解过程。
需要说明的是，该平均值可以是算术平均值、平方平均值或加权平均值。进一步地，还可以以状态变量区间值的中位数作为初值。
其中，该基于雅可比矩阵的线性规划方程如下：
Figure BDA0003344670170000072
Figure BDA0003344670170000073
式中，H为量测点与状态变量之间的关系函数hi(x)对应的雅可比矩阵，ai为系数，对应矩阵H-1中的第i行，Δzn对应矩阵Δz的第n行，Δxi 、
Figure BDA0003344670170000074
为修正值。
采用基于迭代运算的线性规划方法求解式(1)及式(2)，可以得到该修正值Δx、
Figure BDA0003344670170000081
因此，状态变量区间的上限值、下限值的近似解可以表示为：
Figure BDA0003344670170000082
其中，可以根据具体情况设置相应的收敛准则，本发明实施例不限定于此。
步骤S3，根据状态变量可行域及潮流约束建立约束关系式，以量测点在一定置信水平上的相对偏差最小化作为目标函数，根据所述约束关系式和所述目标函数建立状态估计模型。
在一种能够实现的方式中，所述状态估计模型为：
Figure BDA0003344670170000083
s.t.g(x)＝0
Figure BDA0003344670170000084
式中，di表示量测点在一定置信水平上的相对偏差，N为量测点数量，x表示状态变量，
Figure BDA0003344670170000085
为状态变量可行域，g(x)＝0为潮流约束。
步骤S4，求解所述状态估计模型，得到状态估计结果。
其中，步骤S4中状态估计优化问题的求解，可引入库恩塔克条件，再使用牛顿-拉夫逊法求解优化问题，也可以用其他优化问题求解方法来求解。
在一种实施方式中，可以引入松弛变量sxl和sxu，将状态估计模型中的不等式约束转化为下式所示的等式约束：
Figure BDA0003344670170000086
使用内点法求解上述优化问题，得到状态变量的最终值。
本发明还提供了一种电力系统区间抗差状态估计系统。
请参阅图2，图2示出了本发明实施例提供的一种电力系统区间抗差状态估计系统的结构连接示意图。
所述系统包括：
区间约束模型建立模块1，用于建立区间约束模型；
区间约束模型求解模块2，用于对所述区间约束模型进行求解，其中根据得到的状态变量区间值近似解确定初值，当满足收敛准则时结束求解过程，输出状态变量可行域；
状态估计模型建立模块3，用于根据状态变量可行域及潮流约束建立约束关系式，以量测点在一定置信水平上的相对偏差最小化作为目标函数，根据所述约束关系式和所述目标函数建立状态估计模型；
状态估计优化求解模块4，用于求解所述状态估计模型，得到状态估计结果。
在一种能够实现的方式中，所述区间约束模型建立模块1包括：
量测点评价规则建立单元，用于根据量测点在一定置信水平上的相对偏差建立量测点评价规则；
评价函数建立单元，用于根据所述量测点评价规则建立用于判断量测点是否有效的评价函数。
在一种能够实现的方式中，所述量测点评价规则建立单元包括：
相对偏差定义子单元，用于设置所述相对偏差的计算公式为：
Figure BDA0003344670170000091
式中，di表示量测点i在一定置信水平上的相对偏差，hi(x)为量测点i与状态变量x之间的关系函数，zi为量测点i的量测值，Ui为所述一定置信水平对应的量测不确定度；
量测点评价规则设置子单元，用于设置量测点评价规则为：若|di|≤1，则量测点λ为有效量测点；若|di|>1，则量测点i为无效量测点。
在一种能够实现的方式中，所述评价函数的表达式为：
Figure BDA0003344670170000092
式中，参数λ、k的取值根据工程经验设置，λ的取值范围为1～5，k的取值范围为2～4。
在一种能够实现的方式中，所述区间约束模型求解模块2包括：
转化单元，用于将所述评价函数的表达式转化为基于雅可比矩阵的线性规划方程；
求解单元，用于采用基于迭代运算的线性规划系统求解所述线性规划方程，得到状态变量区间值近似解；
迭代单元，用于以状态变量区间值近似解的上限值和下限值的平均值作为初值，不断迭代求解修正所述线性规划方程，直到达到收敛准则时结束求解过程。
在一种能够实现的方式中，所述状态估计模型为：
Figure BDA0003344670170000101
s.t.g(x)＝0
Figure BDA0003344670170000102
式中，di表示量测点在一定置信水平上的相对偏差，N为量测点数量，x表示状态变量，
Figure BDA0003344670170000103
为状态变量可行域，g(x)＝0为潮流约束。
本发明还提供了一种电力系统区间抗差状态估计装置，包括：
存储器，用于存储指令；其中，所述指令为可实现如上任意一项实施例所述的电力系统区间抗差状态估计方法的指令；
处理器，用于执行所述存储器中的指令。
本发明还提供了一种计算机可读存储介质，所述计算机可读存储介质上存储有计算机程序，所述计算机程序被处理器执行时实现如上任意一项实施例所述的电力系统区间抗差状态估计方法。
本发明上述实施例，在对建立的区间约束模型进行求解时，根据得到的状态变量区间值近似解确定初值，并以求解过程结束后输出的状态变量可行域为约束，以量测点在一定置信水平上的相对偏差最小化作为目标函数，建立对应的状态估计模型，进而对状态估计模型进行优化求解，得到状态估计结果；本发明实施例通过取状态变量的可行域为约束，以状态变量可行域设置初值，无需将潮流计算结果作为初值，可有效缩小状态变量的求解范围和区间范围，减少计算量，从而提高计算效率和收敛速度；相比于传统的加权最小二乘状态估计，能够解决不良数据影响下区间估计精度下降的问题，有效抵抗量测坏数据对计算结果的影响，具有更高的鲁棒性和工程实用性。
所属领域的技术人员可以清楚地了解到，为描述的方便和简洁，上述描述的系统，装置和单元的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。
在本申请所提供的几个实施例中，应该理解到，所揭露的系统，装置和方法，可以通过其它的方式实现。例如，以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。
所述集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-OnlyMemory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述，以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (14)

1.一种电力系统区间抗差状态估计方法，其特征在于，包括：
建立区间约束模型；
对所述区间约束模型进行求解，其中根据得到的状态变量区间值近似解确定初值，当满足收敛准则时结束求解过程，输出状态变量可行域；
根据状态变量可行域及潮流约束建立约束关系式，以量测点在一定置信水平上的相对偏差最小化作为目标函数，根据所述约束关系式和所述目标函数建立状态估计模型；
求解所述状态估计模型，得到状态估计结果。
2.根据权利要求1所述的电力系统区间抗差状态估计方法，其特征在于，所述建立区间约束模型，包括：
根据量测点在一定置信水平上的相对偏差建立量测点评价规则；
根据所述量测点评价规则建立用于判断量测点是否有效的评价函数。
3.根据权利要求2所述的电力系统区间抗差状态估计方法，其特征在于，所述根据量测点在一定置信水平上的相对偏差建立量测点评价规则，包括：
设置所述相对偏差的计算公式为：
Figure FDA0003344670160000011
式中，di表示量测点i在一定置信水平上的相对偏差，hi(x)为量测点i与状态变量x之间的关系函数，zi为量测点i的量测值，Ui为所述一定置信水平对应的量测不确定度；
设置量测点评价规则为：若|di|≤1，则量测点i为有效量测点；若|di|>1，则量测点i为无效量测点。
4.根据权利要求3所述的电力系统区间抗差状态估计方法，其特征在于，所述评价函数的表达式为：
Figure FDA0003344670160000012
式中，参数λ、k的取值根据工程经验设置，λ的取值范围为1～5，k的取值范围为2～4。
5.根据权利要求1所述的电力系统区间抗差状态估计方法，其特征在于，所述对所述区间约束模型进行求解，其中根据得到的状态变量区间值近似解确定初值，当满足收敛准则时结束求解过程，包括：
将所述评价函数的表达式转化为基于雅可比矩阵的线性规划方程；
采用基于迭代运算的线性规划方法求解所述线性规划方程，得到状态变量区间值近似解；
以状态变量区间值近似解的上限值和下限值的平均值作为初值，不断迭代求解修正所述线性规划方程，直到达到收敛准则时结束求解过程。
6.根据权利要求1所述的电力系统区间抗差状态估计方法，其特征在于，所述状态估计模型为：
Figure FDA0003344670160000021
s.t.g(x)＝0
Figure FDA0003344670160000022
式中，di表示量测点在一定置信水平上的相对偏差，N为量测点数量，x表示状态变量，
Figure FDA0003344670160000023
为状态变量可行域，g(x)＝0为潮流约束。
7.一种电力系统区间抗差状态估计系统，其特征在于，包括：
区间约束模型建立模块，用于建立区间约束模型；
区间约束模型求解模块，用于对所述区间约束模型进行求解，其中根据得到的状态变量区间值近似解确定初值，当满足收敛准则时结束求解过程，输出状态变量可行域；
状态估计模型建立模块，用于根据状态变量可行域及潮流约束建立约束关系式，以量测点在一定置信水平上的相对偏差最小化作为目标函数，根据所述约束关系式和所述目标函数建立状态估计模型；
状态估计优化求解模块，用于求解所述状态估计模型，得到状态估计结果。
8.根据权利要求7所述的电力系统区间抗差状态估计系统，其特征在于，所述区间约束模型建立模块包括：
量测点评价规则建立单元，用于根据量测点在一定置信水平上的相对偏差建立量测点评价规则；
评价函数建立单元，用于根据所述量测点评价规则建立用于判断量测点是否有效的评价函数。
9.根据权利要求8所述的电力系统区间抗差状态估计系统，其特征在于，所述量测点评价规则建立单元包括：
相对偏差定义子单元，用于设置所述相对偏差的计算公式为：
Figure FDA0003344670160000031
式中，di表示量测点i在一定置信水平上的相对偏差，hi(x)为量测点i与状态变量x之间的关系函数，zi为量测点i的量测值，Ui为所述一定置信水平对应的量测不确定度；
量测点评价规则设置子单元，用于设置量测点评价规则为：若|di|≤1，则量测点i为有效量测点；若|di|>1，则量测点i为无效量测点。
10.根据权利要求9所述的电力系统区间抗差状态估计系统，其特征在于，所述评价函数的表达式为：
Figure FDA0003344670160000032
式中，参数λ、k的取值根据工程经验设置，λ的取值范围为1～5，k的取值范围为2～4。
11.根据权利要求7所述的电力系统区间抗差状态估计系统，其特征在于，所述区间约束模型求解模块包括：
转化单元，用于将所述评价函数的表达式转化为基于雅可比矩阵的线性规划方程；
求解单元，用于采用基于迭代运算的线性规划系统求解所述线性规划方程，得到状态变量区间值近似解；
迭代单元，用于以状态变量区间值近似解的上限值和下限值的平均值作为初值，不断迭代求解修正所述线性规划方程，直到达到收敛准则时结束求解过程。
12.根据权利要求7所述的电力系统区间抗差状态估计系统，其特征在于，所述状态估计模型为：
Figure FDA0003344670160000041
s.t.g(x)＝0
Figure FDA0003344670160000042
式中，di表示量测点在一定置信水平上的相对偏差，N为量测点数量，x表示状态变量，
Figure FDA0003344670160000043
为状态变量可行域，g(x)＝0为潮流约束。
13.一种电力系统区间抗差状态估计装置，其特征在于，包括：
存储器，用于存储指令；其中，所述指令为可实现如权利要求1-6任意一项所述的电力系统区间抗差状态估计方法的指令；
处理器，用于执行所述存储器中的指令。
14.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质上存储有计算机程序，所述计算机程序被处理器执行时实现如权利要求1-6任意一项所述的电力系统区间抗差状态估计方法。
CN202111319373.1A 2021-11-09 2021-11-09 一种电力系统区间抗差状态估计方法、系统及装置 Pending CN113962108A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111319373.1A CN113962108A (zh) 2021-11-09 2021-11-09 一种电力系统区间抗差状态估计方法、系统及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111319373.1A CN113962108A (zh) 2021-11-09 2021-11-09 一种电力系统区间抗差状态估计方法、系统及装置

Publications (1)

Publication Number Publication Date
CN113962108A true CN113962108A (zh) 2022-01-21

Family

ID=79469651

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111319373.1A Pending CN113962108A (zh) 2021-11-09 2021-11-09 一种电力系统区间抗差状态估计方法、系统及装置

Country Status (1)

Country Link
CN (1) CN113962108A (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109711662A (zh) * 2018-11-14 2019-05-03 华北电力大学 一种基于多源数据融合的电网抗差状态估计方法
CN112993989A (zh) * 2021-03-05 2021-06-18 广东电网有限责任公司广州供电局 一种主动配电系统抗差状态估计数据处理方法
CN113139288A (zh) * 2021-04-21 2021-07-20 广东电网有限责任公司广州供电局 基于gpu二次加速的配电网分布式抗差状态估计方法及装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109711662A (zh) * 2018-11-14 2019-05-03 华北电力大学 一种基于多源数据融合的电网抗差状态估计方法
CN112993989A (zh) * 2021-03-05 2021-06-18 广东电网有限责任公司广州供电局 一种主动配电系统抗差状态估计数据处理方法
CN113139288A (zh) * 2021-04-21 2021-07-20 广东电网有限责任公司广州供电局 基于gpu二次加速的配电网分布式抗差状态估计方法及装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
曾顺奇等: "基于动态区间优化方法的主动配电网抗差状态估计", 《南方电网技术》 *
朱华统等: "《GPS坐标系统的变换》", 31 December 1994 *

Similar Documents

Publication Publication Date Title
Judd et al. Asymptotic methods for aggregate growth models
WO2021208079A1 (zh) 获取动力电池寿命数据的方法、装置、计算机设备及介质
CN108089135B (zh) 一种基于极限学习模型的电池状态预测系统及其实现方法
Meng et al. A seventh‐degree cubature Kalman filter
CN106406830B (zh) 一种准周期信号的处理方法、装置和计算机可读存储介质
CN113505490A (zh) 基于gmm的电力系统数字孪生体参数校正方法及装置
CN112818588A (zh) 一种电力系统的最优潮流计算方法、装置及存储介质
CN110783918A (zh) 一种基于线性模型的配电三相区间状态估计求解算法
CN114897277B (zh) 一种基于lstm的突发型滑坡位移预测方法
CN110781626A (zh) 有限差分多重分辨三角函数weno格式的模拟方法
CN114066699A (zh) 碳排放量测算方法、装置及终端设备
CN114239796A (zh) 一种基于扩展卡尔曼滤波的电力系统状态估计方法
CN113962108A (zh) 一种电力系统区间抗差状态估计方法、系统及装置
CN109993374B (zh) 货物量预测方法及装置
CN114372036A (zh) 电力系统的状态估计方法、装置、设备和计算机存储介质
CN104156953A (zh) 一种基于离散正交矩的图像重构方法
CN114117349A (zh) 电力系统随机变量概率分析方法、系统、设备及存储介质
CN113432608A (zh) 适于ins/cns组合导航系统的基于最大相关熵的广义高阶ckf算法
CN107730387B (zh) 一种支付数据的校正方法
CN116029164B (zh) 设备性能退化程度确定方法、系统、电子设备及存储介质
CN113640115B (zh) 适用于准等熵压缩实验数据逆问题求解的优化方法和系统
CN114492816B (zh) 量子态处理方法、计算设备、计算装置及存储介质
CN116029480B (zh) 代理购电测算方法及其系统
CN114580578B (zh) 具有约束的分布式随机优化模型训练方法、装置及终端
CN115511047B (zh) Softmax模型的量化方法、装置、设备及介质

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20220121

RJ01 Rejection of invention patent application after publication