上机实验

已知在特定温度下,煤层中吸附瓦斯的量可用朗格缪尔吸附等温线方程式表示:

$$V_a = \frac{{abp}}{{1 + bp}}$$

式中:

  • \(V_a\) ——单位质量的煤在气体压力为 \(p\) 时吸附瓦斯量在标准状态下的体积,m3
  • \(p\)——气体压力,MPa;
  • \(a\)——单位质量的煤的表面覆盖满单分子层时的吸附量,也称最大吸附量,通常称为朗格缪尔体积,m3/t;
  • \(b\)——吸附系数,MPa。

已知,从某煤矿的一个煤层中取样,在实验室中进行等温吸附实验,得到瓦斯压力 \(p\) 和含量 \(V_a\) 可定义如下:

gas_pressure = [0.05, 0.2, 0.4, 0.6, 0.8, 1, 1.2, 1.4, 1.6, 1.8, 2]  # 单位为 MPa
gas_content = [1.12, 4, 7.53, 10.1, 13, 14.3, 15.9, 17.6, 20.1, 20.4, 22.1]  # 单位为 m**3

请利用 SciPy,用回归分析方法计算吸附常数 \(a\) 和 \(b\)。

预测问题数学建模安全评价方法及其实现