|||
基本上看看这个公式就明白了。
$\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年代做过这方面工作,他完全了解正向模拟的局限,所以他应该斩钉截铁地说“地震的动力学预报根本不可行”。
大气圈的模拟只有比这个更复杂,想一想就知道。所以正向模拟不可能得到二氧化碳的气候敏感性,这个只能用足够多的实际数据来校正——也就是“马后炮”。
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-9-27 12:00
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社