基于Matlab的斷裂帶溫泉水地球化學(xué)特征及地震活動性研究①

伍劍波，張 慧，2，蘇鶴軍，2，李晨樺，2

（1．中國地震局蘭州地震研究所，甘肅 蘭州 730000；2．中國地震局地震預(yù)測研究所蘭州科技創(chuàng)新基地，甘肅 蘭州 730000）

0 引言

地下流體地球化學(xué)方法廣泛應(yīng)用于地震領(lǐng)域，是研究斷裂帶特征以及活動性的一個重要手段。構(gòu)造活動形成的裂隙為溫泉水循環(huán)提供了良好的地質(zhì)通道，流體在上升過程中攜帶了大量的地質(zhì)信息，對構(gòu)造斷裂活動性研究有很重要的地質(zhì)意義。大量的研究表明［1－4］，沿斷裂帶分布的溫泉水其水化特征變化與斷裂活動具有一定的關(guān)聯(lián)性。

隨著地震研究的深入和地下流體學(xué)科的發(fā)展，需要更多計算機軟件技術(shù)的支撐發(fā)展。在地震地下流體研究中，對野外勘查、實驗數(shù)據(jù)分析以及圖形的繪制，需要采用多種分析手段，由于采用的分析方法和手段不同，繪圖軟件就不同，所以研究人員需要掌握多種繪圖軟件。Matlab作為世界三大數(shù)學(xué)軟件之一，具有強大的計算和繪圖能力以及簡單的程序語言，在地震研究方面有重要運用［5］。為了解決地下流體分析軟件分散的局面，我們基于Matlab開發(fā)出了一種圖形用戶界面（GUI）軟件（圖1），該軟件綜合了大部分地球化學(xué)分析軟件的繪圖功能。本文主要探討該軟件水質(zhì)分析功能，并以西秦嶺北緣斷裂帶溫泉水為分析對象進行分析。

圖1 基于Matlab的水質(zhì)分析軟件界面Fig．1 The software interface of water quality analysis based on Matlab．

1 應(yīng)用區(qū)域構(gòu)造概況

西秦嶺北緣斷裂帶位于甘肅東南部。研究區(qū)范圍在北緯33．50°～35．58°，東經(jīng)103．50°～106．50°（圖2）；地震數(shù)據(jù)來自國家地震科學(xué)數(shù)據(jù)共享中心，時間段為1970年1月1日－2012年12月31日；震級ML≥2。

甘東南區(qū)域的構(gòu)造環(huán)境復(fù)雜，是中強震和強震發(fā)生的主要場所，區(qū)域內(nèi)的斷裂構(gòu)造特征嚴(yán)重受到區(qū)域構(gòu)造的影響［6］。西秦嶺北緣斷裂帶為甘東南區(qū)域活動性較強的斷裂之一，晚第四紀(jì)以來，主要以左旋走滑為主，傾滑分量較小，斷裂新活動明顯［7］（圖2）。西秦嶺北緣斷裂帶分為六個次級斷裂，多為張扭性斷層，其派生的小斷層發(fā)育且相互交錯［8］。本文研究的街子溫泉位于天水?dāng)嗔?，武山溫泉位于武山斷裂，而清水溫泉?jù)推測可能位于NW向的通渭－清水隱伏斷裂帶上［9］，不屬于西秦嶺北緣主斷裂系。

2 斷裂帶溫泉水水質(zhì)分析

溫泉受區(qū)域主要斷裂控制，其分布特征往往與地震帶分布一致，通常被認為是確定活動斷裂的一種標(biāo)志［10］。本文把我們采集的溫泉水樣分析數(shù)據(jù)（數(shù)據(jù)III：2011年采集樣品），與前人（數(shù)據(jù)Ⅰ：石雅镠，1987年采集樣品［2］；數(shù)據(jù)Ⅱ：汪萬紅，2006年采集樣品［4］）進行對比研究（表1）。基于 MATLAB繪制矩形水化學(xué)圖、δD－δ18O關(guān)系圖、Na－K－Mg三角圖、以及溫泉水循環(huán)深度三維圖。

圖2 西秦嶺北緣斷裂帶及地震分布Fig．2 Distribution of the north margin fault belt of western Qinling and epicenters of erathquake．

表1 溫泉水樣分析數(shù)據(jù)（單位：mg／L）Table 1 Analysis data of the hot spring water

2．1 溫泉水水化學(xué)類型

矩形水化學(xué)類型圖［11］（圖3）對揭示水文地球化學(xué)規(guī)律起著舉足輕重的作用。由表1幾種主要陰陽離子的濃度數(shù)據(jù)以及圖3顯示結(jié)果，可推測武山和街子溫泉水的總礦化度小于清水溫泉。1987年和2011年的溫泉水化學(xué)組成總體上沒有發(fā)生明顯變化，為低礦化度、淡水。三個溫泉水陽離子含量均以Na＋為主，Ca2＋、Mg2＋較少。但是三者陰離子含量有明顯的區(qū)別，清水溫泉水陰離子以SO2－4為主，而街子和武山溫泉水中陰離子HCO－3含量最多、SO2－4次之。根據(jù)矩形圖分類方法，清水溫泉水化學(xué)類型為Na－SO4·Cl，屬于硫酸型水；街子、武山溫泉水化學(xué)類型為Na－HCO3·SO4·Cl，屬于碳酸型水。

2．2 氫氧同位素分析

最早由Craig［12］提出全球淡水中δD和δ18O具有線性相關(guān)性，并運用全球降水線定義了二者之間的關(guān)系：δD＝8δ18O＋10。Craig大氣降水線只適用全球范圍，由于氣候和地理參數(shù)的變化，局部區(qū)域大氣降水線的斜率和氘截距都與全球線不同。前人經(jīng)過地區(qū)雨水線的建立，求得中國西北地區(qū)的大氣降水線［13］：δD＝7．38δ18O＋7．16。本文以西北地區(qū)的大氣降水線為標(biāo)準(zhǔn)，1987年與2006年的街子、武山、清水三個溫泉水的δD和δ18O線性相關(guān)性如圖4。

圖3 溫泉水矩形水化學(xué)圖Fig．3 Rectangle hydrochemical diagram of the hot spring waters．

圖4 西秦嶺北緣溫泉水δD～δ18 O關(guān)系圖Fig．4 The correlation diagram ofδDandδ18 O of hot spring water samples in the north margin of western Qinling fault zone．

首先，1987年和2006年表現(xiàn)出共同點：溫泉水的氫氧同位素組成均沿著西北地區(qū)大氣降水線附近分布，表明溫泉水補給來源為大氣降水；溫泉水18O總體上相對于降水線向右偏移，即正向氧漂移。由于三個溫泉圍巖的巖性主要為片麻巖、花崗巖，其中街子溫泉的圍巖還含有碳酸鹽巖，地下水與含硅酸鹽礦物或碳酸鹽礦物的巖石之間發(fā)生水—巖反應(yīng)，進行氧同位素平衡交換［14］，即 Si18O2＋ 2H216O ＝Si16O2＋ 2H218O，或 CaC18O3＋H216O ＝ CaC16O3＋H218O，致使地下水18O含量變化。

其次，2006年與1987年相比有明顯的差異性：街子、武山、清水三個溫泉水2006年的δD、δ18O值比1987年更偏向負值一端，說明隨著時間的推移溫泉水發(fā)生了不同程度的氫、氧同位素虧損；2006年的斜率與西北地區(qū)大氣降水線有明顯的偏離。從溫泉水來源上分析，造成這種差異的原因可能是氣候年變化對大氣降水同位素組成的影響。

2．3 水—巖平衡狀態(tài)

Na－K－Mg三角圖解法［15］廣泛應(yīng)用于評價地下熱水的水—巖平衡狀態(tài)和區(qū)分不同類型的水樣。由完全平衡線和部分平衡下限把整個三角圖分為完全平衡、部分平衡和未成熟水三個區(qū)。

由圖5可知，水—巖反應(yīng)在1987與2011年之間有明顯的變化：2011年溫泉水樣比對應(yīng)的1987年樣品更接近完全平衡線。由此推測隨著時間的推移溫泉水的成熟度越高。

三個溫泉均靠近Mg端元，其中武山（WS－Ⅰ、WS－Ⅲ）和街子（JZ－Ⅰ、JZ－Ⅲ）兩個溫泉的水樣接近完全平衡線，說明二者的水—巖反應(yīng)幾乎達到完全平衡狀態(tài)，成熟度較高；而清水溫泉的水樣（QS－Ⅰ、QS－Ⅲ）明顯位于完全平衡線與部分平衡下限之間，屬于部分成熟的水，成熟度小于前兩者。鑒于3個溫泉水具有較高的成熟度，由圖中樣品落點可初步推測溫泉的熱儲溫度：武山為120～160℃、街子為140～180℃，清水為160～180℃。本文研究結(jié)果與石雅謬［2］利用SiO2地?zé)釡貥?biāo)計算得到的熱儲溫度（武山139℃、街子159℃、清水190℃）基本一致。

圖5 溫泉水Na－K－Mg三角圖Fig．5 The triangle diagram of Na－K－Mg of the hot spring water samples．

2．4 溫泉水循環(huán)深度特征

溫泉的熱源如果來源于地?zé)嵩鰷責(zé)?，且熱儲溫度隨著循環(huán)深度的增加而升高，可用下式計算溫泉水循環(huán)深度：

式中D為循環(huán)深度（km）；g為地溫梯度（℃／km）；tR為熱儲溫度（℃）；tcold為當(dāng)?shù)刈罾淅淙臏囟龋ā妫?；h為常溫層厚度（km）。汪萬紅［4］運用式（1）計算了秦嶺北緣斷裂溫泉水循環(huán)深度（表2）。作者把表2數(shù)據(jù)和溫泉位置的經(jīng)緯度導(dǎo)入水質(zhì)分析界面（圖1），繪出三維的秦嶺北緣溫泉水循環(huán)界面（圖6）。圖中顯示了秦嶺北緣斷裂的溫泉水循環(huán)界面。秦嶺北緣的溫泉水循環(huán)深度陜西段最深；甘肅段循環(huán)深度最淺；而青海段與甘肅段相差不大。整個斷裂水循環(huán)界面以甘肅中段為界，呈西高東低，西秦嶺北緣斷裂的溫泉水循環(huán)深度比秦嶺北緣東段小。由表2和圖6可知清水溫泉的循環(huán)深度最深，武山和街子溫泉循環(huán)深度較淺一些。

圖6 秦嶺北緣斷裂帶的溫泉水地下循環(huán)界面Fig．6 The underground circulation interface of the hot spring waters in the north margin of Qinling fault zone．

表2 秦嶺北緣斷裂溫泉水循環(huán)深度［4］Table 2 The circulation depth of the hot spring water in the north margin of Qinling fault zone［4］

3 溫泉水地球化學(xué)特征與斷裂活動性的關(guān)系

本文研究的街子溫泉、武山溫泉分別位于西秦嶺北緣斷裂的天水?dāng)嗔讯魏臀渖綌嗔讯危▓D2），而清水溫泉據(jù)推測可能位于北西向的通渭—清水隱伏斷裂帶上，不屬于西秦嶺北緣斷裂。根據(jù)歷史地震資料和相關(guān)文獻查閱結(jié)果［7－8］，天水?dāng)嗔讯螢?34年天水7級地震的主要發(fā)震斷裂段；武山斷裂可能曾發(fā)生公元前47年隴西63/4級、128年甘谷61/2級和1765年甘谷—武山6 1/2級等地震；通渭—清水隱伏斷裂可能是1718年通渭7 1/2級地震的主要發(fā)震斷裂段。1970年到2012年研究區(qū)震中分布顯示武山段最高、天水段次之、而清水發(fā)震頻數(shù)最少（圖2）。因此從歷史地震和現(xiàn)代地震可初步推測：武山斷裂活動性最強、天水?dāng)嗔汛沃?，通渭－清水隱伏斷裂活動性最弱。

上述溫泉水質(zhì)分析得出清水溫泉的循環(huán)深度最大，街子溫泉次之，武山溫泉最??；而且清水溫泉的水化學(xué)組成特征也明顯不同于武山、街子溫泉。說明溫泉水深循環(huán)對斷裂及其圍巖有強烈的介質(zhì)弱化作用，清水溫泉水地下循環(huán)過程中對圍巖的弱化作用最強，降低斷裂面的有效正壓力，導(dǎo)致孕震時期長，地震活動性低；相反，武山、街子溫泉水對圍巖的弱化作用較小些，有效應(yīng)力和能量積累可能較清水快一些，孕震周期短，所以地震活動性高。

4 結(jié)論與討論

隨著流體學(xué)科研究的深入，Matlab的軟件技術(shù)支撐促進了研究方法的多樣化發(fā)展。矩形水質(zhì)分類法、氫氧同位素分析、水－巖平衡以及循環(huán)深度特征分析在地學(xué)研究中有廣泛的應(yīng)用，本文根據(jù)這些方法開發(fā)的分析軟件，在水質(zhì)研究方面取得較好的應(yīng)用效果。

通過不同年份的樣品對比得出，Ⅰ、Ⅱ、Ⅲ類型水樣的水化學(xué)特征總體上表現(xiàn)一致，溫泉水質(zhì)隨著時間推移變化幅度不大，而三個溫泉之間的水質(zhì)特征有明顯區(qū)別。首先，溫泉水氫氧同位素特征表明三個溫泉水補給來源均為大氣降水，同位素交換作用導(dǎo)致18O正向漂移，氣候變化導(dǎo)致同位素組成的差異；其次，通過三個溫泉水對比分析得出，清水溫泉水化學(xué)特征與武山、街子溫泉也有明顯的區(qū)別：清水溫泉水化學(xué)類型為Na－SO4·Cl，為硫酸型水，水－巖反應(yīng)達到部分平衡，街子、武山溫泉水化學(xué)類型為Na－HCO3·SO4·Cl，為碳酸型水，水－巖反應(yīng)近似處于完全平衡狀態(tài)。

從文中分析可知，清水溫泉水循環(huán)深度最深，熱儲溫度最高，對圍巖弱化作用最強，其附近以及通渭—清水隱伏斷裂發(fā)震頻率較低；而武山溫泉水循環(huán)深度最淺，熱儲溫度最小，西秦嶺北緣斷裂武山段地震活動較為頻繁。溫泉水表現(xiàn)出的這種差異性是否與斷裂活動性大小有相關(guān)性，本文都只做了初步討論，需要后來者更深入的研究。

［1］ R O Fournier，A H Truesdell．Geochemical Indicators of Subsurface Temperature——Part2，Estimation of Temperature and Fraction of Hot Water Mixed with Cold Water［J］．Jour．Research U．S．Geol．Survey，1974，2（3）：263－270．

［2］ 石雅镠．秦嶺北緣大斷裂溫泉的地震化學(xué)研究［D］．蘭州：國家地震局蘭州地震研究所，1988．

SHI Ya－liu．Research on the Seismological Geochemistry in the Northern Margin Fault Zone of Qinling［D］．Lanzhou：Lanzhou Institute of Seismology ，CEA，1988．

［3］ 林元武．紅河斷裂北段溫泉水循環(huán)深度與地震活動性的關(guān)系探討［J］．地震地質(zhì)，1993，15（3）：193－199．

LIN Yuan－wu．A discussion on the Relation of Circulation Depth of Hot Spring Water to Seismic Activity on the Northern Segment of the Honghe Fault Zone［J］．Seismology and geology，1993，15（3）：193－199．

［4］ 汪萬紅，張慧，蘇鶴軍．秦嶺北緣斷裂帶溫泉水循環(huán)深度與地震活動性的關(guān)系研究［J］．西北地震學(xué)報，2008，30（1）：36－41．

WANG Wan－h(huán)ong，ZHANG Hui，SU He－jun．Research on the Relation of Circulation Depth of Hot spring Water with Seismicity in the Northern Margin Fault Zone of Qinling［J］．Northwestern Seismological Journal，2008，30（1）：36－41．

［5］ 劉雙慶，聶永安，高武平，等．Matlab在天津市強震動臺網(wǎng)烈度速報中的應(yīng)用［J］．西北地震學(xué)報，2012，34（4）：405－412．

LIU Shuang－qing，NIE Yong－an，GAO Wu－ping，et al．Application of Matlab Software on Intensity Rapid Report of Tianjin Strong Mation Network［J］．Northwestern Seismological Journal，2012，34（4）：405－412．

［6］ 羅國富，屠泓為，馬禾青，等．甘東南至陜甘寧交界中強震危險區(qū)地震活動能量場時空特征［J］．西北地震學(xué)報，2012，34（2）：132－137．

LUO Guo－fu，TU Hong－wei，MA He－qin，et al．Space－time Distribution Characteristics of Energy Field of Seismic Activity for the Southeastern Part of Gansu Province to Shanxi－Gansu－Ningxia Boundary Region Moderate－strong Earthquake Hazard Area［J］．Northwestern Seismological Journal，2012，34（2）：132－137．

［7］ 滕瑞增，金瑤泉，李西侯，等．西秦嶺北緣斷裂帶新活動特征［J］．西北地震學(xué)報，1994，16（2）：85－80．

TENG Rui－zeng，JIN Yao－quan，LI Xi－h(huán)ou，et al．Recent Activity Characteristics of the Fault Zone at Northern Edge of Western Qinling MT［J］．Northwestern Seismological Journal，1994，16（2）：85－80．

［8］ 邵延秀，袁道陽，王愛國，等．西秦嶺北緣斷裂破裂分段與地震危險性評估［J］．地震地質(zhì)，2011，33（1）：79－80．

SHAO Yan－xiu，YUAN Dao－yang，WANG Ai－guo，et al．The Segmentation of Rupture and Estimate of Earthquake Risk along the North Margin of Western Qinling Fault Zone［J］．Seismology and Geology，2011，33（1）：79－80．

［9］ 張守訓(xùn)，李百祥．天水及其南北地區(qū)溫泉分布的地質(zhì)－地球物理特征［J］．西北地震學(xué)報，2006，28（3）：252－257．

ZHANG Shou－xun，LI Bai－xiang．The Geological and Geophysical Characteristics of Hot Spring Distribution in Tianshui and Its NS Adjacent Areas［J］．Northwestern Seismological Journal，2006，28（3）：252－257．

［10］ 鐘以章，肖秀清．遼東半島溫泉與地震空間分布關(guān)系討論［J］．地震地質(zhì)，1990，12（4）：343－350．

ZHONG Yi－zhang，XIAO Xiu－qing．Preliminary Study on Spatial Relationship between Earthquakes and Hot springs on Liaodong Peninsula［J］．Seismology and Geology，1990，12（4）：343－350．

［11］ Chadha D K．A Proposed New Diagram for Geochemical Classification of Natural Waters and Interpretation of Chemical Dat a［J］．Hydrogeology Journal，1999，7：431－439．

［12］ Harmon Craig．Isotopic Variations in Meteoric Waters［J］．Science，1961，133（3465）：1702－1703．

［13］ 高志發(fā)．西北地區(qū)大氣降水、地表水及地下水同位素組成特征探討［J］．甘肅地質(zhì)學(xué)報，1993，2（2）：94－101．

GAO Zhi－fa．Discussion on Feature of Isotope Component from Atmospheric Water，Ground Water and under Ground Water in Northwest Area，China［J］．Acta Geologica Gansu，1993，2（2）：94－101．

［14］ 桂和榮，陳陸望，宋曉梅．皖北礦區(qū)地下水中氫氧穩(wěn)定同位素的漂移特征［J］．哈爾濱工業(yè)大學(xué)學(xué)報，2005，37（1）：111－114．

GUI He－rong，CHEN Lu－wang，SONG Xiao－mei．Drift Features of Oxygen and Hydrogen Stable Isotopes in Deep Groundwater in Mining area of Northern Anhui［J］．Journal of Harbin Institute of Technology，2005，37（1）：111－114．

［15］ Werner F．Giggenbach．Geothermal solute equilibria：derivation of Na－K－ Mg－Ca geoindicators［J］．Geochimica et Cosmochimica Acta，1988，52：2749－2765．