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

    礦井紅外熱成像遠距離測溫誤差分析與精確測溫方法

    2022-05-27 09:48:34孫繼平范偉強
    煤炭學(xué)報 2022年4期
    關(guān)鍵詞:測量

    孫繼平,范偉強

    (中國礦業(yè)大學(xué)(北京) 機電與信息工程學(xué)院,北京 100083)

    熱成像儀具有非接觸、不干擾被測溫度場、監(jiān)視范圍大、響應(yīng)速度快、探測距離遠等優(yōu)點,已經(jīng)被廣泛應(yīng)用于軍事、航天、農(nóng)業(yè)、環(huán)保、醫(yī)學(xué)、工業(yè)、安全等領(lǐng)域。但現(xiàn)有熱成像儀的直接測溫結(jié)果僅適用于不超過15 m的近距離測溫,當測溫距離增大時,測溫精度將無法滿足工程實際應(yīng)用,需要進行溫度補償(以下簡稱:遠距離測溫)。在遠距離測溫過程中,熱成像儀與目標之間的測量光路受水蒸氣、CO、CH、粉塵、氣溶膠等影響,測量距離越遠,測量誤差越大,難以滿足現(xiàn)場需求。

    礦井火災(zāi)是煤礦重特大事故之一,影響著煤礦安全生產(chǎn)。礦井外因火災(zāi)一般發(fā)生在有電纜、膠帶和機電設(shè)備的巷道、硐室和采掘工作面。礦井外因火災(zāi)監(jiān)測方法主要有溫度、煙霧、氣體等傳感器監(jiān)測,分布式光纖測溫,可見光和紅外圖像等視頻圖像監(jiān)視等。溫度、煙霧、氣體等傳感器監(jiān)測礦井火災(zāi),若礦井全覆蓋布置傳感器,存在著設(shè)備數(shù)量多、成本高、維護工作量大等缺點;若傳感器稀疏布置,遠離監(jiān)測點的火源不能被及時發(fā)現(xiàn)。礦用分布式光纖測溫適用于礦用電纜和膠帶火災(zāi)的監(jiān)測,但存在著光纖易損壞、安裝復(fù)雜、維護困難等缺點?;谝曨l圖像的火災(zāi)監(jiān)測,具有可視化、監(jiān)測范圍廣等優(yōu)點,但可見光和近紅外圖像在井下受礦燈、車燈、巷道燈等光源干擾大。紅外熱成像具有可視化、監(jiān)測范圍廣、對高溫物體敏感、可全面反映監(jiān)視范圍內(nèi)的溫度場信息等優(yōu)點,但遠距離測溫誤差大,難以滿足礦井火災(zāi)監(jiān)測需求。因此,研究熱成像儀的遠距離精確測溫方法,對提高遠紅外攝像機監(jiān)視范圍和火災(zāi)報警準確性,具有重要理論意義和實用價值。

    針對紅外遠距離測溫誤差隨測溫距離變化大的問題,諸多學(xué)者對影響測溫誤差的因素進行了分析,并提出了不同的測溫補償方法。如吳晗平在綜合分析高度修正、傾斜路程以及大氣衰減等因素的基礎(chǔ)上,構(gòu)建了基于工程應(yīng)用的大氣透過率經(jīng)驗公式;部分學(xué)者基于LOWTRAN、MODTRAN等專業(yè)數(shù)據(jù)庫中不同氣體對紅外輻射的吸收特性,對不同路徑下的紅外輻射大氣衰減模型進行了研究;劉丹丹、宋福印等結(jié)合我國實際地理特征,指出不同的大氣環(huán)境有不同的紅外輻射和透射特性,直接采用國外專業(yè)數(shù)據(jù)庫中的大氣透過率模型將導(dǎo)致較大的測溫誤差;陳中偉、杜永成等通過數(shù)值模擬和實驗驗證方法,深入分析了水霧遮蔽對紅外輻射衰減的影響,并得出了不同水霧對紅外測溫距離的衰減規(guī)律;楊楨、張士成等提出的反射溫度補償法、雙溫雙波法、變譜法等測溫模型,有效提高了待測目標的溫度測量精度;孫繼平等分析了礦井巷道中大氣輻射的衰減因素和主要的氣體分子吸收特性;鐘文婷、王紅霞等對影響測溫精度的大氣氣溶膠消光系數(shù)進行了深入研究和測量。但是鮮有關(guān)于礦井紅外測溫距離與測溫誤差的分析和通過紅外熱成像灰度反演的精確測溫方法的研究報道。

    為提高礦井紅外熱成像遠距離測溫精度,實現(xiàn)對煤礦井下高溫?zé)嵩幢O(jiān)測和火災(zāi)預(yù)警,筆者研究了紅外熱成像輻射測溫原理,建立了熱成像儀接收的單波輻射照度與圖像灰度的函數(shù)關(guān)系;分析了礦井紅外輻射路徑中水蒸氣、CO、CH和SO的吸收衰減,粉塵散射衰減對紅外測溫精度的影響,并根據(jù)礦井測量環(huán)境條件,對礦井空氣總透過率進行實時修正;結(jié)合熱成像儀的輻射測溫公式,接收單波輻射照度與圖像灰度之間的函數(shù)關(guān)系,以及礦井氣體和粉塵的衰減模型,提出了紅外熱成像遠距離精確測溫方法;通過選取不同材質(zhì)的灰體目標,在不同表面溫度、不同測溫距離下,采用精確測溫方法對其進行多次測量實驗,驗證了本文所提方法的可行性和有效性。

    1 紅外熱成像輻射測溫原理

    灰體目標表面輻射出射度經(jīng)測量光路傳輸?shù)綗岢上駜x端的單波輻射亮度為

    (1)

    式中,為目標物的單波輻射亮度,W/(m·sr);()為熱成像儀與目標物之間的測溫距離為時,對應(yīng)的大氣光透過率;()為目標物表面溫度為時,在探測器接收光譜區(qū)間內(nèi)的平均法向發(fā)射率;b()為目標物表面溫度為時,對應(yīng)的黑體輻射亮度,W/(m·sr);()為目標物表面溫度為時,在探測器接收光譜區(qū)間內(nèi)的平均法向反射率;u()為環(huán)境溫度為時,對應(yīng)的背景輻射亮度,W/(m·sr);a()為目標與熱成像儀之間的大氣光譜發(fā)射率;a()為大氣溫度為時,對應(yīng)的大氣輻射亮度,W/(m·sr)。

    輻射目標物為灰體時,其發(fā)射率、吸收率、透射率與波長無關(guān),只與溫度有關(guān)。則根據(jù)基爾霍夫定律和式(1)可知,熱成像儀接收到的灰體目標單波輻射照度為

    =π[()()b()+
    ()(1-)u()+a()]

    (2)

    式中,為熱成像儀接收的單波輻射照度,W/m;(),分別為灰體表面溫度為時的表面發(fā)射率和表面吸收率;為熱成像儀最小空間張角所對應(yīng)目標的可視面積,在一定條件下,為常數(shù)。

    所研究灰體均為朗伯體,則根據(jù)朗伯體的輻射特性,被測目標物的單波輻亮度與目標物表面單波輻出度的關(guān)系為

    ()=()π

    (3)

    式中,()為目標物表面的單波輻出度,W/m。

    由斯特藩-玻爾茲曼定律可知,黑體表面在單位面積上所發(fā)射的各種波長的總輻射功率與其熱力學(xué)溫度的4次方成正比:

    (4)

    式中,為黑體表面在單位面積上的總輻射功率,W/m;(,)為黑體表面在單位面積上的單波輻射功率,W/m;為斯蒂芬-玻爾茲曼常數(shù),=5.67×10,W/(m·K)。

    在同一溫度和波長范圍內(nèi),目標物體表面輻射的功率總是小于黑體表面輻射的功率,其灰體表面在單位面積上所發(fā)射的各種波長的總輻射功率與其熱力學(xué)溫度的次方成正比,即灰體單位面積上的輻射功率()為

    ()=()

    (5)

    其中,()為灰體表面溫度為時的發(fā)射率;與探測器的接收光譜范圍有關(guān),當使用不同波段的探測器時,不同,對于8~14 μm的探測器,=4.09;對6~9 μm探測器,=5.33;對3~5 μm探測器,=9.25。

    熱成像儀端接收的輻射照度等效為黑體表面溫度為時的輻射照度,則將式(3),(5)代入式(2),并化簡后得到熱成像儀的輻射測溫公式:

    (6)

    2 礦井氣體及粉塵對紅外熱成像遠距離測溫的影響

    根據(jù)熱成像儀的輻射測溫公式可知,熱成像儀的測溫精度不僅與被測物體輻射溫度、環(huán)境反射溫度、大氣溫度有關(guān),而且與測溫路徑中的空氣透過率緊密相關(guān)。現(xiàn)有熱成像儀可直接用于近距離測溫,但在遠距離測溫時,未考慮空氣透過率對測溫精度的影響,導(dǎo)致熱成像儀測溫準確性和穩(wěn)定性較差。為實現(xiàn)煤礦井下測量環(huán)境中的遠距離精確測溫,分析了礦井氣體對紅外測溫精度的影響,并構(gòu)建了礦井空氣總透過率與輻射測溫距離的衰減模型。

    根據(jù)《煤礦安全規(guī)程》可知,煤礦井下空氣中的CH,CO,CO,HS,SO,NO,NH體積分數(shù)一般較低,其中(CH)≤1.0%,(CO)≤0.5%,(CO)≤0.002 4%,(HS)≤0.000 66%,(SO)≤0.000 5%,(NO)≤0.000 25%,(NH)≤0.004%。在礦井通風(fēng)系統(tǒng)故障時,CH,CO等氣體體積分數(shù)會上升;礦井發(fā)生火災(zāi)時,CO,CO等氣體體積分數(shù)會上升;礦井爆破作業(yè)和膠輪車運輸時,NO,NO,CO,SO等氣體體積分數(shù)會升高;煤與瓦斯突出或瓦斯異常涌出時,CH氣體體積分數(shù)會上升。通過查閱HITRAN光譜數(shù)據(jù)庫可知,CH氣體在8.6~14.0 μm波段基本不具有吸收能力,但是在8.0~8.6 μm波段有較強吸收能力,最強吸收峰8.03 μm處的吸收強度為2.479×10cm/(molecule·cm);CO氣體在8.0~14.0 μm波段不具有吸收能力;CO氣體在8.0~13.0 μm波段具有弱吸收能力,并且在13~14 μm波段有強吸收能力,最強吸收峰13.88 μm處的吸收強度為5.883×10cm/(molecule·cm);SO氣體在8~14 μm波段具有一定吸收能力,在最強吸收峰8.59 μm的吸收強度為6.175×10cm/(molecule·cm);NO氣體在8~14 μm波段具有微弱吸收能力,在最強吸收峰8.09 μm處的吸收強度僅為1.138×10cm/(molecule·cm);NO氣體在8~14 μm波段不具有吸收能力。

    因此,綜合礦井空氣中吸收性分子的光譜吸收能力和體積分數(shù),HS,NO,NO,NH,CO氣體對紅外輻射衰減的影響忽略不計,但是應(yīng)考慮空氣中CH,CO,SO氣體對紅外輻射衰減的影響。此外,煤礦井下有淋水、空氣潮濕,導(dǎo)致礦井空氣中水蒸氣含量高;煤炭采運、巷道掘進、膠輪車尾氣等造成空氣中懸浮大量的煤塵、煙塵、巖塵等粉塵微粒。

    綜上分析可知,在礦井空氣中,造成紅外輻射衰減的主要因素包括:① 空氣中水蒸氣、CO、CH、SO氣體分子吸收;② 煤塵、巖塵和煙塵等散射。假設(shè)空氣對紅外輻射的衰減能力用消光系數(shù)(衰減系數(shù))表示。則由Bougner-Lambert定律得到單波消光系數(shù)()與透過率()的關(guān)系為

    ()=e-()

    (7)

    式中,為熱成像儀與目標物之間的測溫距離, km。

    在礦井巷道中,由各種物質(zhì)引起紅外輻射衰減的透過率相乘得到空氣中的單波透過率:

    ()=()(HO)()(CO)×
    ()(SO)()(CH)()

    (8)

    式中,()(HO)為水蒸氣吸收后的單波透過率;()(CO)為CO吸收后的單波透過率;()(CH)為CH吸收后的單波透過率;()(SO)為SO吸收后的單波透過率;()為空氣散射后的單波透過率。

    2.1 水蒸氣的吸收衰減

    由于水蒸氣對紅外吸收與測溫光路中水分子總數(shù)有關(guān)。在一段大氣光路內(nèi),水蒸氣的含量用凝結(jié)成液態(tài)水柱的高度(可降水分毫米數(shù))表示,并可以通過空氣中溫濕度等參數(shù)計算得到。當大氣壓強和環(huán)境溫度保持不變,在海平面上的輻射距離為1 km時,大氣中水蒸氣凝結(jié)成液態(tài)水柱的高度為

    (9)

    其中,為相對濕度時,空氣中水蒸氣凝結(jié)成液態(tài)水柱的高度;為相對濕度(100%)時,空氣中的水蒸氣凝結(jié)成液態(tài)水柱的高度;為相對濕度時,空氣中的水蒸氣凝結(jié)成液態(tài)水柱的高度。當海平面上的輻射距離為1 km時,空氣中的水蒸氣凝結(jié)成液態(tài)水柱的高度, mm/km。

    在標準大氣壓下,海平面上的水蒸氣含量與平均大氣透過率之間的關(guān)系數(shù)據(jù),見表1。

    通過數(shù)據(jù)擬合可知,空氣中水蒸氣含量與透過率之間的關(guān)系符合指數(shù)變化規(guī)律。其函數(shù)為

    (HO)=0999 3e-0015 04

    (10)

    式中,(HO)為熱成像儀與目標物之間的水蒸氣透過率。

    2.2 CO2的吸收衰減

    在標準大氣壓下,海平面水平路程上的CO體積分數(shù)與平均大氣透過率之間的關(guān)系數(shù)據(jù),見表2。

    表1 海平面上的水蒸氣含量與平均大氣透過率的關(guān)系

    表2 海平面上CO2體積分數(shù)與平均大氣透過率之間的關(guān)系

    通過數(shù)據(jù)擬合得到海平面處CO體積分數(shù)與平均大氣透過率的函數(shù)關(guān)系為

    (CO)=0199e-0526 8(CO)+0800 6e-0028(CO)

    (11)

    式中,(CO)為CO的透過率;(CO)為相對于大氣的CO體積分數(shù)變化。

    2.3 CH4,SO2的吸收衰減

    在實際測量環(huán)境中,CH分子的吸收系數(shù)隨溫度和壓強會發(fā)生變化。因此,將吸收系數(shù)表示為與溫度和壓強有關(guān)的函數(shù):

    ()=()(,)(CH)

    (12)

    式中,()為每單位體積內(nèi)單個分子的譜線強度,()=(),cm/(molecule·cm·atm);(,)為譜線線型函數(shù);為海拔高度處的大氣壓強, atm,1 atm=101 kPa;(CH)為CH氣體體積分數(shù);為氣體分子的吸收譜線強度,cm/(molecule·cm),通過HITRAN數(shù)據(jù)庫可查詢;為氣體分子數(shù)密度;,分別為波數(shù)和中心波數(shù),cm。

    CH氣體吸收譜線線型函數(shù)包括:多普勒線型函數(shù)、洛倫茲線型函數(shù)和Voigt線型函數(shù)。其中,多普勒線型函數(shù)用于低溫度、低氣壓的環(huán)境條件;洛倫茲線型函數(shù)用于分子之間有碰撞展寬的環(huán)境條件。SO和CH具有相同的吸收線型函數(shù),且吸收系數(shù)與溫度和壓強之間的函數(shù)關(guān)系一致。因此,在SO和CH氣體體積分數(shù)測量過程中,根據(jù)不同的環(huán)境溫度和壓強選擇不同的吸收譜線線型函數(shù)。則由式(7)和式(12)得到SO和CH氣體吸收衰減后的積分透過率為

    (13)

    式中,(),()分別為CH和SO氣體的單波消光系數(shù);為積分下限;為積分上限。

    2.4 粉塵微粒的散射衰減

    在工程計算時,散射系數(shù)可利用氣象視程來表示。氣象視程是指目標與背景的對比度隨著距離的增加而減少到2%的距離。氣象能見度表示目標與背景的對比度為1,通過大氣衰減后可感知的對比度為2%時的距離(通常取0.60 μm或者0.55 μm處的目標和背景對比度)。氣象能見度、散射系數(shù)()與目標對比度的函數(shù)關(guān)系為

    ()=(0)exp[-()]

    (14)

    其中,為氣象能見度,km;(0)=1;()=0.02;(λ)為大氣散射系數(shù),km。根據(jù)式(14)得到散射系數(shù)()與氣象視程()的關(guān)系為

    (15)

    式中,通常選擇波長0.55 μm(此波長的大氣實際吸收為0);為視覺對比度閾值。

    對于遠紅外熱輻射,煤礦巷道中粉塵粒徑不大于100 μm,且粉塵粒徑符合Rosin-Rammler分布函數(shù),即粒徑大于0.8 μm的粉塵質(zhì)量占全塵的95%以上。因此,在工程應(yīng)用中,通常忽略瑞利散射的影響。則大氣的散射系數(shù)可表示為

    ()=-

    (16)

    其中,為待定常數(shù);為經(jīng)驗常數(shù)。如果已知0.55 μm處的氣象視程,則根據(jù)式(14),(15)可得該波長下的散射系數(shù),并確定待定常數(shù):

    (17)

    根據(jù)式(16),(17)得到長波紅外譜線的散射系數(shù)為

    (18)

    式中,=0.55 μm;為經(jīng)驗常數(shù),取值為

    由式(18)可得:在8~14 μm波段大氣散射衰減后的積分透過率為

    (19)

    式中,為單波散射消光系數(shù);為大氣散射衰減后的積分透過率。

    2.5 礦井空氣總透過率及透過率參數(shù)修正

    2.5.1 礦井空氣總透過率

    根據(jù)式(10),(11),(13),(19)得到的水蒸氣、CO、CH、SO的吸收透過率和粉塵微粒的散射透過率,代入式(2)可得到8~14 μm遠紅外波段大氣總的積分透過率為

    (20)

    2.5.2 礦井空氣透過率參數(shù)修正

    礦井空氣中水蒸氣、CO等氣體對灰體目標的輻射吸收本領(lǐng),隨環(huán)境溫度和大氣壓強的變化而變化。因此,在不同海拔地區(qū)進行測量過程中,需要對水蒸氣和CO的透過率進行修正。筆者采用距離系數(shù)修正熱成像儀測量現(xiàn)場海拔高度處的水平測溫距離,估計等效海平面水平測溫距離。則可用式(21)近似計算:

    (21)

    式中,為等效海平面水平測溫距離;為海拔高度處的水平測溫距離,km;為海平面處的大氣壓強;為常數(shù),水蒸汽的=0.5,CO的=1.5。

    由式(21)可知,環(huán)境溫度變化對氣體輻射吸收本領(lǐng)的影響很小,可以忽略不計。在實際工程應(yīng)用中,通常采用指數(shù)近似函數(shù)表示

    =e-

    (22)

    式中,為常數(shù),水蒸汽的=0.065 4,CO的=0.19。

    根據(jù)阿伏加德羅定律可知,不同高度的氣體分子密度與溫度和壓強有關(guān),當溫度恒定時,分子密度隨壓強的增加而增大。通常,在海拔2 km范圍內(nèi),高度每升高100 m,大氣壓力下降約1.11 kPa。因此,根據(jù)熱成像儀測溫位置的海拔高度,修正氣體密度變化引起的測量誤差。在工程應(yīng)用中,由高度變化引起的相對濕度測量誤差采用式(23)進行修正:

    =′e-

    (23)

    式中,為海平面處的相對濕度;′為海拔高度處的相對濕度;為恒常數(shù),=0.45,km。

    根據(jù)已有的研究結(jié)論可知,CO在海拔高度處的水平吸收路徑折算成海平面處的水平吸收路徑為

    =′e-

    (24)

    式中,為海平面處的水平吸收路徑;′為海拔高度處的水平吸收路徑,′=;為恒常數(shù),=0.123 km。

    由《煤炭工業(yè)礦井設(shè)計規(guī)范》GB 50215—2015的規(guī)定,礦井通風(fēng)負(正)壓不應(yīng)超過2.94 kPa。礦井通風(fēng)壓力遠小于標準大氣壓101 kPa。因此,在工程應(yīng)用中,由井下巷道通風(fēng)壓力引起的壓強變化可忽略不計。當?shù)V井巷道中水蒸氣含量,以及井下巷道測溫處的實際高度已知,則結(jié)合式(22),(23)得到海平面處等效測溫距離與海平面處水蒸氣相對濕度的乘積為

    =e-0065 4e-045=e-0515 4

    (25)

    同理,當CO體積分數(shù)和井下巷道測溫處的實際高度已知時,則結(jié)合式(22),(23)得到海平面處等效測溫距離與海平面處的水平吸收路徑的乘積為

    ′=e-019e-0123=e-0313

    (26)

    通過式(20),(25),(26)構(gòu)建礦井空氣總透過率與輻射測溫距離的衰減模型:

    (27)

    3 礦井紅外熱成像遠距離精確測溫方法

    由熱成像儀探測器的成像原理可知,探測器中某一探測單元所對應(yīng)的圖像灰度值與熱成像儀光學(xué)鏡頭面積、光學(xué)鏡頭透過率、探測器的光譜響應(yīng)系數(shù)和熱成像儀的灰度標定參數(shù)有關(guān)。此外,光學(xué)鏡頭透過率與波長有關(guān),且不同波長的入射輻射在探測器上的光譜響應(yīng)系數(shù)不同。為此,筆者構(gòu)建了熱成像儀探測器中某一探測單元所對應(yīng)的圖像灰度均值與熱成像儀接收到的單波輻射照度之間的函數(shù)關(guān)系:

    (28)

    式中,為熱成像儀的灰度標定參數(shù);為熱成像儀光學(xué)鏡頭面積,mm;為探測器的光譜響應(yīng)系數(shù)與光學(xué)鏡頭透過率的乘積。

    (29)

    由式(29)得到被測目標物表面的真實溫度為

    (30)

    根據(jù)文獻[27,36]可知,在礦井巷道半徑大于1 m的直巷道中,8~14 μm波段的紅外輻射引起的衰減損耗極小,近似為自由傳播,其輻射衰減主要為礦井水蒸氣、吸收性氣體、粉塵等自損耗媒質(zhì)。因此,根據(jù)構(gòu)建的礦井空氣總透過率與輻射測溫距離的衰減模型,熱成像儀接收到的單波輻射照度與圖像灰度之間的函數(shù)關(guān)系,以及簡化后的熱成像儀輻射測溫公式,提出了礦井紅外熱成像遠距離精確測溫方法。同時,聯(lián)立式(27),(28)和(30)得到精確測溫方法測量后的被測目標物表面的真實溫度。

    (31)

    4 試驗驗證與分析

    為驗證上述模型的正確性,筆者選取噴漆金屬板、硅橡膠板2種材料用于分析易引起礦井火災(zāi)的機電設(shè)備、硅橡膠電纜和運輸皮帶表面溫度變化及測溫距離與測溫精度之間的關(guān)系。

    (1)試驗環(huán)境。本試驗在模擬巷道環(huán)境中進行驗證,試驗環(huán)境中無高于被測物體溫度的干擾源和高體積分數(shù)的CH,SO氣體,并查閱相關(guān)資料可知標準大氣條件下的海平面大氣壓為101.325 kPa,北京市平均海拔為43.71 m,北京平均大氣壓為100.870 kPa。因此,本試驗過程中可忽略大氣壓強、環(huán)境中高溫干擾源和CH,SO氣體吸收衰減引起的溫度測量誤差。

    (2)試驗儀器。熱成像儀采用Fotric 826防火型在線式紅外熱成像儀,其空間分辨率0.32~1.27,探測器響應(yīng)波段8~14 μm,搭配15 mm標準鏡頭;待測溫的灰體目標為噴漆金屬板和硅橡膠板;加熱裝置為硅橡膠加熱面板(內(nèi)部左上、左下、右上、右下嵌入熱電偶);采用伯虎物聯(lián)智能空氣檢測儀測量環(huán)境中溫濕度、粉塵和CO體積分數(shù),具體試驗裝置系統(tǒng)如圖1所示。通常情況下,熱成像儀探測器所對應(yīng)波段的光譜響應(yīng)系數(shù)和光學(xué)鏡頭的透過率均較為平坦,可假設(shè)為常數(shù),本文中=0.9。

    圖1 試驗裝置系統(tǒng)Fig.1 Test device system

    (3)實驗過程。將厚度2 mm的灰體目標粘貼于加熱面板表面,使表面完全黏連,同時保持熱成像儀與待測灰體目標的法線在同一水平線上,且使成像清晰。測溫前對加熱面板預(yù)設(shè)不同溫度,使目標物進行階梯式加熱,為50,70,80,90,100,110,120和130 ℃。

    在標定位置(1 m處),通過熱成像儀測量不同預(yù)設(shè)溫度下噴漆金屬板、硅橡膠板的表面真實溫度;每次保持加熱面板的預(yù)設(shè)溫度穩(wěn)定后讀取其表面真實溫度,并進行20次測量取平均值,測得目標物表面的真實溫度。通過改變加熱面板的預(yù)設(shè)溫度或測溫距離,獲取待測目標物在不同預(yù)設(shè)溫度或測溫距離時,熱成像儀上的熱圖像所對應(yīng)的灰度和熱成像儀自帶算法求得的輻射溫度,每一組數(shù)據(jù)測量20次后,分別記錄輻射溫度、測溫距離和圖像灰度的平均測量值。

    4.1 目標物表面溫度變化對溫度測量的影響

    4.1.1 表面噴漆金屬板實驗

    在環(huán)境溫度21.89 ℃、相對濕度53%、CO體積分數(shù)640×10、粉塵質(zhì)量濃度118 μg/m、空氣透過率0.77的試驗條件下。首先,在熱成像儀標定位置(1 m處)測量不同預(yù)設(shè)溫度下噴漆金屬板表面的真實溫度;接著,通過熱成像儀對測溫距離為35 m處的噴漆金屬板進行測溫,并獲取不同相對應(yīng)的輻射溫度和圖像灰度;最后,記錄實測的,,數(shù)據(jù)和通過精確測溫方法求得的溫度。同時,計算與的絕對誤差(Absolute Error,AE)、相對誤差(Relative Error,RE),以及與的和,試驗數(shù)據(jù)見表3。

    表3 表面噴漆金屬板實驗數(shù)據(jù)

    由表3中的,和,與的和,以及與的和,分別繪制不同預(yù)設(shè)溫度下的真實溫度、輻射溫度和精確測溫方法的溫度的變化曲線(圖2),輻射溫度與真實溫度的絕對誤差直方圖(圖3)和相對誤差折線(圖4),以及精確測溫方法的溫度與真實溫度的絕對誤差直方圖(圖3)和相對誤差折線(圖4)。

    圖2 不同預(yù)設(shè)溫度下真實溫度、輻射溫度、精確測溫方法的測量溫度Fig.2 Real temperature,radiation temperature,accuratemeasurement temperature under different preset temperatures

    圖3 絕對誤差隨測量溫度變化的直方圖Fig.3 Histogram of absolute error with measured temperature

    圖4 相對誤差隨測量溫度變化的折線Fig.4 Line graph of relative error with measured temperature

    根據(jù)表3、圖2~4可知:隨著噴漆金屬板的真實溫度逐漸升高,其與的和與的均趨向變大,與的和與的整體趨于穩(wěn)定。在54.43~135.63 ℃,與的最大達到13.62 ℃,最大達到11.35%,最小為4.79 ℃,最小也達到8.80%,平均達到10.19%;而與的整體較小,其最大為4.87 ℃,最大為3.59%,最小達到1.10 ℃,最小僅為1.80%,平均小于2.31%。采用精確測溫方法求得的平均比紅外熱像儀自帶算法測得的平均減小了7.88%。

    由上可知,在遠距離測溫時,高溫?zé)嵩茨繕藴囟茸兓瘜E的影響較小,但測溫路徑中的紅外輻射衰減對測溫精度有較大影響。熱成像儀自帶算法測得的輻射溫度和精確測溫方法的測量溫度均小于待測目標物的真實溫度,但更接近于真實溫度。進一步驗證了采用本文提出的精確測溫方法更能滿足實際測量要求。

    4.1.2 硅橡膠板實驗

    在環(huán)境溫度21.32 ℃、相對濕度45%、CO體積分數(shù)630×10、粉塵質(zhì)量濃度99 μg/m、空氣透過率0.81的試驗條件下。首先,使用熱成像儀標定位置(1 m處)測量不同預(yù)設(shè)溫度下硅橡膠板表面的真實溫度;接著,通過熱成像儀測量測溫距離為35 m處的硅橡膠板,并獲取不同預(yù)設(shè)溫度相對應(yīng)的輻射溫度和圖像灰度;最后,記錄實測的,,數(shù)據(jù),以及通過精確測溫方法求得的溫度。同時,計算與的和,以及與的和,試驗數(shù)據(jù)見表4。

    表4 硅橡膠板實驗數(shù)據(jù)

    由表4中的,和,與的和,以及與的和,分別繪制不同預(yù)設(shè)溫度下真實溫度、輻射溫度和精確測溫方法的測量溫度的變化曲線(圖5),輻射溫度與真實溫度的絕對誤差直方(圖6)和相對誤差折線(圖7),以及精確測溫方法的測量溫度與真實溫度的絕對誤差直方(圖6)和相對誤差折線(圖7)。

    圖5 不同預(yù)設(shè)溫度下的真實溫度、輻射溫度、精確測溫方法的測量溫度Fig.5 Real temperature,radiation temperature,accuratemeasurement temperature under different preset temperatures

    圖6 絕對誤差隨測量溫度變化的直方圖Fig.6 Histogram of absolute error with measured temperature

    圖7 相對誤差隨測量溫度變化的折線Fig.7 Line graph of relative error with measured temperature

    根據(jù)表4、圖5~7可知:隨著硅橡膠板真實溫度的升高,與的和與的均呈現(xiàn)逐漸增長的趨勢,與的和與的整體趨于穩(wěn)定。在56.68~139.42 ℃,與的,均小于與的,,與的最大達到18.31 ℃,最大達到15.11%,最小達到7.12 ℃,最小也達到12.55%,平均高于14.24%;與的最大為5.42 ℃,最大為4.55%,最小僅為1.92 ℃,最小達到2.62%,平均僅為3.39%。采用精確測溫方法求得的平均比紅外熱像儀自帶算法測得的平均減小了10.85%。

    通過上述分析,進一步驗證了在遠距離測溫時,高溫?zé)嵩茨繕藴囟茸兓瘜Φ挠绊戄^小,但測溫路徑中的紅外輻射衰減對測溫精度有較大影響。在遠距離測溫時,測溫路徑中不同高溫?zé)嵩茨繕说募t外輻射衰減對測溫精度的影響均不容忽視。

    4.2 目標物距離變化對溫度測量的影響

    利用熱成像儀分別采集測溫距離為15,20,25,30,35 m時的噴漆金屬板和硅橡膠板的輻射溫度和對應(yīng)的圖像灰度,同時記錄不同測溫距離下的輻射溫度和圖像平均灰度。

    4.2.1 表面噴漆金屬板實驗

    在環(huán)境溫度21.89 ℃、相對濕度53%、CO體積分數(shù)640×10、粉塵質(zhì)量濃度118 μg/m、表面發(fā)射率0.97,表面真實溫度為106.25 ℃條件下,分別記錄不同測溫距離下的輻射溫度和精確測溫方法的測量溫度,并分別計算與的和,以及與的和,實驗結(jié)果見表5。

    表5 表面噴漆金屬板距離實驗數(shù)據(jù)

    由表5中的與的和,以及與的和,分別繪制不同測溫距離下與的絕對誤差直方圖和相對誤差折線,以及不同測溫距離下與的絕對誤差直方圖和相對誤差折線。絕對誤差直方圖如圖8所示,相對誤差折線如圖9所示。

    圖8 絕對誤差隨測溫距離變化的直方圖Fig.8 Histogram of absolute error with temperaturemeasurement distance

    圖9 相對誤差隨測溫距離變化的折線Fig.9 Line graph of relative error with temperaturemeasurement distance

    根據(jù)表5、圖8~9可知,隨著測溫距離的逐漸變遠,紅外輻射在空氣中的透過率逐漸減小,其與的和與的,以及與的和與的均呈現(xiàn)逐漸增大的趨勢,但與的,值均優(yōu)于與的,,試驗所測結(jié)果符合理論推導(dǎo),空氣透過率變化對紅外熱像儀自帶算法的精度和精確測溫方法的精度均產(chǎn)生較大影響。

    在測溫距離15~35 m內(nèi),熱成像儀自帶算法對應(yīng)的從-4.28%下降至-11.33%,即測溫距離每增遠1 m,平均增加0.35%;精確測溫方法對應(yīng)的從0.25%下降至-1.90%,即測溫距離每增遠1 m,平均增加0.08%。進一步表明,采用礦井空氣透過率修正的遠距離精確測溫方法能夠適用于紅外遠距離精確測溫。

    4.2.2 硅橡膠板實驗

    在環(huán)境溫度21.32 ℃、相對濕度45%、CO體積分數(shù)630×10、粉塵質(zhì)量濃度99 μg/m、表面發(fā)射率0.95、表面真實溫度為107.91 ℃的試驗條件下,分別記錄不同測溫距離下的輻射溫度與精確測溫方法的測量溫度,并分別計算與的和,以及與的和,實驗結(jié)果見表6。

    由表6中的與的和,以及與的和,分別繪制不同測溫距離下與的絕對誤差直方圖和相對誤差折線,以及不同測溫距離下與的絕對誤差直方圖和相對誤差折線。絕對誤差直方圖如圖10所示,相對誤差折線如圖11所示。

    表6 硅橡膠板距離實驗數(shù)據(jù)

    圖10 絕對誤差隨測溫距離變化的直方圖Fig.10 Histogram of absolute error with temperaturemeasurement distance

    圖11 相對誤差隨測溫距離變化的折線Fig.11 Line graph of relative error with temperaturemeasurement distance

    由表6、圖10,11可知,隨著測溫距離的逐漸變遠,紅外輻射在空氣中的透過率逐漸減小,其與的和與的,以及與的和與的均呈現(xiàn)逐漸增大,但與的,均大于與的,,試驗所測結(jié)果符合理論推導(dǎo)。

    在測溫距離15~35 m內(nèi),熱成像儀自帶算法對應(yīng)的從-6.88%下降至-15.08%,即測溫距離每增遠1 m,平均增加0.41%;精確測溫方法對應(yīng)的為-2.30%下降至-3.29%,即測溫距離每增遠1 m,平均增加0.05%。進一步表明,引入礦井空氣透過率修正的遠距離精確測溫方法求得的測溫精度更高,更能滿足實際測溫要求。

    4.3 煤炭燃燒遠距離測溫結(jié)果比較

    根據(jù)4.1和4.2節(jié)的試驗結(jié)論可知,本文提出的遠距離精確測溫方法可用于監(jiān)測礦井外因火災(zāi)。為進一步驗證該方法的有效性,在環(huán)境溫度23.6 ℃、相對濕度91%、表面真實溫度721.4 ℃(采用開普森進口鎧裝熱電偶)、CO體積分數(shù)415×10、粉塵質(zhì)量濃度55 μg/m條件下,分別采集不同測溫距離下的煤炭燃燒狀態(tài)輻射溫度和圖像灰度,計算精確測溫方法的測量溫度。不同測溫距離下的實際測量結(jié)果見表7。

    表7 不同測溫距離下煤炭燃燒的試驗數(shù)據(jù)

    由表7中的與的和,以及與的和,分別繪制不同測溫距離下與的絕對誤差直方圖和相對誤差折線,以及不同測溫距離下與的絕對誤差直方圖和相對誤差折線。絕對誤差直方圖如圖12所示,相對誤差折線如圖13所示。

    圖12 絕對誤差隨測溫距離變化的直方圖Fig.12 Histogram of absolute error with temperaturemeasurement distance

    圖13 相對誤差隨測溫距離變化的折線Fig.13 Line graph of relative error with temperaturemeasurement distance

    根據(jù)表7、圖12,13可知,隨著測溫距離的逐漸變遠,紅外輻射在空氣中的透過率逐漸減小,輻射溫度與真實溫度的,逐漸增大,精確測溫方法的測量溫度與真實溫度的,也逐漸增大,但與的,顯著優(yōu)于與的輻射測溫,。

    在測溫距離35~140 m內(nèi),熱成像儀自帶算法對應(yīng)的從-28.22%下降至-62.52%,即測溫距離每增遠1 m,平均增加0.33%;精確測溫方法對應(yīng)的為-0.08%下降至-14.49%,即測溫距離每增遠1 m,平均僅增加0.14%。進一步表明本文提出的精確測溫方法可用于煤礦井下的遠距離火災(zāi)監(jiān)測,并提高火災(zāi)預(yù)警的準確率。

    5 結(jié) 論

    (1)分析了礦井空氣中主要成分的濃度和紅外吸收強度,得到了礦井紅外輻射衰減的主要影響因素,包括:① 空氣中水蒸氣、CO、CH和SO氣體分子的吸收;② 煤塵、巖塵和煙塵等粉塵顆粒的散射。建立了礦井紅外輻射路徑中水蒸氣、CO、CH、SO氣體的吸收衰減和粉塵顆粒的散射衰減函數(shù),實現(xiàn)了礦井空氣中總透過率的估計。

    (2)分析了礦井環(huán)境中溫度、濕度、氣壓、測溫距離、氣體的吸收衰減和粉塵顆粒的散射衰減等參數(shù)對測溫精度的影響,構(gòu)建了礦井空氣總透過率與輻射測溫距離的衰減模型;分析了目標物體在熱成像儀中的成像灰度,灰度標定系數(shù)、光學(xué)鏡頭面積、探測器響應(yīng)系數(shù)、光學(xué)鏡頭透過率等參數(shù),得到了熱成像儀接收的輻射照度;由礦井空氣總透過率、熱成像儀接收的輻射照度和簡化后的熱成像儀輻射測溫公式,計算出了遠距離目標物體的表面溫度。

    (3)本文提出的精確測溫方法得到的灰體目標表面溫度與目標的真實溫度誤差較小,優(yōu)于熱成像儀自帶測溫算法的測溫精度。在不同預(yù)設(shè)溫度下噴漆金屬板和硅橡膠板試驗中,精確測溫方法的平均相對誤差為2.85%,熱成像儀自帶算法的平均相對誤差為12.22%。在不同測溫距離下噴漆金屬板、硅橡膠板和煤炭燃燒試驗中,熱像儀自帶測溫算法得到的每米相對誤差分別為0.35%,0.41%,0.33%,精確測溫方法得到的每米相對誤差分別為0.08%,0.05%和0.14%,溫度測量誤差大幅下降。

    (4)由試驗結(jié)果可知,熱源目標表面溫度變化對測溫精度的影響較小,但隨著測溫距離的變遠,紅外輻射的總透過率逐漸減小,即測溫路徑中輻射衰減對測量灰體目標溫度的影響更大。此外,本文提出的紅外熱成像遠距離精確測溫方法,不僅可以用于礦井下高溫?zé)嵩春突馂?zāi)的遠距離監(jiān)測,通過修改測溫路徑中的不同氣體吸收參數(shù)后,還可以應(yīng)用于其他場所的大面積高溫?zé)嵩春突馂?zāi)監(jiān)視。

    猜你喜歡
    測量
    測量重量,測量長度……
    把握四個“三” 測量變簡單
    滑動摩擦力的測量和計算
    滑動摩擦力的測量與計算
    測量的樂趣
    二十四節(jié)氣簡易測量
    日出日落的觀察與測量
    滑動摩擦力的測量與計算
    測量
    測量水的多少……
    最新的欧美精品一区二区| 久久狼人影院| 日韩一区二区三区影片| 在线观看国产h片| 一级片免费观看大全| 女性生殖器流出的白浆| 亚洲精品久久成人aⅴ小说| 亚洲精品日本国产第一区| 久久影院123| 亚洲av.av天堂| 亚洲av中文av极速乱| 国产一区二区激情短视频 | 亚洲欧美一区二区三区国产| 亚洲精品av麻豆狂野| 国产一区二区三区av在线| 在线精品无人区一区二区三| 国产精品久久久久久久久免| 久久精品国产综合久久久| 国产精品三级大全| 日本色播在线视频| 国产精品av久久久久免费| 日韩欧美一区视频在线观看| 国产成人欧美| 欧美人与性动交α欧美软件| 久久影院123| 久久久久久伊人网av| 午夜日本视频在线| 九草在线视频观看| 巨乳人妻的诱惑在线观看| 超碰97精品在线观看| 2021少妇久久久久久久久久久| 欧美精品av麻豆av| 精品一区二区免费观看| 国产精品成人在线| 亚洲精品日韩在线中文字幕| 国产不卡av网站在线观看| 中文字幕亚洲精品专区| 我要看黄色一级片免费的| 校园人妻丝袜中文字幕| 午夜激情久久久久久久| 国产极品粉嫩免费观看在线| 精品人妻偷拍中文字幕| 久久国产精品大桥未久av| 捣出白浆h1v1| 日韩视频在线欧美| 久久久精品免费免费高清| 美女午夜性视频免费| videos熟女内射| 热re99久久精品国产66热6| 精品卡一卡二卡四卡免费| www.精华液| 久久这里只有精品19| 青草久久国产| 久久久亚洲精品成人影院| 欧美日韩一区二区视频在线观看视频在线| h视频一区二区三区| 999久久久国产精品视频| 国产精品99久久99久久久不卡 | 18禁裸乳无遮挡动漫免费视频| 国产精品无大码| 边亲边吃奶的免费视频| 久久精品人人爽人人爽视色| 久久毛片免费看一区二区三区| 国产乱人偷精品视频| 另类精品久久| 老司机影院成人| 免费观看无遮挡的男女| av一本久久久久| av片东京热男人的天堂| 日韩不卡一区二区三区视频在线| 成人毛片a级毛片在线播放| 久久久久精品人妻al黑| 成年动漫av网址| 精品亚洲成a人片在线观看| 久久久久久人人人人人| 久久精品久久久久久久性| 亚洲精品日本国产第一区| 精品久久久精品久久久| 高清不卡的av网站| 91久久精品国产一区二区三区| 日韩av免费高清视频| 啦啦啦在线观看免费高清www| 99久久中文字幕三级久久日本| 欧美日韩精品成人综合77777| 1024视频免费在线观看| 制服丝袜香蕉在线| 边亲边吃奶的免费视频| 视频在线观看一区二区三区| 波野结衣二区三区在线| 欧美激情 高清一区二区三区| 一级爰片在线观看| 日韩制服丝袜自拍偷拍| 国产亚洲一区二区精品| 美女国产高潮福利片在线看| 亚洲国产最新在线播放| 看免费av毛片| 在线观看国产h片| 国产亚洲午夜精品一区二区久久| 中国三级夫妇交换| 26uuu在线亚洲综合色| 免费高清在线观看视频在线观看| 免费观看在线日韩| 波野结衣二区三区在线| 成年人免费黄色播放视频| 1024香蕉在线观看| 美女中出高潮动态图| 精品亚洲成a人片在线观看| 三上悠亚av全集在线观看| 三级国产精品片| 9热在线视频观看99| 亚洲婷婷狠狠爱综合网| 国产av国产精品国产| 欧美国产精品一级二级三级| 亚洲情色 制服丝袜| 精品一区在线观看国产| 亚洲精品国产一区二区精华液| 成人国产av品久久久| 我要看黄色一级片免费的| 日本欧美国产在线视频| 国产精品久久久久成人av| 久久 成人 亚洲| 99香蕉大伊视频| 电影成人av| 天天操日日干夜夜撸| 午夜福利网站1000一区二区三区| 一级片免费观看大全| 天天躁夜夜躁狠狠久久av| 免费人妻精品一区二区三区视频| 又大又黄又爽视频免费| 国产精品国产三级专区第一集| 交换朋友夫妻互换小说| 男女免费视频国产| 成年美女黄网站色视频大全免费| 人人妻人人澡人人爽人人夜夜| 日韩,欧美,国产一区二区三区| 亚洲第一青青草原| 熟女av电影| 丝袜美腿诱惑在线| 天天躁夜夜躁狠狠久久av| 日本免费在线观看一区| 人体艺术视频欧美日本| 亚洲人成77777在线视频| 日韩熟女老妇一区二区性免费视频| 久久av网站| 亚洲男人天堂网一区| 黄色一级大片看看| 成年美女黄网站色视频大全免费| 国产免费现黄频在线看| 岛国毛片在线播放| 成人免费观看视频高清| 国产成人精品在线电影| videosex国产| 精品国产一区二区三区久久久樱花| 香蕉丝袜av| 香蕉丝袜av| 久久精品国产亚洲av高清一级| 日韩av在线免费看完整版不卡| 国产精品国产三级专区第一集| 成人手机av| 婷婷色综合大香蕉| 91aial.com中文字幕在线观看| 观看av在线不卡| 久久久久久伊人网av| 国产精品欧美亚洲77777| 国产精品欧美亚洲77777| 国产成人91sexporn| 亚洲精品久久久久久婷婷小说| 交换朋友夫妻互换小说| 人妻一区二区av| 国产不卡av网站在线观看| 性色av一级| 中文字幕亚洲精品专区| 在线观看人妻少妇| 国产伦理片在线播放av一区| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 丝瓜视频免费看黄片| 国产成人欧美| 91成人精品电影| 亚洲第一av免费看| 日韩视频在线欧美| 女人久久www免费人成看片| av网站免费在线观看视频| 伊人亚洲综合成人网| 男女啪啪激烈高潮av片| 亚洲精品久久久久久婷婷小说| 在线观看免费日韩欧美大片| av免费观看日本| 亚洲欧美成人综合另类久久久| 久久久久久伊人网av| 精品久久蜜臀av无| 国产精品免费大片| 有码 亚洲区| 久久影院123| 免费观看无遮挡的男女| 欧美日韩精品网址| 另类精品久久| 欧美av亚洲av综合av国产av | 国产 精品1| 秋霞在线观看毛片| 久久久精品区二区三区| 欧美 日韩 精品 国产| 国产乱人偷精品视频| 丝袜喷水一区| 久热这里只有精品99| 十八禁网站网址无遮挡| 国产高清不卡午夜福利| 交换朋友夫妻互换小说| 亚洲美女黄色视频免费看| 人人妻人人澡人人看| 久久这里只有精品19| 各种免费的搞黄视频| 黄片播放在线免费| 99久久综合免费| 热re99久久国产66热| www日本在线高清视频| 国产日韩欧美视频二区| 国产欧美日韩一区二区三区在线| 久久久a久久爽久久v久久| 国产一区二区三区综合在线观看| 在线天堂中文资源库| 夫妻性生交免费视频一级片| 男女下面插进去视频免费观看| av在线观看视频网站免费| 99热国产这里只有精品6| 国产激情久久老熟女| 国产综合精华液| 亚洲欧美清纯卡通| 亚洲欧美日韩另类电影网站| 亚洲伊人色综图| 国产午夜精品一二区理论片| 日本黄色日本黄色录像| 欧美xxⅹ黑人| 超色免费av| 久久久精品国产亚洲av高清涩受| 亚洲内射少妇av| 啦啦啦视频在线资源免费观看| 中文字幕精品免费在线观看视频| 亚洲av综合色区一区| 亚洲精品一区蜜桃| 亚洲视频免费观看视频| 国产有黄有色有爽视频| 午夜av观看不卡| 美女国产视频在线观看| 免费av中文字幕在线| 亚洲欧美一区二区三区久久| 9热在线视频观看99| 18禁观看日本| 1024香蕉在线观看| 国产福利在线免费观看视频| 国产精品偷伦视频观看了| 久久这里只有精品19| 免费久久久久久久精品成人欧美视频| 我的亚洲天堂| 又黄又粗又硬又大视频| 麻豆乱淫一区二区| 久久精品国产亚洲av天美| 成年女人在线观看亚洲视频| 久久亚洲国产成人精品v| 久久女婷五月综合色啪小说| 熟女电影av网| 久久免费观看电影| 搡女人真爽免费视频火全软件| 亚洲精品国产一区二区精华液| 免费av中文字幕在线| 日韩免费高清中文字幕av| 欧美国产精品va在线观看不卡| 日本欧美国产在线视频| 999久久久国产精品视频| 欧美+日韩+精品| 久久99精品国语久久久| av片东京热男人的天堂| 欧美在线黄色| 国产深夜福利视频在线观看| 久久狼人影院| 国产成人精品在线电影| 欧美精品一区二区免费开放| 一边亲一边摸免费视频| 嫩草影院入口| 国产无遮挡羞羞视频在线观看| 亚洲精华国产精华液的使用体验| 男女下面插进去视频免费观看| 两性夫妻黄色片| 日本猛色少妇xxxxx猛交久久| 少妇猛男粗大的猛烈进出视频| 热re99久久精品国产66热6| 国产成人精品一,二区| 久久鲁丝午夜福利片| 免费高清在线观看视频在线观看| 少妇的丰满在线观看| 在线观看免费高清a一片| 免费看av在线观看网站| 精品亚洲乱码少妇综合久久| 亚洲男人天堂网一区| 国产精品 国内视频| 亚洲天堂av无毛| 又粗又硬又长又爽又黄的视频| av在线播放精品| 黄色视频在线播放观看不卡| av在线app专区| 亚洲精品美女久久av网站| 久久久亚洲精品成人影院| 色吧在线观看| 日日摸夜夜添夜夜爱| 建设人人有责人人尽责人人享有的| 丁香六月天网| 国产成人精品福利久久| 国产极品天堂在线| 精品久久蜜臀av无| 精品少妇内射三级| 老汉色∧v一级毛片| 日韩欧美一区视频在线观看| 国产有黄有色有爽视频| 欧美日韩一级在线毛片| 91久久精品国产一区二区三区| 超碰成人久久| 亚洲美女搞黄在线观看| 91精品三级在线观看| 亚洲激情五月婷婷啪啪| 在线观看三级黄色| 丰满少妇做爰视频| 亚洲精品国产av蜜桃| 久久精品国产亚洲av涩爱| 汤姆久久久久久久影院中文字幕| 午夜久久久在线观看| 欧美人与性动交α欧美软件| 丝瓜视频免费看黄片| 亚洲成人av在线免费| 国产精品秋霞免费鲁丝片| 日日摸夜夜添夜夜爱| 在线观看www视频免费| 成人午夜精彩视频在线观看| 久久精品国产亚洲av高清一级| 亚洲精品乱久久久久久| 精品久久久精品久久久| 高清在线视频一区二区三区| 久久久国产精品麻豆| 欧美精品一区二区免费开放| 一本—道久久a久久精品蜜桃钙片| 性色avwww在线观看| 久久久国产欧美日韩av| 日本免费在线观看一区| 亚洲av男天堂| 黄色毛片三级朝国网站| 熟妇人妻不卡中文字幕| 亚洲精品美女久久av网站| 少妇被粗大猛烈的视频| 日本av手机在线免费观看| 久久久久久免费高清国产稀缺| 日韩在线高清观看一区二区三区| 亚洲国产精品一区二区三区在线| 国产极品粉嫩免费观看在线| 美女xxoo啪啪120秒动态图| 狂野欧美激情性bbbbbb| 亚洲国产看品久久| 国产免费福利视频在线观看| 七月丁香在线播放| 王馨瑶露胸无遮挡在线观看| 晚上一个人看的免费电影| 伦精品一区二区三区| 中文字幕人妻熟女乱码| 激情五月婷婷亚洲| 色婷婷久久久亚洲欧美| 欧美av亚洲av综合av国产av | 精品人妻在线不人妻| 国产av精品麻豆| 成人国语在线视频| 国语对白做爰xxxⅹ性视频网站| 新久久久久国产一级毛片| 久久久久久久久久人人人人人人| av免费在线看不卡| 国产精品久久久久久精品电影小说| 久久这里有精品视频免费| 91午夜精品亚洲一区二区三区| 国产片内射在线| 国产精品免费视频内射| 久久精品国产亚洲av天美| 岛国毛片在线播放| 亚洲国产av新网站| 9191精品国产免费久久| 蜜桃国产av成人99| 精品少妇黑人巨大在线播放| 亚洲综合色惰| 色婷婷久久久亚洲欧美| 久久毛片免费看一区二区三区| 亚洲一码二码三码区别大吗| 自拍欧美九色日韩亚洲蝌蚪91| 免费日韩欧美在线观看| 欧美在线黄色| 国产精品 国内视频| 一级毛片电影观看| 麻豆乱淫一区二区| 深夜精品福利| 中文字幕人妻熟女乱码| 国产成人a∨麻豆精品| 亚洲图色成人| 久久久久久久久久人人人人人人| 五月伊人婷婷丁香| 黑人欧美特级aaaaaa片| 999精品在线视频| 亚洲,欧美,日韩| 久久精品国产自在天天线| 午夜福利视频精品| 亚洲av在线观看美女高潮| 国产精品av久久久久免费| 各种免费的搞黄视频| av女优亚洲男人天堂| 少妇被粗大的猛进出69影院| 日韩av在线免费看完整版不卡| 丰满迷人的少妇在线观看| 亚洲精品一区蜜桃| 90打野战视频偷拍视频| av片东京热男人的天堂| 2021少妇久久久久久久久久久| 精品视频人人做人人爽| 久久99精品国语久久久| 免费观看无遮挡的男女| 麻豆精品久久久久久蜜桃| 99国产精品免费福利视频| 日本wwww免费看| 侵犯人妻中文字幕一二三四区| 巨乳人妻的诱惑在线观看| 精品人妻偷拍中文字幕| 美女xxoo啪啪120秒动态图| 建设人人有责人人尽责人人享有的| 十八禁网站网址无遮挡| 亚洲国产av影院在线观看| 亚洲欧美清纯卡通| 亚洲人成网站在线观看播放| 两性夫妻黄色片| 大片免费播放器 马上看| 国产淫语在线视频| 午夜av观看不卡| 日韩免费高清中文字幕av| 亚洲视频免费观看视频| 日本午夜av视频| 人妻人人澡人人爽人人| 宅男免费午夜| 最新中文字幕久久久久| 久久久久网色| 捣出白浆h1v1| 亚洲av.av天堂| 狂野欧美激情性bbbbbb| 五月开心婷婷网| 国产亚洲欧美精品永久| 国产亚洲一区二区精品| 黑丝袜美女国产一区| 丰满乱子伦码专区| 久久久久久人妻| 91午夜精品亚洲一区二区三区| 中文字幕av电影在线播放| 日本91视频免费播放| 亚洲精品自拍成人| 啦啦啦中文免费视频观看日本| videos熟女内射| 午夜日本视频在线| 久久国产亚洲av麻豆专区| 尾随美女入室| av网站免费在线观看视频| 成人漫画全彩无遮挡| 久久精品国产a三级三级三级| 免费女性裸体啪啪无遮挡网站| 激情五月婷婷亚洲| 成年av动漫网址| 一级片免费观看大全| 男男h啪啪无遮挡| 成人午夜精彩视频在线观看| 少妇人妻 视频| 99热全是精品| 香蕉丝袜av| 天堂8中文在线网| 欧美97在线视频| 最黄视频免费看| 免费日韩欧美在线观看| 看免费成人av毛片| 一级毛片电影观看| 高清在线视频一区二区三区| 欧美日韩成人在线一区二区| 一本久久精品| 日本黄色日本黄色录像| 2018国产大陆天天弄谢| 中文字幕制服av| 亚洲精品日韩在线中文字幕| 国产1区2区3区精品| 一区福利在线观看| 国产麻豆69| 啦啦啦啦在线视频资源| 午夜福利在线观看免费完整高清在| 色视频在线一区二区三区| 欧美在线黄色| 日本午夜av视频| 亚洲欧美中文字幕日韩二区| 国产成人精品福利久久| 男女边吃奶边做爰视频| 久久精品国产亚洲av天美| 婷婷色麻豆天堂久久| 性色avwww在线观看| 国产毛片在线视频| 色吧在线观看| 欧美成人精品欧美一级黄| 精品福利永久在线观看| 精品少妇内射三级| 日韩在线高清观看一区二区三区| 七月丁香在线播放| 另类精品久久| 交换朋友夫妻互换小说| 大片免费播放器 马上看| 国产精品不卡视频一区二区| 爱豆传媒免费全集在线观看| 香蕉精品网在线| 男女啪啪激烈高潮av片| 亚洲国产最新在线播放| 人人妻人人爽人人添夜夜欢视频| 可以免费在线观看a视频的电影网站 | av有码第一页| av天堂久久9| 国产精品麻豆人妻色哟哟久久| 精品国产乱码久久久久久小说| 麻豆乱淫一区二区| 在线看a的网站| 精品卡一卡二卡四卡免费| 黄色一级大片看看| 亚洲欧洲精品一区二区精品久久久 | 视频在线观看一区二区三区| 我要看黄色一级片免费的| 国产免费现黄频在线看| 久久久久久久精品精品| 午夜免费男女啪啪视频观看| 亚洲av福利一区| 丁香六月天网| 少妇被粗大猛烈的视频| 狂野欧美激情性bbbbbb| 亚洲欧美一区二区三区黑人 | 18禁国产床啪视频网站| 国产精品人妻久久久影院| 亚洲图色成人| 波多野结衣一区麻豆| 国产日韩欧美视频二区| 久久人妻熟女aⅴ| 欧美日韩国产mv在线观看视频| 亚洲成av片中文字幕在线观看 | 久久 成人 亚洲| 亚洲精华国产精华液的使用体验| 少妇的丰满在线观看| 欧美成人精品欧美一级黄| 99国产综合亚洲精品| 9色porny在线观看| 久久这里有精品视频免费| 亚洲精品,欧美精品| 卡戴珊不雅视频在线播放| 国产伦理片在线播放av一区| 国产精品 欧美亚洲| 久久国产精品男人的天堂亚洲| 国产一区二区在线观看av| 久久久久久人人人人人| 免费高清在线观看视频在线观看| 麻豆乱淫一区二区| 国产极品天堂在线| 国产精品一二三区在线看| 美女大奶头黄色视频| 免费日韩欧美在线观看| 精品国产乱码久久久久久男人| av卡一久久| 女人久久www免费人成看片| 天美传媒精品一区二区| 国产伦理片在线播放av一区| 国产片内射在线| 成人午夜精彩视频在线观看| 又粗又硬又长又爽又黄的视频| 一个人免费看片子| 99热网站在线观看| 亚洲av在线观看美女高潮| 国产国语露脸激情在线看| 国产精品蜜桃在线观看| 亚洲熟女精品中文字幕| 丰满饥渴人妻一区二区三| 国产黄色免费在线视频| 最近2019中文字幕mv第一页| 欧美精品人与动牲交sv欧美| 美女国产高潮福利片在线看| 满18在线观看网站| 9热在线视频观看99| 午夜福利网站1000一区二区三区| 最近手机中文字幕大全| 久久精品人人爽人人爽视色| 日韩欧美精品免费久久| 精品国产露脸久久av麻豆| 日韩欧美精品免费久久| 色网站视频免费| 一区二区三区乱码不卡18| 叶爱在线成人免费视频播放| 人妻人人澡人人爽人人| 免费高清在线观看日韩| 人妻人人澡人人爽人人| 中文字幕另类日韩欧美亚洲嫩草| 夜夜骑夜夜射夜夜干| 久久人妻熟女aⅴ| 欧美国产精品va在线观看不卡| 久久99一区二区三区| a 毛片基地| 一级,二级,三级黄色视频| 麻豆精品久久久久久蜜桃| 国产精品无大码| 久久久久久人人人人人| 日本wwww免费看| 国产亚洲av片在线观看秒播厂| 成年女人毛片免费观看观看9 | 超碰97精品在线观看| 观看av在线不卡| 国产探花极品一区二区| 欧美日韩视频高清一区二区三区二| 在线精品无人区一区二区三| xxxhd国产人妻xxx| 免费大片黄手机在线观看| 99久久综合免费| 大陆偷拍与自拍|