电力系统潮流计算：那些你必须玩转的标准算例

电力系统潮流计算标准算例的数据从3节点到300节点都齐了。 包含IEEE格式、BPA格式、清华格式同时有各个节点的拓扑图搞电力系统分析的兄弟们应该都懂标准算例就是咱的练功房。从3节点到300节点的各种妖魔鬼怪要是没在标准算例里摸爬滚打过出去实战绝对要跪。今儿咱们就来盘一盘这些藏在数据里的门道。先看IEEE格式这个老大哥它家的数据排列跟强迫症似的。举个9节点系统的例子# IEEE数据读取示例 def parse_ieee_bus(line): bus_num int(line[0:4]) bus_type int(line[5:6]) pd float(line[40:49]) qd float(line[49:59]) # 这切片看得眼疼是吧习惯就好... return {num:bus_num, type:bus_type, pd:pd, qd:qd}注意看那严格的列对齐差一个空格整个解析就崩。特别是节点类型字段bus_type0是PQ节点1是PV节点2是平衡节点这个编码规则比女朋友的心思还难记错。BPA格式就野多了自带方言属性。它的支路数据长这样L 母线A 230. 母线B 230. 0.002 0.014 0.14 400.第一眼看到这格式还以为在看电报码。字母L开头代表线路后面跟着的是电压等级和阻抗参数。处理这种数据建议用正则表达式import re bpa_pattern rL\s(\w)\s(\d\.?)\s(\w)\s(\d\.?)\s([\d.])\s([\d.]) match re.match(bpa_pattern, line) if match: from_bus match.group(1) to_bus match.group(3) r float(match.group(5))注意电压等级单位可能隐含在数值里比如230.可能代表230kV这个坑我当年掉进去过。电力系统潮流计算标准算例的数据从3节点到300节点都齐了。 包含IEEE格式、BPA格式、清华格式同时有各个节点的拓扑图清华格式是国产良心特点是用中文关键字。看到这样的数据别慌节点数量118 [节点数据] 1, 平衡节点, 电压1.06 2, PQ节点, 有功1.2, 无功0.3处理这种格式最爽的是可以直接用CSV解析器import csv with open(thupower.dat) as f: reader csv.DictReader(f, fieldnames[节点号,类型,参数]) for row in reader: if row[类型] 平衡节点: # 这里可以玩点骚操作...但要注意中文字段的匹配曾经有兄弟把平衡节点写成平横节点debug了一整夜别问我是怎么知道的。拓扑图这事必须单说。有个冷知识节点编号不连续的系统用networkx画图会出灵异事件。正确的操作姿势是import networkx as nx G nx.Graph() for line in lines_data: G.add_edge(from_node, to_node, weightimpedance) nx.draw(G, with_labelsTrue, node_size500) # 记得调spring_layout参数不然画出来的像毛线团画300节点拓扑时建议关掉节点标签显示否则你会得到一张毕加索风格的杰作。玩转这些算例有个隐藏技巧准备三套解析器互相验证。当IEEE和BPA的计算结果对不上时用清华格式当裁判。有次我卡在节点相角差过大的问题上最后发现是BPA数据里的并联导纳单位搞反了——这种跨格式校验能救命。最后给个忠告别看300节点数据庞大真正要命的反而是小系统。3节点系统的收敛性问题能把牛顿拉夫逊法逼疯这时候就得祭出PQ分解法或者...直接调求解器参数别死磕才是真谛。