• <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ā)射率在線測量方法
    婷婷丁香在线五月| 国产又色又爽无遮挡免费看| 视频区欧美日本亚洲| 午夜精品在线福利| 亚洲人成网站高清观看| 精品一区二区三区视频在线观看免费| 国产精品一及| 日本 欧美在线| 久久久久国产一级毛片高清牌| 日日爽夜夜爽网站| 18禁观看日本| 哪里可以看免费的av片| 成年人黄色毛片网站| 午夜福利成人在线免费观看| 又黄又爽又免费观看的视频| 国产一区在线观看成人免费| 国产黄色小视频在线观看| 一进一出抽搐动态| 成人欧美大片| 18禁观看日本| 俄罗斯特黄特色一大片| 午夜激情福利司机影院| 99国产精品99久久久久| 波多野结衣高清无吗| 波多野结衣高清作品| 黄片大片在线免费观看| 中文字幕最新亚洲高清| 亚洲专区国产一区二区| 欧美午夜高清在线| 禁无遮挡网站| 欧美最黄视频在线播放免费| 亚洲 国产 在线| 黄色视频不卡| 女警被强在线播放| 搡老熟女国产l中国老女人| 97碰自拍视频| 国产精品av久久久久免费| 99久久国产精品久久久| 亚洲国产日韩欧美精品在线观看 | 国产精品一区二区免费欧美| 成人特级黄色片久久久久久久| 国产高清激情床上av| 69av精品久久久久久| 特大巨黑吊av在线直播| 午夜激情av网站| 亚洲国产看品久久| 99riav亚洲国产免费| 国内久久婷婷六月综合欲色啪| 成年版毛片免费区| 国内揄拍国产精品人妻在线| 18禁裸乳无遮挡免费网站照片| 亚洲第一电影网av| 欧美日韩黄片免| 深夜精品福利| 欧美人与性动交α欧美精品济南到| 日本一本二区三区精品| 亚洲国产高清在线一区二区三| 色综合亚洲欧美另类图片| 日本 欧美在线| 高清在线国产一区| 欧美日韩黄片免| 亚洲精品美女久久久久99蜜臀| 久久性视频一级片| 成人一区二区视频在线观看| 日韩精品青青久久久久久| 九九热线精品视视频播放| 热99re8久久精品国产| 一区二区三区国产精品乱码| 精品电影一区二区在线| 国产三级黄色录像| 少妇的丰满在线观看| 欧美 亚洲 国产 日韩一| av国产免费在线观看| 亚洲一码二码三码区别大吗| 精品第一国产精品| 午夜福利高清视频| 在线观看美女被高潮喷水网站 | 亚洲一区二区三区色噜噜| 亚洲av成人av| 免费在线观看亚洲国产| 不卡av一区二区三区| 国产成人系列免费观看| 亚洲国产欧美一区二区综合| 国产成人av激情在线播放| 欧美三级亚洲精品| 国产一区二区三区在线臀色熟女| 校园春色视频在线观看| 国产午夜福利久久久久久| 男女床上黄色一级片免费看| 日韩有码中文字幕| 欧美一区二区精品小视频在线| 国产成年人精品一区二区| 久久精品成人免费网站| 最新美女视频免费是黄的| 亚洲国产高清在线一区二区三| 波多野结衣高清作品| 精品熟女少妇八av免费久了| 成人永久免费在线观看视频| 1024香蕉在线观看| 成人av在线播放网站| 国产激情偷乱视频一区二区| 亚洲欧美日韩无卡精品| 精品久久久久久久末码| 怎么达到女性高潮| 男人舔奶头视频| 国内毛片毛片毛片毛片毛片| 日本在线视频免费播放| 久久久久久免费高清国产稀缺| 99精品欧美一区二区三区四区| 成人手机av| 精品国内亚洲2022精品成人| 欧美av亚洲av综合av国产av| 中国美女看黄片| 国产又黄又爽又无遮挡在线| 亚洲成人久久性| 亚洲无线在线观看| 人妻久久中文字幕网| 久久精品亚洲精品国产色婷小说| 国产欧美日韩一区二区三| 成人手机av| 国产激情偷乱视频一区二区| 一a级毛片在线观看| 久久精品91无色码中文字幕| 日韩大尺度精品在线看网址| 日韩欧美精品v在线| 久久精品国产亚洲av高清一级| 欧洲精品卡2卡3卡4卡5卡区| 国产精品综合久久久久久久免费| 亚洲av第一区精品v没综合| 国产久久久一区二区三区| 久久精品91蜜桃| 久久久久久国产a免费观看| 夜夜夜夜夜久久久久| 女人爽到高潮嗷嗷叫在线视频| 亚洲美女黄片视频| 给我免费播放毛片高清在线观看| 女人爽到高潮嗷嗷叫在线视频| 18禁黄网站禁片免费观看直播| 老鸭窝网址在线观看| 久久人妻av系列| 一区二区三区高清视频在线| 51午夜福利影视在线观看| www.精华液| 国产亚洲精品久久久久久毛片| 亚洲18禁久久av| 中文字幕精品亚洲无线码一区| 欧美乱色亚洲激情| 免费在线观看完整版高清| 女警被强在线播放| 在线观看66精品国产| 黄片小视频在线播放| 国产真实乱freesex| 国产熟女xx| 婷婷六月久久综合丁香| 99久久国产精品久久久| 女同久久另类99精品国产91| 亚洲狠狠婷婷综合久久图片| 欧美大码av| 亚洲 国产 在线| 成人午夜高清在线视频| 啦啦啦韩国在线观看视频| 亚洲真实伦在线观看| 久久精品国产亚洲av高清一级| 国产aⅴ精品一区二区三区波| 亚洲 国产 在线| 岛国在线免费视频观看| 日日爽夜夜爽网站| 婷婷六月久久综合丁香| 两人在一起打扑克的视频| 手机成人av网站| 波多野结衣高清作品| 久久久久久亚洲精品国产蜜桃av| 国产午夜福利久久久久久| 国产精品自产拍在线观看55亚洲| 欧美久久黑人一区二区| 欧美又色又爽又黄视频| 99riav亚洲国产免费| 精品久久久久久久人妻蜜臀av| 村上凉子中文字幕在线| 麻豆国产av国片精品| 国产午夜精品论理片| 久久精品国产亚洲av香蕉五月| 精品欧美国产一区二区三| 中文资源天堂在线| 美女黄网站色视频| 99久久综合精品五月天人人| 好男人电影高清在线观看| 亚洲欧美一区二区三区黑人| 成人一区二区视频在线观看| 啦啦啦免费观看视频1| 成人手机av| 久久精品亚洲精品国产色婷小说| 19禁男女啪啪无遮挡网站| 性欧美人与动物交配| 国产精品久久久久久久电影 | 亚洲精品美女久久久久99蜜臀| 最新在线观看一区二区三区| 两个人的视频大全免费| a级毛片在线看网站| 国产激情偷乱视频一区二区| 亚洲,欧美精品.| 精品无人区乱码1区二区| 色老头精品视频在线观看| 亚洲精品色激情综合| 久久久久久亚洲精品国产蜜桃av| 亚洲成人久久性| 最近最新中文字幕大全免费视频| 在线国产一区二区在线| 亚洲精品色激情综合| 成年女人毛片免费观看观看9| 午夜两性在线视频| 久久精品国产亚洲av高清一级| 高潮久久久久久久久久久不卡| 在线视频色国产色| 女警被强在线播放| 精品一区二区三区视频在线观看免费| 免费观看人在逋| 50天的宝宝边吃奶边哭怎么回事| 男女下面进入的视频免费午夜| 99精品久久久久人妻精品| 欧美绝顶高潮抽搐喷水| 亚洲熟妇中文字幕五十中出| 伦理电影免费视频| 在线观看一区二区三区| 国产亚洲欧美98| 国产亚洲av高清不卡| 极品教师在线免费播放| 国内精品久久久久精免费| www日本黄色视频网| 别揉我奶头~嗯~啊~动态视频| 欧美日韩瑟瑟在线播放| av在线天堂中文字幕| 亚洲欧洲精品一区二区精品久久久| 成人国产一区最新在线观看| 后天国语完整版免费观看| 超碰成人久久| 欧美午夜高清在线| 桃红色精品国产亚洲av| 性欧美人与动物交配| 搡老熟女国产l中国老女人| 又紧又爽又黄一区二区| 国产精品一及| 在线观看一区二区三区| а√天堂www在线а√下载| 精品日产1卡2卡| 欧美日本视频| 亚洲成人精品中文字幕电影| 舔av片在线| or卡值多少钱| 真人一进一出gif抽搐免费| 成人18禁在线播放| 久久久久免费精品人妻一区二区| 欧美三级亚洲精品| 久久精品影院6| aaaaa片日本免费| 午夜久久久久精精品| 国产片内射在线| av免费在线观看网站| 国产私拍福利视频在线观看| 俺也久久电影网| 国产探花在线观看一区二区| 色老头精品视频在线观看| 久久久久国产精品人妻aⅴ院| 无遮挡黄片免费观看| 少妇裸体淫交视频免费看高清 | 午夜福利成人在线免费观看| 一本综合久久免费| 成人av一区二区三区在线看| 丰满人妻一区二区三区视频av | 久久 成人 亚洲| 午夜福利高清视频| 18禁国产床啪视频网站| 亚洲熟妇中文字幕五十中出| 黄色视频不卡| 国产成人影院久久av| bbb黄色大片| 精品久久久久久久久久久久久| 日韩欧美国产一区二区入口| 午夜福利在线观看吧| 国产av麻豆久久久久久久| 国产主播在线观看一区二区| 欧美乱妇无乱码| 欧美乱码精品一区二区三区| 国产三级黄色录像| 可以在线观看毛片的网站| 国产私拍福利视频在线观看| 国产av一区在线观看免费| 欧美av亚洲av综合av国产av| 欧美丝袜亚洲另类 | 99在线人妻在线中文字幕| 久久久久精品国产欧美久久久| 天天躁狠狠躁夜夜躁狠狠躁| xxxwww97欧美| 99国产极品粉嫩在线观看| 麻豆成人午夜福利视频| 老司机靠b影院| 色播亚洲综合网| 亚洲国产日韩欧美精品在线观看 | 香蕉av资源在线| www国产在线视频色| 嫩草影视91久久| 麻豆国产av国片精品| 国产精品亚洲一级av第二区| a级毛片在线看网站| 亚洲最大成人中文| 两个人视频免费观看高清| 久久久久免费精品人妻一区二区| 国产精品精品国产色婷婷| 不卡av一区二区三区| 99精品欧美一区二区三区四区| 成人18禁在线播放| 99re在线观看精品视频| 精品午夜福利视频在线观看一区| 午夜福利18| 国内精品久久久久久久电影| 黄色 视频免费看| 久久精品成人免费网站| 国产精品亚洲一级av第二区| 国产成人影院久久av| 国语自产精品视频在线第100页| 51午夜福利影视在线观看| 国产熟女xx| 90打野战视频偷拍视频| 日韩欧美在线乱码| 嫩草影院精品99| 男男h啪啪无遮挡| 天堂动漫精品| 亚洲欧美日韩高清在线视频| 十八禁网站免费在线| 黑人欧美特级aaaaaa片| 亚洲av美国av| 亚洲在线自拍视频| www.999成人在线观看| 九色国产91popny在线| 色综合亚洲欧美另类图片| 精品久久久久久成人av| 午夜福利欧美成人| 亚洲国产中文字幕在线视频| 久久精品夜夜夜夜夜久久蜜豆 | 日韩精品青青久久久久久| 在线观看免费午夜福利视频| 嫩草影视91久久| 女警被强在线播放| 日韩精品青青久久久久久| 可以在线观看的亚洲视频| 悠悠久久av| 后天国语完整版免费观看| 999久久久精品免费观看国产| 十八禁人妻一区二区| 国产精品1区2区在线观看.| 国产黄色小视频在线观看| 亚洲专区字幕在线| 亚洲av熟女| 日韩精品中文字幕看吧| 欧美zozozo另类| 欧美激情久久久久久爽电影| 亚洲国产精品成人综合色| 国产精品国产高清国产av| 18禁美女被吸乳视频| 成人永久免费在线观看视频| 狠狠狠狠99中文字幕| 国产亚洲精品一区二区www| av在线天堂中文字幕| 欧美高清成人免费视频www| 两个人免费观看高清视频| 午夜老司机福利片| 午夜影院日韩av| 真人一进一出gif抽搐免费| 亚洲最大成人中文| 成人午夜高清在线视频| av免费在线观看网站| 久久天躁狠狠躁夜夜2o2o| 变态另类成人亚洲欧美熟女| 男插女下体视频免费在线播放| 一区二区三区激情视频| 99精品久久久久人妻精品| 黄色 视频免费看| 别揉我奶头~嗯~啊~动态视频| 日韩欧美 国产精品| 欧美日韩福利视频一区二区| 亚洲精品在线观看二区| 少妇的丰满在线观看| 18禁黄网站禁片免费观看直播| 欧美成人免费av一区二区三区| 国产又色又爽无遮挡免费看| 精品久久久久久久久久久久久| 香蕉av资源在线| 欧美久久黑人一区二区| 女生性感内裤真人,穿戴方法视频| 免费一级毛片在线播放高清视频| 天天添夜夜摸| 搞女人的毛片| 亚洲精华国产精华精| 精品少妇一区二区三区视频日本电影| 18禁裸乳无遮挡免费网站照片| 日韩中文字幕欧美一区二区| 久久久国产成人免费| 亚洲专区中文字幕在线| 两人在一起打扑克的视频| 天天躁狠狠躁夜夜躁狠狠躁| 国产成+人综合+亚洲专区| 久久久精品欧美日韩精品| 国产亚洲精品一区二区www| 黄频高清免费视频| 国产精品免费视频内射| 一本一本综合久久| 97碰自拍视频| 久久久久国产一级毛片高清牌| 国产精品香港三级国产av潘金莲| 亚洲欧洲精品一区二区精品久久久| 国产精品,欧美在线| 欧美一级a爱片免费观看看 | 777久久人妻少妇嫩草av网站| 麻豆av在线久日| 精品国产亚洲在线| 18禁黄网站禁片午夜丰满| 18禁观看日本| 欧美中文日本在线观看视频| 一区福利在线观看| 国产精品 欧美亚洲| 免费电影在线观看免费观看| 国产精品电影一区二区三区| 欧美日韩亚洲综合一区二区三区_| 一本大道久久a久久精品| 国产成人精品无人区| 三级男女做爰猛烈吃奶摸视频| 成人手机av| 亚洲专区中文字幕在线| www国产在线视频色| 欧美黑人精品巨大| 久久久精品大字幕| 亚洲av日韩精品久久久久久密| 真人做人爱边吃奶动态| 搡老妇女老女人老熟妇| 小说图片视频综合网站| 最近在线观看免费完整版| 亚洲专区中文字幕在线| 中出人妻视频一区二区| 国产精品一及| 99国产精品99久久久久| 亚洲七黄色美女视频| 一级黄色大片毛片| 青草久久国产| av天堂在线播放| 久久久久国产一级毛片高清牌| 久久这里只有精品中国| 精品久久久久久久久久免费视频| 九九热线精品视视频播放| 老司机深夜福利视频在线观看| 日日摸夜夜添夜夜添小说| 露出奶头的视频| 亚洲最大成人中文| 亚洲五月婷婷丁香| 精品电影一区二区在线| 国产v大片淫在线免费观看| 色综合亚洲欧美另类图片| 日韩大尺度精品在线看网址| 丝袜美腿诱惑在线| 国产精品1区2区在线观看.| 久久午夜亚洲精品久久| 国产三级中文精品| 热99re8久久精品国产| 亚洲成人久久性| 欧美日韩福利视频一区二区| 丁香欧美五月| 国产v大片淫在线免费观看| 亚洲乱码一区二区免费版| 最新在线观看一区二区三区| 国产精品久久久久久精品电影| 精品一区二区三区av网在线观看| 国产1区2区3区精品| 亚洲狠狠婷婷综合久久图片| 国产黄片美女视频| 99精品久久久久人妻精品| 午夜福利18| 国产片内射在线| 久久香蕉国产精品| 欧美高清成人免费视频www| 欧美zozozo另类| 动漫黄色视频在线观看| e午夜精品久久久久久久| 亚洲美女视频黄频| 亚洲 欧美一区二区三区| 国产精品av久久久久免费| 国产精品99久久99久久久不卡| 精品欧美国产一区二区三| 最好的美女福利视频网| 中文字幕最新亚洲高清| 一个人免费在线观看的高清视频| 国语自产精品视频在线第100页| 好男人在线观看高清免费视频| 亚洲真实伦在线观看| 亚洲国产精品久久男人天堂| 久久香蕉国产精品| 婷婷亚洲欧美| 国产真实乱freesex| 97人妻精品一区二区三区麻豆| 观看免费一级毛片| 99国产综合亚洲精品| 变态另类成人亚洲欧美熟女| 国产成人av教育| 日韩成人在线观看一区二区三区| 正在播放国产对白刺激| 亚洲人成77777在线视频| 亚洲色图av天堂| 午夜精品久久久久久毛片777| 亚洲自拍偷在线| 成人18禁高潮啪啪吃奶动态图| 国产亚洲欧美在线一区二区| 成年版毛片免费区| 免费在线观看黄色视频的| 巨乳人妻的诱惑在线观看| 久久久精品大字幕| 18禁美女被吸乳视频| 妹子高潮喷水视频| 精品久久蜜臀av无| 岛国视频午夜一区免费看| 一本综合久久免费| 给我免费播放毛片高清在线观看| 日韩欧美精品v在线| 国产精品一区二区免费欧美| 少妇裸体淫交视频免费看高清 | 久久亚洲真实| 一区二区三区高清视频在线| 午夜福利成人在线免费观看| 19禁男女啪啪无遮挡网站| 又黄又爽又免费观看的视频| 成人av在线播放网站| 亚洲国产看品久久| 国产av麻豆久久久久久久| 色播亚洲综合网| 精品电影一区二区在线| 日本三级黄在线观看| 国产精品日韩av在线免费观看| 日韩精品中文字幕看吧| 免费在线观看黄色视频的| 日本 av在线| 99国产精品99久久久久| 久久精品91无色码中文字幕| 神马国产精品三级电影在线观看 | 国产一区二区激情短视频| 亚洲午夜精品一区,二区,三区| 欧美人与性动交α欧美精品济南到| www日本黄色视频网| 亚洲 欧美一区二区三区| 亚洲aⅴ乱码一区二区在线播放 | 人妻久久中文字幕网| 黄色视频不卡| av免费在线观看网站| x7x7x7水蜜桃| 丰满人妻一区二区三区视频av | 色综合欧美亚洲国产小说| 99热只有精品国产| 天天一区二区日本电影三级| 在线观看免费午夜福利视频| 婷婷六月久久综合丁香| 91在线观看av| 亚洲五月婷婷丁香| 巨乳人妻的诱惑在线观看| 日日夜夜操网爽| 欧美黑人巨大hd| 人人妻,人人澡人人爽秒播| 日韩欧美国产在线观看| 制服丝袜大香蕉在线| 动漫黄色视频在线观看| 欧美成人免费av一区二区三区| 19禁男女啪啪无遮挡网站| 亚洲精品一卡2卡三卡4卡5卡| 久久久久性生活片| 亚洲精品色激情综合| 757午夜福利合集在线观看| 国产成人aa在线观看| 国产蜜桃级精品一区二区三区| 在线观看舔阴道视频| 亚洲五月婷婷丁香| 欧美在线黄色| 中文字幕av在线有码专区| a在线观看视频网站| 成人三级黄色视频| 黄色毛片三级朝国网站| 久久久久国内视频| 国产精品久久久av美女十八| 国产精品野战在线观看| 亚洲精品久久国产高清桃花| 亚洲国产精品sss在线观看| 国产激情欧美一区二区| 母亲3免费完整高清在线观看| 男人舔奶头视频| 黑人欧美特级aaaaaa片| 欧美zozozo另类| 神马国产精品三级电影在线观看 | 久久中文字幕一级| 每晚都被弄得嗷嗷叫到高潮| 国内久久婷婷六月综合欲色啪| 欧美日韩亚洲综合一区二区三区_| 人妻夜夜爽99麻豆av| 熟女少妇亚洲综合色aaa.| 嫩草影院精品99| 色哟哟哟哟哟哟| 法律面前人人平等表现在哪些方面| 91成年电影在线观看| 久久久久久免费高清国产稀缺| 亚洲精品国产一区二区精华液| 成在线人永久免费视频| 老司机靠b影院| 伦理电影免费视频| 51午夜福利影视在线观看| 欧美+亚洲+日韩+国产| 两个人看的免费小视频| 亚洲国产精品成人综合色| 亚洲真实伦在线观看| 脱女人内裤的视频| 丝袜人妻中文字幕| 午夜福利欧美成人| 国产熟女xx|