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

    Application of response surface methodology for optimization of purge gas recycling to an industrial reactor for conversion of CO2 to methanol

    2017-05-28 10:23:16HadisehKhalilpourmeymandiAzadehMirvakiliMohammadRezaRahimpourAlirezaShariati

    Hadiseh Khalilpourmeymandi,Azadeh Mirvakili,Mohammad Reza Rahimpour ,*,Alireza Shariati

    1 Department of Chemical Engineering,School of Chemical and Petroleum Engineering,Shiraz University,Shiraz 71345,Iran

    2 Department of Chemical Engineering,School of Chemical and Petroleum Engineering,Persian Gulf University,Bushehr 75169,Iran

    3 Shiraz University,Natural Gas Engineering Department,School of Chemical and Petroleum Engineering,Shiraz 71345,Iran

    1.Introduction

    Environmental issues due to emissions of gases induced by petrochemical industries have become worldwide problems.The proportion of gases has increased significantly.Millions of cubic meters of gases are burnt and wasted annually in fuel oil installations in all around the globe[1].This is the beginning of entering millions of tons of carbon dioxide into the atmosphere and therefore the environment is affected by this action.Furthermore,negative changes in climate can highly increase through production of carbon dioxide,black carbon,and the other environmental contaminants.Moreover,there are valuable gases in such flares which are wasted.As a result,they should be utilized in order to lead to a sustainable economic prosperity in developing countries.

    Iran is the third largest country in the world which produces flare gases.The requirement of collecting purge gases and optimum applying of such detrimental gases are quite perceptible due to environmental and economic problems caused by burning and wasting them.In fact,they can be recycled to the unit to pollution reduction and its efficient use.

    The aim of this study is to recycle the purge gas of methanol unit.In order to realize more comprehensive description of methanol production process and modeling,a few mathematical modeling studies on this process are provided.

    There are a vast variety of research projects about the modeling of methanol unit[2–5].In 2004,Rahimpour and Ghader[6]investigated the enhancement of carbon monoxide conversion in a membrane named palladium–silver membrane for methanol synthesis.A novel reactor con figuration with perm-selective Pd–Ag(23 wt%Ag)wall to hydrogen was proposed in their studies.Furthermore,a wide range of important operation conditions were examined.The simulation results demonstrated that carbon monoxide conversion beyond equilibrium value can be enhanced by a membrane reactor.A dynamic model for a membrane dual-type methanol reactor in the presence of long term catalyst deactivation was developed by Rahimpour and Lotfinejad[7]in 2007.The performance of a membrane dual-type methanol reactor in comparison with the conventional one was investigated.Indeed,a favorable temperature pro file and the activity of membrane dual-type reactor relative to single and conventional dual-type reactor systems were shown by simulation results.Accordingly,when a membrane was used in a conventional dual-type methanol reactor,the performance of methanol reactor systems improved.The potential for methanol production from natural gas by direct catalytic partial oxidation was also examined by Edwards and Foster[8]in 1986.

    Purge gas streams of methanol unit contain approximately 63%hydrogen,20%carbon monoxide and carbon dioxide as reactants and 17%nitrogen and methane as inert.For the purpose of recycling the components of purge gas streams,optimal recycle components have been determined by RSM.

    1.1.Response Surface Methodology(RSM)

    The relationships between numerous explanatory variables and one or more response variables are explored by RSM which had been presented by Box and Wilson[9]in 1951.RSM is employed to optimize the recycle components.

    Up to now,different studies have been accomplished on optimization of various processes by RSM[10–14].Samimiet al.[15]has attempted to improve conditions of an industrial methylacetylenepropadiene(MAPD)hydrogenation reactor during the hydrogenation operation with the investigation and modification of different parameters which can directly affect the performance of reactor.The effect of different parameters on the rate of propylene production has been investigated by using RSM.A continuous production of biodiesel from waste cooking oil in a reactive distillation column catalyzed by solid hetero poly acid was done by Noshadiet al.[16]in 2012.RSM based on central composite design(CCD)was utilized to design the experiment and analyzed parameters such as reboiler duty,total feed flow rate,feed temperature and methanol/oil ratio.

    A combination of computational fluid dynamics(CFD)and RSM were utilized by Kumar and Bansal[17]in 2013 in order to model and optimize the photo catalytic degradation of Rhodamine B in an annular photo catalytic reactor.After modeling and simulating the photo catalytic process with CFD,the RSM was utilized for optimization of the process based upon data derived from the CFD simulations.It was indicated that CFD and RSM could be applied for the optimization of photo catalytic processes in order to save time,cost and effort.The catalytic pyrolysis of waste high-density polyethylene to liquid fuel over modified catalysts was optimizedviaRSM by Kumar and Singh[18]in 2014.Countless parameters such as reaction temperature,acidity of modified catalysts and mass ratio between modified catalysts and waste highdensity polyethylene(HDPE)were selected as independent variables.The obtained quadratic model fitted well to predict the response with a high determination coefficient ofR2(0.995).The optimization process with reduced number of costly experiments was successfully achieved by RSM.According to the results which were obtained by RSM,membrane was applied in order to separate gas components.

    1.2.Membrane

    Membrane separation technology is a low-cost process when high purity gas streams are not fundamental[19].Membrane processes for gas separation are gaining a larger acceptance in industries and markets.Adsorption ofCO2,CH4,N2O,and N2on MOF-5,MOF-177,and zeolite 5A was investigated by Sahaetal.[20]in 2010.Both adsorption equilibrium and kinetic data for these gases on all three adsorbents were volumetrically measured at 298 K and gas pressure up to 800 Torr.They understood that both MOF adsorbents have larger adsorption capacities for CO2and CH4in comparison to zeolite 5A at elevated pressure.The influences of gas pressure on diffusivity for various adsorbate–adsorbent systems were also investigated by these authors.Besides,Mfi membranes for carbon dioxide separation were prepared and th

    e separation performance was evaluated by Lindmark and Hedlund[21]in 2010.The permeance of CO2for this membrane was as high as 13 × 10?7mol·m?2·s?1·Pa?1.Their simulation results clearly showed that Mfi membranes were promising candidates for separating carbon dioxide from synthesis gas.Furthermore,Salehiet al.[22]studied the enhancement of carbon monoxide conversion in a novel slurry bubble column reactor for methanol synthesis in 2014.They compared the results of carbon monoxide removal with gas-phase model in the quasi-steady state condition.They finally achieved higher methanol production rate and lower catalyst deactivation rate.Carbon nanotubes(CNTs)with several amines were modified by Suet al.[23]in 2009.These amines included 3-aminopropyl-trimethoxysilyl(APTS),ethylenediamine(EDA)and polyethyleneimine(PEI).The modified CNTs were tested at temperature ranging from 20 to100°C in order to absorb carbon dioxide.The highest capacity of carbon dioxide adsorption was achieved by CNT-APTS at the lowest applied temperature;i.e.55 mg·g?1of CNT-APTS in comparison with 37 mg·g?1of raw CNT.

    1.3.Objective

    There are two purge gas streams in the methanol production unit of domestic petrochemical company that its total flow rate is 1,293,000 mol·h?1.Each purge gas stream contains 63%hydrogen and 20%carbon monoxide and carbon dioxide which are the valuable and beneficial reactants for methanol production reaction.Besides,there is 17%nitrogen and methane as inert gases in the methanol purge gas streams.

    Consequently,the main goal of the current study is to investigate the recycling effect of different components on the methanol production rate.It is undoubted that recycling all components of purge gas leads to production reduction owing to inert gas accumulation.Therefore,the inert gases should be separated from the purge gages and only reactants are recycled to the reactor.A membrane gas separation unit has been considered as a purification unit for separating methane and nitrogen from the purge gas streams.In addition,maximizing methanol production rate and minimizing carbon dioxide production rate have been taken into consideration in this paper.With these purposes in mind,the optimization method of RSM is applied.The optimal value of different parameters like temperature,inlet flow rates of carbon dioxide,carbon monoxide and hydrogen are calculated by RSM.

    2.Process Description

    A schematic diagram of methanol synthesis unit is shown in Fig.1(a).The methanol synthesis unit contains shell and tube heat exchanger which the first reactor(water-cooled reactor)is combined serially with the second one(gas-cooled reactor).Cold feed synthesis gas is preheated inside the tubes of the second reactor(gas-cooled reactor),and subsequently fed into the tube side of the water-cooled reactor where the chemical reaction is initiated by catalysts.Reactions of methanol synthesis are performed over commercial CuO/ZnO/Al2O3catalyst.In this stage,synthesis gas is partly converted into hydrocarbons.The reacting gas leaving the water-cooled reactor is directed into the shell side of the gas-cooled reactor in counter-current mode with synthesis gas flowing through the tube side.The next stage is the purification of the produced methanol.The outlet of reactors is cooled until that methanol is condensed.Purge gases and methanol condensation are separated in a flash drum.The output of the drum unit is divided into two parts:

    ?95%of unconverted syngas is recycled back to the methanol converter to enhance the overall conversion,thereby improving the process economically.

    ?The residual unconverted syngas(approximately 5%)is used as purge gas.

    As clearly observed in Fig.1(b),four branches of purge gas streams are available:

    1-The first stream is zero flaring which means absolute no flaring;closed flare.

    2-The second stream is sent to feed pretreatment as a hydrogen rich stream and mixed with methane gas for feed desulfurization.The flow rate of this stream is 8157 m3·h?1.

    Fig.1.Schematic diagram of(a)a methanol production unit,(b)methanol purge gas streams and(c)purge gas purification unit.

    3-The third stream is sent to the hydrogen recovery by the process of pressure swing adsorption(PSA)which cannot be returned to the reactor.

    4-The fourth steam is sent to the steam reformer and it is used as fuel in the steam reforming unit.The flow rate of this stream is 22,679 m3·h?1.

    There are two returnable purge gas streams in the domestic petrochemical company.One of them is utilized as fuel;however,the other purge gas stream is used as hydrogen rich stream in the desulfurization unit.Although COand CO2are considered as inert gases in the desulfurization unit,they are reactants of methanol synthesis reactors.Therefore,it is suggested that this branch of purge gas stream is recycled to the feed stream.Another stream is used as fuel in the steam reformer.Fuel is much cheaper than purge gas because actually purge gas is an unreacted synthesis gas which is considerably pricey to be produced.As a matter of fact,by recycling all components not only nitrogen and methane as inert gases are accumulated in the reactor during recycling,but also it reduces the reaction rate;moreover,methanol production rate decreases.Therefore,it is necessary to apply a purification unit before recycling these two purge streams in order to remove unwanted components such as nitrogen and methane.

    For the purpose of gas separation,various methods can be introduced which are cryogenic method,absorption with lean oil,adsorption and membrane processes.Membrane processes have been utilized as a purification unit for separating gas components in the current study owing to their simplicity and low maintenance costs.

    2.1.Purification unit

    Hydrogen is the smallest component in the methanol purge gas stream.There is approximately 63%hydrogen in the purge gas.As a result,at first a membrane should be designed for hydrogen separation.There are many membranes for separating hydrogen in the literature.Pd-Ag,a popular membrane for H2separation,is examined in this study.Afterwards,CO and CO2should be separated from the gas mixture;therefore,carbon nanotubes(CNTs)are proffered.

    17%of methanol purge gas contains inert gases like nitrogen and methane.Because of separating such inert gases,a Pd–Ag membrane is applied for hydrogen removal;furthermore,carbon nanotubes(CNTs)are preferred in order to separate carbon dioxide and carbon monoxide.A schematic diagram of this purification unit is illustrated in Fig.1(c).

    Consequently,hydrogen,carbon dioxide and carbon monoxide are recycled to the feed of methanol unit in order to increase the methanol production rate.Besides,methane can be used as fuel in the methanol production unit and nitrogen can be converted into NOxin the steam reforming unit.

    3.Reaction Scheme and Kinetics

    Methanol synthesis includes three overall reactions.They are as follows:

    Reactions(1)–(3)are not independent so that one is a linear combination of the other ones.In the current work,the rate expressions have been selected from Graafet al.[24].The rate equations combining with the equilibrium rate constants provide enough information about kinetics of methanol synthesis[3].The corresponding rate expressions due to the hydrogenation of CO,CO2and reversed water–gas shift reactions over commercial CuO/ZnO/Al2O3catalysts are:

    The reaction rate constants,adsorption equilibrium constants and reaction equilibrium constants in the formulation of kinetic expressions have been tabulated in Table 1,respectively.

    Table 1The reaction rate constants,the adsorption equilibrium constants,and the reaction equilibrium constants for methanol synthesis

    4.Mathematical Model

    4.1.Fixed bed reactors

    A one-dimensional heterogeneous model,which is a conventional model for a catalytic reactor with heat and mass transfer resistances,has been developed for this reactor in order to determine concentration and temperature distributions inside the reactor.The reactor is simulated by MATLAB software.In this model the following assumptions are used:

    The gas mixture is an ideal gas in both catalytic reactors[6,7].

    Both reactors are operated at steady-state conditions.

    Radial variations in both beds are negligible(one-dimensional model).

    Axial diffusion of mass and heat are negligible.

    Bed porosity in axial and radial directions is constant.

    A one-dimensional plug flow in shell and tube sides is considered.

    Chemical reactions are assumed to take place only in the catalyst particles.

    Heat loss to surrounding is neglected.

    4.1.1.Solid phase

    The mass and energy balances for the solid phase are obtained as follows:whereyssandTsare the solid-phase mole fraction and temperature for components(6 components)iin the reactor and η is the effectiveness factor.

    4.1.2.Fluid phase

    The following mass and energy balance equations are developed for the fluid phase:

    whereyigandTgare fluid phase mole fraction and temperature for componentiin the reactor,respectively.The energy equation consists of heat transfer by convection and exchanged heat between the fluid phase and solid phase.Boundary conditions are as follows:

    whereyi0,gT0gandP0gare fluid phase mole fraction,inlet temperature of the reactor,and pressure at the entrance of the reactor,respectively.

    4.2.Puri fi cation unit

    The flux of hydrogen permeating through Pd–Ag membrane is as follows:

    pH2demonstrates the hydrogen partial pressure in Pa.The outer and inner diameters of the Pd/Ag layer are presented by D0and Diin turn.The pre-exponential factorP0above 200°C is reported as 6.33 × 10?8mol·m?2·s?1·Pa?0.5and the activation energy,Ep,is 15.7 kJ·mol?1[25].

    Furthermore,the permeating flux through CNT membrane has been shown in Eq.14.

    where,α is 0.16 barrel for carbon dioxide permeation and the permeation rate constant for carbon monoxide has been assumed approximately near the permeation rate constant of carbon dioxide.

    5.Model Validation

    The formulated model composed of11 ordinary differential equations and the associated boundary conditions lends itself to be an initial value problem.The algebraic equations in the model incorporate the initial conditions,the reaction rates,the ideal gas assumption,as well as aforementioned correlations for the heat and mass transfer coefficients and the physical properties of fluids.These equations along with the discretized ordinary differential equations using backward finite difference form a set of non-linear algebraic equations.The reactor length is then divided into 100 separate sections and the Gauss–Newton method in MATLAB programming environment is used to solve the non-linear algebraic equations in each section.Following this,the results are obtained and validated.The steady state model validation is performed between the pilot plant data reported by domestic petrochemical company and the mathematical modeling of CRC.Results of CRC model and the corresponding observed data of the pilot plant are presented in Table 2.A good agreement is observed between the modeling results and the pilot plant data.Therefore,this mathematical model performs well under industrial conditions.

    Table 2Comparison of mole percent between model results with plant data for fresh catalysts in CRC

    6.Results

    The aim of this study is to recycle purge streams into the methanol production reactor in order that the methanol production increases.According to the investigations,recycling all components from purge streams is not appropriate and effective.At first,the in fluence of recycling different gases from purge streams to the feed of methanol production unit has been studied.

    1.recycling CO,CO2,H2,CH4and N2from purge streams;

    2.recycling just CO,CO2and H2from purge streams;

    3.recycling just CO and H2from purge streams;

    4.recycling just CO and CO2from purge streams;

    5.recycling just CO2and H2from purge streams;

    6.recycling just CO from purge streams;

    7.recycling just CO2from purge streams;

    8.and recycling just H2from purge streams.

    Despite recycling just H2and COis the bestcon figuration,there is almost no membrane to separate CO2from CO.Therefore,an appropriate purification unit has been performed in order that CO2,H2and CO are separated from the other components in the purge streams which is named as Desired Recycle Con figuration(DRC).The gas separation method is described in the process description section.In this study,the modeling results of DRC and CRC have been compared with each other.

    Fig.2 illustrates methanol flow rate according to different recycling con figurations.As shown,the maximum methanol flow rate is achieved in the third con figuration which is related to recycling just COand H2to the feed stream.As described before,however it is the best con figuration,due to limitation in separation technology,recycling of the mixture of CO,CO2and H2is investigated and it is named as DRC.

    Methanol production rate increases by approximately 0.106%in the second con figuration compared to CRC.As mentioned before,the second con figuration describes returning CO,CO2and H2from two purge gas streams to the methanol feed stream named as DRC.

    Fig.3 shows methanol flow rate as a function of reactor length.As obviously, firstly methanol flow rate rises sharply during the first reactor length and then it increases smoothly during the second reactor length which is owing to high reaction rate in the beginning of reaction.It is clearly seen that the flow rate of methanol increases when DRC is applied.In other words,when DRC is applied,the flow rate of methanol is approximately increased by 0.106%in comparison with the one in CRC.Methanol flow rate enhancement demonstrates not only the superiority of DRC to CRC,but also the increase in profit obtained from DRC.

    Hydrogen is one of the main reactants for methanol production.A comparison of hydrogen flow rate in DRC and CRC has been proposed in Fig.4.Generally hydrogen flow rate decreases along the reactor length owing to its fast consumption during the reactions.The flow rate of hydrogen increases considerably in DRC in comparison with the one in CRC.It is due to this fact that there is about 63%hydrogen in each purge stream which is recycled to the reactor feed stream.Obviously,the flow rate of hydrogen in DRC increases by nearly 0.39%compared to the one in CRC.A further advantage of applying DRC is the reduction of catalyst cocking caused by increasing hydrogen flow rate.

    Molar flow rate of methanol,carbon dioxide,carbon monoxide,water and hydrogen at the end of the reactor is depicted in Fig.5(a)–(c)for CRC and DRC.As clearly seen,the proportion of produced methanol in DRC is higher than the one in CRC.

    Fig.2.Schematic diagram of different con figurations comparing methanol flow rate in the outlet.

    Fig.3.Comparison of methanol flow rate as a function of reactor length in CRC and DRC.

    Fig.4.Comparison of hydrogen flow rate as a function of reactor length in CRC and DRC

    The percentage of different component's conversion such as hydrogen,carbon dioxide and carbon monoxide versus the reactor length for both CRC and DRC are demonstrated in Fig.6(a)–(c),respectively.In Fig.6(a),hydrogen conversion increases dramatically in the first two meters of the reactor and afterward it rises gradually.It can be seen from Fig.6(b)that there is a dramatic rise in the conversion of CO2along the reactor length to approximately 12.5%in the first four meters of the reactor.Subsequently,having leveled off until the second 4 m of the reactor,CO2conversion then increases gradually.Moreover,it is obvious from Fig.6(c)that reactor the conversion of CO in the first grows sharply in comparison to the second reactor.Generally,the percentage of these component's conversions has risen in DRC compared to the one in CRC.It occurs because the consumption of hydrogen,carbon dioxide and carbon monoxide in DRC is higher than the one in CRC.

    By applying DRC,the effect of recycling purge streams to the reactor on temperature along the reactor length is negligible.

    The money obtained from methanol sale for DRC and CRC is depicted in Fig.7.The difference between the money which has been obtained from methanol sale in DRC and CRC is approximately 3763.649$·day?1.As clearly shown,the money obtained from methanol sale in DRD is about1.001 times more than the one in CRC.As a matter of fact,obtaining more money from methanol sale in DRC is one of the advantages of applying this novel con figuration which is profitable.

    Fig.8(a)and(b)gives information about the effect of recycling hydrogen and carbon monoxide to the reactor on the methanol production along the reactor length,respectively.It is depicted that as recycling hydrogen and carbon monoxide to the reactor increases,the methanol production enhances as well.Moreover,these figures show that hydrogen and carbon monoxide are the beneficial components which are highly effective in increasing the methanol production.

    The influence of recycling hydrogen on the rate of carbon dioxide production along the reactor length is illustrated in Fig.9.As recycling hydrogen from purge streams to the feed of reactor increases,the flow rate of carbon dioxide reduces.This is the superiority of recycling much more proportion of hydrogen from purge streams to the reactor feed stream.

    7.Optimization

    In the simplest case,an optimization problem consists of maximizing or minimizing a real function by systematically choosing input values from within an allowed set and computing the value of the function.The generalization of optimization theory and techniques to other formulations comprises a large area of applied mathematics.It is imperative to optimize the industrial operations,manufacturing and engineering activities in order to be more efficient.A mathematical model,which predicts the process behavior,is required to optimize the chemical and related processes[26].

    Fig.5.Comparison of molar flow rate of(a)CH3OH and CO2,(b)CO and H2O(c)H2 in the outlet of CRC and DRC.

    Fig.6.Comparison of(a)hydrogen,(b)carbon dioxide and(c)carbon monoxide conversion in CRC and DRC.

    Fig.7.Comparison of money obtained from methanol sale per day in CRC and DRC.

    Fig.8.Effect of(a)hydrogen mole fraction and(b)carbon monoxide mole fraction on methanol production rate and reactor length.

    Fig.9.Effect of hydrogen mole fraction on carbon dioxide production rate and reactor length.

    RSM is utilized in order to optimize the methanol feed stream and operating conditions in the current study.The independent variables in this research are temperature(A),carbon dioxide flow rate(B),carbon monoxide flow rate(C)and hydrogen flow rate(D)which are given in Table 3.Furthermore,maximizing methanol flow rate and minimizing carbon dioxide flow rate have been considered as an objective function.

    Table 3Domain of independent variables in RSM design

    Design of experiment is based on levels of the independent variables resulted in 273 simulation runs.Eqs.15 and 16 present the fitted model that is able to provide a good fitness of the measured methanol and carbon dioxide production data,respectively.As clearly seen,the first parameter in these equations is the intercept and there are various coefficients which represent the effect of each factor as well as factor interactions.

    Predicted values versus observed values for methanol and carbon dioxide are plotted in Fig.10(a)and(b),respectively.As obviously,the simulated values are in agreement with the predicted ones.Moreover,the values which are not simulated by the simulator can be predicted by Eqs.(15)and(16).

    Internally studentized residuals as a function of fitted values for methanol and carbon dioxide are illustrated in Fig.11(a)and(b),respectively.No unusual structure is observed in these figures.These results show that the proposed model is sufficient.

    Analysis of variance(ANOVA)has been applied for each parameter to identify the significance of the effects and interactions between the parameters.Moreover,a p-value less than 0.05 has been statistically considered significant.Design Expert(version 7)has been utilized for ANOVA and regression coefficient calculation.The important terms for the production of methanol and carbon dioxide are tabulated in Tables 4 and 5,respectively.These tables propose that the domain for all variables,which are considered in this model,present maximum methanol production and minimum carbon dioxide production with no significant change in the output and no change in any other variables.

    Fig.10.Normal probability plots of residual for(a)methanol and(b)carbon dioxide.

    Fig.11.Plots of residual vs. fitted values for(a)methanol and(b)carbon dioxide.

    Table 4Model coefficients and the statistical interference of the effect of single variables and their interactions on the mean value of the estimated methanol production rate

    The simultaneous effects of temperature and carbon dioxide flow rate on the production rates of methanol and carbon dioxide are depicted in Fig.12(a)and(b),respectively.Fig.12(a)shows that at first as the temperature increases,methanol production grows slightly and then it decreases.Subsequently,methanol production remains stable with increasing temperature at the end of the reactor.However,the inlet flow rate of carbon dioxide has no effect on the production of methanol.Fig.12(b)illustrates that simultaneous increase in temperature and inlet flow rate of carbon dioxide leads to increasing in the produced carbon dioxide.Therefore,it can be concluded that the effect of temperature on the production of methanoland carbon dioxide is greater than the input amount of carbon dioxide.

    Fig.13(a)and(b)shows the simultaneous effects of hydrogen flow rate and carbon monoxide flow rate on the methanol production rate and carbon dioxide production rate,separately.Fig.13(a)depicts that the increase in the input flow rate of both hydrogen and carbon monoxide causes the increase in the rate of methanol production.But what stands out from the figure is that the impact of the input amount of hydrogen on the methanol production rate is higher than the input amount of carbon monoxide.Fig.13(b)shows that with increasing the input flow rate of hydrogen,carbon dioxide production rises and with the increase in input flow rate of carbon monoxide,carbon dioxide production goes up slightly.Consequently,the inlet amount of hydrogen has a significant in fluence on the production of methanol and carbon dioxide.

    Table 5Model coefficients and the statistical interference of the effect of single variables and their interactions on the mean value of the estimated carbon dioxide production rate

    Fig.12.Simultaneous effects oftemperature and carbon dioxide flowrate on(a)methanol production and(b)carbon dioxide production.

    Fig.13.Simultaneous effects of hydrogen and carbon monoxide flow rates on(a)methanol production and(b)carbon dioxide production.

    Optimized factors which are obtained by RSM for achieving maximum methanol flow rate and minimum carbon dioxide flow rate have been proposed in Table 6.Maximum methanol flow rate of 2083.61 mol·s?1and minimum carbon dioxide flow rate of 1870.46 mol·s?1are achieved when the inlet temperature of reactor is 420.502 K and the inlet flow rates of carbon dioxide,carbon monoxide and hydrogen are 0.01789,43.5661 and 815.11 mol·s?1,in turn.

    Table 6A set of optimized factors obtained by RSM for achieving maximum methanol flow rate and minimum carbon dioxide flow rate

    8.Conclusions

    The purge gas of methanol unit in Iran petrochemical companies is used in two ways which do not have impact on the increase of methanol production rate.However,these purge gas streams contain the unreacted gases of methanol reactor feed stream.Furthermore,a large amount of money has been spent for producing them in the reforming unit.Therefore,the valuable components of such purge gas streams can be utilized as a feed of methanol synthesis rector.There are different con figurations which can be applied for recycling the components of purge streams to the feed stream.The Desired Recycle Con figuration(DRC)has been investigated in this study.Carbon monoxide,carbon dioxide and hydrogen are the optimal components which can be recycled from purge streams to the methanol feed stream for the purpose of achieving maximum methanol flow rate and minimum carbon dioxide flow rate in the outlet of the reactors.For the reduction of accumulation of nitrogen and methane in the reactor,it is suggested that they are separated from the gas mixture,because the presence of these two inert components in the reactor leads to the reduction of methanol production rate.As a result,a separation unit has been proposed in order to separate these components from the other components in the gas mixture of purge streams.The separation unit includes Pd–Ag and carbon nanotube membranes.By recycling carbon monoxide,carbon dioxide and hydrogen,the methanol flow rate increases approximately by 0.106%in DRC compared to CRC.Furthermore,the difference between the cost obtained from the methanol sale in DRC and CRC is about 3763.649USD·d?1.These results show the superiority of DRC to recycling the other components.In addition,the temperature and the molar flow rates of inlet components such as carbon dioxide,carbon monoxide and hydrogen are optimized by RSM which is one of the best and efficient optimization methods.Obtaining maximum value of methanol flow rate as a desired product and also minimum value of carbon dioxide flow rate as a destructive and detrimental component for the environment are the objective functions of RSM.As a matter of fact,10.83%increase in the methanol production rate and 4.078%decrease in the carbon dioxide production rate have been obtained by optimization.

    Nomenclature

    Accross section area of each tube,m2

    avspecific surface area of catalyst pellet,m2·m?3

    Cpspecific heat of gas at constant pressure,J·mol?1·K?1

    ctotal concentration,mol·m?3

    Didiameter of membrane,m

    Epactivation energy of permeability,kJ·kmol?1

    FTotal molar flow rate,mol·s?1

    fifugacity of componenti,Pa

    ΔHf,ienthalpy of formation of componentiJ·mol?1

    ΔHienthalpy of reaction

    hfgas–solid heat transfer coefficient,W·m?2·K?1

    JCNTpermeation rate through the CNT membrane,mol·m?1·s?1

    JPd-Agpermeation rate of hydrogen through the Pd-Ag membrane,mol·m?1·s?1

    Kiadsorption equilibrium constant of speciesiMPa?1

    KP,jequilibrium constant of reactionj

    kgmass transfer coefficient for componenti,m·s?1

    kireaction rate constant of speciesi

    Lreactor length,m

    Ptotal pressure(for exothermic side:Pa;for endothermic side:Pa

    P0the pre-exponential factor,mol·m?2·s?1·Pa?0.5

    ppartial pressure of componenti,Pa

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

    rireaction rate of componenti,mol·kg?1·s?1

    Ttemperature,K

    Uoverall heat transfer coefficient between exothermic and endothermic sides,W·m?2·K?1

    yimole fraction of componenti,mol·mol?1

    zaxial reactor coordinate,m

    α permeation rate constant,mol·m?1·s?1·Pa?1/2

    ρ density of fluid phase,kg·m?3

    η catalyst effectiveness factor

    π 3.1416

    Superscripts

    g In the bulk of gas phase

    s At the surface of catalyst

    Subscripts

    0 inlet conditions

    ichemical species

    [1]M.R.Rahimpour,S.M.Jokar,Feasibility of flare gas reformation to practical energy in Farashband gas re finery:No gas flaring,J.Hazard.Mater.209-210(2012)204–217.

    [2]M.Farniaei,M.Abbasi,A.Rasoolzadeh,M.R.Rahimpour,Enhancement of methanol,DME and hydrogen production via employing hydrogen perm-selective membranes in a novel integrated thermally double-coupled two-membrane reactor,J.Nat.Gas Sci.Eng.14(2013)158–173.

    [3]G.H.Graaf,P.J.J.M.Sijtsema,E.J.Stamhuis,G.E.H.Joosten,Chemical equilibrium in methanol synthesis,Chem.Eng.Sci.J.41(11)(1986)2883–2890.

    [4]N.Kuczynski,M.H.Oyevaar,R.T.Piters,K.R.Westerterp,Methanol synthesis in a counter-current gas–solid–solid trickle flow reactor:an experimental study,Chem.Eng.Sci.J.42(8)(1987)1887–1889.

    [5]N.Park,M.J.Park,Y.J.Lee,K.S.Ha,K.W.Jun,Kinetic modeling of methanol synthesis over commercial catalysts based on three-site adsorption,Fuel Process.Technol.125(2014)139–147.

    [6]M.R.Rahimpour,S.Ghader,Enhancement of CO conversion in a novel Pd–Ag membrane reactor for methanol synthesis,Chem.Eng.Process.43(9)(2004)1181–1188.

    [7]M.R.Rahimpour,M.Lot finejad,Enhancement of methanol production in a membrane dual-type reactor,Chem.Eng.Technol.30(8)(2007)1062–1076.

    [8]J.H.Edwards,N.R.Foster,The potential for methanol production from natural gas by direct catalytic partial oxidation,Fuel Sci.Technol.Int.4(4)(1986)365–390.

    [9]G.E.P.Box,K.B.Wilson,On the experimental attainment of optimum conditions(with discussion),J.R.Stat.Soc.Ser.B13(1)(1951)1–45.

    [10]C.Songa,Y.Kitamurab,S.Lib,Optimization of a novel cryogenic CO2capture process by response surface methodology(RSM),J.Taiwan Inst.Chem.Eng.45(4)(2014)1666–1676.

    [11]S.Toemen,W.A.Wan Abu Bakar,R.Ali,Investigation of Ru/Mn/Ce/Al2O3catalyst for carbon dioxide methanation:Catalytic optimization,physicochemical studies and RSM,J.Taiwan Inst.Chem.Eng.45(5)(2014)2370–2378.

    [12]G.F.Silva,F.L.Camargo,A.L.O.Ferreira,Application of response surface methodology for optimization of biodiesel production by transesterification of soybean oil with ethanol,Fuel Process.Technol.92(3)(2011)3407–3413.

    [13]Y.Wu,S.Zhou,F.Qin,X.Ye,K.Zheng,Modeling physical and oxidative removal properties of Fenton process for treatment of land fill leachate using response surface methodology(RSM),J.Hazard.Mater.180(1–3)(2010)456–465.

    [14]X.Duan,Z.Zhang,C.Srinivasakannan,F.Wang,J.Liang,Regeneration of spent catalyst from vinyl acetate synthesis as porous carbon:Process optimization using RSM,Chem.Eng.Res.Des.92(7)(2014)1249–1256.

    [15]F.Samimi,Z.Khadem Modarresi,O.Dehghani,M.R.Rahimpour,A.Bolhasani,Application of response surface methodology for optimization of an industrial methylacetylene and propadiene hydrogenation reactor,J.Taiwan Inst.Chem.Eng.46(2015)51–64.

    [16]I.Noshadi,N.A.S.Amin,S.Parnas,Continuous production of biodiesel from waste cooking oil in a reactive distillation column catalyzed by solid hetero poly acid:Optimization using response surface methodology(RSM),Fuel94(2012)156–164.

    [17]J.Kumar,A.Bansal,Photo catalytic degradation in annular reactor:Modelization and optimization using computational fluid dynamics(CFD)and response surface methodology(RSM),J.Environ.Chem.Eng.1(2013)398–405.

    [18]S.Kumar,R.K.Singh,Optimization of process parameters by response surface methodology(RSM)for catalytic pyrolysis of waste high-density polyethylene to liquid fuel,J.Environ.Chem.Eng.2(2014)115–122.

    [19]M.Takht Ravanchi,T.Kaghazchi,A.Kargari,Application of membrane separation processes in petrochemical industry:A review,Desalination235(1–3)(2009)199–244.

    [20]D.Saha,Z.Bao,F.Jia,S.Deng,Adsorption ofCO2,CH4,N2O,and N2on MOF-5,MOF-177,and zeolite 5A,Environ.Sci.Technol.44(5)(2010)1820–1826.

    [21]J.Lindmark,J.Hedlund,Carbon dioxide removal from synthesis gas using Mfi membranes,J.Membr.Sci.360(1–2)(2010)284–291.

    [22]K.Salehi,S.M.Jokar,J.Shariati,M.Bahmani,M.A.Sedghamiz,M.R.Rahimpour,Enhancement of CO conversion in a novel slurry bubble column reactor for methanol synthesis,J.Nat.Gas Sci.Eng.21(2014)170–183.

    [23]F.Su,C.Lu,W.Cnen,H.Bai,J.F.Hwang,Capture of CO2from flue gas via multi-walled carbon nanotubes,Sci.Total Environ.407(8)(2009)3017–3023.

    [24]G.H.Graaf,H.Scholtens,E.J.Stamhuis,Intra-particle diffusion limitations in lowpressure methanol synthesis,Chem.Eng.Sci.45(4)(1990)773–783.

    [25]M.Bayat,M.R.Rahimpour,Simultaneous production of ultrapure hydrogen and gasoline in a novel thermally coupled double membrane reactor,Int.J.Energy Res.37(1)(2013)35–48.

    [26]G.P.Rangaiah,Multi-Objective Optimization:Techniques and Applications in Chemical Engineering,National University of Singapore,Singapore,2008.

    国内揄拍国产精品人妻在线| 天堂网av新在线| 一级av片app| 亚洲成a人片在线一区二区| 欧美色欧美亚洲另类二区| 欧美3d第一页| 久久草成人影院| 女的被弄到高潮叫床怎么办| 久久草成人影院| 午夜精品一区二区三区免费看| 搞女人的毛片| 国产三级在线视频| 日本爱情动作片www.在线观看| 在线免费观看不下载黄p国产| 国内精品久久久久精免费| 在线观看66精品国产| 欧美日韩精品成人综合77777| 亚洲av一区综合| 99热全是精品| 91狼人影院| 日产精品乱码卡一卡2卡三| 国产av不卡久久| 夜夜爽天天搞| 亚洲人成网站在线播| 国产探花极品一区二区| 欧美bdsm另类| 欧美激情在线99| 久久草成人影院| 国产免费男女视频| 成人国产麻豆网| 久久6这里有精品| 国产单亲对白刺激| 特级一级黄色大片| 午夜a级毛片| 亚洲七黄色美女视频| 色视频www国产| 午夜免费男女啪啪视频观看| 亚洲成人精品中文字幕电影| 国产高清有码在线观看视频| 国产精品蜜桃在线观看 | 在线观看免费视频日本深夜| 国产精品免费一区二区三区在线| 欧美色视频一区免费| 国产成人a区在线观看| 日本成人三级电影网站| 国产精品一及| 久久九九热精品免费| 中文在线观看免费www的网站| 网址你懂的国产日韩在线| 91aial.com中文字幕在线观看| 亚洲精品国产成人久久av| 一进一出抽搐动态| 久久久久免费精品人妻一区二区| 只有这里有精品99| 亚洲av熟女| 欧美成人免费av一区二区三区| 久久久久久久久久成人| 18禁黄网站禁片免费观看直播| 午夜亚洲福利在线播放| 高清毛片免费看| 老熟妇乱子伦视频在线观看| 丰满人妻一区二区三区视频av| 国产成人影院久久av| 观看免费一级毛片| 国产乱人视频| 精品免费久久久久久久清纯| 在线播放国产精品三级| 国产高潮美女av| АⅤ资源中文在线天堂| 人妻系列 视频| 久久这里有精品视频免费| 精品熟女少妇av免费看| 色综合站精品国产| 国产v大片淫在线免费观看| 亚洲最大成人中文| 99热网站在线观看| 国产亚洲精品久久久com| 精品人妻熟女av久视频| 午夜视频国产福利| 国产真实伦视频高清在线观看| 欧美日韩乱码在线| 99热这里只有精品一区| 亚洲中文字幕一区二区三区有码在线看| 直男gayav资源| 国产高清不卡午夜福利| 久久人人爽人人爽人人片va| 99热6这里只有精品| 大又大粗又爽又黄少妇毛片口| 99久久久亚洲精品蜜臀av| 国产一区二区在线av高清观看| 国产精品美女特级片免费视频播放器| 亚洲av中文字字幕乱码综合| 精品人妻视频免费看| 69人妻影院| 欧美另类亚洲清纯唯美| 欧美三级亚洲精品| 最近视频中文字幕2019在线8| 日本撒尿小便嘘嘘汇集6| 国产成人一区二区在线| 亚洲精品日韩av片在线观看| 免费搜索国产男女视频| 亚洲av成人精品一区久久| 99久久久亚洲精品蜜臀av| 久久久精品94久久精品| 国产精品蜜桃在线观看 | 一级二级三级毛片免费看| 干丝袜人妻中文字幕| 少妇熟女aⅴ在线视频| 久久这里只有精品中国| 久久久色成人| 国产精品人妻久久久久久| 亚洲激情五月婷婷啪啪| 97超视频在线观看视频| 乱人视频在线观看| 狠狠狠狠99中文字幕| 最近手机中文字幕大全| 国产精品伦人一区二区| 国产在视频线在精品| 看黄色毛片网站| 九九在线视频观看精品| 2021天堂中文幕一二区在线观| 99热网站在线观看| 国产黄色小视频在线观看| 久久久国产成人免费| 欧美高清性xxxxhd video| 人人妻人人看人人澡| 寂寞人妻少妇视频99o| 热99在线观看视频| 中国国产av一级| 色视频www国产| 99久久九九国产精品国产免费| 亚洲av.av天堂| 好男人在线观看高清免费视频| 麻豆一二三区av精品| 只有这里有精品99| 亚洲国产精品成人综合色| av黄色大香蕉| av在线观看视频网站免费| 69av精品久久久久久| 日韩欧美国产在线观看| 国产精品久久久久久亚洲av鲁大| 亚洲婷婷狠狠爱综合网| 午夜免费男女啪啪视频观看| 久久精品国产亚洲av香蕉五月| 精品欧美国产一区二区三| 97超碰精品成人国产| eeuss影院久久| 性色avwww在线观看| 99久久精品一区二区三区| 黄片无遮挡物在线观看| 美女cb高潮喷水在线观看| 国产av麻豆久久久久久久| 亚洲乱码一区二区免费版| 91久久精品国产一区二区成人| 国产又黄又爽又无遮挡在线| av在线观看视频网站免费| 91久久精品国产一区二区成人| 久久久久久久亚洲中文字幕| 波多野结衣高清无吗| 两个人视频免费观看高清| 中文字幕人妻熟人妻熟丝袜美| 悠悠久久av| 国产精品女同一区二区软件| 精品久久久噜噜| 亚洲av中文av极速乱| 91在线精品国自产拍蜜月| 日韩欧美国产在线观看| 日韩精品青青久久久久久| 国产高清有码在线观看视频| av专区在线播放| 少妇人妻一区二区三区视频| 日韩亚洲欧美综合| 国产午夜精品论理片| 久久久成人免费电影| 亚洲国产色片| 久久久久久久久久成人| 床上黄色一级片| 国产一级毛片在线| 99久久久亚洲精品蜜臀av| 一夜夜www| 最近最新中文字幕大全电影3| 日韩人妻高清精品专区| 99riav亚洲国产免费| 午夜免费男女啪啪视频观看| 亚洲综合色惰| 婷婷色av中文字幕| 久久久久久九九精品二区国产| 亚洲av第一区精品v没综合| 久久人妻av系列| 国产亚洲5aaaaa淫片| 久久午夜福利片| 国内精品久久久久精免费| 少妇人妻一区二区三区视频| 99久久精品国产国产毛片| 国产精品不卡视频一区二区| 亚洲av成人av| 欧洲精品卡2卡3卡4卡5卡区| 日产精品乱码卡一卡2卡三| 亚洲精品久久国产高清桃花| kizo精华| 99在线人妻在线中文字幕| 男人狂女人下面高潮的视频| 永久网站在线| 观看美女的网站| 午夜久久久久精精品| 免费看a级黄色片| 亚洲精品久久国产高清桃花| 久久久久国产网址| 欧美性猛交╳xxx乱大交人| 亚洲天堂国产精品一区在线| 少妇被粗大猛烈的视频| 99热6这里只有精品| 日本成人三级电影网站| 一区二区三区四区激情视频 | 亚洲欧洲日产国产| 级片在线观看| 在线观看一区二区三区| 国产视频内射| 麻豆久久精品国产亚洲av| 狠狠狠狠99中文字幕| 99热精品在线国产| 亚洲一级一片aⅴ在线观看| 久久精品国产99精品国产亚洲性色| 一级黄色大片毛片| 日本爱情动作片www.在线观看| 亚洲av成人精品一区久久| 免费观看a级毛片全部| 1024手机看黄色片| 成人亚洲欧美一区二区av| 青春草亚洲视频在线观看| 久久久久久久久久久丰满| 国产av在哪里看| 亚洲欧美清纯卡通| 看十八女毛片水多多多| 噜噜噜噜噜久久久久久91| 亚洲天堂国产精品一区在线| 在线观看一区二区三区| 亚洲一区高清亚洲精品| 亚洲精品日韩在线中文字幕 | 亚洲精品乱码久久久久久按摩| 性插视频无遮挡在线免费观看| 精品久久久久久久人妻蜜臀av| 久久久成人免费电影| 国产精品一二三区在线看| 你懂的网址亚洲精品在线观看 | 一边摸一边抽搐一进一小说| 国产av麻豆久久久久久久| eeuss影院久久| 久久国内精品自在自线图片| 午夜视频国产福利| av免费观看日本| av卡一久久| 日本在线视频免费播放| 久久欧美精品欧美久久欧美| 美女cb高潮喷水在线观看| 天天一区二区日本电影三级| 深爱激情五月婷婷| 成人午夜精彩视频在线观看| 国产成年人精品一区二区| 欧美极品一区二区三区四区| 日本与韩国留学比较| 国产淫片久久久久久久久| 两性午夜刺激爽爽歪歪视频在线观看| 国产一区亚洲一区在线观看| 中出人妻视频一区二区| 欧美+亚洲+日韩+国产| 禁无遮挡网站| 插逼视频在线观看| 最近中文字幕高清免费大全6| a级一级毛片免费在线观看| 久久精品国产清高在天天线| 久久精品国产亚洲av天美| 国产精品福利在线免费观看| 欧美色视频一区免费| 婷婷精品国产亚洲av| 偷拍熟女少妇极品色| av专区在线播放| 国产精品免费一区二区三区在线| 国产av一区在线观看免费| 美女大奶头视频| 国产美女午夜福利| 国产爱豆传媒在线观看| 日韩亚洲欧美综合| 国产极品精品免费视频能看的| 亚洲aⅴ乱码一区二区在线播放| 午夜免费男女啪啪视频观看| 欧美一区二区国产精品久久精品| 成年av动漫网址| 欧美日韩一区二区视频在线观看视频在线 | 性色avwww在线观看| 老女人水多毛片| 国产免费一级a男人的天堂| a级毛片免费高清观看在线播放| 三级毛片av免费| 国产午夜精品论理片| 97人妻精品一区二区三区麻豆| 毛片一级片免费看久久久久| 精品久久久久久久久久免费视频| 国产免费一级a男人的天堂| a级毛片a级免费在线| 美女高潮的动态| 日本成人三级电影网站| 亚洲内射少妇av| 免费看av在线观看网站| 婷婷亚洲欧美| 欧美在线一区亚洲| 国产免费男女视频| 国产黄色小视频在线观看| 欧美高清成人免费视频www| 国产亚洲91精品色在线| 成人综合一区亚洲| 黑人高潮一二区| 国内久久婷婷六月综合欲色啪| 婷婷色综合大香蕉| 久久草成人影院| 久久99精品国语久久久| 日本撒尿小便嘘嘘汇集6| 22中文网久久字幕| 成人性生交大片免费视频hd| 边亲边吃奶的免费视频| 亚洲国产精品成人久久小说 | 国产视频内射| 又粗又爽又猛毛片免费看| 老司机影院成人| 夫妻性生交免费视频一级片| 97热精品久久久久久| 黄色配什么色好看| а√天堂www在线а√下载| 尾随美女入室| 久久综合国产亚洲精品| 内射极品少妇av片p| 噜噜噜噜噜久久久久久91| 免费一级毛片在线播放高清视频| 久久久久久久午夜电影| 亚洲人成网站在线观看播放| 成人特级av手机在线观看| 可以在线观看毛片的网站| 欧美精品国产亚洲| 啦啦啦观看免费观看视频高清| 97超视频在线观看视频| 午夜精品一区二区三区免费看| 99久久中文字幕三级久久日本| 一卡2卡三卡四卡精品乱码亚洲| 听说在线观看完整版免费高清| 亚洲人成网站高清观看| 又粗又硬又长又爽又黄的视频 | 日韩中字成人| 欧美另类亚洲清纯唯美| 久久精品国产亚洲av香蕉五月| 蜜桃久久精品国产亚洲av| 精品久久久久久久久亚洲| 国产69精品久久久久777片| 男人舔奶头视频| 最近的中文字幕免费完整| 天天躁夜夜躁狠狠久久av| 在线免费十八禁| 寂寞人妻少妇视频99o| 给我免费播放毛片高清在线观看| 国产免费一级a男人的天堂| 久久久色成人| 男的添女的下面高潮视频| 亚洲国产精品成人久久小说 | 如何舔出高潮| 最新中文字幕久久久久| 少妇熟女aⅴ在线视频| 偷拍熟女少妇极品色| 91狼人影院| 国产精品女同一区二区软件| 亚洲最大成人中文| 六月丁香七月| 九九在线视频观看精品| 给我免费播放毛片高清在线观看| 国产高潮美女av| 成人av在线播放网站| 边亲边吃奶的免费视频| 非洲黑人性xxxx精品又粗又长| 男人舔奶头视频| 免费在线观看成人毛片| 亚洲欧美精品专区久久| 国产精品久久久久久久久免| 男女做爰动态图高潮gif福利片| 国产精品久久电影中文字幕| av免费在线看不卡| 免费观看精品视频网站| 91久久精品国产一区二区三区| 18+在线观看网站| 欧美性猛交╳xxx乱大交人| 国产伦一二天堂av在线观看| 人人妻人人澡欧美一区二区| 国产单亲对白刺激| 国内揄拍国产精品人妻在线| 老师上课跳d突然被开到最大视频| 91麻豆精品激情在线观看国产| 中文字幕人妻熟人妻熟丝袜美| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产精品无大码| 干丝袜人妻中文字幕| 美女大奶头视频| 成人亚洲精品av一区二区| 搡女人真爽免费视频火全软件| 久久精品国产鲁丝片午夜精品| 国产精品一及| 欧美潮喷喷水| 听说在线观看完整版免费高清| 国产伦一二天堂av在线观看| 91久久精品国产一区二区成人| 波野结衣二区三区在线| 99久久精品国产国产毛片| 国产免费一级a男人的天堂| 国产精品美女特级片免费视频播放器| 国产亚洲5aaaaa淫片| 男人狂女人下面高潮的视频| 男女做爰动态图高潮gif福利片| 久久精品国产亚洲av天美| 亚洲av男天堂| 成人欧美大片| 一级毛片久久久久久久久女| 三级经典国产精品| 午夜精品在线福利| 日日啪夜夜撸| 欧美成人免费av一区二区三区| 久久亚洲精品不卡| 成人美女网站在线观看视频| 成人特级av手机在线观看| 国产午夜精品论理片| 日日干狠狠操夜夜爽| 亚洲色图av天堂| 亚洲精品国产av成人精品| 两个人的视频大全免费| 久久99热6这里只有精品| av国产免费在线观看| 免费黄网站久久成人精品| 美女国产视频在线观看| 99热精品在线国产| 久久综合国产亚洲精品| 国产精品一区二区在线观看99 | 长腿黑丝高跟| 色5月婷婷丁香| 长腿黑丝高跟| 又爽又黄a免费视频| 精品久久久久久久久亚洲| 亚洲国产精品成人久久小说 | 91久久精品国产一区二区成人| 日日摸夜夜添夜夜添av毛片| 久久6这里有精品| 在线a可以看的网站| 国产黄色视频一区二区在线观看 | 嫩草影院精品99| 国产成人freesex在线| 一个人观看的视频www高清免费观看| 夫妻性生交免费视频一级片| 在线免费十八禁| 欧美日本亚洲视频在线播放| 亚洲美女搞黄在线观看| 精品午夜福利在线看| 少妇的逼水好多| 日本爱情动作片www.在线观看| 国产精品1区2区在线观看.| 国产精品嫩草影院av在线观看| 天天躁夜夜躁狠狠久久av| 干丝袜人妻中文字幕| 草草在线视频免费看| 亚洲在线观看片| 亚洲最大成人av| 国产精品一区二区三区四区免费观看| 亚洲va在线va天堂va国产| 国产高潮美女av| 久久精品久久久久久久性| 久久这里只有精品中国| 一级毛片aaaaaa免费看小| 免费看美女性在线毛片视频| 国产精品一区二区在线观看99 | 一级毛片电影观看 | 日本-黄色视频高清免费观看| 国产成人aa在线观看| 美女大奶头视频| 天堂影院成人在线观看| 成人特级av手机在线观看| 看免费成人av毛片| 久久午夜亚洲精品久久| 国产亚洲5aaaaa淫片| 亚洲精品乱码久久久v下载方式| 97超视频在线观看视频| 中文资源天堂在线| 成人无遮挡网站| 一区二区三区高清视频在线| 99久国产av精品| 99久久中文字幕三级久久日本| 午夜视频国产福利| 久久韩国三级中文字幕| 91久久精品电影网| 一级av片app| 丰满人妻一区二区三区视频av| 亚洲人成网站在线观看播放| 久久韩国三级中文字幕| 男的添女的下面高潮视频| 91精品国产九色| 97超视频在线观看视频| 国产午夜精品论理片| 能在线免费看毛片的网站| 亚洲激情五月婷婷啪啪| 亚洲不卡免费看| av免费在线看不卡| 亚洲国产日韩欧美精品在线观看| 亚洲婷婷狠狠爱综合网| 亚洲精品乱码久久久久久按摩| 国产91av在线免费观看| 亚洲av成人av| 国产精品乱码一区二三区的特点| 久久久色成人| av国产免费在线观看| 天堂网av新在线| 午夜亚洲福利在线播放| 国产一级毛片七仙女欲春2| 亚洲成a人片在线一区二区| www.色视频.com| 天堂av国产一区二区熟女人妻| 国产黄片视频在线免费观看| 寂寞人妻少妇视频99o| 又爽又黄无遮挡网站| 免费不卡的大黄色大毛片视频在线观看 | 久久久国产成人免费| 91av网一区二区| 天堂网av新在线| 嫩草影院新地址| 国产午夜福利久久久久久| 成人鲁丝片一二三区免费| 在线免费观看的www视频| 国产亚洲精品av在线| 日韩高清综合在线| 欧美日本亚洲视频在线播放| 久久久a久久爽久久v久久| 美女xxoo啪啪120秒动态图| 在线天堂最新版资源| 日本黄色片子视频| 国产一区亚洲一区在线观看| 免费看光身美女| 1024手机看黄色片| 少妇被粗大猛烈的视频| 免费一级毛片在线播放高清视频| 男女视频在线观看网站免费| 国产三级中文精品| 直男gayav资源| 久久久精品欧美日韩精品| 九九在线视频观看精品| 啦啦啦啦在线视频资源| 亚洲精品色激情综合| 国产成人freesex在线| 日韩亚洲欧美综合| 在线观看一区二区三区| 波多野结衣巨乳人妻| videossex国产| 人人妻人人澡欧美一区二区| 99热这里只有是精品50| 日日干狠狠操夜夜爽| 欧美精品一区二区大全| 国产久久久一区二区三区| 日韩强制内射视频| 国产精品久久久久久久电影| 国产精品伦人一区二区| 国产黄a三级三级三级人| 国产黄片美女视频| 麻豆精品久久久久久蜜桃| 免费不卡的大黄色大毛片视频在线观看 | 夜夜夜夜夜久久久久| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 长腿黑丝高跟| 白带黄色成豆腐渣| 中文字幕精品亚洲无线码一区| 精品一区二区免费观看| 国产伦一二天堂av在线观看| 精品欧美国产一区二区三| 国产不卡一卡二| 中文字幕av在线有码专区| 一边摸一边抽搐一进一小说| 欧美日韩在线观看h| 亚洲av成人精品一区久久| 成人午夜高清在线视频| 国产三级中文精品| av在线老鸭窝| 最近中文字幕高清免费大全6| 国内少妇人妻偷人精品xxx网站| 亚洲第一电影网av| 在线观看美女被高潮喷水网站| 精品一区二区三区人妻视频| 欧洲精品卡2卡3卡4卡5卡区| 成人漫画全彩无遮挡| 国产精品一区二区性色av| 高清毛片免费看| 一级二级三级毛片免费看| 麻豆久久精品国产亚洲av| 国产久久久一区二区三区| 一级二级三级毛片免费看| 麻豆久久精品国产亚洲av| 别揉我奶头 嗯啊视频| 乱人视频在线观看| 我的老师免费观看完整版| 高清毛片免费看| 亚洲av男天堂| 亚洲欧美日韩东京热| 最近的中文字幕免费完整| 在线免费观看的www视频| www.av在线官网国产| 非洲黑人性xxxx精品又粗又长| 亚洲av免费在线观看| 大香蕉久久网| 亚洲成a人片在线一区二区| 国产国拍精品亚洲av在线观看| 欧美最新免费一区二区三区| 婷婷色av中文字幕| 欧美一区二区国产精品久久精品| 夜夜看夜夜爽夜夜摸| 成年av动漫网址| 免费看美女性在线毛片视频| 久久久久久九九精品二区国产|