摘要： 活動斷裂帶氣體地球化學(xué)空間分布特征分析是區(qū)域構(gòu)造和斷裂活動性分區(qū)（分段）及其地震危險性分析的重要方法之一。通過臨潭—宕昌斷裂帶9條測線剖面的土壤氣Hg、Rn濃度觀測，結(jié)合斷裂帶分段活動特性與大地?zé)崃鞅尘皵?shù)據(jù)，開展不同分段氣體濃度空間分布特征分析，深入探討土壤氣地球化學(xué)空間分布的特征及成因。結(jié)果表明：斷層土壤氣在臨潭—宕昌斷裂帶具有明顯的分段特征，根據(jù)其活動特性可將斷裂帶分為東、西兩段，且斷層土壤氣地球化學(xué)空間分布與大地?zé)崃骱偷卣鸹顒又g存在良好的相關(guān)關(guān)系;在斷裂西段高熱流區(qū)域，中強地震頻發(fā)，但小震活動較弱，Hg、Rn濃度強度較低;東段低熱流區(qū)域，中強地震較少，但小震活動活躍，Hg、Rn濃度強度較高。研究結(jié)果作為區(qū)域地球化學(xué)背景特征，可為今后地下流體異?；顒臃治鎏峁┲匾罁?jù)，對進一步跟蹤未來可能發(fā)生地震的斷層及其活動段具有重要意義。

關(guān)鍵詞： 斷層土壤氣; 氣體地球化學(xué); 臨潭—宕昌斷裂帶; 地震活動

中圖分類號： P315.724 文獻標(biāo)志碼：A 文章編號： 1000-0844（2024）06-1462-13

DOI：10.20000/j.1000-0844.20240113001

Geochemical spatial distribution characteristics of

soil gas along the Lintan-Tanchang fault zone

LI Ruigang¹， LI Chenhua^{1，2，3}， SU Hejun^{1，2，3}， ZHOU Huiling^{1，2，3}， WANG Yanhong³， WAN Yue¹

（1. Lanzhou Institute of Seismology， CEA， Lanzhou 730000， Gansu， China;

2. Gansu Lanzhou Geophysics National Observation and Research Station， Lanzhou 730000， Gansu， China;

3. Gansu Earthquake Agency， Lanzhou 730000， Gansu， China）

Abstract： The study on geochemical spatial distribution characteristics of gas along the fault zone is an important method for analyzing the regional structure，activity segmentation，and seismic risk of the fault.In this study，we took the Lintan—Tanchang fault zone as study area and observed the soil gas Hg and Rn concentrations in 9 survey lines.We analyzed the spatial distribution characteristics of gas concentration in different fault sections by combining the characteristics of fault activity and the heat flow density.The study show results that the Lintan—Tanchang fault zone has obvious segmentation features.Base on the fault activity characteristics，the Lintan—Tanchang fault could be divided into two segments：east and west.And the spatial distribution of soil gas geochemistry is correlated with heat flow and seismic activity;The western section of the fault is an area of high heat flow，and strong seismic occur frequently in this section，but small seismic activity is weak，with lower Hg and Rn concentration intensity.While in areas with lower heat flow in eastern segment of fault，there are fewer strong seismic but intense small seismic activity，with higher Hg and Rn concentrations intensity.This study provides evidence of subsurface fluid activity in the study area，which is vital in tracing the faults and active segments with seismic risk.

Keywords： fault soil gas;gas geochemistry;Lintan-Tanchang fault;seismic activity

0 引言

地下流體是較為活躍的地球物質(zhì)之一，深大斷裂是地球深部熱能和流體釋放的主要通道之一^［1］。氣體作為一種移動性和壓縮性較強的流體，可沿活動斷裂帶逸出，導(dǎo)致斷裂帶附近土壤氣地球化學(xué)濃度異常^［2-4］。國內(nèi)外學(xué)者研究發(fā)現(xiàn)Hg、H₂、Rn與CO₂等斷層土壤氣在地震發(fā)生前后常出現(xiàn)較為明顯的濃度異常^［5-8］。深入研究發(fā)現(xiàn)這些異常特征與地震強度有關(guān)，通常地震震級越大土壤氣異常幅度越大，震中距越近氣體異常濃度越高，這為斷裂活動性的空間分布特征提供了一種可行的研究方向^［9-12］。20世紀(jì)70年代，Irwin等^［13］通過對北美地區(qū)地震活動構(gòu)造帶附近的泉水進行分析，發(fā)現(xiàn)地震活動性強的區(qū)域存在CO₂濃度異常高值區(qū)，監(jiān)測該區(qū)域CO₂排放濃度變化有助于分析斷裂帶的活動性及地震風(fēng)險。2017年Fu等^［14］在臺灣省北部地區(qū)的活動斷裂帶建立了多個土壤氣體自動監(jiān)測站，經(jīng)歷一年的監(jiān)測后發(fā)現(xiàn)，多數(shù)地震發(fā)生前幾個小時甚至幾天都存在高Rn濃度的異?，F(xiàn)象。近年來，隨著氣體監(jiān)測技術(shù)逐漸成熟，對活動斷裂帶土壤氣體的分析研究成為分析地震危險性的有效手段之一^［15］。

臨潭—宕昌斷裂帶結(jié)構(gòu)復(fù)雜，分支較多，斷裂全長約250 km，總體可分為西段與東段兩個區(qū)域。該斷裂帶經(jīng)過臨潭、岷縣至宕昌，是西秦嶺NWW向斷裂系中的一個重要組成部分^［16］。因其復(fù)雜的構(gòu)造地貌與強烈的活動性，斷裂帶上強震頻發(fā)，如1573年岷縣63/4級、1837年岷縣6級、2003年岷縣5.2級、2013年岷縣—漳縣6.6級、2019年夏河5.7級地震等^［17-19］。這些地震給當(dāng)?shù)厝嗣裆敭a(chǎn)安全造成了嚴(yán)重危害，因此，對該區(qū)域未來地震危險性的判定引起了越來越多的關(guān)注。

近年來甘東南地區(qū)地震頻發(fā)，但過去的研究主要集中在深大斷裂及塊體邊界控制性斷裂的氣體分段特征分析上，對于塊體內(nèi)部臨潭—宕昌斷裂帶氣體分布特征的研究相對較少。臨潭—宕昌斷裂位于甘東南復(fù)雜的構(gòu)造轉(zhuǎn)換區(qū)內(nèi)，其氣體濃度是否存在與相鄰邊界斷層相同的明顯分段特性需要深入研究。因此，本文通過跨斷層測量獲取臨潭—宕昌斷裂土壤氣濃度數(shù)據(jù)，開展Hg與Rn濃度強度及其空間分布特征分析，追蹤斷裂帶不同活動段跨斷層土壤氣地球化學(xué)特征，以期為分析臨潭—宕昌斷裂帶地下流體孕震環(huán)境提供基礎(chǔ)科學(xué)依據(jù)。

1 研究區(qū)地質(zhì)背景

臨潭—宕昌斷裂帶位于青藏高原東北緣東昆侖斷裂帶與西秦嶺北緣斷裂帶過渡地區(qū)^［20］，參與了西秦嶺的形變分配，繼承了東昆侖斷裂帶與西秦嶺北緣斷裂帶的左旋走滑特性。其受到青藏高原東北緣的擠壓以及鄂爾多斯塊體的阻擋作用，產(chǎn)生了明顯的逆沖變形的運動特征^［21］。斷裂帶總體呈NWW—NW向展布，傾向NE，傾角50°～70°，地貌主要呈現(xiàn)為U型谷，部分斷層呈陡坎形態(tài)，形成向南彎曲的弧形構(gòu)造^［22］。臨潭—宕昌斷裂帶由西端的合作開始至岷縣，發(fā)展出多條構(gòu)造相似但規(guī)模不同、錯落分布的次級斷裂，部分次級斷裂較為復(fù)雜，發(fā)育有一條或多條分支。根據(jù)近年的研究結(jié)果，臨潭—宕昌斷裂帶是一條晚第四紀(jì)活動斷裂，活動性質(zhì)自全新世早期以來發(fā)生過改變，是西秦嶺內(nèi)部活動性較強的斷裂。前人研究結(jié)果表明該斷裂帶左旋滑動速率為0.86～1.6 mm/a，垂直滑動速率為0.05～0.1 mm/a，左旋走滑分量大于傾滑分量，總體呈左旋走滑性質(zhì)，個別區(qū)域表現(xiàn)為向南逆沖^［23-24］。斷裂西段曾錯斷地表，區(qū)域構(gòu)造應(yīng)力逐漸積累，為孕震提供了適宜的環(huán)境^［25］。

根據(jù)中國地震臺網(wǎng)地震目錄記錄分析，臨潭—宕昌斷裂帶自全新世以來一直處于高度活躍狀態(tài)，構(gòu)造活動強烈，地震頻發(fā)，2003—2019年就發(fā)生過6次M_S5.0～7.0的中強震及多次不同規(guī)模的小震^［26］。中強震的震源位置與斷裂帶復(fù)雜的構(gòu)造體系密切相關(guān)（圖1）。

2 采樣與分析方法

本研究共布設(shè)了9條測線，每條測線均垂直或大角度相交于臨潭—宕昌斷裂帶主斷裂。其中，斷裂帶西段布設(shè)4條測線，東段布設(shè)5條測線，由西向東分別為合作、貢恰、新城鎮(zhèn)、柰子溝、夭馱、莊子村、麻屆灘、巴沙溝、毛羽村（圖1）。采集不同測線上的斷層土壤氣并進行濃度測量工作，記錄氣體Hg、Rn相關(guān)數(shù)據(jù)。由于氣體在近地表會受到風(fēng)速、風(fēng)向等外界條件干擾，而地下1 m左右深度土壤氣濃度較為穩(wěn)定，故采樣深度選取為地表下80 cm，共監(jiān)測6～12組數(shù)據(jù)。采樣時間為2022年6月每日9：00—18：00，為避免實驗數(shù)據(jù)受到氣象因素影響，選擇天氣狀況穩(wěn)定、良好的時間段進行野外采樣。測量地點選擇植被覆蓋少，土壤下具有原始覆蓋層的區(qū)域。測量過程中時刻監(jiān)視儀器狀態(tài)，保證儀器正常使用且未受到其他因素干擾，每日對所有采樣儀器進行維護以減少誤差，確保數(shù)據(jù)真實準(zhǔn)確。

2.1 Hg測量方法

氣體Hg測量使用加拿大LUMEX公司生產(chǎn)的RA-915M型號便攜式測汞儀。儀器測量范圍：0.5～200 ng/L;測量相對誤差不超過±5%。測線垂直或大角度相交于主斷裂帶，每隔10 m布置一個測點，共計6～12組測點。測量方法采用打釬測量：首先，在選定區(qū)域以錘擊鋼釬鑿出一個深度80 cm的孔洞;然后，將空心取樣器插入其中并填埋周圍土壤，避免地面空氣稀釋土壤氣體;最后，使用橡膠管連接取樣器與儀器，進行測量（圖 2）。儀器使用前需預(yù)熱2～10 min，數(shù)據(jù)無法保存，需人工選取每個測點最大值并記錄。其中，巴沙溝測線測得的Hg濃度低于儀器檢出限，無法獲得準(zhǔn)確數(shù)據(jù)，故下文對該測線的Hg氣體不進行分析。

2.2 Rn測量方法

氣體Rn測量使用德國Genitron公司生產(chǎn)的AlphaGUARD PQ2000便攜式測氡儀。儀器的操作范圍為：2～200 000 Bq/L;溫度：-10～+50 ℃;氣壓：70～110 Pa;相對濕度：0%～95%。測量地點選在Hg測點附近，打孔完畢后使用橡膠管與儀器連接測量，每個測點測量得到10～15組數(shù)據(jù)，最后根據(jù)規(guī)范作數(shù)據(jù)匯總及處理工作。

3 數(shù)據(jù)處理

由于不同斷層或斷層的不同段內(nèi)部破碎程度不同，氣體滲透率不同，導(dǎo)致土壤氣向地表的逸出強度不同，并在空間上表現(xiàn)出明顯的差異性。這種差異性可能會受到土壤結(jié)構(gòu)、測量環(huán)境、氣象條件和人為因素的干擾。為減小場地因素對斷層氣特征的影響，采用單剖面均方差與平均值的計算方法，繪制斷層單剖面土壤氣濃度分布圖（圖3、圖4）。

經(jīng)計算與評估，使用最大值法計算土壤氣濃度強度，使用異常下限確定濃度異常值的量值界限。

最大值法是選取一條測線上所有測點的最大濃度值D_max，將該測線的氣體平均濃度作為數(shù)據(jù)處理的背景值K，二者的比值即為該測點的土壤氣濃度強度S：

S = D_max/K （1）

異常下限是判斷氣體濃度曲線形態(tài)的重要依據(jù)。利用土壤氣背景值K與測線濃度的均方差δ，即可計算出各測線的異常下限X_a［式（2）］。其中異常值可以用背景值加上“n”個標(biāo)準(zhǔn)差（n= 1，2，3，…）來表示，本文中n的取值為1。

X_a=K+δ （2）

上述方法可以較大幅度地減小氣體逸出的環(huán)境因素影響，更直觀地表現(xiàn)各剖面濃度強度的峰值。因此，本文利用此方法，計算得出臨潭—宕昌斷裂帶9條測線土壤氣中Hg、Rn的最大值D_max、背景值K、濃度強度S和異常下限X_a（表1）。

4 分析與討論

4.1 土壤氣空間分布特征與臨潭—宕昌斷裂帶活動性的關(guān)系

臨潭—宕昌斷裂帶是東昆侖斷裂系統(tǒng)和西秦嶺斷裂系統(tǒng)之間一個重要的轉(zhuǎn)換紐帶，其后期構(gòu)造變形受到西秦嶺北緣斷裂向兩側(cè)的擴展與東昆侖斷裂向東北方向的推擠，因此斷裂周邊形成了復(fù)雜的應(yīng)力構(gòu)造環(huán)境。根據(jù)全球定位系統(tǒng)速率場變化情況可知，該斷裂帶呈NWW—NW向展布，且局部區(qū)域發(fā)生構(gòu)造變形，斷裂帶上中強地震的孕育也是受到這種構(gòu)造特征引起的斷裂帶附近的應(yīng)力集中所致。斷裂帶自第四紀(jì)以來活動性質(zhì)以向南逆沖為主，局部陡立或南傾，受到應(yīng)力影響發(fā)生順時針旋轉(zhuǎn)，具有左旋走滑分量。斷裂帶自岷縣向西北方向發(fā)育數(shù)條規(guī)模不同、相互平行或斜接的次級斷裂，其影響寬度范圍為5～10 km。岷縣東南斷裂帶分布較為均一，以主斷裂為主體，至宕昌以南，不同斷裂段上的活動性差異較大。西秦嶺北緣斷裂走滑及向南北兩側(cè)逆沖的“花狀構(gòu)造”是臨潭—宕昌斷裂帶上中強地震頻發(fā)的重要動力因素。靠近西秦嶺北緣斷裂區(qū)域的構(gòu)造應(yīng)力容易集中釋放，并且斷裂帶東段頻繁發(fā)生地震可能會引起西段的應(yīng)力調(diào)整，從而導(dǎo)致中強震的發(fā)生^［27-28］。

根據(jù)斷裂幾何形態(tài)展布特征及其所處的大地?zé)崃鞣植继卣?，以岷縣為界將臨潭—宕昌斷裂帶分為東、西兩段：西段包括合作、貢恰、新城鎮(zhèn)與柰子溝，東段包括夭馱、莊子村、麻屆灘、巴沙溝與毛羽村。由表1可知，臨潭—宕昌斷裂帶Hg與Rn濃度強度范圍分別為1.327～2.878與1.542～8.880;Hg濃度強度高值位于毛羽村測線，Rn濃度強度高值位于麻屆灘測線，均分布于斷裂東段。兩種氣體濃度強度的高值點在斷裂帶上表現(xiàn)出一致的空間分布特征，氣體逸出較為同步，受相似的斷層活動性控制。氣體逸出通道除受斷層活動性影響外，還與各類氣體自身的密度和固有特性有關(guān)^［29-30］。

由圖3、4可知：Hg、Rn土壤氣濃度曲線整體呈現(xiàn)單峰異常狀態(tài)，異常峰值集中分布在斷層周圍10～30 m范圍內(nèi)，推斷所測量斷裂帶區(qū)域近年來活動性較強，地下裂隙連通性穩(wěn)定，為土壤氣向地表逸出提供了較通暢的環(huán)境；少部分測線形態(tài)呈現(xiàn)多峰聚集模式，以斷層為中心的20～50 m內(nèi)Hg與Rn濃度出現(xiàn)多個峰值，在其余范圍內(nèi)濃度變化趨于平緩并形成一個“低谷”。由臨潭—宕昌斷裂帶Hg與Rn濃度強度空間分布特征可知（圖5），東段多處測線濃度強度略高于西段，推測東段可能具有較強的活動性，應(yīng)力積累周期較短且快速釋放，通常表現(xiàn)為小震的形式，導(dǎo)致此區(qū)域常年小震頻發(fā)。在多年地震的影響作用下，斷層核部裂隙閉鎖程度低，裂隙發(fā)育程度高，氣體上升通道穩(wěn)定。由于Hg與Rn氣體密度相對較大，自身遷移速率低，擴散作用不明顯，往往通過載氣作用從地殼深部遷移至地表^［31］。東段長期的斷層滑移導(dǎo)致斷層間相互摩擦作用增強，加大了地下深部氣體擴散率，開辟出更多的裂隙通道，為土壤氣的逸出提供了合適的環(huán)境，經(jīng)過長年不斷的積累，導(dǎo)致東段多處測線氣體濃度異常。斷裂帶西段地震活動相對較少，但M_S≥5.0地震均發(fā)生在此區(qū)域，推測此段構(gòu)造應(yīng)力經(jīng)過長期積累后釋放，但強震離逝時間久，斷層滑動速率低，斷裂逐漸封閉，斷層核部相對閉鎖，氣體向上逸出的通道狹隘，影響土壤氣的滲透性，使得西段部分測點土壤氣濃度強度略低。

4.2 地震活動與大地?zé)崃骺臻g的關(guān)系

大地?zé)崃魇侵傅乇硐虏吭趩挝粫r間內(nèi)以熱傳導(dǎo)方式垂直向上，流經(jīng)單位面積的熱流量，也稱為熱流或大地?zé)崃髅芏?。大地?zé)崃鞯姆植寂c地殼運動、地球內(nèi)部構(gòu)造等地質(zhì)條件之間有顯著的聯(lián)系，熱流既可以為人類提供多種不同用途的地?zé)崮芰抠Y源，也為人類探索地球構(gòu)造提供了重要參考依據(jù)^［32-34］。早在1997年，就有研究者曾在歐洲發(fā)現(xiàn)地表地?zé)崤c地震發(fā)生有顯著的關(guān)系^［35］。我國研究者也提出探討地?zé)釋α私獾卣鸬脑杏兄匾饬x，并通過研究發(fā)現(xiàn)異常地?zé)釁^(qū)域與地震活動在地質(zhì)構(gòu)造上有著一定的吻合^［36］。截至2016年，我國已采集并發(fā)表了1 230個大地?zé)崃鲾?shù)據(jù)^［37］。

受到大地?zé)崃饔绊懙奈镔|(zhì)種類繁雜，相比于長時間與大氣接觸交換的地表氣體，部分地?zé)釟怏w直接來源于地球內(nèi)部，不易受到陸地環(huán)境的影響，可以傳遞地球深部的地質(zhì)環(huán)境信息，使研究人員得以更準(zhǔn)確地分析地下成分^［38-39］。其中，Rn氣體通常來源于地球深部，可溶于水，隨著水流遷移至地表或從地下破碎帶、發(fā)育的裂隙等通道逸散至地表。地?zé)崴鶄鬟f的能量會影響Rn在水中的溶解度，地下溫度越高，Rn在水中的溶解度越低，熱水蒸氣中的Rn含量也越高。這些Rn熱氣最終會通過裂隙通道遷移至地表，在地?zé)釡囟雀摺嗔褞Щ钴S且地下破碎程度高的區(qū)域表現(xiàn)出Rn氣體濃度異常。Hg氣體主要以氣態(tài)方式通過破碎通道從地殼深部逸散至地表，其自然來源有地?zé)峄顒?、火山噴發(fā)、地下含Hg礦物成礦后Hg蒸氣的揮發(fā)等，故Hg氣體也受到地?zé)崤c區(qū)域地質(zhì)構(gòu)造影響^［40］。由此可見，在斷裂帶不同地?zé)釁^(qū)域收集Rn、Hg氣體濃度有助于分析斷裂帶的活動性質(zhì)。

根據(jù)數(shù)據(jù)測量結(jié)果，中國大陸區(qū)域大地?zé)崃髌骄禐?0.4 mW/m²，與全球大陸平均值65 mW/m²接近^［36］。甘肅省位于中國大陸西北地區(qū)，地形呈狹長狀，地貌復(fù)雜多樣，西部以沙漠、戈壁為主要地貌，平均熱流小于60 mW/m²，東部以黃土高原為主，平均熱流大于60 mW/m²，熱流分布總體呈現(xiàn)西部低于東部。結(jié)合中國大陸區(qū)域大地?zé)崃?，可將高于平均熱流的區(qū)域歸為高熱流區(qū)域，將低于平均熱流的區(qū)域歸為低熱流區(qū)域。臨潭—宕昌斷裂帶處于甘東南地區(qū)，近年來地震活動頻發(fā)。如圖6所示，該斷裂帶橫跨多個不同層次的地?zé)釁^(qū)域，總體地?zé)岱秶鸀?8～76 mW/m²。根據(jù)袁伏全等^［41］對青海不同地?zé)釁^(qū)域地震分布的分析，強震分布與斷裂構(gòu)造、大地?zé)崃鞣植加忻芮嘘P(guān)聯(lián)：1970—2016年間，青海共發(fā)生M_S≥6.0強震35次，其中92%的強震發(fā)生在大地?zé)崃鞣秶?0～80 mW/m²的地區(qū);震中位置通常不在地?zé)嶂档闹行奈恢?，往往位于發(fā)生地?zé)嶙兓臋M向梯度帶上。2013年岷縣—漳縣M_S6.6地震前臨潭—宕昌斷裂帶地震b值總體為0.53～1.65^［42］，但存在兩處明顯的低b值（b＜0.6）區(qū)域，即夏河縣與臨潭、岷縣之間的地區(qū)，震中位置位于地?zé)岬臋M向梯度變化帶上。此外，對1990—2010年地震b值進行掃描，發(fā)現(xiàn)岷縣附近地區(qū)b值異常顯著，2013年岷縣—漳縣6.6級地震前岷縣西北方向鄰區(qū)b值明顯降低^［42］。根據(jù)b值大小與強震發(fā)生的可能性關(guān)系來看，低b值異常范圍越大，其強震潛在危險性就越大^［43］，岷縣周圍多次強震恰好印證了此結(jié)論。本文通過臨潭—宕昌斷裂帶所處大地?zé)崃鞣植挤秶cb值異常區(qū)域，將斷裂帶分為東、西兩段：西段位于高地?zé)釁^(qū)域，熱流大于60 mW/m²;東段主要位于低熱流區(qū)域，熱流小于60 mW/m²。1970年以來，該斷裂帶上的M_S≥5.0地震均發(fā)生于熱流范圍61～76 mW/m²的區(qū)域，如高熱流的橫向變化梯度上分別于2003年和2013年發(fā)生了岷縣M_S5.5和岷縣—漳縣M_S6.6地震，高熱流區(qū)域2019年發(fā)生了夏河M_S5.7地震。斷裂自西向東熱流值逐漸遞減，斷裂東段1970—2023年未發(fā)生過M_S≥5.0地震，但頻繁發(fā)生小震。

前人對地?zé)釄雠c地震分布進行研究，發(fā)現(xiàn)我國多數(shù)地震的能量釋放區(qū)域與大地?zé)崃鞯母邿崃鲄^(qū)域相吻合^［44-45］。對1976年唐山M_S7.8地震的研究^［46］表明強震主要發(fā)生于地?zé)釞M向變化較大的梯度帶上，這類區(qū)域具有最大水平剪應(yīng)力與周向應(yīng)力，兩種應(yīng)力共同促進了地震的發(fā)生，而地震所引發(fā)的地殼活動會影響地下破碎帶的破碎程度，進而影響到土壤氣濃度分布。將圖6與表1中所測得的氣體濃度結(jié)合來看，斷裂帶西段為高熱流區(qū)域，其中接近熱流橫向梯度帶變化的區(qū)域發(fā)生過多次強震，應(yīng)力積累周期長。雖然西段高熱流會促進氣體向上遷移，但強震孕育周期長，對破碎帶的影響有限，斷層核部相對閉鎖，Hg、Rn逸出通道受阻，導(dǎo)致地表氣體濃度強度較低。斷裂帶東段為低熱流區(qū)域，內(nèi)部能量穩(wěn)定，地震受熱流的影響較小，應(yīng)力積累周期短，但長期以來小震活躍，導(dǎo)致地下破碎程度較高，氣體逸出通道暢通，所以Hg、Rn濃度強度較高。

此外，地震震中與臨潭—宕昌斷裂帶大地?zé)崃鞣植加兄芮械穆?lián)系。高熱流區(qū)域與熱流發(fā)生梯度變化的區(qū)域發(fā)生強震的風(fēng)險較高，而強震的發(fā)生會導(dǎo)致較強的應(yīng)力釋放，強震之后會有長期的應(yīng)力積累時間，并且會略微降低小震發(fā)生頻次。與小震活躍區(qū)域相比，該區(qū)域?qū)Φ叵缕扑閹У挠绊懴鄬^小，地下氣體向上傳遞能力較弱，氣體濃度強度較低。長期小震活躍的斷裂帶區(qū)域地下裂隙發(fā)育程度高，地下氣體逸出通道穩(wěn)定，并且更容易長期積累氣體，導(dǎo)致土壤氣濃度略高。

4.3 臨潭—宕昌斷裂帶地球化學(xué)特性與地震活動性的關(guān)系

土壤中孔隙分布面積廣，地下氣體在壓力、溫度、濃度差異的驅(qū)動下攜帶地球深部的相關(guān)信息，穿過薄弱的破碎帶，在土壤中擴散、對流，最終釋放至大氣中。在某種程度上，土壤氣濃度強度反映了斷裂帶的構(gòu)造活動性質(zhì)^［47］，而斷裂帶的滑動速率反映了斷層應(yīng)變的積累速率，代表了斷裂長期平均活動水平^［48］。不均勻的滑動速率分布會改變斷裂的走向，影響斷層的形變，而具有一定傾斜角度的斷層有利于能量的積累，易引發(fā)不同強度的地震。結(jié)合土壤氣濃度強度分布特征與斷層滑動速率進行分析，有助于了解不同斷裂帶區(qū)域的地震危險性^［49-50］。前人的研究主要以滑動速率對斷裂造成的構(gòu)造形變作為了解地震特征的主要依據(jù)，現(xiàn)對比大地?zé)崃鞯姆侄翁卣骱蟀l(fā)現(xiàn)，將斷裂的形變特征與大地?zé)崃鲾?shù)據(jù)結(jié)合可以更有效地判斷地震的活動特性，分析不同區(qū)域的土壤氣地球化學(xué)特征。

結(jié)合1970—2023年臨潭—宕昌斷裂帶地震活動分段時序圖（圖7）與地震分段分布表（表2）來看，該斷裂帶53年間發(fā)生過2.0級以上地震376次。其中，西段發(fā)生地震171次，包括5.0～7.0級地震7次，4.0～4.9級地震4次，3.0～3.9級地震35次，2.0～2.9級地震125次;東段發(fā)生地震205次，包括4.0～4.9級地震9次，3.0～3.9級地震30次，2.0～2.9級地震166次。總體上地震活動以中小型為主，5.0級以上地震較少。根據(jù)震源深度分布（圖8）來看，斷裂西段震源深度主要集中在-10～-5 km，小震多為強震的余震，多數(shù)發(fā)生在-10 km以上;斷裂東段震源在-30～-5 km均有分布，M_S4.0～4.9地震主要發(fā)生于-20～-10 km深度，無M_S≥5.0地震發(fā)生，由此也可以看出斷裂帶具有明顯的分段特性。結(jié)合斷裂帶在岷縣西北方向存在多條次級斷裂，在岷縣東南方向以主斷裂為主，推測斷裂帶的應(yīng)力分布在岷縣區(qū)域發(fā)生改變。其中，岷縣以西的斷裂段強震頻發(fā)，7次M_S≥5.0地震主要發(fā)生在西段，同年發(fā)生過多次M_S＜5.0的余震，其余時間段內(nèi)小震活動較少。如圖6所示，多次強震震源位置相近，均發(fā)生在大地?zé)崃鞯臋M向梯度帶，此區(qū)域斷裂形變強度較大。岷縣以東的斷裂段53年間未發(fā)生過M_S≥5.0地震，但小震活動比西段頻繁，地震發(fā)生區(qū)域更廣泛，此段大地?zé)崃髦递^低，斷裂形變強度較小。

將震源位置（圖8）與測線土壤氣濃度強度（圖5）進行對比可以看出，斷裂帶東段M_S2.0～4.9地震頻繁發(fā)生，應(yīng)力積累周期短，地下斷層未完全閉合，地下流體傳輸間隙繁多。在斷層活動和地下流體的相互作用下發(fā)育多條裂隙，再加上地震發(fā)育的成型裂隙通道，加強了裂隙發(fā)育程度，導(dǎo)致流體活動性增強。斷層處Hg與Rn經(jīng)過長期的逸散后濃度較高，結(jié)合小震發(fā)生頻次、滑動速率走向與大地?zé)崃髦?，推測東段斷層活動穩(wěn)定，當(dāng)前產(chǎn)生的應(yīng)變速率低，容易將應(yīng)力快速釋放，形成多次中小地震。強震主要發(fā)生在斷裂帶西段，包括2013年岷縣—漳縣M_S6.6地震和2019年夏河M_S5.7地震，且在主震后發(fā)生了多起余震。斷裂帶西段通過強震將能量釋放后，地下斷層處于應(yīng)力重新積累狀態(tài)，小震活動的頻率降低，地下流體活動強度相對較低，土壤氣濃度強度趨于穩(wěn)定。Hg、Rn受自身特性影響，向上逸出速率小，無法長時間積累，濃度強度相對較低。但西段位于高熱流值區(qū)域，地下熱能較高，又具有左旋走滑的運動分量，應(yīng)力累積速率快，未來仍有發(fā)生強震的風(fēng)險。李晨樺等^［51］對西秦嶺北緣斷裂帶土壤氣進行監(jiān)測分析也發(fā)現(xiàn)：研究區(qū)中段是張性盆地區(qū)，小震頻繁，斷層的形變強度更大，表現(xiàn)出較高的土壤氣濃度強度;東、西兩段小震活動較弱，相對應(yīng)的斷層土壤氣濃度強度也較低。這種分布特征與地下氣體活動性和斷層形變特征密切相關(guān)，說明小震活動對土壤氣濃度的分布具有一定的影響。

綜上所述，臨潭—宕昌斷裂帶斷層深部的流體活動、大地?zé)崃鞯姆植继卣?，以及斷層滑動速率特征均與地震活動有緊密聯(lián)系。高熱流分布區(qū)域與左旋走滑的形變特征有利于強震的孕育，表現(xiàn)為強震多數(shù)發(fā)生于斷裂帶西段，之后強震所釋放的能量需長時間的積累，對地下裂隙的形變影響不穩(wěn)定，深部氣體上升受阻，表現(xiàn)出土壤氣濃度強度略低。斷裂帶東段一直處于較活躍狀態(tài)，熱流范圍較低，無法為強震孕育提供合適的環(huán)境，導(dǎo)致應(yīng)力釋放周期短，小震活動頻繁。由于長期斷層滑動對裂隙通道的影響，此段深部氣體遷移通道穩(wěn)定，表現(xiàn)出土壤氣濃度強度略高。此外，強震的孕育受到更復(fù)雜的空間分布影響，近年來發(fā)生強震區(qū)域的氣體濃度可能產(chǎn)生異常，還需要對臨潭—宕昌斷裂帶進行持續(xù)的觀察與研究。

5 結(jié)論

本文選取臨潭—宕昌斷裂帶西段4條測線、東段5條測線進行土壤氣Hg、Rn濃度的監(jiān)測，并結(jié)合前人的研究資料與野外測量結(jié)果，對該斷裂帶的土壤氣地球化學(xué)特征與空間分布特性進行分析，得出以下結(jié)論：

（1） 臨潭—宕昌斷裂帶在岷縣西北方向存在多條次級斷裂，在岷縣東南方向以主斷裂為主，結(jié)合大地?zé)崃鞣植继卣骺梢詫嗔褞б葬嚎h為界分為東、西兩段，兩段的應(yīng)力積累環(huán)境不同，地震活動也不同。其中，西段多發(fā)生中強震，但地震間隔時間長，斷裂逐漸閉合，影響土壤氣逸出，表現(xiàn)為西段多條測線Hg、Rn濃度略低;東段一直具有較強的活動性，小震活躍，長期的斷層滑動增強斷裂破裂程度，為土壤氣逸出提供穩(wěn)定的通道，表現(xiàn)為東段Hg、Rn濃度相對于西段略高。

（2） 臨潭—宕昌斷裂帶西段位于高地?zé)釁^(qū)域，而高熱流區(qū)域與熱流存在梯度變化的區(qū)域發(fā)生強震的風(fēng)險較大，1970—2023年間M_S≥5.0地震均發(fā)生在此區(qū)域，但小震發(fā)生頻次較低，多處測線的Hg、Rn濃度強度略低。斷裂帶東段熱流偏低且穩(wěn)定，長年未發(fā)生過強震，但小震活動頻繁，此區(qū)域多處測線Hg、Rn濃度強度略高。

（3） 通過對臨潭—宕昌斷裂帶斷裂活動構(gòu)造、斷層滑動速率與歷史地震等方面的綜合分析，發(fā)現(xiàn)該斷裂帶土壤氣濃度強度、大地?zé)崃鞣植紖^(qū)域與斷層滑動速率分布特征有較高的耦合性。具有較高大地?zé)崃鞣植嫉幕顒訕?gòu)造區(qū)域更容易將應(yīng)力長時間積累，形成強震。小震長期活動頻繁的區(qū)域地下裂隙穩(wěn)定，有利于地下流體活動。這為研究區(qū)深部地下流體活動提供了新證據(jù)，對追蹤未來可能發(fā)生地震的斷層和活動段具有重要意義。

致謝：中國地震局蘭州地震研究所張波副研究員在本論文工作測線布設(shè)方面給予了幫助，在此表示感謝！

參考文獻（References）

［1］杜建國，康春麗.地震地下流體發(fā)展概述［J］.地震，2000，20（增刊1）：107-114.

DU Jianguo，KANG Chunli.A brief review on study of earthquake-caused change of underground fluid［J］.Earthquake，2000，20（Suppl01）：107-114.

［2］李營，陳志，胡樂，等.流體地球化學(xué)進展及其在地震預(yù)測研究中的應(yīng)用［J］.科學(xué)通報，2022，67（13）：1404-1420.

LI Ying，CHEN Zhi，HU Le，et al.Advances in seismic fluid geochemistry and its application in earthquake forecasting［J］.Chinese Science Bulletin，2022，67（13）：1404-1420.

［3］ETIOPE G，MARTINELLI G.Migration of carrier and trace gases in the geosphere：an overview［J］.Physics of the Earth and Planetary Interiors，2002，129（3-4）：185-204.

［4］GHOSH D，DEB A，SENGUPTA R.Anomalous radon emission as precursor of earthquake［J］.Journal of Applied Geophysics，2009，69（2）：67-81.

［5］SCIARRA A，CANTUCCI B，COLTORTI M.Learning from soil gas change and isotopic signatures during 2012 Emilia seismic sequence［J］.Scientific Reports，2017，7（1）：14187.

［6］CHEN Z，LI Y，LIU Z，et al.Evidence of multiple sources of soil gas in the Tangshan fault zone，North China［J］.Geofluids，2019（3）：1945450.

［7］SAHOO S K，KATLAMUDI M，BARMAN C，et al.Identification of earthquake precursors in soil radon-222 data of Kutch，Gujarat，India using empirical mode decomposition based Hilbert Huang Transform［J］.Journal of Environmental Radioactivity，2020，222（7）：106353.

［8］CHEN Z，LI Y，MARTINELLI G，et al.Spatial and temporal variations of CO₂ emissions from the active fault zones in the capital area of China［J］.Applied Geochemistry，2020，112：104489.

［9］蘇鶴軍，周慧玲，李晨樺，等.門源M_S6.9地震前嘉峪關(guān)氣氡變化特征分析［J］.地震工程學(xué)報，2022，44（2）：389-397.

SU Hejun，ZHOU Huiling，LI Chenhua，et al.Characteristics of gas radon variation in Jiayuguan before Menyuan M_S6.9 earthquake［J］.China Earthquake Engineering Journal，2022，44（2）：389-397.

［10］ZHANG L，LIU Y W，GUO L S，et al.Isotope geochemistry of mercury and its relation to earthquake in the Wenchuan Earthquake Fault Scientific Drilling Project Hole-1 （WFSD-1）［J］.Tectonophysics，2014，619：79-85.

［11］DI LUCCIO F，VENTURA G，DI GIOVAMBATTISTA R，et al.Normal faults and thrusts reactivated by deep fluids：the 6 April 2009 M_W6.3 L′Aquila earthquake，central Italy［J］.Journal of Geophysical Research （Solid Earth），2010，115（B6）：B06315.

［12］WALIA V，YANG T F，HONG W L，et al.Geochemical variation of soil-gas composition for fault trace and earthquake precursory studies along the Hsincheng fault in NW Taiwan［J］.Applied Radiation and Isotopes，2009，67（10）：1855-1863.

［13］IRWIN W P，BARNES I.Tectonic relations of carbon dioxide discharges and earthquakes［J］.Journal of Geophysical Research：Solid Earth，1980，85（B6）：3115-3121.

［14］FU C C，YANG T F，CHEN C H，et al.Spatial and temporal anomalies of soil gas in northern Taiwan and its tectonic and seismic implications［J］.Journal of Asian Earth Sciences，2017，149：64-77.

［15］ZHOU X C，DU J G，CHENG Z，et al.Advance review of seismic geochemistry［J］.Bulletin of Mineralogy Petrology and Geochemistry，2012，31 （4）：340-346.

［16］李光濤，程理，吳昊，等.臨潭—宕昌主干斷裂南東段晚第四紀(jì)活動的地質(zhì)地貌證據(jù)［J］.地震工程學(xué)報，2020，42（2）：376-383.

LI Guangtao，CHENG Li，WU Hao，et al.Geological and geomorphological evidence of Late Quaternary activity along the southeastern segment of the Lintan—Tanchang major fault，China［J］.China Earthquake Engineering Journal，2020，42（2）：376-383.

［17］顧功敘.中國地震目錄（公元前1831年—公元1969年）［M］.北京：科學(xué)出版社，1983.

GU Gongxu.Catalogue of Chinese earthquakes （1831 BC—1969 AD）［M］.Beijing：Science Press，1983.

［18］何文貴，周志宇，馬爾曼，等.岷縣—卓尼5.0級地震的基本特征和地質(zhì)背景研究［J］.地震研究，2006，29（4）：373-378.

HE Wengui，ZHOU Zhiyu，MA Erman，et al.Basic features and geological background of the Minxian—Zhuoni M5.0 earthquake on Sep. 7，2004［J］.Journal of Seismological Research，2006，29（4）：373-378.

［19］鄭文俊，雷中生，袁道陽，等.1837年甘肅岷縣北6級地震考證與發(fā)震構(gòu)造分析［J］.地震，2007，27（1）：120-130.

ZHENG Wenjun，LEI Zhongsheng，YUAN Daoyang，et al.Structural research on the 1837 northern Minxian M6 earthquake in Gansu Province and its causative structure［J］.Earthquake，2007，27（1）：120-130.

［20］高明星，陳桂華，徐錫偉.地貌參數(shù)指示的臨潭-宕昌斷裂帶最新構(gòu)造隆升差異與地震活動［J］.地震地質(zhì)，2015，37（3）：709-718.

GAO Mingxing，CHEN Guihua，XU Xiwei.Geomorphic indices indicated recent differential tectonic uplift of the Lintan—Dangchang fault and the Minxian—Zhangxian earthquake［J］.Seismology and Geology，2015，37（3）：709-718.

［21］高鵬偉.甘肅東南地區(qū)現(xiàn)今地殼形變及應(yīng)力特征研究［D］.西安：西安科技大學(xué)，2020.

GAO Pengwei.Study on current crustal deformation and stress characteristics in southeast Gansu Province［D］.Xi'an：Xi'an University of Science and Technology，2020.

［22］鄭文俊，閔偉，何文貴，等.2013年甘肅岷縣漳縣6.6級地震震害分布特征及發(fā)震構(gòu)造分析［J］.地震地質(zhì)，2013，35（3）：604-615.

ZHENG Wenjun，MIN Wei，HE Wengui，et al.Distribution of the related disaster and the causative tectonic of the Minxian—Zhanxian M_S6.6 earthquake on July 22，2013，Gansu，China［J］.Seismology and Geology，2013，35（3）：604-615.

［23］張波，田勤儉，王愛國，等.西秦嶺臨潭—宕昌斷裂第四紀(jì)最新活動特征［J］.地震地質(zhì)，2021，43（1）：72-91.

ZHANG Bo，TIANQinjian，WANG Aiguo，et al.Studies on new activity of Lintan—Dangchang fault，West Qinling［J］.Seismology and Geology，2021，43（1）：72-91.

［24］鄭文俊，袁道陽，何文貴，等.甘肅東南地區(qū)構(gòu)造活動與2013年岷縣—漳縣M_S6.6級地震孕震機制［J］.地球物理學(xué)報，2013，56（12）：4058-4071.

ZHENG Wenjun，YUAN Daoyang，HE Wengui，et al.Geometric pattern and active tectonics in Southeastern Gansu Province：discussion on seismogenic mechanism of the Minxian—Zhangxian M_S6.6 earthquake on July 22，2013［J］.Chinese Journal of Geophysics，2013，56（12）：4058-4071.

［25］李光濤，蘇剛，程理，等.臨潭—宕昌斷裂帶遙感解譯及其晚第四紀(jì)活動特征研究［J］.防災(zāi)科技學(xué)院學(xué)報，2019，21（2）：21-28.

LI Guangtao，SU Gang，CHENG Li，et al.Remote sensing interpretation of Lintan—Dangchang fault zone and its Late Quaternary activity characteristics［J］.Journal of Institute of Disaster Prevention，2019，21（2）：21-28.

［26］劉廷.西秦嶺臨潭—宕昌斷裂帶斷層泥及年代學(xué)研究［D］.北京：中國地質(zhì)大學(xué)（北京），2015.

LIU Ting.Study on fault mud and chronology of Lintan—Tanchang fault zone in West Qinling［D］.Beijing：China University of Geosciences，2015.

［27］袁道陽，張培震，劉百篪，等.青藏高原東北緣晚第四紀(jì)活動構(gòu)造的幾何圖像與構(gòu)造轉(zhuǎn)換［J］.地質(zhì)學(xué)報，2004，78（2）：270-278.

YUAN Daoyang，ZHANG Peizhen，LIU Baichi，et al.Geometrical imagery and tectonic transformation of Late Quaternary active tectonics in northeastern margin of Qinghai—Xizang Plateau［J］.Acta Geologica Sinica，2004，78（2）：270-278.

［28］劉雷，邱江濤，莊文泉.臨潭—宕昌斷裂現(xiàn)今閉鎖程度及地震危險性分析［J］.地震地磁觀測與研究，2021，42（2）：145-146.

LIU Lei，QIU Jiangtao，ZHUANG Wenquan.Fault locking and seismic risk analysis of Lintan—Tanchang fault［J］.Seismological and Geomagnetic Observation and Research，2021，42（2）：145-146.

［29］ANNUNZIATELLIS A，BEAUBIEN S E，BIGI S，et al.Gas migration along fault systems and through the vadose zone in the Latera caldera （central Italy）：implications for CO₂ geological storage［J］.International Journal of Greenhouse Gas Control，2008，2（3）：353-372.

［30］張西娟，曾慶利，馬寅生.斷裂帶中的流體活動及其作用［J］.西北地震學(xué)報，2006，28（3）：274-279.

ZHANG Xijuan，ZENG Qingli，MA Yinsheng.Fluid involvement in active fault zone［J］.Northwestern Seismologica Journal，2006，28（3）：274-279.

［31］CIOTOLI G，LOMBARDI S，ANNUNZIATELLIS A.Geostatistical analysis of soil gas data in a high seismic intermontane basin：Fucino Plain，central Italy［J］.Journal of Geophysical Research （Solid Earth），2007，112（B5）：B05407.

［32］POLLACK H N，HURTER S J，JOHNSON J R.Heat flow from the Earth's interior：analysis of the global data set［J］.Reviews of Geophysics，1993，31（3）：267-280.

［33］王丹，魏水建，賈躍瑋，等.地?zé)豳Y源地震勘探方法綜述［J］.物探與化探，2015，39（2）：253-261.

WANG Dan，WEI Shuijian，JIA Yuewei，et al.An overview of methods for geothermal seismic exploration［J］.Geophysical and Geochemical Exploration，2015，39（2）：253-261.

［34］王貴玲，劉彥廣，朱喜，等.中國地?zé)豳Y源現(xiàn)狀及發(fā)展趨勢［J］.地學(xué)前緣，2020，27（1）：1-9.

WANG Guiling，LIU Yanguang，ZHU Xi，et al.The status and development trend of geothermal resources in China［J］.Earth Science Frontiers，2020，27（1）：1-9.

［35］PASQUALE V，VERDOYA M，CHIOZZI P，et al.Rheology and seismotectonic regime in the northern central Mediterranean［J］.Tectonophysics，1997，270（3-4）：239-257.

［36］安鎮(zhèn)文，朱傳鎮(zhèn).地?zé)崤c地震關(guān)系的研究（一）：彈性介質(zhì)中橢圓盤狀裂紋尖端附近的熱應(yīng)力分布［J］.地震學(xué)報，1984，6（2）：211-222.

AN Zhenwen，ZHU Chuanzhen.Study of geothermal variations in relation to large earthquakes （1）：the distribution of thermal stresses in the neighborhood of the sharp end of an elliptical-shaped crack in an infinite elastic medium［J］.Acta Seismologica Sinica，1984，6（2）：211-222.

［37］姜光政，高堋，饒松，等.中國大陸地區(qū)大地?zé)崃鲾?shù)據(jù)匯編（第四版）［J］.地球物理學(xué)報，2016，59（8）：2892-2910.

JIANG Guangzheng，GAO Peng，RAO Song，et al.Compilation of heat flow data in the continental area of China （4th edition）［J］.Chinese Journal of Geophysics，2016，59（8）：2892-2910.

［38］LOWENSTERN J B，BERGFELD D，EVANS W C，et al.Origins of geothermal gases at Yellowstone［J］.Journal of Volcanology and Geothermal Research，2015，302：87-101.

［39］天嬌，龐忠和，李義曼，等.地?zé)釟怏w研究進展［J］.地質(zhì)學(xué)報，2022，96（5）：1752-1766.

TIAN Jiao，PANG Zhonghe，LI Yiman，et al.Research progress on geothermal gas［J］.Acta Geologica Sinica，2022，96（5）：1752-1766.

［40］馮新斌，陳玖斌，付學(xué)吾，等.汞的環(huán)境地球化學(xué)研究進展［J］.礦物巖石地球化學(xué)通報，2013，32（5）：503-530.

FENG Xinbin，CHEN Jiubin，F(xiàn)U Xuewu，et al.Progresses on environmental geochemistry of mercury［J］.Bulletin of Mineralogy，Petrology and Geochemistry，2013，32（5）：503-530.

［41］袁伏全，張超美，孫世瑞，等.青海地區(qū)地震與地?zé)岬姆植继卣鳎跩］.高原地震，2017，29（2）：1-6.

YUAN Fuquan，ZHANG Chaomei，SUN Shirui，et al.Study on spatial relationship between earthquakes and geothermal in Qinghai Province［J］.Plateau Earthquake Research，2017，29（2）：1-6.

［42］馮建剛，張輝，楊萍.2013年岷縣漳縣6.6級地震前地震b值異常特征研究［J］.地震，2016，36（1）：32-37.

FENG Jiangang，ZHANG Hui，YANG Ping.Anormalies of b-value changes before the 2013 Minxian—Zhangxian M_S6.6 earthquake［J］.Earthquake，2016，36（1）：32-37.

［43］易桂喜，聞學(xué)澤，蘇有錦.川滇活動地塊東邊界強震危險性研究［J］.地球物理學(xué)報，2008，51（6）：1719-1725.

YIGuixi，WEN Xueze，SU Youjin.Study on the potential strong-earthquake risk for the eastern boundary of the Sichuan—Yunnan active faulted-block，China［J］.Chinese Journal of Geophysics，2008，51（6）：1719-1725.

［44］WAKITA H，NAKAMURA Y，KITA I，et al.Hydrogen release：new indicator of fault activity［J］.Science，1980，210（4466）：188-190.

［45］吳乾蕃，祖金華，謝毅真，等.云南地區(qū)地?zé)峄咎卣鳎跩］.地震地質(zhì)，1988，10（4）：177-183.

WU Qianfan，ZU Jinhua，XIE Yizhen，et al.Characteristics of geothermal field in Yunnan region［J］.Seismology and geology，1988，10（4）：177-183.

［46］張之立.唐山地震熱力源的研究［J］.地震學(xué)報，1985，7（1）：45-56.

ZHANG Zhili.Investigation on the thermal force sources of 1976 Tangshan earthquake［J］.Acta Seismologica Sinica，1985，7（1）：45-56.

［47］KING C Y，ZHANG W，ZHANG Z C.Earthquake-induced groundwater and gas changes［J］.Pure and Applied Geophysics，2006，163（4）：633-645.

［48］張培震，李傳友，毛鳳英.河流階地演化與走滑斷裂滑動速率［J］.地震地質(zhì)，2008，30（1）：44-57.

ZHANG Peizhen，LI Chuanyou，MAO Fengying.Strath terrace formation and strike-slip faulting［J］.Seismology and Geology，2008，30（1）：44-57.

［49］張培震，徐錫偉，聞學(xué)澤，等.2008年汶川8.0級地震發(fā)震斷裂的滑動速率、復(fù)發(fā)周期和構(gòu)造成因［J］.地球物理學(xué)報，2008，51（4）：1066-1073.

ZHANG Peizhen，XU Xiwei，WEN Xueze，et al.Slip rates and recurrence intervals of the Longmen Shan active fault zone，and tectonic implications for the mechanism of the May 12 Wenchuan earthquake，2008，Sichuan，China［J］.Chinese Journal of Geophysics，2008，51（4）：1066-1073.

［50］CHEN Z，LI Y，LIU Z F，et al.Geochemical and geophysical effects of tectonic activity in faulted areas of the North China Craton［J］.Chemical Geology，2022，609：121048.

［51］李晨樺，張慧，蘇鶴軍，等.西秦嶺北緣斷裂帶深部氣體地球化學(xué)特征與斷層形變空間分布相關(guān)性研究［J］.地震工程學(xué)報，2016，38（6）：955-963.

LI Chenhua，ZHANG Hui，SU Hejun，et al.Study on correlation between geochemical features of deep gas and fault deformation distribution in northern margin fault zone of West Qinling［J］.China Earthquake Engineering Journal，2016，38（6）：955-963.

（本文編輯：趙乘程）