[发明专利]大型受端电网轨迹灵敏度批处理计算方法

权利要求书

1.一种大型受端电网轨迹灵敏度批处理计算方法，其特征在于，基于BPA和MATLAB实现，包括如下步骤：

S1.利用MATLAB根据BPA潮流文件所提供的受端电网的网络数据分别生成故障集文件、动态无功补偿设备安装点集文件和负荷节点电压输出语句集文件，并记录下故障个数n_f和动态无功补偿设备安装点个数n_c；

S2.根据S1生成的故障集文件、动态无功补偿设备安装点集文件和负荷节点电压输出语句集文件，形成原始的.SWI格式稳定程序输入数据文件，并记录下故障参数卡和动态无功补偿设备参数卡的起始行信息；

S3.以批处理的方式遍历n_f×n_c种不同故障搭配不同动态无功补偿设备安装点的组合，得到每种组合下全网所有负荷节点的电压时域分布，所述每种组合下全网所有负荷节点的电压时域分布记录在.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_组合；

S3.4，在MATLAB里利用dos命令调用BPA稳定程序，并将S3.3中生成的稳定程序输入数据文件作为传递参数导入BPA稳定程序中，执行稳定计算；在所遍历的每种组合中执行稳定计算将得到相应的全网所有负荷节点的电压时域分布结果，这些结果以所遍历的组合为单位，记录在不同的稳定程序输出结果文件中；所得稳定程序输出结果文件将按Fault_faultnum_Comp_compnum的命名规则排列；

S4.利用MATLAB读取S3中的稳定程序输出结果文件，得到当节点j，且j＝1,2,...,n,安装动态无功补偿设备后，各负荷节点在不同故障情况下的电压时域分布；并根据权重系数加权求和，得到中间变量R1_j；

S5.不考虑安装动态无功补偿设备，以批处理的方式遍历n_f种不同故障的情况，得到每种情况下全网所有负荷节点的电压时域分布，并记录在.SWX格式的稳定程序输出结果文件中，具体包括如下子步骤：

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的命名规则排列；

S6.利用MATLAB读取S5中的稳定程序输出结果文件，得到不安装动态无功补偿设备时，各负荷节点在不同故障情况下的电压时域分布；并根据权重系数加权求和，得到中间变量R2；

S7.由下式计算得到负荷节点j上的轨迹灵敏度TSI_j：

TSI_j＝R1_j-R2。