地下水井水位及化學(xué)組分的同震差異響應(yīng)特征分析
——以青白江地震為例

顧鴻宇，王東輝，李勝偉，鄭萬模，劉 港，向元英，李 丹，陳能德

（1.中國(guó)地質(zhì)調(diào)查局成都地質(zhì)調(diào)查中心，四川 成都 610081；2.中國(guó)地質(zhì)調(diào)查局西安地質(zhì)調(diào)查中心，陜西 西安 710054；3.四川省工程咨詢研究院，四川 成都 610016；4.四川省開元集團(tuán)有限公司，四川 成都 610081）

近年來，地下水系統(tǒng)對(duì)地震的響應(yīng)機(jī)理引發(fā)了眾多學(xué)者的廣泛關(guān)注。地下水系統(tǒng)對(duì)地震的響應(yīng)，從地震時(shí)間歷程劃分，包括震前響應(yīng)、同震響應(yīng)和震后響應(yīng)3 個(gè)階段。其響應(yīng)方式多種多樣，包括地下水水位、水溫、水化學(xué)組分及地下水流量的改變等[1?2]。

由于地下水水位監(jiān)測(cè)簡(jiǎn)便，前人研究普遍關(guān)注地下水水位對(duì)地震的響應(yīng)規(guī)律[3]。研究形式可歸結(jié)為2 類：（1）單一監(jiān)測(cè)井對(duì)多次地震的監(jiān)測(cè)，此方法有利于消除非地震引發(fā)的水位變動(dòng)以及不同水文地質(zhì)條件的干擾，常用于分析地震引起的應(yīng)力擾動(dòng)[4?5]；（2）多監(jiān)測(cè)井對(duì)同一地震響應(yīng)規(guī)律的監(jiān)測(cè)，監(jiān)測(cè)井的空間分布有利于研究地下水位在不同水文地質(zhì)條件及構(gòu)造條件中對(duì)地震的差異響應(yīng)[6]。由于地震能改變地球的地殼結(jié)構(gòu)，從而導(dǎo)致巖體內(nèi)孔隙壓力變化以及含水層介質(zhì)特性（如滲透率），進(jìn)而產(chǎn)生不同形式的同震水位響應(yīng)。地下水水位的響應(yīng)通常有階躍式、漸變式和振蕩式3 種典型形式以及不同的組合形式，其空間分布與地震震級(jí)及震中距有關(guān)[6]。

此外，由于地下熱流分布廣泛，巖體的熱導(dǎo)率在空間和時(shí)間上變異性較小等特點(diǎn)，地下水溫度也可用于追蹤地震對(duì)地下水系統(tǒng)的影響[7]。溫度剖面和溫度時(shí)間序列數(shù)據(jù)常被用于解釋地下水系統(tǒng)的地震響應(yīng)。地層淺部的溫度通常反映氣候周期的變化規(guī)律（日變化、季節(jié)變化和年變化）。不同深度之間的溫度差異及不同的時(shí)間變化規(guī)律可以造成熱流擴(kuò)散。由于地震改變了巖體滲透性，地下水在流動(dòng)過程中會(huì)對(duì)溫度場(chǎng)進(jìn)行重構(gòu)，改變溫度數(shù)據(jù)的時(shí)間序列特征。因此，對(duì)比地震前后的溫度時(shí)間序列規(guī)律常用于判斷地下水系統(tǒng)特征變化。

地下水水化學(xué)組分的地震響應(yīng)規(guī)律更為復(fù)雜，地震引發(fā)的物理作用（如滲透性改變）和水-巖反應(yīng)均會(huì)造成水化學(xué)組分的改變[8?10]。由于地震可能導(dǎo)致巖體的塑性變形，不同含水層間的隔水層被破壞，不同化學(xué)特征的水源混合，甚至發(fā)生化學(xué)反應(yīng)。同時(shí)，地震對(duì)流體壓力的改變可能導(dǎo)致某些礦物發(fā)生溶解或結(jié)晶。因此，離子本身在場(chǎng)地環(huán)境中的化學(xué)穩(wěn)定性以及濃度是其地震敏感性的重要影響因素。

目前，水位、水溫的聯(lián)合監(jiān)測(cè)較為普遍，而地下水水化學(xué)組分的時(shí)間序列數(shù)據(jù)還較為缺乏。一方面，水化學(xué)數(shù)據(jù)大多是定期采樣檢測(cè)，其采樣周期較長(zhǎng)，難以捕捉到地震期間的短時(shí)響應(yīng)特征。另一方面，由于地下水系統(tǒng)的各個(gè)指標(biāo)對(duì)地震響應(yīng)的敏感性不同，單一指標(biāo)可能會(huì)忽略重要的信息，從而導(dǎo)致對(duì)機(jī)理認(rèn)識(shí)出現(xiàn)偏差。本文采用了水位、水溫、溶解性總固體（TDS）、pH 及氧化還原電位（Eh）多參數(shù)聯(lián)合實(shí)時(shí)監(jiān)測(cè)，獲得了用較高頻率采樣的時(shí)間序列數(shù)據(jù)，并探討了龍泉山地區(qū)地下水系統(tǒng)對(duì)近場(chǎng)地震的差異響應(yīng)機(jī)理。

1 研究區(qū)概況

1.1 地質(zhì)概況

研究區(qū)位于成都市龍泉山（圖1），地貌屬構(gòu)造剝蝕淺切脊?fàn)畹蜕降孛?，地?shì)呈西北高東南低，高差約100 m。龍泉山斷裂總體走向NE20°～30°，斷裂帶北端位于德陽以北，向南經(jīng)過中江縣西側(cè)、金堂鎮(zhèn)、龍泉驛鎮(zhèn)、久隆場(chǎng)鎮(zhèn)和仁壽鎮(zhèn)西側(cè)，直達(dá)樂山新橋鎮(zhèn)附近，全長(zhǎng)大于230 km。研究區(qū)內(nèi)密切相關(guān)的構(gòu)造有金龍寺背斜和四方山斷裂，位于龍泉山背斜北西翼，西距龍泉山斷層約2.5 km[11]。

圖1 研究區(qū)構(gòu)造地質(zhì)概況Fig.1 Tectonic geology in the study area

金龍寺背斜呈北東走向，全長(zhǎng)約16 km，兩翼受四方山斷層的影響呈不對(duì)稱形態(tài)，西翼傾角8°～11°，東翼傾角4°。研究區(qū)內(nèi)出露的地層均為蓬萊鎮(zhèn)組地層，局部出現(xiàn)巖層倒轉(zhuǎn)。

四方山斷裂走向?yàn)楸北蔽飨颍▓D2），地表出露約6 km。斷層兩盤均為蓬萊鎮(zhèn)組地層，兩側(cè)巖層產(chǎn)狀陡峭，局部倒轉(zhuǎn)。斷層帶巖體破碎，裂隙和劈理發(fā)育，此斷層為傾向北西的高角度逆斷層。斷層附近有2 組節(jié)理：25°∠38°，節(jié)理密度7 條/m；2°∠53°，節(jié)理密度0.4 條/m。井位附近地表測(cè)得2 組節(jié)理：15°～30°∠80°～89°和45°～50°∠70°～89°。

圖2 研究區(qū)水文地質(zhì)剖面圖Fig.2 Hydrogeological profile of the study area

研究區(qū)內(nèi)地下水補(bǔ)給主要為大氣降水，地表無自然水體下滲。含水基巖主要為蓬萊鎮(zhèn)組的互層砂泥巖，地下水總體沿順層裂隙自南向北流動(dòng)。

1.2 監(jiān)測(cè)井概況

ZK1 位于垃圾填埋場(chǎng)大壩北滲濾池下游溝谷600 m左側(cè)斜坡上，鉆孔孔口標(biāo)高為576.2 m，深度90.7 m。表層0～8 m 用直徑110 mm 鋼管護(hù)壁，并進(jìn)行水泥永久止水，8 m 以下采用花管。ZK6 位于長(zhǎng)安垃圾場(chǎng)北西20°方位1.2 km 處，鉆孔孔口標(biāo)高為555.5 m，深度61.2 m。表層0～5.6 m 用直徑110 mm 鋼管護(hù)壁，并進(jìn)行水泥永久止水，5.6 m 以下采用花管。ZK6 位于ZK1西北方向下游位置，直線距離約700 m（圖3）。

圖3 監(jiān)測(cè)井位置Fig.3 Location of two monitoring wells

根據(jù)抽水試驗(yàn)結(jié)果，ZK1 在井深7.2～63.7 m 的滲透系數(shù)為9.491×10?5～1.227×10?4cm/s，井深63.7～90.7 m 的滲透系數(shù)為1.632×10?4～1.910×10?4cm/s。ZK6井深7.0～61.2 m 的滲透系數(shù)為4.815×10?4～6.644×10?4cm/s。ZK1 水溫為18.87～20.36 ℃，ZK6 水溫為18.82～19.30 ℃。

根據(jù)鉆孔巖芯、物探電阻率、聲波測(cè)井并結(jié)合水位觀測(cè)等資料，綜合確定了水文地質(zhì)結(jié)構(gòu)，揭露2 層含水層（圖4）：第一含水層，在斜坡區(qū)埋深32.5～46.2 m，在山谷地區(qū)頂板埋深為12.0～21.0 m，底板埋深為31.0～32.2 m，厚度為11.0～16.0 m，巖性以粉砂質(zhì)泥巖為主，中風(fēng)化，表面可見灰綠色粉砂質(zhì)斑點(diǎn)及條帶，發(fā)育5 條斜傾角節(jié)理，傾角為30°～68°，裂面起伏粗糙；發(fā)育一組陡傾角節(jié)理，傾角為85°，裂面起伏粗糙。第二含水層，頂板41.6～48.9 m，底板深度51.6 m 以下（未揭穿），含水巖組（段）厚度3.0～5.5 m，巖性以中-強(qiáng)風(fēng)化粉砂質(zhì)泥巖為主，節(jié)理裂隙較為發(fā)育，巖體破碎，為風(fēng)化帶裂隙孔隙水。

圖4 鉆孔ZK1、ZK7 巖性柱狀圖及含水層位置Fig.4 Lithologic log of wells ZK1 and ZK6 and the depth of two identified aquifers

ZK1 取樣深度為50.0 m，位于第二含水層，其主要化學(xué)指標(biāo)及TDS 值分別為15，16.88，<0.04，431.2 mg/L，水化學(xué)類型為HCO3—Na·Ca·Mg型。ZK6 取樣深度為19.0 m，位于第一含水層，其主要化學(xué)指標(biāo)及TDS 值分別為1 274，68.72，<0.04，1 791.0 mg/L，水化學(xué)類型為HCO3·Cl—Na·Ca·Mg型。結(jié)果表明垃圾滲濾液主要影響第一含水層，對(duì)第二含水層影響較小。

1.3 監(jiān)測(cè)設(shè)備及方法

為了獲得地下水系統(tǒng)的變化規(guī)律，2019年10月25日—2020年7月16日對(duì)2 口井中地下水進(jìn)行自動(dòng)監(jiān)測(cè)，儀器型號(hào)為In-Situ 公司Aqua TROLL@500 型水位水質(zhì)自動(dòng)監(jiān)測(cè)儀，見圖5。監(jiān)測(cè)參數(shù)包括：地下水水位，溫度，溶解性總固體（TDS）、pH、氧化還原電位（Eh）。數(shù)據(jù)采集間隔為12 h，即0 點(diǎn)和12 點(diǎn)整采集數(shù)據(jù)，數(shù)據(jù)發(fā)射頻率為1 次/d。各傳感器在布設(shè)前采用標(biāo)準(zhǔn)液進(jìn)行標(biāo)定。

圖5 監(jiān)測(cè)井及儀器設(shè)備Fig.5 Photo of the monitoring well and equipment

2 結(jié)果

2.1 同震差異響應(yīng)特征

2020年2月3日0 時(shí)5 分，成都市青白江發(fā)生5.1 級(jí)地震，引發(fā)了龍泉山地下水監(jiān)測(cè)井ZK1 和ZK6的水位、水化學(xué)組分及其他參數(shù)的同震異常。兩口監(jiān)測(cè)井距離震中約13 km（圖1）。兩口監(jiān)測(cè)井雖相距僅約700 m，但其響應(yīng)特征明顯不同。差異性響應(yīng)為研究地下水系統(tǒng)的地震響應(yīng)機(jī)理提供了更詳實(shí)的依據(jù)。

（1）水位

震前，兩口井水位未見異常波動(dòng)。地震發(fā)生時(shí)，兩口井的地下水位都發(fā)生了同震階躍形式的上升，但上升幅度不同。ZK1 水位上升了0.80 m，見圖6（a），ZK6 水位上升了0.25 m，見圖7（a）。兩口井的含水層厚度差異可能是導(dǎo)致補(bǔ)給量差異的原因之一。

圖6 ZK1 地下水各類參數(shù)對(duì)青白江地震的響應(yīng)特征Fig.6 Coseismic characteristics of groundwater parameters in ZK1 to the Qingbaijiang earthquake

圖7 ZK6 地下水各類參數(shù)對(duì)青白江地震的響應(yīng)特征Fig.7 Coseismic characteristics of groundwater parameters in ZK6 to the Qingbaijiang earthquake

（2）水化學(xué)組分

Cl?在地下水中化學(xué)性質(zhì)較穩(wěn)定，因此可作為示蹤離子。震前兩口井的Cl?未見異常波動(dòng)，以一定的速率不斷增大。地震發(fā)生時(shí)，ZK1 中Cl?濃度出現(xiàn)了同震階躍形式的下降，濃度下降幅度達(dá)196 mg/L，見圖6（b）。ZK6 中Cl?濃度則出現(xiàn)了相反的變化趨勢(shì)，Cl?濃度同震階躍上升了752 mg/L，見圖7（b）。兩者的差異變化反映了地震時(shí)地下水補(bǔ)給來源發(fā)生改變。

震前，兩口井中的pH 都未發(fā)生異常波動(dòng)。地震發(fā)生時(shí)，ZK1 中的pH 有小幅度升高，但并未偏離其長(zhǎng)期變化趨勢(shì)，見圖6（e）。ZK6 中的pH 出現(xiàn)了較大幅度的增加，從7.27 增加到7.35，見圖7（e）。

ZK1 的Eh為負(fù)值，而ZK6 的Eh為正值，這表明兩口井所處的氧化還原環(huán)境差異明顯，地下水的補(bǔ)給來源不同。地震發(fā)生時(shí)，ZK1 中的Eh變化規(guī)律和變化規(guī)律相似，出現(xiàn)了周期的短暫改變，見圖6（f）。ZK6中的Eh出現(xiàn)了階躍式的下降，從481 mV 下降到404 mV，見圖7（f）。

震前，兩口井的TDS 均保持原有的趨勢(shì)性變化。地震發(fā)生時(shí)，ZK1 中TDS 出現(xiàn)了明顯的下降，從1 435.5 mg/L下降到1 384.9 mg/L，見圖6（g）。ZK6 中TDS 則出現(xiàn)了大幅增長(zhǎng)現(xiàn)象，從1 888.7 mg/L 增加到2 247.3 mg/L，見圖7（g）。地下水補(bǔ)給來源的不同以及不同含水層在地震時(shí)的補(bǔ)給比例差異是導(dǎo)致兩口井TDS 差異變化的因素之一。

地震對(duì)兩口井的溫度并無影響，兩口井溫度都保持在18.7 ℃，見圖6（h）、圖7（h），說明向井補(bǔ)給的含水層無明顯溫度差異，且地震并未導(dǎo)致深部高溫地下水的補(bǔ)給。

2.2 震后差異響應(yīng)特征

兩口觀測(cè)井的震后差異主要表現(xiàn)在各參數(shù)的恢復(fù)特征不同（表1）。ZK1 中的水位和Cl?震后呈現(xiàn)出平移現(xiàn)象，即震后并未出現(xiàn)明顯的恢復(fù)過程，地震后各參數(shù)的變化速率（曲線的斜率）與震前基本一致。TDS 震后也未出現(xiàn)明顯的恢復(fù)過程，但其變化趨勢(shì)相反，TDS 值由逐漸增加變?yōu)榫徛陆?。ZK6 中Cl?、、pH 及TDS 均出現(xiàn)了明顯的恢復(fù)過程，恢復(fù)周期約40 d。水位的恢復(fù)過程較快，周期約15 d，且恢復(fù)曲線特征呈突變的形狀。ZK6 中在震后無恢復(fù)過程，可能是部分不連通裂隙中高濃度地下水在地震作用下短暫開啟并排入井內(nèi)，從而掩蓋了的恢復(fù)過程。

表1 ZK1 和ZK6 水位及水化學(xué)參數(shù)同震變化表Table 1 Coseismic response of water levels and chemical components in ZK1and ZK6

3 地下水系統(tǒng)地震響應(yīng)機(jī)理探討

3.1 地震敏感性分析

地震能量密度代表著地震施加在單位體積介質(zhì)上的最大能量，關(guān)鍵參數(shù)包括震中距以及震級(jí)。Wang 等[12]根據(jù)全球地震數(shù)據(jù)統(tǒng)計(jì)了同震響應(yīng)大小和震中距以及震級(jí)的相關(guān)性，研究表明10?4J/m3是激發(fā)同震響應(yīng)所需的最小能量密度閾值。本文中監(jiān)測(cè)井與震中的距離約13 km，小于單位斷裂長(zhǎng)度。因此，以震源與監(jiān)測(cè)井距離判斷，青白江地震屬于近場(chǎng)地震。根據(jù)Wang 等[13]提出的能量密度計(jì)算公式：

logr=0.48?M?0.33?loge?1.4

式中：r—震中到監(jiān)測(cè)井距離/km；

e—能量密度/（J·m?3）；

M—地震震級(jí)。

本次地震在監(jiān)測(cè)井附近產(chǎn)生的能量密度為4.77×10?1J/m3（圖8），遠(yuǎn)高于同震地下水位響應(yīng)的閾值10?4J/m3。兩口井的位置雖相距較近，但水位變幅卻相差達(dá)到0.55 m。這說明能量密度只能判斷ZK1 對(duì)近場(chǎng)地震的敏感性高于ZK6，可能是由于兩口井結(jié)構(gòu)、含水層厚度、巖體結(jié)構(gòu)以及地形差異等造成的。

圖8 青白江地震在監(jiān)測(cè)井位置的能量密度（修改自文獻(xiàn)[14]）Fig.8 Seismic energy at the monitoring wells triggered by the Qingbaijiang earthquake（modified from Ref.[14]）

3.2 地下水化學(xué)組分及物理參數(shù)分析

離子濃度的變化受多種因素影響，尚無統(tǒng)一的機(jī)理能對(duì)所有現(xiàn)象進(jìn)行解釋。研究表明，地震引起的深部高溫地下水與淺部低溫地下水混合將引起化學(xué)再平衡，最終改變離子濃度及物理參數(shù)。青白江地震前后兩口井的溫度均未發(fā)生變化，這表明此次地震并未連通深部和淺部的含水層。因此，化學(xué)反應(yīng)導(dǎo)致的離子濃度變化不能解釋觀測(cè)數(shù)據(jù)的異常。此外，從觀測(cè)數(shù)據(jù)可以看出，ZK1 和ZK6 的離子濃度耦合變化規(guī)律截然不同，即ZK1 中各離子濃度及參數(shù)變化規(guī)律較為獨(dú)立，而ZK6 中的離子濃度及參數(shù)變化協(xié)同性較好，這表明不同井中離子濃度變化的控制因素存在差異。

地震導(dǎo)致不同水源混合是地下水化學(xué)組分異常的另一可能的原因。從圖9 可知，大氣壓力周期約為128 h，而變化周期依然約為36 h，說明其周期變化并不受大氣壓力控制。推測(cè)可知，影響硝酸根濃度波動(dòng)的因素可能為固體潮，周期性的固體潮使得裂隙反復(fù)地張開和閉合，從而導(dǎo)致不連通裂隙中高濃度間歇性地與井內(nèi)地下水混合，造成濃度周期性波動(dòng)。此外，ZK1 水位及Cl?濃度出現(xiàn)同震階躍變化，其震后未出現(xiàn)緩慢的恢復(fù)過程。這種現(xiàn)象表明地震改變了滲透性。由于同震時(shí)向井排泄的流體壓力遠(yuǎn)高于震后井水流入含水層的壓力，加之裂隙開度在震后減小，使得井水回流到裂隙的體積減少，回流速率急劇降低，最終導(dǎo)致水位及Cl?濃度在震后恢復(fù)過程不顯著。由于鉆孔揭露了兩層含水層，必定存在著其中一層含水層地下水發(fā)生向井內(nèi)的大量補(bǔ)給（圖10）。伴隨著水位上升，Cl?濃度和TDS 顯著下降，說明向井補(bǔ)給的地下水中并未受到垃圾填埋場(chǎng)的污染。同時(shí)，Eh值未發(fā)生同震的顯著異常波動(dòng)，進(jìn)一步說明該主要補(bǔ)給層控制著井內(nèi)Eh，Eh為負(fù)值說明該含水層較為封閉且與地表連通性差。兩層含水層以泥巖作為中間隔水層，隔水性較好，第二含水層理論上應(yīng)具有更低的Eh值。因此，可以推斷地震引發(fā)的深部第二含水層大量向井補(bǔ)給和混合是導(dǎo)致ZK1 中化學(xué)組分同震異常的主要原因。同時(shí)，混合作用并未導(dǎo)致濃度和濃度異常，表明和主要來源于第二含水層，其補(bǔ)給來源主要是大氣降水。

圖9 大氣壓力時(shí)頻譜圖Fig.9 Magnitude Scalogram of the atmospheric pressure

ZK6 中化學(xué)組分的異常較為一致，呈現(xiàn)階躍式變化。除Eh出現(xiàn)階躍下降外，其余組分均出現(xiàn)階躍上升。隨著水位的上升，pH、Cl?濃度和TDS 急劇增加，表明淺部第一含水層（受到垃圾滲濾液的污染）發(fā)生了向井的強(qiáng)烈補(bǔ)給。同時(shí)，和的階躍上升表明深部第二含水層也發(fā)生了向井的補(bǔ)給，但補(bǔ)給比例比第一含水層?。▓D10）。從Cl?濃度和TDS 震后恢復(fù)曲線可以看出，震前和震后恢復(fù)以后其濃度增長(zhǎng)速率未出現(xiàn)顯著變化，表明地震并未改變兩層含水層的滲透性，只是在同震及震后恢復(fù)期間出現(xiàn)了滲透性的短暫變化。井中化學(xué)組分的同震變化是由于地震溝通了含水層中原本不連通的裂隙，使含高濃度離子的地下水補(bǔ)給到井中并與低濃度地下水（賦存于連通裂隙中）混合。此外，ZK6 中并未出現(xiàn)類似于ZK1 中濃度周期性波動(dòng)的現(xiàn)象，這可能與ZK6 附近含水層巖體結(jié)構(gòu)變化有關(guān)。這也反映了巖體結(jié)構(gòu)的差異影響地下水系統(tǒng)的同震響應(yīng)特征。

圖10 ZK1 和ZK6 同震地下水響應(yīng)概念模型Fig.10 Schematic diagrams showing the coseismic differential response of groundwater in wells ZK1 and ZK6

3.3 地下水水位響應(yīng)機(jī)理

根據(jù)前人研究[13?14]，引起地下水位變化的機(jī)理大致可分為4 類：（1）地殼彈性應(yīng)變引起孔隙壓力振蕩；（2）巖土體不排水體積改變，包括不排水剪脹和不排水固結(jié)；（3）地震波引起滲透性增強(qiáng)，包括裂隙張開、新裂隙的形成以及驅(qū)替填充物，通常引起地下水位緩慢且持續(xù)地改變；（4）地震造成井筒-含水層結(jié)構(gòu)共振。

巖土體不排水體積改變被認(rèn)為是引起近場(chǎng)地下水水位階躍變化的主要因素[13]。水位波動(dòng)正比于含水層的體應(yīng)變，但兩者反相，即體應(yīng)變膨脹時(shí)水位下降，壓縮時(shí)水位上升[15]。研究的兩口井均位于砂泥巖互層的巖體中，巖體體積改變首先是結(jié)構(gòu)面的閉合或張開及巖塊本身的變形，其次是巖塊間的滑移。兩者都可能發(fā)生彈性變形或者塑性變形。若地震僅引起彈性變形，則在震前加載過程中，震前一段時(shí)間水位增速應(yīng)隨結(jié)構(gòu)面開度的增加而增加，反之隨開度減小而減小。與此同時(shí)，由于地下水的穩(wěn)定補(bǔ)給，滲流通道的變化會(huì)造成地下水壓力緩慢變大。地震發(fā)生時(shí)，巖體卸荷造成結(jié)構(gòu)面的彈性回跳，地下水在較高的壓力下向井排泄，將造成地下水位的階躍變化。由于彈性變形不改變巖體震前的滲透性，震后地下水位應(yīng)逐漸恢復(fù)到震前水平，或震前水位變化的趨勢(shì)線上。在這個(gè)過程中，恢復(fù)曲線呈指數(shù)型，這是由于在巖體恢復(fù)到震前應(yīng)力狀態(tài)過程中，滲透性也在逐漸恢復(fù)。若地震引起的應(yīng)變超過10?4，將導(dǎo)致巖體的塑性變形，巖體沿裂隙錯(cuò)動(dòng)或產(chǎn)生新的裂隙，將導(dǎo)致滲透性的永久改變[16]。

ZK1 水位出現(xiàn)了階躍上升，震后并未恢復(fù)，而是持續(xù)上升。這說明巖體在不排水條件下發(fā)生的塑性變形導(dǎo)致了裂隙滲透的連通性或者開度的減小，從而降低了滲透性[17?18]。此外，從水位-水化學(xué)組分耦合變化規(guī)律可以看出，第一含水層的滲透性未出現(xiàn)顯著變化，因?yàn)檎鹎昂驼鸷驝l?濃度增長(zhǎng)速率基本一致。第二含水層的滲透性出現(xiàn)了顯著降低，井水在震后無法快速回流到巖體裂隙中，導(dǎo)致了水位無法回落并恢復(fù)到震前水平。由于第二含水層中的濃度波動(dòng)周期并未受到地震影響，可以推斷第二含水層的滲透性降低并不是裂隙連通性減小引起，而是由于裂隙開度的減小引起（圖10）。

ZK6 水位出現(xiàn)了階躍上升，但在約15 d 后恢復(fù)到了震前水平，這種指數(shù)型的恢復(fù)曲線反映了震后滲透性的緩慢恢復(fù)。巖體在地震時(shí)發(fā)生了彈性變形，增加了裂隙的連通程度以及開度，導(dǎo)致了兩層含水層中地下水向井的排泄。地震后巖體內(nèi)應(yīng)力水平遠(yuǎn)低于同震時(shí)巖體內(nèi)的應(yīng)力水平。因此，震后較慢的卸載速率導(dǎo)致了裂隙的回彈速度減慢，進(jìn)而使得地下水位的恢復(fù)出現(xiàn)了滯后的現(xiàn)象（圖10）。

4 結(jié)論

本文利用水化學(xué)組分和水位聯(lián)合監(jiān)測(cè)數(shù)據(jù)分析了龍泉山地區(qū)兩口相鄰地下水監(jiān)測(cè)井的同震差異響應(yīng)特征：

（1）青白江地震在監(jiān)測(cè)井附近產(chǎn)生了4.77×10?1J/m3的能量密度，遠(yuǎn)高于地震觸發(fā)水文響應(yīng)的能量密度閾值10?4J/m3，但相鄰監(jiān)測(cè)井的水位響應(yīng)差異巨大（ZK1 的水位變幅大于ZK6），表明ZK1 的地震敏感性高于ZK6。

（2）不同含水層地下水的混合引起了水化學(xué)組分異常。深部第二含水層大量向井的補(bǔ)給，并與井水混合是導(dǎo)致ZK1 中化學(xué)組分同震異常的主要原因。ZK6水化學(xué)組分異常主要是淺部第一含水層向井補(bǔ)給，并與井水混合造成，第二含水層向井補(bǔ)給的比例較小。

（3）巖體不排水體積改變引起了地下水水位階躍變化。ZK1 水位出現(xiàn)了永久性階躍上升，震后并未恢復(fù)，原因是圍巖在地震作用下發(fā)生了塑性變形，造成第二含水層的裂隙開度減小，從而顯著降低了含水層滲透性，而第一含水層滲透性未生顯著改變。ZK6 水位出現(xiàn)了階躍上升，但震后恢復(fù)到震前水平。這表明地震造成ZK6 圍巖的彈性變形，巖體滲透性雖暫時(shí)改變，但應(yīng)力緩慢恢復(fù)后，滲透性也恢復(fù)到震前水平。