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

    基于多波長激光吸收光譜技術的氣體溫度二維重建

    2016-10-14 01:43:42辛明原宋俊玲
    航天器環(huán)境工程 2016年3期
    關鍵詞:模擬退火吸收光譜譜線

    辛明原,金 星,宋俊玲

    ?

    基于多波長激光吸收光譜技術的氣體溫度二維重建

    辛明原,金 星,宋俊玲

    (裝備學院 激光推進及其應用國家重點實驗室,北京 101416)

    激光吸收光譜技術在航空航天領域具有廣泛的應用前景。文章利用多波長激光吸收光譜技術設計了氣體溫度二維重建系統(tǒng),實現(xiàn)了燃燒場氣體溫度的二維重建,通過增加吸收譜線信息減少了投影個數(shù)。重建算法采用了遺傳模擬退火結合序列二次規(guī)劃(SQP)的混合算法,在實現(xiàn)全局尋優(yōu)的同時提高了運算效率。建立了產(chǎn)物包含水的燃燒模型,利用8條水的近紅外吸收譜線,對燃燒場溫度進行了二維重建。通過在測量信息中添加不同比例的隨機噪聲,驗證了算法的穩(wěn)定性,并與雙波長吸收法進行了對比。數(shù)值仿真實驗結果表明,在低噪聲環(huán)境下,兩種方法的重建結果相當;但對于噪聲較大的環(huán)境而言,多波長吸收法的重建效果更為穩(wěn)定。

    激光吸收光譜;多波長吸收;超光譜層析;遺傳模擬退火算法;二維重建

    0 引言

    基于可調(diào)諧半導體激光吸收光譜技術(TDLAS)的溫度場二維重建,通常采用雙波長吸收譜線方法[1],簡稱雙線法。但是,雙線法二維

    重建對光線投影的角度與數(shù)目有較多的要求[2]:為得到較為理想的結果,至少需要4個投影方向。浙江大學的王飛等人[3]利用4個旋轉平臺及柱面反射器實現(xiàn)了對被測區(qū)域的快速掃描,獲得了足夠的

    譜線信息。然而,航天器內(nèi)部的測量空間有限,很難完成該系統(tǒng)的搭建,而采用固定光路的辦法通常只能獲得2個投影角度。故探索在較少投影角度下的重建辦法有較強的實際意義。

    Ma L等人[4]提出的超光譜層析法(Hyperspectral Tomography, HT),也稱多波長激光吸收光譜法,利用多波長吸收信息減少了投影數(shù)量(共采用10個波長,2個投影角度)。HT法的缺點在于重建時間較長,一次重建需要消耗數(shù)個小時。Cai W等人[5]引入正交分解法縮短了重建時間,2015年又進行了高分辨率重建研究[6]。

    國內(nèi),基于TDLAS技術的溫度場二維重建研究主要集中在雙線法[7-8],對HT法的研究較少。2010年李寧等人[9]利用波分復用技術進行了多波長吸收光譜(4個波長,4個投影方向)二維重建數(shù)值仿真研究,重建結果與給定模型符合較好;但是在光線布局上沒有達到2個投影角度的目標,且溫度重建結果受噪聲影響較大。

    本文設計了基于多波長激光吸收光譜技術的氣體溫度重建系統(tǒng),模擬了產(chǎn)物包含水汽的燃燒場,采用8個吸收波長、2個投影角度進行重建。求解算法采用遺傳模擬退火結合序列二次規(guī)劃(Sequential Quadratic Programming, SQP)的復合式算法,旨在提高計算效率;目標函數(shù)引入了平滑正則化函數(shù),旨在提高算法的抗噪聲能力。

    1 重建理論

    1.1 激光吸收光譜

    根據(jù)Beer-Lambert定律,激光穿過均勻氣體介質(zhì)時透射光和入射光滿足關系式

    其中:入射光的強度為0;透射光強度為t;為激光頻率;為壓強;為吸收氣體組分濃度;()為線強度;φ為線型函數(shù);為光穿過氣體介質(zhì)的長度。

    線強度為

    其中:()為配分函數(shù);0=296K;(0)和(0)分別是參考溫度下的線強度和配分函數(shù);為普朗克常數(shù);為光速;為玻耳茲曼常數(shù);為下態(tài)能級;0為吸收線的中心頻率。()可用多項式進行近似代替[10],、(0)及(0)可在HITRAN光譜數(shù)據(jù)庫[11]中查到。

    線型函數(shù)具有歸一化的性質(zhì),對式(1)在頻域上進行積分可得

    其中為積分吸光度,也稱為光線的投影,可通過實驗測得。

    1.2 重建的數(shù)學模型

    將待測氣體場離散為×的網(wǎng)格,假定在每個網(wǎng)格內(nèi)部氣體參數(shù)均勻分布,壓強在整個待測氣體場范圍內(nèi)均勻分布。為方便表示,按圖1所示對網(wǎng)格進行編號:(1), (2),×××, (×),其中是網(wǎng)格的行數(shù)、是網(wǎng)格的列數(shù)。

    圖1 網(wǎng)格編號方式及離散化示意圖

    Tχ分別表示()編號網(wǎng)格中的溫度與組分濃度。激光頻率用ν表示(其中=1, 2,×××,,為總頻率個數(shù));光線的位置用p來表示(其中=1, 2,×××,,為總光線個數(shù));(p,)表示位置在p的光線穿過編號為()的網(wǎng)格的長度(圖1中的藍色線段);(p,ν)表示對位置為p、頻率為ν的光線的積分吸光度。基于式(3),(p,ν)可表示為

    對氣體的溫度重建問題等價于求解式(4)的方程組。但是該方程組因非線性因子()而具有高

    度的非線性特征,且未知數(shù)個數(shù)較多,直接求解的難度很大[4]。可將方程組求解問題轉換為最小化問題,其優(yōu)化目標函數(shù)為

    式中:(p,ν)m和(p,ν)c分別表示測量與重建的積分吸光度;為測量積分吸光度與重建積分吸光度的歸一化偏差總和。

    2 重建設計

    2.1 搜索算法設計

    針對式(5)的最小化問題,采用遺傳模擬退火算法結合SQP的復合式算法進行求解。算法的設計思路是:通過遺傳模擬退火算法找出全局最優(yōu)解所在區(qū)域,之后利用SQP算法進行深度搜索進而得出全局最優(yōu)解。算法程序流程如圖2所示。

    圖2 程序流程圖

    算法的具體設計如下:

    1)算法及參數(shù)的初始化,其具體參數(shù)如表1所示;

    2)計算種群中每個個體的適應度;

    3)進化次數(shù)加1;

    4)遺傳算法操作,包括選擇、交叉、變異;

    5)判斷進化次數(shù)能否被10整除,若不能則跳轉到第8步;

    6) 對種群進行模擬退火操作,降低退火溫度;

    7)重復步驟6直至退火溫度低于設定值,恢復退火溫度并退出模擬退火程序;

    8)計算群體中每個個體的適應度;

    9)保優(yōu)操作,將上一代種群中適應度最高的個體直接存入下一代種群;

    10)判斷是否滿足終止條件,即適應度變化是否小于設定值,或進化次數(shù)是否達到最大進化次數(shù)設定,如不滿足終止條件,則跳轉至第3步;

    11)以遺傳模擬退火算法得出的最優(yōu)個體作為初值進行SQP搜索,反復迭代10次得到最終重建結果。

    表1 初始化參數(shù)

    2.2 譜線選擇

    重建的氣體場離散為8×8的網(wǎng)格,光線共采用16條,2個投影角度,水平方向和豎直方向上各8條,如圖3所示。

    圖3 氣體場離散及光線布局

    H2O的吸收譜線共選擇了8條,其篩選原則可總結為:選取不同下態(tài)能級、吸收較強且與相鄰譜線干擾較小的譜線[12]。具體的譜線參數(shù)如表2所示。

    表2中最后一行為雙中心波長的重疊譜線,實驗時采用線強度為1.12×10-3 cm-2?atm-1,其中atm指單位大氣壓,下態(tài)能級為1790cm-1。

    2.3 平滑正則化函數(shù)

    式(5)的最佳逼近問題具有大量的局部最優(yōu)解,為改良問題的病態(tài)特性使重建結果穩(wěn)定,在式(5)的基礎上增加了平滑正則化函數(shù)(前提是假設氣體溫度分布是平滑分布)進行修正,即為

    式中T, n表示第行、第列網(wǎng)格中的溫度(或為0時T, n取0)。改進后的目標函數(shù)為

    式中T為正則化參數(shù),用于調(diào)整正則化函數(shù)的幅值,本文中采用T=1×10-7。

    3 數(shù)值仿真實驗

    3.1 溫度模型建立

    溫度分布采用雙高斯峰疊加拋物面的分布,模擬實際燃燒中非對稱、多峰值的溫度分布情況。溫度范圍設定為600~1400K的燃燒溫度,如圖4所示。

    圖4 溫度分布模型

    3.2 重建結果與討論

    目標函數(shù)隨迭代次數(shù)的變化如圖5所示。

    圖5 目標函數(shù)隨計算迭代次數(shù)變化曲線

    during the algorithm

    圖5中虛線的左半部分為遺傳模擬退火(GSA)運算過程,右半部分為SQP運算的過程。GSA計算過程中,遺傳算法共執(zhí)行50次,每10次執(zhí)行

    1次模擬退火運算,因此第10次迭代的目標函數(shù)下降幅度較大。算法執(zhí)行20次后,的下降趨勢變緩,說明種群趨于集中,算法已找到最優(yōu)解所在區(qū)域。在執(zhí)行50次之后采用SQP算法進行局部最優(yōu)解搜索,其搜索效率明顯提高,迅速下降到1×10-3的數(shù)量級。

    為比較搜索效率,圖6給出了相同迭代次數(shù)下遺傳模擬退火算法與單一遺傳算法(GA)及加入SQP算法前后的搜索效果對比。從圖6中可看出,GSA與GA算法在前10次運算中效率相當;之后,GSA引入的模擬退火過程提高了搜索效率;40次運算后,SQP算法的引入進一步提高了搜索效率,算法快速達到了收斂狀態(tài)。

    圖6 GA, GSA與SQP算法效果對比

    加入SQP算法前后的氣體溫度重建結果及偏差如圖7所示,從圖中可看出,由于GSA算法的搜索深度不足,重建結果僅僅在整體上接近了模型分布;加入SQP算法后,重建結果有了顯著改善。

    (a) SQP加入前溫度重建結果 (b) SQP加入前溫度重建偏差

    (c) SQP加入后溫度重建結果 (d) SQP加入后溫度重建偏差

    采用歸一化平均絕對誤差對重建效果進行評價,其表達式為

    為驗證方法的可行性及SQP算法加入前后的重建效果,構建不同溫度范圍的溫度模型并按上述方法進行仿真驗證實驗,仿真結果如表3所示。

    表3 不同溫度模型下SQP算法加入前后的重建結果

    3.3 重建穩(wěn)定性分析

    因?qū)嶋H測量中存在鏡面污垢、平臺振動、擬合偏差等因素,使得測量結果難免會有誤差;在此對測量積分吸光度中加入不同比例的隨機噪聲,以考察重建的穩(wěn)定性。圖8為HT重建與雙線法重建在0~20%噪聲下的溫度重建對比,其中雙線法采用了4個投影角度,共32條光線。圖8顯示,在無噪聲的情況下,雙線法的重建結果與HT法相當,甚至優(yōu)于HT法;但隨著噪聲的增加,雙線法的穩(wěn)定性逐漸降低,重建誤差逐步偏大。

    圖 8 不同比例噪聲下的重建結果

    HT法在不同比例噪聲下的重建結果較為穩(wěn)定,分析其原因為:1)HT法采用多條吸收譜線對同一光線路徑進行測量,隨機噪聲的影響因存在相互抵消而減??;2)平滑正則化函數(shù)使重建結果趨于平滑分布,進一步降低了噪聲帶來的影響。

    4 結束語

    本文設計了基于多波長激光吸收光譜技術的氣體溫度二維重建系統(tǒng),數(shù)值模擬了8×8的氣體溫度場并進行了重建驗證。重建算法上結合了遺傳模擬退火算法與序列二次規(guī)劃算法,提高了搜索效率。通過在測量信息中加入不同比例的噪聲信號,驗證了算法的穩(wěn)定性,并與雙線法重建結果進行了對比。實驗表明,雙線法在較低噪聲情況下的重建結果與HT法相當,但隨著噪聲的增加,其重建誤差逐漸變大;相反HT法則給出了較穩(wěn)定的重建結果。平滑正則化函數(shù)的引入對重建結果有著兩方面的影響:1)改良了HT問題的病態(tài)特性,降低了重建難度,減小了隨機噪聲對重建結果的影響;2) 降低了在峰值處的重建結果,加大了重建誤差。

    (References)

    [[1]] 洪延姬, 宋俊玲, 王廣宇, 等. 激光吸收光譜技術測量非均勻燃燒流場研究進展[J]. 航空學報, 2015, 36(3): 724-736

    Hong Yanji, Song Junling, Wang Guangyu, et al. Review on the research of non-uniform combustion field measurement using laser absorption spectroscopy technique[J]. Acta Aeronautica et Astronautica Sinica, 2015, 36(3): 724-736.

    [2] 宋俊玲, 洪延姬, 王廣宇, 等. 基于激光吸收光譜技術的燃燒場氣體溫度和濃度二維分布重建研究[J]. 物理學報, 2012, 61(24): 240702

    Song Junling, Hong Yanji, Wang Guangyu, et al. Two-dimensional reconstructions of gas temperature and concentration in combustion based on tunable diode laser absorption spectroscopy[J]. Acta Phys Sin, 2012, 61(24): 240702

    [3] Wang F, Cen K F, Li N, et al. Two-dimensional tomography for gas concentration and temperature distributions based on tunable diode laser absorption

    spectroscopy[J]. Measurement Science & Technology, 2010, 21(4): 45301-45310

    [4] Ma L, Cai W. Numerical investigation of hyperspectral tomography for simultaneous temperature and concentration imaging[J].Applied Optics, 2008, 47(21): 3751-3759

    [5] Cai W, Ma L. Hyperspectral tomography based on proper orthogonal decomposition as motivated by imaging diagnostics of unsteady reactive flows[J]. Applied Optics, 2010, 49(4): 601-610

    [6]Cai W, Kaminski C F. A numerical investigation of high-resolution multispectral absorption tomography for flow thermometry[J]. Applied Physics B, 2015, 119(1): 29-35

    [7] 焦萌. 調(diào)諧二極管激光吸收層析成像技術研究進展[J]. 激光雜志, 2014(2): 1-3

    Jiao Meng. Review of tunable diode laser absorption tomography[J]. Laser Journal, 2014(2): 1-3

    [8]夏暉暉, 劉建國, 許振宇, 等. 基于代數(shù)迭代算法的燃燒火焰溫度場和氣體濃度場重建研究[J]. 光譜學與光譜分析, 2015(10): 2697-2702

    Xia Huihui, Liu Jianguo, Xu Zhenyu, et al. Reconstruction research for gas concentration and temperature of flame based on algebraic reconstruction technique[J]. Spectroscopy and Spectral Analysis, 2015(10): 2697-2702

    [9]李寧, 翁春生. 基于多波長激光吸收光譜技術的氣體濃度與溫度二維分布遺傳模擬退火重建研究[J]. 物理學報, 2010, 59(10): 6914-6920

    Li Ning, Weng Chunsheng. Gas concentration and temperature reconstruction by genetic simulated annealing algorithm based on multi-wavelengths diode laser absorption spectroscopy[J]. Acta Phys Sin, 2010, 59(10): 6914-6920

    [10] Gamache R R, Kennedy S, Hawkins R, et al. Total internal partition sums for molecules in the terrestrial atmosphere[J]. Journal of Molecular Structure, 2000, 517: 407-425

    [11] Rothman L S, Gordon I E, Barbe A, et al. The HITRAN, 2008 molecular spectroscopic database[J]. Journal of Quantitative Spectroscopy & Radiative Transfer, 2013, 130(11): 4-50

    [12] Zhou X, Liu X, Jeffries J B, et al. Development of a sensor for temperature and water concentration in combustion gases using a single tunable diode laser[J]. Measurement Science and Technology, 2003, 14(8): 1459

    (編輯:肖福根)

    Two-dimensional reconstruction of gas temperature based on multi-wavelength laser absorption spectroscopy

    Xin Mingyuan, Jin Xing, Song Junling

    (State Key Laboratory of Laser Propulsion and Applications, Equipment Academy, Beijing 101416, China)

    The laser absorption spectroscopy has a broad application prospect in the field of Aeronautics and Astronautics. In this paper, a two-dimensional gas temperature is reconstructed by the multi-wavelengths diode laser absorption spectroscopy. A measurement system for the gas temperature is designed, and the number of projections is reduced by increasing the information of the absorption spectrum. A hybrid algorithm which combines the genetic simulated annealing with the sequential quadratic programming (SQP) is used to solve the problem, which improves the efficiency of the search. A combustion model of the product including water is established, and the temperature of the combustion field is reconstructed by using 8 near-infrared absorption spectrum lines of water. The stability of the algorithm is verified by adding different proportions of random noise in the measurement information, and the results are compared with the two wavelength absorption method. The experimental results show that the reconstructions of the two methods are equivalent in the low-noise environment, but with the increase of the noise, the multi-wavelength absorption method performs more stable.

    laser absorption spectroscopy; multi-wavelengths absorption; hyperspectral tomography; genetic simulated annealing algorithm; two-dimensional reconstruction

    TN247;O433.1

    A

    1673-1379(2016)03-0269-06

    10.3969/j.issn.1673-1379.2016.03.007

    辛明原(1991—),男,碩士研究生,研究方向為激光流場診斷技術;E-mail: ybgaxmy@163.com。指導教師:金 星(1962—),男,研究員,研究領域為激光航天應用技術。

    2016-03-15;

    2016-05-06

    國家自然科學青年基金項目(編號:61505263)

    http://www.bisee.ac.cn E-mail: htqhjgc@126.com Tel: (010)68116407, 68116408, 68116544

    猜你喜歡
    模擬退火吸收光譜譜線
    基于HITRAN光譜數(shù)據(jù)庫的合并譜線測溫仿真研究
    模擬退火遺傳算法在機械臂路徑規(guī)劃中的應用
    測控技術(2018年3期)2018-11-25 09:45:08
    鐵合金光譜譜線分離實驗研究
    電子測試(2018年11期)2018-06-26 05:56:00
    原子吸收光譜分析的干擾與消除應用研究
    淺析原子吸收光譜法在土壤環(huán)境監(jiān)測中的應用
    鍶原子光鐘鐘躍遷譜線探測中的程序控制
    基于模糊自適應模擬退火遺傳算法的配電網(wǎng)故障定位
    茶油氧化過程中紫外吸收光譜特性
    SOA結合模擬退火算法優(yōu)化電容器配置研究
    電源技術(2015年5期)2015-08-22 11:18:24
    基于遺傳-模擬退火算法的城市軌道交通快慢車停站方案
    日韩不卡一区二区三区视频在线| 黄网站色视频无遮挡免费观看| 免费播放大片免费观看视频在线观看| 亚洲成色77777| 黄片播放在线免费| 1024视频免费在线观看| 久久久久久久久久人人人人人人| 亚洲四区av| 久久影院123| 2018国产大陆天天弄谢| 美女国产视频在线观看| 如日韩欧美国产精品一区二区三区| 国产免费福利视频在线观看| 高清毛片免费看| 麻豆乱淫一区二区| 97精品久久久久久久久久精品| 美国免费a级毛片| 成人国产av品久久久| 亚洲av男天堂| 熟女电影av网| 久久精品aⅴ一区二区三区四区 | 91精品伊人久久大香线蕉| 久久久久久伊人网av| 永久免费av网站大全| 久久久久国产网址| 久久久精品免费免费高清| 国产精品麻豆人妻色哟哟久久| 美女国产视频在线观看| 伦理电影大哥的女人| 免费人成在线观看视频色| tube8黄色片| 精品一区二区三区视频在线| 亚洲色图 男人天堂 中文字幕 | 嫩草影院入口| 黑人高潮一二区| 91aial.com中文字幕在线观看| 国产在线免费精品| 考比视频在线观看| 国产在线一区二区三区精| 夜夜骑夜夜射夜夜干| 一区二区三区乱码不卡18| 各种免费的搞黄视频| 免费观看在线日韩| 80岁老熟妇乱子伦牲交| 韩国高清视频一区二区三区| 亚洲av电影在线观看一区二区三区| 午夜日本视频在线| 成年人免费黄色播放视频| 久久精品国产亚洲av天美| 美女内射精品一级片tv| 日韩大片免费观看网站| 日韩av不卡免费在线播放| 精品人妻熟女毛片av久久网站| 亚洲图色成人| 婷婷色综合www| av一本久久久久| 少妇被粗大猛烈的视频| 久久久久精品人妻al黑| 午夜福利视频在线观看免费| 熟女人妻精品中文字幕| av国产久精品久网站免费入址| av免费在线看不卡| 婷婷色麻豆天堂久久| 深夜精品福利| 精品视频人人做人人爽| 又黄又粗又硬又大视频| 亚洲欧美成人综合另类久久久| 亚洲av中文av极速乱| 亚洲av国产av综合av卡| 人成视频在线观看免费观看| 欧美精品高潮呻吟av久久| av黄色大香蕉| 日韩电影二区| 久久这里有精品视频免费| 国产成人午夜福利电影在线观看| 亚洲精品一区蜜桃| 亚洲高清免费不卡视频| 久久免费观看电影| 国产高清不卡午夜福利| 人体艺术视频欧美日本| 久久女婷五月综合色啪小说| 免费av不卡在线播放| 亚洲第一av免费看| 国产一区有黄有色的免费视频| 欧美xxxx性猛交bbbb| 日韩av免费高清视频| 男女高潮啪啪啪动态图| 人人妻人人添人人爽欧美一区卜| 免费大片18禁| 亚洲综合色网址| 这个男人来自地球电影免费观看 | 日韩制服骚丝袜av| 在线观看三级黄色| 在线观看国产h片| 综合色丁香网| 亚洲欧美成人综合另类久久久| 亚洲精品自拍成人| 免费女性裸体啪啪无遮挡网站| 久久久国产精品麻豆| 亚洲精品美女久久久久99蜜臀 | 青春草国产在线视频| 韩国高清视频一区二区三区| 99久久综合免费| 国产高清国产精品国产三级| 欧美精品国产亚洲| 久久这里只有精品19| freevideosex欧美| 亚洲国产精品一区二区三区在线| 精品酒店卫生间| 欧美老熟妇乱子伦牲交| 久久狼人影院| 亚洲精品一区蜜桃| 18禁动态无遮挡网站| 久久久精品区二区三区| 亚洲欧美中文字幕日韩二区| 亚洲精品av麻豆狂野| 中文精品一卡2卡3卡4更新| 国产女主播在线喷水免费视频网站| 久久鲁丝午夜福利片| 欧美老熟妇乱子伦牲交| 80岁老熟妇乱子伦牲交| 精品一区二区三卡| av在线app专区| 在线亚洲精品国产二区图片欧美| 久久久精品区二区三区| 五月天丁香电影| 热99国产精品久久久久久7| 水蜜桃什么品种好| 国产无遮挡羞羞视频在线观看| 高清毛片免费看| 免费黄频网站在线观看国产| 建设人人有责人人尽责人人享有的| 曰老女人黄片| 最近手机中文字幕大全| 欧美xxxx性猛交bbbb| 亚洲少妇的诱惑av| 久久久欧美国产精品| 深夜精品福利| 欧美人与善性xxx| 少妇的逼水好多| 99国产精品免费福利视频| 欧美人与善性xxx| 日本-黄色视频高清免费观看| 亚洲av电影在线进入| 80岁老熟妇乱子伦牲交| 伊人久久国产一区二区| 国产成人欧美| 高清毛片免费看| 精品人妻在线不人妻| 免费看不卡的av| 色5月婷婷丁香| 日韩欧美一区视频在线观看| 国产精品人妻久久久久久| 久久这里有精品视频免费| 久久久a久久爽久久v久久| 国产乱人偷精品视频| 久久久久久久亚洲中文字幕| 丁香六月天网| 欧美日韩av久久| 人人妻人人添人人爽欧美一区卜| 亚洲精品国产av成人精品| 韩国高清视频一区二区三区| 高清在线视频一区二区三区| 黄网站色视频无遮挡免费观看| 色婷婷久久久亚洲欧美| 国产黄色视频一区二区在线观看| 日韩 亚洲 欧美在线| 在线免费观看不下载黄p国产| 少妇熟女欧美另类| 在线免费观看不下载黄p国产| 少妇熟女欧美另类| 少妇 在线观看| videossex国产| 三级国产精品片| 精品一区二区三卡| 天天躁夜夜躁狠狠躁躁| 国产精品一二三区在线看| 最后的刺客免费高清国语| 久久午夜福利片| 好男人视频免费观看在线| 国产亚洲精品久久久com| 久久韩国三级中文字幕| 精品少妇内射三级| 国产精品熟女久久久久浪| 亚洲精品成人av观看孕妇| 王馨瑶露胸无遮挡在线观看| 久久青草综合色| 黑人高潮一二区| 日韩中文字幕视频在线看片| 久久97久久精品| 国产日韩欧美视频二区| 国产精品嫩草影院av在线观看| 久久精品夜色国产| 咕卡用的链子| 视频中文字幕在线观看| 看免费av毛片| 曰老女人黄片| 亚洲伊人色综图| 18禁国产床啪视频网站| 午夜福利视频在线观看免费| 亚洲精品久久午夜乱码| 成人漫画全彩无遮挡| 日本爱情动作片www.在线观看| 国产精品一区www在线观看| 激情视频va一区二区三区| 精品一品国产午夜福利视频| 看十八女毛片水多多多| 日本免费在线观看一区| 欧美人与善性xxx| 国产av码专区亚洲av| 在线天堂最新版资源| 国产福利在线免费观看视频| 免费人妻精品一区二区三区视频| 午夜激情久久久久久久| 嫩草影院入口| av女优亚洲男人天堂| 一边摸一边做爽爽视频免费| 欧美3d第一页| 久久久久久久国产电影| 国产欧美日韩一区二区三区在线| 亚洲国产色片| 黄网站色视频无遮挡免费观看| 国产xxxxx性猛交| 亚洲人成77777在线视频| 免费高清在线观看日韩| 岛国毛片在线播放| 免费av不卡在线播放| 男女免费视频国产| 久久99精品国语久久久| 精品一区二区免费观看| 欧美国产精品一级二级三级| 日日爽夜夜爽网站| av电影中文网址| a级毛色黄片| 美女视频免费永久观看网站| 秋霞伦理黄片| 成人亚洲欧美一区二区av| 只有这里有精品99| 日韩伦理黄色片| 我的女老师完整版在线观看| 欧美另类一区| 亚洲精品av麻豆狂野| 国内精品宾馆在线| 国产乱来视频区| 日韩制服丝袜自拍偷拍| 日韩精品免费视频一区二区三区 | 18禁国产床啪视频网站| 国产白丝娇喘喷水9色精品| 狠狠婷婷综合久久久久久88av| 成年美女黄网站色视频大全免费| 国产黄色免费在线视频| 一区二区av电影网| 九色成人免费人妻av| 日本爱情动作片www.在线观看| 国产在视频线精品| 亚洲人与动物交配视频| 国产一级毛片在线| 男女午夜视频在线观看 | 少妇的丰满在线观看| 黑人欧美特级aaaaaa片| 国产成人a∨麻豆精品| 久久久久精品性色| 我要看黄色一级片免费的| 99热全是精品| 亚洲三级黄色毛片| 亚洲第一区二区三区不卡| 七月丁香在线播放| 久久久久久伊人网av| 2018国产大陆天天弄谢| 黄色怎么调成土黄色| 五月天丁香电影| 久久精品国产鲁丝片午夜精品| 亚洲精品美女久久av网站| 成人毛片a级毛片在线播放| 又大又黄又爽视频免费| 女人精品久久久久毛片| 国产成人av激情在线播放| 久久人人爽人人片av| 另类亚洲欧美激情| 伊人亚洲综合成人网| 一本色道久久久久久精品综合| 国产免费又黄又爽又色| 色5月婷婷丁香| 边亲边吃奶的免费视频| 国产一区二区三区av在线| 国产男人的电影天堂91| 免费久久久久久久精品成人欧美视频 | 99热这里只有是精品在线观看| 国产xxxxx性猛交| 国产成人精品在线电影| 日本wwww免费看| 欧美日韩视频高清一区二区三区二| 免费久久久久久久精品成人欧美视频 | 晚上一个人看的免费电影| 成人二区视频| 激情视频va一区二区三区| 日韩精品免费视频一区二区三区 | 丝袜喷水一区| 久久久久久久久久人人人人人人| 26uuu在线亚洲综合色| 五月玫瑰六月丁香| 蜜桃在线观看..| 美女xxoo啪啪120秒动态图| 精品第一国产精品| 丰满少妇做爰视频| 黄色视频在线播放观看不卡| 久久久久久人妻| 美女中出高潮动态图| 在线精品无人区一区二区三| av视频免费观看在线观看| 欧美最新免费一区二区三区| 免费人妻精品一区二区三区视频| 一边亲一边摸免费视频| 最近中文字幕高清免费大全6| 亚洲精品色激情综合| 国产精品欧美亚洲77777| 久久久久久久精品精品| 男女午夜视频在线观看 | 黑人高潮一二区| 午夜老司机福利剧场| 亚洲人成网站在线观看播放| 成人毛片a级毛片在线播放| 中文字幕最新亚洲高清| 纯流量卡能插随身wifi吗| 插逼视频在线观看| 久久久精品94久久精品| 丰满饥渴人妻一区二区三| 久久精品国产亚洲av涩爱| 69精品国产乱码久久久| 久久99热这里只频精品6学生| 久久久久久久大尺度免费视频| 久久鲁丝午夜福利片| 日韩精品有码人妻一区| 国产成人免费无遮挡视频| 精品99又大又爽又粗少妇毛片| 亚洲精品美女久久久久99蜜臀 | 国产福利在线免费观看视频| 国产一区二区在线观看av| 国产成人a∨麻豆精品| 2021少妇久久久久久久久久久| 亚洲成人一二三区av| 亚洲精品美女久久av网站| 在线观看三级黄色| 99九九在线精品视频| www日本在线高清视频| 免费黄频网站在线观看国产| 久久久欧美国产精品| 交换朋友夫妻互换小说| 成人黄色视频免费在线看| 精品少妇久久久久久888优播| 国产在视频线精品| 日韩制服丝袜自拍偷拍| 国产精品不卡视频一区二区| 日产精品乱码卡一卡2卡三| 日本猛色少妇xxxxx猛交久久| 国产 精品1| 熟女av电影| 婷婷色综合大香蕉| 亚洲欧美日韩卡通动漫| 国产精品久久久久久久久免| 国产在线视频一区二区| 日韩精品有码人妻一区| 美女福利国产在线| 久久精品国产a三级三级三级| 五月开心婷婷网| av免费观看日本| 人人妻人人添人人爽欧美一区卜| 欧美3d第一页| 日本黄色日本黄色录像| 99热网站在线观看| 亚洲成人手机| 丁香六月天网| 大陆偷拍与自拍| 欧美国产精品一级二级三级| 亚洲精品国产av蜜桃| 22中文网久久字幕| 最近中文字幕高清免费大全6| 国产又爽黄色视频| 久久99蜜桃精品久久| 欧美人与性动交α欧美软件 | 在线天堂最新版资源| 日本猛色少妇xxxxx猛交久久| 欧美 日韩 精品 国产| 欧美激情国产日韩精品一区| 亚洲精品色激情综合| 七月丁香在线播放| 精品少妇内射三级| 亚洲图色成人| 人人妻人人澡人人爽人人夜夜| 精品一区二区免费观看| 最后的刺客免费高清国语| 国产深夜福利视频在线观看| 久久久久国产网址| 国产精品久久久久久久电影| 十分钟在线观看高清视频www| 国产成人a∨麻豆精品| 午夜精品国产一区二区电影| 成年美女黄网站色视频大全免费| 激情五月婷婷亚洲| 青春草国产在线视频| 国产白丝娇喘喷水9色精品| 日韩中文字幕视频在线看片| 9热在线视频观看99| 人妻一区二区av| 成人手机av| 国精品久久久久久国模美| 亚洲图色成人| 人成视频在线观看免费观看| 丝袜美足系列| 日韩视频在线欧美| 免费播放大片免费观看视频在线观看| 91精品国产国语对白视频| 天天躁夜夜躁狠狠躁躁| 人体艺术视频欧美日本| 七月丁香在线播放| 精品一区二区三区视频在线| 久久国产精品男人的天堂亚洲 | 欧美精品av麻豆av| 一边摸一边做爽爽视频免费| 中国美白少妇内射xxxbb| av播播在线观看一区| 国精品久久久久久国模美| 高清在线视频一区二区三区| 天堂中文最新版在线下载| 亚洲综合色网址| 最近2019中文字幕mv第一页| a级片在线免费高清观看视频| 日产精品乱码卡一卡2卡三| 天堂8中文在线网| 黄色 视频免费看| 国产av码专区亚洲av| 国产一区亚洲一区在线观看| 黄色视频在线播放观看不卡| 看十八女毛片水多多多| 在线观看国产h片| 少妇的逼水好多| 国产成人91sexporn| 校园人妻丝袜中文字幕| 乱码一卡2卡4卡精品| 91午夜精品亚洲一区二区三区| 黑人猛操日本美女一级片| 麻豆精品久久久久久蜜桃| 欧美日本中文国产一区发布| 国产精品人妻久久久影院| 狠狠精品人妻久久久久久综合| 日韩一本色道免费dvd| 久久精品国产综合久久久 | 美女国产视频在线观看| 日韩一区二区三区影片| 亚洲经典国产精华液单| 桃花免费在线播放| 人妻人人澡人人爽人人| 亚洲人成网站在线观看播放| 免费不卡的大黄色大毛片视频在线观看| 亚洲成人一二三区av| www.av在线官网国产| 啦啦啦在线观看免费高清www| 国产69精品久久久久777片| 国产伦理片在线播放av一区| www.色视频.com| 少妇 在线观看| 免费高清在线观看日韩| 91久久精品国产一区二区三区| 久久久久久久大尺度免费视频| 成年美女黄网站色视频大全免费| 国产精品熟女久久久久浪| 国产免费一区二区三区四区乱码| 免费少妇av软件| freevideosex欧美| 嫩草影院入口| 丁香六月天网| 2022亚洲国产成人精品| 黄色毛片三级朝国网站| 国产又爽黄色视频| 最黄视频免费看| 夜夜骑夜夜射夜夜干| h视频一区二区三区| 看非洲黑人一级黄片| 日韩欧美精品免费久久| 一级爰片在线观看| 亚洲精品视频女| 亚洲欧美日韩卡通动漫| 夜夜爽夜夜爽视频| 男人操女人黄网站| 亚洲欧美清纯卡通| 欧美+日韩+精品| 亚洲av男天堂| 韩国av在线不卡| 插逼视频在线观看| 男女国产视频网站| 国产精品蜜桃在线观看| 成年美女黄网站色视频大全免费| 久久ye,这里只有精品| 亚洲国产成人一精品久久久| 成年美女黄网站色视频大全免费| 中国国产av一级| av在线app专区| 亚洲色图 男人天堂 中文字幕 | 少妇熟女欧美另类| 波野结衣二区三区在线| 高清av免费在线| 亚洲精品456在线播放app| 国产av精品麻豆| 日本爱情动作片www.在线观看| 国产爽快片一区二区三区| 亚洲av日韩在线播放| 视频中文字幕在线观看| av在线老鸭窝| 亚洲精品自拍成人| 在线观看免费高清a一片| 成年动漫av网址| 国产极品粉嫩免费观看在线| 国产精品99久久99久久久不卡 | 韩国精品一区二区三区 | 亚洲高清免费不卡视频| 另类亚洲欧美激情| 亚洲精品乱久久久久久| 热re99久久精品国产66热6| 日本免费在线观看一区| av在线播放精品| 午夜福利乱码中文字幕| 精品国产一区二区久久| 最近的中文字幕免费完整| 亚洲成色77777| 观看美女的网站| av黄色大香蕉| 国产精品无大码| 久久久久网色| 伦理电影大哥的女人| 男女国产视频网站| 日本91视频免费播放| 两性夫妻黄色片 | 亚洲美女搞黄在线观看| 国产黄色免费在线视频| 久久国内精品自在自线图片| av有码第一页| 色哟哟·www| 久久久久久人妻| 免费人妻精品一区二区三区视频| 中文字幕精品免费在线观看视频 | 自拍欧美九色日韩亚洲蝌蚪91| 久久久久久人妻| 久久99蜜桃精品久久| 熟女人妻精品中文字幕| 亚洲第一区二区三区不卡| 日日爽夜夜爽网站| 中国美白少妇内射xxxbb| 高清av免费在线| 精品人妻一区二区三区麻豆| 26uuu在线亚洲综合色| 啦啦啦啦在线视频资源| 亚洲av日韩在线播放| 精品卡一卡二卡四卡免费| 久久精品国产亚洲av涩爱| 97在线视频观看| 涩涩av久久男人的天堂| 亚洲欧美成人综合另类久久久| 国产精品国产三级国产专区5o| 交换朋友夫妻互换小说| 欧美日韩国产mv在线观看视频| 少妇高潮的动态图| 夫妻午夜视频| 国产一区有黄有色的免费视频| 国产精品秋霞免费鲁丝片| 国产成人欧美| av网站免费在线观看视频| 大片免费播放器 马上看| 亚洲成av片中文字幕在线观看 | 插逼视频在线观看| 精品少妇内射三级| 欧美日本中文国产一区发布| 美女xxoo啪啪120秒动态图| 国产午夜精品一二区理论片| 一级片'在线观看视频| 中文字幕亚洲精品专区| 国产极品粉嫩免费观看在线| 内地一区二区视频在线| 日韩不卡一区二区三区视频在线| 亚洲成人一二三区av| 精品人妻偷拍中文字幕| 国产 精品1| 亚洲精品aⅴ在线观看| 久久久久网色| 欧美xxⅹ黑人| 亚洲人与动物交配视频| 精品福利永久在线观看| 一本大道久久a久久精品| 欧美最新免费一区二区三区| 久久97久久精品| 久久精品国产亚洲av涩爱| av天堂久久9| 久久这里有精品视频免费| 国产国拍精品亚洲av在线观看| 91成人精品电影| 国产成人91sexporn| 天天躁夜夜躁狠狠躁躁| 精品国产乱码久久久久久小说| 亚洲国产成人一精品久久久| 免费高清在线观看视频在线观看| 黄色视频在线播放观看不卡| 永久免费av网站大全| 精品久久久精品久久久| 亚洲成人av在线免费| 国产免费一区二区三区四区乱码| 伦理电影大哥的女人| 一区二区三区四区激情视频| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 一区二区日韩欧美中文字幕 | 91精品伊人久久大香线蕉| 女人久久www免费人成看片| 国产精品女同一区二区软件| 国产成人一区二区在线| 伦理电影免费视频| 人人妻人人爽人人添夜夜欢视频| 亚洲国产精品一区三区|