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

    An Extraction Process for Optimal Utilization of Naphtha Based on Molecule Management

    2016-03-22 08:05:02
    中國(guó)煉油與石油化工 2016年3期

    (State Key Laboratory of Chemical Engineering, East China University of Science and Technology, Shanghai 200237)

    An Extraction Process for Optimal Utilization of Naphtha Based on Molecule Management

    Wang Tianxiao; Shen Benxian; Sun Hui

    (State Key Laboratory of Chemical Engineering, East China University of Science and Technology, Shanghai 200237)

    In this paper, the separation of aromatics from light naphtha by using extraction process was investigated for improving the utilization effciency of naphtha. It is indicated that, using a mixture of propylene carbonate-diethylene glycol as the solvent, the optimal extraction conditions cover: a volume fraction of propylene carbonate in the mixed solvent of 0.3, a solvent to feed ratio of 8, and an extraction temperature of 308 K. Through the extraction process, the aromatics mass fraction increases from 10.05% in naphtha to 27.74% in extract oil. It is found that the aromatics yield of extract oil,RA, reaches 92.11%. As a result, in comparison with naphtha, the potential aromatics content of extract oil increases impressively by 18.03%. Meanwhile, the aromatics content of raffnate oil decreases to 1.33%, and the normal paraffn yield of raffnate oil,Rp, is 76.61%. Accordingly, higher total olefns yields can be obtained when using raffnate oil as the raw material for steam cracking. The present results show that the utilization effciency of naphtha is improved through extraction process.

    optimal utilization; naphtha; extraction; aromatics separation

    1 Introduction

    Naphtha is the main feed for the steam cracking process and the catalytic reforming process to produce basic chemicals for the industry. As the demand for naphtha in both processes increases, the contradiction emerges. Hence, the optimal utilization of naphtha attracts extensive research interests in the petrochemical industry.

    The ineffcient utilization of naphtha is being referred to frequently because naphthenes which comprise the potential aromatic content in naphtha are preferred in the catalytic reforming process, while aromatic hydrocarbons can lead to coking during the steam cracking process[1-2]. Therefore, separating aromatic components from naphtha, which is used as a feed for catalytic reforming process, is necessary for the optimal utilization of naphtha fraction. However, aromatics can hardly be separated from naphtha by traditional distillation because of the presence of azeotropes and components with close boiling points[3-4]. The molecule management strategy, which aims at separating the components from a mixture according to different usage requirements of products in molecule, exhibits the superiority in improving the utilization effciency of naphtha[2,5]. As normal paraffins are the preferred components in steam cracking process, the adsorption technology using 5 A zeolite for separating normal paraffins from naphtha is an effcient way for optimal utilization of naphtha[1-2,6-8]. However, the presence of olefns in naphtha leads to coke deposition on 5 A zeolite, which could largely reduce the performance of separation process[9]. Extraction could be considered to be a simple and practical way for separation of naphtha containing olefns. Upon applying an extraction process by using a solvent with high capacity in the separation of aromatics from naphtha, the extract oil with high potential aromatic content can be separated from naphtha as a good feed for catalytic reforming process, meanwhile, the remaining raffnate oil with low aromatics content and relatively highn-paraffin content is a good feed for steam cracking process.

    The aromatics extraction process for naphtha has been widely studied[4,10-12], followed by distillation for solvent recovery. A number of polar solvents, such as sulfolane (SUL)[13], diethylene glycol (DEG)[14], triethylene glycol[15], tetraethylene glycol[16-17], propylene carbonate(PC)[18-20], ionic liquids[21], and other solvents[22], have been studied for separating aromatics from naphtha. The conventional extraction processes, such as the liquid-liquid extraction, the extractive distillation and the azeotropic distillation, have been used for the separation of aromatic hydrocarbons from the mixture with aromatics content exceeding 20%[12]. However, it is a great challenge to separate aromatics from the naphtha which contains less than 20% of aromatics by using the extraction process.

    The present work is aimed at the optimal utilization of light naphtha via liquid-liquid extraction process. Assisted by the Hansen solubility parameter, propylene carbonate (PC) was selected as a suitable solvent for extraction. Then the extraction conditions for PC were optimized. Thereafter, diethylene glycol (DEG) was selected as the cosolvent, and the optimal conditions for the mixed solvent (PC+DEG) were studied through the orthogonal test.

    2 Experimental

    2.1 Materials

    The light naphtha was obtained from the Shandong Jingbo Petrochemical Co., Ltd., with itsn-paraffins mainly consisting ofn-propane ton-decane and its boiling points ranging from 313 K to 433 K. The contents ofn-paraffins, iso-paraffins, olefins, cycloparaffins and aromatics analyzed by gas-chromatograph were 25.90%, 33.73%, 0.54%, 29.78% and 10.05%, respectively, which could change slightly due to the environment and storage conditions.N,N-dimethyl formamide (DMF),N-methylpyrrolidone (NMP), diethylene glycol (DEG), triethylene glycol (TEG) and dimethyl sulfoxide (DMSO) were supplied by the Lingfeng Chemical Reagent Co., Ltd. in Shanghai, with the quality of all materials meeting the chemical reagent purity. Sulfolane (SUL) and propylene carbonate (PC) were supplied by the Sinopharm Chemical Reagent Co., Ltd. in Shanghai, with the quality of all materials meeting the chemical reagent purity. Di-(2-ethylhexyl) sebacate (dioctyl sebacate, DOC) was supplied by the Huakai Resin Co., Ltd. in Jining, with its quality reaching the technically pure grade. Naphtha and solvents were used without further purifcation.

    2.2 Extraction experiments

    The flow diagram of extraction process is shown in Figure 1 (b). The extraction process includes an extraction column, a washing tower and a distillation column. The washing tower and distillation column were used to remove the solvent from the raffinate phase and the extract phase, respectively.

    Figure 1 Flow diagrams of three-stage counter-current extraction experiment (a) and the extraction process (b)

    The flow diagram of the three-stage counter-current extraction experiment[23]is shown in Figure 1 (a), which is used to simulate a three-stage counter-current extraction column. The single-stage extraction experiment and the three-stage counter-current extraction experiment were conducted by using liquid separating funnels and thermostatic bath. Each extraction was followed by shaking for at least 5 min and then settling for at least 10 min to form two clear phases under specifed conditions. The recovery of solvent was carried out for handling both the extract phase and the raffnate phase, respectively, to obtain the extract oil (E) and the raffnate oil (R), respectively. Due to the strong polarity which solvents possess, a small quantity of solvent existing in the raffnate phase can be removed via washing with water. Washing was conducted at a water/oil ratio of 1.0 in order to obtain raffnate oil from the raffnate phase. In consideration of the great difference in boiling point between oil and solvent, with the exception of DMF, the extract oil could be separated from the extract phase through distillation. DMF was removed from the extract phase by washing at a water/oil ratio of 1.0 owing to its similar boiling point to naphtha.

    The water to oil ratio and solvent to feed ratio (S/F) are expressed in volume ratio. The PC content of mixed solvent is expressed in volume fraction.

    2.3 Composition analysis

    The composition of extract oil and raffnate oil was analyzed by a GC-920 gas chromatograph (Shanghai Haixin Chromatograph Instrument Co., Ltd.) equipped with a flame ionization detector (FID) and a PONA column (0.3 mm i.d.×20 m). The hydrocarbon group compositions of naphtha and products were determined by employing the RIPP-Gasoline Composition Analyzer Software V2.01 developed by the Research Institute of Petroleum Processing.

    2.4 Hansen solubility parameter (HSP)

    HSP was used as a convenient way to predict the miscibility of solvents by solubility parameter “distance”(Ra). Hansen[24-27]divides the solubility parameter[28-29]into several individual parts. The relationship between total solubility parameter and HSP is:

    where the subscript d represents dispersion, p represents polar, and hb represents hydrogen bonding.

    Based on HSP, the solubility parameter “distance”between solvent and solute,Ra, is calculated by equation (2)[26]. To use HSP, a detailed database can be found in the literature[24,30-31].

    2.5 Orthogonal analysis

    The optimal solvent composition and extraction conditions for mixed solvent were determined by orthogonal test. Three indexes were selected for analysis of experimental results. Both direct analysis and variance analysis were used to evaluate the results. To make general score for each scheme of conditions, weights for evaluation indexes were computed subjectively by the Analytic Hierarchy Process[32](AHP) and objectively by the entropy method[33]through MATLAB. To make the full use of information, the optimal weight was calculated by considering both subjective weight and objective weight based on the sum of square of deviations[34].

    Variance analysis was carried out for evaluating the indexes to distinguish whether the results of orthogonal test were caused by factors or error. The sum of square of deviations,SS, for factor A is calculated by equation (3).

    whereKis the sum of indexes,n=16 is the experimental time.

    2.6 Product evaluation

    To investigate the performance of extraction process for optimal utilization of naphtha, the performance of naphtha and extract oil for catalytic reforming process were evaluated by the potential aromatics content,P, which is defned as the aromatics content of product when all of the cycloparaffin components are converted into aromatics. In this study, the potential aromatics contents of C6, C7and C8hydrocarbons were calculated for performance evaluation.

    whereMiis the molecular weight of cycloparaffn componenti,Niis the mass fraction of cycloparaffn componenti, andAiis the mass fraction of aromatic componenti.

    3 Results and Discussion

    3.1 Solvent selection

    The selection of solvent was performed by means of the Hansen solubility parameter. Suitable solvents were selected atT=313 K andS/F=3 through single-stage extraction.

    Asn-paraffins are the preferred components for steam cracking process, they possess low reactivity in catalytic reforming process. To preventn-paraffns from being extracted into the extract oil, the HSP ofn-paraffns is calculated to fnd the solvents having overwhelming aromatics solubility overn-paraffns solubility. In case of more than two components in a mixture, the HSP for the mixture can be calculated by equation (5)[24]. The HSP for the mixture ofn-paraffns is listed in Table 1.

    whereφiis the volume fraction of componentiin the mixture.

    Table 1 HSP of n-paraffins[24]

    To determine the appropriate solubility parameter “distance”Ra, DMF, NMP, SUL, DEG and DOC were used in the preliminary study. The result shows that NMP and DOC are good solvents for naphtha, which can even be miscible with naphtha atS/F=1; DMF becomes miscible with naphtha at the stated conditions, and can extract about 30% of naphtha atS/F=1; and DEG and SUL have certain solubility in naphtha.

    A suitable solvent for extractive separation of naphtha needs a considerable solubility for aromatics and a poor solubility forn-paraffins, and a certain solubility (not too much) for naphtha to balance the products. In this context, we have to find a solvent with HSP similar to the aromatic ones (with a small amount of polar force and hydrogen-bonding force), but unlike then-paraffnic ones (without polar force or hydrogen-bonding force). As regards the naphtha in this study, DMF is considered as a good solvent at the boundary surface of HSP solubility sphere with its radius ton-paraffnRaDMF=18.3 MPa1/2calculated by equation (2). The less theRais, the closer the position of solvent to the position of naphtha would be. In other words, solvents inside the sphere are considered as good solvents, and the bad ones are outside or close to it. To exclude good solvents for naphtha, theRaof appropriate solvents is determined to be larger thanRaDMF.

    In order to display the solubility of solvent clearly, the Hansen three-dimensional solubility parameter sphere for the naphtha-solvent system is established for 16 kinds of conventional solvents. The coordinates of the solubility volume center are the HSP ofn-paraffin (δd=15.1,δp=0.0,δhb=0.0) withRa=18.3 MPa1/2. As shown in Figure 2, NMP, DOC and DMF are inside the sphere indicating their high solubility to naphtha, while SUL, DEG are outside the sphere, which agrees well with our preliminary results.

    Figure 2 Hansen three-dimensional solubility sphere for then-paraf fi n-solvent system

    To prevent the solvent from getting miscible with naphtha, DMF, DEG, SUL, PC, TEG and DMSO were obtained as solvents for extraction. Then the single-stage extraction experiments were performed by using 6 solvents mentioned above.

    The yield of aromatics in the extract oil,RA, is defned as the percentage of aromatic components extracted from naphtha into extract oil.

    At the same time, then-paraffn yield of raffnate oil,Rp, is defned by equation (7). The aromatics recovery refers toRA,and then-paraffn recovery refers toRp.

    The yields of aromatics and extract oil are shown in Figure 3. The results indicate that PC, DMF and DMSO have better capacity for extractive separation of naphtha, which means a lowerS/Fratio for aromatics recovery. The yield of aromatics decreases in the following order: DMF>DMSO>PC>SUL>DEG>TEG.

    Figure 3 Yields of aromatics and extract oil. Extraction runs were operated atT=313 K andS/F=3 for DEG, TEG, PC, SUL, DMSO and atT=313 K andS/F=1 for DMF

    Asn-paraffin components are good feedstocks for steam cracking process, solvents should have the selectivity for the separation of aromatics andn-paraffns. The selectivity of the solvent,S, is defned by equation (8). The selectivity of solvents is shown in Figure 4, which decreases in the following order: SUL>DMSO>TEG>PC>DEG>DMF.

    By taking solubility and selectivity into consideration, it is indicated that DMF, DMSO, PC and SUL show advantages in extractive separation of naphtha. DMF, PC and DMSO are found to have high solubility, while DMF possesses low selectivity. SUL and DMSO have high selectivity, but low thermal stability. As a result, PC is selected as a suitable solvent for extractive separation of light naphtha.

    3.2 Extraction conditions optimization for PC

    3.2.1 Effect of temperature

    Once the solvent was selected, the effect of temperature (T) and solvent to feed ratio (S/F) on extraction were studied through the three-stage counter-current extraction experiment. Via the three-stage counter-current extraction, the performance of PC in the multi-stage countercurrent extraction was also tested.

    Figure 4 Selectivity of solvents. Extraction runs took place atT=313 K andS/F=3 for DEG, TEG, PC, SUL, DMSO and atT=313 K andS/F=1 for DMF

    Firstly, the effect of temperature on the composition of products was studied at a temperature of 303, 308, 313, 318 and 323 K, respectively, and aS/Fratio of 6. In Figure 5, it can be seen that the aromatics content of raffinate oil is reduced to about 1.0%, which indicates that PC has a satisfactory capacity for separating aromatics from naphtha in the multi-stage counter-current extraction. With the temperature increasing from 303 K to 323 K, the aromatics content of extract oil decreases significantly because of the increasing solubility for other components; then-paraffns content of raffnate oil increases and reaches a plateau at 318 K, then decreases due to the increasing solubility ofn-paraffns in solvent at a higher temperature. Upon considering the composition of extraction products, the optimal extraction temperature is 313 K.

    3.2.2 Effect ofS/Fratio

    The effect of solvent to feed ratio (S/F)on the composition of extraction products was studied at aS/Fratio of 2, 4, 6, 8, and 10, respectively, and an extraction temperature of 313 K. The results are shown in Figure 6. The aromatics content of raffnate oil is reduced to about 1.0% when theS/Fratio is higher than 6. With the increase ofS/Fratio, the aromatics content of extract oil decreases signifcantly while then-paraffins content of extract oil increases. Meanwhile, then-paraffins content of raffinate oil increases, and iso-paraffns content of raffnate oil increases when theS/Fratio increases from 2 to 6, then decreases with a further increase in the S/F ratio. To separate aro-matics fromn-paraffns, the optimal solvent to feed ratio is 6.

    Figure 5 Effect of temperature on the composition of raf fi nate oil (a) and extract oil (b) at aS/Fratio of 6

    By using an optimal conditions ofS/F=6 andT=313 K throughout the three-stage counter-current extraction process, pure PC can increase the aromatics content from 10.05% in naphtha to 17.58% in extract oil and reduce the aromatics content of raffnate oil to 0.86%. The yield of extract oil is 45.39%. Hence, PC shows high aromatics solubility under optimal conditions. However, then-paraffns content of extract oil reaches 21.38%, which indicates that the solubility of PC ton-paraffins is also relatively high.

    3.3 Analysis of orthogonal test

    Sincen-paraffn components can lead to high total olefns yield in steam cracking process but low aromatics yield in catalytic reforming process, they should remain in raffnate oil based on the idea of molecule management. To improve the recovery ofn-paraffins in raffinate oil, diethylene glycol (DEG) has been selected as a cosolvent, which has lower solubility ton-paraffins and a close boiling point (518 K) to PC (515 K). Orthogonal test was designed to determine the optimal conditions suitable to the mixed solvent.

    Figure 6 Effect ofS/Fratio on the composition of raf fi nate oil (a) and extract oil (b) atT=313 K

    The major factors considered in the orthogonal test covered: the volume fraction of PC in the mixed solvent (x), the solvent to feed volume ratio (S/F), and the extraction temperature (T). Other conditions, such as the pressure, were limited in orthogonal test to reduce error. The [L16(45)] orthogonal table was chosen in this study. As shown in Table 2, each factor has four levels in equal interval, respectively. The aromatics yield of extract oil (RA) and the distribution coefficient ofn-paraffins (Kp) were used to evaluate the extraction performance. Meanwhile, the solubility (defined as the yield of extract oil, %) was also considered in order to obtain the comparative yields of extract oil and raffnate oil.

    According to the properties of indexes, the evaluation matrix was normalized by equation (10).RAandKpare beneft types, which means that they are positively related to the evaluation. And the solubility is a quantitative type with the desired value of 0.5.

    To calculate the subjective weightwAof evaluation indexes, the judgment matrixAof AHP is defined according to the importance of indexes. To improve the recovery ofn-paraffins, the distribution coefficient ofn-paraffins (Kp) is considered as the most important evaluation index for orthogonal test. Then the aromatics yield of extract oil (RA) is in the second place. To keep the balance of products, solubility is also considered. The consistency ratio (CR) of matrixAis 0.015 8. Hence, the AHP weightwAcalculated by matrixAis available.

    By using MATLAB, the subjective weightwAand the objective weightwSof evaluation indexes were calculated. As a result, the optimized combination weight vectorWcwas obtained to calculate the general scores of orthogonal test.

    3.3.1 Direct analysis

    By using the weight vectorWcobtained above, scores were calculated for the normalized matrixZ. The result of direct analysis is shown in Table 2. The range of scores for each factor,R=(knm)max-(knm)min, was calculated to investigate the impact of factors on the results. Larger impact of the factor on scores is made with largerR.

    Table 2 Results of orthogonal test

    Through the direct analysis, the effect of factors on score follows the order: C>B>A. The result indicates that the best scheme of conditions for orthogonal test is A1B3C2, which impliesx=0.2,S/F=8, andT=308 K.

    3.3.2 Variance analysis

    The result of variance analysis is shown in Table 3. The degree of freedomfandFvalue were calculated by the statistical method. By comparingFvalue with the critical value, the factor withFvalue being larger than critical value (α=0.01) is marked as (* *); the factor withFvalue being larger than critical value (α=0.05) is marked as (*).

    If the factor’s sum of squares of deviations is smaller than the errorQ, then the factor will be considered as erroneous, in other words, the new errorQ*, the degree of freedomf, and theFvalues of factors should be recalculated. The 4thand 5throws of [L16(45)] in the orthogonal table were used to calculate the original errorQ. In Table 3, the results show that factors A and B have great significance toRAwithFvalues being larger than the critical value (α=0.01), while the factor C has no significance to the index. Factors A and B have a 99% signifcance toKp, and the factor C has a 95% signifcance to it.Kpis the most important index, which weighs 52.07% in the score. Hence, these factors are all major factors of extraction and evaluation. The variance analysis result of solubility is similar to theRAone. Factors A and B have a 99% significance to the solubility. However, the temperature shows no significance to the solubility. This result might be ascribed to the limited variation of temperature investigated in this orthogonal test.

    Table 3 Results of variance analysis on RA, Kp, and solubility

    It can be concluded that all the factors have great signifcance to the extraction operation, therefore, the test error can be ignored, and the result of orthogonal test is reliable.

    3.3.3 Three-stage counter-current extraction for mixed solvent

    To investigate the performance of multi-stage extraction by using the mixed solvent, the three-stage countercurrent extraction experiment was operated at the best scheme of conditions for orthogonal test. The result is shown in Table 4.

    Table 4 Result of three-stage counter-current extraction by using mixed solvent

    The three-stage counter-current extraction by using mixed solvent achieves a relatively highKp, andRpreaches 82.48%, whileRAis below 90%. To realize the separation of aromatics from naphtha, the satisfedRAshould be over 90%. As a result of variance analysis,RAincreases significantly with the increase ofxandS/Fratio. The load of extraction process will increase if theS/Fratio increases. Hence, it is an efficient way to increase the volume fraction of PC in the mixed solvent to improveRA. The influence ofxon extraction is shown in Figure 7. It is indicated that the solubility of mixed solvent to aromatics increases with the increase ofx, while the solubility ton-paraffins also increases in some degree, which leads to the decrease ofKpandRp. It can be noticed that the general score andKpdecrease signifcantly whenxincreases from 0.2 to 0.4. To ensure the recovery ofn-paraffns in the raffnate oil,xshould not be larger than 0.4. Therefore, further experiment was carried out under conditions leading to increasingxslightly from 0.2 to 0.3, viz.:x=0.3,S/F=8, andT=308 K, which is between the two investigated schemes of conditions in orthogonal test. The result is shown in Table 4. With the increase ofxfrom 0.2 to 0.3,RAincreases to a satisfactory value of 92.11% as expected, andRpdecreases slightly from 82.48% to 76.61%. The yield of extract oil is 33.37%. To achieve high aromatics yield of extract oil and relatively highn-paraffins yield of raffnate oil, the optimal conditions suitable to the mixed solvent throughout the three-stage countercurrent extraction are:x=0.3,S/F=8, andT=308 K.

    Figure 7 The in fl uence of factorxon extraction

    Table 5 Composition of naphtha and extraction products %

    The composition of naphtha and extraction products obtained under optimal conditions is shown in Table 5. It is indicated that the aromatic content of extract oil increases to 27.74%. Meanwhile, the aromatics content of raffnate oil decreases significantly to 1.33%, and then-paraffin content of raffinate oil increases to 29.78%. It is found that the extraction process by using the mixed solvent can extract most of aromatics into extract oil, and an acceptablen-paraffn yield of raffnate oil is also obtained.

    The performance of naphtha and extraction products was evaluated according to the demand of processes. For catalytic reforming process, the potential aromatics content of raw material can be calculated by equation (4); for steam cracking process, the total olefins (ethylene, propylene and butadiene) yields (yt) of raw material can be estimated by equation (13)[35], which is an empirical equation regressed by raw materials with similar composition. As compared with naphtha, the potential aromatic content of extract oil increases by 18.03%. Meanwhile, the total olefins yields of raffinate oil may increase by 2% due to the increase inn-paraffn content. As a result, the application of PC-DEG mixed solvent in the extraction process boosts a greater capacity in improving the utilization effciency of naphtha.

    whereXnis the correctedn-paraffins content excluding the aromatics content (except for ethyl benzene).

    4 Conclusions

    Based on the molecule management, the extract oil with high potential aromatic content has been separated from light naphtha by using extraction process. A mixture of propylene carbonate and diethylene glycol has been found as a suitable solvent. In addition, the optimal extraction conditions include: a volume fraction of PC in mixed solvent of 0.3, a solvent to feed ratio of 8, and an extraction temperature of 308 K. By employing the extraction process, the aromatics content increases from 10.05% in naphtha to 27.74% in the extract oil. Furthermore, the aromatics yield of extract oil reaches 92.11%. As a result, the potential aromatics content increases impressively from 30.13% in naphtha to 48.16% in the extract oil. Meanwhile, the aromatics content of raffinate oil decreases to 1.33% while then-paraffin content increases from 25.90% in naphtha to 29.78% in the raffnate oil. Then-paraffn yield of raffnate oil can reach 76.61%. Consequently, as compared with naphtha, a higher total olefins yield can be obtained when using raffinate oil as the raw material for steam cracking. Hence, the current extraction process achieves an effcient separation of aromatics from naphtha, and the utilizationeffciency of naphtha is improved largely.

    Acknowledgements: This work was financially supported by the Natural Science Foundation of Shanghai, China (16ZR1408100) and the Fundamental Research Funds for the Central Universities of China (22A201514010).

    [1] Shen B, Liu J, Chen H, et al. Method for optimizing naphtha utilization by molecular sieve adsorption: China, CN1710030A[P]

    [2] Liu J, Shen B, Sun H. Adsorption behaviour of normal paraffns in a fxed bed adsorber containing 5 ? molecular sieves[J]. Adsorption Science and Technology, 2006, 24(4): 311-320

    [3] Meindersma G W, Podt A, Klaren M B, et al. Separation of aromatic and aliphatic hydrocarbons with ionic liquids[J]. Chemical Engineering Communications, 2006, 193(11): 1384-1396

    [4] Rawat B S, Gulati I B. Liquid-liquid equilibrium studies for separation of aromatics[J]. Journal of Applied Chemistry and Biotechnology, 1976, 26(8): 425-35

    [5] Sun H, Shen B, Liu J.n-Paraffns adsorption with 5A zeolites: The effect of binder on adsorption equilibria[J]. Separation and Purification Technology, 2008, 64(1): 135-139

    [6] Miano F. Adsorption of hydrocarbon vapour mixtures onto zeolite 5A[J]. Colloids and Surfaces A—Physicochemical and Engineering Aspects, 1996, 110(1): 95-104

    [7] Barcia P S, Silva J A C, Rodrigues A E. Octane upgrading of C5/C6light naphtha by layered pressure swing adsorption[J]. Energy Fuels, 2010, 24(9): 5116-5130

    [8] Cao J, Shen B X, Liu J C. Optimal operation of simulated moving bed technology on utilization of naphtha resource[J]. Separation Science and Technology, 2012, 48(2): 246-253

    [9] Silva J A C, Mata V G, Dias M M, et al. Effect of coke in the equilibrium and kinetics of sorption on 5 A molecular sieve zeolites[J]. Industrial and Engineering Chemistry Research Fundamentals, 2000, 39(4): 1030-1034

    [10] Hamid S H, Ali M A. Comparative study of solvents for the extraction of aromatics from naphtha[J]. Energy Sources, Part A: Recvoery, Utilization and Environmental Effects, 1996, 18(1): 65-84

    [11] Weissermel K, Arpe H-J. Aromatics-Production and Conversion[M]//Industrial Organic Chemistry, Forth Edition. Weinheim: Wiley-VCH Verlag GmbH, 2008

    [12] Meindersma G W, Hansmeier A R, De Haan A B. Ionic liquids for aromatics extraction. Present status and future outlook[J]. Industrial and Engineering Chemistry Research, 2010, 49(16): 7530-7540

    [13] Gaile A A, Zalishchevskii G D, Gafur N N, et al. Removal of aromatic hydrocarbons from reforming naphtha. Combined extraction-extractive-azeotropic distillation process[J]. Chemistry and Technology of Fuels and Oils, 2004, 40(4): 215-221

    [14] Rodrigues M F M, Pinheiro R S, De Sant’ana H B, et al. Removal of aromatic hydrocarbons from hydrocarbon mixture using glycols at 303.15 K and 333.15 K and atmospheric pressure: Experimental and calculated data by NRTL and UNIQUAC models[J]. Fluid Phase Equilibria, 2015, 387: 135-142

    [15] Banerjee T, Sahoo R K, Rath S S, et al. Multicomponent liquid-liquid equilibria prediction for aromatic extraction systems using COSMO-RS[J]. Industrial and Engineering Chemistry Research, 2007, 46(4): 1292-1304

    [16] Saha M, Rawat B S, Khanna M K, et al. Liquid-liquid equilibrium studies on toluene+heptane+solvent[J]. Journal of Chemical and Engineering Data, 1998, 43(3): 422-426

    [17] Wang W, Gou Z, Zhu S. Liquid-liquid equilibria for aromatics extraction systems with tetraethylene glycol[J]. Journal of Chemical and Engineering Data, 1998, 43(1): 81-83

    [18] Ali S H, Lababidi H M S, Merchant S Q, et al. Extraction of aromatics from naphtha reformate using propylene carbonate[J]. Fluid Phase Equilibria, 2003, 214(1): 25-38

    [19] Fahim M A, Merchant S Q. Liquid-liquid equilibria of systems containing propylene carbonate and some hydrocarbons[J]. Journal of Chemical and Engineering Data, 1998, 43(5): 884-888

    [20] Lababidi H M S, Ali S H, Fahim M A. Optimization of aromatics extraction of naphtha reformate by propylene carbonate/diethylene glycol mixed solvent[J]. Industrial and Engineering Chemistry Research, 2006, 45(14): 5086-5097.

    [21] Meindersma G W, De Haan A B. Conceptual process design for aromatic/aliphatic separation with ionic liquids[J]. Chemical Engineering Research and Design,2008, 86(7): 745-752.

    [22] Rodriguez N R, Requejo P F, Kroon M C. Aliphatic-aromatic separation using deep eutectic solvents as extracting agents[J]. Industrial and Engineering Chemistry Research, 2015, 54(45): 11404-11412.

    [23] Deng Xiu, Wu Junsheng. Chemical Separation Engineering[M]. Beijing: Science Press, 2013. (in Chinese)

    [24] Hansen C M. Hansen Solubility Parameters: A User’s Handbook[M]. Boca Raton, FL: CRC Press LLC, 2007.

    [25] Hansen C M. Three-dimensional solubility parameterkey to paint-component affnities: I. Solvents, plasticizers, polymers, and resins[J]. Journal of Paint Technology, 1967, 39(505): 104-117.

    [26] Hansen C M. Three-dimensional solubility parameter-key to paint component affnities. II. Dyes, emulsifers, mutual solubility and compatibility, and pigments[J]. Journal of Paint Technology, 1967, 39(511): 505-510.

    [27] Hansen C M, Skaarup K. Three-dimensional solubility parameter-key to paint component affnities. III. Independent calculation of the parameter components[J]. Journal of Paint Technology, 1967, 39(511): 511-514.

    [28] Hildebrand J H. The Solubility of Non-Electrolytes[M]. Reinhold, NY: Reinhold Pub. Corp., 1950.

    [29] Hildebrand J H, Scott R L. Regular Solutions[M]. Upper Saddle River, NJ: Prentice-Hall, 1962.

    [30] Barton A F M. CRC Handbook of Solubility Parameters and Other Cohesion Parameters[M]. Boca Raton, FL: CRC Press, Inc., 1983.

    [31] Heremisinoff N P. Industrial Solvents Handbook, Second Edition[M]. New York: Marcel Dekker, 2003.

    [32] Saaty T L. How to make a decision: The analytic hierarchy process[J]. European Journal of Operational Research, 1990, 48(1): 9-26

    [33] Shannon C E. A mathematical theory of communication[J]. The Bell System Technical Journal, 1948, 27(3): 379-423

    [34] Chen Wei, Xia Jianhua. An optimal weights combination method considering both subjective and objective weight information[J]. Mathematics in Practice and Theory, 2007, 37(1): 17-22 (in Chinese)

    [35] Wang D Q. Recombination for making good use of naphtha by solvent extraction technology[D]. East China University of Science and Technology, Shanghai, 2008 (in Chinese)

    Received date: 2016-07-07; Accepted date: 2016-07-24.

    Dr. Sun Hui, Telephone: +86-21-64252916; E-mail: sunhui@ecust.edu.cn.

    久久97久久精品| 精品国产超薄肉色丝袜足j| 亚洲精品第二区| 男女午夜视频在线观看| 99国产精品免费福利视频| 日韩av不卡免费在线播放| 久久精品久久久久久噜噜老黄| 黄片播放在线免费| 麻豆乱淫一区二区| 精品一品国产午夜福利视频| 精品国产露脸久久av麻豆| 日韩伦理黄色片| 国产人伦9x9x在线观看| 亚洲欧美激情在线| 欧美黄色片欧美黄色片| 制服诱惑二区| 啦啦啦中文免费视频观看日本| 麻豆乱淫一区二区| 在线精品无人区一区二区三| 成人黄色视频免费在线看| 亚洲欧洲精品一区二区精品久久久 | 成年女人毛片免费观看观看9 | 国产精品99久久99久久久不卡 | 我的亚洲天堂| 天堂俺去俺来也www色官网| 赤兔流量卡办理| 男女之事视频高清在线观看 | 美女扒开内裤让男人捅视频| 老司机在亚洲福利影院| 午夜福利一区二区在线看| 亚洲一码二码三码区别大吗| 欧美日韩亚洲国产一区二区在线观看 | 91老司机精品| 亚洲精品国产区一区二| www.av在线官网国产| 国产精品一国产av| 丁香六月欧美| 伊人久久国产一区二区| 99精国产麻豆久久婷婷| 啦啦啦在线观看免费高清www| 午夜福利影视在线免费观看| 精品一区二区三卡| 大话2 男鬼变身卡| 成人影院久久| 亚洲国产精品成人久久小说| 赤兔流量卡办理| 日韩一本色道免费dvd| 免费久久久久久久精品成人欧美视频| 狂野欧美激情性bbbbbb| 考比视频在线观看| 国产一区二区三区av在线| 国产一区二区 视频在线| 一边摸一边做爽爽视频免费| 欧美黑人欧美精品刺激| 99热国产这里只有精品6| 丁香六月欧美| 无遮挡黄片免费观看| 久久精品人人爽人人爽视色| 男女边摸边吃奶| 永久免费av网站大全| 亚洲国产精品一区二区三区在线| 麻豆av在线久日| 欧美中文综合在线视频| 超碰97精品在线观看| 免费在线观看完整版高清| 欧美亚洲 丝袜 人妻 在线| 国产一区二区 视频在线| 天堂俺去俺来也www色官网| 深夜精品福利| 久久性视频一级片| 久久国产亚洲av麻豆专区| 久久久久久久国产电影| 国产熟女欧美一区二区| 亚洲免费av在线视频| 一级毛片我不卡| 亚洲精品一区蜜桃| 最近最新中文字幕大全免费视频 | 国产片特级美女逼逼视频| 男女无遮挡免费网站观看| 男女下面插进去视频免费观看| 精品午夜福利在线看| 国产精品熟女久久久久浪| 亚洲国产看品久久| 亚洲激情五月婷婷啪啪| 黄片播放在线免费| av在线app专区| 国产熟女欧美一区二区| 国产片内射在线| 日韩中文字幕欧美一区二区 | av免费观看日本| 日日爽夜夜爽网站| 亚洲国产毛片av蜜桃av| 亚洲精品自拍成人| 亚洲,一卡二卡三卡| 午夜精品国产一区二区电影| 肉色欧美久久久久久久蜜桃| 精品一区在线观看国产| 国产欧美日韩一区二区三区在线| 久久女婷五月综合色啪小说| 亚洲av福利一区| 在线 av 中文字幕| 久久亚洲国产成人精品v| 久久性视频一级片| 国产成人av激情在线播放| 女的被弄到高潮叫床怎么办| 免费观看人在逋| 啦啦啦在线观看免费高清www| 99久国产av精品国产电影| 美女午夜性视频免费| av网站免费在线观看视频| 国产视频首页在线观看| 国产深夜福利视频在线观看| 中文字幕人妻丝袜制服| 男女高潮啪啪啪动态图| 嫩草影视91久久| 肉色欧美久久久久久久蜜桃| 日本猛色少妇xxxxx猛交久久| 欧美成人午夜精品| 国产免费现黄频在线看| 大片电影免费在线观看免费| 国产无遮挡羞羞视频在线观看| 久久久久精品人妻al黑| 亚洲伊人久久精品综合| 亚洲欧美精品综合一区二区三区| 精品人妻熟女毛片av久久网站| 亚洲美女搞黄在线观看| 51午夜福利影视在线观看| 美女扒开内裤让男人捅视频| 中文字幕制服av| 欧美日韩亚洲综合一区二区三区_| 成年女人毛片免费观看观看9 | 国产精品一区二区精品视频观看| 99re6热这里在线精品视频| 免费在线观看黄色视频的| 女人爽到高潮嗷嗷叫在线视频| 在线观看一区二区三区激情| 免费人妻精品一区二区三区视频| 久久99热这里只频精品6学生| 免费黄网站久久成人精品| 国产有黄有色有爽视频| 夫妻午夜视频| 亚洲五月色婷婷综合| 亚洲第一青青草原| 国产成人免费观看mmmm| 久久精品久久久久久久性| 亚洲国产日韩一区二区| 亚洲伊人久久精品综合| 国产乱来视频区| 久久鲁丝午夜福利片| 少妇人妻久久综合中文| 侵犯人妻中文字幕一二三四区| 深夜精品福利| 少妇人妻久久综合中文| 女性被躁到高潮视频| 热99久久久久精品小说推荐| 国产精品秋霞免费鲁丝片| 在现免费观看毛片| 精品酒店卫生间| 国产精品亚洲av一区麻豆 | 69精品国产乱码久久久| 国产精品熟女久久久久浪| 亚洲欧美成人精品一区二区| 中国三级夫妇交换| 一边摸一边做爽爽视频免费| 欧美激情高清一区二区三区 | 久久久久久久久久久免费av| 久久久精品94久久精品| 亚洲国产精品国产精品| 十八禁人妻一区二区| 黄色毛片三级朝国网站| 久热这里只有精品99| 黄片无遮挡物在线观看| 久久精品亚洲熟妇少妇任你| 亚洲精品美女久久av网站| 18禁裸乳无遮挡动漫免费视频| 侵犯人妻中文字幕一二三四区| 啦啦啦 在线观看视频| 色吧在线观看| 婷婷成人精品国产| 男的添女的下面高潮视频| 久久99热这里只频精品6学生| 亚洲国产精品国产精品| 欧美激情 高清一区二区三区| 国产精品国产av在线观看| 欧美av亚洲av综合av国产av | 最近最新中文字幕免费大全7| 色综合欧美亚洲国产小说| 欧美精品av麻豆av| 中文字幕最新亚洲高清| 9色porny在线观看| 我的亚洲天堂| av在线播放精品| 国产1区2区3区精品| 卡戴珊不雅视频在线播放| 啦啦啦视频在线资源免费观看| 国产99久久九九免费精品| 久久韩国三级中文字幕| 一本久久精品| 9热在线视频观看99| 精品少妇一区二区三区视频日本电影 | 美国免费a级毛片| 久久久久网色| 亚洲三区欧美一区| 久久久久久免费高清国产稀缺| 国产精品蜜桃在线观看| 黑人巨大精品欧美一区二区蜜桃| 亚洲伊人色综图| 天天影视国产精品| 久久久精品免费免费高清| 欧美老熟妇乱子伦牲交| 岛国毛片在线播放| 精品国产一区二区三区久久久樱花| 免费女性裸体啪啪无遮挡网站| 在线天堂最新版资源| 男女下面插进去视频免费观看| 日韩大码丰满熟妇| 久久久久精品久久久久真实原创| 七月丁香在线播放| 国产亚洲午夜精品一区二区久久| √禁漫天堂资源中文www| 在线观看国产h片| 啦啦啦视频在线资源免费观看| 国产在视频线精品| 亚洲欧美精品综合一区二区三区| 999精品在线视频| 新久久久久国产一级毛片| 男人添女人高潮全过程视频| 久久av网站| 国产精品一区二区在线观看99| 肉色欧美久久久久久久蜜桃| 婷婷色麻豆天堂久久| 别揉我奶头~嗯~啊~动态视频 | 97精品久久久久久久久久精品| 国产一级毛片在线| 精品亚洲成a人片在线观看| 无遮挡黄片免费观看| 国产不卡av网站在线观看| 在线精品无人区一区二区三| 欧美日韩综合久久久久久| 交换朋友夫妻互换小说| av片东京热男人的天堂| 国产精品熟女久久久久浪| 伊人久久大香线蕉亚洲五| 日韩,欧美,国产一区二区三区| 国产女主播在线喷水免费视频网站| 亚洲精品国产一区二区精华液| av国产久精品久网站免费入址| 亚洲国产精品国产精品| 久久国产精品大桥未久av| 亚洲伊人久久精品综合| 嫩草影院入口| av在线老鸭窝| 在线看a的网站| 久久久久精品国产欧美久久久 | 精品人妻熟女毛片av久久网站| 成人毛片60女人毛片免费| 国产激情久久老熟女| 97精品久久久久久久久久精品| 不卡视频在线观看欧美| 日本一区二区免费在线视频| 免费日韩欧美在线观看| 午夜福利一区二区在线看| 国产极品粉嫩免费观看在线| 欧美日韩av久久| 日韩,欧美,国产一区二区三区| 一本—道久久a久久精品蜜桃钙片| 一级a爱视频在线免费观看| 久久久欧美国产精品| 亚洲欧美精品自产自拍| 成人漫画全彩无遮挡| 日韩制服丝袜自拍偷拍| 成年美女黄网站色视频大全免费| 青春草亚洲视频在线观看| 久久久久久久国产电影| 纵有疾风起免费观看全集完整版| 亚洲综合色网址| 国产精品免费视频内射| 国产成人午夜福利电影在线观看| 亚洲一码二码三码区别大吗| 桃花免费在线播放| 制服人妻中文乱码| 啦啦啦啦在线视频资源| 日韩,欧美,国产一区二区三区| 美女主播在线视频| 在线观看免费高清a一片| 国产在线一区二区三区精| 性高湖久久久久久久久免费观看| 亚洲欧洲国产日韩| 婷婷色综合www| 精品视频人人做人人爽| 老鸭窝网址在线观看| 久久久久人妻精品一区果冻| 欧美亚洲日本最大视频资源| 国产精品一区二区在线观看99| 国产人伦9x9x在线观看| 久久久精品免费免费高清| 国产精品人妻久久久影院| av线在线观看网站| 久久久久精品久久久久真实原创| 国产精品蜜桃在线观看| 国产熟女午夜一区二区三区| 免费看av在线观看网站| 欧美精品一区二区免费开放| 国产亚洲一区二区精品| 十八禁人妻一区二区| 老司机在亚洲福利影院| 最近最新中文字幕大全免费视频 | 亚洲欧美激情在线| 亚洲精品久久久久久婷婷小说| 久久 成人 亚洲| 少妇精品久久久久久久| 操美女的视频在线观看| 天天躁夜夜躁狠狠躁躁| 亚洲av欧美aⅴ国产| 久久97久久精品| 午夜福利,免费看| 亚洲精品国产色婷婷电影| 国产xxxxx性猛交| 一级毛片 在线播放| www.自偷自拍.com| 免费在线观看黄色视频的| 国产无遮挡羞羞视频在线观看| 亚洲精品日韩在线中文字幕| 最近2019中文字幕mv第一页| 亚洲精品视频女| 国产 精品1| 亚洲av日韩在线播放| 最近手机中文字幕大全| 欧美日韩亚洲综合一区二区三区_| 美女福利国产在线| 色94色欧美一区二区| 久久97久久精品| 亚洲成av片中文字幕在线观看| 丰满乱子伦码专区| 国产一卡二卡三卡精品 | av在线老鸭窝| 国产一区亚洲一区在线观看| 美女午夜性视频免费| 亚洲成国产人片在线观看| 国产深夜福利视频在线观看| 久久久精品国产亚洲av高清涩受| 国产欧美日韩一区二区三区在线| 免费在线观看黄色视频的| 丝袜在线中文字幕| 欧美日韩福利视频一区二区| 伊人久久大香线蕉亚洲五| 观看av在线不卡| 蜜桃国产av成人99| 99精品久久久久人妻精品| 精品国产超薄肉色丝袜足j| 视频在线观看一区二区三区| 在线免费观看不下载黄p国产| 午夜福利网站1000一区二区三区| 国产 精品1| 纵有疾风起免费观看全集完整版| 水蜜桃什么品种好| 婷婷色麻豆天堂久久| 久久久久国产精品人妻一区二区| 男女边摸边吃奶| 一二三四在线观看免费中文在| 国产福利在线免费观看视频| a级毛片黄视频| 久久99一区二区三区| 亚洲自偷自拍图片 自拍| 免费女性裸体啪啪无遮挡网站| 午夜福利影视在线免费观看| 国产成人精品久久久久久| 国产在视频线精品| 国产成人精品无人区| 亚洲四区av| 男女免费视频国产| 欧美日韩精品网址| 这个男人来自地球电影免费观看 | 日韩精品有码人妻一区| 嫩草影视91久久| 精品国产一区二区久久| 精品亚洲成国产av| 少妇的丰满在线观看| 丝袜喷水一区| 啦啦啦在线观看免费高清www| 国产又爽黄色视频| 美女扒开内裤让男人捅视频| 免费看不卡的av| 久久久久精品国产欧美久久久 | 妹子高潮喷水视频| 波多野结衣一区麻豆| 晚上一个人看的免费电影| 别揉我奶头~嗯~啊~动态视频 | 啦啦啦啦在线视频资源| 少妇被粗大猛烈的视频| 亚洲美女搞黄在线观看| 成人国产av品久久久| 国产精品成人在线| 亚洲精品中文字幕在线视频| 日韩伦理黄色片| 午夜精品国产一区二区电影| 99久久精品国产亚洲精品| 精品亚洲乱码少妇综合久久| 精品久久蜜臀av无| 国产av一区二区精品久久| 中文字幕最新亚洲高清| 十八禁网站网址无遮挡| 卡戴珊不雅视频在线播放| 亚洲人成77777在线视频| 欧美乱码精品一区二区三区| 国产黄色视频一区二区在线观看| 日韩av在线免费看完整版不卡| 精品少妇一区二区三区视频日本电影 | 欧美中文综合在线视频| 欧美激情 高清一区二区三区| 精品一区二区三区av网在线观看 | 搡老乐熟女国产| 国产一卡二卡三卡精品 | 欧美在线黄色| 国产精品欧美亚洲77777| 国产又色又爽无遮挡免| xxxhd国产人妻xxx| 90打野战视频偷拍视频| 欧美国产精品va在线观看不卡| 久久99一区二区三区| 亚洲欧美精品综合一区二区三区| 中文字幕色久视频| 男女无遮挡免费网站观看| 免费av中文字幕在线| 18禁国产床啪视频网站| 看免费av毛片| 99精国产麻豆久久婷婷| 亚洲欧美中文字幕日韩二区| 国产极品天堂在线| 三上悠亚av全集在线观看| 精品久久久久久电影网| 久久国产精品男人的天堂亚洲| av网站在线播放免费| 国产精品久久久人人做人人爽| 久久午夜综合久久蜜桃| 国产精品.久久久| av国产久精品久网站免费入址| 天堂中文最新版在线下载| a 毛片基地| 这个男人来自地球电影免费观看 | 成人国语在线视频| 成人午夜精彩视频在线观看| 看免费av毛片| 国产有黄有色有爽视频| 97精品久久久久久久久久精品| 成人手机av| 亚洲国产成人一精品久久久| 一本色道久久久久久精品综合| 亚洲国产看品久久| 久久综合国产亚洲精品| 久久久国产精品麻豆| 免费在线观看视频国产中文字幕亚洲 | 欧美97在线视频| 欧美少妇被猛烈插入视频| 国产伦理片在线播放av一区| 丝袜美足系列| 国产女主播在线喷水免费视频网站| 99热全是精品| 亚洲欧洲国产日韩| 无遮挡黄片免费观看| √禁漫天堂资源中文www| 欧美人与善性xxx| 黄色毛片三级朝国网站| 国产免费福利视频在线观看| 亚洲三区欧美一区| 久久这里只有精品19| 国产又爽黄色视频| 黑人欧美特级aaaaaa片| 国产男女超爽视频在线观看| a级毛片黄视频| 久久99精品国语久久久| 免费黄网站久久成人精品| 可以免费在线观看a视频的电影网站 | 成人国产av品久久久| 亚洲欧洲日产国产| 欧美日本中文国产一区发布| 精品国产国语对白av| 欧美精品一区二区大全| 老司机靠b影院| 国产成人免费观看mmmm| 国产又色又爽无遮挡免| 美女福利国产在线| 人成视频在线观看免费观看| 国产无遮挡羞羞视频在线观看| 菩萨蛮人人尽说江南好唐韦庄| 日韩电影二区| 久久99热这里只频精品6学生| 99国产综合亚洲精品| 欧美日韩av久久| 男人操女人黄网站| 欧美av亚洲av综合av国产av | 午夜福利在线免费观看网站| 一边摸一边抽搐一进一出视频| 中文欧美无线码| 女人高潮潮喷娇喘18禁视频| 欧美中文综合在线视频| 欧美激情极品国产一区二区三区| 国产精品熟女久久久久浪| 国产精品偷伦视频观看了| 国产在线一区二区三区精| 亚洲国产中文字幕在线视频| www.av在线官网国产| 热99国产精品久久久久久7| 国产黄频视频在线观看| 日日啪夜夜爽| 亚洲精品国产色婷婷电影| 成人免费观看视频高清| 欧美精品亚洲一区二区| 视频区图区小说| 欧美97在线视频| 色视频在线一区二区三区| 丝袜美腿诱惑在线| 精品少妇一区二区三区视频日本电影 | 亚洲精品av麻豆狂野| 亚洲欧美精品自产自拍| 婷婷成人精品国产| 中文字幕色久视频| 亚洲久久久国产精品| 女人爽到高潮嗷嗷叫在线视频| 视频在线观看一区二区三区| 国产在线一区二区三区精| 免费av中文字幕在线| 性高湖久久久久久久久免费观看| 水蜜桃什么品种好| 日韩不卡一区二区三区视频在线| 日韩电影二区| 国产成人免费无遮挡视频| 纵有疾风起免费观看全集完整版| 国产精品蜜桃在线观看| 高清av免费在线| 看免费成人av毛片| 下体分泌物呈黄色| 人人妻人人澡人人看| 久久精品熟女亚洲av麻豆精品| 国产av国产精品国产| 天天操日日干夜夜撸| 国产欧美日韩一区二区三区在线| 午夜久久久在线观看| 黄色一级大片看看| 欧美精品亚洲一区二区| 日本wwww免费看| 赤兔流量卡办理| 国产免费一区二区三区四区乱码| 婷婷色麻豆天堂久久| 高清在线视频一区二区三区| 成人影院久久| 一级毛片黄色毛片免费观看视频| 国产精品一二三区在线看| 日韩大片免费观看网站| 99国产精品免费福利视频| 一本大道久久a久久精品| 日本av免费视频播放| 国产xxxxx性猛交| 一区二区三区四区激情视频| av网站免费在线观看视频| 久久人人97超碰香蕉20202| 中文精品一卡2卡3卡4更新| 精品人妻熟女毛片av久久网站| 色婷婷久久久亚洲欧美| 观看av在线不卡| 秋霞伦理黄片| 久久韩国三级中文字幕| 亚洲精品日韩在线中文字幕| 在线观看人妻少妇| 青青草视频在线视频观看| 在线看a的网站| 午夜精品国产一区二区电影| 激情五月婷婷亚洲| 国产熟女欧美一区二区| 美女视频免费永久观看网站| 午夜激情久久久久久久| 国产免费现黄频在线看| 国产日韩欧美视频二区| 最新在线观看一区二区三区 | 一级毛片黄色毛片免费观看视频| 亚洲精品中文字幕在线视频| 2018国产大陆天天弄谢| 亚洲av电影在线进入| 亚洲五月色婷婷综合| 日韩大片免费观看网站| 久久鲁丝午夜福利片| 亚洲欧美精品综合一区二区三区| 免费不卡黄色视频| 91精品三级在线观看| 丝袜脚勾引网站| 性少妇av在线| 在线精品无人区一区二区三| 9热在线视频观看99| 日本vs欧美在线观看视频| 男女床上黄色一级片免费看| 久久久久精品性色| 久久久精品国产亚洲av高清涩受| 欧美老熟妇乱子伦牲交| 欧美国产精品va在线观看不卡| 日韩伦理黄色片| 1024香蕉在线观看| 日本欧美国产在线视频| 看非洲黑人一级黄片| 蜜桃在线观看..| av不卡在线播放| 69精品国产乱码久久久| 一区二区日韩欧美中文字幕| 美国免费a级毛片| 99香蕉大伊视频| 午夜福利一区二区在线看| 日韩免费高清中文字幕av| 如日韩欧美国产精品一区二区三区| 亚洲精品久久成人aⅴ小说| 欧美乱码精品一区二区三区| av免费观看日本| 久久天堂一区二区三区四区| 热99国产精品久久久久久7| 亚洲欧洲国产日韩| 人人妻人人爽人人添夜夜欢视频| 9热在线视频观看99|