一種輸電線路的雷擊定位方法
【專利摘要】本發(fā)明涉及一種輸電線路的雷擊定位方法,屬于高電壓試驗設(shè)備及測量【技術(shù)領(lǐng)域】。本方法通過查詢雷電流在線監(jiān)測系統(tǒng)監(jiān)測數(shù)據(jù),獲取擊中桿塔的雷電流經(jīng)桿塔和輸電線路分流后,通過該監(jiān)測系統(tǒng)的雷電流傳感器的電流幅值、電流波形、上升時間、半波時間;通過查詢雷電定位系統(tǒng)監(jiān)測數(shù)據(jù),取得同時刻同區(qū)域的雷擊時間、雷擊區(qū)域、雷電流幅值、雷電回?fù)舸螖?shù)數(shù)據(jù);結(jié)合仿真計算得出的雷電波形經(jīng)過桿塔后的衰減關(guān)系,及雷擊點距離桿塔不同位置的電流分流關(guān)系,考慮雷電定位系統(tǒng)的誤差,得出雷擊點的位置區(qū)間。本方法定位精度高,有利于線路雷擊損壞的快速查找和修復(fù);并保證了數(shù)據(jù)有效性和可靠性,節(jié)省了輸電線路雷擊故障的查找和線路修復(fù)時間。
【專利說明】一種輸電線路的雷擊定位方法
【技術(shù)領(lǐng)域】
[0001]本發(fā)明涉及一種輸電線路的雷擊定位方法,屬于高電壓試驗設(shè)備及測量【技術(shù)領(lǐng)域】。
【背景技術(shù)】
[0002]據(jù)不完全統(tǒng)計,我國每年因雷擊造成人員傷亡達3000?4000人,財產(chǎn)損失更是達到驚人的200億元,雷電災(zāi)害已成為破壞性日趨嚴(yán)重的氣象災(zāi)害之一。雷電帶來的電磁暫態(tài)現(xiàn)象屬于電磁污染范圍,而電磁污染是世界公認(rèn)的繼水質(zhì)污染、大氣污染、噪聲污染之后的第四大污染。隨著工農(nóng)業(yè)生產(chǎn)和人民生活的現(xiàn)代化,對鐵路、航空、金融、電力、電信、電視、網(wǎng)絡(luò)等服務(wù)系統(tǒng)和設(shè)施可靠性及服務(wù)質(zhì)量要求也越來越高,人類活動對這些公用事業(yè)的依賴性也越來越大,對雷電防護提出了更高的要求,因此對雷電流準(zhǔn)確測量顯得越來越重要。
[0003]從電力系統(tǒng)的角度出發(fā),雷電已經(jīng)成為了電力系統(tǒng)故障的最主要原因,無論是線路、變電站還是用電設(shè)備,如何使它們能夠更好地抵御雷電和雷電感應(yīng)產(chǎn)生的電磁暫態(tài)過程的沖擊是電力系統(tǒng)雷電防護中最為關(guān)心的問題,這需要建立在對雷電的放電過程和雷電流的各項參數(shù)有充分認(rèn)識和了解的基礎(chǔ)上。
[0004]目前對雷電流特征參數(shù)的研究大都依靠模擬和仿真的方法,而對自然界的雷電流特征參數(shù),包括波形、幅值、上升時間、持續(xù)時間等都缺少實際數(shù)據(jù)。如何能夠獲得真實的雷電流數(shù)據(jù),這對于雷電放電過程的認(rèn)識和雷電流及其空間場的研究都十分有意義。而獲得自然界雷電參數(shù)最直接的方法就是測量,以現(xiàn)在最常見的人工火箭引雷或者自然雷直擊高大物體為例,都是將測量裝置的探頭置于雷電流入地通道上,而將主機置于通道附近或者室內(nèi)進行測量和記錄的。但由于雷電流幅值往往很高,基本為千安培量級,且上升沿往往只有幾百納秒到幾微秒,所以選擇什么樣的測量探頭能夠滿足大的動態(tài)測量范圍和寬的頻率響應(yīng)是最為關(guān)心的問題,羅氏線圈是目前比較通用的選擇,它可以滿足上述兩點的要求,而且它的性能已經(jīng)為很多實驗所證實。
[0005]因為雷電是一種隨機程度很高的自然現(xiàn)象,所以要大量的測量雷電流參數(shù),最好的方法就是大量的安裝測量裝置,但是受到各種實際條件的制約,這樣的做法顯然是不合適的,所以在測量裝置安裝地點的選擇上首先要考慮以往的雷電活動規(guī)律,特別是雷擊點的分布規(guī)律和雷擊次數(shù)的統(tǒng)計規(guī)律。根據(jù)不同的雷電活動強度可以將測量地點劃分為多雷擊區(qū)域和少雷擊區(qū)域,比如通信基站的信號塔,廣播電視發(fā)射塔這些獨立、高大的建筑物,以及容易遭受雷擊的輸電線路的局部和某些變電站,這些都屬于高雷擊區(qū)域。所以在這些地點安裝相對較復(fù)雜,數(shù)據(jù)記錄能力較強的測量裝置。而在其它的少雷擊區(qū)域可以選擇記錄功能較單一(只記錄幅值),結(jié)構(gòu)較為簡單,造價較低的測量裝置。這樣可以充分利用可用的測量資源,最大限度地測量和收集自然界雷電流的參數(shù)信息。
[0006]當(dāng)雷電擊中桿塔或輸電線路時,雷擊點所處的線路中,會有幅值很大的雷電流通過,可能造成線路損壞。同時,雷擊點臨近的絕緣子處可能會發(fā)生閃絡(luò),發(fā)生閃絡(luò)后,雷電流流過閃絡(luò)區(qū)域,會發(fā)出熱量燒蝕絕緣子,影響絕緣子的強度和絕緣特性。如能通過技術(shù)手段對雷擊桿塔或輸電線路進行定位,可以實現(xiàn)對由雷擊損壞的線路進行快速查找和修復(fù)。
【發(fā)明內(nèi)容】
[0007]本發(fā)明目的是提出一種輸電線路的雷擊定位方法,依據(jù)上述雷電定位系統(tǒng)和雷電流在線監(jiān)測系統(tǒng)給出的測量數(shù)據(jù),計算雷擊線路故障發(fā)生的雷擊位置,為電力系統(tǒng)由雷擊造成的線路損壞的快速查找和修復(fù)提供位置依據(jù)。
[0008]本發(fā)明提出的輸電線路雷擊定位方法,包括以下步驟:
[0009](I)從雷電流在線監(jiān)測系統(tǒng)中獲取雷電流在線監(jiān)測系統(tǒng)所安裝的桿塔的位置Pp雷擊時間^、以及擊中桿塔的雷電流經(jīng)桿塔和輸電線路分流后通過在線監(jiān)測系統(tǒng)中雷電流傳感器的雷電流特征參數(shù),雷電流特征參數(shù)包括雷電流的電流幅值I1、雷電流波形、上升時間tupl和半波時間thl,并從雷電流波形中,得到波峰次數(shù)N1 ;
[0010](2)從雷電定位系統(tǒng)中獲取監(jiān)測數(shù)據(jù),監(jiān)測數(shù)據(jù)包括同時刻同區(qū)域的雷擊時間t2、雷電流擊中區(qū)域P2、雷電流幅值I2和雷電回?fù)舸螖?shù)N2 ;
[0011](3)將上述雷電流在線監(jiān)測系統(tǒng)中獲取的雷擊時間、雷電流在線監(jiān)測系統(tǒng)所安裝的桿塔的位置P1、雷電流特征參數(shù)中的電流幅值I1和波形中的波峰次數(shù)N1,與雷電定位系統(tǒng)監(jiān)測數(shù)據(jù)中的雷擊時間t2、雷電流擊中區(qū)域P2、雷電流幅值I2、雷電回?fù)舸螖?shù)N2進行比對,若同時滿足h與t2相差I(lǐng)毫秒 以內(nèi)、P1位于區(qū)域P2內(nèi)、N1與N2相同,且I1小于I2,則判定雷電流在線監(jiān)測系統(tǒng)記錄的數(shù)據(jù)和雷電定位系統(tǒng)記錄的數(shù)據(jù)由同一雷電引起;
[0012](4)根據(jù)雷電擊中輸電線路時,輸電線路中,地線上的雷電流幅值1、上升時間tup和半波時間th與雷擊點和觀測點之間的桿塔數(shù)量及桿塔間線路的段數(shù)的理論關(guān)系表和理論關(guān)系曲線,將上述雷電流幅值I1與雷電流幅值I2的比值I1A2、上升時間tupl和半波時間thl代入上述理論關(guān)系表和理論關(guān)系曲線中,檢索到雷電流在線監(jiān)測系統(tǒng)所在位置與雷電流擊中輸電線路的位置之間的桿塔數(shù)量G和線路的段數(shù)D ;
[0013](5)通過對雷擊線路的仿真計算,得到雷電流沿線路傳播時,經(jīng)過G個桿塔和D段線路廣生幅值裳減后的雷電流幅值占裳減如的雷電流幅值的百分比Q,將上述雷電流在線監(jiān)測系統(tǒng)測量的雷電流幅值I1除以Q,得出衰減前距雷電流擊中線路位置最近處桿塔的雷電流幅值I。:10=Ii/Q ;
[0014](6)根據(jù)雷電擊中輸電線路桿塔之間的地線時,雷電流沿地線向兩側(cè)分流比例與雷擊點距兩側(cè)桿塔的距離之間的理論關(guān)系表,將上述雷電流幅值IjP I2的比例IcZI2代入該理論關(guān)系表中,檢索到雷電流擊中位置與最近桿塔的距離L,若L為0,則判定雷電流直接擊中桿塔,若L不為0,則判定雷電流擊中桿塔之間的線路;
[0015](7)利用下式,計算得到雷電流擊中輸電線路位置與雷電流在線監(jiān)測系統(tǒng)位置之間的距離S:S=GXA+L,其中,A為從電力系統(tǒng)的輸電線路數(shù)據(jù)中得到雷電流在線監(jiān)測系統(tǒng)所在位置處的桿塔之間的平均距離;
[0016](8)對上述雷電流幅值的測量值I2分別增加50%和減少50%后,分別重復(fù)步驟(4)至步驟(7),從得到的雷電流擊中輸電線路位置與雷電流在線監(jiān)測系統(tǒng)位置之間的距離S中,將最大值Smax和最小值Smin分別定義為雷電流擊中輸電線路位置與雷電流在線監(jiān)測系統(tǒng)之間距離的上限和下限;[0017](9)根據(jù)上述步驟(8)的最小值Smin和最大值Smax,得到雷電流擊中輸電線路位置所在的線段[SMIN,Smax],記該線段的中點為SAV,再根據(jù)上述雷電流在線監(jiān)測系統(tǒng)所安裝的桿塔的位置P1,得到雷電流擊中輸電線路的位置Ptl =Pci=PJSav,將Ptl與雷電定位系統(tǒng)監(jiān)測的雷擊區(qū)間P2進行比較,若Ptl處于區(qū)域P2內(nèi),則將Ptl作為雷擊線路的位置,用于雷擊故障修復(fù)和線路維護,若Ptl處于區(qū)域P2外,則結(jié)束定位計算。
[0018]本發(fā)明提出的輸電線路的雷擊定位方法,其優(yōu)點是:
[0019](I)本發(fā)明綜合了雷電定位系統(tǒng)和雷電流在線監(jiān)測系統(tǒng)給出的測量數(shù)據(jù),計算雷擊線路的位置和該雷擊的雷電流波形參數(shù),為電力系統(tǒng)由雷擊造成的線路損壞的快速查找和修復(fù)提供位置依據(jù),其定位精度高于單獨由雷電定位系統(tǒng)或雷電流在線監(jiān)測系統(tǒng)給出的雷電定位結(jié)果,定位誤差小于一至兩個桿塔的距離,對線路雷擊損壞的快速查詢和修復(fù),具有參考意義。
[0020](2)本發(fā)明提出的定位方法,考慮了雷電定位系統(tǒng)的測量誤差,通過雷電定位系統(tǒng)和雷電流在線監(jiān)測系統(tǒng)互相比對,減小定位誤差,得出雷電定位的距離區(qū)間,保證了數(shù)據(jù)的有效性和可靠性,節(jié)省了輸電線路雷擊故障的查找和線路修復(fù)時間。
【專利附圖】
【附圖說明】
[0021]圖1為雷電流擊中桿塔頂部或桿塔之間示意圖。
[0022]圖2為雷電流擊中桿塔之間線路不同位置示意圖。
[0023]圖3為輸電線路中地線上的雷電流幅值I與雷擊點和觀測點間的桿塔數(shù)量的關(guān)系O
[0024]圖4為輸電線路中地線上的雷電流上升時間tup與雷擊點和觀測點間的桿塔數(shù)量的關(guān)系。
[0025]圖5為輸電線路中地線上的雷電流半波時間th與雷擊點和觀測點間的桿塔數(shù)量的關(guān)系。
[0026]圖6為雷電流沿地線向一側(cè)分流比例與雷擊點相對該側(cè)桿塔距離的關(guān)系。
[0027]圖7為考慮雷電定位系統(tǒng)誤差后根據(jù)雷電流幅值判斷雷擊點的位置示意圖。
【具體實施方式】
[0028]本發(fā)明提出的輸電線路雷擊定位方法,包括以下步驟:
[0029](I)從雷電流在線監(jiān)測系統(tǒng)中獲取雷電流在線監(jiān)測系統(tǒng)所安裝的桿塔的位置Pp雷擊時間^、以及擊中桿塔的雷電流經(jīng)桿塔和輸電線路分流后通過在線監(jiān)測系統(tǒng)中雷電流傳感器的雷電流特征參數(shù),雷電流特征參數(shù)包括雷電流的電流幅值I1、雷電流波形、上升時間tupl和半波時間thl,并從雷電流波形中,得到波峰次數(shù)N1 ;
[0030](2)從雷電定位系統(tǒng)中獲取監(jiān)測數(shù)據(jù),監(jiān)測數(shù)據(jù)包括同時刻同區(qū)域的雷擊時間t2、雷電流擊中區(qū)域P2、雷電流幅值I2和雷電回?fù)舸螖?shù)N2 ;
[0031](3)將上述雷電流在線監(jiān)測系統(tǒng)中獲取的雷擊時間h、雷電流在線監(jiān)測系統(tǒng)所安裝的桿塔的位置P1、雷電流特征參數(shù)中的電流幅值I1和波形中的波峰次數(shù)N1,與雷電定位系統(tǒng)監(jiān)測數(shù)據(jù)中的雷擊時間t2、雷電流擊中區(qū)域P2、雷電流幅值I2、雷電回?fù)舸螖?shù)N2進行比對,若同時滿足h與t2相差I(lǐng)毫秒以內(nèi)、P1位于區(qū)域P2內(nèi)、N1與N2相同,且I1小于I2,則判定雷電流在線監(jiān)測系統(tǒng)記錄的數(shù)據(jù)和雷電定位系統(tǒng)記錄的數(shù)據(jù)由同一雷電引起;
[0032](4)根據(jù)雷電擊中輸電線路時,輸電線路中,地線上的雷電流幅值1、上升時間tup和半波時間th與雷擊點和觀測點之間的桿塔數(shù)量及桿塔間線路的段數(shù)的理論關(guān)系表和理論關(guān)系曲線,將上述雷電流幅值I1與雷電流幅值I2的比值I1A2、上升時間tupl和半波時間thl代入上述理論關(guān)系表和理論關(guān)系曲線中,檢索到雷電流在線監(jiān)測系統(tǒng)所在位置與雷電流擊中輸電線路的位置之間的桿塔數(shù)量G和線路的段數(shù)D ;
[0033](5)通過對雷擊線路的仿真計算,得到雷電流沿線路傳播時,經(jīng)過G個桿塔和D段線路廣生幅值裳減后的雷電流幅值占裳減如的雷電流幅值的百分比Q,將上述雷電流在線監(jiān)測系統(tǒng)測量的雷電流幅值I1除以Q,得出衰減前距雷電流擊中線路位置最近處桿塔的雷電流幅值I。:10=Ii/Q ;
[0034](6)根據(jù)雷電擊中輸電線路桿塔之間的地線時,雷電流沿地線向兩側(cè)分流比例與雷擊點距兩側(cè)桿塔的距離之間的理論關(guān)系表,將上述雷電流幅值IjP I2的比例IcZI2代入該理論關(guān)系表中,檢索到雷電流擊中位置與最近桿塔的距離L,若L為0,則判定雷電流直接擊中桿塔,若L不為0,則判定雷電流擊中桿塔之間的線路;
[0035](7)利用下式,計算得到雷電流擊中輸電線路位置與雷電流在線監(jiān)測系統(tǒng)位置之間的距離S:S=GXA+L,其中,A為從電力系統(tǒng)的輸電線路數(shù)據(jù)中得到雷電流在線監(jiān)測系統(tǒng)所在位置處的桿塔之間的平均距離;
[0036](8)對上述雷電流幅值的測量值I2分別增加50%和減少50%后,分別重復(fù)步驟(4)至步驟(7),從得到的雷電流擊中輸電線路位置與雷電流在線監(jiān)測系統(tǒng)位置之間的距離S中,將最大值Smax和最小值Smin分別定義為雷電流擊中輸電線路位置與雷電流在線監(jiān)測系統(tǒng)之間距離的上限和下限;
[0037](9)根據(jù)上述步驟(8)的最小值Smin和最大值Smax,得到雷電流擊中輸電線路位置所在的線段[SMIN,Smax],記該線段的中點為SAV,再根據(jù)上述雷電流在線監(jiān)測系統(tǒng)所安裝的桿塔的位置P1,得到雷電流擊中輸電線路的位置Ptl =Pci=PJSav,將Ptl與雷電定位系統(tǒng)監(jiān)測的雷擊區(qū)間P2進行比較,若Ptl處于區(qū)域P2內(nèi),則將Ptl作為雷擊線路的位置,用于雷擊故障修復(fù)和線路維護,若Ptl處于區(qū)域P2外,則結(jié)束定位計算。
[0038]本發(fā)明定位方法中,輸電線路中地線上的雷電流幅值、上升時間和半波時間與雷擊點與觀測點的距離的關(guān)系為,隨著雷電流經(jīng)過桿塔數(shù)量的增加,雷電波波形上升時間逐漸變短,幅值逐漸變小,半波時間變長,具體計算結(jié)果如下:
[0039]使用國家電網(wǎng)公司出版的《輸電線路設(shè)計標(biāo)準(zhǔn)手冊》作為仿真的設(shè)計標(biāo)準(zhǔn),其中桿塔采用二分裂式酒杯塔,整個輸電線路包含11個典型桿塔,線路末端各設(shè)置為IOOOkm (即線路為無窮長)。檔距均為400m。三相導(dǎo)線采用LGJ-50,水平排列,均距7.5m,兩塔頭同三相導(dǎo)線的垂直距離為5.7m,三相導(dǎo)線距離地面20m,導(dǎo)線為二分裂,分裂間距0.4572m,導(dǎo)線半徑0.02m,電阻率0.03984 Ω /km,垂馳12m ;接地線半徑0.0057m,電阻率為2.8645 Ω /km,垂馳7m。
[0040]雷電流擊中桿塔頂部或桿塔之間的示意圖和雷擊桿塔之間線路不同位置的示意圖分別如圖1和圖2所示。
[0041]雷電流采用標(biāo)準(zhǔn)雷電流波形進行計算,得出當(dāng)雷擊塔頂或桿塔之間的線路時,如圖1所示,通過仿真計算,得到通過臨近幾個桿塔的雷電流仿真波形,其總體規(guī)律如表I所/Jn ο
[0042]表1雷電波形經(jīng)過桿塔后衰減數(shù)據(jù)
[0043]
【權(quán)利要求】
1.一種輸電線路雷擊定位方法,其特征在于該方法包括以下步驟: (1)從雷電流在線監(jiān)測系統(tǒng)中獲取雷電流在線監(jiān)測系統(tǒng)所安裝的桿塔的位置P1、雷擊時間、以及擊中桿塔的雷電流經(jīng)桿塔和輸電線路分流后通過在線監(jiān)測系統(tǒng)中雷電流傳感器的雷電流特征參數(shù),雷電流特征參數(shù)包括雷電流的電流幅值I1、雷電流波形、上升時間tupl和半波時間thl,并從雷電流波形中,得到波峰次數(shù)N1 ; (2)從雷電定位系統(tǒng)中獲取監(jiān)測數(shù)據(jù),監(jiān)測數(shù)據(jù)包括同時刻同區(qū)域的雷擊時間t2、雷電流擊中區(qū)域P2、雷電流幅值I2和雷電回?fù)舸螖?shù)N2 ; (3)將上述雷電流在線監(jiān)測系統(tǒng)中獲取的雷擊時間h、雷電流在線監(jiān)測系統(tǒng)所安裝的桿塔的位置P1、雷電流特征參數(shù)中的電流幅值I1和波形中的波峰次數(shù)N1,與雷電定位系統(tǒng)監(jiān)測數(shù)據(jù)中的雷擊時間&、雷電流擊中區(qū)域P2、雷電流幅值I2、雷電回?fù)舸螖?shù)N2進行比對,若同時滿足h與t2相差I(lǐng)毫秒以內(nèi)、P1位于區(qū)域P2內(nèi)、N1與N2相同,且I1小于I2,則判定雷電流在線監(jiān)測系統(tǒng)記錄的數(shù)據(jù)和雷電定位系統(tǒng)記錄的數(shù)據(jù)由同一雷電引起; (4)根據(jù)雷電擊中輸電線路時,輸電線路中,地線上的雷電流幅值1、上升時間tup和半波時間th與雷擊點和觀測點之間的桿塔數(shù)量及桿塔間線路的段數(shù)的理論關(guān)系表和理論關(guān)系曲線,將上述雷電流幅值I1與雷電流幅值I2的比值I1A2、上升時間tupl和半波時間thlR入上述理論關(guān)系表和理論關(guān)系曲線中,檢索到雷電流在線監(jiān)測系統(tǒng)所在位置與雷電流擊中輸電線路的位置之間的桿塔數(shù)量G和線路的段數(shù)D ; (5)通過對雷擊線路的仿真計算,得到雷電流沿線路傳播時,經(jīng)過G個桿塔和D段線路產(chǎn)生幅值衰減后的雷電流幅值占衰減前的雷電流幅值的百分比Q,將上述雷電流在線監(jiān)測系統(tǒng)測量的雷電流幅值I1除以Q,得出衰減前距雷電流擊中線路位置最近處桿塔的雷電流幅值 10 =10=VQ ; (6)根據(jù)雷電擊中輸電線路桿塔之間的地線時,雷電流沿地線向兩側(cè)分流比例與雷擊點距兩側(cè)桿塔的距離之間的理論關(guān)系表,將上述雷電流幅值Itl和I2的比例IcZI2代入該理論關(guān)系表中,檢索到雷電流擊中位置與最近桿塔的距離L,若L為O,則判定雷電流直接擊中桿塔,若L不為O,則判定雷電流擊中桿塔之間的線路; (7)利用下式,計算得到雷電流擊中輸電線路位置與雷電流在線監(jiān)測系統(tǒng)位置之間的距離S:S=GXA+L,其中,A為從電力系統(tǒng)的輸電線路數(shù)據(jù)中得到雷電流在線監(jiān)測系統(tǒng)所在位置處的桿塔之間的平均距離; (8)對上述雷電流幅值的測量值I2分別增加50%和減少50%后,分別重復(fù)步驟(4)至步驟(7),從得到的雷電流擊中輸電線路位置與雷電流在線監(jiān)測系統(tǒng)位置之間的距離S中,將最大值Smax和最小值Smin分別定義為雷電流擊中輸電線路位置與雷電流在線監(jiān)測系統(tǒng)之間距離的上限和下限; (9)根據(jù)上述步驟(8)的最小值Smin和最大值Smax,得到雷電流擊中輸電線路位置所在的線段[SMIN,Smax],記該線段的中點為SAV,再根據(jù)上述雷電流在線監(jiān)測系統(tǒng)所安裝的桿塔的位置P1,得到雷電流擊中輸電線路的位置Ptl =Pci=PJSav,將Ptl與雷電定位系統(tǒng)監(jiān)測的雷擊區(qū)間P2進行比較,若Ptl處于區(qū)域P2內(nèi),則將Ptl作為雷擊線路的位置,用于雷擊故障修復(fù)和線路維護,若Ptl處于區(qū)域P2外,則結(jié)束定位計算。
【文檔編號】G01R31/08GK103543386SQ201310483383
【公開日】2014年1月29日 申請日期:2013年10月16日 優(yōu)先權(quán)日:2013年10月16日
【發(fā)明者】余占清, 王運超, 王澤眾, 彭瓏, 李 雨, 趙媛, 陳水明, 何金良, 曾嶸, 端木林楠, 牛曉民, 王萍 申請人:清華大學(xué), 國家電網(wǎng)公司, 華北電力科學(xué)研究院有限責(zé)任公司, 華北電網(wǎng)有限公司