行星防禦進階：β 係數、重力鎖孔與時間槓桿的精密力學

一顆撞不中的小行星，為什麼還是讓人睡不著？

想像有一顆直徑 370 公尺的小行星，雷達與光學觀測都告訴你：它在 2029 年不會撞上地球。它會從距離地表約 3.2 萬公里處呼嘯而過——比同步衛星還近，肉眼可見。聽起來是虛驚一場。但對行星防禦的工程師而言，這次「擦身而過」反而是最讓人緊張的劇本之一。因為這顆名為 Apophis（阿波菲斯，99942）的小行星，會在掠過時被地球重力「捏」一把，軌道被微微改寫。問題是：它會不會剛好穿過某個極狹窄的「重力鎖孔（gravitational keyhole）」，導致幾年後的某次回歸真的撞上來？

入門篇談過我們如何「看見」小行星、如何用動能撞擊器（kinetic impactor）把它推開。但真正困難的從來不是「推」，而是「推多少、推得準不準、什麼時候推」。這篇進階篇要把這幾個問句拆成可以計算的物理量：動量轉移的放大係數 $\beta$、鎖孔與共振回歸的幾何、以及偏轉量如何隨「提前期」非線性放大。讀完你會明白，行星防禦本質上是一場與槓桿、誤差與時間賽跑的精密力學。

動量轉移不是只有「撞」那麼簡單：β 係數

入門篇說動能撞擊器靠動量守恆把小行星推開。但如果你真的只用「撞擊器質量 × 速度」去算，你會嚴重低估效果——有時低估到兩三倍。原因藏在一個叫做 $\beta$（beta factor）的放大係數裡。

當撞擊器以高速撞上小行星表面，它不只把自己的動量交給目標。撞擊會炸出大量碎屑（ejecta），這些碎屑以相當高的速度從撞擊坑向後噴出。根據動量守恆，噴出的碎屑就像火箭尾氣，會給小行星一個額外的反作用推力。於是小行星實際獲得的動量，比撞擊器本身攜帶的還要大。

我們用一個無量綱數 $\beta$ 來描述這個放大效果：

$$ m_{\text{ast}}\,\Delta v_{\text{ast}} = \beta \, m_{\text{imp}}\,v_{\text{imp}} $$

其中 $m_{\text{ast}}$ 是小行星質量，$\Delta v_{\text{ast}}$ 是它獲得的速度變化，$m_{\text{imp}}$ 與 $v_{\text{imp}}$ 是撞擊器的質量與相對撞擊速度。

• 若 $\beta = 1$：完全沒有碎屑貢獻，純粹是「撞擊器埋進去」的理想塑性碰撞。
• 若 $\beta > 1$：碎屑反衝提供額外推力。$\beta = 3$ 代表小行星實際得到的動量是撞擊器本身的三倍。

$\beta$ 的大小取決於目標的物理性質：表面是堅硬岩石還是鬆散的碎石堆（rubble pile）？孔隙率多高？這也正是為什麼我們不能只靠紙上計算——必須真的去撞一次看看。

看一個例子：DART 量到的 β

2022 年 9 月，NASA 的 DART（Double Asteroid Redirection Test）任務撞上了小行星 Didymos 的衛星 Dimorphos。這是人類第一次真正測量 $\beta$。

我們無法直接量 Dimorphos 的 $\Delta v$，但可以量它繞 Didymos 的「公轉週期變化」。撞擊前 Dimorphos 繞行週期約 $11\text{ h }55\text{ min}$。科學家原本設定:只要週期縮短 73 秒就算成功。結果呢?週期縮短了約 33 分鐘——遠超預期。

從週期縮短回推，Dimorphos 沿軌道方向獲得的速度變化約為 $\Delta v \approx 2.7\ \text{mm/s}$。把這個值代回上面的公式，並用 DART 的質量（約 $580\ \text{kg}$）、撞擊速度（約 $6.1\ \text{km/s}$）與 Dimorphos 的估計質量，得到：

$$ \beta \approx 3.6 $$

換句話說，碎屑反衝讓動量轉移放大了約三倍多。這一個數字，把行星防禦從「理論上可行」推進到「我們知道實際效率」。它也告訴我們：對碎石堆型小行星，撞擊可能比想像中更有效——但前提是你撞的地方真的會噴出大量碎屑。

鎖孔與共振回歸：為什麼「沒撞中」不等於「安全」

回到開頭的 Apophis。要理解鎖孔，先要理解一個概念：小行星與地球的軌道常處於「平均運動共振（mean-motion resonance）」附近。

假設某顆小行星的公轉週期讓它每繞太陽 $7$ 圈，地球剛好繞 $8$ 圈，那麼大約每隔一段固定時間，它就會回到地球附近——這叫 $7{:}8$ 共振回歸。如果第一次接近時，小行星穿過了某個位置極為精確的小區域，地球重力會把它的軌道週期改寫到「下一次回歸時剛好撞上」的數值。這個致命的小區域，就是鎖孔（keyhole）。

鎖孔有多小？對 Apophis 的某些回歸情境，鎖孔的寬度只有幾百公尺到不到一公里。在 3 萬公里的接近距離尺度下，這意味著小行星的位置只要差幾百公尺，未來命運就天差地別。這也是為什麼天文學家要用雷達把 Apophis 的軌道量到極致——所幸 2021 年的精密觀測已排除了它在 2068 年前撞擊地球的可能。

鎖孔概念對「偏轉策略」有個反直覺的啟示：我們不需要把小行星推離地球幾千公里，只需要讓它錯過那個幾百公尺寬的鎖孔。 在共振回歸發生前去推它，所需的 $\Delta v$ 小到驚人。但若錯過鎖孔的時機、等到它已經鎖定撞擊軌道才動手，需要的推力會暴增好幾個數量級。

偏轉量如何隨時間放大：提前期才是真正的槓桿

這是整個行星防禦中最重要、也最常被誤解的一點。很多人以為「偏轉就是把小行星推開一段距離」，於是直覺認為推力越大越好。但真正的槓桿不是力，而是時間。

考慮一個簡化模型。假設我們對小行星施加一個沿軌道切線方向的速度增量 $\Delta v$。這個微小的速度改變不會讓小行星「瞬間跳開」，而是會緩慢地改變它的軌道週期，使它在軌道上的位置一點一點地「提前」或「落後」。經過時間 $t$ 後，這個位置偏移（沿軌道方向）大約是：

$$ \Delta x \approx \frac{3 \, \Delta v}{a}\, \cdot v_{\text{orb}}\, \cdot t \;\;\sim\;\; \Delta v \cdot t $$

關鍵在於最後那個比例關係：位置偏移與「施力後經過的時間」成正比。 你越早動手，同樣的 $\Delta v$ 能累積出越大的偏移。提前 20 年動手，效果是提前 2 年動手的約 10 倍。

換個比喻：這就像射飛鏢時手抖了一下。如果在出手瞬間（離靶很近）抖，飛鏢只偏一點；但若那個微小角度誤差從很遠的地方就開始累積，到了靶上就會差很多。行星防禦正是利用這個「遠距槓桿」——用極小的早期推力，換取多年後地球軌道尺度的巨大偏移。

動手算一下：偏一個地球半徑要多少 Δv？

我們來估算：若想讓一顆小行星在 $t = 10$ 年後相對撞擊點偏移一個地球半徑（$R_\oplus \approx 6.37\times10^6\ \text{m}$），大約需要多大的沿軌道 $\Delta v$？

用上面的比例關係的數量級版本 $\Delta x \sim 3\,\Delta v \cdot t$（這個係數來自軌道力學中速度擾動對沿軌道漂移的放大），先把時間換成秒：

$$ t = 10\ \text{年} \approx 3.15\times10^{8}\ \text{s} $$

要求 $\Delta x \approx R_\oplus = 6.37\times10^{6}\ \text{m}$，解出：

$$ \Delta v \sim \frac{\Delta x}{3\,t} = \frac{6.37\times10^{6}}{3 \times 3.15\times10^{8}} \approx 6.7\times10^{-3}\ \text{m/s} $$

也就是大約 $\mathbf{7\ mm/s}$——比走路速度的千分之一還小。注意這正好和 DART 在 Dimorphos 上實際造成的 $\Delta v$（約 $2.7\ \text{mm/s}$）同一個數量級！這告訴我們一件令人振奮的事：只要提前十幾二十年，人類現有的單次動能撞擊技術，原則上就足以偏轉一顆中型小行星。 時間，是我們最強大的武器。

（提醒：這是數量級估算，真實任務還要考慮 $\beta$ 不確定性、撞擊方向誤差、小行星自轉與形狀等因素，所以實務上會留下數倍的安全餘裕。）

為什麼大小決定一切：撞擊能量的尺度律

最後我們補上一個入門篇沒細談的維度：不同大小的小行星，威脅完全不在一個數量級。

撞擊釋放的動能是

$$ E = \tfrac{1}{2} m v^2 = \tfrac{1}{2}\left(\tfrac{\pi}{6}\rho D^3\right) v^2 $$

質量隨直徑的三次方成長。所以直徑增加 10 倍，能量增加 $10^3 = 1000$ 倍。

• 直徑約 $20\ \text{m}$（如 2013 年車里雅賓斯克）：空爆能量約 $0.5\ \text{Mt}$ TNT，震碎窗戶、造成局部傷害。
• 直徑約 $140\ \text{m}$：這是 NASA 訂為「需追蹤」的門檻，撞擊足以摧毀一個都會區。
• 直徑約 $1\ \text{km}$：撞擊能量達數萬 Mt，會造成跨大陸氣候衝擊，屬於「全球文明級」威脅。
• 直徑約 $10\ \text{km}$（如導致恐龍滅絕的希克蘇魯伯撞擊體）：能量約 $10^{8}\ \text{Mt}$，是物種滅絕級事件。

這個三次方尺度律解釋了行星防禦的資源分配邏輯：小天體很多但傷害有限且常在大氣中燒毀；真正要嚴防的是 $140\ \text{m}$ 以上、數量相對少但破壞力呈立方放大的那一群。偵測完整度（survey completeness）的目標，正是針對這個尺寸帶設定的。

重點回顧

• $\beta$ 係數把動能撞擊的效果放大：碎屑反衝像火箭尾氣提供額外推力，DART 實測 $\beta \approx 3.6$，證明對碎石堆型小行星撞擊比理想塑性碰撞更有效。
• 鎖孔（keyhole）是軌道上幾百公尺寬的致命區域；穿過它會讓小行星在未來某次共振回歸時撞上地球。偏轉的真正目標常常是「錯過鎖孔」而非「推離地球」。
• 提前期是最大的槓桿：沿軌道位置偏移與施力後經過的時間成正比，越早動手，同樣的 $\Delta v$ 效果越大。
• 數量級估算顯示，提前 10 年只需約 $7\ \text{mm/s}$ 的 $\Delta v$ 就能偏移一個地球半徑——與 DART 實測同量級，代表現有技術原則上已足夠。
• 撞擊能量隨直徑三次方成長，這條尺度律決定了為何要優先嚴防 $140\ \text{m}$ 以上的小行星。

深入探討（研究所視角）

若要把上述直覺推進到研究層級，核心是把偏轉問題寫成軌道動力學的線性化擾動。在開普勒軌道上施加切向速度增量 $\Delta v_t$，最敏感的響應是半長軸 $a$ 的改變。由活力公式（vis-viva）微分可得 $\delta a = \frac{2a^2 v}{\mu}\,\Delta v_t$，其中 $\mu = GM_\odot$。半長軸改變會透過開普勒第三定律改寫週期，於是沿軌道相位（mean anomaly）以正比於時間的速率漂移——這正是「提前期槓桿」的數學根源，且切向施力遠比徑向施力有效率。嚴謹的任務設計會用 B-plane（目標平面）座標：把小行星在接近地球時的軌跡投影到一個垂直於相對速度向量的平面上，並在此平面上量化偏移量與不確定橢圓（uncertainty ellipse），鎖孔正是 B-plane 上一系列離散的小區域。

更前沿的議題包括：（1）$\beta$ 的物理建模——碎屑動量取決於撞擊角、孔隙率與孔隙坍縮（pore crushing），DART 之後的 ESA Hera 任務正前往 Dimorphos 量測撞擊坑形態、密度與精確質量，以反演 $\beta$ 與目標內部結構；（2）慢推（slow-push）方法的數學框架，如重力拖船（gravity tractor）利用太空船自身重力施加 $\sim\!\text{nN}$ 級連續力，需以連續推力下的軌道積分而非脈衝近似求解，其優勢在於不依賴未知的 $\beta$、可即時調整；（3）Yarkovsky 效應——小行星因不對稱熱輻射產生的微小持續加速度，長期可顯著改寫軌道，是高精度撞擊預測（與偏轉殘差估計）中不可忽略的非重力項。把這些放在一起，行星防禦本質上是一個多尺度、含模型不確定性的最佳控制問題：在偵測完整度、軌道測定精度、$\beta$ 不確定性與發射時程之間，尋找使撞擊機率降到可接受水準的最小代價解。