GRB131231A是EP和双追踪的特别风暴。为了进一步分析这场风暴的能谱特征，我们使用贝叶斯方法对这场风暴进行了详细的时间分辨光谱分析，发现24个分辨光谱中有16个具有明显的热成分，主要分布在脉冲的早期和峰值附近，支持伽玛暴喷流从物质主导型向磁性主导型的转变。加入热分量后，参数EP和变得坚硬，而峰值流量Fp几乎没有变化。利用光球层成分计算相对论污水的特征参数后，发现污水的特征参数，如初始半径r0、饱和层rs半径、光球半径rph和污水滑行过程中的洛伦兹系数，都发生了明显的变化。通过转动幂律模型也可以拟合其温度随时间的演变，上升和下降阶段的幂律指数分别为a=7.653.48和a=7.653.48，B=-2.103.48，下降指数B接近-2/3，与之前的结果有较大偏差。