亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        2 840 MHz太陽射電輻射流量計(jì)的升級(jí)改造

        2016-07-11 01:35:48顏毅華徐之材殷興輝譚程明陳志軍
        天文研究與技術(shù) 2016年1期
        關(guān)鍵詞:輻射計(jì)

        杜 靜,顏毅華,徐之材,殷興輝,王 威,譚程明,陳志軍

        (1. 中國科學(xué)院國家天文臺(tái)太陽活動(dòng)重點(diǎn)實(shí)驗(yàn)室,北京 100012; 2. 中國科學(xué)院大學(xué),北京 100049;

        ?

        2 840 MHz太陽射電輻射流量計(jì)的升級(jí)改造

        杜靜1,2,顏毅華1,徐之材3,殷興輝4,王威1,譚程明1,陳志軍1

        (1. 中國科學(xué)院國家天文臺(tái)太陽活動(dòng)重點(diǎn)實(shí)驗(yàn)室,北京100012; 2. 中國科學(xué)院大學(xué),北京100049;

        3. 南京射電天文技術(shù)研究所,江蘇 南京210024; 4. 河海大學(xué),江蘇 南京211100)

        摘要:敘述了國家天文臺(tái)升級(jí)改造后的2 840 MHz太陽射電輻射流量計(jì)的設(shè)計(jì)特點(diǎn)、性能、結(jié)構(gòu)及觀測結(jié)果。升級(jí)改造后的太陽射電輻射流量計(jì)將在較高時(shí)間分辨率上實(shí)時(shí)得到2 840 MHz頻率上的太陽流量,為太陽物理研究積累豐富的觀測數(shù)據(jù),是太陽活動(dòng)監(jiān)測和預(yù)報(bào)的重要參數(shù)之一。

        關(guān)鍵詞:射電天文;太陽射電輻射流量計(jì);流量觀測;輻射計(jì)

        太陽微波輻射流量與黑子相對(duì)數(shù)、耀斑爆發(fā)、X射線爆發(fā)以及高能電子質(zhì)子事件均有一定的關(guān)系[1],而2 840 MHz(波長為10.7 cm)的太陽射電輻射流量與黑子相對(duì)數(shù)及高能電子和質(zhì)子事件的關(guān)系最為密切,是地球空間物理的重要輸入?yún)?shù)。因此,2 840 MHz的射電輻射流量監(jiān)測對(duì)研究太陽活動(dòng)規(guī)律、太陽長期演化行為,建立更加完善的太陽活動(dòng)預(yù)報(bào)模型以及太陽射電天文學(xué)的發(fā)展等具有極其重要的意義。國家天文臺(tái)懷柔2 840 MHz太陽射電輻射流量計(jì)長期參與國際聯(lián)網(wǎng)觀測,連續(xù)三十多年為國際太陽和地球物理數(shù)據(jù)中心提供觀測數(shù)據(jù),為國際太陽活動(dòng)預(yù)報(bào)提供可靠的數(shù)據(jù)支持。文[2]利用懷柔2 840 MHz太陽射電輻射流量計(jì)的長期觀測數(shù)據(jù)研究太陽活動(dòng)周的相似性行為,并預(yù)測第24活動(dòng)周的開始時(shí)間和最大活動(dòng)強(qiáng)度。文[3]利用2 840 MHz射電輻射流量的長期觀測數(shù)據(jù)研究了太陽長期活動(dòng)行為,發(fā)現(xiàn)太陽活動(dòng)的中短時(shí)間尺度的周期性與太陽系行星運(yùn)動(dòng)之間有較強(qiáng)的相關(guān)性。同時(shí),2 840 MHz太陽射電輻射流量計(jì)也是中國空間環(huán)境監(jiān)測網(wǎng)的重要組成部分,為國家航空航天任務(wù)提供空間環(huán)境保障發(fā)揮了重要作用。

        2 840 MHz太陽射電輻射流量計(jì)位于國家天文臺(tái)懷柔太陽觀測基地懷柔水庫北岸,距離北京市中心約60 km,北緯40.3°,東經(jīng)116.6°,海拔高度64 m。附近有較強(qiáng)的機(jī)場雷達(dá)信號(hào)干擾,但干擾信號(hào)特征明顯,數(shù)據(jù)處理過程中,可通過成熟的數(shù)據(jù)處理方法予以剔除。21世紀(jì)初,該設(shè)備由于長年運(yùn)行,多數(shù)元器件老化停止了觀測,考慮到2 840 MHz太陽射電輻射流量計(jì)是國內(nèi)唯一的單頻射電輻射流量計(jì)及其在太陽活動(dòng)研究方面的重要意義,經(jīng)過充分調(diào)研,對(duì)2 840 MHz太陽射電輻射流量計(jì)進(jìn)行了升級(jí)改造。新設(shè)備對(duì)包括微波前端(從低噪放開始)、模擬接收系統(tǒng)、數(shù)字采集系統(tǒng)、存儲(chǔ)及顯示系統(tǒng)都進(jìn)行了升級(jí)改造,天線采用2 m拋物面反射面天線,赤道式裝置。同時(shí),采用了K因子輻射計(jì)和新的定標(biāo)技術(shù),使得定標(biāo)更加高效,定標(biāo)結(jié)果更加精確可信[4]。目前,該設(shè)備已獲取了一些良好的觀測結(jié)果。

        1系統(tǒng)指標(biāo)

        太陽射電望遠(yuǎn)鏡系統(tǒng)是自動(dòng)化程度較高的空間探測裝備,具有定標(biāo)系統(tǒng)、自動(dòng)采集觀測數(shù)據(jù)、自動(dòng)生成和傳送設(shè)備,同時(shí)具有較高的可靠性、穩(wěn)定性,能全天候長期、連續(xù)地運(yùn)行。系統(tǒng)指標(biāo)如表1。

        2系統(tǒng)結(jié)構(gòu)

        2 840 MHz太陽射電輻射流量主要用于太陽活動(dòng)長期行為的研究和太陽預(yù)報(bào),對(duì)測量精度和數(shù)據(jù)的穩(wěn)定性要求很高,因此,在本次升級(jí)改造過程中采用鏡像抑制的超外差式接收機(jī)方案,可很好地抑制干擾,選用K因子輻射計(jì)保證儀器的靈敏度和穩(wěn)定性。同時(shí),近年來關(guān)于太陽活動(dòng)的短期行為也不斷受到太陽物理學(xué)家的高度重視,而2 840 MHz太陽射電流量有最靈敏的響應(yīng),因此,為了滿足這方面的科學(xué)需求,采用了衰減器和對(duì)數(shù)放大器2種方法以提高數(shù)據(jù)的動(dòng)態(tài)范圍。系統(tǒng)結(jié)構(gòu)如圖1。

        3系統(tǒng)特點(diǎn)

        3.1良好的抗干擾設(shè)計(jì)

        采用鏡像抑制的超外差式接收機(jī),在混頻器前加2 840 MHz的聲表面帶通濾波器,選用較高中頻fo=225 MHz,中頻帶寬BWIF=10 MHz,采用BW/fo=5%的聲表面帶通濾波器,矩形系數(shù)接近1,帶外衰減大,對(duì)鏡頻可抑制20~30 dB,從而有效抑制鏡頻干擾和帶外干擾,保證觀測數(shù)據(jù)的可靠性。同時(shí),采用合理的微波放大器,有效抑制交調(diào)互調(diào)干擾信號(hào)的產(chǎn)生,采用良好的電磁兼容性設(shè)計(jì),防止干擾信號(hào)寄生耦合到信號(hào)通道。模數(shù)轉(zhuǎn)換器放在接收機(jī)機(jī)箱內(nèi),與低放輸出短線相接,減小50 Hz等低頻干擾。

        圖12 840 MHz太陽射電輻射流量計(jì)方框圖

        Fig.1The block diagram of the upgraded Solar Radio

        Flux Telescope working at 2840MHz

        3.2大動(dòng)態(tài)范圍

        采用衰減器和對(duì)數(shù)放大器2種方法擴(kuò)大測量動(dòng)態(tài)范圍的技術(shù),構(gòu)成2個(gè)并行輸出的信號(hào)。衰減器法采用1 dB步進(jìn)的最大衰減量為31 dB的數(shù)控步進(jìn)衰減器,用它調(diào)整中頻輸出信號(hào)電平,數(shù)控步進(jìn)衰減器由計(jì)算機(jī)自動(dòng)控制[5]。衰減器法測寧靜太陽和小爆發(fā)的精度較高,但觀測大爆發(fā)時(shí)由于衰減量變化,產(chǎn)生觀測數(shù)據(jù)的跳變,雖然計(jì)算處理后能恢復(fù)連續(xù)變化特性,但精度有所下降。對(duì)數(shù)放大器法可獲得連續(xù)變化的記錄,特別適用于記錄大爆發(fā),精度可達(dá)0.2 dB。

        3.3采用K因子輻射計(jì)

        K因子輻射計(jì)是根據(jù)德國天文學(xué)家Keen提出的K因子測量原理而發(fā)展的新型輻射計(jì),其信號(hào)路工作在全功率計(jì)狀態(tài),靈敏度高。采用天線溫度定標(biāo)的噪聲源與K因子輻射計(jì)調(diào)制式噪聲源合一的方案,通過測量定標(biāo)常數(shù)Cs對(duì)系統(tǒng)進(jìn)行定標(biāo),實(shí)現(xiàn)觀測與溫度標(biāo)定同時(shí)自動(dòng)進(jìn)行,提高了標(biāo)定精度[6-7]。升級(jí)改造后的2 840 MHz太陽射電輻射流量計(jì)采用K因子式輻射計(jì),其原理圖如圖2,圖中對(duì)應(yīng)部件代碼和符號(hào)見表2。

        3.3.1K因子輻射計(jì)的基本關(guān)系

        采用理想平方律檢波器,且放大器工作在線性態(tài)時(shí),VDC和VAC兩路輸出與Ta、Tr和ΔT的關(guān)系[8-9]:

        (1)

        (2)

        K因子定義:K=VDC/VAC,將(1)、 (2)代入上述關(guān)系式,得:K=(Ta+Tr+ΔT)/C,式中:

        故Ta=CK-Tr-ΔT。

        3.3.2用K因子輻射計(jì)觀測射電源

        (1)對(duì)準(zhǔn)射電源時(shí):得到的K因子為Ks,此時(shí),Ta=Ts+Tb,Ts+Tb=CKs-Tr-ΔT。

        (2)對(duì)準(zhǔn)射電源背景時(shí):得到的K因子為Kb,此時(shí),Ta=Tb,Tb=CKb-Tr-ΔT。

        以上兩式相減得到:

        (3)

        Ts為天線溫度中射電源分量;Tb為天線溫度中背景噪聲溫度分量。

        3.3.3定標(biāo)常數(shù)Cs:

        圖2K因子輻射計(jì)原理圖

        Fig.2The block diagram of the radiometer

        for the K-factor measurement

        圖3K因子輻射計(jì)天線口面的噪聲

        分量和時(shí)變狀態(tài)圖

        Fig.3The diagram of the noise components and time

        varying state for the antenna of the radiometer

        with the K-factor measurement

        當(dāng)天線的波束寬度遠(yuǎn)大于射電源視角徑時(shí),Ts=AeS/2k,

        (4)

        Ae為天線的有效接收面積;k為玻爾茲曼常數(shù)。

        結(jié)合(3)、 (4)兩式:C(Ks-Kb)=AeS/2k,

        由此可得到:S=Cs(Ks-Kb),

        (5)

        (6)

        實(shí)際操作中,首先采用其它射電望遠(yuǎn)鏡獲得太陽等射電源流量S,然后分別對(duì)太陽和背景進(jìn)行觀測獲得Ks和Kb,相對(duì)測量的望遠(yuǎn)鏡與它比對(duì)后,依據(jù)(6)式可獲得定標(biāo)常數(shù)Cs。 將Cs數(shù)值輸入定標(biāo)軟件中作為配置參數(shù)使用,之后分別對(duì)背景和太陽進(jìn)行觀測,依據(jù)(5)式即可獲得太陽流量值,相對(duì)測量轉(zhuǎn)換為絕對(duì)測量。

        4系統(tǒng)性能

        4.1靈敏度

        在射電天文觀測中,靈敏度是射電望遠(yuǎn)鏡的重要指標(biāo),決定了接收系統(tǒng)檢測弱信號(hào)的能力,定義為 “最低可測” 的輻射流量密度Smin[1,10]。

        寧靜太陽射電輻射流量:S=80 sfu*http://www.spaceweather.ca/solarflux/sx-5-flux-eng.php?year=2011

        天線有效面積:Ae=ηπd2/4=2.67m2,式中,天線直徑d=2 m;天線效率η=85%

        天線溫度的太陽分量:Ts=AeS/2k=2.67×80×10-22/(2×1.38×10-23)=774K

        輻射計(jì)的工作參數(shù)見表3。

        K因子輻射計(jì)理論靈敏度 (S/N=1):總插損Li=0.89

        天線溫度:Ta=Ts+Tb=804K

        系統(tǒng)噪聲溫度:

        Tsys=LiTa+(1-Li)To+Tr+ΔT=0.89×804+(1-0.89)×290+290+290=1 327K

        噪聲起伏:ΔTmin=Tsys/(2Bwτ)1/2=0.297K/Hz

        系統(tǒng)靈敏度:ΔSmin=2kTmin/Ae=2×1.38×10-23×0.297/2.67=0.031sfu

        4.2計(jì)算機(jī)系統(tǒng)及軟件功能

        (1)用4路A/D采集監(jiān)測電壓Vm、直流輸出電壓Vdc、交流輸出電壓Vac和對(duì)數(shù)放大器輸出電壓Vlog。根據(jù)Vm值控制步進(jìn)衰減器,令接收機(jī)處于合理工作狀態(tài)。根據(jù)人工輸入的 ‘S’(觀測太陽)和 ‘B’(觀測背景)確定Vdc和K因子是Vdcs、Ks還是Vdcb、Kb。進(jìn)行數(shù)字積分時(shí),標(biāo)準(zhǔn)的積分時(shí)間為1 ms和1 s兩種,兩種數(shù)據(jù)可同時(shí)顯示和記錄。升級(jí)改造后的2 840 MHz太陽射電輻射流量計(jì)的觀測軟件界面如圖4。

        (2)進(jìn)行預(yù)處理,如:

        ① 步進(jìn)衰減器的跳變式記錄變?yōu)檫B續(xù)漸變式記錄;

        ② 顯示實(shí)時(shí)的太陽流量;

        ③ 顯示全天太陽流量變化;

        ④ 對(duì)數(shù)值變換為比值顯示。

        (3)數(shù)據(jù)存儲(chǔ):可同時(shí)存儲(chǔ)積分時(shí)間為1 ms(.nwm)和1 s(.nws)的觀測數(shù)據(jù)。采用整型數(shù)據(jù)存儲(chǔ),包括64字節(jié)文件頭和觀測數(shù)據(jù),每2字節(jié)表示一個(gè)采樣點(diǎn)(觀測數(shù)據(jù)點(diǎn))。毫秒級(jí)文件是每10 s存一個(gè)文件,每個(gè)文件20 KB,每天的數(shù)據(jù)量約為56 M,秒級(jí)文件每天的數(shù)據(jù)量約為60 KB (按每天觀測8 h計(jì)算)。

        (4)數(shù)據(jù)發(fā)布:望遠(yuǎn)鏡從1977年開始觀測,2010年進(jìn)行設(shè)備的升級(jí)改造,其數(shù)據(jù)已連續(xù)積累三十多年,定期發(fā)布在國家天文臺(tái)太陽數(shù)據(jù)中心(http://csrh.bao.ac.cn/data/),用戶可自由下載。

        圖4觀測軟件界面

        Fig.4The interface of the observation software

        5觀測結(jié)果

        利用經(jīng)過上述升級(jí)改造后的2 840 MHz太陽射電輻射流量計(jì)對(duì)太陽進(jìn)行觀測,成功觀測到2011年2月15日X2.2級(jí)大耀斑爆發(fā)(圖5)。該事件是一個(gè)典型的長持續(xù)(近1個(gè)小時(shí))復(fù)雜型(多峰)射電爆發(fā)事件。復(fù)雜的多峰結(jié)構(gòu)對(duì)應(yīng)同頻率的頻譜觀測中有斑馬紋、脈動(dòng)、慢漂移等豐富的精細(xì)結(jié)構(gòu)[11]。

        6結(jié)論

        本文介紹了國家天文臺(tái)懷柔觀測站2 840 MHz太陽輻射流量計(jì)的升級(jí)改造。新設(shè)備采用了K因子輻射計(jì)及新的定標(biāo)技術(shù),對(duì)微波前端、模擬接收系統(tǒng)、數(shù)字采集系統(tǒng)、存儲(chǔ)顯示系統(tǒng)也進(jìn)行了升級(jí)改造。通過設(shè)備的升級(jí)改造,改進(jìn)了天線的觀測性能,其數(shù)據(jù)對(duì)研究耀斑物理動(dòng)力學(xué)過程、能量釋放、太陽活動(dòng)規(guī)律及太陽長期演化行為等具有重要意義,為相關(guān)研究提供了豐富的觀測依據(jù)。

        圖52011.2.15 X2.2級(jí)耀斑-新系統(tǒng)

        Fig.5The X2.2 solar flareon Feb.15, 2011

        observed with the new system

        致謝:衷心感謝設(shè)備安裝、調(diào)試過程中南京射電天文技術(shù)研究所吳楓、河海大學(xué)沈亞明等的大力支持和幫助。

        參考文獻(xiàn):

        [1]趙仁揚(yáng), 金聲震, 傅其駿. 太陽射電微波爆發(fā)[M]. 北京: 科學(xué)出版社, 1997.

        [2]Wang Jialong, Zong Weiguo, Le Guiming, et al. Predicting the start and maximum amplitude of solar cycle 24 using similar phases and a cycle grouping[J]. Research in Astronomy and Astrophysics, 2009, 9(2): 133-136.

        [3]Tan Baolin, Cheng Zhuo. The mid-term and long-term solar quasi-periodic cycles and the possible relationship with planetary motions[J]. Astrophysics & Space Science, 2013, 343(2): 511-521.

        [4]彭樹生. 提高微波輻射計(jì)絕對(duì)測量精度的研究[D]. 南京: 中國科學(xué)院紫金山天文臺(tái), 1995.

        [5]董亮, 何樂生, 施碩彪. 太陽射電800~975 MHz波段高動(dòng)態(tài)范圍模擬接收機(jī)的研制[J]. 天文研究與技術(shù)——國家天文臺(tái)臺(tái)刊, 2014, 11(2): 111-117.

        Dong Liang, He Lesheng, Shi Shuobiao. Design of a high dynamic range analog receiver for solar radio observation in 800-975MHz[J]. Astronomical Research & Technology——Publications of National Astronomical Observatories of China, 2014, 11(2): 111-117.

        [6]楊精一. 微波輻射計(jì)定標(biāo)方程研究[J]. 長春光學(xué)精密機(jī)械學(xué)院學(xué)報(bào), 2001, 24(3): 21-24.

        Yang Jingyi. Research of the establish standard of equation with the microwave radiator[J]. Journal of Changchun Institute of Optics and Fine Mechanics, 2001, 24(3): 21-24.

        [7]劉麗佳, 劉彬, 董亮, 等. 云南天文臺(tái)40米射電望遠(yuǎn)鏡S波段總流量觀測系統(tǒng)研制和測試[J]. 天文研究與技術(shù)——國家天文臺(tái)臺(tái)刊, 2013, 10(2): 134-141.Liu Lijia, Liu Bin, Dong Liang, et al. Development and testing of an S-band flux measurement system for the 40m radio telescope of the Yunnan Observatory[J]. Astronomical Research & Technology——Publications of National Astronomical Observatories of China, 2013, 10(2): 134-141.

        [8]彭樹生, 殷興輝, 徐之材. 五毫米大氣遙感微波輻射計(jì)[J]. 遙感技術(shù)與應(yīng)用, 1995, 10(2): 38-41.

        Peng Shusheng, Yin Xinghui, Xu Zhicai. 5mm microwave radiometer for atmosphere remote sensing[J]. Remote Sensing Technology and Application, 1995, 10(2): 38-41.

        [9]吳禮. 毫米波輻射計(jì)特性參數(shù)與性能測試研究[D]. 南京: 南京理工大學(xué), 2004.

        [10]王威, 竇玉江, 顏毅華, 等. CSRH靈敏度分析[J]. 天文研究與技術(shù)——國家天文臺(tái)臺(tái)刊, 2013, 10(1): 22-25.

        Wang Wei, Dou Yujiang, Yan Yihua, et al. Analysis of the sensitivity of the CSRH[J]. Astronomical Research & Technology——Publications of National Astronomical Observatories of China, 2013, 10(1): 22-25.

        [11]Tan Baolin, Yan Yihua, Tan Chengming, et al. Microwave zebra pattern structures in the X2.2 solar flare on 2011 February 15[J]. Astrophysical Journal, 2012, 744(2): 1-10.

        CN 53-1189/PISSN 1672-7673

        Upgrading the Solar Radio Flux Telescope with the Frequency of 2.84GHz

        Du Jing1,2, Yan Yihua1, Xu Zhicai3, Yin Xinghui4, Wang Wei1, Tan Chengming1, Chen Zhijun1

        (1. Key Laboratory of Solar Activity, National Astronomical Observatories, Chinese Academy of Sciences, Beijing 100012, China,Email: jdu@nao.cas.cn; 2. University of Chinese Academy of Sciences, Beijing 100049, China; 3. Nanjing Institute of Radio Astronomical Technology, Nanjing 210024, China; 4. Hehai University, Nanjing 211100, China)

        Key words:Radio astronomy; Solar radio flux telescope; Flux measurement; Radiometer

        Abstract:The solar radio flux telescope with the frequency of 2.84GHz is an important instrument to monitor and forecast the solar activities by collecting and analyzing the solarradio fluxes. At the beginning of the 21st century, the regular monitoring by using such an instrument was stopped due to ageing of most components and difficult to repair during the long running. Now we have updated the microwave front-end, the analog receiving system, the digital acquisition system, the storage and display system, and the 2m parabolic antenna to the system. In the new system we adopted a radiometer for k-factor measurement and a new self-calibration technique which makes the calibration more efficient and accurate. At present, the upgraded solar radio flux telescope working at 2.84GHz has achieved some good results.

        基金項(xiàng)目:國家自然科學(xué)基金 (11103044, 11273030, 11221063, 11433006);科技部科學(xué)基礎(chǔ)性工作專項(xiàng) (2014FY120300) 和國家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃 (973計(jì)劃) (2011CB11401) 資助.

        收稿日期:2015-03-28;

        修訂日期:2015-04-16

        作者簡介:杜靜, 女, 博士. 研究方向:太陽射電. Email: jdu@nao.cas.cn

        中圖分類號(hào):P111.44

        文獻(xiàn)標(biāo)識(shí)碼:A

        文章編號(hào):1672-7673(2016)01-0052-06

        猜你喜歡
        輻射計(jì)
        HY-2B掃描微波輻射計(jì)海溫產(chǎn)品在典型區(qū)域的精度驗(yàn)證和分析
        基于微波輻射計(jì)的張掖地區(qū)水汽、液態(tài)水變化特征分析
        風(fēng)云四號(hào)A星多通道掃描成像輻射計(jì)第一幅彩色合成圖像
        上海航天(2021年3期)2021-07-21 14:02:12
        基于CLEAN算法對(duì)一維綜合孔徑輻射計(jì)成像誤差的校正
        MP-3000 型地基微波輻射計(jì)探測性能分析
        毫米波被動(dòng)探測系統(tǒng)中彈目交會(huì)模型建模與仿真
        微波輻射計(jì)對(duì)昆明霧天的監(jiān)測及應(yīng)用
        高靈敏度、高穩(wěn)定度微波輻射計(jì)技術(shù)研究與實(shí)驗(yàn)驗(yàn)證
        太陽漫反射板衰減監(jiān)測輻射計(jì)的設(shè)計(jì)及性能評(píng)估
        用于目標(biāo)輻射特性測量的W波段輻射計(jì)
        成年无码aⅴ片在线观看| 国产自拍在线视频91| 天堂在线资源中文在线8| av潮喷大喷水系列无码| AV教师一区高清| 国产精品一区又黄又粗又猛又爽| 亚洲精品视频在线一区二区| 中文www新版资源在线 | 亚洲欧美日韩国产一区二区精品| 亚洲人妻av在线播放| 亚洲av精二区三区日韩| 欧美bbw极品另类| 亚洲高清在线不卡中文字幕网| 蜜桃人妻午夜精品一区二区三区| 国产成人a级毛片| 午夜一区欧美二区高清三区| 国产三级精品美女三级| 人妻在线有码中文字幕| 正在播放老肥熟妇露脸| 国产精品原创巨作AV女教师 | 久久精品中文少妇内射| 男女一边摸一边做爽爽的免费阅读| 亚洲精品aⅴ无码精品丝袜足| 国产精品美女一区二区av| 欧美激情一区二区三区成人 | av在线色| 亚洲av一二三四又爽又色又色| 亚洲国产精品久久久av| 亚洲精品黑牛一区二区三区| 日韩成人大屁股内射喷水| 亚洲国产精品悠悠久久琪琪| 女人天堂国产精品资源麻豆| 亚洲欧美国产精品久久| 秒播无码国产在线观看| 一区二区三区国产亚洲网站| 欧美成人家庭影院| 亚洲色无码播放| 亚洲av一区二区国产精品| 国产熟女一区二区三区不卡| 无码国模国产在线观看| 国产清品夜色一区二区三区不卡|