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

        ?

        數(shù)字化能譜獲取中梯形成形研究

        2017-10-13 12:55:02劉寅宇王玉東周榮楊朝文
        核技術(shù) 2017年2期
        關(guān)鍵詞:虛線(xiàn)能譜梯形

        劉寅宇 王玉東 周榮 楊朝文

        ?

        數(shù)字化能譜獲取中梯形成形研究

        劉寅宇 王玉東 周榮 楊朝文

        (四川大學(xué)物理科學(xué)與技術(shù)學(xué)院輻射物理與技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室 成都 610064)

        在數(shù)字化能譜采集系統(tǒng)中,往往需要添加適當(dāng)?shù)拿}沖成形算法對(duì)數(shù)字化核信號(hào)進(jìn)行處理,以提升系統(tǒng)能量分辨率等指標(biāo)。梯形成形算法因具有實(shí)現(xiàn)簡(jiǎn)單、成形脈沖窄等諸多優(yōu)勢(shì)而得到廣泛的應(yīng)用。從探測(cè)器輸出的脈沖信號(hào)進(jìn)行數(shù)學(xué)模型建立及簡(jiǎn)化出發(fā),對(duì)信號(hào)模型的梯形成形算法進(jìn)行分析討論,并利用Simulink在現(xiàn)場(chǎng)可編程邏輯門(mén)陣列(Field-Programmable Gate Array, FPGA)上實(shí)現(xiàn)并改進(jìn)梯形成形算法,通過(guò)實(shí)驗(yàn)測(cè)試了成形算法的可行性及其在能譜測(cè)量過(guò)程中對(duì)能量分辨率的提升。

        數(shù)字譜儀,脈沖模型,梯形成形,現(xiàn)場(chǎng)可編程邏輯門(mén)陣列實(shí)現(xiàn),能量分辨率

        20世紀(jì)90年代以來(lái),隨著各種數(shù)字化器件制造及應(yīng)用技術(shù)的快速發(fā)展,核信號(hào)能譜測(cè)量技術(shù)逐步向數(shù)字化方向發(fā)展[1?4]。相對(duì)于傳統(tǒng)模擬成形能譜測(cè)量多道,數(shù)字化能譜獲取處理方法可以取消線(xiàn)性放大器、采樣保持、成形濾波等傳統(tǒng)的模擬電路系統(tǒng),而直接通過(guò)高采樣率模擬數(shù)字轉(zhuǎn)換器(Analog-to-digital Converter, ADC)對(duì)光電倍增管(Photomultiplier Tubes, PMT)的輸出信號(hào)進(jìn)行實(shí)時(shí)采樣,并利用現(xiàn)場(chǎng)可編程門(mén)陣列(Field- Programmable Gate Array, FPGA)等數(shù)字信號(hào)處理芯片對(duì)數(shù)字化的核信號(hào)進(jìn)行實(shí)時(shí)處理,以獲取能譜,這種方案在價(jià)格和體積等方面有著顯著的優(yōu)勢(shì)[5?8]。

        數(shù)字化能譜獲取方案中采用的ADC器件的有效位數(shù)(Effective Number Of Bits, ENOB)對(duì)系統(tǒng)能量分辨率有直接的影響:假設(shè)ADC的有效位數(shù)為,參考電壓為,則其量化誤差為/2,因此有效位數(shù)越高量化誤差越小。理論上ADC有效位數(shù)越高越有利于能量分辨率提升,然而ADC的有效位數(shù)越高價(jià)格越貴,因此在核信號(hào)處理中涌現(xiàn)出各種降低ADC有效位數(shù)需求的處理方法[9?12],例如高斯成形、Sallen-key濾波成形和梯形成形等。

        梯形成形算法是其中一種較常用的方法,對(duì)比其他成形算法,該成形算法具有在FPGA上實(shí)現(xiàn)更加簡(jiǎn)單、脈寬窄以及參數(shù)調(diào)節(jié)更加靈活等優(yōu)點(diǎn)。本文從探測(cè)器輸出信號(hào)數(shù)學(xué)模型建立出發(fā),對(duì)不同情況下的探測(cè)器輸出信號(hào)進(jìn)行建模分析,對(duì)不同輸出信號(hào)模型的梯形成形算法進(jìn)行討論分析,并利用MATLAB的Simulink工具在FPGA進(jìn)行梯形成形算法實(shí)現(xiàn),最終通過(guò)實(shí)驗(yàn)進(jìn)行驗(yàn)證測(cè)試。

        1 探測(cè)器輸出脈沖信號(hào)模型

        閃爍事件在光電倍增管陽(yáng)極產(chǎn)生的電壓脈沖形狀取決于陽(yáng)極電路的時(shí)間常數(shù),對(duì)于圖1所示的經(jīng)典的光電倍增管輸出電路,其輸出信號(hào)的數(shù)學(xué)模型如下[13]:

        式中:是碘化鈉晶體的發(fā)光衰減時(shí)間;00為電路放電時(shí)間常數(shù)。

        圖1 光電倍增管信號(hào)輸出電路

        Fig.1 PMT signal output circuit.

        當(dāng)00>>時(shí):

        當(dāng)00<<時(shí):

        在進(jìn)行數(shù)字化核信號(hào)處理時(shí),我們更關(guān)心脈沖信號(hào)的時(shí)間特性,因此可對(duì)式(2)、(3)統(tǒng)一簡(jiǎn)化得到雙指數(shù)衰減信號(hào)模型式(4):

        當(dāng)滿(mǎn)足1>>2的條件時(shí),可對(duì)于上述模型進(jìn)一步簡(jiǎn)化得到理想單指數(shù)衰減信號(hào)模型式(5):

        利用MATLAB對(duì)兩種情況下的輸出脈沖信號(hào)進(jìn)行模擬,得到如圖2所示結(jié)果。其中實(shí)線(xiàn)是雙指數(shù)衰減脈沖(Double Exponential, DE)信號(hào),虛線(xiàn)是單指數(shù)衰減脈沖(Single Exponential, SE)信號(hào)。模擬中為了簡(jiǎn)化模擬過(guò)程,假設(shè)輸出相對(duì)幅度相等。

        圖2 MATLAB模擬探測(cè)器輸出信號(hào)

        2 梯形成形算法建立

        對(duì)于圖3所示的理想等腰梯形信號(hào),其時(shí)域表達(dá)式為:

        其中:

        將式(7)?(10)代入式(6)中,并將式(6)離散化后做變換得:

        式中:n=t/、n=t/、n=t/,為ADC采樣周期。

        圖3 等腰梯形信號(hào)

        Fig.3 Isosceles trapezoid signal.

        對(duì)光電倍增管輸出信號(hào)式(4)、(5)作時(shí)域上的離散化處理,并進(jìn)行變換得:

        對(duì)變換后的式(11)、(13)作商,可得單指數(shù)脈沖梯形(Single Exponential Trapezoidal, SET)成形算法:

        SET型算法級(jí)聯(lián)結(jié)構(gòu)如圖4所示,級(jí)聯(lián)構(gòu)造中取消了除法運(yùn)算1/n,以便于算法在FPGA上實(shí)現(xiàn),而將其移至成形算法后。

        圖4 單指數(shù)脈沖梯形成形算法級(jí)聯(lián)結(jié)構(gòu)

        Fig.4 SET filter cascade ladder.

        同理,針對(duì)雙指數(shù)衰減信號(hào),利用式(11)、(12)得雙指數(shù)脈沖梯形(Double Exponential Trapezoidal, DET)成形算法:

        DET型算法級(jí)聯(lián)結(jié)構(gòu)如圖5所示,結(jié)構(gòu)中同樣取消了1/[(1?2)n]的除法運(yùn)算。

        圖5 雙指數(shù)脈沖梯形成形算法級(jí)聯(lián)結(jié)構(gòu)

        Fig.5 DET filter cascade ladder.

        為了驗(yàn)證梯形成形算法的成形效果,利用MATLAB對(duì)以上兩種算法進(jìn)行模擬,得到圖6所示效果:對(duì)于單指數(shù)衰減信號(hào)(圖6(a)中帶點(diǎn)虛線(xiàn)),SET型濾波算法能夠成形出理想的梯形(圖6(a)中短虛線(xiàn));對(duì)于雙指數(shù)衰減信號(hào)(圖6(a)中實(shí)線(xiàn)),SET型濾波成形算法無(wú)法對(duì)其成形為理想梯形(圖6(a)中長(zhǎng)虛線(xiàn));而DET型成形算法能夠?qū)﹄p指數(shù)衰減信號(hào)成形為理想的梯形(圖6(b)中虛線(xiàn))。

        圖6 SET (a)和DET (b)算法成形效果模擬

        SET和DET成形算法的建立均立足于信號(hào)幅度不變化這一基本條件之上,因此粒子在探測(cè)器中沉積能量與成形后脈沖信號(hào)的幅度呈正比關(guān)系。

        3 梯形成形算法的幅頻響應(yīng)

        雖然成形算法的設(shè)計(jì)著眼于信號(hào)的時(shí)域特征,但是兩種成形算法處理在頻域中都相當(dāng)于對(duì)信號(hào)進(jìn)行濾波操作。為了進(jìn)一步分析討論兩種濾波算法在頻域中的濾波效果,對(duì)式(5)、(6)作傅里葉變換,得SET型算法的頻率響應(yīng)函數(shù):

        式中:為分量信號(hào)角頻率。對(duì)式(16)整理取模得:

        同理,對(duì)式(4)、(6)做傅里葉變換,并對(duì)得到的DET型頻率響應(yīng)函數(shù)整理取模得:

        圖7(a)和(c)分別為SET和DET型算法頻率響應(yīng)曲線(xiàn),其中:t=1 μs,t=2 μs,1=40 μs,2=100 μs;圖7(b)為t=2 μs、t=3 μs時(shí)的SET型算法頻率響應(yīng)曲線(xiàn)。

        由圖7可知,SET型成形算法在高頻噪聲抑制方面效果比DET型算法優(yōu)秀;t取值越大,振蕩周期越短,信號(hào)低頻成分相對(duì)增加,對(duì)高頻噪聲抑制效果越好,但是脈沖變寬,增加堆積概率,實(shí)際應(yīng)用中應(yīng)根據(jù)具體情況權(quán)衡選擇。

        針對(duì)上述影響壓裂效果的因素,提出優(yōu)化地質(zhì)選井選層、優(yōu)化壓裂施工規(guī)模、優(yōu)化射孔方式、優(yōu)化壓裂液、加強(qiáng)現(xiàn)場(chǎng)施工管理等措施,減少壓裂波及鄰井情況。2017年優(yōu)化工藝以來(lái)施工30口井78層,均未出現(xiàn)壓竄鄰井現(xiàn)象,工藝成功率100%。

        圖7 SET (a,b)和DET (c)成形算頻率響應(yīng)曲線(xiàn)

        4 成形算法的FPGA實(shí)現(xiàn)

        為了將成形算法應(yīng)用到實(shí)際工程中,利用MATLAB的Simulink工具中的Altera DSP Builder系列模塊搭建梯形成形算法如圖8所示。其中黑色粗線(xiàn)框內(nèi)為SET型算法,在SET型輸出后添加一個(gè)1?2?1參數(shù)級(jí)構(gòu)成完整的DET型算法。

        1?1?1中1是一個(gè)小于1的小數(shù),如果直接采用乘法器對(duì)1?1進(jìn)行運(yùn)算,會(huì)導(dǎo)致大量的FPGA資源浪費(fèi)和性能消耗。為了避免以上情況出現(xiàn),本文統(tǒng)一對(duì)參數(shù)1倍增4096倍(可根據(jù)需求選擇不同倍增因子):用兩個(gè)乘法器Product 1和Product 2對(duì)Adder 4兩邊輸入同時(shí)倍增4096倍;采用兩個(gè)乘法器保證了Adder 4兩端輸入時(shí)序上的同步。本方案避免了除法運(yùn)算和浮點(diǎn)數(shù)運(yùn)算,并嚴(yán)格保證算法內(nèi)部時(shí)序上同步。因此該方案的運(yùn)行時(shí)鐘頻率得到提升,該方案可在Altera公司經(jīng)濟(jì)型FPGA芯片EP4CE30F23C8N上及其100 MHz時(shí)鐘頻率下穩(wěn)定運(yùn)行。

        完成以上算法搭建后,利用Simulink中Signal Compiler工具編譯生成知識(shí)產(chǎn)權(quán)(Intellectual property, IP)核;將生成的IP核添加到Quartus II編寫(xiě)的工程中;最終由Quartus II完成程序的綜合、映射、布線(xiàn)及下載工作。

        圖8 DET成形算法Simulink框圖

        5 實(shí)驗(yàn)驗(yàn)證測(cè)量

        為了驗(yàn)證梯形成形在能譜測(cè)量中對(duì)能量分辨率的影響,本文利用圖9所示的系統(tǒng)對(duì)137Cs放射源進(jìn)行能譜測(cè)量:系統(tǒng)中探測(cè)晶體采用常用的碘化鈉晶體;適配電路(Adapter Circuit)部分采用亞德諾半導(dǎo)體(Analog Devices Inc, ADI)公司的差分放大器AD8352;ADC部分采用ADI公司的AD9434(有效位數(shù)為10 bit)芯片;數(shù)字信號(hào)處理部分采用Altera公司的EP4CE30F23C8N經(jīng)濟(jì)型FPGA芯片作為數(shù)字信號(hào)處理核心芯片。

        對(duì)于NaI(TI)探測(cè)器,利用性能優(yōu)越的NIM能譜分析系統(tǒng)測(cè)量137Cs能譜時(shí),其全能峰分辨率可達(dá)到8%以?xún)?nèi)。圖9是系統(tǒng)沒(méi)有NIM能譜分析系統(tǒng)中的前放、主放及濾波成形等信號(hào)處理電路,因此在直接進(jìn)行能譜測(cè)量時(shí)其能量分辨率較差。

        圖9 實(shí)驗(yàn)系統(tǒng)框圖

        在實(shí)驗(yàn)中抓取了梯形成形前后脈沖信號(hào)如圖10所示。圖10(a)中實(shí)線(xiàn)為高速ADC直接采樣后的脈沖信號(hào)(即原始核脈沖);虛線(xiàn)為經(jīng)過(guò)SET型成形算法成形后的脈沖信號(hào);圖10(b)中虛線(xiàn)為經(jīng)過(guò)DET型成形算法成形后的脈沖信號(hào)。由圖10可知,經(jīng)過(guò)DET型成形算法成形后的脈沖比SET型的脈沖更加接近梯形,其頂部更加平坦,而SET型算法成形后脈沖更加光滑。

        圖10 SET (a)和DET (b)算法對(duì)實(shí)際信號(hào)成形效果

        未經(jīng)過(guò)成形處理的原始脈沖信號(hào)測(cè)量能譜(圖11實(shí)線(xiàn))的137Cs全能峰分辨率為10.78%,經(jīng)過(guò)SET型成形算法成形后的脈沖信號(hào)測(cè)量能譜(圖11虛線(xiàn))的137Cs全能峰分辨率為7.65%,而經(jīng)過(guò)DET型成形算法成形后測(cè)量能譜的137Cs全能峰分辨率為8.34%。

        圖11 梯形成形前后脈沖能譜

        兩種梯形成形算法對(duì)能量分辨率都有顯著的提升。在沒(méi)有其他降噪算法處理時(shí),SET型成形算法的提升效果優(yōu)于DET型的提升效果,這是因?yàn)镾ET型成形算法對(duì)高頻白化噪聲有更加優(yōu)秀的抑制和濾波效果,這與頻率響應(yīng)分析結(jié)果相吻合。

        6 結(jié)語(yǔ)

        通過(guò)改進(jìn)算法在FPGA上的實(shí)現(xiàn)方式,提升了算法運(yùn)行時(shí)鐘頻率,這對(duì)一些采樣率要求高的應(yīng)用場(chǎng)合是一利好消息。本文探討的成形方法對(duì)于后續(xù)更加復(fù)雜的數(shù)字化核脈沖信號(hào)處理具有參考和借鑒價(jià)值。

        1 王海濤, 周建斌, 談?chuàng)P寧, 等. 基于FPGA的伽馬射線(xiàn)脈沖數(shù)字高斯成形算法研究[J]. 核電子學(xué)與探測(cè)技術(shù), 2013, 33(11): 1334?1337. WANG Haitao, ZHOU Jianbin, TAN Yangning,. Research on digital gauss shaping algorithm of gamma-ray pulse based on FPGA[J]. Nuclear Electronics & Detection Technology, 2013, 33(11): 1334?1337.

        2 濮亞男, 曾云, 趙豫斌, 等. 閃爍體中子探測(cè)器讀出ASIC-MaPMT_v10性能測(cè)試的設(shè)計(jì)與實(shí)現(xiàn)[J]. 核技術(shù), 2015, 38(2): 020402. DOI: 10.11889/j.0253-3219.2015. hjs.38.020402. PU Yanan, ZENG Yun, ZHAO Yubin,. Design and implementation of testing system of the readout ASIC-MaPMT_v10 for scintillator neutron detector[J]. Nuclear Techniques, 2015, 38(2): 020402. DOI: 10.11889/j.0253-3219.2015.hjs.38.020402.

        3 楊振雷, 王曉輝, 蘇弘, 等. 一種Flash型FPGA單粒子效應(yīng)測(cè)試方法設(shè)計(jì)及驗(yàn)證[J]. 核技術(shù), 2015, 38(2): 020404. DOI: 10.11889/j.0253-3219.2015.hjs.38.020404. YANG Zhenlei, WANG Xiaohui, SU Hong,. Design and verification of test method for the single event effect in flash-based FPGA[J]. Nuclear Techniques, 2015, 38(2): 020404. DOI: 10.11889/j.0253-3219.2015.hjs.38.020404.

        4 竇玉玲, 張國(guó)光, 王國(guó)寶. 用于脈沖γ劑量測(cè)量裝置數(shù)據(jù)采集系統(tǒng)的設(shè)計(jì)與研制[J]. 核技術(shù), 2014, 37(6): 060402. DOI: 10.11889/j.0253-3219.2014.hjs.37.060402. DOU Yuling, ZHANG Guoguang, WANG Guobao. Design and study of data acquisition system for pulsed γ dose measurement device[J]. Nuclear Techniques, 2014, 37(6): 060402. DOI: 10.11889/j.0253-3219.2014.hjs.37. 060402.

        5 朱玉玉, 焦波, 楊民. 數(shù)字多道中梯形成形算法的FPGA實(shí)現(xiàn)[J]. 核電子學(xué)與探測(cè)技術(shù), 2014, 34(1): 78?82. ZHU Yuyu, JIAO Bo, YANG Min. FPGA-based trapezoidal shaping algorithm in the digital multi-channel system[J]. Nuclear Electronics & Detection Technology, 2014, 34(1): 78?82.

        6 陳小猛, 張羽中, 艾憲蕓, 等. 數(shù)字多道脈沖幅度分析器中數(shù)字Dither加入方法[J]. 核技術(shù), 2015, 38(5): 050401. DOI: 10.11889/j.0253-3219.2015.hjs.38.050401. CHEN Xiaomeng, ZHANG Yuzhong, AI Xianyun,. Method of digital Dither adding for digital multi-channel pulse amplitude analyzer[J]. Nuclear Techniques, 2015, 38(5): 050401. DOI: 10.11889/j.0253-3219.2015.hjs.38. 050401.

        7 張懷強(qiáng), 湯彬, 吳和喜. 數(shù)字核譜儀系統(tǒng)中脈沖幅度提取方法研究[J]. 核技術(shù), 2013, 36(5): 050401. DOI: 10.11889/j.0253-3219.2013.hjs.36.050401. ZHANG Huaiqiang, TANG Bin, WU Hexi. Pulse amplitude extraction in digital nuclear spectrometer system[J]. Nuclear Techniques, 2013, 36(5): 050401. DOI: 10.11889/ j.0253-3219.2013.hjs.36.050401.

        8 李偉男, 楊朝文, 周榮. 基于FPGA脈沖幅度分析器的數(shù)字化基線(xiàn)估計(jì)方法[J]. 核技術(shù), 2015, 38(6): 060403. DOI: 10.11889/j.0253-3219.2015.hjs.38.060403. LI Weinan, YANG Chaowen, ZHOU Rong. Baseline estimation method of digital multi-channel pulse height analyzer based on FPGA[J]. Nuclear Techniques, 2015, 38(6): 060403. DOI: 10.11889/j.0253-3219.2015.hjs.38. 060403.

        9 Yang H T. An improvement in differential non-linearity of analog-digital conversion[J]. Nuclear Electronics & Detection Technology, 2001,21(6): 423?426.

        10 Hu C Y, Kang X, Qu G P,. Theoretical derivation, simulation and experimental research on impulse response function of trapezoidal shaping[J]. Atomic Energy Science & Technology, 2014, 48(2): 341?346.

        11 陳亮, 魏義祥, 屈建石. 便攜式γ譜儀中的核素識(shí)別算法[J]. 清華大學(xué)學(xué)報(bào)(自然科學(xué)版), 2009, 49(5): 635?638. CHEN Liang, WEI Yixiang, QU Jianshi. Nuclide identification algorithm for a portable gamma spectrometer[J]. Journal of Tsinghua University (Science and Technology), 2009, 49(5): 635?638.

        12 肖無(wú)云, 魏義祥, 艾憲蕓. 數(shù)字化多道脈沖幅度分析中的梯形成形算法[J]. 清華大學(xué)學(xué)報(bào)(自然科學(xué)版), 2005, 45(6): 810?812. XIAO Wuyun, WEI Yixiang, AI Xianyun. Trapezoidal shaping algorithm for digital multi-channel pulse height analysis[J]. Journal of Tsinghua University (Science and Technology), 2005, 45(6): 810?812.

        13 Glenn K F. Radiation detection and measurement[M]. 4th ed. John Wiley & Sons, 1989:148?157.

        Trapezoidal filter for digital spectrum acquire

        LIU Yinyu WANG Yudong ZHOU Rong YANG Chaowen

        (College of Physical Science and Technology, Key Laboratory of Radiation Physics and Technology, Ministry of Education, Sichuan University, Chengdu 610064, China)

        Background: Since the 1990s, with the rapid development of various digital devices, nuclear signal digital processing has been widely used. Purpose: To enhance energy resolution of digital spectrum system of nuclear radiation measurement, appropriate filter and pulse shaping algorithm is necessary. Trapezoidal filter is widely used for its advantage of easy to use and narrow pulse width,. Methods: The trapezoidal model and two kinds of trapezoidal filter algorithm based on simplified mathematical model of detector’s output pulse signal are established. Results: The energy resolution of137Cs spectrum has improved to 7.65% from 10.78%. Conclusion: Trapezoid shaping is necessary in nuclear signal digital processing.

        Digital spectroscopy, Pulse model, Trapezoidal shaping, FPGA implementation, Energy resolution

        LIU Yinyu, male, born in1992,graduated from Sichuan University in 2014, master student, focusing on nuclear signal acquisition and processing

        ZHOU Rong, E-mail: zhourong@scu.edu.cn

        2016-11-14, accepted date: 2016-12-06

        TL82

        10.11889/j.0253-3219.2017.hjs.40.020402

        劉寅宇,男,1992年出生,2014年畢業(yè)于四川大學(xué),現(xiàn)為碩士研究生,研究領(lǐng)域?yàn)楹诵盘?hào)獲取與處理

        周榮,E-mail:zhourong@scu.edu.cn

        2016-11-14,

        2016-12-06

        Supported by National Natural Science Foundation of China (No.11205108, No.11475121), Excellent Youth Fund of Sichuan University (No.2016SCU04A13)

        國(guó)家自然科學(xué)基金(No.11205108、No.11475121)、四川大學(xué)優(yōu)秀青年基金(No.2016SCU04A13)資助

        猜你喜歡
        虛線(xiàn)能譜梯形
        玩轉(zhuǎn)梯形
        幼兒100(2023年37期)2023-10-23 11:38:58
        能譜CT在術(shù)前預(yù)測(cè)胰腺癌淋巴結(jié)轉(zhuǎn)移的價(jià)值
        梯形達(dá)人
        一類(lèi)變延遲中立型微分方程梯形方法的漸近估計(jì)
        大牛
        梯形
        啟蒙(3-7歲)(2017年6期)2017-11-27 09:34:55
        M87的多波段輻射過(guò)程及其能譜擬合
        電子材料分析中的能譜干擾峰
        能譜CT和MRI小腸造影的護(hù)理配合
        亚洲人成网站免费播放| 亚洲美女又黄又爽在线观看| 中文字幕精品一二三四五六七八 | 日本久久久久| 国产精品毛片av一区二区三区| 永久免费av无码网站性色av| 亚洲色在线视频| 中日韩字幕中文字幕一区| 精品人妻一区二区三区不卡毛片| 中文字幕影片免费人妻少妇| 蜜桃成熟时在线观看免费视频| 变态另类人妖一区二区三区| 国产乱人偷精品人妻a片| 欧美bbw极品另类| 国产成人精品午夜福利免费APP | 99久久免费看少妇高潮a片特黄| 精品国产18禁久久久久久久| 漂亮的小少妇诱惑内射系列| 亚洲av乱码国产精品观看麻豆| 成熟妇女毛茸茸性视频| 国产猛男猛女超爽免费视频| 亚洲 另类 日韩 制服 无码| 无码一区二区三区老色鬼| 久久久久亚洲AV成人网毛片| 国产高清丝袜美腿视频在线观看| 亚洲精品456在线播放狼人 | 亚洲成av人综合在线观看| 国产欧美日韩一区二区三区| 亚洲人成绝费网站色www| 国产成人综合色在线观看网站| 亚洲七七久久综合桃花| 欧美亚洲日韩国产人成在线播放 | 国产成人综合亚洲国产| 精品女厕偷拍视频一区二区| 欧美老肥婆牲交videos| 无码一区二区三区在线| 欧美视频在线观看一区二区| 久久精品国产av大片| 亚洲午夜精品第一区二区| 激情五月婷婷一区二区| 久热国产vs视频在线观看|