2021年玛多M_S7.4...山东井水位同震响应特征分析_刘凯.pdf

文章编号:1001-8956(2023)02-0210-07中图分类号:P3157文献标识码:A2021 年玛多 MS74 和 2022 年门源 MS69地震引起的山东井水位同震响应特征分析刘 凯1，陈其峰1，张 军2，孙 豪3，宋 磊4(1山东省地震局聊城地震监测中心站，山东 聊城 252000;2山东省地震局菏泽地震监测中心站，山东 菏泽 274000;3山东省地震局烟台地震监测中心站，山东 烟台 264000;4枣庄市应急管理局，山东 枣庄 277000)摘要:2021 年 5 月 22 日青海玛多 MS74 地震和 2022 年 1 月 8 日青海门源 MS69 地震引起栖霞鲁 07井、枣庄鲁 15 井和菏泽鲁 27 井水位不同程度的同震响应。基于秒数据，对比分析 3 口井水位同震变化形态、幅度与分钟值记录的差异，并从含水层渗透性变化、地震能量密度等方面进行同震响应机理探讨。结果表明，井水位秒数据能够更加完整地记录水震波信息，更加精确展现水位同震变化形态、幅度;井水位对于远场大震的震荡形态主要受含水层水文地质条件的影响。在正常应力背景下，远场大震引起枣庄鲁 15 井水位同震响应的地震能量密度阈值大约是 154104Jm3。关键词:水位;同震响应;秒数据doi:1016256/jissn1001-8956202302012地震孕育、发生过程中产生的静态或动态应力作用于含水层，体现在观测井的水位变化上1，可直接、有效地揭示含水层介质对应力应变的响应过程2，其响应范围和变化幅度是传统的抽水、注水等研究手段所达不到的3。由于井水位同震响应分布的广区域性和响应特征的易识别性4，受到许多学者的关注。目前大多数的研究主要是从一井多震、一震多井两个角度开展5。针对不同的地震，记录同一口井的同震响应变化特征，有助于排除井孔水文地质条件的影响，突出地震因素;记录一次地震所引起的不同区域多口井的同震响应特征，有助于分析不同井孔水文地质条件对地震的响应3。已有的研究结果表明井水位同震响应形态主要为震荡和阶变6。同震响应幅度与井孔条件、含水层的导水系数、储水系数、含水层孔隙度7 和地震波类型、周期、振幅8 等因素有关，其受地震面波影响最为明显9。井水位同震变化的观测记录取决于水位仪的采样率。频率越高，越能真实还原水位变化。地震波和水震波都属于高频信号，周期远小于 1 min，相较于分钟值和时值，秒数据具有绝对优势，能更清晰完整地记录到水震波10。本文中利用山东省地下流体观测台网中的栖霞鲁 07 井、枣庄鲁 15 井和菏泽鲁 27 井的水位秒数据(水位均指水位埋深)，结合 2021 年 5 月 22 日青海玛多 MS74 地震和 2022 年 1 月 8 日青海门源 MS69 地震，进行井水位的同震响应特征分析。第37卷第2期2023年6月内陆地震INLANDEATHQUAKEVol37No2Jun2023收稿日期:2022-07-14;修回日期:2022-09-05作者简介:刘凯(1988)，男，工程师，2009 年毕业于山东师范大学地理科学专业，主要从事地震地下流体观测工作 E-mail:LK8822 126com1同震响应变化栖霞鲁 07 井、枣庄鲁 15 井和菏泽鲁 27 井是山东省地下流体观测台网中的骨干观测点，三口井水位日变化动态清晰，含水层承压性好，均有明显的潮汐变化11。栖霞鲁 07 井井深600 m，套管深度 7355 m，主要观测段为 7355 600 m 的花岗岩裂隙承压水，井孔水温140，pH 值 836。枣庄鲁 15 井井深约 501 m，下套管 23962 m，揭露含水层为砂岩孔隙裂隙承压水。菏泽鲁 27 井井深 2 000 m，井孔上层为松散砂土黏土沉积层，1 138 m 以下为基岩，观测水层为深层奥陶系裂隙岩溶水(表 1)。表 1观测井基本信息序号井孔名称井深/m观测含水层岩性断裂带观测仪器型号投入观测时间(年-月-日)1栖霞鲁 07 井600花岗岩刘家亭断裂DSC-2018-07-182枣庄鲁 15 井501砂岩韩庄断裂SWY-2017-12-073菏泽鲁 27 井2 000灰岩聊考断裂带DSW-2018-03-032021 年 5 月 22 日 2 时 4 分，青海果洛州玛多县发生 MS74 地震，震中位于 3459N，9834E。2022 年 1 月 8 日 1 时 45 分，青海海北州门源县发生 MS69 地震，震中位于 3777N，10126E。对于这 2 次地震，栖霞鲁 07 井、枣庄鲁 15 井和菏泽鲁 27 井的水位秒数据记录值和分钟记录值呈现出不同的同震响应变化。玛多 MS74 地震引起的栖霞鲁 07 井水位秒数据呈现震荡形态，向下波动最大幅度 165 cm，向上波动最大幅度 95 cm，向下幅度是向上幅度的174 倍;水位分钟值呈现突降形态，下降幅度 17 cm，向上波动最大幅度 08 cm，向下幅度是向上幅度的 213 倍。枣庄鲁 15 井水位秒数据呈现震荡形态，向下波动最大幅度 21 cm，向上波动最大幅度 15 cm，向下幅度是向上幅度的 14 倍;水位分钟值呈现突降形态，下降幅度 06cm，向上无明显变化。菏泽鲁27 井水位秒数据呈现震荡形态，向下波动最大幅度49 cm，向上波动最大幅度 58 cm，向下幅度是向上幅度的 084 倍;水位分钟值也呈现震荡形态，下降幅度18 cm，上升幅度 09 cm，向下幅度是向上幅度的 20 倍(表 2 和图 1)。表 2玛多 MS74 地震引起的井水位同震变化参数观测井井震距/km水位分钟值/m均值最高点最低点水位秒数据/m均值最高点最低点栖霞鲁 07 井2 055696869606985696868737133枣庄鲁 15 井1 742453094530945315453094529445330菏泽鲁 27 井1 547605936058460611605936053560642门源 MS69 地震引起的栖霞鲁 07 井水位秒数据呈现震荡形态，向下波动最大幅度78 cm，向上波动最大幅度 47 cm，向下幅度是向上幅度的 166 倍;水位分钟值呈现突降形态，突降幅度59 cm，突升幅度03 cm，向下幅度是向上幅度的1967 倍。枣庄鲁15 井水位秒数据1122 期刘 凯等:2021 年玛多 MS74 和 2022 年门源 MS69 地震引起的山东井水位同震响应特征分析波动频率增大，但是变幅不明显;水位分钟值无明显响应变化。菏泽鲁 27 井水位秒数据呈现震荡形态，向下波动最大幅度 06 cm，向上波动最大幅度 06 cm，上下波动对称;水位分钟值也是呈现震荡形态，突降幅度 03 cm，突升幅度 02 cm，向下幅度是向上幅度的 15 倍(表 3 和图2)。图 1玛多 MS74 地震引起的井水位同震变化曲线图(a)栖霞鲁 07 井水位秒数据(b)栖霞鲁 07 井水位分钟值(c)枣庄鲁 15 井水位秒数据(d)枣庄鲁 15 井水位分钟值(e)菏泽鲁 27 井水位秒数据(f)菏泽鲁 27 井水位分钟值Fig1Coseismic variation curves of well water level caused by Maduo MS74 Earthquake表 3门源 MS69 地震引起的井水位同震变化参数观测井井震距/km水位分钟值/m均值最高点最低点水位秒数据/m均值最高点最低点栖霞鲁 07 井1 733533553325394533652895414枣庄鲁 15 井1 48143832438324382743831菏泽鲁 27 井1 281613016129961304613026129661308212内陆地震37 卷图 2门源 MS69 地震引起的井水位同震变化曲线图(a)栖霞鲁 07 井水位秒数据(b)栖霞鲁 07 井水位分钟值(c)枣庄鲁 15 井水位秒数据(d)枣庄鲁 15 井水位分钟值(e)菏泽鲁 27 井水位秒数据(f)菏泽鲁 27 井水位分钟值Fig2Coseismic variation curve of well water level caused by Menyuan MS69 Earthquake2结果与分析远场大震的地震波作用于井含水层，在动态应力作用下，井含水层系统原有的水位动态平衡被打破，含水层介质孔隙压力、渗透系数及层间连通性的改变都有可能引起井水位的变化12。一般来说，同一口井对不同远震的同震响应形态一致13，不受震源机制与震中方位的影响。同震响应幅度与震级、井震距存在一定比例关系14。当响应过程发生改变时，有可能预示着含水层所处区域的应力状态发生变化15。21变化形态远场大震地震波产生的动态应力作用于井含水层系统，含水层介质发生有规律的弹性变化，使得井水在井孔含水层之间流动，出现水位变化。震后随着地震波能量的衰减，水位慢慢恢复至震前水平16。水位震荡是含水层对地震波的弹性响应。玛多 MS74 地震和门源MS69 地震引起的三口井水位秒数据同震变化均为震荡形态。与分钟值相比，秒数据记录的水震波震荡形态更加对称，波形更完整。对比分析栖霞鲁 07 井和菏泽鲁 27 井水位秒数据变化形态可以发现，2 次地震引起的菏3122 期刘 凯等:2021 年玛多 MS74 和 2022 年门源 MS69 地震引起的山东井水位同震响应特征分析泽鲁 27 井水位上下波动幅度比较接近，说明含水层进行了有规律的弹性压缩与拉张;而栖霞鲁 07 井均是突降幅度明显大于上升幅度。井水位的上升或下降，是井孔与含水层交互作用的结果17。井含水层系统既有水平向的井孔含水层间的水流交换，也有垂向的含水层弱透水层间的水流交换18。在地震波动态应力作用下，含水层交替发生压缩与拉张19。在极短时间内，栖霞鲁 07 井含水层受到拉张作用，介质孔隙增大，渗透性变强，水从井孔净流向含水层，再交替受到压应力作用，含水层受到压缩，介质孔隙减小，渗透性降低，水从含水层回流向井孔的速度减弱，导致水位上升幅度变小。在地震波的动态应力消失之后，含水层水文地质条件恢复，井水位便恢复至震前水平。因此，两口井水位呈现出的不同震荡形态与当地的水文地质条件有关，远场地震波只是起到触发作用20。22变化幅度远场大震引起的井水位同震响应幅度与震中距、震级有一定关系。震中距越小、震级越大，引起的井水位同震响应幅度越大21。部分学者据此提出地震能量密度参数22，表示地震波在传播过程中作用在单位体积地层介质上的最大能量值，与地震震级和震中距密切相关23。一般来说，对于同一口井，在相对稳定的井孔结构和含水层水文地质条件的前提下，远场大震能够引起井水位同震响应的地震能量密度存在一个阈值24。当地震能量密度高于阈值时，井水位产生同震响应;低于阈值时，井水位不受影响。如果地震能量密度小于阈值，却能引起井水位同震响应，在一定程度上有可能表明观测井含水层受张应力影响较大，含水层介质孔隙变大，水流渗透性增强，较小的地震能量密度就能引起井水在井含水层之间流动。当含水层受到压应力作用时，含水层介质孔隙率变小，渗透性降低，需要较大的地震能量密度才能引起井水位的同震响应。因此，井水位的记震能力可能会出现地震能量密度超过阈值无法激发同震响应，或者低于阈值时却能记录到水位同震响应的情况，即引起同震响应的地震能量密度阈值发生变化。井水位记震能力的变化在一定程度上可以反应区域应力场变化，对于研究地震的孕育和发展起到一定的指示作用。枣庄鲁 15 井水位秒数据对于玛多 MS74 地震，有明显的同震变化;但是对于门源 MS69地震，只是水位震荡频率增大，波动幅度基本不变。结合近几年井孔周围 100 km 范围内无MS20 地震发生，区域应力场相对稳定，基本可以推断门源 MS69 地震产生的地震能量密度，非常接近引起枣庄鲁 15 井水位同震响应的阈值。对于震级 M，距离为 r 处的地震能量密度 e 可以用经验关系式表示22，logr=048M033loge14(1)式中:r 为井震距，单位 km;e 为地震能