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

    單脈沖搜索在CRAFTS巡天數(shù)據(jù)中的研究與應用

    2022-07-18 00:56:50于徐紅陳宗浩劉志杰游善平
    天文研究與技術 2022年4期
    關鍵詞:單脈沖脈沖星峰度

    張 彬,梁 楠,2,于徐紅,陳宗浩,劉志杰,游善平*

    (1. 貴州師范大學貴州省信息與計算科學重點實驗室,貴州 貴陽 550001;2. 中國科學院云南天文臺,云南 昆明 650216)

    脈沖星是一種快速旋轉、高磁化的中子星,是理解超密度物質的關鍵,可以用于各種基礎物理實驗[1]。目前,脈沖星搜索方法主要分為周期性搜索和單脈沖搜索兩大類。其中,周期性搜索通過應用快速傅里葉變換(Fast Fourier Transform, FFT)將時間序列轉化到頻域,以識別周期性信號,然后在確定的周期內對原始時間序列數(shù)據(jù)進行折疊,以提高周期性信號的信噪比[2]。而單脈沖搜索尋找強的、非周期的脈沖,并未使用快速傅里葉變換和折疊。單脈沖搜索非常適合發(fā)現(xiàn)周期性搜索中無法發(fā)現(xiàn)的孤立爆發(fā)。研究人員應用單脈沖搜索發(fā)現(xiàn)了旋轉射電暫現(xiàn)源(Rotating Radio Transients, RRAT)[3]和快速射電暴的(Fast Radio Bursts, FRB)[4]。

    2003年,文[5]首次提出了一個理論框架,運用單脈沖搜索方法探測射電天文觀測數(shù)據(jù)中的脈沖星信號。2009年,文[6]把單脈沖搜索應用于阿雷西博L波段饋源陣列脈沖星巡天(Pulsar Arecibo L-band Feed Array, PALFA)七波束觀測數(shù)據(jù),發(fā)現(xiàn)了7顆脈沖星。2010年,文[7]使用單脈沖搜索方法對帕克斯多波束脈沖星巡天(Parkes Multibeam Pulsar Survey, PMPS)數(shù)據(jù)進行重新分析,發(fā)現(xiàn)10顆旋轉射電暫現(xiàn)源。2015年,文[8]在單脈沖理論的基礎上,設計針對旋轉射電暫現(xiàn)源的搜索工具RRATtrap,應用在綠岸射電望遠鏡(Green Bank Telescope, GBT)350 MHz觀測數(shù)據(jù),最終探測到18顆旋轉射電暫現(xiàn)源。2016年,文[9]首次將單脈沖搜索與機器學習相結合,提出遞歸峰值識別算法(Recursive Algorithm for Peak IDentification, RAPID),實現(xiàn)了自動化篩選脈沖星候選體,并探測到6顆脈沖星。

    2016年9月,位于貴州喀斯特山區(qū)的500 m口徑球面射電望遠鏡進入科學運行階段,望遠鏡先后配置了超寬帶接收機和L波段19波束接收機[10]。為充分利用FAST優(yōu)異的探測能力,2018年,文[11]設計并實施了FAST多科學目標同時掃描巡天,即同時使用多個數(shù)字終端采集脈沖星、中子氫、分子譜線、旋轉射電暫現(xiàn)源、快速射電暴等多個科學目標的觀測數(shù)據(jù)。鑒于單脈沖搜索方法在各大射電望遠鏡的成功應用,本文對CRAFTS超寬帶數(shù)據(jù)文件的單脈沖搜索結果進行研究,發(fā)現(xiàn)CRAFTS單脈沖搜索篩選得到脈沖星候選體存在數(shù)以萬計的假陽性樣本。因此,我們根據(jù)脈沖星信號在色散量信噪比曲線大致是高斯曲線以及信噪比數(shù)據(jù)分布是高斯分布的假設,構造3個顯著區(qū)分脈沖星信號和干擾的特征,旨在緩解出現(xiàn)大量假陽性樣本的問題。試驗結果表明,添加特征值判斷的單脈沖篩選工具與傳統(tǒng)方法相比,假陽性樣本數(shù)量減少了20%。

    1 CRAFTS超寬帶數(shù)據(jù)單脈沖搜索實驗

    在射電天文觀測數(shù)據(jù)中發(fā)現(xiàn)脈沖星通常分為4個階段:收集、消色散、周期性搜索或單脈沖搜索和人工檢查[1]。第1階段,原始數(shù)據(jù)以電壓時間序列的形式由射電望遠鏡收集;第2階段,脈沖輻射通過星際介質(Interstellar Medium, ISM)的色散效應,導致較低頻率的脈沖比較高頻率晚到達[12],消色散可以去除這些與頻率有關的延遲效應的影響;第3階段,使用周期性搜索或單脈沖搜索找出脈沖星候選體;第4階段,對判斷為脈沖星候選體的數(shù)據(jù)進行人工檢查。

    我們先使用PRESTO(Pulsar Exploration and Search Toolkit)脈沖星搜索工具處理2017年8月至2018年5月CRAFTS的超寬帶數(shù)據(jù),處理流程包括去干擾、消色散和單脈沖搜索等步驟;然后,應用單脈沖篩選工具RRATtrap從PRESTO的Single_pulse_search.py輸出文件中篩選出脈沖星候選體。

    單脈沖篩選工具RRATtrap根據(jù)脈沖星信號的特性區(qū)分脈沖星信號和干擾。(1)脈沖星信號出現(xiàn)在一定色散量(Dispersion Measure, DM)范圍內,在最佳色散量時檢測到峰值信噪比,而在該色散量上下信噪比下降,這是由于脈沖以不準確的色散量消色散時導致脈沖展寬造成的。(2)由于信號在最佳色散量檢測到峰值信噪比,預計來自地面的信號(即射頻干擾)在色散量為0 pc·cm-3時達到峰值。同時射頻干擾(Radio Frequency Interference, RFI)不受色散效應的影響,因此,脈沖星信號出現(xiàn)在一個非常大的色散量范圍內[8]。

    圖1是單脈沖篩選工具RRATtrap探測到已知脈沖星PSR B0540+23的診斷圖。左側子圖縱坐標為試驗色散量范圍,橫坐標為觀測時間,CRAFTS超寬帶數(shù)據(jù)每個觀測文件為52 s,右側子圖描述試驗色散量對應的信噪比大小。從圖1可以看到呈紡錘狀的單脈沖事件組(在相鄰的色散量與時間窗口內所有信噪比大于5的事件)在最佳色散量78 pc·cm-3時取得峰值信噪比68,而在該色散量上下,信噪比下降。其中,強的脈沖星信號標記為Excellent;弱的脈沖星信號表標記為Very good;不規(guī)則的脈沖星信號標記為Good或Ok;射頻干擾標記為RFI;宇宙噪聲標記為Noise。在圖2中,我們看到在較寬的色散量范圍上極強的射頻干擾的實例。此外,在t=20~40 s,極強的射頻干擾錯誤標記為脈沖星信號。

    圖1 RRATtrap探測到脈沖星PSR B0540+23的診斷圖

    圖2 RRATtrap篩選出假陽性實例的診斷圖

    本文試驗選取2017年8月至2018年5月共計約32萬個超寬帶CRAFTS數(shù)據(jù)文件[1]進行單脈沖搜索試驗。由于處理的數(shù)據(jù)規(guī)模較大,對所有數(shù)據(jù)進行全面的人工檢查工作量非常大。因此,我們先使用RRATtrap初步篩選脈沖星候選體。然后,我們對數(shù)據(jù)文件中標記為Excellent,Very good,Good和Ok共計約5萬顆脈沖星候選體的診斷圖進行手工檢查,發(fā)現(xiàn)僅有772個診斷圖真正具有脈沖星信號,對應101顆已知脈沖星。通過

    (1)

    計算的假陽性率(False Positive Rate, FPR)達到98.5%,其中,F(xiàn)P是沒有脈沖星信號的數(shù)據(jù)文件標記為候選脈沖星的數(shù)量;TP是含有脈沖星信號的數(shù)據(jù)文件標記為候選脈沖星的文件數(shù)量。每100個脈沖星候選體中,僅有1~2個包含脈沖星信號(已知或新發(fā)現(xiàn)的)。由于CRAFTS超寬帶數(shù)據(jù)文件中包含的具體脈沖星數(shù)量未知,本文未對假陰性情況進行分析。

    2 區(qū)分脈沖星信號和干擾的特征

    CRAFTS超寬帶數(shù)據(jù)單脈沖搜索結果表明,脈沖星候選體僅有1.5%真正具有脈沖星信號,這主要是未能成功區(qū)分脈沖星信號和干擾造成的,如圖2。本節(jié)首先探討不準確的色散量對信噪比的影響。然后,構造3個顯著區(qū)分脈沖星信號和干擾的特征,并在探測到的101顆脈沖星樣本中選取同時具有脈沖星信號、射頻干擾和宇宙噪聲的79顆脈沖星作為樣本,計算它們在3個特征的取值情況。最后,根據(jù)脈沖星信號、射頻干擾和宇宙噪聲在3個特征取值分布的差異,提出合理閾值應用在單脈沖篩選工具,對脈沖星候選體進行進一步篩選。

    2.1 色散量與信噪比曲線

    文[5]探討了試驗與真實色散量之間的偏差對信噪比的影響。通常,對一個特定色散量的時間序列進行多次下采樣并重新搜索,當有效的采樣時間最接近脈沖寬度時,得到的信噪比最高。在不同的色散量信道中,隨著試驗色散量與真實色散量的偏差越大,信噪比越小,如圖1。測量信噪比S(δDM)與真實信噪比S的比值與色散量的偏差δDM滿足[5]

    (2)

    這里,

    (3)

    其中,erf為誤差函數(shù);δDM為試驗色散量與真實色散量的偏差;Δv為總帶寬;v為中心頻率;W為脈沖寬度,單位ms。

    利用(2)式和(3)式計算天體物理脈沖的信噪比和寬度,我們可以計算預期的色散量偏差δDM。圖3展示了在中心頻率546 MHz,時間分辨率為100 μs和匹配濾波使用30的PSR B2000+40一組單脈沖事件預期和擬合信噪比下降曲線。在這兩種情況下,信噪比峰值處的色散量是真實的色散量。我們計算預期的信噪比下降,利用峰值信噪比和匹配濾波得到脈沖寬度,和實際脈沖寬度很接近。首先,我們使用非線性最小二乘法(Non-linear Least Squares, NLS)對原始數(shù)據(jù)進行回歸,得到擬合的峰值信噪比和脈沖寬度。然后,根據(jù)擬合結果,代入(2)式得到擬合的色散與信噪比曲線。由圖3可以看出,預期和擬合的色散量與信噪比曲線大致是高斯曲線。使用觀測到的峰值信噪比和匹配濾波得到脈沖寬度,預期的信噪比下降與擬合值在峰值右側非常接近,但在左側出現(xiàn)偏差。這是匹配濾波得到的脈沖寬度與實際脈沖寬度的偏差造成的。

    2.2 3個顯著特征

    文[13]指出,脈沖星的輪廓可以簡化為高斯型(對大多數(shù)脈沖星來說是一個合理的近似)。按照文[5]提出的單脈沖搜索理論,天體物理脈沖的色散量和信噪比曲線大致是高斯曲線,如圖3。由于高斯曲線通常是對稱的,文[14]提出了兩種對稱特征(SIDM和SIS/N)分別表征單脈沖事件組的色散量和信噪比的對稱性。SIDM的計算公式為

    (4)

    其中,DMpeak為信噪比峰值對應的色散量;DMmax為單脈沖事件組最大的色散量;DMmin為單脈沖事件組最小的色散量。SIS/N的計算公式為

    (5)

    其中,∑S/Nleft為信噪比峰值左側所有單脈沖事件的信噪比之和;∑S/Nright為信噪比峰值右側所有單脈沖事件的信噪比之和。

    由(4)式和(5)式定義的SIDM和SIS/N的取值范圍在0~1之間,這兩個特征值越高,色散量與信噪比曲線越對稱,相反不遵循(2)式描述規(guī)律的射頻干擾,通常在色散量與信噪比空間是單調遞減(或遞增)。所以,它們通常具有接近于0的對稱值,可以與脈沖星信號進行區(qū)分。圖4和圖5分別展示了脈沖星樣本SIDM和SIS/N的對稱值,并對脈沖星信號、射頻干擾和噪聲進行了對比。從圖4和圖5可以看出,脈沖星信號通常具有較大的對稱值,而射頻干擾和噪聲的對稱值相對較小,這表明脈沖星信號的色散量和信噪比曲線比干擾更對稱。此外,圖4和圖5標記一些對稱值比較大的噪聲,這些通常是高斯噪聲。

    圖4 脈沖星樣本SIDM對稱值

    圖5 脈沖星樣本SIS/N對稱值

    文[5]預測在沒有任何宇宙噪聲和射頻干擾的情況下,信噪比的數(shù)據(jù)分布是高斯分布。峰度是描述總體數(shù)據(jù)分布與高斯分布陡緩程度的統(tǒng)計量。為了判斷脈沖星樣本信噪比的數(shù)據(jù)分布是否符合高斯分布,我們使用峰度值表征單脈沖事件組內信噪比的數(shù)據(jù)分布與高斯分布的符合程度。峰度值的計算公式為

    (6)

    其中,S/N為單脈沖事件的信噪比;μ為單脈沖事件組內平均信噪比;σ為單脈沖事件組內信噪比的方差。

    當(6)式定義的峰度值是3時,信噪比的數(shù)據(jù)分布服從高斯分布,隨著峰度值與3的差距越大,其分布形態(tài)的陡緩程度與高斯分布的差異程度越大。圖6用箱線圖展示了脈沖星樣本的峰度值分布情況,并對脈沖星信號、射頻干擾和噪聲進行了對比。從圖6可以看出,脈沖星信號峰度值集中在3附近,而射頻干擾和噪聲的峰度值分散,且與3差距較大,表明脈沖星信號信噪比分布近似服從高斯分布,而大部分干擾與噪聲不具有此規(guī)律。

    圖4~圖6可以明顯看出脈沖星信號、射頻干擾和宇宙噪聲在3個特征分布的差異。因此,我們根據(jù)脈沖星信號和干擾特征值的分布區(qū)間(表1),選擇合適的特征閾值,剔除不滿足閾值的脈沖星候">選體,從而達到降低假陽性率的目的(比如設置脈沖星信號特征SIDM最小值0.1作為閾值,對低于閾值的候選體視為干擾,那么,特征值在0.02~0.1區(qū)間的射頻干擾將不會錯誤標記為脈沖星候選體)。同時,我們選取不同閾值進行試驗,并統(tǒng)計其對應的假陽性率和遺漏脈沖星信號的結果。由圖7各個特征在不同閾值假陽性及遺漏脈沖星的情況可以看出,隨著閾值變大或者閾值區(qū)間縮小,假陽性率不斷減小,與此同時,遺漏發(fā)現(xiàn)的脈沖星數(shù)量在不斷增加。脈沖星搜尋的首要前提是保證觀測數(shù)據(jù)中所有脈沖星信號不遺漏。因此,我們選取表1中脈沖星信號特征SIDM最小值0.1和SIS/N最小值0.29作為閾值,特征峰度-1.54~27.11作為閾值區(qū)間,對低于閾值或不在閾值區(qū)間的脈沖星候選體視為射頻干擾或宇宙噪聲。

    圖6 脈沖星樣本峰度值

    表1 3個特征在3種類別的分布區(qū)間

    圖7 各個特征在不同閾值假陽性和遺漏脈沖星的情況,其中紅色豎線代表所選擇的閾值

    此外,為了證實所述特征是否在其他觀測數(shù)據(jù)有效,我們對帕克斯多波束數(shù)據(jù)進行了試驗。根據(jù)文[15]公開的帕克斯單脈沖數(shù)據(jù)庫,我們構建了一個PMPS(Parkes Multibeam Pulsar Survey)數(shù)據(jù)庫,其中包含帕克斯望遠鏡發(fā)現(xiàn)的部分脈沖星信號。我們計算得到PMPS數(shù)據(jù)集中脈沖星信號在3個特征值的分布情況,并與CRAFTS進行對比,如表1。由表1中CRAFTS及PMPS的脈沖星信號在3個特征值分布情況可以發(fā)現(xiàn),它們的特征值分布十分接近,盡管它們在消色散網(wǎng)格、周圍射頻干擾的環(huán)境等存在顯著差異。因此,本文所述的特征閾值同樣能夠區(qū)分帕克斯數(shù)據(jù)庫中脈沖星信號和干擾。

    我們將上述閾值和閾值區(qū)間應用于單脈沖篩選工具,對脈沖星候選體進行重新篩選,試驗結果如表2。由表2可以看出,添加單個特征閾值判斷的單脈沖篩選工具能在一定程度上降低假陽性率,并且3個特征閾值組合判斷假陽性率從98.5%降低到78.4%,表明添加特征閾值判斷的單脈沖篩選工具提高了脈沖星搜索效率。

    表2 添加特征后假陽性率結果

    3 總結與展望

    本文使用PRESTO的單脈沖搜索方法對CRAFTS超寬帶數(shù)據(jù)文件進行的試驗表明,現(xiàn)有的單脈沖篩選工具難以區(qū)分真實信號與噪聲或射頻干擾信號,造成數(shù)以萬計的假陽性樣本出現(xiàn),顯著增加人工篩選單脈沖候選體的時間開銷以及候選數(shù)據(jù)存儲壓力。因此,我們提出3個顯著區(qū)分脈沖星信號和干擾的特征,并選取同時具有脈沖星信號、射頻干擾和宇宙噪聲的79顆脈沖星樣本,計算它們在3個特征取值情況。然后,根據(jù)脈沖星信號、射頻干擾和宇宙噪聲在3個特征取值分布的差異,提出合理閾值應用于單脈沖篩選工具,用于對脈沖星候選體進行嚴格的判斷。最后,使用添加特征閾值判斷的篩選工具對脈沖星候選體進行重新試驗。結果表明,在保證所有脈沖星信號不遺漏的情況下,假陽性率從98.5%降低到78.4%。因此,本文所述3個特征具有實用性和有效性,有助于單脈沖搜索在CRAFTS巡天數(shù)據(jù)的應用。

    致謝:本文在500 m口徑球面射電望遠鏡(FAST)數(shù)據(jù)基礎上完成。FAST是由中國科學院國家天文臺運行和管理的國家大科學裝置。感謝中國科學院天文大科學研究中心FAST重大成果培育項目對本文工作的資助。

    猜你喜歡
    單脈沖脈沖星峰度
    “中國天眼”已發(fā)現(xiàn)740余顆新脈沖星
    軍事文摘(2023年12期)2023-06-12 07:51:00
    發(fā)現(xiàn)脈沖星的女天文學家——貝爾
    科學(2022年4期)2022-10-25 02:43:42
    擴散峰度成像技術檢測急性期癲癇大鼠模型的成像改變
    FDA對比幅法單脈沖測向的角度欺騙
    磁共振擴散峰度成像在肝臟病變中的研究進展
    一類帶有慢變參數(shù)的sine-Gordon方程的單脈沖異宿軌道
    基于自動反相校正和峰度值比較的探地雷達回波信號去噪方法
    雷達學報(2018年3期)2018-07-18 02:41:16
    基于虛擬觀測值的X射線單脈沖星星光組合導航
    一種帶寬展寬的毫米波波導縫隙陣列單脈沖天線
    制導與引信(2016年3期)2016-03-20 16:02:00
    磁共振擴散峰度成像MK值、FA值在鑒別高級別膠質瘤與轉移瘤的價值分析
    亚洲片人在线观看| 国产久久久一区二区三区| 韩国av一区二区三区四区| 亚洲中文字幕日韩| 久久中文看片网| 欧美在线一区亚洲| 久久久久国产一级毛片高清牌| 一区二区三区国产精品乱码| 精品电影一区二区在线| 国产av不卡久久| 午夜免费激情av| 99热这里只有精品一区 | 久久久久国产一级毛片高清牌| 一区二区三区国产精品乱码| 日韩精品青青久久久久久| 免费高清视频大片| 亚洲五月婷婷丁香| 在线播放国产精品三级| 日本熟妇午夜| 成人三级黄色视频| 久久中文字幕一级| 亚洲第一欧美日韩一区二区三区| 两性午夜刺激爽爽歪歪视频在线观看 | 黄频高清免费视频| 成人国语在线视频| 亚洲欧美日韩东京热| 国产男靠女视频免费网站| 亚洲精品国产精品久久久不卡| 亚洲精品久久国产高清桃花| svipshipincom国产片| 亚洲欧美日韩东京热| 亚洲欧美日韩无卡精品| 一夜夜www| 99热6这里只有精品| 可以在线观看毛片的网站| 精品国产超薄肉色丝袜足j| 午夜免费观看网址| 国产av又大| 日本一本二区三区精品| 免费在线观看影片大全网站| 国产黄a三级三级三级人| 欧美成狂野欧美在线观看| 制服人妻中文乱码| 日韩中文字幕欧美一区二区| 亚洲熟妇中文字幕五十中出| 国产熟女午夜一区二区三区| 日日摸夜夜添夜夜添小说| 色哟哟哟哟哟哟| 亚洲成人国产一区在线观看| 亚洲欧美日韩高清专用| 99国产精品99久久久久| 女人高潮潮喷娇喘18禁视频| 搡老妇女老女人老熟妇| 亚洲国产中文字幕在线视频| 熟女电影av网| 亚洲av成人av| 中文字幕精品亚洲无线码一区| 国产av麻豆久久久久久久| 老司机靠b影院| 日本精品一区二区三区蜜桃| 可以在线观看的亚洲视频| 怎么达到女性高潮| 久久精品人妻少妇| 午夜亚洲福利在线播放| 亚洲国产精品久久男人天堂| 久久婷婷成人综合色麻豆| 99精品欧美一区二区三区四区| 久久久久九九精品影院| 制服人妻中文乱码| 亚洲全国av大片| 亚洲精品国产精品久久久不卡| 可以在线观看的亚洲视频| 精品高清国产在线一区| 人人妻,人人澡人人爽秒播| 琪琪午夜伦伦电影理论片6080| 日本一区二区免费在线视频| 搞女人的毛片| 久久婷婷人人爽人人干人人爱| 男女视频在线观看网站免费 | 国产成年人精品一区二区| 桃红色精品国产亚洲av| 亚洲中文日韩欧美视频| 亚洲国产精品999在线| 国产亚洲精品久久久久5区| 老司机午夜福利在线观看视频| 又黄又粗又硬又大视频| 婷婷六月久久综合丁香| 高清毛片免费观看视频网站| 99国产精品99久久久久| 国产精品久久久av美女十八| 日本 欧美在线| 亚洲第一电影网av| 国产高清视频在线播放一区| 久久中文字幕一级| 天堂影院成人在线观看| 久久久久久人人人人人| 免费看美女性在线毛片视频| 国产片内射在线| 夜夜看夜夜爽夜夜摸| 国产成年人精品一区二区| 亚洲一卡2卡3卡4卡5卡精品中文| 国产1区2区3区精品| 此物有八面人人有两片| 人成视频在线观看免费观看| 怎么达到女性高潮| 成人三级黄色视频| 青草久久国产| 亚洲一区二区三区不卡视频| 日韩精品免费视频一区二区三区| 真人一进一出gif抽搐免费| 99精品在免费线老司机午夜| 美女免费视频网站| 给我免费播放毛片高清在线观看| 91在线观看av| 国产蜜桃级精品一区二区三区| 国产成人影院久久av| 国产成人av激情在线播放| 亚洲专区中文字幕在线| 国产成人精品无人区| 国产久久久一区二区三区| 亚洲 欧美 日韩 在线 免费| 最近视频中文字幕2019在线8| x7x7x7水蜜桃| 国产精品乱码一区二三区的特点| 18禁国产床啪视频网站| 99热6这里只有精品| 欧美 亚洲 国产 日韩一| 男女下面进入的视频免费午夜| 丝袜人妻中文字幕| 99久久99久久久精品蜜桃| 天天躁狠狠躁夜夜躁狠狠躁| 俺也久久电影网| 窝窝影院91人妻| av免费在线观看网站| 国产亚洲精品第一综合不卡| a级毛片在线看网站| 日本一区二区免费在线视频| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美日韩一级在线毛片| 色综合婷婷激情| 一级a爱片免费观看的视频| 母亲3免费完整高清在线观看| 制服人妻中文乱码| 亚洲美女黄片视频| 亚洲aⅴ乱码一区二区在线播放 | 男女之事视频高清在线观看| 黄色丝袜av网址大全| 一进一出抽搐gif免费好疼| 久久久精品欧美日韩精品| 狂野欧美激情性xxxx| 国产精品99久久99久久久不卡| 99精品欧美一区二区三区四区| 午夜激情福利司机影院| 手机成人av网站| 美女高潮喷水抽搐中文字幕| 精品午夜福利视频在线观看一区| 一区二区三区激情视频| 91av网站免费观看| xxx96com| 成人三级黄色视频| 亚洲五月婷婷丁香| 日本熟妇午夜| 亚洲精品色激情综合| 国产成人一区二区三区免费视频网站| 老熟妇乱子伦视频在线观看| 国产一区二区激情短视频| 性欧美人与动物交配| 国产精品98久久久久久宅男小说| 国产精品爽爽va在线观看网站| 精品久久蜜臀av无| 12—13女人毛片做爰片一| 免费观看人在逋| 国产av一区在线观看免费| 国产精品综合久久久久久久免费| 亚洲av第一区精品v没综合| 中文在线观看免费www的网站 | 色av中文字幕| 日韩免费av在线播放| 一个人免费在线观看电影 | 一区二区三区激情视频| 亚洲精品中文字幕一二三四区| 午夜久久久久精精品| 嫩草影视91久久| 久久久久久久精品吃奶| 国产探花在线观看一区二区| 婷婷六月久久综合丁香| 可以在线观看的亚洲视频| 成人高潮视频无遮挡免费网站| 18禁国产床啪视频网站| 精品久久蜜臀av无| 欧美+亚洲+日韩+国产| 亚洲一码二码三码区别大吗| 麻豆一二三区av精品| 巨乳人妻的诱惑在线观看| 国产精华一区二区三区| 欧美精品啪啪一区二区三区| 国产亚洲欧美在线一区二区| 中文在线观看免费www的网站 | 色尼玛亚洲综合影院| 日韩欧美 国产精品| 两个人视频免费观看高清| 久久久久久久久免费视频了| 亚洲国产精品sss在线观看| www.www免费av| 人人妻人人看人人澡| 日韩精品免费视频一区二区三区| xxxwww97欧美| 老汉色∧v一级毛片| 性色av乱码一区二区三区2| 久久国产精品人妻蜜桃| 大型黄色视频在线免费观看| 欧美zozozo另类| 亚洲 欧美一区二区三区| 一a级毛片在线观看| 麻豆成人午夜福利视频| 小说图片视频综合网站| 午夜亚洲福利在线播放| 国产成人精品久久二区二区免费| 黄片小视频在线播放| 国产亚洲欧美98| 亚洲成人免费电影在线观看| 亚洲av第一区精品v没综合| 啦啦啦观看免费观看视频高清| 日韩欧美精品v在线| 国产精华一区二区三区| 国产精品1区2区在线观看.| 亚洲中文av在线| 亚洲国产中文字幕在线视频| 国产精品久久久久久亚洲av鲁大| 午夜福利在线观看吧| 国产亚洲精品av在线| 91麻豆精品激情在线观看国产| 国产日本99.免费观看| 小说图片视频综合网站| 国产精品av久久久久免费| 男女那种视频在线观看| 国内揄拍国产精品人妻在线| 久久久国产成人精品二区| 亚洲国产高清在线一区二区三| 久久中文字幕一级| 精品久久久久久久久久免费视频| 国产不卡一卡二| 可以在线观看的亚洲视频| 黄色成人免费大全| 亚洲七黄色美女视频| 国产激情欧美一区二区| 亚洲国产高清在线一区二区三| 日韩欧美一区二区三区在线观看| 最新在线观看一区二区三区| 两个人的视频大全免费| 18禁裸乳无遮挡免费网站照片| 亚洲av片天天在线观看| 天堂√8在线中文| 欧美另类亚洲清纯唯美| 国产成人系列免费观看| e午夜精品久久久久久久| 不卡一级毛片| a级毛片在线看网站| 国产亚洲精品一区二区www| 男女床上黄色一级片免费看| 国产人伦9x9x在线观看| 日韩精品免费视频一区二区三区| 精品一区二区三区四区五区乱码| 黄色 视频免费看| avwww免费| 非洲黑人性xxxx精品又粗又长| 国产精品98久久久久久宅男小说| 亚洲一码二码三码区别大吗| 老司机深夜福利视频在线观看| 天堂动漫精品| 欧美大码av| 无遮挡黄片免费观看| 麻豆国产av国片精品| 国产精品一区二区精品视频观看| 亚洲一码二码三码区别大吗| netflix在线观看网站| 美女高潮喷水抽搐中文字幕| 美女扒开内裤让男人捅视频| 神马国产精品三级电影在线观看 | 精品一区二区三区视频在线观看免费| 亚洲黑人精品在线| 不卡av一区二区三区| 琪琪午夜伦伦电影理论片6080| 九色国产91popny在线| 五月伊人婷婷丁香| 成年女人毛片免费观看观看9| 国产高清视频在线播放一区| 两性夫妻黄色片| 成年版毛片免费区| 久久久久久久午夜电影| 一级毛片女人18水好多| 日韩欧美三级三区| 观看免费一级毛片| 丰满的人妻完整版| 国产久久久一区二区三区| 在线视频色国产色| 久久精品国产99精品国产亚洲性色| 国产爱豆传媒在线观看 | 神马国产精品三级电影在线观看 | 在线视频色国产色| 女人被狂操c到高潮| 男人舔女人的私密视频| 日本一区二区免费在线视频| 日韩免费av在线播放| 两性夫妻黄色片| 日本成人三级电影网站| 亚洲国产欧美一区二区综合| 欧美不卡视频在线免费观看 | 午夜两性在线视频| 日韩高清综合在线| 黑人巨大精品欧美一区二区mp4| 大型av网站在线播放| 亚洲一区二区三区不卡视频| 久久人妻福利社区极品人妻图片| 亚洲国产日韩欧美精品在线观看 | 深夜精品福利| 精品一区二区三区av网在线观看| 在线免费观看的www视频| 日本三级黄在线观看| 五月伊人婷婷丁香| 色av中文字幕| 国产精品香港三级国产av潘金莲| 亚洲avbb在线观看| 男人舔女人下体高潮全视频| 亚洲国产中文字幕在线视频| 亚洲美女视频黄频| 法律面前人人平等表现在哪些方面| 精品久久久久久,| 非洲黑人性xxxx精品又粗又长| 成人av一区二区三区在线看| 1024手机看黄色片| 亚洲成人久久性| 桃红色精品国产亚洲av| 在线观看www视频免费| 久久精品91无色码中文字幕| 黑人巨大精品欧美一区二区mp4| 亚洲黑人精品在线| 亚洲精品一卡2卡三卡4卡5卡| 欧美在线一区亚洲| av免费在线观看网站| 国产精品久久久久久人妻精品电影| 欧美绝顶高潮抽搐喷水| 久久精品国产亚洲av高清一级| 国产伦一二天堂av在线观看| 一夜夜www| 欧美乱色亚洲激情| 午夜亚洲福利在线播放| 国内精品一区二区在线观看| 成人午夜高清在线视频| 亚洲欧美日韩东京热| 国产伦人伦偷精品视频| 少妇熟女aⅴ在线视频| 免费观看人在逋| 一本综合久久免费| 51午夜福利影视在线观看| 欧美日韩国产亚洲二区| 亚洲精品中文字幕一二三四区| 熟妇人妻久久中文字幕3abv| 岛国在线免费视频观看| 精品福利观看| 免费看十八禁软件| 在线观看www视频免费| 久久性视频一级片| 不卡一级毛片| 久99久视频精品免费| 99久久久亚洲精品蜜臀av| 亚洲成av人片免费观看| 日本黄色视频三级网站网址| 制服诱惑二区| 日韩欧美 国产精品| 国产真实乱freesex| 悠悠久久av| 好男人在线观看高清免费视频| 身体一侧抽搐| 黄色片一级片一级黄色片| 老熟妇仑乱视频hdxx| 99在线视频只有这里精品首页| 身体一侧抽搐| 热99re8久久精品国产| 波多野结衣高清无吗| 三级毛片av免费| 91字幕亚洲| 女警被强在线播放| 亚洲乱码一区二区免费版| 久久久久久亚洲精品国产蜜桃av| 国产真人三级小视频在线观看| 成人三级做爰电影| 国产精品久久久人人做人人爽| 精品久久久久久久久久久久久| 国产精品久久久人人做人人爽| 亚洲天堂国产精品一区在线| 美女高潮喷水抽搐中文字幕| 国产精品野战在线观看| 精品久久久久久久毛片微露脸| 亚洲av成人一区二区三| 51午夜福利影视在线观看| 又黄又爽又免费观看的视频| 最近最新免费中文字幕在线| 亚洲人成网站在线播放欧美日韩| 毛片女人毛片| 大型黄色视频在线免费观看| 欧美3d第一页| 亚洲国产看品久久| 正在播放国产对白刺激| АⅤ资源中文在线天堂| 看黄色毛片网站| 国产精品av久久久久免费| 欧美绝顶高潮抽搐喷水| 亚洲五月天丁香| 又紧又爽又黄一区二区| 成年人黄色毛片网站| 三级国产精品欧美在线观看 | 精品高清国产在线一区| 91av网站免费观看| 777久久人妻少妇嫩草av网站| 中国美女看黄片| 国产精品免费视频内射| 一本一本综合久久| 神马国产精品三级电影在线观看 | 免费电影在线观看免费观看| 看免费av毛片| 99国产综合亚洲精品| 欧美乱色亚洲激情| 最近最新免费中文字幕在线| 欧美日韩黄片免| 在线观看免费午夜福利视频| 国产亚洲精品第一综合不卡| 亚洲男人的天堂狠狠| 香蕉av资源在线| 高潮久久久久久久久久久不卡| 国内揄拍国产精品人妻在线| 丰满人妻熟妇乱又伦精品不卡| 国产私拍福利视频在线观看| 成人一区二区视频在线观看| 叶爱在线成人免费视频播放| 又黄又爽又免费观看的视频| 男插女下体视频免费在线播放| 悠悠久久av| 老司机午夜十八禁免费视频| 男女下面进入的视频免费午夜| 级片在线观看| 99在线人妻在线中文字幕| 99国产精品99久久久久| 日本a在线网址| 黄色成人免费大全| 99精品久久久久人妻精品| 国产av在哪里看| 亚洲欧美精品综合久久99| 亚洲美女视频黄频| 看免费av毛片| 制服丝袜大香蕉在线| 欧美精品亚洲一区二区| 久9热在线精品视频| 亚洲熟妇熟女久久| 好看av亚洲va欧美ⅴa在| 亚洲一码二码三码区别大吗| 动漫黄色视频在线观看| 很黄的视频免费| 18禁黄网站禁片免费观看直播| 日韩欧美在线二视频| 五月玫瑰六月丁香| 国产亚洲精品av在线| 99国产精品一区二区蜜桃av| 夜夜看夜夜爽夜夜摸| bbb黄色大片| 搡老岳熟女国产| 国产人伦9x9x在线观看| 色综合站精品国产| 日本熟妇午夜| 久久九九热精品免费| 精品电影一区二区在线| 亚洲欧美精品综合久久99| 91麻豆精品激情在线观看国产| 日日爽夜夜爽网站| 琪琪午夜伦伦电影理论片6080| 怎么达到女性高潮| 亚洲精品国产精品久久久不卡| 国产伦在线观看视频一区| 俄罗斯特黄特色一大片| netflix在线观看网站| 亚洲精品国产精品久久久不卡| 一本大道久久a久久精品| 一级片免费观看大全| 久久久久久免费高清国产稀缺| 精品无人区乱码1区二区| 国产又色又爽无遮挡免费看| 亚洲av成人精品一区久久| 亚洲狠狠婷婷综合久久图片| 老熟妇乱子伦视频在线观看| 久久香蕉激情| 五月伊人婷婷丁香| 久久久精品大字幕| 日本免费a在线| 精品不卡国产一区二区三区| 非洲黑人性xxxx精品又粗又长| 一a级毛片在线观看| 久久久精品大字幕| 午夜亚洲福利在线播放| 久久精品91蜜桃| xxxwww97欧美| 1024香蕉在线观看| 亚洲av成人一区二区三| 1024视频免费在线观看| 91字幕亚洲| 精品久久久久久成人av| 国产黄a三级三级三级人| 日本熟妇午夜| 欧美色视频一区免费| 2021天堂中文幕一二区在线观| 亚洲第一电影网av| 亚洲av电影在线进入| 欧美精品亚洲一区二区| 黄色片一级片一级黄色片| 男女下面进入的视频免费午夜| 怎么达到女性高潮| 男女之事视频高清在线观看| 日本黄大片高清| 成人国语在线视频| 国产成年人精品一区二区| 一本综合久久免费| 久久性视频一级片| 亚洲国产精品sss在线观看| xxx96com| 国产黄色小视频在线观看| 午夜两性在线视频| 欧美丝袜亚洲另类 | 男插女下体视频免费在线播放| 午夜福利在线在线| 亚洲天堂国产精品一区在线| 国产精品精品国产色婷婷| 国产视频一区二区在线看| 热99re8久久精品国产| 婷婷六月久久综合丁香| 午夜日韩欧美国产| 国产精品久久视频播放| 97碰自拍视频| 国产av一区在线观看免费| or卡值多少钱| 亚洲成人免费电影在线观看| 亚洲一区高清亚洲精品| 99热只有精品国产| 搞女人的毛片| 国产欧美日韩一区二区精品| 一级a爱片免费观看的视频| 亚洲欧美日韩高清在线视频| 亚洲av第一区精品v没综合| 亚洲精品中文字幕在线视频| 精品日产1卡2卡| 亚洲精品美女久久久久99蜜臀| 亚洲自偷自拍图片 自拍| 亚洲色图 男人天堂 中文字幕| 日韩成人在线观看一区二区三区| 国产在线精品亚洲第一网站| 日本免费a在线| 久久婷婷成人综合色麻豆| 色av中文字幕| videosex国产| 国产精品亚洲一级av第二区| 18禁黄网站禁片午夜丰满| 99久久精品热视频| 91九色精品人成在线观看| 久久性视频一级片| 国产一区在线观看成人免费| 啦啦啦韩国在线观看视频| 日韩av在线大香蕉| 18禁观看日本| 日本在线视频免费播放| a级毛片a级免费在线| 熟女电影av网| 老司机午夜福利在线观看视频| 哪里可以看免费的av片| 国产精品1区2区在线观看.| 国产日本99.免费观看| 在线免费观看的www视频| 日本熟妇午夜| 亚洲人与动物交配视频| 少妇被粗大的猛进出69影院| 久久精品国产综合久久久| 又大又爽又粗| 午夜福利18| 精品国产乱子伦一区二区三区| 亚洲国产精品成人综合色| www.自偷自拍.com| xxx96com| cao死你这个sao货| 婷婷精品国产亚洲av在线| 99热这里只有是精品50| 久久精品国产亚洲av高清一级| 亚洲专区中文字幕在线| 久久午夜综合久久蜜桃| 国产亚洲精品久久久久5区| 日韩欧美 国产精品| 丰满人妻熟妇乱又伦精品不卡| 欧美日韩亚洲综合一区二区三区_| 麻豆久久精品国产亚洲av| √禁漫天堂资源中文www| 成在线人永久免费视频| 在线观看www视频免费| 欧美另类亚洲清纯唯美| 舔av片在线| 人妻夜夜爽99麻豆av| 亚洲精品国产一区二区精华液| 精品久久久久久久毛片微露脸| 亚洲中文字幕一区二区三区有码在线看 | 91大片在线观看| 岛国在线免费视频观看| 久久久久亚洲av毛片大全| 国产在线精品亚洲第一网站| 日韩av在线大香蕉| 免费看a级黄色片| 日日夜夜操网爽| 又大又爽又粗| 国产不卡一卡二| 欧美乱码精品一区二区三区| 麻豆国产av国片精品|