The State of the In-situ Stress and Fault Slide Evaluation of Gonghe Basin, Qinghai Provice
-
摘要: 青海共和盆地地质构造条件复杂,断层十分发育,区域内既有地热、太阳能、矿产等资源丰富,存在龙羊峡水库诱发地震环境背景。因此,分析青海共和盆地地应力特征与地质结构易滑性对于青海东部地区防震减灾具有重要意义。对盆地及附近区域内19个钻孔、65条水压致裂数据和44条应力解除实测数据进行统计分析,并基于断层摩擦强度理论、Byerlee-Anderson理论等断层力学相关理论讨论了研究区域断层易滑性与地震危险性。研究结果表明:研究区域内应力状态在深度350 m左右由逆冲型转换为走滑型,与区域内分布北北西右旋高角度逆断层相吻合;区域内最大水平主应力优势方位为N45°~E60°;地应力场初步结果反演表明研究区域应力场以逆冲型为主,局部地区兼走滑特征,与震源机制解反演结果一致;区域内断层平均摩擦系数为0.41,断层处于稳定状态,即断层易滑性较低,侧压力系数与应力积累指标插值分析结果同样表明断层整体易滑性较低,局部浅部断层带应力积累水平较高,综合分析推测断层易沿NW-SE向滑动。Abstract: The Gonghe Basin has complex geological structural conditions, well-developed faults, abundant geothermal, solar, and mineral resources in the area, as well as the Longyangxia Reservoir-induced earthquake environment background. Therefore, analyzing the characteristics of the in-situ stress and the slidability of the geological structure in the region is very important for earthquake prevention and disaster reduction in the eastern part of Qinghai. Based on the analysis of 65 hydraulic fracturing data and 44 stress relief measured data from 19 boreholes in the basin and nearby areas, and based on the theory of fault friction strength, Byerlee-Anderson theory and other related fault mechanics theories, the fault slidability and seismic risk of the study area are discussed. The results show that the stress state in the study area changes from thrust to strike-slip stress state at a depth of about 350 m, which is consistent with the north-north-west right-lateral high-angle reverse faults in the area; The dominant position of the maximum horizontal principal stress in the area is N45°~E60°; The preliminary inversion of the in-situ stress field also shows that the stress field in the study area is dominated by thrusting, and the local area is also characterized by strike-slip, which is consistent with the inversion result of the focal mechanism solution; The average friction coefficient of the fault in the area is 0.41, and the fault is in a stable state, that is, the fault is relatively slippery; At the same time, the results of interpolation analysis of the lateral pressure coefficient and the stress accumulation index found that the overall slidability of faults in the Gonghe Basin is low, and the level of stress accumulation in local shallow fault zones is high. Comprehensive analysis predicts that it is easy to slide along the NW-SE direction.
-
Key words:
- Gonghe basin /
- in-situ stress /
- focal mechanism solution /
- stress state /
- slip tendency
-
引言
城市地球物理探测与传统的地球物理勘查差别很大,城市地质环境是被人类活动改造过的自然环境,实际情况会更加复杂。地下除地层、构造等天然地质体以外,还有地下管网系统、地下轨道交通、地下储存库等地下综合体。此外,城市地球物理场是动态演化的,在城市化前后甚至还会发生更加显著的变化(陈颙等,2003)。因此,城市地球物理探测具有探测深度浅、探测精度要求高、干扰因素多、场地受限等特点(李万伦等,2018)。目前国内外对目标层埋深仅有几十米(甚至十几米)的超浅层地震勘探经验不多,理论研究也不够深入。从地表至地下几十米的深度范围是各种人工地震勘探的盲区,但它对断层定位和活动性研究意义重大,近年来由于工程建设和防灾减灾的需要,已经开展了不少浅层地震勘探实践和技术研究(朱金芳等,2005)。在城市活断层探测中,用于第四纪覆盖区的活断层探测方法可分为地质学方法和地球物理方法,浅层或超浅层地震方法是地球物理方法调查城市活动断层研究中的重要手段之一(邓起东,2002;邓起东等,2003)。由于浅层地质条件以及新构造运动的特点,浅层或超浅层地震方法与常规地震勘探中采用的方法不完全相同。合理的数据采集方案是获得高分辨率地震剖面的重要基础;对浅层地震波传播速度的研究为数据处理中速度参考值的确定提供了理论依据。近年来,国内外浅层反射地震方法在城市活断层调查中的研究及其实例也相继发表,前人的研究成果反映了浅层地震方法的研究水平正在不断提高(Williams等,1995;潘纪顺等,2002;方盛明等,2002;刘保金等,2002;徐明才等,2005;段生全等,2005;杨晓平等,2007)。在利用浅层地震勘探方法确定隐伏活断层位置及其上断点埋深等工作中,一些地球物理工作者从地震波的激发、数据采集方法、室内资料处理与解释方法等方面进行了有益的探索,并取得了许多宝贵的经验,使得城市活断层浅层地震探测剖面质量和断层定位精度不断提高(何正勤等,2010;李大虎等,2010;杨歧焱等, 2015, 2016;许汉刚等,2016)。
青岛市地处胶东半岛西南部、黄海之滨,新构造运动以隆升为主,缺失第三系和第四系下部。第四纪晚期,受构造作用和海平面上升影响,接受十几米至二十余米厚以砂粘土为主的晚更新世—全新世河、海堆积物,直接覆盖于风化基岩(多为花岗岩)之上。东南部主要为花岗岩组成的低山丘陵,基岩在长期风化作用下形成了一定厚度的风化壳,并在山麓等和缓地带堆积有风化-坡积物,沿海港湾有第四纪晚期河海沉积;西北部主要是第四纪晚期的河海堆积平原或海湾,沉积有10—30m厚的第四系。因此,青岛市地震勘探属于超浅层地震勘探范畴。区域内发育若干北北东向断裂,活动性不甚强烈,地质地貌表现不很清楚,许多段落隐伏于第四系之下,断层的准确位置和活动性都不十分可靠(图 1)。与其他地区的有利探测环境相比,青岛市第四系很薄,横向物性变化大,垂向分层性差。对于不同的地震地质条件和探测目的,必须采用合理的野外采集方法和采集参数,因此,在一个新工区进行正式数据采集之前,应认真进行数据采集参数试验(何正勤等,2007)。本文通过在青岛王哥庄断裂上进行不同观测系统和采集参数的对比试验,获得了几米至二十几米深度范围内高分辨率的地震剖面,为隐伏活动断层位置及断层上断点的确定、钻孔剖面的布设和断层活动性的判定提供了可靠的科学依据。
1. 青岛地区地质构造概况
青岛及邻近地区在大地构造上属华北地台胶辽断隆,基底为元古代的变质岩系,经吕梁运动固结形成;古生代开始为地台发育时期,长期处于构造隆起状态,缺失古生代至三叠纪地层,致使结晶基底大面积出露。总体而言,青岛及邻近地区山势不高,以低山、丘陵和平原为主,海拔较高的山峰多为孤立分布,大致以沧口断裂为界划分为西北、东南2个不同的地貌分区。西北部地势较低,多为海湾、第四系冲洪积平原。由于新构造时期鲁东地区以隆升运动为主,青岛及邻近地区第四纪地层发育不良。许多地区基岩出露,部分地区沉积有晚更新世及全新世堆积物,局部地区发育有中更新世晚期的沉积物。市区东南部主要为花岗岩组成的低山丘陵,在山麓等和缓地带堆积有较厚的风化、坡积物,沿海港湾有第四纪晚期河海堆积;西北部是中生代地层组成的准平原和第四纪河海堆积平原(或海湾),沉积较厚的第四系,最大厚度近30米。
2. 超浅层地震勘探
2.1 青岛地区超浅层地震勘探的特点
(1)目的层比较浅,有效反射波大多在70ms以内,在反射波出现的时窗内,各种干扰波、特别是声波和面波很强,在接收反射波的后几道,出现折射波和反射波混合的情况,严重影响了有效反射波。
(2)在近震源处,震源干扰波条带内干扰严重,很难获取有效波,不同波阻抗界面形成的反射信号幅度远远小于干扰波振幅。
(3)青岛及附近地区沉积结构和地层厚度横向变化大,反射波速度和频率受此影响较大。
(4)青岛及附近地区覆盖层较薄,第四系底界面通常只有一个较强的风化界面,由于风化地层的风化程度是渐变的,风化基岩之间波阻抗差较小,第四系内部垂向分层差。基岩顶部常有数米厚的风化壳,有些基岩面甚至更浅,风化基岩与上覆地层之间的波阻抗差也不是很明显。
2.2 观测系统的设计和采集参数的选取
针对上述地震勘探特点,本文选取王哥庄断层(图 1中F4断裂)为实验对象,为了获得王哥庄断层从近地表至第四系底界面(深约30m)深度范围内的地震剖面,对同一深度范围按照探测精度由低到高的工作思路,在同一地震测线上采用不同的方法技术和探测参数,不同探测参数的地震剖面分别命名为TEST-a和TEST-b。地震数据采集使用德国DMT公司生产的SUMMIT数字地震仪,地震波激发采用了人工夯击和人工锤击,针对同一探测深度、不同精度的探测要求,采用的观测系统参数、地震波激发和数据采集参数见表 1。
表 1 观测系统参数和地震采集参数表Table 1. Survey parameters and seismic data acquisition details剖面编号 道间距/m 偏移距/m 接收道数 覆盖次数 激发震源 检波器固有频率/Hz 采样间隔/ms 记录长度/ms TEST-a 3 6 24 6 人工夯击 60 0.25 256 TEST-b 1 7 24 12 人工锤击 100 0.25 256 在不同道间距的超浅层地震记录上(图 2),图 2(a)、图 2(b)分别是TEST-a、TEST-b对应的超浅层地震共炮点记录。对比图 2(a)和图 2(b)可以看出,采用较大道间距、较小偏移距以及能量较大的震源可以获得一组很清楚的反射波组,但是近震源受干扰波影响严重,在超浅层不能获取有效反射波;采用较小道间距、较大偏移距以及能量较小震源,可以获取超浅层有效反射波,也可以减轻震源干扰波的影响,但是能量较弱,反射波信号相对变弱。
采用较小道间距和较弱震源,可以获取超浅层有效反射波,且垂向分层较好。由于采用小排列接收,道间距越小反射波能量衰减也越慢,因此受震源干扰波影响也变的越小,但是面波比较发育。由此可见,较小道间距和较小偏移距有利于保护反射波,延缓反射波的衰减,但是也带来了面波发育的问题;在保证探测深度要求的前提下,较弱的震源有利于避开震源干扰波的影响。
2.3 超浅层地震数据处理中的方法
采用Promax地震反射处理系统处理地震数据。在数据处理过程中,由于获得的原始记录的反射波近震源干扰严重,远震源反射波能量较强、干扰能量较弱等特点,因此选取记录信噪比高的远震反射波,可最大限度地提高地震资料的分辨率。其数据处理流程和方法主要包括:道编辑、静校正、真振幅恢复、地表一致性振幅处理、地表一致性反褶积、二维随机噪音衰减、f-k滤波、正常时差校正(NMO)、共中心点(CMP)叠加、叠后频率带通滤波和叠后剖面去噪等。
本文利用不同的滤波和处理方法来提高地震数据信噪比,图 3为同一记录数据在不同情况下进行f-k滤波获得的对比分析结果和频谱分析结果。其中,图 3(a)为原始数据和相对应的谱分析结果,可以看出在超浅层地震勘探中,地震反射波受浅层折射、面波和声波的影响,反射波信号相对较弱(图中红色圆圈),信噪比不高,通过频谱分析则显示地震数据主频有2个频段:50Hz和100Hz;图 3(b)为对原始数据应用f-k滤波后的结果及其频谱分析,由图可见应用f-k滤波后,数据中面波和声波得到有效减弱,反射波信号得到加强,信噪比也得到提高,频谱分析显示数据主频为100Hz,同时结果还显示面波频率范围与有效地震反射波接近,滤波困难,数据信噪比低;图 3(c)通过切除数据前几道信号,避免了近道面波对反射波的干扰,再进行f-k滤波,反射波信号明显增强,信噪比也得到进一步提高。对比结果表明,在浅层或超浅层地震勘探中,在取得信噪比较高的情况下,再对地震数据进行f-k滤波可以有效去除干扰波,提高地震反射波的信噪比。
2.4 超浅层地震剖面特征
图 4、图 5分别是TEST-a、TEST-b测线经过数据处理后的地震反射时间及地质解释剖面(图中方向左西右东)。在图 4地震反射时间剖面中,丰富的地层界面反射波组出现在剖面100ms以上的时间范围内,存在1组起伏变化的反射波组,由2—3个强相位组成;在图 5地震反射时间剖面中,丰富的地层界面反射波组出现在剖面60ms以上的时间范围内,存在2组起伏变化的反射波组,分别由2—3个强相位组成。根据剖面上反射波同相轴的横向变化特征和上下反射波组关系,在该测线剖面上解释了1个断点,在图 4和图 5中标识为F。断层F位于TEST-a测线CDP423#和TEST-b测线CDP87#处,反射波同相轴出现了错动,断层为倾向西、高倾角逆断层,根据测区地质资料,判断该断层为王哥庄断裂。
对比TEST-a、TEST-b测线的地震反射波叠加剖面发现,TEST-b测线在超浅层部分,特别是30ms以上部分有清晰的反射同相轴,说明小道距在超浅层地震勘探中,能揭示更多的断层结构细节。
3. 钻孔地质剖面及综合解
在王哥庄北TEST-a测线CDP423#和TEST-b测线CDP87#处的王哥庄断裂断点附近,布设了1组深19.1—23.5m、相邻孔距5—10m的钻孔(ZK1—ZK7),发现在ZK2与ZK6之间基岩顶面和第四系底部被断层错断4m左右,西升东降(图 6)。由图 6可见,除底部层位外,第四系各个层位均连续分布,未受断层构造变动,根据前期第四纪地层划分结果(包括测年结果),这套地层应为上更新统。因此,该处王哥庄断裂在晚更新世早期仍在活动。
4. 结论与讨论
(1)本文从城市活断层探测的角度出发,利用地质学、地球物理学等学科中的相关理论和方法,探讨超浅层地震在复杂地质构造背景下获得有效探测结果的前提条件。综合分析青岛地区地震勘探资料发现:①要依据扩展排列试验并参考观测系统参数选择原则,选好最小和最大偏移距,在超浅层探测中,近震源道记录中反射波与干扰波很难分辨,需要一定的偏移距才能将其分辨清楚;②最大偏移距不能过大,由于目标层埋深很浅,若排列长度过大,记录的多为折射波,所以需要采用短排列、小道间距进行数据采集。
(2)在对超浅层地震数据采集技术、资料处理与解释方法进行系统研究的基础上,以可清晰分辨出断层位置的典型剖面为例,证明在选择合适的观测系统和数据处理方法基础上,超浅层地震反射波法可以获取深度仅有十几米的地层反射信号,且大部分反射剖面都可较清楚地揭示出超浅部断层位置和断层特征。青岛王哥庄断裂地震勘探研究结果显示,超浅层地震反射方法在城市活断层探测中,特别是沿海覆盖层较浅的区域具有良好的应用前景。
(3)在对青岛王哥庄断裂超浅层地震反射时间剖面开展研究和解释的基础上,结合典型断点的钻孔联合剖面探测资料,研究了王哥庄断裂的准确位置、几何结构和活动性。研究表明,王哥庄断裂晚更新世以来可能也是活动的。
-
表 1 研究区域实测地应力数据
Table 1. Measured in-situ stress data in the study area
编号 测段号 测量
深度/m应力值/MPa 印痕方向 最大水平应力
系数KH,max最小水平应力
系数Kh,min侧应力
系数Kav应力积累
指标$ \mu _{\rm{m}}$最大水平
正应力σH最小水平
正应力σh垂直
应力σvGH_ZK1 1 49.46 7.70 5.99 4.49 NE37° 1.71 1.33 1.52 0.26 2 54.57 8.84 8.04 — NE10° 5.97 5.43 5.70 0.71 3 59.57 7.92 6.09 4.59 NW10° 1.73 1.33 1.53 0.27 4 69.81 5.81 4.00 3.00 — 1.94 1.33 1.64 0.32 5 78.57 7.02 4.70 3.90 — 1.80 1.21 1.51 0.29 6 81.33 8.97 5.30 4.50 NW54° 1.99 1.18 1.59 0.33 7 93.11 14.87 8.90 7.00 NW2° 2.12 1.27 1.70 0.36 8 96.05 14.34 8.50 5.50 NW1° 2.61 1.55 2.08 0.45 9 87.35 13.15 7.94 2.26 — 5.82 3.51 4.67 0.71 10 94.45 9.99 7.20 2.44 NW64° 4.09 2.95 3.52 0.61 11 84.95 12.75 7.30 — — 5.54 3.17 4.36 0.69 12 92.58 7.58 5.40 4.90 — 1.55 1.10 1.33 0.21 13 98.59 3.72 3.50 2.60 — 1.43 1.35 1.39 0.18 GH_ZK2 1 45.00 4.90 4.10 1.23 — 3.98 3.33 3.66 0.60 2 59.00 5.50 3.90 1.61 NE61° 3.41 2.42 2.92 0.55 3 85.00 6.50 5.20 2.32 NE21° 2.80 2.24 2.52 0.47 4 113.00 15.30 9.00 3.09 NE41° 4.95 2.91 3.93 0.66 5 183.00 32.90 16.80 5.00 NE53° 6.58 3.36 4.97 0.74 GH_ZK3 1 45.00 11.90 7.70 1.23 NE17° 9.67 6.26 7.97 0.81 2 165.00 22.30 11.70 4.51 NW6° 4.94 2.59 3.77 0.66 GH_ZK4 1 84.00 4.40 3.30 2.30 NE25° 1.92 1.44 1.68 0.31 2 162.00 6.20 4.10 4.43 NE30° 1.40 0.93 1.17 0.17 3 164.00 6.20 4.20 4.48 NE33° 1.38 0.94 1.16 0.16 GH_ZK5 1 147.00 7.50 5.00 4.02 NE36° 1.87 1.24 1.56 0.30 2 266.00 9.70 7.20 7.27 NE40° 1.33 0.99 1.16 0.14 3 284.00 8.40 6.30 7.76 NE42° 1.08 0.81 0.95 0.04 GH_ZK6 1 238.00 17.80 10.40 6.51 NE37° 2.74 1.60 2.17 0.46 2 335.00 17.20 10.40 9.16 NE32° 1.88 1.14 1.51 0.31 3 367.00 22.00 13.20 10.03 NE43° 2.19 1.32 1.76 0.37 GH_ZK7 1 41.00 5.50 3.90 1.12 NE61° 4.91 3.48 4.20 0.66 2 67.00 6.50 5.20 1.83 NE21° 3.55 2.84 3.20 0.56 3 94.00 15.30 9.00 2.57 NE41° 5.95 3.50 4.73 0.71 4 165.00 32.90 16.80 4.51 NE53° 7.29 3.72 5.51 0.76 GH_ZK8 1 45.00 11.10 7.70 1.23 NE17° 9.02 6.26 7.64 0.80 2 165.00 22.30 11.70 4.51 NW6° 4.94 2.59 3.77 0.66 GH_ZK9 1 212.00 10.10 8.10 5.80 NE54° 1.74 1.40 1.57 0.26 2 220.00 10.50 8.70 6.01 NE70° 1.75 1.45 1.60 0.27 3 235.00 11.10 8.80 6.42 NE48° 1.73 1.37 1.55 0.27 GH_ZK10 1 138.00 6.90 4.50 3.77 NW75° 1.83 1.19 1.51 0.29 2 140.00 7.50 5.70 3.83 NW86° 1.96 1.49 1.73 0.32 3 144.00 8.50 5.60 3.94 NE68° 2.16 1.42 1.79 0.37 GH_ZK11 1 186.00 7.30 6.30 5.09 NW64° 1.44 1.24 1.34 0.18 2 190.00 7.80 6.40 5.19 NW81° 1.50 1.23 1.37 0.20 3 198.00 9.00 5.40 5.41 NW70° 1.66 1.00 1.33 0.25 GH_ZK12 1 244.00 8.80 7.40 6.67 NW72° 1.32 1.11 1.22 0.14 2 252.00 9.60 7.40 6.89 NW76° 1.39 1.07 1.23 0.16 3 254.00 7.10 5.20 6.94 NE80° 1.02 0.75 0.89 0.01 GH_ZK13 1 147.00 7.50 4.50 4.02 NW0° 1.87 1.12 1.50 0.30 2 163.00 9.70 6.40 4.46 NW85° 2.18 1.44 1.81 0.37 3 179.00 11.60 6.60 4.89 NW80° 2.37 1.35 1.86 0.41 GH_ZK14 1 148.00 8.10 6.10 4.05 NW65° 2.00 1.51 1.76 0.33 2 150.00 11.00 7.10 4.10 NW82° 2.68 1.73 2.21 0.46 3 158.00 10.40 7.20 4.32 NW69° 2.41 1.67 2.04 0.41 GH_ZK15 1 306.00 12.60 7.30 8.37 NW70° 1.51 0.87 1.19 0.20 2 312.00 13.20 7.90 8.53 NW83° 1.55 0.93 1.24 0.21 3 322.00 12.60 8.20 8.80 NE83° 1.43 0.93 1.18 0.18 GH_ZK16 1 290.00 9.10 5.50 7.93 NW64° 1.15 0.69 0.92 0.07 2 294.00 8.40 5.60 8.04 NW82° 1.05 0.70 0.88 0.02 3 300.00 10.60 6.60 8.20 NW73° 1.29 0.80 1.05 0.13 GH_ZK17 1 250.00 11.00 6.50 6.84 NW74° 1.61 0.95 1.28 0.23 2 252.00 11.80 6.80 6.89 NW85° 1.71 0.99 1.35 0.26 3 258.00 12.00 6.90 7.05 NE80° 1.70 0.98 1.34 0.26 GH_ZK18 1 480.00 23.60 13.00 13.12 NE44° 1.80 0.99 1.40 0.29 2 492.00 24.80 13.40 13.45 NE42° 1.84 1.00 1.42 0.30 GH_ZK19 1 264.00 11.90 6.60 7.22 NE15° 1.65 0.91 1.28 0.24 2 267.00 13.10 7.10 7.30 — 1.79 0.97 1.38 0.28 表 2 300 m深度处测点K、μ值
Table 2. K value and friction coefficient at a depth of 300 m
测点 KH,max Kh,min 有效正应力σ/MPa 剪应力
τ/MPaμ GH_ZK 1 2.80 1.22 6.21 3.38 0.54 GH_ZK 2 4.05 1.94 10.67 5.32 0.50 GH_ZK 3 1.72 0.67 4.02 1.84 0.46 GH_ZK 4 1.82 0.86 4.31 1.80 0.42 GH_ZK 5 1.41 0.98 4.18 0.89 0.21 GH_ZK 6 2.41 1.31 6.41 2.54 0.40 GH_ZK 7 4.00 1.98 10.81 5.23 0.48 GH_ZK 8 1.90 0.67 4.14 2.16 0.52 GH_ZK 9 2.01 1.29 6.03 1.81 0.30 GH_ZK 10 2.13 1.08 5.11 2.21 0.43 GH_ZK 11 1.83 1.05 4.72 1.65 0.35 GH_ZK 12 1.67 1.00 4.39 1.39 0.32 GH_ZK 13 2.27 1.10 5.28 2.47 0.47 GH_ZK 14 2.43 1.32 6.50 2.56 0.39 GH_ZK 15 1.88 1.01 4.57 1.81 0.40 GH_ZK 16 1.64 0.86 4.20 1.43 0.34 GH_ZK 17 1.98 1.00 4.60 2.02 0.44 GH_ZK 18 2.22 1.15 5.48 2.30 0.42 GH_ZK 19 2.13 1.02 4.77 2.29 0.48 -
[1] 安其美, 丁立丰, 王海忠等, 2004. 龙门山断裂带的性质与活动性研究. 大地测量与地球动力学, 24(2): 115—119.An Q. M., Ding L. F., Wang H. Z., et al., 2004. Research of property and activity of Longmen mountain fault zone. Journal of Geodesy and Geodynamics, 24(2): 115—119. (in Chinese) [2] 陈群策, 安其美, 孙东生等, 2010. 山西盆地现今地应力状态与地震危险性分析. 地球学报, 31(4): 541—548.Chen Q. C., An Q. M., Sun D. S., et al., 2010. Current in-situ stress state of Shanxi basin and analysis of earthquake risk. Acta Geoscientia Sinica, 31(4): 541—548. (in Chinese) [3] 丁仨平, 2008. 西秦岭—祁连造山带(东段)交接部位早古生代构造格架及构造演化. 西安: 长安大学.Ding S. P., 2008. Early Palaeozoic tectonic framework and evolution in the junction of Western Qinling orogenic belt and Qilian Orogenic belt. Xi’an: Chang’an University. (in Chinese) [4] 董治平, 姚政生, 雷芳, 1992. 青海东部活断层与现代构造应力场. 地壳形变与地震, 12(4): 64—71.Dong Z. P., Yao Z. S., Lei F., 1992. Active faults and modern tectonic stress field in the region of eastern Qinghai. Crustal Deformation and Earthquake, 12(4): 64—71. (in Chinese) [5] 都昌庭, 2001. 共和地震震源机制解特征. 高原地震, 13(4): 1—5.Du C. T., 2001. Characteristics of the mechanism solutions of Gonghe earthquakes. Plateau Earthquake Research, 13(4): 1—5. (in Chinese) [6] 黄禄渊, 杨树新, 崔效锋等, 2013. 华北地区实测应力特征与断层稳定性分析. 岩土力学, 34(S1): 204—213.Huang L. Y., Yang S. X., Cui X. F., et al., 2013. Analysis of characteristics of measured stress and stability of faults in North China. Rock and Soil Mechanics, 34(S1): 204—213. (in Chinese) [7] 李兵, 丁立丰, 王建新等, 2019. 山东蓬莱近海岸的地应力状态及断层稳定性评价. 地质力学学报, 25(4): 459—466. doi: 10.12090/j.issn.1006-6616.2019.25.04.043Li B., Ding L. F., Wang J. X., et al., 2019. The state of the in-situ stress and fault stability evaluation of the Penglai coast. Journal of Geomechanics, 25(4): 459—466. (in Chinese) doi: 10.12090/j.issn.1006-6616.2019.25.04.043 [8] 刘卓岩, 王成虎, 徐鑫等, 2017. 基于地应力实测数据分析郯庐断裂带中段滑动趋势. 现代地质, 31(4): 869—876. doi: 10.3969/j.issn.1000-8527.2017.04.021Liu Z. Y., Wang C. H., Xu X., et al., 2017. Slip tendency analysis of the Mid-segment of Tan-Lu fault belt based on stress measurements. Geoscience, 31(4): 869—876. (in Chinese) doi: 10.3969/j.issn.1000-8527.2017.04.021 [9] 任海东, 王涛, 2017. 东昆仑—西秦岭造山带对接处三叠纪花岗质岩石时空演化、物源特征对比及其大地构造意义. 地球学报, 38(S1): 59—63. doi: 10.3975/cagsb.2017.s1.16Ren H. D., Wang T., 2017. Temporal-spatial variations, sources and tectonic significances of the Triassic granitic rocks in the junction part of the east Kunlun and West Qinling orogen, Central China. Acta Geoscientia Sinica, 38(S1): 59—63. (in Chinese) doi: 10.3975/cagsb.2017.s1.16 [10] 王成虎, 丁立丰, 李方全等, 2012. 川西北跨度23 a的原地应力实测数据特征及其地壳动力学意义分析. 岩石力学与工程学报, 31(11): 2171—2181. doi: 10.3969/j.issn.1000-6915.2012.11.004Wang C. H., Ding L. F., Li F. Q., et al., 2012. Characteristics of in-situ stress measurement in northwest Sichuan basin with timespan of 23 years and its crustal dynamics significance. Chinese Journal of Rock Mechanics and Engineering, 31(11): 2171—2181. (in Chinese) doi: 10.3969/j.issn.1000-6915.2012.11.004 [11] 王成虎, 宋成科, 郭启良等, 2014. 利用原地应力实测资料分析芦山地震震前浅部地壳应力积累. 地球物理学报, 57(1): 102—114. doi: 10.6038/cjg20140110Wang C. H., Song C. K., Guo Q. L., et al., 2014. Stress build-up in the shallow crust before the Lushan Earthquake based on the in-situ stress measurements. Chinese Journal of Geophysics, 57(1): 102—114. (in Chinese) doi: 10.6038/cjg20140110 [12] 谢富仁, 陈群策, 崔效锋等, 2007. 中国大陆地壳应力环境基础数据库. 地球物理学进展, 22(1): 131—136. doi: 10.3969/j.issn.1004-2903.2007.01.018Xie F. R., Chen Q. C., Cui X. F., et al., 2007. Fundamental database of crustal stress environment in continental China. Progress in Geophysics, 22(1): 131—136. (in Chinese) doi: 10.3969/j.issn.1004-2903.2007.01.018 [13] 徐纪人, 赵志新, 2006. 青藏高原及其周围地区区域应力场与构造运动特征. 中国地质, 33(2): 275—285. doi: 10.3969/j.issn.1000-3657.2006.02.005Xu J. R., Zhao Z. X., 2006. Characteristics of the regional stress field and tectonic movement on the Qinghai-Tibet Plateau and in its surrounding areas. Geology in China, 33(2): 275—285. (in Chinese) doi: 10.3969/j.issn.1000-3657.2006.02.005 [14] 许忠淮, 汪素云, 黄雨蕊等, 1987. 由多个小震推断的青甘和川滇地区地壳应力场的方向特征. 地球物理学报, 30(5): 476—486.Xu Z. H., Wang S. Y., Huang Y. R., et al., 1978. Directions of mean stress axes in southwestern China deduced from microearthquake data. Chinese Journal of Geophysics, 30(5): 476—486. (in Chinese) [15] 张盛生, 张磊, 田成成等, 2019. 青海共和盆地干热岩赋存地质特征及开发潜力. 地质力学学报, 25(4): 501—508. doi: 10.12090/j.issn.1006-6616.2019.25.04.048Zhang S. S., Zhang L., Tian C. C., et al., 2019. Occurrence geological characteristics and development potential of hot dry rocks in Qinghai Gonghe basin. Journal of Geomechanics, 25(4): 501—508. (in Chinese) doi: 10.12090/j.issn.1006-6616.2019.25.04.048 [16] Anderson E. M., 1905. The dynamics of faulting. Transactions of the Edinburgh Geological Society, 8(3): 387—402. [17] Bieniawski Z. T., 1984. Rock mechanics design in mining and tunneling. Rotterdam: A. A. Balkemn. [18] Brown E. T., Hoek E., 1978. Trends in relationships between measured in-situ stresses and depth. International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts, 15(4): 211—215. [19] Byerlee J., 1978. Friction of rocks. Pure and Applied Geophysics, 116(4): 615—626. [20] Herget G., 1987. Stress assumptions for underground excavations in the Canadian shield. International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts, 24(1): 95—97. [21] Jamison D. B., Cook N. G. W., 1980. Note on measured values for the state of stress in the Earth's crust. Journal of Geophysical Research: Solid Earth, 85(B4): 1833—1838. doi: 10.1029/JB085iB04p01833 [22] Lee J. B., Chang C. D., 2009. Slip tendency of Quaternary faults in southeast Korea under current state of stress. Geosciences Journal, 13(4): 353—361. doi: 10.1007/s12303-009-0033-1 [23] Paterson M. S., Wong T. F., 1978. Experimental rock deformation: the brittle field. New York: Springer. [24] Rummel F., 1986. Stresses and tectonics of the upper continental crust-a review. In: Proceedings of the ISRM International Symposium. Stockholm, Sweden: ISRM, 86—177. [25] Savage W. Z., Swolfs H. S., Amadei B., 1992. On the state of stress in the near-surface of the earth's crust. Pure and Applied Geophysics, 138(2): 207—228. doi: 10.1007/BF00878896 [26] Sheorey P. R., 1994. A theory for In Situ stresses in isotropic and transverseley isotropic rock. International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts, 31(1): 23—34. [27] Townend J., Zoback M. D., 2000. How faulting keeps the crust strong. Geology, 28(5): 399—402. doi: 10.1130/0091-7613(2000)28<399:HFKTCS>2.0.CO;2 [28] van Heerden W. L., 1976. Practical application of the CSIR triaxial strain cell for rock stress measurements. In: Symposium I. S. R. M., ed., Investigation of Stress in Rock: Advances in Stress Measurement. Barton, ACT: Institution of Engineers, Australia. [29] Vavryčuk V., 2014. Iterative joint inversion for stress and fault orientations from focal mechanisms. Geophysical Journal International, 199(1): 69—77. doi: 10.1093/gji/ggu224 [30] Zoback M. D., Townend J., Grollimund B., 2002. Steady-state failure equilibrium and deformation of intraplate lithosphere. International Geology Review, 44(5): 383—401. doi: 10.2747/0020-6814.44.5.383 -