CN108631297A - 大型受端电网轨迹灵敏度批处理计算方法 - Google Patents
大型受端电网轨迹灵敏度批处理计算方法 Download PDFInfo
- Publication number
- CN108631297A CN108631297A CN201810273074.0A CN201810273074A CN108631297A CN 108631297 A CN108631297 A CN 108631297A CN 201810273074 A CN201810273074 A CN 201810273074A CN 108631297 A CN108631297 A CN 108631297A
- Authority
- CN
- China
- Prior art keywords
- fault
- file
- stability
- card
- passive compensation
- 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
- 230000035945 sensitivity Effects 0.000 title claims abstract description 36
- 238000004364 calculation method Methods 0.000 title claims abstract description 28
- 238000012545 processing Methods 0.000 title claims abstract description 28
- 238000009434 installation Methods 0.000 claims description 31
- 238000000205 computational method Methods 0.000 claims description 17
- 238000012546 transfer Methods 0.000 claims description 6
- 230000006641 stabilisation Effects 0.000 claims description 4
- 230000005284 excitation Effects 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 238000011105 stabilization Methods 0.000 claims description 3
- 239000003381 stabilizer Substances 0.000 claims description 3
- 238000004519 manufacturing process Methods 0.000 claims description 2
- 238000013515 script Methods 0.000 abstract description 6
- 230000008859 change Effects 0.000 abstract description 4
- 230000000087 stabilizing effect Effects 0.000 abstract 1
- 238000000034 method Methods 0.000 description 10
- 230000008901 benefit Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 4
- 230000005611 electricity Effects 0.000 description 3
- 238000007726 management method Methods 0.000 description 3
- 241001269238 Data Species 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 238000012217 deletion Methods 0.000 description 1
- 230000037430 deletion Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 230000008676 import Effects 0.000 description 1
- 230000035800 maturation Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 230000001052 transient effect 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
-
- 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]
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明提供了一种大型受端电网轨迹灵敏度批处理计算方法，利用MATLAB脚本遍历每种故障与负荷节点的组合，自动修改BPA稳定程序所用到的故障参数卡片与动态无功补偿设备卡片，从而进行稳定计算；并自动对输出数据文件按规律进行命名管理，实现对每种组合电压变化曲线的自动读取，提高计算效率，最终得到受端电网每一个负荷节点上的轨迹灵敏度。现有的手动方法计算大型受端电网所有负荷节点的轨迹灵敏度指标涉及到很大的工程量，而采用本发明的好处在于提供的方法，所有过程均可由MATLAB脚本实现，自动调用BPA稳定程序进行批处理计算，避免重复而枯燥的手动操作，不仅大大增加了计算效率，还避免了出错。
Description
技术领域
本发明涉及电力系统分析技术领域，具体地，涉及一种基于BPA和MATLAB的大型受端电网轨迹灵敏度批处理计算方法。
背景技术
轨迹灵敏度指标可用于衡量电力系统某节点上安装动态无功补偿设备后系统暂态电压稳定的提升程度。对于大型受端电网中的某一负荷节点j，它的轨迹灵敏度TSI可采用下列计算式：
式中：Nl为考察故障总数；WFl为故障Fl的权重，Fl(l＝1,2,...,Nl)为考察故障；Wdi为母线i的权重，n为负荷节点总数；Nk为总的时间采样点数；Wk为tk时刻的权重；Vi,l(tk,Qj0+△Qj(COMP))为节点j安装额定容量为SCOMP的动态无功补偿设备后，在故障Fl下，节点电压Vi在t＝tk时刻的值；△Qj(COMP)为在节点j安装动态无功补偿设备后，该设备释放的无功功率；Vi,l(tk,Qj0)为未安装动态无功补偿设备情况下，在故障Fl下，Vl在t＝tk时刻的值。
可以看出，计算j节点的轨迹灵敏度，关键是需要获得整个电力系统中所有负荷节点在j节点是否安装动态无功补偿设备两种情况下，每种故障发生后的节点电压时域分布。
对于电网节点故障电压的计算和仿真，国内外有不少成熟的软件。中国电力科学院开发的机电暂态软件包PSD-BPA(简称BPA)主要由潮流程序和稳定程序构成，具有计算规模大、计算速度快、数值稳定性好、功能强等特点，已在我国电力系统规划、调度、生产运行及科研部门得到了广泛应用。
利用BPA对大型受端电网的负荷节点进行轨迹灵敏度计算时，需要首先根据网络数据通过潮流程序得到系统所处运行方式下的潮流结果文件(.BSE格式)，然后在稳定程序中导入该潮流结果文件，并手动填写故障参数卡片和动态无功补偿设备参数卡片。在计算受端电网所有负荷节点的轨迹灵敏度从而寻找全网最优动态无功补偿设备安装点时，需要遍历每一种典型故障与补偿设备安装点(也就是负荷节点)的组合，利用BPA稳定程序对每一种组合进行稳定计算，并在每个组合的计算中得到所有负荷节点上的电压变化时域分布。对于大型受端电网，其典型故障数与负荷节点数都具有一定规模，手动修改故障参数卡片和动态无功补偿设备参数卡片将大大增加工作量，降低计算效率；并且稳定程序的输出数据分散在不同的.SWX格式结果文件中，没有统一的命名管理，不利于统计分析。
发明内容
针对现有技术中存在的上述不足，本发明的目的是提供一种基于BPA和MATLAB的大型受端电网轨迹灵敏度批处理计算方法，该方法利用MATLAB脚本遍历每种故障与负荷节点的组合，自动修改BPA稳定程序所用到的故障参数卡片与动态无功补偿设备卡片，从而进行稳定计算；并自动对输出数据文件按规律进行命名管理，实现对每种组合电压变化曲线的自动读取，提高计算效率，最终得到受端电网每一个负荷节点上的轨迹灵敏度。
为了达到上述目的，本发明采用如下技术方案。
一种大型受端电网轨迹灵敏度批处理计算方法，基于BPA和MATLAB实现，包括如下步骤：
S1.利用MATLAB根据BPA潮流文件所提供的受端电网的网络数据分别生成故障集文件、动态无功补偿设备安装点集文件和负荷节点电压输出语句集文件，并记录下故障个数nf和动态无功补偿设备安装点个数nc。
优选地，上述文件均为txt格式。其中，故障集文件包括所有典型故障对应的BPA稳定程序卡片；动态无功补偿设备安装点集文件包括在所要计算轨迹灵敏度指标的所有节点上安装动态无功补偿设备的BPA稳定程序卡片；负荷节点电压输出语句集文件包括全网所有负荷节点上电压时域分布信息的输出指令语句。
S2.根据S1生成的三个文件(故障集文件、动态无功补偿设备安装点集文件和负荷节点电压输出语句集文件)，形成原始的.SWI格式稳定程序输入数据文件，并记录下故障参数卡和动态无功补偿设备参数卡的起始行信息。
优选地，上述原始的.SWI格式稳定程序输入数据文件中，既包括所要遍历所有组合间相同的部分，如控制数据卡、发电机模型卡、励磁系统模型卡、电力系统稳定器(PSS)模型卡、调速器和原动机模型卡、直流控制系统模型卡、负荷模型卡、输出控制卡等；也包括不同组合间产生变化的部分，如故障操作模型卡、动态无功补偿(电力电子)设备模型卡。所有故障操作模型卡和动态无功补偿设备模型卡均需要打上注释符号，以便求解到某种组合时能以去掉注释符号的方式形成相应的.SWI格式稳定程序输入数据文件。
S3.以批处理的方式遍历nf×nc种不同故障搭配不同动态无功补偿设备安装点的组合，根据稳定程序输出结果文件得到每种组合下全网所有负荷节点的电压时域分布(每种组合下全网所有负荷节点的电压时域分布记录在.SWX格式的稳定程序输出结果文件中)。
优选地，包括如下子步骤：
S3.1，定义故障编号变量faultnum和动态无功补偿设备安装点编号变量compnum，利用MATLAB的for循环功能遍历不同故障与安装点的组合。其中，变量faultnum和compnum的格式均为整数。
S3.2，利用MATLAB读取原始的.SWI格式稳定程序输入数据文件，根据故障编号变量faultnum、动态无功补偿设备安装点编号变量compnum以及之前记录下的故障参数卡、补偿设备参数卡起始行信息，删除相应故障操作模型卡和动态无功补偿设备模型卡首列的注释符号，使故障操作模型卡和动态无功补偿设备模型卡生效。
S3.3，利用MATLAB生成每种组合中稳定计算需用到的.SWI格式稳定程序输入数据文件，并统一按照Fault_faultnum_Comp_compnum的规则命名。文件名Fault_faultnum_Comp_compnum的格式为字符串，其中，变量faultnum和compnum的格式由整数通过num2str函数转换成字符串，再与固定字符串Fault_和_Comp_组合。
按照上述规则统一命名有两个好处：一是方便之后调用稳定程序时传递稳定程序输入数据文件(.SWI文件)；二是由于稳定程序输出结果文件(.SWX文件)的名称与稳定程序输入数据文件相同，所以统一命名能方便读取稳定计算结果。
S3.4，在MATLAB里利用dos命令调用BPA稳定程序(swnt.exe)，并将S3.3中生成的稳定程序输入数据文件(.SWI文件)作为传递参数导入BPA稳定程序中，执行稳定计算。在所遍历的每种组合中执行稳定计算将得到相应的全网所有负荷节点的电压时域分布结果，这些结果以所遍历的组合为单位，记录在不同的稳定程序输出结果文件(.SWX文件)中。由于S3.3的操作，所得稳定程序输出结果文件将按Fault_faultnum_Comp_compnum的命名规则排列。
S4.利用MATLAB读取S3中的稳定程序输出结果文件(.SWX文件)，得到当节点j(j＝1,2,...,n)安装动态无功补偿设备后，各负荷节点在不同故障情况下的电压时域分布。并根据权重系数加权求和，得到中间变量R1j。R1j的计算式如下式所示：
式中，Nl为考察故障总数；Nk为总的时间采样点数；WFl为故障Fl的权重；Wdi为母线i的权重；Wk为tk时刻的权重；△Qj(COMP)为在节点j安装动态无功补偿设备后，该设备释放的无功功率；Vi,l(tk,Qj0+△Qj(COMP))为节点j安装额定容量为SCOMP的动态无功补偿设备后，在故障Fl下，节点电压Vi在t＝tk时刻的值。
S5.不考虑安装动态无功补偿设备，以批处理的方式遍历nf种不同故障的情况，得到每种情况下全网所有负荷节点的电压时域分布，并记录在稳定程序输出结果文件中。
优选地，包括如下子步骤：
S5.1，根据S3.1中定义的故障编号变量faultnum，利用MATLAB的for循环功能遍历不同故障情况。
S5.2，利用MATLAB读取原始的.SWI格式稳定程序输入数据文件，根据faultnum以及之前记录下的故障参数卡起始行信息，删除相应故障操作模型卡首列的注释符号，使故障操作模型卡生效。
S5.3，利用MATLAB生成每种故障情况中稳定计算需用到的.SWI格式稳定程序输入数据文件，并统一按照No_Comp_Fault_faultnum的规则命名。文件名No_Comp_Fault_faultnum的格式为字符串，其中，变量faultnum的格式由整数通过num2str函数转换成字符串，再与固定字符串No_Comp_Fault_组合。
S5.4，在MATLAB里利用dos命令调用BPA稳定程序(swnt.exe)，并将S5.3生成的稳定程序输入数据文件(.SWI文件)作为传递参数导入BPA稳定程序中，执行稳定计算。在所遍历的每种故障情况中执行稳定计算将得到相应的全网所有负荷节点的电压时域分布结果，这些结果以所遍历的故障情况为单位，记录在不同的稳定程序输出结果文件(.SWX文件)中。由于S5.3的操作，所得稳定程序输出结果文件将按No_Comp_Fault_faultnum的命名规则排列。
S6.利用MATLAB读取S5中的稳定程序输出结果文件(.SWX文件)，得到不安装动态无功补偿设备时，各负荷节点在不同故障情况下的电压时域分布。并根据权重系数加权求和，得到中间变量R2。R2的计算式如下式所示：
式中，Vi,l(tk,Qj0)为未安装动态无功补偿设备情况下，在故障Fl下，节点电压Vl在t＝tk时刻的值。
S7.由下式计算得到负荷节点j上的轨迹灵敏度TSIj：
TSIj＝R1j-R2。
与现有技术相比，本发明具有如下有益效果：
现有的手动计算大型受端电网所有负荷节点的轨迹灵敏度指标涉及到很大的工程量，而采用本发明提供的方法，所有过程均可由MATLAB脚本实现，自动调用BPA稳定程序进行批处理计算，避免重复而枯燥的手动操作，不仅大大增加了计算效率，还避免了出错。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述，本发明的其它特征、目的和优点将会变得更明显：
图1是基于BPA和MATLAB的大型受端电网轨迹灵敏度批处理计算方法流程图；
图2是MATLAB利用dos命令调用BPA稳定程序，并传递输入数据文件示意图。
具体实施方式
下面对本发明的实施例作详细说明：本实施例在以本发明技术方案为前提下进行实施，给出了详细的实施方式和具体的操作过程。应当指出的是，对本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。
实施例
本实施例提供了一种基于BPA和MATLAB的大型受端电网轨迹灵敏度批处理计算方法，其流程如图1所示，包括下列步骤：
第S1步.利用MATLAB根据潮流文件所提供的受端电网的网络数据生成故障集文件、动态无功补偿设备安装点集文件和负荷节点电压输出语句集文件，并记录下故障个数nf和动态无功补偿设备安装点个数nc。
具体地：
第S1.1步，在MATLAB中读取.dat格式的BPA潮流文件，并识别文件中的支路数据卡，提取支路首末端节点名称及其基准电压信息，据此生成FLT支路三相短路接地故障卡。将所有故障卡依次排列，组成.txt格式的故障集文件。
第S1.2步，在MATLAB中读取.dat格式的BPA潮流文件，并识别文件中的节点数据卡，提取负荷节点的节点名称及其基准电压信息，据此生成动态无功补偿设备安装卡。以STATCOM为例，可生成多张不同安装点的VG/VG+卡。将其依次排列，组成.txt格式动态无功补偿设备安装点集文件。
第S1.3步，根据第S1.2步中提取的负荷节点的节点名称及其基准电压信息，生成BH/B母线输出卡。将其依次排列，组成.txt格式负荷节点电压输出语句集文件。
第S2步.根据第1步生成的三个文件，形成原始的.SWI格式稳定程序输入数据文件，并记录下故障(参数)卡和动态无功补偿设备(参数)卡的起始行。
具体地：利用MATLAB读写文件的功能，将第S1步生成的三个文件中的卡片(故障操作模型卡和动态无功补偿设备模型卡)统一在列首打上注释符号，写入同一文件中，并加上所要遍历的不同组合间相同的部分，如控制数据卡、发电机模型卡、励磁系统模型卡、电力系统稳定器(PSS)模型卡、调速器和原动机模型卡、直流控制系统模型卡、负荷模型卡、输出控制卡等，生成原始的.SWI格式稳定程序输入数据文件。记录下所生成的.SWI格式文件中，FLT卡的起始行startline_F和VG/VG+卡的起始行startline_C。
第S3步.以批处理的方式遍历nf×nc种不同故障搭配不同动态无功补偿设备安装点的组合，得到每种组合下全网所有负荷节点的电压时域分布。
具体地：
第S3.1步，定义故障编号变量“faultnum”和动态无功补偿设备安装点编号变量“compnum”，利用MATLAB的for循环功能遍历faultnum×compnum种不同故障与安装点的组合。
第S3.2步，在每一个循环中，利用MATLAB读取原始的.SWI格式的稳定程序输入数据文件，删除原始文件中第startline_F+faultnum-1行和第startline_C+(compnum-1)*2行首列的注释符号，使相应卡片生效。
第S3.3步，将第S3.2步中修改后的.SWI格式文件统一按照“Fault_faultnum_Comp_compnum”的规则重命名，生成每种组合中稳定计算需用到的.SWI格式稳定程序输入数据文件。由于文件名称为字符串格式，而faultnum和compnum的格式为整数，故需要通过num2str函数将数值变量faultnum和compnum转换成的字符串，再与固定字符串Fault_和_Comp_组合成文件名。
“faultnum”和“compnum”在最初定义的时候为数值类型的变量(数据格式为整数)，用来给故障和补偿设备安装点编号。本步骤中，“Fault_faultnum_Comp_compnum”为字符串格式，在命名的时候，需要将整数型的数值变量(faultnum和compnum)转换成字符串，于是用到num2str函数将数值变量faultnum和compnum转换成字符串。num2str函数是matlab函数库里边现成的函数，直接调用即可。
按照上述规则统一命名有两个好处：一是方便之后调用稳定程序时可根据当前循环中faultnum和compnum的值传递输入数据文件(.SWI文件)；二是由于稳定程序输出结果文件(.SWX文件)的名称与输入数据文件(.SWI文件)相同，所以统一命名能方便读取稳定计算结果。
第S3.4步，在MATLAB里利用dos命令调用BPA稳定程序(swnt.exe)，并将第S3.3步中生成的稳定程序输入数据文件(.SWI文件)作为传递参数导入稳定程序中，执行稳定计算，如图2所示。图中，path_swnt为字符串格式的BPA稳定程序的路径；path_workspace为字符串格式的当前MATLAB脚本所在文件夹路径；
filename为第S3.3步中所生成的.SWI格式稳定程序输入数据文件名。在所遍历的每种组合中执行稳定计算将得到相应的全网负荷节点的电压时域分布结果，这些结果以所遍历的组合为单位，记录在不同的.SWX稳定计算结果文件(稳定程序输出结果文件)中。由于第S3.3步的操作，所得计算结果文件将按“Fault_faultnum_Comp_compnum”的命名规则排列。
在步骤S3中，根据稳定程序输出结果文件得到每种组合下全网所有负荷节点的电压时域分布。在稳定程序输入数据文件(.SWI文件)填好后执行稳定计算，系统会生成稳定程序输出结果文件(.SWX文件)，该.SWX文件中具有当前条件(组合)下全网所有负荷节点的电压时域分布。
第S4步.利用MATLAB读取第S3.4步中的BPA稳定程序输出结果文件(.SWX文件)，得到当节点j(j＝1,2,...,n)安装动态无功补偿设备后，各负荷节点在不同故障情况下的电压时域分布。并根据权重系数加权求和，得到中间变量R1j。R1j的计算式如下式所示：
第S5步.不考虑安装动态无功补偿设备，以批处理的方式遍历nf种不同故障的情况，得到每种情况下全网所有负荷节点的电压时域分布。
具体地：
第S5.1步，根据第S3步中定义的故障编号变量“faultnum”，利用MATLAB的for循环功能遍历不同故障情况。
第S5.2步，在每一个循环中，利用MATLAB读取原始的.SWI格式稳定程序输入数据文件，删除原始文件中第startline_F+faultnum-1行首列的注释符号，使相应故障卡生效。
第S5.3步，将第S5.2步中修改后的.SWI格式文件统一按照“No_Comp_Fault_faultnum”的规则重命名，生成每种故障情况中稳定计算需用到的.SWI格式稳定程序输入数据文件。由于文件名称为字符串格式，而faultnum的格式为整数，故需要通过num2str函数将数值变量faultnum转换成的字符串，再与固定字符串No_Comp_Fault_组合成文件名。
第S5.4步，在MATLAB里利用dos命令调用BPA稳定程序(swnt.exe)，并将第S5.3步生成的稳定程序输入数据文件(.SWI文件)作为传递参数导入稳定程序中，执行稳定计算。调用过程除了filename为第S5.3步中所生成的.SWI格式稳定程序输入数据文件名，其余与第S3.4步相似。在所遍历的每种故障情况中执行稳定计算将得到相应的全网负荷节点的电压时域分布结果，这些结果以所遍历的故障情况为单位，记录在不同的.SWX稳定计算结果文件(稳定程序输出结果文件)中。由于第S5.3步的操作，所得计算结果文件将按“No_Comp_Fault_faultnum”的命名规则排列。
第S6步.利用MATLAB读取第S5.4步中的BPA稳定程序输出结果文件(.SWX文件)，得到不安装动态无功补偿设备时，各负荷节点在不同故障情况下的电压时域分布。并根据权重系数加权求和，得到中间变量R2。R2的计算式如下式所示：
第S7步.由下式计算得到负荷节点j上的轨迹灵敏度TSIj：
TSIj＝R1j-R2。
由于大型受端电网中不止一个负荷节点，所以最终可得到一个n维数组TSI[n]，用来记录大型受端电网中所有负荷节点上的轨迹灵敏度。
需要说明的是：上述对大型受端电网轨迹灵敏度批处理计算方法的描述，不受执行步骤的限制，根据实际应用情况，步骤顺序可以调整。
本实施例提供的大型受端电网轨迹灵敏度批处理计算方法，利用MATLAB脚本遍历每种故障与负荷节点的组合，自动修改BPA稳定程序所用到的故障参数卡片与动态无功补偿设备卡片，从而进行稳定计算；并自动对输出数据文件按规律进行命名管理，实现对每种组合电压变化曲线的自动读取，提高计算效率，最终得到受端电网每一个负荷节点上的轨迹灵敏度。现有的手动方法计算大型受端电网所有负荷节点的轨迹灵敏度指标涉及到很大的工程量，而采用本实施例的好处在于提供的方法，所有过程均可由MATLAB脚本实现，自动调用BPA稳定程序进行批处理计算，避免重复而枯燥的手动操作，不仅大大增加了计算效率，还避免了出错。
以上对本发明的具体实施例进行了描述。需要理解的是，本发明并不局限于上述特定实施方式，本领域技术人员可以在权利要求的范围内做出各种变形或修改，这并不影响本发明的实质内容。
Claims (10)
1.一种大型受端电网轨迹灵敏度批处理计算方法，其特征在于，基于BPA和MATLAB实现，包括如下步骤：
S1.利用MATLAB根据BPA潮流文件所提供的受端电网的网络数据分别生成故障集文件、动态无功补偿设备安装点集文件和负荷节点电压输出语句集文件，并记录下故障个数nf和动态无功补偿设备安装点个数nc；
S2.根据S1生成的故障集文件、动态无功补偿设备安装点集文件和负荷节点电压输出语句集文件，形成原始的.SWI格式稳定程序输入数据文件，并记录下故障参数卡和动态无功补偿设备参数卡的起始行信息；
S3.以批处理的方式遍历nf×nc种不同故障搭配不同动态无功补偿设备安装点的组合，得到每种组合下全网所有负荷节点的电压时域分布，所述每种组合下全网所有负荷节点的电压时域分布记录在.SWX格式的稳定程序输出结果文件中；
S4.利用MATLAB读取S3中的稳定程序输出结果文件，得到当节点j(j＝1,2,...,n)安装动态无功补偿设备后，各负荷节点在不同故障情况下的电压时域分布；并根据权重系数加权求和，得到中间变量R1j；
S5.不考虑安装动态无功补偿设备，以批处理的方式遍历nf种不同故障的情况，得到每种情况下全网所有负荷节点的电压时域分布，并记录在.SWX格式的稳定程序输出结果文件中；
S6.利用MATLAB读取S5中的稳定程序输出结果文件，得到不安装动态无功补偿设备时，各负荷节点在不同故障情况下的电压时域分布；并根据权重系数加权求和，得到中间变量R2；
S7.由下式计算得到负荷节点j上的轨迹灵敏度TSIj：
TSIj＝R1j-R2。
2.根据权利要求1所述的大型受端电网轨迹灵敏度批处理计算方法，其特征在于，所述故障集文件、动态无功补偿设备安装点集文件和负荷节点电压输出语句集文件均为txt格式。
3.根据权利要求1或2所述的大型受端电网轨迹灵敏度批处理计算方法，其特征在于，所述故障集文件包括所有典型故障对应的BPA稳定程序卡片；所述动态无功补偿设备安装点集文件包括在所要计算轨迹灵敏度指标的所有节点上安装动态无功补偿设备的BPA稳定程序卡片；所述负荷节点电压输出语句集文件包括全网所有负荷节点上电压时域分布信息的输出指令语句。
4.根据权利要求1所述的大型受端电网轨迹灵敏度批处理计算方法，其特征在于，所述原始的.SWI格式稳定程序输入数据文件中，包括所要遍历所有组合间相同的部分以及不同组合间产生变化的部分。
5.根据权利要求4所述的大型受端电网轨迹灵敏度批处理计算方法，其特征在于，所述所要遍历所有组合间相同的部分包括：控制数据卡、发电机模型卡、励磁系统模型卡、电力系统稳定器模型卡、调速器和原动机模型卡、直流控制系统模型卡、负荷模型卡以及输出控制卡。
6.根据权利要求4所述的大型受端电网轨迹灵敏度批处理计算方法，其特征在于，所述不同组合间产生变化的部分包括：故障操作模型卡和动态无功补偿设备模型卡。
7.根据权利要求6所述的大型受端电网轨迹灵敏度批处理计算方法，其特征在于，所有故障操作模型卡和动态无功补偿设备模型卡均需要打上注释符号，以便求解到某种组合时能以去掉注释符号的方式形成相应的稳定程序输入数据文件。
8.根据权利要求7所述的大型受端电网轨迹灵敏度批处理计算方法，其特征在于，所述S3包括如下子步骤：
S3.1，定义故障编号变量faultnum和动态无功补偿设备安装点编号变量compnum，利用MATLAB的for循环功能遍历不同故障与安装点的组合；其中，变量faultnum和compnum的格式均为整数；
S3.2，利用MATLAB读取原始的.SWI格式稳定程序输入数据文件，根据故障编号变量faultnum、动态无功补偿设备安装点编号变量compnum以及之前记录下的故障参数卡、补偿设备参数卡起始行信息，删除相应故障操作模型卡和动态无功补偿设备模型卡首列的注释符号，使故障操作模型卡和动态无功补偿设备模型卡生效；
S3.3，利用MATLAB生成每种组合中稳定计算需用到的.SWI格式稳定程序输入数据文件，并统一按照Fault_faultnum_Comp_compnum的规则命名；文件名Fault_faultnum_Comp_compnum的格式为字符串，其中，变量faultnum和compnum的格式由整数通过num2str函数转换成字符串，再与固定字符串Fault_和_Comp_组合；
S3.4，在MATLAB里利用dos命令调用BPA稳定程序，并将S3.3中生成的稳定程序输入数据文件作为传递参数导入BPA稳定程序中，执行稳定计算；在所遍历的每种组合中执行稳定计算将得到相应的全网所有负荷节点的电压时域分布结果，这些结果以所遍历的组合为单位，记录在不同的稳定程序输出结果文件中；所得稳定程序输出结果文件将按Fault_faultnum_Comp_compnum的命名规则排列。
9.根据权利要求8所述的大型受端电网轨迹灵敏度批处理计算方法，其特征在于，所述S5包括如下子步骤：
S5.1，根据S3.1中定义的故障编号变量faultnum，利用MATLAB的for循环功能遍历不同故障情况；
S5.2，利用MATLAB读取原始的.SWI格式稳定程序输入数据文件，根据faultnum以及之前记录下的故障参数卡起始行信息，删除相应故障操作模型卡首列的注释符号，使故障操作模型卡生效；
S5.3，利用MATLAB生成每种故障情况中稳定计算需用到的稳定程序输入数据文件，并统一按照No_Comp_Fault_faultnum的规则命名；文件名No_Comp_Fault_faultnum的格式为字符串，其中，变量faultnum的格式由整数通过num2str函数转换成字符串，再与固定字符串No_Comp_Fault_组合；
S5.4，在MATLAB里利用dos命令调用BPA稳定程序，并将S5.3生成的稳定程序输入数据文件作为传递参数导入BPA稳定程序中，执行稳定计算；在所遍历的每种故障情况中执行稳定计算将得到相应的全网所有负荷节点的电压时域分布结果，这些结果以所遍历的故障情况为单位，记录在不同的稳定程序输出结果文件中；所得稳定程序输出结果文件将按No_Comp_Fault_faultnum的命名规则排列。
10.根据权利要求1所述的大型受端电网轨迹灵敏度批处理计算方法，其特征在于，
R1j的计算式如下式所示：
R2的计算式如下式所示：
式中，Nl为考察故障总数；Nk为总的时间采样点数；WFl为故障Fl的权重；Wdi为母线i的权重；Wk为tk时刻的权重；△Qj(COMP)为在节点j安装动态无功补偿设备后，该设备释放的无功功率；Vi,l(tk,Qj0+△Qj(COMP))为节点j安装额定容量为SCOMP的动态无功补偿设备后，在故障Fl下，节点电压Vi在t＝tk时刻的值；Vi,l(tk,Qj0)为未安装动态无功补偿设备情况下，在故障Fl下，节点电压Vl在t＝tk时刻的值。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810273074.0A CN108631297B (zh) | 2018-03-29 | 2018-03-29 | 大型受端电网轨迹灵敏度批处理计算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810273074.0A CN108631297B (zh) | 2018-03-29 | 2018-03-29 | 大型受端电网轨迹灵敏度批处理计算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108631297A true CN108631297A (zh) | 2018-10-09 |
| CN108631297B CN108631297B (zh) | 2021-05-04 |
Family
ID=63696432
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810273074.0A Active CN108631297B (zh) | 2018-03-29 | 2018-03-29 | 大型受端电网轨迹灵敏度批处理计算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108631297B (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110311384A (zh) * | 2019-07-19 | 2019-10-08 | 南方电网科学研究院有限责任公司 | 电力系统动态无功补偿方法、装置、设备及存储介质 |
| CN111507474A (zh) * | 2020-06-18 | 2020-08-07 | 四川大学 | 一种动态调整Batch-size的神经网络分布式训练方法 |
| CN112600196A (zh) * | 2020-12-02 | 2021-04-02 | 国网安徽省电力有限公司 | 基于bpa节点注入功率的支路功率快速灵敏度计算方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103457279A (zh) * | 2013-09-23 | 2013-12-18 | 广东电网公司电力调度控制中心 | 大规模电网的分区d-statcom综合优化配置方案的确定方法 |
| CN103701140A (zh) * | 2014-01-06 | 2014-04-02 | 国家电网公司 | 提高交直流电网暂态电压稳定的动态无功备用优化方法 |
-
2018
- 2018-03-29 CN CN201810273074.0A patent/CN108631297B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103457279A (zh) * | 2013-09-23 | 2013-12-18 | 广东电网公司电力调度控制中心 | 大规模电网的分区d-statcom综合优化配置方案的确定方法 |
| CN103701140A (zh) * | 2014-01-06 | 2014-04-02 | 国家电网公司 | 提高交直流电网暂态电压稳定的动态无功备用优化方法 |
Non-Patent Citations (1)
| Title |
|---|
| 范丁中: "受端电网动态无功装置优化配置方法研究", 《中国优秀硕士学位论文全文数据库》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110311384A (zh) * | 2019-07-19 | 2019-10-08 | 南方电网科学研究院有限责任公司 | 电力系统动态无功补偿方法、装置、设备及存储介质 |
| CN111507474A (zh) * | 2020-06-18 | 2020-08-07 | 四川大学 | 一种动态调整Batch-size的神经网络分布式训练方法 |
| CN111507474B (zh) * | 2020-06-18 | 2022-07-01 | 四川大学 | 一种动态调整Batch-size的神经网络分布式训练方法 |
| CN112600196A (zh) * | 2020-12-02 | 2021-04-02 | 国网安徽省电力有限公司 | 基于bpa节点注入功率的支路功率快速灵敏度计算方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108631297B (zh) | 2021-05-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108631297A (zh) | 大型受端电网轨迹灵敏度批处理计算方法 | |
| Xu et al. | A decomposition-based practical approach to transient stability-constrained unit commitment | |
| CN105224511A (zh) | 一种智能变电站scd文件的校验方法 | |
| CN107861868A (zh) | 一种提取自动化测试对象的方法及系统 | |
| CN104391692A (zh) | 用于读写psd-bpa数据文件的面向对象编程方法 | |
| CN106777508A (zh) | 一种psd‑bpa的电力系统仿真辅助应用系统 | |
| CN105186502A (zh) | 基于安全域的含双馈风机电力系统暂态稳定性分析方法 | |
| CN105426350A (zh) | Scd文件中提取电压等级、间隔及ied信息的方法 | |
| CN103186366A (zh) | 基于cuda并行计算实现电力系统电磁暂态实时仿真测试方法 | |
| Duckwitz | Power system inertia | |
| CN105373008A (zh) | 一种电力系统用全物理模型仿真方法 | |
| CN111130135A (zh) | 一种适用于高比例新能源接入下的电力系统惯量计算方法 | |
| CN113987781A (zh) | 面向电力系统暂态仿真的设备符号化建模及其应用方法 | |
| CN103942322A (zh) | 一种基于xml文件的配网模型差异数据合并方法 | |
| CN104217003A (zh) | 一种数据建模系统 | |
| CN103886419A (zh) | 一种电网运行方式规划校核分析系统 | |
| CN105279297B (zh) | 一种基于工程应用的pss/e电力系统等值方法 | |
| CN112072692A (zh) | 一种新能源发电场站的阻抗等值方法及装置 | |
| CN104463699A (zh) | 一种智能变电站站控层信息的提取方法及装置 | |
| CN112865071B (zh) | 直流受端电网分布式光伏高渗透下频率偏移预测方法 | |
| CN111553040B (zh) | 一种基于gpu加速的电网拓扑分析高性能计算方法及装置 | |
| Xie et al. | A static equivalent method considering dfig grid-connected control strategy | |
| CN106407534A (zh) | 一种电磁暂态仿真数据的处理方法及装置 | |
| CN112198452A (zh) | 一种适用于工程实用化的新能源短路电流表达式的构建方法 | |
| CN111695332A (zh) | 一种基于变电站监控信息表自动生成cime文件的方法 |
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 |