井水位的同震响应，是指地震引起的井水位振荡或阶变现象，其动态变化能够较直接地反映地震作用于井—含水层系统的过程，对其进行科学研究是揭示地壳介质对应力—应变过程响应直接、有效的手段之一（中国地震局监测预报司，2005）。近几十年来，许多研究者基于水位分钟值数据从不同方面对井水位的同震响应特征和机理进行了研究（Cooper et al，1965；汪成民等，1988；张子广等，1998；刘耀炜等，2005；车用太等，2006；付虹等，2007；杨竹转等，2010）。但是，基于水位秒钟值数据的同震响应研究相对较少。近年来，向阳等（2017）利用水位秒钟值数据计算了单口井水震波的时频特征，并与地震波时频特征进行对比，得到了井水位观测数据所反映的地震波优势周期，计算了含水层水文参数的变化来解释单口井的同震响应机理。对于多口井水位秒钟值数据同震响应变化的研究则少见报道。

北京时间2021年5月21日21:48:34，云南大理漾濞县发生M_S 6.4地震（25.67°N，99.87°E），震源深度8 km。地震引起云南地区部分数字化井水位观测数据产生同震响应。本文收集整理此次地震云南现有在库正常观测的井水位观测数据及其同震响应记录情况，分析了井水位同震响应空间分布和同震水位变化形态特征，从部分观测井水位分钟值与秒钟值数据所反映的同震响应变化形态分析了观测仪器采样率对井水位映震能力的影响，并对水位同震响应机理进行了初步分析。

通过查询云南地球物理台网“十五”数据库和“十项措施”数据库发现，2021年以来，云南地下流体数字化水位观测井网共有57口井（静水位36口井，动水位21口井）正常在网观测，其空间分布如图 1所示。水位观测仪器主要有LN-3A型数字水位观测仪、SWY-Ⅱ型数字化水位观测仪、ZKGD-3000型水位水温综合观测仪，采样率为1次/min，入库保存为分钟值数据。SWY-Ⅱ型数字化水位观测仪还可以进行秒采样，会自动生成秒钟值数据文件，但秒钟值数据保存时间较短且不能直接入库使用，需要手动下载保存和转换格式后才能使用，非常不便。因此，选用水位分钟值观测数据对漾濞M_S 6.4地震引起云南数字化井水位同震响应特征进行分析。

图 1(Fig.1)

图 1 云南省数字化水位观测井分布 F1为小江断裂带；F9为澜沧江断裂带；F10为怒江断裂带；F11为大盈江断裂带；F12为龙陵—瑞丽断裂带；F15为无量山—营盘山断裂；F17为红河断裂带；F18为楚雄—建水断裂 Fig.1 The distribution of digital water level observation wells in Yunnan province

2021年5月21日漾濞M_S 6.4地震发生前后，云南有57口井的水位观测仪器工作正常，观测数据连续可靠，水位动态变化背景相对稳定。通过对这57口井水位分钟值观测数据进行分析，发现漾濞M_S 6.4地震时有13口井记录到明显的同震响应变化，其水位同震响应形态清晰可辨，变化幅度均大于0.005 m。经统计，记录到同震响应的观测井占总观测井的22.8%，将其水位同震响应变化特征整理归纳列于表 1，同震响应观测井分布如图 1所示。

表 1(Table 1

表 1 漾濞MS 6.4地震引起云南数字化井水位同震响应特征Table 1 Co-seismic response characteristics of digital well water level in Yunnan Province caused by Yangbi MS 6.4 earthquake 井点名称 观测类型 井震距/km 同震响应形态 变化幅度/m 下关井 静水位 39 阶升 0.243 大理月溪井 静水位 43 阶升 0.351 洱源井 静水位 50 阶升—阶降复原 0.016 南华井 静水位 153 阶降 0.009 腾冲洞山井 静水位 152 阶升—复原 0.005 腾冲台井 静水位 153 阶降—复原 0.005 牟定井 静水位 175 缓降 0.114 德宏井 静水位 195 阶升—阶降复原 0.015 隆阳井 动水位 93 阶升 0.028 剑川井 动水位 97 上升 0.012 丽江井 动水位 140 阶降—复原 0.005 楚雄井 动水位 180 阶升 0.066 普洱大寨井 动水位 347 阶升 0.032

表 1 漾濞MS 6.4地震引起云南数字化井水位同震响应特征 Table 1 Co-seismic response characteristics of digital well water level in Yunnan Province caused by Yangbi MS 6.4 earthquake

由图 1和表 1可知，漾濞M_S 6.4地震引起水位同震响应的观测井主要分布在震中周围200 km范围内的2个区域。一是沿红河断裂带北侧呈NW走向分布在丽江—剑川—洱源—大理—南华—楚雄一带，共8口井记录到水位同震响应；二是分布在小滇西地区，小滇西地区位于漾濞M_S 6.4地震震中西南方向，共有保山井、隆阳井、腾冲台井、腾冲洞山井、德宏井、施甸井、昌宁井等7口水位井正常观测，其中，隆阳井、腾冲台井、腾冲洞山井、德宏井等4口井记录到水位同震响应。另外，200 km范围外，只有滇西南地区的普洱大寨井记录到水位同震响应。

由表 1可知，对于漾濞M_S 6.4地震，云南井水位产生同震响应的观测井其同震响应形态基本都是阶变型变化，根据其动态变化大致可分为阶升（上升）型、阶降（或缓降）型、阶变（阶升、阶降）—复原型等3种变化类型（图 2）。

图 2(Fig.2)

图 2 2021年5月21日（牟定井为21—22日）水位同震响应动态变化 Fig.2 Co-seismic response dynamic curve of water level

阶升（上升）型同震响应变化表现为当地震发生时，井水位台阶上升或快速上升，之后一直持续在高水位变化，短时间内不会恢复到震前的水位动态背景位置[图 2（a）、2（b）]。漾濞M_S 6.4地震发生时，下关井、大理月溪井、隆阳井、楚雄井、普洱大寨井等5口井表现为阶升变化，升幅最小的是隆阳井，阶升0.028 m，升幅最大的是大理月溪井，阶升0.351 m；剑川井表现为快速上升变化，升幅0.012 m。下关井、大理月溪井、隆阳井震中距小于100 km，楚雄井、普洱大寨井大于100 km，下面以楚雄井为例具体分析。

楚雄井为动水位井，位于川滇菱形块体中南部的滇中块体内楚雄盆地西部边缘，距楚雄—建水断裂较近，距漾濞地震震中180 km，井孔深度为502.80 m，井内50 m处作永久性止水。观测段含水层位于地下57—376 m，岩性主要为含钙质泥岩夹泥灰岩及粉砂质泥岩，地层为白垩系上统江底河组下段（K₂J₁），含裂隙溶孔承压水。2021年5月21日21:48:34漾濞M_S 6.4地震时，楚雄井水位21:50开始从0.563 m台阶式上升至0.629 m，升幅0.066 m，之后一直持续在高水位动态变化，震后残留水位阶。

阶降（或缓降）型同震响应变化表现为当地震发生时，水位台阶下降或缓慢下降，之后一直持续在低水位变化，短时间内不会恢复到震前的水位动态背景位置[图 2（c）、2（d）]。漾濞M_S 6.4地震发生时，南华井表现为阶降变化，降幅0.009 m；牟定井表现为缓慢下降变化，降幅0.114 m，下面以南华井为例具体分析。

南华井为静水位井，位于南华—楚雄断裂带北段，距漾濞地震震中153 km，井深212.5 m，水位埋深约2.3m。观测段含水层埋深75.60 m，岩性为石英砂岩，地层为中生代白垩西系下统高丰寺组（K₁g₃），含基岩裂隙脉状承压水。水位年变动态受降雨影响明显，为春夏低—秋冬高型变化，日变动态特征为双峰—双谷型变化。漾濞M_S 6.4地震时，南华井水位在日变低谷值时产生同震阶降变化，21:50从2.883 m下降到2.892 m，降幅0.009 m，之后一直持续低于水位动态背景值，震后残留水位阶，直至6月中旬进入雨季后才开始转折上升。

阶变（阶升、阶降）—复原型同震响应变化表现为当地震发生时，井水位产生同震阶变（阶升或阶降），之后快速恢复到震前的水位动态背景位置，不会残留水位阶[图（e）、2（f）]。漾濞M_S 6.4地震发生时，洱源井、丽江井、德宏井、腾冲洞山井、腾冲台井等5口井表现为阶变（阶升、阶降）—复原型同震响应变化。震中距最小的是洱源井，相距50 km，其阶变幅度0.016 m。其他4口井震中距都达到或超过140 km，其中，震中距最大的是德宏井，其阶变幅度0.015 m，丽江井、腾冲洞山井、腾冲台井等3口井阶变幅度都是0.005 m。下面以德宏井为例具体分析。

德宏井为静水位井，位于龙陵—瑞丽压扭性大断裂下盘芒市断陷盆地东南缘、龙陵—畹町断裂北部，距漾濞地震震中195 km，井深200.17 m，水位埋深约2.15 m。观测段含水层埋深92.6 m，岩性呈中厚层状泥灰岩，地层年代为中生代侏罗系中统（J₂m），含承压岩溶水。水位年变动态受降雨影响明显，为冬春低—夏秋冬高型变化，日变动态特征为双峰—双谷型变化。漾濞M_S 6.4地震时，德宏井水位在日变低谷值时产生同震阶升—阶降—复原形态变化，21:52从2.409 m阶升至2.398 m，然后阶降至2.413 m，再阶升至2.410 m恢复正常动态，变幅为0.015 m。地震后，水位恢复动态背景值变化，未残留震后水位阶。

通过以上分析可知，漾濞M_S 6.4地震对云南井水位观测的影响主要集中在震中周围200 km范围内。从变化形态和分布上，井水位的同震响应都表现出下列特征：①阶升型同震响应幅度相对大于阶变（阶升、阶降）—复原型和阶降（缓降）型同震响应；②距震中100 km内，观测井无论是在红河断裂带一线还是在小滇西地区，井水位同震响应均主要表现为阶升变化，只有洱源井是阶升—阶降复原变化；③距震中140 km外，井水位同震响应主要表现为阶变（阶升、阶降）—复原变化和阶降（缓降）变化，只有楚雄井和震中距347 km的普洱大寨井是阶升变化。

廖丽霞等（2013）通过对不同采样率的水位同震响应能力及其特征分析认为，不同采样率的仪器记录水位同震响应存在较大差异，不但形态特征上明显不同，响应能力及响应幅度也差别较大。目前，云南有部分井孔水位观测采用SWY-Ⅱ型数字水位仪，具备秒采样功能，但数据收集不方便。使用云南地震台地球物理研究室下载保存的云南省水位秒钟值数据，收集到漾濞M_S 6.4地震时楚雄井、德宏井、洱源井、建水井、剑川井、昆明台井、泸西井、勐腊井、孟连井、弥渡井、姚安井、永胜井、元谋井、麒麟井、昭通井、师宗井等16口井的水位秒采样数据。

漾濞M_S 6.4地震时，楚雄井、德宏井、洱源井、剑川井等4口井水位分钟值、秒钟值都明显记录到同震响应变化。图 3为分钟值、秒钟值记录同震响应形态对比。如图 3（a）所示，楚雄井分钟值水位阶升，秒钟值同样水位阶升，只是在阶升前期记录到短暂的振荡变化，表明楚雄井水位同震响应阶升是可信的。洱源井、德宏井分钟值水位阶变—复原变化，而其秒钟值则表现为振荡—复原变化，表明其水位同震阶变—复原响应变化可能主要是观测仪器采样率不足、数据部分缺失造成的。基于此，漾濞地震水位同震阶变—复原型变化可大致认为是时间较短的振荡—复原变化。

图 3(Fig.3)

图 3 2021年5月21日分钟值、秒钟值记录同震响应形态对比 蓝色为分钟值；黑色为秒钟值 （a）分钟值与秒钟值都记录同震响应的观测井；（b）只有秒钟值记录到同震响应的观测井 Fig.3 Contrast chart of co-seismic response form by minute and second values recorded

其他12口井在漾濞M_S 6.4地震时未记录到水位分钟值同震响应变化或记录不明显，而明显记录到秒钟值同震响应变化[图 3（b）]，但同震响应持续时间只有几秒至几十秒。弥渡井、永胜井、姚安井、元谋井等4口井井震距小于200 km，只有秒钟值明显记录到同震响应。其他8口井井震距均大于200 km，最远的泸西井井震距415 km，同样只有秒钟值记录到同震响应。由图 3可知，秒钟值水位动态变化更为完整，同震响应幅度也相对较大；而分钟值水位动态变化则存在部分信息缺失，甚至未记录到同震响应变化。因此，一定范围内的观测井水位分钟值未记录到同震响应的主要原因是采样率不够，而秒钟值记录同震响应的时间较短可能与观测井—含水层系统、井孔所在位置的水文地质条件和地质构造条件有关。

Manga等（2007）分析认为，地下水的各种应力响应可分为断层错位产生的静应力和来自地震波的动应力。几个震源体范围内的承压含水层的水位阶变是源于多孔弹性介质对静态应力场的响应（Wakita，1975；Roeloffs，1998；Brodsky et al，2003）。Zhang等（2011）认为，在震中距小于1.5倍地震烈度长轴的范围内主要由静态应变引起水位同震阶变，超过该范围的水位同震阶变响应无法用静态应力场的变化来解释。漾濞M_S 6.4地震烈度的Ⅵ度区等震线长轴呈NNW走向，长度为106 km，估算1.5倍烈度长轴为159 km。

分析认为，地震所导致的地表破裂根据其发震断层的走向具有一定的方向性，地震烈度区域随着地表破裂的程度而展布，由此绘制的地震烈度图大多不是圆形区域构成，若只依据烈度长轴估算静态应力（或应变）的影响范围，则估算影响范围为圆形区域，可能与实际相差较大。根据漾濞M_S 6.4地震时周边井水位的同震响应变化形态和分布特征，距震中159 km内同震响应井较多，变化形态各异，也不能完全用静态应力（或应变）的影响来解释。在此基础上，增加地震烈度短轴的影响，漾濞地震烈度短轴长76 km，估算1.5倍烈度短轴为114 km，由长轴和短轴划定的椭圆区域可能更接近静态应力（或应变）的影响范围（图 4）。由此可见，可以依据漾濞M_S 6.4地震时各观测井所处位置的应力状态分析水位同震响应产生的力学机理。

图 4(Fig.4)

图 4 漾濞地震静态应力影响范围示意图 F1为小江断裂带；F9为澜沧江断裂带；F10为怒江断裂带；F11为大盈江断裂带；F12为龙陵—瑞丽断裂带；F15为无量山—营盘山断裂；F17为红河断裂带；F18为楚雄—建水断裂 Fig.4 Sketch map of influence area of static stress of Yangbi earthquake

由图 4可知，漾濞M_S 6.4地震静态应力影响范围内有下关井、大理月溪井、洱源井、剑川井、隆阳井等5口井产生同震响应，除洱源井外，其他4口井水位同震响应均为阶升变化。在漾濞地震前，这些井孔所在区域共同构成了孕育漾濞地震的区域构造应力（静态应力）场，区域压应力逐渐增强，至漾濞地震时，震源应力场的向外推挤作用叠加使这些井孔所在区域含水层受到瞬时增大的压应力，从而导致含水层介质内水瞬间流入井孔造成水位阶升，之后压应力会持续存在一段时间，直至产生新的应力调整，因此震后一段时间内这些井孔水位残留震后水位阶。

在静态应力影响范围以外，有8口井产生同震响应。其中，德宏井、腾冲洞山井、腾冲台井、丽江井阶变—复原型变化结合其秒钟值动态变化实际上是短暂的振荡—复原变化，可认为是地震波作用于含水层导致含水层介质弹性变形引起的。楚雄井水位阶升变化其秒钟值数据在水位阶升前期记录到了振荡变化，表明含水层介质在水位阶升前期经历了弹性变形。分析认为，楚雄井水位阶升可能是受漾濞地震地震波的触发作用，前期产生振荡变化，之后触发区域压应力进一步增强，产生塑性变形，从而表现为阶升变化，并且阶升后残留水位阶，一直持续高水位动态变化。杨竹转等（2005）通过对普洱大寨井水位同震阶变的研究发现，无论地震在哪一个方位发生，无论地震力学机制如何，普洱大寨井水位同震阶变总是上升。南华井位于受静态应力、动态应力影响的交界处，水位可能受静态应力、动态应力的综合作用表现为阶降变化；牟定井水位缓降变化可能是受地震波作用的影响，井—含水层介质产生塑性拉张变形，但拉张变形的程度不大，表现为缓慢下降变化。

综上所述，从井孔所处区域受力影响的状态可以初步分析各种同震响应变化形态的形成机理。但是，对于每一口井水位同震响应变化机理的深入研究，需要结合各观测井所处位置的区域应力状态、同震水位变化形态、观测井的井—含水层系统、井区地质构造条件、水文地质条件等进行综合分析，才能得出更科学合理的解释。

通过对漾濞M_S 6.4地震云南数字化井水位分钟值数据同震响应特征、观测仪器采样率对井水位映震能力的影响进行分析，得到了以下几点认识。

（1）漾濞M_S 6.4地震时有13口井记录到明显的水位同震响应变化，占云南总数字化水位观测井的22.8%。这些水位同震响应观测井主要分布在震中周围200 km范围内的2个区域。一是沿红河断裂带北侧呈NW走向分布在丽江—剑川—洱源—大理—南华—楚雄一带；二是分布在小滇西地区，位于漾濞M_S 6.4地震震中西南方向。在200 km范围以外，只有滇西南地区的普洱大寨井记录到水位同震响应变化。表明漾濞地震引起云南数字化井水位产生同震响应的能力总体不强，影响范围不大。

（2）对于漾濞M_S 6.4地震，云南数字化井水位同震响应形态基本都是阶变型变化，根据其动态变化大致又可细分为阶升（上升）型、阶降（或缓降）型、阶变（阶升、阶降）—复原型等3种变化类型。

（3）对于漾濞M_S 6.4地震，云南数字化井水位的同震响应从变化形态和分布上都表现出下列特征：①阶升型同震响应幅度相对大于阶变（阶升、阶降）—复原型和阶降（缓降）型同震响应；②在静态应力影响范围内，观测井无论是在红河断裂带一线还是在小滇西地区，井水位同震响应主要表现为阶升变化；③在静态应力影响范围以外，井水位同震响应形态除了以阶变（阶升、阶降）—复原变化为主外，还有阶降（缓降）变化和阶升变化。

（4）对于井水位同震响应机理的探讨，在Zhang等（2011）研究的基础上，引入地震烈度短轴的影响，由地震烈度长轴和短轴共同划定的椭圆区域来划分漾濞M_S 6.4地震静态应力、动态应力的影响范围，依据漾濞地震时各观测井所处位置的应力状态，分析其水位同震响应产生的力学机理来解释漾濞地震时各观测井水位的同震响应变化形态可能更接近实际情况。因此，这种划分应该可以作为水位同震响应机理分析的一种参考分析依据。

（5）观测仪器采样率对井水位映震能力的影响较大，漾濞M_S 6.4地震时云南部分观测井没有记录到分钟值数据同震响应，主要原因是仪器采样率不够，而秒钟值记录同震响应的时间较短可能与观测井—含水层系统、井孔所在位置的水文地质条件、地质构造条件等有关。

水位秒钟值数据的同震响应变化形态更加完整，可为井水位同震响应机理研究提供更客观的资料。但是，目前水位秒钟值数据的收集和使用都存在不便，SWY-Ⅱ型数字水位仪具备秒采样功能，但没有自动下载、数据入库功能。如果需要秒钟值数据，需要手动下载保存。若不及时收取，容易造成数据缺失。秒钟值文件数据格式也不是常用软件可以直接使用的格式，需要转换才能使用。建议仪器生产厂家优化水位秒钟值数据下载保存和格式转化使用功能，以便发挥水位秒钟值数据更大的应用研究价值。