-
目前使用地震波资料求取地壳构造应力场的方法主要有两种:一是使用地震的震源机制解来反演应力场;二是使用P波初动符号计算综合震源机制解的方法获得应力场方向[1]。使用P波初动求取单个地震震源机制解需要较多的P波初动信息,并且记录地震波的台站需要合理的覆盖,否则计算出的震源机制解会有较大的误差。再者有许多小微震的震源机制解由于数据原因并不能准确求取,这样通过震源机制解求解应力场的数据量就会大大减少,降低了解的可靠性。而通过P波初动计算综合震源机制解方法则可以利用众多小微震的P波初动信息,增大了输入的信息量,提高了结果的置信度。
利用大量小震P波初动符号推断应力场的方法最早是由Aki[2]提出的,该方法充分运用大量不能确定震源机制解的小微震的P波初动信息,因此获得较广泛的应用。李钦祖等[3]运用该法对单个台站的小震资料分别求取了红山台和沙城台的区域应力场方向,并给出了其置信范围;许忠淮等[4]使用该方法运用多个台站和多个地震资料求取了区域应力场的方向,并通过合成数据的方法进行验证,确定了该方法能得到稳定的应力场;薛宏运等[5]运用该方法求解了鄂尔多斯周围地区分区现代地壳应力场并进行了分析;汪素云等[6]运用该法求取了中国东部大陆分区的地震构造应力场并进行了分析;万永革等[7]运用不同地震距对网格点应力场产生影响的P波综合震源机制法获得了台湾地区三维应力结构的初步结果;卜玉菲等[1]运用P波数据反演了甘东南地区构造应力场,并分析了成因机制;盛书中等[8]运用该法求解了鄂尔多斯周缘地区的应力场方向。
本文所研究的鲁中地区指山东地区的中部,在构造上包括济阳坳陷、鲁西隆起和郯庐断裂带的山东段。在以往的研究中,专门针对鲁中地区构造应力场的研究较少,更多的是包括鲁中地区的更大范围的应力场研究。涉及鲁中地区的应力场数据可能只有几个,是鲁中地区应力场的综合表现。但该区构造复杂,构造应力场可能不是统一的状态,以往研究无法展现该区分辨率更大的应力场数据。本文针对目前鲁中地区构造应力场的研究现状,基于鲁中地区的近震数据资料,利用万永革等[7]发展的综合震源机制解法计算得到了鲁中地区较为精细的应力场数据,弥补了以往大区域应力场研究对局部应力场研究的不足,揭示了研究区小区域应力场的细部特征,以期为该区构造应力场及相关方面的研究提供参考。
-
本研究确定区域应力场方向是按照万永革等[7]发展的方法,综合了P波初动分配到每个网格点的数量以及研究方法的适用性。为保证计算得到的相邻各网格点的应力场数据具有较好的连续性,将鲁中地区(116°~120°E,34°~39°N)划分成0.25°×0.25°的网格。由于记录地震的各台站位置距离网格点的远近不一样,所以给各台站记录到的地震P波初动信息根据距离大小赋予权重,计算方位使用Shen等[9]的大地测量数据计算应变公式
$$W = {e^{ - {r^2}/{D^2}}}{\text{。}}$$ (1) 式中:D为距离衰减常数;r为折合距离,其计算参考公式
$$ r=\sqrt{{\left[111.199\times \left(x-{x}_{0}\right)\cos\frac{(y+{y}_{0})}{2}\right]}^{2}+ {\left[111.199\times (y-{y}_{0})\right]}^{2}+{\left(z-{z}_{0}\right)}^{2}}{\text{。}} $$ (2) 式中:x、y、z为地震震中的经纬度和震源深度;x0、y0、z0为网格点的经纬度和深度。其中,网格点处的深度取水平面以下10 km,即计算10 km深度处区域的应力场方向。如果D值较小,则计算中将使用距离网格点较近的P波初动数据,各网格点的应力场方向较为独立;如果D值较大,则计算中就可以使用网格点周边网格内的P波初动数据,对各网格点的数据具有平滑作用。参考前人的研究经验,本研究取D值为50 km。通过计算可知,当地震台站距离网格点50 km时,台站记录P波初动信息的权重为36.8%;当距离网格点100 km时,P波初动权重为1.83%。因此,我们选取距离网格点的P波初动的最大距离为100 km。
我们采取1°×1°×1°的间隔搜索P、B、T轴的方位,并选择矛盾比(与综合震源机制解不符的P波初动数和总P波初动数的比值)最小的综合机制解为最优解,这样做是为了确定与输入的P波初动符号拟合最好的综合震源机制解,然后分析P、T轴的方位,得到的最优解可视为该点附近主应力方位。只要网格点处计算所运用的P波数据足够多,并且分布方位覆盖较好,则可认为运用多个地震的P波初动求取的应力场方向是最优可靠的。
-
双差地震定位[10]是一种相对定位方法,该方法利用信号的走时差反演震源位置,能够有效地消除震源至台站共同传播路径效应,对地壳速度模型的依赖性相对较小,是一种定位精度较高的重定位方法。本文利用1 402个原始地震数据,通过双差定位法获得了820个地震的重定位数据(图1),双差定位误差如图2所示。由图2可以看到,经度方向的平均误差为0.289 7 km,多数地震的误差在0.7 km以内;纬度方向的平均误差为0.278 1 km,多数地震的误差在0.6 km以内;震源深度的平均误差为0.371 2 km,多数地震的误差在1 km以内。
使用双差定位程序对研究区内地震进行了重新定位,令地震震源位置的经纬度数据更加精确,震源深度也更加合理。本次研究使用了2008年3月15日到2019年8月23日的重新定位后的820个地震(震级范围ML1.0~4.3)的5 048个P波和Pn波的初动数据,初动符号来自于山东地震台网的观测报告。
-
崔鑫等[11]基于Hyposat和HypoDD的定位方法,根据残差结果极小值获得了山东地区最优的P波速度模型,并通过了多重验证。本文计算P波射线方位角和离源角所使用的地壳速度模型采用崔鑫等[11]的研究成果(表1)。
表 1 鲁中地区地壳速度模型
层顶埋深/km P波速度/(km/s) 0 6.13 22 6.88 33 7.93 -
在应用综合震源机制解法求解应力场的过程中,将应力的相对大小设置为固定值0.5,共反演了215个网格点处的应力场数据,扫描得到的最优解为一个或两个,最大矛盾比不超过0.3。
本研究通过计算得到了研究区较为精细的应力场结果,该结果能更好地体现鲁中地区应力场方向的变化趋势。主压应力轴方向(图3)在空间上呈现明显不均匀性,总体上倾角较小且方向上呈现逆时针旋转特性,从研究区北部的NWW-SEE向往南逐渐变为NEE-SWW向,这种趋势在研究区的东部更为明显。整体上看主压应力轴从东向西呈现放射状分布。在研究区的东部及大部主张应力轴方位(图4)由北往南从NNE-SSW逐渐过渡到NNW-SSE,过渡的连续性较好。研究区西南部主张应力轴呈NEE-SWW方位分布,与周边主张应力轴方位有较大差别,方位一致性较差,并且该区主张应力轴的倾角较大。
研究参照世界应力图的划分原则[12],根据震源机制解3个应力轴倾角的大小,将震源机制类型(即其反映的应力状态)分为6种(表2):正断型(NF)、正走滑型(NS)、走滑型(SS)、逆走滑型(TS)、逆断型(TF)和不确定型(U)。根据分类标准,将不同类型的震源机制用不同的颜色表示(图5)。由图4可知,研究区内的应力状态主要为走滑型和逆断型,还有少量的正断型、正走滑型和不确定型。
表 2 震源机制解类型划分表
类型 P 轴倾角 B 轴倾角 T 轴倾角 NF ≥52° <35° NS 40°≤倾角<52° <20° SS <40° >45° <20° <20° >45° <40° TS <20° 40°≤倾角<52° TF <35° ≥52° U 上述类型之外的震源机制解 -
山东地区位于华北地区东部,根据谢富仁等[13]对中国大陆的应力场分区,山东地区属于一级应力场分区中的中国大陆东部应力区,位于二级应力场分区的东北-华北应力区和三级应力场分区的华北应力区、四级应力场分区的华北平原应力区的南部。对于山东地区的应力场,前人已经做了较多研究,但是对于局部应力场的细化研究还是较少。张玲等[14]运用格点尝试法计算了山东地区的应力场,得到山东地区的现代应力场总体特征为NEE(80°)方向主压和NNW向(350°)主张;郑建常等[15]使用小震的震源机制解反演了山东地区的背景应力场,得到了P轴的优势方位为NEE-SWW方向,T轴的优势方位为NNW-SSE方向;李鹏等[16]运用地应力实测数据计算了山东地区的应力场方向,其中最大水平主应力优势方位整体表现为NWW-SEE向,但NEE-SWW向也占有较大比例。另外,根据李钦祖等[17]、许忠淮等[4]、刘丽等[18]对华北地区应力场的研究结果,均认为华北地区的主压应力轴方位在NEE-SWW向附近,而主张应力轴方位在NNW-SSE向附近。由于本文的研究区域较小、数据量相对丰富,且网格划分也较小,因此得到的应力场结果更为精细,并且填补了以往研究没有覆盖的区域,与前文介绍的多位研究者得到的山东地区、华北地区的应力场有部分区域一致性较好,但也有区域一致性较差。
-
研究区内广饶-齐河断裂(F4)以北为济阳坳陷,该区主压应力轴的总体方位为NWW-SEE向,主要受到西北侧的鄂尔多斯地块SEE方向的推挤[15],或是太行山地块SEE方向的推挤作用造成[19],而主张应力轴的总体方位为NNE-SSW向,解秋红[20]对该区应力场的有限元模拟结果与本文基本一致,但与郑建常等[15]、许忠淮等[4]在该区主压应力轴NEE向、主张应力轴NWW向的研究结果不一致,考虑可能是由于该区震源机制解个数较少并且应力场取了大范围的平均值导致,而本文则是较为精细化的计算结果。区内近NS向如白桥断裂、惠东断裂、临邑断裂以及由较小断裂组成的断裂带均为走滑断层,而区内近EW向(包括NEE向和NWW向)的以陈南断裂、义南断裂、埕南断裂等十多条断裂为代表的断裂为正断层断裂,这些正断层断裂与走滑断裂存在调节、派生、继承的关系[21],其性质与本文的应力场结果有较好的一致性,而且对该区的单震震源机制解主应力轴的统计分析结果也显示与本文的应力场结果有较好的吻合[20]。该区综合震源机制解类型为走滑型,这可能与该区东西两侧的两大断裂系(郯庐断裂和兰聊-盐山断裂)的右旋走滑机制有关。
-
研究区内广饶-齐河断裂(F4)以南、郯庐断裂带以西为鲁西隆起,该区主压应力轴的总体方位为NEE-SWW向,与郑建常等[15]、许忠淮等[4]、刘峡等[22]的研究结果一致,但也有地方不同。如该区西北边缘的主压应力轴仍为NWW向,可能受西北侧推挤力影响依然较大,另外在该区南部边缘主压应力轴显示为NWW-SEE向,具体原因有待进一步研究。总体上鲁西隆起区主要受到秦岭大别构造带NE向的挤压力[15]。该区主张应力轴方位主要为NWW-SSE向,与郑建常等[15],许忠淮等[4]、刘峡等[22]在该区的研究结果一致,但在F23、F24、F25断裂包围的区域内主张应力轴方位为NEE-SWW向,由于该区域较小,暂没有与之一致的研究成果。
在图5中,济南、淄博及周边区域的综合震源机制解类型为正走滑型和不确定型,而在不确定型中则又有多数含有正断型分量,因此总体上该区域的应力类型为拉张型兼有走滑机制。该区域内F5泰山西麓断裂为正断型,F6文祖断裂为正断型,F7禹王山断裂、F9王母山断裂、F11上五井断裂为右旋正断型,F8张店-仁河断裂、F10益都断裂为左旋正断型。以上断裂性质与该区的应力场类型一致性较好。
在研究区的西南部(除由F23、F24、F25断裂包围的区域)综合震源机制解类型为逆断型,更为详细的计算数据显示,该区域应力状态中也含有一定的走滑分量(由于本文采用的应力类型划分标准可能不够细致,因此该区应力类型划为了逆断型)。从图3、图4的主压应力轴和主张应力轴方位也可看出,这与鲁西地区的地震呈走滑或近走滑型应力状态有一定的相近[23]。由F23、F24、F25断裂包围的区域综合震源机制解类型也为逆断型,但近似纯逆断型,这可能是由于该区及周边复杂的地质构造和应力场造成的。该区主张应力轴倾角较大(大于52°),这与崔效锋等[24]运用格点尝试法得到的鲁西南地区主张应力轴的平均倾角7°存在较大差别。但济宁西北侧的赵楼井田的原岩地应力实测结果和布格重力异常分布均显示该区应力场状态为扭压状态[25],大地测量的结果显示鲁中、南地区的应变场存在压缩与膨胀相间的情况[26]。吕朋菊等[27]认为鲁西地区晚第三纪受到NNW-SSE向的挤压,使得NW向断裂以右行压扭为主、NEE向断裂以左行压扭为主,第四纪以来该区仍然受到晚第三纪挤压应力背景的影响。总之,由F23、F24、F25断裂包围的区域综合震源机制解类型为逆断型的计算结果有支持的证据,也有不支持的证据,因此该区域的应力场状态有待进一步的研究确定。
总体来看,鲁西隆起区内的应力状态均含有走滑分量,包括走滑型、正走滑型、含走滑分量的逆断型,这与该区东西两侧的两大断裂系(郯庐断裂和兰聊-盐山断裂)的右旋走滑机制有关。
-
研究区内郯庐断裂带南北两段的应力轴方位分别与上文中与其邻近的鲁西隆起和济阳坳陷的情况基本一致。由图5知,该区大部分区域的综合震源机制解类型为走滑型,这与郯庐断裂带右旋走滑机制相符,但在该区潍坊以东附近区域呈正断型。该处附近的钻孔地应力测量数据显示,在距地面约372.5 m以下的应力结构为正断型,这是安丘–莒县断裂的北段区域部分继承了早第三纪的伸展活动方式所致[28]。
-
基于鲁中地区的近震数据资料,利用万永革等发展的综合震源机制解法,计算了鲁中地区的应力场。计算结果在与先前研究结果保持一定一致性的情况下,也揭示了局部应力场的细部特征。P轴总体上倾角较小且方向上呈现逆时针旋转特性,从研究区北部的NWW-SEE向往南逐渐变为NEE-SWW向;而T轴方位由北往南从NNE-SSW逐渐过渡到NNW-SSE,并有部分区域T轴方位一致性较差且倾角较大。该区应力类型有走滑型、正走滑型、正断型和逆断型(部分含走滑分量)。总体上看,研究区内应力类型都含有走滑分量,这表明郯庐断裂和兰聊-盐山断裂的右旋走滑机制对研究区影响显著,个别地区的应力场和断裂还有继承第三纪应力状态和活动方式的情况。
致谢 万永革教授在应力场计算方面的帮助和指导,本文图件使用Generic Mapping Tools (GMT)[29]绘制。
Tectonic Stress Field in Central Area of Shandong Province with Composite Focal Mechanism Method
-
摘要: 利用双差定位法对鲁中地区2008年3月15日至2019年8月23日发生的1 402个地震进行重定位,获得了820个地震的精确位置,重定位后地震位置空间误差显著减小;采用重定位地震的5 048个P(Pn)波初动符号,运用综合震源机制解法计算了鲁中地区0.25°×0.25°网格的构造应力场,结果显示:P轴总体上倾角较小且方向上呈现逆时针旋转特性,从研究区北部的NWW-SEE向往南逐渐变为NEE-SWW向;而T轴方位由北往南从NNE-SSW逐渐过渡到NNW-SSE,并有部分区域T轴方位一致性较差且倾角较大;最后分别对研究区内济阳坳陷、鲁西隆起、郯庐断裂带的应力场特征进行了分析。Abstract: Using the double difference positioning method, 1 402 earthquakes occurred on March 15, 2008 solstice and August 23, 2019 in the middle Shandong area are relocated, and the exact locations of 820 earthquakes were obtained. The tectonic stress field of 0.25°×0.25° grid in Luzhong area was calculated by using the composite focal mechanism method on basis of 5 048 P (Pn) wave initial motion symbols. The results show that the overall dip of the P-axis is small and it’s direction shows counterclockwise rotation characteristics: from north to south it gradually changes to NEE-SWW direction, while the T-axis direction gradually transitioned from NNE-SSW to NNW-SSE and the T-axis azimuth consistency in some areas is poor and the dip is large. Finally, the stress field characteristics of Jiyang Depression, Western Shandong Rise and Tan-Lu fault zone are analyzed.
-
图 3 P轴方位和网格点使用P波初动符号数分布图
注:F1 埕子口断裂;F2 义南断裂;F3 临邑断裂;F4 广饶-齐河断裂;F5 泰山西麓断裂;F6 文祖断裂;F7 禹王山断裂;F8张店-仁河断裂;F9 王母山断裂;F10 益都-无棣断裂;F11 上五井断裂;F12 鄌郚-葛沟断裂;F13 安丘-莒县断裂;F14 肥城断裂;F15 铜冶店-孙祖断裂;F16 沂源-沂水断裂;F17 沂水-汤头断裂;F18 昌邑-大店断裂;F19 山相家-郝戈庄断裂;F20 日照断裂;F21 新泰-蒙阴断裂;F22 巨野断裂;F23 嘉祥断裂;F24 郓城断裂;F25 苍尼断裂;F26 蒙山山前断裂;F27 峄城断裂
表 1 鲁中地区地壳速度模型
层顶埋深/km P波速度/(km/s) 0 6.13 22 6.88 33 7.93 表 2 震源机制解类型划分表
类型 P 轴倾角 B 轴倾角 T 轴倾角 NF ≥52° <35° NS 40°≤倾角<52° <20° SS <40° >45° <20° <20° >45° <40° TS <20° 40°≤倾角<52° TF <35° ≥52° U 上述类型之外的震源机制解 -
[1] 卜玉菲, 张元生, 万永革, 等. P波极性揭示的甘东南地区构造应力场特征[J]. 地震工程学报, 2013, 35(1): 160-165. doi: 10.3969/j.issn.1000-0844.2013.01.0160 [2] Aki K. Earthquake generating stress in Japan for the years 1961 to 1963 obtained by smoothing the first motion radiation patterns[J]. Bulletin of Entomological Research, 1966, 44(2): 447-471. [3] 李钦祖, 王泽皋, 贾云年, 等. 由单台小地震资料所得两个区域的应力场[J]. 地球物理学报, 1973, 16(1): 49-61. [4] 许忠淮, 阎明, 赵仲和. 由多个小地震推断的华北地区构造应力场的方向[J]. 地震学报, 1983, 5(3): 268-279. [5] 薛宏运, 鄢家全. 鄂尔多斯地块周围的现代地壳应力场[J]. 地球物理学报, 1984, 27(2): 144-152. doi: 10.3321/j.issn:0001-5733.1984.02.004 [6] 汪素云, 许忠淮. 中国东部大陆的地震构造应力场[J]. 地震学报, 1985, 7(1): 17-32. [7] 万永革, 吴逸民, 盛书中, 等. P波极性数据所揭示的台湾地区三维应力结构的初步结果[J]. 地球物理学报, 2011, 54(11): 2809-2818. doi: 10.3969/j.issn.0001-5733.2011.11.011 [8] 盛书中, 万永革, 黄骥超, 等. 应用综合震源机制解法推断鄂尔多斯块体周缘现今地壳应力场的初步结果[J]. 地球物理学报, 2015, 58(2): 436-452. doi: 10.6038/cjg20150208 [9] Shen Z K, Jackson D D, Ge B X. Crustal deformation across and beyond the Los Angeles basin from geodetic measurements[J]. Journal of Geophysical Research: Solid Earth, 1996, 101(B12): 27957-27980. doi: 10.1029/96JB02544 [10] Waldhauser F, Ellsworth W L. A double-difference earthquake location algorithm: method and application to the northern Hayward Fault, California[J]. Bulletin of the Seismological Society of America, 2000, 90(6): 1353-1368. doi: 10.1785/0120000006 [11] 崔鑫, 李铂, 朱元清, 等. 基于Hyposat和HypoDD的山东地区地震波一维速度模型[J]. 吉林大学学报: 地球科学版, 2017, 47(5): 1572-1581. [12] Zoback M L. First-and second-order patterns of stress in the lithosphere: the world stress map project[J]. Journal of Geophysical Research: Solid Earth, 1992, 97(B8): 11703-11728. [13] 谢富仁, 崔效锋, 赵建涛, 等. 中国大陆及邻区现代构造应力场分区[J]. 地球物理学报, 2004, 47(4): 654-662. doi: 10.3321/j.issn:0001-5733.2004.04.016 [14] 张玲, 周翠英, 王锋吉, 等. 山东地区各分区地震应力场特征分析[J]. 华北地震科学, 2004, 22(4): 12-15. doi: 10.3969/j.issn.1003-1375.2004.04.004 [15] 郑建常, 王鹏, 李冬梅, 等. 使用小震震源机制解研究山东地区背景应力场[J]. 地震学报, 2013, 35(6): 773-784. doi: 10.3969/j.issn.0253-3782.2013.06.001 [16] 李鹏, 郭奇峰, 刘洪涛, 等. 山东地区现今地应力场特征与应力积累水平分析[J]. 岩石力学与工程学报, 2017, 36(9): 2220-2231. [17] 李钦祖, 靳雅敏, 于新昌. 华北地区的震源机制与地壳应力场[J]. 地震学报, 1982, 4(1): 55-61. [18] 刘丽, 刘檀, 董一兵, 等. 由震源机制解分析华北地区构造应力场[J]. 华北地震科学, 2017, 35(4): 38-43. doi: 10.3969/j.issn.1003-1375.2017.04.007 [19] 范俊喜. 鄂尔多斯地块运动特征研究[D]. 北京: 中国地震局地质研究所, 2002. [20] 解秋红. 渤海及邻区地应力场分析及构造稳定性评价[D]. 青岛: 中国海洋大学, 2011. [21] 赵利. 鲁西隆起—济阳坳陷晚中生代以来断层伸展与走滑关系研究[D]. 青岛: 中国石油大学(华东), 2015. [22] 刘峡, 傅容珊, 杨国华, 等. 用GPS资料研究华北地区形变场和构造应力场[J]. 大地测量与地球动力学, 2006, 26(3): 33-39. [23] 周翠英, 华爱军, 蒋海昆, 等. 以格点尝试法求取的山东地区现代中小地震震源机制解[J]. 东北地震研究, 2003, 19(1): 1-11. [24] 崔效锋, 谢富仁, 李瑞莎, 等. 华北地区构造应力场非均匀特征与煤田深部应力状态[J]. 岩石力学与工程学报, 2010, 29(S1): 2755-2761. [25] 姜春露, 姜振泉, 杨伟峰, 等. 赵楼井田地应力特征及地质构造形成机制[J]. 煤炭学报, 2011, 36(4): 583-587. [26] 赵金瑞. 鲁中、南应变场和地壳的最新运动[J]. 地震学报, 1981, 3(2): 126-134. [27] 吕朋菊, 朱兴珊. 鲁西中、新生代构造应力场的更迭[J]. 山东矿业学院学报, 1989, 8(4): 18-25. [28] 丰成君, 张鹏, 戚帮申, 等. 郯庐断裂带附近地壳浅层现今构造应力场[J]. 现代地质, 2017, 31(1): 46-70. doi: 10.3969/j.issn.1000-8527.2017.01.005 [29] Wessel P, Luis F, Uieda L, et al. The generic mapping tools version 6[J]. Geochemistry, Geophysics, Geosystems, 2019, 20(11): 5556-5564. doi: 10.1029/2019GC008515 -