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

    內(nèi)蒙古東勝地區(qū)砂巖型鈾礦賦礦地層方解石膠結(jié)物C、O同位素特征和成因模型

    2021-07-28 12:02:58劉漢彬金貴善張建鋒韓娟李軍杰張佳石曉徐可
    地質(zhì)論評(píng) 2021年4期
    關(guān)鍵詞:直羅東勝方解石

    劉漢彬,金貴善,張建鋒,韓娟,李軍杰,張佳,石曉,徐可

    核工業(yè)北京地質(zhì)研究院,北京,100029

    內(nèi)容提要: 砂巖型鈾礦賦礦地層方解石膠結(jié)物C、O同位素組成對(duì)研究成礦機(jī)制具有重要的理論和實(shí)踐意義。前人在成礦流體來(lái)源定性判斷方面做了大量的研究工作,但缺乏對(duì)成礦混合流體的組分特征及相關(guān)參數(shù)的定量認(rèn)識(shí)。根據(jù)內(nèi)蒙古東勝地區(qū)東南部皂火壕和西北部納嶺溝鈾礦床賦礦地層直羅組方解石膠結(jié)物C、O同位素組成特征,建立了滲出熱鹵水和滲入地層水兩種不同流體混合生成方解石膠結(jié)物C、O同位素組成定量成因模型,并研究了兩種流體來(lái)源、溶解碳的濃度比、流體比例以及熱液溫度等綜合因素。除少數(shù)樣品方解石膠結(jié)物為地層沉積時(shí)形成外,大部分方解石膠結(jié)物是熱鹵水與地層水兩種流體混合作用的結(jié)果,古生界有機(jī)酸脫羧作用導(dǎo)致熱鹵水中富含CO2。東勝地區(qū)砂巖型鈾礦成礦流體地層水與熱鹵水比例為0.5~0.9,地層水所占的比例較大;熱鹵水與地層水溶解碳濃度比主要范圍為1.5~5.0,部分大于10.0,熱鹵水中溶解碳濃度較高,是富含CO2的流體。混合流體溫度分為兩個(gè)主要范圍55~80℃、90~140℃,結(jié)合古鹽度和鹽度指數(shù)反演,表明成礦流體為有機(jī)與無(wú)機(jī)混合成因的低溫?zé)嵋毫黧w。用砂巖型鈾礦賦礦地層方解石膠結(jié)物C、O同位素組成,可定量模擬兩種流體溶解碳濃度比、流體比例和方解石形成時(shí)溫度等流體成礦條件,解釋碳酸鹽膠結(jié)物成因,以便從成礦流體角度更好地理解成礦作用過(guò)程。

    鄂爾多斯盆地是一個(gè)重要的石油、天然氣、煤炭能源基地。自從20世紀(jì)90年代末,在該盆地東北部東勝地區(qū)先后發(fā)現(xiàn)并探明了皂火壕、納嶺溝、大營(yíng)等砂巖型鈾礦床(張金帶,2016)。鈾礦床成因模式為疊合鈾成礦模式,即在古層間氧化作用的基礎(chǔ)上,疊加后期還原性低溫?zé)崃黧w改造而成礦(李子穎等,2006;秦明寬等,2017;彭云彪等,2017;易超等,2018)。

    目前,砂巖型鈾礦賦礦地層方解石膠結(jié)物C、O同位素組成已經(jīng)在東勝地區(qū)成礦作用研究中得到很好的利用,研究手段主要是根據(jù)不同成因C同位素組成分析、δ13CV-PDB及δ18OV-SMOW圖解投圖、用方解石膠結(jié)物和流體之間的氧同位素分餾系數(shù)計(jì)算流體氧同位素組成等方法,對(duì)其成因進(jìn)行判斷分析,獲得對(duì)流體來(lái)源的認(rèn)識(shí)(馮喬等,2016;吳柏林等,2016;湯超等,2016;張賓等,2020):C同位素組成主要受有機(jī)酸脫羥基作用的控制,O同位素組成揭示其來(lái)源于地層水與低溫深層熱鹵水的混合。但模式比較法只能進(jìn)行定性分析,利用分餾系數(shù)計(jì)算流體O同位素組成時(shí),只把溫度作為限定條件,沒(méi)有限定兩種流體混合比例以及各自的同位素組成等條件,制約了碳酸鹽膠結(jié)物C、O同位素組成在砂巖型鈾礦床成礦作用過(guò)程中的示蹤應(yīng)用。

    在對(duì)東勝地區(qū)東南部皂火壕鈾礦床賦礦地層直羅組方解石膠結(jié)物C、O同位素特征進(jìn)行分析基礎(chǔ)上,試圖根據(jù)已有的研究資料,假定地層水和低溫?zé)猁u水的C、O同位素組成,用兩種不同流體混合生成方解石C、O同位素組成計(jì)算方法(鄭永飛等,2000),從滲入地層水、滲出低溫?zé)猁u水兩種流體來(lái)源、溶解碳的濃度比、流體比例以及熱液溫度等因素綜合考慮,建立定量成因模型,分析并解釋碳酸鹽膠結(jié)物成因,以便從成礦流體角度更好地理解成礦作用過(guò)程。

    1 地質(zhì)概況

    鄂爾多斯盆地是在華北克拉通基礎(chǔ)上,受古亞洲、特提斯和環(huán)太平洋三大構(gòu)造域的作用交切、疊加與復(fù)合而成(李子穎等,2019)。盆地劃分為伊盟隆起、渭北隆起、晉西撓褶帶、陜北斜坡、天環(huán)坳陷和西緣沖斷帶6個(gè)一級(jí)構(gòu)造單元(李宏濤等,2007)。東勝地區(qū)砂巖型鈾礦位于伊盟隆起的東北部。盆地由結(jié)晶基底、穩(wěn)定型地臺(tái)沉積蓋層和活動(dòng)型盆地沉積蓋層組成。其中,結(jié)晶基底為太古宇—古元古界,地臺(tái)沉積蓋層為中元古界—三疊系,盆地蓋層為侏羅系—第四系(圖1)。

    圖1 鄂爾多斯盆地東勝地區(qū)砂巖型鈾礦地質(zhì)簡(jiǎn)圖(據(jù)陳路路等,2017修改)Fig. 1 Geological map of sandstone type uranium deposit in Dongsheng area, Ordos Basin (modified from Chen Lulu et al., 2017&)

    侏羅系直羅組為東勝地區(qū)砂巖型鈾礦賦礦地層,直羅組分為上段和下段,上段為雜色的中細(xì)砂巖、泥巖。下段分為上、下兩個(gè)亞段(焦養(yǎng)泉等,2005),上亞段主要為綠色砂巖、泥巖;含礦主巖為直羅組下段的下亞段砂巖,巖石顏色以灰色為主,結(jié)構(gòu)較疏松。直羅組砂巖總體成熟度較低,砂巖中碎屑磨圓度和分選性均較差,分為粗砂巖、中粒砂巖、細(xì)砂巖。

    直羅組砂巖氧化還原分帶具有其獨(dú)特特征,在垂向上具有后生蝕變分帶(李子穎等,2007;肖新建等,2004a),直羅組下段上亞段綠色砂巖為氧化帶,是發(fā)生后生蝕變作用的主要分帶。直羅組下段下亞段灰色砂巖為還原帶,綠色砂巖與灰色砂巖之間為氧化—還原過(guò)渡帶,鈾礦化位于過(guò)渡帶內(nèi)。

    2 方解石膠結(jié)物同位素組成

    2.1 方解石礦物特征

    東勝地區(qū)砂巖型鈾礦一個(gè)重要的成礦特征是發(fā)生較強(qiáng)的碳酸鹽化蝕變作用,碳酸鹽膠結(jié)物含量一般小于0.5%,因膠結(jié)物含量較少,碎屑顆粒相互接觸,膠結(jié)物僅存在于碎屑顆粒的接觸處,以接觸膠結(jié)為主,孔隙式膠結(jié)為輔。部分砂巖碳酸鹽化強(qiáng)烈,碳酸鹽含量可達(dá)10%~20%或更高,砂巖基質(zhì)和膠結(jié)物含量相對(duì)較高,碎屑顆粒分布于其中,呈基底式膠結(jié)。據(jù)X-射線衍射分析,碳酸鹽幾乎全為方解石,有少量菱鐵礦(李子穎等,2019)。

    在賦礦地層直羅組砂巖特別是礦化段砂巖中,碳酸鹽化較強(qiáng),膠結(jié)物為亮晶方解石或粗晶方解石。碎屑顆粒主要有石英、長(zhǎng)石和巖屑等,結(jié)構(gòu)成熟度及成分成熟度為中等。碎屑物中少部分石英顆粒出現(xiàn)次生加大邊現(xiàn)象(圖2a)。填隙物主要以方解石膠結(jié)物為主,另外,方解石常與黃鐵礦伴生(圖2)。

    圖2 東勝地區(qū)皂火壕鈾礦床直羅組砂巖特征: (a)方解石(Cal)膠結(jié)物膠結(jié)石英(Qtz)、長(zhǎng)石(Or)顆粒(正交偏光);(b)黃鐵礦(Py)與方解石(Cal)(反射光);(c)粒狀黃鐵礦(Py)充填顆??紫叮浇馐?Cal)膠結(jié)物膠結(jié)石英(Qtz)顆粒,鈦鈾礦(Bra)(反射光);(d)浸染狀黃鐵礦(Py)與方解石(Cal)膠結(jié)物(反射光)Fig. 2 Sandstone characteristics of Zhiluo Formation in Zaohuohao uranium deposit,Dongsheng area: (a) Calcitecement cemented quartz and feldspar particles ( magnification , crossed-polarized light); (b) pyriteand calcite ( magnification, reflected light); (c) granular pyrite filled in particle pores, calcite cement colloidal quartz particles, uraninite (× 5 magnification , reflected light); (d) disseminated pyrite and calcite cement (× 5 magnification, reflected light) 方解石—Cal; 石英—Qtz; 長(zhǎng)石—Or; 黃鐵礦— Py; 鈦鈾礦—Bra Calcite—Cal; Quartz—Qtz; Feldspar—Or; Pyrite— Py; Uraninite—Bra

    為研究東勝地區(qū)成礦流體作用,特別是成礦流體來(lái)源及組成、成礦溫度及氧化還原不同分帶差異等,在研究區(qū)東南部皂火壕礦床進(jìn)行賦礦地層直羅組砂巖取樣,在對(duì)砂巖方解石膠結(jié)物C、O同位素組成分析測(cè)試基礎(chǔ)上,進(jìn)行了同位素特征和成因模型研究。

    2.2 分析過(guò)程

    分析測(cè)試過(guò)程主要為,先將樣品粉碎至180目,在真空、25℃恒溫條件下,將適量100%的磷酸加入盛有2 g砂巖樣品的反應(yīng)試管中,反應(yīng)時(shí)間4 h以上, 用冷凍法分離、除去生成的水,收集由方解石膠結(jié)物與磷酸反應(yīng)生成的純凈CO2氣體。在中國(guó)地質(zhì)科學(xué)院礦產(chǎn)資源研究所用MAT251氣體同位素質(zhì)譜儀進(jìn)行雙路法測(cè)量,測(cè)得C、O同位素組成數(shù)據(jù)(表1)。測(cè)量結(jié)果以V-PDB和V-SMOW為標(biāo)準(zhǔn), 分別記為δ13CV-PDB和δ18OV-PDB、δ18OV-SMOW。分析精度優(yōu)于±0. 2‰。碳酸鹽參考標(biāo)準(zhǔn)為GBW04417方解石, 其δ18OV-PDB和δ13CV-PDB分別是-24.12‰ 和-6.06‰(丁悌平,2002)。δ18OV-PDB與δ18OV-SMOW兩種表示之間的轉(zhuǎn)化關(guān)系為,δ18OV-SMOW=1.03091×δ18OV-PDB+30.91。

    2.3 C、O同位素組成特征

    對(duì)分析數(shù)據(jù)進(jìn)行統(tǒng)計(jì)(表1),東勝地區(qū)皂火壕鈾礦床方解石膠結(jié)物C、O同位素組成范圍:

    表1 東勝地區(qū)皂火壕鈾礦床方解石膠結(jié)物C、O同位素組成Table 1 Carbon and oxygen isotopic compositions of calcite cements in Zaohuohao uranium deposit, Dongsheng area

    (1)灰色及綠色含礦砂巖的δ13CV-PDB范圍為-12.5‰~-1.5‰,平均值為-7.8‰;δ18OV-PDB范圍為-12.9‰~-9.7‰,平均為-11.3‰,δ18OV-SMOW范圍為17.5‰ ~22.0‰,平均值為19.3‰。

    (2)灰色砂巖的δ13CV-PDB范圍為-12.1‰~-1.5‰,平均值為-7.3‰;δ18OV-PDB范圍為-14.0‰~-5.0‰,平均值為-10.6‰,δ18OV-SMOW范圍為16.5‰ ~25.7‰,平均值為20.0‰。

    (3)綠色砂巖的δ13CV-PDB范圍為-17.8‰~-1.6‰,平均值為-6.8‰;δ18OV-PDB范圍為-13.6‰~-7.5‰,平均值為-10.4‰,δ18OV-SMOW范圍為16.9‰~23.1‰ ,平均值為20.2‰。全部樣品的δ13CV-PDB范圍為-17.8‰~-1.5‰,平均值為-7.3‰;δ18OV-PDB范圍為-14.0‰~-5.0‰,平均值為-10.7‰,δ18OV-SMOW范圍為16.5‰ ~25.7‰,平均值為19.9‰。

    可以看出,除個(gè)別綠色砂巖的方解石膠結(jié)物δ13CV-PDB值相對(duì)較低,個(gè)別灰色砂巖的方解石膠結(jié)物 O同位素值相對(duì)較低外,三種不同砂巖的方解石膠結(jié)物C、O同位素組成范圍基本一致。

    2.4 流體來(lái)源

    根據(jù)方解石膠結(jié)物C、O同位素組成,可以判斷砂巖型鈾礦成礦流體的來(lái)源(楊曉勇等,2009;薛春紀(jì)等,2010;馮喬等,2016;吳柏林等,2016;湯超等,2016;李子穎等,2019;張賓等,2020)。本次測(cè)量的砂巖方解石膠結(jié)物δ13CV-PDB范圍為-17.8‰~-1.5‰,總體上比海相碳酸鹽C同位素組成(δ13CV-PDB=-4.0‰~4.0‰; Bath et al.,1987)偏負(fù),也比大氣降水中C同位素組成(δ13CV-PDB=-7.0‰~-3.5‰; Schmidt et al.,2005)偏負(fù),但比有機(jī)質(zhì)的C同位素組成(δ13CV-PDB=-33.0‰~-18.0‰; Schmidt et al.,2005)偏正,與上古生界氣田中CO2碳同位素組成范圍(δ13CV-PDB=-20‰~-4‰; 張龍,2017)大致相當(dāng),且約有五分之一樣品的δ13CV-PDB<-10.0‰。因此,可以推測(cè)盆地?zé)猁u水中溶解碳參與了砂巖膠結(jié)物方解石形成。

    在東勝地區(qū)很多學(xué)者對(duì)賦礦地層直羅組方解石膠結(jié)物進(jìn)行了C、O同位素組成特征研究(表2),馮喬等(2016)對(duì)東勝鈾礦區(qū)方解石C、O同位素研究表明,方解石的δ13CV-PDB值范圍為-19.6‰ ~-1.11‰,δ18OV-PDB值范圍-17.13‰~-9.00‰,由于地層滲入水和滲出熱鹵水影響程度不同,導(dǎo)致δ13C值的變化范圍較寬,方解石是由低溫?zé)嵋撼傻V作用形成的。楊曉勇等(2009)對(duì)碳酸鹽膠結(jié)物的C、O同位素測(cè)定結(jié)果表明,δ13CV-PDB值變化范圍為-23.62‰ ~-4. 24‰,δ18OV-SMOW范圍為10. 32‰ ~19. 52‰,且碳酸鹽含量與δ13CV-PDB值呈負(fù)相關(guān)關(guān)系,說(shuō)明有機(jī)碳加入使砂巖中碳酸鹽膠結(jié)物含量增加,碳同位素組成δ13CV-PDB反而變負(fù),證明了賦礦地層直羅組碳酸鹽膠結(jié)物方解石部分來(lái)源于有機(jī)質(zhì)脫羧作用。吳柏林等(2016)對(duì)方解石樣品的C、O 同位素組成特征進(jìn)行了研究,δ13CV-PDB范圍為- 0.3‰~-14.59‰,平均值為-6.19‰,δ18OV-PDB范圍為-9.6‰~-12.31‰,平均值為-11.34‰,其碳的來(lái)源具有明顯有機(jī)碳特征。薛春紀(jì)等(2010)研究表明,方解石的δ13CV-PDB范圍為-14. 0‰~ 2.7‰,δ18OV-SMOW范圍為18.4‰ ~24.6‰,認(rèn)為碳同位素組成指示鈾成礦流體中的大量碳來(lái)自烴類的氧化。湯超等(2016)研究結(jié)果表明,方解石膠結(jié)物δ13CV-PDB范圍為-15.7‰~-1.6‰,平均值為-9.08‰,δ18OV-PDB范圍為-15.6‰~-10‰,平均-12.4‰,方解石的碳來(lái)源為有機(jī)碳,與有機(jī)酸脫羧作用有關(guān)。

    表2 東勝地區(qū)砂巖型鈾礦直羅組碳酸鹽膠結(jié)物 方解石C、O同位素組成范圍

    另外,還有學(xué)者針對(duì)東勝地區(qū)砂巖型鈾礦成礦具有“疊合成因”的特殊性,按砂巖不同顏色分類進(jìn)行了研究。李子穎等(2019)研究了東勝地區(qū)賦礦地層直羅組方解石膠結(jié)物C、O同位素組成,其中,灰色砂巖的δ13CV-PDB范圍為-21.9‰~-7.9‰,δ18OV-PDB范圍為-12.9‰~-9.7‰;灰綠色砂巖的δ13CV-PDB范圍為-31.0‰~-0.4‰,δ18OV-PDB范圍為-13.8‰~-9.9‰;褐紅色砂巖的δ13CV-PDB范圍為-20.9‰~-1.9‰,δ18OV-PDB范圍為-12.8‰~-8.9‰??傮w上,砂巖的δ13CV-PDB范圍為-31.0‰~-0.4‰,表現(xiàn)出差別較大的特征;δ18OV-PDB范圍為-13.8‰~-8.9‰,具有相對(duì)集中的特征,說(shuō)明成礦流體來(lái)源于油氣、地層熱鹵水的雙重參與,具有有機(jī)、無(wú)機(jī)流體雙重成因。

    張賓等(2020)研究了東勝地區(qū)西北部納嶺溝鈾礦床直羅組碳酸鹽膠結(jié)物C、O同位素組成特征,含礦砂巖δ13CV-PDB為-14.0‰~-8.7‰,平均值為-11.12‰,δ18OV-PDB為-13.4‰~-12.4‰,平均值為-12.8‰;灰色砂巖δ13CV-PDB為-14.3‰~-3.8‰,平均值為-8. 44‰,δ18OV-PDB為-14. 3‰~-7. 9‰,平均值為-11. 41‰;綠色砂巖δ13CV-PDB為-20. 6‰~-5.2‰,平均值為-10. 28‰,δ18OV-PDB為-16.3‰~-7.0‰,平均值為-11.61‰。有機(jī)酸脫羧作用產(chǎn)生的方解石δ13CV-PDB以高負(fù)值為特征,大部分樣品方解石膠結(jié)物中的碳主要是有機(jī)質(zhì)脫羧作用形成(圖3)。湯超等(2016)研究東勝地區(qū)砂巖型鈾礦膠結(jié)物方解石C、O同位素結(jié)果表明,方解石的形成與有機(jī)質(zhì)脫羧作用有關(guān),碳來(lái)源為有機(jī)碳(圖3)。

    圖3 方解石膠結(jié)物δ 13CV-PDB-δ18OV-SMOW圖解(據(jù)劉建明等,1997)Fig. 3 Diagram of δ 13CV-PDB-δ18OV-SMOW in calcite cement (modified from Liu Jianming et al., 1997&)

    綜上所述,東勝地區(qū)砂巖型鈾礦賦礦地層直羅組方解石膠結(jié)物C、O同位素研究表明,除個(gè)別樣品呈較為偏負(fù)的C同位素組成外,大部分樣品的δ13CV-PDB范圍在-20.6‰~0.3‰之間,δ18OV-PDB范圍在-20.0‰~-6.1‰,δ18OV-SMOW范圍在10.3‰~24.6‰,本次測(cè)量數(shù)據(jù)范圍與前人分析數(shù)據(jù)范圍基本一致。

    在方解石膠結(jié)物δ13CV-PDB-δ18OV-SMOW圖解上(圖3),樣品點(diǎn)主要分布在沉積有機(jī)物脫羧基作用形成的碳酸鹽分布范圍內(nèi),有10個(gè)樣品分布在海相碳酸鹽范圍內(nèi)的左下部。在方解石膠結(jié)物碳氧同位素組成成因圖解上(圖4),樣品點(diǎn)主要分布在與有機(jī)質(zhì)有關(guān)的碳酸鹽及成巖碳酸鹽分布范圍內(nèi)。這說(shuō)明碳酸鹽膠結(jié)物方解石的碳主要來(lái)源是有機(jī)酸脫羧基作用提供的,還有部分方解石(C、O同位素組成相對(duì)較大)主要是在地層沉積時(shí)形成的。

    圖4 方解石膠結(jié)物C、O同位素組成成因圖解(據(jù)王大銳,2000)Fig. 4 Origin diagram of carbon and oxygen isotopic composition of calcite cement (modified from Wang Darui,2000&)

    3 成因模型及模型判斷

    3.1 模型建立

    如果流體混合作用導(dǎo)致方解石沉淀,那么可以根據(jù)A、B兩種流體的C、O同位素組成,計(jì)算方解石C、O同位素組成(鄭永飛等,2000):

    (1)

    δ18O方解石=δ18Ob+103lnα方解石—H2O+

    XA(δ18Oa-δ18Ob)

    (2)

    式中:XA為混合流體中流體A的摩爾分?jǐn)?shù);P為流體B與流體A之間全部溶解碳的濃度比;a、b分別表示流體A、流體B;δ13Ca為流體A的碳同位素組成;δ13Cb為流體B的碳同位素組成;αm—n為組分m和n之間的分餾系數(shù),1000 lnαm—n為組分m和n之間的富集系數(shù)(Ohmoto and Rye, 1979)。

    表3東勝地區(qū)砂巖型鈾礦賦礦地層直羅組方解石膠結(jié)物C、O同位素組成模擬參數(shù)Table 3 Simulation parameters of carbon and oxygen isotopic compositions of calcite cements in Zhiluo Formation in sandstone type uranium deposit in Dongsheng area

    1000 lnα方解石—CO2=2.962-11.346×103/T+5.358×106/T2-0.388×109/T3

    方解石與H2O之間的氧同位素平衡分餾系數(shù)由下列方程表示:

    1000lnα方解石—H2O=2.78×106/T2-2.89

    (3)

    因此,給定XA、P、T參數(shù)條件,XA=0、T=250 ℃;XA=0.25、T=200℃;XA=0.5、T=150℃;XA=0.75、T=100℃;XA=1、T=50℃時(shí),分別對(duì)應(yīng)取值P=0.2、0.5、1.0、2.0、5.0、10.0,另外,根據(jù)已有地質(zhì)資料,可假設(shè)或根據(jù)已有研究資料確定兩種流體的C、O同位素組成,用上述公式分別計(jì)算出方解石C、O同位素組成(表3),這樣制作成δ13CV-PDB—δ18OV-SMOW成因圖解,通過(guò)圖解能綜合判斷方解石形成的流體參數(shù)條件,進(jìn)而對(duì)其成因進(jìn)行分析。

    薛春紀(jì)等(2010)通過(guò)盆地內(nèi)流體流動(dòng)的數(shù)值模擬研究表明,在侏羅紀(jì)—白堊紀(jì)時(shí),鄂爾多斯盆地邊緣相對(duì)抬升導(dǎo)致地層發(fā)生輕微變傾,產(chǎn)生滲出流和滲入流兩種流體流動(dòng)系統(tǒng)。滲出流系統(tǒng)是盆地下部受超壓驅(qū)動(dòng)流向盆地邊緣淺部、富含烴、碳的熱鹵水,滲入流系統(tǒng)是受重力驅(qū)動(dòng)沿盆邊地層向下流的地層水。

    吳柏林等(2006)對(duì)東勝地區(qū)蝕源區(qū)水的H、O同位素組成進(jìn)行研究,其δDV-SMOW為-54. 6‰,δ18OV-SMOW為-0.15‰。研究區(qū)白堊系地層地下水溶解無(wú)機(jī)碳δ13CV-PDB實(shí)測(cè)值范圍-19.3‰~-4.5‰(吳春勇等,2011),盆地北部上古生界氣田中CO2碳同位素組成δ13CV-PDB主要分布在-20‰~-4‰(張龍,2017)。湯超等(2016)研究成礦期混合流體同位素組成,混合流體δ18OV-SMOW的平均值為7.19‰。根據(jù)上述資料綜合,假定滲入流體的C、O同位素組成分別為δ18OV-SMOW=0.0‰;δ13CV-PDB=0.0‰;滲出熱鹵水流體的C、O同位素組成分別為δ18OV-SMOW=7.0‰;δ13CV-PDB=-20.0‰(表3)。

    在此基礎(chǔ)上進(jìn)行模擬(圖5)并對(duì)成礦條件加以判斷。在公式計(jì)算和模擬圖中,XA代表東勝地區(qū)砂巖型鈾礦賦礦地層直羅組混合流體中地層水的摩爾分?jǐn)?shù);P代表熱鹵水與地層水之間全部溶解碳的濃度比;a、b分別代表地層水、熱鹵水。

    圖5 東勝地區(qū)砂巖型鈾礦賦礦地層直羅組方解石膠結(jié)物C、O同位素組成模型Fig. 5 Carbon and oxygen isotopic compositions model of calcite cements in Zhiluo Formation in sandstone type uranium deposit in Dongsheng area 假定滲入地層水的C、O同位素組成分別為δ18OV-SMOW=0.0‰;δ13CV-PDB=0.0‰;滲出熱鹵水流體的C、O同位素組成分別為δ18OV-SMOW=7.0‰;δ13CV-PDB=-20.0‰。XA代表混合流體中地層水的摩爾分?jǐn)?shù);P代表熱鹵水與地層水之間全部溶解碳的濃度比 It is assumed that the isotopic compositions of carbon and oxygen in the seeping formation water are that δ18OV-SMOW is 0.0‰ and δ13CV-PDB is 0.0‰; the isotopic compositions of carbon and oxygen in the seeping hot brine fluid are that δ18OV-SMOW is 7.0‰; δ13CV-PDB is -20.0‰. XA represents the mole fraction of formation water in the mixed fluid; P represents the concentration ratio of all dissolved carbon between hot brine and formation water

    從公式(2)、(3)可以看出,方解石C、O同位素組成是兩種同位素在方解石與混合流體之間分餾作用的綜合體現(xiàn),受兩種流體來(lái)源、溶解碳濃度比、流體比例以及熱液溫度等多個(gè)因素影響,不僅僅受盆地?zé)猁u水中碳同位素組成控制,因此,需根據(jù)成礦條件綜合考慮,建立并利用成因模型,從更全面的對(duì)賦礦地層方解石膠結(jié)物成因進(jìn)行深入分析。

    從模擬圖可以看出(圖5),隨混合流體中地層水的摩爾分?jǐn)?shù)XA即地層水的比例增大,方解石的δ18OV-SMOW增大;混合流體溫度降低,方解石的δ18OV-SMOW增大;隨熱鹵水與地層水之間全部溶解碳的濃度比P增大,方解石C同位素組成δ13CV-PDB減小。從熱鹵水中溶解碳含量濃度比P的變化情況,也可對(duì)熱鹵水上升期次做出粗略判斷。方解石膠結(jié)物O同位素組成δ18OV-SMOW的變化對(duì)溫度較靈敏。

    3.2 模型判斷

    將皂火壕鈾礦床賦礦地層直羅組方解石膠結(jié)物C、O同位素組成在建立的模型圖上按灰色及綠色含礦砂巖、灰色砂巖、綠色砂巖、全部砂巖分別投圖,并進(jìn)行討論。

    所有含礦砂巖測(cè)試樣品點(diǎn)都落在兩種流體混合模擬曲線范圍內(nèi),混合流體中地層水的摩爾分?jǐn)?shù)約為0.6~0.9,即地層水的比例約占2/3以上,總體上,地層水在混合流體中所占比例較大。熱鹵水與地層水之間全部溶解碳的濃度比P為1.5~5.0,即熱鹵水中溶解碳的濃度是地層水中溶解碳的濃度的1.5至5倍左右。方解石形成溫度具有2個(gè)溫度段,分別為60~80℃和100~110℃(圖6)。形成溫度主要在90~110℃范圍內(nèi),與灰色砂巖、綠色砂巖相比,溫度不是最高,可能這與兩種流體混合后溫度會(huì)相應(yīng)降低有關(guān)。

    圖6 東勝地區(qū)皂火壕砂巖型鈾礦直羅組含礦砂巖方解石膠結(jié)物C、O同位素組成模型Fig. 6 Carbon and oxygen isotopic compositionsmodel of calcite cements in Zhiluo Formation in Zaohuohaosandstone type uranium deposit in Dongsheng area

    灰色砂巖測(cè)試樣品點(diǎn)大部分落在兩種流體混合模擬曲線范圍內(nèi),混合流體中地層水的摩爾分?jǐn)?shù)約為0.5~0.9,即地層水的比例約占一半以上,地層水在混合流體中所占比例相對(duì)較大。熱鹵水與地層水之間全部溶解碳的濃度比P為1.5~5.0,部分樣品在10.0以上。熱鹵水中溶解碳的濃度是地層水中溶解碳的濃度的1.5至5倍左右,部分樣品落在曲線以外,說(shuō)明熱鹵水中溶解碳的濃度是地層水中溶解碳的濃度的10倍以上,熱鹵水中含碳濃度大大高于地層水。方解石形成溫度也具有2個(gè)溫度段,分別為50~80℃和90~140℃(圖7)。形成溫度主要在90~120℃范圍內(nèi),與含礦砂巖、綠色砂巖相比,個(gè)別樣品溫度達(dá)到最高,最高溫度約為135℃,可能是地層水所占比例較小(0.5左右),熱鹵水所占比例相對(duì)較大,兩種流體混合后溫度會(huì)比原地層水的溫度相應(yīng)增高。

    圖7 東勝地區(qū)皂火壕砂巖型鈾礦直羅組灰色砂巖方解石膠結(jié)物C、O同位素組成模型Fig. 7 Carbon and oxygen isotopic compositions model of calcite cement in gray sandstone in Zhiluo Formation in Zaohuohao sandstone type uranium deposit in Dongsheng

    綠色砂巖樣品基本落在兩種流體混合模擬曲線范圍內(nèi),混合流體中地層水的摩爾分?jǐn)?shù)約為0.6~0.9,即地層水的比例約占2/3以上,地層水在混合流體中所占比例也較大。熱鹵水與地層水之間全部溶解碳的濃度比P主要范圍為2.0~10.0,即熱鹵水中溶解碳的濃度是地層水中溶解碳的濃度的2至10倍左右。方解石形成溫度具有2個(gè)溫度段,分別為較低溫度段55~80℃和較高溫度段90~130℃(圖8),最高溫度約為130℃,最高溫度比含礦砂巖溫度高,比灰色砂巖低。

    圖8 皂火壕砂巖型鈾礦直羅組綠色砂巖方解石膠結(jié)物C、O同位素組成模型Fig. 8 Carbon and oxygen isotopic compositions model of calcite cement in green sandstone in Zhiluo Formation in Zaohuohao sandstone type uranium deposit

    皂火壕砂巖型鈾礦直羅組砂巖方解石膠結(jié)物C、O同位素組成成因(圖9)說(shuō)明,該礦床混合流體中地層水的摩爾分?jǐn)?shù)約為0.5~0.9,即地層水的比例約占一半以上,地層水在混合流體中所占比例較大,但在含礦砂巖、灰色砂巖、綠色砂巖之間所占比例有所差異。熱鹵水與地層水之間全部溶解碳的濃度比P主要范圍為1.5~5.0,即熱鹵水中溶解碳的濃度是地層水中溶解碳的濃度的2至5倍左右,個(gè)別樣品達(dá)到了10倍左右。方解石形成溫度具有2個(gè)溫度段,分別為較低溫度段55~80℃和較高溫度段90~140℃(表4)。

    圖9 納嶺溝砂巖型鈾礦直羅組膠結(jié)物方解石C、O同位素組成成因Fig. 9 Carbon and oxygen isotopic composition genesis of calcite in cement in Zhiluo Formation in Nalinggou sandstone type uranium deposit

    表4 直羅組方解石膠結(jié)物C、O同位素模擬 流體成分及溫度參數(shù)

    通過(guò)上述模擬判斷綜合說(shuō)明,賦礦地層中含礦砂巖、灰色砂巖、綠色砂巖都存在熱鹵水及地層水兩種流體的混合作用,只不過(guò)兩種流體的比例有所差別。地層中方解石膠結(jié)物的成因不能簡(jiǎn)單分為有機(jī)成因還是無(wú)機(jī)成因,事實(shí)上是兩種流體共同作用的結(jié)果,只是兩種流體在不同部位對(duì)方解石的形成所起的作用不同,成礦類型為雙重成因。

    需要說(shuō)明的是,灰色砂巖、綠色砂巖中都有部分樣品熱鹵水與地層水流體溶解碳濃度比在10.0以上,這可能是成礦后期直羅組砂巖經(jīng)受了一期更加富含CO2熱鹵水的蝕變作用。

    3.3 納嶺溝鈾礦床模型 判斷

    為說(shuō)明整個(gè)東勝地區(qū)區(qū)域研究的代表性,根據(jù)張賓等(2020)在東勝地區(qū)西北部納嶺溝鈾礦床賦礦地層直羅組方解石膠結(jié)物C、O同位素組成分析數(shù)據(jù),利用建立的模型對(duì)該礦床進(jìn)行有關(guān)參數(shù)研究(圖9)。

    該礦床含礦砂巖、灰色砂巖、綠色砂巖測(cè)試樣品點(diǎn)大部分落在兩種流體混合模擬曲線范圍內(nèi),與皂火壕礦床的分布趨勢(shì)大致相同(圖5)?;旌狭黧w中地層水的摩爾分?jǐn)?shù)主要約為0.5~0.9,主要在0.6左右,即地層水的比例約占一半以上,總體上,地層水在混合流體中所占比例較大。熱鹵水與地層水之間全部溶解碳的濃度比P為1.5~5.0,部分樣品在10.0以上。熱鹵水中溶解碳的濃度是地層水中溶解碳的濃度的1.5至5倍左右,部分樣品落在曲線以外,說(shuō)明熱鹵水中溶解碳的濃度是地層水中溶解碳的濃度的10倍以上,熱鹵水中含碳濃度比地層水明顯增高。方解石形成溫度也具有2個(gè)溫度段,分別為60~80℃和90~140℃(圖9)。形成溫度主要在100~120℃范圍內(nèi)。部分方解石膠結(jié)物在較低溫度、地層水比例較高且熱鹵水中碳含量比例高條件下形成,分布在模擬曲線的右下位置,可能是由一期含溶解碳較高的少量熱鹵水與大量的低溫地層水混合后形成方解石。樣品數(shù)據(jù)點(diǎn)在模擬圖上總體分布特征與皂火壕礦床相近,說(shuō)明東勝地區(qū)不同鈾礦床(礦化點(diǎn))成礦類型具有相似性,礦床類型相同。

    4 成因討論

    4.1 流體溫度

    通過(guò)東勝地區(qū)方解石膠結(jié)物C、O同位素組成混合流體模擬,可知方解石膠結(jié)物形成環(huán)境具有明顯的兩個(gè)溫度段:較低溫度段55~80℃和較高溫度段90~140℃。

    肖新建等(2004b)對(duì)東勝地區(qū)鈾礦床礦石中方解石膠結(jié)物的流體包裹體研究表明,所測(cè)均一溫度區(qū)間為58~ 176℃,平均值為 114.9℃,且有兩個(gè)比較明顯的區(qū)間,第一個(gè)區(qū)間主峰在90~100℃,第二個(gè)區(qū)間主峰為 140~150℃。吳柏林等(2006)對(duì)方解石包裹體均一法測(cè)溫,流體溫度主要在40~60℃和100~125℃兩個(gè)溫度段,平均是103℃。通過(guò)模型所得溫度數(shù)據(jù)與實(shí)際測(cè)量的包裹體測(cè)溫?cái)?shù)據(jù)范圍大致對(duì)應(yīng)。

    總體上,流體包裹體溫度范圍表明流體是一種低溫?zé)嵋毫黧w,部分流體溫度相對(duì)較高。包裹體溫度大致分為兩個(gè)溫度段,單憑溫度數(shù)據(jù),表明賦礦地層可能至少經(jīng)歷過(guò)兩期以上的熱流體活動(dòng),或流體比例不同引起了混合流體溫度差異現(xiàn)象??傊摰貐^(qū)鈾成礦作用具有低溫?zé)嵋鹤饔锰卣?,即發(fā)生熱鹵水沿盆地邊緣滲出,在賦礦地層直羅組與地層水滲入流體相互作用而成礦。

    4.2 流體古鹽度與鹽度指數(shù)

    利用皂火壕礦床碳酸鹽方解石膠結(jié)物C、O同位素組成和古鹽度(S)、鹽度指數(shù)(Z)之間的關(guān)系,S=δ18OV-PDB+21.2/0.61,Z=2.048(δ13CV-PDB+50)+0.498(δ18OV-PDB+50)(Keith and Weber,1964),換算出成礦時(shí)流體的古鹽度和鹽度指數(shù),可以探討方解石膠結(jié)物形成時(shí)的流體環(huán)境(湯超等,2016)。全部50件樣品方解石膠結(jié)物流體古鹽度的范圍為20.8%~29.8%,平均值為24.0%。其中,含礦砂巖方解石膠結(jié)物流體古鹽度的范圍為21.8%~26.2%,平均值為23.4%;灰色砂巖方解石膠結(jié)物流體古鹽度的范圍為20.8%~29.8%,平均值為24.1%;綠色砂巖方解石膠結(jié)物流體古鹽度的范圍為21.2%~27.3%,平均值為24.3%(表1)。

    同樣也統(tǒng)計(jì)了全部50件樣品方解石膠結(jié)物流體鹽度指數(shù),其范圍為86.4~120.5,平均值為107.0。其中,含礦砂巖方解石膠結(jié)物流體鹽度指數(shù)范圍為95.6~119.1,平均值為105.8;灰色砂巖方解石膠結(jié)物流體鹽度指數(shù)范圍為96.5~120.5,平均值為107.0;綠色砂巖方解石膠結(jié)物流體鹽度指數(shù)范圍為84.6~119.9,平均值為108.1(表1)。

    李平等(2017)對(duì)包裹體鹽度研究結(jié)果表明,石英次生包裹體鹽度(NaCl)為25.2%~25.3%,平均25.3%。李榮西等(2011)對(duì)砂巖方解石膠結(jié)物包裹體鹽度進(jìn)行了研究,結(jié)果表明包裹體流體鹽度一般為0.8%~15. 0%,其中氣液兩相包裹體的鹽度從0.5%~9.0%,而含NaCl子礦物的包裹體鹽度較高,最高達(dá)27.7%,且具有隨著溫度增高鹽度也相應(yīng)增大趨勢(shì),流體屬于一種較高密度的流體。

    因此,賦礦地層直羅組砂巖方解石膠結(jié)物流體古鹽度整體較高,古鹽度指數(shù)接近120,說(shuō)明該區(qū)鈾成礦形成時(shí)其鹽度與近似湖水咸化環(huán)境相近,有含鹽較高的熱鹵水參與。富含烴、碳的熱鹵水起源于盆地下部, 從盆地中心下部流向盆地邊緣淺部(薛春紀(jì)等,2010)。因此,膠結(jié)物方解石具有深部熱液成因特征。

    4.3 流體來(lái)源與盆地構(gòu)造演化

    中生代以來(lái)三疊紀(jì)末、晚侏羅世和早白堊世末,鄂爾多斯盆地北部經(jīng)歷了三期主要強(qiáng)烈構(gòu)造運(yùn)動(dòng),其中,晚侏羅世在盆地東北部邊緣東勝地區(qū)發(fā)生構(gòu)造變形和隆升作用,早白堊世末(100 ~20 Ma)發(fā)生快速抬升冷卻(張龍,2017)。盆地邊緣地層抬升和地層構(gòu)造變形控制著成礦流體的形成和礦體賦存空間的定位,侏羅紀(jì)—白堊紀(jì)時(shí)期發(fā)育滲出流和滲入流,兩種流體相互混合,礦體定位空間即鈾礦化地層恰好是上侏羅統(tǒng)(薛春紀(jì)等,2010),這是東勝地區(qū)砂巖型鈾礦成礦顯著的特征。

    根據(jù)東勝地區(qū)砂巖型鈾礦賦礦地層直羅組方解石膠結(jié)物C、O同位素組成模型解釋結(jié)果(圖6 ~圖10),在東勝地區(qū)砂巖型鈾礦成礦作用過(guò)程中,地層水滲入流體所占混合的比例較大,一般為2/3以上,吳柏林等(2006)研究認(rèn)為成礦流體主要來(lái)自大氣降水。本次研究根據(jù)大量分析數(shù)據(jù)定量給出地層滲入水在混合流體所占比例范圍,進(jìn)一步深化了對(duì)東勝地區(qū)砂巖型鈾礦成礦理論認(rèn)識(shí)。

    從熱鹵水與地層水之間全部溶解碳的濃度比P的變化,可以劃分出P的三個(gè)主要分布范圍:1.0~2.0、2.0~5.0、5.0~10.0,P變化范圍較大,說(shuō)明熱鹵水流體來(lái)源具有較大差異,推測(cè)低溫?zé)嵋夯顒?dòng)有多個(gè)期次,至少在三個(gè)期次以上。

    5 結(jié)論

    通過(guò)鄂爾多斯盆地東北部東勝地區(qū)砂巖型鈾礦賦礦地層直羅組含礦砂巖、灰色砂巖及綠色砂巖中方解石膠結(jié)物C、O同位素組成分析,建立了滲入地層水和滲出熱鹵水兩種流體混合C、O同位素組成定量模型,并根據(jù)該模型反演兩種流體中溶解碳的濃度比、流體比例以及熱液溫度等成礦因素,得到以下主要結(jié)論:

    (1)方解石膠結(jié)物C、O同位素組成特征表明,除部分樣品方解石膠結(jié)物是地層沉積形成外,大部分方解石膠結(jié)物形成與有機(jī)酸脫羧基作用相關(guān)。其形成的主要過(guò)程為,盆地深部巖石圈熱活動(dòng)增強(qiáng),古生界有機(jī)酸脫羧基作用形成CO2,晚侏羅世至早白堊世末盆地邊緣抬升和盆地地層變形,伴隨構(gòu)造運(yùn)動(dòng)富含CO2的熱鹵水與地層水兩種流體混合后形成含較高鹽度的成礦流體,伴隨有大量的方解石膠結(jié)物沉淀。

    (2)根據(jù)地層水和地下熱鹵水C、O同位素組成已有基礎(chǔ)資料,假定成礦作用發(fā)生時(shí)滲入地層水和滲出熱鹵水兩種流體C、O同位素組成參數(shù),用賦礦地層方解石膠結(jié)物C、O同位素組成建立了東勝地區(qū)砂巖型鈾礦模型,可定量模擬出方解石形成時(shí)溫度、兩種流體溶解碳濃度比及流體比例等流體成礦條件。

    (3)東勝地區(qū)砂巖型鈾礦成礦流體地層水與熱鹵水比例為0.5~0.9,地層水所占混合流體的比例較大,一般在2/3以上;熱鹵水與地層水溶解碳濃度比為主要范圍為1.5~5.0,部分大于10.0,熱鹵水中溶解碳濃度較高,是富含CO2的流體?;旌狭黧w溫度分為兩個(gè)主要范圍55~80℃、90~140℃,結(jié)合古鹽度和鹽度指數(shù)反演,說(shuō)明成礦流體為有機(jī)與無(wú)機(jī)混合成因的低溫?zé)嵋毫黧w。

    致謝:對(duì)核工業(yè)二〇八大隊(duì)彭云彪研究員、苗愛(ài)生研究員、楊建新研究員在野外工作提供的熱情幫助以及審稿專家和責(zé)任編輯在文章修改方面付出的辛勞表示誠(chéng)摯的感謝。

    參 考 文 獻(xiàn)/References

    (The literature whose publishing year followed by a “&” is in Chinese with English abstract; The literature whose publishing year followed by a “#” is in Chinese without English abstract)

    陳路路, 馮曉曦, 司馬獻(xiàn)章, 李建國(guó), 郭虎, 陳印, 趙華雷, 湯超, 王貴, 劉忠仁, 李曙光. 2017. 鄂爾多斯盆地納嶺溝地區(qū)鈾礦物賦存形式研究及其地質(zhì)意義. 地質(zhì)與勘探, 53(4): 632~642.

    丁悌平. 2002.穩(wěn)定同位素測(cè)試技術(shù)與參考物質(zhì)研究現(xiàn)狀及發(fā)展趨勢(shì). 巖礦測(cè)試, 21(4): 291~300.

    馮喬, 秦宇, 付鎖堂, 柳益群, 周鼎武. 2016.東勝鈾礦砂巖中方解石富集及鈾礦成因. 高校地質(zhì)學(xué)報(bào), 22(1): 53~59.

    焦養(yǎng)泉, 陳安平, 王敏芳, 吳立群, 原海濤, 楊琴, 張承澤, 徐志誠(chéng). 2005.鄂爾多斯盆地東北部直羅組底部砂體成因分析——砂巖型鈾礦床預(yù)測(cè)的空間定位基礎(chǔ). 沉積學(xué)報(bào), 23(3): 371~379.

    李宏濤, 蔡春芳, 羅曉容, 李開(kāi)開(kāi). 2007. 鄂爾多斯北部直羅組中烴類包裹體地球化學(xué)特征及來(lái)源分析.沉積學(xué)報(bào), 25(3):467~473.

    李平, 劉紅旭, 張字龍, 丁波, 宋昊. 2017. 鄂爾多斯盆地納嶺溝鈾礦流體包裹體研究. 第八屆全國(guó)成礦理論與找礦方法學(xué)術(shù)討論會(huì)論文摘要文集, 483.

    李榮西, 段立志, 陳寶赟, 張少妮. 2011. 東勝砂巖型鈾礦氧化酸性流體與還原堿性熱液流體過(guò)渡界面蝕變帶成礦作用研究. 大地構(gòu)造與成礦學(xué), 35(4): 525~532.

    李延河, 段超, 趙悅, 裴浩翔, 任順利. 2016. 氧化還原障在熱液鈾礦成礦中的作用. 地質(zhì)學(xué)報(bào), 90(2): 201~218.

    李子穎, 陳安平, 方錫珩, 歐光習(xí), 張珂, 焦養(yǎng)泉, 夏毓亮, 陳法正, 周文斌, 劉忠厚, 吳仁貴, 肖新建, 孫曄. 2006. 鄂爾多斯盆地東北部砂巖型鈾礦成礦機(jī)理和疊合成礦模式. 礦床地質(zhì), 25(增刊1): 245~248.

    李子穎, 方錫珩, 陳安平, 歐光習(xí), 肖新建, 孫曄, 劉池洋, 王毅. 2007. 鄂爾多斯盆地北部砂巖型鈾礦目標(biāo)層灰綠色砂巖成因. 中國(guó)科學(xué), 37(增刊1): 139~146.

    李子穎, 方錫衍, 秦明寬, 陳安平, 蔡煜琦, 夏毓亮, 孫曄, 張珂, 焦養(yǎng)泉, 歐光習(xí), 陳法正. 2019. 鄂爾多斯盆地北部砂巖鈾成礦作用. 北京:地質(zhì)出版社:175~210.

    劉漢彬, 夏毓亮, 田時(shí)豐. 2007. 東勝地區(qū)砂巖型鈾礦成礦年代學(xué)及成礦鈾源研究. 鈾礦地質(zhì), 23(1): 23~29.

    劉建明, 劉家軍, 顧雪祥. 1997. 沉積盆地中的流體活動(dòng)及其成礦作用. 巖石礦物學(xué)雜志, 16(4): 341~352.

    彭云彪, 陳安平, 楊建新. 2017. 砂巖型鈾礦理論創(chuàng)新與找礦新突破. 中國(guó)核工業(yè), 11: 40~41.

    秦明寬, 李子穎, 劉章月. 2017. 鈾資源勘查技術(shù)新深度. 中國(guó)核工業(yè), 11: 32~34.

    湯超, 司馬獻(xiàn)章, 朱強(qiáng), 馮曉曦, 陳印, 陳路路, 劉曉雪. 2016. 東勝地區(qū)直羅組含鈾砂巖中方解石的碳氧同位素組成及鈾成礦意義. 高校地質(zhì)學(xué)報(bào), 22(4): 698~706.

    王大銳. 2000. 油氣穩(wěn)定同位素地球化學(xué). 北京: 石油工業(yè)出版社:1~295.

    吳柏林, 劉池陽(yáng), 張復(fù)新, 方錫衍, 劉雄. 2006. 東勝砂巖型鈾礦后生蝕變地球化學(xué)性質(zhì)及其成礦意義. 地質(zhì)學(xué)報(bào), 80(5): 740~747.

    吳柏林, 魏安軍, 胡亮, 寸小妮, 孫莉, 羅晶晶, 宋子升, 張龍, 程相虎, 張婉瑩, 王建強(qiáng), 張東東. 2016.內(nèi)蒙古東勝鈾礦區(qū)后生蝕變的穩(wěn)定同位素特征及其地質(zhì)意義. 地質(zhì)通報(bào), 35(12): 2133~2145.

    吳春勇, 蘇小四, 郭金淼, 董維紅. 2011. 鄂爾多斯沙漠高原白堊系地下水水化學(xué)演化的多元統(tǒng)計(jì)分析. 世界地質(zhì), (2): 244~253.

    夏毓亮, 劉漢彬. 2005. 鄂爾多斯盆地東勝地區(qū)直羅組砂體鈾的預(yù)富集與鈾成礦. 世界核地質(zhì)科學(xué), 22(4): 187~191.

    肖新建, 李子穎, 陳安平. 2004a. 東勝地區(qū)砂巖型鈾礦床后生蝕變礦物分帶特征初步研究.鈾礦地質(zhì), 20(3): 136~140.

    肖新建, 李子穎, 方錫珩, 歐光習(xí), 孫曄, 陳安平. 2004b. 東勝砂巖型鈾礦床低溫?zé)嵋毫黧w的證據(jù)及意義.礦物巖石地球化學(xué)通報(bào), 23(4): 301~304.

    徐浩. 2017. 鄂爾多斯盆地煤系礦產(chǎn)資源賦存規(guī)律的構(gòu)造控制研究. 導(dǎo)師:曹代勇. 北京: 中國(guó)礦業(yè)大學(xué)博士學(xué)位論文: 1~130.

    薛春紀(jì), 池國(guó)祥, 薛偉. 2010. 鄂爾多斯盆地砂巖型鈾成礦中兩種流體系統(tǒng)相互作用——地球化學(xué)證據(jù)和流體動(dòng)力學(xué)模擬. 礦床地質(zhì), 29(1): 134~142.

    楊曉勇, 凌明星, 賴小東. 2009. 鄂爾多斯盆地東勝地區(qū)地浸砂巖型鈾礦成礦模型. 地學(xué)前緣, 16(2): 239~249.

    易超, 王貴, 李西得, 張康, 王永君. 2018. 鄂爾多斯盆地東北部直羅組鈾富集特征及鈾成礦模式探討. 礦床地質(zhì), 37(4): 835~852.

    張賓, 劉紅旭, 丁波, 易超, 張艷. 2020. 鄂爾多斯盆地納嶺溝鈾礦床巖石學(xué)、礦物學(xué)特征及其對(duì)鈾成礦作用的指示意義. 地質(zhì)論評(píng), 66(2): 410~424.

    張金帶. 2016. 我國(guó)砂巖型鈾礦成礦理論的創(chuàng)新和發(fā)展. 鈾礦地質(zhì), 32(6): 321~332.

    張龍. 2017. 鄂爾多斯盆地北部天然氣逸散與鈾成礦效應(yīng).導(dǎo)師:劉池陽(yáng).西安: 西北大學(xué)博士論文: 1~167.

    張祖還, 趙懿英, 章邦桐, 沈渭洲.1984. 鈾地球化學(xué).北京:原子能出版社: 160~166.

    鄭永飛, 陳江峰. 2000. 穩(wěn)定同位素地球化學(xué).北京: 科學(xué)出版社: 281~313.

    Bath A H, Milodowski A E, Piro B. 1987. Diagenesis of Wessex and East Yorkshire—Lincolnshire basins, UK: a stable isotope study. Geological Society London Special Publications, 36(1): 173~190.

    Chen Lulu, Feng Xiaoxi, Sima Xianzhang, Li Jianguo, Guo Hu, Chen Yin, Zhao Hualei, Tang Chao, Wang Gui, Liu Zhongren, Li Shuguang. 2017&. Occurrence forms of the uranium minerals in the Nalinggou are of the Ordos Basin and geological implications. Geology and Exploration, 53(4): 632~642.

    Ding Tiping. 2002&. Present status and prospect of analytical techniques and reference materials for stable isotopes. Rock and Mineral Analysis, 21(4): 291~300.

    Feng Qiao, Qin Yu, Fu Suotang, Liu Yiqun, Zhou Dingwu.2016&. Enrichment of calcite of the genesis of uranium deposits in Dongsheng uranium sandstone. Geological Journal of China Universities, 22(1): 53~59.

    Jiao Yangquan, Chen Anping, Wang Minfang, Wu Liqun, Yuan Haitao, Yang Qin, Zhang Chengze, Xu Zhicheng. 2005&. Genetic analysis of the bottom sandstone of Zhiluo Formation, Northeastern Ordos Basin: Predictive base of spatial orientation of sandstone-type uranium deposit. Acta Sedimentologica Sinica, 23(3): 371~378.

    Keith M L, Weber J N. 1964. Carbon and oxygen isotopic composition of selected limestone and fossils. Geochimica et Cosmochimica Acta, 28(10~11): 1787~1816.

    Li Hongtao,Cai Chunfang, Luo Xiaorong, Li Kaikai. 2007&. Geochemical characteristics and source of hydrocarbon inclusions in Zhiluo Formation sandstone in Dongsheng area, northern Ordos Basin. Acta Sedimentologica Sinica, 25(3): 467~473.

    Li Ping, Liu Hongxu, Zhang Zilong, Ding Bo, Song Hao.2017#. Study on fluid inclusions of Nalinggou uranium deposit in Ordos Basin. Abstracts of the 8th National Symposium on Metallogenic Theory and Prospecting Methods: 483.

    Li Rongxi, He Ying, Li Jinbao, Li Jihong, Li Xin. 2006&. Isotope composition of fluid inclusion and tracing to the source of minelization fluid in Dongsheng uranium deposit, north Ordos Basin. Acta Geologica Sinica, 80(5): 753~760.

    Li Rongxi,Duan Lizhi, Chen Baoyun, Zhang Shaoni.2011&.Alteration and metallogeny on the oxic acid/anoxic alkali interface of the Dongsheng uranium deposit in northern Ordos Basin. Geotectonica et Metallogenia,35(4):525~532.

    Li Yanhe, Duan Chao, Zhao Yue, Pei Haoxiang, Ren Shunli.2016&. The role of oxidizing reducing barrier in mineralization of hydrothermal uranium ore. Acta Geologica Sinica, 90(2): 201~218.

    Li Ziying, Chen Anping, Fang Xiheng, Ou Guangxi, Zhang Ke, Jiao Yangquan, Xia Yuliang, Chen Fazheng, Zhou Wenbin, Liu Zhonghou, Wu Rengui, Xiao Xinjian, Sun Ye. 2006&. The metallogenic mechanism and superimposed metallogenic model of sandstone-type uranium ore in the northeast of Ordos Basin. Geological Deposits, 25(Supp. 1): 245~248.

    Li Ziying, Fang Xiheng, Chen Anping, Ou Guangxi, Xiao Xinjian, Sun Ye, Liu Chiyang, Wang Yi. 2007#. Formation of gray-green sandstone in target strata of sandstone-type uranium deposits in northern Ordos Basin. Science China, 37(Supp. 1): 139~146.

    Li Ziying, Fang Xiheng, Qin Mingkuan, Chen Anping, Cai Yuqi, Xia Yuliang, Sun Ye, Zhang Ke, Jiao Yangquan, Ou Guangxi, Chen Fazheng. 2019&. Sandstone Uranium Mineralization in Northern Ordos Basin. Bejing: Geological Press: 175~210.

    Liu Hanbin, Xia Yuliang, Tian Shifeng. 2007&. Study on geochronology and uranium source of sandstone-type uranium deposit in Dongsheng area. Uranium Geology, 23(1): 23~29.

    Liu Jianming, Liu Jiajun, GuXuexiang. 1997&.Basin fluids and their related ore deposits. Acta Petrologica et Mineralogica, 16(4): 341~352.

    Ohmoto and Rye R O. 1979. Isotopes of sulfur and carbon. In: Barnes H L. ed. Geochemistry of hydrothermal ore deposits (2nd Edition). New York:John Wiley and Sons: 509~567.

    Peng Yunbiao, Chen Anping, Yang Jianxin. 2017&. Theoretical innovation and prospecting breakthrough of sandstone-type uranium deposits. China Nuclear Industry, 11: 40~41.

    Qin Mingkuan, Li Ziying, Liu Zhangyue.2017&. New depth of uranium exploration technology. China Nuclear Industry, 11: 32~34.

    Schmidt M, Xeflide S, Botz R, Mann S. 2005. Oxygen isotope fractionation during synthesis of Ca—Mg-carbonate and implications for sedimentary dolomite formation. Geochimica et Cosmochimica Acta, 69(1): 4665~4674.

    Tang Chao, Sima Xianzhang, Zhu Qiang, Feng Xiaoxi, Chen Yin, Chen Lulu, Liu Xiaoxue. 2016&. Carbon and oxygen isotopic composition and uranium mineralization significance of calcite of Zhiluo Formation uranium-bearing sandstones in Dongsheng area. Geological Journal of China Universities, 22(4): 698~706.

    Wang Darui. 2000&.Stable Isotope Geochemistry of Oil and Gas. Beijing: Petroleum Industry Press: 1~295.

    Wu Bolin, LiuChiyang, Zhang Fuxin, Fang Xiheng, Liu Xiong. 2006&. Geochemical characteristics of epigenetic alternation in Dongsheng sandstone-type uranium deposit and its metallogenic signification. Acta Geologica Sinica, 80(5): 740~747.

    Wu Bolin, Wei Anjun, Hu Liang, Cun Xiaoni, Sun Li, Luo Jingjing, Song Zishen, Zhang Long, Chen Xianghu, Zhang Wanying, Wang Jianqiang, Zhang Dongdong. 2016&. The epigenetic alternation stable isotope characteristics of the Dongsheng uranium ore district and their geological implications. Geological Bulletin of China, 35(12): 2133~2145.

    Wu Chunyong, Su Xiaosi, Guo Jinmiao, Dong Weihong. 2011&. Multivariate statistical analysis of hydrogeochemical evolution of groundwater in Cretaceous aquifer Ordos desert plateau. World Geology, (2): 244~253.

    Xia Yuliang, Liu Hanbin. 2005&.Pre-enrichment and metallogeny of uranium in Zhiluo Formation sand bodies of Dongsheng area, Ordos Basin. World Nuclear Geoscience, 22(4): 187~191.

    Xiao Xinjian, Li Ziying, Chen Anping. 2004a&. Preliminary study on features of mineralogical zoning of epigenetic alteration at sandstone-type uranium deposit, Dongsheng area, Ordos basin. Uranium Geology, 20(3): 136~140.

    Xiao Xinjian, Li Ziying, Fang Xiheng, Ou Guangxi, Sun Ye, Chen Anping. 2004b&. The evidences and significances of epithermal mineralization fluid in the Dongsheng sandstone type uranium deposit. Bulletin of Mineralogy, Petrology and Geochemistry,23(4): 301~304.

    Xu Hao. 2017&. Occurrence Characteristics and Tectonic Controls of Coal Series Mineral Resources in Ordos Basin. Tutor: Cao Daiyong. Ph. D. Thesis of China University of Mining and Technology: 1~130.

    Xue Chunji, Chi Guoxiang, Xue Wei. 2010&. Interaction of two fluid systems in formation of sandstone-type uranium deposits in Ordos Basin: Geochemical evidence and hydrodynamic modeling. Mineral Deposits, 29(1): 134~142.

    Yang Xiaoyong, Ling Mingxing, Lai Xiaodong. 2009&. Metallogenic model of the Dongsheng in-situ leaching sandstone-type uranium deposit in the Ordos Basin. Earth Science Frontiers, 16(2): 239~249.

    Yi Chao, WangGui, Li Xide, Zhang Kang, Wang Yongjun. 2018&. A tentative discussion on uranium enrichment characteristics and metallogenic model in Zhiluo Formation, northeastern Ordos Basin. Mineral Deposits, 37(4): 835~852.

    Zhang Bin, Liu Hongxu,Ding Bo, Yi Chao, Zhang Yan. 2020&. Petrological and mineralogical characteristics of Nalinggou uranium deposit in Ordos Basin and its indicative significance for uranium mineralization. Geological Review, 66(2): 410~424.

    Zhang Jindai. 2016&. Innovation anddevelopment of metallogennic theory for sandstone type uranium deposit in China. Uranium Geology, 32(6): 321~332.

    Zhang Long. 2017&.Natural Gas Dissipation and Uranium Mineralization in the Northern Ordos Basin. Tutor: Liuchiyang. Ph. D. Thesis of Northwest University: 1~167.

    Zhang Zuhuan, Zhao Yiying, Zhang Bangtong, Shen Weizhou. 1987#. Uranium Geochemistry. Beijing: Atomic Energy Press: 160~166.

    Zheng Yongfei, Chen Jiangfeng&. 2000.Stable Isotope Geochemistry. Beijing: Science Press: 281~313.

    猜你喜歡
    直羅東勝方解石
    巨厚充水含水層疏放水層位優(yōu)化研究
    鄂爾多斯盆地北部直羅組含水層研究進(jìn)展與水害防治建議
    F-在方解石表面的吸附及其對(duì)方解石表面性質(zhì)的影響
    明初東勝衛(wèi)置廢考辨
    氯化鈣和碳酸鈉對(duì)方解石浮選的影響及其機(jī)理研究
    貴州重晶石與方解石常溫浮選分離試驗(yàn)研究
    陸相頁(yè)巖氣資源評(píng)價(jià)初探:以延長(zhǎng)直羅——下寺灣區(qū)中生界長(zhǎng)7段為例
    直羅油田三疊系延長(zhǎng)組長(zhǎng)8段油氣成藏條件分析
    螢石與方解石、重晶石等鹽類礦物浮選分離現(xiàn)狀
    杭錦旗地區(qū)東勝氣田構(gòu)造及斷裂特征研究
    河南科技(2014年7期)2014-02-27 14:11:29
    亚洲精华国产精华液的使用体验| 久久影院123| 春色校园在线视频观看| 乱人伦中国视频| 中国美白少妇内射xxxbb| 人成视频在线观看免费观看| 欧美国产精品一级二级三级| 少妇 在线观看| 日本91视频免费播放| 亚洲经典国产精华液单| 久久久久精品久久久久真实原创| 91精品伊人久久大香线蕉| 亚洲精品日本国产第一区| 欧美精品亚洲一区二区| av国产久精品久网站免费入址| 久久午夜福利片| 夜夜骑夜夜射夜夜干| 亚洲精品日本国产第一区| 在线免费观看不下载黄p国产| 男人添女人高潮全过程视频| 99视频精品全部免费 在线| 久久久久精品久久久久真实原创| 久久午夜福利片| 十分钟在线观看高清视频www| 久久韩国三级中文字幕| 美女视频免费永久观看网站| 99国产精品免费福利视频| 只有这里有精品99| 精品少妇久久久久久888优播| 国产日韩欧美视频二区| 日韩视频在线欧美| 狠狠婷婷综合久久久久久88av| 男的添女的下面高潮视频| 婷婷色综合大香蕉| 草草在线视频免费看| 欧美国产精品一级二级三级| 黄网站色视频无遮挡免费观看| av视频免费观看在线观看| 午夜福利网站1000一区二区三区| 久久人人爽人人爽人人片va| 97精品久久久久久久久久精品| 日韩一区二区视频免费看| 男女免费视频国产| 久久久久人妻精品一区果冻| a级毛色黄片| 蜜桃在线观看..| 欧美精品人与动牲交sv欧美| 又粗又硬又长又爽又黄的视频| av不卡在线播放| 男女午夜视频在线观看 | 欧美少妇被猛烈插入视频| 春色校园在线视频观看| 日本欧美视频一区| 在线看a的网站| 99热网站在线观看| 婷婷色综合www| 国产乱人偷精品视频| 日韩,欧美,国产一区二区三区| 免费高清在线观看视频在线观看| 久久热在线av| 国产乱来视频区| 精品一区二区三卡| 伊人亚洲综合成人网| 亚洲国产欧美在线一区| 亚洲国产色片| 少妇的逼水好多| 纵有疾风起免费观看全集完整版| 激情五月婷婷亚洲| 久久人人爽人人片av| 国产又爽黄色视频| 亚洲精品aⅴ在线观看| 成人国语在线视频| www日本在线高清视频| 涩涩av久久男人的天堂| 精品熟女少妇av免费看| 高清av免费在线| av不卡在线播放| 91成人精品电影| 黑丝袜美女国产一区| 亚洲欧美成人综合另类久久久| 亚洲,一卡二卡三卡| 丰满饥渴人妻一区二区三| 少妇被粗大的猛进出69影院 | 欧美国产精品一级二级三级| 免费黄色在线免费观看| 久久久久精品人妻al黑| 街头女战士在线观看网站| 久久ye,这里只有精品| 免费黄色在线免费观看| 精品少妇黑人巨大在线播放| 人人妻人人澡人人看| 热re99久久国产66热| av黄色大香蕉| 亚洲av日韩在线播放| 九色成人免费人妻av| 少妇人妻精品综合一区二区| 亚洲国产日韩一区二区| 日韩精品有码人妻一区| 老司机亚洲免费影院| 午夜福利乱码中文字幕| 日本av免费视频播放| 最近中文字幕高清免费大全6| 久久久久久久国产电影| 成年美女黄网站色视频大全免费| 9色porny在线观看| 久久久久久久国产电影| 亚洲精品aⅴ在线观看| 伦理电影免费视频| 天天躁夜夜躁狠狠躁躁| 满18在线观看网站| 丰满迷人的少妇在线观看| 亚洲精品中文字幕在线视频| 国产有黄有色有爽视频| 国产视频首页在线观看| 乱人伦中国视频| 寂寞人妻少妇视频99o| 免费av中文字幕在线| 国产乱来视频区| 国产一区二区三区综合在线观看 | 男人添女人高潮全过程视频| 波多野结衣一区麻豆| 99久久综合免费| tube8黄色片| 亚洲一码二码三码区别大吗| 天堂8中文在线网| 天天躁夜夜躁狠狠躁躁| 男女下面插进去视频免费观看 | 久久久国产精品麻豆| 国产精品秋霞免费鲁丝片| 我要看黄色一级片免费的| 精品一区二区三卡| 欧美日韩精品成人综合77777| 久久午夜福利片| 亚洲国产看品久久| 99香蕉大伊视频| 午夜老司机福利剧场| 国产色婷婷99| 中文欧美无线码| 美女福利国产在线| 亚洲人成网站在线观看播放| 亚洲国产av新网站| 最近的中文字幕免费完整| 国产探花极品一区二区| 自拍欧美九色日韩亚洲蝌蚪91| 99热全是精品| 性高湖久久久久久久久免费观看| 免费久久久久久久精品成人欧美视频 | 99久久人妻综合| 日本免费在线观看一区| 欧美性感艳星| 亚洲av在线观看美女高潮| 亚洲欧美一区二区三区黑人 | 三级国产精品片| 大香蕉久久网| 九色成人免费人妻av| 97人妻天天添夜夜摸| 国产精品久久久av美女十八| 国产精品国产三级国产av玫瑰| 自拍欧美九色日韩亚洲蝌蚪91| 精品卡一卡二卡四卡免费| 久久精品国产亚洲av天美| 午夜福利,免费看| 午夜日本视频在线| 一本色道久久久久久精品综合| 欧美日韩视频高清一区二区三区二| 中文字幕精品免费在线观看视频 | 亚洲激情五月婷婷啪啪| 一本色道久久久久久精品综合| 久久 成人 亚洲| 国产69精品久久久久777片| 人妻 亚洲 视频| 天堂中文最新版在线下载| 精品国产乱码久久久久久小说| 黄片播放在线免费| 国产一区有黄有色的免费视频| 丝袜人妻中文字幕| 欧美人与性动交α欧美软件 | 美女福利国产在线| 欧美老熟妇乱子伦牲交| av又黄又爽大尺度在线免费看| 亚洲精品成人av观看孕妇| 看免费成人av毛片| 伦理电影免费视频| 乱码一卡2卡4卡精品| freevideosex欧美| av电影中文网址| 国产成人午夜福利电影在线观看| 一级毛片 在线播放| 久久久久国产精品人妻一区二区| 成人亚洲精品一区在线观看| 中文字幕最新亚洲高清| 一级爰片在线观看| 最黄视频免费看| 国产精品久久久久久av不卡| 色5月婷婷丁香| 美女国产视频在线观看| 黄色 视频免费看| 精品熟女少妇av免费看| 午夜福利视频精品| 国产精品不卡视频一区二区| 大码成人一级视频| 如何舔出高潮| 日韩av在线免费看完整版不卡| 亚洲精品av麻豆狂野| 天天影视国产精品| 视频区图区小说| 亚洲国产色片| 五月开心婷婷网| 国产男人的电影天堂91| 亚洲成人一二三区av| 免费观看在线日韩| 欧美性感艳星| 两个人免费观看高清视频| 亚洲天堂av无毛| 国产老妇伦熟女老妇高清| 十八禁网站网址无遮挡| 男女国产视频网站| 久久精品国产鲁丝片午夜精品| 日韩av在线免费看完整版不卡| 男男h啪啪无遮挡| 成人午夜精彩视频在线观看| 97精品久久久久久久久久精品| 日韩视频在线欧美| 超色免费av| 一级毛片我不卡| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 日韩大片免费观看网站| 黄色 视频免费看| 欧美日韩一区二区视频在线观看视频在线| 18禁裸乳无遮挡动漫免费视频| 久久久亚洲精品成人影院| 国产探花极品一区二区| 色吧在线观看| 丝袜脚勾引网站| 精品人妻熟女毛片av久久网站| 蜜臀久久99精品久久宅男| 国产欧美另类精品又又久久亚洲欧美| 精品少妇内射三级| 青春草视频在线免费观看| 欧美精品一区二区大全| 9色porny在线观看| 乱码一卡2卡4卡精品| 男女边摸边吃奶| 亚洲欧美日韩卡通动漫| 夜夜爽夜夜爽视频| av播播在线观看一区| 热re99久久国产66热| 国产淫语在线视频| 这个男人来自地球电影免费观看 | 国产欧美亚洲国产| 91国产中文字幕| 亚洲国产精品一区三区| 女人被躁到高潮嗷嗷叫费观| 国产在线免费精品| 精品一区二区免费观看| 一本久久精品| 人妻 亚洲 视频| 国产成人免费观看mmmm| 中文字幕制服av| 91午夜精品亚洲一区二区三区| 国产免费一区二区三区四区乱码| 国产在线免费精品| 久久99热这里只频精品6学生| 免费人成在线观看视频色| 精品久久久久久电影网| 久久久久久人妻| 日产精品乱码卡一卡2卡三| 热99国产精品久久久久久7| 欧美激情极品国产一区二区三区 | av有码第一页| 成人18禁高潮啪啪吃奶动态图| 巨乳人妻的诱惑在线观看| 又大又黄又爽视频免费| 亚洲精品中文字幕在线视频| 大香蕉久久成人网| 少妇 在线观看| 永久网站在线| 人妻人人澡人人爽人人| 成人亚洲欧美一区二区av| 综合色丁香网| 国产精品不卡视频一区二区| 国产片内射在线| 久久国产精品男人的天堂亚洲 | 一区二区日韩欧美中文字幕 | 亚洲美女搞黄在线观看| 国产熟女午夜一区二区三区| 最近的中文字幕免费完整| 久久人妻熟女aⅴ| 91精品伊人久久大香线蕉| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲精品国产色婷婷电影| 老女人水多毛片| 亚洲av成人精品一二三区| 欧美亚洲日本最大视频资源| 亚洲精品,欧美精品| 欧美精品亚洲一区二区| 22中文网久久字幕| 人妻少妇偷人精品九色| 纵有疾风起免费观看全集完整版| 五月天丁香电影| 亚洲成人av在线免费| 久久99热这里只频精品6学生| 久热这里只有精品99| 少妇的丰满在线观看| 精品亚洲乱码少妇综合久久| 久久这里只有精品19| 飞空精品影院首页| 国产精品成人在线| 国产综合精华液| 黄网站色视频无遮挡免费观看| 欧美国产精品va在线观看不卡| 国产爽快片一区二区三区| 欧美激情极品国产一区二区三区 | 毛片一级片免费看久久久久| videos熟女内射| 男的添女的下面高潮视频| 综合色丁香网| 国产日韩欧美亚洲二区| 91午夜精品亚洲一区二区三区| 啦啦啦中文免费视频观看日本| 极品少妇高潮喷水抽搐| 久久久久久久久久久久大奶| 久久青草综合色| 午夜福利视频在线观看免费| 国产极品天堂在线| 美女福利国产在线| 亚洲丝袜综合中文字幕| 哪个播放器可以免费观看大片| 咕卡用的链子| 久久毛片免费看一区二区三区| 日本wwww免费看| 国产精品一区www在线观看| 亚洲久久久国产精品| 在线观看三级黄色| 久久狼人影院| 日本猛色少妇xxxxx猛交久久| 免费少妇av软件| 七月丁香在线播放| 国产亚洲精品第一综合不卡 | 精品亚洲乱码少妇综合久久| av片东京热男人的天堂| 丁香六月天网| 日韩中文字幕视频在线看片| 久久狼人影院| 久久久精品94久久精品| 亚洲精品一二三| 男女国产视频网站| 日韩成人伦理影院| 久久热在线av| 国产在线视频一区二区| 在线观看免费高清a一片| 日本黄色日本黄色录像| 啦啦啦中文免费视频观看日本| 91精品三级在线观看| 日韩成人av中文字幕在线观看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 9色porny在线观看| 高清视频免费观看一区二区| 久久精品久久精品一区二区三区| 国产在线免费精品| 9191精品国产免费久久| 免费观看在线日韩| 久久久久精品性色| 久久久亚洲精品成人影院| 国产黄色视频一区二区在线观看| 亚洲精品久久久久久婷婷小说| 韩国av在线不卡| 视频在线观看一区二区三区| 性色av一级| 香蕉国产在线看| 侵犯人妻中文字幕一二三四区| 一级,二级,三级黄色视频| 久久精品夜色国产| 成人手机av| 妹子高潮喷水视频| av免费观看日本| 久久久久久久大尺度免费视频| 69精品国产乱码久久久| 内地一区二区视频在线| 色吧在线观看| 久久精品熟女亚洲av麻豆精品| 欧美精品国产亚洲| 国产午夜精品一二区理论片| 久久午夜福利片| 啦啦啦视频在线资源免费观看| 亚洲第一区二区三区不卡| 久久青草综合色| 成年女人在线观看亚洲视频| 日韩不卡一区二区三区视频在线| 免费不卡的大黄色大毛片视频在线观看| 国产精品久久久久久久久免| av播播在线观看一区| 亚洲四区av| 五月开心婷婷网| 亚洲欧美日韩卡通动漫| 久久精品国产亚洲av涩爱| 日韩av不卡免费在线播放| a级毛色黄片| 亚洲av日韩在线播放| 久久青草综合色| 9热在线视频观看99| 亚洲成人av在线免费| videossex国产| 亚洲精品美女久久av网站| 久久久精品94久久精品| 亚洲欧美成人精品一区二区| 天堂8中文在线网| 国产伦理片在线播放av一区| 国产精品一区二区在线不卡| 国产av一区二区精品久久| 又粗又硬又长又爽又黄的视频| 欧美精品高潮呻吟av久久| 一级毛片 在线播放| 免费黄频网站在线观看国产| 亚洲综合色网址| 国产淫语在线视频| 亚洲欧洲精品一区二区精品久久久 | 亚洲精品乱码久久久久久按摩| 五月天丁香电影| 午夜福利网站1000一区二区三区| 最近的中文字幕免费完整| 国产黄色免费在线视频| 另类精品久久| 免费播放大片免费观看视频在线观看| 精品亚洲乱码少妇综合久久| 国产一区亚洲一区在线观看| 久热这里只有精品99| 国产片内射在线| 久久久欧美国产精品| 在线观看国产h片| 欧美人与善性xxx| 久久精品久久精品一区二区三区| 久久精品人人爽人人爽视色| 国产在视频线精品| 成人毛片a级毛片在线播放| 亚洲欧美一区二区三区国产| 国产黄频视频在线观看| 中文字幕制服av| 99热国产这里只有精品6| 日韩制服丝袜自拍偷拍| 各种免费的搞黄视频| 青春草视频在线免费观看| 97在线视频观看| 亚洲,欧美精品.| 亚洲精品久久午夜乱码| 亚洲av电影在线进入| 好男人视频免费观看在线| 你懂的网址亚洲精品在线观看| 欧美日韩av久久| 国产精品国产三级专区第一集| 另类精品久久| 久久久久久伊人网av| 天天影视国产精品| 中文字幕人妻熟女乱码| 熟妇人妻不卡中文字幕| 一本色道久久久久久精品综合| 国产精品久久久av美女十八| 一级毛片电影观看| 久久精品国产亚洲av涩爱| 亚洲高清免费不卡视频| 亚洲av日韩在线播放| 欧美成人午夜精品| 免费在线观看黄色视频的| 国产深夜福利视频在线观看| 在线观看三级黄色| 免费看不卡的av| 免费观看无遮挡的男女| 亚洲欧美一区二区三区黑人 | 80岁老熟妇乱子伦牲交| 观看美女的网站| 久久久国产精品麻豆| 最近手机中文字幕大全| 一级,二级,三级黄色视频| 免费看av在线观看网站| 菩萨蛮人人尽说江南好唐韦庄| 乱人伦中国视频| 久久久久精品久久久久真实原创| 国产精品国产三级国产专区5o| 亚洲欧美一区二区三区黑人 | 大香蕉97超碰在线| 免费av不卡在线播放| 精品少妇久久久久久888优播| 人妻 亚洲 视频| av又黄又爽大尺度在线免费看| 中文字幕制服av| 日韩制服骚丝袜av| 男女无遮挡免费网站观看| 久久人人爽av亚洲精品天堂| 亚洲人成77777在线视频| 在现免费观看毛片| 久久99精品国语久久久| 十八禁高潮呻吟视频| 国产国拍精品亚洲av在线观看| 18禁国产床啪视频网站| 毛片一级片免费看久久久久| 国产成人精品久久久久久| 丰满饥渴人妻一区二区三| 日韩不卡一区二区三区视频在线| 极品少妇高潮喷水抽搐| 国产一区二区在线观看日韩| 久久99热6这里只有精品| a级毛片在线看网站| 精品久久久精品久久久| av有码第一页| 精品少妇久久久久久888优播| 精品国产露脸久久av麻豆| 国产精品99久久99久久久不卡 | 美女大奶头黄色视频| 亚洲成人一二三区av| 最近最新中文字幕免费大全7| 国产av一区二区精品久久| 黄色毛片三级朝国网站| 最近最新中文字幕大全免费视频 | 精品福利永久在线观看| 中文字幕精品免费在线观看视频 | 日韩免费高清中文字幕av| 性高湖久久久久久久久免费观看| 在线观看美女被高潮喷水网站| 午夜免费观看性视频| 国产高清国产精品国产三级| videosex国产| 黄色视频在线播放观看不卡| 又大又黄又爽视频免费| 成人综合一区亚洲| 亚洲欧洲日产国产| 日韩电影二区| 中文精品一卡2卡3卡4更新| 久久久a久久爽久久v久久| 五月伊人婷婷丁香| 国产精品久久久av美女十八| 校园人妻丝袜中文字幕| 日本-黄色视频高清免费观看| 日韩一区二区三区影片| 少妇人妻久久综合中文| 熟女av电影| 最近手机中文字幕大全| 日本av免费视频播放| 国产亚洲av片在线观看秒播厂| 91久久精品国产一区二区三区| 五月开心婷婷网| av不卡在线播放| 国产免费一级a男人的天堂| 欧美精品亚洲一区二区| 亚洲熟女精品中文字幕| 成人综合一区亚洲| 老熟女久久久| 国产激情久久老熟女| 在线天堂中文资源库| 色婷婷久久久亚洲欧美| 波多野结衣一区麻豆| 亚洲三级黄色毛片| 卡戴珊不雅视频在线播放| 亚洲av欧美aⅴ国产| 亚洲精品国产av蜜桃| 亚洲内射少妇av| 黄色一级大片看看| a 毛片基地| 亚洲国产av新网站| 热99久久久久精品小说推荐| 王馨瑶露胸无遮挡在线观看| 这个男人来自地球电影免费观看 | 久久久久国产精品人妻一区二区| 精品国产一区二区久久| 热99久久久久精品小说推荐| 国产在视频线精品| 亚洲天堂av无毛| 日韩视频在线欧美| 国产av国产精品国产| 天天操日日干夜夜撸| 中文字幕免费在线视频6| 色哟哟·www| 久久久国产精品麻豆| 人妻系列 视频| 成人国产av品久久久| 只有这里有精品99| 狠狠精品人妻久久久久久综合| 一边亲一边摸免费视频| 伦理电影大哥的女人| 亚洲人与动物交配视频| 亚洲国产av影院在线观看| 中国美白少妇内射xxxbb| 亚洲精品国产av成人精品| 欧美精品一区二区免费开放| 天天影视国产精品| 国产成人a∨麻豆精品| 蜜臀久久99精品久久宅男| 精品亚洲成国产av| 色视频在线一区二区三区| 寂寞人妻少妇视频99o| 中文字幕制服av| 欧美少妇被猛烈插入视频| 99国产精品免费福利视频| 91精品国产国语对白视频| 国产熟女欧美一区二区| 久久久久网色| 男女啪啪激烈高潮av片| 热re99久久精品国产66热6| 欧美日韩精品成人综合77777| 日本91视频免费播放| 日韩成人av中文字幕在线观看| 男女下面插进去视频免费观看 | 少妇人妻久久综合中文| 免费人成在线观看视频色| 久久久久人妻精品一区果冻| 久久精品国产亚洲av涩爱| 精品卡一卡二卡四卡免费| 亚洲欧美精品自产自拍| 飞空精品影院首页| 多毛熟女@视频| 美女脱内裤让男人舔精品视频| 另类精品久久| 97在线视频观看| 亚洲第一区二区三区不卡| 内地一区二区视频在线| 国产男人的电影天堂91|