数值模拟中网格方向的确定

在达西定律的基础上，通过数学分析得出结论：各向异性储集层在各个方向的渗透率数值虽不相等，但存在某种相关性，三维和二维渗透率张量分别是三阶和二阶对称矩阵，因而模拟模型网格设计的关键是求解主渗透率的方向余弦，以主渗透率方向作为模型网格的坐标轴方向。对于二维模拟来说，可以由４口井（Ａ，Ｂ，Ｃ，Ｄ）组成的典型井组为基本单元，先根据地质研究确定ｘ方向的主渗透率（Ｋｘ），再求解Ａ井、Ｂ井连线与Ｋｘ方向的夹角（θ）和渗透率方向系数（ｍ），根据Ｋｘ、θ和ｍ即可求出Ａ井与Ｂ井、Ｃ井、Ｄ井连线方向的渗透率，得出该基本单元的主渗透率初值；求解出若干基本单元后，进行某种加权平均，可得到模拟层及整个模型的网格方向。以长庆气田北区陕８１井组４口生产井２个主力生产层为计算实例。该方法在今后的数值模拟研究中有一定借鉴意义。图２表１参２（陈志宏摘）