[发明专利]一种船舶电能扰动时间定位与类型识别的弧长差分序列方法

摘要

本发明提出一种船舶电能质量扰动时间定位与类型识别的弧长差分序列方法。该方法基于弧长差分序列对形态滤波器滤波后的电能信号进行处理，可以实现船舶电力系统中常见的5种单一扰动和3种复合扰动的时间定位和类型识别。该方法首先采用一种改进的交替混合形态滤波器，对船舶电能信号进行滤波处理，然后提出一种弧长差分序列的定义，并用于精确时间定位扰动时间，提取出扰动时频域特征‑幅值变化特征和波峰特征，最后将特征量和扰动时间输入到定义了扰动类型识别规则的类型识别器，实现扰动类型识别。本发明的方法可用于快速类型识别实际船舶电力系统中的扰动，适用于实时检测和类型识别。

主权项

一种船舶电能质量扰动时间定位与类型识别的弧长差分序列方法，定义弧长差分序列如下：设信号中某采样点为a点，a点之后的每四分之一信号周期，或二分之一信号周期，或一个信号周期的点依次记为b点、c点、d点等，并称a点为弧长序列的起始点，b点、c点、d点等为弧长序列按四分之一，或二分之一，或一个电周期选取的周期点，这些点两两间曲线的弧长可组成弧长序列为则弧长差分序列可以定义如下式(1)，将Δs中的非零项依次称为第一非零元素、第二非零元素、第三非零元素等：设按照每个电周期(0.02s)采样k个点对电能信号进行采样，则四分之一电周期(0.005s)的采样点数为X0＝{x1,x2,…,xn}为采样信号，且x1为过零点；其特征在于船舶电能质量扰动时间定位与类型识别的弧长差分序列方法包括以下步骤：步骤1)取信号的幅值A为基准值，对原始采样信号取标准化。按下式(2)对采样信号X0进行标准化，得到标准化信号X1：步骤2)采用改进的交替混合形态滤波器对步骤1)中处理后的信号进行滤波，改进的交替混合形态滤波器如下式(3)所示：almix(n)＝((f)OC(g)+(f)CO(g))(n)/2   (3)其中：f为采样信号，g(m)为结构元素，Df1＝{p,…,q}，p＝[N/1000]，q＝N‑[N/1000]；Dg＝{0,1,2,…,M}，M为结构元素g(m)中元素的个数；N和M均为整数，N≥M；步骤3)对扰动时间进行精确定位，分为三步：a.按二分之一电周期计算步骤2)中的滤波信号的弧长差分序列，定位扰动时间区间；b.按四分之一电周期计算步骤2)中的滤波信号的弧长差分序列，定位扰动区间中的精确时刻；c.综合扰动时间区间和扰动区间中的精确时刻，完成扰动时间精确定位，具体过程如下：步骤3.1)按二分之一电周期计算步骤2)中的滤波信号的弧长差分序列，定位扰动时间区间：步骤3.1.1)首先，以过零点x1为起始点，按四分之一电周期取周期点。计算两周期点间的弧长，得到弧长序列为X3，最后再如式(1)所示计算得到弧长差分序列X4；步骤3.1.2)提取步骤弧长差分序列中的第一非零元素，查找弧长序列中与第一非零元素相关的弧，此弧所在的时间区间[t1,t1+0.005]就包含扰动发生时刻；步骤3.2)按四分之一电周期计算步骤2)中的滤波信号的弧长差分序列，定位扰动区间中的精确时刻：步骤3.2.1)对步骤2)中滤波后信号的第一个四分之一周期内的点按照用户选定的采样频率进行二次采样；步骤3.2.2)将第一个四分之一电周期内的采样点依次作为起始点，按二分之一电周期取周期点，计算得到弧长差分序列组；步骤3.2.3)提取步骤3.2.2)中的弧长差分序列组中的第一非零元素，组成第一非零元素集合；步骤3.2.4)对步骤3.2.3)中的第一非零元素集合中的元素取绝对值，并记为X5，查找这些元素中的最大值，记此最大值在集合中的序号为j；步骤3.3)综合扰动时间区间和扰动区间中的精确时刻，完成扰动时间精确定位，按下式(10)计算可得扰动时刻T，单位为秒：步骤4)根据步骤3)中的扰动时间定位结果和步骤2)中的滤波后信号，提取扰动发生阶段的弧长差分序列X4'和扰动发生阶段的滤波后信号X2'，进行扰动类型识别，分为两步：a.对X4'进行特征处理，可以识别暂升、暂降、暂态振荡和诊断扰动；b.对X2'进行特征处理，可以识别中断、暂态振荡、暂升/暂降和谐波，记滤波后信号X2'中有Z个采样点，具体的扰动类型识别步骤如下；步骤4.1)对扰动发生阶段的弧长差分序列X4'进行特征处理：步骤4.1.1)对扰动发生阶段的弧长差分序列X4'进行中值滤波，滤波的模板长度为且取为偶数，得到扰动发生后波形幅值特征量S；步骤4.1.2)判断波形幅值特征量S是否大于0，若是，则转步骤4.3)；若否，则转步骤4.4)；步骤4.1.3)判断S的第一个元素与最后一个元素是否大于0.1，即Sstart‑Send＞0.1，若是，则判定扰动为暂态振荡扰动；若否，则判定为暂升扰动；步骤4.1.4)判断S是否大于‑0.9且小于‑0.1，若是，则判定为暂降扰动；若否，则判定为中断；步骤4.2)取对扰动发生阶段的滤波后信号X2'进行特征处理：步骤4.2.1)对扰动发生阶段的滤波后信号X2'进行波峰个数检测，结果记为波动频率特征P；步骤4.2.2)判断P是否等于零，若是，则判定为中断；否则转步骤4.7)；步骤4.2.3)判断M≤P≤M+1是否成立，若是则判定为暂升或者暂降扰动，否则转步骤4.8)；步骤4.2.4)判断P≥N是否成立，若是则判定为暂态振荡，否则判断为谐波扰动，最高次谐波为[P/[z/k]]次谐波。