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

    L-band RFI Measurement and Mitigation at the NSRT Site?

    2019-12-10 11:58:56LIUQiLIUYeCHENMaozhengWANGYueLIUFengYANHaoCAOLiangSUXiaoming
    天文學(xué)報(bào) 2019年6期

    LIU Qi LIU Ye CHEN Mao-zheng WANG Yue LIU Feng YAN Hao CAO Liang SU Xiao-ming

    (1 Xinjiang Astronomical Observatory, Chinese Academy of Sciences, Urumqi 830011)

    (2 University of Chinese Academy of Sciences, Beijing 100049)

    ABSTRACT Radio environment degrades significantly at Nan Shan 26 m Radio Telescope(NSRT)base due to new and upgraded electrical devices on-site and more establishment of surrounding radio communication services.How do we mitigate them effectively? Firstly, a quasi-real time Radio Frequency Interference (RFI) measurement method is employed for RFI detection and spectral analysis based on the layout for potential on-site interference areas and off-site radio communication services.Our investigation at the frequency bands for pulsar observations (1380–1700 MHz) indicates that transient interferences pose the greatest impact on observations.Then we use a portable RFI measurement system for interferences hunting at the NSRT, and the typical characteristics of interferences are analyzed.Secondly, based on the RFI characteristics, main interference sources in the observation room are shielded with shielding cabinets.The results from our shielding efficiency(SE)measurement and evaluation show that these methods work effectively.Furthermore, we propose a Radio Quiet Zone (RQZ)and the regulations for NSRT to mitigate off-site RFIs.

    Key words radio environment, RFI characteristic, RFI mitigation, shielding efficiency

    1 Introduction

    The NSRT base, as shown in the Fig.1, is located at the northern foot of Tianshan Mountain in the center of Eurasian Continent with the longitude and latitude of87?10.67′E and 43?28.27′N, respectively.It is situated at an altitude of 2080 m, and more than 70 km away from Urumqi city.The NSRT is equipped with four cryogenic receivers operating at L, S/X, C, and K bands (see Table 1, in which ST is the System Temperature of the receiver, LCP and RCP are the Left-hand Circular Polarization and Right-hand Circular Polarization, respectively).The unique geographical location makes the telescope to be an important node of Very Long Baseline Interferometry(VLBI)observations for both domestic and international observations.In addition, the NSRT contributes the fundamental research in radio astronomy as well,such as pulsars,molecular lines, and active galactic nuclei.

    Fig.1 The NSRT site

    Table 1 Specifications of the receivers on NSRT

    Radio environment of NSRT is deteriorated for several reasons as follows.Firstly,the recent upgrades on the NSRT include the installment of various new electronic equipments, such as monitoring and controlling electronics, digital backends, safety equipments, and optical telescopes, which introduce potential on-site RFI sources.The Fig.2(a) shows the distribution of the on-site potential RFI sources at the NSRT base.However, shield protection and effective equipment management are ignored during this process[1].Secondly, the range of frequency occupancy is getting broader with the development of technology.The Fig.2(b) shows the distribution of off-site potential RFI sources, which indicates that many transmitters are operating around the NSRT base.Therefore, it is a challenge to identify the RFI sources and their associated signal features for the RFI mitigation.

    Fig.2 Overview of the potential on-site (a) RFI sources, about 600 m range around the NSRT, and off-site (b) RFI sources, where a circular region with a radius of 5 km is indicated by the yellow boundary, and the spots in cyan are mobile phone base stations.

    2 RFI detection and characteristic analysis

    2.1 RFI measurement system

    A portable RFI measurement system is employed for this work.Fig.3 shows the diagram of this system,including a Schwarz-beck 9143 directional antenna,a broadband Low Noise Amplifier (LNA) with maximum gain of 30 dB, a Personal Computer (PC)controlling the R & S FSW13 spectrum analyzer with the General Purpose Interface Bus (GPIB) cable, and a standard noise source with the 24 V Direct Current (VDC)power supply interface.We develop a program on the Visual Studio 2013 platform for the control of spectrum analyzer and data storage.The performance of this system is tested by the Y-factor method[2]showing in Fig.4,indicating that the system sensitivity is good enough for strong RFIs detection.

    2.1.1 Measurement methods

    (1) The RFI detection system is installed on the NSRT at a height close to the feed aperture for both horizontal and vertical polarizations RFI measurements.The RFI levels obtained at this position can be used to analyze the RFI impacting on observations effectively.

    Fig.3 Diagram of the RFI measurement system

    Fig.4 System Gain (SG) and system noise temperature

    (2) Key parameters set for spectrum analyzer are as follows: Resolution Band Width (RBW)=30 kHz, Video Band Width (VBW)=300 kHz, we obtain ten points at every RBW, and the sweep time at one RBW is 100μs; Sample detector and the average trace mode are applied for the measurement,the average time is set to 200 times to improve the signal-to-noise ratio.From the above parameters, the sweep time for one span can be calculated for testing.

    (3) In order to improve the RFI data accuracy, system gain and system noise temperature are measured beforehand.We use the system gain for data calibration;and we can also determine the system reliability according to the system noise temperature.

    (4) The 3 dB beam width of the measuring antenna is about 45?, therefore, RFIs are tested by controlling the NSRT with 45?intervals at a fixed elevation of 10?,which ensures no obstructions impacting on RFI measurement.It is assumed that the northward direction is 0?and the southward direction is 180?.

    (5) A quasi-real-time method is used for daytime RFI measurements from UTC(Coordinated Universal Time) +8 (Beijing time) 10:00 to 22:00 with 1.5 h intervals for one polarization[3], which enables us to analyze the RFI variations against time.

    2.1.2 Data processing and analyzing

    Data processing: we calculate the average values of every 10 points at each RBW band as the final levels corresponding to the channel frequency.

    Data calibration: assuming RFI level at antenna aperture isP, in unit of dBm,which can be calculated according to Eq.(1),

    wherePSAis the signal level obtained by the spectrum analyzer, in unit of dBm;GAis the gain of the testing antenna;GSis the gain of the RFI measurement system except the testing antenna.

    ThenPis further converted into field strengthE,in unit of dBμV·m?1,based on Eq.(2),

    where AF is the antenna factor.

    Interference signals of the spectra from all directions at one time period are counted.The statistical graph is presented in the Fig.5, which shows that narrow signals from 180?,225?,and 270?are greater than those from other directions.Based on the analysis of the distribution of the potential on-site RFI sources as mentioned in the introduction section, most of the RFIs are suspected to come from on-site electronics.

    Fig.5 Statistic of the interference signals from all directions.The black color is the background noise,and the white lines are the signals identified.

    The Frequency Channel Occupancies (FCOs) for pulsar observation band (1400–1700 MHz)are counted as well.The bandwidth of each channel in the digital filter-band is 512 kHz, while the RFI frequency resolution is 30 kHz.So during the statistical process, if one 512 kHz bandwidth in radio environment spectra existing in interference signals, it can be seen that this channel is occupied.The FCO is defined as the percentage of channels occupied vs.the total channels in this paper.The FCO of spectra from all directions and time periods are counted, as presented in the Fig.6[3].This statistical graph shows the FCO varies greatly in the same direction, indicating that transient RFI has great impact on observations.

    2.1.3 Interference source analysis

    In order to verify the RFI sources, measurement steps are as follows: (1) the RFI testing system is installed on the NSRT close to the feed aperture, pointing at the equipment areas (180?).And the measuring antenna is about 30 m away from the equipment areas; then the spectra from both the on-site equipment area and off-site radio environment are tested.(2) The RFI testing system is installed on the top of the most external building at the same height crossing to the equipment areas, pointing to the same direction, as shown in the Fig.7.Then off-site radio environment can be measured in the same direction.A comparison of the spectra between the equipment area and off-site radio environment is presented in Fig.8, indicating that most of the narrow RFIs come from the on-site equipment areas.

    Fig.6 This plot presents frequency channel occupancy vs.time and direction, the frequency channel occupancy of all spectra is marked on the plot for clearance.

    Fig.7 RFI source verification

    Furthermore, the primary RFI sources in equipment areas are tested using a portable on-site RFI measurement system[4].The distance for measurement is set to 1 m; the sample detector and average trace mode are employed for radiate emission testing.Measurement results indicate that the primary RFI sources include servos,network equipment, backends, digital switches, industrial computers and so on.Typical RFI sources and their characteristics are presented in Fig.9,which shows that their emission characteristics are always broadband below 3 GHz; meanwhile emissions from some digital devices are extremely noisy compared to the environment noise(Fig.9(b)).

    Fig.8 Two-state measurement spectra between equipment areas and RFI environment in the direction of 180?.(a) and (b) are the vertical and horizontal polarization measured results, respectively.The black line is the spectra tested on the telescope, and the dotted and gray lines are the RFI environment spectra measured on the top of the building.

    Fig.9 Radiate emissions for the typical digital equipment, (a) shows the emission characteristic of the antenna servo, (b) presents the emission characteristic of the typical backend, where the gray line is the environment noise tested by shutting down all electrical devices.

    2.2 RFI measurements using the NSRT

    Due to NSRT is more sensitive than the RFI measurement system, in order to further hunt and analyze the characteristics of RFI signal, the NSRT is employed for RFI detection at L band.The Fig.10 shows the measurement system links, and the intermediate frequency (IF) spectra from the L-band receiver obtained by the R & S FSW13 spectrum analyzer.For practical testing, the spectra are obtained at 15?intervals in the case of elevation=10?.The measurements analyzing shows the main interference comes from the off-site radio communication services, such as broadcasting satellite,Global Positioning System (GPS), Global Navigation Satellite System (GLONASS),Inmarsat, cell phone stations, military services, and on-site electronics according to the International Telecommunication Union (ITU) radio regulations1http://search.itu.int/history/HistoryDigitalCollectionDocLibrary/1.43.48.en.101.pdfand our on-site interference measurements.

    From the detected RFIs and the analyzed results shown above, we select several strong RFI sources for further analysis of their characteristics, especially for the 1440–1495 MHz band, which have great impacts on the pulsar observation.

    (1) The extremely strong RFI sources come from the azimuth at 255?, as shown in the Fig.11, this broadband interference can be identified in the azimuthal direction from 225?to 270?with the NSRT, and the frequency bands are 1460–1480 MHz,1485–1495 MHz, and 1630–1650 MHz.The interference also makes the L-band receiver saturation when NSRT operating at the azimuthal direction between 255?and 245?,and at the elevation range from 8?to 13?.For this interference, we have sought for its origin several times in the mountains along the direction of 255?.However, we can hardly confirm its source because of the complex terrain.But the azimuthal range that leads to receiver saturation is coded in the pulsar observation software to avoid the strong RFI, the efficiency of pulsar observation is thus increased.

    Fig.10 RFI measurement diagram using the NSRT, the IF signals from L-band receiver were sent to spectrum analyzer.

    Fig.11 Strong RFI in the direction of 255?

    (2) Another strong RFI sources come from the 155?, which can be identified with NSRT in the azimuthal direction ranging from 130?to 190?, covering from 1467 to 1492 MHz,including five sub-bands(1467–1471 MHz,1471–1476 MHz,1476–1479 MHz,1479–1484 MHz, 1484–1492 MHz), as shown in the Fig.12.According to the broadcasting-satellite service[5], this interference comes from satellite digital audio radio uplink service.However, the location of the transmitter could not be determined due to little information provided by local government and the complex terrain factor.Therefore, we should further communicate with the local government for coordinating on this RFI in the future.

    (3)When using the NSRT for RFI measurements, we notice that a strong interference appears at the L band between 1400 and 1500 MHz,which does not vary with AZ.Obviously, this interference impacts on observations greatly.From our experience, the interference may come from the telescope itself.In order to search for this RFI source,the portable interference testing system is used to measure the possible equipment area of the telescope.And we find that the RFI source comes from the newly installed IF switching system in the receiver cabinet.When this system is shut down, the interference signal in the L-band spectrum disappears, as shown in the Fig.13.Now, the IF switching system is turned off during pulsar observations,and we are planning to shield this system in the future.

    Fig.12 RFI characteristics of the broadcasting satellite

    Fig.13 Spectra of the IF Switch during on and off

    (4) Other strong on-site RFI characteristics are presented in the Fig.14(a), which shows that the emissions can enter the telescope receive system directly.The Fig.14(a)shows a lot of narrow signals in the pulsar observing band with an interval of 1.25 MHz from 1580 to 1700 MHz in the azimuthal direction of 195?.Based on the RFI characteristics and the measurement results, the interference source is located coming from the newly installed Fast Radio Burst (FRB) backend close to the window in NRST observing room.The Fig.14(a) indicates that the interference disappeared when the system is shut down.Furthermore we shield this interference source with a shielding cabinet.

    In addition, we find that this interference can be detected obviously when NSRT operating at the azimuth between 120?and 270?, and with the Elevation of 10?, as shown in Fig.14(b).We speculate that this RFI may enter the telescope through the antenna side lobe and the sub-reflector reflections.In order to figure out how the interferences are coupled into the telescope, we need to investigate this issue by further measurement and simulation analysis.

    Fig.14 RFI characteristics of the FRB backend, Fig.14(a) presents the spectra during FRB backend on and off, Fig.14(b) shows the spectra vs.different azimuth.

    3 RFI mitigation

    3.1 On-site RFI mitigation

    Based on the detected RFIs and their characteristics at NSRT site,it is not easy to shield all interference sources in a short time.So we are trying to shield the main RFI sources in the NSRT observing room within our funds available.Shielding efficiency requirements for the RFI source analysis are as follows: (1) the RFI mitigation is for pulsar observations, so the shielded frequency band we cared about is 1.4–1.7 GHz.(2) Interference level limit is?175 dBm for the continuum observations at 1.4 GHz[6];(3) The maximum radiate emissions of the digital equipment counted is?75 dBm at the L band; (4) the NSRT observing room is about 37 m away from the telescope feed aperture, and the path loss calculated with free-space propagation model is about 64 dB; (5) the antenna gain calculated according to the parabolic antenna gain model is about?5 dBi[7?9].The shielding efficiency (SE)SErequirements can be calculated according to Eq.(3),

    whereSEis the SE requirements, TL is the interference level limit at feed aperture of telescope,Gsidelobeis the antenna gain,GLossis the path loss from the observing room to telescope,PEis the measured radiate emission of the equipment,Uis measurement uncertainty.By assuming the measurement uncertainty is 6 dB, the largest calculated SE requirement is 37 dB.

    We propose a shield scheme for the primary RFI sources (including Intra-Day Variability(IDV),FRB,Digital Filter Bank(DFB),Analog Filter Bank(AFB)backend;time and frequency system;the controlling computers for Antenna Control Unit(ACU),Field System (FS), VLBI system, and others) in the observing room using shielding boxes, as shown in Fig.15.And we customize the power line and signal cable filters that installing on interface boards for every cabinet to avoid conducted emission over cables.In addition, in order to realize the network supplying and shield the ethernet devices effectively,the optical fiber links are employed to replace the previously ethernet cables.

    Fig.15 Shielding schematics using RFI cabinets, the dashed line box represents shielding cabinets, the fibers from the switch to the shielding box for ethernet supply.

    In order to verify the effectiveness of electromagnetic shielding enclosures, the SE of the RFI cabinets is measured on-site according to the IEEE standard method[10],and the SE is calculated on the basis of Eq.(4),

    whereP1is the reference power intensity without the enclosure;P2is power intensity with the enclosure in place, in unit of dBm.

    SE measurement results for the doors and interface boards of the shielding cabinets are presented in Fig.16.The figure indicates that the SE of the shielding box are more than 60 dB below 3 GHz, so testing results meet the shield requirements.

    In order to further verify the effectiveness of this shielding work, radiate emissions of the typical interference are tested between the shielding doors opening and closing statuses when the interference sources are installed in the shielding cabinets.And the differences of spectral statuses are presented in the Fig.17, which shows that the emissions can be mitigated more than 40 dB between 1000–2000 MHz.Due to the radiate emission levels of RFI sources and the sensitivity of RFI measurement system,this verification method is not the standard SE measurement method as well, which leads to the SE measurement differences between Fig.16 and Fig.17.However, all measurement results indicate that the SE is enough based on the requirement evaluation analysis above.

    Fig.16 SE measured for shielding cabinet at the cable connected status

    Fig.17 Spectra measured when the shield door is open vs.close.Fig.17(a) shows the shielding cabinet of the time and frequency system testing results; Fig.17(b) presents the shielding cabinet of pulsar observational backend.

    The antenna servo and control system are strong RFI sources,as shown in Fig.9(a).However,it is a challenge for us to shield them due to its high power and various cables coming into and out of the servo modules.If we shield them using traditional method(shield and filter), the cost will be extremely high.From the analysis above, we are planning to move the servo system from the NSRT observing room to the basement under the telescope such that its radiate emissions can be mitigated by the reinforced concrete of the basement.After this work, we make an evaluation for the impact of the servo emission on observations, the results are as follows: (1) the interference level limit for pulsar observation of the NSRT is?175 dBm; (2) the propagation loss from servo to the L-band receiver feed aperture is 58 dB; (3) the antenna gain calculated according to parabolic antenna gain model is about?8 dBi[8]; (4) The maximum RFI levels at the L band is?95 dBm; (5) Reinforced concrete SE measurement result is more than 20 dB,and the SE measurement for servo cabinet door is about 10 dB.From the calculated and measured results above, the servo emissions propagating to the feed aperture of the NSRT is?191 dBm, which is obviously below the pulsar observation limit, indicating that this RFI protection work is effective.

    3.2 Off-site RFI protection

    In recent years, increasing number of new radio communication facilities were constructed in the off-site area around the NSRT,especially the mobile base stations,which have great impacts on pulsar observations.Therefore, a circular region of Radio Quiet Zone (RQZ)with a radius of 10 km is proposed as shown in the Fig.18(a).The off-site RFI protection is conducted with the following requirements: (1)No new transmitter is allowed in the circular region within a radius of 5 km from the NSRT;(2)All operating transmitters within the view of the telescope, and located within a radius of 5 km need to be removed;(3)No new transmitter within the view of the telescope is allowed in the range of 5 to 10 km region; (4) Local government should provide accurate information for all transmitters available in the RQZ for the interference regulation in the future.

    Fig.18 The RQZ for the NSRT is shown in map (a), and mobile stations are marked with signs in blue and purple in map (b).The purple ones will be removed.

    We are proposing the above requirements to the local governments, and to ensure that we can negotiate with the authorities and the regulatory department responsible for the radio transmitters regulation.Meanwhile, the mobile stations within the view of the telescope, and located within a radius of 5 km, as shown in Fig.18(b), should be removed first.

    4 Conclusions

    A portable RFI measurement system and the NSRT are employed for our RFI environment testing and interference hunting;meanwhile we analyze the characteristics of RFIs.The conclusions are as follows: (1) Transient RFIs have great impacts on observations; (2) Both on-site and off-site RFIs can enter the receiver from different azimuths through side lobe and sub-reflector reflections; (3) Dispite the AZ and EL changes,RFIs can always be coupled into the receiver system of telescope,which greatly impact on the observations.

    Additionally, we are trying to shield the primary RFI sources in the NSRT observing room.The results from the SE measurement and emission evaluation indicate that these mitigation works are effective.

    Although we do some work for the RFI measurements and mitigation, many problems remain to be solved in the future.So our work plans are as follows: (1) It is extremely important to enforce on-site electrical equipment management policy so that the RFIs from new electronics can be controlled effectively; (2) There are still many RFI sources at the NSRT site, so the mitigation work will continue; (3) It is necessary for us to figure out the coupling principle between telescope and interference for the further RFI hunting and mitigating; (4) Off-site RFI coordination with the local government is proceeding; (5) A real-time RFI detection system will be developed for further transient interferences analysis.

    五月天丁香电影| 国产成人aa在线观看| 国产成人精品福利久久| 亚洲精品中文字幕在线视频 | 大片电影免费在线观看免费| 波野结衣二区三区在线| 韩国高清视频一区二区三区| 亚洲色图综合在线观看| 免费av毛片视频| 久久久久久久久大av| 菩萨蛮人人尽说江南好唐韦庄| 日本熟妇午夜| 欧美老熟妇乱子伦牲交| 在线免费十八禁| 超碰av人人做人人爽久久| 97人妻精品一区二区三区麻豆| 欧美精品人与动牲交sv欧美| 天堂中文最新版在线下载 | 又黄又爽又刺激的免费视频.| 久久精品久久久久久久性| 午夜爱爱视频在线播放| 18禁裸乳无遮挡免费网站照片| 中文精品一卡2卡3卡4更新| 岛国毛片在线播放| 国产精品蜜桃在线观看| 69人妻影院| 只有这里有精品99| 国产精品一区二区在线观看99| 国产日韩欧美在线精品| 久久精品夜色国产| 一区二区三区乱码不卡18| 色吧在线观看| 制服丝袜香蕉在线| 亚洲无线观看免费| 国产一区有黄有色的免费视频| 下体分泌物呈黄色| 各种免费的搞黄视频| 国产精品偷伦视频观看了| av.在线天堂| 免费人成在线观看视频色| tube8黄色片| 乱系列少妇在线播放| 69人妻影院| 免费观看性生交大片5| 啦啦啦在线观看免费高清www| 深夜a级毛片| 午夜福利网站1000一区二区三区| 91午夜精品亚洲一区二区三区| 色综合色国产| 成年免费大片在线观看| 中文精品一卡2卡3卡4更新| 国产成人福利小说| 欧美最新免费一区二区三区| 99热这里只有精品一区| 99久久人妻综合| 丰满乱子伦码专区| 婷婷色麻豆天堂久久| 99精国产麻豆久久婷婷| 久久精品国产亚洲av涩爱| 国产色爽女视频免费观看| 国产精品蜜桃在线观看| 熟妇人妻不卡中文字幕| 尤物成人国产欧美一区二区三区| 中文字幕亚洲精品专区| 中文字幕亚洲精品专区| av在线app专区| 亚洲成人一二三区av| 97人妻精品一区二区三区麻豆| 免费看av在线观看网站| 国产亚洲一区二区精品| 久久精品国产鲁丝片午夜精品| 国产欧美日韩精品一区二区| 免费看光身美女| 边亲边吃奶的免费视频| 日本色播在线视频| 精品人妻视频免费看| 日韩不卡一区二区三区视频在线| 爱豆传媒免费全集在线观看| 久久影院123| 在线观看免费高清a一片| 国产精品99久久99久久久不卡 | 国产爱豆传媒在线观看| 99re6热这里在线精品视频| 国产成年人精品一区二区| 色吧在线观看| 亚洲av福利一区| 免费电影在线观看免费观看| 大片电影免费在线观看免费| freevideosex欧美| 内射极品少妇av片p| 久久久久久久久久人人人人人人| 日日摸夜夜添夜夜添av毛片| 免费看av在线观看网站| 亚洲av不卡在线观看| 亚洲国产色片| 2021少妇久久久久久久久久久| 久久韩国三级中文字幕| 一区二区三区乱码不卡18| 在线观看免费高清a一片| 麻豆国产97在线/欧美| 国产精品秋霞免费鲁丝片| 国产午夜精品久久久久久一区二区三区| 最近中文字幕2019免费版| 少妇的逼水好多| 丝瓜视频免费看黄片| 婷婷色av中文字幕| 99热6这里只有精品| 99九九线精品视频在线观看视频| 国产高清国产精品国产三级 | 久久精品久久久久久久性| 春色校园在线视频观看| 在线观看免费高清a一片| 日本午夜av视频| 色网站视频免费| 网址你懂的国产日韩在线| 纵有疾风起免费观看全集完整版| 亚洲一区二区三区欧美精品 | 一区二区三区四区激情视频| 一级二级三级毛片免费看| 日韩欧美精品免费久久| 亚洲怡红院男人天堂| 国产免费一级a男人的天堂| www.av在线官网国产| 久久久欧美国产精品| 日韩在线高清观看一区二区三区| 大香蕉久久网| 国产精品.久久久| 搞女人的毛片| 高清毛片免费看| 日韩欧美精品免费久久| av国产久精品久网站免费入址| 亚洲四区av| 亚洲成人久久爱视频| 久久久久久久大尺度免费视频| 一级av片app| 欧美极品一区二区三区四区| 亚洲av成人精品一区久久| 精品视频人人做人人爽| 久久久欧美国产精品| 日韩国内少妇激情av| av天堂中文字幕网| 新久久久久国产一级毛片| 亚洲精品一区蜜桃| 午夜福利视频精品| 美女主播在线视频| 九草在线视频观看| 看黄色毛片网站| 久热久热在线精品观看| 91精品伊人久久大香线蕉| 一区二区三区精品91| av在线app专区| 成人综合一区亚洲| 免费黄频网站在线观看国产| 亚洲人成网站高清观看| 国产 一区 欧美 日韩| av专区在线播放| 亚洲性久久影院| 成人高潮视频无遮挡免费网站| 黄色欧美视频在线观看| 欧美 日韩 精品 国产| 春色校园在线视频观看| 国产老妇伦熟女老妇高清| 亚洲四区av| 美女cb高潮喷水在线观看| 亚洲精品国产色婷婷电影| 一本色道久久久久久精品综合| 国产精品无大码| 九九在线视频观看精品| 亚洲人成网站在线观看播放| 狠狠精品人妻久久久久久综合| 亚洲欧美一区二区三区国产| 69av精品久久久久久| 嫩草影院新地址| xxx大片免费视频| 日韩伦理黄色片| 最新中文字幕久久久久| 深夜a级毛片| 人妻夜夜爽99麻豆av| 国产欧美亚洲国产| 国产精品.久久久| 成人亚洲精品一区在线观看 | 欧美+日韩+精品| 三级经典国产精品| 99视频精品全部免费 在线| 国产熟女欧美一区二区| 成年女人看的毛片在线观看| 免费看a级黄色片| 自拍偷自拍亚洲精品老妇| 女人久久www免费人成看片| 久久精品国产亚洲av涩爱| 亚洲av中文字字幕乱码综合| 久久鲁丝午夜福利片| 波野结衣二区三区在线| 免费av不卡在线播放| 免费av不卡在线播放| 国产成人freesex在线| 午夜福利在线观看免费完整高清在| 国产黄片美女视频| 制服丝袜香蕉在线| 永久免费av网站大全| 日产精品乱码卡一卡2卡三| 免费播放大片免费观看视频在线观看| 亚洲天堂av无毛| 欧美日韩精品成人综合77777| 又大又黄又爽视频免费| av国产免费在线观看| 精品国产三级普通话版| 啦啦啦在线观看免费高清www| 国产乱人视频| 寂寞人妻少妇视频99o| 午夜精品一区二区三区免费看| 人妻系列 视频| 国产精品99久久99久久久不卡 | 在线播放无遮挡| 欧美成人精品欧美一级黄| 涩涩av久久男人的天堂| 真实男女啪啪啪动态图| 日韩视频在线欧美| 天天一区二区日本电影三级| 久久热精品热| 高清日韩中文字幕在线| 亚洲av电影在线观看一区二区三区 | 免费看日本二区| 国产乱来视频区| 国产精品久久久久久精品古装| 国产91av在线免费观看| 久久久a久久爽久久v久久| 日韩视频在线欧美| 欧美+日韩+精品| 国产精品精品国产色婷婷| 久久精品久久精品一区二区三区| 久久久久网色| 国产一级毛片在线| 高清在线视频一区二区三区| 精品一区二区三卡| 国产毛片在线视频| 国产老妇女一区| 亚洲天堂av无毛| xxx大片免费视频| 18禁在线无遮挡免费观看视频| 日韩 亚洲 欧美在线| 久久精品国产亚洲av天美| av在线亚洲专区| 欧美精品人与动牲交sv欧美| 亚洲精品日韩av片在线观看| 亚洲怡红院男人天堂| 亚洲欧洲日产国产| 午夜老司机福利剧场| 久久精品国产a三级三级三级| 欧美变态另类bdsm刘玥| 寂寞人妻少妇视频99o| 男女下面进入的视频免费午夜| xxx大片免费视频| 男男h啪啪无遮挡| 色视频www国产| 看黄色毛片网站| 国产精品一区二区性色av| 国内揄拍国产精品人妻在线| 天美传媒精品一区二区| 在线免费观看不下载黄p国产| 国产色爽女视频免费观看| 色综合色国产| freevideosex欧美| 国产国拍精品亚洲av在线观看| 能在线免费看毛片的网站| 嘟嘟电影网在线观看| 国产一区有黄有色的免费视频| 亚洲色图av天堂| 在线观看av片永久免费下载| 亚洲av.av天堂| 亚洲综合色惰| 亚洲美女视频黄频| 欧美日韩精品成人综合77777| 国产成人免费观看mmmm| av免费观看日本| 日韩伦理黄色片| 丰满人妻一区二区三区视频av| 我的女老师完整版在线观看| www.色视频.com| 国产成人免费观看mmmm| 九九久久精品国产亚洲av麻豆| 亚洲综合色惰| av在线播放精品| 久久精品久久精品一区二区三区| 国产黄色视频一区二区在线观看| 成人综合一区亚洲| 夫妻性生交免费视频一级片| 亚洲成人精品中文字幕电影| av一本久久久久| 免费黄色在线免费观看| 联通29元200g的流量卡| 久久久久久久久久人人人人人人| 3wmmmm亚洲av在线观看| 91在线精品国自产拍蜜月| 一级黄片播放器| 看黄色毛片网站| 国产有黄有色有爽视频| 交换朋友夫妻互换小说| 久久国产乱子免费精品| 高清日韩中文字幕在线| 99热网站在线观看| 国产黄a三级三级三级人| av免费观看日本| 亚洲最大成人av| 久久韩国三级中文字幕| 在线免费观看不下载黄p国产| 亚洲不卡免费看| 人人妻人人看人人澡| 人体艺术视频欧美日本| 日本与韩国留学比较| 自拍欧美九色日韩亚洲蝌蚪91 | 91精品伊人久久大香线蕉| 午夜激情福利司机影院| 天美传媒精品一区二区| 久久人人爽人人片av| 在线观看人妻少妇| 亚洲人成网站在线观看播放| 岛国毛片在线播放| 寂寞人妻少妇视频99o| 免费播放大片免费观看视频在线观看| 一级二级三级毛片免费看| 久久精品国产亚洲网站| 亚洲国产成人一精品久久久| 别揉我奶头 嗯啊视频| 黄色一级大片看看| 亚洲欧美日韩东京热| 亚洲成人中文字幕在线播放| 日本av手机在线免费观看| 晚上一个人看的免费电影| 亚洲欧美日韩东京热| 久久国产乱子免费精品| 一二三四中文在线观看免费高清| 18禁动态无遮挡网站| 国产成人精品久久久久久| 99久久精品热视频| 少妇猛男粗大的猛烈进出视频 | 国产伦精品一区二区三区四那| 最近中文字幕2019免费版| 久久国产乱子免费精品| 精品午夜福利在线看| 亚洲色图av天堂| 国产乱人视频| 色吧在线观看| 国产精品久久久久久av不卡| 啦啦啦啦在线视频资源| 亚洲性久久影院| 91久久精品电影网| 中国三级夫妇交换| 最近中文字幕高清免费大全6| 亚洲欧美日韩另类电影网站 | 亚洲欧美日韩卡通动漫| 麻豆乱淫一区二区| av专区在线播放| 三级男女做爰猛烈吃奶摸视频| av免费观看日本| 免费观看性生交大片5| 亚洲人与动物交配视频| 一级毛片aaaaaa免费看小| 波野结衣二区三区在线| 亚洲自偷自拍三级| 97热精品久久久久久| 中文字幕免费在线视频6| 日本熟妇午夜| 哪个播放器可以免费观看大片| 国产精品久久久久久av不卡| 国产淫片久久久久久久久| 搞女人的毛片| 各种免费的搞黄视频| 51国产日韩欧美| 亚洲精华国产精华液的使用体验| 日韩欧美精品v在线| 在线观看美女被高潮喷水网站| 99久久人妻综合| 肉色欧美久久久久久久蜜桃 | 韩国高清视频一区二区三区| 中国美白少妇内射xxxbb| 男人舔奶头视频| 大片电影免费在线观看免费| 18禁在线播放成人免费| 国产大屁股一区二区在线视频| 欧美精品人与动牲交sv欧美| 久久久国产一区二区| 又爽又黄无遮挡网站| 白带黄色成豆腐渣| 婷婷色综合大香蕉| 能在线免费看毛片的网站| 免费人成在线观看视频色| 男人舔奶头视频| 少妇猛男粗大的猛烈进出视频 | 亚洲欧美精品专区久久| 丰满乱子伦码专区| 乱系列少妇在线播放| 内地一区二区视频在线| 日韩,欧美,国产一区二区三区| 在线看a的网站| 高清av免费在线| 好男人视频免费观看在线| 男女那种视频在线观看| 成人免费观看视频高清| 精品人妻偷拍中文字幕| 97超视频在线观看视频| av国产精品久久久久影院| 国产精品一二三区在线看| 欧美一级a爱片免费观看看| 亚洲精品视频女| 有码 亚洲区| 中文字幕制服av| 国产高清不卡午夜福利| 亚洲内射少妇av| 女人被狂操c到高潮| 网址你懂的国产日韩在线| 成人黄色视频免费在线看| 国语对白做爰xxxⅹ性视频网站| 最近最新中文字幕大全电影3| 中文天堂在线官网| 亚洲aⅴ乱码一区二区在线播放| 国产一区有黄有色的免费视频| 国产午夜精品久久久久久一区二区三区| 欧美成人一区二区免费高清观看| 青春草视频在线免费观看| 美女主播在线视频| 99久久精品国产国产毛片| 国产又色又爽无遮挡免| 亚洲精华国产精华液的使用体验| 最新中文字幕久久久久| 亚洲四区av| 亚洲无线观看免费| 午夜激情久久久久久久| 建设人人有责人人尽责人人享有的 | 午夜福利在线观看免费完整高清在| 天美传媒精品一区二区| 午夜精品一区二区三区免费看| 天天躁夜夜躁狠狠久久av| 国产av国产精品国产| 亚洲在线观看片| 国产 一区精品| 春色校园在线视频观看| 久久久精品免费免费高清| 国产真实伦视频高清在线观看| 亚洲高清免费不卡视频| 午夜精品国产一区二区电影 | 肉色欧美久久久久久久蜜桃 | 久久精品久久久久久久性| 黄色配什么色好看| 国内少妇人妻偷人精品xxx网站| 听说在线观看完整版免费高清| 国产精品伦人一区二区| 国产午夜精品一二区理论片| 哪个播放器可以免费观看大片| 亚洲va在线va天堂va国产| 欧美极品一区二区三区四区| 肉色欧美久久久久久久蜜桃 | 尤物成人国产欧美一区二区三区| 日韩欧美 国产精品| 成人漫画全彩无遮挡| av黄色大香蕉| 国产爱豆传媒在线观看| 舔av片在线| 午夜精品国产一区二区电影 | 美女被艹到高潮喷水动态| 午夜福利高清视频| 亚洲高清免费不卡视频| 一级爰片在线观看| 亚洲欧美日韩无卡精品| 毛片女人毛片| 日韩国内少妇激情av| 国产黄a三级三级三级人| 欧美日韩亚洲高清精品| 国产永久视频网站| 日日摸夜夜添夜夜爱| 亚洲欧洲国产日韩| 一级黄片播放器| 亚洲av日韩在线播放| 一级毛片电影观看| 久久精品夜色国产| 欧美最新免费一区二区三区| 亚洲一级一片aⅴ在线观看| 亚洲最大成人手机在线| 可以在线观看毛片的网站| 亚洲精品自拍成人| 免费大片18禁| 国产精品99久久99久久久不卡 | 精品人妻视频免费看| 亚洲av免费高清在线观看| 亚洲最大成人中文| 免费不卡的大黄色大毛片视频在线观看| 日韩 亚洲 欧美在线| 九九在线视频观看精品| 日本av手机在线免费观看| 老司机影院成人| 国产精品熟女久久久久浪| 中文字幕久久专区| 中文字幕av成人在线电影| 亚洲综合精品二区| 在线精品无人区一区二区三 | 日韩国内少妇激情av| 大片电影免费在线观看免费| 女的被弄到高潮叫床怎么办| 熟女av电影| 国产成人免费观看mmmm| 美女高潮的动态| 久热这里只有精品99| 观看美女的网站| 新久久久久国产一级毛片| 3wmmmm亚洲av在线观看| 麻豆国产97在线/欧美| 中文在线观看免费www的网站| 一区二区三区免费毛片| 色婷婷久久久亚洲欧美| 免费看光身美女| 99久久精品热视频| 亚洲成人久久爱视频| 狠狠精品人妻久久久久久综合| 99re6热这里在线精品视频| 亚洲欧美日韩东京热| 少妇猛男粗大的猛烈进出视频 | 高清av免费在线| 婷婷色av中文字幕| 最近中文字幕2019免费版| 国产成人精品婷婷| 久久精品国产a三级三级三级| 久久女婷五月综合色啪小说 | 欧美三级亚洲精品| 欧美激情在线99| 99热这里只有精品一区| 欧美三级亚洲精品| 国产精品一及| 最近2019中文字幕mv第一页| 色5月婷婷丁香| 日韩大片免费观看网站| 日韩成人av中文字幕在线观看| 联通29元200g的流量卡| av在线观看视频网站免费| 成人国产av品久久久| 日韩av在线免费看完整版不卡| 草草在线视频免费看| 国产成人免费无遮挡视频| 一个人观看的视频www高清免费观看| 成人二区视频| 免费av不卡在线播放| 国产成人免费观看mmmm| 少妇熟女欧美另类| 日本熟妇午夜| 在线观看三级黄色| 天天躁日日操中文字幕| 亚洲av不卡在线观看| 黄片wwwwww| 国产精品一区www在线观看| 亚洲国产高清在线一区二区三| 一级毛片久久久久久久久女| 少妇的逼水好多| 少妇丰满av| 中文字幕免费在线视频6| 波野结衣二区三区在线| 免费观看性生交大片5| 日韩 亚洲 欧美在线| 少妇的逼水好多| 一本色道久久久久久精品综合| 啦啦啦啦在线视频资源| 人人妻人人爽人人添夜夜欢视频 | 亚洲av中文av极速乱| 国产亚洲av嫩草精品影院| 国产精品国产av在线观看| 精品人妻视频免费看| 青春草国产在线视频| 在线播放无遮挡| 日韩伦理黄色片| 九色成人免费人妻av| 亚洲精品乱码久久久v下载方式| 两个人的视频大全免费| a级毛片免费高清观看在线播放| 激情 狠狠 欧美| tube8黄色片| 熟女电影av网| 六月丁香七月| 欧美变态另类bdsm刘玥| 久久人人爽人人爽人人片va| 超碰97精品在线观看| 日韩av在线免费看完整版不卡| 一级毛片久久久久久久久女| 午夜视频国产福利| 国产乱人视频| 丝袜脚勾引网站| 国产伦精品一区二区三区四那| 日韩制服骚丝袜av| 国产亚洲最大av| 一级爰片在线观看| av在线蜜桃| 午夜激情久久久久久久| 一级黄片播放器| 国产精品.久久久| 黄色怎么调成土黄色| 在线亚洲精品国产二区图片欧美 | 日韩免费高清中文字幕av| 白带黄色成豆腐渣| 美女内射精品一级片tv| 岛国毛片在线播放| 日日摸夜夜添夜夜爱| 国产色婷婷99| 国产极品天堂在线| 黄色视频在线播放观看不卡| 日韩人妻高清精品专区| 亚洲欧美日韩另类电影网站 | 日本猛色少妇xxxxx猛交久久| 新久久久久国产一级毛片| 精品国产一区二区三区久久久樱花 | 一级片'在线观看视频| 国产精品久久久久久久久免| 岛国毛片在线播放| 亚洲美女视频黄频| 亚洲久久久久久中文字幕| 亚洲欧美一区二区三区黑人 | 中文字幕制服av| 人体艺术视频欧美日本| 国精品久久久久久国模美|