发明名称 一种二维雷电先导发展模型在山区地形下电位的获取方法
摘要 一种二维雷电先导发展模型在山区地形下电位的获取方法,包括:通过地理信息系统提取输电线路的实际地形数据;对输电线路的实际地形进行拟合,使用获得的实际地形数据求出对应的保角变换函数中的与地形相关的参数,使用保角变换函数将山区地形变换到平原地形,山顶地形下采用的保角变换的函数为:<img file="DDA0000084661570000011.tif" wi="290" he="77" />山谷地形下采用的保角变换的函数为,<img file="DDA0000084661570000012.tif" wi="290" he="77" />其中z=x+i(y+d1),z是保角变换前复平面点的坐标;w=u+i(v+d1),w是保角变换后复平面点的坐标,d和d1是所述与地形相关的参数;保角变换后先做先导通道的镜像,然后根据镜像法计算空间任意一点的电位,空间任意一点S(x,y)的电位的计算公式如下:<img file="DDA0000084661570000013.tif" wi="1251" he="114" /></maths>本发明解决了山区地形下先导发展模型中电位的计算问题。
申请公布号 CN102955055B 申请公布日期 2014.11.26
申请号 CN201110238999.X 申请日期 2011.08.19
申请人 湖北省电力公司电力科学研究院;武汉大学;国家电网公司 发明人 姚尧;阮羚;周文俊;喻剑辉;唐泽洋;李涵
分类号 G01R19/00(2006.01)I 主分类号 G01R19/00(2006.01)I
代理机构 武汉楚天专利事务所 42113 代理人 孔敏
主权项 一种二维雷电先导发展模型在山区地形下电位的获取方法,其特征在于:包括步骤10:通过地理信息系统提取输电线路的实际地形数据;步骤20:对输电线路的实际地形进行拟合,使用步骤10中获得的实际地形数据求出对应的保角变换函数中的与地形相关的参数,然后使用保角变换函数将山区地形变换到平原地形,山顶地形下采用的保角变换的函数为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>w</mi><mo>=</mo><msqrt><msup><mi>z</mi><mn>2</mn></msup><mo>+</mo><msup><mi>d</mi><mn>2</mn></msup></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000543553670000011.GIF" wi="1348" he="101" /></maths>山谷地形下采用的保角变换的函数为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>w</mi><mo>=</mo><msqrt><msup><mi>z</mi><mn>2</mn></msup><mo>-</mo><msup><mi>d</mi><mn>2</mn></msup></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000543553670000012.GIF" wi="1346" he="98" /></maths>式(1)和(2)中z=x+i(y+d1),z是保角变换前复平面点的坐标;w=u+i(v+d1),w是保角变换后复平面点的坐标,d和d1是所述与地形相关的参数;步骤30:保角变换后先做先导通道的镜像,然后根据镜像法计算空间任意一点的电位,空间任意一点S(x,y)的电位的计算公式如下所示:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>V</mi><mi>S</mi></msub><mo>=</mo><mfrac><mn>1</mn><msub><mrow><mn>4</mn><mi>&pi;&epsiv;</mi></mrow><mn>0</mn></msub></mfrac><mrow><mo>(</mo><msubsup><mo>&Integral;</mo><msub><mi>y</mi><mn>10</mn></msub><msub><mi>y</mi><mn>1</mn></msub></msubsup><mi>TH</mi><mrow><mo>(</mo><msub><mi>z</mi><mn>0</mn></msub><mo>)</mo></mrow><msub><mi>r</mi><mn>1</mn></msub><mi>dy</mi><mo>+</mo><msubsup><mo>&Integral;</mo><msub><mi>y</mi><mn>20</mn></msub><msub><mi>y</mi><mn>2</mn></msub></msubsup><mi>TH</mi><mrow><mo>(</mo><msub><mi>z</mi><mn>0</mn></msub><mo>)</mo></mrow><msub><mi>r</mi><mn>1</mn></msub><mi>dy</mi><mo>+</mo><mi>k</mi><mo>+</mo><msubsup><mo>&Integral;</mo><msub><mi>y</mi><mrow><mi>n</mi><mn>0</mn></mrow></msub><msub><mi>y</mi><mi>n</mi></msub></msubsup><mi>TH</mi><mrow><mo>(</mo><msub><mi>z</mi><mn>0</mn></msub><mo>)</mo></mrow><msub><mi>r</mi><mn>1</mn></msub><mi>dy</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>12</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000543553670000013.GIF" wi="1585" he="140" /></maths>式(12)中ξ<sub>0</sub>是真空介电常数,y<sub>n0</sub>为先导通道第n段的末端纵坐标,y<sub>n</sub>为先导通道第n段的起始纵坐标,TH(z<sub>0</sub>)=0.3α+0.7β,<img file="FDA0000543553670000014.GIF" wi="593" he="133" />H为雷云底部距离地面的高度,z<sub>0</sub>为先导头部高度(在公式中z<sub>0</sub>&gt;10),<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>r</mi><mn>1</mn></msub><mo>=</mo><mi>ln</mi><mfrac><msqrt><msup><mrow><mo>(</mo><mi>u</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>-</mo><mi>u</mi><mrow><mo>(</mo><msub><mi>x</mi><mn>0</mn></msub><mo>,</mo><msub><mi>y</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mi>v</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>+</mo><mi>v</mi><mrow><mo>(</mo><msub><mi>x</mi><mn>0</mn></msub><mo>,</mo><msub><mi>y</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup></msqrt><msqrt><msup><mrow><mo>(</mo><mi>u</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>-</mo><mi>u</mi><mrow><mo>(</mo><msub><mi>x</mi><mn>0</mn></msub><mo>,</mo><msub><mi>y</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mi>v</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>-</mo><mi>v</mi><mrow><mo>(</mo><msub><mi>x</mi><mn>0</mn></msub><mo>,</mo><msub><mi>y</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mfrac></mrow>]]></math><img file="FDA0000543553670000021.GIF" wi="1045" he="183" /></maths><maths num="0005" id="cmaths0005"><math><![CDATA[<mfenced open='' close=''><mtable><mtr><mtd><mi>T</mi><mo>=</mo><msub><mi>a</mi><mn>0</mn></msub><mfrac><mrow><msub><mi>L</mi><mn>1</mn></msub><mo>+</mo><msub><mi>L</mi><mn>2</mn></msub><mo>+</mo><mi>k</mi><mo>+</mo><msub><mi>L</mi><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><msqrt><mn>1</mn><mo>+</mo><mfrac><mn>1</mn><msup><mi>k</mi><mn>2</mn></msup></mfrac></msqrt><mrow><mo>(</mo><msub><mi>y</mi><mi>n</mi></msub><mo>-</mo><mi>y</mi><mo>)</mo></mrow></mrow><mrow><mi>H</mi><mrow><mo>(</mo><msub><mi>z</mi><mn>0</mn></msub><mo>)</mo></mrow></mrow></mfrac><msub><mi>I</mi><mn>0</mn></msub><mo>+</mo></mtd></mtr><mtr><mtd><mfrac><mrow><msub><mi>I</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>a</mi><mo>+</mo><mi>b</mi><mrow><mo>(</mo><msub><mi>L</mi><mi>Z</mi></msub><mo>-</mo><mrow><mo>(</mo><msub><mi>L</mi><mn>1</mn></msub><mo>+</mo><msub><mi>L</mi><mn>2</mn></msub><mo>+</mo><mi>k</mi><mo>+</mo><msub><mi>L</mi><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mo>-</mo><msqrt><mn>1</mn><mo>+</mo><mfrac><mn>1</mn><msup><mi>k</mi><mn>2</mn></msup></mfrac></msqrt><mrow><mo>(</mo><msub><mi>y</mi><mi>n</mi></msub><mo>-</mo><mi>y</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>)</mo></mrow></mrow><mrow><mn>1</mn><mo>+</mo><mi>c</mi><mrow><mo>(</mo><msub><mi>L</mi><mi>Z</mi></msub><mo>-</mo><mrow><mo>(</mo><msub><mi>L</mi><mn>1</mn></msub><mo>+</mo><msub><mi>L</mi><mn>2</mn></msub><mo>+</mo><mi>k</mi><mo>+</mo><msub><mi>L</mi><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mo>-</mo><msqrt><mn>1</mn><mo>+</mo><mfrac><mn>1</mn><msup><mi>k</mi><mn>2</mn></msup></mfrac></msqrt><mrow><mo>(</mo><msub><mi>y</mi><mi>n</mi></msub><mo>-</mo><mi>y</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>+</mo><mi>d</mi><msup><mrow><mo>(</mo><msub><mi>L</mi><mi>Z</mi></msub><mo>-</mo><mrow><mo>(</mo><msub><mi>L</mi><mn>1</mn></msub><mo>+</mo><msub><mi>L</mi><mn>2</mn></msub><mo>+</mo><mi>k</mi><mo>+</mo><msub><mi>L</mi><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mo>-</mo><msqrt><mn>1</mn><mo>+</mo><mfrac><mn>1</mn><msup><mi>k</mi><mn>2</mn></msup></mfrac></msqrt><mrow><mo>(</mo><msub><mi>y</mi><mi>n</mi></msub><mo>-</mo><mi>y</mi><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup></mrow></mfrac></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000543553670000022.GIF" wi="1668" he="441" /></maths>其中,I<sub>0</sub>为回击电流峰值,a<sub>0</sub>=1.476×10<sup>‑5</sup>,a=4.857×10<sup>‑5</sup>,b=3.9097×10<sup>‑6</sup>,c=0.522,d=3.73×10<sup>‑3</sup>,k为保角变换前直线的斜率,L<sub>1</sub>、L<sub>2</sub>、L<sub>n‑1</sub>分别表示的第一段、第二段、第n‑1段先导通道的长度。
地址 430077 湖北省武汉市武昌区徐东大街361号