一種基于ADS-B信號的二次雷達精度評估方法

姚 淵

（民航中南空管局技術保障中心，廣東 廣州510403）

0 引言

當前對二次雷達精度評估方法大體分三類：第一類為官方評估工具，廣泛使用的歐控的SASS-C和FAA的RBAT；第二類為商業公司評估工具，權威工具是Intersoft-Electronics公司開發涵蓋雷達系統所有軟硬件分析的RASS系列工具；第三類為各雷達廠家專用評估工具，比如Leonardo公司的ARTESAES，該類工具可對其系統進行更深入的分析。上述工具都是對二次雷達數據重構得出基準數據進行相應評估。本文介紹的評估方法是利用ADS-B數據重構基準數據，得出雷達相應精度評估值。

1 監視數據精度

1.1 二次雷達和精度

二次雷達通過向空中輻射詢問信號，接收應答機應答信號來進行目標定位和監視。提供數據包括目標識別（Mode A）、高度（Mode C）、斜距（Range）和方位（Azimuth）。實際精度評估主要關注距離和方位精度。精度的系統誤差主要來源于時鐘固有誤差、正北校準誤差等，隨機誤差主要來源于方位抖動、時鐘抖動等。民航局二次雷達規范[1]中要求精度如表1所示。

表1 中國民航局標準

1.2 ADS-B和精度

ADS-B通過航空器廣播位置、速度、識別信息等數據，地面站接收處理這些信息實現目標的被動探測、定位。其位置精度取決于機載設備提供的位置精度，包括數據源精度和報告精度。位置數據源通常包括：（1）慣導系統，根據傳感器數據更新位置，精度隨時間降級，通常沒有數據源完整性報告；（2）GPS系統，提供位置數據的同時會給出GPS精度和GPS完整性報告。GPS精度由接收機報告的HFOM表示，DO260A和DO260B應答機使用該值報告導航精度類別NAC；GPS完整性由接收機報告的HPL表示，DO260A和DO260B應答機使用該值來報告導航完整性類別NIC，DO260應答機使用該值來報告導航不確定性類別NUC。DO260應答機沒有給出GPS精度報告，但是Boeing公司研究表明其精度可由NUC值推算得出：Accuracy＜HPL/4。EASA在AMC20-24[2]中給出NUC和NACp之間對應關系如表2所示。

表2 NUC和NACp之間的相互轉換關系

2 精度評估的方法

精度評估過程包含5個步驟：（1）目標錄像和解碼；（2）目標坐標轉換；（3）目標篩選；（4）目標重建；（5）實施評估。評估軟件基于微軟.NET3.5框架，由C#語言實現。錄像部分采用了SharpPcap開源程序包（https：//sourceforge.net/projects/sharppcap/）。

2.1 目標錄像和解碼

二次雷達和ADS-B數據使用歐控ASTERIX數據格式，程序通過WinPcap組件構建標準Pcap錄像。解碼參考標準ASTERIX格式，二次雷達使用CAT034/CAT048包，ADS-B數據使用CAT021包。

2.2 目標坐標轉換

二次雷達和ADS-B數據基于不同坐標系，二次雷達為本地球極坐標，ADS-B為WGS84地表笛卡爾坐標 TA（Lat，Lng，Ha），因此進行精度評估時需進行坐標轉換。見圖1。

圖1 雷達球極坐標和地表笛卡爾坐標

WGS-84框架包含一個基準橢球體、一個標準坐標系以及一個標準大地水平面。其中，基準橢球體中心為地球質心。Z軸指向國際時間服務機構BIH定義[3]的CTP方向；X軸指基準子午線（0°）與CTP定義赤道面的交點；Y軸指向赤道面內右手定理確定的X軸東面90°。基準橢球體參數如表3所示。

表3 WGS-84橢球體參數

坐標轉換需要三個步驟：（1）地表笛卡爾坐標到地心笛卡爾坐標；（2）地心笛卡爾坐標到本地笛卡爾坐標；（3）本地笛卡爾坐標到本地球極坐標。

假設ADS-B目標地表笛卡爾坐標表示為TA（Lat，Lng，Ha），將其轉換為地心笛卡爾坐標為 TAC（Xc，Yc，Zc）。轉換公式為：

由TAC（Xc，Yc，Zc）轉換為本地笛卡爾坐標TRL（Xl，Yl，Zl）。轉換公式為：

其中地表笛卡爾坐標為 R（L，G，h）。Sl為地心笛卡爾坐標到地表笛卡爾坐標的轉換矩陣。Tl為雷達地表笛卡爾坐標到地心笛卡爾坐標的轉換。

由TRL（Xl，Yl，Zl）轉換為雷達本地球極坐標目標TRP（η，θ，ψ）。轉換公式為：

經過上述轉換，ADS-B目標 TA（Lat，Lng，Ha）轉換成了雷達目標TRP（η，θ，ψ）。

2.3 目標篩選

在目標解碼和坐標轉換后需要對用以評估的目標進行篩選，篩選過程如圖2所示。

圖2 目標篩選過程圖

（1）目標區域選擇：用于精度評估的數據在空間上必須是一致的，在目標區域選擇時，可使用門限TA（Rho，Theta，FL）進行篩選；

（2）目標質量選擇：用于精度評估的ADS-B目標必須是高質量的，在目標質量選擇時，可使用門限TQ（Function（NAC，NIC，NUCp））進行篩選；

（3）關聯航跡選擇：在質量篩選完成后，進行關聯目標的選擇。選擇優先級如下：航空器地址AA，航空器識別ID，識別碼（Mode A）。使用關聯表達式A（AA，AId，ModA）。

經過上述篩選和關聯后，對目標進行重建以及精度評估。

2.4 目標重建

因兩者數據更新率不一樣，評估前選擇插值法進行重建。同時民航航班速度小、機動少，航跡點密集，采用一次線性插值。一次線性插值描述：已知點P0（R0，T0，FL0）和P1（R1，T1，FL1）點，使用P0和P1構造直線方程求得 P（R，T，FL）。方程如下：

其中，R表示斜距Rho，T表示方位Theta，FL表示高度層。

數據插值步驟為：（1）選取相關維度。通過分析數據Rho值、Theta值和FL值取值范圍，選取一個維度，如分析Theta精度可選取Rho維度進行重建；（2）重建相關維度。根據評估所需，按一定步進Step進行重建。重建后目標將在此維度上完全一致；（3）構建其他維度。使用前文一次插值方法對其他維度進行重建。

2.5 實施評估

通過選取一致維度，計算其他維度重建數據差值均方根E實現評估。表達式如下：

其中，E為總體差值均方根，ei為重建數據第i個點的差值。

雷達數據單位為Rho-海里、Theta-度，實際誤差ei非常小，平方后會損失精度。因此實際評估過程中使用 Rho-（米）、Theta-0.01（度）。

3 精度評估實驗

評估實驗使用了中南某雷達站S模式二次雷達和ADS-B地面站1小時錄像。篩選參數Rho[0，256]、Theta[315，360]、FL[0，350]、T（NUCp）≥ 8，AA=00B1F0 相 關 ； 選 取 Rho[241.48，244.03]，Rho_Step=0.05海里重建，得出目標Theta緯度E值0.04度。如圖3所示。

圖3 距離重建結果方位偏差

選取Theta[344.64，348.27]，Theta_Step=0.07 度重建，得出目標Rho緯度E值15.22 m。如圖4所示：

圖4 方位重建結果距離偏差

4 總結

本文介紹了使用ADS-B作為基準數據對二次雷達進行精度評估的方法以及相關實驗。該方法可有效實現對二次雷達精度的評估。隨著新監視技術的發展和使用，可以納入更多監視數據源進行分析，也可進行各監視源的互相驗證分析等。