物理化学综合反演

上一课从一个静态液柱反演界面参数。最终项目面对更复杂的可逆反应批次：平衡转化率、早期转化曲线、温升和搅拌依赖同时出现。任务是用最少参数解释这些观测，并用未参与拟合的三十分钟数据判断模型是否真的能预测。

在二百九十八开，一摩尔归一化初始 A 的平衡转化率为 $x_{\mathrm{eq}}=0.800$；十分钟测得 $x=0.400$。先假设传质足够快、温度在动力学阶段受控，宏观逼近平衡可写

$$\frac{dx}{dt}=k(x_{\rm eq}-x),\qquad
x(t)=x_{\rm eq}[1-e^{-kt}].$$

第一阶段，由十分钟点

$$k=-\frac1{10}\ln\left(1-\frac{0.400}{0.800}\right)
=0.06931\ {\rm min^{-1}}.$$

第二阶段，在三十分钟做不再调参的预测：

$$x(30)=0.800[1-e^{-0.06931(30)}]
=0.7000.$$

若留出实测为零点六九零，绝对残差零点零一零、相对平衡转化约百分之一点二五，当前模型可暂时接受。若不同搅拌速率下早期曲线改变，原先“传质足够快”的假设失效，应加入传质系数而非硬改 $k$。

第三阶段，三百一十八开拟合 $k=0.120\mathrm{mi}{\mathrm{n}}^{-1}$，两温度 Arrhenius 估计

$$E_a=R\frac{\ln(0.120/0.06931)}
{1/298-1/318}
=21.62\ {\rm kJ/mol}.$$

另一次近绝热终态中，总热容一点五千焦每开、温升一点六开、反应进度零点零八摩尔，给 $\Delta H\approx -30.0kJ/mol$。动力学活化能与反应焓是不同物理量，数值不能互相替代。

综合诊断按观测来源画残差：组成残差随时间、温度残差随反应进度、搅拌对照随转速。若单一参数只能改善一类数据却恶化另一类，说明结构缺项。留出不同温度和不同搅拌条件，比在同一曲线上增加更多点更能检验迁移能力。

同时检查参数相关矩阵，防止两参数互相补偿。

量化反馈标准：训练与留出转化率误差均低于百分之三；能量闭合低于百分之五；参数置信区间不跨越零且不同批次重叠；残差无时间或温度趋势；每项改进能对应某类残差变化。独立变式：保持 $x_{\mathrm{eq}}=0.800$ 与 $k=0.06931\mathrm{mi}{\mathrm{n}}^{-1}$，求半程到达 $x=0.400$ 的时间。答案恰为十分钟；到达平衡值百分之九十需 $t=\ln 10/k\approx 33.22$ 分钟。

本课完成物理化学项目闭环。下一轮应更换初始浓度与反应器尺度，让同一参数集预测新批次；若必须逐批重拟合，模型仍未抓住可迁移的化学机制。