• <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久久热| 少妇猛男粗大的猛烈进出视频| 韩国高清视频一区二区三区| 国产91精品成人一区二区三区 | 欧美日本中文国产一区发布| 久久九九热精品免费| 91麻豆av在线| 老司机影院成人| 国产深夜福利视频在线观看| 天天添夜夜摸| 欧美黄色淫秽网站| 热99国产精品久久久久久7| 精品亚洲乱码少妇综合久久| 人人妻人人澡人人爽人人夜夜| 国产色视频综合| 夫妻午夜视频| 欧美精品啪啪一区二区三区 | 日日夜夜操网爽| 欧美激情高清一区二区三区| 欧美日韩亚洲综合一区二区三区_| 视频区欧美日本亚洲| 久久女婷五月综合色啪小说| 看十八女毛片水多多多| 国产亚洲精品久久久久5区| 日韩人妻精品一区2区三区| av福利片在线| 成人国产av品久久久| 丰满饥渴人妻一区二区三| 日韩制服丝袜自拍偷拍| 免费少妇av软件| 波多野结衣av一区二区av| 精品一区在线观看国产| 日韩制服骚丝袜av| 2018国产大陆天天弄谢| 99热网站在线观看| 2018国产大陆天天弄谢| 国产高清videossex| 国产亚洲欧美精品永久| 欧美少妇被猛烈插入视频| 人人澡人人妻人| 18禁黄网站禁片午夜丰满| 极品人妻少妇av视频| 男的添女的下面高潮视频| 1024视频免费在线观看| 欧美日韩视频精品一区| 国产精品秋霞免费鲁丝片| 一区二区三区四区激情视频| 亚洲av在线观看美女高潮| 国产三级黄色录像| 七月丁香在线播放| 亚洲午夜精品一区,二区,三区| 国产又爽黄色视频| 亚洲激情五月婷婷啪啪| 校园人妻丝袜中文字幕| 亚洲av在线观看美女高潮| 色94色欧美一区二区| 久久这里只有精品19| 精品高清国产在线一区| 精品亚洲成a人片在线观看| 亚洲人成77777在线视频| 日本一区二区免费在线视频| 亚洲国产最新在线播放| 亚洲国产毛片av蜜桃av| 男的添女的下面高潮视频| 国产成人啪精品午夜网站| 欧美日韩精品网址| 欧美av亚洲av综合av国产av| 一本一本久久a久久精品综合妖精| 国产精品偷伦视频观看了| 精品一区在线观看国产| 免费在线观看影片大全网站 | 99香蕉大伊视频| 亚洲色图综合在线观看| 新久久久久国产一级毛片| 黑丝袜美女国产一区| 极品人妻少妇av视频| 人妻人人澡人人爽人人| 精品视频人人做人人爽| 久久精品亚洲av国产电影网| av国产精品久久久久影院| 少妇粗大呻吟视频| 黑丝袜美女国产一区| 日韩大片免费观看网站| 狂野欧美激情性bbbbbb| 91精品三级在线观看| 亚洲色图 男人天堂 中文字幕| 50天的宝宝边吃奶边哭怎么回事| 中文字幕制服av| 国产又爽黄色视频| 人妻一区二区av| 视频区欧美日本亚洲| 欧美精品人与动牲交sv欧美| 最新的欧美精品一区二区| 制服诱惑二区| 亚洲五月婷婷丁香| 午夜老司机福利片| 一级毛片电影观看| 久久青草综合色| 国产老妇伦熟女老妇高清| 秋霞在线观看毛片| 亚洲精品一卡2卡三卡4卡5卡 | 亚洲精品av麻豆狂野| 久久精品亚洲熟妇少妇任你| av在线app专区| 男女床上黄色一级片免费看| 黄色一级大片看看| 久久久国产欧美日韩av| 午夜影院在线不卡| 五月开心婷婷网| 日本a在线网址| 国产一区二区三区综合在线观看| 人成视频在线观看免费观看| 欧美日韩成人在线一区二区| 久久精品成人免费网站| 麻豆国产av国片精品| 日韩一卡2卡3卡4卡2021年| 久久综合国产亚洲精品| 美女国产高潮福利片在线看| 亚洲精品一区蜜桃| 久久久欧美国产精品| 国产亚洲精品久久久久5区| 国产欧美日韩综合在线一区二区| 欧美国产精品一级二级三级| 精品久久久精品久久久| 亚洲成人免费电影在线观看 | 久久性视频一级片| 久久久国产一区二区| 天天躁日日躁夜夜躁夜夜| 性色av一级| 久久人人97超碰香蕉20202| 久热这里只有精品99| 国产免费现黄频在线看| 国产精品偷伦视频观看了| 亚洲欧洲精品一区二区精品久久久| 又黄又粗又硬又大视频| 欧美国产精品va在线观看不卡| 一级黄色大片毛片| 精品人妻一区二区三区麻豆| a级片在线免费高清观看视频| 天堂中文最新版在线下载| 69精品国产乱码久久久| 青春草亚洲视频在线观看| 人人妻人人添人人爽欧美一区卜| 国产老妇伦熟女老妇高清| 黄色 视频免费看| 亚洲图色成人| 宅男免费午夜| bbb黄色大片| 久久精品亚洲熟妇少妇任你| 在线天堂中文资源库| 免费不卡黄色视频| 日韩一区二区三区影片| 99精国产麻豆久久婷婷| 搡老乐熟女国产| 亚洲国产欧美日韩在线播放| 免费女性裸体啪啪无遮挡网站| 夫妻午夜视频| 女人高潮潮喷娇喘18禁视频| 亚洲av欧美aⅴ国产| 中文精品一卡2卡3卡4更新| 国精品久久久久久国模美| 国语对白做爰xxxⅹ性视频网站| 一级a爱视频在线免费观看| 丝瓜视频免费看黄片| 亚洲成av片中文字幕在线观看| 亚洲欧美精品自产自拍| 成人国产av品久久久| 欧美+亚洲+日韩+国产| 青春草视频在线免费观看| 久久亚洲国产成人精品v| 国产精品 国内视频| 国产成人91sexporn| 日韩中文字幕视频在线看片| 国产高清videossex| 高清欧美精品videossex| 人妻 亚洲 视频| 如日韩欧美国产精品一区二区三区| 国产成人精品久久二区二区免费| 巨乳人妻的诱惑在线观看| 久久99精品国语久久久| 又黄又粗又硬又大视频| 午夜福利一区二区在线看| 99热国产这里只有精品6| 日本欧美国产在线视频| 亚洲中文av在线| 又粗又硬又长又爽又黄的视频| av电影中文网址| 日本欧美视频一区| 又大又黄又爽视频免费| 久久人人爽人人片av| 视频区欧美日本亚洲| 极品人妻少妇av视频| 国产精品99久久99久久久不卡| 日本黄色日本黄色录像| 国产片内射在线| 日韩制服丝袜自拍偷拍| 汤姆久久久久久久影院中文字幕| 9191精品国产免费久久| 老汉色av国产亚洲站长工具| 男女国产视频网站| 免费少妇av软件| 亚洲色图 男人天堂 中文字幕| 一本大道久久a久久精品| 黄片小视频在线播放| 菩萨蛮人人尽说江南好唐韦庄| 一级a爱视频在线免费观看| 亚洲av日韩精品久久久久久密 | 久久久久久久大尺度免费视频| 狠狠婷婷综合久久久久久88av| 久久人人爽av亚洲精品天堂| h视频一区二区三区| 一区在线观看完整版| 不卡av一区二区三区| 在线精品无人区一区二区三| 黑人猛操日本美女一级片| 久久久久久久国产电影| 日韩大码丰满熟妇| 欧美精品啪啪一区二区三区 | 精品免费久久久久久久清纯 | 精品国产乱码久久久久久男人| 精品第一国产精品| 午夜影院在线不卡| 亚洲伊人久久精品综合| 亚洲精品中文字幕在线视频| 男女边摸边吃奶| 真人做人爱边吃奶动态| 亚洲精品久久成人aⅴ小说| 美女主播在线视频| 亚洲午夜精品一区,二区,三区| 国产三级黄色录像| 久久人人爽av亚洲精品天堂| 精品人妻在线不人妻| 又大又黄又爽视频免费| 蜜桃国产av成人99| 国产黄色免费在线视频| 在线 av 中文字幕| 精品国产一区二区三区久久久樱花| 成人影院久久| 国产成人系列免费观看| 老司机午夜十八禁免费视频| 亚洲精品中文字幕在线视频| 最新的欧美精品一区二区| 无限看片的www在线观看| 一二三四社区在线视频社区8| 日韩中文字幕视频在线看片| 精品人妻一区二区三区麻豆| 日本一区二区免费在线视频| 美女视频免费永久观看网站| 亚洲激情五月婷婷啪啪| 一边摸一边抽搐一进一出视频| 视频区欧美日本亚洲| 在线观看国产h片| 91成人精品电影| 如日韩欧美国产精品一区二区三区| 国产老妇伦熟女老妇高清| videosex国产| 高潮久久久久久久久久久不卡| 青春草视频在线免费观看| 国产成人一区二区三区免费视频网站 | 欧美日韩av久久| 成人免费观看视频高清| 50天的宝宝边吃奶边哭怎么回事| 久久天堂一区二区三区四区| 国产成人啪精品午夜网站| 涩涩av久久男人的天堂| 久9热在线精品视频| a 毛片基地| 成在线人永久免费视频| 免费一级毛片在线播放高清视频 | 日本欧美视频一区| 日韩大片免费观看网站| 精品人妻一区二区三区麻豆| 久久这里只有精品19| 在线观看www视频免费| 男女下面插进去视频免费观看| 精品卡一卡二卡四卡免费| 午夜福利,免费看| 亚洲国产中文字幕在线视频| 亚洲第一av免费看| 一本色道久久久久久精品综合| 亚洲精品国产av蜜桃| 视频在线观看一区二区三区| 欧美日韩av久久| 一级毛片黄色毛片免费观看视频| 国产免费一区二区三区四区乱码| 成人免费观看视频高清| 我的亚洲天堂| 人人妻人人爽人人添夜夜欢视频| 国产成人欧美在线观看 | 手机成人av网站| 美女午夜性视频免费| 国产成人精品久久二区二区91| 国产片内射在线| 黄色 视频免费看| 亚洲九九香蕉| 成人三级做爰电影| 久久国产精品影院| 国产成人一区二区在线| 久久久精品国产亚洲av高清涩受| 欧美国产精品va在线观看不卡| 热99久久久久精品小说推荐| 日本色播在线视频| 人妻人人澡人人爽人人| 久久精品aⅴ一区二区三区四区| 欧美精品一区二区免费开放| 亚洲中文日韩欧美视频| 国产一区有黄有色的免费视频| 国产精品久久久久久人妻精品电影 | 一边摸一边做爽爽视频免费| h视频一区二区三区| 亚洲天堂av无毛| 99国产精品一区二区三区| 久久国产精品男人的天堂亚洲| 欧美在线一区亚洲| 久久久久久久久免费视频了| 亚洲精品一区蜜桃| 亚洲av欧美aⅴ国产| 9色porny在线观看| 国产在线一区二区三区精| 大话2 男鬼变身卡| 成人18禁高潮啪啪吃奶动态图| 少妇粗大呻吟视频| 久久人妻熟女aⅴ| 丰满迷人的少妇在线观看| 日韩制服丝袜自拍偷拍| 母亲3免费完整高清在线观看| xxxhd国产人妻xxx| 免费黄频网站在线观看国产| 亚洲七黄色美女视频| 国产片内射在线| 免费在线观看影片大全网站 | 一本色道久久久久久精品综合| 久久久久国产精品人妻一区二区| 日本黄色日本黄色录像| 欧美国产精品一级二级三级| 性色av一级| 欧美激情高清一区二区三区| 欧美精品av麻豆av| 99久久人妻综合| 欧美精品人与动牲交sv欧美| 脱女人内裤的视频| 天天躁狠狠躁夜夜躁狠狠躁| 十八禁人妻一区二区| 丁香六月欧美| 丁香六月天网| 欧美人与性动交α欧美精品济南到| 丰满饥渴人妻一区二区三| 丝袜人妻中文字幕| 热re99久久精品国产66热6| 亚洲精品日韩在线中文字幕| 亚洲一区中文字幕在线| 五月开心婷婷网| 欧美中文综合在线视频| 99热全是精品| 久久久久精品人妻al黑| 亚洲av男天堂| 国产视频一区二区在线看| 久久久亚洲精品成人影院| 手机成人av网站| 交换朋友夫妻互换小说| 丝袜脚勾引网站| 久久精品国产亚洲av高清一级| 一区二区三区乱码不卡18| 人人妻人人澡人人爽人人夜夜| 亚洲中文日韩欧美视频| 丝袜美腿诱惑在线| 精品久久久精品久久久| 久久中文字幕一级| 五月天丁香电影| 99国产精品一区二区蜜桃av | 亚洲精品自拍成人| 国产亚洲一区二区精品| 99精品久久久久人妻精品| 老汉色∧v一级毛片| 欧美人与性动交α欧美精品济南到| 亚洲国产欧美日韩在线播放| 一区二区av电影网| 亚洲精品一区蜜桃| 精品久久蜜臀av无| 两个人免费观看高清视频| www.熟女人妻精品国产| 蜜桃在线观看..| 国产精品一区二区在线不卡| 人人妻人人爽人人添夜夜欢视频| 啦啦啦在线免费观看视频4| 美女高潮到喷水免费观看| 日日摸夜夜添夜夜爱| 亚洲精品美女久久av网站| 丁香六月欧美| 亚洲九九香蕉| a级毛片黄视频| a级毛片在线看网站| 亚洲精品第二区| 大香蕉久久成人网| 在线亚洲精品国产二区图片欧美| 日韩一卡2卡3卡4卡2021年| 精品亚洲成国产av| 成年女人毛片免费观看观看9 | 别揉我奶头~嗯~啊~动态视频 | 国产男女超爽视频在线观看| 午夜福利影视在线免费观看| 欧美在线一区亚洲| 建设人人有责人人尽责人人享有的| 久久精品亚洲av国产电影网| 三上悠亚av全集在线观看| 精品国产一区二区久久| 在线观看一区二区三区激情| 多毛熟女@视频| 国产一卡二卡三卡精品| 国产av国产精品国产| 亚洲成色77777| 亚洲欧美一区二区三区久久| 国产成人精品无人区| 亚洲欧美精品综合一区二区三区| 欧美精品一区二区免费开放| 韩国精品一区二区三区| 久久久久视频综合| 亚洲第一av免费看| 国产高清国产精品国产三级| 亚洲国产日韩一区二区| 国产1区2区3区精品| 欧美日本中文国产一区发布| 激情五月婷婷亚洲| 精品亚洲成国产av| 欧美亚洲 丝袜 人妻 在线| 一本久久精品| 视频区欧美日本亚洲| 欧美乱码精品一区二区三区| 亚洲黑人精品在线| 欧美成人午夜精品| 女警被强在线播放| 国产一区有黄有色的免费视频| 日本猛色少妇xxxxx猛交久久| 精品少妇一区二区三区视频日本电影| 精品高清国产在线一区| 免费看不卡的av| 中文精品一卡2卡3卡4更新| 乱人伦中国视频| av福利片在线| 亚洲精品国产区一区二| 最新的欧美精品一区二区| 在现免费观看毛片| 天堂中文最新版在线下载| 亚洲国产精品一区三区| 精品高清国产在线一区| 久久亚洲国产成人精品v| 2018国产大陆天天弄谢| 97在线人人人人妻| bbb黄色大片| 日韩中文字幕欧美一区二区 | 国产精品欧美亚洲77777| 乱人伦中国视频| 91国产中文字幕| 午夜两性在线视频| 日韩伦理黄色片| 极品少妇高潮喷水抽搐| 欧美国产精品一级二级三级| 看免费成人av毛片| 亚洲自偷自拍图片 自拍| av网站免费在线观看视频| 亚洲成人免费电影在线观看 | 精品国产一区二区三区四区第35| 精品第一国产精品| 久久久国产精品麻豆| www日本在线高清视频| 色精品久久人妻99蜜桃| 女人被躁到高潮嗷嗷叫费观| 中文字幕色久视频| 欧美激情 高清一区二区三区| 亚洲精品av麻豆狂野| 免费观看av网站的网址| 亚洲三区欧美一区| 亚洲国产av新网站| h视频一区二区三区| 欧美久久黑人一区二区| 午夜福利在线免费观看网站| 欧美精品啪啪一区二区三区 | 婷婷色综合大香蕉| 另类精品久久| 欧美日韩成人在线一区二区| 水蜜桃什么品种好| 亚洲黑人精品在线| 国产片特级美女逼逼视频| 中文字幕另类日韩欧美亚洲嫩草| 脱女人内裤的视频| 1024香蕉在线观看| 精品少妇内射三级| 亚洲国产欧美一区二区综合| 午夜精品国产一区二区电影| av又黄又爽大尺度在线免费看| 亚洲第一青青草原| 99久久99久久久精品蜜桃| 欧美精品一区二区大全| 丰满饥渴人妻一区二区三| 国产一区二区在线观看av| avwww免费| 国产99久久九九免费精品| 人体艺术视频欧美日本| 亚洲精品第二区| 少妇人妻 视频| 国产精品香港三级国产av潘金莲 | 亚洲九九香蕉| 亚洲欧美中文字幕日韩二区| 久久精品国产亚洲av高清一级| 夜夜骑夜夜射夜夜干| 国产有黄有色有爽视频| 一本综合久久免费| 免费观看av网站的网址| 亚洲自偷自拍图片 自拍| 久久99一区二区三区| 狂野欧美激情性bbbbbb| 日本午夜av视频| 亚洲欧美中文字幕日韩二区| 美女大奶头黄色视频| a级毛片黄视频| 国产女主播在线喷水免费视频网站| 欧美乱码精品一区二区三区| 肉色欧美久久久久久久蜜桃| 国产在视频线精品| 老司机影院毛片| 午夜日韩欧美国产| 精品免费久久久久久久清纯 | 亚洲视频免费观看视频| 男女边吃奶边做爰视频| 日韩 欧美 亚洲 中文字幕| 日本猛色少妇xxxxx猛交久久| 精品亚洲成国产av| 性少妇av在线| 亚洲av成人不卡在线观看播放网 | 七月丁香在线播放| 精品人妻一区二区三区麻豆| 十分钟在线观看高清视频www| 亚洲国产毛片av蜜桃av| 啦啦啦在线观看免费高清www| 欧美在线黄色| 国产男女内射视频| 亚洲精品一二三| 丁香六月欧美| 久久人妻福利社区极品人妻图片 | 久久精品国产亚洲av高清一级| 亚洲欧美中文字幕日韩二区| 十八禁人妻一区二区| 老司机靠b影院| 日本wwww免费看| 天天躁狠狠躁夜夜躁狠狠躁| 日韩中文字幕视频在线看片| 国产av一区二区精品久久| 亚洲欧洲精品一区二区精品久久久| 99国产精品免费福利视频| 久久 成人 亚洲| 精品国产一区二区三区久久久樱花| 久久精品国产综合久久久| 丝袜美腿诱惑在线| 亚洲精品日本国产第一区| av网站免费在线观看视频| 又紧又爽又黄一区二区| 十八禁人妻一区二区| 91字幕亚洲| 各种免费的搞黄视频| 久久天堂一区二区三区四区| 黄色毛片三级朝国网站| 波野结衣二区三区在线| 中文字幕精品免费在线观看视频| 国产精品久久久久久人妻精品电影 | 巨乳人妻的诱惑在线观看| 国产激情久久老熟女| 欧美老熟妇乱子伦牲交| 国产成人免费无遮挡视频| 一级毛片黄色毛片免费观看视频| 97人妻天天添夜夜摸| 日本色播在线视频| 亚洲少妇的诱惑av| 国产激情久久老熟女| 国产淫语在线视频| 久久国产精品人妻蜜桃| 国产精品.久久久| 亚洲精品av麻豆狂野| 少妇裸体淫交视频免费看高清 | 国产高清视频在线播放一区 | 免费人妻精品一区二区三区视频| 性色av一级| bbb黄色大片| 巨乳人妻的诱惑在线观看| av天堂久久9| 国产欧美日韩精品亚洲av| 国产成人av激情在线播放| 国产精品二区激情视频| 精品第一国产精品| 午夜福利,免费看| 国产精品一区二区免费欧美 | av欧美777| 国产人伦9x9x在线观看| 岛国毛片在线播放| 老司机深夜福利视频在线观看 | 视频区图区小说| 国产日韩欧美在线精品| 亚洲国产中文字幕在线视频| 午夜免费男女啪啪视频观看| 亚洲 欧美一区二区三区| 色94色欧美一区二区| 亚洲欧美清纯卡通| 伊人亚洲综合成人网| 亚洲五月色婷婷综合| 亚洲熟女精品中文字幕| 性色av一级| 黄色片一级片一级黄色片| 日本91视频免费播放| 中文字幕最新亚洲高清| 一二三四社区在线视频社区8| 国语对白做爰xxxⅹ性视频网站| 一级,二级,三级黄色视频| 在线亚洲精品国产二区图片欧美| 777米奇影视久久|