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

    X射線熒光光譜結(jié)合判別分析識(shí)別進(jìn)口鐵礦石產(chǎn)地及品牌

    2020-08-08 07:39:58安雅睿朱志秀
    光譜學(xué)與光譜分析 2020年8期
    關(guān)鍵詞:判別函數(shù)鐵礦石正確率

    張 博,閔 紅,劉 曙*,安雅睿,李 晨,朱志秀

    1. 上海理工大學(xué)理學(xué)院化學(xué)系,上海 200093 2. 上海出入境檢驗(yàn)檢疫局工業(yè)品與原材料檢測(cè)技術(shù)中心,上海 200135

    引 言

    鐵礦石是鋼鐵工業(yè)的重要原材料,不同產(chǎn)地來(lái)源的鐵礦石由于地質(zhì)成因差異,主次元素含量存在一定區(qū)域特征。中國(guó)是全球最大的鐵礦石進(jìn)口國(guó),2017年進(jìn)口量超過(guò)全世界海運(yùn)鐵礦石貿(mào)易量的75%。進(jìn)口鐵礦石中不泛存在摻雜、 摻假、 以次充好的現(xiàn)象,雖然集中于個(gè)案,但對(duì)我國(guó)國(guó)門(mén)安全,經(jīng)濟(jì)安全的危害不容小覷。澳大利亞、 巴西、 南非是全球鐵礦石最主要的出口國(guó),涉及國(guó)際大型礦業(yè)集團(tuán)數(shù)十種品牌鐵礦石。品牌鐵礦石批次多、 數(shù)量大、 質(zhì)量相對(duì)穩(wěn)定。品牌鐵礦石的識(shí)別,可支撐進(jìn)口鐵礦石的風(fēng)險(xiǎn)監(jiān)管,保障貿(mào)易便利化。

    可見(jiàn)光-近紅外光譜、 微波介電光譜、 激光誘導(dǎo)擊穿光譜[1-3]結(jié)合化學(xué)計(jì)量學(xué)或機(jī)器學(xué)習(xí),可實(shí)現(xiàn)不同種類(lèi)鐵礦石的識(shí)別。波長(zhǎng)色散-X射線熒光光譜具有制樣簡(jiǎn)單、 無(wú)損分析、 穩(wěn)定性好、 靈敏度高等優(yōu)點(diǎn),能實(shí)現(xiàn)鐵礦石中主次元素的快速測(cè)定,在海關(guān)系統(tǒng)應(yīng)用非常廣泛,有利于進(jìn)口鐵礦石的快速通關(guān)。孟海東[4]、 Navid Khajehzadeh[5]等分別應(yīng)用X射線熒光光譜結(jié)合神經(jīng)網(wǎng)絡(luò),實(shí)現(xiàn)鐵礦石與銅礦石、 赤鐵礦與磁鐵礦的識(shí)別。判別分析是一種多變量統(tǒng)計(jì)分析方法,逐步判別能實(shí)現(xiàn)有效變量篩選,特別是變量間相關(guān)系數(shù)較大時(shí),能剔除不合適的變量,從而提高判別準(zhǔn)確率。武素茹等[6]以67個(gè)已知國(guó)別鐵礦石樣本X射線熒光光譜無(wú)標(biāo)樣分析數(shù)據(jù)為基礎(chǔ),采用逐步判別法篩選出CaO,MgO,Al2O3,CuO,V2O5五個(gè)特征變量,利用非參數(shù)判別方法建立進(jìn)口國(guó)別的判別模型,準(zhǔn)確率為74.6%。至目前為止,尚沒(méi)有參數(shù)判別分析方法在不同鐵礦石識(shí)別中的報(bào)道。

    澳大利亞、 巴西、 南非作為全球鐵礦石最主要的出口國(guó),主要鐵礦產(chǎn)區(qū)相對(duì)集中。如:澳大利亞90%的鐵礦石資源量和產(chǎn)量都來(lái)自于西澳洲皮爾巴拉克拉通的哈默斯利成礦省,主要鐵礦床包括芒特維爾貝克、 湯姆普萊斯山、 帕拉伯杜等,它們均產(chǎn)于元古宙早期布羅克曼BIF型含鐵建造中[7],這些礦床產(chǎn)出了PB粉、 PB塊、 楊迪粉、 紐塊、 紐粉、 津布巴粉、 國(guó)王粉(見(jiàn)表1)等知名鐵礦石品牌。由于地質(zhì)成因相似,與產(chǎn)地國(guó)別的識(shí)別相比,同一國(guó)家不同品牌鐵礦石的識(shí)別更有難度。目前亦未有進(jìn)口品牌鐵礦石識(shí)別方法的報(bào)道。

    表1 鐵礦石樣品信息Table 1 The information of iron ore samples

    在全國(guó)主要鐵礦石進(jìn)口口岸采集了來(lái)自澳大利亞、 巴西、 南非的14種品牌鐵礦石,236批進(jìn)口鐵礦石代表性樣品(見(jiàn)表1)。采用波長(zhǎng)色散-X射線熒光光譜無(wú)標(biāo)樣分析法共計(jì)檢出24種元素,選擇236批樣品全部檢出的12種元素含量用于判別分析,逐步判別法篩選出其中10種元素含量作為特征變量,采用 Fisher判別分析建立了針對(duì)進(jìn)口鐵礦石產(chǎn)地和品牌的判別模型,討論了不同品牌鐵礦石的化學(xué)成分差異,通過(guò)建模樣品驗(yàn)證、 交叉驗(yàn)證、 測(cè)試樣品驗(yàn)證,確證了模型的準(zhǔn)確性和適用性。

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

    1.1 樣品收集

    根據(jù)GB/T 10322.1—2014《鐵礦石取樣和制樣方法》,從我國(guó)主要的鐵礦石進(jìn)口口岸采集并制備來(lái)自澳大利亞、 南非、 巴西3個(gè)國(guó)家的進(jìn)口鐵礦石化學(xué)分析樣品,包含14個(gè)品牌的共計(jì)236批次樣品。樣品容量大、 種類(lèi)豐富,有一定的獨(dú)立性、 代表性,包含了我國(guó)進(jìn)口鐵礦石主要來(lái)源國(guó)及主流品牌礦種。樣品信息如表1所示,所在礦區(qū)及位置如圖1所示。

    圖1 品牌鐵礦石分布圖Fig.1 Distribution map of brand iron ore

    1.2 方法

    將樣品分裝到干燥瓶中于105 ℃下烘干4 h。采用壓片機(jī)對(duì)烘干樣品進(jìn)行壓片,壓片前用乙醇清洗模具,使用聚乙烯環(huán)使粉末樣品聚攏,壓制樣品在30 t壓力下維持30 s。檢查壓制樣品表面均勻且無(wú)裂紋、 脫落現(xiàn)象,測(cè)量前用洗耳球吹凈樣品表面。

    使用德國(guó)布魯克公司S4 Pioneer波長(zhǎng)色散-X射線熒光光譜儀中的無(wú)標(biāo)樣分析方法檢測(cè)鐵礦石中元素的含量。無(wú)標(biāo)樣分析也稱(chēng)半定量分析,它的基本思路是由儀器和軟件制造商測(cè)定校準(zhǔn)樣品,儲(chǔ)存元素譜線強(qiáng)度和校準(zhǔn)曲線,然后將這些數(shù)據(jù)轉(zhuǎn)到用戶的X射線熒光分析系統(tǒng)中,并用參考樣品校正儀器的漂移,無(wú)標(biāo)樣分析法最大的優(yōu)點(diǎn)是快速。檢測(cè)中使用銠靶光管,四個(gè)分析儀晶體(LiF200,XS-55,PET和Ge),流氣計(jì)數(shù)器(FC),閃爍計(jì)數(shù)器(SC)等元件。表2列出了儀器的部分測(cè)量條件。

    表2 儀器測(cè)量條件Table 2 Instrument measurement conditions

    1.3 逐步判別-Fisher判別分析

    逐步判別分析屬于有監(jiān)督的分類(lèi)方式,先對(duì)已知的樣品進(jìn)行分類(lèi)來(lái)建立模型,再對(duì)未知樣品來(lái)進(jìn)行預(yù)測(cè)分類(lèi)。其中一個(gè)重要的思想就是逐步引入變量,每次只引入一個(gè)變量,同時(shí)也檢驗(yàn)先前引入的變量;如果引入的新變量導(dǎo)致之前變量的判別能力不再顯著,就將先前引入的變量從判別式中移除,篩選至判別式中的變量都很顯著[8]。

    Fisher判別的基本原理是投影,將高維數(shù)據(jù)投影到某個(gè)方向,使得組與組之間區(qū)別最大,組內(nèi)的區(qū)別最小,其判別函數(shù)的建立利用了方差分析的思想[9]。根據(jù)該原則確定判別式,將一未知樣品的變量代入判別式。根據(jù)判別函數(shù)和組質(zhì)心處坐標(biāo)函數(shù),計(jì)算每個(gè)樣品坐標(biāo)與質(zhì)心的距離,與哪個(gè)類(lèi)別的質(zhì)心最近,該樣品就判定為哪個(gè)類(lèi)別。

    本文分析來(lái)自我國(guó)主要鐵礦石進(jìn)口口岸的14個(gè)品牌共計(jì)236個(gè)鐵礦石樣品,建模過(guò)程中選取191個(gè)樣品作為訓(xùn)練集,45個(gè)樣品作為驗(yàn)證集檢驗(yàn)?zāi)P偷臏?zhǔn)確性。訓(xùn)練樣品及測(cè)試樣品的選取如表1所示。

    2 結(jié)果與討論

    2.1 元素的選擇

    針對(duì)采集的236個(gè)鐵礦石樣品,采用波長(zhǎng)色散-X射線熒光光譜無(wú)標(biāo)樣分析共計(jì)檢出Fe,O,Si,Ca,Al,Mn,Tb,Ti,Mg,P,K,S,Cr,Na,Sr,Zr,Zn,V,Cu,Gd,Ba,Cl,Ni和Co共24種元素,其中K,Cu,Zr,Zn,Na,Cl,V,Sr,Gd,Ni,Ba和Co共12個(gè)元素含量存在未檢出的情況,未檢出比例分別為18.20%,50.00%,51.00%,69.90%,70.30%,73.30%,78.00%,83.90%,84.30%,91.50%,92.80%和97.00%,建立鐵礦產(chǎn)地及品牌的識(shí)別模型,在滿足實(shí)際應(yīng)用的前提下,應(yīng)選擇鐵礦樣品檢出比例盡量高的元素,本文選取236個(gè)樣品全部檢出的Fe,O,Si,Ca,Al,Mn,Tb,Ti,Mg,P,Cr和S共12種元素含量用于后續(xù)分析。

    針對(duì)不同進(jìn)口國(guó)家、 品牌鐵礦石的模式識(shí)別,采用逐步判別分析對(duì)Fe,O,Si,Ca,Al,Mn,Tb,Ti,Mg,P,Cr和S共12個(gè)元素含量進(jìn)行變量篩選,變量能否進(jìn)入模型主要取決于協(xié)方差分析的F檢驗(yàn)的顯著性水平,當(dāng)F值大于指定值時(shí)保留該變量,而F值小于指定值時(shí),該變量從模型中剔除。選取合適的F值可以用最少的變量達(dá)到最佳的判別效果。本文選取的F值為3.84,經(jīng)過(guò)逐步判別分析,F(xiàn)e,O,Si,Ca,Al,Mn,Ti,Mg,P和S共10個(gè)元素保留在了模型中,Tb與Cr因未通過(guò)F檢驗(yàn)(F值<3.84)而從模型中剔除[10],最終10個(gè)元素用于建立識(shí)別模型。

    14種品牌鐵礦石12個(gè)元素(Fe,O,Si,Ca,Al,Mn,Tb,Ti,Mg,P,Cr,S)含量的均值對(duì)比分析(圖2)表明:澳精粉Fe,Si和O含量與其他類(lèi)別有顯著差異,巴西粉和混合粉的Si,Mn含量與其他類(lèi)別有顯著差異,南非精粉Ca,Ti,Mg,P和S的含量明顯高于其他類(lèi)別,不同類(lèi)別鐵礦石之間的Al,Mn,Mg,P和S的含量也存在明顯的差異。因此可以利用不同元素的含量組合建立線性判別模型,對(duì)鐵礦石進(jìn)口國(guó)別、 品牌進(jìn)行識(shí)別。不同品牌鐵礦石Tb和Cr含量的平均值在直方圖上差異性不大,這也解釋了逐級(jí)判別分析將這兩個(gè)元素剔除的原因。

    圖2 鐵礦石樣本的元素平均含量條形圖(a):Fe,O;(b):Si,Al;(c):Ca,Mn;(d):Mg,Ti;(e):P,Tb;(f):S,CrFig.2 Bar chart of element average content in iron ore samples(a):Fe,O;(b):Si,Al;(c):Ca,Mn;(d):Mg,Ti;(e):P,Tb;(f):S,Cr

    2.2 識(shí)別進(jìn)口鐵礦石國(guó)別

    使用Fe,O,Si,Ca,Al,Mn,Ti,Mg,P和S共10個(gè)元素含量建立Fisher判別模型,得到2組判別函數(shù)和相應(yīng)的組質(zhì)心處的坐標(biāo)。判別函數(shù)與各變量之間的相關(guān)性如圖3所示,橫坐標(biāo)為函數(shù)1(F1)與各變量的相關(guān)系數(shù),縱坐標(biāo)為函數(shù)2(F2)與各變量的相關(guān)系數(shù),系數(shù)為正表示正相關(guān),系數(shù)為負(fù)表示負(fù)相關(guān),絕對(duì)值越大相關(guān)性越高。Ca,O,Ti,Mg,P各元素含量與函數(shù)1(F1)的相關(guān)系數(shù)分別為0.277,-0.236,0.234,0.230,0.193,Mn,O,F(xiàn)e,S各元素含量與函數(shù)2(F2)的相關(guān)系數(shù)分別為-0.279,0.268,-0.181,0.160,是相關(guān)性相對(duì)較大的元素。

    圖3 變量與判別函數(shù)間的相關(guān)性Fig.3 Correlation between variables and discriminant functions

    判別函數(shù):

    F1=0.525X1-0.598X2+1.4X3+32.627X4+0.654X5-3.936X6+37.01X7-29.4X8-58.953X9-24.002X10-16.337

    F2=0.569X1+0.855X2+0.122X3+7.559X4+1.23X5-4.789X6-9.846X7+4.281X8-128.56X9+147.622X10-61.555

    式中X1—X10分別代表Fe,O,Si,Ca,Al,Mn,Ti,Mg,P和S的含量。

    用函數(shù)1和函數(shù)2的判別得分作散點(diǎn)圖(圖4),橫坐標(biāo)為函數(shù)1得分,縱坐標(biāo)為函數(shù)2得分,可以看出模型對(duì)南非鐵礦石和澳大利亞,巴西兩個(gè)國(guó)家的鐵礦石區(qū)分明顯,澳大利亞與巴西散點(diǎn)存在重疊交叉的現(xiàn)象。

    圖4 判別函數(shù)得分散點(diǎn)圖Fig.4 Discriminant function decentralized point map

    建模樣品為構(gòu)建模型所用的樣品,可以回代到模型,驗(yàn)證模型識(shí)別的準(zhǔn)確性。交叉驗(yàn)證是一種重要的判別效果驗(yàn)證方法,該法可以非常有效地避免強(qiáng)影響點(diǎn)的干擾。本文采用留一交叉驗(yàn)證法對(duì)建模所用的樣品進(jìn)行驗(yàn)證,每次留出一個(gè)作為驗(yàn)證,其余用來(lái)建模,最后計(jì)算平均準(zhǔn)確率作為對(duì)模型的評(píng)價(jià)。測(cè)試樣品為建模過(guò)程中預(yù)留的用于測(cè)試模型識(shí)別正確率的樣品。該模型對(duì)建模樣品分類(lèi)正確率為97.40%,對(duì)南非的鐵礦石樣品識(shí)別正確率為100%,對(duì)澳大利亞、 巴西鐵礦石樣品存在識(shí)別錯(cuò)誤的情況,正確率分別為97.40%,91.70%。模型交叉驗(yàn)證正確率為95.30%,交叉驗(yàn)證的正確率高于80%,說(shuō)明該模型有很好的分類(lèi)準(zhǔn)確度。為了確定模型是否可以對(duì)未包含在模型中的樣品進(jìn)行識(shí)別,分析了建模時(shí)選擇的45個(gè)作為測(cè)試樣品的鐵礦石樣品,模型對(duì)測(cè)試樣品識(shí)別正確率達(dá)到95.50%,其中對(duì)南非和巴西樣品識(shí)別正確度都達(dá)到100%,說(shuō)明此模型可以對(duì)鐵礦石的國(guó)別進(jìn)行很好的識(shí)別[11]。

    2.3 識(shí)別進(jìn)口鐵礦石品牌

    與進(jìn)口鐵礦石國(guó)別的識(shí)別相比較,進(jìn)口鐵礦石品牌的識(shí)別更加具有難度,因?yàn)椴煌放频蔫F礦石可能會(huì)來(lái)自相同國(guó)家相同的礦區(qū),它們物相結(jié)構(gòu)一致,元素含量的差異性也有可能不大。在對(duì)鐵礦石國(guó)別已能進(jìn)行很好識(shí)別的基礎(chǔ)上,嘗試對(duì)鐵礦石品牌做進(jìn)一步識(shí)別。采用Fe,O,Si,Ca,Al,Mn,Ti,Mg,P,S10個(gè)元素含量建立Fisher判別模型,得到10個(gè)判別函數(shù)和相應(yīng)的組質(zhì)心處的坐標(biāo)。

    判別函數(shù):

    F1=0.286X1-0.372X2+0.404X3-9.8X4-0.575X5-2.551X6+54.081X7+18.203X8-95.378X9-17.295X10-12.937

    F2=0.726X1-1.244X2+2.718X3-2.715X4-0.808X5-5.169X6-23.321X7+15.246X8-160.116X9-18.187X10-2.608

    F3=0.587X1-1.372X2+1.45X3+8.922X4+3.122X5+3.139X6+9.551X7-15.179X8+180.272X9+21.586X10-4.251

    F4=0.257X1-0.263X2+0.896X3+33.291X4-1.573X5+1.075X6+37.848X7-34.022X8-62.815X9+25.649X10-6.865

    F5=-0.4X1-0.453X2+1.121X3-13.49X4-1.417X5+11.016X6+9.278X7+0.466X8+19.539X9+87.158X10+35.004

    F6=-0.012X1+0.213X2+1.171X3+6.748X4+2.843X5-6.586X6+12.703X7-2.634X8-138.53X9-59.512X10-9.303

    F7=-0.189X1+0.221X2+2.048X3-4.649X4-3.679X5-0.974X6+10.841X7-7.985X8+221.911X9-69.058X10-0.463

    F8=0.337X1+0.494X2+0.457X3+7.023X4-0.625X5-1.537X6-17.619X7+2.087X8+6.88X9+229.064X10-38.458

    F9=0.518X1-0.035X2+0.885X3-14.995X4-0.257X5-5.636X6+26.099X7-4.9X8-30.377X9+159.188X10-31.76

    F10=1.526X1+1.558X2+1.681X3-1.042X4-0.188X5+4.851X6+6.292X7+0.682X8+7.995X9-75.431X10-145.479

    式中X1—X10分別代表Fe,O,Si,Ca,Al,Mn,Ti,Mg,P,S的含量。

    前3個(gè)判別函數(shù)(F1,F(xiàn)2,F(xiàn)3)分別解釋了總信息的90.6%,5.7%,2.0%,累計(jì)解釋98.40% ,用前3個(gè)函數(shù)建立判別模型,并用判別得分來(lái)繪制三維散點(diǎn)圖(圖5),F(xiàn)1,F(xiàn)2,F(xiàn)3分別為判別函數(shù)F1,F(xiàn)2,F(xiàn)3的得分。三維散點(diǎn)圖分布表明,14個(gè)品牌的鐵礦石可明顯地被劃分為四個(gè)區(qū)域,澳精粉和南非精粉與其他類(lèi)別區(qū)分最為明顯。從圖中還可以看出PB塊與紐曼塊分類(lèi)略有重疊,楊迪粉與國(guó)王粉的集群非常接近。所建立判別模型的分類(lèi)正確率如表3所示。結(jié)果表明:模型對(duì)澳精粉與南非精粉識(shí)別完全正確,因?yàn)榕c其他類(lèi)別的鐵礦石相比,這兩類(lèi)的Ti與Mg的含量與其他類(lèi)別有明顯不同。模型對(duì)于PB塊,楊迪粉,紐塊會(huì)存在識(shí)別錯(cuò)誤的情況。這三個(gè)品牌的鐵礦石都產(chǎn)于澳大利亞皮爾巴拉地區(qū)的哈默斯利鐵礦帶,礦石成因類(lèi)似,元素含量比較接近,因此相對(duì)于其他類(lèi)別更難以區(qū)分。

    圖5 判別函數(shù)得分三維散點(diǎn)圖Fig.5 Discriminant function score three-dimensional scatter plot

    為追求更高的識(shí)別準(zhǔn)確率,選擇使用全部10個(gè)函數(shù)來(lái)建立判別模型。所建立判別模型的分類(lèi)正確率如表3所示。結(jié)果表明:當(dāng)使用全部10個(gè)函數(shù)建模時(shí),模型對(duì)測(cè)試樣品識(shí)別的準(zhǔn)確率有明顯提高,正確率達(dá)到了100%,所建立識(shí)別模型具有很好識(shí)別效果。

    表3 判別模型的預(yù)測(cè)正確率(%)Table 3 Prediction accuracy of discriminant model(%)

    3 結(jié) 論

    利用波長(zhǎng)色散-X射線熒光光譜無(wú)標(biāo)樣分析法測(cè)定澳大利亞、 南非、 巴西3個(gè)國(guó)家14個(gè)品牌236份鐵礦石樣品的元素含量,選擇191個(gè)樣品作為訓(xùn)練樣本,45個(gè)樣品作為測(cè)試樣本,采用逐步判別分析篩選出Fe,O,Si,Ca,Al,Mn,Ti,Mg,P和S共10種元素含量作為特征變量,建立了識(shí)別鐵礦石產(chǎn)地、 品牌的Fisher判別模型。該模型為鐵礦石品牌與元素之間的關(guān)系提供了基礎(chǔ)數(shù)據(jù)與理論依據(jù),僅通過(guò)無(wú)標(biāo)樣分析法測(cè)量鐵礦石樣品10種元素含量建立判別模型,便可以快速進(jìn)行產(chǎn)地、 品牌的識(shí)別。當(dāng)然,模型樣品產(chǎn)地和品牌的確證和樣本數(shù)量是限制模型正確率的關(guān)鍵要素,當(dāng)樣本數(shù)量達(dá)到一定數(shù)量級(jí)之后,所建立識(shí)別模型的準(zhǔn)確率和普適性將得到進(jìn)一步的提升。

    猜你喜歡
    判別函數(shù)鐵礦石正確率
    6.58億噸
    門(mén)診分診服務(wù)態(tài)度與正確率對(duì)護(hù)患關(guān)系的影響
    游樂(lè)設(shè)施事故與危險(xiǎn)量化判別函數(shù)的構(gòu)建
    探究上市公司財(cái)務(wù)預(yù)警的數(shù)學(xué)模型
    生意
    品管圈活動(dòng)在提高介入手術(shù)安全核查正確率中的應(yīng)用
    生意
    鐵礦石中全鐵測(cè)定的探討
    基于Fisher判別函數(shù)的酒店員工離職預(yù)警研究
    中國(guó)鐵礦石需求增長(zhǎng)潛力有限
    金屬礦山(2013年7期)2013-03-11 16:54:15
    久久精品久久久久久噜噜老黄| 日韩伦理黄色片| av.在线天堂| 宅男免费午夜| 久久精品国产a三级三级三级| 中文字幕精品免费在线观看视频| 久久久国产一区二区| 少妇人妻 视频| 精品久久久久久电影网| 91精品三级在线观看| 中文乱码字字幕精品一区二区三区| 黄片播放在线免费| 中国国产av一级| 成人三级做爰电影| 777米奇影视久久| 一本色道久久久久久精品综合| 午夜老司机福利片| 亚洲美女视频黄频| 亚洲国产精品一区三区| 热re99久久国产66热| 免费日韩欧美在线观看| 青春草视频在线免费观看| 在线天堂最新版资源| 久久久久精品久久久久真实原创| 午夜福利视频精品| 美女视频免费永久观看网站| 中文字幕人妻丝袜一区二区 | 一级毛片我不卡| 午夜精品国产一区二区电影| 久久久国产精品麻豆| 哪个播放器可以免费观看大片| 国产视频首页在线观看| 最近2019中文字幕mv第一页| 男人添女人高潮全过程视频| 天天操日日干夜夜撸| 18禁裸乳无遮挡动漫免费视频| 亚洲国产最新在线播放| 亚洲国产毛片av蜜桃av| 午夜福利视频在线观看免费| 母亲3免费完整高清在线观看| 色网站视频免费| 国产在线视频一区二区| 久久精品熟女亚洲av麻豆精品| 在线观看一区二区三区激情| 极品少妇高潮喷水抽搐| 伊人久久国产一区二区| av线在线观看网站| 精品人妻一区二区三区麻豆| 午夜福利视频在线观看免费| 日韩欧美一区视频在线观看| 亚洲欧美色中文字幕在线| 成人午夜精彩视频在线观看| 18禁国产床啪视频网站| 一级毛片黄色毛片免费观看视频| 免费高清在线观看视频在线观看| 丰满迷人的少妇在线观看| 少妇 在线观看| 国产精品国产av在线观看| 日韩一本色道免费dvd| 国产激情久久老熟女| 老司机亚洲免费影院| 搡老岳熟女国产| 肉色欧美久久久久久久蜜桃| 熟女少妇亚洲综合色aaa.| 我要看黄色一级片免费的| 日本91视频免费播放| 99精国产麻豆久久婷婷| 亚洲精品国产一区二区精华液| 超碰成人久久| 欧美av亚洲av综合av国产av | 一区二区三区四区激情视频| 久久精品国产a三级三级三级| 亚洲成色77777| 亚洲欧洲国产日韩| 日韩免费高清中文字幕av| 国产深夜福利视频在线观看| 黄频高清免费视频| 国产精品人妻久久久影院| 夜夜骑夜夜射夜夜干| 新久久久久国产一级毛片| 亚洲成人av在线免费| 性少妇av在线| 国产 精品1| 午夜福利乱码中文字幕| 国产一区二区三区av在线| 黑人欧美特级aaaaaa片| 亚洲三区欧美一区| 五月开心婷婷网| 亚洲欧洲国产日韩| 老司机影院成人| 美女中出高潮动态图| 亚洲精品aⅴ在线观看| 一级片'在线观看视频| 老司机亚洲免费影院| 亚洲精品久久成人aⅴ小说| av网站免费在线观看视频| 久久久国产一区二区| 国产精品一区二区精品视频观看| 热99久久久久精品小说推荐| 91aial.com中文字幕在线观看| 91老司机精品| 国产乱人偷精品视频| 最近中文字幕2019免费版| 成年美女黄网站色视频大全免费| av一本久久久久| 日韩一区二区三区影片| 久久久国产精品麻豆| 一级毛片 在线播放| 国产在视频线精品| 亚洲人成网站在线观看播放| 欧美日韩福利视频一区二区| av不卡在线播放| 国产色婷婷99| 街头女战士在线观看网站| 人人妻人人添人人爽欧美一区卜| 国产在线一区二区三区精| 叶爱在线成人免费视频播放| 欧美国产精品一级二级三级| 在线免费观看不下载黄p国产| 蜜桃在线观看..| 国产精品成人在线| 国产免费一区二区三区四区乱码| 青春草视频在线免费观看| 日韩av在线免费看完整版不卡| 国产一区二区 视频在线| 妹子高潮喷水视频| 丝袜美腿诱惑在线| www.精华液| 最近的中文字幕免费完整| 深夜精品福利| 视频在线观看一区二区三区| www.熟女人妻精品国产| 天天躁夜夜躁狠狠躁躁| 亚洲一码二码三码区别大吗| 成人18禁高潮啪啪吃奶动态图| 晚上一个人看的免费电影| 国产精品av久久久久免费| 中文精品一卡2卡3卡4更新| 欧美日韩精品网址| 久久精品国产亚洲av涩爱| 精品一区二区三区四区五区乱码 | 女性被躁到高潮视频| 精品人妻在线不人妻| 亚洲精品久久久久久婷婷小说| 黄片播放在线免费| 国产深夜福利视频在线观看| 十八禁网站网址无遮挡| 成年动漫av网址| 在线观看免费视频网站a站| 日韩一区二区视频免费看| 九草在线视频观看| 色精品久久人妻99蜜桃| 99热国产这里只有精品6| 男女高潮啪啪啪动态图| 久久国产亚洲av麻豆专区| 亚洲四区av| 性色av一级| 国产乱人偷精品视频| 一级爰片在线观看| 老汉色av国产亚洲站长工具| 精品国产乱码久久久久久男人| 亚洲色图 男人天堂 中文字幕| 亚洲色图 男人天堂 中文字幕| 成人影院久久| 久久久久久人妻| 大话2 男鬼变身卡| xxxhd国产人妻xxx| 亚洲精品久久成人aⅴ小说| 在线观看免费日韩欧美大片| av不卡在线播放| 国产又色又爽无遮挡免| 少妇被粗大猛烈的视频| 国产精品99久久99久久久不卡 | 久久 成人 亚洲| 亚洲成人免费av在线播放| 男的添女的下面高潮视频| 国产高清国产精品国产三级| 男女午夜视频在线观看| 青草久久国产| 老司机亚洲免费影院| 国产毛片在线视频| 2018国产大陆天天弄谢| 飞空精品影院首页| 欧美精品一区二区大全| 热re99久久国产66热| 日韩一区二区视频免费看| 中文字幕人妻熟女乱码| 一级a爱视频在线免费观看| 在线精品无人区一区二区三| 大片电影免费在线观看免费| 乱人伦中国视频| 母亲3免费完整高清在线观看| xxxhd国产人妻xxx| 精品人妻在线不人妻| av福利片在线| 一边摸一边抽搐一进一出视频| 高清视频免费观看一区二区| 又粗又硬又长又爽又黄的视频| 丰满饥渴人妻一区二区三| 国产精品一区二区在线观看99| 大片电影免费在线观看免费| 中文字幕亚洲精品专区| 日韩av不卡免费在线播放| 国产1区2区3区精品| 亚洲欧美色中文字幕在线| 亚洲美女黄色视频免费看| 欧美日韩国产mv在线观看视频| 国产精品久久久人人做人人爽| 99热国产这里只有精品6| 9热在线视频观看99| 欧美日本中文国产一区发布| 免费av中文字幕在线| 亚洲美女搞黄在线观看| 欧美精品一区二区免费开放| 九色亚洲精品在线播放| 国产有黄有色有爽视频| 91老司机精品| 80岁老熟妇乱子伦牲交| av又黄又爽大尺度在线免费看| 青春草国产在线视频| 国产高清不卡午夜福利| 黄色怎么调成土黄色| 一级爰片在线观看| 久久久久精品久久久久真实原创| 啦啦啦啦在线视频资源| 亚洲欧美一区二区三区黑人| 美女福利国产在线| 国产毛片在线视频| 2018国产大陆天天弄谢| 午夜福利视频在线观看免费| 日本一区二区免费在线视频| 黄网站色视频无遮挡免费观看| 精品国产乱码久久久久久男人| 黄片播放在线免费| 欧美变态另类bdsm刘玥| 色综合欧美亚洲国产小说| 婷婷成人精品国产| 麻豆av在线久日| 大片免费播放器 马上看| 亚洲欧美成人精品一区二区| 色综合欧美亚洲国产小说| 国产精品久久久久久久久免| 欧美日韩成人在线一区二区| 国产 精品1| 黑丝袜美女国产一区| 搡老乐熟女国产| 丁香六月天网| 亚洲国产av新网站| 午夜精品国产一区二区电影| 日本av免费视频播放| 亚洲av日韩精品久久久久久密 | 搡老乐熟女国产| tube8黄色片| 69精品国产乱码久久久| 麻豆精品久久久久久蜜桃| 51午夜福利影视在线观看| av在线播放精品| 另类亚洲欧美激情| 国产97色在线日韩免费| 丝瓜视频免费看黄片| 日韩一本色道免费dvd| 人妻人人澡人人爽人人| 女人久久www免费人成看片| 久久久久精品国产欧美久久久 | 欧美国产精品va在线观看不卡| 日韩,欧美,国产一区二区三区| 亚洲欧洲日产国产| 精品亚洲乱码少妇综合久久| av在线观看视频网站免费| 高清欧美精品videossex| 国产一区二区三区综合在线观看| 各种免费的搞黄视频| 亚洲精品第二区| 精品免费久久久久久久清纯 | 啦啦啦在线免费观看视频4| 久热这里只有精品99| 国产一区二区在线观看av| 亚洲综合精品二区| 精品视频人人做人人爽| 菩萨蛮人人尽说江南好唐韦庄| 亚洲成人手机| 欧美黑人精品巨大| 2018国产大陆天天弄谢| 日本欧美视频一区| 在线观看免费午夜福利视频| 亚洲成人一二三区av| av网站免费在线观看视频| 美国免费a级毛片| www.av在线官网国产| 在线观看国产h片| 一级爰片在线观看| 一本色道久久久久久精品综合| 狠狠精品人妻久久久久久综合| 丝袜人妻中文字幕| 老司机影院成人| 天天添夜夜摸| 美女福利国产在线| 亚洲精品美女久久av网站| av天堂久久9| 国产xxxxx性猛交| 亚洲欧美成人综合另类久久久| 天堂俺去俺来也www色官网| 国产福利在线免费观看视频| 午夜福利影视在线免费观看| 亚洲人成77777在线视频| 亚洲精品久久午夜乱码| 高清视频免费观看一区二区| 两个人看的免费小视频| 欧美激情 高清一区二区三区| 高清av免费在线| 涩涩av久久男人的天堂| 操美女的视频在线观看| 极品少妇高潮喷水抽搐| 七月丁香在线播放| 日日啪夜夜爽| 精品一区二区三区四区五区乱码 | 国产精品麻豆人妻色哟哟久久| bbb黄色大片| 婷婷成人精品国产| 一区在线观看完整版| 亚洲熟女毛片儿| 永久免费av网站大全| 9191精品国产免费久久| 亚洲精品美女久久av网站| 黑人欧美特级aaaaaa片| 男女床上黄色一级片免费看| 制服人妻中文乱码| 久久av网站| 国产精品欧美亚洲77777| 国产日韩欧美在线精品| 国产成人欧美在线观看 | 99精国产麻豆久久婷婷| 一区二区三区乱码不卡18| 成人亚洲精品一区在线观看| 99久久综合免费| 女人久久www免费人成看片| 日本爱情动作片www.在线观看| 欧美老熟妇乱子伦牲交| 久久久久人妻精品一区果冻| 尾随美女入室| 国产精品女同一区二区软件| 十八禁网站网址无遮挡| 一区二区三区四区激情视频| av有码第一页| 久久韩国三级中文字幕| 亚洲人成77777在线视频| 卡戴珊不雅视频在线播放| 亚洲成av片中文字幕在线观看| 丰满乱子伦码专区| 亚洲精品久久成人aⅴ小说| 又粗又硬又长又爽又黄的视频| 国产一区有黄有色的免费视频| 日本vs欧美在线观看视频| 亚洲图色成人| 99热网站在线观看| 国产精品麻豆人妻色哟哟久久| 国产97色在线日韩免费| 欧美精品高潮呻吟av久久| 免费不卡黄色视频| 欧美av亚洲av综合av国产av | 亚洲视频免费观看视频| 国产亚洲最大av| 日韩欧美一区视频在线观看| 久久ye,这里只有精品| 热99久久久久精品小说推荐| 两个人看的免费小视频| 99国产精品免费福利视频| 国产成人av激情在线播放| 亚洲视频免费观看视频| 青青草视频在线视频观看| 亚洲男人天堂网一区| 制服人妻中文乱码| 中文字幕人妻丝袜一区二区 | 国产熟女欧美一区二区| 超碰成人久久| 久久久久久久久久久久大奶| 亚洲精品第二区| 这个男人来自地球电影免费观看 | 午夜老司机福利片| 丝瓜视频免费看黄片| 少妇人妻 视频| 亚洲图色成人| 最新在线观看一区二区三区 | 在线观看人妻少妇| 啦啦啦在线观看免费高清www| 少妇精品久久久久久久| 中文精品一卡2卡3卡4更新| 青春草视频在线免费观看| 日本猛色少妇xxxxx猛交久久| avwww免费| 黑人猛操日本美女一级片| 两个人免费观看高清视频| 一边亲一边摸免费视频| 亚洲精华国产精华液的使用体验| 精品一区在线观看国产| 高清视频免费观看一区二区| 一二三四中文在线观看免费高清| 最近中文字幕2019免费版| 一本久久精品| 在线看a的网站| 中文字幕精品免费在线观看视频| 亚洲精品aⅴ在线观看| 亚洲欧美一区二区三区国产| 国产黄频视频在线观看| 飞空精品影院首页| 性少妇av在线| 欧美精品高潮呻吟av久久| 中文精品一卡2卡3卡4更新| 午夜福利乱码中文字幕| 久久影院123| 99re6热这里在线精品视频| 国产成人av激情在线播放| 日本黄色日本黄色录像| 99九九在线精品视频| av.在线天堂| 国产老妇伦熟女老妇高清| a级毛片在线看网站| 亚洲欧美激情在线| 别揉我奶头~嗯~啊~动态视频 | 亚洲精品第二区| 91国产中文字幕| 韩国高清视频一区二区三区| 国产极品粉嫩免费观看在线| 无遮挡黄片免费观看| 午夜日韩欧美国产| 香蕉丝袜av| 国产片特级美女逼逼视频| 精品福利永久在线观看| 在线观看www视频免费| 欧美日韩视频精品一区| 成年av动漫网址| 女人精品久久久久毛片| 成年动漫av网址| 欧美精品av麻豆av| 成年人免费黄色播放视频| av卡一久久| 美女脱内裤让男人舔精品视频| 老司机影院成人| 午夜日本视频在线| 国产精品蜜桃在线观看| 精品国产一区二区久久| 久久国产精品男人的天堂亚洲| 一级片'在线观看视频| 亚洲国产精品成人久久小说| 日韩一本色道免费dvd| 精品国产超薄肉色丝袜足j| 日本色播在线视频| 亚洲图色成人| 欧美日韩亚洲国产一区二区在线观看 | av天堂久久9| 欧美老熟妇乱子伦牲交| 在线观看免费高清a一片| 亚洲成人av在线免费| 亚洲自偷自拍图片 自拍| 老司机影院毛片| 一二三四中文在线观看免费高清| 国产av码专区亚洲av| a 毛片基地| 亚洲精品,欧美精品| 三上悠亚av全集在线观看| 色吧在线观看| 人妻 亚洲 视频| 亚洲国产精品一区三区| 亚洲成色77777| 老司机影院成人| 日韩大码丰满熟妇| 日韩制服丝袜自拍偷拍| 亚洲国产欧美在线一区| kizo精华| 卡戴珊不雅视频在线播放| 久久人妻熟女aⅴ| 久热爱精品视频在线9| 波多野结衣一区麻豆| 极品少妇高潮喷水抽搐| 黑人猛操日本美女一级片| 国产一区二区三区综合在线观看| 中文字幕精品免费在线观看视频| 这个男人来自地球电影免费观看 | 国产一区亚洲一区在线观看| av视频免费观看在线观看| 亚洲 欧美一区二区三区| 亚洲成人av在线免费| 欧美日韩av久久| 国产一区二区三区av在线| av在线app专区| 又粗又硬又长又爽又黄的视频| 熟女av电影| 久久韩国三级中文字幕| 亚洲精品aⅴ在线观看| 免费看av在线观看网站| 亚洲欧洲国产日韩| 日韩电影二区| 黄色怎么调成土黄色| 国产有黄有色有爽视频| 国产亚洲av高清不卡| 视频区图区小说| 男人添女人高潮全过程视频| 人人妻人人添人人爽欧美一区卜| 免费观看人在逋| 少妇人妻精品综合一区二区| 国产国语露脸激情在线看| 精品一区二区三卡| 亚洲精品久久午夜乱码| 成年美女黄网站色视频大全免费| 久久久精品区二区三区| 中国三级夫妇交换| 蜜桃在线观看..| 99久久综合免费| 男的添女的下面高潮视频| 如何舔出高潮| 日韩中文字幕视频在线看片| 黄片小视频在线播放| 午夜日韩欧美国产| 日韩av在线免费看完整版不卡| 纯流量卡能插随身wifi吗| 日韩中文字幕视频在线看片| 国产成人一区二区在线| 日韩欧美精品免费久久| 日韩一本色道免费dvd| 亚洲三区欧美一区| 两性夫妻黄色片| 国产一区亚洲一区在线观看| 国产精品一区二区精品视频观看| 丝袜美腿诱惑在线| 欧美乱码精品一区二区三区| 波野结衣二区三区在线| 亚洲成人手机| 免费高清在线观看日韩| 99久久综合免费| 欧美日韩亚洲综合一区二区三区_| 99九九在线精品视频| 国产欧美亚洲国产| 在线观看免费高清a一片| 女性生殖器流出的白浆| 中文字幕高清在线视频| 亚洲中文av在线| 亚洲国产欧美网| 99精国产麻豆久久婷婷| 最近手机中文字幕大全| 国产免费福利视频在线观看| 欧美国产精品va在线观看不卡| 国产女主播在线喷水免费视频网站| 最近的中文字幕免费完整| 亚洲美女黄色视频免费看| 美女视频免费永久观看网站| 观看av在线不卡| 亚洲精品中文字幕在线视频| 少妇被粗大猛烈的视频| 国产在线一区二区三区精| 成人国产av品久久久| 国产精品久久久久久久久免| 欧美日韩一级在线毛片| 黑人巨大精品欧美一区二区蜜桃| 精品一区二区免费观看| 免费黄频网站在线观看国产| 中文字幕最新亚洲高清| 国产有黄有色有爽视频| 日韩一区二区三区影片| 亚洲欧美成人综合另类久久久| 国产日韩欧美亚洲二区| av网站在线播放免费| 黄色毛片三级朝国网站| 人妻 亚洲 视频| 男女高潮啪啪啪动态图| 免费在线观看完整版高清| 亚洲欧洲精品一区二区精品久久久 | 欧美精品亚洲一区二区| 久久久亚洲精品成人影院| 韩国av在线不卡| 一区在线观看完整版| 国产深夜福利视频在线观看| 一本色道久久久久久精品综合| 青春草视频在线免费观看| av天堂久久9| 亚洲国产日韩一区二区| 嫩草影院入口| 香蕉丝袜av| 欧美日韩精品网址| 久久99精品国语久久久| 纯流量卡能插随身wifi吗| 极品人妻少妇av视频| 日本欧美视频一区| 国产亚洲最大av| 亚洲精品美女久久av网站| 久久精品国产综合久久久| 男女免费视频国产| 免费黄色在线免费观看| av免费观看日本| 精品亚洲成a人片在线观看| 亚洲av福利一区| 亚洲,欧美精品.| 久久精品亚洲熟妇少妇任你| 又大又黄又爽视频免费| 宅男免费午夜| 五月天丁香电影| 久久久精品国产亚洲av高清涩受| 久久av网站| 亚洲国产精品成人久久小说| 少妇人妻久久综合中文| 中文字幕人妻丝袜一区二区 | 午夜av观看不卡| 美女中出高潮动态图| 好男人视频免费观看在线| 国产一区二区三区av在线| 一区二区日韩欧美中文字幕| 日韩欧美精品免费久久| 婷婷色综合大香蕉| 亚洲视频免费观看视频| 在现免费观看毛片| 国产精品秋霞免费鲁丝片| 18禁国产床啪视频网站| 久久精品亚洲熟妇少妇任你| 国产精品秋霞免费鲁丝片| 中文字幕人妻丝袜一区二区 | 男人操女人黄网站|