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

    云南地區(qū)重力異常的小波分解與解釋*

    2014-12-25 05:34:16龐衛(wèi)東張建國沈婭宏
    地震研究 2014年4期
    關(guān)鍵詞:重力云南細節(jié)

    龐衛(wèi)東,張建國,沈婭宏

    (云南省地震局,云南昆明650224)

    0 引言

    地質(zhì)歷史各期構(gòu)造作用,伴隨的巖漿活動在造成大陸巖石圈結(jié)構(gòu)橫向非均質(zhì)性的同時,必然會造成地殼內(nèi)部、上地幔物質(zhì)分布不均勻及其密度等物性差異,并表現(xiàn)為相應的重力異常分布。對重力異常場的研究可以幫助了解區(qū)域尺度范圍內(nèi)的地殼、上地幔結(jié)構(gòu),尤其是基底性質(zhì)、莫霍面起伏等。

    云南及攀西 (攀枝花—西昌)地區(qū)處于印度板塊向中國大陸擠壓作用地帶,區(qū)內(nèi)構(gòu)造復雜,深大斷裂發(fā)育,地殼較為破碎,地震活動十分頻繁,為中國大陸構(gòu)造活動最為強烈的地區(qū)之一。自20世紀80年代以來,很多學者已對該地區(qū)的地殼結(jié)構(gòu)進行了探測和研究工作 (闞榮舉等,1977,1986;熊紹柏等,1986;胡鴻翔等,1986;孫潔等,1989;王椿鏞等,2002)。前人也對區(qū)內(nèi)重力資料進行過大量的研究 (闞榮舉,林中洋,1986;何翔,吳國華,1986;蔣航,劉福生,1987;蔣福珍,方劍,2001),得到了莫霍面深度及地殼與上地幔密度變化等很多成果。但由于資料和研究方法的限制,由重力資料得到的一些結(jié)果與地震測深等研究結(jié)果不太一致,因此有必要作進一步地研究。本文運用小波分析方法對云南地區(qū)的重力資料進行處理,以實現(xiàn)位場的有效分解,并結(jié)合區(qū)域地質(zhì)資料及其它地球物理資料對研究區(qū)的地質(zhì)構(gòu)造、地殼結(jié)構(gòu)和地震活動性進行分析研究。

    1 方法

    小波分析作為傅立葉分析思想的發(fā)展,克服了傅立葉分析不能獲得信號局部特性的缺點,使其在信號處理、圖像分析、地球物理等很多學科領(lǐng)域得到越來越廣泛的應用。小波變換應用其多尺度分析原理 (Mallat,1989)可將信號f(x)分解成各種不同的頻道和頻率成分,并且通過伸縮、平移聚焦到f(x)的任一細節(jié)加以分析 (李世雄等,1994)。利用小波變換的特點能將重力異常場精細地分解到多個不同的尺度上,可反映不同尺度和深度的異常,侯遵澤等 (1997,1998)、高德章等(2000)、楊文采等 (2001)先后把小波分析方法應用于重力場的分離和處理,并詳盡闡述了小波變換及多尺度分析的計算方法及原理,本文只給出重力異常多尺度分解后的表達式。

    假設二維重力異常場為 Δg(x,y)=f(x,y),利用小波多尺度分析原理對其分解后的重力異常表達式為

    式中,D1G-DNG為一階至N階小波變換異常細節(jié);ANG為N階小波變換異常逼近。

    在利用小波分析方法處理數(shù)據(jù)時,其基函數(shù)有多種選擇,對于位場數(shù)據(jù),需要考慮小波基函數(shù)的一些重要性質(zhì),李健等 (2001)對此進行了分析與試驗。小波函數(shù)的完備正交性是多尺度分析的基本要求,正交小波可以保證多尺度分析的低階細節(jié)不變性。小波函數(shù)的對稱性直接影響信號的重構(gòu),如果具有對稱性,重構(gòu)算法失真小,重構(gòu)信號能給出原始信號的一個很好的逼近。小波函數(shù)的消失矩決定了小波函數(shù)逼近光滑函數(shù)時的收斂率,較高的消失矩能夠研究函數(shù)的高階變化。由于不存在正交對稱小波,在處理中考慮到研究區(qū)位場的變化較為平緩,選擇小波基函數(shù)時,傾向于選用振蕩較為平緩和光滑的小波基函數(shù),此外選用的小波基函數(shù)還應有較好對稱性、緊支撐性和較高的消失矩。

    圖1 云南布格重力異常圖 (重力單位:mGal)Fig.1 Bougue gravity anomalies in Yunnan area

    2 資料與數(shù)據(jù)處理

    數(shù)據(jù)資料選用云南省地礦局編制的1:100萬布格重力異常數(shù)據(jù),其數(shù)據(jù)范圍為 (22.5°~29.4°N,98.6°~105.5°E)。圖1顯示區(qū)域布格重力異常值全為負,多在 (-80~-400) ×10-5m·s-2變化,宏觀上從滇東南到滇西北隨地勢升高逐漸降低,等值線呈向南東凹入,由南北轉(zhuǎn)東西的弧形展布。重力異常特點表現(xiàn)為:以紅河斷裂為界,斷裂帶兩側(cè)重力場特征迥異。滇西地區(qū)重力場總的呈南高北低,但又沿北東方向呈高低相間的起伏變化,異常以北西向為主,兼有南北向及北東向的異常表現(xiàn);滇西北則表現(xiàn)為北東向梯級帶上疊加近南北向異常。滇中、滇東地區(qū)重力場總趨勢由南東向北西逐步降低,重力高、低異常呈多個中心分布,異常有較大的升降幅度,但帶狀延展不如滇西明顯,異常帶比較寬,異常數(shù)目比滇西少,在26°N以南變化平緩。滇東地區(qū)南盤江以南,高低異常帶呈東西向展布,表現(xiàn)特殊 (云南省地質(zhì)礦產(chǎn)局,1990)。

    利用小波分析方法對布格重力異常數(shù)據(jù)進行處理,參考李健等 (2001)研究成果,選用雙正交bior3.5小波對重力異常數(shù)據(jù)進行6階分解,并對各階小波細節(jié)進行徑向?qū)?shù)功率譜分析,求取其對應的場源深度 (Spector,Grant,1970;楊文采等,2001),結(jié)果如圖2所示。本文僅展示4~6階的細節(jié)和6階逼近圖像。從統(tǒng)計意義上看,重力異常的徑向平均對數(shù)功率譜曲線可以反映異常場源的平均深度。功率譜計算結(jié)果表明,4階細節(jié)的場源埋深為8 km,反映了上地殼密度不均勻體;5階小波細節(jié)場源埋深為15 km,6階小波細節(jié)場源埋深為25 km,5階和6階小波細節(jié)反映了上地殼底部和中地殼的不均勻性;6階小波逼近場源埋深48 km,則反映了莫霍面起伏的趨勢。

    3 布格重力異常場結(jié)果分析及解釋

    對云南地區(qū)布格重力異常場處理結(jié)果進行分析,并結(jié)合研究區(qū)內(nèi)的地質(zhì)及其它地球物理資料進行相應的地質(zhì)解釋,對區(qū)域地質(zhì)構(gòu)造和地震活動性做初步探討。

    3.1 小波分解結(jié)果

    筆者對云南及周邊地區(qū)布格重力異常小波分解各階結(jié)果進行分析。圖2a為布格重力異常4階小波細節(jié),等值線形態(tài)十分復雜,為上地殼淺部小尺度構(gòu)造及巖體的反映。圖2b是5階小波細節(jié),由圖可見異常帶清晰連續(xù),對區(qū)內(nèi)斷裂帶反映的突出,表現(xiàn)為地下物質(zhì)密度的小波細節(jié)異常的位置與地表斷裂是基本吻合的。同時也存在不完全吻合的地方,一方面這也許是可以對地表斷裂做出相應修正的地方,因為地球物理場具有連續(xù)性,有利于確定斷裂的位置,而地表巖石或斷層的露頭出現(xiàn)比較孤立;另一方面考慮到在此深度上的(5階細節(jié)場源深度為15 km)這種不吻合,也許是斷裂向下延伸時傾角發(fā)生變化而造成。在4階、5階細節(jié)圖中除在局部和沿斷裂帶部位有正的重力異常區(qū) (推測正異常是地幔高密度物質(zhì)沿斷裂上涌至地殼淺部所致)外,在這一深度層次上 (4階細節(jié)8 km,5階細節(jié)15 km)以重力負異常為背景,尤其是區(qū)內(nèi)主要的坳陷和湖盆區(qū),如滇中坳陷內(nèi)的楚雄—南華地區(qū)、鎮(zhèn)沅—思茅弧后盆地及昆明新生代湖盆區(qū)等,異常分布與表面地質(zhì)構(gòu)造相關(guān)。這與王椿鏞等 (2002)利用地震記錄獲得的P、S波三維速度結(jié)構(gòu)在地殼10 km深度上區(qū)域內(nèi)速度異常值相對較小,存在楚雄—南華地區(qū)等負異常區(qū)的結(jié)果相吻合。

    云南地區(qū)6階小波細節(jié)如圖2c所示,計算場源深度為25 km,反映的是中地殼密度不均勻分布,由于重力資料對深部場源的分辨率較低,所以該階小波細節(jié)相較5階細節(jié)對斷裂構(gòu)造的反映沒有那么清晰,但該圖像清楚得顯示深部大尺度不均勻塊體。研究區(qū)中部的西昌—攀枝花—楚雄直至紅河斷裂南北向展布的重力正異常區(qū)與攀西裂谷和滇中坳陷區(qū)位置相對應,滕吉文 (1987),熊紹伯等 (1986)研究指出滇中坳陷與攀西裂谷具有相似的殼幔結(jié)構(gòu)和地球物理特征,表明攀西中生代古裂谷實際上已南延至滇中 (闞榮舉,1989),滇中坳陷具備了裂谷的特征。中生代揚子古陸裂陷形成滇中坳陷,發(fā)育了巨厚的中生代沉積 (云南省地礦局,1990),對應4階和5階細節(jié)重力負異常區(qū),而在巨厚的中生代蓋層之下,中地殼內(nèi)存在鎂鐵質(zhì)巖體和基性巖墻群,此應為沿裂谷中的巖漿通道上涌的高密度物質(zhì),對應西昌—攀枝花—楚雄直至紅河斷裂的南北向展布的重力正異常。在該正異常條帶東側(cè)存在著與其平行展布的重力負異常帶,在地質(zhì)構(gòu)造上對應著康滇古隆起;再向東側(cè)出現(xiàn)南北向的異常等值線寬疏的正異常區(qū),小江斷裂東側(cè)的正異常應是牛首山古陸的反映。圖2d顯示了6階小波逼近,功率譜計算場源深度為48 km,反映了莫霍面起伏趨勢,其形態(tài)與根據(jù)地震測深資料編制的地殼等厚線形態(tài)十分相似,并與胡家富等 (2003)利用S波地殼結(jié)構(gòu)反演所得地殼厚度結(jié)果相吻合。重力異常6階逼近圖像中 (圖2d),可以看到一個向南東突出的舌形,在前舌和后舌之間呈現(xiàn)出“一隆兩坳”的特點 (闞榮舉,林中洋,1986),即攀枝花—楚雄重力高對應的滇中幔隆、以東川為中心的重力低對應的滇東幔坳及滇西重力低對應的滇西復式幔坳 (任治磯,1996)。若扣除前舌與后舌隆坳相間的影響,則研究區(qū)內(nèi)地殼厚度變化趨勢將表現(xiàn)出一個更為平滑的向南東突出的舌形,為此,在六階小波逼近的基礎上進一步做7階小波分解,所得結(jié)果如圖3b所示。7階小波細節(jié) (圖3a)則更清楚的反映出研究區(qū)內(nèi)上地幔隆坳相間展布的特點,與“一隆兩坳”的特征相對應,只是其隆坳中心均發(fā)生了不同程度的北向偏移。云南地區(qū)布格重力異常場處理結(jié)果表明云南及周邊地區(qū)地殼、上地幔物質(zhì)存在明顯的橫向不均勻性。

    圖3 云南重力異常7階細節(jié) (a)和7階逼近圖 (b)(重力單位:mGal)Fig.3 The 7th order wavelet details(a)and the 7th order approximation(b)of the gravity anomalies in Yunnan area

    3.2 地質(zhì)構(gòu)造分析

    為了便于分析,筆者把不同階次的小波細節(jié)數(shù)據(jù)組合起來以獲得不同深度的目標場,將5階和6階小波細節(jié)組合得到重力異常5~6階細節(jié)組合,如圖4所示。組合所得結(jié)果可以避免由于斷裂在不同深度層次上因發(fā)育程度的不同而造成在某一階次圖像中不能完整顯示的缺點,同時也剔除了上部密度不均勻體的干擾。從研究區(qū)內(nèi)深大斷裂在組合圖中投影的結(jié)果看,重力異常的轉(zhuǎn)換帶、梯度帶或突變帶多與斷裂展布相吻合。如程海斷裂為重力正異常區(qū)的東邊界,瀾滄江斷裂南段、綠汁江斷裂北段,小江斷裂帶均為重力異常梯級帶,紅河斷裂帶為正負異常轉(zhuǎn)換帶。研究區(qū)內(nèi)異常條帶展布的另一個顯著的特點是大部分斷裂都投影在異常條帶的正上方或與其緊鄰,這表明區(qū)內(nèi)斷裂多具陡傾的水平走滑性質(zhì)。作為“板塊縫合線”的金沙江—紅河斷裂,金沙江斷裂自滇西北進入本研究區(qū),其重力異常顯示,呈現(xiàn)局部相當狹窄的北西向負異常軸,與斷裂的地面位置符合的很好,特征指示該斷裂產(chǎn)狀相當陡立,與地震測深資料震相分析獲得的斷層產(chǎn)狀 (傾角87.5°)相吻合 (闞榮舉,韓源,1992)。南部北西走向的紅河斷裂表現(xiàn)為梯級帶,并在其西南側(cè)伴有蒼山—哀牢山負異常帶,與地質(zhì)構(gòu)造上的蒼山—哀牢山變質(zhì)帶相對應,針對該負異常,結(jié)合Pg波低速異常等現(xiàn)象認為該區(qū)有逆沖推覆構(gòu)造,存在哀牢山推覆體 (陳元坤,1987)。紅河斷裂北東側(cè)是重力正異常區(qū)與高速異常區(qū)對應,這一正異常影響至與其平行的楚雄—建水斷裂。圖4中紅河斷裂與楚雄—建水斷裂異常條帶平行展布且中間夾持正異常軸帶的特征,可能說明兩條斷裂帶實為深部同一大斷裂在地表附近的不同分支,而中間的異常軸對應著深部大斷裂,又聯(lián)系到地表出露的玄武巖(云南省地礦局,1990),認為重力正異常為地幔、下地殼巖漿順著該深部斷裂侵位到中上地殼,并沿著此平行斷裂系向上傳播的反映。本文分析結(jié)果與白志明和王椿鏞 (2003,2004)對地殼結(jié)構(gòu)層析成像研究顯示的紅河大斷裂在淺部傾向南西,并分支為紅河斷裂、建水斷裂和曲江斷裂,在深部倒轉(zhuǎn),傾向北東切穿地殼的結(jié)果吻合。在南瀾滄江與昌寧—雙江斷裂之間存在一規(guī)模較大的負異常區(qū),應該反映了華力西、印支時代北瀾滄江—昌寧雙江斷裂帶向東俯沖形成的具有低壓高溫和高壓低溫的雙變質(zhì)巖系特征的板塊縫合帶,西側(cè)靠近昌寧雙江斷裂的瀾滄群變質(zhì)巖 (劉增池,1985)及靠近瀾滄江斷裂東側(cè)為一個有著復雜形成歷史的混合—改造型臨滄花崗巖體巖基 (熊家鏞,1983),則表現(xiàn)為重力正異常。

    圖4 云南重力異常5~6階細節(jié)組合(重力單位:mGal)Fig.4 The combination of 5th and 6th order wavelet details of the gravity anomalies in Yunnan area

    3.3 地震活動性分析

    云南地區(qū)是強震多發(fā)區(qū),其地震發(fā)生多受活動斷裂控制。對比云南地區(qū)強震分布 (圖4)和地震帶劃分情況 (秦嘉政等,2005),注意到研究區(qū)重力異常條帶和復雜異常區(qū)與地震分布具有一定的相關(guān)性,大多數(shù)強震震中集中在這些區(qū)域。斷裂兩側(cè)異常突變帶、梯度帶反映了斷裂兩側(cè)塊體物性的差異,復雜異常區(qū)則反映不同斷裂條帶的交匯或交錯造成的塊體破碎,而這些構(gòu)造部位往往利于應力集中、能量積聚,為地震孕育和發(fā)生提供了場所。如建水斷裂和曲江斷裂附近的眾多中強地震形成的通海—石屏地震帶,就是斷裂交匯或交錯的部位;小江斷裂異常條帶對應的小江地震帶則是不同活動塊體的連接邊界,該帶上曾發(fā)生過有記載以來云南最大的1833年嵩明8.0級地震;昌寧—雙江斷裂為1988年瀾滄耿馬7.6級地震的發(fā)震構(gòu)造,在異常圖上呈現(xiàn)出的正負異常轉(zhuǎn)換帶對應著耿馬—瀾滄地震帶,從異常等值線形態(tài)可以注意到沿北西向瀾滄—耿馬地震帶直至騰沖—龍陵地震帶為一個大的異常梯度帶,推測可能為一條隱伏構(gòu)造 (圖4中紅色虛線標注)的反映;中甸—麗江—大理區(qū)域異常條帶表明川滇菱形塊體內(nèi)部劍川斷裂、麗江—小金河斷裂為界兩側(cè)塊體 (北西側(cè)為川西北次級塊體、南東側(cè)為滇中次級塊體)具有活動差異性,與中甸—大理地震帶相吻合。

    4 結(jié)論

    利用小波分析方法,分解云南及周邊地區(qū)重力異常得到的局部異常 (細節(jié)成分)和區(qū)域異常(逼近成分)是合理的。通過對重力異常結(jié)果的分析及地質(zhì)解釋,得到了以下認識:局部異常較好地反映了區(qū)內(nèi)主要的斷裂及主要的坳陷去和湖盆區(qū)的構(gòu)造特征;區(qū)域異常則反映出研究區(qū)內(nèi)地殼、上地幔物質(zhì)存在明顯的橫向不均勻性,主要表現(xiàn)為整體向南東突出的舌形,莫霍面存在南北向的隆起和坳陷,呈“一隆兩坳”的特點,在深部隆坳中心有北向偏移的現(xiàn)象。云南作為強震多發(fā)區(qū),其地震發(fā)生多受活動斷裂控制,區(qū)內(nèi)重力異常條帶和復雜異常區(qū)與地震分布具有一定的相關(guān)性,斷裂帶兩側(cè)的異常反映了斷裂兩側(cè)存在明顯的物性、結(jié)構(gòu)差異,而這些差異部位往往有利于應力的相對集中、能量的積累及地震的孕育和發(fā)生。

    白志明,王椿鏞.2003.云南地區(qū)上部地殼結(jié)構(gòu)和地震構(gòu)造環(huán)境的層析成像研究[J].地震學報,23(2):117-127.

    白志明,王椿鏞.2004.云南遮放—賓川和孟連—馬龍寬角地震剖面的層析成像研究[J].地球物理學報,47(2):257-267.

    陳元坤.1987.哀牢山—點蒼山推覆構(gòu)造的認識依據(jù)和討論[J].云南地質(zhì),16(4):291-297.

    高德章,侯遵澤,唐建.2000.東海及鄰區(qū)重力異常多尺度分解[J].地球物理學報,43(6):842-849.

    何翔,吳國華.1986.云南地區(qū)地殼厚度及重力均衡補償作用[J].地震研究,9(4):487-500.

    侯遵澤,楊文采,劉家琦.1998.中國大陸地殼密度差異多尺度反演[J].地球物理學報,41(5):643-651.

    侯遵澤,楊文采.1997.中國重力異常的小波變換與多尺度分析[J].地球物理學報,40(1):85-95.

    胡鴻翔,陸涵行,王椿鏞,等.1986.滇西地區(qū)地殼結(jié)構(gòu)的爆破研究[J].地球物理學報,29(2):133-144.

    胡家富,蘇有錦,朱雄關(guān),等.2003.云南地殼S波速度與泊松比結(jié)構(gòu)及意義[J].中國科學(D 輯),33(8):714-722.

    蔣福珍,方劍.2001.康滇地區(qū)重力場分離,密度反演與地殼構(gòu)造[J].地震學報,23(4):391-397.

    蔣航,劉福生.1987.從重磁場特征探討攀西裂谷構(gòu)造[A]//張云湘,袁學誠.中國攀西裂谷文集(3)[C].北京:地質(zhì)出版社,76-88.

    闞榮舉,韓源.1992.云南遮放至馬龍地學斷面(說明書)[Z].北京:地震出版社.

    闞榮舉,林中洋.1986.云南地殼上地幔構(gòu)造的初步研究[J].中國地震,2(4):50-61.

    闞榮舉,張四昌,晏風桐,等.1977.我國西南地區(qū)現(xiàn)代構(gòu)造應力場與現(xiàn)在構(gòu)造特征的探討[J].地球物理學報,20(2):96-109.

    闞榮舉.1989.攀西古裂谷及南延部分的地球物理場特征[A]//八十年代地區(qū)中國物理學進展[C].北京:學術(shù)書刊出版社,275-283.

    李健,周云軒,許惠平.2001.重力場數(shù)據(jù)處理中小波母函數(shù)的選擇[J].物探與化探,25(6):410-416.

    劉增池.1985.瀾滄群火山巖特征及其構(gòu)造背景的判別[J].云南地質(zhì),4(4):36-45.

    秦嘉政,皇甫崗,錢曉東.2005.云南強震活動與預測方法研究[M].昆明:云南科技出版社.

    任治磯.1996.云南地體構(gòu)造與成礦作用[M].北京:地質(zhì)出版社.

    孫潔,徐常芳,江釗,等.1989.滇西地區(qū)地殼上地幔電性結(jié)構(gòu)與地殼構(gòu)造運動的關(guān)系[J].地震地質(zhì),11(1):35-45.

    滕吉文.1987.攀枝花—西昌古裂谷與“活化”的地球物理場特征[J].地球物理學報,30(6):581-593.

    王椿鏞,Mooney W.D.,王溪莉,等.2002.川滇地區(qū)地殼上地幔三維速度結(jié)構(gòu)研究[J].地震學報,24(1):1-16.

    熊家鏞.1983.臨滄混合雜巖的基本特征成因探討[A]//青藏高原地質(zhì)文集(13)[C].北京:地質(zhì)出版社.

    熊紹柏,滕吉文,尹周勛.1986.攀西構(gòu)造帶南部地殼與上地幔結(jié)構(gòu)的爆炸地震研究[J].地球物理學報,29(3):235-244.

    楊文采,施志群,侯遵澤,等.2001.離散小波變換與重力異常多重分解[J].地球物理學報,44(4):534-541.

    云南省地質(zhì)礦產(chǎn)局.1990.云南省區(qū)域地質(zhì)志[M].北京:地質(zhì)出版社.

    Mallat S..1989.A theory for multiresolution signal decomposition:the wavelet representation[J].IEEE Pattern Anal.and Machine Intell.,11(7):674-693.

    Spector A..Grant F.S..1970.Statistical models for interpreting aeromagnetic data[J].Geophysics,35(2):293-302.

    猜你喜歡
    重力云南細節(jié)
    瘋狂過山車——重力是什么
    科學大眾(2022年23期)2023-01-30 07:04:16
    云南茶,1200年的發(fā)現(xiàn)
    云南畫報(2021年11期)2022-01-18 03:15:40
    云南邀您來“吸氧”
    云南畫報(2020年12期)2021-01-18 07:19:20
    以細節(jié)取勝 Cambridge Audio AXR100/ FOCAL ARIA 906
    云南是你避暑的最佳選擇
    云南畫報(2020年9期)2020-10-27 02:03:16
    一圖讀懂云南兩新黨建
    留心細節(jié)處處美——《收集東·收集西》
    仰斜式重力擋土墻穩(wěn)定計算復核
    細節(jié)取勝
    Coco薇(2016年10期)2016-11-29 19:59:58
    一張紙的承重力有多大?
    精品第一国产精品| 天天影视国产精品| 嫩草影视91久久| 色在线成人网| 老熟女久久久| 男人舔女人的私密视频| 精品一区二区三区视频在线观看免费 | 91大片在线观看| 在线播放国产精品三级| 高清av免费在线| 伦理电影免费视频| 久久久久视频综合| 国产单亲对白刺激| 高清av免费在线| 久久久久精品国产欧美久久久| 久久ye,这里只有精品| 午夜激情久久久久久久| svipshipincom国产片| 欧美另类亚洲清纯唯美| www.精华液| 久久国产精品人妻蜜桃| 黄片播放在线免费| 精品一区二区三区av网在线观看 | 精品一区二区三区视频在线观看免费 | 欧美精品高潮呻吟av久久| 免费看a级黄色片| 搡老岳熟女国产| 亚洲精品美女久久av网站| 精品国产一区二区三区久久久樱花| 丝袜美腿诱惑在线| 一夜夜www| 国产免费福利视频在线观看| 麻豆av在线久日| 久久精品成人免费网站| 交换朋友夫妻互换小说| 夜夜爽天天搞| 午夜福利在线观看吧| 男女边摸边吃奶| av超薄肉色丝袜交足视频| 久久亚洲真实| av网站免费在线观看视频| a级毛片在线看网站| 极品教师在线免费播放| 女性生殖器流出的白浆| 亚洲av日韩精品久久久久久密| 日本黄色日本黄色录像| 精品人妻在线不人妻| 欧美老熟妇乱子伦牲交| 99re6热这里在线精品视频| 在线亚洲精品国产二区图片欧美| 飞空精品影院首页| 狠狠精品人妻久久久久久综合| 亚洲国产看品久久| 国产亚洲欧美在线一区二区| 午夜福利在线免费观看网站| 国产高清videossex| 搡老乐熟女国产| 精品卡一卡二卡四卡免费| 欧美人与性动交α欧美精品济南到| 老司机午夜福利在线观看视频 | 国产一区二区 视频在线| 少妇猛男粗大的猛烈进出视频| 亚洲专区字幕在线| 亚洲一区中文字幕在线| 久久久久久人人人人人| 岛国在线观看网站| 欧美日韩福利视频一区二区| 色尼玛亚洲综合影院| 久久青草综合色| 午夜91福利影院| 国产免费视频播放在线视频| 在线观看舔阴道视频| 9热在线视频观看99| 菩萨蛮人人尽说江南好唐韦庄| 少妇的丰满在线观看| 最新美女视频免费是黄的| 国产1区2区3区精品| 亚洲午夜理论影院| 午夜免费鲁丝| 中文字幕人妻丝袜制服| 精品亚洲成a人片在线观看| 国产精品亚洲一级av第二区| 麻豆乱淫一区二区| 丝袜在线中文字幕| 国产在线观看jvid| 一进一出好大好爽视频| 国产欧美亚洲国产| 老司机影院毛片| 亚洲一码二码三码区别大吗| 18在线观看网站| bbb黄色大片| 高清av免费在线| 亚洲色图综合在线观看| 亚洲自偷自拍图片 自拍| 亚洲精品一卡2卡三卡4卡5卡| 自线自在国产av| 欧美午夜高清在线| 一个人免费在线观看的高清视频| 亚洲av日韩在线播放| 国产av又大| 成人三级做爰电影| 亚洲av成人一区二区三| 久久国产精品影院| 成人av一区二区三区在线看| 999久久久国产精品视频| 亚洲少妇的诱惑av| 两个人免费观看高清视频| kizo精华| 人人妻人人爽人人添夜夜欢视频| 亚洲av日韩精品久久久久久密| 最近最新中文字幕大全免费视频| 桃红色精品国产亚洲av| 人人妻人人澡人人看| 夫妻午夜视频| 99精品在免费线老司机午夜| 下体分泌物呈黄色| 久久久久久久国产电影| 亚洲美女黄片视频| 高清视频免费观看一区二区| 免费女性裸体啪啪无遮挡网站| 十八禁人妻一区二区| av片东京热男人的天堂| 国产成人啪精品午夜网站| 亚洲欧美日韩高清在线视频 | 无遮挡黄片免费观看| 热99久久久久精品小说推荐| 成人亚洲精品一区在线观看| 狂野欧美激情性xxxx| 久久人妻av系列| 成人永久免费在线观看视频 | 免费女性裸体啪啪无遮挡网站| 久久久久久久久久久久大奶| 日本av免费视频播放| 久久久久久亚洲精品国产蜜桃av| 日本av手机在线免费观看| 久久ye,这里只有精品| 亚洲欧美激情在线| 精品亚洲成国产av| 欧美午夜高清在线| 丝袜美腿诱惑在线| 亚洲人成伊人成综合网2020| 两个人看的免费小视频| 欧美乱码精品一区二区三区| 成年人午夜在线观看视频| 久久久久久免费高清国产稀缺| 人人妻,人人澡人人爽秒播| 久久天堂一区二区三区四区| 国产男女超爽视频在线观看| 国产精品秋霞免费鲁丝片| 色综合婷婷激情| 日本五十路高清| 中文字幕另类日韩欧美亚洲嫩草| 在线av久久热| 免费日韩欧美在线观看| 成年动漫av网址| 黄频高清免费视频| 男女免费视频国产| 国产成人精品无人区| 久久亚洲真实| 久久精品91无色码中文字幕| 久久午夜亚洲精品久久| 香蕉国产在线看| tocl精华| 国产单亲对白刺激| 制服人妻中文乱码| 久久青草综合色| 精品国产国语对白av| 免费日韩欧美在线观看| 999久久久国产精品视频| 美女午夜性视频免费| 亚洲久久久国产精品| 色精品久久人妻99蜜桃| 国产成人免费观看mmmm| 两个人看的免费小视频| 91九色精品人成在线观看| 成人特级黄色片久久久久久久 | 免费在线观看完整版高清| 国产伦人伦偷精品视频| 国产精品偷伦视频观看了| 国产成人欧美在线观看 | 黑人猛操日本美女一级片| 欧美久久黑人一区二区| 美国免费a级毛片| 国产精品亚洲一级av第二区| 精品视频人人做人人爽| 亚洲精品国产精品久久久不卡| 欧美精品人与动牲交sv欧美| 一区二区av电影网| 91大片在线观看| 女人爽到高潮嗷嗷叫在线视频| 国产精品1区2区在线观看. | 9191精品国产免费久久| 国产亚洲av高清不卡| 黑人巨大精品欧美一区二区蜜桃| h视频一区二区三区| 久久av网站| 91av网站免费观看| 国产成人精品无人区| 亚洲一卡2卡3卡4卡5卡精品中文| 黄色片一级片一级黄色片| 国产一区二区 视频在线| 色婷婷久久久亚洲欧美| 国产在线观看jvid| 亚洲黑人精品在线| 老司机福利观看| 成年版毛片免费区| 精品视频人人做人人爽| 国产精品 国内视频| 2018国产大陆天天弄谢| 欧美日韩精品网址| 亚洲视频免费观看视频| 亚洲人成电影观看| 天堂8中文在线网| 日韩中文字幕欧美一区二区| 亚洲国产看品久久| 91成年电影在线观看| 欧美日韩视频精品一区| 欧美中文综合在线视频| av电影中文网址| 伊人久久大香线蕉亚洲五| 天堂动漫精品| 久久久久久久久久久久大奶| 在线av久久热| 亚洲午夜理论影院| 在线亚洲精品国产二区图片欧美| 日韩欧美一区二区三区在线观看 | 免费黄频网站在线观看国产| 国产男女超爽视频在线观看| e午夜精品久久久久久久| 高清欧美精品videossex| 啦啦啦免费观看视频1| 人成视频在线观看免费观看| www.精华液| 超碰97精品在线观看| 十八禁网站网址无遮挡| 亚洲免费av在线视频| 亚洲一区中文字幕在线| 欧美午夜高清在线| 亚洲一卡2卡3卡4卡5卡精品中文| 91大片在线观看| 国产精品国产高清国产av | avwww免费| 啪啪无遮挡十八禁网站| 亚洲国产中文字幕在线视频| 国产精品亚洲一级av第二区| 亚洲成人国产一区在线观看| 国产av又大| 亚洲色图综合在线观看| 嫩草影视91久久| 叶爱在线成人免费视频播放| 欧美另类亚洲清纯唯美| 少妇 在线观看| 99精品欧美一区二区三区四区| 男人操女人黄网站| 日日夜夜操网爽| av天堂在线播放| 一夜夜www| 成人影院久久| 老司机靠b影院| 大陆偷拍与自拍| 2018国产大陆天天弄谢| 中文亚洲av片在线观看爽 | 国产在线免费精品| 狂野欧美激情性xxxx| 一级黄色大片毛片| 精品久久久久久久毛片微露脸| 国产一区二区 视频在线| 日韩视频在线欧美| 久久久久视频综合| 国产野战对白在线观看| 中文字幕人妻熟女乱码| 黄色 视频免费看| 午夜福利在线观看吧| 超碰成人久久| 99国产综合亚洲精品| 精品一区二区三区av网在线观看 | 女人精品久久久久毛片| 国产深夜福利视频在线观看| 中文欧美无线码| 欧美激情久久久久久爽电影 | 成年版毛片免费区| 狠狠狠狠99中文字幕| 午夜福利视频精品| 国产一区二区激情短视频| 久久久久精品国产欧美久久久| 亚洲欧美色中文字幕在线| 三级毛片av免费| 久久久久久免费高清国产稀缺| 夜夜爽天天搞| 老熟妇乱子伦视频在线观看| 极品教师在线免费播放| 男女高潮啪啪啪动态图| 亚洲 国产 在线| 69av精品久久久久久 | 精品视频人人做人人爽| 欧美亚洲 丝袜 人妻 在线| 丁香六月欧美| 国产精品麻豆人妻色哟哟久久| 汤姆久久久久久久影院中文字幕| 怎么达到女性高潮| 国产精品久久久av美女十八| 色精品久久人妻99蜜桃| 久久国产亚洲av麻豆专区| 亚洲欧美色中文字幕在线| 汤姆久久久久久久影院中文字幕| 成人精品一区二区免费| 国产又色又爽无遮挡免费看| 香蕉国产在线看| 日韩免费高清中文字幕av| 18禁国产床啪视频网站| 三级毛片av免费| 久久久久网色| 亚洲成av片中文字幕在线观看| 性少妇av在线| a级毛片黄视频| 90打野战视频偷拍视频| 欧美另类亚洲清纯唯美| 欧美日韩亚洲高清精品| 99国产精品一区二区蜜桃av | av不卡在线播放| 99re6热这里在线精品视频| 狂野欧美激情性xxxx| 建设人人有责人人尽责人人享有的| 精品久久久久久久毛片微露脸| 可以免费在线观看a视频的电影网站| 国产精品亚洲一级av第二区| 色94色欧美一区二区| 法律面前人人平等表现在哪些方面| 国产精品久久久久久精品古装| 国内毛片毛片毛片毛片毛片| 免费av中文字幕在线| 97在线人人人人妻| 汤姆久久久久久久影院中文字幕| 色综合婷婷激情| 久久毛片免费看一区二区三区| 极品人妻少妇av视频| 国产精品98久久久久久宅男小说| 欧美激情极品国产一区二区三区| 国产亚洲欧美精品永久| 老汉色∧v一级毛片| 69av精品久久久久久 | 成年人免费黄色播放视频| av网站在线播放免费| 成人av一区二区三区在线看| 18禁美女被吸乳视频| 国精品久久久久久国模美| 国产在线一区二区三区精| 国产成人影院久久av| 香蕉久久夜色| 高清视频免费观看一区二区| 亚洲国产欧美在线一区| 久久久久久久大尺度免费视频| 亚洲av欧美aⅴ国产| 黑人操中国人逼视频| 久久av网站| 大片电影免费在线观看免费| 国产在线一区二区三区精| 久久午夜亚洲精品久久| 亚洲一区中文字幕在线| 午夜福利免费观看在线| 91老司机精品| 女人高潮潮喷娇喘18禁视频| 久久久久国内视频| 亚洲欧美一区二区三区久久| 中文字幕人妻丝袜一区二区| 婷婷丁香在线五月| 女人高潮潮喷娇喘18禁视频| 欧美+亚洲+日韩+国产| 亚洲成av片中文字幕在线观看| 首页视频小说图片口味搜索| 超碰97精品在线观看| 捣出白浆h1v1| 人妻久久中文字幕网| 午夜两性在线视频| 亚洲九九香蕉| 国产精品一区二区在线不卡| 国产野战对白在线观看| 免费人妻精品一区二区三区视频| 美国免费a级毛片| 免费在线观看影片大全网站| 久久免费观看电影| 大片免费播放器 马上看| 考比视频在线观看| 夜夜夜夜夜久久久久| 久久久精品区二区三区| 亚洲欧美一区二区三区黑人| 精品国内亚洲2022精品成人 | 日韩大片免费观看网站| 免费在线观看影片大全网站| 69av精品久久久久久 | 免费高清在线观看日韩| 久久人人爽av亚洲精品天堂| 久久精品熟女亚洲av麻豆精品| bbb黄色大片| 黄网站色视频无遮挡免费观看| 亚洲精华国产精华精| 99九九在线精品视频| 午夜福利,免费看| 亚洲国产中文字幕在线视频| 国产野战对白在线观看| 99精品久久久久人妻精品| 一级a爱视频在线免费观看| 欧美激情高清一区二区三区| 亚洲色图综合在线观看| 蜜桃在线观看..| 成人18禁在线播放| 在线观看人妻少妇| 亚洲天堂av无毛| 麻豆成人av在线观看| 久久精品国产亚洲av高清一级| 午夜福利视频精品| 久久av网站| 色尼玛亚洲综合影院| 国产精品影院久久| 免费观看a级毛片全部| 日韩一区二区三区影片| 肉色欧美久久久久久久蜜桃| 王馨瑶露胸无遮挡在线观看| 一二三四社区在线视频社区8| 在线播放国产精品三级| 美女午夜性视频免费| 国产成人免费观看mmmm| 成人国产一区最新在线观看| 一区在线观看完整版| 国产日韩欧美视频二区| 亚洲黑人精品在线| 又黄又粗又硬又大视频| 90打野战视频偷拍视频| 久久青草综合色| 免费在线观看完整版高清| 夫妻午夜视频| 在线观看人妻少妇| 亚洲伊人久久精品综合| 午夜精品国产一区二区电影| 如日韩欧美国产精品一区二区三区| 久热爱精品视频在线9| 亚洲av日韩在线播放| 一本久久精品| 一个人免费在线观看的高清视频| xxxhd国产人妻xxx| 国产精品免费一区二区三区在线 | 日本撒尿小便嘘嘘汇集6| 亚洲国产欧美日韩在线播放| 亚洲人成电影观看| 亚洲va日本ⅴa欧美va伊人久久| 极品教师在线免费播放| 美女福利国产在线| 大香蕉久久网| 窝窝影院91人妻| 一本一本久久a久久精品综合妖精| 男女高潮啪啪啪动态图| 日韩 欧美 亚洲 中文字幕| 在线av久久热| 男女午夜视频在线观看| 最黄视频免费看| 婷婷成人精品国产| 一区在线观看完整版| 美女视频免费永久观看网站| 大香蕉久久成人网| 日韩中文字幕欧美一区二区| 亚洲人成电影观看| 精品少妇一区二区三区视频日本电影| 亚洲av欧美aⅴ国产| 女性被躁到高潮视频| 高清黄色对白视频在线免费看| 中文字幕精品免费在线观看视频| 大陆偷拍与自拍| 国产亚洲精品一区二区www | 一边摸一边抽搐一进一小说 | 在线亚洲精品国产二区图片欧美| a在线观看视频网站| 女同久久另类99精品国产91| 黄色片一级片一级黄色片| 丁香欧美五月| 久久久久久久久免费视频了| 人人妻人人澡人人爽人人夜夜| 色婷婷久久久亚洲欧美| 久久九九热精品免费| 99国产极品粉嫩在线观看| 美女扒开内裤让男人捅视频| 捣出白浆h1v1| 亚洲va日本ⅴa欧美va伊人久久| 99riav亚洲国产免费| 天天躁夜夜躁狠狠躁躁| 欧美日韩亚洲国产一区二区在线观看 | 男人舔女人的私密视频| 国产亚洲欧美精品永久| 国产欧美日韩精品亚洲av| 精品国产国语对白av| 久久亚洲真实| 午夜精品国产一区二区电影| 大片电影免费在线观看免费| 日本精品一区二区三区蜜桃| 亚洲专区字幕在线| 久久精品亚洲av国产电影网| 日本五十路高清| 欧美国产精品一级二级三级| 国产精品九九99| www日本在线高清视频| 高清毛片免费观看视频网站 | 香蕉久久夜色| a在线观看视频网站| av线在线观看网站| 黄色丝袜av网址大全| 欧美日韩中文字幕国产精品一区二区三区 | 午夜福利视频在线观看免费| 高清在线国产一区| 露出奶头的视频| 国产精品一区二区精品视频观看| 日本一区二区免费在线视频| 久久亚洲精品不卡| 啦啦啦在线免费观看视频4| 久久精品国产亚洲av高清一级| 久久久久久免费高清国产稀缺| 国产区一区二久久| 蜜桃在线观看..| 9191精品国产免费久久| 一边摸一边抽搐一进一小说 | 极品教师在线免费播放| 亚洲欧美激情在线| 国产区一区二久久| 国产伦理片在线播放av一区| 日韩免费高清中文字幕av| 国产成人影院久久av| 精品人妻1区二区| 亚洲欧洲精品一区二区精品久久久| 90打野战视频偷拍视频| 麻豆成人av在线观看| 精品久久久久久电影网| kizo精华| 一本综合久久免费| 中文欧美无线码| 亚洲欧美色中文字幕在线| 我的亚洲天堂| 亚洲欧美一区二区三区黑人| 国产成人欧美| 精品午夜福利视频在线观看一区 | 久久ye,这里只有精品| 天天添夜夜摸| 精品国产亚洲在线| 日韩中文字幕视频在线看片| 国产一区二区三区综合在线观看| 激情在线观看视频在线高清 | 美女福利国产在线| 1024视频免费在线观看| 99国产极品粉嫩在线观看| 久久ye,这里只有精品| 久久中文字幕人妻熟女| 动漫黄色视频在线观看| 色婷婷av一区二区三区视频| 午夜福利,免费看| 啦啦啦视频在线资源免费观看| 国产av精品麻豆| 国产精品久久久久久精品电影小说| 麻豆av在线久日| 精品国产超薄肉色丝袜足j| 亚洲五月婷婷丁香| 国产欧美日韩一区二区三区在线| 夜夜夜夜夜久久久久| 啦啦啦免费观看视频1| 成人18禁在线播放| 精品少妇黑人巨大在线播放| 亚洲精品av麻豆狂野| 大陆偷拍与自拍| 午夜激情av网站| 大型av网站在线播放| 757午夜福利合集在线观看| 肉色欧美久久久久久久蜜桃| 老司机在亚洲福利影院| 国产熟女午夜一区二区三区| 涩涩av久久男人的天堂| 18禁黄网站禁片午夜丰满| 欧美亚洲日本最大视频资源| 国产精品美女特级片免费视频播放器 | 91大片在线观看| 国产一区二区三区综合在线观看| 热re99久久国产66热| 久久精品国产亚洲av高清一级| 在线观看66精品国产| 亚洲熟女精品中文字幕| 97在线人人人人妻| 2018国产大陆天天弄谢| 亚洲av美国av| 91精品国产国语对白视频| 国产亚洲精品久久久久5区| 亚洲精品一卡2卡三卡4卡5卡| 国产av又大| 色在线成人网| 亚洲欧洲精品一区二区精品久久久| 亚洲精品美女久久av网站| 女人爽到高潮嗷嗷叫在线视频| 欧美在线一区亚洲| 热re99久久精品国产66热6| 99久久精品国产亚洲精品| 久久精品熟女亚洲av麻豆精品| 亚洲情色 制服丝袜| 午夜福利免费观看在线| 一边摸一边抽搐一进一出视频| 这个男人来自地球电影免费观看| 女人被躁到高潮嗷嗷叫费观| 亚洲免费av在线视频| 欧美变态另类bdsm刘玥| 久久99热这里只频精品6学生| 99国产综合亚洲精品| 九色亚洲精品在线播放| 一区在线观看完整版| 精品亚洲成国产av| 亚洲精品国产色婷婷电影| 亚洲中文av在线| 国产精品久久久久成人av| 汤姆久久久久久久影院中文字幕| 大型黄色视频在线免费观看| 婷婷丁香在线五月| 亚洲成av片中文字幕在线观看| 国产在视频线精品|