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

    玄武巖分相Sm-Nd內(nèi)部等時線定年方法流程

    2014-07-19 12:13:35張利國段桂玲楊紅梅譚娟娟段瑞春邱嘯飛劉重芃
    巖礦測試 2014年5期
    關(guān)鍵詞:王水玄武巖同位素

    張利國, 段桂玲, 楊紅梅, 楊 梅, 譚娟娟, 段瑞春, 邱嘯飛, 劉重芃

    (中國地質(zhì)調(diào)查局武漢地質(zhì)調(diào)查中心同位素地球化學研究室, 湖北 武漢 430205)

    玄武巖分相Sm-Nd內(nèi)部等時線定年方法流程

    張利國, 段桂玲, 楊紅梅, 楊 梅, 譚娟娟, 段瑞春, 邱嘯飛, 劉重芃

    (中國地質(zhì)調(diào)查局武漢地質(zhì)調(diào)查中心同位素地球化學研究室, 湖北 武漢 430205)

    長期以來,對玄武巖精確測年一直是困擾地質(zhì)學家的重大科學問題。玄武巖結(jié)構(gòu)和組成特殊,巖石中礦物組成單一、鋯石十分稀少,顆粒很細,采用物理方法挑選單礦物和鋯石十分困難,很難應用內(nèi)部等時線法和鋯石U-Pb法研究其成巖時代。而全巖樣品間因巖漿分異產(chǎn)生的147Sm/144Nd比值差別很小,等時線年齡相對誤差較大;Rb含量很低,Rb/Sr比值很小,全巖Sm-Nd法、Rb-Sr法常常不能給出正確可信的年齡。根據(jù)內(nèi)部等時線法原理,本文通過化學方法,采用王水和氫氟酸-硝酸對玄武巖樣品進行分步溶解,分別對同一件樣品的王水溶解相、王水不溶相和全巖開展Sm-Nd同位素組成分析。結(jié)果表明:通過不同酸介質(zhì)分步溶解,可提取玄武巖中石英、透輝石、長石等礦物組合,該組合與其全巖具有相同的εNd(t)和一致的Nd同位素模式年齡;礦物與全巖構(gòu)筑的內(nèi)部等時線中,147Sm/144Nd比值的變化由全巖之間的0.005擴大到0.11,143Nd/144Nd值的變化由全巖的0.512500~0.512547擴大到0.512500~0.513145。通過該方法獲得了與已有鋯石U-Pb年齡在誤差范圍內(nèi)一致的Sm-Nd等時線年齡:t=(991±21) Ma(MSWD=2.1)。通過對比研究,本文認為:玄武巖分相Sm-Nd內(nèi)部等時線定年方法,適用于前寒武紀及更古老的玄武巖樣品的年齡測定。該方法的建立不僅有效提高玄武巖Sm-Nd等時線定年成功率,也為其他隱晶質(zhì)且不易挑出單礦物樣品的年齡測定提供了新的思路。

    玄武巖; 分步溶樣; 分相; Sm-Nd同位素; 年代學; 礦物內(nèi)部等時線

    玄武巖是地質(zhì)學研究的重要對象,其巖石地球化學及同位素特征與巖漿源區(qū)性質(zhì)、古大地構(gòu)造環(huán)境和深部地球動力學過程有著密切的聯(lián)系[1-4]。根據(jù)玄武巖的地球化學特征可以判別其形成的構(gòu)造環(huán)境[5-7],利用微量元素和稀土元素比值以及Sr-Nd-Pb同位素特征可成功識別出各種地幔端員,可用于反演巖漿作用的動力學過程[3,8-10]。然而,玄武巖一般結(jié)構(gòu)致密、組成均一,通常沒有巖漿結(jié)晶鋯石或含量甚微,不適于鋯石U-Pb定年;而其40Ar/39Ar法年齡測定[11],由于可能存在K和Ar的自然擴散作用或后期變質(zhì)、變形等多因素的擾動作用,很難測出早期地質(zhì)事件的年齡[12],因此,對玄武巖樣品年齡的精確測定一直是困擾地質(zhì)學家的科學難題[13]。

    相對于U-Pb、Rb-Sr及K-Ar體系,Sm、Nd同屬輕稀土元素,其化學性質(zhì)很相近,母體(147Sm)衰變形成的子體(143Nd)易在礦物晶格中保存下來,礦物或巖石中的Sm-Nd同位素體系易保持封閉,具有較強的抗蝕變和變質(zhì)作用的能力。因此,相對于其他同位素體系,Sm-Nd同位素等時線法測定基性巖、超基性巖的形成時代具有其顯著優(yōu)勢[12,14]。而由于Sm-Nd等時線法不僅可以獲得巖石形成的時間,而且可以獲得源區(qū)Nd同位素的初始值[12],長期以來,Sm-Nd同位素體系成功應用于隕石、月球巖石、古老的基性、超基性巖和前寒武紀老地層等領(lǐng)域方面的研究[12,14],是地球科學和宇宙科學研究中重要的工具之一。

    對于Sm-Nd同位素體系,用于制約等時線的147Sm/144Nd范圍越寬,該同位素體系能辨別的年齡越年輕[15]。如果巖石中的最大△147Sm/144Nd為0.05,在此情況下Sm-Nd同位素體系確定年齡絕對誤差約為40 Ma;如果巖石中最大△147Sm/144Nd達到0.1,則其確定的年齡絕對誤差可降至20 Ma左右[16]。但由于Sm和Nd均為輕稀土元素,離子半徑相似,化學性質(zhì)很相近,故自然界中Sm/Nd分餾通常不明顯[16-17];而且147Sm半衰期很長(1.06×1011a),衰變產(chǎn)生的143Nd積累有限,從而使不同來源樣品中Nd同位素組成的變異甚微(<2‰)[18]。這使得Sm-Nd同位素定年的應用在很大程度上受到限制。

    而Sm-Nd內(nèi)部等時線法彌補了全巖等時線的上述缺陷,如Lugmair等[19]從Angra Dos Reis(ADOR)輝石無球粒隕石里分選了3個磷酸鹽礦物和4個輝石類礦物,構(gòu)成了一條十分理想和精確的Sm-Nd內(nèi)部等時線。由于不同單礦物在結(jié)晶時其稀土元素存在分餾差異,因此同一樣品中單礦物的147Sm/144Nd變化范圍較全巖樣品大很多[12],從而獲得高成功率及準確性的Sm-Nd同位素等時線年齡。

    Sm-Nd內(nèi)部等時線法,大多是對采用物理方法分選出的巖石樣品中的單礦物進行分析。但是,對于玄武巖等隱晶結(jié)構(gòu)、組成均一的樣品,物理方法很難分選出其中的單礦物,使Sm-Nd法對玄武巖的測年受到限制。因此,本文嘗試通過化學方法,從玄武巖中分離出部分礦物與全巖構(gòu)筑內(nèi)部等時線。該方法先用王水溶解樣品,不溶部分再經(jīng)氫氟酸-硝酸溶解;然后,對王水溶解相、王水不溶相及全巖,分別進行Sm、Nd同位素分析,組成Sm-Nd礦物內(nèi)部等時線,以準確測定玄武巖年齡。

    1 實驗過程

    1.1 儀器及主要試劑

    固體熱表面電離質(zhì)譜儀(TRITON TI,美國Thermo公司)。

    實驗器皿:聚四氟乙烯坩堝、密封溶樣罐及高壓溶樣彈,石英玻璃交換柱。器皿在使用前,均先用去污粉或洗滌劑將其洗至內(nèi)、外表面不掛水珠或僅有少量細小水珠,再用一次去離子水(經(jīng)Milli-Q水純化系統(tǒng)純化)沖洗,裝入干凈的大燒杯中,依次用王水、硝酸、二次去離子水(一次去離子水再次經(jīng)過水純化系統(tǒng)純化)在電爐上煮至沸騰后保持其微沸4~6 h,于超凈化實驗室風干待用。

    試劑:鹽酸、硝酸、高氯酸、氫氟酸,均為市售優(yōu)級純(經(jīng)亞沸蒸餾)。

    1.2 樣品溶解

    稱取約300 mg粉末樣品(200目),加入適量6 mol/L王水,用低壓密封溶樣罐于110℃加熱48 h,待溶液冷卻至室溫,離心分離。取出上清液蒸干,用2.5 mol/L鹽酸提取,稱為王水溶解相①。殘余粉末置于高壓溶樣彈,加適量氫氟酸和硝酸(體積比10∶1原酸),密封于200℃溶樣48h,蒸干,再加入少許高氯酸加熱直至白煙冒盡,用2.5mol/L鹽酸提取,稱為王水不溶相②。

    1.3 稀土及Nd的分離純化

    稀土分離:采用AG50-8陽離子交換樹脂,交換柱下部樹脂床內(nèi)徑6mm、高80mm,上部盛液管內(nèi)徑為20mm,高100mm,總?cè)莘e40mL。樹脂用6mol/L鹽酸和亞沸蒸餾水交替再生3次。將王水溶解相①、王水不溶相②兩相溶液各分取約1/4(體積比),分別加入一定量145Nd-149Sm混合稀釋劑,混勻上柱,待溶液流干,加入6mL2.5mol/L鹽酸沖洗交換柱盛液管內(nèi)壁,待鹽酸滴干,再加入40mL2.5mol/L鹽酸淋洗Ca、Mg、Fe等基體元素,用15mL6mol/L鹽酸解吸稀土,蒸干解吸液,待質(zhì)譜測定,計算147Sm/144Nd值。

    Nd的純化:王水溶解相①、王水不溶相②兩相未加稀釋劑的溶液,經(jīng)稀土分離后,用P507樹脂(2-乙基己基膦酸單-2-乙基己基酯)分離純化Nd。交換柱規(guī)格與一次分離所用相同,樹脂床高60mm,用6mol/L鹽酸再生,亞沸蒸餾水淋洗至中性。待一次分離解吸液蒸干后,用0.1mol/L鹽酸提取、上柱。加入0.1mol/L鹽酸15mL淋洗La和Ce,用0.2mol/L鹽酸6mL解吸Nd。蒸干,待質(zhì)譜測定143Nd/144Nd值。

    1.4 質(zhì)譜測定

    質(zhì)譜測定采用雙燈絲熱表面電離。將錸燈絲用無水酒精清洗干凈,用點焊機將其點焊在燈絲支架上,將已點焊好燈絲的支架插裝在離子源轉(zhuǎn)盤上,并裝入燒帶裝置中,待真空抽至n×10-5Pa后按預定程序給燈絲供電,在3~5 A電流下燒15 min,以除去燈絲表面雜質(zhì)[20]。將樣品用10 μL去離子水轉(zhuǎn)移到燈絲中央,在1 A電流下蒸干,再緩慢增加電流至燈絲微紅,裝入離子源啟動真空系統(tǒng)。待質(zhì)譜儀的離子源真空達到要求(n×10-7Pa)后,打開通往分析管道的隔離閥,給電離帶燈絲加上電流,緩慢升溫,當燈絲溫度達到1800℃后,再緩慢加熱蒸發(fā)帶燈絲,注意尋找145Nd的離子流,并小心調(diào)節(jié)加到蒸發(fā)帶上的電流,使離子流達到足夠的強度(0.3~1 V)并保持穩(wěn)定,啟動自動測量程序,分別采集Sm、Nd同位素數(shù)據(jù)。本實驗采用靜態(tài)多接收采集數(shù)據(jù),其杯結(jié)構(gòu)如表1所示。

    表 1 靜態(tài)多接收質(zhì)譜儀(TRITON)杯結(jié)構(gòu)

    Table 1 Cup structure of the static multi-collectors mass spectrometer (TRITON)

    杯位置質(zhì)量數(shù)同位素杯位置質(zhì)量數(shù)同位素L3142142Ce+142NdH1146146NdL2143143NdH2147147SmL1144144Nd+144SmH3149149SmC145145Nd

    2 結(jié)果與討論

    2.1 Sm-Nd同位素分析結(jié)果

    2.1.1 標準及樣品

    BCR-2:玄武巖標準樣品(美國地質(zhì)調(diào)查局USGS),由美國地質(zhì)調(diào)查局采自美國俄勒岡州波特蘭以東地區(qū)的Bridal Veil Flow采石場;為現(xiàn)代噴發(fā)火山巖,無斑隱晶質(zhì)結(jié)構(gòu),形成時代為15 Ma[11]。

    EQB:峨眉山玄武巖樣品(La-Ce標準樣品[21]),采自中國四川峨眉山清音電站剖面,是大陸裂谷環(huán)境下的噴溢產(chǎn)物,其結(jié)構(gòu)為微晶玄武巖,形成于晚二疊世早期(259~262 Ma[22-23])。

    SNJ:玄武巖樣品,采自神農(nóng)架地區(qū)東北部姜家橋橋頭附近鄭家埡組地層中,具有間粒間隱結(jié)構(gòu),在該實驗開展之前,沒有年代學數(shù)據(jù)報道;而互層接觸的安山-英安質(zhì)熔結(jié)凝灰?guī)r年齡為1100~900 Ma[24]。

    2.1.2 分相礦物組成

    對玄武巖樣品進行分步溶解時,王水不溶相與全巖的礦物成分在中國地質(zhì)大學(武漢)地質(zhì)過程與礦產(chǎn)資源國家重點實驗室采用X射線粉晶衍射方法完成。分析結(jié)果見表2和圖1,結(jié)果表明:王水不溶相主要為石英、透輝石、長石的礦物組合;由全巖和王水不溶相估算,王水溶解相主要為綠泥石和少量高嶺石(表2)。

    2.1.3 分相Nd同位素組成

    本實驗在中國地質(zhì)調(diào)查局武漢地質(zhì)調(diào)查中心同位素地球化學研究室完成。所選樣品分別為:BCR-2、EQB和SNJ25、SNJ27、SNJ29、SNJ35,均具隱晶結(jié)構(gòu),礦物組成均一。

    質(zhì)譜測定過程中,每個試樣采集14組數(shù)據(jù),每組數(shù)據(jù)取10次掃描,單次掃描積分時間4.194 s,由聯(lián)機計算機取14組數(shù)據(jù)平均值。樣品147Sm/144Nd值的計算:分別采集149Sm/147Sm和145Nd/146Nd由稀釋法原理離線計算得出。143Nd/144Nd值測定:采用146Nd/144Nd=0.7219歸一化。所有樣品Nd同位素組成分析結(jié)果見表3。

    表 2 玄武巖王水不溶相及全巖X射線粉晶衍射半定量分析結(jié)果

    Table 2 Semi-quantitive results of X-ray powder diffraction analysis of the insoluble phase in aqua regia and whole-rock basalt

    樣品對應圖1的序號礦物含量(%)石英透輝石長石綠泥石高嶺石EQB(王水不溶相)A132491900SNJ27(王水不溶相)B116731100SNJ29(王水不溶相)C144371900SNJ35(王水不溶相)D134442200EQB(全巖)A220561905SNJ27(全巖)B212458350SNJ29(全巖)C2273315250SNJ35(全巖)D2223221205

    圖 1 玄武巖王水不溶相和全巖粉晶衍射分析圖譜Fig.1 The X-ray powder diffraction diagrams of undissolved basaltic phase in aqua regia and whole-rock basaltA1、B1、C1、D1分別為EQB,SNJ27、SNJ29、SNJ35王水不溶相圖譜;A2、B2、C2、D2為其對應全巖圖譜。

    表 3 玄武巖樣品及BCR-2標樣的Sm-Nd同位素分析結(jié)果
    Table 3 The Sm-Nd isotopic composition of basalt powder and BCR-2 standard samples dissolved with different acid media

    樣號類型147Sm/144Nd143Nd/144Nd(±2σ)tDM(Ga)t2DM(Ga)εNd(t)SNJ25-1f王水不溶相0.24520.513132±13-0.091.432.81SNJ25-2f王水不溶相0.24720.513145±20-0.021.432.79SNJ35-1f王水不溶相0.26660.513254±150.301.472.18SNJ35-2f王水不溶相0.26920.513264±210.311.492.01SNJ25全巖0.15160.512520±91.541.344.08SNJ25s全巖0.15210.512522±91.551.353.99SNJ35全巖0.15450.512547±71.551.344.18SNJ27全巖0.15680.512534±81.651.383.60SNJ29全巖0.15160.512504±121.581.373.75SNJ07-26-1全巖?0.15480.512525±21.611.373.7SNJ07-26-2全巖?0.15380.512500±11.651.393.4SNJ25-1c王水可溶相0.14660.512486±131.511.293.33SNJ25-2c王水可溶相0.14710.512483±131.521.303.21SNJ35-1c王水可溶相0.15090.512511±91.551.303.26SNJ35-2c王水可溶相0.15000.512519±121.511.313.13EQB全巖0.12100.512509±61.050.910EQB-c王水可溶相0.11610.512501±101.010.910.003EQB-c王水可溶相0.11600.512501±91.010.920.007EQB-f王水不溶相0.14400.512548±81.310.750.013EQB-f王水不溶相0.14220.512544±101.290.72-0.04BCR-2標準推薦值0.13820.512636±2??1.040.761.08BCR-2本文測定值0.13760.512633±41.040.731.15BCR-2-1c王水可溶相0.13060.512620±70.970.760.91BCR-2-2c王水可溶相0.13110.512641±20.940.761.31BCR-2-1f王水不溶相0.15340.512637±121.301.110.81BCR-2-2f王水不溶相0.15500.512644±141.311.120.92

    注:*snj07-26-1值引自文獻[24];**BCR-2值引自GeoRem (http://georem.mpch-mainz.gwdg.de),其147Sm/144Nd值根據(jù)Sm、Nd濃度為計算得出。tDM年齡采用文獻[25]線性模式計算,t2DM年齡采用李獻華等[26]兩階段模式計算。

    2.1.4 全巖及分相Sm-Nd內(nèi)部等時線年齡及誤差分析

    (1) 全巖及分相Sm-Nd內(nèi)部等時線

    根據(jù)“真假Sm-Nd等時線判別經(jīng)驗”:參與等時線年齡計算的各樣品間,有相同εNd(t)值或Nd同位素模式年齡,即可視為滿足Sm-Nd等時線同源、同時和封閉環(huán)境的條件,能參與等時線年齡計算[27]。由表3可以看出:分步溶解后,王水可溶相、不溶相及全巖具有一致的Nd同位素特征εNd(t),且王水可溶相、全巖和礦物組合的等時線(圖2c)與礦物內(nèi)部等時線(圖2b)結(jié)果完全一致。且前人大量研究表明:在巖石中,稀土元素能保持相對穩(wěn)定[12,14]。由此我們認為,玄武巖樣品通過分步溶解,不同相(王水可溶相、不溶相)與全巖之間構(gòu)筑的Sm-Nd礦物內(nèi)部等時線,符合Sm-Nd法測年的理論要求,能真實反映巖石的形成時代。

    圖 2 玄武巖分相Sm-Nd等時線Fig.2 Sm-Nd isochrones of fractional dissolved samplesa—SNJ全巖等時線;b—全巖+王水不溶相;c—全巖+王水不溶相+王水溶解相。

    等時線結(jié)果采用Ludwig(2001)[28]的Isoplot程序處理。根據(jù)經(jīng)驗,實驗室147Sm/144Nd和143Nd/144Nd的系統(tǒng)誤差分別為0.5%和0.003%。

    (2)Sm-Nd等時線年齡誤差分析

    由于Sm和Nd同屬稀土元素,在地質(zhì)作用過程中,Sm和Nd具有相似的地球化學特征和行為,從而導致Sm/Nd比值在巖石中的變化范圍很小[12,16]。用△147Sm/144Nd表示等時線各樣品點147Sm/144Nd值間的最大差值。Zhou等(2011)[22]報道的SNJ07-26-1、SNJ07-26-2和本文實測的SNJ25、SNJ27、SNJ29、SNJ35共6個全巖樣品之間△147Sm/144Nd值僅0.005,143Nd/144Nd在誤差范圍內(nèi)基本一致。

    但不同礦物在其形成階段,稀土元素內(nèi)部可發(fā)生強烈的分餾作用[29-33]。本文通過化學方法,對玄武巖樣品進行分步溶解后,得到石英、透輝石、長石等礦物組合的147Sm/144Nd值變化范圍較全巖樣品大,其△147Sm/144Nd值擴大到0.11。

    根據(jù)Sm-Nd年齡計算方程:

    143Nd/144Nd=143Nd/144Ndi+147Sm/144Nd(eλt- 1)

    (1)

    由最小二乘法原理可以得出如下等時線年齡誤差估算公式:

    △t=±σ143Nd/144Nd/λ△147Sm/144Nd

    (2)

    式中:λ為147Sm衰變常數(shù)=6.54×10-12a-1,σ為143Nd/144Nd比值測定精度(目前質(zhì)譜儀器對143Nd/144Nd比值測定外部精度可以達到±0.000015)[27]。由公式(2)可以估算:當全巖樣品間Δ147Sm/144Nd只有0.005時,Δt將達到±400 Ma,無法準確得到Sm-Nd等時線年齡(圖2a);而通過分步溶解,全巖和王水不溶相(石英、透輝石、長石的礦物組合)的147Sm/144Nd比值變化在0.1516~0.2692范圍內(nèi),Δ147Sm/144Nd絕對差值提高到0.11,由此估算的年齡誤差約為±20 Ma,實際得到的內(nèi)部等時線年齡誤差分別為±21 Ma和±19 Ma(圖2b,c)。

    由此可以看出:玄武巖樣品采用分步溶解,通過化學方法提取物理手段無法挑選的石英、透輝石、長石的礦物組合,與全巖之間構(gòu)成的Sm-Nd礦物內(nèi)部等時線,能有效地擴大Δ147Sm/144Nd值,減小年齡誤差,從而提高玄武巖Sm-Nd同位素年齡測定成功率。

    2.2 實驗標準與流程本底

    本文選取USGS玄武巖標樣BCR-2全巖,作為全流程準確性檢驗。其143Nd/144Nd=0.512633±4,與推薦值(見表3)完全一致;147Sm/144Nd =0.1376與推薦值0.1382誤差為4.3‰,小于1%,符合系統(tǒng)誤差限制。BCR-2全巖測定結(jié)果與推薦值一致,表明本實驗流程準確可靠。SNJ25、SNJ 35及BCR-2等平行雙份分析結(jié)果在誤差范圍內(nèi)基本一致(見表3),表明分步溶解實驗過程具有可重現(xiàn)性。

    質(zhì)譜儀器采用實驗室標準溶液JMC(溶液)監(jiān)控。經(jīng)過2年109次標樣監(jiān)測,其平均值143Nd/144Nd =0.511559±10(2σ,見圖3),標準偏差1×10-5;對標準物質(zhì)JNdi-1(金屬氧化物溶液)分析結(jié)果為143Nd/144Nd =0.512117±9,與Tanaka等(2000)[34]報道值0.512115±7完全一致。

    本實驗全流程空白:CNd≤1×10-10g/g,CSm≤0.7×10-10g/g。

    圖 3 JMC標準溶液測定結(jié)果Fig.3 The determination results of JMC standard solution

    2.3 Nd同位素測定的同質(zhì)異位素干擾

    在進行同位素測定時,通常會有同質(zhì)異位素的干擾。在Nd同位素分析過程中存在的同質(zhì)異位素有:142Ce與142Nd,144Sm與144Nd。

    本文分別采集149Sm/147Sm和145Nd/146Nd兩對同位素比值,通過稀釋法,離線計算樣品中147Sm、144Nd濃度,進而得出147Sm/144Nd值。在此過程中,質(zhì)量數(shù)142不參與計算,避免了142Ce與142Nd疊加的干擾。而143Nd/144Nd值,通過儀器直接測定143Nd、144Nd信號,采用146Nd/144Nd≡0.7219歸一化分餾校正,直接得到結(jié)果,質(zhì)量數(shù)為142的離子信號對143Nd/144Nd值的測定結(jié)果不產(chǎn)生任何影響。同時,楊岳衡等[35]研究工作表明大量Ce的存在并不影響Nd同位素測定。因此,142Ce對142Nd的干擾,在本文所用的147Sm/144Nd與143Nd/144Nd同位素體系中,不影響實驗結(jié)果。

    而144Sm會對143Nd/144Nd值測定產(chǎn)生直接的干擾。在質(zhì)譜測定過程中,通過監(jiān)測147Sm特征峰的信號,觀察144Sm對144Nd的影響。本實驗中,采用二次化學處理,分離Nd與Sm,以消除144Sm對144Nd的干擾。實驗結(jié)果表明:在進行143Nd/144Nd值測定過程中,147Sm/144Nd≤2×10-5(如表4所示),甚至為負值,而自然界中147Sm/144Sm=4.8827[36],可得144Sm/144Nd≤0.4×10-5,基本消除了144Sm對144Nd的干擾。

    表 4143Nd/144Nd值測定時144Sm干擾監(jiān)控值

    Table 4143Nd/144Nd values determined with144Sm disturbance monitoring

    掃描次數(shù)145Nd/144Nd143Nd/144Nd142Nd/144Nd147Sm/144Nd10.3484940.5124941.186942.06673×10-520.3484480.524641.17460-2.38687×10-630.3484570.5125231.170361.22418×10-540.3484410.5125511.166716.49628×10-550.3484570.5125361.163781.23097×10-560.3484460.5124681.161271.45124×10-570.3484400.5125641.159368.56490×10-680.3484520.5125191.157711.43829×10-590.3484240.5124961.156291.44144×10-6100.3484510.5125201.154912.77971×10-6110.3484320.5125671.154011.37760×10-6120.3484360.5125641.153049.42705×10-6130.3484300.5125091.152221.96576×10-5140.3484120.5125101.151372.12260×10-5平均值0.3484400.5125201.159661.01927×10-5

    3 實驗方法的適用性

    3.1 Sm-Nd同位素體系特點

    根據(jù)Sm-Nd等時線原理,酸性、基性到超基性的火成巖、變質(zhì)巖和沉積巖類均可用于Sm-Nd法年齡的測定和研究。由于全巖的Sm-Nd比值通常變化窄小[12,16-17],不易于獲得理想的Sm-Nd全巖等時線年齡結(jié)果,僅可以用于模式年齡測定。因此,斜方輝石、單斜輝石、斜長石和磷酸鹽礦物等是Sm-Nd內(nèi)部等時線年齡測定的常用對象[12]。

    根據(jù)Sm-Nd同位素放射性原理,由公式(2)可以估算:一組具有相同初始143Nd/144Nd比值的樣品,當△147Sm/144Nd值為0.1與0.01時,分別需要母體同位素147Sm經(jīng)過20 Ma和200 Ma的放射積累,△143Nd/144Nd值才能勉強超出儀器測量誤差(±0.000015)。因此,Sm-Nd法更適用于古老巖石的測年與成因研究[12]。

    3.2 分相Sm-Nd同位素方法適用范圍

    為比較不同形成時代樣品對分相實驗的影響,在本文分別選取USGS玄武巖標樣BCR-2(15 Ma)[11]、EQB(峨眉山玄武巖262~259 Ma)[22-23]和SNJ25、SNJ27、SNJ29、SNJ35(1100~900 Ma)[24]三組玄武巖。

    由實驗結(jié)果可以看出:BCR-2與EQB分相后△147Sm/144Nd值僅為0.024、0.026,其估算年齡誤差達到100 Ma。而BCR-2形成時代只有15 Ma,其測定年齡誤差遠大于其自身形成時代,因此,不適合Sm-Nd同位素年齡測定。EQB樣品(260 Ma)分步溶解后,根據(jù)等時線計算原理,其年齡誤差也達到82 Ma(圖4),相對誤差達到32%;所以,形成時代為260 Ma的EQB樣品,同樣不適用于該方法的測定。而SNJ25、SNJ27、SNJ29、SNJ35樣品,其形成時代較老,且分相后△147Sm/144Nd值可達到0.11,通過誤差理論估算,其等時線年齡誤差可降低至±20 Ma,可以得到比較準確的等時線年齡(圖2b)。

    圖 4 EQB分相等時線圖Fig.4 Sm-Nd isochrons of fractional dissolved samples

    通過分析和對比,本文分步溶解實驗方法對260 Ma及更年輕玄武巖樣品,無法得到準確可靠的等時線年齡;而對980 Ma的SNJ玄武巖,得到了良好的Sm-Nd內(nèi)部等時線年齡。其原因大致有:①樣品自身形成年齡小,以致放射成因143Nd累積有限,分相前后樣品的143Nd/144Nd值變化不明顯;②實驗方法,通過本實驗化學流程,無法實現(xiàn)單一礦物的分離,無法完整體現(xiàn)單礦物之間Sm/Nd分餾差別,使樣品分相前后其147Sm/144Nd及143Nd/144Nd變化較小。

    綜上所述,我們認為:分步溶解分相Sm-Nd內(nèi)部等時線測定方法,適用于前寒武紀和更古老的玄武巖樣品的年齡測定。

    4 結(jié)語

    傳統(tǒng)Sm-Nd內(nèi)部等時線法,采用物理方法分選出巖石樣品中的單礦物,分別對礦物和全巖開展Sm-Nd同位素分析。然而,對于玄武巖等隱晶結(jié)構(gòu)、組成均一的樣品,物理方法很難分選出其中的單礦物,使Sm-Nd內(nèi)部等時線法對玄武巖的測年受到限制。本文通過化學方法,對玄武巖樣品采用王水和氫氟酸-硝酸分步溶解,分離出玄武巖樣品中物理方法難以挑選的礦物,實驗結(jié)果令人滿意。

    (1)與全巖Sm-Nd同位素分析結(jié)果相比,該方法有效擴大了等時線中147Sm/144Nd比值差別,分步溶解得到的礦物組合與全巖組成的內(nèi)部等時線中,147Sm/144Nd比值差別由全巖的0.005擴大到0.11,143Nd/144Nd值由全巖的0.512500~0.512547擴大到0.512500~0.513145,成功測定了SNJ玄武巖樣品的同位素年齡:t=(991±21) Ma(MSWD=2.1)。

    (2)通過對形成時代不同的三組樣品的對比研究,我們認為:本實驗為前寒武紀及更古老的玄武巖樣品同位素年齡測定提供了新的方法,也為其他隱晶質(zhì)且不易挑出單礦物的樣品年齡測定提供了新的思路。

    (3)本實驗在方法上,通過本實驗化學流程,無法實現(xiàn)單一礦物的分離,無法完整體現(xiàn)單礦物之間Sm/Nd分餾差別。因此,如何通過化學方法分選出巖石中更多的單一礦物,有待進一步的深入研究。

    [1] Momme P, Tegner C C, Broks C K, Keays R R.Two melting regimes during Paleocene flood basalt generation in East Greenland: Combined REE and PGE modeling[J].ContributionstoMineralogyandPetrology,2006,15(1):88-100.

    [2] Sigrnarssona O, Steinthórsson S.Origin of Icelandic basalts:A review of their petrology and geochemistry [J].JournalofGeodynarnics, 2007,43(1):87-100.

    [3] Kinmana W S, Neal C R, Davidson J P, Font L.The dynamics of Kerguelen Plateau magma evolution:New insights from major element, trace element and Sr isotope microanalysis of plagioclase hosted in Elan Bank bssalts [J].ChemicalGeology,2009,264:247-265.

    [4] Fan W M, Zhang C H, Wang Y J, Guo F, Peng T P.Geochronology and geochemistry of Permian basalts in western Guangxi Province Southwest China:Evidence for plume-lithosphere interaction[J].Lithos,2009,102:218-236.

    [5] Pearce J A.Role of the sub-continental lithosphere in magma genesis at active continental margins[C]// Hawkesworth, Norry.Continental Basalts and Mantle Xenoliths.Nantwich: Shiva,1983:230-249.

    [6] 張樹明,王方正.玄武巖在研究巖石圈深部過程及構(gòu)造背景中的應用[J].地球科學進展,2002,17(5):685-692.

    [7] 厲子龍,楊樹鋒,陳漢林, Langmiur C H,余星,林秀斌,勵音騏.塔西南玄武巖年代學和地球化學特征及其對二疊紀地幔柱巖漿演化的制約[J].巖石學報,2008,24(5):959-970.

    [8] 肖龍,徐義剛,梅厚鈞,沙紹禮.云南賓川地區(qū)峨眉山玄武巖的地球化學特征及時空演化[J].地質(zhì)科學,2003,38(4):478-494.

    [9] Xu Y G, He B, Hnang X L, Luo Z Y, Cheng S L, Xiao L, Zhu D, Shao H, Fan W M, Xu J F, Wang Y J. Identification of mantle plumes in Emeishan Large Igneous Province[J].Episode,2007,30(1):32-42.

    [10] 勵音騏,厲子龍,毛建仁,駱學全,徐貽贛,汪惠惠.江西會昌崠背玄武巖地球化學和Sr-Nd同位素特征及其地質(zhì)意義[J].礦床地質(zhì),2009, 28(5):696-706.

    [11] Brueseke M E, Heizler M T, Hart W K, Mertzman S A.Distribution and geochronology of Oregon Plateau (U.S.A.) flood basalt volcanism: The Steens Basalt revisited[J].JournalofVolcanologyandGeothermalResearch,2007,161:187-214.

    [12] 陳文,萬渝生,李華芹,張宗清,戴橦謨,施澤恩,孫敬博.同位素地質(zhì)年齡測定技術(shù)及應用[J].地質(zhì)學報,2011,85(11):1918-1947.

    [13] 羅照華,莫宣學,萬渝生,李莉,魏陽.青藏高原最年輕堿性玄武巖SHRIMP年齡的地質(zhì)意義[J].巖石學報,2006,22(3):578-584.

    [14] 彭建堂,符亞洲,袁順達,沈能平,張東亮.熱液礦床中含鈣礦物的Sm-Nd同位素定年[J].地質(zhì)論評,2006, 52(5):662-667.

    [15] Farmer G L, DePaolo D J. Source of hydrothermal com-ponents: Heavy Isotopes[M]//Barnes H L. Geochemistry of Hydrothermalre Deposits (Third Edition).New Jersey:John Wiley & Sons, 1997:31-41.

    [16] DePaolo D J. Neodymium Isotope Geochemistry: An In-troduction [M].Berlin Heidelberg: Springer Verlag,1988.

    [17] Pimentel M M, Charnley N.Intracrustal REE fractiona-tion and implications for Sm-Nd model age calculations in late-stage granitic rocks: An example from central Brazil[J].ChemicalGeology,1991,86:123-138.

    [18] 毛存孝,朱乃娟,朱丙泉,施澤恩,黃榮生.Sm-Nd同位素體系測定技術(shù)與應用討論[J].地球化學,1989(1):36-42.

    [19] Lugmair G W, Marti K.Sm-Nd-Pu Timepieces in the Angra dos Reis meteorite[J].EarthandPlanetaryScienceLetters,1977,35:273-284.

    [20] 楊紅梅,路遠發(fā),段桂玲,呂紅,劉 焰,章麗娟,梅玉萍,馬麗艷.茶葉中鉛同位素比值的測定方法[J].地球化學,2005,34(4):373-378.

    [21] 楊紅梅,段桂玲,凌文黎.La-Ce法巖石標準物質(zhì)和Ce同位素標準溶液研制[J].地球化學,2009(2):179-186.

    [22] Zhou M F, Malpas J, Song X Y, Robinson P T, Sun M, Kennedy A K, Lesher C M, Keays R R.A temporal link between the Emeishan large igneous province (SW China) and the end-Guagalupian mass extinction[J].EarthandPlanetaryScienceLetters,2002,196:113-122.[23] Guo F, Fan W M, Wang Y J, Li C W.When did the Emeishan mantle plume activity start? Geochronological and geochemical evidence from ultramafic-mafic dikes in southwestern China[J].InternationalGeologyReview, 2004,46:226-234.

    [24] Qiu X F, Ling W L, Liu X M , Kusky T, Berkana W, Zhang Y H, Gao Y J, Lu S S, Kuang H, Liu C X. Recognition of Grenvillian volcanic suite in the Shennongjia region and its tectonic significance for the South China Craton[J].PrecambrianResearch,2011,191:101-119.

    [25] Goldstein D J, O′Nions R K, Hamillton P J.A Sm-Nd isotopic study of atmospheric dusts and particulates from major river system [J].EarthandPlanetaryScienceLetters,1984,70(1):221-236.

    [26] 李獻華,趙振華,桂訓唐.華南前寒武紀地殼形成時代的Sm-Nd和鋯石U-Pb同位素制約[J].地球化學,1991(3):253-263.

    [27] 李志昌,路遠發(fā),黃圭成.放射性同位素地質(zhì)學方法與進展[M].武漢:中國地質(zhì)大學出版社,2004.

    [28] Ludwig K R.Isoplot/Ex, rev.2.49.A Geochronological Toolkit for Microsoft Excel[M].BerkeleyGeochronologyCenter,2001:56.

    [29] 李華芹,劉家齊,杜國民,魏林.內(nèi)生金屬礦床成礦作用年代學研究——以西華山鎢礦為例[J].科學通報,1992, 37(12):109-1112.

    [30] Chesley J T, Halliday A N, Kyser T K, Spry P G.Direct dating of MVT mineralization: Use of Sm-Nd in fluorite [J].EconomicGeology,1994,89:1192-1199.

    [31] Brugger J, Mass R, Lahaye Y, McRae C, Ghaderi M, Costa S, Lam-bert D, Bateman R, Prince K.Origins of Nd-Sr-Pb isotopic variations in single scheelite grains from Archaean gold deposits, Western Australia[J].ChemicalGeology,2002,182:203-225.

    [32] 彭建堂,胡瑞忠,林源賢,趙軍紅.錫礦山銻礦床熱液方解石的Sm-Nd同位素定年[J].科學通報, 2002,47(10):789-792.

    [33] 彭建堂,胡瑞忠,蔣國豪.螢石Sm-Nd同位素體系對晴隆銻礦床成礦時代和物源的制約[J].巖石學報, 2003,19(4):785-791.

    [34] Tanaka T, Togashi S, Kamioka H, Amakawa H, Kagami H, Hamamoto T, Yuhara M, Orihashi Y J, Yoneda S, Shimizu H, Kunimaru T, Takahashi K, Yanagi T, Nakano T, Fujimaki H, Shinjo R, Asahara Y, Tanimizu M, Dragusanu C.JNdi-1: A neodymium isotopic reference in consistency with LaJolla neodymium [J].ChemicalGeology,2000,168:279-281.

    [35] 楊岳衡,張宏福,謝烈文,吳福元.多接收器電感耦合等離子體質(zhì)譜精確測定釹同位素組成[J].分析化學,2007,35(1):71-74.

    [36] O′Nions R K, Hamilton P J, Evensen N M.Variations in143Nd/144Nd and87Sr/86Sr ratios in ocealts[J].EarthandPlanetaryScienceLetters,1977,34:13-22.

    Procedure for Sm-Nd Internal Isochron Dating Method on Basalt by Phase Separation

    ZHANGLi-guo,DUANGui-ling,YANGHong-mei,YANGMei,TANJuan-juan,DUANRui-chun,

    QIUXiao-fei,LIUChong-peng

    ( Isotopic Geochemistry Laboratory, Wuhan Center of Geological Survey, China Geological Survey, Wuhan 430205, China)

    Accurately dating the age of basalt has been an important scientific issue for geologists. Due to its unique texture and composition, the basaltic minerals are limited to certain types, with scarce zircons in tiny grain sizes. It is difficult to pick pure minerals such as zircon from basalt in physical ways, therefore, studying the formation age of basalt by internal isochron or zircon U-Pb methods is challenging. Since whole rock basalt samples, on the other hand, have a limited range of147Sm/144Nd ratios, low Rb content and low Rb/Sr ratio, the whole rock Sm-Nd method or Rb-Sr method usually cannot yield an accurate age. In this study, the principle of internal isochron and chemical methods, using aqua regia and Hf-HNO3to dissolve the same basalt samples by steps is utilized, and performed by Sm-Nd isotopic analysis on the 3 phases from the same sample: 1) dissolved phase in aqua regia, 2) undissolved phase in aqua regia and 3) whole rock phase. The results shows that, through dissolving the samples in different acid media by steps, quartz-diopside-feldspar mineral assemblage can be extracted from basalt, and this assemblage has the sameεNd(t) and Nd modal age as the whole rock. In the internal isochron constructed by minerals and whole rock, the variation range of147Sm/144Nd ratio increases from 0.005 within the whole rock to 0.11, the143Nd/144Nd ratio range varies from 0.512500-0.512547 within the whole rock to 0.512500-0.513145. The Sm-Nd age yielded with this method is identical to previous zircon U-Pb age within uncertainty:t=(991±21) Ma, MSWD=2.1. By comparing the results from conditional experiments, we propose that Sm-Nd internal isochron dating by phase separation is more suitable for rock samples of Precambrian or older. The establishment of this method can effectively increase the success of Sm-Nd isochron dating for basalt, as well as offer a new estimate for other aphanitic rock samples that cannot provide pure minerals for easy dating.

    basalt; dissolve by steps; phase separation; Sm-Nd isotope; geochronology; internal isochron

    2013-12-04;

    2014-07-15; 接受日期: 2014-07-20

    中國地質(zhì)大調(diào)查項目(1212011120272,1212011121102)

    張利國,工程師,分析化學專業(yè),主要從事Sm-Nd同位素研究。E-mail: whzliguo@163.com。

    0254-5357(2014)05-0640-09

    P588.145; P597.3

    A

    猜你喜歡
    王水玄武巖同位素
    集合中的易錯問題剖析
    玄武巖纖維可用于海水淡化領(lǐng)域
    火山作用對板塊構(gòu)造環(huán)境的判別方法
    2019年全球連續(xù)玄武巖纖維市場產(chǎn)值將達1.047億美元
    捉“妖”記(二)——讓“數(shù)列妖”無處藏身
    冷暖
    冷暖
    短篇小說(2015年10期)2015-11-03 04:22:31
    深空探測用同位素電源的研究進展
    《同位素》(季刊)2015年征訂通知
    同位素(2014年3期)2014-06-13 08:22:28
    硼同位素分離工藝與生產(chǎn)技術(shù)
    同位素(2014年2期)2014-04-16 04:57:15
    欧美成人精品欧美一级黄| 亚洲欧美精品专区久久| 麻豆av噜噜一区二区三区| 最近2019中文字幕mv第一页| 亚洲人与动物交配视频| 亚洲国产精品成人综合色| av在线播放精品| 国产精品久久久久久av不卡| 狂野欧美白嫩少妇大欣赏| 啦啦啦啦在线视频资源| 日韩欧美国产在线观看| 亚洲欧美成人精品一区二区| 成人特级黄色片久久久久久久| 人人妻人人澡人人爽人人夜夜 | 日本一二三区视频观看| 欧美成人精品欧美一级黄| 又黄又爽又刺激的免费视频.| 国产在线男女| 三级经典国产精品| 亚洲精品自拍成人| 国产精品国产三级国产av玫瑰| 亚洲国产欧美人成| 成人午夜高清在线视频| 亚洲国产精品成人久久小说 | 久久久久久伊人网av| 婷婷六月久久综合丁香| 麻豆国产97在线/欧美| 亚洲aⅴ乱码一区二区在线播放| 国产片特级美女逼逼视频| 免费看av在线观看网站| 午夜视频国产福利| 美女脱内裤让男人舔精品视频 | 一本一本综合久久| 一卡2卡三卡四卡精品乱码亚洲| 啦啦啦观看免费观看视频高清| 国产精品一区二区三区四区免费观看| 你懂的网址亚洲精品在线观看 | 人人妻人人澡欧美一区二区| 亚洲18禁久久av| av在线播放精品| 少妇熟女欧美另类| 大又大粗又爽又黄少妇毛片口| 不卡视频在线观看欧美| 亚洲成人av在线免费| 免费观看的影片在线观看| 亚洲精品乱码久久久v下载方式| 男女做爰动态图高潮gif福利片| 国产色爽女视频免费观看| 欧美性感艳星| 国产探花极品一区二区| 国产精品精品国产色婷婷| 一区二区三区高清视频在线| 夫妻性生交免费视频一级片| 欧洲精品卡2卡3卡4卡5卡区| 久99久视频精品免费| 一级毛片久久久久久久久女| 日韩精品有码人妻一区| 亚洲中文字幕一区二区三区有码在线看| 亚洲精华国产精华液的使用体验 | 天堂中文最新版在线下载 | 一夜夜www| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产在线男女| 中文在线观看免费www的网站| 成人二区视频| 桃色一区二区三区在线观看| 国产色爽女视频免费观看| 婷婷精品国产亚洲av| 亚洲成人精品中文字幕电影| 91久久精品国产一区二区三区| 国产亚洲精品久久久com| 亚洲国产精品久久男人天堂| 欧美不卡视频在线免费观看| 18禁黄网站禁片免费观看直播| 国产一区二区三区在线臀色熟女| 午夜爱爱视频在线播放| 国产欧美日韩精品一区二区| 国产精品综合久久久久久久免费| 亚洲熟妇中文字幕五十中出| 午夜精品在线福利| 99热这里只有是精品50| 夜夜看夜夜爽夜夜摸| 国产高清激情床上av| 校园人妻丝袜中文字幕| 欧美色欧美亚洲另类二区| 亚洲三级黄色毛片| 色哟哟哟哟哟哟| 2021天堂中文幕一二区在线观| 精品日产1卡2卡| 精品久久久久久久久亚洲| 在线播放国产精品三级| 高清毛片免费观看视频网站| 91久久精品电影网| 国产午夜福利久久久久久| 亚洲欧美日韩东京热| 国产日韩欧美在线精品| 天堂影院成人在线观看| 亚洲在线观看片| 久久久久免费精品人妻一区二区| 国产精品一区二区在线观看99 | 国产男人的电影天堂91| 青青草视频在线视频观看| av在线天堂中文字幕| av在线亚洲专区| 欧美一级a爱片免费观看看| 欧美另类亚洲清纯唯美| 久久精品国产99精品国产亚洲性色| 久久精品久久久久久久性| 亚洲图色成人| 又黄又爽又刺激的免费视频.| 尾随美女入室| h日本视频在线播放| 天堂中文最新版在线下载 | 人人妻人人澡人人爽人人夜夜 | 久久亚洲精品不卡| 色视频www国产| 久久九九热精品免费| 性色avwww在线观看| 一个人观看的视频www高清免费观看| 色综合站精品国产| 国产一区二区激情短视频| 狂野欧美激情性xxxx在线观看| 两个人视频免费观看高清| 天美传媒精品一区二区| 国产黄片视频在线免费观看| av又黄又爽大尺度在线免费看 | 亚洲欧美中文字幕日韩二区| 亚洲欧美中文字幕日韩二区| 久久久国产成人精品二区| 蜜桃久久精品国产亚洲av| 国语自产精品视频在线第100页| 久久这里只有精品中国| 日韩欧美精品免费久久| 青春草视频在线免费观看| 久久久久久九九精品二区国产| 精品久久久久久久久av| 秋霞在线观看毛片| 久久鲁丝午夜福利片| 91aial.com中文字幕在线观看| 天美传媒精品一区二区| 直男gayav资源| 婷婷精品国产亚洲av| 久久这里只有精品中国| 99久国产av精品国产电影| 青春草视频在线免费观看| 天天一区二区日本电影三级| 白带黄色成豆腐渣| 夜夜看夜夜爽夜夜摸| 亚洲无线观看免费| 在线播放无遮挡| 精品不卡国产一区二区三区| 成人毛片60女人毛片免费| 日韩人妻高清精品专区| 日产精品乱码卡一卡2卡三| 97在线视频观看| 日韩一区二区三区影片| 两性午夜刺激爽爽歪歪视频在线观看| 人妻夜夜爽99麻豆av| 少妇高潮的动态图| 小说图片视频综合网站| 婷婷精品国产亚洲av| 欧美变态另类bdsm刘玥| 欧美另类亚洲清纯唯美| 国产av麻豆久久久久久久| 亚洲欧美精品自产自拍| 国产三级在线视频| 日日撸夜夜添| 亚洲aⅴ乱码一区二区在线播放| 日韩亚洲欧美综合| 99热这里只有是精品在线观看| 午夜久久久久精精品| 久久久久国产网址| 久久久午夜欧美精品| 黄片wwwwww| 婷婷色综合大香蕉| 波多野结衣高清无吗| 精品久久久久久久久亚洲| 91aial.com中文字幕在线观看| 3wmmmm亚洲av在线观看| 九九久久精品国产亚洲av麻豆| 一本久久中文字幕| 午夜视频国产福利| 成人三级黄色视频| 精品久久久久久久末码| 嘟嘟电影网在线观看| 久久久久久久久久久丰满| 青春草视频在线免费观看| 91麻豆精品激情在线观看国产| 日韩欧美一区二区三区在线观看| 中文字幕熟女人妻在线| 国产精品久久视频播放| 成人av在线播放网站| 国产单亲对白刺激| 亚洲欧美日韩高清在线视频| 久久久久久久久久久免费av| av天堂在线播放| 亚洲七黄色美女视频| 久久精品国产亚洲av天美| 亚洲电影在线观看av| 国产一区二区在线av高清观看| 麻豆精品久久久久久蜜桃| 免费看光身美女| 久久久成人免费电影| 国产亚洲av嫩草精品影院| 国产欧美日韩精品一区二区| 日本色播在线视频| 日本一二三区视频观看| 国产伦在线观看视频一区| 搡女人真爽免费视频火全软件| 国产亚洲精品久久久com| 欧美最黄视频在线播放免费| 男女啪啪激烈高潮av片| 波多野结衣高清作品| 中文字幕熟女人妻在线| 亚洲av熟女| 久久国内精品自在自线图片| 国产视频内射| 免费观看精品视频网站| 欧美人与善性xxx| 国产精品1区2区在线观看.| 日韩精品青青久久久久久| 丰满的人妻完整版| 色综合站精品国产| 我要看日韩黄色一级片| 一本久久精品| 中出人妻视频一区二区| 青春草视频在线免费观看| 国产精品.久久久| 日韩精品青青久久久久久| 91久久精品国产一区二区三区| 午夜激情福利司机影院| 亚洲经典国产精华液单| 美女xxoo啪啪120秒动态图| 国产在视频线在精品| 日产精品乱码卡一卡2卡三| 能在线免费看毛片的网站| 久久久欧美国产精品| 久久久成人免费电影| 国产视频内射| 人妻夜夜爽99麻豆av| 欧美激情在线99| 国产在线男女| 亚洲18禁久久av| 最近2019中文字幕mv第一页| 老女人水多毛片| 波多野结衣巨乳人妻| 日韩一本色道免费dvd| 国产精品国产三级国产av玫瑰| 日日干狠狠操夜夜爽| 色视频www国产| 国产av麻豆久久久久久久| 亚洲经典国产精华液单| 看非洲黑人一级黄片| 99国产精品一区二区蜜桃av| 中文字幕制服av| 国产免费男女视频| 欧美bdsm另类| 一级黄片播放器| 中文亚洲av片在线观看爽| 超碰av人人做人人爽久久| 日本一本二区三区精品| 26uuu在线亚洲综合色| 精品人妻熟女av久视频| 小说图片视频综合网站| 欧美激情在线99| 少妇的逼水好多| 久久99热这里只有精品18| 国产精品久久久久久精品电影小说 | 校园人妻丝袜中文字幕| 精品久久久久久成人av| 亚洲成人av在线免费| 欧美bdsm另类| 欧美激情在线99| 啦啦啦韩国在线观看视频| 精品国产三级普通话版| 搡老妇女老女人老熟妇| 成人鲁丝片一二三区免费| 亚洲在久久综合| 精品久久久久久久久亚洲| 69人妻影院| 免费黄网站久久成人精品| 国产三级中文精品| 在线观看午夜福利视频| 可以在线观看毛片的网站| 菩萨蛮人人尽说江南好唐韦庄 | 热99在线观看视频| av在线天堂中文字幕| 日本一二三区视频观看| 一级av片app| 黄色欧美视频在线观看| 国产日韩欧美在线精品| 日日摸夜夜添夜夜添av毛片| 欧美不卡视频在线免费观看| 免费观看在线日韩| 91午夜精品亚洲一区二区三区| 欧美日韩在线观看h| 美女xxoo啪啪120秒动态图| 男人狂女人下面高潮的视频| 国产高清不卡午夜福利| 综合色av麻豆| 一卡2卡三卡四卡精品乱码亚洲| 全区人妻精品视频| 亚洲熟妇中文字幕五十中出| 国产精品久久视频播放| 非洲黑人性xxxx精品又粗又长| 欧洲精品卡2卡3卡4卡5卡区| 毛片一级片免费看久久久久| 欧美xxxx黑人xx丫x性爽| 亚洲va在线va天堂va国产| 国产片特级美女逼逼视频| 欧美成人免费av一区二区三区| 老女人水多毛片| av专区在线播放| 国产精品野战在线观看| 欧美+日韩+精品| 大香蕉久久网| 色哟哟哟哟哟哟| 亚洲真实伦在线观看| 久久久久久伊人网av| 舔av片在线| 午夜福利在线观看吧| 国产片特级美女逼逼视频| 晚上一个人看的免费电影| 可以在线观看毛片的网站| 夫妻性生交免费视频一级片| 日日啪夜夜撸| 九色成人免费人妻av| 国产精品久久视频播放| 日韩 亚洲 欧美在线| 亚洲国产精品久久男人天堂| 国产亚洲av片在线观看秒播厂 | 亚洲av不卡在线观看| 内地一区二区视频在线| 成人欧美大片| 成人亚洲欧美一区二区av| 亚洲一区二区三区色噜噜| 免费av毛片视频| 欧美性感艳星| 能在线免费看毛片的网站| 久久精品国产99精品国产亚洲性色| 久久人妻av系列| av.在线天堂| 午夜福利在线在线| 哪里可以看免费的av片| 欧美日韩乱码在线| 婷婷色av中文字幕| 亚洲国产精品合色在线| 国产精品福利在线免费观看| 人人妻人人看人人澡| 国产伦精品一区二区三区四那| 国产午夜精品一二区理论片| 欧美不卡视频在线免费观看| 菩萨蛮人人尽说江南好唐韦庄 | 欧美又色又爽又黄视频| 精品一区二区免费观看| 亚洲欧美清纯卡通| 国产在视频线在精品| 在线观看66精品国产| 国产精品无大码| 超碰av人人做人人爽久久| 一进一出抽搐动态| 国产成人a区在线观看| 少妇高潮的动态图| 五月伊人婷婷丁香| avwww免费| 欧美一区二区精品小视频在线| 哪个播放器可以免费观看大片| 麻豆国产97在线/欧美| 日韩 亚洲 欧美在线| 久久久久久久久久成人| 亚洲av中文av极速乱| 久久亚洲精品不卡| 一区二区三区免费毛片| 亚洲国产精品久久男人天堂| 神马国产精品三级电影在线观看| 久99久视频精品免费| 网址你懂的国产日韩在线| 欧美在线一区亚洲| 国产午夜精品论理片| 国产一区二区亚洲精品在线观看| 精品久久久久久久久久免费视频| 久久精品国产清高在天天线| 内地一区二区视频在线| av免费观看日本| 欧美成人免费av一区二区三区| 永久网站在线| 国产老妇伦熟女老妇高清| 黄色日韩在线| 免费大片18禁| 久久草成人影院| 成人欧美大片| 黄色日韩在线| 精品久久久久久久久亚洲| 亚洲精品粉嫩美女一区| 一级毛片久久久久久久久女| 国产黄片视频在线免费观看| 国产精品久久久久久av不卡| 午夜激情欧美在线| 久久久欧美国产精品| 99久久久亚洲精品蜜臀av| 免费电影在线观看免费观看| 色播亚洲综合网| 国产精品久久久久久精品电影小说 | 深夜a级毛片| 国内精品宾馆在线| 成人毛片a级毛片在线播放| 一边摸一边抽搐一进一小说| 青春草视频在线免费观看| 日日啪夜夜撸| 99久久无色码亚洲精品果冻| 国产日韩欧美在线精品| 12—13女人毛片做爰片一| 日韩欧美三级三区| av卡一久久| 久久久a久久爽久久v久久| 国产精品一区二区性色av| 久久亚洲精品不卡| 国产视频首页在线观看| 嘟嘟电影网在线观看| 国产午夜精品论理片| 国产精品野战在线观看| 大型黄色视频在线免费观看| 菩萨蛮人人尽说江南好唐韦庄 | 久久久久免费精品人妻一区二区| 网址你懂的国产日韩在线| 日韩欧美三级三区| 久久精品综合一区二区三区| 久久久国产成人免费| 日本在线视频免费播放| 成人亚洲欧美一区二区av| 97超视频在线观看视频| 观看美女的网站| 黄色视频,在线免费观看| 99在线视频只有这里精品首页| 青春草视频在线免费观看| 少妇被粗大猛烈的视频| 久久久久九九精品影院| a级一级毛片免费在线观看| av在线蜜桃| 三级经典国产精品| 99热网站在线观看| 18禁在线无遮挡免费观看视频| 国产高清三级在线| 日本成人三级电影网站| 99久久久亚洲精品蜜臀av| 国产一区二区在线观看日韩| 午夜精品在线福利| 久久精品91蜜桃| 能在线免费看毛片的网站| 色哟哟哟哟哟哟| 可以在线观看的亚洲视频| 网址你懂的国产日韩在线| 欧美成人一区二区免费高清观看| 久久精品国产亚洲av涩爱 | 午夜爱爱视频在线播放| 国产精品乱码一区二三区的特点| 久久久久久久午夜电影| 国产欧美日韩精品一区二区| 一个人免费在线观看电影| 日日摸夜夜添夜夜添av毛片| 毛片一级片免费看久久久久| 亚洲色图av天堂| 久久九九热精品免费| 亚洲人成网站在线播放欧美日韩| 亚洲国产色片| 精品熟女少妇av免费看| 久久精品国产亚洲网站| 国产亚洲5aaaaa淫片| 晚上一个人看的免费电影| 国产午夜福利久久久久久| 18禁在线播放成人免费| 精品久久久久久久久久免费视频| 国产爱豆传媒在线观看| 色播亚洲综合网| 成年版毛片免费区| 青青草视频在线视频观看| 午夜福利成人在线免费观看| 91久久精品电影网| videossex国产| 色综合色国产| 国产一区亚洲一区在线观看| 啦啦啦啦在线视频资源| 男人舔奶头视频| 亚洲aⅴ乱码一区二区在线播放| 十八禁国产超污无遮挡网站| 久久精品久久久久久噜噜老黄 | 日韩强制内射视频| 亚洲欧美日韩高清专用| 国产精品麻豆人妻色哟哟久久 | 卡戴珊不雅视频在线播放| 国产探花极品一区二区| 男女边吃奶边做爰视频| 不卡视频在线观看欧美| 亚洲五月天丁香| .国产精品久久| 国产老妇伦熟女老妇高清| 91久久精品电影网| 国产v大片淫在线免费观看| 18禁黄网站禁片免费观看直播| 国产综合懂色| 亚洲国产欧美在线一区| 熟妇人妻久久中文字幕3abv| 岛国在线免费视频观看| 国产成人精品久久久久久| 国产麻豆成人av免费视频| 蜜臀久久99精品久久宅男| 精品99又大又爽又粗少妇毛片| 亚洲国产精品sss在线观看| 99久久成人亚洲精品观看| 久久久久久久久久黄片| 尾随美女入室| 国产黄片视频在线免费观看| 欧美三级亚洲精品| 麻豆成人午夜福利视频| 国产精品福利在线免费观看| 免费在线观看成人毛片| 国产亚洲5aaaaa淫片| 亚洲一区高清亚洲精品| 五月伊人婷婷丁香| 校园春色视频在线观看| 亚洲自偷自拍三级| 99热这里只有是精品50| 国产一区二区三区在线臀色熟女| 日本熟妇午夜| 久久久精品欧美日韩精品| 夫妻性生交免费视频一级片| 亚洲av免费高清在线观看| 国产精品无大码| 成熟少妇高潮喷水视频| 久久人人爽人人爽人人片va| 99久久精品一区二区三区| 噜噜噜噜噜久久久久久91| 黄色日韩在线| 天天躁日日操中文字幕| ponron亚洲| 久久精品国产亚洲网站| 国产精品一及| 国产成人午夜福利电影在线观看| 成人漫画全彩无遮挡| 日日摸夜夜添夜夜添av毛片| 亚洲av第一区精品v没综合| av福利片在线观看| 免费在线观看成人毛片| 久久国产乱子免费精品| 欧美一级a爱片免费观看看| 我要搜黄色片| 波野结衣二区三区在线| 少妇人妻精品综合一区二区 | 天美传媒精品一区二区| 天堂中文最新版在线下载 | 国产一级毛片在线| 亚洲国产高清在线一区二区三| 美女高潮的动态| 亚洲成人中文字幕在线播放| 午夜精品一区二区三区免费看| 欧美一区二区国产精品久久精品| 亚洲av二区三区四区| 国产精品久久视频播放| av在线观看视频网站免费| 内地一区二区视频在线| 亚洲国产欧美在线一区| www.av在线官网国产| 九九热线精品视视频播放| 免费av不卡在线播放| 一本久久中文字幕| 亚洲国产精品sss在线观看| 蜜桃久久精品国产亚洲av| 九九热线精品视视频播放| 天堂av国产一区二区熟女人妻| 欧美zozozo另类| 国产精品人妻久久久久久| 日本av手机在线免费观看| 简卡轻食公司| 国产精品美女特级片免费视频播放器| 久久精品91蜜桃| 欧美xxxx黑人xx丫x性爽| 1000部很黄的大片| 午夜精品在线福利| 三级经典国产精品| 亚洲四区av| 欧美日韩综合久久久久久| 午夜精品在线福利| 国内精品美女久久久久久| 亚洲七黄色美女视频| 少妇的逼水好多| 黑人高潮一二区| kizo精华| 国产一区二区三区在线臀色熟女| 日韩视频在线欧美| 国产视频首页在线观看| 欧美变态另类bdsm刘玥| 日本免费a在线| 欧美+日韩+精品| 能在线免费看毛片的网站| 色哟哟哟哟哟哟| 99热6这里只有精品| 欧美日韩一区二区视频在线观看视频在线 | 成年免费大片在线观看| 亚洲精品亚洲一区二区| 久久久久九九精品影院| 色5月婷婷丁香| 高清日韩中文字幕在线| 色尼玛亚洲综合影院| 菩萨蛮人人尽说江南好唐韦庄 | 国产91av在线免费观看| 日韩欧美精品v在线| 国产av麻豆久久久久久久| 一区二区三区免费毛片| 久久久久久九九精品二区国产| 精品人妻一区二区三区麻豆| 久久精品影院6| 在线a可以看的网站| 成人性生交大片免费视频hd| 日日撸夜夜添| 亚洲第一电影网av|