---
title: "先拟合响应，再把冰壳—固体退化族算出来"
description: "下载完整Python实验：用合成复响应计算加权残差、壳厚条件剖面和不相容记录，比较液体壳层与全固态退化。"
canonical_url: "https://sandboxphysics.com/zh/stories/ocean-worlds/rerun-a-shell-and-solid-response-fit/"
language: "zh-CN"
page_type: "story"
alternate_languages:
  en: "https://sandboxphysics.com/stories/ocean-worlds/rerun-a-shell-and-solid-response-fit/"
last_modified: "2026-10-28"
---

[从材料记忆走到潮汐加热与海洋世界](https://sandboxphysics.com/zh/stories/ocean-worlds/)

AC12 · Tutorial

# 先拟合响应，再把冰壳—固体退化族算出来

Sandbox Physics · 2026-10-28

这份实验先拟合六个带噪合成响应分量，再构造内部截然不同、拟合却相同的一族天体。把材料解压到新目录，使用Python标准库就能运行。输入、噪声模型、期望输出与失败记录均公开，无须启动浏览器实验，也能亲自检查非唯一性。

## 开始前需要什么

上一章给出可识别性推导；会运行Python文件即可使用复算包，本页解释求解器与假设。

## 下载的是科学复算包，先认清记录来源

压缩包包含分析脚本、两个输入CSV、出处JSON和四份期望输出。主记录的探测频率是合成天体轨道频率的二分之一、一倍和两倍；每行列出重力Love数的实部、虚部及各自独立标准差。这些是本站生成的合成数据，不是Cassini跟踪残差，不是土卫六Love数时间序列，也不是从论文提取的数据集。

生成真值采用五十千米同密度壳、三点五吉帕弹性模量和十的十四点五次方帕秒黏度。高斯扰动使用固定种子12095，实部标准差零点零零二，虚部零点零零一六。保存数值是复算依据；分析时不重新抽噪声，也不暗中改变输入。

使用Python 3.11或更新版本，不需要第三方库。分发目录内有数据和期望输出子目录，脚本不导入网站代码。结果会写出输入的SHA-256哈希；遇到差异时，先与出处记录比较，再判断物理结果。原论文链接提供近似背景，公开包只分发原创教学材料。

对应阅读：Beuthe, Tides on Europa: the membrane paradigm, arXiv:1410.4735v1 · Petricca et al., Titan’s strong tidal dissipation precludes a subsurface ocean, Nature 648 (2025)

## 先运行，再逐个看残差

下载materials.zip并解压，在解压目录运行下面的命令。脚本会打印有效刚度、松弛时间、残差平方和以及条件剖面摘要。拟合分量、每项残差和条件剖面都保存为小型文本表格。把重跑结果写到另一个目录，方便保留原始期望输出。

打开fitted\_response.csv。每行用观测减预测，再除以列出的标准差；把六个标准化残差平方并相加，期望结果约为四点零零四零二九五。拟合使用两个有效参数，在指定独立高斯模型下留下四个残差自由度。自由度计数并不能验证结构假设本身。

几何刚度与弹性模量的拟合乘积约八点二五一七六亿帕，Maxwell时间约八万八千六百零二点五秒。取真实公转频率与零点零零五偏心率，这份拟合响应预测约一百点七四零吉瓦。它与无噪默认例子略有不同，因为拟合来自保存的带噪记录。

```
python3 analyze-response.py --input data --output rerun-output
```

对应阅读：Beuthe, Tides on Europa: the membrane paradigm, arXiv:1410.4735v1

## 看清求解器最小化的对象

两个拟合坐标分别是无量纲有效刚度和Maxwell时间的对数，指数还原后始终为正。脚本先在固定粗网格中找初始点，再用中心差分计算标准化残差的导数，解一个二乘二线性系统得到局部改正；如果平方和没有下降，就缩短步长。

正规方程来自一个局部步骤：在当前参数附近线性展开残差向量，把长度平方展开，对改正量求导并令其为零。它是数值优化工具，不是新加的物理定律。即使最终残差很小，仍可能有多个内部，因为它们映到完全相同的有效坐标。

真实行星分析可能包含相关无线电数据、伴随参数和系统误差，那时需要完整协方差，不能继续使用这里的对角权重。小实验把噪声模型明确写出，才能在不一致时分别追查记录、符号、物理模型或优化步骤。

$$
\boldsymbol r(\boldsymbol p+\Delta\boldsymbol p)\simeq\boldsymbol r+\boldsymbol J\Delta\boldsymbol p,\qquad (\boldsymbol J^{\mathsf T}\boldsymbol J)\Delta\boldsymbol p=-\boldsymbol J^{\mathsf T}\boldsymbol r
$$

每个残差均为观测减预测，再除以分量标准差；Jacobian对这个残差求导。

对应阅读：Beuthe, Tides on Europa: the membrane paradigm, arXiv:1410.4735v1

## 把替代内部逐行构造出来

打开shell\_solid\_family.csv，壳厚从十千米一直列到五百千米全固态。每个厚度先计算论文中的有限壳几何因子，再用拟合有效模量除以它，得到弹性模量；随后乘以拟合Maxwell时间得到黏度。外部预测在每个频率上保持不变。

带噪拟合的全固态行得到约零点二一七吉帕模量、十九点二四万亿帕秒黏度；五十千米行模量约三点五五吉帕。所有行的残差平方和相同。十千米替代体需要约十八点八吉帕，超出交互控件范围，也没有被声称为合理冰材料。独立材料信息的价值，正在这里出现。

再看conditional\_thickness\_profile.csv。严格固定三点五吉帕模量后，脚本在每个四分之一千米网格点对黏度最优化。最佳网格点为五十点五千米，局部剖面增量为一的区间为四十九点五至五十一点七五千米。读这个窄区间时，也要读它旁边的条件；放开模量后，恢复的是整个自由族，而不只是稍宽一点的误差棒。

对应阅读：Beuthe, Tides on Europa: the membrane paradigm, arXiv:1410.4735v1 · Petricca et al., Titan’s strong tidal dissipation precludes a subsurface ocean, Nature 648 (2025)

## 用失败记录检查符号与模型

第二份输入的虚部为正，并远超声明误差；采用正时间虚指数时，被动Maxwell天体无法产生它。脚本拒绝这份故意不相容的记录，并在summary.json写明原因。它不会裁掉响应、为了拟合翻转符号，也不会把它重新解释成炽热内部。

另一篇论文中的正虚数并不自动表示主动材料。作者可能使用相反时间指数，也可能把滞后大小报成正数。比较前应先翻译约定；2025年土卫六论文报告的是正耗散系数，导读会说明转换。这里的失败记录是在符号约定已经固定之后才不相容。

迁移练习可以先把两列误差都加倍再运行：自由族仍严格退化，有条件剖面则会放宽。随后只翻转虚部符号：被动模型应失败，而不是发现另一种海洋。这两次改动分别检验误差和相位究竟提供了什么信息。

对应阅读：Beuthe, Tides on Europa: the membrane paradigm, arXiv:1410.4735v1 · Petricca et al., Titan’s strong tidal dissipation precludes a subsurface ocean, Nature 648 (2025)

## [把推断带回物理场景](https://sandboxphysics.com/zh/labs/tidal-locking/)

1. 导出默认壳与重标度全固态预测，检查Love值，再观察可见内部边界的差别。
2. 把无噪预测与教程拟合输出比较，解释教程功率为什么接近默认读数却不完全相等。

## 练习与解释性答案

1. 模量与黏度都自由时，为什么再增加精确探测频率也不能唯一识别壳厚？

   展开解释

   所有频率依赖相同的两个有效组合，缩放变换在全频域保留它们。新增频率能检验Maxwell频率规律，却不能区分由这条严格对称变换相连的内部。

2. 哪些输出能区分带噪拟合记录与无噪实验默认值？

   展开解释

   拟合响应CSV有非零标准化残差，有效模量、时间和功率也与无噪值略有不同，出处标签不同；对每份选定拟合响应而言，其退化族仍然严格成立。

## 公开复算材料

- [下载完整实验压缩包](https://sandboxphysics.com/stories/ocean-worlds/materials/materials.zip)
- [阅读独立目录运行说明与范围](https://sandboxphysics.com/stories/ocean-worlds/materials/README.md)
- [Python 标准库分析脚本](https://sandboxphysics.com/stories/ocean-worlds/materials/analyze-response.py)
- [合成记录出处、误差与单位](https://sandboxphysics.com/stories/ocean-worlds/materials/provenance.json)
- [拟合响应与条件剖面的期望输出](https://sandboxphysics.com/stories/ocean-worlds/materials/summary.json)

[上一篇：一个复Love数，究竟能识别哪些内部信息](https://sandboxphysics.com/zh/stories/ocean-worlds/what-a-complex-love-number-can-identify/) [下一篇：沿土卫六潮汐研究读懂测量、内部与新实验](https://sandboxphysics.com/zh/stories/ocean-worlds/read-titan-dissipation-and-dragonfly-forecasts/)

## 来源与计算范围

1. [Beuthe, Tides on Europa: the membrane paradigm, arXiv:1410.4735v1](https://arxiv.org/html/1410.4735v1)

   固定为2014年10月17日预印本，Icarus正式发表DOI10.1016/j.icarus.2014.10.027。附录A式105–116给出有限厚度均匀不可压缩壳；附录C讨论Maxwell对应原理，附录E推导功与功率，附录G规定偏心率和光学天平动潮汐。图8–9、12–13展示近似与热模型依赖。核读相关正文和附录；本站的同密度Maxwell例子比真实卫星简单得多。

2. [Petricca et al., Titan’s strong tidal dissipation precludes a subsurface ocean, Nature 648 (2025)](https://www.nature.com/articles/s41586-025-09818-x)

   正式版2025年12月17日，DOI10.1038/s41586-025-09818-x，CC BY 4.0。核读出版方存档全文PMC12711566及Methods。图1给出重力响应估计，图2比较有条件的内部后验，图3是推断内部的示意图。图1误差棒为三倍标准差，图2椭圆为一倍标准差。原文正虚部表示其耗散约定；本站用正时间虚指数，被动Love响应的虚部为负。