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

    基于XRF和Vis-NIR光譜數(shù)據(jù)融合的土壤鎘含量定量分析法

    2020-11-09 10:54:02王清亞李福生江曉宇鄔書良謝濤鋒黃溫鋼
    分析測(cè)試學(xué)報(bào) 2020年11期
    關(guān)鍵詞:融合模型

    王清亞,李福生*,江曉宇,鄔書良,謝濤鋒,黃溫鋼

    (1.東華理工大學(xué) 核技術(shù)應(yīng)用教育部工程研究中心,江西 南昌 330013;2.東華理工大學(xué) 核資源與環(huán)境國(guó)家重點(diǎn)實(shí)驗(yàn)室,江西 南昌 330013)

    隨著工業(yè)進(jìn)程的加快,鎘(Cd)的大量利用造成了嚴(yán)重的環(huán)境污染,全世界每年約向環(huán)境中釋放40萬噸的Cd,其中59%直接進(jìn)入土壤[1]。土壤中Cd含量的傳統(tǒng)檢測(cè)方法存在效率低、成本高、易造成環(huán)境污染等問題,無法滿足較大尺度區(qū)域的土壤鎘含量的調(diào)查需求[2]。近二十年來,近端土壤傳感技術(shù),如X射線熒光光譜(X-ray fluorescence spectroscopy,XRF)[3]和可見光-近紅外(Visible and near-infrared spectra,Vis-NIR)光譜分析技術(shù)[4],因可以實(shí)現(xiàn)土壤中重金屬含量的無損、高效和無污染快速測(cè)定[5],引起了研究人員的關(guān)注。

    但不同類型的土壤基質(zhì)差異大,元素種類繁多,基質(zhì)效應(yīng)明顯[6],單一采用XRF和Vis-NIR時(shí)準(zhǔn)確度[3,7]難以滿足國(guó)標(biāo)對(duì)Cd含量測(cè)定的要求[8]。多源數(shù)據(jù)融合可以將多個(gè)信息源或傳感器獲得的同一個(gè)目標(biāo)的不同信息融合到一起,利用計(jì)算機(jī)技術(shù)對(duì)這些信息進(jìn)行分析、支配和使用,做出決策,往往較單一信息更具準(zhǔn)確性[9]。目前多源數(shù)據(jù)融合已經(jīng)廣泛應(yīng)用于光譜分析領(lǐng)域。武中臣等[10]采用兩種數(shù)據(jù)融合(等權(quán)融合和累加融合)對(duì)硅酸鹽、硫酸鹽、碳酸鹽等4種礦物進(jìn)行分類研究,發(fā)現(xiàn)融合后的預(yù)測(cè)結(jié)果均好于單一檢測(cè)方法。O’rourke等[11]將可見-近紅外與中紅外光譜等權(quán)融合后,測(cè)試了愛爾蘭國(guó)家土壤數(shù)據(jù)庫中的樣品,結(jié)果顯示,相比單一檢測(cè)技術(shù),多傳感器的等權(quán)融合技術(shù)提高了檢測(cè)土壤多項(xiàng)指標(biāo)的準(zhǔn)確性。Stenberg等[12]將Vis-NIR和XRF光譜等權(quán)融合,分別獲得了土壤有機(jī)質(zhì)(SOM)礦物、水分、pH值、銅、鉛含量等重要土壤屬性。史舟等[5]通過數(shù)據(jù)外積融合測(cè)定了土壤中鉻的含量,相比單一數(shù)據(jù)預(yù)測(cè)模型,融合后的模型預(yù)測(cè)精度(ρc=0.88,RMSEP=6.8,RPD=2.30)更佳。綜合以上技術(shù)在土壤相關(guān)參數(shù)定量檢測(cè)方面的研究可知,基于數(shù)據(jù)融合測(cè)定土壤Cd含量是可行的,但目前相關(guān)研究較少。

    本研究以鄱陽湖南磯山為研究區(qū),進(jìn)行野外土壤采樣、室內(nèi)理化分析、光譜采集與處理等工作,并基于XRF和Vis-NIR研究了不同數(shù)據(jù)融合算法對(duì)土壤Cd含量預(yù)測(cè)模型穩(wěn)定性和精度的影響,探索基于兩種光譜數(shù)據(jù)融合的Cd含量的預(yù)測(cè)方法,以為該區(qū)域的重金屬污染調(diào)查提供參考。

    1 實(shí)驗(yàn)部分

    1.1 研究區(qū)域與樣品取樣

    研究的區(qū)域位于江西鄱陽湖西南岸(116°10′24″~116°23′50″E,28°52′21″~29°06′46″N),共采集371個(gè)樣品,土壤類型覆蓋潮土、水稻土、黃棕壤,采用十字法取樣,采樣深度為0~20 cm,采樣時(shí)剔除侵入體,混勻后用四分法收集土樣。經(jīng)風(fēng)干、研磨,過2 mm孔徑篩后,將每份樣品分成3份,分別用于室內(nèi)理化分析和XRF、Vis-NIR光譜采集。

    1.2 實(shí)驗(yàn)室內(nèi)化學(xué)分析

    按照《石墨爐原子吸收分光光度法GB/T 17141-1997》進(jìn)行樣品消解和分析,用HCl-HNO3-HF-HClO4(四酸)破壞礦物晶格,以石墨爐原子吸收分光光度計(jì)測(cè)定試劑中的Cd[13]。約20%的土壤被分析兩次,用以評(píng)估方法的重復(fù)性。設(shè)置95%的置信區(qū)間,剔除超出置信區(qū)間的樣品,最終保留356個(gè)樣品進(jìn)行后續(xù)研究。選用K-S(Kennard-Stone)算法[14]計(jì)算各個(gè)樣品的歐式距離,選用60%的樣本用于建模,40%的樣品用于預(yù)測(cè)和驗(yàn)證(表1)。

    表1 土壤樣品鎘含量統(tǒng)計(jì)特征Table 1 Statistical characteristics of soil Cd content

    1.3 Vis-NIR光譜與XRF光譜測(cè)量

    采用美國(guó)ASD公司的便攜式FieldSpec3高光譜分辨率地物波譜儀測(cè)定土壤樣品的Vis-NIR光譜,其波段值為350~2 500 nm,采樣間隔為1.4 nm(350~1 000 nm)和2 nm(1 000~2 500 nm),重采樣間隔為1 nm。將樣品放置在培養(yǎng)基中(直徑10 cm,高度1 cm),隨機(jī)掃描10次,取信噪比最佳的3次測(cè)定結(jié)果進(jìn)行平均,作為樣品光譜。

    XRF測(cè)量采用泰克松德公司的TS-XH4000-G分析儀,選用土壤模式,每次測(cè)試90 s,同一個(gè)樣品測(cè)試3次。每測(cè)30個(gè)樣品校正1次,最終取3次測(cè)試結(jié)果的平均作為該樣品的光譜。TS-XH4000-G分析儀對(duì)Cd的檢出限為2 mg·kg-1。測(cè)試結(jié)束后,對(duì)測(cè)試結(jié)果進(jìn)行分析,在371個(gè)樣品中只有14個(gè)樣品的Cd含量高于儀器檢出限。大部分樣品的Cd含量低于檢出限。文中儀器對(duì)Cd含量的分析結(jié)果不作為建立模型的考慮范圍。

    1.4 數(shù)據(jù)分析及評(píng)價(jià)指標(biāo)

    最小二乘支持向量機(jī)(Least squares support vector machine,LS-SVM)是基于支持向量機(jī)(SVM)的一種改進(jìn)算法,通過構(gòu)造損失函數(shù)將原SVM中算法的二次尋優(yōu)變?yōu)榍蠼饩€性方程,其求解速度較快,在光譜分析領(lǐng)域中得到了廣泛應(yīng)用[15]。LS-SVM方法共有4種核函數(shù)(線性核、多項(xiàng)式核、徑向基函數(shù)核、高斯核)可供選擇,徑向基函數(shù)核(Radial basis function,RBF)非常適用于光譜分析中的非線性問題處理[16],因此本文使用RBF函數(shù)來實(shí)現(xiàn)LS-SVM建模。

    模型預(yù)測(cè)精度以決定系數(shù)(Determination coefficients,R2)、校正均方根誤差(Root mean squared error of calibration,RMSEC)、預(yù)測(cè)均方根誤差(Root mean square error of prediction,RMSEP)、相對(duì)分析誤差(Relative percent deviation,RPD)4個(gè)參數(shù)衡量。R2反映模型建立和驗(yàn)證的穩(wěn)定性,RMSEC和RMSEP用來檢驗(yàn)?zāi)P偷念A(yù)報(bào)能力[17],RPD是樣本標(biāo)準(zhǔn)差與RMSEP的比值,用以判斷模型的預(yù)測(cè)能力,RPD<1.4時(shí),模型無法對(duì)樣品進(jìn)行預(yù)測(cè);1.4≤RPD<2.0時(shí),認(rèn)為模型效果一般,可用來對(duì)樣品進(jìn)行粗略評(píng)估;RPD≥2.0時(shí),模型具有極好的預(yù)測(cè)能力[18]。

    圖1 不同Cd含量土壤的原始XRF光譜(A)和Vis-NIR光譜(B)Fig.1 Theoriginal XRF spectra (A) and Vis-NIR spectra(B) of soil with different Cd contents

    2 結(jié)果與討論

    2.1 土壤樣品的XRF光譜與Vis-NIR光譜特征

    選取土壤樣品中Cd含量最大值、中間值、最小值(Cd含量分別為1.50、0.76、0.21 mg·kg-1)所對(duì)應(yīng)的3條原始光譜曲線進(jìn)行分析(圖1)。

    根據(jù)莫塞萊定律,土壤中Cd產(chǎn)生的K系和L系特征射線的能量是唯一的。圖1A中23.17 keV處為Cd的Kα特征射線,其能量計(jì)數(shù)與Cd含量呈正相關(guān)。但Cd的L系特征射線(3.13 keV)處于傳感器探測(cè)范圍的邊緣,靈敏度較差,不可作為定量模型的輸入量,故許多定量?jī)x器只將Cd的Kα特征射線計(jì)數(shù)作為輸入量。而復(fù)雜的原級(jí)X射線與土壤樣品的作用過程和元素間特征射線的相互影響(例如Pb的23.82 keV Kβ和Cd的23.17 keV Kα射線相互影響)導(dǎo)致目標(biāo)元素的照射量率與含量之間不再是單純的正比關(guān)系,這也可從圖1A中的特征射線計(jì)數(shù)看出來。而廣泛用于土壤重金屬(如Pb、Cu、Zn、Cr等)檢測(cè)的XRF儀器利用元素特征射線與元素含量之間的正相關(guān)關(guān)系進(jìn)行定量分析,從而導(dǎo)致對(duì)Cd類亞mg·kg-1級(jí)含量的重金屬的檢測(cè)效果不佳。

    不同Cd含量土壤樣品的Vis-NIR光譜(圖1B)具有以下特征:(1)光譜形態(tài)大致相同,近似平行,均呈上凸趨勢(shì),總體單調(diào)遞增;(2)在可見光波段的反射率小于近紅外波段,不同光譜間的差別也較??;(3)3條光譜曲線的差異主要在近紅外波段,出現(xiàn)特征吸收峰的位置大致相同,吸收深度和吸收面積有差別;(4)在400~600 nm區(qū)間斜率較大,和土壤中鐵的含量有關(guān);(5)在1 400、1 900、2 200 nm處存在明顯的水分吸收谷,1 400 nm附近為羥基(OH—)譜帶,1 900 nm附近為H2O譜帶,2 200 nm附近為羥基伸縮振動(dòng)與Al—OH和Mg—OH彎曲振動(dòng)的合譜帶[19]。

    2.2 光譜處理與變換

    為消除噪音和低能射線的影響,首先在光譜中截取數(shù)據(jù)較為豐富的一段:Vis-NIR反射光譜截取450~2 450 nm段(2 000個(gè)通道),并將反射光譜轉(zhuǎn)換成吸收光譜;XRF光譜選取0.405~42.105 keV段(2 000通道)。隨后,對(duì)兩種光譜進(jìn)行基線校正、Sacitzky-Golay(S-G)平滑[20](窗口大小為15,2階多項(xiàng)式)以及標(biāo)準(zhǔn)矢量歸一化(SNV)[21]處理。最后,對(duì)處理后的Vis-NIR和XRF光譜重新采樣(保留200個(gè)通道,分辨率分別為10 nm,0.20 keV),Vis-NIR光譜的數(shù)值范圍為-3.42~4.29,XRF光譜的數(shù)值范圍為-1.16~2.37,解決了數(shù)據(jù)融合中量綱不一致的問題。

    2.3 數(shù)據(jù)融合

    Moros等[22]提出了3種數(shù)據(jù)融合算法:累加融合(Coadditionfusion,CF)、等權(quán)融合(Equal rights fusion,ERF)、外積融合(Outer product fusion,OPF)。累加融合是將XRF和Vis-NIR光譜相關(guān)通道上的數(shù)據(jù)按照一定權(quán)重進(jìn)行加和,得到的新矩陣的行數(shù)等于樣品數(shù)量,列數(shù)不變。本研究中將對(duì)XRF和Vis-NIR光譜進(jìn)行加權(quán)(XRF∶Vis-NIR=3∶1)累加融合得到的定量模型記為CF-LSSVM;等權(quán)融合是將XRF和Vis-NIR光譜數(shù)據(jù)首尾相連,得到的新矩陣行數(shù)為樣品數(shù)量,列數(shù)為兩個(gè)光譜數(shù)據(jù)的列數(shù)之和,定量模型記為ERF-LSSVM;外積融合則通過求取XRF和Vis-NIR光譜數(shù)據(jù)的克羅內(nèi)克積(Kronecker product)[23]進(jìn)行融合:如果A是m×n矩陣,B是p×q矩陣,則求取的克羅內(nèi)克積為mp×nq的分塊矩陣[24],見公式(2)。應(yīng)用克羅內(nèi)克積的形式可以有效增強(qiáng)不同物理量光譜的共同演化性能[23,25-27]。

    (2)

    如圖2所示,外積融合的具體方法是:第m個(gè)被測(cè)樣品可以得到一組具有r個(gè)變量的XRF光譜(1×r)和c個(gè)變量的Vis-NIR光譜(1×c),分別記為rm、cm(公式3),將2個(gè)矩陣進(jìn)行相乘求取克羅克內(nèi)積,得到c×r矩陣(公式4),該矩陣為第m個(gè)樣品的外積融合光譜。將得到的矩陣按照1行的形式展開(Unfold)得到(1×r·c)矩陣,以此矩陣作為建立模型的輸入量,定量模型記為OPF-LSSVM。

    圖2 外積融合算法示意圖Fig.2 Process of data fusion by OPF

    rm=(r1,m,…,ri,m)
    cm=(c1,m,…,cj,m)

    (3)

    (4)

    (5)

    其中rm、cm分別為第m個(gè)樣品的XRF光譜和Vis-NIR光譜,Am是該土壤樣品的外積融合光譜,ri,m是XRF光譜在i能量處的通道計(jì)數(shù),cj,m代表同一個(gè)樣品在j波長(zhǎng)處的吸光度。Ψ為Am矩陣的展開,為1行c×r列矩陣。

    以1號(hào)土壤樣品為例,經(jīng)XRF光譜、Vis-NIR光譜融合算法后,構(gòu)建了CF光譜、ERF光譜、OPF光譜,如圖3所示。從圖中可以看出,CF光譜走勢(shì)和XRF光譜類似,保留了XRF光譜和Vis-NIR中的一些特征峰;ERF光譜圖比CF光譜圖含有更多的信息,包括XRF和Vis-NIR光譜中的所有特征峰;OPF光譜圖所含數(shù)據(jù)量最多,達(dá)到40 000個(gè)變量,其數(shù)值的變化也較大。

    圖3 1號(hào)土壤樣品的XRF光譜圖(A)、Vis-NIR光譜圖(B)、累加融合光譜圖(C)、等權(quán)融合光譜圖(D)和外積融合光譜圖(E)Fig.3 XRF spectra(A),Vis-NIR spectra(B) and their coaddition fusion spectra(C),equal rights fusion spectra (D) andouter product fusion spectra(E) of No.1 soil sample

    2.4 LS-SVM模型

    在建立LS-SVM模型時(shí),首先要選擇合適的核參數(shù)σ2和懲罰系數(shù)γ,通過網(wǎng)格搜索法和留一法(Leave-one-out)進(jìn)行全局尋優(yōu)。數(shù)據(jù)處理中γ的搜索范圍為1~10 000,σ2的搜索范圍為0.01~1 000,確定的最優(yōu)γ和σ2見表2。

    表2 LS-SVM模型中γ和σ2最優(yōu)值Table 2 The optimal values of γ and σ2 in LS-SVM model

    2.4.1 單一光譜預(yù)測(cè)模型使用XRF和Vis-NIR光譜作為輸入變量,建立LS-SVM模型,具體預(yù)測(cè)效果如圖4所示。從圖中可以看出,以單一光譜作為輸入量的模型中,XRF-LS-SVM模型的R2(0.63)略低于Vis-NIR-LS-SVM模型(R2=0.69),其RMSEP(0.13)則高于Vis-NIR-LS-SVM(RMSEP=0.10),表明Vis-NIR-LS-SVM模型具有很好的穩(wěn)定性和精度,略優(yōu)于XRF-LS-SVM的預(yù)測(cè)效果,但兩種模型的預(yù)測(cè)效果均不理想,RPD介于1.40~2.0之間。

    XRF儀器[28]對(duì)Cd類亞mg·kg-1級(jí)含量的重金屬的檢測(cè)效果不佳,使用XRF全光譜作為模型的輸入,雖可以提高痕量元素的檢出限[29],但從實(shí)驗(yàn)結(jié)果看,預(yù)測(cè)精度和穩(wěn)定性仍然無法滿足實(shí)際檢測(cè)要求。

    土壤中的重金屬常與土壤有機(jī)質(zhì)、粘土礦物結(jié)合,這些物質(zhì)在Vis-NIR光譜上具有強(qiáng)烈的特征吸收[30],可基于此使用Vis-NIR對(duì)土壤中的重金屬進(jìn)行預(yù)測(cè)。使用Vis-NIR光譜預(yù)測(cè)土壤Cd含量時(shí),效果略優(yōu)于XRF光譜預(yù)測(cè)模型,這與文獻(xiàn)研究結(jié)果類似[31]。

    表3 不同方法預(yù)測(cè)精度的比較Table 3 Comparison of prediction accuracy of different methods

    2.4.2 數(shù)據(jù)融合預(yù)測(cè)模型在光譜融合的預(yù)測(cè)模型中,CF-LS-SVM的RPD=1.24,基本沒有預(yù)測(cè)能力,從1號(hào)土壤樣品CF光譜圖中可以看出(圖3C),融合后光譜集合了XRF和Vis-NIR光譜的一些信息,但兩個(gè)光譜值相加,導(dǎo)致一些光譜信息消失,增強(qiáng)了融合后光譜的非線性。除CF-LS-SVM外,融合后的模型預(yù)測(cè)效果均優(yōu)于單一光譜模型,表明數(shù)據(jù)融合可以為Cd含量預(yù)測(cè)提供有效方法(表3)。ERF光譜保留了XRF和Vis-NIR的所有光譜信息。ERF-LS-SVM模型的RMSEP=0.09,比兩個(gè)單一光譜及CF-LS-SVM模型效果更佳,但其RPD=1.92也表明ERF-LS-SVM預(yù)測(cè)能力一般,僅可對(duì)土壤中Cd含量進(jìn)行粗略預(yù)測(cè)。OPF-LSSVM在所有模型中效果最好,相比ERF-LS-SVM模型,其RMSEP從0.09降低到0.06,預(yù)測(cè)精度更高,RPD=2.41則表明該模型具有極好的預(yù)測(cè)能力。

    上述結(jié)果說明數(shù)據(jù)融合相較于單一光譜包含了更多的土壤信息,其中外積融合通過求取兩個(gè)光譜的克羅克內(nèi)積,增強(qiáng)了光譜的區(qū)域性,尤其是將不同物理量的XRF和Vis-NIR光譜外積融合后,可以充分利用光譜的不同性質(zhì)和互補(bǔ)信息,更利于土壤參數(shù)的預(yù)測(cè)建模。但在具體運(yùn)算過程中,外積融合的光譜變量是單一光譜變量的幾百倍,建模計(jì)算所消耗的時(shí)間也相對(duì)增加。

    3 結(jié) 論

    本文基于LS-SVM算法對(duì)比了XRF和Vis-NIR單一光譜模型及其等權(quán)融合模型、累加融合模型、外積融合模型的預(yù)測(cè)能力,結(jié)果顯示外積融合模型的預(yù)測(cè)精度和模型穩(wěn)定性最佳(R2=0.85,RMSEP=0.06,RPD=2.41),符合土壤調(diào)查中對(duì)鎘含量的預(yù)測(cè)要求。本研究為土壤重金屬含量的調(diào)查,尤其是Cd、Hg類亞mg·kg-1級(jí)重金屬含量的調(diào)查提供了借鑒方法,對(duì)開發(fā)土壤重金屬分析儀具有重要意義。

    猜你喜歡
    融合模型
    一半模型
    一次函數(shù)“四融合”
    村企黨建聯(lián)建融合共贏
    融合菜
    從創(chuàng)新出發(fā),與高考數(shù)列相遇、融合
    重要模型『一線三等角』
    寬窄融合便攜箱IPFS500
    《融合》
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    3D打印中的模型分割與打包
    日本欧美国产在线视频| 赤兔流量卡办理| 免费观看在线日韩| 亚洲一级一片aⅴ在线观看| 久久久久久久久免费视频了| 一区二区三区激情视频| 亚洲欧美成人精品一区二区| 国产精品偷伦视频观看了| 七月丁香在线播放| av一本久久久久| 精品第一国产精品| 午夜日韩欧美国产| www日本在线高清视频| 中文乱码字字幕精品一区二区三区| 丰满饥渴人妻一区二区三| www.自偷自拍.com| 男女无遮挡免费网站观看| 成年女人在线观看亚洲视频| 亚洲第一青青草原| 午夜福利在线免费观看网站| 欧美激情高清一区二区三区 | 老鸭窝网址在线观看| 国产精品无大码| 日本vs欧美在线观看视频| 欧美另类一区| 91精品伊人久久大香线蕉| 美女国产高潮福利片在线看| 国产一级毛片在线| 伊人久久国产一区二区| 国产日韩欧美在线精品| 亚洲精品一二三| 欧美最新免费一区二区三区| av网站免费在线观看视频| 一区在线观看完整版| 可以免费在线观看a视频的电影网站 | 大陆偷拍与自拍| 交换朋友夫妻互换小说| 欧美日韩视频精品一区| 日本vs欧美在线观看视频| 激情五月婷婷亚洲| av卡一久久| 亚洲欧美一区二区三区国产| 亚洲美女搞黄在线观看| 成年女人在线观看亚洲视频| 女人久久www免费人成看片| 欧美激情极品国产一区二区三区| 香蕉丝袜av| 2022亚洲国产成人精品| 观看av在线不卡| 亚洲综合色网址| 欧美日韩视频高清一区二区三区二| 欧美人与性动交α欧美软件| 精品国产超薄肉色丝袜足j| 狠狠精品人妻久久久久久综合| 久久久久久人妻| 男女无遮挡免费网站观看| 欧美xxⅹ黑人| 欧美少妇被猛烈插入视频| 五月天丁香电影| 久久狼人影院| 国产精品国产三级专区第一集| 青春草亚洲视频在线观看| 麻豆精品久久久久久蜜桃| 黄片小视频在线播放| 欧美亚洲 丝袜 人妻 在线| 搡老乐熟女国产| 亚洲精品日本国产第一区| 男人舔女人的私密视频| 肉色欧美久久久久久久蜜桃| 男女高潮啪啪啪动态图| 丝袜在线中文字幕| 最近中文字幕2019免费版| 色吧在线观看| 狂野欧美激情性bbbbbb| 国产在线视频一区二区| 十分钟在线观看高清视频www| 婷婷色综合www| 日本欧美视频一区| 欧美日韩成人在线一区二区| 久久人人爽av亚洲精品天堂| 国产精品久久久久久久久免| 考比视频在线观看| 成人亚洲欧美一区二区av| 一级,二级,三级黄色视频| 亚洲国产av新网站| 多毛熟女@视频| 卡戴珊不雅视频在线播放| a 毛片基地| 久久精品国产a三级三级三级| 欧美黄色片欧美黄色片| 中文字幕另类日韩欧美亚洲嫩草| 街头女战士在线观看网站| 国产男人的电影天堂91| 最近中文字幕高清免费大全6| 男人爽女人下面视频在线观看| 如何舔出高潮| 午夜免费男女啪啪视频观看| 免费观看在线日韩| 少妇人妻 视频| 黄片播放在线免费| 在现免费观看毛片| 女人被躁到高潮嗷嗷叫费观| 80岁老熟妇乱子伦牲交| 精品国产乱码久久久久久小说| 午夜福利视频精品| 秋霞伦理黄片| av福利片在线| 2021少妇久久久久久久久久久| 日本欧美视频一区| 欧美精品一区二区大全| 看非洲黑人一级黄片| 精品国产一区二区久久| 久久精品久久久久久噜噜老黄| 天天躁夜夜躁狠狠久久av| 亚洲精品自拍成人| 日本午夜av视频| www.精华液| 久久韩国三级中文字幕| 免费黄频网站在线观看国产| 欧美亚洲 丝袜 人妻 在线| 久久99热这里只频精品6学生| 在线亚洲精品国产二区图片欧美| 午夜福利网站1000一区二区三区| 青青草视频在线视频观看| 成年av动漫网址| 国产有黄有色有爽视频| 有码 亚洲区| 韩国高清视频一区二区三区| 久久久久久久久久久免费av| 一本—道久久a久久精品蜜桃钙片| 亚洲欧美清纯卡通| 人妻一区二区av| 看非洲黑人一级黄片| 日韩欧美精品免费久久| a级片在线免费高清观看视频| 69精品国产乱码久久久| 亚洲精品日本国产第一区| 日韩在线高清观看一区二区三区| 日韩一卡2卡3卡4卡2021年| av不卡在线播放| 午夜免费观看性视频| 成人手机av| 啦啦啦视频在线资源免费观看| 久久久a久久爽久久v久久| 母亲3免费完整高清在线观看 | 熟妇人妻不卡中文字幕| 久久99一区二区三区| 国产高清国产精品国产三级| 女人高潮潮喷娇喘18禁视频| av在线观看视频网站免费| 一二三四中文在线观看免费高清| 午夜福利乱码中文字幕| 日韩不卡一区二区三区视频在线| 久久久久国产精品人妻一区二区| 亚洲成人av在线免费| 欧美精品国产亚洲| 久久久久久免费高清国产稀缺| 女性生殖器流出的白浆| av一本久久久久| 国产欧美日韩一区二区三区在线| 欧美xxⅹ黑人| av在线播放精品| 男女午夜视频在线观看| 黄片无遮挡物在线观看| 日韩欧美精品免费久久| 两个人免费观看高清视频| 欧美精品一区二区免费开放| 欧美国产精品va在线观看不卡| 少妇被粗大的猛进出69影院| 日本vs欧美在线观看视频| 777久久人妻少妇嫩草av网站| 欧美人与善性xxx| 日韩一卡2卡3卡4卡2021年| 97人妻天天添夜夜摸| 国产一区二区三区综合在线观看| 熟女少妇亚洲综合色aaa.| 一区二区三区激情视频| 免费女性裸体啪啪无遮挡网站| 搡老乐熟女国产| 久久精品国产亚洲av天美| 午夜福利网站1000一区二区三区| 精品少妇内射三级| 欧美激情 高清一区二区三区| 久久久久精品久久久久真实原创| 一级黄片播放器| 久久精品人人爽人人爽视色| 精品第一国产精品| 免费黄频网站在线观看国产| 欧美日韩精品成人综合77777| 超碰成人久久| 尾随美女入室| 在线观看三级黄色| 伊人久久大香线蕉亚洲五| 欧美精品人与动牲交sv欧美| 精品久久久精品久久久| 人妻一区二区av| 中国三级夫妇交换| 狠狠精品人妻久久久久久综合| 在线看a的网站| 欧美日韩视频高清一区二区三区二| 国产野战对白在线观看| 日本午夜av视频| 日本欧美国产在线视频| 国产精品一区二区在线观看99| 男女边吃奶边做爰视频| 国产精品麻豆人妻色哟哟久久| 青草久久国产| 人妻少妇偷人精品九色| 波野结衣二区三区在线| 久久综合国产亚洲精品| 亚洲成人av在线免费| 国产精品人妻久久久影院| 亚洲精品中文字幕在线视频| 美女大奶头黄色视频| 久久国产亚洲av麻豆专区| 男女无遮挡免费网站观看| 26uuu在线亚洲综合色| 国产 精品1| 女性生殖器流出的白浆| 人人妻人人添人人爽欧美一区卜| 国产av一区二区精品久久| 日本wwww免费看| 久久精品夜色国产| 满18在线观看网站| 免费观看无遮挡的男女| 伦精品一区二区三区| 久久久久国产一级毛片高清牌| 女人被躁到高潮嗷嗷叫费观| 在线观看免费日韩欧美大片| 人妻少妇偷人精品九色| 狠狠婷婷综合久久久久久88av| 综合色丁香网| 男女无遮挡免费网站观看| 欧美日韩一级在线毛片| 国产精品国产三级国产专区5o| 国产有黄有色有爽视频| 少妇猛男粗大的猛烈进出视频| 人妻 亚洲 视频| 日本wwww免费看| 建设人人有责人人尽责人人享有的| 在线天堂最新版资源| 天美传媒精品一区二区| 国产不卡av网站在线观看| 18禁裸乳无遮挡动漫免费视频| 黄色 视频免费看| 一级片免费观看大全| 大片免费播放器 马上看| 18+在线观看网站| 十分钟在线观看高清视频www| 999精品在线视频| 男女免费视频国产| 丰满少妇做爰视频| 欧美激情极品国产一区二区三区| 色视频在线一区二区三区| 一区福利在线观看| 一级毛片 在线播放| 美女视频免费永久观看网站| 欧美日韩成人在线一区二区| 久久精品夜色国产| 伊人亚洲综合成人网| 亚洲国产欧美网| 深夜精品福利| 亚洲av电影在线观看一区二区三区| 一个人免费看片子| 女人精品久久久久毛片| 亚洲伊人色综图| 成人手机av| 欧美日韩av久久| 美女脱内裤让男人舔精品视频| 欧美日韩视频高清一区二区三区二| 成年女人在线观看亚洲视频| 大码成人一级视频| 波多野结衣av一区二区av| 99香蕉大伊视频| 亚洲五月色婷婷综合| av.在线天堂| 久久国产亚洲av麻豆专区| 欧美激情 高清一区二区三区| 大片电影免费在线观看免费| 亚洲国产毛片av蜜桃av| 99热网站在线观看| 久热久热在线精品观看| 免费人妻精品一区二区三区视频| 免费黄色在线免费观看| 麻豆av在线久日| 在线 av 中文字幕| 黑丝袜美女国产一区| 黄色怎么调成土黄色| 少妇的逼水好多| 亚洲精品美女久久av网站| 桃花免费在线播放| 亚洲精品一区蜜桃| 亚洲国产av影院在线观看| 老汉色∧v一级毛片| 亚洲精品日韩在线中文字幕| 国产 精品1| 亚洲第一区二区三区不卡| 亚洲成人手机| 波野结衣二区三区在线| 大片免费播放器 马上看| 日本黄色日本黄色录像| 观看美女的网站| 免费看不卡的av| 亚洲精品国产一区二区精华液| 老女人水多毛片| 99热国产这里只有精品6| 亚洲成人一二三区av| 久久久亚洲精品成人影院| 在现免费观看毛片| 久久99一区二区三区| av电影中文网址| 久久女婷五月综合色啪小说| 国产成人a∨麻豆精品| 不卡视频在线观看欧美| 精品少妇久久久久久888优播| 18禁国产床啪视频网站| 国产日韩一区二区三区精品不卡| 宅男免费午夜| 青春草亚洲视频在线观看| 人妻 亚洲 视频| av网站在线播放免费| 一个人免费看片子| 精品久久蜜臀av无| 丝袜喷水一区| 色吧在线观看| 欧美精品高潮呻吟av久久| 亚洲欧美成人综合另类久久久| 老汉色∧v一级毛片| 在线观看免费日韩欧美大片| 热re99久久精品国产66热6| 日韩av不卡免费在线播放| av女优亚洲男人天堂| 久久人妻熟女aⅴ| 一本久久精品| 精品国产乱码久久久久久男人| 亚洲婷婷狠狠爱综合网| 国产精品.久久久| 亚洲情色 制服丝袜| 日韩在线高清观看一区二区三区| 免费av中文字幕在线| 秋霞伦理黄片| 一级片'在线观看视频| 丁香六月天网| 国产成人免费观看mmmm| 波多野结衣av一区二区av| av片东京热男人的天堂| 亚洲美女搞黄在线观看| av片东京热男人的天堂| 亚洲国产毛片av蜜桃av| av网站在线播放免费| 久久久a久久爽久久v久久| 七月丁香在线播放| 国产日韩欧美在线精品| 国产精品免费大片| 赤兔流量卡办理| 欧美少妇被猛烈插入视频| 久久这里有精品视频免费| 欧美中文综合在线视频| 80岁老熟妇乱子伦牲交| 青春草亚洲视频在线观看| 成人午夜精彩视频在线观看| 亚洲欧洲精品一区二区精品久久久 | 另类精品久久| av福利片在线| 久久精品国产鲁丝片午夜精品| 视频区图区小说| 成年av动漫网址| 久久久久精品久久久久真实原创| 亚洲精品一区蜜桃| 日韩三级伦理在线观看| 秋霞伦理黄片| 国产日韩一区二区三区精品不卡| 十八禁网站网址无遮挡| 一二三四中文在线观看免费高清| 久久久久久久久久人人人人人人| 欧美人与性动交α欧美软件| av网站免费在线观看视频| 国产精品99久久99久久久不卡 | 久久久国产一区二区| 美女中出高潮动态图| 啦啦啦在线观看免费高清www| 亚洲欧美色中文字幕在线| 自拍欧美九色日韩亚洲蝌蚪91| 国产日韩欧美在线精品| 欧美另类一区| 亚洲成av片中文字幕在线观看 | 亚洲,欧美精品.| 在线观看国产h片| 久久av网站| 久久综合国产亚洲精品| 久久久久精品人妻al黑| 一本—道久久a久久精品蜜桃钙片| 妹子高潮喷水视频| 久久精品国产自在天天线| 狠狠精品人妻久久久久久综合| 亚洲欧洲国产日韩| 丝袜喷水一区| 免费av中文字幕在线| 亚洲精华国产精华液的使用体验| 在线观看美女被高潮喷水网站| 国产一区二区三区综合在线观看| 王馨瑶露胸无遮挡在线观看| av.在线天堂| 国产精品嫩草影院av在线观看| 在线天堂最新版资源| 亚洲欧美清纯卡通| 久久精品人人爽人人爽视色| 午夜免费鲁丝| 午夜免费男女啪啪视频观看| 国产成人免费观看mmmm| 水蜜桃什么品种好| 精品国产乱码久久久久久男人| 成人影院久久| 久久精品aⅴ一区二区三区四区 | 欧美97在线视频| 欧美国产精品一级二级三级| 中文字幕人妻丝袜一区二区 | 欧美成人午夜免费资源| 爱豆传媒免费全集在线观看| av视频免费观看在线观看| 欧美精品av麻豆av| 欧美激情极品国产一区二区三区| 嫩草影院入口| 亚洲一码二码三码区别大吗| 欧美成人午夜精品| 欧美bdsm另类| 欧美av亚洲av综合av国产av | 免费日韩欧美在线观看| 18在线观看网站| 免费看av在线观看网站| 精品一区二区免费观看| 久久精品熟女亚洲av麻豆精品| 色哟哟·www| 宅男免费午夜| 美女福利国产在线| 十八禁网站网址无遮挡| 国产视频首页在线观看| 秋霞伦理黄片| 一级片'在线观看视频| av在线老鸭窝| av片东京热男人的天堂| 亚洲图色成人| 成人国产麻豆网| 黄色 视频免费看| 黄片无遮挡物在线观看| 大片免费播放器 马上看| 久久久国产精品麻豆| 天天躁日日躁夜夜躁夜夜| 成人毛片60女人毛片免费| 一级a爱视频在线免费观看| 一区二区三区精品91| 国产欧美日韩一区二区三区在线| 国产97色在线日韩免费| 另类亚洲欧美激情| 欧美激情极品国产一区二区三区| 热re99久久国产66热| 精品国产乱码久久久久久男人| 美女国产高潮福利片在线看| 久久99一区二区三区| 国产精品熟女久久久久浪| 成人亚洲欧美一区二区av| 午夜福利视频精品| 久久久久国产网址| 国产成人av激情在线播放| 亚洲国产精品999| 91久久精品国产一区二区三区| 日韩电影二区| 免费不卡的大黄色大毛片视频在线观看| 伊人久久大香线蕉亚洲五| 18在线观看网站| 老司机亚洲免费影院| 香蕉丝袜av| 一区二区三区乱码不卡18| 国产成人精品在线电影| 91精品国产国语对白视频| 欧美激情极品国产一区二区三区| 青青草视频在线视频观看| 免费在线观看完整版高清| 男女边吃奶边做爰视频| 免费少妇av软件| 亚洲美女视频黄频| 午夜精品国产一区二区电影| 欧美日韩精品成人综合77777| 韩国高清视频一区二区三区| 久久99精品国语久久久| 三级国产精品片| 亚洲精品视频女| 成人二区视频| 少妇的丰满在线观看| 国产午夜精品一二区理论片| 少妇猛男粗大的猛烈进出视频| 大话2 男鬼变身卡| 美女视频免费永久观看网站| 少妇的丰满在线观看| 在线观看免费视频网站a站| 日本爱情动作片www.在线观看| 九草在线视频观看| 免费日韩欧美在线观看| 欧美日韩视频精品一区| 在线观看一区二区三区激情| 9色porny在线观看| 99精国产麻豆久久婷婷| 在线免费观看不下载黄p国产| 国产熟女午夜一区二区三区| 一级爰片在线观看| 人人妻人人添人人爽欧美一区卜| 狠狠婷婷综合久久久久久88av| 国产精品国产三级国产专区5o| 精品一区二区三区四区五区乱码 | 九色亚洲精品在线播放| 日韩不卡一区二区三区视频在线| 日韩精品免费视频一区二区三区| 日本爱情动作片www.在线观看| 久久人妻熟女aⅴ| 精品国产超薄肉色丝袜足j| 亚洲国产av影院在线观看| 成人午夜精彩视频在线观看| 热re99久久国产66热| 国产精品 国内视频| 久久久久久久国产电影| 久久久国产一区二区| 久久鲁丝午夜福利片| 色婷婷久久久亚洲欧美| 日韩中字成人| 青春草视频在线免费观看| av国产精品久久久久影院| 黄频高清免费视频| 丝袜在线中文字幕| 2018国产大陆天天弄谢| 在线天堂中文资源库| 男人爽女人下面视频在线观看| 一区在线观看完整版| 日韩视频在线欧美| 免费观看无遮挡的男女| 日韩人妻精品一区2区三区| 18禁观看日本| 免费av中文字幕在线| 在线天堂中文资源库| 美女国产高潮福利片在线看| 免费在线观看黄色视频的| 日韩一区二区视频免费看| 91精品国产国语对白视频| 国产精品 国内视频| 丝瓜视频免费看黄片| 激情五月婷婷亚洲| 精品国产一区二区三区久久久樱花| 久久人人爽人人片av| 桃花免费在线播放| 在线精品无人区一区二区三| 在线 av 中文字幕| 视频区图区小说| 亚洲人成网站在线观看播放| 一级毛片我不卡| 欧美精品一区二区免费开放| 欧美97在线视频| 一二三四中文在线观看免费高清| 久久免费观看电影| 极品少妇高潮喷水抽搐| 在线观看美女被高潮喷水网站| 大话2 男鬼变身卡| 天天躁夜夜躁狠狠躁躁| 国产片内射在线| 黑丝袜美女国产一区| 丰满饥渴人妻一区二区三| 精品午夜福利在线看| 国产成人精品福利久久| 侵犯人妻中文字幕一二三四区| 女人精品久久久久毛片| 十八禁网站网址无遮挡| 水蜜桃什么品种好| 亚洲av男天堂| 十八禁网站网址无遮挡| 777久久人妻少妇嫩草av网站| 极品人妻少妇av视频| 免费在线观看完整版高清| 欧美国产精品va在线观看不卡| 欧美bdsm另类| kizo精华| 欧美国产精品va在线观看不卡| 99热国产这里只有精品6| 国产又色又爽无遮挡免| 亚洲欧美清纯卡通| 99香蕉大伊视频| 久久毛片免费看一区二区三区| 国产在视频线精品| 婷婷色综合大香蕉| av视频免费观看在线观看| 亚洲综合精品二区| 欧美日韩视频精品一区| 久久精品aⅴ一区二区三区四区 | 国产精品.久久久| 国产 精品1| 女人精品久久久久毛片| 久久鲁丝午夜福利片| 少妇 在线观看| 各种免费的搞黄视频| 国产成人精品福利久久| 欧美成人午夜精品| 自拍欧美九色日韩亚洲蝌蚪91| 久久精品人人爽人人爽视色| 激情视频va一区二区三区| 免费久久久久久久精品成人欧美视频| 亚洲成色77777| 婷婷色av中文字幕| 亚洲欧洲精品一区二区精品久久久 | 国产毛片在线视频| 国产精品久久久久久精品电影小说| 精品一区二区免费观看| 国产探花极品一区二区| 亚洲国产色片| 日本色播在线视频| 精品亚洲成a人片在线观看| 久久久精品94久久精品| 国产野战对白在线观看| 久久久亚洲精品成人影院|