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

    The New Wuqing 70m Radio Telescope and Measurements of Main Electronic Properties in the X-band

    2022-05-24 14:21:02DeQingKongChunLaiLiHongBoZhangYanSuJianJunLiuXinYingZhuJunDuoLiXiPingXueandChenLi

    De-Qing Kong ,Chun-Lai Li,* ,Hong-Bo Zhang ,Yan Su ,Jian-Jun Liu ,Xin-Ying Zhu ,Jun-Duo Li,Xi-Ping Xue,and Chen Li

    1 Key Laboratory of Lunar and Deep Space Exploration,National Astronomical Observatories,Chinese Academy of Sciences,Beijing 100101,China;licl@bao.ac.cn

    2 University of Chinese Academy of Sciences,Beijing 100049,China

    AbstractThe new Wuqing 70 m radio telescope is first used for the downlink data reception in the first Mars exploration mission of China,and will be used for the other deep space communications and radio astronomical observations in the future.The main specifications and measurement results of some properties in the X-band are introduced in this paper,such as pointing calibration,gain and efficiency,system noise temperature,system equivalent flux density,and variations with elevation.The 23 parameters pointing calibration model considering the atmospheric refraction correction in real time is presented in the telescope,and the pointing accuracy reaches in azimuth direction and in elevation direction for different weather conditions.More than 62% efficiencies are achieved at full elevation range,and more than 70%in the mid-elevation.The system equivalent flux density of the X-band in the mid-elevation reaches 26 Jy.

    Key words:telescopes–astronomical instrumentation,methods and techniques–atmospheric effects

    1.Introduction

    The Wuqing 70 m radio telescope (WRT70),70 m in diameter,is located at the Wuqing Station of the National Astronomical Observatories,longitude 117.10 E,latitude 39.54 N,and elevation 8 m.The WRT70 is a turntable Cassegrain antenna with an Azimuth-Zenith mounting support-structure,and is currently the largest fully movable single dish antenna in Asia.The construction of the WRT70 started in 2018 October,started trial operation in 2020 October and was passed the final reception test in 2021 February.

    The WRT70 is first used for the downlink data reception in the first Mars exploration mission of China,and will be used for other deep space communications and radio astronomical observations in the future.The frequency of the downlink signal is X-band in the first Mars exploration mission of China.In the subsequent missions,Ka-band will be used for the downloading of the scientific data.So in the design of the WRT70,the position of the Ka-band is reserved,and the accuracy of the mount and panel has been also considered.The WRT70 has three feeds working at S/X,X and Ku-bands now,and Ka and K or C-bands will be installed later.The feeds can be rotated quickly to the right position automatically to change the observing frequency.In the Mars and following deep space probes,the WRT70 is arrayed with the Miyun 50 m,Miyun 40 m and Kunming 40 m radio telescopes,which makes the equivalent aperture of the synthesized antenna reach 103 m.

    It is necessary to test the main electronic properties for a new telescope,such as pointing calibration,gain and efficiency,system noise temperature,system equivalent flux density,and variations with elevation.For large telescopes,pointing calibration needs to be carried out first(Meeks et al.1968;Yuan et al.1986;Himwich 1993).The linear calibration model is applied in the pointing calibration of mostly large radio telescopes(Meeks et al.1968;Yuan et al.1986;Zhang et al.2009).However,the linear model does not take the structure nonlinear errors into considerations.The pointing model based on Least-square support vector machines was presented in 2007 (Zhang et al.2007),and the model based on generalized extended interpolation was presented and used in the Miyun 50 m radio telescope in 2008(Kong et al.2008),which can partially solve the nonlinear problem.The model considering the nonlinear error of azimuthaxis was presented and used in the Kunming 40 m radio telescope in 2014(Kong et al.2014).Due to the improvement of pointing accuracy of large antennas,higher requirements are put forward for the atmospheric refraction model.

    Calibrations and measurements of the mainly properties of the WRT70 in the X-band were performed from 2020 November to 2021 January.In Section 2,the results show that the WRT70 fully meets the design specifications and has excellent performance compared with large-diameter antennas.The pointing calibration method and result will be introduced in Section 3.The antenna gain and efficiency are measured at full elevation by radio sources.The results of efficiency and gain will be introduced in Section 4.The results of system temperature and its variations with elevation will be introduced in Section 5.The results of the system equivalent flux density and variations with elevation using the Y-factor method by radio source will be introduced in Section 6.

    2.The Main Specifications of the WRT70

    The parabolic reflector is made from solid plate in all the 70 m region.The accuracy of the main surface is 0.35 mm(elevation 48°,by actual measurement),and the accuracy of the secondary reflecting surface is 0.16 mm which is relatively easy to ensure.The main reflective surface dispersion support structure with umbrella-shaped reflector back frame is adopted.The influence of uniformity improves the accuracy of the reflection surface.In the accuracy measurement and adjustment of the main reflection surface,the innovative combination of photographic method and laser scanning method is used,which not only improves the accuracy and speed of surface measurement,but also makes it easy to achieve surface measurement at different elevation angles.The photogrammetry is used in the pre-adjustment of the antenna panel,which can only be carried out at night and takes a long time.The advantage of the laser scanner is that the measurement time is as short as a few minutes.The laser scanner is used to measure the main and sub-reflector surface deformation at different elevation angles,which is mainly used for the preliminary adjustment curves of the sub-surface and the research of the measurement of sunshine temperature surface deformation.

    The main and sub-reflector surface shaped design technology is adopted,which is to optimize the antenna parameters,the radiation level at the edge of the sub-reflector surface,and the aperture field distribution function of the primary reflecting surface.The sub-reflector surface follow-up and the best anastomosis technologies are adopted in the WRT70 design,which ensuring the accuracy of the main reflector surface in the full elevation.The sub-reflector has three adjustment degrees of freedom which are Y and Z-axis translation,and rotation around X-axis.The repetitive positioning accuracy of Y,Z and X-axis is 0.08 mm,0.08 mm and 5″respectively,then the pointing error caused by the sub-reflector repetitive accuracy is less thanAs a result,the calculated antenna efficiency of the X-band is greater than 70% near the panel adjusting positions,and more than 60% within the full elevation.

    Because of the huge structure,the long control line and the very heavy wind load,the servo control equipment is susceptible to interference and the control accuracy is not easy to guarantee.The full digital control loop technology and multi-motor (azimuth eight motors,pitch four motors)electric anti-gap and the Linear Quadratic Gaussian anti-gust disturbance control technologies are adopted,and the high control accuracy of 3″both azimuth and elevation is achieved.Figure 1 shows the picture of the WRT70 radio telescope.

    Figure 1.The WRT70 Radio Telescope.

    Figure 2.Partial temperature and stress sensors distribution.(a)Temperature sensors on the alidade,(b)temperature sensors on the main reflector,(c)temperature and stress sensors on the sub-reflector leg.

    Because large antennas require high pointing accuracy and are easily affected by environmental factors,such as atmospheric refraction,wind and sunshine temperature,temperature sensors,inclinometers and weather meters are placed on key parts of the reflector and alidade.Two inclinometers are respectively installed near the elevation bearing of the alidade.52 temperature sensors and eight stress sensors are placed on the telescope,which will be used for long-term monitoring and correction experiments on structural deformation.The stress sensors are used for the research on wind load and gravity deformation,and the temperature sensors are used for the research on sunshine temperature deformation.Figure 2 shows the temperature and stress sensors distribution.A 23 parameters correction model considering nonlinearity and the atmospheric refraction correction in real time improves the pointing accuracy of the WRT70,which can be used in real scientific observation.

    Table 1 gives the main specifications of the WRT70.

    Table 1 The Main Design Specifications of the WRT70

    3.The Pointing Calibration

    Pointing accuracy is one of the most important performance indicators in large radio telescopes.System pointing calibration is the first step for a new telescope.The widely used linear pointing calibration model is expressed as (Zhang et al.2009)

    Here,δA1and δE1are the azimuth and elevation calibration functions,the correction factor C1~C8can be identified by the least squares method.

    For large radio telescopes,it is necessary to adopt a more complex correction model,because of the nonlinear errors produced by the antenna structure,the influence of sunshine and wind load.For the WQT70,based on the pointing model considering the nonlinear error of azimuth-axis of Kunming 40 m radio telescope (Kong et al.2014),we add the azimuth and elevation encoder scale error terms to correct the errors due to coaxiality of encoder installation.The pointing calibration model we used is expressed as

    The parameters a1~a12and e1~e11in Equation (2) are correction factors,and can be identified by the least squares method,and a6and e6are respectively the azimuth and elevation encoder scale differences.The traditional atmospheric refraction correction model is only for one atmospheric condition without considering the influence of real-time changing factors such as atmospheric temperature,atmospheric pressure,and relative humidity.The troposphere atmospheric refraction model in radio band given by Berman describes the refraction angle at different elevation angles (Berman &Rockwell 1976).The atmospheric refraction angle of the entire troposphere to electromagnetic waves is calculated using the temperature,atmospheric pressure and relative humidity measured in real time on the ground.So the following elevation pointing calibration model considering the atmospheric refraction correction in real time is used in the WRT70

    Here (Berman &Rockwell 1976)

    Here R is the angle of refraction,Z is the actual zenith angle,E is the elevation angle,P is the local atmospheric pressure,T is the local temperature,RH is the relative humidity,and other parameters are constant terms.

    The observation was taken in the X-band,since the Ku-band receiver is not installed at that time.Sources selected for calibration observation are strong intensity,small angular diameter,and uniform distribution in the sky.The selected radio sources are 3C 273B,3C 286,3C 380,3C 279,3C 380,3C 345,DR21,3C 84,etc.

    The pointing errors in total 149 positions,were observed in 2021 November in the first round pointing correction.The pointing calibration results of the WRT70 are shown in Table 2,where e1is the residual pointing error after the calibration using the linear pointing calibration model(Equation (2)),and e2is the error using the calibration model(Equation(3)).It can be seen from Table 2 that e2is better than e1.Particularly,the elevation-axis pointing accuracy increasesfromand improves by 12.4%.The overall pointing accuracy increases fromand improves by 8.7%.Figure 3 shows the sky coverage of the pointing calibration observations and the distribution of the pointing calibration residues of the first round pointing calibration.

    Table 2 Comparison of the First Round Pointing Residuals (Rms) of the WRT70

    In order to be suitable for different weather conditions,we carried out the second round of measurements considering the correction of atmospheric refraction.The pointing errors in total 199 positions were observed in 2021 January.In order to ensure the pointing correction has a wide range of adaptability,observation time covered different weather conditions such as sunny,cloudy,moderate rain and snow.The influenced by the uneven temperature in sunny daytime was not considered in pointing calibration.The pointing calibration results of the WRT70 is shown in Table 3,where e1is the residual pointing error after the calibration using the linear pointing calibration model(Equation(1)),e2is the error using the calibration model(Equation (2)),and e3is the error using the calibration model(Equation(3))with the real-time atmospheric refraction model.Figure 4 shows the pointing sky coverage,the calibration residues,temperature,relative humidity and atmospheric pressure distributions of the second round calibration.

    Figure 3.The sky coverage of the pointing calibration observations and the distribution of the pointing calibration residues of the first round pointing calibration.

    Table 4 shows the rms values of the two models(Equations(2)and(3)).The overall pointing accuracy increases fromand improves by 9.8% by using the realtime atmospheric refraction model.In particular,in the observation of low elevation angles below 20°,the pointing correction residual reduces by 19.3%,23.1% lower than 15°,and 27.5%lower than 10°.The new model makes the elevation correction residuals at different elevation more at the same level,the low elevation angle pointing accuracy deterioration caused by atmospheric refraction is basically eliminated.

    Table 3 Comparison of the Second Round Pointing Residuals (Rms) of the WRT70

    Table 4 Pointing Calibration Residuals Under Certain Elevations

    4.Gain and Gain Variation with Elevation

    The gain curve with elevation was measured from 10°to 80°elevation of WRT70.In the measurement of the Miyun 50 m radio telescope,the radio source Cyg A was used as the measurement beacon,and the effect of Cyg A broadening was removed (Zhang et al.2009).Because the beamwidth of WRT70 is only 130″,it is not suitable for the measurement of the 70 m antenna.We selected 3C 84 as the beacon,since it is almost the strongest radio point source in X-band.The Y-factor of 3C 84 in the X band of WRT70 exceeds 3 dB,so higher gain measurement accuracy can be obtained.

    Azimuth and elevation planes are scanned respectively.Normally the beam shape of large antennas at the X band can be presented well by a Gaussian function.Considering the influence of bottom noise and the change of ground noise,the constant term and the primary term are added in the antenna pattern fitting.So the fitting equation was taken as

    Here x is the distance to the radio source center,p1,p2,...,p5is the fitting parameters.Then the half and one-tenth beamwidth in the two planes can be measured by

    Figure 5 shows an example of the antenna pattern fitting of the scan in elevation direction of 3C 84 in the X-band.Figure 5(a) shows the raw observed pattern and fitting curve.Figure 5(b) shows the normalized observed pat-tern,fitting curve,and the residue.

    Figure 4.The second round pointing calibration with real-time atmospheric refraction.(a)The sky coverage of the pointing calibration observations,(b)the distribution of the pointing calibration residues,(c)pointing calibration residuals using or not using the real-time atmospheric refraction model,(d)residuals lower than certain elevation(rms),(e)number of positions in the 5°elevation range,(f)temperature distribution,(g)relative humidity distribution,(h)atmospheric pressure distribution.

    Figure 5.An example of the antenna pattern fitting of the scan in elevation direction of 3C 84 in the X-band,the continue line is the observed curve and the dashed line is the fitting curve.(a) The raw observed pattern and fitting curve,(b) normalized observed pattern and fitting curve,and the lower line is the residue.

    The antenna gain was then calculated by the following equation (Zhang et al.2009)

    Here,θ3dBAzand θ3dBElare the azimuth and elevation halfpower beamwidth,θ10dBAzand θ10dBElare the azimuth and elevation one-tenth power beamwidth,ΔGFis the feed loss,is the surface accuracy loss,λ is the wavelength and σ is the antenna panel accuracy.

    The method of the above gain variation measurement depends on the beams in the two perpendicular planes at different elevations.So the effect on gain measurements from atmosphere attenuation and system noise temperature error is very small.The observations were carried out in 2020 November,which were dry and clear days.

    The efficiency and gain results are shown in Figure 6.The efficiency is from 62.4% to 74.0%,and correspondingly the gain is 73.7–74.5 dBi.The maximum efficiency and gain are at 46°.9 elevation,which is close to the pre-adjusted position 48°elevation of the main surface.

    5.The System Noise Temperature

    The system temperature of a radio telescope comes from receiver noise,sky back-ground emission,ground emission,atmosphere emission,and antenna structure itself.For the case where both the polarizer and the field amplifier are cooled,system temperature can be express as the following equation

    Here Tris the equivalent temperature of the receiver at the LNA input port,TDis the Dewar temperature,T0is the normal environment temperature,Tskyis the equivalent temperature contributed by the sky back-ground emission,ground emission,atmosphere emission,and antenna structure at the feed output,Lpis the polarizer attenuation coefficient,and Lfis the feed attenuation coefficient without polarizer.

    Because the polarizer is cooled,the waveguide port is not easy to disassemble.The Y-factor measurement method by covering the normal temperature load and pointing to the cold air at elevation 90° can obtain high accuracy (Wang et al.2015).The measurement calculation equation is as follows

    We performed two sets of measurements.One set is a measure point at every 10°elevation intervals,and the other is a fast uniform rotation of the antenna from elevation 7°to 90°.The results are shown in Figure 7.The system temperature is from 48.8 to 25.8 K with the change of the elevation,and the two sets of measurements are basically the same.

    Figure 6.The X-band efficiency and gain of WRT70.(a) efficiency variation with elevation,(b) gain variation with elevation.

    Figure 7.The X-band system temperatures of WRT70.

    6.System Equivalent Flux Density

    System equivalent flux density (SEFD,Jy) is commonly used as the system sensitivity evaluation index in radio astronomy.The calculation formula can be expressed as(Wang et al.2015)

    Here k is the Boltzmann constant,η is antenna efficiency,D is the antenna aperture and Suis unit flow density(Jy).According to the calculation formula of antenna gain,SEFD can be also expressed as the function of the quality factor G/T

    Here λ is the wavelength.It can be seen that SEFD is determined if the antenna gain and system noise temperature are determined.By measuring the system noise temperature and antenna gain respectively,the measurement of SEFD can be realized.However,two measurements will increase error.So we used the on-off Y-factor method of radio source to measure SEFD directly.The major source of error in the measured SEFD is the uncertainty of the radio source flux density.The measured SEFD value can be expressed as

    Here S is the flow of point radio source(Jy),Y is the on-off Yfactor,K1is the atmospheric attenuation correction factor and K2is the radio source expansion correction factor.We used the radio source DR21 as the beacon signal,which is almost the strongest point source of known accurate flux density.The flux of DR21 in the X band has the form (Dent 1972)

    The measurement results of the X-band in right-hand-circular polarization are shown in Figure 8,which were measured on the evening of 2021 January,in dry and clear days.The minimum SEFD reached 25.8 Jy,less than 29.1 Jy above 30°elevations,and the maximum is 50.0 Jy at 7°.3 elevation.It can be seen from Figures 7 and 8 that the change of SEFD with elevation is mainly affected by the system noise temperature.

    Figure 8.The X-band SEFD of WRT70.

    7.Summary

    The 23 parameters pointing calibration model considering the atmospheric refraction correction in real time is presented and used in the telescope,and the pointing accuracy reachesin azimuth direction andin elevation direction respectively for different weather conditions.The elevation overall pointing accuracy increases fromand improves by 9.8%by using the real-time atmospheric refraction model.In particular,in the observation of low elevation angles below 20°,the pointing correction residual reduced by 19.3%,23.1% lower than 15°,and 27.5% lower than 10°.

    The efficiency is from 62.4%to 74.0%,and correspondingly the gain is 73.7 to 74.5 dBi variations with elevation from 10°to 80°elevation.The system temperature is from 48.8 to 25.8 K variations with elevation from 7° to 90° elevation.The minimum SEFD reached 25.8 Jy,less than 29.1 Jy above 30°elevation,and maximum 50.0 Jy at 7.3° elevation.It can be seen that the WRT70 has high performance,and has been successfully applied in the first Mars exploration mission of China.

    Acknowledgments

    Authors thank Mr.Fan Wang and other operators for their assistance on the observations.This work was funded by the Astronomical Joint Fund of the National Natural Science Foundation of China and Chinese Academy of Sciences(U1831114).

    岛国在线免费视频观看| 波多野结衣巨乳人妻| 午夜精品一区二区三区免费看| 色综合欧美亚洲国产小说| 亚洲第一欧美日韩一区二区三区| 可以免费在线观看a视频的电影网站| 精品国内亚洲2022精品成人| 一边摸一边做爽爽视频免费| 亚洲美女视频黄频| 亚洲精品中文字幕一二三四区| aaaaa片日本免费| 亚洲av中文字字幕乱码综合| 午夜免费成人在线视频| or卡值多少钱| 精品熟女少妇八av免费久了| 欧美久久黑人一区二区| 麻豆久久精品国产亚洲av| 婷婷六月久久综合丁香| 这个男人来自地球电影免费观看| 国产aⅴ精品一区二区三区波| 成年版毛片免费区| 男女做爰动态图高潮gif福利片| 亚洲精品中文字幕一二三四区| 久久久久久亚洲精品国产蜜桃av| 国产成+人综合+亚洲专区| 日韩欧美国产一区二区入口| 99国产综合亚洲精品| 三级毛片av免费| 可以在线观看毛片的网站| 不卡一级毛片| 国产一区二区激情短视频| 香蕉国产在线看| 国产人伦9x9x在线观看| 日本成人三级电影网站| 亚洲 欧美 日韩 在线 免费| 免费看十八禁软件| 天堂√8在线中文| 亚洲精品国产一区二区精华液| 国产区一区二久久| 国产成人一区二区三区免费视频网站| 可以免费在线观看a视频的电影网站| www国产在线视频色| 精品久久蜜臀av无| 中文字幕最新亚洲高清| 精品久久久久久久人妻蜜臀av| 超碰成人久久| 欧美三级亚洲精品| 国内少妇人妻偷人精品xxx网站 | 亚洲第一电影网av| 99国产精品99久久久久| 一二三四在线观看免费中文在| 国产男靠女视频免费网站| 长腿黑丝高跟| 18美女黄网站色大片免费观看| xxxwww97欧美| 国产又黄又爽又无遮挡在线| 国产亚洲精品综合一区在线观看 | 男女之事视频高清在线观看| 在线观看免费午夜福利视频| 国产精品久久电影中文字幕| 成人18禁高潮啪啪吃奶动态图| 欧美乱码精品一区二区三区| 少妇粗大呻吟视频| 男女做爰动态图高潮gif福利片| 国产伦人伦偷精品视频| 午夜精品在线福利| 日本三级黄在线观看| e午夜精品久久久久久久| 成人18禁在线播放| 精品久久久久久久人妻蜜臀av| 成人欧美大片| 国产午夜精品论理片| 精品不卡国产一区二区三区| 亚洲精品在线美女| 五月伊人婷婷丁香| 欧美午夜高清在线| 久久香蕉激情| 搡老熟女国产l中国老女人| 最好的美女福利视频网| 夜夜爽天天搞| 淫妇啪啪啪对白视频| 不卡一级毛片| 天天添夜夜摸| 欧美另类亚洲清纯唯美| 久久精品91蜜桃| 美女高潮喷水抽搐中文字幕| 黄色片一级片一级黄色片| 国产亚洲精品综合一区在线观看 | 老司机福利观看| 色老头精品视频在线观看| 久久精品国产亚洲av香蕉五月| 久久天躁狠狠躁夜夜2o2o| 日本三级黄在线观看| 久久草成人影院| 91九色精品人成在线观看| 女人高潮潮喷娇喘18禁视频| 久久亚洲真实| 成人18禁高潮啪啪吃奶动态图| 妹子高潮喷水视频| 午夜福利免费观看在线| 久久久久国产一级毛片高清牌| 国内精品一区二区在线观看| 一进一出好大好爽视频| 国产成人aa在线观看| 亚洲 欧美一区二区三区| 久久久久久人人人人人| 天天添夜夜摸| 高清在线国产一区| 国产一区二区三区视频了| 日韩三级视频一区二区三区| 国产亚洲欧美98| 欧美乱码精品一区二区三区| 熟女电影av网| 久久欧美精品欧美久久欧美| 国产精品av视频在线免费观看| 一二三四社区在线视频社区8| 日韩欧美三级三区| 国产私拍福利视频在线观看| 黄色女人牲交| 国产精品香港三级国产av潘金莲| 丰满的人妻完整版| 成人午夜高清在线视频| 午夜老司机福利片| 亚洲欧美精品综合久久99| 啦啦啦观看免费观看视频高清| 成人特级黄色片久久久久久久| 1024香蕉在线观看| 成年免费大片在线观看| 变态另类成人亚洲欧美熟女| 亚洲av美国av| 亚洲人成电影免费在线| 亚洲专区字幕在线| 国产精品免费一区二区三区在线| 久99久视频精品免费| 欧美日本视频| 丰满的人妻完整版| 欧美色视频一区免费| 精品久久久久久久久久久久久| 国产免费男女视频| 12—13女人毛片做爰片一| 又紧又爽又黄一区二区| 一级a爱片免费观看的视频| 亚洲男人天堂网一区| 亚洲国产精品成人综合色| 操出白浆在线播放| 一卡2卡三卡四卡精品乱码亚洲| 脱女人内裤的视频| 女人爽到高潮嗷嗷叫在线视频| 成人av一区二区三区在线看| 精品一区二区三区视频在线观看免费| 午夜免费成人在线视频| 日韩高清综合在线| 欧美乱码精品一区二区三区| 欧美午夜高清在线| 两人在一起打扑克的视频| 亚洲最大成人中文| 日韩三级视频一区二区三区| 少妇的丰满在线观看| 国产精品98久久久久久宅男小说| 亚洲av电影在线进入| 成人三级做爰电影| 伊人久久大香线蕉亚洲五| 两个人视频免费观看高清| www.熟女人妻精品国产| 久久久久久九九精品二区国产 | 亚洲va日本ⅴa欧美va伊人久久| 久久精品国产清高在天天线| 两人在一起打扑克的视频| 国产精品99久久99久久久不卡| 亚洲精品美女久久久久99蜜臀| 在线看三级毛片| 777久久人妻少妇嫩草av网站| 美女高潮喷水抽搐中文字幕| 亚洲欧美激情综合另类| 午夜激情av网站| 长腿黑丝高跟| svipshipincom国产片| 亚洲精品中文字幕在线视频| 日日干狠狠操夜夜爽| 国产人伦9x9x在线观看| 又粗又爽又猛毛片免费看| 成人午夜高清在线视频| 色综合婷婷激情| 精品久久久久久,| 一本久久中文字幕| 少妇熟女aⅴ在线视频| 欧美黄色片欧美黄色片| 国产成人欧美在线观看| 一进一出抽搐gif免费好疼| 搡老妇女老女人老熟妇| 午夜精品一区二区三区免费看| 久久人人精品亚洲av| 免费看十八禁软件| 不卡av一区二区三区| 免费在线观看完整版高清| 午夜激情福利司机影院| 在线观看www视频免费| 舔av片在线| 国产精品久久电影中文字幕| 人人妻人人澡欧美一区二区| 1024手机看黄色片| 国内精品一区二区在线观看| 精品久久久久久久毛片微露脸| 亚洲精品国产一区二区精华液| 午夜视频精品福利| 久久精品91无色码中文字幕| 国产午夜精品论理片| 宅男免费午夜| 欧美中文日本在线观看视频| 亚洲一区二区三区色噜噜| 淫秽高清视频在线观看| 国产精品永久免费网站| 欧美又色又爽又黄视频| 丝袜人妻中文字幕| 老司机深夜福利视频在线观看| 日韩精品免费视频一区二区三区| 日韩高清综合在线| 欧美最黄视频在线播放免费| 欧美在线黄色| 久久香蕉激情| 国语自产精品视频在线第100页| 亚洲av美国av| 99久久综合精品五月天人人| 亚洲成人久久性| 亚洲九九香蕉| 91老司机精品| 亚洲精品一区av在线观看| 男女视频在线观看网站免费 | 88av欧美| 免费看日本二区| 麻豆成人av在线观看| 国产一区二区激情短视频| 99久久综合精品五月天人人| 欧美日韩精品网址| 国产精品av视频在线免费观看| 亚洲男人天堂网一区| 精品免费久久久久久久清纯| 制服诱惑二区| 国产精品av久久久久免费| 亚洲精品久久国产高清桃花| 国产精品久久电影中文字幕| 欧美国产日韩亚洲一区| av视频在线观看入口| 欧美成人免费av一区二区三区| 黑人操中国人逼视频| 亚洲精品一区av在线观看| 91成年电影在线观看| 国产单亲对白刺激| 欧美成人午夜精品| 久久久久久久久久黄片| 禁无遮挡网站| 窝窝影院91人妻| 国产精品影院久久| 国产伦人伦偷精品视频| 亚洲avbb在线观看| 亚洲专区国产一区二区| 亚洲成人国产一区在线观看| 精品欧美国产一区二区三| 国产精品电影一区二区三区| 男人舔奶头视频| 久久中文字幕人妻熟女| 中出人妻视频一区二区| 欧美精品啪啪一区二区三区| 草草在线视频免费看| 久久久久九九精品影院| 黄色a级毛片大全视频| 国产一区二区三区视频了| 国产一区二区激情短视频| 欧美一级毛片孕妇| 熟女电影av网| 亚洲av成人av| 午夜老司机福利片| 欧美 亚洲 国产 日韩一| 久久精品夜夜夜夜夜久久蜜豆 | 99久久精品国产亚洲精品| 国产免费男女视频| 男女做爰动态图高潮gif福利片| 午夜免费激情av| 欧美最黄视频在线播放免费| 99国产极品粉嫩在线观看| 亚洲精品中文字幕一二三四区| www.熟女人妻精品国产| 身体一侧抽搐| 国产亚洲精品久久久久5区| 首页视频小说图片口味搜索| 久久久久国产精品人妻aⅴ院| 久久人人精品亚洲av| 亚洲国产精品成人综合色| 国产高清videossex| 小说图片视频综合网站| 国产欧美日韩一区二区三| 欧美成人一区二区免费高清观看 | 观看免费一级毛片| 成人av一区二区三区在线看| 国产亚洲欧美在线一区二区| 国产日本99.免费观看| 久久久久国产精品人妻aⅴ院| 国产精品久久久人人做人人爽| 高清毛片免费观看视频网站| 日日摸夜夜添夜夜添小说| 丁香欧美五月| 日韩大码丰满熟妇| 精品福利观看| 亚洲国产精品成人综合色| 久久久久久亚洲精品国产蜜桃av| 悠悠久久av| 18禁裸乳无遮挡免费网站照片| 日本熟妇午夜| 狠狠狠狠99中文字幕| 久久久久免费精品人妻一区二区| 黄色视频,在线免费观看| 国产精品电影一区二区三区| 黄色视频不卡| 脱女人内裤的视频| 99在线人妻在线中文字幕| 国产aⅴ精品一区二区三区波| 亚洲片人在线观看| 欧美中文综合在线视频| 国产高清激情床上av| 亚洲第一欧美日韩一区二区三区| 色在线成人网| 免费看日本二区| 人人妻人人澡欧美一区二区| 亚洲免费av在线视频| 大型黄色视频在线免费观看| 欧美日韩一级在线毛片| 国产高清视频在线播放一区| 精品熟女少妇八av免费久了| 男人舔女人的私密视频| 亚洲成人久久爱视频| 全区人妻精品视频| 桃红色精品国产亚洲av| 国产午夜福利久久久久久| 色精品久久人妻99蜜桃| 亚洲精品一卡2卡三卡4卡5卡| 久久精品国产亚洲av高清一级| 国产亚洲精品av在线| 国产午夜精品久久久久久| 欧美久久黑人一区二区| 亚洲国产欧洲综合997久久,| 一级毛片女人18水好多| 久久精品夜夜夜夜夜久久蜜豆 | 精品一区二区三区四区五区乱码| 夜夜夜夜夜久久久久| 在线观看www视频免费| 久久香蕉国产精品| 国产av一区在线观看免费| 国产精品久久久久久亚洲av鲁大| a在线观看视频网站| 成人永久免费在线观看视频| 99热只有精品国产| 成年女人毛片免费观看观看9| 天天躁夜夜躁狠狠躁躁| 男插女下体视频免费在线播放| 日本熟妇午夜| 90打野战视频偷拍视频| 成人18禁高潮啪啪吃奶动态图| 国产男靠女视频免费网站| 久久久久九九精品影院| 超碰成人久久| 99在线人妻在线中文字幕| 亚洲欧美精品综合久久99| 亚洲欧美日韩无卡精品| 免费看日本二区| 国产成人系列免费观看| 精品午夜福利视频在线观看一区| 亚洲精品久久成人aⅴ小说| 午夜福利成人在线免费观看| 久久精品国产亚洲av香蕉五月| 日韩欧美在线二视频| 久久这里只有精品中国| 亚洲欧洲精品一区二区精品久久久| 99久久精品热视频| xxx96com| 麻豆国产97在线/欧美 | 人人妻人人澡欧美一区二区| 精品欧美一区二区三区在线| 男女下面进入的视频免费午夜| 99热这里只有精品一区 | 亚洲精品美女久久久久99蜜臀| 夜夜躁狠狠躁天天躁| 一二三四社区在线视频社区8| 免费观看精品视频网站| 欧美日韩中文字幕国产精品一区二区三区| 精品国产美女av久久久久小说| 欧美一级毛片孕妇| 国产主播在线观看一区二区| 天堂影院成人在线观看| 十八禁网站免费在线| 日韩大尺度精品在线看网址| 国产一区二区激情短视频| 成人精品一区二区免费| 久久中文字幕人妻熟女| 成人国产综合亚洲| 免费看十八禁软件| 久久婷婷成人综合色麻豆| 99精品久久久久人妻精品| 成人av一区二区三区在线看| 久久久久性生活片| 一级a爱片免费观看的视频| 色av中文字幕| 18禁黄网站禁片免费观看直播| 亚洲av电影不卡..在线观看| 午夜免费观看网址| 国产99久久九九免费精品| 国产午夜精品久久久久久| 中文字幕av在线有码专区| av在线天堂中文字幕| 不卡av一区二区三区| 亚洲黑人精品在线| 亚洲 欧美 日韩 在线 免费| 日本一区二区免费在线视频| 亚洲欧美日韩无卡精品| 精品久久久久久成人av| 香蕉av资源在线| xxx96com| 亚洲av日韩精品久久久久久密| 亚洲国产日韩欧美精品在线观看 | videosex国产| 两性午夜刺激爽爽歪歪视频在线观看 | 中亚洲国语对白在线视频| 婷婷精品国产亚洲av| 高清毛片免费观看视频网站| 日本免费一区二区三区高清不卡| 热99re8久久精品国产| 午夜老司机福利片| 成人国产一区最新在线观看| 免费看日本二区| 男人舔女人下体高潮全视频| 欧美久久黑人一区二区| 久久精品综合一区二区三区| 国产精品亚洲美女久久久| 69av精品久久久久久| 国语自产精品视频在线第100页| 18禁裸乳无遮挡免费网站照片| 亚洲一区二区三区色噜噜| 狂野欧美白嫩少妇大欣赏| 无遮挡黄片免费观看| 一卡2卡三卡四卡精品乱码亚洲| 国产欧美日韩一区二区三| 午夜精品一区二区三区免费看| 国产精品日韩av在线免费观看| 国产亚洲精品av在线| 最近最新中文字幕大全免费视频| 欧美成人午夜精品| bbb黄色大片| 国产亚洲精品第一综合不卡| √禁漫天堂资源中文www| 宅男免费午夜| 波多野结衣巨乳人妻| 黄频高清免费视频| 三级男女做爰猛烈吃奶摸视频| 午夜精品在线福利| 日韩精品青青久久久久久| 亚洲成人久久爱视频| 国产精品乱码一区二三区的特点| 亚洲精品在线美女| 久久久久久久精品吃奶| 国产午夜精品论理片| 黄色视频不卡| 午夜久久久久精精品| 精华霜和精华液先用哪个| 一进一出抽搐动态| 亚洲欧美日韩东京热| 亚洲精品久久成人aⅴ小说| 久久久久国内视频| 久久久精品大字幕| 亚洲欧美日韩无卡精品| 最近在线观看免费完整版| 搞女人的毛片| 亚洲精品美女久久久久99蜜臀| 国产精品 欧美亚洲| 悠悠久久av| 午夜福利视频1000在线观看| 亚洲自偷自拍图片 自拍| 男人舔女人下体高潮全视频| av片东京热男人的天堂| 嫩草影视91久久| 麻豆国产97在线/欧美 | 成人18禁在线播放| 亚洲成人中文字幕在线播放| 亚洲免费av在线视频| 久久婷婷成人综合色麻豆| 在线观看日韩欧美| 国产高清视频在线播放一区| 无人区码免费观看不卡| 国产精品影院久久| 国产成人欧美在线观看| 2021天堂中文幕一二区在线观| 欧美人与性动交α欧美精品济南到| 日本黄色视频三级网站网址| 18美女黄网站色大片免费观看| 婷婷精品国产亚洲av| 麻豆成人午夜福利视频| 日韩成人在线观看一区二区三区| 中文在线观看免费www的网站 | 亚洲国产精品合色在线| 精品欧美国产一区二区三| 久久人人精品亚洲av| 天堂影院成人在线观看| 中文字幕熟女人妻在线| 757午夜福利合集在线观看| 国产精华一区二区三区| 99国产极品粉嫩在线观看| 国产久久久一区二区三区| 一二三四社区在线视频社区8| 成人永久免费在线观看视频| 桃红色精品国产亚洲av| 亚洲乱码一区二区免费版| 国产片内射在线| 久久香蕉国产精品| 手机成人av网站| 国产av又大| 亚洲国产精品999在线| 成人18禁高潮啪啪吃奶动态图| 亚洲狠狠婷婷综合久久图片| 国产精品免费视频内射| 中文字幕精品亚洲无线码一区| 国产主播在线观看一区二区| 国产激情偷乱视频一区二区| 99久久无色码亚洲精品果冻| 国产爱豆传媒在线观看 | 免费高清视频大片| 亚洲自偷自拍图片 自拍| 国产精华一区二区三区| 日韩成人在线观看一区二区三区| 黑人欧美特级aaaaaa片| 日韩欧美免费精品| 日本黄色视频三级网站网址| 欧美3d第一页| 天堂动漫精品| www日本在线高清视频| 1024手机看黄色片| 99re在线观看精品视频| 成人欧美大片| 日本 欧美在线| 97碰自拍视频| 中文字幕最新亚洲高清| 欧美日韩亚洲综合一区二区三区_| 熟女电影av网| 久久久国产成人精品二区| 深夜精品福利| 在线十欧美十亚洲十日本专区| 制服诱惑二区| 国产不卡一卡二| 国产激情欧美一区二区| 欧美在线一区亚洲| 嫩草影视91久久| 精品久久久久久久久久久久久| 嫩草影院精品99| 精品少妇一区二区三区视频日本电影| 亚洲人成伊人成综合网2020| 免费电影在线观看免费观看| 欧美日韩国产亚洲二区| 亚洲色图 男人天堂 中文字幕| 欧美在线黄色| 精品久久久久久久人妻蜜臀av| 久久精品91无色码中文字幕| 2021天堂中文幕一二区在线观| 两个人看的免费小视频| 亚洲精品中文字幕一二三四区| 亚洲成人久久性| 亚洲国产中文字幕在线视频| 日韩三级视频一区二区三区| av免费在线观看网站| aaaaa片日本免费| 免费在线观看完整版高清| 不卡av一区二区三区| 欧美人与性动交α欧美精品济南到| 男女下面进入的视频免费午夜| 波多野结衣巨乳人妻| 日韩高清综合在线| 成人三级黄色视频| 日日摸夜夜添夜夜添小说| 国产精品久久视频播放| 国产野战对白在线观看| 精品国产乱子伦一区二区三区| 岛国在线免费视频观看| 欧美+亚洲+日韩+国产| 岛国在线观看网站| 别揉我奶头~嗯~啊~动态视频| 久久这里只有精品中国| 国产熟女xx| 国产欧美日韩一区二区精品| 久久精品aⅴ一区二区三区四区| 亚洲av成人一区二区三| 毛片女人毛片| 一区福利在线观看| 在线观看舔阴道视频| 禁无遮挡网站| 国产精品香港三级国产av潘金莲| 午夜日韩欧美国产| 久久九九热精品免费| 国产精品香港三级国产av潘金莲| 亚洲av成人一区二区三| 制服诱惑二区| 国产精品av视频在线免费观看| 香蕉久久夜色| 别揉我奶头~嗯~啊~动态视频| 国产精品久久久久久精品电影| 国产精品亚洲美女久久久| 在线观看日韩欧美| 1024香蕉在线观看| 岛国在线观看网站| aaaaa片日本免费| 亚洲熟妇中文字幕五十中出| 欧美精品亚洲一区二区| 欧美黄色淫秽网站| 亚洲国产欧美人成| 日本黄色视频三级网站网址| 亚洲avbb在线观看| 久久中文字幕人妻熟女| 老熟妇仑乱视频hdxx| 亚洲天堂国产精品一区在线| av视频在线观看入口|