Shallow S-wave velocity structure in Chongqing area
-
摘要: 基于重庆台网及邻区固定台站自2011年1月至2018年10月观测到的远震波形记录,利用接收函数直达P波振幅约束浅层结构的方法计算了重庆及周边地区台站下方浅层S波速度结构,结果表明:重庆浅层S波速度结构与盆山构造明显相关,盆地内表现为低速异常,与沉积层特征相对应,大巴山、大娄山区域则表现为相对高速的特征;华蓥山断裂S波速度高于断裂两侧的沉积层速度,表现出四川盆地沉积层中间薄两边厚的特点;华蓥山以东的川东地区多为套滑脱层构造,其浅层S波速度表现为低速异常。最后讨论了2010年以来重庆5次显著地震的孕震环境:川东滑脱构造地区的垫江MS4.4地震和石柱MS4.5地震与该区明显的低速异常有关;荣昌MS4.7和MS4.8地震的震源区无明显的高低速特征,可能与注水相关;武隆MS5.0地震发生在高低速交界部位的有利于积累应变的高速体一侧。Abstract: Based on the teleseismic waveform data recorded at permanent stations in Chongqing and its adjacent area from January 2011 to October 2018, the shallow S-wave velocity structure beneath the stations of Chongqing and its surrounding area is calculated by the method to constrain shallow structure based on direct P-wave amplitude in receiver functions. The results indicated that, the velocity structure of shallow S-wave in Chongqing is obviously related to the basin mountain structure, the low velocity anomaly in the basin corresponds to the characteristics of sedimentary layer, and the Daba mountain and Dalou mountain regions show relatively high velocity anomalies. The S-wave velocity of Huayingshan fault is higher than that of sedimentary layers on both sides of the fault, implying that the sedimentary layers in Sichuan basin are thin in the middle and thick on both sides. Finally, we discussed the seismogenic environment of five earthquakes in Chongqing since 2010. The Dianjiang MS4.4 earthquake and the Shizhu MS4.5 earthquake located in the east Sichuan detachment structure area are related to the obvious low velocity anomaly. The Rongchang MS4.7 and MS4.8 earthquakes have no obvious characteristics of high or low velocities, which may be related to water injection. The Wulong MS5.0 earthquake occurred on the high velocity body side that is conducive to strain accumulation of the high and low velocity junction.
-
引言
重庆地区位于四川盆地东南缘,华南地震区的长江中游地震带内(图1)。重庆区域内的断裂主要以基底断裂为主,如华蓥山断裂、长寿—遵义断裂、七曜山—金佛山断裂、彭水断裂和城口断裂,这些基底断裂对重庆的地震活动具有明显的控制作用。此外,四川盆地与其周缘的龙门山、米仓山—大巴山、七曜山—金佛山、大娄山及大凉山形成了典型的盆山体系。这些盆山构造非常复杂,并且对盆地内油气的分布有较大的控制作用。四川盆地油气资源丰富,其东南缘的地震活动可能与所在地区的盐岩和油气的存在相关,几个典型的注水诱发地震序列的震源集中在三叠系以下的沉积岩地层中(雷兴林等,2014)。特殊的地质构造和地理位置造成了重庆地区地震强度不高、震源浅的特点(丁仁杰,李克昌,2004)。2010年以来重庆及周边地区共发生6次MS4.5以上地震,分别为2010年重庆荣昌MS4.7地震、2013年重庆石柱MS4.5地震、2013年湖北巴东MS5.1地震、2014年湖北秭归MS4.7地震、2016年重庆荣昌MS4.8地震和2017年重庆武隆MS5.0地震。这些地震的震源深度均较浅,其中荣昌MS4.7地震震源深度仅1.85 km左右(王小龙等,2012),最深的武隆MS5.0地震震源深度也仅10 km左右(李翠平等,2019)。
地震的孕育与介质性质及结构密切相关,重庆地区的地震震源较浅,因此有必要深入研究重庆地区的浅层速度结构。并且地壳的浅层速度结构可为盆地基底、微小震定位和地壳上地幔深部结构成像提供重要的研究基础,亦可为地震工程场地评估和防震减灾提供重要参考。
目前,地壳浅层S波速度结构研究主要有地震勘探和天然源等方法。地震勘探虽然具有信噪比高、精度高的特点,但是探测成本较高、对环境具有一定的破坏性。天然源则是利用背景噪声和天然地震事件等数据研究浅层结构,该方法利用地震台站的地震波形资料,节约了成本。例如,基于背景噪声瑞雷面波频散曲线研究地表至浅层数百米深度的S波速度结构(张宝龙等,2016;王娟娟等,2018)。
此外,接收函数直达P波垂向与径向的幅度比与浅层S波速度有关(Ammon,1991)。近年来,基于接收函数直达P波振幅约束浅层结构的新方法取得了新进展。Julià (2007)经过半空间接收函数公式推导,得出了接收函数直达P波的理论振幅公式,证明直达P波垂向与径向的幅度比与射线参数和浅层S波速度有关;并认为随着接收函数有限频率的增加,入射波的波长变短,S波速度值对应的深度范围变浅。钱银苹等(2018)基于接收函数直达P波理论振幅公式进行了数值实验,结果表明该公式能够对近地表S波速度结构进行有效约束,之后又基于不同速度模型对近地表S波速度对应的深度进行了数值试验,最后将该方法应用到青藏高原东北缘。王旭等(2019)同样针对该公式进行了理论模拟,并利用不同频率的接收函数直达P波振幅计算了青藏高原东北缘流动地震台阵下方的浅层S波速度结构。
本文将采用接收函数直达P波约束浅层S波速度的方法研究重庆及周边地区台站下方浅层S波速度结构,探讨该区域浅层速度结构与地质构造的关系,为研究区孕震环境提供科学依据。
1. 数据与方法
1.1 数据资料获取
本文收集整理了重庆台网及邻区固定台站自2011年1月至2018年10月观测到的远震波形记录,选取了MS≥5.5、震中距为30°—90°的远震波形(图2)。对所选波形作如下处理:① 截取P波到时前20 s和后80 s的远震波形数据;② 将截取后的N,E,Z分量波形旋转至R,T,Z分量;③ 选取多个高斯滤波因子α (共6个滤波因子:1.0,1.5,2.5,3,4,5,分别对应不同的通带频率,也即不同深度信息),利用时间域反褶积方法,计算P波接收函数。经过自动及人工挑选,平均每个滤波因子可挑选出震相清晰的接收函数1万2 500多条。图3为华蓥山断裂附近的ROC台6个滤波因子的R向接收函数。
1.2 台站下方浅层S波速度及误差分析
在P波初动附近(±2 s)提取径向接收函数和垂向接收函数的波峰极大值,用垂向分量最大振幅值对径向最大振幅值作归一处理(Ammon,1990),得到的结果即为接收函数直达P波绝对幅度的观测值。
根据入射P波引起的自由表面位移方程(Aki,Richards,2002)可计算直达P波绝对幅度的理论值,即
$$ \frac{{r}_{0}}{{{\textit{z}}}_{0}}=\frac{2p{\eta }_{{v}_{{\rm{S}}0}}}{{\dfrac{1}{v^2_{{\rm{S}}0}}}-2{p}^{2}} \text{,} $$ (1) 式中,
${\eta }_{{v}_{{\rm{S0}}}}=\sqrt {v^2_{\rm{S0}} - 2p^2}$ 为S波的垂直慢度,p为水平慢度(射线参数),$ {{v}}_{\rm{S}0} $ 为浅层S波速度。式(1)表示在低频时震相振幅与平均速度之间的关系(钱银苹等,2018)。绝对幅度的理论值与观测值之间的差异用Misfit (
$ {{v}}_{{\rm{S}}0} $ )表示,即$$ \mathrm{M}\mathrm{i}\mathrm{s}\mathrm{f}\mathrm{i}\mathrm{t} ( {v}_{{\rm{S}}0} ) =\sqrt{{\sum\limits _{i}{ ( {A}^{{\rm{obs}}}-{A}_{i}^{{\rm{syn}}} ) }^{2}}}\;\;\text{,} $$ (2) 式中,
$ {A}^{\mathrm{o}\mathrm{b}\mathrm{s}} $ 为绝对幅度的观测值,$ {A}_{i}^{\mathrm{s}\mathrm{y}\mathrm{n}} $ 为不同$ {{v}}_{{\rm{S}}0} $ 计算的理论值。Misfit最小值所对应的${{v}}_{\rm{S}0}$ 即为vS0的最优值,每条接收函数可确定一个$ {{v}}_{{{\rm{S}}}0} $ 。实际处理中,对于一个台站,按1°的震中距叠加接收函数,利用叠加的接收函数确定
$ {{v}}_{\rm{S}0} $ ,它们的平均值即为这个台站的浅层S波速度。因每叠加1°,接收函数都可以确定一个vS值,为了评估计算结果的质量,以均方误差ε来表示计算误差。当误差大于0.4时,舍弃该台站的结果。图4为ROC台滤波因子为5时确定的
$ {{v}}_{{\rm{S}}0} $ 值,计算结果为(2.68±0.13) km/s。2. 结果
根据上述流程,对重庆及邻区固定台站的不同高斯滤波因子(1.0,1.5,2.5,3,4,5)的接收函数进行了处理,共得到6个不同频率接收函数的浅层S波速度。本文参考了钱银苹等(2018)利用数值试验估计出的不同频率S波速度对应的深度范围,其中:α为1时表示10—11 km深度;α为1.5时表示7 km左右深度;α为2.5时表示4—5 km深度;α为3和4时表示3—4 km深度;α为5时表示2—3 km深度。根据不同滤波因子的接收函数计算的重庆及邻区台站的vS,绘制了S波速度分布图(图5)。从图中可以看出,重庆境内长江沿线的近地表S波速度在2—3 km左右深度处较低(图5a)。七曜山—金佛山断裂为四川盆地的东南边界,在2—3 km深度处四川盆地的S波低速异常与四川盆地沉积特征相符。华蓥山断裂附近台站(ROC,HYS,XHA)的S波速度结构与断裂两侧的台站下方相比表现为高速异常,与王小龙等(2013)周期5 s的相速度结果一致。城口断裂以北的大巴山地区(CHK,HCB,WUX,SHL,SXH)S波速度较高。彭水台网(GOT,LUJ,LAX,XIT,ANZ)S波速度较低,应为水库渗透影响,导致该地区速度出现相对低值。而水库西北方的武隆地区呈现明显高速,速度为2.9 km/s。α=2.5—4表示3—5 km左右深度,该深度范围与α=5对应深度的速度分布特征较一致,台站的速度变化小于0.2 km/s,并且速度值随深度增大。α=1.5和1对应的深度范围,研究区的S波速度进一步增大。武隆、仙女山、涪陵地区高值区范围增大,该区S波平均速度约为3.33 km/s,而四川盆地在α=1.5和1时对应的深度范围仍为明显的低速异常,说明8—10 km左右深度,四川盆地仍为沉积层。
3. 讨论
3.1 重庆及周边浅层S波速度结构
除彭水台网外,重庆浅层S波速度结构与盆山构造明显相关,盆地内低速异常与沉积层特征对应,大巴山、大娄山区域则表现相对高速特征。王小龙等(2013)及王小龙(2017)利用地震背景噪声成像技术反演三峡库区5—40 s的瑞雷面波相速度分布模型及三维剪切波速度模型,结果也显示四川盆地周边的龙门山、米仓山、大巴山、七曜山和大娄山等地质构造与地壳S波速度结构具有明显的对应关系。α=1.5和1对应的深度范围,四川盆地仍为明显的低速异常,说明8—10 km左右深度,四川盆地仍为沉积层;宋晓东等(2015)通过对面波层析成像的三维S波速度模型的研究认为,四川盆地内部上地壳绝大部分都表现为低速,反映出盆地内存在很厚的沉积层;李建有等(2018)利用低频和高频接收函数的两步反演地壳S波速度结构,认为四川盆地4—10 km呈现相对低速区,S波速度介于2.8—3.2 km/s之间,基本反映了四川盆地绝大部分为沉积岩或新生代松软覆盖层的低速分布特征(范军等,2015);华蓥山断裂S波速度(HYS,ROC,XHA)均较其两侧沉积层速度高,反映了四川盆地沉积层中间薄两侧厚的特点;四川盆地的深地震反射剖面显示四川盆地上地壳以沉积层为主,从西北向东南逐渐变薄,在龙门山前沉积层厚度超过15 km,在华蓥山下沉积层减薄到8 km并发生了褶皱变形,形成了华蓥山薄皮褶皱冲断带(王海燕等,2017),与本文计算的S波特征相符。华蓥山以东的川东地区(XIM,FUP,LIJ,HOT,WAZ,SHG,LOH,SHZ,YAL,WUX,JIP)的浅层S波速度表现为低速异常。川东地区是华蓥山与七曜山之间的褶皱变形区,浅层发育一系列NE-NNE向的隔档式褶皱,并且该区发育以寒武系和三叠系膏盐岩为主的多套滑脱层,挤压变形后形成多套滑脱层构造,滑脱构造下覆寒武系底达8 km左右(王志勇等,2018)。
3.2 孕震环境
本文讨论2010年以来重庆5次显著地震,分别为2010年重庆荣昌MS4.7地震(图5a)、2013年重庆石柱MS4.5地震(图5c)、2016年垫江MS4.4地震(图5a)、2016年重庆荣昌MS4.8地震(图5b)和2017年重庆武隆MS5.0地震(图5f)。5次地震的震中位置来自国家地震台网中心正式测定和精定位,震源深度来自CAP方法和精定位。
2010年9月10日荣昌MS4.7地震震源机制解最佳深度为2.93 km (王志伟等,2018),2016年12月27日MS4.8地震最佳拟合深度为3.8 km,同样为浅源地震。α=4和5对应深度处ROC台S波速度分别为2.69 km/s和2.68 km/s,与重庆地区平均速度2.66 km/s较一致,没有明显的高低速特征。许多研究表明荣昌地区地震活动与天然气开采过程中的注水密切相关,注水引起的流体压力扩散导致断层带孔隙压力增加,使得断层带强度降低,触发断层活动(王志伟等,2018),而震源深度与主要储藏层及注水井有关(王小龙等,2012)。
2016年垫江MS4.4地震的波形具有明显的浅源特征,CAP方法反演出的震源深度约为1—2 km (黄世源等,2020)。地震位于川东滑脱构造地区,黄世源等(2020)认为在区域应力作用下,震源区沉积盖层内部的软弱地层受挤压变形,导致地层内斜交的裂隙发生切层失稳滑动,从而诱发了垫江MS4.4地震,而震中附近(XIM台)较低速度的S波为该地震的发生提供了条件。
石柱MS4.5柱地震Rg面波发育明显,震源深度约为4 km (王同军等,2015),地震同样发生在低速异常的川东构造区内,与垫江MS4.4地震不同的是,石柱MS4.5地震震中附近有方斗山断裂。
武隆MS5.0地震发生在七曜山—金佛山断裂与彭水断裂之间,震源深度约为10 km (李翠平等,2019)(图5f),震源区附近S波速度显示出明显的横向不均匀分布特征,武隆地区S波速度明显高于川东滑脱构造区和彭水台网。其中α=1—1.5对应的深度范围,WUL台出现明显的高速异常,S波速度为3.22—3.44 km/s。武隆MS5.0地震发生在高低速交界部位的有利于积累应变的高速体一侧,许多研究表明地壳高低速过渡区域可能是中强地震孕育和发生的关键部位。
4. 结论
本文基于重庆台网及邻区固定台站自2011年1月至2018年10月观测到的远震波形记录,利用接收函数直达P波振幅约束浅层结构的方法计算了重庆及周边地区台站下方浅层S波速度结构,该方法不依赖初始模型,不需要反演,能有效抑制速度结果的不唯一性。最终形成以下几点认识:
重庆浅层S波速度结构与盆山构造明显相关,盆地内低速异常与沉积层特征对应,大巴山、大娄山区域均表现为相对高速特征;华蓥山断裂S波速度均较其两侧沉积层速度高,表现出四川盆地沉积层中间薄两侧厚的特点;华蓥山以东的川东地区的浅层S波速度表现为低速异常,该区存在多套滑脱层构造。
此外,本文还讨论了2010年以来重庆5次显著地震的孕震环境:川东滑脱构造地区的垫江MS4.4地震和石柱MS4.5地震两个浅源地震与该区明显的低速异常有关;荣昌MS4.7和MS4.8地震震源区没有明显的高低速特征,可能与天然气开采过程中的注水相关;武隆MS5.0地震发生在高低速交界部位的有利于积累应变的高速体的一侧。
-
-
丁仁杰, 李克昌. 2004. 重庆地震研究[M]. 北京: 地震出版社: 107. Ding R J, Li K C. 2004. Earthquake Research of Chongqing[M]. Beijing: Seismological Press: 107 (in Chinese).
范军,朱介寿,江晓涛,吴朋. 2015. 利用接收函数方法研究四川地区地壳结构[J]. 地震,35(1):65–76. doi: 10.3969/j.issn.1000-3274.2015.01.008 Fan J,Zhu J S,Jiang X T,Wu P. 2015. Crustal structure of Sichuan from receiver functions[J]. Earthquake,35(1):65–76 (in Chinese).
黄世源,李翠平,王赞军,高见,郭卫英. 2020. 重庆垫江MS4.4地震:一个可能发生在浅地表的地震[J]. 地震学报,42(1):34–43. doi: 10.11939/jass.20190094 Huang S Y,Li C P,Wang Z J,Gao J,Guo W Y. 2020. Dianjiang,Chongqing,MS4.4 earthquake:A possible shallow surface earthquake[J]. Acta Seismologica Sinica,42(1):34–43 (in Chinese).
雷兴林,李霞颖,李琦,马胜利,付碧宏,崔银祥. 2014. 沉积岩储藏系统小断层在油气田注水诱发地震中的作用:以四川盆地为例[J]. 地震地质,36(3):625–643. doi: 10.3969/j.issn.0253-4967.2014.03.007 Lei X L,Li X Y,Li Q,Ma S L,Fu B H,Cui Y X. 2014. Role of immature faults in injection-induced seismicity in oil/gas reservoirs:A case study of the Sichuan basin,China[J]. Seismology and Geology,36(3):625–643 (in Chinese).
李翠平,唐茂云,郭卫英,黄世源,王小龙,高见. 2019. 2017年11月23日重庆武隆MS5.0地震序列重定位及发震断层分析[J]. 地震地质,41(3):603–618. doi: 10.3969/j.issn.0253-4967.2019.03.005 Li C P,Tang M Y,Guo W Y,Huang S Y,Wang X L,Gao J. 2019. Relocation of the 23 November 2017 Wulong MS5.0 earthquake sequence and analysis of its seismogenic fault[J]. Seismology and Geology,41(3):603–618 (in Chinese).
李建有,石宝文,徐晓雅,胡家富. 2018. 利用远震接收函数探测四川盆地及周边地区的地壳结构[J]. 地球物理学报,61(7):2719–2735. doi: 10.6038/cjg2018L0652 Li J Y,Shi B W,Xu X Y,Hu J F. 2018. Crustal structure beneath the Sichuan basin and adjacent regions revealed by teleseismic receiver functions[J]. Chinese Journal of Geophysics,61(7):2719–2735 (in Chinese).
钱银苹,沈旭章,李翠芹,梅秀苹,焦煜媛. 2018. 利用接收函数方法确定青藏高原东北缘近地表S波速度[J]. 地球物理学报,61(10):3951–3963. Qian Y P,Shen X Z,Li C Q,Mei X P,Jiao Y Y. 2018. Constraining the sub-surface S-wave velocity of the northeastern margin of Tibetan Plateau with receiver functions[J]. Chinese Journal of Geophysics,61(10):3951–3963.
宋晓东,李江涛,鲍学伟,李思田,王良书,任建业. 2015. 中国西部大型盆地的深部结构及对盆地形成和演化的意义[J]. 地学前缘,22(1):126–136. doi: 10.13745/j.esf.2015.01.011 Song X D,Li J T,Bao X W,Li S T,Wang L S,Ren J Y. 2015. Deep structure of major basins in western China and implications for basin formation and evolution[J]. Earth Science Frontiers,22(1):126–136 (in Chinese).
王海燕,高锐,卢占武,李文辉,郭华,熊小松,酆少英,李洪强,赵玉莲. 2017. 四川盆地深部地壳结构:深地震反射剖面探测[J]. 地球物理学报,60(8):2913–2923. doi: 10.6038/cjg20170801 Wang H Y,Gao R,Lu Z W,Li W H,Guo H,Xiong X S,Feng S Y,Li H Q,Zhao Y L. 2017. Deep crustal structure in Sichuan basin:Deep seismic reflection profiling[J]. Chinese Journal of Geophysics,60(8):2913–2923 (in Chinese).
王娟娟,姚华建,王伟涛,王宝善,李成,魏斌,冯磊. 2018. 基于背景噪声成像方法的新疆呼图壁储气库地区近地表速度结构研究[J]. 地球物理学报,61(11):4436–4447. doi: 10.6038/cjg2018M0025 Wang J J,Yao H J,Wang W T,Wang B S,Li C,Wei B,Feng L. 2018. Study of the near-surface velocity structure of the Hutubi gas storage area in Xinjiang from ambient noise tomography[J]. Chinese Journal of Geophysics,61(11):4436–4447 (in Chinese).
王同军,蒋璀,贺曼秋,王洪超,陈凯. 2015. 2013年7月18日重庆石柱M4.5地震震源深度探讨[J]. 地震地磁观测与研究,36(3):27–31. doi: 10.3969/j.issn.1003-3246.2015.03.005 Wang T J,Jiang C,He M Q,Wang H C,Chen K. 2015. About focal depth of Shizhu M4.5 earthquake on 18 July,2013[J]. Seismological and Geomagnetic Observation and Research,36(3):27–31 (in Chinese).
王旭,陈凌,凌媛,高一帆,张建勇,姚华建. 2019. 基于接收函数直达P波振幅研究地壳浅层S波速度结构新方法及在青藏高原东北缘的应用[J]. 中国科学:地球科学,49(11):1788–1800. Wang X,Chen L,Ling Y,Gao Y F,Zhang J Y,Yao H J. 2019. A new method to constrain shallow crustal S-wave velocities based on direct P-wave amplitudes in receiver functions and its application in northeastern Tibet[J]. Science China Earth Sciences,62(11):1819–1831. doi: 10.1007/s11430-018-9443-6
王小龙,马胜利,雷兴林,郭欣,王强,余国政,勾宪斌,桑原保人,今西和俊,蒋霞东. 2012. 重庆荣昌诱发地震区精细速度结构及2010年ML5.1地震序列精确定位[J]. 地震地质,34(2):348–358. doi: 10.3969/j.issn.0253-4967.2012.02.013 Wang X L,Ma S L,Lei X L,Guo X,Wang Q,Yu G Z,Gou X B,Kuwahara Y,Imanishi K,Jiang X D. 2012. Fine velocity structure and relocation of the 2010 ML5.1 earthquake sequence in Rongchang gas field[J]. Seismology and Geology,34(2):348–358 (in Chinese).
王小龙,马胜利,郭志,雷兴林,夏英杰,郭欣,余国政,勾宪斌,蒋霞东. 2013. 利用地震背景噪声成像技术反演三峡库区及邻近地区地壳剪切波速度结构[J]. 地球物理学报,56(12):4113–4124. doi: 10.6038/cjg20131216 Wang X L,Ma S L,Guo Z,Lei X L,Xia Y J,Guo X,Yu G Z,Gou X B,Jiang X D. 2013. S-wave velocity of the crust in Three Gorges Reservoir and the adjacent region inverted from seismic ambient noise tomography[J]. Chinese Journal of Geophysics,56(12):4113–4124 (in Chinese).
王小龙. 2017. 利用地震波研究三峡重庆地区及邻区地壳速度结构[M]. 北京: 地震出版社: 53–54. Wang X L. 2017. Study on Crustal Velocity Structure in Chongqing Area of the Three Gorges and its Adjacent Areas by Using Seismic Waves[M]. Beijing: Seismological Press: 53–54 (in Chinese).
王志伟,王小龙,马胜利,雷兴林,刘培洵. 2018. 重庆荣昌地区注水诱发地震的时空分布特征[J]. 地震地质,40(3):523–538. doi: 10.3969/j.issn.0253-4967.2018.03.002 Wang Z W,Wang X L,Ma S L,Lei X L,Liu P X. 2018. Detailed temporal-spatial distribution of induced earthquakes by water injection in Rongchang,Chongqing[J]. Seismology and Geology,40(3):523–538 (in Chinese).
王志勇,康南昌,李明杰,臧殿光,张怡,杨昕. 2018. 四川盆地川东地区滑脱构造特征[J]. 石油地球物理勘探,53(增刊):276–286. doi: 10.13810/j.cnki.issn.1000-7210.2018.S1.044 Wang Z Y,Kang N C,Li M J,Zang D G,Zhang Y,Yang X. 2018. Detached structure features in East Sichuan basin[J]. Oil Geophysical Prospecting,53(S1):276–286 (in Chinese).
张宝龙,李志伟,包丰,邓阳,游庆瑜,张森琦. 2016. 基于微动方法研究五大连池火山区尾山火山锥浅层剪切波速度结构[J]. 地球物理学报,59(10):3662–3673. doi: 10.6038/cjg20161013 Zhang B L,Li Z W,Bao F,Deng Y,You Q Y,Zhang S Q. 2016. Shallow shear-wave velocity structures under the Weishan volcanic cone in Wudalianchi volcano field by microtremor survey[J]. Chinese Journal of Geophysics,59(10):3662–3673 (in Chinese).
Aki K, Richards P G. 2002. Quantitative Seismology[M]. 2nd ed. Sausalito CA, USA: University Science Books: 146-202.
Ammon C J,Randall G E,Zandt G. 1990. On the nonuniqueness of receiver function inversions[J]. J Geophys Res:Solid Earth,95(B10):15303–15318. doi: 10.1029/JB095iB10p15303
Ammon C J. 1991. The isolation of receiver effects from teleseismic P waveforms[J]. Bull Seismol Soc Am,81(6):2504–2510. doi: 10.1785/BSSA0810062504
Julià J. 2007. Constraining velocity and density contrasts across the crust-mantle boundary with receiver function amplitudes[J]. Geophys J Int,171(1):286–301. doi: 10.1111/j.1365-2966.2007.03502.x