文安地震前后地殼形變的時空演化特征

0構造及地震活動與異常2006年7月4日11時56分，河北省文安縣發生了m5.1級地震。首都圈的北京、天津、河北、山西、山東部分地區均有不同程度的震感。研究發現文安地震前重力測網內出現了以上升為主的異常變化,在文安地震震中,短期內形成了重力高梯度帶。同時有專家研究發現體應變、跨斷層位移測點群體趨勢轉折,鉆孔應力應變測點先后出現了中短期異常。上述文獻揭示了地震前前兆信號存在的可能性,因此本文綜合利用洞體應變和鉆孔應變資料,采取不同的方法研究文安地震前后地殼形變的變化特征。洞體應變和鉆孔應變觀測作為重要的地震前兆監測手段,在地震監測預報研究中發揮著重要的作用,尤其在地震短臨預報研究中取得了許多研究成果,為地震預報研究積累了豐富的觀測資料。張晶等綜合分析認為,文安地震前局部應變場的異常變化在200km范圍內,異常持續3個月左右,池順良觀測到了汶川和蘆山強震孕震過程中四川姑咱地震臺的潮汐畸變應變異常現象。周文斌認為在俄羅斯東海岸6.2級地震前敦化臺鉆孔應變資料記錄到了顯著的短期和臨陣異常。山丹臺體應變在汶川8.0級地震前46h、8h、37min時分別記錄到了明顯的短臨階躍異常。本文綜合利用2種應變觀測資料的時間序列法和圖像法,研究文安地震前后應變演化特征及與地震和構造的關系,并進行定性分析。1易縣地震臺—臺站概況在首都圈地區,距離文安震中200km范圍內運行的應變觀測臺站共計10個(圖1),但2006年前正常運行并獲得數據的只有4個臺站,因此選取這4個臺的應變觀測資料進行分析。為對比距離震中較遠的形變臺站的變化情況,本文特別分析了距震中約245km的張家口臺站(表1),這些應變連續觀測站的布設使首都圈地區的應變觀測達到了一定的密度。圖1顯示該次地震發生在牛東斷層和高陽-博野斷裂帶附近。離震中最近的臺站為天堂河地震臺和易縣地震臺。大興天堂河臺位于北京市大興區天堂河醫院院內,距離震中約89km,是首都圈Sacks體應變綜合觀測臺網的一個子站,南苑-通縣斷裂在附近通過。井下探頭使用的是美國華盛頓卡內基研究所研制的Sacks-Evertson鉆孔應變儀,能記錄到清晰的固體潮變化,穩定性好。2001年9月完成數字化改造,實現分鐘值采樣和數據自動傳輸。易縣地震臺位于保定市西北部,在易縣縣城西20km處,距離震中約95km,該臺為形變綜合臺,其中,SS-Y型伸縮儀和TJ-II型鉆孔體應變儀自2001年起正式運行以來,儀器工作正常,產出觀測資料連續率良好,穩定可靠。位于河北省易縣縣城西20km處,太行山山前大斷裂和西部紫荊關大斷裂之間,其附近有次級斷裂構造(大紅門斷層)。昌平地震臺位于北京市昌平區十三陵風景區臥虎山腳下,距離震中150km,構造上位于北東向南口-山前活動斷裂和北西向南口-孫河活動斷裂交匯部位以東7.5km處,屬于應力變化敏感區。其中,RZB-1分量應變儀在2008年6月以前儀器各分量運行狀態基本正常,TJ-II型體應變儀于2005年3月18日安裝并啟用。小辛莊觀測山洞為部隊戰備山洞,位于薊縣城東15km,薊縣臺臺站東約4.5km,距離震中約19km,山洞處于陰山緯向構造體系東部的東西向壓性斷裂帶的許家臺-薊縣-遵化山前斷裂區,與新華夏系華北沉降帶東北邊緣相接。1983年底小辛莊山洞改造建設竣工并投入使用,架設了SSY-Ⅱ石英伸縮儀。張家口地震臺地處京西北、晉冀蒙交界地區,位于陰山東西向構造帶與山西地震帶北端、張家口-渤海地震帶西北端的交匯部位,距離震中最遠,約為245km。該臺的TJ-Ⅰ型鉆孔體應變儀和SS-Y型伸縮儀工作正常,產出觀測資料連續率良好,穩定可靠。2觀測數據的分析本文對上述5個臺站應變觀測資料的日均值和分鐘值進行了初步分析,并對易縣觀測站和張家口觀測站的伸縮應變參數進行計算分析。2.1構造及力學分析從表1所示的5個臺站8套應變儀器來看,200km范圍內的多數臺站存在同震效應。圖2顯示了4個體應變臺的臨震和震后的分鐘值變化曲線,可以看出在震前4天昌平、易縣、張家口臺的應變幅度在7月1—4日均呈下降趨勢,在震后的7月5日達到最低,從7月6日開始應變幅度呈上升趨勢;僅天堂河臺的應變幅度在7月1—4日呈上升趨勢,在7月6日達到最高。同時,4個臺的峰谷幅度差值逐漸降低,在地震當天達到最低,震后峰谷差值逐漸增加。圖3為易縣(a)、張家口(b)、小辛莊(c)的洞體應變原始分鐘值變化曲線圖,從圖上可以清晰地看出,3臺的各分量曲線在1—10日出現上升趨勢,其峰谷差值也呈現上升的趨勢,而地震恰恰發生在峰谷差值最低點的7月4日。圖4顯示了昌平臺鉆孔應變的NS分量應變幅度在7月6日出現轉折,由升高轉為下降,而NE和NW分量應變幅度在7月6日變化速率加快,但總體變化趨勢一致。圖5為小辛莊臺2006年各分量的日均值曲線變化圖,可以看出自2006年3月底至6月底期間,先是在3月底至4月初出現輕微的膨脹變化,然后出現急劇的壓性下降狀態,直到地震前出現階躍變化。圖6為易縣臺伸縮儀三分量觀測曲線,NS向觀測在2001—2005年呈下降的壓性變化,2006年初由下降轉為上升變化,在上升和下降過程中發生了文安5.1級地震,2007年恢復壓性下降狀態,2010年初再次轉為向上的膨脹變化,2012年恢復。EW向觀測在2001—2004年呈下降的壓性變化,2005年變化趨勢逐漸穩定,2011年再次轉為下降的壓性變化。NE向于2005年下半年轉為向下的壓性變化,2006年出現短暫的平緩變化,2007年恢復壓性下降狀態,2009年以后曲線呈相對緩慢的壓性下降狀態,2011年出現轉折。上述臺站均屬華北北部、首都圈地區,在文安地震前,從短臨的分鐘值曲線和長期的日均值曲線來看,鉆孔應變、體應變和洞體應變在臨震前出現了大范圍同步轉折變化,其現象表明,在區域應力、震源應力的共同作用下,大范圍地殼不同層位巖石的應力積累達到較高水平。有研究表明,斷層活動的強弱是區域主壓應力強弱的表現,亦是孕育地震強弱的一個標志。2.2斜邊觀測的觀測在力學中,應變狀態是以主應變、剪切應變、主應變方向等參數來描述的。在應力場作用下,剪應變的變化更能體現巖石破裂過程中的作用,雖然剪應變不可直接觀測到,但可通過伸縮儀3個方向的應變組合觀測,計算出平面應變參數:剪應變、主應變、最大剪應變以及主應變方向等,從而分析區域應變場的變化,進一步研究地表平面應變狀態與地震和構造的關系。目前的應變臺網的定點伸縮儀觀測站具有3個方向觀測,除了2個正交方向的觀測外,還有一個斜邊觀測。根據測得的3個方向的應變值,按照主應變、剪應變計算公式,即可計算出主應變、剪應變、主應力方向。圖7所示的觀測方向為目前大部分觀測站的測線方向,以張家口和易縣臺的觀測方向為例,張家口觀測站斜邊觀測為NW向,易縣觀測站斜邊觀測為NE向。為探索文安地震前的應變狀態變化,本文特對距離震中較近的易縣地震臺和較遠的張家口地震臺的洞體應變進行了應變各參數(最大主應變、最小主應變、最大剪應變和主應變方向)進行了計算。圖8為易縣臺平面應變參數的計算結果,可以看出,最大主應變、最小主應變和最大剪應變曲線在文安地震前后都出現了轉折。在2010年初也發生同步轉折;主應變方向在2006年7月出現大幅度的轉折變化,其后恢復,至2012年12月,未再出現明顯的異常變化。圖9為張家口臺的各種應變參數曲線。可以看出,最大主應變、最小主應變和最大剪應變在2001—2008年變化較為平穩,2009年出現同步轉折,加速上升或下降,2011年再一次發生轉折。主應變方向在2005年出現向下的加速變化,2009年出現了轉折,其后變化較穩定。以上所述的異常變化在時間上可對應2005年8月31日蔚縣4.0級地震,該地震距張家口臺約100km,而2006年7月4日文安地震前該臺沒有異常顯示。3多臺及臺內不同應變參數同步偏轉的可能結果本文分析了文安地震前200km范圍內觀測站及張家口地震臺的洞體應變、鉆孔應變和體應變的變化形態,從震前的分鐘值變化,中短期的日均值變化到應變的各參數變化,在華北北部、特別是首都圈地區應變資料出現了大范圍、同步轉折變化,部分臺站出現畸變變化。其現象表明,文安地震前在區域應力、震源應力的共同作用下,大范圍地殼不同層位巖石的應力積累達到較高水平。離震中最遠的張家口NS向、NW向觀測曲線和最大主應變、最小主應變、最大剪應變曲線在2009年同步出現了加速上升或下降,且變化幅度很大,2011年初又同步發生轉折。而離震中較近的易縣NS向、EW向觀測曲線2011年發生同步轉折,最大主應變和最大剪應變在2010年也出現了同步加速上升變化。以上2個臺及臺內不同應變參數之間發生的同步轉折,可以認為這種變化不是偶然的,表明了首都圈西部區域應變場在2011年前后發生了變化。但首都圈及附近地區近年來除文安地震外,并沒有中強以上地震發生,經分析認為這可能反映了引發周邊特大地震的板塊構造應力場變化,在時間上恰好與2011年3月11日日本本州東海岸附近海域發生的M9.0地震對應。日本9級地震發生在日本島弧的日本海溝地震帶,該地震帶屬于西太平洋地震帶,日本島弧位于太平洋板塊、菲律賓板塊、歐亞板塊與北美板塊相沖的地方。日本島弧與中國東部地塊構造及應力場具有統一性,這無論是從板塊構造角度還是地質力學的角度都能解釋。在日本海溝地區,由于太平洋板塊的斜向俯沖運動,產生持續形變、積累應力,并孕育大地震。隨著太平洋板塊俯沖運動的繼續,在板塊邊界發生變形的同時,陸地一側的應變也不斷地增加。可以認為,太平洋板塊的俯沖運動為歐亞大陸地震的發生提供了力源。華北地區位于歐亞板塊東面,它受到太平洋板塊向西俯沖的力較為直接。因此,多個臺及各個臺內不同應變參數之間發生的同步轉折現象,可能反映的是引發周邊特大地震的板塊構造應力場的變化。通過上述分析,可以得到以下幾點認識:(1)在文安地震前,各種分析證明,距震中200km范圍內的臺站部分伸縮儀的原始觀測曲線和計算的應變參數曲線在震前均出現張性或者壓性轉折變化,或者階躍變化。(2)地表潮汐平面應變與地震