Contemporary tectonic stress field on boundaries of active tectonic blocks in and around western China
-
摘要: 利用哈佛全球矩心矩张量解数据和许忠淮认为1920mdash;1999年可靠的中国大陆震源机制解数据, 反演了中国西部及邻区活动地块边界带上现代构造应力场.通过对FMSI反演程序多次的输入和检验, 得到了边界带上的应力场.边界带上最大主压应力sigma;1轴绝大多数近水平. 在90deg;E以西的中国西部大陆及邻区, sigma;1轴水平方向基本上为近SN向;在青藏高原的东北部, sigma;1轴水平方向基本上为近NE向;在青藏高原的东南部, sigma;1轴水平方向绕喜马拉雅构造东端顺时针方向旋转.最小主压应力sigma;3轴倾角呈两极分布,西域地块区内活动地块边界带和青藏地块区内东北缘部分段sigma;3轴倾角较陡, 而青藏地块区内sigma;3轴倾角近水平, 所以西域地块区和青藏地块区内东北部相对于其它大部分青藏地块区, 有更多的逆冲地震.应力场在同一个边界带具有非均匀性. 北天山带、南天山带、西秦岭mdash;德令哈带、岷山mdash;龙门山带和安宁河mdash;小江带的非均匀性相对要小一些, 西昆仑带、海原mdash;祁连带、东昆仑带、玛尼mdash;玉树带、澜沧江带和滇西西边界带的非均匀性相对要大, 而喀喇昆仑mdash;嘉黎带和喜马拉雅带的非均匀性最显著.由于震源机制解数据的限制, 本文给出的是边界带上部分段的应力场.
-
关键词:
- 活动地块边界带 构造应力场 震源机制解
Abstract: Contemporary tectonic stress field on boundaries of active tectonic blocks in and around western China mainland is obtained from the data of Harvard CMT solution and reliable focal mechanism solutions from 1920 to 1999 selected by Xu Zhonghuai. We obtained the stress field by testing input and checking FMSI code(Gephart) many times. Most of the maximum principal stress sigma;1 axes on active tectonic boundaries are horizontal. To the west of 90deg;E in and around western China mainland, the horizontal sigma;1 axes are nearly in N-S direction. In the northeastern area of Qingzang plateau, the horizontal sigma;1 axes are nearly NE-SW. In the southeastern area of Qingzang plateau, the horizontal sigma;1 axes turn around clockwise around eastern end of Himalaya arc. The minimum principal stress sigma;3 axes show two dominant dip angles. In Xiyu active tectonic block region and on some segments of the northeastern edge of Qingzang plateau, the sigma;3 axes are steep, while in the Qingzang block region the sigma;3 axes are horizontal. So there are more thrust earthquakes in the Xiyu block region and on some segments of the northeastern edge of Qingzang plateau than in the Qingzang plateau block region. The stress field on the same active tectonic boundary is heterogeneous. On north Tianshan belt, south Tianshan belt, west Qinling-Delingha belt, Minshan-Longmenshan belt and Anninghe-Xiaojiang belt, the heterogeneity is relatively weak. On west Kunlun belt, Haiyuan-Qilian belt, east Kunlun belt, Mani-Yushu belt, Lancangjiang belt and western belt of west Yunnan block, the heterogeneity is relatively strong. The heterogeneity on Kalakunlun-Jiali belt and Himalaya belt is more notable. Due to limited focal mechanism data, this paper only gave the stress field on some segments of the active tectonic block boundaries. -
引言
地电阻率前兆观测已在我国开展了40多年,观测结果表明大震前在震源区及其附近区域通常会出现视电阻率的变化,因此该方法成为普遍认同的地震前兆观测手段(张国民等,2001;杜学彬,2010).然而,我国台站的地电阻率观测大部分采用对称四极装置在地表进行观测,得到的视电阻率观测数据是观测范围内不同深度处各地层电阻率的综合反映.因受降雨、地下潜水位变化和农业灌溉等因素的影响,浅部地层(地下数米至数十米)的电阻率会发生周期性的变化,这种变化显然与地震孕育无关.而赵和云和钱家栋(1988)的研究表明,当震源强度与震中距相同时,不同力学结构下介质中的体应变差别明显,这将导致浅部松软覆盖层的电阻率在地震孕育过程中变化最小,而弹性模量较高的深部介质的电阻率变化较大,从而使深部电阻率的变化具有了地震前兆意义.
为获得深部电阻率的变化,前人根据直流电法中的电测深理论发展了地震地电阻率多极距观测方法(钱家栋等,1985;钱家栋,赵和云,1988;薛顺章等,1994),其目的是通过观测不同极距的地电阻率变化,反演地下介质中不同深度的真电阻率的变化过程,为研究该变化与地震或非地震物理关联提供依据.因该方法具有较为可靠的物理基础,美国、前苏联和日本均已先后开展了地电阻率多极距观测的实验研究(Mazzella,Morrison,1974),我国亦开展了大量的理论与实验研究(赵和云,钱家栋,1987;钱家栋,赵和云,1988;薛顺章等,1994;冯锐等,2001;汪雪泉等,2002;李艳东等,2004;Lu et al,2004;毛先进等,2008;王兰炜等,2011).
以往的反演研究均视地下介质的电阻率为一维分布,但由于地壳介质导电性的空间不均匀性,台址下的电阻率结构可能与一维分布存在不同程度的差别,因此往往需要考虑二维结构的情况.本文先研究了一维结构的反演结果,然后考虑更符合实际情况的二维结构,即对二维结构的多极距观测结果进行正演,模拟出多极距观测视电阻率的时间序列,然后对模拟结果作一维反演,再将反演结果与实际电阻率的变化进行对比研究,以便正确地认识和使用地电阻率多极距观测的反演结果.
1. 一维电性结构下地电阻率的反演
3层结构的一维地电阻率模型中,第1,2层厚度分别为7.5 m和20 m.考虑到各层电阻率可能随时间发生周期性变化,因此设定各层电阻率为
(1) 式中: ρi和ρi(0)分别为各层介质的电阻率及其初始值;αi为各层介质电阻率的变化幅度;t为观测时刻,t=1,2,…,24;T为变化周期,T=24.对称四极布置的多极距观测系统的布极参数见表 1.
表 1 一维电性结构下多极距观测系统的布极参数Table 1. The electrodes arrangement parameters of multi-separation array system in one-dimensionAB/m MN/m 30 2 100 16 270 30 500 100 1000 300 研究思路为:首先通过正演获得各测道多极距的模拟观测数据(真电阻率),然后对模拟观测数据进行反演得到各层的电阻率,最后分析各层反演得到的电阻率与真电阻率之间的符合程度.一维反演采用影响系数法(汪雪泉等,2002),该方法的特点是利用两次连续观测数据的差值进行反演.一维地电阻率正演方法已很成熟,这里不再赘述;二维地电阻率正演方法,目前主要有有限元法和积分方程法两类,而本文则采用笔者提出的边界积分方程法(毛先进,鲍光淑,1998).
根据模型所设各层参数及表 1中多极距观测系统的布极参数,通过正演得到每个极距装置系统对应的地表视电阻率曲线,然后利用影响系数法(汪雪泉等,2002)反演得到各层的电阻率(图 1—4中虚线),再与各层的真电阻率(图 1—4中实线)进行比较,两者的误差反映了一维正、反演结果的准确程度.本文主要分以下4种情况进行讨论.
图 3 第1,3层电阻率ρ1和ρ3变化,第2层电阻率ρ2不变时各层的多极距电阻率一维反演结果(虚线)与真电阻率(实线)的比较Figure 3. Comparision of the resistivity by inversion with multi-separation array (dashed line) and the actual resistivity (solid line) for each layer when resistivity of the first-and third-layer ρ1 and ρ3 vary with time, but resistivity of the second-layer ρ2 remains unchanged1) 设定α1=0.03,α2=0,α3=0,T=24,即ρ1变化(变化幅度为3%),ρ2和ρ3不变. 图 1给出了第1—3层电阻率的一维反演结果与真电阻率的时程曲线.可以看出,反演得到的第1—3层电阻率的变化幅度分别为3.10%,0.17%和0.08%,与实际结果非常接近,同时,反演得到的各层电阻率与真电阻率的最大相对误差为1.7/1000,两者的符合度较高.
2) 设定α1=0,α2=0,α3=0.03,T=24,即ρ1和ρ2不变,ρ3变化(变化幅度为3%). 图 2给出了第1—3层电阻率的一维反演结果与真电阻率的时程曲线.可以看出,反演得到的第1—3层的电阻率变化幅度分别为0,0.13%和2.83%,且各层反演得到的电阻率与真电阻率的最大相对误差为2.9/1000,两者的符合度亦较高.
3) 设定α1=0.1,α2=0,α3=0.02,T=24,即ρ1和ρ3变化,变化幅度分别为10%,2%,ρ2不变. 图 3给出了第1—3层电阻率的一维反演结果与真电阻率的时程曲线.可知,反演得到的第1—3层的电阻率变化幅度分别为9.99%,0.25%和1.92%,且各层反演得到的电阻率与真电阻率的最大相对误差为5.4/1000,两者符合程度高.
4) 设定α1=0.3,α2=0.03,α3=0.02,T=24,即ρ1,ρ2,ρ3均变化,变化幅度分别为30%,3%,2%. 图 4给出了第1—3层电阻率的一维反演结果与真电阻率的时程曲线.可以看出,反演得到的第1—3层的电阻率变化幅度分别为31.90%,2.92%和1.95%,各层反演得到的电阻率与真电阻率的最大相对误差为6.7/1000,反演得到的电阻率与真电阻率的符合程度总体较好,并且在第1层真电阻率变化幅度较大(30%),而第2,3层变化幅度较小(分别为3%,2%)的情况下,反演所得结果能够准确地分辨第2,3层电阻率的变化.
以上几种情况下对一维电性结构的地电阻率的正、反演模拟结果表明,对于3层电阻率结构,电阻率发生变化的地层为1—3层不等,但各层反演得到的电阻率与真电阻率的符合程度均较高,这说明在地下一维电性结构下,多极距地电阻率观测的确可以达到区分不同地层电阻率变化的目的.
2. 二维电性结构下地电阻率的反演
目前通常采用一维反演的方法对地电阻率多极距观测视电阻率时间序列进行解释,而多极距观测场地下的电阻率结构可能与一维结构有所差别,换言之,对地电阻率多极距观测结果进行反演解释的理论模型与实际的电性结构可能存在差别.因此,本文考虑台址下更接近实际的二维电阻率结构,先对二维结构进行多极距观测数据的正演,得到模拟的多极距观测的视电阻率时间序列,然后对模拟结果进行一维反演,研究给定的二维电阻率结构与一维反演结果之间的关系.
2.1 二维正演方法的准确性检验
本文采用的二维正演数值模拟方法为作者提出的边界积分方程法(毛先进,鲍光淑,1998),该方法应用于二维电性结构问题的准确性在理论和数值模拟两方面均已得到验证.本文中的准确性检验针对的是一维电阻率结构,分别采用一维方法(解析)与二维边界积分方程法(数值模拟)计算视电阻率,然后考察两者的一致性. 图 5给出了用于二维正演方法准确性检验的一维电阻率结构模型,该模型共分为3层,第1,2层厚度分别为50 m和100 m,第1—3层的电阻率分别为30,60,150 Ω·m. 5套多极距观测系统的电极以图 5所示剖面的地表中心点对称布设,布极参数及一维、二维正演方法得到的地电阻率计算结果列于表 2.
表 2 二维电性结构下多极距观测系统布极参数和针对图 5所示模型利用一维与二维正演方法计算得到的视电阻率结果对比Table 2. The electrodes arrangement parameters of multi-separation array in two dimension and the comparision of apparent resistivities calculated from one-dimensional and two-dimensional forward modeling based on the model of Fig. 5装置编号 AB/m MN/m 一维视电阻率
ρ1D/(Ω·m)二维视电阻率
ρ2D/(Ω·m)相对误差 C1 200 50 38.38 38.20 -0.47% C2 400 150 50.49 50.18 -0.61% C3 600 150 64.24 64.64 0.62% C4 800 250 73.58 74.08 0.68% C5 1000 350 81.46 81.02 -0.32% 注:相对误差=(ρ2D-ρ1D)/ρ1D. 从表 2可见,对于同一模型,二维边界积分方程法与一维正演方法计算得到的视电阻率的最大相对误差的绝对值为0.68%,表明二维边界积分方程法具有足够的精度.
2.2 二维结构下一维反演算例与结果分析
针对不同的二维结构以及不同地层介质电阻率的变化情况,给出几个典型二维电性结构下一维地电阻率反演的算例及其结果分析.
2.2.1 算例1
图 6为根据图 5所示的一维结构模型演化得到的二维电阻率模型,通过改变第1,2层部分单元的电阻率形成层边界有起伏的二维结构.第1,2层厚度分别为50,100 m;第1层下边界的非起伏区域,厚度为50 m,而起伏区厚度为100 m(部分单元伸入第2层),宽度为250 m;第2层下边界的起伏和非起伏区域,厚度均为100 m,起伏区的宽度为250 m(部分单元伸入第3层).
为了模拟前兆观测中普遍存在的地电阻率年变现象,设置浅层(第1,2层)电阻率以不同的幅度作周期性变化,即ρ1(0)=30 Ω·m, ρ2(0)=60 Ω·m, α1=0.03, α2=0.02, 根据式(1) 计算ρ1和ρ2. ρ3恒定为150 Ω·m.经二维正演得到不同时刻的多极距模拟观测数据,之后对这些数据作一维反演.
从图 7a中可以看出,在第1—3层介质真电阻率的变化幅度分别为3%,2%,0时,反演得到的第1—3层介质电阻率的变化幅度分别为3.10%,1.87%,0.25%,各层反演得到的电阻率与真电阻率的变化十分吻合.
2.2.2 算例2
模型与算例1相同,如图 6所示,但浅层电阻率随时间的变化幅度增大,即设ρ1(0)=30 Ω·m, ρ2(0)=60 Ω·m, α1=0.1, α2=0.05, 根据式(1) 计算ρ1和ρ2,ρ3仍设为150 Ω·m.模型第1—3层电阻率的一维反演结果与真电阻率的时程曲线如图 7b所示.
从图 7b可以看出,在第1—3层介质真电阻率的变化幅度分别为10%,5%,0时,反演得到的第1—3层介质电阻率的变化幅度分别为10.30%,5.21%,0.31%,各层反演得到的电阻率与真电阻率的变化情况吻合得较好.
值得注意的是,在算例1和算例2中, 第3层介质的真电阻率是不变的,而反演结果中两个算例第3层介质电阻率的变化幅度分别为0.25%和0.31%,虽然与实际值的偏离均比较小,但存在随浅层(第1,2层)真电阻率变化幅度增大而增大的趋势.因此推测,第1,2层介质的电阻率变化越大,对第3层反演结果的影响也越大.对一定的电阻率结构,当这种变化达到一定程度时,将会影响对第3层(即深部)介质电阻率变化情况的正确判断.
2.2.3 算例3
地电阻率模型如图 8所示,浅层电阻率的参数设置为ρ1(0)=30 Ω·m, ρ2(0)=60 Ω·m, α1=0.3, α2=0.1, 根据式(1) 计算ρ1和ρ2,ρ3仍设为150 Ω·m.
与图 6所示模型不同的是,图 8所示模型的边界起伏只发生在第1,2层之间,第1,2层厚度分别为50,100 m;第1层下边界的非起伏区域层厚为50 m,而起伏区层厚为150 m,达到了第2,3层的交界面上,其宽度为100 m.
图 9为基于图 8所示模型而得到的一维电阻率反演结果与真电阻率时程曲线.可以看出:第1,2层介质的真电阻率变化幅度与算例1,2相比明显增大,变化幅度分别为30%,10%,第3层介质的真电阻率不变;反演得到的第1—3层电阻率的变化幅度分别为31.9%,11.5%和0.52%,反演得到的各层电阻率与真电阻率的变化情况吻合得较好,但反演结果与各层真电阻率的差异均大于前述两个算例,这可能与第1,2层介质的真电阻率变化幅度较大有关,也进一步表明,第1,2层介质的电阻率变化幅度越大,对第3层反演结果的影响越大.
上述3个算列表明,如果地下电性结构并非一维分布,则用多极距观测数据作一维反演获得的不同地层的电阻率与实际电阻率存在差异,且浅层介质电阻率的变化幅度越大,这种差异也越大.
3. 讨论与结论
地电阻率多极距观测的目的,是借助对观测数据的反演获得地下介质中不同地层真电阻率的变化.数据反演是多极距观测方法中的一个重要环节,其准确性和可靠性对地电阻率多极距观测数据的解释与应用十分重要.到目前为止研究人员主要采用一维反演方法对观测数据进行反演处理,由于地下介质的横向不均性,地电阻率观测区域内电阻率结构往往与一维结构存在不同程度的偏离,因此实际地电结构与反演所用的理论模型并不相符,本文初步分析了这种“偏离”或“不符合”对反演结果的影响,得到如下认识:
1) 在一维地下电性结构下,用多极距电阻率观测数据作一维反演可以达到区分不同地层电阻率变化的目的,实践中应尽可能选择地下电阻率结构充分接近一维的场地开展地电阻率多极距观测;
2) 若地下电性结构不满足一维分布,则用多极距观测数据进行一维反演所获得的深部介质的电阻率与实际的深部介质电阻率也会存在差异,且浅层介质的电阻率变化幅度越大,这种差异也越大,反之亦然.
可以认为,定点连续观测中两次连续观测之间的时间间隔越短,则在这段时间内浅部介质的真电阻率变化就越小,因此在观测区的电阻率分布与一维结构存在明显差异的情况下,为了减小反演结果与实际的差异,则应尽量减小两次观测之间的时间间隔.
需要指出的是,本文只是在给定模型下讨论一维结构的多极距观测的一维反演问题,只涉及浅层电阻率周期性变化对反演结果,尤其是研究人员更为关注的深部电阻率反演结果的影响,但实际情况下浅部与深部电阻率的变化组合是十分复杂的,这些复杂情况下的反演结果对深部电阻率变化判断的影响,还有待更为详细地、深入地研究.
-
刁桂苓,徐锡伟, 陈于高, 黄柏寿, 王晓山, 冯向东, 杨雅琼. 2011. 汶川MW7.9和集集MW7.6地震前应力场转换现象及其可能的前兆意义[J]. 地球物理学报, 54(1): 128——136.
杜兴信,邵辉成. 1999.由震源机制解反演中国大陆现代构造应力场[J].地震学报,21(4): 354——360.
李钦祖, 王泽皋, 贾云年, 靳雅敏. 1973. 由单台小地震资料所得两个区域的应力场[J]. 地球物理学报, 16(1): 49——61.
李延兴, 杨国华, 李智, 郭良迁, 黄珹, 朱文耀, 符养, 王琪, 江在森, 王敏. 2003. 中国大陆活动地块的运动与应变状态[J]. 中国科学: D辑,33(刊): 65——81.
马宏生,张国民,刘杰,李丽, 陈化然. 2006. 中国大陆活动地块边界带强震活动特征的研究[J]. 地震地质,28(1):48——60.
康英,杨选,陈杏,陈贵美, 郑斯华.2008.广东及邻区的应力场反演[J]. 地震学报, 30(1): 59——66.
汪素云, 许忠淮. 1985. 中国东部大陆的地震构造应力场[J]. 7(1): 17——22.
谢富仁, 刘光勋. 1989. 阿尔金断裂带中段区域新构造应力场分析[J]. 中国地震, 5(3): 26——36.
谢富仁, 祝景忠, 粱海庆, 刘光勋. 1993. 中国西南地区现代构造应力场基本特征[J].地震学报,15(4): 407——417.
谢富仁, 陈群策, 崔效锋, 李宏, 杨树新, 陈连旺. 2003. 中国大陆地壳应力环境研究[M]. 北京: 地质出版社: 3——12.
谢富仁,崔效锋,赵建涛, 陈群策, 李宏. 2004.中国大陆及邻区现代构造应力场分区[J].地球物理学报,47(4): 654——662.
许忠淮, 阎明,赵仲和. 1983. 由多个小震推断的华北地区构造应力场的方向[J]. 地震学报, 5(3): 268——279.
许忠淮, 戈澍谟. 1984. 用滑动方向拟合法反演富蕴地震断裂带应力场[J]. 地震学报, 8(5): 395——404.
许忠淮, 汪素云, 黄雨蕊, 高阿甲, 金小锋, 常向东. 1987. 由多个小震推断的青甘和川滇地区地壳应力场的方向特征[J]. 地球物理学报, 30(5): 476——486.
许忠淮, 汪素云, 黄雨蕊, 高阿甲. 1989. 由大量地震资料推断我国大陆构造应力场[J]. 地球物理学报, 35(5): 560——572.
许忠淮. 2001. 东亚地区现今构造应力图的编制[J]. 地震学报,23(5):492——501.
薛宏远, 鄢家全. 1984. 鄂尔多斯地块周围的现代地壳应力场[J]. 地球物理学报, 27(2): 144——152.
张国民,张培震.2000.“大陆强震机理与预测”中期学术进展[J].中国基础科学,(10): 4——10.
张国民,马宏生,王辉,王新岭.2005.中国大陆活动地块边界带与强震活动[J]. 地球物理学报,48(3):138——146.
张培震,邓起东,张国民,马瑾, 甘卫军, 闵伟, 毛凤英, 王琪.2003.中国大陆的强震活动与活动地块[J].中国科学: D辑,33(增刊):12——20.
Gephart J W. 1990. FMSI: A fortran program for inverting fault/slickenside and earthquake focal mechanism data to obtain the original stress tensor[J]. Comput Geosci, 16(7): 953——989.
Gephart J W, Forsyth D W. 1984. An improved method for determining the regional stress tensor using focal mechanism data: Application to the San Fernando earthquake sequence[J]. J Geophys Res, 89(B11): 9305——9320.
Wan Y G. 2010. Contemporary tectonic stress field in China[J]. Earthq Sci, 23(4): 377——386.
Zoback M L. 1992. First——and second——order patterns of stress in the lithosphere: The world stress map project[J]. J Geophys Res, 97(B8): 11703——11728.
刁桂苓,徐锡伟, 陈于高, 黄柏寿, 王晓山, 冯向东, 杨雅琼. 2011. 汶川MW7.9和集集MW7.6地震前应力场转换现象及其可能的前兆意义[J]. 地球物理学报, 54(1): 128——136.
杜兴信,邵辉成. 1999.由震源机制解反演中国大陆现代构造应力场[J].地震学报,21(4): 354——360.
李钦祖, 王泽皋, 贾云年, 靳雅敏. 1973. 由单台小地震资料所得两个区域的应力场[J]. 地球物理学报, 16(1): 49——61.
李延兴, 杨国华, 李智, 郭良迁, 黄珹, 朱文耀, 符养, 王琪, 江在森, 王敏. 2003. 中国大陆活动地块的运动与应变状态[J]. 中国科学: D辑,33(刊): 65——81.
马宏生,张国民,刘杰,李丽, 陈化然. 2006. 中国大陆活动地块边界带强震活动特征的研究[J]. 地震地质,28(1):48——60.
康英,杨选,陈杏,陈贵美, 郑斯华.2008.广东及邻区的应力场反演[J]. 地震学报, 30(1): 59——66.
汪素云, 许忠淮. 1985. 中国东部大陆的地震构造应力场[J]. 7(1): 17——22.
谢富仁, 刘光勋. 1989. 阿尔金断裂带中段区域新构造应力场分析[J]. 中国地震, 5(3): 26——36.
谢富仁, 祝景忠, 粱海庆, 刘光勋. 1993. 中国西南地区现代构造应力场基本特征[J].地震学报,15(4): 407——417.
谢富仁, 陈群策, 崔效锋, 李宏, 杨树新, 陈连旺. 2003. 中国大陆地壳应力环境研究[M]. 北京: 地质出版社: 3——12.
谢富仁,崔效锋,赵建涛, 陈群策, 李宏. 2004.中国大陆及邻区现代构造应力场分区[J].地球物理学报,47(4): 654——662.
许忠淮, 阎明,赵仲和. 1983. 由多个小震推断的华北地区构造应力场的方向[J]. 地震学报, 5(3): 268——279.
许忠淮, 戈澍谟. 1984. 用滑动方向拟合法反演富蕴地震断裂带应力场[J]. 地震学报, 8(5): 395——404.
许忠淮, 汪素云, 黄雨蕊, 高阿甲, 金小锋, 常向东. 1987. 由多个小震推断的青甘和川滇地区地壳应力场的方向特征[J]. 地球物理学报, 30(5): 476——486.
许忠淮, 汪素云, 黄雨蕊, 高阿甲. 1989. 由大量地震资料推断我国大陆构造应力场[J]. 地球物理学报, 35(5): 560——572.
许忠淮. 2001. 东亚地区现今构造应力图的编制[J]. 地震学报,23(5):492——501.
薛宏远, 鄢家全. 1984. 鄂尔多斯地块周围的现代地壳应力场[J]. 地球物理学报, 27(2): 144——152.
张国民,张培震.2000.“大陆强震机理与预测”中期学术进展[J].中国基础科学,(10): 4——10.
张国民,马宏生,王辉,王新岭.2005.中国大陆活动地块边界带与强震活动[J]. 地球物理学报,48(3):138——146.
张培震,邓起东,张国民,马瑾, 甘卫军, 闵伟, 毛凤英, 王琪.2003.中国大陆的强震活动与活动地块[J].中国科学: D辑,33(增刊):12——20.
Gephart J W. 1990. FMSI: A fortran program for inverting fault/slickenside and earthquake focal mechanism data to obtain the original stress tensor[J]. Comput Geosci, 16(7): 953——989.
Gephart J W, Forsyth D W. 1984. An improved method for determining the regional stress tensor using focal mechanism data: Application to the San Fernando earthquake sequence[J]. J Geophys Res, 89(B11): 9305——9320.
Wan Y G. 2010. Contemporary tectonic stress field in China[J]. Earthq Sci, 23(4): 377——386.
Zoback M L. 1992. First——and second——order patterns of stress in the lithosphere: The world stress map project[J]. J Geophys Res, 97(B8): 11703——11728.
-
期刊类型引用(2)
1. 叶青,王晓,杜学彬,解滔,范晔,周振贵,刘高川. 中国地震井下地电阻率研究进展. 吉林大学学报(地球科学版). 2022(03): 669-683 . 百度学术
2. 朱涛. 地震电阻率实验研究新进展及展望. 地球与行星物理论评. 2021(01): 61-75 . 百度学术
其他类型引用(3)
计量
- 文章访问数: 1955
- HTML全文浏览量: 1182
- PDF下载量: 325
- 被引次数: 5