• <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
    欧美日韩成人在线一区二区| 女性被躁到高潮视频| 精品国产超薄肉色丝袜足j| 国产xxxxx性猛交| 亚洲性夜色夜夜综合| 久久国产精品人妻蜜桃| 少妇的丰满在线观看| 精品熟女少妇八av免费久了| 一本久久精品| 脱女人内裤的视频| 精品第一国产精品| 欧美人与性动交α欧美软件| 国产免费一区二区三区四区乱码| 9191精品国产免费久久| 9热在线视频观看99| 大片免费播放器 马上看| 黄色怎么调成土黄色| 国产成人免费无遮挡视频| 50天的宝宝边吃奶边哭怎么回事| 日日夜夜操网爽| 黑人操中国人逼视频| 手机成人av网站| 国产视频一区二区在线看| 老司机靠b影院| 人人妻,人人澡人人爽秒播| 国产又爽黄色视频| 91麻豆精品激情在线观看国产 | 国产精品av久久久久免费| 久久性视频一级片| 天堂8中文在线网| 丰满饥渴人妻一区二区三| 少妇被粗大的猛进出69影院| 欧美亚洲 丝袜 人妻 在线| 日日夜夜操网爽| 国产精品熟女久久久久浪| 极品少妇高潮喷水抽搐| 日本a在线网址| 交换朋友夫妻互换小说| 久久这里只有精品19| 人人妻人人添人人爽欧美一区卜| 丰满饥渴人妻一区二区三| 亚洲国产精品999| 亚洲国产欧美一区二区综合| 国产一级毛片在线| 国产一卡二卡三卡精品| 欧美日韩一级在线毛片| 精品少妇内射三级| 国产深夜福利视频在线观看| videosex国产| 久久这里只有精品19| 99国产精品99久久久久| 国产日韩欧美亚洲二区| 日本a在线网址| 91麻豆av在线| 曰老女人黄片| 一区福利在线观看| 满18在线观看网站| 欧美大码av| 亚洲专区国产一区二区| 亚洲成av片中文字幕在线观看| 亚洲人成电影观看| xxxhd国产人妻xxx| 男女午夜视频在线观看| 亚洲精品乱久久久久久| 久久久久久久久免费视频了| 久久久久久久国产电影| 美女脱内裤让男人舔精品视频| 极品人妻少妇av视频| 中文字幕高清在线视频| 欧美精品高潮呻吟av久久| 天堂中文最新版在线下载| 国产成人欧美在线观看 | 日本av手机在线免费观看| 另类精品久久| 国产精品自产拍在线观看55亚洲 | 成人黄色视频免费在线看| 脱女人内裤的视频| 精品亚洲成国产av| 两个人看的免费小视频| 亚洲精品第二区| 日韩欧美一区二区三区在线观看 | 久久天堂一区二区三区四区| 汤姆久久久久久久影院中文字幕| 丁香六月欧美| 欧美av亚洲av综合av国产av| 老鸭窝网址在线观看| 免费av中文字幕在线| 少妇 在线观看| 国产不卡av网站在线观看| 久久久精品免费免费高清| 老汉色∧v一级毛片| 亚洲久久久国产精品| 色94色欧美一区二区| 免费观看av网站的网址| 国产欧美亚洲国产| 伊人亚洲综合成人网| 人妻 亚洲 视频| 成人亚洲精品一区在线观看| 中亚洲国语对白在线视频| 18禁黄网站禁片午夜丰满| 亚洲精品国产色婷婷电影| 99九九在线精品视频| 国产亚洲av高清不卡| 80岁老熟妇乱子伦牲交| 天堂中文最新版在线下载| 国产麻豆69| 久久 成人 亚洲| 亚洲三区欧美一区| 日韩视频在线欧美| 老司机影院毛片| 如日韩欧美国产精品一区二区三区| 成人三级做爰电影| 一本久久精品| 99国产极品粉嫩在线观看| 在线观看免费视频网站a站| 岛国毛片在线播放| 母亲3免费完整高清在线观看| 三级毛片av免费| 91大片在线观看| 亚洲精品一卡2卡三卡4卡5卡 | 成人亚洲精品一区在线观看| 十八禁人妻一区二区| 啦啦啦中文免费视频观看日本| 久久精品久久久久久噜噜老黄| 亚洲欧洲日产国产| 精品国产国语对白av| 婷婷成人精品国产| 国产精品香港三级国产av潘金莲| 十八禁网站网址无遮挡| 欧美老熟妇乱子伦牲交| 狠狠狠狠99中文字幕| 亚洲天堂av无毛| 亚洲中文日韩欧美视频| 国产免费视频播放在线视频| 亚洲中文字幕日韩| 欧美黑人精品巨大| 久久中文字幕一级| 91国产中文字幕| 国产精品欧美亚洲77777| 美国免费a级毛片| 亚洲精品中文字幕在线视频| 午夜福利一区二区在线看| 美女脱内裤让男人舔精品视频| 一进一出抽搐动态| 久久综合国产亚洲精品| 午夜福利视频在线观看免费| 久久久久久久久久久久大奶| 后天国语完整版免费观看| 真人做人爱边吃奶动态| 日韩欧美一区二区三区在线观看 | 久久精品aⅴ一区二区三区四区| 母亲3免费完整高清在线观看| 国产日韩一区二区三区精品不卡| 亚洲中文字幕日韩| 国产精品av久久久久免费| www.精华液| 动漫黄色视频在线观看| 久久av网站| 精品熟女少妇八av免费久了| 黄色视频在线播放观看不卡| 日韩视频一区二区在线观看| 久久综合国产亚洲精品| 国产亚洲精品第一综合不卡| 国产黄频视频在线观看| 久久精品国产综合久久久| 精品国产一区二区三区四区第35| 18禁黄网站禁片午夜丰满| 国产精品熟女久久久久浪| 宅男免费午夜| 久久久久国产一级毛片高清牌| 97精品久久久久久久久久精品| www.av在线官网国产| 女人精品久久久久毛片| 亚洲成人免费电影在线观看| 精品国产一区二区三区久久久樱花| 亚洲国产欧美在线一区| 97人妻天天添夜夜摸| 男女床上黄色一级片免费看| 在线av久久热| 丝袜在线中文字幕| 亚洲av日韩在线播放| 午夜免费成人在线视频| 电影成人av| 久久青草综合色| 精品人妻一区二区三区麻豆| a在线观看视频网站| 亚洲av片天天在线观看| 日本精品一区二区三区蜜桃| 午夜久久久在线观看| 国产一级毛片在线| 亚洲激情五月婷婷啪啪| 国产亚洲精品一区二区www | 日韩欧美一区视频在线观看| 国产成人av教育| 夜夜夜夜夜久久久久| 高清av免费在线| 国产亚洲一区二区精品| 国产精品香港三级国产av潘金莲| 女性生殖器流出的白浆| 国产精品一区二区在线观看99| 亚洲一区二区三区欧美精品| 久久人妻福利社区极品人妻图片| 黄色片一级片一级黄色片| 国产av国产精品国产| 飞空精品影院首页| 日本精品一区二区三区蜜桃| 国产色视频综合| 国产精品一区二区在线不卡| 97精品久久久久久久久久精品| 少妇被粗大的猛进出69影院| 麻豆av在线久日| 制服人妻中文乱码| 一区二区日韩欧美中文字幕| 免费看十八禁软件| 女性被躁到高潮视频| 成年人免费黄色播放视频| 亚洲第一欧美日韩一区二区三区 | 激情视频va一区二区三区| 国产精品自产拍在线观看55亚洲 | 国产精品免费大片| 国产亚洲精品第一综合不卡| 亚洲人成电影免费在线| 久久人人爽人人片av| 99国产综合亚洲精品| 黄片播放在线免费| 日本欧美视频一区| 日韩制服骚丝袜av| 国产一级毛片在线| 欧美日韩成人在线一区二区| 欧美日韩精品网址| 精品国产乱码久久久久久男人| 秋霞在线观看毛片| 久久精品人人爽人人爽视色| 淫妇啪啪啪对白视频 | 久久这里只有精品19| 欧美成人午夜精品| 91av网站免费观看| 精品卡一卡二卡四卡免费| 久久中文看片网| 久久久久视频综合| 97在线人人人人妻| 亚洲天堂av无毛| 久久久水蜜桃国产精品网| 国产精品秋霞免费鲁丝片| 99热全是精品| 欧美在线一区亚洲| 欧美人与性动交α欧美软件| 久久国产精品大桥未久av| 青青草视频在线视频观看| 男女无遮挡免费网站观看| 午夜老司机福利片| 国产欧美日韩一区二区三区在线| 免费人妻精品一区二区三区视频| www日本在线高清视频| 亚洲国产毛片av蜜桃av| 亚洲精品美女久久av网站| 香蕉国产在线看| 欧美日韩国产mv在线观看视频| 一个人免费看片子| 国内毛片毛片毛片毛片毛片| 国产成人av教育| 精品国产国语对白av| 免费不卡黄色视频| 国产高清videossex| 精品高清国产在线一区| 欧美av亚洲av综合av国产av| 啪啪无遮挡十八禁网站| 久久影院123| 国产无遮挡羞羞视频在线观看| 亚洲av美国av| 菩萨蛮人人尽说江南好唐韦庄| 免费在线观看视频国产中文字幕亚洲 | 亚洲精品乱久久久久久| 欧美日韩国产mv在线观看视频| 人人妻人人爽人人添夜夜欢视频| 电影成人av| 99热网站在线观看| 国产免费一区二区三区四区乱码| 国产精品国产三级国产专区5o| 国产精品国产av在线观看| 国产精品 国内视频| 国产日韩一区二区三区精品不卡| 久久久久久人人人人人| 国产精品一区二区在线观看99| 中文字幕色久视频| 老熟妇乱子伦视频在线观看 | 老鸭窝网址在线观看| 色视频在线一区二区三区| 亚洲av成人不卡在线观看播放网 | 超碰成人久久| 一区二区日韩欧美中文字幕| 最近最新免费中文字幕在线| 狂野欧美激情性xxxx| 丝袜脚勾引网站| 久久女婷五月综合色啪小说| 亚洲色图综合在线观看| 欧美日本中文国产一区发布| 免费在线观看日本一区| 亚洲成人国产一区在线观看| 国产人伦9x9x在线观看| 80岁老熟妇乱子伦牲交| 欧美日韩亚洲国产一区二区在线观看 | 老司机影院毛片| 亚洲熟女精品中文字幕| 亚洲色图 男人天堂 中文字幕| 国产av又大| 不卡av一区二区三区| 女人久久www免费人成看片| 欧美另类亚洲清纯唯美| 国产无遮挡羞羞视频在线观看| 成年美女黄网站色视频大全免费| 国产成人啪精品午夜网站| 亚洲国产成人一精品久久久| 日韩欧美免费精品| 青春草亚洲视频在线观看| 嫩草影视91久久| 1024视频免费在线观看| 少妇的丰满在线观看| 9191精品国产免费久久| 80岁老熟妇乱子伦牲交| 别揉我奶头~嗯~啊~动态视频 | 久久久欧美国产精品| 日韩制服骚丝袜av| 老司机亚洲免费影院| av欧美777| 午夜精品久久久久久毛片777| 亚洲精品国产av成人精品| 美女午夜性视频免费| 成年人黄色毛片网站| 一本综合久久免费| 在线永久观看黄色视频| 最黄视频免费看| 亚洲专区字幕在线| 天天躁夜夜躁狠狠躁躁| 人人妻人人澡人人爽人人夜夜| 免费不卡黄色视频| 久久国产亚洲av麻豆专区| 国产精品香港三级国产av潘金莲| 男人操女人黄网站| 无限看片的www在线观看| 老司机影院毛片| 午夜福利在线观看吧| 国产成人免费无遮挡视频| 亚洲第一av免费看| 欧美成狂野欧美在线观看| 一级毛片女人18水好多| 精品国产乱子伦一区二区三区 | 欧美日韩中文字幕国产精品一区二区三区 | 91精品三级在线观看| 麻豆av在线久日| 高清av免费在线| 精品人妻1区二区| 99久久99久久久精品蜜桃| 女人被躁到高潮嗷嗷叫费观| 国产主播在线观看一区二区| 两人在一起打扑克的视频| 国产99久久九九免费精品| 国产日韩欧美视频二区| 人妻 亚洲 视频| 亚洲一码二码三码区别大吗| 最新在线观看一区二区三区| 国产成人a∨麻豆精品| 成人av一区二区三区在线看 | 欧美黄色淫秽网站| 丝袜美足系列| 纯流量卡能插随身wifi吗| 欧美老熟妇乱子伦牲交| 麻豆av在线久日| 欧美少妇被猛烈插入视频| 色婷婷久久久亚洲欧美| 国产在线观看jvid| 免费女性裸体啪啪无遮挡网站| e午夜精品久久久久久久| 99国产精品99久久久久| 亚洲成av片中文字幕在线观看| 90打野战视频偷拍视频| 麻豆国产av国片精品| 人人澡人人妻人| 秋霞在线观看毛片| 国产精品久久久人人做人人爽| 久久久久久久久久久久大奶| 精品欧美一区二区三区在线| 午夜91福利影院| 最近中文字幕2019免费版| 一区在线观看完整版| 久久久久久久精品精品| 丝袜在线中文字幕| 99re6热这里在线精品视频| 欧美成狂野欧美在线观看| 国产成人av激情在线播放| 久久久久久久久久久久大奶| 国产片内射在线| 丰满饥渴人妻一区二区三| 久久av网站| 亚洲熟女精品中文字幕| 午夜免费鲁丝| 亚洲精品国产区一区二| 日本精品一区二区三区蜜桃| 欧美日韩视频精品一区| 脱女人内裤的视频| 亚洲国产成人一精品久久久| 欧美激情极品国产一区二区三区| 9色porny在线观看| 我要看黄色一级片免费的| 搡老岳熟女国产| 黄频高清免费视频| 色综合欧美亚洲国产小说| 欧美国产精品va在线观看不卡| 国产一区有黄有色的免费视频| 欧美xxⅹ黑人| 国产精品免费大片| 国产精品免费视频内射| 亚洲精品一二三| 欧美成狂野欧美在线观看| 久久99热这里只频精品6学生| 免费在线观看影片大全网站| 欧美激情极品国产一区二区三区| 视频区图区小说| 天天躁日日躁夜夜躁夜夜| 午夜影院在线不卡| 精品福利永久在线观看| 欧美另类一区| 亚洲精品在线美女| 国产在线视频一区二区| 午夜老司机福利片| 一本大道久久a久久精品| 男女无遮挡免费网站观看| 亚洲一区二区三区欧美精品| 9色porny在线观看| 在线亚洲精品国产二区图片欧美| 国产一区二区激情短视频 | 国产无遮挡羞羞视频在线观看| 亚洲国产日韩一区二区| 日韩制服骚丝袜av| 在线十欧美十亚洲十日本专区| 一区二区日韩欧美中文字幕| 亚洲五月婷婷丁香| 视频区欧美日本亚洲| 日本一区二区免费在线视频| 国产成+人综合+亚洲专区| 国产真人三级小视频在线观看| 亚洲激情五月婷婷啪啪| 亚洲欧美日韩高清在线视频 | 少妇粗大呻吟视频| 咕卡用的链子| 手机成人av网站| 人人妻人人爽人人添夜夜欢视频| 在线观看免费午夜福利视频| 99九九在线精品视频| 日韩电影二区| 十八禁人妻一区二区| 久久 成人 亚洲| 亚洲伊人久久精品综合| 久久久久久久精品精品| 男女午夜视频在线观看| 亚洲黑人精品在线| 黄色a级毛片大全视频| 午夜免费成人在线视频| 老鸭窝网址在线观看| 欧美亚洲日本最大视频资源| 亚洲国产精品999| 国产一区二区在线观看av| 成年人午夜在线观看视频| 国产伦理片在线播放av一区| 99热国产这里只有精品6| 欧美精品啪啪一区二区三区 | 亚洲国产精品一区二区三区在线| 国产高清国产精品国产三级| 在线亚洲精品国产二区图片欧美| 男女之事视频高清在线观看| 中国国产av一级| av天堂在线播放| 亚洲精品自拍成人| 黄色毛片三级朝国网站| 欧美在线黄色| 久9热在线精品视频| 夜夜骑夜夜射夜夜干| 黄片小视频在线播放| 两人在一起打扑克的视频| 欧美97在线视频| 久久精品国产亚洲av高清一级| 正在播放国产对白刺激| 黑人巨大精品欧美一区二区蜜桃| 国产免费福利视频在线观看| 久久人人爽av亚洲精品天堂| 人人妻,人人澡人人爽秒播| 亚洲一区二区三区欧美精品| 精品少妇内射三级| 两个人免费观看高清视频| 99久久国产精品久久久| 国产精品1区2区在线观看. | 精品国产乱子伦一区二区三区 | 国产一区二区在线观看av| 一二三四在线观看免费中文在| 亚洲国产日韩一区二区| 久久久国产一区二区| 在线亚洲精品国产二区图片欧美| 国产成人精品无人区| 日日爽夜夜爽网站| 久久精品国产亚洲av香蕉五月 | 18禁观看日本| 中亚洲国语对白在线视频| 亚洲av片天天在线观看| 91麻豆精品激情在线观看国产 | 国产精品成人在线| 国产成人精品久久二区二区免费| 欧美变态另类bdsm刘玥| 免费不卡黄色视频| 亚洲成人免费电影在线观看| 午夜福利在线观看吧| netflix在线观看网站| 亚洲视频免费观看视频| 色婷婷av一区二区三区视频| 99热全是精品| 男女之事视频高清在线观看| 亚洲精品久久午夜乱码| 男人添女人高潮全过程视频| 女警被强在线播放| 国产视频一区二区在线看| 少妇裸体淫交视频免费看高清 | 国产欧美亚洲国产| 欧美国产精品一级二级三级| 韩国精品一区二区三区| 汤姆久久久久久久影院中文字幕| 一本大道久久a久久精品| 成年人免费黄色播放视频| 丰满人妻熟妇乱又伦精品不卡| 精品人妻一区二区三区麻豆| 制服诱惑二区| 欧美日韩精品网址| 99国产综合亚洲精品| 脱女人内裤的视频| 亚洲视频免费观看视频| 亚洲国产欧美在线一区| 国产深夜福利视频在线观看| 久久性视频一级片| 777久久人妻少妇嫩草av网站| 午夜福利免费观看在线| 久久久久久久久免费视频了| 国产精品二区激情视频| 久久久国产一区二区| av在线app专区| 日韩有码中文字幕| 久热爱精品视频在线9| 中文字幕人妻丝袜制服| 亚洲 欧美一区二区三区| 国产不卡av网站在线观看| 老司机福利观看| 国产日韩欧美亚洲二区| 精品国产乱码久久久久久小说| 秋霞在线观看毛片| 大型av网站在线播放| 波多野结衣av一区二区av| 波多野结衣一区麻豆| 亚洲一区中文字幕在线| 天天添夜夜摸| 一区二区三区四区激情视频| 欧美激情 高清一区二区三区| 欧美另类亚洲清纯唯美| 少妇精品久久久久久久| 亚洲人成电影免费在线| 999久久久精品免费观看国产| 1024视频免费在线观看| 色老头精品视频在线观看| 免费观看人在逋| 午夜免费鲁丝| 免费在线观看完整版高清| 亚洲avbb在线观看| 国产一区有黄有色的免费视频| 菩萨蛮人人尽说江南好唐韦庄| 中文字幕制服av| 男人操女人黄网站| 国产亚洲一区二区精品| 这个男人来自地球电影免费观看| 手机成人av网站| 亚洲一码二码三码区别大吗| 秋霞在线观看毛片| 男女免费视频国产| 国产精品久久久久成人av| 99精国产麻豆久久婷婷| 热re99久久国产66热| 女人久久www免费人成看片| 亚洲 欧美一区二区三区| 精品少妇黑人巨大在线播放| 欧美少妇被猛烈插入视频| 午夜福利乱码中文字幕| 久久久久国产一级毛片高清牌| 精品人妻一区二区三区麻豆| 久久这里只有精品19| 韩国精品一区二区三区| 精品少妇久久久久久888优播| 亚洲av美国av| 欧美久久黑人一区二区| 国产伦理片在线播放av一区| av天堂久久9| 黄片小视频在线播放| 久久国产精品人妻蜜桃| 亚洲第一青青草原| 精品第一国产精品| 久久av网站| 人人妻人人添人人爽欧美一区卜| 人人妻人人爽人人添夜夜欢视频| 成人18禁高潮啪啪吃奶动态图| 欧美午夜高清在线| 成人18禁高潮啪啪吃奶动态图| 久久中文看片网| 久久综合国产亚洲精品| 一本一本久久a久久精品综合妖精| 国产日韩一区二区三区精品不卡| 黄片小视频在线播放| 91老司机精品| 欧美日本中文国产一区发布| 国产av国产精品国产| 成人免费观看视频高清| 精品卡一卡二卡四卡免费|