CN111181153B - 一种基于支路响应信息确定戴维南等值参数的方法及*** - Google Patents

一种基于支路响应信息确定戴维南等值参数的方法及*** Download PDF

Info

Publication number
CN111181153B
CN111181153B CN201911141762.2A CN201911141762A CN111181153B CN 111181153 B CN111181153 B CN 111181153B CN 201911141762 A CN201911141762 A CN 201911141762A CN 111181153 B CN111181153 B CN 111181153B
Authority
CN
China
Prior art keywords
power system
thevenin
internal potential
reactance
iteration
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
Application number
CN201911141762.2A
Other languages
English (en)
Other versions
CN111181153A (zh
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.)
China Electric Power Research Institute Co Ltd CEPRI
State Grid Hunan Electric Power Co Ltd
Economic and Technological Research Institute of State Grid Hunan Electric Power Co Ltd
Original Assignee
China Electric Power Research Institute Co Ltd CEPRI
State Grid Hunan Electric Power Co Ltd
Economic and Technological Research Institute of State Grid Hunan Electric Power 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 China Electric Power Research Institute Co Ltd CEPRI, State Grid Hunan Electric Power Co Ltd, Economic and Technological Research Institute of State Grid Hunan Electric Power Co Ltd filed Critical China Electric Power Research Institute Co Ltd CEPRI
Priority to CN201911141762.2A priority Critical patent/CN111181153B/zh
Publication of CN111181153A publication Critical patent/CN111181153A/zh
Application granted granted Critical
Publication of CN111181153B publication Critical patent/CN111181153B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/12—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明公开了一种基于支路响应信息确定戴维南等值参数的方法及***，属于电力***安全稳定的电压稳定在线分析和监测技术领域。本发明方法，包括：选择电力***内多条母线中的任意一条母线作为分界，将母线以内的部分作为待测目标；获取待测目标母线节点的电压幅值、有功功率和无功功率；获取迭代轮数t及t+1轮的电力***戴维南内电势幅值、戴维南内电势相角和戴维南电抗；确定迭代轮数t及t+1轮的电力***戴维南内电势幅值、戴维南内电势相角和戴维南电抗的差值的绝对值是否小于预设值；确定迭代轮数t+1轮的电力***戴维南内电势幅值、戴维南内电势相角和戴维南电抗为电力***戴维南等值参数。本发明保障电力***的电压稳定。

Description

一种基于支路响应信息确定戴维南等值参数的方法及***
技术领域
本发明涉及电力***安全稳定的电压稳定在线分析和监测技术领域，并且更具体地，涉及一种基于支路响应信息确定戴维南等值参数的方法及***。
背景技术
我国东部地区的大部分省级电网均为受端电网。为了确保受端电网安全可靠运行不发生电压崩溃，需要进行在线电压稳定监测与分析。电力***戴维南等值是实现在线电压稳定监测与分析的基础技术。
戴维南等值将电力***中任意选定的某子网络等值为内电势与电抗串联的形式。戴维南等值既可以有效地简化复杂电力***，也可以利用其得到的等值内电势与内阻抗来精准地指示电力***的电压稳定性与稳定裕度，从而提升调度人员对电网安全态势的有效感知，防止***出现电压失稳。戴维南等值参数与电网结构和***运行方式紧密相关并随之变化，而戴维南等值参数的准确性直接决定调度人员了对电网安全特性认识的准确度，因此实现各种工况条件下戴维南等值参数的准确确定关乎对电网电压稳定的正确判断，是一项极为重要的工作。
传统的戴维南等值参数确定方法适用于如下场景：在若干个连续采样的相邻时步间，选定的某子网络的戴维南内电势幅值、内电势相角及戴维南电抗均保持不变；与此同时，该子网络以外的外部网络的有功、无功、电压幅值在这些相邻时步间需要有明显的变化。
事实上，在电力***真实的运行工况常常与上述场景差别较大。一个典型的场景是：在对电力***某子网络进行戴维南等值的多个时步内，为调节***电压水平使其保持在合理范围内，除开该子网络以外的外部网络中经常存在无功设备(电感或电容)的投切操作。在这种场景下，该子网络的戴维南等值内电势会显著变化，传统方法存在较大的误差，不再适用。
传统方法仅在戴维南内电势幅值、内电势相角及戴维南电抗在若干个相邻时步上均保持不变的场景下才能使用。在戴维南内电势幅值发生变化的工况下，传统方法得到的戴维南等值参数精度差，难以支撑准确的电压稳定判别。
发明内容
针对上述问题，本发明提出了一种基于支路响应信息确定戴维南等值参数的方法，包括：
选择电力***内多条母线中的任意一条母线作为分界，将母线以内的部分作为待测目标；
根据预设时步，获取待测目标母线节点的电压幅值、待测目标母线线路的有功功率和待测目标母线线路的无功功率；
根据所述待电压幅值、有功功率和无功功率，对电力***戴维南内电势幅值、内电势相角和电抗进行迭代，获取迭代轮数t及t+1轮的电力***戴维南内电势幅值、内电势相角和电抗，所述t大于或等于1；
根据迭代获取迭代轮数t及t+1轮的电力***戴维南内电势幅值、内电势相角和电抗，并确定迭代轮数t及t+1轮的电力***戴维南内电势幅值的差值的绝对值、电力***戴维南内电势相角差值的绝对值和电力***戴维南电抗的差值的绝对值是否小于预设值；
当迭代轮数t及t+1轮数的电力***戴维南内电势幅值的差值的绝对值、内电势相角差值的绝对值和电抗的差值的绝对值均小于预设值时，确定迭代轮数t+1轮的电力***戴维南内电势幅值、内电势相角和电抗为电力***戴维南等值参数。
可选的，预设时步大于或等于2。
可选的，预设值为0.001。
可选的，获取迭代轮数t及t+1轮的戴维南内电势幅值、戴维南内电势相角和电力***戴维南电抗，具体为：
令迭代轮数1轮电力***戴维南内电势幅值为1、迭代轮数1轮电力***戴维南内电势相角为0和迭代轮数1轮电力***戴维南电抗为0进行迭代，获取迭代轮数t轮的电力***戴维南内电势幅值、电力***戴维南内电势相角和电力***戴维南电抗；
根据迭代轮数t轮的电力***戴维南内电势幅值、电力***戴维南内电势相角和电力***戴维南电抗，获取迭代轮数t+1轮电力***戴维南内电势幅值、电力***戴维南内电势相角和电力***戴维南电抗，公式如下：
Figure BDA0002281128000000031
Figure BDA0002281128000000032
Figure BDA0002281128000000033
Figure BDA0002281128000000034
Figure BDA0002281128000000035
其中，η取0.001～0.01、
Figure BDA0002281128000000036
为迭代轮数t+1轮电力***戴维南内电势幅值、δth(t+1)迭代轮数t+1轮电力***戴维南内电势相角、Xth(t+1)迭代轮数t+1轮电力***戴维南电抗、
Figure BDA0002281128000000037
为迭代轮数t轮电力***戴维南内电势幅值、δth(t)迭代轮数t轮电力***戴维南内电势相角、Xth(t)迭代轮数t轮电力***戴维南电抗、Uk为待测目标母线节点的电压幅值、Pk为待测目标母线线路的有功功率、Qk为待测目标母线线路的无功功率、K为预设时步、k为任意一个时步、fk为有功残差和ek为无功残差。
本发明还提供了一种基于支路响应信息确定戴维南等值参数的***，包括：
目标确定模块，选择电力***内多条母线中的任意一条母线作为分界，将母线以内的部分作为待测目标；
参数获取模块，根据预设时步，获取待测目标母线节点的电压幅值、待测目标母线线路的有功功率和待测目标母线线路的无功功率；
迭代模块，根据所述待电压幅值、有功功率和无功功率，对电力***戴维南内电势幅值、内电势相角和电抗进行迭代，获取迭代轮数t及t+1轮的电力***戴维南内电势幅值、内电势相角和电抗，所述t大于或等于1；
比对模块，根据迭代获取迭代轮数t及t+1轮的电力***戴维南内电势幅值、内电势相角和电抗，并确定迭代轮数t及t+1轮的电力***戴维南内电势幅值的差值的绝对值、电力***戴维南内电势相角差值的绝对值和电力***戴维南电抗的差值的绝对值是否小于预设值；
判断模块，确定当迭代轮数t及t+1轮数的电力***戴维南内电势幅值的差值的绝对值、内电势相角差值的绝对值和电抗的差值的绝对值均小于预设值时，确定迭代轮数t+1轮的电力***戴维南内电势幅值、内电势相角和电抗为电力***戴维南等值参数。
可选的，预设时步大于或等于2。
可选的，预设值为0.001。
可选的，迭代模块用于：
令迭代轮数1轮电力***戴维南内电势幅值为1、迭代轮数1轮电力***戴维南内电势相角为0和迭代轮数1轮电力***戴维南电抗为0进行迭代，获取迭代轮数t轮的电力***戴维南内电势幅值、电力***戴维南内电势相角和电力***戴维南电抗；
根据迭代轮数t轮的电力***戴维南内电势幅值、电力***戴维南内电势相角和电力***戴维南电抗，获取迭代轮数t+1轮电力***戴维南内电势幅值、电力***戴维南内电势相角和电力***戴维南电抗，公式如下：
Figure BDA0002281128000000051
Figure BDA0002281128000000052
Figure BDA0002281128000000053
Figure BDA0002281128000000054
Figure BDA0002281128000000055
其中，η取0.001～0.01、
Figure BDA0002281128000000056
为迭代轮数t+1轮电力***戴维南内电势幅值、δth(t+1)迭代轮数t+1轮电力***戴维南内电势相角、Xth(t+1)迭代轮数t+1轮电力***戴维南电抗、
Figure BDA0002281128000000057
为迭代轮数t轮电力***戴维南内电势幅值、δth(t)迭代轮数t轮电力***戴维南内电势相角、Xth(t)迭代轮数t轮电力***戴维南电抗、Uk为待测目标母线节点的电压幅值、Pk为待测目标母线线路的有功功率、Qk为待测目标母线线路的无功功率、K为预设时步、k为任意一个时步、fk为有功残差和ek为无功残差。
本发明可以求得准确的戴维南等值参数，从而支撑准确的电压稳定判别，提升调度人员对电网安全态势的有效感知，防止电力***出现电压失稳，保障电力***的电压稳定。
本发明仅利用母线及其出线上的局部量测数据进行戴维南等值参数的计算，不依赖***模型和参数，对于戴维南内电势幅值变化的工况，本方法得到的戴维南等值参数相比于传统方法更加精确。
附图说明
图1为本发明一种基于支路响应信息确定戴维南等值参数的方法流程图；
图2为本发明一种基于支路响应信息确定戴维南等值参数的方法电力***戴维南等值示意图；
图3为本发明一种基于支路响应信息确定戴维南等值参数的方法待测目标母线线路的有功功率和待测目标母线线路的无功功率曲线图；
图4为本发明一种基于支路响应信息确定戴维南等值参数的方法待测目标母线节点的电压幅值曲线图；
图5为本发明一种基于支路响应信息确定戴维南等值参数的方法电力***戴维南内电势幅值曲线图；
图6为本发明一种基于支路响应信息确定戴维南等值参数的方法电力***戴维南内电势相角曲线图；
图7为本发明一种基于支路响应信息确定戴维南等值参数的方法力***戴维南电抗曲线图；
图8为本发明一种基于支路响应信息确定戴维南等值参数的***结构图。
具体实施方式
现在参考附图介绍本发明的示例性实施方式，然而，本发明可以用许多不同的形式来实施，并且不局限于此处描述的实施例，提供这些实施例是为了详尽地且完全地公开本发明，并且向所属技术领域的技术人员充分传达本发明的范围。对于表示在附图中的示例性实施方式中的术语并不是对本发明的限定。在附图中，相同的单元/元件使用相同的附图标记。
除非另有说明，此处使用的术语(包括科技术语)对所属技术领域的技术人员具有通常的理解含义。另外，可以理解的是，以通常使用的词典限定的术语，应当被理解为与其相关领域的语境具有一致的含义，而不应该被理解为理想化的或过于正式的意义。
本发明提出了一种基于支路响应信息确定戴维南等值参数的方法，如图1所示，包括：
电力***戴维南等值示意图，如图2所示，选择电力***内多条母线中的任意一条母线作为分界，将母线以内的部分作为待测目标，待测目标为子***，子***分界外部作为外部***，Uk为待测目标母线节点的电压幅值、Pk为待测目标母线线路的有功功率、Qk为待测目标母线线路的无功功率、
Figure BDA0002281128000000071
为电力***戴维南内电势幅值、δth电力***戴维南内电势相角和Xth为轮电力***戴维南电抗；
根据预设时步，预设时步大于或等于2；
获取待测目标母线节点的电压幅值、待测目标母线线路的有功功率，待测目标母线节点的电压幅值、待测目标母线线路的有功功率曲线图如图3所示；
获取待测目标母线线路的无功功率，待测目标母线线路的无功功率曲线图如图4所示；
根据所述待测目标母线节点的电压幅值、待测目标母线线路的有功功率和待测目标母线线路的无功功率，对电力***戴维南内电势幅值、电力***戴维南内电势相角和电力***戴维南电抗进行迭代，获取迭代轮数t及t+1轮的电力***戴维南内电势幅值、电力***戴维南内电势相角和电力***戴维南电抗，所述t大于或等于1；
获取迭代轮数t及t+1轮的戴维南内电势幅值、戴维南内电势相角和电力***戴维南电抗，具体为：
令迭代轮数1轮电力***戴维南内电势幅值为1、迭代轮数1轮电力***戴维南内电势相角为0和迭代轮数1轮电力***戴维南电抗为0进行迭代，获取迭代轮数t轮的电力***戴维南内电势幅值、电力***戴维南内电势相角和电力***戴维南电抗；
根据迭代轮数t轮的电力***戴维南内电势幅值、电力***戴维南内电势相角和电力***戴维南电抗，获取迭代轮数t+1轮电力***戴维南内电势幅值、电力***戴维南内电势相角和电力***戴维南电抗，公式如下：
Figure BDA0002281128000000081
Figure BDA0002281128000000082
Figure BDA0002281128000000083
Figure BDA0002281128000000084
Figure BDA0002281128000000085
其中，η取0.001～0.01、
Figure BDA0002281128000000086
为迭代轮数t+1轮电力***戴维南内电势幅值、δth(t+1)迭代轮数t+1轮电力***戴维南内电势相角、Xth(t+1)迭代轮数t+1轮电力***戴维南电抗、
Figure BDA0002281128000000087
为迭代轮数t轮电力***戴维南内电势幅值、δth(t)迭代轮数t轮电力***戴维南内电势相角、Xth(t)迭代轮数t轮电力***戴维南电抗、Uk为待测目标母线节点的电压幅值、Pk为待测目标母线线路的有功功率、Qk为待测目标母线线路的无功功率、K为预设时步、k为任意一个时步、fk为有功残差和ek为无功残差。
根据迭代获取迭代轮数t及t+1轮的电力***戴维南内电势幅值、电力***戴维南内电势相角和电力***戴维南电抗，并确定迭代轮数t及t+1轮的电力***戴维南内电势幅值的差值的绝对值、电力***戴维南内电势相角差值的绝对值和电力***戴维南电抗的差值的绝对值是否小于预设值；
当迭代轮数t及t+1轮数的电力***戴维南内电势幅值的差值的绝对值、电力***戴维南内电势相角差值的绝对值和电力***戴维南电抗的差值的绝对值均小于预设值时，预设值为0.001，确定迭代轮数t+1轮的电力***戴维南内电势幅值、电力***戴维南内电势相角和电力***戴维南电抗为电力***戴维南等值参数；
电力***戴维南内电势幅值、电力***戴维南内电势相角和电力***戴维南电抗的曲线图，分别如图5、6和7所示。
本发明还提供了一种基于支路响应信息确定戴维南等值参数的***200，如图8所示，包括：
目标确定模块201，选择电力***内多条母线中的任意一条母线作为分界，将母线以内的部分作为待测目标；
参数获取模块202，根据预设时步，预设时步大于或等于2，获取待测目标母线节点的电压幅值、待测目标母线线路的有功功率和待测目标母线线路的无功功率；
迭代模块203，根据所述待测目标母线节点的电压幅值、待测目标母线线路的有功功率和待测目标母线线路的无功功率，对电力***戴维南内电势幅值、电力***戴维南内电势相角和电力***戴维南电抗进行迭代，获取迭代轮数t及t+1轮的电力***戴维南内电势幅值、电力***戴维南内电势相角和电力***戴维南电抗，所述t大于或等于1；
迭代模块203用于：
令迭代轮数1轮电力***戴维南内电势幅值为1、迭代轮数1轮电力***戴维南内电势相角为0和迭代轮数1轮电力***戴维南电抗为0进行迭代，获取迭代轮数t轮的电力***戴维南内电势幅值、电力***戴维南内电势相角和电力***戴维南电抗；
根据迭代轮数t轮的电力***戴维南内电势幅值、电力***戴维南内电势相角和电力***戴维南电抗，获取迭代轮数t+1轮电力***戴维南内电势幅值、电力***戴维南内电势相角和电力***戴维南电抗，公式如下：
Figure BDA0002281128000000101
Figure BDA0002281128000000102
Figure BDA0002281128000000103
Figure BDA0002281128000000104
Figure BDA0002281128000000105
其中，η取0.001～0.01、
Figure BDA0002281128000000106
为迭代轮数t+1轮电力***戴维南内电势幅值、δth(t+1)迭代轮数t+1轮电力***戴维南内电势相角、Xth(t+1)迭代轮数t+1轮电力***戴维南电抗、
Figure BDA0002281128000000107
为迭代轮数t轮电力***戴维南内电势幅值、δth(t)迭代轮数t轮电力***戴维南内电势相角、Xth(t)迭代轮数t轮电力***戴维南电抗、Uk为待测目标母线节点的电压幅值、Pk为待测目标母线线路的有功功率、Qk为待测目标母线线路的无功功率、K为预设时步、k为任意一个时步、fk为有功残差和ek为无功残差。
比对模块204，根据迭代获取迭代轮数t及t+1轮的电力***戴维南内电势幅值、电力***戴维南内电势相角和电力***戴维南电抗，并确定迭代轮数t及t+1轮的电力***戴维南内电势幅值的差值的绝对值、电力***戴维南内电势相角差值的绝对值和电力***戴维南电抗的差值的绝对值是否小于预设值；
判断模块205，确定当迭代轮数t及t+1轮数的电力***戴维南内电势幅值的差值的绝对值、电力***戴维南内电势相角差值的绝对值和电力***戴维南电抗的差值的绝对值均小于预设值时，预设值为0.001，确定迭代轮数t+1轮的电力***戴维南内电势幅值、电力***戴维南内电势相角和电力***戴维南电抗为电力***戴维南等值参数。
本发明可以求得准确的戴维南等值参数，从而支撑准确的电压稳定判别，提升调度人员对电网安全态势的有效感知，防止***出现电压失稳，保障电力***的电压稳定。
本发明仅利用母线及其出线上的局部量测数据进行戴维南等值参数的计算，不依赖***模型和参数，对于戴维南内电势幅值变化的工况，本方法得到的戴维南等值参数相比于传统方法更加精确。
本领域内的技术人员应明白，本申请的实施例可提供为方法、***、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(***)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
最后应当说明的是：以上实施例仅用以说明本发明的技术方案而非对其限制，尽管参照上述实施例对本发明进行了详细的说明，所属领域的普通技术人员应当理解：依然可以对本发明的具体实施方式进行修改或者等同替换，而未脱离本发明精神和范围的任何修改或者等同替换，其均应涵盖在本发明的权利要求保护范围之内。

Claims (6)

1.一种基于支路响应信息确定戴维南等值参数的方法，所述方法包括：
选择电力***内多条母线中的任意一条母线作为分界，将母线以内的部分作为待测目标；
根据预设时步，获取待测目标母线节点的电压幅值、待测目标母线线路的有功功率和待测目标母线线路的无功功率；
根据所述待测目标母线节点的电压幅值、有功功率和无功功率，对电力***戴维南内电势幅值、内电势相角和电抗进行迭代，获取迭代轮数t及t+1轮的电力***戴维南内电势幅值、内电势相角和电抗，所述t大于或等于1；
所述获取迭代轮数t及t+1轮的电力***戴维南内电势幅值、内电势相角和电抗，具体为：
令迭代轮数1轮电力***戴维南内电势幅值为1、迭代轮数1轮电力***戴维南内电势相角为0和迭代轮数1轮电力***戴维南电抗为0进行迭代，获取迭代轮数t轮的电力***戴维南内电势幅值、电力***戴维南内电势相角和电力***戴维南电抗；
根据迭代轮数t轮的电力***戴维南内电势幅值、电力***戴维南内电势相角和电力***戴维南电抗，获取迭代轮数t+1轮电力***戴维南内电势幅值、电力***戴维南内电势相角和电力***戴维南电抗，公式如下：
Figure FDA0003820231810000021
Figure FDA0003820231810000022
Figure FDA0003820231810000023
Figure FDA0003820231810000024
Figure FDA0003820231810000025
其中，η取0.001～0.01、
Figure FDA0003820231810000026
为迭代轮数t+1轮电力***戴维南内电势幅值、δth(t+1)迭代轮数t+1轮电力***戴维南内电势相角、Xth(t+1)迭代轮数t+1轮电力***戴维南电抗、
Figure FDA0003820231810000027
为迭代轮数t轮电力***戴维南内电势幅值、δth(t)迭代轮数t轮电力***戴维南内电势相角、Xth(t)迭代轮数t轮电力***戴维南电抗、Uk为待测目标母线节点的电压幅值、Pk为待测目标母线线路的有功功率、Qk为待测目标母线线路的无功功率、K为预设时步、k为任意一个时步、fk为有功残差和ek为无功残差；
根据迭代获取迭代轮数t及t+1轮的电力***戴维南内电势幅值、内电势相角和电抗，并确定迭代轮数t及t+1轮的电力***戴维南内电势幅值的差值的绝对值、电力***戴维南内电势相角差值的绝对值和电力***戴维南电抗的差值的绝对值是否小于预设值；
当迭代轮数t及t+1轮数的电力***戴维南内电势幅值的差值的绝对值、内电势相角差值的绝对值和电抗的差值的绝对值均小于预设值时，确定迭代轮数t+1轮的电力***戴维南内电势幅值、内电势相角和电抗为电力***戴维南等值参数。
2.根据权利要求1所述的方法，所述预设时步大于或等于2。
3.根据权利要求1所述的方法，所述预设值为0.001。
4.一种基于支路响应信息确定戴维南等值参数的***，所述***包括：
目标确定模块，选择电力***内多条母线中的任意一条母线作为分界，将母线以内的部分作为待测目标；
参数获取模块，根据预设时步，获取待测目标母线节点的电压幅值、待测目标母线线路的有功功率和待测目标母线线路的无功功率；
迭代模块，根据所述待测目标母线节点的电压幅值、有功功率和无功功率，对电力***戴维南内电势幅值、内电势相角和电抗进行迭代，获取迭代轮数t及t+1轮的电力***戴维南内电势幅值、内电势相角和电抗，所述t大于或等于1；
所述获取迭代轮数t及t+1轮的电力***戴维南内电势幅值、内电势相角和电抗，具体为：
令迭代轮数1轮电力***戴维南内电势幅值为1、迭代轮数1轮电力***戴维南内电势相角为0和迭代轮数1轮电力***戴维南电抗为0进行迭代，获取迭代轮数t轮的电力***戴维南内电势幅值、电力***戴维南内电势相角和电力***戴维南电抗；
根据迭代轮数t轮的电力***戴维南内电势幅值、电力***戴维南内电势相角和电力***戴维南电抗，获取迭代轮数t+1轮电力***戴维南内电势幅值、电力***戴维南内电势相角和电力***戴维南电抗，公式如下：
Figure FDA0003820231810000031
Figure FDA0003820231810000032
Figure FDA0003820231810000033
Figure FDA0003820231810000034
Figure FDA0003820231810000035
其中，η取0.001～0.01、
Figure FDA0003820231810000036
为迭代轮数t+1轮电力***戴维南内电势幅值、δth(t+1)迭代轮数t+1轮电力***戴维南内电势相角、Xth(t+1)迭代轮数t+1轮电力***戴维南电抗、
Figure FDA0003820231810000041
为迭代轮数t轮电力***戴维南内电势幅值、δth(t)迭代轮数t轮电力***戴维南内电势相角、Xth(t)迭代轮数t轮电力***戴维南电抗、Uk为待测目标母线节点的电压幅值、Pk为待测目标母线线路的有功功率、Qk为待测目标母线线路的无功功率、K为预设时步、k为任意一个时步、fk为有功残差和ek为无功残差；
比对模块，根据迭代获取迭代轮数t及t+1轮的电力***戴维南内电势幅值、内电势相角和电抗，并确定迭代轮数t及t+1轮的电力***戴维南内电势幅值的差值的绝对值、电力***戴维南内电势相角差值的绝对值和电力***戴维南电抗的差值的绝对值是否小于预设值；
判断模块，确定当迭代轮数t及t+1轮数的电力***戴维南内电势幅值的差值的绝对值、内电势相角差值的绝对值和电抗的差值的绝对值均小于预设值时，确定迭代轮数t+1轮的电力***戴维南内电势幅值、内电势相角和电抗为电力***戴维南等值参数。
5.根据权利要求4所述的***，所述预设时步大于或等于2。
6.根据权利要求4所述的***，所述预设值为0.001。
CN201911141762.2A 2019-11-20 2019-11-20 一种基于支路响应信息确定戴维南等值参数的方法及*** Active CN111181153B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911141762.2A CN111181153B (zh) 2019-11-20 2019-11-20 一种基于支路响应信息确定戴维南等值参数的方法及***

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911141762.2A CN111181153B (zh) 2019-11-20 2019-11-20 一种基于支路响应信息确定戴维南等值参数的方法及***

Publications (2)

Publication Number Publication Date
CN111181153A CN111181153A (zh) 2020-05-19
CN111181153B true CN111181153B (zh) 2022-10-11

Family

ID=70651889

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911141762.2A Active CN111181153B (zh) 2019-11-20 2019-11-20 一种基于支路响应信息确定戴维南等值参数的方法及***

Country Status (1)

Country Link
CN (1) CN111181153B (zh)

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102522742B (zh) * 2011-10-26 2014-02-26 重庆大学 基于单点量测信息的外网戴维南等值参数的估计方法
CN102420427B (zh) * 2011-12-02 2014-05-14 余杭供电局 一种考虑外网等值的区域电网电压稳定裕度计算方法
CN103872679A (zh) * 2014-03-07 2014-06-18 国家电网公司 微弱信号下的电网戴维南等值模型参数辨识方法
CN104778367B (zh) * 2015-04-20 2017-08-08 山东大学 基于单一状态断面的广域戴维南等值参数在线计算方法
GB2567460B (en) * 2017-10-12 2020-02-12 Univ Cape Town Method and apparatus for obtaining Thévenin equivalent parameters of a power network
CN107958120B (zh) * 2017-11-30 2021-10-08 上海电力学院 一种基于幂级数展开的***戴维南等值参数计算方法
CN108336730B (zh) * 2018-03-14 2021-07-13 哈尔滨工业大学 一种基于降阶自灵敏度的戴维南等值参数辨识方法

Also Published As

Publication number Publication date
CN111181153A (zh) 2020-05-19

Similar Documents

Publication Publication Date Title
CN110579682B (zh) 一种故障录波数据暂态同源对比方法及装置
CN106505557B (zh) 一种遥测错误辨识方法及装置
CN109633366B (zh) 一种配电网单相接地故障选相方法及选相装置
CN106709147B (zh) 一种基于bpa的短路故障极限切除时间批量自动求取方法
CN109284933B (zh) 一种基于数理统计的电子式互感器状态评估***及方法
CN111130074A (zh) 基于电压首极值时间的直流输电线超高速保护***及方法
CN112989587A (zh) 一种电容式电压互感器劣化成因在线分析方法及***
CN109829246B (zh) 一种基于参数综合可疑度的线路参数辨识方法
CN113484804B (zh) 一种应用于馈线终端装置的错接线识别方法
CN111181153B (zh) 一种基于支路响应信息确定戴维南等值参数的方法及***
CN109901003B (zh) 一种逆变器功率故障检测方法及***
CN110190617B (zh) 多馈入直流电力***的评估方法、***、设备及存储介质
CN114065118B (zh) 一种基于指数函数的电力***鲁棒状态估计方法
CN104849684A (zh) 一种示波器、校正装置及其自动校正水平中心的方法
CN114386510A (zh) 一种电力***量测错误辨识方法及***
CN115879040A (zh) 电容器的故障预测方法、装置、电子设备及存储介质
CN113051741B (zh) 一种高压直流输电***元件可靠性参数求解方法及装置
CN112231981B (zh) 大规模电磁暂态仿真算例的建立方法
CN113036716B (zh) 发电机绕组故障保护方法及装置
CN111612232B (zh) 基于梯度下降的配电网线路重跳概率预测优化方法和装置
CN112180314A (zh) 一种抗干扰自校正gis电子式互感器采集方法和装置
CN112421616A (zh) 一种电网主网中设备参数的辨识方法及***
CN114123164B (zh) 一种电力***短路容量计算方法、装置及终端设备
CN110661261A (zh) 一种用于低频振荡在线分析的pmu数据处理方法及装置
CN104836228B (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