• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

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

    2016-12-02 05:24:11杜學(xué)彬張元生張麗峰
    地震學(xué)報 2016年2期
    關(guān)鍵詞:塘沽亮溫臺站

    章 鑫 杜學(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), 獲得了爆炸點(diǎn)周邊的地球物理場擾動信息. 結(jié)果表明: ① 事件點(diǎn)附近地電阻率對地面爆炸的響應(yīng)較為明顯, 事件發(fā)生后地電阻率觀測曲線幅值急劇變小, 且持續(xù)很長時間; ② 事件發(fā)生的前后兩天, 地電場觀測曲線出現(xiàn)了小幅升降異常, 且與空間電子流量微小幅度下降的時間相同, 而相應(yīng)時段的地磁場則無相關(guān)擾動; ③ 衛(wèi)星熱紅外亮溫在事件發(fā)生后明顯增強(qiáng), 持續(xù)數(shù)天后恢復(fù)至事件前水平; ④ 事件點(diǎn)周圍的地下水位和形變測量未見明顯響應(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ì)的運(yùn)動和變形, 在長期觀測中也記錄了與地震有關(guān)的前兆異常和大量的同震變化以及潮汐響應(yīng)(薄萬舉等, 1994; 孫伶俐等, 2013). 縱觀上述多學(xué)科手段的觀測事實, 無論前兆異常、 同震變化還是潮汐響應(yīng)等均為地下介質(zhì)運(yùn)動和變形的直接測量或由此激發(fā)的地球物理現(xiàn)象. 從這個角度來講, 天津地面爆炸事件也可能引起多學(xué)科觀測手段的地球物理響應(yīng), 因此, 本文利用距離此次爆炸事件點(diǎn)鄰近的京津冀地區(qū)的地面電磁觀測、 衛(wèi)星熱紅外、 地下流體、 地形變觀測數(shù)據(jù)以及空間天氣觀測數(shù)據(jù), 分析和研究該區(qū)域地球物理觀測對該事件的響應(yīng).

    1 觀測方法和背景資料

    本文使用京津冀地區(qū)多種地球物理觀測數(shù)據(jù), 包括2個地電場臺站、 3個地電阻率臺站、 3個地磁場臺站、 3個地下流體臺站(點(diǎn))和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ù)點(diǎn),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)作黑體時的輻射溫度. 熱輻射強(qiá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ā)射, 定點(diǎn)于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, 星下點(diǎn)分辨率為5 km, 一次觀測范圍約為地球表面積的1/3(張元生等, 2011). 相鄰兩次觀測的時間間隔為30分鐘或1小時, 每天至少進(jìn)行24次觀測. 為了減小太陽的影響, 選擇當(dāng)日23:00至次日4:00的多次觀測數(shù)據(jù), 并用補(bǔ)窗法去除或減弱云層的影響, 對相同像元計算其日平均值, 構(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)變和地傾斜等定點(diǎn)形變測量指標(biāo)能直接反映地殼介質(zhì)的微動態(tài)變化, 捕捉到地殼介質(zhì)破裂前的力學(xué)變化信息. 觀測項目包括地傾斜、 洞體應(yīng)變、 鉆孔應(yīng)變和臺站重力等; 觀測儀器有DSQ 型水管傾斜儀、 VS 型垂直擺傾斜儀和DZW 型微伽重力儀等(孫伶俐等, 2013). 本文所用數(shù)據(jù)為2015年8月12—13日天津薊縣小辛莊山洞形變觀測點(diǎn)的數(shù)據(jù), 該臺站距離事件點(diǎn)約120 km, 每日產(chǎn)出1440個數(shù)據(jù), 其測量頻率為1次/min.

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

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

    塘沽地電阻率臺站距事件點(diǎn)約18 km. 圖2a給出了2015年8月10—15日塘沽臺地電阻率觀測值的原始曲線. 通過對比圖2b--d可以看出, 事件發(fā)生前數(shù)天地電阻率具有較大的波動幅度, NS道達(dá)3 Ω·m, EW道較小, 為1 Ω·m. 地電阻率曲線存在明顯的白天波動大、 夜間波動小的特征, 其大、 小幅值變化的時間轉(zhuǎn)換點(diǎn)與輕軌列車開始運(yùn)行和最后收車的時間點(diǎn)相吻合, 因此該現(xiàn)象是輕軌列車運(yù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倍. 但是, 距事件點(diǎn)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)象.

    地電場觀測臺站距事件點(diǎn)較遠(yuǎn), 其中靜海臺距事件點(diǎn)約65 km, 寶坻臺離事件點(diǎn)約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ā)生時刻

    寶坻臺位于事件點(diǎn)正北偏西方向, 靜海臺位于事件點(diǎn)S80°W附近(圖1). 值得指出的是, 事件發(fā)生后位于事件點(diǎn)正北附近的寶坻臺NS和EW兩個方向均接收到了這一突增信號(圖4b), 而位于正西的靜海臺只有EW向測道接收到了該信號(圖4a), 這一現(xiàn)象說明該突增信號具有明顯的方向性, 即EW向強(qiáng), NS向弱. 另外, 8月13日對應(yīng)時間點(diǎn)未出現(xiàn)相似的突增信號, 且位于事件點(diǎn)約N45°E方向的昌黎臺NS和EW兩個方向在相應(yīng)時間點(diǎn)也未出現(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的中強(qiáng)磁暴, 與該事件無關(guān); 8月12日23:00—24:00時段電子流量E2曲線出現(xiàn)微小幅度下降, 時間點(diǎn)為事件發(fā)生時刻附近; 8月13日凌晨2:00—3:00時段兩曲線均出現(xiàn)小幅向下波動,E1曲線波動幅度大于E2曲線, 時間節(jié)點(diǎn)約為事件發(fā)生后2.5小時, 可能與地面爆炸沖擊電離層進(jìn)而影響電子流量有關(guān). 與之對應(yīng)的鄰近事件點(diǎn)的地磁觀測并未出現(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)強(qiáng)烈的熱輻射, 但衛(wèi)星只探測到輻射的一部分, 亮溫圖出現(xiàn)異常; ③ 地面事件的輻射開始變?nèi)酰?衛(wèi)星探測到大部分強(qiáng)烈輻射, 出現(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ā)生點(diǎn)約120 km. 圖8為其觀測結(jié)果的原始時序曲線. 可以看出: 8月12日水平擺傾斜的兩個方向均未出現(xiàn)明顯的異常; 14日0時, 兩個方向的曲線均較12日0時呈小幅下降, 且NS向下降更為明顯(圖8a). 垂直擺傾斜測項也出現(xiàn)類似的變化形態(tài), 其中NS向的測量值下降較大, 而EW向測量值則無明顯的下降趨勢(圖8b), 這可能是由于事件引起的地面震動主要體現(xiàn)在垂直擺NS方向上, 而小辛莊測量點(diǎn)正好位于事件點(diǎn)正北方向略偏西(圖1), 所以垂直擺EW方向測量值未出現(xiàn)明顯響應(yīng). 距離事件發(fā)生點(diǎn)更近的地區(qū)沒有形變測量點(diǎn), 而小辛莊測點(diǎn)距離事件發(fā)生點(diǎn)較遠(yuǎn), 所以尚難明確是否觀測到事件的形變響應(yīng). 由于這兩個測道的時序觀測曲線在事件發(fā)生時間點(diǎn)前后數(shù)十分鐘均無明顯的跳動出現(xiàn), 曲線保持平滑, 故不能確定NS向測值下降是否與該事件有關(guān).

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

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

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

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

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

    3 討論與結(jié)論

    本文通過分析天津港爆炸事件點(diǎn)鄰近地區(qū)的地面電磁、 衛(wèi)星熱紅外、 地下流體和地形變臺站在事件前后的觀測記錄, 得到了不同的地球物理觀測對本次事件的響應(yīng): ① 從事件發(fā)生時刻起, 距離事件點(diǎn)18 km的塘沽臺地電阻率曲線每日幅值減小, 其自然電位差和地電阻率均方差的幅值和均值均大幅下降, 但距離事件點(diǎn)60 km外的青光臺和寶坻臺并未觀測到類似現(xiàn)象; ② 在事件發(fā)生后幾分鐘內(nèi), 距事件點(diǎn)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時開始增強(qiáng), 13—14日亮溫持續(xù)增強(qiáng), 直到16日后亮溫開始衰減, 18日后亮溫恢復(fù)至事件發(fā)生前水平; ⑤ 地下流體觀測未出現(xiàn)明確的與此次事件對應(yīng)的擾動變化, 形變觀測點(diǎn)距事件點(diǎn)超過120 km, 也未出現(xiàn)異常變化. 這些地球物理響應(yīng)的可能原因如下:

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

    2) 塘沽臺地電阻率每日幅值和均方差減小的原因可能有兩方面: 其一是近距離的地鐵津濱線停止運(yùn)行, 由該線路所引起的干擾消失; 其二是爆炸事件的沖擊可能造成了該臺站地下介質(zhì)微裂隙的變化. 在塘沽臺地電阻率觀測中, 每日5時附近地電阻率及其均方差和自然電場開始出現(xiàn)大幅度的變化, 此時間段與每日地鐵津濱線運(yùn)行時間是一致的. 根據(jù)張世中等(2013)和沈紅會等(2014)的研究可知, 每日在地鐵津濱線運(yùn)營時段內(nèi)疊加在塘沽臺地電阻率日均值上的跳動變化屬于地鐵干擾, 是由地鐵運(yùn)行時線路供電所引起的對地漏電所致. 爆炸事件發(fā)生后地鐵津濱線隨即停運(yùn), 相應(yīng)的塘沽臺地電阻率每日幅值和均方差均明顯減小(圖2), 證明了從爆炸事件開始的這些變化是由地鐵停運(yùn)造成的. 但是, 在事件發(fā)生72小時后, 塘沽臺地電阻率觀測中的地鐵干擾仍未完全消失, 僅干擾幅度明顯減小(圖2—3), 其減小但仍未消失的原因可能是該臺站除了受最近距離(6 km)的地鐵津濱線影響外, 還受較遠(yuǎn)距離的其它地鐵運(yù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ā)育和擴(kuò)張, 在宏觀上微裂隙優(yōu)勢取向, 同時低阻水進(jìn)入微裂隙, 改變了地電阻率對每日潮汐的響應(yīng), 進(jìn)而導(dǎo)致每日地電阻率幅值減小等現(xiàn)象. 在地鐵津濱線恢復(fù)試運(yùn)行后, 該臺站的干擾影響又再次增加, 但干擾幅度卻明顯小于爆炸事件之前, 因此也不完全排除介質(zhì)微裂隙變化和低阻水進(jìn)入微裂隙的影響.

    3) 地下流體臺站(包括塘沽臺)未觀測到對應(yīng)此次事件的明確響應(yīng), 其原因可能是爆炸事件本身未造成地下流體的“變異”現(xiàn)象. 地下水位等流體觀測與固體潮汐和大氣壓的相關(guān)性高(晏銳等, 2012; 劉紹華等, 2014), 地面爆炸的瞬時能量遠(yuǎn)大于正常大氣壓值(林大超等, 2001; 胡宏偉等, 2014). 南方周末(2015)報道該事件造成了附近2 km范圍內(nèi)部分建筑物的門窗玻璃破碎、 墻壁受損等, 說明爆炸事件激發(fā)了強(qiáng)烈的沖擊能量. 沖擊能量由地面向四周傳播的同時也向地下傳播, 但是包括近距離的塘沽臺(不超過20 km)在內(nèi)的地下流體臺站卻未出現(xiàn)地下水位等響應(yīng), 由此推測其可能原因是該地區(qū)地下介質(zhì)對沖擊能量的衰減較強(qiáng)(李洪濤等, 2010; 張先武等, 2014), 此次爆炸的沖擊能量未能激發(fā)出地下水位等響應(yīng), 這與電磁和熱紅外對本次事件的響應(yīng)機(jī)制不同. 再者, 除塘沽臺外其它地下流體臺站均位于60 km以外, 地下水位等也無異常. 地形變測量也有類似原因, 小辛莊臺距事件點(diǎn)約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. 強(qiáng)地震附近視電阻率各向異性變化的原因[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).

    李洪濤, 盧文波, 舒大強(qiáng), 楊興國, 易長平. 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ī)理研究[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級地震強(qiáng)余震觀測的電磁同震效應(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)用獨(dú)立分量分析方法提取和剝離地磁觀測中的軌道交通干擾[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

    猜你喜歡
    塘沽亮溫臺站
    中國科學(xué)院野外臺站檔案工作回顧
    霰譜分布特征對強(qiáng)對流云高頻微波亮溫影響的模擬研究
    氣象基層臺站建設(shè)
    西藏科技(2021年12期)2022-01-17 08:46:38
    基于南太平洋的AMSR2 L1R亮溫數(shù)據(jù)質(zhì)量評估
    海洋通報(2020年2期)2020-09-04 09:22:34
    2014年2月12日新疆于田MS7.3地震熱紅外亮溫異常分析
    中國地震(2015年1期)2015-11-08 11:11:32
    基層臺站綜合觀測業(yè)務(wù)管理之我見
    西藏科技(2015年6期)2015-09-26 12:12:13
    關(guān)于塘沽
    齊魯周刊(2015年33期)2015-09-06 00:39:39
    MDOS平臺臺站級使用方法及技巧
    2011年盈江5.8級地震熱輻射亮溫異常分析
    地震研究(2014年3期)2014-02-27 09:30:49
    書 畫
    天津教育(2014年1期)2014-01-18 06:29:52
    精品久久久久久,| 中出人妻视频一区二区| 黑人操中国人逼视频| 精品国产超薄肉色丝袜足j| 日韩欧美一区视频在线观看| av视频在线观看入口| 桃色一区二区三区在线观看| 中文资源天堂在线| 亚洲av成人不卡在线观看播放网| 久久久精品欧美日韩精品| 中文字幕人妻熟女乱码| 可以在线观看毛片的网站| 亚洲成人免费电影在线观看| 岛国视频午夜一区免费看| 欧美黑人巨大hd| 十分钟在线观看高清视频www| 欧美最黄视频在线播放免费| 日本免费一区二区三区高清不卡| 琪琪午夜伦伦电影理论片6080| 一本一本综合久久| 国产99久久九九免费精品| 国产精品98久久久久久宅男小说| 国产蜜桃级精品一区二区三区| 狠狠狠狠99中文字幕| 亚洲免费av在线视频| 看黄色毛片网站| 国产精品综合久久久久久久免费| 久久精品影院6| 久久性视频一级片| 日韩大尺度精品在线看网址| 亚洲中文字幕一区二区三区有码在线看 | 久久久久国产一级毛片高清牌| netflix在线观看网站| 亚洲专区中文字幕在线| 亚洲av成人av| 免费在线观看黄色视频的| 久久青草综合色| 欧美日韩黄片免| 国产成年人精品一区二区| 久久久久九九精品影院| 人妻丰满熟妇av一区二区三区| 一区二区三区激情视频| 黄片小视频在线播放| a在线观看视频网站| 丰满人妻熟妇乱又伦精品不卡| 国产精品国产高清国产av| 在线视频色国产色| 欧美在线一区亚洲| 两个人看的免费小视频| 一区二区三区高清视频在线| 午夜福利免费观看在线| 久久久国产成人精品二区| 国产成+人综合+亚洲专区| 国产男靠女视频免费网站| 日韩视频一区二区在线观看| 欧美一级毛片孕妇| 在线观看66精品国产| 免费av毛片视频| 视频在线观看一区二区三区| 高潮久久久久久久久久久不卡| 妹子高潮喷水视频| 亚洲精品国产精品久久久不卡| 久久九九热精品免费| 九色国产91popny在线| 日本黄色视频三级网站网址| 亚洲熟妇熟女久久| 国产日本99.免费观看| 久久久久久久久中文| 国产视频一区二区在线看| 国产成人系列免费观看| 满18在线观看网站| 中文字幕最新亚洲高清| 久久精品91蜜桃| 亚洲男人天堂网一区| 国产色视频综合| 国产精品一区二区免费欧美| 18禁观看日本| 国产亚洲精品久久久久5区| 国产精品亚洲av一区麻豆| 欧美乱色亚洲激情| 国产精品久久久久久精品电影 | 天天躁夜夜躁狠狠躁躁| 国产精品av久久久久免费| 国产亚洲欧美98| 亚洲av片天天在线观看| 女人高潮潮喷娇喘18禁视频| 亚洲黑人精品在线| 精品久久久久久久久久免费视频| 此物有八面人人有两片| 国内毛片毛片毛片毛片毛片| 欧美国产日韩亚洲一区| 悠悠久久av| 亚洲中文av在线| 中文字幕人妻熟女乱码| 99国产精品一区二区三区| 啦啦啦免费观看视频1| 色播亚洲综合网| 欧美国产日韩亚洲一区| 欧美一区二区精品小视频在线| 久久久久久久午夜电影| 一卡2卡三卡四卡精品乱码亚洲| 亚洲一区中文字幕在线| 国产精品久久久人人做人人爽| 99riav亚洲国产免费| 淫妇啪啪啪对白视频| 黄色 视频免费看| 亚洲中文av在线| 亚洲成a人片在线一区二区| 久久精品亚洲精品国产色婷小说| 日韩欧美一区视频在线观看| 欧美日韩一级在线毛片| 天天一区二区日本电影三级| 久久久久久久久久黄片| 欧美黄色淫秽网站| 国产精品亚洲一级av第二区| 哪里可以看免费的av片| 97超级碰碰碰精品色视频在线观看| 欧美绝顶高潮抽搐喷水| 亚洲激情在线av| 亚洲成av人片免费观看| 热re99久久国产66热| 国产亚洲精品第一综合不卡| 亚洲五月天丁香| 欧美久久黑人一区二区| 国产av在哪里看| 欧美zozozo另类| 国产av一区在线观看免费| 在线观看www视频免费| 国产午夜精品久久久久久| 久99久视频精品免费| 欧美日本视频| 午夜免费观看网址| 午夜激情福利司机影院| 国产乱人伦免费视频| 久久久久久久久久黄片| 国产精品亚洲av一区麻豆| 国产黄片美女视频| 国产亚洲欧美在线一区二区| 十分钟在线观看高清视频www| 激情在线观看视频在线高清| АⅤ资源中文在线天堂| 特大巨黑吊av在线直播 | 在线观看日韩欧美| 日韩精品中文字幕看吧| 哪里可以看免费的av片| 亚洲男人的天堂狠狠| 精品国产亚洲在线| 免费看美女性在线毛片视频| 欧美久久黑人一区二区| 搡老岳熟女国产| 午夜福利在线观看吧| 啦啦啦韩国在线观看视频| 精品人妻1区二区| 69av精品久久久久久| 九色国产91popny在线| 自线自在国产av| 黄频高清免费视频| 脱女人内裤的视频| 村上凉子中文字幕在线| 真人一进一出gif抽搐免费| 精品卡一卡二卡四卡免费| 啦啦啦免费观看视频1| 午夜免费激情av| 一区二区三区国产精品乱码| av欧美777| 两性夫妻黄色片| 久久久久久久久中文| 午夜福利欧美成人| 国产熟女午夜一区二区三区| 婷婷亚洲欧美| videosex国产| 久久午夜综合久久蜜桃| www.自偷自拍.com| 国产成人影院久久av| 真人做人爱边吃奶动态| 免费av毛片视频| 一进一出抽搐动态| 婷婷亚洲欧美| 久久久精品欧美日韩精品| 法律面前人人平等表现在哪些方面| 精品国产亚洲在线| 69av精品久久久久久| 淫妇啪啪啪对白视频| 国产三级在线视频| 制服诱惑二区| 黄色成人免费大全| 国内毛片毛片毛片毛片毛片| 国产精品电影一区二区三区| 一级a爱片免费观看的视频| 一本久久中文字幕| 少妇熟女aⅴ在线视频| 日韩欧美三级三区| 校园春色视频在线观看| 女人高潮潮喷娇喘18禁视频| 久9热在线精品视频| 亚洲va日本ⅴa欧美va伊人久久| 99热只有精品国产| 色综合欧美亚洲国产小说| 一区二区三区激情视频| 亚洲av熟女| 日韩大尺度精品在线看网址| 搡老岳熟女国产| 一本一本综合久久| 成人国产综合亚洲| 亚洲av电影在线进入| 色尼玛亚洲综合影院| 亚洲 国产 在线| а√天堂www在线а√下载| 91老司机精品| 国产精品久久久av美女十八| 制服诱惑二区| 亚洲五月色婷婷综合| 日韩视频一区二区在线观看| 免费观看精品视频网站| 啦啦啦免费观看视频1| 欧美人与性动交α欧美精品济南到| 欧美国产精品va在线观看不卡| 欧美黄色片欧美黄色片| 久久精品国产99精品国产亚洲性色| 亚洲 欧美 日韩 在线 免费| 国产区一区二久久| 老熟妇仑乱视频hdxx| 国产成人精品久久二区二区91| 一级毛片女人18水好多| 草草在线视频免费看| 日本一区二区免费在线视频| 欧美精品亚洲一区二区| 少妇裸体淫交视频免费看高清 | 好看av亚洲va欧美ⅴa在| 亚洲第一青青草原| 国产精品久久电影中文字幕| 此物有八面人人有两片| 麻豆成人午夜福利视频| 搞女人的毛片| 国产精品爽爽va在线观看网站 | 午夜福利欧美成人| 国产区一区二久久| 丝袜在线中文字幕| 亚洲国产欧美一区二区综合| 熟妇人妻久久中文字幕3abv| 淫秽高清视频在线观看| 看免费av毛片| 国产精品免费一区二区三区在线| 99热这里只有精品一区 | 国产熟女xx| 日韩欧美三级三区| 黄色丝袜av网址大全| 最好的美女福利视频网| 久99久视频精品免费| 变态另类成人亚洲欧美熟女| 满18在线观看网站| 91九色精品人成在线观看| 亚洲男人天堂网一区| 特大巨黑吊av在线直播 | 国产亚洲欧美精品永久| 亚洲真实伦在线观看| 久久九九热精品免费| 久久久久国产精品人妻aⅴ院| 久久久国产精品麻豆| 一本精品99久久精品77| 日韩欧美一区视频在线观看| 色综合欧美亚洲国产小说| 大型黄色视频在线免费观看| tocl精华| 精品国产亚洲在线| 国产精品久久久久久人妻精品电影| 亚洲av中文字字幕乱码综合 | 精品久久久久久久人妻蜜臀av| 日韩欧美一区二区三区在线观看| 精品午夜福利视频在线观看一区| 国产精品影院久久| 欧美日韩福利视频一区二区| 亚洲中文字幕一区二区三区有码在线看 | 97人妻精品一区二区三区麻豆 | 亚洲电影在线观看av| 亚洲av第一区精品v没综合| 麻豆一二三区av精品| 18禁国产床啪视频网站| 日本三级黄在线观看| 无遮挡黄片免费观看| 麻豆成人午夜福利视频| 999精品在线视频| 男女视频在线观看网站免费 | 免费看a级黄色片| 男人操女人黄网站| 精品一区二区三区视频在线观看免费| 18禁裸乳无遮挡免费网站照片 | 嫩草影院精品99| 国产成人av激情在线播放| 日韩欧美一区视频在线观看| 首页视频小说图片口味搜索| 日本五十路高清| 久9热在线精品视频| 亚洲人成网站在线播放欧美日韩| 母亲3免费完整高清在线观看| 91在线观看av| 亚洲无线在线观看| 成人精品一区二区免费| 色综合站精品国产| 日日爽夜夜爽网站| 99在线人妻在线中文字幕| 一本精品99久久精品77| 91麻豆精品激情在线观看国产| 国产v大片淫在线免费观看| 日韩欧美国产一区二区入口| 18禁黄网站禁片午夜丰满| 国产一级毛片七仙女欲春2 | 日本五十路高清| 欧美日韩福利视频一区二区| 国产高清视频在线播放一区| 欧美色视频一区免费| 在线观看免费视频日本深夜| 久久精品91无色码中文字幕| 超碰成人久久| 少妇裸体淫交视频免费看高清 | 国产精品99久久99久久久不卡| 免费女性裸体啪啪无遮挡网站| 国产99久久九九免费精品| 精品乱码久久久久久99久播| 最近最新中文字幕大全电影3 | 欧美黄色淫秽网站| 男人舔奶头视频| 中文字幕人妻熟女乱码| 国产爱豆传媒在线观看 | 50天的宝宝边吃奶边哭怎么回事| 国产精品二区激情视频| 精品久久久久久久久久免费视频| 国产精品一区二区三区四区久久 | 老汉色∧v一级毛片| 18禁裸乳无遮挡免费网站照片 | 久久精品影院6| 黑人欧美特级aaaaaa片| 99热只有精品国产| 欧美一级a爱片免费观看看 | 免费观看人在逋| 一区二区三区国产精品乱码| 亚洲人成77777在线视频| 成人欧美大片| 免费高清视频大片| 怎么达到女性高潮| 国产黄a三级三级三级人| 精品久久久久久久久久免费视频| 波多野结衣高清作品| 欧美久久黑人一区二区| 国产精品一区二区三区四区久久 | 一区二区日韩欧美中文字幕| svipshipincom国产片| 一区福利在线观看| 午夜成年电影在线免费观看| 国产蜜桃级精品一区二区三区| 国产成人啪精品午夜网站| 男人舔女人下体高潮全视频| 日本在线视频免费播放| 宅男免费午夜| 国内精品久久久久精免费| 99热只有精品国产| 免费观看精品视频网站| 欧美日韩中文字幕国产精品一区二区三区| 成人精品一区二区免费| 亚洲性夜色夜夜综合| 成人精品一区二区免费| 极品教师在线免费播放| 熟妇人妻久久中文字幕3abv| 大型av网站在线播放| 国产精品乱码一区二三区的特点| 亚洲黑人精品在线| 精品少妇一区二区三区视频日本电影| a在线观看视频网站| 一边摸一边做爽爽视频免费| 香蕉国产在线看| 精品一区二区三区视频在线观看免费| 国产色视频综合| 黄网站色视频无遮挡免费观看| 嫁个100分男人电影在线观看| 一边摸一边做爽爽视频免费| 亚洲色图 男人天堂 中文字幕| 精品国产超薄肉色丝袜足j| 国产色视频综合| 亚洲精品国产区一区二| 一进一出抽搐动态| 亚洲精品在线观看二区| 最新在线观看一区二区三区| 精品国产一区二区三区四区第35| 久久伊人香网站| 久久99热这里只有精品18| 国产精品永久免费网站| 久热这里只有精品99| 午夜福利成人在线免费观看| 99久久久亚洲精品蜜臀av| 俺也久久电影网| 91老司机精品| 欧美一区二区精品小视频在线| 搞女人的毛片| 99re在线观看精品视频| 日韩有码中文字幕| 女人爽到高潮嗷嗷叫在线视频| 亚洲av美国av| 99久久综合精品五月天人人| 国产一级毛片七仙女欲春2 | 操出白浆在线播放| 久久久久国产精品人妻aⅴ院| 十八禁网站免费在线| 欧美丝袜亚洲另类 | 黄色视频,在线免费观看| 日本在线视频免费播放| 亚洲专区中文字幕在线| 亚洲一码二码三码区别大吗| 国产精品野战在线观看| 久热这里只有精品99| 午夜老司机福利片| 曰老女人黄片| 香蕉av资源在线| e午夜精品久久久久久久| 身体一侧抽搐| 在线观看日韩欧美| 91成人精品电影| 可以免费在线观看a视频的电影网站| 国产精品久久久av美女十八| 高清在线国产一区| 亚洲男人天堂网一区| 日日夜夜操网爽| 亚洲国产欧洲综合997久久, | 香蕉丝袜av| 国产精品乱码一区二三区的特点| 1024视频免费在线观看| 身体一侧抽搐| e午夜精品久久久久久久| 美女免费视频网站| 日韩欧美国产在线观看| 啦啦啦 在线观看视频| 中文字幕高清在线视频| 欧美激情高清一区二区三区| 久久久久国产精品人妻aⅴ院| 中亚洲国语对白在线视频| 成年版毛片免费区| 一区二区三区精品91| 国产一级毛片七仙女欲春2 | 午夜a级毛片| 一本精品99久久精品77| АⅤ资源中文在线天堂| 国产亚洲精品一区二区www| 免费人成视频x8x8入口观看| 日本熟妇午夜| www.www免费av| 免费在线观看影片大全网站| 亚洲免费av在线视频| 精华霜和精华液先用哪个| 丰满人妻熟妇乱又伦精品不卡| 岛国视频午夜一区免费看| 国产精品野战在线观看| 51午夜福利影视在线观看| 黄色成人免费大全| av福利片在线| 国产精品久久久av美女十八| 一本精品99久久精品77| 国产精品日韩av在线免费观看| 国产熟女xx| 色播亚洲综合网| 动漫黄色视频在线观看| aaaaa片日本免费| 搡老岳熟女国产| 国产人伦9x9x在线观看| 在线天堂中文资源库| 韩国av一区二区三区四区| 99国产精品一区二区蜜桃av| 精品一区二区三区视频在线观看免费| 日本撒尿小便嘘嘘汇集6| 国产又色又爽无遮挡免费看| 淫妇啪啪啪对白视频| 露出奶头的视频| 不卡av一区二区三区| 18美女黄网站色大片免费观看| 12—13女人毛片做爰片一| 亚洲免费av在线视频| 午夜福利一区二区在线看| 国内精品久久久久精免费| 黄色视频,在线免费观看| 精品福利观看| 欧美日韩一级在线毛片| 精华霜和精华液先用哪个| 国产aⅴ精品一区二区三区波| 国产黄a三级三级三级人| 中文字幕精品免费在线观看视频| 亚洲精品久久国产高清桃花| 亚洲 国产 在线| 精品日产1卡2卡| 美女高潮喷水抽搐中文字幕| 亚洲国产日韩欧美精品在线观看 | 身体一侧抽搐| 性欧美人与动物交配| 亚洲激情在线av| 成人精品一区二区免费| 免费搜索国产男女视频| 久久久久久人人人人人| 欧美成人免费av一区二区三区| 亚洲,欧美精品.| 成人精品一区二区免费| 国产精品久久久久久人妻精品电影| 身体一侧抽搐| 中出人妻视频一区二区| 亚洲久久久国产精品| 亚洲一区二区三区不卡视频| 十八禁人妻一区二区| www.自偷自拍.com| 日本三级黄在线观看| 亚洲精品国产精品久久久不卡| 天天一区二区日本电影三级| 一本一本综合久久| 亚洲国产精品999在线| xxx96com| 制服丝袜大香蕉在线| 在线观看66精品国产| 18禁黄网站禁片午夜丰满| 国产精品香港三级国产av潘金莲| 免费在线观看日本一区| 成人午夜高清在线视频 | 亚洲一区高清亚洲精品| 在线免费观看的www视频| 91九色精品人成在线观看| 亚洲精品久久国产高清桃花| 午夜免费成人在线视频| 久久午夜综合久久蜜桃| 亚洲九九香蕉| 精品国产超薄肉色丝袜足j| 十分钟在线观看高清视频www| 91麻豆精品激情在线观看国产| 精品卡一卡二卡四卡免费| 最近最新免费中文字幕在线| 97碰自拍视频| 国产高清videossex| 成人亚洲精品一区在线观看| 变态另类成人亚洲欧美熟女| 免费在线观看亚洲国产| 国产91精品成人一区二区三区| 久久精品国产综合久久久| 视频在线观看一区二区三区| 国产av在哪里看| 91大片在线观看| 天天添夜夜摸| 亚洲 欧美 日韩 在线 免费| 国产精品亚洲av一区麻豆| 国产成年人精品一区二区| 亚洲自偷自拍图片 自拍| tocl精华| www日本黄色视频网| 精品国内亚洲2022精品成人| 成人午夜高清在线视频 | 精品午夜福利视频在线观看一区| 老汉色av国产亚洲站长工具| 岛国在线观看网站| 久99久视频精品免费| 草草在线视频免费看| 亚洲av片天天在线观看| 俺也久久电影网| 久久久国产欧美日韩av| 欧美一级a爱片免费观看看 | 日本三级黄在线观看| 露出奶头的视频| 久久香蕉精品热| xxxwww97欧美| 亚洲男人天堂网一区| 亚洲最大成人中文| 在线观看免费午夜福利视频| 黄色女人牲交| 搡老熟女国产l中国老女人| 国产色视频综合| 巨乳人妻的诱惑在线观看| 窝窝影院91人妻| 最近最新免费中文字幕在线| 女人高潮潮喷娇喘18禁视频| 91老司机精品| cao死你这个sao货| 日韩欧美国产一区二区入口| 国产又黄又爽又无遮挡在线| 久久人妻福利社区极品人妻图片| 亚洲精品国产区一区二| 精品久久久久久成人av| 久久久国产精品麻豆| 欧美成人午夜精品| 首页视频小说图片口味搜索| 99在线人妻在线中文字幕| 国语自产精品视频在线第100页| 久久久精品国产亚洲av高清涩受| 亚洲第一青青草原| www日本在线高清视频| 国产麻豆成人av免费视频| 香蕉av资源在线| 欧美 亚洲 国产 日韩一| 久99久视频精品免费| 香蕉av资源在线| 日本熟妇午夜| 热99re8久久精品国产| 国产91精品成人一区二区三区| 两性午夜刺激爽爽歪歪视频在线观看 | 久久久精品国产亚洲av高清涩受| 精品高清国产在线一区| a在线观看视频网站| 18禁黄网站禁片免费观看直播| 亚洲 欧美一区二区三区| 1024视频免费在线观看| 日韩成人在线观看一区二区三区| 日韩欧美国产在线观看| 亚洲七黄色美女视频| 亚洲 国产 在线| xxx96com| 99久久无色码亚洲精品果冻| 日本五十路高清| 精品国产乱子伦一区二区三区| 老鸭窝网址在线观看| 色在线成人网| 国产激情偷乱视频一区二区| 久久久久九九精品影院| 国产精品香港三级国产av潘金莲| 在线观看舔阴道视频| 国产真实乱freesex|