• <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ī)制分析
    亚洲人成网站在线观看播放| 一级毛片久久久久久久久女| 日韩中字成人| 人妻制服诱惑在线中文字幕| 日韩成人伦理影院| 97超碰精品成人国产| 国产免费一区二区三区四区乱码| 欧美高清性xxxxhd video| 男女无遮挡免费网站观看| 欧美人与善性xxx| 丰满少妇做爰视频| 搡老乐熟女国产| 大码成人一级视频| 国产成人福利小说| 午夜激情久久久久久久| 一级a做视频免费观看| 亚洲精华国产精华液的使用体验| 免费黄频网站在线观看国产| 日本黄色片子视频| 精品一区二区三区视频在线| 精品久久久久久久末码| 国产成人免费无遮挡视频| 99九九线精品视频在线观看视频| av一本久久久久| 亚洲国产欧美在线一区| 日本wwww免费看| 看非洲黑人一级黄片| 视频区图区小说| 少妇丰满av| 精品国产三级普通话版| 精品久久久久久久久亚洲| 久久99精品国语久久久| 国产高潮美女av| 国产爱豆传媒在线观看| 亚洲欧美精品自产自拍| 国产男女超爽视频在线观看| 精品酒店卫生间| 高清视频免费观看一区二区| 国产伦理片在线播放av一区| 久久这里有精品视频免费| a级一级毛片免费在线观看| 欧美一级a爱片免费观看看| 国产 一区 欧美 日韩| 久久久久性生活片| 亚洲最大成人手机在线| 午夜福利视频精品| 美女内射精品一级片tv| 亚洲精品影视一区二区三区av| 亚洲精品影视一区二区三区av| 在线观看免费高清a一片| 最后的刺客免费高清国语| 丰满少妇做爰视频| 国内精品宾馆在线| 精品人妻熟女av久视频| 亚洲国产色片| 大香蕉97超碰在线| 乱系列少妇在线播放| 亚洲欧美一区二区三区国产| 国产 精品1| 久久久久性生活片| 韩国高清视频一区二区三区| 性色av一级| 韩国高清视频一区二区三区| 蜜桃亚洲精品一区二区三区| 免费在线观看成人毛片| 久热久热在线精品观看| 水蜜桃什么品种好| 亚洲欧美一区二区三区国产| 久久久久国产精品人妻一区二区| 毛片女人毛片| 婷婷色综合www| 能在线免费看毛片的网站| av国产精品久久久久影院| 欧美激情在线99| 国模一区二区三区四区视频| 国模一区二区三区四区视频| 性色avwww在线观看| 男人爽女人下面视频在线观看| 亚洲精品第二区| 久久久久国产精品人妻一区二区| 国产精品久久久久久精品电影小说 | 2021少妇久久久久久久久久久| 2021少妇久久久久久久久久久| 国产精品99久久99久久久不卡 | 国产一区二区亚洲精品在线观看| 亚洲精品国产成人久久av| 91在线精品国自产拍蜜月| 亚洲精品国产av成人精品| 国产精品一区www在线观看| 国产精品久久久久久精品古装| 久久精品国产a三级三级三级| 又爽又黄a免费视频| 男人添女人高潮全过程视频| 内射极品少妇av片p| 联通29元200g的流量卡| 九草在线视频观看| av福利片在线观看| 日韩人妻高清精品专区| 精华霜和精华液先用哪个| 又大又黄又爽视频免费| 成人国产麻豆网| 久久久久久久久大av| 精品久久久久久久久av| 久久久久久久精品精品| 丝袜美腿在线中文| 国产老妇伦熟女老妇高清| 青春草国产在线视频| 中文字幕人妻熟人妻熟丝袜美| 成年女人在线观看亚洲视频 | 亚洲一级一片aⅴ在线观看| 亚洲av一区综合| 国产精品秋霞免费鲁丝片| 啦啦啦在线观看免费高清www| av卡一久久| 一级黄片播放器| 日日啪夜夜爽| 在线播放无遮挡| 熟女电影av网| 国产精品国产三级专区第一集| 最近最新中文字幕免费大全7| 欧美少妇被猛烈插入视频| 一个人观看的视频www高清免费观看| 日韩亚洲欧美综合| 亚洲精品456在线播放app| 国产 一区 欧美 日韩| 久久鲁丝午夜福利片| 麻豆久久精品国产亚洲av| 婷婷色综合大香蕉| 高清av免费在线| 国精品久久久久久国模美| 国产片特级美女逼逼视频| 一本一本综合久久| 久久久久久久大尺度免费视频| 国产日韩欧美在线精品| 欧美高清成人免费视频www| 亚洲成色77777| 人妻夜夜爽99麻豆av| 亚洲av成人精品一二三区| 激情 狠狠 欧美| 国产亚洲av嫩草精品影院| 日本三级黄在线观看| 在线观看一区二区三区激情| 国产日韩欧美在线精品| 国产毛片a区久久久久| 在线观看av片永久免费下载| 亚洲欧美成人综合另类久久久| 国产亚洲最大av| 大片免费播放器 马上看| 国产精品国产三级国产专区5o| 久久久久久久亚洲中文字幕| 伊人久久精品亚洲午夜| 日韩欧美 国产精品| 日本午夜av视频| 免费不卡的大黄色大毛片视频在线观看| 欧美性猛交╳xxx乱大交人| 国语对白做爰xxxⅹ性视频网站| 青春草国产在线视频| av在线老鸭窝| 亚洲性久久影院| 肉色欧美久久久久久久蜜桃 | 最后的刺客免费高清国语| 久久久久性生活片| 日韩电影二区| 免费av不卡在线播放| 国产精品国产三级专区第一集| 亚洲自偷自拍三级| av国产精品久久久久影院| 国产美女午夜福利| 中文字幕亚洲精品专区| 丝袜喷水一区| 色5月婷婷丁香| 日韩人妻高清精品专区| 中国三级夫妇交换| 少妇裸体淫交视频免费看高清| 麻豆成人午夜福利视频| 亚洲国产成人一精品久久久| 午夜精品一区二区三区免费看| 久久久久久久久久成人| 国产精品一区二区性色av| 日本wwww免费看| 欧美性猛交╳xxx乱大交人| 成年人午夜在线观看视频| 国产91av在线免费观看| 中文资源天堂在线| 成年女人看的毛片在线观看| 搡女人真爽免费视频火全软件| 亚洲精品国产成人久久av| 97热精品久久久久久| 性色av一级| 秋霞伦理黄片| av天堂中文字幕网| 亚洲欧美日韩无卡精品| 国产精品国产三级国产av玫瑰| 日韩av免费高清视频| 日本爱情动作片www.在线观看| 亚洲国产色片| 亚洲图色成人| 伦精品一区二区三区| 男的添女的下面高潮视频| 中国三级夫妇交换| 国产精品久久久久久av不卡| 一级av片app| 免费少妇av软件| 岛国毛片在线播放| 亚洲精品久久久久久婷婷小说| 亚洲av男天堂| 久久99热这里只频精品6学生| 亚洲四区av| 成年免费大片在线观看| 欧美潮喷喷水| av福利片在线观看| 3wmmmm亚洲av在线观看| 色哟哟·www| 麻豆成人午夜福利视频| 老师上课跳d突然被开到最大视频| 成人高潮视频无遮挡免费网站| 精品人妻视频免费看| 欧美国产精品一级二级三级 | www.av在线官网国产| 蜜臀久久99精品久久宅男| 国产男人的电影天堂91| 欧美xxxx黑人xx丫x性爽| 在线亚洲精品国产二区图片欧美 | 国产精品一区二区三区四区免费观看| 久久精品人妻少妇| 最近最新中文字幕免费大全7| 亚洲国产最新在线播放| 亚洲久久久久久中文字幕| 亚洲色图av天堂| 国产一区二区在线观看日韩| 国精品久久久久久国模美| 女的被弄到高潮叫床怎么办| 国产成人freesex在线| 亚洲天堂av无毛| 欧美性感艳星| av女优亚洲男人天堂| 小蜜桃在线观看免费完整版高清| av网站免费在线观看视频| 亚洲内射少妇av| 精品人妻熟女av久视频| 日韩电影二区| 国内少妇人妻偷人精品xxx网站| 丰满人妻一区二区三区视频av| 成人无遮挡网站| 亚洲美女搞黄在线观看| 黄色一级大片看看| 街头女战士在线观看网站| 校园人妻丝袜中文字幕| 国产一级毛片在线| 国产毛片a区久久久久| 亚洲av成人精品一二三区| 色播亚洲综合网| 熟妇人妻不卡中文字幕| 99re6热这里在线精品视频| 亚洲真实伦在线观看| 久久人人爽av亚洲精品天堂 | 亚洲精品成人av观看孕妇| 国产伦在线观看视频一区| 国产爱豆传媒在线观看| 国产成人精品婷婷| 欧美日韩国产mv在线观看视频 | 午夜精品国产一区二区电影 | 国产爱豆传媒在线观看| 一级片'在线观看视频| 五月天丁香电影| 可以在线观看毛片的网站| 亚洲精品aⅴ在线观看| 麻豆国产97在线/欧美| 免费观看的影片在线观看| 国产日韩欧美在线精品| 国内少妇人妻偷人精品xxx网站| 国产一区二区亚洲精品在线观看| 亚洲人与动物交配视频| 久久精品熟女亚洲av麻豆精品| 边亲边吃奶的免费视频| 国产男女内射视频| 国产在线男女| 汤姆久久久久久久影院中文字幕| 国产精品av视频在线免费观看| 一区二区三区免费毛片| 久久精品国产鲁丝片午夜精品| 在线天堂最新版资源| 少妇的逼好多水| 免费在线观看成人毛片| 国产探花极品一区二区| 少妇丰满av| 最近手机中文字幕大全| 国产色爽女视频免费观看| 亚洲欧美精品自产自拍| 欧美另类一区| 成年免费大片在线观看| 大香蕉久久网| 国产成人91sexporn| 99久久精品国产国产毛片| 午夜免费观看性视频| 在线观看av片永久免费下载| 黄色怎么调成土黄色| 伊人久久国产一区二区| 一级毛片黄色毛片免费观看视频| 啦啦啦在线观看免费高清www| 韩国av在线不卡| 国产成人freesex在线| 草草在线视频免费看| 春色校园在线视频观看| 可以在线观看毛片的网站| 日产精品乱码卡一卡2卡三| 69人妻影院| 别揉我奶头 嗯啊视频| 国产精品久久久久久久电影| 久久久久久久午夜电影| 五月天丁香电影| 精品久久久久久久人妻蜜臀av| 少妇人妻一区二区三区视频| 91久久精品国产一区二区三区| 性色avwww在线观看| 日韩欧美一区视频在线观看 | 日韩亚洲欧美综合| 在线播放无遮挡| 亚洲av成人精品一二三区| 国语对白做爰xxxⅹ性视频网站| 一级毛片aaaaaa免费看小| 日韩欧美精品v在线| 国产成人免费无遮挡视频| 久久久a久久爽久久v久久| 国产永久视频网站| 国产黄片视频在线免费观看| 免费黄网站久久成人精品| 国产av码专区亚洲av| 亚洲内射少妇av| 美女高潮的动态| 人人妻人人澡人人爽人人夜夜| 大片电影免费在线观看免费| 超碰97精品在线观看| 亚洲成人av在线免费| 又粗又硬又长又爽又黄的视频| 久久精品国产亚洲网站| 欧美激情国产日韩精品一区| www.色视频.com| 熟女人妻精品中文字幕| 国产毛片在线视频| 午夜福利在线在线| 欧美最新免费一区二区三区| 26uuu在线亚洲综合色| 热99国产精品久久久久久7| 午夜免费男女啪啪视频观看| 久久99蜜桃精品久久| 日日摸夜夜添夜夜爱| 久久久久久久亚洲中文字幕| 免费在线观看成人毛片| 黄色一级大片看看| 国产一区亚洲一区在线观看| 久久99热6这里只有精品| 国产白丝娇喘喷水9色精品| 国产精品偷伦视频观看了| 欧美日韩国产mv在线观看视频 | 国产亚洲av片在线观看秒播厂| 亚洲av欧美aⅴ国产| 成人漫画全彩无遮挡| 日本三级黄在线观看| 国内少妇人妻偷人精品xxx网站| 观看美女的网站| 欧美高清性xxxxhd video| 三级国产精品欧美在线观看| 亚洲av一区综合| 国产日韩欧美亚洲二区| 亚洲精品自拍成人| 国产成人a∨麻豆精品| 秋霞伦理黄片| 国产 一区精品| 男人爽女人下面视频在线观看| 欧美高清成人免费视频www| 久久精品国产自在天天线| 五月天丁香电影| 亚洲欧美日韩东京热| 国产精品国产三级国产av玫瑰| 秋霞在线观看毛片| 白带黄色成豆腐渣| 国产亚洲91精品色在线| 成年免费大片在线观看| 国产精品久久久久久久电影| 亚洲av电影在线观看一区二区三区 | 亚洲精品456在线播放app| 永久网站在线| 超碰av人人做人人爽久久| 亚洲天堂国产精品一区在线| 欧美亚洲 丝袜 人妻 在线| av在线播放精品| 青春草视频在线免费观看| 亚洲欧美精品自产自拍| 日韩强制内射视频| 亚洲精品国产av蜜桃| 永久免费av网站大全| 69人妻影院| 午夜视频国产福利| 26uuu在线亚洲综合色| 国产 精品1| 熟女电影av网| 免费大片18禁| 真实男女啪啪啪动态图| av在线观看视频网站免费| 亚洲欧美清纯卡通| 国产伦理片在线播放av一区| 免费电影在线观看免费观看| 黄色视频在线播放观看不卡| 看十八女毛片水多多多| 综合色av麻豆| 真实男女啪啪啪动态图| 亚洲成人中文字幕在线播放| 久久久久国产精品人妻一区二区| 三级男女做爰猛烈吃奶摸视频| 久久ye,这里只有精品| 亚洲美女搞黄在线观看| 久久精品人妻少妇| 一本久久精品| www.av在线官网国产| 男女国产视频网站| 别揉我奶头 嗯啊视频| 国产精品成人在线| 国内精品宾馆在线| 免费大片黄手机在线观看| 精品久久国产蜜桃| 直男gayav资源| 国产午夜精品久久久久久一区二区三区| 最近手机中文字幕大全| 人人妻人人爽人人添夜夜欢视频 | 天堂俺去俺来也www色官网| 网址你懂的国产日韩在线| 欧美性感艳星| 国产伦精品一区二区三区四那| 99久久精品热视频| 亚洲精品成人av观看孕妇| 国产av码专区亚洲av| 免费观看无遮挡的男女| 熟女av电影| 国产成人福利小说| 精品亚洲乱码少妇综合久久| 国产精品偷伦视频观看了| 日韩亚洲欧美综合| www.av在线官网国产| 嫩草影院精品99| 国产免费又黄又爽又色| 交换朋友夫妻互换小说| 91精品伊人久久大香线蕉| 欧美日本视频| 亚洲精品色激情综合| 亚洲美女搞黄在线观看| 欧美日韩视频精品一区| 国产爽快片一区二区三区| 亚洲图色成人| 搞女人的毛片| 日韩av在线免费看完整版不卡| 日本爱情动作片www.在线观看| 99久久精品一区二区三区| 国产一区亚洲一区在线观看| 国产欧美日韩一区二区三区在线 | 99热网站在线观看| 欧美成人精品欧美一级黄| 男人舔奶头视频| 日韩欧美精品免费久久| 亚洲最大成人中文| 中国国产av一级| 夜夜爽夜夜爽视频| 亚洲一级一片aⅴ在线观看| 国产精品伦人一区二区| av在线观看视频网站免费| 黄色视频在线播放观看不卡| 国产熟女欧美一区二区| 国产成人a∨麻豆精品| 99re6热这里在线精品视频| 天美传媒精品一区二区| av专区在线播放| 大片电影免费在线观看免费| 久久久a久久爽久久v久久| 日韩不卡一区二区三区视频在线| 国产一级毛片在线| 在线看a的网站| 亚洲va在线va天堂va国产| 成人亚洲欧美一区二区av| 色播亚洲综合网| 岛国毛片在线播放| 免费观看性生交大片5| 欧美+日韩+精品| 激情 狠狠 欧美| 国产在线男女| 日本猛色少妇xxxxx猛交久久| 国产欧美日韩精品一区二区| 国产成人福利小说| 日韩电影二区| 九色成人免费人妻av| 久久精品熟女亚洲av麻豆精品| 久久人人爽人人爽人人片va| 免费少妇av软件| 国产一区有黄有色的免费视频| 伦精品一区二区三区| 欧美性猛交╳xxx乱大交人| 不卡视频在线观看欧美| 国精品久久久久久国模美| 精品视频人人做人人爽| 久久人人爽人人爽人人片va| 精品久久久久久久人妻蜜臀av| 国产成人福利小说| 久久这里有精品视频免费| 成人特级av手机在线观看| 国产色婷婷99| 欧美少妇被猛烈插入视频| a级一级毛片免费在线观看| 精品久久久久久久久av| 看黄色毛片网站| 偷拍熟女少妇极品色| 夜夜爽夜夜爽视频| 久久久国产一区二区| 亚洲av中文字字幕乱码综合| 久久ye,这里只有精品| 国产精品人妻久久久久久| 美女xxoo啪啪120秒动态图| 18+在线观看网站| a级毛片免费高清观看在线播放| 国产色婷婷99| 有码 亚洲区| 人人妻人人看人人澡| 亚洲av电影在线观看一区二区三区 | 日韩成人伦理影院| 久久综合国产亚洲精品| 国产一区二区三区综合在线观看 | 午夜老司机福利剧场| 亚洲av中文av极速乱| av线在线观看网站| 久久精品综合一区二区三区| 一个人看的www免费观看视频| 国产精品久久久久久精品电影小说 | 哪个播放器可以免费观看大片| 国产精品麻豆人妻色哟哟久久| 久久人人爽av亚洲精品天堂 | 国产亚洲91精品色在线| 日本熟妇午夜| 国产老妇伦熟女老妇高清| 亚洲精品成人久久久久久| 男人添女人高潮全过程视频| 免费观看无遮挡的男女| 天天躁夜夜躁狠狠久久av| 精品99又大又爽又粗少妇毛片| 一级a做视频免费观看| 白带黄色成豆腐渣| 久久久久国产精品人妻一区二区| 久久99热这里只有精品18| kizo精华| 亚洲高清免费不卡视频| 婷婷色麻豆天堂久久| 新久久久久国产一级毛片| 日韩强制内射视频| 你懂的网址亚洲精品在线观看| 免费看av在线观看网站| 777米奇影视久久| 欧美激情久久久久久爽电影| 精品人妻熟女av久视频| 菩萨蛮人人尽说江南好唐韦庄| 亚洲精品一二三| 男插女下体视频免费在线播放| 亚洲精品,欧美精品| 国产一区有黄有色的免费视频| 久久久久久久精品精品| 男的添女的下面高潮视频| 一级a做视频免费观看| 欧美变态另类bdsm刘玥| 精品久久久精品久久久| 精品人妻视频免费看| 久久久久久久亚洲中文字幕| av在线播放精品| 日韩av不卡免费在线播放| 久久精品国产a三级三级三级| 久久影院123| 国产午夜福利久久久久久| 又爽又黄无遮挡网站| 啦啦啦在线观看免费高清www| 女人久久www免费人成看片| 国产中年淑女户外野战色| 最近最新中文字幕免费大全7| 精品少妇黑人巨大在线播放| 午夜爱爱视频在线播放| 精品一区二区免费观看| 伦精品一区二区三区| 婷婷色综合www| 精品久久久久久久久亚洲| 成人欧美大片| 日韩欧美精品免费久久| 2021天堂中文幕一二区在线观| 男女那种视频在线观看| 欧美zozozo另类| 精品人妻偷拍中文字幕| 欧美高清性xxxxhd video| 久久影院123| 大话2 男鬼变身卡| 亚洲av不卡在线观看| 男女国产视频网站| 精品熟女少妇av免费看| 如何舔出高潮| 国产一区二区亚洲精品在线观看| 精品国产露脸久久av麻豆| 国产v大片淫在线免费观看| 国产黄a三级三级三级人| 大陆偷拍与自拍| 香蕉精品网在线| 综合色av麻豆| 国产免费福利视频在线观看| 亚洲精品乱码久久久v下载方式| 偷拍熟女少妇极品色| 交换朋友夫妻互换小说| 国产乱人偷精品视频| 欧美性感艳星| 身体一侧抽搐| 亚洲精品国产av蜜桃| 欧美性猛交╳xxx乱大交人| 国产免费福利视频在线观看| 亚洲图色成人| 校园人妻丝袜中文字幕|