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

    基于多波長數(shù)字鎖相檢測的光學拓撲成像方法

    2013-12-06 12:10:58朱蘋蘋張麗敏
    關鍵詞:鎖相異質光學

    高 峰 ,朱蘋蘋,易 茜,張 偉,陳 琛,張麗敏

    (1.天津大學精密儀器與光電子工程學院,天津300072;2. 天津市生物醫(yī)學檢測技術與儀器重點實驗室,天津300072)

    光學拓撲成像是指在組織體表面利用反射測量給出組織體表面下淺層的光學參數(shù)二維變化.光學拓撲成像主要應用于腦功能成像研究,即監(jiān)測生理或心理活動引起的大腦皮層內各調控區(qū)的應激反應過程[1],其作為一種近紅外光譜測量方法,具有完全無創(chuàng)性,高時間分辨率,以及可在自然條件下長時間測量等諸多優(yōu)點.目前它已發(fā)展成為臨床診斷和監(jiān)護中的有效手段,成為現(xiàn)有腦功能成像模態(tài)(功能磁共振成像、腦電圖等)的一種極為有益的補充.

    利用反射光強信息,光學拓撲成像可以重建出組織體內部重要生理參數(shù)的變化分布圖[2],例如氧合血紅蛋白和脫氧血紅蛋白的濃度變化分布.雖然此成像方式不能測量確切的血容量濃度,但是可以測量分子物質的濃度變化[3],可實時連續(xù)無創(chuàng)檢測血氧的相對濃度變化,提供腦組織中血氧的動態(tài)變化情況,已成為臨床診斷和監(jiān)護的重要手段之一.

    考慮到光學拓撲成像方法的優(yōu)點以及應用前景,多家公司已經研制出光學拓撲成像系統(tǒng),但其對探測器要求很高,成本昂貴難以普及.筆者提出的基于數(shù)字鎖相檢測技術的光學拓撲成像系統(tǒng),采用光開關簡化了系統(tǒng)對光源的編碼過程,并且利用數(shù)字鎖相檢測技術提高了系統(tǒng)信噪比,實現(xiàn)了多波長同步動態(tài)測量,降低了系統(tǒng)對探測器的要求.

    光學拓撲成像可以選擇時間分辨測量、頻域測量和連續(xù)波測量3 種模式.經比較,連續(xù)波測量有更實用的優(yōu)點:高信噪比,高時間分辨率,動態(tài)范圍大,簡單性,成本低,臨床方便等[2,4].故本文選擇連續(xù)波測量模式實現(xiàn)光學拓撲成像.假定每個測量區(qū)域的散射系數(shù)是均勻的,且在相鄰的2 次生理活動前后不發(fā)生變化[5],在此基礎上以血紅蛋白的吸收來測量組織中的氧合情況.

    光入射到組織體后,其強度會以指數(shù)形式衰減,而數(shù)字鎖相檢測技術就是一種可以實現(xiàn)微弱光信號動態(tài)檢測的有效技術.相對于其他數(shù)字檢測技術,數(shù)字鎖相檢測技術可以實現(xiàn)多路調制信號的并行檢測,在相同的采樣時間內有更好的抗噪性能[6-7].增大源-探測光纖間距,使得多波長耦合信號穿透表面組織到達目標位置時,依然可以利用此檢測技術于噪聲中提取到有用信息.

    本文詳細介紹了每個組成部分,設計出符合要求的源-探測光纖排布以及測量支架,并且模擬組織體光學參數(shù)進行了液態(tài)仿體的實驗研究.光學拓撲成像實驗過程中,系統(tǒng)將兩路不同反射波長的激光信號耦合后作為入射光源,并以數(shù)字鎖相技術同步檢測反射光強.以修正的Beer-Lamber 定律為模型計算出每對鄰近的源-探測光纖處生理或心理活動前后血紅蛋白的濃度變化量,利用此信息大致反映出大腦皮層內部的氧合狀況.

    1 光學拓撲成像系統(tǒng)

    如圖1 所示,基于數(shù)字鎖相檢測技術的光學拓撲成像系統(tǒng)主要由光源系統(tǒng)、組織體與測量支架、檢測處理系統(tǒng)組成.此系統(tǒng)可實現(xiàn)多路調制信號入射組織體,同步信號監(jiān)測、數(shù)據(jù)采集以及并行數(shù)字鎖相處理等功能.光源系統(tǒng)由單片機控制的兩路調制信號源、激光器控制模塊、激光器輸出尾纖、波分復用器以及光開關組成,可以提供多路入射光源.檢測處理系統(tǒng)由光電探測器、數(shù)據(jù)采集卡以及PC 機組成,可以并行檢測、采集組織體表面四路反射光強,并將采集的信息輸入到PC 機中進行鎖相處理、數(shù)值計算.根據(jù)修正后Beer-Lamber 定律,可以由處理后獲得的幅值信息計算出采樣點處的 aμΔ ,進行線性插值后即可得到重建的圖像.

    圖1 系統(tǒng)組成Fig.1 System configuration

    為更好地實現(xiàn)本文設計的光學拓撲成像系統(tǒng),選擇了符合設計需求的元器件和設備,得到該系統(tǒng)的實物連接圖,見圖1(b).下面詳細介紹具體的選擇情況.

    1.1 光源系統(tǒng)

    為了更好地區(qū)分氧合血紅蛋白和脫氧血紅蛋白[8],本文采用波長分別為660,nm 和830,nm 的溫控半導體激光器系統(tǒng)作為輻射光源,其由溫度-電流控制模塊(LTC100-B,美國Thorlabs Inc)和激光尾纖(分別為LPS-PM660-FC 以及LPS-PM830-FC,美國Thorlabs Inc)兩部分組成.此外為了實驗過程的簡易可行,采用波分復用器(WDMs,加拿大OZ Optics Ltd)耦合兩路光信號,同時采用光開關實現(xiàn)輻射光源入射位置的切換.

    為提取反射出的微弱有用信號,需要以一定頻率調制激光器.如此數(shù)字鎖相檢測技術才能更好地去除背景噪聲,分辨出多波長耦合信號中的調制信號[7].本文采用單片機接收LABVIEW 發(fā)送的信息,控制程控放大器(LT6912,美國Linear Technologies Inc)產生正弦調制信號的特定幅度,同時控制直接數(shù)字頻率合成器(direct digital synthesizer,DDS,AD9831,美國Analog Device Inc)產生正弦調制信號的特定頻率.以此為基礎,經過設計、制板與調試,得到多路可調頻可調幅的調制信號源[9].

    1.2 組織體與測量支架

    采用固態(tài)仿體和液態(tài)仿體兩類進行實驗.系統(tǒng)可行性驗證階段,選用的是固態(tài)標準仿體(吸收系數(shù)μa=0.01mm-1,約化散射系數(shù)μs′=1mm-1).在光學拓撲成像驗證階段,暫不能進行人體實驗,故選用的是單獨配置的液態(tài)仿體.

    仿體實驗時選用不同的測量支架搭配進行.以固態(tài)仿體進行實驗時,采用的是平移支架,具體實現(xiàn)過程在后續(xù)實驗中介紹.液態(tài)仿體是目前驗證階段的主體部分,在此詳細介紹液態(tài)仿體的配置以及實驗中用到的測量支架.

    根據(jù)朗伯定律

    式中:I0為初始光強;I 為測量得到的光強;μa為吸收系數(shù);e 為常數(shù)2.718;ε(λ)為摩爾吸光系數(shù),不隨濃度C 和光程的改變而改變.由式(1)和式(2)可知,直接計算出 Δμa,即可反映出血紅蛋白濃度的相對濃度變化信息.

    為了相對精確地測量 aμΔ ,實驗過程中使用液態(tài)仿體實現(xiàn)吸收系數(shù)的變化.將Intralipid 溶液與印度墨水以相應比例混合后即可配置出所需光學參數(shù)的液態(tài)仿體[10],其中混合比例需根據(jù)所需光學參數(shù)來改變.同時,為了改善拓撲成像效果、提高空間分辨率,還要考慮源-探測光纖的分布情況以及異質體的放置位置[11].圖 2(a)展示了設定的光纖排布情況.它滿足以下規(guī)則:①每對源入射光纖和探測光纖之間的距離為20,mm,此距離直接影響光程;②每對源入射光纖和探測光纖的中間點定義為采樣點,且探測光纖測得的信息反映的即是中間點處的 aμΔ ;③每個采樣點和相鄰位置處的采樣點等距.

    在40,mm×40,mm 的源-探測光纖分布區(qū)域,共有 12 個采樣點,相鄰位置的采樣點間隔為14.14,mm.為獲得相對高的空間分辨率,將異質體放置在采樣點位置處.圖2(b)為筆者設計的測量支架,能夠放置源光纖、探測光纖和塑料試管,其中塑料試管直徑10,mm,用來放入異質體.小孔和大孔分別為光纖和異質體放置位置.

    圖2 光纖排布Fig.2 Optode array in the phantom experiment

    1.3 檢測處理系統(tǒng)

    信號檢測部分采用PIN Si光電二極管(PDA36A,美國Thorlabs Inc)實現(xiàn)微弱光信號探測[12],數(shù)據(jù)采集卡(PCI-1712L,臺灣研華Advantech Ltd)實現(xiàn)多通道同步數(shù)據(jù)采集.考慮到設定的源-探測光纖間距以及光電探測器自身靈敏度的限制,經探測器放大之后的信號依然比較微弱時,可以對此時的信號進行簡單的濾波放大處理,以期提高數(shù)據(jù)采集的有效性.

    采集到PC 機內的離散信號為

    式中:A0為直流分量;Ak為2 個頻率 fk(k = 1,2)的諧波幅度;fs為采樣頻率;fk為調制頻率.它們之間滿足

    式中 Ns為數(shù)據(jù)采集卡設定的采樣點數(shù).本文采用的調制信號為10,kHz、15,kHz,采樣頻率為400,kHz,采樣點 Ns為2 000.

    本文需要對采集到的數(shù)據(jù)進行多波長數(shù)字鎖相檢測處理,此過程中需要兩路離散的參考信號:同相參考信號和正交參考信號,此兩路信號要有90°的相位差,其公式為

    將式(3)中的 ()M n 與式(5)和式(6)中的兩路參考信號 ()C n 、 ()S n 分別相乘,并通過低通濾波器濾除交流分量,此過程可以通過離散傅里葉變換的(discrete Fourier transform,DFT)矩陣來實現(xiàn),即

    于是,可得被測信號的幅值信號為

    根據(jù)修正后的Beer-Lamber 定律,被測信號的幅值信息與血紅蛋白濃度的關系為式中:Coxy和 Cdeoxy分別為含氧血紅蛋白與脫氧血紅蛋白的濃度;Aio(λi)為波長λi的入射光強;Ai(λi)為靜止狀態(tài)時波長λi的輸出光強;L 為設定的源-探測器光纖間距;B 為光程因子,本文引用其經驗值;ε為血紅蛋白的摩爾吸光系數(shù),在特定波長處為常數(shù)[13].

    腦活動狀態(tài)時,被測信號的幅值信息與血紅蛋白濃度的關系為

    式中 Ai′ (λi)為血紅蛋白濃度變化后(活動狀態(tài))波長λi的輸出光強.因此,血紅蛋白的相對濃度變化量如式(11)所示:

    由式(11)即可計算出血紅蛋白濃度的相對變化量.目前系統(tǒng)暫時處于模擬驗證階段,本文以吸收系數(shù)的變化來反映血紅蛋白濃度的變化情況,故后續(xù)實驗過程中,探測的是吸收系數(shù)的相對變化量.

    2 實 驗

    2.1 光學拓撲成像方法的可行性驗證

    實驗過程中簡化系統(tǒng),以驗證光學拓撲成像方法的可行性,即光源系統(tǒng)不采用光開關,直接將波分復用器耦合后的光入射仿體,同時檢測處理系統(tǒng)只采用一個光電探測器以及一個采樣通道.本文采用固態(tài)標準仿體模擬組織體,其光學參數(shù):μa=0.01,mm-1,μs′=1,mm-1(800,nm 波長下測得),并采用2 個固定支架固定源光纖和探測光纖.

    實驗分為以下兩類.

    (1) 以10,kHz 正弦信號調制發(fā)射波長為830,nm的激光器,15,kHz 正弦信號調制發(fā)射波長為660,nm的激光器.固定源光纖的支架,以0.5,mm 為步長移動探測光纖支架,實現(xiàn)不同的源-探測光纖間距下幅值信息的檢測.將經過數(shù)字鎖相處理并歸一化后的信號,與外推條件下距離源光纖一定水平距離處反射光強歸一化后的解析解相比較,以此來判斷系統(tǒng)所測信號幅值的有效性.

    圖3中測量值與解析解并不完全一致.這是因為不同波長下,光學參數(shù)不同,基于半無限平面模型所計算出的反射光強不同,因此正常情況下測量的幅值與解析解會有一定的差異.由以上結果可見,入射光經過仿體后,系統(tǒng)測得的幅值隨著源-探測光纖間距的改變而改變,其趨勢與理論上解析解的變化趨勢大體一致.

    (2) 理論上,輸入的調制信號幅值與多波長數(shù)字鎖相檢測處理后得到的信號幅值成正比關系.為相對準確地判斷系統(tǒng)的測量結果,對于兩路調制信號,固定其中一路調制信號的幅值,同時有規(guī)律地改變另一路調制信號的幅值,以此實現(xiàn)較準確的驗證目的.

    圖3 實驗測量值與解析解的對比Fig.3 Comparison between the phantom experimental result and the analytical solution

    圖4所示分別為維持15,kHz 正弦調制信號幅值不變,將10,kHz 正弦調制信號幅值從600,mV 以50,mV 為步長依次降低至100,mV 時的測量幅值曲線,以及維持10,kHz 正弦調制信號幅值不變,將15,kHz 正弦調制信號幅值從600,mV 以50,mV 為步長降低至100,mV 時的測量幅值曲線.

    圖4 激發(fā)光波長和調制頻率分別為830,nm、10,kHz 和660,nm、15,kHz 時測量幅值與調制信號峰峰值的關系曲線 Fig.4 Normalized amplitudes measured with the input single changing using the excitation light of different wavelengths and modulation frequencies:830,nm,10,kHz;660,nm,15,kHz

    由以上結果可知,調制信號的輸入幅值與數(shù)字鎖相檢測處理后得到的幅值成正比關系.測量中可忽略手動調節(jié)的實驗誤差和相同調制信號峰峰值時不同波長帶來的吸收系數(shù)差異.綜合以上兩類實驗結果可知,光學拓撲成像方法是可行的,系統(tǒng)測量的幅值信息是有效的.

    2.2 光學拓撲成像方法的實驗驗證

    2.2.1 單異質體

    背景液態(tài)仿體的光學參數(shù)設定為:aμ=0.005 mm-1,λ2=660,nm,f2=15,kHz,s1μ′ = mm-1,異質體的參數(shù)為:μa=0.03mm-1,μs′ = 1mm-1(以上光學參數(shù)以660,nm 波長處測得的吸光度來配置),重建圖像見圖5.圖5(a)、(b)分別為830,nm、660,nm 波長處Δμa的重建圖像,異質體位置處求得的 Δμa分別為:Δμa(830) =0.003 9,mm-1和 Δμa(660)=0.005,9,mm-1.

    在660,nm 波長處,理論吸收系數(shù)的相對變化量為: Δμ(a6 60) = 0.025mm-1,測 量 結 果 的 誤 差 為76.4%.造成該誤差的原因分析如下:首先,在實驗過程中不可避免地會受到外界雜散光的影響,且系統(tǒng)檢測到的反射光強度十分微弱,造成了較低的信噪比;其次,配置的液態(tài)仿體的吸收系數(shù)和散射系數(shù)分別由光譜儀測量和經驗值計算所得,且液態(tài)仿體因重力影響導致縱向密度分布不均勻;此外,實驗過程中還存在人為因素造成的誤差.綜上所述,定性的結果誤差很大,但是由成像圖可知,系統(tǒng)可以有效地分辨異質體且定位較準確,異質體對背景其他位置的影響不大,可用來定性分析、分辨吸收系數(shù)的變化.

    圖5 激發(fā)光波長分別為830 和660,nm 時單異質體的Δμa重建圖像 Fig.5 Reconstructed-images of the single inhomogeneous target using the excitation light of different wavelengths(830,nm,660,nm)

    2.2.2 不同吸收系數(shù)的雙異質體

    背景液態(tài)仿體的光學參數(shù)設定為:μa=0 .02mm-1,約化散射系數(shù)μs′ = 1mm-1,異質體光學參數(shù)為:μa= 0.04mm-1和μa= 0.06mm-1,μs′ = 1mm-1(以上光學參數(shù)是以波長660,nm 處的吸光度來配置),重建圖像如圖6 所示.圖6(a)為830,nm 波長處 Δμa的重建圖像,上、下側異質體位置處測得的 Δμa分別為:Δμa′(8 30) =0.006,1,mm-1和 Δμa(8 30) =0.003,1,mm-1. 圖6(b)為660,nm 波長處 Δμa的重建圖像,上、下側異質體位置處測得的 Δμa分別為:Δμa′(6 60)= 0.007,2,mm-1和 Δμa(6 60)= 0.003,8,mm-1.

    圖6 激發(fā)光波長分別為830 nm 和660 nm 時不同吸收系數(shù)的雙異質體的Δμa 重建圖像 Fig.6 Reconstructed Δμa images of two inhomogeneous tar- gets with different concentrations using the excitation light of different wavelengths(830,nm,660,nm)

    在660,nm 波長處,理論上吸收系數(shù)的相對變化量分別為:Δμa(6 60)= 0.02,mm-1、Δμa′ (660) = 0.04mm-1,故誤差分別為81%、82%,雙異質體的吸收系數(shù)比例的測量誤差為5.25%.測量誤差比較大,其原因與單異質體的情況相同,可以通過優(yōu)化液態(tài)仿體光學參數(shù)的配置以及改善遮光條件等方法降低測量誤差.另外由吸收系數(shù)變化的成像圖可知,系統(tǒng)可以有效地分辨異質體且定位較準確,雙異質體的吸收系數(shù)比例與理論值基本一致,即系統(tǒng)可以對吸收系數(shù)的變化趨勢正確成像.

    [5]Gao F,Zhao H J,Tanikawa Y,et al. Optical tomographic mapping of cerebral haemodynamics by means of time-domain detection:Methodology and phantom validation[J]. Physics in Medicine and Biology,2004,49:1055-1078.

    [6]Lasker J M,Masciotti J M,Li Y,et al. Dynamic optical tomographic imager with optimized digital lock-in filtering[J]. SPIE-OSA Biomedical Optics , 2007 ,6629(3):1-11.

    [7]Masciotti J M,Lasker J M,Hielscher A H. Digital lock-in detection for discriminating multiple modulation frequencies with high accuracy and computational efficiency[J]. IEEE Transactions on Instrumentation and Measurement,2008,57(1):182-189.

    [8]Everdell N L,Gibson A P,Tullis I D C. A frequency multiplexed near-infrared topography system of imaging functional activation in the brain[J]. Review of Scientific Instruments,2005,76(9):093705-1-093705-5.

    [9]焦玉婷. 動態(tài)擴散光學層析成像中多波長數(shù)字鎖相光電檢測技術的研究[D]. 天津:天津大學精密儀器與光電子工程學院,2010.Jiao Yuting. Research on Multi-Wavelength Digital Lock-In Photoelectric Detection in Dynamic Diffuse Optical Tomography[D]. Tianjin:School of Precision Instrument and Opto-Electronics Engineering,Tianjin University,2010(in Chinese).

    [10]van Staveren H J,Moes C J M,van Marie,et al. Light scattering in Intralipid-10% in the wavelength range of 400-1100,nm[J]. Applied Optics , 1991 , 30(31):4507-4514.

    [11]Yamamoto T,Maki A,Kadoya T,et al. Arranging optical fibres for the spatial resolution improvement of topographical images[J]. Physics in Medicine and Biology,2002,47:3429-3440.

    [12]安毓英,曾曉東. 光電探測原理[M]. 西安:西安電子科技大學出版社,2009.An Yuying,Zeng Xiaodong. Photoelectrical Detection Theory[M]. Xi’an:Xi’an University of Electronic Science and Technology Press,2009(in Chinese).

    [13]Matcher S J,Cooper C E,Cope M,et al. Performance comparison of several published tissue near-infrared spectroscopy algorithms[J]. Analytical Biochemistry,1995,227:54-68.

    猜你喜歡
    鎖相異質光學
    滑輪組的裝配
    光學常見考題逐個擊破
    基于雙dq鎖相技術的三電平SVG研究
    電測與儀表(2016年6期)2016-04-11 12:06:08
    隨機與異質網(wǎng)絡共存的SIS傳染病模型的定性分析
    Ag2CO3/Ag2O異質p-n結光催化劑的制備及其可見光光催化性能
    寬頻鎖相的一種實現(xiàn)方法
    電測與儀表(2015年8期)2015-04-09 11:50:10
    MoS2/ZnO異質結的光電特性
    物理實驗(2015年10期)2015-02-28 17:36:52
    基于LabVIEW的鎖相放大器的設計與測量
    物理實驗(2015年9期)2015-02-28 17:36:50
    光學遙感壓縮成像技術
    CDSC在單相電網(wǎng)同步鎖相技術中的應用
    国产精品一区二区在线不卡| 激情视频va一区二区三区| 一级作爱视频免费观看| av有码第一页| 国产精品久久久av美女十八| 午夜激情av网站| 亚洲欧洲精品一区二区精品久久久| 亚洲七黄色美女视频| 美女大奶头视频| av免费在线观看网站| 亚洲在线自拍视频| 精品熟女少妇八av免费久了| 国产精品九九99| 日韩欧美国产一区二区入口| 国语自产精品视频在线第100页| 国产精品一区二区在线不卡| 久久性视频一级片| videosex国产| 丁香欧美五月| 亚洲欧美日韩另类电影网站| www.熟女人妻精品国产| 老鸭窝网址在线观看| 国产蜜桃级精品一区二区三区| 亚洲色图av天堂| 正在播放国产对白刺激| 亚洲人成电影观看| 国产精品一区二区三区四区久久 | 免费不卡黄色视频| 亚洲人成77777在线视频| 成人三级黄色视频| 一级黄色大片毛片| 在线播放国产精品三级| 在线观看免费日韩欧美大片| 亚洲精品美女久久av网站| 久久人人97超碰香蕉20202| 欧美绝顶高潮抽搐喷水| 精品国内亚洲2022精品成人| 欧美日韩一级在线毛片| 国产精品久久久久久人妻精品电影| 欧美激情 高清一区二区三区| 每晚都被弄得嗷嗷叫到高潮| 亚洲av片天天在线观看| 色综合欧美亚洲国产小说| 午夜免费激情av| 亚洲av成人av| 91成人精品电影| 久热爱精品视频在线9| 人人澡人人妻人| 欧美成狂野欧美在线观看| 99re在线观看精品视频| 国产亚洲精品久久久久久毛片| 免费在线观看日本一区| 欧美人与性动交α欧美精品济南到| 国产成人精品在线电影| 色在线成人网| 亚洲自偷自拍图片 自拍| 91麻豆av在线| 亚洲一码二码三码区别大吗| 啦啦啦免费观看视频1| 欧美日韩精品网址| 在线播放国产精品三级| 久久久国产成人免费| 亚洲av日韩精品久久久久久密| 97碰自拍视频| 九色国产91popny在线| 一个人免费在线观看的高清视频| 又黄又爽又免费观看的视频| 国产精品 欧美亚洲| 18禁美女被吸乳视频| 亚洲片人在线观看| 每晚都被弄得嗷嗷叫到高潮| 日韩国内少妇激情av| 黄色丝袜av网址大全| 99riav亚洲国产免费| 亚洲无线在线观看| 一级a爱片免费观看的视频| 国产野战对白在线观看| 国产av一区二区精品久久| 天天一区二区日本电影三级 | 久久久久久久久中文| 一区二区三区精品91| 香蕉国产在线看| 如日韩欧美国产精品一区二区三区| 久久久久国产一级毛片高清牌| 午夜两性在线视频| 欧美另类亚洲清纯唯美| 成年人黄色毛片网站| 男女之事视频高清在线观看| 99精品欧美一区二区三区四区| 波多野结衣av一区二区av| 久久精品91蜜桃| 亚洲va日本ⅴa欧美va伊人久久| 中文字幕久久专区| 女人被狂操c到高潮| 禁无遮挡网站| 久久天躁狠狠躁夜夜2o2o| 国产伦人伦偷精品视频| 久久香蕉国产精品| 丁香欧美五月| 日韩免费av在线播放| 女同久久另类99精品国产91| 老司机深夜福利视频在线观看| 午夜久久久久精精品| 久久精品国产99精品国产亚洲性色 | 老司机在亚洲福利影院| 亚洲第一电影网av| 久久久久久亚洲精品国产蜜桃av| 精品欧美国产一区二区三| a级毛片在线看网站| 老汉色av国产亚洲站长工具| 日韩中文字幕欧美一区二区| 久久精品国产综合久久久| 日韩有码中文字幕| 国产精品精品国产色婷婷| 成人特级黄色片久久久久久久| 亚洲午夜理论影院| 亚洲在线自拍视频| 精品午夜福利视频在线观看一区| 国产色视频综合| 免费人成视频x8x8入口观看| 长腿黑丝高跟| 久热这里只有精品99| 神马国产精品三级电影在线观看 | 久久久水蜜桃国产精品网| 在线观看免费视频网站a站| 成人三级黄色视频| 精品不卡国产一区二区三区| 色尼玛亚洲综合影院| 成人国语在线视频| 亚洲精品久久国产高清桃花| 91麻豆精品激情在线观看国产| 男人舔女人的私密视频| 99国产精品99久久久久| 精品午夜福利视频在线观看一区| 国产成年人精品一区二区| 久久欧美精品欧美久久欧美| av免费在线观看网站| 国内久久婷婷六月综合欲色啪| 后天国语完整版免费观看| 色综合站精品国产| 一级a爱片免费观看的视频| 手机成人av网站| 日本三级黄在线观看| 大香蕉久久成人网| 男女床上黄色一级片免费看| 女人高潮潮喷娇喘18禁视频| 熟妇人妻久久中文字幕3abv| 9191精品国产免费久久| 精品第一国产精品| 国产成人av激情在线播放| 国产高清videossex| 亚洲国产精品久久男人天堂| 性欧美人与动物交配| 一进一出抽搐动态| 欧美人与性动交α欧美精品济南到| 免费搜索国产男女视频| 亚洲欧美一区二区三区黑人| 国产色视频综合| 一本久久中文字幕| 18禁美女被吸乳视频| or卡值多少钱| 国产午夜精品久久久久久| av视频免费观看在线观看| 啦啦啦免费观看视频1| 美女大奶头视频| www.自偷自拍.com| 黄色 视频免费看| 国产一级毛片七仙女欲春2 | 丰满人妻熟妇乱又伦精品不卡| 精品久久蜜臀av无| 久久久国产成人精品二区| 1024视频免费在线观看| avwww免费| 亚洲中文字幕日韩| 亚洲精品国产一区二区精华液| 久久人妻福利社区极品人妻图片| 免费看a级黄色片| 在线观看舔阴道视频| 国产精品秋霞免费鲁丝片| 国产三级黄色录像| 亚洲中文字幕一区二区三区有码在线看 | 后天国语完整版免费观看| 欧美午夜高清在线| 在线观看www视频免费| 少妇 在线观看| 欧美日韩精品网址| 又紧又爽又黄一区二区| 动漫黄色视频在线观看| 夜夜看夜夜爽夜夜摸| 亚洲五月色婷婷综合| 一二三四在线观看免费中文在| 国语自产精品视频在线第100页| 亚洲 欧美 日韩 在线 免费| 欧美亚洲日本最大视频资源| 欧美日韩黄片免| 久久久久久人人人人人| 午夜福利18| 午夜免费观看网址| 国产一区二区在线av高清观看| 在线天堂中文资源库| 国产午夜精品久久久久久| 亚洲人成77777在线视频| 国产亚洲精品久久久久5区| 一级片免费观看大全| 午夜精品在线福利| 日本撒尿小便嘘嘘汇集6| 亚洲男人的天堂狠狠| 91大片在线观看| 国产99久久九九免费精品| 中文字幕精品免费在线观看视频| 一级a爱视频在线免费观看| 久久久久久免费高清国产稀缺| 亚洲av成人不卡在线观看播放网| 日本欧美视频一区| 午夜福利一区二区在线看| 黄片小视频在线播放| 天堂动漫精品| 99国产精品免费福利视频| 久久香蕉国产精品| 亚洲av熟女| 国产极品粉嫩免费观看在线| 成熟少妇高潮喷水视频| 久久人人爽av亚洲精品天堂| 色综合婷婷激情| 中文字幕精品免费在线观看视频| 久热爱精品视频在线9| 亚洲 国产 在线| 国产又爽黄色视频| 免费高清视频大片| 叶爱在线成人免费视频播放| 色播亚洲综合网| 男人操女人黄网站| 亚洲一区高清亚洲精品| 亚洲在线自拍视频| 在线十欧美十亚洲十日本专区| 亚洲 欧美 日韩 在线 免费| 男女午夜视频在线观看| 一区二区日韩欧美中文字幕| 国产精品九九99| 亚洲精华国产精华精| 国产熟女午夜一区二区三区| 亚洲欧美精品综合久久99| 黄色a级毛片大全视频| 搞女人的毛片| 老司机靠b影院| avwww免费| 亚洲人成77777在线视频| 日本 欧美在线| 99精品久久久久人妻精品| 一卡2卡三卡四卡精品乱码亚洲| 丁香六月欧美| 激情在线观看视频在线高清| 精品国产亚洲在线| 国产在线观看jvid| 99热只有精品国产| 午夜精品在线福利| 最新美女视频免费是黄的| 香蕉久久夜色| 亚洲国产欧美日韩在线播放| 日韩成人在线观看一区二区三区| 精品乱码久久久久久99久播| 久久人妻熟女aⅴ| 黄片大片在线免费观看| 精品欧美一区二区三区在线| 脱女人内裤的视频| 九色亚洲精品在线播放| 18禁观看日本| 在线永久观看黄色视频| 亚洲人成77777在线视频| 国产色视频综合| 国产成人精品久久二区二区免费| 欧美日韩福利视频一区二区| 国产成人欧美在线观看| 啦啦啦观看免费观看视频高清 | 久久久久久国产a免费观看| 1024视频免费在线观看| 成人国语在线视频| 亚洲精品av麻豆狂野| 激情在线观看视频在线高清| 成人国产一区最新在线观看| 一进一出好大好爽视频| 91av网站免费观看| 欧美人与性动交α欧美精品济南到| 成人亚洲精品av一区二区| 丝袜美足系列| 国产视频一区二区在线看| 人人妻人人澡欧美一区二区 | 欧美黑人欧美精品刺激| 我的亚洲天堂| 日本 av在线| 成人亚洲精品一区在线观看| 国产欧美日韩一区二区三区在线| www.自偷自拍.com| 午夜福利一区二区在线看| 免费高清在线观看日韩| 久久久久九九精品影院| 美女扒开内裤让男人捅视频| 成人欧美大片| 久久香蕉精品热| 午夜福利18| 欧美激情极品国产一区二区三区| 一级片免费观看大全| 国产麻豆成人av免费视频| 黑人欧美特级aaaaaa片| 九色亚洲精品在线播放| 亚洲美女黄片视频| 国产欧美日韩精品亚洲av| 18美女黄网站色大片免费观看| 首页视频小说图片口味搜索| 日日摸夜夜添夜夜添小说| 中文亚洲av片在线观看爽| 亚洲国产看品久久| 黄色视频,在线免费观看| 国产亚洲av嫩草精品影院| 亚洲一区高清亚洲精品| 精品久久蜜臀av无| 黄片播放在线免费| 天天一区二区日本电影三级 | 日日干狠狠操夜夜爽| x7x7x7水蜜桃| 欧美久久黑人一区二区| 亚洲一区中文字幕在线| 男女午夜视频在线观看| 90打野战视频偷拍视频| 日本在线视频免费播放| netflix在线观看网站| 国产成人欧美| 免费在线观看完整版高清| 亚洲欧美一区二区三区黑人| 性少妇av在线| 久久国产亚洲av麻豆专区| 99久久久亚洲精品蜜臀av| 宅男免费午夜| 啪啪无遮挡十八禁网站| 亚洲免费av在线视频| 日韩欧美国产在线观看| 在线国产一区二区在线| 伊人久久大香线蕉亚洲五| 精品乱码久久久久久99久播| 国产高清有码在线观看视频 | 久久久精品国产亚洲av高清涩受| 制服人妻中文乱码| 一边摸一边做爽爽视频免费| 国产不卡一卡二| 精品国产乱码久久久久久男人| 性少妇av在线| 国产亚洲精品av在线| 欧美一区二区精品小视频在线| 精品国产美女av久久久久小说| 18美女黄网站色大片免费观看| 久久久国产成人精品二区| 好看av亚洲va欧美ⅴa在| www.自偷自拍.com| 久久香蕉国产精品| tocl精华| 成人亚洲精品av一区二区| 午夜日韩欧美国产| 国产精品久久久久久亚洲av鲁大| 18禁美女被吸乳视频| 色综合欧美亚洲国产小说| 精品高清国产在线一区| 侵犯人妻中文字幕一二三四区| 韩国精品一区二区三区| 在线免费观看的www视频| 黄片大片在线免费观看| 黑人操中国人逼视频| 巨乳人妻的诱惑在线观看| 亚洲狠狠婷婷综合久久图片| 日本一区二区免费在线视频| 午夜成年电影在线免费观看| 少妇熟女aⅴ在线视频| ponron亚洲| 国产精品久久久久久精品电影 | 亚洲专区国产一区二区| 精品人妻在线不人妻| 久久人妻福利社区极品人妻图片| 欧美色欧美亚洲另类二区 | 精品一品国产午夜福利视频| 成在线人永久免费视频| 后天国语完整版免费观看| 成人亚洲精品av一区二区| 国产1区2区3区精品| 夜夜看夜夜爽夜夜摸| 国产精品久久久久久精品电影 | 国产精品香港三级国产av潘金莲| or卡值多少钱| 久久这里只有精品19| 欧美最黄视频在线播放免费| 一区二区三区高清视频在线| 国产一区二区三区综合在线观看| 后天国语完整版免费观看| 国内精品久久久久精免费| 国产高清视频在线播放一区| 视频区欧美日本亚洲| 欧美黄色淫秽网站| 国产精品乱码一区二三区的特点 | 可以在线观看的亚洲视频| 可以在线观看毛片的网站| 亚洲全国av大片| 亚洲性夜色夜夜综合| bbb黄色大片| 欧美乱码精品一区二区三区| 亚洲国产精品sss在线观看| 国产成人精品无人区| 久久伊人香网站| 亚洲精品国产区一区二| 亚洲av成人av| 亚洲欧美激情综合另类| 亚洲一区中文字幕在线| 亚洲av美国av| 日本 欧美在线| 两人在一起打扑克的视频| 少妇熟女aⅴ在线视频| 成人三级黄色视频| 久久人妻熟女aⅴ| 91国产中文字幕| 国产激情欧美一区二区| 女同久久另类99精品国产91| 国产又色又爽无遮挡免费看| 国产高清有码在线观看视频 | 一级毛片高清免费大全| 亚洲精品在线观看二区| av视频在线观看入口| www.精华液| 动漫黄色视频在线观看| 欧美成人性av电影在线观看| 黄片播放在线免费| 午夜精品国产一区二区电影| 久久久国产精品麻豆| 婷婷丁香在线五月| 欧美日韩瑟瑟在线播放| 国产精品,欧美在线| 一区在线观看完整版| 咕卡用的链子| 久久狼人影院| 久久性视频一级片| 亚洲国产日韩欧美精品在线观看 | 一进一出抽搐动态| 丝袜美足系列| 亚洲片人在线观看| 嫁个100分男人电影在线观看| 人人澡人人妻人| 国产精品98久久久久久宅男小说| 后天国语完整版免费观看| 一级毛片女人18水好多| 美女大奶头视频| 久久久久国内视频| 日日夜夜操网爽| 中亚洲国语对白在线视频| 久久九九热精品免费| 好看av亚洲va欧美ⅴa在| 午夜福利一区二区在线看| 国产精品爽爽va在线观看网站 | 99精品欧美一区二区三区四区| 热99re8久久精品国产| 亚洲最大成人中文| 日韩欧美三级三区| 久久久久久亚洲精品国产蜜桃av| 国语自产精品视频在线第100页| 国产精品一区二区三区四区久久 | 亚洲熟妇熟女久久| 狂野欧美激情性xxxx| 国产视频一区二区在线看| 精品久久久久久久毛片微露脸| 丰满人妻熟妇乱又伦精品不卡| 黑人巨大精品欧美一区二区蜜桃| 免费人成视频x8x8入口观看| 99精品在免费线老司机午夜| 国产不卡一卡二| 男人操女人黄网站| 欧美日韩瑟瑟在线播放| 色综合站精品国产| 9色porny在线观看| 日本免费一区二区三区高清不卡 | 国产成人免费无遮挡视频| 人人妻人人澡欧美一区二区 | 不卡av一区二区三区| 88av欧美| 精品日产1卡2卡| 每晚都被弄得嗷嗷叫到高潮| 亚洲欧美一区二区三区黑人| 怎么达到女性高潮| 在线观看免费视频网站a站| 黄色片一级片一级黄色片| 97人妻天天添夜夜摸| 国产一区二区在线av高清观看| 美女扒开内裤让男人捅视频| 欧美一区二区精品小视频在线| 日韩av在线大香蕉| 国产激情欧美一区二区| 国产精品香港三级国产av潘金莲| 欧美色视频一区免费| 欧美精品啪啪一区二区三区| 国产精品九九99| 国产成人啪精品午夜网站| 亚洲全国av大片| 日本 欧美在线| 精品高清国产在线一区| 日韩欧美一区视频在线观看| 极品人妻少妇av视频| 老汉色∧v一级毛片| 亚洲av第一区精品v没综合| 夜夜躁狠狠躁天天躁| 亚洲成人精品中文字幕电影| 亚洲一卡2卡3卡4卡5卡精品中文| 日韩精品中文字幕看吧| 国产精品九九99| 久久久久国产一级毛片高清牌| 老熟妇仑乱视频hdxx| 久久精品成人免费网站| 免费在线观看视频国产中文字幕亚洲| 亚洲成av片中文字幕在线观看| 欧美另类亚洲清纯唯美| 免费看a级黄色片| 精品欧美国产一区二区三| 侵犯人妻中文字幕一二三四区| 黄频高清免费视频| 成人精品一区二区免费| 亚洲,欧美精品.| 99久久综合精品五月天人人| 亚洲 国产 在线| 香蕉国产在线看| 国产精品98久久久久久宅男小说| 极品教师在线免费播放| 久久 成人 亚洲| 男女床上黄色一级片免费看| 国产精品 欧美亚洲| 最好的美女福利视频网| 麻豆国产av国片精品| 在线观看66精品国产| 搡老妇女老女人老熟妇| 国产亚洲精品综合一区在线观看 | 久久午夜综合久久蜜桃| 国产不卡一卡二| 国产亚洲欧美98| 亚洲全国av大片| 亚洲一卡2卡3卡4卡5卡精品中文| 黑人巨大精品欧美一区二区mp4| 亚洲七黄色美女视频| 国产亚洲精品av在线| av有码第一页| 国产成人av激情在线播放| 欧美午夜高清在线| 亚洲第一欧美日韩一区二区三区| 国产精品久久久久久亚洲av鲁大| 久久香蕉激情| 国产午夜福利久久久久久| 1024视频免费在线观看| 女性被躁到高潮视频| 91av网站免费观看| 精品卡一卡二卡四卡免费| 亚洲精品美女久久av网站| 久久热在线av| 91麻豆精品激情在线观看国产| 亚洲欧美精品综合久久99| 男女午夜视频在线观看| 国产成人av教育| 欧美一区二区精品小视频在线| 欧美成狂野欧美在线观看| 亚洲人成网站在线播放欧美日韩| 人成视频在线观看免费观看| 又黄又粗又硬又大视频| 久久中文字幕人妻熟女| 午夜日韩欧美国产| 淫秽高清视频在线观看| 丁香六月欧美| 久久香蕉国产精品| 亚洲 欧美一区二区三区| 中文字幕人妻丝袜一区二区| 欧美大码av| 他把我摸到了高潮在线观看| 午夜久久久在线观看| 国产不卡一卡二| 欧美成人免费av一区二区三区| 婷婷精品国产亚洲av在线| 国产精品美女特级片免费视频播放器 | 大码成人一级视频| 中文字幕高清在线视频| 亚洲专区中文字幕在线| 亚洲av第一区精品v没综合| 搡老熟女国产l中国老女人| 成在线人永久免费视频| 精品卡一卡二卡四卡免费| 大陆偷拍与自拍| 国产熟女xx| 久久中文字幕人妻熟女| 亚洲一区二区三区不卡视频| 免费观看人在逋| 超碰成人久久| 91av网站免费观看| 超碰成人久久| 免费人成视频x8x8入口观看| 国产精品一区二区精品视频观看| 色精品久久人妻99蜜桃| 一边摸一边做爽爽视频免费| 国产色视频综合| 每晚都被弄得嗷嗷叫到高潮| 午夜精品在线福利| 亚洲专区中文字幕在线| 欧美日韩福利视频一区二区| 悠悠久久av| 在线观看午夜福利视频| 国产精品亚洲一级av第二区| 久久性视频一级片| 少妇粗大呻吟视频| 非洲黑人性xxxx精品又粗又长| 亚洲精品一区av在线观看| www国产在线视频色| 一a级毛片在线观看| 日本三级黄在线观看| 国产免费av片在线观看野外av| 在线播放国产精品三级| av天堂久久9| 国产欧美日韩一区二区三| 国产人伦9x9x在线观看|