卫星退休后去哪？3.6 万公里高空的墓地轨道

先别往下翻

（在心里选一个，再往下看。）

3.6 万公里的"车位"为什么值钱

地球同步轨道（GEO）高度约 35,786 公里。在这个高度，卫星绕地球一圈的周期恰好等于地球自转周期（约 23 小时 56 分），于是它看上去"悬停"在地球某一点上方。一颗 GEO 卫星服务的区域永远不变——电视、广播、气象、通信，都靠这个性质。

但 GEO 不是一条线，而是一个"带"。按 NASA 轨道碎片首席科学家 Nicholas L. Johnson 的定义：GEO 保护区 = 35,786 km ± 200 km 高度、南北纬 ±15° 围成的空间。位置是稀缺的：两颗卫星不能占用同一个"车位"。这就是为什么说——轨道是资产。

退役的卫星不能留在原地，也不能随便掉下来

寿命到了，卫星不能留在原地：它会慢慢漂移、晃动，威胁旁边还在工作的"车位"。也不能像低轨卫星那样掉进大气烧掉：GEO 距地面 3.6 万公里，没有大气可以依靠，把它拖回来需要的燃料比推上去多得多，再入还不可控。

所以国际共识是第三条路：把它推到更高的轨道——一条比 GEO 高几百公里的"墓地轨道"（graveyard orbit），也叫超同步处置轨道。规则由两个国际组织给出：

（对照：低轨卫星走的是另一条规则——寿命终了后 25 年内离轨。）

抬升要花多少燃料？两个数就能算出来

这里只需要两个常数：地球引力常数 μ ≈ 398600 km³/s²、地球半径 R ≈ 6378 km。算出来的 GEO 圆轨道速度是约 3.0747 km/s。

要把卫星从 GEO 抬升到处置轨道，用的是霍曼转移——两次点火。脚本对三种抬升量算出总 Δv：

看懂这张表，关键在于：每秒约 3.07 公里的飞行速度下，几十米每秒只是"轻轻推一下"——但这点推力足以把卫星抬到一条永不回落到 GEO 保护区的轨道。燃料成本主要不在这次抬升，而在卫星十几年寿命里"稳在地位"的长期消耗；如果前期把燃料省没了，末了想推也推不动。

IADC 公式按典型通信卫星取值（光压系数 C_R=1.4、面质比 A/m=0.02 m²/kg）算出来是 235 + 1000×1.4×0.02 = 263 km。这正是现实中卫星做出来的数量级。

一个做对了的样本：GSAT-12

2023 年 3 月，印度空间研究组织（ISRO）给服役超过 11 年的通信卫星 GSAT-12 做"后任务处置"：7 次点火，把卫星抬升到比 GEO 高约 400 公里的超同步圆轨道，最低近地点抬升 261 km、偏心率控制在 0.003 以内——与 IADC 公式算出的 263 km 几乎一致。ISRO 声明这次处置完全符合 IADC 与联合国的减缓指南。

一个做错了的样本：EchoStar-7 和全球首张"太空垃圾罚单"

做对不难，做错也容易。美国 DISH 公司的 EchoStar-7 卫星 2022 年退役时，2012 年获批的处置计划要求把它抬到比 GEO 高 300 公里的墓地轨道；结果因为推进剂余量不足，只抬高了 122 公里——比要求低了 178 公里。

2023 年 10 月，美国联邦通信委员会（FCC）开出全球第一张太空碎片罚单：罚款 15 万美元，DISH 同意和解并承认责任。罚款金额不大，但它把"处置承诺"变成了有成本的事——这比道德呼吁管用。

三个问题

1. 卫星退役为什么不能留在原轨道？ 报废卫星会在轨道带上慢慢漂移、晃动，威胁旁边的"车位"。GEO 是稀缺资产，报废物体赖着不走，等于把位置占死——这也是碎片减缓规则存在的意义。

2. 为什么 GEO 卫星不坠入大气烧毁？ GEO 距地面约 3.6 万公里，没有大气可以依靠；把它拖回地球需要的燃料比推上去多得多，再入还不可控。所以唯一的可行处置是向上推——推到更高的墓地轨道。

3. 墓地轨道能一直用下去吗？ 能放很久——离 GEO 保护区足够远、几乎不受扰动，规则设计就是要让报废卫星离保护区越远越好。真正的难题不是技术，而是"愿不愿意在寿命末期留够燃料"——EchoStar-7 就是燃料没留够的代价。

你可以自己验

到这里的所有数字，都不是"我们说的"，而是公开文献算出来的：

• GEO 保护区定义（35786 ± 200 km、±15° 纬度）逐字取自 NASA 轨道碎片首席科学家 Nicholas L. Johnson 的论文（NTRS 20100007939）；
• IADC 抬升公式 235 + 1000×C_R×A/m、ITU 的 300 km、低轨 25 年规则，逐字取自联合国 UNOOSA 文件 A/AC.105/820（2003）；
• GSAT-12 处置来自 ISRO 官网新闻稿；EchoStar-7 罚单来自 FCC 执法局 Consent Decree。

脚本读取冻结快照 → 用 μ 与 R 重算 GEO 速度、周期、三种抬升量的 Δv 与近地点高度、IADC 公式验算 → 输出与正文完全一致的表格与哈希。你（或任何 AI agent）都可以重跑核对，见附录。

附录：可复现验证

python3 verify_graveyard.py                          # 读默认快照，重跑并核对
python3 verify_graveyard.py graveyard_orbit.json     # 指定快照文件

脚本会做三件事：

1. 读取冻结快照（GEO 保护区定义 + IADC/ITU 规则 + 工作实例 + 罚单 + 常数）；
2. 重算派生数字（GEO 速度 3.0747 km/s、周期约 23 小时 56 分、三种抬升量 Δv 8.53/10.88/14.48 m/s、处置轨道近地点高度 36021/36086/36186 km、IADC 公式 263 km）；
3. 输出与正文完全一致的表格、验证层级与快照 sha256。

说明：四个 primary 锚点逐字取自公开文献——GEO 保护区定义（NASA NTRS 20100007939）与 IADC 抬升公式、ITU 300 km、低轨 25 年（均出自联合国 UNOOSA 文件 A/AC.105/820，2003-12-15，印度对 IADC 建议的意见）。每个 primary 锚点内含原文引用句 quote 与出处定位 locator，可逐字复核；GSAT-12 处置（ISRO 新闻稿）与 EchoStar-7 罚单（FCC Consent Decree DA-23-888A1）为 secondary——归因清楚的官方二手口径。正文表格与派生数字由脚本从快照重算，冻结快照可精确复现。

原始材料（可直接复制）

下面的脚本与数据快照全文，供你（或任何 AI agent）直接复制使用，无需下载。

复现材料（含 sha256 锚点，随嵌入文件自动更新）：
verify_graveyard.py（复现脚本） · sha256 d566c19f22ff4b3a15d4ef6c6f6238d240f9ba69af951d38a77b15900498a819
graveyard_orbit.json（数据快照） · sha256 550e7c9af0b3bfc9dafd6816df73b1cbfd82fb3df89b238786ccc01410ec5275

verify_graveyard.py — 复现脚本全文

#!/usr/bin/env python3
# verify_graveyard.py — 科普附录：验证墓地轨道（GEO 处置）派生数字（零依赖）
# 输入：graveyard_orbit.json 冻结快照 v2（4 个 primary 锚点 + 常数 + 工作实例 + 罚单）
# 输出：GEO 基本量 + 霍曼转移 Δv + IADC 公式验算 + 证据层级 + 快照 sha256
# 用法：python3 verify_graveyard.py            （默认读 graveyard_orbit.json）
#       python3 verify_graveyard.py FILE       （指定快照文件）
# 退出码：0 成功；2 解析失败；3 证据契约违规（primary 缺 quote/locator）
import hashlib
import json
import math
import sys

TIERS = ("primary", "secondary", "media")


def check_schema(snap):
    """证据门禁：primary 档必须带原文引用句 quote + 出处定位 locator（可复现≠可验证）。"""
    errors = []
    for a in snap.get("anchors", []):
        tier = a.get("tier")
        if tier not in TIERS:
            errors.append(f"锚点 {a.get('label', '?')}: tier 缺失或非法（{tier}）")
            continue
        if tier == "primary":
            if not a.get("quote", "").strip():
                errors.append(f"锚点 {a.get('label', '?')}: primary 档缺 quote（原文引用句）")
            if not a.get("locator", "").strip():
                errors.append(f"锚点 {a.get('label', '?')}: primary 档缺 locator（出处定位）")
    for key, item in snap.get("background", {}).items():
        if isinstance(item, dict) and "tier" in item:
            tier = item.get("tier")
            if tier not in TIERS:
                errors.append(f"background.{key}: tier 非法（{tier}）")
            elif tier == "primary":
                if not item.get("quote", "").strip():
                    errors.append(f"background.{key}: primary 档缺 quote")
                if not item.get("locator", "").strip():
                    errors.append(f"background.{key}: primary 档缺 locator")
    return errors


def hohmann_dv_mps(r1_km, r2_km, mu_km3_s2):
    """圆轨道 r1 → r2 的霍曼转移总 Δv（两次点火之和），单位 m/s。"""
    v1 = math.sqrt(mu_km3_s2 / r1_km)
    v2 = math.sqrt(mu_km3_s2 / r2_km)
    a = (r1_km + r2_km) / 2.0
    v_p = math.sqrt(mu_km3_s2 * (2.0 / r1_km - 1.0 / a))
    v_a = math.sqrt(mu_km3_s2 * (2.0 / r2_km - 1.0 / a))
    return ((v_p - v1) + (v2 - v_a)) * 1000.0


def fmt_hms(seconds):
    h = int(seconds // 3600)
    m = int((seconds % 3600) // 60)
    s = int(seconds % 60)
    return f"{h}小时{m:02d}分{s:02d}秒"


def main():
    path = sys.argv[1] if len(sys.argv) > 1 else "graveyard_orbit.json"
    with open(path, "rb") as f:
        raw = f.read()
    sha = hashlib.sha256(raw).hexdigest()

    try:
        snap = json.loads(raw.decode("utf-8"))
        anchors = snap["anchors"]
        c = snap["constants"]
        ex = snap["worked_example"]
        en = snap["enforcement"]
    except Exception as e:
        print(f"ERROR: 快照解析失败: {e}", file=sys.stderr)
        return 2

    schema_errors = check_schema(snap)
    if schema_errors:
        print("ERROR: 证据契约违规（可复现≠可验证）：", file=sys.stderr)
        for e in schema_errors:
            print(f"  - {e}", file=sys.stderr)
        return 3

    mu = c["mu_km3_s2"]
    R_E = c["R_E_km"]
    GEO_alt = c["GEO_alt_km"]
    C_R = c["C_R"]
    A_m = c["A_over_m_m2_kg"]

    r_GEO = R_E + GEO_alt
    v_GEO = math.sqrt(mu / r_GEO)          # km/s
    T_GEO = 2 * math.pi * math.sqrt(r_GEO ** 3 / mu)  # s

    print("=== 墓地轨道（GEO 处置）派生数字 ===")
    print(f"GEO 地心半径 r = {R_E} + {GEO_alt} = {r_GEO:.1f} km")
    print(f"GEO 圆轨道速度 v = √(μ/r) = {v_GEO:.4f} km/s")
    print(f"GEO 轨道周期 T = 2π√(r³/μ) = {T_GEO:.0f} s ≈ {fmt_hms(T_GEO)}（约一个恒星日；35786 km 为取整值，周期会有一两秒之差）")

    print("-" * 52)
    print("霍曼转移总 Δv（从 GEO 抬升到处置轨道，两次点火之和）:")
    raises = [235, 300, 400]
    for rk in raises:
        r2 = r_GEO + rk
        dv = hohmann_dv_mps(r_GEO, r2, mu)
        print(f"  +{rk:>3} km → Δv = {dv:.2f} m/s（处置轨道近地点高度 {r2 - R_E:.0f} km）")

    print("-" * 52)
    iadc_raise = 235 + 1000 * C_R * A_m
    print(f"IADC 公式验算：235 + 1000×C_R×A/m = 235 + {1000 * C_R * A_m:.0f} = {iadc_raise:.0f} km")
    print(f"  （C_R={C_R}、A/m={A_m} m²/kg，即典型通信卫星取值）")
    print(f"  处置轨道近地点高度 = {r_GEO + iadc_raise - R_E:.0f} km（在 GEO 之上 {iadc_raise:.0f} km）")
    print(f"  GSAT-12 实际最低近地点抬升 {ex['min_perigee_raise_km']} km，与上式 {iadc_raise:.0f} km 差 {abs(iadc_raise - ex['min_perigee_raise_km']):.0f} km")

    print("-" * 52)
    print(f"工作实例 {ex['date']}: {ex['title']}（{ex['tier']}）")
    print(f"  {ex['note']}；最终轨道：{ex['final_orbit_note']}；偏心度要求 ≤ {ex['ecc_requirement']}，共 {ex['burns']} 次点火")
    print(f"罚单 {en['date']}: {en['title']}（{en['tier']}）")
    print(f"  {en['note']}；罚款 ${en['penalty_usd']:,}")

    print("-" * 52)
    print("验证层级汇总（可复现≠可验证，quote 为外部正确性逐字证据）:")
    print("  primary（主源逐字）: GEO 保护区定义（NTRS）+ IADC 公式（A/AC.105/820）+ ITU 300km（A/AC.105/820）+ 低轨 25 年（A/AC.105/820）")
    print("  secondary（官方二手）: GSAT-12 处置（ISRO 新闻稿）+ EchoStar-7 罚单（FCC Consent Decree）")
    print(f"快照格式: {snap.get('format', '?')} · 锚点行数 {snap.get('rows', '?')}")
    print(f"快照 sha256: {sha}")
    print("验证方式：重跑本脚本，与文章正文表格与派生数字逐项核对（Δv 两次点火之和、近地点高度、IADC 公式）。")
    return 0


if __name__ == "__main__":
    sys.exit(main())

graveyard_orbit.json — 数据快照全文

{
  "format": "chenji-graveyard-orbit-v2",
  "schema": {
    "version": 2,
    "evidence_tiers": ["primary", "secondary", "media"],
    "rule": "primary 档锚点必须带原文引用句(quote)+出处定位(locator)；secondary 为归因清楚的官方二手口径；media 为媒体汇编。可复现≠可验证：quote 是外部正确性的逐字证据，让 AI agent 无需依赖权威即可核验。"
  },
  "generated": "2026-08-16",
  "source": "联合国和平利用外层空间委员会 UNOOSA A/AC.105/820（2003-12-15）+ NASA NTRS 20100007939（Nicholas L. Johnson, IAC-10-A6.4.1）公开文献原文",
  "rows": 4,
  "constants": {
    "mu_km3_s2": 398600.4418,
    "R_E_km": 6378.137,
    "GEO_alt_km": 35786,
    "C_R": 1.4,
    "A_over_m_m2_kg": 0.02,
    "note": "μ（地球引力常数）与 R_E（赤道半径）为 WGS84 标准值；GEO 标称高度 35786 km 为国际惯用值。"
  },
  "anchors": [
    {
      "asof": "2010",
      "label": "GEO 保护区域定义",
      "tier": "primary",
      "quote": "geosynchronous orbit (GEO), which includes the volume of space encompassed by 35,786 km +/- 200 km in altitude and +/- 15 degrees in latitude",
      "locator": "NTRS 20100007939《Medium Earth Orbits: Is There a Need for a Third Protected Region?》引言第 1 段（第 1 页）",
      "ref": "NASA NTRS 20100007939，Nicholas L. Johnson，IAC-10-A6.4.1"
    },
    {
      "asof": "2003",
      "label": "IADC 墓地轨道抬升公式（GEO 处置）",
      "tier": "primary",
      "quote": "235 km + (1,000.CR.A/m)",
      "locator": "UNOOSA A/AC.105/820 第 2 页（印度意见第 1 段，引 IADC 空间碎片减缓指南建议 5.3.1）",
      "ref": "联合国 UNOOSA 文件 A/AC.105/820（2003-12-15），IADC 空间碎片减缓建议：会员国意见汇编"
    },
    {
      "asof": "2003",
      "label": "ITU 建议最低抬升 300 km",
      "tier": "primary",
      "quote": "the satellite orbit should be raised to no less than 300 km fixed altitude above the geostationary orbit",
      "locator": "UNOOSA A/AC.105/820 第 2 页（印度意见第 3 段，引 ITU-R 建议 S.1003《对地静止轨道环境保护》）",
      "ref": "联合国 UNOOSA 文件 A/AC.105/820 + 国际电信联盟 ITU-R S.1003"
    },
    {
      "asof": "2003",
      "label": "低轨 25 年离轨规则",
      "tier": "primary",
      "quote": "a space system at the end of its mission should be left in low-Earth orbit so that the atmospheric drag would restrict the orbital lifetime to a limited number of years. Twenty-five years is recommended as a reasonable and appropriate lifetime limit for such disposal.",
      "locator": "UNOOSA A/AC.105/820 第 2-3 页（印度意见第 4-5 段，引 IADC 指南建议 5.3.2）",
      "ref": "联合国 UNOOSA 文件 A/AC.105/820（2003-12-15）"
    }
  ],
  "worked_example": {
    "title": "GSAT-12 后任务处置（做对了的样本）",
    "date": "2023-03-23",
    "agency": "ISRO",
    "launch_date": "2011-07-15",
    "min_perigee_raise_km": 261,
    "ecc_requirement": 0.003,
    "burns": 7,
    "final_orbit_note": "近 400 km 超同步圆轨道（GEO 标称高度 35786 km 之上）",
    "tier": "secondary",
    "locator": "ISRO 官网新闻稿《Post Mission Disposal of GSAT-12》（2023-04）",
    "note": "最低近地点抬升 261 km——与 IADC 公式按 C_R=1.4、A/m=0.02 m²/kg 算得的 263 km 几乎一致；PMD 完全符合 IADC 与联合国减缓指南"
  },
  "enforcement": {
    "title": "EchoStar-7 处置违规——全球首张太空碎片罚单",
    "agency": "FCC（美国联邦通信委员会）",
    "penalty_usd": 150000,
    "date": "2023-10-02",
    "retired": "2022-05",
    "required_raise_km": 300,
    "actual_raise_km": 122,
    "shortfall_km": 178,
    "tier": "secondary",
    "locator": "FCC 执法局 Consent Decree DA-23-888A1（2023-10-02）",
    "note": "DISH Operating 同意和解并承认责任；因推进剂余量不足，EchoStar-7 只抬升到距对地静止弧 122 km，比 2012 年获批计划要求的 300 km 低 178 km"
  },
  "background": {
    "geo_region_note": {
      "tier": "primary",
      "quote": "the GEO protected region is only a small portion of the GEO region, as defined previously, since the protected region only includes the region between 15 degrees north latitude and 15 degrees south latitude",
      "locator": "NTRS 20100007939 第 III 节（第 2-3 页）",
      "note": "GEO 保护区只是 GEO 轨道带（35786±200 km）与南北纬 15° 围成的一小片"
    },
    "graveyard_rationale": {
      "tier": "primary",
      "quote": "GEO space structures be maneuvered to long-term disposal orbits sufficiently far above the GEO region",
      "locator": "NTRS 20100007939 第 III 节（第 3 页）",
      "note": "IADC 指南建议：GEO 航天器寿命终了时机动到充分高于 GEO 区域的长期处置轨道"
    }
  }
}

数据来源：联合国 UNOOSA 文件 A/AC.105/820（2003-12-15）—— IADC 空间碎片减缓建议（抬升公式 235+1000·C_R·A/m、ITU-R S.1003 不低于 300 km、低轨 25 年规则）；NASA NTRS 20100007939（Nicholas L. Johnson）—— GEO 保护区定义（35786±200 km、±15° 纬度）。
历史背景：ISRO 官网新闻稿 —— GSAT-12 后任务处置（2023-03-23，7 次点火，最低近地点抬升 261 km，最终约 400 km 超同步轨道）；FCC 执法局 Consent Decree DA-23-888A1（2023-10-02）—— DISH/EchoStar-7 处置违规，罚款 $150,000。