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

    基于原子磁力計的穿戴式腦磁圖動態(tài)測量研究

    2022-09-09 08:26:42陳春巧郭清乾徐佳玉馮曉宇楊曉冬
    波譜學雜志 2022年3期
    關(guān)鍵詞:磁力計被試者梯度

    陳春巧,張 欣,郭清乾,徐佳玉,馮曉宇,常 嚴,胡 濤#,楊曉冬*

    基于原子磁力計的穿戴式腦磁圖動態(tài)測量研究

    陳春巧1,2,3,張 欣2,3,郭清乾2,3,徐佳玉1,2,馮曉宇2,3,常 嚴2,3,胡 濤2,3#,楊曉冬1,2,3*

    1. 長春理工大學 電子信息工程學院,吉林 長春 130022;2. 中國科學院蘇州生物醫(yī)學工程技術(shù)研究所,江蘇 蘇州 215163;3. 季華實驗室,廣東 佛山 528200

    腦磁圖作為一種無創(chuàng)的腦功能成像技術(shù),依靠超高的時間及空間溯源分辨率,在腦科學研究和臨床應用領(lǐng)域中有著極其重要的價值.本文介紹了自主搭建的基于原子磁力計的穿戴式腦磁圖系統(tǒng),通過設計勻場補償線圈組并結(jié)合參考傳感器陣列,實現(xiàn)被試頭部運動區(qū)域內(nèi)剩磁在±1 nT以內(nèi),保證動態(tài)測量過程中傳感器輸出維持在動態(tài)范圍以內(nèi);同時提出了一種虛擬合成梯度去噪方法,顯著抑制了環(huán)境共模噪聲;最終在被試者頭部自然運動狀態(tài)下,成功檢測到高信噪比的節(jié)律信號與聽覺誘發(fā)磁場信號,證實了該系統(tǒng)的有效性,為穿戴式腦磁圖應用推廣提供更多的可能性.

    穿戴式腦磁圖;原子磁力計;動態(tài)測量;勻場補償;噪聲抑制

    引 言

    近年來,一種新型的量子弱磁傳感器——無自旋交換弛豫(spin-exchange relaxation-free,SERF)原子磁力計(optically-pumped magnetometer,OPM)得到了快速發(fā)展[13-15],其利用堿金屬原子光泵浦自旋極化和零場磁共振檢測原理,通過抑制自旋交換弛豫,它無需工作在低溫環(huán)境就能實現(xiàn)和SQUID相當?shù)奶綔y靈敏度,并且具有輕便、可小型化、探測距離近等優(yōu)勢[16].通過陣列形式將小型原子磁力計排布在腦磁測量帽插槽內(nèi)可實現(xiàn)穿戴式腦磁圖——OPM-MEG,并可在被試者自然運動狀態(tài)下進行測量,大大提升了系統(tǒng)的靈活性和便利性[17,18].2018年,英國諾丁漢大學團隊首次利用OPM-MEG實現(xiàn)了被試者運動過程中波(13~30 Hz)腦磁信號測量[19].2019年,Barry等[20]利用OPM-MEG實現(xiàn)了過程中海馬體波(4~8 Hz)信號測量動態(tài).2021年,Seymour等[21]利用OPM-MEG開展了動態(tài)過程中聽覺誘發(fā)信號的檢測.輕量化穿戴式腦磁圖動態(tài)測量的實現(xiàn)在提升被試者體驗感的同時擴大了其臨床應用范圍,尤其是對癲癇、多動癥、帕金森病等無法控制自身活動的患者,彌補了傳統(tǒng)腦磁圖技術(shù)的不足.

    本文詳細介紹了自主搭建的基于原子磁力計的穿戴式腦磁圖系統(tǒng),為實現(xiàn)動態(tài)腦磁信號測量,設計了勻場補償線圈組并結(jié)合參考傳感器陣列,用于被試頭部運動區(qū)域內(nèi)剩磁補償,保證動態(tài)測量過程中傳感器輸出維持在動態(tài)范圍內(nèi).此外,研究了一種虛擬合成梯度去噪方法,用于抑制環(huán)境共模噪聲.最后,設計了兩組動態(tài)腦磁信號測量實驗,用于驗證系統(tǒng)的可靠性和有效性.

    1 穿戴式腦磁圖系統(tǒng)

    1.1 系統(tǒng)整體介紹

    本文自主搭建的基于原子磁力計的穿戴式腦磁圖系統(tǒng)基本結(jié)構(gòu)如圖1(a)所示,其中OPM傳感器均使用美國QuSpin公司二代原子磁力計(QZFM Gen-2.0),該OPM最多支持雙通道(軸與軸)腦磁信號采集,其測量方向如圖1(b)所示.系統(tǒng)主體在磁屏蔽室中,磁屏蔽室內(nèi)部尺寸為 1.75 m×1.95 m×2.23 m,可將地磁場屏蔽至10 nT以內(nèi),使傳感器處于正常工作狀態(tài),同時利用左右兩側(cè)勻場補償線圈組進一步補償屏蔽房中心區(qū)域剩磁.被試頭部佩戴柔性腦磁測量帽,表面共排布32個插槽用于放置OPM探測傳感器陣列,測量帽可根據(jù)任意頭型調(diào)整大小,并使傳感器緊貼頭皮. 4個OPMs作為參考傳感器陣列固定在左右兩側(cè)支架上,用于測量屏蔽房內(nèi)部剩磁及環(huán)境噪聲. OPM參考傳感器陣列和探測傳感器陣列信號經(jīng)屏蔽房外的OPM電子學系統(tǒng),輸出到控制采集系統(tǒng),進行數(shù)據(jù)處理.被試接受腦磁誘發(fā)刺激模塊產(chǎn)生的聽覺、視覺或觸覺等相關(guān)刺激[22],并由腦磁誘發(fā)刺激控制系統(tǒng)所控制.控制采集系統(tǒng)采用C++語言編寫,以1 024 Hz采樣率記錄原始數(shù)據(jù)并與刺激序列保持同步.

    圖1 (a)自主搭建的基于原子磁力計(OPM)的穿戴式腦磁圖儀系統(tǒng)整機圖;(b) OPM傳感器雙軸測量方向

    1.2 勻場補償

    由于磁屏蔽室內(nèi)部殘余剩磁的大小約為10 nT,且在空間中存在較大的梯度場,在OPM-MEG動態(tài)實驗過程中,頭部的輕微晃動都有可能超出OPM傳感器自身的動態(tài)范圍(±1.5 nT),使其無法正確探測腦磁信號.因此,如何使用勻場補償線圈組進一步補償磁屏蔽室內(nèi)部剩磁及梯度場至關(guān)重要[19,23]. 本文基于平面型載流面[24]自主設計了7組對稱補償線圈組(3組勻場線圈B、BB,4組一階梯度線圈dB/d、dB/d、dB/d和dB/d),補償中心區(qū)域為40 cm×40 cm×40 cm,7組補償線圈張疊加粘貼在左右兩側(cè)固定支架上,線圈尺寸均為165 cm×140 cm,左右兩側(cè)線圈間距120 cm.7組補償線圈系數(shù)實測值如表1所示.

    表1 補償線圈系數(shù)實測值

    采用4個OPMs作為參考傳感器陣列擺放在被試頭部左右兩側(cè),兩側(cè)間距為40 cm.具體位置如圖2(a)所示,4個OPMs參考傳感器圍成一個30 cm×30 cm×40 cm區(qū)域,可涵蓋被試頭部自然運動范圍.OPMs參考傳感器均工作在三軸直流測量(field-zero)模式,實時探測所在位置處三軸方向剩磁,參考傳感器位置擺放包含各個方向梯度磁場(例如1號與2號參考傳感器作差可得到方向的梯度場),便于最優(yōu)化補償.勻場補償需建立優(yōu)化目標函數(shù),目標函數(shù)定義為:

    為測試勻場補償效果,我們將一個OPM探測傳感器插入腦磁測量帽中,傳感器工作在單軸弱磁測量模式.被試頭部在30 cm×30 cm×40 cm區(qū)域內(nèi)各個方向自由晃動,測量補償前后剩磁變化,結(jié)果如圖2(b)所示.可以看出,未施加勻場補償情況下,被試頭動極易導致被測磁場(細虛線)超出傳感器動態(tài)范圍(±1 500 pT)達到飽和,無法正常進行腦磁測量;而在完成勻場補償后(粗實線),被試頭部任意晃動,被測磁場仍在±1 000 pT(即±1 nT)以內(nèi),傳感器可正常工作.

    圖2 (a) OPM參考傳感器位置擺放示意圖;(b)被試自然頭動狀態(tài)下,探測傳感器所測磁場變化(細虛線:未施加勻場補償;粗實線:施加勻場補償)

    1.3 虛擬合成梯度去噪

    環(huán)境噪聲的抑制將有利于獲取高信噪比的腦磁信號,本文提出一種虛擬合成梯度去噪方法用于環(huán)境共模噪聲抑制.該系統(tǒng)配置由腦磁測量帽上的OPMs探測傳感器陣列和遠離頭皮的參考傳感器陣列組成,參考傳感器陣列擺放位置與圖2(a)保持一致.在完成勻場補償后,4個參考傳感器工作在單軸弱磁測量模式,虛擬合成梯度去噪具體方法如下:

    圖3 (a)頭模照片;(b)空屏蔽房中,3通道傳感器的虛擬合成梯度去噪結(jié)果

    2 腦磁圖動態(tài)測量

    為檢驗穿戴式腦磁圖系統(tǒng)的可靠性和實用性,我們開展了兩組被試者自然頭動狀態(tài)下的腦磁測量實驗,分別是節(jié)律光刺激和聽覺誘發(fā)刺激,每組實驗被試者4人,實驗設計均相同,節(jié)律實驗重復三次,聽覺實驗重復兩次,本文展示其中一名被試實驗結(jié)果[其余三名被試者實驗結(jié)果見圖S1和圖S2(掃描文章首頁二維碼或在論文網(wǎng)頁版查看)].在兩組實驗中,被試頭部在30 cm×30 cm×40 cm區(qū)域內(nèi)自由擺動.腦磁圖動態(tài)測量實驗開始前,參考OPMs傳感器陣列在三軸直流測量模式下工作用于勻場補償.補償完畢后,補償線圈組電流保持恒定;此時,參考傳感器陣列切換到單軸弱磁測量模式(軸)測量環(huán)境噪聲,用于梯度去噪,梯度去噪覆蓋整個動態(tài)腦磁測量過程并與探測傳感器保持同步.

    對于節(jié)律光刺激實驗,5通道OPMs探測傳感器配置在被試者的枕骨區(qū)域,因為節(jié)律(8~13 Hz范圍)信號源主要位于枕骨區(qū)[4,25,26].光刺激任務設計如圖4所示,被試者通過空氣耳機指示在整個實驗中每5 s交替睜眼和閉眼.傳感器測量得到的全頻段時域信號經(jīng)過虛擬合成梯度去噪及8~13 Hz帶通濾波后,結(jié)果如圖5(a)所示,灰色和黃色區(qū)域代表睜眼和閉眼的狀態(tài).從圖中可以看出5通道時域數(shù)據(jù)均有明顯的強弱交替的節(jié)律信號變化,當被試者睜眼時,信號幅度明顯減弱.各通道間信號幅度強弱則取決于OPM傳感器與顱內(nèi)信號源相對距離,越靠近信號源則幅度越強.

    圖4 a節(jié)律光刺激實驗設計

    聽覺信號源主要位于大腦皮層顳葉,因此聽覺誘發(fā)磁場信號測量的實驗中將10通道OPMs探測傳感器放置在被試者頭部的左側(cè)和右側(cè)顳葉區(qū)域,兩側(cè)各5通道.被試者佩戴空氣耳機聽取來自刺激裝置產(chǎn)生的聲音刺激,實驗總共施加400次聽覺刺激,持續(xù)時間為0.3 s,刺激時間間隔為1.7 s.同時,聽覺刺激頻率隨機呈現(xiàn)1 000 Hz或1 200 Hz,以保證被試不易產(chǎn)生聽覺疲勞,其中1 000 Hz刺激的數(shù)量占刺激總數(shù)的80%.聽覺誘發(fā)磁場信號所處波段為2~40 Hz,我們同樣采用虛擬合成梯度去噪方法對原始信號進行噪聲抑制并使用帶通濾波截取2~40 Hz范圍內(nèi)信號,對400個聽覺刺激信號做疊加取平均,得到如圖5(b)所示的聽覺誘發(fā)磁場時域信號,正負信號反映兩側(cè)腦磁信號方向相對傳感器測量方向相反.從圖中可以看出,在刺激產(chǎn)生大約100 ms后,10通道OPM傳感器均出現(xiàn)明顯的磁場強度峰值,該峰值即為M100信號.M100是一個典型的聽覺刺激誘發(fā)反應尖峰信號,產(chǎn)生于聽覺刺激發(fā)生100 ms以后,該反應是一個瞬態(tài)過程,后續(xù)不再有明顯反應[9,27].

    圖5 (a)被試者自然頭動狀態(tài)下5通道a節(jié)律信號時域圖;(b)被試自然頭動狀態(tài)下10通道聽覺誘發(fā)磁場信號時域圖

    兩組實驗均表明該系統(tǒng)可實現(xiàn)高信噪比的腦磁圖動態(tài)測量,證實了該系統(tǒng)的有效性和實用性.

    3 結(jié)論

    本文詳細描述了自主搭建的基于原子磁力計的穿戴式腦磁圖動態(tài)測量系統(tǒng),通過配置參考傳感器陣列和勻場補償線圈組,結(jié)合相應的補償算法,可實現(xiàn)30 cm×30 cm×40 cm區(qū)域內(nèi)剩磁在±1 nT以內(nèi),保證動態(tài)測量過程中傳感器輸出維持在動態(tài)范圍以內(nèi).同時本文提出了一種虛擬合成梯度去噪方法,實現(xiàn)了環(huán)境共模噪聲的有效抑制,使信號本底噪聲可降低至傳感器靈敏度附近.最后設計了兩組實驗用于驗證穿戴式腦磁圖系統(tǒng)動態(tài)測量性能,結(jié)果表明被試在頭部自然運動狀態(tài)下,可以成功檢測到高信噪比的節(jié)律信號與聽覺誘發(fā)磁場信號,證實了該系統(tǒng)的有效性.該研究為穿戴式腦磁圖在癲癇、多動癥、帕金森病等患者動態(tài)測量的應用推廣提供了更多的可能性.

    圖S1 其余三名被試自然頭動狀態(tài)下5通道節(jié)律信號時域圖

    圖S2 其余三名被試自然頭動狀態(tài)下10通道聽覺誘發(fā)磁場信號時域圖

    [1] YANG L Q, LIN F C, LEI H. Resting state functional connectivity in brain studied by fMRI approach[J]. Chinese J Magn Reson, 2010, 27(3): 326-340.

    楊麗琴, 林富春, 雷皓. 靜息狀態(tài)下腦功能連接的磁共振成像研究[J]. 波譜學雜志, 2010, 27(3): 326-340.

    [2] CHENG L W, WANG L L, ZHONG K. Application of fMRI in transcranial direct current stimulation researches[J]. Chinese J Magn Reson, 2020, 37(4): 533-546.

    程力維, 王璐璐, 鐘凱. fMRI在經(jīng)顱直流電刺激研究中的應用進展[J]. 波譜學雜志, 2020, 37(4): 533-546.

    [3] VARONE G, HUSSAIN Z, SHEIKH Z, et al. Real-time artifacts reduction during TMS-EEG co-registration: A comprehensive review on technologies and procedures[J]. Sensors (Basel, Switzerland), 21(2): 637.

    [4] COHEN D. Magnetoencephalography: detection of the brain's electrical activity with a superconducting magnetometer[J]. Science, 1972, 175(4022): 664-666.

    [5] COHEN. Magnetoencephalography: evidence of magnetic fields produced by alpha-rhythm currents[J]. Science, 1968, 161(3843): 784-786.

    [6] GRATTA C D, PIZZELLA V, TECCHIO F, et al. Magnetoencephalography - a noninvasive brain imaging method with 1 ms time resolution[J]. Rep Prog Phys, 2001, 64(12): 1759-1814.

    [7] VRBA J. Multichannel SQUID biomagnetic systems. applications of superconductivity[M]. Springer Netherlands, 2000. 61-138.

    [8] CHEN M J, LIAO S X,YANG H C, et al. Nuclear magnetic resonance and imaging of hyperpolarized3He using high-c superconducting quantum interference device in microtesla magnetic fields[J]. Chinese J Magn Reson, 2010. 27(3): 386-395.

    陳名杰, 廖書賢, 楊鴻昌, 等. 采用超導量子干涉組件在微特斯拉磁場下獲取超極化3He的核磁共振波譜和影像[J]. 波譜學雜志, 2010, 27(3): 386-395.

    [9] H?M?L?INEN M, HARI R, ILMONIEMI R J, et al. Magnetoencephalography—theory, instrumentation, and applications to noninvasive studies of the working human brain[J]. Rev Mod Phys, 1993, 65(2): 413-497.

    [10] WEINSTOCK H. SQUID sensors: fundamentals, fabrication and applications[J]. Springer Netherlands, 1996.

    [11] VRBA J. Magnetoencephalography: The art of finding a needle in a haystack[J]. J Psychophysiol, 2003, 17(4): 237-237.

    [12] MUKAMEL R, GELBARD H, ARIELI A, et al. Coupling between neuronal firing, field potentials, and FMRI in human auditory cortex[J]. Science, 2005, 309(5736): 951-954.

    [13] WANG X F, SUN X P, ZHAO X C, et al. Progress in biomagnetic signal measurements with ultra-sensitive atomic magnetometers[J]. Chinese Journal of Lasers, 2018, 45(2): 0207012.

    王曉飛, 孫獻平, 趙修超等. 超靈敏原子磁力計在生物磁應用中的研究進展[J]. 中國激光, 2018, 45(2): 0207012

    [14] MA H F, WU Y T, ZHAO W, et al. Research progress of magnetoencephalography in the functional mechanism of bilingual brain[J]. Chinese Journal of Biomedical Engineering, 2021, 40(4): 477-484.

    馬恒芬, 吳云濤, 趙文, 等. 雙語腦功能機制的腦磁圖研究進展[J]. 中國生物醫(yī)學工程學報, 2021, 40(4): 477-484.

    [15] XU W J, JIANG M, PENG X H. Study on ultra-low-field nuclear magnetic resonance spectroscopy based on high-sensitivity atomic magnetometer[J]. Journal of University of Science and Technology of China, 2020, 50(8): 1138-1143.

    徐文杰, 江敏, 彭新華. 基于高靈敏度原子磁力計的超低場核磁共振譜學研究[J]. 中國科學技術(shù)大學學報, 2020, 50(8): 1138-1143.

    [16] ZHANG S L, CAO N. A synthetic optically pumped gradiometer for magnetocardiography measurements[J]. Chinese Physics B, 2020, 29(4).

    [17] BOTO E, MEYER S S, SHAH V, et al. A new generation of magnetoencephalography: Room temperature measurements using optically-pumped magnetometers[J]. Neuroimage, 2017, 149: 404-414.

    [18] ZHANG X, CHEN C Q, ZHANG M K, et al. Detection and analysis of MEG signals in occipital region with double-channel OPM sensors[J]. J Neurosci Methods, 2020, 346: 108948.

    [19] BOTO E, HOLMES N, LEGGETT J, et al. Moving magnetoencephalography towards real-world applications with a wearable system[J]. Nature, 2018, 555(7698): 657-661.

    [20] BARRY D N, TIERNEY T M, HOLMES N, et al. Imaging the human hippocampus with optically-pumped magnetoencephalography[J]. Neuroimage, 2019, 203: 116192.

    [21] SEYMOUR R A, ALEXANDER N, MELLOR S, et al. Using OPMs to measure neural activity in standing, mobile participants[J]. Neuroimage, 2021, 244: 118604.

    [22] SUN W, WANG H, ZHANG Y, et al. Optimal design for quantification of gas concentration based olfactory stimulator[J]. Chinese J Magn Reson, 2021, 38(1): 12-21.

    孫韋, 王慧, 張寅, 等. 基于氣體濃度定量的嗅覺刺激器優(yōu)化設計[J]. 波譜學雜志, 2021, 38(1): 12-21.

    [23] IIVANAINEN J, ZETTER R, GRON M, et al. On-scalp MEG system utilizing an actively shielded array of optically-pumped magnetometers[J]. Neuroimage, 2019, 194: 244-258.

    [24] HOLMES N, LEGGETT J, BOTO E, et al. A bi-planar coil system for nulling background magnetic fields in scalp mounted magnetoencephalography[J]. Neuroimage, 2018, 181: 760-774.

    [25] KAWABATA N. Nonstationary power spectrum analysis of the photic alpha blocking[J]. Kybernetik, 1972, 12(1): 40-44.

    [26] BRENNER D, WILLIAMSON S J, KAUFMAN L. Visually evoked magnetic fields of the human brain[J]. Science, 1975, 190(4213): 480-482.

    [27] JOHNSON C N, SCHWINDT P D, WEISEND M. Multi-sensor magnetoencephalography with atomic magnetometers[J]. Phys Med Biol, 2013, 58(17): 6065-6077.

    Moving Wearable Magnetoencephalography Measurement Study Based on Optically-pumped Magnetometer

    1,2,3,2,3,2,3,1,2,2,3,2,3,2,3#,1,2,3*

    1.School of Electronic and Information Engineering, Changchun University of Science and Technology, Changchun 130022, China; 2. Suzhou Institute of Biomedical Engineering and Technology, Chinese Academy of Sciences, Suzhou 215163, China; 3. Jihua Laboratory, Foshan 528200, China

    Magnetoencephalography is a non-invasive technology for brain function imaging, which is of enormous value to brain science research and clinical application due to its ultra-high temporal and spatial trace resolution. In this paper, we introduce a self-built and atomic magnetometer based wearable magnetoencephalography system. By designing bi-planar coils system and combining with reference sensor array, the residual magnetic field in the subject’s head movement area is controlled to be within±1 nT, which ensures the sensors are maintained within their dynamic range during the moving measurement. At the same time, a virtual gradiometer-based noise reduction method is proposed to suppress the common-mode magnetic-field noise. Finally, the alpha rhythm and auditory evoked magnetic field signals with high signal-to-noise ratio are successfully detected under the subject’s natural head movement and the effectiveness of the system is confirmed. This study could provide more possibilities for the application and promotion of moving wearable magnetoencephalography.

    wearablemagnetoencephalography, optically-pumped magnetometer, moving measurement, field nulling, noise reduction

    O441.5;Q64

    A

    10.11938/cjmr20222975

    2022-02-15;

    2022-03-25

    蘇州市基礎(chǔ)研究試點項目(SJC2021024);季華實驗室項目(X190131TD190);江蘇省自然科學基金青年項目(BK20200215).

    # Tel: 17625326300, E-mail: hutao@sibet.ac.cn;

    * Tel: 18900616030, E-mail: xiaodong.yang@sibet.ac.cn.

    猜你喜歡
    磁力計被試者梯度
    一個改進的WYL型三項共軛梯度法
    基于遞推最小二乘法加速度計信息輔助的磁力計標定方法?
    六軸IMU補償?shù)拇帕τ媱討B(tài)穩(wěn)定校準
    一種自適應Dai-Liao共軛梯度法
    德西效應
    基于低功耗藍牙的磁力計在線校準系統(tǒng)設計
    德西效應
    一類扭積形式的梯度近Ricci孤立子
    兒童花生過敏治療取得突破
    為什么我們會忍不住發(fā)朋友圈?
    奧秘(2016年3期)2016-03-23 21:58:57
    这个男人来自地球电影免费观看| 久久欧美精品欧美久久欧美| 丰满的人妻完整版| 老鸭窝网址在线观看| 免费不卡黄色视频| 国产伦一二天堂av在线观看| 欧美av亚洲av综合av国产av| 18禁美女被吸乳视频| 在线观看一区二区三区| 国产精品精品国产色婷婷| av福利片在线| 男男h啪啪无遮挡| 麻豆久久精品国产亚洲av| 在线观看免费视频网站a站| av有码第一页| 亚洲国产精品成人综合色| 老鸭窝网址在线观看| 亚洲全国av大片| 色尼玛亚洲综合影院| 久久久久久久久免费视频了| 午夜福利成人在线免费观看| 在线观看日韩欧美| 国产av又大| 午夜亚洲福利在线播放| 麻豆av在线久日| 夜夜爽天天搞| av欧美777| 成人三级黄色视频| 久久人人爽av亚洲精品天堂| 高潮久久久久久久久久久不卡| 人人妻人人澡人人看| 男女之事视频高清在线观看| 夜夜爽天天搞| 一级,二级,三级黄色视频| 黑丝袜美女国产一区| 亚洲片人在线观看| 一级a爱片免费观看的视频| 免费在线观看完整版高清| 午夜福利一区二区在线看| АⅤ资源中文在线天堂| 亚洲男人天堂网一区| 黄色 视频免费看| 美女大奶头视频| 欧美午夜高清在线| 波多野结衣高清无吗| 亚洲中文字幕日韩| 国产亚洲精品av在线| 久久伊人香网站| 久久 成人 亚洲| 一边摸一边抽搐一进一出视频| 成人18禁高潮啪啪吃奶动态图| 成人欧美大片| 久久性视频一级片| 99riav亚洲国产免费| 高清黄色对白视频在线免费看| 91国产中文字幕| 黄色视频不卡| 99在线人妻在线中文字幕| 十八禁网站免费在线| 欧美激情 高清一区二区三区| 久久精品国产亚洲av香蕉五月| 国产一区在线观看成人免费| 国产成人精品在线电影| 久久精品国产99精品国产亚洲性色 | 日韩欧美一区视频在线观看| 高清黄色对白视频在线免费看| 婷婷丁香在线五月| 精品卡一卡二卡四卡免费| 老司机在亚洲福利影院| 村上凉子中文字幕在线| 日韩欧美免费精品| 日韩有码中文字幕| 精品国产乱子伦一区二区三区| 色哟哟哟哟哟哟| 亚洲精品久久成人aⅴ小说| 一进一出抽搐gif免费好疼| 91字幕亚洲| 久热爱精品视频在线9| 久久欧美精品欧美久久欧美| 波多野结衣av一区二区av| 精品乱码久久久久久99久播| 黑丝袜美女国产一区| 国产成人精品久久二区二区91| 变态另类丝袜制服| 国产私拍福利视频在线观看| 国产精品二区激情视频| 性色av乱码一区二区三区2| 亚洲人成77777在线视频| 嫁个100分男人电影在线观看| 啪啪无遮挡十八禁网站| 桃色一区二区三区在线观看| 久久精品国产99精品国产亚洲性色 | 午夜福利高清视频| 他把我摸到了高潮在线观看| 亚洲男人天堂网一区| 日韩大码丰满熟妇| 国产成人免费无遮挡视频| 国产午夜福利久久久久久| АⅤ资源中文在线天堂| 两人在一起打扑克的视频| 精品久久久精品久久久| 欧美+亚洲+日韩+国产| 国产男靠女视频免费网站| 天堂动漫精品| 欧美不卡视频在线免费观看 | 欧美av亚洲av综合av国产av| 啦啦啦免费观看视频1| 久久天堂一区二区三区四区| 国产成人精品无人区| 欧洲精品卡2卡3卡4卡5卡区| 欧美激情久久久久久爽电影 | 可以免费在线观看a视频的电影网站| 亚洲一区二区三区不卡视频| 亚洲av成人不卡在线观看播放网| 涩涩av久久男人的天堂| 无遮挡黄片免费观看| 一区二区日韩欧美中文字幕| 亚洲片人在线观看| 午夜免费观看网址| 欧美黄色片欧美黄色片| 亚洲一区二区三区不卡视频| 九色国产91popny在线| 两个人看的免费小视频| 国产激情久久老熟女| 亚洲一区高清亚洲精品| 校园春色视频在线观看| 色尼玛亚洲综合影院| 日日夜夜操网爽| 国产99久久九九免费精品| 国产亚洲精品综合一区在线观看 | 少妇粗大呻吟视频| 十八禁网站免费在线| 国产成人影院久久av| 日韩欧美在线二视频| 欧美绝顶高潮抽搐喷水| 美女大奶头视频| 午夜免费鲁丝| 免费在线观看视频国产中文字幕亚洲| 亚洲国产毛片av蜜桃av| 大型av网站在线播放| 国产成人欧美在线观看| 国产乱人伦免费视频| 天天躁夜夜躁狠狠躁躁| 亚洲狠狠婷婷综合久久图片| 十分钟在线观看高清视频www| 精品久久久精品久久久| 在线视频色国产色| 成在线人永久免费视频| 亚洲免费av在线视频| 精品国产国语对白av| 婷婷精品国产亚洲av在线| 亚洲无线在线观看| 成年女人毛片免费观看观看9| 看黄色毛片网站| 此物有八面人人有两片| 久久中文字幕人妻熟女| 国产欧美日韩综合在线一区二区| 久久久久久久久免费视频了| 不卡一级毛片| 亚洲中文av在线| 99香蕉大伊视频| 韩国av一区二区三区四区| 久久精品aⅴ一区二区三区四区| 亚洲国产精品合色在线| 精品卡一卡二卡四卡免费| 亚洲国产毛片av蜜桃av| 久久九九热精品免费| 色综合婷婷激情| 午夜a级毛片| 波多野结衣一区麻豆| 无限看片的www在线观看| 亚洲激情在线av| 黄色片一级片一级黄色片| 国产精品野战在线观看| 欧美亚洲日本最大视频资源| 国产91精品成人一区二区三区| 欧美人与性动交α欧美精品济南到| 国产亚洲精品一区二区www| av视频免费观看在线观看| 可以在线观看的亚洲视频| 国产精品乱码一区二三区的特点 | 久久伊人香网站| 午夜精品国产一区二区电影| 1024视频免费在线观看| 91精品国产国语对白视频| 国产一区在线观看成人免费| 人人妻人人澡人人看| 免费无遮挡裸体视频| 变态另类丝袜制服| 级片在线观看| 久久国产精品男人的天堂亚洲| 禁无遮挡网站| 欧美日韩精品网址| 国内毛片毛片毛片毛片毛片| 国产激情欧美一区二区| 亚洲激情在线av| aaaaa片日本免费| 午夜激情av网站| 国产午夜精品久久久久久| 91麻豆av在线| 亚洲伊人色综图| 1024视频免费在线观看| 在线国产一区二区在线| 深夜精品福利| 色尼玛亚洲综合影院| 欧美+亚洲+日韩+国产| 老鸭窝网址在线观看| 99国产精品99久久久久| 欧美激情 高清一区二区三区| e午夜精品久久久久久久| 亚洲人成网站在线播放欧美日韩| 一夜夜www| 亚洲欧美日韩高清在线视频| 日韩欧美一区视频在线观看| 欧美中文日本在线观看视频| 久久人人97超碰香蕉20202| 夜夜爽天天搞| 欧美乱色亚洲激情| 国产精品日韩av在线免费观看 | 国产欧美日韩精品亚洲av| 久久久久国产一级毛片高清牌| 老汉色∧v一级毛片| 最近最新中文字幕大全免费视频| 91在线观看av| 欧美激情极品国产一区二区三区| 精品第一国产精品| 一区二区日韩欧美中文字幕| 少妇的丰满在线观看| 不卡一级毛片| www.www免费av| 国产精品亚洲美女久久久| 男女床上黄色一级片免费看| 成人18禁在线播放| 欧美日本亚洲视频在线播放| 久9热在线精品视频| 午夜激情av网站| 免费高清在线观看日韩| 啦啦啦韩国在线观看视频| 久久国产精品男人的天堂亚洲| 别揉我奶头~嗯~啊~动态视频| 亚洲色图av天堂| 国产三级黄色录像| 精品一区二区三区av网在线观看| 国产精品久久久久久精品电影 | 欧美 亚洲 国产 日韩一| www.999成人在线观看| 精品久久蜜臀av无| 人人妻人人澡欧美一区二区 | 91成人精品电影| 精品欧美国产一区二区三| 老熟妇乱子伦视频在线观看| 欧美日韩乱码在线| 在线观看舔阴道视频| 日本精品一区二区三区蜜桃| 午夜久久久久精精品| 91成人精品电影| 麻豆av在线久日| 国产成+人综合+亚洲专区| 91精品三级在线观看| 午夜久久久在线观看| 婷婷精品国产亚洲av在线| 亚洲成av人片免费观看| 国产亚洲精品久久久久久毛片| 搡老妇女老女人老熟妇| 一本综合久久免费| 黑丝袜美女国产一区| 人人澡人人妻人| 美女免费视频网站| 一级作爱视频免费观看| 日韩中文字幕欧美一区二区| 18禁美女被吸乳视频| 成人精品一区二区免费| 天堂动漫精品| 国产精品一区二区免费欧美| 国产私拍福利视频在线观看| 成人三级做爰电影| 搡老岳熟女国产| 99久久国产精品久久久| 精品电影一区二区在线| 中文字幕人妻熟女乱码| 丝袜在线中文字幕| 久久久久九九精品影院| 一二三四在线观看免费中文在| 欧美日韩亚洲综合一区二区三区_| 国产99久久九九免费精品| 在线永久观看黄色视频| 国产精品久久久久久人妻精品电影| 免费一级毛片在线播放高清视频 | 成人精品一区二区免费| 国产免费av片在线观看野外av| 亚洲成a人片在线一区二区| 一进一出抽搐gif免费好疼| 午夜福利视频1000在线观看 | 欧美日韩一级在线毛片| 每晚都被弄得嗷嗷叫到高潮| 波多野结衣一区麻豆| 精品少妇一区二区三区视频日本电影| a级毛片在线看网站| 久久精品亚洲熟妇少妇任你| 黄频高清免费视频| 欧美绝顶高潮抽搐喷水| 国产精品自产拍在线观看55亚洲| 9191精品国产免费久久| 久久久精品欧美日韩精品| 亚洲 欧美一区二区三区| 国产私拍福利视频在线观看| 午夜久久久久精精品| 国产精品香港三级国产av潘金莲| 色av中文字幕| 男女下面进入的视频免费午夜 | 国产精品久久电影中文字幕| av在线播放免费不卡| 99久久精品国产亚洲精品| 欧美亚洲日本最大视频资源| 好男人在线观看高清免费视频 | 他把我摸到了高潮在线观看| 丰满人妻熟妇乱又伦精品不卡| 国产xxxxx性猛交| 精品国产一区二区久久| 国产精品久久久久久亚洲av鲁大| 人人妻人人澡人人看| 欧美日韩黄片免| 日韩大码丰满熟妇| 一区二区三区激情视频| 日韩大码丰满熟妇| 国产亚洲av高清不卡| cao死你这个sao货| 老司机在亚洲福利影院| 精品久久久久久久毛片微露脸| 国产激情欧美一区二区| 亚洲熟妇熟女久久| 国产1区2区3区精品| 天堂√8在线中文| 久久人人爽av亚洲精品天堂| 亚洲熟妇中文字幕五十中出| 啦啦啦韩国在线观看视频| 国产精品二区激情视频| 精品一区二区三区av网在线观看| 99国产精品免费福利视频| 两个人看的免费小视频| 91精品三级在线观看| 久久久久久久久免费视频了| 长腿黑丝高跟| 国产成人免费无遮挡视频| 在线天堂中文资源库| 日韩大尺度精品在线看网址 | 亚洲国产日韩欧美精品在线观看 | 国产精品亚洲av一区麻豆| 精品人妻1区二区| 美女 人体艺术 gogo| 亚洲 国产 在线| bbb黄色大片| 亚洲专区中文字幕在线| 老司机午夜十八禁免费视频| av中文乱码字幕在线| 国内精品久久久久久久电影| 午夜影院日韩av| 亚洲免费av在线视频| av中文乱码字幕在线| 老司机午夜十八禁免费视频| 精品久久久久久,| 欧美丝袜亚洲另类 | √禁漫天堂资源中文www| 午夜a级毛片| 亚洲男人的天堂狠狠| avwww免费| 午夜久久久在线观看| 在线观看免费视频网站a站| 熟女少妇亚洲综合色aaa.| 男人舔女人下体高潮全视频| 在线播放国产精品三级| 亚洲男人天堂网一区| 人人澡人人妻人| 国产精品二区激情视频| 黄色毛片三级朝国网站| 国产av精品麻豆| 久久精品亚洲精品国产色婷小说| av视频在线观看入口| 欧美激情久久久久久爽电影 | 12—13女人毛片做爰片一| 69精品国产乱码久久久| 日韩欧美国产一区二区入口| 国产精品亚洲av一区麻豆| 天天躁狠狠躁夜夜躁狠狠躁| 精品久久久久久成人av| 1024香蕉在线观看| 久久青草综合色| 女生性感内裤真人,穿戴方法视频| 免费一级毛片在线播放高清视频 | 激情在线观看视频在线高清| 日本 欧美在线| 精品少妇一区二区三区视频日本电影| 久久九九热精品免费| 999精品在线视频| 在线观看日韩欧美| 亚洲少妇的诱惑av| 成人三级做爰电影| 久久精品国产清高在天天线| 中文字幕久久专区| 成人免费观看视频高清| 中文字幕精品免费在线观看视频| 亚洲黑人精品在线| 一区二区三区激情视频| 操出白浆在线播放| 国产亚洲精品第一综合不卡| 亚洲av电影不卡..在线观看| 1024香蕉在线观看| 亚洲精品av麻豆狂野| 露出奶头的视频| av电影中文网址| 久久精品国产亚洲av高清一级| 一边摸一边抽搐一进一出视频| 久久天躁狠狠躁夜夜2o2o| 亚洲欧美激情在线| 免费在线观看黄色视频的| 国产精品久久久久久亚洲av鲁大| 狂野欧美激情性xxxx| 久久狼人影院| 久久热在线av| 亚洲精品国产精品久久久不卡| 国产亚洲精品综合一区在线观看 | 亚洲无线在线观看| 不卡av一区二区三区| 电影成人av| 一区二区三区激情视频| 如日韩欧美国产精品一区二区三区| 一区二区三区高清视频在线| 国产97色在线日韩免费| 午夜a级毛片| 91大片在线观看| 免费av毛片视频| 夜夜躁狠狠躁天天躁| 国产精品日韩av在线免费观看 | 久久人妻熟女aⅴ| 精品福利观看| 男人操女人黄网站| 不卡一级毛片| 久久这里只有精品19| 神马国产精品三级电影在线观看 | 精品无人区乱码1区二区| 亚洲自拍偷在线| 欧洲精品卡2卡3卡4卡5卡区| 99国产精品免费福利视频| 99国产精品99久久久久| 狠狠狠狠99中文字幕| 久久国产精品影院| 精品熟女少妇八av免费久了| 日韩成人在线观看一区二区三区| 亚洲专区国产一区二区| 给我免费播放毛片高清在线观看| 老熟妇仑乱视频hdxx| 国产亚洲欧美98| 亚洲精品久久成人aⅴ小说| 国产精品一区二区免费欧美| or卡值多少钱| 一边摸一边抽搐一进一出视频| 亚洲人成电影观看| www.自偷自拍.com| 欧美日韩亚洲国产一区二区在线观看| 深夜精品福利| 好看av亚洲va欧美ⅴa在| 日韩av在线大香蕉| 身体一侧抽搐| 日韩成人在线观看一区二区三区| 黄色丝袜av网址大全| 十八禁网站免费在线| 91在线观看av| 色精品久久人妻99蜜桃| ponron亚洲| 国产单亲对白刺激| 久久中文看片网| 99国产精品99久久久久| 欧美在线一区亚洲| 老汉色∧v一级毛片| av片东京热男人的天堂| av视频免费观看在线观看| 欧美最黄视频在线播放免费| 免费女性裸体啪啪无遮挡网站| 变态另类成人亚洲欧美熟女 | 亚洲第一av免费看| 9色porny在线观看| av在线天堂中文字幕| 色老头精品视频在线观看| 91国产中文字幕| 亚洲精品在线观看二区| 色尼玛亚洲综合影院| 日韩三级视频一区二区三区| 欧美国产精品va在线观看不卡| 嫁个100分男人电影在线观看| 自线自在国产av| 欧美日韩亚洲国产一区二区在线观看| 天堂动漫精品| 亚洲自偷自拍图片 自拍| 国产成人影院久久av| 伊人久久大香线蕉亚洲五| 成人亚洲精品av一区二区| 亚洲成av片中文字幕在线观看| 欧美日本视频| 视频在线观看一区二区三区| 日韩精品青青久久久久久| 美女扒开内裤让男人捅视频| 国产黄a三级三级三级人| 搞女人的毛片| 欧美成狂野欧美在线观看| 神马国产精品三级电影在线观看 | 亚洲一码二码三码区别大吗| 亚洲五月色婷婷综合| 日本一区二区免费在线视频| 亚洲片人在线观看| 黄色毛片三级朝国网站| 欧美在线一区亚洲| ponron亚洲| 9热在线视频观看99| 麻豆国产av国片精品| 搡老岳熟女国产| 久久人人爽av亚洲精品天堂| 少妇裸体淫交视频免费看高清 | 国内精品久久久久精免费| 中文字幕另类日韩欧美亚洲嫩草| 精品国产亚洲在线| 国产单亲对白刺激| www.999成人在线观看| 亚洲自偷自拍图片 自拍| 一区二区三区精品91| 国产亚洲精品第一综合不卡| xxx96com| 久久人妻熟女aⅴ| 久久久久久免费高清国产稀缺| e午夜精品久久久久久久| 国产一卡二卡三卡精品| 精品久久久精品久久久| 成人手机av| 亚洲精品一区av在线观看| 亚洲第一电影网av| 天天添夜夜摸| 黄网站色视频无遮挡免费观看| 日韩一卡2卡3卡4卡2021年| 精品国产美女av久久久久小说| 精品卡一卡二卡四卡免费| 欧美成人一区二区免费高清观看 | 亚洲第一电影网av| videosex国产| 亚洲欧美日韩另类电影网站| x7x7x7水蜜桃| 欧美激情 高清一区二区三区| 美国免费a级毛片| 久久性视频一级片| 一级a爱片免费观看的视频| 啦啦啦观看免费观看视频高清 | 动漫黄色视频在线观看| 久久草成人影院| 天堂√8在线中文| 色老头精品视频在线观看| 操出白浆在线播放| 国产精品久久久人人做人人爽| 99国产精品99久久久久| 欧美成人一区二区免费高清观看 | 欧美中文综合在线视频| 90打野战视频偷拍视频| 久久香蕉精品热| 国产精品亚洲av一区麻豆| 十八禁网站免费在线| 成人国产综合亚洲| 国产成人一区二区三区免费视频网站| 久久天堂一区二区三区四区| 亚洲专区字幕在线| 波多野结衣一区麻豆| 日韩精品免费视频一区二区三区| 精品高清国产在线一区| 99riav亚洲国产免费| 无遮挡黄片免费观看| 日韩大尺度精品在线看网址 | 99久久久亚洲精品蜜臀av| 美女扒开内裤让男人捅视频| 国产精品爽爽va在线观看网站 | 日韩欧美国产一区二区入口| 国产日韩一区二区三区精品不卡| 亚洲av第一区精品v没综合| ponron亚洲| 亚洲精品国产色婷婷电影| 91av网站免费观看| 亚洲国产欧美日韩在线播放| √禁漫天堂资源中文www| 国产精品自产拍在线观看55亚洲| 亚洲精品在线美女| 人妻久久中文字幕网| 午夜福利免费观看在线| 日本黄色视频三级网站网址| 老司机福利观看| 久久久久久大精品| 久久香蕉精品热| 欧美一级毛片孕妇| 亚洲无线在线观看| 亚洲精品在线美女| 免费久久久久久久精品成人欧美视频| 免费高清在线观看日韩| 亚洲精品国产区一区二| 多毛熟女@视频| 日韩欧美免费精品| 少妇粗大呻吟视频| 黄色丝袜av网址大全| 国产伦人伦偷精品视频| 亚洲专区中文字幕在线| 国产日韩一区二区三区精品不卡| 最近最新中文字幕大全免费视频| 此物有八面人人有两片| 精品久久久久久久久久免费视频| av视频在线观看入口| 亚洲精品在线观看二区| 国内精品久久久久精免费| 一级,二级,三级黄色视频| 波多野结衣av一区二区av| 免费一级毛片在线播放高清视频 | 亚洲精品一区av在线观看| 脱女人内裤的视频|