摘要：以岩石实验中矿物的几何形态及空间分布为建模依据,以实验条件及单矿物电导率的测量结果为约束条件,用有限元数值方法模拟了不同微观结构的斜长石、辉石混合物在施加电压后电势及电流的分布情况,并计算了混合模型在不同温度条件下的电导率.研究结果显示,数值模型网格数及矿物颗粒数的选取对电导率计算结果的精度有较大影响,在体导电情况下,模型电导率因矿物比例含量和排列结构而异.当斜长石及辉石随机分布时,随着辉石含量的增加,混合模型电导率在不同温度下均有所增加,且温度越高,增加幅度越大,电导率的有限元模拟计算结果接近于有效介质渗透理论模型,且位于并、串联模型之间以及HS模型的上、下边界范围内;在斜长石及辉石含量一定的情况下,各矿物的排列分布对电导率计算结果也有一定的影响,当矿物颗粒大小接近且分布均匀时,模型中电势沿电流传导方向变化较为均匀,模拟计算得出的电导率相对较高,当矿物颗粒大小差别较大及分布不均匀时,电势分布受到一定的扰动,电导率计算结果也较低.将混合模型电导率有限元计算结果与辉长岩、辉绿岩及玄武岩实验测量结果进行比较,显示这3种岩石样品电导率与温度变化关系的斜率均与混合模型计算结果的斜率相接近,表明这些岩石在所选温度段导电机制与斜长石、辉石混合模型相似,用斜长石、辉石混合模型的电导率研究玄武岩、辉长岩及辉绿岩的导电性具有适用性.将混合模型有限元计算结果与玄武岩、辉长岩、辉绿岩覆盖区地壳大地电磁实测结果对比,发现大地电磁电导率结果位于混合模型计算结果范围内,用斜长石、辉石混合模型模拟玄武岩、辉长岩等岩石地壳具有一定的可行性.

关键词：有限元数值模拟 矿物混合模型 电导率 

单位：中国地震局地球物理研究所; 北京100081; 中国地震局地震观测与地球物理成像重点实验室; 北京100081; 中国地质大学(北京)地球物理与信息技术学院; 北京100083; 天津华北地质勘查局; 天津300170; 北京市劳动保护科学研究所; 北京100054