Analysis of the Jingle Well Water-level Anomaly Based on Numerical Simulation and Aquifer Stress Inversion
-
摘要: 利用数值模拟与含水层垂向应力反演的方法,对降雨与河流荷载作用下静乐井区应力进行定量计算,并从区域应力场与同一构造区其他前兆测项变化特征对2种方法计算结果进行分析。数值模拟结果显示,降雨与河流荷载的共同作用使静乐井所在区域产生的最大垂向位移约为0.7mm,最大垂向应力约为30hPa;同时段静乐井-含水层垂向应力反演结果显示静乐井-含水层最大垂向应力约为79hPa;数值模拟与静乐井-含水层垂向应力反演结果均表明,静乐井水位高值异常与同时段降雨量增多、河流荷载效应增强关系密切,但2种结果存在一定差异。通过可反映井区应力场变化特征的山西地震带3级地震缺震和b值、穿过静乐井的GPS基线、地震矩释放、同一构造单元同时段其他前兆测项与静乐井水位高值异常对比等分析结果可知,2017年静乐井水位高值异常期间井区应力场存在显著异常,推测静乐井水位高值异常除受降雨与河流荷载作用的影响外,也可能受构造活动增强的影响。Abstract: Using the methods of numerical simulation and aquifer stress inversion, we calculate the stress values in Jingle well area under rainfall and river load. The calculation results of the two methods are analyzed from the regional stress field and other precursor observation values in the same tectonic area. The numerical simulation results show that the maximum vertical displacement and the maximum vertical stress caused by the combined load of river and rainfall are about 0.7mm and 30 hPa. At the same time, the inversion results of the vertical stress of the interval well-aquifer show that the maximum vertical stress of Jingle well-aquifer is about 79 hPa. Both numerical simulation and vertical stress inversion results of well-aquifer show that the abnormal high water level of Jingle well is closely related to the increase of rainfall and the enhancement of river load effect in the same period, but there are certain differences between the two results. Through the comparative analysis of the magnitude 3 earthquake deficiency and earthquake b value in Shanxi seismic belt, GPS baseline passing through jingle well, earthquake moment release, other precursor observation values of the same structural unit and the anomaly of jingle well water level high value, which can reflect the variation characteristics of the well stress field, it is found that there are significant anomalies in the well stress field during the anomaly of jingle well water level high value in 2017. It is speculated that the anomaly of jingle well water level high value not only is influenced by river and rainfall load, but also contains the factors of structural activity enhancement.
-
表 1 静乐井-含水层多种参数计算结果
Table 1. Calculation results of various parameters of Jingle well-aquifer
时间 气压系数/(mm·hPa−1) 潮汐效率/(mm·10−9) 体积压缩系数/pa−1 孔隙度/% 固体体积压缩系数/pa−1 水位变化量/m 含水层垂向应力/hPa 2016-01 4.618 0.191 69.285 1.806 5.117 0.013 1.4019 2016-02 2.522 0.330 9.766 2.442 0.751 0.007 0.5703 2016-03 3.267 0.308 167.63 1.008 6.112 −0.004 −0.8453 2016-04 5.997 0.180 960.22 0.434 18.121 −0.014 −5.9814 2016-05 4.328 0.342 530.25 0.516 10.970 0.005 1.8874 2016-06 5.032 0.264 4926.8 0.053 10.884 −0.0009 −3.2224 2016-07 5.742 0.316 221.25 0.657 6.260 −0.003 −0.8655 2016-08 9.897 0.238 592.66 0.203 5.651 −0.017 −15.5497 2016-09 3.960 0.167 170.08 1.215 8.059 −0.006 −0.9238 2016-10 4.642 0.148 90.599 1.785 6.621 −0.007 −0.8333 2016-11 2.459 0.269 162.6 1.366 6.875 0.002 0.2584 2016-12 3.767 0.147 1709.5 0.106 6.9188 0.008 14.8921 2017-01 4.007 0.154 123.43 1.52 7.348 0.007 0.8916 2017-02 2.861 0.214 254.43 1.089 9.405 0.001 0.1522 2017-03 2.784 0.227 1454.1 0.071 3.449 0.001 1.3134 2017-04 3.808 0.292 148.14 0.997 5.684 −0.007 −1.3365 2017-05 6.093 0.206 600.47 0.494 12.952 0.011 4.2036 2017-06 4.756 0.285 84.983 1.233 4.318 0.005 0.8134 2017-07 6.225 0.249 114.4 1.032 5.173 −0.004 −0.7851 2017-08 18.400 0.191 4233.5 0.018 36.619 −0.03 −323.3068 2017-09 15.148 0.164 751.41 0.388 14.207 −0.015 −7.5280 2017-10 6.703 0.269 3324.2 0.055 8.1277 −0.012 −42.1179 2017-11 3.473 0.161 813.31 0.676 20.424 0.010 2.9897 2017-12 — 0.246 — — — — — -
车用太, 鱼金子, 刘五洲, 1999.华北北部地区3次强震前地下流体异常场及其形成与演化机理.中国地震, 15(2):139-150. http://www.cqvip.com/Main/Detail.aspx?id=3579726 车用太, 鱼金子, 2004.地下流体典型异常的调查与研究.北京:气象出版社. 丁风和, 哈媛媛, 王勇等, 2015.基于数字化水位的张渤带地区构造应力场时序特征分析.地震, 35(2):133-138. doi: 10.3969/j.issn.1000-3274.2015.02.014 黄辅琼, 晏锐, 陈颙等, 2004.利用深井地下水位动态研究大华北地区现今构造应力场状态.地震, 24(1):112-118. http://d.old.wanfangdata.com.cn/Periodical/diz200401017 李春洪, 陈益惠, 田竹君, 1990.井-含水层系统对固体潮的动态响应及其影响因素.中国地震, 6(2):37-45. http://www.cqvip.com/Main/Detail.aspx?id=410642 刘雁冰, 裴顺平, 2017.汶川地震前后b值的时空变化及构造意义.地球物理学报, 60(6):2104-2112. http://www.cqvip.com/QK/94718X/20176/672463573.html 刘峡, 马瑾, 占伟等, 2013.汶川地震前后山西断陷带的地壳运动.大地测量与地球动力学, 33(3):5-10. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dkxbydz201303002 刘瑞春, 李自红, 赵文星等, 2014.汶川M8.0地震前后山西地震带水平形变场变化特征研究.地震工程学报, 36(3):634-638. http://d.old.wanfangdata.com.cn/Periodical/xbdzxb201403034 孙小龙, 刘耀炜, 晏锐, 2011.利用水位资料反演华北地区构造应力场变化.地震, 31(2):42-49. http://d.old.wanfangdata.com.cn/Periodical/diz201102005 王爱英, 张淑亮, 杨占山, 1998.山西静乐井水位高值异常及其与华北地震活动的关系.地震, 18(2):201-205. doi: 10.3321/j.issn:0253-3782.1998.02.012 吴小平, 1990.b值物理机制的再探讨.西北地震学报, 12(3):1-13. http://d.old.wanfangdata.com.cn/Thesis/Y065065 鱼金子, 车用太, 张大维等, 1994.降雨对汤坑群井水位动态影响的剖析.地震研究, 17(3):273-280. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK199400058071 张淑亮, 刘巍, 王爱英, 1997.从地下流体附加力源看晋5-1井的水位变化.山西地震, (4):27-31. https://www.zhangqiaokeyan.com/academic-journal-cn_earthquake-research-shanxi_thesis/0201254234742.html 张淑亮, 李冬梅, 范雪芳, 2005.井水位前驱波与气压、风、降雨及强震关系的分析.地震, 25(3):69-77. http://d.old.wanfangdata.com.cn/Periodical/diz200503010 张昭栋, 张广诚, 1987.利用水位阶变资料反演震时应力的调整变化.地震研究, 10(6):693-702. http://www.cnki.com.cn/Article/CJFDTotal-DZYJ198706004.htm 张昭栋, 郑金涵, 冯初刚, 1989.水井水位的气压效率和降水荷载效率之间的定量关系.地震, (6):38-44. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK000000213971 张昭栋, 郑金涵, 冯初刚等, 1990.深井水位的地表荷载效应及其与含水层物理力学参数关系初探.地震学刊, (1):27-32. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK000000378638 张昭栋, 郑金涵, 冯初刚等, 1993.根据井水位阶变反演含水层应力变化的三种方法的对比.华北地震科学, 11(1):39-44. https://www.ixueshu.com/document/946e317b94c487682ce096f8b8e9acd0318947a18e7f9386.html 朱艾斓, 徐锡伟, 胡平等, 2005.首都圈地区b值随震源深度的变化:对地震成核的意义.科学通报, 50(8):788-792. doi: 10.3321/j.issn:0023-074X.2005.08.012 朱艾斓, 解朝娣, 徐锡伟等, 2010.鄂尔多斯块体周缘地区近期地震活动性与汶川地震应力触发作用的关系.地学前缘, 17(5):206-214. http://d.old.wanfangdata.com.cn/Periodical/dxqy201005018 Bredehoeft J. D., 1967.Response of well-aquifer systems to Earth tides. Journal of Geophysical Research(72):3075-3087. doi: 10.1029-JZ072i012p03075/