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

    微量生物碳酸鹽ICP-OES元素分析中的儀器優(yōu)化

    2011-01-11 14:41:26于心科常鳳鳴仇曉華南青云
    海洋科學 2011年6期
    關鍵詞:碳酸鹽譜線比值

    路 波, 于心科 常鳳鳴 仇曉華, 南青云

    (1. 中國科學院 海洋研究所 海洋地質與環(huán)境重點實驗室, 山東 青島 266071; 2. 中國科學院 研究生院, 北京 100049)

    微量生物碳酸鹽ICP-OES元素分析中的儀器優(yōu)化

    路 波1,2, 于心科1, 常鳳鳴1, 仇曉華1,2, 南青云1

    (1. 中國科學院 海洋研究所 海洋地質與環(huán)境重點實驗室, 山東 青島 266071; 2. 中國科學院 研究生院, 北京 100049)

    為優(yōu)化電感耦合等離子體發(fā)射光譜儀(ICP-OES)測試微量生物碳酸鹽中元素含量的工作條件, 在Thermo-Fisher公司生產的iCAP 6300 Radial型ICP-OES上, 對1 050組儀器條件、32種不同組成的標準溶液、36個譜線對和兩種儀器校正工作曲線進行了測試; 借助開源科學計算軟件Scilab數據插值和可視化技術對測試數據進行了比較分析。結果表明, 該儀器測試微量生物碳酸鹽中Mg/Ca和Sr/Ca比值的最佳條件為: 霧化器, 0.1 mL/min同心霧化器; 冷卻氣流量, 缺省設置12 L/min; 輔助氣流量, 0.5 L/min; 霧化氣壓力, 0.24 MPa; 蠕動泵速, 16 r/min; RF功率, 1 000 W; 觀察高度, 15 mm; 譜線, Ca 373.690 nm, Mg 279.553 nm {121}, Sr 407.771 nm; 儀器校準工作曲線為強度-濃度工作曲線。

    ICP-OES; Scilab; 生物碳酸鹽; 儀器優(yōu)化

    生物碳酸鹽(珊瑚、貝殼、有孔蟲、鈣質超微化石和介形蟲等)廣泛存在于海洋沉積物中, 其主要化學成分為CaCO3, 并含有微量Mg、Sr等元素。近年來, 生物碳酸鹽中的元素比值Mg/Ca、Sr/Ca已被廣泛應用于古海洋環(huán)境的重建[1-10]。在生物碳酸鹽微量元素分析中, 電感耦合等離子體發(fā)射光譜儀(ICPOES)因成本相對低廉、維護簡單、操作方便和可對多個元素快速同時測定等特點而被廣泛應用[2,5,11-14]。通常, 在對生物碳酸鹽進行 ICP-OES分析時, 所獲得的樣品量一般很少, 使得分析溶液中微量元素的濃度很低, 這對分析靈敏度提出了很高的要求; 同時由于生物碳酸鹽的樣品量常變化很大(以有孔蟲為例,常見固體樣品量約10~103μg), 使得分析溶液的濃度也變化不一, 形成明顯的Ca元素基體效應[2,15]。雖然de Villiers等[15]提出可通過使用強度比-濃度比工作曲線的方式消減, 但對此方法的實驗論證還很少,需進行更多的評估。

    除此之外, 由于 ICP-OES儀器間性能的差異和選擇測試項目的不同, 分析之前常需優(yōu)化儀器參數(霧化氣壓力、蠕動泵速、RF功率和觀察高度, 譜線選擇等)以求取得最佳分析結果[16-17]。但是, 針對生物碳酸鹽 ICP-OES分析, 如何優(yōu)化分析儀器參數的報道尚不多見。本文就儀器參數變化對分析性能的影響進行了實驗, 并利用開源科學運算軟件 Scilab三維數據插值和可視化技術進行了直觀展示, 確定了最佳儀器參數組合。同時, 對影響分析結果的譜線對組合、工作曲線類型也進行了考察。實驗結果對生物碳酸鹽的微量元素分析有重要幫助, 對其他類型樣品的分析也有重要參考價值。

    1 材料和方法

    1.1 試劑和溶液

    存放溶液的容器皆為 Nalgene公司生產的高密度聚乙烯(HDPE)瓶。使用前用 1 mol/L分析純硝酸浸泡1周以上, 重蒸餾水沖洗5遍, 去離子水(由重蒸餾水經 Pall?Purelab Classic型純水機制得, 電阻率≥18.2 MΩ)沖洗一遍, 60℃烘干或凍干備用。

    標準溶液均按重量法配置。先將高純固體藥品CaCO3、SrCO3、MnCO3、Fe2O3、Mg(NO3)2·6H2O、Fe(NO3)3·9H2O 和 Al(NO3)3·9H2O(美國 Sigma 公司生產)溶解于高純硝酸(質譜純級, 江蘇晶玻試劑有限公司)溶液中配置成高濃度單元素標液, 再按需要配置成不同組成的溶液系列(表1)。

    表1 實驗中用到的標準溶液Tab. 1 The standard solutions in experiments

    1.2 儀器介紹

    實驗中所用ICP-OES系賽默飛世(Thermo-Fisher)公司生產的iCAP 6300 Radial型。儀器為徑向觀測方式, 配備中階梯光柵和電荷注入式固體檢測器(CID),檢測波長為166 ~ 847 nm。27.12 MHz高頻固態(tài)發(fā)生器的功率在750 ~ 1 500 W范圍內連續(xù)可調, 蠕動泵在0 ~ 125 r /min范圍內連續(xù)可調。分析時選用Ar氣為工作氣體, 冷卻氣固定為12 L/min, 輔助氣0.5 L/min, 霧化器為GE公司生產的0.1 mL/min的微量玻璃同心霧化器。

    1.3 Scilab介紹

    Scilab軟件由法國國立信息與自動化研究院(INRIA)和法國國立橋梁學院(ENPC)開發(fā), 于 1994年開始向外界提供完全開放的源代碼, 于2001年起由中法聯(lián)合實驗室在中國推廣, 具有矩陣運算、數值運算、程序設計和豐富的繪圖功能等。本研究使用該軟件完成了儀器參數優(yōu)化環(huán)節(jié)中對測試數據的數據插值和三維圖示。

    2 結果與討論

    2.1 儀器參數優(yōu)化

    在利用 ICP-OES分析時, 霧化氣壓力(NP)、觀察高度(VH)、高頻發(fā)生器功率(RF)和蠕動泵速(PS)的變化對分析性能影響顯著, 可以用譜線強度(INT)和信背比(IBR)來考察。

    INT指譜線的強度, INT越高, 說明對分析元素濃度變化的響應越靈敏, 檢測限也低; 但 INT太高也意味著對其他譜線的干擾增強, 有時過高將出現(xiàn)溢出檢測器檢測能力而使分析失敗的現(xiàn)象, 即信號過飽和問題。IBR為譜線強度和背景強度的比值, 值越高, 則背景相對影響越小, 越有利于分析的準確性和精確性; 反之, 則分析的準確性和精確性越低。

    霧化氣在霧化器處和分析溶液一起形成氣溶膠,隨后經霧室、矩管進入等離子體。當NP增大時, 單位時間內被攜帶進入等離子體的氣溶膠增多, 分析元素相應譜線的INT增強, IBR升高。但當NP過高時, 氣溶膠將很快穿過等離子體, 會因停留時間過短來不及被充分激發(fā); 另一方面則因氣溶膠的冷卻作用進一步降低被激發(fā)的比例, 從而造成譜線 INT升高到一定峰值后下降, 但是IBR卻會一直增大。

    VH為分析譜線在等離子體上的采集位置。隨VH增加, 譜線INT將減小, 但IBR增大。

    高頻發(fā)生器為等離子體提供能量。RF越大, 元素被激發(fā)的比例越高, 譜線 INT越高; 但是增高到一定值后, INT隨RF增加而增加的幅度將逐漸減小,直至幾無變化。同時由于干擾譜線強度增加、背景輻射增強等原因, IBR將持續(xù)減小。

    蠕動泵將分析溶液提升至霧化器處, 當PS增大時, 溶液提升量增加, 氣溶膠中液滴濃度將增加, 此時的譜線INT增強, IBR升高。但PS過大, 霧化效率降低, 氣溶膠對等離子體的冷卻作用增強。此時,INT、IBR隨PS升高的幅度將減小。如果繼續(xù)增大PS, INT、IBR將不增反降。

    為考查NP、VH、RF、PS對INT和INT、IBR的影響, 并針對生物碳酸鹽的元素組成特點, 配置了 Ca濃度為16.20 μg/g、Mg濃度19.39 ng/g、Sr濃度45.04 ng/g的標準溶液(33)(表1)考察了不同NP、VH、RF、PS條件下INT和IBR的變化情況。實驗中選擇了譜線Mg27951(279.553 nm, 121級次)、Sr4077(407.771 nm)[11]和 Ca3736(373.693 nm)為分析譜線, 共考察了NP在0.16~0.24 MPa之間、PS在8 ~ 24 r/min之間、RF在1 000~1 300 W之間、VH在9~19 mm之間的儀器參數組合, 共計1 050(5×5×7×6)種(表2)。

    表2 ICP-OES優(yōu)化中部分參數考察范圍Tab. 2 Range of experimental parameters in ICP-OES Optimization

    對應每一條譜線, 得到1 050個INT(IBR)值, 先按NP取值分為5組, 每組包括210(5×7×6)個PS、RF和VH組合下對應的INT(IBR)值。在Scilab中以PS、RF、VH為坐標軸構建三維空間, 將NP分組中每組210個INT(IBR)值在此空間中插值求出一定間隔等值點的 PS、RF、VH坐標, 并連接成面得到在PS、RF、VH三維空間中INT(IBR)等值面的分布(源代碼修改自 Roberts[18])(圖 1)。INT等值面以紅色曲面表示, 紅色越深, INT越高。IBR等值面以藍色曲面表示, 顏色越深, IBR越高。

    圖1 IBR和INT等值面隨NP、PS、RF和VH變化的情況Fig. 1 Isosurfaces of IBR and INT relative to varying NP, PS, RF and VH

    Ca3736、Mg27951和Sr4077的INT、IBR等值面(圖 1)與 PS、RF、VH 三個坐標軸都相交, 說明INT(IBR)受多個儀器參數綜合作用的影響。但是Ca3736、Mg27951和Sr4077的INT(IBR)等值面形狀和位置各不相同, 說明NP、PS、RF和VH對不同元素譜線的影響也不相同。

    由于生物碳酸鹽樣品以碳酸鈣為主要組成, 故元素Ca成為分析時的基體元素, 其INT和IBR值總是很高。為降低其基體效應, 應選擇 Ca譜線的INT較低的儀器參數組合。而 Mg、Sr等元素含量低, INT和IBR值總是很低。為了分析結果的準確和精確, 在Mg、Sr譜線的INT高于一定值時, IBR應越大越好。此外, 由于實際分析中生物碳酸鹽樣品量微少, 應控制分析時溶液消耗量, 即PS值越低越好。對元素組成接近生物碳酸鹽的標準溶液(33)來說, 當NP為0.24 MPa, VH為15 mm, RF為1 000 W, PS為16 r/min時, 連續(xù)3次曝光測量的溶液消耗不超過1 mL。同時Ca3736譜線INT為13 179.1 cts/S,IBR為18.7; Mg27951譜線INT為412.6 cts/s, IBR為6.6; Sr407譜線INT為4018.8 cts/s, IBR為6.4, 滿足微量碳酸鹽樣品的分析要求, 因此確定為最佳儀器參數。

    2.2 譜線選擇

    由于生物碳酸鹽樣品的組成和樣品量多變, 進行ICP-OES分析時的溶液組成常常變化不一。為了考察溶液組成變化對不同譜線的影響, 選擇出受組成變化影響最小的譜線組合, 對組成不同的32種溶液(Blank、溶液(1)~(31), 見表1)進行了測定和比較。其中, (1)~(11)Mg/Ca、Sr/Ca 不變, 但 Ca、Mg 和 Sr濃度連續(xù)增加; (12)~(21)Ca濃度變化不大, Mg、Sr濃度逐漸增加; (22)~(31)Mg、Sr濃度變化不大, Ca濃度逐漸增加。

    實際應用中, 分析結果多表達為元素與Ca的比值(如珊瑚 Sr/Ca, 有孔蟲、介形蟲 Mg/Ca)。為了便于比較, 本實驗選擇了以譜線對靈敏度比值為衡量指標, 并對其進行了標準化(表3)。

    表3 標準化的譜線對靈敏度比值Tab. 3 Normalized sensitivity ratios of element pairs

    譜線靈敏度為元素譜線強度與溶液濃度比值,譜線對靈敏度比值為相應譜線的靈敏度比值。以接近生物碳酸鹽樣品平均組成的標準溶液(5)為基準,將其余溶液的譜線對靈敏度比值與溶液(5)作商, 即得標準化后的譜線對靈敏度比值。分別考察了譜線Ca3158、Ca3179(317.933 nm)、Ca3706(370.603 nm)、Ca4226(422.673 nm)、Ca4454(445.478 nm)、Mg27950(279.553 nm, 120級次)、Mg27951、Mg2802(280.270 nm)、Mg2852(285.213 nm)、Sr4077 和 Sr4215(421.552 nm)。對溶液(1)~(31)的統(tǒng)計顯示(表 3), 譜線對Mg27951/Ca3736、Sr4077/Ca3736靈敏度比值隨溶液組成變化的相對標準偏差最小, 分別為 0.95 % (1σ)和 2.47 % (1σ), 說明譜線 Ca3736、Mg27951和Sr4077受溶液變化的差異影響最小, 可以選擇為分析用譜線。

    2.3 工作曲線比較

    使用ICP-OES分析碳酸鹽元素組成時, Ca的基體效應一直受到關注[2,15]。de Villers等[15]提出可以通過建立強度比-濃度比工作曲線(即以元素濃度比值為橫坐標、以相應譜線對的強度比值為縱坐標在笛卡爾坐標系中通過最小二乘法擬合建立的工作曲線)的方式消減。為了評估其適用性, 選用了有孔蟲Mg/Ca、Sr/Ca分析中工作溶液(62)前后 4個月的測試數據, 共計12批次, 1021個測量值, 將強度比-濃度比工作曲線與傳統(tǒng)強度-濃度工作曲線(即以元素濃度為橫坐標、以相應譜線的強度為縱坐標在笛卡爾坐標系中通過最小二乘法擬合建立的工作曲線)下的Mg/Ca、Sr/Ca計算結果進行了比較。

    其中, 強度-濃度工作曲線根據標準溶液(52)、(63)、(54)~(56)建立, 強度比-濃度比工作曲線根據標準溶液(58)~(61)建立。對于 Mg/Ca, 強度-濃度工作曲線下計算結果的相對誤差為 0.604 %, 精確度為0.515 %(1σ RSD); 強度比-濃度比工作曲線下相對誤差為-0.937 %, 精確度1.035 %(1σ RSD)。對于Sr/Ca分析, 強度-濃度工作曲線下計算結果的相對誤差為-0.924 %, 精確度為 0.818 %(1σ RSD); 強度比-濃度比工作曲線下相對誤差為-1.680 %, 精確度0.898 %(1σ RSD)。這說明在 iCAP 6300 Radial平臺上, 建立強度比-濃度比工作曲線并不能提高分析的準確度和精確度。

    3 結論

    在iCAP 6300 Radial型ICP-OES上, 通過考察不同儀器參數、不同譜線和工作曲線對分析結果的影響, 可以得出如下結論:

    (1) 元素譜線強度和信背比受多個儀器參數綜合影響; 選擇儀器參數時, 可利用Scilab的數據插值和可視化技術方便地確定儀器參數組合。利用iCAP 6300 Radial型ICP-OES分析生物碳酸鹽的Mg/Ca、Sr/Ca時, 最佳儀器參數為: 霧化器, 0.1 mL/min同心霧化器; 冷卻氣流量, 缺省設置12 L/min; 輔助氣流量, 0.5 L/min; 霧化氣壓力, 0.24 MPa; 蠕動泵速, 16 r/min, 高頻發(fā)生器功率, 1 000 W; 觀察高度, 15 mm。

    (2) 進行元素比值分析時, 可利用標準化后的元素譜線對靈敏度比值衡量溶液組成變化對分析結果的影響, 確定最佳元素譜線組合。利用iCAP 6300 Radial型ICP-OES分析生物碳酸鹽的Mg/Ca、Sr/Ca時, 譜線Ca 373.690 nm、Mg 279.553 nm {121}和Sr 407.771 nm是最佳譜線組合。

    (3) 進行元素比值分析時, 強度比-濃度比工作曲線下計算結果并不總是優(yōu)于傳統(tǒng)強度-濃度工作曲線下計算結果。利用iCAP 6300 Radial型ICP-OES分析生物碳酸鹽的Mg/Ca、Sr/Ca時, 最佳工作曲線為強度-濃度工作曲線。

    [1]Lea D W, Pak D K, Spero H J. Climate impact of late quaternary equatorial Pacific sea surface temperature variations[J]. Science, 2000, 289(5485): 1719-1724.

    [2]Schrag D P. Rapid analysis of high-precision Sr/Ca ratios in corals and other marine carbonates[J].Paleoceanography, 1999, 14(2): 97-102.

    [3]Mucci A. Influence of temperature on the composition of magnesian calcite overgrowths precipitated from seawater[J]. Geochimica et Cosmochimica Acta, 1987,51(7): 1977-1984.

    [4]Katz A. The interaction of magnesium with calcite during crystal growth at 25-90°C and one atmosphere[J]. Geochimica et Cosmochimica Acta,1973, 37(6): 1563-1578.

    [5]Mitsuguchi T, Matsumoto E, Abe O, et al. Mg/Ca thermometry in coral-skeletons[J]. Science, 1996, 274(5289): 961-963.

    [6]Correge T. Preliminary results of paleotemperature reconstruction using the magnesium to calcium ratio of deep-sea ostracode shells from the late Quaternary of site 822, leg 133 (western Coral Sea)[J]. 1993, 133:175-180.

    [7]Dwyer G S, Cronin T M, Baker P A, et al. North Atlantic Deepwater Temperature Change During Late Pliocene and Late Quaternary Climatic Cycles[J]. Science,1995, 270(5240): 1347-1351.

    [8]Nürnberg D, Bijma J, Hemleben C. Assessing the reliability of magnesium in foraminiferal calcite as a proxy for water mass temperatures[J]. Geochimica et Cosmochimica Acta, 1996, 60(5): 803-814.

    [9]Lea D W, Mashiotta T A, Spero H J. Controls on magnesium and strontium uptake in planktonic foraminifera determined by live culturing[J]. Geochimica Et Cosmochimica Acta, 1999, 63(16): 2369-2379.

    [10]Elderfield H, Ganssen G. Past temperature and delta O-18 of surface ocean waters inferred from foraminiferal Mg/Ca ratios[J]. Nature, 2000, 405(6785):442-445.

    [11]黃超, 路波, 于心科. ICP-OES直接測定海洋沉積物中有孔蟲的 Mg和 Ca[J]. 光譜實驗室, 2008, 25(4):727-731.

    [12]Andreasen D H, Sosdian S, Perron-Cashman S, et al.Fidelity of radially viewed ICP-OES and magnetic-sector ICP-MS measurement of Mg/Ca and Sr/Ca ratios in marine biogenic carbonates: Are they trustworthy together?[J]. Geochemistry Geophysics Geosystems, 2006, 7.

    [13]韋剛健, 余克服, 李獻華, 等. 南海北部珊瑚 Sr/Ca和Mg/Ca溫度計及高分辨率SST記錄重建嘗試[J].第四紀研究, 2004, 24(3): 325-331.

    [14]鄧文峰, 韋剛健, 李獻華. 有孔蟲的高精度 Mg/Ca比值的 ICP—AES分析[J]. 海洋地質與第四紀地質,2005, 25(2): 147 - 151.

    [15]de Villiers S, Greaves M, Elderfield H. An intensity ratio calibration method for the accurate determination of Mg/Ca and Sr/Ca of marine carbonates by ICP-AES[J]. Geochemistry Geophysics Geosystems,2002, 3

    [16]Sun Y L, Sun M. Inductively coupled plasma atomic emission spectrometry: An efficient tool for precise determination of Mg/Ca and Sr/Ca ratios in biocarbonates[J]. Applied Spectroscopy, 2003, 57(6):711-714.

    [17]Green D R H, Cooper M J, German C R, et al.Optimization of an inductively coupled plasma-optical emission spectrometry method for the rapid determination of high-precision Mg/Ca and Sr/Ca in foraminiferal calcite[J]. Geochemistry Geophysics Geosystems, 2003, 4

    [18]Roberts T. IMAGE TOOLS-contour3d.sci. www.scilab.org/contrib/index_contrib.php?page=display-Contributio n&fileID=1236, 2009.

    Optimization of ICP-OES for analysis of micro-mass biogenic carbonate

    LU Bo1,2, YU Xin-ke1, CHANG Feng-ming1, QIU Xiao-hua1,2, NAN Qing-yun1
    (1. Key Laboratory of Marine Geology and Environment, Institute of Oceanology, Chinese Academy of Sciences, Qingdao 266071, China; 2. Graduate University of Chinese Academy of Sciences, Beijing 100049,China)

    Jun., 10, 2010

    ICP-OES; Scilab; Biogenic carbonate; Instrument optimization

    To optimize Inductively Coupled Plasma-Optical Emission Spectrometer (ICP-OES) for elemental analysis of micro-mass biogenic carbonate, 1050 groups of instrument parameters, 32 standard solutions with different constituents, 36 analyte pairs and two instrument calibration methods were tested on a Thermo-Fisher iCAP 6300 Radial ICP-OES. The data were analyzed by means of conventional statistics and Scilab. The optimized instrument conditions of Thermo-fisher iCAP 6300 Radial ICP-OES for analysis of Mg/Ca and Sr/Ca ratios in micro-mass biogenic carbonate were enlisted as the follows: glass concentric nebulizer, 0.1 ml/min; coolant gas, default 12 L/min; nebulizer gas pressure, 0.24 MPa; peristaltic pump speed, 16 rpm; radio frequency generator power,1000 w; view height, 15 mm; analyte lines, Ca 373.690 nm, Mg 279.553 nm {121}, Sr 407.771 nm; instrument calibration method, conventional intensity-concentration method.

    P736.4

    A

    1000-3096(2011)06-0001-07

    2010-06-10;

    2010-12-18

    國家重點基礎研究發(fā)展計劃“973”項目(2007CB815903);中國科學院知識創(chuàng)新工程項目(KZCX2-YW-221);中國科學院海洋研究所知識創(chuàng)新工程青年人才領域前沿項目(No.2007-10);國家自然科學基金青年科學基金項目(40906038); 海洋地質國家重點實驗室開放基金項目(MG0903)

    路波(1982-), 男, 江蘇沛縣人, 博士研究生, 主要從事地球化學和古海洋學研究 ,電話 : 0411-84783433, E-mail: luboemail@gmail.com; 于心科, 通信作者, 電話: 0532-82898533, E-mail:xyu@ms.qdio.ac.cn

    致謝:感謝孫晗杰、陳雙喜、熊志方博士在實驗過程中的熱情幫助, 感謝Thermo-Fisher公司工程師一直以來的技術支持和有益建議。文章修改過程中, 幾位審稿人提出了熱情中肯的寶貴意見, 在此一并感謝。

    康亦兼)

    猜你喜歡
    碳酸鹽譜線比值
    基于HITRAN光譜數據庫的合并譜線測溫仿真研究
    鐵合金光譜譜線分離實驗研究
    電子測試(2018年11期)2018-06-26 05:56:00
    比值遙感蝕變信息提取及閾值確定(插圖)
    河北遙感(2017年2期)2017-08-07 14:49:00
    鍶原子光鐘鐘躍遷譜線探測中的程序控制
    不同應變率比值計算方法在甲狀腺惡性腫瘤診斷中的應用
    海相碳酸鹽烴源巖生烴潛力模糊評價方法
    藥芯焊絲GMAW電弧光譜的研究
    機械工程師(2014年6期)2014-12-23 06:44:22
    超高壓均質聯(lián)合二甲基二碳酸鹽對荔枝汁中污染菌及其微生物貨架期的影響
    碳酸鹽型滑石物相分析方法研究
    雙電機比值聯(lián)動控制系統(tǒng)
    特大巨黑吊av在线直播| 人妻制服诱惑在线中文字幕| 免费av不卡在线播放| 在线观看66精品国产| 1000部很黄的大片| 少妇的逼好多水| 午夜福利视频1000在线观看| 十八禁国产超污无遮挡网站| 午夜两性在线视频| 香蕉av资源在线| АⅤ资源中文在线天堂| 黄色视频,在线免费观看| 久9热在线精品视频| 久久精品国产亚洲av天美| 国产欧美日韩精品一区二区| 麻豆久久精品国产亚洲av| 亚洲综合色惰| 久久久久九九精品影院| 一进一出抽搐动态| 日韩欧美在线乱码| 最好的美女福利视频网| 女生性感内裤真人,穿戴方法视频| 色哟哟·www| 国产国拍精品亚洲av在线观看| 亚洲国产欧美人成| 亚洲性夜色夜夜综合| 日韩欧美免费精品| 丝袜美腿在线中文| 日日干狠狠操夜夜爽| 欧美成人一区二区免费高清观看| 给我免费播放毛片高清在线观看| 国产中年淑女户外野战色| 欧美性感艳星| 岛国在线免费视频观看| 男插女下体视频免费在线播放| 欧美日韩福利视频一区二区| 97热精品久久久久久| 国产精品美女特级片免费视频播放器| 男人狂女人下面高潮的视频| 真实男女啪啪啪动态图| 国产午夜精品论理片| 看十八女毛片水多多多| 日韩欧美一区二区三区在线观看| 国产亚洲av嫩草精品影院| 日韩欧美国产在线观看| 国产免费一级a男人的天堂| 97超视频在线观看视频| 亚洲最大成人手机在线| 麻豆一二三区av精品| 12—13女人毛片做爰片一| 精品午夜福利在线看| 禁无遮挡网站| 国产色婷婷99| 97碰自拍视频| 日本免费a在线| 久久99热6这里只有精品| 成人精品一区二区免费| 日本与韩国留学比较| 欧美性感艳星| 精品午夜福利视频在线观看一区| 在线a可以看的网站| 久久精品夜夜夜夜夜久久蜜豆| 91在线精品国自产拍蜜月| 老司机午夜福利在线观看视频| 成人永久免费在线观看视频| 日韩欧美精品v在线| 国产精品女同一区二区软件 | 桃色一区二区三区在线观看| 老司机深夜福利视频在线观看| 久久国产乱子伦精品免费另类| 免费看美女性在线毛片视频| 亚洲男人的天堂狠狠| 黄色一级大片看看| 国产欧美日韩一区二区三| 欧美xxxx黑人xx丫x性爽| 亚洲第一电影网av| 两个人的视频大全免费| 极品教师在线免费播放| 欧美另类亚洲清纯唯美| 国产白丝娇喘喷水9色精品| 搡女人真爽免费视频火全软件 | 日韩欧美一区二区三区在线观看| 久久久国产成人免费| 免费人成在线观看视频色| 在线观看av片永久免费下载| 一个人免费在线观看电影| 色视频www国产| 成人永久免费在线观看视频| 日韩欧美一区二区三区在线观看| 简卡轻食公司| 亚洲男人的天堂狠狠| 亚洲专区国产一区二区| 九色成人免费人妻av| 能在线免费观看的黄片| 99久久99久久久精品蜜桃| 国产男靠女视频免费网站| 欧美性猛交╳xxx乱大交人| 99久久精品国产亚洲精品| 日日摸夜夜添夜夜添av毛片 | 男女之事视频高清在线观看| 国产av一区在线观看免费| 午夜两性在线视频| 亚洲最大成人av| 久久久国产成人精品二区| 少妇熟女aⅴ在线视频| 精品国产三级普通话版| 别揉我奶头~嗯~啊~动态视频| 精品久久久久久久人妻蜜臀av| 此物有八面人人有两片| 观看美女的网站| 欧美成人免费av一区二区三区| 成人国产综合亚洲| 亚洲av成人av| 国产成人欧美在线观看| 久久国产精品影院| 欧美精品国产亚洲| 亚洲成a人片在线一区二区| 熟妇人妻久久中文字幕3abv| 色吧在线观看| 最近在线观看免费完整版| 亚洲av中文字字幕乱码综合| 午夜福利欧美成人| 好男人在线观看高清免费视频| 成人欧美大片| 18禁黄网站禁片免费观看直播| 日韩av在线大香蕉| 国产成人a区在线观看| 日韩欧美 国产精品| 九九热线精品视视频播放| 欧美3d第一页| 亚洲av二区三区四区| 在线免费观看的www视频| 亚洲国产色片| 搡老熟女国产l中国老女人| 亚洲成人久久性| 99视频精品全部免费 在线| 精品国内亚洲2022精品成人| 国产黄片美女视频| 此物有八面人人有两片| 精品乱码久久久久久99久播| 午夜免费男女啪啪视频观看 | 91麻豆精品激情在线观看国产| 女生性感内裤真人,穿戴方法视频| 最后的刺客免费高清国语| 亚洲av电影不卡..在线观看| 欧美高清成人免费视频www| 日本黄大片高清| 久久久久久久久大av| 日本a在线网址| 能在线免费观看的黄片| 淫妇啪啪啪对白视频| .国产精品久久| 18禁黄网站禁片免费观看直播| 极品教师在线免费播放| aaaaa片日本免费| 日日摸夜夜添夜夜添小说| 人妻夜夜爽99麻豆av| 日韩欧美一区二区三区在线观看| 黄片小视频在线播放| 毛片女人毛片| 欧美绝顶高潮抽搐喷水| 国产91精品成人一区二区三区| 国产人妻一区二区三区在| 久久久久免费精品人妻一区二区| 久久久久久久久久黄片| 成年女人毛片免费观看观看9| av天堂中文字幕网| 极品教师在线免费播放| 又黄又爽又刺激的免费视频.| 神马国产精品三级电影在线观看| 欧美成人一区二区免费高清观看| 成人永久免费在线观看视频| 国产av不卡久久| 夜夜夜夜夜久久久久| 国产黄片美女视频| 欧美一区二区亚洲| 人妻丰满熟妇av一区二区三区| 欧美又色又爽又黄视频| 久久性视频一级片| 国产亚洲精品av在线| 久久亚洲精品不卡| 99久久精品国产亚洲精品| 波野结衣二区三区在线| 国产精品自产拍在线观看55亚洲| 又爽又黄a免费视频| 日本免费一区二区三区高清不卡| 亚洲五月天丁香| 99热精品在线国产| www.色视频.com| 99精品久久久久人妻精品| 午夜日韩欧美国产| 日韩欧美三级三区| 亚洲男人的天堂狠狠| 国产高清视频在线播放一区| 欧美成人a在线观看| 国产精品电影一区二区三区| 精品一区二区免费观看| 又爽又黄无遮挡网站| 夜夜看夜夜爽夜夜摸| 乱人视频在线观看| 国产色婷婷99| 久久久成人免费电影| 在线免费观看不下载黄p国产 | 九九久久精品国产亚洲av麻豆| 精品一区二区三区视频在线| 婷婷精品国产亚洲av| 国产亚洲欧美在线一区二区| 九色国产91popny在线| 亚洲乱码一区二区免费版| 亚洲欧美清纯卡通| 亚洲国产日韩欧美精品在线观看| 男女床上黄色一级片免费看| 亚洲久久久久久中文字幕| 麻豆久久精品国产亚洲av| а√天堂www在线а√下载| 亚洲在线观看片| av国产免费在线观看| 精品欧美国产一区二区三| 99国产精品一区二区三区| 亚洲精品色激情综合| 国产单亲对白刺激| 一区二区三区激情视频| 88av欧美| 久久久久久久精品吃奶| 黄色女人牲交| 午夜福利视频1000在线观看| 人人妻人人看人人澡| 美女免费视频网站| 久久天躁狠狠躁夜夜2o2o| 尤物成人国产欧美一区二区三区| 成人国产综合亚洲| 免费电影在线观看免费观看| 国产人妻一区二区三区在| 搞女人的毛片| 久久久久久久久大av| 可以在线观看的亚洲视频| 搡老熟女国产l中国老女人| 五月伊人婷婷丁香| 18禁裸乳无遮挡免费网站照片| 亚洲欧美日韩高清在线视频| 此物有八面人人有两片| 国产精品女同一区二区软件 | 69av精品久久久久久| 内地一区二区视频在线| 三级男女做爰猛烈吃奶摸视频| 免费av毛片视频| 男插女下体视频免费在线播放| 午夜精品一区二区三区免费看| 此物有八面人人有两片| 两个人视频免费观看高清| 色尼玛亚洲综合影院| 亚洲国产高清在线一区二区三| 久久精品综合一区二区三区| 国产一区二区在线av高清观看| 天天躁日日操中文字幕| 国产综合懂色| 国产三级在线视频| 男女之事视频高清在线观看| 99热这里只有是精品50| 午夜免费成人在线视频| 亚洲av一区综合| 尤物成人国产欧美一区二区三区| 免费在线观看影片大全网站| 俄罗斯特黄特色一大片| 少妇人妻精品综合一区二区 | 国产久久久一区二区三区| 精品久久久久久久久久免费视频| 成人毛片a级毛片在线播放| 精品久久久久久久久av| 免费在线观看亚洲国产| 9191精品国产免费久久| 国产av在哪里看| 欧美zozozo另类| 国产aⅴ精品一区二区三区波| 国产精品一区二区性色av| 九九热线精品视视频播放| 欧美中文日本在线观看视频| 日韩国内少妇激情av| 搡女人真爽免费视频火全软件 | 伊人久久精品亚洲午夜| 国产精品av视频在线免费观看| 久久天躁狠狠躁夜夜2o2o| 天堂网av新在线| 欧美日本亚洲视频在线播放| 免费搜索国产男女视频| 日韩大尺度精品在线看网址| 亚洲成人精品中文字幕电影| 99久久无色码亚洲精品果冻| 9191精品国产免费久久| 欧美一区二区国产精品久久精品| 欧美潮喷喷水| 天堂动漫精品| 99热这里只有是精品50| 99热只有精品国产| 国产在视频线在精品| 亚洲精品日韩av片在线观看| 中文字幕精品亚洲无线码一区| 欧美三级亚洲精品| 午夜老司机福利剧场| 天天一区二区日本电影三级| 搡老妇女老女人老熟妇| 韩国av一区二区三区四区| 每晚都被弄得嗷嗷叫到高潮| 亚洲成av人片免费观看| 永久网站在线| 男插女下体视频免费在线播放| 老熟妇乱子伦视频在线观看| 色播亚洲综合网| 熟妇人妻久久中文字幕3abv| 观看美女的网站| 18+在线观看网站| www.999成人在线观看| 亚洲国产精品久久男人天堂| 搡女人真爽免费视频火全软件 | 国产精品久久久久久人妻精品电影| 亚洲午夜理论影院| 久久久久久大精品| 看黄色毛片网站| 少妇的逼好多水| 99热这里只有是精品50| 亚洲国产欧洲综合997久久,| 人妻丰满熟妇av一区二区三区| 国产精品一及| av在线老鸭窝| 久久99热这里只有精品18| 又爽又黄无遮挡网站| av中文乱码字幕在线| 精品熟女少妇八av免费久了| 日韩人妻高清精品专区| 亚洲成人久久爱视频| 1000部很黄的大片| 精品熟女少妇八av免费久了| 国产乱人视频| 精品人妻熟女av久视频| 亚洲精品亚洲一区二区| 精品无人区乱码1区二区| 12—13女人毛片做爰片一| 成年人黄色毛片网站| 小说图片视频综合网站| 一本久久中文字幕| 999久久久精品免费观看国产| 午夜老司机福利剧场| 亚洲精品一卡2卡三卡4卡5卡| 欧美成人性av电影在线观看| 757午夜福利合集在线观看| 免费看美女性在线毛片视频| 久久精品国产亚洲av涩爱 | 欧美日韩综合久久久久久 | 狂野欧美白嫩少妇大欣赏| 麻豆久久精品国产亚洲av| 给我免费播放毛片高清在线观看| 夜夜看夜夜爽夜夜摸| 最后的刺客免费高清国语| 久久久久国产精品人妻aⅴ院| 免费看美女性在线毛片视频| 伊人久久精品亚洲午夜| 亚洲国产精品成人综合色| 熟妇人妻久久中文字幕3abv| 亚洲av熟女| 伊人久久精品亚洲午夜| 色播亚洲综合网| 欧美不卡视频在线免费观看| 欧美成人一区二区免费高清观看| 亚洲国产高清在线一区二区三| 99riav亚洲国产免费| 精品不卡国产一区二区三区| 午夜免费激情av| 欧美日韩综合久久久久久 | 51国产日韩欧美| 成人毛片a级毛片在线播放| 99精品久久久久人妻精品| 日韩欧美 国产精品| 成人av在线播放网站| 麻豆久久精品国产亚洲av| 真人一进一出gif抽搐免费| 动漫黄色视频在线观看| 99在线视频只有这里精品首页| 人妻夜夜爽99麻豆av| 麻豆成人午夜福利视频| 久久精品人妻少妇| 一级黄片播放器| av在线观看视频网站免费| 老司机午夜十八禁免费视频| 精品久久久久久,| 欧美黄色淫秽网站| 亚洲成人中文字幕在线播放| 亚洲美女搞黄在线观看 | 免费看a级黄色片| 亚洲成a人片在线一区二区| 最好的美女福利视频网| 极品教师在线视频| 草草在线视频免费看| 色av中文字幕| 久久久成人免费电影| 精品人妻一区二区三区麻豆 | 麻豆成人av在线观看| 高清毛片免费观看视频网站| 999久久久精品免费观看国产| 日韩亚洲欧美综合| 国产精品久久久久久亚洲av鲁大| xxxwww97欧美| 日韩欧美精品v在线| 男人狂女人下面高潮的视频| 国产精品av视频在线免费观看| 精品人妻一区二区三区麻豆 | 麻豆久久精品国产亚洲av| 国产人妻一区二区三区在| 97碰自拍视频| 国产成人欧美在线观看| 久久精品国产亚洲av香蕉五月| 天堂影院成人在线观看| 精品午夜福利在线看| 久久久精品大字幕| 一个人免费在线观看的高清视频| 久久久久久久久久成人| 欧美在线一区亚洲| 亚洲中文字幕日韩| 国产淫片久久久久久久久 | 午夜福利视频1000在线观看| ponron亚洲| 18禁黄网站禁片午夜丰满| 色哟哟哟哟哟哟| 简卡轻食公司| 色在线成人网| 国产精品久久久久久久久免 | 成人av在线播放网站| 美女高潮的动态| 日本免费一区二区三区高清不卡| 成年免费大片在线观看| 国产主播在线观看一区二区| 国产高清三级在线| 欧美日韩中文字幕国产精品一区二区三区| 亚洲无线观看免费| 国产精品久久久久久人妻精品电影| 在线a可以看的网站| 亚洲av免费高清在线观看| 大型黄色视频在线免费观看| 国产精品久久久久久久久免 | 亚洲一区高清亚洲精品| 国产高清激情床上av| 久久国产精品人妻蜜桃| 色综合站精品国产| 欧美精品啪啪一区二区三区| 高清毛片免费观看视频网站| 好男人电影高清在线观看| xxxwww97欧美| 中文字幕人成人乱码亚洲影| 亚洲天堂国产精品一区在线| 国产一区二区亚洲精品在线观看| 欧美激情国产日韩精品一区| 久久午夜福利片| 国产精品一区二区性色av| x7x7x7水蜜桃| 日韩中字成人| 精品久久久久久久久久免费视频| or卡值多少钱| 国产乱人伦免费视频| 男插女下体视频免费在线播放| 国内少妇人妻偷人精品xxx网站| 99热这里只有是精品在线观看 | 看免费av毛片| 美女xxoo啪啪120秒动态图 | 夜夜爽天天搞| 天天一区二区日本电影三级| 一个人免费在线观看电影| 在线十欧美十亚洲十日本专区| 国内精品一区二区在线观看| 国产高潮美女av| 18禁黄网站禁片午夜丰满| 亚洲av成人不卡在线观看播放网| 成人性生交大片免费视频hd| 三级国产精品欧美在线观看| 欧美不卡视频在线免费观看| 亚洲无线在线观看| 亚洲在线自拍视频| 丰满的人妻完整版| 欧美另类亚洲清纯唯美| 最新中文字幕久久久久| 欧美高清性xxxxhd video| 国内揄拍国产精品人妻在线| 极品教师在线视频| 亚洲精品在线观看二区| 一进一出抽搐动态| 大型黄色视频在线免费观看| 99久久无色码亚洲精品果冻| 日本成人三级电影网站| 性色avwww在线观看| 欧美三级亚洲精品| 国产成人影院久久av| 女同久久另类99精品国产91| 一个人免费在线观看电影| 99国产精品一区二区三区| 极品教师在线免费播放| 小蜜桃在线观看免费完整版高清| 久久久成人免费电影| 中文字幕免费在线视频6| 少妇被粗大猛烈的视频| 国产真实乱freesex| 三级国产精品欧美在线观看| 少妇人妻精品综合一区二区 | 天堂影院成人在线观看| 久久99热这里只有精品18| 精品99又大又爽又粗少妇毛片 | 欧美精品国产亚洲| 99热这里只有是精品50| www.999成人在线观看| 欧美最黄视频在线播放免费| 亚洲狠狠婷婷综合久久图片| 久久久久久久久中文| 欧美性猛交╳xxx乱大交人| 久久久久久久久久成人| 免费观看人在逋| 久久精品国产亚洲av涩爱 | 亚洲avbb在线观看| 久久亚洲精品不卡| 国产成人啪精品午夜网站| 色吧在线观看| 在线播放国产精品三级| 黄色女人牲交| 中文字幕人妻熟人妻熟丝袜美| 国产主播在线观看一区二区| 亚洲人成网站在线播| 黄色女人牲交| 尤物成人国产欧美一区二区三区| 91麻豆精品激情在线观看国产| 久久久精品大字幕| 欧美精品国产亚洲| 国产伦精品一区二区三区四那| 久久精品综合一区二区三区| 亚洲三级黄色毛片| 有码 亚洲区| 亚洲成av人片在线播放无| 黄色视频,在线免费观看| 熟妇人妻久久中文字幕3abv| 两个人视频免费观看高清| 成人国产综合亚洲| 国产精品久久久久久人妻精品电影| 在线观看免费视频日本深夜| 国产精华一区二区三区| 免费在线观看影片大全网站| bbb黄色大片| 亚洲精品亚洲一区二区| 欧美黄色淫秽网站| 亚洲美女黄片视频| 欧美极品一区二区三区四区| 亚洲国产精品sss在线观看| 久久久久精品国产欧美久久久| 黄片小视频在线播放| 午夜精品一区二区三区免费看| 3wmmmm亚洲av在线观看| 亚洲美女搞黄在线观看 | 国产精品一区二区免费欧美| 久久精品综合一区二区三区| 成人高潮视频无遮挡免费网站| 91久久精品国产一区二区成人| 成人国产综合亚洲| 国产黄a三级三级三级人| 99久国产av精品| 色吧在线观看| 亚洲久久久久久中文字幕| av中文乱码字幕在线| 中出人妻视频一区二区| 欧美bdsm另类| 精品午夜福利视频在线观看一区| 欧美午夜高清在线| 美女高潮的动态| 精品国产三级普通话版| 国产在线精品亚洲第一网站| 国产欧美日韩精品亚洲av| 欧美高清成人免费视频www| 亚洲欧美日韩卡通动漫| 亚洲专区国产一区二区| 国产淫片久久久久久久久 | 麻豆国产97在线/欧美| 午夜免费男女啪啪视频观看 | av在线蜜桃| 欧美成人一区二区免费高清观看| 三级国产精品欧美在线观看| 97超级碰碰碰精品色视频在线观看| 久久精品国产自在天天线| 亚洲天堂国产精品一区在线| 男女床上黄色一级片免费看| 日日摸夜夜添夜夜添小说| 观看美女的网站| 在线观看舔阴道视频| 一边摸一边抽搐一进一小说| 国产亚洲精品久久久久久毛片| 久久久色成人| 此物有八面人人有两片| 中文字幕免费在线视频6| 色综合亚洲欧美另类图片| 特大巨黑吊av在线直播| 久久久久久久久大av| 老女人水多毛片| 亚洲熟妇中文字幕五十中出| 淫妇啪啪啪对白视频| 久久久精品欧美日韩精品| 91在线观看av| 亚洲国产精品成人综合色| 欧美区成人在线视频| 好男人在线观看高清免费视频| 日本五十路高清| 欧美另类亚洲清纯唯美| 精品午夜福利在线看| 欧美日本视频| 欧美成人一区二区免费高清观看| 男人舔奶头视频| 直男gayav资源| 床上黄色一级片| 成人永久免费在线观看视频| 嫩草影院入口| 久久热精品热| 成人亚洲精品av一区二区| 免费电影在线观看免费观看| 国产在线男女|