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

    北京云蒙山變質(zhì)核雜巖大水峪韌性剪切帶的應變特征及構(gòu)造意義*s

    2021-09-24 05:41:58趙騰格侯泉林石夢巖丁林
    巖石學報 2021年8期
    關鍵詞:東段大水渦度

    趙騰格 侯泉林 石夢巖 丁林

    1. 中國科學院青藏高原研究所大陸碰撞與高原隆升重點實驗室,北京 100101

    2. 中國科學院大學地球與行星科學學院;中國科學院計算地球動力學重點實驗室,北京 100049

    3. 中國科學院青藏高原地球科學卓越創(chuàng)新中心,北京 100101

    4. 河南理工大學資源環(huán)境學院,焦作 454150

    華北克拉通在中生代以來發(fā)生了強烈的破壞與巨量的減薄,并在地殼淺部發(fā)育了一系列的伸展構(gòu)造與變質(zhì)核雜巖。云蒙山變質(zhì)核雜巖位于華北克拉通北部,是我國最早發(fā)現(xiàn)的變質(zhì)核雜巖。20世紀80年代,王玉芳(1989)在云蒙山水峪地區(qū)發(fā)現(xiàn)了我國第一條低角度正斷層。隨后,Davis與鄭亞東等人的合作研究將其確立為我國第一個變質(zhì)核雜巖(Davisetal., 1996)。三十多年來,不同的學者對云蒙山變質(zhì)核雜巖進行了包括發(fā)育機制、巖石類型、地球化學以及年代學等多方面的研究工作,但仍存在一定的爭議與討論空間。

    爭議的焦點之一是對云蒙山變質(zhì)核雜巖發(fā)育機制的認識,即是伸展作用主導還是巖漿作用主導?Zhengetal.(1988)首先提出了熱隆-滑覆-推覆的構(gòu)造模式,把河防口低角度正斷層與云蒙山北部逆沖斷層的形成歸因為云蒙山核部花崗巖基的熱力隆起,之后又將其修正為先發(fā)生南北向的縮短作用形成四合堂剪切帶,再發(fā)生地殼熔融致使花崗閃長巖隆起進而側(cè)向伸展形成大水峪韌性剪切帶。Davisetal.(1996)則認為云蒙山變質(zhì)核雜巖的發(fā)育機制是類似科迪勒拉型變質(zhì)核雜巖的簡單伸展模式。此模式假設北部的推覆作用較早發(fā)生,之后才發(fā)生了云蒙山巖體的隆升以及拆離斷層的發(fā)育,變質(zhì)核雜巖整體受伸展作用所控制。Passchieretal.(2005)通過研究云蒙山花崗閃長巖體侵入于周邊變閃長巖體中同構(gòu)造花崗巖脈體的運動學與地球化學特征,推斷巖漿熱隆作用是云蒙山變質(zhì)核雜巖形成的主要原因并將四合堂剪切帶與大水峪剪切帶解釋為同期的伸展剪切帶。陳印(2014)則認為云蒙山變質(zhì)核雜巖的形成更接近滾動樞紐模式,主張變質(zhì)核雜巖形成的主導因素不是云蒙山巖體的侵位而是區(qū)域性的拉張作用。云蒙山野外實地勘探的結(jié)果表明四合堂剪切帶的上盤是向南倒轉(zhuǎn)的大型背斜,并且存在大量的“褶疊帶”與地層重復現(xiàn)象,所以部分觀點將云蒙山地區(qū)解釋為簡單剪切成因有所不妥。此外,大水峪韌性剪切帶與其上覆的河防口低角度正斷層向NE延伸,并不環(huán)繞包圍云蒙山巖基,這與科迪勒拉型變質(zhì)核雜巖所具有的“斷層包圍的內(nèi)部變質(zhì)核”的特征不符。因此云蒙山變質(zhì)核雜巖的形成機制還有待進一步探討。

    爭議的焦點之二是大水峪韌性剪切帶內(nèi)應變類型的變化及其與變質(zhì)核雜巖形成的關系。夏浩然(2011)認為伸展作用早期的韌性變形以純剪切為主,其后簡單剪切組分逐漸增加,在伸展變形后期過渡為簡單剪切。李建波(2012)認為大水峪韌性剪切帶的形成是地殼伸展與巖漿底辟聯(lián)合作用的產(chǎn)物。剪切帶在發(fā)育的早期經(jīng)歷了簡單剪切作用,晚期則逐漸過渡為以純剪切為主的變形。陳印(2014)對大水峪韌性剪切帶內(nèi)不同位置的5塊樣品進行了運動學渦度計算,推斷大水峪韌性剪切帶各處均發(fā)生以簡單剪切為主的變形。即便如此,這種數(shù)量的樣品數(shù)據(jù)似乎并不能很好地反應剪切帶內(nèi)的運動學渦度特征。

    無論是云蒙山變質(zhì)核雜巖的形成機制還是大水峪韌性剪切帶內(nèi)的應變特征變化,其核心的問題在于云蒙山地區(qū)構(gòu)造演化過程中巖漿作用與伸展作用的主次關系。巖漿作用所造成侵位隆升為變形過程提供了純剪切組分,伸展作用所造成的韌性剪切為變形過程提供了簡單剪切組分。厘清大水峪韌性剪切帶形成過程中巖漿作用與伸展作用的關系,能更好地了解云蒙山地區(qū)的構(gòu)造演化過程。

    通過對前人的研究方法與數(shù)據(jù)進行分析,覺得仍存在一些有待商榷的問題,比如大水峪韌性剪切帶內(nèi)部的運動學渦度值。渦度本身是一個瞬時概念,指變形過程中某一階段的剪切狀態(tài)。除非剪切過程是穩(wěn)定遞進變形,否則變形過程中的運動學渦度值并不恒定。在運動學渦度的研究過程中,有不同的學者提出了多個研究方法,不同的方法得到的計算結(jié)果存在些許差異。因為在實際操作中,研究者觀察到的是變形的最終產(chǎn)物,并未涉及到變形事件的全過程,所以不同方法所計算出的運動學渦度值屬于變形的不同時期。盡管如此,運動學渦度值的計算仍不失其意義,通過多種方法之間的對比我們還是能夠?qū)^(qū)域內(nèi)的變形事件有所了解。在大水峪韌性剪切帶中,夏浩然(2011)和李建波(2012)所測量的剖面與傾向以較大角度斜交,無法判斷其變化是沿剪切方向的變化還是沿走向的不均性變化。因此這樣測量出的結(jié)果難以代表大水峪韌性剪切帶內(nèi)的應變情況。陳印(2014)采用石英C軸組構(gòu)法測量得到河防口處的運動學渦度值為0.98。由最大有效力矩準則可知:最大主應力軸σ1與伸展褶劈理(C’)的夾角為55°,因此σ1與剪切帶法線的夾角ξ不會大于35°,在公式Wk=sin2ξ的限制下,運動學渦度Wk的最大值不會大于0.94。換言之,因巖石內(nèi)聚力的存在,在同一構(gòu)造應力場作用下,剪切帶必然遭受減薄,Wk<0.94。因此大水峪韌性剪切帶內(nèi)的應變特征還存在一定的討論空間。

    本文以大水峪韌性剪切帶為研究對象,在剪切帶的不同部位沿傾向方向進行樣品采集。在鏡下對采集到的定向樣品進行有限應變測量,使用石英長石溫度計與EBSD石英組構(gòu)分析來確定剪切帶的變形溫度區(qū)間,采用極摩爾圓法、伸展褶劈理法以及Rs/θ法對剪切帶的運動學渦度進行計算。以此判斷巖漿作用與伸展作用在剪切帶發(fā)育過程中的影響,為云蒙山變質(zhì)核雜巖的構(gòu)造演化提供新的證據(jù)。

    1 區(qū)域地質(zhì)背景

    云蒙山位于華北克拉通北部,燕山構(gòu)造帶西段,地處北京北部的密云水庫以西(圖1)。研究區(qū)主要出露有太古代的變質(zhì)基底、中新元古界沉積巖蓋層以及侵入其中的云蒙山花崗巖基以及石城閃長巖體(張慧等, 2018)。

    圖1 云蒙山地質(zhì)簡圖(據(jù)陳印等, 2014修改)Fig.1 Geological sketch map of Yunmengshan (modified after Chen et al., 2014)

    區(qū)域內(nèi)的構(gòu)造主體為NE-SW向的云蒙山背形。背形的北緣為四合堂韌性剪切帶,東南緣為大水峪韌性剪切帶,北西側(cè)為崎峰茶高角度正斷層,南東側(cè)為河防口低角度正斷層(圖1)。云蒙山背形、大水峪韌性剪切帶和河防口低角度正斷層共同組成了云蒙山變質(zhì)核雜巖。

    云蒙山背形長約60km,寬約15km。其主要組成部分為晚侏羅-早白堊世的云蒙山花崗巖基以及兩翼的元古界、寒武系和中侏羅世的未變形的火山碎屑巖。云蒙山巖基出露面積約300km2,以中、細粒二長花崗巖為主(姬廣義等, 2004)。巖基南側(cè)與長園巖體相接,東側(cè)為密云閃長巖和沙坨子花崗巖。云蒙山巖基在大水峪剪切帶的形成過程中發(fā)育了不同程度的韌性與脆性變形構(gòu)造,由巖基向外逐漸由片麻巖過渡到糜棱巖,再由糜棱巖過渡到斷層角礫巖和斷層泥。巖體中的脈體與面理呈不同角度相交。前人的研究資料表明,云蒙山巖體侵位的時間為146~137Ma(Davisetal., 1998; Shietal., 2009; Wangetal., 2012; 劉翠等, 2004; 陳印等, 2014)。

    大水峪韌性剪切帶位于云蒙山變質(zhì)核雜巖的東南緣。巖石組成主要為片麻狀與糜棱狀的花崗巖和閃長巖,普遍發(fā)育有旋轉(zhuǎn)碎斑、S-C組構(gòu)、伸展褶劈理、里德爾破裂等韌性與脆性變形構(gòu)造且普遍具有清晰的透入性面理和線理。這些現(xiàn)象均指示著上盤向SE的剪切作用。剪切帶的走向在梨樹溝附近由NE變?yōu)镹NE,因此本文將大水峪韌性剪切帶分為北東段和南西段兩個部分進行研究。北東段韌性剪切帶走向NNE,沿梨樹溝-石城一線展布,長約10km,寬約為1km,傾角20°~45°,上盤向SE運動。南西段韌性剪切帶走向NE,沿河防口-白道峪-梨樹溝一線展布,長約15km,寬約3km,傾角20°~50°,上盤向SE運動。在石城附近,剪切帶切穿云蒙山花崗巖基與石城閃長巖體。帶內(nèi)可見多期次的花崗巖侵入脈體,受變形作用影響與面理呈不同角度相交。年代學研究表明大水峪剪切帶的活動時間為135~114Ma(Davisetal., 1996; Shietal., 2009; Wangetal., 2012; Zhuetal., 2015)。

    河防口低角度正斷層疊加在大水峪韌性剪切帶之上,總長約60km,走向NE-NNE,傾向SE,傾角15°~45°。關于河防口低角度正斷層的成因,目前尚存在一定的爭議,陳印(2014)認為斷層和大水峪韌性剪切帶都受控于同樣的伸展機制,與云蒙山巖體侵位無關,是135~100Ma地殼隆升過程中由塑性變形轉(zhuǎn)變?yōu)榇嘈詳嗔训慕Y(jié)果。Zhuetal.(2015)認為河防口低角度正斷層是受到均衡回彈影響而發(fā)育的地殼淺層次伸展構(gòu)造,持續(xù)時間從114Ma到早白堊世末期。

    2 大水峪韌性剪切帶的變形特征

    2.1 宏觀變形特征

    大水峪韌性剪切帶的南西段切穿長園花崗巖體,北東段穿過云蒙山花崗巖基與石城閃長巖體。剪切帶的上盤為侏羅系、元古宇以及太古宇的碎裂和角礫狀巖石,下盤為角礫狀太古宇片麻巖及侵入其中的花崗質(zhì)片麻巖。

    大水峪韌性剪切帶的不同部位發(fā)育有不同變形程度的構(gòu)造巖,由巖基向斷層方向依次出露片麻巖、糜棱巖、初糜棱巖以及斷層角礫巖,剪切帶南西段的青龍峽附近可見下寒武統(tǒng)君山組豹皮灰?guī)r向南逆沖推覆在中侏羅統(tǒng)的火山碎屑巖之上,斷層面北傾,傾角約15°,上盤為一系列高角度的正斷層(圖2a)。在白道峪村村口,可見花崗質(zhì)糜棱巖的露頭,其上發(fā)育有良好的眼球狀旋轉(zhuǎn)碎斑和S-C組構(gòu),指示向SE的剪切作用(圖2b)。剪切帶北東段的云蒙峽剖面內(nèi)可見同樣指示SE方向剪切作用的伸展褶劈理(圖2c)。此外,剖面內(nèi)還發(fā)現(xiàn)了鞘褶皺的存在(圖2d),這指示著云蒙山早期南北向的構(gòu)造運動。鞘褶皺的形成有多種方式。當剪應變大于10時,直接形成沿剪切方向的拉長不對稱褶皺;當剪應變較小時,鞘褶皺的形成有兩種可能:一是先期褶皺在剪切作用下樞紐被彎曲,形成翼間角很小的刀鞘狀褶皺,這是疊加變形的結(jié)果;二是先期褶皺中存在原始層厚度不等或褶皺的原始層面與剪切方向斜交以及其他的局部不均一性,在遞進的剪切作用下發(fā)育為鞘褶皺(張逸鵬等, 2012)。有限應變測量結(jié)果顯示大水峪韌性剪切帶內(nèi)的剪應變均不大于5,屬于剪應變較小的情況。前人的研究結(jié)果表明云蒙山地區(qū)早先存在過一次由北向南的逆沖推覆,因此鞘褶皺可能是在這個階段初步發(fā)育,之后在經(jīng)歷了東西向的剪切作用后形成的。另一方面,鞘褶皺的發(fā)育說明剪切帶在形成過程中發(fā)生了韌性減薄。

    2.2 顯微構(gòu)造特征

    鏡下觀察可以發(fā)現(xiàn),大水峪韌性剪切帶中同樣保留了不同特征的韌性與脆性變形構(gòu)造。在靠近河防口斷層的部位,韌性變形的S-C組構(gòu)(圖3a)與脆性變形的多米諾構(gòu)造(圖3b)均指示向SE的剪切。此外,大水峪韌性剪切帶不同位置的石英與長石的重結(jié)晶機制存在差異。在靠近巖基的部位,通常為花崗質(zhì)片麻巖或花崗質(zhì)糜棱巖,鏡下觀察可見長石膨凸重結(jié)晶與亞顆粒旋轉(zhuǎn)重結(jié)晶現(xiàn)象共存(圖3c),斜長石的邊部發(fā)育明顯的蠕英結(jié)構(gòu)(圖3d),石英則均為顆粒邊界遷移重結(jié)晶(圖3e)。在剪切帶中間部位,北東段因為云蒙山巖基與石城閃長巖體侵入接觸,主要為花崗質(zhì)初糜棱巖與閃長質(zhì)初糜棱巖互層。長石的重結(jié)晶機制變?yōu)榕蛲怪亟Y(jié)晶(圖3f),蠕英結(jié)構(gòu)的數(shù)量減少,石英則過渡為顆粒邊界遷移重結(jié)晶與亞顆粒旋轉(zhuǎn)重結(jié)晶共存(圖3g)。南西段主要為花崗質(zhì)與花崗閃長質(zhì)的初糜棱巖,石英與長石重結(jié)晶機制基本保持不變。在靠近河防口斷層的部位,北東段位于石城閃長巖體內(nèi),長石顆粒均為膨凸重結(jié)晶,石英多為亞顆粒旋轉(zhuǎn)重結(jié)晶(圖3h),但個別樣品中的石英仍可觀察到顆粒邊界遷移重結(jié)晶現(xiàn)象。南西段的長石顆粒同樣為膨凸重結(jié)晶,石英則過渡為亞顆粒旋轉(zhuǎn)重結(jié)晶。

    圖3 大水峪韌性剪切帶鏡下顯微構(gòu)造(a)長石的S-C組構(gòu);(b)斜長石的多米諾構(gòu)造;(c)鉀長石的亞顆粒旋轉(zhuǎn)重結(jié)晶現(xiàn)象;(d)斜長石內(nèi)部的蠕英結(jié)構(gòu);(e)石英的顆粒邊界遷移重結(jié)晶;(f)長石的膨凸重結(jié)晶;(g)石英的亞顆粒旋轉(zhuǎn)重結(jié)晶(SGR)與顆粒邊界遷移重結(jié)晶(GBM)共存;(h)石英的亞顆粒旋轉(zhuǎn)重結(jié)晶.Pl-斜長石;Kfs-鉀長石;Qtz-石英;Hbl-角閃石;Mus-白云母;Bt-黑云母;SGR-亞顆粒旋轉(zhuǎn)重結(jié)晶;GBM-顆粒邊界遷移重結(jié)晶Fig.3 Microstructures of rocks of the Dashuiyu ductile shear zone(a) S-C fabric of feldspar; (b) domino structure of plagioclase; (c) subgrain rotational recrystallization of K-feldspar; (d) myrmekitic texture within plagioclase; (e) grain boundary migration recrystallization of quartz; (f) bulging recrystallization of feldspar; (g) coexistence of subgrain rotation recrystallization (SGR) and grain boundary migration recrystallization (GBM) in quartz; (h) grain boundary migration recrystallization of quartz. Pl-plagioclase; Kfs-potash feldspar; Qtz-quartz; Hbl-amphibole; Mus-muscovite; Bt-biotite; SGR-subgrain rotational recrystallization; GBM-grain boundary migration recrystallization

    前人的研究結(jié)果表明,石英與長石在不同的溫度區(qū)間內(nèi)存在著不同的重結(jié)晶樣式(向必偉等, 2007)。通過觀察長石與石英在顯微鏡下的動態(tài)重結(jié)晶形式,可以限定大水峪韌性剪切帶的變形溫度區(qū)間。因此,初步推斷大水峪韌性剪切帶的構(gòu)造巖相在高綠片巖相與低角閃巖相之間變化。變形溫度區(qū)間為480~650℃(大部分樣品的變形溫度低于600℃)。這與張慧等(2018)對云蒙峽剖面使用二長石溫度計估算的結(jié)果(450~630℃)基本吻合。

    3 有限應變測量與運動學渦度估算

    3.1 有限應變測量

    發(fā)生過構(gòu)造運動的區(qū)域,其間的巖石會在一定程度上發(fā)生變形,產(chǎn)生應變。通過測量巖石中某些標志物體(礫石、構(gòu)造角礫、化石、鮞粒、氣孔以及枕狀熔巖等)的應變狀態(tài),可以了解巖石的變形現(xiàn)象,從而進一步了解研究區(qū)的構(gòu)造屬性與經(jīng)歷的構(gòu)造事件。這種定量分析研究區(qū)域應變狀況的手段稱為“有限應變測量”。

    有限應變測量首先要確定所測量的應變標志體。對所采集的樣品進行了鏡下觀察研究發(fā)現(xiàn),剪切帶內(nèi)的石英多數(shù)發(fā)生不同程度的重結(jié)晶作用,如亞顆粒旋轉(zhuǎn)重結(jié)晶以及顆粒邊界遷移重結(jié)晶,不適合作為有限應變測量的標志體。剪切帶內(nèi)的長石多發(fā)生脆性破裂和低溫膨凸重結(jié)晶。即使在出現(xiàn)長石亞顆粒旋轉(zhuǎn)重結(jié)晶的樣品中,也可觀察到部分膨凸重結(jié)晶的長石顆粒。因此可將長石在變形過程中的變化視為剛體旋轉(zhuǎn)。長石在不同的變質(zhì)條件下具有不同的變形現(xiàn)象:在低級變質(zhì)條件(<400℃)下主要為脆性破裂,在中低級變質(zhì)條件(400~500℃)和中級變質(zhì)條件下(450~600℃)出現(xiàn)位錯攀移與亞顆粒旋轉(zhuǎn)動態(tài)重結(jié)晶作用,在高級變質(zhì)條件(>600℃)下則以亞顆粒旋轉(zhuǎn)重結(jié)晶為主(胡玲等, 2009)。由前文可知,大水峪韌性剪切帶的變形溫度區(qū)間為480~650℃,且大多數(shù)樣品位于中低級和中級變質(zhì)條件下。Bailey and Eyster (2003)的研究結(jié)果表明:在中低級韌性變形過程中,長石顆粒間的滑動較少,近于剛體旋轉(zhuǎn)。此外,大水峪韌性剪切帶內(nèi)主要發(fā)育糜棱巖和初糜棱巖,基質(zhì)含量約在50%左右。鏡下觀察也可發(fā)現(xiàn)長石顆粒較多而石英多為裂隙充填。綜合以上原因,本研究選取鏡下較為完整的長石顆粒作為應變測量標志體。

    有限應變測量的基礎是要將變形巖石中標志體的變形狀態(tài)與變形前進行比較。通過測量應變標志體的主應變軸的方位和大小等變形數(shù)據(jù),算出研究區(qū)域內(nèi)各樣點的應變狀態(tài),求取各點應變橢球體主軸的方位、長短軸比以及Flinn指數(shù)(k)等參數(shù),進而推算研究區(qū)域的減薄量和韌性位移量等數(shù)據(jù)。

    因此,應變測量的首要任務是要確定主應變軸的方位及其大小。一般情況下,主應變軸的方位由糜棱巖中的面理和拉伸線理的方位來確定。規(guī)定拉伸線理方向為X軸,平行面理垂直線理方向為Y軸,垂直面理方向為Z軸。XYZ三個軸共同組成了應變分析所需要的XY面(最大壓扁面)、XZ面和YZ面。通過測量三個面中任意兩個面上的長短軸比,將其標準化到體積為4π/3的橢圓上,即可得到XYZ三個軸上的應變參數(shù)。

    本次研究在大水峪韌性剪切帶的5條剖面內(nèi)共采集41塊定向樣品進行研究。平行于XZ和YZ面切制定向薄片。使用Straindesk軟件對鏡下照片中的長石顆粒輪廓進行勾勒,即可得到長石顆粒的長短軸比以及長軸與剪切帶邊界夾角等應變測量數(shù)據(jù)(李志勇和曾佐勛, 2006)。每個樣品采集30個以上的長石顆粒數(shù)據(jù),并對測量得到的應變數(shù)據(jù)計算調(diào)和平均值來減少誤差。詳細測量結(jié)果見表1。

    表1 大水峪韌性剪切帶有限應變測量、運動學渦度分析數(shù)據(jù)表Table 1 Data table for finite strain measurement and kinematic vorticity analysis of the Dashuiyu ductile shear zone

    3.2 應變特征分析

    許氏圖解和Flinn圖解是分析地質(zhì)體應變類型的常用方法。通過將應變測量得到的數(shù)據(jù)投影到圖上可以直觀地表達應變橢球的形態(tài)以及其所代表的應變類型。

    剪切帶北東段的Flinn圖解投圖結(jié)果顯示云蒙峽剖面和南石城剖面的點位較為分散,北石城剖面內(nèi)的點則較為集中,但Flinn指數(shù)k的平均值相差不大,分別為0.23、0.28和0.24,均指示壓扁應變(圖4a)。許氏圖解中三條剖面的Lode參數(shù)ν均位于0到1之間,結(jié)果與Flinn圖解的結(jié)果相同,仍指示壓扁變形的橢球(圖4b)。剪切帶南西段的Flinn圖解中白道峪剖面與小水峪剖面的投點位置更靠近k=0.1的軸,k的平均值分別為0.13和0.15,這指示著南西段的應變狀態(tài)接近軸對稱壓縮(圖4c)。許氏圖解的Lode參數(shù)ν同樣位于0到1之間,指示壓扁應變(圖4d)。

    圖4 大水峪韌性剪切帶的Flinn圖解與許氏圖解Fig.4 Flinn and Hsǖ diagram of the Dashuiyu ductile shear zone

    綜合分析表明大水峪韌性剪切帶在形成過程中的主要應變類型為壓扁應變,其由巖漿作用所主導的巖基侵位隆升所提供,且南西段的壓扁作用要大于北東段。此外,應變測量結(jié)果顯示南西段與北東段的X軸與Y軸均發(fā)生了拉長。剪切帶南西段的寬度(2~3km)大于北東段(1km)。前人研究表明,云蒙山東北部巖體(馮家峪巖體、西白蓮峪巖體和大東溝巖體)中未觀察到明顯的構(gòu)造變形現(xiàn)象(陳印, 2014)。因此,推測大水峪韌性剪切帶的北東段的伸展拆離過程受到云蒙山東北部巖體(西白蓮峪巖體、馮家峪巖體和大東溝巖體)的阻擋,進而沿剪切帶走向方向發(fā)生了一定程度上的擠壓作用。綜合上述討論,可知剪切帶在發(fā)育過程受巖漿作用所主導的純剪切的影響較大,伸展拆離作用所提供的簡單剪切作用于剪切帶形成晚期。

    3.3 韌性位移量與減薄率計算

    表2 大水峪韌性剪切帶中各剖面的韌性位移量Table 2 Ductile displacement of each section in the Dashuiyu ductile shear zone

    剪切帶在形成過程中受到垂直和水平方向的力,必然會出現(xiàn)一定程度上的體積變化,韌性的減薄與增厚不可避免。Simpson and De Paor (1993)提出韌性剪切帶在形成過程中受到次簡單剪切作用影響而發(fā)生厚度變化,剪切帶在其長度增加、體積變小和第三位方向伸展的情況下可發(fā)生韌性減薄。Wallis (1995)認為垂直剪切帶邊界方向的減薄率(減薄量與原厚度的比值)與運動學渦度Wk以及XZ面的軸比Rs

    分別采用公式法和極摩爾圓圖解法對云蒙峽剖面和小水峪剖面的XZ面和YZ面上的減薄率進行計算(表3),發(fā)現(xiàn)XZ面的減薄率高達60%~90%,YZ面的減薄率也在60%以上。這同樣指示著大水峪韌性剪切帶的形成過程中受壓扁作用的影響較大,與Flinn圖解與許氏圖解的投圖結(jié)果一致。

    表3 大水峪韌性剪切帶的韌性減薄率Table 3 Ductile thinning rate of the Dashuiyu ductile shear zone

    4 運動學渦度

    渦度(vorticity)原本是流體力學中描述流體旋轉(zhuǎn)程度的概念,地質(zhì)學上把它作為描述巖石變形過程中內(nèi)部旋轉(zhuǎn)程度的向量。為了度量簡單剪切與純剪切各組分在變形過程中所占的比例,地質(zhì)學家又引入了運動學渦度(Wk)的概念。運動學渦度可以更加清晰地度量剪切帶中簡單剪切與純剪切各組分在變形過程中所占的比例,進而更好地解釋剪切帶的運動學性質(zhì)。當0

    Bobyarchick(1986)把運動學渦度(Wk)的表達式簡化為Wk=cosα。式中,α為兩特征方向間的夾角。純剪切的運動學渦度值為0,簡單剪切的運動學渦度值為1,一般剪切的運動學渦度值介于0和1之間(鄭亞東等, 2008)。

    運動學渦度的研究與使用迄今已有40年,計算運動學渦度的方法也在不斷完善和進步。目前,估算運動學渦度的主流方法有極摩爾圓法(Zhang and Zheng, 1997)、伸展褶劈理法(鄭亞東等, 2007)以及Rs/θ圖解法(Tikoff and Fossen, 1993)。本文結(jié)合上述三種方法對大水峪韌性剪切帶的運動學渦度值(Wk)進行估算。

    4.1 極摩爾圓法

    Simpson and De Paor (1993)最先提出了極摩爾圓的構(gòu)建方法。Zhang and Zheng (1997)對其進行了改進,提出了更為實用的利用極摩爾圓計算運動學渦度的方法。本文采用這個方法,對大水峪韌性剪切帶內(nèi)5條剖面上的定向樣品進行運動學渦度的計算。利用有限應變測量得到的XZ面長短軸比RS-XZ以及其長軸與剪切帶邊界的夾角θ來構(gòu)建極摩爾圓,從而得到兩特征向量的夾角α,再利用公式Wk=cosα計算運動學渦度(表1)。將計算得到的運動學渦度值結(jié)合采樣位置投影在剖面圖上,可以得到大水峪韌性剪切帶走向方向的渦度變化趨勢。

    剪切帶北東段的運動學渦度值變化幅度較大。云蒙峽剖面的運動學渦度為0.08~0.84,平均值為0.42;南石城剖面的運動學渦度為0.06~0.5,平均值為0.25;北石城剖面的運動學渦度為0.11~0.65,平均值為0.44(表1)。結(jié)合剖面在地質(zhì)圖上的位置(圖1)分析得到大水峪剪切帶北東段的運動學渦度變化趨勢由巖基向斷層方向先減小后增大,剪切類型總體仍為純剪切(圖5)。

    圖5 大水峪韌性剪切帶北東段采樣剖面與運動學渦度值變化Fig.5 Sampling locations and kinematic vorticity variations in the northeast section of the Dashuiyu ductile shear zone

    剪切帶南西段的渦度變化范圍較小。白道峪剖面的渦度為0.44~0.66,平均值為0.54;小水峪剖面的渦度為0.60~0.85,平均值為0.77。結(jié)合剖面地理位置分析得到大水峪剪切帶南西段由北向南運動學渦度值逐漸增大,純剪切組分減少,剪切類型逐漸由純剪切過渡到簡單剪切(圖6)。

    圖6 大水峪韌性剪切帶南西段采樣剖面及運動學渦度值變化Fig.6 Sampling locations and kinematic vorticity variations in the southwest section of the Dashuiyu ductile shear zone

    4.2 伸展褶劈理(C’)法

    伸展褶劈理(C’)形成于剪切帶晚期的低溫階段,代表的是剪切帶形成晚期的運動學渦度(鄭亞東等, 2008)。鄭亞東(1999)提出的“最大有效力矩準則”可用于伸展褶劈理發(fā)育時運動學渦度的計算。C′面理與C面理間的夾角為0°~35°,C′與最大主應力(σ1)間的夾角為55°。測量糜棱巖中形成的伸展褶劈理(C’)與剪切帶面理(C)的夾角,可得最大主應力軸(σ1)與剪切帶法線的夾角(ξ),通過公式 可以計算得到運動學渦度Wk(鄭亞東等, 2008)。

    本文針對白道峪、云蒙峽、南石城和北石城4條剖面中選取能觀察到伸展褶劈理發(fā)育的XZ面薄片進行測量計算。為了減小誤差,每張薄片測量30組數(shù)據(jù)以上并計算調(diào)和平均值(表1)。白道峪剖面的運動學渦度值為0.27~0.47,平均值為0.39。云蒙峽剖面的運動學渦度值為0.10~0.59,平均值為0.40;南石城剖面的運動學渦度值為0.15~0.42,平均值為0.30;北石城剖面的運動學渦度值為0.22~0.62,平均值為0.44。四條剖面的伸展褶劈理法計算結(jié)果均指示純剪切為主的變形作用。結(jié)合極摩爾圓法計算得到的運動學渦度數(shù)據(jù)與剖面采樣位置(圖5、圖6)分析,說明剪切帶形成的后期純剪切組分占比較高且北東段的運動學渦度存在由巖基向外先減小后增大的趨勢。

    4.3 RS/θ曲線法

    在某些剪切帶中,有限應變橢球體可以由垂直于剪切帶和平行于剪切方向的截面中的變形標志體計算。Fossen and Tikoff (1993)的研究表明如果應變橢圓的軸率RS和應變橢圓長軸與高應變帶邊界夾角θ已知,則可通過關系圖解估算運動學渦度值的范圍與地質(zhì)體所經(jīng)歷的擠壓或拉張作用。劉江等(2012)基于前人的研究,繪制了平均運動學渦度Wm關于RS與θ的等值線圖。本文主要通過測量S面理與C面理間的夾角來獲得θ的值,結(jié)合特征向量間的夾角α以及長短軸比RS進行投圖(圖7)。結(jié)果顯示,剪切帶北東段與南西段均發(fā)生減薄的一般剪切作用,且北東段內(nèi)的渦度變化范圍要大于南西段。

    圖7 大水峪韌性剪切帶RS/θ法圖解Fig.7 RS/θ diagrams of the Dashuiyu ductile shear zone

    此外,Tikoff and Fossen (1993)還提出了轉(zhuǎn)換擠壓和轉(zhuǎn)換拉張下運動學渦度(Wk)分別與θ和α的關系圖。將θ和α投影在圖上,發(fā)現(xiàn)大水峪韌性剪切帶的南西段和北東段均為以純剪切為主的變形(圖8a, b),且均指示以純剪切為主的轉(zhuǎn)換擠壓作用(圖8c, d)。白道峪剖面的純剪切組分多于小水峪剖面,這與極摩爾圓法測得的結(jié)果一致。盡管剪切帶北東段剖面的長度要小于剪切帶南西段剖面的長度,但北東段剖面內(nèi)的渦度變化范圍要大于南西段剖面的渦度變化范圍。結(jié)合云蒙山地區(qū)的地質(zhì)背景,這說明大水峪韌性剪切帶的北東段的相較南西段發(fā)生了更為復雜的構(gòu)造演化事件。

    圖8 轉(zhuǎn)換擠壓和轉(zhuǎn)換拉張下運動學渦度(Wk)分別與θ和α的關系Fig.8 The relationship between kinematic vorticity (Wk) and θ and α, under transpression and transtension

    5 EBSD石英組構(gòu)分析

    在構(gòu)造地質(zhì)學研究中,巖石組構(gòu)分析是研究巖石物理性質(zhì)、變形機制以及運動學規(guī)律的重要方向和手段。僅在露頭尺度上對巖石的變形現(xiàn)象進行研究不僅不精確而且難以定量分析。因此,在巖石組構(gòu)的研究中常常需要借助一些技術手段。

    電子背散射衍射(Electron backscatter diffraction,簡稱EBSD)技術是一種近年來在地質(zhì)領域很流行的技術。其原理是結(jié)合波譜和能譜成分分析等設備使用高分辨率的掃描電子顯微鏡來確定晶體的粒度、邊界、取向差、物相鑒定與應變等信息(黃學猛等, 2016)。通過EBSD技術分析能夠把掃描得到的晶體參數(shù)以不同的圖件表示,其中就包括礦物的CPO樣式。

    樣片表面晶體的優(yōu)選方向稱為CPO(crystal preferred orientation)。不同的CPO樣式可以指示不同的運動方向、確定變形溫度區(qū)間和度量應變大小。CPO樣式在共軸變形與非共軸變形中呈現(xiàn)不同的對稱關系,前者為對稱關系,后者為單斜對稱關系(Langilleetal., 2010)。此外,CPO樣式也能指示不同的變形溫度區(qū)間。以石英為例,變形溫度由低到高依次為底面滑移、菱面滑移、柱面滑移和柱面滑移。

    本研究采用石英作為研究載體,實驗在中國地質(zhì)科學院地質(zhì)研究所大陸動力學實驗室完成,測試方法及流程見文獻(劉俊來等, 2008)。根據(jù)所測試的5個樣品呈現(xiàn)出不同的CPO樣式,并將其結(jié)合采樣位置繪制在地質(zhì)圖上(圖9)。由圖9可知,5塊樣品的石英C軸組構(gòu)圖均指示向SE的剪切作用,但環(huán)帶類型和極密位置存在差異。剪切帶兩端靠近巖基的樣品(YMX19554)和靠近斷層的樣品(YMX121401和SC694)呈現(xiàn)為靠近Y軸極密的柱面滑移,指示高溫變形。剪切帶中間部位的樣品(DX682和SC692)則為Z軸與YZ軸之間的極密,為指示中低溫變形的菱面滑移和底面滑移。這說明剪切帶中間部位存在不同溫度巖片疊加的現(xiàn)象并且變形溫度由巖基向外呈現(xiàn)先降低后升高的變化趨勢。

    圖9 大水峪韌性剪切帶北東段石英EBSD組構(gòu)及采樣位置Fig.9 EBSD fabric and sampling location of quartz in the northeast part of the Dashuiyu ductile shear zone

    6 討論

    6.1 變形條件

    糜棱巖通常形成于動力變質(zhì)條件下。Tullisetal.(2000)指出:韌性剪切帶不僅是巖石的變形帶,同時也是變質(zhì)帶。在不同的變質(zhì)環(huán)境下,糜棱巖呈現(xiàn)出不同的礦物組合和變形特征。這些現(xiàn)象同時也可以反映出巖石形成時的物理環(huán)境。大水峪韌性剪切帶北東段和南西段的礦物組成略有不同。北東段的三條剖面中可以觀察到綠泥石、綠簾石和白云母等綠片巖相環(huán)境下的典型礦物。此外還發(fā)現(xiàn)了蠕英結(jié)構(gòu)這種中高溫變質(zhì)條件下的產(chǎn)物。蠕英結(jié)構(gòu)的成分由斜長石和石英組成。董彥龍等(2019)的研究表明,變形花崗巖中蠕英結(jié)構(gòu)的形成在一定程度上會受到剪切帶發(fā)育過程中力的作用。因此,蠕英結(jié)構(gòu)很可能形成于大水峪韌性剪切帶的發(fā)育過程之中。綜合以上特征,結(jié)合長石的不同重結(jié)晶現(xiàn)象,可以推斷大水峪韌性剪切帶的北東段經(jīng)歷了從綠片巖相到角閃巖相的變化。與北東段不同的是,南西段的兩條剖面中均普遍發(fā)育蠕英結(jié)構(gòu)以及石英顆粒邊界遷移重結(jié)晶等中高溫條件下的變形現(xiàn)象。因此推測大水峪剪切帶南西段的巖石大部分形成于角閃巖相條件下。

    由石英長石礦物溫度計以及EBSD石英組構(gòu)分析可知:大水峪韌性剪切帶北東段的變形溫度區(qū)間變化較大,變形主要發(fā)生在中低溫的綠片巖相條件下,少數(shù)發(fā)生在中高溫的角閃巖相中。剪切帶內(nèi)的溫度由巖基向斷層方向先降低后升高。前人的研究表明:云蒙山東側(cè)的冷卻速率為28.6℃/Myr(500~300℃)和50.0℃/Myr(300~200℃)(陳印等, 2014)。剪切帶整體的變形溫度在480~650℃之間,幾百米之內(nèi)出現(xiàn)如此大的溫度變化可能是因為伸展作用與走滑擠壓作用將不同變形溫度的巖片堆疊到了一起。

    6.2 變形特征分析

    Flinn圖解和許氏圖解均指示剪切帶內(nèi)發(fā)生以壓扁作用為主的應變,但北東段的壓扁程度要小于南西段。Lode參數(shù)ν和Flinn指數(shù)k的值指示北東段的應變類型為一般壓縮,南西段的應變類型接近軸對稱壓縮。減薄率的計算結(jié)果表明剪切帶北東段和南西段均存在60%以上的減薄,因此大水峪剪切帶在發(fā)生變形時不僅僅是X軸發(fā)生了伸長,同時Y軸也有所伸長。由應變分析的結(jié)果可知,北東段和南西段的Y軸伸長原因有所不同:剪切帶南西段是因為壓扁作用較大導致了Z軸縮短,與此同時X軸與Y軸方向同時發(fā)生了加長;剪切帶北東段則是因為剪切帶在伸展拆離過程中受到了東側(cè)巖體的阻擋,進而沿剪切帶走向方向發(fā)生了擠壓造成的Y軸拉長。

    6.3 運動學渦度變化的構(gòu)造意義

    極摩爾圓法測量運動學渦度需要使用剪切方向與拉伸方向之間的夾角以及拉伸方向的長短軸比。因此極摩爾圓法獲得的是拉伸線理和糜棱面理形成時,也就是變形早期階段的運動學渦度值 (倪金龍等, 2019)。由圖5和圖6可知,極摩爾圓法計算出的Wk值在南西段自北向南逐漸增大,由純剪切過渡到簡單剪切;在北東段由巖基向斷層方向先減小后增大,除靠近巖基和斷層方向的個別位置為簡單剪切外,剖面內(nèi)整體仍以純剪切為主。

    北東段與南西段運動學渦度值變化的差異與變形過程的不同有關。剪切帶的南西段早期受到自北向南的逆沖推覆作用影響,積累了一定程度的簡單剪切組分。剪切帶活動階段又發(fā)生了向SE的拆離,兩次向南的構(gòu)造事件疊加造成了剪切帶南西段的運動學渦度值自北向南逐漸增大。北東段剖面的走向基本與逆沖推覆作用的方向垂直,所以北東段剖面的運動學渦度值不受南北向逆沖推覆作用的影響,剖面內(nèi)的簡單剪切組分均由向SE的伸展拆離作用產(chǎn)生。此外,云蒙山背斜的樞紐為NE-SW方向(Davisetal., 1996),這導致了北東段的伸展拆離過程要比南西段更為順利,其在拆離階段所積累的簡單剪切組分也要多于南西段。應變特征分析的結(jié)果暗示了北東段在拆離過程中受東側(cè)巖體的阻擋,進而沿剪切帶走向方向發(fā)生了Y軸方向的伸長。因此北東段的運動學渦度存在由巖基向斷層方向先減小后增大的變化趨勢。夏浩然(2011)在云蒙峽內(nèi)做的EBSD石英組構(gòu)分析結(jié)果顯示云蒙峽剖面內(nèi)1200m的范圍內(nèi)出現(xiàn)了由單環(huán)帶過渡到共軛的雙環(huán)帶最后再變?yōu)閱苇h(huán)帶的變化過程。Leeetal.(1987)以及Malavieille (1993)的研究表明:在拆離斷層中,沿上盤運動方向,應變由小變大,石英的C軸組構(gòu)圖像由雙環(huán)帶過渡為單環(huán)帶(張進江等, 1998)。這指示著云蒙峽剖面內(nèi)存在著“簡單剪切-純剪切-簡單剪切”的變化趨勢,與極摩爾圓法得到的結(jié)論相符。綜上,剪切帶北東段在發(fā)育過程中經(jīng)歷了巖漿作用與伸展作用,受巖漿作用控制;南西段除了受到巖漿作用以及伸展作用的影響外,其運動學渦度還被云蒙山早期的逆沖推覆作用所影響,最終自北向南由純剪切過渡到簡單剪切。因此,大水峪韌性剪切帶在形成的早期受巖漿作用所主導的純剪切所控制。

    伸展褶劈理(C’)法是鄭亞東教授在最大有效力矩準則的基礎上提出的計算運動學渦度的方法。C’面理是先存面理受到剪應力作用發(fā)生旋轉(zhuǎn)而形成的面理,因此伸展褶劈理(C’)法所計算出的運動學渦度值反映的是較為晚期的剪切作用(王新社等, 2002)。伸展褶劈理只能在剪切帶中靠近斷層的部位觀察到,白道峪、云蒙峽、南石城和北石城剖面中的運動學渦度均指示純剪切組分占主導地位(0.1

    RS/θ法使用的θ的值為S面理與C面理的夾角,代表了晚期的運動學渦度。北東段與南西段均指示純剪切為主的轉(zhuǎn)換擠壓作用,同樣說明巖漿作用在剪切帶發(fā)育的后期占據(jù)了主導地位。

    綜上所述,極摩爾圓法計算的是早期的運動學渦度,其計算出的值相較其他兩種方法更大,說明伸展作用主要作用于剪切帶形成的早期。即使這樣,極摩爾圓法所計算出的渦度值也多小于0.75,指示剪切帶發(fā)育早期仍受純剪切控制。伸展褶劈理(C’)法與RS/θ法所指示的是晚期的運動學渦度,這同樣說明巖漿作用在剪切帶形成過程中所占的組分更大。

    6.4 變形時間與動力學機制

    大水峪韌性剪切帶的變形主要發(fā)育于云蒙山巖基和石城閃長巖體內(nèi),對它們內(nèi)部的花崗巖、閃長巖以及同構(gòu)造脈體進行實驗定年,能夠有效地限定大水峪韌性剪切帶發(fā)育不同階段的時間節(jié)點。

    云蒙山巖基的鋯石U-Pb定年結(jié)果為143~141Ma(Davisetal., 1996)、144±4Ma(劉翠等, 2004)、145~141Ma(Shietal., 2009)、137±3Ma、141±2Ma(Wangetal., 2012)以及146~141Ma(陳印, 2014)。綜合前人的數(shù)據(jù),推測云蒙山花崗巖基的侵位時間在146~137Ma之間。

    石城巖體的鋯石U-Pb定年結(jié)果為159±1Ma(Davisetal., 1996)、162~151Ma(Shietal., 2009)和158Ma(Wangetal., 2012)。侵入到石城閃長巖體中的花崗巖脈的捕獲鋯石的U-Pb年齡為158~153Ma(Wangetal., 2012; 陳印, 2014; 陳印等, 2014)。因此,石城巖體的鋯石U-Pb年齡在162~151Ma之間,屬晚侏羅世巖體。

    大水峪韌性剪切帶切穿云蒙山巖基與石城巖體。剪切帶的黑云母40Ar-39Ar測年結(jié)果為113~106Ma(王玉芳等, 1989);鉀長石40Ar-39Ar的測年結(jié)果為118Ma(Davisetal., 2002);剪切帶內(nèi)變形巖體的鉀長石、黑云母和角閃石40Ar-39Ar測年結(jié)果分別為114Ma、120Ma和125Ma(Wangetal., 2012);剪切帶中變形的云蒙山花崗巖和侵入到剪切帶中的同構(gòu)造花崗巖脈的鋯石U-Pb同位素測年結(jié)果為135~126Ma,剪切帶中的糜棱巖黑云母40Ar-39Ar測年結(jié)果為120~116Ma(陳印, 2014)。綜上,推測大水峪韌性剪切帶的活動時間為135~126Ma。前文中提到云蒙山巖基的侵位年齡為146~137Ma,這說明剪切帶并不是在巖基侵位剛發(fā)生時就開始發(fā)育的。此外,大水峪韌性剪切帶的變形溫度為480~650℃,要大于黑云母、角閃石、鉀長石K-Ar同位素封閉體系的溫度(300±50℃),因此125~114Ma為大水峪剪切帶的抬升冷卻時間。

    在云蒙山東北部的3個巖體中均未觀察到明顯的構(gòu)造變形現(xiàn)象。前人同樣對這3個巖體進行了鋯石U-Pb測年:大東溝巖體的年齡為127~128Ma(Davisetal., 1996)和129±2Ma(孫會一等, 2016);西白蓮峪巖體的年齡為124±5Ma(Wangetal., 2012)、130.4±1.3Ma(陳印, 2014)和130±1Ma(孫會一等, 2016);馮家峪巖體的年齡為128±5Ma (Wangetal., 2012)、127±1Ma、130.7±0.72Ma(陳印, 2014)以及127±1Ma(孫會一等, 2016)。綜合分析前人數(shù)據(jù),云蒙山北部3個小巖體的侵位時間應該在130~127Ma。這個時間早于剪切帶的快速抬升時間,說明東側(cè)巖體的侵位隆升確實可能對剪切帶北東段的伸展拆離過程構(gòu)成影響。

    6.5 大水峪韌性剪切帶的構(gòu)造演化

    由前文可知,大水峪韌性剪切帶的南西段與北東段經(jīng)歷了不同的構(gòu)造演化事件。剪切帶南西段早先受由北向南逆沖推覆作用的影響,積累了一定的簡單剪切組分。在135~126Ma,大水峪韌性剪切帶開始活動,簡單剪切組分進一步增多,與逆沖推覆作用的簡單剪切組分疊加,導致了運動學渦度值自北向南逐漸增大。125~114Ma,大水峪韌性剪切帶發(fā)生了快速隆升,純剪切組分迅速增加。最終,大水峪韌性剪切帶的南西段整體受一般剪切控制。

    剪切帶北東段的構(gòu)造演化過程較為復雜。162~151Ma,石城閃長巖體首先發(fā)生侵位(圖10a)。146~137Ma,云蒙山巖基發(fā)生侵位,與石城巖體發(fā)生侵入接觸(圖10b)。根據(jù)前文的研究,大水峪韌性剪切帶的變形溫度在480~650℃之間,這說明變質(zhì)核雜巖的顯著隆升發(fā)生在云蒙山巖基侵位之后,侵位深度約7km(許繼峰等, 1994)。135~126Ma,大水峪韌性剪切帶開始活動,對云蒙山巖基和石城閃長巖體產(chǎn)生剪切作用,為變形的過程提供簡單剪切組分。在此期間,石城巖體中更靠近巖基方位的高溫變形部分被向SE的剪切作用搬運到更遠的東部(圖10c)。韌性位移的距離約為1500m。130~127Ma,在大水峪韌性剪切帶活動的期間,云蒙山東部的馮家峪、西白蓮峪和大東溝巖體發(fā)生侵位。巖體內(nèi)未觀察到明顯的構(gòu)造變形,說明其未受大水峪韌性剪切帶的影響。結(jié)合韌性位移量的計算結(jié)果與東部巖體(馮家峪巖體、西白蓮峪巖體和大東溝巖體)的地理位置,推斷其阻礙了大水峪韌性剪切帶北東段的伸展拆離過程。因為剪切帶的主要形成動力是巖體侵位所引發(fā)的垂向應力,不是區(qū)域應力。在區(qū)域應力導致拆離作用發(fā)生時,東部巖體已經(jīng)成巖,因此會對向東的拆離作用產(chǎn)生阻擋以削弱剪切帶的伸展作用。因為成巖作用已經(jīng)完成,所以東部巖體內(nèi)并未發(fā)生明顯的變形現(xiàn)象。剪切帶的伸展受阻擋作用的影響,從而發(fā)生Y軸方向上的拉長。125~114Ma,大水峪韌性剪切帶發(fā)生了快速隆升,在抬升的過程中大水峪韌性剪切帶主要受到純剪切組分的影響,但其過程并未伴隨有巖漿作用。在隆升的過程中,剪切帶由韌性域過渡到脆性域,野外觀察也發(fā)現(xiàn)大水峪韌性剪切帶中云蒙山巖基與石城巖體的交界處可觀察到花崗質(zhì)初糜棱巖與閃長質(zhì)初糜棱巖互層,這可能指示著不同變形溫度的構(gòu)造巖片在拆離過程中發(fā)生了堆疊,高溫的巖片被插入到了相對低溫的巖片之中,二者互相影響造成了石英CPO的樣式中高溫變形與低溫變形疊加的現(xiàn)象(圖10d)。由于云蒙山東部西白蓮峪等巖體的阻擋,使大水峪韌性剪切帶的北東段無法向東繼續(xù)伸展拆離,進而沿剪切帶走向方向發(fā)生了走滑作用。這也導致了剪切帶中心部位的純剪切組分增大,所以大水峪剪切帶內(nèi)的運動學渦度存在著“中間低兩端高”的現(xiàn)象。隨后的113~100Ma,大水峪韌性剪切帶緩慢隆升至地表(圖10e)。

    圖10 大水峪韌性剪切帶北東段演化模式(a) 162~151Ma(石城巖體侵位);(b) 146~137Ma(云蒙山巖基侵位);(c) 135~126Ma(大水峪韌性剪切帶活動);(d) 125~114Ma(快速抬升);(e) 113~100Ma(緩慢隆升)Fig.10 Evolution model of northeast segment of the Dashuiyu ductile shear zone(a) 162~151Ma (the Shicheng rock mass intrusion); (b) 146~137Ma (the Yunmengshan batholith intrusion); (c) 135~126Ma (the Dashuiyu ductile shear zone activity); (d) 125~114Ma (rapid uplift); (e) 113~100Ma (slow uplift)

    7 結(jié)論

    (1)北京云蒙山變質(zhì)核雜巖東南緣的大水峪韌性剪切帶的變形溫度區(qū)間為480~650℃(大多數(shù)樣品的變形溫度區(qū)間位于600℃以下)。南西段的變形溫度自北向南逐漸降低,北東段的變形溫度由巖基向河防口斷層方向先降低后升高。

    (2)采用長石作為應變測量標志體,對大水峪韌性剪切帶進行應變測量與運動學渦度的估算。應變測量的結(jié)果指示剪切帶在形成過程中主要受壓扁作用的影響。南西段的運動學渦度自北向南逐漸增大,北東段的運動學渦度由巖基向外先減小后增大,二者均以純剪切組分居多且南西段的壓扁程度與韌性位移量均大于北東段。巖漿作用導致的巖基侵位隆升主要提供純剪切組分,伸展拆離作用主要提供簡單剪切組分。因此推斷云蒙山變質(zhì)核雜巖發(fā)育過程中巖漿作用居于主導地位,并且剪切帶北東段向SE方向的拆離過程受到了云蒙山東部巖體(馮家峪巖體、西白蓮峪巖體以及大東溝巖體)的阻擋。

    致謝中國地質(zhì)科學院地質(zhì)研究所梁鳳華老師在EBSD實驗過程中給予了切實的指導;中國科學院大學研究生胡俊成在野外工作中給予了大量幫助;審稿人提出了寶貴的修改意見和建議;作者在此向他們致以誠摯的謝意。

    在李繼亮先生逝世一周年之際,謹以此文寄托哀思,先生千古!

    猜你喜歡
    東段大水渦度
    英屬印度“科學邊疆”擴張戰(zhàn)略與中印邊界問題東段的形成
    嘩啦嘩啦大水來啦!
    幼兒畫刊(2019年3期)2019-11-04 02:41:36
    南海相對渦度的時空變化特征分析*
    色爾騰山山前斷裂東段與西段轉(zhuǎn)折處構(gòu)造演化特征
    斜壓渦度的變化與臺風暴雨的關系研究
    奇怪森林
    2017年11月16日至17日博州降雪天氣分析
    科技資訊(2018年30期)2018-03-11 09:30:58
    變大水“漫灌”為精準“滴灌”
    青海湖高原
    中國詩歌(2015年3期)2015-11-15 00:33:15
    一次東北南部特大暴雪的動力機制分析
    国产亚洲欧美在线一区二区| 日本一区二区免费在线视频| 少妇 在线观看| 婷婷六月久久综合丁香| 精品国产亚洲在线| 久久精品国产综合久久久| 两个人看的免费小视频| aaaaa片日本免费| 18禁黄网站禁片午夜丰满| 欧美黑人巨大hd| 国产精品爽爽va在线观看网站 | 国产精品影院久久| 两性午夜刺激爽爽歪歪视频在线观看 | 日韩欧美国产一区二区入口| 国产精品亚洲一级av第二区| 精品国产乱子伦一区二区三区| 正在播放国产对白刺激| 狂野欧美激情性xxxx| 在线观看免费午夜福利视频| 日韩大码丰满熟妇| 十八禁人妻一区二区| 日韩高清综合在线| av福利片在线| tocl精华| 人人妻人人澡人人看| 欧美人与性动交α欧美精品济南到| 亚洲精华国产精华精| 欧美激情高清一区二区三区| 日韩视频一区二区在线观看| 琪琪午夜伦伦电影理论片6080| 男女床上黄色一级片免费看| 在线永久观看黄色视频| 午夜福利一区二区在线看| 亚洲精品在线观看二区| 欧美在线黄色| 男女下面进入的视频免费午夜 | 午夜亚洲福利在线播放| 亚洲国产精品成人综合色| av视频在线观看入口| av欧美777| 久久久久久九九精品二区国产 | 日韩欧美一区二区三区在线观看| 视频区欧美日本亚洲| 草草在线视频免费看| 亚洲国产看品久久| 国产精品久久久久久精品电影 | 俄罗斯特黄特色一大片| 国产高清有码在线观看视频 | 国产亚洲av高清不卡| 欧美三级亚洲精品| 欧美三级亚洲精品| 欧美亚洲日本最大视频资源| 日本黄色视频三级网站网址| 国产午夜福利久久久久久| 在线观看www视频免费| 久久久久久九九精品二区国产 | 99国产极品粉嫩在线观看| 999久久久国产精品视频| 搡老熟女国产l中国老女人| 成年版毛片免费区| 999久久久精品免费观看国产| 天天躁狠狠躁夜夜躁狠狠躁| 欧美黑人精品巨大| 男人舔女人下体高潮全视频| 国产又色又爽无遮挡免费看| 婷婷精品国产亚洲av在线| 男人舔女人的私密视频| 国产亚洲精品一区二区www| 亚洲av成人一区二区三| 热re99久久国产66热| 麻豆一二三区av精品| 757午夜福利合集在线观看| 一级a爱视频在线免费观看| 午夜免费成人在线视频| 2021天堂中文幕一二区在线观 | 久99久视频精品免费| 国产精品二区激情视频| 精品久久久久久久毛片微露脸| 亚洲精品久久国产高清桃花| 一级作爱视频免费观看| 久久精品国产99精品国产亚洲性色| 麻豆国产av国片精品| 亚洲欧美激情综合另类| 国产成人啪精品午夜网站| 成人18禁高潮啪啪吃奶动态图| 欧美最黄视频在线播放免费| 91九色精品人成在线观看| 在线观看舔阴道视频| 在线十欧美十亚洲十日本专区| 久久国产精品影院| 久久久久久人人人人人| 国产三级在线视频| 男女下面进入的视频免费午夜 | 99热只有精品国产| 欧美午夜高清在线| 国产蜜桃级精品一区二区三区| 亚洲精品久久成人aⅴ小说| 国产又色又爽无遮挡免费看| 国产男靠女视频免费网站| 精品一区二区三区视频在线观看免费| 久久 成人 亚洲| 精品日产1卡2卡| 久久香蕉激情| 女警被强在线播放| 亚洲人成电影免费在线| 一级作爱视频免费观看| 国产麻豆成人av免费视频| 狂野欧美激情性xxxx| 日韩有码中文字幕| 国产极品粉嫩免费观看在线| 国产国语露脸激情在线看| 老司机午夜十八禁免费视频| 亚洲第一欧美日韩一区二区三区| 免费看日本二区| 精品第一国产精品| 日日摸夜夜添夜夜添小说| 中文资源天堂在线| 曰老女人黄片| 国产午夜福利久久久久久| 悠悠久久av| 亚洲精品一区av在线观看| 女人高潮潮喷娇喘18禁视频| 亚洲第一av免费看| 久久精品影院6| 日本免费一区二区三区高清不卡| 黄网站色视频无遮挡免费观看| 国产亚洲精品久久久久久毛片| 神马国产精品三级电影在线观看 | 熟女电影av网| 亚洲精品国产区一区二| 黑人巨大精品欧美一区二区mp4| 亚洲一区二区三区色噜噜| 一进一出抽搐gif免费好疼| 波多野结衣av一区二区av| 午夜a级毛片| 亚洲精品久久成人aⅴ小说| 91麻豆av在线| 国产熟女午夜一区二区三区| 国产极品粉嫩免费观看在线| 女警被强在线播放| 热99re8久久精品国产| 成年女人毛片免费观看观看9| 精品一区二区三区av网在线观看| 久久国产亚洲av麻豆专区| 岛国在线观看网站| 麻豆成人av在线观看| 国产一级毛片七仙女欲春2 | 精品乱码久久久久久99久播| av有码第一页| 男女之事视频高清在线观看| 久久中文字幕一级| 久久中文字幕人妻熟女| 国产成人精品久久二区二区91| 免费在线观看完整版高清| 老汉色av国产亚洲站长工具| 性色av乱码一区二区三区2| 日本 欧美在线| tocl精华| 首页视频小说图片口味搜索| 日韩欧美国产在线观看| 在线国产一区二区在线| 欧美不卡视频在线免费观看 | 午夜两性在线视频| 亚洲专区国产一区二区| 精品不卡国产一区二区三区| 亚洲熟女毛片儿| 久久欧美精品欧美久久欧美| 国产真人三级小视频在线观看| 少妇被粗大的猛进出69影院| 国产欧美日韩精品亚洲av| 久久久国产精品麻豆| 欧美黑人巨大hd| 99热6这里只有精品| 日韩欧美在线二视频| 中文资源天堂在线| 一进一出抽搐gif免费好疼| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲中文av在线| 青草久久国产| 国产单亲对白刺激| 亚洲人成网站高清观看| 国产激情欧美一区二区| 亚洲自偷自拍图片 自拍| 一区二区三区激情视频| 欧美一区二区精品小视频在线| 波多野结衣高清无吗| 欧美在线黄色| 亚洲久久久国产精品| 琪琪午夜伦伦电影理论片6080| 国产成人精品久久二区二区91| 青草久久国产| 久久精品国产99精品国产亚洲性色| 中文字幕最新亚洲高清| 淫妇啪啪啪对白视频| 又黄又粗又硬又大视频| 亚洲三区欧美一区| 最好的美女福利视频网| 午夜免费观看网址| 男女之事视频高清在线观看| 久久九九热精品免费| 欧美日韩亚洲国产一区二区在线观看| xxx96com| 欧美精品亚洲一区二区| 12—13女人毛片做爰片一| 国产又爽黄色视频| 国产精品免费一区二区三区在线| 黄片小视频在线播放| 叶爱在线成人免费视频播放| 亚洲精品久久国产高清桃花| aaaaa片日本免费| 中文字幕最新亚洲高清| 韩国av一区二区三区四区| 国产精品98久久久久久宅男小说| 午夜视频精品福利| 天天躁狠狠躁夜夜躁狠狠躁| 91麻豆精品激情在线观看国产| 俄罗斯特黄特色一大片| 精品欧美国产一区二区三| 91麻豆av在线| 丁香六月欧美| 国语自产精品视频在线第100页| 12—13女人毛片做爰片一| 国产精品久久久久久精品电影 | 免费看日本二区| 亚洲在线自拍视频| 亚洲国产精品合色在线| 少妇粗大呻吟视频| 中文字幕人成人乱码亚洲影| 日韩三级视频一区二区三区| 最近最新免费中文字幕在线| 美女 人体艺术 gogo| 美女扒开内裤让男人捅视频| 亚洲 欧美一区二区三区| 美女国产高潮福利片在线看| 午夜视频精品福利| 久久99热这里只有精品18| 成在线人永久免费视频| 亚洲五月婷婷丁香| 女人被狂操c到高潮| 怎么达到女性高潮| 亚洲一区中文字幕在线| 大型黄色视频在线免费观看| 午夜福利一区二区在线看| 久久精品国产亚洲av高清一级| 欧美精品啪啪一区二区三区| 悠悠久久av| 日韩欧美一区二区三区在线观看| 国产精品精品国产色婷婷| 一本一本综合久久| 国内久久婷婷六月综合欲色啪| 欧美乱码精品一区二区三区| 久久香蕉精品热| 日韩精品免费视频一区二区三区| 女性生殖器流出的白浆| 青草久久国产| 国产精品一区二区免费欧美| 欧美日韩黄片免| 美女高潮喷水抽搐中文字幕| 欧美成人性av电影在线观看| 首页视频小说图片口味搜索| 国产伦在线观看视频一区| av中文乱码字幕在线| 美女高潮到喷水免费观看| 国产欧美日韩一区二区精品| 精品一区二区三区av网在线观看| 日韩中文字幕欧美一区二区| 国产欧美日韩一区二区三| 免费搜索国产男女视频| 欧美日韩亚洲国产一区二区在线观看| 少妇熟女aⅴ在线视频| 久久国产乱子伦精品免费另类| 久久久精品欧美日韩精品| 给我免费播放毛片高清在线观看| 国产久久久一区二区三区| 亚洲午夜精品一区,二区,三区| 国产又色又爽无遮挡免费看| 免费高清在线观看日韩| 精品熟女少妇八av免费久了| 亚洲人成网站高清观看| 黄网站色视频无遮挡免费观看| 久久天堂一区二区三区四区| 在线观看午夜福利视频| 亚洲人成网站高清观看| 日韩中文字幕欧美一区二区| 亚洲男人天堂网一区| 亚洲成av人片免费观看| 人人澡人人妻人| 动漫黄色视频在线观看| 日韩欧美三级三区| 久久久国产成人免费| 我的亚洲天堂| 身体一侧抽搐| 国产黄a三级三级三级人| 欧美色视频一区免费| 亚洲 欧美一区二区三区| 亚洲中文字幕一区二区三区有码在线看 | 国产精品av久久久久免费| 成年免费大片在线观看| 国产激情偷乱视频一区二区| 不卡一级毛片| 变态另类丝袜制服| 一本久久中文字幕| 淫妇啪啪啪对白视频| 亚洲三区欧美一区| 久久精品国产清高在天天线| 午夜福利在线观看吧| 亚洲avbb在线观看| 成人国产综合亚洲| 男人舔女人下体高潮全视频| 国产成人精品无人区| 欧美日韩亚洲国产一区二区在线观看| 久久精品亚洲精品国产色婷小说| 亚洲欧美日韩无卡精品| 欧美日韩亚洲综合一区二区三区_| 99久久精品国产亚洲精品| 国产黄片美女视频| 美女高潮喷水抽搐中文字幕| 亚洲一卡2卡3卡4卡5卡精品中文| 九色国产91popny在线| 在线观看免费视频日本深夜| 亚洲精品粉嫩美女一区| 日本一区二区免费在线视频| 波多野结衣巨乳人妻| 91成人精品电影| 欧美激情高清一区二区三区| 在线av久久热| 欧美日韩瑟瑟在线播放| 亚洲欧美日韩无卡精品| 亚洲精品中文字幕一二三四区| 欧美中文综合在线视频| 国产精品一区二区精品视频观看| 黄色女人牲交| 亚洲成av人片免费观看| 国产1区2区3区精品| 午夜视频精品福利| 亚洲无线在线观看| 久99久视频精品免费| 看免费av毛片| 午夜视频精品福利| 亚洲精品国产精品久久久不卡| 精品欧美国产一区二区三| 成年女人毛片免费观看观看9| 国产伦一二天堂av在线观看| av中文乱码字幕在线| 性色av乱码一区二区三区2| www.熟女人妻精品国产| 午夜视频精品福利| 亚洲aⅴ乱码一区二区在线播放 | 视频在线观看一区二区三区| 免费观看精品视频网站| 一级a爱片免费观看的视频| 黄片大片在线免费观看| 久久久久免费精品人妻一区二区 | 1024视频免费在线观看| 免费观看精品视频网站| 国产主播在线观看一区二区| 免费搜索国产男女视频| 91在线观看av| 免费高清视频大片| 无人区码免费观看不卡| av电影中文网址| 午夜福利欧美成人| 中国美女看黄片| 91老司机精品| 久久精品人妻少妇| 日韩欧美免费精品| www.精华液| 夜夜看夜夜爽夜夜摸| 亚洲精品美女久久av网站| 国产成人一区二区三区免费视频网站| 欧美乱码精品一区二区三区| 成人三级黄色视频| 色精品久久人妻99蜜桃| 看片在线看免费视频| 12—13女人毛片做爰片一| 极品教师在线免费播放| 亚洲avbb在线观看| 女同久久另类99精品国产91| 国产三级黄色录像| 亚洲精品国产一区二区精华液| 免费在线观看影片大全网站| 亚洲人成77777在线视频| 欧美成人免费av一区二区三区| 精品国产国语对白av| 国产爱豆传媒在线观看 | 欧美在线黄色| 91av网站免费观看| 亚洲成a人片在线一区二区| 少妇熟女aⅴ在线视频| 国内揄拍国产精品人妻在线 | 最近最新中文字幕大全电影3 | 精品乱码久久久久久99久播| 精品久久蜜臀av无| 此物有八面人人有两片| 亚洲精华国产精华精| 精品卡一卡二卡四卡免费| 国产高清激情床上av| 日本一区二区免费在线视频| 国产成人系列免费观看| 国产三级黄色录像| 日本撒尿小便嘘嘘汇集6| av免费在线观看网站| 亚洲无线在线观看| 91大片在线观看| 成人国产一区最新在线观看| 久久精品国产亚洲av高清一级| 国产人伦9x9x在线观看| 97碰自拍视频| 久久精品国产亚洲av高清一级| 精品第一国产精品| 97人妻精品一区二区三区麻豆 | 亚洲天堂国产精品一区在线| 女人高潮潮喷娇喘18禁视频| 欧美在线一区亚洲| 最近最新中文字幕大全电影3 | 女生性感内裤真人,穿戴方法视频| 欧美日韩亚洲国产一区二区在线观看| 日韩免费av在线播放| 精品久久久久久久久久久久久 | 亚洲av成人一区二区三| 老司机靠b影院| 两个人视频免费观看高清| 亚洲精品国产一区二区精华液| 亚洲全国av大片| 黄色视频不卡| 亚洲第一电影网av| 又紧又爽又黄一区二区| 成在线人永久免费视频| 亚洲天堂国产精品一区在线| 亚洲成人精品中文字幕电影| 色播亚洲综合网| www.精华液| 日韩欧美国产一区二区入口| 国产精品免费一区二区三区在线| 久久天躁狠狠躁夜夜2o2o| 亚洲 欧美一区二区三区| 波多野结衣av一区二区av| 这个男人来自地球电影免费观看| 久久 成人 亚洲| 美女午夜性视频免费| 久久久精品国产亚洲av高清涩受| 日韩精品中文字幕看吧| 十八禁网站免费在线| 国内精品久久久久久久电影| 侵犯人妻中文字幕一二三四区| 欧美日韩中文字幕国产精品一区二区三区| 亚洲一区高清亚洲精品| 老司机靠b影院| 在线观看免费日韩欧美大片| 亚洲国产精品合色在线| videosex国产| 日韩欧美在线二视频| 每晚都被弄得嗷嗷叫到高潮| 欧美日韩乱码在线| 欧美久久黑人一区二区| 国产亚洲精品第一综合不卡| 午夜福利在线观看吧| 老司机深夜福利视频在线观看| 99在线人妻在线中文字幕| 99久久综合精品五月天人人| 亚洲电影在线观看av| 人成视频在线观看免费观看| 夜夜躁狠狠躁天天躁| 久久亚洲精品不卡| 国产欧美日韩一区二区精品| 欧美国产精品va在线观看不卡| 十八禁网站免费在线| 91麻豆av在线| 国产精品久久久久久亚洲av鲁大| 亚洲aⅴ乱码一区二区在线播放 | 啦啦啦韩国在线观看视频| 成年版毛片免费区| 午夜福利免费观看在线| 成人18禁高潮啪啪吃奶动态图| 自线自在国产av| 欧美成人一区二区免费高清观看 | 啦啦啦 在线观看视频| 欧美中文综合在线视频| 91字幕亚洲| 丝袜人妻中文字幕| 国产熟女xx| 波多野结衣高清作品| 人妻久久中文字幕网| 免费在线观看影片大全网站| 身体一侧抽搐| 欧美zozozo另类| 久久久久久久久中文| 精品国产美女av久久久久小说| 久久久精品国产亚洲av高清涩受| 亚洲专区国产一区二区| 女性生殖器流出的白浆| 丰满人妻熟妇乱又伦精品不卡| 女人高潮潮喷娇喘18禁视频| 国产亚洲av嫩草精品影院| 欧美国产精品va在线观看不卡| 天堂√8在线中文| 人妻久久中文字幕网| 伊人久久大香线蕉亚洲五| 男人的好看免费观看在线视频 | 黄色毛片三级朝国网站| 久热这里只有精品99| 夜夜夜夜夜久久久久| 亚洲成人久久性| 夜夜爽天天搞| 欧美乱码精品一区二区三区| 岛国视频午夜一区免费看| 国产日本99.免费观看| 韩国av一区二区三区四区| 十八禁网站免费在线| 国产久久久一区二区三区| 国产三级在线视频| 一级作爱视频免费观看| 麻豆成人av在线观看| 国产高清有码在线观看视频 | 在线观看免费午夜福利视频| 亚洲国产高清在线一区二区三 | 久久人妻av系列| 久久精品国产亚洲av香蕉五月| 国产成人欧美在线观看| 好男人电影高清在线观看| 久久香蕉激情| 熟女电影av网| av超薄肉色丝袜交足视频| 日本熟妇午夜| 久久人妻福利社区极品人妻图片| 精品国内亚洲2022精品成人| 亚洲av美国av| 成人18禁高潮啪啪吃奶动态图| 男女下面进入的视频免费午夜 | 好看av亚洲va欧美ⅴa在| 久久久久久久午夜电影| 淫妇啪啪啪对白视频| 99国产精品99久久久久| 岛国视频午夜一区免费看| svipshipincom国产片| 久久婷婷人人爽人人干人人爱| 国产久久久一区二区三区| 国产私拍福利视频在线观看| 91大片在线观看| 亚洲av美国av| 免费观看精品视频网站| 欧美黑人精品巨大| 老司机福利观看| 999久久久国产精品视频| 国产精品美女特级片免费视频播放器 | 天堂动漫精品| 亚洲天堂国产精品一区在线| 精品国产美女av久久久久小说| 国产精品久久电影中文字幕| 午夜福利欧美成人| 亚洲avbb在线观看| 日韩 欧美 亚洲 中文字幕| 婷婷丁香在线五月| 无遮挡黄片免费观看| 久久久精品国产亚洲av高清涩受| 日韩精品青青久久久久久| 99久久国产精品久久久| 亚洲色图 男人天堂 中文字幕| 99久久久亚洲精品蜜臀av| 丝袜人妻中文字幕| 国产又黄又爽又无遮挡在线| 久久精品人妻少妇| 亚洲中文av在线| 97超级碰碰碰精品色视频在线观看| 日本成人三级电影网站| 最近最新中文字幕大全免费视频| 国产欧美日韩一区二区精品| 成年版毛片免费区| 欧美成人性av电影在线观看| 99国产精品99久久久久| 女生性感内裤真人,穿戴方法视频| 日韩大尺度精品在线看网址| 日韩 欧美 亚洲 中文字幕| 一进一出抽搐动态| 妹子高潮喷水视频| 狠狠狠狠99中文字幕| 色综合婷婷激情| 久久中文字幕人妻熟女| 亚洲avbb在线观看| 丁香欧美五月| 一卡2卡三卡四卡精品乱码亚洲| 亚洲精品在线观看二区| 久久久精品国产亚洲av高清涩受| 好男人电影高清在线观看| 人人妻人人澡人人看| 欧美日韩黄片免| 最近在线观看免费完整版| 亚洲精品国产区一区二| 国产欧美日韩一区二区精品| 久久精品91蜜桃| 最近最新中文字幕大全免费视频| 午夜日韩欧美国产| av在线天堂中文字幕| 成人永久免费在线观看视频| 久久久久久久午夜电影| 热re99久久国产66热| 少妇被粗大的猛进出69影院| 老司机福利观看| 18禁国产床啪视频网站| 日本三级黄在线观看| 无遮挡黄片免费观看| 午夜免费激情av| 亚洲国产精品合色在线| 一个人免费在线观看的高清视频| 欧美日韩亚洲国产一区二区在线观看| 男女午夜视频在线观看| 非洲黑人性xxxx精品又粗又长| 好男人电影高清在线观看| x7x7x7水蜜桃| 午夜福利18| 欧美人与性动交α欧美精品济南到| 国产一区二区三区视频了| 一级作爱视频免费观看| 国内少妇人妻偷人精品xxx网站 | 日本一本二区三区精品|