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

    山西靈丘縣刁泉銀銅礦流體包裹體特征及成礦流體演化*

    2014-04-13 04:09:28牛斯達(dá)李勝榮張聚全王志華
    巖石學(xué)報(bào) 2014年11期
    關(guān)鍵詞:子石矽卡巖鹽度

    牛斯達(dá) 李勝榮** 張聚全,2 王志華

    NIU SiDa1,LI ShengRong1**,ZHANG JuQuan1,2 and WANG ZhiHua3

    1. 地質(zhì)過(guò)程與礦產(chǎn)資源國(guó)家重點(diǎn)實(shí)驗(yàn)室,中國(guó)地質(zhì)大學(xué)地球科學(xué)與資源學(xué)院,北京 100083

    2. 石家莊經(jīng)濟(jì)學(xué)院寶石與材料工藝學(xué)院,石家莊 050031

    3. 西安地質(zhì)調(diào)查中心,西安 710054

    1. State Key Laboratory of Geological Processes and Mineral Resources,School of the Earth Sciences and Resources,China University of Geosciences,Beijing 100083,China

    2. School of Gemology and Materials Technology,Shijiazhuang University of Economics,Shijiazhuang 050031,China

    3. Xi’an Center of Geological Survey,CGS,Xi’an 710054,China

    2014-02-27 收稿,2014-06-18 改回.

    刁泉銀銅礦位于山西省東北部靈丘縣,地處太行山北段,是一個(gè)以銀銅為主的矽卡巖型多金屬礦床。前人對(duì)刁泉銀銅礦開(kāi)展了許多研究,周利霞和唐耀林(1997)初步建立了其成礦模式;李兆龍和張連營(yíng)(1999)對(duì)其成因機(jī)制進(jìn)行了探索;李喜榮等(1999)對(duì)金銀賦存狀態(tài)進(jìn)行了研究;李曉剛等(2004)初步建立了找礦模式;王永明等(2006)在礦物學(xué)特征方面、李愛(ài)生和韓志安(2001)在化探方面、李惠等(2008)在構(gòu)造疊加暈方面也取得了一些成果。

    許多學(xué)者對(duì)矽卡巖型銀銅金礦床從地球化學(xué)、流體包裹體、成礦機(jī)制等方面進(jìn)行過(guò)較深入研究(Morrison et al.,1991;Johnson and Meinert,1994;Zhao et al.,1999)。流體包裹體能夠很好地指示流體成礦的全過(guò)程(范宏瑞等,1998;Shepherd et al.,1985),但目前刁泉銀銅礦床的成礦流體特征研究程度較低,在很大程度上制約了成礦模型的建立。本文在描述刁泉銀銅礦地質(zhì)特征的基礎(chǔ)上,利用流體包裹體顯微測(cè)溫和激光拉曼光譜分析技術(shù)對(duì)該礦床的流體特征和演化過(guò)程進(jìn)行了研究,進(jìn)一步揭示了礦床的形成機(jī)制。

    1 礦床地質(zhì)特征

    圖1 華北克拉通位置圖(a,據(jù)湯艷杰等,2006;羅軍燕,2009)和刁泉銀銅礦區(qū)域地質(zhì)圖(b,據(jù)周利霞和唐耀林,1997;易洪波,2012)1-第四系;2-白堊系;3-中侏羅系火山巖;4-二疊系頁(yè)巖、砂巖;5-石炭系頁(yè)巖;6-奧陶系灰?guī)r、頁(yè)巖;7-寒武系灰?guī)r、頁(yè)巖、白云巖;8-震旦系砂巖;9-太古界五臺(tái)群變質(zhì)巖系;10-燕山期花崗巖類;11-燕山期閃長(zhǎng)巖類;12-元古代閃長(zhǎng)巖類;13-太古宙花崗巖類;14-向斜軸;15-隱伏斷裂;16-斷裂;17-礦區(qū)位置Fig.1 Location of the North China Craton (a,modified after Tang et al.,2006;Luo,2009)and regional geological map of the Diaoquan Ag-Cu deposit (b,modified after Zhou and Tang,1997;Yi,2012)1-Quaternary;2-Cretaceous;3-Middle Jurassic volcanic rocks;4-Permian shale and sandstone;5-Carboniferous shale;6-Ordovician limestone and shale;7-Cambrian limestone,shale and dolomite;8-Sinian sandstone;9-metamorphic rocks of Achaean Wutai Group;10-Yanshanian granite;11-Yanshanian diorite;12-Proterozoic diorite;13-Archaean granite;14-synclinal axis;15-buried fault;16-fault;17-location of the deposit

    圖2 刁泉銀銅礦礦區(qū)地質(zhì)圖(據(jù)周利霞和唐耀林,1997;易洪波,2012)1-凝灰角礫巖夾流紋質(zhì)熔巖;2-安山質(zhì)角礫巖及凝灰?guī)r;3-凝灰角礫巖;4-中、上奧陶統(tǒng)大理巖;5-中、上寒武統(tǒng)大理巖;6-礦體及編號(hào);7-黑云母石英二長(zhǎng)巖;8-花崗斑巖;9-輝長(zhǎng)巖;10-閃長(zhǎng)玢巖脈;11-安山巖;12-細(xì)晶巖脈;13-煌斑巖脈;14-輝綠玢巖脈;15-矽卡巖;16-背斜軸;17-斷裂;18-推測(cè)巖體邊界或環(huán)形斷裂;19-巖層產(chǎn)狀Fig. 2 Geological map of the Diaoquan Ag-Cu deposit(modified after Zhou and Tang,1997;Yi,2012)1-tuff breccia with rhyolitic lava;2-andesitic breccia and tuff;3-tuff breccia;4-the Upper and Middle Ordovician marble;5-the Upper and Middle Cambrian marble;6-ore bodies and their numbers;7-biotite quartz monzonite;8-granite porphyry;9-gabbro;10-diorite porphyrite dyke;11-andesite;12-aplite dyke;13-lamprophyre dyke;14-diabase-porphyrite dike;15-skarn;16-anticline axis;17-fault;18-presumable rock boundaries or ring fault;19-occurrence of layers

    刁泉銀銅礦處于華北克拉通中北部,區(qū)域上為一中生代斷陷火山盆地(周利霞和唐耀林,1997)。該礦床位于燕山造山帶與太行山造山帶的交切部位(Deng et al.,2007;羅照華等,1999),區(qū)內(nèi)出露地層主要為五臺(tái)群中淺變質(zhì)的碎屑巖、中酸性火山巖建造,其次為震旦系砂巖,奧陶系灰?guī)r、頁(yè)巖和寒武系灰?guī)r、頁(yè)巖、白云巖等(圖1)。

    區(qū)內(nèi)巖漿活動(dòng)以五臺(tái)期和燕山期為主,與該區(qū)成礦作用密切相關(guān)的是強(qiáng)烈而頻繁的燕山期巖漿活動(dòng)(周利霞和唐耀林,1997)。中生代早中期,區(qū)域巖石圈處于歐亞板塊碰撞后強(qiáng)烈加厚狀態(tài),區(qū)域性隆升和剝蝕顯著(劉紅濤等,2002)。華北中部巖石圈在造山帶基礎(chǔ)上開(kāi)始減薄過(guò)程,主要始于130 ~110Ma 的拆沉減薄(羅照華等,2006;Li and Santosh,2014)。巖石圈減薄以及軟流圈物質(zhì)上涌導(dǎo)致該區(qū)的殼型、殼?;旌闲椭兴嵝?酸性巖漿活動(dòng),并侵入形成大規(guī)?;◢弾r類(Li et al.,2013,2014)。

    礦區(qū)內(nèi)地層由老到新自南東向北西依次出露有寒武系中上統(tǒng)、奧陶系、侏羅系及第四系(圖2),巖層走向30 ~50°,傾向NW,傾角10 ~50°,其中寒武系中上統(tǒng)泥質(zhì)白云巖、條帶狀灰?guī)r和奧陶系下統(tǒng)冶里組、亮甲山組泥晶灰?guī)r、細(xì)晶白云巖,中統(tǒng)馬家溝組角礫狀灰?guī)r,是刁泉復(fù)式巖體和銀銅礦體的圍巖(李兆龍和張連營(yíng),1999)。刁泉礦區(qū)的斷裂主要有北西向的刁泉斷裂和北東向的槍頭嶺斷裂,為成礦提供了有利的運(yùn)移通道和賦存空間。

    與礦化有關(guān)的中酸性雜巖體刁泉巖體呈巖株?duì)町a(chǎn)出,侵位于兩組斷裂構(gòu)造的交會(huì)部位,直徑約900m,面積0.7km2,巖性以花崗斑巖為主(圖4c),并有中基性脈巖發(fā)育。李兆龍和張連營(yíng)(1999)測(cè)定刁泉巖體花崗斑巖和石英斑巖全巖Rb-Sr 等時(shí)線年齡為130.9 ±1.4Ma,黑云母石英二長(zhǎng)巖的黑云母K-Ar 同位素年齡為130.5Ma,易洪波(2012)用單顆粒鋯石U-Pb 法測(cè)定黑云母石英二長(zhǎng)巖及花崗斑巖的年齡分別為163.7 ±1.8Ma 和154.5 ±1.1Ma,均屬燕山期產(chǎn)物。巖體主要巖石屬于鋁過(guò)飽和型,呈鈣堿性,稀土元素分布屬于輕稀土富集型,成因類型屬殼?;煸葱?周利霞和唐耀林,1997;李兆龍和張連營(yíng),1999)。易洪波(2012)對(duì)礦體中石英細(xì)脈型輝鉬礦礦石進(jìn)行了Re-Os 同位素測(cè)定,得到礦床形成年齡為132.3 ±6.1Ma,而張聚全對(duì)其測(cè)定的Re-Os 模式年齡為130.9 ±1.6Ma(未發(fā)表數(shù)據(jù))。

    刁泉巖體與圍巖的環(huán)形接觸帶是本區(qū)礦體的賦存部位,Ⅰ礦帶為礦區(qū)主要礦帶,礦體數(shù)量多達(dá)33 個(gè),賦存于矽卡巖帶及其外接觸帶中,具有工業(yè)意義的礦體主要有3 號(hào)、8 號(hào)及25 號(hào)三個(gè)礦體。其中25 號(hào)礦體規(guī)模最大,延長(zhǎng)大于200m,延深100 ~300m,厚度一般1 ~15m,最大厚度40 余米。礦體在平面上主要呈脈狀、似層狀,透鏡狀、囊狀,在剖面上以彎月?tīng)?、鉤狀產(chǎn)出,分枝復(fù)合、尖滅再現(xiàn)現(xiàn)象明顯,總體產(chǎn)狀從上到下先內(nèi)傾而后外傾,垂向上顯示上富Ag 下富Cu、Au 的礦化規(guī)律(李兆龍和張連營(yíng),1999)。

    蝕變分帶從巖體到圍巖可分為:Ⅰ內(nèi)蝕變帶,巖體不同程度地發(fā)生絹云母化、鈉長(zhǎng)石化、硅化、綠泥石化、矽卡巖化和碳酸巖化,0 ~20m 寬,主要金屬礦化為輝鉬礦化,局部成為工業(yè)小礦體;Ⅱ矽卡巖帶,主要由透輝石-鈣鐵榴石矽卡巖、綠簾石矽卡巖、鈣鋁榴石矽卡巖組成,2 ~60m 寬,銀銅礦體主要賦存于此帶中;Ⅲ矽卡巖化大理巖帶,0 ~50m,巖體發(fā)生透輝石化、鈣鐵榴石化、鈣鋁榴石化,局部有銀銅礦體賦存;Ⅳ大理巖、角巖帶,100 ~500m 寬。

    巖漿作用伴隨而來(lái)的熱液活動(dòng)使該區(qū)圍巖遭受不同程度的變質(zhì)作用,具體可分為兩個(gè)類型:一是與熱變質(zhì)有關(guān)的大理巖化和角巖化;二是與熱液交代有關(guān)的絹云母化、鈉長(zhǎng)石化、鉀長(zhǎng)石化、硅化、矽卡巖化和碳酸鹽化等,其中矽卡巖化與本區(qū)成礦關(guān)系最為密切(王永明等,2006)。

    礦石類型以脈狀(圖3b)、浸染狀(圖3c,d)、塊狀為主。礦石結(jié)構(gòu)主要有有他形粒狀結(jié)構(gòu)(圖3d)、填隙結(jié)構(gòu)(圖4d)、自形-半自形粒狀結(jié)構(gòu)(圖4d)、乳滴狀固熔體分離結(jié)構(gòu)(圖4f)、文象結(jié)構(gòu)(圖4f)和共邊結(jié)構(gòu)等。礦石礦物主要為輝銀礦、碲銀礦、自然銀、金銀礦、黃銅礦、斑銅礦、輝銅礦、輝鉬礦等,脈石礦物主要為石榴子石、方解石、白云石等。

    圖3 刁泉銀銅礦典型手標(biāo)本特征(a)-矽卡巖,石榴子石晶形較好,發(fā)育輕微綠泥石化,方解石沿裂隙充填;(b)-石英細(xì)脈穿切花崗巖,伴隨輕微輝鉬礦化;(c)-浸染狀斑銅礦膠結(jié)石榴子石矽卡巖角礫;(d)-斑銅礦和輝銅礦呈稠密浸染狀與方解石伴生,方解石粒度較大. Grnt-花崗巖;Sk-矽卡巖;Grt-石榴子石;Phl-金云母;Qz-石英;Cal-方解石;Mo-輝鉬礦;Bn-斑銅礦;Cc-輝銅礦Fig. 3 Typical hand specimens of the Diaoquan Ag-Cu deposit(a)-crystal shapes of garnets are relatively complete in skarn with slight chloritization,and calcite distributes along the fissures;(b)-a quartz vein with slight molybdenite cuts through the granite;(c)-the garnet skarn breccias are cemented with disseminated bornite;(d)-dense disseminated bornite and chalcocite are born associating with coarse grained calcite. Grnt-granite;Sk-skarn;Grt-garnet;Phlphlogopite;Qz-quartz;Cal-calcite;Mo-molybdenite;Bn-bornite;Cc-chalcocite

    2 成礦期次劃分和主要礦物特征

    2.1 成礦期次劃分

    根據(jù)野外脈體穿插關(guān)系、礦物組合、空間分布及結(jié)構(gòu)構(gòu)造特點(diǎn),將刁泉銀銅礦的成礦作用從早到晚劃分為3 個(gè)成礦期5 個(gè)成礦階段:Ⅰ矽卡巖期,可分為兩個(gè)階段,Ⅰ1 早矽卡巖階段和Ⅰ2 晚矽卡巖階段,Ⅰ1 階段最常見(jiàn)礦物為石榴子石,輝石、硅灰石等無(wú)水硅酸鹽,其中石榴子石含量較高且自形程度好,透輝石含量相對(duì)較少,Ⅰ2 階段主要形成符山石、透閃石、陽(yáng)起石、綠簾石等含水硅酸鹽礦物;Ⅱ石英-硫化物期,Ⅱ1 鉬礦化階段主要金屬礦物是輝鉬礦,Ⅱ2 銅銀礦化階段以石英、綠泥石等非金屬礦物和斑銅礦、黃銅礦、黃鐵礦等金屬硫化物為主,并伴隨著銀礦物和金礦物的沉淀;Ⅲ碳酸鹽期,這一時(shí)期主要發(fā)育碳酸鹽化,并伴隨銅礦化和銀礦化。

    2.2 主要矽卡巖礦物

    圖4 刁泉銀銅礦床巖相學(xué)特征(a)-具明顯環(huán)帶的石榴子石間隙被方解石充填(單偏光);(b)-早期黑云母發(fā)育綠泥石化(單偏光);(c)-花崗斑巖,斑晶為石英和鉀長(zhǎng)石(正交偏光);(d)-晚期的黃銅礦、黃鐵礦、方解石細(xì)脈充填石榴子石間隙(反射光);(e)-方鉛礦、閃鋅礦、黃鐵礦共生(反射光);(f)-斑銅礦與輝銅礦呈文象結(jié)構(gòu),閃鋅礦中分布星點(diǎn)狀黃銅礦(反射光);(g、h)-黃銅礦中分布輝銀礦、碲銀礦等(背散射圖像). Bi-黑云母;Chl-綠泥石;Or-鉀長(zhǎng)石;Sh-白鎢礦;Ccp-黃銅礦;Py-黃鐵礦;Gn-方鉛礦;Sp-閃鋅礦;Arg-輝銀礦;Hes-碲銀礦Fig.4 Petrographical characteristics of samples from the Diaoquan Ag-Cu deposit(a)-garnets with clear zoning are filled with calcite (under transmitted light);(b)-biotite formed in an early stage is partly altered with choritization (under transmitted light);(c)-graniteporphyry with quartz and K-feldspar phenocrysts (under transmitted light);(d)-later chalcopyrite,pyrite and calcite veinlet fill the gaps of earlier garnets (under reflective light);(e)-galena and sphalerite are associated with pyrite (under reflective light);(f)-graphic structure formed by bornite and chalcocite,star-like chalcopyrite distributing in sphalerite (reflected light);(g,h)-chalcopyrite with small argentite and hessite (BSE image). Bi-biotite;Chl-chlorite;Or-orthoclase;Tr-tremolite;Sh-scheelite;Ccp-chalcopyrite;Pypyrite;Gn-galena;Sp-sphalerite;Arg-argentite;Hes-hessite

    刁泉銀銅礦的石榴子石多呈淺黃綠色(圖3a),單偏光下常為淺黃褐色,極正高突起(圖4a),正交偏光下異常干涉色可達(dá)一級(jí)灰白,自形-他形,為鈣鐵榴石-鈣鋁榴石系列,環(huán)帶結(jié)構(gòu)發(fā)育。黑云母呈褐色,多色性明顯,常見(jiàn)綠泥石化(圖4b);金云母多呈深綠色(圖3d),薄片中無(wú)色至淺棕色,近平行消光。巖體中的云母以黑云母為主而矽卡巖中則以金云母為主,黑云母處于鎂質(zhì)黑云母與鐵質(zhì)黑云母的過(guò)渡位置。

    綠簾石呈黃綠色,多色性弱,二級(jí)干涉色鮮艷而明亮;輝石以透輝石為主,多呈綠色或暗綠色,短柱狀,正交偏光下可見(jiàn)鮮艷的二級(jí)干涉色;角閃石呈草綠色,可見(jiàn)典型的閃石式解理。

    2.3 主要金屬礦物

    斑銅礦在刁泉銀銅礦的分布十分廣泛,多與輝銅礦、黃鐵礦、黃銅礦等硫化物共生,礦石構(gòu)造以稠密浸染狀(圖3c,d)、塊狀為主,其次為團(tuán)塊狀和網(wǎng)脈狀構(gòu)造。多呈他形-半自形結(jié)構(gòu),固溶體分離結(jié)構(gòu)、共邊結(jié)構(gòu)、交代殘余結(jié)構(gòu)十分發(fā)育,固溶體分離結(jié)構(gòu)主要表現(xiàn)為斑銅礦呈乳滴狀或葉片狀產(chǎn)于輝銅礦(圖4f)。

    刁泉礦區(qū)黃鐵礦并不十分發(fā)育,呈半自形-自形(圖4e),常見(jiàn)填隙結(jié)構(gòu)或呈脈狀充填于石榴子石等早期形成的矽卡巖礦物之間。黃銅礦呈銅黃色,他形不規(guī)則塊狀,呈填隙結(jié)構(gòu)充填于早期形成的礦物之間(圖4d),與閃鋅礦呈固溶體分離結(jié)構(gòu),薄片中常呈星點(diǎn)狀(圖4f)。方鉛礦主要呈半自形-他形,閃鋅礦主要呈褐紅色調(diào),常與黃銅礦呈固溶體分離結(jié)構(gòu),與黃鐵礦、黃銅礦等其他硫化物共生(圖4e,f)。

    3 流體包裹體研究

    3.1 樣品及測(cè)試方法

    用于研究的樣品兼顧了與成礦有密切關(guān)系的各類地質(zhì)體,主要針對(duì)形態(tài)較為規(guī)則的原生流體包裹體進(jìn)行研究。主要巖性包括石榴子石矽卡巖,含石英脈花崗巖、蝕變花崗斑巖、黑云母石英二長(zhǎng)巖以及礦石,包裹體寄主礦物主要是石榴子石、石英和方解石,分別來(lái)自石榴子石矽卡巖,石英脈和碳酸鹽脈(圖3a,b,d)。野外共采集樣品33 件,根據(jù)產(chǎn)出特征和手標(biāo)本觀察結(jié)果從中挑選17 件樣品磨制成0.3mm厚雙面拋光的包裹體光薄片,在進(jìn)行流體包裹體巖相學(xué)觀察后,挑選出合適的測(cè)溫樣品,用丙酮浸泡并清洗薄片,然后進(jìn)行均一法-冷凍法測(cè)溫研究及成分分析。流體包裹體測(cè)溫工作在中國(guó)地質(zhì)大學(xué)(北京)礦物標(biāo)型實(shí)驗(yàn)室進(jìn)行,所用儀器為德國(guó)Leica DM2500P 偏光顯微鏡,英國(guó)Linkam THMSG600 冷熱臺(tái),其可測(cè)溫度范圍為- 196 ~600℃,測(cè)試精度為±1.0℃。加熱、冷凍過(guò)程中設(shè)置的初始控溫速率為30℃/min,開(kāi)始觀察到變化時(shí)降低速率,相變點(diǎn)附近速率一般為1 ~2℃/min。

    單個(gè)流體包裹體激光拉曼顯微探針測(cè)試在核工業(yè)北京地質(zhì)研究院分析測(cè)試研究中心完成。測(cè)試儀器為L(zhǎng)ABHRVIS LabRAM HR800 研究級(jí)顯微激光拉曼光譜儀,使用Yag晶體倍頻固體激光器,波長(zhǎng)為532nm,掃描范圍100 ~4200(cm-1)。

    3.2 流體包裹體類型

    通過(guò)對(duì)刁泉銀銅礦中流體包裹體進(jìn)行巖相學(xué)觀察,發(fā)現(xiàn)在三種被選作主要研究對(duì)象的寄主礦物中石英和方解石包裹體均較發(fā)育,石榴子石中包裹體數(shù)量相對(duì)較少。根據(jù)包裹體在室溫下的相態(tài)可以將其分為3 類5 個(gè)亞類(圖5),具體分述如下。

    圖5 刁泉銀銅礦流體包裹體顯微照片(a、b)-石榴子石中的1 型包裹體,(b)中可見(jiàn)晶體的生長(zhǎng)紋;(c)-方解石中呈負(fù)晶形的1 型包裹體;(d、e)-石英中1 型包裹體,(e)為原生包裹體集群分布;(f)-石英中1 型包裹體和2B 型包裹體密切共生,顯示流體的沸騰特征;(g)-石英中的2A 型包裹體;(h)-石英中的2B 型包裹體;(i)-石英中的3A 型包裹體;(j)-石榴子石中的3B 型包裹體;(k)-石榴子石中含有兩個(gè)子晶的3B 型包裹體Fig.5 Micrographs of fluid inclusions in the Diaoquan Ag-Cu deposit(a,b)-type 1 inclusions in garnet,growth striation can be seen in(b);(c)-a negative crystal shaped type 1 inclusion in calcite;(d,e)-type 1 inclusions in quartz. (e)presents a primary inclusion group;(f)-type 1 inclusions and type 2B inclusions are in close association within quartz,displaying boiling characteristics of the fluid;(g)-a type 2A inclusion in quartz;(h)-type 2B inclusion in quartz;(i)-a type 3A inclusion in quartz;(j)-type 3B inclusions in garnet;(k)-a type 3B inclusion containing two daughter minerals (a halite and a sylvite)

    圖6 刁泉銀銅礦床流體包裹體激光拉曼圖譜Fig. 6 Laser Raman spectra of fluid inclusions of the Diaoquan Ag-Cu deposit

    1 型富液相包裹體:室溫下氣相充填度小于50%,加熱后均一至液相。1 型包裹體在石榴子石、石英、方解石中均占重要地位,是該礦區(qū)最為發(fā)育的包裹體類型,石榴子石中的許多該類包裹體表現(xiàn)出臨界均一或近臨界均一行為。

    2 型富氣相包裹體:室溫下氣相充填度大于50%,只在石英中發(fā)現(xiàn),數(shù)量較少。依據(jù)其均一方式可以將該類型包裹體劃分為兩個(gè)亞類:2A 型,加熱后均一為液相,許多該類包裹體表現(xiàn)出近臨界均一;2B 型,加熱后均一為氣相。

    3 型含子礦物包裹體:室溫下存在一個(gè)或一個(gè)以上子礦物的三相包裹體,集群或星散狀出現(xiàn)。子礦物類型主要為石鹽和鉀鹽,石鹽子礦物多呈立方體,鉀鹽子礦物呈立方體或渾圓狀,均無(wú)色透明。由于溫度升高時(shí)鉀鹽溶解度增加的速度遠(yuǎn)大于石鹽,因此鉀鹽子礦物先于石鹽子礦物溶解。根據(jù)其子礦物的消失溫度與氣液均一溫度的關(guān)系將其分為:3A型,氣液均一溫度高于子礦物溶化溫度,只出現(xiàn)在石榴子石;3B 型,氣液均一溫度低于子礦物熔化溫度,只出現(xiàn)于石英,多呈負(fù)晶形。

    本文對(duì)石榴子石和石英流體包裹體進(jìn)行了激光拉曼顯微探針?lè)治?圖6)。方解石由于其較大的雙折率,未對(duì)其進(jìn)行激光拉曼測(cè)試。石榴子石流體包裹體氣相的測(cè)點(diǎn)除寄主礦物的特征峰以外,未檢測(cè)到其它成分,可能由于含量太低所致。石英中流體包裹體氣相主要成分為H2O,少量CO2,液相成分僅檢測(cè)到寬緩的水峰。因此將刁泉銀銅礦流體歸為NaCl-H2O 體系進(jìn)行計(jì)算和討論較為合理。

    3.3 不同礦化階段流體包裹體特征及均一溫度

    本次研究選取矽卡巖期石榴子石、石英-硫化物期石英和碳酸鹽期方解石3 種易于觀察且包裹體發(fā)育較為理想的主礦物,對(duì)各個(gè)類型的包裹體進(jìn)行了顯微測(cè)溫,結(jié)果見(jiàn)表1和圖7。

    3.3.1 矽卡巖期

    矽卡巖期的流體包裹體主要賦存于石榴子石,形成于早期矽卡巖階段,其中包裹體的捕獲溫度可以較好地代表這一階段的形成溫度。石榴子石中包裹體類型為1 型和3A 型,其中1 型包裹體最為發(fā)育,不同區(qū)域包裹體豐度變化較大,大小以3 ~10μm 為主。

    這一階段的1 型包裹體均一溫度變化范圍為339 ~587℃,均值為464℃。許多包裹體表現(xiàn)出臨界均一或接近臨界均一,即氣泡跳動(dòng)不明顯且臨近均一時(shí)氣液邊界消失,而非氣泡逐漸變小至消失而達(dá)到均一。含子礦物3A 型包裹體升溫后石鹽先于氣泡消失,完全均一溫度的變化范圍481 ~589℃,平均值為528℃,這類包裹體被認(rèn)為是從原始不飽和溶液中捕獲的(盧煥章等,2004)。

    3.3.2 石英-硫化物期

    石英-硫化物期的包裹體主要賦存于石英之中,包裹體種類多,且豐度較大,主要類型為1 型,其次為3B 型,少量2A、2B 型。大小變化較大,一般在3 ~15μm,形狀各異,有負(fù)晶形、渾圓形等等。1 型包裹體均一溫度變化范圍較大,為170 ~510℃,平均280℃,峰值出現(xiàn)于200 ~240℃區(qū)間。2A和2B 型包裹體統(tǒng)計(jì)情況類似,都分布于360 ~420℃區(qū)間。觀察到的多數(shù)2 型包裹體表現(xiàn)出接近臨界均一,且大多出現(xiàn)于含輝鉬礦的石英脈。1 型和2 型包裹體均一溫度均覆蓋了360 ~420℃區(qū)間,暗示了沸騰的發(fā)生。3B 型包裹體升溫后氣泡先消失,通過(guò)石鹽溶解達(dá)到均一狀態(tài),這類包裹體捕獲于原始過(guò)飽和溶液(盧煥章等,2004),最終均一溫度為259~483℃,平均值為351.4℃。因此,流體在360 ~420℃區(qū)間流體發(fā)生了沸騰導(dǎo)致礦化,此為Ⅱ1 輝鉬礦化階段,而Ⅱ2 銀銅礦化階段發(fā)生在200 ~240℃區(qū)間(圖7)。

    表1 刁泉銀銅礦床流體包裹體數(shù)據(jù)表Table 1 Microthermometric data for fluid inclusions from the Diaoquan Ag-Cu deposit

    圖7 刁泉銀銅礦不同類型包裹體均一溫度直方圖Fig.7 Histogram showing homogenization temperature distribution for different types of fluid inclusions in the Diaoquan Ag-Cu deposit

    3.3.3 碳酸鹽期

    碳酸鹽期仍然伴隨銅銀礦化,大量方解石從含豐富碳酸鹽組份的溶液中結(jié)晶出來(lái),包裹體主要賦存其中。包裹體種類單一,僅發(fā)現(xiàn)1 型,大小多在4 ~12μm 之間,個(gè)體較小者多集群而較大者星散狀出現(xiàn),均一溫度分布于133 ~323℃,平均值為226℃,其均一溫度頻數(shù)的峰值出現(xiàn)區(qū)間與石英中1型包裹體基本吻合,這印證了銀銅礦化發(fā)生的溫度區(qū)間為200 ~240℃。

    3.4 鹽度

    3 型含子礦物包裹體均含有石鹽子晶,僅極個(gè)別包裹體觀察到鉀鹽子晶,因此按照NaCl-H2O 體系進(jìn)行鹽度計(jì)算。1型包裹體的鹽度利用冰點(diǎn)溫度求得,3 型包裹體的鹽度利用子礦物融化溫度求得(Hall et al.,1988)。2 型包裹體由于冷凍現(xiàn)象不明顯,其冰點(diǎn)溫度未測(cè)出。1 型和3 型包裹體計(jì)算結(jié)果見(jiàn)表1 和圖8。

    該礦床流體包裹體的鹽度(%NaCleqv)整體變化范圍很大,介于0.2% ~64.0%,1 型包裹體為低鹽度而3 型包裹體為中高鹽度。石榴子石和石英中包裹體的鹽度可以分為兩個(gè)區(qū)間,這可能是由于子礦物處于亞穩(wěn)態(tài)造成(Shepherd et al.,1985),也可能是流體相分離的結(jié)果(Heinrich,2005),而方解石則相對(duì)集中于低鹽度。1 型包裹體中,石榴子石鹽度最高,集中于10% ~20%,石英鹽度變化范圍最大,方解石鹽度最低(<10%)。對(duì)于3 型包裹體,從石榴子石中的3A型到石英中的3B 型穿越了NaCl 飽和曲線,并大致沿NaCl飽和曲線演化。3A 型包裹體的鹽度分布于50% ~60%,而3B 型包裹體的鹽度集中在30% ~50%。

    圖8 刁泉銀銅礦流體包裹體均一溫度-鹽度雙變量圖解虛線指在一定溫度和鹽度下NaCl 溶液的的蒸氣壓力(MPa)(據(jù)Calagari,2004)Fig.8 Homogenization temperature-salinity phase diagram of fluid inclusions in the Diaoquan Ag-Cu depositDashed lines referring to vapor pressures of NaCl solutions at the indicated temperatures and salinity (after Calagari,2004)

    表2 刁泉銀銅礦流體包裹體密度、壓力和捕獲深度估計(jì)Table 2 Density,pressure and the estimation of depth when captured

    圖9 刁泉銀銅礦流體包裹體密度分布直方圖Fig. 9 Histogram showing density distribution of fluid inclusions in the Diaoquan Ag-Cu deposit

    從石榴子石到石英、方解石,其中包裹體的溫度和鹽度不斷下降,反映了成礦過(guò)程中的由高溫到低溫,自高鹽度向低鹽度的流體演化趨勢(shì)。

    3.5 密度

    流體包裹體密度的計(jì)算根據(jù)氣液均一溫度和鹽度算出(劉斌和段光賢,1987;劉斌,2001)。統(tǒng)計(jì)結(jié)果(圖9 和表2)表明,石榴子石中流體包裹體密度分布于兩個(gè)區(qū)間,分別為0.4 ~0.6g/cm3和0.7 ~1.0g/cm3,而石英密度分布范圍較廣,從0.3 ~1.2g/cm3均有分布,峰值出現(xiàn)在0.4 ~0.5g/cm3,0.6 ~0.8g/cm3和1.0 ~1.1g/cm3和0.7 ~1.0g/cm3。方解石密度分布最為集中,介于0.6 ~0.8g/cm3,峰值出現(xiàn)在0.7 ~0.8g/cm3,總體密度介于石榴子石和石英之間。石英的密度高于方解石可能與石英-硫化物階段流體含礦性較高有關(guān),而方解石中雖仍伴隨礦化,但在碳酸鹽期流體中的礦質(zhì)逐漸沉淀,流體活動(dòng)趨于結(jié)束,因此方解石中流體包裹體的密度分布較為集中。

    圖10 刁泉銀銅礦包裹體在NaCl-H2O 體系P-X-T 投影圖(據(jù)Shinohara and Hedenquist,1997;Hedenquist et al.,1998;Wilkinson,2001)Fig.10 P-X-T phase diagram of fluid inclusions in NaCl-H2O system with projections of fluid inclusions from the Diaoquan Ag-Cu deposit on it (after Shinohara and Hedenquist,1997;Hedenquist et al.,1998;Wilkinson,2001)

    3.6 壓力和深度

    進(jìn)行壓力測(cè)算的包裹體均屬于NaCl-H2O 體系,因此在NaCl-H2O 體系P-X-T 圖解(圖10)上求得了各階段包裹體的最低捕獲壓力和深度,并參考了其他計(jì)算方法(Driesner and Heinrich,2007)。求得矽卡巖期流體包裹體最低捕獲壓力為23 ~66MPa,石英-硫化物期5 ~46MPa,碳酸鹽期<15MPa。若采用10MPa/km(Zoback and Zoback,1989)的靜水壓力和25MPa/km 的靜巖壓力(Grandi et al.,2002)計(jì)算各階段流體的深度,矽卡巖期按照靜巖壓力計(jì)算對(duì)應(yīng)的流體深度為0.9~2.6km,石英-硫化物期按照靜水壓力計(jì)算的流體深度為0.5 ~4.6km,Mo 礦化主要發(fā)生在1 ~3km,而Cu-Ag 礦化主要發(fā)生在0.5 ~1.5km,碳酸鹽期靜水壓力捕獲深度<1.5km。

    4 討論

    華北克拉通在燕山期突然發(fā)生強(qiáng)烈而廣泛的巖漿-構(gòu)造-成礦事件(鄧晉福等,1999)。刁泉銀銅礦床形成于燕山晚期,在此時(shí)期,巖石圈減薄,軟流圈物質(zhì)上涌,殼型、殼幔混合型中酸性-酸性巖漿沿基底斷裂上升侵位(Li and Santosh,2014;Li et al.,2013,2014)。在巖漿侵位的最后階段,結(jié)晶作用使流體中開(kāi)始出現(xiàn)高溫氣液(Candela and Piccoli,1995)。

    刁泉礦床受火山口構(gòu)造影響,加上斷裂構(gòu)造為熱液提供了運(yùn)移通道和上升空間,富含成礦物質(zhì)的流體向聚礦構(gòu)造快速運(yùn)移,于巖漿結(jié)晶的最后階段從淺部出溶的流體在矽卡巖期主礦物石榴子石中形成1 型富液相和3 型含子礦物包裹體。出溶水溶液的鹽度越高其攜帶成礦金屬元素能力越強(qiáng)(于淼等,2013),因此從熔體相中分離出來(lái)的高溫高鹽度流體為主要載礦流體。

    隨著流體演化至石英-硫化物期,在大致400℃左右,封閉體系逐漸變?yōu)殚_(kāi)放體系,壓力降低導(dǎo)致原始均勻流體減壓沸騰,造成大量酸性氣體組分逸出,從而使流體pH 值增大,有利于成礦物質(zhì)沉淀(Skinner,1979)。這與激光拉曼研究中石榴子石流體包裹體沒(méi)有檢測(cè)到CO2而石英中流體包裹體氣相成分具有CO2特征峰的結(jié)果一致。流體溫度壓力條件發(fā)生突變導(dǎo)致許多絡(luò)合物不再穩(wěn)定,也致使大量成礦元素沉淀并富集成礦(芮宗瑤等,2003)。石英中的3 型含子礦物包裹體穿越石鹽飽和曲線,溫度和壓力條件突變導(dǎo)致流體發(fā)生沸騰,此時(shí)在石英中出現(xiàn)1 型富氣相與2 型富液相包裹體共生的現(xiàn)象,這是在石英-硫化物期的早期形成鉬礦化的直接動(dòng)因。熱液系統(tǒng)在350℃以下尤其在300 ~180℃之間是金屬硫化物的溶解度從高到低的轉(zhuǎn)變階段(Helgeson,1969),因此石英-硫化物期后期和碳酸鹽期銅和銀的硫化物大量生成。刁泉銀銅礦床的Mo 礦化主要形成于1 ~3km,Cu-Ag 礦化主要形成于0.5 ~1.5km,與矽卡巖礦床一般形成深度1 ~4.5km 之間(趙一鳴和林文蔚,1990)的結(jié)論基本吻合。

    5 結(jié)論

    (1)刁泉銀銅礦成礦過(guò)程從早到晚可以劃分為3 個(gè)成礦期5 個(gè)成礦階段,可見(jiàn)典型的矽卡巖礦物組合。Ⅰ矽卡巖期,包括Ⅰ1 早矽卡巖階段和Ⅰ2 晚矽卡巖階段;Ⅱ石英-硫化物期,可分為Ⅱ1 鉬礦化階段和Ⅱ2 銅銀礦化階段;Ⅲ碳酸鹽期。

    (2)刁泉銀銅礦流體包裹體類型主要包括富液相、富氣相和含子晶包裹體,成礦流體為NaCl-H2O 體系。淺成條件下出溶的巖漿發(fā)生不混溶作用,矽卡巖期的流體由高溫(481~587℃)高鹽度(51% ~64% NaCleqv)和較高溫(339 ~587℃)中鹽度(11% ~18%NaCleqv)的流體組成。石英-硫化物期在鉬礦化階段發(fā)生減壓沸騰作用,形成富氣相與富液相包裹體共生,流體的溫度(170 ~510℃)鹽度(35% ~58%NaCleqv、2% ~20%NaCleqv)與矽卡巖期相比有明顯下降,隨后伴隨銅銀礦物沉淀。碳酸鹽期流體溫度(132 ~323℃)與鹽度(<8%NaCleqv)繼續(xù)降低且分布集中。

    (3)石榴子石中流體包裹體密度分布于0.4 ~0.6g/cm3和0.7 ~1.0g/cm3,石英密度從0.3 ~1.2g/cm3均有分布,方解石密度分布最為集中,介于0.6 ~0.8g/cm3。矽卡巖期流體包裹體最低捕獲壓力為23 ~66MPa,按照靜巖壓力計(jì)算對(duì)應(yīng)的流體深度為0.9 ~2.6km;石英-硫化物期5 ~46MPa,按照靜水壓力計(jì)算的流體深度為0.5 ~4.6km,Mo 礦化主要發(fā)生在1 ~3km,而Cu-Ag 礦化主要發(fā)生在0.5 ~1.5km;碳酸鹽期<15MPa,靜水壓力捕獲深度<1.5km。刁泉銀銅礦為一與淺成低溫?zé)嵋河嘘P(guān)的矽卡巖型礦床。

    致謝 野外工作得到山西紫金礦業(yè)公司各位領(lǐng)導(dǎo)和技術(shù)人員的大力支持,課題組張雪鴿、賈利賓、李增達(dá)給予了無(wú)私的幫助;實(shí)驗(yàn)室工作得到張秀寶老師、歐光習(xí)老師的指導(dǎo)和幫助;羅照華教授對(duì)初稿進(jìn)行了審閱,匿名審稿專家對(duì)本文提出許多寶貴意見(jiàn);在此一并表示誠(chéng)摯感謝。鄧晉福教授長(zhǎng)期給予通訊作者多方面的教誨和學(xué)術(shù)熏陶,本文特別為先生八十壽辰所作,祝先生健康長(zhǎng)壽,繼續(xù)指引我們的團(tuán)隊(duì)奮勇前行。

    Calagari AA. 2004. Fluid inclusion studies in quartz veinlets in the porphyry copper deposit at Sungun,East-Azarbaidjan,Iran. Journal of Asian Earth Sciences,23(2):179 -189

    Candela PA and Piccoli PL. 1995. Model ore-metal partitioning from melt into vapor and vapor/brine mixtures. In:Thompson JFH (ed.).Magmas,F(xiàn)luid and Ore Deposits. Mineralogical Association of Canada Short Course Handbook,23:101 -127

    Deng JF,Mo XX,Zhao HL,Luo ZH,Zhao GC and Dai SQ. 1999. The Yanshanian lithosphere-asthenosphere catastrophe and metallogenic environment in East China. Mineral Deposits,18(4):309 -315(in Chinese with English abstract)

    Deng JF,Su SG,Niu YL,Liu C,Zhao GC,Zhao XG,Zhou S and Wu ZX. 2007. A possible model for the lithospheric thinning of North China Craton:Evidence from the Yanshanian (Jura-Cretaceous)magmatism and tectonism. Lithos,96(1):22 -35

    Driesner T and Heinrich CA. 2007. The system H2O-NaCl. Part I:Correlation formulae for phase relations in temperature-pressurecomposition space from 0 to 1000℃,0 to 5000bar,and 0 to 1 XNaCl. Geochimica et Cosmochimica Acta,71(20):4880 -4901

    Fan HR,Xie YH and Wang YL. 1998. Fluid-rock interaction during mineralization of the Shanggong structure controlled alteration-type gold deposit in western Henan Province,central China. Acta Petrologica Sinica,14(4):529 - 541 (in Chinese with English abstract)

    Grandi S,Rao RVN and Toksoz MN. 2002. Geomechanical modeling of in-situ stresses around a borehole. Massachusetts Institute of Technology. Earth Resources Laboratory,1 -16

    Hall DL,Sterner SM and Bodnar RJ. 1988. Freezing point depression of NaCl-KCl-H2O solutions. Economic Geology,83(1):197 -202

    Hedenquist JW,Arribas A and Reynolds TJ. 1998. Evolution of an intrusion-centered hydrothermal system: Far Southeast-Lepanto porphyry and epithermal Cu-Au deposits,Philippines. Economic Geology,93(4):373 -404

    Heinrich CA. 2005. The physical and chemical evolution of low-salinity magmatic fluids at the porphyry to epithermal transition: A thermodynamic study. Mineralium Deposita,39(8):864 -889

    Helgeson HC. 1969. Thermodynamics of hydrothermal systems at elevated and pressures. American Journal of Science,367(7):729 -804

    Johnson TW and Meinert LD. 1994. Au-Cu-Ag skarn and replacement mineralization in the McLaren deposit,New World District,Park County,Montana. Economic Geology,89(5):969 -993

    Li AS and Han ZA. 2001. The role of geochemical exploration in the discovery of a large Ag-Au-Cu deposit at Diaoquan,Shanxi. Geology in China,28(11):29 -34 (in Chinese with English abstract)

    Li H,Zhang GY,Li DL,Yu B,Wang ZY,Xu LY and Chen JF. 2008.The structural superimposed halo model for the Diaoquan skarn silver-copper deposit in Lingqiu County, Shanxi Province.Geophysical and Geochemical Exploration,32(5):529 -533 (in Chinese with English abstract)

    Li SR,Santosh M,Zhang HF,Shen JF,Dong GC,Wang JZ and Zhang JQ. 2013. Inhomogeneous lithospheric thinning in the central North China Craton:Zircon U-Pb and S-He-Ar isotopic record from magmatism and metallogeny in the Taihang Mountains. Gondwana Reserch,23:141 -160

    Li SR and Santosh M. 2014. Metallogeny and craton destruction:Records from the North China Craton. Ore Geology Reviews,56:376 -414

    Li SR,Santosh M,Zhang HF,Luo JY and Zhang JQ. 2014. Metallogeny in response to lithospheric thinning and craton destruction:Geochemistry and U-Pb zircon chronology of the Yixingzhai gold deposit,central North China Craton. Ore Geology Reviews,56:457-471

    Li XG,Ge MR,Liu BT,Liu HF and Guo JL. 2004. Integrated prospecting model of Diaoquan copper-silver deposit. 2004. Journal of Taiyuan University of Technology,35(6):703 -706 (in Chinese with English abstract)

    Li XR,Liu XJ and Wang PY. 1999. A research on occurrence of gold and silver in ore of Diaoquan silver copper deposit. Contributions to Geology and Mineral Resources Research,14(2):55 - 62 (in Chinese with English abstract)

    Li ZL and Zhang LY. 1999. Geological characteristics and genetic mechanism of the Diaoquan silver copper deposit,Shanxi Province.Mineral Deposits,18 (1):11 - 21 (in Chinese with English abstract)

    Liu B and Duan GX. 1987. The density and isochoric formulae for NaCl-H2O fluid inclusions (salinity ≤25%)and their applications. Acta Mineralogica Sinica,7(4):345 - 352 (in Chinese with English abstract)

    Liu B. 2001. Density and isochoric formulae for NaCl-H2O inclusions with medium and high salinity and their applications. Geological Review,47(6):617 -622 (in Chinese with English abstract)

    Liu HT,Zhai MG,Liu JM and Sun SH. 2002. The Mesozoic granitoids in the northern marginal region of North China Craton:Evolution from post-collisional to anorogenic settings. Acta Petrologica Sinica,18(4):433 -448 (in Chinese with English abstract)

    Lu HZ,F(xiàn)ang HR,Ni P,Ou GX,Shen K and Zhang WH. 2004. Fluid Inclusion. Beijing:Science Press,1 -487 (in Chinese)

    Luo JY. 2009. Genetic mineralogy and metallogenic prognosis of Yixingzhai gold deposit in Fanshi County,Shanxi Province. Ph. D.Dissertation. Beijing:China University of Geosciences,1 -138 (in Chinese with English summary)

    Luo ZH,Deng JF and Han XQ. 1999. On Characteristics of Magmatic Activities and Orogenic Process of Taihangshan Intraplate Orogeny.Beijing:Geological Publishing House,1 -124 (in Chinese)

    Luo ZH,Wei Y,Xin HT,Ke S,Li WT,Li DD and Huang JX. 2006.The Mesozoic intraplate orogeny of the Taihang Mountains and the thinning of the continental lithosphere in North China. Earth Science Frontiers,13(6):52 -63 (in Chinese with English abstract)

    Morrison GW,Rose WJ and Jaireth S. 1991. Geological and geochemical controls on the silver content (fineness)of gold in gold-silver deposits. Ore Geology Reviews,6(4):333 -364

    Rui ZY,Zhao YM,Wang LS and Wang YT. 2003. Role of volatile components in formation of skarn and porphyry deposits. Mineral Deposits,22(1):141 -148 (in Chinese with English abstract)

    Shepherd T,Rankin AH and Alderton DHM. 1985. A Practical Guide to Fluid Inclusion Studies. Glasgow:Blackie & Son,239

    Shinohara H and Hedenquist JW. 1997. Constraints on magma degassing beneath the Far Southeast porphyry Cu-Au deposit,Philippines.Journal of Petrology,38(12):1741 -1752

    Skinner BJ. 1979. The many origins of hydrothermal mineral deposits.In:Barnes HL (ed.). Geochemistry of Hydrothermal Ore Deposit.New York:John Wiley and Sons,1 -12

    Tang YJ,Zhang HF,Ying JF and Zhang J. 2006. Source characteristics and temporal-spatial evolution of Mesozoic and Cenozoic basaltic magmatism in the Taihang Mountains. Acta Petrologica Sinica,22(6):1657 -1664 (in Chinese with English abstract)

    Wang YM,Yuan LD and Wen CG. 2006. Mineralogical characteristics of the Diaoquan contact metasomatic-medium to low temperature hydrothermal Ag-Cu deposit. Huabei Land and Resources,(1):4-7 (in Chinese)

    Wilkinson JJ. 2001. Fluid inclusions in hydrothermal ore deposits.Lithos,55(1):229 -272

    Yi HB. 2012. Geochemical characteristics of the Diaoquan Ag-Cu deposit and their constrains on genesis of the deposit in Lingqiu,Shanxi Province. Master Degree Thesis. Changsha: Central South University,1 -68 (in Chinese with English summary)

    Yu M,F(xiàn)eng CY,Xiao Y,Liu JN,Li DX,Li GC and Ma SC. 2013.Features and evolution of metallogenic fluid in Jiadanggen porphyry copper deposit of Gonghe Country,Qinghai Province. Mineral Deposits,32(1):133 -147 (in Chinese with English abstract)

    Zhao YM and Lin WW. 1990. Skarn Deposit of China. Beijing:Geological Publishing House,161 -180 (in Chinese)

    Zhao YM,Zhang YN and Bi CS. 1999. Geology of gold-bearing skarn deposits in the Middle and Lower Yangtze River Valley and adjacent regions. Ore Geology Reviews,14(3):227 -249

    Zhou LX and Tang YL. 1997. Geological characteristics and genesis of the Diaoquan silver-copper deposit,Shanxi Province. Jour. Geol. &Min. Res. North China,12(2):122 - 136 (in Chinese with English abstract)

    Zoback MD and Zoback ML. 1989. Stress in the earth’s lithosphere. In:Fairbridge RW (ed.). Encyclopedia of Earth Sciences Series. New York:Nostrand Reinhold Co.,1221 -1232

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

    鄧晉福,莫宣學(xué),趙海玲,羅照華,趙國(guó)春,戴圣潛. 1999. 中國(guó)東部燕山期巖石圈-軟流圈系統(tǒng)大災(zāi)變與成礦環(huán)境. 礦床地質(zhì),18(4):309 -315

    范宏瑞,謝奕漢,王英蘭. 1998. 豫西上宮構(gòu)造蝕變巖型金礦成礦過(guò)程中的流體-巖石反應(yīng). 巖石學(xué)報(bào),14(4):529 -541

    李愛(ài)生,韓志安. 2001. 化探在山西省刁泉大型銀金銅礦床發(fā)現(xiàn)中的作用. 中國(guó)地質(zhì),28(11):29 -34

    李惠,張國(guó)義,李德亮,禹斌,王志遠(yuǎn),許麗云,陳建飛. 2008. 山西靈丘刁泉矽卡巖型銀銅礦床構(gòu)造疊加暈?zāi)P? 物探與化探,32(5):529 -533

    李曉剛,葛民榮,劉邦濤,劉鴻福,郭俊麗. 2004. 刁泉銅銀礦床綜合找礦模式. 太原理工大學(xué)學(xué)報(bào),35(6):703 -706

    李喜榮,劉新江,王培英. 1999. 刁泉銀銅礦床金銀賦存狀態(tài)的研究. 地質(zhì)找礦論叢,14(2):55 -62

    李兆龍,張連營(yíng). 1999. 山西省刁泉銀銅礦床地質(zhì)特征及成因機(jī)制.礦床地質(zhì),18(1):11 -21

    劉斌,段光賢. 1987. NaCl-H2O 溶液包裹體的密度式和等容式及其應(yīng)用. 礦物學(xué)報(bào),7(4):345 -352

    劉斌. 2001. 中高鹽度NaCl-H2O 包裹體的密度式和等容式及其應(yīng)用. 地質(zhì)論評(píng),47(6):617 -622

    劉紅濤,翟明國(guó),劉建明,孫世華. 2002. 華北克拉通北緣中生代花崗巖:從碰撞后到非造山. 巖石學(xué)報(bào),18(4):433 -448

    盧煥章,范宏瑞,倪培,歐光習(xí),沈坤,張文淮. 2004. 流體包裹體.北京:科學(xué)出版社,1 -487

    羅軍燕. 2009. 山西省繁峙縣義興寨金礦床成因礦物學(xué)研究與成礦預(yù)測(cè). 博士學(xué)位論文. 北京:中國(guó)地質(zhì)大學(xué),1 -138

    羅照華,鄧晉福,韓秀卿. 1999. 太行山造山帶巖漿活動(dòng)及其造山過(guò)程反演. 北京:地質(zhì)出版社,1 -124

    羅照華,魏陽(yáng),辛后田,柯珊,李文韜,李德東,黃金香. 2006. 太行山中生代板內(nèi)造山作用與華北大陸巖石圈巨大減薄. 地學(xué)前緣,13(6):52 -63

    芮宗瑤,趙一鳴,王龍生,王義天. 2003. 揮發(fā)份在矽卡巖型和斑巖型礦床形成中的作用. 礦床地質(zhì),22(1):141 -148

    湯艷杰,張宏福,英基豐,張瑾. 2006. 太行山地區(qū)中、新生代玄武質(zhì)巖漿的源區(qū)特征與時(shí)空演化. 巖石學(xué)報(bào),22(6):1657 -1664

    王永明,袁立東,溫常貴. 2006. 刁泉接觸交代-中低溫?zé)嵋盒豌y銅礦床礦物學(xué)特征. 華北國(guó)土資源,(1):4 -7

    易洪波. 2012. 山西省靈丘縣刁泉銀銅礦地球化學(xué)特征及其對(duì)礦床成因的約束. 碩士學(xué)位論文. 長(zhǎng)沙:中南大學(xué),1 -68

    于淼,豐成友,肖曄,劉建楠,李大新,李國(guó)臣,馬圣鈔. 2013. 青海共和縣加當(dāng)根斑巖銅礦床成礦流體特征及演化. 礦床地質(zhì),32(1):133 -147

    趙一鳴,林文蔚. 1990. 中國(guó)矽卡巖礦床. 北京:地質(zhì)出版社,298 -317

    周利霞,唐耀林. 1997. 山西刁泉銀銅礦地質(zhì)特征及礦床成因. 華北地質(zhì)礦產(chǎn)雜志,12(2):122 -136

    猜你喜歡
    子石矽卡巖鹽度
    激電聯(lián)合剖面在判斷矽卡巖型礦床礦體產(chǎn)狀中的應(yīng)用
    鹽度和pH對(duì)細(xì)角螺耗氧率和排氨率的影響
    鹽度脅迫對(duì)入侵生物福壽螺的急性毒性效應(yīng)
    Geochemical and geochronological studies of the Aketas granite from Fuyun County,Xinjiang:the implications of the petrogenesis and tectonic setting
    Acta Geochimica(2015年4期)2015-10-25 02:03:16
    適用于高鹽度和致密巖層驅(qū)油的表面活性劑
    打括號(hào)的警察之墓
    牡丹(2014年10期)2014-10-17 09:31:59
    石榴子石Lu-Hf年代學(xué)研究進(jìn)展
    WorldView—Ⅱ在青海東昆侖求勉雷克塔格地區(qū)矽卡巖型礦床找礦中的應(yīng)用
    河北遙感(2014年3期)2014-07-10 13:16:47
    選礦尾礦中的石榴子石可選性試驗(yàn)研究
    膠州灣夏季鹽度長(zhǎng)期輸運(yùn)機(jī)制分析
    中文字幕免费在线视频6| 国产欧美日韩精品亚洲av| 日韩人妻高清精品专区| 国产黄片美女视频| 亚洲不卡免费看| 97碰自拍视频| 国产伦一二天堂av在线观看| 国产色婷婷99| 99精品在免费线老司机午夜| www.熟女人妻精品国产| 久久午夜福利片| 精品午夜福利在线看| 91麻豆精品激情在线观看国产| 免费电影在线观看免费观看| 日本精品一区二区三区蜜桃| 亚洲无线在线观看| 制服丝袜大香蕉在线| 日韩高清综合在线| 一级毛片久久久久久久久女| 白带黄色成豆腐渣| 国产视频一区二区在线看| 99久久久亚洲精品蜜臀av| 精品久久久久久,| 久久久成人免费电影| 国产一区二区亚洲精品在线观看| 亚洲人成网站在线播| 国产亚洲av嫩草精品影院| 亚洲美女视频黄频| 男人舔奶头视频| 国内少妇人妻偷人精品xxx网站| 欧美又色又爽又黄视频| 中文字幕av成人在线电影| 欧美一区二区国产精品久久精品| 老女人水多毛片| 99riav亚洲国产免费| 99riav亚洲国产免费| 亚洲国产色片| av视频在线观看入口| 国产精品三级大全| 我的女老师完整版在线观看| 国产精品久久久久久亚洲av鲁大| 亚洲中文字幕一区二区三区有码在线看| 欧美日韩国产亚洲二区| 熟妇人妻久久中文字幕3abv| 中文资源天堂在线| 乱码一卡2卡4卡精品| 中文在线观看免费www的网站| 在线播放国产精品三级| 日本撒尿小便嘘嘘汇集6| 精品欧美国产一区二区三| 国产av不卡久久| 神马国产精品三级电影在线观看| 老司机深夜福利视频在线观看| 在线观看舔阴道视频| 亚洲中文日韩欧美视频| 老司机午夜福利在线观看视频| 99热这里只有精品一区| 高潮久久久久久久久久久不卡| 亚洲成人久久爱视频| 99久久无色码亚洲精品果冻| 他把我摸到了高潮在线观看| 天堂动漫精品| 免费大片18禁| 高清毛片免费观看视频网站| 国产精品美女特级片免费视频播放器| 又爽又黄无遮挡网站| 国产在视频线在精品| 午夜视频国产福利| 午夜视频国产福利| 三级男女做爰猛烈吃奶摸视频| 一进一出抽搐gif免费好疼| 97碰自拍视频| 一进一出好大好爽视频| 欧美日韩中文字幕国产精品一区二区三区| 国产黄色小视频在线观看| 一本一本综合久久| 国产单亲对白刺激| 天堂网av新在线| 精品人妻一区二区三区麻豆 | 亚洲色图av天堂| 91久久精品电影网| 丁香欧美五月| 99热这里只有精品一区| 一本一本综合久久| 嫁个100分男人电影在线观看| 看片在线看免费视频| 成年女人毛片免费观看观看9| 桃色一区二区三区在线观看| 一边摸一边抽搐一进一小说| 久久天躁狠狠躁夜夜2o2o| 国产真实乱freesex| 特大巨黑吊av在线直播| 亚洲人成网站高清观看| 国产免费男女视频| 亚洲av二区三区四区| 久久久久精品国产欧美久久久| 黄片小视频在线播放| 精品久久久久久久久亚洲 | 一边摸一边抽搐一进一小说| 我的女老师完整版在线观看| 亚洲乱码一区二区免费版| 亚洲精品一区av在线观看| 黄色女人牲交| 国产三级黄色录像| www.www免费av| 两人在一起打扑克的视频| 男人的好看免费观看在线视频| 亚洲成人中文字幕在线播放| 3wmmmm亚洲av在线观看| 欧美+亚洲+日韩+国产| 搡老岳熟女国产| 欧美黄色片欧美黄色片| 岛国在线免费视频观看| 在线免费观看不下载黄p国产 | 麻豆成人午夜福利视频| 欧美3d第一页| 免费观看的影片在线观看| 午夜老司机福利剧场| 久久精品国产清高在天天线| 亚洲美女视频黄频| 午夜福利免费观看在线| 亚洲aⅴ乱码一区二区在线播放| 一个人观看的视频www高清免费观看| 久久人人爽人人爽人人片va | 国产欧美日韩一区二区精品| 性色av乱码一区二区三区2| 欧美性感艳星| 毛片女人毛片| 毛片女人毛片| bbb黄色大片| 别揉我奶头 嗯啊视频| 女人被狂操c到高潮| 狂野欧美白嫩少妇大欣赏| 国产精品自产拍在线观看55亚洲| 午夜福利欧美成人| 国产伦精品一区二区三区视频9| 人人妻人人澡欧美一区二区| 少妇被粗大猛烈的视频| www.999成人在线观看| 久久天躁狠狠躁夜夜2o2o| 精品人妻一区二区三区麻豆 | 国产亚洲精品久久久久久毛片| 国产乱人视频| www.色视频.com| 国内揄拍国产精品人妻在线| 国产精品影院久久| 婷婷精品国产亚洲av在线| ponron亚洲| 波多野结衣高清无吗| 国产伦一二天堂av在线观看| 1000部很黄的大片| 国产精品影院久久| 欧美三级亚洲精品| 欧美潮喷喷水| 又黄又爽又免费观看的视频| 亚洲一区二区三区不卡视频| 少妇高潮的动态图| 色综合婷婷激情| 国产大屁股一区二区在线视频| 可以在线观看毛片的网站| 51国产日韩欧美| 亚洲五月婷婷丁香| 一边摸一边抽搐一进一小说| bbb黄色大片| 亚洲黑人精品在线| 亚洲精华国产精华精| 黄色视频,在线免费观看| 国产精品日韩av在线免费观看| 91麻豆精品激情在线观看国产| 精品免费久久久久久久清纯| 午夜激情福利司机影院| 精品久久久久久久久久久久久| 精品久久久久久久久久久久久| 国产精品乱码一区二三区的特点| 校园春色视频在线观看| 亚洲精品影视一区二区三区av| 久久精品影院6| 欧美成人a在线观看| 日本撒尿小便嘘嘘汇集6| 国产欧美日韩精品一区二区| 波多野结衣巨乳人妻| 日韩亚洲欧美综合| 在线看三级毛片| 他把我摸到了高潮在线观看| 精品国产三级普通话版| 他把我摸到了高潮在线观看| 最近视频中文字幕2019在线8| 最近最新中文字幕大全电影3| 国产精品亚洲av一区麻豆| 久久久久精品国产欧美久久久| 国产精品1区2区在线观看.| 亚洲av成人不卡在线观看播放网| 亚洲美女黄片视频| 精品国产三级普通话版| 国产av不卡久久| 国产黄色小视频在线观看| 禁无遮挡网站| 国产亚洲欧美98| 99久久精品一区二区三区| 网址你懂的国产日韩在线| 亚洲人成网站在线播放欧美日韩| 此物有八面人人有两片| 波多野结衣高清作品| 毛片一级片免费看久久久久 | 国产乱人伦免费视频| 特级一级黄色大片| 午夜两性在线视频| 久久精品国产亚洲av香蕉五月| 波多野结衣高清作品| 日本一二三区视频观看| 亚洲av电影不卡..在线观看| 久久久久免费精品人妻一区二区| 午夜福利高清视频| 十八禁人妻一区二区| 久久6这里有精品| 久久6这里有精品| 成人性生交大片免费视频hd| aaaaa片日本免费| 看片在线看免费视频| 1000部很黄的大片| 国产高清激情床上av| 最近中文字幕高清免费大全6 | 欧美一区二区国产精品久久精品| av女优亚洲男人天堂| 精品久久久久久久久久免费视频| 亚洲精品一卡2卡三卡4卡5卡| 国产伦精品一区二区三区视频9| 国模一区二区三区四区视频| 免费黄网站久久成人精品 | АⅤ资源中文在线天堂| 桃红色精品国产亚洲av| 桃色一区二区三区在线观看| 精品人妻一区二区三区麻豆 | 亚洲av一区综合| 在线免费观看不下载黄p国产 | 国产成+人综合+亚洲专区| 自拍偷自拍亚洲精品老妇| 久久久精品欧美日韩精品| aaaaa片日本免费| 搡老熟女国产l中国老女人| 午夜福利18| 国产黄色小视频在线观看| 蜜桃亚洲精品一区二区三区| 欧美+日韩+精品| 亚洲精品日韩av片在线观看| 日本免费a在线| 男人狂女人下面高潮的视频| 看十八女毛片水多多多| 免费看美女性在线毛片视频| 日日夜夜操网爽| 少妇被粗大猛烈的视频| 国产成人aa在线观看| 欧美最新免费一区二区三区 | 老女人水多毛片| 村上凉子中文字幕在线| 51国产日韩欧美| 日本成人三级电影网站| 女人十人毛片免费观看3o分钟| 美女被艹到高潮喷水动态| 草草在线视频免费看| 免费在线观看成人毛片| 欧美三级亚洲精品| 精品久久久久久久末码| 国产精品不卡视频一区二区 | 99久久99久久久精品蜜桃| 国产91精品成人一区二区三区| 日本撒尿小便嘘嘘汇集6| 亚洲 欧美 日韩 在线 免费| 色播亚洲综合网| 国产精品人妻久久久久久| 国产高清激情床上av| 一级黄片播放器| 三级男女做爰猛烈吃奶摸视频| 黄色视频,在线免费观看| 亚洲人成网站高清观看| 国产精品久久久久久精品电影| 日本撒尿小便嘘嘘汇集6| 成年女人看的毛片在线观看| 国模一区二区三区四区视频| 在线十欧美十亚洲十日本专区| 哪里可以看免费的av片| 亚洲欧美精品综合久久99| 99国产精品一区二区三区| 亚洲成人中文字幕在线播放| 亚洲五月天丁香| 亚洲av成人精品一区久久| 中文字幕av成人在线电影| 亚洲av成人av| 亚洲精品影视一区二区三区av| 亚洲久久久久久中文字幕| 动漫黄色视频在线观看| 国产成年人精品一区二区| 好男人在线观看高清免费视频| 欧美一级a爱片免费观看看| 亚洲专区中文字幕在线| 欧美黄色片欧美黄色片| 1024手机看黄色片| 久久久久久久亚洲中文字幕 | 少妇裸体淫交视频免费看高清| 最近在线观看免费完整版| 亚洲在线自拍视频| 在线a可以看的网站| 大型黄色视频在线免费观看| 国产探花极品一区二区| 亚洲成人久久性| 国产精品,欧美在线| 午夜老司机福利剧场| 国产在视频线在精品| 一区二区三区激情视频| 亚洲熟妇中文字幕五十中出| 少妇人妻精品综合一区二区 | 亚洲男人的天堂狠狠| 在线观看免费视频日本深夜| 欧美精品国产亚洲| 成人性生交大片免费视频hd| 亚洲无线观看免费| 精品久久久久久,| 久久久久久久精品吃奶| 国产高清三级在线| 欧美+亚洲+日韩+国产| 一区二区三区四区激情视频 | 国产乱人伦免费视频| 成年女人永久免费观看视频| 最新在线观看一区二区三区| 亚洲av五月六月丁香网| 熟妇人妻久久中文字幕3abv| 午夜影院日韩av| 搡老妇女老女人老熟妇| 欧美国产日韩亚洲一区| 国产免费一级a男人的天堂| 一二三四社区在线视频社区8| 18禁黄网站禁片午夜丰满| 波野结衣二区三区在线| 国产精品日韩av在线免费观看| 免费黄网站久久成人精品 | 99riav亚洲国产免费| 韩国av一区二区三区四区| 日韩欧美一区二区三区在线观看| 99热这里只有是精品在线观看 | 成人永久免费在线观看视频| 亚洲av.av天堂| 欧美丝袜亚洲另类 | 波野结衣二区三区在线| 国产免费av片在线观看野外av| 欧美色欧美亚洲另类二区| 婷婷精品国产亚洲av| www.色视频.com| 国产人妻一区二区三区在| 国产探花在线观看一区二区| 一区二区三区高清视频在线| 亚洲av成人av| 久久精品国产99精品国产亚洲性色| 在线观看66精品国产| 久久久久久国产a免费观看| 最新在线观看一区二区三区| av福利片在线观看| 亚洲国产欧美人成| 亚洲av不卡在线观看| 国产老妇女一区| 成人特级av手机在线观看| 成熟少妇高潮喷水视频| 亚洲 欧美 日韩 在线 免费| 亚洲性夜色夜夜综合| 丰满人妻熟妇乱又伦精品不卡| 999久久久精品免费观看国产| a级毛片a级免费在线| 亚洲欧美精品综合久久99| 在线天堂最新版资源| 免费看美女性在线毛片视频| 9191精品国产免费久久| 久久精品综合一区二区三区| 国产精品一区二区免费欧美| 国产蜜桃级精品一区二区三区| 成年免费大片在线观看| 神马国产精品三级电影在线观看| 99在线视频只有这里精品首页| 日本成人三级电影网站| 天天躁日日操中文字幕| 少妇高潮的动态图| 岛国在线免费视频观看| 日本熟妇午夜| 欧美一区二区亚洲| 免费av观看视频| 亚洲无线在线观看| 波多野结衣巨乳人妻| 亚洲18禁久久av| 天堂影院成人在线观看| 久久精品国产亚洲av香蕉五月| 欧美成人免费av一区二区三区| 亚洲精品一区av在线观看| 日韩有码中文字幕| av在线蜜桃| 男人的好看免费观看在线视频| h日本视频在线播放| 日本a在线网址| 午夜视频国产福利| 国产一区二区亚洲精品在线观看| 可以在线观看毛片的网站| 亚洲美女视频黄频| 免费黄网站久久成人精品 | 热99re8久久精品国产| 日本免费一区二区三区高清不卡| 亚洲人成电影免费在线| 亚洲精品在线美女| 禁无遮挡网站| 三级国产精品欧美在线观看| 久久亚洲真实| av福利片在线观看| 欧美一区二区亚洲| 脱女人内裤的视频| 麻豆国产97在线/欧美| 中文字幕免费在线视频6| 亚洲成人中文字幕在线播放| 免费在线观看亚洲国产| 精品人妻1区二区| 国产成人aa在线观看| 国产一区二区三区视频了| 中国美女看黄片| 日本成人三级电影网站| 成人午夜高清在线视频| 国产免费男女视频| 性色av乱码一区二区三区2| 18+在线观看网站| 亚洲国产欧美人成| 真人做人爱边吃奶动态| 久久久久久久午夜电影| 色综合婷婷激情| 色综合站精品国产| 看片在线看免费视频| 成人特级av手机在线观看| 亚洲中文日韩欧美视频| 偷拍熟女少妇极品色| 搞女人的毛片| 久久久久久久精品吃奶| h日本视频在线播放| 色精品久久人妻99蜜桃| 欧美区成人在线视频| 天天躁日日操中文字幕| 免费av不卡在线播放| 日韩欧美精品免费久久 | 91狼人影院| 一边摸一边抽搐一进一小说| 免费人成在线观看视频色| 婷婷精品国产亚洲av在线| 久久天躁狠狠躁夜夜2o2o| 久久久久久久久中文| 日韩欧美三级三区| 精品一区二区三区av网在线观看| 久久国产乱子免费精品| 别揉我奶头~嗯~啊~动态视频| 久久久久国产精品人妻aⅴ院| 国产欧美日韩精品亚洲av| 亚洲精品456在线播放app | 好男人在线观看高清免费视频| 国产熟女xx| 久久精品久久久久久噜噜老黄 | 美女cb高潮喷水在线观看| 国产精品人妻久久久久久| 精品人妻视频免费看| 热99re8久久精品国产| 国产精品99久久久久久久久| 亚洲自偷自拍三级| 欧美另类亚洲清纯唯美| 欧美+亚洲+日韩+国产| a在线观看视频网站| 成人永久免费在线观看视频| 国产黄a三级三级三级人| 天堂影院成人在线观看| 亚洲乱码一区二区免费版| 精品午夜福利视频在线观看一区| 99久久久亚洲精品蜜臀av| 亚洲精华国产精华精| 一区二区三区激情视频| 非洲黑人性xxxx精品又粗又长| 成人三级黄色视频| 直男gayav资源| 长腿黑丝高跟| 成人永久免费在线观看视频| 欧美成人一区二区免费高清观看| 午夜福利高清视频| 亚洲欧美精品综合久久99| 91麻豆精品激情在线观看国产| 欧美最黄视频在线播放免费| 偷拍熟女少妇极品色| 精品福利观看| 久9热在线精品视频| 久久久久亚洲av毛片大全| 亚洲第一欧美日韩一区二区三区| 国产精品,欧美在线| 天堂av国产一区二区熟女人妻| 听说在线观看完整版免费高清| 成人高潮视频无遮挡免费网站| 最近最新免费中文字幕在线| 久久久久久九九精品二区国产| 国产视频内射| 色综合亚洲欧美另类图片| 可以在线观看毛片的网站| 日本免费a在线| 国产老妇女一区| 亚洲成人久久爱视频| 久久国产乱子免费精品| 男女那种视频在线观看| 99国产综合亚洲精品| 91在线观看av| 极品教师在线免费播放| 国产成人a区在线观看| 一二三四社区在线视频社区8| 精品99又大又爽又粗少妇毛片 | 99久久久亚洲精品蜜臀av| 亚洲国产精品成人综合色| 国产精品亚洲美女久久久| 麻豆国产97在线/欧美| 怎么达到女性高潮| 国产精品免费一区二区三区在线| 成年女人毛片免费观看观看9| 91午夜精品亚洲一区二区三区 | 一边摸一边抽搐一进一小说| 日韩精品中文字幕看吧| 久久久久久国产a免费观看| 亚洲精品456在线播放app | 此物有八面人人有两片| 国产精华一区二区三区| 欧美日韩国产亚洲二区| 久久精品国产亚洲av香蕉五月| 久久久久性生活片| 九九在线视频观看精品| 亚洲18禁久久av| 亚洲av电影不卡..在线观看| 国产三级中文精品| 久99久视频精品免费| 99视频精品全部免费 在线| 毛片一级片免费看久久久久 | 欧美日韩国产亚洲二区| 性色avwww在线观看| 波多野结衣高清无吗| 欧美极品一区二区三区四区| 精品午夜福利在线看| 亚洲一区二区三区不卡视频| 一边摸一边抽搐一进一小说| 午夜福利在线在线| av黄色大香蕉| 淫秽高清视频在线观看| 欧美最新免费一区二区三区 | 久久香蕉精品热| 免费在线观看成人毛片| 亚洲激情在线av| 人人妻人人看人人澡| 国产精品久久久久久人妻精品电影| 久久香蕉精品热| 久久99热6这里只有精品| 久久精品国产清高在天天线| 色综合婷婷激情| 午夜福利在线观看吧| 国产爱豆传媒在线观看| 亚洲美女搞黄在线观看 | 国产主播在线观看一区二区| 国产伦精品一区二区三区四那| 免费人成视频x8x8入口观看| 三级国产精品欧美在线观看| 欧美成人免费av一区二区三区| АⅤ资源中文在线天堂| 欧美乱色亚洲激情| 在现免费观看毛片| 精品久久久久久久久亚洲 | 不卡一级毛片| 亚洲欧美精品综合久久99| 哪里可以看免费的av片| 又黄又爽又免费观看的视频| 国产成人aa在线观看| 热99在线观看视频| 日韩欧美在线乱码| a在线观看视频网站| 午夜影院日韩av| 精品一区二区三区av网在线观看| 日本a在线网址| 日韩欧美 国产精品| 国产精品久久久久久久久免 | 国产单亲对白刺激| 88av欧美| 欧美激情久久久久久爽电影| 亚洲国产色片| 国产真实伦视频高清在线观看 | 欧美日韩黄片免| 成熟少妇高潮喷水视频| 国产精品日韩av在线免费观看| 中文资源天堂在线| 制服丝袜大香蕉在线| 婷婷丁香在线五月| 国产亚洲精品av在线| АⅤ资源中文在线天堂| 日本一本二区三区精品| 亚洲综合色惰| 久久人妻av系列| 欧美黑人欧美精品刺激| 丁香六月欧美| 综合色av麻豆| 欧美日韩瑟瑟在线播放| 波野结衣二区三区在线| 男女之事视频高清在线观看| 亚洲,欧美,日韩| 岛国在线免费视频观看| 国语自产精品视频在线第100页| 亚州av有码| 国产综合懂色| АⅤ资源中文在线天堂| 91在线观看av| 岛国在线免费视频观看| 久久香蕉精品热| 成人av在线播放网站| 午夜福利在线在线| 欧美另类亚洲清纯唯美| 精品乱码久久久久久99久播| 精品久久久久久久人妻蜜臀av| 国产精品女同一区二区软件 | 午夜视频国产福利|