基于观测环境的洞体应变压性异常综合研究.pdf

第4 3卷 第9期2 0 2 3年9月大 地 测 量 与 地 球 动 力 学J o u r n a l o fG e o d e s ya n dG e o d y n a m i c sV o l.4 3N o.9S e p t.,2 0 2 3收稿日期:2 0 2 2-1 1-0 4项目来源:中国地震局震情跟踪定向工作任务(2 0 2 2 0 1 0 2 1 9);安徽省地震局科研创新团队。第一作者简介:刘莉,工程师,主要从事地形变监测、预测研究,E-m a i l:l i u l i 3 1 6 5 4 4 11 2 6.c o m。D O I:1 0.1 4 0 7 5/j.j g g.2 0 2 3.0 9.0 1 8文章编号:1 6 7 1-5 9 4 2(2 0 2 3)0 9-0 9 8 0-0 6基于观测环境的洞体应变压性异常综合研究刘 莉1 郑海刚1 李军辉1 倪红玉1 周冬瑞11 安徽省地震局,合肥市长江西路5 5 8号,2 3 0 0 3 1摘 要:从观测环境角度出发,通过对金寨台气象因素进行回归分析发现,气温、气压与洞体应变之间虽存在线性关系,但拟合结果与洞体应变压性下降实际观测应变变化量相差较大;洞体应变压性下降起始时间与梅山水库开始泄洪、下游史河水位陡升时间存在较强的相关性。三维不规则载荷模型数值模拟结果显示,洞体应变压性异常的前期快速下降变化与同期梅山水库载荷变化的相关性较小,与史河水位变化有较大相关性;洞体应变压性异常发生前后研究区构造应力场未发生显著变动。由此推定,金寨台洞体应变压性下降异常可能是周边场地环境载荷变化的影响所致。关键词:洞体应变;压性异常;气象;载荷;构造应力场中图分类号:P 3 1 5 文献标识码:A 定点地形变观测具有精度高、频带宽和连续性好等优点,是研究构造运动、断层滑移和地震前兆等诸多动力学问题的重要手段之一1。目前,一般采用S S-Y型伸缩仪进行定点台站洞体应变固体潮观测,获得的线应变固体潮观测资料可进行潮汐主应变、主方向和剪应变计算2-3,以分析观测场地的平面应变状态及区域应变场变化和地震前兆异常。就地震前兆研究而言,全面、深入地分析震前形变异常特征及其成因机制,对地震分析预测不仅重要也极具必要性。金寨台位于华北地块、下扬子地块与秦岭-大别断褶带的交汇处,附近有多条断裂带经过,区域构造 较 为 复 杂。台 站S S-Y伸 缩 仪EW分 量2 0 2 0-0 5-0 7开始趋势陡转,出现大幅持续下降形态的压性异常变化,至0 8-0 6(共计9 2d)下降幅度达1 27 8 41 0-1 0。虽然对台站观测系统和周边环境进行了多次调查,但至今仍未厘清该异常的性质及成因机制。本文根据金寨台区域水文、构造特征及地震活动性等,基于观测环境,拟从气象因素(气温、气压、降雨、洞温)变化、附近水体(水库、河流)载荷情况及区域地壳应力场变化等可能的成因机制角度,探讨此次异常的最佳成因,以期为该区域地震前兆研究提供助益。1 台站概况及区域构造背景金寨台形变观测山洞全长1 6 4.6 4 2m,引洞(近N S向)长5 8.8 4 2m,主洞正N S向长5 7.9 0 0m、正EW向长4 7.9 0 0m,主要有地倾斜和洞体应变2种观测项目,其中洞体应变观测仪器S S-Y型伸缩仪放置在主洞2个方向上(图1)。图1 金寨台山洞形变仪器布置分布F i g.1 L a y o u t a n dd i s t r i b u t i o no f c a v ed e f o r m a t i o ni n s t r u m e n t sa t J i n z h a i s t a t i o n 第4 3卷第9期刘 莉等:基于观测环境的洞体应变压性异常综合研究金寨台洞体应变测量原理是测量地表面两点间基线长度的相对变化量,即=L-LL=LL(1)式中,L为原地壳表面两点间的距离,称为基线长;L 为变化后地壳表面两点间的距离;L为基线的绝对变化量;为应变量,即单位长度的相对变化量。根据约定,压缩为负、拉张为正,即0时L拉张。2 EW分量压性异常特征判断金寨洞体应变自2 0 1 6-1 1开始观测以来EW分量除了出现过2次大幅下降的压性异常变化外,其他时段年变规律、趋势变化相对稳定(图2)。第1次压性异常变化为2 0 1 7-0 4-1 00 5-2 8(共计4 9d),EW分量变化幅度1 37 3 11 0-1 0,经核实认定为自然环境和场地环境干扰综合事件,因降雨、水库水位涨落、水库泄洪、史河水位上涨等因素综合影响引起的快速变化;第2次压性异常变化自2 0 2 0-0 5-0 7开始,EW分量趋势陡转,出现大幅持续下降,至0 8-0 6(共计9 2d)EW分量数据下降幅度达1 27 8 41 0-1 0;0 8-0 7开始EW分量在保持缓慢下降的趋势背景下出现3次张性鼓包现象,2 0 2 1-0 1底异常结束,数据曲线转入正常年变形态。2次变化均为快速下降压性变化,2 0 1 7年的下降变化结束后立即恢复正常年变,2 0 2 0-0 8-0 6下降变化结束后并未回归正常年变,而是出现了“三上三下”的应变调整形态(图2),并于2 0 2 1-0 1底回归正常年变。图2 金寨台洞体应变EW分量观测整点值曲线F i g.2 P o i n tv a l u ec u r v eo fEWc o m p o n e n to b s e r v a t i o no f t u n n e l b o d ys t r a i na t J i n z h a i s t a t i o n2次压性异常变化期间,观测系统运行正常,无检修仪器等人为干扰情况。为进一步分析金寨洞体应变EW分量大幅压性异常变化期间数据的可靠性,计算2 0 1 72 0 2 2年洞体应变N S分量和EW分量M 2波潮汐因子,结果见图3。结果表明,N S分量和EW分量的M 2波潮汐因子相对稳定。N S分量在2 0 1 8-0 91 0因多次检修更换传感器,导致潮汐因子变化较大,其他时间潮汐因子均在0.6 1 6左右;EW分量潮汐因子比N S分量潮汐因子大,但均在1.0 4 2左右。除受到仪器检修影响外,2个分量潮汐因子相对稳定,故认为仪器观测数据质量可靠,压性异常排除观测系统存在问题的可能性。对比金寨台其他形变观测项目及周边泾县台、淮北台、泗县台、麻城台洞体应变同时段观测数据,均未发现类似压性异常变化。图3 金寨台洞体应变N S、EW分量M 2波潮汐因子F i g.3 T i d a l f a c t o ro fM 2w a v eo fN S,EWc o m p o n e n t so f t u n n e l v o l u m es t r a i na t J i n z h a i s t a t i o n综上可知,在排除观测系统故障及人为干扰等因素后,金寨台2 0 2 0年EW分量的压性下降变化与2 0 1 7年的压性下降变化在形态上存在相似性和差异性,亟需从其他角度分析判断变化原因。本文拟从气象因素影响、周边场地环境变化、区域地壳应力场变动等方面分析此次压性下降变化的原因。3 EW分量压性异常3种可能的成因机制3.1 气象因素触发机制分析金寨台除了地倾斜和洞体应变观测外,还配有气象三要素辅助观测,收集历年来的气象观测资料及洞温观测资料。本文收集整理洞体应变观测以来气温、气压、洞温、降雨量、实测应变等变化量,建立回归分析模型,得到2 0 2 0年压性下降变化期间因气温、气压、洞温及降雨量变化引起的理论应变变化量,并与金寨台同时段实测洞体应变变化量进行对比分析。首先,针对各气象因素与应变变化之间进行一元线性回归分析,由于本文收集到的气象因素变化导致的洞体应变变化量与2 0 2 0年洞体应变189大 地 测 量 与 地 球 动 力 学2 0 2 3年9月EW向压性下降变化量相差较大,作图时存在压制现象,因此对所有气象因素变化量、洞体应变变化量作对数处理,再进行一元线性回归分析。如图4所示,气温及气压均与洞体应变存在明显的线性关系,但2 0 2 0-0 5-0 70 8-0 6实际观测应变变化量(图4中红色圆圈)与一元线性回归拟合曲线相差较远;降雨量及洞温均与洞体应变无明显线性关系。图4 应变变化回归分析F i g.4 R e g r e s s i o na n a l y s i so f s t r a i nc h a n g e其次,运用气温、气压、降雨量及洞温与洞体应变建立多元线性回归模型,表达式为:y=9.5 2 2+1.5 1 2x1+1 5.7 0 5x2+0.1 8 2x3+5 0.9 2 6x4(2)式中,y为洞体应变变化量,x1为气温变化量,x2为气压变化量,x3为降雨量,x4为洞温变化量。如图5所示,多元线性回归拟合曲线与洞体应变变化曲线相差较大,拟合优度为0.5 4 3,表明洞体应变变化的5 4.3%可由模型确定;误差方差为8 9.6 1 5,模型误差较大。将2 0 2 0-0 5-0 70 8-0 6气温、气压、降雨量、洞温变化值代入式(2),得到模型理论应变变化量为3 9 41 0-1 0,与同时段EW分量应变变化幅度1 27 8 41 0-1 0相差较大。图5 气温、气压、降雨、洞温与洞体应变多元线性回归分析F i g.5 M u l t i p l e l i n e a r r e g r e s s i o na n a l y s i so f a i rt e m p e r a t u r e,a i rp r e s s u r e,r a i n f a l l,c a v et e m p e r a t u r ea n dc a v eb o d ys t r a i n综上可知,平常情况下,金寨台洞体应变EW分量在气温、气压与应变之间存在线性影响。针对2 0 2 0-0 5-0 70 8-0 6洞体应变EW分量压性下降变化,无论是气温、气压的一元线性回归拟合结果,还是气温、气压、降雨量及洞温的多元线性回归拟合结果,均与实际观测应变变化量相差较大,说明此次大幅压性下降变化不是气象因素所致。3.2 水体载荷变化分析梅山水库位于金寨台西南方向约2k m,总库容2 3.3 7 0亿m3,常年蓄水约1 0亿m3。史河从梅山水库发源,从台站西北方向流过,最近距离约5 0 0m。收集梅山水库和史河的相关水文数据可以发现,2 0 2 0-0 5-0 7金寨台洞体应变EW分量发生压性下降时,梅山水库开始泄洪、下游史河水位陡升,时间上存在较强的相关性(图6)。图6 梅山水库库容变化、史河水位变化与洞体应变EW分量对比F i g.6 C o m p a r i s o no fEWc o m p o n e n t so f r e s e r v o i rc a p a c i t yc h a n g e,S h i h er i v e rw a t e r l e v e l c h a n g ea n dt u n n e lb o d ys t r a i n i nM e i s h a nr e s e r v o i r建立三维不规则形状载荷模型4推算压性下降变化时段内水库库容变化及史河水位变化造成的应变变化量。针对三维不规则形状载荷模型,可将该模型的总作用力F进行网格化,分别求取各网格内作用力f对台站观测点造成的垂向和水平N S、EW向位移,然后求所有网格块对观测点作用的矢量和,即可得到观测点垂向和水平N S、EW向位移。各网格内作用力f对台站观测点造成的水平EW、N S向位移为5:u=f x z4 G R3-f x4(G+)R(R+z)(3)v=f y z4 G R3-f y4(G+)R(R+z)(4)式中,u为x轴(即水平EW向)位移;v为y轴(即水平N S向)位移;f为各网格内载荷作用力;R=x2+y2+z2为各网格到台站观测点的距离;和G为拉梅常数,对于弹性模量E=4 3.3 2 8G P a,泊松比=0.1 3 2,G=E2(1+)1.9 1 4289 第4 3卷第9期刘 莉等:基于观测环境的洞体应变压性异常综合研究1 01 0P a,=E(1+)(1-2)6.8 6 51 09P a。首先建立梅山水库的三维不规则形状载荷模型,将其简略划分为5 8个截面为10 0 0m10 0 0m的网格块,如图7(a)所示,运用式(3)和式(4)分别计算不同水库库容时这5 8个网格块对洞体应变观测点造成的位移矢量和,结果如图8所示。根据模型计算结果,2 0 2 0-0 5-0 70 8-0 6金寨台洞体应变EW分量受水库载荷变化引起的理论位移变化量为-2.0 3 0 mm。金 寨台洞体应 变仪EW分量基线长L=2 8.2 8m,2 0 2 0-0 5-0 70 8-0 6观测值总变化量为-1 27 8 41 0-1 0,运用式(1)可换算得到台站洞体应变观测基线的绝对变化量L为-0.0 3 6mm,与数值模拟的同时段理论位移变 化 量-2.0 3 0 mm相 差2个 数 量 级,说 明2 0 2 0年金寨台洞体应变压性异常变化与同时期梅山水库载荷变化的相关性较小。其次,建立史河的三维不规则形状载荷模型,将其简略划分为1 4个截面为5 0 0m5 0 0m的网格块,如图7(b)所示,运用式(2)和式(3)分别计算史河不同水位时1 4个网格块对洞体应变观测点造成的位移矢量和,结果如图9所示。根据模型计算结果,2 0 2 0-0 5-0 70 8-0 6金寨台洞体应变EW分量受史河水位变化(水位最大变幅为2.4 6m)引起的理论位移变化量为-0.0 2 3mm,与洞体应变实测位移变化量-0.0 3 6mm相差不大,考虑与三维不规则载荷模型的地形近似、网格近似划分、未考虑沉积层影响等近似条件的设置有关,因此认为2 0 2 0年压性下降的前期快速下降变化与史河水位变化有较大相关性。图7 梅山水库、史河不规则形状载荷网格化示意图F i g.7 M e s h i n gd i a g r a mo f i r r e g u l a r s h a p e l o a do fM e i s h a nr e s e r v o i ra n dS h i h er i v e r图8 梅山水库三维不规则载荷模型水平方向位移解F i g.8 H o r i z o n t a l d i s p l a c e m e n t s o l u t i o no f 3 Di r r e g u l a r l o a dm o d e l o fM e i s h a nr e s e r v o i r图9 史河三维不规则载荷模型水平方向位移解F i g.9 H o r i z o n t a l d i s p l a c e m e n t s o l u t i o no f 3 Di r r e g u l a r l o a dm o d e l o fS h i h er i v e r389大 地 测 量 与 地 球 动 力 学2 0 2 3年9月 结合图8和图9还可以看出,无论是水库库容变化还是史河水位变化,模型中N S分量位移量均大于EW分量,说明梅山水库库容变化或史河水位变化对金寨台周边区域位移的影响N S向大于EW向。3.3 区域地壳应力场变动分析为进一步分析2 0 2 0年压性异常变化时段金寨台所在区域(3 0.2 3 3.2 N,1 1 3.8 1 1 7.8 E)地壳应力场的变化情况,采用S n o k e方法对正常时段(2 0 1 8-0 1-0 12 0 2 0-0 5-0 6)和洞体应变压性下降异常时段(2 0 2 0-0 5-0 72 0 2 1-0 1-3 1)研究区ML2.3地震的震源机制解进行解算,共得到3 0个较可靠的震源机制解,结果如图1 0所示。之后,利用力轴张量法进行反演,获取该区域地壳应力场的时空变化特征(图1 1)。图1 0 压性异常前后研究区域内3 0次ML2.3地震的震源机制解F i g.1 0 F o c a lm e c h a n i s ms o l u t i o n so f 3 0ML2.3e a r t h q u a k e s i nt h es t u d i e da r e ab e f o r ea n da f t e rt h ec o m p r e s s i v ea n o m a l y结果显示,正常时段该区域主压应力1的方位角为8 1.4、倾角为1 9.2,呈N E E向、低倾角;主张应力3的方位角为1 7 3、倾角为4.4,呈S S E向、近水平,表明正常时段该区域处于近EW向低倾角压应力和近N S向近水平张应力的区域构造应力作用。而2 0 2 0年压性异常时段,该区域主压应力1的方位角为1 0 4.5、倾角为6.6,呈S E E向、应力轴接近水平;主张应力3的方位角为1 9 7.3、倾角为2 3,呈S SW向、低倾角,可见虽然也处于近空心圆、空心三角形分别代表P轴、T轴,实心圆、实心菱形和实心正方形分别代表主压应力1、中间应力2、主张应力3图1 1 震源机制解P、T轴和主应力轴1、2、3在下半球吴尔夫网上的投影F i g.1 1 P-a x i s,T-a x i s,a n dp r i n c i p a l s t r e s sa x i so ff o c a lm e c h a n i s ms o l u t i o n1,2,3p r o j e c t i o no nt h eW o l f fn e t w o r ki nt h e l o w e rh e m i s p h e r eEW向压应力和近N S向张应力的区域构造应力作用中,但相较于正常时段而言,主压应力方位角发生顺时针偏转,偏转幅度约2 4。倪红玉等6利用阻尼区域应力反演方法研究表明,安徽地区受近EW向水平挤压和近N S向水平拉张的现代构造应力场作用,最大主压应力方向在空间上从北到南也存在一定程度的偏转。金寨台洞体应变EW分量压性异常发生前后,研究区主压应力轴1和主张应力轴3的方位角总体趋势较为一致,存在一定程度的顺时针偏转,这与倪红玉等6的研究结果一致。由此说明,金寨台洞体应变EW分量压性异常发生前后该区域构造应力场未发生显著变动,2 0 2 0年金寨台洞体应变压性异常变化由区域构造应力场变化引起的可能性较小。4 结 语针对金寨台洞体应变EW分量在2 0 2 0年出现的压性下降异常,在排除观测系统故障、人为干扰等因素后,提出并分析了3种可能的成因机制,得到以下初步认识:1)通过气温、气压、降雨、洞温的一元或多元线性回归分析发现,气温变化、气压变化与应变变化之间虽然存在线性影响,但拟合结果与实际观测应变变化量相差较大,说明此次大幅压性下降变化不是气象因素所致。2)2 0 2 0-0 5-0 7金寨台洞体应变EW分量发生压性下降时,梅山水库开始泄洪、下游史河水位陡升,时间上存在较强的相关性;三维不规则载荷模型数值模拟结果显示,洞体应变压性的前期快速下降变化与同时期梅山水库载荷变化的相关性较489 第4 3卷第9期刘 莉等:基于观测环境的洞体应变压性异常综合研究小,与史河水位变化有较大相关性,说明此次大幅压性下降变化由梅山水库泄洪触发,之后受下游史河蓄水量较大影响所致。3)洞体应变EW分量压性异常发生前后,研究区主压应力轴1和主张应力轴3的方位角总体趋势较为一致,进一步说明异常前后该区域构造应力场未发生显著变动,因此洞体应变压性异常变化由区域构造应力场变化引起的可能性较小。初步推定,金寨台洞体应变EW分量压性下降异常很可能是受周边场地环境载荷变化的影响所致。受水文资料限制,本文对地形变化、地下介质分布、物性参数、相关模型参数等进行了简化和近似,数值模拟得到的结果与实际位移变化量在数量级上基本一致,可为今后定点形变出现大幅变化后的判定提供定性与定量分析的思路。下一步考虑结合同台其他定点形变观测资料作进一步对比分析,为该区域震情跟踪提供助益。参考文献1 杨小林,王军,王希彬,等.陕西韩城台巨幅地倾斜异常的成因分析J.地震学报,2 0 1 8,4 0(6):7 6 0-7 7 3(Y a n gX i a-o l i n,W a n gJ u n,W a n gX i b i n,e ta l.Wh a tC a u s e st h eR e-m a r k a b l eT i l tA n o m a l i e sa t t h eH a n c h e n gG e o d y n a m i cO b-s e r v a t o r y i nS h a a n x i P r o v i n c e?J.A c t aS e i s m o l o g i c aS i n-i c a,2 0 1 8,4 0(6):7 6 0-7 7 3)2 刘序俨,李平,张雁滨.地表的面应变和体应变固体潮理论值计 算 及 其 调 和 分 析 J.地 壳 形 变 与 地 震,1 9 8 8,8(4):3 5 4-3 5 8(L i uX u y a n,L i P i n g,Z h a n gY a n b i n.C a l c u-l a t i o no fT h e o r e t i c a lV a l u eo fC u b i cE x p a n s i o nT i d ea n dI t sH a r m o n i cA n a l y s i sJ.C r u s t a lD e f o r m a t i o na n d E a r t h-q u a k e,1 9 8 8,8(4):3 5 4-3 5 8)3 蒋骏,张雁滨.潮汐线应变组合观测的物理力学实质及其调和分析J.地球物理学报,1 9 9 4,3 7(增2):2 0 4-2 1 2(J i a n gJ u n,Z h a n g Y a n b i n.T h eP h y s i c s-M e c h a n i c a lS i g n i f i c a n c eo fT i d a l S t r a i nC o m b i n a t i o na n dI t sH a r m o n i cA n a l y s i sJ.C h i-n e s e J o u r n a l o fG e o p h y s i c s,1 9 9 4,3 7(S 2):2 0 4-2 1 2)4 闫伟,牛安福,陈国琴.不规则载荷影响地表近场垂向位移和倾斜观测的定量计算J.中国地震,2 0 1 5,3 1(2):2 9 7-3 0 7(Y a nW e i,N i uA n f u,C h e nG u o q i n.C a l c u l a t i o no fG r o u n dV e r t i c a lD i s p l a c e m e n t s,T i l t sa n dT h e i rD i s t r i b u-t i o nC h a r a c t e r i s t i c sa r o u n da nI r r e g u l a r l yL o a dJ.E a r t h-q u a k eR e s e a r c h i nC h i n a,2 0 1 5,3 1(2):2 9 7-3 0 7)5 邱泽华.钻孔应变观测点到地面载荷干扰源最小“安静”距离的理论分析J.岩石力学与工程学报,2 0 0 4,2 3(2 3):40 6 3-40 6 7(Q i uZ e h u a.T h e o r e t i cS t u d i e so nt h eS h o r t e s tQ u i e tD i s t a n c eb e t w e e nB o r e h o l eS t r a i nO b s e r v a t o r yS i t ea n dS u r f i c i a lM e c h a n i c a lD i s t u r b a n c eS o u r c e sJ.C h i n e s eJ o u r n a l o f R o c k M e c h a n i c s a n d E n g i n e e r i n g,2 0 0 4,2 3(2 3):40 6 3-40 6 7)6 倪红玉,洪德全,刘泽民,等.安徽地区震源机制解与应力场特征研究J.大地测量与地球动力学,2 0 1 7,3 7(增2):7 6-8 0(N iH o n g y u,H o n gD e q u a n,L i uZ e m i n,e ta l.R e-s e a r c ho nF o c a lM e c h a n i s mS o l u t i o n s a n dC h a r a c t e r i s t i c so fS t r e s sF i e l do fA n h u iR e g i o nJ.J o u r n a lo fG e o d e s ya n dG e o d y n a m i c s,2 0 1 7,3 7(S 2):7 6-8 0)C o m p r e h e n s i v eS t u d yo nC o m p r e s s i v eA n o m a l yo fC a v e-S t r a i nB a s e do nO b s e r v a t i o nE n v i r o n m e n tL I UL i1 ZHENGH a i g a n g1 L IJ u n h u i1 N IH o n g y u1 ZHO UD o n g r u i11 A n h u iE a r t h q u a k eA g e n c y,5 5 8W e s t-C h a n g j i a n gR o a d,H e f e i 2 3 0 0 3 1,C h i n aA b s t r a c t:F r o mt h ep e r s p e c t i v eo fo b s e r v a t i o ne n v i r o n m e n t,t h r o u g ht h er e g r e s s i o na n a l y s i so fm e t e-o r o l o g i c a l f a c t o r s,w e f i n d t h a t a l t h o u g ht h e r e i s a l i n e a r e f f e c tb e t w e e nt e m p e r a t u r e c h a n g e,p r e s s u r ec h a n g ea n ds t r a i nc h a n g e,t h e f i t t i n gr e s u l t s a r eq u i t ed i f f e r e n t f r o mt h ea c t u a l o b s e r v e ds t r a i nc h a n g ew h e nt h e s t r a i np r e s s u r ed r o p s.T h e r e i s a s t r o n gc o r r e l a t i o nb e t w e e n t h e s t a r t t i m eo f t h ed e c l i n e a n dt h e t i m ew h e nt h e M e i s h a nr e s e r v o i rs t a r t st od i s c h a r g ea n dt h et i m ew h e nt h ew a t e rl e v e lo ft h ed o w n s t r e a mS h i h e r i v e r r i s e s s h a r p l y.T h en u m e r i c a l s i m u l a t i o n r e s u l t s o f t h e t h r e e-d i m e n s i o n a l i r r e g-u l a r l o a dm o d e l s h o wt h a t t h ee a r l yr a p i dd e c l i n eo f t h ea n o m a l o u ss t r a i nc o m p r e s s i b i l i t yo f t h ec a v eb o d yh a s l i t t l e c o r r e l a t i o nw i t ht h e l o a dc h a n g eo f t h eM e i s h a nr e s e r v o i r i nt h e s a m ep e r i o d,a n dh a s al a r g ec o r r e l a t i o nw i t ht h ew a t e r l e v e l c h a n g eo f t h eS h i h er i v e r.T h e t e c t o n i cs t r e s s f i e l d i nt h es t u d ya r e ad i dn o t c h a n g e s i g n i f i c a n t l yb e f o r e a n da f t e r t h eo c c u r r e n c eo f c a v e-s t r a i na n dc o m p r e s s i v e a n o m a-l y.T h e r e f o r e,w ep r e l i m i n a r i l y i n f e r t h a t t h ea b n o r m a ld r o po fs t r a i na n dp r e s s u r e i nJ i n z h a i s t a t i o nc a v em a yb ec a u s e db yt h ee n v i r o n m e n t a l l o a dc h a n g e s i nt h es u r r o u n d i n gs i t e s.K e yw o r d s:c a v e-s t r a i n;c o m p r e s s i v ea n o m a l y;m e t e o r o l o g y;l o a d;t e c t o n i cs t r e s s f i e l dF o u n d a t i o ns u p p o r t:T h eE a r t h q u a k eT r a c k i n gT a s ko fC E A,N o.2 0 2 2 0 1 0 2 1 9;S c i e n t i f i cR e s e a r c ha n d I n n o v a t i o nT e a mo fA n h u iE a r t h-q u a k eA g e n c y.A b o u t t h e f i r s t a u t h o r:L I UL i,e n g i n e e r,m a j o r s i nt e r r a i nc h a n g em o n i t o r i n ga n df o r e c a s t i n g,E-m a i l:l i u l i 3 1 6 5 4 4 11 2 6.c o m.589