文章描述了用于解决质量平衡方程（2.25.1）的Crank-Nicholson解法。方程的完全显式分解在（2.25.15）中给出，用于求解水位线上方和下方的气体和液相浓度。

Crank-Nicholson解法结合了完全显式和隐式表示，形成了方程（2.25.16）。该方程通过标准的三对角线解法求解，如（2.25.17）所示。

文章还提到，在每个时间步长上执行两个甲烷平衡检查，以确保严格遵守甲烷分子的守恒（除了产生和消耗之外）和碳原子。