準天頂衛星 サブメータ級補強機能の性能評価 電子情報通信学会 宇宙・航行エレクトロニクス研究会 July 29, 2005 準天頂衛星 サブメータ級補強機能の性能評価 坂井 丈泰、福島 荘之介、新井 直樹、 伊藤 憲 電子航法研究所

Introduction 我が国は、2008年頃の衛星打上げを目指して準天頂衛星システム（QZSS）を計画中。 July 2005 Sakai, ENRI Introduction 我が国は、2008年頃の衛星打上げを目指して準天頂衛星システム（QZSS）を計画中。 ８の字軌道で高仰角から測位・放送サービスを提供。測位ミッションについては、基本的にGPSの補完あるいは補強を行う。 補完＝GPS同様の信号を追加 補強＝ディファレンシャル補正情報（＋インテグリティ情報） 各省庁・研究機関が研究開発を分担。電子航法研究所は、サブメータ級補強情報サービスの開発を担当。 現在、サブメータ級補強信号の内容について、SBAS（静止衛星によるGPS補強システム）をベースとして検討中。プロトタイプを作成し、性能評価を試みた。

準天頂衛星の軌道（１） 楕円軌道 低軌道衛星（LEO） NNSS（極軌道） 傾斜地球同期軌道（IGSO） ＝８の字軌道 中軌道衛星（MEO） July 2005 Sakai, ENRI 準天頂衛星の軌道（１） 楕円軌道 低軌道衛星（LEO） NNSS（極軌道） 傾斜地球同期軌道（IGSO） ＝８の字軌道 中軌道衛星（MEO） GPS, GLONASS, Galileo 静止軌道（GEO） 高度35,786km SBAS 地球 90 180 -60 60 Longitude, deg Latitude, deg 東経140度を中心に配置した場合の地上軌跡の例 離心率0、軌道傾斜角55度

準天頂衛星の軌道（２） 離心率を変えて、円→楕円軌道とした例。 軌道傾斜角55度、遠地点は北半球側。北半球ではゆっくり移動する。 July 2005 Sakai, ENRI 準天頂衛星の軌道（２） -180 -90 90 180 -60 60 Longitude, deg Latitude, deg e=0 e=0.2 e=0.4 e=0.6 長楕円軌道 ８の字軌道 離心率を変えて、円→楕円軌道とした例。 軌道傾斜角55度、遠地点は北半球側。北半球ではゆっくり移動する。 e=0.1～0.2程度が有望。

アベイラビリティ計算例（マスク20度、要求精度＝水平10m） July 2005 Sakai, ENRI アベイラビリティの向上 1 0.90 0.95 1 0.999 0.9995 125 130 135 140 145 150 25 30 35 40 45 50 Longitude, E Latitude, N 125 130 135 140 145 150 25 30 35 40 45 50 Longitude, E Latitude, N 平均 衛星数 5.3 平均 衛星数 9.3 GPS衛星（24機）のみ GPS＋MTSAT 2機＋準天頂3機 アベイラビリティ計算例（マスク20度、要求精度＝水平10m）

サブメータ級補強信号 我が国全域を対象としたディファレンシャル補正情報（広域ディファレンシャルGPS＝WADGPS）。 July 2005 Sakai, ENRI サブメータ級補強信号 我が国全域を対象としたディファレンシャル補正情報（広域ディファレンシャルGPS＝WADGPS）。 GPS L1周波数で放送。同一のアンテナで受信できる。 補強対象：GPS、準天頂衛星、（ガリレオ）。目標精度＝1m。 ベクトル補正方式：衛星軌道、クロック、電離層伝搬遅延をそれぞれ別々に補正する。 インテグリティ情報あり：信頼性の高い位置情報を提供。 すでに実用化されているSBAS（静止衛星による補強システム）をベースとして開発。 主要な課題：電離層伝搬遅延補正、インテグリティ情報

広域ディファレンシャル補正方式 衛星軌道情報の誤差 対流圏 電離層 電離層遅延（～100m） 対流圏遅延（～20m） 衛星クロック誤差 July 2005 Sakai, ENRI 広域ディファレンシャル補正方式 衛星軌道情報の誤差 対流圏 電離層 電離層遅延（～100m） 対流圏遅延（～20m） 衛星クロック誤差 ユーザ位置の関数ではない すべてのユーザに対して 同じ寄与 SA ONなら速い変動 寄与の程度はユーザ位置による （視線方向成分が問題） 変動の周期は数10分以上 ユーザ位置の関数 垂直構造は薄膜で近似など ユーザ位置（特に高度）の関数 モデルによる補正が有効

補強メッセージ検討 目的：日本全域に対して有効な広域補強情報をリアルタイムにユーザに伝送し得るメッセージ構成であること。 July 2005 Sakai, ENRI 補強メッセージ検討 目的：日本全域に対して有効な広域補強情報をリアルタイムにユーザに伝送し得るメッセージ構成であること。 制約：データレートは 250bps 以内。 すでに実用化されているSBAS（静止衛星による補強システム）メッセージをベースとする。 静止衛星依存部分は準天頂衛星向けに変更する。 電離層・対流圏伝搬遅延については、高精度化が可能な拡張メッセージを検討する。 実現可能な測位精度について予備検討が必要（量子化単位は12.5cmでよいか？）。 今回の課題：実際にSBASベースの補強メッセージで所要の性能が達成できるか、プロトタイプにより評価する。

補正情報の具体的内容（ＳＢＡＳ） 衛星クロック補正：エフェメリスから得た衛星クロックに対する補正量。 July 2005 Sakai, ENRI 補正情報の具体的内容（ＳＢＡＳ） 衛星クロック補正：エフェメリスから得た衛星クロックに対する補正量。 衛星軌道補正：エフェメリスから得た衛星のECEF直交座標値に対する補正量。Δx、Δy、Δzの３成分。 電離層遅延補正：経緯度で５度毎に設定された格子点（IGP）における垂直遅延量として放送。 ユーザ側ではあらかじめ決められた内挿法（双一次線形補間）により必要な位置（IPP）における遅延量を得る。 垂直構造は薄膜で近似。衛星の仰角により垂直→傾斜変換。 対流圏遅延補正：モデルにより遅延量を計算する。 衛星の仰角により垂直→傾斜変換する。 拡張メッセージについては今回は考慮しない。

補正情報の生成手順 モニタ局における測定データを収集 対流圏遅延補正 電離層遅延補正 クロック／軌道誤差を分離 前処理など July 2005 Sakai, ENRI 補正情報の生成手順 モニタ局における測定データを収集 対流圏遅延補正 電離層遅延補正 クロック／軌道誤差を分離 前処理など 電離層遅延補正データを作成 補正情報の品質チェック ユーザに向けて放送 2周波の観測データ 時刻あわせ 整合性チェック モデルで補正する ２周波データにより補正 データの整合性チェック 実際に補正情報を適用してみる 標準フォーマットで放送

クロック／軌道誤差分離 対流圏遅延・電離層遅延を除いた測距誤差は、衛星クロック誤差と軌道情報の誤差の線形結合。 July 2005 Sakai, ENRI クロック／軌道誤差分離 対流圏遅延・電離層遅延を除いた測距誤差は、衛星クロック誤差と軌道情報の誤差の線形結合。 クロックと軌道誤差の性質の違いを利用して、これらを分離： 衛星クロック誤差：すべてのユーザに一様な誤差となる。 軌道誤差：ユーザ位置によって影響が異なる。それぞれのユーザにとっての視線方向成分が誤差として現れる。 真の位置 エフェメリスによる位置 ユーザ２ ユーザ１ ユーザ１にとっての軌道誤差 ユーザ２にとっての軌道誤差 衛星位置の誤差

クロック／軌道誤差分離 マルチパス誤差などを避けるため、比較的長い時定数を持たせてカルマンフィルタで処理する。 July 2005 Sakai, ENRI クロック／軌道誤差分離 マルチパス誤差などを避けるため、比較的長い時定数を持たせてカルマンフィルタで処理する。 多数のモニタ局が幾何学的にまんべんなく分布しているほうが分離しやすい： ただし、ディファレンシャル補正を行うことが目的なので、必ずしも完全に分離する必要はない。 モニタ局はサービスエリア内でなるべく広い面積をカバーするように配置するとよい。 精密軌道暦がリアルタイムに取得できる場合は、それを利用してもよい。 衛星クロック誤差については別途推定したほうがよい。

電離層伝搬遅延補正 今回のプロトタイプでは、MSAS/WAASが採用している プレーナフィット（planar fit）方式を利用： July 2005 Sakai, ENRI 電離層伝搬遅延補正 今回のプロトタイプでは、MSAS/WAASが採用している プレーナフィット（planar fit）方式を利用： 平面モデル：IGP周辺の電離層垂直 遅延量が経緯度の一次関数になる ものとして、IGP位置における遅延量 を求める。 周波数間バイアス：電離層 遅延量の測定で問題となる 周波数間バイアス（IFB）は、 リアルタイムに推定・除去。 推定処理は1分毎に実行 （IFB推定は10分毎） 電離層遅延量 Rmax IGP IGPにおける垂直遅延量の推定値

プロトタイプの実行条件 ＰＣ上で動作するプロトタイプシステムを作成した。 July 2005 Sakai, ENRI プロトタイプの実行条件 ＰＣ上で動作するプロトタイプシステムを作成した。 準天頂システムの5モニタ局、あるいはMSASと同じ６モニタ局で、電離層活動の違いにより３つの時期について計算を実行。 時期I 静穏 2004年6月22～24日 時期II 磁気嵐 2004年7月22～24日 時期III 強い磁気嵐 2004年11月8～10日 利用するのは二周波の擬似距離のみ。搬送波位相は使わない。データソースは国土地理院GEONET（30秒サンプル）。 量子化はすべて 12.5cm 単位（SBASと同じ）。 電離層伝搬遅延補正：IGPは5度毎に配置（SBASと同じ） 対流圏伝搬遅延補正：SBASのモデルによる。 完全なSBASメッセージを毎秒生成、ファイルに記録。 NovAtelフォーマット（$FRMAレコード）を流用

July 2005 Sakai, ENRI 補強メッセージのサンプル $FRMA,272,86403.130,134,80811EA4,250,53081FFDFFDFFDFFDFFDFFDFFDFFDFFDFFDFFDFFDFFFBBBBBBBBBBBBAC1CD280*7C $FRMA,272,86404.130,134,80811EA4,250,9A0C1FFDFFDFFDFFDFFDFFDFFDFFDFFDFFDFFDFFDFFFBBBBBBBBBBBBB7E76F80*0F $FRMA,272,86405.130,134,80811EA4,250,C661FFDFFDFFDFFDFFDFFFBBBBB8800000000000000000000000000036CD8A40*70 $FRMA,272,86406.130,134,80811EA4,250,5306FFBFFFF8000000000000000000000000000000000000000000002B963FC0*0D $FRMA,272,86407.130,134,80811EA4,250,9A091FFDFFDFFDFFDFFDFFDFFDFFDFFDFFDFFDFFDFFFBBBBBBBBBBBB806D3340*77 $FRMA,272,86408.130,134,80811EA4,250,C60D1FFDFFDFFDFFDFFDFFDFFDFFDFFDFFDFFDFFDFFFBBBBBBBBBBBB924AAE40*08 $FRMA,272,86409.130,134,80811EA4,250,5361FFDFFDFFDFFDFFDFFFBBBBB89000000000000000000000000000021FE640*73 $FRMA,272,86410.130,134,80811EA4,250,9A61FFDFFDFFDFFDFFDFFFBBBBB8A00000000000000000000000000039994D00*05 $FRMA,272,86411.130,134,80811EA4,250,C60A1FFDFFDFFDFFDFFDFFDFFDFFDFFDFFDFFDFFDFFFBBBBBBBBBBBBA6BE8CC0*03 $FRMA,272,86412.130,134,80811EA4,250,530E1FFDFFDFFDFFDFFDFFDFFDFFDFFDFFDFFDFFDFFFBBBBBBBBBBBBA99E5040*0A メッセージタイプID （左６ビット） プリアンブル メッセージ長 衛星PRN番号 （WAAS-PORと同じ） 時刻 CRC

ユーザ受信機シミュレータ ユーザ受信機シミュレータも、ＰＣ上で動作する計算機プログラム。 July 2005 Sakai, ENRI ユーザ受信機シミュレータ ユーザ受信機シミュレータも、ＰＣ上で動作する計算機プログラム。 ユーザ局における観測データを入力として、測位計算を実行： SBAS補強メッセージも入力された場合、これを適用する。 測位結果のほか、プロテクションレベルも計算。 補強メッセージは、NovAtelフォーマット（$FRMAレコード）。SBASとの区別はまったくない。 一周波の擬似距離のみを利用。搬送波位相は使わないが、キャリアスムージングを適用。データソースは国土地理院GEONET（30秒サンプル）。 航法ユーザを想定し、出力フィルタはなし。

モニタ局およびユーザ局の配置 MSAS 6局 QZSS 5局 JAXAが設置するモニタ局の候補地 July 2005 Sakai, ENRI モニタ局およびユーザ局の配置 MSAS 6局 QZSS 5局 ● モニタ局 ● ユーザ局 （南鳥島→父島に変更） JAXAが設置するモニタ局の候補地 実際には、補強情報生成用にはさらにGEONETを利用できる MSASモニタ局の設置場所 さらにハワイとオーストラリアに標定局がある

ユーザ測位誤差（御前崎：水平） RMS誤差 1.842m (S/A) 0.423m (QZSS) 0.423m (MSAS) July 2005 Sakai, ENRI ユーザ測位誤差（御前崎：水平） RMS誤差 1.842m (S/A) 0.423m (QZSS) 0.423m (MSAS)

ユーザ測位誤差（御前崎：垂直） RMS誤差 2.958m (S/A) 0.833m (QZSS) 0.692m (MSAS) July 2005 Sakai, ENRI ユーザ測位誤差（御前崎：垂直） RMS誤差 2.958m (S/A) 0.833m (QZSS) 0.692m (MSAS)

July 2005 Sakai, ENRI 補正量のサンプル

July 2005 Sakai, ENRI 測位精度の比較 RMS値、単位 [m] 1.566 3.003 0.688 0.854 0.674 0.744 1.909 3.083 0.471 0.822 0.437 0.746 1.813 3.300 0.510 0.775 0.462 0.693 1.627 3.020 0.627 0.920 0.450 0.727 1.863 3.349 0.657 0.744 0.691 0.748 1.842 2.958 0.423 0.833 0.423 0.692 1.929 3.305 0.470 0.719 0.453 0.640 1.837 3.058 0.709 0.885 0.479 0.653 3.151 4.867 2.958 2.561 2.008 2.561 3.585 5.624 1.641 2.046 1.375 1.659 3.125 5.342 1.756 2.119 1.257 1.624 3.171 5.143 1.552 2.455 1.265 1.815 単独 QZSS 5局 MSAS 6局 時期I （静穏） 時期II （磁気嵐） 時期III （強い磁気嵐） 男鹿 水平　　垂直 御前崎 高山 高知 条件 （一周波、30秒サンプル、キャリアスムージングあり）

補正性能の傾向 強い磁気嵐が発生していない限り、QZSS 5局のモニタ局配置で、水平・垂直とも測位誤差を1m以内に抑えられる。 July 2005 Sakai, ENRI 補正性能の傾向 強い磁気嵐が発生していない限り、QZSS 5局のモニタ局配置で、水平・垂直とも測位誤差を1m以内に抑えられる。 ユーザ側は一周波のコード擬似距離のみを使用。 量子化単位は 12.5cm で、SBASと同じまま。 水平0.4～0.7m、垂直0.6～0.9m程度の測位誤差（RMS値）。 強い磁気嵐のもとでは、補正情報に改善の余地がある。 モニタ局の配置については、MSAS配置のほうが若干有利。 QZSSに比べてコンパクトにモニタ局が配置されているため。 ただし、サービスエリア端部での測位精度は要検証。 QZSSは、GEONETを利用してモニタ局を追加できる。 実際には１秒サンプルの観測データが利用でき、キャリアスムージングが有効に作用して測位精度は若干改善する。

July 2005 Sakai, ENRI プロテクションレベル（水平） プロテクションレベル 現実の測位誤差

July 2005 Sakai, ENRI プロテクションレベル（垂直） プロテクションレベル 現実の測位誤差

July 2005 Sakai, ENRI Conclusion 準天頂衛星（測位ミッション）が放送するサブメータ級補強信号の検討のため、プロトタイプシステムを作成してユーザ測位誤差を評価した。 SBAS方式補強メッセージにより、所要の測位精度は達成可能： ユーザ側は一周波のコード擬似距離のみを使用。 強い磁気嵐が発生していない限り、水平0.4～0.7m、垂直0.6～0.9m程度の測位誤差（RMS値）。 強い磁気嵐のもとでは、補正情報に改善の余地がある。 今後の検討課題： １秒サンプルによる評価（キャリアスムージングが有効に働く） モニタ局を追加（GEONETの利用を想定） 磁気嵐の際の補強性能改善