利用python调用PSS/E进行电力系统仿真

PSS/E是西门子公司（PTI）开发的电力系统仿真软件，可以实现潮流计算、短路计算、机电暂态等功能，但是在科研或者项目研究时可能需要将仿真数据导出，在MATLAB、python等平台上进行进一步处理，进行更进一步的研究，换句话说，仿真计算只是项目研究的一个部分，仿真计算仅仅是数据来源。对于**“如何将仿真数据导出”**这一问题，PSS/E开发了psspy模块，实现了PSS/E与python2.7的交互，可以通过在python编辑器上调用psspy等模块实现对PSS/E的控制，进行仿真计算。

写在前面

psspy模块是32位的，因此要调用psspy必须使用python2.7-bit32位，64位的python是无法实现的。如果不确定python版本，可以通过cmd输入python查看：

为了简便起见，小编安装了anaconda-2.1.0-x86安装包

环境配置

首先将psspy导入，psspy的地址一般在PSSE安装地址中的PSSEBIN文件夹中。然后配置环境变量，所谓环境变量就是python可执行文件所在的位置。配置时把psspy模块所在文件夹添加到sys.environ[‘PATH’]中即可。最后导入psspy包

import os, sys
PYTHONPATH = r'C:\Program Files (x86)\PTI\PSSE33\PSSBIN'
sys.path.append(PYTHONPATH)
os.environ['PATH'] += ';' + PYTHONPATH
import psspy

调用PSS/E

导入psspy后就可以通过调用psspy模块进行潮流计算等仿真工作了，像在PSS/E图像界面操作一样，首先打开PSS/E进行初始化，打开文件，进行潮流计算……对应的代码如下：

if __name__ == '__main__':CASE = r"C:\Program Files (x86)\PTI\PSSE33\EXAMPLE\savnw.sav"psspy.psseinit(buses=2000) psspy.case(CASE) b=psspy.fnsl()

以PSS/E自带的savnw.sav为例，CASE保存了该算例的绝对地址。
因为导入的是整个模块，所以调用模块里的函数时要添加模块名。
psseinit（buses）函数用于初始化PSS/E，buses参数是PSS/E母线数目的最大值；
case(CASE)函数用于打开算例文件，参数是算例所在路径；
fnsl()为牛拉法潮流计算函数，共有8个参数，具体可以参考PSS/E的API手册。不给参数是按照缺省值运行。

运行结果：

如图，自上而下一次是打开PSS/E，读取文件，潮流计算结果。经过6次迭代，计算结果收敛，

全部代码

# Set environment
import os, sys
PYTHONPATH = r'C:\Program Files (x86)\PTI\PSSE33\PSSBIN'#
sys.path.append(PYTHONPATH)#告诉python模块导入的位置
os.environ['PATH'] += ';' + PYTHONPATH #增加环境变量 import psspyif __name__ == '__main__':CASE = r"C:\Program Files (x86)\PTI\PSSE33\EXAMPLE\savnw.sav"psspy.psseinit(buses=2000) #初始化PSSE，buses是母线数量psspy.case(CASE) #打开文件b=psspy.fnsl()

参考资料

Run PSSE from Python and not the other way around - Python for Power Systems