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

    2013年蘆山MS7.0地震前甘孜臺地電阻率變化分析1

    2013-09-09 06:37:30滔1軍1任越霞2趙民淵3
    震災防御技術 2013年4期
    關鍵詞:甘孜蘆山電阻率

    解 滔1) 盧 軍1) 任越霞2) 趙民淵3)

    ?

    2013年蘆山7.0地震前甘孜臺地電阻率變化分析

    解 滔盧 軍任越霞趙民淵

    1)中國地震臺網中心,北京100045 2)四川省地震局,成都 610041 3)甘孜地震臺,甘孜 626700

    甘孜臺地電阻率N30oE測道觀測資料自2011年7月開始出現趨勢上升變化,N60oW測道地電阻率則從2012年出現趨勢上升變化。2011年測區(qū)原317國道實施擴建工程,N30oE測道測量電極分別向供電極方向移動10m。采用甘孜臺電測深曲線以水平層狀模型反演了測區(qū)的電性結構,理論計算表明,測量電極的移動將會引起N30oE測道4W·m的上升變化,扣除這部分變化后,甘孜臺兩測道觀測值于2012年同步上升。以水平層狀模型計算了甘孜臺兩測道各層介質的影響系數,兩測道淺層兩層介質影響系數均為負,能合理地解釋甘孜臺地電阻率在雨季降水量增加時觀測值上升、旱季降水量減少時觀測值下降這一“夏高冬低”的年變現象。建立三維有限元模型計算了317國道拓寬部分對觀測的影響,計算結果表明,拓寬部分僅能引起N60oW測道約0.15W·m的下降變化和N30oE測道約0.1W·m的上升變化,其對觀測的影響非常小。同時2013年1月甘孜臺兩測道年變低值顯著高于2008年以來各年的年變低值,在蘆山地震前呈同步的上升變化,但是與汶川地震前的下降變化相反,因此甘孜臺自2012年的趨勢上升變化是不是蘆山地震的前兆異常還難以確定。

    蘆山地震 甘孜臺 地電阻率 趨勢上升 年變化

    引言

    2013年4月20日在四川省雅安市蘆山縣發(fā)生7.0地震,造成196人遇難,20余人失蹤和愈萬人受傷。此次地震是繼2008年汶川8.0地震后發(fā)生在龍門山斷裂帶的又一次巨大破壞性地震。據中國地震臺網測定,震中位于東經102o57′,北緯30o19′,震源深度為13km(http://www.ceic.ac.cn/),主震震中距汶川地震震中約87km。余震分布區(qū)沿長軸(龍門山斷裂走向)約45km,沿短軸(北西方向)約20km,余震分布區(qū)域距汶川地震余震區(qū)約37km(杜方等,2013;劉杰等,2013),震源機制研究結果表明,此次蘆山地震為純逆沖型(呂堅等,2013;曾祥方等,2013)。

    我國將地電阻率方法應用于地震中短期預測已有近50年,在長期的觀測實踐中多次記錄到了大震、中強地震前突出的異常現象(錢復業(yè)等,1982;1990;杜學彬,2010;張學民等,2009;桂燮泰等,1989;錢家棟等,2013;Lu等,1999)。關于地震前地電阻率異常機理的研究也相繼展開和不斷深入。從巖石介質微觀導電機理上,DD模式和IPE模式給出含水巖層在受壓應力作用時,介質電阻率呈下降變化,受張應力作用時電阻率呈上升變化(錢家棟,1985)。我國地電阻率觀測采用對稱四極觀測裝置,供電極距一般為AB=1000m左右,且認為孕震應力主要作用于測區(qū)深部巖層。深部巖層對地電阻率觀測的影響系數為正且大于表層介質的影響系數,因而理論上地電阻率觀測能記錄到與深部巖層電阻率受應力作用時同步的變化(趙和云等,1987;錢家棟等,1988;1998;Lu等,2004)。對巖石樣品的應力加載實驗也表明,隨著應力的增加巖石電阻率呈下降變化,同時在相同應力加載作用下,巖石不同方向的電阻率變化亦不同,呈現出各向異性變化(陸陽泉等,1990;陳峰等,2003;2013)。野外原地實驗結果顯示,巖層受壓應力作用時地電阻率觀測值呈下降變化,壓應力卸載和受張應力作用時地電阻率觀測值上升(趙玉林等,1983;1990)。

    近年來地電阻率觀測受到的干擾日趨嚴重,造成有些臺站資料失去了地震分析預測效能。蘆山地震震中約300km范圍內有成都、冕寧、西昌小廟和甘孜4個地電阻率臺站,其中成都臺和冕寧臺受干擾嚴重、西昌小廟臺于2012年年底停測,資料難以分析。甘孜臺觀測環(huán)境相對較好,但是觀測極距僅為AB=500m,同時由于測區(qū)下有一高阻層,使得觀測資料較易受來自地表的干擾,且呈現“夏高冬低”年變形態(tài)。甘孜臺地電阻率N30oE測道觀測值自2011年7月開始出現趨勢上升變化,N60oW測道則從2012年出現趨勢上升變化,2011年同期穿過測區(qū)的317國道實施擴建工程。本文將結合甘孜臺電測深資料討論其地電阻率變化的基本特征,并建立測區(qū)空間三維有限元模型分析國道擴建對觀測造成的影響,進而為分析蘆山地震前甘孜臺地電阻率的變化提供依據。

    1 甘孜地電阻率臺站簡介

    甘孜臺位于四川省甘孜縣城郊以東約4km處,臺站所在區(qū)域為鮮水河斷裂帶西北端、甘孜-馬尼干戈斷裂東端的斯俄盆地,臺站所在區(qū)域斷層構造如圖1所示。測區(qū)地貌平坦,臺站附近有三疊系變質砂巖、板巖出露。盆地第四紀沖積覆蓋層厚度可達40m以上,主要為含碎石粘土、亞粘土等高階地沖積物(李有才,1984)。臺站所在區(qū)域的鮮水河斷裂是我國地震非常活躍的一條地震帶,北西自甘孜東谷,向南經爐霍、道孚至康定磨西附近,全長約350km,整體呈左旋走滑特征,并發(fā)生多次7級以上地震(王貴宣,1995;張國民等,2005;Wen等,2008)。

    甘孜臺地電阻率觀測始于1970年,后因原測區(qū)受到金屬水管干擾,于1987年搬遷至現今位置(王德富,1990)。臺站布設N30oE和N60oW兩相互正交測道,采用對稱四極觀測裝置,并沿布極中心對稱分布(圖2)。兩測道供電極距AB=500m,原測量極距均為MN=120m,2011年317國道擴建后N30oE測道測量極距變更為140m,目前臺站使用ZD8BI數字地電儀進行觀測。自觀測以來在臺站周圍發(fā)生了多次中強地震,甘孜臺地電阻率記錄了較為突出的前兆異常。其中在1975年九龍6.2地震、1979年玉樹6.6地震、1981年道孚6.9地震和1982年甘孜6.0地震前N60oW測道出現了年尺度的下降異常(李有才,1984;王德富,1985);在1989年巴塘6.7雙震前兩測道同步出現了半年尺度的下降異常(王德富,1990);在2008年汶川8.0地震前兩測道出現約兩年的趨勢下降變化(張學民等,2009;杜學彬,2010;錢家棟等,2013)。

    圖1 甘孜臺地區(qū)斷層構造

    圖2 甘孜臺地電阻率布極示意圖

    2 甘孜地電阻率變化基本特征

    甘孜臺2007年至2013年6月的地電阻率和降雨量觀測資料如圖3所示。表層介質電阻率隨季節(jié)變化顯著,夏季降水量較多且表層水為液態(tài),介質電阻率相對較低,冬季降水較少且甘孜臺表層土壤處于凍結狀態(tài),介質電阻率較高。甘孜臺兩測道地電阻率在夏季降水量增加,淺層介質電阻率值降低時觀測值上升,而在冬季降水量減少和地表土層處于封凍,淺層介質電阻率值上升時觀測值下降,年變形態(tài)表現為“夏高冬低”,與我國大多數臺站地電阻率“夏低冬高”年變形態(tài)相反,同時年變幅度也較大,約為觀測值的8%。甘孜臺沿N30oE和N60oW兩方向的電測深曲線如圖4所示,兩方向電測深曲線均為K型,形態(tài)在淺層較為相近,中間含有高阻層。在深部電測深曲線顯示介質N60oW方向電阻率顯著低于N30oE,表現出電阻率各向異性。表1是依據圖4中電測深曲線以水平層狀模型計算并參考地電阻率觀測資料得到的測區(qū)地層電性結構。

    表1 甘孜臺水平層狀電性結構

    圖3 (a)甘孜臺N30oE測道地電阻率觀測值;(b)甘孜臺N60oW測道地電阻率觀測值;(c)甘孜臺月降雨量

    圖4 (a)甘孜臺N30oE方向電測深曲線;(b)甘孜臺N60oW方向電測深曲線

    如果將地電阻率測區(qū)劃分為任意的塊區(qū)域,每一塊區(qū)域介質電阻率為,=1,2,…,。在測區(qū)電性結構確定、觀測裝置和極距以及布極位置確定時,地電阻率是各分區(qū)介質電阻率的函數(錢家棟等,1985;1998;趙和云等,1987;Park等,1991):

    在多數情況下,各分區(qū)介質電阻率在一定時間內的相對變化非常小,,因此

    將(1)式進行Taylor級數展開,二階及高階項遠遠小于一階項,可忽略不計。地電阻率相對變化可以簡單地表示為各分區(qū)介質電阻率相對變化的加權和:

    采用表1層狀電性結構,計算N30oE和N60oW測道各層介質影響系數如圖5所示,甘孜臺兩測道在供電極距AB/2=250m時地表兩層介質對地電阻率的影響系數為負。因此在淺層介質電阻率受降雨影響季節(jié)性變化時,甘孜臺觀測的地電阻率表現出與之相反的年變化。同時兩測道地表兩層介質的影響系數也顯著大于其他臺站,表層介質相對變化量對地電阻率變化的貢獻較大,因而甘孜臺地電阻率的年變化幅度較大。

    3 317國道擴建對觀測的影響

    3.1 移動測量電極產生的影響

    原N30oE測道測量極距MN=120m,2011年317國道擴建,原測量電極位置被覆蓋,兩新測量電極位置各自向供電電極方向移動10m,移動后測量極距MN=140m。以表1中N30oE測道電性結構,采用二極裝置濾波器算法計算N30oE測道在測量電極移動前后地電阻率變化(O’Niell等,1984),計算結果如表2所示。計算結果表明,測量電極移動可以引起該測道地電阻率4W·m的階躍上升變化,將其后觀測資料扣除這部分變化后,N30oE測道在2011年的上升變化并不顯著。

    表2 測量電極移動產生的影響(W·m)

    3.2 國道擴建部分產生的影響

    317國道的擴建部分位于兩測道測線和電極附近,是地表地電阻率觀測較為敏感的區(qū)域,通常該區(qū)域介質電性結構的改變對觀測的影響不容忽視。公路施工時表層路基被壓實,這部分介質電阻率值應升高,同時路面為高阻混凝土結構,因此我們將這兩部分簡化為一厚度為2m的高阻體,建立三維有限元模型計算其對甘孜臺兩測道地電阻率觀測的影響。

    3.2.1 穩(wěn)恒電流場有限元方法

    地電阻率定點臺站觀測采用對稱四極裝置,觀測時在供電電極A、B輸入直流電流,在測量電極M、N測量電勢差,此問題可視為穩(wěn)恒電流場計算,電流場遵守Maxwell方程組和電荷守恒定律,電位分布滿足Possion方程:

    式中,是由電流源產生的電位;是介質電導率;是Dirac delta函數。

    應用虛功原理可得到穩(wěn)恒電流場Possion方程的有限元弱解形式:

    式中,W為計算區(qū)域;為任意的虛位移函數,在滿足Dirichlet邊界條件的邊界上,虛位移函數=0。

    由于電位在電流源處存在奇異性,數值求解方程(5)在電流源附近得到的結果誤差較大。常用的去除奇異性的方法是將電流源產生的電位分解為由供電電流在均勻半空間介質()中產生的一次場電位V和非均勻介質()引起的二次場電位V,且滿足=V+V和。一次場電位V可以通過解析表達式得到,二次場電位V通過單元離散化后數值計算得到(Lowry等,1989;Zhao等,1996)。

    地電阻率觀測在地表自然地滿足Neumann邊界條件,在水平方向和垂直方向(深度)可視為無窮遠邊界,可以施加Dirichlet邊界條件(=0),也可以施加Neumann邊界條件(Coggon,1971)。但是建立的模型在水平和垂直方向上的尺度不可能是無限的,對于一固定尺寸的模型,在供電極距AB大于一定的值后,對無窮遠邊界施加Dirichlet邊界條件時計算得到的地電阻率值將小于實際值,而對無窮遠邊界施加Neumann邊界條件時計算得到的地電阻率值將大于實際值(Li等,2005;Dey等,1979)。對固定的供電極距AB,模型尺寸越大,邊界對計算結果的影響就越小,但是模型越大,計算量也就越大,因此需要合理地選擇模型水平方向的尺寸和最底層厚度。

    模型經單元離散化、施加電流源和邊界條件后可對單元節(jié)點上的自由度(電位)進行有限元數值求解,求解出電位分布后可以獲得測量電極間的電勢差,進而依據對稱四極裝置系數計算地電阻率。

    3.2.2 模型建立

    甘孜臺測區(qū)淺層介質主要為沖積沉積物,N30oE和N60oW兩方向電測深曲線在淺層時形態(tài)相近,反映淺層介質電阻率各向異性不顯著。深部巖層由于長期受構造應力作用,微裂隙發(fā)育可能不均勻,造成深部巖層電阻率在不同方向有差異,兩方向電測深曲線顯示深部巖層N30oE方向電阻率大于N60oW方向。模型采用水平層狀結構,假定淺層三層介質電阻率均勻,僅在模型最底層采用不同的電阻率值計算兩方向的地電阻率值,模型參數如表1所示,圖6為模型示意圖。

    甘孜臺兩個測道的供電極距相同,且都沿布極中心對稱分布,因而模型在N60oW和N30oE方向采用相同的尺寸。圖7a是在模型最底層厚度為1000m時N60oW測道地電阻率計算值隨模型水平尺寸的變化,在>3000m之后N60oW測道地電阻率計算值基本上已不隨的增加而變化,所以模型中取=4000m。圖7b則是模型水平尺寸=4000m時N60oW測道地電阻率計算值隨模型最底層厚度的變化,在>1000m之后N60oW測道地電阻率計算值基本上也不隨的增加而變化,所以在模型中取=1500m。在計算甘孜臺兩測道地電阻率受國道擴建影響的最終模型大小為4000×4000×1551m,在建立模型計算時臺站觀測裝置布極中心始終位于模型水平面中心。

    3.2.3 計算結果

    計算時將國道擴建部分和路基壓實部分簡化為一厚度為2m的高阻體,電阻率值設定為10000W·m。先將高阻體部分電阻率設定為表層介質電阻率值,計算未受干擾時的地電阻率,而后將高電阻率值賦予之,計算受干擾后的地電阻率,從而計算出國道擴建部分對兩測道地電阻率的影響,計算結果如表3所示。317國道擴建部分僅能引起N30oE測道約0.1W·m的上升變化和N60oW測道約0.15W·m的下降變化變化。由于甘孜臺測區(qū)表層介質電阻率值較高,在高阻介質中嵌入高阻體,加之擴建部分的規(guī)模不大(寬度為4m),因而其對地電阻率觀測值的影響很小。N30oE測道在2011年8月17日移動電極,隨后測區(qū)內路基工程開始實施,直至9月完成,期間兩測道整點觀測值如圖8所示,從圖中可以看出施工期間甘孜臺兩測道地電阻率變化連續(xù)平穩(wěn),施工對觀測并沒有造成顯著的影響。

    表3 國道擴建部分產生的影響(W·m)

    甘孜臺兩測道地電阻率觀測值在扣除電極移動和317國道擴建部分產生的影響后,于2012年6月開始同步大幅度快速上升。2012年電磁學科組現場落實認為,當年6月、7月和9月三個月份降雨量顯著大于2008—2011年同月份降雨量,甘孜臺兩測道地電阻率同步快速上升與降雨量增加有關(電磁學科赴甘孜臺異常落實工作小組,2012)。直至2013年1月,測區(qū)幾無降雨,表層處于凍結時兩測道地電阻率年變低值仍高于往年,因而甘孜臺地電阻率在排除干擾后存在趨勢上升變化。

    4 討論

    甘孜臺測區(qū)淺層由沖積沉積物構成的地層電阻率較高,依據電測深曲線采用水平層狀模型反演的電性結構顯示第三層為高阻層,其下則為低阻巖層。高阻巖層的存在使得供電電流更多地集中在淺層介質中,加之臺站供電極距較小,這使得淺層介質對觀測到的地電阻率相對變化貢獻較多。從兩測道各層介質的影響系數可以看出,第三層介質影響系數最大,且大于1,第二層次之。一般認為由孕震引起的應力變化更多地影響深部巖層電阻率變化,甘孜臺深部巖層影響系數盡與表層介質相當,這對捕捉孕震信息是不利的。由于表層介質影響系數較大,甘孜臺的年變幅度也相應的較大,同時容易受到來自地表的干擾。甘孜臺為電磁學科在鮮水河上唯一的地電阻率臺站,對鮮水河上震情判定具有重要意義。從圖5可以看出,若將供電極距擴展至AB=1000m,最底層介質影響系數最大,且表層介質影響系數將降低一個數量級。測區(qū)覆蓋層厚度較薄,鉆孔將電極直接放置于第四層,實施井下觀測也是較好的選擇。

    從相互垂直的兩條電測深曲線(圖4)可以看出,甘孜臺電性結構在第四層存在顯著的各向異性特征,文中采用的是依據各條電測深曲線按一維層狀均勻介質模型反演的電性結構,它是對實際結構的一種簡化和近似。單一電測深曲線是對地下非均勻(或各向異性)結構的綜合反映,而以一維層狀均勻介質模型反演的電性結構是對地表觀測結果的一種解譯。甘孜臺地電阻率觀測測量極距的移動和國道擴建均在地表,因而分別采用一維模型予以計算的結果是可信的。

    甘孜臺距2008年汶川地震震中約330km,距蘆山地震震中約300km,汶川地震前甘孜臺地電阻率表現為中期下降變化異常(張學民等,2009),在蘆山地震前則表現為趨勢上升變化。汶川地震和蘆山地震都發(fā)生在龍門山斷裂帶上,發(fā)震動力源為巴彥喀拉地塊相對四川盆地的逆沖作用。震源機制解表明,汶川地震在破裂區(qū)中南段以逆沖為主,在北段以走滑為主,主震以逆沖型為主(陳運泰等,2008),蘆山地震發(fā)生在汶川地震以西南約85km處,主震為逆沖型(呂堅等,2013;曾祥方等,2013)。同時甘孜臺與兩次地震震中的距離幾乎一樣,兩次地震與甘孜臺的方位角也幾乎一樣,雖然震級上有差別,但兩次地震相對甘孜臺而言構成了大地震事件的重現性。然而甘孜臺地電阻率資料在汶川地震前為下降變化,而在蘆山地震前卻為上升變化,這與我們通常理解的物理現象重現性存在矛盾。地電阻率觀測值是局部范圍內測區(qū)地下介質電阻率的綜合反映,在排除測區(qū)干擾和年變產生的影響后,地電阻率變化被認為是深部巖層受應力作用時巖層電阻率變化引起的。理論模型(DD模式和IPE模型)以及實驗資料表明,含水巖層在壓應力積累階段電阻率總體呈下降變化,在壓應力卸載或張應力積累作用下電阻率呈上升變化。因此結合測區(qū)所在區(qū)域應力場資料來分析地電阻率變化更為客觀,甘孜臺有一體應變觀測,但觀測質量很差,無法獲取汶川地震和蘆山地震前測區(qū)深部巖層應變相對變化狀態(tài)。

    5 結論

    甘孜臺地電阻率測區(qū)下有一高阻層,使得甘孜臺在現有布極方式下淺層介質的影響系數為負,在降雨增加淺層介質電阻率降低時觀測值上升,降雨減少和淺層電阻率升高時觀測值降低,表現出“夏高冬低”的年變形態(tài)。同時由于下覆高阻層和觀測極距較短,淺層介質影響系數較大,從而年變幅度較大,且觀測易受來自地表的干擾。317國道擴建造成N30oE測道測量極距變化,從而引起該測道約4W·m的階躍上升變化。將擴建部分簡化為在地表嵌入的高阻體,有限元模型計算表明,擴建部分僅能引起N30oE測道約0.1W·m的上升變化和引起N60oW測道約0.15W·m的下降變化??鄢蓴_變化后,甘孜臺兩測道地電阻率資料自2012年6月開始同步出現快速上升變化,且年變峰值顯著高于前幾年,其中包含了2012年降雨量顯著大于前幾年引起的上升變化部分。在2013年1月幾無降雨、表層處于凍結時年變谷值也高于前幾年,因此甘孜臺兩測道存在趨勢上升變化。汶川地震前甘孜臺兩測道資料出現中期的下降變化,而蘆山地震前卻表現為上升變化,同時甘孜臺距震中距離約為300km,因此甘孜臺自2012年的趨勢上變化是不是蘆山地震的前兆異?,F在還難以予以確定。

    致謝:四川省地震局張茂軍和甘孜臺工作人員一同參與了甘孜臺的科學考察工作,并提供了一些重要的資料,作者在此一并表示衷心的感謝。

    陳峰,安金珍,廖椿庭,2003.原始電阻率各向異性巖石的電阻率變化各向異.地球物理學報,46(2):271—280.

    陳峰,馬麥寧,安金珍,2013.承壓介質電阻率變化的方向性與主應力的關系.地震學報,35(1):84—93.

    陳運泰,許力生,張勇,2008. 2008年5月12日汶川特大地震震源特性分析報告.

    電磁學科赴甘孜臺異常落實工作小組,2012.2012年10月7—12日四川甘孜臺地電阻率異常落實報告.震情研究,94:107—111.

    杜方,龍鋒,阮祥等,2013.四川蘆山7.0級地震及其與汶川8.0級地震的關系.地球物理學報,56(5):1772—1783.

    杜學彬,2010.在地震預報中的兩類視電阻率變化.中國科學(D輯),40(10):1321—1330.

    桂燮泰,關華平,戴經安,1989.唐山、松潘地震前視電阻率短臨異常圖像重現性.西北地震學報,11(4):71—75.

    李有才,1984.強震前甘孜臺視電阻率的前兆異常.地震,2:62—64.

    劉杰,易桂喜,張致偉等,2013.2013年4月20日四川蘆山7.0級地震介紹.地球物理學報,56(4):1404—1407.

    呂堅,王曉山,蘇金蓉等,2013.蘆山7.0級地震序列的震源位置與震源機制解特征.地球物理學報, 56(5):1753—1763.

    陸陽泉,錢家棟,劉建毅,1990.大型花崗巖標本緩慢膨脹破裂過程中電阻率和聲發(fā)射前兆特征的實驗研究.西北地震學報,12(2):35—41.

    錢復業(yè),趙玉林,于謀明,1982.地震前地電阻率異常變化.中國科學(B輯),(9):831—839.

    錢復業(yè),趙玉林,劉婕,1990.唐山7.8級地震地電阻率臨震功率譜異常.地震,(3):33—38.

    錢家棟,陳有發(fā),金安忠,1985.地電阻率法在地震預報中的應用.北京:地震出版社,103—107,187—225.

    錢家棟,趙和云,1988.地電阻率的數值模擬和多極距觀測系統(tǒng).地震學報,10(1):77—89.

    錢家棟,曹愛民,1998.1976年唐山7.8級地震地電阻率和地下水前兆綜合物理機制研究.地震,18(增刊):1—9.

    錢家棟,馬欽忠,李劭秾,2013.汶川8.0地震前成都臺NE測線地電阻率異常的進一步研究.地震學報,35(1):4—17.

    王德富,李志明,1985.四川地電阻率預報地震的效能.四川地震,1:31—36.

    王德富,李志明,1990.甘孜臺地電測區(qū)改造和巴塘6.7級地震的預報.四川地震,6:38—42.

    王貴宣,鄭大林,張肇誠等,1995.鮮水河斷裂帶地震活動特征及強震發(fā)生隨時間增長概率.地震研究,18(3):221—226.

    曾祥方,羅艷,韓立波等,2013.2013年4月20日四川蘆山Ms7.0地震:一個高角度逆沖地震.地球物理學報,56(4):1418—1424.

    張國民,馬宏生,王輝等,2005.中國大陸活動地塊邊界帶與強震活動.地球物理學報,48(3):602—610.

    張學民,李美,關華平,2009.汶川8.0級地震前的地電阻率異常分析.地震,29(1):108—115.

    趙和云,錢家棟,1987.郫縣臺多極距觀測資料的反演和分析.中國地震,3(增刊):73—78.

    趙玉林,錢復業(yè),楊休成等,1983.原地電阻率變化的實驗.地震學報,5(2):217—225.

    趙玉林,錢復業(yè),許同春,1990.巖土層受力時電阻率變化與應變的關系.地震學報,12(1):87—93.

    Coggon J.H., 1971. Electromagnetic and electrical modeling by the finite element method. Geophysics, 36: 132—155.

    Dey A., Morrison H.F., 1979. Resistivity modeling for arbitrary shaped three-dimensional structures. Geophysics, 44 (4): 753—780.

    Li Yuguo, Spitzer K., 2005. Finite element resistivity modeling for three-dimensional structures with arbitrary anisotropy. Pepi., 150: 15—27.

    Lowry T., Allen M.B., Shive P.N., 1989. Singularity removal: a refinement of resistivity modeling techniques. Geophysics, 54: 766—774.

    Lu Jun, Qian Fuye, Zhao Yulin, 1999. Sensitivity analysis of the Schlumberger monitoring array: application to changes of resistivity prior to the 1976 earthquake in Tangshan, China. Tectonophy., 307: 397—405.

    Lu Jun, Xue Shunzhang, Qian Fuye et al., 2004. Unexpected changes in resistivity monitoring for earthquakes of the Longmen Shan in Sichuan, China, with a fixed Schlumberger sounding array. Phys. Earth Planet In., 145 (1-4): 87—97.

    O’Neill D.J., Merrick N.P., 1984. A digital linear filter for resistivity sounding with a generalized electrode array. Geophys. Prospect., 32 (1): 105—123.

    Park S.K., Van G.P., 1991. Inversion of pole-pole data for 3-D resistivity structure beneath arrays of electrodes. Geophysics, 56: 951—960.

    Seigel H.O., 1959. Mathematical formulation and type curves for induced polarization. Geophysics, 24: 547—565.

    Wen Xueze, Ma Shilin, Xu Xiwei et al., 2008. Historical pattern and behavior of earthquake ruptures along the eastern boundary of the Sichuan-Yunnan faulted-block, southwestern China. Physics of the Earth and Planetary Interiors, 168 (1-2): 16—36.

    Zhao S., Yedlin M., 1996. Some refinement on the finite-difference method for 3-D dc resistivity modeling. Geophysics, 61: 1301—1307.

    Analysis on Apparent Resistivity Variations of Ganzi Station before 2013 Lushan7.0 Earthquake

    Xie Tao, Lu Jun, Ren Yuexiaand Zhao Minyuan

    1) China Earthquake Networks Center, Beijing 100045, China 2) Sichuan Earthquake Administration, Chengdu 610041, China 3) Ganzi Earthquake Station, Ganzi 626700, China

    Apparent resistivity data of N30oE monitoring direction observed at Ganzi station began to increase since July 2011. And data of N60oW direction also increased since 2012. In 2011 the potential electrodes of N30oE direction moved 10m towards the current electrodes due to the expansion construction of national road No. 317. By using a horizontal layered model our analysis indicates that the potential electrodes shift can cause 4raiseto measurement of N30oE direction. Therefore, apparent resistivity of the two directions increased in the same time of 2012 after offset the effect from electrodes shift. Sensitive coefficients of two directions were also obtained using the model. Sensitive coefficients of shallow layers were minus for the two directions and the unexpected annual variations of Ganzi station were well explained. In order to quantitatively analyze effects on observation from the national road, we constructed an finite element model based on the electric structure. Analysis results also suggested that, the national road expansion construction can only cause 0.15decrease on N60oW monitoring direction and 0.1increase on N30oE direction. Additionally, the valley values of annual variation of 2013 were distinctively higher than other years from 2008, meaning that there existed abnormal rise of apparent resistivity of the two directions at Ganzi station before the Lushan earthquake. However, the rise was contrary to the decline variation before Wenchuan earthquake. Therefore, it was still unconfirmed that whether or not the rise variation was related to Lushan earthquake.

    Lushan earthquake; Ganzi station; Apparent resistivity; Trend increase; Annual variation

    國家科技支撐計劃專題(2012BAK19B02-03)和國家自然科學基金(41204057)共同資助

    2013-07-19

    解滔,男,生于1986年。碩士。主要從事地震電磁學和熱紅外衛(wèi)星遙感應用研究。E-mail:xtaolake@163.com

    解滔,盧軍,任越霞,趙民淵,2013. 2013年蘆山7.0地震前甘孜臺地電阻率變化分析. 震災防御技術,8(4):347—358.

    猜你喜歡
    甘孜蘆山電阻率
    丁真的甘孜,到底有多極致?
    金橋(2021年1期)2021-05-21 08:22:08
    川藏高原甘孜
    三維電阻率成像與高聚物注漿在水閘加固中的應用
    甘孜藏區(qū)中小學體育與健康教育課程教學模式探索
    體育科技(2016年2期)2016-02-28 17:06:14
    隨鉆電阻率測井的固定探測深度合成方法
    海洋可控源電磁場視電阻率計算方法
    春回蘆山
    粉煤灰摻量對水泥漿體電阻率與自收縮的影響
    蘆山地震公路地質災害調查及評估
    災害學(2014年1期)2014-03-01 02:25:58
    蘆山7.0級地震前后巖石圈磁場異常變化研究
    地震研究(2014年1期)2014-02-27 09:29:41
    天美传媒精品一区二区| 美女cb高潮喷水在线观看| 99热只有精品国产| 精品一区二区三区人妻视频| 国产激情偷乱视频一区二区| 三级毛片av免费| 亚洲欧美日韩卡通动漫| 国产精品久久久人人做人人爽| 中文亚洲av片在线观看爽| 免费在线观看日本一区| 成年免费大片在线观看| 日本a在线网址| 日韩欧美三级三区| 久久久久亚洲av毛片大全| 国产探花极品一区二区| 尤物成人国产欧美一区二区三区| 麻豆成人午夜福利视频| 亚洲av一区综合| 欧美午夜高清在线| 中文字幕av成人在线电影| 亚洲最大成人中文| 国内精品久久久久精免费| 免费看日本二区| 精品欧美国产一区二区三| 国产黄色小视频在线观看| 成人国产一区最新在线观看| 亚洲最大成人中文| 国产精品野战在线观看| 法律面前人人平等表现在哪些方面| 亚洲人成电影免费在线| 亚洲精华国产精华精| 露出奶头的视频| 香蕉av资源在线| 少妇丰满av| 亚洲成人免费电影在线观看| 国产亚洲精品一区二区www| 国产高清三级在线| 又爽又黄无遮挡网站| 最近在线观看免费完整版| 丝袜美腿在线中文| 日韩av在线大香蕉| 岛国在线免费视频观看| 欧美成人a在线观看| 国产亚洲欧美在线一区二区| 国内揄拍国产精品人妻在线| 12—13女人毛片做爰片一| 欧美乱色亚洲激情| 国产精品一区二区免费欧美| 成年免费大片在线观看| 给我免费播放毛片高清在线观看| 俄罗斯特黄特色一大片| 国产蜜桃级精品一区二区三区| 精品电影一区二区在线| 国产精品国产高清国产av| 丁香六月欧美| 国产精品永久免费网站| av片东京热男人的天堂| 久久九九热精品免费| 制服丝袜大香蕉在线| 国产欧美日韩精品亚洲av| 啦啦啦免费观看视频1| 国产高潮美女av| 欧美最新免费一区二区三区 | 亚洲av日韩精品久久久久久密| 久久婷婷人人爽人人干人人爱| 国产乱人视频| 午夜两性在线视频| 久久久久免费精品人妻一区二区| 久久人妻av系列| 少妇的丰满在线观看| 亚洲五月婷婷丁香| 精品免费久久久久久久清纯| 岛国视频午夜一区免费看| eeuss影院久久| 欧美一区二区精品小视频在线| 欧美日韩国产亚洲二区| 夜夜躁狠狠躁天天躁| 熟女少妇亚洲综合色aaa.| 18禁黄网站禁片午夜丰满| 性色avwww在线观看| 男女之事视频高清在线观看| 长腿黑丝高跟| 亚洲片人在线观看| 日本五十路高清| 欧美日韩黄片免| 成年版毛片免费区| 久久国产精品人妻蜜桃| 老司机午夜福利在线观看视频| 欧美日韩一级在线毛片| 欧美成人性av电影在线观看| 久久精品综合一区二区三区| 他把我摸到了高潮在线观看| 91av网一区二区| 欧美+亚洲+日韩+国产| 精品一区二区三区视频在线观看免费| 国产精品乱码一区二三区的特点| a级一级毛片免费在线观看| 精品久久久久久久久久久久久| 久久草成人影院| x7x7x7水蜜桃| 国产精品嫩草影院av在线观看 | 久久草成人影院| x7x7x7水蜜桃| 欧美日韩一级在线毛片| 男插女下体视频免费在线播放| 亚洲黑人精品在线| 欧美性感艳星| 90打野战视频偷拍视频| 国产国拍精品亚洲av在线观看 | 夜夜爽天天搞| 国产高清有码在线观看视频| 亚洲精品美女久久久久99蜜臀| 国产主播在线观看一区二区| av国产免费在线观看| 女人被狂操c到高潮| 最近在线观看免费完整版| 午夜免费激情av| 久久天躁狠狠躁夜夜2o2o| 国产高清有码在线观看视频| 国内毛片毛片毛片毛片毛片| 欧美日韩中文字幕国产精品一区二区三区| 88av欧美| 亚洲av美国av| 香蕉av资源在线| 热99在线观看视频| 国产真人三级小视频在线观看| 村上凉子中文字幕在线| 亚洲精品亚洲一区二区| 两个人的视频大全免费| 日韩中文字幕欧美一区二区| 一级黄色大片毛片| 日韩精品中文字幕看吧| 99国产精品一区二区三区| 国产伦一二天堂av在线观看| 亚洲国产日韩欧美精品在线观看 | 高清日韩中文字幕在线| 久久精品国产亚洲av香蕉五月| 国产精品香港三级国产av潘金莲| 精品久久久久久,| 亚洲美女黄片视频| 波多野结衣高清无吗| 国产爱豆传媒在线观看| 91久久精品电影网| 午夜福利高清视频| 国产免费一级a男人的天堂| 亚洲精品在线美女| 精品不卡国产一区二区三区| 亚洲精品一卡2卡三卡4卡5卡| 欧美激情久久久久久爽电影| 久久久久免费精品人妻一区二区| 韩国av一区二区三区四区| 中文字幕高清在线视频| 国产激情偷乱视频一区二区| 一级毛片高清免费大全| 美女被艹到高潮喷水动态| 脱女人内裤的视频| 成人一区二区视频在线观看| 人妻丰满熟妇av一区二区三区| 国产精品永久免费网站| 国产精品永久免费网站| 国内久久婷婷六月综合欲色啪| 日韩欧美在线二视频| 色尼玛亚洲综合影院| 国产高清videossex| 内地一区二区视频在线| 国产精品久久电影中文字幕| 亚洲精品乱码久久久v下载方式 | ponron亚洲| 99久久久亚洲精品蜜臀av| 国产爱豆传媒在线观看| 好男人电影高清在线观看| 好男人电影高清在线观看| 日韩欧美国产一区二区入口| 久久精品人妻少妇| 神马国产精品三级电影在线观看| 好男人电影高清在线观看| 一进一出抽搐gif免费好疼| 午夜福利高清视频| 国产男靠女视频免费网站| 此物有八面人人有两片| 2021天堂中文幕一二区在线观| 国产男靠女视频免费网站| 亚洲精品日韩av片在线观看 | 女生性感内裤真人,穿戴方法视频| 最近视频中文字幕2019在线8| 欧美xxxx黑人xx丫x性爽| 黄色片一级片一级黄色片| 亚洲第一电影网av| 在线国产一区二区在线| 一本综合久久免费| 欧美不卡视频在线免费观看| 国产精品98久久久久久宅男小说| 不卡一级毛片| 国产男靠女视频免费网站| 一级a爱片免费观看的视频| 久久亚洲精品不卡| 国产欧美日韩精品亚洲av| 欧美一区二区亚洲| 一个人免费在线观看电影| 欧美成人a在线观看| 午夜精品久久久久久毛片777| 成年免费大片在线观看| 国产精品久久电影中文字幕| 国产在线精品亚洲第一网站| 欧美高清成人免费视频www| 91九色精品人成在线观看| 亚洲国产精品合色在线| 又紧又爽又黄一区二区| 亚洲 欧美 日韩 在线 免费| 三级国产精品欧美在线观看| 国产三级中文精品| 美女免费视频网站| 国产精品av视频在线免费观看| 他把我摸到了高潮在线观看| 少妇高潮的动态图| 九九在线视频观看精品| 亚洲国产高清在线一区二区三| 国产高清视频在线播放一区| 国产精品影院久久| 精华霜和精华液先用哪个| 国产一区二区在线观看日韩 | 国产精品三级大全| 18禁在线播放成人免费| 51午夜福利影视在线观看| 老鸭窝网址在线观看| 亚洲成av人片免费观看| 日韩有码中文字幕| 久久精品国产99精品国产亚洲性色| 精品日产1卡2卡| 男女那种视频在线观看| 99久久精品一区二区三区| 久久久久免费精品人妻一区二区| 日韩精品中文字幕看吧| 日本黄色视频三级网站网址| 大型黄色视频在线免费观看| 精品久久久久久,| 亚洲成人精品中文字幕电影| 久久精品国产亚洲av香蕉五月| 亚洲精品美女久久久久99蜜臀| 精品久久久久久久毛片微露脸| 亚洲七黄色美女视频| 久久久久国内视频| 亚洲人与动物交配视频| 亚洲第一欧美日韩一区二区三区| 国产精品永久免费网站| av专区在线播放| 久久精品国产亚洲av涩爱 | 男女之事视频高清在线观看| 国产精品久久视频播放| 一级作爱视频免费观看| www国产在线视频色| 中文字幕人妻熟人妻熟丝袜美 | 在线免费观看的www视频| 亚洲无线在线观看| 夜夜爽天天搞| 久99久视频精品免费| 成熟少妇高潮喷水视频| 午夜精品一区二区三区免费看| 一本综合久久免费| 亚洲人成网站在线播放欧美日韩| 在线观看av片永久免费下载| 美女大奶头视频| 美女高潮喷水抽搐中文字幕| 天天添夜夜摸| 成人无遮挡网站| 美女 人体艺术 gogo| 日韩欧美在线二视频| 免费一级毛片在线播放高清视频| 色老头精品视频在线观看| av在线蜜桃| 18禁黄网站禁片免费观看直播| 亚洲av免费在线观看| 99国产综合亚洲精品| 久久久久九九精品影院| 99久久精品国产亚洲精品| 午夜久久久久精精品| 亚洲国产欧美网| 久久国产精品影院| 午夜福利在线在线| 国产91精品成人一区二区三区| 最近视频中文字幕2019在线8| 一边摸一边抽搐一进一小说| 国产日本99.免费观看| 欧美大码av| 成人永久免费在线观看视频| 亚洲精品在线美女| 国产一区在线观看成人免费| 欧美黑人欧美精品刺激| 亚洲欧美一区二区三区黑人| 国产淫片久久久久久久久 | 国产精品女同一区二区软件 | 又紧又爽又黄一区二区| 亚洲精品一区av在线观看| 国产精品久久久人人做人人爽| 久久久精品欧美日韩精品| 十八禁人妻一区二区| 欧美中文日本在线观看视频| 看免费av毛片| 国语自产精品视频在线第100页| 91麻豆精品激情在线观看国产| 俺也久久电影网| 99久久99久久久精品蜜桃| 久久久成人免费电影| 国产精品精品国产色婷婷| 18禁裸乳无遮挡免费网站照片| 成年免费大片在线观看| 天堂网av新在线| 国产黄片美女视频| av天堂在线播放| 天堂动漫精品| 老鸭窝网址在线观看| 成年免费大片在线观看| 免费人成在线观看视频色| 久久国产精品影院| 美女免费视频网站| 嫩草影院精品99| 成人av在线播放网站| 亚洲国产精品久久男人天堂| 久久精品亚洲精品国产色婷小说| 日韩欧美 国产精品| 国产日本99.免费观看| 看免费av毛片| 麻豆成人av在线观看| 亚洲乱码一区二区免费版| 亚洲熟妇中文字幕五十中出| 亚洲五月天丁香| 亚洲第一欧美日韩一区二区三区| 天堂av国产一区二区熟女人妻| 99视频精品全部免费 在线| 3wmmmm亚洲av在线观看| 成人特级av手机在线观看| 亚洲成av人片免费观看| 丰满人妻熟妇乱又伦精品不卡| 无限看片的www在线观看| 日韩人妻高清精品专区| 美女高潮喷水抽搐中文字幕| 欧美黑人巨大hd| av片东京热男人的天堂| 国产欧美日韩精品一区二区| 国产精品久久久久久久电影 | 欧美色欧美亚洲另类二区| 脱女人内裤的视频| 亚洲 欧美 日韩 在线 免费| 久久这里只有精品中国| 日日干狠狠操夜夜爽| av片东京热男人的天堂| 天堂av国产一区二区熟女人妻| 国产99白浆流出| 久久精品国产清高在天天线| 欧美在线一区亚洲| 久久草成人影院| 一本综合久久免费| 国产av在哪里看| 亚洲av电影在线进入| 99精品久久久久人妻精品| 淫妇啪啪啪对白视频| 夜夜看夜夜爽夜夜摸| 他把我摸到了高潮在线观看| 欧美zozozo另类| 欧美成人免费av一区二区三区| 国产伦在线观看视频一区| 久久久久亚洲av毛片大全| 欧美日韩中文字幕国产精品一区二区三区| 精品欧美国产一区二区三| 日本一本二区三区精品| h日本视频在线播放| 欧美bdsm另类| 女人十人毛片免费观看3o分钟| 少妇的丰满在线观看| 欧美日韩国产亚洲二区| 国产av在哪里看| 最近视频中文字幕2019在线8| 国产麻豆成人av免费视频| 国产精品电影一区二区三区| av在线蜜桃| 久久精品影院6| 好男人在线观看高清免费视频| 国产精品一区二区免费欧美| 亚洲成人免费电影在线观看| 他把我摸到了高潮在线观看| 中国美女看黄片| 国产高清videossex| 在线观看66精品国产| 国产伦精品一区二区三区视频9 | 色视频www国产| 久久久久免费精品人妻一区二区| 女人高潮潮喷娇喘18禁视频| 欧美性猛交╳xxx乱大交人| 国产亚洲欧美98| 又黄又粗又硬又大视频| 亚洲国产精品999在线| 香蕉丝袜av| 蜜桃亚洲精品一区二区三区| 伊人久久精品亚洲午夜| av国产免费在线观看| 欧美zozozo另类| 亚洲精品乱码久久久v下载方式 | 国产黄色小视频在线观看| 国产v大片淫在线免费观看| 亚洲国产精品合色在线| 亚洲片人在线观看| 久久久久国内视频| 日本一本二区三区精品| 1000部很黄的大片| 色综合站精品国产| 亚洲成a人片在线一区二区| 欧美乱色亚洲激情| 美女免费视频网站| 国产精品 国内视频| 制服人妻中文乱码| 麻豆久久精品国产亚洲av| 琪琪午夜伦伦电影理论片6080| 99热6这里只有精品| 在线播放无遮挡| 一级a爱片免费观看的视频| 亚洲国产精品成人综合色| 尤物成人国产欧美一区二区三区| 一区二区三区激情视频| 老汉色av国产亚洲站长工具| 成人性生交大片免费视频hd| 99久久99久久久精品蜜桃| 久久久久久九九精品二区国产| 国产一区二区三区在线臀色熟女| 亚洲一区高清亚洲精品| www日本在线高清视频| 色综合亚洲欧美另类图片| 亚洲成av人片免费观看| 免费一级毛片在线播放高清视频| 人人妻,人人澡人人爽秒播| 男女之事视频高清在线观看| 国产精品久久电影中文字幕| 中文亚洲av片在线观看爽| 国产国拍精品亚洲av在线观看 | 两个人看的免费小视频| 国产午夜精品论理片| 97人妻精品一区二区三区麻豆| 欧美成人a在线观看| 日本一本二区三区精品| 又粗又爽又猛毛片免费看| 国产中年淑女户外野战色| 精品国内亚洲2022精品成人| 人人妻,人人澡人人爽秒播| 亚洲av熟女| 国产高清三级在线| 欧美+亚洲+日韩+国产| 午夜亚洲福利在线播放| 久久伊人香网站| 精品一区二区三区人妻视频| 国产精品爽爽va在线观看网站| 99久久综合精品五月天人人| 一级黄色大片毛片| 麻豆国产97在线/欧美| 首页视频小说图片口味搜索| 黄色日韩在线| 国产av不卡久久| 国产单亲对白刺激| 麻豆成人av在线观看| 国产一区二区亚洲精品在线观看| 首页视频小说图片口味搜索| 一级黄色大片毛片| 国产精品日韩av在线免费观看| 亚洲精品456在线播放app | 18禁美女被吸乳视频| 欧美大码av| 午夜福利在线在线| 欧美成人a在线观看| 可以在线观看毛片的网站| 美女大奶头视频| 午夜福利视频1000在线观看| 18+在线观看网站| 久久婷婷人人爽人人干人人爱| 国产真实乱freesex| 成人高潮视频无遮挡免费网站| 美女高潮喷水抽搐中文字幕| aaaaa片日本免费| 变态另类成人亚洲欧美熟女| 国产亚洲欧美在线一区二区| 欧美另类亚洲清纯唯美| 两个人的视频大全免费| а√天堂www在线а√下载| 国产成人av激情在线播放| 国产午夜精品久久久久久一区二区三区 | 色综合欧美亚洲国产小说| 亚洲一区二区三区色噜噜| 制服丝袜大香蕉在线| 国产在线精品亚洲第一网站| 亚洲美女黄片视频| 国产亚洲精品久久久com| 欧美+日韩+精品| 中文字幕人妻丝袜一区二区| 久久精品国产99精品国产亚洲性色| 成人特级黄色片久久久久久久| 极品教师在线免费播放| 观看免费一级毛片| 国产成人系列免费观看| 男人和女人高潮做爰伦理| 亚洲欧美日韩东京热| 国产成人福利小说| 亚洲国产中文字幕在线视频| 精品人妻偷拍中文字幕| 精品久久久久久久末码| 欧美日本视频| 黄片小视频在线播放| 日本黄色片子视频| 亚洲中文字幕一区二区三区有码在线看| 国产精品一区二区三区四区久久| 青草久久国产| www日本在线高清视频| 亚洲精品成人久久久久久| 夜夜爽天天搞| 久久精品国产综合久久久| 国产精品,欧美在线| 女警被强在线播放| 国产老妇女一区| 少妇裸体淫交视频免费看高清| 高潮久久久久久久久久久不卡| 国产不卡一卡二| 免费在线观看亚洲国产| 亚洲精品在线美女| 91久久精品国产一区二区成人 | 亚洲 欧美 日韩 在线 免费| 国产伦人伦偷精品视频| av专区在线播放| 91麻豆精品激情在线观看国产| 欧美黄色淫秽网站| 久久久久久人人人人人| 久久久久国内视频| 天天一区二区日本电影三级| 久久久久久久午夜电影| 观看免费一级毛片| 国产精品香港三级国产av潘金莲| 亚洲人成伊人成综合网2020| 亚洲 欧美 日韩 在线 免费| 日韩av在线大香蕉| 亚洲一区二区三区不卡视频| 国产黄a三级三级三级人| 天堂影院成人在线观看| 欧美bdsm另类| 夜夜夜夜夜久久久久| 久久久久久九九精品二区国产| a在线观看视频网站| 少妇高潮的动态图| 免费大片18禁| 天堂av国产一区二区熟女人妻| 黄片小视频在线播放| 亚洲精品乱码久久久v下载方式 | 51国产日韩欧美| 亚洲一区二区三区色噜噜| 国产精品久久视频播放| 欧美精品啪啪一区二区三区| 亚洲狠狠婷婷综合久久图片| 日日夜夜操网爽| 国产精品免费一区二区三区在线| 精品久久久久久成人av| 亚洲欧美日韩无卡精品| x7x7x7水蜜桃| 亚洲专区中文字幕在线| 久久午夜亚洲精品久久| 精品不卡国产一区二区三区| 欧美在线一区亚洲| 一个人看视频在线观看www免费 | 精品国内亚洲2022精品成人| 日韩欧美三级三区| 级片在线观看| 久久这里只有精品中国| 婷婷丁香在线五月| 亚洲成人精品中文字幕电影| 免费在线观看亚洲国产| 国产乱人视频| 在线观看日韩欧美| 麻豆久久精品国产亚洲av| 亚洲五月婷婷丁香| 欧美区成人在线视频| 中亚洲国语对白在线视频| 手机成人av网站| 欧美激情在线99| 亚洲国产精品合色在线| 国产黄色小视频在线观看| 国产亚洲精品一区二区www| 亚洲av电影在线进入| 一区二区三区国产精品乱码| 国产精品三级大全| 人妻丰满熟妇av一区二区三区| 搞女人的毛片| 极品教师在线免费播放| 国产亚洲精品综合一区在线观看| 丰满乱子伦码专区| 麻豆一二三区av精品| 国产一区二区激情短视频| 窝窝影院91人妻| 麻豆国产97在线/欧美| 啦啦啦韩国在线观看视频| 五月玫瑰六月丁香| 国产成人a区在线观看| 亚洲av美国av| 国产成人av激情在线播放| 中文字幕人妻丝袜一区二区| 国产高清三级在线| 日本精品一区二区三区蜜桃| 亚洲在线观看片| 天堂√8在线中文| 国产激情欧美一区二区| 有码 亚洲区| 亚洲中文日韩欧美视频| 综合色av麻豆| 韩国av一区二区三区四区| 深爱激情五月婷婷| 久久精品国产亚洲av涩爱 | 波野结衣二区三区在线 | 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 免费无遮挡裸体视频| 男女视频在线观看网站免费| 波野结衣二区三区在线 | 婷婷丁香在线五月| 一区福利在线观看|