基本上看看这个公式就明白了。

$\frac{\partial \left(\phi\rho s_i \right )}{\partial t} =\nabla\left[\phi\rho s_i\frac{\kappa_i}{\mu_i} \left(\nabla p - \rho\vec{g} \right )\right ] +\nabla\left[ D_i \nabla\left (\phi\rho s_i \right )\right] +r_i$

这是质量守恒决定的油气连续性方程，目的是求取岩石孔隙流体的密度 $\rho$ （或者压力 $p$ ）和任一流体成分 $i$ 的饱和度 $s$ 随时间空间的变化。其它几乎每个参数背后的都由另外一些方程决定。

我们不可能确切知道压实的参数，所以不可能精确模拟孔隙度 $\phi$ ;

由于无法准确了解地球各个圈层的热生成和热传递，不可能精确模拟温度对时空的变化；

而其它参数（渗透率 $\kappa$ 、粘度 $\mu$ 、扩散系数 $D$ 、反应速度 $r$ ）不但依赖于空间，而且依赖于温度，所以更不可能准确求得；

除此之外，反应速度依赖于基元化学反应及其动力学参数，这个不可能准确知道；

所以流体成分不可能准确算出；

而渗透率、粘度和扩散系数又依赖于成分；

初始条件很粗略；

边界条件更任意，不可能准确设定边界压力（这个涉及到毛细管力，又是空间不均匀而且依赖于成分和温度），不可能准确设定破裂压力；

所有这些不确定性决定了，虽然模型可以大致算出油气的量，但是不可能算出一公里范围内某一深度的油气饱和度s和压力p，更不要说每种成分的含量了。

石耀霖院士1990年代做过这方面工作，他完全了解正向模拟的局限，所以他应该斩钉截铁地说“地震的动力学预报根本不可行”。

大气圈的模拟只有比这个更复杂，想一想就知道。所以正向模拟不可能得到二氧化碳的气候敏感性，这个只能用足够多的实际数据来校正——也就是“马后炮”。