--- title: "从一份GOES文件重算耀斑曲线" description: "读取固定NOAA原始文件，解码UTC与质量标志，重算两个X射线波段的峰值和条件通量积分。" canonical_url: "https://sandboxphysics.com/zh/stories/solar-magnetic-reconnection/one-goes-file-one-flare/" language: "zh-CN" page_type: "story" alternate_languages: en: "https://sandboxphysics.com/stories/solar-magnetic-reconnection/one-goes-file-one-flare/" last_modified: "2026-11-03" --- [从磁应力走向磁重联](https://sandboxphysics.com/zh/stories/solar-magnetic-reconnection/) L3 · 磁重联与太阳耀斑信号 # 从一份GOES文件重算耀斑曲线 Sandbox Physics · 2026-11-03 固定的一份NOAA文件，就能重算2024年5月10日事件的两个X射线波段峰值，以及定义清楚的通量积分。先认原始字节，再读时间与质量字段，最后才选耀斑区间；计算从头到尾使用观测记录。 ## 开始前需要什么 需要会用Python数组、文件哈希，并分清辐照度与能量。离线分析需要Python3.12、NumPy和h5py；科学包给出固定依赖，画图另需Matplotlib。 ## 先确认是哪一份原始文件 选用GOES16科学质量XRS一分钟产品，日期为2024年5月10日，版本为v2-2-2。这次重处理于2026年10月发布。包中保留原始NetCDF4/HDF5文件，共301570字节，元数据也保持原样。 读取数值前，程序同时检查长度与完整SHA-256。文件名也是身份的一部分。名称相近的业务产品或旧版本，校准和处理可能不同，不能因为曲线相似就混用。 辅助程序的计划与验证模式都不联网。明确要求重新获取时，它只请求这个具名URL，并拒绝改变或过大的输入。所选文件允许自由使用和再分发，包内保留NOAA NCEI太阳辐照团队署名。 ``` python scripts/astronomy/ac28/reacquire_goes.py --plan python scripts/astronomy/ac28/reacquire_goes.py --verify docs/astronomy/ac28/research/sci_xrsf-l2-avg1m_g16_d20240510_v2-2-2.nc ``` | 输入身份 | 固定值 | | --- | --- | | 字节 | 301570 | | SHA-256 | 93146eb0a920a8c0cfe5f74897c90978e17006f8280dcc73a5f99aeafec8ba16 | | 产品 | sci\_xrsf-l2-avg1m\_g16\_d20240510\_v2-2-2.nc | 来源：[NOAA NCEI, GOES16 XRS science-quality one-minute means, 2024-05-10, v2-2-2](https://data.ngdc.noaa.gov/platforms/solar-space-observing-satellites/goes/goes16/l2/data/xrsf-l2-avg1m_science/2024/05/sci_xrsf-l2-avg1m_g16_d20240510_v2-2-2.nc) · [NOAA, GOES-R XRS L2 Data ReadMe, release v2-2-2, 6 October 2026](https://data.ngdc.noaa.gov/platforms/solar-space-observing-satellites/goes/goes16/l2/docs/GOES-R_XRS_L2_Data_ReadMe.pdf) · [NOAA NCEI Open Data Policy, signed 2023, §4a; selected file’s unrestricted-reuse statement](https://www.ncei.noaa.gov/sites/default/files/2023-12/NCEI%20PD-10-2-02%20-%20Open%20Data%20Policy%20Signed.pdf) ## 读清波段和真正的时钟 短波段XRS-A覆盖 $0.05\text{–}0.4\,\mathrm{nm}$，XRS-B覆盖 $0.1\text{–}0.8\,\mathrm{nm}$。两列都给出观测者处的波段辐照度，单位为 $\mathrm{W\,m^{-2}}$。 文件有1440个按分钟起点标记的均值。时间坐标从2000年1月1日12:00:00 UTC起计秒，忽略闰秒，不是天文学里以TT定义的J2000时刻。程序用两种独立日历转换相互比较，保留源文件的约定。 写成06:54的峰，指从06:54开始的那一分钟。比较峰值或更快的变化时，要一直保留这个一分钟平均的时间分辨率。 来源：[NOAA, GOES-R XRS L2 Data Users Guide, §§2.1–2.4 and quality-flag tables](https://data.ngdc.noaa.gov/platforms/solar-space-observing-satellites/goes/goes16/l2/docs/GOES-R_XRS_L2_Data_Users_Guide.pdf) · [NOAA NCEI, GOES16 XRS science-quality one-minute means, 2024-05-10, v2-2-2](https://data.ngdc.noaa.gov/platforms/solar-space-observing-satellites/goes/goes16/l2/data/xrsf-l2-avg1m_science/2024/05/sci_xrsf-l2-avg1m_g16_d20240510_v2-2-2.nc) ## 先筛分钟，再数实际参与的秒 读取主质量标志、参与的一秒样本数和电子修正标志。这里要求主标志为零，有效样本在54到60之间，辐照度有限且高于源下限；电子修正还须有效，没有插值、衰减或显著主波段修正。计数字段的填充值255要拒绝。 五十四个样本是本分析选择的九成覆盖门槛，不会补出缺失的秒。主产品标记为插值的分钟也保留缺口，不能仅因画折线方便而连过去。 06:00至08:00 UTC内，两个波段各有120个分钟通过筛选。但A波段四个分钟没有采满，实际7195秒，最少的一个分钟只有58秒；B波段则有完整7200个一秒样本。分钟通过率与完整曝光，是不同的量。 | 所选事件 | XRS-A | XRS-B | | --- | --- | --- | | 通过的分钟数 | 120 | 120 | | 参与的秒数 | 7195 | 7200 | | 未采满的分钟数 | 4 | 0 | | 单分钟最少秒数 | 58 | 60 | 来源：[NOAA, GOES-R XRS L2 Data Users Guide, §§2.1–2.4 and quality-flag tables](https://data.ngdc.noaa.gov/platforms/solar-space-observing-satellites/goes/goes16/l2/docs/GOES-R_XRS_L2_Data_Users_Guide.pdf) · [NOAA NCEI, GOES16 XRS science-quality one-minute means, 2024-05-10, v2-2-2](https://data.ngdc.noaa.gov/platforms/solar-space-observing-satellites/goes/goes16/l2/data/xrsf-l2-avg1m_science/2024/05/sci_xrsf-l2-avg1m_g16_d20240510_v2-2-2.nc) ## 重算两个峰值 运行下面的分析命令，会写出结果JSON和原始元数据报告。画图是另外一条可选命令，不装绘图库也能检查数值。 A波段峰所在的分钟从06:52 UTC开始，均值为 $1.23114\times10^{-4}\,\mathrm{W\,m^{-2}}$；B波段从06:54开始，均值为 $3.92869\times10^{-4}\,\mathrm{W\,m^{-2}}$。把后者除以名义X级阈值，可得约3.93。 可以把它叫作科学产品中约X3.93级的一分钟均值。业务耀斑目录可能采用不同产品、时间分辨率与校准。两个所选分钟相差两分钟，也不足以直接量出精密的气体冷却时间。 ``` python -W error scripts/astronomy/ac28/analyze_goes.py python -W error scripts/astronomy/ac28/make_figures.py ``` ![GOES16实测全天波段辐照度和选中的两小时事件区间。](https://sandboxphysics.com/stories/solar-magnetic-reconnection/materials/goes16-two-band-flare.png) 两个波段来自同一原始产品。质量筛选拒绝的分钟保留缺口，均值仍用源文件的分钟起点时标。 来源：[NOAA NCEI, GOES16 XRS science-quality one-minute means, 2024-05-10, v2-2-2](https://data.ngdc.noaa.gov/platforms/solar-space-observing-satellites/goes/goes16/l2/data/xrsf-l2-avg1m_science/2024/05/sci_xrsf-l2-avg1m_g16_d20240510_v2-2-2.nc) · [NOAA, GOES-R XRS L2 Data Users Guide, §§2.1–2.4 and quality-flag tables](https://data.ngdc.noaa.gov/platforms/solar-space-observing-satellites/goes/goes16/l2/docs/GOES-R_XRS_L2_Data_Users_Guide.pdf) · [NOAA, GOES-R XRS L2 Data ReadMe, release v2-2-2, 6 October 2026](https://data.ngdc.noaa.gov/platforms/solar-space-observing-satellites/goes/goes16/l2/docs/GOES-R_XRS_L2_Data_ReadMe.pdf) ## 说清积分到底估计什么 辐照度乘秒，得到单位面积上的能量。把分钟均值逐一乘六十再相加，估计保留区间的整分钟通量积分。若六十个样本齐全，使用的就是完整分钟均值；只有54到59个样本时，则假设现有秒样本能代表整分钟。 A波段的整分钟估计为 $0.14607970\,\mathrm{J\,m^{-2}}$，按实际样本数加权则得另一个代理量 $0.14595377\,\mathrm{J\,m^{-2}}$。两者之差不是缺失五秒的误差上界。B波段整分钟估计为 $0.65148857\,\mathrm{J\,m^{-2}}$。 再指定背景区间：05:30至05:50，或05:40至06:00，各用均值与中位数。这四种选择给出B波段带符号超额通量约每平方米0.606到0.610焦耳。程序另列只积正值的结果，因为截去负值会改变统计量。 源文件的au\_factor还给出单独命名的一天文单位修正结果，原生辐照度并未预先乘上它。这里也没有套用较早GOES8至15系列的历史比例因子。四种背景与这些修正选择，都属于敏感性比较。 $$ \widehat{\mathcal F}_{60}=\sum_{i\in\mathrm{good}}\overline F_i(60\,\mathrm{s}),\qquad\mathcal F_{\mathrm{available}}=\sum_{i\in\mathrm{good}}\overline F_i(n_i\,\mathrm{s}) $$ 未采满的分钟里，第一项把现有均值外推到整分钟；第二项按可用秒数加权。被拒绝的分钟没有补值。 $$ \widehat{\mathcal F}_{\mathrm{excess}}=\sum_{i\in\mathrm{good}}(\overline F_i-F_{\mathrm{background}})(60\,\mathrm{s}) $$ 带符号超额量允许单个分钟为负；将负值截成零，就是另一个估计。 来源：[NOAA, GOES-R XRS L2 Data Users Guide, §§2.1–2.4 and quality-flag tables](https://data.ngdc.noaa.gov/platforms/solar-space-observing-satellites/goes/goes16/l2/docs/GOES-R_XRS_L2_Data_Users_Guide.pdf) · [NOAA NCEI, GOES16 XRS science-quality one-minute means, 2024-05-10, v2-2-2](https://data.ngdc.noaa.gov/platforms/solar-space-observing-satellites/goes/goes16/l2/data/xrsf-l2-avg1m_science/2024/05/sci_xrsf-l2-avg1m_g16_d20240510_v2-2-2.nc) ## 让定义跟着物理结论走 得到的积分属于这些波段、这个观察方向、选定区间与分钟代表性假设。它不是整场耀斑的总辐射能，更不是释放的磁能。这里也没有用两个波段反演经过校准的温度和发射量。 若要与MHD耗散热定量比较，还需补上密度、热力学、热传导、辐射与仪器响应。眼下，一边提供磁场与流体的能量转换，另一边提供真实X射线信号，各自都能沿原始输入重算到结果。 换下一份任务数据时，仍按这个次序读：输入与版本、时钟、波段、质量标志、实际样本、估计方法、背景和物理解释。引用数值时，把这些选择一起带走。 来源：[NOAA, GOES-R XRS L2 Data Users Guide, §§2.1–2.4 and quality-flag tables](https://data.ngdc.noaa.gov/platforms/solar-space-observing-satellites/goes/goes16/l2/docs/GOES-R_XRS_L2_Data_Users_Guide.pdf) · [NOAA, GOES-R XRS L2 Data ReadMe, release v2-2-2, 6 October 2026](https://data.ngdc.noaa.gov/platforms/solar-space-observing-satellites/goes/goes16/l2/docs/GOES-R_XRS_L2_Data_ReadMe.pdf) ## 独立重跑科学计算 小型材料包含独立Python、精确NOAA输入、自算MHD状态，以及来源、单位和失败说明。固定依赖与离线命令都在README中。 - [solar-reconnection-science.zip](https://sandboxphysics.com/stories/solar-magnetic-reconnection/materials/solar-reconnection-science.zip) 6.33 MB - [include-manifest.json](https://sandboxphysics.com/stories/solar-magnetic-reconnection/materials/include-manifest.json) 0.01 MB - [README-science.md](https://sandboxphysics.com/stories/solar-magnetic-reconnection/materials/README-science.md) 0.01 MB - [runtime-domain-reference.json](https://sandboxphysics.com/stories/solar-magnetic-reconnection/materials/runtime-domain-reference.json) 0.21 MB - [goes-reference-results.json](https://sandboxphysics.com/stories/solar-magnetic-reconnection/materials/goes-reference-results.json) 0.01 MB - [retained-counterexamples.json](https://sandboxphysics.com/stories/solar-magnetic-reconnection/materials/retained-counterexamples.json) 0.00 MB 精确文件与使用许可 solar-reconnection-science.zip 6325678 bytes SHA-256: 4f680fc833816eb4ae27fd1bce178ee61af985c0f800097e86e11cd489710f07 include-manifest.json 10996 bytes SHA-256: 363ff93f8e96ae524583d067b560f865ff526021cb45bc725a21efd41161e405 README-science.md 9685 bytes SHA-256: f0aee6da22e4e5ff73aa326ee9b8b3028601bed9b1e2a9974056588db417d32a runtime-domain-reference.json 208128 bytes SHA-256: 33cdb0c1deec2889759846f572b8d32ce1da3c2de9278e549a068bece7391f5d goes-reference-results.json 5924 bytes SHA-256: cc4e4785d66571658b1824d9c4aa5fdacfdba5f0643c2f126245836771a5c197 retained-counterexamples.json 1826 bytes SHA-256: 1e89efaf39a01b431ea1aeb75c8b9a108f508a4ab1972351e47df68266db561e 所选NOAA原始输入允许自由使用与再分发。新写独立Python采用MIT，自算数据、图与文档采用CC BY4.0并保留来源署名。原论文PDF和网站Engine代码不在包中。 ## 自己推一步 一个分钟只有58个一秒样本。其均值乘六十，是否就是完整实测的六十秒积分？ 它是在现有样本能代表缺失秒数的假设下，估计整分钟积分。乘五十八给出按可用秒数加权的代理量，两种运算都无法告诉我们缺失两秒发生了什么。 为什么四种背景结果不构成一个置信区间？ 它们只是四种指定分析情景，没有给出背景、校准或缺失样本的概率分布，因此这个范围没有置信水平的含义。 能否把MHD热曲线与GOES曲线对到最好看，从而测出磁扩散率？ 模型尚无从无量纲热量到这些波段的密度、温度、辐射及仪器映射。按形状对齐会留下未定义的辐射与仪器转换。识别磁扩散率，需要先建立模型状态到实测波段的物理映射。 ### [辨认从耗散到发光还缺什么](https://sandboxphysics.com/zh/labs/magnetic-reconnection/) 1. 暂停一个电流片状态，写清模型实际记录哪一种热量。 2. 列出预测某个X射线波段还需要的物理场。 3. 展示比较时，给真实GOES图与合成耗散热图各自命名。 ## 原始来源 1. [NOAA NCEI, GOES16 XRS science-quality one-minute means, 2024-05-10, v2-2-2](https://data.ngdc.noaa.gov/platforms/solar-space-observing-satellites/goes/goes16/l2/data/xrsf-l2-avg1m_science/2024/05/sci_xrsf-l2-avg1m_g16_d20240510_v2-2-2.nc) 2. [NOAA, GOES-R XRS L2 Data Users Guide, §§2.1–2.4 and quality-flag tables](https://data.ngdc.noaa.gov/platforms/solar-space-observing-satellites/goes/goes16/l2/docs/GOES-R_XRS_L2_Data_Users_Guide.pdf) 3. [NOAA, GOES-R XRS L2 Data ReadMe, release v2-2-2, 6 October 2026](https://data.ngdc.noaa.gov/platforms/solar-space-observing-satellites/goes/goes16/l2/docs/GOES-R_XRS_L2_Data_ReadMe.pdf) 4. [NOAA NCEI Open Data Policy, signed 2023, §4a; selected file’s unrestricted-reuse statement](https://www.ncei.noaa.gov/sites/default/files/2023-12/NCEI%20PD-10-2-02%20-%20Open%20Data%20Policy%20Signed.pdf) [上一篇：沿着重联计算追踪能量](https://sandboxphysics.com/zh/stories/solar-magnetic-reconnection/follow-the-converted-energy/) [下一篇：重联的流体图景在哪里不再够用](https://sandboxphysics.com/zh/stories/solar-magnetic-reconnection/where-the-fluid-picture-runs-out/)