Featured image of post 人工衛星メガコンステレーションと軌道力学:数万機の衛星が織りなす低軌道通信網の数理と現実

人工衛星メガコンステレーションと軌道力学:数万機の衛星が織りなす低軌道通信網の数理と現実

Starlinkが切り拓いたメガコンステレーション。ウォーカー・コンステレーション、ケプラー軌道要素、大気ドラッグと軌道離脱、ケスラーシンドロームの危機。

はじめに:人類の空を覆うメガコンステレーションの幕開け

21世紀の宇宙開発において、最も野心的かつ劇的な変革をもたらしているのが「人工衛星メガコンステレーション(Mega-Constellation)」である。SpaceXが主導するStarlinkをはじめ、OneWeb、AmazonのProject Kuiperなど、数千から数万機というかつてない規模の人工衛星群が地球の低軌道を覆い尽くそうとしている。これは単なる通信技術の進化にとどまらず、人類が宇宙空間という三次元のキャンバスを、数学と物理学の極限の精度でデザインし、統制下に置こうとする壮大な試みである。 本稿では、メガコンステレーションを成立させている根底の理論である「軌道力学(Orbital Mechanics)」、宇宙通信の基礎である電磁波伝搬、ハードウェアの推進系、そして宇宙環境の持続可能性を脅かすデブリ問題に至るまで、極めて詳細な数理的アプローチを用いて徹底的に解き明かしていく。

第1章:静止軌道(GEO)の限界と低軌道(LEO)メガコンステレーションへのパラダイムシフト

1.1 静止軌道(GEO)通信の物理的制約

宇宙空間を利用した通信システムは、長らく赤道上空約35,786 kmに位置する静止軌道(Geostationary Earth Orbit: GEO)が主役であった。GEOは地球の自転周期(恒星日:約23時間56分4秒)と衛星の公転周期が完全に一致するため、地上からは常に同じ方角に静止しているように見える。この特性により、地上のアンテナは固定されたまま追尾機構を必要とせず、広大な地域を1機の衛星でカバーできるという絶大なメリットがあった。

$$ t_{delay} = \frac{4 \times 35,786,000}{3 \times 10^8} \approx 0.477 \text{ s} = 477 \text{ ms} $$

実際の通信プロトコルにおける処理遅延、エラー訂正符号(FEC)の計算時間、および地上ネットワークのルーティング遅延を加味すると、往復遅延時間は容易に600msから800msに達する。これは、リアルタイム性が要求される現代のオンラインゲーム、高頻度取引(HFT)、遠隔医療、あるいはスムーズなビデオ会議において致命的な遅延である。TCP/IPプロトコルにおけるウィンドウサイズの制約(BDP: Bandwidth-Delay Product)も相まって、GEOでは広帯域であってもスループットが劇的に低下する問題が存在した。

1.2 低軌道(LEO)へのパラダイムシフト

$$ t_{LEO\_delay} = \frac{4 \times 550,000}{3 \times 10^8} \approx 0.0073 \text{ s} = 7.3 \text{ ms} $$

地上ネットワークのルーティング遅延を含めても20〜30msという、地上光ファイバー網に匹敵する、あるいは長距離通信においてはそれを凌駕する超低遅延インターネットが実現可能となる。ガラスファイバー内の光の屈折率は約1.5であり、光速の約2/3(約 $2 \times 10^8$ m/s)に低下する。一方、真空の宇宙空間では光速は維持されるため、大陸間通信のような数千キロの距離では、LEOを経由した方が物理的により早くデータが到達するのだ。

1.3 フリスの伝搬公式(Friis Transmission Equation)とリンクバジェット解析

$$ P_r = P_t G_t G_r \left( \frac{\lambda}{4 \pi d} \right)^2 \frac{1}{L_a L_s} $$$$ L_{FSPL} = \left( \frac{4 \pi d}{\lambda} \right)^2 $$$$ L_{FSPL}(dB) = 20 \log_{10}(d) + 20 \log_{10}(f) + 20 \log_{10}\left(\frac{4 \pi}{c}\right) $$$$ \Delta L = 20 \log_{10}\left(\frac{36000}{550}\right) \approx 20 \log_{10}(65.45) \approx 36.3 \text{ dB} $$

すなわち、LEO衛星はGEO衛星と比較して、同じ周波数帯において約36.3dB(電力比で約4,200倍)も電波の減衰が少ない。これにより、衛星側の送信電力(EIRP)を大幅に抑えつつ、ユーザー端末のアンテナ開口面を小型化できる。この強力なリンクバジェットが、直径わずか50cm程度の家庭用アンテナでのブロードバンド通信を可能にした。

第2章:軌道力学の数理と摂動論

数万機の衛星が互いに衝突することなく、かつ地球上のあらゆる地点をシームレスにカバーし続けるためには、精密な数理モデルが不可欠である。ここでは二体問題から摂動論までを詳細に追う。

2.1 ケプラーの法則と二体問題方程式

$$ m \frac{d^2\mathbf{r}}{dt^2} = -G \frac{Mm}{r^3} \mathbf{r} $$$$ \ddot{\mathbf{r}} + \frac{\mu}{r^3} \mathbf{r} = 0 $$$$ v^2 = \mu \left( \frac{2}{r} - \frac{1}{a} \right) $$

高度550 kmの円軌道($a = 6371 + 550 = 6921$ km, $r=a$)の場合、速度は $v = \sqrt{\mu/a} \approx 7.59 \text{ km/s}$(時速約27,300 km)となる。この猛烈な速度が、後述するドップラーシフトとハンドオーバーの極端な頻度を生み出す。

2.2 ケプラーの軌道6要素 (Keplerian Elements)

三次元空間において衛星の軌道と位置を完全に特定するためには、6つの独立したパラメータが必要である。

  1. 軌道長半径 (Semi-major axis, $a$): 軌道のエネルギーと周期を決定する。
  2. 離心率 (Eccentricity, $e$): 軌道の形状(円軌道なら $e=0$)。通信品質を一定に保つため、LEOコンステレーションは $e \approx 0.0001$ 程度の極めて真円に近い軌道をとる。
  3. 軌道傾斜角 (Inclination, $i$): 赤道面と軌道面がなす角度。Starlinkでは53度、70度、97.6度などが採用される。
  4. 昇交点赤経 (Right Ascension of the Ascending Node, $\Omega$): 春分点(Vernal Equinox)方向から昇交点(衛星が赤道を南から北へ横切る点)までの赤道面上の角度。
  5. 近地点引数 (Argument of Perigee, $\omega$): 昇交点から近地点までの軌道面上の角度。
  6. 真近点角 (True Anomaly, $\nu$): 近地点から測った現在の衛星の位置を表す角度。

2.3 地球の扁平率に伴う重力ポテンシャル $J_2$ 摂動

$$ U = \frac{\mu}{r} \left[ 1 - \sum_{n=2}^{\infty} J_n \left(\frac{R_e}{r}\right)^n P_n(\sin \phi) \right] $$

ここで、$R_e$ は地球の赤道半径(6378.137 km)、$P_n$ はルジャンドル多項式、$\phi$ は地心緯度である。最も影響が大きいのが、赤道の膨らみを表す2次の帯状調和係数 $J_2 \approx 1.08263 \times 10^{-3}$ である。

$$ \dot{\Omega} = -\frac{3}{2} J_2 \left(\frac{R_e}{p}\right)^2 n \cos i $$$$ \dot{\omega} = \frac{3}{4} J_2 \left(\frac{R_e}{p}\right)^2 n (5 \cos^2 i - 1) $$

ここで、$p = a(1-e^2)$ は半直弦、$n = \sqrt{\mu/a^3}$ は平均運動である。

軌道傾斜角 $i$ が90度より小さい場合(順行軌道)、$\dot{\Omega}$ は負となり、軌道面は地球の自転と逆行する形で西へ回転する(ノード退行:Nodal Regression)。高度550 km、傾斜角53度の場合、$\dot{\Omega}$ は約 $-5.2^\circ / \text{day}$ となる。メガコンステレーションでは、すべての衛星が同じ高度と傾斜角を持つように精密に制御される。これにより、$J_2$ 摂動による $\dot{\Omega}$ の変化率が全プレーンで同一となり、コンステレーションの網目構造が長期間にわたって相対的な形を崩さずに維持されるのである。

2.4 太陽同期軌道(SSO)の設計原理

$$ 0.9856^\circ/\text{day} = -\frac{3}{2} J_2 \left(\frac{R_e}{a}\right)^2 n \cos i $$

これを満たすためには $\cos i < 0$、すなわち軌道傾斜角 $i > 90^\circ$ の逆行軌道でなければならない。高度550 kmの場合、$i \approx 97.6^\circ$ となる。Starlinkの一部のシェルは極地(北極・南極周辺)をカバーするためにSSOに近い極軌道を採用している。

第3章:ウォーカー・コンステレーション幾何学

数千もの衛星で地球全体を隙間なくカバーする幾何学的な最適解が「ウォーカー・コンステレーション(Walker Constellation)」である。

3.1 Walker-Delta配置 $i: T/P/F$ の数学的定義

ジョン・G・ウォーカーが考案したWalker-Deltaパターンは、$i: T/P/F$ という記法で完全に定義される。

  • $i$: 軌道傾斜角 (Inclination)
  • $T$: コンステレーションを構成する総衛星数
  • $P$: 軌道面の数 (Number of orbital Planes)
  • $F$: 隣接軌道面間での衛星の位相差パラメータ ($0 \le F \le P-1$ の整数)

各軌道面には $S = T/P$ 機の衛星が均等に配置される。軌道面内の衛星間隔は $\Delta \nu = 360^\circ / S$ である。 昇交点赤経 $\Omega$ は赤道上で均等に分割され、隣接軌道面との間隔は $\Delta \Omega = 360^\circ / P$ となる。 さらに、隣接する東側の軌道面における衛星の真近点角のズレ(位相差)は $\Delta \Phi = F \times (360^\circ / T)$ で与えられる。

例えば、Starlinkの第一世代における代表的なシェル(Shell 1)は、高度550 km、傾斜角53度で、$T=1584, P=72$ という巨大なWalker配置を採用している(各面に $S=22$ 機)。隣接プレーン間での位相差 $F$ を最適化することで、軌道が最も密集する最高緯度(北緯・南緯53度付近)における衛星同士の衝突リスクを最小化しつつ、地上から見上げた際に常に1機以上の衛星が仰角25度以上に見える連続カバレッジ(Continuous Coverage)を保証している。

3.2 衛星間光リンク(ISL)による宇宙メッシュネットワーク

第一世代のメガコンステレーションは、衛星が地上のユーザー端末とゲートウェイ局(地上局)の両方と同時に通信できる範囲(ベントパイプ接続)でしかインターネットを提供できなかった。これでは大洋の真ん中や極地にサービスを提供できない。

この限界を突破するのが、レーザー通信を用いた衛星間光リンク(Inter-Satellite Link: ISL)である。真空の宇宙空間では大気による光の減衰やシンチレーション(大気の揺らぎ)が存在しないため、波長1.55 $\mu$m 帯(Cバンド)のレーザーによる数Gbps〜数十Gbpsの大容量・低遅延通信が可能となる。 各衛星は4基の光通信ターミナルを搭載し、同じ軌道面内の前後2機(Intra-plane ISL)と、隣接する軌道面の左右2機(Inter-plane ISL)とレーザーリンクを確立する。

3.3 最短経路ダイクストラ法とトポロジー動的更新

ISLによって形成されるネットワークは、ノード(衛星)が秒速約7.5 kmで移動するため、ネットワークトポロジーが秒単位で激しく変化する。特に極地方に向かうにつれて軌道面が交差するため、隣接プレーンの衛星とのレーザーリンク(Inter-plane ISL)は定期的に切断と再接続(ハンドオーバー)を繰り返す。

$$ C_{ij} = \alpha \cdot d_{ij} + \beta \cdot Q_{ij} + \gamma \cdot L_{ij} $$

ここで、$d_{ij}$ は物理的距離(遅延)、$Q_{ij}$ はキュー長(輻輳)、$L_{ij}$ はリンクの残り維持可能時間である。 データパケットは、宇宙空間を光速のまま直線的にホップしていく。地球の曲率に沿って光ファイバーを這わせる地上ネットワークと比較して、経路長が短く、屈折率による遅延(ファイバー中は $c/1.5$)がないため、ニューヨーク・ロンドン間のような超長距離通信では、ISLを経由した方が理論的に高速になるのである。

第4章:Starlink衛星のハードウェアと推進系

メガコンステレーションの成立条件は、衛星の大量生産と劇的なコストダウンである。

4.1 クリプトン/アルゴンホールスラスタの比推力と推進薬質量計算

衛星は軌道投入後、自力で運用軌道まで上昇し、運用中は大気ドラッグを補償し、寿命末期には軌道離脱(Deorbit)を行う必要がある。これらのマヌーバに必要な速度増分 $\Delta V$ を達成するため、ホールスラスタ(Hall-effect Thruster)と呼ばれる電気推進システムが採用されている。

$$ m_p = m_0 \left( 1 - e^{-\frac{\Delta V}{I_{sp} g_0}} \right) $$

ここで、$I_{sp}$ は比推力、$g_0$ は標準重力加速度、$m_0$ は初期質量である。 従来の電気推進は高価なキセノン(Xenon)を使用していたが、SpaceXは第1世代でクリプトン(Krypton)、第2世代(V2 Mini)ではアルゴン(Argon)を採用した。アルゴンは大気中に豊富に存在し極めて安価だが、イオン化エネルギーが高く推力効率が落ちる。しかし、磁場トポロジーの最適化により、アルゴンホールスラスタは比推力 2500秒 ≒ 24.5 km/s の排気速度を達成し、大量打ち上げ時の推進薬コストを革命的に引き下げた。

4.2 フェーズドアレイアンテナのビームフォーミング数理

地上端末との通信には、機械的な可動部を持たずに電波のビーム方向を瞬時に変えられるフェーズドアレイアンテナ(Phased Array Antenna)が使用される。 アンテナ素子を格子状に配列し、各素子から送信される電波の位相(Phase)を制御することで、干渉効果によって特定の方向に強い電波ビームを形成する。

$$ \Delta \Phi_{mn} = -\frac{2\pi}{\lambda} (x_m \sin\theta \cos\phi + y_n \sin\theta \sin\phi) $$

Starlink衛星とユーザーターミナルは、高度なビームフォーマーICを搭載し、毎秒数千回という速度で位相ウェイト行列を再計算している。これにより、上空を高速移動する衛星を電子的かつシームレスに追尾し続けることができる。

4.3 自動衝突回避アルゴリズム

$$ P_c = \frac{1}{2\pi |C_p|^{1/2}} \iint_{A} \exp\left( -\frac{1}{2} \mathbf{r}^T C_p^{-1} \mathbf{r} \right) dx dy $$

ここで $C_p$ は投影された共分散行列、$A$ は衝突断面積である。$P_c$ が $10^{-5}$(10万分の1)を超えた場合、衛星は自律的にホールスラスタを点火し回避マヌーバを実行する。AIと予測制御(MPC)の融合により、人手を介さない安全確保が実現されている。

第5章:宇宙デブリ問題とケスラーシンドロームの恐怖

5.1 衝突確率のポアソン過程モデル

$$ d\lambda = \rho \cdot A \cdot v_{rel} \cdot dt $$$$ P_c = 1 - e^{-\int_0^T \rho A v_{rel} dt} $$

低軌道での正面衝突では相対速度が約 $10 \sim 15 \text{ km/s}$ に達する。わずか1cmのアルミニウム破片でも、その運動エネルギーは手榴弾に匹敵し、衛星を完全に粉砕する。

5.2 大気ドラッグによる自然落下メカニズム

$$ \mathbf{a}_{drag} = -\frac{1}{2} \rho_{atm} \frac{C_D A}{m} v_{rel}^2 \frac{\mathbf{v}_{rel}}{v_{rel}} $$

大気密度 $\rho_{atm}$ は高度が下がるほど指数関数的に増加し、太陽活動による極端紫外線(EUV)で熱圏が膨張するとさらに上昇する。Starlinkが高度550 kmを選んだ最大の理由はここにある。万が一制御不能になっても、数年以内(通常1〜5年)に大気ドラッグで自然に高度が低下し、大気圏に突入して燃え尽きる「自己浄化軌道」だからである。高度1000 km以上では数百年滞留してしまう。

5.3 軌道上衝突による連鎖破壊:ケスラーシンドローム

NASAのドナルド・ケスラーが1978年に提唱した「ケスラーシンドローム(Kessler Syndrome)」は最悪のシナリオである。 大型物体が衝突すると数千個のデブリ雲が発生し、それが他の衛星に衝突する確率を急激に引き上げる。連鎖的に衝突が発生し、破片が指数関数的に増殖する。 臨界密度を超えると、新たな打ち上げを行わずとも自己増殖が止まらなくなり、特定の軌道領域(例えば高度700〜1,000 km)が数百年から数千年にわたって使用不能になる。FCCが「運用終了後5年以内の軌道離脱」を義務付けた背景には、この破滅的連鎖を未然に防ぐ強い危機感がある。

第6章:天文学への光害問題と宇宙の持続可能性

6.1 衛星の反射光度と光学望遠鏡への影響

$$ m_1 - m_2 = -2.5 \log_{10} \left( \frac{F_1}{F_2} \right) $$

初期のStarlink衛星は視等級で $+3$ 〜 $+5$ に達し、ルビン天文台のような超高感度広視野望遠鏡のCCDセンサーを飽和(サチュレーション)させ、深刻なクロストークを引き起こした。これにより、地球接近小惑星(NEO)の探査や宇宙論的観測に致命的な影響が出た。

6.2 遮光対策:バイザーサットと誘電体ミラーフィルム

SpaceXと天文学界は協働で対策に乗り出した。

  1. ダークサット (DarkSat): 表面を黒く塗装したが、太陽熱を吸収し熱設計が破綻した。
  2. バイザーサット (VisorSat): 展開式サンバイザーで影を作ったが、レーザー通信機器と干渉し、大気抵抗も増大した。
  3. 誘電体ミラーフィルム: 第2世代(V2 Mini)では、光を地上ではなく宇宙空間へ鏡面反射させる特殊なブラッグ反射膜(Dielectric Mirror Film)と、黒色塗料を組み合わせた高度な熱・光学制御を採用。これにより、肉眼で見えない $+7$ 等級以下まで暗くすることに成功しつつある。

6.3 宇宙交通管理(STM)の未来

数万機の人工衛星が織りなす低軌道メガコンステレーションは、全人類にブロードバンドをもたらす革命である。しかしそれは同時に、軌道力学の過酷な数理と、ケスラーシンドロームという宇宙環境の限界に対する、人類のモラルの試金石でもある。 現在、国連COPUOSを中心に、海洋における「航行の自由」や「COLREG」に相当する宇宙交通管理(Space Traffic Management: STM)の枠組み作りが急ピッチで進められている。持続可能な宇宙開発(Space Sustainability)の実現こそが、次世代に対する我々の最大の責務である。

付録:メガコンステレーションの通信容量モデリング

$$ C_{beam} = B \log_2 \left( 1 + \text{SINR} \right) $$

ここで、$B$ は帯域幅(例えばKuバンドにおける 250 MHz などのチャネル幅)、SINR(Signal-to-Interference-plus-Noise Ratio)は信号対干渉雑音比である。

$$ \text{SINR} = \frac{P_r}{N_0 B + \sum I_{intra} + \sum I_{inter}} $$
  • $P_r$: 受信電力(フリスの伝搬公式から算出)
  • $N_0$: 雑音電力密度($N_0 = k T_{sys}$、$k$ はボルツマン定数、$T_{sys}$ はシステム雑音温度)
  • $\sum I_{intra}$: 同一システム内の他のビームまたは他の衛星からの自己干渉(Intra-system interference)
  • $\sum I_{inter}$: OneWebやKuiperなど、他社のコンステレーションやGEO衛星からの干渉(Inter-system interference)
$$ C_{sat} = \sum_{i=1}^{N_{beam}} B_i \log_2 \left( 1 + \text{SINR}_i \right) $$$$ C_{eff} = C_{total} \times \eta_{land} \times \eta_{pop} \times \eta_{utilization} $$

ここで $\eta_{utilization}$ はネットワークの稼働率およびルーティング効率である。 この数式から明らかなように、メガコンステレーションの経済的自立性を高めるためには、海洋上の航空機や船舶へのサービス提供、あるいはISLを用いた長距離バックホール通信によって、本来無駄になるはずの「海洋上空の衛星キャパシティ」をいかにマネタイズするかが、ビジネスモデル上の極めて重要な命題となるのである。

comments powered by Disqus