基于时空序列的大电网静态稳定边界特征提取方法及系统技术方案

本发明专利技术提供一种基于时空序列的大电网静态稳定边界特征提取方法及系统，属于大电网安全稳定运行领域。本发明专利技术方法包括：采集静态稳定边界全息样本数据；选取静态稳定边界样本数据；映射静态稳定边界样本数据为三张样本图片；利用基于尺度不变特性变换法提取源样本图片、网样本图片和荷样本图片三张样本图片特征，计算三张样本图片高斯差分尺度空间；对边界荷样本图片和边界网样本图片进行特征匹配；提取样本荷‑网的静态稳定边界特征。本发明专利技术避免了实际运行数据少且工况覆盖面小的缺点，得出的静态稳定边界特征更具说服力。

【技术实现步骤摘要】
基于时空序列的大电网静态稳定边界特征提取方法及系统
本专利技术涉及大电网安全稳定运行领域，具体地涉及一种基于时空序列的大电网静态稳定边界特征提取方法及系统。
技术介绍
随着经济的发展，用电负荷不断增长，电网不断逼近其输电能力极限，电力系统的安全稳定运行受到不断挑战；新能源、电力电子设备大规模接入电网，其不确定性严重影响了电网安全稳定运行能力，全球范围内发生多起大规模停电事故。这些都对电力系统边界特征提取提出了迫切的技术要求。传统静态稳定评估研究或稳定边界研究大多仅给出了静态稳定裕度或边界范围，没有直接根据稳定边界作出研究，如基于电力系统静态电压稳定域(Staticvoltagestabilityregion,SVSR)边界拓扑特性，实现高精度SVSR的快速检索；通过改进传统的直接延拓法，提出了一种二维参数稳定边界追踪方法。其研究方法大多针对物理模型和高强度仿真计算，传统模型难以适应当前电力系统，导致得出的结论与电力系统实际运行状态存在差异，其在线工程适用性无法保证。静态稳定边界的研究大多根据电压特性展开研究，分析的全面性无法保证。目前对于静态稳定的判定存在多种指标，但大多数指标只能体现电网的某一个侧面(如电压稳定指标仅反映电网的电压水平，功角稳定指标仅反映发电机的运行状态)，其无法从宏观上审视电网的整体状态。由于广域测量系统(wideareameasurementsystem，WAMS)在电力系统中广泛应用，数据驱动、信息驱动的电力系统分析理论框架被广泛提出，通过大数据分析技术实现大电网全景安全防御和稳定态势量化评估。电网存在大量仿真或实测数据，利用成熟的特征提取算法为电网数据特征的提取提供可能。有关特征提取的研究大多应用于图像处理领域，识别图像特征，可根据特征匹配算法实现图像特征的提取、匹配以及精准定位。
技术实现思路
本专利技术的目的在于直接挖掘大电网静态稳定临界状态下全息信息的关联程度，全面体现静稳极限时的运行状态，改进大电网静态稳定态势评估指标，提出了一种时空序列的大电网静态稳定边界特征提取方法及系统，所述方法包括：设定电网运行初始状态，利用连续潮流法，穷举多种负荷增长和发电机出力增加的运行工况，采集静态稳定边界全息样本数据；将全息样本数据按采样时间顺序排列，形成矩阵，公式如下：XN*T＝[x(t1),x(t2),.....,x(ti),...]其中，x(ti)为在采样时刻ti的全息样本数据，每个全息样本的量测数据构成的列向量，在全息样本数据中添加高斯白噪声，并对含噪音的全息样本数据进行平滑处理，使用实时分离窗技术，设置实时分离窗宽度TW，在ti时刻获取样本矩阵为：XN*Tw(ti)＝[x(ti-Tw+1),x(ti-Tw+2),.....,x(ti)]取数据处理后的X矩阵的最后一列作为静态稳定边界样本数据；将静态稳定边界样本数据按源、网和荷分类，将边界源样本、边界网样本和边界荷样本合成三个矩阵，具体方法为：以所有边界样本为行，电网源数据变量为列构成源矩阵，以此类推构成网矩阵和荷矩阵，将三个矩阵中的数据进行归一化处理，形成像素矩阵，映射为源样本图片、网样本图片和荷样本图片三张样本图片；利用基于尺度不变特性变换法提取源样本图片、网样本图片和荷样本图片三张样本图片特征；计算三张样本图片高斯差分尺度空间D(x,y,σ)，公式如下：D(x,y,σ)＝(G(x,y,kσ)-G(x,y,σ))*I(x,y)其中(x，y)是空间坐标，σ是尺度坐标，I(x,y)为图片像素矩阵，G(x,y,σ)是尺度坐标为σ下的尺度可变高斯函数，G(x,y,kσ)是尺度坐标为kσ下的尺度可变高斯函数，G(x,y,kσ)表达式如下：在高斯差分尺度空间上，针对三张样本图片的每一个采样点和图像域与尺度域的相邻点相比较，得到尺度空间上的极值点，并依据下式去除边缘响应：其中H为由高斯尺度空间构成的海森矩阵，Tr(H)为海森矩阵的迹，Det(H)为海森矩阵的行列式，γ为该矩阵最大特征值与最小特征值间的倍数；将高斯差分尺度空间上的关键点周围的采样点带入梯度公式，确定梯度主方向，以关键点为中心取16个8×8的窗口，计算梯度方向直方图，绘制每个梯度方向的累加值，形成一个128维的描述子；对边界荷样本图片和边界网样本图片进行特征匹配；计算边界荷样本图片和边界网样本图片特征描述子的欧式距离，若计算结果小于设定的阈值，则判定匹配成功；若大于阈值，则判定不成功，对下一个特征点描述子进行匹配。欧式距离D的计算方法如下式所示：D＝[arccos(des1*des2)]^1/2其中，des1为荷样本图片特征描述子,des2为网样本图片特征描述子；将荷样本图片与网样本图片匹配到的特征数与网样本图片的总特征数相除得到静稳极限状态下的样本荷-网的静态稳定边界特征Khw，Khw公式如下：其中，Fm1为荷样本图片与网样本图片匹配到的特征数，Fw为网样本图片的总特征数。可选的，全息样本数据包括发电机节点电压、流出有功功率、无功功率、母线电压、流出有功功率、无功功率、负荷有功功率和无功功率。可选的，电网运行初始状态包括网络拓扑结构和初始负荷状态。可选的，发电机出力增加的运行工况包括单负荷增长、多负荷增长、全负荷增长、单台发电机出力增长和多台发电机出力增长。可选的，ti取值范围从电网运行初始时刻到电网静态稳定极限状态运行时刻。可选的，静态稳定边界特征Khw根据不同的功率增长百分比行成静态稳定边界特征Khw的实时变化曲线。可选的，三个矩阵的合成具体流程把包括：以所有边界样本为行，电网源数据变量为列合成源矩阵、网矩阵和荷矩阵。本专利技术系统包括：采集模块：设定电网运行初始状态，利用连续潮流法，穷举多种负荷增长和发电机出力增加的运行工况，采集静态稳定边界全息样本数据；数据处理模块：将全息样本数据按采样时间顺序排列，形成矩阵，公式如下：XN*T＝[x(t1),x(t2),.....,x(ti),...]其中，x(ti)为在采样时刻ti的全息样本数据，每个全息样本的量测数据构成的列向量，在全息样本数据中添加高斯白噪声，并对含噪音的全息样本数据进行平滑处理，使用实时分离窗技术，设置实时分离窗宽度TW，在ti时刻获取样本矩阵为：XN*Tw(ti)＝[x(ti-Tw+1),x(ti-Tw+2),.....,x(ti)]取数据处理后的X矩阵的最后一列作为静态稳定边界样本数据；映射模块：将静态稳定边界样本数据按源、网和荷分类，将边界源样本、边界网样本和边界荷样本合成三个矩阵，将三个矩阵中的数据进行归一化处理，形成像素矩阵，映射为源样本图片、网样本图片和荷样本图片三张样本图片；计算模块：利用基于尺度不变特性变换法提取源样本图片、网样本图片和荷样本图片三张样本图片特征，计算三张样本图片高斯差分尺度空间D(x,y,σ)，公式如下：D(x,y,σ)＝(G(x,y,kσ)-G(x,y,σ))*I(x,y)其中(x，y)是空间坐标，σ是尺度坐标，I(x,y)为图片像素矩阵，G(x,y,σ)是尺度坐标为σ下的尺度可变高斯函数，G(x,y,kσ)是尺度坐标为kσ下的尺度可变高斯函数，G(x,y,kσ)表达式如下：在高斯差分尺度空间上，针对三张样本图片的每一个采样点和图像域与尺度域的相邻点相比较，得到尺度空间上的本文档来自技高网...

【技术保护点】
1.一种基于时空序列的大电网静态稳定边界特征提取方法，其特征在于，包括：设定电网运行初始状态，利用连续潮流法，穷举多种负荷增长和发电机出力增加的运行工况，采集静态稳定边界全息样本数据；将全息样本数据按采样时间顺序排列形成矩阵公式如下：XN*T＝[x(t1),x(t2),.....,x(ti),...]    (1)其中，x(ti)为在采样时刻ti的全息样本数据，每个全息样本的量测数据构成的列向量，在全息样本数据中添加高斯白噪声，并对含噪音的全息样本数据进行平滑处理，使用实时分离窗技术，设置实时分离窗宽度TW，在ti时刻获取样本矩阵为：XN*Tw(ti)＝[x(ti‑Tw+1),x(ti‑Tw+2),.....,x(ti)]    (2)取数据处理后的X矩阵的最后一列作为静态稳定边界样本数据；将静态稳定边界样本数据按源、网和荷分类，将边界源样本、边界网样本和边界荷样本合成三个矩阵，将三个矩阵中的数据进行归一化处理，形成像素矩阵，映射为源样本图片、网样本图片和荷样本图片三张样本图片；利用基于尺度不变特性变换法提取源样本图片、网样本图片和荷样本图片三张样本图片特征，计算三张样本图片高斯差分尺度空间D(x,y,σ)，公式如下：D(x,y,σ)＝(G(x,y,kσ)‑G(x,y,σ))*I(x,y)    (3)其中(x，y)是空间坐标，σ是尺度坐标，I(x,y)为图片像素矩阵，G(x,y,σ)是尺度坐标为σ下的尺度可变高斯函数，G(x,y,kσ)是尺度坐标为kσ下的尺度可变高斯函数，G(x,y,kσ)表达式如下：...

【技术特征摘要】
1.一种基于时空序列的大电网静态稳定边界特征提取方法，其特征在于，包括：设定电网运行初始状态，利用连续潮流法，穷举多种负荷增长和发电机出力增加的运行工况，采集静态稳定边界全息样本数据；将全息样本数据按采样时间顺序排列形成矩阵公式如下：XN*T＝[x(t1),x(t2),.....,x(ti),...](1)其中，x(ti)为在采样时刻ti的全息样本数据，每个全息样本的量测数据构成的列向量，在全息样本数据中添加高斯白噪声，并对含噪音的全息样本数据进行平滑处理，使用实时分离窗技术，设置实时分离窗宽度TW，在ti时刻获取样本矩阵为：XN*Tw(ti)＝[x(ti-Tw+1),x(ti-Tw+2),.....,x(ti)](2)取数据处理后的X矩阵的最后一列作为静态稳定边界样本数据；将静态稳定边界样本数据按源、网和荷分类，将边界源样本、边界网样本和边界荷样本合成三个矩阵，将三个矩阵中的数据进行归一化处理，形成像素矩阵，映射为源样本图片、网样本图片和荷样本图片三张样本图片；利用基于尺度不变特性变换法提取源样本图片、网样本图片和荷样本图片三张样本图片特征，计算三张样本图片高斯差分尺度空间D(x,y,σ)，公式如下：D(x,y,σ)＝(G(x,y,kσ)-G(x,y,σ))*I(x,y)(3)其中(x，y)是空间坐标，σ是尺度坐标，I(x,y)为图片像素矩阵，G(x,y,σ)是尺度坐标为σ下的尺度可变高斯函数，G(x,y,kσ)是尺度坐标为kσ下的尺度可变高斯函数，G(x,y,kσ)表达式如下：在高斯差分尺度空间上，针对三张样本图片的每一个采样点和图像域与尺度域的相邻点相比较，得到尺度空间上的极值点去除边缘响应，去除边缘响应公式如下：其中H为由高斯尺度空间构成的海森矩阵，Tr(H)为海森矩阵的迹，Det(H)为海森矩阵的行列式，γ为该矩阵最大特征值与最小特征值间的倍数；将高斯差分尺度空间上的关键点周围的采样点带入梯度公式，确定梯度主方向，以关键点为中心取16个8×8的窗口，计算梯度方向直方图，绘制每个梯度方向的累加值，形成一个128维的描述子；对边界荷样本图片和边界网样本图片进行特征匹配；计算边界荷样本图片和边界网样本图片特征描述子的欧式距离，若计算结果小于设定的阈值，则判定匹配成功；若大于阈值，则判定不成功，对下一个特征点描述子进行匹配，欧式距离D的计算方法如下式所示：D＝[arccos(des1*des2)]^1/2(6)其中，des1为荷样本图片特征描述子,des2为网样本图片特征描述子；将荷样本图片与网样本图片匹配到的特征数与网样本图片的总特征数相除得到静稳极限状态下的样本荷-网的静态稳定边界特征Khw，Khw公式如下：其中，Fm1为荷样本图片与网样本图片匹配到的特征数，Fw为网样本图片的总特征数。2.根据权利要求1所述的方法，所述的静态稳定边界特征Khw根据不同的功率增长百分比行成静态稳定边界特征Khw的实时曲线。3.根据权利要求1所述的方法，所述的全息样本数据包括发电机节点电压、流出有功功率、无功功率、母线电压、流出有功功率、无功功率、负荷有功功率和无功功率。4.根据权利要求1所述的方法，所述的电网运行初始状态包括网络拓扑结构和初始负荷状态。5.根据权利要求1所述的方法...

【专利技术属性】
技术研发人员：刘道伟，田一童，李柏青，王克非，李丹，高德宾，侯凯元，夏得明，杨红英，章锐，李宗翰，赵高尚，郑恒峰，李京，杨少波，
申请(专利权)人：中国电力科学研究院有限公司，国家电网有限公司，国家电网公司东北分部，
类型：发明
国别省市：北京,11