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

    導(dǎo)波式雷達(dá)物位計(jì)測(cè)量誤差分析與軟件修正

    2015-10-25 02:16:48徐科軍
    計(jì)量學(xué)報(bào) 2015年5期
    關(guān)鍵詞:物位導(dǎo)波寬度

    劉 錚, 魏 萌, 徐科軍

    (合肥工業(yè)大學(xué)電氣與自動(dòng)化工程學(xué)院,安徽 合肥 230009)

    導(dǎo)波式雷達(dá)物位計(jì)測(cè)量誤差分析與軟件修正

    劉 錚, 魏 萌, 徐科軍

    (合肥工業(yè)大學(xué)電氣與自動(dòng)化工程學(xué)院,安徽 合肥 230009)

    分析了導(dǎo)波雷達(dá)物位計(jì)測(cè)量誤差產(chǎn)生的原因,提出軟件修正方法。在硬件方面分析了發(fā)射脈沖寬度和等效時(shí)間采樣間隔的影響;在軟件方面分析了信號(hào)處理過(guò)程產(chǎn)生的粗大誤差、隨機(jī)誤差和系統(tǒng)誤差?;诂F(xiàn)有的硬件平臺(tái),提出一整套處理方法以減小測(cè)量誤差,包括結(jié)合回波特征參數(shù)去判斷物位回波,9點(diǎn)移動(dòng)平均濾波預(yù)處理,選取最小斜率點(diǎn)作為回波定位點(diǎn),對(duì)計(jì)算值進(jìn)行中位值平均濾波。水位測(cè)量實(shí)驗(yàn)結(jié)果表明,在理論誤差為2~4 cm情況下,實(shí)際測(cè)量誤差小于1 cm,這驗(yàn)證了軟件修正方法的有效性。

    計(jì)量學(xué);導(dǎo)波雷達(dá)物位計(jì);測(cè)量誤差分析;軟件修正

    1 引 言

    導(dǎo)波式雷達(dá)物位計(jì)作為新型物位測(cè)量?jī)x表,廣泛應(yīng)用于石油、化工領(lǐng)域中水類(lèi)、油類(lèi)和固料類(lèi)等的接觸式測(cè)量[1-3]。由于導(dǎo)波桿對(duì)發(fā)射脈沖的匯聚導(dǎo)向作用,使電磁波能量更大程度地導(dǎo)向測(cè)量介質(zhì),因此,具有適用于小介電常數(shù)液位測(cè)量、受粉塵和泡沫等影響小的特點(diǎn)。根據(jù)準(zhǔn)確度要求,雷達(dá)物位計(jì)可劃分為過(guò)程控制級(jí)和貿(mào)易結(jié)算級(jí),導(dǎo)波式雷達(dá)物位計(jì)屬于過(guò)程控制級(jí)。

    目前,國(guó)內(nèi)文獻(xiàn)大都側(cè)重于介紹導(dǎo)波雷達(dá)物位計(jì)的工作原理和實(shí)際應(yīng)用,研究測(cè)量準(zhǔn)確度問(wèn)題的文獻(xiàn)較少。國(guó)外某些文獻(xiàn)提出了提高測(cè)量準(zhǔn)確度方法,產(chǎn)品的測(cè)量誤差最小可至1 mm,而國(guó)內(nèi)產(chǎn)品測(cè)量誤差最小為1 cm。例如,使用不完全外部介電傳輸探測(cè)器的特制導(dǎo)波桿[4],結(jié)合參考曲線的熵濾波法[5],安裝反射器[6]等。然而,國(guó)外的方法需要加裝反射器和特定的導(dǎo)波桿或需要空罐參考曲線。

    本文首先從硬件和軟件兩方面分析導(dǎo)播雷達(dá)物位計(jì)測(cè)量誤差的來(lái)源,然后,從信號(hào)處理方法角度,提出修正導(dǎo)波式雷達(dá)物位計(jì)測(cè)量誤差的方法,以提高其測(cè)量準(zhǔn)確度。

    2 誤差分析

    2.1硬件誤差

    導(dǎo)波雷達(dá)物位計(jì)依據(jù)時(shí)域反射原理測(cè)量物位。雷達(dá)主機(jī)的時(shí)基電路產(chǎn)生脈沖波形。該信號(hào)分為兩路。一路經(jīng)前置觸發(fā)電路放大整形,變成幅度較大、前沿較抖的脈沖信號(hào)送至脈沖發(fā)射器;脈沖發(fā)射器通過(guò)發(fā)射天線沿導(dǎo)波桿發(fā)送周期性窄脈沖電磁波;電磁波在不同介質(zhì)的分界面上反射后形成一系列回波信號(hào),由接收天線接收后送往等效時(shí)間采樣電路。另一路到采樣脈沖發(fā)生器,產(chǎn)生采樣脈沖,以控制等效時(shí)間采樣電路。等效時(shí)間采樣電路對(duì)回波信號(hào)進(jìn)行采樣,再將其變成中低頻信號(hào),送至單片機(jī)[7]。

    導(dǎo)波雷達(dá)物位計(jì)硬件由脈沖發(fā)射電路、等效時(shí)間采樣電路和信號(hào)處理電路3個(gè)部分組成,如圖1所示。下面分析脈沖發(fā)射電路和等效時(shí)間采樣電路所產(chǎn)生的誤差。

    圖1 硬件原理框圖

    (1)脈沖發(fā)射電路

    根據(jù)雷達(dá)原理中雷達(dá)距離分辨率的定義可知,脈沖式雷達(dá)可分辨兩目標(biāo)間的最小距離為:

    式中:ΔR為距離分辨率;c為電磁波在真空中的傳播速度;τ為脈沖收發(fā)單元產(chǎn)生的發(fā)射脈沖的寬度。

    可見(jiàn),脈沖式雷達(dá)能分辨目標(biāo)間的最小距離與發(fā)射脈沖的寬度有關(guān)。例如,某國(guó)產(chǎn)導(dǎo)波雷達(dá)物位計(jì)發(fā)射脈沖的寬度為15 ns,若空氣中電磁波傳播速度按光速計(jì)算,則ΔR=2.25 m。即在雷達(dá)脈沖波的傳輸途徑中,間距小于2.25 m的兩目標(biāo)的反射回波將會(huì)產(chǎn)生混疊。接收的反射回波混疊程度越高,則對(duì)目標(biāo)回波的分辨造成的影響也就越大,從而影響了測(cè)量的準(zhǔn)確度。相關(guān)資料表明,德國(guó)E+H公司導(dǎo)波雷達(dá)物位計(jì)的發(fā)射脈沖寬度僅為0.8 ns。所以,發(fā)射脈沖的寬度會(huì)影響測(cè)量準(zhǔn)確度[8]。通過(guò)減小發(fā)射脈沖寬度可提高距離分辨率,進(jìn)而減弱回波信號(hào)的混疊程度,有利于測(cè)量準(zhǔn)確度的提高。發(fā)射脈沖寬度越小,所包含的頻率成分越多,因此反射回波波峰越陡,也有利于測(cè)量準(zhǔn)確度的提高。

    (2)等效時(shí)間采樣電路

    由于在測(cè)量過(guò)程中,脈沖信號(hào)的行程時(shí)間為ns量級(jí),如果通過(guò)實(shí)時(shí)采樣來(lái)直接測(cè)量時(shí)差,普通的單片機(jī)無(wú)法滿(mǎn)足要求,所以,一般依據(jù)等效時(shí)間采樣原理先對(duì)回波信號(hào)進(jìn)行采樣[9],將高速信號(hào)在時(shí)域上放大,擴(kuò)展成低速信號(hào)。等效時(shí)間采樣方法的示意圖如圖2所示。

    圖2 等效時(shí)間采樣原理示意圖

    設(shè)接收脈沖重復(fù)周期為T(mén),兩路脈沖間隔微小時(shí)間差為Δt,單片機(jī)采樣頻率為Fs,則等效時(shí)間采樣放大倍數(shù)K為[10]:

    單片機(jī)采集經(jīng)過(guò)等效時(shí)間采樣電路轉(zhuǎn)換和放大后的回波信號(hào),其相鄰采樣點(diǎn)之間所對(duì)應(yīng)的物位距離d為:

    可見(jiàn),提高等效時(shí)間采樣電路的放大倍數(shù)K或者提高單片機(jī)采樣頻率Fs可以提高采樣精度。由式(2)可知,要提高放大倍數(shù)K就要減小延遲時(shí)間Δt。例如,某國(guó)產(chǎn)導(dǎo)波雷達(dá)物位計(jì)采用447 kHz的發(fā)射脈沖頻率,Δt為10 ps,F(xiàn)s為50 kHz,所以,

    另外,Δt的穩(wěn)定性也會(huì)影響測(cè)量準(zhǔn)確度。若Δt不夠穩(wěn)定,會(huì)導(dǎo)致等效放大倍數(shù)K的變化,從而導(dǎo)致在時(shí)域上拓展的回波波形的傳播時(shí)間發(fā)生變化,這樣會(huì)影響回波時(shí)間的計(jì)算準(zhǔn)確度,從而影響測(cè)量準(zhǔn)確度。

    2.2軟件誤差

    (1)粗大誤差

    信號(hào)處理中的粗大誤差來(lái)源于物位回波判斷錯(cuò)誤和失波。前者主要是由信號(hào)處理方法造成,后者由硬件條件、環(huán)境等造成。

    物位回波為介質(zhì)表面反射回波。在采集的回波信號(hào)中,除物位回波外,還有由其它因素產(chǎn)生的反射回波,這些統(tǒng)稱(chēng)為噪聲。以物位回波信號(hào)的頻率作為參照,噪聲可分為高頻噪聲和物位回波同頻噪聲兩類(lèi)。其中,高頻噪聲對(duì)應(yīng)于白噪聲的高頻部分;因操作條件、電源干擾、儀器不穩(wěn)定、障礙物產(chǎn)生的尖峰或突變信號(hào)等,被看成信號(hào)同頻噪聲??梢?jiàn),若無(wú)法從各類(lèi)噪聲中區(qū)分物位回波,將造成對(duì)真實(shí)回波的誤判,從而導(dǎo)致粗大誤差。

    傳統(tǒng)方法以回波峰值點(diǎn)幅值最大[11,12]作為條件進(jìn)行判斷。然而,由于實(shí)際信號(hào)的復(fù)雜性,此方法條件單一,在出現(xiàn)干擾回波幅值較大的情況下會(huì)發(fā)生誤判,造成粗大誤差。

    由于射頻電路的不穩(wěn)定、介質(zhì)的強(qiáng)烈揮發(fā)或液面晃動(dòng)等造成反射回波波形特別弱或消失,即所謂的“失波”,將導(dǎo)致無(wú)法確定物位回波或錯(cuò)誤判斷。

    (2)隨機(jī)誤差

    在導(dǎo)波雷達(dá)物位計(jì)的信號(hào)處理中可以看出,即使同一液位的測(cè)量值也并不穩(wěn)定,呈現(xiàn)上下波動(dòng)的特點(diǎn),故存在隨機(jī)誤差。此外,在實(shí)驗(yàn)測(cè)量過(guò)程中,通過(guò)鋼尺或米尺讀取標(biāo)準(zhǔn)值,由于鋼尺或米尺的精度限制,使讀數(shù)中不可避免地存在隨機(jī)誤差。

    (3)系統(tǒng)誤差

    信號(hào)處理中的系統(tǒng)誤差主要存在于傳播時(shí)間計(jì)算過(guò)程中。其中,回波定位點(diǎn)的選取是傳播時(shí)間計(jì)算的關(guān)鍵。物位回波波形類(lèi)似鐘形曲線,具有一定的寬度,有峰值點(diǎn)、上升沿、下降沿。應(yīng)當(dāng)選擇幅度相對(duì)穩(wěn)定且隨距離變化敏感的位置作為回波定位點(diǎn)。通常選擇回波峰值點(diǎn)作為定位點(diǎn)。然而,由于信號(hào)本身分辨率不高,會(huì)導(dǎo)致即使是同一液位,峰值點(diǎn)也會(huì)出現(xiàn)1~2點(diǎn)的波動(dòng)。根據(jù)式(5)可知,將導(dǎo)致1.34~2.68 cm的誤差。此外,實(shí)際中會(huì)出現(xiàn)平頂峰波形情況,真實(shí)峰值點(diǎn)無(wú)法確定,導(dǎo)致在一段液位內(nèi)確定的峰值點(diǎn)相同,也造成誤差。

    標(biāo)準(zhǔn)值若通過(guò)鋼尺或米尺獲得,則測(cè)量環(huán)境的溫度會(huì)對(duì)尺子造成熱脹冷縮,帶來(lái)誤差。此外,讀數(shù)人員的讀數(shù)習(xí)慣也會(huì)導(dǎo)致系統(tǒng)誤差。

    2.3其他因素誤差

    電路本身的設(shè)計(jì)和元器件的質(zhì)量對(duì)輸出高頻信號(hào)的信噪比有直接影響。測(cè)量介質(zhì)的介電常數(shù)直接影響回波波形。機(jī)械安裝的影響主要在于障礙物干擾,障礙物產(chǎn)生反射信號(hào),容易導(dǎo)致對(duì)被測(cè)介質(zhì)反射信號(hào)的誤判,從而增大錯(cuò)誤測(cè)量的概率。

    3 誤差修正

    從軟件修正出發(fā),研究合適的信號(hào)處理方法,來(lái)提高導(dǎo)波雷達(dá)物位計(jì)的測(cè)量準(zhǔn)確度。

    3.1信號(hào)預(yù)處理

    回波信號(hào)中含有高頻噪聲,需要進(jìn)行濾波預(yù)處理。去除高頻噪聲的濾波器有以巴特沃斯低通濾波為代表的頻域?yàn)V波器和以數(shù)字平滑算法為代表的時(shí)域?yàn)V波器,例如:移動(dòng)平均濾波器[13-14]、中值濾波器等。與頻域?yàn)V波器相比,時(shí)域?yàn)V波具有算法簡(jiǎn)單、運(yùn)行速度快的優(yōu)點(diǎn),因而,特別適合信號(hào)的實(shí)時(shí)處理。為此,選用時(shí)域?yàn)V波器進(jìn)行預(yù)處理。對(duì)于移動(dòng)平均濾波器,濾波結(jié)果是移動(dòng)窗口內(nèi)數(shù)據(jù)加權(quán)求和的結(jié)果,且所有數(shù)據(jù)分配了相同的加權(quán)因子。不同窗口寬度具有不同的效果:窗口越寬,平滑效果越好,但也會(huì)使波形變化率變小,信號(hào)幅值降低,信號(hào)畸變大,因此,應(yīng)當(dāng)合理選擇移動(dòng)窗口寬度。通常選擇比半峰寬略小的寬度,由于物位回波的寬度不固定,從10點(diǎn)至20點(diǎn)都存在,因此,對(duì)5點(diǎn)、7點(diǎn)、9點(diǎn)移動(dòng)平均濾波效果進(jìn)行了對(duì)比。

    隨機(jī)選擇若干個(gè)液位;隨機(jī)保留各液位下的一組信號(hào),用來(lái)計(jì)算傳播時(shí)間。其中,傳播時(shí)間的計(jì)算以傳統(tǒng)方法的峰值點(diǎn)為準(zhǔn)。對(duì)液位-傳播時(shí)間數(shù)據(jù)進(jìn)行最小二乘擬合,并用均方根誤差評(píng)價(jià)其線性度,即作為濾波性能的評(píng)價(jià)指標(biāo)。式中:σ為均方根誤差;Xi為樣本值;為樣本均值;n為樣本數(shù)。計(jì)算結(jié)果如表1所示。

    表1 不同窗寬移動(dòng)平均濾波比較

    經(jīng)過(guò)移動(dòng)平均處理后,均方根誤差明顯減??;5點(diǎn)移動(dòng)平均與9點(diǎn)移動(dòng)平均濾波效果相近,略低于7點(diǎn)移動(dòng)平均。此外,考慮到越多點(diǎn)數(shù)的移動(dòng)平均濾波改善平頂峰的效果越好,如圖3和4所示,將有利于后續(xù)特征參數(shù)的計(jì)算、回波定位點(diǎn)的確定,因此,選擇9點(diǎn)移動(dòng)平均濾波。

    3.2物位回波識(shí)別

    為了避免粗大誤差,對(duì)物位回波識(shí)別方法進(jìn)行改進(jìn),提出基于多特征參數(shù)識(shí)別的物位回波處理方法,利用回波的多個(gè)參數(shù),通過(guò)不同方式的結(jié)合構(gòu)成多個(gè)判斷條件,增強(qiáng)了回波判斷的可靠性。

    圖3 9點(diǎn)移動(dòng)平均濾波效果

    圖4 平頂峰濾波后效果

    基于在實(shí)驗(yàn)室條件下采集的數(shù)據(jù),對(duì)表征回波特征的參數(shù)進(jìn)行統(tǒng)計(jì)分析,包括峰值點(diǎn)幅值、寬度、變化率等,結(jié)果表明,多數(shù)情況下物位回波具有最大幅值,寬度基本穩(wěn)定,但回波疊加會(huì)導(dǎo)致寬度變大;變化率較大,但隨時(shí)間沒(méi)有呈現(xiàn)規(guī)律性變化。不同干擾回波有不同特點(diǎn),難以通過(guò)單一條件來(lái)確定物位回波。為此,根據(jù)干擾回波的特點(diǎn)將其分類(lèi),構(gòu)成不同條件進(jìn)行排除。

    在回波判斷時(shí),第一步是通過(guò)設(shè)置回波寬度、變化率、下降比的下限來(lái)排除具有一個(gè)或多個(gè)參數(shù)與物位回波明顯不同的虛假回波。這些回波稱(chēng)之為第一類(lèi)回波,如電子干擾產(chǎn)生的隨機(jī)窄尖峰脈沖、傳感器本身結(jié)構(gòu)反射回波和寬度幅值很小的噪聲。第二步再采用綜合多個(gè)特征參數(shù)計(jì)算置信度的方法予以排除一個(gè)或多個(gè)特征參數(shù)和物位回波較為接近的干擾波,如二次或多次回波、干擾物產(chǎn)生的回波及寬度較大的噪聲。對(duì)特征與真實(shí)回波接近的干擾波特性進(jìn)行統(tǒng)計(jì)分析可知:寬度越小、峰值點(diǎn)幅值越大,變化率越大的回波為物位回波的可能性越大。因此,根據(jù)變化率、峰值點(diǎn)幅值、寬度來(lái)設(shè)置一個(gè)置信度來(lái)確定真實(shí)回波,其經(jīng)驗(yàn)公式為:置信度=0.6×變化率+0.1×峰值點(diǎn)幅值-0.3×寬度。其中的0.6、0.1和0.3是權(quán)值,可以根據(jù)實(shí)際情況進(jìn)行調(diào)整。

    該方法可以排除虛假回波,準(zhǔn)確判斷幅值較小的物位回波,分別以不同液位(15 cm和76 cm)的水的回波信號(hào)為例,說(shuō)明該方法的有效性,如圖5和6所示。圖中,一組△和○對(duì)應(yīng)一個(gè)經(jīng)過(guò)去除第一類(lèi)回波后的剩余回波。其中,△代表起點(diǎn),○代表結(jié)束點(diǎn);豎直虛線之間標(biāo)記了經(jīng)過(guò)置信度計(jì)算后最終確定的物位回波范圍。在圖5中,干擾波的特征與回波明顯不同,通過(guò)判斷回波寬度、變化率和下降比等特征就能將其排除,確定出真實(shí)物位回波。在圖6中,經(jīng)過(guò)第一步判斷,可以找到3個(gè)波,即圖中箭頭所指的物位回波和兩個(gè)干擾回波。其中,第一個(gè)干擾回波比物位回波幅值略高,變化率小,寬度大;第二個(gè)干擾回波比物位回波幅值低,變化率小,寬度大。將這3個(gè)波的特征參數(shù)(幅值、變化率、寬度)代入上述經(jīng)驗(yàn)公式計(jì)算各自置信度,物位回波的置信度最大,從而確定其為真實(shí)回波。

    圖5 15 cm水位回波信號(hào)

    3.3回波定位點(diǎn)選擇

    經(jīng)過(guò)研究,回波的特征點(diǎn)有峰值點(diǎn)、起點(diǎn)、終點(diǎn)、最大斜率點(diǎn)、最小斜率點(diǎn)等,其中,起點(diǎn)和終點(diǎn)往往變動(dòng)較大,已經(jīng)不能真實(shí)反映回波隨液位變化情況,可不考慮。

    在實(shí)驗(yàn)室條件下,選用5種回波定位點(diǎn)進(jìn)行比較,分別是:峰值點(diǎn)[15]、經(jīng)2階Lagrange插值的峰值點(diǎn)、最大斜率點(diǎn)、最小斜率點(diǎn),以及以這兩點(diǎn)作為邊界的區(qū)域中心點(diǎn)進(jìn)行了比較。中心點(diǎn)計(jì)算公式為:

    圖6 76 cm水位回波信號(hào)

    式中:xmid為區(qū)域中心點(diǎn)橫坐標(biāo);xL為左邊界點(diǎn);xR為右邊界點(diǎn);y為回波曲線;x為橫坐標(biāo)。

    先對(duì)20~100 cm液位范圍內(nèi)隨機(jī)抽取了27個(gè)液位,并在各個(gè)液位分別采集一組信號(hào),對(duì)信號(hào)進(jìn)行9點(diǎn)移動(dòng)平均濾波預(yù)處理后,分別用峰值點(diǎn)、峰值點(diǎn)2階Lagrange插值點(diǎn)、最大斜率點(diǎn)、最小斜率點(diǎn)、中心點(diǎn)作為回波定位點(diǎn),以各回波定位點(diǎn)確定的(物位-時(shí)間)線性模型均方根誤差作為衡量指標(biāo),如表2所示。

    表2 不同回波定位點(diǎn)比較

    由表2數(shù)據(jù)可知,用最小斜率點(diǎn)作為回波定位點(diǎn)進(jìn)行最小二乘擬合的均方根誤差最小,即線性度最好,用于計(jì)算物位的準(zhǔn)確度最高。

    3.4計(jì)算結(jié)果后處理

    在實(shí)際工作現(xiàn)場(chǎng)中,由于各種干擾的存在,即使物位恒定,傳播時(shí)間計(jì)算值仍不斷變化。在同一液位,連續(xù)采集5 000組的傳播時(shí)間計(jì)算值,如圖7所示??傻脗鞑r(shí)間的波動(dòng)區(qū)間為[233.441 9,234.605 9],最大偏差為1.164 0 ns,對(duì)應(yīng)換算成距離為1.560 9 cm,即同一物位測(cè)量值的最大偏差為1.560 9 cm。

    根據(jù)概率密度的計(jì)算方法,將該液位下傳播時(shí)間的樣本波動(dòng)區(qū)間分為30段,對(duì)樣本進(jìn)行排序后計(jì)算在區(qū)間上的概率密度,如圖8所示。

    圖7 傳播時(shí)間樣本序列

    圖8 傳播時(shí)間樣本概率密度函數(shù)

    由圖8可知,傳播時(shí)間樣本的概率密度函數(shù)不滿(mǎn)足高斯分布。故對(duì)傳播時(shí)間采用平均值濾波的方法不能有效減小噪聲。為此,選擇中值濾波[16]。對(duì)傳播時(shí)間的樣本序列進(jìn)行中值濾波后再計(jì)算概率密度函數(shù),并對(duì)其進(jìn)行高斯擬合,其結(jié)果見(jiàn)圖9。

    圖9 傳播時(shí)間樣本概率密度函數(shù)高斯擬合

    由圖9可看出,經(jīng)中值濾波后的樣本概率密度符合高斯分布,將濾波后的樣本再進(jìn)行平均濾波,可以有效減小隨機(jī)誤差。隨機(jī)選擇30組液位記錄液位和對(duì)應(yīng)傳播時(shí)間計(jì)算值,所得數(shù)據(jù)通過(guò)最小二乘擬合確定均方根誤差,如表3所示。通過(guò)對(duì)比表2和表3可知,進(jìn)行兩級(jí)濾波有助于減小誤差,提高計(jì)算準(zhǔn)確度。

    表3 兩級(jí)濾波后物位-傳播時(shí)間模型的均方根誤差

    4 實(shí)驗(yàn)驗(yàn)證

    為了驗(yàn)證本文所提方法的有效性,在實(shí)驗(yàn)室進(jìn)行了水位測(cè)量實(shí)驗(yàn)。導(dǎo)波桿長(zhǎng)103 cm,采樣頻率為50 kHz。分別以峰值點(diǎn)的2階Lagrange插值點(diǎn)和最小斜率點(diǎn)作為回波定位點(diǎn)進(jìn)行對(duì)比實(shí)驗(yàn)。為了考核的嚴(yán)格性,選取能夠反映最大測(cè)量誤差的部分測(cè)量結(jié)果,如表4所示。

    表4 測(cè)量結(jié)果比較cm

    可見(jiàn),以峰值點(diǎn)的2階Lagrange插值點(diǎn)為回波定位點(diǎn)的測(cè)量誤差大于1 cm,以最小斜率點(diǎn)為回波定位點(diǎn)的測(cè)量誤差在1 cm以?xún)?nèi),與理論誤差相比均有較大減小。

    5 結(jié) 論

    為了提高導(dǎo)波雷達(dá)物位計(jì)測(cè)量準(zhǔn)確度,從硬件和軟件兩方面對(duì)其誤差來(lái)源進(jìn)行分析,在信號(hào)處理方法上對(duì)影響準(zhǔn)確度的因素進(jìn)行修正。

    (1)在硬件方面,發(fā)射脈沖寬度大、等效時(shí)間采樣電路的時(shí)間間隔大、單片機(jī)的采樣頻率低是導(dǎo)致導(dǎo)波式雷達(dá)物位計(jì)測(cè)量準(zhǔn)確度不高的原因。需要改進(jìn)硬件平臺(tái),如減小發(fā)射脈沖寬度、提高采樣頻率、減小等效時(shí)間采樣電路的時(shí)間間隔等,以提高測(cè)量的準(zhǔn)確度。

    (2)在信號(hào)處理過(guò)程中,僅根據(jù)最大峰值來(lái)判斷有效回波信號(hào)是造成導(dǎo)波雷達(dá)物位計(jì)測(cè)量不可靠、測(cè)量盲區(qū)大的主要原因之一。為此,基于現(xiàn)有硬件條件的限制,提出結(jié)合回波特征參數(shù),構(gòu)成多個(gè)條件來(lái)排除虛假回波、確定物位回波的方法,提高物位回波判斷的可靠性,有效地排除干擾,減小粗大誤差。

    (3)針對(duì)不可避免的噪聲干擾導(dǎo)致信號(hào)抖動(dòng)、毛刺現(xiàn)象等,提出采用9點(diǎn)移動(dòng)平均濾波對(duì)信號(hào)進(jìn)行預(yù)處理,以及將最小斜率點(diǎn)作為回波定位點(diǎn)進(jìn)行傳播時(shí)間計(jì)算;針對(duì)計(jì)算值的波動(dòng)問(wèn)題,根據(jù)其不符合高斯分布的特點(diǎn),提出采用中位值平均濾波減小隨機(jī)誤差。實(shí)驗(yàn)結(jié)果表明,在理論誤差為2~4 cm的情況下,實(shí)際測(cè)量誤差小于1 cm,有效地提高了測(cè)量準(zhǔn)確度。

    [1] Cataldo A,Tarricone L,Attivissimo F,et al.A TDR method for real-time monitoring of liquids[J].IEEE Trans Instrum Meas,2007,56(5):1616-1625.

    [2] 張寧,黃楊,李建雙.關(guān)于物位計(jì)測(cè)量方法的研究[J].計(jì)量學(xué)報(bào),2008,29(z1):104-106.

    [3] 張寧,黃楊.關(guān)于雷達(dá)物位計(jì)實(shí)際應(yīng)用及其測(cè)量方法的探討[J].計(jì)量學(xué)報(bào),2009,30(z1):111-114.

    [4] 羅斯蒙特儲(chǔ)罐雷達(dá)股份公司.通過(guò)多模式傳播的具有介電常數(shù)補(bǔ)償?shù)膶?dǎo)波雷達(dá)物位計(jì)系統(tǒng):中國(guó),CN201120366374.7[P].2012-06-06.

    [5] Wall G D.Method for an entropy filter echo processing in time-of-flightorlevelmeasurementsystems:U.S,6935177B2[P].2005.

    [6] Delin H,Edvardsson O.System and method for filling level determination:U.S,7525476 B1[P].2009.

    [7] 蘭茂亮,趙輝,劉偉佳.導(dǎo)波雷達(dá)液位回波信號(hào)處理電路的設(shè)計(jì)[J].計(jì)量技術(shù),2010,(10):14-17.

    [8] 張龍浩,祝懷標(biāo),譙艷.CPLD和MSP430單片機(jī)在導(dǎo)波雷達(dá)物位計(jì)中的應(yīng)用[J].電子設(shè)計(jì)工程,2012,(21):141-144.

    [9] 劉偉佳,趙輝,唐美芹.一種導(dǎo)波雷達(dá)液位計(jì)的設(shè)計(jì)與實(shí)現(xiàn)[J].自動(dòng)化與儀表,2009,(11):14-17.

    [10] 姜波,陰麗娟,陳祥光,等.等效時(shí)間采樣在雷達(dá)液位儀中的應(yīng)用[J].儀表技術(shù)與傳感器,2005,(6):50-54.

    [11] 趙建文,侯媛彬.信號(hào)距離度模型及其波形復(fù)雜信號(hào)辨識(shí)[J].電子測(cè)量與儀器學(xué)報(bào),2013,27(6):485-491.

    [12] 李存龍,陳偉民,章鵬.正弦調(diào)制型微波測(cè)距系統(tǒng)的不確定度分析[J].儀器儀表學(xué)報(bào),2014,35(1):23 -29.

    [13] 王穎,金志軍.常用數(shù)字濾波算法[J].中國(guó)計(jì)量,2012,(3):99-100.

    [14] 李聚剛.雷達(dá)物位計(jì)的應(yīng)用研究及算法改進(jìn)[D].上海:上海交通大學(xué),2008.

    [15] 馬遠(yuǎn),魏萌,徐科軍.基于反向回波的導(dǎo)波雷達(dá)物位計(jì)信號(hào)處理方法[J].化工自動(dòng)化及儀表,2013,(12):1473-1476.

    [16] 李建安,楊學(xué)志,高雋,等.兩類(lèi)自適應(yīng)非線性濾波器設(shè)計(jì)方法研究[J].儀器儀表學(xué)報(bào),2005,26(8):374-377.

    Measurement Error Analysis and Software Correction of Guided Wave Radar Level Gauge

    LIU Zheng, WEI Meng, XU Ke-Jun
    (School of Electrical and Automation Engineering,Hefei University of Technology,Hefei,Anhui 230009,China)

    The causes of the error of the guided wave radar level gauge(GWRLG)are analyzed,and a software correction method is proposed.The impacts of transmit pulse width and equivalent time sampling interval are studied in hardware.The gross error,the random error and the systematic error generated in the signal processing are researched in software.Based on the existing hardware platform,a set of signal processing method is put forward to reduce the measurement error,i.e.,combining feature parameters to recognize the actual level echo,performing 9 point moving average to preprocess the signal,selecting the minimum slope point as the echo positing point,and conducting median average filtering for the calculated values.The experimental results of water level measurement show that the actual measurement error is reduced to within 1 cm when the theoretical error is 2~4 cm,which validates the effectiveness of the proposed software correction method.

    metrology;guided wave radar level gauge;measurement error analysis;software correction

    TB938.1

    A

    1000-1158(2015)05-0506-07

    10.3969/j.issn.1000-1158.2015.05.13

    2014-08-27;

    2015-03-16

    劉錚(1990-),男,陜西西安人,合肥工業(yè)大學(xué)碩士研究生,研究方向?yàn)榍度胧较到y(tǒng)。292231698@qq.com徐科軍為通訊作者。dsplab@hfut.edu.cn

    猜你喜歡
    物位導(dǎo)波寬度
    超聲導(dǎo)波技術(shù)在長(zhǎng)輸管道跨越段腐蝕檢測(cè)中的應(yīng)用
    一種改進(jìn)的寬波束線陣?yán)走_(dá)物位測(cè)量算法
    深圳計(jì)為針對(duì)不同工況為用戶(hù)提供相應(yīng)物位測(cè)量解決方案
    卷簧缺陷檢測(cè)的超聲導(dǎo)波傳感器研制
    電子制作(2019年9期)2019-05-30 09:42:00
    馬屁股的寬度
    紅細(xì)胞分布寬度與血栓的關(guān)系
    孩子成長(zhǎng)中,對(duì)寬度的追求更重要
    人生十六七(2015年5期)2015-02-28 13:08:24
    磁致伸縮導(dǎo)波激勵(lì)傳感器模型及輸出特性
    Endress+Hauser Micropilot FMR5X雷達(dá)物位儀
    20~1000kHz超聲導(dǎo)波用于錨桿錨固質(zhì)量檢測(cè)的探索
    金屬礦山(2013年11期)2013-03-11 16:55:12
    成人特级黄色片久久久久久久| 天堂动漫精品| 女人被躁到高潮嗷嗷叫费观| 亚洲国产欧美日韩在线播放| 欧美黑人欧美精品刺激| 国产精品九九99| 男人的好看免费观看在线视频 | 国产黄a三级三级三级人| 人成视频在线观看免费观看| 国产精品 欧美亚洲| 国产熟女xx| 国产一区二区三区视频了| 国产蜜桃级精品一区二区三区| 90打野战视频偷拍视频| 久9热在线精品视频| 级片在线观看| 亚洲熟女毛片儿| 久久精品aⅴ一区二区三区四区| 国产精品日韩av在线免费观看 | 国产成人欧美| 国产精品免费一区二区三区在线| 两性夫妻黄色片| www.精华液| 好看av亚洲va欧美ⅴa在| 日韩免费av在线播放| 久热这里只有精品99| 免费观看人在逋| 国产精品秋霞免费鲁丝片| 一级毛片女人18水好多| 午夜福利高清视频| 国产亚洲欧美在线一区二区| 中文字幕人妻熟女乱码| 男男h啪啪无遮挡| 欧美亚洲日本最大视频资源| 国产一区二区在线av高清观看| 黄片大片在线免费观看| 日韩大尺度精品在线看网址 | 美女 人体艺术 gogo| 50天的宝宝边吃奶边哭怎么回事| 国产精品1区2区在线观看.| 美女免费视频网站| 久久亚洲真实| 国产aⅴ精品一区二区三区波| 国产乱人伦免费视频| 99re在线观看精品视频| 午夜福利在线观看吧| 欧美另类亚洲清纯唯美| 亚洲欧美精品综合一区二区三区| 亚洲成人国产一区在线观看| 丝袜人妻中文字幕| 法律面前人人平等表现在哪些方面| 久久草成人影院| 丁香欧美五月| 国产精品久久久av美女十八| 身体一侧抽搐| 欧美在线黄色| 深夜精品福利| 久久中文字幕人妻熟女| 久久久久久久精品吃奶| 18美女黄网站色大片免费观看| 后天国语完整版免费观看| 国产麻豆69| 亚洲欧美激情综合另类| 可以免费在线观看a视频的电影网站| 亚洲国产中文字幕在线视频| 国产精品野战在线观看| 女人爽到高潮嗷嗷叫在线视频| 欧美激情久久久久久爽电影 | 麻豆久久精品国产亚洲av| 99在线视频只有这里精品首页| 少妇的丰满在线观看| 亚洲avbb在线观看| 怎么达到女性高潮| 搡老岳熟女国产| 午夜视频精品福利| 亚洲美女黄片视频| 美女大奶头视频| 成人三级黄色视频| 中文字幕久久专区| 亚洲色图综合在线观看| 伦理电影免费视频| 天堂√8在线中文| 俄罗斯特黄特色一大片| 中亚洲国语对白在线视频| av在线播放免费不卡| 自线自在国产av| 一进一出好大好爽视频| 亚洲五月天丁香| 日本vs欧美在线观看视频| 久久久久久久久免费视频了| 精品国产国语对白av| 亚洲三区欧美一区| 激情视频va一区二区三区| 男人舔女人的私密视频| 国产精品一区二区三区四区久久 | 国产亚洲精品av在线| 久久久久久免费高清国产稀缺| 他把我摸到了高潮在线观看| 国产成人精品久久二区二区91| 成人国语在线视频| 啦啦啦韩国在线观看视频| 欧美乱码精品一区二区三区| 曰老女人黄片| 曰老女人黄片| 国产男靠女视频免费网站| 国产精品九九99| 欧美日韩乱码在线| 国产亚洲欧美98| 国内久久婷婷六月综合欲色啪| 国产亚洲精品一区二区www| 精品国产超薄肉色丝袜足j| 国产精品国产高清国产av| 日韩成人在线观看一区二区三区| 亚洲熟妇熟女久久| 午夜亚洲福利在线播放| 一进一出抽搐动态| 9色porny在线观看| 黄色成人免费大全| 亚洲欧美激情在线| 国产精品综合久久久久久久免费 | 国产精品日韩av在线免费观看 | 亚洲精品中文字幕在线视频| 国产又色又爽无遮挡免费看| 国产一区二区激情短视频| 亚洲av成人av| 亚洲人成77777在线视频| 一区二区三区精品91| 成人18禁在线播放| 少妇熟女aⅴ在线视频| 日韩精品青青久久久久久| 动漫黄色视频在线观看| 波多野结衣av一区二区av| 久久人人精品亚洲av| 久久天堂一区二区三区四区| 黑丝袜美女国产一区| 一级作爱视频免费观看| 国产av一区在线观看免费| 搞女人的毛片| 亚洲 欧美 日韩 在线 免费| 国产xxxxx性猛交| 国产精品九九99| 国产国语露脸激情在线看| 国产人伦9x9x在线观看| 97人妻天天添夜夜摸| 日韩 欧美 亚洲 中文字幕| 天堂动漫精品| 免费高清在线观看日韩| 久久久久国内视频| 日韩一卡2卡3卡4卡2021年| 咕卡用的链子| 99国产极品粉嫩在线观看| 亚洲精品国产精品久久久不卡| 午夜福利免费观看在线| 亚洲精品中文字幕一二三四区| bbb黄色大片| 淫妇啪啪啪对白视频| 欧美激情久久久久久爽电影 | 国产精品久久视频播放| 一卡2卡三卡四卡精品乱码亚洲| 又黄又爽又免费观看的视频| 美女免费视频网站| 90打野战视频偷拍视频| 日韩av在线大香蕉| 国产麻豆69| 丁香六月欧美| 午夜精品国产一区二区电影| 国产一区二区三区综合在线观看| 大香蕉久久成人网| 色综合欧美亚洲国产小说| 亚洲国产欧美网| 国产精品99久久99久久久不卡| 亚洲精品在线观看二区| videosex国产| 日韩一卡2卡3卡4卡2021年| 亚洲中文日韩欧美视频| 又紧又爽又黄一区二区| 亚洲美女黄片视频| av有码第一页| 国产精品精品国产色婷婷| 极品人妻少妇av视频| 亚洲精品av麻豆狂野| 日韩免费av在线播放| 国产一区二区在线av高清观看| 国产成人av教育| 亚洲av片天天在线观看| 老司机靠b影院| 国产一级毛片七仙女欲春2 | 日本欧美视频一区| 黄频高清免费视频| 欧美丝袜亚洲另类 | 欧美在线一区亚洲| 色播亚洲综合网| 少妇裸体淫交视频免费看高清 | 99精品在免费线老司机午夜| 午夜视频精品福利| 国产精品二区激情视频| av超薄肉色丝袜交足视频| 可以免费在线观看a视频的电影网站| 亚洲一区二区三区不卡视频| 99国产极品粉嫩在线观看| 日韩欧美一区二区三区在线观看| 黑丝袜美女国产一区| 最近最新中文字幕大全免费视频| 久久性视频一级片| 精品久久久久久成人av| 香蕉久久夜色| 久久狼人影院| 亚洲五月婷婷丁香| 亚洲在线自拍视频| 神马国产精品三级电影在线观看 | 999精品在线视频| 久久久水蜜桃国产精品网| 无人区码免费观看不卡| 日日爽夜夜爽网站| 亚洲全国av大片| 亚洲精品美女久久久久99蜜臀| tocl精华| 亚洲天堂国产精品一区在线| 亚洲国产欧美日韩在线播放| 国产欧美日韩一区二区三| 十八禁人妻一区二区| 啦啦啦观看免费观看视频高清 | 亚洲五月婷婷丁香| 黄片播放在线免费| 日本免费一区二区三区高清不卡 | 黑人操中国人逼视频| 免费在线观看影片大全网站| 99国产精品免费福利视频| 中文字幕人妻熟女乱码| 男男h啪啪无遮挡| 国产成人欧美| 亚洲男人天堂网一区| 欧美色视频一区免费| 日日夜夜操网爽| 少妇 在线观看| 精品人妻在线不人妻| 日本a在线网址| 国产精品亚洲av一区麻豆| 亚洲色图综合在线观看| 久久国产精品影院| 午夜福利,免费看| 69精品国产乱码久久久| 精品一区二区三区av网在线观看| 国产成人欧美| 18禁国产床啪视频网站| 午夜精品在线福利| 精品一区二区三区av网在线观看| 999久久久精品免费观看国产| 欧美激情极品国产一区二区三区| 好男人在线观看高清免费视频 | 丝袜美足系列| 女人高潮潮喷娇喘18禁视频| 美女高潮到喷水免费观看| 三级毛片av免费| netflix在线观看网站| 亚洲成人免费电影在线观看| 欧美黑人精品巨大| 色精品久久人妻99蜜桃| 天堂影院成人在线观看| 精品福利观看| 国产欧美日韩一区二区三| 桃红色精品国产亚洲av| 午夜免费鲁丝| 嫩草影视91久久| 大型av网站在线播放| 午夜福利成人在线免费观看| 美女国产高潮福利片在线看| 黑丝袜美女国产一区| 久久人人爽av亚洲精品天堂| 午夜a级毛片| 黄色片一级片一级黄色片| 又紧又爽又黄一区二区| 久久亚洲真实| 韩国av一区二区三区四区| 欧美黑人欧美精品刺激| 天堂影院成人在线观看| 午夜a级毛片| 桃色一区二区三区在线观看| 999久久久精品免费观看国产| 亚洲aⅴ乱码一区二区在线播放 | 日日摸夜夜添夜夜添小说| 色综合亚洲欧美另类图片| 亚洲成a人片在线一区二区| 999久久久精品免费观看国产| 亚洲黑人精品在线| 又黄又粗又硬又大视频| 精品久久久久久久毛片微露脸| 91老司机精品| 男女做爰动态图高潮gif福利片 | 亚洲伊人色综图| 禁无遮挡网站| 成年版毛片免费区| 亚洲七黄色美女视频| 中文字幕高清在线视频| 91精品三级在线观看| www.精华液| 久久久国产欧美日韩av| 极品人妻少妇av视频| 久久久久久免费高清国产稀缺| 69av精品久久久久久| 丝袜美足系列| 久久狼人影院| 高清毛片免费观看视频网站| 国产熟女xx| 国产成年人精品一区二区| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲 欧美 日韩 在线 免费| 嫩草影院精品99| 国产激情久久老熟女| 久久精品亚洲熟妇少妇任你| 午夜福利欧美成人| 亚洲五月天丁香| 国内精品久久久久精免费| 好男人电影高清在线观看| 亚洲精品久久成人aⅴ小说| 亚洲精品久久国产高清桃花| 伦理电影免费视频| 国产精品国产高清国产av| 免费高清视频大片| 亚洲一区高清亚洲精品| 国产精品一区二区精品视频观看| 99热只有精品国产| 99久久国产精品久久久| 老熟妇乱子伦视频在线观看| 99精品欧美一区二区三区四区| 国产成人av教育| 老熟妇乱子伦视频在线观看| 黄色丝袜av网址大全| 精品一区二区三区视频在线观看免费| 极品教师在线免费播放| 国产亚洲欧美精品永久| 亚洲中文字幕日韩| 视频区欧美日本亚洲| 精品午夜福利视频在线观看一区| 九色亚洲精品在线播放| 国产欧美日韩精品亚洲av| 99久久精品国产亚洲精品| 侵犯人妻中文字幕一二三四区| 国产成年人精品一区二区| 性少妇av在线| 国产99久久九九免费精品| 制服人妻中文乱码| 少妇粗大呻吟视频| 一夜夜www| 久久久久精品国产欧美久久久| 精品国产超薄肉色丝袜足j| 国产一区二区在线av高清观看| 精品久久蜜臀av无| 在线国产一区二区在线| av欧美777| 精品一区二区三区四区五区乱码| 一卡2卡三卡四卡精品乱码亚洲| 九色亚洲精品在线播放| 午夜福利免费观看在线| 免费久久久久久久精品成人欧美视频| 国语自产精品视频在线第100页| 久久香蕉国产精品| 日韩欧美在线二视频| 18禁观看日本| 91av网站免费观看| 最新美女视频免费是黄的| 不卡一级毛片| 狠狠狠狠99中文字幕| 中文字幕久久专区| 男女下面进入的视频免费午夜 | 男女做爰动态图高潮gif福利片 | www.精华液| 成年版毛片免费区| 美女高潮到喷水免费观看| 精品午夜福利视频在线观看一区| 亚洲欧美激情在线| 亚洲性夜色夜夜综合| 法律面前人人平等表现在哪些方面| av免费在线观看网站| 国产精品爽爽va在线观看网站 | 亚洲国产欧美一区二区综合| 欧美久久黑人一区二区| 神马国产精品三级电影在线观看 | 老司机靠b影院| 亚洲精品中文字幕一二三四区| 国产99久久九九免费精品| 一区二区日韩欧美中文字幕| 国产成人一区二区三区免费视频网站| 国产亚洲精品久久久久久毛片| 18禁美女被吸乳视频| 久久 成人 亚洲| 亚洲成人国产一区在线观看| 极品人妻少妇av视频| 亚洲熟女毛片儿| 一夜夜www| 夜夜躁狠狠躁天天躁| 在线视频色国产色| 免费看十八禁软件| 久久青草综合色| av欧美777| 人人妻人人澡人人看| 久久国产精品影院| 叶爱在线成人免费视频播放| 久久久久久久久免费视频了| 国产成人系列免费观看| 国产成人av教育| 国产精品影院久久| 久久草成人影院| 精品国产一区二区三区四区第35| 多毛熟女@视频| 国产av一区二区精品久久| 日本在线视频免费播放| 91麻豆精品激情在线观看国产| 91大片在线观看| 亚洲av美国av| 日日夜夜操网爽| 视频区欧美日本亚洲| 午夜日韩欧美国产| 搡老妇女老女人老熟妇| 色老头精品视频在线观看| 老鸭窝网址在线观看| 黄色成人免费大全| 成人手机av| 久久久久久免费高清国产稀缺| 好男人在线观看高清免费视频 | 国产精品久久久久久精品电影 | 亚洲国产高清在线一区二区三 | 国产97色在线日韩免费| 日本欧美视频一区| 亚洲男人的天堂狠狠| 少妇被粗大的猛进出69影院| 怎么达到女性高潮| 色综合站精品国产| 国产亚洲精品一区二区www| 男女之事视频高清在线观看| 亚洲国产精品合色在线| 欧美成人一区二区免费高清观看 | 日本撒尿小便嘘嘘汇集6| 人人妻人人澡欧美一区二区 | 日韩欧美三级三区| 日日夜夜操网爽| 亚洲人成77777在线视频| 国产亚洲精品av在线| 男女之事视频高清在线观看| 一二三四在线观看免费中文在| 桃红色精品国产亚洲av| 99久久久亚洲精品蜜臀av| 操出白浆在线播放| 色在线成人网| 日本三级黄在线观看| 精品乱码久久久久久99久播| 精品高清国产在线一区| 18禁美女被吸乳视频| 777久久人妻少妇嫩草av网站| 免费看十八禁软件| 国产精品综合久久久久久久免费 | 悠悠久久av| 亚洲人成电影免费在线| 亚洲欧美日韩无卡精品| 欧美中文综合在线视频| 亚洲久久久国产精品| 黄色毛片三级朝国网站| 久久中文字幕一级| 黄色成人免费大全| 99香蕉大伊视频| 久久久久亚洲av毛片大全| 女人被狂操c到高潮| 日韩欧美一区视频在线观看| 亚洲情色 制服丝袜| 美女高潮到喷水免费观看| 香蕉久久夜色| 色精品久久人妻99蜜桃| 国产一区二区三区视频了| 三级毛片av免费| 美女大奶头视频| 国内久久婷婷六月综合欲色啪| 男女床上黄色一级片免费看| netflix在线观看网站| 999久久久精品免费观看国产| 嫩草影院精品99| 欧美国产精品va在线观看不卡| av网站免费在线观看视频| 亚洲 欧美 日韩 在线 免费| 日本一区二区免费在线视频| 欧美在线黄色| 国产在线精品亚洲第一网站| 90打野战视频偷拍视频| 手机成人av网站| 午夜福利视频1000在线观看 | 可以免费在线观看a视频的电影网站| 男女做爰动态图高潮gif福利片 | 久久精品国产清高在天天线| 欧美最黄视频在线播放免费| 精品免费久久久久久久清纯| 日韩中文字幕欧美一区二区| 黄色片一级片一级黄色片| 国产精品久久久av美女十八| 亚洲情色 制服丝袜| 国产成人精品在线电影| 免费无遮挡裸体视频| 后天国语完整版免费观看| 午夜免费成人在线视频| 多毛熟女@视频| 长腿黑丝高跟| 无限看片的www在线观看| 1024视频免费在线观看| 欧美一区二区精品小视频在线| 国产精品亚洲av一区麻豆| 国产欧美日韩一区二区三| 黄色毛片三级朝国网站| 最近最新中文字幕大全免费视频| 国产一区二区三区视频了| 99久久精品国产亚洲精品| 久久亚洲精品不卡| 亚洲欧美日韩高清在线视频| 免费观看人在逋| 国产精品亚洲一级av第二区| 欧美色视频一区免费| 亚洲成人精品中文字幕电影| 中文字幕最新亚洲高清| 国产成人精品在线电影| 九色亚洲精品在线播放| 国产国语露脸激情在线看| 中文字幕人成人乱码亚洲影| 怎么达到女性高潮| 精品免费久久久久久久清纯| 精品国产超薄肉色丝袜足j| 操美女的视频在线观看| 久久精品国产亚洲av香蕉五月| 日日摸夜夜添夜夜添小说| 欧美成人性av电影在线观看| 999精品在线视频| 欧美中文日本在线观看视频| 男女床上黄色一级片免费看| 欧美乱色亚洲激情| 禁无遮挡网站| 亚洲欧美日韩无卡精品| 黑人欧美特级aaaaaa片| 日韩欧美免费精品| 亚洲 欧美 日韩 在线 免费| 国产区一区二久久| 1024香蕉在线观看| 亚洲一区二区三区不卡视频| 大香蕉久久成人网| 无遮挡黄片免费观看| 91大片在线观看| 国产精品精品国产色婷婷| 日韩精品中文字幕看吧| 亚洲一卡2卡3卡4卡5卡精品中文| 中文字幕人成人乱码亚洲影| 精品国内亚洲2022精品成人| 日日干狠狠操夜夜爽| 亚洲中文字幕一区二区三区有码在线看 | 亚洲av电影不卡..在线观看| 亚洲av成人av| 欧美激情久久久久久爽电影 | 亚洲国产精品sss在线观看| 中文字幕最新亚洲高清| 国产av精品麻豆| 国产精品亚洲一级av第二区| 夜夜夜夜夜久久久久| 99热只有精品国产| 国产免费男女视频| 男人的好看免费观看在线视频 | 日本免费一区二区三区高清不卡 | 久久久水蜜桃国产精品网| 久久久久国产精品人妻aⅴ院| 一区在线观看完整版| 国产真人三级小视频在线观看| 精品国内亚洲2022精品成人| 色播亚洲综合网| 好男人电影高清在线观看| svipshipincom国产片| or卡值多少钱| 黄网站色视频无遮挡免费观看| 欧美丝袜亚洲另类 | 九色国产91popny在线| 久久中文字幕一级| 国产午夜福利久久久久久| 大码成人一级视频| 国产一区二区在线av高清观看| 老司机在亚洲福利影院| 非洲黑人性xxxx精品又粗又长| 国产精品久久久人人做人人爽| 咕卡用的链子| 两性夫妻黄色片| 中文字幕人成人乱码亚洲影| 免费一级毛片在线播放高清视频 | 一区二区三区高清视频在线| 色av中文字幕| 琪琪午夜伦伦电影理论片6080| 亚洲五月色婷婷综合| av片东京热男人的天堂| 亚洲成av片中文字幕在线观看| 色老头精品视频在线观看| 午夜老司机福利片| 黄片小视频在线播放| 亚洲va日本ⅴa欧美va伊人久久| 亚洲一卡2卡3卡4卡5卡精品中文| 婷婷六月久久综合丁香| 精品卡一卡二卡四卡免费| 午夜福利视频1000在线观看 | 亚洲情色 制服丝袜| 久久精品国产亚洲av高清一级| 午夜精品国产一区二区电影| 在线播放国产精品三级| 99国产综合亚洲精品| 麻豆久久精品国产亚洲av| 法律面前人人平等表现在哪些方面| 成人免费观看视频高清| 黄网站色视频无遮挡免费观看| 精品欧美一区二区三区在线| 亚洲中文av在线| 又紧又爽又黄一区二区| 久久久国产成人精品二区| 19禁男女啪啪无遮挡网站| 亚洲欧美日韩另类电影网站| www国产在线视频色| 国产三级黄色录像| 亚洲人成网站在线播放欧美日韩|