基于實測ADS-B數(shù)據(jù)的地面站覆蓋性分析

白宇駿

（中國民用航空華東地區(qū)空中交通管理局，上海 200335）

基于實測ADS-B數(shù)據(jù)的地面站覆蓋性分析

白宇駿

（中國民用航空華東地區(qū)空中交通管理局，上海 200335）

根據(jù)ADS-B地面站覆蓋范圍的需要，考慮地形遮蔽、地球曲率、頂空盲區(qū)及設(shè)備性能等因素，計算了ADSB地面站的理論覆蓋范圍?；贏DS-B數(shù)據(jù)，采用Lambert投影方法，提出了一種凹包生成算法；依據(jù)實測數(shù)據(jù)計算不同高度層的ADS-B覆蓋圖，提出一種ADS-B連續(xù)性計算方法，并計算漏點率和連續(xù)性等級。結(jié)果表明，該方法能夠提供地形遮蔽影響下的理論最大覆蓋范圍，可用于輔助ADS-B地面站的選址決策。

ADS-B；覆蓋范圍；Lambert投影方法；凹包生成算法

廣播式自動相關(guān)監(jiān)視技術(shù)（ADS-B，automatic dependent surveillance-broadcast）是一種能夠自動對外界廣播飛行信息的新一代監(jiān)視系統(tǒng)。相比傳統(tǒng)的監(jiān)視設(shè)備，ADS-B設(shè)備重量輕、安裝方便、維護費用低、電磁污染少、更新速率快，不僅可以使空中交通的安全等級得到提高，運行效率得到優(yōu)化，空域容量得到增加，還能夠使監(jiān)視空域范圍得以擴大。作為安全、經(jīng)濟、高效的新一代監(jiān)視設(shè)備，有極大的發(fā)展前景。中國民用航空ADS-B實施規(guī)劃提出，2016—2020年要實現(xiàn)中國范圍內(nèi)全空域ADS-B OUT的全面運行。并計劃在這5年內(nèi)完成386個ADS-B地面站新建工程，實現(xiàn)航線、航路及運輸機場的監(jiān)視覆蓋。而ADS-B地面站的覆蓋范圍需考慮地形遮蔽、地球曲率、頂空盲區(qū)及設(shè)備性能等一系列因素，以實現(xiàn)最優(yōu)覆蓋效果。因此，為實現(xiàn)良好的空域覆蓋并滿足監(jiān)視性能的相關(guān)標(biāo)準(zhǔn)，對ADS-B地面站覆蓋性的研究分析十分必要。

國內(nèi)外學(xué)者對ADS-B地面站的理論覆蓋范圍進行了研究。焦衛(wèi)東等[1]通過設(shè)計ADS-B規(guī)范標(biāo)準(zhǔn)的天氣模型，獲得了天線空間輻射場的增益數(shù)據(jù)，在此基礎(chǔ)上，根據(jù)信號自由空間傳輸距離、視距傳輸距離、遮蔽角導(dǎo)致的視距截止距離的計算結(jié)果，并結(jié)合地理高程數(shù)據(jù)，得出基站的信號覆蓋情況，表明地形對ADSB基站信號的傳播影響非常明顯。劉文評等[2]提出一種針對不規(guī)則地形的ADS-B地面站信息覆蓋分析方法，并通過實驗表明，不規(guī)則地形是影響地面站覆蓋的關(guān)鍵因素，Longley-Rice模型可以準(zhǔn)確模擬信號衰減。沈笑云等[3]提出一種針對開闊區(qū)域ADS-B地面站信號覆蓋分析方法，并做了相關(guān)實驗，結(jié)果表明，開闊區(qū)域地表反射是影響地面站覆蓋的關(guān)鍵因素，二徑模型可以準(zhǔn)確模擬信號衰減。駱敏華等[4]借鑒國內(nèi)外相關(guān)的ADS-B測試評估方案，對上海浦東國際機場及其周邊各個航路的不同高度層和特定區(qū)域的ADS-B信號進行了大量測試，根據(jù)實際接收的信號質(zhì)量綜合決策現(xiàn)場有效覆蓋區(qū)初步評估結(jié)果。Kostic等[5]介紹了虛擬三維地理信息系統(tǒng)（GIS）中的雷達覆蓋預(yù)測與可視化。Boci等[6]概述了用于設(shè)計和構(gòu)建無線電站點布局的射頻（RF）覆蓋分析方法，以支持美國聯(lián)邦航空管理局（FAA）的ADS-B空中交通現(xiàn)代化計劃。

基于ADS-B數(shù)據(jù)，采用Lambert投影方法，將采集的ADS-B數(shù)據(jù)中的航跡點等信息投影到以ADS-B地面站為原點的平面坐標(biāo)系中，利用凹包生成算法，依據(jù)給定高度層上的實測ADS-B航跡計算不同高度層的ADS-B覆蓋圖。在給定高度層下，通過ADS-B連續(xù)性計算，依據(jù)漏點率結(jié)果確定連續(xù)性等級，得到不同連續(xù)性級別的ADS-B信號在航段上的分布覆蓋結(jié)果。

1 ADS-B地面站理論覆蓋范圍計算

ADS-B地面站的覆蓋半徑不僅受到機載ADS-B發(fā)射機信號最大輻射距離的限制，還受地球曲率的影響。ADS-B信號以空間波形式傳播，即視距傳播：電波沿視線直接傳播到接收點。因此在建立ADS-B地面站覆蓋模型時有必要考慮地形遮蔽對ADS-B信號的影響?；贏DS-B通訊鏈的傳播特點，當(dāng)有障礙物（如山峰、丘陵等）存在于電波的傳輸路徑上時，機載ADS-B發(fā)射器所發(fā)出的信號就無法被地面站接收到，形成一個信號盲區(qū)。如圖1所示，沿ADS-B地面站接收天線與障礙物頂端做一條射線，該射線與水平面夾角為θ2，這時在山峰后面會形成一個盲區(qū)，當(dāng)飛機的下滑角小于θ2時，其機載ADS-B發(fā)射器所發(fā)出的信號無法被地面站接收到。

圖1 仰角對覆蓋盲區(qū)的影響Fig.1 Influence of elevation angle on coverage blind zone

遮蔽角是指障礙物頂點和地面站天線之間的連線與地平線之間的夾角，即

其中：ht為地表某點的海拔高度；ds為地表某點與地面站的斜距；hr為地面站天線高度；Re為等效地球半徑（這里取Re=8 496 km）。將{θi}的最大值作為基準(zhǔn)遮蔽角θ0。實際應(yīng)用中需對遮蔽角加入修正因子，即

其中：λ為ADS-B系統(tǒng)工作波長；dsmax為基準(zhǔn)遮蔽角對應(yīng)的斜距。將遮蔽角的相應(yīng)DEM網(wǎng)格點可視信息儲存起來，作為視點與目標(biāo)點是否通視的判定依據(jù)之一。

所得修正遮蔽角可用于計算不同方位ADS-B的視線截止距離，即

進一步計算該遮蔽角下某高度層H0上的ADS-B信號的視線截止距離，如圖2所示，即

圖2 遮蔽角計算模型Fig.2 Calculation model of defilade angle

根據(jù)視線截止距離所形成的一系列視線截止點，將所有方向上的視線截止點擬合成一條封閉曲線，則該曲線所形成的區(qū)域即為ADS-B的實際覆蓋范圍。

2 ADS-B地面站實際覆蓋范圍計算

通過采集不同高度的ADS-B航跡點，并將航路點采用Lambert投影方式投影至平面，計算投影平面上航路點的最小凹包，由此得到ADS-B地面站在不同高度層上的覆蓋范圍。

2.1 ADS-B航跡投影方法

如圖3所示，Lambert投影是一種等角圓錐投影，設(shè)想用一個正圓錐割于球面兩標(biāo)準(zhǔn)緯線，應(yīng)用等角條件將地球面投影到圓錐面上，然后沿一母線展開，即為Lambert投影平面。Lambert等角投影后緯線為同心圓弧，經(jīng)線為同心圓半徑。

圖3 Lambert投影Fig.3 Lambert projection

以圖幅的原點經(jīng)線L0（一般是圖幅中央經(jīng)線）作縱坐標(biāo)x軸（向北為正），原點經(jīng)線L0與原點緯線B0（一般是圖幅最南端緯線）的交點作為原點O，過此點的切線作為橫坐標(biāo)y軸（向東為正），構(gòu)成Lambert平面直角坐標(biāo)系，如圖4所示。

圖4 平面直角系坐標(biāo)Fig.4 Lambert coordinate system

將地理坐標(biāo)（Bd，Ld）轉(zhuǎn)化為Lambert投影下的平面坐標(biāo)（xd，yd），Lambert等角投影正解公式為

縱坐標(biāo)為

橫坐標(biāo)為

其中：a為地球長半軸；e為第一偏心率；B1為第一標(biāo)準(zhǔn)緯線緯度；B2為第二標(biāo)準(zhǔn)緯線緯度；r0為原點維度處的 r值；rd為 Bd緯度處的 r值；mB1、mB2分別為標(biāo)準(zhǔn)緯線B1和 B2處的m 值；tB1、tB2分別為標(biāo)準(zhǔn)緯線 B1和 B2處的t值。

2.2 ADS-B航跡覆蓋范圍凹包

依據(jù)給定高度層上的實測ADS-B航跡計算不同高度層的ADS-B覆蓋圖。即基于步驟2中得到的給定高度層上所有ADS-B航跡點在平面上的投影坐標(biāo)，計算信號覆蓋凹包，如圖5所示。

凹包計算步驟如下：

步驟1首先以ADS-B地面站為原點O（xd，yd），將水平面分為n個扇區(qū)，因此每個扇區(qū)的夾角為2π/n。

圖5 凹包生成算法示意圖Fig.5 Diagram of concave hull algorithm

步驟2對于每個ADS-B航跡點Pi在平面上的投影坐標(biāo)（xd，yd），計算向量 OPi與 x軸的夾角，即

步驟3計算第k個扇區(qū)中包含的航跡點集合Sk，對于每個ADS-B航跡點Pi，判斷向量OPi是否位于扇區(qū) k，判斷依據(jù)為（k+1）2π/n ＞ αi≥ k2π/ni。

步驟4若第k個扇區(qū)中包含的航跡點集合Sk不為空，則計算點集Sk中距離ADS-B地面站最遠的航跡點 Pk，即任意反之，若Sk為空，則令該扇區(qū)離ADS-B地面站最遠的航跡點為P0。

步驟5將每個扇區(qū)離ADS-B地面站最遠的航跡點構(gòu)成多邊形，所得凹多邊形即為覆蓋所有航跡點的凹包。

3 ADS-B信號連續(xù)性在航段上的分布

3.1 ADS-B信號連續(xù)性的計算

在正常情況下為了保證縮小航空器間隔和運行的安全性，ADS-B的更新時間間隔應(yīng)為0.5～2 s。但無論如何，在計算正確的情況下，每1 s都應(yīng)該有位置報文信息的存在?；谶@樣的原理可以判定，只要每1 s有位置報文，則認為報文不存在漏點，假若某1 s沒有位置報文，則判定報文數(shù)據(jù)丟失，認為報文存在1次漏點，連續(xù)多秒不存在報文數(shù)據(jù)，則判定ADS-B多個報文數(shù)據(jù)丟失，判定為多個漏點。

首先將ADS-B的漏點數(shù)初始化為0，然后在一段時間間隔內(nèi)提取某架飛機的ADS-B報文數(shù)據(jù)，如果連續(xù)m個系統(tǒng)周期內(nèi)未收到航跡點，則認為ADSB探測出現(xiàn)漏點，并為ADS-B漏點數(shù)加m，否則繼續(xù)提取數(shù)據(jù)，直到時間間隔結(jié)束。最后根據(jù)情況有

3.2 ADS-B信號連續(xù)性在航段上的分布

給定高度層下，通過ADS-B信號連續(xù)性計算，依據(jù)漏點率結(jié)果確定連續(xù)性等級，得到不同連續(xù)性級別的ADS-B信號在航段上的分布覆蓋結(jié)果。按民航劃分的飛行高度層，將不同高度層上航班的ADS-B實測數(shù)據(jù)進行分類，并根據(jù)距離誤差和角度誤差準(zhǔn)則將ADS-B實測數(shù)據(jù)匹配到對應(yīng)的航路。某航空器的ADS-B實測數(shù)據(jù)輸出位置為Pi（xi，yi），直線航段兩端的坐標(biāo)分別為 A（xa，ya）、B（xb，yb）?？汕蟮?ADS-B 實測數(shù)據(jù)到航段的距離為

假設(shè)ADS-B實測數(shù)據(jù)的航跡角為φi，則ADS-B實測數(shù)據(jù)與航段的夾角為

當(dāng)ADS-B實測數(shù)據(jù)到航段的距離di≤δ，且ADS-B實測數(shù)據(jù)與航段的夾角Δθi≤α?xí)r，該數(shù)據(jù)匹配到該航段，否則不能匹配到該航段，其中δ與α為給定的閾值。

對上述分段航路，統(tǒng)計在一段時間間隔內(nèi)，該航段上所有航空器出現(xiàn)的漏點數(shù)，依據(jù)ADS-B信號連續(xù)性計算原理，通過漏點率結(jié)果確定連續(xù)性等級。選定任意一架航空器，連續(xù)m秒內(nèi)未收到航跡點，則判斷為ADS-B的航跡漏點，上述航段上所有航空器出現(xiàn)的漏點數(shù)，根據(jù)漏點率計算公式得到丟失點百分比計算結(jié)果。依據(jù)事先設(shè)置的連續(xù)性等級辨別表，判斷不同航段上漏點率所處的范圍，確定該航段的連續(xù)性等級，得到給定高度層下不同連續(xù)性級別的航段在航路上的分布，如表1所示。

表1 連續(xù)性等級判別表Tab.1 Continuity grade

4 案例研究

為驗證本文所提方法的正確性，以南京某ADS-B地面站獲取的數(shù)據(jù)為例進行分析，該ADS-B地面站天線所處位置為：東經(jīng)118.786 111°、北緯31.939 167°、高度50 m。采集2017年3月22日上午9：00—12：00的航跡數(shù)據(jù)，在NASA World Wind地理信息系統(tǒng)平臺上開發(fā)了ADS-B地面站覆蓋性分析軟件，結(jié)果如圖6～圖8所示。

圖6 ADS-B地面站5 400 m高度理論覆蓋范圍Fig.6 Theoretical coverage of ADS-B station atheight of 5 400 m

圖7 ADS-B地面站5 400 m高度實際覆蓋范圍凹包及ADS-B航跡連續(xù)性分布圖Fig.7 Concave hull of actual coverage of ADS-B station at height of 5400 m and ADS-B track continuity

圖8 ADS-B地面站5 400 m高度理論覆蓋范圍與實際覆蓋范圍凹包疊加結(jié)果Fig.8 Concave hull overlap of theoretical coverage and actual coverage of ADS-B station at height of 5 400 m

圖6所示為ADS-B地面站5 400 m高度理論覆蓋范圍，其中圓形區(qū)域為僅考慮地球曲率影響下的理論最大覆蓋范圍，而多邊形區(qū)域為考慮地形遮蔽影響下的理論最大覆蓋范圍。該結(jié)果可用于輔助ADS-B地面站的選址決策。

圖7所示為將空域分為60個扇區(qū)后ADS-B地面站5 400 m高度實際覆蓋范圍凹包，該凹包可用于ADS-B運行階段，為判斷航空器是否即將飛出ADS-B地面站覆蓋范圍提供依據(jù)。

圖8所示為ADS-B地面站5 400 m高度考慮地形遮蔽因素時理論覆蓋范圍與實際覆蓋范圍凹包疊加結(jié)果。從圖中可以看出，ADS-B實際數(shù)據(jù)的最遠距離驗證了ADS-B理論覆蓋模型的結(jié)果，但由于實際ADS-B數(shù)據(jù)較少，二者覆蓋的區(qū)域差別較大，需要進一步增加實測數(shù)據(jù)來提高實際覆蓋范圍的精確性。

5 結(jié)語

考慮地形遮蔽、地球曲率及設(shè)備性能對ADS-B地面站覆蓋范圍的影響，以南京某ADS-B地面站為例計算了ADS-B地面站的理論覆蓋范圍，測得不同高度層下ADS-B的實測覆蓋圖，計算并匹配得到不同連續(xù)性級別的ADS-B信號在航段上的分布覆蓋結(jié)果。該研究可為ADS-B地面站的選址和運行輔助決策提供有效的依據(jù)和支持。

[1]焦衛(wèi)東,周 波,沈笑云.ADS-B地面站信號覆蓋分析及選址應(yīng)用[J].計算機仿真,2014,31(4):51-55.

[2]劉文評,曹 博,劉志剛,等.不規(guī)則地形對ADS-B地面站信號覆蓋影響分析[J].電光與控制,2016(6):84-89.

[3]沈笑云,曹 博,張思遠,等.開闊區(qū)域ADS-B地面站信號覆蓋仿真分析[J].計算機仿真,2015,32(3):94-99.

[4]駱敏華,周祿華,毛 建.浦東ADS-B覆蓋測試與分析方法研究[J].網(wǎng)絡(luò)安全技術(shù)與應(yīng)用,2015(6):115-116.

[5]KOSTIC A,RANCIC D．Radar Coverage Analysis in Virtual GIS Environment[C]//TELSIKS 6th International Conference on Telecommunications in Modern Satellite,Cable and Broadcasting Service,IEEE,Serbia and Montenegro,2003:721-724．

[6]BOCI E,SARKANI S,MAZZUCHI T．Optimizing ADS-B RF Coverage[C]//Integrated Communications,Navigation&Surveillance Conference,Arlington,VA,USA,2009:1-10.

（責(zé)任編輯：劉佩佩）

Coverage analysis on ADS-B station based on measured data

BAI Yujun
（East China Regional Air Traffic Management Bureau of Civil Aviation of China,Shanghai 200335,China）

According to the factors such as topography,earth curvature,headspace blind zone and equipment performance,theoretical coverage of ADS-B ground station is calculated.Based on ADS-B data,Lambert projection method is used to project the generation algorithm of concave hull.Based on measured data,the coverage of ADS-B under different layers of height is calculated.An ADS-B continuity calculation method is proposed to calculate the leakage rate and continuity level.Results show that the proposed method can provide theoretical maximum coverage under the influence of terrain shelter and can be used to assist ADS-B ground station’s location decision.

ADS-B;coverage;Lambert projection method;generation algorithm of concave hull

V243

：A

：1674－5590（2017）04－0011－05

2017-04-15；

：2017-05-23

：中國民用航空局空中交通管理局科技項目（KJ1506）

白宇駿（1975—），男，廣東清遠人，工程師，工學(xué)博士，研究方向為空管通信導(dǎo)航監(jiān)視設(shè)備規(guī)劃與運行管理.