Shear wave splitting of local earthquakes in Weiyuan area,Sichuan
-
摘要: 基于布设在四川威远地区的流动地震台站2019年11月至2020年5月记录的近震波形数据,利用剪切波分裂分析方法获得了研究区内23个地震台站的剪切波分裂参数。结果显示:威远地区地震台站快波偏振方向多数为北西向,与区域主压应力方向一致;有个别台站可能由于受到局部构造环境的影响,快波偏振方向为北东向。而在地壳速度变化明显的威远背斜附近,有6个台站结果显示为两个快波偏振优势方向,表明该地区快波偏振方向是区域应力场和局部构造环境共同作用的结果,威远地区慢波延迟时间均值为4.43 ms/km,威远背斜南侧区域的慢波延迟时间值普遍大于其北部区域,表明在威远地区南部区域的各向异性强度强于北部区域。Abstract: In this paper, the shear wave splitting parameters of 23 seismic stations in Weiyuan area, Sichuan are obtained using the shear wave splitting analysis method, base on the local earthquake waveform data recorded by the temporary seismic stations from November 2019 to May 2020 in the study area. The results show that the polarization directions of fast shear wave at most stations in Weiyuan area are northwest, which is consistent with the direction of regional principal compressive stress. And the polarization directions of a few stations are northeast, which may be affected by the local tectonic environment. Near Weiyuan anticline with obvious differences in crustal velocity, there are six stations with two dominant polarization directions of fast wave, indicating that the polarization direction of fast wave in this area is influenced by the combined effects of regional stress field and local tectonic environment. The average delay time of slow wave in Weiyuan area is 4.43 ms/km. And the delay time of slow wave in the south of Weiyuan anticline is generally greater than that in the north, suggesting that the anisotropy intensity in the south of Weiyuan area is stronger than that in the north.
-
引言
龙门山断裂带位于青藏高原东缘、四川盆地西北缘,长约500 km,宽约30—50 km,沿NE−SW向展布,断层滑动以逆冲为主兼具右旋走滑分量。如图1所示,龙门山断裂带的主要断裂自西向东分布有后山断裂(汶川—茂县断裂)、中央断裂(映秀—北川断裂)、前山断裂(江油—灌县断裂)(唐荣昌等,1991;邓起东等,1994;贾东等,2003)。在该断裂带上相继发生了2008年汶川MS8.0地震和2013年芦山MS7.0地震,两次地震的地震成因和该断裂带的后续地震危险性研究备受关注。构建三维地质模型能够实现浅层地质调查和深部构造资料的综合,对岩石圈(地壳)结构进行直观描述,是进行孕震环境和地震机制分析研究的前提。
图 1 龙门山断裂带地质构造背景及剖面分布Ⅰ :竹巴龙—资中剖面;Ⅱ :花石峡—简阳剖面;Ⅲ :黑水—三台剖面;Ⅳ :唐克—蒲江—阆中三角剖面. F1-1:汶川—茂汶断裂;F1-2:平武—青川断裂;F2:映秀—北川断裂;F3:灌县—安县断裂Figure 1. The geological structure of Longmenshan fault zone and profile positionⅠ :Zhubalong-Zizhong profile;Ⅱ :Huashixia-Jianyang profile;Ⅲ :Heishui-Santai profile;Ⅳ :Tangke-Pujiang-Langzhong profile. F1-1:Wenchuan-Maowen fault;F1-2:Pingwu-Qingchuan fault; F2:Yingxiu-Beichuan fault;F3:Guanxian-Anxian fault近二三十年关于龙门山断裂带的研究已经累积了丰富的地球物理勘探资料,如深地震测深剖面(陈学波等,1988;韩渭宾,蒋国芳,2000;Wang et al,2007;滕吉文等,2008;徐朝繁等,2008;朱介寿,2008;Li et al,2009;Jia et al,2014;王夫运等,2015;张新彦等,2017)、基于接收函数反演或层析成像的速度结构模型构建(王椿镛等,2002,2003;刘启元等,2009;郭飚等,2009;雷建设等,2009;吴建平等,2009;Li et al,2009;邓文泽等,2014)、基于重磁数据的密度模型反演(邓起东,1994;滕吉文等,2008;张季生等,2009;唐新功等,2012;张恩会等,2015;杨光亮等,2015)以及大地电磁探测(李立,金国元,1987;孙洁等,2003;陈高等,2006;徐朝繁等,2008;王绪本等,2008,2009)。
龙门山区域的研究成果虽然丰富,但呈现单一分散的状态,针对这些已有成果进行客观整理的报道甚少。为此,本文拟将龙门山区域具有代表性的二维剖面作为数据源,通过专业的三维建模软件构建该区域深达70 km的三维地质模型,并通过实测的布格重力异常数据检验模型,继而通过重力资料约束反演对该模型进一步予以修正,为有效预测龙门山未知区域提供参考。
1. 构造背景
龙门山断裂带处于华北地块、华南地块与羌塘地块的结合部位,西北侧是松潘—甘孜古残留洋盆,东南侧是稳定的扬子克拉通(Jia et al,2014)。龙门山断裂带及邻近区域的速度结构反演结果表明:该区域的速度结构从青藏高原东南部的山地到扬子地台具有强烈的不均匀性;四川盆地所处沉积层较厚,且其沉积环境自晚古生代以来变化不大(王椿镛等,2002),地壳结构整体稳定,地壳及上地幔的平均速度均明显偏高(王椿镛等,2015),其中上地壳呈现低速异常,表明该地区处于坳陷状态(胥颐,2009);松潘—甘孜地块地壳因介质岩性变化而呈增厚状态(王椿镛等,2015),其上地壳速度较高,中下地壳速度偏低,上地幔也表现为明显低速异常(胥颐,2009;王椿镛等,2015)。
2. 数据与方法
三维地质建模是一门通过空间信息管理、分析与地穴统计技术构造三维地质体,并进行地质解释的多学科教程技术(潘懋等,2007)。地质单元具有多尺度的非均匀性和连续性,根据地质单元的连续性可以从测得的部分信息推算其相邻区域的趋势。但是,这种推算方法的主观性较强,建模过程会掺杂建模人员的主观判断,因此建模的结果具有多解性,即模型的随机性。
建模是基于一定的算法将二维资料推演成三维空间信息的过程。首先,在模型表面上确定剖面与断裂带的位置,在模型表面约束剖面的地表迹线、断层的走向及影响范围;然后,由剖面迹线建立对应的剖面,对剖面信息进行数字化,剖面信息包含断层在剖面位置向下延展轨迹数据、层位在该剖面的具体埋深及形态信息;最后,依据适合的算法,由二维剖面生成三维模型。实际建模过程中需要人为干预,协调各剖面信息以确保数据的一致性。对模型赋予物性参数并进行正演,与实测的重力数据进行对比,从而验证模型的可信度。
本文使用GeoModeller软件(4.0.5版)进行三维地质建模,该软件支持输入不同类型的数据,如二维剖面、地质图、钻孔数据等,支持对数据的可视化和编辑,定义不同地质单元的接触关系后,使用势场插值的方式构建三维模型(Calcagno et al,2011)。本文取研究区范围为(102°E—107°E,29°N—34°N),建模内容包括断层模型和地层模型。用于比对的布格重力原数据是(98°E—108°E,26°N—34°N)范围内的采样间隔为5 km的川滇地区格网数据。
2.1 断层参数
用于三维建模的断层参数包括地表迹线和断裂倾向、倾角。断层地表迹线根据模型分辨率进行了适当的简化。
本文只考虑前山断裂带、中央断裂带和后山断裂带,忽略其它断层的影响。龙门山后山断裂包括北段的平武—青川断裂(F1-2)和南段的汶川—茂汶断裂(F1-1)。平武—青川断裂走向N60°—70°E (贾东等,2003);汶川—茂汶断裂中段走向为N25°—45°E,大约10 km深度以上产状的倾向为NW,倾角约为55°,在5 km深度处开始变缓,倾角约为35°,在北端与平武—青川断裂斜交(唐荣昌等,1991;马保起等,2005;冯杨洋,2016),汶川—茂汶断裂南段的走向为N45°E (徐锡伟等,2008)。
龙门山中央断裂,即映秀—北川断裂,分为北段和南段两部分。北川—映秀断裂北段(F2-2)的走向为N45°E,倾向NW,整体倾角约为50°—70°,在平均19 km深度以上倾向NW,倾角约为55°,在6 km深度处开始变缓;北川—映秀断裂南段(F2-1)的走向为N35°—45°E,倾向为NW,平均17 km深度以上倾向NW,倾角约为55°,在9 km深度处开始变缓(邓起东等,1994;赵小麟等,1994;冯杨洋,2016)。
龙门山前山断裂(F3)又称灌县—安县断裂,总体走向呈N35°—45°E,断面倾向NW,倾角为50°—70°,由南到北依次为F3-1,F3-2,F3-3,F3-4:南段的F3-1和F3-2走向为N43°E,倾向为NW,倾角45°—65°不等,在10 km深度以上倾向NW,倾角约为50°,在4 km处开始变缓,倾角为25°;中段F3-3 (灌县—江油断裂)主要发育在中生代地层中,倾向NW,倾角较陡,产状为北东陡南西缓;北段F3-4在平均14 km深度以上倾向NW,倾角约为45°,在8 km处开始变缓(邓起东等,1994;冯杨洋,2016)。
2.2 地层参数
地层参数包括地层界面形态和地层物性参数。地层密度根据密度与速度的转化公式(冯锐,1985)和实测剖面的地震波速度推算,结合重力反演最终确定,详见表1。地层界面形态根据收集到的剖面获取,经过筛选,最终使用以下四个剖面作为建模数据源:唐克—蒲江—阆中三角测线DSS剖面Ⅳ(陈学波等,1988;邓阳凡等,2011)、黑水—三台剖面Ⅲ(邓阳凡等,2011)、花石峡—简阳剖面Ⅱ(王有学等,2005)和竹巴龙—资中剖面Ⅰ(王椿镛,2002,2003;Jia et al,2014),各剖面的地表迹线见图1,图2为经过修改和数字化的剖面。
表 1 地层密度Table 1. Formation density地层 密度/(g·cm−3) 沉积层 2.36 上地壳 2.67 中地壳 2.9 下地壳 3.14 地幔 3.35 图 2 本文采用的五个数字化剖面(a) 竹巴龙—资中剖面(改自Wang et al,2007);(b) 唐克—蒲江剖面(改自韩渭宾,蒋国芳,2000);(c) 蒲江—阆中剖面(改自韩渭宾,蒋国芳,2000);(d) 阆中—唐克剖面(改自韩渭宾,蒋国芳,2000);(e) 花石峡—简阳剖面(改自王有学等,2005);(f) 黑水—三台剖面(改自王成,2012)Figure 2. Five digital profiles used in this paper(a) Zhubalong-Zizhong profile (revised from Wang et al,2007);(b) Tangke-Pujian profile (revised from Han and Jiang,2000);(c) Pujiang-Langzhong profile (revised from Han and Jiang,2000);(d) Langzhong-Tangke profile (revised from Han and Jiang,2000);(e) Huashixia-Jianyang profile (revised from Wang et al,2005);(f) Heishui-Santai profile (revised from Wang,2012)2.3 重力数据
本文重力数据购买于全国地质资料馆。中国地质调查局发展研究中心于2000年对原有的各行业间的数据进行整合,按照统一的重力基本网(1985)、坐标系(北京54)和国家高程标准(1985)、正常重力场公式、地改半径和中间层密度,建立了规范的、安全的全国范围重力数据库。该数据库收录原地矿部共计255个1:20万图幅的1:20万数据、84个1:20万图幅的1:50万数据、27个1:100万图幅的1:100万数据以及1999年以来 “国土资源大调查” 专项部署完成的青藏高原6个1:100万图幅和遍布全国的78个1:20万图幅的区域重力调查成果数据。
2.4 方法
本文利用Geomodeller进行建模。该软件支持多种地质数据的导入、可视化与编辑,在定义了描述不同地层单元之间关系的地层柱状图后,利用势场插值法计算三维模型。Geomodeller基于势场理论进行三维插值,使用协同克里金算法计算接触点(地质、断层)和方向数据,并将其视作三维标量场的等势面增量和导数(Lajaunie et al,1997)。根据不同接触关系定义地层层序,定义曲面为不同场的特殊等值面,根据地层层序的定义生成三维地质模型(Calcagno et al,2011)。
地球物理正演是指在地球物理资料解释理论中,由地质体的形状、产状、空间位置等赋存状态和密度、磁性、电性、弹性、速度等物性参数计算该地质体引起的场异常过程。正演模型是利用计算机模拟方法来获得已知地质体的地震响应特征,从而指导地震资料的构造和岩性解释。与反演相比,正演的优势在于改变模型参数方便,选用的计算方法灵活,模型制作简单。
正演算法采用多线程空间卷积算法:给定分辨率,将连续的地质模型离散化为多个正则或半正则三维体元;根据已知的观测值计算每个体元的权重,并将其作为卷积算子;对模型进行空间卷积,统计不同体元的权重总和,得到模型各点的重力异常;评估每一次正演结果的残差,迭代上述计算,当残差值达到最小即停止迭代(Gibson et al,2013)。
3. 讨论与分析
3.1 模型及重力异常残差
最终建成的地层模型和断层模型分别如图3a和图3b所示。由图3a可见,沉积层整体较薄,受分辨率的影响并未表现出较大的变化;地壳层以龙门山断裂带为界,两侧有明显的厚度差异,其中下地壳贡献的厚度变化值最大,上地壳次之,中地壳厚度变化最不明显;本模型的地幔层主要反映上地幔,在龙门山断裂带西北侧的松潘—甘孜地块埋深较大,经过龙门山断裂带的抬升,在东南侧的四川盆地上地幔约有10—20 km的抬升。
实测剖面(王成,2012)显示,位于四川盆地的龙泉驿、成都、彭州至都江堰段莫霍面埋深处于45—48 km之间,郸县、双流、资阳、简阳、成都至汉川段莫霍面埋深介于42—44 km之间;位于龙门山断裂带的耿达(埋深48 km)处莫霍面陡降至达维的68 km,经小金、丹巴缓慢上升至64 km,越过鲜水河断裂,道孚附近的莫霍面变浅至60 km,在茂县(45 km)陡降至黑水的58 km,莫霍面出现了13 km的陡降带,随后,经壤口、红原、唐克、麦溪缓慢上升至50 km。模型整体趋势与实测数据吻合。
对模型进行正演以验证其可靠性,参数见表1。布格重力异常原件的重力值范围为−76.17×10−5—−435.15×10−5 m/s2 (图4a),计算的重力值范围为−172.88×10−5—−424.23×10−5 m/s2 (图4b)。根据总体形态,计算场与观测场具有一致性,能反映断裂沿东北—西南对角线分布,但根据数值分析,残差较大(图4c)。如图所示,观测值西侧边缘表现为带状的重力低值不规则圈闭,计算值表现为存在于西北角的小面积重力低值圈闭,根据竹巴龙—资中测线(Wang et al,2007)的解释剖面,本文推断西南角的重力低值与埋深17—20 km处存在的延伸至龙门山断裂带的低速异常体有关;南侧中部偏北的重力低值与花石峡—简阳剖面(王有学等,2005)上埋深25—40 km处的低速异常体、阿坝—双流测线(朱介寿,2008)上埋深20—30 km处低速异常体有关。至于龙门山断裂带,观测值中显示出一条明显的重力梯度带,与计算值形态基本吻合,但断裂带西南端和东北端的重力梯度带表现为明显的向内侧收缩,出现这种差异我们推测原因主要有两点,其一是由于断裂带两端缺少实测剖面数据约束,使得曲面插值算法在模型边缘产生了误差,其二是模型正演在边缘处存在一定的误差。
根据实际建模中遇到的误差问题,总结如下:① 剖面数据由人工反演获取,反演的算法与人工解释的过程均存在一定的误差;② 在人工地震测深等试验下,剖面位置均不是规则的线,后期处理中经过一定程度的加工;③ 模型建造中需要将不同剖面进行匹配,为了保证整个模型的一致,需要人为修改部分剖面上的信息。
从整体来看,模型正演得到的重力异常展布与观测值基本一致,即整体重力值呈现为自西北至东南单调上升的重力异常梯度带,龙门山断裂带体现为高梯度变化带,这表明速度剖面构造的正演模型能很好地反映重力场的变化趋势。
3.2 各地层正演结果
通过调整物性参数的设置,排除其它各层的影响从而实现各地层的单独正演,结果如图5所示。
1) 沉积层。沉积层引起的布格重力异常整体表现为−45.53×10−5—−444.63×10−5 m/s2范围内的负异常,西北高、东南低,分区明显。沿龙门山断裂带形成了一条重力梯度带,该梯度带将整个区域划分为东西两侧,松潘—甘孜盆地沉积层薄、四川盆地沉积层厚。
2) 中上地壳。上地壳和中地壳皆显示为负异常值,前者变化范围较大,在−49.29×10−5—−352.03×10−5 m/s2之间,后者变化范围较小,在−89.14×10−5—−117.68×10−5 m/s2;两者形态差异较大,上地壳大致呈现为自西北至东南的梯度带,龙门山断裂带表现为明显且粗糙的高梯度变化带,而中地壳呈现为多个无规律的圈闭,或重力高圈闭,或重力低圈闭。
上地壳西北角表现为低重力异常,与奔子栏—唐克埋深在10—25 km范围内速度为5.7—5.9 km/s的低速异常体有关(王椿镛等,2003);东侧中部的阿坝至理县(朱介寿,2008)、花石峡附近(王有学等,2005;邓阳凡等,2011)及黑水附近(蔡学林,2007)的中地壳存在低速异常体,与正演结果所显示的重力低值圈闭相对应。
3) 下地壳。较中上地壳和沉积层,下地壳引起的布格重力异常变化趋于光滑,整体为正异常,西北高、东南低,范围为26.37×10−5—52.73×10−5 m/s2,变化较小,龙门山断裂带表现为光滑的异常梯度带。下地壳异常变化小,说明厚度变化小,整体厚度较大,平均厚度达30 km。
4) 莫霍面。莫霍面的重力异常逐渐趋于简单,局部重力异常现象消失,反映了上地幔速度结构的简单化。重力异常沿龙门山断裂带表现为一条重力梯度带,梯度值介于75.08×10−5—257.82×10−5 m/s2。模型内布格重力异常数据呈现为规则的自西北至东南单调上升的重力异常梯度带。地幔层的正演结果与利用重力数据反演获取的4阶小波逼近图(深度等同于地幔层)的趋势一致(姜文亮,张景发,2011)。
龙门山断裂带西北侧的松潘—甘孜地块夹持于昆仑断裂带与龙门山断裂带之间,莫霍面在整体上呈凹陷形态,且地壳自西北向东南方向逐渐变薄,表现为台阶状抬升,莫霍面深度为60—70 km,东南侧四川盆地的莫霍面深度经龙门山的抬升,其埋深陡增至40—50 km (王有学等,2005)。
3.3 模型反演
为降低残差,在重力数据约束下对模型进行反演,根据残差计算(图6a−c)可知反演结果良好,图6d−f为反演后的模型展示,沿断层走向莫霍面的抬升情况很明显,与上述分析一致。
4. 讨论与结论
本文挑选整合了龙门山区域的经典文献研究资料,从中挑选了部分二维剖面,经过数字化、检校编辑、地层定义等步骤,构建了龙门山区域的三维地质模型,利用该区域实测重力数据验证,结果证明了模型的可用性;之后通过重力数据反演,修正了模型的误差,形成了龙门山区域基本可用的三维地质模型。
模型显示,龙门山区域的莫霍面埋深约40—66 km,分区明显,西低东高、北低南高;下地壳厚度处于9—31 km的变化范围,西北端厚,向东南端逐渐变薄;中地壳厚度变化无明显规律,存在明显的隆起与凹陷,整体厚度值在7—19 km;上地壳厚度约为22 km,自南至北呈增厚趋势,其中部分区域存在隆起与凹陷;沉积层埋深呈现东北浅、西南深的分布状态,埋深范围为1—9 km。其中中上地壳存在的隆起、凹陷与前人研究中所揭示的低速异常体吻合,地幔的变化趋势与前人反演的自东北向西南单调上升重力梯度带吻合。
本文通过建模对龙门山区域的地球物理勘探成果进行了有选择的统一整理,试图综合管理现有的成果并以三维的形式展示,是对地球物理数据统一管理与共享方式的有益探索。模型的精度受原数据精度限制,未能充分还原地壳结构的细节,建模过程忽略了一些特殊地质体,若需对特殊地质体进行建模需要高精度数据支撑;建模过程中主要存在不同剖面数据间的不一致性,受建模人员的主观影响较大,这是三维地质建模下一步需要解决的问题。
-
图 2 威远地区构造背景、台站和2019年11月—2020年5月地震分布
F1:牛颈垇断层;F2:大垭口断层;F3:东兴场断层;F4:高桥断层;F5:资中断层;F6:墨林场断层;F7:重滩断层;F8:黄桷坡断层;F9:长岩断层,下同
Figure 2. Tectonic settings,the distribution of seismic stations and the epicenters between November 2019 and May 2020 in Weiyuan region
F1:Niujing’ao fault;F2:Dayakou fault;F3:Dongxingchang fault;F4:Gaoqiao fault;F5:Zizhong fault;F6:Molinchang fault;F7:Chongtan fault; F8:Huangjuepo fault;F9:Changyan fault,the same below
图 3 地震波形及偏振分析示例
(a) 台站SC38记录的2020-12-06 17:05:08发生的ML2.1地震波形;(b) 台站SC12记录的2021-05-11 21:40:46发生的ML1.5地震波形;(c) 台站SC06记录的2021-04-07 11:26:42发生的ML1.9地震波形。左图为三分量地震波形,t1,t2,t3为用于偏振分析的剪切波时间段;右图为三个时段剪切波质点偏振图,S1和S2分别为快剪切波和慢剪切波的到时
Figure 3. Seismic waveforms and diagrams of the polarization analysis
(a) Seismic waves of the ML2.1 earthquake at 17:05:08 on December 6,2020 recorded by the station SC38;(b) Seismic waves of the ML1.5 earthquake at 21:40:46 on May 11,2021 recorded by the station SC12;;(c) Seismic waves of the ML1.9 earthquake at 11:26:42 on April 7,2021 recorded by the station SC06 ;and the left panels are the three-component records of seismic waveform. t1,t2 and t3 are the time periods of shear wave for polarization analysis,respectively. The right panels are the particle trails of the shear waves in three time periods. S1 and S2 are the arrival times of fast and slow shear waves,respectively.
表 1 威远地区各台站剪切波分裂参数结果
Table 1 The results of shear wave splitting parameters at the seismic stations in Weiyuan area
台站
代码有效事件
个数优势方向
个数快波偏振
方向/°慢波延迟时间
/(ms·km−1)台站
代码有效事件
个数优势方向
个数快波偏振
方向/°慢波延迟时间
/(ms·km−1)SC01 34 1 −39±17.8 4.85±1.21 SC15 78 1 −42±18.3 5.23±1.72 SC02 32 1 −34±30.1 5.12±1.04 SC16 43 1 −46±15.7 4.66±1.38 SC03 5 1 −55±4.0 3.74±0.85 46 2 51±13.2 4.73±1.19 SC05 17 1 −46±20.8 3.98±1.01 SC17 35 1 −41±16.8 4.49±1.54 SC06 38 1 −49±10.9 4.13±1.28 SC19 15 1 −17±23.7 3.31±0.36 11 2 37±18.3 3.99±1.06 SC20 7 1 −39±11.7 3.41±0.71 SC07 36 1 20±14.6 4.35±1.13 SC21 6 1 −36±12.9 3.48±0.76 SC08 41 1 −46±15.3 4.76±1.42 SC22 31 1 −33±25.6 3.65±1.03 12 2 15±9.9 4.21±1.16 SC27 67 1 −48±14.1 4.87±1.20 SC09 16 1 −39±11.6 4.40±1.02 SC28 48 1 −44±18.2 4.34±0.68 7 2 25±17.5 3.89±0.88 SC30 5 1 −57±13.5 4.92±1.38 SC10 13 1 −43±16.1 3.37±1.18 SC38 53 1 −45±14.7 4.52±1.15 SC11 15 1 −31±20.5 4.24±1.61 14 2 56±11.8 4.33±1.07 SC12 24 1 −42±15.8 5.06±1.80 SC39 47 1 −36±22.3 4.37±1.06 17 2 44±13.2 4.62±1.21 -
白卓立,张学刚,薛二峰,张希,李瑞莎,秦珊兰. 2020. 基于GPS速度场的四川地区1999—2018年应变场特征演化分析[J]. 地质科学,55(4):1315–1326. doi: 10.12017/dzkx.2020.079 Bai Z L,Zhang X G,Xue E F,Zhang X,Li R S,Qin S L. 2020. The characteristic evolution analysis of scale strain field in Sichuan area from 1999−2018 based on GPS velocity field[J]. Chinese Journal of Geology,55(4):1315–1326 (in Chinese).
常利军,王椿镛,丁志峰. 2008. 四川及邻区上地幔各向异性研究[J]. 中国科学:D辑,28(12):1589–1599. Chang L J,Wang C Y,Ding Z F. 2008. Study on upper mantle anisotropy in Sichuan and its adjacent areas[J]. Science in China:Series D,28(12):1589–1599 (in Chinese).
陈安国,高原,石玉涛. 2019. 龙门山断裂带域上地壳各向异性及其变化[J]. 地球物理学报,62(8):2959–2981. doi: 10.6038/cjg2019M0179 Chen A G,Gao Y,Shi Y T. 2019. Seismic anisotropy and its variation in the upper crust beneath the Longmenshan fault zone[J]. Chinese Journal of Geophysics,62(8):2959–2981 (in Chinese).
董大忠,高世葵,黄金亮,管全中,王淑芳,王玉满. 2014. 论四川盆地页岩气资源勘探开发前景[J]. 天然气工业,34(12):1–15. Dong D Z,Gao S K,Huang J L,Guan Q Z,Wang S F,Wang Y M. 2014. A discussion on the shale gas exploration & development prospect in the Sichuan basin[J]. Natural Gas Industry,34(12):1–15 (in Chinese).
杜广宝,吴庆举,张雪梅,张瑞青. 2021. 四川威远及邻区中小地震活动特征及地壳精细结构研究[J]. 地球物理学报,64(11):3983–3996. doi: 10.6038/cjg2021O0452 Du G B,Wu Q J,Zhang X M,Zhang R Q. 2021. Activity characteristics of moderate and small earthquakes and fine crustal feature beneath Weiyuan,Sichuan and adjacent areas[J]. Chinese Journal of Geophysics,64(11):3983–3996 (in Chinese).
郭铁龙,高原. 2020. 剪切波分裂揭示的青藏高原上地壳地震各向异性基本特征[J]. 地球物理学报,63(3):1085–1103. doi: 10.6038/cjg2020N0156 Guo T L,Gao Y. 2020. Seismic anisotropy in the upper crust within Tibetan Plateau revealed by shear-wave splitting[J]. Chinese Journal of Geophysics,63(3):1085–1103 (in Chinese).
何超红,谭孝成,黄成智. 2015. 威远县地质灾害发育规律及成因分析[J]. 防灾科技学院学报,17(1):7–11. doi: 10.3969/j.issn.1673-8047.2015.01.002 He C H,Tan X C,Huang C Z. 2015. Development patterns and formation mechanisms of geological hazards in Weiyuan county[J]. Journal of Institute of Disaster Prevention,17(1):7–11 (in Chinese).
黄周传,吉聪,吴寒婷,石宇通,耿嘉琪,徐弥坚,韩存瑞,徐鸣洁,王良书. 2021. 青藏高原东南缘地壳结构与变形机制研究进展[J]. 地球与行星物理论评,52(3):291–307. doi: 10.19975/j.dqyxx.2021-005 Huang Z C,Ji C,Wu H T,Shi Y T,Geng J Q,Xu M J,Han C R,Xu M J,Wang L S. 2021. Review on the crustal structures and deformations in the southeastern margin of the Tibetan Plateau[J]. Reviews of Geophysics and Planetary Physics,52(3):291–307 (in Chinese).
李长军,甘卫军,秦姗兰,郝明,梁诗明,杨帆. 2019. 青藏高原东南缘南段现今变形特征研究[J]. 地球物理学报,62(12):4540–4553. doi: 10.6038/cjg2019M0692 Li Z J,Gan W J,Qin S L,Hao M,Liang S M,Yang F. 2019. Present-day deformation characteristics of the southeast borderland of the Tibetan Plateau[J]. Chinese Journal of Geophysics,62(12):4540–4553 (in Chinese).
刘莎,吴朋. 2015. 紫坪铺水库水位变化对剪切波分裂参数的影响[J]. 地球物理学报,58(11):4106–4114. doi: 10.6038/cjg20151118 Liu S,Wu P. 2015. The effect of water level changes in Zipingpu reservoir on the parameters of shear wave splitting[J]. Chinese Journal of Geophysics,58(11):4106–4114 (in Chinese).
刘莎,吴朋,杨建思,苏金蓉. 2015. 2012年6月24日宁蒗—盐源MS5.7地震的剪切波分裂探讨[J]. 地震学报,37(5):787–795. doi: 10.11939/jass.2015.05.007 Liu S,Wu P,Yang J S,Su J R. 2015. Discussion on shear-wave splitting of the June 24,2012 Ninglang-Yanyuan MS5.7 earthquake[J]. Acta Seismologica Sinica,37(5):787–795 (in Chinese).
刘顺,罗志立,熊荣国,宋颖辉. 2000. 从实验岩石变形过程探讨四川盆地加里东古隆起的形成机制[J]. 成都理工学院学报,27(4):343–347. Liu S,Luo Z L,Xiong R G,Song Y H. 2000. On formation mechanism of Caledonian Paleoswell in Sichuan basin from experimental rocks deformational process[J]. Journal of Chengdu University of Technology,27(4):343–347 (in Chinese).
刘顺. 2001. 四川盆地威远背斜的形成时代及形成机制[J]. 成都理工学院学报,28(4):340–343. Liu S. 2001. Formation time and formation mechanism of the Weiyuan anticline in Sichuan basin[J]. Journal of Chengdu University of Technology,28(4):340–343 (in Chinese).
石强,陈鹏,王秀芹,刘凤新. 2017. 页岩气水平井高产层段判识方法及其应用:以四川盆地威远页岩气示范区下志留统龙马溪组为例[J]. 天然气工业,37(1):60–65. doi: 10.3787/j.issn.1000-0976.2017.01.007 Shi Q,Chen P,Wang X Q,Liu F X. 2017. A method for identifying high-productivity intervals in a horizontal shale gas well and its application:A case study of the Lower Silurian Longmaxi Fm in Weiyuan-shale gas demonstration area,Sichuan basin[J]. Natural Gas Industry,37(1):60–65 (in Chinese).
石玉涛,高原,赵翠萍,姚志祥,太龄雪,张永久. 2009. 汶川地震余震序列的地震各向异性[J]. 地球物理学报,52(2):398–407. Shi Y T,Gao Y,Zhao C P,Yao Z X,Tai L X,Zhang Y J. 2009. A study of seismic anisotropy of Wenchuan earthquake sequence[J]. Chinese Journal of Geophysics,52(2):398–407 (in Chinese).
石玉涛,高原,张永久,王辉,姚志祥. 2013. 松潘—甘孜地块东部、川滇地块北部与四川盆地西部的地壳剪切波分裂[J]. 地球物理学报,56(2):481–494. Shi Y T,Gao Y,Zhang Y J,Wang H,Yao Z X. 2013. Shear-wave splitting in the crust in eastern Songpan-Garzê block,Sichuan-Yunnan block and western Sichuan basin[J]. Chinese Journal of Geophysics,56(2):481–494 (in Chinese).
唐荣昌,黄祖智,马声浩,龚宇,周荣军. 1995. 四川活动断裂带的基本特征[J]. 地震地质,17(4):390–396. Tang R C,Huang Z Z,Ma S H,Gong Y,Zhou R J. 1995. Basic characteristics of active fault zones in Sichuan Province[J]. Seismology and Geology,17(4):390–396 (in Chinese).
王雪鹤,李永华,呼楠. 2021. 腾冲地区近震S波分裂研究[J]. 地球物理学报,64(1):131–145. doi: 10.6038/cjg2021O0098 Wang X H,Li Y H,Hu N. 2021. Shear wave splitting analysis of local earthquakes at Tengchong[J]. Chinese Journal of Geophysics,64(1):131–145 (in Chinese).
吴晶,高原,陈运泰,黄金莉. 2007. 首都圈西北部地区地壳介质地震各向异性特征初步研究[J]. 地球物理学报,50(1):209–220. Wu J,Gao Y,Chen Y T,Huang J L. 2007. Seismic anisotropy in the crust in northwestern capital area of China[J]. Chinese Journal of Geophysics,50(1):209–220 (in Chinese).
谢军. 2018. 长宁—威远国家级页岩气示范区建设实践与成效[J]. 天然气工业,38(2):1–7. doi: 10.3787/j.issn.1000-0976.2018.02.001 Xie J. 2018. Practices and achievements of the Changning-Weiyuan shale gas national demonstration project construction[J]. Natural Gas Industry,38(2):1–7 (in Chinese).
徐炳高. 1997. 四川省威远地区燕山早期的构造反转及威远背斜雏形的形成[J]. 四川地质学报,17(2):96–102. Xu B G. 1997. Tectonic inversion in Weiyuan region,Sichuan in early Yanshanian epoch and formation of Weiyuan anticline embryo[J]. Acta Geologica Sichuan,17(2):96–102 (in Chinese).
杨妍. 2019. 川滇地区地壳速度结构及方位各向异性的联合反演研究[D]. 合肥: 中国科学技术大学: 29–47. Yang Y. 2019. Joint Inversion for Crustal Structure and Azimuthal Anisotropy: Application to the Sichuan−Yunnan Region, SW China[D]. Hefei: University of Science and Technology of China: 29–47 (in Chinese).
易桂喜,龙锋,梁明剑,赵敏,王思维. 2020. 四川盆地荣县—威远—资中地区发震构造几何结构与构造变形特征:基于震源机制解的认识和启示[J]. 地球物理学报,63(9):3275–3291. doi: 10.6038/cjg2020O0095 Yi G X,Long F,Liang M J,Zhao M,Wang S W. 2020. Geometry and tectonic deformation of seismogenic structures in the Rongxian-Weiyuan-Zizhong region,Sichuan basin:Insights from focal mechanism solutions[J]. Chinese Journal of Geophysics,63(9):3275–3291 (in Chinese).
易桂喜,龙锋,梁明剑,赵敏,王思维,宫悦,乔慧珍,苏金蓉. 2019. 2019年6月17日四川长宁MS6.0地震序列震源机制解与发震构造分析[J]. 地球物理学报,62(9):3432–3447. Yi G X,Long F,Liang M J,Zhao M,Wang S W,Gong Y,Qiao H Z,Su J R. 2019. Focal mechanism solutions and seismogenic structure of the 17 June 2019 MS6.0 Sichuan Changning earthquake sequence[J]. Chinese Journal of Geophysics,62(9):3432–3447 (in Chinese).
喻颐,李忠权,杨渊宇,应文峰,杭文艳. 2013. 四川盆地威远地区构造演化特征及其对下古生界油气富集的控制作用[J]. 天然气勘探与开发,36(2):1–4. doi: 10.3969/j.issn.1673-3177.2013.02.001 Yu Y,Li Z Q,Yang Y Y,Ying W F,Hang W Y. 2013. Characteristics of tectonic evolution of Weiyuan area in Sichuan basin and its effect on lower Paleozoic oil and gas reservoirs[J]. Natural Gas Exploration and Development,36(2):1–4 (in Chinese).
曾求,储日升,盛敏汉,危自根. 2020. 基于地震背景噪声的四川威远地区浅层速度结构成像研究[J]. 地球物理学报,63(3):944–955. Zeng Q,Chu R S,Sheng M H,Wei Z G. 2020. Seismic ambient noise tomography for shallow velocity structures beneath Weiyuan,Sichuan[J]. Chinese Journal of Geophysics,63(3):944–955 (in Chinese).
张国苓,杨晓松,陈建业,闫小兵. 2010. 中下地壳岩石弹性波各向异性的影响因素[J]. 地震地质,32(2):327–337. Zhang G L,Yang X S,Chen J Y,Yan X B. 2010. The influencing factor of elastic anisotropy in middle to lower continental crust[J]. Seismology and Geology,32(2):327–337 (in Chinese).
张永久,高原,石玉涛,程万正. 2008. 四川区域地震台网的剪切波分裂研究[J]. 地震学报,30(2):123–134. Zhang Y J,Gao Y,Shi Y T,Cheng W Z. 2008. A study on S wave splitting using waveform data from Sichuan regional seismic network[J]. Acta Seismologica Sinica,30(2):123–134 (in Chinese).
张永久,高原,石玉涛,太龄雪. 2010. 四川紫坪铺水库库区地震剪切波分裂研究[J]. 地球物理学报,53(9):2091–2101. Zhang Y J,Gao Y,Shi Y T,Tai L X. 2010. The shear-wave splitting study of Sichuan Zipingpu reservoir region[J]. Chinese Journal of Geophysics,53(9):2091–2101 (in Chinese).
朱瑜,桑琴,吴昌龙,杨通水. 2010. 威远气田震旦系灯影组储层特征研究[J]. 重庆科技学院学报(自然科学版),12(5):11–13. Zhu Y,Sang Q,Wu C L,Yang T S. 2010. Reservoir characteristics of the Sinian Dengying formation in the Weiyuan gasfield,Sichuan[J]. Journal of Chongqing University of Science and Technology (Natural Sciences Edition),12(5):11–13 (in Chinese).
Crampin S. 1981. A review of wave motion in anisotropic and cracked elastic-media[J]. Wave Motion,3(4):343–391. doi: 10.1016/0165-2125(81)90026-3
Crampin S,Peacock S. 2005. A review of shear-wave splitting in the compliant crack-critical anisotropic earth[J]. Wave Motion,41(1):59–77. doi: 10.1016/j.wavemoti.2004.05.006
Gao Y,Crampin S. 2004. Observations of stress relaxation before earthquakes[J]. Geophys J Int,157(2):578–582. doi: 10.1111/j.1365-246X.2004.02207.x
Gao Y,Crampin S. 2008. Shear wave splitting and earthquake forecasting[J]. Terra Nova,20(6):440–448. doi: 10.1111/j.1365-3121.2008.00836.x
Gao Y,Wu J,Fukao Y,Shi Y T,Zhu A L. 2011. Shear wave splitting in the crust in North China:Stress,faults and tectonic implications[J]. Geophys J Int,187(2):642–654. doi: 10.1111/j.1365-246X.2011.05200.x
Han C R,Huang Z C,Xu M J,Wang L S,Mi N,Yu D Y,Li H. 2019. Focal mechanism and stress field in the northeastern Tibetan Plateau:Insight into layered crustal deformations[J]. Geophys J Int,218(3):2066–2078. doi: 10.1093/gji/ggz267
Lei X L,Wang Z W,Su J R. 2019a. The December 2018 ML5.7 and January 2019 ML5.3 earthquakes in south Sichuan basin induced by shale gas hydraulic fracturing[J]. Seismol Res Lett,90(3):1099–1110. doi: 10.1785/0220190029
Lei X L,Wang Z W,Su J R. 2019b. Possible link between long-term and short-term water injections and earthquakes in salt mine and shale gas site in Changning,south Sichuan basin,China[J]. Earth Planet Phys,3(6):510–525. doi: 10.26464/epp2019052
Vinnik L P,Kind R,Makeyeva L I,Kosarev G L. 1989. Azimuthal anisotropy in the lithosphere from observations of long-period S-waves[J]. Geophys J Int,99(3):549–559. doi: 10.1111/j.1365-246X.1989.tb02039.x
Yang H F,Zhou P C,Fang N,Zhu G H,Xu W B,Su J R,Meng F B,Chu R S. 2020. A shallow shock:The 25 February 2019 ML4.9 earthquake in the Weiyuan shale gas field in Sichuan,China[J]. Seismol Res Lett,91(6):3182–3194. doi: 10.1785/0220200202
Zheng Y,Ge C,Xie Z J,Yang Y J,Xiong X,Hsu H T. 2013. Crustal and upper mantle structure and the deep seismogenic environment in the source regions of the Lushan earthquake and the Wenchuan earthquake[J]. Science China Earth Sciences,56(7):1158–1168. doi: 10.1007/s11430-013-4641-2
-
期刊类型引用(9)
1. 石玉燕,张春鹏,颜启,苗庆杰. 2022-01-15汤加火山喷发的地震学震级测定. 大地测量与地球动力学. 2023(02): 209-214 . 百度学术
2. 韩光洁,刘奕君,席楠. 2023年山东平原M_S 5.5地震宽频带面波震级和近场地震动反应谱空间分布特征分析. 地震地磁观测与研究. 2023(06): 13-19 . 百度学术
3. Qicheng Li,Jingwen Sun,Guimei Xi,Jing Liu. The Doppler effect induced by earthquakes: A case study of the Wenchuan M_S8.0 earthquake. Geodesy and Geodynamics. 2022(05): 435-444 . 必应学术
4. 孙冬军,刘芳,毕波. 芦山地震和九寨沟地震台站实测震级差异性分析. 中国地震. 2022(01): 112-119 . 百度学术
5. 石玉燕,张春鹏,颜启,苗庆杰. 2022年1月8日门源6.9级地震山东地震台网测定面波震级对比分析. 地震地磁观测与研究. 2022(04): 1-7 . 百度学术
6. 孔韩东,刘瑞丰,边银菊,李赞,王子博,胡岩松. 地震辐射能量测定方法研究及其在汶川8.0级地震中的应用. 地球物理学报. 2022(12): 4775-4788 . 百度学术
7. 陈建龙. 爆破地震波传播过程中的多普勒效应研究. 爆破器材. 2021(01): 48-52 . 百度学术
8. 李启成,闵也,何书耕. 用小波变换讨论地震中的多普勒现象. 吉林大学学报(地球科学版). 2021(03): 909-918 . 百度学术
9. 邓文泽. 2021年3月5日新西兰克马德克群岛M7.8面波震级测定差异分析. 地震科学进展. 2021(09): 395-399 . 百度学术
其他类型引用(4)