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

    基于GPS資料研究云南地區(qū)地殼形變動(dòng)態(tài)特征

    2018-11-26 09:32王巖洪敏邵德勝汪志民
    地震研究 2018年3期

    王巖 洪敏 邵德勝 汪志民

    摘要:基于1999—2007年和2011—2013年2期云南地區(qū)的GPS速度場(chǎng)資料,結(jié)合最小二乘配置方法給出了2個(gè)時(shí)段應(yīng)變參數(shù)和GPS速度剖面結(jié)果,綜合分析了汶川8.0級(jí)地震前后云南地區(qū)地殼變形動(dòng)態(tài)演化特征。GPS 速度場(chǎng)動(dòng)態(tài)演化結(jié)果表明:整體上GPS速度場(chǎng)方向從川滇塊體南部的南向逐漸轉(zhuǎn)為滇南塊體的南東向,塊體邊界斷裂帶附近存在明顯的相對(duì)運(yùn)動(dòng)特征。相比于1999—2007年,點(diǎn)位更加密集的2011—2013年GPS速度場(chǎng)結(jié)果還反映了斷裂帶分段變形特征;在滇南塊體及周邊地區(qū),南東向的運(yùn)動(dòng)有顯著的增加。應(yīng)變參數(shù)動(dòng)態(tài)演化結(jié)果表明:整體上呈現(xiàn)“中部拉張、兩端擠壓”的特點(diǎn);區(qū)域變形(剪切、拉張或擠壓)的分布特征與區(qū)域主要斷裂帶的背景運(yùn)動(dòng)及變形特征密切相關(guān);相比于1999—2007年,2011—2013年的拉張變形區(qū)有所擴(kuò)大,并且有向東和向南擴(kuò)展的趨勢(shì);川滇塊體東邊界(小江斷裂帶)整體的擠壓變形已經(jīng)不再顯著,主要集中在小江斷裂帶北段靠北的區(qū)域??鐢嗔褞PS速度剖面結(jié)果顯示:小江斷裂帶左旋滑動(dòng)速率從北向南逐漸減?。ㄓ?0 mm/a減至5 mm/a),北段變形寬度較寬;紅河斷裂帶右旋滑動(dòng)速率約為4 mm/a,變形寬度較寬;動(dòng)態(tài)結(jié)果表明,汶川地震對(duì)2條斷裂帶變形模式影響較小。

    關(guān)鍵詞:最小二乘配置;GPS速度場(chǎng);應(yīng)變參數(shù);GPS速度剖面

    中圖分類(lèi)號(hào):P315.725 文獻(xiàn)標(biāo)識(shí)碼:A 文章編號(hào):1000-0666(2018)03-0368-07

    0 引言

    云南地區(qū)地處青藏高原東南緣地區(qū),是我國(guó)地震活動(dòng)最為顯著的地區(qū)之一,該區(qū)域活動(dòng)構(gòu)造豐富且復(fù)雜(張俊昌等,1980;王義昭等,1988;李西,2015),這種特殊的構(gòu)造環(huán)境也造就了該區(qū)域特殊的地殼變形特征。本區(qū)域位于南北地震南段,多位研究學(xué)者針對(duì)南北地震帶地區(qū)地殼變形的特征開(kāi)展了大量研究(江在森等,2001;方穎等,2008;楊國(guó)華等,2009;丁開(kāi)華等,2013;張希等,2003;武艷強(qiáng)等,2012;魏文薪等,2012;季靈運(yùn)等,2015),另外,也有學(xué)者針對(duì)云南地區(qū)的地殼變形特征開(kāi)展了相關(guān)的研究(楊國(guó)華等,2003;洪敏等,2014;王伶俐等,2015,周海濤等,2017),但針對(duì)云南地區(qū)的地殼變形的動(dòng)態(tài)演化特征研究還略顯不足。2008年汶川8.0級(jí)地震發(fā)生對(duì)云南地區(qū)地殼變形的影響如何,該地震前后區(qū)域的變形是否有顯著變化等問(wèn)題值得進(jìn)一步研究。因此,本文基于1999—2007年和2011—2013年云南地區(qū)的GPS速度場(chǎng)資料,結(jié)合最小二乘配置方法給出2個(gè)時(shí)段應(yīng)變參數(shù)結(jié)果和GPS速度剖面結(jié)果,綜合分析了汶川地震前后云南地區(qū)地殼變形動(dòng)態(tài)演化特征。

    1 GPS資料概況及數(shù)據(jù)處理

    我國(guó)從20世紀(jì)80年代末逐步開(kāi)始GPS觀測(cè),較大規(guī)模、高精度的GPS觀測(cè)從 “中國(guó)地殼運(yùn)動(dòng)觀測(cè)網(wǎng)絡(luò)”建設(shè)完成后開(kāi)始,該網(wǎng)絡(luò)分布于中國(guó)大陸的GPS區(qū)域站約1 000個(gè)。在中國(guó)地殼運(yùn)動(dòng)觀測(cè)網(wǎng)絡(luò)基礎(chǔ)上,2007年進(jìn)行二期建設(shè)的“中國(guó)大陸構(gòu)造環(huán)境監(jiān)測(cè)網(wǎng)絡(luò)”(簡(jiǎn)稱(chēng)陸態(tài)網(wǎng)絡(luò))把區(qū)域網(wǎng)點(diǎn)增加到2 000個(gè)。為研究云南地區(qū)地殼變形的動(dòng)態(tài)演化特征,本文主要利用中國(guó)地殼運(yùn)動(dòng)觀測(cè)網(wǎng)絡(luò)建成后的GPS復(fù)測(cè)資料進(jìn)行研究,這些資料包括:中國(guó)地殼運(yùn)動(dòng)觀測(cè)網(wǎng)絡(luò)GPS區(qū)域網(wǎng)1999,2001,2004,2007年復(fù)測(cè)資料,中國(guó)大陸構(gòu)造環(huán)境監(jiān)測(cè)網(wǎng)絡(luò)GPS區(qū)域網(wǎng)2011,2013年復(fù)測(cè)資料。

    GPS區(qū)域網(wǎng)速度場(chǎng)精密解算工作與GPS時(shí)間序列解算類(lèi)似,分2步進(jìn)行。第一步為GAMIT解算,由于同步觀測(cè)的測(cè)點(diǎn)較多,考慮到GAMIT單次解算的測(cè)點(diǎn)有限,一般采用分區(qū)解算方式,此步的主要產(chǎn)出為單日松弛解hfile。第二步解算以單日松弛解hfile作為輸入,融合中國(guó)大陸GPS連續(xù)站和IGS單日松弛解結(jié)果,采用GLOBK軟件或QOCA軟件(Dong et al,1998)進(jìn)行GPS速度解算。由于GPS速度解算至少需要2期GPS觀測(cè),因此統(tǒng)一2期(多期)解算結(jié)果的參考基準(zhǔn)就顯得尤為重要,此時(shí)參考框架點(diǎn)的選取原則與時(shí)間序列類(lèi)似。另外一個(gè)需要關(guān)注的問(wèn)題是同震位移的扣除,如果在地震后進(jìn)行過(guò)GPS觀測(cè),則可以利用觀測(cè)數(shù)據(jù)建立同震位移分布模型,對(duì)同震影響進(jìn)行扣除。

    基于上述處理原則,采用GAMIT/GLOBK軟件解算,得出了1999—2007年共4期和2011—2013年共2期云南地區(qū)GPS速度場(chǎng)結(jié)果??紤]到南北地震帶中南部西側(cè)為中國(guó)大陸最大的相對(duì)穩(wěn)定的華南地塊,在云南地區(qū)的GPS站實(shí)測(cè)速度中,扣除華南地塊的歐拉剛體運(yùn)動(dòng)參數(shù)確定的理論值,獲得該區(qū)域相對(duì)于華南地塊的地殼運(yùn)動(dòng)速度圖,既能夠反映云南地區(qū)屬于青藏地塊區(qū)的強(qiáng)烈變形與穩(wěn)定的華南地塊的差異,也能夠兼顧反映云南地區(qū)主要塊體邊界帶的相對(duì)運(yùn)動(dòng)。筆者還針對(duì)制約多期GPS速度場(chǎng)動(dòng)態(tài)演化特征分析的關(guān)鍵問(wèn)題,利用改進(jìn)擬準(zhǔn)檢定(QUAD)初選指標(biāo)的方法,有效解決多期速度場(chǎng)的基準(zhǔn)偏移對(duì)粗差剔除的影響,給出了相對(duì)華南基準(zhǔn)的速度場(chǎng)結(jié)果(鄒鎮(zhèn)宇等,2015)。圖1給出了由上述GPS觀測(cè)資料解算獲得的1999—2007年及2011—2013年云南地區(qū)相對(duì)于華南地塊的GPS速度場(chǎng)結(jié)果。

    2 地殼相對(duì)運(yùn)動(dòng)及其動(dòng)態(tài)變化特征 從云南地區(qū)總體地殼相對(duì)運(yùn)動(dòng)與變形特征來(lái)看(圖1),青藏高原東部物質(zhì)向邊緣流動(dòng)受到邊界較穩(wěn)定塊體的阻擋,造成在接近垂直青藏地塊邊界帶方向地殼縮短與沿地塊邊界帶方向地殼伸展,兼有扭動(dòng)與剪切變形的特征。根據(jù)速度衰減或方向改變,大致為南東向或南向。作為本區(qū)域重要的塊體邊界帶的小江斷裂帶GPS站速度矢量大小方向在空間上的不斷變化,顯示有較強(qiáng)的相對(duì)運(yùn)動(dòng)與變形。這表明由于其特殊的構(gòu)造環(huán)境及其東側(cè)被穩(wěn)定的地塊阻擋,其對(duì)青藏地塊向東運(yùn)動(dòng)及物質(zhì)東移的吸收十分顯著。云南地區(qū)以近南北向地殼縮短變形、近東西向地殼伸展為主。在最南部南北向800 km以?xún)?nèi)最大縮短量約8 mm/a,而在最南端近東西向地殼伸張量為近10 mm/a。該區(qū)地殼相對(duì)運(yùn)動(dòng)與變形以川滇地塊向南擠出滑移最為顯著,且川滇地塊中部運(yùn)動(dòng)速度最大。南部整體上呈現(xiàn)順時(shí)針扭轉(zhuǎn)運(yùn)動(dòng)。GPS站速度從西到東及從北到南不斷向南、南部偏轉(zhuǎn),使站速度從近東西向運(yùn)動(dòng)偏轉(zhuǎn)到川滇地塊南側(cè)地區(qū)為向南或南西運(yùn)動(dòng)。在南北地震帶南部東邊緣,即川滇地塊與華南地塊邊界帶的GPS站速度在平行地塊邊界帶的分量自西向東衰減顯著,表現(xiàn)為明顯的左旋變形特征,總體左旋扭動(dòng)量約為10~12 mm/a。

    由圖1a可見(jiàn),1999—2007年 GPS速度場(chǎng)特征主要體現(xiàn)在:(1)GPS速度場(chǎng)方向從川滇塊體南部南向逐漸轉(zhuǎn)為在滇南塊體的南西向,反映出了繞喜馬拉雅東構(gòu)造結(jié)做順時(shí)針旋轉(zhuǎn)的變形特征;(2)在川滇塊體東邊界的南段存在顯著的左旋剪切運(yùn)動(dòng),滇西北地區(qū)存在顯著的拉張變形;(3)青藏高原物質(zhì)東流,受穩(wěn)定的華南塊體的阻擋,使得云南地區(qū)形成現(xiàn)有的變形特征,雖然區(qū)域的構(gòu)造十分復(fù)雜,但主要的塊體邊界斷裂帶的相對(duì)運(yùn)動(dòng)特征仍能在速度場(chǎng)中清晰的顯示出來(lái)。由圖1b可見(jiàn),2011—2013年GPS點(diǎn)位的密度顯著增加,但是由于觀測(cè)時(shí)間較短,誤差大于1999—2007年的結(jié)果;速度場(chǎng)結(jié)果不僅反映了塊體邊界帶的相對(duì)運(yùn)動(dòng)特征,也可以反映斷裂帶分段變形

    特征(例如小江斷裂帶的北段和中段的左旋滑動(dòng)速率顯著的大于南段的滑動(dòng)速率);并且,相比于1999—2007年速度場(chǎng)結(jié)果,在滇南塊體及周邊地區(qū),2011—2013年南東向的運(yùn)動(dòng)有顯著的增加。

    3 應(yīng)變率場(chǎng)動(dòng)態(tài)變化特征

    雖然采用相對(duì)運(yùn)動(dòng)速度場(chǎng)可以較直觀地反映地殼變形特征,但就對(duì)地殼變形的定量描述方面來(lái)看,應(yīng)變場(chǎng)的分析更有優(yōu)勢(shì)。應(yīng)變場(chǎng)各種參數(shù)能夠全面表達(dá)變形的不同性質(zhì)與強(qiáng)度,從一般概念講,由于應(yīng)變量是由觀測(cè)點(diǎn)之間的相對(duì)運(yùn)動(dòng)來(lái)決定的,它不依賴(lài)于參考基準(zhǔn)。但目前應(yīng)變計(jì)算有多種方法,不同研究者采用不同的方法用相同來(lái)源的GPS資料所給出的應(yīng)變場(chǎng)圖像也有較明顯差異(石耀霖,朱守彪,2006)。Wu 等(2011)借助模擬實(shí)驗(yàn)數(shù)據(jù)對(duì)最小二乘配置應(yīng)變場(chǎng)解算方法結(jié)果與三角單元法、多面函數(shù)擬合和球諧函數(shù)展開(kāi)這幾種應(yīng)變場(chǎng)方法解算結(jié)果進(jìn)行了比較分析,根據(jù)其與理論值的逼近程度,客觀估計(jì)不同方法模型的優(yōu)缺點(diǎn)及其抗差性與實(shí)用性。結(jié)果表明,最小二乘方法最優(yōu)。因此本文利用最小二乘配置球面應(yīng)變率計(jì)算方法對(duì)研究區(qū)域的應(yīng)變率場(chǎng)特征進(jìn)行分析(江在森,劉經(jīng)南,2010)。

    圖2給出了2個(gè)時(shí)段的面應(yīng)變率及最大主應(yīng)變率的分布圖,1999—2007年的結(jié)果顯示:(1)拉張變形主要集中在紅河斷裂帶、程海斷裂帶、麗江小金河斷裂帶南段等區(qū)域,擠壓變形主要集中在滇西及川滇塊體東邊界南段的區(qū)域,并且整體上呈現(xiàn)“中部拉張、兩端擠壓”的變形特征;(2)拉張變形的區(qū)域與地質(zhì)結(jié)果給出的斷裂帶變形特征具有一致性(劉光勛等,1986;黃小巾等,2014),擠壓變形特征可能與此區(qū)域的動(dòng)力源有關(guān),即青藏高原物質(zhì)東流受到穩(wěn)定的華南塊體阻擋,在川滇塊體東邊界南段(小江斷裂帶)形成擠壓區(qū);(3)在滇南的邊界區(qū)域形成了擠壓變形的高值區(qū),由于最小二乘配置方法在計(jì)算的過(guò)程中存在邊界效應(yīng),因此,高擠壓變形區(qū)可信度不高。2011—2013年的結(jié)果一方面繼承了上一期整體變形特征,另一方面還存在著顯著的差異性,主要表現(xiàn)為:拉張變形區(qū)有所擴(kuò)大,并且有向東和向南擴(kuò)展的趨勢(shì);在靠近則木河斷裂帶附近存在顯著的擠壓變形的高值區(qū),川滇塊體東邊界(小江斷裂帶)整體的擠壓變形已經(jīng)不再顯著,主要集中在小江斷裂帶北段靠北的區(qū)域。

    圖3給出了最大剪應(yīng)變率的結(jié)果,1999—2007年的空間分布最顯著的特征是應(yīng)變率高值主要沿川滇地塊東邊界帶南段(小江斷裂帶)分布。并且以程海斷裂帶為界,形成了最大剪應(yīng)變率的高值與低值的過(guò)渡區(qū)域,可能預(yù)示著該區(qū)域剪切應(yīng)力積累水平較高,需關(guān)注該區(qū)域的強(qiáng)震危險(xiǎn)性。另外,區(qū)域剪切變形的分布特征與區(qū)域主要斷裂帶的背景運(yùn)動(dòng)及變形特征密切相關(guān)(例如左旋走滑運(yùn)動(dòng)為主的小江斷裂帶處于最大剪應(yīng)變率的高值區(qū),右旋走滑運(yùn)動(dòng)為主的紅河斷裂帶處于最大剪應(yīng)變率的低值區(qū))。從2011—2013年的空間分布特征來(lái)看,與上一時(shí)段的一致性較好,這也表明了地殼變形隨著時(shí)間的推移,具有較好的繼承性特征;另外,相比于上一時(shí)段,程海斷裂帶附近的最大剪應(yīng)變的高值區(qū)有向北擴(kuò)展的趨勢(shì);與其同步的還有紅河斷裂帶的最大剪應(yīng)變的高值區(qū)也有向北擴(kuò)展的趨勢(shì)。

    4 主要斷裂帶運(yùn)動(dòng)特性研究

    由于最小二乘配置應(yīng)變率計(jì)算方法基于連續(xù)變形假設(shè),因此其結(jié)果在識(shí)別具體斷裂帶的具體變形模式方面不具優(yōu)勢(shì)。GPS速度剖面圖像可直觀展現(xiàn)位移分布與斷層的關(guān)系,經(jīng)常被應(yīng)用于斷層變形特征的分析中(張培震等,2005,2008;Maurin et al,2010;Wu et al,2015)。因此,本文以小江斷裂帶和紅河斷裂帶為例分析斷裂帶運(yùn)動(dòng)變形的動(dòng)態(tài)演化特征。

    4.1 小江斷裂帶

    小江斷裂帶北起巧家以北,南至建水東南,全長(zhǎng)約400 km,總體呈南北走向。根據(jù)其內(nèi)部結(jié)構(gòu)分成3段,分別為北段、中段和南段。在長(zhǎng)期的活動(dòng)過(guò)程中,曾經(jīng)歷壓、張、扭不同力學(xué)性質(zhì)的轉(zhuǎn)化,晚第四紀(jì)以來(lái),斷裂以強(qiáng)烈左旋走滑為特征(宋方敏,汪一鵬,1998)。按照上述的分段原則,以斷裂帶為中心,左右兩側(cè)分別圈定一定范圍,利用該范圍內(nèi)的GPS實(shí)測(cè)速度場(chǎng)數(shù)據(jù)進(jìn)行高斯投影,求得各GPS點(diǎn)的沿?cái)嗔褞ё呦蚍较蚣按怪睌嗔褞ё呦蚍较虻乃俾?;并按照走滑斷層的震間變形模式對(duì)剖面進(jìn)行擬合,最后給出了2個(gè)時(shí)間段平行斷裂帶GPS速度剖面和擬合結(jié)果,如圖4所示。

    由圖4可見(jiàn):(1)小江斷裂帶北段2個(gè)時(shí)段的滑動(dòng)速率變化不大,均為10 mm/a,且表現(xiàn)為一定的應(yīng)變積累的變形特征,變形寬度較寬;(2)小江斷裂帶中段:相比于1999—2007年,2011—2013年的遠(yuǎn)場(chǎng)滑動(dòng)速率有所增加(從7 mm/a增加到10 mm/a),小江斷裂帶從中段開(kāi)始分為2支,現(xiàn)有的GPS點(diǎn)位密度還不能識(shí)別出東西2支的變形特征,所得到的滑動(dòng)速率結(jié)果應(yīng)為2支滑動(dòng)速率之和。雖然滑動(dòng)速率有所增加,并且從剖面的形態(tài)來(lái)看仍然表現(xiàn)為一定的應(yīng)變積累特征,但2個(gè)時(shí)段的變形寬度變化不大;相比于北段,中段的變形寬度較窄;(3)小江斷裂帶南段:相比于北段,該段的滑動(dòng)速率有了明顯的減少,約為5 mm/a,這可能與曲江石屏斷裂帶吸收了部分變形有關(guān),從2個(gè)時(shí)段的動(dòng)態(tài)演化特征來(lái)看,變形寬度沒(méi)有顯著的變化,且變形寬度比小江斷裂帶北段窄。綜合上述分析可以看出,小江斷裂帶以左旋走滑運(yùn)動(dòng)為主,不同段落變形有著明顯的差異,滑動(dòng)速率也有所不同;根據(jù)走滑斷層的震間期的變形模型推測(cè)小江斷裂帶的應(yīng)變積累水平較高;另外,從各段落的動(dòng)態(tài)演化特征來(lái)看,小江斷裂帶不同段落在不同時(shí)段變化并不明顯,汶川地震同震及震后調(diào)整對(duì)此斷裂帶影響較小。

    4.2 紅河斷裂帶中段

    根據(jù)紅河斷裂帶幾何結(jié)構(gòu)特征,可以分為3個(gè)變形區(qū),并且北西段和南東段存在多個(gè)次級(jí)斷裂(向宏發(fā)等,2004),顧及到實(shí)際GPS點(diǎn)位分布特征,本文僅對(duì)紅河斷裂的中段進(jìn)行GPS剖面分析。由于紅河斷裂帶是以右旋走滑為主,因此,圖5給出了平行斷裂帶的GPS速度剖面和擬合結(jié)果。

    從GPS速度剖面及擬合結(jié)果可以看出,紅河斷裂帶中段主要以右旋走滑運(yùn)動(dòng)為主,滑動(dòng)速率約為4 mm/a,2個(gè)時(shí)段的滑動(dòng)速率變化不大;從斷裂帶的變形形態(tài)來(lái)看,存在著一定的應(yīng)變積累特征,以現(xiàn)有的GPS點(diǎn)位做剖面擬合可以看出,2個(gè)時(shí)段的變形寬度變化不大,但變形寬度較寬,可能預(yù)示著該段落可能處于孕震晚期;從以上分析也可看出,汶川地震的發(fā)生對(duì)此斷裂的變形影響較小。

    5 討論與結(jié)論

    基于1999—2007年和2011—2013年的2期GPS速度場(chǎng)資料,并結(jié)合最小二乘配置方法給出的2個(gè)時(shí)段應(yīng)變參數(shù)結(jié)果,綜合分析了汶川地震前后云南地區(qū)地殼變形動(dòng)態(tài)演化特征。獲得了以下主要認(rèn)識(shí):

    (1)2個(gè)時(shí)段GPS速度場(chǎng)的動(dòng)態(tài)結(jié)果表明:整體上GPS速度場(chǎng)方向從川滇塊體南部南向逐漸轉(zhuǎn)為在滇南塊體的南西向,反映出了繞喜馬拉雅東構(gòu)造結(jié)做順時(shí)針旋轉(zhuǎn)的變形特征。塊體邊界斷裂帶附近存在明顯的相對(duì)運(yùn)動(dòng)特征。2011—2013年點(diǎn)位的密度顯著增加,速度場(chǎng)結(jié)果不僅反映了塊體邊界帶的相對(duì)運(yùn)動(dòng)特征,也可以反映斷裂帶分段變形特征;相比于1999—2007年的GPS速度場(chǎng)結(jié)果,在滇南塊體及周邊地區(qū),南東向的運(yùn)動(dòng)有顯著的增加。

    (2)2個(gè)時(shí)段的應(yīng)變率場(chǎng)動(dòng)態(tài)結(jié)果表明,整體上“中部拉張、兩端擠壓”;區(qū)域變形(剪切、拉張或擠壓)的分布特征與區(qū)域主要斷裂帶的背景運(yùn)動(dòng)及變形特征密切相關(guān);相比于1999—2007年,2011—2013年給出的拉張變形區(qū)有所擴(kuò)大,并且有向東和向南擴(kuò)展的趨勢(shì);川滇塊體東邊界(小江斷裂帶)整體的擠壓變形已經(jīng)不再顯著,主要集中在小江斷裂帶北段靠北的區(qū)域;程海斷裂帶附近的最大剪應(yīng)變的高值區(qū)有向北擴(kuò)展的趨勢(shì);與其同步的還有紅河斷裂帶的最大剪應(yīng)變的低值區(qū)也有向北擴(kuò)展的趨勢(shì)。

    (3)跨斷裂帶GPS速度剖面結(jié)果顯示:小江斷裂帶左旋滑動(dòng)速率從北向南逐漸減?。ㄓ?0 mm/a減至5 mm/a),北段變形寬度較寬;紅河斷裂帶右旋滑動(dòng)速率約為4 mm/a,變形寬度較寬;動(dòng)態(tài)結(jié)果表明,汶川地震對(duì)2條斷裂帶變形模式影響較小。

    GPS 數(shù)據(jù)處理得到中國(guó)地震局第一監(jiān)測(cè)中心武艷強(qiáng)研究員數(shù)據(jù)解算團(tuán)隊(duì)的技術(shù)支持,GPS的變形分析得到了中國(guó)地震局地震預(yù)測(cè)研究所魏文薪博士的有益討論,在此一并表示感謝。

    參考文獻(xiàn):

    丁開(kāi)華,許才軍,鄒蓉,等.2013.利用GPS分析川滇地區(qū)活動(dòng)地塊運(yùn)動(dòng)與應(yīng)變模型[J].武漢大學(xué)學(xué)報(bào)(信息科學(xué)版),38(7):822-827.

    方穎,江在森,顧國(guó)華.2008.昆侖山MS8.1地震對(duì)川滇地區(qū)地殼運(yùn)動(dòng)影響分析[J].大地測(cè)量與地球動(dòng)力學(xué),28(4):25-30.

    洪敏,張勇,邵德盛,等.2014.云南地區(qū)近期地殼活動(dòng)特征[J].地震研究,37(3):367-372.

    黃小巾,吳中海,李家存,等.2014.滇西北裂陷帶的構(gòu)造地貌特征與第四紀(jì)構(gòu)造活動(dòng)性[J].地質(zhì)通報(bào),33(4):578-593.

    季靈運(yùn),劉立煒,郝明.2015.利用InSAR技術(shù)研究滇西南鎮(zhèn)康—永德地區(qū)現(xiàn)今地殼形變特征[J].地震研究,38(1):84-89.

    江在森,劉經(jīng)南.2010.應(yīng)用最小二乘配置建立地殼運(yùn)動(dòng)速度場(chǎng)與應(yīng)變場(chǎng)的方法[J].地球物理學(xué)報(bào),53(5):1109-1117.

    江在森,張希,崔篤信,等.2001.青藏塊體東北緣近期水平運(yùn)動(dòng)與變形[J].地球物理學(xué)報(bào),44(5):636-644.

    李西.2015.川滇地塊云南地區(qū)不同發(fā)育階段邊界斷裂破裂特征研究[D].北京:中國(guó)地震局地質(zhì)研究所.

    劉光勛,李方全,李桂榮.1986.我國(guó)滇西北地震活動(dòng)區(qū)的活動(dòng)構(gòu)造與應(yīng)力狀態(tài)[J].地震地質(zhì),8(1):1-14.

    石耀霖,朱守彪.2006.用GPS位移資料計(jì)算應(yīng)變方法的討論[J].大地測(cè)量與地球動(dòng)力學(xué),26(1):1-8.

    宋方敏,汪一鵬.1998.中國(guó)活斷層研究專(zhuān)輯——小江活動(dòng)斷裂帶[M].北京:地震出版社,257.

    王伶俐,邵德盛,施發(fā)奇.2015.基于GPS的云南地區(qū)活動(dòng)地塊現(xiàn)今運(yùn)動(dòng)及應(yīng)變特征分析[J].地震地磁觀測(cè)與研究,36(1):28-36.

    王義昭,熊家鏞,林堯明.1988.云南地質(zhì)構(gòu)造的若干特點(diǎn)[J].云南地質(zhì),7(2):3-9.

    魏文薪,江在森,武艷強(qiáng),等.2012.利用GPS數(shù)據(jù)研究川滇塊體東邊界主要斷裂帶運(yùn)動(dòng)特性[J].武漢大學(xué)學(xué)報(bào)(信息科學(xué)版),37(9):1041-1044.

    武艷強(qiáng),江在森,楊國(guó)華,等.2012.南北地震帶北段近期地殼變形特征研究[J].武漢大學(xué)學(xué)報(bào)(信息科學(xué)版),37(9):1045-1048.

    向宏發(fā),韓竹軍,虢順民,等.2004.紅河斷裂帶大型右旋走滑運(yùn)動(dòng)與伴生構(gòu)造地貌變形[J].地震地質(zhì),26(4):597-610.

    楊國(guó)華,韓月萍,楊博.2009.川滇地區(qū)地殼水平運(yùn)動(dòng)與變形場(chǎng)的演化特征及其機(jī)制討論[J].地震研究,32(3):275-282.

    楊國(guó)華,王琪,王敏,等.2003.云南地區(qū)現(xiàn)今地殼水平運(yùn)動(dòng)與變形特征[J].大地測(cè)量與地球動(dòng)力學(xué),23(2):7-14.

    張俊昌,姚孝執(zhí),吳伯黔.1980.云南主要構(gòu)造體系與強(qiáng)震分布規(guī)律的初步分析[J].地震研究,3(1):47-57.

    張培震,甘衛(wèi)軍,沈正康,等.2005.中國(guó)大陸現(xiàn)今構(gòu)造作用的地塊運(yùn)動(dòng)和連續(xù)變形耦合模型[J].地質(zhì)學(xué)報(bào),79(6):748-756.

    張培震,徐錫偉,聞學(xué)澤,等.2008.2008年汶川8.0級(jí)地震發(fā)震斷裂的滑動(dòng)速率、復(fù)發(fā)周期和構(gòu)造成因[J].地球物理學(xué)報(bào),51(4):1066-1073.

    張希,江在森,王琪,等.2003.川滇地區(qū)地殼水平運(yùn)動(dòng)特征與強(qiáng)震關(guān)系研究[J].大地測(cè)量與地球動(dòng)力學(xué),23(3):35-41.

    鄒鎮(zhèn)宇,江在森,武艷強(qiáng),等.2015.基于GPS速度場(chǎng)變化結(jié)果研究汶川地震前后南北地震帶地殼運(yùn)動(dòng)動(dòng)態(tài)特征[J].地球物理學(xué)報(bào),58(5):1597-1609.

    周海濤,陳兆輝,朱爽,等.2017.川滇地區(qū)部分連續(xù)GPS測(cè)點(diǎn)形變異常分析[J].地震研究,40(1):94-100.

    Dong D,Hering T A,King R W.1998.Estimating regional deformation from a combination of space and terrestrial geodetic data[J].J Geophys Res,72(4):200-214.

    Maurin T,Masson F,Rangin C,et al.2010.First global positioning system results in northern Myanmar:Constant and localized slip rate along the Sagaing fault[J].Geology,38(7):591-594.

    Wu Y Q,Jiang Z S,Yang G H.et al.2011.Comparison of GPS strain rate computing methods and their reliability[J].Geophys J Int,185(2):703-717,doi:10.1111/j.1365-246X.2011.04976.x.

    Wu Y Q,Jiang Z S,Zhao J,et al.2015.Crustal Deformation before the 2008 Wenchuan MS8.0 Earthquake studied using GPS Data[J].J Geodyn,(85):11-23,doi:10.1016/j.jog.2014.12.002.

    Abstract Based on the GPS velocity field data of 1999—2007 and 2011—2013,we used the least squares configuration method and GPS velocity profile results to synthetically analyze the dynamic evolution characteristics of crustal deformation in Yunnan area before and after the Wenchuan earthquake. The dynamic evolution of GPS velocity field shows that the direction is gradually changed from the south of the southern part of the Sichuan-Yunnan block to the south-west of the southern Yunnan block and there is a clear relative motion characteristic near the block boundary fault zone. Compared with the GPS velocity of 1999—2007,the results of 2011—2013 also reflect segmental deformation characteristics of the block boundary fault zone. South-east movement has a significant increase. This may be related to crustal deformation adjustment after the Wenchuan earthquake. The dynamic evolution of strain parameters shows that ‘middle tension,squeeze at both ends is in the whole area and the distribution of deformation(shear,tension or extrusion)is closely related to the background motion and deformation characteristics of the main fault zone. Compared with results during the period of 1999—2007,the tensile zone in 2011—2013 has been expanded in the eastward and southward direction. The extrusion deformation of the eastern boundary(Xiaojiang fault zone)of the Sichuan-Yunnan block is no longer significant,mainly concentrated in the northern section of the Xiaojiang fault zone,which may be related to the post-seismic deformation adjustment of the Wenchuan earthquake. The GPS velocity profile results show that the left-lateral slip velocity of the Xiaojiang fault zone reduced gradually from north to south(from 10 mm/a to 5 mm/a),and the width of the northern section is wider. The right-hand slip rate of the Honghe fault zone is about 4 mm/a,and the deformation width is wider. The dynamic results show that the Wenchuan earthquake has little effect on the deformation modes of the two fault zones.

    Keywords:least squares configuration;GPS velocity field;strain parameter;GPS velocity profile

    99香蕉大伊视频| 丝袜脚勾引网站| 免费av中文字幕在线| 大码成人一级视频| 搡女人真爽免费视频火全软件| 亚洲伊人久久精品综合| 国产白丝娇喘喷水9色精品| 国产午夜精品一二区理论片| av在线app专区| 久久精品久久久久久久性| 侵犯人妻中文字幕一二三四区| 亚洲精品av麻豆狂野| 久久热在线av| 最近中文字幕高清免费大全6| 啦啦啦中文免费视频观看日本| 国产一区亚洲一区在线观看| 亚洲,欧美精品.| 国产精品女同一区二区软件| 晚上一个人看的免费电影| 午夜激情久久久久久久| 国产一区二区三区综合在线观看 | 国产欧美亚洲国产| 国产精品麻豆人妻色哟哟久久| 亚洲图色成人| 亚洲成人一二三区av| 韩国高清视频一区二区三区| 色婷婷久久久亚洲欧美| 亚洲色图综合在线观看| 狂野欧美激情性xxxx在线观看| 99久久综合免费| 大片电影免费在线观看免费| 日韩制服骚丝袜av| 亚洲欧美成人综合另类久久久| av免费观看日本| 欧美bdsm另类| 少妇被粗大猛烈的视频| 精品国产乱码久久久久久小说| av在线观看视频网站免费| a级毛片在线看网站| 亚洲精品自拍成人| 最新的欧美精品一区二区| 一级,二级,三级黄色视频| 夫妻性生交免费视频一级片| 国产不卡av网站在线观看| 热99久久久久精品小说推荐| 精品一区二区三卡| 久久狼人影院| 亚洲三级黄色毛片| 男人添女人高潮全过程视频| 91成人精品电影| 一个人免费看片子| 精品国产一区二区三区四区第35| 精品亚洲成国产av| 午夜精品国产一区二区电影| 黄网站色视频无遮挡免费观看| av女优亚洲男人天堂| 国产老妇伦熟女老妇高清| 女性被躁到高潮视频| 少妇人妻精品综合一区二区| 曰老女人黄片| 国产精品99久久99久久久不卡 | 欧美日本中文国产一区发布| 黄色配什么色好看| 色94色欧美一区二区| 久久久久久久久久久免费av| 亚洲欧美成人精品一区二区| 精品人妻在线不人妻| 成人午夜精彩视频在线观看| 久久久精品免费免费高清| 亚洲三级黄色毛片| 免费大片黄手机在线观看| 亚洲国产看品久久| 少妇的逼水好多| 亚洲精品一二三| 日韩,欧美,国产一区二区三区| 人妻一区二区av| 伦理电影大哥的女人| 久久精品人人爽人人爽视色| 国产精品久久久久久久电影| 精品熟女少妇av免费看| 男女啪啪激烈高潮av片| 欧美亚洲 丝袜 人妻 在线| 只有这里有精品99| 欧美亚洲 丝袜 人妻 在线| 日韩中字成人| 一级毛片 在线播放| 色婷婷av一区二区三区视频| 色吧在线观看| 日韩中文字幕视频在线看片| 午夜激情久久久久久久| 中文字幕最新亚洲高清| 在线观看国产h片| 啦啦啦中文免费视频观看日本| 一二三四中文在线观看免费高清| 亚洲综合精品二区| 欧美日韩av久久| 美女国产高潮福利片在线看| 黄色视频在线播放观看不卡| 侵犯人妻中文字幕一二三四区| 国产精品久久久久久精品古装| 啦啦啦啦在线视频资源| 内地一区二区视频在线| 一个人免费看片子| 国产又色又爽无遮挡免| 成人亚洲欧美一区二区av| 亚洲激情五月婷婷啪啪| 国产熟女午夜一区二区三区| 超色免费av| 少妇被粗大猛烈的视频| 777米奇影视久久| 人人澡人人妻人| av卡一久久| 国产在视频线精品| 亚洲国产色片| 97精品久久久久久久久久精品| 人妻系列 视频| 在线亚洲精品国产二区图片欧美| 国产爽快片一区二区三区| tube8黄色片| 午夜福利,免费看| 国产男女超爽视频在线观看| 亚洲欧美精品自产自拍| 欧美精品一区二区免费开放| 亚洲av男天堂| 国产精品一区二区在线不卡| 高清不卡的av网站| 欧美人与性动交α欧美软件 | 中文字幕精品免费在线观看视频 | 十八禁网站网址无遮挡| 欧美变态另类bdsm刘玥| 亚洲久久久国产精品| 一级毛片黄色毛片免费观看视频| 天天躁夜夜躁狠狠久久av| 哪个播放器可以免费观看大片| 我的女老师完整版在线观看| 欧美精品国产亚洲| 综合色丁香网| 精品一区二区三区视频在线| 亚洲成人一二三区av| 丝袜喷水一区| 亚洲精品一二三| 免费观看在线日韩| 日产精品乱码卡一卡2卡三| 成人国产麻豆网| 大片电影免费在线观看免费| xxxhd国产人妻xxx| 老司机影院成人| 国产一级毛片在线| 巨乳人妻的诱惑在线观看| 欧美老熟妇乱子伦牲交| 人体艺术视频欧美日本| 免费观看在线日韩| 亚洲精品自拍成人| 成人手机av| 日本wwww免费看| 男女边吃奶边做爰视频| 国产av国产精品国产| 欧美国产精品va在线观看不卡| 国产69精品久久久久777片| 日本91视频免费播放| 大香蕉久久网| 91精品三级在线观看| 国产午夜精品一二区理论片| 亚洲欧美成人精品一区二区| 人体艺术视频欧美日本| 免费高清在线观看视频在线观看| 爱豆传媒免费全集在线观看| 国产精品欧美亚洲77777| 久久人人爽av亚洲精品天堂| 精品少妇内射三级| 曰老女人黄片| 久热久热在线精品观看| 午夜福利视频在线观看免费| 97超碰精品成人国产| 99国产精品免费福利视频| 一级毛片黄色毛片免费观看视频| 这个男人来自地球电影免费观看 | 97精品久久久久久久久久精品| 香蕉精品网在线| 一级爰片在线观看| 欧美日本中文国产一区发布| 免费观看av网站的网址| 日韩人妻精品一区2区三区| 国产黄频视频在线观看| 色吧在线观看| 国产av国产精品国产| 亚洲欧美一区二区三区国产| 自线自在国产av| 久久国产精品男人的天堂亚洲 | 国产片内射在线| 热re99久久国产66热| 欧美精品亚洲一区二区| 亚洲美女搞黄在线观看| 日韩在线高清观看一区二区三区| 人体艺术视频欧美日本| 一级,二级,三级黄色视频| 丝袜人妻中文字幕| 国产成人免费观看mmmm| 国产成人免费无遮挡视频| 99热这里只有是精品在线观看| 午夜影院在线不卡| 国产在线免费精品| 久久99蜜桃精品久久| 亚洲欧洲日产国产| 国产在视频线精品| 欧美 日韩 精品 国产| 国产欧美亚洲国产| 中文欧美无线码| 又粗又硬又长又爽又黄的视频| 久久ye,这里只有精品| av卡一久久| 亚洲成人一二三区av| 人妻 亚洲 视频| 日本-黄色视频高清免费观看| 成人黄色视频免费在线看| 亚洲成国产人片在线观看| 亚洲精品av麻豆狂野| 一级黄片播放器| 国产乱人偷精品视频| 欧美激情 高清一区二区三区| 欧美最新免费一区二区三区| 亚洲精品,欧美精品| 一本—道久久a久久精品蜜桃钙片| 最近中文字幕高清免费大全6| 在线 av 中文字幕| 亚洲精品成人av观看孕妇| 国产精品麻豆人妻色哟哟久久| 久久鲁丝午夜福利片| 国产福利在线免费观看视频| av福利片在线| 成人影院久久| 国产黄色视频一区二区在线观看| 欧美精品高潮呻吟av久久| 国产不卡av网站在线观看| 日韩,欧美,国产一区二区三区| 亚洲精品视频女| av又黄又爽大尺度在线免费看| av卡一久久| 另类精品久久| 久久久精品免费免费高清| 免费女性裸体啪啪无遮挡网站| 久久精品久久久久久噜噜老黄| 人人澡人人妻人| 久久人人爽人人片av| 亚洲丝袜综合中文字幕| 一级毛片黄色毛片免费观看视频| 亚洲欧洲国产日韩| 亚洲精品久久午夜乱码| 国产黄色视频一区二区在线观看| 午夜影院在线不卡| 国产精品久久久久久久电影| 亚洲欧洲国产日韩| 一区二区三区四区激情视频| 成人二区视频| 日韩在线高清观看一区二区三区| 男的添女的下面高潮视频| 777米奇影视久久| 欧美日韩国产mv在线观看视频| 超碰97精品在线观看| 国产精品人妻久久久久久| 午夜日本视频在线| 国产视频首页在线观看| 51国产日韩欧美| 久热这里只有精品99| 狠狠婷婷综合久久久久久88av| 99久久中文字幕三级久久日本| 中文字幕av电影在线播放| 啦啦啦啦在线视频资源| 日韩在线高清观看一区二区三区| 2022亚洲国产成人精品| 精品一区在线观看国产| 久久精品夜色国产| 99久久人妻综合| 国产成人精品在线电影| 看非洲黑人一级黄片| 国产免费视频播放在线视频| 成年av动漫网址| 欧美激情极品国产一区二区三区 | 久久精品夜色国产| 久久毛片免费看一区二区三区| 亚洲激情五月婷婷啪啪| √禁漫天堂资源中文www| 午夜激情av网站| 精品一品国产午夜福利视频| 大话2 男鬼变身卡| 精品午夜福利在线看| 中文字幕精品免费在线观看视频 | www.色视频.com| 人妻人人澡人人爽人人| 乱码一卡2卡4卡精品| 美女脱内裤让男人舔精品视频| 七月丁香在线播放| 欧美日韩国产mv在线观看视频| 高清在线视频一区二区三区| 欧美日韩一区二区视频在线观看视频在线| 一本—道久久a久久精品蜜桃钙片| 婷婷色综合www| 国产精品久久久av美女十八| 亚洲欧美精品自产自拍| 久久人人爽人人爽人人片va| 最近的中文字幕免费完整| 久久毛片免费看一区二区三区| 肉色欧美久久久久久久蜜桃| 99热网站在线观看| 欧美日韩国产mv在线观看视频| 伊人久久国产一区二区| 国产成人精品在线电影| 七月丁香在线播放| 日韩免费高清中文字幕av| 午夜福利视频在线观看免费| 精品国产一区二区三区久久久樱花| 一区二区av电影网| 亚洲精品av麻豆狂野| 国产在线一区二区三区精| 插逼视频在线观看| 菩萨蛮人人尽说江南好唐韦庄| 天天影视国产精品| 中文字幕另类日韩欧美亚洲嫩草| 精品人妻在线不人妻| 精品视频人人做人人爽| 免费观看av网站的网址| 又粗又硬又长又爽又黄的视频| 欧美国产精品va在线观看不卡| 在线观看一区二区三区激情| 久久女婷五月综合色啪小说| 日韩大片免费观看网站| 少妇 在线观看| 亚洲,欧美精品.| 大香蕉久久成人网| 午夜免费观看性视频| 午夜福利在线观看免费完整高清在| 亚洲av福利一区| 熟女人妻精品中文字幕| 久热这里只有精品99| 久久99蜜桃精品久久| 久久精品久久精品一区二区三区| 丝袜在线中文字幕| 两个人免费观看高清视频| 国产免费现黄频在线看| 精品一区二区免费观看| 最后的刺客免费高清国语| 国产精品久久久久久精品古装| 五月开心婷婷网| 一本色道久久久久久精品综合| 久久久久网色| 欧美老熟妇乱子伦牲交| 99香蕉大伊视频| 又黄又粗又硬又大视频| 一区二区日韩欧美中文字幕 | 日本与韩国留学比较| 高清不卡的av网站| 男女午夜视频在线观看 | 亚洲av成人精品一二三区| 久久午夜福利片| 国产亚洲av片在线观看秒播厂| 肉色欧美久久久久久久蜜桃| 男女边吃奶边做爰视频| 韩国高清视频一区二区三区| 久久精品国产亚洲av涩爱| 日韩不卡一区二区三区视频在线| 性色avwww在线观看| 一个人免费看片子| 久久免费观看电影| 制服人妻中文乱码| 亚洲国产成人一精品久久久| 国产av精品麻豆| 日韩中文字幕视频在线看片| av女优亚洲男人天堂| 一区二区日韩欧美中文字幕 | 久久影院123| 欧美成人精品欧美一级黄| 黑丝袜美女国产一区| 啦啦啦在线观看免费高清www| a 毛片基地| 97超碰精品成人国产| 日本wwww免费看| 欧美 亚洲 国产 日韩一| 97超碰精品成人国产| av播播在线观看一区| 色视频在线一区二区三区| 亚洲美女搞黄在线观看| 免费观看在线日韩| 亚洲精品国产色婷婷电影| 久热久热在线精品观看| 免费高清在线观看日韩| 亚洲av男天堂| 熟女人妻精品中文字幕| 欧美激情国产日韩精品一区| 久久久久久伊人网av| 少妇人妻精品综合一区二区| 国产在线免费精品| 777米奇影视久久| 亚洲精品中文字幕在线视频| 午夜福利在线观看免费完整高清在| 日韩三级伦理在线观看| 建设人人有责人人尽责人人享有的| 亚洲经典国产精华液单| √禁漫天堂资源中文www| 99热这里只有是精品在线观看| 久久女婷五月综合色啪小说| 国产亚洲午夜精品一区二区久久| 国产亚洲一区二区精品| 久久久国产精品麻豆| 中文精品一卡2卡3卡4更新| 美女视频免费永久观看网站| 亚洲经典国产精华液单| 一边亲一边摸免费视频| 欧美成人午夜精品| 人人妻人人澡人人爽人人夜夜| 激情五月婷婷亚洲| 久久精品久久久久久久性| 老司机影院毛片| 欧美国产精品一级二级三级| 晚上一个人看的免费电影| 欧美xxⅹ黑人| 久久青草综合色| 国产淫语在线视频| 观看av在线不卡| 亚洲高清免费不卡视频| 午夜福利乱码中文字幕| 国产熟女午夜一区二区三区| 美女国产高潮福利片在线看| 青春草视频在线免费观看| 高清黄色对白视频在线免费看| 最黄视频免费看| 1024视频免费在线观看| www日本在线高清视频| 青青草视频在线视频观看| 少妇人妻 视频| 日本av免费视频播放| 亚洲av电影在线观看一区二区三区| 欧美成人精品欧美一级黄| 全区人妻精品视频| 国产综合精华液| tube8黄色片| 免费人妻精品一区二区三区视频| 成人无遮挡网站| 国产伦理片在线播放av一区| 亚洲国产av影院在线观看| 少妇 在线观看| 99精国产麻豆久久婷婷| 综合色丁香网| 免费人成在线观看视频色| 精品一区二区免费观看| 国产乱来视频区| 永久网站在线| 国产成人免费无遮挡视频| 日韩在线高清观看一区二区三区| 国产精品国产三级国产av玫瑰| 在线亚洲精品国产二区图片欧美| 欧美精品一区二区大全| 亚洲成人手机| 午夜福利,免费看| 精品午夜福利在线看| 亚洲精品国产av蜜桃| 亚洲av中文av极速乱| av电影中文网址| 日本黄色日本黄色录像| 日本欧美国产在线视频| 18禁动态无遮挡网站| 亚洲一码二码三码区别大吗| 午夜福利视频精品| 日韩在线高清观看一区二区三区| 精品午夜福利在线看| 亚洲图色成人| 两性夫妻黄色片 | 极品少妇高潮喷水抽搐| 欧美日韩成人在线一区二区| 国产精品久久久久久精品电影小说| 波野结衣二区三区在线| 人人妻人人添人人爽欧美一区卜| 久久精品aⅴ一区二区三区四区 | 亚洲精品aⅴ在线观看| 精品少妇久久久久久888优播| 国产淫语在线视频| 美女国产视频在线观看| 亚洲精品456在线播放app| 伊人亚洲综合成人网| 午夜老司机福利剧场| 免费观看a级毛片全部| 少妇熟女欧美另类| 母亲3免费完整高清在线观看 | 极品人妻少妇av视频| 欧美成人午夜精品| 成人毛片60女人毛片免费| 嫩草影院入口| 精品少妇黑人巨大在线播放| 亚洲欧洲精品一区二区精品久久久 | 亚洲av免费高清在线观看| 久久久久国产网址| 国产亚洲最大av| 国产一区有黄有色的免费视频| 又大又黄又爽视频免费| 超碰97精品在线观看| 男女免费视频国产| 国产成人a∨麻豆精品| 亚洲精品美女久久久久99蜜臀 | 丰满乱子伦码专区| 免费少妇av软件| 少妇人妻 视频| 亚洲精品国产av蜜桃| 国产精品 国内视频| 欧美国产精品va在线观看不卡| 免费高清在线观看日韩| 一区二区三区四区激情视频| 人人妻人人爽人人添夜夜欢视频| 丰满乱子伦码专区| 在线观看免费视频网站a站| 一级毛片 在线播放| 九九爱精品视频在线观看| 国产 一区精品| 啦啦啦啦在线视频资源| 美女大奶头黄色视频| 中文字幕精品免费在线观看视频 | 卡戴珊不雅视频在线播放| 国产成人一区二区在线| av一本久久久久| 考比视频在线观看| 一级a做视频免费观看| 午夜影院在线不卡| 女人精品久久久久毛片| 国产在线一区二区三区精| 美女主播在线视频| 亚洲av综合色区一区| freevideosex欧美| 狂野欧美激情性bbbbbb| 天天操日日干夜夜撸| 亚洲成人手机| videosex国产| 精品少妇久久久久久888优播| 色5月婷婷丁香| 国产综合精华液| 99re6热这里在线精品视频| av天堂久久9| 热99久久久久精品小说推荐| 国产男女内射视频| 久久精品熟女亚洲av麻豆精品| 国产欧美另类精品又又久久亚洲欧美| 男女午夜视频在线观看 | 天天躁夜夜躁狠狠躁躁| 哪个播放器可以免费观看大片| 国产亚洲一区二区精品| 大香蕉久久成人网| 韩国av在线不卡| 只有这里有精品99| 青春草国产在线视频| 亚洲av国产av综合av卡| 国产精品99久久99久久久不卡 | a级毛片黄视频| 亚洲国产精品一区二区三区在线| 免费日韩欧美在线观看| 亚洲欧美一区二区三区黑人 | 国内精品宾馆在线| 90打野战视频偷拍视频| 成人无遮挡网站| 肉色欧美久久久久久久蜜桃| av不卡在线播放| 最近中文字幕2019免费版| 欧美日韩国产mv在线观看视频| 亚洲av电影在线进入| 午夜91福利影院| 一级毛片黄色毛片免费观看视频| 人妻人人澡人人爽人人| 免费观看av网站的网址| 亚洲av免费高清在线观看| 在线看a的网站| 五月天丁香电影| 欧美精品一区二区大全| 水蜜桃什么品种好| 永久免费av网站大全| av女优亚洲男人天堂| 日韩一本色道免费dvd| 秋霞在线观看毛片| 国产黄频视频在线观看| 亚洲精品中文字幕在线视频| 男人操女人黄网站| 国产欧美日韩一区二区三区在线| 日韩伦理黄色片| 乱人伦中国视频| 久久久久久伊人网av| 一级爰片在线观看| 高清av免费在线| 成人午夜精彩视频在线观看| 国产精品三级大全| 少妇被粗大的猛进出69影院 | 国产 精品1| 丝袜人妻中文字幕| 欧美亚洲日本最大视频资源| 日韩精品有码人妻一区| 国产在线一区二区三区精| 爱豆传媒免费全集在线观看| 精品视频人人做人人爽| 丰满乱子伦码专区| 久久久久久久久久成人| 久久97久久精品| 亚洲成人av在线免费| 国产精品久久久久成人av| 国产成人精品无人区| 亚洲欧美日韩另类电影网站| 这个男人来自地球电影免费观看 | 菩萨蛮人人尽说江南好唐韦庄| 男女午夜视频在线观看 | 国产一区有黄有色的免费视频| 大片电影免费在线观看免费| 少妇被粗大的猛进出69影院 | 亚洲国产欧美日韩在线播放| 欧美日韩视频精品一区| 国产片特级美女逼逼视频| 在线天堂最新版资源| 丝袜脚勾引网站| 一区二区三区精品91| 亚洲精品456在线播放app| 日韩伦理黄色片| 精品亚洲成a人片在线观看| 欧美精品国产亚洲|