岩石宏观流变性质的数值模拟-岩石学报.PDFVIP

１００００５６９／２００８／０２４（０６）１４１７２４ ＡｃｔａＰｅｔｒｏｌｏｇｉｃａＳｉｎｉｃａ　岩石学报 岩石宏观流变性质的数值模拟 崔晓佳　张怀　石耀霖 ＣＵＩＸｉａｏＪｉａ，ＺＨＡＮＧＨｕａｉａｎｄＳＨＩＹａｏＬｉｎ 中国科学院研究生院计算地球动力学实验室，北京　１０００４９ ＬａｂｏｒａｔｏｒｙｏｆＣｏｍｐｕｔａｔｉｏｎａｌＧｅｏｄｙｎａｍｉｃｓ，ＧｒａｄｕａｔｅＵｎｉｖｅｒｓｉｔｙｏｆＣｈｉｎｅｓｅＡｃａｄｅｍｙｏｆＳｃｉｅｎｃｅｓ，Ｂｅｉｊｉｎｇ１０００４９，Ｃｈｉｎａ ２００７０６０４收稿，２００７０９２５改回． ＣｕｉＸＪ，ＺｈａｎｇＨｕａｉａｎｄＳｈｉＹＬ．２００８．Ｎｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎｏｎｒｈｅｏｌｏｇｙｏｆｍｕｌｔｉｍｉｎｅｒａｌｒｏｃｋｓ．ＡｃｔａＰｅｔｒｏｌｏｇｉｃａＳｉｎｉｃａ， ２４（６）：１４１７－１４２４ Ａｂｓｔｒａｃｔ　　Ｔｈｅｆｉｎｉｔｅｅｌｅｍｅｎｔｍｅｔｈｏｄｗａｓｕｓｅｄｔｏｓｔｕｄｙｔｈｅｒｈｅｏｌｏｇｉｃａｌｐｒｏｐｅｒｔｙｏｆｒｏｃｋｓｃｏｍｐｏｓｅｄｏｆｍｕｌｔｉｐｌｅｍｉｎｅｒａｌ，ｇｉｖｅｎｔｈｅ ｒｈｅｏｌｏｇｉｃａｌｐｒｏｐｅｒｔｉｅｓｏｆｍｏｎｏｍｉｎｅｒａｌｓｕｎｄｅｒｃｅｒｔａｉｎｔｅｍｐｅｒａｔｕｒｅａｎｄｐｒｅｓｓｕｒｅ．Ｗｅｒａｎｄｏｍｌｙｐａｒｔｉｔｉｏｎｅｄｔｈｅｒｏｃｋｖｏｌｕｍｅ，ｄｅｓｉｇｎａｔｅｄ ｄｉｆｆｅｒｅｎｔｍｉｎｅｒａｌｓｔｏｓｕｂｐａｒｔｉｔｉｏｎｓ，ａｎｄｅｓｔａｂｌｉｓｈｅｄｔｈｅ３Ｄｆｉｎｉｔｅｅｌｅｍｅｎｔｍｏｄｅｌ．Ｔｈｅｎｗｅｃｏｍｐｕｔｅｔｈｅｓｔｒｅｓｓａｎｄｓｔｒａｉｎｒａｔｅｕｎｄｅｒ ｐｒｅｓｓｕｒｅｔｏｇｅｔｔｈｅｒｏｃｋ＇ｓｒｈｅｏｌｏｇｉｃａｌｐｒｏｐｅｒｔｙ，ｗｉｔｈｎｏｎｅｅｄｔｏｄｏｒｅａｌｒｏｃｋｅｘｐｅｒｉｍｅｎｔｓｕｎｄｅｒｈｉｇｈｐｒｅｓｓｕｒｅａｎｄｈｉｇｈｔｅｍｐｅｒａｔｕｒｅ．Ｗｉｔｈ ｔｈｅａｓｓｕｍｐｔｉｏｎｔｈａｔｉｎａｍｕｌｔｉｍｉｎｅｒａｌｒｏｃｋｓｉｎｇｌｅｍｉｎｅｒａｌｓｄｏｎ＇ｔｒｅａｃｔｗｉｔｈｅａｃｈｏｔｈｅｒａｎｄｂｅｈａｖｅｔｈｅｓａｍｅｗａｙａｓｔｈｅｙａｒｅｓｅｐｅｒａｔｅｄ， ｏｕｒｍｅｔｈｏｄｃａｎｂｅｅｆｆｉｃｉｅｎｔａｎｄｅｆｆｅｃｔｉｖｅ．ＩｎｏｕｒｒｅｓｅａｒｃｈｗｅｄｉｄｔｈｅｒｈｅｏｌｏｇｙｅｘｐｅｒｉｍｅｎｔｏｎｒｏｃｋｓｃｏｍｐｏｓｅｄｏｆＧｅｏｌｉｖｉｎｅａｎｄＧｅ ｓｐｉｎｅｌ．Ｆｏｒａｃｅｒｔａｉｎｍｉｎｅｒａｌｃｏｍｐｏｎｅｎｔｒａｔｉｏ，ｗｅｇｅｎｅｒａｔｅｒａｎｄｏｍｌｙｍａｎｙｍｏｄｅｌｓ，ａｎｄａｎａｌｙｓｉｓｓｔａｔｉｓｔｉｃａｌｌｙ．Ｗｅｆｉｎｄｔｈａｔｔｈｅ ｒｈｅｏｌｏｇｉｃａｌｐａｒａｍｅｔｅｒｓｏｆｔｈｅｒｏｃｋａｒｅｒｅｇｕｌａｒｌｙｄｉｓｔｒｉｂｕｔｅｄｂｅｔｗｅｅｎｔｈｅｃｏｒｒｅｓｐｏｎｄｉｎｇｐａｒａｍｅｔｅｒｓｏｆｍｏｎｏｍｉｎｅｒａｌｓ，ｂｕｔｔｈｅｒｅａｒｅ ｄｉｆｆｅｒｅｎｃｅｓｄｕｅｔｏｔｈｅｄｉｓｔｒｉｂｕｔｉｏｎｏｆｍｉｎｅｒａｌｓ．Ｓｏｍｅｔｉｍｅｓｔｈｅｄｉｆｆｅｒｅｎｃｅｉｓｌａｒｇｅ，ａｎｄｔｈａｔｉｓｒｅｌａｔｅｄｗｉｔｈｔｈｅｒａｔｉｏｏｆｍｉｎｅｒａｌｓ． Ｋｅｙｗｏｒｄｓ Ｆｉｎｉｔｅｅｌｅｍｅｎｔ，Ｍｕｌｔｉｍｉｎｅｒａｌｒｏｃｋ，Ｒｈｅｏｌｏｇｙ，Ｓｉｍｕｌａｔｉｏｎ 摘　要　　本文用有限元方法研究多矿物岩石的流变性质，即在一定温度压力范围内单矿物流变性质已知的情况下，采用随 机网格划分和材料性质指定的方法，建立三维多矿物岩石的有限元模型，在计算机上计算该岩石模型在不同温度压力下的变 形，根据这种数值模拟实验的结果，确定岩石整体的流变性质，而不必再进行真实的多矿物岩石的高温高压实验。用有限元 方法进行多矿物岩石流变性质的研究可以更为高效和节省，但其前提是多矿物间没有发生化学反应，变形机制仍与单矿物实 验时的主导机制一样。本研究中进行了由锗橄榄石和锗尖晶石组成的双矿物岩石的流变实验，鉴于单个的随机生成的模型 不一定居有代表性，因此对每一种矿物比例均随机生成了多个岩石模型，进行大量ＭｏｎｔｅＣａｒｌｏ实验和统计分析。计算结果表 明，由两种矿物组成的双矿物岩石整体流变性质介于两者之间且按一定规律分布，与矿物比例有关；但即使同等矿物比例混 合后的岩石会由于内部矿物分布结构的不同仍会出现一定差异，特别在试件尺度小时这种效应更为明显，因此多矿物岩石实 验试件应该有足够大的尺度。 关键词　　有限单元；多矿物岩石；流变性；模拟 中图法分类号　　Ｐ５８９．１；Ｐ５