收稿日期: 2016-01-18
修回日期: 2016-03-15
网络出版日期: 2016-07-15
基金资助
国家自然科学基金重大项目（41086965）资助
Calculations of rheological properties of multi-mineral foliated rocks using finite element method based on the rheological properties of minerals
Received date: 2016-01-18
Revised date: 2016-03-15
Online published: 2016-07-15
高温高压下岩石表现出流变性质。流变参数是数值模拟地球动力过程必需的重要参数。获得流变参数需要借助于室内的高温高压蠕变实验以及野外的地球物理资料观测的反演。对单矿物的高温高压实验比较多，多矿物组成的岩石的实验则较难较少。本文探讨在单矿物流变参数已知的条件下，能否运用有限元法，计算多矿物岩石整体的流变参数。研究表明，计算结果与高温高压实验获得的参数可以对应。即使矿物力学性质可以近似为各向同性时，如果矿物排列有优势方向（例如存在面理的花岗质糜棱岩），计算表明岩石整体会显示各向异性。当加载外力与矿物面理方向垂直或平行时，岩石的等效黏滞系数较大；当所成角度为30°左右时，等效黏滞系数最小。
蒋龙腾 , 刘贵 , 石耀霖 . 利用矿物流变性质对面理岩石流变性质的有限元方法计算[J]. 中国科学院大学学报, 2016 , 33(4) : 513 -518 . DOI: 10.7523/j.issn.2095-6134.2016.04.012
Rocks show rheological properties under high temperature and high pressure.Rheological parameter is an important parameter for numerical simulation of the dynamic process of the earth.The rheological parameters can be obtained by means of the high temperature and high pressure creep experiments in the laboratory and by the inversion of the field geophysical data.There have been some high temperature and high pressure experiments on single minerals.However, it is difficult to carry out multi-mineral experiments.In this study, we discuss the possibility of using the finite element method to calculate the rheological parameters of multi-mineral rocks under the condition that the rheological parameters of minerals are known.The results show that the calculated results match the parameters obtained in the high temperature and high pressure experiments.The results show that the rock as a whole will show anisotropy if the mineral arrangement has the advantage of direction (such as the existence of the granitic mylonite foliation), even though mineral mechanical properties are approximatly isotropic.When loading force is perpendicular or parallel to mineral foliation, equivalent viscous coefficient of rock is larger; and when the angle is 30 degrees, equivalent viscous coefficient is minimized.
[1] 石耀霖,张贝,张斯奇.地震数值预报[J].物理,2013,42(4): 237-255.
[2] Liu G, Zhou Y, He C h, et al.An experimental study of effect of pre-existing fabric on deformation of foliated mylonite at high temperature and pressure[J].Geological Journal, 2014: 51(1):92-112.
[3] Burgmann R, Dresen G.Rheology of the lower crust and upper mantle: evidence from rock mechanics, geodesy, and field observations[J].Annual Review of Earth and Planetary Sciences, 2008,36: 531-567.
[4] 张晁军,曹建玲,石耀霖.从震后形变探讨青藏高原下地壳黏滞系数[J].中国科学 D辑:地球科学,2008,38(10):1 250-1 257.
[5] 刘善琪,李永兵,田会全,等.含湿孔隙岩石有效热导率的数值分析[J].地球物理学报,2012,55(12):4 239- 4 248.
[6] Tullis T, Horowitz F G, Tullis J.Flow laws of polyphase aggregates from end-member flow laws[J].Journal of Geophysical Research, 1991, 96(B5): 8 081-8 096.
[7] 崔晓佳,张怀,石耀霖.岩石宏观流变性质的数值模拟[J].岩石学报,2008,24(6):1 417-1 424.
[8] Jing H, Zhang H, Li H, et al.Parallel numerical analysis on the rheology of the martian ice-rock mixture[J].Journal of Earth Science, 2011, 22(2): 176-181.
[9] Offerhaus L J, Wirth R, Dresen G.High-temperature creep of polycrystalline albite[J].Deformation Mechanisms, Rheology and Tectonics, 2001, 124:107-131.
[10] Rutter E H, Brodie K H.Experimental intracrystalline plastic flow in hot-pressed synthetic quartzite prepared from Brazilian quartz crystals[J].Journal of Structural Geology, 2004, 26(2):259-270.
[11] Mariani E, Brodie K H, Rutter E H.Experimental deformation of muscovite shear zones at high temperatures under hydrothermal conditions and the strength of phyllosilicate-bearing faults in nature[J].Journal of Structural Geology, 2006, 28:1 569-1 587.
[12] 刘贵,石耀霖,周永胜.各向异性花岗岩的强度和变形机制实验研究[J].吉林大学学报:地球科学版,2015,45(7):1 506-1 514.
[13] 石耀霖,曹建玲.中国大陆岩石圈等效粘滞系数的计算和讨论[J].地学前缘,2008,15(3):82-95.
[14] 刘贵,周永胜,石耀霖.先存组构对各向异性岩石流变强度的影响[J].地震地质,2014,36(3):918-928.
/
| 〈 |
|
〉 |
