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

    膠東郭家?guī)X巖體固結(jié)冷卻軌跡與隆升剝蝕*

    2015-04-13 04:25:02豆敬兆付順張華鋒
    巖石學(xué)報(bào) 2015年8期
    關(guān)鍵詞:隆升角閃石斜長(zhǎng)石

    豆敬兆 付順 張華鋒**

    DOU JingZhao1,F(xiàn)U Shun2 and ZHANG HuaFeng1**

    1. 中國(guó)地質(zhì)大學(xué)地球科學(xué)與資源學(xué)院,北京 100083

    2. 成都理工大學(xué),數(shù)學(xué)地質(zhì)四川省重點(diǎn)實(shí)驗(yàn)室,成都 610059

    1. School of Earth Sciences and Resources,China University of Geosciences,Beijing 100083,China

    2. Sichuan Mathematical Geology Key Laboratory,Chengdu University of Technology,Chengdu 610059,China

    2014-09-07 收稿,2015-01-31 改回.

    1 引言

    研究地殼隆升剝蝕過(guò)程對(duì)了解造山帶的演化及找礦具有重要意義。地殼隆升剝蝕的研究方法主要有熱年代學(xué)、盆地碎屑鋯石年齡統(tǒng)計(jì)分析及礦物壓力計(jì)。不同礦物的同位素體系封閉溫度不同而導(dǎo)致記錄的同位素年齡也存在差異,這一特征能夠反映地質(zhì)體的冷卻歷史(Hart,1964)。Dodson(1973)在此基礎(chǔ)上提出了礦物封閉溫度的概念,建立了礦物同位素冷卻年齡理論,并得到了良好應(yīng)用(Chen et al.,2013;Lin et al.,2011)。此方法前提是隆升速率等于地表剝蝕速率。盆地碎屑鋯石年齡統(tǒng)計(jì)分析是通過(guò)對(duì)沉積物中碎屑鋯石的年齡與區(qū)域地質(zhì)體同位素年齡進(jìn)行對(duì)比,通過(guò)同位素示蹤判斷沉積物源區(qū),結(jié)合造山時(shí)限,進(jìn)而反演造山帶的隆升剝蝕過(guò)程(Yang et al.,2006,2014c)。礦物壓力計(jì)則是利用巖體中某些特定礦物的化學(xué)成分估算巖體侵位固結(jié)壓力,該方法簡(jiǎn)單有效而得到應(yīng)用(馬昌前等,1995;Sial et al.,1999;Helmy et al.,2004;Zhang et al.,2006)。近年來(lái)鋯石、磷灰石(U-Th)/He 熱年代學(xué)開(kāi)始在礦床剝露史示蹤中得到應(yīng)用(Li et al.,2014;Liu et al.,2014)。本文采用礦物溫壓計(jì)對(duì)膠東郭家?guī)X巖體的冷卻固結(jié)條件進(jìn)行分析,結(jié)合前人發(fā)表的同位素?cái)?shù)據(jù),探討其冷卻抬升歷史。

    郭家?guī)X巖體位于膠東半島西北部,玲瓏金礦區(qū)及招掖成礦帶的北部。該巖體侵入膠東巖群及侏羅紀(jì)玲瓏巖體,形成時(shí)間在130 ~126Ma(關(guān)康等,1998;Wang et al.,1998,2014;Yang et al.,2012,2014b)。前人通過(guò)礦物學(xué)及巖石學(xué)的研究,提出巖漿混合(陳光遠(yuǎn)等,1993;曲曉明和王鶴年,1997;Hou et al.,2007)、高Ba-Sr 花崗巖(Wang et al.,2014)及下地殼鎂鐵質(zhì)巖石脫水部分熔融(楊進(jìn)輝等,2003)等不同認(rèn)識(shí)。

    對(duì)于郭家?guī)X巖體形成物理化學(xué)條件,前人研究頗少。林文蔚和殷秀蘭(1998)通過(guò)Q-Ab-Or-H2O 體系實(shí)驗(yàn)相圖、礦物溫壓計(jì)等對(duì)郭家?guī)X巖體形成的溫度、壓力及氧逸度等物理化學(xué)條件進(jìn)行了研究,認(rèn)為郭家?guī)X巖體形成過(guò)程中富水的巖漿流體捕獲了金,熱液的濃縮作用促進(jìn)了金礦的形成。陳光遠(yuǎn)等(1993)通過(guò)詳細(xì)的成因礦物學(xué)提出,由東向西,角閃石由短柱狀至長(zhǎng)柱狀的變化反映的是郭家?guī)X巖體由東至西巖體侵位逐漸變淺。而角閃石-斜長(zhǎng)石溫壓計(jì)獲得的結(jié)果則顯示由東向西溫度及壓力逐漸升高,這種現(xiàn)象被認(rèn)為是膠東西部差異隆升的結(jié)果(陸麗娜等,2011)。盡管前人對(duì)郭家?guī)X巖體形成的物理化學(xué)條件有過(guò)相應(yīng)研究,并對(duì)膠東西部100Ma 年以來(lái)地殼隆升剝蝕速率進(jìn)行過(guò)研究(柳振江等,2010),但對(duì)于郭家?guī)X巖體形成后的冷卻史、隆升剝蝕及其機(jī)制缺乏詳細(xì)研究,對(duì)巖體形成與金礦的關(guān)系也存在爭(zhēng)議。本文通過(guò)郭家?guī)X巖體物理化學(xué)條件的詳細(xì)研究,結(jié)合已有資料針對(duì)上述問(wèn)題進(jìn)行深入探討。

    2 區(qū)域地質(zhì)和巖體地質(zhì)

    膠東半島位于中國(guó)東部,西界為郯廬斷裂,南側(cè)為蘇魯-大別造山帶(圖1a),由經(jīng)歷超高壓變質(zhì)作用的揚(yáng)子陸塊和具有華北克拉通屬性的膠北隆起組成,是金礦床的主要富集區(qū)(Li and Santosh,2014;Fan et al.,2014)。兩大構(gòu)造單元?jiǎng)t以桃村和米山斷裂圍限的昆崳山混雜帶分隔(Zhai et al.,2000)。區(qū)域上出露的地層主要有太古宙膠東巖群、古元古界荊山群和粉子山群變質(zhì)巖以及少量震旦亞界蓬萊群(陳光遠(yuǎn)等,1993)。巖漿活動(dòng)主要為中生代花崗巖類,伴生少量中基性脈巖,出露的花崗巖體主要為侏羅紀(jì)玲瓏、灤家河及白堊紀(jì)郭家?guī)X和艾山巖體(徐金方和沈步云,1989;孫豐月等,1995)。玲瓏巖體位于新城-焦家斷裂帶和招遠(yuǎn)-平度斷裂帶之間,呈NNE 帶狀展布,巖性為片麻狀黑云母花崗巖和花崗閃長(zhǎng)巖,區(qū)域上侵入太古宙膠東巖群;鋯石U-Pb 年齡顯示形成時(shí)代在163 ~149Ma(Jiang et al.,2012;Yang et al.,2014b),是膠西北金礦床主要的容礦圍巖。賦存于玲瓏巖體內(nèi)的金礦類型主要為石英脈型金礦(楊立強(qiáng)等,2014;Wen et al.,2014;Song et al.,2015)。研究區(qū)內(nèi)發(fā)育多條NE、NNE 向斷裂帶,是中生代金礦的主要控礦斷裂(鄧軍等,1996,1999;呂古賢等,2006;宋明春等,2009)。

    郭家?guī)X巖體是膠東西北部金礦床的主要控礦圍巖,呈NE 向帶狀分布(圖1b),自東向西依次為郭家?guī)X、叢家、北截、上莊、三山島等巖體(陳光遠(yuǎn)等,1993)。郭家?guī)X巖體侵入玲瓏巖體,并在東部郭家?guī)X地區(qū)被艾山巖體侵入。賦存于郭家?guī)X巖體中的金礦為破碎蝕變巖型金礦,其礦化時(shí)間在126 ~120Ma(Wang et al.,1998;Qiu et al.,2002;Hu et al.,2013)。巖體內(nèi)廣泛發(fā)育韌性剪切變形,同時(shí)可見(jiàn)大量暗色微粒包體以及暗色礦物為主的流動(dòng)狀條帶(圖2a)。

    本次野外樣品采自郭家?guī)X63(37°29'18″N,120°36'01″E),叢家15B(37°32'25″N,120°33'15″E),北截11B(37°29'17″N,120°23'28″E)及上莊07(37°24'54″N,120°11'40″E)等巖體。巖性主要為花崗閃長(zhǎng)巖、二長(zhǎng)花崗巖及石英閃長(zhǎng)巖,普遍為似斑狀結(jié)構(gòu),塊狀構(gòu)造,斑晶主要為斜長(zhǎng)石(40% ~50%),鉀長(zhǎng)石(25% ~30%),石英(20% ~25%),角閃石(4% ±),黑云母(5% ±),另有副礦物榍石、褐簾石、綠簾石、磷灰石及磁鐵礦等(圖2d,f)。

    斜長(zhǎng)石自形-半自形,普遍發(fā)育聚片雙晶,大小約0.5 ×1mm,并見(jiàn)蠕蟲(chóng)結(jié)構(gòu)(圖2d)。斜長(zhǎng)石存于礦物顆粒之間及鉀長(zhǎng)石中,表明鉀長(zhǎng)石晚于斜長(zhǎng)石結(jié)晶。鉀長(zhǎng)石,自形-半自形,斑晶顆粒較大,約4 ×3cm,基質(zhì)鉀長(zhǎng)石中發(fā)育自形斜長(zhǎng)石及黑云母等(圖2b,e)。

    角閃石較自形,斑晶大小約1.5 ×2mm,不同的巖體中角閃石含量略有變化,約4% ~5%。

    圖1 膠東西部區(qū)域地質(zhì)簡(jiǎn)圖(a,據(jù)Zhao et al.,1999 修改;b,據(jù)山東省地勘局,1999①山東省地勘局.1999.1∶150 萬(wàn)山東省地質(zhì)圖修改)Fig.1 Geological sketch map of the Guojialing granodiorites in western Jiaodong Peninsula,China (a,after Zhao et al.,1999)

    黑云母形態(tài)上明顯的分為兩類,一類自形,存于鉀長(zhǎng)石中及礦物顆粒間;另一類彎曲變形,發(fā)育于礦物顆粒間,可能受后期的構(gòu)造作用所致(圖2c,d)。兩類黑云母多色性較弱,淺綠或淺褐色-褐色,局部綠泥石化,并沿解理析出磁鐵礦。石英他形充填于礦物顆粒間,見(jiàn)波狀消光。褐簾石自形柱狀,邊緣為綠簾石,并顯示出弱的成分環(huán)帶,無(wú)次生蝕變現(xiàn)象,表明為巖漿成因(圖2f)。

    巖相學(xué)特征顯示,基質(zhì)中角閃石及斜長(zhǎng)石自形,鉀長(zhǎng)石自形-半自形,且礦物內(nèi)部包含有自形的黑云母及斜長(zhǎng)石,說(shuō)明基質(zhì)中的鉀長(zhǎng)石晚于斜長(zhǎng)石及黑云母結(jié)晶,黑云母發(fā)育于自形板狀的斜長(zhǎng)石礦物顆粒之間,暗示黑云母結(jié)晶于斜長(zhǎng)石之后,而他形的石英則表明其最晚結(jié)晶。

    3 礦物化學(xué)

    我們對(duì)黑云母、角閃石、鉀長(zhǎng)石及斜長(zhǎng)石進(jìn)行礦物化學(xué)分析,采用角閃石全鋁壓力計(jì)(Anderson and Smith,1995),獲得巖體侵位壓力;利用角閃石-斜長(zhǎng)石溫度計(jì)(Holland and Blundy,1994)計(jì)算其平衡溫度;采用黑云母-角閃石溫壓計(jì)(Симонова,1979)、黑云母溫度計(jì)(Henry et al.,2005)及二長(zhǎng)石溫度計(jì)(Putirka,2008)進(jìn)行礦物平衡溫度估算。在中國(guó)地質(zhì)科學(xué)院電子探針實(shí)驗(yàn)室JXA-8230 型電子探針測(cè)試儀器上完成,加速電壓15kV,電流20nA,束斑直徑5μm。測(cè)試結(jié)果見(jiàn)表1。

    3.1 黑云母

    礦物化學(xué)分析顯示,兩類黑云母在化學(xué)成分上基本相同。黑云母-角閃石平衡溫壓反映的是共生的角閃石和黑云母達(dá)到平衡時(shí)的溫壓條件,故邊部成分更接近平衡時(shí)的狀態(tài)。為采用黑云母-角閃石平衡溫壓計(jì),我們選擇黑云母干凈且無(wú)后期蝕變的邊部進(jìn)行礦物化學(xué)的分析。其氧化物含量SiO2=34.94% ~37.69%,F(xiàn)eO =15.70% ~17.92%,MgO=12.01% ~14.39%,TiO2=1.81% ~2.43%,Mg/(Mg+Fe)=0.55 ~0.63;黑云母分類圖解中屬于鎂黑云母(圖3a),與前人的結(jié)果一致(陳光遠(yuǎn)等,1993)。在成因分類中,落入再平衡區(qū)域(圖3c)。黑云母鏡下多色性較弱,呈淺綠色,并見(jiàn)含鈦磁鐵礦沿黑云母解理及黑云母與鉀長(zhǎng)石接觸邊緣發(fā)育,這與前人對(duì)于再平衡黑云母特征的描述一致,可能與巖漿后期流體作用有關(guān)(Nachit et al.,2005),而在形態(tài)上發(fā)生變形彎曲的黑云母,與巖體固結(jié)后的韌性剪切作用有關(guān)。在黑云母MgO-FeOT-Al2O3構(gòu)造環(huán)境判別圖中,顯示郭家?guī)X花崗閃長(zhǎng)巖為造山帶鈣堿性花崗巖(圖3d),與前人利用全巖化學(xué)成分獲得的結(jié)果一致(陳光遠(yuǎn)等,1993)。

    圖2 膠東郭家?guī)X巖體野外及巖相學(xué)特征(a、b)郭家?guī)X巖體中閃長(zhǎng)質(zhì)包體及鉀長(zhǎng)石斑晶;(c)自形及彎曲變形的兩類黑云母;(d、e)黑云母與石英、榍石、角閃石、斜長(zhǎng)石及鉀長(zhǎng)石共生;(f)自形褐簾石及粒狀角閃石. 礦物縮寫(xiě)據(jù)Whitney and Evans (2010):Am-角閃石;Bi-黑云母;Kfs-鉀長(zhǎng)石;Pl-斜長(zhǎng)石;Qz-石英;Spn-榍石;Aln-褐簾石;Ep-綠簾石Fig.2 Field and micrographical photos for the Guojialing granodiorites in Jiaodong Peninsula

    3.2 角閃石

    角閃石斑晶在花崗閃長(zhǎng)巖中以短柱狀、長(zhǎng)柱狀產(chǎn)出,半自形-自形。礦物化學(xué)顯示,角閃石SiO2= 42.54% ~44.76%,Al2O3=8.91% ~10.18%,MgO、CaO 含量較高,分別為10.64% ~11.99%和11.11% ~11.52%,其CaB>1.5,(Na+K)A<0.5 屬于鈣質(zhì)角閃石亞族中的鎂角閃石(圖3b)。

    3.3 長(zhǎng)石

    與角閃石共生的斜長(zhǎng)石礦物邊部化學(xué)分析表明,斜長(zhǎng)石SiO2= 62.13% ~62.53%,Na2O 含量較高,為8.97% ~9.51%,Al2O3=22.81% ~23.51%;Ab =0.79 ~0.82,An =0.17 ~0.21。

    鉀長(zhǎng)石產(chǎn)于斑晶及基質(zhì)中,電子探針結(jié)果顯示:SiO2=60.05% ~64.48%,Al2O3= 18.71% ~18.78%,K2O =13.69% ~14.26%;Ab=0.14 ~0.18,Or=0.82 ~0.86。

    4 巖漿冷卻固結(jié)條件

    4.1 溫壓條件

    根據(jù)前人數(shù)據(jù)(Zhang et al.,2010;陸麗娜等,2011),據(jù)Watson and Harrison (1983)提供的花崗巖鋯飽和溫度計(jì)算方法,獲得郭家?guī)X巖體的鋯飽和溫度為726 ~800℃,平均為755℃,代表巖漿侵位溫度。結(jié)合已有的主量數(shù)據(jù)(楊進(jìn)輝等,2003;Zhang et al.,2010;陸麗娜等,2011),采用QAb-Or 標(biāo)準(zhǔn)礦物的共結(jié)壓力圖解(Huang and Wyllie,1975),獲得巖體固結(jié)壓力集中在4 ~5kbar 之間。少量點(diǎn)偏離壓力線而處于5kbar 之下,這可能與巖體局部鉀長(zhǎng)石斑晶較多有關(guān),反映巖體局部的鉀長(zhǎng)石斑晶先結(jié)晶的特點(diǎn)(圖4a)。

    表1 膠東郭家?guī)X巖體黑云母、角閃石、斜長(zhǎng)石及鉀長(zhǎng)石礦物化學(xué)成分表(wt%)Table1 Chemicalcomposition ofbiotite, amphibole, plagioclaseand potassiumfeldsparfromGuojialingpluton in JiaodongPeninsula, China(wt%)

    圖3 郭家?guī)X巖體角閃石和黑云母礦物化學(xué)分類圖解(a,b)黑云母、角閃石分類圖解(底圖分別據(jù)Foster,1960;Leake et al. ,1997);(c)黑云母10TiO2-FeOT-MgO 成因類型圖解(Nachit et al. ,2005),其中A、B、C 分別代表原生、再平衡、新生黑云母區(qū)域;(d)黑云母MgO-FeOT-Al2O3 相關(guān)圖(Abdel-Rahman,1994),A、C、P 分別代表非造山帶堿性花崗巖、造山帶鈣堿性花崗巖、過(guò)鋁質(zhì)花崗巖. 空心圓圈為本人數(shù)據(jù),實(shí)心圓圈引自陳光遠(yuǎn)等(1993)Fig.3 Chemical composition classification diagrams of amphibole and biotite for the Guojialing granodiorites(a,b)Classification diagrams of biotite and calcic amphibole (after Foster,1960;Leake et al. ,1997,respectively);(c)10TiO2-FeOT-MgO diagram of biotite (after Nachit et al. ,2005),A,B and C correspond to the domains of primary magmatic biotites,reequilibrated biotites and the neoformed biotites,respectively;(d)MgO-FeOT-Al2O3 biotite discriminant diagram for biotite in anorogenic alkaline suites(field A),biotite in calcalkaline orogenic suites (field C)and biotite in peralumino us suites (field P)(after Abdel-Rahman,1994). Data represented by open circles from this paper and filled circles after Chen et al.(1993)

    Holland and Blundy (1994)基于淺閃石-透閃石(A)和淺閃石-鈉透閃石(B)反應(yīng)平衡的研究,分別提出兩個(gè)角閃石-斜長(zhǎng)石溫度計(jì)(A、B),其方程式如下:

    因?yàn)闇囟扔?jì)B 所獲得溫度能夠被其他礦物溫度計(jì)所再現(xiàn),故溫度計(jì)B 獲得的結(jié)果更為可信(Anderson,1996)。采用溫度計(jì)B 計(jì)算時(shí),需要代入相應(yīng)的壓力值,故先利用Schmit (1992)提到的角閃石全鋁壓力計(jì)獲得壓力結(jié)果,將獲得的壓力值代入溫度計(jì)B 中,得到平衡溫度為657 ~717℃(±40),集中在693 ~717℃(±40)。采用角閃石全鋁壓力計(jì)(Anderson and Smith,1995),并利用溫度計(jì)B 所獲得的結(jié)果對(duì)所獲壓力進(jìn)行校正,獲得的壓力為4.8 ±0.6kbar。

    圖4 郭家?guī)X巖體固結(jié)冷卻溫壓條件圖解(a)花崗巖Q-Ab-Or 共結(jié)壓力圖解(底圖據(jù)Huang and Wyllie,1975),數(shù)據(jù)引自Zhang et al. (2010);陸麗娜等(2011);(b、c)角閃石-黑云母礦物對(duì)Mg/(Ti+Mg+Fe+Mn)Bi-Mg/(Ti+Mg+Fe+Mn)Am溫度及Al/(Ti+Mg+Fe+Mn)Bi-Mg/(Ti+Mg+Fe+Mn)Am壓力圖解(底圖據(jù)Симонова,1979);(d)黑云母Ti-Mg/(Mg+Fe)溫度圖解(底圖據(jù)Henry et al. ,2005)Fig.4 Estimated pressures and temperatures for the Guojialing granodiorites(a)Q-Ab-Or diagram for the Guojialing granodiorites in Jiaodong Peninsula (after Huang and Wyllie,1975),data from Zhang et al. (2010),Lu et al. (2011);(b,c)temperature diagram of Mg/(Ti+Mg+Fe+Mn)Bi vs. Mg/(Ti+Mg+Fe+Mn)Am,and pressure diagram of Al/(Ti+Mg+Fe+Mn)Bi vs. Mg/(Ti+Mg+Fe+Mn)Am for amphibole and biotite (after Симонова,1979;(d)Ti-Mg/(Mg +Fe)geothermometry diagram for biotite (after Henry et al. ,2005)

    在角 閃 石-黑 云 母 壓 力 與 溫 度 圖 解 中(Симонова,1979),獲得二者平衡條件為4 ~4.4kbar,650 ~700℃(圖4b,c),與黑云母Ti-Mg/(Mg +Fe)投圖(Henry et al.,2005)獲得的溫度結(jié)果相近(圖4d)。利用Putirka (2008)提出的二長(zhǎng)石溫度計(jì),其方程式為:

    所獲得的平衡溫度略低,僅為570 ~580℃。

    4.2 氧逸度

    郭家?guī)X巖體中礦物組合多見(jiàn)石英+榍石+磁鐵礦組合,而該組合可以用于計(jì)算礦物結(jié)晶時(shí)的氧逸度條件(Wones,1989),其公式為:logfO2= -30930/T +14.98 +0.142(P -1)/T。根據(jù)角閃石-斜長(zhǎng)石溫度計(jì)所獲溫度(717 ~657℃)及角閃石全鋁壓力計(jì)獲得的壓力結(jié)果(4.8 ±0.6kbar),代入公式獲得的氧逸度為-18.3 ~-16.3。在云母氧緩沖劑圖解中(Wones and Eugster,1965),黑云母投點(diǎn)落在NB 氧緩沖劑反應(yīng)線附近(圖5a),結(jié)合黑云母結(jié)晶溫度,在不同緩沖劑的溫度-氧逸度圖解中(圖5c),黑云母結(jié)晶時(shí)的氧逸度條件大約在-18 ~-17,該結(jié)果與黑云母的Fe3+/(Fe3++Fe2+)與氧逸度關(guān)系圖中所獲得的氧逸度條件相近(圖5b)。根據(jù)角閃石-黑云母溫壓計(jì)估算結(jié)果,利用Wones(1989)氧逸度計(jì)算公式獲得的氧逸度值為-18.5 ~-16.8,與上述兩種方法所得結(jié)果一致。

    圖5 黑云母氧逸度及角閃石溫壓條件圖解(a)黑云母組成與氧緩沖劑相關(guān)圖(據(jù)Wones and Eugster,1965),MH,NB,F(xiàn)MQ 及WM 分別為磁鐵礦-赤鐵礦,自然鎳-綠鎳礦,鐵橄欖石-磁鐵礦-石英,方鐵礦-磁鐵礦緩沖劑組合;(b)黑云母100 ×Fe3+ /(Fe3+ +Fe2+)、結(jié)晶溫度與氧逸度相關(guān)圖(ΠeрчукЛЛ,1973);(c)不同氧緩沖劑溫度-氧逸度關(guān)系圖(Eugster and Wones,1962);(d)角閃石壓力計(jì)及角閃石-斜長(zhǎng)石溫度計(jì)計(jì)算所獲溫壓范圍Fig.5 Oxygen fugacity of biotite and P-T conditions for amphibole(a)the correlative diagram between biotite composition and oxygen buffer-reagent (after Wones and Eugster,1965);(b)100 × Fe3+ /(Fe3+ +Fe2+)crystal temperature and oxygen fugacity of biotite (ΠeрчукЛЛ et al. ,1973);(c)the different oxygen buffer-reagent temperature-oxygen fugacity diagram (after Eugster and Wones,1962);(d)the result of amphibole barometer and amphibole-plagioclase thermometer calculation

    5 討論

    5.1 巖體固結(jié)冷卻及隆升剝蝕

    花崗巖Q-Ab-Or 標(biāo)準(zhǔn)礦物結(jié)晶壓力圖解、角閃石全鋁壓力計(jì)及角閃石-黑云母壓力計(jì)所獲壓力基本一致,結(jié)果表明巖體侵位壓力在4 ~5kbar 之間;鏡下觀察到,巖石中發(fā)育大量自形的褐簾石和綠簾石(圖2f),且無(wú)溶蝕現(xiàn)象,表明這些綠簾石及褐簾石為巖漿成因,它的出現(xiàn)暗示巖體侵位壓力在3 ~5kbar(Schmidt and Thompson,1996)。綜上所述,我們認(rèn)為采用角閃石全鋁壓力計(jì)(Anderson et al.,1995)計(jì)算獲得的平衡壓力(4.8 ±0.6kbar)比較合理。按照上地殼平均比重(2.7g/cm3)和平均地壓梯度2.7km/kbar 進(jìn)行計(jì)算,郭家?guī)X巖體侵位深度約13 ±1.6km。

    郭家?guī)X巖體的鋯飽和溫度約在755℃,代表其巖漿侵位溫度。礦物平衡溫度則反映結(jié)晶冷卻過(guò)程。我們的結(jié)果表明,角閃石-斜長(zhǎng)石平衡溫度為657 ~717℃,角閃石-黑云母平衡溫度為650 ~700℃,而二長(zhǎng)石平衡溫度在570 ~580℃之間,這與巖相學(xué)觀察到的礦物結(jié)晶順序一致(基質(zhì)中礦物結(jié)晶順序?yàn)榻情W石-斜長(zhǎng)石-黑云母-鉀長(zhǎng)石-石英)?;|(zhì)中的鉀長(zhǎng)石盡管為晚期結(jié)晶礦物,但二長(zhǎng)石記錄的平衡溫度明顯較飽和水花崗巖固相線低了近70 ~100℃。這很可能與巖體侵位深度大,巖體固結(jié)冷卻后與圍巖溫度相差不大而相互平衡有關(guān)。我們所測(cè)二長(zhǎng)石均為基質(zhì)中新鮮潔凈顆粒,不存在后期蝕變改造可能。如果膠東地區(qū)早白堊世的地溫梯度按照35℃/km 和40℃/km 分別計(jì)算,郭家?guī)X巖體侵位固結(jié)后的圍巖溫度分別為455 ±56℃和520 ±64℃。后者與二長(zhǎng)石570 ~580℃的平衡溫度比較吻合。因此,我們將二長(zhǎng)石平衡溫度解釋為巖體與圍巖平衡的結(jié)果是合理的。另外,推算的本區(qū)早白堊世地溫梯度在40℃/km 左右比較合理。當(dāng)時(shí)的高地溫梯度可能與其構(gòu)造背景及區(qū)域巖漿作用關(guān)系密切。

    郭家?guī)X巖體鋯石年齡為130 ~126Ma(關(guān)康等,1998;Wang et al.,1998;Yang et al.,2012,2014b),而角閃石和黑云母的40Ar-39Ar 年齡分別為130 ±2Ma 和124 ±1Ma,二者理論上的封閉溫度分別在500 ~400℃和300 ~200℃之間(Reiners and Brandon,2006;Li et al.,2003)。它們的封閉溫度與礦物冷卻速率、深部地質(zhì)體隆升剝蝕速率關(guān)系密切,礦物冷卻速率和地質(zhì)體抬升速率越快,其封閉溫度則越高(Reiners and Brandon,2006)。如果我們假設(shè)郭家?guī)X巖體固結(jié)后發(fā)生快速冷卻抬升,角閃石與黑云母封閉溫度分別為最高的500℃和300℃。那么,巖體從130 ~126Ma(固結(jié)溫度650℃)至124Ma 的300℃,平均冷卻速率約60 ~175℃/Myr。假設(shè)角閃石冷卻到500℃獲得40Ar-39Ar 封閉年齡是130Ma,冷卻到黑云母記錄125 ~123Ma 時(shí)的封閉溫度為300℃,則巖體的平均冷卻速率為40 ~70℃/Myr。

    角閃石與黑云母40Ar-39Ar 開(kāi)始冷卻封閉的溫度不代表巖體當(dāng)時(shí)所處深度上圍巖溫度。但是,對(duì)于深成巖體根據(jù)地溫梯度可以估算礦物40Ar-39Ar 封閉時(shí)的最大深度。因此根據(jù)地溫梯度(40℃/km)計(jì)算,郭家?guī)X巖體與圍巖溫度平衡后經(jīng)歷了快速隆升過(guò)程。黑云母40Ar-39Ar 體系封閉溫度在300~200℃,對(duì)應(yīng)的最大深度約7.5 ~5km。由于郭家?guī)X巖體黑云母化學(xué)成分上顯示為再平衡結(jié)果,巖石學(xué)特征上也顯示大量黑云母經(jīng)歷過(guò)剪切變形作用。所以,我們認(rèn)為黑云母40Ar-39Ar冷卻封閉年齡很可能是區(qū)域韌性剪切作用后的冷卻年齡,即韌性剪切時(shí)代發(fā)生在126 ~124Ma,前人測(cè)得郭家?guī)X韌性剪切帶糜棱巖中角閃石和黑云母的40Ar-39Ar 年齡分別為124 ±1.8Ma、123 ±1.5Ma(Charles et al.,2013),也能說(shuō)明這一點(diǎn)。韌性剪切作用發(fā)生的深度至少在10km 左右。這意味著韌性剪切作用之后,郭家?guī)X巖體發(fā)生了快速隆升剝蝕。侵入郭家?guī)X巖體的艾山巖體鋯石年齡顯示為116 ±2Ma(Goss et al.,2010),侵位壓力在1kbar 左右(孫豐月等,1995),即深度僅2.7km 左右。所以,郭家?guī)X巖體黑云母40Ar-39Ar 記錄124 ±1Ma 時(shí)的合理深度約7.5 ~2.7km。本區(qū)金礦流體包裹體研究顯示金礦成礦壓力變化在1 ~3kbar,深度相當(dāng)于2.7~8.1km(Fan et al.,2003)。根據(jù)地溫梯度計(jì)算,成礦深度范圍內(nèi)地殼的溫度在100 ~320℃之間。這與郭家?guī)X巖體黑云母40Ar-39Ar 記錄124Ma 的封閉溫度吻合。

    綜上,郭家?guī)X巖體在130Ma 形成后到艾山巖體侵位(116Ma),地殼總體隆升剝蝕達(dá)10 ±1.6km 左右,其冷卻和隆升路徑見(jiàn)圖6。大量同位素年齡揭示出膠東金礦集中爆發(fā)在120 ±5Ma(Yang and Zhou,2000,2001;陳衍景等,2004;Li et al.,2008;于學(xué)峰等,2012;Guo et al.,2013;Hu et al.,2013),巖體快速隆升剝蝕發(fā)生在韌性剪切作用之后,約124 ~116Ma 之間。這說(shuō)明,成礦作用發(fā)生在韌性剪切作用后的地殼快速隆升剝蝕過(guò)程中。從時(shí)間和郭家?guī)X隆升剝蝕過(guò)程判斷,郭家?guī)X巖體在金礦爆發(fā)期間的快速隆升與金礦作用相關(guān)的巖漿-構(gòu)造熱事件緊密相關(guān)。

    圖6 郭家?guī)X巖體抬升軌跡圖其中TZr 為鋯飽和溫度(數(shù)據(jù)引自Zhang et al.,2010;陸麗娜等,2011);方法據(jù)Watson and Harrison (1983);壓力計(jì)算方法據(jù)花崗巖Q-Ab-Or 標(biāo)準(zhǔn)礦物結(jié)晶壓力圖解(數(shù)據(jù)引自楊進(jìn)輝等,2003;Zhang et al.,2010;陸麗娜等,2011). TAm-Pl為角閃石-斜長(zhǎng)石平衡溫壓,計(jì)算方法據(jù)Holland and Blundy (1994)及Anderson and Smith (1995).TAm-Bi為角閃石-黑云母平衡溫壓,計(jì)算方法據(jù)Симонова(1979).TPl-Kf為二長(zhǎng)平衡石溫度,計(jì)算方法據(jù)Putirka(2008). 角閃石40 Ar-39 Ar 年齡及封閉溫度引自Reiners and Brandon (2006). 黑云母40 Ar-39 Ar 年齡及封閉溫度引自Li et al. (2003)及Reiners and Brandon(2006);侵位深度按2.7km/kbar 的地壓梯度計(jì)算;艾山巖體侵位時(shí)的鋯石年齡及侵位壓力分別引自Goss et al. (2010)、孫豐月等(1995),對(duì)應(yīng)的郭家?guī)X巖體此時(shí)的溫度根據(jù)1kbar 的侵位壓力和2.7km/kbar 的地壓梯度計(jì)算,深度約在2.7km,以40℃/km 的地溫梯度計(jì)算,得到郭家?guī)X巖體與圍巖平衡溫度約100℃. 飽和水花崗巖固相線引自Luth et al. (1964). 巖體冷卻隆升路徑沿A-B-C 為理想結(jié)果,即巖體抬升過(guò)程中總是處于與圍巖溫度均衡狀態(tài)下,依據(jù)地溫梯度和礦物封閉溫度獲得的溫壓變化路徑;A-B’-C 則是極端的情況下,巖體韌性剪切作用后快速隆升到艾山巖體侵位時(shí)的1kbar 深度后,快速冷卻路徑. 因此,A-B-B’-C 組成的三角形區(qū)域內(nèi)為郭家?guī)X巖體冷卻抬升剝蝕的可能路徑區(qū)域Fig.6 Cooling and exhumation for the Guojialing granodiorites

    不同深度上形成的礦體目前處于同一地殼層次上不僅說(shuō)明成礦作用是多階段的,同時(shí)也進(jìn)一步佐證了成礦期間地殼發(fā)生過(guò)快速的隆升剝蝕作用。呂古賢(1997)曾對(duì)本區(qū)玲瓏和焦家金礦成礦深度和階段進(jìn)行過(guò)詳細(xì)研究,結(jié)果顯示本區(qū)金礦存在多階段且深度不同的礦體。盡管其采用的壓力校正方法獲得了相對(duì)較淺的成礦深度(約3.5 ~0.7km),但不同成礦階段的成礦深度存在明顯差距的特征有力的說(shuō)明了成礦期地殼發(fā)生過(guò)快速隆升剝蝕。

    5.2 地殼隆升機(jī)制與成礦作用關(guān)系

    郭家?guī)X巖體在10Myr 內(nèi)發(fā)生了大規(guī)模的隆升剝蝕(平均隆升剝蝕速率約1 ±0.1km/Myr),其機(jī)制可能與華北東部中生代巖石圈減薄導(dǎo)致的地殼淺部隆升、伸展及拆離有關(guān)(周新華等,2002;Zhai et al.,2002;Yang et al.,2003;Li and Santosh,2014)。膠西北金礦田多發(fā)育脆性斷裂帶中,其中賦存大量金礦體。它們的形成與地殼拆離和成礦流體上升并沿著這些脆性斷裂帶發(fā)生成礦作用緊密相關(guān)(范宏瑞等,2005;Yang et al.,2014a)。在伸展背景下,地殼的大規(guī)模水平拆離作用能夠造成上地殼巖石在極短的時(shí)間內(nèi)發(fā)生快速位移,并造成下覆地質(zhì)體快速減壓隆升和剝蝕作用。

    另外,本文對(duì)郭家?guī)X巖體的氧逸度估算結(jié)果顯示,氧逸度值-18 ~-16 之間,而氧逸度的高低是控制Au 沉淀的重要因素,氧逸度的降低更有利于Au 的沉淀(Williams-Jones et al.,2009;Li et al.,2013)。金礦沉淀大致發(fā)生在-30 ~-23 之間(Gibert et al.,1998;Li et al.,2013)。因此,我們推測(cè),在郭家?guī)X巖體固結(jié)晚期,排出的流體如果含礦則會(huì)在上覆地層中發(fā)生沉淀成礦。

    6 結(jié)論

    (1)郭家?guī)X巖體在約130Ma 侵位,侵位溫度在726 ~800℃,固結(jié)溫度650 ~700℃,就位壓力在4 ~5kbar;固結(jié)后巖體冷卻至570 ~580℃,此時(shí)二長(zhǎng)石溫度與圍巖溫度達(dá)到平衡;之后溫度降至500 ±50℃,隨后巖體冷卻到124Ma 的300±30℃左右;

    (2)郭家?guī)X巖體在126 ~116Ma 期間發(fā)生了快速隆升剝蝕,隆升剝蝕量在10km 左右;其發(fā)生時(shí)間與中國(guó)東部中生代巖石圈減薄的峰期時(shí)限耦合,暗示膠東西部地殼快速隆升剝蝕是深部地質(zhì)過(guò)程的淺部響應(yīng)。

    致謝 感謝中國(guó)地質(zhì)科學(xué)院地質(zhì)力學(xué)研究所呂古賢研究員和中國(guó)科學(xué)院地質(zhì)與地球物理所劉玄博士提出的修改建議;電子探針測(cè)試分析在中國(guó)地質(zhì)科學(xué)院電子探針實(shí)驗(yàn)室完成,對(duì)陳振宇老師的大力幫助和指導(dǎo),表示感謝。

    Abdel-Rahman AFM. 1994. Nature of biotites from alkaline,calcalkaline,and peraluminous magmas. Journal of Petrology,35(2):525 -541

    Anderson JL and Smith DR. 1995. The effects of temperature and fO2on the Al-in-hornblende barometer. American Mineralogist,80:549-559

    Anderson JL. 1996. Status of thermobarometry in granitic batholiths.Geological Society of America Special Papers,315:125 -138

    Charles N,Augier R,Gumiaux C,Monié P,Chen Y,F(xiàn)aure M and Zhu RX. 2013. Timing,duration and role of magmatism in wide rift systems:Insights from the Jiaodong Peninsula (China,East Asia).Gondwana Research,24(1):412 -428

    Chen G,Ding C,Xu LM,Zhang HR,Hu YX,Yang F,Li N,Mao XN and Li Y. 2013. Analysis on the thermal history and uplift process of Zijinshan intrusive complex in the eastern Ordos basin. Chinese Journal of Geophysics,56(1):78 -90

    Chen GY,Sun DS,Zhou XR,Shao W,Gong RT and Shao Y. 1993.Genetic Mineralogy and Gold Mineralization of Guojialing Granodiorite in Jiaodong Region. Wuhan: Chinese University Geosciences Press,8 -219 (in Chinese)

    Chen YJ,Pirajno F,Lai Y and Li C. 2004. Metallogenic time and tectonic setting of the Jiaodong gold province,eastern China. Acta Petrologica Sinica,20(4):907 - 922 (in Chinese with English abstract)

    Deng J,Xu SL,F(xiàn)ang Y,Zhou XQ and Wai L. 1996. The Tectonic Systems and Gold Metallogenic Dynamics in the Northwestern Jiaodong,China. Beijing:Geological Publishing House,1 -98 (in Chinese)

    Deng J,Zhai YS,Yang LQ,Xiao RG and Sun ZS. 1999. Tectonic evolution and dynamics of metallogenic system-an example from the gold ore deposits-concentrated area in Jiaodong,Shandong,China.Earth Science Frontiers,6(2):315 -323 (in Chinese with English abstract)

    Dodson MH. 1973. Closure temperature in cooling geochronological and petrological systems. Contributions to Mineralogy and Petrology,40(3):259 -274

    Eugster HP and Wones DR. 1962. Stability relations of the ferruginous biotite,annite. Journal of Petrology,3(1):82 -125

    Fan HR,Zhai MG,Xie YH and Yang JH. 2003. Ore-forming fluids associated with granite-hosted gold mineralization at the Sanshandao deposit,Jiaodong gold province,China. Mineralium Deposita,38(6):739 -750

    Fan HR,Hu FF,Yang JH,Shen K and Zhai MG. 2005. Fluid evolution and large-scale gold metallogeny during Mesozoic tectonic transition in the eastern Shandong province. Acta Petrologica Sinica,21(5):1317 -1328 (in Chinese with English abstract)

    Fan HR,Hu FF,Yang KF,Wen BJ and Liu X. 2014. Ore-forming fluids and ore genesis in the world-class Mesozoic gold province,Jiaodong Peninsula,eastern China. Acta Geologica Sinica,88(Suppl.2):1089 -1091

    Foster MD. 1960. Interpretation of the Composition of Trioctahedral Micas. New York:US Government Printing Office,11 -39

    Gibert F,Pascal ML and Pichavant M. 1998. Gold solubility and speciation in hydrothermal solutions:Experimental study of the stability of hydrosulphide complex of gold (AuHS)at 350 to 450℃and 500bars. Geochimica et Cosmochimica Acta,62(17):2931-2947

    Goss SC,Wilde SA,Wu FY and Yang JH. 2010. The age,isotopic signature and significance of the youngest Mesozoic granitoids in the Jiaodong Terrane,Shandong Province,North China Craton. Lithos,120(3 -4):309 -326

    Guan K,Luo ZK,Miao LC and Huang JZ. 1998. SHRIMP in zircon chronology for Guojialing suite granite in Zhaoye,Jiaodong district.Scientia Geologica Sinica,33(3):318 - 328 (in Chinese with English abstract)

    Guo P,Santosh M and Li SR. 2013. Geodynamics of gold metallogeny in the Shandong Province,NE China:An integrated geological,geophysical and geochemical perspective. Gondwana Research,24(3 -4):1172 -1202

    Hart SR. 1964. The petrology and isotopic-mineral age relations of a contact zone in the Front Range,Colorado. The Journal of Geology,72(5):493 -525

    Helmy HM,Ahmed AF,El Mahallawi MM and Ail SM. 2004. Pressure,temperature and oxygen fugacity conditions of calc-alkaline granitoids,Eastern Desert of Egypt,and tectonic implications.Journal of African Earth Sciences,38(3):255 -268

    Henry DJ,Guidotti CV and Thomson JA. 2005. The Ti-saturation surface for low-to-medium pressure metapelitic biotites:Implications for geothermometry and Ti-substitution mechanisms. American Mineralogist,90(2 -3):316 -328

    Holland T and Blundy J. 1994. Non-ideal interactions in calcic amphiboles and their bearing on amphibole-plagioclase thermometry.Contributions to Mineralogy and Petrology,116(4):433 -447

    Hou ML,Jiang YH,Jiang SY,Ling HF and Zhao KD. 2007. Contrasting origins of Late Mesozoic adakitic granitoids from the northwestern Jiaodong Peninsula,East China:Implications for crustal thickening to delamination. Geological Magazine,144(4):619 -631

    Hu FF,F(xiàn)an HR,Jiang XH,Li XC,Yang KF and Mernagh T. 2013.Fluid inclusions at different depths in the Sanshandao gold deposit,Jiaodong Peninsula,China. Geofluids,13(4):528 -541

    Huang WL and Wyllie PJ. 1975. Melting reactions in the system to 35 kilobars,dry and with excess water. The Journal of Geology,83(6):737 -748

    Jiang N,Chen JZ,Guo JH and Chang GH. 2012. In situ zircon U-Pb,oxygen and hafnium isotopic compositions of Jurassic granites from the North China craton:Evidence for Triassic subduction of continental crust and subsequent metamorphism-related18O depletion. Lithos,142 -143:84 -94

    Leake BE,Woolley AR,Arps CES et al. 1997. Nomenclature of amphiboles:Report of the Subcommittee on Amphiboles of the International Mineralogical Association Commission on New Minerals and Mineral Names. The Canadian Mineralogist,35:219 -246

    Li GM,Cao MJ,Qin KZ,Evans JN,McInnes BIA and Liu YS. 2014.Thermal-tectonic history of the Baogutu porphyry Cu deposit,West Junggar as constrained from zircon U-Pb,biotite Ar/Ar and zircon/apatite (U-Th)/He dating. Journal of Asian Earth Sciences,79:741 -758

    Li JW,Vasconcelos PM,Zhang J,Zhou MF,Zhang XJ and Yang FH.2003.40Ar/39Ar constraints on a temporal link between gold mineralization,magmatism,and continental margin transtension in the Jiaodong gold province,eastern China. The Journal of Geology,111(6):741 -751

    Li QL,Chen FK,Yang JH and Fan HR. 2008. Single grain pyrite Rb-Sr dating of the Linglong gold deposit,eastern China. Ore Geology Reviews,34(3):263 -270

    Li SR and Santosh M. 2014. Metallogeny and craton destruction:Records from the North China Craton. Ore Geology Reviews,56:376 -414 Li XC,F(xiàn)an HR,Santosh M,Hu FF,Yang KF and Lan TG. 2013.Hydrothermal alteration associated with Mesozoic granite-hosted gold mineralization at the Sanshandao deposit,Jiaodong Gold Province,China. Ore Geology Reviews,53:403 -421

    Lin W,Monié P,F(xiàn)aure M,Sch?rer U,Shi YH,Breton NL and Wang QC. 2011. Cooling paths of the NE China crust during the Mesozoic extensional tectonics:Example from the south-Liaodong peninsula metamorphic core complex. Journal of Asian Earth Sciences,42(5):1048 -1065

    Lin WW and Yin XL. 1998. The forming physicochemical conditions of Linglong granitic complex and its geological significance. Acta Geoscientia Sinica,19(1):40 - 49 (in Chinese with English abstract)

    Liu X,F(xiàn)an HR,Evans NJ,Batt GE,McInnes BIA,Yang KF and Qin KZ. 2014. Cooling and exhumation of the Mid-Jurassic porphyry copper systems in Dexing City,SE China:Insights from geo-and thermochronology. Mineralium Deposita,49(7):809 -819

    Liu ZJ,Wang JP,Zheng DW,Liu JJ,Liu J and Fu C. 2010.Exploration prospect and post-ore denudation in the northwestern Jiaodong Gold Province,China:Evidence from apatite fission track thermochronology. Acta Petrologica Sinica,26(12):3597 -3611(in Chinese with English abstract)

    Lu LN,F(xiàn)an HR,Hu FF,Yang KF and Lan TG. 2011. Emplacement depth of the Guojialing granodiorites from the northwestern Jiaodong Peninsula, eastern China: Evidences from hornblende thermobarometry and fluid inclusions. Acta Petrologica Sinica,27(5):1521 -1532 (in Chinese with English abstract)

    Luth WC,Jahns RH and Tuttle OF. 1964. The granite system at pressures of 4 to 10 kilobars. Journal of Geophysical Research,69(4):759 -773

    Lü GX. 1997. Study and estimate of depths of the formation of the Linglong and Jiaojia gold deposits,Shandong. Science China(Series D),27(4):337 -342 (in Chinese)

    Lü GX,Guo T,Shu B,Shen YK,Liu DJ and Zhou GF. 2006.Geological characteristics of rock-controlling and ore-controlling structures in the Jiaodong gold ore concentration area. Acta Geoscientica Sinica,27(5):471 -478 (in Chinese with English abstract)

    Ma CQ,Yang KG,Tang ZH,Long Y,Ehlers C and Lindroos A. 1995.Formation and differential rock uplift-exhumation of high-pressure metamorphic terrane in Dabie mountains,central China:Evidence from igneous rocks. Earth Science,20(5):515 -520 (in Chinese with English abstract)

    Nachit H,Ibhi A,Abia EH and Ohoud MB. 2005. Discrimination between primary magmatic biotites,reequilibrated biotites and neoformed biotites. Comptes Rendus Geoscience,337(16):1415 -1420

    Putirka KD. 2008. Thermometers and barometers for volcanic systems.Reviews in Mineralogy and Geochemistry,69(1):61 -120

    Qiu YM,Groves DI,McNaughton NJ,Wang LG and Zhou TH. 2002.Nature,age,and tectonic setting of granitoid-hosted,orogenic gold deposits of the Jiaodong Peninsula,eastern North China craton,China. Mineralium Deposita,37(3 -4):283 -305

    Qu XM and Wang HN. 1997. Dynamic study on the crustal-mantle magma mixing and emplacement mechanism of Guojialing granite.Scientia Geologica Sinica,32(4):445 - 454 (in Chinese with English abstract)

    Reiners PW and Brandon MT. 2006. Using thermochronology to understand orogenic erosion. Annual Review Earth Planetary Sciences,34:419 -466

    Schmit MW. 1992. Amphibole composition in tonalite as a function of pressure:An experimental calibration of the Al-in-hornblende barometer. Contribution to Mineralogy and Petrology,110(2 -3):304 -310

    Schmidt MW and Thompson AB. 1996. Epidote in calc-alkaline magmas:An experimental study of stability,phase relationships,and the role of epidote in magmatic evolution. American Mineralogist,81:462 -474

    Sial AN,Toselli AJ,Saavedra J,Parada MA and Ferreira VP. 1999.Emplacement,petrological and magnetic susceptibility characteristics of diverse magmatic epidote-bearing granitoid rocks in Brazil,Argentina and Chile. Lithos,46(3):367 -392

    Song MC,Xu JX and Wang PC. 2009. Tectonic Framework and Tectonic Evolution of the Shandong Province. Beijing:Geological Publishing House,1 -272 (in Chinese)

    Song MC,Li SZ,Santosh M et al. 2015. Types,characteristics and metallogenesis of gold deposits in the Jiaodong Peninsula,Eastern North China Craton. Ore Geology Reviews,65:612 -625

    Sun FY,Shi ZL and Feng BZ. 1995. Gold Ore Geology,Lithogenesis and Mantle-Derived C-H-O Fluids in Jiaodong Penisula,Eastern China. Changchun:Jilin People’s Press,32 -56 (in Chinese)

    Wang LG,Qiu YM,McNaughton NJ,Groves DI,Luo ZK,Huang JZ,Miao LC and Liu YK. 1998. Constraints on crustal evolution and gold metallogeny in the northwestern Jiaodong Peninsula,China,from SHRIMP U-Pb zircon studies of granitoids. Ore Geology Reviews,13(1 -5):275 -291

    Wang ZL,Yang LQ,Deng J et al. 2014. Gold-hosting high Ba-Sr granitoids in the Xincheng gold deposit,Jiaodong Peninsula,East China:Petrogenesis and tectonic setting. Journal of Asian Earth Sciences,95:274 -299

    Watson EB and Harrison TM. 1983. Zircon saturation revisited:Temperature and composition effects in a variety of crustal magma types. Earth and Planetary Science Letters,64(2):295 -304

    Wen BJ,F(xiàn)an HR,Santosh M,Hu FF,Pirajno F and Yang KF. 2014.Genesis of two different types of gold mineralization in the Linglong gold field,China:Constrains from geology,fluid inclusions and stable isotope. Ore Geology Reviews,65(3):643 -658

    Whitney DL and Evans BW. 2010. Abbreviations for names of rockforming minerals. American Mineralogist,95(1):185 -187

    Williams-Jones AE,Bowell RJ and Migdisov AA. 2009. Gold in solution. Elements,5(5):281 -287

    Wones DR and Eugster HP. 1965. Stability of biotite:Experiment,theory,and application. American Mineralogist,50 (9):1228-1272

    Wones DR. 1989. Significance of the assemblage titanite + magnetite +quartz in granitic rocks. American Mineralogist,74(7 -8):744-749

    Xu JF and Shen BY. 1989. On the granitoid related to gold mineralization in Jiaodong block. Geology of Shandong,5(2):1 - 125 (in Chinese with English abstract)

    Yang JH and Zhou XH. 2000. The Rb-Sr isochron of ore and pyrite subsamples from Linglong gold deposit,Jiaodong Peninsula,eastern China and their geological significance. Chinese Science Bulletin,45(24):2272 -2277

    Yang JH and Zhou XH. 2001. Rb-Sr,Sm-Nd,and Pb isotope systematics of pyrite:Implications for the age and genesis of lode gold deposits. Geology,29(8):711 -714

    Yang JH,Wu FY and Wilde SA. 2003. A review of the geodynamic setting of large-scale Late Mesozoic gold mineralization in the North China Craton:An association with lithospheric thinning. Ore Geology Reviews,23(3 -4):125 -152

    Yang JH,Zhu MF,Liu W and Zhai MG. 2003. Geochemistry and petrogenesis of Guojialing granodiorites from the northwestern Jiaodong Peninsula,eastern China. Acta Petrologica Sinica,19(4):692 -700 (in Chinese with English abstract)

    Yang JH,Wu FY,Shao JA,Wide SA,Xie LW and Liu XM. 2006.Constraints on the timing of uplift of the Yanshan Fold and Thrust Belt,North China. Earth and Planetary Science Letters,246(3):336 -352

    Yang KF,F(xiàn)an HR,Santosh M,Hu FF,Wilde SA,Lan TG,Lu LN and Liu YS. 2012. Reactivation of the Archean lower crust:Implications for zircon geochronology, elemental and Sr-Nd-Hf isotopic geochemistry of Late Mesozoic granitoids from northwestern Jiaodong Terrane,the North China Craton. Lithos,146 -147:112 -127 Yang LQ,Deng J,Goldfarb RJ,Zhang J,Gao BF and Wang ZL. 2014a.40Ar/39Ar geochronological constraints on the formation of the Dayingezhuang gold deposit:New implications for timing and duration of hydrothermal activity in the Jiaodong gold province,China. Gondwana Research,25(4):1469 -1483

    Yang LQ,Deng J,Wang ZL,Zhang L,Guo LN,Song MC and Zheng XL. 2014. Mesozoic gold metallogenic system of the Jiaodong gold province,eastern China. Acta Petrologica Sinica,30(9):2447 -2467 (in Chinese with English abstract)

    Yang QY,Santosh M,Shen JF and Li SR. 2014b. Juvenile vs. recycled crust in NE China:Zircon U-Pb geochronology,Hf isotopes and an integrated model for Mesozoic gold mineralization in the Jiaodong Peninsula. Gondwana Research,25(4):1445 -1468

    Yang WT,Yang JH,Wang XF and Du YS. 2014c. Uplift-denudation history of the Qinling orogen:Constrained from the detrital-zircon UPb geochronology. Journal of Asian Earth Sciences,89:54 -65

    Yu XF,Li HK and Shan W. 2012. Study on coupling between Yanshannian tectonic thermal events and gold mineralization in Jiaodong ore concentrating area in Shandong Province. Acta Geologica Sinica,86(12):1946 -1956 (in Chinese with English abstract)

    Zhai MG,Cong BL,Guo JH,Liu WJ,Li YG and Wang QC. 2000. Sm-Nd geochronology and petrography of garnet pyroxene granulites in the northern Sulu region of China and their geotectonic implication.Lithos,52(1):23 -33

    Zhai MG,Yang JH,F(xiàn)an HR,Miao LC and Li YG. 2002. A large-scale cluster of gold deposits and metallogenesis in the eastern North China Craton. International Geology Review,44(5):458 -476

    Zhang J,Zhao ZF,Zheng YF and Deng MN. 2010. Postcollisional magmatism:Geochemical constraints on the petrogenesis of Mesozoic granitoids in the Sulu orogen,China. Lithos,119(3):512 -536

    Zhang SH,Zhao Y and Song B. 2006. Hornblende thermobarometry of the Carboniferous granitoids from the Inner Mongolia Paleo-uplift:Implications for the tectonic evolution of the northern margin of North China block. Mineralogy and Petrology,87(1 -2):123 -141

    Zhao GC,Wilde SA,Cawood PA and Lu LZ. 1999. Thermal evolution of two types of mafic granulites from the North China Craton:Implications for both mantle plume and collisional tectonics.Geological Magazine,136(3):223 -240

    Zhou XH,Yang JH and Zhang LC. 2003. Metallogenesis of superlarge gold deposits in Jiaodong region and deep processes of subcontinental lithosphere beneath North China Craton in Mesozoic. Science in China (Series D),46(Suppl.1):14 -25

    Симонова Л И. 1979. Темноцветные минералы-индикаторы глубины кристаллизации гранитоидов (на примере интрузивов Кураминской зоны,Средняя Азия). Геохимия,(9):1307-1322

    ΠeрчукЛЛ HykM. 1973. Tepmoдинамигeckий peжимг Убиного петрогенезиса. Hayka

    附中文參考文獻(xiàn)

    陳光遠(yuǎn),孫岱生,周珣若,邵偉,宮潤(rùn)潭,邵岳. 1993. 膠東郭家?guī)X花崗閃長(zhǎng)巖成因礦物學(xué)與金礦化. 武漢:中國(guó)地質(zhì)大學(xué)出版社,8 -219

    陳衍景,Pirajno F,賴勇,李超. 2004. 膠東礦集區(qū)大規(guī)模成礦時(shí)間和構(gòu)造環(huán)境. 巖石學(xué)報(bào),20(4):907 -922

    鄧軍,徐守禮,方云,周顯強(qiáng),萬(wàn)麗. 1996. 膠東西北部構(gòu)造體系及金成礦動(dòng)力學(xué). 北京:地質(zhì)出版社,1 -98

    鄧軍,翟裕生,楊立強(qiáng),肖榮閣,孫忠實(shí). 1999. 構(gòu)造演化與成礦系統(tǒng)動(dòng)力學(xué)——以膠東金礦集中區(qū)為例. 地學(xué)前緣,6(2):315-323

    范宏瑞,胡芳芳,楊進(jìn)輝,沈昆,翟明國(guó). 2005. 膠東中生代構(gòu)造體制轉(zhuǎn)折過(guò)程中流體演化和金的大規(guī)模成礦. 巖石學(xué)報(bào),21(5):1317 -1328

    關(guān)康,羅鎮(zhèn)寬,苗來(lái)成,黃佳展. 1998. 膠東招掖郭家?guī)X型花崗巖鋯石SHRIMP 年代學(xué)研究. 地質(zhì)科學(xué),33(3):318 -328

    林文蔚,殷秀蘭. 1998. 玲瓏花崗質(zhì)雜巖體形成的物理化學(xué)條件及其地質(zhì)意義. 地球?qū)W報(bào),19(1):40 -49

    柳振江,王建平,鄭德文,劉家軍,劉俊,付超. 2010. 膠東西北部金礦剝蝕程度及找礦潛力和方向——來(lái)自磷灰石裂變徑跡熱年代學(xué)的證據(jù). 巖石學(xué)報(bào),26(12):3597 -3611

    陸麗娜,范宏瑞,胡芳芳,楊奎鋒,藍(lán)廷廣. 2011. 膠西北郭家?guī)X花崗閃長(zhǎng)巖侵位深度:來(lái)自角閃石溫壓計(jì)和流體包裹體的證據(jù).巖石學(xué)報(bào),27(5):1521 -1532

    呂古賢. 1997. 山東玲瓏金礦田和焦家金礦田成礦深度的測(cè)算與研究方法. 中國(guó)科學(xué)(D 輯),27(4):337 -342

    呂古賢,郭濤,舒斌,申玉科,劉杜鵑,周?chē)?guó)發(fā). 2006. 膠東金礦集中區(qū)構(gòu)造控巖控礦地質(zhì)特征研究. 地球?qū)W報(bào),27(5):471 -478

    馬昌前,楊坤光,唐仲華,龍昱,Ehlers C,Lindroos A. 1995. 華中大別山高壓變質(zhì)地體的形成和差異巖石隆升-剝蝕:來(lái)自火成巖的證據(jù). 地球科學(xué),20(5):515 -520

    曲曉明,王鶴年. 1997. 郭家?guī)X巖體殼慢巖漿混合作用與侵位機(jī)制的動(dòng)力學(xué)研究. 地質(zhì)科學(xué),32(4):445 -454

    宋明春,徐軍祥,王沛成. 2009. 山東省大地構(gòu)造格局和地質(zhì)構(gòu)造演化. 北京:地質(zhì)出版社,1 -272

    孫豐月,石準(zhǔn)立,馮本智. 1995. 膠東金礦地質(zhì)及幔源C-H-O 流體分異成巖成礦. 長(zhǎng)春:吉林人民出版社,32 -56

    徐金方,沈步云. 1989. 膠北地塊與金礦有關(guān)的花崗巖類的研究. 山東地質(zhì),5(2):1 -125

    楊進(jìn)輝,朱美妃,劉偉,翟明國(guó). 2003. 膠東地區(qū)郭家?guī)X花崗閃長(zhǎng)巖的地球化學(xué)特征及成因. 巖石學(xué)報(bào),19(4):692 -700

    楊立強(qiáng),鄧軍,王中亮,張良,郭林楠,宋明春,鄭小禮. 2014. 膠東中生代金成礦系統(tǒng). 巖石學(xué)報(bào),30(9):2447 -2467

    于學(xué)峰,李洪奎,單偉. 2012. 山東膠東礦集區(qū)燕山期構(gòu)造熱事件與金礦成礦耦合探討. 地質(zhì)學(xué)報(bào),86(12):1946 -1956

    周新華,楊進(jìn)輝,張連昌. 2002. 膠東大型金礦的形成與中生代華北大陸巖石圈深部過(guò)程. 中國(guó)科學(xué)(D 輯),32(增刊):11 -20

    猜你喜歡
    隆升角閃石斜長(zhǎng)石
    角閃石礦物化學(xué)研究進(jìn)展
    角閃石高溫高壓實(shí)驗(yàn)研究進(jìn)展及其地球物理意義
    桂北新寨角閃花崗巖角閃石礦物化學(xué)特征及其指示意義
    內(nèi)蒙古獲各琦地區(qū)海西期石英閃長(zhǎng)巖中角閃石成因礦物學(xué)研究
    科科斯脊玄武巖斜長(zhǎng)石礦物化學(xué)及地質(zhì)意義
    南北構(gòu)造帶北部“古脊梁”演化過(guò)程探討
    西南印度洋中脊63.9°E斜長(zhǎng)石超斑狀玄武巖對(duì)超慢速擴(kuò)張洋脊巖漿過(guò)程的指示*
    不同濃度水楊酸對(duì)斜長(zhǎng)石溶解能力的研究
    南迦巴瓦峰第四紀(jì)隆升期次劃分的熱年代學(xué)證據(jù)
    測(cè)定、花崗巖類巖石中長(zhǎng)石成分的分析方法
    地球(2016年10期)2016-04-14 21:07:54
    国产高清videossex| av在线天堂中文字幕| 亚洲欧美日韩东京热| 91在线精品国自产拍蜜月 | 日本精品一区二区三区蜜桃| 久久久久久人人人人人| 国产激情偷乱视频一区二区| 两人在一起打扑克的视频| 首页视频小说图片口味搜索| 亚洲成人精品中文字幕电影| 欧美高清成人免费视频www| 啦啦啦观看免费观看视频高清| 久久久久久国产a免费观看| av在线蜜桃| 日韩欧美在线乱码| 高潮久久久久久久久久久不卡| 国产在视频线在精品| 三级毛片av免费| 国产精品久久视频播放| 久9热在线精品视频| 老司机福利观看| 精品国内亚洲2022精品成人| 精品熟女少妇八av免费久了| 亚洲成av人片免费观看| 啪啪无遮挡十八禁网站| 欧美午夜高清在线| 国产免费av片在线观看野外av| 亚洲在线观看片| 午夜福利在线观看吧| 一区福利在线观看| 一区福利在线观看| 久久99热这里只有精品18| 人人妻,人人澡人人爽秒播| 亚洲国产精品合色在线| 人人妻人人看人人澡| 99久久精品国产亚洲精品| 亚洲精品影视一区二区三区av| 一夜夜www| 日日摸夜夜添夜夜添小说| av女优亚洲男人天堂| 国产精品免费一区二区三区在线| 日韩成人在线观看一区二区三区| 国产精品影院久久| 男人舔奶头视频| 亚洲av成人精品一区久久| 精品无人区乱码1区二区| 丝袜美腿在线中文| 亚洲av日韩精品久久久久久密| 国产精品香港三级国产av潘金莲| 欧美三级亚洲精品| 国内精品一区二区在线观看| 国产高清videossex| 欧美乱妇无乱码| 91在线精品国自产拍蜜月 | 国产伦人伦偷精品视频| 丰满人妻一区二区三区视频av | 他把我摸到了高潮在线观看| 亚洲精品亚洲一区二区| 免费无遮挡裸体视频| 欧美日韩瑟瑟在线播放| 特大巨黑吊av在线直播| 午夜福利在线观看吧| 噜噜噜噜噜久久久久久91| 国产精品久久久人人做人人爽| 又粗又爽又猛毛片免费看| 波野结衣二区三区在线 | 久久香蕉精品热| 一本综合久久免费| 久久伊人香网站| 色精品久久人妻99蜜桃| 黄片小视频在线播放| 久久久成人免费电影| 搡女人真爽免费视频火全软件 | 一本久久中文字幕| 国产精品自产拍在线观看55亚洲| 亚洲av第一区精品v没综合| 日本 欧美在线| 国产精品精品国产色婷婷| 日韩欧美三级三区| www.999成人在线观看| 欧美性感艳星| 免费大片18禁| 一区福利在线观看| 又黄又爽又免费观看的视频| bbb黄色大片| bbb黄色大片| 国产97色在线日韩免费| 亚洲av不卡在线观看| 亚洲人成电影免费在线| 久久久色成人| 男女之事视频高清在线观看| 色吧在线观看| www.999成人在线观看| 久久久久久久午夜电影| a级毛片a级免费在线| 亚洲成人久久性| 搞女人的毛片| 中文字幕久久专区| 在线播放无遮挡| 99久久精品一区二区三区| 成人av在线播放网站| 欧美激情在线99| 日韩欧美免费精品| 成人特级av手机在线观看| 男女床上黄色一级片免费看| 日韩欧美一区二区三区在线观看| eeuss影院久久| 精品一区二区三区人妻视频| 日韩欧美 国产精品| 每晚都被弄得嗷嗷叫到高潮| 亚洲国产欧美网| 亚洲av成人精品一区久久| 日韩欧美在线二视频| 亚洲精品456在线播放app | 欧美激情久久久久久爽电影| 国产亚洲精品av在线| 免费在线观看日本一区| 国产精品 欧美亚洲| 色视频www国产| 久99久视频精品免费| 男插女下体视频免费在线播放| 欧美激情在线99| 在线观看美女被高潮喷水网站 | 久久久久性生活片| 亚洲成av人片免费观看| 色精品久久人妻99蜜桃| 狠狠狠狠99中文字幕| a级毛片a级免费在线| 久久精品影院6| 五月玫瑰六月丁香| 18禁国产床啪视频网站| 成人三级黄色视频| 搡女人真爽免费视频火全软件 | aaaaa片日本免费| 两个人视频免费观看高清| 亚洲av熟女| 日本一二三区视频观看| 99久久精品国产亚洲精品| 又爽又黄无遮挡网站| 国产v大片淫在线免费观看| 日韩欧美在线二视频| 久久人妻av系列| 丁香六月欧美| 国产精品电影一区二区三区| 日本黄色片子视频| 国产亚洲精品久久久com| 亚洲内射少妇av| www.色视频.com| 岛国视频午夜一区免费看| 久久久精品欧美日韩精品| 最后的刺客免费高清国语| 男插女下体视频免费在线播放| 午夜老司机福利剧场| 亚洲av免费高清在线观看| 天天躁日日操中文字幕| a在线观看视频网站| 18禁黄网站禁片午夜丰满| 亚洲男人的天堂狠狠| 国产精品电影一区二区三区| 看黄色毛片网站| 老司机午夜福利在线观看视频| 首页视频小说图片口味搜索| 老司机午夜福利在线观看视频| 制服人妻中文乱码| 法律面前人人平等表现在哪些方面| 国产v大片淫在线免费观看| 99热这里只有精品一区| 婷婷丁香在线五月| 欧洲精品卡2卡3卡4卡5卡区| www.熟女人妻精品国产| 午夜久久久久精精品| 欧美最黄视频在线播放免费| 免费一级毛片在线播放高清视频| 波多野结衣巨乳人妻| 午夜免费男女啪啪视频观看 | 叶爱在线成人免费视频播放| 嫩草影院精品99| 听说在线观看完整版免费高清| or卡值多少钱| 18禁在线播放成人免费| 观看美女的网站| 叶爱在线成人免费视频播放| 观看美女的网站| 女人十人毛片免费观看3o分钟| 免费人成视频x8x8入口观看| 母亲3免费完整高清在线观看| 国产精品久久视频播放| 夜夜爽天天搞| or卡值多少钱| 18禁裸乳无遮挡免费网站照片| 亚洲欧美激情综合另类| 中出人妻视频一区二区| 免费看日本二区| 少妇的丰满在线观看| 中出人妻视频一区二区| 成人一区二区视频在线观看| 欧美日韩瑟瑟在线播放| 国产一区二区三区视频了| 又粗又爽又猛毛片免费看| 搡老岳熟女国产| 久久香蕉国产精品| 最近最新免费中文字幕在线| 午夜福利视频1000在线观看| 一区二区三区免费毛片| 免费观看的影片在线观看| 精品久久久久久成人av| 日韩精品中文字幕看吧| 欧美bdsm另类| 亚洲熟妇熟女久久| 一区二区三区免费毛片| 久久草成人影院| 日日摸夜夜添夜夜添小说| 亚洲av五月六月丁香网| 身体一侧抽搐| 少妇人妻一区二区三区视频| 一进一出好大好爽视频| 亚洲av成人av| 俄罗斯特黄特色一大片| 欧美乱妇无乱码| svipshipincom国产片| 18+在线观看网站| 黄色成人免费大全| 中文资源天堂在线| 国产极品精品免费视频能看的| 亚洲精品日韩av片在线观看 | 中文字幕人妻丝袜一区二区| 国产一区二区三区在线臀色熟女| 最近最新中文字幕大全免费视频| 亚洲五月婷婷丁香| 美女高潮的动态| 真实男女啪啪啪动态图| 人人妻人人看人人澡| 免费av不卡在线播放| 久久久久亚洲av毛片大全| 精品国产美女av久久久久小说| 国产精品,欧美在线| 亚洲av熟女| 一本一本综合久久| 日韩有码中文字幕| 欧美中文日本在线观看视频| 级片在线观看| 少妇的逼水好多| 99久久九九国产精品国产免费| 在线观看免费视频日本深夜| 国产精品综合久久久久久久免费| 免费看a级黄色片| 丰满的人妻完整版| 久久精品综合一区二区三区| 观看美女的网站| 亚洲欧美精品综合久久99| 国内精品久久久久精免费| 成年版毛片免费区| 国产精品香港三级国产av潘金莲| 亚洲av一区综合| 免费看十八禁软件| 国产成人欧美在线观看| 亚洲av中文字字幕乱码综合| 又紧又爽又黄一区二区| 亚洲精品成人久久久久久| 一进一出好大好爽视频| 国产精品爽爽va在线观看网站| 真人一进一出gif抽搐免费| 国产v大片淫在线免费观看| 欧美+日韩+精品| 神马国产精品三级电影在线观看| 九色成人免费人妻av| 国模一区二区三区四区视频| 老司机午夜十八禁免费视频| 亚洲人成网站高清观看| 国产精品亚洲一级av第二区| 亚洲美女视频黄频| 欧美zozozo另类| 可以在线观看的亚洲视频| 人妻久久中文字幕网| 精品久久久久久久久久免费视频| 久久久久免费精品人妻一区二区| 亚洲电影在线观看av| 在线观看一区二区三区| 亚洲av电影在线进入| 久久九九热精品免费| 亚洲av不卡在线观看| 日韩欧美在线乱码| 国产一级毛片七仙女欲春2| 久久欧美精品欧美久久欧美| 久久精品国产自在天天线| 久久99热这里只有精品18| 久久亚洲真实| 黄色丝袜av网址大全| 熟妇人妻久久中文字幕3abv| 久久精品91蜜桃| 最好的美女福利视频网| 国产成人福利小说| 成人性生交大片免费视频hd| 窝窝影院91人妻| 麻豆久久精品国产亚洲av| 久久这里只有精品中国| 在线观看美女被高潮喷水网站 | 国产不卡一卡二| 一个人免费在线观看的高清视频| 亚洲欧美日韩高清专用| 国内久久婷婷六月综合欲色啪| 老司机深夜福利视频在线观看| 欧美+亚洲+日韩+国产| 色吧在线观看| 欧美乱妇无乱码| 观看免费一级毛片| 国内精品一区二区在线观看| 波多野结衣高清无吗| 欧美黑人欧美精品刺激| 色播亚洲综合网| xxxwww97欧美| 日日摸夜夜添夜夜添小说| 亚洲成人免费电影在线观看| 国产一区在线观看成人免费| 国产精品日韩av在线免费观看| 嫁个100分男人电影在线观看| 欧美极品一区二区三区四区| 在线国产一区二区在线| 麻豆久久精品国产亚洲av| 男女做爰动态图高潮gif福利片| 亚洲av一区综合| 国产中年淑女户外野战色| 亚洲欧美日韩高清在线视频| 又紧又爽又黄一区二区| 精华霜和精华液先用哪个| 亚洲性夜色夜夜综合| 精品欧美国产一区二区三| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | а√天堂www在线а√下载| 成年女人毛片免费观看观看9| 亚洲在线观看片| 午夜两性在线视频| 精品国产美女av久久久久小说| 两个人视频免费观看高清| 精品日产1卡2卡| 色综合亚洲欧美另类图片| 99国产精品一区二区蜜桃av| 夜夜夜夜夜久久久久| 99久久99久久久精品蜜桃| 婷婷精品国产亚洲av在线| 久久久久久久午夜电影| 日韩有码中文字幕| 啦啦啦韩国在线观看视频| 亚洲精品一卡2卡三卡4卡5卡| 97人妻精品一区二区三区麻豆| 国产一区二区三区在线臀色熟女| 黄色片一级片一级黄色片| 制服人妻中文乱码| 国产视频一区二区在线看| 久久精品国产清高在天天线| 99久久无色码亚洲精品果冻| 亚洲av不卡在线观看| 三级男女做爰猛烈吃奶摸视频| 国产欧美日韩一区二区三| 日韩av在线大香蕉| 色综合欧美亚洲国产小说| 亚洲色图av天堂| 亚洲激情在线av| 毛片女人毛片| 午夜精品一区二区三区免费看| 中文亚洲av片在线观看爽| 国内揄拍国产精品人妻在线| 亚洲国产精品成人综合色| 亚洲久久久久久中文字幕| av天堂中文字幕网| 免费电影在线观看免费观看| 性色avwww在线观看| 亚洲在线自拍视频| 1000部很黄的大片| 俄罗斯特黄特色一大片| 黄色丝袜av网址大全| 成年版毛片免费区| 中文字幕人成人乱码亚洲影| 日韩av在线大香蕉| 黄色丝袜av网址大全| 一区二区三区激情视频| 88av欧美| 欧美日韩一级在线毛片| 免费av观看视频| 国产av在哪里看| 在线观看美女被高潮喷水网站 | 国产主播在线观看一区二区| 母亲3免费完整高清在线观看| 欧美另类亚洲清纯唯美| 成人国产综合亚洲| 草草在线视频免费看| 日本一二三区视频观看| 成人欧美大片| 又粗又爽又猛毛片免费看| 51午夜福利影视在线观看| 亚洲不卡免费看| 亚洲中文字幕日韩| 国产精品日韩av在线免费观看| 久久午夜亚洲精品久久| 精品一区二区三区av网在线观看| 在线观看av片永久免费下载| 亚洲内射少妇av| 桃色一区二区三区在线观看| 老司机在亚洲福利影院| 国产欧美日韩一区二区三| 麻豆成人午夜福利视频| 怎么达到女性高潮| 男人舔奶头视频| 丝袜美腿在线中文| 精品久久久久久久久久免费视频| 国产三级中文精品| 亚洲av熟女| 亚洲第一欧美日韩一区二区三区| 成人永久免费在线观看视频| 午夜影院日韩av| 免费看美女性在线毛片视频| 色综合欧美亚洲国产小说| 成人特级黄色片久久久久久久| 色综合站精品国产| 精品午夜福利视频在线观看一区| av福利片在线观看| 哪里可以看免费的av片| 免费看日本二区| 色老头精品视频在线观看| 在线免费观看的www视频| 成年人黄色毛片网站| 国产精品99久久久久久久久| 欧美黄色片欧美黄色片| 欧美中文日本在线观看视频| 精品不卡国产一区二区三区| 久久香蕉国产精品| 男女视频在线观看网站免费| 国产在线精品亚洲第一网站| 午夜福利在线在线| 一个人观看的视频www高清免费观看| 啦啦啦观看免费观看视频高清| 变态另类成人亚洲欧美熟女| 亚洲熟妇中文字幕五十中出| 夜夜夜夜夜久久久久| 一个人看的www免费观看视频| 国产成年人精品一区二区| 热99re8久久精品国产| 成人国产综合亚洲| 哪里可以看免费的av片| www国产在线视频色| 在线观看免费视频日本深夜| 欧美一区二区国产精品久久精品| 精品日产1卡2卡| 身体一侧抽搐| 亚洲七黄色美女视频| 观看免费一级毛片| 午夜福利18| 亚洲av免费高清在线观看| 国产亚洲精品综合一区在线观看| 3wmmmm亚洲av在线观看| 男插女下体视频免费在线播放| 免费一级毛片在线播放高清视频| 在线国产一区二区在线| 91久久精品电影网| 亚洲国产欧美人成| 欧美国产日韩亚洲一区| 亚洲 国产 在线| 日韩av在线大香蕉| 99热6这里只有精品| 成人特级黄色片久久久久久久| 久久精品91无色码中文字幕| 国内精品美女久久久久久| 一级黄片播放器| 97人妻精品一区二区三区麻豆| 亚洲av电影不卡..在线观看| 亚洲av美国av| 91在线观看av| 免费搜索国产男女视频| 网址你懂的国产日韩在线| 一区二区三区免费毛片| 三级国产精品欧美在线观看| 老司机午夜十八禁免费视频| 国产主播在线观看一区二区| 一本精品99久久精品77| 国产成人福利小说| 香蕉久久夜色| 18美女黄网站色大片免费观看| 色哟哟哟哟哟哟| 麻豆成人午夜福利视频| 欧美av亚洲av综合av国产av| 制服人妻中文乱码| 噜噜噜噜噜久久久久久91| 激情在线观看视频在线高清| 真人一进一出gif抽搐免费| 男女午夜视频在线观看| 精品人妻一区二区三区麻豆 | 99在线视频只有这里精品首页| 欧美日本亚洲视频在线播放| 国产高清videossex| 波多野结衣高清作品| 中文资源天堂在线| 嫩草影院入口| 日本免费一区二区三区高清不卡| 成人亚洲精品av一区二区| 国产欧美日韩一区二区三| 麻豆一二三区av精品| 一夜夜www| 亚洲电影在线观看av| 美女 人体艺术 gogo| 午夜a级毛片| 一级毛片女人18水好多| 香蕉丝袜av| 国产亚洲精品久久久久久毛片| 午夜激情欧美在线| 啦啦啦免费观看视频1| 99热这里只有是精品50| 国产欧美日韩精品亚洲av| 免费av观看视频| 亚洲精华国产精华精| 国产爱豆传媒在线观看| 校园春色视频在线观看| 国产成人福利小说| 成人国产综合亚洲| 在线观看免费午夜福利视频| 男女下面进入的视频免费午夜| 精品久久久久久久久久免费视频| 欧美日韩乱码在线| 女同久久另类99精品国产91| 男人的好看免费观看在线视频| 久久久久亚洲av毛片大全| 国产精品 欧美亚洲| 久久久久久久午夜电影| 黄色日韩在线| 午夜激情欧美在线| 日韩欧美一区二区三区在线观看| 18+在线观看网站| 久久天躁狠狠躁夜夜2o2o| ponron亚洲| 99热这里只有精品一区| 久久久久久久午夜电影| 波野结衣二区三区在线 | 熟女电影av网| e午夜精品久久久久久久| 久久精品国产清高在天天线| 搡老熟女国产l中国老女人| 国产成人影院久久av| 人人妻人人澡欧美一区二区| 九色国产91popny在线| 叶爱在线成人免费视频播放| 精品国产亚洲在线| 老司机午夜福利在线观看视频| 国产精品一区二区三区四区久久| 亚洲欧美日韩东京热| av女优亚洲男人天堂| 一a级毛片在线观看| 中文字幕久久专区| 91久久精品国产一区二区成人 | 免费观看的影片在线观看| 久久国产乱子伦精品免费另类| 亚洲av成人精品一区久久| 男女之事视频高清在线观看| 色哟哟哟哟哟哟| 亚洲av电影不卡..在线观看| 国产精品日韩av在线免费观看| 国产黄片美女视频| 一级作爱视频免费观看| 亚洲美女视频黄频| 91在线观看av| 亚洲最大成人手机在线| 亚洲精品成人久久久久久| 窝窝影院91人妻| 天天躁日日操中文字幕| 午夜免费成人在线视频| 国产在线精品亚洲第一网站| 日日干狠狠操夜夜爽| 日韩免费av在线播放| 亚洲精品乱码久久久v下载方式 | 亚洲专区国产一区二区| 亚洲成a人片在线一区二区| 亚洲精品色激情综合| 国产免费男女视频| av女优亚洲男人天堂| 亚洲美女视频黄频| 欧美乱码精品一区二区三区| 久久久久九九精品影院| 亚洲成人中文字幕在线播放| 日日摸夜夜添夜夜添小说| 精品久久久久久久末码| 村上凉子中文字幕在线| 国内精品久久久久精免费| 色播亚洲综合网| 午夜福利在线在线| 国产成人a区在线观看| 91在线精品国自产拍蜜月 | 91麻豆av在线| 欧美日韩中文字幕国产精品一区二区三区| 欧美日韩综合久久久久久 | 天美传媒精品一区二区| 日韩欧美 国产精品| 成人18禁在线播放| 免费观看人在逋| 久久久久久九九精品二区国产| 黄色丝袜av网址大全| 欧美日韩乱码在线| 午夜福利高清视频| 制服人妻中文乱码| 岛国在线免费视频观看| 国产精品久久久久久亚洲av鲁大| 日本一本二区三区精品| 91久久精品国产一区二区成人 | 成人性生交大片免费视频hd| 丰满的人妻完整版| 国产探花极品一区二区| 女警被强在线播放| 婷婷六月久久综合丁香| 久久久久免费精品人妻一区二区| 久9热在线精品视频| 国产精品久久久久久精品电影| 国产免费男女视频| 亚洲人与动物交配视频| 久久久久久久久久黄片| 欧美不卡视频在线免费观看| 97超级碰碰碰精品色视频在线观看| 一二三四社区在线视频社区8| av天堂中文字幕网|