部分熔融岩石波速的有限元数值模拟.pdfVIP

·296· 中国地球物理2011 (8)地球内部结构及其动力学 部分熔融岩石波速的有限元数值模拟 NumericalSimulationforthe of MoltenRocks FEM VelocityPartially Using 孟祥翠’陈祖安 ChenZuan MengXiangcui 中国科学院地质与地球物理研究所北京100029 地震波探测发现在大陆地壳内部．特别是在当今构造活动带的不同深度上断续地存在壳内低速层． 例如，在青藏高原的15～35 kin／ S和3．2～3．5 km的区 km／s：在南纬2l。～23。的安第斯山脉的中段，自西科迪勒拉到东科迪勒拉长约350 值为特征。很多学者认为部分熔融和脱水作用是引起岩石地震波速度降低的主要原因(Ito，1990：宋茂双 等，1996；高平等，1994；林传勇等，2001)。 新近的地震探测揭示出青藏高原地壳深部(20～50km)具有广泛和显著的各向异性，其各向异性的平均 et et et 强度达到10％，甚至更高(Ozacaral，2004；Sherringtona1．，2004；Shapiroal，2004)。理论分析和实验 结果均表明，地壳深部岩石的各向异性可以由矿物定向分布所引起，但是由LPO引起的地壳深处各类岩石 的平均各向异性强度通常不超过5％，远不足以解释在青藏高原地壳中所观测到的弹性波各向异性的强度。 地质和地球物理证据均表明部分熔融在青藏高原地壳深部可能广泛存在。熔体对构造变形响应非常明显。 在构造应力作用下，熔体以定向分布的扁平状熔体囊或熔体条带形式存在，熔体囊长轴或熔体条带的延伸方 向与剪切面的夹角为15。～25。【l~21。熔体的定向分布有可能显著地改变岩石原来的对称性，从而改变岩石 的弹性性质(弹性波速和各向异性)，影响地震测深的结果。 中外学者【3-41都有通过高温高压实对部分熔融岩石波速进行过研究，结果表明岩石中出现熔体时会导 致波速出现突然变，而且波速随着部分熔融岩石中熔体含量的增加而持续降低。但是定量评价部分熔融程 度对于弹性波波速的影响相当困难，因为波速同时受部分熔融岩石中熔体的形态、分布的影响，而在实验 过程中对于熔体的分数和熔体分布形态很难人为控制。 作者希望采用数值模拟的方法对部分熔融岩石中熔体对于波速的影响进行定量的研究。思路是通过热 弹性有限元方法来计算部分熔融岩石在受力作用下的应力应变情况，然后根据虚功原理计算出平均应力和 平均应变，继而计算出杨氏模量和泊松比，而波速与杨氏模量和泊松比存在一定的关系，所以可以计算出 部分熔融岩石的波速。我们可以通过改变部分熔融岩石中熔体的比例，研究波速随熔体含量的变化情况： 改变部分熔融岩石中熔体的分布形态，研究波速随熔体的分布形态的变化情况。当熔体在构造应力作用下 而成条带聚集时，研究对于波速的影响。根据计算结果得到的波速随熔体含量和温度变化的拟合曲线，可 以用来粗略的反推出地下熔体的比例。 初步计算结果表明岩石总体的杨氏模量随着熔体含量的增加而减小，而波速随着岩石中熔体含量的增 加而减小，这与高温高压实验得到的数据基本符合。熔体的形态对于波速也有比较大的影响。当熔体在岩 石中呈条带聚集分布时，岩石的波速出现大幅度的减小，而且岩石表现为各向异性的特征。 本研究由国家自然科学基助。 参考文献 D mantle RevEarthPlanet M molten Sci