999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

基于近紅外光譜技術(shù)的魚雷燃料組分含量分析系統(tǒng)

2012-05-28 02:17:06郁俊江
水下無人系統(tǒng)學(xué)報 2012年3期
關(guān)鍵詞:分析系統(tǒng)

李 偉, 郁俊江, 鄧 鵬

?

基于近紅外光譜技術(shù)的魚雷燃料組分含量分析系統(tǒng)

李 偉1, 郁俊江2, 鄧 鵬1

(1. 海軍潛艇學(xué)院 導(dǎo)彈兵器系, 山東 青島, 266042; 2. 海軍91663部隊, 山東 青島, 266071)

分析了現(xiàn)行魚雷燃料成分測定方法在使用中存在的問題, 基于相關(guān)的化學(xué)計量學(xué)和光譜理論研究, 提出了將近紅外光譜(NIR)技術(shù)運用到魚雷燃料組分含量檢測分析的方法。基于選擇建模樣品集、恢復(fù)近紅外光譜弱信號以及建立和優(yōu)化近紅外數(shù)學(xué)模型, 建立了軟硬件結(jié)合的魚雷動力燃料組分含量NIR分析系統(tǒng), 依據(jù)該系統(tǒng)和相關(guān)方法, 對不同批次和儲存狀態(tài)的魚雷燃料樣本進(jìn)行了取樣化驗和成分測定, 并與現(xiàn)行方法進(jìn)行了試驗數(shù)據(jù)對比分析, 驗證了本系統(tǒng)的準(zhǔn)確性和適用性, 實現(xiàn)了遠(yuǎn)程在線分析、車載移動分析和現(xiàn)場應(yīng)急分析。

魚雷燃料; 近紅外光譜; 組分含量

0 引言

為了保證魚雷的航速和航程等與作戰(zhàn)使用密切相關(guān)的動力性能, 魚雷燃料組分含量范圍要求嚴(yán)格控制, 因此需要定時進(jìn)行化學(xué)成分的分析。

美國的OTTO-Ⅱ燃料是目前世界上較為先進(jìn)和使用范圍較廣的熱動力魚雷用推進(jìn)劑[1]。該推進(jìn)劑能量密度大、毒性小、安全使用性能好。本文以與OTTO-Ⅱ化學(xué)成分相同或相近的水下動力推進(jìn)用燃料為研究對象, 采用新技術(shù)和新方法, 對其組分含量進(jìn)行分析。

1 現(xiàn)有分析方法存在的問題

常用的魚雷動力燃料主要由1,2-丙二醇二硝酸酯(PGDN)、2-硝基二苯胺 (2-NDPA)、癸二酸二丁酯(DBS)3種成分組成, 目前對這3種成分含量的測定基本沿用20世紀(jì)70年代美國軍用規(guī)范OTTOⅡ燃料MIL-O-87672(OS)的做法[2], 即利用氮量法來測試燃料中PGDN的百分含量, 利用分光光度法或容量溴化法來測試燃料中2-NDPA的含量, 燃料中DBS的含量則用減量法測定, 即由100%減去PGDN和2-NDPA的百分含量求得。PGDN測定要求2次平行測定的允許差值不大于0.3%, 2-NDPA測定要求2次平行誤差應(yīng)不大于0.05%, 兩者均取其算術(shù)平均值為測定結(jié)果。

現(xiàn)有測定方法在使用中存在的問題具體表現(xiàn)在以下幾個方面。

1) 化驗所用的汞是一種具有沉積毒性的重金屬, 可以透過皮膚侵入人體, 還容易污染環(huán)境。所用濃硫酸具有強(qiáng)腐蝕性和強(qiáng)氧化性, 化學(xué)反應(yīng)過程中生成的氧化氮及二氧化氮等氣體容易對人體呼吸道造成傷害, 對環(huán)境危害極大。

2) 對化驗場所環(huán)境要求高。當(dāng)化驗環(huán)境溫度過高時, 反應(yīng)將過快而無法控制, 易導(dǎo)致化驗失敗。同時, 化驗場所要配備通風(fēng)、防爆等設(shè)施。因此, 現(xiàn)行方法對燃料進(jìn)行化驗通常要在冬春季進(jìn)行, 不利于全年實時監(jiān)測燃料技術(shù)狀態(tài)。

3) 化驗分析時間長, 工作強(qiáng)度大, 人員操作技能要求高。在操作過程中, 必須嚴(yán)格控制所稱樣品質(zhì)量和反應(yīng)速度, 在反應(yīng)過程中要長時間搖晃裝有大量汞的氮素反應(yīng)器, 存在安全隱患。

4) 機(jī)動檢測適應(yīng)性不強(qiáng)。在試驗過程中要用到大量的汞和濃硫酸、冰醋酸、乙醇等溶劑, 且化驗器材大多為玻璃材質(zhì), 運輸中易碎, 運輸難度大, 不利于機(jī)動檢測。

以上存在的不足限制了魚雷動力燃料安全快速檢測, 需要依托新技術(shù)探索更為有效的方法和手段。

2 光譜分析理論

近紅外光譜(near infrared, NIR)[3]主要是含X-H(X=C、N、O)基團(tuán)的樣品在中紅外區(qū)域基頻振動的倍頻和合頻吸收。魚雷動力燃料中, 1,2-丙二醇二硝酸酯(PGDN)、2-硝基二苯胺 (2-NDPA)、癸二酸二丁酯 (DBS)這3種主要成分都含有C-H或N-H基團(tuán), 這就為采用NIR技術(shù)定量分析魚雷燃料中3種組分含量奠定了化學(xué)信息基礎(chǔ)。

雖然NIR理論上非常適合用于碳?xì)溆袡C(jī)物質(zhì)的組成性質(zhì)測量[4], 但是在該區(qū)域內(nèi), 含氫基團(tuán)化學(xué)鍵振動的倍頻與合頻吸收強(qiáng)度很弱, 靈敏度相對較低, 吸收帶較寬且重疊嚴(yán)重, 因此, 依靠傳統(tǒng)建立工作曲線的方法進(jìn)行定量分析是十分困難的, 這也是早期影響NIR分析技術(shù)發(fā)展的致命原因所在。而化學(xué)計量學(xué)的發(fā)展為這一問題的解決奠定了數(shù)學(xué)基礎(chǔ)。

化學(xué)計量學(xué)是綜合使用數(shù)學(xué)、統(tǒng)計學(xué)和計算機(jī)科學(xué)等方法從化學(xué)測量數(shù)據(jù)中提取信息的一門新興的交叉學(xué)科。大量化學(xué)計量學(xué)方法被編成軟件, 并成為分析儀器(尤其是NIR分析儀)的重要組成部分。其中的數(shù)學(xué)處理方法主要有[4]: 多元線性回歸(multiple linear regression, MLR)、逐步回歸(stepwise regression, SWR)、主組分分析(principal component analysis, PCA)、主組分回歸(principal component regression, PCR)、偏最小二乘法(partial least squares, PLS)、人工神經(jīng)網(wǎng)絡(luò)(artificial neural network, ANN)和拓?fù)?Toplogical)等。MLR和SMR方法在分析樣品時只用了一些特征波長點的光譜信息, 其他點的信息被丟失, 易產(chǎn)生模型的過適應(yīng)性(Over fitting)。PCR和PLS方法的顯著特點就是利用了全部光譜信息, 可以壓縮所需樣品數(shù)量, 將高度相關(guān)的波長點歸于一個獨立變量中, 根據(jù)為數(shù)不多的獨立變量建立回歸方程, 通過內(nèi)部檢驗來防止過模型現(xiàn)象, 相比MLR和SMR分析精度提高。

3 魚雷燃料的NIR分析方法

NIR分析是近十年來發(fā)展最為迅速的高新分析技術(shù)之一。NIR系統(tǒng)的一個重要特點就是技術(shù)本身的成套性, 即近紅外光譜儀、化學(xué)計量學(xué)軟件和應(yīng)用模型的三位一體性, 性能優(yōu)異的近紅外光譜儀是該技術(shù)的基礎(chǔ)和前提。在軍工應(yīng)用領(lǐng)域, NIR已被用于無機(jī)物紅煙硝酸密度、四氧化二氮、磷酸、氫氟酸和水分含量等5項指標(biāo)的快速測定, 并測定了液體推進(jìn)劑燃料混胺中的三乙胺、二甲苯胺和水分含量以及密度等; 測定了發(fā)射藥生產(chǎn)過程中固體含能材料結(jié)構(gòu)內(nèi)揮發(fā)和結(jié)構(gòu)外揮發(fā)的含量[5]。

本文首次提出了將NIR分析技術(shù)運用到魚雷動力燃料組分含量檢測中, 提取燃料光譜特征, 建立并優(yōu)化數(shù)學(xué)模型, 研制軟硬件結(jié)合的NIR分析系統(tǒng), 可以大大提高部隊魚雷燃料的檢測效率, 減輕勞動強(qiáng)度, 滿足使用要求。圖1是NIR分析系統(tǒng)的工作流程圖。

梨友開始還勸秀容月明為他們效力,苦口婆心,不厭其煩,見秀容月明絲毫不為所動,就叫大將馬刺把他押往北國,讓單于處置,自己繼續(xù)率兵向南挺進(jìn)。

圖1 NIR分析系統(tǒng)工作流程圖

3.1 建模樣品集的選擇

數(shù)學(xué)模型預(yù)測樣品的效果決定于建模所用數(shù)據(jù)以及(用算法)對建模數(shù)據(jù)中信息的充分提取。NIR分析大致有一半誤差來自于建模數(shù)據(jù), 因此, 運用關(guān)聯(lián)算法建立數(shù)學(xué)模型前, 需要運用恢復(fù)、壓縮技術(shù)對數(shù)據(jù)進(jìn)行前處理, 依據(jù)“少而精”的原則, 提高建模數(shù)據(jù)的信息有效率。

為了保證NIR分析系統(tǒng)的穩(wěn)定性, 必須全面考慮魚雷燃料樣品的復(fù)雜背景, 這就需要大量收集燃料樣品, 并在規(guī)范的條件下測定這些樣品的NIR和組分含量。根據(jù)魚雷燃料樣品的不同類型(組分含量), 運用算法從譜庫中選擇代表性樣品作為建模樣品集。對收集配制的魚雷燃料樣品, 測定其主組分含量后進(jìn)行挑選。在挑選過程中盡可能包含那些處于兩端(即最高、最低范圍內(nèi))的樣品。由于大多數(shù)自然產(chǎn)物化學(xué)值呈正態(tài)分布, 因此選擇的魚雷燃料建模樣品集的組分含量也應(yīng)該呈正態(tài)分布, 且覆蓋范圍較寬, 從而保證組分含量NIR分析系統(tǒng)對大多數(shù)魚雷燃料樣品的適應(yīng)性。

對于大多數(shù)對象來說, 一般所挑選的建模樣品數(shù)應(yīng)不少于50個, 但也不能過多, 樣品數(shù)過多會帶來更多的干擾因素, 因為模型中每增加一個樣品, 在增加信息量的同時也增加了干擾信息, 當(dāng)建模樣品的數(shù)量太多時干擾信息會掩蓋大部分有用信息, 使所建模型性能大大降低。由此選擇了70個具有代表性的均勻且主組分含量呈自然正態(tài)分布的魚雷燃料樣品作為建模樣品集。對于異常樣品, 其中有些是光譜或組分含量誤差較大, 這些異常樣品的有效信息率低, 應(yīng)加以剔除; 有些異常樣品的光譜或組分含量是正確的, 只是由于其類型(背景)特殊, 這些異常樣品有的保留在校正樣品集中, 這樣雖然有可能降低預(yù)測的準(zhǔn)確度, 但卻可以提高所建模型的適配范圍, 因此要考慮實際情況, 根據(jù)整個NIR分析系統(tǒng)的可靠性和適配范圍來決定取舍。

3.2 NIR集弱信號的恢復(fù)

NIR有效信息的強(qiáng)度較低, 屬于弱信號; 另外, 每個光譜數(shù)據(jù)除了包含有效信息以外還包含測量誤差等不確定的、難以校正的、干擾測定的無效信息或稱干擾信息, 這些干擾信息的強(qiáng)度一般較大。弱信號的處理需要放大, 但是伴隨弱信號的放大, 噪聲也隨著放大, 所以處理弱信號首先要通過數(shù)學(xué)方法, 降低疊加在弱信號上的噪聲與系統(tǒng)干擾(誤差), 增強(qiáng)有效信息的相對強(qiáng)度, 以恢復(fù)弱信號的原有特征, 這就是弱信號的恢復(fù)。需要用校正技術(shù), 即消除或降低信號中系統(tǒng)誤差和隨機(jī)誤差的干擾, 來處理恢復(fù)NIR的弱信號。NIR的誤差主要有隨機(jī)誤差和系統(tǒng)誤差2種, 樣品集光譜的前處理就是為了降低這兩類誤差, 可以通過多次測定取平均、平滑、傅立葉濾波等方法降低隨機(jī)誤差; 用導(dǎo)數(shù)光譜等方法來降低或消除某些系統(tǒng)誤差, 可以用各種算法如各種求平均、平滑的方法。采用傅立葉技術(shù)和小波變換來提高近紅外信號的信噪比, 增強(qiáng)有效信息的相對強(qiáng)度, 從而恢復(fù)NIR集的弱信號。

3.3 NIR數(shù)學(xué)模型的建立和優(yōu)化

魚雷燃料組分含量NIR分析系統(tǒng)的建立過程主要是近紅外數(shù)學(xué)模型的建立和優(yōu)化過程, 也就是經(jīng)過前處理的魚雷燃料樣品的光譜數(shù)據(jù)與樣品主組分含量之間進(jìn)行數(shù)學(xué)關(guān)聯(lián)的過程。所謂關(guān)聯(lián)是指建立兩組量(燃料樣品NIR和與之對應(yīng)的主組分含量)的數(shù)學(xué)關(guān)系(即數(shù)學(xué)模型), 通過數(shù)學(xué)關(guān)系可以從一個量得到另外一個量。關(guān)聯(lián)的方法可以通過解方程組, 得到方程確定的解, 也可以通過統(tǒng)計擬合的方法(回歸方法)建立模型。通過特殊的算法(PCR, PLS, ANN等)來選擇適當(dāng)形式的光譜函數(shù)和定標(biāo)方程的項數(shù), 選擇最佳波長組合并確定相應(yīng)的方程常數(shù)項和各項的系數(shù), 建立多元濃度定標(biāo)方程, 即數(shù)學(xué)模型。通過對定標(biāo)方程求解和確立, 來對模型進(jìn)行創(chuàng)建和優(yōu)化。

本文采用的全譜區(qū)PLS用組分分解法, 將原光譜中上千個光譜數(shù)據(jù)點最終壓縮成幾個組分進(jìn)行選取, 再通過線性回歸或非線性回歸與待測數(shù)據(jù)關(guān)聯(lián)建立數(shù)學(xué)模型, 以提高建模光譜數(shù)據(jù)的有效信息率。

圖2為優(yōu)化處理后典型的燃料樣品NIR(吸收光譜度)圖。

圖2 魚雷燃料樣品NIR圖

3.4 分析系統(tǒng)的工程實現(xiàn)

對數(shù)學(xué)模型進(jìn)行轉(zhuǎn)移、修正, 建立魚雷燃料組分含量NIR分析系統(tǒng), 用待測樣品對NIR分析系統(tǒng)進(jìn)行評價。

該課題研究的主要成果為一套軟硬件結(jié)合的“魚雷動力燃料組分含量NIR分析系統(tǒng)”, 系統(tǒng)組成框圖如圖3所示。

4 試驗結(jié)果與分析

本系統(tǒng)研制成功后, 對魚雷燃料樣本進(jìn)行了取樣化驗和成分測定。使用中對3類樣本進(jìn)行采樣: 第1類是從雷上放出的燃料在專用容器中儲存約一周后進(jìn)行化驗(簡稱“雷上放出”); 第2類是直接對從雷上取出的燃料進(jìn)行化驗(簡稱“雷體直取”); 第3類是對還沒有裝雷的儲藏罐內(nèi)的燃料進(jìn)行化驗(簡稱“罐裝直取”)。3類情況下共進(jìn)行了98次試驗, 不同燃料批次、儲存狀態(tài)和測定時間的10次典型試驗結(jié)果匯總列表見表1。

圖3 分析系統(tǒng)組成框圖

表1 魚雷燃料成分試驗測定數(shù)據(jù)表

Table 1 Detection data of torpedo fuel contents

注: 表中單位(m/m,%)為質(zhì)量百分比。

表中的“×”、“√”分別表示魚雷燃料組分含量為不合格、合格, 合格判據(jù)為:

1) 1,2-丙二醇二硝酸酯(PGDN)含量97.8~ 98.2 , (m/m, %);

2) 2-硝基二苯胺 (2-NDPA) 含量3.6~3.8 , (m/m, %);

3) 水分含量≤2.3, (m/m, %)。

基本性能試驗中, 用本系統(tǒng)NIR模型預(yù)測魚雷燃料樣品的組分含量, PGDN預(yù)測值與GJB XX X-XX檢測方法測定值的最大相對誤差僅為研制要求的56%, 2-NDPA預(yù)測值最大相對誤差為研制要求的59%。

準(zhǔn)確度試驗中, 對于PGDN, 預(yù)測均方差(root mean square error of prediction, RMSEP)在相對誤差平均數(shù)為0.25時為0.31%, 小于0.5%的研制要求; 對于2-NDPA , RMSEP在相對誤差平均數(shù)為0.00195時為1.17%, 小于5%的研制要求。

精密度試驗中, 在30 min內(nèi)對魚雷燃料母液樣品進(jìn)行25次重復(fù)測定, 每次間隔1 min, 極差僅為0.04, RMSEP為0.05, 遠(yuǎn)小于GJB要求的平行誤差0.3。

目前, 本系統(tǒng)已經(jīng)替代原有設(shè)備, 用于裝備質(zhì)量監(jiān)測單位魚雷燃料的定期檢查。使用結(jié)果表明, 本系統(tǒng)操作中所需儀器少, 測量過程簡單, 不需要經(jīng)常標(biāo)定, 可以在復(fù)雜環(huán)境中進(jìn)行測試操作, 實現(xiàn)了遠(yuǎn)程在線分析、車載移動分析和現(xiàn)場應(yīng)急分析; 另外, 通過建立加載模型還可以實現(xiàn)一臺儀器對多種產(chǎn)品的測定, 提高了儀器的集成度和綜合使用效能, 節(jié)省了人力物力。

5 結(jié)束語

本文介紹的基于NIR分析技術(shù)的魚雷燃料組分含量分析系統(tǒng), 在高效性、快捷性、易用性、擴(kuò)展性和環(huán)保性等方面都明顯優(yōu)于現(xiàn)行檢測方法, 推動了NIR技術(shù)在武器液體推進(jìn)劑領(lǐng)域的應(yīng)用, 借助擴(kuò)展積分球等配件支持, 還可以應(yīng)用于固體推進(jìn)劑的化學(xué)分析。

[1] 于俊衛(wèi). HAP三組元推進(jìn)劑熱動力系統(tǒng)研究[J]. 魚雷技術(shù), 2001, 9(4): 33-37.

[2] 趙小鋒. 國外水下液體化學(xué)推進(jìn)劑的研究現(xiàn)狀[J]. 火炸藥學(xué)報, 2009, 32(9): 62-65.

Zhao Xiao-feng. Current Situation of Underwater Chemical Propellan[J]. Chinese Journal of Explosives and Propellants, 2009, 32(9): 62-65.

[3] 張卉, 宋研, 冷靜, 等. 近紅外光譜分析技術(shù)[J]. 光譜實驗室, 2007, 24(3): 388-395.Zhang Hui, Song Yan, Leng Jing, et al. Near Infrared Spectroscopic Analytical Technique[J]. Journal of Pharmaceutical and Biomedical Analysis, 2007, 24(5): 388-395.

[4] 解國玲, 任芊, 董守龍, 等. 近紅外光譜技術(shù)在含能材料成分分析中的建模研究[J]. 火炸藥學(xué)報, 2003, 26(11): 78-80.Xie Guo-ling, Ren Qian, Dong Shou-long, et al. Study on Modeling of Multicomponent Analysis in Propellants Products with NIR Spectroscopy[J]. Chinese Journal of Explosives & Propellants, 2003, 26(11): 78-80.

[5] 褚小立, 袁洪福, 陸婉珍. 近年來我國近紅外光譜分析技術(shù)的研究與應(yīng)用進(jìn)展[J]. 分析儀器, 2006, 37(2): 1-9. Chu Xiao-li, Yuan Hong-fu, Lu Wan-zhen. Research and Applications of Near Infrared Spectroscopy in China in Recent Years[J]. Analytical Instrumentation, 2006, 37(2): 1-9.

Detection System of Component Content for Torpedo Fuel Based on NIR

LI Wei1, JU Jun-jiang2, DENG Peng1

(1. Department of Missile and Weaponry Engineering, Navy Submarine Academy, Qingdao 266042, China; 2. 91663thUnit, The People′s Liberation Army of China, Qingdao 266071, China)

This paper analyzes the problems in detecting component content of torpedo fuel with the existing methods. Based on chemometrics and spectral theory, a method for detecting component content of torpedo fuel with near infrared (NIR) spectroscopy is proposed. By selecting modeling sample set, restoring weak NIR signal, and establishing and optimizing near infrared mathematical model, an NIR analysis system of component content for torpedo fuel is developed. With this system and related method, the torpedo fuel samples of different batches and storage states are examined and the components are detected, and the results are compared with the test data from the existing methods. The comparison verifies the accuracy and feasibility of the proposed detection system. Thus, long-range on-line analysis, mobile analysis on vehicle and emergency analysis on-the-spot are realized with this method.

torpedo fuel; near infrared (NIR) spectroscopy; component content

TJ630.34;E925.23

A

1673-1948(2012)03-0206-04

2011-12-21;

2012-02-15.

李 偉(1967-), 男, 碩士, 副教授, 主要從事魚雷作戰(zhàn)使用及技術(shù)保障教學(xué)工作.

(責(zé)任編輯: 陳 曦)

猜你喜歡
分析系統(tǒng)
Smartflower POP 一體式光伏系統(tǒng)
WJ-700無人機(jī)系統(tǒng)
隱蔽失效適航要求符合性驗證分析
ZC系列無人機(jī)遙感系統(tǒng)
北京測繪(2020年12期)2020-12-29 01:33:58
基于PowerPC+FPGA顯示系統(tǒng)
半沸制皂系統(tǒng)(下)
電力系統(tǒng)不平衡分析
電子制作(2018年18期)2018-11-14 01:48:24
連通與提升系統(tǒng)的最后一塊拼圖 Audiolab 傲立 M-DAC mini
電力系統(tǒng)及其自動化發(fā)展趨勢分析
中西醫(yī)結(jié)合治療抑郁癥100例分析
主站蜘蛛池模板: 免费一级毛片完整版在线看| 91在线精品麻豆欧美在线| 中文字幕永久视频| 激情综合网激情综合| 99在线视频免费| 中文字幕在线一区二区在线| 亚卅精品无码久久毛片乌克兰| 精品综合久久久久久97超人| 国产在线无码av完整版在线观看| 夜夜操国产| 天天摸天天操免费播放小视频| 丁香婷婷激情综合激情| 国产色网站| 国产精品视频观看裸模| 四虎亚洲精品| 免费无码AV片在线观看中文| 久久福利片| 日韩毛片视频| 亚洲三级网站| 另类重口100页在线播放| 欧美性久久久久| 麻豆精品视频在线原创| 日本一区二区不卡视频| 久久精品国产亚洲麻豆| 91破解版在线亚洲| 亚洲熟女中文字幕男人总站| 无码精品一区二区久久久| 成人精品在线观看| 亚洲精品无码高潮喷水A| 欧美精品1区2区| 欧美在线精品一区二区三区| 99福利视频导航| 香蕉国产精品视频| 美女国产在线| 国产精品免费入口视频| 亚洲欧美另类久久久精品播放的| 日本一区二区三区精品国产| 四虎永久在线精品国产免费| 成人国产精品2021| 中文纯内无码H| 国产成人艳妇AA视频在线| AV不卡无码免费一区二区三区| 欧美特黄一免在线观看| 亚洲欧美天堂网| 国产成人夜色91| 性激烈欧美三级在线播放| 国产网站一区二区三区| 欧美第一页在线| 国产精品一区二区无码免费看片| 中文成人无码国产亚洲| 欧美日韩中文字幕在线| 亚洲人在线| 免费AV在线播放观看18禁强制| 国产精品一区二区在线播放| 朝桐光一区二区| 另类欧美日韩| 国产亚洲欧美日韩在线观看一区二区| 在线观看91精品国产剧情免费| 亚洲狠狠婷婷综合久久久久| 国产91全国探花系列在线播放| 国产区免费| 69精品在线观看| 亚洲色图欧美视频| 片在线无码观看| 国内a级毛片| 日本妇乱子伦视频| 成人av手机在线观看| 一级香蕉视频在线观看| 欧美午夜在线视频| 欧美精品另类| 国产性爱网站| 日本精品一在线观看视频| 在线中文字幕网| 午夜国产精品视频| 国产女人18水真多毛片18精品| 青草视频久久| 看av免费毛片手机播放| 亚洲第一国产综合| 1769国产精品免费视频| 欧美日韩精品一区二区视频| 午夜三级在线| 久久国产热|