章 鑫 杜學(xué)彬 張元生 張麗峰

1) 中國蘭州730000中國地震局蘭州地震研究所 2) 中國蘭州730000蘭州地球物理國家野外科學(xué)觀測研究站

?

京津冀地區(qū)多種地球物理觀測對天津爆炸事件的響應(yīng)

章 鑫1)杜學(xué)彬1,2),張元生1,2)張麗峰1)

1) 中國蘭州730000中國地震局蘭州地震研究所 2) 中國蘭州730000蘭州地球物理國家野外科學(xué)觀測研究站

本文利用京津冀地區(qū)地球物理觀測數(shù)據(jù)和衛(wèi)星熱紅外遙感數(shù)據(jù)， 研究了2015年8月12日天津港地面爆炸時的地球物理響應(yīng)， 獲得了爆炸點周邊的地球物理場擾動信息． 結(jié)果表明： ① 事件點附近地電阻率對地面爆炸的響應(yīng)較為明顯， 事件發(fā)生后地電阻率觀測曲線幅值急劇變小， 且持續(xù)很長時間； ② 事件發(fā)生的前后兩天， 地電場觀測曲線出現(xiàn)了小幅升降異常， 且與空間電子流量微小幅度下降的時間相同， 而相應(yīng)時段的地磁場則無相關(guān)擾動； ③ 衛(wèi)星熱紅外亮溫在事件發(fā)生后明顯增強， 持續(xù)數(shù)天后恢復(fù)至事件前水平； ④ 事件點周圍的地下水位和形變測量未見明顯響應(yīng)． 上述結(jié)果表明不同的地球物理觀測對該事件的響應(yīng)有差異， 這對于分析和研判與地震等相關(guān)的異常具有參考意義， 同時本文討論的觀測結(jié)果也有利于深入認(rèn)知天津爆炸事件的過程．

地面爆炸 地電場 衛(wèi)星紅外遙感 地下流體 地電阻率 形變

引言

2015年8月12日23時34分， 天津市濱海區(qū)發(fā)生了地面爆炸事件(以下簡稱"事件")(中國發(fā)展觀察， 2015)， 該事件激發(fā)的震動波及到整個京津冀地區(qū)． 中國地震局在京津冀地區(qū)布設(shè)了大量的測震、 地面電磁(地電場、 地電阻率、 地磁場)、 地下流體和地形變等地球物理觀測臺站， 獲得了多種地球物理觀測數(shù)據(jù)； 同時中國氣象局靜止氣象衛(wèi)星(FY-2系列)在事件發(fā)生前后也產(chǎn)出了熱紅外觀測數(shù)據(jù)． 在事件發(fā)生時， 測震臺站的記錄明顯響應(yīng)該事件， 但其它地球物理觀測是否也響應(yīng)了該事件， 是值得研究的問題， 這將有助于認(rèn)知事件本身及其所激發(fā)的地球物理場變化和不同地球物理手段對地震等災(zāi)害事件的監(jiān)測能力．

在以地震監(jiān)測預(yù)報為目的地面電磁觀測和衛(wèi)星熱紅外觀測中， 已記錄了豐富的與地震發(fā)生、 火山噴發(fā)有關(guān)的前兆異常(屈春燕等， 2006； 杜學(xué)彬， 2010； 譚大誠等， 2010， 2014； 張元生等， 2010， 2011)， 同時在地面和電離層還觀測到了同震和震前電磁參量的擾動(湯吉等， 2008； 張學(xué)民等， 2010； 閆相相等， 2013)． 在地下流體的水位觀測中， 也得到了大量的與地震有關(guān)的前兆異常、 同震變化以及水位對潮汐的響應(yīng)(車用太等， 2006， 2011； 晏銳等， 2012). 地形變觀測直接測量地殼介質(zhì)的運動和變形， 在長期觀測中也記錄了與地震有關(guān)的前兆異常和大量的同震變化以及潮汐響應(yīng)(薄萬舉等， 1994； 孫伶俐等， 2013)． 縱觀上述多學(xué)科手段的觀測事實， 無論前兆異常、 同震變化還是潮汐響應(yīng)等均為地下介質(zhì)運動和變形的直接測量或由此激發(fā)的地球物理現(xiàn)象． 從這個角度來講， 天津地面爆炸事件也可能引起多學(xué)科觀測手段的地球物理響應(yīng)， 因此， 本文利用距離此次爆炸事件點鄰近的京津冀地區(qū)的地面電磁觀測、 衛(wèi)星熱紅外、 地下流體、 地形變觀測數(shù)據(jù)以及空間天氣觀測數(shù)據(jù)， 分析和研究該區(qū)域地球物理觀測對該事件的響應(yīng)．

1 觀測方法和背景資料

本文使用京津冀地區(qū)多種地球物理觀測數(shù)據(jù)， 包括2個地電場臺站、 3個地電阻率臺站、 3個地磁場臺站、 3個地下流體臺站(點)和1個形變觀測臺站記錄的數(shù)據(jù)以及8天連續(xù)的衛(wèi)星熱紅外遙感數(shù)據(jù)和1周空間天氣監(jiān)測數(shù)據(jù). 研究區(qū)各觀測臺站分布如圖1所示， 箭頭所指為實際爆炸場景(中國發(fā)展觀察， 2015). 另外， 使用2014年相應(yīng)農(nóng)歷時段的地電阻率觀測值作時間參考， 使用與該事件相距1000 km外的杭州地磁臺2015年8月6—20日觀測值作遠(yuǎn)空間參考．

圖1 研究區(qū)觀測臺站分布圖Fig.1 Distribution of stations in studied region

1.1 地電和地磁觀測

在地電場和地電阻率觀測中引入自然電場測量和對稱四極測量探測方法， 監(jiān)測地下介質(zhì)電性結(jié)構(gòu)隨時間的變化．

地電場臺站在地表布設(shè)NS和EW兩個水平正交測量方向和一個斜交測量方向， 每個方向有長、 短兩種極距， 長極距一般為200—400 m， 短極距約為長極距的一半； 測量儀器頻帶為0—0.005 Hz， 儀器分辨率為10 μV/km， 數(shù)據(jù)產(chǎn)出不低于1次/(min·測道)， 每日數(shù)據(jù)量為1440個． 近年來增加了地電場測量頻率， 采樣率為1 Hz， 能夠記錄到更加豐富的地電場變化信息．

地電阻率觀測采用對稱四極的布極方式， 觀測臺站在地表布設(shè)NS和EW兩個水平正交測項， 供電極距一般約為1000 m． 測量儀器分辨率為0.01 Ω·m， 數(shù)據(jù)產(chǎn)出為1次/h， 每日產(chǎn)出24個地電阻率數(shù)據(jù)． 數(shù)據(jù)處理時對地電阻率曲線取幅值計算， 表達(dá)式為

(1)

式中，S(i)為采樣序列曲線，i為數(shù)據(jù)點，k為區(qū)間[m,n]內(nèi)取最值(最大值和最小值)的個數(shù)．

地磁場觀測為水平分量H、 磁偏角D和地磁場總量F觀測， 其記錄頻率為1 Hz， 每日產(chǎn)出8萬6400個數(shù)據(jù)． 地磁觀測采樣率較高， 記錄到的磁變化信息豐富， 并且與空間電離層擾動和電流系的關(guān)系密切．

1.2 衛(wèi)星熱紅外遙感探測

衛(wèi)星熱紅外遙感儀器能直接觀測地表物體的熱輻射度， 經(jīng)定標(biāo)處理和幾何校正后的數(shù)據(jù)為相應(yīng)通道(波段)的輻射亮度． 應(yīng)用普朗克輻射定理的黑體輻射公式能計算獲得輻射物的輻射溫度(張元生等， 2010)． 為了與物體的真實溫度區(qū)別， 稱之為亮溫， 代表把輻射物體當(dāng)作黑體時的輻射溫度． 熱輻射強度除了與物體溫度和比輻射率有關(guān)， 還與物體表面曲率有關(guān)． 短時間內(nèi)， 同一物體在同一方向的熱輻射變化主要取決于物體本身的溫度．

本文選取靜止氣象衛(wèi)星FY-2C的遠(yuǎn)熱紅外亮溫數(shù)據(jù)(國家衛(wèi)星氣象中心， 2015)． FY-2C衛(wèi)星于2004年10月發(fā)射， 定點于105°E上空， 距地面約3.5萬千米， C衛(wèi)星的有效觀測范圍為(45°E—165°E， 60°S—60°N)． 靜止衛(wèi)星自旋掃描輻射儀熱紅外1， 2通道所對應(yīng)的波段為10.3—11.3 μm和11.5—12.5 μm， 星下點分辨率為5 km， 一次觀測范圍約為地球表面積的1/3(張元生等， 2011)． 相鄰兩次觀測的時間間隔為30分鐘或1小時， 每天至少進(jìn)行24次觀測． 為了減小太陽的影響， 選擇當(dāng)日23:00至次日4:00的多次觀測數(shù)據(jù)， 并用補窗法去除或減弱云層的影響， 對相同像元計算其日平均值， 構(gòu)成最終的亮溫日值．

1.3 地下流體和形變觀測

地下流體觀測包括地下靜水位、 逸出氣氡濃度、 水溫、 電導(dǎo)率和水中氡濃度變化等測項， 一般在井下觀測， 受地表人為活動的干擾較少． 這些數(shù)據(jù)中含有豐富的地殼動力信息， 特別是井水位動態(tài)數(shù)據(jù)中含有地球固體潮、 大氣壓力變化、 地震波、 前驅(qū)波、 地表水體荷載作用和斷層蠕動等信息， 為地球科學(xué)研究提供了有價值的信息資源(車用太等， 2006)． 本文采用塘沽臺、 靜海臺和王3井的流體觀測數(shù)據(jù)， 流體觀測采樣頻率為1次/min， 單測道每日產(chǎn)出1440個數(shù)據(jù). 處理數(shù)據(jù)時應(yīng)盡可能尊重原始數(shù)據(jù)， 對時序曲線僅作必要的去突跳處理．

地應(yīng)變和地傾斜等定點形變測量指標(biāo)能直接反映地殼介質(zhì)的微動態(tài)變化， 捕捉到地殼介質(zhì)破裂前的力學(xué)變化信息． 觀測項目包括地傾斜、 洞體應(yīng)變、 鉆孔應(yīng)變和臺站重力等; 觀測儀器有DSQ 型水管傾斜儀、 VS 型垂直擺傾斜儀和DZW 型微伽重力儀等(孫伶俐等， 2013)． 本文所用數(shù)據(jù)為2015年8月12—13日天津薊縣小辛莊山洞形變觀測點的數(shù)據(jù)， 該臺站距離事件點約120 km， 每日產(chǎn)出1440個數(shù)據(jù)， 其測量頻率為1次/min．

2 地球物理觀測對事件的響應(yīng)

2.1 地電阻率和地電場響應(yīng)

塘沽地電阻率臺站距事件點約18 km． 圖2a給出了2015年8月10—15日塘沽臺地電阻率觀測值的原始曲線． 通過對比圖2b--d可以看出， 事件發(fā)生前數(shù)天地電阻率具有較大的波動幅度， NS道達(dá)3 Ω·m， EW道較小， 為1 Ω·m. 地電阻率曲線存在明顯的白天波動大、 夜間波動小的特征， 其大、 小幅值變化的時間轉(zhuǎn)換點與輕軌列車開始運行和最后收車的時間點相吻合， 因此該現(xiàn)象是輕軌列車運營通電(或漏電)的響應(yīng)(謝凡等， 2011； 張宇等， 2014). 總體來看， 在該事件發(fā)生后的數(shù)日內(nèi)， 地電阻率幅值大幅度減小， 而其平均值與事件發(fā)生前數(shù)日的平均值相近．

表1為2015年8月6—20日期間每日的地電阻率幅值(1天計算1次幅值)統(tǒng)計結(jié)果． 可以看出， 12—13日塘沽臺所有測向的幅值均發(fā)生大幅度的變化， 13日NS測道地電阻率幅值僅為12日幅值的1/2， 幅值減小近1.1 Ω·m． EW測道幅值減小更大， 最大的為自然電位差和EW向地電阻率均方差， 12日為13日的近3倍． 但是， 距事件點60 km的青光臺并未出現(xiàn)幅值減小的現(xiàn)象， 其地電阻率觀測值和地電阻率均方差均保持平穩(wěn)過渡或略有增大．

圖3給出了2015年8月6—20日及2014年7月18日—8月1日(時間參考段)地電阻率和地電阻率均方差幅值變化情況. 可以看出： 從事件發(fā)生的12日夜間， 塘沽臺地電阻率幅值大幅減小， 12日前和13日后幅值均處于平穩(wěn)水平； 地電阻率均方差由于本身數(shù)值較小， 幅值減小的現(xiàn)象不夠明顯， 但其減小倍數(shù)與地電阻率減小倍數(shù)(表1加黑行)近似． 青光臺地電阻率及其均方差幅值曲線并未出現(xiàn)幅值減小的現(xiàn)象．

地電場觀測臺站距事件點較遠(yuǎn)， 其中靜海臺距事件點約65 km， 寶坻臺離事件點約80 km． 圖4給出了事件發(fā)生前幾十分鐘內(nèi)靜海臺和寶坻臺測得的地電場響應(yīng)情況． 可以看出， EW測道地電場曲線出現(xiàn)了增大波動(簡稱“突增”)， 但與事件本身無關(guān)， 作為地表突發(fā)事件， 不可能出現(xiàn)地電場前兆． 張世中等(2013)和沈紅會等(2014)的研究指出， 由于受到城鐵輕軌漏電的干擾， 地電場觀測曲線在6:00—23:00時段內(nèi)疊加了高頻成分， 如圖4a中靜海臺EW測道的地電場觀測曲線在6:00—23:00時段內(nèi)疊加了大量的高頻信號， 使得天然場源的高頻部分難以被識別． 圖4中除了靜海臺NS測道外， 兩個臺站的其它測道在事件發(fā)生前幾十分鐘內(nèi)均疊加了該信號， 地電場曲線從下降的趨勢上突增， 形成一個較大的不明信號源(圖4中8月13日前)．

圖2 2015年8月10—15日地電阻率臺站記錄的原始曲線， 圖中虛線表示天津爆炸事件發(fā)生時刻

寶坻臺位于事件點正北偏西方向， 靜海臺位于事件點S80°W附近(圖1)． 值得指出的是， 事件發(fā)生后位于事件點正北附近的寶坻臺NS和EW兩個方向均接收到了這一突增信號(圖4b)， 而位于正西的靜海臺只有EW向測道接收到了該信號(圖4a)， 這一現(xiàn)象說明該突增信號具有明顯的方向性， 即EW向強， NS向弱． 另外， 8月13日對應(yīng)時間點未出現(xiàn)相似的突增信號， 且位于事件點約N45°E方向的昌黎臺NS和EW兩個方向在相應(yīng)時間點也未出現(xiàn)突增信號(圖中未給出)． 由于地電場和地磁場的變化成分很大程度上與空間電流體系的擾動相關(guān)(孫正江， 王華俊， 1984)， 故有必要參考相應(yīng)時段的空間天氣情況．

表1 2015年8月6—20日塘沽臺和青光臺地電阻率幅值統(tǒng)計結(jié)果

注：ρ為地電阻率，αρ為地電阻率均方差，Vsp為自然電位差．

圖3 2015年8月6—20日塘沽臺(a)、 青光臺(b)和2014年7月18日—8月1日塘沽臺(c)地電阻率ρ及其均方差αρ幅值變化曲線

圖4 2015年8月12—13日靜海臺(a)和寶坻臺(b)記錄的地電場觀測曲線

2.2 空間電離層和地磁響應(yīng)

圖5給出了國家空間天氣監(jiān)測預(yù)警中心(2015)發(fā)布的風(fēng)云衛(wèi)星監(jiān)測周報， 包括電子流量和地磁場Kp指數(shù)監(jiān)測圖． 可以看出：Kp指數(shù)在事件發(fā)生時間附近均未出現(xiàn)明顯的異常， 8月15日夜間至16日凌晨發(fā)生Kp=7的中強磁暴， 與該事件無關(guān)； 8月12日23:00—24:00時段電子流量E2曲線出現(xiàn)微小幅度下降， 時間點為事件發(fā)生時刻附近； 8月13日凌晨2:00—3:00時段兩曲線均出現(xiàn)小幅向下波動，E1曲線波動幅度大于E2曲線， 時間節(jié)點約為事件發(fā)生后2.5小時， 可能與地面爆炸沖擊電離層進(jìn)而影響電子流量有關(guān)． 與之對應(yīng)的鄰近事件點的地磁觀測并未出現(xiàn)磁擾(圖6)． 8月13日12時左右， 寶坻臺、 徐莊子臺和遠(yuǎn)距離的杭州臺H分量均出現(xiàn)小型擾動， 但該時段Kp指數(shù)僅為4， 由此可見該擾動的原因并不明確； 圖4中靜海臺和寶坻臺地電場觀測曲線在13日12時左右也出現(xiàn)了小幅上升， 相同的變化在12日12時左右并未出現(xiàn)， 故推斷此次磁擾與事件本身無關(guān)．

圖5 2015年8月10—16日空間天氣情況(引自國家空間天氣監(jiān)測預(yù)警中心， 2015)

圖6 2015年8月12—13日地磁變化時序曲線

2.3 衛(wèi)星熱紅外響應(yīng)

圖7為2015年8月11—18日研究區(qū)衛(wèi)星熱紅外遙感亮溫異常分布情況. 為了消除太陽的影響， 對當(dāng)日23:00至次日4:00的數(shù)據(jù)取平均作為該日的亮溫值， 故圖中所示的亮溫均為已去掉背景值的亮溫異常． 可以看出： 8月11日無亮溫異常； 12日出現(xiàn)了小幅亮溫異常； 13日和14日的亮溫異常范圍增加、 幅度增大； 15日以后異常范圍和幅度逐步減??； 至18日幾乎恢復(fù)至事件發(fā)生前11日的亮溫水平． 綜上可知， 異常出現(xiàn)的演變過程為“無(11日)→出現(xiàn)(12日)→增加/穩(wěn)定(13—15日)→衰減(16—17日)→基本消失(18日)”． 由于熱輻射的波段不同， 在大氣中傳播的速度也不相等， 所以熱輻射的最大能量并不在第一時間出現(xiàn)在亮溫異常上， 事件發(fā)生時間與異常出現(xiàn)時間之間有一小段時間滯后． 該亮溫異常過程可以理解為： ① 事件發(fā)生前幾乎無熱紅外亮溫異常； ② 事件發(fā)生后， 地面出現(xiàn)強烈的熱輻射， 但衛(wèi)星只探測到輻射的一部分， 亮溫圖出現(xiàn)異常； ③ 地面事件的輻射開始變?nèi)酰?衛(wèi)星探測到大部分強烈輻射， 出現(xiàn)了較大幅度和范圍的熱紅外異常； ④ 在事件發(fā)生3—4天后， 地面事件的熱輻射極弱， 衛(wèi)星探測到的異常逐步衰減； ⑤ 事件發(fā)生7天后， 衛(wèi)星幾乎探測不到熱紅外亮溫， 亮溫異常消失．

圖7 研究區(qū)2015年8月11—18日熱紅外亮溫異常分布

事件引起熱輻射的原因可能有兩種： 一是化學(xué)品爆炸直接產(chǎn)生的熱量， 包括化學(xué)品和廠房等其它物質(zhì)劇烈燃燒生成的大量熱； 二是未爆炸的化學(xué)品持續(xù)遇水反應(yīng)生成的熱， 化學(xué)反應(yīng)的能量沒有爆炸的能量大， 但其持續(xù)時間極長， 對采樣率較低的衛(wèi)星熱紅外亮溫貢獻(xiàn)較大．

2.4 地下形變響應(yīng)

天津薊縣小辛莊臺的形變觀測包括水平擺傾斜、 垂直擺傾斜、 水管傾斜和洞體形變等測項， 該臺距離事件發(fā)生點約120 km． 圖8為其觀測結(jié)果的原始時序曲線. 可以看出: 8月12日水平擺傾斜的兩個方向均未出現(xiàn)明顯的異常； 14日0時， 兩個方向的曲線均較12日0時呈小幅下降， 且NS向下降更為明顯(圖8a)． 垂直擺傾斜測項也出現(xiàn)類似的變化形態(tài)， 其中NS向的測量值下降較大， 而EW向測量值則無明顯的下降趨勢(圖8b)， 這可能是由于事件引起的地面震動主要體現(xiàn)在垂直擺NS方向上， 而小辛莊測量點正好位于事件點正北方向略偏西(圖1)， 所以垂直擺EW方向測量值未出現(xiàn)明顯響應(yīng)． 距離事件發(fā)生點更近的地區(qū)沒有形變測量點， 而小辛莊測點距離事件發(fā)生點較遠(yuǎn)， 所以尚難明確是否觀測到事件的形變響應(yīng)． 由于這兩個測道的時序觀測曲線在事件發(fā)生時間點前后數(shù)十分鐘均無明顯的跳動出現(xiàn)， 曲線保持平滑， 故不能確定NS向測值下降是否與該事件有關(guān)．

圖8 2015年8月12—13日小辛莊臺形變觀測時序曲線

由圖8c可以看出， 13日水管傾斜曲線幅度大于12日， 但前后數(shù)日的趨勢均是如此， 這可能是由儀器零點漂移等因素引起的， 與本次事件無關(guān)． 由圖8d可以看出， 洞體應(yīng)變觀測結(jié)果中未出現(xiàn)明顯的突變異常， 13日EW方向觀測值與12日相比呈小幅降低， 但由于引起洞體應(yīng)變變化的因素很多， 諸如大氣壓、 地震等， 所以不能確定該形變測量的變化是否與此次事件有關(guān)． 形變測量的所有測項中均不能明確看出與該事件相關(guān)的應(yīng)變響應(yīng)， 這可能是由于臺站距離事件發(fā)生點太遠(yuǎn)， 處于事件能量波及范圍以外．

2.5 地下流體響應(yīng)

在爆炸事件發(fā)生前后， 京津冀地區(qū)的地下流體觀測并未出現(xiàn)明顯的異常， 如圖9所示． 塘沽地下流體觀測點位于事件附近， 自7月下旬靜水位處于持續(xù)下降狀態(tài)， 其響應(yīng)情況類似于小辛莊臺形變觀測的水平擺和垂直擺傾斜測項的NS測道， 14日0時測值較12日0時有所下降. 靜海臺氣壓也出現(xiàn)同樣的變化規(guī)律， 但其下降趨勢并非事件后才出現(xiàn)， 而是持續(xù)時間較長． 其余測點的靜水位則比較平穩(wěn)， 如王3井． 塘沽臺靜水位下降的原因可能是儀器零點漂移或局部氣壓的影響， 與此次事件無關(guān)． 值得指出的是， 小辛莊臺垂直擺傾斜NS測道觀測曲線與塘沽臺靜水位曲線具有相同的趨勢， 該現(xiàn)象似乎表明靜水位測項的變化趨勢也具有方向性， 但與圖4中地電場觀測到的方向變化是否具有統(tǒng)一的機制， 還有待進(jìn)一步深入研究．

圖9 2015年8月6—20日各臺站流體觀測曲線. 右圖為左圖虛線部分的放大圖

3 討論與結(jié)論

本文通過分析天津港爆炸事件點鄰近地區(qū)的地面電磁、 衛(wèi)星熱紅外、 地下流體和地形變臺站在事件前后的觀測記錄， 得到了不同的地球物理觀測對本次事件的響應(yīng)： ① 從事件發(fā)生時刻起， 距離事件點18 km的塘沽臺地電阻率曲線每日幅值減小， 其自然電位差和地電阻率均方差的幅值和均值均大幅下降， 但距離事件點60 km外的青光臺和寶坻臺并未觀測到類似現(xiàn)象； ② 在事件發(fā)生后幾分鐘內(nèi)， 距事件點65—80 km的靜海臺、 寶坻臺和較遠(yuǎn)的昌黎臺等地電場觀測曲線同步出現(xiàn)了瞬間脈沖變化， 在事件發(fā)生前幾十分鐘內(nèi)， 靜海臺和寶坻臺EW測道地電場觀測曲線出現(xiàn)了先小幅度增大(突增)后又恢復(fù)的波動變化， 寶坻臺NS測道也出現(xiàn)類似信號； ③ 在上述靜海臺等地電場觀測曲線出現(xiàn)擾動時， 地磁場觀測曲線并未出現(xiàn)擾動響應(yīng)， 但在事件發(fā)生后的2.5小時， 電離層電子流量E2曲線出現(xiàn)了小幅度下降， 4小時候后E1和E2曲線均出現(xiàn)中等幅度的下降擾動， 此時地磁場和地電場均未觀測到擾動； ④ 衛(wèi)星熱紅外亮溫于12日11時開始增強， 13—14日亮溫持續(xù)增強， 直到16日后亮溫開始衰減， 18日后亮溫恢復(fù)至事件發(fā)生前水平； ⑤ 地下流體觀測未出現(xiàn)明確的與此次事件對應(yīng)的擾動變化， 形變觀測點距事件點超過120 km， 也未出現(xiàn)異常變化． 這些地球物理響應(yīng)的可能原因如下：

1) 電離層高能電子流量的波動與衛(wèi)星熱紅外亮溫異常的機理是吻合的． 地面大能量的爆炸沖擊以及因注水形成的物質(zhì)間的持續(xù)化學(xué)反應(yīng)放熱， 其能量擴散到空間引起衛(wèi)星熱紅外亮溫增加和電離層電子流量的微小幅度下降(林大超等， 2001； 閆相相等， 2013)． 空間電子流量微小幅度下降的影響是瞬時和小范圍的， 可能未影響到地磁場或者引起的地磁場變化不在觀測頻段內(nèi)． 首先， 在事件發(fā)生后幾分鐘內(nèi)， 事件點周圍的靜海、 寶坻和昌黎臺記錄到同步地電場脈沖變化， 而在遠(yuǎn)距離的其它臺站則未觀測到同步脈沖變化， 顯然此脈沖變化不是由空間電磁擾動引起的； 其次， 靜海、 寶坻和昌黎臺均記錄到了同步脈沖， 說明不是由臺站觀測環(huán)境所引起的． 因此， 上述3個臺站出現(xiàn)的同步脈沖變化是緣于檢測到了由本次爆炸事件激發(fā)出的局部空間范圍的電磁效應(yīng)． 在事件發(fā)生前幾十分鐘內(nèi)， 靜海臺和寶坻臺的地電場觀測曲線出現(xiàn)小幅波動變化， 其原因可能是直接觀測到了爆炸事件前其它事件對地電場的影響， 或者是地電場對電離層電子流量微小幅度下降的響應(yīng)． 在13日電子流量大幅度下降時， 地電場上升變化幅度沒有地磁場下降幅度明顯， 因此可以排除后一種可能性． 對爆炸事件后亮溫增加、 電子流量下降的另一種解釋為， 與地磁暴類似， 空間電子流量下降擾動引起了地磁場水平分量H大幅下降， 對應(yīng)的地電場垂直分量受影響最大， 而對地電場水平分量影響則較?。?衛(wèi)星熱紅外亮溫增強能夠證實地面事件的沖擊確實存在熱輻射(馬曉靜等， 2009； 張元生等， 2010)， 這將引起空間電子流量波動， 所以認(rèn)為該事件引起了強烈的熱輻射， 此后是持續(xù)的化學(xué)反應(yīng)放熱， 二者疊加使熱輻射量增加， 導(dǎo)致熱紅外亮溫異常于14—15日達(dá)到最強， 至16日后逐漸消散．

2) 塘沽臺地電阻率每日幅值和均方差減小的原因可能有兩方面： 其一是近距離的地鐵津濱線停止運行， 由該線路所引起的干擾消失； 其二是爆炸事件的沖擊可能造成了該臺站地下介質(zhì)微裂隙的變化． 在塘沽臺地電阻率觀測中， 每日5時附近地電阻率及其均方差和自然電場開始出現(xiàn)大幅度的變化， 此時間段與每日地鐵津濱線運行時間是一致的． 根據(jù)張世中等(2013)和沈紅會等(2014)的研究可知， 每日在地鐵津濱線運營時段內(nèi)疊加在塘沽臺地電阻率日均值上的跳動變化屬于地鐵干擾， 是由地鐵運行時線路供電所引起的對地漏電所致． 爆炸事件發(fā)生后地鐵津濱線隨即停運， 相應(yīng)的塘沽臺地電阻率每日幅值和均方差均明顯減小(圖2)， 證明了從爆炸事件開始的這些變化是由地鐵停運造成的． 但是， 在事件發(fā)生72小時后， 塘沽臺地電阻率觀測中的地鐵干擾仍未完全消失， 僅干擾幅度明顯減小(圖2—3)， 其減小但仍未消失的原因可能是該臺站除了受最近距離(6 km)的地鐵津濱線影響外， 還受較遠(yuǎn)距離的其它地鐵運營的影響， 但因其距離較遠(yuǎn)， 所以造成的干擾幅度較?。?由此看來， 地電阻率觀測中避開或抑制地鐵干擾是相當(dāng)重要的． 其次， 錢書清等(2003)根據(jù)巖石實驗推斷， 若地下相對穩(wěn)定的介質(zhì)裂隙受到地震等大的不可逆應(yīng)力作用后， 裂隙將受力貫通， 宏觀上地電阻率減小， 且長時間保持低地電阻率水平． 據(jù)地下各向異性介質(zhì)中電阻率變化與微裂隙和裂隙水的本構(gòu)關(guān)系可知， 介質(zhì)中微裂隙貫通和水進(jìn)入微裂隙是引起視電阻率減小的直接原因(杜學(xué)彬等， 2007； 杜學(xué)彬， 2010)． 天津濱海地區(qū)東北部沿海岸向南沉積物粒度逐漸變細(xì)， 從砂和粉砂質(zhì)砂到砂質(zhì)粉砂， 再到黏土質(zhì)粉砂(岳軍等， 2011)， 具有富水條件， 水中含鹽量也高(即低阻水)， 造成了介質(zhì)具有良好的導(dǎo)電性． 塘沽臺地下介質(zhì)就屬于這種情況， 該臺地電阻率觀測供電極距AB=1500 m， 探測深度/范圍較大(杜學(xué)彬等， 2007)， 探測層介質(zhì)是水飽和或趨于飽和的， 兩個測道的實測地電阻率值在10—14 Ω·m之間， 并存在明顯的方向性． 由阿契(Archie)定律可知， 當(dāng)介質(zhì)飽和或趨于飽和的情況下， 再增加水含量對巖石電阻率的影響并不大． 在這種情況下， 根據(jù)杜學(xué)彬等(2007)和杜學(xué)彬(2010)給出的本構(gòu)關(guān)系， 介質(zhì)電阻率變化主要取決于介質(zhì)微裂隙發(fā)育和優(yōu)勢取向． 所以本文認(rèn)為爆炸事件所引起的震動波及到近距離的塘沽臺， 可能導(dǎo)致了該臺站地下介質(zhì)中微裂隙的發(fā)育和擴張， 在宏觀上微裂隙優(yōu)勢取向， 同時低阻水進(jìn)入微裂隙， 改變了地電阻率對每日潮汐的響應(yīng)， 進(jìn)而導(dǎo)致每日地電阻率幅值減小等現(xiàn)象． 在地鐵津濱線恢復(fù)試運行后， 該臺站的干擾影響又再次增加， 但干擾幅度卻明顯小于爆炸事件之前， 因此也不完全排除介質(zhì)微裂隙變化和低阻水進(jìn)入微裂隙的影響．

3) 地下流體臺站(包括塘沽臺)未觀測到對應(yīng)此次事件的明確響應(yīng)， 其原因可能是爆炸事件本身未造成地下流體的“變異”現(xiàn)象． 地下水位等流體觀測與固體潮汐和大氣壓的相關(guān)性高(晏銳等， 2012； 劉紹華等， 2014)， 地面爆炸的瞬時能量遠(yuǎn)大于正常大氣壓值(林大超等， 2001； 胡宏偉等， 2014)． 南方周末(2015)報道該事件造成了附近2 km范圍內(nèi)部分建筑物的門窗玻璃破碎、 墻壁受損等， 說明爆炸事件激發(fā)了強烈的沖擊能量． 沖擊能量由地面向四周傳播的同時也向地下傳播， 但是包括近距離的塘沽臺(不超過20 km)在內(nèi)的地下流體臺站卻未出現(xiàn)地下水位等響應(yīng)， 由此推測其可能原因是該地區(qū)地下介質(zhì)對沖擊能量的衰減較強(李洪濤等， 2010； 張先武等， 2014)， 此次爆炸的沖擊能量未能激發(fā)出地下水位等響應(yīng)， 這與電磁和熱紅外對本次事件的響應(yīng)機制不同． 再者， 除塘沽臺外其它地下流體臺站均位于60 km以外， 地下水位等也無異常． 地形變測量也有類似原因， 小辛莊臺距事件點約120 km， 通過地下介質(zhì)的衰減作用， 事件的沖擊能量未能傳播至很遠(yuǎn)距離． 對本次事件， 包括近距離的塘沽臺在內(nèi)的靜水位觀測未出現(xiàn)對該事件的響應(yīng)， 其原因目前尚不明確． 實際上， 在地震監(jiān)測中僅從地磁場觀測曲線上識別出同震變化的情況也不多， 這與地電場是不同的．

另外， 天津爆炸事件發(fā)生前幾十分鐘靜海臺和寶坻臺均出現(xiàn)了地電場小幅“突增”現(xiàn)象， 是否與事件發(fā)生前的其它事件有關(guān)， 還有待進(jìn)一步深入研究．

中國地震前兆臺網(wǎng)提供了相關(guān)前兆數(shù)據(jù)， 國家空間天氣監(jiān)測預(yù)警中心提供了空間天氣周報， 審稿專家提出了寶貴意見， 作者在此一并表示感謝．

薄萬舉, 謝覺民, 楊炳順. 1994. 地殼形變測量用于地震預(yù)報的研究[J]. 地殼形變與地震, 14(2): 13--19.

Bo W J, Xie J M, Yang B S. 1994. Study on earthquake prediction by use of the data of crustal deformation[J].CrustalDeformationandEarthquake, 14(2): 13--19 (in Chinese).

車用太, 趙文忠, 魚金子, 劉成龍. 2006. 京津冀地區(qū)井水位的數(shù)字化觀測及其地震前兆監(jiān)測效能評估[J]. 地震, 26(4): 103--112.

Che Y T ， Zhao W Z， Yu J Z， Liu C L. 2006. Digitalized well-water-level observation and monitoring efficiency evaluation of earthquake precursor in the Beijing-Tianjin-Hebei region[J].Earthquake, 26(4): 103--112 (in Chinese).

車用太, 魚金子, 劉成龍, 徐桂明, 鄭益銘. 2011. 判別地下水異常的干擾性與前兆性的原則及其應(yīng)用實例[J]. 地震學(xué)報, 33(6): 800--808.

Che Y T, Yu J Z, Liu C L, Xu G M, Zheng Y M. 2011. Principles on distinguishing interference from seismic precursor of underground water variation and its application[J].ActaSeismologicaSinica, 33(6): 800--808 (in Chinese).

杜學(xué)彬， 李寧， 葉青， 馬占虎， 閆睿． 2007． 強地震附近視電阻率各向異性變化的原因[J]. 地球物理學(xué)報， 50(6)： 1802--1810.

Du X B， Li N， Ye Q， Ma Z H， Yan R. 2007. A possible reason for the anisotropic changes in apparent resistivity near the focal region of strong earthquake[J].ChineseJournalofGeophysics， 50(6)： 1802--1810 (in Chinese).

杜學(xué)彬. 2010. 在地震預(yù)報中的兩類視電阻率變化[J]. 中國科學(xué): 地球科學(xué), 40(10): 1321--1330.

Du X B. 2011. Two types of changes in apparent resistivity in earthquake prediction[J].ScienceChina:EarthSciences, 54(1): 145--156.

國家空間天氣監(jiān)測預(yù)警中心. 2015. 空間天氣周報[EB/OL]. [2015-08-17]. http:∥img.nsmc.org.cn/spacedata/SWWR/SEVP_NSMC_SWWR_NCSW_ESWE_AGLB_LNO_P9_20150818010016800.pdf.

National Center for Space Weather. 2015. Space weather weekly report[EB/OL]. [2015-08-17]. http:∥img.nsmc.org.cn/spacedata/SWWR/SEVP_NSMC_SWWR_NCSW_ESWE_AGLB_LNO_P9_20150818010016800.pdf (in Chinese).

國家衛(wèi)星氣象中心. 2015. 衛(wèi)星紅外亮溫數(shù)據(jù)[EB/OL]. [2015-08-18]. http:∥satellite.cma.gov.cn/PortalSite/Sup/User/LoginUser.aspx?url=http:∥satellite.cma.gov.cn/PortalSite/Data/ShoppingCart.aspx.

National Satellite Meteorological Center. 2015. Satellite infrared brightness temperature data[EB/OL]. [2015-08-18]. http:∥satellite.cma.gov.cn/PortalSite/Sup/User/LoginUser.aspx?url=http:∥satellite.cma.gov.cn/PortalSite/Data/ShoppingCart.aspx (in Chinese).

胡宏偉, 宋浦, 郭煒, 馮海云, 張立建. 2014. 地面爆炸沖擊波的相互作用[J]. 高壓物理學(xué)報, 28(3): 353--357.

Hu H W, Song P, Guo W, Feng H Y, Zhang L J. 2014. Interaction of shock waves in ground burst[J].ChineseJournalofHighPressurePhysics, 28(3): 353--357 (in Chinese).

李洪濤, 盧文波, 舒大強, 楊興國, 易長平. 2010. 爆破地震波的能量衰減規(guī)律研究[J]. 巖石力學(xué)與工程學(xué)報, 29(增刊1): 3364--3369.

Li H T, Lu W B, Shu D Q, Yang X G, Yi C P. 2010. Study of energy attenuation law of blast-induced seismic wave[J].ChineseJournalofRockMechanicsandEngineering, 29(Suppl 1): 3364--3369 (in Chinese).

林大超, 白春華, 張奇. 2001. 爆炸地震地面豎向振動的幅值特性[J]. 北京理工大學(xué)學(xué)報, 21(2): 173--176.

Lin D C, Bai C H, Zhang Q. 2001. Amplitude characteristics of vertical ground surface vibration caused by explosion seism[J].JournalofBeijingInstituteofTechnology, 21(2): 173--176 (in Chinese).

劉紹華, 梁杏, 朱常坤, 馬斌, 徐敏, 高業(yè)新. 2014. 衡水地區(qū)淺層高頻地下水位動態(tài)的氣壓響應(yīng)及其特征[J]. 安全與環(huán)境工程, 21(1): 1--8.

Liu S H, Liang X, Zhu C K, Ma B, Xu M, Gao Y X. 2014. Atmospheric pressure response and characteristics of the high frequency level dynamic of shallow groundwater in Hengshui region[J].SafetyandEnvironmentalEnginee-ring, 21(1): 1--8 (in Chinese).

馬曉靜, 鄧志輝, 陳梅花, 楊竹轉(zhuǎn), 高祥林. 2009. 從衛(wèi)星紅外亮溫與大地?zé)崃鞯年P(guān)系看地震前的熱紅外異常[J]. 地球物理學(xué)報, 52(11): 2746--2751.

Ma X J, Deng Z H, Chen M H, Yang Z Z, Gao X L. 2009. A perspective to thermal infrared anomalies before earthquakes from the relationship between satellite infrared brightness temperature and terrestrial heat flow[J].ChineseJournalofGeophysics, 52(11): 2746--2751 (in Chinese).

南方周末. 2015. 國家級環(huán)境應(yīng)急戰(zhàn)三十天: 修復(fù)天津港[EB/OL]. [2015-09-10]. http:∥www.infzm.com/content/111706.

Southern Weekend. 2015. National environmental emergency war continued thirty days: Repair Tianjin Port[EB/OL]. [2015-09-10]. http:∥www.infzm.com/content/111706 (in Chinese).

錢書清, 郝錦綺, 周建國, 高金田. 2003. 巖石受壓破裂的ULF和LF電磁前兆信號[J]. 中國地震, 19(2): 109--116.

Qian S Q, Hao J Q, Zhou J G, Gao J T. 2003. Precursory electric and magnetic signals at ULF and LF bands during the fracture of rocks under pressure[J].EarthquakeResearchinChina, 19(2): 109--116 (in Chinese).

屈春燕, 單新建, 馬瑾. 2006. 衛(wèi)星熱紅外遙感在火山活動性監(jiān)測中的應(yīng)用[J]. 地震地質(zhì), 28(1): 99--110.

Qu C Y, Shan X J, Ma J. 2006. Application of satellite thermal infrared remote sensing in detection of volcano activity[J].SeismologyandGeology, 28(1): 99--110 (in Chinese).

沈紅會, 李鴻宇, 袁慎杰, 繆阿麗. 2014. 江寧臺地電場地地鐵干擾試驗結(jié)果分析[J]. 地震, 34(3): 125--131.

Shen H H, Li H Y, Yuan S J, Miu A L. 2014. Testing results of subway interferences to the Jiangning electromagnetic observation site[J].Earthquake, 34(3): 125--131 (in Chinese).

孫伶俐, 李明, 蔣玲霞, 羅俊秋, 戴苗, 馬武剛. 2013. 湖北省潮汐形變觀測異常及干擾識別[J]. 大地測量與地球動力學(xué), 33(增刊I): 36--40.

Sun L L, Li M, Jiang L X, Luo J Q, Dai M, Ma W G. 2013. Anomaly recognition of tidal deformation and disturbance factors in Hubei Province[J].JournalofGeodesyandGeodynamics, 33(Suppl I): 36--40 (in Chinese).

孫正江, 王華俊. 1984. 地電概論[M]. 北京: 地震出版社: 3--95.

Sun Z J, Wang H J. 1984.IntroductiontoGeoelectricity[M]. Beijing: Seismological Press: 3--95 (in Chinese).

譚大誠, 趙家騮, 席繼樓, 杜學(xué)彬, 徐建明. 2010. 潮汐地電場特征及機理研究[J]. 地球物理學(xué)報, 53(3): 544--555.

Tan D C, Zhao J L, Xi J L, Du X B, Xu J M. 2010. A study on feature and mechanism of the tidal geoelectrical field[J].ChineseJournalofGeophysics, 53(3): 544--555 (in Chinese).

譚大誠, 趙家騮, 劉小鳳, 范瑩瑩, 劉君, 陳軍營. 2014. 自然電場的區(qū)域性變化特征[J]. 地球物理學(xué)報, 57(5): 1588--1598.

Tan D C, Zhao J L, Li X F, Fan Y Y, Liu J, Chen J Y. 2014. Features of regional variations of the spontaneous field[J].ChineseJournalofGeophysics, 57(5): 1588--1598 (in Chinese).

湯吉, 詹艷, 王立鳳, 徐建郞, 趙國澤, 陳小斌, 董澤義, 肖騎彬, 王繼軍, 蔡軍濤, 徐光晶. 2008. 5月12日汶川8.0級地震強余震觀測的電磁同震效應(yīng)[J]. 地震地質(zhì), 30(3): 739--744.

Tang J, Zhan Y, Wang L F, Xu J L, Zhao G Z, Chen X B, Dong Z Y, Xiao Q B, Wang J J, Cai J T, Xu G J. 2008. Coseismic signal associated with aftershock of theMS8.0 Wenchuan earthquake[J].SeismologyandGeology, 30(3): 739--744 (in Chinese).

謝凡, 滕云田, 徐平. 2011. 應(yīng)用獨立分量分析方法提取和剝離地磁觀測中的軌道交通干擾[J]. 地球物理學(xué)進(jìn)展, 26(5): 1824--1831.

Xie F, Teng Y T, Xu P. 2011. Removal of EM interference generated by urban railway transit from geomagnetic observation by ICA method[J].ProgressinGeophysics, 26(5): 1824--1831 (in Chinese).

晏銳, 張立, 簡春林. 2012. 云南曲靖井水位潮汐動態(tài)特征分析[J]. 地震學(xué)報, 34(3): 363--373.

Yan R, Zhang L, Jian C L. 2012. Dynamic analysis of well water level tides in Qujing well, Yunnan Province[J].ActaSeismologicaSinica, 34(3): 363--373 (in Chinese).

閆相相, 單新建, 曹晉濱, 湯吉, 王飛飛. 2013. 日本MW9.0級特大地震前電離層擾動初步分析[J]. 地球物理學(xué)進(jìn)展, 28(1): 155--164.

Yan X X, Shan X J, Cao J B, Tang J, Wang F F. 2013. Preliminary study of the seismoionospheric perturbation before Tohoku-OkiMW9.0 earthquake[J].ProgressinGeophysics, 28(1): 155--164 (in Chinese).

岳軍, Dong Y, 張寶華, 牟林, 王國明, 陳安蜀, 袁寶印, 劉景蘭, 魏俊浩. 2011. 渤海灣西岸幾種地球化學(xué)的環(huán)境指標(biāo)[J]. 地質(zhì)學(xué)報, 85(7): 1239--1250.

Yue J, Dong Y, Zhang B H, Mu L, Wang G M, Chen A S, Yuan B Y, Liu J L, Wei J H. 2011. Several geochemical indicators of the west coast of Bohai Sea[J].ActaGeologicaSinica, 85(7): 1239--1250 (in Chinese).

張世中, 石航, 王蘭煒, 胡哲, 劉大鵬, 魏連生, 鞠永. 2013. 地電臺站受城市軌道交通干擾的測試分析與抗干擾措施研究[J]. 地震學(xué)報, 35(1): 117--124.

Zhang S Z, Shi H, Wang L W, Hu Z, Liu D P, Wei L S, Ju Y. 2013. Test analysis on disturbances caused by urban rail transit at geoelectric stations and measures to reduce its influence[J].ActaSeismologicaSinica, 35(1): 117--124 (in Chinese).

張學(xué)民, 劉靜, 申旭輝, Parrot M, 錢家棟, 歐陽新艷, 趙庶凡, 黃建平. 2010. 2005年3月28日蘇門答臘8.6級地震前的電離層擾動[J]. 地球物理學(xué)報, 53(3): 567--575.

Zhang X M, Liu J, Shen X H, Parrot M, Qian J D, Ouyang X Y, Zhao S F, Huang J P. 2010. Ionospheric perturbations associated with theM8.6 Sumatra earthquake on 28 March 2005[J].ChineseJournalofGeophysics, 53(3): 567--575 (in Chinese).

張先武, 高云澤, 方廣有. 2014. 消除探地雷達(dá)數(shù)據(jù)的子波衰減和頻散的反濾波方法[J]. 地球物理學(xué)報, 57(3): 932--938.

Zhang X W, Gao Y Z, Fang G Y. 2014. An inverse filtering method for removing the wavelet attenuation and dispersion of Ground Penetrating Radar data[J].ChineseJournalofGeophysics, 57(3): 932--938 (in Chinese).

張宇, 王蘭煒, 張興國, 朱旭, 劉大鵬, 顏蕊. 2014. 相關(guān)檢測技術(shù)在低頻交流地電阻率觀測中的應(yīng)用[J]. 地球物理學(xué)進(jìn)展, 29(4): 1973--1979.

Zhang Y, Wang L W, Zhang X G, Zhu X, Liu D P, Yan R. 2014. Application of correlation detection technology in low-frequency AC geo-resistivity observation[J].ProgressinGeophysics, 29(4): 1973--1979 (in Chinese).

張元生, 郭曉, 鐘美嬌, 沈文榮, 李穩(wěn), 何斌. 2010. 汶川地震衛(wèi)星熱紅外亮溫變化[J]. 科學(xué)通報, 55(10)： 904--910.

Zhang Y S, Guo X, Zhong M J, Shen W R, Li W, He B. 2010. Wenchuan earthquake: Brightness temperature changes from satellite infrared information[J].ChineseScienceBulletin, 55(18): 1917--1924.

張元生, 郭曉, 魏從信, 沈文榮, 惠少興. 2011. 日本9級和緬甸7.2級地震熱輻射表現(xiàn)特征[J]. 地球物理學(xué)報, 54(10): 2575--2580.

Zhang Y S， Guo X， Wei C X， Shen W R， Hui S X. 2011. The characteristics of seismic thermal radiation of JapanMS9.0 and MyanmarMS7.2 earthquake[J].ChineseJournalofGeophysics, 54(10): 2575--2580 (in Chinese).

中國發(fā)展觀察. 2015. 天津塘沽發(fā)生爆炸事故[EB/OL]. [2015-08-20]. http:∥www.chinado.cn/?p=2931.

China Development Observation. 2015. Explosion occurred in Tanggu, Tianjin[EB/OL]. [2015-08-20]. http:∥www.chinado.cn/?p=2931 (in Chinese).

Multiple geophysical observation response to the Tianjin explosion in Beijing-Tianjin-Hebei region

Zhang Xin1)Du Xuebin1,2),Zhang Yuansheng1,2)Zhang Lifeng1)

1)LanzhouInstituteofSeismology,ChinaEarthquakeAdministration,Lanzhou730000,China2)LanzhouNationalObservatoryofGeophysics,Lanzhou730000,China

The surface explosion occurred in Tanggu district of Tianjin on August 12, 2015 caused wide attention. This paper presented the geophysical observation response to this explosion based on the data from a number of geophysical observation stations in Beijing-Tianjin-Hebei region as well as the satellite infrared remote sensing data, and obtained the geophysical disturbances information around the point. The results showed that the response of georesistivity to the explosion near the event site is apparent, that is, the georesistivity curve declined rapidly after this event and lasted for a long time. It was also observed that the geoelectric field anomalies appeared with a small amplitude two days before and after the event, and the space electronic flow declined slightly synchronously, however, no disturbances of geomagnetic field appeared at the same period. Meanwhile, the satellite infrared brightness temperature significantly increased after the event, and lasted for several days. In addition, there was no obvious response of underground water level and deformation to the event. The results indicate that different types of observations have different responses to this event, it will be valuable for analyzing and determing the anomalies related to earthquake in the future. What is more, it is also useful for a deeper cognition for the process of Tianjin explosion.

surface explosion; geoelectric field; satellite infrared remote sensing; underground fluid; georesistivity; deformation

章鑫， 杜學(xué)彬， 張元生， 張麗峰. 2016. 京津冀地區(qū)多種地球物理觀測對天津爆炸事件的響應(yīng). 地震學(xué)報, 38(2): 283--297. doi:10.11939/jass.2016.02.012.

Zhang X, Du X B, Zhang Y S, Zhang L F. 2016. Multiple geophysical observation response to the Tianjin explosion in Beijing-Tianjin-Hebei region.ActaSeismologicaSinica, 38(2): 283--297. doi:10.11939/jass.2016.02.012.

國家自然科學(xué)基金(41374080)資助

2015-10-14收到初稿， 2016-02-24決定采用修改稿.

e-mail: duxb@163.com

10.11939/jass.2016.02.012

P315.72

A