• <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)同步鎖相技術中的應用
    欧美激情极品国产一区二区三区| 国产欧美日韩一区二区三 | 黄色毛片三级朝国网站| 欧美日韩亚洲国产一区二区在线观看 | 成人三级做爰电影| 菩萨蛮人人尽说江南好唐韦庄| 久9热在线精品视频| 中文字幕人妻熟女乱码| 伊人亚洲综合成人网| 欧美xxⅹ黑人| 99精国产麻豆久久婷婷| 国产熟女欧美一区二区| 精品少妇黑人巨大在线播放| 亚洲男人天堂网一区| 欧美黑人精品巨大| 久热这里只有精品99| 无遮挡黄片免费观看| 亚洲精品一区蜜桃| 亚洲熟女毛片儿| 国产视频一区二区在线看| 妹子高潮喷水视频| 成年人午夜在线观看视频| 精品久久久精品久久久| 咕卡用的链子| 国产免费一区二区三区四区乱码| 久久久久久久精品精品| 国产精品一区二区在线不卡| 精品一区在线观看国产| av又黄又爽大尺度在线免费看| 亚洲伊人久久精品综合| 欧美日韩一级在线毛片| 久久性视频一级片| 成人黄色视频免费在线看| 两性夫妻黄色片| 香蕉国产在线看| 亚洲av成人精品一二三区| 少妇人妻久久综合中文| 日本色播在线视频| 国产一区二区三区av在线| 亚洲一卡2卡3卡4卡5卡精品中文| 久久九九热精品免费| 亚洲一卡2卡3卡4卡5卡精品中文| 99精品久久久久人妻精品| 波多野结衣一区麻豆| 国产精品偷伦视频观看了| 只有这里有精品99| 午夜激情av网站| 亚洲五月色婷婷综合| 99热全是精品| 国产亚洲av高清不卡| 久久久久久免费高清国产稀缺| 亚洲精品一卡2卡三卡4卡5卡 | 91字幕亚洲| 黄色视频在线播放观看不卡| 青春草亚洲视频在线观看| e午夜精品久久久久久久| 日韩av不卡免费在线播放| a级毛片在线看网站| 日韩人妻精品一区2区三区| 一级片免费观看大全| 欧美xxⅹ黑人| 99国产精品一区二区蜜桃av | 亚洲九九香蕉| 黄频高清免费视频| 男女免费视频国产| 亚洲欧美精品自产自拍| 99国产综合亚洲精品| 深夜精品福利| 欧美黑人欧美精品刺激| 久久女婷五月综合色啪小说| 男的添女的下面高潮视频| 又大又爽又粗| 日本一区二区免费在线视频| 又大又黄又爽视频免费| 69精品国产乱码久久久| 美女主播在线视频| 别揉我奶头~嗯~啊~动态视频 | 成年美女黄网站色视频大全免费| 免费高清在线观看日韩| 精品国产一区二区三区久久久樱花| 视频区欧美日本亚洲| 五月开心婷婷网| 看免费成人av毛片| 亚洲久久久国产精品| 国产欧美亚洲国产| 亚洲久久久国产精品| videosex国产| 国产亚洲精品久久久久5区| 午夜福利免费观看在线| 国产熟女欧美一区二区| 亚洲精品自拍成人| 国产熟女午夜一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| 色婷婷av一区二区三区视频| 可以免费在线观看a视频的电影网站| 国产亚洲欧美精品永久| 欧美大码av| 黄色毛片三级朝国网站| 欧美精品啪啪一区二区三区 | 人人澡人人妻人| 国产av国产精品国产| 天天躁日日躁夜夜躁夜夜| 欧美日韩精品网址| 亚洲精品中文字幕在线视频| 午夜免费成人在线视频| 国产在线观看jvid| 国产淫语在线视频| 久久久欧美国产精品| 久久精品人人爽人人爽视色| 国产女主播在线喷水免费视频网站| 久久99精品国语久久久| 欧美人与性动交α欧美精品济南到| 国产一区二区三区av在线| 成人国语在线视频| 免费av中文字幕在线| 看免费av毛片| 欧美 亚洲 国产 日韩一| 黄网站色视频无遮挡免费观看| 亚洲专区中文字幕在线| 无遮挡黄片免费观看| 亚洲,欧美,日韩| 国产男女超爽视频在线观看| 人成视频在线观看免费观看| 最新的欧美精品一区二区| 日本午夜av视频| 啦啦啦在线观看免费高清www| 国产亚洲欧美精品永久| 久久热在线av| 人人妻人人添人人爽欧美一区卜| 在线观看免费高清a一片| 777久久人妻少妇嫩草av网站| 性色av一级| 丁香六月天网| 18禁国产床啪视频网站| 美女脱内裤让男人舔精品视频| 婷婷色麻豆天堂久久| 国产97色在线日韩免费| 伦理电影免费视频| 亚洲国产欧美在线一区| 久久影院123| 久久精品国产综合久久久| 亚洲欧美日韩另类电影网站| 国产在线一区二区三区精| 久久国产精品男人的天堂亚洲| 久热这里只有精品99| 欧美激情极品国产一区二区三区| 精品免费久久久久久久清纯 | 在线精品无人区一区二区三| 男女床上黄色一级片免费看| 成人免费观看视频高清| 亚洲欧美一区二区三区黑人| 国产色视频综合| 国产精品欧美亚洲77777| 91精品伊人久久大香线蕉| 夜夜骑夜夜射夜夜干| 99热全是精品| 精品人妻1区二区| 久久毛片免费看一区二区三区| 成人亚洲精品一区在线观看| 日韩av在线免费看完整版不卡| 69精品国产乱码久久久| av又黄又爽大尺度在线免费看| 18禁裸乳无遮挡动漫免费视频| 国产亚洲欧美精品永久| 99国产精品99久久久久| 一边摸一边抽搐一进一出视频| 丰满饥渴人妻一区二区三| 日韩 欧美 亚洲 中文字幕| 日韩,欧美,国产一区二区三区| 国产亚洲精品久久久久5区| 亚洲一卡2卡3卡4卡5卡精品中文| 国产激情久久老熟女| 亚洲情色 制服丝袜| 亚洲欧美清纯卡通| 精品国产超薄肉色丝袜足j| 午夜影院在线不卡| 国产91精品成人一区二区三区 | 50天的宝宝边吃奶边哭怎么回事| 成年人午夜在线观看视频| 久久人人97超碰香蕉20202| 黄片小视频在线播放| 欧美xxⅹ黑人| 中文字幕人妻熟女乱码| 搡老乐熟女国产| 午夜激情av网站| 国产日韩欧美亚洲二区| 国产又色又爽无遮挡免| 天天添夜夜摸| 交换朋友夫妻互换小说| 如日韩欧美国产精品一区二区三区| 啦啦啦在线观看免费高清www| 久久热在线av| 曰老女人黄片| 久久精品国产a三级三级三级| 蜜桃在线观看..| 大型av网站在线播放| 国产欧美日韩综合在线一区二区| 国产成人啪精品午夜网站| 50天的宝宝边吃奶边哭怎么回事| 国产老妇伦熟女老妇高清| 秋霞在线观看毛片| 国产日韩欧美亚洲二区| 久久久欧美国产精品| 成在线人永久免费视频| 久久精品成人免费网站| 人人妻人人澡人人看| 久久热在线av| 老司机亚洲免费影院| 十八禁高潮呻吟视频| 免费在线观看完整版高清| 免费观看av网站的网址| 国产熟女午夜一区二区三区| 黑人猛操日本美女一级片| 精品久久久精品久久久| 亚洲欧美日韩另类电影网站| 美女午夜性视频免费| 黄色视频不卡| 亚洲精品久久久久久婷婷小说| 一级毛片 在线播放| 黑人欧美特级aaaaaa片| 免费日韩欧美在线观看| 韩国精品一区二区三区| 深夜精品福利| 麻豆av在线久日| 亚洲免费av在线视频| 亚洲,一卡二卡三卡| 久久久久网色| 国产爽快片一区二区三区| 欧美人与善性xxx| av网站免费在线观看视频| 免费观看a级毛片全部| 自拍欧美九色日韩亚洲蝌蚪91| 国产免费一区二区三区四区乱码| 国产亚洲一区二区精品| 国产一区有黄有色的免费视频| 国产一区二区激情短视频 | 天天添夜夜摸| 观看av在线不卡| 爱豆传媒免费全集在线观看| 国产成人免费无遮挡视频| 亚洲精品第二区| 两个人免费观看高清视频| 国产亚洲精品久久久久5区| 亚洲天堂av无毛| 亚洲欧洲日产国产| 国产片特级美女逼逼视频| 少妇精品久久久久久久| av有码第一页| 久久99精品国语久久久| 亚洲欧洲日产国产| 天天躁日日躁夜夜躁夜夜| 亚洲 国产 在线| 99国产精品免费福利视频| 欧美精品人与动牲交sv欧美| 国产成人免费观看mmmm| 国产免费视频播放在线视频| 999精品在线视频| 久久天堂一区二区三区四区| 中文字幕另类日韩欧美亚洲嫩草| 日韩大码丰满熟妇| 免费av中文字幕在线| 手机成人av网站| 一级a爱视频在线免费观看| 午夜影院在线不卡| 99精国产麻豆久久婷婷| 精品国产一区二区久久| 黄色毛片三级朝国网站| 三上悠亚av全集在线观看| 国产精品久久久人人做人人爽| 国产欧美日韩精品亚洲av| 中文欧美无线码| 青春草亚洲视频在线观看| 国产精品99久久99久久久不卡| 亚洲av成人不卡在线观看播放网 | 国产精品久久久人人做人人爽| 满18在线观看网站| 午夜久久久在线观看| 热99国产精品久久久久久7| 免费黄频网站在线观看国产| 亚洲人成77777在线视频| 少妇人妻久久综合中文| 国精品久久久久久国模美| 国产精品久久久久久精品古装| 国产午夜精品一二区理论片| 天天操日日干夜夜撸| 国产成人免费观看mmmm| 丰满少妇做爰视频| 人人妻,人人澡人人爽秒播 | 最新在线观看一区二区三区 | 亚洲精品一区蜜桃| 老司机影院毛片| 婷婷色综合大香蕉| 黄色怎么调成土黄色| 亚洲国产欧美日韩在线播放| 精品少妇黑人巨大在线播放| 一区在线观看完整版| 看免费av毛片| 纯流量卡能插随身wifi吗| 色94色欧美一区二区| 又大又黄又爽视频免费| 国产成人a∨麻豆精品| 精品亚洲乱码少妇综合久久| 人人妻人人澡人人爽人人夜夜| 女警被强在线播放| 亚洲色图综合在线观看| 999精品在线视频| 成年人免费黄色播放视频| 69精品国产乱码久久久| 亚洲 欧美一区二区三区| 国产成人欧美| 啦啦啦 在线观看视频| av福利片在线| 亚洲av片天天在线观看| 亚洲欧美色中文字幕在线| 国产成人精品久久二区二区91| 中文精品一卡2卡3卡4更新| 一边摸一边抽搐一进一出视频| 国产主播在线观看一区二区 | 色婷婷久久久亚洲欧美| 国产真人三级小视频在线观看| 国产免费福利视频在线观看| 大片电影免费在线观看免费| 国产亚洲av高清不卡| 久久亚洲国产成人精品v| 日韩av不卡免费在线播放| 国产视频首页在线观看| 水蜜桃什么品种好| 欧美黑人欧美精品刺激| 日本色播在线视频| 国产成人a∨麻豆精品| 在线观看人妻少妇| 国产伦人伦偷精品视频| 一级毛片 在线播放| 一级毛片我不卡| 日韩视频在线欧美| 欧美在线一区亚洲| www.精华液| 精品国产超薄肉色丝袜足j| 国产欧美日韩精品亚洲av| 国产极品粉嫩免费观看在线| 91精品伊人久久大香线蕉| 国产野战对白在线观看| 欧美xxⅹ黑人| 激情五月婷婷亚洲| 妹子高潮喷水视频| www.av在线官网国产| 国产成人91sexporn| 又黄又粗又硬又大视频| av视频免费观看在线观看| 久久性视频一级片| 久久精品久久精品一区二区三区| 亚洲少妇的诱惑av| 高清黄色对白视频在线免费看| 欧美国产精品一级二级三级| 亚洲国产欧美日韩在线播放| 夫妻性生交免费视频一级片| 国产精品人妻久久久影院| 欧美日韩av久久| 黄片小视频在线播放| 69精品国产乱码久久久| 亚洲精品国产av成人精品| 在线观看免费视频网站a站| 一区二区三区精品91| 真人做人爱边吃奶动态| 亚洲av在线观看美女高潮| 亚洲欧美精品自产自拍| 丰满迷人的少妇在线观看| 成人国产一区最新在线观看 | 亚洲精品av麻豆狂野| 涩涩av久久男人的天堂| 日韩人妻精品一区2区三区| 欧美精品av麻豆av| 欧美成人精品欧美一级黄| 久久人妻福利社区极品人妻图片 | 久久精品熟女亚洲av麻豆精品| 麻豆av在线久日| 黑人欧美特级aaaaaa片| 国产在线一区二区三区精| 伊人亚洲综合成人网| 欧美日韩视频高清一区二区三区二| 亚洲国产精品国产精品| 亚洲黑人精品在线| 看免费av毛片| 极品少妇高潮喷水抽搐| 日韩av在线免费看完整版不卡| 在线观看免费高清a一片| 国产又爽黄色视频| 婷婷色麻豆天堂久久| 精品少妇内射三级| 美女中出高潮动态图| 青草久久国产| 亚洲欧美清纯卡通| 国产又色又爽无遮挡免| 99精品久久久久人妻精品| 老司机亚洲免费影院| 国产亚洲精品第一综合不卡| 母亲3免费完整高清在线观看| 天堂8中文在线网| 男女午夜视频在线观看| 亚洲精品日韩在线中文字幕| 啦啦啦在线免费观看视频4| 无遮挡黄片免费观看| 在线观看免费日韩欧美大片| 久久国产亚洲av麻豆专区| 国产亚洲精品第一综合不卡| 亚洲精品日韩在线中文字幕| 丁香六月欧美| 考比视频在线观看| 少妇人妻 视频| 欧美日韩精品网址| 啦啦啦 在线观看视频| 成人国语在线视频| 亚洲成人免费电影在线观看 | 久久热在线av| av在线app专区| 视频在线观看一区二区三区| 啦啦啦中文免费视频观看日本| 丝袜喷水一区| 欧美亚洲 丝袜 人妻 在线| 中文字幕高清在线视频| 国产精品香港三级国产av潘金莲 | 精品国产超薄肉色丝袜足j| 亚洲国产日韩一区二区| 久9热在线精品视频| 中文精品一卡2卡3卡4更新| 十分钟在线观看高清视频www| 日本vs欧美在线观看视频| 两人在一起打扑克的视频| bbb黄色大片| 欧美日韩福利视频一区二区| 精品欧美一区二区三区在线| 超碰成人久久| 熟女av电影| 午夜老司机福利片| 狂野欧美激情性bbbbbb| 国产精品.久久久| 国产片内射在线| 一个人免费看片子| 久久精品久久久久久噜噜老黄| 美女国产高潮福利片在线看| 岛国毛片在线播放| www.熟女人妻精品国产| 乱人伦中国视频| 国产色视频综合| 老司机亚洲免费影院| 19禁男女啪啪无遮挡网站| 国产成人啪精品午夜网站| 日韩制服丝袜自拍偷拍| 日韩人妻精品一区2区三区| 成人手机av| 国产精品一区二区免费欧美 | 欧美+亚洲+日韩+国产| 欧美精品一区二区大全| 久久国产亚洲av麻豆专区| 考比视频在线观看| 菩萨蛮人人尽说江南好唐韦庄| 亚洲三区欧美一区| 成年美女黄网站色视频大全免费| 久久99热这里只频精品6学生| 国产淫语在线视频| 精品国产乱码久久久久久小说| 一级a爱视频在线免费观看| 国产男女超爽视频在线观看| 国产日韩一区二区三区精品不卡| 亚洲精品一区蜜桃| 观看av在线不卡| 一边摸一边做爽爽视频免费| www.av在线官网国产| 成人影院久久| 青春草视频在线免费观看| 国产精品亚洲av一区麻豆| 免费看av在线观看网站| 欧美在线一区亚洲| 国语对白做爰xxxⅹ性视频网站| 久久毛片免费看一区二区三区| 亚洲av男天堂| 老汉色∧v一级毛片| 日韩av不卡免费在线播放| 99久久99久久久精品蜜桃| 国产精品国产av在线观看| 成人国产av品久久久| 男女之事视频高清在线观看 | 男人舔女人的私密视频| 久久99一区二区三区| 啦啦啦中文免费视频观看日本| 国产精品久久久av美女十八| 日本av免费视频播放| 亚洲精品成人av观看孕妇| 男女之事视频高清在线观看 | 在线观看www视频免费| 免费观看av网站的网址| 校园人妻丝袜中文字幕| 国产精品国产av在线观看| 最新在线观看一区二区三区 | 九草在线视频观看| 久久ye,这里只有精品| 丁香六月天网| 精品国产一区二区久久| 精品少妇一区二区三区视频日本电影| 亚洲av电影在线进入| 高清不卡的av网站| 热99久久久久精品小说推荐| 成年女人毛片免费观看观看9 | 国产免费一区二区三区四区乱码| 人人妻人人爽人人添夜夜欢视频| 精品久久蜜臀av无| 亚洲图色成人| 男女边摸边吃奶| 无限看片的www在线观看| 国产精品九九99| 男女边吃奶边做爰视频| 少妇精品久久久久久久| 日本一区二区免费在线视频| 亚洲五月色婷婷综合| 纵有疾风起免费观看全集完整版| 精品卡一卡二卡四卡免费| 亚洲男人天堂网一区| 多毛熟女@视频| 亚洲成人免费av在线播放| 少妇人妻久久综合中文| 久久久久久人人人人人| 中文欧美无线码| 丝袜脚勾引网站| 日本a在线网址| 国产免费一区二区三区四区乱码| 国产精品.久久久| 一区二区三区四区激情视频| 亚洲第一青青草原| 国产激情久久老熟女| 亚洲国产最新在线播放| 日韩熟女老妇一区二区性免费视频| 欧美大码av| 亚洲精品一二三| 久久久久精品国产欧美久久久 | 老司机影院成人| 色网站视频免费| 一本久久精品| 美女高潮到喷水免费观看| 久久青草综合色| 亚洲精品美女久久久久99蜜臀 | 日韩欧美一区视频在线观看| 成人国语在线视频| 天天躁夜夜躁狠狠躁躁| 一区在线观看完整版| 一区二区三区激情视频| 国产黄色免费在线视频| 侵犯人妻中文字幕一二三四区| 国产成人精品久久久久久| 91九色精品人成在线观看| 99国产综合亚洲精品| av国产精品久久久久影院| 十八禁网站网址无遮挡| av在线播放精品| 十八禁网站网址无遮挡| 青青草视频在线视频观看| 日韩精品免费视频一区二区三区| 亚洲国产欧美在线一区| 国产在线观看jvid| 成人国产av品久久久| 午夜福利在线免费观看网站| 久久ye,这里只有精品| 日本91视频免费播放| e午夜精品久久久久久久| 伊人久久大香线蕉亚洲五| 日本一区二区免费在线视频| 在线亚洲精品国产二区图片欧美| 制服诱惑二区| 女人被躁到高潮嗷嗷叫费观| 别揉我奶头~嗯~啊~动态视频 | 女人被躁到高潮嗷嗷叫费观| 日本色播在线视频| 成年av动漫网址| 国产成人一区二区三区免费视频网站 | 2021少妇久久久久久久久久久| 九草在线视频观看| 啦啦啦在线免费观看视频4| 在线天堂中文资源库| 亚洲欧美日韩高清在线视频 | 国产精品熟女久久久久浪| 最近最新中文字幕大全免费视频 | 99国产综合亚洲精品| 精品久久久精品久久久| 国产精品.久久久| 亚洲国产最新在线播放| 岛国毛片在线播放| 91精品伊人久久大香线蕉| 大陆偷拍与自拍| 高清欧美精品videossex| 肉色欧美久久久久久久蜜桃| 下体分泌物呈黄色| 日韩av在线免费看完整版不卡| 秋霞在线观看毛片| 欧美人与性动交α欧美精品济南到| 亚洲精品国产av蜜桃| 亚洲欧洲国产日韩| 在线观看一区二区三区激情| 一级毛片电影观看| 亚洲精品美女久久久久99蜜臀 | 欧美变态另类bdsm刘玥| 亚洲天堂av无毛| 中文字幕制服av| 精品久久久久久电影网| 午夜福利视频在线观看免费| 亚洲精品乱久久久久久| 涩涩av久久男人的天堂| 国产男人的电影天堂91| 男女无遮挡免费网站观看| 两个人看的免费小视频| 亚洲精品美女久久av网站| 一区二区三区乱码不卡18| 精品国产乱码久久久久久男人| 国产精品久久久久成人av| 精品福利永久在线观看| 亚洲欧美日韩另类电影网站| 亚洲一卡2卡3卡4卡5卡精品中文|