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

    基于背景噪聲研究華北克拉通中部Rayleigh波相速度和方位各向異性

    2021-08-03 10:56:48黃翔丁志峰寧杰遠(yuǎn)常利軍
    地球物理學(xué)報(bào) 2021年8期
    關(guān)鍵詞:克拉通華北方位

    黃翔, 丁志峰 , 寧杰遠(yuǎn), 常利軍

    1 中國(guó)地震局地球物理研究所, 北京 100081 2 北京大學(xué)地球與空間科學(xué)學(xué)院, 北京 100871

    0 引言

    華北克拉通是世界上最古老的太古代克拉通之一,與其他穩(wěn)定且擁有巨厚巖石圈的典型克拉通不同,其東部的華北盆地,在顯生宙經(jīng)歷了顯著的熱構(gòu)造活化和巖石圈破壞,其巖石圈地幔以及地殼都受到了不同程度的改造和減薄(Menzies et al., 1993; Griffin et al., 1998; Zheng et al., 2006),其演化過(guò)程和破壞機(jī)制一直是地球科學(xué)界關(guān)注的熱點(diǎn)問(wèn)題.前人研究認(rèn)為在1.8 Ga前華北克拉通由其東部塊體和西部塊體碰撞合并而成,在其中部形成了華北中部造山帶(Trans-North China Orogen)(Zhao et al., 2005),主要由太行造山帶、呂梁造山帶以及山西地塹等構(gòu)造單元組成.一般認(rèn)為華北克拉通西部的鄂爾多斯塊體仍然保留了其克拉通特性,內(nèi)部比較穩(wěn)定,而位于鄂爾多斯塊體與華北盆地中間的中部造山帶則表現(xiàn)出更為復(fù)雜的構(gòu)造特征(Ai et al., 2019).前人在該區(qū)域進(jìn)行的層析成像和各向異性研究大都基于固定臺(tái)站(Li et al., 2011; Bao et al., 2013; Li S L et al., 2018; Lü, 2019; Ai et al., 2020; 顧勤平等, 2020a)或只能覆蓋研究區(qū)域局部范圍的臺(tái)陣(Chen et al., 2015; Zheng et al., 2019),受限于臺(tái)站覆蓋范圍及密度不足等原因,目前對(duì)于整個(gè)華北克拉通中部造山帶高分辨率的三維速度結(jié)構(gòu)以及復(fù)雜的構(gòu)造活動(dòng)和變形特征的認(rèn)識(shí)還比較有限.

    山西地塹位于太行造山帶與呂梁造山帶之間,北起燕山造山帶,南連渭河地塹,是華北克拉通中部最重要的組成單元之一,其速度結(jié)構(gòu)和構(gòu)造運(yùn)動(dòng)歷史都非常復(fù)雜,歷來(lái)都是地球物理學(xué)以及地質(zhì)學(xué)等地球科學(xué)領(lǐng)域的重點(diǎn)研究區(qū).前人基于地質(zhì)調(diào)查研究認(rèn)為渭河—山西裂谷系統(tǒng)是在漸新世至上新世期間,從渭河裂谷開(kāi)始逐步發(fā)展到山西裂谷中部至北部,即裂谷作用是由南向北進(jìn)行的(Zhang et al., 1998).然而,該系統(tǒng)的演化發(fā)展機(jī)制還存有爭(zhēng)議,如對(duì)于造成華北克拉通的伸展變形的原因,有的學(xué)者認(rèn)為主要與亞洲大陸東緣下方太平洋板塊的西向俯沖有關(guān)(Northrup et al., 1995),另一部分學(xué)者則認(rèn)為印度—亞洲板塊碰撞的遠(yuǎn)場(chǎng)效應(yīng)是引起華北克拉通的大陸內(nèi)部裂谷作用的主要原因(Molnar and Tapponnier, 1977; Liu et al., 2004),也有學(xué)者將兩種原因結(jié)合討論其對(duì)華北克拉通伸展作用的影響(Zhang et al., 1998; Xu et al., 2004; Xu, 2007).Ai等(2019)基于背景噪聲和地震面波反演得到的S波速度模型認(rèn)為渭河—山西裂谷系統(tǒng)南北段的驅(qū)動(dòng)機(jī)制不同,其中北段目前主要是由地幔物質(zhì)上涌驅(qū)動(dòng)的,而南段則主要受青藏高原持續(xù)隆升的影響.

    大同火山群位于山西地塹的最北端,是華北克拉通最大的第四紀(jì)板內(nèi)火山區(qū)之一(陳文寄等, 1992),在早更新世(0.74 Ma)首次噴發(fā),巖漿活動(dòng)持續(xù)到晚更新世(陳文寄等, 1992; Xu et al., 2005).雖然巖漿規(guī)模相對(duì)較小(Xu and Ma, 1992),但作為板內(nèi)火山,大同火山活動(dòng)與山西地塹的演化發(fā)展密切相關(guān),是華北克拉通內(nèi)備受關(guān)注的構(gòu)造單元之一.前人對(duì)大同火山的研究大都針對(duì)其巖漿的深部來(lái)源.如Lei (2012)基于遠(yuǎn)震P波到時(shí)層析成像結(jié)果認(rèn)為大同火山的巖漿可能來(lái)源于下地幔深部上涌的地幔柱;Li S L等(2018)基于背景噪聲和遠(yuǎn)震面波層析成像結(jié)果認(rèn)為大同火山的巖漿可能來(lái)源于青藏高原東北緣下的軟流圈.然而大同火山巖漿在地殼內(nèi)的上涌通道以及巖漿分布范圍目前還不清楚.因此,開(kāi)展該區(qū)域高分辨率的背景噪聲成像研究對(duì)于揭示大同火山群下方巖漿活動(dòng)機(jī)制具有重要意義.

    地震各向異性是研究巖石圈變形和構(gòu)造單元演化過(guò)程的重要工具.前人利用遠(yuǎn)震SKS分裂、接收函數(shù)Ps轉(zhuǎn)換波、Pn波走時(shí)成像等方法對(duì)該區(qū)域地殼及上地幔各向異性進(jìn)行了大量的研究(常利軍等, 2011; Li et al., 2011; Lü, 2019; Zheng et al., 2019).常利軍等(2011)基于SKS分裂結(jié)果與GPS觀測(cè)的運(yùn)動(dòng)場(chǎng)認(rèn)為鄂爾多斯塊體周緣殼幔形變符合垂直連貫變形模式;Zheng等(2019)基于Ps接收函數(shù)各向異性結(jié)果與GPS觀測(cè)速度場(chǎng)結(jié)果和SKS分裂結(jié)果對(duì)比,認(rèn)為華北克拉通中部造山帶的西北部表現(xiàn)為殼幔耦合變形模式,而東部則表現(xiàn)為殼幔解耦的變形模式;高原等(2010)基于近震剪切波分裂結(jié)果與遠(yuǎn)震SKS分裂結(jié)果對(duì)比,提出華北地區(qū)地殼與上地幔不是簡(jiǎn)單的殼幔解耦型,也不是強(qiáng)耦合型.因此,對(duì)于華北克拉通區(qū)域殼幔變形的耦合情況目前仍然存在爭(zhēng)議.此外,由于SKS分裂得到的各向異性主要來(lái)自于上地幔,接收函數(shù)Ps轉(zhuǎn)換波的結(jié)果是整個(gè)地殼的平均結(jié)果,Pn波成像只反映上地幔頂部的各向異性,很難得到地殼尺度內(nèi)隨深度的變化的各向異性特征.而不同周期面波的方位各向異性可以反映不同深度的速度結(jié)構(gòu)方位各向異性(顧勤平等, 2020b),有著更好的垂向分辨率,因此利用面波層析成像獲得華北克拉通區(qū)域的速度結(jié)構(gòu)及其方位各向異性對(duì)于解釋該區(qū)域的構(gòu)造變形和殼幔耦合特征及其地球動(dòng)力學(xué)機(jī)制具有重要意義.Chen等(2015)利用Rayleigh波層析成像得到了華北克拉通東北部的各向異性結(jié)果,然而由于臺(tái)站覆蓋范圍有限,對(duì)整個(gè)華北克拉通中部造山帶的研究不足.本文基于華北克拉通中部306個(gè)寬頻帶流動(dòng)地震臺(tái)站的波形記錄,利用背景噪聲互相關(guān)技術(shù)獲得了研究區(qū)8~35 s的Rayleigh波相速度及其方位各向異性,得到了具有高分辨率的成像結(jié)果,進(jìn)而對(duì)山西地塹南北段不同的速度結(jié)構(gòu)及變形特征、大同火山群下地幔熱物質(zhì)上涌通道、區(qū)域內(nèi)的地殼與地幔耦合關(guān)系等問(wèn)題進(jìn)行了探討.

    1 數(shù)據(jù)和方法

    本文使用了中國(guó)地震科學(xué)臺(tái)陣探測(cè)項(xiàng)目三期一階段布設(shè)于華北克拉通中部的306個(gè)臺(tái)站從2017年1月1日至2017年7月28日觀測(cè)的垂向連續(xù)波形數(shù)據(jù),臺(tái)站統(tǒng)一應(yīng)用了型號(hào)為REFTEK數(shù)據(jù)采集器和CMG-3T地震計(jì)(頻帶范圍為50 Hz~120 s),臺(tái)間距~30 km,位置如圖1所示.臺(tái)站的分布主要覆蓋了華北克拉通中部造山帶及其鄰近區(qū)域,包括鄂爾多斯塊體的東半部分和華北盆地的西部局部地區(qū).

    圖1 研究區(qū)域構(gòu)造背景圖(a)和研究所用臺(tái)站分布圖(b)(a) 圖中黑實(shí)線表示華北克拉通構(gòu)造線(Zhao et al., 2005),紅色實(shí)線代表斷層分布,紅色三角形代表大同火山群; (b) 圖中藍(lán)色三角形代表本研究使用的臺(tái)站.Fig.1 Tectonic settings of the study area (a) and the location of seismic stations for this study (b)(a) The black solid lines in the figure represent the tectonic lines of North China Craton (Zhao et al., 2005), the red solid lines represent the locations of faults, and the red triangle represents the Datong Volcanos; (b) The blue triangles represent the stations used in this study.

    單臺(tái)數(shù)據(jù)預(yù)處理過(guò)程主要參照Bensen等(2007)提出的處理方法.首先將原始數(shù)據(jù)重采樣至1 Hz,然后進(jìn)行去儀器響應(yīng)、帶通濾波(5~50 s)、去線性趨勢(shì)和去均值處理,圖2展示了部分臺(tái)站一天的連續(xù)記錄經(jīng)過(guò)上述處理過(guò)后的波形,除了少數(shù)波形存在局部間斷外,絕大部分觀測(cè)數(shù)據(jù)質(zhì)量都比較好.接著采用滑動(dòng)絕對(duì)平均法對(duì)波形進(jìn)行時(shí)域歸一化,并進(jìn)行了頻域譜白化.

    完成數(shù)據(jù)預(yù)處理之后,對(duì)所有的臺(tái)站對(duì)每天的記錄進(jìn)行了波形互相關(guān)運(yùn)算,為了提高波形的信噪比,我們采用兩步法對(duì)每天的互相關(guān)波形進(jìn)行疊加.首先,把每天的互相關(guān)結(jié)果按時(shí)間排序并對(duì)每5天的記錄進(jìn)行線性疊加,然后利用Li G L等(2018)提出的時(shí)頻相位加權(quán)疊加法(tf-PWS)對(duì)5天疊加結(jié)果進(jìn)行進(jìn)一步疊加得到最終的臺(tái)站對(duì)間的互相關(guān)波形.相對(duì)于傳統(tǒng)的線性疊加,該方法被證明能夠有效地提升波形信噪比同時(shí)保持波形的頻散特性.其中,波形的信噪比是參照Bensen等(2007)進(jìn)行定義的,即信號(hào)窗內(nèi)的最大振幅與信號(hào)窗之后500 s的噪聲窗內(nèi)振幅的均方根之比.為了進(jìn)一步提高信噪比,我們將互相關(guān)波形的因果信號(hào)與果因信號(hào)進(jìn)行疊加得到最后的對(duì)稱分量進(jìn)行頻散曲線的提取,以臺(tái)站13801為例,圖3展示了最后的互相關(guān)對(duì)稱分量波形.

    圖2 經(jīng)過(guò)去儀器響應(yīng)、帶通濾波(5~50 s)、去線性趨勢(shì)和去均值處理的連續(xù)波形Fig.2 Continuous waveforms after removing instrument response, bandpass filtering (5~50 s), de-linear trend and de-mean

    圖3 以臺(tái)站13801為中心的互相關(guān)對(duì)稱分量波形Fig.3 Symmetric component of the cross-correlations between station 13801 and other stations

    圖4 華北克拉通中部Rayleigh波層析成像分辨率測(cè)試結(jié)果Fig.4 Resolution test results of Rayleigh wave tomography in the central North China Craton

    圖5 不同周期Rayleigh波相速度和方位各向異性分布圖(a) YM:陰山山脈;LM:呂梁山脈;TM:太行山脈;HG:河套地塹;OB:鄂爾多斯盆地;DB:大同盆地;TB:太原盆地;YB:運(yùn)城盆地;SR:山西地塹;WR:渭河地塹;NCB:華北盆地. (g)紅色五角星標(biāo)注了約35 km的深震位置.Fig.5 Maps for phase velocity and azimuthal anisotropy of Rayleigh waves at different periods(a) YM: Yinshan Mountains; LM: Luliang Mountains; TM: Taihang Mountains; HG: Hetao Graben; OB: Ordos Basin; DB: Datong Basin; TB: Taiyuan Basin; YB: Yuncheng Basin; SR: Shanxi Rift; WR: Weihe Rift; NCB: North China Basin. (g)The red star marks the location of the deep earthquake at ~35 km depth.

    圖6 (a) 研究區(qū)域的平均S波速度模型(紅色),藍(lán)色線條是AK135速度模型(Kennett et al., 1995); (b) 不同周期Rayleigh波相速度對(duì)S波速度的敏感核Fig.6 (a) The average S-wave velocity model (red line) of the study area. The blue line is the AK135 earth model (Kennett et al., 1995); (b) Sensitivity kernels to S-wave velocity for Rayleigh wave phase velocities at different periods

    本文采用Levshin和Ritzwoller(2001)發(fā)展的自動(dòng)頻時(shí)分析方法(FTAN)來(lái)獲取每個(gè)臺(tái)站對(duì)之間8~35 s的Rayleigh波相速度頻散曲線.為滿足遠(yuǎn)場(chǎng)近似條件,我們要求臺(tái)站對(duì)之間的距離大于三倍波長(zhǎng),同時(shí)為了進(jìn)一步提高Rayleigh波相速度的可靠性,只采用信噪比大于20的互相關(guān)結(jié)果.

    經(jīng)過(guò)對(duì)臺(tái)站對(duì)間的頻散曲線進(jìn)行初步篩選后,我們采用了Barmin等(2001)發(fā)展的一種基于線性射線理論的阻尼最小二乘法反演方法來(lái)反演二維Rayleigh波相速度及其方位各向異性.對(duì)于弱各向異性介質(zhì),面波相速度(或群速度)可以表示為(Smith and Dahlen, 1973)

    c(T)=c0(T)+c1(T)cos(2θ)+c2(T)sin(2θ)

    +c3(T)cos(4θ)+c4(T)sin(4θ),

    (1)

    其中T為面波周期,θ為方位角,c0為各向同性分量,c1,c2為2θ項(xiàng)各向異性分量,c3,c4為4θ項(xiàng)各向異性分量.考慮到4θ項(xiàng)通常可以被忽略(Montagner and Nataf, 1986; Debayle and Sambridge, 2004),本文只反演2θ項(xiàng),各向異性強(qiáng)度和快波方位角分別可由(2)、(3)式得到

    (2)

    (3)

    為了進(jìn)一步剔除頻散曲線測(cè)量中可能的異常值,我們采用了兩步反演法來(lái)獲取最終的二維Rayleigh波相速度模型.首先利用初步篩選的頻散結(jié)果進(jìn)行第一次反演,接著參考Guo等(2016)的標(biāo)準(zhǔn)對(duì)反演的走時(shí)殘差進(jìn)行評(píng)估,剔除了走時(shí)殘差大于6 s的測(cè)量值,然后將篩選過(guò)后的測(cè)量值再次進(jìn)行反演,最終得到了二維Rayleigh波相速度模型.在反演過(guò)程中,發(fā)現(xiàn)各向異性的絕對(duì)強(qiáng)度受到選擇的平滑參數(shù)的影響,而相對(duì)各向異性強(qiáng)度和快波方向基本不受影響,因此本文只對(duì)各向異性強(qiáng)度的相對(duì)變化和快波方向進(jìn)行討論.

    反演采用0.25°×0.25°的網(wǎng)格進(jìn)行插值,為檢測(cè)結(jié)果的可靠性,同時(shí)進(jìn)行了分辨率測(cè)試.圖4 展示了不同周期的分辨率測(cè)試結(jié)果,測(cè)試顯示研究區(qū)域內(nèi)的大部分地區(qū)都能得到小于25 km的分辨率.隨著周期增加,研究區(qū)域的邊緣部分分辨率有所下降,而研究區(qū)域中心的大部分都保持了高分辨率,保證了反演結(jié)果的可靠性.

    2 反演結(jié)果

    基于上述反演過(guò)程,我們得到了8~35 s的Rayleigh波各向同性相速度及其方位各向異性結(jié)果.圖5展示了不同周期的Rayleigh波相速度及其方位各向異性反演結(jié)果,黑色短棒代表方位各向異性強(qiáng)度和快波方向.

    為評(píng)估不同周期相速度對(duì)不同深度范圍速度結(jié)構(gòu)的反映,我們利用研究區(qū)域內(nèi)的平均S波速度模型計(jì)算了不同周期Rayleigh波相速度對(duì)S波速度的深度敏感核函數(shù)(Sensitivity Kernel),如圖6所示.本文將利用不同周期Rayleigh波相速度及其方位各向異性來(lái)討論對(duì)應(yīng)的研究區(qū)域上、中、下地殼及上地幔頂部的速度結(jié)構(gòu)和方位各向異性.8~12 s的相速度主要對(duì)上地殼的剪切波速度敏感,隨著周期增大,相速度的敏感深度范圍也逐漸增加,因此短周期相速度對(duì)剪切波速度結(jié)構(gòu)分辨率較高,而長(zhǎng)周期相速度的分辨率有所降低.16~20 s 的相速度主要反映中地殼的速度結(jié)構(gòu),而24~30 s 的相速度主要反映的是下地殼速度結(jié)構(gòu),35 s則主要反映上地幔頂部的速度結(jié)構(gòu).

    8~12 s的相速度異常主要與地表的結(jié)晶基底和沉積層的分布有關(guān).在華北盆地、鄂爾多斯盆地以及河套地塹、渭河地塹、山西地塹中的斷陷盆地等區(qū)域表現(xiàn)出明顯的低速異常,其中沿著渭河、山西地塹的斷陷盆地主要有運(yùn)城盆地、太原盆地以及大同盆地,區(qū)域內(nèi)低速異常反映了這些區(qū)域擁有較厚的沉積層.陰山山脈、呂梁山脈以及太行山脈等造山帶地區(qū)則表現(xiàn)出顯著的高速異常,這與地表的地質(zhì)構(gòu)造背景也基本吻合.該周期段的面波方位各向異性也主要與地表的地質(zhì)構(gòu)造背景有關(guān),在盆地和地塹與造山帶的交匯處,各向異性強(qiáng)度相對(duì)較大且快波方向主要沿著交匯邊界的走向分布,如華北盆地與太行造山帶的分界線附近,以及河套地塹與陰山造山帶的分界線附近,都表現(xiàn)出該特性,這與高原和吳晶(2008)基于直達(dá)剪切波分裂的結(jié)果在太行造山帶東緣表現(xiàn)出與邊界大致平行的快波方向相符合.而在構(gòu)造單元內(nèi)部,各向異性強(qiáng)度則相對(duì)較小且快波方向主要沿著構(gòu)造單元的走向分布,如在陰山造山帶內(nèi)部,快波方向以E-W向?yàn)橹鳎趨瘟?、太行造山帶?nèi)部,快波方向則以NNE-SSW方向?yàn)橹?,這與前人的面波各向異性研究的結(jié)果(Shen et al., 2016;胥鴻睿,2018)大致符合.值得注意的是鄂爾多斯塊體內(nèi)部大部分區(qū)域都表現(xiàn)出較弱的低速異常和較弱的方位各向異性,而其東北角鄰近大同火山群的部分地區(qū)則表現(xiàn)為極大的高速異常以及較強(qiáng)的NNE-SSW向各向異性.

    16~20 s 的相速度異常主要反映了研究區(qū)域中部地殼的速度結(jié)構(gòu),隨著周期增大,相速度異常表現(xiàn)出與較短周期不同的特性.如華北盆地、河套地塹、渭河地塹區(qū)域的低速異常逐漸減弱,陰山、呂梁、太行造山帶地區(qū)的高速異常也同樣減弱.而大同火山群附近的低速異常區(qū)域逐漸增大,其范圍擴(kuò)展至太行造山帶北段、陰山造山帶東部,向西與河套地塹低速異常相連接,向南延伸至山西地塹北段,且低速異常逐漸增強(qiáng),可能與大同火山群的巖漿活動(dòng)有關(guān).隨著周期增加,山西地塹南段的低速異常逐漸減小,轉(zhuǎn)變?yōu)楦咚佼惓?,即沿著山西地塹,其南段與北段表現(xiàn)出完全不同的速度異常,可能代表著沿著山西地塹在中地殼存在明顯的南北差異.還有一個(gè)值得注意的變化是華北盆地的低速異常隨著周期增大開(kāi)始向高速異常轉(zhuǎn)變,可能代表在該區(qū)域已經(jīng)逐漸接近下地殼.該周期段的方位各向異性強(qiáng)度相對(duì)短周期有所減小,但總體的快波方向基本保持不變,即沿山西地塹整體表現(xiàn)為NNE-SSW方向?yàn)橹?,而在陰山造山帶表現(xiàn)為E-W向?yàn)橹?,與前人基于直達(dá)剪切波分裂的結(jié)果(高原和吳晶,2008; 趙博等,2011)大致吻合.與Ps波接收函數(shù)的結(jié)果相比,陰山造山帶近E-W向的快波方向與楊妍等(2018)的結(jié)果一致,而大同盆地附近區(qū)域NNE-SSW向的快波方向與Zheng等(2019)的結(jié)果相符合.其中,楊妍等(2018)和Zheng等(2019)的結(jié)果在大同盆地區(qū)域存在差異,我們的結(jié)果與Zheng等(2019)的結(jié)果更為接近.在鄂爾多斯塊體內(nèi)部以及華北盆地方位各向異性強(qiáng)度趨近于0,在大同火山群附近表現(xiàn)出較弱的NNE-SSW向快波方向.值得注意的是在20 s 周期,山西地塹南段的快波方向逐漸開(kāi)始向NEE-SWW向轉(zhuǎn)變.

    24~30 s 的相速度異常主要代表了下地殼的速度結(jié)構(gòu)特性,在該周期段內(nèi),研究區(qū)域內(nèi)最顯著的低速異常主要集中在大同火山群及其鄰近的鄂爾多斯東北角并向西延伸至河套地塹.隨著周期增大,低速異常的范圍以及強(qiáng)度也逐漸變大,該低速異常的中心表現(xiàn)出向西遷移的特性,可能反映了巖漿上升的通道.在鄂爾多斯內(nèi)部整體依然表現(xiàn)為高速異常,代表其仍然保留著克拉通的特性.而山西地塹中南段則完全轉(zhuǎn)變?yōu)楦咚佼惓#c西邊的鄂爾多斯東南部和東邊的太行造山帶的高速異常連為一體,表現(xiàn)出與山西地塹北段完全不同的速度結(jié)構(gòu)特征.華北盆地則完全轉(zhuǎn)變?yōu)楦咚佼惓?,代表該區(qū)域可能已經(jīng)到達(dá)了上地幔頂部.該周期段方位各向異性的分布與較短周期也有所不同,在陰山造山帶和山西地塹最北端靠近大同火山群低速異常體的區(qū)域出現(xiàn)了圍繞該低速異常體的環(huán)狀快波方向分布,即在該低速異常體的西北邊緣表現(xiàn)為近NE-SW或NEE-SWW向,而東北邊緣表現(xiàn)為較弱的NWW-SEE向.另一個(gè)值得注意的變化是在山西地塹中段,即太原盆地附近,快波方向由NE-SW向轉(zhuǎn)變?yōu)镹EE-SWW向,逐漸與山西地塹南段表現(xiàn)一致.

    35 s的相速度異常主要反映了上地幔頂部的速度結(jié)構(gòu).在該周期段內(nèi),研究區(qū)域內(nèi)的低速異常分布與下地殼較為一致,大同火山群區(qū)域的低速異常體中心逐漸遷移至其西邊150 km處,低速異常的范圍基本保持不變.值得注意的變化是在該低速異常的南端出現(xiàn)了向東南方向擴(kuò)展的趨勢(shì),即太原盆地東邊的太行造山帶中部逐漸轉(zhuǎn)變?yōu)榈退佼惓?鄂爾多斯塊體、山西地塹中南段依舊表現(xiàn)為高速異常,而華北盆地的高速異常有所減弱.整體的方位各向異性強(qiáng)度相對(duì)于下地殼有所增強(qiáng),圍繞大同火山群附近的低速異常中心環(huán)狀的各向異性更為明顯,可能代表了上地幔頂部軟流圈物質(zhì)的上涌位置.另一個(gè)較為顯著的現(xiàn)象是在山西地塹中南段附近各向異性強(qiáng)度較大,且快波方向表現(xiàn)為近E-W向,這與山西地塹北段NNE-SSW向的快波方向有顯著差異.

    3 討論

    本研究得到的Rayleigh波相速度各向同性分量和方位各向異性分量反映了研究區(qū)域地殼至上地幔頂部50 km深度左右顯著的水平不均勻性和復(fù)雜的構(gòu)造變形特征.其中相速度主要對(duì)介質(zhì)的剪切波速度敏感,同時(shí)也與介質(zhì)的溫度、組成成分和介質(zhì)內(nèi)的流體性質(zhì)有關(guān),而造成相速度方位各向異性的因素卻非常復(fù)雜.地殼中的各向異性通常被歸因于裂隙的走向或礦物的定向排列(Rabbel and Mooney, 1996).其中在上地殼淺部,一般認(rèn)為地震波的方位各向異性可能是由區(qū)域地殼應(yīng)力場(chǎng)導(dǎo)致的地層裂隙和構(gòu)造單元的走向引起的(Crampin, 1994),因此上地殼的方位各向異性主要與地表的構(gòu)造特征有關(guān),一般表現(xiàn)為與斷層走向一致.而在更深處的下地殼和上地幔中的各向異性則一般認(rèn)為是形變導(dǎo)致的各向異性礦物定向排列造成的,其中下地殼和上地幔方位各向異性分別主要與角閃石和橄欖巖的定向排列有關(guān)(Liu et al., 2016).巖石圈地幔主要由橄欖石組成,其各向異性可以在構(gòu)造期內(nèi)長(zhǎng)期保存,因此可以用來(lái)研究構(gòu)造變形歷史以及與相鄰構(gòu)造塊體的相互作用過(guò)程,而軟流圈地幔的各向異性可能與地幔流的流動(dòng)方向有關(guān)(Long and Becker, 2010).考慮到本研究獲得的相速度最大周期可能無(wú)法有效反映軟流圈地幔的各向異性信息,我們將主要討論華北克拉通中部地殼及上地幔頂部的速度結(jié)構(gòu)和構(gòu)造變形特征.

    3.1 山西地塹

    本研究得到的Rayleigh波相速度以及方位各向異性結(jié)果都顯示出山西地塹南北段在中地殼至上地幔頂部存在顯著差異.在上地殼對(duì)應(yīng)周期段(8~12 s),相速度的結(jié)果顯示山西地塹南北段包括大同盆地、太原盆地、運(yùn)城盆地等區(qū)域都存在低速異常,推測(cè)可能與這些區(qū)域內(nèi)的新生代沉積層有關(guān)(Zhang et al., 1998).南北段的方位各向異性表現(xiàn)出相似特性,即主要沿NNE-SSW方向,與該區(qū)域斷層以及構(gòu)造單元的走向大致吻合,推測(cè)山西地塹上地殼的相速度方位各向異性主要來(lái)自于斷層裂隙和地表構(gòu)造單元走向.在中下地殼以及上地幔對(duì)應(yīng)周期段(16 s以上),山西地塹北段的低速異常逐漸增強(qiáng)且范圍逐漸增大,主要與該區(qū)域第四紀(jì)大同火山群的巖漿活動(dòng)有關(guān)(唐有彩等, 2011; Li S L et al., 2018).而在山西地塹的南段,隨著周期增加,低速異常逐漸消失,取而代之的是與西側(cè)鄂爾多斯塊體和東側(cè)太行造山帶一致的高速異常,表明山西地塹的巖漿活動(dòng)目前只在北段較為活躍(Bao et al., 2013).在該周期段內(nèi),相速度方位各向異性也表現(xiàn)出不同特征,北段的各向異性強(qiáng)度較弱,且主要沿NNE-SSW或NE-SW向,而中南段的各向異性強(qiáng)度較大且快波方向逐漸由NNE-SSW向順時(shí)針旋轉(zhuǎn)至E-W向.南北段的各向異性強(qiáng)度差異可能主要受到巖漿活動(dòng)的影響,即北段由于地幔熱物質(zhì)上涌使得下地殼以及上地幔頂部的礦物排列方向受到影響導(dǎo)致方位各向異性強(qiáng)度相對(duì)較弱.Li等(2015)利用山西地塹的震源機(jī)制解得到了該區(qū)域的應(yīng)力場(chǎng)分布,結(jié)果顯示山西地塹的最大壓應(yīng)力方向在北段主要呈NE-SW向,而在南段則主要呈NEE-SWW向,這與我們得到的地殼對(duì)應(yīng)周期的Rayleigh波快波方向大致吻合,同時(shí)也與Qu等(2014)基于GPS觀測(cè)結(jié)果得到的區(qū)域應(yīng)變場(chǎng)的方向比較對(duì)應(yīng).Ai等(2020)結(jié)合臨汾盆地附近區(qū)域的震源深度達(dá)到35 km以及表現(xiàn)出較高的速度異常,認(rèn)為在該區(qū)域中下地殼表現(xiàn)出較強(qiáng)的脆性變形,這與本文結(jié)果顯示該區(qū)域地殼內(nèi)快波方向與最大壓應(yīng)力方向大致平行比較符合.因此我們認(rèn)為華北克拉通中部地殼內(nèi)的面波方位各向異性主要受應(yīng)力場(chǎng)的影響,但北段還受到巖漿活動(dòng)影響.結(jié)合Li等(2015)應(yīng)力場(chǎng)的研究結(jié)果,山西地塹北段主要表現(xiàn)為NW-SE向的拉張,南段則表現(xiàn)為NNW-SSE向的拉張,可能與GPS觀測(cè)到的鄂爾多斯塊體逆時(shí)針旋轉(zhuǎn)有關(guān)(李延興等, 2005),即青藏高原隆升對(duì)鄂爾多斯西南角的持續(xù)東北向擠壓導(dǎo)致了其逆時(shí)針旋轉(zhuǎn),進(jìn)而造成山西地塹南段形成了NNW-SSE方向的簡(jiǎn)單剪切,而在北段則形成NW-SE方向的純剪切.在上地幔頂部對(duì)應(yīng)的周期(35 s),山西地塹中南段的快波方向呈E-W向,這與常利軍等(2021)得到的SKS分裂快波方向比較吻合,可能與地幔物質(zhì)流動(dòng)有關(guān).綜合以上討論,我們的結(jié)果傾向于支持山西地塹的演化發(fā)展主要受印度—亞洲板塊碰撞的遠(yuǎn)場(chǎng)效應(yīng)影響.

    3.2 大同火山群

    大同火山群區(qū)域的低速異常是本研究中非常突出的一個(gè)相速度異常特征,該低速異常體在前人的研究觀測(cè)中也可以被觀察到(Tang et al., 2013; Chen et al., 2015; Li S L et al., 2018).由于本研究使用的臺(tái)站覆蓋更為密集,因此結(jié)果分辨率更高,可以得到地殼及上地幔頂部更精細(xì)的結(jié)果.本文結(jié)果顯示該低速異常體沿深度方向可以被連續(xù)觀測(cè)到且其特性表現(xiàn)出很強(qiáng)的垂向變化,包括低速異常強(qiáng)度、范圍、位置等.在上地殼對(duì)應(yīng)周期(8~12 s),該區(qū)域低速異常較弱,且范圍較小,中心位于大同盆地,推測(cè)應(yīng)該是由盆地淺部的低速沉積層造成.其方位各向異性較強(qiáng)且主要沿NEE-SWW和NE-SW向,與該區(qū)域斷層走向一致,推測(cè)是由地表斷層構(gòu)造產(chǎn)生的.在中下地殼以及上地幔頂部對(duì)應(yīng)周期段(16 s以上),其低速異常強(qiáng)度隨深度(周期)有所增加,且范圍逐漸增大,低速體的中心逐漸向西偏移至150 km處.中下地殼以及上地幔頂部的低速異常主要與大同火山群的巖漿活動(dòng)有關(guān),該結(jié)果也和前人觀測(cè)到的新生代地幔源巖漿巖出露相一致(Xu et al., 2005),Zhang等(2016)利用大地電磁觀測(cè)到大同火山群下存在的低阻體同樣反映了地殼中的巖漿活動(dòng).該低速異常隨深度表現(xiàn)出強(qiáng)度增加以及范圍擴(kuò)大可能代表了自下而上的熱物質(zhì)上涌對(duì)地殼的加熱甚至可能導(dǎo)致了部分熔融.而其中心向西偏移可能代表了地幔熱物質(zhì)上涌的通道,即地幔熱物質(zhì)從鄂爾多斯東北角之下的軟流圈經(jīng)鄂爾多斯東北角與中部造山帶之間向西傾斜的剪切帶上涌,形成了第四紀(jì)大同火山群的玄武質(zhì)巖漿噴發(fā),如圖7所示.

    圖7 大同火山群下方巖漿上涌示意圖Fig.7 Schematic diagram of magma upwelling under DatongVolcanos

    下地殼至上地幔頂部對(duì)應(yīng)周期(24 s以上)的相速度方位各向異性圍繞該低速異常體的環(huán)狀快波方向大致勾勒出地幔熱物質(zhì)在不同深度的侵入位置和范圍,與相速度低速異常相吻合.這種環(huán)狀快波方向分布在常利軍等(2021)在該區(qū)域利用SKS分裂得到的結(jié)果中也被觀測(cè)到,證實(shí)該環(huán)狀分布的快波方向延伸到了上地幔.低速異常體中心的各向異性較弱,而邊緣的各向異性較強(qiáng),類似的特性也被前人在該區(qū)域的面波各向異性研究觀察到(Chen et al., 2015),推測(cè)可能是地幔熱物質(zhì)上涌導(dǎo)致了礦物排列方向的垂向排列,形成了徑向各向異性,而方位各向異性被減弱(Iyer and Hirahara, 1993; Chen et al., 2015).其邊緣較強(qiáng)的環(huán)狀快波方向則可能是由熱物質(zhì)上涌對(duì)地殼的加熱甚至部分熔融導(dǎo)致該低速體邊緣較高的速度梯度所造成.因此,本研究結(jié)合面波速度異常及其方位各向異性對(duì)大同火山群下地幔熱物質(zhì)在地殼中上涌的通道有了很好的限制,而熱物質(zhì)更深處的來(lái)源還有待進(jìn)一步的研究.

    3.3 殼幔耦合分析

    殼幔之間的耦合關(guān)系與巖石圈的歷史動(dòng)力學(xué)過(guò)程密切相關(guān),同時(shí)也受到現(xiàn)今板塊運(yùn)動(dòng)的巨大影響.華北克拉通的殼幔耦合關(guān)系對(duì)于研究華北地區(qū)深部動(dòng)力學(xué)機(jī)制和運(yùn)動(dòng)模型具有非常重要的意義.本研究使用背景噪聲互相關(guān)得到的相速度及方位各向異性的周期范圍反映的速度結(jié)構(gòu)及變形特征限于地殼以及上地幔頂部,而一般認(rèn)為SKS分裂得到的各向異性主要來(lái)自上地幔(常利軍等, 2008, 2011),因此結(jié)合本研究得到的方位各向異性以及前人的SKS分裂結(jié)果,可以對(duì)研究區(qū)域的地殼與地幔運(yùn)動(dòng)變形的耦合問(wèn)題進(jìn)行探討.圖8顯示了本研究得到的下地殼對(duì)應(yīng)周期的Rayleigh波相速度快波方向與常利軍等(2011, 2012, 2021)得到的SKS分裂結(jié)果的對(duì)比.可以觀察到相速度快波方向與SKS分裂快波方向在局部區(qū)域存在著顯著的差異.如在山西地塹北段大同火山群鄰近區(qū)域,下地殼的面波快波方向主要表現(xiàn)為NE-SW和NNE-SSW方向,該結(jié)果與Zheng等(2019)利用Ps波接收函數(shù)和張暉等(2020)利用直達(dá)剪切波分裂在該區(qū)域得到的地殼快波方向比較吻合,而SKS分裂的快波方向則表現(xiàn)出與相速度快波方向近乎垂直的NW-SE和NWW-SEE方向,這與Chen等(2015)利用面波各向異性得到的該區(qū)域在30 km和80 km快波方向不同比較對(duì)應(yīng).殼幔不一致的快波方向似乎表現(xiàn)出殼幔解耦的變形特征,然而考慮到該區(qū)域受巖漿活動(dòng)影響劇烈,面波的各向異性來(lái)源非常復(fù)雜,可能無(wú)法準(zhǔn)確反映地殼變形特性(Yao et al., 2010).常利軍等(2011)基于GPS觀測(cè)的地表運(yùn)動(dòng)速度場(chǎng)與SKS分裂快波方向推測(cè)鄂爾多斯周緣的殼幔變形符合垂直連貫變形模式,而高原等(2010)根據(jù)華北克拉通東北部近場(chǎng)資料得到的剪切波分裂快波方向與遠(yuǎn)震SKS分裂快波方向存在差異,提出華北地區(qū)地殼與上地幔不是簡(jiǎn)單的殼幔解耦型,也不是強(qiáng)耦合型.考慮到GPS一般被認(rèn)為可以代表上地殼的形變特征(魯來(lái)玉等, 2014),而中下地殼變形模式可能非常復(fù)雜,同時(shí)該區(qū)域還受到巖漿活動(dòng)的影響,因此該區(qū)域的殼幔耦合情況還有待進(jìn)一步的研究.而在山西地塹中南段及其鄰近區(qū)域,下地殼以及上地幔頂部的相速度表現(xiàn)為高速異常,其快波方向主要表現(xiàn)為NEE-SWW和E-W向,這與常利軍等(2011)SKS分裂在該區(qū)域表現(xiàn)出沿E-W、NEE-SWW和NWW-SEE的快波方向大致吻合,因此推測(cè)在山西地塹中南段地殼與地幔的運(yùn)動(dòng)變形更符合耦合特性,即該區(qū)域表現(xiàn)出殼幔的垂直連貫變形特征.考慮到地殼和地幔中的方位各向異性的來(lái)源非常復(fù)雜,可能受到多種因素的影響,因此該區(qū)域準(zhǔn)確的殼幔變形模式還需要利用不同資料和方法進(jìn)一步進(jìn)行探測(cè).

    圖8 Rayleigh波30 s和35 s的方位各向異性快波方向(黑色短棒)與SKS分裂(常利軍等, 2021)快波方向(紅色短棒)對(duì)比圖Fig.8 Comparison of azimuthal anisotropic fast wave directions from Rayleigh wave at 30 s and 35 s (black bars) and SKS splitting (red bars) (Chang et al., 2021)

    4 結(jié)論

    本文利用背景噪聲層析成像獲得了華北克拉通中部區(qū)域8~35 s的Rayleigh波相速度及其方位各向異性分布.結(jié)果顯示,短周期面波在山西地塹南北段都表現(xiàn)出低速異常,可能與山西地塹內(nèi)斷陷盆地中的新生代沉積層有關(guān),其方位各向異性的快波方向主要與區(qū)域內(nèi)斷層走向以及構(gòu)造單元走向一致.中長(zhǎng)周期面波代表的山西地塹南北段的速度結(jié)構(gòu)以及方位各向異性都存在較大的差異.隨著周期增大,北段的大同火山群低速異常逐漸增強(qiáng)且范圍逐漸變大,其中心逐漸向西遷移至大同火山群西部150 km處,可能代表了地幔熱物質(zhì)在地殼內(nèi)沿鄂爾多斯東北角與中部造山帶之間向西傾斜的剪切帶上涌的通道,而中南段的低速異常逐漸消失,取而代之的是垂向上連續(xù)而穩(wěn)定的高速異常,一直延伸到上地幔頂部,可能代表著目前山西地塹的巖漿活動(dòng)只在北段比較活躍.南北段方位各向異性強(qiáng)度的差異可能主要與北段大同火山群第四紀(jì)巖漿活動(dòng)有關(guān),北段受地幔熱物質(zhì)上涌影響導(dǎo)致方位各向異性減小,快波方向圍繞上涌的熱物質(zhì)呈環(huán)狀分布,與低速異常表現(xiàn)的地幔熱物質(zhì)上涌通道吻合.地殼對(duì)應(yīng)周期的面波快波方向與區(qū)域最大主壓應(yīng)力方向大致符合,推測(cè)地殼內(nèi)的方位各向異性主要受地殼應(yīng)力場(chǎng)的影響,但北段還受到巖漿活動(dòng)影響.中南段的下地殼至上地幔頂部表現(xiàn)為高速異常,其快波方向逐漸與前人SKS分裂結(jié)果一致,可能代表了該區(qū)域的殼幔耦合變形模式.而北段的下地殼和上地幔頂部的快波方向與SKS分裂結(jié)果不一致,似乎表現(xiàn)出殼幔解耦特征,但考慮到北段受到巖漿活動(dòng)影響,其殼幔變形模型有待進(jìn)一步研究.

    致謝感謝中國(guó)地震局地球物理研究所“中國(guó)地震科學(xué)探測(cè)臺(tái)陣數(shù)據(jù)中心”為本研究提供地震波形數(shù)據(jù).感謝Rice University的Li G L博士為我們提供了數(shù)據(jù)處理程序.

    猜你喜歡
    克拉通華北方位
    巖石圈地幔分層性對(duì)克拉通穩(wěn)定性的影響
    認(rèn)方位
    幼兒園(2021年12期)2021-11-06 05:10:20
    華北玉米市場(chǎng)將進(jìn)入筑底期
    有關(guān)克拉通破壞及其成因的綜述
    華北克拉通重力剖面重力點(diǎn)位GPS測(cè)量精度分析
    Literature Review on Context Translation Mode
    借助方位法的拆字
    說(shuō)方位
    幼兒100(2016年28期)2016-02-28 21:26:17
    拉張槽對(duì)四川盆地海相油氣分布的控制作用
    基于TMS320C6678的SAR方位向預(yù)濾波器的并行實(shí)現(xiàn)
    av在线蜜桃| 久久久久国内视频| 精品国产超薄肉色丝袜足j| 精品久久久久久久久久免费视频| 美女 人体艺术 gogo| 久久久久九九精品影院| 黄频高清免费视频| 亚洲av熟女| 99国产精品一区二区蜜桃av| 天天躁狠狠躁夜夜躁狠狠躁| 老司机在亚洲福利影院| 亚洲成av人片在线播放无| 一区二区三区国产精品乱码| 欧美大码av| 中文字幕最新亚洲高清| 国产v大片淫在线免费观看| 国产成人av教育| 中文字幕熟女人妻在线| 中文字幕人妻丝袜一区二区| 国产精品一及| 精品不卡国产一区二区三区| 香蕉久久夜色| 久久久久久久精品吃奶| 男女之事视频高清在线观看| 综合色av麻豆| 免费电影在线观看免费观看| 久久亚洲精品不卡| 美女大奶头视频| 丁香欧美五月| 在线视频色国产色| 日韩欧美精品v在线| ponron亚洲| 日韩中文字幕欧美一区二区| 两性午夜刺激爽爽歪歪视频在线观看| 国产激情偷乱视频一区二区| 亚洲国产中文字幕在线视频| 特大巨黑吊av在线直播| 国产精品美女特级片免费视频播放器 | 国产高清videossex| 国产成人av教育| 一边摸一边抽搐一进一小说| 亚洲国产精品成人综合色| 亚洲九九香蕉| 88av欧美| 国产私拍福利视频在线观看| 在线观看一区二区三区| 国产野战对白在线观看| 国产男靠女视频免费网站| 国产精品亚洲一级av第二区| 不卡一级毛片| 国产三级中文精品| 99国产极品粉嫩在线观看| 免费看a级黄色片| 99re在线观看精品视频| 村上凉子中文字幕在线| 婷婷精品国产亚洲av在线| 国产精品女同一区二区软件 | 级片在线观看| 熟女人妻精品中文字幕| 美女 人体艺术 gogo| 久久久色成人| 国产成人精品久久二区二区91| 夜夜爽天天搞| 禁无遮挡网站| 欧美日韩一级在线毛片| 国产美女午夜福利| 热99re8久久精品国产| 国产亚洲欧美98| 99精品久久久久人妻精品| 在线观看日韩欧美| 免费在线观看亚洲国产| 美女cb高潮喷水在线观看 | 麻豆一二三区av精品| 国产激情欧美一区二区| 9191精品国产免费久久| 熟妇人妻久久中文字幕3abv| 老司机在亚洲福利影院| 一区二区三区激情视频| 九九在线视频观看精品| 日本免费一区二区三区高清不卡| 国产高潮美女av| 亚洲专区中文字幕在线| 亚洲国产精品999在线| 日韩精品青青久久久久久| 99久久99久久久精品蜜桃| 亚洲最大成人中文| 757午夜福利合集在线观看| 亚洲熟女毛片儿| 99久久久亚洲精品蜜臀av| 国产av在哪里看| 天天躁狠狠躁夜夜躁狠狠躁| 日本五十路高清| 久久婷婷人人爽人人干人人爱| 国产av不卡久久| 美女黄网站色视频| 久久久国产精品麻豆| 国产av不卡久久| 天天添夜夜摸| 免费搜索国产男女视频| 中国美女看黄片| 给我免费播放毛片高清在线观看| 一进一出抽搐动态| 美女高潮的动态| 欧洲精品卡2卡3卡4卡5卡区| 欧美中文综合在线视频| 国产探花在线观看一区二区| 久久久久免费精品人妻一区二区| 老司机在亚洲福利影院| 1000部很黄的大片| 村上凉子中文字幕在线| 99re在线观看精品视频| 久久精品国产清高在天天线| 亚洲av日韩精品久久久久久密| 成人午夜高清在线视频| 亚洲成av人片免费观看| 最新在线观看一区二区三区| 99久久综合精品五月天人人| www.熟女人妻精品国产| av黄色大香蕉| 中文字幕高清在线视频| 搡老熟女国产l中国老女人| 国产精品一区二区精品视频观看| 日本黄色视频三级网站网址| 国产极品精品免费视频能看的| 久久久国产欧美日韩av| 午夜福利18| 国产激情久久老熟女| 日本撒尿小便嘘嘘汇集6| 亚洲国产中文字幕在线视频| 真人做人爱边吃奶动态| 两性午夜刺激爽爽歪歪视频在线观看| 在线观看免费视频日本深夜| www.自偷自拍.com| 国内揄拍国产精品人妻在线| 变态另类丝袜制服| 亚洲精品美女久久av网站| 久久久久久大精品| 两个人视频免费观看高清| 国产淫片久久久久久久久 | 观看免费一级毛片| 国产 一区 欧美 日韩| 一级毛片高清免费大全| 曰老女人黄片| 亚洲熟女毛片儿| 欧美又色又爽又黄视频| 亚洲成人精品中文字幕电影| 啦啦啦免费观看视频1| 黄色女人牲交| 国产精品自产拍在线观看55亚洲| 最近最新免费中文字幕在线| www.www免费av| 午夜激情福利司机影院| 国产v大片淫在线免费观看| 看免费av毛片| 亚洲自偷自拍图片 自拍| 欧美不卡视频在线免费观看| 男人舔女人的私密视频| 国产高潮美女av| 18禁黄网站禁片免费观看直播| 精品乱码久久久久久99久播| 国产精品精品国产色婷婷| 久久九九热精品免费| 欧美一级毛片孕妇| 亚洲欧美日韩高清在线视频| 老司机在亚洲福利影院| 精品一区二区三区视频在线观看免费| 好男人电影高清在线观看| 一a级毛片在线观看| 久久精品国产亚洲av香蕉五月| 在线观看免费视频日本深夜| 男人舔奶头视频| 午夜视频精品福利| 搡老岳熟女国产| 亚洲精品美女久久久久99蜜臀| 国产精品日韩av在线免费观看| 日本精品一区二区三区蜜桃| 午夜亚洲福利在线播放| 老鸭窝网址在线观看| 俺也久久电影网| 国产日本99.免费观看| 国内久久婷婷六月综合欲色啪| 亚洲成av人片在线播放无| 国内精品美女久久久久久| 曰老女人黄片| 亚洲熟女毛片儿| 久久亚洲精品不卡| 久久中文字幕一级| 亚洲国产欧美一区二区综合| 别揉我奶头~嗯~啊~动态视频| 色av中文字幕| 真人做人爱边吃奶动态| 中文亚洲av片在线观看爽| 国产男靠女视频免费网站| 亚洲精品色激情综合| av视频在线观看入口| 亚洲欧美激情综合另类| 日韩欧美国产在线观看| 国产毛片a区久久久久| 一级a爱片免费观看的视频| 97超视频在线观看视频| 日日干狠狠操夜夜爽| 国产欧美日韩一区二区精品| 亚洲精品456在线播放app | 90打野战视频偷拍视频| 午夜亚洲福利在线播放| 制服丝袜大香蕉在线| 国产高清有码在线观看视频| 免费看美女性在线毛片视频| 高清毛片免费观看视频网站| 欧美另类亚洲清纯唯美| 日韩大尺度精品在线看网址| 色老头精品视频在线观看| 国产av在哪里看| 久久久水蜜桃国产精品网| 女人高潮潮喷娇喘18禁视频| 老熟妇仑乱视频hdxx| 巨乳人妻的诱惑在线观看| 亚洲精品美女久久av网站| 最近最新中文字幕大全免费视频| 操出白浆在线播放| 国产aⅴ精品一区二区三区波| av国产免费在线观看| 国产高清视频在线播放一区| 国产精品久久久av美女十八| 黄色成人免费大全| 亚洲av中文字字幕乱码综合| 夜夜看夜夜爽夜夜摸| 久久精品夜夜夜夜夜久久蜜豆| 欧美日韩乱码在线| 免费在线观看成人毛片| 免费在线观看日本一区| 欧美日韩综合久久久久久 | 97超级碰碰碰精品色视频在线观看| 国产精品永久免费网站| 搞女人的毛片| 每晚都被弄得嗷嗷叫到高潮| 一进一出好大好爽视频| 国内毛片毛片毛片毛片毛片| 亚洲人成伊人成综合网2020| 一级毛片精品| 在线a可以看的网站| 12—13女人毛片做爰片一| 国产免费av片在线观看野外av| 欧美日韩乱码在线| 欧美丝袜亚洲另类 | 国产精品 欧美亚洲| 色综合婷婷激情| 两个人的视频大全免费| 久久婷婷人人爽人人干人人爱| 性色av乱码一区二区三区2| 国产麻豆成人av免费视频| 精品一区二区三区av网在线观看| 色播亚洲综合网| 熟女人妻精品中文字幕| 真实男女啪啪啪动态图| 日本熟妇午夜| 国产野战对白在线观看| or卡值多少钱| 中文字幕精品亚洲无线码一区| 中文字幕av在线有码专区| 国产精品99久久99久久久不卡| 日韩欧美 国产精品| www.999成人在线观看| 美女cb高潮喷水在线观看 | 国产av在哪里看| or卡值多少钱| 国产成人一区二区三区免费视频网站| 久久人妻av系列| 日本熟妇午夜| 一级a爱片免费观看的视频| 在线观看美女被高潮喷水网站 | 亚洲av中文字字幕乱码综合| 亚洲专区国产一区二区| 少妇丰满av| 亚洲欧美日韩卡通动漫| 婷婷丁香在线五月| 无限看片的www在线观看| 婷婷精品国产亚洲av在线| 国产精品久久久久久精品电影| 久久久久久大精品| 精品一区二区三区视频在线 | 欧洲精品卡2卡3卡4卡5卡区| 好看av亚洲va欧美ⅴa在| 国产精品综合久久久久久久免费| 国产精品,欧美在线| 岛国在线免费视频观看| 很黄的视频免费| 国产伦精品一区二区三区四那| 精品一区二区三区视频在线观看免费| 欧美黄色片欧美黄色片| 一本久久中文字幕| 欧美中文日本在线观看视频| 亚洲av片天天在线观看| 欧美又色又爽又黄视频| 欧美黑人欧美精品刺激| 久久亚洲真实| 热99re8久久精品国产| 免费在线观看视频国产中文字幕亚洲| 色哟哟哟哟哟哟| 亚洲激情在线av| 欧美在线一区亚洲| 成人亚洲精品av一区二区| 无限看片的www在线观看| 我的老师免费观看完整版| 日韩欧美一区二区三区在线观看| 亚洲欧美日韩卡通动漫| 女警被强在线播放| 亚洲人成网站高清观看| 亚洲精品在线美女| 99精品欧美一区二区三区四区| 这个男人来自地球电影免费观看| 亚洲国产精品久久男人天堂| 天天一区二区日本电影三级| 亚洲成人精品中文字幕电影| 极品教师在线免费播放| 国产视频内射| 国产激情偷乱视频一区二区| 日本黄色视频三级网站网址| 午夜免费观看网址| 国产成人影院久久av| 中文在线观看免费www的网站| 三级男女做爰猛烈吃奶摸视频| 色吧在线观看| 国模一区二区三区四区视频 | 色哟哟哟哟哟哟| 久久久水蜜桃国产精品网| 免费搜索国产男女视频| 久久精品国产99精品国产亚洲性色| 麻豆一二三区av精品| 日韩大尺度精品在线看网址| 一边摸一边抽搐一进一小说| 久久精品夜夜夜夜夜久久蜜豆| 欧美绝顶高潮抽搐喷水| 老司机午夜福利在线观看视频| 男人和女人高潮做爰伦理| 国产激情久久老熟女| 黄色丝袜av网址大全| 观看免费一级毛片| 曰老女人黄片| 色噜噜av男人的天堂激情| 亚洲精华国产精华精| 国产aⅴ精品一区二区三区波| 精品不卡国产一区二区三区| 免费av毛片视频| 女警被强在线播放| 亚洲无线观看免费| 精品久久久久久,| 又大又爽又粗| 日本一本二区三区精品| 精品国产亚洲在线| 国产精品日韩av在线免费观看| 精品福利观看| 1000部很黄的大片| 琪琪午夜伦伦电影理论片6080| 久久久久性生活片| 国产激情偷乱视频一区二区| 国产av在哪里看| 亚洲av熟女| 国产又黄又爽又无遮挡在线| 亚洲人成网站在线播放欧美日韩| 丰满的人妻完整版| 又大又爽又粗| 亚洲av中文字字幕乱码综合| 人人妻,人人澡人人爽秒播| 两个人视频免费观看高清| 精品久久久久久成人av| 欧美日韩综合久久久久久 | 亚洲 欧美 日韩 在线 免费| 天堂动漫精品| www国产在线视频色| 欧美日韩综合久久久久久 | 国产高清视频在线播放一区| 99国产综合亚洲精品| 午夜福利成人在线免费观看| 国产免费男女视频| av欧美777| 毛片女人毛片| 欧美大码av| 久久精品人妻少妇| 91在线精品国自产拍蜜月 | 日韩免费av在线播放| 黄色视频,在线免费观看| 精品午夜福利视频在线观看一区| 午夜福利在线观看免费完整高清在 | 国产黄片美女视频| 男女做爰动态图高潮gif福利片| 欧美色视频一区免费| 天天一区二区日本电影三级| 在线永久观看黄色视频| 日本 欧美在线| 成人欧美大片| 国产野战对白在线观看| 色播亚洲综合网| 国产探花在线观看一区二区| 精品久久久久久,| 91久久精品国产一区二区成人 | 一本一本综合久久| 美女被艹到高潮喷水动态| 亚洲国产看品久久| 99在线视频只有这里精品首页| 亚洲国产日韩欧美精品在线观看 | 91av网一区二区| 熟女人妻精品中文字幕| 亚洲av片天天在线观看| 国产不卡一卡二| 麻豆成人午夜福利视频| 欧美精品啪啪一区二区三区| 久久精品影院6| 国产精品永久免费网站| 男人舔奶头视频| 欧美xxxx黑人xx丫x性爽| 亚洲中文日韩欧美视频| 色噜噜av男人的天堂激情| 别揉我奶头~嗯~啊~动态视频| 日本熟妇午夜| 久久精品影院6| 性欧美人与动物交配| 日本成人三级电影网站| 少妇裸体淫交视频免费看高清| 丰满人妻熟妇乱又伦精品不卡| 我的老师免费观看完整版| 久久香蕉国产精品| 蜜桃久久精品国产亚洲av| 97碰自拍视频| 亚洲欧美日韩卡通动漫| 99久久99久久久精品蜜桃| 88av欧美| 免费看日本二区| 日本黄色视频三级网站网址| 亚洲,欧美精品.| 精品一区二区三区四区五区乱码| 国产精品亚洲一级av第二区| 欧美性猛交黑人性爽| 亚洲avbb在线观看| 女同久久另类99精品国产91| 听说在线观看完整版免费高清| 日本三级黄在线观看| 午夜福利在线观看免费完整高清在 | 国产一区二区在线av高清观看| 亚洲成a人片在线一区二区| 欧美日韩黄片免| 一夜夜www| 看片在线看免费视频| 午夜a级毛片| 老鸭窝网址在线观看| 欧美一级a爱片免费观看看| 亚洲激情在线av| 午夜福利视频1000在线观看| 免费在线观看视频国产中文字幕亚洲| 老汉色av国产亚洲站长工具| 99久久成人亚洲精品观看| 国产高潮美女av| 亚洲五月婷婷丁香| 国产精品乱码一区二三区的特点| 美女黄网站色视频| 欧美一级毛片孕妇| 精品一区二区三区四区五区乱码| 禁无遮挡网站| 亚洲最大成人中文| 精品久久久久久久久久久久久| 日韩中文字幕欧美一区二区| 一夜夜www| 日韩欧美在线乱码| 欧美午夜高清在线| 欧美色欧美亚洲另类二区| 身体一侧抽搐| 夜夜爽天天搞| 人人妻人人看人人澡| 国产精品免费一区二区三区在线| 日本黄色片子视频| 国产一区二区三区视频了| 啪啪无遮挡十八禁网站| 中文资源天堂在线| 亚洲精品色激情综合| svipshipincom国产片| 亚洲成人精品中文字幕电影| 99riav亚洲国产免费| 成年版毛片免费区| 成在线人永久免费视频| 人妻丰满熟妇av一区二区三区| 草草在线视频免费看| 久久精品国产99精品国产亚洲性色| 国产高清视频在线播放一区| 久久亚洲真实| 欧美乱色亚洲激情| 19禁男女啪啪无遮挡网站| 男女那种视频在线观看| 少妇人妻一区二区三区视频| 日韩有码中文字幕| 1000部很黄的大片| 一个人观看的视频www高清免费观看 | 午夜成年电影在线免费观看| 人人妻,人人澡人人爽秒播| 久久午夜亚洲精品久久| 网址你懂的国产日韩在线| 极品教师在线免费播放| 一级毛片高清免费大全| 一二三四社区在线视频社区8| 成人18禁在线播放| 日韩成人在线观看一区二区三区| 亚洲人与动物交配视频| 动漫黄色视频在线观看| 午夜福利在线观看吧| 99国产精品一区二区三区| 舔av片在线| 精品国内亚洲2022精品成人| 亚洲av美国av| 国产精品久久电影中文字幕| 色在线成人网| 91av网站免费观看| 亚洲中文日韩欧美视频| 99在线人妻在线中文字幕| 性欧美人与动物交配| 伊人久久大香线蕉亚洲五| 1024手机看黄色片| 在线观看免费视频日本深夜| 欧美在线一区亚洲| 日韩精品中文字幕看吧| 日本免费a在线| 日韩欧美在线乱码| 亚洲av电影在线进入| 啦啦啦韩国在线观看视频| 免费电影在线观看免费观看| www.熟女人妻精品国产| 色噜噜av男人的天堂激情| 99久久精品国产亚洲精品| 成人国产一区最新在线观看| 成人av在线播放网站| 亚洲 欧美 日韩 在线 免费| a级毛片在线看网站| 欧美日韩瑟瑟在线播放| 成在线人永久免费视频| 一区二区三区国产精品乱码| 88av欧美| 久久人妻av系列| 久久久成人免费电影| 国产成人欧美在线观看| 亚洲国产精品sss在线观看| 美女 人体艺术 gogo| 国产高潮美女av| 欧美日韩福利视频一区二区| 三级国产精品欧美在线观看 | 国产乱人伦免费视频| 国产伦人伦偷精品视频| 中文资源天堂在线| 舔av片在线| 少妇的逼水好多| 天堂动漫精品| xxx96com| 日韩欧美三级三区| 国产人伦9x9x在线观看| 国产一区在线观看成人免费| 精品久久蜜臀av无| 91av网一区二区| 少妇的丰满在线观看| 亚洲美女黄片视频| 成人鲁丝片一二三区免费| 女人高潮潮喷娇喘18禁视频| 日本撒尿小便嘘嘘汇集6| 久久久久久人人人人人| 免费观看的影片在线观看| 一二三四在线观看免费中文在| 色播亚洲综合网| 老司机在亚洲福利影院| 欧美黑人巨大hd| 白带黄色成豆腐渣| 国产黄色小视频在线观看| 在线国产一区二区在线| 女生性感内裤真人,穿戴方法视频| 国产探花在线观看一区二区| 成年女人永久免费观看视频| 中文字幕高清在线视频| 成年女人永久免费观看视频| 中文字幕高清在线视频| 亚洲狠狠婷婷综合久久图片| 91老司机精品| 欧美黑人欧美精品刺激| 中文字幕久久专区| 男女床上黄色一级片免费看| 免费观看精品视频网站| 国产高清视频在线观看网站| 免费一级毛片在线播放高清视频| 久久精品影院6| 窝窝影院91人妻| 国产一区二区三区在线臀色熟女| 亚洲av片天天在线观看| 久久香蕉国产精品| 99久久99久久久精品蜜桃| 久久久精品欧美日韩精品| 国产成人精品久久二区二区免费| 舔av片在线| 搞女人的毛片| 国产高清有码在线观看视频| 99久久精品一区二区三区| 国内毛片毛片毛片毛片毛片| 国产亚洲精品综合一区在线观看| 叶爱在线成人免费视频播放| 色老头精品视频在线观看| 免费搜索国产男女视频| 欧美不卡视频在线免费观看| 老鸭窝网址在线观看| av福利片在线观看| 性色avwww在线观看| 动漫黄色视频在线观看| 国产真实乱freesex| 午夜免费成人在线视频| 黄色片一级片一级黄色片| 日本与韩国留学比较| 久久久久国产一级毛片高清牌| 亚洲国产精品999在线| 欧洲精品卡2卡3卡4卡5卡区| 黄色女人牲交| 制服人妻中文乱码| 亚洲精品粉嫩美女一区| 人人妻人人看人人澡|