Measurement results of in-situ stress in Guyuan area of Ningxia on the southwest margin of Ordos block and its causation analysis
-
摘要:
鄂尔多斯地块西南缘新构造活动强烈、地震频发,具有复杂多样的构造变形模式和活动特征。为了解鄂尔多斯地块西南缘地壳浅部地应力分布规律及断层稳定性,利用宁夏固原地区的水压致裂地应力测量数据,结合其他实测及震源机制解资料,分析了鄂尔多斯地块西南缘构造应力场特征。结果表明:①研究区2个钻孔的主应力关系整体表现为SH > Sv > Sh,水平应力起主导作用,属于走滑型应力状态,钻孔附近最大水平主应力方位平均为N59°W,与震源机制解获得的青藏高原东北缘主压应力方位有差异,推断鄂尔多斯地块西南缘现今NWW向走滑剪切应力环境的形成可能主要受到海原断裂带和六盘山断裂带的影响,应为局部构造和区域构造应力场共同作用的结果。②利用Mohr-Coulomb准则及Byerlee定律,摩擦系数取0.6~1.0,对研究区的现今地应力状态分析后发现,鄂尔多斯地块西南缘海原断裂带和六盘山断裂带的地应力大小未达到地壳浅部断层产生滑动失稳的临界条件,处于较稳定的应力状态。该研究成果为鄂尔多斯地块关键构造部位的断裂活动性分析和地质环境安全评价提供了参考依据。
Abstract:The southwestern margin of the Ordos block has experienced many great earthquakes with complicated tectonic deformation and tectonic activities.In order to understand in-situ stress environment of southwestern margin of the Ordos block and to evaluate the seismic risk, the authors carried out the in-situ stress measurement by hydraulic fracturing in two boreholes in the Guyuan area of Ningxia.Combined with other measured data and focal mechanism solutions, the authors discussed the tectonic stress field characteristics of the southwestern margin of the Ordos block.Some conclusions have been reached:(1) The results show that the relationship of principal stress of the two boreholes in the study area is SH > Sv > Sh, belonging to the strike-slip stress state.The azimuth of maximum horizontal principal stress near the borehole is N59°W; however, it is different from the P-axis orientation of focal mechanism solutions on the northeastern margin of the Tibetan Plateau.It is considered that the formation of the NWW strike-slip shear stress environment on the southwestern margin of the Ordos block is mainly affected by the Haiyuan fault and Liupanshan fault zone.The current stress environment may be the result of the combined effect of local structure and regional tectonic stress fields.(2) The crustal activity in the study area is discussed using the Mohr-Coulomb criterion and Byerlee's law under the premise that the friction coefficient is 0.6~1.0.The in-situ stress state of the southwestern margin of the Ordos block has not reached the sliding critical condition of the shallow fault and is in a relatively stable crustal stress state.The results have great significance for analysis of active faults as well as assessment of the regional geological environment and geological disasters prevention.
-
自Shand 1927年[1]提出过铝质花岗岩的概念以来,许多学者[2-18]就其岩石学特征、岩石地球化学特征、构造环境、成因等进行过研究。随着大陆动力学研究的深入开展,运用过铝质花岗岩探讨大地构造背景成为研究热点。
内蒙古东乌珠穆沁旗(以下简称东乌旗)地处华北板块与西伯利亚板块之间的中亚-蒙古造山带中东段的兴蒙褶皱带,属于二连-贺根山基性-超基性岩带(即二连-贺根山板块对接带)和中蒙古-得尔布干深断裂之间西伯利亚板块东南缘的晚古生代安第斯型陆缘增生带[19]。该陆缘增生带内岩浆岩广布,构成二连-东乌旗晚古生代巨型花岗岩带,即查干敖包-奥尤特-朝不楞构造-岩浆岩带[20],其演化记录了华北板块与西伯利亚板块汇聚和古亚洲洋闭合的历史,是研究两大板块拼合及古亚洲洋最终关闭时限证据的重要载体[21]。近些年来,众多学者在内蒙古东乌旗一带开展研究工作,在原厘定的燕山期复合型岩体中发现众多华力西期岩体,对其开展了岩浆岩与大地构造背景讨论,对晚古生代西伯利亚板块与华北板块沿二连-贺根山缝合带的碰撞造山事件进行了时间上的约束[22-32]。前人研究认为,区内晚石炭世—早二叠世发育大量具有后造山特征的I型花岗岩类、I-S型混合特征花岗岩类、A型花岗岩类及碱性花岗岩类,少有报道晚石炭世—早二叠世具有典型后碰撞强过铝质S型花岗岩类。笔者在内蒙古东乌旗开展区域地质填图过程中,在区内发现查干哈达音亨嘎岩体,岩性主要为含石榴子石二云母二长花岗岩,具有强过铝质花岗岩特征。前人对该岩体未开展过相关研究。为此,本文通过野外和室内岩石学、岩石地球化学及同位素年代学的研究,确定其岩浆作用时代和成因,为深入认识东乌旗地区地壳演化历史和过程提供依据。
1. 地质背景及岩体地质
研究区内地层发育为晚泥盆世安格尔音乌拉组和上新世宝格达乌拉组(图 1)。前者为一套砂岩、粉砂岩和硅质岩夹灰岩透镜体组合的滨浅海相碎屑沉积地层,变形复杂,褶皱发育。研究区内发育北西向、北东向2组断裂,且多被后期不同时代的酸性侵入岩类侵入破坏,发生角岩化,多处可见残留顶垂体;后者为一套红色和绿色泥岩组合的湖相沉积,不整合于其形成之前的所有地质体之上。
图 1 内蒙古东乌珠穆沁旗地区大地构造位置和地质简图(角图据参考文献[33])Figure 1. Simplified geological map of the study area and tectonic position of the East Ujimqin Banner, Inner Mongolia研究区内侵入岩极发育,但岩性较简单,均为华里西期酸性花岗岩类。极少量闪长岩类、细粒花岗岩类呈脉状产出,按照接触关系及粒度可以分解为中细粒似斑状含黑云母二长花岗岩单元、细粒二长花岗岩单元和中细粒(似斑状)含石榴子石二云母二长花岗岩单元。
查干哈达音亨嘎岩体为中细粒似斑状含石榴子石二云母二长花岗岩单元,分布于内蒙古东乌旗查干哈达音享嘎一带,岩体出露面积约16.35km2,独立产出,周边多被第四系或上新统不整合覆盖,仅在西南侧可见极少量泥盆系安格尔音乌拉组二段(D3a2)呈顶垂体侵入接触,在接触带附近地层中变质粉砂岩具角岩化现象,未见同其他岩体单元接触界线。岩体内部岩性极简单,无岩相变化,以中细粒花岗岩为主,部分可见较大钾长石和石英斑晶,包体、脉体不发育。
2. 样品概况及测试方法
2.1 样品概况
本文同位素年龄样的采样坐标为北纬46°02′01″、东经117°21′30″,化学全分析样品的采样具体位置见图 1。样品岩性均为中细粒(似斑状)含石榴子石二云母二长花岗岩,呈浅灰色-灰白色,中细粒花岗结构,少许似斑状结构,块状构造(图 2),主要组成矿物为斜长石、钾长石、石英、黑云母、白云母,极少量石榴子石。斜长石呈半自形板状,杂乱分布,粒径一般为2.0~3.0mm,部分为0.2~2.0mm,轻绢云母化、高岭土化及局部白云母化,可见环带构造,聚片双晶发育,少数被钾长石呈净边状交代,约占45%;钾长石呈他形粒状,少数呈半自形宽板状,杂乱分布,粒径一般2.0~5.0mm,部分0.2~2.0mm,具高岭土化,晶内嵌布少量斜长石、石英等小包体,局部交代斜长石,占25%~30%;石英呈他形粒状,不均匀分布,粒径一般为0.2~2.0mm,少数为2.0~3.0mm,可见轻波状消光,集合体似堆状分布,约占25%;黑云母、白云母呈鳞片状-叶片状,星散状分布,粒径一般为2.0~3.0mm,部分为0.1~2.0mm,黑云母稍绿泥石化,白云母交代黑云母,占1%~5%;石榴子石呈近等轴粒状,星散状分布,粒度一般为0.1~0.5mm,云母沿其裂纹发生交代。岩石内还见少量铁质等充填的裂纹、裂隙。主要副矿物为锆石、石榴子石、独居石,另可见1~10粒毒砂、黄铁矿、方铅矿、磁铁矿颗粒。
2.2 测试方法
同位素年龄样品为(似斑状)含石榴子石二云母二长花岗岩,年代测试采用锆石U-Pb激光烧蚀同位素测年法,室内对样品进行切割去除风化表面后用粉碎机粉碎至120目,用于锆石分选。锆石的分选在河北省区域地质矿产调查研究所完成。对选取的锆石样品进行制靶、抛光和反射光、透射光、阴极发光照片拍摄,选择具有典型晶体特征、无包体、无裂纹的锆石用于U-Pb同位素年龄测试。上述工作均在西北大学大陆动力学国家重点实验室完成,所用阴极发光仪器为配有英国Gatan公司的MonoCL3+型阴极荧光探头的美国FEI公司Quan ta 400 FEG扫描电镜,U-Pb测年系统中激光剥蚀系统为GeoLas2005,ICP- MS为Agilent 7500a。年龄测试激光剥蚀过程中载气为氦气,激光频率为10Hz,强度80MJ,剥蚀束斑直径30μm,每个时间分辨分析数据为30s的空白信号和50s的样品信号。定年过程中采用NIST610硅玻璃标准优化仪器,锆石标准91500为外标进行同位素分馏校正,用锆石标准GJ-1观察仪器的状态及验证测试结果的精确度。测试分析中每分析5个样品点校正2次91500标样,并测试1次GJ-1标样,每10个样品标准化1次NIST610。最终测试数据的离线处理采用软件ICPMASDataCal[34]完成,U-Pb年龄谐和图绘制和年龄权重平均计算均采用Isoplot3.0[35]完成。
主量和微量元素含量测试选择天然露头下不同地区的新鲜样品,具体采样位置见图 1。室内对样品进行切割去除风化表面后用粉碎机粉碎至200目,采用酸溶法制成溶液用于岩石地球化学测试。测试分析在核工业北京地质研究院分析测试中心完成,主量元素使用XRF法测试(二价和三价铁由化学法测定),所用仪器为飞利浦PW2404X射线荧光光谱仪,测试误差优于5%;微量元素采用酸溶法,将制备好的样品溶液在ICPMS上测试,所用仪器为德国Finnigan MAT公司制造的HR-ICP-MS(Element I),工作温度、相对湿度分别为20℃和30%,微量元素含量大于10×10-6时的相对误差小于5%,小于10×10-6时的相对误差小于10%。
3. 测试结果
3.1 同位素年龄
用于测试的锆石在双目镜下颜色均为粉黄色,透明、金刚光泽,半自形柱状,晶体表面较光滑,晶棱晶面较平直完整,微有铁染,锆石结晶后改造痕迹不明显,但可见锥柱不对称的歪晶及固相黑包体,揭示锆石结晶时的介质环境不适宜按照理想形态生长。锆石粒径以0.01~0.1mm为主,个别为0.1~0.3mm,伸长系数以1.5~3为主,少量为3~5。全部锆石晶群集中,类型均为(100)型复柱岩浆锆石,具有较晚较低温度结晶的同源岩浆产物特征。锆石阴极发光图像(图 3)显示岩浆型锆石的振荡韵律环带结构或明暗相间的条带结构,且一般为较窄的岩浆环带,这种较窄的岩浆环带一般为低温条件下微量元素的扩散速度慢形成的[36],该结果与锆石的结晶类型特征一致。锆石的Th、U含量分别为158×10-6~ 563×10-6和250×10-6~2461×10-6,Th/U值为0.17~0.84,表现为典型的岩浆成因锆石[36-38]。
U-Th-Pb法年龄测定可以同时获得4个年龄值,如果这4个值较接近,其算术平均值即为一致年龄,代表矿物结晶年龄。由于U、Pb的活动性较强,而Th4+的地球化学性质与U4+相似,已形成的岩石和矿物难免受到后期地质作用的影响,造成母、子体核素不同程度丢失(或获得),破坏了体系的封闭性,导致测定的4个年龄数据不一致,而经常存在t208<t206<t207<t206/207的顺序。引起不一致年龄的原因主要是不同子体的丢失程度不同,这时t206/207年龄最接近矿物结晶年龄。因为207Pb和206Pb化学性质相似,故丢失率也较一致,这一年龄值可消除因Pb丢失产生的误差。
为了排除由于矿物中子体同位素丢失引起的测年误差,U-Pb谐和曲线提供了较好的解决方法,U-Pb谐和曲线方程如下:
206Pb∗238U=(eλ238t−1)207Pb∗235U=(eλ235t−1) 在以206Pb*/238U为纵坐标和207Pb*/235U为横坐标的图中,对一个给定的年龄值,可得出相应的206Pb*/238U和207Pb*/235U值。通过选取不同的年龄,求出一条U-Pb体系的理论曲线,该曲线称之为U-Pb谐和曲线。
锆石年龄测定结果见表 1。共计测试样品点15个,除6、10、11三个测试点的U-Pb年龄谐和性较差外,绝大多数样品的谐和度大于90%。12个样品数据点位于谐和线上或其附近(图 4),均具有较好的谐和性,12个测试点中,12和14号测试点的206Pb/238U年龄分别为319.4±5.0Ma和318.3±2.8Ma,其可能为稍早期捕获锆石的年龄。其他锆石206Pb/238U年龄在297.0±5.4~306.2±3.5Ma之间,年龄加权平均值为299.2±2.2Ma(n=10,MSWD=2.7)。该年龄可代表灰白色中细粒含石榴子石二云母二长花岗岩单元锆石结晶年龄,即早二叠世早期侵位年龄。
表 1 含石榴子石二云母二长花岗岩单颗粒锆石U-Th-Pb激光烧蚀法测年数据(D4029-TW1)Table 1. LA-ICP-MS zircon U-Th-Pb data of garnet-bearing two-mica monzonitic granites测点编号 含量/10-6 Th/U 同位素比值 同位素年龄/Ma 206Pb 232Th 238U 207Pb/206Pb 1σ 206Pb/238U 1σ 207Pb/235U 1σ 208Pb/232Th 1σ 207Pb/206Pb 1σ 206Pb/238U 1σ 207Pb/235U 1σ 208Pb/232Th 1σ 谐和度 D4029- TW1-01 161 464 842 0.55 0.0524 0.0015 0.0475 0.0007 0.346 0.012 0.0152 0.0003 305.6 53.7 299.4 4.2 301.5 9.1 304.9 6.3 99% D4029- TW1-02 99 336 515 0.65 0.0528 0.0014 0.0475 0.0004 0.344 0.010 0.0152 0.0002 305.6 61.1 298.9 2.7 300.1 7.4 304.9 4.3 99% D4029- TW1-03 251 434 1323 0.33 0.0523 0.0011 0.0477 0.0006 0.347 0.010 0.0146 0.0003 283.4 37.0 300.4 3.6 302.5 7.4 292.0 5.1 99% D4029- TW1-04 395 562 1993 0.28 0.0494 0.0011 0.0486 0.0006 0.333 0.009 0.0147 0.0003 189.0 33.3 306.2 3.5 292.1 7.1 295.2 6.5 95% D4029- TW1-05 48 157 250 0.63 0.0588 0.0032 0.0485 0.0009 0.390 0.022 0.0149 0.0006 479.7 75.0 305.2 5.5 334.0 16.0 299.9 11.1 90% D4029- TW1-06 162 414 847 0.49 0.0651 0.0016 0.0473 0.0007 0.429 0.015 0.0171 0.0003 772.2 55.6 297.9 4.3 362.5 10.3 343.5 6.8 80% D4029- TW1-07 144 619 740 0.84 0.0549 0.0016 0.0472 0.0005 0.360 0.012 0.0151 0.0003 409.3 44.4 297.4 0 3.3 312.1 8.8 303.6 5.8 95% D4029- TW1-08 462 416 2460 0.17 0.0572 0.0014 0.0471 0.0009 0.374 0.013 0.0214 0.0004 498.2 25.9 297.0 5.4 322.9 9.3 427.9 8.5 91% D4029- TW1-09 210 406 1079 0.38 0.0531 0.0012 0.0472 0.0006 0.347 0.009 0.0156 0.0004 344.5 40.7 297.6 3.5 302.7 7.0 312.4 7.1 98% D4029- TW1-10 182 502 969 0.52 0.0632 0.0020 0.0465 0.0004 0.398 0.012 0.0168 0.0002 716.7 166.7 293.3 2.7 340.1 8.4 336.3 4.9 85% D4029- TW1-11 229 445 1128 0.39 0.0818 0.0022 0.0494 0.0007 0.545 0.014 0.0281 0.0010 1328. 7 58.3 310.6 4.4 441.8 9.2 560.4 19.3 65% D4029- TW1-12 166 334 763 0.44 0.0564 0.0017 0.0508 0.0008 0.395 0.013 0.0159 0.0003 450.1 52.8 319.4 5.0 338.1 9.7 318.8 6.4 94% D4029- TW1-13 403 448 1966 0.23 0.0533 0.0007 0.0469 0.0005 0.346 0.006 0.0158 0.0002 346.4 26.9 295.2 2.8 301.6 4.9 317.0 4.5 97% D4029- TW1-14 164 378 733 0.52 0.0539 0.0013 0.0506 0.0005 0.378 0.010 0.0161 0.0002 361.2 56.5 318.3 2.8 325.4 7.5 323.4 4.9 97% D4029- TW1-15 253 537 1168 0.46 0.0584 0.0015 0.0485 0.0006 0.395 0.013 0.0156 0.0003 538.9 51.9 305.4 3.9 338.3 9.3 313.4 5.2 90% 3.2 岩石地球化学测试
含石榴子石二云母二长花岗岩单元主量、微量元素测试样品共取7件,具体样品测试结果见表 2。SiO2含量为74.11% ~75.69%,Al2O3含量为13.59%~14.00%,CaO含量为0.37%~0.51%,Al2O3含量为13.52%~14%,Fe2O3为1.1%~1.73%,FeO为0.8% ~1.15%,MnO为0.073% ~1.12%,MgO为0.061% ~0.11%,P2O5为0.092% ~0.25%,K2O+Na2O为7.87% ~8.56%,K2O/Na2O为0.96~1.32,K2O含量略高于Na2O含量,具有富硅、铝,偏碱,而贫钙、镁、铁的特征。里特曼指数σ为1.93~2.35,A/CNK为1.112~1.196,在A/CNK- A/NK图解(图 5-a)中,样品点均位于S型强过铝质花岗岩区内,在SiO2-K2O图解(图 5-b)中,则全部为高钾钙碱性系列。
表 2 主量、微量和稀土元素含量及特征参数Table 2. Major, trace element and REE composition for garnet-bearing two-mica monzonitic granites含量 D4029-H1 D4142-H1 D4140-H1 D4139-H1 D4137-H1 D4136-H1 D4135-H1 SiO2 75.57 74.73 74.24 74.8 74.18 75.11 75.69 TiO2 0.051 0.033 0.042 0.044 0.053 0.061 0.039 Al2O3 13.59 14 13.78 13.8 13.69 13.52 13.65 Fe2O3 1.29 1.59 1.7 1.42 1.73 1.72 1.1 FeO 0.9 1.15 1.2 0.8 0.8 0.85 1 MnO 0.11 0.12 0.081 0.078 0.022 0.126 0.073 MgO 0.09 0.061 0.083 0.099 0.11 0.103 0.067 CaO 0.42 0.37 0.36 0.46 0.42 0.51 0.37 Na2O 4.11 4.18 3.9 4.08 3.68 3.66 3.76 K2O 4.34 4.06 4.35 3.94 4.88 4.21 4.21 P2O5 0.12 0.12 0.25 0.12 0.092 0.12 0.129 烧失量 0.22 0.5 1.01 0.93 0.94 0.69 0.9 总量 100.8 100.9 100.9 100.5 100.5 100.6 100.9 K2O+Na2O 8.45 8.24 8.25 8.02 8.56 7.87 7.97 K2O/Na2O 1.06 0.97 1.12 0.97 1.33 1.15 1.12 A/CNK 1.112 1.172 1.17 1.168 1.131 1.175 1.196 刚玉(C) 1.65 2.33 2.6 2.28 1.82 2.3 2.54 SI 0.84 0.55 0.74 0.96 0.99 0.98 0.66 AR 4.04 3.69 3.8 3.57 4.08 3.56 3.63 σ43 2.2 2.14 2.18 2.02 2.35 1.93 1.94 DI 93.62 92.52 92.65 92.72 93 91.93 93.18 Sc 4.58 5.72 4.18 4.43 3.05 3.58 3.96 Cr 17 15.5 20.1 23.9 23.6 16.2 12.6 Co 0.796 1.04 1.89 1.34 1.52 1.59 0.841 Ni 2.19 3.7 4.76 3.78 3.99 4.97 2.77 Rb 312 348 290 255 208 229 304 Nb 13 14.9 10.4 12.6 6.34 6.5 8.19 Ta 2.68 1.56 1.56 1.79 0.561 1.26 1.9 Th 6.63 8.42 7.73 7.78 6.29 5.58 6.07 Ba 35.6 19.6 60.4 62.9 63.3 87.1 35.4 Sr 19.1 10.9 25.9 35.1 36.5 36.5 14.8 V 2.75 3.9 5.52 7.05 8.36 17.6 7.23 Zr 44.5 53.1 62.1 61.9 51.6 56.8 43.4 Hf 2.52 3.58 3.16 3.1 2.5 2.66 2.63 K 36012.77 33689.36 36095.74 32693.62 40493.62 34934.04 34934.04 P 523.94 523.94 1091.55 523.94 401.69 523.94 563.24 Ti 306 198 252 264 318 366 234 Rb/Sr 16.335 31.927 11.197 7.265 5.699 6.274 20.541 Rb/Ba 8.764 17.755 4.801 4.054 3.286 2.629 8.588 Sr/Ba 0.537 0.556 0.429 0.558 0.577 0.419 0.418 Zr/Hf 17.659 14.832 19.652 19.968 20.64 21.353 16.502 Zr/Th 6.712 6.306 8.034 7.956 8.203 10.179 7.15 La 4.52 3.45 6.19 5.94 6.15 6.88 3.83 Ce 10.8 9.01 14.9 12.5 13.4 15.4 9.28 Pr 1.38 1.24 1.86 1.66 1.66 1.96 1.19 Nd 5.65 5.22 7.41 7.1 6.27 6.56 4.61 Sm 2.05 2.25 2.71 2.77 1.81 1.92 1.51 Eu 0.093 0.021 0.115 0.269 0.128 0.213 0.077 Gd 1.95 2.17 2.2 3.74 1.81 1.69 1.25 Tb 0.483 0.611 0.58 0.992 0.457 0.367 0.259 Dy 2.73 3.29 3.32 5.83 2.83 2.3 1.69 Ho 0.441 0.521 0.559 1.15 0.514 0.404 0.243 Er 1.2 1.38 1.74 3.71 1.6 1.21 0.651 Tm 0.217 0.248 0.336 0.67 0.308 0.252 0.137 Yb 1.43 1.59 2.5 4.81 2.34 1.68 1.08 Lu 0.206 0.211 0.367 0.747 0.345 0.251 0.142 Y 15.2 17.2 16 31.6 15 12.3 7.21 ΣREE 33.15 31.21 44.79 51.89 39.62 41.09 25.95 LREE 24.49 21.19 33.19 30.24 29.42 32.93 20.5 HREE 8.66 10.02 11.6 21.65 10.2 8.15 5.45 LREE/HREE 2.83 2.11 2.86 1.4 2.88 4.04 3.76 (La/Yb)N 3.16 2.17 2.48 1.23 2.63 4.1 3.55 δEu 0.14 0.03 0.14 0.26 0.21 0.35 0.17 δCe 1.03 1.05 1.05 0.94 0.99 1 1.04 注:主量元素含量单位为%,微量和稀土元素含量为10-6 稀土元素总量∑REE为25.95×10-6~51.89×10-6,含量偏低,平均值为38.24×10-6,远低于中国花岗岩平均值和天山-兴安造山系花岗岩平均值[41]。轻稀土元素(LREE)为20.50×10-6~33.19×10-6,平均值为27.42×10-6,重稀土元素(HREE)为5.45×10-6~ 21.65×10-6,平均值为10.82×10-6,LREE/HREE值为1.40~4.04,(La/Yb)N值为1.23~4.10,轻、重稀土元素基本无分馏,δEu=0.14~0.35,平均值为0.19,具强烈的Eu负异常。稀土元素配分曲线表现为Eu强烈负异常的平坦“雁式”型式(图 6-a)。
微量元素除Rb、Th平均含量较高外,其余元素均接近或低于世界花岗岩平均值而略高于天山-兴安造山系花岗岩平均值[41]。在微量元素原始地幔标准化蛛网图(图 6-b)上,强烈富集大离子亲石元素Rb、K和高场强元素Th,强烈亏损大离子亲石元素Ba、Sr和高场强元素Ti。
4. 讨论
4.1 岩石成因类型与岩浆源岩
Sylvester[14]对强过铝质S型花岗岩进行了系统的阐述,指出典型的S型花岗岩指含铝黑云母及其他含铝矿物,如白云母、堇青石、石榴子石等矿物的强过铝花岗岩类岩石,A/CNK > 1.1,刚玉标准分子大于1%[5]。查干哈达音亨嘎岩体花岗岩中出现数量不等的原生白云母、黑云母矿物,为(含)二云母花岗岩类,岩石副矿物中含大量石榴子石、独居石,未见包体组分;岩石中SiO2、Al2O3含量高,CaO、P2O5含量低,呈现富硅、铝而贫钙、磷特征,标准矿物计算中均出现刚玉,且含量均大于1%,A/CNK值均大于1.1,为过铝质高钾钙碱性系列岩石;微量元素中富集Rb、Th、U而强烈亏损Sr、Ba、Ti,表现出典型的强过铝质S型花岗岩特征。这与Zr -TiO2判别图解(图 7-a)和SiO2-P2O5判别图解(图 7-b)的结论一致。
尽管过铝质花岗岩类型和成因多种多样,Barbarin[8]将最主要的、最常见的过铝质花岗岩划分为2种类型:含白云母过铝质花岗岩类(MPG)和含堇青石过铝质花岗岩类(CPG),且代表了2种不同的成因。查干哈达音亨嘎岩体花岗岩颗粒粗,原生白云母矿物出现,且具清晰的自形形态等特征,副矿物中包含大量石榴子石矿物,还表现为典型的MPG类花岗岩。MPG花岗岩在造山带中往往呈巨大的深成侵入体或岩基形式产出,如喜马拉雅的马钠斯卢峰深成岩体、西欧海西造山带或苏格兰加里东造山带中的几个深成岩体群[44]。然而,查干哈达音亨嘎岩体是中亚造山带内部散布的少而孤立的深成岩体,明显不同于上述岩体,而具有与澳大利亚拉克伦褶皱带花岗岩一致的产出特征[6, 45]。因此,其可能产在地壳加厚的横切造山带的横推断层或逆掩型韧性剪切带中,由地壳加厚引起的深熔作用形成,且形成于地壳岩石“湿”的深熔作用和岩浆的结晶分离作用[44]。
CaO/Na2O值是判断源区成分一个极其重要的指标,在贫长石、富粘土的源区产生的过铝质花岗岩的熔融物中,该值较低(一般小于0.3)[44]。查干哈达音亨嘎岩体花岗岩CaO/Na2O值为0.09~0.14,反映其源区可能与变质泥岩部分熔融有关。在Rb/Sr-Rb/Ba图解[14] (图 8)上,样品也显示形成于变质泥岩的部分熔融。由实验得知,对于已知源区的成分和压力,较热的、大规模的部分熔融比相对冷的、小规模的部分熔融物Al2O3/TiO2值要低[44]。查干哈达音亨嘎岩体花岗岩Al2O3/TiO2值为221.64~424.24,反映其岩浆形成时温度可能较低,且源岩部分熔融程度也较低,这与锆石饱和温度揭示的岩浆形成温度729~757℃一致,也与锆石阴极发光图像特征一致。
图 8 S型花岗岩Rb/Sr-Rb/Ba图解[14]Figure 8. Diagram of Rb/Sr versus Rb/Ba of S-type granitioids4.2 构造背景环境分析
华北板块北缘和西伯利亚板块南缘缝合带上碰撞花岗岩的侵位年代应该是对碰撞缝合时间可靠的限制,碰撞花岗岩代表西伯利亚南缘和华北板块北缘两大板块最终碰撞缝合之后,由于陆壳加厚而发生重熔的产物,碰撞花岗岩的侵位时间代表缝合时间的上限,弧花岗岩的侵位时间代表缝合时间的下限[46]。最新研究表明,在中亚北造山带地区发育大量490~422Ma与俯冲相关的弧岩浆[47-50]和423~337Ma的碰撞花岗岩[20, 24, 31, 48, 51]。发育更多的是在碰撞造山之后具有后造山特征的碱性、过碱性花岗岩或具有造山后特征的A型花岗岩类[21, 23, 30, 52-54]。贺根山蛇绿岩中获取的微晶辉长岩和斜长花岗岩的年龄分别为354Ma和333Ma [50],在二连浩特地区获得的蛇绿岩年龄为354.2~344.8Ma [55]。说明在晚古生代,查干哈达音亨嘎岩体强过铝质花岗岩形成之前,二连—东乌旗一带已经发生了碰撞造山作用,之后为后碰撞-造山后的伸展扩张时期。区域上,在后碰撞-造山后时期,同期具有陆相沉积特征的宝力高庙组碎屑-火山岩不整合沉积于晚泥盆世安格尔音乌拉组海相地层之上,而在同碰撞造山期间,则缺失早石炭世沉积物,也较好地约束了这一结论。
据Sylvester [14]研究,强过铝质花岗岩可划分为高压型和高温型后碰撞型花岗岩类。其中,澳大利亚东南的拉克伦褶皱带中花岗岩属于高温型。本区强过铝质花岗岩表现为MPG类花岗岩,类似于澳大利亚东南的拉克伦褶皱带中的花岗岩,如发育规模较小,伴生大量的同构造-后构造钙碱性(I型)及SP型(S型)花岗岩侵位[21, 23-24, 27-28]。因此,查干哈达音亨嘎岩体强过铝质花岗岩还属于高温型强过铝质花岗岩类。该类花岗岩岩浆的形成往往缺乏阿尔卑斯和喜马拉雅山特有的高压变质作用和极端的地壳加厚作用及抬升作用,但具有板块汇聚作用及碰撞作用的特点,如早期岛弧岩石及薄皮逆冲断层[44],可与上述区域岩浆及构造特征对应。
5. 结论
(1) 通过岩石地球化学分析,该岩体具有富硅、铝,偏碱,而贫钙、镁、铁的特征,里特曼指数σ为1.93~2.35,A/CNK值为1.112~1.196,轻、重稀土元素基本无分馏,强烈的负Eu异常,强烈富集大离子亲石元素Rb、K和高场强元素Th,强烈亏损大离子亲石元素Ba、Sr和高场强元素Ti,为高钾钙碱性S型强过铝质花岗岩。
(2) 根据CaO/Na2O(0.09~0.14)值判断,查干哈达音亨嘎岩体花岗岩源区成分与变质泥岩部分熔融有关。
(3) 通过锆石U-Pb同位素定年,206Pb/238U年龄加权平均值为299.2±2.2Ma(n=10,MSWD=2.7);对比板块不同位置花岗岩的年龄限制,认为晚古生代,查干哈达音亨嘎岩体强过铝质花岗岩形成之前,二连—东乌旗一带已经发生了碰撞造山作用,之后为后碰撞-造山后的伸展扩张时期。
(4) 结合地球化学特征、同位素年龄及构造环境判别,认为区域内岩体为具有典型后碰撞特征的晚石炭世—早二叠世高温型强过铝质S型花岗岩类。
-
图 1 鄂尔多斯地块地形地貌、活动断裂、地震分布特征(a)和研究区地质构造简图[29](b)
F1—海原断裂带;F2—马东山褶皱逆冲带;F3—小关山断裂;F4—六盘山东麓断裂;F5—六盘山西麓断裂;F6—陇县-岐山-马召断裂;F7—固关-虢镇断裂。图中地震目录数据(MS≥3.0)来自USGS,震源机制解数据(MS≥4.5)来自GCMT
Figure 1. Landform, active tectonic characteristics and earthquake distribution of the Ordos block(a)and simplified geological map of the research region(b)
表 1 鄂尔多斯西南缘水压致裂地应力测量结果
Table 1 Results of in-situ stress measurement on the southwest margin of Ordos block
钻孔 深度/m 压裂参数/MPa 应力值/MPa 破裂方位 Pb Pr Ps PH P0 Sh SH Sv ZK7 128.50 5.77 4.09 2.91 1.29 0.64 2.91 3.99 3.41 157.50 10.44 4.70 3.36 1.58 0.93 3.36 4.44 4.17 181.20 7.11 5.37 3.96 1.81 1.16 3.96 5.35 4.80 N62°W 202.00 14.72 6.98 4.52 2.02 1.37 4.52 5.21 5.35 N56°W 213.50 13.97 10.74 6.56 2.14 1.49 6.56 7.44 5.66 231.50 17.70 9.53 6.25 2.32 1.67 6.25 7.54 6.13 N45°W 250.80 15.08 8.57 6.11 2.51 1.86 6.11 7.90 6.65 ZK10 206.00 6.77 4.63 4.29 2.06 1.09 4.29 7.15 5.46 248.00 9.91 5.36 4.10 2.48 1.51 4.10 5.43 6.57 N75°W 250.00 9.69 6.95 6.15 2.50 1.53 6.15 9.97 6.63 258.10 9.90 5.28 5.00 2.58 1.61 5.00 8.11 6.84 N65°W 269.50 11.94 5.67 5.25 2.70 1.72 5.25 8.35 7.14 N52°W 282.40 7.90 6.69 5.26 2.82 1.85 5.26 7.25 7.48 注:Pb为岩石原地破裂压力;Pr为破裂面重张压力;Ps为破裂面瞬时关闭压力;PH为静水压力;P0为孔隙压力;SH为最大水平主应力;Sh为最小水平主应力;Sv为根据上覆岩层重量计算的垂向主应力(岩石密度取2.65 g/cm3) 表 2 不同地区SH、Sh随深度变化情况
Table 2 Variation of SH and Sh with the depth in different regions
-
盛书中, 万永革, 黄骥超, 等.应用综合震源机制解法推断鄂尔多斯块体周缘现今地壳应力场的初步结果[J].地球物理学报, 2015, 58(2):436-452. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dqwlxb201502008 郭祥云, 蒋长胜, 王晓山, 等.鄂尔多斯块体周缘中小地震震源机制及应力场特征[J].大地测量与地球动力学, 2017, 37(7):675-685. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dkxbydz201707003 薛宏运, 鄢家全.鄂尔多斯地块周围的现代地壳应力场[J].地球物理学报, 1984, 27(2):144-152. doi: 10.3321/j.issn:0001-5733.1984.02.004 张培震.中国大陆岩石圈最新构造变动与地震灾害[J].第四纪研究, 1999, 19(5):404-413. doi: 10.3321/j.issn:1001-7410.1999.05.003 邓起东, 尤惠川.鄂尔多斯周缘断陷盆地带的构造活动特征及其形成机制, 现代地壳运动研究[M].北京:地震出版社, 1985:58-78. 徐锡伟, 程国良, 马杏垣, 等.华北及其邻区块体转动模式和动力来源[J].地球科学, 1994, 19(2):129-138. doi: 10.3321/j.issn:1000-2383.1994.02.005 陈强森, 鲍学伟, 徐树斌, 等.利用背景噪声反演鄂尔多斯块体及其南缘地区地壳速度结构[J].高校地质学报, 2013, 19(3):504-512. doi: 10.3969/j.issn.1006-7493.2013.03.012 张培震, 邓起东, 张竹琪, 等.中国大陆的活动断裂、地震灾害及其动力过程[J].中国科学:地球科学, 2013, 43(10):1607-1620. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=zgkx-cd201310005 张培震, 张会平, 郑文俊, 等.东亚大陆新生代构造演化[J].地震地质, 2014, 36(3):574-585. doi: 10.3969/j.issn.0253-4967.2014.03.003 Wang C Y, Sandvol E, Zhu L, et al.Lateral variation of crustal structure in the Ordos block and surrounding regions, North China, and its tectonic implications[J].Earth and Planetary Science Letters, 2014, 387:198-211. doi: 10.1016/j.epsl.2013.11.033
丰成君, 陈群策, 李国歧, 等.青藏高原东南缘丽江-剑川地区地应力测量与地震危险性[J].地质通报, 2014, 33(4):524-534. doi: 10.3969/j.issn.1671-2552.2014.04.009 李孟銮, 赵知军.宁夏地区现代构造应力场及其与地震活动的关系[J].地震研究, 1986, 9(3):33-48. http://www.cnki.com.cn/Article/CJFDTotal-DZYJ198603003.htm 单修政, 徐世芳, 段峰.鄂尔多斯周缘地带未来强震发生地区初探[J].华北地震科学, 2002, 20(1):10-15. doi: 10.3969/j.issn.1003-1375.2002.01.002 赵知军, 刘秀景, 康凌燕.宁夏及邻近地区震源机制解特征[J].西北地震学报, 2002, 24(2):162-166. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=xbdzxb200202011 范俊喜, 马瑾, 刁桂苓.由小震震源机制解得到的鄂尔多斯周边构造应力场[J].地震地质, 2003, 25(1):88-99. doi: 10.3969/j.issn.0253-4967.2003.01.009 张辉, 刘小凤, 赵凌云.由多个震源机制解分析甘肃及边邻地区应力场特征[J].华南地震, 2007, 27(2):33-39. doi: 10.3969/j.issn.1001-8662.2007.02.005 刘芳, 王晓山, 杨雅琼.内蒙古中西部地区小震震源机制解分析[J].大地测量与地球动力学, 2010, 30(S1):7-11. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dkxbydz2010z1002 Zoback M D, Healy J H.Friction, faulting and in-situ stress[J].Annales Geophsicae, 1984, 2:689-698.
Byerlee J.Friction of rocks[J].Pure and Applied Geophysics, 1978, 116(4/5):615-626. doi: 10.1007/BF00876528
王锋, 赵红格.鄂尔多斯地块西部构造研究中几个值得注意的问题[J].中国地质, 2007, 34(3):392-399. doi: 10.3969/j.issn.1000-3657.2007.03.004 邓起东, 程绍平, 闵伟, 等.鄂尔多斯块体新生代构造活动和动力学的讨论[J].地质力学学报, 1999, 5(3):13-21. doi: 10.3969/j.issn.1006-6616.1999.03.003 张岳桥, 廖昌珍, 施炜, 等.鄂尔多斯盆地周边地带新构造演化及其区域动力学背景[J].高校地质学报, 2006, 12(3):285-297. doi: 10.3969/j.issn.1006-7493.2006.03.001 Tapponnier P, Peltzer G, Le Dain A Y, et al.Propagating extrusion tectonics in asia:new insights from simple experiments with plasticine[J].Geology, 1982, 10(12):611. doi: 10.1130/0091-7613(1982)10<611:PETIAN>2.0.CO;2
Peltzer G, Tapponnier P, Zhang Z T, et al.Neogene and quaternary faulting in and along the Qinling Shan[J].Nature, 1985, 317(6037):500-505. doi: 10.1038/317500a0
Zhang Y Q, Vergely P, Mercier J L.Pliocene-Quaternary faulting pattern and left-slip propagation tectonics in North China[J].Episodes, 1999, 22(2):84-88. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dzlxxb199903007
Uyeda S, Kanamori H.Back-arc opening and the mode of subduction[J].Journal of Geophysical Research Solid Earth, 1979, 84(B3):1049-1061. doi: 10.1029/JB084iB03p01049
Northrup C J, Royden L H, Burchfiel B C.Motion of the Pacific plate relative to Eurasia and its potential relation to Cenozoic extension along the eastern margin of Eurasia[J].Geology, 1995, 23(8):719-722. https://ui.adsabs.harvard.edu/abs/1995Geo....23..719N/abstract
Tian X, Teng J, Zhang H, et al.Structure of crust and upper mantle beneath the ordos block and the yinshan mountains revealed by receiver function analysis[J].Physics of the Earth & Planetary Interiors, 2011, 184(3):186-193. https://www.sciencedirect.com/science/article/abs/pii/S0031920110002475
刘兴旺, 袁道阳, 吴赵, 等.六盘山断裂带活动性差异及其对六盘山隆升的影响[J].第四纪研究, 2016, 36(4):898-906. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dsjyj201604010 李新男.鄂尔多斯西南缘活动构造几何图像、运动特征及构造变形模式[J].国际地震动态, 2018, 473(5):46-48. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=gjdzdt201805009 徐锡伟, 于贵华, 陈桂华, 等.青藏高原北部大型走滑断裂带近地表地质变形带特征分析[J].地震地质, 2007, 29(2):201-207. doi: 10.3969/j.issn.0253-4967.2007.02.002 李清河, 郭守年, 吕德徽.鄂尔多斯西缘与西南缘深部结构与构造[M].北京:地震出版社, 1999:1-257. 孟庆筱, 景鹏旭, 何申海, 等.GPS约束下陇西地区断裂带现今滑动速率的非连续接触模拟研究[J].大地测量与地球动力学, 2018, 38(12):1227-1231. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dkxbydz201812003 王海燕, 高锐, 尹安, 等.深地震反射剖面揭示的海原断裂带深部几何形态与地壳形变[J].地球物理学报, 2012, 55(12):3902-3909. doi: 10.6038/j.issn.0001-5733.2012.12.003 孙赫, 徐晶, 柳皓元.基于InSAR的广义海原断裂带中东段现今深部运动特征[J].大地测量与地球动力学, 2017, 37(11):1141-1145. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dkxbydz201711009 许英才, 高原, 石玉涛, 等.鄂尔多斯块体西缘地壳介质各向异性:从银川地堑到海原断裂带[J].地球物理学报, 2019, 62(11):4239-4258. doi: 10.6038/cjg2019M0309 杜方, 闻学泽, 冯建刚, 等.六盘山断裂带的地震构造特征与强震危险背景[J].地球物理学报, 2018, 61(2):545-559. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dqwlxb201802014 邓起东, 张维岐.海原走滑断裂带及其尾端挤压构造[J].地震地质, 1989, 11(1):1-14. http://www.cnki.com.cn/Article/CJFDTotal-DZDZ198901000.htm Haimson B C, Cornet F H.ISRM suggested methods for rockstress estimation-Part 3:hydraulic fracturing(HF)and/orhydraulic testing of pre-existing fractures(HTPF)[J].International Journal of Rock Mechanics and Mining Sciences, 2003, 40(7/8):1011-1020.
Haimson B C, Rummel F.Hydrofracturing stress measurements in the Iceland research drilling project drill hole at Reydarfjordur, Iceland[J].Journal of Geophysical Research:Solid Earth, 1982, 87(B8):6631-6649. doi: 10.1029/JB087iB08p06631
Zoback M D, Apel R, Baumgärtner J, et al.Upper-crustal strength inferred from stress measurements to 6 km depth in the KTB borehole[J].Nature 1993, 365:633-635. doi: 10.1038/365633a0
Zhao X G, Wang J, Qin X H, et al.In-situ stress measurements and regional stress field assessment in the Xinjiang candidate area for China's HLW disposal[J].Engineering geology, 2015, 197:42-56. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=bc90eb6ae03f1f459cb05dbf5fabc79d
陈群策, 丰成君, 孟文, 等.5·12汶川地震后龙门山断裂带东北段现今地应力测量结果分析[J].地球物理学报, 2012, 55(12):3923-3932. doi: 10.6038/j.issn.0001-5733.2012.12.005 Aamodt L, Kuriyagawa M.Measurement of instantaneous shut in pressure in crystalline rock[C]//Monterey C A.Presented at the Workshop on Hydraulic Fracturing Stress Measurements.1981, 218(4): 715-716.
Hayashi K, Haimson B C.Characteristics of shut-in curves in hydraulic fracturing stress measurements and determination of in situ minimum compressive stress[J].Journal of Geophysical Research, 1991, 96(B11):18311-18321. doi: 10.1029/91JB01867
Hayashi K, Sakurai I.Interpretation of hydraulic fracturing shut-in curves for tectonic stress measurements[J].International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts, 1989, 26(6):477-482. doi: 10.1016/0148-9062(89)91424-1
丰成君, 陈群策, 吴满路, 等.水压致裂应力测量数据分析——对瞬时关闭压力Ps的常用判读方法讨论[J].岩土力学, 2012, 33(7):2149-2159. doi: 10.3969/j.issn.1000-7598.2012.07.035 Meng W, Chen Q C, Zhao Z, et al.Characteristics and implications of the stress state in the Longmen Shan fault zone, eastern margin of the Tibetan Plateau[J].Tectonophysics, 2015, 656:1-19. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=f45b5d5ec9bbdb7e940bf677268a190f
Tan C, Wang R, Sun Y, et al.Numerical modelling estimation of the'tectonic stress plane'(TSP)beneath topography with quasi-U-shaped valleys[J].International Journal of Rock Mechanics and Mining Sciences, 2004, 2(41):303-310. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=bd3311626cfabbaeea908a4ef99baf9c
黄禄渊, 杨树新, 崔效锋, 等.华北地区实测应力特征与断层稳定性分析[J].岩土力学, 2013, 34(S1):204-213. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=ytlx2013z1032 杨树新, 姚瑞, 崔效锋, 等.中国大陆与各活动地块、南北地震带实测应力特征分析[J].地球物理学报, 2012, 55(12):4207-4217. doi: 10.6038/j.issn.0001-5733.2012.12.032 王艳华, 崔效锋, 胡幸平, 等.基于原地应力测量数据的中国大陆地壳上部应力状态研究[J].地球物理学报, 2012, 55(9):3016-3027. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dqwlxb201209020 景锋, 盛谦, 张勇慧, 等.中国大陆浅层地壳实测地应力分布规律研究[J].岩石力学与工程学报, 2007, 26(10):2057-2062. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=yslxygcxb200710014 Anderson E M.The dynamics of faulting[J].Transactions of the Edinburgh Geological Society, 1905, 8(3):387-402. https://www.researchgate.net/publication/238338902_The_Dynamics_of_Faulting
Wang Q.Present-Day Crustal Deformation in China Constrained by Global Positioning System Measurements[J].Science, 2001, 294(5542):574-577. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=zgkx-ed200210001
徐婉桢, 孟国杰, 苏小宁.基于GPS观测的六盘山断裂震间闭锁特征研究[J].地震, 2016, 36(3):14-24. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=diz201603002 谢富仁, 陈群策, 崔效锋, 等.中国大陆地壳应力环境基础数据库[J].地球物理学进展, 2007, 22(1):131-136. doi: 10.3969/j.issn.1004-2903.2007.01.018 牛琳琳, 丰成君, 张鹏, 等.鄂尔多斯地块南缘地应力测量研究[J].地质力学学报, 2018, 24(1):25-34. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dzlxxb201801003 孟召平, 蓝强, 刘翠丽, 等.鄂尔多斯盆地东南缘地应力, 储层压力及其耦合关系[J].煤炭学报, 2013, (1):122-128. http://www.cnki.com.cn/Article/CJFDTotal-MTXB201301019.htm 陈小斌, 臧绍先, 刘永岗, 等.鄂尔多斯地块的现今水平运动状态及其与周缘地块的相互作用[J].中国科学院大学学报, 2005, 22(3):309-314. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=zgkxyyjsyxb200503008 Su S, Stephansson O.Effect of a fault on in situ stresses studied by the distinct element method[J].International Journal of Rock Mechanics and Mining Sciences, 1999, 8(36):1051-1056. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=c57c0d657e5062ed4a1ae2543a2c2197
Lin W, Yeh E C, Ito H, et al.Current stress state and principal stress rotations in the vicinity of the Chelungpu fault induced by the 1999 Chi-Chi, Taiwan, earthquake[J].Geophysical Research Letters, 2007, 34(16). http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=74831bd2876d3862fdff632d3e3a6d8f
Hudson J A, Harrison J P.Engineering rock mechanics:an introduction to the principles[M].Netherlands:Elsevier Science and Technology, 2000:1-456.
Hickman S, Zoback M.Stress orientations and magnitudes in the SAFOD pilot hole[J].Geophysical Research Letters, 2004, 31(15):1-4. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=10.1029/2004GL020043
许英才, 曾宪伟, 许文俊, 等.基于台阵的青藏高原东北缘海原-六盘山断裂带及邻区地壳结构研究[J].中国地震, 2018, 34(3):101-114. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=zgdz201803010 秦向辉, 陈群策, 谭成轩, 等.龙门山断裂带西南段现今地应力状态与地震危险性分[J].岩石力学与工程学报, 2013, 32(1):2870-2876. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=yslxygcxb2013z1038 张鹏, 秦向辉, 丰成君, 等.郯庐断裂带山东段深孔地应力测量及其现今活动性分析[J].岩土力学, 2013, 34(8):2329-2335. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=ytlx201308042 Zoback M D, Healy J H.In situ stress measurements to 3.5 km depth in the Cajon Pass Scientific Research Borehole:Implications for the mechanics of crustal faulting[J].Journal of Geophysical Research, 1992, 97(B4):5039-5057. doi: 10.1029/91JB02175
崔笃信, 郝明, 李煜航, 等.鄂尔多斯块体周缘地区现今地壳水平运动与应变[J].地球物理学报, 2016, 59(10):3646-3661. doi: 10.6038/cjg20161012 Liao C T.Stress change near the Kunlun fault before and after the Ms 8.1 Kunlun earthquake[J].Geophysical Research Letters, 2003, 30(20):2027. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=10.1029/2003GL018106