王良俊 黃仁桂 朱國偉 趙愛平 肖孟仁 李雨澤 袁 俊
1)中國江西332006 江西九江揚(yáng)子塊體東部地球動力學(xué)野外科學(xué)觀測研究站
2)中國南昌 330039 江西省地震局
氡濃度觀測是國際上公認(rèn)的地震監(jiān)測手段之一,也是我國地震地下流體觀測臺網(wǎng)中重要測項(xiàng),在地震趨勢分析與短臨震情研判中發(fā)揮著重要作用(Wang et al,2014,2018;Zhou et al,2020;Alam et al,2021;Muto et al,2021;Zhao et al,2021)。測氡儀對深層地下水(或溫泉)中溶解氣、逸出氣及斷裂帶土壤氣體中的氡氣濃度進(jìn)行連續(xù)觀測,從而捕捉地震信息(劉仕錦等,2019;劉春國等,2021)。氡被巖石、土壤顆粒表面吸附與解附,或溶解于地下水并隨地下水遷移。在地震孕育與發(fā)生過程中,巖石應(yīng)力與熱力狀態(tài)的變化及深部物質(zhì)運(yùn)移,可能會導(dǎo)致巖石氡射氣系數(shù)改變、氡溶解度變化。但不可否定,氣氡觀測中也存在許多影響觀測資料質(zhì)量的因素(趙冬等,2018;劉仕錦等,2019),如井水流量、脫氣系統(tǒng)、觀測室溫度、濕度、氣壓等環(huán)境因素,以及水井周邊環(huán)境的干擾因素等(徐長銀等,2017)。因此,利用氣氡濃度資料進(jìn)行地震分析預(yù)測時(shí),必須對異常信息開展充分的核實(shí),以判斷異常是地震前兆還是干擾所致(劉磊等,2017)。本文以九江2 井氣氡濃度異常為例,對非深部地下介質(zhì)變化氣氡濃度異常信息進(jìn)行分析,以期為氣氡濃度異常核實(shí)過程提供實(shí)例,為其他地震地下流體異常核實(shí)拓寬思路。
九江中心站位于江西廬山西北側(cè),地理位置為29.65°N、116.01°E,海拔110 m,構(gòu)造上屬于揚(yáng)子陸塊下?lián)P子地塊中部,區(qū)域構(gòu)造位于廬山西北緣邊緣剪切帶上的夏家-威家?guī)X左行走滑斷裂帶,該斷裂帶呈NE 走向。九江中心站內(nèi)的九江1 井主要用于物理量觀測,九江2 井主要用于化學(xué)量觀測(趙愛平等,2015)。九江2 井成井于2008 年,觀測含水層為構(gòu)造裂隙水,屬于承壓自流井,詳細(xì)的井孔剖面特征如圖1 所示。該井深71 m,套管11.9 m,地表至6.2 m 為第四系聯(lián)圩組亞黏土、砂礫石,6.2—11.9 m 為震旦系皮園村組強(qiáng)風(fēng)化碳質(zhì)灰?guī)r,11.9—71.0 m 均為下元古代碳質(zhì)灰?guī)r,含水層分別在 17.0—22.5 m、55.0—63.5 m,目前日流量約為300 t(圖1)。距該井約1 km處有天花井水庫,50 m 處有1 條溪流,溪流源頭為廬山馬尾水泉水。
圖1 九江2 號觀測井井孔柱狀圖Fig.1 Column diagram of hole of Jiujiang No.2 observation well
近年來,九江2 井氣氡濃度3 次出現(xiàn)異常,出現(xiàn)異常后均對脫氣裝置、氣路等進(jìn)行排查處理,3 次排查均為氡濃度上升異常:①2014 年7 月22 日氣氡濃度異常。九江2 井氣氡濃度從7 月21 日11 時(shí)開始增大,之前背景值為30—40 Bq/L,15 時(shí)增大到170 Bq/L 以上,之后持續(xù)上升,7 月22 日16 時(shí)達(dá)512.4 Bq/L,7 月27 日氣氡濃度回至40 Bq/L 左右。②2015 年6 月24 日氣氡濃度異常。九江2 井氣氡濃度從6 月22 日11 時(shí)開始增大,23 日由均值360 Bq/L 升高到480 Bq/L。③2021 年4 月5 日氣氡異常。4 月5 日起九江2 井氣氡濃度持續(xù)升高,日變化幅度約為10 Bq/L,6 日內(nèi)累積上升84 Bq/L,氣氡濃度上升幅度達(dá)43%。
2014年7月21日九江2井氣氡濃度異常如圖2所示。九江中心站各觀測儀器運(yùn)行正常,觀測環(huán)境未發(fā)生變化,對同測點(diǎn)其他流體測項(xiàng)和同臺形變測項(xiàng)等資料進(jìn)行綜合分析,未發(fā)現(xiàn)與氣氡濃度同步異常測項(xiàng),基本可排除是構(gòu)造運(yùn)動引起的氣氡濃度異常,即異常屬干擾類型。由圖2 可見,2014 年7 月15—16 日,九江中心站區(qū)域出現(xiàn)強(qiáng)降雨,降雨量達(dá)78 mm,7 月24 日再次降雨52 mm,九江靜水位明顯上升,達(dá)80—100 mm。九江井氣氡脫氣裝置為傳統(tǒng)的濺落式脫氣裝置,容積約1.25 L,被置放于下沉式池內(nèi),靠近主井管,極易受井管噴溢出水的影響。由于區(qū)域強(qiáng)降雨使水位上升,井管出水口泄流速率增大,增強(qiáng)了井管溢出水噴濺能力,井管溢出水強(qiáng)烈的噴濺打破了SD-3A 型測氡儀脫氣裝置的平穩(wěn)狀態(tài),改變了脫氣裝置內(nèi)部水—?dú)猸h(huán)境的平衡(高小其等,2021)。受此影響,7 月21 日起,氡氣濃度出現(xiàn)大幅異常變化,異常持續(xù)4 天左右后開始下降,異常期間水位處于較高水平,即水位較大幅度的變化對氣氡濃度有影響。此次異常主要影響機(jī)制為水位大幅升高導(dǎo)致主井管溢流速率增大,溢出井水直接濺落到脫氣裝置,從而影響脫氣裝置內(nèi)水—?dú)馄胶?,進(jìn)而影響氣氡濃度變化。為改變脫氣裝置內(nèi)水—?dú)馄胶猓?015 年1 月1 日九江中心站采用自然吸氣脫氣裝置進(jìn)行鼓泡脫出氡氣。
圖2 2014 年九江降雨量(a)、靜水位(b)、2 井氣氡濃度異常數(shù)據(jù)(c)Fig.2 Fluid observation data of Jiujiang in 2014 (a) rainfall,(b) static water level,and (c) abnormal radon concentration of well 2
2015 年6 月24 日九江2 井氣氡濃度異常如圖3 所示。在異常期間,中心站觀測儀器、觀測環(huán)境正常,流體、形變、測震等觀測資料均未出現(xiàn)同步變化,即此次異常為非深部地下介質(zhì)變化引起的氣氡濃度異常,屬于干擾類型。由圖3 可見,2015 年6 月20—21 日九江地區(qū)發(fā)生強(qiáng)降雨,降雨量達(dá)74 mm,水位明顯上升,井水噴出井口,氣氡濃度升高滯后于水位上升,水位上升滯后于降雨。分析認(rèn)為,主要影響機(jī)制如下:①大量降雨增加含水層荷載,擠壓含水層孔隙、裂隙,使其中地下深處含水層部分氡組分隨水進(jìn)入井孔;②水位上升增加井孔水壓,進(jìn)而增加水中氡氣的溶解度;③井口噴水增大出水流速,含水層內(nèi)水流速度在一定程度上也增大,水流沖刷含水層能力增強(qiáng),將含水層中更多的氡成分帶入井孔;④井水自井口濺落,觀測室內(nèi)混合氣體中的氡氣濃度升高,自然吸氣脫氣裝置吸收該氣體鼓泡,使氡氣濃度測值進(jìn)一步增大。九江中心站氣氡濃度此次在原背景值水平下升高幅度達(dá)30%,其主要干擾因素為強(qiáng)降雨直接導(dǎo)致氣氡濃度的升高。
圖3 2015 年九江降雨量(a)、靜水位(b)、2 井氣氡濃度異常數(shù)據(jù)(c)Fig.3 Fluid observation data of Jiujiang in 2015 (a) rainfall,(b) static water level,(c) abnormal radon concentration of well 2
2021 年4 月5 日九江2 井氣氡濃度異常如圖4 所示。對九江中心站九江1 井、中心站周邊環(huán)境進(jìn)行調(diào)查,未發(fā)現(xiàn)較大的施工工程開工、打井、抽水等明顯干擾源。此次氣氡濃度異??赡苌婕暗降牡叵滤瘜W(xué)組分、水位埋深、地傾斜、應(yīng)變等均沒有明顯變化。由圖4 可見,九江2 井氣氡濃度與觀測室室內(nèi)溫度、室外溫度之間存在緩慢同步上升現(xiàn)象,特別是4 月28 至5 月10 日室內(nèi)溫度與氣氡濃度出現(xiàn)間斷性同步突跳變化現(xiàn)象。分析認(rèn)為,氣氡濃度變化的另一主要影響因素為觀測系統(tǒng)的影響。圖5 為九江2 井觀測系統(tǒng)供水管道圖及示意圖。由于氡氣觀測儀器檢測平臺管道(3 路)與氣氡供水管道從主井觀測室出發(fā),沿室外地埋溝(管槽)一直進(jìn)入到氡氣觀測儀器檢測平臺室內(nèi)(觀測房),4 路管道相鄰安置。此次氣氡濃度異常之前6 個(gè)月,檢測平臺一直處于滿管實(shí)驗(yàn)狀態(tài),2021 年3 月下旬起停止實(shí)驗(yàn),把檢測平臺水路管道和恒流水箱水排空,4 月初氣氡濃度出現(xiàn)異常。隨著夏季到來,氣溫不斷升高,管道出現(xiàn)熱脹冷縮,流水管道內(nèi)形成部分間隙,同時(shí)氡氣濃度隨溫度升高會產(chǎn)生更多的逸出氣體,逸出氣氡混合到地下水的量緩慢達(dá)到最高,而后緩慢降低,恢復(fù)到正常的狀態(tài),即溫度升高,水路管道擠壓導(dǎo)致壓力變化,氣氡濃度緩慢變化。
圖4 2021 年九江室外溫度(a)、室內(nèi)溫度(b)、2 井氣氡濃度異常數(shù)據(jù)(c)Fig.4 Fluid observation data of Jiujiang in 2021 (a) outdoor temperature,(b) indoor temperature,(c) abnormal radon concentration of well 2
圖5 九江2 井觀測系統(tǒng)供水管道Fig.5 Water supply pipeline diagram of Jiujiang Well 2 observation system
表1 是九江中心站7 次降雨、水位變化、氣氡濃度異常統(tǒng)計(jì)表。在一定范圍內(nèi)九江中心站氣氡濃度與水流量呈正相關(guān),氣氡濃度高值異常由水位變化、降雨量所引起。如2015 年6 月21—22 日(B 時(shí)段)氣氡濃度異常與降雨量、水位上升間為直接相關(guān),氣氡濃度變化率為43.68%;2021 年4 月1—14 日(G 時(shí)段)降雨量75 mm,水位變化0.089 m,氡值變化率為43.68%,與2015 年6 月21—22 日(B 時(shí)段)類似,從九江中心站氣氡值變化率、水位和降雨對應(yīng)圖3 來看,相關(guān)性較好。但從A、C、D、E、F 時(shí)間段看,這些時(shí)段的降雨量和水位變化均大于B、G 時(shí)間段,而氣氡濃度變化率較低。這說明,氣氡濃度異常變化與九江中心站降雨量、水位變化間不是簡單的線性關(guān)系,需要分別作具體分析,才能搞清楚三者間的關(guān)聯(lián)性。
表1 九江中心站總降雨量、水位變化、氣氡濃度異常Table 1 The rainfall,change of the water level and radon concentration anomaly in Jiujiang Central Station
2021年4月1—14日(G深度)氣氡濃度異常期間,采集九江中心站3個(gè)位置(1—井口;2—氡檢測平臺;3—站內(nèi)廢井)的3 次水樣送至中國地震局地震預(yù)測研究所進(jìn)行水化組分測試,測試結(jié)果如圖6 所示。由圖6 可見,Piper 圖顯示常量離子的含量較接近,說明補(bǔ)給源一致。Schoeller 圖清晰地顯示主要離子成分的含量基本接近,說明是同一補(bǔ)給來源(張磊等,2016,2019),與Piper 圖的結(jié)果一致。同時(shí)可見,除Cl-離子有微小變化外,其他基本沒有變化,這進(jìn)一步證實(shí)異常是由觀測系統(tǒng)所致。
在實(shí)際觀測中,氣氡濃度受水位、降雨量、溫度、水溫、壓力、礦化度等多種因素的影響,在核實(shí)氣氡濃度3 次異常變化時(shí)發(fā)現(xiàn),九江中心站其他地球物理測項(xiàng)均無異常,可排除相關(guān)影響。研究結(jié)論如下:①水位上升,脫氣裝置進(jìn)水流量增大,單位時(shí)間內(nèi)吸入的空氣也隨之增加,進(jìn)而提升水中溶解氡的脫氣效率,同時(shí),水位上升導(dǎo)致的主井管地下水溢流可直接影響脫氣裝置的水—?dú)馄胶?,也?dǎo)致氣氡濃度上升;②大量降雨增加九江2 井含水層荷載,擠壓含水層孔隙、裂隙,含水層內(nèi)水流速度在一定程度上也增大,水流沖刷含水層能力增強(qiáng),將含水層中更多的氣氡帶入井孔,多重因素的共同作用導(dǎo)致氣氡濃度增大;③氣氡觀測系統(tǒng)受溫度的影響,水路管道熱脹冷縮,水路擠壓形成間隙,導(dǎo)致脫氣壓力緩慢變化,氣氡濃度隨之出現(xiàn)緩慢上升。
為進(jìn)一步分析九江2 井氣氡濃度3 次出現(xiàn)異常的原因,后續(xù)將開展加密實(shí)驗(yàn),建立完整的水位、水流量與氣氡濃度整體觀測系統(tǒng),分析氣氡濃度與水位、降雨量、水流量間的關(guān)系及其變化特征,為后期對氣氡濃度觀測資料進(jìn)行分析、干擾排查、脫氣裝置改造等提供參考。必要時(shí)還應(yīng)充分利用異常測項(xiàng)本身的物理、化學(xué)性質(zhì)及水井的基礎(chǔ)資料進(jìn)行綜合判定,以期為氣氡濃度異常核實(shí)跟蹤提供實(shí)驗(yàn)支撐,從而提高九江中心站氣氡濃度觀測效能。