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

    鄂爾多斯地區(qū)上地幔巖石圈三維速度結(jié)構(gòu)面波反演研究

    2012-12-18 05:29:06周仕勇陳永順馮永革
    地球物理學(xué)報 2012年5期
    關(guān)鍵詞:平面波克拉通面波

    李 多,周仕勇,陳永順,馮永革,李 鵬

    北京大學(xué)地球與空間科學(xué)學(xué)院,北京 100871

    鄂爾多斯地區(qū)上地幔巖石圈三維速度結(jié)構(gòu)面波反演研究

    李 多,周仕勇*,陳永順,馮永革,李 鵬

    北京大學(xué)地球與空間科學(xué)學(xué)院,北京 100871

    雙平面波擬合法是一種新的面波成像方法,反演中考慮地震波場中的非平面波成分,提高反演的分辨率.本文利用雙平面波擬合法,反演獲得鄂爾多斯地區(qū)上地幔巖石圈的速度結(jié)構(gòu).所用資料為國家數(shù)字地震臺網(wǎng)69個寬頻帶地震儀和北京大學(xué)34個流動數(shù)字地震臺觀測到的地震波面波資料.首先從面波記錄中提取了研究區(qū)域20~125s瑞利波相速度頻散曲線,進(jìn)而得到各個周期瑞利波相速度異常分布圖.結(jié)果顯示,短周期瑞利波相速度異常與地表的構(gòu)造特征吻合較好,中長周期的瑞利波相速度可以反映出上地幔巖石圈的速度異常分布以及構(gòu)造特征.由研究區(qū)20~125s的瑞利波相速度分布圖可以反演得到地表到地下200km范圍內(nèi)的三維剪切波速度結(jié)構(gòu).結(jié)果顯示,鄂爾多斯塊體內(nèi)部穩(wěn)定均一,活化或改造的痕跡不明顯;鄂爾多斯塊體西南緣受到青藏高原的強(qiáng)烈作用,有大量地幔物質(zhì)流動的痕跡存在;中央轉(zhuǎn)換帶下超過200km深度存在地幔物質(zhì)上涌,可能與太平洋板塊的俯沖和青藏高原板塊的擠壓有關(guān).

    雙平面波擬合法,面波成像,鄂爾多斯,華北克拉通,巖石圈

    1 引 言

    鄂爾多斯塊體屬于華北克拉通的一部分.華北克拉通可分為東部塊體、西部塊體和中間轉(zhuǎn)換帶三部分,其中西部塊體包括陰山塊體和鄂爾多斯塊體.地質(zhì)學(xué)研究表明,華北克拉通的東、西塊體在早元古代至中元古代相互分離,各自經(jīng)歷了不同的發(fā)育歷史.根據(jù)鄂爾多斯塊體北緣的孔慈巖證據(jù)顯示,鄂爾多斯塊體在約2.0~1.9Ga與北邊的陰山板塊碰撞拼合,在約1.85Ga與東部板塊碰撞拼合,并形成了現(xiàn)在的中央轉(zhuǎn)換帶[1-3].晚中生代到新生代時期,伴隨著巖石圈減薄、斷陷盆地形成、地下熱流、火山活動等,發(fā)生了華北克拉通的活化[4-7].

    最近十幾年來,中國東部巖石圈的減薄和華北克拉通的活化成了國內(nèi)外學(xué)者的熱門研究問題,相關(guān)的研究成果層出不窮.地震波層析成像是用于研究地球內(nèi)部結(jié)構(gòu)的最有效手段,在地球深部結(jié)構(gòu)探測研究中具有重要作用.已有很多地球物理學(xué)家利用寬頻帶地震臺網(wǎng)記錄,運(yùn)用各種地震波層析成像技術(shù)來探測鄂爾多斯地區(qū)的地殼和上地幔結(jié)構(gòu).已有的瑞利面波成像和P波成像得到了華北地區(qū)的三維地殼和上地幔結(jié)構(gòu),結(jié)果顯示鄂爾多斯塊體整體堅固穩(wěn)定[8-10].接收函數(shù)的結(jié)果顯示出華北克拉通西部板塊和東部板塊的殼-幔過渡帶存在明顯的結(jié)構(gòu)差異[11-12].橫波分裂的結(jié)果顯示山西斷陷帶上地幔的各向異性方向與地表斷陷帶的走向一致,這與拉張區(qū)地幔物質(zhì)流動有關(guān)[13].利用P波和S波層析成像技術(shù)得到的華北克拉通速度異常分布結(jié)果顯示,中央轉(zhuǎn)換帶下超過500km深的低速異常是由于地幔物質(zhì)上涌引起的,這與華北克拉通東部板塊的活化密切相關(guān)[14].

    利用面波相速度的頻散特征反演地下速度結(jié)構(gòu),是地殼及上地幔結(jié)構(gòu)探測的主要地震學(xué)方法之一.目前,單臺法和雙臺法仍然是國內(nèi)地震學(xué)家提取面波頻散特性的主要方法[15].單臺法測量的是源臺間面波的頻散特性,在較近的臺站上計算得到的面波相速度會存在較大的誤差,因此不可能用于區(qū)域小尺度速度結(jié)構(gòu)精細(xì)反演.而雙臺法則需要兩個接收臺站和震源位于同一大圓路徑上,利用滿足此條件的雙臺地震資料可以求得雙臺間的面波頻散特性,然而兩個接收臺站和震源嚴(yán)格位于同一大圓路徑上的記錄并不很多,使雙臺法的應(yīng)用同樣受到局限.

    傳統(tǒng)的單臺法及雙臺法提取面波頻散特性的理論基礎(chǔ)也面臨挑戰(zhàn).Wielandt和Friederich自1993年來[16-19]發(fā)表的系列理論研究表明單臺法及雙臺法提取面波頻散特性的理論基礎(chǔ)所隱含的平面波假設(shè)在實(shí)際中很難成立.Wielandt[19]研究指出在地震面波的非平面波成分不可忽視的時候,傳統(tǒng)面波相速度測量方法所得到的只是波場的動態(tài)相速度.動態(tài)相速度除了包含傳播介質(zhì)的物理性質(zhì)信息之外,還包含了面波波場振幅局部不均勻分布的信息;因此不能直接反映研究區(qū)域介質(zhì)的物理性質(zhì).只有正確估計面波波場局部的結(jié)構(gòu),才可能得到真正與介質(zhì)相關(guān)的面波結(jié)構(gòu)相速度.為解決這一問題我們擬應(yīng)用Forsyth和Li[20]提出的一種通過用雙平面面波擬合地震臺陣記錄的非均勻面波波場,提取臺陣覆蓋區(qū)介質(zhì)的面波相速度頻散曲線的方法,這樣因?yàn)闊o平面波假設(shè),減少了“結(jié)構(gòu)性誤差”,理論更為嚴(yán)謹(jǐn).而且面波波場擬合法同時利用了臺站記錄的面波振幅和相位信息,且不受大圓路徑的局限,利用的資料更豐富.

    Jiang等[21]曾經(jīng)利用雙平面面波擬合法開展了藏南地區(qū)面波精細(xì)成像并發(fā)現(xiàn)藏南底下存在約40km厚的軟弱下地殼層,為青藏高原下地殼流模型提供了支持.本文作者李鵬曾經(jīng)利用雙平面面波擬合法開展了鄂爾多斯塊體和中央轉(zhuǎn)換帶地區(qū)地下結(jié)構(gòu)的面波反演研究[22],探測到穩(wěn)定的鄂爾多斯塊體巖石圈可能延伸至底下120~140km深度,而山西斷陷區(qū)及華北盆地的巖石圈可能減薄至70~80km深度,給出了華北克拉通現(xiàn)代處于活化過程中及其可能動力學(xué)圖像的面波探測證據(jù),同時也驗(yàn)證了雙平面面波擬合法的有效性.但因?yàn)楫?dāng)時沒有收集到銀川、烏海等地鄂爾多斯西北緣的幾個中國國家臺網(wǎng)臺站的記錄(圖1),不能對西緣—青藏板塊與鄂爾多斯塊體交匯區(qū)底下結(jié)構(gòu)有清晰成像,并影響研究區(qū)成像的整體分辨率.本研究我們將在李鵬等[22]的工作的基礎(chǔ)上,補(bǔ)充銀川臺等23個鄂爾多斯西北緣的中國國家臺網(wǎng)臺站的記錄,用雙平面波擬合法反演開展鄂爾多斯塊體和中央轉(zhuǎn)換帶地區(qū)瑞利波相速度二維分布圖像的反演,期望得到研究區(qū)更高分辨率的面波成像結(jié)果,為認(rèn)識華北克拉通的形成和演化提供更可靠的地震學(xué)證據(jù).

    圖1 所用臺站分布圖及研究區(qū)域構(gòu)造地質(zhì)構(gòu)造圖.紅色三角為國家地震局固定臺站,綠色三角為北京大學(xué)流動臺站.白色方框?yàn)檠芯繀^(qū)域.臺站主要分布在山西斷陷帶和渭河斷陷帶附近Fig.1 Distribution of the stations and tomography and tectonics of the study region.Red triangles denote the stations of CENC;Green triangles denote the stations of PKU.White pane denotes the study region.Stations distribute densely along Shanxi Rift and Weihe Rift

    2 數(shù)據(jù)處理

    本文的數(shù)據(jù)來源為兩部分:中國地震局臺網(wǎng)中心在33°N—42°N和105°E—115°E的區(qū)域范圍內(nèi)的69個固定臺站的地震波形數(shù)據(jù)[23].記錄時間為2007年8月—2009年11月;北京大學(xué)在山西斷陷帶布設(shè)的34個二維流動數(shù)字臺網(wǎng)的地震波形數(shù)據(jù),記錄時間為2007年1—6月.在數(shù)據(jù)選取時,只選取了瑞利面波的垂直向記錄,同時根據(jù)信噪比,篩選了196個震級≥6.0級、震中距為30°~120°、方位角分布較為均勻的遠(yuǎn)震事件(圖2).

    處理數(shù)據(jù)過程中,將所有臺站的波形信號統(tǒng)一校正為CMG-3ESPC的拾震器的儀器響應(yīng).選取帶寬為10mHz的窄帶濾波器(巴特沃斯型濾波器,級數(shù)4,通道數(shù)2).從0.008Hz到0.05Hz選取了12個中心頻率點(diǎn),對地震波進(jìn)行濾波,從而獲取了12個窄頻段的波形數(shù)據(jù).對于每個地震事件,反演方程中包含6個波場參數(shù)(每個平面波需要用振幅、初相位與傳播方位角3個參數(shù)描述)和2個格點(diǎn)的相速度參數(shù),因此,我們保留射線記錄數(shù)大于9的地震事件,共得到約27000條地震射線,各頻段(周期)的記錄數(shù)條形統(tǒng)計圖如圖3.圖4顯示了50s周期的地震射線分布,可以看出射線覆蓋較好的地區(qū)為南部和東部,北部和西部由于受到臺站分布和地震事件分布的局限,射線分布較為稀疏.

    圖4 周期為50s的地震射線覆蓋分布圖.研究區(qū)域東部和南部射線覆蓋密集Fig.4 Great circle ray path coverage at the period of 50s.The density on east and south is good

    3 方法概述

    3.1 動態(tài)波數(shù)與結(jié)構(gòu)波數(shù)的區(qū)別

    考慮一列沿地球表面?zhèn)鞑サ牡卣鸩ǎ娌ǎ?,其子波頻率為ω的波場函數(shù)可以寫為

    (x,y)地表上某觀測點(diǎn)坐標(biāo),波數(shù)矢量Ko=(kx,ky),ko=|Ko|,令

    則有

    由此可得到(x,y)點(diǎn)上頻率為ω的相速度:

    注意到(3),(4)式中的波數(shù)和相速度是直接由測量面波波場在空間的分布而得到,我們稱之為波場動態(tài)波數(shù)|相速度(Dynamic wave wavenumber|phase velocity).

    由動力學(xué)方程出發(fā),我們可以推導(dǎo)出(1)式表達(dá)的波場函數(shù)滿足Helmholtz方程,即有

    (5)式中k=ω/c.由于(5)式直接由動力學(xué)方程推導(dǎo),k或c由波的傳播介質(zhì)的彈性參數(shù)唯一確定,我們將k,c稱為結(jié)構(gòu)波數(shù)和結(jié)構(gòu)相速度(structural wavenumber|phase velocity).由(5)式有

    顯然(6)式與(3)式所表達(dá)的波數(shù)在物理上并不完全一樣,數(shù)學(xué)表達(dá)上也存在差別.

    由于波場函數(shù)必須滿足Helmholtz方程,我們將(2)式代入(6)式有

    將(3)式代入(7)則有

    由(8)式可清楚地看出:當(dāng)面波不是平面波時,即臺陣記錄的同一面波因記錄臺位置不同,其振幅不一樣,(1)式中的A亦即(8)式中a在空間上不是常量,則觀測的波場動態(tài)波數(shù)k0與結(jié)構(gòu)波數(shù)k是存在系統(tǒng)偏差的.傳統(tǒng)的臺源法和雙臺法測量的面波相速度頻散特性只是觀測的面波波場動態(tài)相速度頻散特性,是一種平面波假設(shè)下的近似測量.

    3.2 雙平面波擬合法

    為得到精確的面波反演結(jié)果,本研究擬采用Forsyth與Li[20]方法原理,其基本原理概述如下:

    對頻率為ω沿水平方向傳播的非平面波可以用兩組同頻率、不同水平方向傳播、具有不同相位和振幅的平面波的干涉效應(yīng)表達(dá),即非平面波垂直向位移波場可以表達(dá)為[20]

    我們將研究區(qū)分成N個等尺度網(wǎng)格,一般情形下,介質(zhì)速度具有各向異性,因此不同地震的面波通過j網(wǎng)格,由于射線穿過的方位不同,速度不同.因此i地震面波在j網(wǎng)絡(luò)傳播的速度為

    我們將研究區(qū)域劃分為0.5°×0.5°的網(wǎng)格,忽略每個格點(diǎn)的各向異性成分,即令B1,B2為0,每個格點(diǎn)的瑞利波相速度為均一值.對每個格點(diǎn)加入高斯平滑因子,則每個格點(diǎn)的慢度表示為

    圖5 研究區(qū)域建立的區(qū)域直角坐標(biāo)系.參考臺站位于坐標(biāo)原點(diǎn),(r,ψ)是與直角坐標(biāo)系(x,y)相對應(yīng)的極坐標(biāo)系坐標(biāo),橫軸為射線大圓弧路徑方向,平面波以與橫軸夾角為θ的方向入射.虛線方框?yàn)檠芯繀^(qū)域,相對到時零點(diǎn)(τ=0)為沿大圓弧路徑傳播的面波波前最先接觸到研究區(qū)域的時刻[24]Fig.5 Local coordinate system of two-plane-wave method.Reference station is at(0,0),(r,ψ)is the position of polar coordinate as(x,y)in rectangular coordinate.x-coordinate is along the great circle path,while y axis is perpendicular to xaxis.The incoming plane wave is along the direction of propagationθ.The dashed square frame denotes the study region.The time that the wavefront along the great circle path arrives atthe first point of the study region is 0(τ=0)[24]

    其中,i表示地震編號,j表示格點(diǎn)編號,高斯權(quán)重函數(shù)為

    其中,特征長度Lw直接影響高斯權(quán)重函數(shù)的形狀,可以控制反演得到的速度結(jié)構(gòu)的光滑程度,Lw過小會導(dǎo)致空間小尺度內(nèi)出現(xiàn)突變;過大則會使得合理的速度差異不夠顯著.所以,在反演過程中需要調(diào)整特征長度來調(diào)節(jié)結(jié)果的分辨率和標(biāo)準(zhǔn)差范圍[20].經(jīng)過不同的嘗試,以獲取研究區(qū)最佳分辨率為標(biāo)準(zhǔn),我們將Lw取為50km.

    將6個波場參數(shù)和每個格點(diǎn)的慢度參數(shù)帶入反演,分別利用模擬退火法[25]和廣義最小二乘法反演[26],可以得到每個格點(diǎn)的瑞利波相速度.

    4 反演結(jié)果與分析

    4.1 檢測板試驗(yàn)

    檢測板是驗(yàn)證結(jié)果的分辨率和可靠性的重要依據(jù)[27].在本文的研究中,我們將輸入的速度模型的異常范圍設(shè)定為±5%之間,異常尺度設(shè)置為2°×2°,即空間尺度約為200km×200km,這一尺度基本可以分辨出研究區(qū)域的構(gòu)造單元的速度異常.

    通過將反演結(jié)果與輸入模型進(jìn)行比對(圖6),可以看出,短周期時,僅有山西斷陷帶地區(qū)臺站密集處可以得到反演結(jié)果,這是因?yàn)楦哳l信號衰減快、信噪比低.中長周期時,鄂爾多斯塊體、山西斷陷帶和南部的渭河斷陷帶處結(jié)果較為清晰準(zhǔn)確,而北部,尤其是西北部,由于射線分布稀疏,反演可信度不高.由于本文的目的就是研究鄂爾多斯塊體和山西斷陷帶的上地幔巖石圈結(jié)果,西北部和東北部的低分辨檢測結(jié)果對我們后面的研究結(jié)論影響不大.

    4.2 瑞利波相速度

    首先將研究區(qū)域內(nèi)各點(diǎn)瑞利波相速度設(shè)為相同,反演得到一維瑞利波相速度頻散曲線(圖7).在短周期時,該地區(qū)相速度基本吻合AK135模型[28],而在中長周期時,該地區(qū)相速度低于AK135模型[28],反映了華北克拉通的活化而引起的巖石圈整體的溫度升高.

    以一維平均瑞利波速度為初始模型,反演得到二維瑞利波相速度異常分布圖(圖8).圖中紅色為低速異常,藍(lán)色為高速異常.短周期時(20~29s)的瑞利波反映地表到莫霍面的速度異常分布,與地表的構(gòu)造特征基本吻合,如太原盆地為新生代的斷陷盆地,沉積層較厚,表現(xiàn)為瑞利波的低速異常,而太行山、呂梁山等為高速異常.

    中長周期(33~125s)反映了上地幔巖石圈的速度結(jié)構(gòu).這里可以看到鄂爾多斯塊體整體為明顯的高速異常體,內(nèi)部速度差異較小,與這里塊體內(nèi)部少地震的事實(shí)相吻合.以山西斷陷帶為界,東西兩側(cè)的高低速異常分割明顯.低速異常沿中央轉(zhuǎn)換帶呈南北條帶分布.

    鄂爾多斯塊體是中國大陸最古老的克拉通,自形成以來,沒有經(jīng)歷過活化和改造,被認(rèn)為是最堅固穩(wěn)定的大陸巖石圈[29].晚太古代至元古代時期,華北克拉通的西部板塊與東部板塊碰撞拼合,拼合部位在中央轉(zhuǎn)換帶[30].晚中生代到新生代,華北克拉通再次活化,華北地區(qū)巖石圈拉張減薄,形成了一系列的斷陷盆地[31-32].巖石圈的熱結(jié)構(gòu)研究[33]顯示,山西斷陷帶存在較高的熱值,對應(yīng)這里的低速異常,很有可能與地幔物質(zhì)流動有關(guān).

    4.3 剪切波速度結(jié)構(gòu)

    圖6 二維瑞利波相速度檢測板試驗(yàn).輸入速度異常模型為-5%~5%之間Fig.6 Checkerboard test for the 2-D Rayleigh wave isotropic phase velocity.The input velocity model is with anomalies within-5%~5%

    剪切波速度相對于物質(zhì)的流變性變化敏感.不同深度剪切波速度對于不同頻率的瑞利面波相速度的敏感程度不同,大體上為面波波長的1/3處為峰值,向兩側(cè)逐漸遞減[34].并且越深處的剪切波對各頻率相速度的敏感程度差別減弱.本文所用研究方法是Weerarantne等[35]和Li等[27]在Saito[36]的方法基礎(chǔ)上而發(fā)展起來的.參考模型選擇了AK135模型[28],初始地殼模型為三層(圖9).反演結(jié)果受到剪切波速度阻尼系數(shù)和地殼厚度阻尼系數(shù)的共同控制.阻尼系數(shù)過大會導(dǎo)致誤差較大,無法得到最優(yōu)化解;阻尼系數(shù)過小則會出現(xiàn)約束不夠強(qiáng),導(dǎo)致解的不穩(wěn)定.所以阻尼系數(shù)的選擇要權(quán)衡兩方面的因素,反演過程中要根據(jù)結(jié)果調(diào)整參數(shù),以得到較好的結(jié)果.本文中根據(jù)多次反演結(jié)果的比較,設(shè)定剪切波速度阻尼系數(shù)和地殼厚度阻尼系數(shù)為0.05和2.另外,由于莫霍面深度對于反演結(jié)果影響較大,所以反演過程不限定莫霍面深度,這樣保證了結(jié)果的穩(wěn)定性和可靠性[37].

    圖7 一維瑞利波相速度頻散曲線.實(shí)線為AK135參考模型,長虛線為輸入模型,短虛線黑線為反演結(jié)果Fig.7 Average Rayleigh-wave phase velocities in Ordos region(black circles with error bar)compared with those of global AK135model(solid black line).The long dashed line is the input velocities for 12different periods

    圖8 二維瑞利波相速度異常分布圖Fig.8 Maps of Rayleigh wave isotropic phase velocity anomalies

    一維平均剪切波速度曲線圖(圖9)顯示,該地區(qū)整體上地幔巖石圈速度低于AK135全球平均模型[28].須指出的是,圖7顯示鄂爾多斯的相速度頻散曲線在幾乎整個地殼(50s以內(nèi))內(nèi)與AK135非常一致,但是圖9顯示的反演得到的一維剪切波速度結(jié)構(gòu)則顯示鄂爾多斯內(nèi)地殼速度明顯低于AK135,這盡管可能與我們觀測的頻散曲線缺乏高頻段信息(短周期截止點(diǎn)在20s),對地殼速度的約束較弱有關(guān),但長周期(大于50s)部分低相速度,也可能造成對地殼速度的反演結(jié)果系統(tǒng)偏低.鄂爾多斯內(nèi)上地幔巖石圈速度明顯偏低的事實(shí)是清晰的,這可能與晚中生代到新生代的華北克拉通的活化而引起的熱物質(zhì)的流動有關(guān),具體分析將在后面詳細(xì)說明.

    二維剪切波速度差異分布圖(圖10)顯示出了不同深度處剪切波速度結(jié)構(gòu).淺層的剪切波速度與地表的地質(zhì)構(gòu)造特征吻合較好.莫霍面兩側(cè)的速度差異反映出了地幔巖石圈的結(jié)構(gòu)變化特征.可以看出鄂爾多斯塊體在莫霍面以下超過200km深處均為明顯的高速異常.

    鄂爾多斯塊體西南緣出現(xiàn)低速異常,深度超過200km.這里是青藏板塊擠壓處,隆起形成了六盤山.青藏板塊受印度板塊的擠壓作用,印度板塊仍以每年6cm的速度向北運(yùn)動,由此擠壓青藏板塊,青藏板塊遇到堅硬的鄂爾多斯板塊,引起地殼和地幔物質(zhì)沿邊界向東溢出.橫波分裂結(jié)果顯示這里的上地幔各向異性強(qiáng)烈,證明了地幔存在大量的物質(zhì)溢出[38-39].接收函數(shù)的結(jié)果顯示,青藏板塊與鄂爾多斯板塊的過渡帶地殼變形強(qiáng)烈,說明這里強(qiáng)烈的擠壓作用[40].我們的觀測結(jié)果顯示,這里的低速異常區(qū)深度約為200km,這預(yù)示著上地幔較大范圍內(nèi)可能存在熱物質(zhì)的流動.

    圖9 一維剪切波速度結(jié)構(gòu).實(shí)線為AK135參考模型,虛線為反演結(jié)果Fig.9 1-D average shear wave velocities in Ordos region(black dashed line).The solid line is the AK135 global model with three crustal layers.The thickness of the crust is 40km from the shear wave inversion

    中央轉(zhuǎn)換帶下的低速異常明顯,深度超過200km,印證了中央轉(zhuǎn)換帶內(nèi)新生代以來的活化導(dǎo)致了山西斷陷帶中的一系列斷陷盆地的形成[41-42]這一結(jié)論.橫波分裂的結(jié)果[13-14,43]顯示,中央轉(zhuǎn)換帶的各向異性方向與斷陷帶的走向一致,即為北東-南西向.而我們的結(jié)果給出了上地幔巖石圈的低速異常范圍,很好地補(bǔ)充了前人的研究.

    值得注意的是,中央轉(zhuǎn)換帶下的低速異常并不是連續(xù)的,在中段下,存在低速異常的間斷.這個間斷在莫霍面以下的各個圖像上均存在,超過200km深.這個間斷面的形成可能與整個華北克拉通的活化與減薄的成因有關(guān).

    垂向的剪切波速度分布顯示(圖11),鄂爾多斯地區(qū)整體的平均莫霍面深度為40km,山西斷陷帶和渭河斷陷帶的地殼出現(xiàn)較明顯的抬升和減薄,這符合拉張區(qū)的構(gòu)造特征.接收函數(shù)結(jié)果顯示,鄂爾多斯地區(qū)的莫霍面深度為40~43km,山西斷陷帶莫霍面比鄂爾多斯地塊薄5km以上[44-45].面波方法對莫霍面的深度約束并不強(qiáng),只能一定程度上反映莫霍面的變化范圍,我們的結(jié)果與接收函數(shù)的結(jié)果較為一致.六盤山、太行山下的地殼厚度較大,符合造山帶的特點(diǎn).太行山北段和華北平原的西緣下存在明顯的低速異常帶,這兩處異常帶都處在中央轉(zhuǎn)換帶上,華北克拉通的活化可能導(dǎo)致原本的板塊拼合帶出現(xiàn)了拉張減薄,引起地幔物質(zhì)上涌.

    圖10 三維剪切波速度異常分布圖.圖中數(shù)字表示出所在深度范圍.黑色實(shí)線表示縱向剖面圖Fig.10 Maps of shear wave velocity anomaly.The depth range of each map has been marked.The thick black lines are vertical profiles in further study

    圖11 剪切波速度垂直剖面.藍(lán)線表示地表地勢.黑線表示莫霍面位置.剖面位置在圖10中表示出.剖面AA′和BB′為經(jīng)向剖面,分別沿山西斷陷帶的東西兩側(cè).剖面CC′和DD′為緯向剖面,其中CC′剖面經(jīng)過渭河斷陷帶,DD′剖面橫貫鄂爾多斯塊體內(nèi)部Fig.11 Vertical profiles of shear wave velocity.Blue lines show the topography.Black lines show the positions of Moho.Positions of these profiles are shown in Fig.10.Profile AA′and BB′are along the western and eastern edges of Shanxi Rift.Profile CC′is along the Weihe Rift,and DD′passes through Ordos block.

    5 結(jié) 論

    我們利用雙平面波擬合法,得到了鄂爾多斯地區(qū)的瑞利波相速度結(jié)構(gòu),進(jìn)而反演得到剪切波三維速度結(jié)構(gòu).從雙平面波波場擬合方法的原理來看,該方法很大程度上還原了波場中的非平面波能量成分,使得反演結(jié)果更為準(zhǔn)確.在資料處理方面,拋棄了臺源大圓弧路徑這一條件,使得每一個地震事件在研究區(qū)域內(nèi)的全部臺站資料都可以得到有效利用,大大提高了反演的效率.

    短周期瑞利波相速度異常分布圖(圖8)很好地反映了地表的地質(zhì)構(gòu)造特征,太原盆地、臨汾盆地在新生代接受了大量沉積物,表現(xiàn)為低速異常,太行山、呂梁山為堅固穩(wěn)定的塊體,表現(xiàn)為高速異常.中長周期瑞利波相速度反映出上地幔巖石圈的構(gòu)造特征,為分析華北克拉通的地質(zhì)構(gòu)造背景提供了重要的依據(jù).

    瑞利波相速度結(jié)果和剪切波速度結(jié)果均顯示鄂爾多斯塊體具有超過200km的克拉通根基,堅固穩(wěn)定,沒有活化和改造的跡象.其西南緣上,青藏高原本身的地殼增厚和向東北方向的推擠,使得這里產(chǎn)生了沿著鄂爾多斯塊體邊界向東的地幔物質(zhì)流動.

    中部轉(zhuǎn)換帶下存在超過200km的低速異常并不是呈連續(xù)條帶分布,說明這里的地幔物質(zhì)上涌并不是連續(xù)貫穿整個中央轉(zhuǎn)換帶的,這種現(xiàn)象并不是單一因素造成的,很可能要受到太平洋板塊俯沖下插和青藏高原板塊的俯沖推覆的雙重作用.而這兩種作用各自的機(jī)制以及兩者之間的關(guān)聯(lián),還需要進(jìn)一步的觀測資料才能解釋.

    致 謝 感謝中國地震局地球物理研究所“國家數(shù)字測震臺網(wǎng)數(shù)據(jù)備份中心”為本研究提供地震波形數(shù)據(jù).感謝參與北京大學(xué)二維流動臺站布設(shè)的所有老師和同學(xué).感謝中國科學(xué)院地球物理研究所姜明明博士對本研究提出的意見和幫助.

    (References)

    [1] Zhao G,Sun M,Wilde S A,et al.Late Archean to Paleoproterozoic evolution of the North China Craton:key issues revisited.Precambrian Research,2005,136(2):177-202.

    [2] Wu C,Zhong C T.The Paleoproterozoic SW-NE collision model for the central North China Craton.Progress in Precambrian Research,1998,21(3):28-50.

    [3] Lu L,Xu X,Liu F.Early Precambrian Khondalite Series in North China.Changchun:Changchun Publishing House,1996:219-234.

    [4] Griffin W L,Andi Z,O'Reilly S Y,et al.Phanerozoic evolution of the lithosphere beneath the Sino-Korean Craton.∥Flower M F J,Chung S L,Lo C H,Lee T Y eds.Mantle Dynamics and Plate Interactions in East Asia.Washington:American Geophysical Union,1998:107-126.

    [5] Menzies M A,Xu Y G.Geodynamics of the North China Craton.∥Flower M F J,Chung S L,Lo C H,Lee T Y eds.Mantle Dynamics and Plate Interactions in East Asia.Washington:American Geophysical Union,1998:155-165.

    [6] 吳福元,葛文春,孫德有等.中國東部巖石圈減薄研究中的幾個問題.地學(xué)前緣,2003,10(3):51-60.Wu F Y,Ge W C,Sun D Y,et al.Discussions on the lithospheric thinning in Eastern China.Earth Science Frontier(in Chinese),2003,10(3):51-60.

    [7] Wu F,Xu Y,Gao S,et al.Lithospheric thinning and destruction of the North China Craton.Acta Petrologica Sinica,2008,24(6):1145-1174.

    [8] 陳國英,宋仲和,安昌強(qiáng)等.華北地區(qū)三維地殼上地幔結(jié)構(gòu).地球物理學(xué)報,1991,34(2):172-181.Chen G Y,Song Z H,An C Q,et al.Three-dimensional crust and upper mantle structure of the North China region.Chinese J.Geophys.(in Chinese),1991,34(2):172-181.

    [9] 易桂喜,姚華建,朱介壽等.中國大陸及鄰區(qū)Rayleigh面波相速度分布特征.地球物理學(xué)報,2008,51(2):402-411.Yi G X,Yao H J,Zhu J S,et al.Rayleigh-wave phase velocity distribution in China continent and its adjacent regions.Chinese J.Geophys.(in Chinese),2008,51(2):402-411.

    [10] Li C,van der Hilst R D,Toks?z M N.Constraining P-wave velocity variations in the upper mantle beneath Southeast Asia.Physics of the Earth and Planetary Interiors,2006,154(2):180-195.

    [11] Zheng T Y,Zhao L,Zhu R X.Insight into the geodynamics of cratonic reactivation from seismic analysis of the crustmantle boundary.Geophys.Res.Lett.2008,35(8):L08303,doi:10.1029/2008GL033439.

    [12] Zheng T Y,Chen L,Zhao L,et al.Crust-mantle structure difference across the gravity gradient zone in North China Craton:Seismic image of the thinned continental crust.Physics of the Earth and Planetary Interiors,2006,159(1-2):43-58.

    [13] Zhao L,Zheng T.Complex upper-mantle deformation beneath the North China Craton:implications for lithospheric thinning.Geophysical Journal International,2007,170(3):1095-1099.

    [14] Zhao L,Allen R M,Zheng T,et al.Reactivation of an Archean craton:Constraints from P-and S-wave tomography in North China.Geophys Res Lett,2009,36(17):L17306,doi:1011029/2009GL039781.

    [15] 姚華建,徐果明,肖翔等.基于圖像分析的雙臺面波相速度頻散曲線快速提取方法.地震地磁觀測與研究,2004,25(1):1-8.Yao H J,Xu G M,Xiao X,et al.A quick tracing method based on image analysis technique for the determination of dual stations phase velocities dispersion curve of surface wave.Seismol.Geomagn.Observ.Res.(in Chinese),2004,25(1):1-8.

    [16] Friederich W.Wave-theoretical inversion of teleseismic surface waves in a regional network:phase-velocity maps and a three-dimensionalupper-mantle shear-wave-velocity model for southern Germany.Geophysical Journal International,1998,132(1):203-225.

    [17] Friederich W,Wielandt E,Stange S.Multiple forward scattering of surface waves:Comparison with an exact solution and Born single-scattering methods.Geophysical Journal International,1993,112(2):264-275.

    [18] Friederich W,Wielandt E.Interpretation of seismic surfacewaves in regional networks-joint estimation of wave-field geometry and local phase-velocity-method and numerical tests.Geophysical Journal International,1995,120(3):731-744.

    [19] Wielandt E.Propagation and structural interpretation of nonplane waves.Geophysical Journal International,1993,113(1):45-53.

    [20] Forsyth D W,Li A.Array analysis of two-dimensional variations in surface wave phase velocity and azimuthal anisotropy in the presence of multipathing interference.//Levander A,Nolet G eds.Seismic Earth:Array Analysis of Broad-band Seismograms.Washington:American Geophysical Union.2005:81-97.

    [21] Jiang M M,Zhou S Y,et al.3-D lithospheric structure beneath southern Tibet from Rayleigh-wave tomography with a 2-D seismic array.Geophysical Journal International,2010,185(2):593-608.

    [22] 李鵬,周仕勇,陳永順等.利用雙平面波干涉面波層析成像方法研究山西斷陷盆地及鄂爾多斯地臺三維速度結(jié)構(gòu).CT理論與應(yīng)用研究,2010,19(3):47-60.Li P,Zhou S Y,Chen Y S,et al.3Dvelocity structure in Shanxi graben and Ordos from two plane waves method.CT Theory and Applications(in Chinese),2010,19(3):47-60.

    [23] 鄭秀芬,歐陽飚,張東寧等.“國家數(shù)字測震臺網(wǎng)數(shù)據(jù)備份中心”技術(shù)系統(tǒng)建設(shè)及其對汶川大地震研究的數(shù)據(jù)支撐.地球物理學(xué)報,2009,52(5):1412-1417.Zheng X F,Ouyang B,Zhang D N,et al.Technical system construction of Data Backup Centre for China Seismograph Network and the data support to researches on the Wenchuan earthquake.Chinese J.Geophys.(in Chinese),2009,52(5):1412-1417.

    [24] 姜明明.藏南巖石圈結(jié)構(gòu)研究及地震活動性統(tǒng)計建模[博士論文].北京:北京大學(xué)地球物理系,2008.Jiang M M.Studies on the lithospheric stucture in souther Tibet and modeling the seismicity of north China by statistics[Ph.D.thesis](in Chinese).Beijing:Geophysics Department of Peking University,2008.

    [25] Press W H.Numerical Recipes in FORTRAN:the art of Scientific Computing.Cambridge:Cambridge Univ Pr,1992.

    [26] Tarantola A,Valette B.Generalized nonlinear inverse problems solved using the least squares criterion.Rev.Geophys.Space Phys.,1982,20(2):219-232.

    [27] Li A,F(xiàn)orsyth D W,F(xiàn)ischer K M.Shear velocity structure and azimuthal anisotropy beneath eastern North America from Rayleigh wave inversion.J.Geophys.Res.,2003,108(B8):2362.

    [28] Kennett B L N,Engdahl E R,Buland R.Constraints on seismic velocities in the Earth from traveltimes.Geophysical Journal International,1995,122(1):108-124.

    [29] Liu D Y,Nutman A P,Compston W,et al.Remnants of≥3800Ma crust in the Chinese part of the Sino-Korean craton.Geology,1992,20(4):339-342.

    [30] Zhao G,Wilde S A,Cawood P A,et al.Archean blocks and their boundaries in the North China Craton:lithological,geochemical,structural and PT path constraints and tectonic evolution.Precambrian Research,2001,107(1-2):45-73.

    [31] Ren J,Tamaki K,Li S,et al.Late Mesozoic and Cenozoic rifting and its dynamic setting in Eastern China and adjacent areas.Tectonophysics,2002,344(3-4):175-205.

    [32] Zhao L,Zheng T.Seismic structure of the Bohai Bay Basin,northern China:Implications for basin evolution.Earthand PlanetaryScienceLetters,2005,231(1-2):9-22.

    [33] 臧紹先,劉永剛,寧杰遠(yuǎn).華北地區(qū)巖石圈熱結(jié)構(gòu)的研究.地球物理學(xué)報,2002,45(1):56-66.Zang S X,Liu Y G,Ning J Y.Thermal structure of the lithosphere in North China.Chinese J.Geophys.(in Chinese),2002,45(1):56-66.

    [34] Forsyth D W,Webb S C,Dorman L M,et al.Phase velocities of Rayleigh waves in the MELT experiment on the East Pacific Rise.Science,1998,280(5367):1235-1238.

    [35] Weeraratne D S,F(xiàn)orsyth D W,F(xiàn)ischer K M,et al.Evidence for an upper mantle plume beneath the Tanzanian craton from Rayleigh wave tomography.J.Geophys.Res.,2003,108(B9):2427.

    [36] Saito M.DISPER80:A subroutine package for the calculation of seismic normal-mode solution.//Doornbos D J eds.Seismological Algorithms:Computational Methods and Computure Programs.New York:Elsevier.1988:293-319.

    [37] Jiang M M,Zhou S Y,Sandvol E,et al.3-D lithospheric structure beneath southern Tibet from Rayleigh-wave tomography with a 2-D seismic array.Geophysical Journal International,2011,185(2):593-608.

    [38] Yin A,Harrison T M.Geologic evolution of the Himalayan-Tibetan orogen.Annual Review of Earth and Planetary Sciences,2000,28(1):211-280.

    [39] Bilham R,Larson K,F(xiàn)reymueller J,et al.GPS measurements of present-day convergence across the Nepal Himalaya.Nature,1997,386(6620):61-64.

    [40] 陳九輝,劉啟元,李順成等.青藏高原東北緣—鄂爾多斯地塊地殼上地幔S波速度結(jié)構(gòu).地球物理學(xué)報,2005,48(2):333-342.Chen J H,Liu Q Y,Li S C,et al.Crust and upper mantle S 2wave velocity structure across Northeastern Tibetan Plateau and Ordos block.Chinese J.Geophys.(in Chinese),2005,48(2):333-342.

    [41] Zhang Y Q,Mercier J L,Vergély P.Extension in the graben systems around the Ordos(China),and its contribution to the extrusion tectonics of south China with respect to Gobi-Mongolia.Tectonophysics,1998,285(1):41-75.

    [42] Zhang Y,Ma Y,Yang N,et al.Cenozoic extensional stress evolution in North China.Journal of Geodynamics,2003,36(5):591-613.

    [43] Zhao L,Zheng T.Using shear wave splitting measurements to investigate the upper mantle anisotropy beneath the North China Craton:Distinct variation from east to west.Geophys.Res.Lett.,2005,32(10):L10309,doi:10.1029/2005GL022585.

    [44] 唐有彩,馮永革,陳永順等.山西斷陷帶地殼結(jié)構(gòu)的接收函數(shù)研究.地球物理學(xué)報,2010,53(9):2102-2109.Tang Y C,F(xiàn)eng Y G,Chen Y S,et al.Receiver function analysis at Shanxi Rift.Chinese J.Geophys.(in Chinese),2010,53(9):2102-2109.

    [45] 田寶峰,李娟,王為民等.華北太行山區(qū)地殼各向異性的接收函數(shù)證據(jù).地球物理學(xué)報,2008,51(5):1459-1467.Tian B F,Li J,Wang W M,et al.Crust anisotropy of Taihangshan mountain range in north China inferred from receiver functions.Chinese J.Geophys.(in Chinese),2008,51(5):1459-1467.

    3-D lithospheric structure of upper mantle beneath Ordos region from Rayleigh-wave tomography

    LI Duo,ZHOU Shi-Yong*,CHEN Yong-Shun,F(xiàn)ENG Yong-Ge,LI Peng
    Institute of Theoretical and Applied Geophysics,School of Earth and Space Science,Peking University,Beijing100871,China

    We apply two-plane-wave tomography,which takes the influences of the non-plane wavefield into consideration,to study the structure of velocity of the upper mantle lithosphere of Ordos region.The source of data consists of two parts,one is from 69broad-band seismographs of China Earthquake Networks Center,and the other is from 34mobile digital broad-band seismographs of PKU.At short periods most high and low velocity anomalies correlate well with surface geological features.The difference of the structure of upper mantle lithosphere is revealed by long-period surface waves.We extracted the 3-D structure of shear velocity anomalies of 200 km depth from the Rayleigh wave phase velocities.The results reveal that the Ordos block has the high velocity beyond 200km depth,and no evidence for reactivation.There is fierce interaction between Ordos block and Tibet block at the southwestern edge of Ordos block,which causes theupper mantle flow there.Upwelling beneath the Central zone may be caused by both the subduction of Pacific plate and India plate.

    Two-plane-wave method,Surface wave tomography,Ordos,North China Craton,Lithosphere

    10.6038/j.issn.0001-5733.2012.05.019

    P315

    2012-01-14,2012-03-17收修定稿

    地震行業(yè)科研專項(xiàng)課題(201208009)與國家自然科學(xué)基金(41074030,90814002)聯(lián)合資助.

    李多,女,1988年生,碩士研究生在讀,主要從事地震學(xué)方面的研究.E-mail:liduoduo07@gmail.com

    *通訊作者周仕勇,男,教授,主要從事震源過程、定位及地震活動性研究.E-mail:zsy@pku.edu.cn

    李多,周仕勇,陳永順等.鄂爾多斯地區(qū)上地幔巖石圈三維速度結(jié)構(gòu)面波反演研究.地球物理學(xué)報,2012,55(5):1613-1623,

    10.6038/j.issn.0001-5733.2012.05.019.

    Li D,Zhou S Y,Chen Y S,et al.3-D lithospheric structure of upper mantle beneath Ordos region from Rayleigh-wave tomography.Chinese J.Geophys.(in Chinese),2012,55(5):1613-1623,doi:10.6038/j.issn.0001-5733.2012.05.019.

    (本文編輯 胡素芳)

    猜你喜歡
    平面波克拉通面波
    巖石圈地幔分層性對克拉通穩(wěn)定性的影響
    Landau-Lifshitz方程平面波解的全局光滑性
    5G OTA測量寬帶平面波模擬器的高效優(yōu)化方法與應(yīng)用
    gPhone重力儀的面波頻段響應(yīng)實(shí)測研究
    地震研究(2021年1期)2021-04-13 01:04:56
    自適應(yīng)相減和Curvelet變換組合壓制面波
    有關(guān)克拉通破壞及其成因的綜述
    華北克拉通重力剖面重力點(diǎn)位GPS測量精度分析
    基于GPU并行運(yùn)算的超聲平面波成像仿真
    電子制作(2016年11期)2016-11-07 08:43:45
    拉張槽對四川盆地海相油氣分布的控制作用
    淺析工程勘探的面波勘探方法
    河南科技(2014年8期)2014-02-27 14:07:40
    搡老岳熟女国产| 又紧又爽又黄一区二区| 亚洲五月婷婷丁香| 精品久久久久久久人妻蜜臀av| 精品一区二区三区人妻视频| 欧美乱码精品一区二区三区| 日本一本二区三区精品| 嫁个100分男人电影在线观看| 免费观看人在逋| 亚洲国产欧美人成| 一区福利在线观看| 五月伊人婷婷丁香| 一a级毛片在线观看| 成人国产综合亚洲| 女人被狂操c到高潮| 99国产精品一区二区蜜桃av| 亚洲电影在线观看av| 亚洲中文日韩欧美视频| 韩国av一区二区三区四区| 99国产精品一区二区三区| 欧美区成人在线视频| 婷婷精品国产亚洲av在线| 国产免费av片在线观看野外av| 午夜a级毛片| 精品国产美女av久久久久小说| 成人国产一区最新在线观看| 此物有八面人人有两片| 欧美bdsm另类| 久久久成人免费电影| 五月玫瑰六月丁香| 长腿黑丝高跟| 99精品欧美一区二区三区四区| 久久国产乱子伦精品免费另类| 亚洲va日本ⅴa欧美va伊人久久| 欧美绝顶高潮抽搐喷水| 欧美不卡视频在线免费观看| xxxwww97欧美| 别揉我奶头~嗯~啊~动态视频| 亚洲人成伊人成综合网2020| 欧美一级a爱片免费观看看| 国产成人a区在线观看| 搞女人的毛片| 麻豆久久精品国产亚洲av| 村上凉子中文字幕在线| 美女高潮喷水抽搐中文字幕| 日日干狠狠操夜夜爽| 中亚洲国语对白在线视频| 少妇人妻精品综合一区二区 | 久久精品亚洲精品国产色婷小说| 亚洲av成人精品一区久久| 国产熟女xx| 欧美三级亚洲精品| 国产三级在线视频| 97碰自拍视频| 日韩国内少妇激情av| 天堂网av新在线| 亚洲欧美日韩高清在线视频| 91在线观看av| 韩国av一区二区三区四区| 国内精品一区二区在线观看| 成人高潮视频无遮挡免费网站| 一个人看的www免费观看视频| 日本熟妇午夜| 亚洲人成网站在线播放欧美日韩| 午夜福利高清视频| 两个人的视频大全免费| 国产一区二区三区视频了| 村上凉子中文字幕在线| 99久久久亚洲精品蜜臀av| 国产色爽女视频免费观看| 国产v大片淫在线免费观看| 中文字幕精品亚洲无线码一区| 欧美激情在线99| 一个人免费在线观看的高清视频| 一夜夜www| 免费观看的影片在线观看| 日韩有码中文字幕| 国产精品永久免费网站| 亚洲国产欧美网| 欧美绝顶高潮抽搐喷水| 在线观看66精品国产| 国产真实伦视频高清在线观看 | 九色成人免费人妻av| 人人妻人人澡欧美一区二区| 91久久精品电影网| 日本免费一区二区三区高清不卡| 人人妻人人澡欧美一区二区| 久久精品夜夜夜夜夜久久蜜豆| 亚洲av成人av| 国内毛片毛片毛片毛片毛片| 午夜亚洲福利在线播放| 一边摸一边抽搐一进一小说| 嫁个100分男人电影在线观看| 午夜福利视频1000在线观看| 青草久久国产| 午夜激情欧美在线| 国产激情偷乱视频一区二区| 中文资源天堂在线| 精品国产亚洲在线| 97人妻精品一区二区三区麻豆| 91麻豆精品激情在线观看国产| 男女之事视频高清在线观看| 国产精品美女特级片免费视频播放器| 国产久久久一区二区三区| 成人永久免费在线观看视频| 两个人的视频大全免费| 日本精品一区二区三区蜜桃| 男女视频在线观看网站免费| 狂野欧美激情性xxxx| 在线播放无遮挡| 国产午夜福利久久久久久| 成人三级黄色视频| 欧美黄色片欧美黄色片| 日本精品一区二区三区蜜桃| 一二三四社区在线视频社区8| 久久香蕉国产精品| 我要搜黄色片| 精品欧美国产一区二区三| 中文字幕人妻熟人妻熟丝袜美 | 欧美又色又爽又黄视频| 黄色女人牲交| 欧美性猛交黑人性爽| 免费人成视频x8x8入口观看| 天天躁日日操中文字幕| 国产视频内射| 法律面前人人平等表现在哪些方面| 高潮久久久久久久久久久不卡| 特级一级黄色大片| 亚洲中文字幕一区二区三区有码在线看| 亚洲欧美精品综合久久99| 特大巨黑吊av在线直播| 亚洲国产日韩欧美精品在线观看 | 国产真实乱freesex| 国产色婷婷99| 精品国产超薄肉色丝袜足j| 美女高潮喷水抽搐中文字幕| a在线观看视频网站| 午夜激情福利司机影院| 国产精品永久免费网站| 亚洲va日本ⅴa欧美va伊人久久| 嫩草影视91久久| 久久久久久国产a免费观看| 国产美女午夜福利| 高潮久久久久久久久久久不卡| 黄色日韩在线| 成人三级黄色视频| 久久久久久大精品| 真人一进一出gif抽搐免费| 国内久久婷婷六月综合欲色啪| 成人永久免费在线观看视频| 国内少妇人妻偷人精品xxx网站| 国产精品免费一区二区三区在线| 日日干狠狠操夜夜爽| 欧美丝袜亚洲另类 | 美女高潮喷水抽搐中文字幕| 母亲3免费完整高清在线观看| 精品午夜福利视频在线观看一区| 欧美乱码精品一区二区三区| 好男人电影高清在线观看| 免费看美女性在线毛片视频| 日本一本二区三区精品| 婷婷精品国产亚洲av在线| 亚洲国产精品成人综合色| 亚洲在线观看片| 2021天堂中文幕一二区在线观| 欧美黑人巨大hd| 国产色婷婷99| 日本 av在线| 真人一进一出gif抽搐免费| 在线观看舔阴道视频| 热99在线观看视频| 免费在线观看亚洲国产| 精品国产亚洲在线| 欧美zozozo另类| 精品一区二区三区av网在线观看| www日本黄色视频网| 美女cb高潮喷水在线观看| 亚洲午夜理论影院| 欧美黄色片欧美黄色片| 午夜福利成人在线免费观看| 国产乱人视频| 国模一区二区三区四区视频| 欧美日韩国产亚洲二区| 亚洲18禁久久av| 欧美黄色淫秽网站| 国产精品亚洲一级av第二区| 舔av片在线| 国产精品永久免费网站| 色在线成人网| 精品欧美国产一区二区三| 无限看片的www在线观看| 免费电影在线观看免费观看| 99久久综合精品五月天人人| 91字幕亚洲| 女人十人毛片免费观看3o分钟| 夜夜躁狠狠躁天天躁| 午夜影院日韩av| 一区二区三区免费毛片| 国产探花极品一区二区| 有码 亚洲区| 亚洲人成网站在线播| 精品人妻1区二区| 99国产综合亚洲精品| 黄色丝袜av网址大全| 亚洲 国产 在线| 夜夜爽天天搞| 中文字幕高清在线视频| 国产高清视频在线观看网站| 丰满人妻熟妇乱又伦精品不卡| 91字幕亚洲| 老鸭窝网址在线观看| 麻豆一二三区av精品| 国产精品香港三级国产av潘金莲| x7x7x7水蜜桃| 精品国产超薄肉色丝袜足j| 制服丝袜大香蕉在线| 国产成人a区在线观看| 精品久久久久久,| 久久性视频一级片| 网址你懂的国产日韩在线| 欧美日韩中文字幕国产精品一区二区三区| 久久久色成人| 一个人看的www免费观看视频| 亚洲无线在线观看| 亚洲av成人精品一区久久| 亚洲精品日韩av片在线观看 | 老熟妇乱子伦视频在线观看| 18禁在线播放成人免费| 两个人看的免费小视频| 五月伊人婷婷丁香| 国产97色在线日韩免费| 嫩草影视91久久| www.www免费av| 好看av亚洲va欧美ⅴa在| 亚洲精品色激情综合| 国产精品亚洲av一区麻豆| 欧美日韩黄片免| 国产高清有码在线观看视频| 香蕉久久夜色| 亚洲精品粉嫩美女一区| 日本在线视频免费播放| 成人三级黄色视频| 亚洲成a人片在线一区二区| 日韩中文字幕欧美一区二区| 亚洲黑人精品在线| 又黄又粗又硬又大视频| 嫩草影院精品99| 久久久精品大字幕| 久久久久久久久大av| 两人在一起打扑克的视频| 久久精品91蜜桃| 欧美激情在线99| 久久久久久久久中文| 宅男免费午夜| 色av中文字幕| 国产野战对白在线观看| 亚洲人成伊人成综合网2020| 精品人妻偷拍中文字幕| 真人一进一出gif抽搐免费| 老司机在亚洲福利影院| 成人午夜高清在线视频| 国产亚洲精品av在线| www.999成人在线观看| 亚洲不卡免费看| 人人妻人人澡欧美一区二区| 亚洲人与动物交配视频| 国产免费一级a男人的天堂| 久久亚洲真实| 天天一区二区日本电影三级| 99国产综合亚洲精品| 中文字幕久久专区| 免费av毛片视频| 狂野欧美激情性xxxx| 中文字幕av成人在线电影| 成人无遮挡网站| 国产探花极品一区二区| 香蕉丝袜av| 成人一区二区视频在线观看| 深夜精品福利| 两个人看的免费小视频| 色吧在线观看| 国产高清videossex| 国产成人系列免费观看| 国产视频内射| 日本五十路高清| tocl精华| 999久久久精品免费观看国产| 久久亚洲真实| 精品福利观看| 免费人成视频x8x8入口观看| 日韩精品中文字幕看吧| 亚洲欧美一区二区三区黑人| 内射极品少妇av片p| 欧美一区二区国产精品久久精品| 成人一区二区视频在线观看| АⅤ资源中文在线天堂| 99热这里只有是精品50| 高清日韩中文字幕在线| 国产在视频线在精品| 国产真人三级小视频在线观看| 国产日本99.免费观看| 人妻夜夜爽99麻豆av| 中国美女看黄片| 国产毛片a区久久久久| 超碰av人人做人人爽久久 | 久久精品国产自在天天线| 久久久久性生活片| 亚洲精品在线观看二区| 少妇人妻精品综合一区二区 | 一级a爱片免费观看的视频| 热99re8久久精品国产| 日本黄大片高清| 国产精品久久视频播放| 国语自产精品视频在线第100页| 18禁美女被吸乳视频| 免费av观看视频| 亚洲内射少妇av| 亚洲人成网站高清观看| 五月玫瑰六月丁香| 怎么达到女性高潮| 免费高清视频大片| 免费观看精品视频网站| 亚洲熟妇熟女久久| 国产激情欧美一区二区| 床上黄色一级片| 久久国产精品人妻蜜桃| 色综合亚洲欧美另类图片| 欧美极品一区二区三区四区| 中文字幕人成人乱码亚洲影| 最新在线观看一区二区三区| 日韩成人在线观看一区二区三区| 欧美3d第一页| 欧美av亚洲av综合av国产av| 久久草成人影院| 亚洲美女黄片视频| 成年女人毛片免费观看观看9| 女人高潮潮喷娇喘18禁视频| 国产69精品久久久久777片| 亚洲精品在线观看二区| 国产97色在线日韩免费| 亚洲国产精品999在线| 女警被强在线播放| 91字幕亚洲| 熟女人妻精品中文字幕| 97碰自拍视频| e午夜精品久久久久久久| 此物有八面人人有两片| 国产视频内射| 亚洲精品日韩av片在线观看 | 久久精品人妻少妇| x7x7x7水蜜桃| 午夜福利18| av中文乱码字幕在线| 精品熟女少妇八av免费久了| 嫩草影视91久久| 国产麻豆成人av免费视频| 日本免费一区二区三区高清不卡| 中文字幕精品亚洲无线码一区| 国内少妇人妻偷人精品xxx网站| 一本一本综合久久| 国产一区二区三区视频了| 毛片女人毛片| 三级国产精品欧美在线观看| 午夜精品久久久久久毛片777| 悠悠久久av| 1024手机看黄色片| 91久久精品电影网| 亚洲人成伊人成综合网2020| 色在线成人网| 欧美另类亚洲清纯唯美| 九九久久精品国产亚洲av麻豆| 一卡2卡三卡四卡精品乱码亚洲| av国产免费在线观看| 国产日本99.免费观看| 午夜免费成人在线视频| 老司机午夜福利在线观看视频| 国产真实伦视频高清在线观看 | 女警被强在线播放| 亚洲真实伦在线观看| aaaaa片日本免费| 国产精品久久久久久亚洲av鲁大| 夜夜夜夜夜久久久久| 国产野战对白在线观看| 老司机福利观看| 在线观看一区二区三区| 色尼玛亚洲综合影院| 国内揄拍国产精品人妻在线| 午夜久久久久精精品| 亚洲,欧美精品.| 亚洲av成人精品一区久久| 十八禁网站免费在线| 欧美最黄视频在线播放免费| 丰满乱子伦码专区| 国产精品亚洲美女久久久| 午夜亚洲福利在线播放| 精品福利观看| 国产亚洲精品久久久com| 亚洲一区二区三区色噜噜| 国产精华一区二区三区| 91久久精品电影网| e午夜精品久久久久久久| 在线视频色国产色| 久久伊人香网站| 久久久久久久久大av| xxx96com| 亚洲精品国产精品久久久不卡| a级一级毛片免费在线观看| 国产亚洲欧美在线一区二区| 香蕉久久夜色| 99久久无色码亚洲精品果冻| 欧美区成人在线视频| 成年人黄色毛片网站| 老汉色∧v一级毛片| 一级作爱视频免费观看| 人人妻人人澡欧美一区二区| 欧美另类亚洲清纯唯美| 黄片小视频在线播放| xxx96com| 久久久国产精品麻豆| 国产熟女xx| 乱人视频在线观看| 国产精品野战在线观看| 亚洲国产欧洲综合997久久,| 在线国产一区二区在线| 久久久国产精品麻豆| 国产在线精品亚洲第一网站| 97超视频在线观看视频| 午夜视频国产福利| 少妇裸体淫交视频免费看高清| 亚洲欧美日韩东京热| 亚洲精品影视一区二区三区av| 全区人妻精品视频| 国产美女午夜福利| 97人妻精品一区二区三区麻豆| 亚洲国产中文字幕在线视频| 男女床上黄色一级片免费看| 中文字幕高清在线视频| 欧美成人性av电影在线观看| 麻豆成人av在线观看| tocl精华| 男女午夜视频在线观看| 国产久久久一区二区三区| 97超视频在线观看视频| 色综合婷婷激情| 18禁美女被吸乳视频| 亚洲专区国产一区二区| 久9热在线精品视频| 日韩欧美国产一区二区入口| 欧美极品一区二区三区四区| 美女高潮喷水抽搐中文字幕| 欧美高清成人免费视频www| 天美传媒精品一区二区| 亚洲人成网站在线播放欧美日韩| 国产美女午夜福利| 欧美在线黄色| 真实男女啪啪啪动态图| 精品久久久久久久毛片微露脸| 日本熟妇午夜| 美女cb高潮喷水在线观看| 日本免费一区二区三区高清不卡| 欧美中文日本在线观看视频| 国产欧美日韩一区二区三| 国模一区二区三区四区视频| 免费在线观看影片大全网站| 日本撒尿小便嘘嘘汇集6| 一a级毛片在线观看| 国产色爽女视频免费观看| 久久精品综合一区二区三区| 天天一区二区日本电影三级| 午夜福利视频1000在线观看| 欧美又色又爽又黄视频| 亚洲不卡免费看| 一二三四社区在线视频社区8| 亚洲国产日韩欧美精品在线观看 | 悠悠久久av| 桃色一区二区三区在线观看| 男女做爰动态图高潮gif福利片| 99热这里只有是精品50| 国产精品自产拍在线观看55亚洲| 黄色日韩在线| 别揉我奶头~嗯~啊~动态视频| 亚洲美女黄片视频| 中文在线观看免费www的网站| 国产一区二区在线av高清观看| 欧美+日韩+精品| 一边摸一边抽搐一进一小说| 欧美+亚洲+日韩+国产| 欧美成人免费av一区二区三区| 亚洲,欧美精品.| 在线国产一区二区在线| 国产精品 国内视频| 亚洲一区二区三区色噜噜| 97超级碰碰碰精品色视频在线观看| eeuss影院久久| 精品一区二区三区av网在线观看| 叶爱在线成人免费视频播放| 久久国产精品影院| 一级黄片播放器| 亚洲精品日韩av片在线观看 | 国产精品香港三级国产av潘金莲| 一级毛片女人18水好多| а√天堂www在线а√下载| 国产视频一区二区在线看| 51国产日韩欧美| 午夜a级毛片| 夜夜躁狠狠躁天天躁| 免费看美女性在线毛片视频| 午夜两性在线视频| 桃色一区二区三区在线观看| 18禁国产床啪视频网站| 中文资源天堂在线| 俺也久久电影网| 99国产极品粉嫩在线观看| 法律面前人人平等表现在哪些方面| 少妇的逼好多水| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲最大成人手机在线| 午夜福利18| 真实男女啪啪啪动态图| 超碰av人人做人人爽久久 | 亚洲aⅴ乱码一区二区在线播放| 国产精品久久久久久久电影 | 香蕉av资源在线| 日韩欧美在线乱码| 日韩欧美精品免费久久 | 欧美性感艳星| 日本黄色视频三级网站网址| 性欧美人与动物交配| 久久久久久久久久黄片| 搞女人的毛片| 亚洲欧美日韩高清专用| 国产真人三级小视频在线观看| 日韩大尺度精品在线看网址| 亚洲七黄色美女视频| 我的老师免费观看完整版| 成人鲁丝片一二三区免费| 国产精品,欧美在线| 在线观看舔阴道视频| 十八禁人妻一区二区| 国产三级在线视频| 99热这里只有是精品50| 51午夜福利影视在线观看| 欧美日韩黄片免| 99久国产av精品| 国产精品免费一区二区三区在线| 天堂网av新在线| 亚洲国产欧美人成| 国产午夜精品论理片| 亚洲av熟女| 一级黄色大片毛片| 18禁裸乳无遮挡免费网站照片| 悠悠久久av| 国产视频内射| av国产免费在线观看| 日韩欧美在线二视频| 宅男免费午夜| 脱女人内裤的视频| 久久久久免费精品人妻一区二区| 久久久精品欧美日韩精品| 18+在线观看网站| 日本a在线网址| 久久国产精品影院| 精品欧美国产一区二区三| 国产乱人视频| 成人鲁丝片一二三区免费| 91麻豆精品激情在线观看国产| 国产亚洲精品av在线| 免费观看的影片在线观看| 18禁黄网站禁片免费观看直播| 国产成人av激情在线播放| 女人十人毛片免费观看3o分钟| 老汉色∧v一级毛片| www国产在线视频色| 99在线人妻在线中文字幕| 精品福利观看| xxxwww97欧美| www.999成人在线观看| 99久久九九国产精品国产免费| 午夜激情福利司机影院| 国产精品,欧美在线| 国产精品乱码一区二三区的特点| 在线观看免费午夜福利视频| 丰满的人妻完整版| 欧美乱色亚洲激情| 成人性生交大片免费视频hd| 久久精品国产亚洲av香蕉五月| 在线国产一区二区在线| 免费观看的影片在线观看| 国产一区二区在线观看日韩 | 久久这里只有精品中国| 18+在线观看网站| 给我免费播放毛片高清在线观看| 欧美bdsm另类| 久久精品国产综合久久久| 午夜精品在线福利| 国产主播在线观看一区二区| 国产亚洲精品综合一区在线观看| 青草久久国产| 国内揄拍国产精品人妻在线| 亚洲国产日韩欧美精品在线观看 | 日韩欧美精品v在线| 窝窝影院91人妻| 99riav亚洲国产免费| 精华霜和精华液先用哪个| 波多野结衣巨乳人妻| a在线观看视频网站| 3wmmmm亚洲av在线观看| 中文字幕av在线有码专区| 青草久久国产| 国产高清激情床上av| 成人永久免费在线观看视频| 操出白浆在线播放| 在线免费观看的www视频| 99久久九九国产精品国产免费| av片东京热男人的天堂|