亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        動態(tài)心電監(jiān)測系統(tǒng)介紹及關(guān)鍵技術(shù)水平分析

        2017-01-18 07:57:32王吉鳴呂穎瑩包濤樓理綱馮靖祎
        中國醫(yī)療設(shè)備 2016年10期
        關(guān)鍵詞:心電電信號特征提取

        王吉鳴,呂穎瑩,包濤,樓理綱,馮靖祎

        浙江大學(xué)醫(yī)學(xué)院附屬第一醫(yī)院 醫(yī)學(xué)工程部,浙江 杭州 310003

        動態(tài)心電監(jiān)測系統(tǒng)介紹及關(guān)鍵技術(shù)水平分析

        王吉鳴,呂穎瑩,包濤,樓理綱,馮靖祎

        浙江大學(xué)醫(yī)學(xué)院附屬第一醫(yī)院 醫(yī)學(xué)工程部,浙江 杭州 310003

        心臟系統(tǒng)疾病的防治和診斷己成為當(dāng)今醫(yī)學(xué)界面臨的重要問題,動態(tài)心電信號的監(jiān)測是解決該問題的主要手段。本文首先在心電信號分析理論的基礎(chǔ)上,介紹動態(tài)心電系統(tǒng)的硬件組成,對心電電極及硬件電路研究結(jié)果進(jìn)行總結(jié);其次,重點闡述了現(xiàn)階段QRS波自動檢測技術(shù),舉例說明運動偽跡干擾的濾波技術(shù);再次,綜述了QRS波特征提取的國內(nèi)外相關(guān)研究成果;最后,對比了基于心電信號數(shù)據(jù)和心電信號特征參數(shù)的質(zhì)量評估方法,并列舉了國內(nèi)外的心電信號質(zhì)量評估案例。

        動態(tài)心電圖;電極;QRS檢波;運動偽跡;質(zhì)量評估

        心血管疾病是嚴(yán)重威脅人類健康的常見疾病,心臟疾病的預(yù)防和治療問題亟待解決[1]。目前主要的心電監(jiān)測系統(tǒng)包括靜態(tài)心電監(jiān)測系統(tǒng)[2]、心率失常檢測、運動心電系統(tǒng)[3]及動態(tài)心電監(jiān)測系統(tǒng)(Holter系統(tǒng))。前3種方法的應(yīng)用受患者的運動狀態(tài)以及關(guān)鍵參數(shù)(波形振幅、頻率等)的限制,缺乏靈活性及全面性。動態(tài)心電監(jiān)測系統(tǒng)通過24小時監(jiān)測患者日常生活中心電的變化情況,可以全面直觀的了解患者心電變化,所以被廣泛應(yīng)用于目前的心電圖監(jiān)測中[4]。該系統(tǒng)根據(jù)心電信號的生理形成原因,可通過檢測P波、QRS波、T波的不同波段以及間期組成,分析異常特性診斷患者的病理[5-6]。

        動態(tài)心電信號具備信號幅值小、頻率低、差異性大、抗干擾能力弱等特點,其極易受到肌電噪聲、工頻噪聲、環(huán)境噪聲等因素影響[7]。所以設(shè)計的智能心電監(jiān)測系統(tǒng)需要克服以上不良影響,得到最真實的心電信號,并在相關(guān)專業(yè)軟件基礎(chǔ)上分析心電信號的關(guān)鍵點和異常點。智能心電監(jiān)測系統(tǒng)主要由硬件和軟件組成,硬件主要負(fù)責(zé)心電信號的采集處理與記錄,軟件主要負(fù)責(zé)心電信號的分析與結(jié)果比對。整個動態(tài)心電監(jiān)測流程,見圖1。本文總結(jié)了智能心電監(jiān)測系統(tǒng)在軟硬件研發(fā)過程中涉及到的相關(guān)技術(shù),對相關(guān)技術(shù)進(jìn)行分析及展望。

        1 動態(tài)心電監(jiān)測系統(tǒng)硬件介紹

        圖1 動態(tài)心電監(jiān)測流程

        動態(tài)心電監(jiān)測系統(tǒng)的硬件主要由兩個部分組成:心電傳感器和信號采集模塊。心電傳感器即心電電極是本系統(tǒng)的重要組成部分,其性能直接影響心電源信號的質(zhì)量?,F(xiàn)階段廣泛應(yīng)用的心電電極是一次性銀-氯化銀濕性電極,因銀-氯化銀材料的導(dǎo)電性優(yōu)越,且信號在該材料上損耗小,該電極在獲取心電源信號上具有獨特優(yōu)勢。但長期接觸銀-氯化銀會導(dǎo)致毒性滲入皮膚,隨著動態(tài)心電監(jiān)測系統(tǒng)品質(zhì)要求的提升,逐漸涌現(xiàn)出一批新型材料心電電極。干電極的出現(xiàn)克服了濕電極的部分缺點,北卡羅萊納州立大學(xué)的朱勇等人發(fā)明了一種基于銀納米線的干電極,該電極較傳統(tǒng)干電極的優(yōu)勢明顯,其信號采集質(zhì)量與濕電極性能相當(dāng),且其不采用凝膠層結(jié)構(gòu),提升了患者佩戴舒適感[8]。由于濕電極的性能優(yōu)越,研究者并沒有放棄該類型電極的研究,Lobodzinski等[9]應(yīng)用生物電位光纖傳感器解決了濕電極與人體長期接觸引起的皮膚敏感問題。由于材料技術(shù)的發(fā)展,出現(xiàn)了電子織物心電電極,其將導(dǎo)電電極紡線編織入衣物內(nèi),患者只要穿著該類型衣物,便可以通過鑲嵌在衣物中的電極實時采集心電信號。較早完整地實現(xiàn)該技術(shù)的是韓國高科技研究院,他們提出的時尚平面電路板P-FCB,將電極集成到衣物上進(jìn)行動態(tài)心電監(jiān)測[10]。

        動態(tài)心電監(jiān)測系統(tǒng)的信號采集模塊,見圖2。其基本原理是通過電極采集心電源信號,由于動態(tài)心電監(jiān)測信號的復(fù)雜特性,采集到的源信號幅值?。? mV左右)、頻率低(0.1~100 Hz)、噪聲大。噪聲的來源主要有環(huán)境熱噪聲、工頻干擾等,另外由于人體處于運動狀態(tài),動作會造成電極的變形、移位、擠壓,種種因素均會導(dǎo)致采集到的心電源信號不純凈,工程上稱之為運動偽跡干擾。所以需要對源信號進(jìn)行信號調(diào)理,主要的調(diào)理步驟包括前端放大、低通濾波、工頻陷波、后端放大。采用多級放大電路主要為了克服噪聲信號的放大;低通濾波器將不屬于心電信號的噪聲濾除,一般以二階及以上低通濾波為宜;工頻陷波電路可以有效濾除50 Hz的工頻信號,一般采用雙T型陷波電路。將調(diào)理后的信號傳輸至微處理器,現(xiàn)階段常用的微處理器主要有MSP430、ARM芯片等,該類處理器的功耗低,在電池系統(tǒng)中表現(xiàn)出良好的功耗性能,且處理器自身帶有12位的ADC模塊,既能滿足心電信號的采集精度要求,又能節(jié)約硬件成本。由于動態(tài)心電監(jiān)測信號需要記錄24 h的心電信號,系統(tǒng)的數(shù)據(jù)量大,一般數(shù)據(jù)存儲模塊會采用SD卡,其存儲容量可以達(dá)到GB級別,且具備體積小、數(shù)據(jù)讀取速度快、數(shù)據(jù)不易丟失等優(yōu)點。

        圖2 動態(tài)心電監(jiān)測系統(tǒng)信號采集模塊

        采集到的心電信號可以通過有線和無線兩種形式傳送至上位機(jī)進(jìn)行數(shù)據(jù)分析。有線數(shù)據(jù)傳輸?shù)木窒拊谟谛枰颊咴谂宕鲃討B(tài)心電監(jiān)測系統(tǒng)24 h后再將所有數(shù)據(jù)導(dǎo)入上位機(jī),耗費的時間較長,但是其數(shù)據(jù)傳輸可靠性高,較少出現(xiàn)數(shù)據(jù)丟失的情況;無線數(shù)據(jù)傳輸可以克服空間時間的限制,做到實時數(shù)據(jù)上傳顯示,但無線傳輸可靠性較難保證,存在數(shù)據(jù)傳輸丟包現(xiàn)象,而心電信號對數(shù)據(jù)完整性要求較高,所以對無線數(shù)據(jù)傳輸提出了較高要求。2010年,Teo等人研發(fā)了基于SoC芯片的無線傳感器節(jié)點,將其應(yīng)用于動態(tài)心電監(jiān)測系統(tǒng)的數(shù)據(jù)傳輸上,較好的實現(xiàn)了心電信號的監(jiān)測[11]。同年,我國的張兢等人也利用無線傳感模塊,研發(fā)了一套功耗低、數(shù)據(jù)傳輸穩(wěn)定性良好的動態(tài)心電監(jiān)測系統(tǒng),并在實際應(yīng)用中證實了裝置的可靠性。

        2 QRS波自動檢測技術(shù)的介紹

        心電信號監(jiān)測需要軟件配合將心電信號中的有效信號提取出來,常用的抑制噪聲干擾的方法有數(shù)字濾波器法、小波濾波、自適應(yīng)濾波等。1996年,Levkov等提出了采用分段處理信號的方法將采集到的信號進(jìn)行線性和非線性段濾波[12]。國內(nèi)也逐步涌現(xiàn)出一批有效的濾波算法,Zhang等[13]提出了基于多級、多尺度形態(tài)學(xué)的濾波器,該方法對心電信號的去干擾效果明顯。目前,由于肌電干擾帶來的心電信號噪聲仍然沒有特殊的解決方法,所以心電信號去干擾研究仍然十分具有前景。

        2.1 運動偽跡干擾處理技術(shù)

        在對經(jīng)硬件調(diào)理后的心電信號進(jìn)行特征檢測之前還需要完成另外一項重要的工作,即運動偽跡去除。運動偽跡給心電信號帶來的干擾圖,見圖3。

        運動偽跡干擾的心電圖與普通異常的心電圖相似,若不加以區(qū)分極有可能給患者帶來錯誤的診斷。各學(xué)者在對運動偽跡信號干擾方面進(jìn)行了大量深入研究,提出了各種信號處理技術(shù),其中較成熟的信號處理包括自適應(yīng)濾波、卡爾曼濾波以及小波變換濾波。自適應(yīng)濾波最初是由Widrow B提出的,他將帶干擾的信號與參考信號進(jìn)行比較,將得到的比較結(jié)果作為依據(jù)改變?yōu)V波器的相關(guān)參數(shù),改進(jìn)濾波器對應(yīng)特定輸入信號的濾波特性,提高信號的信噪比[14]。另外也有相關(guān)學(xué)者應(yīng)用加速度傳感器記錄人體運動特征,作為最終的數(shù)據(jù)處理中對應(yīng)時間段的心電信號噪聲參考源[15]。

        圖3 帶運動偽跡干擾的心電聯(lián)導(dǎo)信號

        2.2 心電信號特征提取技術(shù)

        在完成心電信號的濾波之后,需要對信號波形的特征進(jìn)行提取,主要是提取并分析信號中的QRS波群,然后對QRS波峰峰值、面積、寬度、橫向偏移量、縱向偏移量以及信噪比等進(jìn)行檢測。對QRS波群的特征檢測效果的好壞是決定動態(tài)心電監(jiān)測系統(tǒng)產(chǎn)品關(guān)鍵,其直接影響患者的診斷結(jié)果,所以國內(nèi)外學(xué)者提出了眾多對QRS波群的特征檢測方法,其中主要包括基于數(shù)學(xué)模型的QRS波特征提取、基于濾波和閾值檢測的QRS波特特征提取、基于模板匹配的QRS波特特征提取、基于小波變換的QRS波特特征提取等。特征提取的主要流程,見圖4。

        圖4 QRS波特征提取的流程

        基于濾波和閾值檢測的方法復(fù)雜性較低且應(yīng)用簡單,在早期的動態(tài)心電監(jiān)測產(chǎn)品中應(yīng)用較為廣泛,但是由于方法的局限性以及信號的復(fù)雜性,對閾值的設(shè)定需要人為根據(jù)經(jīng)驗法設(shè)定,其信號特征提取效果有限。當(dāng)然早期學(xué)者也在此基礎(chǔ)上提出了較多的改進(jìn)閾值檢測方法,例如Suppappola等[16]提出了改進(jìn)閾值檢測方法,GritZali[17]提出了用閾值方法去分隔幅值,提取有效的信號特征。

        基于數(shù)學(xué)模型的QRS波特征提取則將QRS波群進(jìn)行數(shù)學(xué)模型的建立,根據(jù)數(shù)學(xué)模型來區(qū)分有效的QRS波,雖然效果顯著但是建模的步驟計算量大,而且模型建立對高頻噪聲特別敏感,在實際應(yīng)用存在缺陷[18-19]。

        根據(jù)小波變換方法展開的QRS波特征提取方法是在信號進(jìn)行時域和頻域分析,將心電信號分解為多個尺度進(jìn)行多維分析,其效果較好。國內(nèi)的學(xué)者在基于小波變換的波群特征提取方面取得較多的成果。例如梁崴巍等人應(yīng)用樣條小波的方法捕捉到了QRS波的始末點,針對確定區(qū)域內(nèi)的QRS波形分析對病患的診斷更具說服力[20]。李翠微等人也是應(yīng)用樣條小波的原理,但是他們在此基礎(chǔ)上分割出小波分析將噪聲信號與QRS波信號尺度,根據(jù)不同尺度的分析可以區(qū)分QRS波的始末點。也有學(xué)者將QRS波的特征進(jìn)行聚類分析,得到相應(yīng)的分析結(jié)果供病患研究[21]。

        隨著數(shù)學(xué)計算方法與處理器硬件的發(fā)展,基于模板匹配的QRS波特征提取被越來越多的應(yīng)用于實際心電監(jiān)測產(chǎn)品,該方法的基本原理就是處理后的心電信號與模板庫中的信號進(jìn)行比對,得到對應(yīng)的特征參數(shù),模板庫的大小決定了該方法的有效性,現(xiàn)階段的處理器內(nèi)存技術(shù)已經(jīng)可以完全滿足模板庫要求。

        3 動態(tài)心電信號質(zhì)量評估技術(shù)

        心電信號用不同的方式進(jìn)行硬件信號采集調(diào)理、軟件去噪、運動偽跡去除以及波形特征提取后,得到的信號也是千差萬別,要使研發(fā)的產(chǎn)品性能好、診斷正確率高,需要做心電信號質(zhì)量評估。質(zhì)量評估方案是近幾年提出的,其對心電監(jiān)測系統(tǒng)產(chǎn)品質(zhì)量的提高帶來了巨大的好處。

        心電信號的質(zhì)量評估方法總體上可以分為兩種:① 基于整體或者局部的信號數(shù)據(jù)進(jìn)行質(zhì)量評估;② 基于信號特征參數(shù)提取進(jìn)行質(zhì)量評估。前者主要涉及到數(shù)學(xué)計算方面的均方誤差、標(biāo)準(zhǔn)差、信噪比、峰值誤差等原理,通過對原始數(shù)據(jù)的計算,以及現(xiàn)有的質(zhì)量評估標(biāo)準(zhǔn),完成心電信號質(zhì)量的分類。后者主要在完成QRS波特征提取之后,對得到的波形特征參數(shù)R波檢測匹配度,并且對Q波、R波、T波重合率等進(jìn)行質(zhì)量評估。前者應(yīng)用數(shù)學(xué)方法可以實現(xiàn)簡單的心電信號質(zhì)量分類,但并不能分辨出異常信號或者干擾信號處于整體信號的位置;后者雖然能分別異常信號與干擾信號,但同時也受制于噪聲干擾,若前端的QRS波特征提取效果不佳,將會直接影響信號質(zhì)量的評估。

        部分學(xué)者提出了一種直接的方法,在完成QRS波的特征提取后,將得到的波形面積進(jìn)行分析,順序比較QRS波兩兩之間的面積,以直方圖直觀的顯示心電信號質(zhì)量問題。該方法雖然直觀形象,但是僅采用波形面積一項參數(shù)作為質(zhì)量評估基準(zhǔn),缺乏科學(xué)性。后面又出現(xiàn)了以心電信號的振幅時間、基線漂移振幅、最大振幅等參數(shù)作為評估基準(zhǔn),結(jié)合Mason-Likar導(dǎo)聯(lián)方法,綜合評價信號質(zhì)量,使得評估結(jié)果更具說服力[22]。國內(nèi)學(xué)者易曉霖[23]提出了一套動態(tài)心電信號質(zhì)量評估方法,跟上述方法不同的是兩者提出的依據(jù)參數(shù)不同,其將QRS波群特征提取后得到的R波檢測匹配度也作為參數(shù),結(jié)合QRS波群的功率譜密度比值,劃分不同質(zhì)量等級的心電信號,該方法可行性較高且效果良好。

        4 總結(jié)

        我國是人口大國,心血管疾病患者及潛在患者數(shù)量眾多,這對動態(tài)心電監(jiān)測系統(tǒng)提出很高的要求。心電監(jiān)測系統(tǒng)在未來應(yīng)該向全信息、多導(dǎo)聯(lián)、多功能方向發(fā)展,提高心電信號的監(jiān)測性能和運動偽跡干擾的濾除,進(jìn)一步提升系統(tǒng)的自動分析能力。

        [1] Hossain P,Kawar B,Nahas M.Obesity and diabetes in the developing world a growing challenge[J].N Engl J Med,2007, 356(3):213-215.

        [2] Pan J,Tompkins WJ.A real-time QRS detection algorithm[J]. IEEE Trans Biomed Eng,1985,32(3):230-236.

        [3] 白家蓮.動態(tài)心電信號分析系統(tǒng)[D].長沙:中南大學(xué),2009.

        [4] 于南詳,陳東義,夏侯士戟.可穿戴計算技術(shù)及其應(yīng)用的新發(fā)現(xiàn)[J].數(shù)字通信,2012,39(4):13-20.

        [5] Huigen E,Peper A,Grimbergen CA.Investigation into the origin of the noise of surface electrodes[J].Med and Biol Eng and Comput,2002,40(3):332-338.

        [6] 解鳳英,范秀華,梁朝霞,等.動態(tài)心電圖在小二心律失常與心率變異性診斷中的意義[J].現(xiàn)代生物醫(yī)學(xué)進(jìn)展,2015,15(36):7144-7146.

        [7] 郭繼鴻,張萍.動態(tài)心電圖學(xué)[M].北京:人民衛(wèi)生出版社,2003: 10-36.

        [8] 李肅義,楊美玲,丁梅,等.動態(tài)心電采集記錄系統(tǒng)的發(fā)展與展望[J].生物醫(yī)學(xué)工程雜志,2012,29(1):175-178.

        [9] Lobodzinski SM,Laks MM.Biopotential fiber sensor[J]. J Electrocardiol,2006,39(4):S41-S46.

        [10] Yoo J,Yan L,Lee S,et al.A wearable ECG acquisition system with compact planar-fashionable circuit oard-based shirt information technology in bio-medicine[C].IEEE Transaction on 2009,13(6):897-902.

        [11] Teo HT,Qian XB,Gopalakrishnan PK,et al.A 700-wireless sensor node SoC for continuous real-time health monitoring[J]. IEEE J Solid-State Circuits,2010,45(11):2292-2299.

        [12] Dotsinsky IA,Daskalov IK.Accuracy of 50Hz interference subtraction from an electrocardiogram[J].Med Biol Eng Comput, 1996,34(6):489-494.

        [13] Zhang F, Lian Y.QRS detection based on multi-scale mathematical morphology foe wearable ECG devices in body area networks[J]. IEEE Trans Biomed Circuits Syst.2009,3(4):220-228.

        [14] Romero,Inaki B,Torfinn B,et al.Motion artifact reduction in ambulatory ECG monitoring: an integrated system approach[A]. Proceedings of the 2nd Conference on Wireless Health, ACM[C]. 2011:1124-1135.

        [15] Sung WY,Se DM,Yong HY,et al.Adaptive motion artifacts reduction using 3-axis accelerometer in e-textile ECG measurement system[J]. J Med Syst,2008,32(2):101-106.

        [16] Suppappola S,Sun Y.Nonlinear transforms of ECG signals for digital QRS detection: a quantitative analysis[J].IEEE Trans on Biomed Eng,1994,41(4):197-400.

        [17] Gritzali F.Towards a generalized scheme for QRS detection in ECG waveforms[J].Signal processing,1988,15(2):183-192.

        [18] 楊凱,叢林,胡文東,等.基于BMD101的嵌入式無線心電監(jiān)測系統(tǒng)[J].電子技術(shù)應(yīng)用,2014,40(1):122-124.

        [19] 凌振寶,李妹妍,張銘,等.基于ZigBee技術(shù)的無導(dǎo)聯(lián)線動態(tài)心電監(jiān)測儀設(shè)計[J].電子技術(shù)應(yīng)用,2015,41(8):66-71.

        [20] 董春桃,隋偉,喻洪流.基于改進(jìn)的差分閾值法心臟遠(yuǎn)程監(jiān)測的實時與回顧性分析[J].中國醫(yī)學(xué)裝備,2014,(11):51-54.

        [21] 晏明軍.基于小波變換和FPGA的心電監(jiān)測系統(tǒng)設(shè)計[J].電子測試,2014,(20):35-39.

        [22] Langley P, Di Marco, LY,et al.An algorithm for assessment of Quality of ECGs acquired via Mobile telephones[J].Comput Cardiol,2011,38:281-284.

        [23] 易曉霖.可穿戴心電信號質(zhì)量綜合評估及節(jié)律分析系統(tǒng)[D].上海:東華大學(xué),2012.

        Introduction of Dynamic ECG Monitoring System and Analysis of Key Technical Level

        WANG Ji-ming, LV Ying-ying, BAO Tao, LOU Li-gang, FENG Jing-yi
        Department of Clinical Engineering, The First Affiliated Hospital of Medical School of Zhejiang University, Hangzhou Zhejiang 310003, China

        Heart system disease prevention and diagnosis has become the top issue in medical field. Dynamic electrocardiography (ECG) signal monitor was the main mean to solve this problem. On the basis of ECG signal analysis theory, this paper introduced the hardware components of the dynamic ECG monitoring system, investigated the status of ECG electrodes and hardware circuit at home and abroad. In the algorithm aspects, it also elaborated QRS wave automatic detection technology. The paper illustrated the f lter technology of motion artifact interference, studied the related research results of QRS potter character extraction. Finally, it compared different quality evaluation method based on ECG signal data and ECG signal characteristic parameters, illustrated different cases of ECG signal quality evaluation.

        dynamic ECG; electrode; QRS wave detection; motion artifact; quality evaluation

        R541

        A

        10.3969/j.issn.1674-1633.2016.10.021

        1674-1633(2016)10-0071-04

        2016-03-31

        2016-05-07

        浙江省醫(yī)藥衛(wèi)生科技計劃項目(201463043)。

        馮靖祎,研究員,主要研究方向為臨床醫(yī)學(xué)工程技術(shù)與管理。

        通訊作者郵箱:casper_feng@163.com

        猜你喜歡
        心電電信號特征提取
        基于聯(lián)合聚類分析的單通道腹部心電信號的胎心率提取
        心電向量圖診斷高血壓病左心室異常的臨床應(yīng)用
        基于非接觸式電極的心電監(jiān)測系統(tǒng)
        電子制作(2019年19期)2019-11-23 08:41:40
        穿戴式心電:發(fā)展歷程、核心技術(shù)與未來挑戰(zhàn)
        基于Code Composer Studio3.3完成對心電信號的去噪
        科技傳播(2019年24期)2019-06-15 09:29:28
        更正啟事
        基于Daubechies(dbN)的飛行器音頻特征提取
        電子制作(2018年19期)2018-11-14 02:37:08
        基于隨機(jī)森林的航天器電信號多分類識別方法
        Bagging RCSP腦電特征提取算法
        基于MED和循環(huán)域解調(diào)的多故障特征提取
        久久人人爽av亚洲精品| 1000部精品久久久久久久久| 日本一极品久久99精品| 成人aaa片一区国产精品| 人妖精品视频在线观看| 亚洲色图专区在线观看| 久久频精品99香蕉国产| 日射精情感性色视频| 国产成年女人特黄特色毛片免| 久久精品国内一区二区三区| 天天色天天操天天日天天射| 亚洲欧美中文v日韩v在线| 国产av无码专区亚洲av毛网站| 高清少妇一区二区三区| 亚洲av永久无码精品秋霞电影影院| 亚洲2022国产成人精品无码区| 亚洲一区不卡在线导航| 午夜精品久久久久久久久久久久| 国产亚洲精品品视频在线| 欧美最大胆的西西人体44| 亚洲一区二区刺激的视频| 爆乳午夜福利视频精品| 狠狠色噜噜狠狠狠狠米奇777| 亚洲一区二区三区自拍麻豆| 成人精品一级毛片| 日本a级片免费网站观看| 无码三级国产三级在线电影| 正在播放国产对白孕妇作爱| 无码a级毛片免费视频内谢| 91免费国产| 看国产黄大片在线观看| 日韩精品熟妇一区二区三区| 日韩久久免费精品视频| 亚洲欧洲高潮| 大肉大捧一进一出好爽视频| 一本色道久久综合亚洲精品不| 免费国产一级片内射老| 色狠狠一区二区三区香蕉| 精品熟女av中文字幕| 久久99国产精一区二区三区| 亚洲av日韩av综合aⅴxxx|