Effect Earthquake Response of Seismic Underground Fluid from the Zhaogezhuang Observation Well
-
摘要: 赵各庄井位于首都圈地区,是夏垫断裂带北端的地震观测井。其地下流体同震响应的统计数据显示,在对6级以上远场地震的响应次数上,赵各庄井的水位较水温更显著;在对MS 7.0井水位的响应幅度变化范围为3.0mm—770.0mm,典型水温的响应幅度为0.0129℃,井水位响应幅度明显强于水温,响应形态以振荡型为主。赵各庄井水温、水位的异常变化和夏垫断裂带的活动性密切相关,同时也会影响断裂带的活动性。综合分析认为,在对赵各庄井水温和水位两大测项进行观测时,应以水位为主,还应关注夏垫断裂带的活动性。同时,为了监测夏垫断裂带的活动性,应对赵各庄井水位和水温进行长期监测,以保证首都圈地区的人民生产安全。Abstract: The Zhaogezhuang well is loacated in the capital zone and it's an earthquake observation well along the Xiadian Fault zone. In this poper we conducted sfatistical study of the co-seismic corresponding of and we fand that respcmse of groundwater level is more remarkable than that of temperature for MS ≥ 6.0 earthquakes on responding times and MS ≥ 7.0 earthquakes on responding minute as well; Responding amplitude ranges from 3.0-770.0mm for groundwater level and 0.0129℃ for typical groundwater temperature, and the variation of groundwater level is more obvious than that of temperature. Oscillation mode is the main responsive form. The anomaly of groundwater level and temperature is related to the activity of the Xiadian fault zone in consideration of its location and maybe impact the fault activity. In conclusion, groundwater level of the Zhaogezhuang well should be paid more attention in order to monitor the activity of Xiadian fault zone. Meanwhile, in order to monitor the Xiadian fault zone activity, the observation should be conducted in long-time of the guarantee of production safety of capital circle people.
-
Key words:
- Xiadian active fault /
- Zhaogezhuang well /
- Water temperature /
- Water level
-
引言
赵各庄井位于首都圈地区,该井自2001年起进行数字化水位和水温观测。武安绪等(2010)发现赵各庄井流体对汶川大地震有很好的同震响应。赵各庄井是位于夏垫断裂带北端的地震观测井(图 1),由于第四纪的活动性和曾发生的三河—平谷8级地震(烈度Ⅺ),其所处的断裂带成为首都圈地区重要的地震断层之一,引起人们广泛重视。前人在该断裂带上开展了一系列研究工作,诸如沉积环境(杨晓平等,2012;张超等,2014)、年代地层(李梁等,2011;李梁,2012)、地貌(江娃利,1999;毛昌伟等,2010;丁锐等,2014)、地球物理(邓前辉等, 2000, 2001;何付兵等,2013;田优平等,2014;段美芳等,2018)、浅层探测(张先康等,2002;贾辉等,2008;冉志杰等,2013)、人工探槽(江娃利等,2000)、钻探(徐锡伟等,2000)、遥感(王雷等,2014)和地球化学(陆丽娜等, 2016, 2018)等,并积累了大量资料。目前,该项工作中仍缺乏详细的数字化水位、水温的地震流体研究。大地震产生的应力能够瞬间改变地壳介质状态,引起含水层中地下流体的流动(Wang等,2014)。近年来,越来越多的证据表明水文地质学与地震动力学存在着密切的关系(Telesca等,2015),尤其是井水位和井水温方面(兰双双等,2011;孙小龙等, 2011, 2013;陈亮等,2012;Sun等,2012;Gu等,2013;车用太等,2014;廖欣等,2014;Shi等,2015;宋洋等,2016;向阳等, 2017a, 2017b;Kim等,2018;Petitta等,2018;Liu等,2018)。井水位和井水温的同震效应是揭示地壳介质对应力—应变过程响应的有效手段,能够反映地壳动力作用下地下介质应变和孔隙压的变化特征,从而便于深入了解地下介质的动力学过程(Roeloffs,1998;Montgomery等,2003)。
图 1 观测井构造背景示意图(据高战武,2001)Figure 1. Tectonic background diagram of the observation well为了更好地监测和研究夏垫断裂带的地震活动性以及赵各庄井的同震变化情况,本文以地震地下流体两大测项——水温和水位为研究对象,针对赵各庄井水位和水温的同震响应进行分析,同时探讨赵各庄井孔水位变化对夏垫断裂带的影响。
1. 研究背景
1.1 观测井基本信息及地质概况
赵各庄井位于北京市平谷区城关镇赵各庄村内,地理坐标为117.01°E、40.13°N,海拔26.2m,属于平谷区地震办公室管理的地方台观测井。该井位处夏垫断裂带北端西侧。夏垫断裂带位于华北平原区北部,为一东南盘下降、西北盘上升的正断兼具右旋走滑分量的断裂带,其总体走向NE45°,倾向SE,倾角50°—70°,全长45km,新生代断距大(可达3—4km)。该断裂被认为是1679年三河—平谷MS 8.0地震的重要发震断裂(图 1)(韩晓昆等,2013;万永魁等,2014)。断裂带上盘为黄褐色黏土、黄色粉砂、灰黄色亚黏土以及灰黑色黏土,可见零星钙核;下盘则为灰黑色黏土、灰褐色亚黏土、灰色粉砂、灰褐色亚黏土、黄褐色黏土以及黄色亚砂土(江娃利等,2000)。
赵各庄井井孔柱状图如图 2所示,井区周围地表多被第四系松散层覆盖,观测含水层属第四系孔隙承压水,主要含水层有:①第四系卵石、砾石混黏砂含水层,位于井下171.00—254.00m;②第四系石英砂岩层,位于井下556.81—557.26m。其成井后下设无缝钢管至182.10m,其下安装同径滤水管至234.60m,井的观测含水性为182.10—25400m上层冷水和556.81m以下层热水的混合水。该观测井257.59—479.00m为卵石、砾石混黏岩,479.00—556.81m为卵石、漂石混黏岩,透水性较好。
图 2 赵各庄测井井孔柱状图(据韩孔艳等,2016)Figure 2. The stratigraphy of Zhaogezhuang wellhole赵各庄井建于1989年,井深557.26m,第四系厚度达556m,下伏蓟县系砂岩。该自2001年9月开始进行数字化水位、水温观测,使用LN—3型水位仪和SZW—1A型水温仪(表 1)。温度探头放置在180m处,且自安装以来位置基本不变。2011年5月水温仪增加协议转换器,并将LN—3型水位仪更换为SWY—1A型水位仪进行观测,2015年更换为SWY—2型水位仪。多年的观测资料显示该井对于北京地区4级以上地震有较好的映震能力(韩孔艳等,2016)。
表 1 观测仪器设备参数Table 1. Parameters of observation instrument仪器名称 型号 生产厂家 检测线 观测精度 分辨率 探头埋深 测项背景值 水位仪 LN-3 中国地震局分析预报中心 传感器 ±0.2%FS 优于1mm 20.6m 13.7m 数采 ±(0.05%R±2) 0.1mv 水温仪 SZW-1A 地壳所 0—100℃ ±0.05℃ ±0.0001℃ 180.0m 16.46℃ 注:数据源引自北京市地震局(2006)。 1.2 观测资料概况
赵各庄井自1990年投入使用以来,具有多年的水位、水温观测资料,观测精度和观测方式均为分钟值。本次研究水温观测资料数据时段为2006.01.01—2015.07.10,日均值资料缺数51天,连续性较好;水位观测资料选取数据时段为2003.01.01—2015.07.10,日均值资料缺数77天,资料连续性较好,水位数据较水温数据突变少,稳定性也强于水温数据。
2. 井水位、水温动态协调变化特征
选取的静水位和水温数据起止时间为2006.01.01—2015.07.10。刘耀炜等(2010)研究认为井孔水位、水温的协调变化特征主要有同步同向协调变化、同步反向协调变化及无规律变化3种类型。2012年之前,赵各庄井水位和水温在形态上总体呈同步反向协调变化,即水位下降,而水温上升(图 3(a));2012年之后,两者呈同步缓升的趋势(图 3(c))。井水位日值年动态曲线反映出每年冬季及次年春季静水位呈现出下降的趋势,而在夏季及秋季则出现逐渐上升。水温的变化规律则表现出2种情况:2012年之前,每年冬季及次年春季呈上升趋势;2012年之后,每年冬季及次年春季呈下降趋势。
从统计期内水位、水温旬均值关系(图 3(b)、(d))不难看出,水位和水温之间的相关系数在2012年之前为0.54494,2012年之后为0.54872,两者相关性较明显,这表明地震孕育过程对两者的作用程度是有差别的。
3. 井水位、水温同震响应特征
自Cooper等(1965)对开口井孔中水位波动对地震波的响应研究和Mogi等(1989)对日本伊豆半岛某温泉水温突变对于强震的响应研究以来,国内外开展了对水温和水位同震响应的大量研究。研究充分表明,地震流体观测中水位和水温对远场大震的同震响应是地震波作用于井—含水层系统最直接的体现(孙小龙等,2008a)。井水位、水温的同震响应特征是指井水位和水温在地震波的作用下所表现出来的异常升降(宋洋等,2016)。为提高研究的准确性,选取可能引起井水位和水温同震响应的所有地震目录(表 2),其中井水位为静水位,地震目录来源于国家地震科学数据共享中心1(刘瑞丰等,2007;),具体选取标准为:70°E以东,10°S—60°N地区的MS 4.0以上地震(震级统一用MS 标度,无MS 标度时利用经验公式Mb=0.63MS +2.5(Gutenberg等,1956;陈运泰等,2004)和MS =1.13ML-1.08(郭履灿,1971)进行换算)。以发震时刻为基点,甄别井水位和水温动态数据及同震响应的地震,分析井水位响应类型,计算震中距,筛选及计算结果见表 3。
表 2 赵各庄井地震及响应次数统计Table 2. Numbers of earthquakes and response times of observation well测项 数据时段 次数 MS 4.0—4.9 5.0—5.9 6.0—6.9 7.0—7.9 ≥8.0 静水位 2003-01-01—2015-07-10 发震 2 182 542 91 12 响应 0 0 6 13 4 水温 2006-01-01—2015-07-10 发震 1 126 385 72 9 响应 0 0 0 0 1 表 3 赵各庄井水位及水温同震响应特征Table 3. The co-seismic responsive characteristic of groundwater level and temperature of Zhaogezhuang well发震日期 发震时刻 水位开始变化时间 北纬/° 东经/° MS 震源深度
/km震中距/km 参考地点 水位变化幅度
/m同震响应形态类型 水温变化幅度/
℃滞后时间
/min2014-02-12 17:19:48 17:37 36.14 82.51 7.3 10 3031 新疆于田 0.0150 振荡—脉冲型 无 18 2008-05-12 14:27:59 14:37 31.01 103.42 8.0 14 1590 四川汶川县 0.2075 脉冲型 0.0129 9 2008-03-21 06:32:59 06:54 35.64 81.54 7.3 21 3131 新疆于田县 0.0240 振荡型 无 21 2008-01-09 16:26:44 16:44 32.39 85.27 6.6 10 2955 西藏改则县 0.0030 振荡型 无 18 2007-04-02 04:39:55 05:12 -8.5 157.00 7.7 10 6788 所罗门群岛 0.0535 振荡—脉冲型 无 33 2006-12-26 20:26:19 20:46 21.86 120.6 7.3 10 2061 中国南海 0.0150 振荡型 无 20 2006-12-26 20:34:10 20:46 21.87 120.73 6.8 10 2060 中国南海 0.0090 振荡型 无 12 2006-04-21 07:24:59 07:44 60.97 167.48 7.9 32 4103 勘查加半岛 0.0590 振荡型 无 19 2005-07-24 23:42:04 00:03 7.83 92.19 7.4 19 4355 尼科巴群岛 0.0280 振荡型 缺数据 21 2005-03-29 00:09:35 00:31 2.03 97.05 8.4 34 4690 印尼苏门答腊 0.2335 振荡型 缺数据 22 2004-12-26 08:58:51 09:25 3.15 95.79 8.8 40 4633 印尼苏门答腊 0.7700 振荡型 缺数据 27 2004-11-29 02:32:13 02:42 43.08 145.14 6.9 41 2351 日本北海道 0.0220 脉冲型 缺数据 10 2004-11-12 05:26:31 05:50 -9.16 125.41 7.2 10 5549 印尼帝汶 0.0620 振荡型 缺数据 24 2004-10-23 16:55:58 17:05 37.10 139.03 6.7 28 1937 日本 0.0050 振荡型 缺数据 9 2004-10-15 12:08:49 12:18 24.51 122.83 6.3 107 1820 中国台湾宜兰以东海中 0.0070 脉冲型 缺数据 10 2004-09-05 18:07:08 18:17 33.02 136.94 7.0 35 1942 日本本州 0.0140 脉冲型 缺数据 10 2004-09-05 22:57:16 23:08 33.15 137.16 7.3 15 1952 日本本州南 0.0230 振荡型 缺数据 11 2004-02-07 10:42:33 11:07 -3.88 135.31 7.3 11 5238 印尼伊里安岛 0.0125 振荡型 缺数据 25 2003-12-10 12:38:11 12:57 23.01 121.51 6.7 13 1950 中国台湾台东近海 0.0095 脉冲型 缺数据 19 2003-11-17 14:43:06 15:05 51.23 178.57 7.5 30 4789 拉特群岛 0.0605 脉冲型 缺数据 22 2003-10-01 09:03:25 09:18 50.13 87.82 7.2 9 2524 俄、蒙、中交界 0.0065 振荡型 缺数据 15 2003-09-27 19:33:25 19:46 50.02 87.87 7.7 16 2517 中蒙边界 0.0550 振荡型 缺数据 13 2003-09-26 03:50:05 04:00 41.79 143.71 8.2 16 2241 日本北海道 0.5455 振荡—脉冲型 缺数据 10 3.1 同震响应地震次数及分布
对于赵各庄井静水位,选取地震829次,其中具有同震响应的地震次数为23次,占全部地震的2.77%;对于水温,选取地震593次,呈现同震响应的地震为1次,占全部地震的0.16%。静水位与水温对汶川MS 8.0地震均有显著且一致的同震响应,如图 4所示,由图可见当汶川地震发震时,水位呈现突降,同时水温也呈现突降,二者表现一致。而对于其它地震,二者较少表现出明显的一致性。对于6级以上的远场地震,水位较水温有更显著的响应,对于5级及以下的地震,二者响应能力相当。不同地震发震地点不同,震级不同,释放能量不同,从而造成同一观测井对于不同地震的响应差异(巩浩波等,2015)。对于5级及以下地震,由于震级较小,地震释放的能量较弱,当地震波传播到测井时,能量削弱,不足以引起水位的变化,即未响应;而5级以上地震,震级较大,释放的能量较强,综合考虑震中距的大小,以及地震波传播过程中所经过的地质构造和岩性的不同,会引起测井水位显著变化或无明显变化,即不全有响应。
3.2 响应时间
从响应特征筛选和计算结果(表 3)可知,赵各庄井静水位同震响应时间最短的是2008年汶川MS 8.0地震和2004年10月23日日本MS 6.7地震,震后9分钟井水位即响应;而响应时间最长的2007年4月2日所罗门群岛MS 7.7地震,震后33分钟后才响应。2004年9月5日,日本本州发生了MS 7.0和MS 7.3地震,两者相隔时间为290分钟,分别在震后10分钟和11分钟出现了响应。而2006年12月26日,中国南海发生了2次地震,震级分别是MS 7.3和MS 6.8,水位的响应时间分别为20分钟和12分钟。整体来看,对于7—8级地震水位具有显著的响应能力,一般响应时间为10—20分钟。根据震中距与响应时间关系拟合图(图 5),获得拟合曲线斜率为3.28,推测在水位同震响应过程中地震面波起到了主要作用(牛安福等,2005)。此外,由图 5可以看出,震中距越大响应时间越长,震中距越小响应时间越短,当震中距小于3000km时,曲线的拟合性最好,数据最集中,说明在此震中距内,水位的同震响应更规律,效果也更好。当震中距相同或相近时,如2004年12月26日苏门答腊MS 8.8地震和2005年3月29日苏门答腊MS 8.4地震,2次地震的井水位响应时间也基本一致。
3.3 响应幅度
赵各庄井水位对各地震的同震响应幅度不同,持续时间也不等。其中,2004年12月26日印尼苏门答腊MS 8.7地震响应幅度最大,震后水位振荡幅度达770.0mm,2008年1月9日西藏改则县MS 6.9地震响应幅度最小,振荡幅度为3.0mm。此外,具有同震响应特征的地震震级多集中在MS 6.0—8.0之间(图 6(a))。
(1)2003年以来,日本附近发生的具有同震响应的地震共有5次,发震地点与赵各庄井的距离比较相近,震级分别为MS 6.9、MS 6.7、MS 7.0、MS 7.3和MS 8.2,其井水位响应幅度也不一致。震级最大的地震为2003年9月26日MS 8.2地震,其井水位响应幅度为545.5mm;震级较小的2次地震为2004年10月23日MS 6.7地震和2004年12月29日MS 6.9地震,其井水位响应幅度分别为5.0mm和22.0mm(表 3)。由此可见,当震中距相近时,震级越大,相应的井水位响应幅度也越大。
(2)2003年以来,统计范围内共发生12次MS 8.0级以上地震,赵各庄井对其中的4次地震有明显的响应(表 3):①2008年汶川MS 8.0地震距赵各庄井1590km,震源深度14km,在这12次地震中距赵各庄井震中距最小,井水位响应幅度为207.5mm;②2004年12月26日印尼苏门答腊MS 8.8地震距赵各庄井4633km,震源深度为40km,其井水位响应幅度最大,为770.0mm。
(3)2008年汶川MS 8.0地震和日本附近的多次地震,其发震地点与赵各庄井的距离在1500—2500km之间。然而响应幅度差异较大,汶川MS 8.0地震引起的水位响应幅度为207.5mm,日本附近的地震按表 3中的顺序响应幅度依次为22.0mm、5.0mm、14.0mm、23.0mm、545.5mm,对应的震级依次为MS 6.9、MS 6.7、MS 7.0、MS 7.3、MS 8.2。可以看出,在震中距相近的情况下,震级越大,同震响应幅度越大。
(4)具有响应的中国大陆型地震共计5次(表 3),震级分别为MS 7.3、MS 8.0、MS 7.3、MS 6.6、MS 7.7,其响应幅度依次为15.0mm、207.5mm、24.0mm、3.0mm、55.0mm;国外非大陆型地震共计7次,震级分别为MS 7.7、MS 7.9、MS 7.4、MS 8.4、MS 8.8、MS 7.3、MS 7.5,其响应幅度依次为53.5mm、59.0mm、28.0mm、233.5mm、770.0mm、12.5mm、60.5mm。国外非大陆型远场大震引起赵各庄井水位响应幅度明显高于中国大陆型地震,特别是2004年12月26日印尼苏门答腊MS 8.8地震其引起的水位响应幅度高达770.0mm。
(5)2008年7月5日鄂霍次克海MS 7.8地震和2015年5月30日日本志摩市MS 8.3地震,其震中距较小,赵各庄井水位无同震响应,可能是震源较深的原因,鄂霍次克海地震和日本志摩市地震震源深度分别为636km和680km。此外,通过统计也发现,当震级相同时,水位响应幅度并未随震中距场大呈减小趋势,而是不一而同,对于日本北海道MS 8.0地震与汶川MS 8.0地震,前者震中距大于后者,两者震源深度都为33km,但井水位响应幅度前者却比后者大。
对井水温同震响应统计(图 6(b))可以看出,赵各庄井仅有1次同震响应的记录,其响应幅度明显弱于井水位。有响应记录的地震为2008年5月12日的汶川地震,水温响应幅度为0.0129℃。而对于其它地震,部分没有数据记录,部分没有响应记录。整体而言,赵各庄井水位响应能力强于水温。
3.4 响应形态
响应特征统计结果(表 3)和典型水位同震响应曲线(图 7)显示,赵各庄井水位同震响应形态为振荡型、振荡—脉冲型和脉冲型3种。其中,图 7(a)、7(b)、7(c)和7(h)中的水位变化为振荡型,图 7(e)和7(f)为典型的脉冲型,图 7(d)和7(g)则是兼有振荡和脉冲变化的振荡—脉冲型。图 7(e)和7(f)中的水位响应形态虽同为脉冲型,但图 7(e)为向上脉冲型,图 7(f)则为向下脉冲型,二者皆有显著的脉冲峰。振荡型也分为等振幅和不等振幅2种类型,其中图 7(b)是典型的等幅振荡型。振荡型为水位响应形态的主要类型。图 7(a)、7(b)、7(c)和7(f)的水位变化幅度分别为770.0mm、233.5mm、545.5mm和207.5mm,变化幅度最显著,而其响应时间分别是27min、22min、10min和9min,结合这些台站到地震震中的震中距,分析认为该井水位变化对于这几次地震的响应最快。
4. 异常成因探讨
4.1 水位变化与震中距、震级之间的关系
由震中距、震级和井水位实测变化值,得出三者之间的空间分布情况(图 8)。据前人研究(Roeloffs,1998;杨竹转等,2005;陈大庆等, 2007a, 2007b),水位的变化幅度与地震震级、震中距之间存在一定的数量关系:
$${\rm{lg}}\Delta {h_i} = {\omega _1}M + {\omega _2}{\rm{lg}}D + {\omega _3} $$ (1) 式中:ω1、ω2、ω3为常数;Δhi为水位变化幅度(cm);D为震中距(km);M为MS ,表示面波震级。
Roeloffs(1998)通过对美国加利福尼亚洲BV井的研究得出以下关系:
$$ {\rm{lg}}\Delta h = 0.89M - 1.63{\rm{lg}}D - 0.69 $$ (2) 当Δh≥5cm时,有;
$$ M \ge 1.55 + 1.82{\rm{lg}}D $$ (3) 即当震级与震中距满足(3)式时,会产生超过5cm的水位上升。
针对赵各庄井,根据式(1),不考虑水位上升或下降的问题,对所记录到的水位变化幅度与震级、震中距(表 3)之间做二元回归计算分析,得到以下关系式:
$$ {\rm{lg}}\Delta {h_i} = 1.923M + 1.084{\rm{lg}}D - 19.350 $$ (4) 当Δh≥0.7cm时,有:
$$ M \ge 9.98 - 0.563{\rm{lg}}D $$ (5) 即当震级与震中距满足式(5)时,井水位会产生超过0.7cm的水位变化。
经过计算,在所选震例中共有4次地震(MS 7.7、MS 8.2、MS 8.4、MS 8.8)符合式(5),4次地震的实测震级大于预测震级。其它地震的实测震级均小于预测震级,且在这些地震中,某些地震震中距很大,反而引起了大于0.7cm的水位变幅,因此,赵各庄井的井水位变化显然不受震级与震中距的严格约束。根据式(4),绘制出震级、震中距与水位变化幅度关系的三维拟合图(图 8),从中可以看出引起井水位变化的地震分布情况,同时反映出震级在7级以上的地震其水位变化幅度最为显著,震中距对水位变化的影响无明显规律。
4.2 井水位变幅与含水层岩性的关系
赵各庄井周围被大量第四纪堆积物所覆盖(图 2),且堆积厚度大,多为卵石、砾石混黏岩,其最深处为第四系石英砂岩层。根据史浙明(2015)对岩性与同震水位响应之间关系的研究可知,岩性为火成岩、砂岩以及灰岩并不是引起井水位变幅的主要因素。而赵各庄井并非这种情况,其第四纪堆积物并未成岩,孔隙较大,导水性较好,当远场中强震的地震波通过这种介质时,产生的动态应力会使得含水层渗透性进一步增大或减小,渗透性产生的变化导致水位发生振荡,这与响应地震的相应形态类型基本相符。
5. 讨论
5.1 井孔水位、水温映震效应特征
井孔水位、水温对于地震的映震响应是客观存在的,其形态类型复杂多样,但并非无规律可循。根据前人研究,对于同震响应的水位变化,近场地震以同震阶变为主(付虹等,2002;王芳等,2012),中场地震多出现缓变,而水位振荡则主要出现于远场大震(陈大庆等,2007a;尹宏伟等,2016);而水温对于远场地震,其特征表现为震时水温突降—缓慢恢复(陈大庆等,2007a;孙小龙等,2008b)。通常震前以阶变类异常为主(车用太等,1996),两者之间常出现同震水位振荡且水温下降的特征(石耀霖等,2007)。这些同震响应的方式除与地震特征有关外,还主要取决于观测井局部的地质构造、水文地质条件以及观测井周边的环境等。如在含水层埋深大致相同时,水位同震响应的幅度取决于观测井含水层的岩性,一般灰岩 > 变质岩 > 砂岩 > 第四系砂砾石,含水层顶板埋深亦对其有影响。
赵各庄井水位对一些远场大震同震响应形态以振荡为主,同时也出现脉冲、振荡—脉冲的形态。对于2008年汶川MS 8.0地震,其震后9分钟出现阶降,而后水位并未恢复到原来的水位高度,但有所上升,达207.5mm。
对于远场大震,其震后响应时间为几分钟到十几分钟,有时达到20分钟以上,当地震波传到观测井后,会导致观测井—含水层系统介质发生一定形变,从而促使其孔隙压力发生改变,进而引起水流速度发生变化,当孔隙压力增大时,含水层中水会向井孔流动,井孔水位因此上升;当孔隙压力减小时,井孔中水会向含水层流动,导致井孔水位下降。地震波的这种持续影响,导致观测井孔水位持续上升或持续下降,震后缓慢恢复或者难以恢复;井孔水位上升与下降交替转变,持续一段时间后停止。
5.2 井孔水位变化与断裂
地下流体观测井井孔均尽可能地布设在断层及其邻近地区,这些观测井的水位、水温及一些化学量的异常变化都有可能与观测井含水层相关断层的活动性有直接的关系(王博等,2008)。地下水位的异常变化除了受地壳形变的影响外,还主要受降雨的影响,并且降雨也可以导致地壳的形变(陆明勇等,2005)。降雨可以在地表形成地表积水,使得地面的荷载作用加强,其下岩土中的很大范围内产生较大的附加应力。此外,降雨也能增加地下含水层中的含水量,使得孔隙压力增大,导致含水层中的水向井孔流动,水位上升;反之,可导致水位下降,水位的这种变化将产生动水压力。在地震孕育过程中的中短、短临阶段,水位上升或下降,都会对断层活动性产生一定的影响。此外,黄辅琼等(2005)通过研究断层的变形行为与降雨及地下水的关系,分析认为断层的年变形量与年降雨量相关,较少或者较多的不规则降雨都会导致断层的变形行为发生畸变,降雨通过改变断层带上的地下水位(或断层的孔隙压力)影响着其断层的变形行为,浅部断层的流体可以通过影响断层的变形行为而参与或者影响区域构造应力应变的调整。因此,为监测夏垫断裂带的活动性,应对赵各庄井水位进行长期监测,尤其是降雨量发生突变的季节,更应重点监测其前后变化。
对于布设在夏垫断裂带上的观测井——赵各庄井,其井水位下降,可反应该地区地下水位整体呈下降趋势。在下降初期,虽然含水层孔隙压力降低,但是水位变化引起的动水压力却增强,依然会增强断层面上的总压力及摩擦力,故在下降的初期会导致断层的滑动。但是随着水位下降的速度变缓慢,孔隙压力继续降低,水位变化幅度变小,反而阻止断层活动。
6. 结论
地震次数和水温、水位同震响应次数的统计数据显示,对于6级以上的远场地震,赵各庄井水位较水温具有更显著的响应,对于5级及以下的地震二者响应能力相当。而在响应时间上,对于7—8级地震水位具有显著的响应能力,响应时间一般为10—20min。在响应幅度上,具有水位同震响应的地震多集中在6—8级,幅度变化范围为3.0—770.0mm;而水温的典型响应地震为2008年的汶川地震,其响应幅度为0.0129 ℃,其它地震几乎无明显响应。水位响应形态可分为振荡型、脉冲型和振荡—脉冲型,且以振荡型为主,其中变化幅度最大的4次地震,其响应时间也是最快速的。赵各庄井周围的第四纪堆积物未成岩,空隙较大,导水性较好,有利于地震波的传播,进而对水位的振荡产生影响。
除上述因素外,观测井的地质构造、水文地质条件等也决定了同震响应的特征。由于赵各庄井位于夏垫断裂带的北端,其水温、水位的异常变化和夏垫断裂带的活动性相关,且水温、水位的异常又会影响夏垫断裂带的活动性。今后对赵各庄井水温、水位两大测项进行观测时,应以水位为主,同时,还需关注夏垫断裂带的活动性。同时为了监测夏垫断裂带的活动性,应对赵各庄井水位和水温进行长期监测,尤其是降雨量发生突变的季节,更应重点监测其前后变化,以保证首都圈地区的人民生产安全。
致谢: 防灾科技学院宋洋老师、姜纪沂教授协助收集论文资料,贾建鹏老师对于数据的获取给予指导,在此表示感谢! -
图 1 观测井构造背景示意图(据高战武,2001)
Figure 1. Tectonic background diagram of the observation well
图 2 赵各庄测井井孔柱状图(据韩孔艳等,2016)
Figure 2. The stratigraphy of Zhaogezhuang wellhole
表 1 观测仪器设备参数
Table 1. Parameters of observation instrument
仪器名称 型号 生产厂家 检测线 观测精度 分辨率 探头埋深 测项背景值 水位仪 LN-3 中国地震局分析预报中心 传感器 ±0.2%FS 优于1mm 20.6m 13.7m 数采 ±(0.05%R±2) 0.1mv 水温仪 SZW-1A 地壳所 0—100℃ ±0.05℃ ±0.0001℃ 180.0m 16.46℃ 注:数据源引自北京市地震局(2006)。 表 2 赵各庄井地震及响应次数统计
Table 2. Numbers of earthquakes and response times of observation well
测项 数据时段 次数 MS 4.0—4.9 5.0—5.9 6.0—6.9 7.0—7.9 ≥8.0 静水位 2003-01-01—2015-07-10 发震 2 182 542 91 12 响应 0 0 6 13 4 水温 2006-01-01—2015-07-10 发震 1 126 385 72 9 响应 0 0 0 0 1 表 3 赵各庄井水位及水温同震响应特征
Table 3. The co-seismic responsive characteristic of groundwater level and temperature of Zhaogezhuang well
发震日期 发震时刻 水位开始变化时间 北纬/° 东经/° MS 震源深度
/km震中距/km 参考地点 水位变化幅度
/m同震响应形态类型 水温变化幅度/
℃滞后时间
/min2014-02-12 17:19:48 17:37 36.14 82.51 7.3 10 3031 新疆于田 0.0150 振荡—脉冲型 无 18 2008-05-12 14:27:59 14:37 31.01 103.42 8.0 14 1590 四川汶川县 0.2075 脉冲型 0.0129 9 2008-03-21 06:32:59 06:54 35.64 81.54 7.3 21 3131 新疆于田县 0.0240 振荡型 无 21 2008-01-09 16:26:44 16:44 32.39 85.27 6.6 10 2955 西藏改则县 0.0030 振荡型 无 18 2007-04-02 04:39:55 05:12 -8.5 157.00 7.7 10 6788 所罗门群岛 0.0535 振荡—脉冲型 无 33 2006-12-26 20:26:19 20:46 21.86 120.6 7.3 10 2061 中国南海 0.0150 振荡型 无 20 2006-12-26 20:34:10 20:46 21.87 120.73 6.8 10 2060 中国南海 0.0090 振荡型 无 12 2006-04-21 07:24:59 07:44 60.97 167.48 7.9 32 4103 勘查加半岛 0.0590 振荡型 无 19 2005-07-24 23:42:04 00:03 7.83 92.19 7.4 19 4355 尼科巴群岛 0.0280 振荡型 缺数据 21 2005-03-29 00:09:35 00:31 2.03 97.05 8.4 34 4690 印尼苏门答腊 0.2335 振荡型 缺数据 22 2004-12-26 08:58:51 09:25 3.15 95.79 8.8 40 4633 印尼苏门答腊 0.7700 振荡型 缺数据 27 2004-11-29 02:32:13 02:42 43.08 145.14 6.9 41 2351 日本北海道 0.0220 脉冲型 缺数据 10 2004-11-12 05:26:31 05:50 -9.16 125.41 7.2 10 5549 印尼帝汶 0.0620 振荡型 缺数据 24 2004-10-23 16:55:58 17:05 37.10 139.03 6.7 28 1937 日本 0.0050 振荡型 缺数据 9 2004-10-15 12:08:49 12:18 24.51 122.83 6.3 107 1820 中国台湾宜兰以东海中 0.0070 脉冲型 缺数据 10 2004-09-05 18:07:08 18:17 33.02 136.94 7.0 35 1942 日本本州 0.0140 脉冲型 缺数据 10 2004-09-05 22:57:16 23:08 33.15 137.16 7.3 15 1952 日本本州南 0.0230 振荡型 缺数据 11 2004-02-07 10:42:33 11:07 -3.88 135.31 7.3 11 5238 印尼伊里安岛 0.0125 振荡型 缺数据 25 2003-12-10 12:38:11 12:57 23.01 121.51 6.7 13 1950 中国台湾台东近海 0.0095 脉冲型 缺数据 19 2003-11-17 14:43:06 15:05 51.23 178.57 7.5 30 4789 拉特群岛 0.0605 脉冲型 缺数据 22 2003-10-01 09:03:25 09:18 50.13 87.82 7.2 9 2524 俄、蒙、中交界 0.0065 振荡型 缺数据 15 2003-09-27 19:33:25 19:46 50.02 87.87 7.7 16 2517 中蒙边界 0.0550 振荡型 缺数据 13 2003-09-26 03:50:05 04:00 41.79 143.71 8.2 16 2241 日本北海道 0.5455 振荡—脉冲型 缺数据 10 -
北京市地震局, 2006.北京市地震监测志.北京:地震出版社. 车用太, 鱼金子, 刘春国, 1996.我国地震地下水温度动态观测与研究.水文地质工程地质, (4):34-37. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK199600494353 车用太, 何案华, 鱼金子, 2014.水温微动态形成的水热动力学与地热动力学机制.地震学报, 36(1):106-117. doi: 10.3969/j.issn.0253-3782.2014.01.009 陈大庆, 刘耀炜, 杨选辉等, 2007a.远场大震的水位、水温同震响应及其机理研究.地震地质, 29(1):122-132. http://d.old.wanfangdata.com.cn/Periodical/dzdz200701011 陈大庆, 万永芳, 赵明利, 2007b.北京良乡、赵各庄井水位、水温同震响应对比研究.华南地震, 27(2):69-75. http://d.old.wanfangdata.com.cn/Periodical/hndz200702010 陈亮, 丁富雄, 张蕾等, 2012.建水地下流体监测井水位骤降异常分析.地震研究, 35(4):471-476. doi: 10.3969/j.issn.1000-0666.2012.04.005 陈运泰, 刘瑞丰, 2004.地震的震级.地震地磁观测与研究, 25(6):1-12. doi: 10.3969/j.issn.1003-3246.2004.06.001 邓前辉, 王继军, 汤吉等, 2000.三河-平谷8级大震区地壳上地幔电性结构特征.见: 中国地震学会第八次学术大会论文摘要集.北京: 地震学报出版社, 282. http://cpfd.cnki.com.cn/Article/CPFDTOTAL-ZGEM200000001282.htm 邓前辉, 王继军, 汤吉等, 2001.三河-平谷8级大震区地壳上地幔电性结构特征研究.地震地质, 23(2):178-185. doi: 10.3969/j.issn.0253-4967.2001.02.006 丁锐, 张世民, 毛昌伟等, 2014.夏垫断裂断层陡坎的微地貌分析.见: 2014年中国地球科学联合学术年会——专题15: 活动断层长期滑动习性、深部结构与地震论文集.北京: 中国地球物理学会, 全国岩石学与地球动力学研讨会组委会, 中国地质学会构造地质学与地球动力学专业委员会, 中国地质学会区域地质与成矿专业委员会. http://cpfd.cnki.com.cn/Article/CPFDTOTAL-ZGDW201410020016.htm 段美芳, 陆丽娜, 李莹甄等, 2018.张家口-渤海断裂带西段及中西段b值时空扫描.震灾防御技术, 13(1):138-148. http://zzfy.eq-j.cn/zzfyjs/ch/reader/view_abstract.aspx?flag=1&file_no=20180112&journal_id=zzfyjs 付虹, 刘丽芳, 王世芹等, 2002.地方震及近震地下水同震震后效应研究.地震, 22(4):55-66. http://d.old.wanfangdata.com.cn/Periodical/diz200204010 高战武, 2001.张家口-蓬莱断裂带地震地质特征研究.北京: 中国地震局地质研究所. 巩浩波, 郭卫英, 李光科等, 2015.重庆井网水位水温同震响应特征分析.震灾防御技术, 10(S1):794-804. doi: 10.11899/zzfy2015s111 郭履灿, 1971.华北地区的地方性震级ML和面波震级MS 经验关系.见: 全国地震工作会议.1-10. 韩孔艳, 赵健, 2016.赵各庄地震观测井观测现状分析及改造方案探讨.煤炭与化工, 39(7):129-132. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=hebhg201607043 韩晓昆, 李营, 杜建国等, 2013.夏垫断裂中南段土壤气体地球化学特征.物探与化探, 37(6):976-982. http://d.old.wanfangdata.com.cn/Periodical/wtyht201306003 何付兵, 白凌燕, 王继明等, 2013.夏垫断裂带深部构造特征与第四纪活动性讨论.地震地质, 35(3):490-505. doi: 10.3969/j.issn.0253-4967.2013.03.004 黄辅琼, 陈顒, 白长清等, 2005.八宝山断层的变形行为与降雨及地下水的关系.地震学报, 27(6):637-646. doi: 10.3321/j.issn:0253-3782.2005.06.008 贾辉, 何正勤, 叶太兰, 2008.用瞬态多道瑞利波法研究夏垫隐伏断裂附近的浅层速度结构变化.地震地磁观测与研究, 29(1):15-19. doi: 10.3969/j.issn.1003-3246.2008.01.003 江娃利, 1999.北京平谷地区地表陡坎的成因识别.地震地质, 21(4):309-315. doi: 10.3969/j.issn.0253-4967.1999.04.003 江娃利, 侯治华, 肖振敏等, 2000.北京平原夏垫断裂齐心庄探槽古地震事件分析.地震地质, 22(4):413-422. doi: 10.3969/j.issn.0253-4967.2000.04.010 兰双双, 迟宝明, 姜纪沂, 2011.地下水位对近震和远震异常响应的比较——以汶川地震和苏门答腊地震为例.吉林大学学报(地球科学版), 41(1):145-152. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=cckjdxxb201101017 李梁, 张世红, 蔡向民等, 2011.北京市平原区夏垫断裂活动性的磁学研究.见: 中国地球物理学会第二十七届年会论文集.长沙: 中国地球物理学会, 161. http://cpfd.cnki.com.cn/Article/CPFDTOTAL-ZGDW201110001130.htm 李梁, 2012.北京市平原区夏垫断裂活动性的磁性地层学研究.北京: 中国地质大学(北京). http://cdmd.cnki.com.cn/Article/CDMD-11415-1012365104.htm 廖欣, 刘春平, 石云等, 2014.川06井水位固体潮效应变化初探.地震学报, 36(2):299-305. doi: 10.3969/j.issn.0253-3782.2014.02.014 刘瑞丰, 蔡晋安, 彭克银等, 2007.地震科学数据共享工程.地震, 27(2):9-16. doi: 10.3969/j.issn.1000-3274.2007.02.002 刘耀炜, 陆明勇, 付虹等, 2010.地下流体动态信息提取与强震预测技术研究.北京:地震出版社. 陆丽娜, 李静, 沈军等, 2016.夏垫活动断层土壤氡地球化学特征.震灾防御技术, 11(4):736-746. http://zzfy.eq-j.cn/zzfyjs/ch/reader/view_abstract.aspx?flag=1&file_no=20160404&journal_id=zzfyjs 陆丽娜, 杨明, 李静等, 2018.土壤气汞探测在夏垫断裂带的应用研究.地质与勘探, 54(1):112-120. doi: 10.3969/j.issn.0495-5331.2018.01.012 陆明勇, 黄辅琼, 刘善华等, 2005.地壳变形与地下水相互作用及其异常关系初探.地震, 25(1):67-73. http://d.old.wanfangdata.com.cn/Periodical/diz200501009 毛昌伟, 丁锐, 龚正等, 2010.1679年三河-平谷8级地震地表断层陡坎的GPS测量.地壳构造与地壳应力文集, (22):11-18. http://www.cnki.com.cn/Article/CJFDTotal-SEIS201000003.htm 牛安福, 张晶, 吉平, 2005.强地震引起的同震形变响应.内陆地震, 19(1):1-7. doi: 10.3969/j.issn.1001-8956.2005.01.001 冉志杰, 李皓, 吕国军等, 2013.夏垫断裂夏垫段浅部构造特征地震探测.地震工程学报, 35(3):656-663. doi: 10.3969/j.issn.1000-0844.2013.03.0651 石耀霖, 曹建玲, 马丽等, 2007.唐山井水温的同震变化及其物理解释.地震学报, 29(3):265-273. doi: 10.3321/j.issn:0253-3782.2007.03.005 史浙明, 2015.地下水位同震响应特征及机理研究.北京: 中国地质大学(北京). http://cdmd.cnki.com.cn/Article/CDMD-11415-1015391754.htm 宋洋, 谷洪彪, 李海君等, 2016.断裂两盘井水位同震响应特征对比分析——以北京八宝山断裂带中段大灰厂两观测井为例.吉林大学学报(地球科学版), 46(6):1815-1822. http://d.old.wanfangdata.com.cn/Periodical/cckjdxxb201606018 孙小龙, 刘耀炜, 2008a.苏门答腊8.5级地震引起的水温响应变化.华北地震科学, 26(1):35-40. http://d.old.wanfangdata.com.cn/Periodical/hbdzkx200801008 孙小龙, 刘耀炜, 2008b.塔院井水位和水温的同震响应特征及其机理探讨.中国地震, 24(2):105-115. http://d.old.wanfangdata.com.cn/Periodical/zgdz200802002 孙小龙, 刘耀炜, 宴锐, 2011.利用水位资料反演华北地区构造应力场变化.地震, 31(2):42-49. http://d.old.wanfangdata.com.cn/Periodical/diz201102005 孙小龙, 刘耀炜, 晏锐, 2013.云南姚安井2009年10月后水位下降的成因分析.地震学报, 35(3):410-420. doi: 10.3969/j.issn.0253-3782.2013.03.012 田优平, 余达远, 万永革等, 2014.三河-平谷地震区地球物理特征研究.地球物理学进展, 29(4):1563-1572. http://www.cnki.com.cn/Article/CJFDTotal-DQWJ201404010.htm 万永魁, 沈军, 于晓辉等, 2014.北京平原区夏垫活动断裂滑动速率及古地震复发间隔.防灾科技学院学报, 16(3):38-45. doi: 10.3969/j.issn.1673-8047.2014.03.007 王博, 刘耀炜, 孙小龙, 2008.地下流体与断裂活动关系的研究综述.地震研究, 31(3):296-302. doi: 10.3969/j.issn.1000-0666.2008.03.018 王芳, 黄辅琼, 2012.安溪一号井水位的同震阶变响应特征.地震地磁观测与研究, 33(5):169-177. doi: 10.3969/j.issn.1003-3246.2012.05/06.031 王雷, 沈军, 林玲玲等, 2014.北京平原区夏垫断裂的多源遥感影像特征.防灾科技学院学报, 16(4):33-39. doi: 10.3969/j.issn.1673-8047.2014.04.006 武安绪, 李国江, 李平安等, 2010.北京地区流体观测方法记录到的汶川大地震同震变化.地震地磁观测与研究, 31(1):71-76. doi: 10.3969/j.issn.1003-3246.2010.01.014 向阳, 孙小龙, 高小其等, 2017a.新10井水位对九寨沟MS 7.0精河MS 6.6地震同震响应.中国地震, 33(4):563-574. http://d.old.wanfangdata.com.cn/Periodical/zgdz201704012 向阳, 孙小龙, 杨朋涛, 2017b.新疆阿克陶MS 6.7地震引起的新10井水位同震响应研究.地震学报, 39(6):899-909. http://d.old.wanfangdata.com.cn/Periodical/dizhen201706008 徐锡伟, 计凤桔, 于贵华等, 2000.用钻孔地层剖面记录恢复古地震序列:河北夏垫断裂古地震研究.地震地质, 22(1):9-19. doi: 10.3969/j.issn.0253-4967.2000.01.002 杨晓平, 曹景虎, 陈献程, 2012.夏垫活动断裂两盘岩心氧化铁变化.地震地质, 34(4):659-671. doi: 10.3969/j.issn.0253-4967.2012.04.010 杨竹转, 邓志辉, 赵云旭等, 2005.云南思茅大寨井水位同震阶变的初步研究.地震学报, 27(5):569-574. doi: 10.3321/j.issn:0253-3782.2005.05.012 尹宏伟, 梁丽环, 韩文英等, 2016.黄骅井水位映震效能分析.地震地磁观测与研究, 37(3):23-29. http://d.old.wanfangdata.com.cn/Periodical/dzdcgcyyj201603004 张超, 沈军, 赵江涛等, 2014.北京平原夏垫断裂沉积旋回研究.防灾科技学院学报, 16(2):49-57. doi: 10.3969/j.issn.1673-8047.2014.02.009 张先康, 赵金仁, 刘国华等, 2002.三河-平谷8.0级大震区震源细结构的深地震反射探测研究.中国地震, 18(4):326-336. doi: 10.3969/j.issn.1001-4683.2002.04.002 Cooper Jr. H. H., Bredehoeft J. D., Papadopulos I. S., et al., 1965. The response of well-aquifer systems to seismic waves. Journal of Geophysical Research, 70 (16):3915-3926. doi: 10.1029/JZ070i016p03915 Gu S. Y., Liu B. W., Zhang H., et al., 2013. Numerical simulation of coseismic response and its mechanism of well water temperature to farfield strong earthquakes. Earthquake Research in China, 27 (3):316-330. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=zgdzyj-e201303004 Gutenberg B., Richter C. F., 1956. Magnitude and energy of earthquakes. Annals of Geophysics, 9 (1):1-15. http://d.old.wanfangdata.com.cn/OAPaper/oai_doaj-articles_e788e8839342702b331c76926f6e3f47 Kim S. W., Choi H. S., Park D. U., et al., 2018. Water level response measurement in a steel cylindrical liquid storage tank using image filter processing under seismic excitation. Mechanical Systems and Signal Processing, 101:274-291. doi: 10.1016/j.ymssp.2017.08.035 Liu C. Y., Chia Y. P., Chuang P. Y., et al., 2018. Impacts of hydrogeological characteristics on groundwater-level changes induced by earthquakes. Hydrogeology Journal, 26 (2):451-465. http://www.springerlink.com/openurl.asp?id=doi:10.1007/s10040-017-1684-z Mogi K., Mochizuki H., Kurokawa Y., 1989. Temperature changes in an artesian spring at Usami in the Izu Peninsula (Japan) and their relation to earthquakes. Tectonophysics, 159 (1-2):95-108. doi: 10.1016/0040-1951(89)90172-8 Montgomery D. R., Manga M., 2003. Streamflow and water well responses to earthquakes. Science, 300 (5628):2047-2049. doi: 10.1126/science.1082980 Petitta M., Mastrorillo L., Preziosi E., et al., 2018. Water-table and discharge changes associated with the 2016-2017 seismic sequence in central Italy:hydrogeological data and a conceptual model for fractured carbonate aquifers. Hydrogeology Journal, 3 (4):1009-1026. Roeloffs E. A., 1998. Persistent water level changes in a well near Parkfield, California, due to local and distant earthquakes. Journal of Geophysical Research, 103 (B1):869-889. doi: 10.1029/97JB02335 Shi Z. M., Wang G. C., Manga M., et al., 2015. Mechanism of co-seismic water level change following four great earthquakes-insights from co-seismic responses throughout the Chinese mainland. Earth and Planetary Science Letters, 430:66-74. doi: 10.1016/j.epsl.2015.08.012 Sun X. L., Liu Y. W., 2012. Changes in groundwater level and temperature induced by distant earthquakes. Geosciences Journal, 16 (3):327-337. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=96c834f61c6f5f38625cb44f017942a1 Telesca L., do Nascimento A. F., Bezerra F. H. R., et al., 2015. Analysis of the cross-correlation between water level and seismicity at ACU reservoir (Brazil). Tectonophysics, 658:151-158. doi: 10.1016/j.tecto.2015.07.017 Wang C. Y., Manga M., 2014. Earthquakes and water. Encyclopedia of Complexity and Systems Science. New York: Springer, 1-18. 期刊类型引用(4)
1. 陆丽娜,俞美艳,李静,杜宁,袁有恩,贾建鹏. 川05井动水位对典型地震的映震效能分析. 地震. 2025(01): 16-31 . 百度学术
2. 刘阁,姬霄鹤,郭少峰,李志涛. 范县井水位对远场大震的同震响应特征. 华北地震科学. 2023(03): 74-79 . 百度学术
3. 也尔扎提·巴合提汉. 乌什地震台静水位观测一致性及干扰因素分析. 内陆地震. 2021(04): 399-405 . 百度学术
4. 孔凡亮,邢喜民. 巴彦岱1号井静水位周期变化主要影响因素研究. 震灾防御技术. 2020(03): 643-650 . 本站查看
其他类型引用(1)
-