利用虚震源反射成像法研究江南造山带东段上地壳结构

肖尧, 唐启家, 王林松

肖尧, 唐启家, 王林松. 2020: 利用虚震源反射成像法研究江南造山带东段上地壳结构. 地震学报, 42(5): 567-579. DOI: 10.11939/jass.20190193
引用本文: 肖尧, 唐启家, 王林松. 2020: 利用虚震源反射成像法研究江南造山带东段上地壳结构. 地震学报, 42(5): 567-579. DOI: 10.11939/jass.20190193
Xiao Yao, Tang Qijia, Wang Linsong. 2020: Crustal structure across the Jiangnan orogen from teleseismic virtual-source reflection method. Acta Seismologica Sinica, 42(5): 567-579. DOI: 10.11939/jass.20190193
Citation: Xiao Yao, Tang Qijia, Wang Linsong. 2020: Crustal structure across the Jiangnan orogen from teleseismic virtual-source reflection method. Acta Seismologica Sinica, 42(5): 567-579. DOI: 10.11939/jass.20190193

利用虚震源反射成像法研究江南造山带东段上地壳结构

基金项目: 中央高校基本科研业务费专项(CUGCJ1707)和中国地质大学(武汉)地质过程与矿产资源国家重点实验室专项(MSFGPMR01-5)联合资助
详细信息
    通讯作者:

    唐启家: e-mail:iori89724@hotmail.com

  • 中图分类号: P315.2

Crustal structure across the Jiangnan orogen from teleseismic virtual-source reflection method

  • 摘要: 江南造山带是华南地区扬子地块与华夏地块碰撞的产物,其地壳构造记录了两地块的碰撞过程,研究江南造山带的地壳构造有助于重建扬子地块与华夏地块的碰撞过程。本研究在江南造山带上布设了两条流动地震台阵,利用虚震源反射法提取其所记录的远震事件初至P波在地表的反射波(PPdp)波形,重构了研究区内两条测线下方的上地壳结构。结果显示:江绍断裂两侧上地壳沉积层的厚度变化明显,推断该断裂是扬子地块与华夏地块的东边界;相较华夏地块,江南造山带与扬子地块的层位连续性更强,符合江南造山带先与扬子地块合为整体后再与华夏地块碰撞的多期构造过程及其对应产生的亲扬子地块属性;江绍断裂西北侧的地层不整合以及赣东北断裂区域的断陷构造,可为了解古华夏洋向扬子陆块俯冲及碰撞和随后的构造运动过程提供参考依据。
    Abstract: The Jiangnan orogen belt (JOB) was formed by the collision between the Yangtze block (YB) and the Cathaysia block (CB). The crustal structure of the JOB recorded the collisional process between the two blocks. Hence, studying crust of the JOB can illuminate the process between the YB and the CB, and thus better understand the tectonic evolution of the South China Plate (SCP). In this study, we applied the teleseismic virtual reflection method to derive the phase of the reflection waves (PPdp) of large teleseismic earthquakes recorded by two temporary seismic arrays that were deployed around the JOB and obtained the upper crust structures. The results suggest that the thickness of the upper crust sediments on two sides of the Jiangshao fault changes obviously. It is also inferred that the Jiangshao fault is the eastern boundary between the Yangtze block and the Cathaysia block. The JOB and the YB have more apparent stratigraphic continuity than that of the CB. This result agrees with a multi-stage tectonic model which suggests that the JOB first bounded the YB and then collided with the CB. Hence the characteristics of the JOB are closer to that of the YB. The stratigraphic unconformities in the northwest of the Jiangshao fault and the fault structure in the northeast Jiangxi region provides a basic profile for understanding the subducting and colliding processes between the ancient Cathaysian ocean and the YB as well.
  • 井水位观测是地震地下流体前兆监测的重要手段之一(刘耀炜,2006),井水位的异常变化能够客观地、灵敏地反映地壳介质应力-应变变化。在地震孕育、发生的过程中,局部应力的加卸载作用或断层的大尺度活动均会导致岩体发生变形,使含水层孔隙压力产生变化,从而导致地壳介质渗流场的变化,井孔水位也随之动态变化(田竹君,谷圆珠,1985Montgomery,Manga,2003)。付虹等(1997)Sato等(2004)车用太等(2008)的观测实践已表明,地震前井水位是存在异常变化的。研究地震孕育、发生过程中井水位的各种异常变化,对揭示地震孕育、发生和发展过程具有重要的物理意义,在地震前兆监测与地震预测实践中发挥着重要作用(车用太等,2003)。

    云南省会泽井的井水位自2012年数字化观测以来,该井周边及附近地区发生的中强地震震前井水位均出现了显著的上升,显示出较强的映震能力。本文拟以会泽井水位为研究对象,系统地分析会泽井水位在几次中强地震前出现的显著异常特征,并依据会泽井水位及其周边的地震活动情况,利用Molchan图表法对该井水位观测资料进行预测效能定量化评估与检验,提取其具有前兆指示意义的地震预测预报指标,进而对会泽井水位上升异常的机理进行解释,以期为本地区的震情跟踪及地震形势判定提供可能的判定依据。

    图  1  会泽井的水文地质条件(a)和井孔柱状图(b)
    Figure  1.  Hydrogeological conditions (a) and the bore hole histogram (b) of Huize well

    会泽井位于云南省曲靖市会泽县城以西18 km处的娜姑湖盆地边缘,地理坐标为(26.52°N,103.15°E),海拔为2 005 m,所属区域属于云贵高原中部盆岭地貌,地势北高南低;位于川滇菱形地块构造的东南部,小江断裂东侧8 km处,白雾街不对称双弧形构造的交会部位。会泽井深103.15 m,套管深度为87.80 m,水位埋深约30.00 m,其中34.06—87.80 m为滤水管,观测段为34.06—103.15 m。会泽井的水文地质条件如图1a所示,地下水为裂隙承压水,水温为16.0℃左右,受区域活动断裂控制和影响,该井的松散堆积层孔隙水的富水性弱,单井涌水量小于100 m3/d,地表以下11.19 m深处为0.15 m厚的棕色黏土隔水层,其隔水效果较好。井孔周边岩性如图1b所示:地表为第四系湖相沉积,观测含水层岩性为第四系风化玄武岩、砾石、碎石、黏土含细砂层;下伏地层为二叠系灰色玄武岩。

    会泽井水位采用LN-3型数字化水位仪,探头放置于水下4.925 m,自2012年采用数字化观测以来,该井水位的动态特征表现为持续下降趋势,无年变,不受降雨影响,下降速率为0.3 m/a。如图2a所示,在表1所示的4次地震前,该井水位在多年缓降趋势的正常动态背景上,均表现出显著的突升异常现象,异常幅度较大,表明其具有一定的映震能力。

    图  2  会泽井静水位与该地区的降雨观测曲线(a)及4次地震前会泽井静水位的显著异常变化(b)
    Figure  2.  Curves of water level in Huizi well and the rainfall in the area (a) as well as abnormal changes of water level of Huize well before four earthquakes (b)
    表  1  会泽井水位出现异常变化及对应的地震
    Table  1.  Abnormal changes of water level and corresponding earthquakes
    发震日期
    年-月-日
    震中位置 地点 MS 井震距
    /km
    水位变化
    幅值/m
    水位变化开始时间
    年-月-日
    水位变化距
    发震时间/d
    北纬/° 东经/°
    2012−09−07 27.50 103.95 云南彝良 5.7 136 0.27 2012−08−13 25
    2013−04−20 30.30 103.00 四川芦山 7.0 421 0.21 2013−04−06 14
    2014−05−07 25.48 101.92 云南元谋 4.7 161 0.38 2014−04−18 19
    2014−08−03 27.10 103.40 云南鲁甸 6.5 71 0.33 2014−07−20 14
    下载: 导出CSV 
    | 显示表格

    为进一步分析中强地震前会泽井水位的异常变化特征,分别将图2a中的4次地震前井水位出现的显著异常时段经时间平移后,叠加异常形态曲线,得到各异常动态曲线,如图2b所示,可见4次地震前会泽井水位的异常形态特征较为相似,均为短时间内的上升异常变化,变化幅值存在差异,介于0.21—0.38 m之间,详见表1。其中井震距最短的地震为2014年8月3日鲁甸MS 6.5地震,震前14天井水位呈明显上升变化,幅值为0.33 m,并在地震发生后又出现了一次幅值为0.32 m的明显上升变化(此为震后效应)。几次地震前会泽井水位异常变化距离发震的时间均在1个月以内,说明该井水位具有显著的地震短临异常特征。

    从会泽井水位异常的变化特征来看,当井水位出现上升异常变化之后,在一定时间内均会发生中强地震,说明该井水位异常与地震活动存在一定的关系,据此可为该区域的地震形势判定提供可能的参考依据。

    会泽井水位出现上升异常与本区或周边地区中强以上地震发生的对应关系,显示出会泽井水位在多次中强震前都有较好的映震能力,但并非该区域内所有地震震前均存在井水位的异常变化。因此,本文将采用Molchan图表法对其预测效能进行定量化检验。Molchan图表法是目前国际上“地震可预测性合作研究”计划中采用的6种统计检验方法之一,该方法不仅能客观地和科学地进行地震预测评估,还能解决固定研究区域内地震时间的预测问题,并能给出相应的概率解释(Jordan,2006),已被广泛应用于常规的确定性和概率性预测的统计检验和效能评估中(Zechar,Jordan,2008)。

    研究区范围为(21°—32°N,97°—110°E),地震目录引自中国地震台网中心(2017),震源机制解引自USGS (2017)。依据震级与震中距一般对应关系:MS4.0—4.9地震的震中距小于200 km;MS5.0—5.9地震的震中距小于250 km;MS6.0—6.9地震的震中距小于300 km;MS7.0以上地震的震中距小于500 km,筛选出2012—2015年会泽井周围发生的地震,共31次,其中MS4.0—4.9地震19次,MS5.0—5.9地震8次,MS6.0—6.9地震3次,MS7.0—7.9地震1次。利用实际观测数据进行地震预测分析时,常需要对数据进行预处理。由于会泽井水位无年变特征,水位的异常特征表现为高值异常,对所积累的震例进行日常跟踪分析可知,地震一般发生在井水位异常升高之后,因此,本文选取去趋势处理和差分分析两种方法分别对会泽井水位的原始数据进行预处理,再利用Molchan图表法对水位异常期相应的地震依次进行检验,并对不同预处理方法的检验结果加以分析,最终得到会泽井的预测效能和优势对应异常时间的定量结果。

    Molchan图表法主要是针对预测值与目标地震差异度的检验(Molchan,1990),该方法既能直观地反映观测资料的整体预测效能,又能定量地分析异常,提取最佳阈值所对应的异常判定指标。图表法所计算的变量为:漏报率v,即预测无震而实际发震的数量与总的实际发震数之比;异常时空占有率τ,即以不同的阈值提取异常的时空范围与总的时空范围之比。

    本文利用此方法对会泽井水位观测数据的预测效果进行检验时,需考虑时间占有率,即所有异常时间及其有效预测期所累积的去除二者重复时段后的总时间长度除以被检验数据的总时间长度。例如:被检验数据的起始时间为2010年1月1日,结束时间为2010年12月31日,总时间长度为365天。在这一年中出现了测值高于某一阈值的异常3次,每次异常的持续时间为15天,预测期为30天,并且3次异常的持续时间和预测期中共计有20天是重复的。那么,其时间占有率为[3×(15+30)−20]/365 ≈ 0.315。通过不断降低预测的“警报”阈值,分别计算异常在时间上的占有率τ和相应的漏报率v,得到Molchan图表中的τ-v曲线。曲线与图表的边界线所包围的面积代表检验的预测效果,面积越小,预测效果越好;此外,还需参考概率增益G。在Molchan图表法中G被定义为(Molchan,1991Zechar,Jordan,2008

    $G {\text{=}} \frac{{1 - \nu }}{\tau }{\text{,}}$

    (1)

    G越大,预测效果越好,若τ-v曲线接近于直线G=1,表示无统计意义。实际检验的具体过程为给定阈值后,超出阈值的数据为异常值,地震发生在异常值所在时段及其有效预测期之外时,称为漏报。

    图  3  会泽井周边地震的空间分布(a)及水位的Molchan检验结果(b)
    Figure  3.  The spatial distribution of earthquakes around the Huize well (a) and the test results of water level by Molchan chart method (b)

    2012—2015年会泽井周围发生的31次地震的空间分布如图3a所示,利用Molchan检验法对会泽井水位去趋势后数据进行预测效能检验,阈值由大到小滑动,根据阈值确定时间占有率τ和漏投率v,得到τ-v曲线结果如图3b所示。图中点(1.0,0)表示地震全都报准,但其时间占有率也最大,占据整个数据的时间段;点(0,1.0)表示地震全部漏报,其时间占有率最小,相当于未作出预测。因此,评估观测资料的预测效能需要从漏报率和时间占有率这两方面综合判定。

    图  4  川滇菱形地块及附近地震的空间分布(a)及水位Molchan检验结果(b)
    Figure  4.  The spatial distribution of earthquakes in and around Sichuan-Yunnan block (a) and the test results of water level by Molchan chart method (b)

    对会泽井水位观测数据经去趋势处理后的Molchan检验结果显示,该井水位对MS≥4.0地震与MS≥5.0地震的预测效果基本相当,且G约为2,只是MS≥5.0地震的时间占有率略小,说明从去趋势曲线来看,会泽井水位对于其周围或附近地区、整个大区域不同的震级档未显示出明显的优势预测效果。

    会泽井位于川滇菱形地块的东南部,处于多条断裂的交会部位,其所在地区的地震活动受区域活动断裂控制和影响,因此,针对MS≥5.0地震按会泽井所在川滇菱形地块以东和以西进行分选,地震震中分布如图4a所示;分别对川滇菱形地块以东和以西的地震预测效能进行Molchan检验,结果如图4b所示,可以明显地看出:会泽井水位去趋势曲线对地块以东的地震预测效能较好,检验线与横纵坐标包围的面积较小,大部分数据点G>2;而对地块以西的地震预测效能并不理想,基本都在随机线G=1附近,且G<2,说明当会泽井水位出现异常变化时,中强地震发生的优势预测区域为川滇菱形地块以东,这可以为将来地震发生的地点和强度提供一定的参考依据。

    图  5  会泽井水位5日差分曲线(a,c)及Molchan检验结果(b,d)
    Figure  5.  Five-day difference water level curves (a,c) and the test results of water level by Molchan chart method (b,d)

    以上分析和检验是基于对原始水位数据进行去趋势预处理后的数据,考虑到会泽井水位震前多表现为上升型高值异常,且不同的预处理方法可能会得到不同的结果,因此为了更加突出会泽井水位的高值异常,分析不同的预处理方法对检验结果的影响,本文还尝试用差分法对其数据进行预处理,结果如图5a图5c所示,可见绝大部分地震前都出现了井水位高值现象,但并不是所有高值之后均有地震发生。

    利用Molchan图表法对差分处理后的会泽井水位进行检验,结果如图5b图5d所示。由图5b可直观地看出,MS≥5.0地震的红色检验线包围的面积较MS≥4.0地震的绿色检验线包围面积更小,且MS≥5.0地震的概率增益更大,说明会泽井水位5日差分曲线对其附近或周围地区MS≥5.0地震的预测效能更好;同样,如图5d所示,川滇菱形地块以东地震的红色检验线所包围的面积较地块以西地震的绿色检验线所包围面积更小,其概率增益更大,说明会泽井水位5日差分曲线对川滇菱形地块以东地震的预测效能更好。由此可见,利用去趋势法和差分法对会泽井水位进行预处理后,二者的检验效果与已有认识基本一致,即会泽井水位的上升高值异常对周边及川滇菱形地块以东的中强地震具有一定的预测效能。

    图  6  会泽井水位去趋势(a,b)和5日差分值(c,d)处理的预测效能
    Figure  6.  Prediction efficiency of water level in Huize well based on the de-trended data (a,b) and five-day difference data (c,d)

    基于上述分析,分别对会泽井去趋势曲线和5日差分曲线进行预测效能的检验,结果如图6a图6c所示。可以看出,会泽井水位在去趋势和5日差分处理后的整体预测效能都很好,即概率增益均大于1,而且有效预测时间的预测效能均在0.5以上,尤其是5日差分处理后的会泽井水位绝大多数预测时间的预测效能均在0.6以上,部分在0.8以上,说明会泽井水位用5日差分方法进行预处理后,预测效能更好。

    图6b图6d分别为两种方法预测效能的优势地震对应时段,可见,会泽井水位去趋势处理后,优势地震对应时间段为3个月以内,预测效能介于0.6—0.8之间,3个月之后预测效能降低至同一水平(图6b);会泽井水位5日差分处理后,优势地震对应时间段也为3个月以内,预测效能介于0.7—0.8之间,3个月之后预测效能逐渐降低,直至同一水平(图6d),说明会泽井水位在1—3个月以内的短期预测效果较好。

    综上所述,从空间预测角度来看,会泽井水位的上升高值异常对周边及川滇菱形地块以东的中强地震具有一定的预测效能,从时间预测角度来看,其预测期多在1—3个月,且1个月左右的预测效能最佳。定量检验的结果与水位原始曲线(经去趋势处理)的映震效果基本一致,如表1中统计的4次显著震例,其异常出现时间距发震时间均在30天以内。

    国内外研究人员已经开展了大量的地下水位异常变化与地震关系的机理研究,提出了多种模式。20世纪70年代初的扩容-扩散模式、裂隙串通或雪崩模式和微裂-顶位移模式,均认为地震孕育到发生的过程中,震源附近应力积累使得地下介质变化,引起地震震中位置或周围地下水位的异常(国家地震局预测预防司,1997)。车用太和杨会年(1985)总结了国外研究关于震前地下水动态异常的3种机理:滑动机理、破裂机理和变形锋传播机理,归纳出震前断裂的预滑或微破裂的产生与发展可形成震中压缩区和拉张区,从而导致井水位上升或下降异常。车用太(1990)进一步研究认为震前地下水位异常与震源区域应力场活动相关。至20世纪80年代,地下水前兆机理研究从“震源效应”转向“大范围应力场效应”,认为井水位上升、下降异常与地壳应力变化—含水层变形—水动力条件改变等作用过程有关,是区域应力场和震源应力场共同作用的结果(车用太,杨会年,1985王吉易等,2002)。

    2012年以来会泽井水位在多次中强地震前均出现了上升高值异常变化,本文利用Molchan图表法的检验结果显示,会泽井水位的高值异常对川滇菱形地块以东地震的预测效能较高,这可能与会泽井所处的构造位置及井孔的水文条件密切相关。从构造区域环境来看,会泽井位于则木河—小江断裂和西鱼河—昭通断裂的交会区附近(图7a),两断裂位于川滇菱形地块与华南地块之间的边界带上,属于活动地块边界带的一部分(张培震等,2003),由一系列大规模逆冲断裂系组成,具有显著逆冲分量的右旋走滑性质,朝东南向推覆(闻学泽等,2013),在青藏高原隆升、川滇地块朝南东运动的区域构造作用下,会泽井所处的区域应力场呈挤压状态,其井水位在本文分析的多次中强地震前表现为上升高值异常变化(图7b)。

    从区域水文环境来看,由于处于构造灵敏点,会泽井的观测含水层对区域构造作用具有很好的反映能力。2012年云南彝良MS5.7地震和2014年云南鲁甸MS6.5地震发生在离会泽井较近的西鱼河—昭通断裂附近,从水位观测曲线上也能看出,这两次地震前出现的水位异常变化也是最显著的。但是,并不是所有地震前均会出现类似异常变化,例如2014年云南永善MS5.3地震,一方面此次地震发生在离会泽井较远的ENE向蓬峰断裂附近(近200 km),另一方面可能由于此次地震的孕育过程中会泽井所处位置未受到区域构造作用的影响,即永善MS5.3地震的震源过程与会泽井附近区域的断裂带并无直接相关性。结合先前会泽井异常变化与川滇菱形地块构造活动密切相关的认识,会泽井位于构造灵敏部位,对川滇菱形地块以东区域的整体构造作用较为敏感,但对距其较远的局部应力调整无显著反映。

    此外,观测井自身的水文条件也直接影响到水位对构造活动作用的灵敏程度。距会泽井约120 km处的昭通井、昭阳井,虽同处西鱼河—昭通断裂附近(图7a),但会泽井水位在3次地震前均有显著异常出现(图7b);昭通井在2012年云南彝良MS5.7地震和2014年云南鲁甸MS6.5地震前出现了上升高值异常,在2014年云南元谋MS4.7地震前却未出现;而昭阳井仅在2014年云南鲁甸MS6.5地震前出现了上升高值异常变化。可见,同样是水位观测井,会泽井由于处于构造灵敏点,其对地震孕育过程或构造活动的响应能力较好,即使同处一个构造区(如昭通井和昭阳井),不同的井-含水层系统条件也直接影响其对地震的响应能力。从图7b可知,昭通井水位深35 m左右,昭阳井水位埋深4.5 m左右,其观测含水层明显不同,进而对附近区域的地震响应能力也不完全相同。

    图  7  会泽及周边观测井的位置及其水位与地震的关系
    F1:莲峰断裂;F2:小江断裂;F3:西鱼河—昭通断裂;F4:磨盘山—绿汁江断裂
    Figure  7.  Location of wells and comparision of earthquakes and water levels nearby by Huize well
    F1:Lianfeng fault;F2:Xiaojiang fault;F3:Xiyuhe-Zhaotong fault;F4:Mopanshan-Lüzhijiang fault

    综上分析,井水位异常的变化机理较为复杂,既与区域构造活动有关,也与井孔自身的水文条件密切相关。会泽井位于则木河—小江断裂和西鱼河—昭通断裂的交会区附近,处于构造敏感部位,其水位对区域构造作用及地震孕育过程较为灵敏,因而具有较好的预测效能,本文定量检验的结果与其一致。

    本文通过对会泽井水位在几次中强地震前出现的显著异常特征的分析,认为会泽井水位对其周围一定范围内的中强地震具有较高的映震能力。利用Molchan图表法对会泽井水位的预测效能进行了定量检验,得到了会泽井的水位优势预测时间,最后从不同的角度对会泽井水位上升异常机理进行了初步分析与解释,结论如下:

    1) 会泽井水位在几次显著的上升异常变化后,一个月以内周围一定范围内均有中强地震发生,说明会泽井水位异常与地震活动存在一定的关系,且均为短临异常,这对周边区域的地震预测探索有重要意义。

    2) 会泽井水位异常的整体预测效能较好,对地震发生的时间和地点具有一定的指示意义。优势地震对应时间段为3个月以内,优势预测地区为川滇菱形地块以东区域。

    3) 处于断裂带交会部位的会泽井,对与震源过程相关的区域构造活动响应较为灵敏,其水位在川滇地区几次中强地震前均出现显著上升异常,与区域构造作用密切相关。

  • 图  1   江南造山带及其邻近地区地质图

    Figure  1.   Geological map of the Jiangnan orogen and its surrounding region

    图  2   远震地震射线PPdp与PPmp传播示意图

    Figure  2.   A diagrammatic view of the ray path of the teleseismic PPdp and PPmp phases

    图  3   正演模型与反演结果

    (a) 正演使用的两层水平模型;(b) 运用虚震源反射方法反演得到的时间-距离剖面

    Figure  3.   Forward model and inversion results

    (a) The two-layer horizontal model;(b) Time-distance profile obtained by the teleseismic virtual-source reflection (TVR) inversion

    图  4   事件②波形在AA′剖面台站的处理结果

    (a) 去仪器响应并通过0.1—2 Hz的带通滤波后的波形,右图为平均各道所得到的震源子波,背景中细灰色线代表用于叠加的各道波形,倒三角形标出了用于图6的台站位置;(b) 事件②虚震源反射成像剖面;(c) 事件①,②,③的虚震源反射叠加剖面;(d) 叠加后对各道进行三道圆滑的结果

    Figure  4.   Waveform processing results of event ② along the AA′ profile

    (a) Waveform data after removing the instrument response from the raw data and it was filtered by a 0.1−2 Hz band-pass filter。The averaged seismic trace is shown on the right,and the gray background curves are the overlapped traces of all channels。The inverse triangle marks the location where the station is taken for Fig. 6;(b) The original TVR profile from event ②;(c) The stacked TVR profile for teleseismic events ①,② and ③;(d) The smoothed TVR profile

    图  5   事件⑤波形在BB′剖面台站的处理结果

    (a) 去仪器响应并通过0.1—2 Hz的带通滤波后的波形,右图为平均各道所得到的震源子波,背景中细灰色线代表用于叠加的各道波形,倒三角形标出了用于图6的台站位置;(b) 事件⑤虚震源反射成像剖面;(c) 叠加事件④—⑧的虚震源反射剖面;(d) 叠加后对各道进行三道圆滑的结果

    Figure  5.   Waveform processing results of event ⑤ at the BB′ profile

    (a) Waveform data after removing the instrument response from the raw data, it was filtered by a 0.1−2 Hz band-pass filter. The averaged seismic trace is shown on the right,and the gray background curves are the overlapped traces of all channels. The inverse triangle marks the locations where the station is taken for Fig. 6;(b) The original TVR profile from event ⑤;(c) The stacked TVR profile for teleseismic events ④−⑧;(d) The smoothed TVR profile

    图  6   远震虚源反射法与传统接收函数结果对比

    (a) 事件②在台站 (图4a中倒三角形)处的虚源反射剖面(上图)与合成信号(下图)的对比,箭头代表识别的震相,合成信号所用参数为vP为5.2 km/s,vS为3.1 km/s,反射层深度为6.5 km;(b) 事件②在台站(图4a中倒三角形)处Ps转换波接收函数(箭头)(上图)与人工合成信号(下图)的对比,合成信号所用参数与图(a)相同;(c) 与图(a)相似,事件⑤在台站(图5a中倒三角形)处的对比,合成信号所用参数为vP为5.2 km/s,vS为3.1 km/s,反射层深度为5 km;(d) 与图(b)相似,事件⑤在台站(图5a中倒三角形)处的对比,合成信号所用参数与图(c)相同

    Figure  6.   Example of TVR profile (top) for events ② and ⑤ at the marked station (inverse triangle in Figs. 4a and 5a) versus synthetic seismogram (bottom)

    (a) Example of TVR profile (top) for the event ② at the marked station (inverse triangle in Fig. 4a) versus synthetic seismogram (bottom). The arrow represents the recognized seismic phase. The synthetics were calculated to fit the data using vP of 5.2 km/s and vS of 3.1 km/s and the depth of the seismic layer generating reflection phase is 6.5 km. (b) Example of conventional Ps receiver function (top) at the marked station (inverse triangle in Fig. 4a) compared with the synthetic seismogram (bottom). The parameters used to synthesize the signal are the same as those in Fig. (a). (c) Example of TVR profile (top) for the event ⑤ at the marked station (inverse triangle in Fig. 5aversus synthetic seismogram (bottom). The arrow represents the recognized seismic phase. The synthetics were calculated to fit the data using vP of 5.2 km/s and vS of 3.1 km/s and the depth of the seismic layer generating reflection phase is 5 km. (d) Example of conventional Ps receiver function (top) at the marked station (inverse triangle in Fig. 5a) compared with the synthetic seismogram (bottom). The parameters used to synthesize the signal are the same as those in Fig. (c)

    图  7   结合已知的地层信息与断层对成像剖面AA′ (a)和 BB′ (b)的解释

    两剖面以纵波速度5.2 km/s进行时深转换结果。绿色线表示沉积层,蓝色和红色线代表上地壳内两段反射面

    Figure  7.   Interpretation of profiles AA′ (a) and BB′ (b) from the interfaces and known faults

    The two profiles are converted using average P-wave velocity of 5.2 km/s. The green lines represent sedimentary layers,and the blue and red lines represent two reflections in the shallow crust

    表  1   本研究所用地震事件的参数

    事件
    编号
    发震日期 MS 地理位置 震源深度
    /km
    震中距
     年-月-日 时:分:秒
    2019-04-03 05:35:33 6.4 (52.25°N,178.00°E) 20 50.711
    2019-04-07 05:55:00 6.1 (6.90°S,125.07°E) 540 36.746
    2019-04-12 22:51:32 5.8 (6.43°S,148.57°E) 20 47.157
    2019-06-03 14:04:36 5.5 (0.35°N,97.72°E) 20 35.244
    2019-06-15 04:10:52 5.6 (5.86°S,130.75°E) 120 37.696
    2019-06-16 06:55:00 7.2 (30.80°S,178.10°W) 20 85.704
    2019-06-16 13:17:15 6.2 (30.97°S,178.06°W) 30 85.836
    2019-06-17 14:02:04 5.9 (30.85°S,177.60°W) 20 86.071
    注:震源参数引自中国地震台网中心
    下载: 导出CSV
  • 黄汲清,任纪舜,姜春发,张之孟,徐志琴. 1977. 中国大地构造基本轮廓[J]. 地质学报,51(2):117–135.

    Huang J Q,Ren J S,Jiang C F,Zhang Z M,Xu Z Q. 1977. An outline of the characteristics of China[J]. Acta Geology Sinica,51(2):117–135 (in Chinese).

    李吉均,张林源,邓养鑫,周尚哲. 1983. 庐山第四纪环境演变和地貌发育问题[J]. 中国科学:B辑,13(8):734–745.

    Li J J,Zhang L Y,Deng Y X,Zhou S Z. 1983. Problems on Quaternary environmental evolution and geomorphic development in the Lushan Mountain[J]. Science in China:Series B,13(8):734–745 (in Chinese).

    梁锋,吕庆田,严加永,刘振东. 2014. 长江中下游宁芜火山岩盆地深部结构特征:来自反射地震的认识[J]. 岩石学报,30(4):941–956.

    Liang F,Lü Q T,Yan J Y,Liu Z D. 2014. Deep structure of Ningwu volcanic basin in the middle and lower reaches of Yangtze River:Insights from reflection seismic data[J]. Acta Petrologica Sinica,30(4):941–956 (in Chinese).

    梁兴,叶舟,吴根耀,郑华平,徐克定,张廷山,刘家铎. 2006. 鄱阳盆地构造-沉积特征及其演化史[J]. 地质科学,41(3):404–429.

    Liang X,Ye Z,Wu G Y,Zheng H P,Xu K D,Zhang T S,Liu J D. 2006. Sedimento-tectonic features and geological evolution of the Poyang basin[J]. Chinese Journal of Geology,41(3):404–429 (in Chinese).

    舒良树,周新民,邓平,余心起. 2006. 南岭构造带的基本地质特征[J]. 地质论评,52(2):251–265.

    Shu L S,Zhou X M,Deng P,Yu X Q. 2006. Principal geological features of Nanling tectonic belt,South China[J]. Geological Review,52(2):251–265 (in Chinese).

    舒良树. 2012. 华南构造演化的基本特征[J]. 地质通报,31(7):1035–1053.

    Shu L S. 2012. An analysis of principal features of tectonic evolution in South China block[J]. Geological Bulletin of China,31(7):1035–1053 (in Chinese).

    王孝磊,周金城,陈昕,张凤凤,孙梓铭. 2017. 江南造山带的形成与演化[J]. 矿物岩石地球化学通报,36(5):714–735.

    Wang X L,Zhou J C,Chen X,Zhang F F,Sun Z M. 2017. Formation and evolution of the Jiangnan Orogen[J]. Bulletin of Mineralogy,Petrology and Geochemistry,36(5):714–735 (in Chinese).

    徐涛,张忠杰,田小波,刘宝峰,白志明,吕庆田,滕吉文. 2014. 长江中下游成矿带及邻区地壳速度结构:来自利辛—宜兴宽角地震资料的约束[J]. 岩石学报,30(4):918–930.

    Xu T,Zhang Z J,Tian X B,Liu B F,Bai Z M,Lü Q T,Teng J W. 2014. Crustal structure beneath the middle-lower Yangtze metallogenic belt and its surrounding areas:Constraints from active source seismic experiment along the Lixin to Yixing profile in East China[J]. Acta Petrologica Sinica,30(4):918–930 (in Chinese).

    杨晓东,吴中海,张海军. 2016. 鄱阳湖盆地的地质演化、新构造运动及其成因机制探讨[J]. 地质力学学报,22(3):667–684.

    Yang X D,Wu Z H,Zhang H J. 2016. Geological evolution,neotectonics and genetic mechanism of the Poyang lake basin[J]. Journal of Geomechanics,22(3):667–684 (in Chinese).

    张培震,邓起东,张国民,马瑾,甘卫军,闵伟,毛凤英,王琪. 2003. 中国大陆的强震活动与活动地块[J]. 中国科学:D辑,33(增刊1):12–20.

    Zhang P Z,Deng Q D,Zhang G M,Ma J,Gan W J,Min W,Mao F Y,Wang Q. 2003. Active tectonic blocks and strong earthquakes in the continent of China[J]. Science in China:Series D,46(2):13–24.

    周新民,朱云鹤. 1992. 江绍断裂带的岩浆混合作用及其两侧的前寒武纪地质[J]. 中国科学:B辑,22(3):296–303.

    Zhou X M,Zhu Y H. 1992. Magma mixing in the Jiang-Shao fault zone and the Precambrian geology of both side of the Jiang-Shao zone[J]. Science in China:Series B,22(3):296–303 (in Chinese).

    Chen M,Niu F L,Liu Q Y,Tromp J,Zheng X F. 2015. Multiparameter adjoint tomography of the crust and upper mantle beneath East Asia:1. Model construction and comparisons[J]. J Geophys Res:Solid Earth,120(3):1762–1786. doi: 10.1002/2014JB011638

    Guo L H,Gao R,Shi L,Huang Z R,Ma Y W. 2019. Crustal thickness and Poisson’s ratios of South China revealed from joint inversion of receiver function and gravity data[J]. Earth Planet Sci Lett,510:142–152. doi: 10.1016/j.jpgl.2018.12.039

    Guo L Z,Shi Y S,Lu H F,Ma R S,Dong H G,Yang S F. 1989. The pre-Devonian tectonic patterns and evolution of South China[J]. J Southeast Asian Earth Sci,3(1/2/3/4):87–93.

    Guo Z,Gao X,Li T,Wang W. 2018. Crustal and uppermost mantle structures of the South China from joint analysis of receiver functions and Rayleigh wave dispersions[J]. Phys Earth Planet Inter,278:16–25. doi: 10.1016/j.pepi.2018.03.001

    Hsü K J,Li J L,Chen H H,Wang Q C,Sun S,Şengör A M C. 1990. Tectonics of South China:Key to understanding West Pacific geology[J]. Tectonophysics,183(1/2/3/4):9–39.

    Huang J L,Zhao D P. 2006. High-resolution mantle tomography of China and surrounding regions[J]. J Geophys Res,111(B9):B09305.

    Huang R,Xu Y X,Zhu L P,He K. 2015. Detailed Moho geometry beneath southeastern China and its implications on thinning of continental crust[J]. J Asian Earth Sci,112:42–48. doi: 10.1016/j.jseaes.2015.09.002

    Lü Q T,Yan J Y,Shi D N,Dong S W,Tang J T,Wu M G,Chang Y F. 2013. Reflection seismic imaging of the Lujiang-Zongyang volcanic basin,Yangtze Metallogenic Belt:An insight into the crustal structure and geodynamics of an ore district[J]. Tectonophysics,606:60–77. doi: 10.1016/j.tecto.2013.04.006

    Li H Y,Song X D,Lü Q T,Yang X Y,Deng Y F,Ouyang L B,Li J P,Li X F,Jiang G M. 2018. Seismic imaging of lithosphere structure and upper mantle deformation beneath east-central China and their tectonic implications[J]. J Geophys Res,123(4):2856–2870. doi: 10.1002/2017JB014992

    Li X Q,Nábĕlek J L. 1999. Deconvolution of teleseismic body waves for enhancing structure beneath a seismometer array[J]. Bull Seismol Soc Am,89(1):190–201.

    Ligorría J P,Ammon C J. 1999. Iterative deconvolution and receiver-function estimation[J]. Bull Seismol Soc Am,89(5):1395–1400.

    Lin S F,Xing G F,Davis D W,Yin C Q,Wu M L,Li L M,Jiang Y,Chen Z H. 2018. Appalachian-style multi-terrane Wilson cycle model for the assembly of South China[J]. Geology,46(4):319–322. doi: 10.1130/G39806.1

    Luo S,Yao H J,Li Q S,Wang W T,Wan K S,Meng Y F,Liu B. 2019. High-resolution 3D crustal S-wave velocity structure of the Middle-Lower Yangtze River Metallogenic Belt and implications for its deep geodynamic setting[J]. Science China Earth Sciences,62(9):1361–1378. doi: 10.1007/s11430-018-9352-9

    Luo Y H,Xu Y X,Yang Y J. 2012. Crustal structure beneath the Dabie orogenic belt from ambient noise tomography[J]. Earth Planet Sci Lett,313/314:12–22. doi: 10.1016/j.jpgl.2011.11.004

    Shen W S,Ritzwoller M H,Kang D,Kim Y H,Lin F C,Ning J Y,Wang W T,Zheng Y,Zhou L Q. 2016. A seismic reference model for the crust and uppermost mantle beneath China from surface wave dispersion[J]. Geophys J Int,206(2):954–979. doi: 10.1093/gji/ggw175

    Shu L S,Faure M,Yu J H,Jahn B M. 2011. Geochronological and geochemical features of the Cathaysia block (South China):New evidence for the Neoproterozoic breakup of Rodinia[J]. Precambrian Res,187(3/4):263–276.

    Tong P,Chen C W,Komatitsch D,Basini P,Liu Q Y. 2014. High-resolution seismic array imaging based on an SEM-FK hybrid method[J]. Geophys J Int,197(1):369–395. doi: 10.1093/gji/ggt508

    Tseng T L,Chen W P. 2006. Probing the southern Indian shield with P-wave receiver-function profiles[J]. Bull Seismol Soc Am,96(1):328–333. doi: 10.1785/0120050074

    Tseng T L,Chen W P,Nowack R L. 2009. Northward thinning of Tibetan crust revealed by virtual seismic profiles[J]. Geophys Res Lett,36(24):L24304. doi: 10.1029/2009GL040457

    Wang D Z,Shu L S. 2012. Late Mesozoic basin and range tectonics and related magmatism in Southeast China[J]. Geosci Front,3(2):109–124. doi: 10.1016/j.gsf.2011.11.007

    Yang Z H,Sheehan A F,Yeck W L,Miller K C,Erslev E A,Worthington L L,Harder S H. 2012. Imaging basin structure with teleseismic virtual source reflection profiles[J]. Geophys Res Lett,39(2):L02303.

    Yu J H, Schuster G T. 2001. Crosscorrelogram migration of IVSPWD Data[C]//2001 SEG Annual International Meeting. San Antonio, Texas: Society of Exploration Geophysicists: 456–459.

    Zhao G C. 2015. Jiangnan Orogen in South China:Developing from divergent double subduction[J]. Gondwana Res,27(3):1173–1180. doi: 10.1016/j.gr.2014.09.004

    Zhou X M,Li W X. 2000. Origin of Late Mesozoic igneous rocks in Southeastern China:Implications for lithosphere subduction and underplating of mafic magmas[J]. Tectonophysics,326(3/4):269–287.

    Zou Z H,Zhou H W,Gurrola H. 2014. Teleseismic virtual‐source imaging of the basin structures in the Three Gorges region,China[J]. Bull Seismol Soc Am,104(4):2142–2147. doi: 10.1785/0120130110

图(7)  /  表(1)
计量
  • 文章访问数:  1008
  • HTML全文浏览量:  557
  • PDF下载量:  54
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-01-01
  • 修回日期:  2020-06-07
  • 网络出版日期:  2020-12-07
  • 发布日期:  2020-09-14

目录

/

返回文章
返回