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

    噴氣燃料近紅外光譜分析的整體不確定度評定

    2023-07-03 02:24:00邢志娜
    兵器裝備工程學(xué)報 2023年6期
    關(guān)鍵詞:紅外光譜儀噴氣分量

    邢志娜,袁 林

    (1.海軍航空大學(xué), 山東 煙臺 264001;2. 92154部隊, 山東 煙臺 264002)

    0 引言

    噴氣燃料主要用作噴氣式發(fā)動機燃料,廣泛用于軍用、民航飛機和軍事武器裝備,定期對其進行質(zhì)量指標(biāo)檢測,是噴氣式發(fā)動機正常工作的保證。由于噴氣燃料質(zhì)量指標(biāo)分析檢測方法操作大多較復(fù)雜和耗時,不適合要求及時給出質(zhì)量分析報告的工作場合,例如軍用燃料應(yīng)急分析檢測、石化煉油廠噴氣燃料生產(chǎn)在線檢驗等[1]。近紅外光譜技術(shù)(NIR)憑借綜合運用計算機技術(shù)、光譜技術(shù)和化學(xué)計量學(xué),僅通過采集樣品的一張光譜結(jié)合多個質(zhì)量指標(biāo)定量分析模型就可以同時獲得這些質(zhì)量指標(biāo)的分析數(shù)據(jù),近年來廣泛應(yīng)用于石油產(chǎn)品的分析應(yīng)用中[2-4]。

    由于近紅外光譜分析(以下簡稱NIR)是儀器本身與分析模型兩部分共同作用的一種間接校正分析,特別是分析模型依賴于標(biāo)準(zhǔn)數(shù)據(jù)和建模過程,因此對近紅外光譜分析方法誤差的評定不能單純地依靠以往對測量結(jié)果的誤差分析。測量不確定度是在誤差理論的基礎(chǔ)之上延伸出的新的表征測量結(jié)果的概念,被定義為“一個與測定方法有關(guān)的、用來表征由被測變量分散程度的參數(shù)”[5]。它表示了重復(fù)性測量條件下得到的測量結(jié)果的分散性,它的數(shù)值大小跟測量方法、數(shù)據(jù)處理方法等有關(guān)系[6]。目前,測量不確定度應(yīng)用于石油產(chǎn)品的近紅外光譜分析的研究還較少,僅限于單個質(zhì)量指標(biāo)的化學(xué)分析檢測[7-13],但薛忠、徐冰[14-17]、J.Mantanus[18-20]等已將不確定度評估應(yīng)用到在醫(yī)藥學(xué)方面的中藥近紅外光譜定量分析中。

    噴氣燃料的運動粘度大小對燃料的霧化、供油量、燃料泵潤滑等都有重要影響,是一項重要質(zhì)量指標(biāo)。同時運動粘度信息與噴氣燃料的主要成分大小相關(guān),能夠利用噴氣燃料的近紅外光譜信息獲得其運動粘度的指標(biāo)大小。因此,本文中以噴氣燃料為分析對象,對其20 ℃運動粘度指標(biāo)進行NIR分析。從整體測量不確定度的角度,充分考慮近紅外光譜分析方法中各種不確定度來源,對噴氣燃料近紅外光譜分析進行全面的不確定度評估研究,以期能夠科學(xué)的評價近紅外光譜分析方法的誤差。

    1 方法原理

    1.1 標(biāo)準(zhǔn)數(shù)據(jù)測定

    1.1.1工作原理

    依據(jù)GB/T265—1988《石油產(chǎn)品運動粘度測定法和動力粘度計算法》,測定噴氣燃料20 ℃運動粘度。工作原理為用相對法測定一定體積的液體在重力作用下流經(jīng)毛細(xì)管所需時間,以求得液體的運動粘度,計算公式見式(1)。

    νt=Cτt

    (1)

    式中,νt為溫度t℃時試樣的運動粘度,mm2/s;C為毛細(xì)管粘度計常數(shù),mm2/s2;τt為溫度t℃時試樣的平均流動時間,s。

    1.1.2儀器與試劑

    HCR1043石油產(chǎn)品運動粘度測定儀,河南海克爾儀器儀表有限公司;0~100 ℃、分格為1 ℃,最大允許誤差±0.1 ℃的玻璃水銀溫度計;精確到0.01 s的電子秒表;內(nèi)徑為0.4 mm、常數(shù)C=0.004 113 mm2/s2的毛細(xì)血管粘度計;石油醚;3號噴氣燃料(68個樣品)。

    1.2 近紅外光譜分析模型建立

    1.2.1分析儀器

    NIR-2000近紅外光譜儀配有NIR分析軟件,北京英賢儀器有限公司;5 cm玻璃樣品池,CCD/2048像素檢測器,光譜范圍700~1 100 nm,波長準(zhǔn)確度≤1.5 nm,波長重復(fù)性<0.05 nm(10次連續(xù)掃描),吸光度重復(fù)性≤0.000 5 AU。

    1.2.2分析過程

    將樣品裝于5 cm的玻璃樣品池中放入近紅外光譜儀中,采用NIR透射方式采集光譜,如圖1所示。

    圖1 3號噴氣燃料近紅外光譜圖Fig.1 NIR curve of No.3 jet fuel

    將68個樣品的光譜與對應(yīng)的運動粘度測定標(biāo)準(zhǔn)數(shù)值結(jié)合建立樣品集。利用K-S分類方法將樣品集分為校正集和驗證集,其中,校正集58個樣品,用于校正分析模型的建立和模型的完全交互驗證;驗證集10個樣品,不參與建模,用于模型的外部驗證。對校正集樣本進行光譜預(yù)處理,然后確定參與建模的光譜波段,利用多元校正方法建立分析模型,通過比較模型評價4個判據(jù)的數(shù)值(校正相關(guān)系數(shù)、校正標(biāo)準(zhǔn)偏差SEC、預(yù)測相關(guān)系數(shù)和預(yù)測標(biāo)準(zhǔn)偏差SEP)來評價模型性能,最終通過驗證集的預(yù)測標(biāo)準(zhǔn)偏差(SEP)值和配對t檢驗對校正模型進行驗證,從而建立穩(wěn)健準(zhǔn)確的近紅外分析模型。

    2 NIR定量分析模型不確定度評定的技術(shù)原理

    對于近紅外光譜定量分析模型的建立,由于測定儀器和測定原理復(fù)雜,不確定度來源較多,因此,為了更簡單準(zhǔn)確地進行定量分析模型的不確定度評定,采用基于驗證數(shù)據(jù)的不確定度評定法——β-容度容忍區(qū)間(β-content tolerance interval,β-CTI)。β-CTI是在β-ETI(β-期望容忍區(qū)間)的基礎(chǔ)上改進而來的,它是一種將樣品的指標(biāo)水平與不確定度聯(lián)系起來,以便進行方法驗證和不確定度計算的方法。2004年,Feinberg等[21]將不確定度與物質(zhì)的指標(biāo)水平聯(lián)系起來,率先在驗證數(shù)據(jù)的基礎(chǔ)上,采用如lnU=a+blnY等常用數(shù)學(xué)模型擬合,將物質(zhì)的指標(biāo)水平和不確定度聯(lián)系起來,不僅便于不確定度的計算,更直觀地反映出模型在實際應(yīng)用中的預(yù)測能力。β-CTI評定法通過合理的實驗設(shè)計獲取數(shù)據(jù),經(jīng)過數(shù)據(jù)處理和不同數(shù)學(xué)模型對指標(biāo)水平與不確定度進行曲線擬合,選取最佳擬合曲線就可用以預(yù)測不同指標(biāo)水平對應(yīng)的不確定度。當(dāng)指標(biāo)水平在指標(biāo)要求水平范圍內(nèi)(即設(shè)定容忍區(qū)間內(nèi))時,所得最佳擬合曲線可始終用于預(yù)測該指標(biāo)不同水平下的不確定度,免除了許多重復(fù)實驗,可操作性強,在單個實驗室就可以完成,是一個科學(xué)、實用、便利的不確定度評定方法。

    根據(jù)ISO/DTS 21748指南[22],在某一指標(biāo)水平C下,樣品測量值Y的不確定度評估模型為式(2):

    (2)

    Feinberg在忽略其他因素造成的不確定度的基礎(chǔ)上,建立了容忍區(qū)間[L,U]和不確定度u(Y)之間的關(guān)系,關(guān)系如式(3)所示:

    (3)

    式中,L為容忍區(qū)間下限;U為容忍區(qū)間上限;Qt為自由度為ν時的t分布分位數(shù)。

    由上述關(guān)系可計算定量分析不確定度u(Y),見式(4):

    (4)

    (5)

    (6)

    3 NIR分析不確定度分量來源分析

    不確定度來源分析是不確定度評定的第一步,是合理評定測量不確定度的先決條件。分析噴氣燃料運動粘度測定產(chǎn)生不確定度的主要因素并制作相關(guān)因果圖,見圖2。

    圖2 噴氣燃料運動粘度NIR分析的不確定度來源圖Fig.2 Uncertainty source chart of kinetic viscosity NIR measurement of jet fuel

    根據(jù)近紅外光譜分析的整個過程,其測量不確定度來源主要有以下幾個方面:

    1) 近紅外光譜儀引入。近紅外光譜儀是近紅外光譜分析用到的主要儀器,儀器本身是經(jīng)過校驗檢定的,并附有主要技術(shù)指標(biāo)規(guī)格和儀器的驗收報告。因此,近紅外光譜儀引入的不確定度屬于測量儀器引入的不確定度;

    2) 標(biāo)準(zhǔn)數(shù)據(jù)測定引入。噴氣燃料的質(zhì)量指標(biāo)測定結(jié)果是近紅外光譜定量分析建模的標(biāo)準(zhǔn)數(shù)據(jù),其測定結(jié)果本身具有不確定度,因此將此不確定度分量引入到近紅外光譜分析的不確定度中,屬于測量方法引入的不確定度;

    3) NIR定量分析模型建立引入。NIR定量分析模型建立是近紅外光譜分析最重要的環(huán)節(jié),是根據(jù)樣品的質(zhì)量指標(biāo)數(shù)據(jù)和光譜數(shù)據(jù)的相關(guān)性建立數(shù)學(xué)模型的過程,實際上就是一個近似的擬合過程,存在不可避免的不確定度,對預(yù)測結(jié)果的準(zhǔn)確性和可靠性都有影響。因此,NIR定量分析模型建立引入的不確定度屬于測量方法引入的不確定度。

    4 NIR不確定度評定

    不確定度通常分為兩類:A類不確定度和B類不確定度,其中A類指可以根據(jù)一系列測量值的統(tǒng)計分布進行評定的分量;B類指可以根據(jù)經(jīng)驗或其他信息獲得的概率密度函數(shù)進行評定的分量。

    4.1 近紅外光譜儀引入的不確定度分量

    由光譜儀使用說明書和實驗條件可知,儀器標(biāo)定的標(biāo)準(zhǔn)波長為807.5 nm。近紅外光譜儀引入的不確定度分量主要包括近紅外光譜儀波長準(zhǔn)確度和波長重復(fù)性引入的不確定度,二者相互獨立。近紅外光譜儀的波長準(zhǔn)確度不大于1.5 nm,即波長測量允許最大誤差值a=1.5 nm。近紅外光譜儀的波長重復(fù)性不大于0.05 nm,即對標(biāo)準(zhǔn)波長為807.5 nm 的波進行10次測量,測得波長誤差最大平均值為 0.05 nm。因此,近紅外光譜儀引入的不確定度分量計算結(jié)果見表1。

    表1 近紅外光譜儀引入的不確定度分量Table 1 Effect of the influential factors and their uncertainty contribution of near-infrared spectrograph

    4.2 標(biāo)準(zhǔn)數(shù)據(jù)測定引入的不確定度分量

    由標(biāo)準(zhǔn)數(shù)據(jù)測定實際操作過程分析可知,運動粘度標(biāo)準(zhǔn)數(shù)據(jù)測定引入的不確定度分量包括重復(fù)性實驗、恒溫浴溫度波動、時間測量、毛細(xì)管粘度計、毛細(xì)管粘度計放置位置和溫度計引入的不確定度,上述不確定度分量互相獨立。

    在重復(fù)性條件下,使用粘度計常數(shù)C為0.004 113 mm2/s2的粘度計對同一噴氣燃料試樣重復(fù)測量10次,測量結(jié)果見表2,從而可以計算得到重復(fù)性實驗帶來的A類不確定度。

    表2 噴氣燃料樣品20 ℃時運動粘度的重復(fù)測定結(jié)果Table 2 Repeated determination results of kinetic viscosity of a jet fuel sample

    由于需要測定20 ℃下的運動粘度,選取19~21 ℃溫度范圍作為實驗區(qū)間,以0.2 ℃為溫度間隔,分別測定每個溫度下的運動粘度,且每個溫度下測定3次,取其算術(shù)平均值作為每個溫度的運動粘度,并進行線性直線擬合,如圖3所示,可得運動粘度溫度變化率即該不確定度靈敏系數(shù)為0.004 068 mm2/(s·℃)。

    圖3 運動粘度與恒溫浴溫度的擬合直線Fig.3 Fitting curve of kinetic viscosity and thermostatic bath temperature

    根據(jù)相關(guān)JJG檢定規(guī)程,內(nèi)徑為0.4 mm、常數(shù)C=0.004 113 mm2/s2的毛細(xì)血管粘度計的粘度計常數(shù)的相對擴展不確定度U=0.6%(包含因子k=2)[23]。實驗選用的電子秒表規(guī)定其在10 min內(nèi)的最大允許誤差為0.07 s[24]。在相關(guān)技術(shù)資料中可以查到,平氏粘度計的傾斜誤差不超過±0.3%,半寬為0.3%[25]。0~100 ℃、分格為0.1 ℃的玻璃水銀溫度計的相對擴展不確定度U=0.05 ℃[26]。運動粘度標(biāo)準(zhǔn)數(shù)據(jù)測定引入的不確定度分量見表3。

    表3 運動粘度標(biāo)準(zhǔn)數(shù)據(jù)測定引入的不確定度分量Table 3 Effect of the influential factors and their uncertainty contribution of standard data of kinetic viscosity

    4.3 NIR定量分析模型引入的不確定度分量

    選取“3×3×6”的實驗方案,即將實驗分3個間隔時間相同的時間段(季度)進行,每隔相同時間段(季度)進行1次實驗,每次實驗包含6個運動粘度水平,每個運動粘度水平進行3次平行重復(fù)實驗,共獲得54個外部驗證集,用于定量分析的不確定度(B類)評定。

    首先,根據(jù)所建立的近紅外定量分析模型對外部驗證集樣本進行性質(zhì)預(yù)測,根據(jù)式(7)計算出每個驗證樣本的相對預(yù)測偏差:

    (7)

    為了便于在同一度量下,對不同運動粘度水平下的不確定度進行比較,采用相對預(yù)測偏差進行計算。由基于驗證數(shù)據(jù)的不確定度計算原理,采用β-CTI方法,根據(jù)式(2)—式(6),可依次計算出不同運動粘度水平下的近紅外光譜定量分析的β-CTI容忍區(qū)間和相對不確定度,如表4所示。

    表4 不同運動粘度水平的相對不確定度Table 4 Relative extension uncertainty of different kinetic viscosity levels

    為預(yù)測不同運動粘度水平下的不確定度,將相對不確定度與運動粘度水平聯(lián)系起來,通過曲線擬合,建立兩者之間的對應(yīng)關(guān)系,擬合曲線如圖4所示,可用于一定范圍內(nèi)任意運動粘度水平的不確定度預(yù)測。

    圖4 運動粘度與相對不確定度的最佳擬合曲線Fig.4 Best fitting curve between kinetic viscosity and relative uncertainty

    因此,由擬合曲線可知,當(dāng)運動粘度為1.645 mm2/s時,u(M)為0.683%。

    5 NIR分析不確定度合成

    各不確定度分量的影響權(quán)重見圖5。

    圖5 各不確定度分量對合成不確定度的影響權(quán)重Fig.5 Influence weight of the different components on the synthetic uncertainty

    取置信水平為95%,包含因子k=2,得到相對擴展不確定度U=urel×2 = 0.78%×2=1.56%,則該噴氣燃料樣品20 ℃運動粘度近紅外光譜分析結(jié)果表示為1.645±0.026 mm2/s。

    6 結(jié)論

    本文中完成了噴氣燃料近紅外光譜分析整體不確定度的評定過程。其中,近紅外光譜儀、運動粘度標(biāo)準(zhǔn)數(shù)據(jù)測定、NIR定量分析模型等對試驗結(jié)果都有一定影響,得到的近紅外光譜儀、標(biāo)準(zhǔn)數(shù)據(jù)測定和NIR定量分析模型引入的不確定度分量分別為0.11%、0.36%和0.683%,各自在整體不確定度中所占權(quán)重分別為10%、31%和59%。最后得到噴氣燃料近紅外光譜分析的整體相對合成不確定度為1.56%。從不確定度分量的大小可以看出,NIR定量分析方法的準(zhǔn)確性在很大程度上依賴于所建分析模型。因此,在分析過程中應(yīng)嚴(yán)格按照標(biāo)準(zhǔn)進行試驗,并重點進行分析模型的構(gòu)建和優(yōu)化。噴氣燃料的近紅外光譜分析整體不確定度的評定可有效評價該類分析方法的準(zhǔn)確性和可靠性,并可為其他物質(zhì)近紅外光譜分析測定的不確定度評估提供借鑒。

    猜你喜歡
    紅外光譜儀噴氣分量
    帽子的分量
    一物千斤
    智族GQ(2019年9期)2019-10-28 08:16:21
    MEMS近紅外光譜儀微鏡驅(qū)動控制系統(tǒng)研發(fā)
    論《哈姆雷特》中良心的分量
    分量
    顯微紅外光譜儀不同檢測模式的探索
    噴氣的板車
    小布老虎(2016年10期)2016-12-01 05:46:41
    噴氣的木板車
    The United States Set Out to Build a New Generation of Unmanned Drones SR-71
    科技傳播(2013年22期)2013-10-17 11:16:36
    噴氣滑板
    亚洲精品影视一区二区三区av| 国产 一区精品| 欧美日韩综合久久久久久| 欧美日韩一区二区视频在线观看视频在线 | 亚洲aⅴ乱码一区二区在线播放| 亚洲最大成人av| 淫秽高清视频在线观看| 日本猛色少妇xxxxx猛交久久| 国产一区二区亚洲精品在线观看| 久久久久久久午夜电影| 精品久久国产蜜桃| 成人特级av手机在线观看| 亚洲经典国产精华液单| 九九爱精品视频在线观看| 91精品国产九色| 精品久久久久久成人av| 久久久久久久久久久丰满| 国产乱人偷精品视频| 午夜久久久久精精品| 99热这里只有是精品50| 麻豆一二三区av精品| 精品国内亚洲2022精品成人| 国产激情偷乱视频一区二区| 亚洲欧美成人精品一区二区| 在线a可以看的网站| 亚洲在久久综合| 久久精品熟女亚洲av麻豆精品 | 亚洲av成人av| 亚洲精品色激情综合| 日韩三级伦理在线观看| 99在线人妻在线中文字幕| 成人特级av手机在线观看| 听说在线观看完整版免费高清| 校园人妻丝袜中文字幕| 免费搜索国产男女视频| 能在线免费看毛片的网站| 麻豆成人午夜福利视频| 久久精品久久久久久久性| 亚洲成av人片在线播放无| 91午夜精品亚洲一区二区三区| 国产成人福利小说| 日韩一区二区三区影片| 国产精品不卡视频一区二区| 深夜a级毛片| 久久综合国产亚洲精品| 国产91av在线免费观看| 三级男女做爰猛烈吃奶摸视频| 久久精品夜夜夜夜夜久久蜜豆| 久久精品久久精品一区二区三区| 亚洲精品一区蜜桃| 日本黄色视频三级网站网址| 日韩大片免费观看网站 | 久久久久久久久久成人| 看非洲黑人一级黄片| 一夜夜www| 婷婷色av中文字幕| 亚洲色图av天堂| 久久精品夜夜夜夜夜久久蜜豆| 免费搜索国产男女视频| 男女那种视频在线观看| 在线观看av片永久免费下载| 国产精品日韩av在线免费观看| 乱系列少妇在线播放| 免费人成在线观看视频色| 人妻夜夜爽99麻豆av| 亚洲国产欧美在线一区| 热99re8久久精品国产| 精品久久久久久久人妻蜜臀av| 日韩av不卡免费在线播放| 久久久久久国产a免费观看| 国产人妻一区二区三区在| 欧美潮喷喷水| 自拍偷自拍亚洲精品老妇| 欧美性感艳星| 国产精品久久久久久久久免| 三级男女做爰猛烈吃奶摸视频| 精品久久久久久电影网 | 亚洲精品aⅴ在线观看| 久久人人爽人人爽人人片va| 日韩一区二区三区影片| 免费无遮挡裸体视频| 亚洲欧美成人综合另类久久久 | 午夜福利网站1000一区二区三区| 中文字幕熟女人妻在线| 亚洲国产成人一精品久久久| 日本黄色视频三级网站网址| 亚洲国产欧美在线一区| 国产亚洲一区二区精品| 久久久久久久久久久免费av| 亚洲av成人av| 国产精品伦人一区二区| 国内精品一区二区在线观看| 最近最新中文字幕免费大全7| 国产精华一区二区三区| 国产探花在线观看一区二区| 人人妻人人澡欧美一区二区| 欧美成人a在线观看| 又爽又黄a免费视频| 中文亚洲av片在线观看爽| 日韩三级伦理在线观看| 国产精品国产三级国产专区5o | 亚洲欧美清纯卡通| 精品久久久噜噜| ponron亚洲| 日韩强制内射视频| 成年av动漫网址| 国产高清不卡午夜福利| 麻豆国产97在线/欧美| 国产美女午夜福利| 少妇熟女aⅴ在线视频| 最新中文字幕久久久久| 麻豆成人午夜福利视频| 伊人久久精品亚洲午夜| 国产精品久久久久久精品电影| 国产亚洲av嫩草精品影院| 日本一本二区三区精品| 男的添女的下面高潮视频| 丰满乱子伦码专区| 别揉我奶头 嗯啊视频| 你懂的网址亚洲精品在线观看 | 伊人久久精品亚洲午夜| 色综合亚洲欧美另类图片| 国产亚洲精品久久久com| 久久精品熟女亚洲av麻豆精品 | 国内揄拍国产精品人妻在线| 成人av在线播放网站| 欧美日本视频| 国产成人免费观看mmmm| 少妇熟女欧美另类| av天堂中文字幕网| 久久久精品欧美日韩精品| 国产单亲对白刺激| 高清午夜精品一区二区三区| 我要看日韩黄色一级片| 一级av片app| 尾随美女入室| 亚洲婷婷狠狠爱综合网| 波野结衣二区三区在线| 可以在线观看毛片的网站| 人人妻人人看人人澡| 国产精品乱码一区二三区的特点| 亚洲色图av天堂| 午夜福利网站1000一区二区三区| 草草在线视频免费看| 综合色丁香网| 水蜜桃什么品种好| 丝袜喷水一区| 国产白丝娇喘喷水9色精品| 国产精品福利在线免费观看| 国产一区二区亚洲精品在线观看| 午夜福利视频1000在线观看| 国产v大片淫在线免费观看| 欧美日本视频| 亚洲激情五月婷婷啪啪| 国产亚洲91精品色在线| 精品久久久久久久末码| 99久久精品一区二区三区| 国产精品福利在线免费观看| 欧美成人一区二区免费高清观看| 中文资源天堂在线| av又黄又爽大尺度在线免费看 | 精品久久久久久久人妻蜜臀av| 亚洲精品日韩av片在线观看| 国产片特级美女逼逼视频| 久久久久免费精品人妻一区二区| 99热网站在线观看| 国产精品国产三级国产专区5o | 亚洲欧美中文字幕日韩二区| 国产精品.久久久| 国产亚洲91精品色在线| 九九爱精品视频在线观看| 1000部很黄的大片| 卡戴珊不雅视频在线播放| 国产成人精品一,二区| 午夜日本视频在线| 色播亚洲综合网| 亚洲三级黄色毛片| 亚洲av成人精品一二三区| 在线免费十八禁| 伊人久久精品亚洲午夜| 男女边吃奶边做爰视频| 精品人妻一区二区三区麻豆| 男女啪啪激烈高潮av片| 午夜老司机福利剧场| 午夜视频国产福利| 免费电影在线观看免费观看| 91精品国产九色| 日本-黄色视频高清免费观看| 噜噜噜噜噜久久久久久91| 精品久久久久久成人av| 日本五十路高清| 91精品伊人久久大香线蕉| 欧美97在线视频| 精品人妻视频免费看| 亚洲综合精品二区| 亚洲国产精品合色在线| 激情 狠狠 欧美| 久久精品国产亚洲av天美| 久久久久久九九精品二区国产| 搡老妇女老女人老熟妇| 亚洲av中文字字幕乱码综合| 亚洲精品乱久久久久久| 成年av动漫网址| 男的添女的下面高潮视频| 天美传媒精品一区二区| 天堂网av新在线| 91aial.com中文字幕在线观看| 插阴视频在线观看视频| 免费看光身美女| 国产亚洲91精品色在线| 高清在线视频一区二区三区 | 2021少妇久久久久久久久久久| 国产伦理片在线播放av一区| 国产成人精品婷婷| 在线播放国产精品三级| 国产 一区精品| 老司机影院成人| 国产免费男女视频| 久久精品影院6| 午夜福利成人在线免费观看| 日日啪夜夜撸| 欧美xxxx黑人xx丫x性爽| 免费av观看视频| 国产私拍福利视频在线观看| 精品99又大又爽又粗少妇毛片| 边亲边吃奶的免费视频| 成人无遮挡网站| 国产精品人妻久久久久久| 精品一区二区免费观看| 久热久热在线精品观看| av在线观看视频网站免费| 一级毛片我不卡| 久久亚洲国产成人精品v| 国产私拍福利视频在线观看| 老司机福利观看| 天天躁夜夜躁狠狠久久av| 国产成人a∨麻豆精品| 99久久精品一区二区三区| 成人国产麻豆网| 狂野欧美白嫩少妇大欣赏| 大香蕉久久网| 丰满乱子伦码专区| 黄色欧美视频在线观看| 我的女老师完整版在线观看| 欧美极品一区二区三区四区| 久久99蜜桃精品久久| 亚洲欧美清纯卡通| 久久久久久久久久久免费av| 日本av手机在线免费观看| 欧美97在线视频| 亚洲av成人精品一二三区| 国产精品国产三级国产专区5o | 秋霞在线观看毛片| 纵有疾风起免费观看全集完整版 | av女优亚洲男人天堂| 欧美不卡视频在线免费观看| 日本av手机在线免费观看| 国产欧美日韩精品一区二区| 亚州av有码| 亚洲人成网站在线播| 欧美又色又爽又黄视频| 国产美女午夜福利| 午夜福利网站1000一区二区三区| 禁无遮挡网站| 国产熟女欧美一区二区| 欧美xxxx黑人xx丫x性爽| 国产伦一二天堂av在线观看| or卡值多少钱| 成人欧美大片| 国产国拍精品亚洲av在线观看| 蜜桃久久精品国产亚洲av| 免费看a级黄色片| 日本欧美国产在线视频| 人体艺术视频欧美日本| 男女视频在线观看网站免费| 久久鲁丝午夜福利片| 亚洲国产欧美人成| 青春草视频在线免费观看| 亚洲国产精品专区欧美| 欧美高清性xxxxhd video| 国产精品福利在线免费观看| 久久久久性生活片| 亚洲精品亚洲一区二区| 亚洲欧美一区二区三区国产| 麻豆久久精品国产亚洲av| 欧美成人a在线观看| 亚洲,欧美,日韩| 成年女人永久免费观看视频| 国产淫语在线视频| 久久人人爽人人爽人人片va| 亚洲在线观看片| 国内揄拍国产精品人妻在线| 亚洲av日韩在线播放| 一区二区三区免费毛片| 女人久久www免费人成看片 | 亚洲av日韩在线播放| 久久午夜福利片| 欧美+日韩+精品| 国产成人免费观看mmmm| 午夜爱爱视频在线播放| 国产淫语在线视频| 成人三级黄色视频| 亚洲av电影在线观看一区二区三区 | av线在线观看网站| 欧美日韩综合久久久久久| 国产成人福利小说| 免费看av在线观看网站| 免费无遮挡裸体视频| 国产欧美另类精品又又久久亚洲欧美| 51国产日韩欧美| 两性午夜刺激爽爽歪歪视频在线观看| 精品一区二区免费观看| 国产成人aa在线观看| 熟妇人妻久久中文字幕3abv| 国产精品久久久久久精品电影| 精品一区二区三区视频在线| 欧美日本视频| 国产精品av视频在线免费观看| 国产人妻一区二区三区在| 小说图片视频综合网站| 能在线免费看毛片的网站| 亚洲色图av天堂| 国产黄片视频在线免费观看| 免费看光身美女| 偷拍熟女少妇极品色| 亚洲三级黄色毛片| 精品久久久久久久久av| av专区在线播放| 国产欧美另类精品又又久久亚洲欧美| 亚洲在线自拍视频| 黑人高潮一二区| 又爽又黄a免费视频| 热99在线观看视频| 国产av一区在线观看免费| 成人国产麻豆网| 欧美成人免费av一区二区三区| 2021天堂中文幕一二区在线观| 午夜激情福利司机影院| 日产精品乱码卡一卡2卡三| 欧美xxxx性猛交bbbb| 亚洲av日韩在线播放| 国产亚洲av嫩草精品影院| 国产av码专区亚洲av| 18禁在线无遮挡免费观看视频| 少妇高潮的动态图| 好男人视频免费观看在线| 日本熟妇午夜| 日日摸夜夜添夜夜爱| 亚洲国产日韩欧美精品在线观看| 免费大片18禁| 久久鲁丝午夜福利片| 日本爱情动作片www.在线观看| 国产伦在线观看视频一区| 中国国产av一级| 18+在线观看网站| 久久久久久大精品| 91狼人影院| 嫩草影院精品99| 日韩欧美在线乱码| 深爱激情五月婷婷| 不卡视频在线观看欧美| 久久99热这里只频精品6学生 | 床上黄色一级片| 成人一区二区视频在线观看| 男女国产视频网站| 亚洲性久久影院| 视频中文字幕在线观看| 亚洲欧美日韩东京热| 日本爱情动作片www.在线观看| 久久久久久久久久成人| 久久6这里有精品| 国产精品乱码一区二三区的特点| 卡戴珊不雅视频在线播放| 成年版毛片免费区| 97热精品久久久久久| 成人亚洲欧美一区二区av| 成人无遮挡网站| 熟女电影av网| av女优亚洲男人天堂| 成人毛片a级毛片在线播放| 国产片特级美女逼逼视频| 午夜福利网站1000一区二区三区| 国产麻豆成人av免费视频| 男人舔奶头视频| 免费在线观看成人毛片| 日韩欧美精品v在线| 男人狂女人下面高潮的视频| 国内揄拍国产精品人妻在线| 久久精品国产亚洲av天美| 国产精品久久久久久精品电影| 精品人妻一区二区三区麻豆| 亚洲国产欧美在线一区| 小蜜桃在线观看免费完整版高清| 国产精品女同一区二区软件| 九草在线视频观看| 日韩高清综合在线| 永久免费av网站大全| 亚洲国产精品成人久久小说| 亚洲精品日韩在线中文字幕| 女人被狂操c到高潮| 97人妻精品一区二区三区麻豆| 女的被弄到高潮叫床怎么办| 久久99热这里只有精品18| videos熟女内射| 亚洲美女视频黄频| 麻豆精品久久久久久蜜桃| 波野结衣二区三区在线| 成人特级av手机在线观看| 亚洲av不卡在线观看| 蜜臀久久99精品久久宅男| 在线免费观看不下载黄p国产| 亚洲第一区二区三区不卡| 久久综合国产亚洲精品| 老女人水多毛片| 午夜a级毛片| 亚洲国产精品久久男人天堂| 国产一区亚洲一区在线观看| 波多野结衣巨乳人妻| 99九九线精品视频在线观看视频| 久久久久国产网址| 最后的刺客免费高清国语| 亚洲中文字幕一区二区三区有码在线看| av女优亚洲男人天堂| 午夜日本视频在线| 特大巨黑吊av在线直播| 国产精品一及| 久久精品熟女亚洲av麻豆精品 | 欧美激情在线99| 成人亚洲欧美一区二区av| 免费观看的影片在线观看| 成人毛片a级毛片在线播放| 99热6这里只有精品| 亚洲国产精品sss在线观看| 精品久久久久久成人av| 国产欧美另类精品又又久久亚洲欧美| 亚洲人成网站在线观看播放| 精品免费久久久久久久清纯| .国产精品久久| 日本欧美国产在线视频| 国产精品国产三级专区第一集| 国产一级毛片在线| 国产成人免费观看mmmm| 久久这里有精品视频免费| 深爱激情五月婷婷| 国产精品一区二区在线观看99 | 国产中年淑女户外野战色| 黄片无遮挡物在线观看| 欧美zozozo另类| 国产毛片a区久久久久| 国模一区二区三区四区视频| 欧美精品一区二区大全| 国产69精品久久久久777片| 久久精品夜夜夜夜夜久久蜜豆| 最近最新中文字幕大全电影3| 有码 亚洲区| 亚洲国产日韩欧美精品在线观看| 亚洲国产精品成人综合色| 国产熟女欧美一区二区| 狂野欧美白嫩少妇大欣赏| 少妇人妻一区二区三区视频| 高清在线视频一区二区三区 | 又爽又黄a免费视频| 久久精品影院6| 国产伦理片在线播放av一区| 亚洲五月天丁香| 亚洲欧美精品自产自拍| 青春草国产在线视频| 久久精品国产自在天天线| 国产亚洲av片在线观看秒播厂 | 久久6这里有精品| 熟女电影av网| 亚洲性久久影院| av卡一久久| 国产成人精品久久久久久| 狂野欧美白嫩少妇大欣赏| 91av网一区二区| 级片在线观看| 久久久久久大精品| 亚洲色图av天堂| 日韩亚洲欧美综合| 国产一区二区三区av在线| 国产乱人偷精品视频| 国产成人福利小说| 成人高潮视频无遮挡免费网站| 国语对白做爰xxxⅹ性视频网站| ponron亚洲| 蜜桃亚洲精品一区二区三区| 爱豆传媒免费全集在线观看| 中国美白少妇内射xxxbb| 久久韩国三级中文字幕| 国产高潮美女av| 久久精品国产99精品国产亚洲性色| 亚洲伊人久久精品综合 | 日韩精品有码人妻一区| 久久久精品大字幕| 午夜亚洲福利在线播放| 成人亚洲欧美一区二区av| 国产成人aa在线观看| av免费观看日本| 国产精品麻豆人妻色哟哟久久 | 美女高潮的动态| 菩萨蛮人人尽说江南好唐韦庄 | 日韩欧美三级三区| 中国美白少妇内射xxxbb| 亚洲婷婷狠狠爱综合网| 91午夜精品亚洲一区二区三区| 亚洲国产成人一精品久久久| 天堂√8在线中文| 搞女人的毛片| 国内精品宾馆在线| 精品人妻视频免费看| 国产 一区 欧美 日韩| 国产高清国产精品国产三级 | 午夜爱爱视频在线播放| eeuss影院久久| 亚洲在线观看片| 少妇裸体淫交视频免费看高清| 人体艺术视频欧美日本| 中文字幕精品亚洲无线码一区| 久久精品久久久久久噜噜老黄 | 精品免费久久久久久久清纯| 日本-黄色视频高清免费观看| 久久99热6这里只有精品| 菩萨蛮人人尽说江南好唐韦庄 | 国产精品人妻久久久影院| 高清视频免费观看一区二区 | 国产色爽女视频免费观看| 乱人视频在线观看| 最近中文字幕2019免费版| 一级黄片播放器| 国产成人91sexporn| 国产精品人妻久久久久久| 欧美xxxx黑人xx丫x性爽| 男人舔奶头视频| 国产乱人偷精品视频| 中文字幕熟女人妻在线| 最近中文字幕2019免费版| 久久久久性生活片| 最近最新中文字幕免费大全7| 久久亚洲国产成人精品v| 日韩高清综合在线| 亚洲欧美精品专区久久| 一级爰片在线观看| 最近中文字幕高清免费大全6| 精品一区二区三区人妻视频| 亚洲av不卡在线观看| 国产黄a三级三级三级人| 岛国毛片在线播放| 日韩 亚洲 欧美在线| 九九久久精品国产亚洲av麻豆| 亚洲人成网站在线观看播放| 欧美成人一区二区免费高清观看| 久久99精品国语久久久| 国产色爽女视频免费观看| 免费观看人在逋| 国产午夜福利久久久久久| 国产一级毛片在线| 国产黄a三级三级三级人| 亚洲色图av天堂| 少妇人妻精品综合一区二区| 亚洲av男天堂| 久久久久久久久久成人| 好男人在线观看高清免费视频| 日韩中字成人| 午夜视频国产福利| 一级毛片我不卡| 欧美日韩在线观看h| 久久久精品大字幕| 一二三四中文在线观看免费高清| 久久久久久九九精品二区国产| 高清毛片免费看| or卡值多少钱| 欧美三级亚洲精品| 18禁动态无遮挡网站| 亚洲欧洲日产国产| 精品国产一区二区三区久久久樱花 | 国产在线男女| av国产久精品久网站免费入址| 2021天堂中文幕一二区在线观| 99热这里只有是精品50| 欧美激情在线99| 亚洲av中文av极速乱| 国产又黄又爽又无遮挡在线| 精品免费久久久久久久清纯| 久久久久久九九精品二区国产| av线在线观看网站| 偷拍熟女少妇极品色| 我的老师免费观看完整版| 日本三级黄在线观看| 18禁在线播放成人免费| 麻豆久久精品国产亚洲av| 高清毛片免费看| 美女xxoo啪啪120秒动态图| 大话2 男鬼变身卡| 亚洲精品乱码久久久v下载方式| 舔av片在线| 成年女人看的毛片在线观看| 女的被弄到高潮叫床怎么办| 九九热线精品视视频播放| 久久精品91蜜桃| 中文字幕制服av| 综合色av麻豆| 欧美日韩一区二区视频在线观看视频在线 | 草草在线视频免费看| 国产伦理片在线播放av一区| 1024手机看黄色片| 菩萨蛮人人尽说江南好唐韦庄 | a级毛片免费高清观看在线播放| 插逼视频在线观看| 女的被弄到高潮叫床怎么办| 国产视频首页在线观看| av.在线天堂| 又爽又黄无遮挡网站| 国产精品av视频在线免费观看| 午夜精品一区二区三区免费看|