顏玉聰 劉峰立 郭麗爽 周曉成 歐陽(yáng)澍培 李靜超 王萬(wàn)麗

摘要：對(duì)地震活動(dòng)強(qiáng)烈地區(qū)進(jìn)行水文地球化學(xué)監(jiān)測(cè)可以很好地獲取與地震相關(guān)的深部流體地球化學(xué)信息。通過(guò)研究龍門山斷裂帶及其周圍13個(gè)溫泉及觀測(cè)井的水文地球化學(xué)特征，建立該斷裂帶溫泉水文循環(huán)模型，揭示其水化學(xué)變化與地震活動(dòng)的關(guān)系。結(jié)果表明：①龍門山斷裂帶溫泉水主要來(lái)自其周圍0.8～3.2 km高山的大氣降水;②溫泉水化學(xué)類型沿龍門山斷裂從北到南、從西到東依次為重碳酸型、硫酸型、氯化物型，溫泉水的循環(huán)深度、水巖反應(yīng)程度及微量元素富集因子不斷增加;③在汶川MS8.0地震、蘆山MS7.0地震發(fā)生后，震中距200 km以內(nèi)的溫泉的離子組分呈下降趨勢(shì)，這可能是在震后愈合過(guò)程中，龍門山斷裂內(nèi)部滲透性減弱、水巖反應(yīng)程度衰減所致。

關(guān)鍵詞：溫泉;水文地球化學(xué);地震;龍門山斷裂帶

中圖分類號(hào)：P314.1?? 文獻(xiàn)標(biāo)識(shí)碼：A?? 文章編號(hào)：1000-0666（2021）02-0170-15

0 引言

溫泉多沿?cái)嗔褞С雎?，是溝通地殼深部與淺部的窗口，和深部地下水系統(tǒng)有著密切的關(guān)系（Zhang et al，2015;Zhou et al，2017）。研究活動(dòng)斷層的溫泉流體地球化學(xué)特征有助于深入了解該斷層的活動(dòng)方式、活動(dòng)強(qiáng)度以及切割深度等斷層運(yùn)動(dòng)特性，對(duì)探索斷層帶深部、特別是震源區(qū)附近介質(zhì)的物理化學(xué)環(huán)境、研究各類地球化學(xué)組分的地震前兆機(jī)制、尋找特定地區(qū)的前兆觀測(cè)組分有重要意義（Lambrakis，Kallergis，2005;Liu et al，2009）。因此，在地震易發(fā)區(qū)對(duì)溫泉的水文地球化學(xué)特征進(jìn)行監(jiān)測(cè)，可捕捉來(lái)自地殼深部活動(dòng)的信息（蒲小武等，2013;王云等，2014;Walraevens et al，2018;Xu et al，2019）。

近年來(lái)，國(guó)內(nèi)外許多專家對(duì)溫泉?jiǎng)討B(tài)與地震活動(dòng)的關(guān)系進(jìn)行了研究（Tsunogai，Wakita，1995;Shi et al，2020）。Skelton等（2014）對(duì)冰島2012年10月以及2013年4月2次5級(jí)地震之間以及地震前后4～6個(gè)月的地下水中穩(wěn)定同位素比率和鈉、鈣等溶質(zhì)的濃度變化規(guī)律進(jìn)行研究，獲得了該地區(qū)地震前后地下水水文地球化學(xué)變化的規(guī)律;Okuyama等（2016）研究發(fā)現(xiàn)1965—1968年日本松代地震群前后水中某些溶解組分隨井深呈線性增長(zhǎng);Pérez 等（2008）研究1995年西班牙加利西亞4.6級(jí)地震以及1997年5.3級(jí)地震發(fā)現(xiàn)，在震中距為90 km的地下水觀測(cè)井中Cl-濃度明顯增加。另外，Shi等（2020）的研究顯示，氫氧同位素以及微量元素的變化對(duì)地震的響應(yīng)也較為明顯，2018年通海MS5.0地震后，江川觀測(cè)井的水化學(xué)以及氫氧同位素明顯上升，微量元素（Li+、Pb2+、As2+等元素）的濃度呈明顯下降。

龍門山斷裂帶是青藏高原東緣山脈與揚(yáng)子板塊的結(jié)合部位，地質(zhì)構(gòu)造復(fù)雜，地震活動(dòng)劇烈（Lin et al，2018）。自2008年汶川MS8.0地震后，我國(guó)在該地區(qū)建立了大量的流動(dòng)觀測(cè)網(wǎng)、溫泉和熱水井固定觀測(cè)點(diǎn)，豐富了該地區(qū)的流體地球化學(xué)背景（程萬(wàn)正等，2013）。目前，關(guān)于龍門山斷裂帶地下流體的研究主要集中在地震前后溫泉水溫、水化學(xué)組分、水氡以及溫泉逸出氣異常等方面（Chen et al，2015;晏銳等，2015;Zhou et al，2015;楊耀等，2019）。Chen等（2015）對(duì)2008—2010年川西地區(qū)32個(gè)溫泉的水化學(xué)變化進(jìn)行了研究，發(fā)現(xiàn)了汶川地震后溫泉的離子濃度呈下降變化;Zhou等（2015）對(duì)該地區(qū)32個(gè)溫泉?dú)膺M(jìn)行了3次采樣，結(jié)果表明汶川MS8.0地震后，四川西部地區(qū)有更多地幔流體進(jìn)入地殼，地幔成分的流體上涌可能與大地震的發(fā)生有關(guān)系;楊耀等（2019）研究了川西地區(qū)地下流體觀測(cè)井的地下水的水化學(xué)特征，確定了四川西部地區(qū)主要觀測(cè)井的水文地球化學(xué)背景。

前人對(duì)龍門山斷裂帶溫泉水文地球化學(xué)特征的研究主要集中于川西地區(qū)附近，沒(méi)有對(duì)整個(gè)龍門山斷裂帶展開(kāi)更為細(xì)致的分析，鑒于此，本文對(duì)龍門山斷裂帶13個(gè)溫泉及觀測(cè)井的水文地球化學(xué)特征進(jìn)行分析，針對(duì)該斷裂帶不同斷裂計(jì)算溫泉水的熱儲(chǔ)溫度與循環(huán)深度，確定溫泉水的補(bǔ)給來(lái)源和補(bǔ)給高程，判斷其水文地球化學(xué)的控制因素和影響因素，建立龍門山斷裂帶溫泉水循環(huán)模型。通過(guò)探討2008—2016年龍門山斷裂帶溫泉水文地球化學(xué)變化與汶川、蘆山地震活動(dòng)的關(guān)系，為該地區(qū)震情跟蹤、異常落實(shí)和地震預(yù)測(cè)提供基礎(chǔ)研究。

1 研究區(qū)地質(zhì)概況

龍門山斷裂帶是青藏高原東緣的邊界山脈，北起廣元、南至天全，長(zhǎng)約500 km、寬約30～50 km，呈NE—SW向展布，形成于印支晚期，受多期復(fù)雜構(gòu)造運(yùn)動(dòng)影響（鄧起東等，1994;Clark，Royden，2000;Godard et al，2009;Kouketsu et al，2017）。該構(gòu)造帶是由一系列大致平行的逆沖、逆沖兼走滑的斷裂構(gòu)成，自西向東發(fā)育有茂汶—汶川斷裂（F1）、北川—映秀斷裂（F2）、安縣—灌縣斷裂（F3）和廣元—大邑隱伏斷裂帶（F4）以及一系列龍門山山前斷裂（圖1）（王二七等，2001;徐錫偉等，2008;付碧宏等，2008）。

龍門山地區(qū)位于高原氣候區(qū)與季風(fēng)氣候區(qū)交界處，受大氣降水和高山冰雪融水的共同補(bǔ)給，水系較為發(fā)育（李?yuàn)^生等，2015），從西南到東北依次發(fā)育青衣江、岷山、沱江、涪江和嘉陵江等一系列斜交于龍門山斷裂帶的河流（王志等，2010）。龍門山斷裂帶由多期構(gòu)造變動(dòng)形成，水文地質(zhì)條件復(fù)雜，溫泉多沿?cái)嗔褞Оl(fā)育且呈帶狀分布。研究區(qū)主要的含水層包括第四系松散堆積層孔隙含水層、侏羅系—白堊系碎屑巖裂隙孔隙—基巖裂隙含水層、震旦系—三疊系碳酸鹽巖巖溶裂隙—基巖裂隙水以及前震旦系基底含水層4大含水層，且在褶皺和斷裂構(gòu)造共同作用下，各含水層之間相互聯(lián)系（陶廣斌，2019）。前人已收集到的地震數(shù)據(jù)顯示龍門山地區(qū)最早記錄的地震為1161年北川MS4.45地震，根據(jù)中國(guó)地震局監(jiān)測(cè)預(yù)報(bào)司統(tǒng)計(jì)，自公元1161年至2020年7月，本區(qū)共發(fā)生MS>43/4地震196次，其中MS5.0～5.9地震112次，MS6.0～6.9地震21次，MS>7.0地震5次，因此，該區(qū)具有重要的研究?jī)r(jià)值（張?jiān)罉虻龋?008）。

2 研究方法

筆者對(duì)2008—2016年龍門山斷裂帶13個(gè)溫泉（表1）進(jìn)行了9次調(diào)查，共59個(gè)溫泉水樣，每個(gè)溫泉水樣采集4瓶用于測(cè)定其主量元素、微量元素、氫氧同位素及SiO2含量。野外采樣時(shí)，采樣容器為250 mL的聚乙烯塑料瓶，將聚乙烯瓶置于水中潤(rùn)洗2～3次后裝樣，取樣時(shí)將瓶?jī)?nèi)空氣排凈，避免水樣受空氣影響，并現(xiàn)場(chǎng)測(cè)定溫泉水的水溫、pH、電導(dǎo)率，并向測(cè)定微量元素的樣品中加入1.00 mL的濃硝酸（pH<2），在室溫條件下保存。

溫泉水樣品中常量元素的測(cè)定在中國(guó)地震局地震預(yù)測(cè)重點(diǎn)實(shí)驗(yàn)室完成，使用儀器為Dionex ICS-900離子色譜儀及AS40自動(dòng)進(jìn)樣器，檢測(cè)限為0.01 mg/L。溫泉水中的HCO-3和CO2-3濃度使用ZDJ-3D型電位滴定儀通過(guò)標(biāo)準(zhǔn)滴定程序進(jìn)行測(cè)量，滴定所用鹽酸濃度為0.05 mol/L。為了校準(zhǔn)色譜，在測(cè)量每批水樣品之前需測(cè)量標(biāo)準(zhǔn)樣品，測(cè)量偏差在±2‰以內(nèi)。離子平衡（ib）計(jì)算公式如下：

ib=∑陽(yáng)離子-∑陰離子（∑陽(yáng)離子+∑陰離子）×0.5×100%（1）

式中：陰陽(yáng)離子的測(cè)量誤差均小于5%。

微量元素的含量分析在中國(guó)地震局地殼應(yīng)力研究所采用Element型電感耦合等離子體質(zhì)譜儀ICP-MS完成，樣品的分析精度優(yōu)于5%。氫氧同位素分析采用氣體同位素質(zhì)譜儀MAT253，以國(guó)際標(biāo)準(zhǔn)即標(biāo)準(zhǔn)平均海洋水（SMOW）為標(biāo)準(zhǔn)，分別記為δD和δ18O，氫氧同位素參考標(biāo)準(zhǔn)物質(zhì)為國(guó)家標(biāo)準(zhǔn)物質(zhì)（GBW04401，04402），樣品的分析精度為δD<±1‰和δ18O<±0.2‰。利用電感耦合等離子體發(fā)射光譜儀Optima-5300DV檢測(cè)SiO2含量。

3 溫泉水樣分析結(jié)果

對(duì)所采集的溫泉水樣進(jìn)行分析，結(jié)果見(jiàn)表2。研究區(qū)溫泉水樣水溫變化范圍為10.8℃～64.0 ℃;pH變化范圍為6.86～9.50，平均值為7.96，均為典型的弱堿性水;電導(dǎo)率變化范圍為154.00～132 400.00 μs/cm，TDS的變化范圍為0.08～201.37 g/L，其中WJJ為鹵水，LFS、ZGS為鹽水，HSW、TC、PL為咸水，GYG、SZZ、JYG為淡水，LX、WC、ZBS、ZLJ為微咸水溫泉。從表2還可看出，水樣中的主量元素相對(duì)豐度：陽(yáng)離子主要為Na+、Mg2+、Ca2+;陰離子主要為SO2-4、HCO-3、Cl-;水樣中δD和δ18O測(cè)值的分布范圍分別為-6.18‰～-17.80‰和-51.50‰～-119.20‰。微量元素含量較低，其豐度結(jié)果見(jiàn)表3，SiO2含量的分布范圍為0.81～29.70 mg/L。

4 討論

4.1 溫泉水來(lái)源

Craig（1961）通過(guò)總結(jié)全球各地區(qū)大氣降水中的δD、δ18O數(shù)值變化范圍及二者之間的關(guān)系，最先提出了用于判斷地下水現(xiàn)代補(bǔ)給來(lái)源的全球大氣降水線（GMLW）。在此之后，不斷有學(xué)者提出了一系列降水線方程，本文選取了Xu等（2015）對(duì)龍門山地區(qū)建立的大氣降水線方程（LMLW）：δD=7.6δ18O+11.3，繪制龍門山斷裂帶溫泉δD-δ18O關(guān)系圖。從圖2可觀察到，龍門山斷裂帶的溫泉水同位素?cái)?shù)據(jù)主要分布在大氣降水線兩側(cè)，表明研究區(qū)溫泉水及觀測(cè)井水主要來(lái)源于大氣降水。

研究區(qū)水樣δD和δ18O測(cè)值的分布范圍分別為-6.18‰～-17.80‰和-51.50‰～-119.20‰，由于龍門山地區(qū)海拔變化較大，氫氧同位素會(huì)受到高程效應(yīng)的影響，因此可以利用同位素的高程效應(yīng)估算溫泉水的補(bǔ)給高程（Clark，Royden，2000;于津生，1997;Zhou et al，2010;劉成龍等，2020）。

我國(guó)西部降水的δD值與高程H的關(guān)系為（于津生，1997）：

δD=-0.026H-30.2（2）

中國(guó)大氣降水δD的高程效應(yīng)可以表示為（Zhou et al，2010）：

δD=-0.03HALT -27（3）

式中：HALT為海拔高程值。

根據(jù)這些公式計(jì)算研究區(qū)溫泉補(bǔ)給高程再取平均值（表4）。由表4可見(jiàn)，研究區(qū)溫泉的補(bǔ)給高程大致為0.8～3.2 km;斷裂帶南西段溫泉的補(bǔ)給高程為1.4～3.2 km，補(bǔ)給高程較高;北東段溫泉的補(bǔ)給高程為0.8～1.9 km，相對(duì)較低，這與龍門山斷裂帶總體呈現(xiàn)由南向北、由西向東傾斜的趨勢(shì)一致。

4.2 水化學(xué)成分特征

4.2.1 常量元素

本文使用了舒卡列夫分類法（毫克當(dāng)量百分比分別大于25%的陰陽(yáng)離子參與命名）將龍門山斷裂帶的溫泉水劃分為Na-HCO3、Ca·Na-HCO3、Na-HCO3·SO4、Ca-HCO3·SO4、Ca·Mg-SO4、Na-HCO3·Cl、Na-SO4·Cl、Na-Cl八種水化學(xué)類型，利用origin軟件繪制Piper三線圖，觀測(cè)離子分配規(guī)律（圖3），從空間上看，該地區(qū)水化學(xué)類型復(fù)雜，從時(shí)間上看，該地區(qū)溫泉水化學(xué)類型在采樣期間沒(méi)有明顯改變。

位于茂縣—汶川斷裂帶的GRG、JYG溫泉水化學(xué)類型為Na-HCO3、Na-HCO3·SO4，礦化度較低，其含水層為分別為侏羅系蓮花口組鈣質(zhì)砂巖、粉砂巖以及震旦系燈影組白云巖、頁(yè)巖。溫泉中Na+主要來(lái)源于鈉長(zhǎng)石、鉀長(zhǎng)石、鈣長(zhǎng)石等礦物發(fā)生的溶解作用和陽(yáng)離子吸附的交替作用。HCO-3主要來(lái)源于灰?guī)r、白云巖等碳酸鹽巖的溶解，一般來(lái)說(shuō)，隨著地下水入滲深度的增加，溶濾作用增強(qiáng)，溫泉水中的SO2-4含量也會(huì)增加。WC、LX溫泉的水化學(xué)類型為Ca·Mg-SO4，出露于志留系-泥盆系泥質(zhì)灰?guī)r夾黃鐵礦巖層中，含水介質(zhì)中黃鐵礦的氧化作用可加快灰?guī)r及白云巖的溶解，從而使WC、LX溫泉富集Ca2+、Mg2+及SO2-4。

位于龍門山地區(qū)北川—映秀斷裂、安縣—灌縣斷裂及廣元—大邑隱伏斷裂帶南段附近的WJJ、ZGS、HSW溫泉，靠近四川盆地，礦化度高，為鹵水和咸水，水化學(xué)類型為Na-Cl型及Ma-SO2-4型，而LFS、ZBS、ZLJ溫泉水化學(xué)類型為Na-HCO3·Cl型。這可能是由于四川盆地從震旦系到白堊系各個(gè)地層均賦存有地下鹵水（周訓(xùn)等，2015），而在構(gòu)造應(yīng)力相對(duì)集中帶、斷裂破碎帶和裂隙發(fā)育帶是鹵水的主要富集地帶（林耀庭，2001;楊耀等，2019），因此WJJ、ZGS、HSW溫泉的形成可能起源于古海水，在水體循環(huán)過(guò)程中以濃縮作用為主，最終形成高礦化度的深部鹵水（曹琴等，2015）。而LFS、ZBS、ZLJ溫泉可能是地下水體循環(huán)過(guò)程中經(jīng)過(guò)含鹽地層，水體溶濾大量氯化物巖類（周訓(xùn)，1993），從而導(dǎo)致地下水的礦化度增高。研究區(qū)中段SZZ溫泉水化學(xué)類型為Ca·Na-HCO3型，而北段PL、TC溫泉的水化學(xué)類型為Ca-SO4·HCO3，都出露于三疊系嘉陵江組灰?guī)r、白云巖。而北段溫泉中SO2-4含量增加是由該溫泉含水介質(zhì)中三疊系嘉陵江組巖層中夾厚層硬石膏所致。

4.2.2 微量元素

深循環(huán)地下熱水中微量元素的含量特征與區(qū)域地球化學(xué)背景值、元素性質(zhì)、地下水的物理化學(xué)條件、水-巖相互作用等密切相關(guān)（Zhang et al，2016）。研究區(qū)溫泉圍巖以白云巖、泥質(zhì)灰?guī)r、碎屑巖及砂巖為主，溫泉水循環(huán)于碳酸鹽巖中，微量元素富集因子的大小可以定性判斷溫泉的水巖反應(yīng)程度（Soto-Jiménez，Páez-Osuna，2001），其計(jì)算方法為：

EFi=（Ci/CR）w /（Ci/CR）r（4）

式中：CR為選定的參比元素含量;Ci為樣品中元素含量;w為水樣中元素濃度;r為巖石中元素濃度。因此將龍門山地區(qū)方解石及石英巖脈的微量元素含量（薛鈞月，2009）作為標(biāo)準(zhǔn)，使用Ti濃度對(duì)微量元素?cái)?shù)據(jù)進(jìn)行歸一化，計(jì)算研究區(qū)微量元素富集因子EF（圖4）。從圖4中可以明顯看到微量元素的富集狀況，結(jié)合微量元素含量，對(duì)研究區(qū)微量元素地球化學(xué)特征進(jìn)行分析，發(fā)現(xiàn)斷裂帶南西段溫泉微量元素相對(duì)富集，且微量元素富集因子比北東段大。

4.3 水巖反應(yīng)程度

Na-K-Mg三角圖能夠反應(yīng)溫泉的水巖平衡狀態(tài)和平衡溫度（Giggenbach，1988）。從研究區(qū)的水巖反應(yīng)圖（圖5）可以看出，WJJ溫泉落在完全平衡區(qū)。林耀庭等（1997）和Zhou等（2015）根據(jù)四川盆地地下鹵水的δD和δ18O值將其劃分為4種類型：大氣水淋濾型鹵水、海相沉積型鹵水、海相沉積-大氣降水疊加型鹵水以及海相沉積-巖漿水疊加型鹵水。WJJ溫泉礦化度較高，其δD和δ18O落于大氣降水線上部，表明其成因與大氣降水有關(guān)，為大氣淋濾型鹵水，且水-巖反應(yīng)達(dá)到離子平衡狀態(tài)。HSW06、ZLJ06、LFA06、ZBS06、ZGS06、ZGS07為部分平衡水，其它溫泉為未成熟水在圖中Mg1/2的一側(cè)，表明水巖反應(yīng)程度低，并受到一些淺層冷水混合作用。

4.4 熱儲(chǔ)溫度與循環(huán)深度

目前常用的地?zé)釡貥?biāo)主要分為兩大類，一類是陽(yáng)離子溫標(biāo)，另一類為SiO2溫標(biāo)。陽(yáng)離子溫標(biāo)方法一般用于估算部分平衡水及平衡水的熱儲(chǔ)溫度（Fournier，Rowe，1966），陽(yáng)離子溫標(biāo)法中K-Mg溫標(biāo)適用于中低溫地?zé)崴?。研究區(qū)WJJ平衡水溫泉及HSW、ZLJ、LFA、ZBS、ZGS溫泉等部分平衡水為中低溫地?zé)崴到y(tǒng)，因此采用K-Mg溫標(biāo)計(jì)算熱儲(chǔ)溫度更適合（王逸凌等，2020）。而SiO2相對(duì)于其他礦物而言，穩(wěn)定性較高，能夠指示未成熟水的熱儲(chǔ)溫度（龐忠和等，1990），因此該地區(qū)未成熟水溫泉采用SiO2溫標(biāo)。通過(guò)計(jì)算得到龍門山地區(qū)熱儲(chǔ)溫度范圍為20.1 ℃～123.4 ℃（表5）。

研究區(qū)溫泉及觀測(cè)井為大氣降水入滲補(bǔ)給通過(guò)深循環(huán)獲得深部熱流加熱之后在升流至地表形成，其熱水循環(huán)深度計(jì)算公式如下：

H=（T-T0）/g+h（5）

式中：H為循環(huán)深度（km）;T為地?zé)崴臒醿?chǔ)溫度（℃）;T0為研究區(qū)的平均氣溫（℃），研究區(qū)內(nèi)平均氣溫取17 ℃;g為地溫梯度，為24.5 ℃/km;h為常溫帶深度，為30 m（崔希林，2014）。

由表5得到龍門山斷裂帶內(nèi)熱水循環(huán)深度范圍為0.1～4.4 km（表5）。其北東段TC及PL溫泉循環(huán)深度為1.1～1.5 km，中部GRG等溫泉的循環(huán)深度為0.7～2.5 km，其中以Cl-為主的LFS、ZBS、ZLJ溫泉，其循環(huán)深度為0.9～2.8 km，北部HSW、WJJ、ZGS溫泉循環(huán)深度為4.0～4.4 km。結(jié)合溫泉常量元素及微量元素特征可以發(fā)現(xiàn)，以HCO-3為主的溫泉，循環(huán)深度較淺，SO2-4次之，而以Cl-為主要陰離子的溫泉循環(huán)深度較深。結(jié)合Na-K-Mg三角圖（圖4）、水化學(xué)組分特征以及循環(huán)深度綜合分析，龍門山斷裂帶溫泉的水化學(xué)類型沿茂縣—汶川斷裂帶、北川—映秀斷裂、安縣—灌縣斷裂及廣元—大邑隱伏斷裂帶從北到南、從西到東依次為重碳酸型、硫酸型、氯化物型溫泉水，其水巖反應(yīng)程度、水化學(xué)組分含量、循環(huán)深度沿?cái)嗔逊植疾粩嗉由睢?/p>

4.5 溫泉水化學(xué)變化與地震的響應(yīng)

構(gòu)造活動(dòng)和地震活動(dòng)會(huì)打破地殼中流體原有的水-巖平衡狀態(tài)，使溫泉循環(huán)深度，熱儲(chǔ)溫度，深部熱水循環(huán)程度，水-巖相互作用改變，造成離子組分，氫氧同位素等發(fā)生變化（Skelton et al，2014）。通過(guò)龍門山地區(qū)溫泉的各離子濃度繪制主要離子相對(duì)變化率（（測(cè)量值-平均值）/平均值×100%）隨時(shí)間變化的折線圖（圖6），根據(jù)采集樣品離子含量的方差占平均值的百分比來(lái)判斷研究區(qū)異常值的臨界。

通過(guò)對(duì)比距離汶川MS8.0地震震中108 km以內(nèi)的GRG、JYG、SZZ、LX、WC溫泉在震后1個(gè)月及5個(gè)月采集的兩次溫泉水樣（表6），發(fā)現(xiàn)隨著時(shí)間的變化，溫泉中的主要離子組分（Na+、Cl-、SO2-4）及微量元素含量（Li+、Pb2+、As2+）呈下降趨勢(shì)，同時(shí)2008年6月至2010年6月采集水樣的δD和δ18O也回歸至大氣降水線附近。從δD-δ18O關(guān)系圖（圖2）發(fā)現(xiàn)，2008年6月和10月采集的GRG01、GRG02、JYG01、LX01、SZZ01、WC02的氧同位素偏離于大氣降水線的右下方。這可能是由于區(qū)域構(gòu)造活動(dòng)的改變使斷層封閉性改變，增加了溫泉水與礦物的接觸，使其水巖反應(yīng)增強(qiáng)，水化學(xué)離子含量增加。

2013年蘆山MS7.0地震后，在2013年4月及2014年8月采集了震中距為65 km以內(nèi)HSW和WJJ溫泉以及震中距為195 km以內(nèi)的SZZ、WC的水樣，其化學(xué)組分部分呈下降趨勢(shì)，但變化不明顯，也有少量元素呈增加的趨勢(shì)（表6）。這可能與蘆山地震后研究區(qū)多次發(fā)生5級(jí)以上余震，且在2014年11月發(fā)生了康定MS6.5地震（見(jiàn)圖6中紅色框），其深部熱水供應(yīng)變化不大有關(guān)。從δD-δ18O關(guān)系圖可以發(fā)現(xiàn)，2013年5月采集的JYG06、LFS06、PL06、SZZ06、WC06、WJJ06、ZBS06、ZGS06、ZLJ06氫氧同位素偏離于大氣降水線的左上方。因此，推測(cè)2013年蘆山地震時(shí)各溫泉離子濃度發(fā)生變化的原因?yàn)榈卣鹨鸬母骱畬娱g滲透性增強(qiáng)，地下水混合所致。

地震活動(dòng)還會(huì)改變應(yīng)力場(chǎng)變化，斷裂帶裂隙發(fā)生變化，引起滲透率發(fā)生變化，從而導(dǎo)致地震震中區(qū)周圍的流體行為改變（Zhou et al，2017）。周永勝等（2008）、宋娟和周永勝（2013）研究發(fā)現(xiàn)川滇地區(qū)的龍門山斷裂以及紅河斷裂帶發(fā)生內(nèi)部流體孔隙壓力周期性變化以及斷層帶脆-塑性轉(zhuǎn)化、裂縫張開(kāi)與愈合等現(xiàn)象。強(qiáng)震發(fā)生后，在巖石塑性變形、壓溶、靜態(tài)重結(jié)晶、水巖反應(yīng)等作用下，斷層內(nèi)部由地震形成的斷層和微裂縫逐漸被愈合，把斷層深部的流體密封。Luo等（2017）研究發(fā)現(xiàn)地下熱水的氧同位素向右偏移的原因?yàn)榈叵聼崴谏钛h(huán)過(guò)程中與富集δ18O圍巖發(fā)生水巖作用，伴隨氧同位素的交換，導(dǎo)致地下熱水δ18O的升高，這與2008年氫氧同位素的變化表現(xiàn)一致。2013年5月，氫氧同位素偏離大氣降水線的左上方，其原因?yàn)闇厝c深部來(lái)源的CO2相互溶解的過(guò)程中發(fā)生了δ18O同位素交換（Benavente et al，2016），這與2013年氫氧同位素的變化表現(xiàn)一致。

地震的孕育和發(fā)生過(guò)程中也會(huì)改變含水層及其周圍巖土環(huán)境，引起滲透率變化或水體混合，從而引起水文地球化學(xué)變化。地震后，斷層愈合，各離子濃度又大幅度回落，這可能是由地震后斷裂內(nèi)部滲透性減弱、深部流體供應(yīng)減少、水巖反應(yīng)程度衰減以及當(dāng)?shù)亟涤炅康脑龆?、大量淺層水混入所致。

4.6 溫泉成因模式

結(jié)合已有的地球物理資料和本文的研究成果建立模型，龍門山斷裂帶周圍3.2 km的高山為溫泉的補(bǔ)給區(qū)，大氣降水入滲沿?cái)嗔褞R聚補(bǔ)給地下水，地下水沿?cái)鄬悠扑閹Щ蛄严稁Ы?jīng)歷深循環(huán)，獲得來(lái)自深部熱源加熱，此時(shí)水體的溫度增加至123 ℃，循環(huán)深度增大至地下4 km左右，與周圍的巖石發(fā)生水巖反應(yīng)，并在水頭差驅(qū)動(dòng)下上升，在地形較低的山谷和河谷處出露地表（圖7）。而龍門山斷裂帶兩側(cè)中上地殼存在低速高導(dǎo)層，為龍門山中上地殼內(nèi)淺源地震的孕育提供了深部構(gòu)造條件。地震的孕育及發(fā)生可以改變地下應(yīng)力應(yīng)變狀態(tài)，觸發(fā)斷裂帶內(nèi)部含水層破碎，導(dǎo)致深部流體以及淺層冷水混入從而影響地下熱水的水巖反應(yīng)程度，并改變研究區(qū)溫泉的水化學(xué)組分及同位素特征。本文通過(guò)以上研究建立了龍門山地區(qū)溫泉水化學(xué)模型，結(jié)合研究區(qū)溫泉水化學(xué)變化對(duì)地震的響應(yīng)，為該地區(qū)的震情跟蹤、異常落實(shí)和地震預(yù)測(cè)提供基礎(chǔ)研究。

5 結(jié)論

本文對(duì)龍門山斷裂帶13個(gè)溫泉及觀測(cè)井的水文地球化學(xué)特征進(jìn)行分析研究，得到以下結(jié)論：

（1）龍門山斷裂帶溫泉水的主要來(lái)源為周圍3.2 km高山的大氣降水。

（2）根據(jù)Na-K-Mg三角圖、水化學(xué)組分特征以及循環(huán)深度綜合分析，龍門山斷裂帶溫泉的水化學(xué)類型沿龍門山斷裂帶從北到南、從西到東依次為重碳酸型、硫酸型、氯化物型，水巖反應(yīng)程度、水化學(xué)組分含量、循環(huán)深度沿?cái)嗔逊植疾粩嘣黾?。其中茂縣—汶川斷裂帶溫泉為未成熟水，水巖反應(yīng)程度較弱，循環(huán)深度達(dá)到0.7～2.5 km;北川—映秀斷裂、安縣—灌縣斷裂及廣元—大邑隱伏斷裂帶北部為未成熟水，水巖反應(yīng)程度較弱，循環(huán)深度達(dá)到1.1～1.5 km;中部為部分平衡水，水巖反應(yīng)程度較強(qiáng)，循環(huán)深度達(dá)到0.9～2.8 km;南部為平衡水及部分平衡水，水巖反應(yīng)程度強(qiáng)，循環(huán)深度達(dá)到4.0～4.4 km。

（3）對(duì)比汶川MS8.0地震、蘆山MS7.0地震震后幾個(gè)月內(nèi)溫泉水化學(xué)組分（Na+、Cl-、SO42-、Li、Sr、As、δ18O及δD）的變化，發(fā)現(xiàn)震中距200 km以內(nèi)的溫泉水化學(xué)組分呈下降趨勢(shì)，結(jié)合氫氧同位素的變化推測(cè)這可能是由震后斷層愈合、斷層滲透系數(shù)減弱、水巖反應(yīng)程度衰減所致。以上研究可為該地區(qū)震情跟蹤、異常落實(shí)和地震預(yù)測(cè)提供基礎(chǔ)。

感謝中國(guó)科學(xué)院地質(zhì)與地球物理研究所蘭州油氣資源研究中心與中國(guó)地震局地殼應(yīng)力研究所的相關(guān)工作人員在樣品分析過(guò)程中給予的幫助，在此表示衷心的感謝！

參考文獻(xiàn)：

曹琴，周訓(xùn)，張歡，等.2015.四川盆地臥龍河儲(chǔ)鹵構(gòu)造地下鹵水的水化學(xué)特征及成因[J].地質(zhì)通報(bào)，34（5）：990-997.

陳志.2014.汶川地震后川西地區(qū)溫泉水地球化學(xué)研究[D].合肥：中國(guó)科學(xué)技術(shù)大學(xué).

程萬(wàn)正，官致君，李軍.2013.對(duì)汶川8.0級(jí)地震前四川地區(qū)地下流體觀測(cè)異常的研究[J].四川地震，（2）：1-8.

崔希林.2014.四川盆地中西部溫泉成因模式研究[D].成都：成都理工大學(xué).

鄧起東，陳社發(fā)，趙小麟.1994.龍門山及其鄰區(qū)的構(gòu)造和地震活動(dòng)及動(dòng)力學(xué)[J].地震地質(zhì)，16（4）：389-403.

付碧宏，時(shí)丕龍，張之武.2008.四川汶川MS8.0大地震地表破裂帶的遙感影像解析[J].地質(zhì)學(xué)報(bào)，82（12）：1679-1687.

柯斌.2014.四川茂縣吉魚(yú)溫泉成因模式及開(kāi)發(fā)遠(yuǎn)景分析[D].成都：成都理工大學(xué).

李?yuàn)^生，趙國(guó)華，李勇，等.2015.龍門山地區(qū)水系發(fā)育特征及其對(duì)青藏高原東緣隆升的指示[J].地質(zhì)論評(píng)，61（2）：345-355.

李生紅.2013.川東褶皺帶中、下三疊統(tǒng)水文地質(zhì)特征研究[D].成都：成都理工大學(xué).

林耀庭，唐慶，熊淑君，等.1997.四川盆地鹵水的氫、氧同位素地球化學(xué)特征及其成因分類研究[J].地質(zhì)地球化學(xué)，（4）：20-26.

林耀庭.2001.四川盆地三疊紀(jì)鹵水成藏條件[J].化工礦產(chǎn)地質(zhì)，23（1）：19-24.

劉成龍，王廣才，史浙明，等.2020.云南硫磺洞溫泉水文地球化學(xué)特征和成因分析[J].地震研究，43（2）：278-286.

龐忠和，汪集旸，樊志成.1990.利用SiO2混合模型計(jì)算漳州地?zé)崽餆醿?chǔ)溫度[J].科學(xué)通報(bào)，（1）：57-59.

蒲小武，武銀，狄國(guó)榮，等.2013.甘肅清水溫泉井與臨夏井水溫在汶川地震前后的異常變化特征[J].地震研究，36（3）：269-274.

宋娟，周永勝.2013.斷層帶流體對(duì)斷層強(qiáng)度和強(qiáng)震孕育的影響[J].國(guó)際地震動(dòng)態(tài)，420（12）：5-16.

陶廣斌.2019.川西淡礦化溫泉地球化學(xué)特征及成因研究[D].成都：成都理工大學(xué).

王二七，孟慶任，陳智樑，等.2001.龍門山斷裂帶印支期左旋走滑運(yùn)動(dòng)及其大地構(gòu)造成因[J].地學(xué)前緣，8（2）：165-174.

王逸凌，穆文清，史浙明.2020.鮮水河斷裂南段溫泉水文地球化學(xué)特征[J].地震研究，43（2）：287-295.

王云，趙慈平，劉峰，等.2014.小江斷裂帶及鄰近地區(qū)溫泉地球化學(xué)特征與地震活動(dòng)關(guān)系研究[J].地震研究，37（2）：228-243.

王志，趙琳娜，張國(guó)平，等.2010.汶川地震災(zāi)區(qū)堰塞湖流域面雨量計(jì)算方法研究[J].氣象，36（6）：7-12.

徐錫偉，聞學(xué)澤，葉建青，等.2008.汶川MS8.0地震地表破裂帶及其發(fā)震構(gòu)造[J].地震地質(zhì)，30（3）：597-629.

薛鈞月.2009.龍門山構(gòu)造帶中—北段構(gòu)造流體地球化學(xué)特征及其與成藏關(guān)系的探討[D].成都：成都理工大學(xué).

晏銳，官致君，劉耀煒.2015.川西溫泉水溫觀測(cè)及其在蘆山MS7.0地震前的異?，F(xiàn)象[J].地震學(xué)報(bào)，37（2）：347-356.

楊耀，周曉成，官致君，等.2019.川西地下流體觀測(cè)井水文地球化學(xué)特征[J].礦物巖石地球化學(xué)通報(bào)，38（5）：966-976.

于津生.1997.中國(guó)同位素地球化學(xué)研究[M].北京：科學(xué)出版社.

張?jiān)罉颍瑮钷r(nóng)，施煒，等.2008.青藏高原東緣新構(gòu)造及其對(duì)汶川地震的控制作用[J].地質(zhì)學(xué)報(bào)，82（12）：1668-1678.

周曉成，王萬(wàn)麗，李立武，等.2020.金沙江—紅河斷裂帶溫泉?dú)怏w地球化學(xué)特征[J].巖石學(xué)報(bào)，36（7）：2197-2214.

周訓(xùn)，曹琴，尹菲，等.2015.四川盆地東部高褶帶三疊系地層鹵水和溫泉的地球化學(xué)特征及成因[J].地質(zhì)學(xué)報(bào)，89（11）：1908-1920.

周訓(xùn).1993.四川盆地龍女寺儲(chǔ)鹵構(gòu)造深層地下鹵水的水文地球化學(xué)特征及成因[J].現(xiàn)代地質(zhì)，7（1）：83-92.

周永勝，何昌榮，楊曉松.2008.中地殼韌性剪切帶中的水與變形機(jī)制[J].中國(guó)科學(xué)：地球科學(xué)，38（7）：819-832.

Benavente O，Tassi F，Reich M，et al.2016.Chemical and isotopic features of cold and thermal fluids discharged in the Southern Volcanic Zone between 32.5°S and 36°S：Insights into the physical and chemical processes controlling fluid geochemistry in geothermal systems of Central Chile[J].Chemical Geology，420（1）：97-113.

Chen Z，Zhou X C，Du J G，et al.2015.Hydrochemical characteristics of hot spring waters in the Kangding district related to the Lushan MS= 7.0 earthquake in Sichuan，China[J].Natural Hazards & Earth System Ences Discussions，2（6）：1149-1156.

Clark M K，Royden L H.2000.Topographic ooze：Building the eastern margin of Tibet by lower crustal flow[J].Geology，28（8）：703-706.

Craig H.1961.Isotopic Variations in Meteoric Waters[J].Science，133（3465）：1702-1703.

Fournier R O，Rowe J J.1966.The deposition of silica in hot springs[J].29（1）：585-587.

Giggenbach W F.1988.Geothermal solute equilibria.Derivation of Na-K-Mg-Ca geoindicators[J].Geochimica Et Cosmochimica Acta，52（12）：2749-2765.

Godard V，Pik R，Cattin R，et al.2009.Late Cenozoic evolution of the central Longmen Shan，eastern Tibet：Insight from（U-Th）/He thermochronometry[J].Tectonics，28（5）：TC5009.

Kouketsu Y，Shimizu I，Wang Y，et al.2017.Raman spectra of carbonaceous materials in a fault zone in the Longmenshan thrust belt，China;comparisons with those of sedimentary and metamorphic rocks[J].Tectonophysics，699（1）：129-145.

Lambrakis N，Kallergis G.2005.Contribution to the study of Greek thermal springs：hydrogeological and hydrochemical characteristics and origin of thermal waters[J].Hydrogeology Journal，13（3）：506-521.

Lin X，Dreger D，Ge H，et al.2018.Spatial and Temporal Variations in the Moment Tensor Solutions of the 2008 Wenchuan Earthquake Aftershocks and Their Tectonic Implications[J].Tectonics，37（3）：989-1005.

Liu Y，Ren H，Wang B.2009.Application of environmental isotopes and tracer techniques to seismic subsurface fluids[J].Earth Science Frontiers，16（1）：369-377.

Luo J，Pang Z H，Kong Y K，et al.2017.Geothermal potential evaluation and development prioritization based on geochemistry of geothermal waters from Kangding area，western Sichuan，China[J].Environmental Earth Ences，76（9）：343.

Mokadem N，Demdoum A，Hamed Y，et al.2016.Hydrogeochemical and stable isotope data of groundwater of a multi-aquifer system：Northern Gafsa basin- Central Tunisia[J].Journal of African Earthences，114（1）：174-191.

Pérez N M，Hernández P A，Igarashi G，et al.2008.Searching and detecting earthquake geochemical precursors in CO2-rich groundwaters from Galicia，Spain[J].Geochemical Journal，42（1）：75-83.

Okuyama Y，F(xiàn)unatsu T，F(xiàn)ujii T，et al.2016.Mid-crustal fluid related to the Matsushiro earthquake swarm（1965-1967）in northern Central Japan：geochemical reproduction，Tectonophysics，679（3）：61-72.

Shi Z，Zhang H，Wang G.2020.Groundwater trace elements change induced by MS5.0 earthquake in Yunnan[J].Journal of Hydrology，581：124424.

Skelton A，Andren M，Kristmannsdottir H，et al.2014.Changes in groundwater chemistry before two consecutive earthquakes in Iceland[J].Nature Geoscience，7（10）：752-756.

Soto-Jiménez M F，Páez-Osuna F.2001.Distribution and Normalization of Heavy Metal Concentrations in Mangrove and Lagoonal Sediments from Mazatlán Harbor（SE Gulf of California）[J].Estuarine Coastal & Shelf Ence，53（3）：259-274.

Tsunogai U，Wakita H.1995.Precursory Chemical Changes in Ground Water：Kobe Earthquake，Japan[J].Science，5220（269）：61-63.

Walraevens K，Bakundukize C，Mtoni Y E，et al.2018.Understanding the hydrogeochemical evolution of groundwater in Precambrian basement aquifers：A case study of Bugesera region in Burundi[J].Journal of Geochemical Exploration，188（3）：24-42.

Xu Q，Hoke G D，Jing L Z，et al.2015.Stable isotopes of surface water across the Longmenshan margin of the eastern Tibetan Plateau[J].Geochemistry Geophysics Geosystems，15（8）：3416-3429.

Xu Z，Zhen F，Li SX，et al.2019.The 2018 MS5.9 Mojiang Earthquake：Source model and intensity based on near-field seismic recordings[J].Earth and Planetary Physics，3（3）：88-101.

Zhang Y F，Tan H B，Zhang W J，et al.2015.A New Geochemical Perspective on Hydrochemical Evolution of The Tibetan Geothermal System[J].Geochemistry International，53（12）：1090-1106.

Zhang Y F，Tan H B，Zhang W J，et al.2016.Geochemical constraint on origin and evolution of solutes in geothermal springs in western Yunnan，China[J].Chemie der Erde-Geochemistry，76（1）：63-75.

Zhou X，Jiang C，Zhao J，et al.2015.Occurrence and resource evaluation of the subsurface high-K brines in the Pingluoba brine-bearing structure in western Sichuan Basin[J].Environmental Earth Sciences，73（12）：8565-8574.

Zhou X C，Li C，Ju X M，et al.2010.Origin of Subsurface Brines in the Sichuan Basin[J].Groundwater，35（1）：53-58.

Zhou X C，Liu L，Chen Z，et al.2017.Gas geochemistry of the hot spring in the Litang fault zone，Southeast Tibetan Plateau[J].Applied Geochemistry，79（1）：17-26.

Hydrogeochemical Characteristics of the Hot Springsin the Longmenshan Fault Zone

YAN Yucong1，LIU Fengli1，GUO Lishuang2，ZHOU Xiaocheng1，OUYANG Shupei1，LI Jingchao1，WANG Wanli1

（1.Institute of Earthquake Forecasting，China Earthquake Administration，Beijing 100036，China）

（2.National Institute of Natural Hazards，Ministry of Emergency Management，Key Laboratory ofCrustal Dynamics，China Earthquake Administration，Beijing 100085，China）

Abstract

Hydrogeochemical Monitoring in the regions with strong earthquake activities would facilitate the obtaining of the geochemical information about deep fluids correlated to seismic anomalies.In this paper the hydrogeochemical characteristics of 13 hot springs（wells）in the Longmenshan Fault Zone were investigated，so as to establish the hydrogeochemical cycle model of the hot springs and reveal the relationship between hydrogeochemical variations and seismic activity in the Fault Zone.The results show that：①The hot springs in the Longmenshan Fault Zone were recharged by the meteoric water at elevations ranging from 0.8 to 3.2 km.②The chemical types of the hot springs from north to south，and from west to east along the Longmenshan Fault Zone were respectively carbonic acid，sulfuric acid and chloride.Likewise，the circulation depth of spring water deepened;the water-rock reaction intensified;and enrichment factors of the trace element of the hot springs continuously increased.③After the Wenchuan MS8.0 earthquake and the Lushan MS7.0 earthquake，the dissolved hydrochemical composition got declining in hot springs within the epicenter of 200 kilometers，which could be the result of the permeability decrease and the attenuation of water-rock reaction inside the fault in the healing process of the Longmenshan Fault after earthquake，

Keywords：hot spring;hydrogeochemistry;earthquake;the Longmenshan Fault Zone

收稿日期：2021-01-11.

基金項(xiàng)目：國(guó)家重點(diǎn)研發(fā)計(jì)劃（2017YFC1500501）、國(guó)家自然科學(xué)基金面上項(xiàng)目（41673106）和中國(guó)地震局地震預(yù)測(cè)研究所基本科研業(yè)務(wù)費(fèi)（2017IES010205，2016IES010304，2018IEF010104，2020IEF0604，2020IEF0703）聯(lián)合資助.

第一作者簡(jiǎn)介：顏玉聰（1996-），女，在讀碩士研究生，主要從事地震地球化學(xué)研究.E-mail：yanyucong2020@163.com.

通訊作者簡(jiǎn)介：郭麗爽（1983-），高級(jí)工程師，主要從事地下流體地球化學(xué)方面的研究.E-mail：guolsh02@163.com.