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

    Improvement potential analysis for integrated fractionating and heat exchange processes in delayed coking units☆

    2016-06-01 03:00:58YangLeiDanlinZengGuanghuiWang

    Yang Lei,Danlin Zeng,Guanghui Wang

    School of Chemical Engineering and Technology,Wuhan University of Science and Technology,Wuhan 430081,China

    1.Introduction

    Fractionating and heatexchange processes are commonly integrated system in re fineries,such as complex fractionator and heat exchanger network(HEN)systems in delayed coking units(DCUs)[1].Analysis and optimization of these energy integrated systems are mainly performed with mathematical modeling and thermodynamic analysis.

    Currently,exergy analysis is conducted to identify and quantify the parts of systems with excessive irreversibility and determine the thermodynamic efficiency of fractionating or heat exchange processes[2-4].It describes the maximum available work when a type of energy is converted reversibly to a reference system in equilibrium with environmental conditions[5]and is a powerful tool for reducing the energy demand by locating inefficiencies from exergy losses of the system[6].Recently,pinch analysis and exergy analysis are combined to identify related thermodynamic imperfection with environmental concerns based on graph methods.Ishida and Ohno[7]presented a diagram(energy-direction factor diagram)to analyze exergy destructions based on the conception of a reversible distillation column,which is only acceptable to binary mixtures where the composition is uniquely determined by temperature.A new method for graphic exergy analysis of multicomponent distillation and for visualization of rate characteristics of a distillation column was proposed by using an energy utilization diagram(EUD)[8,9].Taprap and Ishida[10],by using a mass transfer model coupled with a mass and heat balance,developed the EUD for a distillation column in relation to thex-ydiagram in a logarithmic scale.The diagrams identified the amount of energy transformation and exergy loss of individual process,the heating and mixing of the liquid phase,cooling and mixing of the vapor phase and phase changes.By using a graphicalexergy analysis,Jin and Ishida[11]introduced EUDs to analyze advanced thermal power systems.Suggestions for improving the systems were described based on EUDs.Jinet al.[12]also used the graphical exergy methodology based on EUDs to analyze gas turbine cycle,based on the combination of humid air turbine cycle with an externally fired heater.Bouapetchet al.[13]extended EUDs to power availability for exergy loss analysis.A Modified subsystem power availability diagram was proposed to present the graphical exergy distribution.Using material utilization diagram,Budimanet al.[14]analyzed separation performance and exergy loss phenomena for a retro fitted distillation column.

    The exergy loss pro file,characterized by its straightforwardness,has been widely used for a distillation column.For example,Le Goffet al.[15]introduced exergy loss pro files to design and optimize a distillation column.Budiman and Ishida[16]presented the exergy utilization diagram to indentify the targets for reducing energy and exergy consumptions.Khoaet al.[17]reported a novel method to analyze the exergy loss of the effect of design and operating parameters,which is graphically visualized by three dimensional exergy analysis curves,which revealthe correlations between exergy loss,design and operating parameters in a distillation column.

    Feng and Zhu[18]used the strength of pinch analysis,in which system information is represented using simple diagrams(e.g.composite curves,grand composite curves),and obtained targets for the system under consideration prior to design.By combining the strengths of pinch and exergy analysis methods,the method could represent a whole system,including individualunitson one diagram,which ishelpful to quickly screen the promising modifications for improving a base case design.In order to address the energy integration,Anantharamanet al.[19]had Modified and extended the concept of energy level proposed by Feng and Zhu.As pressure,temperature and composition changes were taken into account when developing the theory for this method,energy level composite curves(ELCCs)could be applied to a wide range of processes and in particular to energy intensive chemical plants.

    Exergy analysis identifies the system components with the highest thermodynamic inefficiencies and related processes,but only a part of the exergy destructions in the system could be avoided in practice.Based on the concepts of avoidable/unavoidable exergy destructions and investment costs,Weiet al.[20]presented an exergy analysis and an exergoeconomic evaluation to identify the energy-saving potentials in distillation processes.

    Previous researches mainly concentrated on the exergy analysis for individual processes[7-10,14-17,20],such as fractionating or heat exchange,rarely for integrated systems[1].In this study,extending the concept of traditional ELCC,an advanced ELCC considering temperature dependence of heat capacity and integrating fractionating and heat exchange processes is constructed by using simulation results and a stepwise procedure.The concepts of avoidable/unavoidable exergy destructions are utilized to identify practically and economically possible improvements in the integrated system,which would give the maximum potential for process improvements.The application benefits of the method are illustrated with a DCU.

    2.Advanced ELCC and Exergy Analysis

    2.1.Advanced ELCC

    In order to evaluate the quality of the energy fora system,the energy level for heat is de fined as[18]

    whereT0is the temperature in ambient conditions,andTis the temperature for streams.

    For a steady-state- flow system,the heat exergy can be expressed as

    The ε-Qdiagram is the basis of the ELCC,with the energy level(ε)and enthalpy(Qor ΔH)as the ordinate and abscissa,respectively,as shown in Fig.1.For a complex fractionator with superheated vapor feed and pump-around heat removal,the ε-Qdiagram for heat exchange and fractionating processes is shown in Fig.1(b).

    When the thermal energy is concerned,the energy levelis expressed as Eq.(1),and ΔHis equivalent to the exchanged heatQ.The heat exergy is deduced to

    The ε-Qdiagram is an upgraded version of temperature-enthalpy(T-H)diagram.On the premise of no heat dissipating,the energy composite curves(so-calledT-Hdiagram)for a simple heat exchange process can be converted into the energy level composite curve(so-calledε-Qdiagram).The shaded area represents exergy destruction in the heat transfer or fractionating process as shown in Fig.1.

    The complex fractionating and heat exchange processes can be revealed by ELCC in Fig.2(a)and(b),respectively.The two processes can be integrated in one figure as shown in Fig.2(c).To decrease the exergy destruction for the fractionating process,some modifications can be implemented for improving the energy-use efficiency.For example,as shown in Fig.2(d),as the return temperature of the pumparound stream increases,the curve for heat removal is changed.The new exergy distribution for the fractionating process is illustrated by a dotted line.The curve for pump-around heat removal,which is the hot stream in the heat exchange process,approaches to the curve for feed.Therefore,some parameters are varied for the heat exchange process,including the outlet temperature of the hot and cold streams in the heat exchanger.The new exergy distribution for the heat exchange process is illustrated in Fig.2(e)by a dotted line.These two improved processes can be integrated in one figure as shown in Fig.2(f).It reveals that a part of the decreased exergy destruction in the fractionating process is diverted into the exergy destruction in the heat exchange process.More exergy can be recovered.

    Hence,in the ELCC for the integrated system,the exergy destruction distribution is clearly indicated in one figure.More accurate practical results can be obtained by considering temperature dependence of heat capacity.The interactions between different processes are identified when the processes are improved.

    Energy is generally transferred from fractionating processes at higher energy levels to heat exchange processes at lower energy levels.The interaction between fractionating processes(energy donating)and heatexchange processes(energy accepting)can be properly revealed in the advanced ELCC.Therefore,both energy and exergy balances for a whole system can be represented simultaneously.

    Astep wise procedure forconstruction of advanced ELCC isdescribed as follows.

    Step 1:To simplify the flow sheet and de fine the boundary of the chosen integrated system,its envelope should be de fined first containing the fractionator and all heat exchangers involving air and water coolers.Then the simulation for the system is set up in order to obtain the physical properties of involved streams.

    Step 2:The curves in the ELCC are subdivided into two categories:one represents materials(e.g.feed or product),another represents energy flow(e.g.heat removal).The final temperature of the material is de fined as 288.15 K,which is different with energy flows.The final temperature of heat removal is the return temperature.Especially,the cold utilities cannot be plotted in the ELCC such as cooling water and air,which move discharged exergy from selected system to the surroundings.

    Step 3:Forevery processin the integrated system,the temperature and corresponding enthalpy formaterialor energy flows are obtained based on the simulation.The temperature is converted to the energy levelby calculating Eq.(1).In the ELCC,the ordinate and abscissa represent the energy level and enthalpy difference,respectively.To obtain more accurate fitting results,as many data as possible are de fined to plot a curve for presenting material or energy flows.

    Step 4:The curves for the energy donating/accepting terms of the fractionator can be moved horizontally,relative to one another such that the higher energy level curve is to the left of the lowerenergy levelcurve and cold streamsin the heatexchange process are moved horizontally until they are covered completely by corresponding hot streams.

    Fig.1.ε-Q diagram.(a)Heat exchange process;(b)complex fractionating process.

    Step 5:The least squares method is used to fitting the equation of each curve.These equations are similar to

    whereaandbare the corresponding polynomial fitting coefficients.Equations with a higher degree can ensure the R-square more than 0.9999.In this equation,temperature dependence of heat capacity is considered.

    Step 6:The exergy destruction distribution of the integrated system is determined by calculating the integral areas in the advanced ELCC.

    Similar to the composite curves of pinch analysis,there is no limit to move the curves horizontally.This graphical representation of energy donating and accepting composite curves gives physical insight to an engineer regarding energy degradation and interaction between fractionating and heat exchange processes.The principle of energy transfer from a higher value of ε to a lower value is valid as higher temperature streams will have a higher ε from Eq.(1).Hence,this graphical representation,developed to the analogous pinch representation for heat recovery,can effectively reveal the performance of individual process units indicating the exergy donating/accepting trends.It is worth noting that the heat capacity of each stream is not considered as a constant in the advanced ELCC.Therefore,the calculated results of exergy destructions are more accurate to suggest the potential for improvements.

    2.2.Exergy analysis based on avoidable/unavoidable exergy destructions

    The exergy destruction rate˙ED,kassociated with thekth component of the system consists of one avoidable(superscript AV)and another unavoidable(superscript UN)part.An exergy balance and the total exergy destruction in a process can be described as[18]

    where˙EF,k,˙EP,kand˙EL,kdenote the fuel,product,and exergy loss of thekth component,respectively.

    Fig.2.ELCC for an integrated system with process improvement.(a)Fractionating process;(b)heatexchange process;(c)integrated system;(d)fractionating process with improvement;(e)heat exchange process with improvement;(f)integrated system with process improvement.

    The thermodynamic inefficiency of a component consists exclusively of exergy loss(i.e.=0)in this study.Hence,exergy efficiency,ηk,related to avoidable exergy destruction,indicates the comparison of exergy efficiency of different components,de fined as

    and the potential for improvements for an individual component is measured by the ratio

    Practical and valuable improvements can be found by the two evaluation indexes through calculating Eqs.(7)and(8)for the integrated system.The accurate estimation of unavoidable exergy destruction for the heat exchange and fractionating processes is crucial to the quantitative evaluation and identification for potential improvements.

    The exergy destructionis calculated by the advanced ELCC through the difference of integral area between hot and cold stream curves for an adiabatic heat exchange process.

    Based on the concept of avoidable/unavoidable exergy destructions,for an adiabatic heat exchange process is given by[18]

    whereis the heat duty of a heat exchanger,Thais the mean temperature of the hot stream,Tcais the mean temperature of the cold stream,Ris the universal gas constant,pis the pressure,and Δpis the pressure drop.For simplicity,the pressure drop for hot and cold streams is ignored in the heat exchanger(Δph= Δpc=0).

    The optimal temperature approach ΔToptis applied to calculate(Tha-Tca)opt[21]

    Under the current economic and technical conditions(tradeoff between the costs of fuel and heat exchanger investment),some optimal temperature differences are presented for different conditions[21].In detail,the optimal temperature difference is de fined as 15 K for the heat exchangers,in which there is no phase transition for streams.Otherwise,the temperature approach is de fined as 10 K.

    Thus,the unavoidable exergy destruction rateof a heat exchanger is calculated from Eqs.(9)and(10).Furthermore,we can analyze and evaluate the energy-saving potential and efficiency of energy use for different heat exchangers by Eqs.(7)and(8).

    For an adiabatic fractionating process,the total exergy destruction is

    whereanddenote the exergy destruction of an irreversible separation process and the thermodynamic exergy difference,respectively.

    The minimum or unavoidable exergy destruction is determined by the reversible fractionating process.Under ideal conditions,when fractionating is a reversible process for mixing,the minimumexergy destruction is the thermodynamic exergy difference for fractionating,which can be calculated by[22]whereNjdenotes the mole ratio between fractionjand allfractions,andR=8.314 J·mol-1·K-1.

    is calculated using Eqs.(11)and(12),in which energy-saving potential of process improvements is de fined for the fractionator.On the basis of the above analysis,andcorrespond to avoidable and unavoidable exergy destructions for the fractionating process,respectively.

    Thus far,the avoidable exergy destruction is usually applied to analyze the process inefficiency for single processes,while the interaction between fractionating and heat exchange processes is rarely considered.In general,energy efficiency improvement will be effectively arrived by considering interactions between different processes.

    3.Improving Potential Analysis for a Delayed Coking Unit

    A DCU in a Chinese re finery is considered as a case study to reveal the exergy analysis for an integrated system based on the advanced ELCC,which is shown in Fig.3.The simulation strategy about the complex fractionator can be found in Ref.[23].A steady state simulation is set up using the simulation system PRO/II.The data for the improving potential analysis are from the simulation results,which are valid compared with the operational data.

    3.1.Integrated system and its ELCC

    AconventionalDCUgenerally consists of three sections,i.e.coking,oil and gas separation,and heat exchange network.Envelope I is chosen as the boundary of the flow sheet for exergy analysis in Fig.3,which represents the fractionating and heat exchange processes.It includes a complex fractionator,10 heat exchangers,3 water coolers and 4 air coolers.The heat sources include pump-around heat removal and products donated by the complex fractionator.The heatsinks contain feeds,reboiling streams in the absorption and stabilization system and feed water for steam generation.The integrated system looks simple,but the interaction between the components is complex.Forexample,the heatremoval distribution or the operating parameters of the fractionator may be varied in the fractionating process optimization.The fractionator is optimized in isolation,but the HEN may be influenced by those changes so that the utility cost increases obviously.For clarity,the system flow sheet is simplified to Fig.4,where solid lines represent material flows,dotted lines represent energy flows,and black boxes refer to several sections of the column,in which the mass transfer is not considered.

    It is assumed that the system is adiabatic and the pressure exergy is ignored.Then the constructing method is also used to generate the ELCC for the case study as shown in Fig.5.

    Fitting equations are de fined for hot and cold stream curves by Eq.(4)and the results are listed in Table 1.For this integrated system,the simplified exergy destruction calculation results without considering temperature dependence of heatcapacity are compared with the results from the advanced ELCC in Fig.6.The results of some components are different,so it is necessary to consider non-linear heat capacities varied with temperature for some material and energy flows.

    Moreover,basic information for calculating the unavoidable exergy destruction of the fractionator is listed in Table 2.According to Eq.(12),the unavoidable exergy destruction of the fractionator in this case study is given by

    3.2.Exergy analysis and evaluation

    Fig.3.Flow diagram for a DCU.

    Fig.4.Schematic diagram of an integrated system.

    Fig.5.Advanced ELCC for the case study.

    Table 3 reveals the unavoidable exergy destruction rate,avoidable exergy destruction rate,the ratio ofand,and the exergy efficiency ηkwhich are calculated for each component.Corresponding to Fig.4,avoidable and unavoidable destruction distributions of the integrated fractionating and heat exchange processes are shown in Fig.7.The rateand ratioare the thermodynamic measures for the potential of improving the performance of thekth component in isolation,which are shown in Fig.8.

    3.3.Discussion

    The main contributors to the exergy destruction in the integrated system are in descending order of importance:HE05,HE08,fractionator,HE10,HE06,HE04,HE09,HE07,HE03,HE01,and HE02.However,indicates the potential of improvements when the component is considered in isolation,and represents the percentage of the incurred exergy destruction in a component that could be avoided.Moreover,this variable isde fined as the energy-saving potential.anddemonstrate a dominant part ofwhich can be reduced by using the best technology in some components of the integrated system.In Fig.8,the value ofis more than 90%in HE05 and HE08.A significant part of the exergy destruction in these high temperature heat exchangers is avoidable.The avoidable exergy destruction is mainly caused by ahigher temperature approach between hot and cold streams.For example,the cold stream ofHE08 is the deaerated water,which is heated to generate steam at 1.0 MPa by the intermediate pump-around.Because the energy level for the cold stream in HE03 is higher than that in HE08,we suggest increasing the area of HE03 so as to decrease the temperature of the hot stream in HE08.

    Table 1Basic information and calculated results of heat exchange processes

    Fig.6.Comparison of the results for the simplified calculation method and the advanced ELCC.

    Table 2Basic information for calculating the unavoidable exergy destruction of the fractionator

    Table 3Results from the exergy analysis including avoidable and unavoidable exergy destruction rates and evaluation for selected components

    Meanwhile,exergy efficiencies without considering avoidable/unavoidable exergy destructions for the components are ranked in descending order as follows:HE02,fractionator,HE03,HE04,HE01,HE07,HE08,HE09,HE06,HE10,and HE05.The results of conventional exergy analysis can be strongly supplemented by the exergy analysis when considering avoidable/unavoidable exergy destructions.According to the irorigins,irreversibilities identified in the conventionalexergy analysis can be split.Only part of the irreversibilities thatcan be avoided should be considered for the improvement of a system with the integrated processes.The values of ηkfor the components are ranked in descending order as follows:HE02,fractionator,HE01,HE03,HE04,HE07,HE09,HE08,HE10,HE06,and HE05.In a conventional exergy analysis,the lower exergy efficiency is generally considered to have a greater retro fit potential.It unreasonably indicates that the fractionator has a great energy-saving potential.However,due to the fractionating requirements,the valuesand ηkfor the fractionator are 38.1%and 97.3%,respectively.Compared with other components,it is obvious that the fractionator is not the most promising component for improvement.

    Fig.7.Avoidable and unavoidable destruction distributions of the integrated fractionating and heat exchange processes.

    Fig.8.Exergy evaluation for the case study.

    Based on the temperature pro file of the integrated system as shown in Fig.5,some improvements can be effectively evaluated.The fractionation and heat exchange processes in an integrated system tightly interact with each other.When the exergy destruction of fractionator is decreased,the total exergy destruction of the integrated system may be decreased less than that of the fractionator.In order to improve the energy efficiency of the fractionation process,the exergy destruction should be decreased.According to Fig.5,the curves representing energy accepting terms of the fractionator should be closer to the ones representing energy donating terms of the fractionator.For the heat exchange process,the curves representing the cold streams in HENmay be close to the ones representing energy accepting terms of the fractionator.For example,when the return temperature of the pump-around is increased and heatremovalduty is unchanged,the exergy destruction of fractionator is decreased.However,by the effect of heat exchange process,the total exergy destruction of the integrated system may not have the same reduction with the fractionator.Similar to the pinch method,the integration opportunities between the system and the streams outside the system could be appointed in the ELCC,in which the energy level is considered.

    4.Conclusions

    By applying the advanced ELCC and concepts of avoidable/unavoidable exergy destructions,a novel method of exergy analysis is introduced for improving integrated fractionating and heat exchange processes.In the advanced ELCC,temperature dependence of heat capacity is considered and fractionating and heat exchange processes are integrated in one figure.More accurate results of exergy destructions are obtained.

    An integrated system for fractionating and heat exchange processes in a delayed coking unit is analyzed and evaluated to demonstrate the effectiveness of the proposed method.Through the calculation for avoidable/unavoidable exergy destructions,the major cause of inef ficiency in the system is identified.The exergy destruction results with the advanced ELCC are different from that with simplified exergy destruction calculation.More accurate results can be obtained.Furthermore,exergy efficiencies without considering avoidable/unavoidable exergy destructions for the components are different.The results reveal that the heat exchange between gas oil and deethanization gasoline has the highestpotentialfor improvement,and the lowestexergy efficiency of the component in the integrated system is only 29.4%.The improvement potential and exergy efficiency for the fractionator are 38.1%and 97.3%,respectively.Compared with other components,it is obvious that the fractionator is not the most promising component for improvement.

    Nomenclature

    aincidence coefficient

    bincidence coefficient

    ˙Eexergy rate,MW

    ppressure,MPa

    Qheat duty,MW

    Runiversal gas constant,J·mol-1·K-1

    Ttemperature,K

    Tcamean temperature of cold stream,K

    Thamean temperature of hot stream,K

    ε energy level

    η efficiency,%

    Superscripts

    AV avoidable

    UN unavoidable

    Subscripts

    a average

    c cold stream

    D destruction

    F feed

    Frac fractionator

    h hot stream

    in/out inlet stream/outlet stream

    jpetroleum fractions

    kcomponents in the system

    opt optimization

    P product

    0 ambient conditions

    [1]Y.Lei,B.J.Zhang,Q.L.Chen,Thermodynamic analysis and modeling for typical feedpreheating and fractionating processes in delayed cokers,Comput.Aided Chem.Eng.31(2012)655-659.

    [2]Y.Demirel,Thermodynamic analysis of separation systems,Sep.Sci.Technol.39(2004)3897-3942.

    [3]Y.Demirel,Assessment of thermodynamic performances for distillation columns,Int.J.Exergy3(2006)345-361.

    [4]J.M.Dealy,M.E.Weber,Thermodynamic analysis of process efficiency,Appl.Energy6(1980)177-201.

    [5]R.Rivero,Application of the exergy concept in the petroleum re fining and petrochemical industry,Energy Convers.Manag.43(2002)1199-1220.

    [6]Y.Lara,P.Lisbona,A.Martínez,L.M.Romeo,Design and analysis of heat exchanger networks for integrated Ca-looping systems,Appl.Energy111(2013)690-700.

    [7]M.Ishida,T.Ohno,Application of energy direction factor diagram for exergy analysis of a distillation column,J.Chem.Eng.Jpn.16(4)(1983)281-287.

    [8]M.Ishida,R.Taprap,Application of energy-utilization diagram for graphic exergy analysis of multicomponent distillation column,J.Chem.Eng.Jpn.25(4)(1992)396-402.

    [9]M.Ishida,R.Taprap,Introduction of individual energy level for exergy analysis of process systems with multiple components,J.Chem.Eng.Jpn.26(4)(1993)437-439.

    [10]R.Taprap,M.Ishida,Graphic exergy analysis of processes in distillation column by energy-utilization diagrams,AIChE J.42(6)(1996)1633-1641.

    [11]H.G.Jin,M.Ishida,Graphical exergy analysis of complex cycles,Energy18(6)(1993)615-625.

    [12]H.G.Jin,H.B.Zhao,Z.L.Liu,R.X.Cai,A novel EFHAT system and exergy analysis with energy utilization diagram,Energy29(12)(2004)1983-1991.

    [13]W.Bouapetch,P.Tungwungwiwat,T.Srinophakun,Extension of energy-utilisation diagram to power availability for exergy loss analysis,Int.J.Exergy12(3)(2013)323-343.

    [14]A.Budiman,Sutijan,D.R.Sawitri,Graphical exergy analysis of retro fitted distillation column,Int.J.Exergy8(4)(2011)477-493.

    [15]P.Le Goff,T.Cachot,R.Rivero,Exergy analysis of distillation processes,Chem.Eng.Technol.19(6)(1996)478-485.

    [16]A.Budiman,M.Ishida,Three-dimensional graphical exergy analysis of a distillation column,J.Chem.Eng.Jpn29(4)(1996)662-668.

    [17]T.D.Khoa,M.Shuhaimi,H.Hashim,M.H.Panjeshahi,Optimal design of distillation column using three dimensional exergy analysis curves,Energy35(12)(2010)5309-5319.

    [18]X.Feng,X.X.Zhu,Combining pinch and exergy analysis for process modifications,Appl.Therm.Eng.17(3)(1997)246-261.

    [19]R.Anantharaman,O.S.Abbas,T.Gundersen,Energy level composite curves—A new graphical methodology for the integration of energy intensive processes,Appl.Therm.Eng.26(13)(2006)1378-1384.

    [20]Z.Q.Wei,B.J.Zhang,S.Y.Wu,Q.L.Chen,G.Tsatsaronis,Energy-use analysis and evaluation of distillation systems through avoidable exergy destruction and investment costs,Energy42(1)(2012)424-433.

    [21]I.C.Kemp,Pinch analysis and process integration,Butterworth-Heinemann,Amsterdam,2007.

    [22]B.Hua,Estimation of exergy difference in the separating process of petroleum and its distillates,Acta Pet.Sin.(Pet.Process.Sect.)2(3)(1986)102-107.

    [23]Y.Lei,B.J.Zhang,X.Q.Hou,Q.L.Chen,A novel strategy for simulating the main fractionator of delayed cokers by separating the de-superheating process,Chin.J.Chem.Eng.21(3)(2013)285-294.

    伊人亚洲综合成人网| 五月天丁香电影| 亚洲五月色婷婷综合| 亚洲欧美一区二区三区久久| 一级黄色大片毛片| 91麻豆精品激情在线观看国产 | 一本—道久久a久久精品蜜桃钙片| 黄色a级毛片大全视频| 在线 av 中文字幕| 久久人人爽人人片av| 国产精品国产三级国产专区5o| 亚洲国产欧美日韩在线播放| 精品人妻在线不人妻| 欧美xxⅹ黑人| 亚洲伊人久久精品综合| 亚洲精品av麻豆狂野| 国产精品免费视频内射| 国产不卡av网站在线观看| 在线看a的网站| bbb黄色大片| xxxhd国产人妻xxx| 伊人久久大香线蕉亚洲五| 男女床上黄色一级片免费看| 桃红色精品国产亚洲av| 色视频在线一区二区三区| 欧美成人午夜精品| 亚洲成av片中文字幕在线观看| 男女免费视频国产| 久久久久视频综合| 亚洲va日本ⅴa欧美va伊人久久 | 国产熟女午夜一区二区三区| av国产精品久久久久影院| 国产精品一二三区在线看| 日韩欧美一区视频在线观看| 又大又爽又粗| 天堂8中文在线网| 国产视频一区二区在线看| 午夜免费观看性视频| 建设人人有责人人尽责人人享有的| 天天躁狠狠躁夜夜躁狠狠躁| 性色av乱码一区二区三区2| 99热网站在线观看| 热99国产精品久久久久久7| av天堂在线播放| 美女高潮到喷水免费观看| 麻豆国产av国片精品| 久久热在线av| 老熟妇仑乱视频hdxx| 十八禁网站网址无遮挡| bbb黄色大片| 黄色毛片三级朝国网站| 天堂俺去俺来也www色官网| 久久毛片免费看一区二区三区| 日本a在线网址| 国产精品亚洲av一区麻豆| 91九色精品人成在线观看| 日韩制服丝袜自拍偷拍| 国产精品av久久久久免费| 久久这里只有精品19| 国产成人精品久久二区二区免费| 精品亚洲成a人片在线观看| 免费高清在线观看日韩| 各种免费的搞黄视频| 一区二区三区乱码不卡18| av天堂在线播放| 亚洲性夜色夜夜综合| 国产精品秋霞免费鲁丝片| 国产1区2区3区精品| 老熟妇仑乱视频hdxx| 国产欧美亚洲国产| 夜夜夜夜夜久久久久| 欧美日韩国产mv在线观看视频| www.自偷自拍.com| 亚洲欧美精品自产自拍| 亚洲精品国产色婷婷电影| 国产精品1区2区在线观看. | 超碰97精品在线观看| 一本一本久久a久久精品综合妖精| 18在线观看网站| 精品第一国产精品| 欧美+亚洲+日韩+国产| 欧美日韩中文字幕国产精品一区二区三区 | 老熟妇仑乱视频hdxx| 十八禁网站免费在线| 超色免费av| 久久影院123| 黄网站色视频无遮挡免费观看| 久久久欧美国产精品| 又大又爽又粗| 啦啦啦中文免费视频观看日本| 亚洲精品在线美女| 午夜福利免费观看在线| 18禁裸乳无遮挡动漫免费视频| 免费在线观看视频国产中文字幕亚洲 | 亚洲色图 男人天堂 中文字幕| 淫妇啪啪啪对白视频 | 亚洲一码二码三码区别大吗| 精品久久久久久电影网| 亚洲av成人一区二区三| 久久狼人影院| 国产av又大| 亚洲 欧美一区二区三区| 久久久久久久大尺度免费视频| 大香蕉久久网| 大码成人一级视频| 999久久久精品免费观看国产| 亚洲精品粉嫩美女一区| 亚洲精品乱久久久久久| 亚洲欧美激情在线| 90打野战视频偷拍视频| 久久久国产一区二区| 九色亚洲精品在线播放| 国产在线视频一区二区| 大码成人一级视频| 不卡一级毛片| 亚洲av日韩精品久久久久久密| 日本撒尿小便嘘嘘汇集6| 欧美人与性动交α欧美软件| 亚洲第一青青草原| 午夜老司机福利片| 岛国在线观看网站| 亚洲第一青青草原| 亚洲免费av在线视频| www日本在线高清视频| 亚洲精品久久成人aⅴ小说| 国产黄色免费在线视频| 99九九在线精品视频| 国产在线观看jvid| 亚洲人成电影观看| 王馨瑶露胸无遮挡在线观看| 免费少妇av软件| av网站免费在线观看视频| 中文字幕最新亚洲高清| 丝瓜视频免费看黄片| 免费观看人在逋| 多毛熟女@视频| 国产日韩欧美视频二区| 亚洲精品美女久久av网站| 亚洲欧美成人综合另类久久久| 一本综合久久免费| 1024香蕉在线观看| 欧美日韩一级在线毛片| 美女国产高潮福利片在线看| 岛国毛片在线播放| 国产男女内射视频| 亚洲 国产 在线| 中文字幕人妻丝袜制服| 亚洲欧美日韩高清在线视频 | 精品久久久久久久毛片微露脸 | 乱人伦中国视频| 亚洲精品中文字幕一二三四区 | 亚洲熟女毛片儿| 国产在线免费精品| 亚洲欧洲日产国产| www.自偷自拍.com| 最新的欧美精品一区二区| 亚洲精品久久午夜乱码| 啦啦啦免费观看视频1| 十八禁人妻一区二区| av国产精品久久久久影院| 国产亚洲精品第一综合不卡| 嫁个100分男人电影在线观看| 亚洲精品av麻豆狂野| 亚洲国产毛片av蜜桃av| 成人手机av| 热re99久久国产66热| 亚洲欧美色中文字幕在线| 午夜福利在线观看吧| 窝窝影院91人妻| 丝袜喷水一区| 69精品国产乱码久久久| 9色porny在线观看| 国产亚洲精品第一综合不卡| 91国产中文字幕| 国产免费一区二区三区四区乱码| 宅男免费午夜| tocl精华| 在线观看免费高清a一片| 久久国产精品男人的天堂亚洲| 国产成人影院久久av| av国产精品久久久久影院| 免费黄频网站在线观看国产| 如日韩欧美国产精品一区二区三区| 亚洲国产欧美网| 欧美97在线视频| 在线观看免费日韩欧美大片| 中文字幕精品免费在线观看视频| 一个人免费看片子| 大陆偷拍与自拍| 91国产中文字幕| av在线播放精品| 精品亚洲成a人片在线观看| 午夜福利在线观看吧| 午夜福利视频精品| kizo精华| 两人在一起打扑克的视频| 精品少妇一区二区三区视频日本电影| 国产精品自产拍在线观看55亚洲 | 久久久精品国产亚洲av高清涩受| 久久人人爽人人片av| 国产一区二区三区在线臀色熟女 | 十八禁网站网址无遮挡| 国产成人啪精品午夜网站| 国产一级毛片在线| 色精品久久人妻99蜜桃| 国产亚洲精品一区二区www | 国产91精品成人一区二区三区 | 一本大道久久a久久精品| 热99久久久久精品小说推荐| 18禁裸乳无遮挡动漫免费视频| 亚洲中文字幕日韩| 亚洲av欧美aⅴ国产| 在线观看免费高清a一片| 一区二区日韩欧美中文字幕| 亚洲国产av新网站| 不卡一级毛片| av超薄肉色丝袜交足视频| 黄色毛片三级朝国网站| 久久av网站| 各种免费的搞黄视频| 不卡一级毛片| 这个男人来自地球电影免费观看| 岛国在线观看网站| 美女脱内裤让男人舔精品视频| 中文字幕最新亚洲高清| 国产又爽黄色视频| 两人在一起打扑克的视频| 建设人人有责人人尽责人人享有的| 久久中文字幕一级| 国精品久久久久久国模美| 日韩电影二区| 日韩大片免费观看网站| 纯流量卡能插随身wifi吗| 69av精品久久久久久 | 国产深夜福利视频在线观看| 欧美精品一区二区免费开放| 少妇 在线观看| 老司机影院毛片| 91大片在线观看| 国产亚洲精品第一综合不卡| bbb黄色大片| 精品少妇一区二区三区视频日本电影| 国产成人精品久久二区二区91| 一级毛片精品| 亚洲男人天堂网一区| 777米奇影视久久| 久久久久久久国产电影| 91国产中文字幕| 五月天丁香电影| 久久久国产精品麻豆| 正在播放国产对白刺激| 精品视频人人做人人爽| 少妇的丰满在线观看| 999精品在线视频| 伊人亚洲综合成人网| 国产免费一区二区三区四区乱码| 国产在线免费精品| 精品国产乱子伦一区二区三区 | 国产成人av激情在线播放| 亚洲国产看品久久| 精品久久蜜臀av无| 91成人精品电影| 国产色视频综合| 一本大道久久a久久精品| 啦啦啦 在线观看视频| 正在播放国产对白刺激| 欧美黑人精品巨大| 2018国产大陆天天弄谢| av在线老鸭窝| 久久久精品国产亚洲av高清涩受| 亚洲国产成人一精品久久久| 亚洲av国产av综合av卡| 美女国产高潮福利片在线看| 天堂中文最新版在线下载| 黄色怎么调成土黄色| 一区二区日韩欧美中文字幕| 亚洲精品国产一区二区精华液| av天堂在线播放| 青草久久国产| 狂野欧美激情性bbbbbb| 精品福利永久在线观看| 欧美性长视频在线观看| 日韩欧美免费精品| 午夜91福利影院| 人妻一区二区av| 国产成人啪精品午夜网站| 丁香六月天网| 日本vs欧美在线观看视频| 超碰97精品在线观看| 亚洲黑人精品在线| 又黄又粗又硬又大视频| 国产亚洲欧美在线一区二区| 精品国内亚洲2022精品成人 | 久久久久网色| 成在线人永久免费视频| 王馨瑶露胸无遮挡在线观看| 国产三级黄色录像| 亚洲色图综合在线观看| 亚洲国产日韩一区二区| 淫妇啪啪啪对白视频 | 女人精品久久久久毛片| 久久久久久人人人人人| 亚洲人成电影观看| 一二三四在线观看免费中文在| 国产免费现黄频在线看| 亚洲 欧美一区二区三区| 他把我摸到了高潮在线观看 | 不卡av一区二区三区| 午夜福利免费观看在线| 日日摸夜夜添夜夜添小说| 免费在线观看完整版高清| 午夜成年电影在线免费观看| 久久精品aⅴ一区二区三区四区| 狂野欧美激情性xxxx| 欧美+亚洲+日韩+国产| 一区在线观看完整版| 亚洲三区欧美一区| 十八禁人妻一区二区| 国产精品99久久99久久久不卡| 亚洲欧美一区二区三区黑人| 不卡av一区二区三区| 精品国产一区二区三区久久久樱花| 国产成人免费无遮挡视频| 午夜福利,免费看| 国产片内射在线| 精品一品国产午夜福利视频| 国产激情久久老熟女| 亚洲欧美色中文字幕在线| av片东京热男人的天堂| 国产亚洲精品第一综合不卡| 一本—道久久a久久精品蜜桃钙片| 亚洲av男天堂| 777久久人妻少妇嫩草av网站| 蜜桃在线观看..| 麻豆av在线久日| 国产免费视频播放在线视频| 老熟女久久久| 欧美精品一区二区大全| 50天的宝宝边吃奶边哭怎么回事| 国产精品影院久久| 母亲3免费完整高清在线观看| 国产一区二区 视频在线| 精品亚洲乱码少妇综合久久| 男女免费视频国产| 一区二区三区四区激情视频| 巨乳人妻的诱惑在线观看| 成人手机av| 亚洲熟女毛片儿| 亚洲av国产av综合av卡| 天堂俺去俺来也www色官网| 这个男人来自地球电影免费观看| 麻豆乱淫一区二区| 免费高清在线观看视频在线观看| 久久国产精品影院| 国产亚洲精品久久久久5区| xxxhd国产人妻xxx| 欧美日本中文国产一区发布| 啦啦啦在线免费观看视频4| 日韩,欧美,国产一区二区三区| 精品少妇久久久久久888优播| 国产伦理片在线播放av一区| 精品少妇久久久久久888优播| 亚洲成av片中文字幕在线观看| 久久久久久久久免费视频了| 久久久久精品人妻al黑| 99国产精品99久久久久| 黑丝袜美女国产一区| 王馨瑶露胸无遮挡在线观看| 日韩欧美免费精品| 免费在线观看视频国产中文字幕亚洲 | 中文字幕最新亚洲高清| 深夜精品福利| 国产av国产精品国产| 精品少妇黑人巨大在线播放| 国产成人av激情在线播放| 欧美黑人欧美精品刺激| 一区二区三区四区激情视频| 一个人免费在线观看的高清视频 | 91精品三级在线观看| 天堂8中文在线网| 91精品三级在线观看| 久久99热这里只频精品6学生| 一二三四社区在线视频社区8| 免费在线观看日本一区| 国产免费现黄频在线看| 一本大道久久a久久精品| 亚洲欧美精品自产自拍| av免费在线观看网站| 电影成人av| 亚洲欧美日韩另类电影网站| 免费人妻精品一区二区三区视频| 12—13女人毛片做爰片一| 最黄视频免费看| 久久久久久人人人人人| 操美女的视频在线观看| 美女午夜性视频免费| av网站在线播放免费| 成年动漫av网址| 欧美中文综合在线视频| 国产伦理片在线播放av一区| 日日爽夜夜爽网站| 久久中文字幕一级| 下体分泌物呈黄色| 国产黄色免费在线视频| 每晚都被弄得嗷嗷叫到高潮| 性少妇av在线| 日韩电影二区| 欧美人与性动交α欧美精品济南到| 人妻人人澡人人爽人人| 国产在线视频一区二区| 欧美精品av麻豆av| 精品亚洲乱码少妇综合久久| 欧美亚洲 丝袜 人妻 在线| 下体分泌物呈黄色| 90打野战视频偷拍视频| 亚洲第一欧美日韩一区二区三区 | av天堂久久9| 国产日韩欧美在线精品| 免费在线观看黄色视频的| 日韩 亚洲 欧美在线| 亚洲av电影在线进入| 高潮久久久久久久久久久不卡| 精品一品国产午夜福利视频| 日韩一卡2卡3卡4卡2021年| 免费在线观看影片大全网站| 伊人亚洲综合成人网| 99久久99久久久精品蜜桃| 91精品三级在线观看| 久久久久精品国产欧美久久久 | 中亚洲国语对白在线视频| 亚洲va日本ⅴa欧美va伊人久久 | 爱豆传媒免费全集在线观看| 久久久久久人人人人人| 嫁个100分男人电影在线观看| 午夜激情av网站| 免费观看人在逋| 下体分泌物呈黄色| 两性夫妻黄色片| 中文字幕制服av| e午夜精品久久久久久久| 巨乳人妻的诱惑在线观看| www.熟女人妻精品国产| 啦啦啦 在线观看视频| 在线观看舔阴道视频| 欧美xxⅹ黑人| 亚洲熟女毛片儿| 欧美人与性动交α欧美精品济南到| 亚洲av成人不卡在线观看播放网 | 亚洲成人国产一区在线观看| 国产av国产精品国产| 最近中文字幕2019免费版| 99九九在线精品视频| www.av在线官网国产| 黄片播放在线免费| 久久女婷五月综合色啪小说| 亚洲一区中文字幕在线| 男女无遮挡免费网站观看| 国产一区二区三区av在线| 久久久精品区二区三区| 久久久精品国产亚洲av高清涩受| 精品国产一区二区三区四区第35| 国产人伦9x9x在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 美女扒开内裤让男人捅视频| 人妻一区二区av| 丁香六月欧美| 亚洲欧美一区二区三区黑人| 日韩中文字幕欧美一区二区| 精品一区二区三区四区五区乱码| 精品一区二区三卡| 久久国产精品人妻蜜桃| 国产免费现黄频在线看| 欧美精品啪啪一区二区三区 | www日本在线高清视频| 视频在线观看一区二区三区| 12—13女人毛片做爰片一| av片东京热男人的天堂| 在线永久观看黄色视频| 久久国产亚洲av麻豆专区| 91av网站免费观看| 亚洲熟女毛片儿| 精品国产一区二区久久| 欧美精品高潮呻吟av久久| 国产精品免费视频内射| 欧美日韩一级在线毛片| 啦啦啦免费观看视频1| 亚洲av电影在线观看一区二区三区| 亚洲国产日韩一区二区| 免费日韩欧美在线观看| 成年人黄色毛片网站| 日本五十路高清| 多毛熟女@视频| 19禁男女啪啪无遮挡网站| 国产精品久久久久成人av| 午夜福利在线观看吧| 午夜精品国产一区二区电影| av网站在线播放免费| 日本精品一区二区三区蜜桃| 人妻人人澡人人爽人人| 丝袜人妻中文字幕| 久久精品国产a三级三级三级| 免费久久久久久久精品成人欧美视频| 亚洲精品美女久久av网站| 高清欧美精品videossex| 久久久久国内视频| 亚洲全国av大片| 欧美日韩亚洲国产一区二区在线观看 | 中文字幕高清在线视频| 国产福利在线免费观看视频| 婷婷成人精品国产| 久久久久国产一级毛片高清牌| 亚洲国产精品999| 久久国产亚洲av麻豆专区| 国产成人精品久久二区二区91| 国产无遮挡羞羞视频在线观看| 国产97色在线日韩免费| 欧美激情 高清一区二区三区| 精品少妇内射三级| 中亚洲国语对白在线视频| 亚洲全国av大片| 国产精品免费大片| 19禁男女啪啪无遮挡网站| www日本在线高清视频| 欧美日韩视频精品一区| av在线老鸭窝| 我的亚洲天堂| 免费女性裸体啪啪无遮挡网站| 又紧又爽又黄一区二区| 色婷婷av一区二区三区视频| 欧美老熟妇乱子伦牲交| 女人爽到高潮嗷嗷叫在线视频| 免费看十八禁软件| 十八禁高潮呻吟视频| av天堂久久9| 国产av精品麻豆| 亚洲精品成人av观看孕妇| 精品国产一区二区久久| 另类亚洲欧美激情| 亚洲成人手机| 久久天躁狠狠躁夜夜2o2o| 老熟女久久久| av天堂久久9| 亚洲九九香蕉| 亚洲欧美日韩另类电影网站| 后天国语完整版免费观看| 亚洲专区中文字幕在线| 一区二区av电影网| 日韩视频在线欧美| 多毛熟女@视频| 91精品三级在线观看| 午夜日韩欧美国产| 香蕉丝袜av| a 毛片基地| 午夜老司机福利片| 国产不卡av网站在线观看| 午夜91福利影院| 日韩大码丰满熟妇| 老汉色∧v一级毛片| 日本一区二区免费在线视频| 国产野战对白在线观看| 一本色道久久久久久精品综合| 亚洲精品久久成人aⅴ小说| 久久久久久久精品精品| cao死你这个sao货| 国产深夜福利视频在线观看| 搡老乐熟女国产| 免费观看人在逋| 国产成人欧美| 法律面前人人平等表现在哪些方面 | 丰满人妻熟妇乱又伦精品不卡| 午夜福利免费观看在线| 亚洲中文字幕日韩| 91麻豆av在线| 欧美日韩中文字幕国产精品一区二区三区 | 一区二区三区乱码不卡18| 久久性视频一级片| 日韩 欧美 亚洲 中文字幕| 亚洲 欧美一区二区三区| 国产精品 国内视频| 午夜免费成人在线视频| 国产高清国产精品国产三级| 狠狠婷婷综合久久久久久88av| 免费在线观看影片大全网站| 国产亚洲午夜精品一区二区久久| 国产精品自产拍在线观看55亚洲 | 国产精品1区2区在线观看. | 满18在线观看网站| 自线自在国产av| 青春草亚洲视频在线观看| 91av网站免费观看| 老熟妇仑乱视频hdxx| 精品一区在线观看国产| 中文字幕另类日韩欧美亚洲嫩草| 在线观看www视频免费| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲精品美女久久久久99蜜臀| 国产老妇伦熟女老妇高清| 人人妻人人澡人人看| 精品乱码久久久久久99久播| 老司机影院成人| 最黄视频免费看| 蜜桃在线观看..| www日本在线高清视频| 国产成人一区二区三区免费视频网站| 久久中文看片网| 亚洲一区中文字幕在线| 日本撒尿小便嘘嘘汇集6| 亚洲精品美女久久久久99蜜臀| 九色亚洲精品在线播放| 欧美黄色淫秽网站| 日本黄色日本黄色录像| 又紧又爽又黄一区二区| 大码成人一级视频| 亚洲av电影在线进入| 欧美亚洲日本最大视频资源|