• <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ī)制分析
    我要看日韩黄色一级片| 日韩欧美 国产精品| 国模一区二区三区四区视频| 看非洲黑人一级黄片| 美女脱内裤让男人舔精品视频| 免费看日本二区| av国产精品久久久久影院| www.色视频.com| h日本视频在线播放| 国产在线一区二区三区精| 国产在视频线精品| 少妇被粗大猛烈的视频| av专区在线播放| 18禁在线无遮挡免费观看视频| 国产亚洲91精品色在线| 免费高清在线观看视频在线观看| 免费看不卡的av| 简卡轻食公司| 一级毛片电影观看| 免费观看a级毛片全部| 亚洲,欧美,日韩| 国产69精品久久久久777片| 欧美zozozo另类| 精品午夜福利在线看| 日本与韩国留学比较| 亚洲久久久国产精品| 亚洲精品成人av观看孕妇| 久久久久人妻精品一区果冻| 九九久久精品国产亚洲av麻豆| 亚洲国产欧美人成| 国产v大片淫在线免费观看| av不卡在线播放| 国产高清有码在线观看视频| 国产成人91sexporn| 国产乱人视频| 看免费成人av毛片| av专区在线播放| 赤兔流量卡办理| 免费黄网站久久成人精品| 卡戴珊不雅视频在线播放| 亚洲伊人久久精品综合| 午夜福利在线观看免费完整高清在| 高清视频免费观看一区二区| 99热这里只有是精品50| 亚洲精品第二区| 国产成人a∨麻豆精品| 男女国产视频网站| 观看av在线不卡| 边亲边吃奶的免费视频| 看免费成人av毛片| 久热这里只有精品99| 最近最新中文字幕大全电影3| 成人二区视频| 国产精品久久久久久精品电影小说 | 夜夜爽夜夜爽视频| 亚洲成人一二三区av| 日本欧美国产在线视频| 在线观看国产h片| 十分钟在线观看高清视频www | 一级黄片播放器| av国产精品久久久久影院| 99九九线精品视频在线观看视频| 777米奇影视久久| 日本色播在线视频| 男的添女的下面高潮视频| 少妇高潮的动态图| 中文资源天堂在线| tube8黄色片| 欧美极品一区二区三区四区| 水蜜桃什么品种好| 久久国内精品自在自线图片| 亚洲欧美日韩卡通动漫| 国产男人的电影天堂91| 少妇被粗大猛烈的视频| 街头女战士在线观看网站| 日韩av不卡免费在线播放| 全区人妻精品视频| 亚洲国产日韩一区二区| 国产精品人妻久久久久久| 夜夜爽夜夜爽视频| 国产大屁股一区二区在线视频| 亚洲图色成人| 美女高潮的动态| 免费在线观看成人毛片| 亚洲国产精品成人久久小说| 一本一本综合久久| 国产精品爽爽va在线观看网站| 18禁裸乳无遮挡动漫免费视频| 草草在线视频免费看| 欧美日韩视频精品一区| 永久免费av网站大全| 网址你懂的国产日韩在线| 亚洲av中文字字幕乱码综合| 老师上课跳d突然被开到最大视频| 建设人人有责人人尽责人人享有的 | 成人黄色视频免费在线看| .国产精品久久| 国产毛片在线视频| 精品久久久噜噜| 久久久久久久国产电影| 美女高潮的动态| 国产人妻一区二区三区在| 18禁在线无遮挡免费观看视频| 亚洲aⅴ乱码一区二区在线播放| 久久ye,这里只有精品| 又粗又硬又长又爽又黄的视频| 精品久久久久久电影网| av免费观看日本| 国产精品国产三级专区第一集| 国产黄色免费在线视频| 亚洲欧美精品自产自拍| 波野结衣二区三区在线| 免费看光身美女| 中文资源天堂在线| 联通29元200g的流量卡| 国产一级毛片在线| 国产亚洲一区二区精品| 日本猛色少妇xxxxx猛交久久| 最近最新中文字幕免费大全7| 99久久中文字幕三级久久日本| 国产成人免费观看mmmm| 一本色道久久久久久精品综合| 亚洲国产欧美在线一区| 高清av免费在线| 伦理电影免费视频| 99久国产av精品国产电影| 久久午夜福利片| 欧美xxxx黑人xx丫x性爽| 欧美日韩国产mv在线观看视频 | 97超碰精品成人国产| 亚洲无线观看免费| 国产伦在线观看视频一区| 亚洲无线观看免费| 亚洲中文av在线| 黄色一级大片看看| 一个人看的www免费观看视频| 成人亚洲欧美一区二区av| 国产av码专区亚洲av| 全区人妻精品视频| 永久免费av网站大全| 国产乱人视频| 深夜a级毛片| 国产精品秋霞免费鲁丝片| 亚洲自偷自拍三级| 色哟哟·www| 在线观看美女被高潮喷水网站| 精品少妇久久久久久888优播| av在线播放精品| 久久韩国三级中文字幕| 午夜老司机福利剧场| 亚洲欧美中文字幕日韩二区| 青春草国产在线视频| 色视频www国产| 国产精品蜜桃在线观看| 亚洲不卡免费看| av在线老鸭窝| 久久久精品免费免费高清| 内射极品少妇av片p| 97超碰精品成人国产| 十八禁网站网址无遮挡 | 中文字幕久久专区| 久久久精品免费免费高清| 91精品一卡2卡3卡4卡| 久久久久国产网址| 人体艺术视频欧美日本| 午夜福利在线观看免费完整高清在| 我的老师免费观看完整版| 国产亚洲5aaaaa淫片| 在线观看一区二区三区| 亚洲av.av天堂| 国产亚洲91精品色在线| 国产亚洲欧美精品永久| 一级二级三级毛片免费看| 国精品久久久久久国模美| 好男人视频免费观看在线| 高清不卡的av网站| 日韩免费高清中文字幕av| 国产精品嫩草影院av在线观看| 亚洲精品日本国产第一区| 国产老妇伦熟女老妇高清| 国产69精品久久久久777片| 丝瓜视频免费看黄片| 涩涩av久久男人的天堂| 少妇的逼水好多| 亚洲性久久影院| av福利片在线观看| 全区人妻精品视频| 男女啪啪激烈高潮av片| av播播在线观看一区| 亚洲美女视频黄频| 国产精品一及| 国产国拍精品亚洲av在线观看| av一本久久久久| 最后的刺客免费高清国语| 亚洲av免费高清在线观看| 婷婷色av中文字幕| 又爽又黄a免费视频| 赤兔流量卡办理| 国产精品一区二区在线观看99| 我的老师免费观看完整版| 国产极品天堂在线| 色婷婷av一区二区三区视频| 2018国产大陆天天弄谢| 高清不卡的av网站| 亚洲人与动物交配视频| 国产久久久一区二区三区| 亚洲av成人精品一区久久| 国产欧美日韩精品一区二区| 久久亚洲国产成人精品v| 日本黄色日本黄色录像| 国产免费一级a男人的天堂| 亚洲一区二区三区欧美精品| 久久久成人免费电影| 男女边吃奶边做爰视频| 亚洲aⅴ乱码一区二区在线播放| 色婷婷久久久亚洲欧美| 免费看av在线观看网站| 乱码一卡2卡4卡精品| 伦精品一区二区三区| 一级毛片黄色毛片免费观看视频| av线在线观看网站| 18+在线观看网站| 欧美精品人与动牲交sv欧美| 中文欧美无线码| 国产精品一区二区在线不卡| 日韩 亚洲 欧美在线| 老女人水多毛片| 久久久久久久久久人人人人人人| 深夜a级毛片| 18禁在线播放成人免费| 少妇熟女欧美另类| 国产欧美日韩精品一区二区| 日本黄大片高清| 黑人高潮一二区| 夜夜骑夜夜射夜夜干| 亚洲图色成人| 亚洲电影在线观看av| 大片电影免费在线观看免费| 日日啪夜夜爽| 久久国产乱子免费精品| a级毛色黄片| 亚洲欧美一区二区三区黑人 | 一级av片app| 纵有疾风起免费观看全集完整版| 美女cb高潮喷水在线观看| 亚洲欧美日韩东京热| 男女边吃奶边做爰视频| 中文字幕久久专区| 亚洲精品乱久久久久久| 极品少妇高潮喷水抽搐| 日韩不卡一区二区三区视频在线| 香蕉精品网在线| 中文字幕av成人在线电影| 一本久久精品| 成人高潮视频无遮挡免费网站| 一边亲一边摸免费视频| 亚洲人成网站在线播| 日产精品乱码卡一卡2卡三| 午夜福利影视在线免费观看| av天堂中文字幕网| 男男h啪啪无遮挡| 尾随美女入室| 夫妻午夜视频| 国产精品麻豆人妻色哟哟久久| 看非洲黑人一级黄片| 久久久久性生活片| 午夜福利高清视频| 亚洲国产色片| 国产免费一级a男人的天堂| 天天躁日日操中文字幕| 欧美激情国产日韩精品一区| 一级二级三级毛片免费看| 男人添女人高潮全过程视频| 成人国产av品久久久| 国产成人免费无遮挡视频| 欧美高清成人免费视频www| 国产精品国产三级国产av玫瑰| 中文字幕制服av| 国产亚洲5aaaaa淫片| 国产片特级美女逼逼视频| 黑丝袜美女国产一区| av国产精品久久久久影院| 国产精品人妻久久久久久| 国产亚洲91精品色在线| 免费少妇av软件| 精品国产三级普通话版| 欧美成人精品欧美一级黄| 91午夜精品亚洲一区二区三区| 午夜福利影视在线免费观看| 大陆偷拍与自拍| 高清不卡的av网站| 欧美激情极品国产一区二区三区 | 看十八女毛片水多多多| 国产黄片视频在线免费观看| 人妻一区二区av| 国产男女内射视频| 成年女人在线观看亚洲视频| 大片电影免费在线观看免费| 男人添女人高潮全过程视频| 国产成人精品婷婷| 免费黄频网站在线观看国产| 美女国产视频在线观看| av在线app专区| av国产免费在线观看| 99热这里只有是精品在线观看| 人妻夜夜爽99麻豆av| 国产在视频线精品| 国产高清不卡午夜福利| av国产久精品久网站免费入址| 国产伦在线观看视频一区| 欧美另类一区| 国产又色又爽无遮挡免| 搡老乐熟女国产| 国产伦精品一区二区三区四那| 五月玫瑰六月丁香| 日韩视频在线欧美| 国产亚洲欧美精品永久| 免费人成在线观看视频色| 日韩中文字幕视频在线看片 | 有码 亚洲区| 国产成人精品一,二区| 亚洲av.av天堂| av线在线观看网站| 色婷婷av一区二区三区视频| 国产日韩欧美亚洲二区| 亚洲av成人精品一二三区| 18禁动态无遮挡网站| 国产精品无大码| a级毛色黄片| 亚洲精品日韩在线中文字幕| 99热6这里只有精品| 交换朋友夫妻互换小说| 国产成人aa在线观看| 亚洲精品aⅴ在线观看| av不卡在线播放| 插逼视频在线观看| 99久久精品国产国产毛片| 九九在线视频观看精品| 国语对白做爰xxxⅹ性视频网站| 久久女婷五月综合色啪小说| 亚洲高清免费不卡视频| 夫妻性生交免费视频一级片| 日韩欧美精品免费久久| 久久人妻熟女aⅴ| 夜夜骑夜夜射夜夜干| 九九在线视频观看精品| 亚洲四区av| 日本-黄色视频高清免费观看| 日韩免费高清中文字幕av| 最近中文字幕2019免费版| 少妇猛男粗大的猛烈进出视频| 人妻制服诱惑在线中文字幕| 国产av国产精品国产| 伦精品一区二区三区| 偷拍熟女少妇极品色| 中国美白少妇内射xxxbb| 狂野欧美激情性xxxx在线观看| 赤兔流量卡办理| 一级毛片电影观看| 另类亚洲欧美激情| 国产精品一区二区在线不卡| 欧美激情国产日韩精品一区| 亚洲欧美日韩卡通动漫| av福利片在线观看| 久久国内精品自在自线图片| 成人毛片60女人毛片免费| 久久久久网色| 亚洲av中文av极速乱| 少妇人妻精品综合一区二区| 亚洲欧美精品专区久久| 天堂俺去俺来也www色官网| 国产精品偷伦视频观看了| 国产精品一二三区在线看| 2022亚洲国产成人精品| 男人狂女人下面高潮的视频| 美女主播在线视频| 亚洲精品乱码久久久v下载方式| 亚洲美女视频黄频| 精品亚洲成国产av| 一级片'在线观看视频| 美女中出高潮动态图| 自拍欧美九色日韩亚洲蝌蚪91 | 男的添女的下面高潮视频| 熟女人妻精品中文字幕| 免费人妻精品一区二区三区视频| 久久精品久久久久久久性| 91aial.com中文字幕在线观看| 成年免费大片在线观看| 久久久久久久久久久丰满| 九九久久精品国产亚洲av麻豆| 高清在线视频一区二区三区| 亚洲国产欧美在线一区| 欧美变态另类bdsm刘玥| 天堂8中文在线网| 免费看光身美女| 一二三四中文在线观看免费高清| 欧美区成人在线视频| av不卡在线播放| 日韩欧美 国产精品| 狂野欧美激情性bbbbbb| 日韩一本色道免费dvd| 少妇人妻精品综合一区二区| 精品人妻视频免费看| 国产欧美日韩一区二区三区在线 | 国产 一区精品| 一级av片app| 男女下面进入的视频免费午夜| 久久精品国产a三级三级三级| 久久久久久伊人网av| 在线观看一区二区三区激情| 欧美少妇被猛烈插入视频| 欧美老熟妇乱子伦牲交| 国产高清有码在线观看视频| 大陆偷拍与自拍| 亚洲av免费高清在线观看| 男的添女的下面高潮视频| 韩国av在线不卡| 少妇人妻一区二区三区视频| 亚洲精品一区蜜桃| 一级片'在线观看视频| 欧美变态另类bdsm刘玥| 久久青草综合色| 伦理电影大哥的女人| 久久午夜福利片| 成年美女黄网站色视频大全免费 | 婷婷色综合www| 永久免费av网站大全| 亚洲欧美日韩另类电影网站 | 最近的中文字幕免费完整| 一区二区三区精品91| 极品教师在线视频| 99久久中文字幕三级久久日本| 中文字幕久久专区| 偷拍熟女少妇极品色| 大香蕉久久网| 日韩免费高清中文字幕av| 人妻夜夜爽99麻豆av| 少妇熟女欧美另类| 最新中文字幕久久久久| 免费黄频网站在线观看国产| 激情五月婷婷亚洲| 精品少妇黑人巨大在线播放| 欧美日韩国产mv在线观看视频 | 一级av片app| 国产精品秋霞免费鲁丝片| 久久久久久久亚洲中文字幕| 搡老乐熟女国产| 九草在线视频观看| av不卡在线播放| 男女啪啪激烈高潮av片| 欧美 日韩 精品 国产| 美女国产视频在线观看| 99久久精品一区二区三区| 久久人人爽av亚洲精品天堂 | 黄色视频在线播放观看不卡| 国产黄频视频在线观看| 亚洲人成网站在线播| 97超视频在线观看视频| 亚洲精华国产精华液的使用体验| av国产精品久久久久影院| 99视频精品全部免费 在线| 简卡轻食公司| 哪个播放器可以免费观看大片| 国产乱人偷精品视频| 高清av免费在线| 人人妻人人爽人人添夜夜欢视频 | 欧美xxxx性猛交bbbb| 国产精品福利在线免费观看| 久久精品国产鲁丝片午夜精品| 国产成人精品婷婷| 精品少妇黑人巨大在线播放| 亚洲精品国产av成人精品| 身体一侧抽搐| 婷婷色av中文字幕| 国产精品国产三级国产av玫瑰| 人妻 亚洲 视频| 久久精品久久精品一区二区三区| 人体艺术视频欧美日本| 国产中年淑女户外野战色| 99热国产这里只有精品6| 国产精品麻豆人妻色哟哟久久| 伦理电影免费视频| 网址你懂的国产日韩在线| 国产老妇伦熟女老妇高清| 老熟女久久久| 亚州av有码| 蜜桃在线观看..| 国产成人a∨麻豆精品| 1000部很黄的大片| 国产精品久久久久久久久免| 亚洲自偷自拍三级| 亚洲,欧美,日韩| 精品久久久精品久久久| 久久久午夜欧美精品| 欧美高清性xxxxhd video| 午夜激情久久久久久久| 欧美人与善性xxx| 午夜老司机福利剧场| 男女边摸边吃奶| 精品99又大又爽又粗少妇毛片| 免费不卡的大黄色大毛片视频在线观看| 免费看不卡的av| 欧美成人a在线观看| av网站免费在线观看视频| 晚上一个人看的免费电影| av又黄又爽大尺度在线免费看| 人人妻人人看人人澡| 国产精品久久久久久av不卡| 天天躁夜夜躁狠狠久久av| 中文字幕免费在线视频6| 欧美精品国产亚洲| 日本黄色日本黄色录像| 色婷婷久久久亚洲欧美| 哪个播放器可以免费观看大片| 少妇精品久久久久久久| 天堂8中文在线网| 亚洲国产高清在线一区二区三| 麻豆精品久久久久久蜜桃| 亚洲综合色惰| 日日摸夜夜添夜夜添av毛片| 午夜福利影视在线免费观看| 在线观看免费日韩欧美大片 | 久久久久久人妻| 亚洲内射少妇av| 九草在线视频观看| 国产 精品1| 免费播放大片免费观看视频在线观看| 精品一区二区三卡| 亚洲精品久久午夜乱码| 国产免费一级a男人的天堂| 晚上一个人看的免费电影| 午夜激情久久久久久久| 欧美激情极品国产一区二区三区 | 婷婷色综合大香蕉| 成人一区二区视频在线观看| 久久国产亚洲av麻豆专区| 国内少妇人妻偷人精品xxx网站| 在线亚洲精品国产二区图片欧美 | 一本色道久久久久久精品综合| 尾随美女入室| 久久国内精品自在自线图片| 亚洲天堂av无毛| 中文字幕久久专区| 国产免费视频播放在线视频| 91午夜精品亚洲一区二区三区| 国产久久久一区二区三区| av一本久久久久| 欧美日韩综合久久久久久| 国产视频内射| 亚洲自偷自拍三级| 黄色配什么色好看| 国产美女午夜福利| 国产成人91sexporn| 尾随美女入室| 国产精品熟女久久久久浪| av国产免费在线观看| h视频一区二区三区| 亚洲aⅴ乱码一区二区在线播放| 国产永久视频网站| 九九爱精品视频在线观看| 国产视频首页在线观看| 国产男女超爽视频在线观看| 99热这里只有是精品在线观看| 亚洲精品乱久久久久久| 亚洲最大成人中文| 亚洲欧美一区二区三区国产| 日本wwww免费看| 日日摸夜夜添夜夜添av毛片| 国内精品宾馆在线| 久久精品国产自在天天线| 一级毛片aaaaaa免费看小| 一二三四中文在线观看免费高清| 亚洲精品视频女| 免费黄频网站在线观看国产| 免费观看性生交大片5| 22中文网久久字幕| 极品少妇高潮喷水抽搐| 日韩成人av中文字幕在线观看| 国产女主播在线喷水免费视频网站| 18禁动态无遮挡网站| 一级毛片电影观看| 高清午夜精品一区二区三区| 极品教师在线视频| a级毛片免费高清观看在线播放| 精品国产露脸久久av麻豆| 亚洲精品456在线播放app| 国产成人精品婷婷| 波野结衣二区三区在线| 欧美丝袜亚洲另类| 欧美精品国产亚洲| 亚洲经典国产精华液单| 日本黄色片子视频| 日韩av免费高清视频| 国产白丝娇喘喷水9色精品| 免费黄色在线免费观看| 亚州av有码| 女性生殖器流出的白浆| 国产视频内射| 一边亲一边摸免费视频| 国产高清不卡午夜福利| 蜜桃亚洲精品一区二区三区| 久久久精品94久久精品| 九九在线视频观看精品| 女的被弄到高潮叫床怎么办| 国产精品爽爽va在线观看网站| 精品熟女少妇av免费看| 熟妇人妻不卡中文字幕| 成人18禁高潮啪啪吃奶动态图 | 小蜜桃在线观看免费完整版高清| 亚洲精品国产av成人精品| 国产免费福利视频在线观看| 欧美日韩亚洲高清精品| 国产精品一区二区在线不卡| 高清黄色对白视频在线免费看 | 老司机影院毛片|