地球物理学报  2012, Vol. 55 Issue (8): 2611-2619   PDF    
新加坡地应力测量
孟文1,2 , 陈群策1,2 , 杜建军1,2 , 丰成君1,2 , 秦向辉1,2 , 安其美3     
1. 中国地质科学院地质力学研究所, 北京 100081;
2. 国土资源部新构造运动与地质灾害重点实验室, 北京 100081;
3. 中国地震局地壳应力研究所, 北京 100085
摘要: 2008年和2009年, 通过和新加坡宜康地质技术有限责任公司合作, 在新加坡18个浅钻孔中进行了水压致裂地应力测量. 通过测量获得了新加坡地壳浅部的现今应力状态, 为新加坡拟建地下管缆隧道设计施工提供了地应力依据, 也为世界应力环境数据库提供了参考. 测试结果表明, 在38.00~56.95 m深度域内最大水平主应力σh值为3.45~9.64 MPa, 最小水平主应力σh为2.49~6.28 MPa. 水平主应力和用上覆岩层重量计算的垂直应力σv之间的关系为σhhv,最大、最小水平主应力平均值分别是垂直应力平均值的5.13倍和3.23倍, 表明新加坡应力场以水平应力为主导. 初步分析认为, 新加坡陆域浅部构造应力场的量值大致呈南、北两端偏低, 中间高的态势, 在北端最大与最小水平主应力值都小于它们各自的平均值, 而南端则与平均值相当. 最大水平主应力优势方向为NE—NEE向, 该方向与前人在其邻区获得的构造应力场方向和由震源机制解反映的主压应力方向较为一致. 新加坡现今构造应力场状态与印度板块、澳大利亚板块、菲律宾板块与欧亚板块的汇聚作用有着密切的联系. 根据库仑摩擦滑动准则, 断层面摩擦系数取0.6~1.0, 分析了研究区内断层产生摩擦滑动的可能性.
关键词: 新加坡      水压致裂法      地应力测量      摩擦滑动     
In-situ stress measurement in Singapore
MENG Wen1,2, CHEN Qun-Ce1,2, DU Jian-Jun1,2, FENG Cheng-Jun1,2, QIN Xiang-Hui1,2, AN Qi-Mei3     
1. Institute of Geomechanics, Chinese Academy of Geological Sciences, Beijing 100081, China;
2. Key Laboratory of Neotectonic Movement and Geohazard, Ministry of Land and Resources, Beijing 100081, China;
3. Institute of Crustal Dynamics, China Earthquake Administration, Beijing 100085, China
Abstract: In the year 2008 and 2009, in-situ stresses were measured in 18 boreholes by the hydraulic fracturing method to determine the present-day stress state of the shallow crust in Singapore, which provide key parameters for the design and construction of underground tunnels and data for the study of world stress environment. The test results show that at depths ranging from 38.00~56.95 m, the maximum horizontal principal stress σh is 3.45~9.64 MPa, and the minimum horizontal principal stress σh is 2.49~6.28 MPa. The relative stress magnitudes are σhhv, and the mean values of maximum horizontal principal stress and minimum horizontal principal stress are respectively 5.13 and 3.23 times the vertical stress, which reveals that the horizontal stress is dominant. The magnitudes of the tectonic stress field in the shallow crust of Singapore exhibit heterogeneity, which are low in the north and south and high in the middle. The maximum and minimum principal stresses are both lower than their mean values in the northern end, and are well matched in the southern end. The dominant direction of the maximum principal tectonic stress is NE—NEE, in agreement with the orientation of the tectonic stress field and that from focal mechanism solutions in adjacent areas from previous work. The current stress state of Singapore has a close relationship with convergence of the India, Australia, Philippine and Eurasia plates. According to Coulomb friction criteria, and taking the friction coefficient as 0.6~1.0, the possible slide on the faults in Singapore has been analyzed in this work..
Key words: Singapore      Hydraulic fracturing      In-situ stress measurement      Frictional sliding     
1 引 言

新加坡位于连接印度洋与太平洋的马六甲海峡,北隔狭窄的柔佛海峡与马来西亚紧邻,南部隔新加坡海峡与印尼相望,是亚洲最重要的金融、服务和海运中心之一.有限的领土已满足不了经济的迅速发展和人口的不断增长,要充分利用有效的空间挖掘领土的潜力,需要实施各种地下工程.因此了解地质构造、地应力状态对地下构筑物围岩稳定性的影响是非常重要的.

地应力状态是评价地壳稳定性不可或缺的重要因素之一,也是地下隧道施工的重要参考因素[1].因而进行地应力测量,获取实测地应力资料是非常关键的.受新加坡宜康地质技术有限责任公司的委托,2008 年和2009 年,我们在新加坡进行了18个钻孔的地应力测量,这18个钻孔自北向南基本贯穿了新加坡国土.钻孔位置见图 1.

图 1 钻孔分布位置及最大主应力方向示意图 Fig. 1 Sketch showing location of drilling boreholes and the direction of the maximum horizontal principal stres

测量工作严格按照国际岩石力学学会颁布的规范进行,测量曲线标准、规范,曲线上各特征值确切.通过测量获得了新加坡地壳浅部主应力的量值范围和相互关系.初步分析了新加坡浅部构造应力场的分布特点,为新加坡地下岩土工程设计提供了翔实的地应力资料.

本次测量,自北向南贯穿新加坡全境,获得的测试结果不仅为其地下工程设计提供了重要依据,并且对世界应力图数据库具有补充价值,为研究板块作用机制、环太平洋地震带构造应力场提供了宝贵的资料.

2 区域地质构造概况

从区域地质构造的角度来看,新加坡属于晚古生代造山带,主要岩性为花岗岩和玄武岩1).新加坡邻近环太平洋火山地震带,印度尼西亚群岛、菲律宾群岛等均处于地震带上.位于其东南方向的加里曼丹岛位于印度-澳大利亚板块、太平洋板块和欧亚板块的斜向汇聚作用带内.新加坡虽未直接处于板块边界上,但印度-澳大利亚板块、太平洋板块和欧亚板块的汇聚作用和周围活动断裂的广泛分布[2-5]定会对其产生一定的影响.新加坡及其周围区域构造背景见图 2.

图 2 新加坡及周围区域构造背景简图(据文献[3]更改) Fig. 2 Simplified map showing tectonic setting of Singapore and adjacent areas (modified from Ref. [3])
3 水压致裂地应力测量方法原理

水压致裂地应力测量方法是国际岩石力学学会试验方法委员会建议的确定岩体应力方法之一[6].该方法利用一对可膨胀的橡胶封隔器在选定的测量深度封隔一段钻孔,然后通过泵入流体对封隔段加压直至在孔壁周围产生裂缝,压裂过程中同时记录压力随时间的变化曲线.由实测曲线可得破裂压力Pb,裂缝的瞬时闭合压力Ps,裂缝重新张开的压力Pr.再根据相关公式可计算得最大、最小水平主应力,垂直应力及岩石抗拉强度[7-9].如果诱发裂缝是竖直的即与垂直压力平行,则可根据印模器得到的裂缝的方向确定最大水平主应力的方向.

(1)

(2)

(3)

(4)

式中Po 为孔隙压力;γ 为测试段上覆岩体的容重;H为测试段上覆岩体的厚度.

由式(1)可见,孔隙压力是确定最大水平主应力的重要参数之一.国内外的研究和实际测量表明,在地壳浅部低渗透率岩石中,孔隙压力与水柱静压力大致相等,在计算最大水平主应力时可以采用水柱静压力代替孔隙压力[10-13].新加坡的地应力测量钻孔深度不足百米,在计算最大水平主应力值时,也采用水柱静压力代替孔隙水压力.

由于水压致裂地应力测量方法操作简单、不需要岩石力学参数参与计算,尤其是可直接确定最小主应力,测量深度理论上不受限制等[8].因而在水电、矿山、交通等工程领域及大陆动力学研究、区域地壳稳定性评价等研究领域得到了广泛应用[14-17].

4 测量结果及分析 4.1 测量结果

根据工程设计安排,18 个测量钻孔的深度一般都在60m 以内.自北向南钻孔编号分别为:PBH-01,PBH-66,PBH-80,PBH-12,PBH-18,PBH-138,PBH-24,PBH-90,PBH-151,PBH-30,PBH-37,PBH-42,PBH-46,PBH-132,PBH-104,PBH-54,PBH-130和PBH-119.

按照新加坡甲方要求,在每个钻孔中只进行一段压裂测量和印模定向.在测量过程中,始终按照国际岩石力学学会颁布的规范进行.测量曲线都比较标准.并根据国际岩石力学学会推荐的5种判读瞬时关闭压力方法(分别为:单切线法、dp/dt法、马斯卡特法、dT/dP法和压力-流量法)[6]中的dp/dt[7, 18-19]和dT/dP[7, 20-23]判读Ps.测量结果见表 1.

表 1 新加坡地应力测量结果 Table 1 Results of in-situ stress measurement in Singapore

前人很多研究表明,水压致裂应力测量结果可信度较高,基本可代表测量区域内的构造应力场状态[24-25].

4.2 测量结果分析讨论 4.2.1 水平主应力大小

图 3 是18 个钻孔水平主应力值的分布情况,由图可见,在测试范围内,水平主应力值呈非均匀分布,大致呈南、北两端偏低,中间偏高的态势.北端的最大水平主应力值一般为4~5 MPa,最小水平主应力值一般在2~3 MPa,均小于它们各自的平均应力值,而中间段除个别点外,最大水平主应力值一般为7~10 MPa,最小水平主应力值一般在4~6 MPa,南端应力值则与平均值相当.现场测量时观察到,各测量钻孔之间的地形、地貌变化不大,测区内局部断裂构造以及岩石风化程度的差异是导致应力场非均匀分布的主要原因.在靠近断裂带边缘、岩层风化较深以及裂隙较发育的钻孔中应力值相对较低.南北两端相比较,南端又高于北端,前人的分析认为,板块间的作用力沿整个苏门答腊-安达曼群岛破裂带由南向北递减[26],新加坡总体南高北低的应力状态可能是这一信息的反映.但限于我们没有确切的新加坡地质背景资料,因而仅作讨论.

图 3 各钻孔应力大小分布示意图 Fig. 3 Schematic illustration of magnitudes of stress in different boreholes

最大、最小水平主应力平均值分别是垂直应力平均值的5.13倍和3.23倍(表 1图 3).显示新加坡地层浅部具有较高水平应力.可以认为,较高水平应力及其较高的差应力是影响地下建筑物围岩稳定的主导因素.

新加坡地层浅部具有较高水平应力的原因比较复杂.分析认为,印度-澳大利亚板块、太平洋板块、菲律宾板块和欧亚板块的汇聚作用是其高应力背景的主因.由图 2可见,新加坡较为靠近板块汇聚作用带,板块间的挤压、俯冲以及地壳深部物质上涌,导致新加坡地壳整体性隆升,在地壳隆升区由于垂直应力不断减小导致水平应力不断加强.新加坡虽然较靠近印度尼西亚等地震多发区,但是几乎没有发生过地震,一方面是由于巨大的苏门答腊断裂带和加里曼丹岛上多条断裂消减了外围强烈构造活动对新加坡的影响程度,使其成为构造活动相对比较稳定、封闭较好的局部地块,另一方面表明新加坡地壳基底的完整性和刚性性质,在这种封闭较好的区块内有利于应变能的积累.

4.2.2 最大水平主应力方向

表 1中18个钻孔测量的最大水平主应力方向大多为NE40°-70°,相对比较离散.考虑原因可能为测量钻孔较浅,而地层浅部的应力方向对局部构造、岩石各向异性的反应比较敏感[27].但从表 1不难看出,新加坡现今构造应力场的主导方向为NE-NEE 向,平均为N47°E.

世界地应力测量中虽未有新加坡的测试结果,但其邻近区域(印度尼西亚)的应力测试结果具有一定的参考价值.由其他资料反映的新加坡邻区的主应力方向[28-30]及由震源机制解得到的主压应力轴方向[31-34]为NNE-NE 向,与新加坡的测试结果大体一致(如图 4).

图 4 新加坡及其邻区最大主应力方向示意图 (A)新加坡邻区(印度尼西亚)最大主应力方向图,截自世界应力图;(B)新加坡邻区震源机制解,改自文献[33] ; (C)印度尼西亚地震带的震源机制解最大主压应力方位玫瑰图,改自文献[34]. Fig. 4 Directions of the maximum principal stress in Singapore and adjacent areas (A) Directions of the maximum principal stress in Singapore and adjacent areas from world stress map; (B) Focal mechanism solutions in adjacent areas of Singapore (modified from Ref. [33]) ;C) Rose diagram of directions of the maximum principal stress in Indonesian seis¬ mic zone from focal mechanism solutions (modified from Ref. [34]).

分析认为,新加坡构造应力场是印度-澳大利亚板块、太平洋板块、菲律宾板块和欧亚板块的汇聚作用以及物质流上涌等多因素综合作用的结果.其主动力源来自印度板块以每年约50 mm 的速度向NNE 方向挤压[35],和澳大利亚板块以每年约60mm的速度向NE30°-35°方向俯冲[36].而菲律宾板块、太平洋板块对本区的影响相对较弱[37].

4.3 断层摩擦滑动的讨论

库仑准则指出,当断层面上的剪应力τ 大于或等于滑动摩擦阻力μσn 时,断层出现滑动.Zoback等[38-39]根据Byerlee的断层滑动条件[40],引入有效应力的概念,认为在可能产生滑动的断层面上最大与最小有效应力之比是摩擦系数μ 的函数,并用(5)式表示:

(5)

式中S1 和S3 分别为断层面上的最大与最小主应力值,P为孔隙压力,μ 为断层面摩擦系数

若最大与最小有效主应力之比小于((μ2 +1)1/2 +μ)2 值,则断层面稳定,若大于或等于此值则断层面在方位合适的面上可能发生活动.这里方位合适的面指断层面的法线方向与最大主应力S1的夹角为φ 的面.前人研究认为φμ 的关系为:φ =0.5(π/2+arctan-1μ)[13].Byerlee综合各种不同岩石类型的实验室摩擦试验资料,结果表明大部分岩石的摩擦系数μ 位于0.6~1.0之间.

根据测得的数据可知主应力之间的关系为:σH >σh >σv.利用式(5),其中S1 =σHS3 =σvP=Poμ 取0.6~1.0.得到断层滑动的临界条件分别为(σH -Po)/(σv-Po)=3.12(μ=0.6);(σH-Po)/(σv -Po)= 5.83 (μ = 1.0).即当(σH -Po)/(σv -Po)<3.12时,断层不会发生瞬间滑动;当3.12< (σH-Po)/(σv-Po)时,断层有滑动的可能.

图 5是18个钻孔的分析结果.各孔获得的(σH -Po)/(σv -Po)值均超过了使断层发生滑动的最小临界值,且多数达到并超过了使断层滑动的最大临界值.值得注意的是,此次测量获得的是地壳浅部的应力状态,一定程度上受地形和局部构造的影响[27, 41-43],而对断层的活动性影响更大的是深部应力状态.该计算结果还表明测区范围内应力值已积累到一定的程度.虽然隧道开挖深度较浅,花岗岩的物理力学指标较高,但是高水平应力以及高差应力也可能导致不良结构面、块体的变形破坏.

图 5 新加坡测区范围内断层活动性示意图 Fig. 5 Schematic illustration of fault activity of in the study area of Singapore
5 结 论

水压致裂地应力测量方法是获得岩体应力的重要方法之一,在科学研究和工程领域都得到了广泛的应用.本文基于新加坡地应力测量数据获得了该区地壳浅部应力状态,并分析讨论了新加坡的构造位置、该区板块运动特征及区域构造应力场的动力学来源,结论如下:

(1) 在新加坡18个测试钻孔中,38.00~56.95m深度域内测量的最大水平主应力值为3.45~9.64MPa,最小水平主应力值为2.49~6.28 MPa.应力值总体趋势为南北两端偏低,中间偏高.

(2) 测量的水平主应力与计算的垂直主应力之间的关系为σH >σh >σv,且最大、最小水平主应力平均值分别是垂直应力平均值的5.13 倍和3.23倍.显示新加坡地层浅部具有相对较高的水平应力,该应力状态有利于逆断层活动.

(3) 测量的最大水平主应力方向比较离散,为NE40°-70°,优势方向为NE-NEE,该方向与前人在其邻区获得的构造应力场方向和由震源机制解反映的主压应力方向较为一致.由前面分析讨论可知,新加坡现今构造应力场状态是印度-澳大利亚板块、太平洋板块、菲律宾板块和欧亚板块的汇聚作用以及物质流上涌,周围活动断裂广泛分布等多因素综合作用的结果.

(4) 根据库仑摩擦滑动准则,断层面摩擦系数取0.6~1.0,分析了区内断层产生摩擦滑动的可能性.各孔(σH -Po)/(σv-Po)比值均已超过了断层摩擦滑动的最小临界值,大多甚至超过了断层滑动的最大临界值,可见新加坡现今构造应力场具有较高的能量累积.较高的水平应力及其较高的差应力是影响地下建筑物围岩稳定的主导因素.

致谢

感谢新加坡宜康地质技术有限责任公司为测量工作提供的支持和帮助.感谢中国地质科学院地质力学研究所范桃园副研究员、孙东生博士为本文提供的建议.

参考文献
[1] Klee G, Rummel F, Williams A. Hydraulic fracturing stress measurements in Hong Kong. International Journal of Rock Mechanics and Mining Sciences , 1999, 36(6): 731-741. DOI:10.1016/S0148-9062(99)00036-4
[2] 赵财胜, 孙丰月, 李碧乐, 等. 马来西亚沙捞越Punda走滑构造及其动力学成因. 吉林大学学报 (地球科学版) , 2003, 33(4): 419–432. Zhao C S, Sun F Y, Li B L, et al. Punda strike-slip faults in Sarawak Malaysia and its dynamic process. Journal of Jilin University (Earth Science Edition) (in Chinese) , 2003, 33(4): 419-432.
[3] Hall R, Blundell D. Reconstructing Cenozoic SE Asia. Geological Society of London Special Publication , 1996, 106(1): 153-184. DOI:10.1144/GSL.SP.1996.106.01.11
[4] 朱章显, 杨振强. 东南亚新生代两类埃达克岩的分布、成因和含矿性. 地质力学学报 , 2008, 14(4): 328–338. Zhu Z X, Yang Z Q. Distribution, origin and mineralization of two types of Cenozoic adakite and adakite-like rock in Southeastern Asia. Journal of Geomechanics (in Chinese) , 2008, 14(4): 328-338.
[5] 郑勇, 傅容珊, 熊熊. 中国大陆及周边地区现代岩石圈演化动力学模拟. 地球物理学报 , 2006, 49(2): 415–427. Zheng Y, Fu R S, Xiong X. Dynamic simulation of lithospheric evolution from the modern China mainland and its surrounding areas. Chinese J. Geophys. (in Chinese) , 2006, 49(2): 415-427.
[6] Haimson B C, Cornet F H. ISRM suggested methods for rock stress estimation—Part 3: hydraulic fracturing (HF) and/or hydraulic testing of pre-existing fractures (HTPF). International Journal of Rock Mechanics & Mining Sciences , 2003, 40(7-8): 1011-1020.
[7] 丰成君, 陈群策, 吴满路, 等. 水压致裂应力测量数据分析——对瞬时关闭压力Ps的常用判读方法讨论. 岩土力学 , 2012, 33(7): 2149–2159. Feng C J, Chen Q C, Wu M L, et al. Analysis of the hydro-fracturing data——The discussion of the methods frequently used to determine the instantaneous shut-in pressure. Rock and Soil Mechanics (in Chinese) , 2012, 33(7): 2149-2159.
[8] Haimson B C. The hydrofracturing stress measuring method and recent field results. Int. J. Rock Mech. Min. Sci. & Geomech. Abstr. , 1978, 15(4): 167-178.
[9] Haimson B C. Measurement of in situ stress. Methods Exp. Phys. , 1987, 24(B): 377-408.
[10] Moos D, Zooback M D. Utilization of observations of well bore failure to constrain the orientation and magnitude of crustal stresses: Application to Continental, Deep Sea Drilling Project, and Ocean Drilling Program Boreholes. J. Geophys. Res. , 1990, 95(B6): 9305-9325. DOI:10.1029/JB095iB06p09305
[11] Haimson B C, Doe T W. State of stress, permeability, and fractures in the Precambrian granite of Northern Illinois. J. Geophys. Res. , 1983, 88(B9): 7355-7371. DOI:10.1029/JB088iB09p07355
[12] Zoback M D, Hickman S. In situ study of the physical mechanisms controlling induced seismicity at Monticello Reservoir, South Carolina. J. Geophys. Res. , 1982, 87(B8): 6959-6974. DOI:10.1029/JB087iB08p06959
[13] 苏恺之, 李方全, 张伯崇, 等. 长江三峡坝区地壳应力与孔隙压力综合研究. 北京: 地震出版社, 1996 . Su K Z, Li F Q, Zhang B C, et al. Comprehensive Study of Crustal Stress and Pore Pressure at the Site of Three Gorges Dam (in Chinese). Beijing: Seismological Press, 1996 .
[14] 安其美, 丁立丰, 王海忠, 等. 龙门山断裂带的性质与活动性研究. 大地测量与地球动力学 , 2004, 24(2): 115–119. An Q M, Ding L F, Wang H Z, et al. Research of property and activity of Longmen Mountain fault zone. Journal of Geodesy and Geodynamics (in Chinese) , 2004, 24(2): 115-119.
[15] 白世伟, 李光煜. 二滩水电站坝区岩体应力场研究. 岩石力学与工程学报 , 1982, 1(1): 45–55. Bai S W, Li G Y. Research on stress field around dam area of Ertan hydropower station. Chinese Journal of Rock Mechanics and Engineering (in Chinese) , 1982, 1(1): 45-55.
[16] 刘允芳. 在单钻孔中水压致裂法的三维地应力测量. 岩石力学与工程学报 , 1999, 18(2): 192–196. Liu Y F. 3-dimensional geostress measurement by hydraulic fracturing technique in one borehole. Chinese Journal of Rock Mechanics and Engineering (in Chinese) , 1999, 18(2): 192-196.
[17] 郭啟良, 王成虎, 马洪生, 等. 汶川Ms8.0级大震前后的水压致裂原地应力测量. 地球物理学报 , 2009, 52(5): 1395–1401. Guo Q L, Wang C H, Ma H S, et al. In-situ hydro-fracture stress measurement before and after the Wenchuan Ms8.0 earthquake of China. Chinese J. Geophys. (in Chinese) , 2009, 52(5): 1395-1401.
[18] Zoback M D, Haimson B C. Status of the hydraulic fracturing method for in situ stress measurements. Proc. U. S. Symp. Rock Mech. , 1982, 23: 143-156.
[19] Lee M Y, Haimson B C. Statistical evaluation of hydraulic fracturing stress measurement parameters. Int. J. Rock Mech. Min. Sci. & Geomech. Abstr. , 1989, 26(6): 447-456.
[20] 陈群策, 毛吉震, 张钧. 水压致裂应力测量数据分析——对关闭压力Ps值的讨论.// 水压致裂裂缝的形成和扩展研究. 北京: 地震出版社, 1999: 191-194. Chen Q C, Mao J Z, Zhang J. Analysis of the Hydro-Fracturing Data—the Discussion of the Techniques used to Determine the Shut-in Pressure.// The Research of the Generation and Extension of the Hydraulic Fracture (in Chinese). Beijing: Seismological Press, 1999: 191-194.
[21] Hayashi K, Sakurai I. Interpretation of hydraulic fracturing shut-in curves for tectonic stress measurements. Int. J. Rock Mech. Min. Sci. & Geomech. Abstr. , 1989, 26(6): 477-482.
[22] Cheung L S, Haimson B C. Laboratory study of hydraulic fracturing pressure data—how valid is their conventional interpretation?. Int. J. Rock Mech. Min. Sci. & Geomech. Abstr. , 1989, 26(6): 595-605.
[23] Hayashi K, Haimson B C. Characteristics of shut-in curves in hydraulic fracturing stress measurements and determination of in situ minimum compressive stress. J. Geophys. Res. , 1991, 96(18): 311-321.
[24] Haimson B C. 水压致裂法应力测量. //丁建民译. 地应力测量与研究.北京: 地震出版社, 1982: 93-122. Haimson B C. Hydrofracturing stress measurements. // Ding J M Trans. In-situ Stress Measurements and Research (in Chinese). Beijing: Seismological Press, 1982: 93-122.
[25] 翟青山, 毛吉震, 张钧, 等. 根据钻孔崩落资料确定剑川地区应力场方向. 地震地质 , 1989, 11(2): 46–52. Zhai Q S, Mao J Z, Zhang J, et al. Contemporary orientation of tectonic stress field in the Jianchuan determined from the borehole. Seismology and Geology (in Chinese) , 1989, 11(2): 46-52.
[26] Ammon C J, Ji C, Thio H K, et al. Rupture process of the 2004 Sumatra-Andaman Earthquake. Science , 2005, 308(5725): 1133-1139. DOI:10.1126/science.1112260
[27] 陈群策, 毛吉震, 侯砚和. 利用地应力实测数据讨论地形对地应力的影响. 岩石力学与工程学报 , 2004, 23(23): 3990–3995. Chen Q C, Mao J Z, Hou Y H. Study on influence of topography on in-situ stress by interpretation of measurement data of in-situ stress. Chinese Journal of Rock Mechanics and Engineering (in Chinese) , 2004, 23(23): 3990-3995.
[28] Zoback M L. First-and second-order patterns of stress in the lithosphere: The world stress map project. J. Geophys. Res. , 1992, 97(B8): 11703-11728. DOI:10.1029/92JB00132
[29] Mueller B, Reinecker J, Heidbach O, et al. The 2000 release of the World Stress Map (available online at www. world-stress-map. org).
[30] 黄玺瑛, 魏东平. 世界应力图2000版 (WSM2000) 介绍及使用说明. 地球物理学进展 , 2003, 18(2): 234–246. Huang X Y, Wei D P. WSM2000: The release 2000 world stress map. Progress in Geophysics (in Chinese) , 2003, 18(2): 234-246.
[31] McCloskey J, Lange D, Tilmann F, et al. The September 2009 Padang Earthquake. Nature Geoscience , 2010, 3(2): 70-71. DOI:10.1038/ngeo753
[32] 周云好, 许力生, 陈运泰. 2000年6月4日印度尼西亚苏门答腊南部MS8.0地震的震源机制. 地震学报 , 2002, 24(5): 462–469. Zhou Y H, Xu L S, Chen Y T. Mechanism of the June 4,2000 southern Sumatra, Indonesia, MS8.0 earthquake. Acta Seismologica Sinica (in Chinese) , 2002, 24(5): 462-469.
[33] 林纪曾, 梁国昭, 赵毅, 等. 东南沿海地区的震源机制与构造应力场. 地震学报 , 1980, 2(3): 245–257. Lin J Z, Liang G Z, Zhao Y, et al. Focal mechanism and tectonic stress field of coastal southeast China. Acta Seismologica Sinica (in Chinese) , 1980, 2(3): 245-257.
[34] 赵小艳, 苏有锦, 付虹, 等. 欧亚地震带现代构造应力场及其分区特征. 地震研究 , 2007, 30(2): 146–151. Zhao X Y, Su Y J, Fu H, et al. Recent tectonic stress field of Eurasian earthquake zone and the characteristics of its subareas. Journal of Seismological Research (in Chinese) , 2007, 30(2): 146-151.
[35] Molnar P, Tapponnier P. Cenozoic tectonics of Asia: effects of a continental collision. Science , 1975, 189(4201): 419-426. DOI:10.1126/science.189.4201.419
[36] 杨志根. 对2004年苏门答腊—安达曼岛地震的测地学机制讨论. 天文学进展 , 2007, 25(3): 270–278. Yang Z G. Discussions on possible geodetic mechanism for the 2004 Sumatra-Andaman islands earthquake. Progress in Astronomy (in Chinese) , 2007, 25(3): 270-278.
[37] 臧绍先, 宁杰远. 菲律宾海板块与欧亚板块的相互作用及其对东亚构造运动的影响. 地球物理学报 , 2002, 45(2): 188–197. Zang S X, Ning J Y. Interaction between Philippine sea plate (PH) and Eurasia (EU) plate and its influence on the movement eastern Asia. Chinese J. Geophys. (in Chinese) , 2002, 45(2): 188-197.
[38] Zoback M D, Healy J H. In-situ stress measurements to 3.5 km depth in Cajon. Pass Scientific Research Bore-hole: implications for the mechanics of crustal faulting. J. Geophys. Res. , 1992, 97(B4): 5069-5057.
[39] Zoback M D, Healy J H. Friction, faulting and in-situ stress. Annales Geophysics , 1984, 2(6): 689-698.
[40] Byerlee J O. Friction of rocks. Pageoph. , 1978, 116(4/5): 615-626.
[41] 朱焕春, 余启华, 赵海滨. 河谷地应力测值的数值检验. 岩石力学与工程学报 , 1997, 16(5): 471–477. Zhu H C, Yu Q H, Zhao H B. Numerical check of in-situ stresses measured at valley. Chinese Journal of Rock Mechanics and Engineering (in Chinese) , 1997, 16(5): 471-477.
[42] 朱焕春, 李浩. 论岩体构造应力. 水利学报 , 2001(9): 81–85. Zhu H C, Li H. On tectonic stress in rock mass. Journal of Hydraulic Engineering (in Chinese) , 2001(9): 81-85.
[43] 谢富仁, 崔效锋, 赵建涛, 等. 中国大陆及邻区现代构造应力场分区. 地球物理学报 , 2004, 47(4): 654–662. Xie F R, Cui X F, Zhao J T, et al. Regional division of the recent tectonic stress field in China and adjacent areas. Chinese J. Geophys. (in Chinese) , 2004, 47(4): 654-662.