• <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.

    麻豆一二三区av精品| 国产精品九九99| 国产精品久久久久久亚洲av鲁大| www日本在线高清视频| 人妻久久中文字幕网| 午夜激情欧美在线| 亚洲五月婷婷丁香| 69av精品久久久久久| 久久精品综合一区二区三区| 色吧在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 高清在线国产一区| 99热6这里只有精品| 久久九九热精品免费| 91久久精品国产一区二区成人 | 国产精品亚洲美女久久久| 国产日本99.免费观看| 麻豆国产av国片精品| 国产蜜桃级精品一区二区三区| 国产av在哪里看| 国产乱人伦免费视频| 欧美最黄视频在线播放免费| 精品久久久久久成人av| 日日干狠狠操夜夜爽| 琪琪午夜伦伦电影理论片6080| 国产免费av片在线观看野外av| 国产在线精品亚洲第一网站| 国产单亲对白刺激| 久久精品91无色码中文字幕| 亚洲av第一区精品v没综合| 精品久久久久久久毛片微露脸| 国产精品一及| 精品国产三级普通话版| 一a级毛片在线观看| 亚洲精品一区av在线观看| 亚洲 国产 在线| 亚洲熟女毛片儿| 一级黄色大片毛片| 欧美精品啪啪一区二区三区| 9191精品国产免费久久| 日日摸夜夜添夜夜添小说| 国产精品永久免费网站| 国产97色在线日韩免费| bbb黄色大片| a在线观看视频网站| 观看美女的网站| 1024香蕉在线观看| 亚洲精品一卡2卡三卡4卡5卡| 老司机深夜福利视频在线观看| 国产精品野战在线观看| 欧美日韩一级在线毛片| 男人和女人高潮做爰伦理| 老鸭窝网址在线观看| 亚洲av第一区精品v没综合| 琪琪午夜伦伦电影理论片6080| 岛国视频午夜一区免费看| 少妇人妻一区二区三区视频| 久久久久亚洲av毛片大全| 精品一区二区三区视频在线观看免费| 亚洲电影在线观看av| 日本黄色视频三级网站网址| 免费观看精品视频网站| 久久午夜综合久久蜜桃| 久久精品影院6| 丰满人妻熟妇乱又伦精品不卡| av片东京热男人的天堂| 在线观看日韩欧美| 禁无遮挡网站| 亚洲无线在线观看| 亚洲中文av在线| 国产黄a三级三级三级人| 亚洲熟妇中文字幕五十中出| 黄色女人牲交| 亚洲精品一区av在线观看| 色在线成人网| 国产男靠女视频免费网站| 不卡一级毛片| 国产亚洲精品久久久久久毛片| 亚洲自偷自拍图片 自拍| 又黄又粗又硬又大视频| 亚洲在线自拍视频| 级片在线观看| 精品一区二区三区av网在线观看| 亚洲欧美日韩高清在线视频| 久久久精品大字幕| 精品久久久久久久人妻蜜臀av| 可以在线观看毛片的网站| 国产亚洲精品综合一区在线观看| 欧美高清成人免费视频www| 黑人操中国人逼视频| www日本在线高清视频| 国产精品98久久久久久宅男小说| 国产精品一区二区精品视频观看| 变态另类成人亚洲欧美熟女| 国产精品一区二区三区四区免费观看 | 性色av乱码一区二区三区2| or卡值多少钱| 欧美性猛交黑人性爽| 亚洲性夜色夜夜综合| 亚洲性夜色夜夜综合| 日韩欧美免费精品| 偷拍熟女少妇极品色| 99国产综合亚洲精品| 午夜亚洲福利在线播放| 狠狠狠狠99中文字幕| 国产精品 国内视频| 久9热在线精品视频| 别揉我奶头~嗯~啊~动态视频| 一进一出抽搐动态| 免费在线观看亚洲国产| e午夜精品久久久久久久| www国产在线视频色| 国产69精品久久久久777片 | 亚洲电影在线观看av| 波多野结衣高清作品| 欧美午夜高清在线| 国内揄拍国产精品人妻在线| 宅男免费午夜| 欧美日韩中文字幕国产精品一区二区三区| 少妇裸体淫交视频免费看高清| 亚洲黑人精品在线| 麻豆国产av国片精品| 久久午夜亚洲精品久久| 丰满的人妻完整版| 国产伦精品一区二区三区四那| 97超级碰碰碰精品色视频在线观看| 欧美激情在线99| 精品乱码久久久久久99久播| 精品国产三级普通话版| 久久久久亚洲av毛片大全| 亚洲av五月六月丁香网| 天堂网av新在线| 成人鲁丝片一二三区免费| 中亚洲国语对白在线视频| 丝袜人妻中文字幕| 中国美女看黄片| 成人性生交大片免费视频hd| 国产麻豆成人av免费视频| 99热6这里只有精品| 可以在线观看的亚洲视频| 欧美丝袜亚洲另类 | 搡老岳熟女国产| 男女下面进入的视频免费午夜| 九九久久精品国产亚洲av麻豆 | 久久久久久久精品吃奶| 亚洲精品久久国产高清桃花| 成人性生交大片免费视频hd| 久久久国产成人精品二区| 国产成人aa在线观看| 97超级碰碰碰精品色视频在线观看| svipshipincom国产片| 99久久精品一区二区三区| 久久久久免费精品人妻一区二区| 成人国产综合亚洲| 不卡一级毛片| 亚洲乱码一区二区免费版| 精品国产三级普通话版| 久久精品影院6| 香蕉丝袜av| 在线观看午夜福利视频| 亚洲七黄色美女视频| 国产成人影院久久av| 精品一区二区三区视频在线 | 亚洲欧美日韩卡通动漫| 欧美黑人巨大hd| 国产亚洲精品久久久com| e午夜精品久久久久久久| 欧美日韩亚洲国产一区二区在线观看| 国产在线精品亚洲第一网站| 一本综合久久免费| 在线观看免费视频日本深夜| 婷婷精品国产亚洲av在线| 国产免费av片在线观看野外av| 琪琪午夜伦伦电影理论片6080| 国产成人精品无人区| 青草久久国产| 精品欧美国产一区二区三| 久久这里只有精品19| 18禁黄网站禁片午夜丰满| 99热这里只有精品一区 | 日韩欧美在线乱码| 免费在线观看视频国产中文字幕亚洲| 国语自产精品视频在线第100页| 1000部很黄的大片| 欧美日韩综合久久久久久 | 日本 av在线| 天堂动漫精品| 最近视频中文字幕2019在线8| 天天躁日日操中文字幕| 熟女人妻精品中文字幕| 99久久综合精品五月天人人| 国产麻豆成人av免费视频| 熟女人妻精品中文字幕| 99在线视频只有这里精品首页| www国产在线视频色| 国产在线精品亚洲第一网站| 中文字幕久久专区| 亚洲最大成人中文| 身体一侧抽搐| 天天躁狠狠躁夜夜躁狠狠躁| 欧美黄色片欧美黄色片| 91av网站免费观看| 国产亚洲精品久久久com| 亚洲精品456在线播放app | 国产成+人综合+亚洲专区| 俄罗斯特黄特色一大片| 精品福利观看| 久久伊人香网站| 亚洲av美国av| 大型黄色视频在线免费观看| 亚洲国产中文字幕在线视频| 狠狠狠狠99中文字幕| 五月伊人婷婷丁香| 欧美另类亚洲清纯唯美| 麻豆国产97在线/欧美| 亚洲av中文字字幕乱码综合| 99热这里只有是精品50| 国产伦精品一区二区三区四那| 精品久久久久久成人av| 制服人妻中文乱码| 黑人操中国人逼视频| 国产激情偷乱视频一区二区| 亚洲第一电影网av| 亚洲国产欧美一区二区综合| 久久精品夜夜夜夜夜久久蜜豆| 日本成人三级电影网站| 国产亚洲精品久久久久久毛片| 日韩免费av在线播放| 精品午夜福利视频在线观看一区| 色尼玛亚洲综合影院| 亚洲avbb在线观看| 成人特级黄色片久久久久久久| 国产爱豆传媒在线观看| 亚洲午夜精品一区,二区,三区| 国产一级毛片七仙女欲春2| 老汉色∧v一级毛片| 久久久久久人人人人人| 国产主播在线观看一区二区| 九九久久精品国产亚洲av麻豆 | 美女 人体艺术 gogo| 午夜激情欧美在线| 中文字幕精品亚洲无线码一区| 午夜免费激情av| 在线观看免费午夜福利视频| 欧美国产日韩亚洲一区| 日本在线视频免费播放| 成人亚洲精品av一区二区| 国产精品野战在线观看| 小蜜桃在线观看免费完整版高清| 日韩欧美在线乱码| av天堂中文字幕网| 国内精品美女久久久久久| 国内精品美女久久久久久| 男女下面进入的视频免费午夜| 国内揄拍国产精品人妻在线| 国产又色又爽无遮挡免费看| 巨乳人妻的诱惑在线观看| 亚洲熟妇中文字幕五十中出| 欧美性猛交╳xxx乱大交人| 黑人欧美特级aaaaaa片| 看黄色毛片网站| 成人国产综合亚洲| 日本成人三级电影网站| 国产高清三级在线| 亚洲精品乱码久久久v下载方式 | 人人妻,人人澡人人爽秒播| 日韩 欧美 亚洲 中文字幕| 草草在线视频免费看| 久久99热这里只有精品18| 久久久久久久午夜电影| 一夜夜www| 国产精品亚洲一级av第二区| 每晚都被弄得嗷嗷叫到高潮| 黄片小视频在线播放| 免费人成视频x8x8入口观看| 国产激情久久老熟女| 成年版毛片免费区| 国产亚洲av高清不卡| 99在线人妻在线中文字幕| 亚洲人成电影免费在线| 欧美国产日韩亚洲一区| 国产精品99久久99久久久不卡| 97超视频在线观看视频| 母亲3免费完整高清在线观看| 国产精品亚洲av一区麻豆| 91老司机精品| 欧美日韩精品网址| 男女下面进入的视频免费午夜| 久久亚洲真实| 在线看三级毛片| 老司机深夜福利视频在线观看| 不卡一级毛片| 成人特级黄色片久久久久久久| 久久久久久大精品| 99国产精品一区二区蜜桃av| 高清在线国产一区| 亚洲午夜精品一区,二区,三区| 国产私拍福利视频在线观看| 成年人黄色毛片网站| 精品国产乱子伦一区二区三区| 国内精品久久久久精免费| 日本成人三级电影网站| 在线免费观看的www视频| 国产精品美女特级片免费视频播放器 | 国产精品久久视频播放| 国产av一区在线观看免费| 亚洲人成网站在线播放欧美日韩| 精华霜和精华液先用哪个| 色精品久久人妻99蜜桃| 国产亚洲精品久久久com| av黄色大香蕉| 国产三级中文精品| 免费人成视频x8x8入口观看| 男女那种视频在线观看| 欧美日韩瑟瑟在线播放| 一级a爱片免费观看的视频| 亚洲黑人精品在线| 观看免费一级毛片| 午夜视频精品福利| 露出奶头的视频| 国产视频内射| 中文资源天堂在线| 人妻夜夜爽99麻豆av| 十八禁网站免费在线| 两个人视频免费观看高清| 午夜免费激情av| 免费一级毛片在线播放高清视频| 亚洲av成人不卡在线观看播放网| 啪啪无遮挡十八禁网站| www.精华液| 日韩欧美 国产精品| 亚洲一区二区三区不卡视频| 国产极品精品免费视频能看的| 中文字幕熟女人妻在线| 国产久久久一区二区三区| 国语自产精品视频在线第100页| 18禁国产床啪视频网站| 99久久无色码亚洲精品果冻| 欧美色欧美亚洲另类二区| 一二三四在线观看免费中文在| 亚洲精品色激情综合| 欧美日韩福利视频一区二区| netflix在线观看网站| 五月伊人婷婷丁香| 久久久国产成人免费| 岛国在线观看网站| 免费电影在线观看免费观看| 成人av在线播放网站| 欧美日韩精品网址| 精品久久久久久,| 欧美乱妇无乱码| 成人一区二区视频在线观看| 真人做人爱边吃奶动态| 亚洲国产精品sss在线观看| cao死你这个sao货| 成人鲁丝片一二三区免费| 在线十欧美十亚洲十日本专区| 婷婷六月久久综合丁香| 国产美女午夜福利| 国产精品99久久久久久久久| 18禁国产床啪视频网站| 国产精品亚洲av一区麻豆| 又紧又爽又黄一区二区| 亚洲精品国产精品久久久不卡| 国产精品久久电影中文字幕| 久久精品91蜜桃| 免费av毛片视频| 亚洲成av人片免费观看| 欧美日韩国产亚洲二区| 88av欧美| 亚洲性夜色夜夜综合| 级片在线观看| 最好的美女福利视频网| 最近最新中文字幕大全免费视频| 免费看日本二区| 一区二区三区高清视频在线| 国产精品免费一区二区三区在线| 国产视频内射| 亚洲中文字幕一区二区三区有码在线看 | 麻豆成人午夜福利视频| 手机成人av网站| 亚洲国产欧美人成| 日韩欧美国产在线观看| 日本与韩国留学比较| 亚洲精华国产精华精| 香蕉久久夜色| 亚洲美女视频黄频| 成年免费大片在线观看| 国产黄片美女视频| 19禁男女啪啪无遮挡网站| 国产成人av教育| 欧美乱码精品一区二区三区| 又粗又爽又猛毛片免费看| 床上黄色一级片| x7x7x7水蜜桃| 9191精品国产免费久久| 欧美三级亚洲精品| 免费av毛片视频| 香蕉丝袜av| 久久热在线av| av欧美777| 狂野欧美白嫩少妇大欣赏| 亚洲午夜精品一区,二区,三区| 美女高潮喷水抽搐中文字幕| 最近最新免费中文字幕在线| 国产久久久一区二区三区| 白带黄色成豆腐渣| 国产精品一区二区三区四区久久| 久久久久九九精品影院| 久久久精品大字幕| 日韩欧美国产一区二区入口| 精品久久久久久久末码| 看片在线看免费视频| 国产真人三级小视频在线观看| 毛片女人毛片| 99精品在免费线老司机午夜| 无遮挡黄片免费观看| 一个人观看的视频www高清免费观看 | 黑人欧美特级aaaaaa片| 少妇人妻一区二区三区视频| 69av精品久久久久久| 日本一本二区三区精品| 五月伊人婷婷丁香| 国产成人av教育| 夜夜看夜夜爽夜夜摸| 色播亚洲综合网| 搡老熟女国产l中国老女人| 国产午夜精品久久久久久| 高清在线国产一区| 美女免费视频网站| 在线免费观看不下载黄p国产 | av天堂在线播放| 又大又爽又粗| 老熟妇乱子伦视频在线观看| 亚洲精品久久国产高清桃花| 免费在线观看影片大全网站| 免费在线观看亚洲国产| 国产三级在线视频| 蜜桃久久精品国产亚洲av| 美女黄网站色视频| 脱女人内裤的视频| or卡值多少钱| 亚洲av片天天在线观看| а√天堂www在线а√下载| 国产视频内射| 成人国产一区最新在线观看| 中文字幕精品亚洲无线码一区| 精品久久久久久,| 久久精品国产清高在天天线| 最新美女视频免费是黄的| 日韩三级视频一区二区三区| 99精品在免费线老司机午夜| 亚洲无线在线观看| 国产免费男女视频| 国产精品久久视频播放| 1000部很黄的大片| av女优亚洲男人天堂 | 无人区码免费观看不卡| 亚洲欧美日韩无卡精品| 国产日本99.免费观看| 成人三级黄色视频| 色综合欧美亚洲国产小说| 美女免费视频网站| 亚洲第一电影网av| 变态另类丝袜制服| 国产精品亚洲美女久久久| 天天躁日日操中文字幕| 亚洲 欧美 日韩 在线 免费| 一级作爱视频免费观看| 亚洲av成人不卡在线观看播放网| 国产成人欧美在线观看| 久久久成人免费电影| 久久午夜综合久久蜜桃| 国产精品九九99| 99久久99久久久精品蜜桃| 亚洲七黄色美女视频| 婷婷精品国产亚洲av| 色精品久久人妻99蜜桃| 国内精品美女久久久久久| 国产成人一区二区三区免费视频网站| 97人妻精品一区二区三区麻豆| 日本三级黄在线观看| 国产欧美日韩精品亚洲av| 好看av亚洲va欧美ⅴa在| 香蕉丝袜av| 国产黄片美女视频| 高清毛片免费观看视频网站| 日本黄色视频三级网站网址| 熟女人妻精品中文字幕| 午夜影院日韩av| 99久久成人亚洲精品观看| 久久中文看片网| 99久久精品一区二区三区| а√天堂www在线а√下载| 夜夜看夜夜爽夜夜摸| 中文字幕精品亚洲无线码一区| 色噜噜av男人的天堂激情| 国产 一区 欧美 日韩| 97超视频在线观看视频| 亚洲国产日韩欧美精品在线观看 | 午夜激情欧美在线| 午夜福利欧美成人| 国产又黄又爽又无遮挡在线| 床上黄色一级片| 日本熟妇午夜| www.www免费av| 精品国产超薄肉色丝袜足j| 在线免费观看不下载黄p国产 | 天堂网av新在线| 国产三级黄色录像| 国产精品国产高清国产av| 久久久色成人| 91字幕亚洲| 熟女人妻精品中文字幕| 久久久久亚洲av毛片大全| xxxwww97欧美| 啦啦啦韩国在线观看视频| 三级国产精品欧美在线观看 | 给我免费播放毛片高清在线观看| 一区二区三区激情视频| 香蕉国产在线看| av在线蜜桃| 亚洲人与动物交配视频| 嫩草影视91久久| 国产精品亚洲av一区麻豆| 男女下面进入的视频免费午夜| www.精华液| 一进一出好大好爽视频| 中文字幕最新亚洲高清| 五月玫瑰六月丁香| 国产探花在线观看一区二区| 丰满人妻一区二区三区视频av | 国产视频内射| 美女午夜性视频免费| 色av中文字幕| 五月伊人婷婷丁香| www.精华液| 丰满人妻一区二区三区视频av | 精品福利观看| 亚洲狠狠婷婷综合久久图片| 免费在线观看日本一区| 男人舔女人下体高潮全视频| av黄色大香蕉| 变态另类丝袜制服| 久久久久久久久中文| 免费观看人在逋| 免费av不卡在线播放| 亚洲欧美日韩高清专用| 老汉色∧v一级毛片| x7x7x7水蜜桃| 久久人人精品亚洲av| 99久国产av精品| 亚洲片人在线观看| 亚洲男人的天堂狠狠| 一二三四在线观看免费中文在| 免费看十八禁软件| 三级国产精品欧美在线观看 | 国产三级中文精品| 99久国产av精品| 丰满人妻一区二区三区视频av | 怎么达到女性高潮| 欧美一区二区国产精品久久精品| 久久久久久久久中文| 国产黄色小视频在线观看| 国产综合懂色| 宅男免费午夜| 日本一本二区三区精品| 天堂√8在线中文| 亚洲中文日韩欧美视频| 听说在线观看完整版免费高清| 免费在线观看亚洲国产| 精品一区二区三区av网在线观看| 精华霜和精华液先用哪个| 一卡2卡三卡四卡精品乱码亚洲| 一个人看的www免费观看视频| a在线观看视频网站| 伊人久久大香线蕉亚洲五| 国产激情久久老熟女| 久久精品影院6| 亚洲自偷自拍图片 自拍| 老鸭窝网址在线观看| 一进一出抽搐gif免费好疼| 色精品久久人妻99蜜桃| 女同久久另类99精品国产91| 叶爱在线成人免费视频播放| 亚洲色图av天堂| 黄片大片在线免费观看| 亚洲精品在线美女| 99精品久久久久人妻精品| 亚洲精品在线观看二区| 欧美一级毛片孕妇| 国产蜜桃级精品一区二区三区| 国产成人av激情在线播放| 久久久久久大精品| 757午夜福利合集在线观看| 日本一二三区视频观看| 国内少妇人妻偷人精品xxx网站 | 又紧又爽又黄一区二区| 麻豆成人午夜福利视频| 免费看a级黄色片| 国产成人啪精品午夜网站| 一二三四社区在线视频社区8| 小蜜桃在线观看免费完整版高清| 制服丝袜大香蕉在线| 免费av不卡在线播放| 麻豆成人午夜福利视频| 亚洲成人精品中文字幕电影| 黄色丝袜av网址大全| 香蕉av资源在线| 可以在线观看毛片的网站| 听说在线观看完整版免费高清| 老熟妇乱子伦视频在线观看| 久久午夜亚洲精品久久| 婷婷精品国产亚洲av在线| 亚洲av电影不卡..在线观看|