CN116581754A - 一种配电网实时潮流计算方法、装置、设备及介质 - Google Patents
一种配电网实时潮流计算方法、装置、设备及介质 Download PDFInfo
- Publication number
- CN116581754A CN116581754A CN202310848301.9A CN202310848301A CN116581754A CN 116581754 A CN116581754 A CN 116581754A CN 202310848301 A CN202310848301 A CN 202310848301A CN 116581754 A CN116581754 A CN 116581754A
- Authority
- CN
- China
- Prior art keywords
- line
- data
- calculating
- node
- distribution network
- 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
- 238000009826 distribution Methods 0.000 title claims abstract description 107
- 238000000034 method Methods 0.000 title claims abstract description 71
- 238000005259 measurement Methods 0.000 claims abstract description 96
- 238000010586 diagram Methods 0.000 claims abstract description 45
- 238000004364 calculation method Methods 0.000 claims abstract description 35
- 230000008878 coupling Effects 0.000 claims abstract description 12
- 238000010168 coupling process Methods 0.000 claims abstract description 12
- 238000005859 coupling reaction Methods 0.000 claims abstract description 12
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 11
- 238000005070 sampling Methods 0.000 claims description 24
- 230000008569 process Effects 0.000 claims description 17
- 238000003860 storage Methods 0.000 claims description 12
- 238000004590 computer program Methods 0.000 claims description 10
- 238000012937 correction Methods 0.000 claims description 7
- 238000010276 construction Methods 0.000 claims description 4
- 238000012795 verification Methods 0.000 claims description 4
- 230000004931 aggregating effect Effects 0.000 claims description 3
- 230000000295 complement effect Effects 0.000 claims description 3
- 238000006467 substitution reaction Methods 0.000 claims description 3
- 239000013589 supplement Substances 0.000 claims description 3
- 238000011144 upstream manufacturing Methods 0.000 claims description 3
- 238000004458 analytical method Methods 0.000 abstract description 4
- 239000010949 copper Substances 0.000 description 8
- 230000006870 function Effects 0.000 description 6
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 5
- 229910052782 aluminium Inorganic materials 0.000 description 5
- 238000004804 winding Methods 0.000 description 5
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 description 4
- 229910052802 copper Inorganic materials 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 239000000463 material Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 229910000838 Al alloy Inorganic materials 0.000 description 1
- 229910000831 Steel Inorganic materials 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 229910003460 diamond Inorganic materials 0.000 description 1
- 239000010432 diamond Substances 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
- 230000009466 transformation Effects 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
-
- 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
- G06F2113/00—Details relating to the application field
- G06F2113/04—Power grid 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
- 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
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Power Engineering (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种配电网实时潮流计算方法、装置、设备及介质，其方法包括获取目标配电网的线路信息；基于线路信息，获取各线路的下属设备台账、拓扑数据以及量测数据；根据拓扑数据构建各线路的拓扑数据图，并通过量测数据对拓扑数据图进行校正，对校正后的拓扑数据图进行简化；根据下属设备台账计算各线路的电气参数，并通过量测数据和简化后的拓扑数据图对电气参数进行校正；对校正后的电气参数、简化后的拓扑数据图以及量测数据进行耦合关联和差异化状态估计；基于差异化状态估计结果，采用前推回代算法计算潮流值；本发明能够实现配电网实时潮流计算，支撑配电网精准规划决策、精益运行分析。
Description
技术领域
本发明涉及一种配电网实时潮流计算方法、装置、设备及介质，属于电力系统技术领域。
背景技术
电力系统潮流计算是分析电力系统稳态运行情况的一种基本电气计算，是电力系统安全、经济分析和实时控制与调度的基础，大电网潮流计算（10kV以上）已伴随调度自动化系统（D5000）的推广而广泛应用，而配电网由于数据分散在多个业务系统、数据质量不高、量测监测点不足等问题，导致配电网实时潮流尚未实现在线工程化计算。
配电网理论潮流计算方法主要包括前推回代法、PQ分解法、牛顿法等，理论方法已相对成熟，但在实际应用中，主要存在以下问题：
1、理论计算研究与配电网实际工程情况相脱离
配电网潮流计算理论研究中，基于PMU(phasor measurement unit 相量测量装置)的配网潮流计算方法，假定配电网拓扑、参数及量测已知基础上，构建各类求解模型，但其假定已知输入实际配电网并无法提供，而在实际中，由于配电网点多面广，加装PMU投资量巨大，难以实现。
2、长期以来配网潮流计算依靠线下建模，难以大规模计算
长期以来，配电网计算程序尚未打通与现有设备管理系统PMS，调度自动化D5000，配电自动化D5200以及用电信息采集系统之间的联系，在部分地区，配电线路建模仍依赖于人工，工作量巨大，难以大规模计算。
3、配电自动化尚未全面覆盖，不能满足计算需求
部分做法参考调度自动化对于主网计算方法，试图依赖配电自动化数据进行配网潮流计算。但由于配电自动化FTU（Feeder Terminal Unit，馈线终端单元）、DTU(Distribution Terminal Unit，配电终端单元)、TTU(Transformer Terminal Unit，变压器终端单元)等设备尚未实现全量覆盖，利用这些设备采集的电气量无法满足配电网可观性要求，导致配电网潮流不可计算。
4、配电网缺乏设备实测电气参数
现有方法假设已知配网设备电气参数，主网设备如主变、输电线路等均有实测的电阻电抗等电气参数，而配电网设备如架空线、电缆、配电变压器等均无实测电气参数，导致配电网潮流无法计算。
5、配电网拓扑模型不准确
现有方法假设拓扑模型准确，或仅存在拓扑联通性等问题，然而配网拓扑还存在线变关系不一致，开关状态错误等拓扑异常情况，导致配网潮流计算失真。
6、采集设备一次采集成功率较低，无法满足实时计算需求
现有方法并未考虑实时计算情形下，配电网采集设备一次采集成功率较低，导致配网实时量测缺失，进而导致配网潮流无法计算。
发明内容
本发明的目的在于克服现有技术中的不足，提供一种配电网实时潮流计算方法、装置、设备及介质，解决配电网实时潮流计算的技术问题，从而支撑配电网精准规划决策、精益运行分析。
为达到上述目的，本发明是采用下述技术方案实现的：
第一方面，本发明提供了一种配电网实时潮流计算方法，包括：
获取目标配电网的线路信息，所述线路信息包括各线路的线路ID和额定电压；
获取各线路的下属设备台账、拓扑数据以及量测数据；
根据拓扑数据构建各线路的拓扑数据图，并通过量测数据对拓扑数据图进行校正，对校正后的拓扑数据图进行简化；
根据下属设备台账计算各线路的电气参数，并通过量测数据和简化后的拓扑数据图对电气参数进行校正；
对校正后的电气参数、简化后的拓扑数据图以及量测数据进行耦合关联和差异化状态估计；
基于差异化状态估计结果，采用前推回代算法计算各线路的潮流值。
可选的，所述量测数据包括采集时刻和历史时段的量测数据，所述历史时段为采集时刻的前一天；若采集时刻的量测数据存在缺失，则根据当前采集时刻的上一采样时刻、下一采样时刻以及历史时段中同一采样时刻的量测值计算当前采集时刻的补足量测数据：
；
；
；
式中，为当前采集时刻的补足量测数据，/>为当前采集时刻的上一时刻、下一时刻的量测数据，/>为历史时段中同一采样时刻的量测值，/>为历史时段中各采样时刻的量测值平均值。
可选的，所述拓扑数据包括各配网设备的设备ID、所属节点或支线、首端关系节点、末端关系节点以及设备标签；所述设备标签包括T节点、负荷节点、开关、支线，所述开关包括出线开关和联络开关；
所述拓扑数据图的构建包括：以出线开关为起始节点，以负荷节点或联络开关为终止节点，以所属节点或支线、首端关系节点、末端关系节点为依据，在起始节点和终止节点之间遍历生成拓扑数据图。
可选的，所述通过量测数据对拓扑数据图进行校正包括：
根据线路的量测数据获取拓扑数据图中各节点的电压曲线；
将起始节点作为基准节点，将除起始节点以外的节点作为校验节点；
基于获取的电压曲线根据电压相似性原理计算基准节点和各校验节点之间的皮尔逊相关系数；
若，则相应的校验节点拓扑校验正确，无需修正；
若，则相应的校验节点拓扑校验错误，需要修正；所述修正为将校验节点从拓扑数据图中剔除。
可选的，所述对校正后的拓扑数据图进行简化包括：
去除拓扑数据图中不存在负荷节点的支线；
将T节点和开关上游的开关、支线、负荷节点剔除；
聚合相邻两个节点之间的支线；
相邻的两个节点之间仅保留一个开关。
可选的，所述下属设备台账包括各配网设备的设备ID及其设备型号；
所述根据下属设备台账计算各线路的电气参数包括：根据各设备型号与预构建的配网设备典型参数库进行匹配，获取各设备ID的典型参数；根据各典型参数计算各设备ID的电气参数。
可选的，所述根据各典型参数计算各设备ID的电气参数包括：
计算配电变压器的等效阻抗和对地导纳：
；
；
；
；
式中，分别为配电变压器的等值电阻和等值电抗，/>分别为对地电导和对地电纳，/>分别为配电变压器的额定电压、额定功率、短路损耗、短路电压百分比、空载损耗功率以及空载电流百分比；
计算架空线和电缆线的等效阻抗和对地电容：
；
式中，分别为架空线的等值电阻、等值电抗和对地电容，分别为架空线的长度、单位长度的等值电阻、等值电抗和对地电容；分别为电缆线的等值电阻、等值电抗和对地电容，/>分别为电缆线的长度、单位长度的等值电阻、等值电抗和对地电容。
可选的，所述通过量测数据和简化后的拓扑数据图对电气参数进行校正包括：
根据线路的量测数据获取简化后的拓扑数据图的各节点的潮流值；
根据各节点的潮流值计算各设备ID对应的电阻参数和电抗参数；
根据各设备ID对应的电阻参数和等值电阻计算电阻相对偏差，若电阻相对偏差大于预设电阻偏差阈值，则等值电阻校验错误，将等值电阻修正为电阻参数；若电阻相对偏差小于等于预设电阻偏差阈值，则等值电阻校验正确；
根据各设备ID对应的电抗参数和等值电抗计算电抗相对偏差，若电抗相对偏差大于预设电抗偏差阈值，则等值电抗校验错误，将等值电抗修正为电抗参数；若电抗相对偏差小于等于预设电抗偏差阈值，则等值电抗校验正确。
可选的，所述差异化状态估计包括：
采样加权最小二乘法计算量测数据的差异化状态估计值：
；
式中，为采样时刻的量测数据，/>为电压幅值相角，/>为非线性量测函数，为量测数据的置信度；
；
式中，、/>为历史时段中各采样时刻的量测值平均值和标准差。
可选的，所述采用前推回代算法计算各线路的潮流值包括：
获取线路末端电压和功率的差异化状态估计，通过逐段推算出首端电压和功率的计算量；
计算首端电压的计算量和差异化状态估计的差值误差，若差值误差小于等于预设误差阈值，则记录首端电压的计算量，根据首端电压的计算量计算各节点的潮流值；
若差值误差大于预设误差阈值，则获取线路首端电压的差异化状态估计，结合首端功率的计算量，通过逐段推算出末端电压和功率的计算量；
将末端电压和功率的计算量作为线路末端电压和功率的差异化状态估计，重复上述步骤。
第二方面，本发明提供了一种配电网实时潮流计算装置，所述装置包括：
信息获取模块，用于获取目标配电网的线路信息，所述线路信息包括各线路的线路ID和额定电压；
数据获取模块，用于获取各线路的下属设备台账、拓扑数据以及量测数据；
拓扑计算模块，用于根据拓扑数据构建各线路的拓扑数据图，并通过量测数据对拓扑数据图进行校正，对校正后的拓扑数据图进行简化；
电气计算模块，用于根据下属设备台账计算各线路的电气参数，并通过量测数据和简化后的拓扑数据图对电气参数进行校正；
耦合估计模块，用于对校正后的电气参数、简化后的拓扑数据图以及量测数据进行耦合关联和差异化状态估计；
潮流计算模块，用于基于差异化状态估计结果，采用前推回代算法计算各线路的潮流值。
第三方面，本发明提供了一种电子设备，包括处理器及存储介质；
所述存储介质用于存储指令；
所述处理器用于根据所述指令进行操作以执行根据上述方法的步骤。
第四方面，本发明提供了一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现上述方法的步骤。
与现有技术相比，本发明所达到的有益效果：
本发明提供的一种配电网实时潮流计算方法、装置、设备及介质，通过获取各线路的下属设备台账、拓扑数据以及量测数据，生成对校正后的电气参数、简化后的拓扑数据图以及量测数据，并进行耦合关联和差异化状态估计；基于差异化状态估计结果，采用前推回代算法计算潮流值；可解决长期以来采集点不足、实时量测缺失、无实测参数、拓扑不准等问题导致的配电网实时潮流不可算及算不准的难题，支撑配电网精准规划决策、精益运行分析。
附图说明
图1是本发明实施例一提供的一种配电网实时潮流计算方法的流程图。
图2是本发明实施例一提供的一种构建好的拓扑数据图示例的示意图。
图3是本发明实施例一提供的图2简化后的拓扑数据图示例的示意图。
具体实施方式
下面结合附图对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案，而不能以此来限制本发明的保护范围。
实施例一：
如图1所示，本发明提供了一种配电网实时潮流计算方法，包括以下步骤：
S1、获取目标配电网的线路信息，线路信息包括各线路的线路ID和额定电压。
S2、获取各线路的下属设备台账、拓扑数据以及量测数据。
S3、根据拓扑数据构建各线路的拓扑数据图，并通过量测数据对拓扑数据图进行校正，对校正后的拓扑数据图进行简化。
S4、根据下属设备台账计算各线路的电气参数，并通过量测数据和简化后的拓扑数据图对电气参数进行校正。
S5、对校正后的电气参数、简化后的拓扑数据图以及量测数据进行耦合关联和差异化状态估计。
S6、基于差异化状态估计结果，采用前推回代算法计算各线路的潮流值。
具体的：
1、量测数据包括采集时刻和历史时段的量测数据，历史时段为采集时刻的前一天；若采集时刻的量测数据存在缺失，则根据当前采集时刻的上一采样时刻、下一采样时刻以及历史时段中同一采样时刻的量测值计算当前采集时刻的补足量测数据：
；
；
；
式中，为当前采集时刻的补足量测数据，/>为当前采集时刻的上一时刻、下一时刻的量测数据，/>为历史时段中同一采样时刻的量测值，/>为历史时段中各采样时刻的量测值平均值。
在本实施例中，量测数据的来源有3个：一来源于调度自动化系统的配网线路出线开关量测，二来源于用电信息采集系统中配电变压器关口电表量测，三来源于配电自动化系统中柱上开关、环网柜、开闭所内开关FTU、DTU等配电自动化装置量测。
部分量测数据如表1-3所示：
表1：配网线路出线开关量测
| 时间 | 有功功率 | 无功功率 | 电流 | 线电压 | A相电压 | B相电压 | C相电压 |
| 2023-2-19 00:00:00 | 0.8555 | -0.1097 | 49.2195 | 10.4061 | 5.9848 | 6.0668 | 6.02 |
| 2023-2-19 00:15:00 | 0.8555 | -0.1097 | 48.868 | 10.3593 | 5.9438 | 6.0258 | 6.0083 |
| 2023-2-19 00:30:00 | 0.8336 | -0.1097 | 47.4617 | 10.3358 | 5.9438 | 6.0258 | 6.0083 |
| 2023-2-19 00:45:00 | 0.7824 | -0.1097 | 43.5944 | 10.3534 | 5.9321 | 6.0317 | 6.0024 |
| 2023-2-19 01:00:00 | 0.7458 | -0.1097 | 43.946 | 10.3534 | 5.9321 | 6.0317 | 6.0024 |
| 2023-2-19 01:15:00 | 0.7093 | -0.1097 | 41.8366 | 10.3651 | 5.9321 | 6.0434 | 5.9907 |
表2：配电变压器关口电表量测
| ID | 时间 | 电流 | A相电压 | B相电压 | C相电压 | 有功功率 | 无功功率 |
| 14000976371200 | 2023-2-19 00:00:00 | 0.177 | 229.3 | 229.5 | 230 | 66.336 | -6.192 |
| 14000976371200 | 2023-2-19 00:15:00 | 0.159 | 227.8 | 227.9 | 228.4 | 66.24 | -5.376 |
| 14000976371200 | 2023-2-19 00:30:00 | 0.152 | 227.7 | 227.9 | 228.3 | 64.704 | -5.952 |
| 14000976371200 | 2023-2-19 00:45:00 | 0.13 | 228 | 228.1 | 228.5 | 56.88 | -5.424 |
| 14000976371200 | 2023-2-19 01:00:00 | 0.132 | 227.8 | 227.9 | 228.4 | 53.04 | -5.376 |
| 14000976371200 | 2023-2-19 01:15:00 | 0.112 | 228.1 | 228.2 | 228.6 | 47.472 | -5.616 |
| 14000976371200 | 2023-2-19 01:30:00 | 0.128 | 228 | 228.2 | 228.6 | 47.952 | -5.232 |
| 14000976371200 | 2023-2-19 01:45:00 | 0.11 | 228.2 | 228.4 | 228.8 | 48.48 | -4.56 |
表3：配电自动化装置量测
| ID | 时间 | 有功功率 | 无功功率 | A相电流 | B相电流 | C相电流 |
| 14000976403846 | 2023-02-19 00:00:00 | 0.36 | -0.024 | 19.800001 | 19.800001 | 19.800001 |
| 14000976403846 | 2023-02-19 00:15:00 | 0.336 | -0.012 | 19.800001 | 19.2 | 18.6 |
| 14000976403846 | 2023-02-19 00:30:00 | 0.36 | -0.024 | 20.400002 | 20.400002 | 19.2 |
| 14000976403846 | 2023-02-19 00:45:00 | 0.324 | -0.012 | 18 | 18 | 18.6 |
| 14000976403846 | 2023-02-19 01:00:00 | 0.3 | -0.024 | 18.6 | 17.400002 | 16.800001 |
| 14000976403846 | 2023-02-19 01:15:00 | 0.288 | -0.012 | 17.400002 | 17.400002 | 16.2 |
| 14000976403846 | 2023-02-19 01:30:00 | 0.288 | 0 | 16.800001 | 15.6 | 15.000001 |
| 14000976403846 | 2023-02-19 01:45:00 | 0.288 | -0.024 | 16.800001 | 15.6 | 15.6 |
| 14000976403846 | 2023-02-19 02:00:00 | 0.264 | -0.024 | 15.000001 | 14.400001 | 15.000001 |
2、拓扑数据包括各配网设备的设备ID、所属节点或支线、首端关系节点、末端关系节点以及设备标签；设备标签包括T节点（母线、分支箱等）、负荷节点（普通负荷）、开关、支线（电缆、终端头、导线、塔杆等），开关包括出线开关和联络开关；通常通过获取线路cim文件地址，读取cim文件内容，解析cim文件，整理拓扑连接结构。
拓扑数据如表4所示：
表4：拓扑数据
| 设备ID | 所属节点或支线 | 首端关系节点 | 末端关系节点 | 设备标签 |
| SBID000000654AA5E41AB148158CC50EC048DE233A | T1303 | - | -09544- | T节点 |
| SBID00000071BBA6450A90442A86E26D362BCF6A8F | 03212 | 09533 | 09534 | 开关 |
| 14000978947835 | 03215 | 09536 | 09537 | 开关 |
| 14000976257656 | 03213 | 09534 | 09535 | 负荷节点 |
| 14000976257672 | 03214 | 09535 | 09536 | 负荷节点 |
| d37e0e20-e79f-4fbb-9441-7b41fbf45320 | T1304 | 09537 | --- | T节点 |
| 14000976257633 | 09533 | T1303 | 03212 | 支线 |
| 14000976257934 | 09534 | 03212 | 03213 | 支线 |
拓扑数据图的构建包括：以出线开关为起始节点，以负荷节点或联络开关为终止节点，以所属节点或支线、首端关系节点、末端关系节点为依据，在起始节点和终止节点之间遍历生成拓扑数据图；
如图2所示，本实施例提供了一种构建好的拓扑数据图示例，图中六边形为出线开关、菱形为联络开关、矩形为其他开关、方形为负荷节点、圆形为T节点，箭头为潮流方向，省略了支线末端的负荷节点。
3、通过量测数据对拓扑数据图进行校正包括：
根据线路的量测数据获取拓扑数据图中各节点的电压曲线；
将起始节点作为基准节点，将除起始节点以外的节点作为校验节点；
基于获取的电压曲线根据电压相似性原理计算基准节点和各校验节点之间的皮尔逊相关系数；
若，则相应的校验节点拓扑校验正确，无需修正；
若，则相应的校验节点拓扑校验错误，需要修正；所述修正为将校验节点从拓扑数据图中剔除；在剔除后，可进一步的将剔除的校验节点放入其他线路的拓扑数据图构建过程中，并再次进行拓扑校验。
4、对校正后的拓扑数据图进行简化包括：
去除拓扑数据图中不存在负荷节点的支线；
将T节点和开关上游的开关、支线、负荷节点剔除；
聚合相邻两个节点之间的支线；
相邻的两个节点之间仅保留一个开关；
如图3所示，为图2简化后的拓扑数据图示例，将主支线上的两个负荷节点和一个开关进行了收缩简化。
5、下属设备台账包括各配网设备的设备ID及其设备型号；
根据下属设备台账计算各线路的电气参数包括：根据各设备型号与预构建的配网设备典型参数库进行匹配，获取各设备ID的典型参数；根据各典型参数计算各设备ID的电气参数。
其中，配网设备典型参数库为根据国标规定为明确设备型号与典型参数的对应关系所构建的库，其中设备涉及配电变压器、架空线、电缆线三类，表5-7分别为配电变压器、架空线、电缆线典型参数库的部分数据：
表5：配电变压器典型参数库示例
| 型号 | 变压器类型 | 调压方式 | 额定容量(MVA) | 空载损耗(kW) | 空载电流(%) | 短路损耗(kW) | 短路阻抗(%) |
| SCB11-630/20(10) | 双绕组 | 无励磁调压 | 0.63 | 1.28 | 0.54 | 7.34 | 6 |
| SGB10-800/10 | 双绕组 | 有载调压 | 0.8 | 0.85 | 1.52 | 6.99 | 6 |
| SCB-630/10 | 双绕组 | 有载调压 | 0.63 | 0.85 | 0.965 | 5.306666667 | 4 |
| ... | |||||||
| SCB13 | 双绕组 | 有载调压 | 1.25 | 1.5 | 0.85 | 8.746666667 | 6 |
| S13-M-630/10 | 双绕组 | 有载调压 | 0.63 | 0.6 | 0.57 | 6.2 | 4.5 |
表6：架空线典型参数库示例
| 型号 | 导线类型 | 导线材质 | 截面积(mm²) | 电阻(Ω) | 电抗(Ω) | 载流量(A) |
| JKLYJ-25 | 架空线路型号 | 铝 | 25 | 1.2634 | 0.4096 | 97 |
| JKLYJ-240/10 | 架空线路型号 | 铝 | 240 | 0.13 | 0.356 | 500 |
| .... | ||||||
| JKLYJ-10-150 | 架空线路型号 | 铝 | 150 | 0.26368 | 0.295839341 | 285 |
| JNRLH3/LBY-300/55 | 架空线路型号 | 铝包钢芯铝合金 | 300 | 0.141184 | 0.272558996 | 565 |
| JKLYJ-150 | 架空线路型号 | 铝 | 150 | 0.26368 | 0.295839341 | 285 |
表7：电缆线典型参数库示例
| 型号 | 导线类型 | 导线材质 | 截面积(mm²) | 电阻(Ω) | 电抗(Ω) | 载流量(A) |
| ZC-YJV22-18/30-3*400 | 电缆线路型号 | 铝 | 400 | 0.078 | 0.06 | 520 |
| ZRYJV22-18/20-3*240 | 电缆线路型号 | 铜 | 240 | 0.0771 | 0.06 | 8.7 |
| .... | ||||||
| ZA-YJV22-8.7/15-3*400 | 电缆线路型号 | 铜 | 400 | 0.0463 | 0.06 | 333 |
| YJV22-3*400 | 电缆线路型号 | 铜 | 400 | 0.06016 | 0.078240474 | 610 |
| YJV22-3*240 | 电缆线路型号 | 铜 | 240 | 0.096512 | 0.082846541 | 500 |
配电变压器、架空线、电缆线台账信息的部分数据，如表8-10所示：
表8：配电变压器台账信息
| ID | 型号 | 额定容量 | 联结组标号 |
| 14000976371200 | SGB10-800/10 | 800 | Dyn11 |
| 14000976375583 | SGB10-800/10 | 800 | Dyn11 |
| .... | |||
| 14000976383990 | SCB-630/10 | 630 | Dyn11 |
| 4cace66e4e8a06deec54f3d0c401554cac22d01054 | S13-M-630/10 | 630 | Dyn11 |
| d9f40ad119ff80808163d0e8100163d9f34e3f4bc0 | SCB-630/10 | 630 | Dyn11 |
| dd14698574ff80808163d0e8100163dd139bb468eb | SCB-630/10 | 630 | Dyn11 |
表9：架空线台账信息
| ID | 型号 | 长度 | 截面积 |
| SBID000000BAF003482C4B4339B63F889F9A3FF136 | JKLYJ-150 | 30 | 150 |
| SBID0000003CBC0E8559C24765957EF734493211ED | JKLYJ-150 | 69.05 | 150 |
| SBID000000C08246A7DCA2452FA5DB3605F4C02298 | JKLYJ-150 | 39.56 | 150 |
| a7c3e44f-c194-4371-92ab-17c367fd77fc | ZS-JKL/LHYJ-150 | 42.5975 | 150 |
| .... | |||
| SBID0000008439E1FCAE9145A3AF2F91125A599F60 | JKLYJ-150 | 51.67 | 150 |
| 2741d66d-63ed-45ff-b207-25d7d53b4460 | JKLYJ-10-150 | 50 | 150 |
| SBID000000C4F94BDA3B6B41FABDE490A82D01ACFD | ZS-JKL/LHYJ-150 | 70 | 150 |
表10：电缆线台账信息
| ID | 型号 | 长度 | 截面积 |
| SBID000000654AA5E41AB148158CC50EC048DE233A | YJV22-3*400 | 11.43 | 400 |
| SBID00000071BBA6450A90442A86E26D362BCF6A8F | YJV22-3*400 | 100 | 400 |
| 14000978947835 | YJV22-3*240 | 375 | 240 |
| 14000976257696 | YJV22-3*70 | 439 | 70 |
| .... | |||
| d37e0e20-e79f-4fbb-9441-7b41fbf45320 | YJV22-3*400 | 122 | 400 |
| 14001000168255 | YJV22-3*95 | 25 | 95 |
匹配获取的各设备ID的典型参数，如表11-13所示：
表11：配电变压器典型参数
| ID | 型号 | 容量 | 空载电流 | 空载损耗 | 短路损耗 | 短路电压百分比 |
| 14000976371200 | SGB10-800/10 | 800 | 0.85 | 1.52 | 6.99 | 6 |
| 14000976375583 | SGB10-800/10 | 800 | 0.85 | 1.52 | 6.99 | 6 |
| .... | ||||||
| 14000976383990 | SCB-630/10 | 630 | 0.85 | 0.965 | 5.306666667 | 4 |
| 4cace66e4e8a06deec54f3d0c401554cac22d01054 | S13-M-630/10 | 630 | 0.6 | 0.57 | 6.2 | 4.5 |
| d9f40ad119ff80808163d0e8100163d9f34e3f4bc0 | SCB-630/10 | 630 | 0.85 | 0.965 | 5.306666667 | 4 |
| dd14698574ff80808163d0e8100163dd139bb468eb | SCB-630/10 | 630 | 0.85 | 0.965 | 5.306666667 | 4 |
表12：架空线典型参数
| ID | 型号 | 每千米电阻(Ω) | 每千米电抗(Ω) | 每千米对地电容(nf) |
| SBID000000BAF003482C4B4339B63F889F9A3FF136 | JKLYJ-150 | 0.26368 | 0.295839341 | 12.44552089 |
| SBID0000003CBC0E8559C24765957EF734493211ED | JKLYJ-150 | 0.26368 | 0.295839341 | 12.44552089 |
| SBID000000C08246A7DCA2452FA5DB3605F4C02298 | JKLYJ-150 | 0.26368 | 0.295839341 | 12.44552089 |
| a7c3e44f-c194-4371-92ab-17c367fd77fc | ZS-JKL/LHYJ-150 | 0.26368 | 0.295839341 | 12.44552089 |
| .... | ||||
| SBID0000008439E1FCAE9145A3AF2F91125A599F60 | JKLYJ-150 | 0.26368 | 0.295839341 | 12.44552089 |
| 2741d66d-63ed-45ff-b207-25d7d53b4460 | JKLYJ-10-150 | 0.26368 | 0.295839341 | 12.44552089 |
| SBID000000C4F94BDA3B6B41FABDE490A82D01ACFD | ZS-JKL/LHYJ-150 | 0.26368 | 0.295839341 | 12.44552089 |
表13：电缆线典型参数
| ID | 型号 | 每千米电阻(Ω) | 每千米电抗(Ω) | 每千米对地电容(nf) |
| SBID000000654AA5E41AB148158CC50EC048DE233A | YJV22-3*400 | 0.06016 | 0.078240474 | 531.1 |
| SBID00000071BBA6450A90442A86E26D362BCF6A8F | YJV22-3*400 | 0.06016 | 0.078240474 | 531.1 |
| 14000978947835 | YJV22-3*240 | 0.096512 | 0.082846541 | 429 |
| 14000976257696 | YJV22-3*70 | 0.34304 | 0.097716721 | 269.3 |
| .... | ||||
| d37e0e20-e79f-4fbb-9441-7b41fbf45320 | YJV22-3*400 | 0.06016 | 0.078240474 | 531.1 |
| 14001000168255 | YJV22-3*95 | 0.24704 | 0.0935021 | 299.9 |
6、根据各典型参数计算各设备ID的电气参数包括：
计算配电变压器的等效阻抗和对地导纳：
；
；
；
；
式中，分别为配电变压器的等值电阻和等值电抗，/>分别为对地电导和对地电纳，/>分别为配电变压器的额定电压、额定功率、短路损耗、短路电压百分比、空载损耗功率以及空载电流百分比；
计算架空线和电缆线的等效阻抗和对地电容：
；
式中，分别为架空线的等值电阻、等值电抗和对地电容，/>分别为架空线的长度、单位长度的等值电阻、等值电抗和对地电容；/>分别为电缆线的等值电阻、等值电抗和对地电容，/>分别为电缆线的长度、单位长度的等值电阻、等值电抗和对地电容。
7、通过量测数据和简化后的拓扑数据图对电气参数进行校正包括：
根据线路的量测数据获取简化后的拓扑数据图的各节点的潮流值；
根据各节点的潮流值计算各设备ID对应的电阻参数和电抗参数；
根据各设备ID对应的电阻参数和等值电阻计算电阻相对偏差，若电阻相对偏差大于预设电阻偏差阈值，则等值电阻校验错误，将等值电阻修正为电阻参数；若电阻相对偏差小于等于预设电阻偏差阈值，则等值电阻校验正确；
根据各设备ID对应的电抗参数和等值电抗计算电抗相对偏差，若电抗相对偏差大于预设电抗偏差阈值，则等值电抗校验错误，将等值电抗修正为电抗参数；若电抗相对偏差小于等于预设电抗偏差阈值，则等值电抗校验正确。
8、差异化状态估计包括：
采样加权最小二乘法计算量测数据的差异化状态估计值：
；
式中，为采样时刻的量测数据，/>为电压幅值相角，/>为非线性量测函数，为量测数据的置信度；
；
式中，、/>为历史时段中各采样时刻的量测值平均值和标准差。
以配电变压器的有功功率为例，表14为差异化状态估计值：
表14：差异化状态估计值
| ID | 类型 | 断面时间 | 状态估计 | 量测评估置信度 | 差异化状态估计值 |
| 14000976371200 | 配电变压器 | 2023-02-20 13:45:00 | 70.7049 | 0.9933 | 70.2312 |
| 14000976386482 | 配电变压器 | 2023-02-20 13:45:00 | 30.8657 | 0.9814 | 30.2916 |
| 14000976375583 | 配电变压器 | 2023-02-20 13:45:00 | 68.6982 | 0.8721 | 59.9117 |
| 14000976387987 | 配电变压器 | 2023-02-20 13:45:00 | 40.6464 | 0.9432 | 38.3377 |
| .... | |||||
| 143064475 | 线路出线开关 | 2023-02-20 13:45:00 | 518.5269 | 0.9436 | 489.282012 |
9、采用前推回代算法计算各线路的潮流值包括：
获取线路末端电压和功率的差异化状态估计，通过逐段推算出首端电压和功率的计算量；
计算首端电压的计算量和差异化状态估计的差值误差，若差值误差小于等于预设误差阈值，则记录首端电压的计算量，根据首端电压的计算量计算各节点的潮流值；
若差值误差大于预设误差阈值，则获取线路首端电压的差异化状态估计，结合首端功率的计算量，通过逐段推算出末端电压和功率的计算量；
将末端电压和功率的计算量作为线路末端电压和功率的差异化状态估计，重复上述步骤。
以配电变压器的有功功率为例，表15为潮流值：
表15：潮流值
| ID | 类型 | 断面时间 | 差异化状态估计 | 潮流值 |
| 14000976371200 | 配电变压器 | 2023-02-20 13:45:00 | 70.2312 | 69.2310 |
| 14000976386482 | 配电变压器 | 2023-02-20 13:45:00 | 30.2916 | 30.1820 |
| 14000976375583 | 配电变压器 | 2023-02-20 13:45:00 | 59.9117 | 58.8901 |
| 14000976387987 | 配电变压器 | 2023-02-20 13:45:00 | 38.3377 | 38.1128 |
| .... | ||||
| 143064475 | 线路出线开关 | 2023-02-20 13:45:00 | 489.2820 | 509.7064 |
实施例二：
本发明实施例提供了一种配电网实时潮流计算装置，所述装置包括：
信息获取模块，用于获取目标配电网的线路信息，所述线路信息包括各线路的线路ID和额定电压；
数据获取模块，用于获取各线路的下属设备台账、拓扑数据以及量测数据；
拓扑计算模块，用于根据拓扑数据构建各线路的拓扑数据图，并通过量测数据对拓扑数据图进行校正，对校正后的拓扑数据图进行简化；
电气计算模块，用于根据下属设备台账计算各线路的电气参数，并通过量测数据和简化后的拓扑数据图对电气参数进行校正；
耦合估计模块，用于对校正后的电气参数、简化后的拓扑数据图以及量测数据进行耦合关联和差异化状态估计；
潮流计算模块，用于基于差异化状态估计结果，采用前推回代算法计算各线路的潮流值。
实施例三：
基于实施例一，本发明实施例提供了一种电子设备，其特征在于，包括处理器及存储介质；
所述存储介质用于存储指令；
所述处理器用于根据所述指令进行操作以执行根据上述方法的步骤。
实施例四：
基于实施例一，本发明实施例提供了一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现上述方法的步骤。
本发明应用于某省份全部41790条配电网线路每15分钟在线自动实时潮流计算，可以通过计算获得全部线路和节点状态的线路占全部线路的比例达95%以上，节点计算值与已有量测值比对偏差在5%以内的比例达到99%。
本领域内的技术人员应明白，本申请的实施例可提供为方法、系统、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质（包括但不限于磁盘存储器、CD-ROM、光学存储器等）上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备（系统）、和计算机程序产品的流程图和／或方框图来描述的。应理解可由计算机程序指令实现流程图和／或方框图中的每一流程和／或方框、以及流程图和／或方框图中的流程和／或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和／或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和／或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和／或方框图一个方框或多个方框中指定的功能的步骤。
以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明技术原理的前提下，还可以做出若干改进和变形，这些改进和变形也应视为本发明的保护范围。
Claims (13)
1.一种配电网实时潮流计算方法，其特征在于，包括：
获取目标配电网的线路信息，所述线路信息包括各线路的线路ID和额定电压；
获取各线路的下属设备台账、拓扑数据以及量测数据；
根据拓扑数据构建各线路的拓扑数据图，并通过量测数据对拓扑数据图进行校正，对校正后的拓扑数据图进行简化；
根据下属设备台账计算各线路的电气参数，并通过量测数据和简化后的拓扑数据图对电气参数进行校正；
对校正后的电气参数、简化后的拓扑数据图以及量测数据进行耦合关联和差异化状态估计；
基于差异化状态估计结果，采用前推回代算法计算各线路的潮流值。
2.根据权利要求1所述的配电网实时潮流计算方法，其特征在于，所述量测数据包括采集时刻和历史时段的量测数据，所述历史时段为采集时刻的前一天；若采集时刻的量测数据存在缺失，则根据当前采集时刻的上一采样时刻、下一采样时刻以及历史时段中同一采样时刻的量测值计算当前采集时刻的补足量测数据：
；
；
；
式中，为当前采集时刻的补足量测数据，/>为当前采集时刻的上一时刻、下一时刻的量测数据，/>为/>的差值绝对值，/>为历史时段中同一采样时刻的量测值，/>为历史时段中各采样时刻的量测值平均值，/>为/>的差值绝对值。
3.根据权利要求1所述的配电网实时潮流计算方法，其特征在于，所述拓扑数据包括各配网设备的设备ID、所属节点或支线、首端关系节点、末端关系节点以及设备标签；所述设备标签包括T节点、负荷节点、开关、支线，所述开关包括出线开关和联络开关；
所述拓扑数据图的构建包括：以出线开关为起始节点，以负荷节点或联络开关为终止节点，以所属节点或支线、首端关系节点、末端关系节点为依据，在起始节点和终止节点之间遍历生成拓扑数据图。
4.根据权利要求3所述的配电网实时潮流计算方法，其特征在于，所述通过量测数据对拓扑数据图进行校正包括：
根据线路的量测数据获取拓扑数据图中各节点的电压曲线；
将起始节点作为基准节点，将除起始节点以外的节点作为校验节点；
基于获取的电压曲线根据电压相似性原理计算基准节点和各校验节点之间的皮尔逊相关系数；
若，则相应的校验节点拓扑校验正确，无需修正；
若，则相应的校验节点拓扑校验错误，需要修正；所述修正为将校验节点从拓扑数据图中剔除。
5.根据权利要求3所述的配电网实时潮流计算方法，其特征在于，所述对校正后的拓扑数据图进行简化包括：
去除拓扑数据图中不存在负荷节点的支线；
将T节点和开关上游的开关、支线、负荷节点剔除；
聚合相邻两个节点之间的支线；
相邻的两个节点之间仅保留一个开关。
6.根据权利要求2所述的配电网实时潮流计算方法，其特征在于，所述下属设备台账包括各配网设备的设备ID及其设备型号；
所述根据下属设备台账计算各线路的电气参数包括：根据各设备型号与预构建的配网设备典型参数库进行匹配，获取各设备ID的典型参数；根据各典型参数计算各设备ID的电气参数。
7.根据权利要求6所述的配电网实时潮流计算方法，其特征在于，所述根据各典型参数计算各设备ID的电气参数包括：
计算配电变压器的等效阻抗和对地导纳：
；
；
；
；
式中，分别为配电变压器的等值电阻和等值电抗，/>分别为对地电导和对地电纳，/>分别为配电变压器的额定电压、额定功率、短路损耗、短路电压百分比、空载损耗功率以及空载电流百分比；
计算架空线和电缆线的等效阻抗和对地电容：
；
式中，分别为架空线的等值电阻、等值电抗和对地电容，/>分别为架空线的长度、单位长度的等值电阻、等值电抗和对地电容；/>分别为电缆线的等值电阻、等值电抗和对地电容，/>分别为电缆线的长度、单位长度的等值电阻、等值电抗和对地电容。
8.根据权利要求7所述的配电网实时潮流计算方法，其特征在于，所述通过量测数据和简化后的拓扑数据图对电气参数进行校正包括：
根据线路的量测数据获取简化后的拓扑数据图的各节点的潮流值；
根据各节点的潮流值计算各设备ID对应的电阻参数和电抗参数；
根据各设备ID对应的电阻参数和等值电阻计算电阻相对偏差，若电阻相对偏差大于预设电阻偏差阈值，则等值电阻校验错误，将等值电阻修正为电阻参数；若电阻相对偏差小于等于预设电阻偏差阈值，则等值电阻校验正确；
根据各设备ID对应的电抗参数和等值电抗计算电抗相对偏差，若电抗相对偏差大于预设电抗偏差阈值，则等值电抗校验错误，将等值电抗修正为电抗参数；若电抗相对偏差小于等于预设电抗偏差阈值，则等值电抗校验正确。
9.根据权利要求1所述的配电网实时潮流计算方法，其特征在于，所述差异化状态估计包括：
采样加权最小二乘法计算量测数据的差异化状态估计值：
；
式中，为采样时刻的量测数据，/>为电压幅值相角，/>为非线性量测函数，/>为量测数据的置信度；
；
式中，、/>为历史时段中各采样时刻的量测值平均值和标准差。
10.根据权利要求1所述的配电网实时潮流计算方法，其特征在于，所述采用前推回代算法计算各线路的潮流值包括：
获取线路末端电压和功率的差异化状态估计，通过逐段推算出首端电压和功率的计算量；
计算首端电压的计算量和差异化状态估计的差值误差，若差值误差小于等于预设误差阈值，则记录首端电压的计算量，根据首端电压的计算量计算各节点的潮流值；
若差值误差大于预设误差阈值，则获取线路首端电压的差异化状态估计，结合首端功率的计算量，通过逐段推算出末端电压和功率的计算量；
将末端电压和功率的计算量作为线路末端电压和功率的差异化状态估计，重复上述步骤。
11.一种配电网实时潮流计算装置，其特征在于，所述装置包括：
信息获取模块，用于获取目标配电网的线路信息，所述线路信息包括各线路的线路ID和额定电压；
数据获取模块，用于获取各线路的下属设备台账、拓扑数据以及量测数据；
拓扑计算模块，用于根据拓扑数据构建各线路的拓扑数据图，并通过量测数据对拓扑数据图进行校正，对校正后的拓扑数据图进行简化；
电气计算模块，用于根据下属设备台账计算各线路的电气参数，并通过量测数据和简化后的拓扑数据图对电气参数进行校正；
耦合估计模块，用于对校正后的电气参数、简化后的拓扑数据图以及量测数据进行耦合关联和差异化状态估计；
潮流计算模块，用于基于差异化状态估计结果，采用前推回代算法计算各线路的潮流值。
12.一种电子设备，其特征在于，包括处理器及存储介质；
所述存储介质用于存储指令；
所述处理器用于根据所述指令进行操作以执行根据权利要求1-10任一项所述方法的步骤。
13.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该程序被处理器执行时实现权利要求1-10任一项所述方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310848301.9A CN116581754B (zh) | 2023-07-12 | 2023-07-12 | 一种配电网实时潮流计算方法、装置、设备及介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310848301.9A CN116581754B (zh) | 2023-07-12 | 2023-07-12 | 一种配电网实时潮流计算方法、装置、设备及介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116581754A true CN116581754A (zh) | 2023-08-11 |
| CN116581754B CN116581754B (zh) | 2023-09-26 |
Family
ID=87543462
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310848301.9A Active CN116581754B (zh) | 2023-07-12 | 2023-07-12 | 一种配电网实时潮流计算方法、装置、设备及介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116581754B (zh) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109167362A (zh) * | 2018-08-08 | 2019-01-08 | 山东大学 | 一种计及电缆热特性的配电网潮流计算方法 |
| CN109980650A (zh) * | 2019-05-21 | 2019-07-05 | 华北水利水电大学 | 一种辐射型配电系统潮流计算方法 |
| WO2023029388A1 (zh) * | 2021-09-06 | 2023-03-09 | 广东电网有限责任公司湛江供电局 | 一种基于cim模型的配电网辅助规划方法及系统 |
-
2023
- 2023-07-12 CN CN202310848301.9A patent/CN116581754B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109167362A (zh) * | 2018-08-08 | 2019-01-08 | 山东大学 | 一种计及电缆热特性的配电网潮流计算方法 |
| CN109980650A (zh) * | 2019-05-21 | 2019-07-05 | 华北水利水电大学 | 一种辐射型配电系统潮流计算方法 |
| WO2023029388A1 (zh) * | 2021-09-06 | 2023-03-09 | 广东电网有限责任公司湛江供电局 | 一种基于cim模型的配电网辅助规划方法及系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116581754B (zh) | 2023-09-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111026927B (zh) | 一种低压台区运行状态智能监测系统 | |
| CN102567603B (zh) | 基于实测拓扑和量测数据自动生成bpa计算文件的方法 | |
| CN109752629B (zh) | 一种电网量测问题智能诊断方法及系统 | |
| CN109034666B (zh) | 一种基于负荷突变的低压电网拓扑大数据辨识系统及方法 | |
| CN104462779B (zh) | 采用快速建模与数据质量分析的理论线损在线计算方法 | |
| CN109286188B (zh) | 一种基于多源数据集的10kV配电网理论线损计算方法 | |
| CN113689119B (zh) | 基于数字孪生的配电网供电可靠性评价方法、设备及介质 | |
| CN112330117B (zh) | 一种配电网规划年供电可靠性评估方法和装置 | |
| CN111680084A (zh) | 一种基于多源数据融合的线损计算方法 | |
| CN105119282A (zh) | 一种电网理论线损在线计算系统及方法 | |
| CN106022641A (zh) | 一种变电站负荷特性的在线分析系统和方法 | |
| CN104834305B (zh) | 基于dms系统的配网自动化终端遥测异常分析系统及方法 | |
| CN115685046A (zh) | 互感器计量异常识别方法、装置、设备及存储介质 | |
| CN109142988B (zh) | 一种基于电能质量监测数据的配网故障定位方法及系统 | |
| CN116581754B (zh) | 一种配电网实时潮流计算方法、装置、设备及介质 | |
| CN111965479A (zh) | 一种基于数据精细化量测比对的线损异常定位方法 | |
| CN111984925A (zh) | 基于回路阻抗的线路异常定位方法、存储介质和计算设备 | |
| CN109494733B (zh) | 一种电力负荷模型的辨识参数优化方法及系统 | |
| CN115712811A (zh) | 基于低压拓扑配电网下的阻抗计算模型 | |
| CN114994401B (zh) | 一种线路线损异常的检测方法及装置 | |
| CN115982902A (zh) | 一种基于地理坐标的电网拓扑关系识别方法及装置 | |
| CN106786498B (zh) | 一种主站-变电站数据协同辨识方法及其装置 | |
| CN104700318A (zh) | 基于设备功率流拓扑的配电网薄弱环节辨识方法 | |
| CN110555244B (zh) | 一种配电网同期线损计算系统及方法 | |
| CN113092936A (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 |