2. 河北省地震动力学重点实验室, 河北三河 065201;
3. 中国地震局地球物理研究所, 北京 100081
2. Hebei Key Laboratory of Earthquake Dynamics, Sanhe 065201, Hebei, China;
3. Institute of Geophysics, China Earthquake Administration, Beijing 100081, China
地震可以引起井水温的变化,井水温的同震-震后响应能够揭示地壳介质对应力应变过程的响应。Shimamura等(1981)较早报道了对井水温度震前异常和同震响应的研究。在国内,付子忠(1988)首次用仪器观测到水温同震变化现象。前人对1999年中国台湾集集7.6级(Wang et al,2012)、2004年苏门答腊8.5级(孙小龙等,2008a)、2008年汶川8.0级(杨竹转等,2008;He et al,2020)、2011年日本9.0级(鱼金子等,2012;马玉川等,2014)、2013年芦山7.0级(张彬等,2013;巩浩波等,2015)、2015年尼泊尔8.1级(张彬等,2015)等地震的井水温同震-震后响应特征做了系统研究。陈大庆等(2007)对远场大震水温同震响应和机理进行研究,提出同震水温下降的气体脱逸模式。
在单井多震的研究方面,如石耀霖等(2007)研究唐山井水温的同震变化,使用有限单元法模型计算分析,认为井水垂直振荡时搅动井水引起的弥散效应是造成同震水温变化的主要原因。杨竹转等(2007)、孙小龙等(2008b)研究北京塔院井水温同震效应,认为井孔中的水体受振荡激发而加速对流与掺混是导致水温先下降的主要原因,水温同震幅度与震级、井震距等有关。杨竹转等(2018)研究丽江党校井水温同震变化特征,认为该井水温同震主要受水位泄流状态和水位变化的影响。在单井多层水温同震-震后响应的研究方面,如张慧等(2013)研究了海口ZK26井三层水温的同震和震后效应,认为观测层位的水文地质条件、背景噪声等的不同有可能造成不同层位井水温对远场大震的响应不同,远场大震井水温震后效应是由地震波导致的井-含水层系统参数的改变引起的。董蕾等(2017)研究发现随观测深度的增加,荣昌华江井三层水温同震响应能力增强。He等(2017)研究了川03井三层水温对2008年汶川8.0级和2011年日本9.0级地震的同震响应,认为水温同震变化与区域应力应变状态变化、含水层分布及特征密切相关,井孔中观测到的不同层的水温变化很可能受水平方向水流的控制。关于井水温微动态的形成机制,前人提出由水流动产生的热对流引起井水温变化的水热动力学机制和井区岩土中大地热流作用或热传导引起井水温变化的地热动力学机制,并指出多数井水温同震响应是水热动力学机制引起的,震后效应多是地热动力学机制引起的(车用太等,2008、2014),同时也开展了一些井水温动态数值模拟方面的研究(杨竹转,2011;顾申宜等,2013)。
重庆荣昌井水温对中、远场大震和近场地震均响应灵敏,积累了丰富的同震和震后观测资料。本研究统计分析了荣昌井水温对2008年1月—2021年9月全球MS≥7.0、川滇地区MS≥6.0、重庆及周边地区20km范围内MS≥4.0共273次地震的同震-震后响应动态特征,给出引起荣昌井水温同震变化的地震波能量密度阈值,对井水温同震优势方向成因及机理进行深入研究,为井水温动态的定量分析和机理解释提供研究思路。
1 构造背景与观测井概况荣昌井位于重庆市荣昌区安富镇(29.373°N,105.459°E),高程326m(图 1)。该井井深251m,其中0~5m深度段为套管,6~69m深度段设花管,70m深度以下为裸孔。井径240~110mm,井口套管直径240mm。井孔位于三叠系上统须家河组下亚组(T3 xjⅠ),岩性主要为中粒长石英砂岩、泥页岩及煤层。井孔地质概况为:0~5.2m为第四系松散岩土层;5.2~37.2m为须家河组灰黑色页岩、煤层,该层裂隙发育;37.2~61.1m为须家河组灰白色细砂岩,该层裂隙不发育;61.1~251.0m为须家河组中粒长石英砂岩,局部裂隙发育,岩层倾角45°,该层为主要观测含水层(图 2)。
该井构造上位于华蓥山基底断裂带西南侧,距断裂带最近距离约1km。华蓥山断裂位于中国西南地区,是一条右旋逆走滑型断裂,该断裂带北起万源,向南西经达川、荣昌至宜宾西南,全长约600km,为川中台拱与川东陷褶束2个三级构造单元的边界断裂。断裂走向N45°E,断面总体倾向SE,倾角30°~70°,在地表由规模不等的若干条断层组成。断裂带上的地表断裂有过多期活动,最新活动时代为距今(223.6±31.6)ka,历史中强震活动集中于断裂带的西南段(杨蓉等,2010)。
荣昌井水温自2008年1月正式观测以来,呈升温漂移型,2008年以来平均水温为23.70℃,年均升温约0.01℃,该井水温在2008年汶川8.0级地震时有明显的同震和震后效应(图 3)。荣昌井观测情况见表 1。
本文统计了2008年1月—2021年9月期间荣昌井水温对全球MS≥7.0、川滇地区MS≥6.0、重庆及周边地区MS≥4.0共273次地震的同震响应,统计结果见表 2。该井水温的同震变化往往伴随着水位的变化,水温的同震响应比(同震个数/地震个数)要远低于同井水位,荣昌井水温与水位同震响应比为1 ︰ 5.7,且荣昌井水温同震震级响应范围、震中距响应范围均小于同井水位,说明水温的同震响应能力弱于水位。
经统计发现,荣昌井水温26次同震响应方向始终为上升,井水温同震响应的地震分布见图 4,井水温典型同震形态见图 5。荣昌井3层不同深度水温(水温探头分别置于井下225m、215m、133m)在2014年3月6日—2016年3月13日同步观测,3层水温探头均位于砂岩主要观测含水层中(61.1~251.0m),且均处于温度正梯度范围内(图 2)。3层水温同步观测期间,探头225m处的同震响应次数最多(8次),其次为215m处(5次),最少为133m处(1次),说明荣昌井水温随观测深度增加,水温记录地震能力增强,在225m处的同震响应能力最强。其中,有2层以上水温同时响应的地震共有5次,分别为2015年4月25日尼泊尔8.1级(井震距2040km)、2015年5月30日日本小笠原8.0级(井震距3411km)、2015年11月14日中国东海海域7.2级(井震距2237km)、2015年12月7日塔吉克斯坦7.4级(井震距3146km)和2016年3月2日印尼苏门答腊岛海域7.8级地震(井震距2967km),且2层水温5次同震方向均为上升;3层水温同时响应的只有1次,即2015年尼泊尔8.1级地震。由图 6、表 3可见,荣昌井水位同震先于水温,水温同震是由深及浅的顺序发生,同震响应幅度随深度的增加而减小,水温的同震响应持续时间较水位长,且不同层位的水温同震方向均为上升。由此说明,单个井水温对不同地震的同震响应存在优势响应方向,其水温同震升降的性质均不因地震的远近、大小、震源机制或方位的变化而发生改变。
统计荣昌井水位和水温同时响应的9次地震同震响应参数(表 4),与同井水位同震变化方向进行比较,发现荣昌井水位-水温对中、远场震的同震多为同向上升正相关关系或振荡—上升,但在2008年汶川8.0级、2010年荣昌4.5级、2021年泸县6.0级地震时,为同震水位下降—水温上升的反相关关系,其可能原因为这3次地震为近场地震,其他地震均为中、远场地震,近场地震和中、远场地震引起的水位同震响应变化机制不一致所致(Lai et al,2021)。由此可见,井水温同震变化虽与水动力作用有关,但并非简单地由水位同震变化引起。
在此引入地震波能量密度e(r),表示地震波在传播过程中作用在单位体积地层介质上的最大地震能量值,其与地震震级M和震中距R存在一定的关系。荣昌井地震波能量密度分布与震级和震中距的关系如图 7所示,图中的直线为满足公式lgR=0.48M-0.33lge(r)-1.4 (Wang,2007)的能量密度等值线,坐标轴内的每个点代表一次地震。引起荣昌井水温同震变化的地震波能量密度e(r)>10-5J/m3,而引起荣昌井水温和水位反向同震变化的地震波能量密度e(r)>1J/m3。
荣昌井水温自观测以来表现出丰富的震后响应变化,一种为上升—下降型正脉冲变化,如2008年汶川8.0级和2010年荣昌4.5级地震;另一种为震后上升,持续数天后恢复到之前的变化水平,如2016年荣昌4.9级和2017年九寨沟7.0级地震。
荣昌井水温在2008年汶川8.0级地震后第5天开始出现上升—下降型正脉冲变化,2010年荣昌4.5级地震后亦如此,震后第11天开始出现上升—下降型正脉冲变化(图 8、表 5)。这2次地震后的荣昌井水温响应最为显著,且形态变化相似,汶川地震后的震后响应持续时间约为荣昌地震后的1倍,且仅有这2次地震及2021年泸县6.0级地震对应的荣昌水位同震响应方向为下降,其他地震对应的荣昌井水位同震均为上升或振荡。将荣昌井水温上升—下降脉冲型震后变化与荣昌4.5级地震震区ML≥1.0余震序列活动进行对比,发现随着荣昌井水温上升—下降脉冲型震后变化逐渐减小直至消失后,余震活动也由起伏逐渐趋于平静,表明荣昌井水温上升—下降脉冲型震后变化与震区的余震有较好的对应关系。
荣昌井水温在2010年玉树7.1级、2011年日本9.0级、2016年苏门答腊海域7.8级、2016年荣昌4.9级、2017年九寨沟7.0级地震后出现加速上升变化,持续数天后恢复到之前的变化水平。其中,2层水温(探头位置分别为井下225m和215m)在2016—2017年并行观测期间共记录到3次同步的上升—持续高值—下降型同震-震后响应变化,如2016年苏门答腊海域7.8级、2016年荣昌4.9级、2017年九寨沟7.0级等地震震后均有记录(表 6、图 9)。震级越大、震中距越小,荣昌井水温震后上升幅度越大,持续时间1~3个月。215m处水温的震后响应幅度和持续时间均大于225m处。
对于引发水温变化的动力源,地震释放的能量主要由摩擦能、破裂能和地震波能三部分组成。摩擦热主要集中在发震断裂附近,破裂能主要集中在震中区附近,地震波能占地震释放能量的比重不大,但可在大范围内释放(马玉川等,2014)。从动力源的角度看,荣昌井水温同震变化多数与地震波的影响有关。而地震波引起的渗透性变化的机制主要有:①堵塞或冲刷裂隙通道,造成颗粒移动(Brodsky et al,2003;Elkhoury et al,2011;Candela et al,2014);②新的微裂隙生成,含水层系统被破坏(Manga et al,2012);③水滴、气泡等的迁移(Manga et al,2012)。地震发生后,由于裂隙里的沉积物被冲刷而新形成的通道会被重新堵塞,使得渗透性能够逐渐恢复到震前水平。渗透率恢复到预估值的时间取决于重新堵塞孔隙的时间和封闭孔隙的地球化学过程所需的时间。井水温的震后响应是由于地震波周期性振荡作用于含水层系统后,一方面激活了孔隙、裂隙中的充填物(如气体、滞水等),另一方面促使局部裂隙的串通,使得各水体间及水体与围岩的热量交换加强或减弱,从而导致井水温变化。
目前,井水温同震上升机理主要有:①在地震波震动作用下,承压含水层下部热水上升引起温度上升;②地震波引起的水体波动与井壁及水温传感器摩擦引起温度上升;③井孔水温梯度为负时,上部热水在震动作用下下渗引起温度上升(毛慧玲等,2012)。井水温同震下降机理主要有:①气体逸出说,即当井水气体释放时,同时释放出井水中的热量,从而降低了井水温度;②热弥散说,井水垂直振荡时搅动井水,引起的弥散效应造成同震水温变化,在一定的条件下形成同震水温降低现象;③冷水下渗说,即井孔含水层周边上层地下水随着振动效应的作用,向下垂直运动的速率有所加大,低温水快速混入观测含水层中,引起温度的快速下降(孙小龙,2008c)。浅循环含水层表明水温将会下降,而深循环含水层水温将会升高。通过浅循环含水层向井孔注入冷水,通过深循环含水层向井孔注入热水。
由地貌、地质特征可以看出,荣昌井区附近的地下水系主要由大气降水补给,在含水裂隙带被切割处或侵蚀基准面附近进行排泄(图 10(a))。该井含水层地层出露在背斜区,地层由平缓向两翼变陡,造成自流斜坡,地下水富集具有一定承压性。砂岩厚度大,受力后产生的张裂隙较多,是主要含水层,下部页岩相对隔水(图 10(b))。同位素测试结果表明,荣昌井水为大气降水成因,piper三线图显示该井水水化学类型为Na-Cl型,Na-K-Mg三角图可以较好地用于解释水的起源(Giggenbach,1988),结果显示荣昌井水为部分平衡水,说明水的循环深度较深,携带一定的深部物质(图 11)。荣昌井区从20世纪50年代后集中开采50多年,采煤层延伸至观测井底部。2016年以来,荣昌井区周边陆续开始进行页岩气开采。页岩气是赋存于有机质泥页岩及其夹层中,以吸附和游离状态为主要存在方式的非常规天然气,成分以甲烷为主。截至2021年10月,荣昌井井面时有气泡上涌,佐证了该井含水层含气的事实。2021年4月在井口架设了实时摄像头,记录到5月22日青海玛多7.4级地震震时气泡量明显增多(图 12),而震时荣昌井水温也出现了上升变化(图 5)。
注:T3xjⅠ:三叠系上统须家河组下亚组,主要为泥页岩、砂岩及煤层;T3xjⅡ:三叠系上统须家河上亚组,主要为砂岩、泥页岩;J1z:侏罗系下统珍珠冲组,主要为泥页岩夹砂岩;Q4:第四系全新统,主要为砾石、砂、黏土。 |
Craig(1961));LMLW为重庆地区大气降水线(据胡菡等(2015))。 | 注:δD-δ18O分布图中GMLW为全球大气降水线(据
气体具有很大的活动性,在地壳内分布广且分散,一般分布在岩石的孔隙、裂隙及溶洞中,有的溶解于水中,有的封存在合适的构造中形成气芷。地震、构造活动、应力的变化、局部气体的形成、承压水的变化和流动均可形成压力梯度,造成气体从压力高的地方向压力低的地方流动或迁移(方震等,2012)。当地震波经过井区含水层时,会造成裂隙串通、井内水体连续振荡,深部含水层孔隙中的气体及溶解于井水中气体就会随地下水一起运移至表面。荣昌井水温同震上升变化可能是地震波的扰动造成井下深部气体大于本底的释放量,大量气泡沿裂隙上升进入井含水层系统而引起的。震后一段时间,气体释放量逐渐减小,水温也逐渐恢复至震前水平。
荣昌井水温的同震-震后响应变化较为特殊,其变化机理也不同于前人所提出的井水温同震上升机理。而根据对西南地区井水温同震响应特征的分析(巩浩波等,2021)可知,并非所有井水温同震优势下降的井位都是浅循环井,如西昌川32井,井深410.2m,观测层128.17~410.20m。水温同震优势方向上升型和下降型井位并无明显共性。水温的同震变化是探头周围的地下水与其周围介质进行能量交换的结果,与水位变化幅度及运动方式、水流速度、井孔周围孔隙度、井孔温度梯度、水温探头放置位置等有关。由此说明井水温同震-震后响应机理的复杂性,今后仍需对更多的井孔进行更加深入的研究。
5 结论重庆荣昌井水温对远大震和近震均响应灵敏,积累了丰富的同震和震后观测资料,且荣昌井(多层)水温同震响应方向均为上升,是研究井水温同震-震后响应动态特征的优质井位。本文统计分析了荣昌井水温对2008年1月—2021年9月期间共273次地震的同震-震后响应动态特征,给出引起该井水温同震变化的地震波能量密度阈值,对井水温同震优势方向成因和机理进行了深入研究,得到以下结论:
(1) 荣昌井水温同震-震后响应能力较好,对近震和远震均可记录到。该井水温同震响应由深及浅顺序发生,同震-震后响应持续时间随观测深度的增加而增加,同震响应幅度随观测深度的增加而减小,且该井水温同震-震后响应持续的时间较水位更长。
(2) 荣昌井多层水温同震响应方向始终为上升,说明单个井水温对不同地震的同震响应存在优势响应方向,其水温同震升降的性质均不因地震的远近、大小、震源机制或方位的变化而改变。荣昌井水温同震响应优势方向上升变化可能是地震波的扰动造成井下深部气体释放,沿裂隙上升进入井含水层观测层而引起的。
(3) 荣昌井水位-水温对中、远场地震的同震多为同向上升正相关关系或振荡—上升,但在汶川8.0级、荣昌4.5级、泸县6.0级地震时为同震水位下降—水温上升的反相关关系,其原因可能为这3次地震为近场地震,其他地震均为中、远场地震,近场地震和中、远场地震引起的水位同震响应变化机制不一致。引起荣昌井水温同震变化的地震波能量密度e(r)>10-5J/m3,而引起荣昌井水温和水位同震反向变化的地震波能量密度e(r)>1J/m3。
车用太、何案华、鱼金子, 2014, 水温微动态形成的水热动力学与地热动力学机制, 地震学报, 36(1): 106-117. DOI:10.3969/j.issn.0253-3782.2014.01.009 |
车用太、刘成龙、鱼金子, 2008, 井水温度微动态及其形成机制, 地震, 28(4): 20-28. DOI:10.3969/j.issn.1000-3274.2008.04.003 |
陈大庆、刘耀炜、杨选辉等, 2007, 远场大震的水位、水温同震响应及其机理研究, 地震地质, 29(1): 122-132. DOI:10.3969/j.issn.0253-4967.2007.01.011 |
董蕾、陈敏、马伟等, 2017, 荣昌华江井水温多层微动态特征, 地震地磁观测与研究, 38(6): 72-78. DOI:10.3969/j.issn.1003-3246.2017.06.013 |
方震、刘耀炜、杨选辉等, 2012, 地震断裂带中气体来源及运移机制研究进展, 地球物理学进展, 27(2): 483-495. DOI:10.6038/j.issn.1004-2903.2012.02.011 |
付子忠, 1988. 地热动态观测与地热前兆. 见: 地壳构造与地壳应力文集(1). 北京: 地震出版社, 1~8.
|
巩浩波、郭卫英、李光科等, 2015, 重庆井网水位水温同震响应特征分析, 震灾防御技术, 10(增刊): 794-804. |
巩浩波、李光科、李翠平等, 2021, 西南地区地震观测井泉水温动态特征研究, 地球物理学进展, 36(4): 1393-1403. |
顾申宜、刘博维、张慧等, 2013, 井孔水温对远场巨震同震响应及其机制的数模研究, 地震, 33(1): 29-39. |
胡菡、王建力, 2015, 重庆市2013年10—12月大气降水中氢氧同位素特征及水汽来源分析, 中国岩溶, 34(3): 247-253. |
马玉川、刘耀炜、马未宇, 2014, 2011年日本M w9.0地震引起的环渤海地区井水温度变化分析, 地震, 34(3): 40-49. DOI:10.3969/j.issn.1000-3274.2014.03.004 |
毛慧玲、秦嘉政、杨正纲, 2012, 日本9.0级和汶川8.0级地震引起云南地区井水温同震响应研究, 高原地震, 24(1): 16-22. DOI:10.3969/j.issn.1005-586X.2012.01.003 |
石耀霖、曹建玲、马丽等, 2007, 唐山井水温的同震变化及其物理解释, 地震学报, 29(3): 265-273. DOI:10.3321/j.issn:0253-3782.2007.03.005 |
孙小龙、刘耀炜, 2008a, 苏门答腊8.5级地震引起的水温响应变化, 华北地震科学, 26(1): 35-40. |
孙小龙、刘耀炜, 2008b, 塔院井水位和水温的同震响应特征及其机理探讨, 中国地震, 24(2): 105-115. |
孙小龙、刘耀炜、王博, 2008c, 宾川井对印尼大震的同震响应特征及其机理解释, 地震, 28(3): 69-78. |
杨蓉、尊珠桑姆、许长海等, 2010, 四川盆地东部华蓥山断裂滑动分析与古应力重建, 内蒙古石油化工, 36(4): 97-100. |
杨竹转, 2011. 地震波引起的井水位水温同震变化及其机理研究. 博士学位论文. 北京: 中国地震局地质研究所.
|
杨竹转、邓志辉、邓盛昌等, 2018, 井水温度动态变化分析——以云南丽江党校井为例, 地震学报, 40(2): 172-184. |
杨竹转、邓志辉、刘春国等, 2008, 中国大陆井水位与水温动态对汶川M s8.0地震的同震响应特征分析, 地震地质, 30(4): 895-905. |
杨竹转、邓志辉、陶京玲等, 2007, 北京塔院井数字化观测水温的同震效应研究, 地震学报, 29(2): 203-213. |
鱼金子、车用太、刘成龙等, 2012, 金沙江水网对日本9.0级地震的同震响应及其特征与机理, 地震, 32(1): 59-69. |
张彬、刘耀炜、高小其等, 2015, 2015年尼泊尔M s8.1地震引起的井水位与井水温同震效应及其相关性分析, 地震学报, 37(4): 533-540. |
张彬、刘耀炜、杨选辉, 2013, 中国大陆井水温对汶川8.0级、玉树7.1级、芦山7.0级和岷县6.6级地震响应特征的对比研究, 地震工程学报, 35(3): 535-541. |
张慧、顾申宜、李志雄等, 2013, 单井多层位水温微动态初步研究, 地震, 33(1): 101-110. |
Brodsky E E, Roeloffs E, Woodcock D, et al, 2003, A mechanism for sustained groundwater pressure changes induced by distant earthquakes, J Geophys Res, 108(B8): 2390. |
Candela T, Brodsky E E, Marone C, et al, 2014, Laboratory evidence for particle mobilization as a mechanism for permeability enhancement via dynamic stressing, Earth Planet Sci Lett, 392: 279-291. |
Craig H, 1961, Isotopic variations in meteoric waters, Science, 133(3465): 1702-1703. |
Elkhoury J E, Niemeijer A, Brodsky E E, et al, 2011, Laboratory observations of permeability enhancement by fluid pressure oscillation of in situ fractured rock, J Geophys Res: Solid Earth, 116(B2): B02311. |
Giggenbach W F, 1988, Geothermal solute equilibria. Derivation of Na-K-Ma-Ca geoindicators, Geochim Cosmochim Acta, 52(12): 2749-2765. |
He A H, Singh R P, 2020, Coseismic groundwater temperature response associated with the Wenchuan earthquake, Pure Appl Geophys, 177(1): 109-120. |
He A H, Zhao G, Sun Z H, et al, 2017, Co-seismic multilayer water temperature and water level changes associated with Wenchuan and Tohoku-Oki earthquakes in the Chuan no. 03 well, China, J Seismol, 21(4): 719-734. |
Lai G J, Lei X L, Jiang C S, et al, 2021, Multiple mechanisms of coseismic water level changes at the Rongchang well in a seismically active area in China, Tectonophysics, 819: 229083. |
Manga M, Beresnev I, Brodsky E E, et al, 2012, Changes in permeability caused by transient stresses: field observations, experiments, and mechanisms, Rev Geophys, 50(2): RG2004. |
Shimamura H, Watanabe H, 1981, Coseismic changes in groundwater temperature of the Usu volcanic region, Nature, 291(5811): 137-138. |
Wang C Y, 2007, Liquefaction beyond the Near Field, Seismol Res Lett, 78(5): 512-517. |
Wang C Y, Manga M, Wang C H, et al, 2012, Transient change in groundwater temperature after earthquakes, Geology, 40(2): 119-122. |