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

    Design and Control of Self-Heat Recuperative Distillation Process for Separation of Close-Boiling Mixtures: n-Butanol and iso-Butanol

    2015-06-22 14:19:10LiLuminLiuYuliangZhaiJianSunLanyiTianYuanyu
    中國煉油與石油化工 2015年4期

    Li Lumin; Liu Yuliang; Zhai Jian; Sun Lanyi; Tian Yuanyu

    (State Key Laboratory of Heavy Oil Processing, China University of Petroleum, Qingdao 266580)

    Design and Control of Self-Heat Recuperative Distillation Process for Separation of Close-Boiling Mixtures: n-Butanol and iso-Butanol

    Li Lumin; Liu Yuliang; Zhai Jian; Sun Lanyi; Tian Yuanyu

    (State Key Laboratory of Heavy Oil Processing, China University of Petroleum, Qingdao 266580)

    In order to explore the advantages of self-heat recuperative distillation (SHRD) process, the design and control of the SHRD process was studied for the separation of n-butanol and iso-butanol mixtures. The economic superiority of SHRD process is presented when a comparison on the total annual cost (TAC) of the conventional distillation process, the vapor recompression distillation process and the SHRD process was made. For the SHRD process, 37.74% and 11.35% savings of TAC can be achieved as compared to the conventional distillation process and vapor recompression distillation process, respectively. The dynamic characteristics of this promising SHRD sequence had been studied, and the dynamic responses demonstrated that 10% changes in both feed flow rate and feed composition can be well handled by the control strategy with dual-temperature control. It is proven that the SHRD system not only can provide economical savings but also can operate normally with good controllability.

    n-butanol; iso-butanol; self-heat recuperative distillation; total annual cost; control

    1 Introduction

    Distillation is one of unit operations widely used in petrochemical processes. However, distillation is also a wellknown process which requires a great deal of energy for the phase change and separation. Moreover, the CO2emissions are strongly related to energy consumption in distillation since the energy is mostly generated through the combustion of fossil fuels. Thus the energy savings of distillation process have been continuously attracting enormous interest from all over the world both in academic and industrial circles. Heat integration in distillation systems is an effective method to save energy, such as the dividing-wall column (DWC)[1-2], the heat pump assisted distillation[3-4]and the heat-integrated distillation column (HIDiC), etc.[5-6]Among these innovative methods, the heat pump assisted distillation process and the heat-integrated distillation process often use a compressor to increase the pressure and temperature of vapor stream, which after compression can be used to heat the liquid stream. However, the traditional heat recovery methods are adopted in these distillation processes without considering the preheating section. Recently, Kansha, et al. developed a self-heat recuperation technology which focused on the use of both latent heat and sensible heat.[7]In this attractive technology, the heat of process streams can be recycled to the process without any heat addition, which can bring about significant energy savings to the process. Then this attractive self-heat recuperation technology was applied to distillation process and the innovative self-heat recuperative distillation (SHRD) was proposed by Kansha, et al.[8]Figure 1(a) shows the flowchart of a SHRD process consisting of the preheating module and the distillation module, and in each module the total enthalpy of feed streams is equal to all enthalpy items of effluent counterparts. Figure 1(b) presents the temperature and heat relationship for the SHRD process, of which Tstdand Tbare the standard temperature and boiling point of stream 1, respectively. Each number in Figure 1(b) stands for the corresponding stream number of Figure 1(a). In this innovative SHRD process the vapor stream from the column overhead will be used as the heat source of feed stream and bottoms, which leads to significant energysavings. It is reported that the energy requirement of SHRD process can be decreased to 1/6—1/8 as compared with the traditional heat-exchanged distillation process[9]. Kansha, et al.[10]applied the self-heat recuperation technology to the flash distillation process, and the innovative process can achieve more than 75% reduction in the energy requirement as compared with the traditional process. On the basis of self-heat recuperation technology, Kansha and his co-workers then studied the HIDiC process, and the results showed that the developed distillation process was effective with further reduction of energy requirement[11]. Matsuda, et al.[12-13]also demonstrated that the distillation and petrochemical processes which applied the self-heat recuperation technology could achieve significant energy savings.

    In the carbonylation process for synthesizing n-butanol, the residues contain n-butanol, iso-butanol and some high boiling point components. However, the separation of these close-boiling mixtures consumes a great deal of energy. To reduce the energy requirements in the separation of n-butanol and iso-butanol mixtures, Gao, et al.[14]studied two types of mechanical vapor recompression heat pump distillation processes. The results indicated that the two processes were more economically favorable in both energy requirements and total annual cost (TAC). Nevertheless, the preheating section of feed stream has not been considered in the case of n-butanol/iso-butanol mixture separation.

    Though most of the previous works reported the design of distillation process that applied the self-heat recuperation technology, the related control system design of the SHRD process has not been fully investigated previously[15]. This work can be considered as an important supplement to the previous studies, with the focus being put on the application of SHRD and control issues. In this work, the SHRD process is used to separate the mixtures of n-butanol and iso-butanol for the first time, and then the TAC of the conventional distillation process, vapor recompression distillation column (VRC) process and SHRD process is calculated systematically. In addition, the dynamic controllability of the SHRD process is examined on the feed flow rate as well as the feed composition changes. Finally, the conclusions are drawn for evaluating the SHRD process.

    Figure 1 (a) Flowchart of SHRD process; (b) Temperature and heat relationship in SHRD process

    2 Steady-State Designs

    As an example of close-boiling mixture separation, the system of n-butanol (with a boiling point of 117.3 ℃ at 100 kPa) and iso-butanol (with a boiling point of 107.3 ℃ at 100 kPa) was studied. Gao, et al.[14]investigated this system without considering the preheating section. However, the preheating of feed is significant as the desired feed temperature is 114.2 ℃ for the separation. In this section, the conventional distillation process, VRC process and SHRD process will be explored for the separation of n-butanol and iso-butanol mixture. Figure 2 gives the conceptual design flowsheet of the three processes, and in each process the feed should be heated and the two products should be cooled. The design of these processes were implemented using the commercial simulator Aspen Plus, and the Wilson model was applied in this binary system[14].

    Figure 2 Design flowsheets of: (a) conventional distillation; (b) vapor recompression distillation; (c) self-heat recuperative distillation

    2.1 Simulations for the proposed three processes

    Figure 3 gives the steady-state flowsheet of the original one-column distillation. As shown in Figure 3, the feed flow rate is 5 000 kg/hr with the mixture composed of 55/45 n-butanol/iso-butanol, and the feed is heated before it is fed into the column. To obtain the minimum TAC, Gao, et al. investigated a range of operating pressures in their work and gave an optimal one of 90 kPa for the VRC process[14], and the operating pressure of 90 kPa was also assumed in this case. The specifications for the purity of two products (n-butanol and iso-butanol) were all set at 99.0%. To obtain a direct comparison, the feed and product conditions of VRC process and SHRD process should be consistent with those of the conventional process.

    Figure 4 shows the steady-state flowsheet of VRC process in detail. In the VRC process, the overhead stream of the column is compressed by a compressor to increase the temperature and pressure. Then the compressed stream is used as the heat source of the reboiler and is simultaneously condensed. As shown in Figure 4, the condensed stream is subsequently cooled by a trim condenser. Finally the overhead stream is divided into the distillate and reflux of the column.

    Figure 2(c) shows the flowsheet of the advanced SHRD process, in which the vapor distillate of the column is split into two parts. Subsequently, the divided streams are compressed by the compressors C1 and C2, respectively. The overhead product stream (compressed by C1) is cooled through exchanging heat with the feed stream. The latent heat of the other stream (compressed by C2)is treated as the reboiler duty, and then the condensed stream is recycled as the reflux to the column. The bottom stream of the column is divided, and the product stream (pure n-butanol) is used to preheat the feed stream, and the other stream functioning as the boilup is vaporized by the heat of compressed stream.

    Figure 3 Process flowsheet of conventional distillation system with details

    Figure 4 Process flowsheet of vapor recompression distillation system with details

    2.2 Economical optimization

    In this section, economical optimization is carried out to evaluate properly the economic feasibility of SHRD. TAC is used as the objective function, which can be described as Eq. (1),

    where the operating cost contains the steam cost for heater and reboiler, the cost of cooling water and the electricity cost for the compressor. The capital investment includes the cost of column trays, column shell, heat exchangers and compressors excluding the small equipment like pumps and valves because of their low costs. 3 years of capital payback period and 8 000 hours/year of available operating time in a plant are assumed. 0.852 kW/(K·m2) and 0.568 kW/(K·m2) are used as the heat transfer coefficients of condensers and reboilers, respectively[16]. The capital cost estimates were suggested by Douglas[17]with updated M&S index of 1 536.5 (Chemical Engineering Plant Cost Index, 4 th Q, 2011).

    The minimum TAC exists when a series of total number of trays (NT) varies under a certain feed location. To achieve the optimum parameters, the TAC of the SHRD process was calculated, and Figure 5 presents graphically the TAC of the SHRD process via changing the number of stages. It can be seen from Figure 5(a) that when the number of stages increases, the operating cost changes slightly,while the capital cost decreases firstly when NT is less than 67, and then it increases when NT is equal to or greater than 67. The minimum TAC of 682.17×103$/year occurs when NT is 67 as revealed in Figure 5(b).

    Figure 5 Effects of the number of stages on (a) capital cost and operating cost; (b) TAC of the SHRD process

    To compare the three processes fairly, the conventional process and VRC process are also optimized in the same way as described in SHRD process, and the detailed simulation results are given in Table 1. As shown in Table 1, the VRC process provides a 29.77% reduction in TAC as compared with the conventional process. Moreover, the SHRD configuration can achieve 11.35% of TAC savings, as compared to the VRC design. Hence the SHRD design is the most economical one when it is compared with the conventional design and VRC design.

    Table 1 Optimum results of the SHRD design, the conventional distillation and VRC distillation designs

    3 Control Structure Design

    Because the SHRD process has a greater advantage in TAC than the other two designs, it would also be interesting to examine the dynamic characteristics and controllability of this advanced process. Some dynamic parameters of the SHRD process should be specified before the simulation is exported to Aspen dynamics as the pressuredriven simulation. The diameter of SHRD column is calculated through the tray sizing function in Aspen Plus. The size of column base is specified based on a 5-min liquid holdup time with a half full level.

    The overhead vapor flow rate is manipulated to control the pressure at the top of column, and the bottoms flow rate is used to control the base level of the column. After the pressure and base level are controlled, some disturbances, which include the throughput changes and feed composition changes, are introduced to investigate the stability of the proposed SHRD design. Figure 6 displays the responses of product purity under 10% feed changes. As shown in Figure 6, the feed disturbance in each case would always result in the purity reduction of one product. As the stage temperature is not controlled to meet the separation requirement, when feed flow rate increases, there is not enough energy to heat the feed stream, which will make more iso-butanol get into the distillation column bottom. On the contrary, the distillate will contain more n-butanol when the throughput decreases. In the case of feed composition disturbances, the changes on iso-butanol will have a direct impact on the product purity. For example, a +10% iso-butanol feed composition disturbance will lead to an increase in iso-butanol product purity and accordingly a decrease in n-butanol product purity.

    As shown in Figure 6, though the SHRD system can achieve stable operation after various feed disturbances are introduced, the products cannot keep at their desired purity level. Therefore the control of the SHRD design will be studied afterwards. The temperature profile along the optimum SHRD column is presented in Figure 7, and the relative large slopes can be observed on the stages 20 and 50. Thus, according to the slope criterion the stages 22 and 50 will be considered as the temperature control points later. Besides the control over column pressure and base level, the feed flowrate must be controlled, and the flow rate of iso-butanol product is used to control the feed temperature. In addition, the reflux rate is manipulated to control the temperature of the stage 22 via adjusting the brake power of the compressor C2. The temperature controllers with a 1-min deadtime are tuned based on relay-feedback tests by the Tyreus-Luyben tuning rule[18]. These proposed controllers used in the loops are conventional PI controllers, and the basic control structure is presented in Figure 8.

    Figure 9 shows the dynamic responses to ±10 % changesin feed stream for this basic control strategy. It is shown that despite some basic controllers used in the common cases, the two products cannot maintain the desired purity after the feed disturbances are introduced, which shows the similar trends as presented in Figure 6. Symmetrical responses with relatively large deviations can be observed in all temperature responses of the stage 22. However, the temperature responses of the stage 50 have not fluctuated similarly as the tendency of the stage 22 (see Figure 9). When the feed flow rate ranges from 5 000 kg/hr to 5 500 kg/hr, the tray temperature especially on the stage 50 shows a large decrease, leading to a great reduction in n-butanol product purity. The effect on the stage 50 is quite small when the feed flowrate is reduced by 10%, resulting in a small offset of n-butanol product purity. For an increase of 10% isobutanol in the feed composition, the tray temperature ofthe stage 50 decreases significantly with large deviations, which will affect the purity of n-butanol product.

    Figure 6 Dynamic responses when the disturbances are: (a) +10% in feed flow rate; (b) -10% in feed flow rate; (c) +10% in iso-butanol composition; (d) -10% in iso-butanol composition

    Figure 7 Column temperature proflle for the optimum SHRD process

    Figure 8 Basic control structure of the SHRD process

    Figure 9 Dynamic responses of the SHRD process for the basic control structure with ±10% disturbances in (a) feed flow rate; (b) feed composition

    To overcome the defects of the basic control structure, this low relative volatility system requires the use of dualtemperature control. An improved control structure for the SHRD design is proposed after a lot of control structures have been investigated following continuously explored efforts. Based on the original control structure, another temperature controller TC3 is added to control the tray temperature of the stage 50 as shown in Figure 10, which is adjusted by the ratio of boilup rate to feed flow rate (B/F). The sequential iterative tuning procedure is used to tune the temperature controllers[19].

    Figure 10 Improved control structure of the SHRD process

    Figure 11 reveals the simulation results of the SHRD process with the improved control structure for ±10% feed disturbances. As shown in Figure 11, the temperature control point of the stage 22 is able to quickly return back to its setpoint. Although the settling time is relatively long when the SHRD system is adapted to a 10% increase in iso-butanol content of the feed composition, the temperature responses of the stage 50 are all able to hold fairly close to their setpoints, which leads to the improvement in n-butanol quality. Judging from the closed-loop responses, it can be seen that the improved control strategy with the B/F ratio can handle different feed disturbances effectively. Moreover, the SHRD system with the improved control scheme can achieve stability faster than with the basic one. Therefore, the improved control structure is suitable for this SHRD process.

    4 Conclusions

    The self-heat recuperative distillation is an attractive technology used in the separation process. In this work, the conventional distillation process, VRC process and SHRD process are investigated for the separation of n-butanol and iso-butanol mixtures. The optimum design is obtained by minimizing the TAC, and the results show that, for the SHRD process, 37.74% and 11.35% reductions in TAC can be obtained compared with the conventional distillation process and VRC process, respectively.

    The dynamic characteristics of this SHRD system are investigated, and two control structures are proposed for the optimum SHRD design. The dynamic responses reveal that the SHRD process with the basic control structure can barely bring the product purity to their setpoints. Therefore, in order to control the system effectively, an improved control structure with dual-temperature control is proposed. It is shown that with the improved control structure, the SHRD system can handle various feed disturbances effectively, and the high-purity products can be obtained.

    Acknowledgements:Financial supports from the National Natural Science Foundation of China (Grant: 21276279 and Grant: 21476261) and the Fundamental Research Funds for the Central Universities (No. 14CX05030A; No. 15CX06042A) are acknowledged with gratitude.

    [1] Halvorsen I J, Dejanovi? I, Skogestad S, et al. Internal configurations for a multi-product dividing wall column[J]. Chemical Engineering Research and Design, 2013, 91(10): 1954-1965

    [2] Wu Y C, Lee H Y, Huang, H. P, et al. Energy-saving dividing-wall column design and control for heterogeneous azeotropic distillation systems[J]. Industrial & Engineering Chemistry Research, 2014, 53(4): 1537-1552

    [3] Kiran B, Jana A K. Introducing vapor recompression mechanism in heat-integrated distillation column: impact of internal energy driven intermediate and bottom reboiler[J]. AIChE Journal, 2015, 61(1): 118-131

    [4] Waheed M A, Oni A O, Adejuyigbe S B, et al. Performance enhancement of vapor recompression heat pump[J]. Applied Energy, 2014, 114: 69-79

    Figure 11 Dynamic responses of the SHRD process for the improved control structure with ±10% disturbances in: (a) feed flow rate; (b) feed composition

    [5] Suphanit B. Optimal heat distribution in the internally heatintegrated distillation column (HIDiC)[J]. Energy, 2011, 36(7): 4171-4181

    [6] Bruinsma O S L, Krikken T, Cot J, et al. The structured heat integrated distillation column[J]. Chemical Engineering Research and Design, 2012, 90(4): 458-470

    [7] Kansha Y, Tsuru N, Sato K, et al. Self-heat recuperation technology for energy saving in chemical processes[J]. Industrial & Engineering Chemistry Research, 2009, 48(16): 7682-7686

    [8] Kansha Y, Tsuru N, Fushimi C, et al. New design methodology based on self-heat recuperation for production by azeotropic distillation[J]. Energy & Fuels, 2010, 24(11): 6099-6102

    [9] Kansha Y, Tsuru N, Fushimi C, et al. Integrated process module for distillation processes based on self-heat recuperation technology[J]. Journal of Chemical Engineering of Japan, 2010, 43(6): 502-507

    [10] Kansha Y, Tsuru N, Fushimi C, et al. An innovative modularity of heat circulation for fractional distillation[J]. Chemical Engineering Science, 2010, 65(1): 330-334

    [11] Kansha Y, Kishimoto A, Tsutsumi A. Process design methodology for high-energy saving HIDiC based on self-heat recuperation[J]. Asia-Pacific Journal of Chemical Engineering, 2011, 6(3): 320-326

    [12] Matsuda K, Kawazuishi K, Hirochi Y, et al. Advanced energy saving in the reaction section of the hydro-desulfurization process with self-heat recuperation technology[J]. Applied Thermal Engineering, 2010, 30(16): 2300-2305

    [13] Matsuda K, Kawazuishi K, Kansha Y, et al. Advanced energy saving in distillation process with self-heat recuperation technology[J]. Energy, 2011, 36(8): 4640-4645

    [14] Gao X, Chen J, Ma Z, et al. Simulation and optimization of distillation processes for separating a close-boiling mixture of n-butanol and isobutanol[J]. Industrial & Engineering Chemistry Research, 2014, 53(37): 14440-14445

    [15] Kansha Y, Kishimoto A, Tsutsumi A. Analysis of dynamic characteristics for self-heat recuperative process[J]. Chemical Engineering Transactions, 2012, 29: 307-312

    [16] Sun L, Wang Q, Li L, et al. Design and control of extractive dividing wall column for separating benzene/cyclohexane mixtures[J]. Industrial & Engineering Chemistry Research, 2014, 53(19): 8120-8131

    [17] Douglas J M. Conceptual Design of Chemical Processes[M]. New York: McGraw-Hill: 1988

    [18] Luyben W L. Plantwide Dynamic Simulators in Chemical Processing and Control[M]. New York: Marcel Dekker: 2002

    [19] Huang H P, Jeng J C, Chiang C H, et al. A direct method for multi-loop PI/PID controller design[J]. Journal of Process Control, 2003, 13(8): 769-786

    Trial Production of Needle Coke at Shanghai Petrochemical Company

    Needle coke is the starting material for manufacture of premium graphite electrode, which at present cannot be produced by SINOPEC and even other domestic refining enterprises so that the carbon producing plants inside China have to import a significant amount of high-quality needle coke every year. The project “Production of premium needle coke and development of engineering techniques” jointly undertaken by the SINOPEC Research Institute of Petroleum Processing and the SINOPEC Shanghai Petrochemical Company Limited (SPCL) after completion of the preparative stages involving distillation of slurry oil and hydrodesulfurization of slurry fractions on July 10, 2015 began to enter into the commercial-scale coking test phase. Till July 18, 2015, the test team has performed four operating cycles, and the first stage of trial coking had manufactured a total of 1800 tons of coke, the properties of which are being investigated. The said commercial tests were carried out in the 1.0 Mt/a delayed coking unit at SPCL. In the course of commercial tests there were a lot of difficulties related with control over the operation of equipment because of the drastic changes in production process conditions, and any malfunction could result in tower flooding or coking inside the heating furnace. The technical personnel of RIPP have been working cooperatively to formulate emergency response measures while making directives and instructions to secure the safe and smooth operation of the process unit.

    date: 2015-08-03; Accepted date: 2015-09-09.

    Prof. Sun Lanyi, Telephone: +86-13854208340; Fax: +86-532-86981787; E-mail: sunlanyi@upc.edu.cn.

    国产免费一级a男人的天堂| 亚洲精品国产av蜜桃| 老熟女久久久| 黄色配什么色好看| 亚洲精品,欧美精品| 夫妻午夜视频| 曰老女人黄片| 久热这里只有精品99| 国产男女超爽视频在线观看| 午夜av观看不卡| 国产亚洲最大av| 亚洲精品国产av成人精品| 三上悠亚av全集在线观看| 中国美白少妇内射xxxbb| 精品国产露脸久久av麻豆| 一级毛片黄色毛片免费观看视频| 免费观看a级毛片全部| av播播在线观看一区| av一本久久久久| 下体分泌物呈黄色| 国产日韩欧美在线精品| 永久免费av网站大全| 一本色道久久久久久精品综合| 国精品久久久久久国模美| 热re99久久精品国产66热6| 午夜视频国产福利| 亚洲激情五月婷婷啪啪| 精品久久国产蜜桃| 少妇 在线观看| 久久这里有精品视频免费| 高清黄色对白视频在线免费看| 爱豆传媒免费全集在线观看| 亚洲成人手机| 黄色视频在线播放观看不卡| 日产精品乱码卡一卡2卡三| 午夜福利在线观看免费完整高清在| 亚洲欧美一区二区三区黑人 | 热99国产精品久久久久久7| 精品亚洲成a人片在线观看| 国产免费现黄频在线看| 亚洲综合色惰| 热99国产精品久久久久久7| 亚洲av成人精品一区久久| 最近中文字幕2019免费版| 国产成人91sexporn| 毛片一级片免费看久久久久| 国产亚洲av片在线观看秒播厂| 国产精品一区二区在线观看99| 国产国语露脸激情在线看| 亚洲欧美日韩卡通动漫| 国产深夜福利视频在线观看| 精品熟女少妇av免费看| 你懂的网址亚洲精品在线观看| 又大又黄又爽视频免费| 丰满少妇做爰视频| 伊人久久国产一区二区| 国产熟女欧美一区二区| 亚洲av福利一区| 在现免费观看毛片| 久久这里有精品视频免费| 22中文网久久字幕| 久久久久人妻精品一区果冻| 热re99久久精品国产66热6| 高清视频免费观看一区二区| 搡女人真爽免费视频火全软件| 热re99久久国产66热| 久久久久久人妻| 黄色一级大片看看| 日本wwww免费看| 亚洲婷婷狠狠爱综合网| tube8黄色片| 日本免费在线观看一区| 国产男女内射视频| 国产成人aa在线观看| 在线观看免费视频网站a站| 母亲3免费完整高清在线观看 | 亚洲精品自拍成人| 视频中文字幕在线观看| 久热久热在线精品观看| 成人午夜精彩视频在线观看| 欧美老熟妇乱子伦牲交| 欧美成人午夜免费资源| 男人添女人高潮全过程视频| 色94色欧美一区二区| 亚洲精品视频女| 最近最新中文字幕免费大全7| 一边亲一边摸免费视频| 日韩成人av中文字幕在线观看| 美女内射精品一级片tv| 熟女人妻精品中文字幕| 高清毛片免费看| 精品人妻在线不人妻| 女人久久www免费人成看片| 天天操日日干夜夜撸| 欧美精品一区二区免费开放| 国产精品国产三级国产av玫瑰| 免费大片18禁| 免费观看的影片在线观看| 99热6这里只有精品| 欧美成人午夜免费资源| 伦精品一区二区三区| 国产精品无大码| 高清不卡的av网站| 久久久久久久久久久丰满| 国产精品久久久久久久电影| 大码成人一级视频| 狂野欧美激情性bbbbbb| 亚洲国产日韩一区二区| 国产欧美另类精品又又久久亚洲欧美| 一边摸一边做爽爽视频免费| 久久久久久久久久久免费av| 午夜久久久在线观看| 一级黄片播放器| 99国产精品免费福利视频| 精品国产乱码久久久久久小说| 99九九在线精品视频| 久久久久久久亚洲中文字幕| 欧美亚洲 丝袜 人妻 在线| 亚洲精品久久午夜乱码| 国产亚洲精品第一综合不卡 | 街头女战士在线观看网站| 久久精品久久精品一区二区三区| 少妇精品久久久久久久| 久久久午夜欧美精品| 国产成人av激情在线播放 | 免费高清在线观看视频在线观看| 国产av国产精品国产| 亚洲av在线观看美女高潮| av.在线天堂| 亚洲内射少妇av| 午夜激情福利司机影院| 久久精品国产亚洲av涩爱| 日韩精品有码人妻一区| 亚洲欧美中文字幕日韩二区| 午夜福利,免费看| 国产一区二区在线观看日韩| 一本一本综合久久| 高清视频免费观看一区二区| 亚州av有码| 亚洲人与动物交配视频| 亚洲经典国产精华液单| 欧美国产精品一级二级三级| 欧美成人午夜免费资源| 亚洲精品乱码久久久久久按摩| 综合色丁香网| 国产成人精品一,二区| 久久免费观看电影| 亚洲精品av麻豆狂野| videos熟女内射| 母亲3免费完整高清在线观看 | 26uuu在线亚洲综合色| 午夜福利影视在线免费观看| 久久久精品免费免费高清| 国产精品成人在线| 亚洲综合精品二区| 人妻一区二区av| 国产成人91sexporn| 下体分泌物呈黄色| 久久精品人人爽人人爽视色| 午夜久久久在线观看| 婷婷色av中文字幕| 欧美xxⅹ黑人| 91在线精品国自产拍蜜月| 国产日韩欧美亚洲二区| av播播在线观看一区| videossex国产| 人妻 亚洲 视频| 亚洲熟女精品中文字幕| 免费日韩欧美在线观看| 色5月婷婷丁香| 亚洲av成人精品一二三区| 在线 av 中文字幕| 久久久久国产网址| 国产欧美另类精品又又久久亚洲欧美| 中文字幕亚洲精品专区| 高清午夜精品一区二区三区| 国产综合精华液| 亚洲精品一区蜜桃| 汤姆久久久久久久影院中文字幕| 最近手机中文字幕大全| 国产一区二区三区综合在线观看 | 国产精品嫩草影院av在线观看| 国产成人午夜福利电影在线观看| 日日撸夜夜添| 国产精品久久久久久精品古装| 国产av精品麻豆| 国国产精品蜜臀av免费| 亚洲少妇的诱惑av| 亚洲国产精品国产精品| 精品人妻偷拍中文字幕| 精品卡一卡二卡四卡免费| 秋霞伦理黄片| 亚洲av.av天堂| 亚洲欧美清纯卡通| 久久国产精品男人的天堂亚洲 | 视频区图区小说| 亚洲婷婷狠狠爱综合网| videos熟女内射| 九九在线视频观看精品| 少妇的逼好多水| 特大巨黑吊av在线直播| 日韩人妻高清精品专区| 久久 成人 亚洲| 美女中出高潮动态图| 亚洲精品日韩av片在线观看| 亚洲人与动物交配视频| 狠狠婷婷综合久久久久久88av| 久久影院123| 国产在线一区二区三区精| 亚洲一区二区三区欧美精品| 一本一本综合久久| 波野结衣二区三区在线| 久久人人爽av亚洲精品天堂| 精品久久久噜噜| 亚洲av中文av极速乱| 99久久综合免费| 99久久精品一区二区三区| 又大又黄又爽视频免费| 丝袜脚勾引网站| 国产淫语在线视频| 男的添女的下面高潮视频| 成人国产麻豆网| 九九在线视频观看精品| 国产一区二区在线观看日韩| 日日摸夜夜添夜夜爱| 亚洲国产精品国产精品| 纵有疾风起免费观看全集完整版| 美女xxoo啪啪120秒动态图| 18禁动态无遮挡网站| 欧美日韩综合久久久久久| 一本色道久久久久久精品综合| 99久国产av精品国产电影| 中文字幕久久专区| 日韩一本色道免费dvd| 久久97久久精品| 午夜影院在线不卡| 人人妻人人澡人人爽人人夜夜| 大码成人一级视频| 国产日韩欧美视频二区| 精品一品国产午夜福利视频| av播播在线观看一区| 成人无遮挡网站| 国产成人精品一,二区| 18在线观看网站| 久久 成人 亚洲| 一级爰片在线观看| 这个男人来自地球电影免费观看 | 热re99久久国产66热| 热re99久久精品国产66热6| 精品国产露脸久久av麻豆| 久久久精品94久久精品| 妹子高潮喷水视频| 午夜视频国产福利| 大香蕉久久成人网| 亚洲国产最新在线播放| 在线观看人妻少妇| 国产黄色免费在线视频| 在线观看美女被高潮喷水网站| 啦啦啦啦在线视频资源| videosex国产| 一级毛片 在线播放| 男人爽女人下面视频在线观看| 国产精品人妻久久久久久| 欧美人与性动交α欧美精品济南到 | 多毛熟女@视频| 国产亚洲一区二区精品| 精品亚洲成国产av| 亚洲欧美精品自产自拍| 国产国拍精品亚洲av在线观看| 少妇人妻精品综合一区二区| 亚洲欧美成人精品一区二区| 日本免费在线观看一区| xxx大片免费视频| 黄色毛片三级朝国网站| 简卡轻食公司| 三上悠亚av全集在线观看| 久热这里只有精品99| 国产成人精品福利久久| 国产成人精品无人区| 69精品国产乱码久久久| 又粗又硬又长又爽又黄的视频| 久久久久久久精品精品| 热99国产精品久久久久久7| 特大巨黑吊av在线直播| 一级毛片aaaaaa免费看小| 国产精品.久久久| 亚洲av国产av综合av卡| 中国三级夫妇交换| 国产综合精华液| 曰老女人黄片| 乱人伦中国视频| 麻豆成人av视频| 交换朋友夫妻互换小说| 又大又黄又爽视频免费| 91久久精品国产一区二区三区| 蜜桃国产av成人99| 久久久久久久久久人人人人人人| 亚洲精品自拍成人| 各种免费的搞黄视频| 日日爽夜夜爽网站| 少妇熟女欧美另类| tube8黄色片| 国产在线免费精品| 午夜日本视频在线| 各种免费的搞黄视频| 久久97久久精品| 秋霞在线观看毛片| 亚洲第一av免费看| 性高湖久久久久久久久免费观看| 一级,二级,三级黄色视频| 最近最新中文字幕免费大全7| 免费观看性生交大片5| 亚洲欧美中文字幕日韩二区| 极品少妇高潮喷水抽搐| 不卡视频在线观看欧美| 亚洲av国产av综合av卡| 中文精品一卡2卡3卡4更新| 日本vs欧美在线观看视频| 人人妻人人添人人爽欧美一区卜| 国产精品嫩草影院av在线观看| 国产深夜福利视频在线观看| 人人妻人人爽人人添夜夜欢视频| 边亲边吃奶的免费视频| 99热网站在线观看| 日日摸夜夜添夜夜爱| 热re99久久精品国产66热6| 一区二区日韩欧美中文字幕 | 九九爱精品视频在线观看| 色5月婷婷丁香| 91午夜精品亚洲一区二区三区| 爱豆传媒免费全集在线观看| 大片免费播放器 马上看| 欧美亚洲日本最大视频资源| 国产午夜精品一二区理论片| 国产免费福利视频在线观看| 久久久久久伊人网av| 亚洲精品av麻豆狂野| 国产69精品久久久久777片| 人妻 亚洲 视频| 永久免费av网站大全| 久久久久视频综合| av国产久精品久网站免费入址| 欧美成人精品欧美一级黄| 亚洲国产精品999| 亚洲精品乱码久久久久久按摩| 自拍欧美九色日韩亚洲蝌蚪91| 国语对白做爰xxxⅹ性视频网站| 久久久久久人妻| 久久精品国产亚洲av涩爱| 免费人妻精品一区二区三区视频| 有码 亚洲区| 女性被躁到高潮视频| 精品亚洲成a人片在线观看| 国产免费现黄频在线看| 少妇被粗大猛烈的视频| 国产成人aa在线观看| .国产精品久久| 日韩三级伦理在线观看| 亚洲国产最新在线播放| av在线app专区| 大陆偷拍与自拍| 搡女人真爽免费视频火全软件| 国产av国产精品国产| 欧美成人精品欧美一级黄| 欧美激情国产日韩精品一区| 汤姆久久久久久久影院中文字幕| 夫妻午夜视频| 日韩 亚洲 欧美在线| 日韩一区二区视频免费看| 人成视频在线观看免费观看| 又粗又硬又长又爽又黄的视频| 老司机亚洲免费影院| 国产色爽女视频免费观看| 在线天堂最新版资源| 成人国语在线视频| 插阴视频在线观看视频| 久久久国产欧美日韩av| 99精国产麻豆久久婷婷| 久久精品国产自在天天线| 亚洲色图综合在线观看| 如何舔出高潮| 久久鲁丝午夜福利片| 99热国产这里只有精品6| 少妇猛男粗大的猛烈进出视频| 午夜精品国产一区二区电影| 大香蕉久久网| 亚洲怡红院男人天堂| 最后的刺客免费高清国语| 91成人精品电影| 夜夜看夜夜爽夜夜摸| 精品午夜福利在线看| 男人爽女人下面视频在线观看| 成年女人在线观看亚洲视频| a级毛片在线看网站| 久久ye,这里只有精品| 国产精品国产三级国产专区5o| 尾随美女入室| 国产熟女午夜一区二区三区 | 97超碰精品成人国产| 亚洲精品,欧美精品| 色视频在线一区二区三区| 欧美 日韩 精品 国产| 久久久久久久久久久免费av| www.色视频.com| 高清欧美精品videossex| 国产免费又黄又爽又色| 22中文网久久字幕| 日韩强制内射视频| 欧美3d第一页| 精品酒店卫生间| 久久久久久伊人网av| 老司机影院成人| 日韩欧美一区视频在线观看| 蜜臀久久99精品久久宅男| 久久久久久久精品精品| 九九爱精品视频在线观看| 在线天堂最新版资源| 2018国产大陆天天弄谢| 久久久久国产精品人妻一区二区| a级片在线免费高清观看视频| 亚洲熟女精品中文字幕| 在线观看免费高清a一片| 国产亚洲精品第一综合不卡 | 亚洲第一区二区三区不卡| 夫妻午夜视频| 亚洲精品456在线播放app| 亚洲人成77777在线视频| 亚洲国产欧美在线一区| 狂野欧美激情性bbbbbb| 99re6热这里在线精品视频| 国产精品秋霞免费鲁丝片| 亚洲av免费高清在线观看| 两个人免费观看高清视频| 日本av免费视频播放| av国产久精品久网站免费入址| 精品一区二区三卡| 成人手机av| av国产久精品久网站免费入址| 大又大粗又爽又黄少妇毛片口| 超色免费av| 美女内射精品一级片tv| 久久狼人影院| 国产午夜精品久久久久久一区二区三区| 亚洲色图综合在线观看| 一个人免费看片子| 天天躁夜夜躁狠狠久久av| 久久青草综合色| 丝袜喷水一区| 午夜日本视频在线| 五月天丁香电影| 日本免费在线观看一区| 99久久综合免费| .国产精品久久| 99久久综合免费| 国产成人精品在线电影| 日日摸夜夜添夜夜爱| 精品酒店卫生间| 久久久精品94久久精品| 一区二区日韩欧美中文字幕 | 国产精品无大码| 亚洲av综合色区一区| 满18在线观看网站| 久久国产亚洲av麻豆专区| 亚洲国产日韩一区二区| 日日摸夜夜添夜夜爱| 一个人看视频在线观看www免费| 我的女老师完整版在线观看| 黑人巨大精品欧美一区二区蜜桃 | 亚洲精品国产色婷婷电影| 亚洲少妇的诱惑av| 日韩中文字幕视频在线看片| 国产成人午夜福利电影在线观看| 成人国语在线视频| 午夜精品国产一区二区电影| 成人黄色视频免费在线看| 久久久久精品性色| 天堂中文最新版在线下载| 一区二区三区乱码不卡18| 久久这里有精品视频免费| 你懂的网址亚洲精品在线观看| 精品一区在线观看国产| 日韩一区二区视频免费看| 一级a做视频免费观看| 日韩电影二区| 国产亚洲欧美精品永久| 丝袜喷水一区| 成人国语在线视频| 国产在线一区二区三区精| 一区二区av电影网| 伊人久久国产一区二区| 日日撸夜夜添| 日本欧美视频一区| 五月玫瑰六月丁香| tube8黄色片| 十分钟在线观看高清视频www| 在线亚洲精品国产二区图片欧美 | 亚洲av福利一区| 精品99又大又爽又粗少妇毛片| 国产一区二区在线观看日韩| 亚洲成人一二三区av| 美女大奶头黄色视频| 一级毛片我不卡| 大陆偷拍与自拍| 亚洲av免费高清在线观看| 插逼视频在线观看| 国产老妇伦熟女老妇高清| 大话2 男鬼变身卡| 久久久久国产精品人妻一区二区| 亚洲欧美一区二区三区黑人 | 欧美日韩视频精品一区| 三级国产精品片| 亚洲精品日韩在线中文字幕| 亚洲精品成人av观看孕妇| 久久99蜜桃精品久久| 丝袜在线中文字幕| av福利片在线| 精品亚洲成a人片在线观看| 精品午夜福利在线看| 最近的中文字幕免费完整| 欧美日本中文国产一区发布| 国产精品女同一区二区软件| 久久久久久久亚洲中文字幕| 国产欧美另类精品又又久久亚洲欧美| 亚洲第一av免费看| 一级a做视频免费观看| 国产午夜精品一二区理论片| 久久午夜综合久久蜜桃| 午夜福利,免费看| xxxhd国产人妻xxx| 女人久久www免费人成看片| 国语对白做爰xxxⅹ性视频网站| 夫妻性生交免费视频一级片| 亚洲综合精品二区| 成人亚洲精品一区在线观看| 日韩欧美精品免费久久| 免费av不卡在线播放| 日韩欧美精品免费久久| 天堂中文最新版在线下载| a 毛片基地| freevideosex欧美| 国产成人aa在线观看| 国产精品偷伦视频观看了| 国产爽快片一区二区三区| 十分钟在线观看高清视频www| 涩涩av久久男人的天堂| 精品久久久久久电影网| 国产亚洲精品第一综合不卡 | 国产伦精品一区二区三区视频9| 国产精品麻豆人妻色哟哟久久| 亚洲内射少妇av| 欧美老熟妇乱子伦牲交| 91久久精品国产一区二区成人| 亚洲人成77777在线视频| 热99国产精品久久久久久7| 久久久午夜欧美精品| 在线亚洲精品国产二区图片欧美 | 欧美成人午夜免费资源| xxx大片免费视频| 精品亚洲成国产av| 有码 亚洲区| 午夜免费男女啪啪视频观看| 观看美女的网站| 国产白丝娇喘喷水9色精品| 99久久精品国产国产毛片| 成人毛片60女人毛片免费| 国产乱来视频区| 日本av免费视频播放| 精品99又大又爽又粗少妇毛片| www.色视频.com| tube8黄色片| 自拍欧美九色日韩亚洲蝌蚪91| 蜜桃在线观看..| 少妇猛男粗大的猛烈进出视频| 国产av码专区亚洲av| 久久久久国产精品人妻一区二区| 日产精品乱码卡一卡2卡三| 精品99又大又爽又粗少妇毛片| 亚洲精品色激情综合| 免费看不卡的av| 国产高清国产精品国产三级| videosex国产| 老司机影院成人| 久久99蜜桃精品久久| 日日啪夜夜爽| 多毛熟女@视频| 十八禁网站网址无遮挡| 亚洲高清免费不卡视频| 欧美精品人与动牲交sv欧美| 在线观看www视频免费| 日本vs欧美在线观看视频| 日韩欧美一区视频在线观看| 国产黄片视频在线免费观看| 99视频精品全部免费 在线| 简卡轻食公司| 久久99蜜桃精品久久| 丝瓜视频免费看黄片| 久久久久视频综合| 在线观看免费视频网站a站| 久久久久久久精品精品| 国产在线免费精品| av卡一久久| 国产成人精品福利久久| 色婷婷久久久亚洲欧美| 亚洲色图 男人天堂 中文字幕 | 久久99热6这里只有精品| 日韩精品有码人妻一区| 美女中出高潮动态图| 久久青草综合色| 国产av码专区亚洲av| 九色亚洲精品在线播放| 免费观看a级毛片全部| 少妇精品久久久久久久| 爱豆传媒免费全集在线观看| 18+在线观看网站| 国产av码专区亚洲av|