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

    Development of a doppler asymmetric spatial heterodyne interferometer for ground-based wind field detection at the 557.7 nm wavelength

    2023-10-07 07:40:50LIUHuanJIANGLunZHANGXiaofeiFUYunSONGYansongTONGShoufengLIUXianzhu
    中國光學 2023年5期

    LIU Huan,JIANG Lun,2,3*,ZHANG Xiao-fei,F(xiàn)U Yun,SONG Yan-song,2,3,TONG Shou-feng,2,3,LIU Xian-zhu,2,3

    (1.School of Opto-Electronic Engineering, Changchun University of Science and Technology,Changchun 130022, China;2.Peng Cheng Laboratory, Shenzhen 518000, China;3.Key Laboratory of Fundamental Science for National Defense of Aero and Ground Laser Communication Technology, Changchun University of Science and Technology, Changchun 130022, China)

    * Corresponding author,E-mail: jlciomp@163.com

    Abstract: A ground-based Doppler Asymmetric Spatial Heterodyne (DASH) interferometer with a high Signal-to-Noise Ratio (SNR) and large etendue (A?) with thermal compensation was developed to detect wind field information in the middle atmosphere.The detailed parameters and index of the DASH interferometer were developed for the 557.7 nm oxygen airglow spectral line.The system was designed with an expanded Field Of View (FOV) and thermal compensation.The half-FOV angle reached 2.815°, the etendue was 0.095 25 cm2sr, and the system’s SNR was approximately 113.75.Through the thermal compensation design,the final optical path difference with temperature variation (dΔd0/dT) was only 2.224×10-7mm/°C.The optical system was designed and optimized according to the corresponding parameters.Image-side telecentric and bilateral telecentric optical system structures were used in the entrance optics and exit optics, respectively,and parameters such as telecentricity and distortion met the detection requirements.To verify the design results, a ground-based DASH interferometer experimental platform was constructed, and indoor and outdoor ground-based experiments were conducted.In the final experiment, clear interference fringes were obtained,which proves that the system design results of the DASH interferometer are reasonable, and the system’s SNR and etendue meet the detection requirements.

    Key words: ground-based DASH interferometer; 557.7 nm oxygen atom airglow spectral line; optical design; signal-to-noise ratio

    1 Introduction

    The wind field in the middle and upper atmosphere affects climate change, the flight trajectory parameters of spacecraft, and the precise guidance of military missiles.Therefore, it is necessary to research the accurate and efficient detection of windfield information.Over the last century, the Fabry-Perot interferometer has been regarded as a typical passive wind-measuring device owing to its high sensitivity and small temperature influence on the optical path difference of the interferometer.However, the Fabry-Perot étalon has high accuracy requirements and a difficult manufacturing process.Additionally, its Field Of View (FOV) is not easily broadened (usually less than 2°)[1], resulting in a low Signal-to-Noise Ratio (SNR).The Michelson interferometer is another device commonly used for passive wind measurements and it tracks the Doppler frequency shift by monitoring the phase shift of multiple interference fringes.FOV broadening technology can be used to expand the FOV.However,dynamic measurement methods reduce the measurement accuracy of the wind speed[2].

    Doppler Asymmetric Spatial Heterodyne (DASH) technology has both grating diffraction and spatial interference capabilities.In the optical path, an FOV prism can be used to increase the FOV of the system to obtain a large etendue and SNR.Moreover, DASH technology has a series of advantages,such as high spectral resolution and simultaneous multi-spectral line static detection.Since 2006,DASH has been used to detect middle and upper atmospheric wind fields[3].The DASH interferometer takes the airglow spectral line in the middle and upper atmosphere as the light source, inverts the windfield parameters via the Doppler effect of the spectral line, and finally obtains the wind speed using the interference absolute phase difference.Owing to its high resolution and simple structure, the DASH system has clear advantages in detecting atmospheric spectral Doppler shifts, which has attracted the attention of many research teams.

    Englert[3]proposed the concept of DASH technology to detect the wind field in the middle and upper atmosphere in 2006.They used a breadboard to detect two lasers that simultaneously changed the voltage and found that they caused the same trend in interference phase changes.The feasibility of DASH technology for multispectral line observation and real-time phase tracking was subsequently verified[3].The following year, with the support of the US Naval Laboratory, the first breadboard structure of a DASH interferometer without an expanded FOV was constructed[4].In 2010, a ground-based (Redline DASH Demonstration, REDDI)[5]and airborne(Redline Interferometer for Doppler Winds, ARROW)[6]principle prototype was successfully developed for the oxygen red line with a wavelength of 630 nm in the atmosphere, and a spectral line with a wavelength of 630.48 nm emitted by an Ne lamp was used for calibration and real-time phase calibration.Launched on the Ionospheric Connection Explorer (ICON) satellite in October 2019, the Michelson Interferometer for Global High-resolution Thermospheric Imaging (MIGHTI) developed by the team can remotely measure atmospheric wind and temperature fields within a range of 90-300 km from Earth by detecting the naturally occurring red(630.0 nm) and green (557.7 nm) lines of oxygen atoms and oxygen molecules (762 nm) near airglow spectral lines[7-8].

    By building the DASH interferometer platform,New York University conducted a wind-field simulation experiment with a 630 nm Ne laser and obtained wind speed data by preliminary inversion[9].Using this, the SWIFT-DASH principle prototype for stratospheric wind-field and infrared ozone detection was developed, which opened a new avenue for expanding infrared ozone detection[10].

    The Anhui Institute of Optics and Fine Mechanics developed a large-offset DASH interferometer (LODI) based on simulations and manual system design.The characteristic emission lines emitted by a potassium lamp with wavelengths of 769.90 nm and 766.49 nm were used for indoor simulation experiments[11].Another DASH system was used to perform ground-based measurement experiments on the Doppler shift of the solar absorption line[12].

    From 2011, the Xi’an Institute of Optics and Precision Mechanics, Chinese Academy of Sciences,started using DASH technology to detect the atmospheric wind field including the single-arm-type structure of the Dual Field Of View (DFOV) Doppler heterodyne interferometer optical system simulation design.Two-FOV-directions interference figure at different levels can be obtained without a field scanning device[13].A Doppler differential interferometer measuring 630 nm oxygen red line was designed[14].

    For the detection of middle atmosphere windfield information, there is a need for middle atmosphere consensus for global network observation.Ground-based detection technology has good local time cover features, and the DASH system is a relatively cheap and sophisticated example of passive wind equipment.The establishment of ground-based DASH interferometers will certainly lay the foundation for spaceborne exploration and space-physics research on a global scale.

    In this paper, a ground-based DASH interferometer based on the 557.7 nm oxygen airglow is proposed to measure the wind field in the middle atmosphere.Owing to the high sensitivity of the DASH interferometer, clouds and air between the airglow and the ground can significantly affect its detection accuracy.Therefore, the system requires a high SNR.In Section 2, the structural principle of the ground-based DASH interferometer is introduced.As the core component of the system, the interference module was designed in detail for its optimal optical path difference, FOV extension, and thermal compensation, and the SNR of the system was calculated.Next, the design process of the entrance optics and exit optics is described.In Section 3, the indoor and outdoor ground-based experiments and the experimental results produced by the DASH system are described.Finally, in Section 4, the DASH system and experimental results are summarized.

    2 Overall parameter design

    2.1 Structural principle

    The DASH interferometer is based on the structure of the Michelson interferometer, where the mirrors at the ends of the two optical paths are replaced by gratings placed at Littrow angles (θL) to the orthogonal plane of the optical axis.After reaching the interferometer, light experiences grating diffraction and returns with different frequencies and the optical axis to form different angles, which separates different frequencies of light in space.The spectral line with a specified frequency is extracted using a Fourier transform, the phase is calculated,and the wind speed is retrieved.The two optical paths of the DASH interferometer are asymmetrically designed, and the light passes through the two optical paths to form an optical path difference of 2Δdto increase interference phase sensitivity.

    As shown in Figure 1, the airglow-emission light with a given Doppler shift is collimated into the interferometer by the entrance optics.The incident wavefront is divided into two optical paths of equal intensity by the beam splitter, and the light reaches the grating through the extended field prism of the two optical paths.After diffraction occurs on the grating surface, the light returns to the beam splitter on the original path and forms interference fringes at the interferometer exit.The interference fringes are scaled using the detector imaging system and imaged on the detector.The asymmetric length ?tof the beam splitter, the refractive index of the beam splitter material, and the length difference of the air spacers between the two optical paths of the interferometer jointly provide the optimal optical path difference of the interferometer.The final interference equation is[11]

    whereB(σ) is the incident spectral radiance,σis the arbitrary wave number in the band,σLis the fundamental frequency wave number of the system,θLis the Littrow angle of the grating,xis the position coordinate of the detector, and ?dis the bias of one optical path of the interferometer.The optimal optical path differenceLopt= 2?d.

    2.2 Interferometer design

    The interferometer is the core device of the system, and its design quality also determines the accuracy of the wind-speed detection.The following is a theoretical analysis of the design of each part of the interferometer.

    2.2.1 Optimal optical path difference

    The target airglow wave numberσ0changes toσ'result from the Doppler shift caused by the wind speedv.That is,σ'=(1-v/c)σ0, wherecis the speed of light.The relationship between the phase changeδ?corresponding to the frequency shift and the wind speedvis as follows:

    Therefore, the phase sensitivity can be increased by increasing the optical path differenceL;however, the detection target is the airglow spectral line with a given linewidth, and the spectral linewidth caused by temperature is[15]

    whereTis the temperature, andMis the molar mass of the emission source particle.Therefore, for the DASH interferometer, the final interferogram light intensity distribution is

    Because the intensity distribution (Equation(4)) of the interferogram generated by the expanded spectral line is an enveloped cosine function, increasing the optical path difference leads to a decrease in the interference modulation regime; that is,the interference signal becomes weaker.Therefore,there must be an optimal optical path difference that achieves the optimal sensitivity and interference modulation regime in the interference phase.

    The system provides optical paths difference with different air thicknesses, different glass thicknesses and different material refractive indexes of the two arm optical paths and the refractive index of the glass.The optical paths difference is satisfied.

    In this study, the oxygen atom green line O(1S)557.7 nm airglow with a wavelength of 557.7 nm located 90-110 km above the ground was used as the detection target light source.Unlike aurorae, airglow has a relatively weak radiation intensity, but is not controlled by the geomagnetic field.In addition,the airglow is only slightly affected by solar activity and is relatively uniformly distributed across the globe.The O(1S) 557.7 nm airglow is the strongest optical radiation in the visible band of night airglow,which is generated when oxygen atoms transition from one energy level to another.The transition time is approximately 0.8 s, the peak height of radiation is approximately 95 km, and the thickness is approximately 15 km[16-17].

    It is known that the temperature of the oxygen atom green line O(1S) 557.7 nm airglow is approximately 200 K at an atmospheric altitude of 95 km.Oxygen has a molar mass of 16 g/mol, and the Lopt of the DASH system is calculated as follows:

    The relationship between the modulation degree, interference phase value, and optical path difference of the 557.7 nm (200 K) airglow line at heights of 90-110 km is shown in Figure 2.

    The interference phase difference increases with an increase in the optical path difference, but the degree of interference modulation decreases with an increase in the optical path difference.The intersection of the two curves is the optimal choice of the two parameters, which is consistent with the calculated optimal optical path difference of 3.5 cm.

    To avoid the influence of another green oxygen atom’s weak emission zone of 150-300 km on the wind-field measurement at the measured height,the optimal optical path difference of the oxygen atom green line at this height was simulated such that the atmospheric temperature at the corresponding height was 800 K.The functional relationship between the modulation degree, interference phase value, and optical path difference of the spectral line is shown in Figure 3.

    The optimal optical path difference corresponding to the weak emission of green oxygen atoms at heights of 150-300 km was approximately 2 cm, and the modulation degree of the 557.7 nm spectral line at heights of 150-300 km corresponding to the optical path difference of 3.5 cm was very small.Although the two emission bands were airglows of the same wavelength, the atmospheric temperatures differed, and the optimal optical path difference required by the measurement system was also different; therefore, the weak emission band would not affect the detection of the spectral line from 90-110 km.

    To verify whether the optimal optical path difference drifts under different wind speeds, the interference signals corresponding to wind speeds of 30 m/s、70 m/s and 150 m/s were subtracted from the interference signal corresponding to a wind speed of 0 m/s, and the difference is shown in Figure 4.The wind speed affected the signal amplitude, but the maximum position of the signal amplitude was always the same; that is, the optimal optical path difference did not change with wind speed.

    2.2.2 Expanded field of view prism

    To achieve the SNR and etendue required for wind measurement, DASH systems usually increase etendue by increasing the aperture or FOV of the optical system.The aperture index is limited by the weight of manufactured lens and the production cost.Therefore, We adopted FOV expansion technology.When the DASH interferometer satisfies the mutual restriction relationship between the solid angle ? and the spectral resolution abilityR[18-20],

    The on-axis and off-axis ray phase difference is less than the difference between the π conditions, interference fringe gap within a stripe distance, and the view angle corresponding to off-axis light at the time for the largest view angle.Two prisms are present in the two optical paths of the interferometer.The prisms have a certain vertex angleαthat can compensate for the phase difference at different incidence angles and obtain a large FOV under conditions with a high-contrast fringe.

    The minimum deviation angle of the incident prism surface also satisfies the following formula:

    By optimizing the value of the prism’s top angle, the gap between the phase difference of the off-axis and on-axis FOVs is reduced, and a larger FOV tolerance is obtained[21].

    The prism top angleα= 8.742 4° is deduced by selecting an appropriate prism material and combining it with Snell’s refraction formula.As the target light source for outdoor ground-based wind measurement is 557.7 nm airglow, and the target light source for indoor measurement is 557.0 nm krypton lamp.When the system measures the wind on the outdoor ground, it also needs a krypton lamp(557.0 nm) to correct the system, which needs to distinguish the target light source of these two wavelengths, requiring the system to meet the wavelength resolution of 1.093 3 cm-1(0.034 nm).hence, the resolution ability of the interferometer is

    A? is the product of the aperture of the optical system and the solid angle of the FOV used to characterize the luminous flux of the system.When the aperture diameter of the system is 40 mm, the etendue value is:

    compared with the spaceborne DASH system[7-8],the ground-based system required a larger etendue and a less stringent system weight.Under the premise that the FOV was not different, the system aperture should be appropriately increased to improve the system etendue.

    2.2.3 Thermal compensation

    The DASH interferometer inverts the wind speed by determining the phase difference.Therefore, it requires high phase sensitivity.However, the detection site is usually outdoors or spaceborne, and the observation environment is harsh.The variation in ambient temperature affects the instruments’characteristics.Furthermore, the thermal effect of the interferometer originates from the beam splitter, the thermal effects of the refractive index’s change in the two optical paths’ FOV prism, and the thermal expansion of the wedge spacers at the air spacers.The above thermal effects shift the light propagation angle in the interferometer and change the optical path difference, thus reducing the wind speed measurement’s accuracy.

    As shown in Figure 5,d1andd2are the thicknesses of the air spacers,ωis the incident angle of the incident surface of the interferometer,n1is the refractive index of the beam splitter,t1is the thickness of the asymmetric edge of the beam splitter.The bias ?dprovided by the difference between the beam splitter and the spacer on the two optical paths is[22]

    whereθ0,θ1, andθ2are the angles of light on each optical surface.According to the law of refraction,

    to eliminate the influence of angleωand broaden the field, the second-order term should be set as zero.

    Therefore,

    For the integrated DASH interferometer, the beam splitter, prism, and grating are in contact with spacers composed of optical glass, so the thermal expansion of the spacers also affects ?d.As a result,

    wheren1is the refractive index of the beam splitter,t1is the asymmetric thickness of the beam splitter,d1andd2are the thicknesses of spacers 1 and 2, respectively, andα1,α2, andα3are the thermal expansion coefficients of the material of the beam splitter and spacers 1 and 2, respectively.In order to achieve the thermal compensation effect, materials with small temperature influences were selected, as described in Table 3.The final heat dissipation difference result Δd0/dT=2.224×10-7mm/°C, resulting in a phase shift of 0.05 rad/°C.Compared with the phase shift of other DASH systems of 4.29 rad /°C[23],this system has better thermal stability.

    2.3 Signal-to-noise ratio

    2.3.1 Simulation of system SNR

    The luminance of the oxygen atom green line at 557.7 nm at night is 5 kR[11]according to the formula,

    whereRis the Rayleigh unit to represent the luminance, that is the oxygen atom green line airglow luminance isR=106×5 000,λrefers to the central wavelength of the spectral band, whose unit is μm, h is Planck’s constant (that is, 6.626 070 15×10-34J·s),andcis the speed of light in a vacuum (that is,3×108m/s).To convert brightness from Rayleigh units to radiant brightness units, the following is used:Ls(σ)=1.418×10-10W/cm2/sr.

    For the DASH interferometer, the photoelectron numbere(x) of the signal received by the detector is[24]

    whereks(σ) is the spectral radiation brightness of the target received by a single pixel of the detector per unit of time, represented by the number of photoelectrons[24-26],

    To improve the SNR, it is necessary to use a detector with a refrigeration function.In this study,the EMCCD of Andor Company was selected as the interferometer detector in the ground-based detection experiment.The definitions of the system parameters are as follows.

    The photon noise can be expressed as

    Using Equation (27), the SNR was calculated to be approximately 113.75.The relationship between the SNR of the interferogram and pixel position is presented in Figure 6 (a).Table 1 is the DASH system parameter list.

    Compared with the SNR of other DASH systems[14], this system improved the SNR of the system without affecting the imaging quality by optimizing the optical system and choosing a smaller F #.

    Tab.1 DASH system parameter list

    2.3.2 Influence of SNR on wind speed accuracy

    Noise is one of the important factors in reducing the quality of initial interferogram data.We conducted a simulation analysis on the variation process of wind speed error with SNR.

    Firstly, an interferogram with a wind speed of 100 m/s was simulated, and different white Gaussian noise images were added to simulate the initial interferogram corresponding to the SNR of different values.Then, the inverse Fourier transform of the corresponding frequency spectrum was extracted by the Fourier transform and window function of the initial interferogram.The final wind speed was obtained by inversion.

    In the inversion process, the Nuttall window function of the Gaussian window type with a small sidelobe effect was adopted.When selecting the width of the window function (FWHM, Full width Half Maximum), the spectrum after the Fourier transform should be considered as the convolution of the actual spectrum and the instrument’s linear function, which itself has a given width.Meanwhile,with a smaller window function width, there is a larger sidelobe in the spectrum and a larger distortion region after the convolution of the interferogram.With a wider window function width, there is more noise and a greater wind speed inversion error.In this paper, FWHM=10 and the wind speed obtained without noise inversion is 99.993 6 m/s.

    The final influence curve of the SNR on wind speed accuracy was shown in Figure 6 (b).When the SNR is 0-60, the wind speed error is less than 10-4m/s.When the SNR is larger than 60, the wind speed error is almost 0.The SNR obtained by simulation in this paper is 113.75, so the interference fringe with high contrast can still be detected in a poor environment without a temperature control system.

    2.4 Optical system design

    For ground-based DASH interferometers, the imaging quality of the optical system has a significant impact on the entire DASH interferometer system.The optical system of the DASH interferometer includes entrance and exit optics.The entrance optics are located in front of the interferometer and are used to focus the source to be measured on the grating of the DASH.Exit optics are a transmission optical system with appropriate magnification and image quality, which can transmit the image of the interferogram to the detector.

    2.4.1 Entrance optics

    The magnification of the image-side telecentric optical path does not change with the position of the image plane, and the illumination of the image plane is uniform.To ensure the interference characteristics of the DASH system and uniformity of detector illuminance, an image-side telecentric optical system with a reasonable length was selected for the entrance optics.A narrowband filter was installed in front of the entrance optics to eliminate the signal in the excess spectral area outside the spectral range,and the signal light was transmitted into the interferometer by the entrance optics to be imaged on the gratings.To ensure that the interferometer can receive the target light flux, the entrance pupil diameter of the entrance optics was initially determined to be 40 mm, and the paraxial image height was -7 mm-7 mm.Based on the indicators above, an appropriate initial structure was selected and optimized.

    The Modulation Transfer Function (MTF)curves of the entrance optics are shown in Figure 7.The MTF curve for each FOV was very close to the diffraction limit, and the imaging quality of the system was good.For the 550-560 nm band, the image-side telecentricity of the entrance optics was less than 0.006°, and the distortion was less than 0.05%.The design results met the design requirements, and the total length of the system was 252.8 mm, which is convenient for indoor and external field observation experiments.

    2.4.2 Exit optics

    The exit optics adopt a bilateral telecentric optical system structure.A feature of the bilateral telecentric optical system structure is that the incident parallel light ensures the system has a sufficiently large depth-of-field range, whereas the parallel light leaving the lens ensures that the working distance can vary greatly within the depth-of-field range and the magnification will not change.

    The detector used in the system was a C-type interface, and the rear intercept should be greater than 25 mm.The image-side detector chip was 13.3 mm × 13.3 mm, and the effective width of the object-side grating was 13.69 mm × 13.69 mm.Then,the magnification was[25]

    whereyandy′ represent the object-side and imageside field sizes, respectively.Applying Equation(28), the system magnification was -0.971 5.

    To ensure the telecentricity of the system, software-related operands were used to control the angle of the deviating main ray from the optical axis.The telecentric degree should be as small as possible according to the design requirements of the telecentric system.In the optimization function setting, entrance pupil was set to 0.500, 0.707, and

    1.000 to optimize the telecentric degree of the system.For a dual-telecentric lens optical system, distortion of the optical system distorts the interference fringe and reduces the measurement accuracy of the wind speed.Therefore, in addition to the above parameters, it is necessary to focus on optimizing the system distortion.

    According to the Nyquist frequency formula,

    As the detector pixel sizedwas 13 μm, the cutoff frequency should be greater than 38.46 lp/mm.At 38.46 lp/mm, the MTF of the exit optics was greater than 0.5[26].The MTF describes the relationship between the contrast of the object and the image at different spatial frequencies after the object passes through the optical system.The MTF is related to both the aberration of the optical system and its diffraction effects.This is an objective and reliable image quality evaluation method.The design results for the exit optics are shown in Figure 8 and Table 2.

    Tab.2 Design indices of the dual telecentric optical system

    In the measured target band, the MTF of the entrance and exit optics was close to the diffraction limit, which offers good imaging quality and satisfies all of the indicators.Two optical systems and interference modules were placed into the same system via software, and the image plane of the entrance optics (the end grating surface of the two optical paths) was successfully matched with the object plane of the exit optics (the end grating surface of the two optical paths).Two-dimensional and three-dimensional images of the system are shown in Figure 9 (a) and 9 (b).The surfaces of all lenses of the two optical systems were spherical,making them easy to process.Considering the thermal difference and economic conditions, reasonable materials were selected to produce the real object.

    3 Experimental verification

    3.1 Verification results

    To verify whether the design results meet the requirements of ground-based detection and whether the ground-based DASH interferometer can receive a sufficient etendue and SNR, an experimental platform for the indoor DASH interferometer was constructed based on the design results.A physical diagram of the interferometer is shown in Figure 10,and the basic parameters are shown in Table 3.The specifications of DASH system are shown in Table 4.One optical path of the grating was glued,and the other optical path required experimental debugging.

    Tab.3 Basic parameters of DASH interferometer

    Krypton lamp emission light contains a 557.0 nm emission line, and the light intensity is weak.Therefore, a krypton lamp could be used to simulate the airglow spectrum after passing through the narrowband filter in the laboratory experiment.The narrowband filter is shown in Figure 11.The narrowband filter diameter was 50 mm, the thickness was 5 mm, the material was fused quartz, the central wavelength was 557.7 nm, the bandwidth was 3 nm,and the peak transmittance was no less than 95%.

    Tab.4 DASH system specifications

    The laboratory experimental platform is shown in Figure 12 (a) and (b).The target light emitted by the krypton lamp was transmitted through the optical fiber into the collimation system.The collimated target light was semi-inverse and semi-permeable to the wind speed simulator owing to beam splitter 1.The wind speed simulator was driven by a motor with a turntable attached with reflective film, so that the incident light was reflected into the entrance optics with a narrow-band filter after the frequency shift.The signal was imaged on the two gratings of the interference module by the entrance optics, and the interference fringe was formed by the interference module.The interference fringe was scaled and imaged on the detector using the exit optics.

    The DASH interferometer uses an absolute phase to invert the wind speed differentially.The Littrow angle is key to determining the spatial frequency and phase of the interferogram.Therefore,an accurate Littrow angle assembly is crucial to the accuracy of wind speed measurements.The krypton lamp was used as the light source, and the grating Littrow angle was finely tuned mechanically.The krypton light source, with a wavelength of 557.0 nm, formed interference fringes through the interferometer that were received by the detector.The interference pattern is shown in Figure 13.Because the krypton lamp was relatively bright compared with the gas glow, the integration time in the laboratory experiment was only 1/60 that of the ground experiment, which was too small and lead to a transverse bright line in the detector that was unrelated to the interferogram.It was preliminarily proven that the interference module and optical system were reasonably designed.

    Laboratory debugging experiments have proven the feasibility of using krypton lamps as the light source for DASH interferometer detection.On this basis, the DASH interferometer was glued to an integrated system with a stable structure and mounted on a tripod with an adjustable pitch angle, as shown in Figure 14.Ground-based detection experiments were conducted using night airglow as the detection target.

    Because the brightness of night airglow is very low, the detection distance is large, and stray light may also affect the experimental results in an environment with as high air transmittance as possible.Therefore, a field with almost no stray light was selected to conduct the ground-based detection experiment under sunny and cloudy weather conditions.The outdoor ground-based experiment was carried out in Changchun, China, on May 14, 2022, and a photograph of the experimental setup is shown in Figure 15.On the day of the experiment, the weather was clear with a full moon.In order to avoid the influence of moonlight as background light, a narrowband filter was still used in front of the optical system to filter out most of the background light.The measurement time was selected to be 9 p.m.Beijing time when the moon had not fully risen to its zenith and rested instead at 45° to that.The device faced away from the moon, toward the darkest part of the sky at an elevation of 70 degrees.During the detection process, the integration time was 60 s.At the same time, an EMCCD with refrigeration and gain functions was selected.The cooling temperature was -70 °C and the gain was 20 to ensure that the influence of background light was minimized.The interference fringe is shown in Figure 16.

    On this basis, we also conducted several experiments, among which the indoor and external interference graphs are obtained on September 27, 2022.They are shown in Fig.17 and Fig.18.

    Fig.1 Diagram of the DASH interferometer structure

    Fig.2 Curves of the modulation and phase difference varying with optical path difference of the 557.7 nm spectral line at 90-110 km

    Fig.3 Curves of the modulation and phase difference varying with optical path difference of the 557.7 nm spectral line at 150-300 km

    Fig.4 Plot of the optical path difference versus the interferogram intensity difference

    Fig.5 Optical path diagram of the interference module

    Fig.6 (a) Relationship between signal to noise ratio and pixel position; (b) graph of variation of wind speed error with SNR

    Fig.7 The MTF curves of the entrance optics

    Fig.8 (a) MTF of the exit optics and (b) dot chart of the exit optics

    Fig.9 (a) Two-dimensional system diagram; (b) three-dimensional system diagram

    Fig.10 Real image of the interferometer

    Fig.11 Real image of the narrow band filter

    Fig.12 (a) Schematic diagram of the laboratory debugging device; (b) photograph of the laboratory debugging device

    Fig.13 Streaks made by krypton lamps in the laboratory

    Fig.14 Integrated ground-based DASH interferometer

    Fig.15 Photograph of the ground experiment setup

    Fig.16 Night airglow interferogram

    Fig.17 Interferogram of krypton lamp on September 27

    Fig.18 Interferogram of night airglow on September 27

    Fig.19 Schematic diagram of wind speed simulation

    However, there were angle differences between the stripes measured in the indoor and outdoor experiments, and the reasons are analyzed as follows:

    (1)The bump in the transportation process would cause a slight deviation of Littrow angle.After the temperature control system is improved,the system will be reinforced to avoid this problem.

    (2) There was a gap between the indoor and outdoor temperatures, and the system was very sensitive to temperature.There was no temperature control system at that moment.To improve the temperature control system, this problem could be avoided.

    The obvious interference fringes shown in Figures 16-18 prove that the system design results of the DASH interferometer are reasonable, and the signal-to-noise ratio and luminous flux of the system meets the detection requirements, verifying the feasibility of ground-based detection of the DASH interferometer.

    On this basis, we will invert the interference fringes measured in the laboratory.The speed of the beam carried by the krypton lamp and the wind speed simulator is:

    whereNsis the rotating speed of the turntable of the wind speed simulator,ris the distance from the round point of the turntable of the wind speed simulator,αris the angle between the beam and the reflecting disk, andrris the angle from the vertical direction.Schematic diagram is shown in Figure 19.

    When the rotational speed is 1 810, 2 021 and 2 203 r/min, the simulated wind speed is 37.527,41.9 and 45.675 m/s, respectively.The corresponding interference diagram is shown in Fig.20.

    Fig.20 The corresponding interferogram is 1 810, 2 021 and 2 203 r/min, respectively

    The measured interference fringe pattern shown in Fig.20 is preprocessed, and the average phase difference obtained by inversion is 0.087,0.148 7 and 0.221 8 rad, respectively.The corresponding wind speed is 66.425, 113.1 and 168.7 m/s.

    Due to the lack of calibration, the accurate phase of fundamental frequency cannot be obtained for the time being.During the experiment, the ambient temperature changes greatly and exceed the ambient temperature required by the experiment(20 °C±0.1 °C).Moreover, the wind speed simulator will also produce errors, resulting in the measured wind speed errors.Although the wind speed error is large, the wind speed obtained by inversion is within a reasonable range and the function of the principle verification level of the interferometer can be verified.

    4 Conclusion

    Based on the characteristics of the 557.7 nm oxygen green atom airglow spectral line, a ground DASH interferometer was developed, its design scheme was formulated, and the parameters of its interference module were calculated.Simultaneously, the optical system of the DASH interferometer was designed.The system entrance and exit optics were optimized with the image telecentric mirror group and double telecentric mirror group as initial structures.In the target band, the image telecentric degree of the entrance optics was less than 0.006°,and the distortion was less than 0.05%.The image square and object square telecentricity of the exit optics were both less than 7.0 × 10-3, the distortion was less than 2.4 × 10-5%, the Nyquist space frequency was 38.46 lp/mm, the MTF was greater than 0.5,and all indicators met the requirements.Compared to other passive wind measuring systems, this system has a large FOV and is designed for thermal dissipation.The half-field angle reached 2.815°, the system SNR was approximately 113.75, and the etendue was approximately 0.095 25 cm2sr.Clear interference fringes were detected in both indoor and outdoor experiments, which proves that the groundbased DASH interferometer with thermal compensation and high SNR and large etendue meet the requirements of ground-based detection, and the feasibility of DASH interferometer ground-based detection is verified.

    色尼玛亚洲综合影院| 岛国视频午夜一区免费看| 国产视频一区二区在线看| 成人精品一区二区免费| 两人在一起打扑克的视频| 狠狠狠狠99中文字幕| 一夜夜www| 1024香蕉在线观看| 91字幕亚洲| 国产精品99久久99久久久不卡| 99精品在免费线老司机午夜| 99热这里只有精品一区 | 精品免费久久久久久久清纯| 午夜福利18| 黄色成人免费大全| 久久青草综合色| 无人区码免费观看不卡| 俄罗斯特黄特色一大片| 亚洲av成人av| 巨乳人妻的诱惑在线观看| 少妇粗大呻吟视频| 国产精品98久久久久久宅男小说| 亚洲欧洲精品一区二区精品久久久| 可以在线观看的亚洲视频| 中文字幕人妻丝袜一区二区| 国产又爽黄色视频| 亚洲全国av大片| 欧美国产精品va在线观看不卡| 最新美女视频免费是黄的| 欧美丝袜亚洲另类 | 亚洲aⅴ乱码一区二区在线播放 | 亚洲自拍偷在线| 国产伦一二天堂av在线观看| 一本精品99久久精品77| 国产v大片淫在线免费观看| 中亚洲国语对白在线视频| 免费女性裸体啪啪无遮挡网站| 国产av不卡久久| 午夜福利一区二区在线看| 精品电影一区二区在线| 国产亚洲精品av在线| 国产又黄又爽又无遮挡在线| 欧美中文日本在线观看视频| 国产国语露脸激情在线看| 午夜免费激情av| 听说在线观看完整版免费高清| 国产精品二区激情视频| 久久久久久九九精品二区国产 | 日日摸夜夜添夜夜添小说| 黑丝袜美女国产一区| 男人舔女人的私密视频| 久久国产精品影院| 日本一本二区三区精品| 亚洲人成伊人成综合网2020| 人妻丰满熟妇av一区二区三区| tocl精华| 最近最新中文字幕大全免费视频| 岛国视频午夜一区免费看| 亚洲自拍偷在线| 国产精品久久久av美女十八| 久久久久久人人人人人| 午夜两性在线视频| 亚洲va日本ⅴa欧美va伊人久久| 精品一区二区三区av网在线观看| aaaaa片日本免费| 久久中文字幕一级| 亚洲熟妇熟女久久| 宅男免费午夜| 美女大奶头视频| 国产成人系列免费观看| 欧美又色又爽又黄视频| 妹子高潮喷水视频| 国产精品免费视频内射| 91在线观看av| netflix在线观看网站| 国产精品久久电影中文字幕| 免费看日本二区| 国产视频内射| 免费在线观看黄色视频的| 久久久久九九精品影院| 搡老妇女老女人老熟妇| 亚洲av电影在线进入| 国产黄a三级三级三级人| 99在线人妻在线中文字幕| 久久久久久久久久黄片| 哪里可以看免费的av片| 99精品欧美一区二区三区四区| 黄片小视频在线播放| 国产免费av片在线观看野外av| 久久婷婷人人爽人人干人人爱| 人人妻人人澡欧美一区二区| 神马国产精品三级电影在线观看 | xxx96com| 窝窝影院91人妻| 午夜免费鲁丝| 色综合婷婷激情| 亚洲 欧美一区二区三区| 日韩欧美国产在线观看| 国产片内射在线| 精品日产1卡2卡| 99热6这里只有精品| 久久久久九九精品影院| 免费看十八禁软件| 一级片免费观看大全| 男女视频在线观看网站免费 | 国产私拍福利视频在线观看| 国产亚洲av高清不卡| 国产精品自产拍在线观看55亚洲| 国产伦人伦偷精品视频| 久久午夜亚洲精品久久| 特大巨黑吊av在线直播 | 黑人操中国人逼视频| 久久久久久久久免费视频了| 国内久久婷婷六月综合欲色啪| 欧美成人免费av一区二区三区| 在线观看免费日韩欧美大片| 久久天堂一区二区三区四区| 男人舔女人的私密视频| 国产精品久久久久久亚洲av鲁大| 日韩免费av在线播放| 午夜福利高清视频| 嫩草影视91久久| www.自偷自拍.com| 美女大奶头视频| 男女下面进入的视频免费午夜 | 日韩大码丰满熟妇| 国产极品粉嫩免费观看在线| 最近在线观看免费完整版| 久久久国产成人免费| 99久久精品国产亚洲精品| 国产精品 欧美亚洲| 国产av一区在线观看免费| 怎么达到女性高潮| 黄色视频,在线免费观看| 精华霜和精华液先用哪个| 午夜福利欧美成人| 国产av一区二区精品久久| 免费在线观看完整版高清| 国产亚洲欧美在线一区二区| 日本免费a在线| 久久中文字幕一级| 在线观看午夜福利视频| 青草久久国产| 亚洲欧美激情综合另类| 国产97色在线日韩免费| 法律面前人人平等表现在哪些方面| 一本精品99久久精品77| 99国产精品99久久久久| 激情在线观看视频在线高清| 侵犯人妻中文字幕一二三四区| 国产99白浆流出| 国产三级在线视频| 波多野结衣av一区二区av| 在线观看一区二区三区| 搡老岳熟女国产| 99re在线观看精品视频| 国产一区在线观看成人免费| 国产av一区二区精品久久| 亚洲av美国av| 黄色视频不卡| 搡老妇女老女人老熟妇| 久久九九热精品免费| 国产色视频综合| 一区二区三区精品91| 97超级碰碰碰精品色视频在线观看| 在线看三级毛片| 日韩欧美在线二视频| 国内精品久久久久精免费| 亚洲一区二区三区不卡视频| 亚洲av五月六月丁香网| 12—13女人毛片做爰片一| 99国产极品粉嫩在线观看| 午夜福利欧美成人| 久久精品影院6| 制服丝袜大香蕉在线| 国产成人影院久久av| 在线观看免费日韩欧美大片| e午夜精品久久久久久久| 精品午夜福利视频在线观看一区| 日韩大尺度精品在线看网址| 欧美国产精品va在线观看不卡| 精品福利观看| 亚洲精品在线美女| 男人舔奶头视频| 最好的美女福利视频网| 一进一出好大好爽视频| 国产精品自产拍在线观看55亚洲| 色精品久久人妻99蜜桃| 别揉我奶头~嗯~啊~动态视频| √禁漫天堂资源中文www| 不卡av一区二区三区| 搡老岳熟女国产| 国产亚洲精品久久久久5区| 99热这里只有精品一区 | 欧美性猛交黑人性爽| 老司机午夜十八禁免费视频| 亚洲中文日韩欧美视频| 在线观看66精品国产| 日韩欧美 国产精品| 一级a爱片免费观看的视频| 欧美黄色片欧美黄色片| 俄罗斯特黄特色一大片| 午夜福利视频1000在线观看| av福利片在线| 日本在线视频免费播放| 波多野结衣巨乳人妻| 黄网站色视频无遮挡免费观看| 久久国产精品男人的天堂亚洲| 长腿黑丝高跟| 亚洲五月天丁香| 精品一区二区三区四区五区乱码| 一本精品99久久精品77| 午夜日韩欧美国产| 少妇 在线观看| 午夜视频精品福利| 精品国产美女av久久久久小说| 女性生殖器流出的白浆| 大型av网站在线播放| 国产成人一区二区三区免费视频网站| 级片在线观看| 啦啦啦韩国在线观看视频| 欧美精品啪啪一区二区三区| 黄色 视频免费看| 女人高潮潮喷娇喘18禁视频| 国产av一区二区精品久久| 欧美日韩亚洲综合一区二区三区_| 18禁裸乳无遮挡免费网站照片 | 91老司机精品| 国产乱人伦免费视频| 免费看a级黄色片| 欧美午夜高清在线| 日本a在线网址| 97超级碰碰碰精品色视频在线观看| 日韩大尺度精品在线看网址| 精品福利观看| 自线自在国产av| 日本精品一区二区三区蜜桃| av超薄肉色丝袜交足视频| 精品一区二区三区四区五区乱码| 久久热在线av| 黑丝袜美女国产一区| 国产一区二区三区视频了| 精品人妻1区二区| 亚洲精品色激情综合| 久久久精品欧美日韩精品| 国产欧美日韩一区二区精品| av片东京热男人的天堂| 91老司机精品| 亚洲专区中文字幕在线| 欧美zozozo另类| 亚洲一卡2卡3卡4卡5卡精品中文| 国产片内射在线| 国产乱人伦免费视频| 国产亚洲精品第一综合不卡| 深夜精品福利| tocl精华| 啦啦啦韩国在线观看视频| 午夜两性在线视频| 亚洲国产精品999在线| 精品久久久久久久毛片微露脸| 麻豆成人午夜福利视频| 日本三级黄在线观看| 美国免费a级毛片| 日本五十路高清| 亚洲中文av在线| 热99re8久久精品国产| 国产精品电影一区二区三区| 午夜视频精品福利| 亚洲av成人av| 国产亚洲av嫩草精品影院| 丰满人妻熟妇乱又伦精品不卡| 手机成人av网站| 成人三级黄色视频| 亚洲狠狠婷婷综合久久图片| 国产国语露脸激情在线看| 看片在线看免费视频| 12—13女人毛片做爰片一| 日韩欧美一区二区三区在线观看| 91av网站免费观看| 很黄的视频免费| 热re99久久国产66热| 丝袜在线中文字幕| 久久中文看片网| 久久久久久久精品吃奶| 亚洲成人久久爱视频| 50天的宝宝边吃奶边哭怎么回事| 久久国产精品影院| 欧美午夜高清在线| 久久精品成人免费网站| 国产亚洲欧美精品永久| 757午夜福利合集在线观看| 无遮挡黄片免费观看| 天堂动漫精品| 禁无遮挡网站| 久久热在线av| 亚洲五月色婷婷综合| 18禁国产床啪视频网站| 欧美成人一区二区免费高清观看 | 亚洲 欧美一区二区三区| 熟妇人妻久久中文字幕3abv| 男人舔奶头视频| 色综合站精品国产| 午夜福利在线观看吧| 熟女少妇亚洲综合色aaa.| 久久精品国产综合久久久| 午夜免费成人在线视频| 丝袜美腿诱惑在线| 成年女人毛片免费观看观看9| 欧美乱码精品一区二区三区| 亚洲av美国av| 韩国av一区二区三区四区| 国产精品98久久久久久宅男小说| 国产真人三级小视频在线观看| 两个人视频免费观看高清| 久久草成人影院| 精品卡一卡二卡四卡免费| 欧美国产精品va在线观看不卡| 亚洲真实伦在线观看| 亚洲,欧美精品.| 俺也久久电影网| 亚洲人成网站高清观看| 女同久久另类99精品国产91| 国产亚洲精品第一综合不卡| 免费高清在线观看日韩| 99热6这里只有精品| 黄色毛片三级朝国网站| 欧美又色又爽又黄视频| 国产av不卡久久| 丝袜人妻中文字幕| 精品日产1卡2卡| 成年版毛片免费区| av欧美777| 视频在线观看一区二区三区| 欧美另类亚洲清纯唯美| 亚洲av五月六月丁香网| 国产私拍福利视频在线观看| 两性夫妻黄色片| 啦啦啦韩国在线观看视频| 波多野结衣高清无吗| 国产亚洲精品久久久久5区| 老司机福利观看| 白带黄色成豆腐渣| 亚洲第一青青草原| 日韩高清综合在线| 精品久久久久久久久久久久久 | 波多野结衣av一区二区av| 亚洲七黄色美女视频| 国产精品 欧美亚洲| 亚洲成人免费电影在线观看| 国产av不卡久久| 18禁黄网站禁片午夜丰满| 欧美性长视频在线观看| 国产私拍福利视频在线观看| 亚洲色图 男人天堂 中文字幕| 久久久久久久久中文| 中文字幕人成人乱码亚洲影| 亚洲精华国产精华精| 九色国产91popny在线| 丝袜在线中文字幕| 久久久久久国产a免费观看| 亚洲人成电影免费在线| 夜夜躁狠狠躁天天躁| 亚洲成人久久性| 成人亚洲精品一区在线观看| 亚洲欧美精品综合久久99| 国产精品1区2区在线观看.| 色综合亚洲欧美另类图片| 在线观看午夜福利视频| 免费观看精品视频网站| 国产爱豆传媒在线观看 | 成人亚洲精品一区在线观看| 久久香蕉激情| 日日摸夜夜添夜夜添小说| 一级a爱片免费观看的视频| 狠狠狠狠99中文字幕| av片东京热男人的天堂| 欧美绝顶高潮抽搐喷水| 免费看a级黄色片| 免费女性裸体啪啪无遮挡网站| xxx96com| 18禁观看日本| 国产精品精品国产色婷婷| 首页视频小说图片口味搜索| 黑人巨大精品欧美一区二区mp4| 国产黄a三级三级三级人| 国产97色在线日韩免费| 欧美日韩精品网址| 在线永久观看黄色视频| 亚洲一卡2卡3卡4卡5卡精品中文| 男女那种视频在线观看| 国产午夜福利久久久久久| 色综合站精品国产| 午夜福利成人在线免费观看| 婷婷精品国产亚洲av| 免费电影在线观看免费观看| 在线观看日韩欧美| 可以在线观看毛片的网站| 97碰自拍视频| 欧美成人性av电影在线观看| 九色国产91popny在线| 男人舔女人的私密视频| 制服丝袜大香蕉在线| 久久午夜综合久久蜜桃| 亚洲国产欧美日韩在线播放| 天堂√8在线中文| 国产日本99.免费观看| 白带黄色成豆腐渣| 国产精品免费一区二区三区在线| 三级毛片av免费| 色婷婷久久久亚洲欧美| 国产91精品成人一区二区三区| 三级毛片av免费| 曰老女人黄片| 十八禁网站免费在线| 很黄的视频免费| 国产精品一区二区精品视频观看| 88av欧美| 精品欧美国产一区二区三| 成人永久免费在线观看视频| 午夜a级毛片| 女警被强在线播放| 激情在线观看视频在线高清| 亚洲久久久国产精品| 成年免费大片在线观看| 国产一区二区三区视频了| 国产亚洲欧美在线一区二区| 成人亚洲精品av一区二区| 老司机靠b影院| 国产精品精品国产色婷婷| 级片在线观看| 999久久久国产精品视频| 69av精品久久久久久| 一级毛片女人18水好多| 久久久久亚洲av毛片大全| a在线观看视频网站| 麻豆一二三区av精品| 亚洲国产看品久久| 国产亚洲av嫩草精品影院| 国产成人精品久久二区二区免费| 精品久久久久久,| 欧美人与性动交α欧美精品济南到| 岛国在线观看网站| 成人亚洲精品一区在线观看| 欧美国产日韩亚洲一区| 日韩欧美国产一区二区入口| 一个人观看的视频www高清免费观看 | 桃色一区二区三区在线观看| 亚洲精品一卡2卡三卡4卡5卡| 国产午夜精品久久久久久| 色综合欧美亚洲国产小说| 丝袜在线中文字幕| 亚洲黑人精品在线| 热re99久久国产66热| 天堂动漫精品| 国产成人啪精品午夜网站| av片东京热男人的天堂| 黑人巨大精品欧美一区二区mp4| 欧美成人性av电影在线观看| 国产成人精品久久二区二区免费| 丰满的人妻完整版| 女人被狂操c到高潮| 嫩草影视91久久| 国产精华一区二区三区| 婷婷精品国产亚洲av| 51午夜福利影视在线观看| 国产亚洲精品综合一区在线观看 | 亚洲人成77777在线视频| 热re99久久国产66热| 狠狠狠狠99中文字幕| 熟女电影av网| 中国美女看黄片| 亚洲色图av天堂| 亚洲精品美女久久av网站| 俄罗斯特黄特色一大片| 久久久久久久久免费视频了| 中文字幕人成人乱码亚洲影| 18禁黄网站禁片午夜丰满| 精品不卡国产一区二区三区| 满18在线观看网站| 99久久99久久久精品蜜桃| 亚洲精品美女久久av网站| 制服人妻中文乱码| 男人舔女人下体高潮全视频| 一本一本综合久久| 夜夜看夜夜爽夜夜摸| 久久精品91无色码中文字幕| 欧美 亚洲 国产 日韩一| 国产精品爽爽va在线观看网站 | 国产精品野战在线观看| АⅤ资源中文在线天堂| 日本一本二区三区精品| 一进一出好大好爽视频| 天堂影院成人在线观看| 久热这里只有精品99| 久久久久久亚洲精品国产蜜桃av| 在线观看www视频免费| cao死你这个sao货| av电影中文网址| 国产精品日韩av在线免费观看| 国产人伦9x9x在线观看| 美女午夜性视频免费| 亚洲免费av在线视频| 999久久久国产精品视频| 91国产中文字幕| 黄色成人免费大全| 精品免费久久久久久久清纯| 色播在线永久视频| 啦啦啦 在线观看视频| 日韩欧美国产在线观看| 丁香欧美五月| 亚洲一区高清亚洲精品| 午夜福利在线在线| 黄色片一级片一级黄色片| 午夜福利高清视频| 久99久视频精品免费| 日韩免费av在线播放| 欧美日韩一级在线毛片| 99久久久亚洲精品蜜臀av| 欧美日韩一级在线毛片| 一进一出抽搐gif免费好疼| av福利片在线| 亚洲专区中文字幕在线| 热99re8久久精品国产| 久久欧美精品欧美久久欧美| 热99re8久久精品国产| 99国产极品粉嫩在线观看| 97超级碰碰碰精品色视频在线观看| 国产三级在线视频| www国产在线视频色| 在线观看一区二区三区| 欧洲精品卡2卡3卡4卡5卡区| 日韩欧美一区二区三区在线观看| 级片在线观看| 18禁国产床啪视频网站| 怎么达到女性高潮| 精品免费久久久久久久清纯| 国产成人欧美| 99热6这里只有精品| 色综合亚洲欧美另类图片| 亚洲国产欧美网| 麻豆av在线久日| 一进一出抽搐gif免费好疼| 国产亚洲欧美精品永久| 黄色视频不卡| 亚洲色图av天堂| 三级毛片av免费| 日韩欧美一区二区三区在线观看| 久久精品国产综合久久久| 日韩高清综合在线| 亚洲真实伦在线观看| 成人18禁高潮啪啪吃奶动态图| 一二三四社区在线视频社区8| 一本综合久久免费| 亚洲色图 男人天堂 中文字幕| www日本黄色视频网| 窝窝影院91人妻| 老鸭窝网址在线观看| 日韩欧美 国产精品| www.精华液| 亚洲无线在线观看| 又紧又爽又黄一区二区| 免费在线观看亚洲国产| 免费av毛片视频| 91九色精品人成在线观看| 午夜福利视频1000在线观看| 97人妻精品一区二区三区麻豆 | 亚洲自拍偷在线| 我的亚洲天堂| 人成视频在线观看免费观看| 自线自在国产av| 久久午夜综合久久蜜桃| 搡老妇女老女人老熟妇| 亚洲精品国产区一区二| 成熟少妇高潮喷水视频| 久久久久久久午夜电影| 在线观看免费视频日本深夜| 十分钟在线观看高清视频www| 亚洲国产欧洲综合997久久, | 后天国语完整版免费观看| 精品国产一区二区三区四区第35| 国产激情偷乱视频一区二区| 最近最新免费中文字幕在线| 婷婷精品国产亚洲av在线| 亚洲人成网站高清观看| 少妇裸体淫交视频免费看高清 | 无遮挡黄片免费观看| 淫妇啪啪啪对白视频| 欧美zozozo另类| 久久欧美精品欧美久久欧美| 欧美久久黑人一区二区| 欧美日韩乱码在线| 老熟妇仑乱视频hdxx| 在线播放国产精品三级| 午夜亚洲福利在线播放| 国产爱豆传媒在线观看 | 欧美成人午夜精品| 亚洲成国产人片在线观看| 国产一区二区在线av高清观看| 国产精品一区二区精品视频观看| 成年免费大片在线观看| 亚洲精品美女久久av网站| 亚洲精品一卡2卡三卡4卡5卡| 夜夜夜夜夜久久久久| 国产欧美日韩精品亚洲av| 亚洲专区国产一区二区| 成年女人毛片免费观看观看9| 亚洲精品国产精品久久久不卡| 在线观看免费视频日本深夜| 99久久综合精品五月天人人| 国产精品日韩av在线免费观看| 亚洲熟女毛片儿| 欧美另类亚洲清纯唯美| 日韩视频一区二区在线观看| 麻豆久久精品国产亚洲av| 成人18禁在线播放|