王漢堃 秦兆軍 金 師

?

呼包鄂地區(qū)閃電定位資料誤差分析及場(chǎng)地誤差優(yōu)化

王漢堃1秦兆軍1金 師2

1.內(nèi)蒙古自治區(qū)雷電預(yù)警防護(hù)中心，內(nèi)蒙古 呼和浩特 010051 2.呼和浩特市氣象局，內(nèi)蒙古 呼和浩特 010051

對(duì)目前常用的磁定向法誤差來源進(jìn)行分析，通過對(duì)現(xiàn)有的呼包鄂地區(qū)2010-2014年閃電定位數(shù)據(jù)進(jìn)行篩選，確定包頭青山、東勝、和林格爾三個(gè)分析站點(diǎn)，建立分析樣本，利用陳明理等基于Orville的本征值技術(shù)發(fā)展完善的場(chǎng)地誤差優(yōu)化模式進(jìn)行優(yōu)化，用對(duì)應(yīng)的雷達(dá)資料和雷災(zāi)資料與訂正后的定位資料進(jìn)行對(duì)比。驗(yàn)證表明：經(jīng)過場(chǎng)地誤差優(yōu)化后原本在強(qiáng)回波區(qū)域外的定位數(shù)據(jù)返回了強(qiáng)回波區(qū)域中，明顯更符合實(shí)際情況，場(chǎng)地優(yōu)化模式可以顯著提高閃電定位資料的準(zhǔn)確性，對(duì)今后的雷電監(jiān)測(cè)預(yù)警、雷災(zāi)調(diào)查等工作都有積極的意義。

閃電定位；誤差；優(yōu)化

呼包鄂地區(qū)位于內(nèi)蒙古自治區(qū)中西部的核心區(qū)域，包括包頭市、呼和浩特市、鄂爾多斯市，呈品字形分布于黃河兩岸，是內(nèi)蒙古最具活力的城市經(jīng)濟(jì)圈，被譽(yù)為內(nèi)蒙古的“金三角”地區(qū)。呼包鄂地區(qū)地型復(fù)雜，北部大青山和東南部蠻漢山為山地地形，黃河從中穿過，南部以平原為主，眾多因素導(dǎo)致了呼包鄂地區(qū)雷電頻發(fā)，年平均雷暴日在25天以上，屬于中雷區(qū)，是全國(guó)防雷減災(zāi)工作的重點(diǎn)地區(qū)之一。近幾十年來，由于對(duì)閃電現(xiàn)象的深入研究以及遙感遙測(cè)等技術(shù)的發(fā)展，在閃電探測(cè)特別是對(duì)地閃的監(jiān)測(cè)定位方面取得很大的進(jìn)展[2]。閃電定位系統(tǒng)的基本原理是依靠接收地閃回?fù)舭l(fā)出的甚低頻信號(hào)（VLF）來實(shí)現(xiàn)定位的。由于安裝場(chǎng)地存在差異使得閃電定位系統(tǒng)在使用中存在較大誤差，對(duì)閃電定位系統(tǒng)在實(shí)際工作中的應(yīng)用產(chǎn)生了較大的影響，這就需要對(duì)這種定位誤差的來源和處理方法進(jìn)行深入研究，得出有效解決閃電定位系統(tǒng)存在探測(cè)誤差的方法。呼包鄂地區(qū)采用的是北京華云東方探測(cè)技術(shù)有限公司開發(fā)的ADTD閃電定位系統(tǒng)，自2009年運(yùn)行以來一直未開展系統(tǒng)的探測(cè)效率、誤差分析等相關(guān)研究，閃電定位資料雖然已廣泛的應(yīng)用于雷擊災(zāi)害風(fēng)險(xiǎn)評(píng)估、雷電災(zāi)害監(jiān)測(cè)預(yù)警、雷災(zāi)技術(shù)調(diào)查等眾多防雷業(yè)務(wù)中，但定位資料一直不夠準(zhǔn)確，給日常的防雷減災(zāi)業(yè)務(wù)帶來了不必要的損失，閃電定位資料的準(zhǔn)確性已成為制約其在相關(guān)業(yè)務(wù)中進(jìn)一步應(yīng)用的重要因素。

1 閃電定位資料誤差分析

目前國(guó)內(nèi)外常用的閃電定位方法有時(shí)差法、磁定向法及混合法，本文僅對(duì)磁定向法的誤差進(jìn)行分析和探討。磁定向法是通過各方位探測(cè)器（DF）測(cè)量指向閃電位置的方位角，進(jìn)而進(jìn)行球面三角交匯來確定閃電位置的。

由于各種隨機(jī)因素及站址周圍場(chǎng)地環(huán)境的影響使實(shí)測(cè)方向具有兩種不同誤差，即隨機(jī)誤差和場(chǎng)地誤差。隨機(jī)誤差只有1°左右。但因站址而引起的場(chǎng)地誤差卻可達(dá)5°～15°，甚至更高。[3]場(chǎng)地誤差可看作兩部分組成，一部分是由于DF本身天線的原因造成的旋轉(zhuǎn)誤差，這部分誤差為常數(shù)；另一部分是由于測(cè)量場(chǎng)地周邊環(huán)境因素的影響造成的誤差，這部分誤差是方位角的函數(shù)[4]。由于旋轉(zhuǎn)誤差為常數(shù)，比較容易訂正，本文將著重圍繞后一部分場(chǎng)地誤差的優(yōu)化方法展開，為方便描述，后文中所有提到場(chǎng)地誤差除特別之外的全部為不包含旋轉(zhuǎn)誤差的場(chǎng)地誤差。

2 基于本征值技術(shù)的場(chǎng)地誤差優(yōu)化模介紹

基于Orville的本征值技術(shù)[5]，1990年陳明理等發(fā)展完善了一套場(chǎng)地誤差優(yōu)化模式，對(duì)北京地區(qū)的閃電定位資料進(jìn)行了優(yōu)化處理，確定出場(chǎng)地誤差曲線并優(yōu)化閃電位置，發(fā)現(xiàn)該模式具有明顯的場(chǎng)地誤差處理功能[6]。

2011年張其林等利用該模式對(duì)甘肅平?jīng)龅貐^(qū)的閃電定位資料進(jìn)行了優(yōu)化和驗(yàn)證，收到良好的效果，證明該模式對(duì)于山地地形也有一定的適用性。但是，由于復(fù)雜地形的影響，不同地區(qū)的閃電定位系統(tǒng)誤差各不相同，該模式在呼包鄂地區(qū)的適用性有待進(jìn)一步驗(yàn)證。此外，本文采用了呼包鄂地區(qū)2011—2014年近200萬條資料，對(duì)于場(chǎng)地誤差的優(yōu)化可靠性更強(qiáng)。[7]

2.1 不考慮場(chǎng)地誤差下的優(yōu)化模式

在不考慮場(chǎng)地誤差的情況下，閃電定位資料只存在隨機(jī)誤差，為便于表示，用i表示探測(cè)站的序號(hào)，用j表示該探測(cè)站監(jiān)測(cè)到的地閃序號(hào)，那么有：

其中：θij表示第探測(cè)器i監(jiān)測(cè)到的第j個(gè)地閃的方位角；αij表示探測(cè)器（DF）i與第j個(gè)地閃真實(shí)的方位角，eij為探測(cè)器（DF）i探測(cè)地閃j產(chǎn)生的隨機(jī)誤差。

記地閃j的經(jīng)緯度坐標(biāo)為（Lj，φj）=（經(jīng)度，緯度），探測(cè)器（DF）i的坐標(biāo)（Ldi，φdi）=（經(jīng)度，緯度），那么αij可表示為：

最終Orville[5]推導(dǎo)出對(duì)應(yīng)第j個(gè)地閃的優(yōu)化位置P（Lj，φj）為：

2.2 考慮場(chǎng)地誤差下的優(yōu)化模式

考慮到場(chǎng)地誤差的情況，探測(cè)角與實(shí)際方位存在如下關(guān)系：

其中：Nf表示地閃樣本數(shù)，表示場(chǎng)地誤差存在情況下探測(cè)器i監(jiān)測(cè)到的第j個(gè)地閃的測(cè)量角，θij表示經(jīng)場(chǎng)地誤差訂正后的角度，為探測(cè)器i監(jiān)測(cè)到的第j個(gè)地閃的場(chǎng)地誤差。

我們假定系統(tǒng)具有下列形式的場(chǎng)地誤差：

Nh為傅里葉基數(shù)項(xiàng)數(shù)。這樣，對(duì)于某一個(gè)的值，由（5）-（9）可知第j個(gè)閃電最小的的值為，于是對(duì)所有閃電樣本可建立目標(biāo)函數(shù)?。

3 場(chǎng)地優(yōu)化模式應(yīng)用

3.1 定位資料來源

基礎(chǔ)數(shù)據(jù)源于內(nèi)蒙古自治區(qū)氣象局ADTD閃電定位系統(tǒng)2011-2014年的閃電定位觀測(cè)資料。通過對(duì)原始數(shù)據(jù)反復(fù)篩選，確定包頭青山、東勝、和林格爾3個(gè)訂正目標(biāo)站點(diǎn)，并篩選出大量同時(shí)性較好的回?fù)魯?shù)據(jù)，在篩選中還發(fā)現(xiàn)2013年8月11日閃電定位數(shù)據(jù)多達(dá)60688條，剔除掉無用數(shù)據(jù)依舊包含大量有效數(shù)據(jù)，適合用于雷達(dá)回波的對(duì)比驗(yàn)證。

根據(jù)歷年研究資料表明，負(fù)地閃發(fā)生的區(qū)域在回波強(qiáng)度大于35 dbz，云頂高度7 km以上，但由于雷達(dá)資料的原因，本文選取了呼和浩特雷達(dá)站回波強(qiáng)度大于35 dbz，云頂高度超過8 km的雷達(dá)資料用于對(duì)比驗(yàn)證。

3.2 閃電定位資料場(chǎng)地誤差的優(yōu)化

通過對(duì)數(shù)據(jù)的篩選和計(jì)算，最終得出包頭青山、東勝、和林格爾三個(gè)站點(diǎn)的場(chǎng)地誤差曲線如下：

圖1 各站點(diǎn)誤差曲線

其中：橫坐標(biāo)為實(shí)測(cè)方位角，縱坐標(biāo)為場(chǎng)地誤差。坐標(biāo)0點(diǎn)從上到下依次對(duì)應(yīng)包頭青山、東勝、和林格爾。

從圖1可以看出，各站點(diǎn)的場(chǎng)地誤差均在10°以內(nèi)。利用圖1的結(jié)果和公式（8）就可以對(duì)每一個(gè)探測(cè)站閃電定位資料場(chǎng)地誤差進(jìn)行優(yōu)化，然后重新交匯就可以求得經(jīng)過訂正后的閃電定位資料。[9]

3.3 雷達(dá)資料對(duì)比驗(yàn)證

由于雷達(dá)回波資料為6分鐘一組，通過閃電定位資料的分析發(fā)現(xiàn)，2013年8月11日00：00—00：08、2013年8月11日5：40—5：48兩個(gè)時(shí)間段的閃電定位樣本數(shù)據(jù)量較大，位置也相對(duì)集中，與雷達(dá)數(shù)據(jù)對(duì)比驗(yàn)證效果比較明顯。[10]因此，確定2013年8月11日00：00—00：08、2013年8月11日5：40—5：48兩個(gè)時(shí)間段的雷達(dá)資料進(jìn)行對(duì)比驗(yàn)證，效果如圖2、圖3：

圖2

圖3

圖2為實(shí)測(cè)數(shù)據(jù)與雷達(dá)回波疊加效果，圖3為經(jīng)場(chǎng)地優(yōu)化模式后數(shù)據(jù)與雷達(dá)回波效果疊加，由圖2與圖3比可以看出，部分原本在雷達(dá)強(qiáng)回波區(qū)域外的探測(cè)點(diǎn)返回了強(qiáng)回波區(qū)域內(nèi)，明顯更接近于實(shí)際情況。對(duì)站點(diǎn)閃電定位資料場(chǎng)地誤差優(yōu)化效果比較明顯。

3.4 雷災(zāi)資料驗(yàn)證

2013年9月9日13：40左右，呼和浩特市海拉爾大街華泰4S店院遭受雷擊，造成部分設(shè)備損壞，直接經(jīng)濟(jì)損失6萬余元。有人目擊稱當(dāng)時(shí)看到院內(nèi)桿塔遭受直接雷擊，根據(jù)現(xiàn)場(chǎng)勘查及相關(guān)數(shù)據(jù)，桿塔接地線實(shí)測(cè)剩磁為2.0T，并在斷接卡處有熔痕，說明桿塔確實(shí)有大電流流過。根據(jù)當(dāng)日閃電定位資料顯示，當(dāng)日在13：30—13：50分在雷擊點(diǎn)附近共監(jiān)測(cè)到三次地閃，坐標(biāo)如表1。

表1 監(jiān)測(cè)到雷擊點(diǎn)附近閃電坐標(biāo)

其中最接近雷擊點(diǎn)的228號(hào)地閃距離雷擊點(diǎn)的距離近2公里，229號(hào)地閃距離雷擊點(diǎn)的距離近4.7公里。在228號(hào)地閃數(shù)據(jù)中包含了本次優(yōu)化站點(diǎn)和林格爾和包頭青山兩個(gè)站點(diǎn)的回?fù)魯?shù)據(jù)，229號(hào)地閃數(shù)據(jù)中包含了和林格爾、包頭青山和東勝三個(gè)站點(diǎn)的回?fù)魯?shù)據(jù)，由于229號(hào)數(shù)據(jù)中和林格爾的回?fù)魯?shù)據(jù)探測(cè)時(shí)間上存在明顯誤差，應(yīng)予以剔除.因此，我們對(duì)228號(hào)數(shù)據(jù)和229號(hào)數(shù)據(jù)中優(yōu)化站點(diǎn)的回?fù)魯?shù)據(jù)利用誤差曲線優(yōu)化后重新交匯，得出新的閃電定位數(shù)據(jù)為：228號(hào)數(shù)據(jù)（經(jīng)度：40.8297，緯度111.613）、229號(hào)數(shù)據(jù)（40.8114，111.6013）。進(jìn)過優(yōu)化后的228號(hào)地閃距離雷擊點(diǎn)的距離為1.1公里，229號(hào)地閃距離雷擊點(diǎn)的距離近3公里，更接近與實(shí)際情況。但由于目前可以用于對(duì)比驗(yàn)證的直擊雷在案例僅有一例，因此雷災(zāi)案例的對(duì)比驗(yàn)證結(jié)果可能存在一定的偶然性，僅能用于參考。

4 結(jié)論與思考

本文主要利用了2011—2014年呼包鄂地區(qū)閃電定位資料，經(jīng)過對(duì)數(shù)據(jù)的反復(fù)篩查保留同時(shí)性較好的數(shù)據(jù)后建立數(shù)據(jù)樣本，確立包頭青山、東勝、和林格爾三個(gè)目標(biāo)站點(diǎn)，通過場(chǎng)地誤差優(yōu)化模式對(duì)數(shù)據(jù)樣本進(jìn)行優(yōu)化，繪制出包頭青山、東勝、和林格爾三個(gè)站點(diǎn)的場(chǎng)地誤差曲線，進(jìn)過誤差曲線對(duì)閃電定位資料進(jìn)行場(chǎng)地誤差訂正后，重新進(jìn)行交匯并與選取的雷達(dá)云頂高度和回波強(qiáng)度進(jìn)行對(duì)比驗(yàn)證。驗(yàn)證表明，該模式在地形復(fù)雜呼包鄂地區(qū)同樣具有一定的適用性。但在驗(yàn)證試驗(yàn)中發(fā)現(xiàn)，當(dāng)天氣過程較弱，閃電分布較分散的個(gè)別情況下，優(yōu)化后的定位資料依舊無法返回雷達(dá)強(qiáng)回波區(qū)，原因一直沒有找到，希望在日后的工作中能夠進(jìn)一步的完善和探討。

[1]許小峰．國(guó)家閃電監(jiān)測(cè)網(wǎng)的建設(shè)與技術(shù)分析[J]．中國(guó)工程科學(xué)，2002（5）：7-13．

[2]張義軍，孟青，馬明，等．閃電探測(cè)技術(shù)發(fā)展和資料應(yīng)用[J]．應(yīng)用氣象學(xué)報(bào)，2006，17（5）：611-620．

[3]張其林，劉曉東，劉明遠(yuǎn)，王靜．甘肅平?jīng)鲩W電定位資料誤差分析及其優(yōu)化[J]．氣象科技，2011，39（6）：808-813．

[4]陳明理，劉欣生，郭昌明，等．確定閃電定位系統(tǒng)場(chǎng)地誤差的參數(shù)化方法[J]．高原氣象，1990，9（3）：207-319．

[5]Orville R E．An analytical solution to obtain the optimum source location usingmultiple direction finders on a spherical suface [J]．J Geophys Res，1987，92（9）；10877-10866．

[6]陳明理，劉欣生，郭昌明．北京地區(qū)閃電定位系統(tǒng)場(chǎng)地誤差及其結(jié)構(gòu)分析[J]．氣象學(xué)報(bào)，1933，51（1）：66-74．

[7]Schutte T．An experimental stydy of the angle correction of the direction findersof the LLP system in Sweden[C]//Institute of High Tension Research，Uppsala，Sweden，UURIE，1984：158-184．

[8]馮民學(xué)，韋海容，焦圣明，等．南京市閃電定位資料的對(duì)比分析[J]．南京氣象學(xué)院學(xué)報(bào)，2008，31（2）：151-157．

[9]田芳，肖穩(wěn)安，馮民學(xué)，等．閃電定位儀觀測(cè)結(jié)果的修訂分析[J]．華東電力，2008，36（6）：38-42．

[10]徐洪澤．江蘇閃電定位網(wǎng)資料的分析與研究[D]．南京：南京信息工程大學(xué)，2007．

Error analysis of lightning location data and optimization of site error in the area of

Wang Hankun1Qin Zhaojun1Jin Shi2

1. Inner Mongolia Autonomous Region lightning protection center，Inner Mongolia Hohhot 010051 2. Hohhot Meteorological Bureau，Inner MongoliaHohhot 010051

The analysis of the sources of error of magnetic orientation method are used，based on the existing hu-bao-e area 2010-2014 lightning location data screening，analysis to determine the three stations in Baotou Dongsheng，Castle Peak，Helingeer，established by Chen Mingli and sample analysis，the eigenvalue of Orville technology development site error optimization model optimization based on the comparison with the location information corresponding to the radar data and lightning disaster data and revised. Results show that after site error of optimized positioning data originally in the strong echo area outside the back of strong echo area，significantly more in line with the actual situation，site optimization model can significantly improve the accuracy of lightning location data，have a positive significance for the future of the lightning monitoring and early warning，disaster investigation work.

lightning location；error；optimization

P427.32；P412

A

1009-6434（2016）11-0153-04