发明名称 基于受扰电压轨迹簇特征的电力系统暂态稳定性预测方法
摘要 本发明公开了一种基于WAMS受扰电压轨迹簇特征的电力系统暂态稳定性预测方法,该方法对电力系统中WAMS系统所采集到的历史信息进行初步筛选,得到故障后所有关键发电机节点电压幅值的时域轨迹簇构成原始数据集;对原始数据计算轨迹簇的27个几何特征;利用Relief算法对所有特征量进行权重计算,选出与系统暂态稳定性强相关的若干个广域故障特征,并作为暂态稳定预测算法的输入数据集;基于输入数据集构造SVM预测模型。本发明不建立电力系统的分析模型,而是根据WAMS系统得到的电力系统的响应信息,直接预测电力系统的暂态稳定性;本发明不但能够快速预测系统稳定性,具有极强的适应性和鲁棒性。
申请公布号 CN104578053B 申请公布日期 2017.01.04
申请号 CN201510013138.X 申请日期 2015.01.09
申请人 北京交通大学 发明人 吴俊勇;冀鲁豫;郝亮亮;周艳真;于之虹;严剑峰;华科
分类号 H02J3/00(2006.01)I 主分类号 H02J3/00(2006.01)I
代理机构 北京正理专利代理有限公司 11257 代理人 张文祎
主权项 基于受扰电压轨迹簇特征的电力系统暂态稳定性预测方法,其特征在于,所述方法的步骤包括:S1、基于WAMS系统对电力系统采集得到的历史数据信息进行初步筛选,得到故障后所有关键发电机节点电压幅值的时域轨迹簇,并构成原始数据集;所述步骤S1中,在某个故障情况下,得到故障清除后m个周波的所有n个发电机节点电压轨迹簇,记为轨迹簇矩阵:{x<sub>ij</sub>}<sub>m×n</sub>,其中,x<sub>ij</sub>即为某一周波下某一轨迹采样点的值,i表示采样点的周波序号,j表示轨迹簇中的某一条轨迹序号;S2、对所述原始数据进行轨迹簇中27个几何特征量的计算;所述步骤S2中所述27个几何特征量包括:特征1:瞬时轨迹质心,<img file="FDA0001054157060000011.GIF" wi="270" he="126" />c<sub>i</sub>为一m行的列向量;特征2:轨迹簇断面离散度,<img file="FDA0001054157060000012.GIF" wi="430" he="157" />特征3:上包络线,ue<sub>i</sub>={max(x<sub>i</sub>)},i=1,2,...,m;特征4:下包络线,le<sub>i</sub>={min(x<sub>i</sub>)},i=1,2,...,m;特征5:瞬时轨迹中心,me<sub>i</sub>={(max(x<sub>i</sub>)+min(x<sub>i</sub>))/2},i=1,2,...,m;特征6:上包络线与质心线之差,dcue<sub>i</sub>={|ue<sub>i</sub>‑c<sub>i</sub>|},i=1,2,...,m;特征7:下包络线与质心线之差,dcle<sub>i</sub>={|le<sub>i</sub>‑c<sub>i</sub>|},i=1,2,...,m;特征8:包络线高度,he<sub>i</sub>={|max(x<sub>i</sub>)‑min(x<sub>i</sub>)|},i=1,2,...,m;特征9:轨迹质心与中心之差,dcm<sub>i</sub>={|c<sub>i</sub>‑me<sub>i</sub>|},i=1,2,...,m;特征10:质心线瞬时变化率,<img file="FDA0001054157060000013.GIF" wi="702" he="135" />其中,h为采样时间间隔;特征11:离散度瞬时变化率,<img file="FDA0001054157060000014.GIF" wi="718" he="135" />特征12:上包络线变化率,<img file="FDA0001054157060000015.GIF" wi="774" he="135" />特征13:下包络线变化率,<img file="FDA0001054157060000016.GIF" wi="741" he="142" />特征14:中心线变化率,<img file="FDA0001054157060000021.GIF" wi="806" he="135" />特征15:包络线高度变化率,<img file="FDA0001054157060000022.GIF" wi="774" he="139" />特征16:轨迹曲率,<img file="FDA0001054157060000023.GIF" wi="1172" he="183" />特征17:质心线曲率,<img file="FDA0001054157060000024.GIF" wi="925" he="135" />特征18:离散度曲率,<img file="FDA0001054157060000025.GIF" wi="942" he="138" />特征19:上包络线曲率,<img file="FDA0001054157060000026.GIF" wi="1014" he="135" />特征20:下包络线曲率,<img file="FDA0001054157060000027.GIF" wi="973" he="134" />特征21:中心线曲率,<img file="FDA0001054157060000028.GIF" wi="1069" he="142" />特征22:质心线变化加速度,<img file="FDA0001054157060000029.GIF" wi="757" he="139" />特征23:离散度变化加速度,<img file="FDA00010541570600000210.GIF" wi="770" he="134" />特征24:上包络线加速度,<img file="FDA00010541570600000211.GIF" wi="830" he="135" />特征25:下包络线加速度,<img file="FDA00010541570600000212.GIF" wi="790" he="135" />特征26:中心线加速度,<img file="FDA00010541570600000213.GIF" wi="862" he="135" />特征27:包络线高度变化加速度,<img file="FDA00010541570600000214.GIF" wi="830" he="139" />其中,轨迹曲线曲率的计算采用了三点二次插值算法;对N个电压轨迹簇样本计算特征量后得到的数据集记为X<sub>N×M</sub>,N为样本个数,M为特征量个数,此处M=27;S3、利用Relief算法对步骤S2中所述的27个几何特征量进行权重w<sub>i</sub>计算,并选出与系统暂态稳定性强相关的若干个广域故障特征作为暂态稳定预测算法的输入数据集;S4、根据所述输入数据集构建支持向量机SVM预测模型,并采取交叉验证的方法得到最佳精度的预测模型;S5,利用步骤S2和S3对WAMS系统采集得到的新电力系统数据信息进行处理,构建预测数据集;S6、利用步骤S4中构建的预测模型对所述预测数据集进行判断,获得暂态稳定性的预测结果。
地址 100044 北京市海淀区上园村3号