多孔介质相对渗透率的耦合模型-凝聚态物理专业论文.docxVIP

华 华 中 科 技 大 学 硕 士 学 位 论 文 PAGE 10 PAGE 10 1 绪 论 1.1 多孔介质 多孔介质——或者叫多孔材料——的定义是内部含有孔隙的固体。在微观层面 上，所有的固体分子之间都极细微的空隙，叫做“分子间隙”；而在宏观层面，对很 多固体材料而言，一些空隙相对其自身尺度不可忽略，这些空隙我们可以称其为“洞 穴”。在本文中，以及大多数讨论多孔介质性质的文章中，孔隙的尺度都是介于“分 子间隙”和“洞穴”之间的。这些孔隙可以是相互连通的也可以是不连通的，流体 可以在相互连通的孔隙中流动，这部分空间和整个多孔介质体积的比叫做多孔介质 的孔隙率。 按照这一定义，自然界中的很多事物都可以归纳为多孔介质来研究。大多数岩 石是具有孔隙的，例如泡沫岩、石灰岩、白云岩等[1]，地下油藏的开采、岩土力学的 研究都离不开多孔介质理论；土壤是很多细小颗粒和小石块组成的，对于研究土壤 孔隙率、渗透特性、透气率等和农业生产密切相关的方面，多孔介质理论显得尤为 重要；生物体内的组织器官等，由于具有盘根错杂的导管系统，和细胞组织间存在 的孔隙，也属于多孔介质。当然，除了天然形成的多孔介质以外，人造材料也有很 多属于多孔介质。例如航天复合材料、建筑保温材料、工程陶瓷等，渗透率、孔隙 率、导热率等，是这些材料在设计和提高性能是的重要参数；另外，广泛用于芯片 设计的半导体材料也属于多孔介质，孔隙率对于半导体材料的导电性能具有很大的 影响。 多孔介质虽然广泛的存在于自然界和生产实际中，但是由于其不规则的几何形 状和复杂的结构，给其研究带来了很大的困难。人们只能通过实验[2][3][4]和数值模拟 [5][6][7]来确定多孔介质的特性和参数。不过随着分形理论的发展，大量报道揭示了多 孔介质的微结构、孔隙分布、颗粒大小满足分形分布规律。从而为我们通过理论分 析和预测多孔介质特性和参数提供了数学工具。已有研究表明，天然形成的多孔介 质（如土壤、油气藏、煤层气结构、地下水资源、沙石等）和部分人工制造的多孔 介质（如硅晶片微结构）、哺乳动物的脏器与循环系统、植物根系、枝叶及其微循环 系统、颗粒型多孔介质、多孔纤维织物等都是分形多孔介质[8][9][10][11][12]。 1.2 渗透率概述 对于多孔介质，我们通常关心气体或流体在多孔介质的固体骨架中的传输特性。 描述多孔介质流动特性的一个重要物理量叫做渗透率，其定义基于最初研究多 孔介质中流动特性的达西定律。1856 年，法国物理学家 H. Darcy 在研究第荣城（Dijon） 的给水问题时，用直立砂柱进行了渗流实验，并得出了他的结论，也就是著名的达 西定律： Q = K ? A ? ΔP  （1-1） μ L 流体通过砂柱横截面的体积流量 Q 与横截面积 A 和水头压力差 ΔP 成正比，而与 砂柱长度 L 和流体粘度 μ 成反比。比值 K 称为介质的渗透率。它与流体性质无关， 而只是多孔介质本身结构的函数，所以我们也称它为多孔介质的绝对渗透率。 多孔介质渗透率的研究方法也是多种多样的，下面我们逐一介绍： 1.2.1 实验直接测量 有些人用实验方法直接测量多孔介质的渗透率，并将实验结果综合成各种经验 关系式，其缺点是每个经验关系式中通常有几个经验常数。如Adler and Thovert[7]报 告了对于玻璃粒子多孔介质，渗透率经验公式为： eK = 0.117r 2φ 4.57 e  （1-2） 上式中 re 为粒子等效半径,φ 为孔隙度。 实验直接测量方法虽然最为直接可靠，但不可避免的存在拟合常数，而且不能 很好地从物理原理上解释所测量到的数据。 1.2.2 数值模拟方法 通过数值模拟方法测量渗透率[7]错误！未找到引用源。的一般采用以下方法来模拟： A. 构造一个周期性的几何结构来模拟该复杂介质； B. 有限元法、有限差分法、自动格子机方法、Monte Carlo、Lattice Boltzmann Method； C. 所得结果通常以图表或经验关系式给出。 数值模拟方法有其固有的缺点，模拟的结果无法与实际多孔介质微结构参数相 匹配。其结果的物理意义往往最不明确。 1.2.3 理论解析方法 对于一些理想结构（如圆形毛细管、矩形孔道）的多孔介质，前人给出了一些 渗透率的分析解[1][14]。对于一类多孔纤维织物，文献[14]解析地推导并获得了渗透率 的通用分析解模型： K = ag H g (H g /12 +  k⊥ / 2) / A  （1-3） 式中 ag 为最大孔隙面积， H g 为孔隙直径， A 为周期单位面积， k⊥ 为纤维束的 横向渗透率。其中没有经验常数，其预测值与渗透率实验结果吻合很好。可此模型 只适合结构规则的多孔介质，对大多数无序多孔介质，此