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

    Particle Image Velocimetry Measurement of the Flow Field in the Play of the Drilling Pump Valve

    2010-03-01 01:46:54YANGGuoanYINXinSONGZhengandHUANGCong

    YANG Guoan, YIN Xin, , SONG Zheng, and HUANG Cong

    1 Diagnosis & Self-recovery Engineering Research Center, Beijing University of Chemical Technology,Beijing 100029, China

    2 Production Operations Center, PetroChina Guangxi Petrochemical Company, Qinzhou 535008, China

    1 Introduction

    The drilling pump valve is one of the most critical components in a drilling pump, having a direct impact on the efficiency and reliability of the whole machine.However, due to the complex working conditions in practical use, the service life of the valve is the shortest among all the wearing parts of a drilling pump, which restricts the development of the drilling pump in the direction of high-pressure, high-stroke and high-output.

    Over the decades, a plenty of studies were carried out on the issue of valve failure, however, due to the effect of a complex system of forces on the valve, there are difficulties in the theoretical research and experimental study. Judging from the literature, it is found most of the studies on the pump valve were directly aimed at the valve itself,concentrating on exploring the mechanism of valve failure and improving valve design. Especially in the discussion of the failure mechanism, many points of view have been presented. For example, JOFFE, et al[1], believed that the main causes of the valve failure were abrasive wear, fatigue wear as well as the combined effect of cavitation and fatigue wear; in Ref. [2], analyses are made respectively on different parts of the valve, and it is believed that the causes of valve failure are different for each part of the valve; With the quantitative calculation of fault tree and an analysis of the damaged valve surface, WANG, et al[3], presented that the slurry erosion and the material of the valve had the greatest impact on the valve failure; differently, as to the work of XIA, et al[4], the study was from the perspective of impact fatigue, putting forward experimental method to measure the impact force. Researches have shown that many cases of valve failure are caused by the erosion of slurry. However, since there is no experimental study special on the flow field of the valve play, it is much lacking in the knowledge of the flow field characteristics,and thus research work in this aspect is hard to get furthered.

    In this paper, a particle image velocimetry(PIV)measurement of the flow field in the model valve was performed to experimentally investigate the characteristics of the flowing fluid in the valve play, which verified and improved the previous computational fluid dynamics(CFD)simulations. This paper firstly introduces the preparation work for the measurement, including the design of the experimental setup, particle selection, and adjustment of PIV devices, such as the camera and the laser. Then,analyses are made on the obtained PIV data for different working conditions and conclusions are drawn from these results. After that, in order to get more precise data, an advanced PIV test with fluorescent particles is introduced as a supplement. Finally in this article, the measuring accuracy is discussed and the conclusions are summarized.

    2 Design of Experimental Setup

    Fig. 1 shows a sketch of the experimental setup. This experimental setup basically consists of a pump, a tank,pipelines and a model valve with three pairs of scaled valve discs and valve seats. Among all these parts, design for the body of the model valve and the valve discs and seats therein are keys to the whole installation, taking the 7#drilling pump valve commonly used in oilfields as a reference and scaled down with proportion of 1:4. These optically transparent models are made of acrylic glass with polished external and internal surfaces, mainly assembled with three parts which are the upper chamber, the lower chamber and the valve disc and seat in between.

    Fig. 1. Sketch of the experimental setup

    Considering the structure of the model and the fluid pressure under the working condition of the experiment,there are two assumptions introduced in the design of the models: (1) stress σrin the radial direction is negligible, as the body model is a thin-walled vessel and the scale of fluid pressure therein is not higher than 1 MPa; (2) the impact stress occurring with valve closing could be neglected,because the shock of the valve disc is tiny in the experimental model. On the basis of the proportion mentioned above and the design criterion for pressure vessels[5], the figure and structure of the designed model are decided as what’s shown in Fig. 2. Fig. 3 shows the model valve discs and seats. As shown in Fig. 4, the external wall of the model is designed as square, which helps to minimize the refraction at the surface.

    3 Experimental Procedures

    A two-dimenional PIV system was used to measure the velocity of the flow within the model. The light sheet was emitted by a integrated high power Nd-YAG duolaser, and the images of the particles were recorded by a chargecoupled device(CCD) camera, Redlake Megaplus II ES2001(resolution 1 600×1 200 pixels, frame rate 31.85 f/s). The operation of the laser and the camera was synchronized by a digital Delay/Pulse Generator, model MicroPulse710.Since the region to be photographed was as small as an area of 20 mm×20 mm and the interrogation window was even narrower and smaller of which the maximal width was less than 5 mm and the flow speed thereof was much higher, a conventional camera lens was not suitable for this shooting.So, a microfocus lens (model Nikon AF-S VR Micro-Nikkor 105 mm f/2.8G IF-ED) was used to get enlarged images.

    Fig. 2. Designed model

    Fig. 3. Photos of model valve discs & seats

    Fig. 4. Cross section of the model

    The particle type selection and the particle seeding are two crucial issues for PIV measurement. To select a suitable type of particles for this test, the scattering performance and the uniformity of the beads dispersing should be taken into account[6–7]. A type of hollow glass beads (density 1.05 g/cm3and diameter 1–5 μm) were chosen as tracers for this test. For the seeding of the particles, the appropriate dosage of the particles should be seeded into the circulating flow. The number of the particles within an interrogation region should be sufficient to describe the flow state, and it should be noted if the number of the particles is too large, speckles may form,which may cover the tracers[7]. However, it is not easy to figure out such a right dosage quantitatively, as many factors, such as laser power, model structure, etc, may influence the quality of PIV images. In the tests, the dosage of the particles was determined in such a way: adding particles into water bit by bit, and series of real-time images of the seeding beads were monitored and displayed on the computer screen to seed the beads, from which the right dosage was decided when the desirable picture was obtained.

    In the model, positional relation between the laser and the camera is shown in Fig. 5.

    Fig. 5. Positional relations among the model, laser and camera

    Due to the light refraction on the conical surfaces of valve seats, if the camera shot vertically into the model,images of the valve play could not be clearly photographed,and the gap between the disc and the seat shown in the recorded pictures would be more narrow than it actually was, which might cause some difficulties in computing velocities in the next step. To avoid these, after repeated trials, the method of off-axis measurement (see Fig. 6) was finally adopted. It has been proved that when the camera and the model’s symmetry axis were placed at an angle of about 20°–25°, the influence of the refraction caused by the conical structure would get minimized, which contributed to get better and clearer images of the valve play.

    Fig. 6. Vertical shot and off-axis shot

    In addition, the region where the flow field is investigated is the valve play, i.e., the annular gap between the conical surface of a valve disc and that of a seat. If the camera were laid horizontally, images of the valve play would be correspondingly angled towards the horizontal position, which might also cause inconvenience to marking the interrogation area. As shown in Fig. 7, the shooting angle of the camera in the vertical plane could be adjusted to make the images horizontal.

    Fig. 7. Angle relationship between the camera and particle image

    4 Processing and Analysis of Experimental Results

    In order to explore the velocity distribution in the flow field for different valve angles (the valve angle refers to the angle between the axis and the generatrix of the valve disc)and for different valve lifts, experiments were conducted for 9 working conditions. For each valve seat and valve disc, tests were performed at lifts of 3 mm, 4 mm and 5 mm,respectively.

    When the particle images were acquired, MicroVec2 software was used to process these raw picture data. Within the software, the speed of the particles was essentially computed on the basis of cross-correlation, and additionally techniques, such as iterative algorithm, erroneous vector correction, were employed to reduce the computing error.By using the software Tecplot, vector data obtained in the computations, including velocity, velocity distributions and vorticity distributions, could be plotted as shown in Fig. 8.Fig. 8 shows the flow field state of 45° valve at the valve lift of 4 mm. As shown in Fig. 8(a), the flow direction changes in the vicinity of position A, i.e., the base angle of the valve disc. Compared with Fig. 8(b) at the same position, it is found that the maximal velocity of the entire field takes place near the corner of the disc.

    Fig. 8. Flow field state of 45° valve at lift of 4 mm

    As illustrated in Fig. 9 in detail, when the water flows upwards and approaches the entry to the valve play where the cross-section of the flow channel is getting narrow gradually, some of the fluids will flow in the direction of the valve play, whereas others will be forced to move straight to the bottom of the valve disc. Because of this, a fluid pressure in the area beneath the bottom of the disc arises and increases gradually, and also the internal energy thereof accumulates, which forces these impeded fluids to rush into the valve play area. As can be seen in Fig. 9, from section B-B' to A-A', width of the flow channel continuously contracts and reaches the narrowest width near the section A-A'. According to the continuity principle,it can be concluded that the maximal speed occurs in the vicinity of section A-A'. The conclusion is applicable to other working conditions.

    Fig. 9. Sketch of flow movement near the entry to valve play

    In further, the rushing flow from under the disc bottom will be squeezed by that of the lower part in the passage and its flow direction will be turned along the valve play.As a result of fluid inertia, this turning of direction occurs mostly in the middle part of the entrance. As shown in Fig.9, on the upper wall of the channel, fluids in the proximity of position C are protected by both the bottom of disc and the exterior faster flow, less affected by the subsequent flow.Hence it can be deduced that substances therein exchange much more slowly with those of the outer liquid, which contributes to the forming of vortexes. So it can be speculated here that, for the practical operating condition of a pump, fluid therein is drilling slurry containing lots of solid granules, therefore eddies at position C will probably lead to the remain of those grains, which may become one of the factors inducing the valve failure.

    Looking through the graphs of flow fields for different conditions (see Fig. 10–Fig. 17), obviously, the color of the middle part of the channel is brighter than that near the two walls. The brighter color indicates the higher mean flow velocities and more uniform flow vectors. However, in the vicinity of the upper and lower walls, the color is darker,which means the lower speed and slower change in flow directions. These results are analogous to those from previous CFD simulations[8]that the highest velocity exists at the midpoint of the cross-section and gradually decreases to both sides. Moreover, from velocity profiles, it can be concluded that the flow velocities near the entry to the valve play are remarkably higher than those near the exit. In practice, this is helpful that the sealing ring shall be fixed on the upside of the valve disc to reduce the erosion.

    Fig. 10. Flow field state of 30° valve at lift of 3 mm

    Fig. 11. Flow field state of 35° valve at lift of 3 mm

    Fig. 12. Flow field state of 45° valve at lift of 3 mm

    Fig. 13. Flow field state of 30° valve at lift of 4 mm

    Fig. 14. Flow field state of 35° valve at lift of 4 mm

    Fig. 15. Flow field state of 30° valve at lift of 5 mm

    Fig. 16. Flow field state of 35° valve at lift of 5 mm

    Fig. 17. Flow field state of 45° valve at lift of 5 mm

    From vorticity profiles (see Fig. 18–Fig. 26), it can be inferred that the absolute values of vorticity in the adjacent areas of both walls of the channel are much higher than those in the middle. Vorticity values near the upper wall are positive while the lower are negative. This indicates that the flow direction changes a lot in the vicinity of the both walls of the passage, usually with an anticlockwise turning near the upper wall or with a clockwise turning near the lower wall. Considering the actual working conditions of a pump, as drilling slurry passes through the valve play, both of the turnings mentioned above will cause a phenomenon that the fluids with plenty of solid granules will strike the surfaces of the valve disc and seat.

    Fig. 18. Vorticity profile of 30° valve at lift of 3 mm (s–1)

    Fig. 19. Vorticity profile of 35° valve at lift of 3 mm (s–1)

    Fig. 20. Vorticity profile of 45° valve at lift of 3 mm (s–1)

    Fig. 21. Vorticity profile of 30° valve at lift of 4 mm (s–1)

    Fig. 22. Vorticity profile of 35° valve at lift of 4 mm (s–1)

    Fig. 24. Vorticity profile of 30° valve at lift of 5 mm (s–1)

    Fig. 25. Vorticity profile of 35° valve at lift of 5 mm (s–1)

    Fig. 26. Vorticity profile of 45° valve at lift of 5 mm (s–1)

    After that, the granules slow down, and some of them congregate near the walls and the amount may increase with time. When the disc impacts the seat with speed of about 20 m/s[9], some of the massing granules may not be washed out of the play immediately. The remained solids might be compressed heavily in between the disc and the seat, damaging both of the surfaces.

    As summarized in Table, the maximal velocities and the mean velocities of each flow field are obtained after a classification of the data. Based on Table, trend charts for both the maximal velocities and the mean velocities are respectively plotted.

    Table. Max. velocities and mean velocities of each flow field

    From Fig. 27(a), it can be seen that when the valve lift is constant, the maximal speed in the flow field shows a tendency of decreasing with the angle of the valve disc.This experimental result proves the CFD software computed result in Ref. [8] which indicates that, when the diameter of the disc bottom is constant (and so is the diameter of the valve hole), velocity in the flow field decreases with the increase in the valve angle. The scale of average velocities could be used to evaluate the rapidness of the flow. As shown in Fig. 28(a), when the valve angle increases, the mean speed also exhibits a downward trend.

    Similarly, in Fig. 27(b), the plot of maximal speed exhibits a trend of reducing when the valve disc rises. It is interesting to note here that, as shown in this chart, the curve of the velocities is somewhat approximate to that of a negative exponential function, with a larger decreasing rate of velocity at shorter valve lifts but a smaller one when the lifts get higher. According to Fig. 28(b), the variation in the average velocity is still consistent with the general trend of decreasing, and the speed of the flow reduces when the valve lifts increase. Hence a conclusion could be drawn that,the flow velocity in the valve play is declining continuously while the valve disc is rising, and the rate of the decreasing is getting smaller gradually. Combined with the valve failure mechanisms, it can be inferred that when a working valve is about to open or to close, its lift is tiny and the flow velocity is extremely high. With such high-speed drilling slurry containing abrasive solid particles, serious erosion damage to valves, especially to sealing rings is ineluctable and the injuries may accelerate the failure of the valve.

    Fig. 27. Variations of maximal velocity

    Fig. 28. Variations of mean velocity

    In order to get more precise data, an improved measurement, in which fluorescent particles were introduced, was carried out for the 30° valve at the lift of 4 mm. The fluorescent particles, of which the main component is Rhodamine B, are usually seeded in Micro-PIV experiments. Unlike those ordinary tracer particles, in the experimental tests with fluorescent particles,when these tracers are illuminated by laser, the fluorescence with a wavelength of around 590 nm can be excited. Meanwhile, when recording the particle images,the camera shoots into the flow with a bandpass color filter fixed in front of its lens, therefore visible light of other wavelengths could be effectively filtered out and its harmful influences on the images would get minimized.

    Results of this measurement are presented in Fig. 29, in which A-A' denotes the section of the entrance to the valve play.

    By analyzing these graphs, the following conclusions can be drawn.

    (1) With the velocity distribution (see Fig. 29(b)), it is evident that when the water flows into the area of the channel(above Section A-A'), its velocity increases rapidly and approximately shows a parabolic layered distribution,which experimentally verifies the CFD result computed in Ref. [8];

    Fig. 29. Results of measurement with fluorescent particles

    (2) Fig. 29(d) presents a distribution of the velocity fluctuation. Obviously the values of velocity fluctuation in the channel area (above section A-A') are higher than those of the outer (below section A-A'), and generally the velocity fluctuation in the area near the wall is higher than that in the middle. Besides, there exists the highest velocity fluctuation in the proximity of position B, and that is the site where the actual valve failure initially appears.

    Compared with experiments using hollow glass beads, it is found that results from this test are consistent with the previous conclusions.

    5 Discussion of Measuring Accuracy

    The accuracy of PIV may be affected by three factors:the precision of the PIV system, the following behavior of tracing particles, and the refraction of light, which are detailed as follows.

    (1) Precision of the PIV system. The precision of the PIV system here mainly refers to the measuring accuracy of the computed results by processing particle images with the PIV software. With an experimental assessing approach presented in Ref. [10], the precision of the MicroVec2 software has been calibrated, and the average relative error is 1%.

    (2) The following behavior of tracing particles. Results computed by the PIV system are actually the movement of those tracing particles, and that is why the following behavior must be considered when selecting the particles.In the work of LI, et al[11], an analysis of the following behavior was presented. And judging from their results, the particles used in this measurement, whose density is 1.05 g/cm3and diameter is 1–5 μm, have very good following feature, which is suitable for the test.

    (3) Light refraction. During the collection of particle images, what the camera records were not the original images of the particles but the refracted ones penetrating the water and the acrylic glass, and that directly affects the measurement results. Since the refractive index of acrylic glass is about 1.49 and that of the water is about 1.33 (after the particles being seeded, the index will be a little higher),the angle of refraction from water to acrylic glass will be very small. However, the angle will be larger when light travels from acrylic glass to the air, and thus an error will be caused. As a solution of this, the external wall of the model is designed as square to minimize this negative influence, which has already been mentioned in Fig. 4.

    6 Conclusions

    (1) Near the base angle of the valve disc, flow direction changes distinctively, with a rapid velocity and serious erosion at this position. However,in the particular region which is marked as C in Fig. 9, an area of slower flow exists there, which may be easy to form eddies and cause the remaining of solid granules.

    (2) The central region of the flow field has larger average flow velocities and more uniform flow vectors.Contrastively, in the vicinity of the both walls of the valve play, flow speeds thereof slow down and the direction changes greatly. In addition, the experimental test has proved that the velocities near the entrance to the play are significantly higher than that near the exit, thus sealing rings of the pump valve should be fixed on the upper part of the valve disc to reduce the erosion.

    (3) The absolute values of the vorticity near the walls are much higher than those in the middle area, and the upper has positive values whilst the lower has negative values.Therefore, in practice, some of solid particles in drilling slurry may gather near both of the walls, and when the pump valve closes, they may be compressed, leading to the failure of the valve.

    (4) A conclusion of previous CFD simulations (in Ref. [8]) has been verified that when the diameter of the disc bottom is constant (and so is the diameter of the valve hole), velocity in the flow field decreases with the increase in the valve angle. Besides, a group of plots have been presented in Fig. 27(b), from which it could be speculated that serious erosion would happen when the disc gets a smaller lift.

    (5) With the results from the improved measurement with fluorescent particles, a parabolic layered distribution of velocities has been plotted, which experimentally verifies the previous CFD results (in Ref. [8]); it has been measured that values of velocity fluctuation in the valve play are higher than that out of it, and the highest fluctuation occurs on the site where the actual valve failure initially appears.

    [1] JOFFE S H D, ALLEN C. The wear of pump valves in fine particle quartzite slurries[J]. Wear, 1990, 135(2): 279–291.

    [2] YE Yongbiao, ZHU Yongyou, WANG Longlong. The failure analysis of the valve in reciprocating pump[J]. Water Conservancy& Electric Power Machinery, 2006, 28(1): 37–38. (in Chinese)

    [3] WANG Zhonghui, LU Yongming, ZHANG Jinzhong. Fault tree analysis for failure of drilling pump valve[J]. Oil Field Equipment,1995, 24(1): 26–27. (in Chinese)

    [4] XIA Yuanbai, LUO Hui. Testing method of the impacting force of drilling pump valve in the course of valve closing[J]. Oil Field Equipment, 2001, 30(6): 30–32. (in Chinese)

    [5] ZHENG Jinyang, DONG Qiwu, SANG Zhifu. Process equipment design[M]. Beijing: Chemical Industry Press, 2001. (in Chinese)

    [6] MEINHART C D, WERELEY S T, SANTIAGO J G. PIV measurements of a microchannel flow[J]. Experiments in Fluids,1999, 27: 414–419.

    [7] WEN Jian, LI Yanzhong, WANG Simin, et al. Investigation of flow patterns in the entrance of heat exchanger by PIV technology[J].Journal of Harbin Institute of Technology, 2008, 40(1): 113–117. (in Chinese)

    [8] YANG Guoan, HUANG Cong, YU Li. Compute and analysis of the flow field in the play of the drilling pump valve based on the simulation by FLUENT[J]. Oil Field Equipment, 2008, 38(3):41–44. (in Chinese)

    [9] YANG Guoan, ZHANG Dong, HUANG Cong. The impact characteristic analysis of drilling pump valve[J]. Journal of Vibration and Shock, 2008, 27(12): 18–22. (in Chinese)

    [10] DONG Mingzhe, WANG Yang, JIANG Ningtao. An experimental method for evaluating the PIV measurement error and its application in optimizing PIV experiment parameters[J]. Journal of Experiments in Fluid Mechanics, 2005, 19(2): 79–83. (in Chinese)

    [11] LI Enbang, LI Zhiping, LI Chun, et al. Numerical analysis of following behaviors of particle tracers in turbulent[J]. Chinese Journal of Scientific Instrument, 2009, 30(2): 225–231. (in Chinese)

    日本精品一区二区三区蜜桃| 久99久视频精品免费| 蜜桃久久精品国产亚洲av| 国产av麻豆久久久久久久| 三级经典国产精品| 亚洲综合色惰| 亚洲av第一区精品v没综合| 精品无人区乱码1区二区| 亚洲在线自拍视频| 成年女人看的毛片在线观看| 久久99热这里只有精品18| 综合色丁香网| 国产v大片淫在线免费观看| 国产在线男女| 99国产极品粉嫩在线观看| 变态另类丝袜制服| 日本一二三区视频观看| 天堂网av新在线| 在线观看66精品国产| 97人妻精品一区二区三区麻豆| 日韩强制内射视频| 日本撒尿小便嘘嘘汇集6| av专区在线播放| 国产午夜福利久久久久久| 美女免费视频网站| 女人十人毛片免费观看3o分钟| 成年免费大片在线观看| 精品少妇黑人巨大在线播放 | 欧洲精品卡2卡3卡4卡5卡区| 三级经典国产精品| 久久久成人免费电影| 日本a在线网址| 久久99热6这里只有精品| 亚洲av中文字字幕乱码综合| 亚洲精品在线观看二区| 在线观看美女被高潮喷水网站| 国产黄色视频一区二区在线观看 | 亚洲图色成人| 日本成人三级电影网站| 国产精品免费一区二区三区在线| 日韩亚洲欧美综合| 99久久精品一区二区三区| 我的老师免费观看完整版| 丰满的人妻完整版| 国产亚洲精品久久久com| 天堂影院成人在线观看| 1000部很黄的大片| 狂野欧美激情性xxxx在线观看| 国产精品美女特级片免费视频播放器| 国产色婷婷99| 亚洲激情五月婷婷啪啪| 欧美性感艳星| 欧美极品一区二区三区四区| 亚洲经典国产精华液单| 色吧在线观看| 成人av在线播放网站| 国产探花在线观看一区二区| 日本五十路高清| 综合色av麻豆| 男人狂女人下面高潮的视频| 亚洲色图av天堂| 男人舔女人下体高潮全视频| 两个人视频免费观看高清| 欧美3d第一页| 久久久久国产网址| 一进一出好大好爽视频| 欧美成人一区二区免费高清观看| 亚洲av成人精品一区久久| 欧美三级亚洲精品| 国产亚洲91精品色在线| 香蕉av资源在线| 毛片一级片免费看久久久久| 国产精品一区二区免费欧美| 色哟哟哟哟哟哟| 五月玫瑰六月丁香| 成人av在线播放网站| 可以在线观看毛片的网站| 12—13女人毛片做爰片一| 久久人人爽人人爽人人片va| 国产片特级美女逼逼视频| 国内精品久久久久精免费| 日本黄大片高清| 欧美精品国产亚洲| 少妇熟女欧美另类| 黄色配什么色好看| 在线免费观看的www视频| 亚洲av熟女| 99久国产av精品国产电影| av女优亚洲男人天堂| 中文字幕精品亚洲无线码一区| 午夜激情欧美在线| 久久精品国产亚洲av天美| 两性午夜刺激爽爽歪歪视频在线观看| 婷婷六月久久综合丁香| 久久午夜亚洲精品久久| 成人一区二区视频在线观看| 欧美日本亚洲视频在线播放| 国产一区二区亚洲精品在线观看| 日韩欧美免费精品| av在线播放精品| 一边摸一边抽搐一进一小说| 精品久久久久久久久久免费视频| 色哟哟哟哟哟哟| 麻豆国产97在线/欧美| 99久国产av精品国产电影| 精品久久久久久成人av| 国产亚洲精品av在线| 寂寞人妻少妇视频99o| 最新中文字幕久久久久| 国产国拍精品亚洲av在线观看| 亚洲国产精品成人综合色| 久久99热6这里只有精品| 国产又黄又爽又无遮挡在线| 插逼视频在线观看| 亚洲精品日韩av片在线观看| 久久人人爽人人爽人人片va| а√天堂www在线а√下载| 国产精品嫩草影院av在线观看| 午夜福利成人在线免费观看| 亚洲av中文字字幕乱码综合| 亚洲欧美精品自产自拍| 亚洲成人中文字幕在线播放| 国产国拍精品亚洲av在线观看| 亚洲国产精品成人综合色| 少妇熟女欧美另类| 久久久久久久亚洲中文字幕| 日韩国内少妇激情av| 国产成人精品久久久久久| 少妇人妻一区二区三区视频| 久久久久久国产a免费观看| 嫩草影院入口| 日本撒尿小便嘘嘘汇集6| 黄色欧美视频在线观看| 精品福利观看| 亚洲自拍偷在线| 欧美zozozo另类| 国产精品一区www在线观看| 亚洲国产精品国产精品| 亚洲最大成人中文| 1000部很黄的大片| 色吧在线观看| aaaaa片日本免费| 一个人看视频在线观看www免费| 婷婷六月久久综合丁香| 精品一区二区三区视频在线| 国产毛片a区久久久久| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲精品成人久久久久久| 亚洲婷婷狠狠爱综合网| 日韩欧美在线乱码| 国内揄拍国产精品人妻在线| 毛片女人毛片| 国产伦一二天堂av在线观看| av免费在线看不卡| 国产视频内射| 国产视频一区二区在线看| 国产成年人精品一区二区| 国产人妻一区二区三区在| 91狼人影院| 国内揄拍国产精品人妻在线| 51国产日韩欧美| 特级一级黄色大片| 亚州av有码| 一级黄片播放器| 看免费成人av毛片| 国产国拍精品亚洲av在线观看| 国产亚洲欧美98| 国内精品一区二区在线观看| 亚洲一级一片aⅴ在线观看| 亚洲国产色片| 又黄又爽又免费观看的视频| 久久精品国产清高在天天线| 秋霞在线观看毛片| 国产老妇女一区| 有码 亚洲区| 中出人妻视频一区二区| 夜夜看夜夜爽夜夜摸| a级毛片a级免费在线| 又爽又黄无遮挡网站| 亚洲经典国产精华液单| 淫秽高清视频在线观看| 精品午夜福利视频在线观看一区| 不卡一级毛片| 欧美性感艳星| 精品久久久噜噜| 国产精品久久久久久亚洲av鲁大| 午夜福利在线观看吧| 欧美一级a爱片免费观看看| 99国产极品粉嫩在线观看| 美女内射精品一级片tv| av在线亚洲专区| 男女做爰动态图高潮gif福利片| 国产亚洲av嫩草精品影院| 看片在线看免费视频| 亚洲乱码一区二区免费版| 夜夜爽天天搞| 波多野结衣巨乳人妻| 国产精品美女特级片免费视频播放器| 一进一出好大好爽视频| 日本免费a在线| 麻豆乱淫一区二区| 久久精品国产鲁丝片午夜精品| 亚洲最大成人中文| 久久人人精品亚洲av| 久久久久久久午夜电影| 亚洲欧美日韩高清专用| 国产精品一区二区免费欧美| 日本成人三级电影网站| 久久99热这里只有精品18| 亚洲高清免费不卡视频| 高清毛片免费看| 婷婷精品国产亚洲av| 国产一区二区亚洲精品在线观看| 久久久a久久爽久久v久久| 2021天堂中文幕一二区在线观| 久久久久国内视频| 黄色一级大片看看| 成人特级av手机在线观看| 久久久久久久久久久丰满| 男人舔奶头视频| 午夜福利在线观看吧| 久久亚洲国产成人精品v| 亚洲一级一片aⅴ在线观看| 国产69精品久久久久777片| 亚洲电影在线观看av| 搡老妇女老女人老熟妇| 久久久久久久久中文| 99久久精品热视频| 欧美成人免费av一区二区三区| 身体一侧抽搐| 性色avwww在线观看| 国产伦一二天堂av在线观看| 91久久精品国产一区二区三区| 欧美国产日韩亚洲一区| eeuss影院久久| 国内精品宾馆在线| 三级男女做爰猛烈吃奶摸视频| 舔av片在线| 天堂影院成人在线观看| 五月伊人婷婷丁香| 国产精品女同一区二区软件| 亚洲成a人片在线一区二区| 欧美日韩综合久久久久久| 我要看日韩黄色一级片| 成人一区二区视频在线观看| 深夜精品福利| 日韩欧美国产在线观看| 午夜a级毛片| 三级毛片av免费| 看十八女毛片水多多多| 在线免费观看的www视频| 午夜激情福利司机影院| 欧美高清性xxxxhd video| 日日啪夜夜撸| 91久久精品国产一区二区成人| 真人做人爱边吃奶动态| 婷婷六月久久综合丁香| 久久草成人影院| 99视频精品全部免费 在线| 成人漫画全彩无遮挡| 色尼玛亚洲综合影院| 国产成人freesex在线 | 久久久国产成人免费| 又黄又爽又免费观看的视频| 色哟哟·www| 天堂影院成人在线观看| av在线观看视频网站免费| 天天躁日日操中文字幕| 小蜜桃在线观看免费完整版高清| 久久精品国产亚洲av香蕉五月| 美女高潮的动态| 亚洲成人av在线免费| 国产黄色小视频在线观看| 99热这里只有是精品50| 日本与韩国留学比较| av卡一久久| 啦啦啦观看免费观看视频高清| 日韩制服骚丝袜av| 简卡轻食公司| 婷婷亚洲欧美| 淫秽高清视频在线观看| 国产色婷婷99| avwww免费| 日韩亚洲欧美综合| 全区人妻精品视频| 美女xxoo啪啪120秒动态图| 国产久久久一区二区三区| 成年女人看的毛片在线观看| 18禁在线无遮挡免费观看视频 | 成人三级黄色视频| 在线免费观看不下载黄p国产| 久久鲁丝午夜福利片| 久久久久久久久久黄片| 国产乱人偷精品视频| 国产乱人视频| 久99久视频精品免费| 色综合色国产| 久久精品综合一区二区三区| 中文字幕免费在线视频6| 亚洲中文字幕一区二区三区有码在线看| a级毛片免费高清观看在线播放| 如何舔出高潮| 夜夜看夜夜爽夜夜摸| 毛片一级片免费看久久久久| 小蜜桃在线观看免费完整版高清| 久久精品国产99精品国产亚洲性色| 我的老师免费观看完整版| 午夜亚洲福利在线播放| 国产在线精品亚洲第一网站| 一级a爱片免费观看的视频| 国产精华一区二区三区| 欧美日韩综合久久久久久| 乱人视频在线观看| 尾随美女入室| 天堂动漫精品| 人妻久久中文字幕网| 国产三级在线视频| 麻豆国产av国片精品| 内射极品少妇av片p| 亚洲国产欧洲综合997久久,| 又爽又黄无遮挡网站| 国产视频一区二区在线看| 免费大片18禁| 中国美白少妇内射xxxbb| 亚洲av美国av| 又爽又黄无遮挡网站| av中文乱码字幕在线| 亚洲va在线va天堂va国产| 午夜日韩欧美国产| 激情 狠狠 欧美| 老女人水多毛片| avwww免费| 少妇熟女欧美另类| 亚洲欧美精品自产自拍| 97热精品久久久久久| 91狼人影院| 亚洲av免费在线观看| 国产av一区在线观看免费| 中文字幕av成人在线电影| 亚洲七黄色美女视频| 真实男女啪啪啪动态图| 深夜a级毛片| 久久久精品欧美日韩精品| 国产大屁股一区二区在线视频| 日产精品乱码卡一卡2卡三| 免费搜索国产男女视频| 少妇熟女aⅴ在线视频| 亚洲综合色惰| 一级黄片播放器| 精品乱码久久久久久99久播| 午夜福利18| 精品久久久久久久久亚洲| 听说在线观看完整版免费高清| 韩国av在线不卡| 欧美不卡视频在线免费观看| 淫秽高清视频在线观看| 丝袜喷水一区| 亚洲乱码一区二区免费版| 一级毛片我不卡| 日韩在线高清观看一区二区三区| 日韩精品有码人妻一区| 少妇人妻精品综合一区二区 | 欧美日韩一区二区视频在线观看视频在线 | 亚洲五月天丁香| 国产色婷婷99| 搡女人真爽免费视频火全软件 | 日韩欧美国产在线观看| 国产高清激情床上av| 精品无人区乱码1区二区| 中文亚洲av片在线观看爽| 亚洲精品久久国产高清桃花| 欧美+日韩+精品| 18禁在线播放成人免费| 国产在线精品亚洲第一网站| 精品午夜福利在线看| 免费看日本二区| 亚洲成人精品中文字幕电影| 精品一区二区三区视频在线| 国产精品一区二区性色av| 欧美日韩乱码在线| 国产成人一区二区在线| 国产成人福利小说| 亚洲av中文字字幕乱码综合| 国产成年人精品一区二区| 免费看av在线观看网站| 亚洲欧美日韩无卡精品| 亚洲成人久久爱视频| 一区二区三区高清视频在线| 午夜老司机福利剧场| 国产激情偷乱视频一区二区| 欧美日韩在线观看h| 久久久色成人| 欧美一级a爱片免费观看看| 少妇丰满av| 麻豆精品久久久久久蜜桃| 精品一区二区三区视频在线观看免费| 午夜福利在线观看吧| 91狼人影院| av在线天堂中文字幕| 啦啦啦观看免费观看视频高清| 免费高清视频大片| 亚洲天堂国产精品一区在线| 国产精品一及| 色尼玛亚洲综合影院| 久久99热这里只有精品18| 国产单亲对白刺激| 国产高潮美女av| 亚洲人与动物交配视频| av.在线天堂| 亚洲最大成人手机在线| 亚洲最大成人av| 欧美潮喷喷水| 国产精品一区二区三区四区久久| 成人综合一区亚洲| 久久中文看片网| 黄色欧美视频在线观看| 亚洲欧美成人综合另类久久久 | 国产午夜精品久久久久久一区二区三区 | 99国产精品一区二区蜜桃av| av天堂中文字幕网| 99国产精品一区二区蜜桃av| 亚洲欧美日韩东京热| av视频在线观看入口| 97超碰精品成人国产| 国产男人的电影天堂91| 内地一区二区视频在线| 国产午夜精品久久久久久一区二区三区 | 国产v大片淫在线免费观看| 免费在线观看影片大全网站| 亚洲aⅴ乱码一区二区在线播放| 国产黄片美女视频| 精品免费久久久久久久清纯| 黄色一级大片看看| 精品免费久久久久久久清纯| 日韩成人伦理影院| 亚洲aⅴ乱码一区二区在线播放| 欧美日本视频| a级毛片a级免费在线| 久久亚洲国产成人精品v| 长腿黑丝高跟| 99热这里只有是精品50| 国产一区二区在线av高清观看| 国模一区二区三区四区视频| 国产白丝娇喘喷水9色精品| 日韩av在线大香蕉| 欧美激情久久久久久爽电影| 美女xxoo啪啪120秒动态图| 免费看光身美女| 亚洲最大成人av| 可以在线观看毛片的网站| 日韩三级伦理在线观看| 免费人成在线观看视频色| 日韩欧美一区二区三区在线观看| 麻豆乱淫一区二区| 日本色播在线视频| 免费看av在线观看网站| 欧美bdsm另类| 亚洲av不卡在线观看| 国产三级在线视频| 99热这里只有精品一区| 在线观看一区二区三区| 亚洲人成网站高清观看| 男女做爰动态图高潮gif福利片| 中文在线观看免费www的网站| 91av网一区二区| 别揉我奶头 嗯啊视频| 91av网一区二区| 亚洲第一电影网av| 成人国产麻豆网| 又黄又爽又刺激的免费视频.| 日韩欧美一区二区三区在线观看| 一级毛片aaaaaa免费看小| 国产精品人妻久久久影院| 午夜精品一区二区三区免费看| 国产真实伦视频高清在线观看| 国产人妻一区二区三区在| 久久久午夜欧美精品| 在线看三级毛片| 两个人的视频大全免费| 国模一区二区三区四区视频| 国产高潮美女av| 别揉我奶头~嗯~啊~动态视频| 日本五十路高清| 日韩 亚洲 欧美在线| 永久网站在线| 秋霞在线观看毛片| 亚洲av中文av极速乱| 最近2019中文字幕mv第一页| 无遮挡黄片免费观看| 亚洲精品456在线播放app| 久久精品国产清高在天天线| 九九爱精品视频在线观看| 黑人高潮一二区| 国产综合懂色| a级毛片免费高清观看在线播放| 国产免费男女视频| 欧美在线一区亚洲| 寂寞人妻少妇视频99o| 久久久色成人| 久久国产乱子免费精品| 性欧美人与动物交配| 国产大屁股一区二区在线视频| 草草在线视频免费看| 蜜桃久久精品国产亚洲av| 国产精品人妻久久久影院| 精品久久久久久久久久免费视频| 中文字幕熟女人妻在线| 特级一级黄色大片| 午夜免费激情av| 久久精品91蜜桃| 麻豆乱淫一区二区| 日本与韩国留学比较| 尤物成人国产欧美一区二区三区| 国产一区二区激情短视频| 日本熟妇午夜| 午夜精品在线福利| 成人三级黄色视频| 国产高清三级在线| 在线观看美女被高潮喷水网站| 波野结衣二区三区在线| 最后的刺客免费高清国语| 久久精品人妻少妇| 97超级碰碰碰精品色视频在线观看| 国产三级在线视频| 老师上课跳d突然被开到最大视频| 久久精品国产亚洲av香蕉五月| 精品国内亚洲2022精品成人| 嫩草影院新地址| 亚洲电影在线观看av| 日本欧美国产在线视频| 国产精品久久久久久久电影| 午夜精品在线福利| 午夜福利高清视频| 日本欧美国产在线视频| 1000部很黄的大片| 精品日产1卡2卡| 最新中文字幕久久久久| 大又大粗又爽又黄少妇毛片口| 天堂√8在线中文| 一个人看视频在线观看www免费| av在线播放精品| 国内精品美女久久久久久| av女优亚洲男人天堂| www日本黄色视频网| 如何舔出高潮| 国产国拍精品亚洲av在线观看| 特级一级黄色大片| 欧美成人a在线观看| 国内少妇人妻偷人精品xxx网站| 桃色一区二区三区在线观看| av国产免费在线观看| 精品久久久噜噜| 少妇的逼好多水| 久久鲁丝午夜福利片| 夜夜爽天天搞| 久久精品人妻少妇| 麻豆乱淫一区二区| 日韩欧美三级三区| 国产精品亚洲美女久久久| 岛国在线免费视频观看| 久久热精品热| 国产视频内射| 国产高清视频在线观看网站| 日日摸夜夜添夜夜添小说| 一区二区三区四区激情视频 | 大又大粗又爽又黄少妇毛片口| 嫩草影院新地址| 国产精品一区二区性色av| 麻豆精品久久久久久蜜桃| 免费不卡的大黄色大毛片视频在线观看 | 精品国内亚洲2022精品成人| а√天堂www在线а√下载| 亚洲专区国产一区二区| 亚洲人与动物交配视频| 亚洲成人精品中文字幕电影| 啦啦啦观看免费观看视频高清| 久久久精品94久久精品| 中文资源天堂在线| 亚洲精品国产av成人精品 | 亚洲最大成人av| 亚洲人成网站在线播放欧美日韩| 精品一区二区三区视频在线| 亚洲美女视频黄频| 国产 一区 欧美 日韩| 免费看美女性在线毛片视频| 国产精品1区2区在线观看.| 久久久久性生活片| 自拍偷自拍亚洲精品老妇| 99热这里只有是精品50| 一级毛片久久久久久久久女| 欧美不卡视频在线免费观看| 麻豆乱淫一区二区| 精品午夜福利在线看| 中文资源天堂在线| 最近在线观看免费完整版| 男女啪啪激烈高潮av片| 又爽又黄a免费视频| 国产黄色小视频在线观看| 变态另类成人亚洲欧美熟女| a级毛色黄片| 精华霜和精华液先用哪个| 啦啦啦观看免费观看视频高清| or卡值多少钱| 国产极品精品免费视频能看的| 国产中年淑女户外野战色| 色吧在线观看| 日韩一区二区视频免费看| 中国美女看黄片| 天天躁夜夜躁狠狠久久av| 99热这里只有是精品在线观看| 亚州av有码| 免费av观看视频| 最后的刺客免费高清国语| 午夜福利高清视频| 久久人妻av系列|