999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

The development, performances and applications of the monochromatic X-rays facilities in(0.218-301)keV at NIM,China

2021-07-02 09:31:58SiMingGuoJinJieWuDongJieHou
Nuclear Science and Techniques 2021年6期

Si-Ming Guo· Jin-Jie Wu · Dong-Jie Hou

Abstract Space scientific exploration is rapidly becoming the primary battlefield for humankind to explore the universe.Countries worldwide have launched numerous space exploration satellites. Accurate calibration of the detectors on the ground is a crucial element for space science satellites to obtain observational results.For the purpose of providing calibration for various satellite-borne detectors,multiple monochromatic X-rays facilities have been built at the National Institute of Metrology, P.R. China (NIM).These facilities mainly pertain to grating diffraction and Bragg diffraction, and the energy range of the produced monochromatic X-rays is 0.218-301 keV. These facilities have a high performance in terms of energy stability,monochromaticity, and flux stability. The monochromaticity was greater than 3.0%.The energy stability of the facility is 0.02% at 25 keV over 8 h, and the flux stability was within 1.0% at 25 keV over 8 h. Calibration experiments on the properties of satellite-borne detectors,such as energy linearity, energy resolution, detection efficiency,and temperature response, can be conducted at the facilities. Thus far, the calibration of two satellites has been completed by the authors, and the work on three other satellites is in progress. This study will contribute to the advancement of X-ray astronomy the development of Chinese space science.

Keywords Monochromatic X-rays · Monochromator ·Bragg diffraction · Energy spectra · Calibration

1 Introduction

The use of rockets or high-altitude balloons and other vehicles began in the 1960s; it was also when collimating scintillation crystals or gas detectors were employed to realize space X-ray detection. This marked the start of space astronomical observation research. X-ray astronomy largely involves the study of high-energy radiation celestial bodies through X-ray observation. The main observation objects include black holes, neutron stars, as well as interstellar high-temperature hot gases.Studies in this field largely focus on physical processes under extreme conditions such as extremely high density, extremely strong magnetic fields, and extremely strong gravitational fields.Between 2000 and 2014, the USA launched a total of 84 space exploration satellites, accounting for 27.6% of the total number of satellites launched worldwide.In 2017,the first Chinese space science satellite Hard X-ray Modulation Telescope (HXMT) was successfully launched, relaying significant research data [1, 2]. Presently, China has numerous space science satellites under development or are ready to be launched; examples include the gravitational wave high-energy electromagnetic counterpart all-sky monitor(GECAM)[3],hard X-ray imager(HXI)[4],space variable objects monitor (SVOM) [5, 6], and the enhanced X-ray timing and polarimetry mission (eXTP) [7]. The detectors utilized by many space science satellites are integral to space observations. Key parameters such as energy resolution, detection efficiency, and uniformity of these detectors must be ensured through ground calibration.Otherwise, the inversion from observation data to real astronomical objects cannot be realized. Hence, research on monochromatic X-ray calibration facilities and calibration methods for detectors is required.There are several methods to produce monochromatic X-rays, namely radionuclides, synchrotron radiation, K fluorescence, and Bragg diffraction. Although radioactive sources are broadly used in detector calibration,their energy values are constant, rare, and uncontrollable. The synchrotron X-ray source is extensively used in a variety of fields, but it has the disadvantages of high cost, inconvenient application,and instability. K fluorescent X-rays are generated by primary X-rays to excite the secondary target. Owing to material restrictions, they are limited to a specific energy.Based on the continuous spectrum from an X-ray machine,monochromatic X-rays can be obtained through Bragg diffraction,whose energy is continuous and can be adjusted using different crystals and by altering the Bragg angle.This method is convenient, effective, and low cost for the application of detector performance calibration.

Internationally, monochromatic X-ray radiation facilities predicated on X-ray machines include the SOLEX device from the Becquerel laboratory in France [9-11],XACT in Italy [12],XRCF in the USA [13]and PANTER in Germany [14]; the main performance comparison between them is shown in Table 1. Satellites calibrated by these facilities include BeppoSAX[15],INTEGRAL/JEMX [16], AXAF [17], Swift [18], and GBM [19, 20]. With regard to China, the HXCF [21] was jointly built by the National Institute of Metrology (NIM) of People’s Republic of China and the Institute of High Energy Physics, Chinese Academy of Sciences.The energy range was 15-168 keV. The said facilities provide on-ground calibrations on energy linearity, energy resolution, detection efficiency, energy response matrix, and effective area of China’s first astronomical satellite,HXMT[22],which laid the foundation for future scientific achievements [1, 23].

Detectors carried by space satellites form the core of space observation. Parameters such as the energy resolution,detection efficiency,and uniformity of these detectors must be calibrated on the ground. Otherwise, it is impossible to realize inversion from observation data to real astronomical phenomena. Therefore, it is necessary to develop monochromatic X-ray calibration devices and detector calibration methods. At present, research on monochromatic X-ray facilities in hard X-ray sections has achieved considerable research results and applications. In view of the shortcomings of monochromatic X-ray measurements and detector calibration in the soft X-ray energy section,combined with the demand for future space science development in China, we study monochromatic X-ray facilities in (0.218-301) keV to meet the demands for calibration of detectors. Owing to the large energy range,four sets of facilities were developed to realize monochromatic X-rays, and measurement research was also carried out.The theoretical calculation of the detection efficiency of the standard detector and experimental calibration under a standard radioactive source were completed.The measurement results were traced to the national primary standard for radionuclide activity. A traceability system of monochromatic X-ray flux in the energy range was established to ensure the accuracy and reliability of the applications. The research results will provide calibration services for space science satellite projects that have been established in China and data traceability support for nextgeneration space science satellites in space science planning in China. In addition, it could provide calibration support for the detectors used in dark matter research andgravitational wave detection research and may promote China’s breakthroughs in the frontier scientific area. This research achievement may support the development of space astronomy in China and enhance Chinese international competitiveness in X-ray astronomy.

Table 1 The emblematic monochromatic X-rays facilities in the world

2 Methods and experiment setup

For building monochromatic X-ray facilities to provide calibration and research conditions for space science satellite-borne detectors and other detectors, much research has been conducted at NIM, and several facilities have been come up. To date, (0.218-301) keV monochromatic X-rays generated by the facilities have been completed. The (0.218-1.6) keV monochromatic X-rays are generated by X-ray machines and diffracted by gratings, whereas a (5-301) keV energy range is produced using the Bragg diffraction of crystals.

2.1 The (0.218-1.6) keV X-ray facility

The (0.218-1.6) keV monochromatic X-ray beam facility predominantly consists of an X-ray machine,a nonharmonic grating monochromator, and a detector. The Bremsstrahlung radiation spectrum generated by the X-ray machine enters the non-harmonic monochromator to achieve dispersion focusing, upon which the adjustable monochromatic X-ray beam is obtained at the exit. The source is a customized multi-target windowless X-ray machine (WORX, model XWT-065-SE) with five target materials. The structure is illustrated in Fig. 1.

Fig. 1 (Color online) Schematic configuration of the soft X-ray monochromator

Considering the convenience and stability provided by the focused light source, the monochromator is directly connected to the source without a front-focusing system.The non-harmonic grating monochromator adopts the CT optical structure and employs two toroidal mirrors to realize collimation and focusing functions. The focus positions are at the entrance and exit of the slit. The incident and diffracted beams of the grating are parallel in the spectral resolution direction. During wavelength scanning,only one element of the grating needs to be rotated,and the focusing conditions can be met without moving the exit slit.

Characteristic X-rays corresponding to different target materials can be generated when the target materials are changed. The main parameters of the X-ray machine and monochromator are listed in Table 2.

The monochromator could produce (218-1600) eV monochromatic X-rays; if the Bragg diffraction crystal is used to generate the monochromatic X-rays, the energy range may be larger. Because the X-tube is multi-target, it could produce stable characteristic peaks of the targets,such as Si, Kα1: 1.74 keV, Kα2: 1.739 keV, Kβ1:1.832 keV,Ti, Kα:4.5 keV, Kβ:4.93 keV,Lα:0.452 keV,Lβ: 0.458 keV, Ag, Lα1: 2.984 keV, Lα2: 2.978 keV, Lβ1:3.151 keV, Lβ2: 3.384 keV, Lγ: 3.519 keV.

2.2 The (5-40) keV X-ray facility

The (5-40) keV monochromatic X-ray beam facility is principally composed of an X-ray machine, a diffraction crystal, and a synchronous rotating device. The monochromator comprises a high-precision rotator and a Bragg diffraction crystal, with monochromatic X-rays producing and passing through the monochromator. The process follows Bragg’s law [24]:

where d is the crystal constant,n is the diffraction series,θ is the Bragg’s angle,and λ is the wavelength.According to Eq. (2), the energy of monochromatic X-rays can be deduced as follows:

The uncertainty of the photon flux introduced by the position uncertainty is substantially large, and a position deviation of 1 mm may cause a 20% flux error; furthermore, the position accuracy is critical for the detection efficiency calibration of the detector. The double-crystal monochromator is advantageous in that there is only a slight translation in the X-ray emission direction when different energies are adjusted, and no angular deflection occurs. This is beneficial to the detection efficiency calibration. In contrast, the single-crystal monochromator has an energy corresponding to an angle.Every time the energy is changed, the angle needs to be adjusted; thus, every angle change alters the position of the detector,introducinga considerable amount of uncertainty. The best method to address this issue involves fixing the direction of the X-ray emissions,with no need to alter the position of the detector when adjusting different energies. To guarantee that the direction of the X-ray emissions remains unchanged, the design scheme adopted is as shown in Fig. 2.

Table 2 The first instrument characteristics

The monochromator was separately placed on a turntable;the monochromator and the tube were also placed on the same platform, and the upper and lower rotations were controlled by two motors. The upper motor controls the rotation of the crystal and the lower motor controls the X-ray machine. When the crystal rotates at an angle θ, the X-ray machine rotates at an angle of 2θ to ensure that the direction of the monochromatic X-rays does not change.The detailed parameters are listed in Table 3

2.3 The (20-161) keV X-ray facility

The (20-161) keV X-ray facility comprises an X-ray tube, monochromator, detector, and control system, as depicted in Fig. 3.The X-ray continuous spectrum emitted by the X-ray tube subsequently impacts the double-crystal monochromator. The crystal alters the angle under the control of a high-precision turntable and T structure, and monochromatic X-rays are produced. Monochromatic X-rays are ultimately detected after passing through the collimator and the beam-limiting aperture.

Fig. 2 (Color online) Schematic configuration of the (5-40) keV monochromatic X-rays beam facility

The double-crystal monochromator is an X-ray diffraction device with high precision and resolution,utilizing the Bragg diffraction law for crystals, and includes a doublecrystal, high-precision rotator, T structure,beam regulator,and aperture. The principle of the structure is shown in Fig. 3. The monochromator employs two parallel crystals as the original dispersion; the first crystal (‘‘crystal I’’)realizes monochromaticity. Under the action of the fixed height difference structure,the second crystal(‘‘crystal II’’)maintains the exit direction and height of the output monochromatic light relative to the incident light, so as to obtain a light spot with a fixed position. Detailed parameters are listed in Table 4.

2.4 The (21-301) keV X-ray facility

The configuration of the (21-301) keV X-ray facility is illustrated in Fig. 4.Owing to the limitation of the structure of the double-crystal monochromator, higher-energy monochromatic X-rays cannot be generated by the diffraction of a double-crystal monochromator. Singlecrystal monochromators have, therefore, become a better choice.The Bragg angle can be adjusted by controlling the rotator, and thus, different monochromatic X-rays can be generated.

The X-ray machine utilizes the Y.MG605 type X-ray machine of YXLON Company, with a rated voltage range of 20-600 kV. The detailed parameters are listed in Table 5.

3 Performance test and results

3.1 Detectors

3.1.1 HPGe

A high-purity germanium (HPGe) detector (GL0110)was chosen as the first standard. Combined with theparameters provided by the manufacturer, the internal structure of the detector is obtained by industrial CT scanning, the geometric model is established to simulate the detection efficiency by Monte Carlo (MC) codes, andthe point source extrapolation experiment was verified before the detection efficiency of the HPGe was finally obtained. A perspective image of the HPGe is shown in Fig. 5.

Table 3 The second instrument characteristics

Fig. 3 (Color online)Schematic configuration of the(20-161) keV monochromatic X-rays beam facility

Table 4 The third instrument characteristics

Fig. 4 (Color online)Schematic configuration of the high-energy facility

Table 5 The fourth instrument characteristics

An experimental calibration was conducted,as shown in Fig. 6. Three radioactive sources,241Am,57Co, and109Cd,were selected for the experiment. Upon calculation of the experimental data,the experimental efficiencies at different energies were obtained.Specific experimental methods are provided in Ref. [25].

On comparing the simulation and experiments [26], the results were determined to fit well after necessary corrections were met. The results are presented in Table 6.

The HPGe detection efficiency curves are shown in Fig. 7.The test data points in the figure were measured by experiments pursuant to the point source extrapolation method. The curve is calculated by establishing a Monte Carlo model based on the CT scan image of the HPGe detector and is used for the beam flux calculation.

3.1.2 SDD

A silicon drift detector (SDD) is another standard detector. Owing to its high performance, the said detector has been extensively used in scientific research. The detection area is 20 mm2, its length is 450 μm, the thickness of incident beryllium window is 8 μm, and the FWHM is greater than 133 eV at 5.9 keV.

Fig. 5 The CT image of the detector

Fig. 6 (Color online) Calibration experiment

The detection efficiency was determined by calculating the response of the detector to different incident energy photons using MC simulation software, and then the reliability of the calculation results is verified using the known activity of the radioactive source for calibration experiments. The structure data of the detector used in the calculation mainly refer to the factory instructions and industrial CT scanning results. The structure used in the simulation is shown in Fig. 8.

The simulation results show that the detection efficiency of the SDD is excellent below 10 keV, and the detection efficiency decreases rapidly with an increase in energy(Fig. 9).

3.1.3 CCD

Two CCD cameras were used.One was PIXIS-XF 2048 B. It has excellent performance in radiation detection owing to the 2048 × 2048 array and the rapid measurement system. The other CCD was from Andor DO934PBN; the imaging array used was 1024 × 1024, the pixels were 13 μm × 13 μm, and the maximum speed was 5 MHz. The detection efficiency of the two cameras is shown in Fig. 10.

Table 6 Results from simulation and experiments

Fig. 7 The detection efficiency obtained from MC for parallel incident photons and experiment

Fig. 8 (Color online) The SDD structure

Fig. 9 Calculated intrinsic efficiency of SDD

3.2 Data and discussion

3.2.1 Measurements of (0.218-1.6) keV photons

The zero-order spot was recorded on the CCD detector when the grating was adjusted to the horizontal position.The position coordinates of the center point of the zeroorder spot were(1318,1092),and the incident angle was α.The first-order spot was recorded while the grating angle was rotated to 949.8 eV (λ = 1.305 nm). The position of the light spot center could be ultimately located in the zeroorder coordinates (1318, 1092) by further rotating the grating angle slightly; at this point, Δα = 0.851°. Pursuant to grating diffraction equation

Thus, from Δα and λ, the actual incident angle of the zero-order α = 87.0715° can be obtained. Subsequently, a first-order spot of 828 eV was collected on the CCD while the grating was rotated to 0.976°, as shown in Fig. 11.

After replacing the target with Ti and rotating the grating to 1.762°, the first-level spot of 458.3 eV was recorded by the CCD. When the grating was rotated, the limit of the rotation angle was 0.5052° and 3.7° for the 1500 lp/mm grating, and the corresponding energies were 1600 eV and 218 eV, respectively. Hence, the working energy range of the monochromator was 218-1600 eV.On the CCD detection surface, the spot length was 103 pixels(approximately 1.39 mm), and the spot height was 39 pixels (approximately 0.53 mm). The total photon count was 1.64 × 106/min, and hence, the count rate was 2.7 × 104cps at 949.8 eV.

Fig. 10 Detection efficiency of PXF:2048B (left) and DO934P-BN (right)

Fig.11 (Color online)Cu target X-ray source 949.8 eV primary spot(left) and 828 eV primary spot (right)

At the best resolution position of the first-level spot at 949.8 eV (λ = 1.305 nm), the grating incident angle was 87.0661°.The number of pixels for broadening the spectral lines, ΔN = 31(FWHM), was acquired from the measured spot. Furthermore, the distance between the center of the CCD array and the focusing lens was L = 800 mm; hence,the angle broadening corresponding to the FWHM of the spectral line was calculated as ΔN × 13.5 μm/L = 5.23 × 10-4. In accordance with the grating equation:

Next, the energy broadening corresponding to the angular broadening ΔE = 20.496 eV was obtained. Additionally, the energy resolution was obtained as E/ΔE≈46 at 949.8 eV (2.2% at 949.8 eV), as shown in Fig. 12. The energy stability was better than 0.2% at 1 h, and the flux stability was within 0.5% at 1 h. Benefitting from the high precision of the vacuum turntable structure and a very stable X-ray generator, its stability was found to be excellent.

Fig. 12 Experimental result of the 949.8 eV X-ray

3.2.2 Measurements of (5-40) keV X-rays

The SDD detector was used to collect data every 30 min, and the energy stability was obtained through long-term measurement, thereby ensuring the stability of the device, as shown in Fig. 13.

The energy stability of the facility was 0.02%at 25 keV over 8 h, and the flux stability thereof was within 1.0% at 25 keV 8 h.

The energy range of the facility was determined experimentally, whereas the energy spectra were obtained using the SDD and HPGe detectors. A portion of the experimental results is shown in Fig. 14.

Photons are scattered after passing through a 3 mm beam limiting diaphragm, resulting in a spot size larger than the aperture of the diaphragm.The CCD detector was moved to the center of the laser position. The coordinates of the edge of the spot were measured as(739,713),(1016,713), (874, 551), and (874, 881). The image obtained by CDD measurement was 13.5 μm per pixel; thus, thehorizontal length of the spot was 3739.5 μm and the vertical length was 4455 μm, as shown in Fig. 15.

Fig. 13 (Color online) The flux stability at 25 keV

Fig. 14 (Color online) The measured monochromatic spectra under LiF220 crystal monochromator

The HPGe detector was used to measure the performance of the device in more detail. The performance parameters, such as flux and monochromaticity, are listed in Table 7.In the present experimental mode,the minimum flux was greater than 5000 cm-2s-1. This value can be altered by adjusting the parameters of the X-ray machine.The monochromaticity was obtained by subtracting the energy resolution of the detector from the measurement results. The monochromaticity of the device in the measured energy range was below 3.24%, which is highly favorable and able to meet the calibration requirements of almost all detectors. The monochromaticity of the facility was then measured. Because the detector has energy resolution, the energy resolution of the measured all-round peak deducting the energy resolution of the detector itself is the monochromaticity of the monochromatic X-ray source.

3.2.3 Measurements of (20-161) keV X-rays

In this study, the double-crystal monochromator could cover the energy range of 20-161 keV using Bragg diffraction crystals Si(220) and Si(551). Several spectra measured by HPGe were selected, and the double-crystal monochromator was rotated so that the maximum diffraction angle could reach the lower limit of the energy,and the minimum diffraction angle could obtain the maximum energy. Owing to the limitation of the structure, the rotation range of the double-crystal monochromator for different crystals would be slightly different, and the energywould also be limited because of the influence of the crystal length and the distance between the two crystals.In theory, although various crystals have a significantly wide energy range,the theoretical range will be reduced because of the limitation of the processing structure. Each crystal had an excellent flux and energy resolution within the appropriate energy range.The measured spectra are shown in Fig. 16.

Fig. 15 (Color online)Measurement of the spot distribution

Table 7 The measurement results from the experiment

The flux stability was measured after the X-ray tube was warmed.Here,the flux changes within 10 h were recorded under certain X-ray tube voltages and currents. Two sets per hour were recorded for each set for 1000 s. The flux stability was greater than 0.8% over 10 h (Fig. 17).

Under normal circumstances, a uniform spot with good monochromaticity was obtained at less than 10 mm. The size of the light plate primarily depends on the size of the beam-limiting diaphragm. A CCD detector was used to measure the light spot under a beam-limiting diaphragm with a diameter of 4 mm, as shown in Fig. 18.

The HPGe detector was used to measure the performance of the device in more detail. The performance parameters, such as flux and monochromaticity, are listed in Table 8.In the present experimental mode,the minimum flux was greater than 2000 cm-2s-1. The monochromaticity was obtained by subtracting the energy resolution of the detector from the measurement results. The monochromaticity of the device in the measured energy range was below 3.55%.

Fig. 16 (Color online) The measured monochromatic spectra using a double-crystal monochro mator

Fig. 17 The flux stability at 127 keV

Fig. 18 (Color online)Measurement of light spot

3.2.4 Measurements of (21-301) keV X-rays

The monochromatic X-rays were measured using an HPGe detector, and the device was studied in detail. Several typical energy spectra are shown in Fig. 19. In fact,three Si crystals were used in the monochromator, and the generated monochromatic X-rays were related to different crystals. Ultimately, (21-301) keV was realized on this facility.

The detailed measurement parameters are listed in Table 9. The monochromaticity was better than 5.74%,whereas the minimum flux was greater than 2000 cm-2s-1. After detailed measurements, this facility was found to have good performance in terms of stability and a linear relationship between flux and cube current.The energy stability was better than 0.4% within 50 h,whereas the flux stability was better than 1.4%within 50 h.The said facility can be used to conduct X-ray diffraction studies and detector calibration studies.

Table 8 The specification of the double-crystal monochromator

Fig. 19 (Color online) The measured monochromatic spectra under Si(551) crystal

4 Applications

In recent years,several space science projects related to X-ray detection have been developed in China. Insight-HXMT was completed in 2017. High-energy detectors are the main components of the HXMT payload. The main detectors and backup detectors thereof have completed calibration of energy linearity, detection efficiency, and energy resolution on the present calibration device,and the energy range covers (20-150) keV [27]. The experimental energy resolution of the HED Z01-25 detector is shown in Fig. 20. The calibrated detection efficiencies of all 24 detectors at the present facility are presented in Fig. 21.The experimental calibration results verified the detection efficiency of the theoretical calculations.

Table 9 The specification of the high-energy monochromator

Fig. 20 (Color online) EC relations for a NaI detector

Fig. 21 (Color online) Detection efficiency experiments and MC calculation

Fig. 22 (Color online) E-FWHM relations of GRDs

Fig. 23 (Color online) Calibrated detection efficiency of some GRD detectors and MC simulation results

GECAM was launched in 2020.Its main load is a GRD detector that can measure photons in the 6 keV-5 MeV energy range. The core of the GRDs is the LaBr3crystal.Thus far,67 detectors (including 50 main detectors and 17 backup detectors) have been calibrated on these monochromatic facilities.The energy range covers(6-160)keV. Four of the detectors completed a fine calibration of the absorption edge (energy interval 0.1 keV). The experimental results are shown in Figs. 22 and 23.

In addition to the two satellite projects described above,the HXI, SVOM, and Gamma Ray Integrated Detectors(GRID) will be calibrated in due course. These projects have also undergone several preliminary calibration experiments.Detector research has also been conducted on these facilities, such as the calibration of a CdTe detector X-ray CCD detector.

5 Conclusion

Monochromatic X-ray facilities predicated on the X-ray tube and diffraction were introduced, and their detailed parameters were evaluated. The standard detectors were constructed through Monte Carlo simulations and experiments with a radioactive source with known activity. The(0.218-1.6) keV monochromatic X-ray beam facility was realized in a vacuum environment. Monochromatic X-rays have good monochromaticity and stability in the measurable energy range. The energy range was still relatively narrow and limited by the structure of the device. It was observed that an improved grating support structure or the use of a Bragg diffraction crystal to replace the grating could achieve(1-10)keV X-rays.For(5-301)keV energy range,several measurement and improvement studies have been conducted to provide stable detector testing and calibration services. The results show that there is a linear relationship between flux stability and the X-ray tube current. The experimental results also implied that the performance of the facilities was stable.

In the future, we will broaden the energy range to include both higher and lower energy as part of more detailed research on monochromatic X-ray facilities below 5 keV, to improve the calibration capability of this energy range. Further, attempts will be made to conduct research on new standard detectors such as XTES to create conditions for providing better calibration services. The authors are confident that an internationally renowned monochromatic X-ray calibration base can be built.

Author contributionsAll authors contributed to the study conception and design. Material preparation, data collection and analysis were performed by Si-Ming Guo, Jin-Jie Wu and Dong-Jie Hou. The first draft of the manuscript was written by Si-Ming Guo and all authors commented on previous versions of the manuscript. All authors read and approved the final manuscript.

主站蜘蛛池模板: 99久久国产自偷自偷免费一区| 日本高清免费不卡视频| 色精品视频| 一级成人欧美一区在线观看| 99精品视频九九精品| 波多野结衣的av一区二区三区| 免费国产高清视频| 日本免费新一区视频| 国产人妖视频一区在线观看| 国产精品吹潮在线观看中文| 日本国产精品| 免费无遮挡AV| 天天视频在线91频| 波多野结衣在线一区二区| 成年人国产网站| 在线观看亚洲天堂| 麻豆精品在线播放| 无码AV高清毛片中国一级毛片| 人人艹人人爽| 亚洲 欧美 偷自乱 图片| 国产欧美视频综合二区| 国产毛片网站| 国产极品粉嫩小泬免费看| 青青青国产精品国产精品美女| 尤物在线观看乱码| 国产欧美视频综合二区| 国产欧美在线观看一区| 久久人人97超碰人人澡爱香蕉 | 久久黄色一级视频| www.youjizz.com久久| 久久国产高潮流白浆免费观看| 免费观看国产小粉嫩喷水| 老司机久久99久久精品播放| 免费毛片网站在线观看| 老司机久久精品视频| 永久在线精品免费视频观看| 精品国产电影久久九九| 少妇人妻无码首页| 在线看国产精品| 国产精品片在线观看手机版| 国产成人乱码一区二区三区在线| 亚洲另类色| 欧美一级在线看| 91尤物国产尤物福利在线| 国产午夜无码专区喷水| 不卡的在线视频免费观看| 久久伊人久久亚洲综合| 91在线激情在线观看| 色偷偷一区| 亚洲成aⅴ人片在线影院八| 在线免费看片a| 久996视频精品免费观看| 国产AV毛片| 中文字幕亚洲精品2页| 久久久噜噜噜| 国产精品成人观看视频国产| 国产亚洲精品97AA片在线播放| 亚洲国产av无码综合原创国产| 国产乱人激情H在线观看| 免费久久一级欧美特大黄| 青青国产成人免费精品视频| 波多野结衣在线一区二区| 亚洲欧美色中文字幕| 亚洲天堂网在线视频| 三区在线视频| 欧美成人看片一区二区三区| 天堂在线视频精品| 国产成熟女人性满足视频| 亚洲精品第1页| 伊人久久大香线蕉aⅴ色| 亚洲有无码中文网| 女同久久精品国产99国| 小蝌蚪亚洲精品国产| 欧美另类视频一区二区三区| 日本精品影院| 亚洲欧洲日产国产无码AV| 国产精品久久久久无码网站| 亚洲色图在线观看| 99这里只有精品免费视频| 国产成人一区在线播放| 亚洲精品第一页不卡| 大乳丰满人妻中文字幕日本|