甘東南4井水溫資料在汶川MS8.0地震前后的異常變化特征分析＊

蒲小武 梅東林 陳彥平 葉媛媛 王軍燕 許可娟

學術論文

甘東南4井水溫資料在汶川MS8.0地震前后的異常變化特征分析＊

蒲小武 梅東林 陳彥平 葉媛媛 王軍燕 許可娟

（甘肅省地震局隴南中心地震臺，甘肅隴南746000）

較為系統地分析了甘東南成縣、樊壩、清水溫泉、兩水4口水溫井自2007年以來的水溫觀測數據，發現在2008年5月12日汶川MS8.0地震前后4口水溫井資料都出現了不同程度的前兆異常變化，但前兆異常信息大部分淹沒在趨勢上升的背景中。為了更好地提取前兆異常信息，文中采用了流體資料處理時常用的差分法、變差率法、從屬函數法、趨勢速率法等前兆異常信息提取方法，對數據進行了處理，并分析了4井水溫前兆異常變化特征。

水溫前兆異常；差分法；變差率；從屬函數；趨勢速率

引言

井水溫度動態觀測是地下流體學科重要的觀測項目之一，中國自開展水溫觀測以來，已積累了大量的水溫地震前兆異常實例。對水溫前兆異常特征［1］、干擾因素［2］、動態類型［3－7］等已有許多專家進行過系統研究。此前總結的水溫地震前兆異常的一般特征是：異常多出現在震前幾天至幾十天時間段內，為短臨異常；異常幅度大，一般為正常起伏度的幾至幾十倍，異常顯著而易識別，因此，在地震短臨監測中，表現出特有的優勢［8－10］。這些研究成果為我國的地震監測與預測，尤其對在短臨預測與地震預報方面提供了重要的數據。

在“十五”期間，甘肅甘東南區域，架設了許多水溫觀測系統，極大地促進了本區域地下流體觀測項目的發展。2008年5月12日在青藏塊體東部邊緣發生了汶川MS8.0地震，這次地震前，位于甘東南區域內的成縣、樊壩、清水溫泉、兩水4口水溫井資料出現了異常變化，由于4井水溫多年動態均屬于升溫型形態，異常信息大部分淹沒在趨勢變化中。為了更好地提取前兆異常，我們對4井水溫資料采用了一階差分法、變差率法、從屬函數法、趨勢速率法進行數據處理，并分析了4井水溫地震前兆異常特征。

1 臺站基本情況

1.1觀測井

（1）成縣井。該井位于成縣縣城以東2km的蘭空農廠，觀測含水層為盆地第三系空隙承壓水，由秦嶺山區降水補給，水頭高出地面18m，含水層透水性差，井深264.48m，井內套管深度255.90m。2007年4月開始水溫觀測，使用的儀器是SZW－1A型數字溫度計。該井封閉性較好，年動態穩定，一般降雨對流量沒有影響。

（2）樊壩井。該井位于隴南區甘泉鄉樊壩村，距隴南市55km，井口地面標高180m。地處秦嶺西北端米倉山山脈，是長江流域和黃河流域的分水嶺，白龍江斷裂通過該處。出露的地層為第四系坡洪積亞粘土及沙粒土，下伏泥質灰巖、硬質泥巖。隴南樊壩水溫觀測是“十五”期間新上的流體觀測項目，該井是為了流體觀測新打的一口井，井深達200m，2007年3月安裝了SZW－1A型數字式水溫計，與水位同井觀測。

（3）清水溫泉井。該井位于縣城北東8km的湯浴河谷內，海拔高程1 445m。地勢北東高，南西低，屬低山丘陵。其構造位置屬于隴西系南段、秦嶺北緣大斷裂的北側，六盤山斷裂帶西南側（圖1）。臺基巖性為通渭—清水斷裂晚古生代花崗巖，該溫泉是深部熱水沿巖石裂隙運動，在一定地質構造條件下出露于地表而形成的。清水溫泉井深443.3m。2007年6月開始水溫觀測，使用的儀器是SZW－1A型數字溫度計。

（4）兩水井。該井位于甘肅省隴南區兩水鎮前村，處于武都弧形構造帶前弧西翼，處在多組構造交匯復合部位。區內分布有武都弧形構造帶，東西向構造帶，北北東向、北東向及北西向構造帶，在本區東南方向四川省境內還分布有大致成N60°E方向展布的龍門山構造帶。井深75.92m，處于志留系石灰巖中，該灰巖中有片巖夾層，井內有裂隙水滲入。隴南臺兩水水溫始測于2001年3月，與水位同井觀測，配備了由地殼所研制的SZW－1A型數字式溫度計。

圖1 區域構造背景及井點分布圖

1.2正常動態

對于任何一個地震前兆觀測項目，只有在了解其正常動態特征的基礎上，才能夠有效地提取地震前兆異常信息。井水溫度的背景值，一般取決于井孔所在地的大地熱流q或地溫梯度dT／dz。趙剛等對全國近277個觀測井的地熱正常動態進行了系統研究，將其劃分為降溫型、升溫型、波動型、穩定型、跳變型、長周期型6種基本類型［11－12］。據統計，在整個地下流體臺網中，測點溫度為降溫型、升溫型和波動型占了70%左右。

本文研究區域內的4口水溫井，水溫多年動態都呈上升趨勢變化（圖2）。其中成縣、樊壩、清水溫泉井水溫探頭位置處于恒溫帶下，所測為深部水溫變化，不受地表溫度影響，無明顯干擾源。成縣井年變幅度達到0.04℃，月變幅約為0.003℃。樊壩井年變幅度達到0.009℃，月變幅約為0.000 8℃。清水溫泉井水溫年變幅度達到0.016℃，月變幅約為0.001 5℃。

圖2 汶川MS8.0地震前后4井水溫整點值、年均值對比圖

兩水井水溫自2001年開始觀測以來，基本呈現上升趨勢變化（圖2），年變幅度達到0.01℃，月變幅約為0.000 9℃。長期動態特征屬于升溫型，但由于井孔只有75m左右，探頭位置雖然處于恒溫帶下，但由于水的對流，水溫變化受水位、氣溫波動變化影響，夏秋測值變化較大。

2 流體異常提取方法簡介

一階差分法是一種常規數據處理方法，方法簡單，本文不做介紹。

2.1變差率法

變差率是表征一條曲線的相對變化幅度，就是定量地確定本年相對于前一年同一月的變化程度，如果出現明顯的變化差，則表明本年度變化程度較大，若變差率較小，則表明本年與前一年的年動態相似。其具體做法如下：

（1）以水溫月均值作為基本的數據序列；

（2）定義變差率為某個月的水溫值與前一年同月的水溫值之差再除以前一年同月的水溫值。其計算式為：

式中，Ri為變差率；Ci和Ci－12分別為水溫月均值和前一年同月的水溫月均值。

2.2從屬函數法

盡管前兆觀測資料的異常變化形態多種多樣，諸如突跳、轉折、加速等，但究其本質都是觀測曲線隨著時間變化呈現出的斜率變化，而從屬函數就表征一條曲線變化的斜率和跳動的程度［13］。由于從屬函數的這種特性，因而對處理趨勢變化，定量分析其曲線斜率變化，效果獨特。

從屬函數μ值公式如下：

ki為觀測值變化斜率；ri為月均值M（t）與時間t的相關系數，它反應了觀測序列內在質量的好壞；α為經驗常數。

2.3趨勢速率法

趨勢變化特征是地下流體中短期異常的主要特征之一。由于不同的流體測點所處的構造、含水層巖性、井（泉）結構等物理化學性質的差異，趨勢性異常就有上升、下降和波動之分。本文研究的4井水溫正常背景變化為趨勢上升型，許多前兆異常信息就會淹沒在趨勢上升的背景里，定量計算趨勢速率（曲線斜率）變化，才可以有效提取地震前兆異常。盡管趨勢速率方法與從屬函數方法都表示某一時段內時間序列斜率大小，但兩種方法得到的結果不完全相同，在資料分析運用中可以作為參考。具體做法是：

（1）以月均值為基本數據序列。

（2）計算趨勢速率值（v）。選擇線性窗Y＝a＋bX擬合模型。

（3）判定異常指標。若趨勢速率值v≥0.5則為異常指標。

3 水溫異常特征分析

3.1原始曲線特征分析

兩水井水溫自2007年年初，成縣、樊壩、清水3井水溫自2007年架設以來，一直處于快速上升狀態（圖2），在汶川8.0級地震前后，沒有明顯的突變、轉折等異常形態出現，但水溫上升變化速率在汶川地震前要快于地震后，2007—2008年的水溫曲線斜率大于2009年以后其他年份的，在2010年以前水溫變化總體呈近勻減速變化，這種變化與對數函數曲線形態很類似。2010年后4井水溫變化速度相對穩定，數據以比較緩慢而相對固定的速度上升，曲線形態呈現近直線斜形上升變化。

3.2一階差分異常特征分析

差分法是一種常用的數據處理方法，在本文中差分法所用數據為月均值數據。

成縣、樊壩、清水3井水溫在汶川地震前，一階差分值有明顯的快速下降變化，這種趨勢變化一直持續到2009年初。2009年至2010年成縣、清水兩井水溫一階差分值呈現緩慢下降變化，速度明顯低于2008年以前的，2011年至2012年變化平穩。樊壩井自2009年起一階差分值圍繞某一基值呈小幅波動變化。

兩水井水溫在2007—2012這6年間，一階差分值呈現年初上升，出現最大峰值后下降，在年底出現最小峰值后又回升，年年輪回，呈現準周期變化。這種變化主要是水位的準周期漲落引起水溫的準周期變化。在汶川地震前后，一階差分值沒有出現明顯異常。

3.3變差率法異常特征分析

4井變差率在汶川地震前都出現了幾乎同步、快速下降變化，變化形態非常相似（圖3）。汶川地震后，這種變化狀態都持續了一段時段，但持續時間4井有差異。樊壩井自2008年7月起結束了快速下降變化，出現緩慢波動下降，2011年出現了轉折上升變化。清水井2008年12月結束了快速下降變化狀態，2009年起變差率出現了緩慢下降，2011年出現轉折上升變化。兩水井在汶川地震前變差率總體以波動下降變化為主，地震后出現快速下降變化，2009結束了快速下降變化，以波動變化為主，2011年起出現大幅波動上升變化。成縣井在汶川地震前后都以快速下降為主，2010年5月下降速率減緩。

3.4從屬函數異常特征及趨勢速率異常特征分析

趨勢速率方法與從屬函數方法都表示某一時段內時間序列斜率大小，都是對線性相關性綜合特征的描述。兩種方法得到的結果總體基本相似（圖3）。汶川地震前4井水溫從屬函數μ值和趨勢速率值呈現明顯的高值變化（圖2），且μ值、v值都大于0.5，總體呈下降趨勢。

圖3 4井水溫經4種方法處理后的曲線圖

成縣井μ值2008年2月前超過0.6，3月急速降至0，汶川地震后又急速回升，數值超過0.6，9月份后從屬函數值為0。樊壩井和清水井從屬函數值變化形態相似，汶川地震前后一段時間μ值超過0.5；地震前μ值快速下降，地震發生后μ值轉平，不久又快速下降，2008年底后μ值再次低于0.5。兩水井μ值變化有一定的特點，幾乎呈正弦波形態變化，每年5～6月出現高峰值，當年12月至次年1～2月出現低峰值；2007年、2008年的μ高峰值都大于0.5，明顯大于以后各年份峰值幅度。

4 結論與討論

（1）4井在汶川地震前水溫測值都出現了明顯的快速上升變化，這種變化是漂移還是異常？兩水水溫是2001年開始觀測，資料積累了很多年，也就在汶川地震前存在快速上升變化，在本文中為了和其他水溫井比較所以選取了2007年以來的數據進行分析對比。其他3井的變化形態總體和兩水水溫測值的變化趨勢相似，2007年以來呈現快速上升，但速率卻呈現逐漸下降的變化。這種趨勢上升變化的水溫資料，異常信息都已淹沒在趨勢上升的背景變化里了。本文定量計算了4井水溫一階差分值、月變差率、從屬函數、趨勢速率，發現4井在汶川地震前后有著幾乎相似或相近的變化形態，其值都是震前處于高值區域，而且測值變化是由高向低快速下降。

（2）但對于水溫異常特征，目前多處于定性研究階段，主要是水溫異常與地震的對應關系的定性分析，定量的研究較少。對于長期漂移上升（下降）型，采用一定的數學方法可以更好地提取地震前兆異常。差分法、月變差率法、從屬函數法、趨勢速率法等數學方法就是非常好的前兆異常提取方法。這些方法，對于我們今后在處理、分析長趨勢漂移上升（下降）型地震前兆觀測數據時，具有極其重要的作用，為今后利用水溫數據進行預報地震探索開辟了一條新路徑。

（3）在地震前兆異常研究中盡管有了一定的認識，但這種認識仍然還是比較零散的。此外，不同震級水平的地震，其震兆特征差異性可能很大。由于地震屬于小概率事件，特大地震發生的概率特別小，而現在研究的地震前兆多是4～6級地震的前兆，對于8級以上地震前兆的研究非常少，所以在面臨地震前兆異常判斷時，可能會不知不覺中以“小”代“大”。由于汶川地震是近幾十年來發生在甘東南毗鄰區最大的一次地震，而這4井水溫觀測時間又不長，汶川地震前，對此4口井的變化難免認識不足，因此，對這些水溫井前兆異常進行再分析、再認識非常有必要。

（作者電子信箱，蒲小武：wdpuxw＠163.com）

［1］王勇，張國紅，代建華，等.昌平地震臺地溫前兆特征分析.西北地震學報，2000，22（3）：236－240

［2］胡敦寬，李淑芳，劉永銘.白家疃觀測井的地熱干擾排除及機理探討.地震學報，1996，18（1）：128－131

［3］車用太，魚金子，劉春國.我國地震地下水溫度動態觀測與研究.水文地質工程地質，1996（4）：34－37

［4］車用太，劉喜蘭，姚寶樹，等.首都圈地區井水溫度的動態類型及其成因分析.地震地質，2003，25（3）：404－419

［5］車用太，魚金子.塔院井水溫的大幅度奇異變化∥中國地震局監測預報司.地下流體典型異常的調查與研究.北京：氣象出版社，2004：414－418

［6］谷元珠，車用太，魚金子，等.塔院井水溫微動態研究.地震，2003，23（1）：102－108

［7］魚金子，車用太，劉伍洲.井水溫度微動態形成的水動力學機制研究.地震，1997，17（4）：389－396

［8］付子忠.地熱動態觀測與地熱前兆.地殼構造與地殼應力，1988，1（1）：1－5

［9］車用太，魚金子，劉五洲，等.張北—尚義地震的地下流體異常場及其成因分析.地震學報，1999，21（2）：194－201

［10］劉耀煒，孫小龍，王世芹，等.井孔水溫異常與2007年寧洱6.4級地震關系分析.地震研究，2008，31（4）：347－353

［11］趙剛，何案華，王軍.地熱正常動態特征的研究.地震，2009，29（3）：109－116

［12］趙剛，何案華，馬文娟，等.不同動態背景的地熱對比觀測研究.地震學報，2011，33（1）：51－61

［13］葉秀薇.粵閩地區地下流體從屬函數異常與地震關系的初步研究.防災減災工程學報，2004，24（2）：195－201

Characteristics analysis on the abnormal changes of the water temperature before and after the Wenchuan MS8.0form 4 wells located in southeast of Gansu

Pu Xiaowu，Mei Donglin，Chen Yanping，Ye Yuanyuan，Wang Junyan，Xu Kejuan
（Longnan Seismic Station，Earthquake Administration of Gansu Province，Gansu Longnan 746000，China）

This paper systematically analyzes water temperature observation data since 2007from 4wells Chenxian，Fanba，Qingshui and Liangshui which are located in southeast of Gansu.It is found in 2008，before and after the May 12wenchuan MS8.0，water temperature of 4wells appeared different degree of precursory anomaly changes，but most abnormal precursory informations submerged in the background of the rising trend.In order to extract precursory anomaly information，methods such as the difference method，the variable rate method，the subordinate function method and trend rate method which are often used in fluid data processing are adopted to exact the abnormal precursory information，and analyze the changes of water temperature anomalies.

water temperature anomaly；difference；variable rate；subordinate function；trend rate

P315.72＋3；

： A；

10.3969／j.issn.0235－4975.2014.02.004

2013－09－09；

2013－09－22。

中國地震局蘭州地震研究所地震科技發展基金野外站基金（2013096）。