P27
第 25 太阳活动周耀斑等待时间分布的独立检验：非平稳泊松模型与幂律尾指数
来源说明：本则出自独立撰写的第二批方案。它与前 20 则同样遵循八块结构与硬门槛要求，但撰写时未做英文文献检索，新颖性边界依据的是历届获奖图谱与既有知识，而非当轮查新。因此其「需核实」条目更多，启动前须自行补一轮英文检索。
1 · 研究问题
第 25 太阳周(2020–2026)GOES 软 X 射线耀斑目录的等待时间分布，其幂律尾指数与已发表的早期周期结果(Wheatland 2000 等报告的约 −2.1 至 −2.4 区间，精确值需核实)是否在 95% 置信区间内一致？偏离能否由非平稳泊松模型(活动率随周期演化)完全解释？
2 · 研究背景与空白
太阳耀斑是磁能突然释放事件,GOES 卫星自 1970 年代起连续记录全日面软 X 射线流量，耀斑目录(时间、等级)是太阳物理最长的均匀事件序列之一。相邻耀斑的"等待时间"分布形态是区分触发机制的探针：若耀斑为速率缓变的泊松过程，分布应为活动率分布的混合指数，其长时尾呈幂律——这是 Wheatland(2000, ApJ)非平稳泊松图像的核心预言；自组织临界(SOC)框架(Aschwanden 等)也对该分布有预期。
已有工作覆盖第 21–23 周，报告长时尾幂律并支持非平稳泊松解释；第 24 周有零星分析。
空白在于时间覆盖：第 25 周(2019 年底起,2024–2025 达极大)刚刚完整，尚无(或极少，需核实：搜索关键词 "solar flare waiting time distribution cycle 25")用同一口径对其完整重做的公开分析。方法成熟、数据公开、指数漂移与否两个结果都有物理含义，是教科书级的"新样本独立检验"型学生课题。
3 · 可检验假设
- H1：第 25 周 M 级及以上耀斑等待时间分布长时尾(> 1 天)幂律指数与本项目在第 23 周复现值的差 ≤ 0.3(置信区间重叠)。
- H2：用当周耀斑发生率曲线驱动的非平稳泊松蒙特卡洛，能在 KS 检验 p > 0.05 水平上再现观测分布；若被拒绝，超出部分(如短时聚集)即为主要结论(同震余震式聚集的证据，标注为待解释现象)。
4 · 量化验收标准
- 方法学校验(硬门槛)：用与已发表工作相同的历史时段(如第 22–23 周)重跑等待时间分布，幂律尾指数与文献值偏差 ≤ 10%，分布形态图与文献图目视逐段一致。此条不过关，后续全部结论无效。
- 统计口径预先写死：幂律尾用最大似然(CSN 方法)估计 + 自助法 95% 置信区间;x_min 由 KS 距离定；不用双对数直线拟合作主结果。事件阈值敏感性:C/M/X 三档阈值分别报告全套结果。
- 数据完整性预算：仪器换代(GOES-8→16 系列)与目录缺口时段列表化，剔除缺口 > 1 天的区间并报告剔除比例(< 5%)。
- H2 蒙特卡洛 ≥ 10⁴ 次合成序列,KS 统计量的零分布由合成序列自举。
- 样本量下限：第 25 周 M 级以上事件数 ≥ 500(预计满足，准确数第 2 周核实)。
- 可复现性：下载脚本 + 清洗规则 + 分析脚本全部开源，从原始目录到全部图表一键重跑。
5 · 数据与工具
| 用途 | 来源 / 工具 |
|---|---|
| 耀斑目录 | NOAA SWPC / NCEI 的 GOES X 射线耀斑事件表(公开,swpc.noaa.gov 与 ngdc.noaa.gov 事件列表；具体文件路径以当前官网为准，需核实)；备选：通过 SunPy 的 Fido/HEK 接口检索耀斑事件(sunpy.org，开源) |
| 覆盖范围 | 1975 至今全日面记录；本项目用 1996–2026(第 23–25 周)；表格级数据 < 100 MB，纯 CPU |
| 分析 | Python:pandas + numpy + powerlaw 包 + 自写蒙特卡洛；全套分析笔记本级 CPU < 1 小时 |
| 理论对照 | Wheatland 2000 非平稳泊松框架与早期周期指数,仅用于校验与对比，不计入本项目的数据贡献 |
6 · 方法路径
- 装 SunPy/pandas 环境，下载并解析全部目录；列出仪器换代与缺口清单。
- 在第 22–23 周历史段复现已发表分布与指数(验收 1)。
- 固化清洗与阈值规则，生成第 25 周等待时间序列。
- CSN 拟合 + 自助法，三阈值敏感性矩阵。
- 构建非平稳泊松蒙特卡洛(速率 = 平滑月发生率)，做 KS 检验(H2)。
- 核实边界:GOES 换代前后目录定标是否一致(读 NCEI 说明文档；不一致则分段分析)。
- 独立交叉校验：用秩序统计(去趋势后的归一化等待时间)复核幂律尾稳健性。
7 · 新颖性边界
- 本课题不声称提出等待时间统计新方法，不声称判定耀斑触发机制；非平稳泊松框架与早期指数为 Wheatland(2000)等已发表结论。
- 已有工作：第 21–23 周分布与指数、SOC 解释(Aschwanden 等綜述)。
- 本项目贡献(主结论)：第 25 周完整数据的同口径独立检验——指数漂移则提示分布依赖周期强度(第 25 周弱于 23 周)，不漂移则普适性获得一次独立支持。两种结果都构成有效结论。属"时间覆盖空白 + 新样本独立检验"。
- 核对获奖清单：本赛道尚无太阳统计课题；最接近的公开数据获奖题为 "Developing an Automated Pipeline for Identifying False Positives Among Released TESS Objects of Interest"(2021 铜)与 "Investigating Physical Conditions across the Center of M82"(2024 铜)，均非太阳物理，说明评委接受档案数据管线类工作。
8 · 决策门槛(go / no-go)
⚠️ 第 2 周末：核实数据入口与事件数。阈值：能取到 1996–2026 连续目录、第 25 周 M+ 事件 ≥ 500。若 NCEI 文件接口变动，改走 SunPy/HEK；若 M+ 不足 500,降级：阈值降至 C 级(事件数增一个数量级)，框架不变。 ⚠️ 第 6 周末：完成验收 1 历史复现，指数偏差 ≤ 10%；不达标则冻结新分析，排查清洗规则差异(常见原因：重复事件、阈值口径)，直到达标。 ⚠️ 第 20 周检查点：若第 25 周分布与非平稳泊松完全一致且与历史无差(最平淡结果)，按预案把"三阈值 + 跨周期指数表 + 完全可复现管线"作为主交付，结论定位为普适性支持——仍是完整论文。