发明名称 一种双极直流输电线路无线电干扰确定方法
摘要 本发明公开了一种双极直流输电线路无线电干扰确定方法,步骤是:1、输电线路模型参数的获取:输电线路导线参数和输电线路各导线所处环境;双极直流输电线路无线电干扰的方程组确定;2、双极直流输电线路无线电干扰确定方法:包括双极直流输电线路无线电干扰的模型建立;双极直流输电线路无线电干扰的方程组确定;各分裂导线的中心未知电荷量确定;各分裂导线内等效电荷位置和电荷量的确定;双极直流输电线路无线电干扰确定。本发明可确定直流输电输电线路无线电干扰,可判断好天气下直流输电输电线路无线电干扰不应大于55dB(μV/m)。可运用于在双极直流输电线路设计、运行和改造工程。
申请公布号 CN102590686A 申请公布日期 2012.07.18
申请号 CN201210060642.1 申请日期 2012.03.08
申请人 中国南方电网有限责任公司超高压输电公司检修试验中心 发明人 邓军;肖遥;吕金壮;陈禾;楚金伟;王竣;陆春玉;骆晓龙;杨方明;张飚
分类号 G01R31/00(2006.01)I 主分类号 G01R31/00(2006.01)I
代理机构 广州科粤专利商标代理有限公司 44001 代理人 黄培智
主权项 1.一种双极直流输电线路无线电干扰确定方法,其特征在于,包括如下步骤:S1、输电线路模型参数的获取:输电线路导线参数和输电线路各导线所处环境;所述输电线路导线参数包含单根导线的分裂数、分裂半径和导线的截面积;两导线之间的距离;所述输电线路导线所处环境包括输电线路对地投影为起点,测量各导线对地高度;S2、双极直流输电线路无线电干扰确定方法S21,双极直流输电线路无线电干扰的模型建立:将测量的输电线路导线参数和输电线路各导线所处环境构建物理模型,绘制输电线路无线电干扰的二维图形模型,包括:双极直流输电线路无线电干扰的模型的二维图形绘制包括根据换单根导线的分裂数、分裂半径和导线的截面积;两导线之间的距离绘制的二维图形;双极直流输电线路无线电干扰的模型的二维图形包括根据测量各导线对地高度绘制各导线相对于地面的镜像;根据双极直流输电线路无线电干扰的模型得到每根分裂导线及其镜像的圆心坐标(x<sub>i</sub>,y<sub>i</sub>)、导线半径r<sub>i</sub>,并将这些信息存储在矩阵A,所述矩阵A每一行代表一根分裂导线的信息,第一列代表分裂导线编号,第二列代表分裂导线在地面的水平距离,第三列代表分裂导线对地高度,第四列代表分裂导线半径,矩阵A的总行数等于分裂导线及镜像的总和N;S22,双极直流输电线路无线电干扰的方程组确定首先设置每根分裂导线作为孤立导体对待,且认为每根子导线上电荷相等并集中在该分裂导线的中心,因此,给定每个分裂导线的电荷初始向量,并保存在数组q<sub>0</sub>,q<sub>0</sub>的大小等于之前所述N;根据双极直流输电线路无线电干扰的模型和信息存储矩阵A得到双极直流输电线路无线电干扰的方程组如式(1);<maths num="0001"><![CDATA[<math><mrow><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>u</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><msub><mi>u</mi><mn>2</mn></msub></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>u</mi><mi>N</mi></msub></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>D</mi><mn>11</mn></msub></mtd><mtd><msub><mi>D</mi><mn>12</mn></msub></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msub><mi>D</mi><mrow><mn>1</mn><mi>N</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>D</mi><mn>21</mn></msub></mtd><mtd><msub><mi>D</mi><mn>22</mn></msub></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>D</mi><mrow><mi>N</mi><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>D</mi><mrow><mi>N</mi><mn>2</mn></mrow></msub></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msub><mi>D</mi><mi>NN</mi></msub></mtd></mtr></mtable></mfenced><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>q</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><msub><mi>q</mi><mn>2</mn></msub></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>q</mi><mn>2</mn></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中向量u为分裂导线运行电压,向量q为各分裂导线的未知电荷量,矩阵D是电位系数矩阵。所述向量u是由u<sub>1</sub>,u<sub>2</sub>…u<sub>N</sub>组成且因分裂导线运行电压已知而该向量已知,向量q是由q<sub>1</sub>,q<sub>2</sub>…q<sub>N</sub>组成且为未知量,分裂导线运行电压是指双极直流输电线路正常运行时导线上的电位值,矩阵D是由D<sub>11</sub>,D<sub>12</sub>…D<sub>NN</sub>组成,其中矩阵D的元素由式(2)计算。<maths num="0002"><![CDATA[<math><mrow><msub><mi>D</mi><mi>ij</mi></msub><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mfrac><mn>1</mn><mrow><mn>2</mn><mi>&pi;</mi><msub><mi>&epsiv;</mi><mn>0</mn></msub></mrow></mfrac><mi>ln</mi><mrow><mo>(</mo><mfrac><mrow><mn>2</mn><msub><mi>h</mi><mi>i</mi></msub></mrow><msub><mi>R</mi><mi>i</mi></msub></mfrac><mo>)</mo></mrow><mi>i</mi><mo>=</mo><mi>j</mi></mtd></mtr><mtr><mtd><mfrac><mn>1</mn><mrow><mn>2</mn><mi>&pi;</mi><msub><mi>&epsiv;</mi><mn>0</mn></msub></mrow></mfrac><mi>ln</mi><mrow><mo>(</mo><mfrac><msubsup><mi>L</mi><mi>ij</mi><mo>&prime;</mo></msubsup><msub><mi>L</mi><mi>ij</mi></msub></mfrac><mo>)</mo></mrow><mi>i</mi><mo>&NotEqual;</mo><mi>j</mi></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中ε<sub>0</sub>是真空的介电常数<img file="FDA0000141805020000023.GIF" wi="236" he="106" />R<sub>i</sub>为分裂导线半径,L′<sub>ij</sub>是分裂导线i和分裂导线j的镜像之间的距离,L<sub>ij</sub>是分裂导线i和分裂导线j之间的距离;S23,各分裂导线的中心未知电荷量确定根据所述双极直流输电线路无线电干扰的方程组,进行如下基本流程:a)根据所述节点初始向量q<sub>0</sub>计算初始残值r<sub>0</sub>如式(3);r<sub>0</sub>=Dq<sub>0</sub>-u    (3)b)利用式(4)求解共轭梯度向量p<sub>k+1</sub>、<img file="FDA0000141805020000024.GIF" wi="75" he="57" />和向量β<sub>k</sub>;<img file="FDA0000141805020000025.GIF" wi="1391" he="669" />c)换流变磁场分析的未知向量C的确定和误差计算如(5);<maths num="0003"><![CDATA[<math><mrow><msub><mi>&alpha;</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>=</mo><mrow><mo>(</mo><msub><mi>z</mi><mi>k</mi></msub><mo>,</mo><msubsup><mi>r</mi><mi>k</mi><mo>*</mo></msubsup><mo>)</mo></mrow><mo>/</mo><mrow><mo>(</mo><msubsup><mi>p</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mo>*</mo></msubsup><mo>,</mo><msub><mi>Dp</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow></mrow></math>]]></maths>q<sub>k+1</sub>=q<sub>k</sub>+α<sub>k+1</sub>p<sub>k+1</sub>    (5)<maths num="0004"><![CDATA[<math><mrow><msub><mi>r</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>=</mo><msub><mi>r</mi><mi>k</mi></msub><mo>-</mo><msub><mi>&alpha;</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><msub><mi>Dp</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>;</mo><msubsup><mi>r</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mo>*</mo></msubsup><mo>=</mo><msubsup><mi>r</mi><mi>k</mi><mo>*</mo></msubsup><mo>-</mo><msub><mi>&alpha;</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><msup><mi>D</mi><mi>T</mi></msup><msubsup><mi>p</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mo>*</mo></msubsup></mrow></math>]]></maths>d)若r<sub>k+1</sub>的值大于10<sup>-7</sup>,q<sub>k+1</sub>作为新的q<sub>0</sub>带入a)步骤,重新计算q<sub>k+1</sub>和r<sub>k+1</sub>,直到r<sub>k+1</sub>的值小于10<sup>-7</sup>;否则,将q<sub>k+1</sub>和r<sub>k+1</sub>作为各分裂导线的中心未知电荷量的未知向量的解和误差;S24,各分裂导线内等效电荷位置和电荷量的确定根据所属各分裂导线的中心未知电荷量确定,假设分裂导线m,其他分裂导线在分裂导线m内的等效电荷位置和电荷量可采用镜像法确定。其他分裂导线在分裂导线m内的等效电荷量等于原电荷,但符号相反;其他分裂导线在分裂导线m内的等效电荷位置位于K<sub>mj</sub>,其计算方法如式(6);<maths num="0005"><![CDATA[<math><mrow><msub><mi>K</mi><mi>mj</mi></msub><mo>=</mo><mfrac><msubsup><mi>R</mi><mi>m</mi><mn>2</mn></msubsup><msub><mi>L</mi><mi>mj</mi></msub></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow></math>]]></maths>式(6)中参数与式(2)定义一样,j是其他分裂导线需要将电荷量等效到分裂导线m内的导线编号;S25)双极直流输电线路无线电干扰确定根据所述各分裂导线内等效电荷位置和电荷量,计算分裂导线几何中心与子导线中心连线方向上的子导线外表面上的电场φ<sub>n</sub>如式(7);<maths num="0006"><![CDATA[<math><mrow><msub><mi>&phi;</mi><mi>n</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mfrac><mn>1</mn><mrow><mn>2</mn><mi>&pi;</mi><msub><mi>&epsiv;</mi><mn>0</mn></msub></mrow></mfrac><mi>ln</mi><mrow><mo>(</mo><mfrac><msubsup><mi>L</mi><mi>ni</mi><mo>&prime;</mo></msubsup><msub><mi>L</mi><mi>ni</mi></msub></mfrac><mo>)</mo></mrow><msub><mi>q</mi><mi>i</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中q<sub>i</sub>表示分裂导线内i点处电荷量,L′<sub>ni</sub>和L<sub>ni</sub>的定义方法如式(2)。根据电场φ<sub>n</sub>计算双极直流输电线路无线电干扰E dBV/m;<maths num="0007"><![CDATA[<math><mrow><mi>E</mi><mo>=</mo><mn>38</mn><mo>+</mo><mn>1.6</mn><mrow><mo>(</mo><msub><mi>&phi;</mi><mi>n</mi></msub><mo>-</mo><mn>24</mn><mo>)</mo></mrow><mo>+</mo><mn>46</mn><mi>lgr</mi><mo>+</mo><mn>54</mn><mi>gn</mi><mo>+</mo><mi>&Delta;</mi><msub><mi>E</mi><mi>f</mi></msub><mo>+</mo><mn>33</mn><mi>lg</mi><mfrac><mn>20</mn><mi>D</mi></mfrac><mo>+</mo><mi>&Delta;</mi><msub><mi>E</mi><mi>W</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中r为分裂导线半径;n为导线分裂数;D为计算点到导线表面的距离;ΔE<sub>W</sub>为气象修正项,好天气取0 dBV/m,坏天气取3 dB V/m;ΔE<sub>f</sub>为干扰频率修正项,ΔE<sub>f</sub>=5[1-2(log10f)<sup>2</sup>],f为干扰频率,适用于0.15~30MHz频段。
地址 510663 广东省广州市萝岗区科学城科学大道181号A4栋检修试验中心