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

    Current sensor based on diamond nitrogen-vacancy color center

    2023-09-05 08:47:56ZiYangShi史子陽WeiGao高偉QiWang王啟HaoGuo郭浩JunTang唐軍ZhongHaoLi李中豪HuanFeiWen溫?zé)wZongMinMa馬宗敏andJunLiu劉俊
    Chinese Physics B 2023年7期
    關(guān)鍵詞:高偉劉俊唐軍

    Zi-Yang Shi(史子陽), Wei Gao(高偉), Qi Wang(王啟), Hao Guo(郭浩), Jun Tang(唐軍),Zhong-Hao Li(李中豪), Huan-Fei Wen(溫?zé)w), Zong-Min Ma(馬宗敏), and Jun Liu(劉俊)

    Shanxi Province Key Laboratory of Quantum Sensing and Precision Measurement,School of Instrument and Electronics,North University of China,Taiyuan 030051,China

    Keywords: current sensor,diamond,high precision,nitrogen-vacancy(NV)color center

    1.Introduction

    As a basic physical quantity, it is extremely significant to accurately measure the current.Physicists and engineers have been exploring various methods to measure current.Precise DC current sensing technology is also the need of many scientific research and industrial applications.[1]For example,current source usually realizes high-precision current control through feedback, so high-precision current measurement is very important for accurate current control.The nuclear magnetic resonance gyroscope realizes the precise measurement of angular velocity through the precise measurement of nuclear spin precession frequency in the magnetic field, which requires a high-precision and stable magnetic field.The highprecision current measurement is of great significance to generate the stable magnetic field in the nuclear magnetic resonance gyroscope; Accurate DC current measurement can be used for the evaluation of pA level current generators to calibrate other instruments.[2]In addition, current sensors also have important applications in other industries, such as automotive power management, solar cell system, smart grid and other application industries that need non-contact accurate current measurement.[3]

    Current sensors have developed to different degrees at home and abroad, which can be roughly divided into direct measurement and indirect measurement, and there are many kinds of non-contact current sensors.Current sensors commonly used in the industrial field based on Ampere loop law include Rogowski coil current sensor, Hall effect current sensor,[4,5]fluxgate current sensor,[6,7]magnetoresistance(AMR, GMR, TMR, CMR, etc) and current transformer,[8]which have their own advantages and disadvantages in all aspects.In other fields,there are current sensors based on NMR,quantum Hall effect, superconducting quantum interference squid and other principles.They have high requirements for the application environment and are expensive.So far, some technologies are not mature and are in the development or improvement stage.A small amount of them are used in laboratory instruments and equipment.In addition, the optical fiber current sensor based on Faraday magneto-optical effect has good performance for measuring AC large current (such as 100 kA),but its performance problem for DC measurement needs to be solved urgently.

    Rogowski coil current sensor and Hall current sensor are relatively traditional current sensors.Their sensitivity is generally low, and their linearity is between 0.05%and 1%.[9,10]New current sensors include optical fiber current sensor based on Faraday magneto-optical effect and GMR current sensor based on giant magnetoresistance effect.They have very high sensitivity, and the linearity of optical fiber type is 0.2%,The linearity of GMR current sensor is between 0.001% and 0.05%.[11–13]It can be seen from the above that for the field of precision detection, the existing current sensors have a common disadvantage that the linearity is not high enough.In addition, optical fiber current sensors are prone to aging.Therefore, a new current sensing method based on diamond nitrogen-vacancy(NV)color center is proposed to solve these problems.It not only has very high linearity and sensitivity, but also has the advantages of low energy consumption and good time stability.At present, there is a large development space for current sensors based on diamond at home and abroad.[14–17]In the field of precision measurement, it has a very broad application prospect.[18,19]

    In this paper,a high-precision DC current sensor based on diamond NV color center is proposed.By using a ring-shaped magnetic concentrating structure with open air gap, the magnetic field generated by the current in the conductor is concentrated to generate a magnetic field that has a stable relationship with the current to be measured.In order to achieve high precision,the diamond NV color center is used to accurately measure the air gap magnetic field,and the corresponding accurate current value is obtained by solving the relevant formula.Finally, we experimentally demonstrate the high accuracy and linearity of the DC current sensor.

    2.Methods and experiment

    Diamond has a strict face centered cubic latticestructure,and there is a carbon atom at any lattice point on the diamond cell.[20,21]NV color center is a stable defect structure, calledC3Vsymmetric structure, formed by capturing a hole nearby after a carbon atom in diamond crystal is replaced by a nitrogen atom.[22]Therefore, there are four possible orientations of NV color centers in diamond.NV color centers in these four directions have equal probability distribution,and the included angle of each NV axis is 109.47?,as shown in Fig.1(a).Its ground state is triplet state (3A2), which isms=0 state,ms=+1 state andms=?1 state respectively.When the external magnetic field is 0,ms=+1 state andms=?1 state are degenerate double states(ms=±1),and there is a natural zero field splitting with a value of 2.87 GHz betweenms=±1 andms=0 state.[23]The zero-field splitting constant is denoted byD.When the applied external magnetic field is not zero,thems=±1 state is separated from each other by energy level splitting.Among them,two metastable singlets1A1and1E exist between the ground state and excited state of diamond NV color center,and their existence greatly affects the fluorescence signal of NV color center.The transition mechanism of NV color center is shown in Fig.1(b).When a 532 nm green laser is used to irradiate the NV color center, the electrons in the ground state of the NV color center will absorb energy and transition to the excited state3E.The NV color center in the excited state3E has two paths back to the ground state.One path is that the electrons in the excited statems=0 will radiate and transition to the ground statems=0 due to instability, releasing red fluorescence signal.The other path does not release fluorescent, as described below.The electrons in the partially excited statems=+1 have a radiative transition back to the ground statems=+1, some have a non radiative transition,decay to the metastable state1A1and1E,and then have a non radiative transition to the ground statems=0.It can be seen that using 532 nm wavelength laser to excite the NV color center can make the NV color center spin polarized to the ground statems=0.[24]The electron spin state of NV color center can be determined by the collected fluorescence signal.

    Fig.1.Schematic diagram of diamond NV color center and energy level structure.(a)The schematics show four possible axes of the NV center in a diamond crystal.(b)The NV center is excited by 532 nm laser and emits photons with a broad spectrum of 637–800 nm.

    By irradiating 532 nm laser and scanning the fluorescence signal of NV color center at a certain frequency, the optical detection magnetic resonance (ODMR) spectrum can be obtained.Under the condition of no magnetic field, only 2.87 GHz formant exists.When an external magnetic field is applied, the degenerate statems=±1 occurs Zeeman splitting,forming two frequency formants.And the magnitude of the external magnetic field projected onto the NV axis can be calculated by the formant frequency of the ODMR spectrum.The ground state Hamiltonian(H)of NV color center is[25]

    whereDis the zero-field splitting constant of NV color center due to electron spin–spin interaction, and its value is 2.87 GHz;Bis the external magnetic induction intensity;Sis the electron spin state of NV color center, and the value is 1;sx,sy, andszare the components of spin angular momentum in Cartesian coordinates;gis the RAND factor;Eis electric field stress;μBis Bohr magneton.The relation between Zeeman splitting degree(or frequency shift)and external magnetic induction intensity is

    Through simplification,we can get

    where ?fis frequency shift;γis the gyromagnetic ratio;his Planck constant;θis the angle between the projection of the NV axis in the magnetic field direction and the NV axis.In this experiment, one of the NV axes of the diamond was adjusted to be consistent with the direction of the magnetic field to be measured,that is,the direction of the NV axis was parallel to the direction of the magnetic field,as shown in Fig.1(a).As a result, the angle between the NV axis and the magnetic field direction is 0 and cosθ=1.

    The magnetic field around the energized wire is not a uniform magnetic field, so it is not suitable for direct measurement.In this paper, a ring-shaped magnetic gathering structure with open air gap(hereinafter referred to as the magnetic gathering ring) is used to gather the magnetic field generated by the current in the wire,[26,27]so that a uniform and parallel magnetic field is formed in the center of the air gap of the magnetic gathering ring.And the magnetic field generated by the electrified wire is in direct proportion to the current flowing.From the accurate measured value of the magnetic field, the corresponding accurate value of the current can be obtained.Therefore, combined with the relationship between the magnetic field calculated by ODMR formant frequency shift and the magnetic field and current in the air gap of the magnetic concentrating ring, the current passing through the conductor can be obtained.

    Fig.2.(a) Structure diagram of current sensor.The dotted box shows the partial top view of the magnetic gathering ring and diamond.(b)Simulation diagram of air gap magnetic field of magnetic gathering ring.The results show that a uniform magnetic field can be formed in the center of the air gap.

    The sensing system for diamond current measurement is shown in Fig.2(a).The corresponding physical image is shown in Fig.S1 in the supplementary materials.In the experiment, first, the green laser is emitted through the MGL-III-532-100 mW laser, then all the green light is reflected through the long wave pass dichroic mirror (Thorlabs DMLP605R),and through the objective lens(OLYMPUS PLN 10×NA=0.25)focus on the diamond,and the NV color center of the diamond will excite a red fluorescent signal after laser irradiation.The fluorescent signal is filtered by the long wave pass filter(Thorlabs FELH0600)and the long wave pass dichroic mirror,and then focused by the short wave pass filter(Thorlabs FESH0800)and the condensing lens(ACL25416UB NA=0.79),and finally collected by the photodetector(Thorlabs PDA100A2).In addition, current source (KEYSIGHT B2962A) is used to provide stable high-precision DC current of different sizes.The swept microwave is transmitted to the diamond by the microwave source (KEYSIGHT N5183B MXG)connected to the microstrip antenna.The output power is?8 dBm and the frequency range is 2.51–3.24 GHz.In this way, by changing the input current through the current source,the change of the magnetic field in the center of the air gap of the magnetic concentrating ring under different current conditions can be measured,which changes the resonance frequency of the ODMR signal, so as to realize the detection of the current by the diamond NV color center.

    The diamond used in the experiment was purchased from ElementSix with a size of 1.5 mm×1.5 mm×1 mm, nitrogen content is about 10?4.Polished surface of diamond(1.5 mm×1.5 mm) is the (100) surface, which is irradiated for 3 h under 10 MeV electron irradiation and then annealed for 2 h under 850?C.[28–33]The NV color center concentration formed is about 10?6.The functions of magnetic gathering ring mainly include magnetic gathering and antiinterference, so it is appropriate to use soft magnetic material and it is required to have high permeability, low coercivity low loss, high resistivity and high saturation magnetic induction.The material properties of permalloy fully meet the above requirements, and permalloy products have been standardized, low material cost.Therefore, after comprehensive consideration,1J85 permalloy is used in this experiment,which has a high magnetic permeability (initial permeability of 30000 mH/m, maximum permeability of 150000 mH/m),lower than 1.6 Hc/A·m?1coercivity.And 0.8 T saturation magnetic induction intensity(much higher than mT level magnetic field intensity generated during the experiment), these excellent performance can satisfy completely the requirements of the experiment.[34]And its structural dimensions are 26 mm inner diameter, 40 mm outer diameter, 7 mm height and 3.5 mm air gap width.COMSOL Multiphysics software is used to simulate the structure of the magnetic gathering ring,as shown in Fig.2(b).The darker the color in the figure, the greater the magnetic field.The simulation results show that the magnetic gathering ring has a good ability to gather the magnetic field around the energized wire,and a uniform magnetic field can be formed in the center of the air gap.Besides,a 4 mm2copper core wire was used in the experiment.The maximum allowable current is 32 A,which exceeds the maximum current 8 A generated by the required current source in the experiment, which can ensure the safety and measurement requirements of the experiment.In addition,there is no current source with large range output that meets the measurement range under laboratory conditions,so the equal Ampere turn method is used for the experiment.A 10 m long wire is wound into a pattern of 5 turns side by side and passes through the center of the magnetic gathering ring,that is,the 5 turns of wire connected side by side into a certain currentIis equivalent to a wire connected with a 5Icurrent,so as to expand the measurement range of the designed sensor.

    3.Results and discussion

    Once the sensing system is built, the magnetic field of the air gap in the magnetic gathering ring is in a stable proportional relationship with the input current in the wire.In the experiment,firstly,the change of the magnetic field of the air gap of the magnetic gathering ring with the current under the input current of the current source of 1–8 A (the equivalent current through the equal Ampere turn method is 5–40 A)is measured by a Gauss meter, and the average value of 10 groups of measurement results is taken (Table S1 in the supplementary materials).The linear relationship between the air gap magnetic field and the input current is obtained by linear fitting of the data.The available formula is

    whereB0is the air gap magnetic field,I0is the input current,andkandbare fitting constants.

    The results show thatk=0.2419 mT/A,b=0.08571 mT,so the linear relationship between the calculated current and the magnetic field calculated by the diamond ODMR resonance frequency solution can be deduced from formula (5),which can be expressed as

    whereI1is the output current andB1is the air gap magnetic field calculated by ODMR resonance frequency solution.

    By changing the input current and observing and recording the formant frequency shift of diamond NV color center ODMR signal, ODMR curves under different current conditions can be obtained, as shown in Fig.3(a).There are two pairs of formants in the figure,of which the outer pair is generated by the NV axis parallel to the direction of the magnetic field, and the middle pair is generated by the other three NV axes together and superimposed together.At the same time,we can also get the relationship between the ODMR signal frequency shift ?fand the input currentI0,which is expressed as

    And it is easy to know that the corresponding coefficient fitting value in the formula is thatk′=13.50475 MHz/A,b′=4.88361 MHz.It can be seen that with the increase of the input current,the distance between the two resonance peak points on the outside of the ODMR curve increases,that is,with the increase of the current,the magnetic field strength formed at the center of the air gap where the magnetic concentrating ring is placed on the diamond gradually increases.In addition,as the fluorescence intensity contrast of ODMR signal will decrease with the increase of magnetic field intensity,the current range measured in the experiment is 5–40 A after considering the comprehensive influence of magnetic field on the experimental results.

    Fig.3.(a)Shift of ODMR formant frequency under different input currents.(b)The relationship between the calculated current and the actual current.

    In order to evaluate the performance of the designed current sensor, we use Keysight B2962A current source, which can be used as a standard current to calibrate the output current calculated by the current sensor.We set several test points from 5 A to 40 A,and each point is the average of the data recorded ten times.In addition, we have simulated and collected 100 sets of data for the regional magnetic field of 2 mm×2 mm in the middle of the air gap,and the percentage error of the magnetic field strength data is less than 0.005%,so the position error of the diamond can be ignored.After placing the diamond in the air gap,the positions of all devices,including the wire, the diamond, the magnetic gathering ring and so on, will be fixed.And from the experimental data, it can be seen that the experimental repeatability is good,so the reading of the current sensor is stable.

    The output current is calculated from the ODMR resonance frequency solution,and the linearity of input current and output current is analyzed and fitted.The results in Fig.3(b).show the relationship between the calculated current and the actual current.It can be seen that the current sensor shows a high linearity in a large working range of 5–40 A, about 33 ppm, which is much higher than the current sensors commonly used in the market.The formula for calculating the linearity is shown in the supplementary materials as Eq.(S1).And the performance of the proposed current sensor compared with other technologies is shown in Table S2 of the supplementary materials.The new current sensor designed in this paper can achieve such excellent linearity thanks to the important role of the magnetic gathering ring and diamond NV color center.On the one hand,the ring enhances the magnetic induction intensity generated by the measured current and improves the measurement sensitivity.It also eliminates the interference of magnetic fields in the environment.On the other hand,thanks to the excellent measurement performance of diamond NV color center ensemble system,quantum operation and readout can be carried out through the spin of NV color center,which can realize quantitative and nondestructive measurement of magnetic properties.Therefore,the current sensor has an approximate and low absolute error in the whole working range.

    When the measuring current range is 5–40 A,the absolute error of the measuring current is less than 51μA,as shown in Fig.4(a).Equation(S2)shown in the supplementary materials is the formula for calculating the absolute error.In the current sensor,the equivalent accuracy of magnetic field measurement is better than 0.56μT.The absolute error of the current source used in the experiment and the systematic error caused by the wire placement will help to improve the performance of the designed current sensor prototype.According to experience,in this experiment, the wire is just completely placed in the inner hole of the sensor without any gap.The wire is placed in the center of the magnetic gathering ring through hole of the sensor to ensure that the position is not biased as much as possible, and the maximum absolute error of the current source used is very small and can be ignored.

    Since the absolute error has no obvious change,the relative error decreases with the increase of the actual current,as shown in Fig.4(b).The relative error at 40 A is 2.42×10?6.Equation(S3)shown in the supplementary materials is the formula for calculating the absolute error.In addition, we also evaluated the performance of the current sensor through Allan variance, as shown in Fig.5.When the current is 10 A, the sampling rate is 250 Hz,and the output of the current sensor is recorded for 1 hour.The results show that the short-term stability of the sensor is lower than 0.0002 A.The value related to the average time of 10 min is consistent with the absolute error data in the preceding paragraph.

    Fig.4.(a) Absolute error and under actual current.(b) Relative error under actual current.

    Fig.5.Allan variance of the current sensor data.The short-term stability of the sensor is lower than 0.0002 A.

    Although the preliminary experiment in this paper has achieved quite good results, there are still many details to be improved, and many other factors need to be considered and tested in future work.For example,temperature is an influential factor that almost any experiment should consider.According to the research of Shaozhuo Lin and others,[35]a modulation and demodulation technology and PID frequency locking technology are used to realize independent detection of two spin symmetry states of NV color center,the common mode cancellation method is used to suppress spin information noise caused by temperature, so as to suppress temperature noise in NV magnetometer and improve accuracy and sensitivity of NV magnetometer.At present,the dual formant system can suppress temperature noise by 10 times and this method still has a lot of room for improvement.In future experiments,this method can be combined with the current sensor based on diamond NV color center designed in this paper,so that the NV quantum sensor can be applied to a wider environment.In addition,the noise in the current should also be considered and filtered to reduce the impact on the experimental results in the future optimization experiments.The traditional method of installing harmonic compensation devices to suppress harmonics and noise can be considered,that is,using LC tuned filters.This method has a simple structure, which can compensate both harmonics and reactive power, and can meet the experimental requirements.Furthermore,we can also optimize the future experimental results from the design of the experimental device itself.For example,in order to further optimize the experimental system, we can make design choices from the shape,material,thickness,air gap size and other aspects of the magnetic focusing structure,so as to obtain better magnetic focusing effect and better experimental results.

    4.Conclusion

    In summary, we have designed and proved a highprecision DC current sensor based on diamond NV color center,which can maintain an ultra-high linearity in a large measurement range, as high as 33 ppm.When the range of input current is 5–40 A, the absolute error of calculated current is less than 51μA,and the relative error is 2.42×10?6at 40 A.In addition,there are other interference effects during the experiment,such as temperature,wire placement,noise of electric equipment, etc., which will have a certain impact on the experimental data.Even so, such a high measurement linearity is enough to show its outstanding performance.Compared with other current sensors,the current sensor based on optical diamond NV color center has outstanding advantages,and will have a very wide application prospect in the field of precision measurement.

    Acknowledgements

    Project supported in part by the National Natural Science Foundation of China (Grant Nos.51922009, 51727808,62175219, 62103385, and 51821003), the Key Laboratory of Shanxi Province(Grant No.201905D121001),and the Shanxi‘1331 Project’Key Subjects Construction.

    猜你喜歡
    高偉劉俊唐軍
    Effect of a static pedestrian as an exit obstacle on evacuation
    High-fidelity resonant tunneling passage in three-waveguide system
    劉俊
    唐軍治療圍絕經(jīng)期失眠經(jīng)驗總結(jié)
    Unpinning the spiral waves by using parameter waves*
    Analysis of asymmetry of the Dα emission spectra under the Zeeman effect in boundary region for D–D experiment on EAST tokamak?
    秋天的風(fēng)
    金山(2020年3期)2020-04-15 03:56:36
    In fluence of Ni/Mn ratio on magnetostructural transformation and magnetocaloric effect in Ni48?x Co2Mn38+x Sn12(x=0,1.0,1.5,2.0, and 2.5)ferromagnetic shape memory alloys?
    我和你打個賭
    唐軍 留守少年的逆襲
    時代人物(2014年12期)2015-01-29 14:10:30
    水蜜桃什么品种好| 中文字幕av成人在线电影| 日韩欧美一区视频在线观看 | 内地一区二区视频在线| 欧美成人午夜免费资源| 大又大粗又爽又黄少妇毛片口| 亚洲av成人精品一区久久| 久久久欧美国产精品| .国产精品久久| 99精国产麻豆久久婷婷| 我的老师免费观看完整版| 偷拍熟女少妇极品色| 少妇丰满av| 欧美国产精品一级二级三级 | 国产免费视频播放在线视频| 99久久综合免费| 激情五月婷婷亚洲| 亚洲美女视频黄频| 久久热精品热| 色综合色国产| 久久精品熟女亚洲av麻豆精品| 亚洲一级一片aⅴ在线观看| 亚洲av二区三区四区| 中国美白少妇内射xxxbb| 小蜜桃在线观看免费完整版高清| av线在线观看网站| 一级av片app| 少妇 在线观看| 丝袜脚勾引网站| 日本欧美国产在线视频| 性色avwww在线观看| 国产免费视频播放在线视频| 亚洲av男天堂| 亚洲精品视频女| 精品久久久噜噜| 97在线人人人人妻| 纯流量卡能插随身wifi吗| 国产片特级美女逼逼视频| 毛片一级片免费看久久久久| 精品少妇久久久久久888优播| 男女国产视频网站| 国产成人精品婷婷| 黑人高潮一二区| 丰满人妻一区二区三区视频av| 日本vs欧美在线观看视频 | 久久99蜜桃精品久久| 又大又黄又爽视频免费| 草草在线视频免费看| 国产永久视频网站| 不卡视频在线观看欧美| 国产成人精品一,二区| 国产精品国产三级国产专区5o| 欧美xxxx黑人xx丫x性爽| 91久久精品国产一区二区三区| 少妇丰满av| 男女啪啪激烈高潮av片| 日韩在线高清观看一区二区三区| 欧美高清性xxxxhd video| 精品久久久精品久久久| 九色成人免费人妻av| 我的女老师完整版在线观看| 精品人妻视频免费看| 在线观看一区二区三区| 久久久亚洲精品成人影院| 亚洲久久久国产精品| 秋霞伦理黄片| 日韩一区二区视频免费看| 免费少妇av软件| 国产一区有黄有色的免费视频| 欧美日韩亚洲高清精品| 亚洲四区av| 亚洲精品国产色婷婷电影| 久久99蜜桃精品久久| 成年美女黄网站色视频大全免费 | 亚洲国产精品专区欧美| 亚洲av中文av极速乱| 亚洲色图av天堂| 国产一区有黄有色的免费视频| 久久精品国产a三级三级三级| 国产精品蜜桃在线观看| 如何舔出高潮| 亚洲成人一二三区av| 日韩一区二区三区影片| 日本黄色片子视频| 欧美丝袜亚洲另类| 性色av一级| 欧美丝袜亚洲另类| 天堂俺去俺来也www色官网| 国产精品蜜桃在线观看| 极品教师在线视频| 国产精品国产av在线观看| 尾随美女入室| 亚洲无线观看免费| 成人免费观看视频高清| 日本免费在线观看一区| 久久久午夜欧美精品| 亚洲精品日韩av片在线观看| 日韩中字成人| 欧美日韩在线观看h| 国产高潮美女av| 国产高潮美女av| 在线观看av片永久免费下载| 在线观看av片永久免费下载| 天天躁夜夜躁狠狠久久av| 欧美精品人与动牲交sv欧美| 91午夜精品亚洲一区二区三区| 国产女主播在线喷水免费视频网站| 26uuu在线亚洲综合色| 伦理电影免费视频| 国产在视频线精品| 亚洲精品乱码久久久v下载方式| 久久久久久久久久久丰满| 国产欧美日韩精品一区二区| 人妻制服诱惑在线中文字幕| 国产精品蜜桃在线观看| 成人午夜精彩视频在线观看| 日本-黄色视频高清免费观看| 黄片无遮挡物在线观看| 欧美精品人与动牲交sv欧美| 亚洲一区二区三区欧美精品| 亚洲美女视频黄频| 免费观看性生交大片5| 欧美3d第一页| 色5月婷婷丁香| 中文欧美无线码| 欧美精品国产亚洲| 久久 成人 亚洲| a级毛色黄片| 国产日韩欧美亚洲二区| 久久这里有精品视频免费| 不卡视频在线观看欧美| 久久久久久久久久成人| 午夜福利网站1000一区二区三区| 国产伦精品一区二区三区视频9| 伦理电影免费视频| 下体分泌物呈黄色| 内射极品少妇av片p| 免费播放大片免费观看视频在线观看| 免费观看的影片在线观看| 一级毛片aaaaaa免费看小| freevideosex欧美| 国模一区二区三区四区视频| 人妻 亚洲 视频| 国产精品福利在线免费观看| 国产亚洲精品久久久com| 国产在线男女| 亚洲欧洲国产日韩| 亚洲人与动物交配视频| 国产成人a区在线观看| 在线观看三级黄色| 简卡轻食公司| 亚洲色图av天堂| 国产精品人妻久久久影院| 最近的中文字幕免费完整| 欧美+日韩+精品| 亚洲国产av新网站| a级毛色黄片| 亚洲av二区三区四区| 又粗又硬又长又爽又黄的视频| 秋霞在线观看毛片| 国产av码专区亚洲av| 免费av中文字幕在线| 如何舔出高潮| 三级国产精品欧美在线观看| 成人特级av手机在线观看| 精品人妻熟女av久视频| 97精品久久久久久久久久精品| 日韩亚洲欧美综合| 国产乱来视频区| 国产乱人偷精品视频| 80岁老熟妇乱子伦牲交| 久久女婷五月综合色啪小说| 亚洲中文av在线| 亚洲精品乱码久久久久久按摩| 中文精品一卡2卡3卡4更新| 性色avwww在线观看| 熟妇人妻不卡中文字幕| 国产免费福利视频在线观看| 亚洲精品久久午夜乱码| 日本欧美视频一区| 国产一区二区在线观看日韩| 狂野欧美白嫩少妇大欣赏| 看免费成人av毛片| 亚洲国产欧美人成| 高清毛片免费看| 国内少妇人妻偷人精品xxx网站| 国精品久久久久久国模美| 国产 一区 欧美 日韩| 天天躁夜夜躁狠狠久久av| 亚洲色图综合在线观看| 国产成人午夜福利电影在线观看| 搡女人真爽免费视频火全软件| 简卡轻食公司| 99精国产麻豆久久婷婷| 日本爱情动作片www.在线观看| 老师上课跳d突然被开到最大视频| 亚洲精品自拍成人| 丰满少妇做爰视频| 成人高潮视频无遮挡免费网站| 国产精品人妻久久久影院| 国产无遮挡羞羞视频在线观看| 久久久久久久久久人人人人人人| 97在线视频观看| 免费黄网站久久成人精品| 两个人的视频大全免费| 丰满少妇做爰视频| 看免费成人av毛片| 伦理电影免费视频| 精品一区二区三卡| 亚洲综合色惰| 午夜精品国产一区二区电影| 久久久午夜欧美精品| 欧美xxxx黑人xx丫x性爽| 国产成人午夜福利电影在线观看| 777米奇影视久久| 18禁动态无遮挡网站| 中文天堂在线官网| 日本一二三区视频观看| 亚洲av二区三区四区| 91久久精品国产一区二区三区| 亚洲精品一区蜜桃| 国产一区二区三区av在线| 日本一二三区视频观看| 成人亚洲精品一区在线观看 | 亚洲精品日本国产第一区| 内射极品少妇av片p| 十分钟在线观看高清视频www | 97在线人人人人妻| 黄片wwwwww| 精品人妻熟女av久视频| 欧美日韩一区二区视频在线观看视频在线| 午夜免费鲁丝| av视频免费观看在线观看| 少妇人妻久久综合中文| 国产精品久久久久久精品古装| 美女脱内裤让男人舔精品视频| 少妇的逼好多水| 国产深夜福利视频在线观看| 国产精品国产三级专区第一集| videossex国产| 色5月婷婷丁香| 自拍欧美九色日韩亚洲蝌蚪91 | .国产精品久久| 在线观看国产h片| 丝袜喷水一区| 久久精品国产自在天天线| 精品人妻熟女av久视频| 午夜福利在线观看免费完整高清在| 中国国产av一级| 乱码一卡2卡4卡精品| 深夜a级毛片| 日韩免费高清中文字幕av| 久久av网站| 久久人人爽人人爽人人片va| av在线老鸭窝| 久久99精品国语久久久| 一本一本综合久久| 久久久久久久久久久免费av| 视频中文字幕在线观看| 蜜桃在线观看..| 成人国产av品久久久| 99久久精品国产国产毛片| 成人高潮视频无遮挡免费网站| 各种免费的搞黄视频| 3wmmmm亚洲av在线观看| 热99国产精品久久久久久7| 99re6热这里在线精品视频| 精品久久国产蜜桃| 国产高清三级在线| www.色视频.com| 韩国高清视频一区二区三区| 日韩av不卡免费在线播放| 18禁在线无遮挡免费观看视频| 亚洲无线观看免费| 中国国产av一级| 国产淫语在线视频| 久久人人爽人人爽人人片va| 国产日韩欧美在线精品| 少妇人妻久久综合中文| 丝瓜视频免费看黄片| 国产黄频视频在线观看| 国产成人aa在线观看| 国产熟女欧美一区二区| 欧美日韩国产mv在线观看视频 | 18禁裸乳无遮挡免费网站照片| 人妻制服诱惑在线中文字幕| 亚洲精品久久午夜乱码| 99国产精品免费福利视频| 激情五月婷婷亚洲| 联通29元200g的流量卡| 国产精品人妻久久久影院| 欧美成人午夜免费资源| 国内揄拍国产精品人妻在线| 纯流量卡能插随身wifi吗| 亚洲真实伦在线观看| 大香蕉97超碰在线| 国产高清不卡午夜福利| 国产免费视频播放在线视频| 国产毛片在线视频| 欧美另类一区| 另类亚洲欧美激情| 夫妻午夜视频| 80岁老熟妇乱子伦牲交| 欧美少妇被猛烈插入视频| 国产日韩欧美在线精品| 免费不卡的大黄色大毛片视频在线观看| 美女视频免费永久观看网站| 日本av免费视频播放| 欧美xxxx性猛交bbbb| 久久人妻熟女aⅴ| 嫩草影院入口| av福利片在线观看| 日韩欧美精品免费久久| 九九在线视频观看精品| 久久韩国三级中文字幕| 五月开心婷婷网| 亚洲av电影在线观看一区二区三区| 欧美成人午夜免费资源| 久久女婷五月综合色啪小说| 男女下面进入的视频免费午夜| 亚洲综合精品二区| 国产精品人妻久久久影院| 国产视频内射| www.av在线官网国产| 97超视频在线观看视频| 亚洲精品成人av观看孕妇| 性高湖久久久久久久久免费观看| av又黄又爽大尺度在线免费看| 少妇人妻精品综合一区二区| 成人亚洲精品一区在线观看 | 黄色视频在线播放观看不卡| 亚洲最大成人中文| 一本—道久久a久久精品蜜桃钙片| 久久99热6这里只有精品| 在线天堂最新版资源| 成人18禁高潮啪啪吃奶动态图 | 国产精品一及| 久久久久久久久久久丰满| 男女下面进入的视频免费午夜| 欧美成人精品欧美一级黄| 亚洲真实伦在线观看| 人妻少妇偷人精品九色| 国产成人精品久久久久久| 中文字幕精品免费在线观看视频 | 最近的中文字幕免费完整| 日韩欧美 国产精品| 久久久色成人| 人体艺术视频欧美日本| 一区二区三区乱码不卡18| av天堂中文字幕网| 精品久久久精品久久久| 亚洲经典国产精华液单| 最近中文字幕高清免费大全6| 亚洲色图综合在线观看| 在线观看免费视频网站a站| 国产精品国产三级国产av玫瑰| 又爽又黄a免费视频| 少妇人妻久久综合中文| 老女人水多毛片| 亚洲av男天堂| videos熟女内射| 高清日韩中文字幕在线| 伊人久久国产一区二区| 日本vs欧美在线观看视频 | 色网站视频免费| 丰满少妇做爰视频| 日本黄大片高清| 久久久色成人| 久久精品人妻少妇| 日本黄色片子视频| 国产一区二区在线观看日韩| 亚洲国产精品国产精品| 七月丁香在线播放| 精品久久久噜噜| 又粗又硬又长又爽又黄的视频| 久久久a久久爽久久v久久| 日日撸夜夜添| 日韩视频在线欧美| 夫妻午夜视频| 欧美日本视频| 在线观看美女被高潮喷水网站| 在线播放无遮挡| 久久ye,这里只有精品| 人妻 亚洲 视频| 亚洲精品国产成人久久av| 国产又色又爽无遮挡免| 午夜免费鲁丝| 国产精品.久久久| 好男人视频免费观看在线| 精品国产一区二区三区久久久樱花 | 国产成人免费无遮挡视频| 伦理电影大哥的女人| 久久6这里有精品| 黄色一级大片看看| 最后的刺客免费高清国语| 大陆偷拍与自拍| 国语对白做爰xxxⅹ性视频网站| 黑丝袜美女国产一区| 欧美xxxx黑人xx丫x性爽| 熟女电影av网| a 毛片基地| 丰满乱子伦码专区| 国产一区有黄有色的免费视频| 777米奇影视久久| 国产欧美日韩精品一区二区| 丰满迷人的少妇在线观看| 内射极品少妇av片p| 亚洲欧美日韩卡通动漫| 波野结衣二区三区在线| 国产精品久久久久成人av| 国产免费一级a男人的天堂| 黄色视频在线播放观看不卡| 国产精品福利在线免费观看| 成年av动漫网址| 国产成人免费观看mmmm| 国产免费又黄又爽又色| 激情五月婷婷亚洲| 国产在线视频一区二区| 亚洲成人手机| 老司机影院毛片| 亚洲va在线va天堂va国产| 午夜福利在线观看免费完整高清在| 麻豆成人午夜福利视频| 亚洲真实伦在线观看| 亚洲精品456在线播放app| 少妇被粗大猛烈的视频| 免费看日本二区| 亚洲精品亚洲一区二区| 免费黄色在线免费观看| 亚洲成人中文字幕在线播放| 日本黄色片子视频| 亚洲精品乱码久久久久久按摩| 乱系列少妇在线播放| 国产成人freesex在线| 亚洲国产欧美人成| 在线精品无人区一区二区三 | 日本黄大片高清| 最近的中文字幕免费完整| 国产午夜精品久久久久久一区二区三区| 18禁动态无遮挡网站| 日韩制服骚丝袜av| 精品视频人人做人人爽| 成人18禁高潮啪啪吃奶动态图 | 国产久久久一区二区三区| 丰满人妻一区二区三区视频av| 成人一区二区视频在线观看| 国产精品女同一区二区软件| .国产精品久久| 老司机影院成人| 日日啪夜夜撸| 中文字幕免费在线视频6| 一级毛片 在线播放| 联通29元200g的流量卡| 在线观看美女被高潮喷水网站| 三级经典国产精品| 男人爽女人下面视频在线观看| 亚洲激情五月婷婷啪啪| 亚洲精品456在线播放app| 亚洲国产精品专区欧美| 国产伦精品一区二区三区四那| 久久综合国产亚洲精品| 黑丝袜美女国产一区| 亚洲国产精品成人久久小说| 国产真实伦视频高清在线观看| 国产精品熟女久久久久浪| 国产精品一区二区性色av| 在线免费观看不下载黄p国产| 街头女战士在线观看网站| 男人狂女人下面高潮的视频| 又大又黄又爽视频免费| 男人添女人高潮全过程视频| 国产综合精华液| 老司机影院毛片| 简卡轻食公司| 五月伊人婷婷丁香| 精品亚洲成国产av| 噜噜噜噜噜久久久久久91| 国产精品一区二区三区四区免费观看| 一本一本综合久久| 街头女战士在线观看网站| 久久精品国产a三级三级三级| 一级片'在线观看视频| 男女无遮挡免费网站观看| 啦啦啦在线观看免费高清www| 国产高清三级在线| 亚洲美女视频黄频| 久久精品久久久久久久性| 91久久精品电影网| 日本一二三区视频观看| 国产av码专区亚洲av| 欧美成人一区二区免费高清观看| 性高湖久久久久久久久免费观看| av国产精品久久久久影院| 亚洲久久久国产精品| 国产综合精华液| 高清视频免费观看一区二区| 丰满乱子伦码专区| 国产伦精品一区二区三区四那| 国产国拍精品亚洲av在线观看| 少妇精品久久久久久久| 欧美日韩视频精品一区| 亚洲av综合色区一区| 久久久午夜欧美精品| 女性被躁到高潮视频| 成年美女黄网站色视频大全免费 | 欧美日韩精品成人综合77777| 嘟嘟电影网在线观看| 国产真实伦视频高清在线观看| 久久久久久人妻| 丰满人妻一区二区三区视频av| av在线老鸭窝| 亚洲国产欧美人成| 成人免费观看视频高清| 水蜜桃什么品种好| 少妇人妻 视频| 久久精品熟女亚洲av麻豆精品| 久久99热这里只频精品6学生| 亚洲电影在线观看av| 水蜜桃什么品种好| 一区二区av电影网| 久久久欧美国产精品| 六月丁香七月| 男女边摸边吃奶| 国产欧美日韩一区二区三区在线 | www.色视频.com| 国产精品国产av在线观看| 六月丁香七月| 这个男人来自地球电影免费观看 | 精品99又大又爽又粗少妇毛片| 欧美精品国产亚洲| 免费大片黄手机在线观看| 在线 av 中文字幕| 精品国产乱码久久久久久小说| 精品视频人人做人人爽| 91久久精品国产一区二区成人| 成人毛片a级毛片在线播放| 国产欧美日韩一区二区三区在线 | 免费不卡的大黄色大毛片视频在线观看| 尾随美女入室| 毛片女人毛片| av国产免费在线观看| 男人添女人高潮全过程视频| 色视频在线一区二区三区| 亚洲av二区三区四区| 少妇人妻精品综合一区二区| 久久久久久伊人网av| 日韩成人伦理影院| 毛片女人毛片| 最黄视频免费看| 中文字幕久久专区| 高清不卡的av网站| 国产永久视频网站| 伊人久久国产一区二区| 偷拍熟女少妇极品色| 亚洲欧美一区二区三区黑人 | 久久亚洲国产成人精品v| 亚洲精品国产av成人精品| 亚洲欧美日韩无卡精品| 欧美97在线视频| 一二三四中文在线观看免费高清| 中文天堂在线官网| 直男gayav资源| 成人漫画全彩无遮挡| 日本免费在线观看一区| 久热久热在线精品观看| 午夜激情久久久久久久| 免费大片黄手机在线观看| 日韩 亚洲 欧美在线| 精品久久久久久电影网| 蜜桃在线观看..| 中文字幕精品免费在线观看视频 | 精品一区二区三卡| 美女福利国产在线 | 亚洲国产色片| 日韩不卡一区二区三区视频在线| 欧美高清成人免费视频www| 婷婷色综合www| 亚洲成人中文字幕在线播放| 99久久精品一区二区三区| 亚洲成人手机| 街头女战士在线观看网站| 亚洲三级黄色毛片| 99久久综合免费| 国产精品成人在线| 国产淫语在线视频| 久久午夜福利片| av在线观看视频网站免费| 久久久a久久爽久久v久久| av国产精品久久久久影院| 男女边摸边吃奶| 国产高清三级在线| 一二三四中文在线观看免费高清| www.色视频.com| 麻豆成人av视频| 老司机影院成人| 亚洲怡红院男人天堂| 国产成人一区二区在线| 国产精品久久久久久精品电影小说 | 国产精品99久久99久久久不卡 | 精品国产三级普通话版| 久久亚洲国产成人精品v| 五月天丁香电影| 久久久久精品久久久久真实原创| 日韩在线高清观看一区二区三区| 国产精品熟女久久久久浪| 亚洲av成人精品一二三区| 亚洲性久久影院| 乱系列少妇在线播放| 成人国产av品久久久| 免费少妇av软件| 91精品伊人久久大香线蕉| 午夜福利在线在线| 夜夜爽夜夜爽视频| 男男h啪啪无遮挡| 亚洲美女黄色视频免费看|