蘆山MS7．0地震前遠、近場形變時空演化特征研究*

牛安福 顧國華 曹景鵬1， 張凌空閆 偉 趙 靜 吉 平

1）中國北京 100045 中國地震臺網(wǎng)中心

2）中國北京 100036 中國地震局地震預測研究所

引言

大地震的短期前兆現(xiàn)象一直是國內外學者關注與爭論的焦點之一．基于擴容擴散和擴容引起的摩擦滑動模型及地震破裂成核過程模型（Nur，1972；Mjachkinetal，1975；Rummeletal，1978；Das，Scholz，1981；Dubrovskiy，Sergeev，2006），地傾斜應變觀測等被認為是監(jiān)測地殼形變和捕捉地震前兆的重要手段．而實際觀測結果表明，震級小于4．3級的小地震前地傾斜應變信號與周圍地震具有一定的相關性（Johnston，Mortensen，1974；McHugh，Johnston，1978），但異常臺站至震中距離很難達到0km或0m（Rikitake，1975），因此，Takemoto（1991）提出了地震發(fā)生的最小前兆距離的概念；對于震級大于6．0級地震，臨震前盡管觀測到一些突出的短期異常，但通常遠離震中區(qū)域（Mogi，1985；牛安福等，2009，2012）；而靠近震源區(qū)域，由于臨震前變形不顯著，同樣引起對以監(jiān)測斷層預滑動預報地震思路的質疑（Takemoto，1991；Johnstonetal，1990；Lindeetal，1992；Bilham，2005；Amoruso，Crescentini，2010；牛安福等，2012）．

臨近蘆山MS7．0地震前的形變前兆分布，為認識該地震孕育過程及其伴隨的前兆現(xiàn)象等提供了不可多得的信息．2013年4月20日發(fā)生的蘆山MS7．0地震位于龍門山斷裂帶南段，其破裂特征與汶川MS8．0地震相似，屬于逆沖型地震．該震震前沿鮮水河斷裂帶一些跨斷層基線觀測到顯著的加速轉折變化，顯示因巴顏喀拉地塊的快速東向運動，打破了鮮水河斷裂長期左旋運動的格局，是突出的場兆變化；沿安寧河、則木河斷裂帶個別場地的跨斷層水準基線呈巨幅變化；沿龍門山斷裂帶一些水準觀測在汶川MS8．0地震后呈持續(xù)調整變化；而鮮水河、龍門山、安寧河3條主要斷裂圍成的三叉口地區(qū)，地傾斜、應變、GPS、重力等連續(xù)形變變化不大，斷層蠕變及斷層水準變化較小，是形變變化最為薄弱的區(qū)域．

本文針對蘆山MS7．0地震臨震前多種形變觀測結果進行綜合分析，提出了震前出現(xiàn)的形變異常“空區(qū)”現(xiàn)象，并對該現(xiàn)象給出了可能的解釋．

1 蘆山MS7．0地震前突出形變異常分布

自20世紀70年代末開始，四川省地震局沿龍門山斷裂帶、鮮水河斷裂帶、安寧河斷裂帶、則木河斷裂帶陸續(xù)開展了跨斷層水準基線測量和斷層蠕變測量，斷層水準基線觀測周期為1—3個月．本文使用的觀測數(shù)據(jù)截止時間為蘆山地震前最后一期，即2013年3月份．該地區(qū)自2003年開始多數(shù)地傾斜和地應變臺站等實現(xiàn)了數(shù)字化觀測，采用的資料截止時間為震前1天．

1．1 鮮水河斷裂帶異常分布

蘆山MS7．0地震前鮮水河斷裂帶較為顯著的異常變化主要是侏倭、格簍、虛墟和龍燈壩跨斷層基線．這些基線長期變化趨勢表現(xiàn)為巴彥喀拉塊體南邊界以左旋運動為主．而2013年1月、3月份的觀測結果卻給出了顯著的轉折變化，短期的活動以右旋為主（圖1）．

圖1 蘆山MS7．0地震前沿鮮水河斷裂帶跨斷層基線變化．LAB和LAC分別表示AB邊和AC邊的邊長Fig．1 The changes in cross fault baseline lengths along Xianshuihe fault prior to the Lushan earthquake，where LABand LACare the length of ABand AClines，respectively

該斷裂帶跨斷層水準變化與斷層基線相比不顯著．沿鮮水河斷裂帶還有蝦拉沱、恰叫、溝普、老乾寧和龍燈壩斷層蠕變測量，除溝普外，多以趨勢異常為主．

1）侏倭基線．該場地距離蘆山MS7．0地震震中280km．該斷層觀測場地自1980年開始觀測，AB邊長168m，AC邊長288m左右．早期每2—3個月測量一期，汶川MS8．0地震后每月測量一期．1982年以來，該處斷層的水平運動呈左旋張性活動．AB邊以平均0．8mm／a左右速率拉張，臨震前有短時異常擾動；AC邊以平均1．0mm／a速率縮短，汶川MS8．0地震后，AC邊變化較小，2013年1月前后加速壓縮，震前最大壓縮量達2．8mm，異常持續(xù)時間約3—4個月．

2）格簍基線．該場地距離蘆山MS7．0地震震中270km．該斷層觀測場地自1982年開始觀測，AB邊長72m，AC邊長96m左右．早期每2—3個月測量一期，2004年起每月測量一期．1982年以來，該處斷層的水平運動呈左旋張性活動，平均活動量在0．1mm／a左右．2001年昆侖山口西MS8．1地震后，年變幅度增大，但年變速率未有大的改變．2012年5月起基線趨勢出現(xiàn)轉折，2013年1月有明顯加速轉向，AB測邊出現(xiàn)大幅度縮短，變化幅度為3．5mm，異常持續(xù)時間近4個月．

3）虛墟水準、基線．該斷層觀測場地距離蘆山MS7．0地震震中260km．該場地基線觀測自1980年開始，AB邊和AC邊基線分別長215m和240m左右，每1—3個月觀測一期．基線LAB和LAC受2001年昆侖山口西MS8．1地震、2004年印尼蘇門答臘MS9．0地震、以及2008年汶川MS8．0地震影響，基線變化趨勢發(fā)生改變．汶川地震后，斷層左旋張性活動，2011年底至2012年10月拉張的趨勢有所減弱，2013年1月AB測邊出現(xiàn)大幅度壓性變化，加速持續(xù)時間近4個月，變化幅度為2．96mm．

4）龍燈壩基線．該斷層觀測場地距離蘆山MS7．0地震震中160km．該場地基線觀測自1985年開始，AB邊和AC邊基線分別長48m和72m左右，每1—3個月觀測一期．龍燈壩場地的水平運動歷來較弱，1991—2002年兩測邊均以0．04mm／a左右速率變化，2003—2006年以0．19mm／a左右速率拉張．2011年11月至2013年1月，AB測邊出現(xiàn)持續(xù)13個月的下降，最大幅度為2．03mm．

5）其它場地基線．老乾寧斷層觀測場地距離蘆山MS7．0地震震中140km．該場地基線觀測自1979年開始，每1—3個月觀測一期．2012年以來兩測邊觀測曲線波動變化，臨震前無明顯的加速斷層活動．折多塘斷層基線（距蘆山地震震中110km）呈現(xiàn)較穩(wěn)定的趨勢，無顯著的短期加速異常變化．

6）溝普水平蠕變．該場地距離蘆山MS7．0地震震中200km．溝普蠕變觀測1—3測邊在2008年汶川MS8．0地震后趨勢發(fā)生轉折，在2010年玉樹地震后出現(xiàn)較大幅度起伏．異常核實報告顯示，未發(fā)現(xiàn)儀器故障、人為干擾等．2011年5—9月持續(xù)上升，幅度為0．59mm，10月轉折下降．2012年1—10月繼續(xù)大幅度上升，幅度為1．1mm（圖2）．

從鮮水河斷層基線及蠕變觀測結果可知，顯著的短期異常變化主要分布在溝普以西地段，而鮮水河東南段，龍燈壩基線、老乾寧及折多塘基線則沒有明顯的短期加速異常．

圖2 溝普水平蠕變1—3測邊觀測曲線Fig．2 The horizontal fault creep of 1—3 lines at Goupu site

1．2 沿龍門山斷裂帶形變異常分布

沿龍門山斷裂帶斷層活動主要表現(xiàn)為汶川MS8．0地震后的調整．周圍較近的觀測場地包括寶興、雙河和灌縣．

1）寶興水準．該場地水準觀測開始于汶川MS8．0地震后，距離蘆山MS7．0地震震中15 km．從2010年以來的觀測結果看，水準變化幅度較小，震前3年變化幅度僅0．6mm．震前有小幅擾動，但也在動態(tài)范圍內（圖3）．

圖3 蘆山MS7．0地震前沿龍門山斷裂帶跨斷層水準變化圖中HXX1，H3—4，H2—1分別表示XX1邊、3—4邊和2—1邊測線高程變化Fig．3 The changes in cross fault leveling along the Longmenshan fault before the Lushan earthquake，where HXX1，H3—4，and H2—1denote the leveling of XX1，3—4，and 2—1lines，respectively

2）雙河水準．雙河場地跨越龍門山前山斷裂帶之安縣—灌縣斷層中段的雙河斷層，距離蘆山MS7．0地震震中90km．1991年以來雙河短水準的變化主要在2mm范圍內波動．汶川MS8．0特大地震前沒有監(jiān)測到明顯的前兆異常，但同震變化非常顯著，上盤相對下盤下降7．76mm．汶川MS8．0地震后，斷層一直處于壓性調整活動階段，年動態(tài)變化在1．5mm范圍內（圖3）．

3）灌縣水準．灌縣水準場地距離蘆山MS7．0地震100km．2000年以來具有較穩(wěn)定的年周期變化特征．受汶川MS8．0地震的影響，產(chǎn)生顯著的同震變形．之后至蘆山地震前，呈現(xiàn)持續(xù)下降，幅度達3mm，是龍門山斷裂帶上水準變化速率相對較大的場地之一（圖3）．

1．3 沿安寧河、則木河斷裂帶形變異常

蘆山MS7．0地震前沿安寧河、則木河斷裂帶突出的形變異常主要有西昌、爾烏水準和湯家坪水準基線．寧南、冕寧、安順場等場地水準存在一定的變化，但不顯著；紫馬跨斷層蠕變在汶川MS8．0地震后出現(xiàn)年變化幅度增大現(xiàn)象，震前年變化特征較為穩(wěn)定．

1）西昌水準．西昌短水準測量場地位于西昌市李金堡鎮(zhèn)，距蘆山MS7．0地震震中270km．該場地布設了A-B-C-A短水準監(jiān)測環(huán)線．其中AB、BC為跨斷層斜交邊．資料顯示，1988—2000年平均速率為0．27mm／a左右，斷層略顯張性活動．2004年開始對該場地進行每年6期的加密觀測，2004—2009年連續(xù)6年的觀測顯示該處斷層無明顯變化，年變化比較清晰．2010年4—6月的觀測顯示，該場地兩測線同時出現(xiàn)明顯的變化，AB測線下降了3．89mm，AC測線上升了4．46mm．2010年7月—2011年3月有所轉折．2011年5月—2012年5月，兩測線又同時出現(xiàn)明顯的變化，AB測線下降了10．6mm，AC測線上升了10．5mm．2012年7—9月的觀測資料顯示曲線有所轉折，該時段內在川滇交界先后發(fā)生了2012年6月24日的鹽源MS5．7地震和2012年9月7日的彝良MS5．6地震，震中距分別為150km和180km．兩次地震發(fā)生在該水準轉折變化過程中（圖4）．2013年1月觀測結果顯示再次發(fā)生轉折，持續(xù)時間4個月左右．

圖4 蘆山MS7．0地震前沿安寧河、則木河斷裂帶跨斷層水準、基線變化HAB，HAC，H1—3分別表示相應測線的水準變化，LAB表示AB測線長度Fig．4 The changes in cross fault leveling and baseline lengths along the Anninghe and Zemuhe faults before the Lushan earthquake．HAB，HAC，H1—3separately represent corresponding leveling changes for the lines AB，AC，and 1—3，and LABis the length of ABline

2）爾烏水準．爾烏短水準在2008年之前的變化顯示此處斷層垂直活動較弱．汶川MS8．0大震后打破了往年同期的上升變化，出現(xiàn)近3mm幅度的異常擾動．2010年玉樹地震后，再次發(fā)生轉折，上升幅度達2．34mm后恢復．2012年開始異常再次出現(xiàn)，上升幅度為2．23mm，臨震前出現(xiàn)轉折（圖4）．

該場地基線長期變化具有10年左右周期．2003—2007年出現(xiàn)較大幅度的異常擾動，2008年汶川MS8．0地震后，趨勢出現(xiàn)轉折．2011年開始加速伸長，意味著發(fā)生了較大幅度的左旋運動，年最大幅度達20mm，鹽源MS5．7和彝良MS5．6兩次中等地震均發(fā)生在基線達到最大值前后，臨近蘆山MS7．0地震前，基線出現(xiàn)轉折（圖4）．

2 蘆山MS7．0地震前近場形變變化特征

應用地震震級估算此次地震破裂尺度約30km．將2倍破裂尺度內定義為近場，而2倍尺度外定義為遠場．近場的范圍為震中附近60km，沿地震破裂方向可接近90km．

蘆山MS7．0地震近震源區(qū)域的形變觀測點包括寶興跨斷層場地，雅安和姑咱等連續(xù)形變觀測臺站．

蘆山MS7．0地震最近的跨斷層觀測場地是寶興場地，震中距離15km．該場地跨斷層水準結果表明，其變化幅度在研究區(qū)域內是最小的，震前3年變化幅度也只有0．6mm，且震前擾動變化未超出2倍動態(tài)范圍（圖3）．

1）雅安臺地傾斜變化．雅安地震臺處于龍門山斷裂帶西南段，位于蘆山MS7．0地震極震中區(qū)域．雅安地震臺地形變觀測儀器為JB型金屬水平擺傾斜儀，2003年完成數(shù)字化改造，儀器工作周期14．2s±0．1s．儀器洞室東西走向，條石拱砌，進深7m，頂層及兩側覆蓋厚度＞20m，儀器墩高60cm，巖性為砂泥巖．儀器洞室年溫度變幅≤1．6℃，日溫度變幅≤0．2℃，相對濕度＞95%．日照和洞內潮濕對觀測會有些干擾，但年周期變化較為穩(wěn)定．

該臺傾斜潮汐記錄數(shù)據(jù)分析結果表明，除個別計算粗差外，M2波潮汐因子及相位較為穩(wěn)定（圖5）．2012年6月24日鹽源MS5．7地震前變化不大．

平時室內要保持空氣新鮮流通，相對濕度60%左右；保持寶寶的呼吸道暢通，根據(jù)病情給予化痰，及時清除分泌物，必要時給予支氣管解痙治療和氧療。

應用傅里葉周期分析方法，可提取相應的年周期成分．其結果表明，該臺地傾斜矢量模在震前沒有出現(xiàn)顯著的異常變化，但年周期成分自2008年開始則出現(xiàn)逐年下降特征（圖6b）．

圖5 蘆山MS7．0地震前雅安臺地傾斜M2波潮汐因子（a）、相位（b）及相應的誤差Fig．5 The changes of tilt tidal factor（a）and phase（b）at Ya′an station before the Lushan earthquake where the bars denote the monthly errors

圖6 蘆山MS7．0地震前雅安臺傾斜矢量模（a）及年周期變化（b）Fig．6 The tilt vector modulus（a）and the annual period changes（b）at Ya′an station before the Lushan earthquake

2）姑咱臺地形變變化．姑咱臺處于北西向的鮮水河斷裂帶、北東向的龍門山斷裂帶和南北向的安寧河斷裂帶復合部位靠北的地段，距離蘆山MS7．0地震震中90km．洞室的巖性為前震旦紀似斑狀黑云母花崗閃長巖，形變洞進深90m，覆蓋厚度大于70m，洞溫18．5℃左右，日溫變幅小于0．03℃，年溫變幅小于0．5℃，相對濕度大于95%；洞室內裝有DZW型重力儀、SQ型水平擺、DSQ型水管傾斜儀、SSY型伸縮儀等，該臺在2007年還安裝了YRY分量鉆孔應變儀．

由該臺多套傾斜儀、應變儀記錄到數(shù)據(jù)計算得到的M2波潮汐因子及相位無顯著的短期異常變化．該臺定點重力結果中的年周期成分在2004—2005年、2010—2011年擾動幅度較大，然后逐年變小，之后分別發(fā)生了汶川MS8．0地震和蘆山MS7．0地震（圖7）．

蘆山MS7．0地震前兩個臺站資料計算得到的年周期成分衰減現(xiàn)象，反映了臨震前孕震區(qū)域能量耗散存在逐年減小的特征，這是開放系統(tǒng)所具有的減熵特征．1850年，克勞修斯提出熱力學第二定律，即熱量可以自發(fā)地從高溫物體轉移到低溫物體，但不能自發(fā)地從低溫物體轉移到高溫物體．這一定律指出了開放的熱力學系統(tǒng)自發(fā)變化的必然方向，并被應用在通訊、生物等諸多領域．因此，臨震前記錄到的各種前兆變化，都是半無限空間內、弱約束條件下孕震系統(tǒng)演化的結果，這也是許多前兆現(xiàn)象難以使用一些有限尺度的試驗結果或一些物理模型來解釋的原因之一．20世紀末，我國地震學者也曾開展了地震孕育過程中的熵減現(xiàn)象研究，取得了一些有意義的成果（朱傳鎮(zhèn)，1989；朱令人等，1990；王海濤，楊馬陵，1990；牛安福，1991）．本文給出的蘆山MS7．0地震前年變幅度減少現(xiàn)象是地震孕育過程中熵減的又一體現(xiàn)，對借助連續(xù)變形觀測進行近震源區(qū)域判定提供了一種重要的判據(jù)．

圖7 蘆山MS7．0地震前姑咱臺重力（a）及年周期成分變化（b）Fig．7 The changes in gravity（a）and the annual periodic component（b）at Guzan station before Lushan earthquake

3 蘆山MS7．0地震前形變異常分布特征討論

蘆山MS7．0地震前出現(xiàn)的一些形變突出異常是否與該地震相關，以及如何解釋近震源區(qū)域出現(xiàn)的形變變化等均是地震前兆研究的基本問題．以往通常按照地域來建立異常與地震的相關性，但發(fā)現(xiàn)異常特征與地震三要素的關系是復雜的．

圖8 遠場異常與地震相關性判別的時間依賴準則Fig．8 Sketch map of the time dependent criterion for constructing the correlativity between far-filed deformation anomalies and earthquakes

對于大地震而言，與其相關的前兆變化可能出現(xiàn)在很遠的區(qū)域．例如，1944年日本Tonankai 8級地震前觀測到顯著的傾斜異常，Kagegawa測點距離該地震震中約300km（Mogi，1985）．為解決遠場形變異常與地震的相關性，作者提出了時間依賴判別準則，即異常越接近地震發(fā)生時刻，其與地震的相關性越高（圖8）．地震波及同震變形與地震震級和地點具有較確定的聯(lián)系，擁有最高的相關系數(shù)為1．0．而中期、短期及短臨異常都低于最高值1．0，且以短臨異常最高，短期異常次之，這體現(xiàn)了地震預測的“最后一分鐘”現(xiàn)象（牛安福，2005）．地震波與同震變形場等與地震參數(shù)都可建立一定的模型聯(lián)系，作為震前的前兆也期待有這種結果．作者曾開展過相關的研究與試驗工作（牛安福，2003）．

考慮到蘆山MS7．0地震前鮮水河斷裂帶出現(xiàn)的斷層基線短期加速變化及觀測周期，估計其異常持續(xù)時間為90—120天；安寧河、則木河斷層水準在巨幅異常結束至蘆山地震前，持續(xù)時間為180—240天．依據(jù)形變短期加速異常持續(xù)時間與地震震級的聯(lián)系（牛安福，2003），估計其相應的震級范圍分別為（6．4—6．6）±0．4和（6．8—7．0）±0．4，因此可認為蘆山地震震前發(fā)生的這些異常變化與蘆山地震具有一定的相關性．此處所指相關只是從時間依賴的角度而言，并非意味著該地區(qū)今后不會發(fā)生較大地震．

對蘆山MS7．0地震前各種形變資料的進一步分析發(fā)現(xiàn)，震前形變異常形成的“空區(qū)”是較明顯的，該空區(qū)屬于變形幅度小或異常變化不顯著的區(qū)域（圖9）．形變異?！翱諈^(qū)”是臨近強地震前出現(xiàn)的一種重要的前兆分布現(xiàn)象．周碩愚等（1997）利用精密水準復測結果和跨斷層資料，研究了唐山地震前震源區(qū)域形變場的演化，并驗證了形變空區(qū)與地震活動區(qū)及斷層閉鎖區(qū)域間的聯(lián)系．牛安福等（1999）利用F檢驗方法，研究了1996年2月3日麗江MS7．0地震前后異常強度的分布，提出了異常強度空區(qū)，并研究了形變“空區(qū)”內外地震活動性的差異．蘆山地震前觀測到的變形空區(qū)明顯受到臺站布局及觀測手段多樣性的影響，這里展示的只是一個大致的分布．

圖9 蘆山MS7．0地震前的形變異常分布及形變異?！翱諈^(qū)”紅色矩形表示蘆山地震前形變變化較大或異常顯著的場地，虛線橢圓表示形變異?！翱諈^(qū)”，紅色實心圓表示2012年以來研究區(qū)域內發(fā)生的MS5．5以上地震，黑色三角形、矩形分別表示連續(xù)形變觀測臺站和斷層觀測場地，紅色線、綠色線分別表示斷層及河流分布Fig．9 The deformation anomaly distribution showing the deformation weak area（denoted by dashed ellipse）Red rectangles denote the sites with significant deformation changes or anomalies，the red circles denote the MS≥5．5earthquakes since 2012，respectively；black triangles and rectangles denote continuous deformation observation stations and cross fault observation sites，the red and blue lines show the faults and rivers

圖10 蘆山MS7．0地震震中區(qū)域水平位移年速率分布虛線橢圓表示水平運動速率較低的區(qū)域Fig．10 Annual velocities of horizontal displacements in the epicentral area of Lushan earthquake The dashed ellipse delineates the area with lower horizontal movement velocities

GPS資料也驗證了蘆山MS7．0地震前形變空區(qū)的存在．利用GPS陸態(tài)網(wǎng)絡和中國地震局地震預測研究所汶川地震科考GPS站2010年6月30日—2013年2月20日數(shù)據(jù)，以東部穩(wěn)定基準計算獲得了各站的水平位移年速率，依據(jù)其速率大小也可勾畫出運動速率較小的臺站所在的區(qū)域（圖10）．該區(qū)域內，靠近震中最近的LS05和LS06臺站距離震中17km左右．其水平位移年速率較周圍臺站是最小的，垂直位移速率也具有相近的特征．

比較圖9與圖10，盡管給出的兩個形變“空區(qū)”在形態(tài)和范圍上有所差別，但基本上都反映出近震源區(qū)域是變形速率最小、異常最少的區(qū)域．

4 結論

對蘆山MS7．0地震前多種形變觀測數(shù)據(jù)的分析表明，強地震前近場與遠場形變變化的特征具有顯著的差別．沿鮮水河斷裂帶一些跨斷層基線自2013年1月開始的短期加速變化，及沿安寧河、則木河斷裂帶個別場地的跨斷層水準基線自2010年4月出現(xiàn)的巨幅異常等與蘆山地震在異常持續(xù)時間和發(fā)震構造上具有一定的聯(lián)系．這些突出的形變異常主要分布在蘆山地震的外圍區(qū)域．在靠近震中附近的鮮水河、龍門山和安寧河3條主要斷裂圍成的三叉口地區(qū)，連續(xù)形變（傾斜、應變、重力）及斷層水準和蠕變觀測臨震前均未有顯著的異常變化，GPS水平、垂直位移年速率最小，是形變變化或形變異常分布的“空區(qū)”．其中，近震源區(qū)域附近臺站的傾斜、重力年變化幅度震前2—3年不是增大，而是減?。?/p>

蘆山MS7．0地震前各種形變資料的分析結果與構造活動引起的變形場也有較大的差別，甚至也有悖于常規(guī)意義上對“以場求源”，即在近震源區(qū)域捕捉失穩(wěn)變化預報地震思路的理解，但這種結果可能體現(xiàn)了強地震應變積累與地震釋放之間矛盾的交替（吳翼麟，1990），即通過大變形釋放地震能量，通過小變形積累起大的地震應變．

蘆山MS7．0地震前觀測到的形變前兆現(xiàn)象特征與汶川MS8．0地震等震前的前兆現(xiàn)象較為接近（Johnstonetal，1990；Lindeetal，1992；Bilham，2005；Amoruso，Crescentini，2010；牛安福等，2012）．因此，蘆山MS7．0地震前近震源區(qū)域觀測到的連續(xù)形變觀測前兆為認識與強地震孕育有關的變形過程及開展地震預報研究等提供了新的思路．

牛安福．1991．陽高6．1級地震前傾斜固體潮振幅因子熵的變化［J］．地殼形變與地震，11（2）：1-7．

牛安福．2003．地傾斜變化的突變性及與地震關系的研究［J］．地震學報，25（4）：441-445．

牛安福．2005．地震短期、短臨預測面臨的“最后一分鐘”［J］．國際地震動態(tài)，（2）：29-31．

牛安福，張凌空，閆偉，吉平．2012．汶川地震近震源區(qū)地形變短期前兆現(xiàn)象的解析［J］．地震，32（2）：52-63．

牛安福，張凌空，閆偉，賈曉東，李曉帆．2009．汶川地震前南北地震帶中北段地形變變化特征的研究［J］．地震，29（1）：100-107．

牛安福，張雁濱，柯麗君，劉祖蔭，鄭文衡．1999．地震前地殼形變異常分布的非均勻性特征［J］．地震，19（2）：149-154．

王海濤，楊馬陵．1990．地震前兆系統(tǒng)信息熵的初步研究［J］．中國地震，6（2）：10-14．

吳翼麟．1990．定點形變前兆預報地震的觀測技術與分析方法［J］．地震，（5）：33-45．

周碩愚，施順英，帥平．1997．唐山地震前后地殼形變場的時空分布、演化特征與機理研究［J］．地震學報，19（6）：559-565．

朱傳鎮(zhèn)．1989．應用熵的原理推進地震活動定量化研究［J］．地震，12（1）：75-76．

朱令人，王海濤，白超英，周仕勇，龔宇清．1990．瀾滄—耿馬7．6級地震前中小地震活動信息熵和分維的研究［J］．內陸地震，4（4）：310-317．

Amoruso A，Crescentini L．2010．Limits on earthquake nucleation and other pre-seismic phenomena from continuous strain in the near field of the 2009L′Aquila earthquake［J］．GeophysResLett，37：L10307，doi：10．1029／2010GL043308，2010．

Bilham R．2005．Co-seismic strain and the transition to surface afterslip recorded by creep-meters near the 2004Parkfield epicenter［J］．SeismResLett，76（1）：49-57．

Das S，Scholz C H．1981．Theory of time-dependent rupture in the Earth［J］．JGeophysRes，86（B7）：6039-6051．

Dubrovskiy V A，Sergeev V N．2006．Short-and medium-term earthquake precursors as evidence of the sliding instability along faults［J］．PhysicsoftheSolidEarth，42（10）：802-808．

Johnston M J S，Linde A T，Gladwin M T．1990．Near-field high precision strain prior to the October 18，1989Loma PrietaML7．1earthquake［J］．GeophysResLett，17：1777-1780．

Johnston M J S，Mortensen C E．1974．Tilt precursors before earthquakes on the San Andreas fault，California［J］．Science，186：1031-1034．

Linde T，Gladwin M T，Johnston M J S．1992．The Loma Prieta earthquake，1989and earth strain tidal amplitudes：An unsuccessful search for associated changes［J］．GeophysResLett，19（3）：317-320．

Mjachkin V I，Brace W F，Sobolev G A，Dieterich J H．1975．Two models for earthquake forerunners［J］．PureAppl Geophys，113（1）：169-181．

McHugh S，Johnston M J S．1978．Dislocation modeling of creep-related tilt changes［J］．BullSeismolSocAm，68：155-168．

Mogi K．1985．Temporal variation of crustal deformation during the days preceding a thrust-type great earthquake：The 1944Tonankai earthquake of magnitude 8．1，Japan［J］．PureApplGeophys，122（6）：765-780．

Nur A．1972．Dilatancy，pore fluids and premonitory variations oftS／tPtravel times［J］．BullSeismolSocAm，62（5）：1217-1222．

Rikitake T．1975．Dilatancy model and empirical formulas for an earthquake area［J］．PureApplGeophys，113（1）：141-146．

Rummel F，Alheid H J，F(xiàn)rohn C．1978．Dilatancy and fracture induced velocity changes in rock and their relation to frictional sliding［J］．PureApplGeophys，116：743-764．

Takemoto S．1991．Some problems on detection of earthquake precursors by means of continuous monitoring of crustal strains and tilts［J］．JGeophysRes，96（B6）：10377-10390．