Featured image of post 人造衛星巨型星座與軌道力學:覆蓋人類天空的序幕

人造衛星巨型星座與軌道力學:覆蓋人類天空的序幕

本文探討人造衛星巨型星座(Mega-Constellation)背後的軌道力學、通訊物理學、硬體推進系統以及太空垃圾等議題,以詳細的數學與物理方法深入解析 Starlink 等系統的運作原理。

前言:覆蓋人類天空的巨型星座之序幕

在21世紀的太空開發中,帶來最具野心且戲劇性變革的,莫過於「人造衛星巨型星座(Mega-Constellation)」。以 SpaceX 主導的 Starlink 為首,加上 OneWeb、Amazon 的 Project Kuiper 等,數千至數萬顆前所未見規模的人造衛星群正準備覆蓋地球的低軌道。這不僅僅是通訊技術的進化,更是人類試圖以數學與物理學的極限精度,來設計並將宇宙空間這個三維畫布置於控制之下的宏大嘗試。 本文將探討促成巨型星座成立的底層理論「軌道力學(Orbital Mechanics)」、太空通訊基礎的電磁波傳播、硬體的推進系統,以及威脅太空環境永續性的太空垃圾問題,透過極其詳細的數理方法進行徹底的解析。

第1章:地球同步軌道(GEO)的極限與低地球軌道(LEO)巨型星座的典範轉移

1.1 地球同步軌道(GEO)通訊的物理限制

利用宇宙空間的通訊系統,長期以來一直以位於赤道上空約 35,786 公里的地球同步軌道(Geostationary Earth Orbit: GEO)為主角。由於 GEO 的地球自轉週期(恆星日:約 23 小時 56 分 4 秒)與衛星的公轉週期完全一致,從地面上看起來,衛星始終靜止在同一個方向。憑藉這項特性,地面天線可以保持固定而不需要追蹤機構,且僅需一顆衛星即可覆蓋廣大地區,具有極大的優勢。

$$ 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} $$

也就是說,與 GEO 衛星相比,LEO 衛星在相同頻段下的電波衰減少了約 36.3 dB(功率比約為 4,200 倍)。如此一來,在大幅降低衛星端發射功率(EIRP)的同時,也能夠縮小使用者終端的天線孔徑面。這種強大的鏈路預算,使得僅靠直徑約 50 公分的家用天線進行寬頻通訊成為可能。

第2章:軌道力學的數理與攝動理論

為了讓數萬顆衛星在互不碰撞的情況下,持續且無縫地覆蓋地球上任何地點,精密的數理模型是不可或缺的。在此,我們將詳細探討從二體問題到攝動理論(Perturbation Theory)。

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)。這樣猛烈的速度,產生了後文即將提到的都卜勒頻移(Doppler Shift)以及極高頻率的換手(Handover)。

2.2 克卜勒軌道六要素 (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$ 是地心緯度。影響最大的,是代表赤道膨脹的二次帶諧係數 $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)$ 是半通徑(Semi-latus rectum),$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 的部分星殼(Shell)為了覆蓋極地(北極、南極周邊),採用了接近 SSO 的極軌道。

第3章:沃克星座幾何學

用數千顆衛星毫無縫隙地覆蓋整個地球的幾何學最佳解,就是「沃克星座(Walker Constellation)」。

3.1 Walker-Delta 配置 $i: T/P/F$ 的數學定義

由約翰·G·沃克(John G. Walker)所設計的 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 度附近)使衛星之間的碰撞風險降至最低,同時保證了從地面仰望時,始終有一顆以上衛星處於仰角 25 度以上的連續覆蓋(Continuous Coverage)。

3.2 透過星際光學鏈路(ISL)建立的太空網狀網路

第一代的巨型星座,只能在衛星能同時與地面使用者終端以及閘道站(地面站)進行通訊的範圍內(彎管連接,Bent-pipe)提供網際網路。這導致無法在海洋中央或極地提供服務。

突破此一限制的,是使用雷射通訊的星際光學鏈路(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}$ 是鏈路剩餘可維持時間。 資料封包以光速在宇宙空間中直線躍點(Hop)。與沿著地球曲率鋪設光纖的地面網路相比,其路徑長度較短,且沒有折射率所帶來的延遲(光纖中為 $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 在第一代採用了氪(Krypton),在第二代(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}$(十萬分之一),衛星便會自主點燃霍爾推進器執行迴避機動。透過結合 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 的唐納德·凱斯勒(Donald Kessler)在 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$,這讓像是魯賓天文台(Rubin Observatory)這類超高感度廣角望遠鏡的 CCD 感測器達到飽和(Saturation),並引起嚴重的串擾。這對近地小行星(NEO)的探測與宇宙學觀測帶來了致命的影響。

6.2 遮光對策:遮陽傘衛星與介電質反射膜

SpaceX 與天文學界展開合作尋求對策。

  1. 黑衛星 (DarkSat): 表面塗黑,但會吸收太陽熱能導致熱設計失效。
  2. 遮陽傘衛星 (VisorSat): 透過展開式遮陽板製造陰影,但會與雷射通訊設備產生干擾,且大氣阻力也隨之增加。
  3. 介電質反射膜: 在第二代(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 進行長距離回程(Backhaul)通訊,來將本來會白白浪費的「海洋上空衛星容量」變現,是商業模式上極其重要的命題。

comments powered by Disqus