石榴子石族矿物状态方程研究进展

范大伟 李博 陈伟 许金贵 匡云倩 叶之琳 周文戈 谢鸿森

范利丹, 徐峰, 余永强, 张志伟, 余雳伟, 周桂杰. 一维静载与循环冲击共同作用下砂岩动态力学特性试验研究[J]. 高压物理学报, 2022, 36(3): 034101. doi: 10.11858/gywlxb.20210879
引用本文: 范大伟, 李博, 陈伟, 许金贵, 匡云倩, 叶之琳, 周文戈, 谢鸿森. 石榴子石族矿物状态方程研究进展[J]. 高压物理学报, 2018, 32(1): 010101. doi: 10.11858/gywlxb.20170597
FAN Lidan, XU Feng, YU Yongqiang, ZHANG Zhiwei, YU Liwei, ZHOU Guijie. Experimental Study of Dynamic Mechanical Characteristics of Sandstone under One-Dimensional Coupled Static-Cyclic Impact Loads[J]. Chinese Journal of High Pressure Physics, 2022, 36(3): 034101. doi: 10.11858/gywlxb.20210879
Citation: FAN Dawei, LI Bo, CHEN Wei, XU Jingui, KUANG Yunqian, YE Zhilin, ZHOU Wenge, XIE Hongsen. Research Progress of the Equation of State for Garnet Minerals[J]. Chinese Journal of High Pressure Physics, 2018, 32(1): 010101. doi: 10.11858/gywlxb.20170597

石榴子石族矿物状态方程研究进展

doi: 10.11858/gywlxb.20170597
基金项目: 

中国科学院战略性先导科技专项(B类) XDB18010401

国家自然科学基金 41274105

国家自然科学基金 41374107

国家自然科学基金 U1632112

中国科学院“西部之光”人才培养引进计划“西部青年学者”A类项目 

中国科学院青年创新促进会专项基金 

详细信息
    作者简介:

    范大伟(1982-), 男, 副研究员, 硕士生导师, 主要从事高温高压矿物物理研究. E-mail:fandawei@vip.gyig.ac.cn

    通讯作者:

    周文戈(1967-), 男, 研究员, 博士生导师, 主要从事高温高压岩石、矿物物性研究. E-mail:zhouwenge@vip.gyig.ac.cn

  • 中图分类号: O521.2

Research Progress of the Equation of State for Garnet Minerals

  • 摘要: 石榴子石是重要的造岩矿物,是上地幔、地幔转换带以及(超)高压变质岩中最重要的造岩矿物之一,研究其状态方程对于约束地球内部物质组成和状态、正确理解大洋岩石圈俯冲板块和地幔动力学过程具有重要意义。文中综述了20世纪70年代以来石榴子石p-V(压强-晶胞体积)和p-V-T(压强-晶胞体积-温度)状态方程的研究进展,重点讨论高温高压条件下石榴子石的稳定性以及组分变化和含水对热弹性参数的影响。最后简略概括石榴子石状态方程研究存在的问题,并指出发展方向。

     

  • 浅部矿物资源的日渐枯竭使资源开采工作向地下更深处发展。矿物采掘工作往往需要多次爆破,加之大型机械设备的连续钻探施作[1-3],以及随深度增加而增大的地应力,导致巷道围岩遭受频繁的爆破冲击和机械振动,围岩体的破坏程度不断加剧,岩爆发生的概率也大大增加,严重威胁煤炭开采人员和机械设备的安全,开采进度难以保证。

    分离式霍普金森压杆(split Hopkinson pressure bar,SHPB)是目前国内外众多学者研究深部岩石动态力学性能和破坏机制的主要工具。金解放等[4-6]通过改造的SHPB动静组合加载装置,对长径比为1.0的砂岩进行了动静组合循环冲击试验,发现冲击试样呈现锥体的破坏模式,表明砂岩的单位体积吸收能随循环冲击次数的增加而增大,致使岩石内部损伤加剧、强度降低,他们还将岩石的动态应力-应变曲线划分为5个发展阶段,即压密阶段、弹性阶段、内部裂纹扩展的加载阶段、第1卸载阶段和第2卸载阶段。王志亮等[7]在岩石的单轴循环冲击试验中,利用静态应力-应变曲线,推导出静态裂纹起裂应力,进而得出动态裂纹起裂应力。田诺成等[8]对不同轴压下的花岗岩进行等幅循环冲击,发现在相同的入射波应力循环冲击下,循环冲击次数随轴压的增大先增加后减少,即轴压使循环冲击荷载作用下的岩石动力学特性呈先强化后劣化的变化趋势,结合声发射技术,发现强度强化是由于预加轴压小于岩石的起裂应力,反之则会造成岩石强度劣化。Li等[9]利用大直径(75 mm)压杆对无轴压和无围压下的花岗岩试样进行循环冲击试验,发现当动载应力不超过岩石静态强度的60%时,循环冲击基本不会对岩石造成明显损伤,该结论在文献[5]中也得到证实。唐礼忠等[10]研究岩石在预加轴压与循环冲击共同作用时发现,岩石的应力-应变曲线在不同轴压下呈现回弹与不回弹两种特性(定义为局部回弹),且岩石的破碎块度随轴压的增大而增大。

    很多理论和试验研究均表明,在静载与动载的共同作用下岩石的力学性能和破坏模式相较于静载条件下有所不同。目前,已有部分学者对砂岩开展了一维或三维动静组合加载下的循环冲击试验[11-13],但鲜有涉及一维动静组合循环加载下的岩石力学性能和能量耗散,以及冲击气压对岩石动态力学性能的影响等问题。对于深部巷道开挖时的围岩体,矿柱在承受一个方向上的静荷载(地应力)时还要承受外部较大且较频繁的动载扰动,明显的应力集中极易造成矿柱失稳。为此,本研究主要探讨砂岩在同一方向的动静组合循环加载下的力学性能、能量耗散及破坏模式,这对地下工程建设的安全评价具有重要的现实意义。

    砂岩岩样取自鹤煤三矿井下800 m巷道工作面,其完整性和均质性较好。依据国际岩石力学学会(ISRM)推荐的试样尺寸[11, 14],将砂岩加工成长径比为0.5(直径为50 mm,长度为25 mm)的圆柱试块。为了确保试验结果的精确度,精细打磨试件,使试件两个端面的不平行度小于0.02 mm,不垂直度小于0.001 rad。取样所选砂岩为天然湿度状态,测得7个试块的单轴抗压强度分别为89.31、123.75、117.47、97.78、120.36、125.94、138.36 MPa,根据离散性要求,舍去最大值和最小值,求得试样的平均单轴抗压强度(fcu)为116.94 MPa,其他静态物理力学参数见表1,其中:ft为抗拉强度,E为弹性模量,μ为泊松比,v为纵波波速。

    表  1  砂岩的静态物理力学参数
    Table  1.  Static physical and mechanical parameters of sandstone
    fcu/MPaft/MPaE/GPaμv/(m·s−1)
    116.945.3031.250.263923
    下载: 导出CSV 
    | 显示表格

    采用河南理工大学冲击实验室改进的SHPB试验装置(见图1):压杆系统由子弹、入射杆、透射杆和吸收杆等部分组成,其中子弹为直径37 mm、长度 400 mm的合金冲头,入射杆为直径50 mm、长度2400 mm的变截面合金钢,透射杆为直径50 mm、长度1200 mm的合金钢,吸收杆为直径37 mm的铝质杆,合金钢杆中的纵波波速为5190 m/s,弹性模量为210 GPa;氮气为动力源,通过调节气压控制器改变子弹冲击速度,子弹速度由红外测速仪测量;试验数据由江苏东华生产的DH8302超动态应变放大仪采集。该试验系统可进行常规的单轴、一维和三维冲击试验,亦可进行轴压和围压的手动或数控设定。

    图  1  SHPB 试验装置
    Figure  1.  SHPB test apparatus

    图2为SHPB试验装置示意图。子弹在冲击气压下以一定速度撞击入射杆,入射杆中产生入射波传递给试样,由于试块与试样的波阻抗存在较大差异,使得部分入射波经试样传给透射杆转化为透射波,一部分经试样反射回入射杆转化为反射波,另有少数被试样吸收,转化为吸收波,用于岩样裂缝的产生与发展。其中,入射波εi(t)和反射波εr(t)均由粘贴在入射杆表面的应变片A1测得,而透射波εt(t)由应变片A2测得。由于试样厚度远小于应力波波长,因此在“应力均匀性假设”的前提下,可认为试样中的应力均匀无衰减,根据一维应力波理论,由关系式εi(t)+εr(t)=εt(t),将“三波法”简化为“二波法”[15],推导出试样的动态应力σs(t)、应变εs(t)、应变率˙εs(t)分别为

    图  2  SHPB试验装置示意图
    Figure  2.  Sketch of the SHPB system
    σs(t)=E0A0Asεt(t)
    (1)
    εs(t)=2c0Lst0εr(t)dt
    (2)
    ˙εs(t)=2c0Lsεr(t)
    (3)

    式中:A0c0E0分别为弹性杆的横截面面积(mm2)、纵波波速(m/s)和弹性模量(GPa),AsLs分别为试样的横截面面积(mm2)和长度(mm)。

    选取最具代表性的砂岩试样应力测量原始数据,得到图3(其中U为电压),对入射波和反射波进行叠加,可见其波形与透射波基本重合,表明试样左右两端受力基本相等,满足应力均衡条件[16]

    图  3  砂岩的动态应力均衡检验
    Figure  3.  Calibration of dynamic stress balance of sandstone

    为探究砂岩在一维静压下的循环冲击动态力学特性,基于试验装置的加载限度,将轴压分别设置为 5、15 和 20 MPa,对应于砂岩单轴抗压强度的 4%、13%和 17%,每组试验对3个试样进行冲击。试验前,为减少压杆与试样之间的摩擦力,保证良好接触,在试样与压杆接触端面涂抹少量黄油用作耦合剂;取边长15 mm、厚0.3 mm的正方形H62黄铜片作为波形整形器贴于入射杆靠近子弹的端面;由同一细直钢杆将子弹送入发射腔内同一位置,尽可能保证相同气压下入射波的峰值相同或相近。一维动静组合岩石加载模型见图4(a),采用等截面积加载方式对试样进行循环冲击,加载的恒应变率半正弦波应力脉冲[9, 11, 17]图4(b)所示 (其中:pt为总载荷;σas为预加轴压;pc为冲击气压,即动载荷),冲击完成后对试样拍照并观察,直至试样产生宏观裂纹并完全破碎。试样的物理参数及加载方案见表2,其中:L为试样长度,D为试样直径,ρs为试样密度。

    图  4  一维动静组合加载模型示意图
    Figure  4.  Sketch of model under one-dimensional coupled static and dynamic loads
    表  2  代表性试样的基本物理参数及加载方案
    Table  2.  Basic parameters of the representative samples and the loading scheme
    Sample No.L/mmD/mmρs/(kg·m−3)σas/MPapc/MPa
    S3-1
    S3-2
    S3-3
    S3-4
    S3-5
    S3-6
    S3-7
    24.97
    25.26
    25.38
    25.26
    25.35
    25.04
    24.94
    49.56
    49.60
    49.53
    49.84
    49.90
    49.88
    49.65
    2591
    2571
    2567
    2561
    2580
    2564
    2628
    5
    15
    20
    15
    20
    15
    20
    0.8
    0.8
    0.8
    1.0
    1.0
    1.2
    1.2
    下载: 导出CSV 
    | 显示表格

    记录每块岩样的累计冲击次数,并绘制动态应力-应变曲线,相同加载及冲击类型下进行多次试验,选取最具代表性的试验数据进行分析。表3列出了试样的累计冲击次数(n)。对于试样S3-7,因试验数据采集不及时造成数据缺失,无法对其各项动态力学性能进行分析评价,但其冲击破坏形态及循环冲击次数等数据仍具有有效性。

    表  3  代表性试样的总循环冲击次数
    Table  3.  Total number of the cyclic impacts of the representative samples
    S3-1S3-2S3-3S3-4S3-5S3-6S3-7
    1510811962
    下载: 导出CSV 
    | 显示表格

    图5显示了循环冲击下试样S3-1的原始波形。可以看出,在轴压与冲击气压不变的情况下,每次冲击的入射波基本重合,说明试样每次受到的冲击荷载基本相同。首次冲击的反射波峰值最大,透射波峰值最小;随着冲击次数的增加,反射波和透射波曲线整体下移;而后,随着冲击的持续进行,反射波与透射波曲线又整体上移,直至试样发生宏观破碎。

    图  5  砂岩试样循环冲击原始波形
    Figure  5.  Original waveforms of cyclic impacts for sandstone samples

    图6为砂岩试样在循环冲击作用下的动态应力-应变曲线。可以看出:所有曲线在应力峰值前的变化趋势呈现一致性,即应力随应变呈递增关系;应力峰值前,曲线斜率逐渐减小,没有压密阶段,仅有弹性阶段和非线性变形阶段。应力峰值后,应力-应变曲线表现出两种变形特征:一是极少数的局部回弹现象,如试样S3-4的第7次冲击,可能的原因是轴压的施加使试样内部提前储存一定的弹性能,导致试样在经历多次冲击没有完全破碎的情况下内部积攒的弹性力大于外部冲击力,而在此瞬间,试样上的质点速度可能高于入射杆端面的质点速度,使得试样在以上两种因素的综合作用下于卸载阶段出现曲线局部回弹现象;二是无局部回弹,即与一般单轴压缩试验得到的应力-应变曲线变化趋势一致,该现象可能是由于试样内部的弹性能过早释放,造成卸载阶段的弹性力始终小于外部加载应力。

    图  6  砂岩试样在循环冲击荷载下的应力-应变曲线
    Figure  6.  Stress-strain curves of sandstone subjected to cyclic impact

    为进一步说明图5中冲击波形曲线的变化特征,对图6所示的砂岩试样在循环冲击作用下的动态应力-应变曲线变化特征进行分析。定义曲线峰值应力(σf)对应的应变为峰值应变(εf),应变的最大值为最大应变(εm)。在试块冲击破坏的整个过程中,动态应力-应变曲线整体出现回弹现象(定义为整体回弹),即动态峰值应力总体呈现先上升后下降的变化趋势,而最大应变则刚好相反(见图7)。在加载初期,由于试块内部具有较多孔隙及细小裂纹,大部分入射波在孔隙或裂纹表面发生反射,小部分发生透射,因此造成反射波峰值较大,透射波峰值较小,导致动态峰值应力逐渐增大,峰值应变逐渐减小,最大应变也相应减小。这说明在前几次冲击过程中,砂岩内部众多微小孔隙(裂隙)在荷载作用下不断变小甚至“愈合”,使得大部分入射波穿过试样,造成试块抵抗外界动态荷载的能力逐渐增强,即动态应力-应变曲线整体回弹。随着冲击的继续进行,砂岩的动态强度呈减小趋势,而峰值应变却增大,说明砂岩试块在持续的损伤软化作用下,内部开始产生细小裂纹,越来越多的入射波发生反射,且随着冲击次数的增加,裂纹不断发展,直至产生宏观破坏。

    图  7  峰值应力及最大应变与冲击次数的关系
    Figure  7.  Variations of peak stress and maximum strain with impact times

    图8显示了各试样每次冲击的峰值应力与循环冲击次数的二次抛物曲线拟合关系,拟合曲线的最高点为冲击回弹转折点。对比试样S3-2和S3-3可以发现,在相同的冲击气压下,峰值应力曲线随轴压的增大而整体下移,试样S3-4和S3-5也有相同的变化趋势,表明峰值应力随轴压的增大而减小,此现象说明轴压对岩石的承载能力起到削弱效果,在一定程度上验证了文献[18-20]得出的“轴压弱化岩石承载能力”这一结论。另外,回弹转折点所对应的循环冲击次数随轴压的增大而减小,即回弹转折点前移(试样S3-1和S3-3、试样S3-4和S3-5)。在相同的轴压(试样S3-2、S3-4和S3-6,以及试样S3-3和S3-5)下,随着冲击气压的增大,峰值应力曲线整体上移,表明峰值应力随冲击气压的增大而增大,而回弹转折点对应的循环冲击次数随冲击气压的增大而减小。同样地,对比图9中试样S3-2和S3-3可以发现,在相同的冲击气压下,最大应变曲线随轴压的增大而整体上移,试样S3-4和S3-5也有相同的变化趋势,表明最大应变随轴压的增大而增大,说明轴压对岩石的承载能力起到削弱效果。另外,回弹转折点所对应的循环冲击次数随轴压的增大而减小,即回弹转折点前移(试样S3-1和S3-3、试样S3-4和S3-5)。在相同的轴压下,随着冲击气压的增大,试样S3-3和S3-5、试样S3-4和S3-6的最大应变曲线具有上移的变化趋势,表明最大应变随冲击气压的增大而增大,而回弹转折点对应的循环冲击次数随冲击气压的增大而减小。

    图  8  峰值应力与循环冲击次数的关系
    Figure  8.  Relationship between peak stress and cyclic impact times
    图  9  最大应变与循环冲击次数的关系
    Figure  9.  Relationship between the maximum strain and cyclic impact times

    普遍而言,多数砂岩试样在经历1~5次冲击后承载能力才会显著下降,说明动态冲击次数不超过某一阈值时,冲击并不会对岩石内部造成实质性损伤,甚至还会提高其动态强度。

    以试样S3-2、S3-4、S3-5及S3-6为例,绘制平均应变率随循环冲击次数的关系曲线,如图10所示。在一维应力波理论下,仅考虑试样与入射杆及透射杆之间应力波的一次透反射,有如下关系式[17]

    图  10  平均应变率与循环冲击次数的关系
    Figure  10.  Relationship between average strain rate and cyclic impact times
    F=ρscsρ0c0ρscs+ρ0c01<F<0
    (4)
    εr(t)=Fεi(t)
    (5)
    εt(t)=(1F2)εi(t)
    (6)

    式中:ρ0c0ρscs分别为压杆和试样的波阻抗,F为试件与入射杆之间的反射波系数。

    分析认为,在砂岩试样的前几次冲击过程中,由于压密作用,岩石内部微粒错动,内部孔隙与裂缝闭合,微裂纹的闭合量总体上大于冲击造成的损伤量[5]。在以压密作用为主导的情况下,试样变得更加致密,造成试样的波阻抗ρscs逐渐增大,反射系数F也逐渐增大,当入射波εi(t)相同时,反射波εr(t)逐渐减小。同时,由式(3)可知,平均应变率˙εs(t)与反射波εr(t)成正比,因此,多数情况下(试样S3-2、S3-4和S3-5),在前几次冲击过程中,试样的平均应变率逐渐减小,且减小幅度呈逐渐减小趋势,而后趋于稳定,说明压密效果显著,试样的力学性能得以增强。当试样经历几次冲击作用后,其内部裂隙达到最大闭合率,动态强度达到最大,继续冲击加载后,岩石内部开始衍生裂纹,损伤软化效果逐渐加剧,此时试样内部由冲击造成的损伤量大于微裂纹的闭合量,岩石的波阻抗减小。同样地,反射系数F减小,引起反射波εr(t)逐渐增强,试样的平均应变率随冲击次数的增加而增大,表明岩石内部正在进行裂缝的生成与扩展,力学性能降低,从而形成了如图10所示的平均应变率先减小后增大的变化趋势。

    对比试样S3-2、S3-4、S3-6可以发现,在相同的轴压(15 MPa)、不同的气压(0.8、1.0、1.2 MPa)下,首次冲击的平均应变率分别为18.7、12.9、11.38 s−1,而后随着冲击气压的增大,砂岩的平均应变率总体呈先减小后增大的趋势。

    图11显示了砂岩承受的累计循环冲击次数与冲击气压的变化关系。当轴压为15、20 MPa时,对不同冲击气压下试样的累计循环冲击次数进行一元二次多项式回归,二者呈上凸的二次抛物线关系:轴压不变时,砂岩承受的累计冲击次数随冲击气压的增大呈先增大后减小的变化趋势;而当冲击气压不变时,砂岩的累计循环冲击次数随轴压的增大而减小,与唐礼忠等[10]得出的结论一致。同样地,在图12所示的不同冲击气压下轴压对砂岩累计循环冲击次数的影响中可以看出,累计循环冲击次数随轴压的增大而减小;而轴压不变时,累计冲击次数随冲击气压的增大先增大后减小。由此可以大致预测砂岩在爆破冲击下的失稳条件,为深部岩石巷道开挖提供理论依据。

    图  11  冲击次数与冲击气压的关系
    Figure  11.  Relationship between impact times and impact pressure
    图  12  冲击次数与预加轴压的关系
    Figure  12.  Relationship between impact times and preloading axial pressure

    根据SHPB一维试验原理,由入射波σi(t)、反射波σr(t)和透射波σt(t)及能量守恒定律,通过下式对能量进行求解

    Ei=A0C0ρst0σ2i(t)dt
    (7)
    Er=A0C0ρst0σ2r(t)dt
    (8)
    Et=A0C0ρst0σ2t(t)dt
    (9)
    EA=EiErEt
    (10)
    eV=EAVs
    (11)

    式中:EiErEtEA分别为入射能、反射能、透射能和吸收能,J;ρs为试样密度,kg/m3eV为单位体积耗能,J/cm3Vs为试样体积,cm3

    图13为相同冲击气压、不同轴压下试样的单位体积耗能与循环冲击次数的关系。当轴压为5 MPa时,随着冲击次数的增加,砂岩由释能区向吸能区方向发展;整个循环冲击过程中,释放的能量先增大后减小,随后向吸收能量过渡。分析认为,前几次冲击时,试样内部微裂纹逐渐闭合,随着冲击次数的增加,微裂纹闭合速度减缓,释放的能量逐渐减小,并达到释放能量峰值;随着冲击的持续进行,原本已经闭合以及将要闭合的微裂纹开始劣化,加之新微裂纹产生,不断消耗试样储存的弹性能,且消耗弹性能的速度随着裂纹扩展速度的加快而加快,最终导致试样破坏;当试样释放的能量接近零时,预加轴压产生的弹性能无法满足试样破碎而弹射的弹性能,试样吸收能量补充弹性能以使碎岩弹出。对于预加轴压分别为15、20 MPa的试样S3-2和S3-3,单位体积耗能与循环冲击次数曲线的变化趋势与试样S3-1具有一致性,但就试样S3-3而言,直到循环冲击结束,始终未出现由释能区向吸能区的过渡现象,原因在于试样破碎所消耗的弹性能较少,无法使碎石发生弹射现象。由砂岩的单位体积耗能规律可知,地下深部矿柱失稳大概率是在静载与动载的相互作用下导致的,当静压一定时,动载的频繁扰动是主要诱导因素。

    图  13  单位体积耗能与循环冲击次数的关系
    Figure  13.  Relationship between energy consumption per volume and cyclic impact times

    砂岩试样在整个循环冲击直至破坏的过程中,大致经历3个发展阶段:压密阶段、内部裂纹发展阶段和加速破坏阶段[21]。(1) 压密阶段:此阶段为前几次冲击过程中岩石内部持续演化过程,由于砂岩实体材料内部微粒较大,孔隙较多,在预加轴压和冲击荷载的作用下,微粒与微粒之间发生错动,通过局部变形调整使孔隙填塞(无法完全填塞),在此阶段试样结构更加致密均匀,力学性能有所改善,抵抗外界冲击荷载的能力逐渐增强。(2) 内部裂纹发展阶段:当试块经历多次冲击后,内部裂隙达到最小,强度达到最大,当继续冲击,伴随着损伤微裂纹的发展,轴向裂纹发展尤为明显,此阶段岩石承受动态荷载的能力逐渐减弱,且趋势不断加快。(3) 加速破坏阶段:岩石内部裂纹发展到一定程度,破坏面沿着试块轴向劈裂加速进行,最终产生宏观破碎形态,试样基本失去承载能力。

    图14显示了试样在循环冲击作用下典型的破坏形态。在单轴单次或循环冲击试验中,黄油的耦合效果很好,有效地减小了端部效应带来的影响,使得试样的破坏模式为沿纵向面的剥离破坏。对于预加轴向静压的一维循环冲击加载,轴压将黄油挤出,试样与各端面之间的黄油所剩无几,端部效应并无明显改善。另外,在轴压加载过程中,试样内部原本的横向裂纹被挤压,轴向产生微裂纹,导致试样尺寸沿横向变大且不均匀,沿纵向变小(变化甚微);由于入/透射杆与试样的材质不同,试样横向尺寸的变化将产生较大的摩擦力,试样横向外侧变形被限制,从而成为端部效应的又一助长因素(由于轴压经透射杆持续施加于试样,因此该效应只对透射界面端部有较大的影响)。由此可以发现,破坏面是沿着试样轴向倾斜发展的,即试样入射面破坏较大,透射面破坏相对较小。

    图  14  循环冲击作用下砂岩的破坏形态
    Figure  14.  Failure morphology of sandstone under cyclic impact loads

    当入射杆的应力波到达入射界面时,入射杆与试样材质的差异以及试样内部孔隙、裂缝的随机分布使入射波发生透射和反射,反射波到达入射界面以及透射波到达透射界面又会进行下一阶段的反射和透射,如此循环往复,对试样造成动态疲劳破坏。试样被冲击压缩时,与空气接触的圆周面因位于临空面而不受横向约束,越往试样内部,横向约束越大,导致横向变形越发困难,从而形成了如图15所示的试样破坏。图15中共轭双曲破坏面形成的主要原因是端部的约束效应,另一个原因则是临空面的无约束作用,当最外侧(图15中区域Ⅰ)破碎剥落时(前一次破坏为后一次破坏提供补偿空间),由于应力波的速度远大于裂纹的张裂速度,导致试样内侧部分因快速压缩而与外侧将要剥落部分发生摩擦,即形成剪切力。因此,砂岩最终的宏观破碎是由张裂破坏与剪切破坏共同导致的。

    图  15  循环冲击作用下砂岩的破坏模式
    Figure  15.  Failure mode of sandstone under cyclic impact loads

    通过对硬岩巷道砂岩开展一维动静组合循环冲击加载试验,分析了砂岩的应力-应变曲线特征、冲击气压和预加轴压对循环冲击次数的影响、能量耗散规律以及砂岩的破坏机制,得出以下主要结论。

    (1) 循环冲击过程中,砂岩的应力-应变曲线出现整体回弹现象,即随着冲击次数的增加,砂岩的峰值应力先增大后减小,而峰值应变、最大应变以及平均应变率则呈现相反的变化趋势。

    (2) 在相同的冲击气压下,砂岩的峰值应力及承受的累计循环冲击次数随轴压的增大而减小;在相同的预加轴压下,砂岩的峰值应力随冲击气压的增大而增大,但承受的累计冲击次数却随冲击气压的增大先增大后减小。

    (3) 根据应力-应变曲线和平均应变率随循环冲击次数的变化规律,可将砂岩的整个循环冲击过程分为压密阶段、内部裂纹发展阶段以及加速破坏阶段,对应于释放能先增后减以及向吸收能过渡的发展趋势。

    (4) 一维动静组合循环冲击试验证明了轴压的施加对微元体的力学性能起到劣化作用,砂岩体会发生由张裂和剪切共同作用的宏观破坏模式。

  • 图  上地幔和过渡带在地幔岩模型(Pyrolite模型和Eclogite模型)中的矿物组成(根据文献[10]修改)

    Figure  1.  Mineral composition of the upper mantle and transition zone in pyrolite models including Pyrolytic model and Eclogite model (modified from Ref.[10])

    图  石榴子石矿物的晶体结构(红球表示氧原子)

    Figure  2.  Crystal structure of garnet (The red sphere is the oxygen atom)

    图  镁铝-铁铝榴石固溶体体弹模量随铁铝榴石摩尔分数的变化[27]

    Figure  3.  Bulk modulus vs.mole fraction of almandine for almandine-pyrope binary system[27]

    图  锰铝-铁铝榴石固溶体体弹模量随铁铝榴石摩尔分数的变化[42]

    Figure  4.  Bulk modulus vs.mole fraction of almandine for spessartine-almandine binary system[42]

    图  钙铝-钙铁榴石固溶体体弹模量随钙铝榴石摩尔分数的变化[44]

    Figure  5.  Bulk modulus vs.mole fraction of grossular for grossular-andradite binary system[44]

    图  镁铝-钙铝榴石固溶体体弹模量随钙铝榴石摩尔分数的变化[45]

    Figure  6.  Bulk modulus vs.mole fraction of grossular for pyrope-grossular binary system[45]

    图  玄武岩中地幔岩包体(a)、金伯利岩中地幔包体和金刚石包体(b)以及(超)高压变质带榴辉岩(c)中石榴子石的端元组成(涉及的参考文献:[52-55](玄武岩包体)、[56-59](金伯利岩和金刚石包体)、[60-63]((超)高压变质带榴辉岩))

    Figure  7.  Composition of garnets from the mantle-derived basalt xenolith (a), kimberlite and diamond xenolith (b) and (ultra) high pressure metamorphic (UHPM) eclogite (c)(Corresponding reference:[52-55](mantle-derived basalt xenolith), [56-59](kimberlite and diamond xenolith), [60-63]((ultra) high pressure metamorphic eclogite).)

    表  1  常温高压条件下不同组分石榴子石的弹性参数

    Table  1.   Elastic parameters of garnets with different chemical compositions at room temperature and high pressure

    Sample Composition V0/nm3 K0/GPa K0 Ref.
    Natural Prp 175 [19]
    Natural Prp 1.5093(3) 173.7(32) 4.0a [31]
    Natural Prp Prp67Alm20Grs11 1.5377(6) 179(3) 4.0a [32]
    Synthetic Prp Prp100 171(3) 1.8(7) [21]
    Synthetic Prp Prp100 1.5034(5) 175(1) 4.5(5) [23]
    Synthetic Prp Prp100 175a 3.3(10) [24]
    Synthetic Prp Prp100 1.5029(3) 171(2) 4.4(2) [25]
    Synthetic Prp Prp100 1.50615(16) 163.7(17) 6.4(4) [28]
    Synthetic Prp Prp100 1.5027 190(6) 5.45a [18]
    Synthetic Alm Alm100 175(7) 1.5(16) [21]
    Synthetic Alm Alm100 1.5286 185(3) 4.2(3) [25]
    Synthetic Alm Alm100 1.53352(1) 172.6(15) 5.8(5) [28]
    Synthetic Alm Alm100 1.5336 168(5) 5.45a [18]
    Natural Prp-Alm Prp60Alm31 1.5292 177(6) 5.45a [18]
    Natural Prp-Alm Prp22Alm72 1.5300 173(6) 5.45a [18]
    Synthetic Prp-Alm Prp83Alm17 1.511(1) 172(4) 4.3a [27]
    Synthetic Prp-Alm Prp54Alm46 1.515(2) 174(2) 4.3a [27]
    Synthetic Prp-Alm Prp30Alm70 1.526(1) 183(2) 4.3a [27]
    Synthetic Prp-Alm Prp60Alm40 1.51632(13) 167.2(17) 5.6(5) [28]
    Natural Spe Spe93 1.5730 171(1) 5.4(2) [33]
    Synthetic Spe Spe100 171.8a 7.4 (10) [24]
    Synthetic Spe Spe100 1.5636 183(4) 5.1(6) [25]
    Natural Grs Grs90 1.6644 173(2) 4.25a [20]
    Natural Grs Grs97 1.6623 175(4) 4.25a [20]
    Synthetic Grs Grs100 168(25) 6.2(4) [29]
    Synthetic Grs Grs100 1.6602 170(4) 5.2(6) [25]
    Natural And And100 1.7476(5) 159(2) 4.0a [32]
    Synthetic And And100 1.7513 162(5) 4.4(7) [25]
    Natural Grs-And Grs14And84 1.6848(3) 166(2) 4.0a [26]
    Natural Grs-And Grs34And64 1.6909(4) 168(3) 4.0a [26]
    Natural Grs-And Grs63And34 1.6992(5) 173(2) 4.0a [26]
    Natural Uva Uva62Grs35 1.6975 160(1) 5.8(1) [33]
    Synthetic Uva Uva100 162a 4.7(7) [24]
    Synthetic Kat Kat100 66(4) 4.1(5) [29]
    Synthetic Maj Maj100 1.5131 161.2 4.0a [30]
    Synthetic Maj 1.5470(3) 164.8(34) 4.0a [31]
    Synthetic My-Na Maj 1.5054(2) 175.1(13) 4.0a [31]
    Synthetic Na-Maj 1.4855(3) 191.5(25) 4.0a [31]
    Synthetic Ca-Maj 1.5246(5) 169.3(34) 4.0a [31]
      Note: The superscript “a” represents the value was fixed in the equation of state fitting.Prp, Alm, Spe, Grs, And, Uva, Kat, and Maj stand for pyrope, almandine, spessartine, grossular, andradite, uvarovite, katoite, and majorite, respectively.
    下载: 导出CSV

    表  2  高温高压下不同组分石榴子石的热弹性参数

    Table  2.   Thermoelastic parameters of garnets with different chemical compositions at high temperature and high pressure

    Sample Composition V0/nm3 K0/GPa K0 (K/T)p/
    (GPa·K-1)
    α0/
    (10-5K-1)
    Ref.
    Synthetic Prp Prp100 1.5031(5) 170(2) 5.0a -0.020(3) 2.30(20) [34]
    Synthetic Prp Prp100 1.5007(19) 164(9) 4.9(12) -0.024(13) 2.97(45) [35]
    Synthetic Alm Alm100 1.53105(7) 179(3) 4.0a -0.043(14) 2.60(50) [39]
    Natural Alm Alm86Prp7Spe7 1.5396(9) 177(2) 4.0a -0.032(16) 3.10(70) [41]
    Synthetic Spe Spe100 1.56496a 171(4) 5.3(8) -0.049(7) 2.46(54) [37]
    Natural Spe-Alm Spe38Alm62 1.5446(6) 180(4) 4.0(4) -0.028(5) 3.16(14) [42]
    Natural Spe-Alm Spe64Alm36 1.5577(9) 176(4) 4.0(5) -0.029(5) 3.04(16) [42]
    Synthetic Grs Grs100 1.6630(10) 159.7(4.0) 5.10(48) -0.021(2) 2.77(24) [36]
    Natural Grs Grs97Alm3 1.66608a 168.2(1.7) 4.0a -0.016(3) 2.78(2) [40]
    Natural And And99 1.75405a 158.0(1.5) 4.0a -0.020(3) 3.16(2) [40]
    Synthetic Grs-And Grs50And50 1.7069(2) 164(2) 4.7(2) -0.018(2) 2.94(7) [44]
    Synthetic Prp-Grs Prp80Grs20 1.5394(2) 159.1(2) 4.4a 2.382(11) [45]
    Synthetic Prp-Grs Prp60Grs40 1.5784(2) 161.8(1) 4.4a 2.425(4) [45]
    Synthetic Prp-Grs Prp40Grs60 1.6120(2) 160.7(1) 4.4a 2.258(1) [45]
    Synthetic Prp-Grs Prp20Grs80 1.6384(2) 158.3(1) 4.4a 2.129(33) [45]
    Synthetic Uva Uva100 1.7368(8) 162(3) 4.3(4) -0.021(4) 2.72(14) [43]
    Synthetic hydrous Prp Prp100 1.5054(3) 162(1) 4.9(2) -0.018(4) 3.20(10) [46]
    Synthetic Na-Maj Na-Maj 1.47588 184(4) 3.8(6) -0.023(5) 3.22(20) [38]
      Note: The superscript “a” represents the value was fixed in the equation of state fitting.
    下载: 导出CSV
  • [1] 谢鸿森.地球深部物质科学导论[M].北京:科学出版社, 1997.
    [2] 谢鸿森, 侯渭, 周文戈, 等.地球深部物质科学——在静高压大腔体实验研究方面的某些进展[J].地学前缘, 2000, 7(1):217-228. http://d.wanfangdata.com.cn/periodical_dxqy200001021.aspx

    XIE H S, HOU W, ZHOU W G, et al.On material science of the Earth's interior[J]. Earth Science Frontiers, 2000, 7(1):217-228. http://d.wanfangdata.com.cn/periodical_dxqy200001021.aspx
    [3] 杨晓志.浅谈高温高压实验地球科学:方法和应用[J].矿物岩石地球化学通报, 2015, 34(3):509-525. http://www.cqvip.com/QK/84215X/201503/665661852.html

    YANG X Z.A brief introduction of high temperature and high pressure experimental geosciences:methods and advances[J]. Bulletin of Mineralogy, Petrology and Geochemistry, 2015, 34(3):509-525. http://www.cqvip.com/QK/84215X/201503/665661852.html
    [4] RINGWOOD A E.Composition and petrology of the Earth's mantle[M]. New York:McGraw-Hill, 1975.
    [5] ANDERSON D L, BASS J D.Mineralogy and composition of the upper mantle[J]. Geophysical Research Letters, 1984, 11(7):637-640. doi: 10.1029/GL011i007p00637
    [6] DUFFY T S, ANDERSON D L.Seismic velocities in mantle minerals and the mineralogy of the upper mantle[J]. Journal of Geophysical Research:Solid Earth, 1989, 94(B2):1895-1912. doi: 10.1029/JB094iB02p01895
    [7] RINGWOOD A E.Phase transformations and their bearing on the constitution and dynamics of the mantle[J]. Geochimica et Cosmochimica Acta, 1991, 55(8):2083-2110. doi: 10.1016/0016-7037(91)90090-R
    [8] ITA J, STIXRUDE L.Petrology, elasticity, and composition of the mantle transition zone[J]. Journal of Geophysical Research:Solid Earth, 1992, 97(B5):6849-6866. doi: 10.1029/92JB00068
    [9] IRIFUNE T, RINGWOOD A E.Phase transformations in subducted oceanic crust and buoyancy relationships at depths of 600-800 km in the mantle[J]. Earth and Planetary Science Letters, 1993, 117(1/2):101-110. https://www.sciencedirect.com/science/article/pii/0012821X9390120X
    [10] WOOD B J, KISEEVA E S, MATZEN A K.Garnet in the Earth's mantle[J]. Elements, 2013, 9(6):421-426. doi: 10.2113/gselements.9.6.421
    [11] MURAKAMI M, SINOGEIKIN S V, LITASOV K, et al.Single-crystal elasticity of iron-bearing majorite to 26 GPa:implications for seismic velocity structure of the mantle transition zone[J]. Earth and Planetary Science Letters, 2008, 274(3/4):339-345. https://www.sciencedirect.com/science/article/pii/S0012821X08004834
    [12] ANDERSON D L.Theory of the Earth[M]. Oxford:Blackwell Scientific Publications, 1989.
    [13] IRIFUNE T, SEKINE T, RINGWOOD A E, et al.The eclogite-garnetite transformation at high pressure and some geophysical implications[J]. Earth and Planetary Science Letters, 1986, 77(2):245-256. doi: 10.1016/0012-821X(86)90165-2
    [14] IRIFUNEA T, RINGWOOD A E. Phase transformations in primitive MORB and pyrolite compositions to 25 GPa and some geophysical implications[M]//MANGHNANI M, SYONO Y. High-Pressure Research in Mineral Physics: A Volume in Honor of Syun-iti Akimoto. Washington D C: American Geophysical Union, 1987: 231-242.
    [15] LIOU J G, ERNST W G, ZHANG R Y, et al.Ultrahigh-pressure minerals and metamorphic terranes-the view from China[J]. Journal of Asian Earth Sciences, 2009, 35(3):199-231. http://www.sciencedirect.com/science/article/pii/S1367912008001648
    [16] MEAGHER E P. Silicate garnet[M]//RABBE P H. Ortho-Silicates: Review in Mineralogy. Washington D C: Mineralogical Society of America, 1982: 25-66.
    [17] 李胜荣.结晶学与矿物学[M].北京:地质出版社, 2008.
    [18] TAKAHASHI T, LIU L G.Compression of ferromagnesian garnets and the effect of solid solutions on the bulk modulus[J]. Journal of Geophysical Research, 1970, 75(29):5757-5766. doi: 10.1029/JB075i029p05757
    [19] DUBA A, OLINGER B.Compression of garnet to 100 kilobars[J]. Journal of Geophysical Research, 1972, 77(14):2496-2499. doi: 10.1029/JB077i014p02496
    [20] WEAVER J S, TAKAHASHI T, BASS J.Isothermal compression of grossular garnets to 250 kbar and the effect of calcium on the bulk modulus[J]. Journal of Geophysical Research, 1976, 81(14):2475-2482. doi: 10.1029/JB081i014p02475
    [21] SATO Y, AKAOGI M, AKIMOTO S I.Hydrostatic compression of the synthetic garnets pyrope and almandine[J]. Journal of Geophysical Research:Solid Earth, 1978, 83(B1):335-338. doi: 10.1029/JB083iB01p00335
    [22] HAZEN R M, FINGER L W.Crystal structures and compressibilities of pyrope and grossular to 60 kbar[J]. American Mineralogist, 1978, 63(3/4):297-303. doi: 10.1007/BF00199500
    [23] LEVIEN L, PREWITT C T, WEIDNER D J.Compression of pyrope[J]. American Mineralogist, 1979, 64(7/8):805-808. http://www.researchgate.net/publication/285076500_Compression_of_pyrope?ev=auth_pub
    [24] LEGER J M, REDON A M, CHATEAU C.Compressions of synthetic pyrope, spessartine and uvarovite garnets up to 25 GPa[J]. Physics and Chemistry of Minerals, 1990, 17(2):161-167. doi: 10.1007/BF00199668
    [25] ZHANG L, AHSBAHS H, KUTOGLU A, et al.Single-crystal hydrostatic compression of synthetic pyrope, almandine, spessartine, grossular and andradite garnets at high pressures[J]. Physics and Chemistry of Minerals, 1999, 27(1):52-58. doi: 10.1007/s002690050240
    [26] FAN D W, WEI S Y, LIU J, et al.High pressure X-ray diffraction study of a grossular-andradite solid solution and the bulk modulus variation along this solid solution[J]. Chinese Physics Letters, 2011, 28(7):076101. doi: 10.1088/0256-307X/28/7/076101
    [27] HUANG S, CHEN J H.Equation of state of pyrope-almandine solid solution measured using a diamond anvil cell and in situ synchrotron X-ray diffraction[J]. Physics of the Earth and Planetary Interiors, 2014, 228:88-91. doi: 10.1016/j.pepi.2014.01.014
    [28] MILANI S, NESTOLA F, ALVARO M, et al.Diamond-garnet geobarometry:the role of garnet compressibility and expansivity[J]. Lithos, 2015, 227:140-147. doi: 10.1016/j.lithos.2015.03.017
    [29] OLIJNYK H, PARIS E, GEIGER C A, et al.Compressional study of katoite[Ca3Al2(O4H4)3] and grossular garnet[J]. Journal of Geophysical Research:Solid Earth, 1991, 96(B9):14313-14318. doi: 10.1029/91JB01180
    [30] YAGI T, UCHIYAMA Y, AKAOGI M, et al.Isothermal compression curve of MgSiO3 tetragonal garnet[J]. Physics of the Earth and Planetary Interiors, 1992, 74(1/2):1-7. https://www.sciencedirect.com/science/article/pii/S0031920197000290
    [31] HAZEN R M, DOWNS R T, CONRAD P G, et al.Comparative compressibilities of majorite-type garnets[J]. Physics and Chemistry of Minerals, 1994, 21(5):344-349. doi: 10.1007/BF00202099.pdf
    [32] HAZEN R M, FINGER L W.High-pressure crystal chemistry of andradite and pyrope:revised procedures for high-pressure diffraction experiments[J]. American Mineralogist, 1989, 74(3/4):352-359. https://hazen.carnegiescience.edu/sites/hazen.gl.ciw.edu/files/133-Andradite-1989.pdf
    [33] DIELLA V, SANI A, LEVY D, et al.High-pressure synchrotron X-ray diffraction study of spessartine and uvarovite:a comparison between different equation of state models[J]. American Mineralogist, 2004, 89(2/3):371-376. https://www.researchgate.net/publication/229430893_High-pressure_synchrotron_X-ray_diffraction_study_of_spessartine_and_uvarovite_A_comparison_between_different_equation_of_state_models
    [34] WANG Y B, WEIDNER D J, ZHANG J Z, et al.Thermal equation of state of garnets along the pyrope-majorite join[J]. Physics of the Earth and Planetary Interiors, 1998, 105(1/2):59-71. https://www.sciencedirect.com/science/article/pii/S0031920197000721
    [35] ZOU Y, GRÉAUX S, IRIFUNE T, et al.Thermal equation of state of Mg3Al2Si3O12 pyrope garnet up to 19 GPa and 1 700 K[J]. Physics and Chemistry of Minerals, 2012, 39(7):589-598. doi: 10.1007/s00269-012-0514-z
    [36] GRÉAUX S, KONO Y, NISHIYAMA N, et al.P-V-T equation of state of Ca3Al2Si3O12 grossular garnet[J]. Physics and Chemistry of Minerals, 2011, 38(2):85-94. doi: 10.1007/s00269-010-0384-1
    [37] GRÉAUX S, YAMADA A.P-V-T equation of state of Mn3Al2Si3O12 spessartine garnet[J]. Physics and Chemistry of Minerals, 2014, 41(2):141-149. doi: 10.1007/s00269-013-0632-2
    [38] DYMSHITS A M, LITASOV K D, SHATSKIY A, et al.P-V-T equation of state of Na-majorite to 21 GPa and 1 673 K[J]. Physics of the Earth and Planetary Interiors, 2014, 227:68-75. doi: 10.1016/j.pepi.2013.11.005
    [39] ARIMOTO T, GRÉAUX S, IRIFUNE T, et al.Sound velocities of Fe3Al2Si3O12 almandine up to 19 GPa and 1 700 K[J]. Physics of the Earth and Planetary Interiors, 2015, 246:1-8. doi: 10.1016/j.pepi.2015.06.004
    [40] PAVESE A, DIELLA V, PISCHEDDA V, et al.Pressure-volume-temperature equation of state of andradite and grossular, by high-pressure and-temperature powder diffraction[J]. Physics and Chemistry of Minerals, 2001, 28(4):242-248. doi: 10.1007/s002690000144
    [41] FAN D W, ZHOU W G, LIU C Q, et al.The thermal equation of state of (Fe0.86Mg0.07Mn0.07)3Al2Si3O12 almandine[J]. Mineralogical Magazine, 2009, 73(1):95-102. doi: 10.1180/minmag.2009.073.1.95
    [42] FAN D W, XU J G, MA M N, et al.P-V-T equation of state of spessartine-almandine solid solution measured using a diamond anvil cell and in situ synchrotron X-ray diffraction[J]. Physics and Chemistry of Minerals, 2015, 42(1):63-72. doi: 10.1007/s00269-014-0700-2
    [43] FAN D W, XU J G, MA M N, et al.P-V-T equation of state of Ca3Cr2Si3O12 uvarovite garnet by using a diamond-anvil cell and in-situ synchrotron X-ray diffraction[J]. American Mineralogist, 2015, 100(2/3):588-597. https://www.researchgate.net/publication/259324214_High-Pressure_and_High-Temperature_Stability_and_Equation_of_State_of_Superhydrous_Phase_B
    [44] FAN D W, KUANG Y Q, XU J G, et al.Thermoelastic properties of grossular-andradite solid solution at high pressures and temperatures[J]. Physics and Chemistry of Minerals, 2017, 44(2):137-147. doi: 10.1007/s00269-016-0843-4
    [45] DU W, CLARK S M, WALKER D.Thermo-compression of pyrope-grossular garnet solid solutions:non-linear compositional dependence[J]. American Mineralogist, 2015, 100(1):215-222. doi: 10.2138/am-2015-4752
    [46] FAN D W, LU C, XU J G, et al.Effects of water on P-V-T equation of state of pyrope[J]. Physics of the Earth and Planetary Interiors, 2017, 267:9-18. doi: 10.1016/j.pepi.2017.03.005
    [47] LIU X, SHIEH S R, FLEET M E, et al.High-pressure study on lead fluorapatite[J]. American Mineralogist, 2008, 93(10):1581-1584. doi: 10.2138/am.2008.2816
    [48] ALLRED A L.Electronegativity values from thermochemical data[J]. Journal of Inorganic and Nuclear Chemistry, 1961, 17(3/4):215-221. https://www.sciencedirect.com/science/article/pii/0022190261801425
    [49] SHANNON R D.Revised effective ionic radii and systematic studies of interatomic distances in halides and chalcogenides[J]. Acta Crystallographica Section A, 1976, 32(5):751-767. doi: 10.1107/S0567739476001551
    [50] SANI A, QUARTIERI S, BOSCHERINI F, et al.Fe2+-O and Mn2+-O bonding and Fe2+-and Mn2+-vibrational properties in synthetic almandine-spessartine solid solutions:an X-ray absorption fine structure study[J]. European Journal of Mineralogy, 2004, 16(5):801-808. doi: 10.1127/0935-1221/2004/0016-0801
    [51] GRIGGS D T, BLACIC J D.Quartz:anomalous weakness of synthetic crystals[J]. Science, 1965, 147(3655):292-295. doi: 10.1126/science.147.3655.292
    [52] CHEN S H, O'REILLY S Y, ZHOU X H, et al.Thermal and petrological structure of the lithosphere beneath Hannuoba, Sino-Korean Craton, China:evidence from xenoliths[J]. Lithos, 2001, 56(4):267-301. doi: 10.1016/S0024-4937(00)00065-7
    [53] XU W, LIU X, WANG Q, et al.Garnet exsolution in garnet clinopyroxenite and clinopyroxenite xenoliths in early Cretaceous intrusions from the Xuzhou region, eastern China[J]. Mineralogical Magazine, 2004, 68(3):443-453. doi: 10.1180/0026461046830198
    [54] IONOV D A, ASHCHEPKOV I, JAGOUTZ E.The provenance of fertile off-craton lithospheric mantle:Sr-Nd isotope and chemical composition of garnet and spinel peridotite xenoliths from Vitim, Siberia[J]. Chemical Geology, 2005, 217(1):41-75. https://www.sciencedirect.com/science/article/pii/S000925410400484X
    [55] HUANG X L, XU Y G, LO C H, et al.Exsolution lamellae in a clinopyroxene megacryst aggregate from Cenozoic basalt, Leizhou Peninsula, South China:petrography and chemical evolution[J]. Contributions to Mineralogy and Petrology, 2007, 154(6):691-705. doi: 10.1007/s00410-007-0218-4
    [56] ZHANG H F, ZHOU M F, SUN M, et al.The origin of Mengyin and Fuxian diamondiferous kimberlites from the North China Craton:implication for Palaeozoic subducted oceanic slab-mantle interaction[J]. Journal of Asian Earth Sciences, 2010, 37(5):425-437. https://www.sciencedirect.com/science/article/pii/S1367912009002478
    [57] ALIFIROVA T A, POKHILENKO L N, KORSAKOV A V.Apatite, SiO2, rutile and orthopyroxene precipitates in minerals of eclogite xenoliths from Yakutian kimberlites, Russia[J]. Lithos, 2015, 226:31-49. doi: 10.1016/j.lithos.2015.01.020
    [58] SPETSIUS Z V, BOGUSH I N, KOVALCHUK O E.FTIR mapping of diamond plates of eclogitic and peridotitic xenoliths from the Nyurbinskaya pipe, Yakutia:genetic implications[J]. Russian Geology and Geophysics, 2015, 56(1/2):344-353. https://www.sciencedirect.com/science/article/pii/S1068797115000267
    [59] RICHES A J V, ICKERT R B, PEARSON D G, et al.In situ oxygen-isotope, major-, and trace-element constraints on the metasomatic modification and crustal origin of a diamondiferous eclogite from Roberts Victor, Kaapvaal Craton[J]. Geochimica et Cosmochimica Acta, 2016, 174:345-359. doi: 10.1016/j.gca.2015.11.028
    [60] RUBATTO D, HERMANN J.Zircon formation during fluid circulation in eclogites (Monviso, Western Alps):implications for Zr and Hf budget in subduction zones[J]. Geochimica et Cosmochimica Acta, 2003, 67(12):2173-2187. doi: 10.1016/S0016-7037(02)01321-2
    [61] BUCHER K, FAZIS Y, CAPITANI C D, et al.Blueschists, eclogites, and decompression assemblages of the Zermatt-Saas ophiolite:high-pressure metamorphism of subducted Tethys lithosphere[J]. American Mineralogist, 2005, 90(5/6):821-835. https://pubs.geoscienceworld.org/msa/ammin/article-abstract/90/5-6/821/44411/blueschists-eclogites-and-decompression
    [62] GLODNY J, RING U, KUHN A, et al.Crystallization and very rapid exhumation of the youngest Alpine eclogites (Tauern Window, Eastern Alps) from Rb/Sr mineral assemblage analysis[J]. Contributions to Mineralogy and Petrology, 2005, 149(6):699-712. doi: 10.1007/s00410-005-0676-5
    [63] LIN W, SHI Y H, WANG Q C.Exhumation tectonics of the HP-UHP orogenic belt in Eastern China:new structural-petrological insights from the Tongcheng massif, Eastern Dabieshan[J]. Lithos, 2009, 109(3):285-303. http://www.sciencedirect.com/science/article/pii/S0024493708002314
    [64] 李晓东, 李晖, 李鹏善.同步辐射高压单晶衍射实验技术[J].物理学报, 2017, 66(3):130-142. http://www.oalib.com/paper/4339226

    LI X D, LI H, LI P S.High pressure single-crystal synchrotron X-ray diffraction technique[J]. Acta Physica Sinica, 2017, 66(3):130-142. http://www.oalib.com/paper/4339226
  • 加载中
图(7) / 表(2)
计量
  • 文章访问数:  9360
  • HTML全文浏览量:  3669
  • PDF下载量:  278
出版历程
  • 收稿日期:  2017-06-23
  • 修回日期:  2017-07-07

目录

/

返回文章
返回