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

    吉林省地電阻率震兆異常分析研究*

    2014-08-02 03:55:03史紅軍趙衛(wèi)星張可佳朱偉楠
    地震學(xué)報(bào) 2014年3期

    史紅軍 趙衛(wèi)星 張可佳 朱偉楠 李 寧

    1) 中國(guó)長(zhǎng)春130400榆樹地震臺(tái) 2) 中國(guó)長(zhǎng)春130012長(zhǎng)春市地震速測(cè)速報(bào)中心

    吉林省地電阻率震兆異常分析研究*

    1) 中國(guó)長(zhǎng)春130400榆樹地震臺(tái) 2) 中國(guó)長(zhǎng)春130012長(zhǎng)春市地震速測(cè)速報(bào)中心

    應(yīng)用吉林省內(nèi)榆樹、 四平兩個(gè)地電阻率臺(tái)站觀測(cè)資料,采用形態(tài)法、 歸一化變化速率法和各向異性度法,分析研究了1999年汪清MS7.0深震以及2006年乾安—前郭MS5.0淺源地震前的地電阻率異常變化. 結(jié)果表明, 震前地電阻率出現(xiàn)長(zhǎng)趨勢(shì)下降變化或短臨異常變化,表現(xiàn)為: ① 汪清MS7.0深震前,榆樹臺(tái)和四平臺(tái)地電阻率月均值曲線上均顯示出震前持續(xù)2年尺度的中期下降異常; ② 汪清MS7.0深震以及乾安—前郭MS5.0淺源中強(qiáng)震前,榆樹臺(tái)和四平臺(tái)地電阻率變化速率均大于異常指標(biāo); ③ 乾安—前郭MS5.0地震前,四平臺(tái)地電阻率各向異性度曲線顯示出明顯的低值異常變化,而榆樹臺(tái)各向異性度S曲線則出現(xiàn)破年變變化. 上述異常變化的原因可能為,太平洋板塊向NW方向擠壓,導(dǎo)致震源區(qū)及附近最大主壓應(yīng)力方位(或近于該方位)NW向的擠壓作用突出,引起介質(zhì)內(nèi)部導(dǎo)電流體快速進(jìn)入或重新分布,從而使地電阻率出現(xiàn)長(zhǎng)趨勢(shì)下降變化和各向異性變化.

    地電阻率 地震 震兆異常 異常形態(tài) 歸一化變化速率 各向異性度

    引言

    從1966年邢臺(tái)地震開始,40多年的地震預(yù)測(cè)預(yù)報(bào)實(shí)踐表明,在多次大地震、 中等地震前均記錄到了地電阻率異常變化. 趙玉林和錢復(fù)業(yè)(1978)、 錢復(fù)業(yè)和趙玉林(1980)報(bào)道1976年唐山MS7.8強(qiáng)震和其它一些中小地震前震中及其附近地區(qū)電阻率下降, 并在震中形成一個(gè)電阻率下降的異常區(qū). Qian等(1996)在唐山震例研究中得到,震前二三年多臺(tái)出現(xiàn)地電阻率同步下降趨勢(shì)異常. 關(guān)華平等(2008)利用四川及其周邊省市的地電阻率觀測(cè)資料分析2008年汶川地震前的地電阻率異?,F(xiàn)象時(shí),發(fā)現(xiàn)這次MS8.0地震前,距震中310 km范圍內(nèi)的5個(gè)地電阻率臺(tái)站中,有3個(gè)臺(tái)站出現(xiàn)2年以上趨勢(shì)下降異常. 杜學(xué)彬等(2001)提出了歸一化月速率方法,并用該方法處理了中等以上地震前近震中區(qū)地電阻率數(shù)據(jù),提取出了地震地電阻率各向異性變化的前兆異常. 毛桐恩等(1999)提出了無(wú)量綱地電阻率各向異性度(S)方法,并將其應(yīng)用于地震預(yù)報(bào)實(shí)踐; 此后,其他學(xué)者也進(jìn)行了大量震前地電阻率各向異性變化方面的研究,積累了不少震例(馮志生等,2004; 阮愛國(guó)等,2004).

    圖1 吉林省地電阻率臺(tái)站與震中分布圖Fig.1 Distribution of earth resistivity observation stations and epicenters of Jilin Province

    本文在考察吉林省地電阻率臺(tái)站的觀測(cè)環(huán)境、 觀測(cè)儀器、 觀測(cè)資料以及干擾源等前提下,在前人研究基礎(chǔ)上,利用榆樹和四平這兩個(gè)地電阻率臺(tái)站的觀測(cè)數(shù)據(jù),采用形態(tài)法、 地電阻率各向異性度法和歸一化變化速率法,研究了區(qū)域(121.60°—131.40°E,40.75°—46.35°N)內(nèi)發(fā)生的MS≥4.3 地震前各臺(tái)站觀測(cè)到的地電阻率震兆異常,期望為今后地電阻率觀測(cè)資料的分析和應(yīng)用,尤其是吉林省地震地電阻率異常的正確判定提供一定的依據(jù).

    1 吉林省地電阻率臺(tái)站概況

    目前,吉林省內(nèi)有3個(gè)地電阻率觀測(cè)臺(tái)站,分別為四平臺(tái)、 白城臺(tái)和榆樹臺(tái). 各臺(tái)站及震中分布如圖1所示. 因文中未使用白城臺(tái)資料,故僅對(duì)四平臺(tái)和榆樹臺(tái)的概況予以介紹.

    1.1 四平臺(tái)概況

    四平臺(tái)位于東北斷塊區(qū)松遼沉降坳陷帶與張廣才嶺—老爺嶺斷塊隆起帶的交會(huì)部位,靠近松遼沉降坳陷帶東部斜坡帶的邊緣. 其附近有NE向伊通—舒蘭斷裂帶經(jīng)過(guò),基底為海西期花崗巖,埋深約500 m. 地表出露第四系黃土層,基巖為白堊系泉頭統(tǒng)砂巖、 粉砂巖. 地下水類型為基巖裂隙水. 地表下11 m為第四系黏土、 亞黏土; 11—96 m為白堊系粉砂巖,泥巖互層; 96 m以下基底為海西期花崗巖. 四平臺(tái)地處吉林省四平—長(zhǎng)春—榆樹中強(qiáng)地震活動(dòng)帶,屬小震活動(dòng)頻發(fā)區(qū)(吉林省地震局,2005).

    四平臺(tái)地電阻率于1977年4月開始觀測(cè),觀測(cè)儀器為ZD8B,布設(shè)N40°E 和N50°W共兩個(gè)測(cè)道,供電極距AB=1100 m,測(cè)量極距MN=300 m. 供電/測(cè)量電極采用厚5 cm、 100 cm×100 cm正方形鉛板電極,電極埋深3.0 m. 供電電極接地電阻為 4.0—26.0 Ω,測(cè)量電極接地電阻為4.0—20.0 Ω. 2003年11月以前,外線路采用架空方式,供電線、 測(cè)量線均使用2×1.5 mm2電纜; 2003年11月以后,外線路采用地埋方式,供電線、 測(cè)量線均使用2×2.5 mm2鍍鋅鎧裝電纜,外線路絕緣≥50 MΩ. 該臺(tái)裝置系統(tǒng)符合DB/T 18.1—2006關(guān)于臺(tái)站建設(shè)的技術(shù)要求(杜學(xué)彬等,2006).

    1.2 榆樹臺(tái)概況

    榆樹臺(tái)所處大地構(gòu)造位置為東北斷塊區(qū)松遼斷陷沉降帶的東部隆起區(qū),位于NE向伊通—舒蘭、 四平—長(zhǎng)春兩個(gè)深大斷裂帶北延部分的中間地帶、 NNW向卡岔河斷裂西側(cè)2.5 km的平原地區(qū). 沿NE向四平—長(zhǎng)春斷裂帶及伊通—舒蘭斷裂帶上曾發(fā)生過(guò)中等強(qiáng)度的地震,現(xiàn)代中小震活動(dòng)明顯. 該臺(tái)站位于NE和NNW向3條活動(dòng)斷裂交會(huì)部位附近. 臺(tái)站東側(cè)的NNW向卡岔河斷裂構(gòu)造比較發(fā)育,臺(tái)址所處的地質(zhì)背景對(duì)捕捉地震前兆信息十分有利. 榆樹臺(tái)臺(tái)址基巖主要是白堊系砂巖, 第四紀(jì)覆蓋大多在70—90 m,巖石結(jié)構(gòu)致密,對(duì)監(jiān)測(cè)明顯地震前兆信息是有利的(史紅軍等,2011).

    榆樹臺(tái)地電阻率于1993年1月正式觀測(cè),觀測(cè)儀器為ZD8BI,布設(shè)N45°E 和N45°W兩個(gè)測(cè)道,供電極距AB=900 m,測(cè)量極距MN=300 m. 供電/測(cè)量電極采用厚6 cm、 100 cm×100 cm正方形鉛板電極,電極埋深2.5 m. 供電電極接地電阻為 6.0—27.0 Ω,測(cè)量電極接地電阻為5.5—27.5 Ω. 2003年11月以前,外線路采用架空方式,供電線使用2×1.5 mm2電纜,測(cè)量線1997年4月以前使用被覆線,1997年5月—2003年10月,使用2×1.5 mm2電纜; 2003年11月以后,外線路采用地埋方式,供電線、 測(cè)量線均使用2×2.5 mm2鍍鋅鎧裝電纜,外線路絕緣≥50 MΩ. 該臺(tái)裝置系統(tǒng)符合DB/T 18.1—2006關(guān)于臺(tái)站建設(shè)的技術(shù)要求(杜學(xué)彬等,2006).

    1.3 臺(tái)站運(yùn)行觀測(cè)情況

    對(duì)四平臺(tái)、 榆樹臺(tái)建臺(tái)以來(lái)數(shù)據(jù)統(tǒng)計(jì)結(jié)果表明,四平臺(tái)、 榆樹臺(tái)的裝置系統(tǒng)穩(wěn)定性均滿足《地震及前兆數(shù)字觀測(cè)技術(shù)規(guī)范——電磁觀測(cè)(試行)》(以下簡(jiǎn)稱《規(guī)范》, 中國(guó)地震局,2001)技術(shù)要求: ① 測(cè)量極接地電阻<100 Ω; ② 供電極接地電阻<30 Ω; ③ 漏電電流影響系數(shù)ε1<0.1%; ④ 漏電電位差影響系數(shù)ε2<0.5%.

    地電阻率觀測(cè)精度由觀測(cè)系統(tǒng)誤差和偶然誤差兩部分組成(中國(guó)地震局,2001). 觀測(cè)系統(tǒng)誤差主要是測(cè)量?jī)x器誤差,一般可通過(guò)對(duì)儀器的定期標(biāo)定予以克服; 偶然誤差則是多種因素的綜合作用使觀測(cè)值產(chǎn)生某些變化. 偶然誤差使測(cè)量值在均值附近波動(dòng),其大小用日均方誤差的月均值來(lái)衡量,一般取其相對(duì)均方根誤差Kσ. 從四平臺(tái)、 榆樹臺(tái)儀器標(biāo)定統(tǒng)計(jì)結(jié)果來(lái)看,建臺(tái)以來(lái),這兩個(gè)臺(tái)站的儀器標(biāo)定結(jié)果均合格. 多年來(lái),這兩個(gè)臺(tái)的地電儀器,尤其是數(shù)字地電儀器產(chǎn)生的地電阻率相對(duì)均方根誤差Kσ的月均值均在0.3%以內(nèi)變化,保證了可靠識(shí)別不低于1%的地電阻率異常.

    2 震兆異常分析

    吉林省屬于少震省份,省內(nèi)自1960年榆樹土橋MS5.8地震、 1966年懷德范家屯MS5.2地震后,直至2006年3月31日乾安—前郭MS5.0地震發(fā)生,5級(jí)以上淺源中強(qiáng)地震發(fā)生的時(shí)間跨度達(dá)40年; 而處于我國(guó)深源地震區(qū)的汪清,于1999年、 2002年先后發(fā)生MS7.0、MS7.2深源地震(史紅軍等,2011). 經(jīng)整理吉林省地電阻率臺(tái)站觀測(cè)數(shù)據(jù),采用形態(tài)法、 歸一化變化速率法和各向異性度法深入分析研究得到的結(jié)果表明,在汪清MS≥7.0深源地震以及乾安—前郭MS5.0等中強(qiáng)地震前,均記錄到了較為明顯的地電阻率震兆異常.

    2.1 地電阻率形態(tài)異常分析

    形態(tài)分析是根據(jù)原始資料的變化形態(tài)進(jìn)行分析的一種常用方法. 形態(tài)法由于沒(méi)有作任何數(shù)據(jù)處理,因此能較好反映震前各個(gè)時(shí)段地電阻率的原始形態(tài),直觀、 簡(jiǎn)單,物理意義明確. 1978年唐山MS7.8地震前河北昌黎臺(tái)觀測(cè)的地電阻率出現(xiàn)大幅度下降性異常變化; 同年8月四川松潘MS7.2地震前甘肅武都臺(tái)觀測(cè)的地電阻率也出現(xiàn)了下降性異常變化,地震發(fā)生時(shí)異常幅度達(dá)到最大,地震后逐漸恢復(fù)(錢家棟,1993). 杜學(xué)彬等(1999,2000)應(yīng)用歸一化變化速率法處理了我國(guó)100多個(gè)地電臺(tái)站資料. 統(tǒng)計(jì)表明,MS7.0—7.8地震附近約150 km范圍內(nèi)出現(xiàn)異常的臺(tái)站全部為下降性異常;MS6.0—6.9地震約50 km范圍內(nèi)出現(xiàn)異常的臺(tái)站也全部為下降性異常. 隨著震中距增加,下降性異常占異常總數(shù)目的比例減小,MS5.0—5.9地震出現(xiàn)的下降性異常仍為多數(shù)(張繼紅等,2010).

    應(yīng)用形態(tài)法分析榆樹臺(tái)和四平臺(tái)地電阻率觀測(cè)數(shù)據(jù),可以看出汪清MS7.0深震前,兩個(gè)臺(tái)站地電阻率月均值曲線均存在長(zhǎng)趨勢(shì)下降. 汪清MS7.0深震前,榆樹臺(tái)兩測(cè)向地電阻率月均值曲線從1997年1月初開始幾乎同步出現(xiàn)長(zhǎng)趨勢(shì)下降(圖2); 震后,兩測(cè)向地電阻率月均值曲線均緩慢回升(史紅軍等,2008). 這種異常變化,與唐山、 松潘以及汶川大地震的震例總體異常特征相似(趙玉林,錢復(fù)業(yè),1978; 錢復(fù)業(yè),趙玉林,1980; 張學(xué)民等,2009). 盡管兩測(cè)向地電阻率的變化幅度均小于1%,但均顯示出2年尺度的中期異常. 四平臺(tái)N50°W向地電阻率月均值曲線從1997年1月初出現(xiàn)長(zhǎng)趨勢(shì)下降,到1998年11月降至最低點(diǎn)后開始轉(zhuǎn)折加速上升,上升至最高點(diǎn)時(shí)發(fā)生MS7.0深震. 該曲線顯現(xiàn)了一個(gè)完整的下降—上升異常變化,與杜學(xué)彬等(2000)的研究成果相一致. 雖然異常變化幅度只有0.47%,但1—2年尺度中期異常較明顯.

    2.2 地電阻率歸一化變化速率異常分析

    歸一化變化速率法為杜學(xué)彬等(2001)提出的一種提取地電阻率中短期至短期短臨異常的方法,杜學(xué)彬(2010)進(jìn)一步完善了該方法. 其數(shù)學(xué)定義為

    (1)

    圖2 榆樹臺(tái)和四平臺(tái)地電阻率月均值曲線Fig.2 The curves of monthly mean values of earth resistivity recorded at Yushu and Siping stations

    (2)

    圖3與圖4分別給出了榆樹臺(tái)和四平臺(tái)1997—2002年及2003—2011年地電阻率歸一化變化速率曲線,滑動(dòng)步長(zhǎng)為8. 由圖可見,在1999年汪清MS7.0、 2002年汪清MS7.2、 2011年琿春MS6.1深震以及2006年乾安—前郭MS5.0、 2009年伊通—公主嶺MS4.3淺源中強(qiáng)地震前,榆樹臺(tái)、 四平臺(tái)地電阻率歸一化變化速率均出現(xiàn)不同程度的異常變化. 對(duì)于MS6.1以上深震而言,地電阻率歸一化月速率值于震前3—28個(gè)月出現(xiàn)異常. 例如,汪清MS7.2深震前,榆樹臺(tái)N45°E測(cè)向地電阻率歸一化月速率值于2002年4—6月出現(xiàn)最大值為6.04的短臨異常,N45°W測(cè)向地電阻率歸一化月速率值于2000年3—6月出現(xiàn)最大值7.29后,于2001年3—6月又出現(xiàn)最小值為-6.77的中期異常; 四平臺(tái)N50°W測(cè)向地電阻率歸一化月速率值于2001年5、 6月出現(xiàn)最小值為-3.28的中期異常. 而對(duì)于淺源中強(qiáng)地震而言,地電阻率歸一化月速率值于震前3—28個(gè)月出現(xiàn)異常. 例如,乾安—前郭MS5.0地震前,四平臺(tái)N40°E測(cè)向地電阻率歸一化月速率值于2005年7—9月出現(xiàn)近1年尺度的高值異常變化后,又于2006年1—2月出現(xiàn)短臨低值異常; 其N50°W測(cè)向地電阻率歸一化月速率值于2003年12月—2004年1月出現(xiàn)2年尺度的低值異常.

    圖3 1997—2002年榆樹臺(tái)和四平臺(tái)地電阻率歸一化變化速率曲線

    圖4 2003—2011年榆樹臺(tái)和四平臺(tái)地電阻率歸一化變化速率曲線

    2.3 地電阻率各向異性度異常分析

    (3)

    式中,ρSN和ρEW為NS向和EW向的地電阻率月均值,n為月數(shù)(n≥7). 在孕震等周圍環(huán)境因素?zé)o變化時(shí),S通常為常數(shù). 若兩個(gè)測(cè)量方向的觀測(cè)值ρSN與ρEW之差增大,表明巖石介質(zhì)的各向異性增強(qiáng); 反之,則表明介質(zhì)的各向異性減弱. 即S值的大小及變化,表征了孕震過(guò)程中巖石介質(zhì)各向異性的強(qiáng)弱及變化,或者巖石破裂程度的變化. 其實(shí)質(zhì)反映了孕震構(gòu)造應(yīng)力場(chǎng)的動(dòng)態(tài)變化(陸陽(yáng)泉等,1998; 錢復(fù)業(yè),趙玉林,1998; 毛桐恩等,1999).

    根據(jù)式(3),分別對(duì)四平臺(tái)、 榆樹臺(tái)近年地電阻率進(jìn)行計(jì)算,得到了四平臺(tái)和榆樹臺(tái)的地電阻率各向異性度S曲線(圖5,6).

    圖5 2004—2009年四平臺(tái)電阻率月均值曲線及各向異性度(S)曲線

    圖6 2004—2009年榆樹臺(tái)電阻率月均值曲線及其各向異性度(S)曲線

    四平臺(tái)地電阻率各向異性度值從2005年6月份開始出現(xiàn)加速下降,至10—11月各向異性度出現(xiàn)最低值,在各向異性度值恢復(fù)上升過(guò)程中,發(fā)生了2006年3月乾安—前郭MS5.0地震. 同樣,其各向異性度值在2008年6月再次開始出現(xiàn)加速下降,至8—9月份各向異性度出現(xiàn)最低值,然后轉(zhuǎn)折上升,在此過(guò)程中,發(fā)生了伊通—公主嶺MS4.3地震. 在這兩次地震前,四平臺(tái)地電阻率各向異性度曲線均顯現(xiàn)出半年尺度的短臨異常.

    由圖6可見,榆樹臺(tái)地電阻率各向異性度曲線變化比較有規(guī)律,呈現(xiàn)類似于正弦波變化,但在2006年3月乾安—前郭MS5.0地震前,這一規(guī)律被打破,各向異性度曲線在2005年8月—2006年5月呈近乎直線變化,從2005年8月份出現(xiàn)破年變變化,到地震發(fā)生,持續(xù)時(shí)間為8個(gè)月; 地震發(fā)生后,各向異性度曲線逐漸恢復(fù)正常.

    3 地電阻率異常變化特征

    研究結(jié)果表明,大震地電阻率異常范圍,其半徑大于200 km、 小于700 km,異常時(shí)間在12—20個(gè)月內(nèi)的發(fā)震概率最高,為97.8%; 強(qiáng)震地電阻率異常范圍大多在300 km以內(nèi),趨勢(shì)異常時(shí)間在12個(gè)月內(nèi)發(fā)震的概率最高,為86.5%; 中強(qiáng)地震地電阻率異常范圍,半徑一般在200 km以內(nèi),趨勢(shì)異常時(shí)間在180天之內(nèi)的發(fā)震概率最高,為95.0%(汪志亮等,1995). 在吉林省境內(nèi)發(fā)生的MS≥6.1深源地震及MS≥4.3淺源中強(qiáng)地震前,地電阻率能夠較好地反映震前異常變化,異常特征及地震的各項(xiàng)參數(shù)如表1所示. 從表1可以清楚地看到,映震范圍為50—550 km; 異常出現(xiàn)至發(fā)震時(shí)間間隔為90—840天; 異常持續(xù)時(shí)間為30—730天; 異常變化幅度為0.47%—0.89%. 異常變化形態(tài)主要以地電阻率月均值曲線呈長(zhǎng)趨勢(shì)下降、 歸一化速率曲線出現(xiàn)高值或低值、 各向異性度曲線出現(xiàn)低值等異常形式出現(xiàn).

    4 地電阻率異常變化的物理解釋

    吉林省的汪清、 琿春至黑龍江省的東寧、 穆棱一帶(129°—132°E,40.5°—44.5°N),是我國(guó)唯一的深源地震區(qū). 張立敏和唐曉明(1983)研究了西太平洋板塊俯沖運(yùn)動(dòng)對(duì)東北深震帶的影響,認(rèn)為東北地區(qū)深震是西太平洋板塊俯沖到中國(guó)大陸之下造成的. 孟憲森等(1996)研究認(rèn)為,深震是地幔流驅(qū)動(dòng)太平洋板塊向歐亞大陸板塊俯沖時(shí)下插板塊發(fā)生斷裂所致, 且深震活動(dòng)與東北淺震活動(dòng)明顯相關(guān), 因此東北淺震的動(dòng)力來(lái)源也是板塊俯沖作用. 由于巖石層俯沖深度大于70 km時(shí), 應(yīng)力狀態(tài)是壓性的, 主壓應(yīng)力軸的方向與俯沖板塊的傾斜方向一致, 因而震源機(jī)制解中的P軸方向即代表俯沖板塊的運(yùn)動(dòng)方向. 呂政(2003)研究吉林省汪清深震時(shí)指出,2002年6月29日地震矩張量解顯示的震源類型是逆斷層性質(zhì),而1999年4月8日地震則是逆斷層兼走滑分量,但兩次地震的壓應(yīng)力軸均為WNW方向,反映了太平洋板塊向NW方向俯沖并與向東漂移的歐亞板塊碰撞的事實(shí).

    上述研究結(jié)果為地電阻率震前出現(xiàn)長(zhǎng)趨勢(shì)下降以及各向異性變化提供了物理解釋. 參照杜學(xué)彬等(2007)關(guān)于地震前地電阻率各向異性變化的研究,作者認(rèn)為,由于榆樹臺(tái)、 四平臺(tái)位于汪清、 琿春深震NW向或近NW向,在孕震擴(kuò)容階段,臺(tái)站所處介質(zhì)內(nèi)部微破裂數(shù)量急劇增大,裂隙走向沿最大主壓應(yīng)力方向(NW向)優(yōu)勢(shì)排列,導(dǎo)電流體快速進(jìn)入或重新分布,導(dǎo)致真電阻率快速、 大幅度下降變化,從而引起地電阻率快速下降變化.

    在強(qiáng)地震孕震晚期階段,由于太平洋板塊向NW方向擠壓,導(dǎo)致震源區(qū)及附近最大主壓應(yīng)力方位(或近于該方位)NW向的擠壓作用突出. 處于NW向的榆樹臺(tái)和四平臺(tái)介質(zhì)內(nèi)部豎向微裂隙發(fā)育、 快速發(fā)展,微裂隙走向沿最大加壓方向優(yōu)勢(shì)取向,導(dǎo)電水溶液快速進(jìn)入或重新分布,產(chǎn)生沿最大加壓方向快速、 大幅度變化為主的真電阻率各向異性變化. 由于真、 視電阻率變化最顯著的方向不同,所以產(chǎn)生了垂直加壓方向快速、 大幅度變化為主的視電阻率各向異性變化. 這一觀點(diǎn)與杜學(xué)彬等(2007)給出的強(qiáng)震前視電阻率各向異性變化的原因相吻合.

    5 結(jié)論

    綜上所述,通過(guò)對(duì)吉林省內(nèi)榆樹和四平兩個(gè)地電阻率臺(tái)站觀測(cè)數(shù)據(jù)的處理與分析研究表明,在中強(qiáng)地震及深源地震前,地電阻率能夠記錄到明顯的異常變化. 這一結(jié)果可以為吉林省地震前期預(yù)測(cè)提供一定的參考依據(jù).

    1) 榆樹臺(tái)和四平臺(tái)觀測(cè)到了2年尺度長(zhǎng)趨勢(shì)的地電阻率下降異?,F(xiàn)象. 在震前,四平臺(tái)地電阻率各向異性度出現(xiàn)下降,而榆樹臺(tái)各向異性度則出現(xiàn)破年變變化. 這兩個(gè)臺(tái)站的地電阻率歸一化變化速率出現(xiàn)短期至2年尺度異常變化.

    2) 榆樹臺(tái)、 四平臺(tái)所觀測(cè)到的地電阻率2 年尺度長(zhǎng)趨勢(shì)下降異常,源于介質(zhì)內(nèi)部微破裂、 微裂隙數(shù)目急劇增大,裂隙走向沿最大主壓應(yīng)力方向(NW向)優(yōu)勢(shì)排列,導(dǎo)電流體快速進(jìn)入或重新分布,導(dǎo)致真電阻率快速、 大幅度下降變化.

    3) 由于太平洋板塊向NW方向擠壓,導(dǎo)致在強(qiáng)地震孕震晚期階段,震源區(qū)及附近最大主壓應(yīng)力方位(或近于該方位)NW向的擠壓作用突出. 處于NW向的榆樹臺(tái)和四平臺(tái)介質(zhì)內(nèi)部微裂隙發(fā)育、 快速發(fā)展,微裂隙走向沿最大加壓方向優(yōu)勢(shì)取向,導(dǎo)電水溶液快速進(jìn)入或重新分布,產(chǎn)生沿最大加壓方向快速、 大幅度變化為主的真電阻率各向異性變化,最終導(dǎo)致地面觀測(cè)的視電阻率各向異性變化.

    地震的孕育和發(fā)生是一個(gè)極其復(fù)雜的過(guò)程,吉林省又屬于少震省份,可供借鑒的震例有限. 上述地電阻率異常變化是否具有普適性,仍有待震例資料的進(jìn)一步積累和檢驗(yàn).

    中國(guó)地震局地震預(yù)測(cè)研究所錢家棟研究員、 趙家騮研究員和中國(guó)地震局地震預(yù)測(cè)研究所蘭州科技創(chuàng)新基地杜學(xué)彬研究員以及上海市地震局馬欽忠研究員對(duì)本研究提出了有益建議; 審稿專家及編輯部老師對(duì)稿件提出了修改意見與建議; 本文資料由吉林省地震局四平地震臺(tái)和榆樹地震臺(tái)提供. 作者在此一并表示感謝!

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

    Du X B. 2010. Two types of changes in apparent resistivity in earthquake prediction[J].ScienceChina:EarthSciences, 40(10): 1321--1330 (in Chinese).

    杜學(xué)彬, 任廣軍, 薛順章. 1999. 中國(guó)大陸多種前兆異常研究及強(qiáng)震的嘗試性預(yù)測(cè)[J]. 西北地震學(xué)報(bào), 21(2): 113--122.

    Du X B, Ren G J, Xue S Z. 1999. Study on many kinds of precursory anomalies and trial prediction of strong earthquakes in the continent of China[J].NorthwesternSeismologicalJournal, 21(2): 113--122 (in Chinese).

    杜學(xué)彬, 薛順章, 郝臻, 張世中. 2000. 地電阻率中短期異常與地震的關(guān)系[J]. 地震學(xué)報(bào), 22(4): 368--376.

    采用累計(jì)預(yù)扣法,到2019年2月,小周應(yīng)繳稅稅款=(10000×2-5000×2-1500×2-1000×2)×3%-75(累計(jì)已預(yù)扣繳稅額)=75元。以此類推。

    Du X B, Xue S Z, Hao Z, Zhang S Z. 2000. On the relation of moderate-short term anomaly of earth resistivity to earthquake[J].ActaSeismologicaSinica, 22(4): 368--376 (in Chinese).

    杜學(xué)彬, 阮愛國(guó), 范世宏, 郝臻. 2001. 強(qiáng)震近震中區(qū)地電阻率變化速率的各向異性[J]. 地震學(xué)報(bào), 23(3): 289--297.

    Du X B, Ruan A G, Fan S H, Hao Z. 2001. Anisotropy of the variation rate of apparent resistivity near the epicentral region of strong earthquakes[J].ActaSeismologicaSinica, 23(3): 289--297 (in Chinese).

    杜學(xué)彬, 趙家騮, 譚大誠(chéng), 席繼樓, 錢家棟, 盧軍, 唐宇雄, 康好林, 陸陽(yáng)泉, 陳有發(fā), 王得志. 2006. 地震臺(tái)站建設(shè)規(guī)范: 地電臺(tái)站(第1部分): 地電阻率臺(tái)站(DB/T 18.1—2006)[S]. 北京: 地震出版社: 2--3, 9--12.

    Du X B, Zhao J L, Tan D C, Xi J L, Qian J D, Lu J, Tang Y X, Kang H L, Lu Y Q, Chen Y F, Wang D Z. 2006.SpecificationfortheConstructionofSeismicStation:GeoelectricalStation,Part1:GeoelctricalResistivityObservatory(DB/T18.1—2006)[S]. Beijing: Seismological Press: 2--3, 9--12 (in Chinese).

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

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

    馮志生, 楊建軍, 梅衛(wèi)萍, 耿杰, 汪雪泉, 劉義高. 2004. 華東地區(qū)地電阻率各向異性度的地震前兆異常特征初步研究[J]. 地震學(xué)報(bào), 26(2): 223--227.

    Feng Z S, Yang J J, Mei W P, Geng J, Wang X Q, Liu Y G. 2004. Preliminary study on characteristics of seismic precursor of earth resistivity anisotropy degree in East China[J].ActaSeismologicaSinica, 26(2): 223--227 (in Chinese).

    關(guān)華平, 張學(xué)民, 李美. 2008. 汶川地震前的地電阻率異常分析[J]. 國(guó)際地震動(dòng)態(tài), (11): 144.

    Guan H P, Zhang X M, Li M. 2008. Study on geo-resistivity anomaly before Wenchuan earthquake[J].RecentDevelopmentsinWorldSeismology, (11): 144 (in Chinese).

    吉林省地震局. 2005. 吉林省地震監(jiān)測(cè)志[M]. 北京: 地震出版社: 75--85, 128--135, 138--148.

    Earthquake Administration of Jilin Province. 2005.SeismicMonitoringAnnalsofJilinProvince[M]. Beijing: Seismological Press: 75--85, 128--135, 138--148 (in Chinese).

    陸陽(yáng)泉, 劉建毅, 梁子斌. 1998. 巖石電學(xué)性質(zhì)的實(shí)驗(yàn)研究概述[G]∥地震地電學(xué)發(fā)展與展望. 蘭州: 蘭州大學(xué)出版社: 32--34.

    Lu Y Q, Liu J Y, Liang Z B. 1998. Overview of experimental studies on the electrical property of rocks[G]∥DevelopmentandProspectofSeismicGeoelectricity. Lanzhou: Lanzhou University Press: 32--34 (in Chinese).

    呂政. 2003. 吉林省汪清兩次深源地震的研究[J]. 東北地震研究, 19(4): 1--8.

    Lü Z. 2003. Research on the two deep earthquakes in Jilin Province[J].SeismologicalResearchofNortheastChina, 19(4): 1--8 (in Chinese).

    毛桐恩, 王鐵城. 1995. 唐山地震地電阻率各向異性變化特征研究[J]. 地震學(xué)報(bào), 17(4): 505--510.

    Mao T E, Wang T C. 1995. Variations of ground resistivity anisotropy during the Tangshan earthquake[J].ActaSeismologicaSinica, 17(4): 505--510 (in Chinese).

    毛桐恩, 陳大元, 劉新美, 余素榮, 宋大捷. 1998. 承載巖石地電阻率各向異性變化的實(shí)驗(yàn)研究與現(xiàn)場(chǎng)觀測(cè)結(jié)果[J]. 地震, 18(S1): 33--38.

    Mao T E, Chen D Y, Liu X M, Yu S R, Song D J. 1998. Experimental research on change of georesistivity anisotropy of pressed rocks and on the spot observation[J].Earthquake, 18(S1): 33--38 (in Chinese).

    毛桐恩, 胥廣銀, 范思源, 趙明, 孫景芳. 1999. 地電阻率各向異性度的動(dòng)態(tài)演化圖像與地震孕育過(guò)程[J]. 地震學(xué)報(bào), 21(2): 180--186.

    Mao T E, Xu G Y, Fan S Y, Zhao M, Sun J F. 1999. Dynamic evolution patterns of ground resistivity anisotropy and the seismogenic process[J].ActaSeismologicaSinica, 21(2): 180--186 (in Chinese).

    孟憲森, 朱景春, 孫文斌, 許學(xué)禮, 李登恒. 1996. 東北地區(qū)淺源中強(qiáng)震及深震與西太平洋板塊俯沖[J]. 東北地震研究, 12(2): 12--23.

    Meng X S, Zhu J C, Sun W B, Xu X L, Li D H. 1996. Shallow moderate and deep earthquake in Northeast China and Western Pacific Plate subduction[J].SeismologicalResearchofNortheastChina, 12(2): 12--23 (in Chinese).

    錢復(fù)業(yè), 趙玉林. 1980. 地震前地電阻率變化十例[J]. 地震學(xué)報(bào), 2(2): 186--197.

    Qian F Y, Zhao Y L. 1980. Ten examples of changes in earth-resistivity prior to strong earthquakes[J].ActaSeismologicaSinica, 2(2): 186--197 (in Chinese).

    錢復(fù)業(yè), 趙玉林. 1998. 地電阻率的各向異性地震前兆及其在探索震中區(qū)域應(yīng)力場(chǎng)分布中的應(yīng)用[G]∥地震地電學(xué)發(fā)展與展望. 蘭州: 蘭州大學(xué)出版社: 64--69.

    Qian F Y, Zhao Y L. 1998. Precursory of anisotropy of earth resistivity and its application to exploring distribution of regional stress field in the epicenter region[G]∥DevelopmentandProspectofSeismicGeoelectricity. Lanzhou: Lanzhou University Press: 64--69 (in Chinese).

    錢家棟. 1993. 與大震孕育過(guò)程有關(guān)的地電阻率變化研究[J]. 中國(guó)地震, 9(4): 341--350.

    Qian J D. 1993. A study on the changes in geoelectrical resistivity associated with preparatory process of great earthquakes in China[J].EarthquakeResearchinChina, 9(4): 341--350 (in Chinese).

    阮愛國(guó), 盧軍, 符瓊玉, 吳振利. 2004. 兩次大震的視電阻率各向異性特征[J]. 華南地震, 24(4): 1--7.

    Ruan A G, Lu J, Fu Q Y, Wu Z L. 2004. Anisotropic characteristics of earth apparent resistivity during the two strong earthquakes[J].SouthChinaJournalofSeismology, 24(4): 1--7 (in Chinese).

    史紅軍, 趙衛(wèi)星, 龐晶源. 2008. 吉林汪清深源地震前榆樹地震臺(tái)地電阻率異常分析[J]. 地震研究, 31(3): 233--237.

    Shi H J, Zhao W X, Pang J Y. 2008. Earth resistivity anomalies observed at Yushu Seismic Station before Wangqing deep earthquakes[J].JournalofSeismologicalResearch, 31(3): 233--237 (in Chinese).

    史紅軍, 趙衛(wèi)星, 張振文, 楊柳. 2011. 榆樹地電阻率臺(tái)觀測(cè)效能評(píng)價(jià)[J]. 防災(zāi)減災(zāi)學(xué)報(bào), 27(1): 57--64.

    Shi H J, Zhao W X, Zhang Z W, Yang L. 2011. Assessment on the observation efficiency of Yushu earth resistivity station[J].JournalofDisasterPreventionandReduction, 27(1): 57--64 (in Chinese).

    汪志亮, 王志賢, 余素榮. 1995. 地電前兆異常特征及瀾滄—耿馬7.6級(jí)等地震預(yù)報(bào)問(wèn)題[J]. 地震, 25(3): 240--245.

    Wang Z L, Wang Z X, Yu S R. 1995. Features of anomalous precursors and prediction of Lancang--Gengma earthquake withM=7.6 and other events[J].Earthquake, 25(3): 240--245 (in Chinese).

    張繼紅, 趙國(guó)澤, 王晉生, 唐廷梅, 齊麗杰. 2010. 2008年汶川MS8.0地震與青島臺(tái)地電阻率異常變化關(guān)系[J]. 地震地質(zhì), 32(3): 409--416.

    Zhang J H, Zhao G Z, Wang J S, Tang T M, Qi L J. 2010. Preliminary study on the relationship between theMS8.0 Wenchuan earthquake of 2008 and the anomalies of earth resistivity observed at Qingdao seismic station[J].SeismologyandGeology, 32(3): 409--416 (in Chinese).

    張立敏, 唐曉明. 1983. 西太平洋板塊俯沖運(yùn)動(dòng)與中國(guó)東北深震帶[J]. 地球物理學(xué)報(bào), 26(4): 331--340.

    Zhang L M, Tang X M. 1983. The underthrusting movement of the Western Pacific Plate and the deep focus earthquake zone of Northeast China[J].ChineseJournalofGeophysics, 26(4): 331--340 (in Chinese).

    張學(xué)民, 李美, 關(guān)華平. 2009. 汶川8.0級(jí)地震前的地電阻率異常分析[J]. 地震, 29(1): 108--115.

    Zhang X M, Li M, Guan H P. 2009. Anomaly analysis of earth resistivity observations before the Wenchuan earthquake[J].Earthquake, 29(1): 108--115 (in Chinese).

    趙玉林, 錢復(fù)業(yè). 1978. 唐山7.8級(jí)強(qiáng)震前震中周圍形變電阻率的下降異常[J]. 地球物理學(xué)報(bào), 21(3): 181--190.

    Zhao Y L, Qian F Y. 1978. Electrical resistivity anomaly observed in and around the epicentral area prior to the Tang-shan earthquake of 1976[J].ActaGeophysicaSinica, 21(3): 181--190 (in Chinese).

    中國(guó)地震局. 2001. 地震及前兆數(shù)字觀測(cè)技術(shù)規(guī)范: 電磁觀測(cè)(試行)[M]. 北京: 地震出版社: 3--6, 19--36.

    China Earthquake Administration. 2001.TechnicalSpecificationforDigitalObservationofEarthquakesandPrecursors:ElectromagneticObservation(ForTrialImplementation)[M]. Beijing: Seismological Press: 3--6, 19--36 (in Chinese).

    Qian J D, Gui X T, Lu Y Q. 1996. Review on the observation and experiment in earthquake prediction by using geoelectrical methods[G]∥TheSelectedPapersofEarthquakePredictioninChina. Beijing: Seismological Press: 170.

    Characteristics of precursor anomalies of earth resistivity at Yushu and Siping seismic stations of Jilin Province

    1)YushuSeismicStation,Changchun130400,China2)EarthquakeFastDeterminationandReportCenterofChangchunCity,Changchun130012,China

    This paper studies the anomalous changes of earth resistivity before 1999 WangqingMS7.0 earthquake and 2006 Qian’an--QianguoMS5.0 earthquake observed at Yushu and Siping stations of Jilin Province. The results show that, at the Yushu and Siping stations a long-trend decreasing or short-impending anomaly of earth resistivity appeared before the two earthquakes. Monthly mean curves of earth resistivity displayed two-year decrease anomalies at Yushu and Siping stations before the Wangqing earthquake. Before Wangqing and Qian’an--Qianguo earthquakes, normalized variation rates of earth resistivity are larger than anomaly index at Yushu and Siping stations. Before Qian’an--Qianguo earthquake, earth resistivity anisotropy degree showed obviously low value anomaly at Siping station, and at Yushu station it broke its annual variation value. Possible reasons for the anomalous changes maybe are related to the maximum principal compressive stress orientation’s (or close to the range) compressing NWwards in the source region and its vicinity resulting from NWwards extrusion of Pacific Plate, which makes the conductive fluid in the media quickly enter or redistribute, so that the earth resistivity shows a long-trend decrease and anisotropic variation.

    earth resistivity; earthquake; earthquake precursor anomaly; shape anomaly; normalized variation rate; anisotropy degree

    10.3969/j.issn.0253-3782.2014.03.011.

    中國(guó)地震局“地震監(jiān)測(cè)、 預(yù)報(bào)、 科研三結(jié)合”課題(201112, 140701)資助.

    2012-12-03收到初稿,2013-03-10決定采用修改稿.

    e-mail: dzshj@sina.com

    10.3969/j.issn.0253-3782.2014.03.011

    P319.3+2

    A

    史紅軍, 趙衛(wèi)星, 張可佳, 朱偉楠, 李寧. 2014. 吉林省地電阻率震兆異常分析研究. 地震學(xué)報(bào), 36(3): 452--463.

    Shi H J, Zhao W X, Zhang K J, Zhu W N, Li N. 2014. Characteristics of precursor anomalies of earth resistivity at Yushu and Siping seismic stations of Jilin Province.ActaSeismologicaSinica, 36(3): 452--463. doi:10.3969/j.issn.0253-3782.2014.03.011.

    日韩中文字幕欧美一区二区 | 亚洲伊人色综图| 国产高清videossex| 久久久久久久久久久久大奶| 精品免费久久久久久久清纯 | 国产色视频综合| 十八禁人妻一区二区| 操出白浆在线播放| 操出白浆在线播放| 亚洲伊人久久精品综合| e午夜精品久久久久久久| 极品人妻少妇av视频| 国产精品 欧美亚洲| 久久影院123| 一本久久精品| 建设人人有责人人尽责人人享有的| 新久久久久国产一级毛片| 国产不卡av网站在线观看| 夫妻性生交免费视频一级片| 美女脱内裤让男人舔精品视频| 天天躁夜夜躁狠狠久久av| 午夜老司机福利片| 国产片内射在线| 一级毛片我不卡| 麻豆乱淫一区二区| 在线观看人妻少妇| 久久天堂一区二区三区四区| 欧美日韩综合久久久久久| 免费看十八禁软件| 国产成人免费观看mmmm| 亚洲av电影在线观看一区二区三区| 国产成人精品久久二区二区免费| 欧美成人午夜精品| 老汉色∧v一级毛片| 人人妻,人人澡人人爽秒播 | 国产成人精品在线电影| 久久人妻熟女aⅴ| 一级,二级,三级黄色视频| 日本av免费视频播放| 观看av在线不卡| 操美女的视频在线观看| 欧美中文综合在线视频| 国产视频首页在线观看| 一级毛片女人18水好多 | av有码第一页| 亚洲欧美一区二区三区黑人| 中国美女看黄片| 日本av免费视频播放| 一区二区av电影网| 亚洲综合色网址| 搡老乐熟女国产| 国语对白做爰xxxⅹ性视频网站| 国产在视频线精品| 男女午夜视频在线观看| 亚洲伊人色综图| 国产成人系列免费观看| 欧美在线黄色| 男女无遮挡免费网站观看| 99re6热这里在线精品视频| 最近中文字幕2019免费版| 成人免费观看视频高清| 男女之事视频高清在线观看 | 亚洲免费av在线视频| 国产熟女午夜一区二区三区| 精品国产乱码久久久久久小说| 性少妇av在线| 首页视频小说图片口味搜索 | 国产一区有黄有色的免费视频| 天天躁夜夜躁狠狠躁躁| av不卡在线播放| 久久午夜综合久久蜜桃| 亚洲欧美成人综合另类久久久| 男女免费视频国产| 日本91视频免费播放| 人人妻人人添人人爽欧美一区卜| 国产成人影院久久av| 欧美老熟妇乱子伦牲交| 免费看av在线观看网站| 飞空精品影院首页| 老司机亚洲免费影院| 亚洲情色 制服丝袜| 中文字幕色久视频| 亚洲一卡2卡3卡4卡5卡精品中文| 满18在线观看网站| 国产成人精品久久二区二区91| 大型av网站在线播放| 男的添女的下面高潮视频| 永久免费av网站大全| 男女边摸边吃奶| 久久久亚洲精品成人影院| 国产精品av久久久久免费| 亚洲男人天堂网一区| 久久狼人影院| 亚洲图色成人| 性色av一级| 97人妻天天添夜夜摸| 丰满少妇做爰视频| 久久青草综合色| 男女国产视频网站| 夫妻午夜视频| 国产成人啪精品午夜网站| 国产精品亚洲av一区麻豆| 亚洲av成人精品一二三区| 美女主播在线视频| 国产精品一国产av| 最新的欧美精品一区二区| 男女下面插进去视频免费观看| 叶爱在线成人免费视频播放| 国产免费现黄频在线看| 亚洲视频免费观看视频| av线在线观看网站| 老熟女久久久| 国产又色又爽无遮挡免| 亚洲成色77777| 亚洲av国产av综合av卡| 国产精品av久久久久免费| 五月开心婷婷网| 国产精品国产av在线观看| 极品少妇高潮喷水抽搐| 国产无遮挡羞羞视频在线观看| 2021少妇久久久久久久久久久| 校园人妻丝袜中文字幕| 2018国产大陆天天弄谢| 久久久精品国产亚洲av高清涩受| 男的添女的下面高潮视频| 最近中文字幕2019免费版| 国产一区二区三区综合在线观看| h视频一区二区三区| av线在线观看网站| 欧美日韩黄片免| 成年人午夜在线观看视频| 国精品久久久久久国模美| 亚洲精品国产色婷婷电影| 国产精品亚洲av一区麻豆| 在线观看免费午夜福利视频| 免费看十八禁软件| 国产福利在线免费观看视频| 亚洲精品美女久久av网站| 国产91精品成人一区二区三区 | 黄色视频不卡| 精品少妇黑人巨大在线播放| 免费看不卡的av| 在线观看免费视频网站a站| 少妇精品久久久久久久| 精品人妻一区二区三区麻豆| 欧美精品高潮呻吟av久久| 欧美在线黄色| 久久毛片免费看一区二区三区| h视频一区二区三区| 99精国产麻豆久久婷婷| 香蕉国产在线看| 多毛熟女@视频| 我的亚洲天堂| 一本久久精品| 每晚都被弄得嗷嗷叫到高潮| 欧美成人午夜精品| 在线观看人妻少妇| 国产精品.久久久| 免费高清在线观看视频在线观看| 久久久久久亚洲精品国产蜜桃av| 50天的宝宝边吃奶边哭怎么回事| 精品国产国语对白av| 国产成人精品久久久久久| 视频在线观看一区二区三区| 国产成人av教育| 男女国产视频网站| 精品国产国语对白av| 亚洲午夜精品一区,二区,三区| 妹子高潮喷水视频| 亚洲国产欧美日韩在线播放| 欧美精品av麻豆av| 亚洲精品久久午夜乱码| 丝袜美足系列| 在现免费观看毛片| 热99国产精品久久久久久7| 亚洲国产精品一区三区| 免费人妻精品一区二区三区视频| 99香蕉大伊视频| 大话2 男鬼变身卡| 交换朋友夫妻互换小说| 只有这里有精品99| 国产成人一区二区三区免费视频网站 | svipshipincom国产片| 首页视频小说图片口味搜索 | 黄片小视频在线播放| 色精品久久人妻99蜜桃| 色婷婷av一区二区三区视频| 50天的宝宝边吃奶边哭怎么回事| 亚洲av电影在线观看一区二区三区| 欧美日韩亚洲综合一区二区三区_| 丝袜喷水一区| 国产一区二区三区综合在线观看| 天天操日日干夜夜撸| 亚洲一卡2卡3卡4卡5卡精品中文| 日韩一卡2卡3卡4卡2021年| 人人澡人人妻人| 亚洲av日韩在线播放| 国产欧美日韩一区二区三 | 国产日韩一区二区三区精品不卡| 午夜福利,免费看| 亚洲精品美女久久久久99蜜臀 | 免费高清在线观看视频在线观看| a级毛片在线看网站| 美女中出高潮动态图| 欧美人与善性xxx| 香蕉丝袜av| 亚洲人成电影观看| 性高湖久久久久久久久免费观看| 国产精品国产av在线观看| 美女脱内裤让男人舔精品视频| 精品熟女少妇八av免费久了| 免费av中文字幕在线| 九色亚洲精品在线播放| 在线观看人妻少妇| 色网站视频免费| 免费在线观看黄色视频的| 多毛熟女@视频| 欧美精品啪啪一区二区三区 | av天堂久久9| 久久久精品区二区三区| 免费看不卡的av| 日韩一本色道免费dvd| 成人18禁高潮啪啪吃奶动态图| 另类精品久久| 美女高潮到喷水免费观看| 国产成人91sexporn| 女性被躁到高潮视频| 欧美黑人精品巨大| 视频在线观看一区二区三区| 人人妻人人添人人爽欧美一区卜| 亚洲欧美激情在线| 19禁男女啪啪无遮挡网站| 性少妇av在线| 一边摸一边抽搐一进一出视频| 欧美变态另类bdsm刘玥| 日韩人妻精品一区2区三区| 又粗又硬又长又爽又黄的视频| 无限看片的www在线观看| 午夜视频精品福利| 99re6热这里在线精品视频| 嫁个100分男人电影在线观看 | 久久亚洲精品不卡| 欧美在线一区亚洲| 亚洲精品国产一区二区精华液| 2021少妇久久久久久久久久久| av线在线观看网站| 老司机影院毛片| 色婷婷av一区二区三区视频| 人妻人人澡人人爽人人| 精品亚洲乱码少妇综合久久| 十八禁网站网址无遮挡| 又大又黄又爽视频免费| 美国免费a级毛片| 天天影视国产精品| 你懂的网址亚洲精品在线观看| 91麻豆av在线| 狂野欧美激情性bbbbbb| 日韩欧美一区视频在线观看| 电影成人av| 精品一区在线观看国产| 激情视频va一区二区三区| 亚洲国产日韩一区二区| 青草久久国产| 国产精品一区二区免费欧美 | 精品少妇一区二区三区视频日本电影| 亚洲成人国产一区在线观看 | 亚洲中文字幕日韩| 久久天堂一区二区三区四区| 电影成人av| 国产日韩一区二区三区精品不卡| 欧美成人午夜精品| 国产又色又爽无遮挡免| 国产女主播在线喷水免费视频网站| 操美女的视频在线观看| 脱女人内裤的视频| 啦啦啦啦在线视频资源| 高清视频免费观看一区二区| 只有这里有精品99| 国产精品亚洲av一区麻豆| 色播在线永久视频| 欧美日韩精品网址| 免费在线观看影片大全网站 | av在线老鸭窝| tube8黄色片| 亚洲第一青青草原| 久久中文字幕一级| 国产一区二区三区av在线| 久久久国产一区二区| 亚洲国产精品一区三区| 亚洲国产精品成人久久小说| 亚洲精品国产区一区二| h视频一区二区三区| 成在线人永久免费视频| 男人操女人黄网站| av网站在线播放免费| netflix在线观看网站| 建设人人有责人人尽责人人享有的| 男人添女人高潮全过程视频| 国产精品.久久久| 国产成人av教育| 看免费av毛片| 婷婷成人精品国产| 成人亚洲欧美一区二区av| 熟女少妇亚洲综合色aaa.| 亚洲国产欧美日韩在线播放| 男女边摸边吃奶| 亚洲色图综合在线观看| 亚洲人成网站在线观看播放| 亚洲精品乱久久久久久| 18禁观看日本| 亚洲精品日韩在线中文字幕| 免费在线观看视频国产中文字幕亚洲 | 人人妻人人爽人人添夜夜欢视频| 大香蕉久久成人网| 亚洲天堂av无毛| 亚洲熟女精品中文字幕| av网站在线播放免费| 咕卡用的链子| 欧美+亚洲+日韩+国产| 亚洲中文av在线| 这个男人来自地球电影免费观看| 久久国产精品大桥未久av| 18禁黄网站禁片午夜丰满| 多毛熟女@视频| 日韩中文字幕欧美一区二区 | 午夜福利在线免费观看网站| 午夜av观看不卡| svipshipincom国产片| 国产在视频线精品| 麻豆乱淫一区二区| 无遮挡黄片免费观看| 亚洲成人国产一区在线观看 | 欧美国产精品一级二级三级| av不卡在线播放| 国产在线观看jvid| 国产精品一区二区在线观看99| 欧美日韩成人在线一区二区| 一区二区av电影网| av片东京热男人的天堂| 91麻豆av在线| 欧美大码av| 久久久久久久久免费视频了| av国产久精品久网站免费入址| 欧美精品亚洲一区二区| 精品国产超薄肉色丝袜足j| 久久久久国产精品人妻一区二区| 欧美日韩黄片免| 国产爽快片一区二区三区| 91成人精品电影| 精品人妻一区二区三区麻豆| 成年动漫av网址| 亚洲第一青青草原| 久久人妻福利社区极品人妻图片 | 日韩一卡2卡3卡4卡2021年| 91麻豆av在线| 国产日韩一区二区三区精品不卡| 热re99久久国产66热| 免费日韩欧美在线观看| 下体分泌物呈黄色| 久热爱精品视频在线9| 国产精品一二三区在线看| 国产精品熟女久久久久浪| 欧美日韩视频精品一区| 高潮久久久久久久久久久不卡| 精品亚洲成a人片在线观看| 亚洲欧美一区二区三区黑人| 一级片'在线观看视频| 亚洲国产精品成人久久小说| 1024香蕉在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 高潮久久久久久久久久久不卡| 久久精品久久久久久久性| 黄色一级大片看看| 极品人妻少妇av视频| 又黄又粗又硬又大视频| 成年动漫av网址| 亚洲精品日韩在线中文字幕| 久久人人爽av亚洲精品天堂| 999久久久国产精品视频| 色网站视频免费| 午夜激情久久久久久久| 亚洲五月色婷婷综合| 免费在线观看完整版高清| 亚洲欧美成人综合另类久久久| 免费一级毛片在线播放高清视频 | 欧美中文综合在线视频| 久久久久久免费高清国产稀缺| 亚洲七黄色美女视频| 人人妻人人添人人爽欧美一区卜| 欧美精品高潮呻吟av久久| 波多野结衣一区麻豆| 我的亚洲天堂| av福利片在线| 王馨瑶露胸无遮挡在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 只有这里有精品99| 超碰成人久久| 午夜影院在线不卡| 亚洲欧美日韩另类电影网站| 欧美变态另类bdsm刘玥| 丝袜人妻中文字幕| 亚洲 欧美一区二区三区| 午夜免费成人在线视频| 国产精品一区二区精品视频观看| 热99国产精品久久久久久7| 亚洲成人手机| 乱人伦中国视频| 性色av乱码一区二区三区2| 亚洲av电影在线观看一区二区三区| 久久久久久久久久久久大奶| 黑人猛操日本美女一级片| 久热爱精品视频在线9| 国产精品.久久久| 国产男人的电影天堂91| 亚洲三区欧美一区| 欧美老熟妇乱子伦牲交| 亚洲伊人色综图| 亚洲av日韩在线播放| 50天的宝宝边吃奶边哭怎么回事| 老司机靠b影院| 久久人妻熟女aⅴ| 国产成人免费无遮挡视频| 五月天丁香电影| 久热爱精品视频在线9| 国产精品二区激情视频| 最新在线观看一区二区三区 | 午夜免费成人在线视频| 国产熟女欧美一区二区| 人妻 亚洲 视频| 亚洲七黄色美女视频| 9色porny在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 日韩大码丰满熟妇| 日韩免费高清中文字幕av| 1024香蕉在线观看| 国产成人欧美在线观看 | 亚洲国产欧美网| 赤兔流量卡办理| 各种免费的搞黄视频| 国产高清视频在线播放一区 | h视频一区二区三区| 亚洲欧美精品自产自拍| 精品久久蜜臀av无| 不卡av一区二区三区| 亚洲精品一区蜜桃| 国产在线免费精品| 两性夫妻黄色片| 亚洲成国产人片在线观看| 国产欧美亚洲国产| 亚洲精品自拍成人| 各种免费的搞黄视频| 99国产精品99久久久久| 欧美亚洲日本最大视频资源| 国产一卡二卡三卡精品| 久久天躁狠狠躁夜夜2o2o | 亚洲成人手机| 欧美另类一区| h视频一区二区三区| 搡老乐熟女国产| 国产一区二区三区综合在线观看| 午夜91福利影院| 久久精品国产亚洲av高清一级| 国产黄色视频一区二区在线观看| 男人操女人黄网站| 男女午夜视频在线观看| 国产亚洲午夜精品一区二区久久| 欧美黄色片欧美黄色片| 大片免费播放器 马上看| 国产伦人伦偷精品视频| 精品少妇一区二区三区视频日本电影| 亚洲天堂av无毛| 国产在线一区二区三区精| 日韩中文字幕视频在线看片| 亚洲人成网站在线观看播放| 亚洲 欧美一区二区三区| 国产av精品麻豆| 女性被躁到高潮视频| 汤姆久久久久久久影院中文字幕| 亚洲专区中文字幕在线| 91麻豆av在线| 亚洲国产精品999| 美女福利国产在线| 大话2 男鬼变身卡| 久久久精品区二区三区| 精品一区二区三区四区五区乱码 | 免费在线观看黄色视频的| videos熟女内射| 亚洲国产av影院在线观看| 欧美精品亚洲一区二区| 韩国高清视频一区二区三区| 国产男女超爽视频在线观看| 色视频在线一区二区三区| 久热这里只有精品99| 少妇人妻久久综合中文| 成年av动漫网址| 99精品久久久久人妻精品| 亚洲精品在线美女| 国产成人影院久久av| 亚洲精品日本国产第一区| 性色av乱码一区二区三区2| 久久国产亚洲av麻豆专区| 一本一本久久a久久精品综合妖精| 免费av中文字幕在线| 久久精品熟女亚洲av麻豆精品| 中文精品一卡2卡3卡4更新| 日本五十路高清| 久久这里只有精品19| 亚洲 国产 在线| 日本午夜av视频| 国产亚洲一区二区精品| 满18在线观看网站| av国产久精品久网站免费入址| 国产一区二区激情短视频 | 亚洲激情五月婷婷啪啪| 久久国产亚洲av麻豆专区| 欧美人与善性xxx| 多毛熟女@视频| 一个人免费看片子| 老司机午夜十八禁免费视频| 国产欧美日韩一区二区三区在线| 国产精品偷伦视频观看了| 秋霞在线观看毛片| 亚洲激情五月婷婷啪啪| 一本久久精品| 黑丝袜美女国产一区| 国产av国产精品国产| 黄片小视频在线播放| 国产精品免费大片| 国产在线观看jvid| 欧美黄色淫秽网站| 精品人妻1区二区| 菩萨蛮人人尽说江南好唐韦庄| 女人被躁到高潮嗷嗷叫费观| 亚洲一区二区三区欧美精品| 亚洲专区国产一区二区| 大话2 男鬼变身卡| 超碰成人久久| 日韩人妻精品一区2区三区| 在线亚洲精品国产二区图片欧美| 又粗又硬又长又爽又黄的视频| 99re6热这里在线精品视频| 中文字幕精品免费在线观看视频| 国产精品偷伦视频观看了| 国产野战对白在线观看| 久久影院123| 精品国产一区二区三区久久久樱花| 一二三四在线观看免费中文在| 18禁黄网站禁片午夜丰满| 成年动漫av网址| 久久 成人 亚洲| 人成视频在线观看免费观看| 久久久久久久国产电影| 国产亚洲欧美精品永久| 国产一区二区激情短视频 | 精品一区二区三卡| 一区二区三区四区激情视频| 国产福利在线免费观看视频| 久久久久久久大尺度免费视频| 又粗又硬又长又爽又黄的视频| 国产一区二区三区综合在线观看| 男人爽女人下面视频在线观看| 中文字幕人妻熟女乱码| 一区福利在线观看| 久久 成人 亚洲| 日韩一本色道免费dvd| 美女视频免费永久观看网站| 又大又黄又爽视频免费| 狠狠精品人妻久久久久久综合| 欧美日韩一级在线毛片| 另类精品久久| 一区在线观看完整版| 亚洲成人免费电影在线观看 | 亚洲欧美日韩另类电影网站| 免费人妻精品一区二区三区视频| 丝袜在线中文字幕| 国产精品秋霞免费鲁丝片| 精品亚洲成国产av| 最新在线观看一区二区三区 | 国产精品久久久久久人妻精品电影 | 中文乱码字字幕精品一区二区三区| 欧美日韩亚洲综合一区二区三区_| 国产日韩一区二区三区精品不卡| 亚洲,欧美精品.| 午夜激情久久久久久久| 国产精品香港三级国产av潘金莲 | 国产成人系列免费观看| 熟女少妇亚洲综合色aaa.| 少妇被粗大的猛进出69影院| 丰满迷人的少妇在线观看| 国产一区二区三区综合在线观看| 亚洲国产欧美网| 久久精品国产亚洲av高清一级| 水蜜桃什么品种好| 欧美亚洲 丝袜 人妻 在线| 亚洲国产欧美日韩在线播放| 国产精品秋霞免费鲁丝片| 精品亚洲成a人片在线观看| av网站在线播放免费| 免费看不卡的av| 看免费av毛片| 国产精品av久久久久免费| 国产老妇伦熟女老妇高清| 中国美女看黄片| 悠悠久久av| 国产成人免费观看mmmm| 久久精品熟女亚洲av麻豆精品| 制服人妻中文乱码| 十八禁网站网址无遮挡| 丝袜脚勾引网站| 亚洲伊人久久精品综合| 欧美人与善性xxx| 国产精品久久久久久精品电影小说| 久久这里只有精品19| 亚洲男人天堂网一区| 男人爽女人下面视频在线观看|