低渗透多孔介质中非线性渗流规律LB模拟-%27--.pdfVIP

塾主塑墨量些煎盆盟 383 低渗透多孔介质中的非线性的渗流规律的LB模拟4 刘日武 周富信 (中国科学院力学研究基础室100060) 摘要t根据低渗透多孔介质孔隙结构特征，选择LatdceBoltzmann方法作为研 究低渗透多孔介质中的非线性渗流规律的模拟手段，对于中高渗多孔介 质，验证了LatticeBoltzmann方法用于多孔介质中渗流问题研究的正确 性，从模拟结果得到了流量与压力梯度呈线性关系、流量与粘度呈反比 关系、流量与渗透率呈线性关系。在此基础上研究了液体以气体在低渗 多孔介质中非线性渗流规律，与实验结果的表征是一致的，说明了孔隙 结构是引起低渗多孔介质中的非线性渗流规律的主要因素。、 关键词多孔介质非线性渗流力学数值模拟La岖ceBoltzmann方法 低渗透多孔介质是指孔隙度较低渗透性较差的孔隙介质材料【1l，如混凝土结构、陶瓷 材料、地下致密的岩石等物质。液体在低渗透多孔介质中的流动称为低渗透多孔介质渗 流，这种渗流有别于中高渗透性多孔介质中的渗流，最本质也是最明显的一点就是其中的 流动规律不再符合经典的渗流规律——D懈y定律，大量的室内模拟实验和矿场试验已经 对此做出了极好的证明，该问题的解决对多孔介质中的低速流动的理论和工程应用都有着 积极的意义。但是如何从理论模拟上对低渗透多孔介质渗流机理作出论证目前国内外尚未 见报导，本文以低渗透多孔介质的孔隙特征为突破点进行了尝试性的研究，并获得了较满 意结果。 1基本理论方法 要从理论上研究低渗透多孔介质的渗流机理，必须构造合理的低渗透多孔介质孔隙结 构模型并选择合适的计算方法。首先根据低渗透多孔介质孔隙大小符合Maxwell分布的特 点，建立了低渗透多孔介质孔隙结构的数值化模型。然后选择了Lam∞Boltzmann方法“l Boltamnn方 作为首选的计算方法。选择LattioBBoltanann方法的理论依据是因为bt曲e 法具有区别于其他数值方法的以下特征：第一，在相空间(或速度空间)的LB方法的对流因 子(或过程)是线性的。这一特征是借自能量理论并且和宏观描述的其他方法在线性对流项 相反。简单的对流项与松弛方法相结合，允许我们通过多重尺度展开处理非线性对流项。 第二，LBM方法在近不可压缩极限条件下可以获得不可压缩的Marie-Stokes方程。用状态 方程可以计算LB方法的压力。相反，在不可压缩的Marie-Stokes方程的直接数值模拟 中，压力满足一个具有类似源作用的速度级数项的Possion方程。为获得压力求解此方程 本研究得到中自博士后基金支持和中饵科学院涪流力学专项基金支恃 384 1堂生士垦谴±星芏查盘金迨塞叁 经常产生需要用叠代或松弛一类的特殊处理才能解决的数值困难。第三：在相空间中 LBM方法使用最小的速度集。在具有Maxwe-BoltTmaann平衡分布的传统的能量理论中， 相空间是完全的函数空间。平均的方法包含来自整个速度相空问的所有信息。因为在 LBM方法中只有使用了一两个速度和很少的运动方向，因此从微观分布函数到宏观量的 变换关系是十分简单的，并且只是由简单的代数运算组成。以下是LatticeBoltzmann方法 的基本方程： Q“x，O 动理演化方程：L(x+e。Ax，t+at)=A(x，tH 碰撞项：n。：一三(毛一，?) 一 f 质量守恒方程：∑la=P 动量守恒方程：∑毛e。=pu 2 Lattice Boltzmann方法用于多孔介质中流体流动的验证 多jL低渗流中最基本的规律是达西定律，其微分形式：q=一二gradp。为了验证达西 “ 定律，主要做了以下三方面的工作：①流量与压力梯度关系的模拟，得到的流量与压力梯 度呈线