坚固体孕震模式下地下流体异常时空演化的数值模拟
-
摘要: 利用固液二相介质平面应变问题的有限元方法,研究了倾滑型地震地下流体场兆、源兆特征及时空演变过程.结果表明:孕震的早期阶段,即弹性积累和非线性阶段,地下流体缓慢变化,异常不明显;中期阶段,即硬化及局部膨胀阶段,震源区出现膨胀,地下流体异常主要集中在震源区,源兆出现得早,场兆出现得晚;进入中短期阶段,即大范围膨胀阶段,震源区异常继续发展,异常区不断扩大,在震源区以外出现新的异常值较大的地区.
-
关键词:
- 地下流体 场兆 源兆 机理 数值模拟
-
引言
张家口—渤海地震活动带(以下简称为张渤地震带)穿北京、天津延至渤海,大致走向为WNW,长约600 km,宽十几至数十千米,是我国东部的一条重要的地震活动带,曾发生1976年唐山MS7.8、1975年海城MS7.3大地震和大量中小地震(图 1).其东段与郯城—庐江断裂带(以下简称为郯庐断裂带)的北段相交.本文研究区范围为113°E—126°E、34°N—44°N,地处华北,是我国经济文化最发达、人口最密集的地区之一. 对张渤地震带及其邻区的地震特性研究可为该地区的地震预测提供深部地球物理依据,同时为该地区的可持续发展提供参考.
图 1 研究区域构造背景以及台站和1970年以来M>4.5地震震中分布Figure 1. Tectonic settings as well as distribution of stations (triangles) and M>4.5 earthquakes (red dots) since the year 1970 in the research areaThe belt extending in NW- -SE direction is Zhangjiakou-Bohai seismic active belt, the other is Tancheng-Lujiang fault zone张渤地震带及其邻区位于华北克拉通中东部地区(图 1).华北克拉通是世界上最古老的克拉通之一(Wilde et al,2001;翟明国,2010),古生代以前保持着典型克拉通的性质(赵国春等,2002;Menzies et al,2007),之后其东部的岩石圈地幔经历了物理和化学性质的根本改变. 相关研究认为华北克拉通东部已经被改造,甚至发生了破坏(邓晋福等,2006; Zheng et al,2007;嵇少丞等,2008;An et al,2009;朱日祥,郑天愉,2009).诸多研究人员对该区域的构造开展了大量的研究工作,包括利用人工地震测深方法对壳幔结构的研究(王夫运等,2005;熊小松等,2011)、体波和面波的层析成像研究(朱介寿等,2002;Huang et al,2003;Liang et al,2004;Sun,Toksöz,2006)、利用接收函数对间断面的研究(Ai,Zheng,2003; Li et al,2008;Chen,Ai,2009; Gao et al,2010)等,对揭示该区域的构造演化具有非常重要的意义,同时也为后续研究提供参考.
地球内部的间断面是探讨地球结构、界定深部构造性质的基础(高原,周蕙兰,1998).莫霍面是地球最重要的界面之一,表征着地震波速度和岩石物性参数的改变,其深度是地壳结构的一个重要参数.地壳平均泊松比可以揭示岩石组分的变化(Ji et al,2009).对莫霍面起伏特征和地壳平均泊松比变化的研究,可为揭示该地区地壳形变和深部构造提供重要信息.
地震图中远震P波通过莫霍面时产生的P- -S转换波非常微弱而难以辨识,而接收函数对于突出这一转换波非常有效,故成为研究地壳结构的理想方法(吴庆举,曾融生,1998;张莹莹,高原,2012).Zandt和Ammon(1995)首先利用接收函数中的转换震相及多次反射转换震相估计了全球地壳厚度和泊松比;Zhu和Kanamori(2000)对该方法作了进一步约束,使其更加稳定.之后该方法在地球物理领域得到了广泛的应用(李永华等,2006;张洪双等,2009;Wang et al,2010;张广成等,2013),也被称为接收函数的H-k叠加搜索方法.
本研究采用接收函数的H-k叠加搜索方法,计算张渤地震带及其邻区台站下方的地壳厚度和泊松比.由于接收函数方法的稳定性依赖于地震事件的反方位角分布,本研究采用较罗艳等(2008)及葛粲等(2011)更长时间的波形记录,以期获得更多的远震事件和更为丰富的方位角分布,从而保证本文结果的稳定性,有助于进一步了解华北地区及张渤地震带和郯庐断裂带北段的深部结构.
2. 数据及方法
本文采用位于山东、河北、辽宁、北京和天津等地区的103个宽频带固定台站(图 1)于2007年8月1日—2011年12月31日(共4年零4个月)记录到的波形资料,选取震中距为30°—90°、M>5.5的地震事件1 000余个(图 2).在此基础上,采用带宽为0.02—0.5 Hz的频率域反褶积方法计算接收函数,从中筛选出转换震相清晰且信噪比高的接收函数,用于接收函数的H-k叠加搜索方法计算.
接收函数的H-k叠加搜索方法利用了接收函数的时间和振幅属性.在台站下方莫霍面起伏和速度横向变化不大的情况下,莫霍面产生的转换波Ps以及多次波PpPs和PpSs+PsPs的走时仅与台站下方的莫霍面深度、P波与S波速度比有关.当地壳平均P波速度(vP)、波速比(k=vP/vS)和地壳厚度(H)的取值与实际地壳情况相符时,转换波和多次波的到时所对应的振幅加权叠加值最大.在给定地壳平均P波速度的情况下,对每个台站的接收函数进行扫描叠加,可以得到每个台站下方的地壳厚度和平均速度比.扫描叠加振幅由震相Ps,PpPs和PpSs+PsPs的振幅加权值给出,权值和为1.在Zhu和Kanamori(2000)的研究中,Ps,PpPs和PpSs+PsPs振幅的权值通常分别采用0.7,0.2和0.1,本文通过多次试算比较认为这样的权值比较合理.根据人工测深结果(刘昌铨,杨健,1982;嘉世旭等,2005),计算中我们取地壳平均P波速度为6.3 km/s,并设定H和k的扫描范围分别为30—50 km和1.6—1.9,相应扫描间隔分别为1 km和0.005,扫描结果的误差根据叠加函数的方差给定.
3. 结果及分析
通过对远震事件波形处理,得到了各个台站的接收函数,并对其中Ps及多次反射震相清晰、信噪比高的99个台站的接收函数进行了地壳厚度和波速比扫描,图 3给出了DSA,BXI,XIT等3个台站的H-k扫描分析结果.利用泊松比与波速比之间的换算关系σ=0.5[1-(k2-1)-1],将叠加结果中的波速比换算为泊松比,结果列于表 1.
图 3 利用H-k方法计算台站下方的地壳厚度(H)和平均波速比(νP/νS)a)H-k搜索结果;(b)接收函数按射线参数叠加的结果,同时给出了台站名称及参与计算的接收函数叠加条数Figure 3. Crustal thickness and the distribution of νP/νS ratio beneath the stations calculated by using H-k method(a)The searching result by H-k method; (b)Receiver functions stacking according to ray parameters. The name of stations and the number of receive functions are also given on the end of receive function waves表 1 各台站H-k扫描结果Table 1. H-k searching results for all stations used in this study台站 H/km k σ N 台站 H/km k σ N JIX 34±1.84 1.75±0.039 0.18 79 HAY 32±1.23 1.78±0.042 0.27 120 BHC 35±1.01 1.62±0.028 0.19 54 HUR 33±1.17 1.77±0.057 0.27 106 RCH 32±0.81 1.64±0.033 0.2 177 JIX 33±1.69 1.79±0.062 0.27 128 HXQ 33±2.07 1.65±0.058 0.21 61 KAB 40±1.23 1.78±0.038 0.27 174 LIS 34±1.43 1.65±0.033 0.21 113 KDN 32±1.34 1.77±0.053 0.27 176 JNX 37±2.04 1.67±0.345 0.22 61 LQU 34±1.01 1.79±0.028 0.27 135 QIL 34±1.22 1.67±0.035 0.22 61 LUQ 33±0.95 1.77±0.027 0.27 105 ZCH 36±1.68 1.66±0.042 0.22 201 LYA 36±2.36 1.7±0.044 0.27 175 JZH 36±1.50 1.68±0.033 0.23 152 LZH 33±1.05 1.77±0.042 0.27 130 XTT 32±1.40 1.68±0.039 0.23 129 MQI 31±1.17 1.77±0.045 0.27 156 DSD 34±1.24 1.7±0.044 0.24 149 NKY 36±1.23 1.78±0.031 0.27 183 KUC 35±1.52 1.71±0.046 0.24 89 QID 31±0.78 1.78±0.023 0.27 109 LOK 34±1.29 1.7±0.045 0.24 196 WED 30±0.9 1.79±0.03 0.27 92 LSH 35±1.15 1.72±0.042 0.24 179 WFD 33±1.05 1.77±0.036 0.27 161 SHS 32±1.56 1.71±0.051 0.24 102 WUL 35±1.63 1.77±0.049 0.27 140 TLK 33±1.28 1.7±0.038 0.24 75 XIL 34±0.92 1.79±0.035 0.27 61 XIT 35±1.30 1.71±0.035 0.24 139 XLD 33±1.06 1.78±0.038 0.27 109 XYN 35±0.96 1.71±0.028 0.24 140 ANS 30±2.31 1.8±0.135 0.28 210 YON 30±2.49 1.71±0.063 0.24 112 BXI 34±1.79 1.81±0.066 0.28 52 YTA 33±1.66 1.72±0.044 0.24 63 CHD 34±1.05 1.8±0.036 0.28 92 CHC 40±1.23 1.74±0.044 0.25 166 CXT 34±1.07 1.82±0.035 0.28 89 GAX 32±1.69 1.74±0.057 0.25 31 DHC 36±1.75 1.8±0.055 0.28 311 HNS 31±1.41 1.73±0.045 0.25 76 GUS 32±2.05 1.81±0.073 0.28 115 HST 34±1.54 1.73±0.048 0.25 222 HSH 31±1.02 1.82±0.043 0.28 90 JCA 35±2.23 1.73±0.042 0.25 256 JIN 33±1.45 1.8±0.055 0.28 179 JUN 33±0.81 1.73±0.027 0.25 164 LHT 31±1.78 1.8±0.04 0.28 91 JUX 33±1.87 1.73±0.036 0.25 83 TCH 32±1.14 1.82±0.037 0.28 75 LOH 36±1.37 1.74±0.028 0.25 148 XFN 30±1.55 1.82±0.041 0.28 25 QYU 33±1.43 1.73±0.228 0.25 62 XIT 33±1.24 1.81±0.048 0.28 68 RSH 32±1.53 1.74±0.067 0.25 139 XMN 30±1.52 1.8±0.035 0.28 113 SNY 32±0.98 1.74±0.042 0.25 320 ZHB 39±1.22 1.82±0.037 0.28 136 XBZ 43±1.33 1.73±0.03 0.25 77 ZHT 36±1.49 1.81±0.046 0.28 232 YKO 30±1.26 1.74±0.037 0.25 86 BBS 34±1.56 1.85±0.038 0.29 41 YUS 34±1.08 1.73±0.033 0.25 159 CHY 36±1.25 1.83±0.036 0.29 19 CLI 32±1.43 1.76±0.062 0.26 58 CSH 31±1.62 1.84±0.064 0.29 191 FEN 37±1.54 1.75±0.029 0.26 137 FKU 30±1.99 1.85±0.039 0.29 94 LIC 33±1.50 1.76±0.047 0.26 191 LAY 30±1.66 1.83±0.038 0.29 75 LYN 33±1.00 1.75±0.034 0.26 187 LBP 35±0.70 1.84±0.022 0.29 118 NLA 34±2.00 1.75±0.075 0.26 142 MIY 34±1.17 1.85±0.039 0.29 195 TIL 31±1.08 1.75±0.037 0.26 240 NAP 32±1.55 1.83±0.06 0.29 145 WAT 34±1.02 1.75±0.027 0.26 154 RZH 33±1.01 1.84±0.034 0.29 110 WEC 37±1.79 1.75±0.055 0.26 134 SSL 34±1.11 1.84±0.035 0.29 212 YSH 35±1.13 1.76±0.038 0.26 118 TIA 30±1.54 1.85±0.039 0.29 106 ANQ 34±1.41 1.78±0.027 0.27 101 BEP 31±1.56 1.87±0.052 0.30 178 CHD 32±0.94 1.78±0.034 0.27 90 BZH 31±1.46 1.86±0.053 0.30 163 DDO 31±1.38 1.78±0.043 0.27 145 QIX 30±1.92 1.86±0.045 0.30 60 DL2 30±1.26 1.77±0.038 0.27 129 SXT 34±0.83 1.87±0.026 0.30 102 FXI 33±1.18 1.77±0.042 0.27 189 ZJK 36±1.35 1.87±0.035 0.30 60 GAN 36±1.66 1.78±0.055 0.27 102 YAY 34±1.23 1.89±0.034 0.31 81 注:H为莫霍面深度,k为波速比,σ为泊松比,N为地震数. 3.1 地壳厚度
由表 1可以看出,张渤地震带及其邻区地壳厚度为30—46 km,变化跨度很大,与王峻等(2009)以及许卫卫和郑天愉(2005)的研究结果在空间分布上基本一致.图 4和图 5分别给出了用插值和按台站位置两种方式展示的地壳厚度和泊松比分布情况.从整体来看(图 4),研究区地壳厚度与地表地形有很好的对应性,山脉隆起区地壳厚度较大,盆地和平原地区地壳厚度较小,这与艾利均衡假说一致,可能暗示了研究区的地壳经历了破坏与均衡共同作用的过程.由于华北平原缺少宽频带地震台,本研究在该地区无法给出计算结果,故该地区的差值结果可靠性不佳. 相关研究(葛粲等,2011;刘琼林等,2011)在该地区也未得到可靠结果.从图 4可以看出,华北平原北部的山脉隆起地区的地壳厚度介于36—43 km之间,具有从造山带到平原或盆地逐渐减薄的变化趋势,这与赵金仁等(2006)的人工源反射剖面研究结果一致.
从局部来看(图 5),燕山隆起区地壳较厚,最厚处位于康保台(KAB)和张北台(ZHB),其下方地壳厚度为40 km左右,地壳厚度从张渤地震带的NW- -SE方向逐渐减薄(图 4),这与嘉世旭等(2005)、许卫卫和郑天愉(2005)和和王峻等(2009)的研究结果一致.从图 4可以看出: 研究区内的华北克拉通中、西部地区平均地壳厚度大于40 km; 燕山隆起区地壳整体较厚,并且西北部地区较东南部地区略厚.本文结果显示,北京地区出现了莫霍面局部高低起伏的现象,这可能与北京地区位于燕山隆起与太行山隆起交汇处的复杂构造有关. 该现象在罗艳等(2008)、王峻等(2009)、 葛粲等(2011)及刘琼林等(2011)的研究结果中并未显示.在华北平原向太行山隆起的过渡地区出现了莫霍面陡变,最厚处井陉台(JNX)下方地壳厚度高达37 km,最薄处永年台(YON)仅为30 km.结果还显示,分布在太行山山前断裂两侧的邻近台站地壳厚度表现出剧烈的变化,例如井陉台(JNX)和鹿泉台(LUQ)对应的莫霍面深度分别为37和33 km,这暗示着太行山山前断裂局部可能深达莫霍面.郯庐断裂带南段的山东地区地壳厚度较厚,而该断裂带北端的辽宁 地区地壳较薄且变化不大,平均厚度为32.3 km,这与葛粲等(2011)的研究结果一致.故整体来看: 郯庐断裂带下方地壳呈南厚北薄的特点;辽东隆起地区山脊处地壳厚度较两侧略厚; 鲁西和鲁东隆起部分地区地壳厚度略厚,若此处的克拉通破坏方式以底侵作用为主(翟明国等,2005),则在地幔岩浆底侵作用下,地壳和软流圈地幔的岩石发生混存而使该地区地壳变厚;位于渤海中间岛屿上的北隍城台(BHC),其下方地壳厚度较大,约35 km.
3.2 泊松比
研究区的泊松比与地形和地壳厚度虽没有明显的对应关系,却呈现出明显的区域分布特征.大陆地壳的泊松比一般为0.25—0.27(Zandt,Ammon,1995; 童蔚蔚等,2007),而位于张渤地震带西北段断层密集区域的台站下方地壳泊松比值均较高(图 6). 例如张家口台(ZJK)、十三陵台(SSL)、八宝山台(BBS)、牛口峪台(NKY)和大灰场台(DHC),其下方泊松比均大于0.27.一般来说较高泊松比意味着岩石破碎、多空裂隙,而低泊松比则说明岩石相对完整、强度大.断层密集区域附近的台站下方地壳出现高泊松比,其原因一是由于断裂带附近岩石比较破碎,二是在岩石的破碎过程中应力释放导致局部应力变低,从而造成地壳岩石的部分熔融所致.太行山山前地区泊松比变化剧烈,这可能是由于太行山山前地区地壳结构比较复杂,地壳物质非均匀分布所致.根据野外调查结果(张瑾等,2005;张锡明等,2007),鲁西地区广泛发育典型的幔源岩浆岩,其中中生代侵入岩分布十分广泛且类型多样,这些幔源岩浆岩大多属于基性、超基性岩,具有较高的泊松比.本文认为鲁西隆起地区出现泊松比高低复杂分布的现象,可能是由于地幔物质的底侵作用 所引起. 底侵作用伴随着局部的火山活动,将高泊松比的地幔物质带入地壳,导致局部地壳泊松比增大.鲁西隆起南部的苍山台(CSH)下方地壳泊松比值高达0.29,东部的泰安台(TIA)和济南台(JIN)泊松比值也分别高达0.29和0.28. 臧绍先等(2002)的大地热流结果表明这些台站位于大地热流高值区,故本文认为该地区地壳可能有局部熔融的现象.郯城台(TCH)下方地壳泊松比值高达0.28,这是由于该台站位于郯庐断裂带起始位置,断层处岩石较破碎所致; 而远离该断层的邹城台(ZCH)和独山岛台(DSD)下方地壳泊松比值分别为0.22和0.24,岩石完整强度较大,所以泊松比值较低;临沭台(LIS)下方地壳泊松比值为0.21. 鲁东地区地表广泛出露榴辉岩和超基性岩(凌贤长等,1996),本文认为鲁东隆起的高泊松比与地壳岩性有关.由图 7可看出,位于渤海中间岛屿上的北隍城台(BHC)下方地壳泊松比值较低,为0.19. 杨婷等(2012)的成像结果显示该区域45 km深度处为低速(但该处为研究区边缘,结果可靠性不佳),据此分析该区泊松比较低可能是热物质上涌所致.
4. 讨论与结论
通过对99个台站的接收函数H-k叠加搜索计算,结果表明张渤地震带及其邻区的地壳厚度在30—46 km之间变化,跨度很大,且地壳厚度与地表地形相关性很好,即山脉隆起区地壳厚度大,盆地区地壳厚度小,地壳厚度在鲁西和鲁东隆起处较大,地壳从燕山、太行山地区到华北平原逐渐变薄;辽东隆起地区地壳厚度则由山脊向两侧山坡逐渐减薄.太行山隆起的东缘附近,地壳厚度变化剧烈,暗示太行山山前断裂很可能是一条深达莫霍面的断裂构造.位于郯庐断裂带北段渤海湾北部的辽宁地区地壳较薄,平均厚度为32.3 km;郯庐断裂带南段的山东地区,地壳厚度较厚. 整体来看,位于研究区内的郯庐断裂带下方地壳厚度呈现南厚北薄的特点,这是由于郯庐断裂带恰好在此处穿过地壳较厚的鲁西隆起地区,而地壳厚度在郯庐断裂带两侧没有明显变化,可能暗示该区域地壳在纵向上没有经历明显的压缩和拉伸作用.
本文得到的研究区泊松比值处于0.18—0.31范围内,其分布具有分区特点.在张渤地震带西北段断层比较密集的区域,部分台站下方地壳泊松比值大于0.27,表现为高值异常,可能是地震活动造成岩石比较破碎所致; 而应力释放导致的地壳岩石部分熔融也可能是引起高泊松比的另一个重要原因.在郯庐带东侧的鲁西隆起出现了高低泊松比复杂分布的现象. 泊松比高值异常出现在火山活动区,对应的地表有大面积侵入岩出露,其原因可能是局部的火山活动将高泊松比的地幔物质带入地壳,从而导致地壳泊松比的增大;泊松比低值异常出现在远离断层的地区,推测岩石较为完整,强度较大.太行山山前地区泊松比变化剧烈,可能是由于该地区地壳结构比较复杂,地壳物质非均匀分布所致.位于渤海中间岛屿上的北隍城台(BHC)地壳厚度较大,为35 km,泊松比较低,为0.19,其原因尚不清楚,有待于今后综合更多的观测资料和借助多种分析技术对该地区开展更为精细的研究.
中国地震局地球物理研究所“国家数字测震台网数据备份中心”为本研究提供了地震波形数据,郑秀芬研究员在资料的收集和处理过程中提供了帮助,钮凤林教授和张洪双博士在数据处理过程中给予了指导,作者在此一并表示衷心的感谢.
-
-
期刊类型引用(14)
1. 关成尧,姜纪沂,袁四化,白相东,刘晓燕,张艳. 密云水库及邻区地貌成因、控盆要素及其区域构造意义. 城市地质. 2025(01): 47-59 . 百度学术
2. 张扬,商琳,王淼,吴海涛,段彬. 南堡凹陷地壳构造稳定性:源自地震活动时空特征的证据. 地质力学学报. 2025(02): 313-324 . 百度学术
3. Mengke Cai,Gongcheng Zhang,Wanyin Wang,Dingding Wang. Study on crustal thickness and the prediction of prolific depressions: the Bohai Basin as an example. Acta Oceanologica Sinica. 2024(04): 92-104 . 必应学术
4. 孙安辉,高原,赵国峰,任超,梁姗姗. 2022年1月8日青海门源6.9级地震的震源区结构特征和b值意义初探. 地球物理学报. 2022(03): 1175-1183 . 百度学术
5. 刘嘉栋,丁志峰,武岩,姜磊. 利用远震接收函数研究华北克拉通北部造山带地壳厚度及泊松比. 地震学报. 2022(03): 357-373 . 本站查看
6. 刘迁,赵瑞,栗宁. 利用远震接收函数探测襄樊台下方地壳厚度及泊松比. 地震科学进展. 2021(05): 230-233 . 百度学术
7. 李婷婷,刘利,范文华,孙业君. 基于远震P波接收函数研究江苏地区地壳厚度和泊松比. 地震研究. 2020(04): 680-688+768 . 百度学术
8. 吴鹏,贾华,张小涛,王雪飞,刘爽. 晋冀蒙交界地区的地壳厚度与泊松比特征. 地震. 2019(04): 63-75 . 百度学术
9. 李静,周晓成,石宏宇,孙凤霞,陆丽娜,陈超,曾令华,陈志. 首都圈西北部主要活动断裂CO_2、Rn、Hg脱气对环境的影响. 环境化学. 2018(05): 931-941 . 百度学术
10. 许英才,曾宪伟,许文俊,马禾青,金涛,任家琪. 基于台阵的青藏高原东北缘海原-六盘山断裂带及邻区地壳结构研究. 中国地震. 2018(03): 484-497 . 百度学术
11. 王喜龙,李营,杜建国,陈志,周晓成,李新艳,崔月菊,王海燕,张志宏. 首都圈地区土壤气Rn, Hg, CO_2地球化学特征及其成因. 地震学报. 2017(01): 85-101+155 . 本站查看
12. Guang-hua Pang,Ji-Kun Feng,Jun Lin. Crustal structure beneath Liaoning province and the Bohai Sea and its adjacent region in China based on ambient noise tomography. Earthquake Science. 2017(01): 1-15 . 必应学术
13. WANG ChunYong,WU QingJu,DUAN YongHong,WANG ZhiShuo,LOU Hai. Crustal and upper mantle structure and deep tectonic genesis of large earthquakes in North China. Science China(Earth Sciences). 2017(05): 821-857 . 必应学术
14. 王椿镛,吴庆举,段永红,王志铄,楼海. 华北地壳上地幔结构及其大地震深部构造成因. 中国科学:地球科学. 2017(06): 684-719 . 百度学术
其他类型引用(9)
计量
- 文章访问数: 1211
- HTML全文浏览量: 18
- PDF下载量: 64
- 被引次数: 23