周期结构混合介质电性谱的三维有限差分模拟

论文摘要

混合介质电性谱的研究对于矿产勘探、材料设计、医学检查和环境监测都有很重要的意义。介电常数和电导率的频散现象广泛地存在于各种天然介质和人工介质中。近来,越来越多的实验发现一些普通混合介质的介电常数在低频域会出现异常大的增强现象,对此人们争议还比较大,还没有统一的解释。本文采用三维有限差分方法求解混合介质内部的电场分布,进而提取混合介质的等效介电常数和电导率,在差分网格的划分中采用了一种近似方法来解决含薄膜层的介质的数值模拟问题,通过经典解析公式的验证,表明本文对膜结构的近似处理是合理的,没有引起大的误差。本文主要从以下几个方面对混合介质的等效电性参数及其频散特性进行了深入研究:（1）对无损耗介质的等效介电常数进行了简单的分析,阐述了其等效介电常数的主要控制因素和变化规律。对无损耗介质的模拟计算表明:在包裹体体积较小时,而且包裹体和宿主的介电常数比值不大的情况下,Maxwell-Garnett公式和FDM的计算结果非常接近。当包裹体体积较大时,Maxwell-Garnett公式就不太适用了,所以其与FDM的计算结果差异也较大。在包裹体介电常数大于宿主介电常数的情况下,嵌入比越大,等效介电常数随包裹体介电常数增长而增长的速度越快;球形包裹体嵌入比在0.4到0.7之间时,等效介电常数随嵌入比增长速度最快。（2）对一般混合介质的电性谱进行了分析,研究了包裹体的几何参数和电性参数对混合介质电性谱的影响规律,并分析了弱介电增强和弱电导率增强现象及其产生机理。对一般混合介质的模拟表明:在某些情况下,等效介电常数会在低频端出现微弱的增强;某些情况下,等效电导率会在高频端出现微弱的增强。等效介电常数的增强和介质各组分电导率的差异有关;等效电导率的增强和介质各组分介电常数的差异有关。等效介电常数的增强是由于等势线在传导电流为主导时（低频端）在局部汇集引起的;等效电导率的增强是由于等势线在位移电流为主导时（高频端）在局部汇集引起的。其中,等效电导率的增强现象尚未见有报道。对于各组分介电常数相同的混合介质,包裹体和宿主的电导率比值控制着等效介电常数在低频端的增强,等效电导率介于包裹体和宿主电导率之间。对于各组分电导率相同的介质,包裹体和宿主的介电常数比值控制着等效电导率在高频端的增强,等效介电常数介于包裹体和宿主介电常数之间。当各组分介电常数同时增大但比值不变时,等效电导率在高频端的增强幅度不变,但是电性谱的驰豫频率会发生平移。当各组分电导率同时增大但比值不变时,等效介电常数在低频端的增强幅度不变,但是电性谱的驰豫频率会发生平移。等厚的两层介质模型的计算表明:对于介电常数相同,电导率不同的两层介质,随着频率的降低,等势线会逐步集中到电导率较小的介质中,同时,在低频端等效介电常数会增大一倍;对于电导率相同,介电常数不同的两层介质,当频率升高时,等势线会逐步集中到介电常数较小的介质中,同时,在高频端等效电率会增大一倍。对于介电常数相等,中间一层电导率较低的三层介质模型,发现中间一层厚度越小,低频端的介电增强作用越强,同时电导率曲线整体向左上方平移。对于含球形包裹体的模型,当包裹体嵌入比R小于0.5时,低频端的等效介电常数随着R的增加而增加;当包裹体嵌入比R大于0.5时,低频端的等效介电常数随着R的增加而降低。所以说R=0.5是一个临界状态,因为当包裹体的嵌入比R大于0.5时,这些包裹体就相互连接起来,如果包裹体是导电的,那么混合体系的导电性就会极大地增强。所以当R大于0.5后,混合介质的等效电参数会有很大的变化。椭球包裹体的倾角、椭球包裹体的纵横比和圆柱包裹体的纵横比都对电性谱有一定影响,其中相似的一个特点是:包裹体在所施加电场方向的延伸越长,低频端的介电增强作用越大。混合介质中的组分越多,其电性谱越复杂,对于含三种不同组分的混合介质,其电性谱会出现两次驰豫。在准静态电场的假设下,等比例缩放模型时,电性谱不发生改变。含八个单元和一个单元的模型的电性谱完全一致,说明基于混合定律对周期结构混合介质的简化处理是正确的。（3）对含薄膜层的介质进行了分析,发现低电导率薄膜会引起等效介电常数在低频端出现异常大的增强。薄膜电导率、薄膜厚度、薄膜面积对低频端的介电增强作用有很大影响。包裹体的几何形状、倾角等因素对介质的电性谱也有一定影响。类似地,低介电常数薄膜会引起等效电导率在高频端出现很大的增强。当膜的厚度大于10-8m时,介电增强作用随着膜厚度的减小而增强;当膜的厚度小于10-8m时,介电增强作用随着膜厚度的减小而减弱。因为当薄膜厚度太薄时,薄膜就起不到阻断传导电流的作用了。随着膜电导率的降低,低频端的介电增强作用不断增强,但是这种增强作用随着膜电导率的降低逐步变弱,最后当膜的电导率小到一定程度后,低频端的介电常数趋向一个饱和值(如:膜电导率σm=10-8S/m和σm=0对应的曲线基本重合)。因为膜电导率越低,等势线越向膜层中汇集,当膜电导率低到一定程度后,等势线就全部汇集到膜层中了,这样就到达一种饱和状态了。随着膜面积（水平层状薄膜）的增大,低频端的介电常数越来越大（介电增强作用越来越强）;当膜层完全将模型上部和下部隔开时,（即:将导电通路完全切断）,这时传导电流很难穿过整个模型,而低频时,以传导电流为主,所以低频端的等效电导率就特别低;高频时以位移电流为主,所以高频端的等效电导率未受影响。含膜的两层介质的计算表明:低频端的介电增强幅度和薄膜参数之间存在着定量关系;同样地,高频端的电导率增强幅度和薄膜参数之间也存在着定量关系。等效介电常数在低频端的增强是电容模型的有效厚度在低频端减小引起的;等效电导率在高频端的增强是电阻模型的有效长度在高频端减小引起的。所以可以认为混合介质等效电参数的频散现象是因为样本的有效尺寸随频率发生改变而引起的。（4）对比了经典频散模型和本文数值计算的电性谱的变化规律,分析了经典公式和数值方法中模型参数对电性谱的控制作用。介电常数和电导率的经典频散模型（Cole-Cole模型）分别用了4个参数来控制介电常数和电导率频散曲线的形态和位置,但是其中3个参数的物理意义不明确。本文的数值模方法通过改变混合介质中各个组分的电性参数和几何参数也可以实现对散曲线形态和位置的控制。数值方法计算频散曲线的速度远远慢于经典解析模型,而且引入的模型参数更多,但是数值方法中采用的参数都有明确的物理意义,更方便分析介质频散特性和介质本身的物理参数之间的关系。（5）对几种岩石模型的电性谱进行了分析,证实在含有机质页岩、含金属颗粒的岩石、亲油岩石等岩石中会出现巨大的介电增强现象。对于饱水纯砂岩模型,在砂岩孔隙相互连通的情况下,等效介电常数和电导率都不会随频率发生变化。对于完全饱油的纯砂岩模型,电性谱会出现频散现象,而且低频端会出现较弱的介电增强现象。饱含低电导率流体的裂缝模型的计算表明:这种介质会出现很大的介电增强现象。含有机质页岩的模型的计算表明:低电导率的有机质层片会引起很大的介电增强现象。对于含金属颗粒的岩石模型,如果金属颗粒表面存在低电导率膜层（氧化壳）,该介质就会在低频端出现很大的介电增强。亲油岩石由于岩石颗粒表面覆盖着油性膜层,所以在低频端也会出现很大的介电增强,但其等效电导率基本不随频率发生变化。

论文目录

• 摘要
• ABSTRACT
• 第一章 绪论
• 1.1 本文的研究目的与意义
• 1.1.1 研究目的
• 1.1.2 研究意义
• 1.2 研究现状综述
• 1.2.1 基于有效介质理论的等效电性参数研究
• 1.2.2 岩石电性参数频散现象的研究
• 1.2.3 介电增强现象
• 1.3 主要研究内容
• 1.4 研究方法与研究成果
• 1.4.1 研究方法
• 1.4.2 拟解决关键问题
• 1.4.3 主要创新点
• 第二章 基于有限差分方法的建模和计算
• 2.1 混合定律与计算模型
• 2.2 基于3D-FDM 的等效电性参数计算
• 2.3 模型参数的定义
• 2.4 3D-FDM方法的验证
• 2.5 FDM 计算精度和效率的评价
• 2.5.1 网格密度的影响
• 2.5.2 终止迭代的误差容限的影响
• 2.5.3 计算耗时的估计
• 2.6 本章小结
• 第三章 一般混合介质的电性谱
• 3.1 无损耗混合介质的等效介电常数
• 3.2 有损耗混合介质的等效介电常数和电导率
• 3.2.1 各组分电导率比值的影响
• 3.2.2 各组分介电常数比值的影响
• 3.2.3 各组分介电常数同时增大但比值不变
• 3.2.4 各组分电导率同时增大但比值不变
• 3.3 层状模型的分析
• 3.3.1 两层模型
• 3.3.2 三层模型
• 3.4 包裹体几何参数的影响
• 3.4.1 球形包裹体半径的影响
• 3.4.2 椭球倾角的影响
• 3.4.3 椭球包裹体纵横比的影响
• 3.4.4 圆柱体纵横比的影响
• 3.5 三种组分的混合介质
• 3.6 模型尺度的影响
• 3.7 混合定律的检验
• 3.8 本章小结
• 第四章 含薄膜结构的混合介质的电性谱
• 4.1 球形带膜包裹体
• 4.1.1 薄膜层的影响
• 4.1.2 薄膜厚度的影响
• 4.1.3 薄膜电导率的影响
• 4.1.4 球形包裹体半径的影响
• 4.1.5 模型尺度的影响
• 4.2 包裹体几何参数的影响
• 4.2.1 包裹体形状的影响
• 4.2.2 包裹体倾角的影响
• 4.2.3 椭球包裹体纵横比的影响
• 4.2.4 圆柱纵横比的影响
• 4.2.5 膜的面积的影响
• 4.3 含膜的两层模型
• 4.3.1 介电常数剧增
• 4.3.2 电导率剧增
• 4.4 含多种包裹体的混合介质
• 4.5 与经典频散模型的对比
• 4.5.1 介电常数的频散模型
• 4.5.2 电导率的频散模型
• 4.5.3 讨论
• 4.6 本章小结
• 第五章 几种岩石模型的计算
• 5.1 一般岩石的简单模型
• 5.1.1 饱水纯砂岩
• 5.1.2 饱油纯砂岩
• 5.2 岩石的介电巨增
• 5.2.1 含低电导率流体的微裂缝模型
• 5.2.2 含有机质的页岩
• 5.2.3 含金属颗粒的砂岩
• 5.2.4 亲油岩石
• 5.3 本章小结
• 第六章 结论与问题
• 6.1 主要结论
• 6.2 存在问题与对后续研究工作的建议
• 致谢
• 参考文献
• 附录 常用包裹体的几何参数设置