发明名称 基于关联矩阵的配电网行波故障定位方法
摘要 本发明涉及电力系统自动化技术领域,是一种基于关联矩阵的配电网行波故障定位方法,包括第一步,数据预处理;第二步,构建关联矩阵;第三步,故障区间判断;第四步,判断向量元素满足相关性及对称性;第五步,坏数据识别;第六步,构建误差矩阵,并对故障进行定位,测距结果输出。本发明基于误差整体分布原则,通过构建误差矩阵,查找全部节点中整体误差最小的故障定位结果,提高了故障定位的整体精度。本发明利用多端数据构建关联矩阵,通过关联矩阵的对称性及相关性确定故障区间并识别异常装置数据,在此基础上进一步完成了故障点精确定位,从而提高了配电网故障定位的可靠性及精度。
申请公布号 CN106443360A 申请公布日期 2017.02.22
申请号 CN201610994603.7 申请日期 2016.11.11
申请人 国网新疆电力公司巴州供电公司;国家电网公司 发明人 徐波;陈旭;崔立忠;刘肖骢;郭宁明;杜向楠
分类号 G01R31/08(2006.01)I 主分类号 G01R31/08(2006.01)I
代理机构 乌鲁木齐合纵专利商标事务所 65105 代理人 董燕;周星莹
主权项 一种基于关联矩阵的配电网行波故障定位方法,其特征在于包括以下步骤:第一步,输入数据,进行数据预处理,包括(1)相模变换,采用Clark变换将三相电压变换获得α、β、0模量,选取α模量进行分析;(2)查找故障初始时刻,利用小波变换检测暂态电压信号突变点,作为各测量节点的故障初始时刻;(3)异常节点筛选,根据测量到的初始故障时刻及幅值对节点进行筛选,排除明显异常节点;第二步,构建关联矩阵,各主干节点构成节点集合N={n1,n2,n3,n4,n5,n6},两两节点关联配对构成n阶关联计算矩阵A,矩阵元素为节点两两配对双端测距计算结果,即:A<sub>12</sub>=(n1,n2),计算公式如下:A<sub>12</sub>=[l<sub>1</sub>-(t’<sub>2</sub>-t’<sub>1</sub>)×v]/2   (3)其中,t’<sub>1</sub>、t’<sub>2</sub>为主干节点的初始波头时刻,l<sub>1</sub>主干节点的间距离,v为定线路波速;当故障点位于节点n<sub>3</sub>、n<sub>4</sub>之间时,D<sub>2</sub>为故障点距主干节点n<sub>2</sub>的距离,A<sub>21</sub>≈0,A<sub>23</sub>≈l<sub>2</sub>,A<sub>24</sub>≈A<sub>25</sub>≈A<sub>26</sub>≈D<sub>2</sub>,矩阵元素小于测距标准差σ,则将矩阵元素设为零,其他矩阵元素计算方法同A<sub>12</sub>计算方法;令对角线元素A<sub>nn</sub>=0,则得出关联矩阵A,对应数值如下:<img file="FDA0001151001290000011.GIF" wi="1791" he="454" />其中关联阵行向量、列向量分别为:<maths num="0001"><math><![CDATA[<mrow><mi>A</mi><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>R</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><msub><mi>R</mi><mn>2</mn></msub></mtd></mtr><mtr><mtd><msub><mi>R</mi><mn>3</mn></msub></mtd></mtr><mtr><mtd><msub><mi>R</mi><mn>4</mn></msub></mtd></mtr><mtr><mtd><msub><mi>R</mi><mn>5</mn></msub></mtd></mtr><mtr><mtd><msub><mi>R</mi><mn>6</mn></msub></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>C</mi><mn>1</mn></msub></mtd><mtd><msub><mi>C</mi><mn>2</mn></msub></mtd><mtd><msub><mi>C</mi><mn>3</mn></msub></mtd><mtd><msub><mi>C</mi><mn>4</mn></msub></mtd><mtd><msub><mi>C</mi><mn>5</mn></msub></mtd><mtd><msub><mi>C</mi><mn>6</mn></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001151001290000012.GIF" wi="1862" he="447" /></maths>当监测线路上节点n数据异常,则关联矩阵中对应行向量R<sub>n</sub>及列向量C<sub>n</sub>异常;第三步,基于关联节点的故障区间判断,由公式4可知,关联矩阵行向量零元素最多两行就对应于故障点相邻节点,且两行向量是相邻的,之后进入第四步;第四步,判断向量元素是否满足相关性及对称性;若向量元素满足相关性及对称性,则进入第六步,若向量元素不满足相关性及对称性,则进入第五步;第五步,坏数据识别,在故障定位计算中,查找线路长度、波速v、初始波头时刻测量节点数据是否存在误差,确认坏数据存在的位置,之后进入第六步;第六步,误差矩阵构建及故障点定位,由关联矩阵A相对位置得出,存在以下相关性:(1)D<sub>1</sub>≈D<sub>2</sub>+l<sub>1</sub>,则定义测距计算结果G<sub>14</sub>=A<sub>14</sub>-l<sub>1</sub>-l<sub>2</sub>,G<sub>14</sub>≈G<sub>15</sub>≈G<sub>16</sub>;(2)D<sub>2</sub>≈D<sub>3</sub>+l<sub>2</sub>,则定义测距计算结果G<sub>24</sub>=A<sub>24</sub>-l<sub>1</sub>-l<sub>2</sub>,G<sub>24</sub>≈G<sub>25</sub>≈G<sub>26</sub>;对D<sub>4</sub>、D<sub>5</sub>、D<sub>6</sub>做(1)、(2)处理可将关联矩阵A转换为测量阵G,由于G<sub>1n</sub>=M<sub>n1</sub>,测量阵G可进行降阶处理,如下式所示:<maths num="0002"><math><![CDATA[<mrow><mi>M</mi><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>M</mi><mn>11</mn></msub></mtd><mtd><msub><mi>M</mi><mn>12</mn></msub></mtd><mtd><msub><mi>M</mi><mn>13</mn></msub></mtd></mtr><mtr><mtd><msub><mi>M</mi><mn>21</mn></msub></mtd><mtd><msub><mi>M</mi><mn>22</mn></msub></mtd><mtd><msub><mi>M</mi><mn>23</mn></msub></mtd></mtr><mtr><mtd><msub><mi>M</mi><mn>31</mn></msub></mtd><mtd><msub><mi>M</mi><mn>32</mn></msub></mtd><mtd><msub><mi>M</mi><mn>33</mn></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001151001290000021.GIF" wi="1579" he="222" /></maths> 其中M<sub>12</sub>=G<sub>14</sub>,即节点n<sub>1</sub>与n<sub>4</sub>配对计算得到的故障点位置,其他元素计算相同;由于误差整体上呈现正态分布原则,在无坏数据干扰情况下,最精确的配对计算结果与其他元素相比较总测量误差应最小,与其他元素相比,M<sub>11</sub>测量误差如下:<maths num="0003"><math><![CDATA[<mrow><msub><mi>E</mi><mn>11</mn></msub><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mo>|</mo><msub><mi>M</mi><mn>12</mn></msub><mo>-</mo><msub><mi>M</mi><mn>11</mn></msub><mo>|</mo></mrow></mtd><mtd><mrow><mo>|</mo><msub><mi>M</mi><mn>13</mn></msub><mo>-</mo><msub><mi>M</mi><mn>11</mn></msub><mo>|</mo></mrow></mtd></mtr><mtr><mtd><mrow><mo>|</mo><msub><mi>M</mi><mn>21</mn></msub><mo>-</mo><msub><mi>M</mi><mn>11</mn></msub><mo>|</mo></mrow></mtd><mtd><mrow><mo>|</mo><msub><mi>M</mi><mn>22</mn></msub><mo>-</mo><msub><mi>M</mi><mn>11</mn></msub><mo>|</mo></mrow></mtd><mtd><mrow><mo>|</mo><msub><mi>M</mi><mn>23</mn></msub><mo>-</mo><msub><mi>M</mi><mn>11</mn></msub><mo>|</mo></mrow></mtd></mtr><mtr><mtd><mrow><mo>|</mo><msub><mi>M</mi><mn>31</mn></msub><mo>-</mo><msub><mi>M</mi><mn>11</mn></msub><mo>|</mo></mrow></mtd><mtd><mrow><mo>|</mo><msub><mi>M</mi><mn>32</mn></msub><mo>-</mo><msub><mi>M</mi><mn>11</mn></msub><mo>|</mo></mrow></mtd><mtd><mrow><mo>|</mo><msub><mi>M</mi><mn>33</mn></msub><mo>-</mo><msub><mi>M</mi><mn>11</mn></msub><mo>|</mo></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001151001290000022.GIF" wi="1317" he="229" /></maths>则总的测量误差为:<img file="FDA0001151001290000023.GIF" wi="486" he="127" />构成误差矩阵:<maths num="0004"><math><![CDATA[<mrow><mi>S</mi><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>S</mi><mn>11</mn></msub></mtd><mtd><msub><mi>S</mi><mn>12</mn></msub></mtd><mtd><msub><mi>S</mi><mn>13</mn></msub></mtd></mtr><mtr><mtd><msub><mi>S</mi><mn>21</mn></msub></mtd><mtd><msub><mi>S</mi><mn>22</mn></msub></mtd><mtd><msub><mi>S</mi><mn>23</mn></msub></mtd></mtr><mtr><mtd><msub><mi>S</mi><mn>31</mn></msub></mtd><mtd><msub><mi>S</mi><mn>32</mn></msub></mtd><mtd><msub><mi>S</mi><mn>33</mn></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001151001290000024.GIF" wi="1373" he="223" /></maths> 误差矩阵中值最小元素即为精确配对结果,该结果为故障点对应线路主干节点的位置,当该结果与节点n<sub>3</sub>、n<sub>4</sub>距离大于σ,则认为故障点位于主干线路上;当故障点距离n<sub>3</sub>、n<sub>4</sub>距离小于σ,则认为故障点位于分支上,则节点n<sub>3</sub>、n<sub>4</sub>配对进行双端定位计算,故障点距离实际终端节点3计算公式如下:L=[l<sub>3</sub>+b<sub>3</sub>+b<sub>4</sub>-(t<sub>4</sub>-t<sub>3</sub>)×v]/2   (11)其中l<sub>3</sub>为主干节点的间距离,b<sub>3</sub>、b<sub>4</sub>分别为各分支线路长度,t<sub>3</sub>、t<sub>4</sub>为分支线路末端的终端节点的初始波头时刻,v为线路波速。
地址 841000 新疆维吾尔自治区巴音郭楞蒙古自治州库尔勒市石化大道电力公司812室
您可能感兴趣的专利