可变形多孔介质渗透系数的测定方法α.PDF

第 13 卷　第 3 期 实　验　力　学 V o l. 13　N o. 3 1998 年 9 月 JOU RNAL O F EXPER IM EN TAL M ECHAN ICS Sep. 1998 　 可变形多孔介质渗透系数的测定方法 徐曾和　徐小荷　许继军 ( 东北大学岩石破裂与失稳研究中心, 沈阳, 110006) 摘要　在B io t 理论基础上给出可变形多孔介质耦合渗流基本方程; 求出小试件 一维定常耦合渗流问题的解答; 表明在一维流固耦合情况下试件内部压力梯度有明显 的非均匀性。因此通过实验确定可变形多孔介质渗透系数在数学上可归结为微分方程 的反问题, 传统的测试渗透系数的方法需要改进。介绍了可变形多孔介质渗透系数的 测试原理和测试方法。对粒状多孔材料实验的结果表明, 传统实验方法得到的渗透系 数误差较大。 关键词　可变形多孔介质　渗透系数　流固耦合 1　前言 　　流体通过多孔介质的渗透是许多工程学科的基础。无论是地下渗透, 还是工程渗流, 所涉 及的多孔介质大多是可变形体。此时多孔介质变形与其中流体流动存在比较强的相互作用, 分 别研究流体渗流与介质变形已无法反映这类问题的特点。可变形多孔介质渗流的最主要特点 ( ) ( ) 是介质渗透率 或渗透系数 是介质变形 或介质应力和孔隙压力 的函数, 在工程上一般采用 实验的方法确定具体多孔介质的渗透系数。 在经典渗流力学中, 由于忽略介质变形与流体渗流的相互作用, 在一维稳定渗流的情况 下, 介质内部的压力梯度是常量。这样可通过小试件一维渗流实验直接确定介质的渗透系数。 对于可变形多孔介质, 通过一维稳定渗流试验确定介质的渗透系数, 仍是最方便的方法。但此 时许多研究者仍假设介质内部压力梯度是均匀的常量, 通过曲线拟合的方法, 由测试值得到介 质渗透系数的经验函数。事实上, 在流固耦合的情况下, 即使是一维稳定渗流, 试件内部的压力 梯度也是非均匀的, 因此由实测值直接拟合试验曲线的方法失效。 由于试件内压力梯度非均匀, 首先需要建立适当的数学模型, 该模型能反映介质渗透系数 与实测值的关系。这样由试验确定渗透系数的问题在数学上归结为微分方程的反问题, 一般情 况下是函数待定的反问题, 在特殊情况下可简化为参数待定的反问题。 本文首先按B io t 建立三维固结理论的基本思想, 对其理论加以扩展, 给出可变形多孔介 本文于 1997 年 7 月 27 日收到第1 稿, 1998 年 4 月 21 日收到修改稿 © 1995-2004 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved. 第 3 期　　　　　　　　　徐曾和等: 可变形多孔介质渗透系数的测定方法　　　　　　　　　　　 315 质耦合渗流的基本方程, 求出一维定常耦合渗流的解答, 然后介绍了可变形多孔介质渗透系数 的测试原理与测试方法, 对粒状多孔材料试验的结果表明, 传统方法测得的渗透系数误差较 大。 2　可变形多孔介质耦合渗流的基本方程 　　 提出的三维固结理论, 最早考察了三向压缩情况下多孔介质变形与孔隙流体流动的 B io t 相互作用。B io t 理论是建立在线性弹性和D arcy 定律的基础上的。即假设多孔介质的变形线 性地依赖于介质应力和孔隙流体压力, 孔隙率的改变也线性地依赖于介质应力与孔隙压力, 且 忽略了孔隙压力对介质剪切变形的影响, 这样按照岩石力学中以压为正以拉为负的惯用规定, 可得本构方程为[ 1 ] = 2G + ( + p ) ( , j = 1, 2, 3) ( ) ij ij ij v