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

    The Simultaneous Three-channel Multicolor CCD Photometric System of the 1.2 m Telescope at Jilin Astronomical Observatory*

    2022-05-24 06:34:04BingLiNiuChengZhiLiuZhenWeiLiZheKangandYouLv

    Bing-Li Niu ,Cheng-Zhi Liu ,Zhen-Wei Li ,Zhe Kang ,and You Lv

    1 Changchun Observatory,National Astronomical Observatories,University of Chinese Academy of Sciences,Changchun 130117,China;lcz@cho.ac.cn

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

    Abstract We will introduce the construction and design of a new simultaneous three-channel multicolor CCD photometer.This photometer has been mounted on the 1.2 m telescope at Jilin Astronomical Observatory,and is applied to study space debris,gamma-ray burst afterglows,asteroids and other rapidly variable objects.As one of the ground follow-up telescopes of the Chinese-French Space-based multi-band astronomical Variable Objects Monitor(SVOM) mission,it appears very essential to evaluate the performance of the photometry system of this photometer,which can achieve simultaneous imaging within a field of view ofandin the Sloan Digital Sky Survey (SDSS) g′,r′,i′ bands,respectively.Photometric calibrations were carried out by using plenty of SDSS standard stars,and the relationship between the photometric system and the Johnson-Bessel filter system was also studied.The results of the performance evaluation from observing open cluster M67 are presented.

    Key words:astronomical instrumentation– photometer– photoelectric photometry

    1.Introduction

    In recent years,multicolor CCD photometric technology has attracted great attention for its wide applications in astronomical observations including variable stars,active galactic nuclei and gamma-ray burst (GRB) afterglows (e.g.,Greiner et al.2008),so that the underlying mechanism and physical origin of the observed properties could be studied (Warner 2006).In addition,photometric measurements of fast-moving objects can determine their surface materials,shapes,rotations and internal structure.It can also be employed to monitor the operating state of space debris (Santoni et al.2013;Sun &Zhao 2013).Currently,the most used photometry system is based on the filter wheel,which changes the filter in turn to achieve photometry at different wave bands.However,it will introduce incorrect results when observing rapid variables and fast moving objects due to the rapid variation of brightness or color.In comparison with the traditional multicolor CCD photometry system,three-channel multicolor CCD photometry has advantages in making observations simultaneously in different filter bands,as well as with higher temporal resolution.This approach is especially suitable for measurements of shortperiod variables,GRBs (Greiner et al.2008) and space debris.

    Nowadays,many countries have carried out simultaneous multicolor CCD photometry with different telescopes,such as the Chilean GROND-a 7-Channel imager which is installed on the MPI/ESO 2.2 m telescope (Greiner et al.2008) and ULTRACAM,which is a portable and high-speed imaging photometer(Dhillon et al.2007)that has been attached to both the 4.2 m William Herschel Telescope on La Palma and the 8.2 m Very Large Telescope in Chile.The Japanese Multicolor Imaging Telescopes for Survey and Monstrous Explosions(MITSuME) have been equipped on two 50 cm telescopes at Akeno and Okayama,respectively (Kotani et al.2005).However,in China there has been little research on the application of a multichannel photometer(Mao et al.2013).We always have strong motivation to do simultaneous multicolor CCD photometry for various variable objects in our observatory to promote the capability of astronomical observation in our country.In our work,a simultaneous three-channel multicolor CCD photometry system was designed and constructed utilizing the Sloan Digital Sky Survey (SDSS) g′,r′,and i′ bands,and was installed on the Cassegrain focus of the 1.2 m large-aperture telescope at Changchun Observatory,National Astronomical Observatories,Chinese Academy of Sciences,China.The 1.2 m telescope is also a member of the ground-based telescope in the Chinese-French Space-based multi-band astronomical Variable Objects Monitor (SVOM)mission,for GRB studies (Paul et al.2011),so it is of great significance to study the performance of this equipment.

    We have characterized the photometry system of the photometer and results are reported in this paper.The telescope and site conditions are briefly introduced in Section 2.In Section 3,we describe the design and structure of the threechannel CCD photometer.The calibration and performance of the photometry system are reported in Section 4 and Section 5,respectively.Finally,we give our conclusion in Section 6.

    2.The Telescope Parameters and Site Conditions

    2.1.The Telescope Parameters

    The 1.2 m telescope that belongs to Changchun Observatory of NAOC was designed and made by Shanghai Astronomical Observatory.A photograph and the optical path of the 1.2 m telescope are shown in Figures 1 and 2,respectively.The 1.2 m telescope has a reflective optical structure and two optical focus systems,one is a prime focus system and the other is a Cassegrain focus system.The switch between these two focus systems is realized by the secondary mirror controlled by a switching mechanism.The simultaneous three-channel multicolor CCD photometry system was installed in the Cassegrain focus system,so we will introduce the Cassegrain focus system in detail,and the performance of the prime focus system will be described in another paper.Both mirrors and mount characteristics of the 1.2 m telescope are listed in Table 1.The Cassegrain focus system was designed with a focal ratio of 8,and both primary and secondary mirrors of the 1.2 m telescope are hyperboloidal mirrors,with effective diameters of 1180 mm and 340mm,respectively,forming a spherical and coma aberration free Ritchey-Chrétien (R-C) optical system.Over the entire field of view (FOV),80% of the energy is concentrated within a 1″ point-spread function,and the transmittance of the optical system is better than 80%.The mount system is an equatorial fork mount with high accuracy friction servo motors.

    Figure 1.The 1.2 m telescope.

    Figure 2.The optical path of the telescope.

    Figure 3.Monthly observational time in each month during the past four years.Left panel:numbers of monthly observational nights;Right panel:numbers of monthly observational hours.

    2.2.The Site Conditions

    The 1.2 m telescope is located at Xiaosuihe Village in Jilin City (126°20′24″ East,43°49′48″ North),with an elevation of about 313 m.We measured the seeing of Jilin Astronomical Observatory by using a differential image motion monitor(DIMM) with the help of Yunnan Observatory.The average value of seeing is aboutand the best value is better thanThe number of clear nights was achieved by using data recorded during the past four years from January 2017 to December 2020 as depicted in Figure 3.The observational time of the past four years was obtained according to the observation log of our Mini Optoelectronic Telescope Array.The left panel of Figure 3 displays the number of nights when observations were taken in each month.Meanwhile,the number of hours when the telescope was being operated during each month is shown in the right panel of Figure 3.As can be seen from Figure 3,during the summer time from June to August,the number of observation nights was minimum each year due to rainy and cloudy weather.In winter,it is necessary to defrost and maintain the dome due to extremely low temperatures,so actually there are more clear nights from November to January of the following year.As the plots demonstrate,during the past four years the number of observation nights in each year ranges from 165 to 223 with an average of 189 nights,and the number of observation hours per year ranges from 1599 hr to 2179 hr with an average of 1856 hr.

    Figure 4.Photograph of the three-channel CCD photometer.

    3.Parameters of the Three-Channel CCD Photometer

    There are many ways to achieve multi-channel imaging simultaneously,such as by using a Philips prism (Berquist &Berquist 1996) or spectroscopic lens.A Simultaneous Photometer for Infrared and Visual light (SPIV),which was constructed by the Steward Observatory Mirror Lab of the University of Arizona,adopts the method of using a spectroscopic lens to achieve simultaneous photometry (Eisendart 1984).There are many advantages of this method,such as simple structure,easy assembly,lower polarization effect and accuracy of components that are easy to guarantee.In addition,it is easier to make the device with a larger size.In order to obtain a larger FOV of about 21′×21′,we adopt a spectroscopic lens to design the simultaneous three-channel multicolor CCD photometer.

    3.1.Physical Structure

    The structure of a three-channel CCD photometer is depicted in Figure 4.

    Figure 5.Optical structure of the three-channel photometer.

    Figure 6.Transmission of the u′ g′ r′ i′ z′ filters (Fukugita et al.1996).

    1.Reflector:this part is composed of three mirrors including M3,M4 and M5 to ensure the alignment between the photometer and the main optical path because of the long focal length of the Cassegrain focus system.Each mirror has three degrees of freedom adjustment,which are pitch,azimuth and optical axis translation for the convenience of installation and alignment.

    2.Field lens:this mirror is located between M5 and the first dichroic lens to collimate the light reflected by M5.

    3.Beam-splitting lens:as the core of the three-channel CCD photometer,this part consists of two dichroic lenses that require special coating (Eisendart 1984) to achieve separation of the three color beams.

    4.Imaging lens group:in order to optimize the beam quality,three lens groups were designed to correct the light path for three channels,respectively.

    5.Filters:the light in the specified band is obtained after filtering by the filter,which is used for imaging in three channels.

    6.CCD detectors:each channel is equipped with the same CCD camera,which is an Andor iKon-L DZ936N-BV,as an imaging terminal.

    3.2.Optical System

    The optical structure of the three-channel CCD photometer is illustrated in Figure 5.The optical wave bands chosen for the three channels are the g′,r′ and i′ bands of the SDSS system(Figure 6),respectively.The bandpasses of the SDSS filter system are wider than those of the Johnson/Bessel filter system,which is more suitable for faint object detection(Fukugita et al.1996).

    To obtain a larger FOV,the three bands of g′,r′ and i′ are separated by using a beam-splitting lens,which divides the light collimated by the field lens into three beams for each channel.The first dichroic lens reflects the g′band beam,while the r′ band beam and the i′ band beam pass through,and the second dichroic lens reflects the r′ band beam and transmits the i′band beam.After filtering by the corresponding filter,the beams of each band enter the CCD detector to complete the operation of photometry and photometric measurement.Many factors are taken into account to evaluate the image quality,including FOV,degradations,full width at half maximum(FWHM),ellipticity,Moffat β and asymmetry.The results are shown in Figures 7,8 and Table 2.

    3.3.CCD Detectors

    The three-channel CCD photometer is equipped with three uniform detectors,which are backside illuminated.Parameters of the detector including readout noise,pixel size,quantum efficiency (QE),dark current and chip scale are under consideration,and the specifications of the CCD camera are listed in Table 3.Due to differences on the chip,the stability of bias,dark current and linearity of the CCD camera in each channel was tested,as displayed in Figures 9,10 and 11,respectively.In addition,the test results are plotted in Table 4.The QE curve of iKon-L DZ936N-BV is depicted in Figure 12,and the specific parameters of each channel are listed in Table 5.

    Table 5 Parameters of three Channels

    Table 6 Atmospheric Extinction Coefficients and Sky Brightness at Jilin Astronomical Observatory

    Figure 11.The linear response curve of the iKon-L DZ936N-BV with scan rate of 1 MHz and gain 1×.

    Figure 12.QE curve of iKon-L DZ936N-BV.

    Table 4 Test Results of the CCD Detectors

    Figure 9.The stability of 100 bias frames from the iKon-L DZ936N-BV at a temperature of -70°C,with scan rate of 1 MHz.

    Figure 10.The dark current of iKon-L DZ936N-BV varies with different temperature with scan rate of 1 MHz and gain 1×.

    Table 2 The Results of Image Quality Evaluation

    Table 3 Specifications of the CCD Camera

    Figure 7.Original images of open cluster M67 in the g′,r′and i′bands.Top left panel:the original image of the g′channel,with an FOV of 215×215;top right panel:the original image of the r′channel,with an FOV of 215×215;bottom panel:the original image of the i′channel,with an FOV of 213×213.The design of the optical path leads to the FOV of i′ band being different from the other two bands.

    Figure 8.Degradations in optical quality along the off angle direction in the three channels.The result shows that the degradations in optical quality caused by aberrations and the possible CCD tilt are within 9%,which would be a side effect of large FOV.The star distortion at the edge of the FOV is inevitable,so we pay more attention to the image quality at the center of the FOV.

    3.4.Software

    The software system of the three-channel CCD photometer consists of three parts:observation support system,data acquisition system and data processing system.Each system communicates with each other through a socket based on TCP/IP(Ge et al.2015).The three cameras of each channel are controlled by three industrial personal computers,which can independently control the camera in this band for acquisition.The whole system runs in the Linux environment,and the three-channel camera is controlled by command for acquisition.In the Windows environment,the simultaneous acquisition of three-channel CCD cameras can be carried out manually.

    1.The observation support system mainly provides a unified and high-precision time reference to the data acquisition and processing system through the time latch.The telescope control software and CCD camera control software adopt the on-demand synchronization mode,while the observation control service software and data processing software adopt the timing synchronization mode to realize the simultaneous observation of the three channel cameras.

    2.The data acquisition system includes trigger source,client,control server and the computer that controls the telescope,as well as three industrial personal computers that control the CCD camera.The software that controls the CCD camera in each channel runs on the Linux platform,and integrates the terminal racer,focusing and auto focusing functions.The system can realize tracking observation and data acquisition flow of space fragments,conventional astronomical targets and emergencies.

    3.The data processing system runs on the Linux platform,and the core algorithm of photometry is based on astronomical tools including Source Extractor (SExtractor,Bertin &Arnouts 1996),Astrometry-net (Lang et al.2010) and Astronomy in Python (AstroPy,Astropy Collaboration et al.2013).

    4.System Calibration

    The three-channel CCD photometer was installed on the 1.2 m telescope on 2020 September 26.In order to checkwhether the optical path of the three-channel CCD photometer was smooth or not after assembly,laser beams were used to test the overall optical path.The actual imaging operation was then performed and the mechanical focusing was carried out to ensure the imaging was displayed in the center of the FOV in each channel.

    The system was calibrated during four observing runs taken in December of 2020,and March,April and November of 2021,respectively.This work was achieved by observing plenty of SDSS standard stars and Landolt standard stars which vary in magnitude and color,covering a wide sky region.The observation night was a standard photometric night,which was clear and moonless with a stable atmosphere.The calibration process includes the determination of the photometric zero-points for the three bands,system efficiency,field curvature and parfocality.

    4.1.System Efficiency

    The three-channel CCD photometer is an independent product with adjustment devices that can ensure the photometer is directly mounted on the telescope for imaging.In order to guarantee highefficiency of the three-channel CCD photometer,the system efficiency is taken into account in the design of its optical system,overall structure layout and electronic control communication system,so as to optimize the position of each mirror.In order to maximize the ability of the three-channel CCD photometer,each channel of the photometer is equipped with the same type of CCD camera as the imaging terminal.

    The method of calculating system efficiency has been introduced in detail by Mao et al.(2013),which is adopted in our work to evaluate the efficiency of the three-channel CCDphotometer.By observing BD+750325 (airmass=1.2)during the system’s calibration,we calculate the efficiency of the three-channel CCD photometer.The result shows that the efficiencies of the g′,r′,i′ channels are 65.6%,68.3% and 63.7%,respectively.Considering the machining error of components and the overall assembly error of the three-channel CCD photometer,the efficiency fully meets the expectation.

    4.2.Limiting Magnitudes

    By observing Vega (α Lyr) at different airmasses,we determined the photometric zero-points of each band.The zeropoints of the g′,r′ and i′ bands are 23.49,22.85 and 21.97,respectively,as affirmed in Figure 13.Meanwhile,the atmospheric extinction coefficients in the g′,r′and i′bands together with night sky brightness (Zhang et al.2013) at Jilin Astronomical Observatory were measured,as listed in Table 6.Based on the measured zero-points,signal-to-noise ratio(SNR)and extinction coefficient,we calculate the limiting magnitudes of the g′,r′and i′ bands,which were 15.26,16.39,15.63,respectively (1 s exposure time and SNR=5).

    Table 7 The Parameters of M67 member stars observed on 2021 November 27

    Table 8 Observing Log of M67

    Figure 13.Zero-points in the g′,r′and i′bands.The photometric zero-points of each band are based on flux from a 1 s exposure of Vega (α Lyr) at different airmasses.Top panel:g′ channel;middle panel:r′ channel;bottom panel:i′channel.

    Figure 14.The relationship between standard magnitude (x) and transformed magnitude (y) is calculated by using transformation equations.Top left panel:g′channel;top right panel:r′ channel;bottom panel:i′ channel.

    4.3.Photometric System Transformations

    To make the system transformation results more reliable,we observed various types of SDSS standard stars.Meanwhile,a wide range of airmass and high SNRs were considered.The calibrated observations were carried out under standard photometric nights on 2020 December 6,2021 March 24,2021 April 26 and 2021 November 25.A number of standard stars were observed,and the position,magnitude and color index of the objects were taken into account.The transformation accuracies from the three-channel photometer system (T)to the SDSS standard system(S)are 2.9%in g′band,2.4%in r′band and 3.1% in i′ band,which are indicated by Figure 14,and the transformation equations are listed as follows(November 25):

    Figure 15.The photometric errors of standard stars with different magnitudes(exposure time=10 s,airmass=1.131).Top left panel:g′channel;top right panel:r′channel;bottom panel:i′ channel.

    Figure 16.Differential photometric results of M67.Top panel:light curves of EX Cnc in three bands;bottom panel:light curves of the comparison star relative to the check star.There is a constant shift in the y-direction in order to optimize the visibility of light curves in the two panels.

    Figure 17.The distribution of photometric error of M67 member stars.

    The relationship between the three-channel photometer system(T)and the Johnson-Bessel filter system (J)are studied by observing a number of Landolt standard stars (Landolt 1992,2013).Transformation equations are listed as follows(November 25):

    The photometric accuracy was evaluated by observing abundant standard stars in SDSS and clusters (M67,Czernik 43).The photometric precision of the system in 9–14 mag is better than 4.3%in g′band and r′band,and 4.7%in i′band,as shown in Figure 15.We also calculated the primary and secondary extinction coefficients,and the sky brightness of Jilin Astronomical Observatory,as displayed in Table 6.The determination of the atmospheric extinction coefficient and the sky brightness takes a long time.Similarly,the system transformation from the three-channel photometer system to SDSS filter system is also a long-term work.In the future,we will improve the method for more comprehensive and systematic measurement,and the results will be described in another paper.

    5.Performance of the Three Channel CCD Photometer

    The simultaneous three-channel photometer is operated in Jilin Astronomical Observatory at present.As one of the ground-assisting telescopes of the SVOM project,it is of great significance to evaluate the performance of the photometer.In this paper,M67 (also known as Messier 67 or NGC 2682,hereinafter referred to as M67) is selected as the observation object,which is an open star cluster located in the constellation of Cancer,with about 500 members (Zhang et al.2005;Pribulla et al.2008;Peng et al.2016;Canton et al.2021).The star field is not dense and the brightness is concentrated in 10–11 mag,which is an ideal photometric calibration area.Our photometric work was performed in this cluster using a photometric standard catalog provided by the American Association of Variable Star Observers (AAVSO) [https://www.aavso.org].

    EX Cnc is a blue straggler with δ Scuti-like variability in the solar-age Galactic open cluster M67 (Zhang et al.2005;Yakut et al.2009),whose period of light variation is 0.04852 days.The variable stars in the open cluster M67 were measured by using the method of differential photometry,according to the following principles:the comparison star and variable star are close in position,brightness and color,so that six standard stars in the cluster were selected as comparisons for eliminating the effects of atmospheric extinction.In addition,the calibration stars were calibrated by six check stars.We conducted overnight tracking observations of the open clusters,and detailed information about the observed objects and observing log of M67 are provided in Tables 7 and 8,respectively.

    The magnitudes of the comparison stars are stable as affirmed in Figure 16,thus it is feasible to use these six stars to evaluate the photometric accuracy and differential magnitudes of EX Cnc.In addition,by calculating the rms of the differential magnitudes between the comparison stars and check stars,we find that the photometric accuracies of the g′,r′and i′ bands are 1.59%,1.48% and 1.62%,respectively.In order to study the distribution of the photometric error,135 member stars of M67 ranging from 10 mag to 16 mag were selected,as depicted in Figure 17.

    The photometric error of most M67 member stars is better than 0.1 mag,indicating that the three-channel CCD photometer has good stability and reliability.There is good coordination and consistency in the work,as a result of the same type of camera.The photometric error of faint stars is relatively large due to the limiting magnitude of each channel,which could be improved by appropriately increasing the exposure time while ensuring a higher SNR.The variation of photometric error with the magnitude of the standard star can be suggested by Figure 15.The results fully meet the expectations and work demands,considering the inevitable machining and assembly errors of the parts of the three-channel CCD photometer.

    6.Summary

    The multicolor CCD photometer installed in Jilin Astronomical Observatory is a three-channel simultaneous photometer with an FOV ofin the SDSS g′,r′,i′bands,respectively.The efficiency of each band is 65.6%,63.7%and 68.3%,respectively.The limiting magnitudes of the g′,r′,i′ bands are 15.46,16.89 and 15.93,and the accuracy of differential photometry of each band in terms of rms is 1.59%,1.48%and 1.62%,respectively.In the future,we will upgrade the system software to do the data processing after receiving the images generated by the CCD camera so as to improve the work efficiency.Moreover,we will try to optimize the image quality of the edge FOV to make good use of the photometer for observation tasks,expecting better performance.

    Acknowledgments

    This work is a part of research done by the Astronomical Joint Fund of the National Natural Science Foundation of China (NSFC,Grant No.U2031129),together with the Youth Innovation Promotion Association of the Chinese Academy of Sciences (2018079),and Project of the National Center for Astronomy,Chinese Academy of Sciences (2019).

    Thanks go to Xiaomeng Lu and Chao Wu at National Astronomical Observatories for their help in this project,and also to Bowen Li,Guangsheng Zhang and Chunlei Guo at Jilin Astronomical Observatory for their hard work.

    美女脱内裤让男人舔精品视频| 在线播放国产精品三级| 国产视频内射| 日韩 亚洲 欧美在线| 日本免费a在线| 校园人妻丝袜中文字幕| 狠狠狠狠99中文字幕| 免费观看性生交大片5| 欧美zozozo另类| 99久久无色码亚洲精品果冻| 精品久久久久久久末码| 午夜免费激情av| 精品欧美国产一区二区三| 亚洲最大成人手机在线| 亚洲aⅴ乱码一区二区在线播放| 国语对白做爰xxxⅹ性视频网站| 1024手机看黄色片| 国国产精品蜜臀av免费| 嫩草影院入口| 久久国内精品自在自线图片| 免费黄网站久久成人精品| 91精品一卡2卡3卡4卡| 女人十人毛片免费观看3o分钟| 日韩强制内射视频| 午夜亚洲福利在线播放| 午夜亚洲福利在线播放| 毛片一级片免费看久久久久| 99久久无色码亚洲精品果冻| 一二三四中文在线观看免费高清| 久久久色成人| 国产成人aa在线观看| 国产乱来视频区| 亚洲成人中文字幕在线播放| 男人舔女人下体高潮全视频| 精品午夜福利在线看| 日本猛色少妇xxxxx猛交久久| 99在线人妻在线中文字幕| 日韩强制内射视频| 视频中文字幕在线观看| 国产精品女同一区二区软件| 青春草亚洲视频在线观看| 免费大片18禁| 久久韩国三级中文字幕| 国产成人a∨麻豆精品| 91久久精品国产一区二区成人| 亚洲国产精品成人综合色| 卡戴珊不雅视频在线播放| 国产精品久久久久久av不卡| 亚洲国产色片| 少妇人妻精品综合一区二区| 亚洲一区高清亚洲精品| 三级国产精品片| 中文字幕人妻熟人妻熟丝袜美| 国产 一区 欧美 日韩| 亚洲在线自拍视频| 麻豆成人av视频| 老司机影院毛片| 国产免费福利视频在线观看| 婷婷色麻豆天堂久久 | 91精品伊人久久大香线蕉| 亚洲av免费高清在线观看| 亚洲人与动物交配视频| 欧美成人免费av一区二区三区| 亚洲在线自拍视频| 欧美日韩国产亚洲二区| 夫妻性生交免费视频一级片| 天堂中文最新版在线下载 | 亚洲欧美日韩卡通动漫| 国产精品久久久久久av不卡| 久久国内精品自在自线图片| 能在线免费观看的黄片| 国产精品麻豆人妻色哟哟久久 | 可以在线观看毛片的网站| 国产成人freesex在线| 国语自产精品视频在线第100页| 亚洲成人av在线免费| 日本黄大片高清| 国产亚洲5aaaaa淫片| 亚洲成人中文字幕在线播放| 高清日韩中文字幕在线| 高清日韩中文字幕在线| 久久久精品94久久精品| 天天躁日日操中文字幕| 精品欧美国产一区二区三| 秋霞在线观看毛片| 晚上一个人看的免费电影| 十八禁国产超污无遮挡网站| 亚洲人成网站在线观看播放| 一个人看的www免费观看视频| 天天躁夜夜躁狠狠久久av| 色综合亚洲欧美另类图片| 青春草亚洲视频在线观看| 欧美激情在线99| 91久久精品国产一区二区成人| 男插女下体视频免费在线播放| 国产v大片淫在线免费观看| 亚洲在线观看片| 国产午夜精品论理片| 99视频精品全部免费 在线| 亚洲成人精品中文字幕电影| 亚洲欧洲国产日韩| 精品久久国产蜜桃| 国产免费视频播放在线视频 | 午夜福利视频1000在线观看| 亚洲精品国产av成人精品| 一级黄片播放器| 亚洲人成网站在线播| 国产免费又黄又爽又色| 91狼人影院| av天堂中文字幕网| 人体艺术视频欧美日本| 久久久亚洲精品成人影院| 国内揄拍国产精品人妻在线| 久久99蜜桃精品久久| 国产精品国产三级国产av玫瑰| 少妇的逼好多水| 两个人视频免费观看高清| 一级毛片我不卡| 99热全是精品| 久久久久久久久久久丰满| 国产激情偷乱视频一区二区| 免费看光身美女| 亚洲欧美精品自产自拍| 国产三级在线视频| 亚洲精品乱久久久久久| 一夜夜www| 日本一二三区视频观看| 国产精品一二三区在线看| 国产又黄又爽又无遮挡在线| 内射极品少妇av片p| ponron亚洲| 亚洲av成人精品一二三区| 国国产精品蜜臀av免费| 国产高清有码在线观看视频| 天天躁日日操中文字幕| 日韩一区二区三区影片| 久久久久久久亚洲中文字幕| 日韩一区二区视频免费看| 小说图片视频综合网站| 中文欧美无线码| 丰满乱子伦码专区| 亚洲综合精品二区| 美女高潮的动态| av黄色大香蕉| 91精品一卡2卡3卡4卡| 身体一侧抽搐| eeuss影院久久| 国产精品日韩av在线免费观看| 久久综合国产亚洲精品| 午夜免费激情av| 中文字幕免费在线视频6| 人人妻人人看人人澡| 国产av码专区亚洲av| 国产三级在线视频| 国产私拍福利视频在线观看| 国产高清国产精品国产三级 | 午夜福利在线观看吧| 青春草国产在线视频| 国产精品永久免费网站| 亚洲成av人片在线播放无| 成人高潮视频无遮挡免费网站| 亚洲中文字幕一区二区三区有码在线看| 亚洲真实伦在线观看| 亚洲熟妇中文字幕五十中出| 久久99热这里只频精品6学生 | 免费av观看视频| 国产极品精品免费视频能看的| 久久99热这里只有精品18| 一级毛片电影观看 | 国产综合懂色| 白带黄色成豆腐渣| 国产女主播在线喷水免费视频网站 | 亚洲av熟女| 国产精品一区www在线观看| 高清在线视频一区二区三区 | 久久精品久久久久久久性| 嘟嘟电影网在线观看| 欧美97在线视频| 亚洲无线观看免费| 国产白丝娇喘喷水9色精品| 91aial.com中文字幕在线观看| 少妇的逼好多水| 色播亚洲综合网| 国产视频内射| 联通29元200g的流量卡| 一级黄色大片毛片| 国产午夜精品久久久久久一区二区三区| 长腿黑丝高跟| 免费看美女性在线毛片视频| 最近2019中文字幕mv第一页| 国产精品麻豆人妻色哟哟久久 | 久久久久久久久久久免费av| 精华霜和精华液先用哪个| 国产精品三级大全| 日韩高清综合在线| 干丝袜人妻中文字幕| 日日摸夜夜添夜夜添av毛片| 麻豆久久精品国产亚洲av| 高清日韩中文字幕在线| 成人午夜高清在线视频| 免费在线观看成人毛片| 国产探花在线观看一区二区| 少妇高潮的动态图| 人人妻人人澡欧美一区二区| 99在线人妻在线中文字幕| 一边摸一边抽搐一进一小说| 女的被弄到高潮叫床怎么办| 日韩成人伦理影院| 身体一侧抽搐| 我的女老师完整版在线观看| 狠狠狠狠99中文字幕| 成人鲁丝片一二三区免费| 亚洲成av人片在线播放无| 神马国产精品三级电影在线观看| 99在线人妻在线中文字幕| 国产精品熟女久久久久浪| 国产一区有黄有色的免费视频 | 真实男女啪啪啪动态图| 特级一级黄色大片| 男人舔奶头视频| 亚洲欧洲日产国产| 禁无遮挡网站| 精品99又大又爽又粗少妇毛片| 99久久九九国产精品国产免费| 国产三级中文精品| 国国产精品蜜臀av免费| 麻豆乱淫一区二区| 亚洲人成网站在线观看播放| 人人妻人人澡欧美一区二区| 舔av片在线| 又爽又黄无遮挡网站| 亚洲最大成人中文| 五月伊人婷婷丁香| 日韩av在线大香蕉| 日本三级黄在线观看| 色综合色国产| 国产不卡一卡二| 亚洲av熟女| 免费黄网站久久成人精品| 久久久久久久久久久免费av| av免费在线看不卡| 国产亚洲av嫩草精品影院| 真实男女啪啪啪动态图| 亚洲第一区二区三区不卡| 欧美精品国产亚洲| 国产精品美女特级片免费视频播放器| 国产精品永久免费网站| 日韩,欧美,国产一区二区三区 | 国产视频内射| 国产伦在线观看视频一区| 欧美日本亚洲视频在线播放| 日本av手机在线免费观看| 女人被狂操c到高潮| 黄色配什么色好看| 国产伦精品一区二区三区视频9| 男女视频在线观看网站免费| 国产久久久一区二区三区| 精品国产一区二区三区久久久樱花 | 欧美高清成人免费视频www| 久久精品人妻少妇| 久久久精品94久久精品| 韩国高清视频一区二区三区| 两性午夜刺激爽爽歪歪视频在线观看| 久久精品夜夜夜夜夜久久蜜豆| 97热精品久久久久久| 天堂网av新在线| 1000部很黄的大片| 国国产精品蜜臀av免费| 99热6这里只有精品| 亚洲国产日韩欧美精品在线观看| 人人妻人人澡欧美一区二区| 一区二区三区四区激情视频| 久久久久久久久久久丰满| 夜夜看夜夜爽夜夜摸| 国产成年人精品一区二区| 午夜爱爱视频在线播放| 精品久久久久久久人妻蜜臀av| 91午夜精品亚洲一区二区三区| 国产高清国产精品国产三级 | 在线播放无遮挡| 中文欧美无线码| 男女视频在线观看网站免费| 直男gayav资源| 色网站视频免费| 国产精品一区二区三区四区久久| 久久久国产成人免费| 国产午夜精品一二区理论片| 国产三级在线视频| 狠狠狠狠99中文字幕| 级片在线观看| 你懂的网址亚洲精品在线观看 | 亚洲国产精品久久男人天堂| 国产精品久久久久久久电影| 亚洲熟妇中文字幕五十中出| 久久久久网色| av在线亚洲专区| 五月玫瑰六月丁香| 禁无遮挡网站| 欧美变态另类bdsm刘玥| 日韩欧美精品v在线| 日韩成人av中文字幕在线观看| 亚洲av不卡在线观看| 欧美不卡视频在线免费观看| 最新中文字幕久久久久| 春色校园在线视频观看| a级一级毛片免费在线观看| 亚洲自拍偷在线| 丝袜美腿在线中文| 又黄又爽又刺激的免费视频.| av在线蜜桃| 国产精品三级大全| 又爽又黄无遮挡网站| 国产毛片a区久久久久| 国产精品国产三级国产av玫瑰| 免费在线观看成人毛片| 亚洲最大成人中文| 亚洲一级一片aⅴ在线观看| 欧美色视频一区免费| 国产伦在线观看视频一区| 久久精品夜夜夜夜夜久久蜜豆| 国语对白做爰xxxⅹ性视频网站| 天堂网av新在线| 亚洲精品,欧美精品| 毛片一级片免费看久久久久| 中文字幕亚洲精品专区| 免费人成在线观看视频色| 日本免费一区二区三区高清不卡| 在线免费观看的www视频| 六月丁香七月| 国产伦精品一区二区三区四那| 亚洲av中文av极速乱| 天美传媒精品一区二区| 高清在线视频一区二区三区 | 高清日韩中文字幕在线| 日本欧美国产在线视频| 欧美潮喷喷水| 村上凉子中文字幕在线| 欧美xxxx黑人xx丫x性爽| 中文字幕久久专区| 日本色播在线视频| 日韩亚洲欧美综合| 美女被艹到高潮喷水动态| 又粗又爽又猛毛片免费看| 搡老妇女老女人老熟妇| 成年av动漫网址| 大话2 男鬼变身卡| 久久99热这里只频精品6学生 | 永久网站在线| a级毛片免费高清观看在线播放| 韩国高清视频一区二区三区| 国产熟女欧美一区二区| 日韩强制内射视频| 日本三级黄在线观看| 天天躁夜夜躁狠狠久久av| av视频在线观看入口| 中文字幕久久专区| 成人三级黄色视频| 国产精品嫩草影院av在线观看| 国产伦一二天堂av在线观看| 日韩欧美三级三区| 色视频www国产| 美女脱内裤让男人舔精品视频| 九九热线精品视视频播放| videos熟女内射| av在线观看视频网站免费| 国产在视频线精品| 中文资源天堂在线| 丰满人妻一区二区三区视频av| 欧美3d第一页| 日本-黄色视频高清免费观看| av在线天堂中文字幕| 十八禁国产超污无遮挡网站| 日本午夜av视频| 久久久国产成人免费| 成年版毛片免费区| 国产高清国产精品国产三级 | 少妇熟女aⅴ在线视频| 久久国产乱子免费精品| 亚洲国产最新在线播放| 亚洲精华国产精华液的使用体验| 国内揄拍国产精品人妻在线| 成人av在线播放网站| 色网站视频免费| 69人妻影院| 国产日韩欧美在线精品| 亚洲精品成人久久久久久| 简卡轻食公司| 亚洲经典国产精华液单| 久久久精品94久久精品| 国产乱人偷精品视频| 欧美又色又爽又黄视频| 最后的刺客免费高清国语| 国产欧美日韩精品一区二区| 日韩精品有码人妻一区| 国产精品美女特级片免费视频播放器| 久久久久久久久久久免费av| 久久精品综合一区二区三区| 蜜桃久久精品国产亚洲av| 波多野结衣高清无吗| 综合色丁香网| 国产成人精品一,二区| 免费看a级黄色片| 国产亚洲av嫩草精品影院| 久久精品国产亚洲av涩爱| 亚洲国产欧美在线一区| 热99re8久久精品国产| 夜夜看夜夜爽夜夜摸| 免费搜索国产男女视频| 麻豆一二三区av精品| 在线观看一区二区三区| 成人性生交大片免费视频hd| 少妇熟女aⅴ在线视频| 亚洲av成人av| 欧美成人午夜免费资源| 日本黄大片高清| 丰满少妇做爰视频| 国产成人freesex在线| 九色成人免费人妻av| 乱码一卡2卡4卡精品| 亚洲精品,欧美精品| 日本欧美国产在线视频| 不卡视频在线观看欧美| 变态另类丝袜制服| 中文在线观看免费www的网站| 亚洲精品亚洲一区二区| 国内精品宾馆在线| 成人高潮视频无遮挡免费网站| 久久久久久久久久黄片| 免费看av在线观看网站| 在线观看66精品国产| 深爱激情五月婷婷| 欧美潮喷喷水| 婷婷色麻豆天堂久久 | 国产精品野战在线观看| 亚洲av.av天堂| 一级av片app| 纵有疾风起免费观看全集完整版 | 精品久久久久久电影网 | 亚洲精品乱码久久久久久按摩| 欧美zozozo另类| 人人妻人人看人人澡| 大香蕉97超碰在线| 麻豆成人av视频| 国产av不卡久久| 亚洲av不卡在线观看| 婷婷色麻豆天堂久久 | 国产熟女欧美一区二区| 91久久精品国产一区二区三区| 亚洲美女搞黄在线观看| 自拍偷自拍亚洲精品老妇| 色吧在线观看| 人人妻人人澡欧美一区二区| 一级毛片我不卡| 免费av观看视频| 国产成人精品婷婷| 尤物成人国产欧美一区二区三区| 成人三级黄色视频| 熟女人妻精品中文字幕| 九九在线视频观看精品| 天美传媒精品一区二区| 在线播放无遮挡| 中文精品一卡2卡3卡4更新| 精品99又大又爽又粗少妇毛片| av在线观看视频网站免费| 精品人妻熟女av久视频| 亚洲国产高清在线一区二区三| 18禁在线无遮挡免费观看视频| 91狼人影院| 国产 一区精品| 国产乱来视频区| 亚洲中文字幕一区二区三区有码在线看| 精品久久久久久久末码| 免费看光身美女| 淫秽高清视频在线观看| 老师上课跳d突然被开到最大视频| 村上凉子中文字幕在线| 日韩亚洲欧美综合| 国产老妇女一区| 69人妻影院| 日韩在线高清观看一区二区三区| 爱豆传媒免费全集在线观看| 成人亚洲精品av一区二区| 国产三级中文精品| 狠狠狠狠99中文字幕| 久久亚洲国产成人精品v| 亚洲18禁久久av| 又爽又黄a免费视频| 精品久久久久久久久亚洲| 免费在线观看成人毛片| 两个人的视频大全免费| 男人和女人高潮做爰伦理| 欧美最新免费一区二区三区| 91在线精品国自产拍蜜月| 亚洲无线观看免费| 听说在线观看完整版免费高清| 日韩一区二区视频免费看| 国产精品国产三级国产av玫瑰| 亚洲图色成人| 久久99蜜桃精品久久| 国产老妇伦熟女老妇高清| 欧美最新免费一区二区三区| 久久99热这里只频精品6学生 | 日韩精品青青久久久久久| 国产淫片久久久久久久久| 中国美白少妇内射xxxbb| 成人特级av手机在线观看| 可以在线观看毛片的网站| 精品午夜福利在线看| 久久人妻av系列| 国产一区二区三区av在线| 免费观看的影片在线观看| 美女脱内裤让男人舔精品视频| or卡值多少钱| 免费人成在线观看视频色| 丰满人妻一区二区三区视频av| 欧美成人免费av一区二区三区| 国产三级中文精品| 国产精品美女特级片免费视频播放器| av在线蜜桃| av女优亚洲男人天堂| 亚洲第一区二区三区不卡| 国产精品熟女久久久久浪| 日本-黄色视频高清免费观看| 国语对白做爰xxxⅹ性视频网站| 国产白丝娇喘喷水9色精品| 欧美性感艳星| 欧美+日韩+精品| 日韩强制内射视频| 久久久久久久久中文| 婷婷色综合大香蕉| 99热全是精品| 国产午夜精品久久久久久一区二区三区| 国产精品美女特级片免费视频播放器| 国产女主播在线喷水免费视频网站 | 亚洲欧美日韩无卡精品| a级毛色黄片| 亚洲精品日韩在线中文字幕| 久久99热6这里只有精品| 三级国产精品欧美在线观看| 男人的好看免费观看在线视频| av.在线天堂| 毛片女人毛片| 1000部很黄的大片| 国产人妻一区二区三区在| 国产三级中文精品| 99热这里只有是精品在线观看| 男人舔奶头视频| 成人鲁丝片一二三区免费| 看非洲黑人一级黄片| 欧美高清性xxxxhd video| 午夜a级毛片| 亚洲国产精品国产精品| 亚洲av男天堂| 成人综合一区亚洲| 久久精品综合一区二区三区| 日韩成人伦理影院| 天美传媒精品一区二区| 国产毛片a区久久久久| 一个人观看的视频www高清免费观看| 亚洲经典国产精华液单| av国产久精品久网站免费入址| 人体艺术视频欧美日本| 国产精华一区二区三区| 永久免费av网站大全| 少妇人妻精品综合一区二区| 久久99精品国语久久久| 麻豆精品久久久久久蜜桃| 久久99精品国语久久久| 如何舔出高潮| 菩萨蛮人人尽说江南好唐韦庄 | av.在线天堂| 国产老妇伦熟女老妇高清| 午夜福利在线在线| 成人二区视频| 日韩中字成人| 午夜福利成人在线免费观看| 国产成人91sexporn| 日本一本二区三区精品| 变态另类丝袜制服| 日韩成人伦理影院| 成人午夜高清在线视频| 黄色日韩在线| av国产久精品久网站免费入址| 天堂中文最新版在线下载 | 一个人看的www免费观看视频| 国语对白做爰xxxⅹ性视频网站| 国产免费福利视频在线观看| 午夜激情欧美在线| 欧美性猛交╳xxx乱大交人| 18+在线观看网站| 午夜亚洲福利在线播放| 美女脱内裤让男人舔精品视频| 午夜精品在线福利| 91午夜精品亚洲一区二区三区| 亚洲一区高清亚洲精品| 国产精品av视频在线免费观看| 国产乱人偷精品视频| 国产精品一区二区三区四区久久| 亚洲丝袜综合中文字幕| 我要搜黄色片| 亚洲色图av天堂| 禁无遮挡网站| 久久欧美精品欧美久久欧美| 精品99又大又爽又粗少妇毛片| 我的女老师完整版在线观看| a级一级毛片免费在线观看| 亚洲av电影在线观看一区二区三区 | 婷婷色综合大香蕉| 亚洲不卡免费看| 亚洲国产欧洲综合997久久,| 少妇的逼水好多| 国产黄色小视频在线观看| 夫妻性生交免费视频一级片| 久久久国产成人免费| 汤姆久久久久久久影院中文字幕 |