2023年土耳其两次强震对周围地区的静态应力影响.pdf

书书书第 卷第 期 年 月防灾科技学院学报 ，年土耳其两次强震对周围地区的静态应力影响靳志同，崔华伟，刘佳璐，冯淦，余海琳（防灾科技学院，河北 三河 ；河北省地震动力学重点实验室，河北 三河 ；山东省地震局，山东 济南 ）摘要：基于 年 月 日土耳其卡赫拉曼马拉 地震和 地震的有限断层破裂模型，计算了两次强震在周围地区产生的位移场、应力场以及同震库伦应力变化对 的余震的应力触发。结果表明：（）两次强震在地表产生了较大的位移，其中，地震产生的最大水平位移约为 ，最大隆升约为 ，最大沉降约为 。地震产生的最大水平位移约为 ，最大隆升约为 ，最大沉降约为 。两次强震后地表的最大水平位移约为 ，最大隆升约为 ，最大沉降约为 。（）地震震中南北两侧为面膨胀区，东西两侧为面压缩区，震中区域的主压应力方向为近南北向，主张应力方向为近东西向。地震震中北东和南西侧为面膨胀区，南东和北西侧为面压缩区，震中区域的主压应力方向为近东西向，主张应力方向为近南北向。（）地震在 地震震中处的库伦应力增加约为 ，说明 地震对 地震的发生是有触发作用的。两次强震的发生使得有限断层周围出现较大范围的应力卸载，库伦应力增加的区域主要出现在两个地震的有限断层的两端，应该注意库伦应力增加区域发生地震的危险性。（）地震对 地震所在断层的南段和北段上 的余震是有促进作用的，深度上的 的余震主要发生 地震所在断层北段与 地震所在断层东段的交叉区域，地震震中西侧的 的余震是由 地震所触发。关键词：土耳其地震；破裂模型；位移场；应力场；库伦应力分布中图分类号：文献标识码：文章编号：（）收稿日期：基金项目：河北省地震动力学重点实验室开放基金（）；河北省高等学校自然科学技术项目（）；国家自 然 科 学 基 金（）；中 央 高 校 科 研 业 务 费 专 项（）；河 北 省 地 震 科 技 星 火 计 划 项 目（、）作者简介：靳志同（），男，博士，副教授，主要从事地震学和地壳应力场方面的研究 引言 年 月 日 时 分（），在土耳其卡赫拉曼马拉发生 地震，震中位于 、。年 月 日 时 分（），在该地震北东方向约 又发生 了 地 震，震 中 位 于 、。截至 年 月 日，土耳其 个省份至少有 人死亡，人受伤，座建筑被毁。在叙利亚西北部，至少有 人死亡，人受伤，座建筑被毁，数千人受损，土耳其戈尔巴西和哈塔伊发生了液化和地面沉降 。土耳其卡赫拉曼马拉 地震是由浅部走滑断 裂 引 起 的，该 事 件 破 裂 为 一 条 走 向 北东 南西的近直立左旋断层，土耳其卡赫拉曼马拉 地震的发震断层为近直立的左旋断层，走向为东西向 。这两次地震地处安纳托利亚板块、阿拉伯板块和非洲板块的交界处附近的东安纳托利亚断裂带或死海转换断裂带上 。死海转换断层是走滑断层，容纳了阿拉伯板块半岛相对于非洲和欧亚板块北向运动，土耳其北部的安纳托利亚断层是右旋走滑结构，容纳了安纳托利亚块相对于欧亚大陆和非洲向西的大部分平移运动 。大量研究表明，通过计算地震引起的库伦应力变化发现，地震之间存在相互作用，这对后续地震 发 生 的 地 点 和 时 间 都 有 一 定 影 响 。等 对 年发生的 地震序列（、）进行了分析，研究发现余震丛集基本上都发生在应力增加区。万永革等 研究发现，中国 个大地震（年邢台地震、年唐山地震、年青海共和地震和 年云南丽江地震）的第一次子破裂事件所产生的静态库伦应力变化均对后续的子破裂事件有触发作用。沈正康等 采用更加接近真实岩石圈状况的黏弹性水平层状介质模型，对 年花石峡 地震、年都兰 地震、年玛尼 和 年玛尼 地震的应力触发进行研究，结果表明东昆仑断裂带上之前发生的地震对其后续地震和余震有触发关系。等 的研究发现，如果已知周围的断层滑动习性，就可以对大的余震的发生地点进行预测。等 研究了 年苏门答腊地震周围断层的库伦应力变化，发现 断层上的库伦应力显著增加，随后便有一个 的大地震发生在该断层上。汶川地震发生后，等 、万永革等 分别给出了 年汶川 地震和于田 地震在周围断层的不同断层段上产生的库伦应力变化，并指出部分断层上的大震的发震被提前或滞后，年 月 日发生在汶川 地震南部的芦山 地震和 年 月 日发生在于田地震震中东北部的 地震就验证了他们的研究成果的合理性。在 年尼泊尔地震后，盛书中等 计算了距尼泊尔首都加德满都约 的尼泊尔 地震对中国大陆地区的静态应力触发情况，发现尼泊尔地震对中国大陆产生的应力变化量值很小，其产生的应力加载主要集中在其邻近的西藏和新疆地区的部分断层上。万永革等 基于破裂模型及均匀弹性半空间模型，分别计算了尼泊尔强震序列（年 月 日至 月 日期间发生的 次 地震）在中国大陆及其附近主要断层上产生的库伦应力变化和水平应力变化，得到中国西藏日喀则地区的定日县 地震、聂拉木县 地 震 是 由 尼 泊 尔 地 震 序 列 触 发 的，年尼泊尔地震序列对喜马拉雅主山前逆冲断裂和青藏高原腹地的大部分拉张构造起到触发作用。计算方法地震是由地下岩石错动产生的，在地震之后的比较短的时间内，弹性能的释放是主要的物理过程，可以把地球中各种固态、液态等介质一阶近似为半无限空间内均匀各向同性的完全弹性体 ，。如果已知断层面的几 何 参 数 和 错 动量，则可严格算出弹性体内部产生的位移场和应变场 。提出了利用正交归一法计算地震应力场格林函数的方法，并在此基础上建立了黏弹松弛分层模型下的地震同震及震后形变的模型，并发展了相应的数值模拟方法 ，该方法同样可以用于成层介质弹性半空间下的静态应力触发。本文使用的破裂模型是基于对远震和大地测量观测的分析而得到（图 ），图 中 、和 分别表示 地震的有限断层破裂模型的 个子断层 ；、和 分别表示 地震的有限断层破裂模型的 个子断层 。考虑到活动断层面上的库伦应力变化受到剪应力、正应力和摩擦系数的影响，库伦应力变化（）可表示为 （）式中，为断层面滑动方向上的剪切应力变化（当 和断层面滑动方向一致时为正）；为断层面上正应力变化（拉伸为正）；为视摩擦系数，由于地下介质可能会包含孔隙流体，那么，如果考虑孔隙流体和断层面上的介质特性，视摩防灾科技学院学报第 卷图 土耳其两次强震的破裂模型 擦系数 通常取 ，本文依据前人的经验 ，取 。表 和 地震所在地壳的速度模型 序号深度 波速度 （）波速度 （）密度 （）如果计算得到该处静态库伦破裂应力变化为正值（负值），则表明该处受之前地震的影响，该处发生地震的时间将提前（推迟），有大量研究表明 ，当应力变化达到 就能有效地影响后续地震的发生，为应力触发的阈值。两次强震产生的位移场为研究土耳其两次地震对周围物质的运动情况造成的影响，基于两次强震的破裂模型 ，和 给出的地壳成层速度模型（表 ），运用 软件包，计算土耳其两次强震在周围地区产生的同震位移场（图 ）。图中黑色箭头表示水平位移场，底色表示垂直位第 期靳志同等：年土耳其两次强震对周围地区的静态应力影响图 土耳其两次强震在地表产生的位移场 移场，从黄色到红色表示隆升，从白色到蓝色表示沉降水平位移取对数后再绘图。地震产生的地表位移场（图 ）具有以下特征：从水平位移场来看，震中周围物质从地震南北两侧向震中汇聚，然后向东西两侧流出，符合 走 滑 型 地 震 产 生 的 位 移 场 的 基 本 特防灾科技学院学报第 卷征 ，最 大 水 平 位 移 约 为 （、）。从垂直位移场来看，震中附近隆升和沉降区域相互交错，沿地震的走向（）所在直线近似为轴对称图形，呈蝴蝶状，具体表现为震中南北两侧由震中向外大致呈现先沉降后隆升，震中东西两侧由震中向外大致呈现先隆升后沉降，最大隆升值约为 ，最大沉降值约为 。地震产生的地表位移场（图 ）具有以下特征：从水平位移场来看，震中周围物质分别从地震北东和南西向震中汇聚，然后向北西和南东向流出，符合走滑型地震产生的位移场的基本特征 ，最 大 水 平 位 移 约 为 （、）。从垂直位移场来看，震中附近隆升和沉降区域相互交错，沿地震的走向（）所在直线近似为轴对称图形，具体表现为震中北东和南西侧由震中向外大致呈现先沉降后隆升，震中北西和南东侧由震中向外大致呈现先隆升后沉降，最大隆升值约为 ，最大沉降值约为 。两次强震共同产生的地表位移场（图 ）具有以下特征：从水平位移场来看，震中周围物质分别从地震北东和南西向震中汇聚，然后向北西和南东向流出，符合走滑型地震产生的位移场的基本特征 ，最大水平位移约为 （、）。从垂直位移场来看，震中附近隆升和沉降区域相互交错，较为复杂，但是整体上还是呈现出震中南北两侧为隆升区域，东西两侧为沉降区域，最大隆升值约为 ，最大沉降值约为 。王永哲等 ，基于均匀弹性半空间位错模型和美国地质调查局震源机制解，分别模拟了两次强震的三维同震形变场，得到 地震产生的最大水平向位移达到 ，垂直向位移达到 （沉降），得到 地震产生的最大水平向位移达到 ，垂直向位移达到 （隆升）。这与本文的结果有差 别：本 文 得到 地震产生的最大水平向位移达到 ，垂直向位移达到 （沉降），得到 地震产生的最大水平向位移达到 ，垂直向位移达到 （隆升）。这种差别主要是因为王永哲等 ，的结果是基于均匀弹性半空间位错模型和美国地质调查局震源机制解，模拟了两次地震的同震形变场，并利用 卫星的 成像几何参数分别计算了此次地震的地表同震位移，而本文使用的两次地震的破裂模型是基于远震和区域地震和大地测量观测，包括远震体波、个强运动加速度计站、个高速率 站和 个静态 偏移的观测结果等分析后给出的破裂模型（该模型更新时间为 ），。两次强震产生的应力场为了直观地反映地震产生的应力变化对附近区域的动力学影响，本节分别讨论两次强震在水平面上产生的应力场。将应力变化张量写成矩阵的形式 ：（）（）槡（）（）槡（）式中，表示最大应力、表示最小应力。水平应力的面应力分量表示为：（）式中，表示水平应力造成的该地区挤压和扩张。最大应力 与 轴的夹角 可以由下式求得：()（）可以看到，如果水平面应力大于零，则表示该地区伸张（膨胀），反之，水平面应力小于零则表示压缩。这样就可以在水平面内把应力变化张量表示出来，并用此来分析地震对周围介质产生的应力变化对该地震附近区域的动力学过程产生的影响 ，下文对应力取对数后再绘图。本节基于土耳其 年 月 日在卡赫拉曼马拉发生的 和 两次强震的破裂模型 ，和 给出的地壳成层速度模型，运用 软件包 ，计算两次强震在周围地区产生的同震应力场（图 ）。图 第 期靳志同等：年土耳其两次强震对周围地区的静态应力影响图 土耳其两次强震在震源深度上产生的应力场 中，绿色箭头和黑色箭头分别表示主张应力（最大主应力）和主压应力（最小主应力），底色表示水平面应力，黄色到红色表示面膨胀区，白色到蓝色表示面压缩区，拉张为正。地震在震源深度（）上的应力场（图 ）具有以下特征：从面应力上来看，震中南北两侧为面膨胀区，膨胀区内主张应力起主要作用，东西两侧为面压缩区，压缩区内主压应力起主要作用，膨胀区和压缩区的分界线与发震断层的两个节面比较一致。从水平主压应力和主张应力来看，震中区域的主压应力为近南北向，主张应力方向为近东西向。土耳其 地震为走滑性质，结合有限断层模型 ，得到断层倾向南西，因此震中东侧为上盘，西侧为下盘。从上文得到的该地震产生的地表的水平位移（图 ）看出，上盘物质从南侧向震中汇聚后又向东侧流出，下盘物质从北侧向震中汇聚后向西侧流出。由各点的位移情况可以得到，地震南侧上盘区域各个点的受力状况应该以张应力为主，而东侧上盘区域各个点的受力状况应该以压应力为主，同时，地震西侧下盘区域各个点的受力状况应该以压应力为主，而北侧下盘区域各个点的受力状况应该以张应力为主。通过各点的运动情况分析得到的应力状态与本节计算的应力场的结果（图 ）是比较一致的。地震在震源深度（）上的应力场（图 ）具有以下特征：从面应力上来看，震中北东和南西侧为面膨胀区，南东和北西侧为面压缩区，膨胀区和压缩区的分界线与发震断层的两个节面比较一致。从水平主压应力和主张应力来看，震中区域的主压应力为近东西向，主张应力方向为近