镜质组反射率测试及其所反映的构造应力场.pdf

蒋建平,高广运,康继武.镜质组反射率测试及其所反映的构造应力场.地球物理学报, 2007 , 501 138～145 Jiang J P , Gao G Y, KangJ W. Tests on vitrinite reflectance of coal and analysis of tectonic stressfield.Chinese J.Geophys. in Chinese , 2007 , 501 138～145 镜质组反射率测试及其所反映的构造应力场 蒋建平 1 ,高广运 2 ,康继武 3 1上海海事大学交通运输学院,上海 200135 2同济大学地下建筑与工程系,上海 200092 3河南理工大学资源与环境工程系,焦作 454000 摘 要 为了利用煤镜质组反射率推断豫北构造应力场,基于构造应力是造成镜质组反射率各向异性VRA的主 要控制因素,对在豫北焦作、 安阳、 鹤壁三个矿区井下取的定向煤样品进行了室内显微镜下反射率测试,对测试数 据按有关公式进行计算后就得出了由最大反射率值、 中间反射率值、 最小反射率值为轴组成的镜质组反射率椭球 体VRI .结果表明,研究区煤样显示二轴晶正光性,它们经过了强烈的构造变形,研究区煤样镜质组反射率各向异 性是由构造应力所引起;由镜质组反射率椭球体VRI得出的豫北构造应力场与钻探、 节理统计、 河南省区域构造 应力场分析等方法得出的应力场基本上是吻合的,但也有一定的区别,主要体现在应力场到安阳、 鹤壁矿区后发生 了向左的偏转.总之,利用煤镜质组反射率各向异性推断构造应力场的方法是可行的、 是有效的. 关键词 镜质组反射率,构造应力场,煤,豫北 文章编号 0001 - 5733200701 - 0138 - 08中图分类号 P313收稿日期 2005 - 11 - 01 , 2006 - 10 - 29收修定稿 基金项目 上海市交通运输规划与管理重点学科建设项目和上海市教委科研项目05FZ03资助. 作者简介 蒋建平,男,1966年生,2002年6月于南京大学获博士学位,2004年5月于同济大学博士后出站.现在上海海事大学任教,主要从事 岩土工程、 港航工程和工程地质方面的教学和研究工作. E2mail jjpwx 1631com Tests on vitrinite reflectance of coal and analysis of tectonic stress field JIANGJian2Ping1, G AO Guang2Yun2, K ANGJi2Wu3 1College of Traffic Transportation,Shanghai Maritime University,Shanghai200135 ,China 2Department of Geotechnical Engineering,Tongji University,Shanghai200092 ,China 3College of Resource and Environment Engineering,Henan University of Science and Technology,Jiaozuo454000 ,China Abstract In order to deduce the tectonic stressfield in North Henan by coal vitrinite reflectance , based on the principle that the tectonic stress is the main controlling factor of the vitrinite reflectance anisotropy VRA , we have pered tests of coal vitrinite reflectance under microscope for directional coal specimen collected from Jiaozuo , Anyang , Hebi miningfield in Henan province , and gained vitrinite reflectance ellipsoids VRI made up of reflectance maximum , reflectance middle value and reflectance minimum through some calculation to testing data. The results show that coal specimen in studied area indicate indicatrix of optical biaxial crystal , which suffered from intense tectonic deation and the vitrinite reflectance anisotropy VRA is induced by tectonic stress. The tectonic stress field of North Henan obtained from vitrinite reflectance ellipsoids VRI is accordant basically with the tectonic stress field acquired from other s , such as drilling , joint statistics and Henan province territorial tectonic stress field analysis , but there is small difference that the tectonic stress field has left deflexion in Anyang and Hebi. In short , the deducing tectonic stress field by vitrinite reflectance anisotropy VRA is feasible and effective. Keywords Vitrinite reflectance , Tectonic stress field , Coal , North Henan 第50卷 第1期 2007年1月 地 球 物 理 学 报 CHINESE JOURNAL OF GEOPHYSICS Vol.50 , No.1 Jan. , 2007 1994-2007 China Academic Journal Electronic Publishing House. All rights reserved. 1 引 言 光、 波等在遇到界面时一般会产生反射或折 射 [1 ,2] .自然界中有些物质的反射率可作为标志性指 标来研究该物质的一些重要性能.如烧结矿中铁酸 钙的反射率与其组成结构、 冶金性能有相关关系 [3] ; 动物有机碎屑如笔石、 胞石、 沥青等的反射率可作 为烃源岩成熟度的指标 [4～7] ;镜质体反射率记录了 地层经历的温度历史,已成为研究盆地热史的主要 手段之一 [8] ;有机物镜质体反射率是目前国际上公 认的有机质成熟度指标 [9] . 煤层中的镜质组反射率是研究煤层变质程度的 有效指标 [10] ,它架起了研究煤变形与变质关系的桥 梁 [11 ,12] .由于煤层相对围岩的软弱性,构造煤镜质组 反射率的演化又受到构造应力环境的影响,它最完 整地记录了构造变形历史中应力作用和应变环境, 是进行构造煤所处构造环境研究的重要标志物. 20 世纪70年代末期以来,煤镜质组反射率各向异性分 析方法作为一种新手段,逐步应用于煤田构造研究, 在研究推覆构造特征、 煤反射率异常与基底断裂的 关系、 有限应变等方面 [13～16] ,已取得一定成果. 本文就是基于煤镜质组反射率对豫北地区的构 造应力场进行分析.构造应力及其应力场 [1～8 ,17～24] 不仅是基础地质研究的需要,也是力学中的重要研 究内容之一,特别是构造应力场是岩土力学或岩土 工程中的一个重要问题,它关系到各类工程的稳定 性和安全性,如在地下采矿工程、 隧道工程、 各类建 构筑物的地基基础工程中.所以本文的研究具有 重要的理论和实际意义. 2 取样与测试 211 取 样 研究区为河南省境内黄河以北、 京广线以西的 地区,包括焦作、 鹤壁、 安阳3个矿区.大地构造位置 处于华北古板块南缘带,是华北晚古生代聚煤盆地 的一部分.该地区地层发育比较齐全,二l煤所在的 煤系软岩位于二叠系下统山西组底部,二l煤平均厚 5 m ,属稳定煤层,煤种为瘦煤～无烟煤. 取样主要在上述的3个矿区二l煤中进行,取样 点分布如图1所示.取样时先用地质罗盘定向,取好 后将医用白胶布条贴于样品上,再在胶布上用箭头 标明方向及度数,同时在胶布上写明样品号及取样 地点,样品的详细情况还应在记录本上写清. 煤样全部采于井下巷道及工作面.采样点选择 在无明显断层和褶曲的区段,并要选在煤体破坏最 严重的部位,以反映取样点煤体镜质组最大反射率. 图1 研究区取样点分布 Fig.1 Locations of sampling sites in the studied area 212 测 试 1制煤砖方法 先将煤标本煮胶后切制成3～4 cm的立方体 煤砖 , 共做了18块煤砖.选择一个面平行于取样 时已定的方向,另两个面与之垂直,将3个面抛光 这样就得到测试光片54块,浸油后可放到显微镜 下测试 , 建立独立直角坐标系,以3个互相垂直截 面的3条交线代表3个坐标轴X,Y,Z 图2a . 2测量方法 用显微光度计分别顺制好的3个光面截面上 测量镜质组的油浸视最大反射率R 0 max、 视最小反 射率R 0 min及视最小反射率与坐标轴的夹角θ , 一 般每个光面上应测大于30个有效测点,然后分别计 算出各自的平均值R 0 max,R 0 min, θ , 这样3个光面可 得9个参数参加下面的计算. 3计算方法 用前面测得的数据R 0 max,R 0 min, θ计算下列参数 图2b .YOZ面、ZOX面与XOY面类似 令 λ ′ x 1 λx , λ ′ y 1 λy , λ ′ z 1 λz , γ ′ γ λ γ λ ′; λ ′ x cos 2 θx R0max 2 xy sin 2 θx R0min 2 xy , λ ′ y cos 2 θy R0max 2 yz sin 2 θy R0min 2 yz , 931 1期蒋建平等镜质组反射率测试及其所反映的构造应力场 1994-2007 China Academic Journal Electronic Publishing House. All rights reserved. 图2 坐标系及煤砖正交截面内VR椭圆示意图 a坐标系及煤砖截面椭圆 ;b XOY反射率椭圆. Fig.2 Coordinates and coal brick VR ellipses a Coordinates and coal brick section ellipse ; bXOYreflectance ellipse. λ ′ z cos 2 θz R0max 2 zx sin 2 θz R0min 2 zx ; γ ′ xy 1 R0max 2 xy - 1 R0min 2 xy sinθxcosθx, γ ′ yz 1 R0max 2 yz - 1 R0min 2 yz sinθycosθy, γ ′ zx 1 R0max 2 zx - 1 R0min 2 zx sinθzcosθz. 式中λ代表平方长度比,γ代表二维应变椭圆的剪 应变. 由这些参数得到一个椭球方程 λ ′ xx 2 λ ′ yy 2 λ ′ zz 2 - 2 γ ′ xyxy- 2 γ ′ yzyz- 2 γ ′ zxzx 1. 解此方程可得到光率体椭球的主半轴长度 1 eiλi 3 3wi ,i 1 ,2 ,3e为伸长度 其中 w1 2acosα, w2 2acosα 2 3π , w3 2acosα- 2 3π ; a 2K1 3 , K1 J 2 1 3 -J2, cos3α 2K2 a 3 , K2J3- 1 3 J1J2 2 J1 3 3 ; J1 λ ′ x λ ′ y λ ′ z, J2 λ ′ x λ ′ y λ ′ y λ ′ z λ ′ z λ ′ x- γ ′ 2 xy- γ ′ 2 yz- γ ′ 2 zx, J3 λ ′ x λ ′ y λ ′ z- 2 γ ′ xy γ ′ yz γ ′ zx- λ ′ x γ ′ 2 yz - λ ′ y γ ′ 2 zx- λ ′ z γ ′ 2 xy. 这样,在3个光面椭圆相容的前提下,煤的反射 率光率体椭球的3个主半轴长度 1 ei,即真反射 率就可求得,即椭球体形态已确定. 213 镜质组反射率测试结果 将测得的研究区各煤样的反射率光率体椭球的 3个 主 半 轴 长 度表 1 转 换 成 有 关 参 数 K ln 1 e1 1 e2 ln 1 e2 1 e3 ,投于Ramsy [25] 改进的弗林 图解上图3 ,可以看到,研究区的煤样点除一点外 都落在K 1直线的上方,即K 1 ,亦即 1 e1 - 1 e2 1 e2 - 1 e3或Ra-RbRb-Rc Ra、Rb、Rc分别为最大、 中间、 最小反射率轴 , 根据前面的分析及Chandra理论 [26] 知,研究区煤样 显示二轴晶正光性,它们经过了强烈的构造变形. 图3 研究区镜质组反射率椭球体弗林图解 Fig.3 Fulin diagram of VRI in the studied area K 1的点,即在K 1线上方,其反射率光率 041 地 球 物 理 学 报Chinese J. Geophys.50卷 1994-2007 China Academic Journal Electronic Publishing House. All rights reserved. 体为长球形椭球体;在K 1线下方的点属于扁球 形椭球体. 由表1中K值大小和双反射率大小可知研究 区煤样在构造应力作用下变形程度的相对大小. 表1 研究区镜质组反射率光率体参数表 Table 1 VRI parameters in the studied area 矿区矿点煤样类型 最大反射率 Rmax 最小反射率 Rmin 双反射率 Rmax-Rmin ln 1 e1 1 e2 ln 1 e2 1 e3 K 1 冯营矿无烟煤316512131811333018840123331794 焦作矿区 2 冯营矿无烟煤313752102611349019550130131173 3 九里山矿无烟煤410232160811415019930128431496 4 韩王矿无烟煤316092131401755016270125421468 安阳矿区5 大众矿贫煤217601165211108015380165901816 鹤壁矿区6 六 矿瘦煤213781141001968016130133911808 3 用镜质组反射率判定构造应力场 图4 镜质组反射率各向异性与应力应变的关系对应图解 Fig.4 Corresponding connection between VRA and stress , strain 311 镜质组反射率与构造应力相关的原理 煤镜质组反射率与煤的分子结构密切相关,反 射率随着芳香稠环缩合程度的增加和稠环中碳原子 数的增多而增大 [10] ,即芳香核叠片的排列由无序向 有序转变,镜质组反射率各向异性VRA明显地反 映了这一转变.镜质组反射率各向异性VRA是指 芳香环叠片逐渐采取定向排列,是多种因素静压 力、 构造应力、 温度、 时间综合作用的结果,但起主 要作用的则是构造应力中的压力因素,构造应力的 作用是控制反射率演化的主要因素 [27] . 研究表明,剪切应力可以促使石墨晶格形成并 加速石墨化,表现为VRA增强;在足够强度的剪应 力作用下,石墨结构甚至可形成于低煤级阶段;在同 一煤层中,构造变形煤分层和非构造煤分层在镜煤 反射率特征方面存在差异,反映了应力促使VRA增 强的基本规律. 定向压力给芳香环叠片发育施加一个优势方 位,而温度只能增加分子排列的随机性,因此,定向 压力是VRA定向性的根本原因. Levine et al. [28] 指 出,正如应力作用下晶体可以重新定向一样,定向压 力促使芳香环叠片重新排列,新的芳环层平行于最 小挤压应力方向发育,导致煤最大反射率Rmax趋于 最小挤压应力方向,煤最小反射率Rmin趋于垂直最 小挤压应力方向,最终形成的镜煤反射率椭球体 VRI型式代表了煤化进程中增量应变的总和. 根据VRA成因的理论,可以建立应力、 应变方 位与主反射率的对应关系图4 . 从而,不计煤化作用进程中镜质组反射率的绝 对增长,仅考虑反射率相对各向异性以VRI三主轴 的轴率表示 , 可以将VRI近似比于有限应变椭球体, 这就为利用VRA开展有限应变分析提供了可能. 在煤化作用过程中,地质应力包括静压力和构 造应力,两者都可以造成煤光率体椭球形态的变化. Levine et al. [28] 通过对美国加利福尼亚州Broad T op 煤田的反射率研究认为,在构造变形很弱或不变形 的地区以水平岩层的静压力作用为主,最大反射率 Ra在层面上,且各个方向近于相等,而最小反射 率Rc与层面垂直.镜质组的光性与矿物的一轴晶 负光性相似Ra-RbRb-Rc , 呈三轴不等的椭 球体,并且最小光率轴Rc与最大压应力方向平 行.在构造变形较弱的地区,如Broad T op煤田东部, 煤的光性显示二轴负光性.由此可看出,煤的光性从 一轴负光性、 二轴负光性、 二轴正光性依次反映煤的 构造变形从很弱或无到逐渐增强的变化趋势. 因此,地质应力可以促使煤中的芳环叠片重新 排列,以致光轴c始终保持着与最大压应力方向一 致的方位.镜质组反射率可以由各向同性向一轴负 光性或二轴光率体变化,二轴光率体的镜质组反射 率为一个三轴不等的椭球体.其轴率及主轴方位与 构造作用密切相关,反映了构造应力的相对大小和 受力的大致方位,表明煤的镜质组反射率椭球体也 是一种可以用来进行有限应变分析的应变标志物. 312 构造应力场的判定 由前面分析知,研究区煤样的变形是由构造应 力引起的.因此,由镜质组反射率各向异性VRA即 镜质组反射率椭球体与应力应变的关系,推出了研 究区的构造应力场表2和图5 .图5中的黑色椭圆 为反射率椭球体按图4得出的水平投影.由图5可 知,研究区煤受到了先期的挤压变形和后期的拉伸 变形,其中焦作矿区和安阳、 鹤壁矿区的有区别. 表2 研究区反射率椭球体三轴的方位 Table 2 Three2axis directions of VRI in the studied area 矿区矿 点Ra方位即最大主应力σ1方位Rc方位即最小主应力σ3方位Rb方位即中间主应力σ2方位 冯营矿32 ∠13125 ∠23221 ∠54 焦作矿区 冯营矿38 ∠11134 ∠15230 ∠45 韩王矿72 ∠21155 ∠16252 ∠57 九里山矿58 ∠17165 ∠8211 ∠59 安阳矿区大众矿328 ∠12234 ∠22140 ∠48 鹤壁矿区六 矿19 ∠8144 ∠19322 ∠73 图5 镜质组反射率椭球体反映的构造应力场 a焦作矿区 ;b 安阳矿区、 鹤壁矿区. Fig.5 Tectonic stress field derived from VRI of coal a Jiaozuo mines; b Anyang , Hebi mines. 图6为作者拍摄的煤砖显微照片 [11] ,照片显示 沿褶曲的轴部有一张性断层,表明该煤砖是先受到 了挤压作用,然后又受到拉伸作用. 河南煤田地质勘探公司 [29] 从节理统计、 河南省 区域构造应力场、 钻探等方面的研究表明,豫北存在 两期构造变形.先期存在向NW逆冲的挤压变形;后 期为伸展变形. 1节理统计表明太行区存在两套共轭剪节理 图7 .第1套共轭剪节理一组NNW - NNE向,左 型剪切;一组NWW向,右型剪切,最大主压应力指 示NW向,平均方位317 ∠3.第2套共轭剪节理 一组NNE向,右型剪切;一组NEE - SEE向,左型剪 切,最大主压应力指示NE向,平均方位52 ∠5. 第1套共轭剪节理表明,太行区在燕山早期发 育一期NW - SE向的挤压构造变形.第2套共轭剪 节理表明,太行区在燕山晚期受到一期NE向的伸 241 地 球 物 理 学 报Chinese J. Geophys.50卷 1994-2007 China Academic Journal Electronic Publishing House. All rights reserved. 图6 煤样显微构造照片放大250倍 Fig.6 Microscope photograph of coal sample 250 图7 节理统计分析图 Fig.7 Joint statistics chart 展构造变形. 2河南省区域构造应力场方面的研究表明,河 南省大地构造演化经历了8个阶段.其中第5阶段 和第6阶段为影响研究区煤层的两个重要阶段,第 7阶段的影响已不明显. 第5阶段即燕山构造运动早期.板块间左型走 滑运动在板内形成NW向具有走滑平移性质的断裂 带,这是区域内NWW - SEE挤压形成的构造形迹. 豫北也受到NWW - SEE向挤压,应力场测量结果为 最大主压应力方向NW ,产状317 ∠3.此阶段为第 1期二1煤层变形期图8 . 第6阶段即燕山构造运动晚期.是陆壳强烈压 缩期,表现为华北板块向南仰冲,华南板块向北俯 冲,河南区域应力场是近SN向挤压.但由于河南板 内动力来源于南部边界,应力传递由南向北衰减,应 力场测量结果为最大主压应力方向NE,产状52 ∠5.此阶段为第2期二1煤层变形期图9 . 3如图10所示,罗1井的钻探表明,中原油田 在聊兰断裂NNE向下存在一个逆冲断裂系,走向 和聊兰断裂平行,倾向SE,向NW向逆冲,称作新罗 逆断层.从断层交接关系看,聊兰断裂切割了新罗逆 图8 河南省燕山早期构造应力场[29] Fig.8 Tectonic stress field of Henan province in the early Yanshan time[29] 图9 河南省燕山晚期构造应力场[29] Fig.9 Tectonic stress field of Henan province in the late Yanshan time[29] 断层,说明在NE向拉伸断裂形成之前,存在一期向 341 1期蒋建平等镜质组反射率测试及其所反映的构造应力场 1994-2007 China Academic Journal Electronic Publishing House. All rights reserved. 图10 罗1井钻探剖面图[29] N -上第三系Tertiary;K-白垩系Cretaceous;O -奥陶系Ordovician; P -二叠系Permian;Q -第四系Quaternary;C-石炭系Carboniferous Fig.10 Cutaway view of Luo 1 drilling well[29] NW向逆冲的挤压变形. 综上所述,本文的研究结果和其他方法得到的 应力场基本上是吻合的,但也有一定的区别.区别主 要体现在安阳、 鹤壁矿区,说明构造应力场在豫北的 安阳、 鹤壁矿区发生了变化,即,在燕山早期,豫北的 焦作矿区受到的NWW - SEE向挤压,到豫北的安 阳、 鹤壁矿区已转化为NEE - SWW向挤压;在燕山 晚期,焦作矿区的NNE - SSW向拉伸,到安阳、 鹤壁 矿区已转化为NNW - SSE向拉伸. 4 结 论 1构造应力是造成镜质组反射率各向异性 VRA的主要控制因素. 2研究区煤样显示二轴晶正光性,它们经过了 强烈的构造变形,研究区煤样镜质组反射率各向异 性是由构造应力所引起. 3由镜质组反射率椭球体VRI得出的豫北构 造应力场与钻探、 节理统计、 河南省区域构造应力场 分析等方法得出的应力场基本上是吻合的,但也有 一定的区别,主要体现在应力场到安阳、 鹤壁矿区后 发生了向左的偏转,即,在燕山早期,豫北的焦作矿 区受到的NWW - SEE向挤压,到豫北的安阳、 鹤壁 矿区已转化为NEE - SWW向挤压;在燕山晚期,焦 作矿区的NNE - SSW向拉伸,到安阳、 鹤壁矿区已 转化为NNW - SSE向拉伸. 4煤镜质组反射率各向异性分析法推断构造 应力场是可行的、 有效的. 参考文献References [ 1 ] 孙成禹,杜世通.平面波在粗糙界面上的反射特征研究.地球 物理学报,2006 ,493 903～907 Sun C Y,Du S T. Reflection features of planar acoustic waves from rough interfaces.Chinese J.Geophys. in Chinese , 2006 ,493 903～907 [ 2 ] 刘福平,李瑞忠,杨长春等.非均匀电磁波在导电媒质界面反 射的横向偏移.地球物理学报,2006 ,492 533～539 Liu F P ,Li R Z, Yang C C, et al. The lateral shift of inhomogeneous reflected electromagnetic wave on the interface of conductive media. Chinese J.Geophys. in Chinese , 2006 ,492 533～539 [ 3 ] 许传智.烧结矿中铁酸钙的反射率与其组成结构、 冶金性能的 关系.武钢技术,1998 ,1 40～44 Xu C Z.Theconnectionbetweenreflectanceofferritein agglomeration and its structure , metallurgy perance.Wisco Technology, 1998 ,1 40～44 [ 4 ] Bertrand R , HeeouxY Chitinozoan.Graptoliteandscolecodont reflectance as an alternative to vitrinite and pyrobitumen reflectance in Ordovician and Silurian strata Anticosti Island , Quebec , Canada. Am.Assoc.Pet.Geol. , 1987 , 2 71～77 [ 5 ] 肖贤明,刘德汉,傅家谟.沥青反射率作为烃源岩成熟度指标 的意义.沉积学报,1991 ,9增 138～146 Xiao X M, Liu D H, Fu J M.The significance of Bitumen reflectanceasamatureparameterofsourcerocks.Acta Sedimentologica Sinicain Chinese , 1991 ,9 Suppl. 138～146 [ 6 ] 祝幼华,李永军.笔石反射率测试中的若干问题.大庆石油地 质与开发,1999 ,181 19～21 Zhu Y H,Li Y J.Some questions in graptolite reflectivity test. Petroleum Geology Oilfield Development in Daqingin Chinese , 1999 ,181 19～21 [ 7 ] 曾凡刚,程克明,吴朝东.应用海相镜质组反射率研究华北地 区下古生界成熟度.地质地球化学,1998 ,263 21～24 Zeng F G,Cheng KM, Wu C D. Maturity of the lower Palaeozoic in North China in terms of reflectivity of marine vitrinites.Geology Geochemistryin Chinese , 1998 ,263 21～24 [ 8 ] 王 玮,周祖翼,于 鹏.镜质体反射率与最高温度及其附近 温度变化率的关系 几种镜质体反射率计算模型的比较. 地球物理学报,2005 ,486 1375～1383 Wang W,Zhou Z Y, Yu P. Relations between vitrinite reflectance , peak temperature and its neighboring temperature variation rate A comparison of s.Chinese J.Geophys. in Chinese , 2005 , 486 1375～1383 [ 9 ] 刘祖发,肖贤明,傅家谟等.海相镜质组反射率用作早古生代 烃源岩成熟度指标研究.地球化学,1999 ,286 580～588 Liu Z F , Xiao XM, Fu J M, et al. Marine vitrinite reflectance as a maturity indicator of lower Palaeozoic hydrocarbon source rocks. Geochemicain Chinese ,1999 ,286 580～588 [10] 杨 起.煤地质学进展.北京科学出版社,1987 441 地 球 物 理 学 报Chinese J. Geophys.50卷 1994-2007 China Academic Journal Electronic Publishing House. All rights reserved. Yang Q. Evolvement of Coal Geology in Chinese . Beijing Science Press , 1987 [11] 蒋建平.豫北大地构造演化与二1煤变形变质研究[硕士论 文].焦作河南理工大学,199916 Jiang J P. Study on North Henan tectonic evolvement and 21coal ation and metamorphose [ Master’s thesis ] in Chinese . Jiaozuo Henan University of Science and Technology , 199616 [12] 蒋建平,罗国煜,康继武.豫北煤田焦作、 鹤壁、 安阳矿区二1煤 层优势断裂构造地球化学特征.地球化学,2000 ,294 367～375 Jiang J P , Luo G Y, Kang J W. Tectonochemical characteristics of preferred fault zones of the 21coal in Jiaozuo , Hebi , Anyang mining area, North Henan province.Geochemicain Chinese ,2000 ,29 4 367～375 [13] Stone IJ , Cock A C. The influence of some tectonic structures upon vitrinite reflectance.J.Geol. , 1979 , 875 479～508 [14] Hower J C, Davis A. Vitrinite reflectance anisotropy as a fabric element.Geology, 1981 , 94 165～168 [15] 曹代勇.安徽淮北煤田推覆构造中煤镜质组反射率各向异性 研究.地质论评,1990 ,364 333～340 Cao D Y. The vitrinite reflectance anisotropy in the nappe structure in theHuaibeicoalfield ,Anhuiprovince.GeologyReview in Chinese , 1990 ,364 333～340 [16] 王文侠.涟源坳陷煤反射率变化及其与深部断裂构造的关系. 煤田地质与勘探,1991 ,192 20～25 Wang W X. Coal reflectance changes and connection with deep faults in Lianyuancoalfield.MiningGeologyandExploration in Chinese , 1991 ,192 20～25 [17] 朱艾斓,徐锡维,周永胜等.川西地区小地震重新定位及其活 动构造意义.地球物理学报,2005 ,483 629～636 Zhu A L ,Xu X W, Zhou YS, et al. Relocation of small earthquake in western Sichuan , China and its implications for active tectonics. Chinese J.Geophys. in Chinese , 2005 ,483 629～636 [18] 李建忠,潘忠习,冯心涛等.聂拉木地区高喜马拉雅岩石磁组 构及其构造含义.地球物理学报,2006 ,492 496～503 Li J Z, Pan Z C, Feng X T, et al. Magnetic fabrics and tectonic implications of the Higher Himalayan rocks in Nyalam, southern Tibe