王卫东,林晨,张永志,赵云峰
(1.长安大学地质工程与测绘学院,西安710054;2.浙江大学理学院,杭州310027)
从20世纪70年代发现钻孔注水诱发地震,到20世纪90年代国际上出现地壳流体研究的热点,人们越来越认识到地壳中流体的存在,并对地下流体,特别是热水在地震孕育与发生过程中的作用进行了深入的研究[1-6],国内也对地热水与地震活动的关系进行了一些研究[7-8]。
西安地区属于国家级地震重点监视防御区,历史上发生过多次强烈地震,过量开采地下承压水已引起西安市严重的地面沉降和地裂缝活动,刘快胜等发现西安市地裂缝活动和地下热水的开采有明显的相关性[9],吴富春等认为西安地区的地下热水开采没有诱发西安地区的地震[10]。我们以往的研究则表明西安地区地下热水开采和区域地震活动密切相关,西安市抽取地下热水的影响范围主要局限于西安及邻区[11-12]。本文在研究西影井水位观测资料的基础上,分析了西安市地热水活动的季节性特征,及其与西安及邻区、关中东部地区和关中地区地震活动的相关性。发现西安市地下热水活动与区域地震活动密切相关,其相关性基本不受地热水活动季节性变化的影响,表明与构造活动相关的区域应力场变化引起地下水系统的变形和岩石孔隙水压力的变化,可能导致地下热水活动的异常。
西安地区的地下热水开采有明显的季节性变化,一二季度用水量大,三四季度用水量小[13],但地下热水开采的季节性变化是否与区域地震活动相关尚缺乏研究。位于西安市南郊的西影井观测段是第三系热水,其周围的地热水开采始于1992年。因此,本研究选取西影井1992-2006年的水位日均值资料,研究地下水位季节性变化与区域地震活动的关系。将观测水位日均值资料,经气压校正后,对校正值进行一阶差分,对1-2季度和3-4季度超出背景阈值的超差次数进行了统计,得到西影井水位1-2季度和3-4季度超差次数的累计曲线图(图1-2)。
图1 1992-2006年1-2季度水位异常曲线
图2 1992-2006年3-4季度水位异常曲线
由图可以看出,该井每年的水位变幅(或超差次数)相差悬殊,具有明显的阶段性特点,可分为3个阶段:第Ⅰ阶段(1992-1993年)水位动态比较平静,第Ⅱ阶段(1994-2000年)水位动态比较活跃,第Ⅲ阶段(2000-2006年)水位动态又趋于平静;该井1-2季度和3-4季度的季节性变化不很明显,但在1997-1998年间略有不同,1997年3-4季度的超差次数明显多于1997年1-2季度的超差次数,可能与1998年1月5日泾阳5.2级地震对地下水的影响有关。
利用1992-2006年间陕西省地震目录,对关中地区(33.5°-35.5°N 、106°-111°E)、关中东部地区(33.5°-35.5°N 、108°-111°E)和西安及邻区(33.8°-34.7°N 、108.5°-109.5°E)的地震活动性进行了分析,得到了累积地震应变能年释放量时间序列图[12](图 3-5)。
图3 关中地区地震应变释放累计曲线
由图可见,1992-2006年,地震活动也具有较为明显的阶段性特征,且每阶段的持续时间、活动水平与西影井水位异常活动有较好的一致性。
图4 西安及邻区地震应变释放累计曲线
图5 关中东部地区地震应变释放累计曲线
对比图1-5可见,三个区域的地震应变释放与西影井水位异常均存在较为一致的变化趋势。为进一步研究西安市地热水活动的季节性特征与区域地震活动相关性的差异,分别对西影井1992-2006年1-2季度和3-4季度水位异常曲线与上述三个区域的地震应变释放曲线进行了相关分析,结果如下(表1)。
表1 西影井水位异常与区域地震活动的相关分析
相关分析的结果表明,西影井水位异常与西安及邻区的地震活动最为密切,与关中东部地区和关中地区的地震活动也有较明显的相关性,基本不受季节性用水的影响。究其原因,由于西影井水位所观测的是1 520 m以下深层热水含水层的孔隙压力,观测层位是具有一定封闭条件的承压系统,并能通过压力传递等方式灵敏地反映地壳中的应力应变状态,其动态反映了地壳中的应力应变状态,使岩体中微小应力应变都能反映到孔隙压力的改变上,因此它的动态是构造活动的深部反映。西影热水井水位变化与关中地区、关中东部地区和西安及邻区的地震活动均相关,反映出它们受同一构造力源的控制,而西影热水井水位变化与西安及邻区的地震活动最为密切的内在原因,可能在于它们都处于长安-临潼活动断裂的控制范围内。而季节性用水对区域构造应力场的影响很小,因而不影响西安市地下热水活动与区域地震活动的相关性。
由上述分析可知,西安地区地下热水活动与西安及邻区的地震活动最为密切,与关中东部地区和关中地区的地震活动也有较明显的相关性,基本不受季节性用水的影响,反映出它们受同一力源构造的控制;西安地区地下热水活动与西安及邻区的地震活动最为密切的原因,可能在于它们都处于长安-临潼活动断裂的影响范围内,都与区域构造活动相关,其活动水平取决于对构造活动的响应程度;研究结果也表明,与构造活动相关的区域应力场变化引起地下水系统的变形和岩石孔隙水压力的变化,可能导致了地下热水活动的异常。
[1] Asteriadis G,Livieratos E.Pre-Seismic Responses of the Underground Water Level and Temperature concerning a 4.8 Magnitude Earthquake in Greece on 20 October 1988[J].Tectonophysics,1989,170:165-169.
[2] Contadakis M E,Asteriadis G.Hydrologic changes as possible earthquake precursors in Greece[J].Natural Hazard,2001,23:29-47.
[3] Hayder J,Al-Shukri,Brian JM.Reduced seismic velocities in the source zoneof New Madrid earthquakes[J].Bull.Seism.Soc.Am,1998,78(4):1491-1509.
[4] Xu Huilong,Yasue Oki,Toshikata Ito.The 1995 Northern Niigata Earthquake with special attention to a geopressured hydrothermal system[J].The Island Arc,1998,7:647-659.
[5] King C Y,Azuma S,Igarashi G,et al.Earthquakerelated Water Level Changes at 16 closely clustered Wells in Tono,Central Japan[J].J Geophys Res.,1999,104:13073-13082.
[6] Oki Y,Hiraga S.Groundwater monitoring for earthquake by an amateur network in Japan[J].Pageoph,1988,126:211-240.
[7] 宋贯一,易立新,宋晓冰.地下热水对断裂活动与地震活动的影响研究[J].地震学报,2000,22(6):632-636.
[8] 张淑亮,李冬梅,马朝辉.山西地震带热流体对强震活动的影响[J].地震,2002,22(1):84-90.
[9] 刘快胜,张家明,米秋霞.西安地裂缝近期活动与热水井水位关系研究[J].华北地震科学,1996,14(2):67-73.
[10] 吴富春,宋立胜,朱兴国,等.西安地区的地热水开采与诱发地震关系的研究[J].地震学报,2001,23(4):407-412.
[11] 王卫东,彭建兵,张永志,等.西安市地热水开采现状及其环境问题[J].水土保持研究,2005,12(5):266-267.
[12] 王卫东,林晨,张永志,等.区域地震活动与西安市地下热水活动相关性的研究[J].地震,2009,29(2):118-122.
[13] 姜规模,韩凤霞.西安城市地下水动态研究[J].城市勘测,2007(2):124-127.