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

    非平表面的紅外熱像測溫修正方法研究

    2021-03-06 02:47:34付萬超范春利
    紅外技術(shù) 2021年2期
    關(guān)鍵詞:熱像儀頂角發(fā)射率

    付萬超,范春利,楊 立

    〈測量技術(shù)〉

    非平表面的紅外熱像測溫修正方法研究

    付萬超,范春利,楊 立

    (海軍工程大學 動力工程學院,湖北 武漢 430033)

    利用紅外熱像儀測溫需先設(shè)定被測表面的法向發(fā)射率,該發(fā)射率通常為定值。而當熱像儀處于被測點的天頂角大于50°的位置范圍時,由于被測點定向發(fā)射率的變化,必造成這些點的測溫誤差。對于非平表面,這樣的點大量存在。因此,必須對其測溫結(jié)果進行修正。本文針對使用單目紅外熱像儀測量非平表面溫度時由于各點定向發(fā)射率的變化引起的測量誤差進行研究,并依據(jù)物體表面定量發(fā)射率的變化規(guī)律,給出了測量點的溫度修正系數(shù)。同時,通過點云三維建模,利用熱像儀的幾何成像原理推導(dǎo)出紅外熱像圖與實際被測表面中點與點的對應(yīng)關(guān)系,給出了通過紅外熱像儀測量非平表面的溫度分布的誤差修正方法。實驗證明了該方法的有效性。

    紅外無損檢測;紅外熱像儀;測溫誤差;溫度修正;非平表面

    0 引言

    在工業(yè)生產(chǎn)中,及時掌握設(shè)備的工作狀態(tài)是十分必要的。好多情況下,物體表面溫度可以間接地反映設(shè)備的工作狀態(tài),因此開展依據(jù)物體表面溫度對工業(yè)設(shè)備運行狀態(tài)進行檢測和監(jiān)測的理論與方法研究具有重要的意義[1-2]。

    對于設(shè)備表面溫度分布的測量,最常用的設(shè)備為紅外熱像儀。通過紅外熱像儀進行測溫,一般要求正對測量表面。然而在狹小空間或由其他原因造成可達性較差的位置進行測量時,如艦船的機艙,常常無法滿足紅外熱成像測溫的要求,因此就需要人們對特殊測量條件下所測得的溫度分布進行修正。造成紅外熱像測溫誤差的因素很多,相關(guān)的研究包括紅外熱像儀的測溫的影響因素分析及誤差的計算[1,3]、太陽、大氣等環(huán)境的影響[1,4-5],周圍高溫物體的影響[6-7]、測量距離及視角的影響[8-9]、探測器溫度的影響[10]等。對于本文討論的非平設(shè)備表面,尤其在熱像儀處于被測點的天頂角的50°范圍之外時,發(fā)射率將發(fā)生明顯的變化。因此,對于設(shè)為恒定表面發(fā)射率的熱像儀,必將會產(chǎn)生顯著的測溫誤差。因此,開展非平表面紅外溫度測量的修正研究具有重要的應(yīng)用價值。

    本文針對利用紅外熱像儀大角度測溫或?qū)Ψ瞧奖砻孢M行測溫因無法或部分無法處于表面各點定向發(fā)射率恒定的天頂角范圍內(nèi)時熱像儀的測溫誤差進行研究,并基于物體表面定向發(fā)射率的變化規(guī)律,通過對所測表面進行3D建模,建立實際表面和熱像圖之間的坐標對應(yīng)關(guān)系,提出熱像儀所測溫度分布的修正方法。方法通過實驗進行驗證。

    1 紅外測溫及修正系數(shù)的計算

    根據(jù)紅外熱像儀溫度計算公式[1]:

    式中:為物體表面發(fā)射率;a為大氣的光譜透射率;r為熱像儀指示的輻射溫度;a為大氣輻射率;u為環(huán)境溫度。在3~5mm,=9.2554;在8~12mm,=3.9889。

    由式(1)可知,溫度測量值與被測物體的發(fā)射率密切相關(guān)。在測量復(fù)雜物體表面時,需在熱像儀上將發(fā)射率設(shè)為一固定值,通常為該表面材料的法向發(fā)射率。天頂角是某個方向和天頂方向的夾角,通常記為(見圖1)。本文中利用天頂角來描述熱像儀測量物體表面中點的溫度時,熱像儀成像中心點相對于物體表面的位置關(guān)系。在一次測量時,熱像儀所處的位置,對于非平表面上不同的點具有不同的天頂角。當角超出一定范圍(如大于50°)時,定向發(fā)射率的變化必然對紅外測溫的準確性產(chǎn)生影響。

    圖1 熱像儀機位與其對應(yīng)的天頂角

    根據(jù)菲涅爾定律[1],光譜方向反射率表示為:

    根據(jù)Snell定律,折射角與入射角滿足:

    根據(jù)Kirchhoff定律:

    據(jù)表1[11]中鋁(Al)、黃銅(Brass)、聚乙烯(PS)、聚苯乙烯(PE)4種材料分別在3~5mm與8~12mm工作波段的復(fù)折射率可計算出其定向發(fā)射率隨天頂角變化的理論值。

    表1 材料各波段復(fù)折射率

    式(4)的計算結(jié)果如圖2所示,可觀察到金屬與非導(dǎo)電體定向發(fā)射率隨天頂角的變化規(guī)律:各類材料定向發(fā)射率在被測點天頂角大于50°附近開始出現(xiàn)明顯的變化。

    在工程測量中,需在測量前設(shè)定被測物體表面的發(fā)射率,通常設(shè)為其法向發(fā)射率^。對于相對平整的表面,可使用^代表所有點的發(fā)射率。當熱像儀處于被測點的天頂角大于50°的位置范圍時,需采用該角度的定向發(fā)射率。對于不同角度的定向發(fā)射率,根據(jù)(1)式可知實際溫度0與測量溫度1的關(guān)系為:

    (a) 金屬材料發(fā)射率隨天頂角變化

    (a) Variation of emissivity of metallic materials with zenith angle

    (b) 非金屬材料發(fā)射率隨天頂角變化

    (b) Variation of emissivity of nonmetallic materials with zenith angle

    圖2 定向發(fā)射率隨天頂角的變化

    Fig.2 Variation of direction emissivity with zenith angle

    定義修正系數(shù)為:

    則可得修正系數(shù)曲線如圖3所示。

    由圖3觀察得出,對于金屬表面溫度測量結(jié)果的修正系數(shù)在天頂角大于50°時減小,在90°附近上升迅速,而非導(dǎo)電體在天頂角大于70°時迅速增加,最大時測溫誤差可達50%。因此,需要對紅外熱像圖中不同位置的溫度值依據(jù)其對應(yīng)實際點的修正系數(shù)進行修正。

    圖3 溫度修正系數(shù)隨天頂角的變化

    2 被測表面各點天頂角的計算

    如前所述,要計算被測表面各個點的修正系數(shù),必須首先計算各個點紅外測溫時熱像儀所處的天頂角。

    2.1 坐標與法向量的獲取

    以一鍋爐為例,繪制出其點云坐標和法向量如圖4所示。

    圖4 某鍋爐點云坐標與法向量示例

    2.2 天頂角的計算

    則可計算天頂角為:

    由空間幾何關(guān)系可知,對于天頂角小于90°的被測點,由于透視關(guān)系無法在熱像儀上成像,因此可得可視部分如圖5所示。

    其中坐標為(-0.75,0.12,0.1)的點為紅外熱像儀中心所在點,黃色部分為鍋爐表面在該位置關(guān)系下熱像儀的可視點,藍色為不可視點。

    圖5 鍋爐的紅外熱像儀可視點

    2.3 溫度修正系數(shù)的分布

    由該位置關(guān)系可計算得出各點對應(yīng)天頂角。由于實際工程應(yīng)用中,通常在金屬表面填涂一層非導(dǎo)電涂料,且非導(dǎo)電體在法向發(fā)射率相較于導(dǎo)電體更大,理論的光譜定向發(fā)射率與實際相差不大,因此以8~13mm的PS( Poly styrene)材料作為被測物體材料。在圖4的位置關(guān)系下,計算得出該鍋爐的修正系數(shù)的頻率分布直方圖如圖6所示。經(jīng)計算可得,該位置關(guān)系下修正系數(shù)大于1.2的點占比44.1%。由此可見,對實際非平表面測溫時的定向發(fā)射率進行修正是非常必要的。

    圖6 不同溫度修正系數(shù)對應(yīng)的被測點所占總被測點的比例

    3 被測點在圖像的映射

    天頂角與成像點的實際位置和其法向量密切相關(guān)。因此,針對紅外熱像儀所測溫度的修正需明確各成像點在世界坐標中的實際位置。

    世界坐標系下被測點w經(jīng)紅外熱像儀投影至平面像素點p=[,]T,在成像過程中經(jīng)過兩次坐標變換,分別是由世界坐標系到熱像儀坐標系的變換和由熱像儀坐標系到像素坐標系的變換[12],兩個變換分別通過外參矩陣和內(nèi)參矩陣來實現(xiàn)。

    3.1 外參矩陣的獲取

    世界坐標w到熱像儀坐標i的轉(zhuǎn)換關(guān)系為:

    其中:

    為熱像儀的外參矩陣,i為熱像儀的基坐標在世界坐標系的坐標,描述熱像儀在世界坐標系中的姿態(tài),為熱像儀中心在世界坐標系中的位置。

    計算可得R=[1,2,3]為熱像儀在世界坐標系旋轉(zhuǎn)到當前姿態(tài)的旋轉(zhuǎn)矩陣。

    當熱像儀角度變化時,可根據(jù)熱像儀的俯角與轉(zhuǎn)角求得熱像儀坐標系軸與水平面的交點:

    3.2 內(nèi)參矩陣的獲取

    在世界坐標轉(zhuǎn)換為熱像儀坐標后,需將熱像儀坐標轉(zhuǎn)換為像素坐標,完成三維到二維的轉(zhuǎn)變。

    熱像儀坐標i到像素坐標p的轉(zhuǎn)換關(guān)系為:

    內(nèi)參矩陣與熱像儀的型號和圖像的分辨率有關(guān),kk為熱像儀與方向的像素焦距,cc為圖片主點偏移,內(nèi)參矩陣均以像素為單位。

    ff為熱像儀與方向的毫米焦距,對于毫米焦距和像素焦距的轉(zhuǎn)換關(guān)系有:

    式中:ccd_size是熱像儀廠商給出的CCD尺寸,以英寸為單位;image_、image_是圖像、方向的像素數(shù),與成像的分辨率有關(guān)。

    3.3 熱像儀標定

    在工程應(yīng)用中,外參矩陣主要通過機位的計算得出,內(nèi)參矩陣經(jīng)熱像儀標定得出。

    對于熱像儀標定,一般使用張氏標定法[13]。通過Matlab軟件自帶的熱像儀標定程序camera calibrator可對單目熱像儀進行標定??梢姽庀鄼C標定需打印棋盤格,對于紅外熱像儀,將棋盤格用金屬板代替。黑色部分使用隔熱膠帶,白色部分使用金屬板或錫紙。

    在程序中輸入標定板在不同角度下的紅外熱像圖,可得到熱像儀中心與標定板相對位置和熱像儀內(nèi)參、外參、圖片尺寸等參數(shù)。本文主要使用標定得到的熱像儀內(nèi)參矩陣。

    4 溫度修正計算流程

    4.1 溫度修正系數(shù)矩陣

    總體的溫度修正計算流程如圖7所示。

    據(jù)式(7)計算出天頂角,查得該天頂角下的發(fā)射率,然后計算式(6)得到世界坐標系下各測溫點的修正系數(shù)(w)。經(jīng)熱像儀成像后,測溫點由世界坐標w映射到成像平面坐標p=[,]T。w到p的轉(zhuǎn)換關(guān)系由3.3節(jié)標定的熱像儀矩陣計算得出。經(jīng)上述轉(zhuǎn)換關(guān)系可得測溫點在二維平面坐標對應(yīng)的修正系數(shù)(p)。

    熱像儀拍攝的溫度圖像將溫度以矩陣的形式存儲,因此,應(yīng)當將p從離散的點云轉(zhuǎn)換為與熱像圖格式相同的p,j,其中,p,j代表像素塊(,)的坐標。

    對于修正系數(shù)在像素塊范圍內(nèi)梯度較小的點,該像素塊的修正系數(shù)取像素塊內(nèi)的平均值。對在熱像圖上像素位置為(,)的點,該像素的修正系數(shù)為:

    當像素內(nèi)修正系數(shù)梯度變化較大的,可采用該像素內(nèi)修正系數(shù)最大值附近點取平均值的方法,在此像素塊內(nèi)凸顯該形狀非光滑點。

    4.2 成像溫度修正

    對于紅外熱像儀,其所測溫度矩陣為0(p)。修正后的溫度矩陣為所測溫度矩陣與溫度修正系數(shù)矩陣做基本積(哈達瑪積):

    所得溫度(p)即修正后的溫度。

    5 實驗驗證

    本文通過一個熱水桶表面的溫度測量實驗來驗證該算法有效性。首先使用Solidworks軟件進行幾何建模,獲得被測物體幾何信息,經(jīng)Meshlab軟件求得熱像儀可視表面的點云及各點處的法向量,再經(jīng)Matlab計算出各點溫度的修正系數(shù)以及各點對應(yīng)的熱像圖中的坐標,最后將紅外熱像儀直接測得的溫度分布信息導(dǎo)入Matlab進行修正得到水桶表面的修正后的溫度。

    5.1 紅外測溫與熱像儀標定

    圖7 溫度修正流程圖

    圖8 熱像儀機位

    圖9 紅外熱像圖與溫度截線

    使用熱像儀拍攝紅外標定板(如圖10所示),獲取標定板的紅外熱像圖,將熱像圖輸入標定程序,獲得熱像儀內(nèi)參矩陣。

    圖10 紅外熱像儀的標定

    5.2 測溫結(jié)果的修正

    通過Solidworks軟件建模并將得到的水桶表面的云數(shù)據(jù)導(dǎo)入Matlab軟件中,可得到水桶表面與熱像儀在空間中的位置關(guān)系,如圖11所示(黃色點為被測表面的可視點,藍色為被測表面的非可視點)。本實驗中,略去了接近桶口的外徑較大的部分。

    圖11 熱像儀與被測表面的位置關(guān)系

    依據(jù)熱像儀與可視點的空間位置關(guān)系,據(jù)式(7)可計算得各可視點的天頂角,據(jù)圖2可得不同天頂角時,所對應(yīng)的發(fā)射率。然后據(jù)式(6)算得每個可視點的修正系數(shù)。如4.1節(jié)所述,對空間中可視點進行熱像儀矩陣的變換,映射至平面可得出熱像圖中各點對應(yīng)的修正系數(shù),見圖12。進而依據(jù)式(15)得到熱像圖中各點修正后的溫度分布。

    圖12 修正系數(shù)分布

    為了更好地顯示本文方法對紅外熱像儀所測得的溫度分布,本節(jié)在所研究的可視表面高度方向的中心,取與圖9中相同位置的截線進行展示,修正系數(shù)如圖13所示。將所得修正系數(shù)與成像截線溫度相乘得到修正后的截線溫度,修正前后溫度的對比如圖14所示。

    由圖14可見,修正前最大誤差近40℃,真實溫度為335℃,相對誤差約為12.4%;修正后,集中在兩側(cè)大角度測量時的誤差不超過5.0℃,相對誤差小于1.5%。從此可見該方法的有效性。值得說明的是圖13中修正系數(shù)分布曲線并不完全對稱,主要原因是該圖的橫坐標為非常微小的像素點,研究過程中像素點的選取沒有完全對稱所致。更多的平行截線結(jié)果均表明,該修正方法較好地消除了非平表面在紅外熱像測溫時部分因為測量角度過大而引起的發(fā)射率變化所造成的誤差。

    圖13 截線上各點的修正系數(shù)

    圖14 截線上溫度分布修正前和修正后的對比

    6 結(jié)論

    針對使用紅外熱像儀測量非平表面溫度時不同被測量點的定向發(fā)射率的變化對測溫結(jié)果的影響,本文依據(jù)物體表面定量發(fā)射率的變化規(guī)律,給出了各個測量點的修正系數(shù)。同時,通過點云三維建模,利用熱像儀通用幾何成像原理推導(dǎo)出紅外熱像圖像素與實際物體的被測表面之間各個測量點的對應(yīng)關(guān)系,給出了紅外熱像測溫過程中非平表面上的大天頂角被測點的溫度測量結(jié)果的誤差修正方法。實驗證明了該方法的有效性。

    [1] 楊立, 楊楨. 紅外熱成像測溫原理與技術(shù)[M]. 北京: 科學出版社,2016.

    YANG Li, YANG Zhen.[M]. Beijing: Beijing Science Press, 2016.

    [2] 范春利. 幾何形狀導(dǎo)熱反問題方法與應(yīng)用[M]. 北京: 科學出版社, 2015.

    FAN Chunli.:[M]. Beijing: Science Press, 2015.

    [3] 廖盼盼, 張佳民. 紅外測溫精度的影響因素及補償方法的研究[J]. 紅外技術(shù), 2017, 39(2): 173-177.

    LIAO Panpan, ZHANG Jiamin. Research on influence factors for measuring and method of correction in infrared thermometer[J]., 2017, 39(2): 173-177.

    [4] 劉慧開, 楊立. 太陽輻射對紅外熱像儀測溫誤差的影響[J]. 紅外技術(shù), 2002, 24(1): 34-37.

    LIU Huikai, YANG Li. Effect of radiation of the sun on infrared temperature measurement[J]., 2002, 24(1): 34-37.

    [5] 呂游, 楊波, 魏仲慧, 等. 大氣與環(huán)境影響分析的紅外比色測溫方法[J]. 紅外與激光工程, 2015, 44(8): 2321.

    LV You, YANG Bo, WEI Zhonghui, et al. Middle infrared colorimetric temperature measurement method considered influence of atmosphere and ambient[J]., 2015, 44(8): 2321.

    [6] 石東平, 吳超, 李孜軍, 等. 基于反射溫度補償及入射溫度補償?shù)募t外測溫影響分析[J]. 紅外與激光工程, 2015, 44(8): 0204006.

    SHI Dongping, WU Chao, LI Zijun, et al. Analysis of the influence of infrared temperature measurement based on reflected temperature compensation and incidence temperature compensation[J]., 2015, 44(8): 0204006.

    [7] 張健, 楊立, 劉慧開. 環(huán)境高溫物體對紅外熱像儀測溫誤差的影響[J].紅外技術(shù), 2005, 27(5): 419-422.

    ZHANG Jian, YANG Li, LIU Huikai. Effect of environmental object on infrared temperature measurement[J]., 2005, 27(5): 419-422.

    [8] 周志成, 魏旭, 謝天喜, 等. 觀測距離及視角對紅外熱輻射檢測的影響研究[J]. 紅外技術(shù), 2017, 39(1): 86-90.

    ZHOU Zhicheng, WEI Xu, XIE Tianxi, et al. Influence of observation distance and angle of view on the detection accuracy of infrared thermal radiation[J]., 2017, 39(1): 86-90.

    [9] 杜玉璽, 胡振琪, 葛運航, 等. 距離對不同強度熱源紅外測溫影響及補償[J]. 紅外技術(shù), 2019,41(10): 976-981.

    DU Yuxi, HU Zhenqi, GE Yunhang, et al. Distance influence and compensation of infrared temperature measurement with different intensity heat sources[J]., 2019, 41(10): 976-981.

    [10] 范春利, 楊立, 華順芳. 熱探測器溫度對紅外熱像儀測溫的影響[J].紅外技術(shù), 2002, 24(5): 22-24.

    FAN Chunli, YANG Li, HUA Shunfang. Effect of temperature of thermal detector on temperature measurement of uncooled infrared thermography[J]., 2002,24(5):22-24.

    [11] Mikhail Polyanskiy. Refractiveindex.info database [DB/OL]. [2020 -04-01] https://refractiveindex.info/2008-2020.

    [12] 馬頌德, 張正友. 計算機視覺: 計算理論與算法基礎(chǔ)[M]. 北京: 科學出版社, 1998.

    MA S, ZHANG Z.[M]. Beijing: Science Press, 1998.

    [13] ZHANG Z. A flexible new technique for camera calibration[J]., 2000, 22(11): 1330-1334.

    Compensation Method for Temperature Distribution Measured by Infrared Thermography for Non-flat Surfaces

    FU Wanchao,F(xiàn)AN Chunli,YANG Li

    (,,430033,)

    When employing an infrared thermal imager to measure surface temperature, the emissivity of the surface to be measured should be set first and kept constant during the measurement process. However, when the infrared imager is placed in the range of more than 50°of the zenith angle of the points to be measured, the emissivities of the points in this angle will vary significantly; hence, temperature measurement errors will occur, especially for points on non-flat surfaces. In view of the measurement error caused by the variation of emissivity of different points in the measured non-flat surface when using a monocular infrared thermal imager, this paper provides a compensation factor based on the variation rules of the emissivity with the measuring angle. In addition, based on 3D modeling technology, the relationship between the positions of the points in the thermographic image and those in the actual surface is determined. The compensation method of the temperature measurements for a non-flat surface is presented. The feasibility of the method was verified through experiments.

    thermographic non-destructive inspection, infrared imager, temperature measurement error, compensation method, non-flat surface

    TK31

    A

    1001-8891(2021)02-0179-06

    2020-07-29;

    2020-08-20.

    付萬超(1996-),男,碩士生,主要從事紅外無損檢測和傳熱反問題等方面的研究。E-mail:1009202412@qq.com。

    范春利(1978-),男,副教授,主要從事紅外無損檢測和傳熱反問題等方面的研究。E-mail:chlfan@163.com。

    國家自然科學基金項目(51479203)。

    猜你喜歡
    熱像儀頂角發(fā)射率
    一般三棱鏡最大頂角與折射率的關(guān)系
    涼亭中的數(shù)學
    氧氣A(O,O)波段氣輝體發(fā)射率和臨邊輻射強度模擬與分析
    頂角為100°的等腰三角形性質(zhì)的應(yīng)用
    低溫狀態(tài)下的材料法向發(fā)射率測量
    空間噪聲等效溫差測試方法研究
    平房倉空調(diào)設(shè)計及熱像儀應(yīng)用報告
    美國菲力爾公司(FLIR Systems)推出全新T系列紅外熱像儀
    塔克拉瑪干沙漠地表發(fā)射率及分布變化特征
    不透明材料波段法向發(fā)射率在線測量方法
    非洲黑人性xxxx精品又粗又长| 欧美区成人在线视频| 伊人久久精品亚洲午夜| 日本免费一区二区三区高清不卡| 成人国产麻豆网| 亚洲18禁久久av| 亚洲av电影在线观看一区二区三区 | 乱码一卡2卡4卡精品| 特大巨黑吊av在线直播| 国产成人精品一,二区| 免费看日本二区| 色视频www国产| 久热久热在线精品观看| av在线老鸭窝| 久久鲁丝午夜福利片| 国产精品久久久久久久电影| 亚洲成人精品中文字幕电影| 国产av不卡久久| 日本av手机在线免费观看| 成人特级av手机在线观看| 国产亚洲最大av| 免费搜索国产男女视频| 午夜激情欧美在线| 国产一区有黄有色的免费视频 | 亚洲av成人精品一区久久| 中国美白少妇内射xxxbb| 一级黄色大片毛片| 亚洲激情五月婷婷啪啪| 国产一区二区三区av在线| 精品午夜福利在线看| 偷拍熟女少妇极品色| 日本一二三区视频观看| 欧美另类亚洲清纯唯美| 国产成人午夜福利电影在线观看| 国产免费福利视频在线观看| 看非洲黑人一级黄片| 一边亲一边摸免费视频| 国产精品久久久久久av不卡| 国产精品国产三级国产专区5o | 3wmmmm亚洲av在线观看| 男女国产视频网站| 2021少妇久久久久久久久久久| 女的被弄到高潮叫床怎么办| 久久久久久久久中文| 亚洲成人av在线免费| 亚洲av中文av极速乱| 精品一区二区免费观看| 日韩成人av中文字幕在线观看| 亚洲av免费高清在线观看| 国产一级毛片在线| 日韩在线高清观看一区二区三区| 日本-黄色视频高清免费观看| 国产 一区精品| 亚洲av成人av| 欧美色视频一区免费| 成年av动漫网址| 97超碰精品成人国产| 全区人妻精品视频| 亚洲美女视频黄频| 久久精品久久精品一区二区三区| 人妻少妇偷人精品九色| 国产乱人偷精品视频| 国产伦一二天堂av在线观看| 国产精品综合久久久久久久免费| 最近最新中文字幕免费大全7| 91在线精品国自产拍蜜月| 最近的中文字幕免费完整| 国产一区有黄有色的免费视频 | 国产综合懂色| 亚洲色图av天堂| 日本五十路高清| 久久人人爽人人爽人人片va| 2021少妇久久久久久久久久久| 99久久人妻综合| 国产午夜精品论理片| 国产免费一级a男人的天堂| 午夜a级毛片| 国产精品女同一区二区软件| 乱系列少妇在线播放| 真实男女啪啪啪动态图| 免费av不卡在线播放| 欧美日韩在线观看h| 成年女人永久免费观看视频| 亚洲欧美中文字幕日韩二区| 成人综合一区亚洲| 在线天堂最新版资源| www日本黄色视频网| 淫秽高清视频在线观看| 国产老妇伦熟女老妇高清| 亚洲天堂国产精品一区在线| 一个人看的www免费观看视频| 久久99精品国语久久久| 九色成人免费人妻av| 日韩欧美精品免费久久| 色尼玛亚洲综合影院| 国产精品99久久久久久久久| 黄片wwwwww| 国产精品女同一区二区软件| 国产乱来视频区| 看黄色毛片网站| 国产毛片a区久久久久| av在线观看视频网站免费| 国产精品日韩av在线免费观看| av卡一久久| 午夜爱爱视频在线播放| 日韩视频在线欧美| 国产女主播在线喷水免费视频网站 | 久久欧美精品欧美久久欧美| 男女边吃奶边做爰视频| 成人性生交大片免费视频hd| 国产伦精品一区二区三区四那| 身体一侧抽搐| 亚洲内射少妇av| 中文字幕精品亚洲无线码一区| 美女被艹到高潮喷水动态| 两个人的视频大全免费| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 搞女人的毛片| 亚洲精品aⅴ在线观看| 有码 亚洲区| 一本一本综合久久| 中文字幕久久专区| 成人毛片a级毛片在线播放| 2021天堂中文幕一二区在线观| 麻豆久久精品国产亚洲av| 在线观看av片永久免费下载| 亚洲欧美日韩无卡精品| 美女国产视频在线观看| 最近视频中文字幕2019在线8| 少妇裸体淫交视频免费看高清| 日本午夜av视频| 久久99热这里只有精品18| av又黄又爽大尺度在线免费看 | 寂寞人妻少妇视频99o| 看免费成人av毛片| 青春草视频在线免费观看| 又爽又黄a免费视频| 中文字幕免费在线视频6| 国产黄色视频一区二区在线观看 | 男人狂女人下面高潮的视频| av视频在线观看入口| 69av精品久久久久久| 最后的刺客免费高清国语| 亚洲欧美中文字幕日韩二区| 久久久久九九精品影院| 午夜福利视频1000在线观看| 91精品伊人久久大香线蕉| 99久久无色码亚洲精品果冻| 国产视频首页在线观看| 国产精品爽爽va在线观看网站| 在线观看美女被高潮喷水网站| 嫩草影院新地址| 高清日韩中文字幕在线| 亚洲熟妇中文字幕五十中出| 亚洲av免费在线观看| 亚州av有码| 日本猛色少妇xxxxx猛交久久| 美女国产视频在线观看| 美女大奶头视频| 成人毛片a级毛片在线播放| 中文欧美无线码| 国产一区二区三区av在线| 中文精品一卡2卡3卡4更新| 综合色av麻豆| 禁无遮挡网站| 亚洲美女搞黄在线观看| 久久这里有精品视频免费| av天堂中文字幕网| 2021天堂中文幕一二区在线观| 久久99热这里只有精品18| 亚洲国产欧美人成| 亚洲av一区综合| 天天躁日日操中文字幕| 午夜爱爱视频在线播放| 久久久久久久国产电影| 99热精品在线国产| 少妇熟女aⅴ在线视频| 日韩人妻高清精品专区| 亚洲av成人精品一二三区| 国产精品久久久久久精品电影| 免费大片18禁| 一本久久精品| 男人狂女人下面高潮的视频| 亚洲三级黄色毛片| 黄色配什么色好看| 波多野结衣高清无吗| 国产精品一及| 午夜亚洲福利在线播放| 中文在线观看免费www的网站| 欧美区成人在线视频| 国产欧美另类精品又又久久亚洲欧美| 久久久色成人| 色5月婷婷丁香| av播播在线观看一区| 日本av手机在线免费观看| 色尼玛亚洲综合影院| 国产熟女欧美一区二区| 亚洲精品国产成人久久av| 久久久国产成人精品二区| 国产精品99久久久久久久久| 最新中文字幕久久久久| 97热精品久久久久久| 国产亚洲av片在线观看秒播厂 | 色5月婷婷丁香| 免费av毛片视频| 国产亚洲精品久久久com| 国产精品野战在线观看| 日本五十路高清| 嫩草影院入口| 夜夜爽夜夜爽视频| 村上凉子中文字幕在线| 国产成人福利小说| 国产亚洲一区二区精品| 国产又色又爽无遮挡免| 亚洲电影在线观看av| 免费观看a级毛片全部| 少妇裸体淫交视频免费看高清| 欧美极品一区二区三区四区| 欧美日本视频| 哪个播放器可以免费观看大片| 成人午夜高清在线视频| 最近最新中文字幕大全电影3| 1000部很黄的大片| 欧美一区二区亚洲| 直男gayav资源| 插逼视频在线观看| 日本一本二区三区精品| 国产淫片久久久久久久久| 中文字幕久久专区| 国产精品伦人一区二区| 99热这里只有是精品在线观看| 久久久久性生活片| 亚洲人成网站在线播| 国产色爽女视频免费观看| 国产精品野战在线观看| 国产免费男女视频| 国内精品美女久久久久久| 亚洲国产精品sss在线观看| 日韩欧美三级三区| 一个人免费在线观看电影| 欧美性猛交╳xxx乱大交人| 国产精华一区二区三区| 晚上一个人看的免费电影| 小蜜桃在线观看免费完整版高清| 亚洲第一区二区三区不卡| 日韩精品有码人妻一区| 长腿黑丝高跟| 丰满少妇做爰视频| 欧美另类亚洲清纯唯美| 久久久久久久久中文| 免费无遮挡裸体视频| 人体艺术视频欧美日本| 亚洲欧美日韩无卡精品| 色吧在线观看| 亚洲av一区综合| 国产片特级美女逼逼视频| 久久99热6这里只有精品| 天堂影院成人在线观看| 日韩一区二区视频免费看| 亚洲人成网站高清观看| 国产成人freesex在线| 哪个播放器可以免费观看大片| 丰满乱子伦码专区| 亚洲美女搞黄在线观看| 日本免费在线观看一区| 色综合站精品国产| 国产v大片淫在线免费观看| 国产精品一区二区在线观看99 | 能在线免费观看的黄片| 精品酒店卫生间| 中文精品一卡2卡3卡4更新| 亚洲国产色片| 99久久九九国产精品国产免费| 日韩三级伦理在线观看| 建设人人有责人人尽责人人享有的 | 天美传媒精品一区二区| 丝袜喷水一区| 日本-黄色视频高清免费观看| 亚洲中文字幕一区二区三区有码在线看| 成人毛片a级毛片在线播放| 成人毛片60女人毛片免费| 国内精品美女久久久久久| 久久久精品94久久精品| 久久久久网色| 欧美激情久久久久久爽电影| 一边亲一边摸免费视频| 国产精品伦人一区二区| 精品人妻偷拍中文字幕| 好男人视频免费观看在线| 久久人妻av系列| 91午夜精品亚洲一区二区三区| 一二三四中文在线观看免费高清| 最近最新中文字幕大全电影3| 午夜日本视频在线| 好男人视频免费观看在线| 亚洲人成网站在线观看播放| 日本wwww免费看| 2021少妇久久久久久久久久久| av专区在线播放| 哪个播放器可以免费观看大片| 男插女下体视频免费在线播放| 午夜免费男女啪啪视频观看| 亚洲欧美成人综合另类久久久 | 可以在线观看毛片的网站| 中文欧美无线码| 一级毛片我不卡| 啦啦啦观看免费观看视频高清| 国产午夜精品久久久久久一区二区三区| 性色avwww在线观看| 午夜视频国产福利| 在线观看一区二区三区| 51国产日韩欧美| 白带黄色成豆腐渣| 欧美性猛交╳xxx乱大交人| 国产精品永久免费网站| av在线观看视频网站免费| 91精品伊人久久大香线蕉| 久久人人爽人人片av| 中文字幕亚洲精品专区| 国产精品野战在线观看| 夜夜看夜夜爽夜夜摸| 国产成人精品婷婷| 国产精华一区二区三区| 国产精品国产三级国产专区5o | h日本视频在线播放| 日本与韩国留学比较| 91久久精品电影网| av国产免费在线观看| 有码 亚洲区| 免费观看的影片在线观看| 欧美成人免费av一区二区三区| 天天一区二区日本电影三级| 日本与韩国留学比较| 嫩草影院新地址| 欧美3d第一页| 不卡视频在线观看欧美| 国产 一区精品| 国产视频首页在线观看| 日本黄色片子视频| 国产成人一区二区在线| 亚洲欧美日韩无卡精品| 视频中文字幕在线观看| 超碰97精品在线观看| 免费看光身美女| 看非洲黑人一级黄片| 男女视频在线观看网站免费| 国产乱来视频区| 欧美+日韩+精品| 神马国产精品三级电影在线观看| 亚洲国产欧美人成| 日韩 亚洲 欧美在线| 久久久久久久久中文| 国产免费一级a男人的天堂| 日韩大片免费观看网站 | 搡老妇女老女人老熟妇| 国产单亲对白刺激| 人妻制服诱惑在线中文字幕| eeuss影院久久| 嫩草影院入口| 又爽又黄a免费视频| 人体艺术视频欧美日本| 国内精品一区二区在线观看| 婷婷六月久久综合丁香| 禁无遮挡网站| 国产大屁股一区二区在线视频| 色网站视频免费| 99久久中文字幕三级久久日本| 午夜老司机福利剧场| 亚洲精品乱码久久久v下载方式| 三级男女做爰猛烈吃奶摸视频| 有码 亚洲区| videos熟女内射| 久久欧美精品欧美久久欧美| 国产精品久久久久久av不卡| 国产伦精品一区二区三区四那| 最新中文字幕久久久久| 插阴视频在线观看视频| 国产成人午夜福利电影在线观看| 亚洲精品乱码久久久久久按摩| 99热网站在线观看| 91在线精品国自产拍蜜月| 一本一本综合久久| 欧美成人午夜免费资源| 国产黄片美女视频| 亚洲美女视频黄频| 超碰av人人做人人爽久久| 2021少妇久久久久久久久久久| 亚洲av成人精品一区久久| 亚洲中文字幕日韩| 欧美一区二区国产精品久久精品| 最近的中文字幕免费完整| 成人一区二区视频在线观看| 日本-黄色视频高清免费观看| 日韩欧美三级三区| 午夜精品在线福利| 少妇的逼水好多| av播播在线观看一区| 欧美成人a在线观看| 中文字幕人妻熟人妻熟丝袜美| 国产黄片美女视频| 国内揄拍国产精品人妻在线| 亚洲av熟女| 久久这里有精品视频免费| 色综合亚洲欧美另类图片| 干丝袜人妻中文字幕| 国产高潮美女av| 色播亚洲综合网| 看黄色毛片网站| 精品午夜福利在线看| 99久久成人亚洲精品观看| 亚洲乱码一区二区免费版| 1024手机看黄色片| av福利片在线观看| 日本爱情动作片www.在线观看| 日韩视频在线欧美| 长腿黑丝高跟| 免费观看精品视频网站| 听说在线观看完整版免费高清| 亚洲av成人精品一区久久| 黄片无遮挡物在线观看| 欧美xxxx黑人xx丫x性爽| 老女人水多毛片| 成人高潮视频无遮挡免费网站| 在线观看一区二区三区| 国产高清有码在线观看视频| 婷婷色av中文字幕| 久久精品国产亚洲av涩爱| 亚洲欧美精品专区久久| 亚洲欧美日韩卡通动漫| 精品人妻一区二区三区麻豆| 波多野结衣高清无吗| av女优亚洲男人天堂| 欧美三级亚洲精品| 午夜精品在线福利| 大话2 男鬼变身卡| 男人的好看免费观看在线视频| 国产女主播在线喷水免费视频网站 | 人妻夜夜爽99麻豆av| 亚洲欧美日韩高清专用| 国内精品美女久久久久久| 大香蕉97超碰在线| 九九在线视频观看精品| 成人无遮挡网站| 黄色日韩在线| 午夜老司机福利剧场| 欧美区成人在线视频| 91久久精品国产一区二区成人| 亚洲国产精品合色在线| 少妇猛男粗大的猛烈进出视频 | 青春草国产在线视频| 国产三级中文精品| 淫秽高清视频在线观看| 青春草视频在线免费观看| 免费看a级黄色片| 久久草成人影院| 亚洲精品亚洲一区二区| 午夜激情欧美在线| 国产精品综合久久久久久久免费| 一区二区三区免费毛片| 国产高清不卡午夜福利| 久久精品久久精品一区二区三区| 欧美97在线视频| 99久久精品热视频| 中文字幕人妻熟人妻熟丝袜美| 中文欧美无线码| 国产老妇伦熟女老妇高清| 国产精品久久久久久久电影| 日本一二三区视频观看| 大话2 男鬼变身卡| 一本久久精品| 精品不卡国产一区二区三区| 最后的刺客免费高清国语| 一级黄色大片毛片| 午夜老司机福利剧场| 亚洲中文字幕日韩| 老司机影院成人| 你懂的网址亚洲精品在线观看 | 又黄又爽又刺激的免费视频.| ponron亚洲| 丰满少妇做爰视频| 网址你懂的国产日韩在线| 国产精品久久久久久av不卡| 色网站视频免费| 禁无遮挡网站| 欧美一区二区精品小视频在线| 国产黄色视频一区二区在线观看 | 99热这里只有是精品在线观看| 久久99精品国语久久久| 女人被狂操c到高潮| av在线蜜桃| 精品免费久久久久久久清纯| 免费无遮挡裸体视频| 性插视频无遮挡在线免费观看| 在线a可以看的网站| 美女内射精品一级片tv| 欧美最新免费一区二区三区| 欧美成人一区二区免费高清观看| 国产精品乱码一区二三区的特点| 久久久久久久午夜电影| 国产激情偷乱视频一区二区| 欧美性猛交黑人性爽| 老师上课跳d突然被开到最大视频| 高清日韩中文字幕在线| 精品久久久久久久久久久久久| 亚洲怡红院男人天堂| 欧美成人午夜免费资源| 久久99热这里只有精品18| 午夜久久久久精精品| 久久婷婷人人爽人人干人人爱| 国产亚洲av片在线观看秒播厂 | 少妇熟女aⅴ在线视频| 欧美xxxx黑人xx丫x性爽| 一区二区三区高清视频在线| 午夜福利成人在线免费观看| 少妇的逼好多水| 少妇熟女aⅴ在线视频| 亚洲av免费高清在线观看| 只有这里有精品99| 国产在线一区二区三区精 | 女人十人毛片免费观看3o分钟| 色5月婷婷丁香| 日本猛色少妇xxxxx猛交久久| 国产乱人偷精品视频| 亚洲欧美成人综合另类久久久 | 国产欧美日韩精品一区二区| 日韩欧美三级三区| 女的被弄到高潮叫床怎么办| 免费看a级黄色片| 免费在线观看成人毛片| 91久久精品电影网| 中文乱码字字幕精品一区二区三区 | 亚洲av中文字字幕乱码综合| 丰满乱子伦码专区| 人人妻人人看人人澡| 综合色av麻豆| 少妇的逼水好多| 亚洲最大成人av| 99久久精品一区二区三区| 亚洲av电影在线观看一区二区三区 | 国产精品国产三级国产av玫瑰| 亚洲av成人精品一区久久| 国产男人的电影天堂91| 成人美女网站在线观看视频| 亚洲一级一片aⅴ在线观看| 亚洲不卡免费看| 国产精品久久久久久精品电影| 丝袜喷水一区| 国产av不卡久久| a级一级毛片免费在线观看| 日韩强制内射视频| 少妇的逼水好多| 国产精品女同一区二区软件| av在线亚洲专区| 中国国产av一级| 国产伦一二天堂av在线观看| 99在线人妻在线中文字幕| 中文在线观看免费www的网站| 亚洲国产精品专区欧美| 蜜桃久久精品国产亚洲av| 免费观看a级毛片全部| 99热6这里只有精品| 91av网一区二区| 18禁动态无遮挡网站| 蜜臀久久99精品久久宅男| 亚洲va在线va天堂va国产| 熟女电影av网| 精品欧美国产一区二区三| 身体一侧抽搐| 久久久色成人| 两个人的视频大全免费| 国产美女午夜福利| 亚洲精品,欧美精品| 色视频www国产| 日本五十路高清| 亚洲精品456在线播放app| 精品人妻视频免费看| 内地一区二区视频在线| 国产成人91sexporn| 一级爰片在线观看| 久久精品综合一区二区三区| 色5月婷婷丁香| 99久久无色码亚洲精品果冻| 亚洲电影在线观看av| 国产精品久久久久久精品电影小说 | 国产精品爽爽va在线观看网站| 欧美区成人在线视频| 一级av片app| 男女边吃奶边做爰视频| 最近中文字幕2019免费版| 99热这里只有是精品在线观看| 国产高清视频在线观看网站| 黄色一级大片看看| 国产黄色视频一区二区在线观看 | 人人妻人人看人人澡| 一夜夜www| 九九在线视频观看精品| 精品人妻熟女av久视频| 亚洲成人精品中文字幕电影| 成年av动漫网址| 亚洲成人精品中文字幕电影| 国产久久久一区二区三区| 99久久精品热视频| 久久久a久久爽久久v久久| 床上黄色一级片| 欧美高清成人免费视频www| 日本免费a在线| 夫妻性生交免费视频一级片| 日韩高清综合在线| 成人毛片a级毛片在线播放| 麻豆乱淫一区二区| www日本黄色视频网| videos熟女内射| 精品久久国产蜜桃| 国产精品不卡视频一区二区| av视频在线观看入口|