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

    Performance optimization of scintillator neutron detectors for EMD in CSNS

    2023-12-02 09:22:22XiaojieCai蔡小杰QianYu于潛ChangHuang黃暢BinTang唐彬ShihuiZhou周詩慧XiaohuWang王小胡XiupingYue岳秀萍andZhijiaSun孫志嘉
    Chinese Physics B 2023年11期

    Xiaojie Cai(蔡小杰), Qian Yu(于潛), Chang Huang(黃暢), Bin Tang(唐彬),4,?, Shihui Zhou(周詩慧),Xiaohu Wang(王小胡), Xiuping Yue(岳秀萍), and Zhijia Sun(孫志嘉),4,§

    1Spallation Neutron Source Science Center,Dongguan 523803,China

    2Southwest University of Science and Technology,Mianyang 621002,China

    3Institute of Modern Physics,Chinese Academy of Sciences,Lanzhou 730000,China

    4Institute of High Energy Physics,Chinese Academy of Sciences,Beijing 100049,China

    5Shenzhen Key Laboratory of Ultraintense Laser and Advanced Material Technology,Center for Advanced Material Diagnostic Technology,and College of Engineering Physics,Shenzhen Technology University,Shenzhen 518118,China

    Keywords: scintillator neutron detector,silicon photomultipliers,wavelength shifting fiber, 6LiF/ZnS scintillator,detection efficiency inhomogeneity

    1.Introduction

    With the rapid development of neutron scattering technology,the scattering technology based on high intensity source is gradually playing an irreplaceable important role in the fields of material science and life science.

    The spallation neutron source produced by acceleratordriven charged particles (such as protons, etc.) can provide high effective neutron flux, superior pulse time structure,and rich high-energy short-wavelength neutron combined with time-of-flight technology.[1]Engineering neutron diffraction can nondestructively and noninvasively probe stress,strain, temperature, and phase evolutions deep within bulk materials.[2]The mechanical behavior of real time can be achieved by applying temperature and tensile environment to the sample through the load platform,and thus provide strong technical support for manufacturing process, failure analysis,and so on.Chinese spallation neutron source is one of the four largest pulse spallation neutron sources in the world.Engineering material diffractometer (EMD), which can realize high-precisionin-situstress and strain detections of materials, will be constructed at the Beamline 8 between 2019 and 2023.It will serve for deep stress detection and dynamic change research of materials.A novel thermal neutron detector has been developed for the requirements of EMD applications.The detector consists of6LiF/ZnS(Ag)scintillation screens,wavelength-shifting fiber(WLSF)arrays silicon photomultipliers(SiPMs), and application specific integrated circuit(ASIC)read-out electronics.There are 14 detector modules with nearly 3 m2covered area has been assembled in two 90?neutron scattering positions in the scattering room of EMD phase I.Before their installation, the performances of all the detector modules need to be tested.[3]Due to the existence of various influencing factors in the production process,the performances between each detector modules and between each channel in one detector are not consistent.The inhomogeneity detection efficiencies between each pixel and each detector are not consistent.[4]The inhomogeneity mainly caused by the inconsistency of the WLSFs in collecting scintillation photons,needs to be mitigated before the detector installation.In this work,the detection efficiency and its inhomogeneity of the detector module has been tested in BL.20 in CSNS.By adjusting the readout electronics threshold of each pixel, the detection efficiency inhomogeneity has been improved.A standard Si powder sample was also used to prove the diffraction peaks achieved by all the optimized detector modules,and it worked correctly.

    2.Working principle and structure of scintillator detector

    The scintillator neutron detectors for neutron spectrometers have been developed in CSNS since 2011.Version I scintillator detectors with6LiF/ZnS(Ag) screens, wavelength shifting fiber (WLSF) arrays, and multi-anode photo multiplier tubes (MA-PMTs) have been successfully installed and operated well in general purpose of powder diffractometer(GPPD)in CSNS.[5]The detector adopts the upper and lower double scintillation screens intermediately sandwiched for two-dimensional wave shifted fiber array structure, as shown in Fig.1.

    Fig.1.Scintillation detector unit structure diagram.

    In such detectors, the6LiF/ZnS(Ag) screens convert the incident neutrons into secondary charged particles which cause the scintillation photons in the following reaction:

    n+6Li→α(2.05 MeV)+3H(2.73 MeV).

    The purpose of this design is to increase the thickness of the scintillation screen and increase the probability of neutrons reacting with6Li nuclide in the scintillation screen to release more photons.[5,6]

    However,the EMD requires higher demand for detection efficiencies.The detection efficiency of the scintillation screen depends on the probability of nuclear reactions between neutrons and6Li nuclide, and photon emission probability in the scintillation screen.Therefore,we designed the detector with oblique incident shutter structure.[7]

    Version II scintillator detector developed for the EMD application achieves higher detection efficiencies.This detector module is composed of inclined6LiF/ZnS(Ag)scintillation screens (scintillator/ND/300-μm thickness with6LiF:ZnS mass ratio 1:2), wavelength-shifting fiber array (WLSF, Kuraray/Y11(200)M/round?1 mm with 200 ppm wave shifting doping), the silicon photomultiparray(SiPM,SENSL/miroFC-30035-SMT), and the ASIC readout electronics.The fluorescing light is collected by the WLSF which absorbs blue light and reemits green light.The green light was transmitted to the end surface of the WLSF,which is coupled to the SiPMs by silicon oil.SiPMs convert the optical signals into pulsed electrical signals.The signals are amplified, shaped, and discriminated by custom-designed read-out ASIC electronics and converts them into digital signals which will be input to the FPGA chip for further processing.[8]Different from the traditional discrete component electronics,the ASIC has the advantages of high integration and low power consumption,which is more suitable for engineering applications of multi-pixel detectors.As shown in Fig.2,the head of the detector adopts oblique incident shutter structure.It could improve the detect efficiency effectively because the opaque6LiF/ZnS(Ag) scintillation screen makes the incident path of the neutron larger, while the outgoing path of the fluoresce light is not changed.The tilt angle of the scintillation screen is designed as 17?to the incident neutron.

    The WLSF arrays distributed longitudinally on the surface of the scintillation screen to collect the fluorescence from the neutron screens.[7]

    Fig.2.The schematic diagram of the detector head.

    In order to verify the fact that the optimal detection efficiency can be achieved when the scintillation screen of the detector is tilted at an angle of 17?.We used GEANT4 to simulate the detector structure of the scintillation screen at different tilt angles.Figure 3 is the neutron conversion efficiency of the scintillation screen under various tilt angles at different wavelengths.The results show that with the decrease of the tilt angle between the scintillation screen and incident neutrons,the detection efficiency of the scintillation screen increases gradually.Considering that the pixel size is 3 mm and the space requirement of the WLSF arrangement, the inclination angle of the scintillation screen of the detector is designed as 17?.The neutron reaction efficiency of the scintillation screen could reach 81.2%@1 ?A neutron wavelength.

    Fig.3.The neutron reaction efficiency of the scintillation screen under various tilt angles at different wavelengths.

    The hade of the detector module is designed as a circular arc with a diameter of 2 m to realize circular neutron detection detector cell with common mechanical chamber for supporting and light shielding.Figure 4 shows the photo of the detector cell.Five detector modules are also fixed to the mechanical support frame to form an arc-shaped construction with the radius of 2 m.Except the neutron incident window,the out-shell of the cell adopts boron-containing aluminum plates to protect it from neutron radiations.Each module in the cell consists of 643 mm×180 mm size pixels.The cell dimensions are 1066 mm (length)×203 mm (height)×450 mm (thickness).All the detector cells are mounted on the bank bracket according to the laser positioning data.

    Fig.4.Photo of a detector cell consist of 5 detector modules.

    3.Performance test of the detector module

    According to the physical design of EMD, the main parameters of the scintillator detector are detection efficiency and detection efficiency uniformity.The consistency of position resolution of the detectors is mainly determined by the accuracy of the mechanical processing and the optical isolation for each pixel.These differences are usually small.However the detection efficiency performances for each of detector modules and for each of pixels are not consistent.It mainly related to the uniformity of tilt angle and the thickness of the scintillation screen, the uniformity of the optical transmission efficiency of each WLSF from its bending and coupling process,and the uniformity of the gain of each SiPM.The uniformity of the thickness of the scintillation screen is less than±5%according to the examination report from the manufacturer.The difference of the tilt angle is also less than 1?(±6%)by the finish machining of the mechanical support.The gain difference of each of SiPMs is less than±5%according to the manufacturer quality report.The uniformity of the optical transmission efficiency of each WLSF from its bending and coupling process is less than±10%with the manufacturing process quality control.All these differences may cause a maximum difference of 1.6 times inhomogeneity of the detect efficiency for each of pixels.So before the detector modules installation,the inhomogeneity needs to be measured and mitigated.

    The typical way for the calibration of the detection efficiency is using isotopic neutron source with 4πemitting uniformly.[9,10]

    It will cost quite long time for all the detector pixel detections to gain a good efficiency uniformity.The time schedule of the EMD installation is not allowed us to do such an experiment.

    The BL.20 neutron beam line of the CSNS is designed specifically for the detector testing.We use high incoherent neutron scattering sample to disperse the neutron beam and the neutron counts of each pixel.By adjusting the readout electronics threshold of each pixel, the inhomogeneity of the neutron detect efficiency can be improved.

    3.1.Experimental setup

    The BL.20 beam line in CSNS uses a decoupled narrowing liquid hydrogen moderator.After several collimations and absorption,the diameter of the neutron beam spots at the outlet is 20 mm.The neutron fluency is about 108n/(cm2·s).The neutron wavelength range is 0.2 ?A–10 ?A and can be changed by two chopper with 5.6-?A bandwidth.The distance from moderator to the neutron outlet spot is about 9 m.Figure 5 shows the structure of the BL.20 beam line in CSNS.

    The measurement setup for measuring relative detection efficiency is shown in Fig.6.The standard3He tube(LND252299/pressure 20 atm,1 atm=1.01×105Pa)/diameter 1 inch,1 inch=2.54 cm)is placed to evaluate the incident neutron flux.

    The outgoing neutrons pass through the beam limiting hole from the BL.20 beam line.The limited neutron beam then incident into the mica monochromator with an incident angle of 45?.Because the distribution of the neutron wavelength from the BL20 is wide and the neutron flux is still too high for the3He tube,the mica(0 0 10)(0 0 8)monochromator is used to further reduce the neutron flux to prevent the count rate saturation of the3He tube and to choose special wavelengths of neutrons (wavelength 1.4 ?A, 2.8 ?A, and 4.8 ?A) for calculating the detect efficiency accurately.in the beam exit.The detect efficiency can be calculated by the following formula:

    Fig.5.The structure of the BL.20 beam line in CSNS.

    Fig.6.Detection efficiency measurement setup diagram.

    The time-of-flight(TOF)spectrum of the monochromatic neutrons is gained by the standard3He tube and the detector respectively in the same position.Both two neutron counts are measured normally by the counts of the beam monitor placed

    whereNSDandN3Heare normalized neutron counts of scintillator detectors and standard3He tube respectively,andη3Heis the detection efficiency of standard3He tube detectors.

    Using high incoherent neutron scattering sample, such as water or polyethylene, to disperse the neutron beam, the neutron counts of each pixel in each detector cell could be obtained.[11]These counts reflect the inhomogeneity of the neutron detection efficiency for each of pixels.The detector cells are placed in an aluminum bracket with 2-m distance from the sample.

    Fig.7.Field experimental test of imitating engineering material spectrometer.

    And the angle between the center of the neutron beam line and that of the detector cell is 23?,which is the maximum angle that the scattering room of the BL.20 beam line can be placed.The inhomogeneous area of the neutron scattering in forward direction for water sample needs to be avoid.The aluminum bracket can place 7 detector cells trying to simulate the actual placement of the detector cells in EMD.The schematic diagram of the experimental set-up is shown in Fig.7.

    A boron-containing rubber plate collimating tube with 1.16-m long,10-cm diameter,and 5-mm thickness was placed at 11.4-m position from the beam hole.The collimating tube reduces the neutron divergence.A slit with 8-mm width was set in front of water sample.Finally, the diffracted neutrons incident into the detectors in 2-m arc.And in order to reduce the spatial neutron background, all the tables were wrapped with boron-containing polyethylene.Before measurements,the uniformity of the moderated neutrons after 2-cm thickness water with 10 cm×10 cm area was simply simulated by GEAN4.A cured ideal neutron detector with 6 m×6 m area was placed on a 2-m arc.The moderated neutrons counting distribution is approximately uniform in a range of angles 10?–25?between the neutron incident direction and the ideal detector.

    Fig.8.GEANT4 simulation results.

    The diffraction peak of the stander sample could be calculated through the data of the detector by changing water to stander Si powder in the later stage.

    3.2.Measurement results and discussion

    Using the mica monochromator, the counts for neutron wavelengths of 2.9 ?A and 1.4 ?A for the standard3He tube and that of the scintillator detector were obtained.And the inner wall is fitted with a 5-mm polyethylene block collimated pipe being placed at 11.4-m position from the beam hole.The collimating tube reduces the neutron divergence.A slit with 8-mm width was set in front of the water sample.Finally, the diffracted neutrons incident into the detectors with 2-m distance.And in order to reduce the spatial neutron background,all the tables was wrapped with boron-containing polyethylene.The diffraction peak of the stander sample could be calculated through the data of the detector by changing the water to stander Si powder in the later stage respectively.

    3.2.1.Detection efficiency

    The incident neutron wavelength can be calculated by the TOF(time of flight)way using the following formula:

    whereLis the neutron flight distance,mnis neutron mass,andhis the Planck constant.The neutron wavelength spectrum can be calculated with this formula.

    The TOF spectrum of the standard3He tube and the particularly pixel of the detector is calculated and is shown in Fig.9(a).We find that the two characteristic peaks with wavelengths 1.4 ?A,2.8 ?A,and 4.8 ?A are clear.However,in the short wavelength part(1 ?A–2 ?A),the spatial neutron background is more serious in the TOF spectrum.The reason is that part of the short wavelength neutrons background cannot shielded effectively by the boron-containing polyethylene.

    By calculating the total number of neutrons in the specific wavelength segment,such as 1.38 ?A–1.42 ?A or 4.78 ?A–4.82 ?A,the detection efficiency of the corresponding wavelength of the detector can be obtained by using formula (1).the detection efficiency of the scintillator detector is 53%@1.4 ?A,77%@2.8 ?A,and 88%@4.8 ?A,respectively.

    Figure 9(b) is the detection efficiency curve for the particularly pixel of the detector.The difference in detection efficiency between the three pixels can be seen from the picture.This will affect the accuracy of the test result of samples in the engineering application of the stress testing.

    Therefore, the detection efficiency of the detector needs to be optimized.In the later stage, the detection efficiency of each pixel of the detector needs to be modified to further reduce the influence of the non-uniformity of detection efficiency on the stress testing result of the sample.

    Fig.9.TOF spectra obtained by the standard 3He tube and the prototype scintillator detector(a).The detection efficiency of the scintillator detector as a function of neutron wavelength(b).

    3.2.2.Detection efficiency uniformity

    According to the former test of the neutron signal charge spectrum(see Fig.10)of the scintillator detector module with the252Cf radioactive source, the signal spectrum is a wide Landau distribution.This is caused by the second particle energy loss distribution of the6Li(n, d)αnuclear reaction, the opacity of the6LiF/ZnS(Ag) scintillator, and the optical collection efficiency difference of the round WLSFs for different incident scintillate photon positions.When the electronics threshold is setting between 1 pc and 5 pc,the neutron signals can be separated from the background noise,and the detection efficiency can maintain in 40%–60%@1.4 ?A.

    Fig.10.The neutron signal spectrum of the scintillator detector.

    This means that by adjusting the electronic thresholds of each channel of the detector,it is possible to optimize the uniformity of the detector efficiency of the scintillator detectors.For stress measurement, the inhomogeneity of detection efficiency of the detector cannot guarantee the accurate measurement of continuous stress of the sample.Only with better uniformity of the detection efficiency can the weak stress variation of samples be detected.[12]

    Before the neutron uniformity measurement,the standard detector MWPC is used to test the spatial uniformity of neutrons.By this testing,the moderated neutrons counting distribution is approximately uniform.Thus the neutrons are generally uniformly distributed.[5]

    Fig.11.A comparison of neutron counts for each channel changed before and after adjusting electronics.

    By adjusting the threshold of each channel to optimize the neutron counts uniformity of the detector,the uniformity of the neutron counts for each pixel are greatly improved.Figure 11 shows the neutron counts distribution of a detector cell before and after the adjusting.The original neutron counts distribution showed that the detect efficiency non-uniformity caused by the device difference is obvious.The way of adjusting the electronics threshold only optimizes the detect efficiency in a limiting range.If the difference caused by the wrong operation in manufacturing process or caused by the bad device performance is too large, the detect efficiency homogeneity of the pixel is hard to be optimized to an acceptable level,such as 3 pixels shown in Fig.11.In general,the efficiency uniformity of the detector can be optimized to be better than 90%.

    Figure 12 shows the neutron counts distribution of all 14 detector cells after electronics threshold adjusting.We can see some bad channels in several cells, such as pixel 128, pixel 192, and pixel 256 still had bad performance.It is mainly caused by the WLSF break in the processing according to our return inspection.

    Fig.12.The uniformity of fourteen detectors after calibration.

    The WLSF breakage will cause most of the light collected from the scintillation screen hardly being transmitted to the SiPM, which will result in a small signal in this channel.It is difficult to improve the neutron detection efficiency of these kinds of channels by adjusting the electronic threshold.

    3.2.3.Standard Si sample testing

    After the detect efficiency optimization tests for 14 detector cells, the standard Si powder samples are used to check the performance of all the detector cells.Figure 13 shows the diffraction spectrum of the standard Si powder sample obtained by the optimized detector.According to the Bragg’s diffraction formula

    whereθis the angle between the incident direction of the neutron and the diffracted neutron.λis the wavelength of the neutron.Theλvalue can be obtained from the time-of-flight spectrum,and thus thedspectrum of the detector pixel can be obtained.

    As can be seen from Fig.13, multiple diffraction peaks appear.By comparing the test results with the diffraction spectrum of the Si powder sample from the scintillator detectors of the GPPD in CSNS,it is shown that the numbers and the positions of the diffraction peaks are consistent.However, the results are interfered by the short wavelength neutrons background from the Si powder incoherent scattering.

    Fig.13.The diffraction spectrum of the standard Si powder sample obtained by the optimized detector.

    4.Conclusion and outlook

    By the testing performance of detectors,the detection efficiency of the detector is increased more than 40% at the wavelength 1.4 ?A,which meets the design of the EMD spectrometer.The detect efficiency uniformity caused by small different fabrication processes of each independent devices of the detector cells is also measured by water sample experiment.The detection efficiency of each pixel channel of the detectors is optimized to be consistent by adjusting the electronics thresholds.The detection efficiency of several channels did not adjust to the better state because of the WLSFs break.After the optimization,the mean square error is less than 5%,and there is no obvious uneven interval.

    Fig.14.The EMD spectrometer real scene picture.

    The performances of the testing detectors show that all the 14 fabricated detectors can be used in engineering applications in EMD.These new scintillator detectors not only meet the urgent needs of CSNS engineering, but also improve the level of neutron detection technology in our country.The EMD has great significance in the material science study.The EMD spectrometer real scence picture is shown in Fig.14.

    The EMD spectrometer can provide domestic highperformance neutron detectors for many advanced neutron scattering experimental platforms in China, and at the same time, it can also serve as the stress measurement for material design and process upgrading and improvement.

    The establishment of engineering material diffraction spectrometer will bring great help to solve the problems encountered by Chinese industry,especially to improve the level of China’s aerospace industry.

    Acknowledgements

    Project supported by the National Natural Science Foundation of China (Grant Nos.11975255 and 11875273) and Guangdong Basic and Applied Basic Research Foundation(Grant No.2020B1515120025).

    各种免费的搞黄视频| 亚洲国产看品久久| 黄网站色视频无遮挡免费观看| 变态另类成人亚洲欧美熟女 | 亚洲一区中文字幕在线| 激情在线观看视频在线高清| 黑人巨大精品欧美一区二区mp4| 亚洲全国av大片| 丰满的人妻完整版| 国产精品自产拍在线观看55亚洲| av网站免费在线观看视频| 亚洲五月婷婷丁香| 久久久久国产一级毛片高清牌| 精品国产乱子伦一区二区三区| 可以在线观看毛片的网站| 国产精品国产高清国产av| 老汉色∧v一级毛片| 欧美在线黄色| 在线天堂中文资源库| 日日干狠狠操夜夜爽| 久久亚洲真实| 人成视频在线观看免费观看| 1024视频免费在线观看| 99riav亚洲国产免费| 欧美人与性动交α欧美精品济南到| 老司机午夜十八禁免费视频| 久久热在线av| 看黄色毛片网站| 手机成人av网站| 在线观看www视频免费| 日韩精品青青久久久久久| 免费少妇av软件| 精品国产美女av久久久久小说| 日本精品一区二区三区蜜桃| 亚洲av电影在线进入| 国产熟女午夜一区二区三区| 日日摸夜夜添夜夜添小说| 999精品在线视频| 另类亚洲欧美激情| 成人手机av| 成人18禁高潮啪啪吃奶动态图| 国产精品免费视频内射| 深夜精品福利| 欧美日韩一级在线毛片| 中文字幕高清在线视频| 精品久久久久久,| 国产av精品麻豆| 国产麻豆69| 热re99久久国产66热| 欧美+亚洲+日韩+国产| 19禁男女啪啪无遮挡网站| 国产熟女午夜一区二区三区| 免费在线观看日本一区| 波多野结衣高清无吗| 午夜福利免费观看在线| 亚洲免费av在线视频| 久久精品国产综合久久久| 亚洲av成人不卡在线观看播放网| 亚洲精品国产精品久久久不卡| 99香蕉大伊视频| 黄色毛片三级朝国网站| 韩国精品一区二区三区| 性色av乱码一区二区三区2| 中文字幕色久视频| 亚洲午夜精品一区,二区,三区| 国产不卡一卡二| 69av精品久久久久久| 人妻久久中文字幕网| 热99国产精品久久久久久7| 美女午夜性视频免费| 高清黄色对白视频在线免费看| 大码成人一级视频| 少妇 在线观看| 他把我摸到了高潮在线观看| 欧美人与性动交α欧美精品济南到| 啦啦啦在线免费观看视频4| 国产欧美日韩一区二区三区在线| 少妇裸体淫交视频免费看高清 | 亚洲av五月六月丁香网| 99热只有精品国产| 国产成人一区二区三区免费视频网站| 亚洲自偷自拍图片 自拍| 他把我摸到了高潮在线观看| 又大又爽又粗| 丁香六月欧美| 婷婷精品国产亚洲av在线| 成人亚洲精品一区在线观看| 久久精品亚洲av国产电影网| av网站在线播放免费| 国产精品1区2区在线观看.| 国产欧美日韩一区二区三| a在线观看视频网站| 最近最新中文字幕大全电影3 | 午夜免费鲁丝| 黑人巨大精品欧美一区二区蜜桃| 日韩av在线大香蕉| 久久久久久久久免费视频了| 日本五十路高清| 免费观看人在逋| 成人三级做爰电影| 国产99白浆流出| 丰满的人妻完整版| 国产麻豆69| 超碰成人久久| 亚洲,欧美精品.| 亚洲色图 男人天堂 中文字幕| 欧美最黄视频在线播放免费 | 欧美乱妇无乱码| 脱女人内裤的视频| 亚洲国产精品999在线| 69av精品久久久久久| 级片在线观看| 波多野结衣高清无吗| 性色av乱码一区二区三区2| 国产激情久久老熟女| 欧美成人午夜精品| 大陆偷拍与自拍| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品一区二区精品视频观看| 99在线人妻在线中文字幕| 国产麻豆69| 免费日韩欧美在线观看| 亚洲视频免费观看视频| 精品一区二区三区四区五区乱码| 高清欧美精品videossex| 午夜a级毛片| 免费看十八禁软件| 精品国产国语对白av| 欧美在线黄色| 免费在线观看黄色视频的| 黑人操中国人逼视频| 巨乳人妻的诱惑在线观看| 亚洲一区二区三区欧美精品| 色播在线永久视频| 国产又色又爽无遮挡免费看| 久久性视频一级片| 99国产精品一区二区三区| 丰满的人妻完整版| av片东京热男人的天堂| 国产精品98久久久久久宅男小说| 久久久久久久午夜电影 | 欧美老熟妇乱子伦牲交| 国产区一区二久久| 曰老女人黄片| 99精品欧美一区二区三区四区| 视频区图区小说| 老司机深夜福利视频在线观看| 黄色女人牲交| x7x7x7水蜜桃| 一区二区三区激情视频| 天堂中文最新版在线下载| 级片在线观看| 久久久国产欧美日韩av| 久久午夜亚洲精品久久| 精品卡一卡二卡四卡免费| 国产真人三级小视频在线观看| 久久久久九九精品影院| 成人永久免费在线观看视频| 中文字幕人妻丝袜一区二区| 人妻久久中文字幕网| 久久精品91蜜桃| 亚洲狠狠婷婷综合久久图片| 丰满人妻熟妇乱又伦精品不卡| 老汉色∧v一级毛片| 香蕉丝袜av| 黄色丝袜av网址大全| 亚洲成人久久性| 欧美激情高清一区二区三区| 怎么达到女性高潮| 午夜a级毛片| 男人操女人黄网站| 一边摸一边抽搐一进一出视频| 日本精品一区二区三区蜜桃| 国产激情久久老熟女| 999精品在线视频| 日韩欧美免费精品| 亚洲精品久久午夜乱码| 久久亚洲精品不卡| 日本撒尿小便嘘嘘汇集6| 老熟妇乱子伦视频在线观看| 免费高清在线观看日韩| 人人澡人人妻人| 免费看a级黄色片| 国产精品久久视频播放| 韩国av一区二区三区四区| 午夜福利一区二区在线看| 99热国产这里只有精品6| 亚洲va日本ⅴa欧美va伊人久久| 午夜免费观看网址| 99久久综合精品五月天人人| 97碰自拍视频| 久久久久国内视频| 国产伦一二天堂av在线观看| 长腿黑丝高跟| 多毛熟女@视频| 好男人电影高清在线观看| 久久99一区二区三区| 人人妻人人添人人爽欧美一区卜| 国产又爽黄色视频| 精品国产一区二区久久| 亚洲精品一卡2卡三卡4卡5卡| 国产一区二区在线av高清观看| 最新在线观看一区二区三区| 午夜视频精品福利| 高清毛片免费观看视频网站 | 久久久水蜜桃国产精品网| 首页视频小说图片口味搜索| 麻豆久久精品国产亚洲av | 天堂俺去俺来也www色官网| 国产精品爽爽va在线观看网站 | 午夜免费激情av| 久久精品国产综合久久久| 国产精品国产av在线观看| 精品无人区乱码1区二区| 国产乱人伦免费视频| 国产精品九九99| 夫妻午夜视频| 人人妻,人人澡人人爽秒播| 麻豆av在线久日| 丝袜美足系列| 老汉色av国产亚洲站长工具| av天堂久久9| 亚洲成人精品中文字幕电影 | 一二三四在线观看免费中文在| 久久久国产精品麻豆| 国产片内射在线| 99国产极品粉嫩在线观看| 在线观看午夜福利视频| 亚洲午夜精品一区,二区,三区| 久久久国产精品麻豆| 午夜久久久在线观看| 男人舔女人的私密视频| 久久九九热精品免费| 日韩国内少妇激情av| 亚洲精品美女久久av网站| 精品午夜福利视频在线观看一区| 日本五十路高清| 中文欧美无线码| 老司机亚洲免费影院| 国产午夜精品久久久久久| 免费久久久久久久精品成人欧美视频| 日韩成人在线观看一区二区三区| 日本一区二区免费在线视频| 亚洲精品成人av观看孕妇| 国产精品自产拍在线观看55亚洲| xxxhd国产人妻xxx| 午夜福利,免费看| 成人三级做爰电影| 久久久久亚洲av毛片大全| 老鸭窝网址在线观看| 亚洲av成人一区二区三| 国产亚洲欧美在线一区二区| 人妻丰满熟妇av一区二区三区| 国产又色又爽无遮挡免费看| 成年女人毛片免费观看观看9| 纯流量卡能插随身wifi吗| 亚洲一区二区三区欧美精品| 亚洲av五月六月丁香网| 18禁美女被吸乳视频| 91精品国产国语对白视频| 18禁国产床啪视频网站| 男女高潮啪啪啪动态图| 午夜福利欧美成人| 少妇的丰满在线观看| 免费在线观看影片大全网站| 久久九九热精品免费| 老司机在亚洲福利影院| 日日干狠狠操夜夜爽| 色综合婷婷激情| 91精品三级在线观看| 日韩欧美在线二视频| 男女午夜视频在线观看| 久久这里只有精品19| 国产99白浆流出| 一级片免费观看大全| 美国免费a级毛片| 久久精品影院6| 在线播放国产精品三级| 欧美黑人精品巨大| 美女扒开内裤让男人捅视频| 国产真人三级小视频在线观看| x7x7x7水蜜桃| 神马国产精品三级电影在线观看 | 最好的美女福利视频网| 久久这里只有精品19| 色哟哟哟哟哟哟| av超薄肉色丝袜交足视频| 日韩免费高清中文字幕av| 美女高潮到喷水免费观看| 久久人妻av系列| 国产高清激情床上av| 女人被狂操c到高潮| 纯流量卡能插随身wifi吗| 国产精华一区二区三区| 一区二区三区激情视频| 男人的好看免费观看在线视频 | 欧美日韩一级在线毛片| 欧美日韩亚洲高清精品| 99在线视频只有这里精品首页| 亚洲专区国产一区二区| 久久久久久久久久久久大奶| 老熟妇乱子伦视频在线观看| 人成视频在线观看免费观看| 亚洲av五月六月丁香网| 日韩av在线大香蕉| 亚洲av五月六月丁香网| 夜夜爽天天搞| 国产精品98久久久久久宅男小说| 人成视频在线观看免费观看| 精品福利观看| 母亲3免费完整高清在线观看| 欧美乱色亚洲激情| 午夜亚洲福利在线播放| 超碰成人久久| 高清黄色对白视频在线免费看| 一本综合久久免费| 欧美日韩乱码在线| 国产伦一二天堂av在线观看| 黄频高清免费视频| 不卡av一区二区三区| 叶爱在线成人免费视频播放| 成人永久免费在线观看视频| 免费女性裸体啪啪无遮挡网站| 色综合婷婷激情| 欧美激情极品国产一区二区三区| 亚洲精品久久午夜乱码| 国产黄a三级三级三级人| 黑人欧美特级aaaaaa片| 国产成年人精品一区二区 | 亚洲av电影在线进入| 国产精品野战在线观看 | 欧美性长视频在线观看| а√天堂www在线а√下载| 搡老熟女国产l中国老女人| 久久久国产精品麻豆| 日日爽夜夜爽网站| 一二三四社区在线视频社区8| 9色porny在线观看| 国产精品免费一区二区三区在线| 天堂影院成人在线观看| 亚洲,欧美精品.| 久久午夜综合久久蜜桃| 99热只有精品国产| 亚洲精品久久午夜乱码| 国产午夜精品久久久久久| 人妻丰满熟妇av一区二区三区| 深夜精品福利| 亚洲 欧美一区二区三区| 国产成年人精品一区二区 | 国产av精品麻豆| 精品卡一卡二卡四卡免费| 视频区欧美日本亚洲| 麻豆一二三区av精品| 欧美日本中文国产一区发布| 国产99久久九九免费精品| 88av欧美| 亚洲精品av麻豆狂野| 99在线人妻在线中文字幕| 久久青草综合色| 一二三四社区在线视频社区8| 黄色a级毛片大全视频| 亚洲第一青青草原| 人人妻人人添人人爽欧美一区卜| 亚洲精品国产精品久久久不卡| av国产精品久久久久影院| 天堂动漫精品| 久久国产亚洲av麻豆专区| 国产一区二区三区综合在线观看| 国产成人av激情在线播放| 亚洲久久久国产精品| 精品久久久久久久久久免费视频 | 国产精品永久免费网站| 亚洲精品一二三| 丰满迷人的少妇在线观看| 18禁国产床啪视频网站| 亚洲精品久久午夜乱码| 久久精品人人爽人人爽视色| 91九色精品人成在线观看| 黄色视频,在线免费观看| 咕卡用的链子| 久久99一区二区三区| 亚洲午夜精品一区,二区,三区| 人成视频在线观看免费观看| 国产精品1区2区在线观看.| 麻豆av在线久日| 亚洲欧美激情在线| 成人黄色视频免费在线看| 老汉色∧v一级毛片| av福利片在线| av天堂在线播放| 成人国语在线视频| 久久狼人影院| 亚洲欧洲精品一区二区精品久久久| 国产精品国产av在线观看| 欧美大码av| 88av欧美| 热re99久久国产66热| 俄罗斯特黄特色一大片| 国产精品国产高清国产av| 电影成人av| 亚洲五月色婷婷综合| 大香蕉久久成人网| 一进一出抽搐gif免费好疼 | av视频免费观看在线观看| www.熟女人妻精品国产| 91精品三级在线观看| 可以在线观看毛片的网站| 久久精品人人爽人人爽视色| 国产精品久久久久成人av| 亚洲国产精品一区二区三区在线| 侵犯人妻中文字幕一二三四区| a在线观看视频网站| 亚洲一区二区三区色噜噜 | 日韩中文字幕欧美一区二区| 99香蕉大伊视频| 人人妻,人人澡人人爽秒播| av片东京热男人的天堂| 国产精品香港三级国产av潘金莲| 午夜成年电影在线免费观看| 岛国视频午夜一区免费看| 亚洲 欧美一区二区三区| videosex国产| 欧美国产精品va在线观看不卡| 免费人成视频x8x8入口观看| 亚洲熟女毛片儿| 波多野结衣一区麻豆| 亚洲全国av大片| 桃色一区二区三区在线观看| 日本 av在线| 精品久久久久久电影网| 国产99白浆流出| www.www免费av| 在线观看66精品国产| 五月开心婷婷网| 岛国视频午夜一区免费看| 婷婷丁香在线五月| 纯流量卡能插随身wifi吗| 91国产中文字幕| 亚洲欧美日韩另类电影网站| 久久久久久久精品吃奶| 午夜日韩欧美国产| 亚洲成av片中文字幕在线观看| 国产精品 欧美亚洲| 亚洲第一av免费看| 性少妇av在线| 999久久久国产精品视频| 亚洲色图 男人天堂 中文字幕| 欧美成人午夜精品| 亚洲国产欧美一区二区综合| 两性午夜刺激爽爽歪歪视频在线观看 | xxxhd国产人妻xxx| 宅男免费午夜| 丰满人妻熟妇乱又伦精品不卡| 欧美在线黄色| 久久久久久亚洲精品国产蜜桃av| 91字幕亚洲| 午夜精品在线福利| 国产男靠女视频免费网站| 老司机福利观看| 亚洲成人国产一区在线观看| 看黄色毛片网站| 亚洲精品中文字幕一二三四区| 亚洲欧美日韩另类电影网站| 日本一区二区免费在线视频| 精品国产一区二区三区四区第35| 精品一区二区三区av网在线观看| 亚洲av电影在线进入| 国产精品国产高清国产av| 男人的好看免费观看在线视频 | 亚洲人成电影免费在线| 亚洲国产精品sss在线观看 | 日韩精品免费视频一区二区三区| 欧美久久黑人一区二区| 日韩视频一区二区在线观看| 9热在线视频观看99| 久久午夜亚洲精品久久| 免费看a级黄色片| 日韩欧美在线二视频| 天堂动漫精品| 一本大道久久a久久精品| 久久天躁狠狠躁夜夜2o2o| 亚洲一区二区三区色噜噜 | 三上悠亚av全集在线观看| 极品人妻少妇av视频| 一a级毛片在线观看| 精品一区二区三区av网在线观看| 99国产综合亚洲精品| 精品福利观看| 国产欧美日韩综合在线一区二区| 成人三级黄色视频| 丰满迷人的少妇在线观看| 精品久久久久久久毛片微露脸| 亚洲av五月六月丁香网| 99精品久久久久人妻精品| 亚洲国产精品999在线| 91大片在线观看| 久久久国产精品麻豆| 激情在线观看视频在线高清| 一级作爱视频免费观看| 精品国产一区二区三区四区第35| 亚洲欧美日韩高清在线视频| 黄频高清免费视频| 999久久久精品免费观看国产| aaaaa片日本免费| 80岁老熟妇乱子伦牲交| 美女福利国产在线| 交换朋友夫妻互换小说| 老熟妇仑乱视频hdxx| 少妇的丰满在线观看| 日韩欧美三级三区| 欧美乱码精品一区二区三区| 日韩高清综合在线| 99国产极品粉嫩在线观看| 国产精品永久免费网站| 淫秽高清视频在线观看| 国产乱人伦免费视频| 精品熟女少妇八av免费久了| 国产又色又爽无遮挡免费看| 亚洲 欧美一区二区三区| 亚洲国产精品一区二区三区在线| 国产成人欧美| 婷婷丁香在线五月| 国产一区二区三区综合在线观看| 欧美日韩瑟瑟在线播放| 久久九九热精品免费| 久久国产乱子伦精品免费另类| 伦理电影免费视频| 国产av精品麻豆| 中文亚洲av片在线观看爽| 日本vs欧美在线观看视频| 国产精品久久久人人做人人爽| 男男h啪啪无遮挡| 电影成人av| www.自偷自拍.com| 国产av一区在线观看免费| 亚洲人成77777在线视频| 一级a爱视频在线免费观看| 无人区码免费观看不卡| 日韩欧美一区二区三区在线观看| 色综合婷婷激情| 中国美女看黄片| 日本五十路高清| 国产日韩一区二区三区精品不卡| 91av网站免费观看| 国产成人欧美| 亚洲男人的天堂狠狠| av有码第一页| 一级黄色大片毛片| 性少妇av在线| 亚洲男人天堂网一区| 日韩免费高清中文字幕av| 久久九九热精品免费| 天天影视国产精品| 啦啦啦在线免费观看视频4| 免费人成视频x8x8入口观看| 好男人电影高清在线观看| 黄色丝袜av网址大全| av在线天堂中文字幕 | 啪啪无遮挡十八禁网站| 亚洲av成人不卡在线观看播放网| 搡老乐熟女国产| 欧美日韩亚洲综合一区二区三区_| 精品国产亚洲在线| 自拍欧美九色日韩亚洲蝌蚪91| 久久香蕉国产精品| 亚洲欧美日韩另类电影网站| 十分钟在线观看高清视频www| 亚洲第一青青草原| 一区二区三区国产精品乱码| 性色av乱码一区二区三区2| 欧美精品亚洲一区二区| 久久久久九九精品影院| 好男人电影高清在线观看| 大码成人一级视频| 他把我摸到了高潮在线观看| 免费观看精品视频网站| 久久精品人人爽人人爽视色| 日韩精品中文字幕看吧| 大香蕉久久成人网| 日韩精品青青久久久久久| 成人亚洲精品av一区二区 | 国产99久久九九免费精品| 国内毛片毛片毛片毛片毛片| 国产精品98久久久久久宅男小说| 一区在线观看完整版| 久久人人精品亚洲av| 国内久久婷婷六月综合欲色啪| 久久精品影院6| 无限看片的www在线观看| 又大又爽又粗| 亚洲av成人不卡在线观看播放网| 久久久国产一区二区| www.999成人在线观看| 成人永久免费在线观看视频| 一边摸一边抽搐一进一小说| 999久久久国产精品视频| 成人三级做爰电影| 天堂动漫精品| 亚洲自偷自拍图片 自拍| 亚洲五月婷婷丁香| 女人高潮潮喷娇喘18禁视频| 久久伊人香网站| 国产精品九九99| 国产精品国产高清国产av| 国产熟女午夜一区二区三区| 丁香欧美五月| 9色porny在线观看| 久久人妻福利社区极品人妻图片| 亚洲欧美精品综合一区二区三区| 在线观看日韩欧美| 免费少妇av软件| 淫妇啪啪啪对白视频| 久久中文字幕人妻熟女| 久热爱精品视频在线9| 美女 人体艺术 gogo|