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

    Integration of energy recovery network including recycling residual pressure energy with pinch technology

    2017-05-28 08:50:16JianqiangDengZhengCaoDongboZhangXiaoFeng

    Jianqiang Deng *,Zheng Cao Dongbo Zhang ,Xiao Feng

    1 School of Chemical Engineering and Technology,Xi'an Jiaotong University,Xi'an 710049,China

    2 Shaanxi Coal and Chemical Technology Institute Co.,Ltd,Xi'an 710070,China

    1.Introduction

    In recent years,energy-saving in the energy and chemical industry has gone through a stage of process integration,i.e.,the entire system is deemed as an organic whole and the work focus is the global optimum design.The two most common forms of energy are heat and pressure.At the beginning of integration technology development,the recovery and integration of heat energy is generally accepted as a key consideration.Many integration technologies,such as mathematical programming technologies and pinch technology,have been proposed for optimizing the heat exchanger network(HEN).In particular,pinch technology is a simple and efficient method for process optimization,and presents a clear overall view of energy consumption efficiency[1].In recent decades,pinch technology has led to the recovery of large amounts of heat energy through the proper design of HEN,and it has also been applied in hydrogen network and water network optimizations,benefiting the entire process with minimum utility consumptions and cost-savings[2,3].

    In addition to thermal power plants,compression or boosting is a major consumer of energy in many chemical plants,including those producing Lique fied Natural Gas(LNG),refineries,air enrichment and ammonia,etc.Some industrial streams need work for compression,while others can export work through expansion.Traditionally,the work produced or required can be converted into electric energy for furtherusage or integration.Hanet al.[4]developed an online optimization system and applied it to the condensing steamturbine network of a chemical plant.An optimization design was formulated by utilizing a developed model to maximize the total electric power recovery from the steam turbine network.Aspelundet al.[5]proposed a graphical heuristic methodology based on extending traditional pinch analysis with exergy calculations to minimize energy requirements(total shaft work)in sub-ambient processes.Similarly,Shinet al.[6]proposed a Mixed Integer Linear Programming(MILP)formulation for optimizing boil-off gas compressor operations to minimize the total average power consumption in an LNG receiving and re-gasification terminal.El-Halwagiet al.[7]introduced a targeting pinch diagram for the integration of heat and power generation process,and the optimal steam flows and process target of cogeneration were rapidly identified.Del Nogalet al.[8,9]reported a model able to obtain an optimal power plant system for utility networks.Appropriate gas/steam turbines as drivers were matched with compressor stages by using mathematical programming methods,and the developed methodology was applied to a LNG case study.

    In fluid process systems,whenever a high-pressure fluid stream(stream HP)is depressurized,a potential for pressure energy recovery exists.However,pressurizing a low-pressure fluid stream(stream LP)always requires compression work.Optimal integration or matching of the required booster and expansion devices may lead to significant savings in energy.Chenget al.[10]proposed the concept of a work exchanger(WE)in 1960,and described the basic working process of a positive displacement-type work exchanger(PDWE).Based on the working principle,work exchangers can be classified as PDWE or turbine-type(centrifugal type)WE.The PDWE is an isobaric energy recovery device characterized by its working principle.It can be classified as piston type pressure exchanger and rotary pressure exchanger.The former uses a pair of stationary cylinder where pistons reciprocate regularly associated with valve motion.The latter uses a cylindrical rotor with a plurality of open-ended ducts,and the rotor spins between two end covers to realize a continuous pressure exchange process[11].Because the hydraulic energy is directly transferred between high pressure stream and feed stream,the efficiency of PDWE is higher than that of turbine-type WE[12].Investigations show that PDWEs have an energy transfer efficiency of 92%–97%[13–15],while turbine-type WEs have an efficiency of 72%–75%[13,16].Currently,work exchange is successfully applied in worldwide seawater reverse osmosis(SWRO)desalination process,where the PDWE contributes to as much as 60%of energy savings in the RO process[17].An RO desalination demonstration project operated in Changdao and in Shengsi suggests the PDWE in such plants can decrease the energy consumption to 4.0 kW·h·m-3for freshwater production[18].

    Based on the similarity of the heat and mass transfer process,Huang and Fan[19]proposed the work exchange network(WEN)in 1996.By applying the basic principles of energy conservation,necessary and sufficient conditions were derived to identify feasible stream matching in work exchangers.Afterwards,Razibet al.[20]modeled the details of compressor and turbine operating curves to identify high-pressure and low-pressure streams that could be matched for work exchange via a turbine-type WE.Considering the characteristics of the PDWE,Liuet al.[21]applied a graphical integration method in the streammatching process.This allowed for the maximum energy recovery and optimal network structure to be simultaneously identified.Pressure energy exists not only in liquid streams but also in gas streams.Denget al.[22]analyzed the performance of gas–gas work exchangers based on the energy conservation law and the real-gas law.It was proposed that a heat-work coupling transfer network could be a more flexible way to recover energy in gas streams.Onishiet al.[23]presented a superstructure for HEN synthesis by adjusting pressure levels of process streams to improve heat integration.The mathematical model was developed using generalized disjunctive programming and reformulated in mixed-integer nonlinear programming.It indicated that the integration of heat exchange and pressure manipulation equipment could decrease the total annualized cost.

    Energy integration on a large scale is the current trend in technology development.However,as indicated by the above investigations,there are very few studies focused on WENs,and especially on the heat work coupling transfer network.Currently,pinch technology has been successfully applied in HENs with simple,visualized and flexible characteristics.Therefore,a systematic approach and a detailed scheme including recycling residual pressure energy with pinch technology are proposed to construct a heat-work coupling transfer network.Because of its high energy transfer efficiency,the PDWE has been selected to build a WEN in this work.For comparison,how to build a WEN based on turbine-type WEs is also discussed.Related issues,such as heat transfer in the device,starting mode,etc.,are analyzed in a practical work transfer process on a coupling network.Finally,as a case study,a typical rectisol process in a coal-water slurry gasification section of synthetic ammonia is chosen to demonstrate the detailed synthesis process.

    2.Problem Statement

    Original streams are a set of liquid streams with different temperatures and pressures.Since it is possible to transfer both heat and work among these original streams,the optimized matching between each stream constitutes a complicated heat-work coupling transfer network.To develop a reasonable matching scheme based on efficient and logical procedures has become the main content of this integration design.The aim of the integration design is to make effective use of heat energy and pressure energy from original streams.As this is an optimization problem,the given parameters of each original stream are fluid properties such as composition,density,viscosity and thermal conductivity,and operation parameters such as the flow rate,pressure and temperature.Besides the required parameters set in the heat exchange process,corresponding parameters in the pressure exchange process have also been previously set.This includes the pressure drop of each work exchanger,the minimum driving pressure differential(ΔPdr),the lowest economical pressure drop(ΔPmin)and the minimum economical indicated work.Relevant variables include the flow rate,pressure and temperature of each stream in the new heat-work network.The optimization target is to find the matching scheme with the highest energy recovery rate.

    3.Systematic Approach and Procedure

    In process design,pinch analysis is applied as an efficient tool to guarantee minimum energy levels in design of HEN.Since work and heat have similar transfer mechanisms and analogous exchange network designs,it is therefore possible to apply pinch technology to the design and optimization of WENs.By analogy with the temperature pinch point in an HEN,the place with the minimum pressure difference in a WEN is defined as the pressure pinch point.Similarly,the pressure pinch point can be obtained from the composite curve of a simple WEN,while the problem table method is more efficient for complex WENs with more streams.Once the pinch pressure is obtained,the matching rules and restrictions involved by pinch position helps to identify the various matching schemes and to select the most appropriate one for integration.The basic systematic approach on a heat-work coupling transfer network can be described as follows:An HEN is constructed first,which commonly determines the main effect of energy recovery.This step can be neglected when reforming an existing exchange network because its HEN is usually already made.It is then necessary to analyze basic data and build a WEN based on pinch technology.Since the WEN always alters the state parameters of the previous HEN to some extent,the final heat-work coupling transfer network can be achieved only by adjusting the operation parameters of the network.What needs to be noted is that there can be different WEN schemes established by using different types of work exchangers.Fig.1 shows the systematic approach to building a heat-work coupling transfer network using pinch analysis,with a detailed explanation below.

    Fig.1.Scheme and steps to construct a heat-work coupling transfer network.

    Step 1:Construction of an HEN

    Heat recovery has played a very important role in energy-saving technology.The working procedures of an HEN optimal design are relatively mature.There are several reasons for constructing an HEN in advance:

    (1)Recovered heat is commonly an important part of the total recovered energy in an energy or chemical system.Thus,it is necessary to analyze the heat recovery first and build an HEN in advance.

    (2)Stream pressure shows little change in the heat transfer process.When an HEN has been constructed,the pressure distribution of the system remains relatively stable,and this provides maximal recovery potential of the pressure energy for a WEN.

    (3)After the HEN design process,all system node temperatures can be deemed as known data.A work exchanger can then be located at the lower temperature position of a stream,and material requirements can be reduced accordingly.This also provides a relatively benign operational environment when operating a WEN.

    (4)The temperature variation of a stream flowing through the PDWE is very small.A newly added WEN has a small impact on the previous HEN.

    (5)Because an HEN usually exists in a practical energy or chemical process,the systematic approach used in this study can be directly used in a renovation project.

    Step 2:Data extraction

    The basic data mainly includes pressure,temperature and flow rate.According to changes in pressure and flow rate during the process,the streams as optimized objects of the WEN are picked out.Among these,the streams providing pressure energy are defined as work sources,while the streams requiring work are defined as work sinks.

    Since liquid streams are assumed to be incompressible,the volume change and temperature effect on the work exchange process are neglected.In this case,the amount of indicated workWindcan be calculated based on the pressure difference ΔPbetween original pressure and target pressure[19],

    whereVis the volumetric flow rate of a stream.

    Step 3:Initial screening

    In practical applications,there are more restrictions in the work transfer process than that in the heat transfer process:

    (1)The original pressure of the work source(PHP,o)must be greater than the target pressure of the work sink(PLP,t),and the target pressure of the work source(PHP,t)must be less than the original pressure of the work sink(PLP,o).This restriction exists mainly in the PDWE,such as for DWEER and PX commercial products.As shown in Fig.2,the in flow streams in the cylinder can be discharged and the transfer process can be accomplished only whenPHP,ois higher thanPLP,tandPHP,tis lower thanPLP,oin each cycle.In this study,the minimum pressure difference required to push the piston in one piston stroke of PDWE is 0.1 MPa[10,20].

    Fig.2.Pressurization process in a cylinder.

    (2)The work transfer process between two streams must occur in the same transfer device or in two interconnected devices with the same shaft.The distance between the two streams is limited to a certain range.This restriction is particularly important for an actual engineering process.

    (3)According to pinch technology,when the driving force near the pinch is tiny,the increased number of work exchangers will not significantly improve the energy recovery efficiency.Based on this,the lowest economical pressure drop near the pinch point ΔPminis set to 0.2 MPa.In addition,a stream has a certain amount of recoverable pressure energy when its flow rate is large,regardless of the pressure level.A WE is deemed economically feasible to be set when a stream with the indicated work is greater than 5 kW.If a stream has little exchangeable pressure energy,it cannot serve as the exchanging object when considering equipment investment and operational costs.

    In view of the above-mentioned restrictions,this work proposes a two-stage screening method for choosing streams.Initial screening is carried out when the pinch point is unknown.According to restriction(3),the streams with small exchangeable pressure energy and the completely un fit pressure ranges should be ruled out,and marked as “N”.A Turbine-type WE can be used in broad work conditions without considering restriction(1).Then streams that only use the turbine-type WE for pressure energy recovery are selected and labeled as “T”.The rest of the streams can be used as work exchange objects for the PDWE,and temporarily marked as“Y”.By singling out the streams that are suitable for the PDWE matching scheme,the matching workload in subsequent steps is decreased.

    The values of work exchanger efficiency and pump efficiency are based on those most commonly found in the literature.The efficiency of a centrifugal pump with a large flow rate( flow rate no less than 50 m3·h?1)is 80%,while that of a smaller pump( flow rate less than 50 m3·h?1)is 65%[24,25].The energy transfer efficiency of a turbinetype WE is 75%[16].Based on a large number of applications data for SWRO,it is reasonable to set the energy transfer efficiency of a PDWE to 94%.As PDWEs are proposed as preferred equipment in this work,it is preferable to build a WEN containing PDWEs.

    Step 4a:Determination of pinch position

    In order to construct a WEN scheme based on PDWEs,the pressure pinch point is needed.In reference to the problem table method of the HEN,the pressure pinch point of the WEN can be obtained by the following procedure:

    Step 4.1:Divide pressure intervals;

    After decreasing work source pressure by half of ΔPmin,and increasing work sink pressure by half of ΔPmin,the average pressure of each stream is obtained.Then by arranging all average pressures in descending order,the pressure region is divided into several intervals.

    Step 4.2:Calculate work balance of each interval;

    The delivered work values of sinks and sources are set to positive and negative respectively,which are added up to calculate the de ficit work of each interval.

    Step 4.3:Calculate cumulative work of each interval without work input;

    The output work of each interval is calculated as the difference between input work and deficit work,and is equal to the input work of the next interval.

    Step 4.4:Determine minimum work input of utility engineering;

    When the output work is less than zero,external work needs to be introduced to ensure the sufficiency of work in each interval.Therefore,the minimum work in put is set to the opposite value of the minimum of the surplus work output.

    Step 4.5:Determine pressure pinch point.

    Step 4b:Matching scheme based on turbine-type WEs

    There are no pressure limits in the turbine-type WE,and even a high-pressure fluid can be pressurized by a low-pressure fluid.In some cases,the PDWE cannot be used in a transfer network and turbine is the only appropriate device,although energy transfer efficiency of turbine equipment is relatively low.

    Planning a WEN scheme based on turbines does not require consideration of the pinch point.The primary concern is to match the amount of energy of the work sources and work sinks to maximize the quantity of recovered work.In addition,there is an upper limit to the matching number of streams when network stability is taken into account.The larger the number of other streams matching with a certain stream,the more susceptible the system is to be instable.In this work,the upper limit of the matching number is set as 2,which means one stream matches no more than two other streams.After establishing a WEN scheme based on turbines,we can proceed directly to Step 8.

    Step 5:Second screening

    The second screening is carried out when the pressure pinch point is obtained.Key streams that supply or need more indicated work are primarily investigated.If the minimum pressurize utility is positive in pinch analysis,the biggest work sources are the key streams.If the minimum pressurize utility is negative,the biggest work sinks are the key streams.Based on the restriction that streams do not match across the pinch,the target pressure of work source across the pinch should be set to pinch pressure.Considering the work exchange efficiency,the amount of practical recovered workWprcan be calculated by multiplyingWindwith energy transfer efficiency η,

    In addition,some streams must be boosted in advance to the pressure that is specified by rule(1)of Step 3 before matching,which is termed the pre-boost process.“Y”can then be subdivided into three situations,namely “N”,“M”,and “P”“N”represents the stream that is unsuitable for matching after identification in the second screening,such as the stream with too little recovered work after the pre-boost process.“P”refers to the stream that can be directly applied by the PDWE.“M”indicates the stream that can be applied by the PDWE after pre-boosting.

    Step 6:WEN matching scheme based on PDWEs

    As shown in Fig.3,work sources are listed on the left and work sinks are on the right,whereas the work exchangers are installed between them.The streams above the dotted line are above the pinch pressure.Due to the matching restrictions that work exchange is not applied across the pinch,the implementation of pinch point can simplify the superstructure of a WEN.Similar to an HEN,there are five principles that should be considered when applying pinch technology to build a WEN based on PDWEs,especially to match key streams:

    Fig.3.Superstructure of a WEN with pressure pinch point.

    (1)Number of streams.In theory,pressure-relief devices should not be set above the pinch.This means that all HP streams achieve de pressurization by transferring work to LP streams,i.e.,NHP≤NLP.Here,NHPis the number of work sources with pressure higher than the pinch.NLPis the number of work sinks above the pinch.Similarly,the booster should not be set below the pinch point,i.e.,NHP≥NLPat the moment.

    (2)Volumetric flow rate.In order to ensure that the minimum pressure difference is at the pinch point,the streams around the pinch should meet the following rules:above the pinch,VHP≥VLP;below the pinch,VHP≤VLP.

    The two principles(1)and(2)play guiding roles only near the pinch,and they are not strict requirements.

    (3)Maximize exchange load.In order to ensure the minimum number of work exchange units,each matching process should not be ended until the energy of one stream is used up.

    (4)Sequence rule.Streams are matched in the order of descending power supplied by work sources.

    (5)Matching numbers.One stream matches no more than two other streams.

    After the key streams are matched,the remaining streams can be operated manually.Then a WEN based on PDWEs is generated.

    Step 7:Construction of the final WEN

    For streams which are weak and not suitable for the PDWE,the turbine-type WE is used to recover the rest pressure energy.Combining the characteristics of turbines and the PDWE,a final WEN can be constructed with the aim of achieving the maximum recovered work.This also can be adjusted based on equipment investment and production cost.

    Step 8:Synthesis of the heat-work coupling transfer network

    Considering the temperature variations of streams caused by heat transfer or throttling that occurs during the work exchange process,the operating conditions of the initial HEN need to be adjusted.The optimized heat-work coupling transfer network can then be formed.

    4.Discussion Regarding the Application Problem

    The two working fluids involving pressure energy transfer can be in liquid–liquid form,liquid–gas form,gas–liquid form or gas–gas form.Among these,the liquid–liquid form has the highest energy transfer efficiency(near 1.0)[26].Moreover,work transfer in the liquid–liquid form abounds in a practical process.Therefore,this section and the subsequent case study mainly focus on application problems in the liquid–liquid form.

    4.1.Equipment selection

    The turbine-type WE is simple,runs smoothly,has few working condition limitations and is already applied in practical projects.However,the turbine-type WE has a lower energy transfer efficiency.In comparison,the PDWE has a higher energy transfer efficiency,even though it has a complex structure and more working condition constraints.Moreover,the successful application of PDWE in SWRO desalination indicates the possibility that the device could be similarly applied in the energy or chemical industry.Therefore,in terms of the higher efficiency and energy savings,a WEN matching scheme based on PDWE is the best choice.In specific cases,a hybrid matching scheme based on PDWE and the turbine-type WE,and even a pure matching scheme based on the turbine-type WE,is an alternative choice.

    4.2.Throttling effects

    For the PDWE,high pressure streams are de pressurized transiently at the beginning of their flow in the cylinder,and the resulting phase change and cooling effects persist for an extremely short period of time.Therefore,no phase change is assumed in the recovery process,and the throttling effect can be ignored when establishing and optimizing a WEN because it has a little effect on the system.

    4.3.Interaction analysis

    4.3.1.Heat transfer effect

    Since HP streams and LP streams with different temperatures flow through the cylinder alternately in the transfer process,the heat transfer phenomena will inevitably occur between the HP stream and the cylinder as well as between the cylinder and the LP stream.In addition,the pushing and mixing between two streams in the cylinder also leads to temperature variation of the streams.Fortunately,the amount of mixing is small according to the existing service experience of the PDWE in the SWRO system.As a result,the influence of temperature variation of the two streams will spread to the entire HEN.Therefore,the whole network must be adjusted.Considering that the heat transfer process in the PDWE is similar to a regenerative heat exchanger(or regenerator),this work uses the heat-transfer model of a regenerator[27]to approximately calculate the temperature change of the streams.

    4.3.2.Pressure resistance of a HE

    The maximum operating pressure of a heat exchanger in a heatwork coupling transfer network can be reduced by means of the work transfer process,or more specifically,by only adjusting its position from the front of a WE to the back.By adjusting the position,the pressures of the two streams in an HE can be closer to each other.The closer the two streampressures are in the heat transfer process,the higher the resulting pressure resistance capacity of the heat exchanger.As a result,the material grade of the heat exchanger can be reduced.

    4.4.Starting mode

    Some mature experiences in the processing industry can be referenced to analyze the starting mode of the WEN.A pre-configured pump provides the original driving force to push a key work-sink stream.The power of the pump should meet the work demand of the stream to pressurize from original pressure to target pressure,and the priority selection of the driving system is a variable frequency motor.

    The starting mode of the process can beforecasted as firstly building liquid circulation to make streams flow in an over all system.The pumps will start corresponding to key work sinks in order to boost the front part of the system.The work exchangers driven by the pressure differential between the front and the back of the system then begin to work.Meanwhile,the related gas streams and work-source streams should be under control when the system pressure rises.When the system pressure reaches the normal state,all the work exchangers will start and the power of pre-configured pumps will simultaneously decrease.Then the system enters normal operating conditions.

    5.Case Study

    The rectisol process involves a gas purification technology that is widely used in the coal chemical industry.This process uses methanol as the sorption solvent and removes acidic gas in crude syngas using the selective absorption characteristics of a low temperature methanol solution.In this study,we analyze a typical rectisol process in a constructed synthetic ammonia plant with an annual output of 300 k tons.Since the HEN already exists,the design of the WEN and the optimization of the heat-work coupling transfer network in this system start from Step 2.

    Step 2:Data extraction

    Fig.4 shows the basic process of the typical five-tower rectisol.It can be observed that the high-pressure methanol stream comes into contact with the process gas in the absorption tower(1)and absorbs acid gas like CO2and H2S during the washing process.After that,one part of the stream flows out at the bottom of the absorption tower(1)and enters into the separator(S)after being depressurized to 1.75 MPa(G)by a pressure-reducing valve.Another part of the stream is extracted at the lower section of the tower,and is depressurized to 1.75 MPa(G)by a pressure-reducing valve.The volumetric flow rates of the two branch streams whose pressures significantly change are both larger than 100 m3·h?1and can be deemed as two main work sources(SR1,SR2).The other two work sources(SR3,SR5)are the two liquid streams separated by the separator(S).The two streams both are depressurized to approximately 0.2 MPa(G)by pressure-reducing valves.The volumetric flow rates of these are relatively large,and a considerable amount of pressure energy can be recovered.

    As shown in Fig.4,the gas generated by the gas–liquid separator and the nitrogen-rich hydrogen gas coming from the liquid nitrogen wash mix together form the first work sink(SK1).Corresponding to the two pumps used to pressurize the extracting solutions from the knockout drum(KD)and the concentration tower(3),two work sinks(SK2,SK3)exist in the cycling system between the absorption tower(1)and the CO2desorption tower(2).The methanol circulating pump corresponding to the largest work sink(SK5)is the largest power consumer,and the methanol solution is pressurized by the pump from 0.05 MPa(G)to 6.4 MPa(G)before entering the absorption tower(1).There are two other work sinks(SK6,SK7)at the bottom and top of the methanol regeneration tower(4),which are used for the tower's liquid transportation and the top re flux,respectively.The basic data of the work sources and work sinks are shown in Table 1,where SR represents a work source and SK represents a work sink.

    Step 3:Initial screening

    Since this work mainly discusses the liquid–liquid matching process,and PDWE is not expected to be suitable for gas streams,SK1 is not taken into accountas a work sink.According to stipulations of the initial screening,SR6,SK2,and SK7 are also not suitable for matching because their pressure energy is too small.Therefore,these four streams are marked as“N”,and the rest of the streams that can be further investigated are marked as“Y”.The initial screening results are displayed in column 7 in Table 1.

    Step 4a:Determination of pinch position

    As stated in Section 3,for the PDWE,the energy transfer efficiency is 0.94 and the pinch point pressure difference ΔPminis set as 0.2 MPa.Based on the specific data analysis of the SR and SK with the problem table method used in pinch technology,the entire system is first divided into 13 intervals as shown in Fig.5.

    Fig.4.Schematic diagram of typical five-tower rectisol process.1—Absorption tower;2—CO2 Desorption tower;3—Concentration tower;4—Methanol Regeneration tower;5—Methanol–Water Separator;C—Compressor;H—Heater;KD—knockout drum;S—Separator.

    Deficit work refers to the work consumption of each interval,and a negative deficit work indicates that the pressure interval has surplus work.The computed deficit work of respective pressure intervals are shown in column 3 of Table 2.After calculating the work balance of each interval,the cumulative work of each pressure interval without input work from the outside are shown in columns 4 and 5 of Table 2.To make the cumulative output work of each pressure interval nonnegative,external work with a minimum of 155.34 kW is introduced.Finally,the calculated results of recounting the output work of each pressure interval are shown in the last two columns of Table 2.Since all the pressure intervals' output work is not less than zero and the eleventh interval output work is zero,the pinch pressure is therefore determined to be 0.15 MPa.

    Step 4b:Matching scheme based on turbine-type WEs

    As previously stated,application of the turbine-type WE is not restricted by pressure range.A WEN based on the turbine-type WE is built mainly as a schematic comparison.According to the matching principle mentioned in Step 4b of Section 3,a turbine-type WEN matching scheme can be obtained as shown in Fig.6.In the diagram,“T”represents a turbine-type work exchanger,while the gray solid circle represents supercharging equipment such as a pump.The efficiency of a turbine-type WE is 0.75,and the amount of recovered work can be calculated based on the basic data in Table 1.The calculated values of recovered work are listed in Table 3.

    Fig.5.Pressure interval division of a specific rectisol process.

    Table 1Basic data and screening results

    Step 5:Second screening

    According to pinch technology,the minimum pressurizing utility is positive(155.34 kW).The two work-source streams with the largestindicated works(SR1 and SR2)are chosen as the key streams.

    A second screening is carried out based on the pressure limit(PHP,o>PLP,t,PLP,o>PHP,t)and pinch point value(0.15 MPa).The pressures of key streams SR1 and SR2 are above the pinch,which are also between the original pressure and the target pressure of SK5(Table 1).Therefore,SR1 and SR2 can be marked as “P”in the last column of Table 1.Since the original pressures of SR3 and SR5 are above the pinch,and the target pressure of SK4 or SK6 is within the pressure ranges of SR3 and SR5,SR3 and SR5 can be marked as“M”when they are matched with SK4 or SK6 with the aid of a certain supporting process.Similarly,the pressure range of SR4 is above the pinch,and it can be matched with SK4 or SK5 through a supporting process.Therefore,SR4 can be marked as “M”.It is found that the pressure interval of SK4 is across the pinch.Because the target pressure meets the matching requirements of pressure and volume flow,SK4 can be marked as“M”by pre-boosting the rest part below the pinch.The pressure range of stream SK5 is the largest(0.05 MPa–6.4 MPa)among the SK streams,and it can be matched with most of the SR streams.Therefore,SK5 is marked as “M”.The original pressure of SK3 is 0.11 MPa which cannot be matched unless it is pressurized to 0.35 MPa in advance.This is because the target pressure of SR streams upon the pinch should be 0.25 MPa,and the driving pressure differential(0.1 MPa)is also required.The remaining energy ofSK3 is only 2.7 kW,and SK3 is therefore marked as “N”.Similarly,the actual energy that can be transferred for SK6 is only 3.47 kW,and SK6 is also marked as“N”.The specific screening results are shown in the last column of Table 1.

    Table 2Computation sheet of problem table method

    Table 3Recovered work from the turbine-type WE matching scheme

    Step 6:WEN matching scheme based on PDWEs

    The matching scheme is divided into two parts that are above and below the pinch.Matching of the both parts proceeds respectively.According to the volumetric flow rate rule(2),the two streams cannot satisfy the volume flow requirement above and below the pinch simultaneously.Because the pinch pressure is small(only 0.15 MPa)without considering of the driving force requirement(0.1 MPa),the matching scheme below the pinch will not be considered in this case.

    Based on the results of the second screening and its five rules,two key streams(SR1 and SR2)can be matched as follows:SR1 can be matched with the branch of SK5(the corresponding flow is 113 m3·h?1,the original pressure and target pressure are 1.95 MPa and 5.05 MPa,respectively);SR2 can be matched with the other part of SK5(the corresponding flow is 110m3·h?1,the original and target pressures are 1.95 MPa and 5.05 MPa,respectively).

    On the basis ofthe above results,there are 4 streams left that need to be matched:SR3,SR4,SR5 and SK4.Among the work sources,the work supplied by SR5 is closest to the energy needed by SK4,so SR5 is chosen to match with SK4.The matching scheme based on PDWEs is shown in Fig.7.In the diagram,“M”indicates a PDWE unit with a pre-boosting process.However,it is observed that the direct matching,namely “P”,is not used widely and does not occur in this case.A gray solid circle in Fig.7 indicates supercharging equipment such as a pump,and a horizontal triangle indicates a pressure-reducing valve.According to the basic data and the energy transfer efficiency of the related equipment type,the amount of recovered work in each exchange unit can be obtained,and the specific values are shown in Table 4.

    Fig.6.Matching scheme based on turbine-type WEs.

    Fig.7.Matching scheme based on PDWEs.

    Table 4Recovered work of the PDWE matching scheme

    Step 7:Construction of the final WEN

    Considering the matching scheme based on the PDWE and characteristics of the turbine-type WE,the remaining SR streams that are not matched can be utilized by using pure turbines(T1,T2,T3 in Table 5)to recover their pressure to generate electricity.The initial WEN can then be achieved,and the required configuration is shown in Fig.8.The specific quantity of recovered work of the initial WEN can be calculated.These results are shown in Table 5.

    Step 8:Synthesis of the heat-work coupling transfer network

    The incoming WEN will change the parameters of the original HEN.To obtain the final heat-work coupling exchange network,the influence of temperature changes caused by the WEN must be considered.Table 6 shows the original pressure,temperature and flow rate data of the streams that pass through the three PDWE exchange units.

    The heat-transfer calculation model in PDWEs is similar to that of the regenerator mentioned in Section 3.The detailed description of the equations is omitted here.The difference in heat transfer temperature differences can be calculated iteratively.The specific calculation results are shown in Tables 7,8 and 9.

    From Table 7,itcan be observed that the most reasonable heat transfer temperature difference(Δtm)of unit 1 is approximately 50.0 °C,and the temperature changes of the two streams are both close to 7.8°C.

    From Table 8,it can be noted that the optimal heat transfer temperature difference(Δtm)of unit 2 is approximately 60.0 °C,and the temperature changes of the two streams are both 9.8°C.The optimaltmof unit 3 is approximately 18.0°C,and the stream temperature changes are nearly 2.2°C,as shown in Table 9.

    As the stream temperature changes are elucidated,we can adjust the operating conditions of the HEN and then generate the final heat-work coupling transfer network.

    Table 5Specific value of recovered work of the initial WEN

    6.Conclusions

    Global energy integration of an energy system has become more sophisticated and is not only a combination of heat energy,but also an integration of pressure energy.The specific integration method to building a heat-work coupling transfer network has been presented in this study.A systematic approach favors an HEN to be built first before a WEN is added to the energy transfer network,based on pinch technology.The final heat-work coupling transfer network is generated through fine-tuning the operation parameters of entire network.With the highest efficiency of energy transfer,the positive displacement type work exchanger is given priority to build a work exchange network.The study results show that pinch technology can be successfully used in a heat-work coupling transfer network.The influence of a positive displacement work exchanger on the energy transfer network in a coupling transfer system was investigated.The results showed that the effect of the phase change,throttle cooling and heat transfer could be neglected in the synthesis process.The theoretical analysis indicated that the approach was conceptually feasible.Finally,a case study of the typical rectisol process and the effect of added energy recovery is presented.The operability of the proposed method on the comprehensive energy transfer network is verified.This work provides a reference for constructing a more complex and larger energy synthesis network.

    Fig.8.Configuration of the initial WEN.

    Table 6Basic data of the three exchange units

    Table 7Temperature change of work exchanger unit 1

    Table 9Temperature change of unit 3

    Nomenclature

    Fflow rate,m3·h?1

    Nnumber of streams above the pinch

    Ppressure,MPa

    Qquantity of heat transfer in a period per cycle,kJ

    Vvolume,m3

    Wwork,kW

    Ttemperature,°C

    Δtmtemperature difference(set value),°C

    Δtm′ temperature difference(calculated value),°C

    η energy transfer efficiency

    Subscripts

    HEN heat exchanger network

    HP high pressure

    HT high temperature

    ind indicated work

    LP low pressure

    LT low temperature

    o original

    PDWE positive displacement-type work exchanger

    pr practical recovered work

    t target

    WE work exchanger

    WEN work exchange network

    Acknowledgments

    This work was supported by the National Natural Science Foundation of China(No.20936004 and No.21376187).

    [1]M.Ebrahim,A.Kawari,Pinch technology:An efficient tool for chemical-plantenergy and capital-cost saving,Appl.Energy65(1)(2000)45–49.

    [2]G.L.Liu,H.Li,X.Feng,C.Deng,Pinch location of the hydrogen network with purification reuse,Chin.J.Chem.Eng.21(12)(2013)1332–1340.

    [3]Z.Y.Liu,Y.M.Li,G.L.Zhang,Y.Z.Yang,Simultaneously designing and targeting for networks with multiple resources of different qualities,Chin.J.Chem.Eng.17(3)(2009)445–453.

    [4]I.S.Han,Y.H.Lee,C.H.Han,Modeling and optimization of the condensing steam turbine network of a chemical plant,Ind.Eng.Chem.Res.45(2)(2006)670–680.

    [5]A.Aspelund,D.O.Berstad,T.Gundersen,An extended pinch analysis and design procedure utilizing pressure based exergy for subambient cooling,Appl.Therm.Eng.27(11)(2007)2633–2649.

    [6]M.W.Shin,D.Shin,S.H.Choi,E.S.Yoon,C.G.Han,Optimization of the operation of boil-off gas compressors at a liquefied natural gas gasification plant,Ind.Eng.Chem.Res.46(20)(2007)6540–6545.

    [7]M.El-Halwagi,D.Harell,H.D.Spriggs,Targeting cogeneration and waste utilization through process integration,Appl.Energy86(6)(2009)880–887.

    [8]F.L.Del Nogal,J.K.Kim,S.Perry,R.Smith,Synthesis of mechanical driver and power generation configurations,Part 1:Optimization framework,AIChE J.56(9)(2010)2356–2376.

    [9]F.L.Del Nogal,J.K.Kim,S.Perry,R.Smith,Synthesis of mechanical driver and power generation configurations,Part 2:LNG applications,AIChE J.56(9)(2010)2377–2389.

    [10]C.Y.Cheng,S.W.Cheng,L.T.Fan,Flow work exchanger,AIChE J.13(3)(1967)438–442.

    [11]O.M.Al-Hawaj,The design aspects of rotary work exchanger for SWRO,Desalin.Water Treat.8(1-3)(2009)131–138.

    [12]J.X.Sun,Y.Wang,S.C.Xu,S.C.Wang,Energy recovery device with a fluid switcher for seawater reverse osmosis system,Chin.J.Chem.Eng.16(2)(2008)329–332.

    [13]O.V.Sallangos,E.Kantilaftis,Operating experience of the Dhekelia seawater desalination plant using an innovative energy recovery system,Desalination173(1)(2004)91–102.

    [14]W.T.Andrews,W.F.Pergande,G.S.McTaggart,Energy performance enhancements of a 950 m3/d seawater reverse osmosis unit in Grand Cayman,Desalination135(1)(2001)195–204.

    [15]B.Pe?ate,J.A.De la Fuente,M.Barreto,Operation of the RO Kinetic?energy recovery system:Description and real experiences,Desalination252(1)(2010)179–185.

    [16]J.Y.Mak,D.Heaven,Syngas purification in gasification-based ammonia/urea plants,Ammonia Plant Saf.Relat.Facil.39(1999)319–331.

    [17]R.L.Stover,Seawater reverse osmosis with isobaric energy recovery devices,Desalination203(1)(2007)168–175.

    [18]Y.W.Tan,B.Tan,Q.Wang,Progress in desalination project in China,Technol.Water Treat.33(1)(2007)1–3(in Chinese).

    [19]Y.L.Huang,L.T.Fan,Analysis of a work exchanger network,Ind.Eng.Chem.Res.35(10)(1996)3528–3538.

    [20]M.S.Razib,M.M.F.Hasan,I.A.Karimi,Preliminary synthesis of work exchange networks,Comput.Chem.Eng.37(2012)262–277.

    [21]G.L.Liu,H.Zhou,R.J.Shen,X.Feng,A graphical method for integrating work exchange network,Appl.Energy114(2014)588–599.

    [22]J.Q.Deng,J.Q.Shi,Z.X.Zhang,X.Feng,Thermodynamic analysis on work transfer process of two gas streams,Ind.Eng.Chem.Res.49(24)(2010)12496–12502.

    [23]V.C.Onishi,M.A.Ravagnani,J.A.Caballero,Simultaneous synthesis of heat exchanger networks with pressure recovery:Optimal integration between heat and work,AIChE J.60(3)(2014)893–908.

    [24]Y.D.Zozulyak,Y.L.Levashov,V.T.Novatsky,Precision setup for fluid pump efficiency testing,Heat Transf.Res.28(7)(1997)499–502.

    [25]I.B.Cameron,R.B.Clemente,SWRO with ERI's PXpressure exchanger device-a global survey,Desalination221(1)(2008)136–142.

    [26]J.Q.Deng,L.Ma,X.Feng,Work recovery efficiency of work exchange in chemical process,J.Chem.Ind.Eng.62(11)(2011)3017–3023.

    [27]H.T.Cui,F.Yang,Thermal-storage Technology and Application,Chemical Industry Press,Beijing,2004(in Chinese).

    美女脱内裤让男人舔精品视频| 亚洲综合色网址| 久久精品人人爽人人爽视色| 久久精品熟女亚洲av麻豆精品| 久久精品夜色国产| 最近手机中文字幕大全| 亚洲色图综合在线观看| 久久综合国产亚洲精品| 2022亚洲国产成人精品| 午夜精品国产一区二区电影| 制服诱惑二区| 亚洲国产精品国产精品| 卡戴珊不雅视频在线播放| 伊人久久精品亚洲午夜| 国产成人91sexporn| 亚洲欧美中文字幕日韩二区| 亚洲色图 男人天堂 中文字幕 | 日日撸夜夜添| 秋霞伦理黄片| 久热久热在线精品观看| 日韩 亚洲 欧美在线| 日韩免费高清中文字幕av| 日本午夜av视频| 日韩视频在线欧美| 看非洲黑人一级黄片| 在线亚洲精品国产二区图片欧美 | 国模一区二区三区四区视频| 婷婷成人精品国产| 九草在线视频观看| 日韩大片免费观看网站| 国产在线一区二区三区精| 亚洲国产精品一区二区三区在线| 黑人猛操日本美女一级片| 国产探花极品一区二区| 欧美少妇被猛烈插入视频| 飞空精品影院首页| 国产精品偷伦视频观看了| 国产亚洲av片在线观看秒播厂| 啦啦啦中文免费视频观看日本| 一级,二级,三级黄色视频| 另类亚洲欧美激情| 99九九在线精品视频| 国产熟女午夜一区二区三区 | 纵有疾风起免费观看全集完整版| 欧美成人精品欧美一级黄| 国产一区二区在线观看av| 国产精品不卡视频一区二区| 91在线精品国自产拍蜜月| 日韩大片免费观看网站| 亚洲欧洲国产日韩| 少妇被粗大猛烈的视频| 最新中文字幕久久久久| 成年人免费黄色播放视频| 在线观看国产h片| 亚洲中文av在线| 亚洲国产成人一精品久久久| 黄色怎么调成土黄色| 少妇人妻 视频| 免费人妻精品一区二区三区视频| 亚洲精品视频女| 高清午夜精品一区二区三区| 成年美女黄网站色视频大全免费 | 久久精品国产亚洲av涩爱| 欧美 日韩 精品 国产| 最黄视频免费看| 两个人免费观看高清视频| 亚洲欧美日韩另类电影网站| 一区二区三区四区激情视频| 最近中文字幕高清免费大全6| 亚洲国产精品一区三区| 国产亚洲午夜精品一区二区久久| a级毛片免费高清观看在线播放| 免费av中文字幕在线| 中文欧美无线码| 我的老师免费观看完整版| 亚洲国产精品一区二区三区在线| 国产熟女欧美一区二区| 国产高清不卡午夜福利| 欧美激情极品国产一区二区三区 | 纯流量卡能插随身wifi吗| 国产无遮挡羞羞视频在线观看| 少妇熟女欧美另类| 九色亚洲精品在线播放| 99热全是精品| 韩国高清视频一区二区三区| 亚洲欧美日韩另类电影网站| 国产日韩欧美视频二区| 久久人人爽人人片av| 91在线精品国自产拍蜜月| 免费观看a级毛片全部| 国产免费又黄又爽又色| 亚洲欧美清纯卡通| 色哟哟·www| 日韩av在线免费看完整版不卡| 人妻制服诱惑在线中文字幕| 国产深夜福利视频在线观看| 啦啦啦中文免费视频观看日本| 这个男人来自地球电影免费观看 | 能在线免费看毛片的网站| 久久人人爽av亚洲精品天堂| a级毛片免费高清观看在线播放| av福利片在线| 在线观看免费日韩欧美大片 | 精品卡一卡二卡四卡免费| 亚洲av.av天堂| 在线观看免费视频网站a站| 久久久精品区二区三区| 日日撸夜夜添| 国产成人精品在线电影| 狂野欧美激情性bbbbbb| 七月丁香在线播放| 日本黄大片高清| av有码第一页| 精品酒店卫生间| 伦精品一区二区三区| 18禁观看日本| 97精品久久久久久久久久精品| 亚洲国产av新网站| 久久99蜜桃精品久久| 国产熟女欧美一区二区| 三上悠亚av全集在线观看| av电影中文网址| 男男h啪啪无遮挡| 午夜91福利影院| 亚洲国产色片| 三级国产精品片| 91精品国产九色| 久久久a久久爽久久v久久| 久久久久久久国产电影| 美女xxoo啪啪120秒动态图| av在线观看视频网站免费| 国产毛片在线视频| 久久久久久久国产电影| 欧美3d第一页| 秋霞在线观看毛片| 婷婷色av中文字幕| 伊人久久国产一区二区| 中文字幕久久专区| 亚洲无线观看免费| 亚洲av成人精品一区久久| 亚洲精品国产av蜜桃| 91精品国产九色| 久久精品国产亚洲网站| 亚洲精品第二区| 免费大片18禁| 亚洲伊人久久精品综合| 色吧在线观看| 3wmmmm亚洲av在线观看| 亚洲综合色惰| 亚洲综合色网址| 亚洲av电影在线观看一区二区三区| av国产精品久久久久影院| 亚洲四区av| 好男人视频免费观看在线| 亚洲精品国产色婷婷电影| 久久久久视频综合| 夫妻午夜视频| 日本爱情动作片www.在线观看| 91精品伊人久久大香线蕉| 国产老妇伦熟女老妇高清| 欧美丝袜亚洲另类| 乱码一卡2卡4卡精品| 午夜视频国产福利| 欧美激情极品国产一区二区三区 | 久久精品国产鲁丝片午夜精品| 国产免费又黄又爽又色| 色5月婷婷丁香| 免费av不卡在线播放| 大话2 男鬼变身卡| 午夜av观看不卡| 美女脱内裤让男人舔精品视频| 免费大片黄手机在线观看| 91精品三级在线观看| 97在线人人人人妻| 国产成人aa在线观看| 亚洲婷婷狠狠爱综合网| 欧美日韩视频高清一区二区三区二| 大片电影免费在线观看免费| 亚洲av福利一区| av.在线天堂| 亚洲精华国产精华液的使用体验| 国产高清不卡午夜福利| 国产av码专区亚洲av| 国产av精品麻豆| 午夜福利在线观看免费完整高清在| 夜夜看夜夜爽夜夜摸| 天天躁夜夜躁狠狠久久av| 国产成人精品在线电影| 你懂的网址亚洲精品在线观看| 80岁老熟妇乱子伦牲交| 久久久国产欧美日韩av| 国产精品一二三区在线看| 最近2019中文字幕mv第一页| av国产精品久久久久影院| 久久久精品94久久精品| 美女大奶头黄色视频| 99re6热这里在线精品视频| 久久久久久久亚洲中文字幕| 在线观看www视频免费| 99热全是精品| 999精品在线视频| 久久精品国产a三级三级三级| 成年av动漫网址| 自线自在国产av| 久久精品国产a三级三级三级| 亚洲国产色片| 亚洲成人一二三区av| 性色av一级| 99九九在线精品视频| 国产在视频线精品| 简卡轻食公司| 特大巨黑吊av在线直播| 日韩一区二区视频免费看| 国产精品.久久久| 亚洲精品日韩在线中文字幕| 桃花免费在线播放| 人人澡人人妻人| 大片电影免费在线观看免费| 精品人妻偷拍中文字幕| 亚洲欧美精品自产自拍| 一级毛片黄色毛片免费观看视频| 高清午夜精品一区二区三区| 边亲边吃奶的免费视频| 成人黄色视频免费在线看| 我的女老师完整版在线观看| 国产成人av激情在线播放 | 99九九在线精品视频| 狂野欧美激情性bbbbbb| 老女人水多毛片| av电影中文网址| 青青草视频在线视频观看| 亚洲精品日本国产第一区| 三上悠亚av全集在线观看| 久久精品人人爽人人爽视色| 亚洲av欧美aⅴ国产| 欧美精品人与动牲交sv欧美| 日韩一本色道免费dvd| 国产精品一国产av| 久久久久久久久大av| 男女边吃奶边做爰视频| 免费观看无遮挡的男女| 亚洲欧美日韩另类电影网站| 国产乱人偷精品视频| 在线精品无人区一区二区三| 亚洲av欧美aⅴ国产| 人人妻人人澡人人看| 中国三级夫妇交换| 天堂8中文在线网| 99久国产av精品国产电影| 国产一级毛片在线| 国产午夜精品一二区理论片| 日本免费在线观看一区| 久久精品久久久久久噜噜老黄| 搡老乐熟女国产| 下体分泌物呈黄色| 一个人免费看片子| 亚洲人成77777在线视频| 国产熟女午夜一区二区三区 | 一级黄片播放器| 久久国产精品男人的天堂亚洲 | av有码第一页| 国产 一区精品| 男女高潮啪啪啪动态图| 国产亚洲精品久久久com| 青春草国产在线视频| 女的被弄到高潮叫床怎么办| 黄色欧美视频在线观看| 欧美日本中文国产一区发布| 在线看a的网站| 精品人妻偷拍中文字幕| 一区二区av电影网| 五月玫瑰六月丁香| 国产伦精品一区二区三区视频9| 天天躁夜夜躁狠狠久久av| 成年人免费黄色播放视频| 国产精品不卡视频一区二区| 中文字幕人妻熟人妻熟丝袜美| 欧美+日韩+精品| 天堂中文最新版在线下载| 亚洲欧洲国产日韩| av一本久久久久| 女性被躁到高潮视频| 极品少妇高潮喷水抽搐| 国产亚洲精品第一综合不卡 | 国产成人av激情在线播放 | 纯流量卡能插随身wifi吗| 亚洲国产成人一精品久久久| 三上悠亚av全集在线观看| 高清视频免费观看一区二区| 亚洲人成77777在线视频| 人体艺术视频欧美日本| 午夜日本视频在线| 最黄视频免费看| 观看av在线不卡| 国产精品 国内视频| 肉色欧美久久久久久久蜜桃| 欧美+日韩+精品| 热99久久久久精品小说推荐| 夫妻性生交免费视频一级片| 欧美精品人与动牲交sv欧美| 成人毛片a级毛片在线播放| 欧美亚洲 丝袜 人妻 在线| 亚洲欧洲国产日韩| 欧美性感艳星| 老司机影院成人| 国产精品99久久99久久久不卡 | 男女无遮挡免费网站观看| 十八禁高潮呻吟视频| 中文字幕亚洲精品专区| 亚洲欧美一区二区三区国产| 男人爽女人下面视频在线观看| 晚上一个人看的免费电影| 成人午夜精彩视频在线观看| av.在线天堂| 日韩制服骚丝袜av| 黄色怎么调成土黄色| 国产精品久久久久久精品电影小说| 亚洲图色成人| 亚洲av综合色区一区| 少妇猛男粗大的猛烈进出视频| 亚洲婷婷狠狠爱综合网| 在线观看国产h片| 国产精品女同一区二区软件| 久久久久久久精品精品| 色婷婷久久久亚洲欧美| 乱人伦中国视频| 日韩欧美一区视频在线观看| 亚洲av二区三区四区| 成人黄色视频免费在线看| 美女主播在线视频| 最近最新中文字幕免费大全7| 只有这里有精品99| 免费人成在线观看视频色| 国产精品国产三级国产专区5o| 亚洲人与动物交配视频| 一个人免费看片子| 亚洲人成网站在线观看播放| 国产成人精品久久久久久| 久久精品国产自在天天线| 一级片'在线观看视频| 亚洲综合精品二区| 在线观看人妻少妇| 亚洲欧美一区二区三区黑人 | 日韩大片免费观看网站| 国产综合精华液| 久久久久久久大尺度免费视频| 大香蕉久久网| 视频区图区小说| 国模一区二区三区四区视频| 女的被弄到高潮叫床怎么办| 成人免费观看视频高清| 亚洲怡红院男人天堂| 91精品伊人久久大香线蕉| 亚洲精品一区蜜桃| 日韩熟女老妇一区二区性免费视频| 少妇被粗大猛烈的视频| 色视频在线一区二区三区| 在线亚洲精品国产二区图片欧美 | 欧美另类一区| 性色av一级| 免费大片黄手机在线观看| 亚洲天堂av无毛| 午夜免费鲁丝| 久久久久久久大尺度免费视频| 亚洲国产精品999| 中文天堂在线官网| 日韩精品有码人妻一区| 搡老乐熟女国产| av女优亚洲男人天堂| 亚洲国产欧美日韩在线播放| 99国产精品免费福利视频| 我的老师免费观看完整版| 日日撸夜夜添| 午夜福利视频精品| 亚洲国产精品一区三区| 精品卡一卡二卡四卡免费| 如何舔出高潮| 精品少妇黑人巨大在线播放| 蜜臀久久99精品久久宅男| 免费av中文字幕在线| 欧美精品亚洲一区二区| 欧美人与善性xxx| 国产国语露脸激情在线看| 精品久久久精品久久久| 男女无遮挡免费网站观看| 日韩av在线免费看完整版不卡| 最近中文字幕高清免费大全6| 女的被弄到高潮叫床怎么办| 久久国产精品男人的天堂亚洲 | 欧美国产精品一级二级三级| 欧美性感艳星| 超碰97精品在线观看| 乱人伦中国视频| 亚洲精品日韩在线中文字幕| 王馨瑶露胸无遮挡在线观看| 亚洲少妇的诱惑av| 天天躁夜夜躁狠狠久久av| 亚洲情色 制服丝袜| 国产精品久久久久久精品古装| 性高湖久久久久久久久免费观看| 五月玫瑰六月丁香| 天堂俺去俺来也www色官网| .国产精品久久| 男女高潮啪啪啪动态图| 美女视频免费永久观看网站| 夜夜看夜夜爽夜夜摸| 狂野欧美激情性xxxx在线观看| 丝袜喷水一区| 国产黄频视频在线观看| 91精品伊人久久大香线蕉| 亚洲国产毛片av蜜桃av| 天天影视国产精品| 欧美精品一区二区大全| 国产成人91sexporn| 亚洲欧美日韩卡通动漫| 免费黄网站久久成人精品| 欧美日韩综合久久久久久| 色哟哟·www| 亚洲,一卡二卡三卡| 韩国高清视频一区二区三区| 国产色婷婷99| 永久网站在线| 亚洲美女视频黄频| videossex国产| 精品久久久噜噜| 99视频精品全部免费 在线| 观看av在线不卡| av国产久精品久网站免费入址| 熟妇人妻不卡中文字幕| 天堂中文最新版在线下载| 国产片内射在线| 韩国av在线不卡| 看免费成人av毛片| 亚洲美女黄色视频免费看| 少妇猛男粗大的猛烈进出视频| 新久久久久国产一级毛片| 两个人的视频大全免费| 视频中文字幕在线观看| 久久精品国产a三级三级三级| 老司机影院毛片| 欧美性感艳星| 免费黄频网站在线观看国产| 国产黄色免费在线视频| 久久午夜综合久久蜜桃| 亚洲av不卡在线观看| 国产男女超爽视频在线观看| 国产男人的电影天堂91| 亚洲精品国产av蜜桃| 看十八女毛片水多多多| 自拍欧美九色日韩亚洲蝌蚪91| 韩国高清视频一区二区三区| 男女啪啪激烈高潮av片| 18禁裸乳无遮挡动漫免费视频| 亚洲高清免费不卡视频| 99久久精品一区二区三区| 国产精品国产三级国产av玫瑰| 一级黄片播放器| 精品人妻一区二区三区麻豆| 亚洲少妇的诱惑av| 亚洲激情五月婷婷啪啪| 精品国产国语对白av| 91精品一卡2卡3卡4卡| 精品午夜福利在线看| 国产一级毛片在线| 人妻夜夜爽99麻豆av| 国产高清国产精品国产三级| 日日啪夜夜爽| 中文欧美无线码| 色婷婷久久久亚洲欧美| av有码第一页| 国产亚洲欧美精品永久| 超碰97精品在线观看| 人人妻人人澡人人爽人人夜夜| 亚洲av男天堂| 亚洲av不卡在线观看| 大香蕉久久成人网| av女优亚洲男人天堂| 国产白丝娇喘喷水9色精品| 制服人妻中文乱码| 91成人精品电影| 日本午夜av视频| 水蜜桃什么品种好| 最近中文字幕2019免费版| 国产亚洲最大av| 免费黄频网站在线观看国产| 成人国产麻豆网| 欧美亚洲日本最大视频资源| 久久免费观看电影| 久久国产亚洲av麻豆专区| 日韩在线高清观看一区二区三区| 久久久久久久久大av| 国产一级毛片在线| 精品人妻偷拍中文字幕| 插阴视频在线观看视频| 看免费成人av毛片| 国产精品一区二区在线不卡| 日韩精品免费视频一区二区三区 | h视频一区二区三区| 草草在线视频免费看| 嘟嘟电影网在线观看| 18在线观看网站| 女性被躁到高潮视频| 99视频精品全部免费 在线| av在线播放精品| 赤兔流量卡办理| 中文精品一卡2卡3卡4更新| 亚洲av在线观看美女高潮| 制服诱惑二区| 秋霞在线观看毛片| 国产成人精品无人区| 亚洲婷婷狠狠爱综合网| 国产一区二区三区综合在线观看 | 精品久久久久久久久av| 久久久久久人妻| 91aial.com中文字幕在线观看| 欧美日韩av久久| 亚洲精品aⅴ在线观看| 精品一区二区三区视频在线| 高清午夜精品一区二区三区| 国产成人a∨麻豆精品| 99国产精品免费福利视频| 丰满乱子伦码专区| 精品久久久噜噜| 女人精品久久久久毛片| 三上悠亚av全集在线观看| 成人手机av| 亚洲精品日本国产第一区| 少妇高潮的动态图| 青青草视频在线视频观看| 男人添女人高潮全过程视频| 色婷婷av一区二区三区视频| 亚洲av综合色区一区| 久久免费观看电影| 精品久久久精品久久久| 99国产精品免费福利视频| 伦理电影免费视频| 乱码一卡2卡4卡精品| 亚洲熟女精品中文字幕| av女优亚洲男人天堂| .国产精品久久| 久久99热这里只频精品6学生| 欧美激情 高清一区二区三区| 大片电影免费在线观看免费| 久久久久久久亚洲中文字幕| 又大又黄又爽视频免费| 亚洲精品国产av蜜桃| av免费观看日本| 狠狠精品人妻久久久久久综合| 国产熟女午夜一区二区三区 | 色5月婷婷丁香| 久久午夜综合久久蜜桃| 18+在线观看网站| 成人影院久久| 久久毛片免费看一区二区三区| 国国产精品蜜臀av免费| 国产日韩欧美在线精品| 97超视频在线观看视频| 少妇精品久久久久久久| 菩萨蛮人人尽说江南好唐韦庄| 内地一区二区视频在线| 免费播放大片免费观看视频在线观看| 啦啦啦啦在线视频资源| 亚洲av中文av极速乱| 夜夜骑夜夜射夜夜干| 五月开心婷婷网| 天美传媒精品一区二区| 国产成人午夜福利电影在线观看| 有码 亚洲区| 欧美精品高潮呻吟av久久| 极品人妻少妇av视频| 如何舔出高潮| 国产午夜精品久久久久久一区二区三区| 在线 av 中文字幕| 寂寞人妻少妇视频99o| 91精品国产国语对白视频| 国产精品人妻久久久久久| 国产成人精品婷婷| 18禁动态无遮挡网站| 亚洲天堂av无毛| 亚洲成人一二三区av| av在线观看视频网站免费| 狂野欧美白嫩少妇大欣赏| 两个人的视频大全免费| 黄色配什么色好看| 精品午夜福利在线看| 一级爰片在线观看| 五月玫瑰六月丁香| 在线亚洲精品国产二区图片欧美 | 天堂8中文在线网| 80岁老熟妇乱子伦牲交| 久久人人爽人人爽人人片va| 美女中出高潮动态图| 日韩av免费高清视频| 精品久久久久久久久亚洲| 在线 av 中文字幕| 亚洲无线观看免费| 有码 亚洲区| 久久婷婷青草| 亚洲欧洲日产国产| 亚洲性久久影院| 日本av手机在线免费观看| 亚洲性久久影院| 久久国产亚洲av麻豆专区| 久久久久久久久久久免费av| 精品国产国语对白av| 夜夜看夜夜爽夜夜摸| 亚洲欧美中文字幕日韩二区| 日韩av在线免费看完整版不卡| 青春草国产在线视频| 亚洲色图 男人天堂 中文字幕 | 国产精品一区二区在线不卡| 色5月婷婷丁香| 亚洲精品久久午夜乱码| 久久精品人人爽人人爽视色| 99九九在线精品视频|