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

    模擬常減壓蒸餾裝置的物性計(jì)算方法和模型結(jié)構(gòu)

    2018-06-29 07:26:50
    石油煉制與化工 2018年7期
    關(guān)鍵詞:減壓塔逸度減壓蒸餾

    鄭 文 剛

    (中國(guó)石油化工股份有限公司煉油事業(yè)部,北京 100728)

    隨著計(jì)算機(jī)軟硬件技術(shù)的高速發(fā)展,流程模擬技術(shù)已逐漸被大專院校、研究設(shè)計(jì)單位、企業(yè)各級(jí)管理部門及工藝工程師所掌握。中國(guó)石油化工股份有限公司自21世紀(jì)初開(kāi)始在各下屬企業(yè)推廣和應(yīng)用該技術(shù),培養(yǎng)了一批能熟練掌握和使用該技術(shù)的工程技術(shù)人員。目前流程模擬技術(shù)已在煉油廠裝置設(shè)計(jì)、工況分析、故障診斷、挖潛增效、加工方案優(yōu)化等方面發(fā)揮了重要作用,取得了顯著的經(jīng)濟(jì)效益。

    1 常減壓蒸餾裝置流程模擬

    長(zhǎng)期以來(lái),在初餾塔、常壓塔的塔盤效率估算,初餾塔、常壓塔和減壓塔進(jìn)料壓力、壓降估算,減壓爐轉(zhuǎn)油線、汽化段、減壓塔填料等板高度估算,重油的氣液平衡預(yù)測(cè)等方面一直存在問(wèn)題,加上現(xiàn)場(chǎng)原油評(píng)價(jià)、物料流量、溫度、壓力、分析數(shù)據(jù)等都存在一定誤差,因此難以準(zhǔn)確模擬常減壓蒸餾裝置。常減壓蒸餾裝置是煉油廠的“龍頭”,其穩(wěn)定和優(yōu)化運(yùn)行對(duì)全廠效益具有重要影響,因此有必要建立較準(zhǔn)確的模型,指導(dǎo)該裝置的設(shè)計(jì)、操作和優(yōu)化。

    2 常減壓蒸餾裝置模型結(jié)構(gòu)

    (1)總體思路。采用Aspen Plus V10.1軟件,應(yīng)用聯(lián)立方程(Equation Oriented,簡(jiǎn)稱EO)建模法建立某煉油廠的常減壓蒸餾裝置模型,然后選用不同的物性計(jì)算方法,建立特殊的減壓部分模型結(jié)構(gòu),調(diào)用動(dòng)態(tài)矩陣求解器(DMO),在數(shù)據(jù)整定(Reconciliation)模式下求優(yōu),使得模擬計(jì)算結(jié)果與標(biāo)定數(shù)據(jù)的誤差平方和最小。把各種物性計(jì)算方法求得的最小誤差平方和按從小到大排列,得出誤差最小的物性計(jì)算方法。該方法的優(yōu)點(diǎn)是能把原始標(biāo)定數(shù)據(jù)的誤差對(duì)模擬計(jì)算值的影響降低到最小程度,缺點(diǎn)是建模過(guò)程較復(fù)雜。

    (2)塔盤效率。為了能使模型中的塔盤和實(shí)際塔盤一一對(duì)應(yīng),初餾塔和常壓塔均采用實(shí)際塔盤進(jìn)行建模,塔盤效率采用Murphree效率,其最優(yōu)估計(jì)值采用EO建模法,用DMO求解器在Reconciliation求優(yōu)模式下求解得到。減壓塔塔盤數(shù)采用填料高度等板高度(HETP)折算成理論塔板數(shù),由于填料的HETP也和氣液相負(fù)荷有關(guān),所以按填料的設(shè)計(jì)HETP計(jì)算出理論塔板數(shù)后,還需要進(jìn)一步設(shè)置理論塔盤的Murphree效率。

    (3)初餾塔、常壓塔和減壓塔塔頂氣。一般情況下,模型預(yù)測(cè)的初餾塔、常壓塔塔頂氣流量均偏小,為了準(zhǔn)確模擬初餾塔和常壓塔塔頂瓦斯流量、回流罐溫度,需要在初餾塔和常壓塔進(jìn)料中各添加一股氣體,其組成和流量采用標(biāo)定數(shù)據(jù)。在EO建模法中,可以對(duì)這兩股氣體流量進(jìn)行數(shù)據(jù)整定,求解出最優(yōu)流量估計(jì)值。

    減壓塔塔頂氣一般流量很小,即使不做任何處理,通過(guò)模型也能準(zhǔn)確模擬減頂氣的溫度和流量。

    (4)初餾塔、常壓塔、減壓塔進(jìn)料壓力、全塔壓降估算。采用與塔盤效率估算相同的方法,即采用EO建模法中的數(shù)據(jù)整定模式求出最優(yōu)估算值。

    (5)減壓部分特殊的模型結(jié)構(gòu)。減壓爐轉(zhuǎn)油線在進(jìn)減壓塔的汽化段前,因?yàn)閴毫Σ粩嘟档停糠忠合嗥鼰嵋约稗D(zhuǎn)油線散熱導(dǎo)致溫度降低。由于轉(zhuǎn)油線直徑較大,液相停留在下部,氣相充滿上部,且兩相速度不一致,因此液相和氣相難以達(dá)到相平衡,造成氣相溫度高,液相溫度低。從減壓爐出口氣液兩相呈相平衡狀態(tài)到減壓塔汽化段入口氣液兩相呈不平衡狀態(tài)的過(guò)程中,存在某個(gè)壓力,在該壓力下,氣相和液相完成相分離,建模時(shí)應(yīng)把氣相單獨(dú)進(jìn)洗滌段,液相進(jìn)汽化段。計(jì)算表明,這種建模方法更準(zhǔn)確,同時(shí)建模時(shí)也需考慮在汽化段閃蒸過(guò)程中存在霧沫夾帶現(xiàn)象。如果不考慮該因素,減四線殘?zhí)亢徒饘俸棵黠@偏低。如減四線殘?zhí)繉?shí)際值為4%~14%,主要集中在7%~10%,而模擬計(jì)算值往往小于5%。建模時(shí),可把減壓塔分成汽化段、汽化段以上、汽化段以下三部分考慮,對(duì)汽化段應(yīng)考慮部分霧沫夾帶[1-4]。

    (6)序貫?zāi)K法(SM)和聯(lián)立方程法(EO)。SM法沿著流程圖走向,對(duì)單元設(shè)備逐個(gè)求解。遇到有循環(huán)的流量,則先設(shè)置物流初值,下游單元完成計(jì)算后,把循環(huán)物流計(jì)算值賦值給上次初始值,如此迭代,直到流程圖全部收斂。EO法把系統(tǒng)的所有單元設(shè)備的物料、能量平衡方程、約束方程、變量連接方程、邊界條件一起放在一個(gè)大型稀疏矩陣中,采用序貫二次型(SQP)法求解非線性方程組,從而得到流程圖中所有單元設(shè)備的物料和能量流股信息。EO法求解速率快,不僅能求解設(shè)計(jì)問(wèn)題,也能求解優(yōu)化問(wèn)題,還能給出參數(shù)的最優(yōu)估計(jì)值。但該方法需要較大的計(jì)算機(jī)內(nèi)存,且模型不收斂時(shí)不易查找問(wèn)題[5]。

    (7)換熱網(wǎng)絡(luò)。本課題目的是給出相對(duì)準(zhǔn)確的物性計(jì)算方法用于常減壓蒸餾裝置模擬,換熱網(wǎng)絡(luò)影響不大,因此簡(jiǎn)化了換熱網(wǎng)絡(luò)。模型中初餾塔、常壓塔和減壓塔的進(jìn)料溫度固定。

    (8)物性計(jì)算方法選用。物性計(jì)算方法是指一組熱力學(xué)、傳遞性質(zhì)計(jì)算模型,不是專指熱力學(xué)狀態(tài)方程(EOS),EOS只是其中的一部分。如Maxwell-Bonnell[6]物性計(jì)算方法,氣液平衡常數(shù)用Maxwell-Bonnell關(guān)聯(lián)式計(jì)算,氣體和液體的焓、熵、Gibbs自由能、氣體和液體的摩爾體積用Soave-Redlich-Kwong[7]熱力學(xué)狀態(tài)方程計(jì)算。在Aspen Plus中物性計(jì)算方法統(tǒng)稱為性質(zhì)模型。由于不涉及換熱網(wǎng)絡(luò)、塔盤水力學(xué)核算、吸附分離、催化反應(yīng)等,因此不討論這些模型。

    對(duì)于常減壓蒸餾裝置,特別是減壓塔,各種物性方法的模擬結(jié)果差別較大。在實(shí)踐中,經(jīng)常遇到減壓渣油的模擬收率和實(shí)際標(biāo)定值有差距,一般偏低。最大的可能是減壓部分沒(méi)有選擇合適的物性計(jì)算方法。以下探討在Aspen Plus中用于常減壓蒸餾裝置模擬的準(zhǔn)確物性計(jì)算方法。

    在Aspen Plus中選用不同的物性計(jì)算方法對(duì)某常減壓蒸餾裝置進(jìn)行模擬。其中初餾塔和常壓塔選用一種方法,減壓部分(含減壓爐、減壓塔)選用另一種方法,并在BK10[8]、Maxwell-Bonnell(簡(jiǎn)稱MX)、Chao-Seader[9](簡(jiǎn)稱CS)、Grayson-Streed[10](簡(jiǎn)稱GS)、PENG-ROB或Peng-Robinson[11](簡(jiǎn)稱PR)、Redlich-Kwong-SOAVE(簡(jiǎn)稱RKS)、Soave-Redlich-Kwong (簡(jiǎn)稱SRK)模型之間進(jìn)行組合。

    BK10和MX法屬于蒸氣壓-液相逸度模型,CS和GS法屬于氣相狀態(tài)方程-液相活度系數(shù)模型,PR、RKS和SRK法屬于狀態(tài)方程模型。這些方法都是一系列模型的組合,比如MX方法使用理想氣體狀態(tài)方程計(jì)算氣相逸度系數(shù),Maxwell-Bonnell關(guān)聯(lián)公式計(jì)算飽和蒸氣壓,從而得到氣液平衡常數(shù)、Soave-Redlich-Kwong方程計(jì)算氣體和液體的焓、熵、Gibbs自由能、氣體和液體的摩爾體積。

    3 EO建模法下Aspen Plus不同物性計(jì)算方法對(duì)常減壓蒸餾裝置模擬結(jié)果的影響

    任意挑選一套某煉油廠的常減壓蒸餾裝置標(biāo)定數(shù)據(jù),含原油評(píng)價(jià)、主要操作參數(shù)、物料平衡和產(chǎn)品分析數(shù)據(jù)。用Aspen Plus的EO建模法,建立了該裝置的模型。

    不僅物性計(jì)算方法會(huì)影響模擬結(jié)果,原油評(píng)價(jià)數(shù)據(jù)、塔盤效率、進(jìn)料壓力、塔壓降、流量誤差、溫度誤差、分析誤差等也會(huì)影響模擬結(jié)果。為了嚴(yán)格比較物性計(jì)算方法的影響,需要剔除干撓因素或?qū)⑵溆绊懣刂圃谧钚》秶鷥?nèi)。因而采用如下方法:選定物性計(jì)算方法后,在數(shù)據(jù)整定模式下對(duì)模擬結(jié)果與標(biāo)定數(shù)據(jù)的誤差平方和求最小值。把多個(gè)物性計(jì)算方法對(duì)應(yīng)的目標(biāo)函數(shù)從小到大排列,就可得出最準(zhǔn)確的方法。

    以初餾塔為例,原油中的水分、進(jìn)料壓力、塔壓降、塔盤效率、塔頂石腦油流量、塔頂氣體流量是自變量。初餾塔塔頂石腦油干點(diǎn)、回流罐溫度、塔底溫度和塔頂回流量是因變量。這些變量除塔盤效率外,都存在測(cè)量值。在EO建模法的Reconciliation求解模式下,把關(guān)鍵變量的實(shí)測(cè)值和模擬值的誤差平方和作為目標(biāo)函數(shù)。并在裝置各種約束范圍內(nèi),求目標(biāo)函數(shù)(OF)的最小值。

    (1)

    式中:Xim為某變量的實(shí)測(cè)值;Xis為該變量的模擬值;i為標(biāo)準(zhǔn)偏差。

    建立的EO模型,自由變量總數(shù)是33個(gè),目標(biāo)函數(shù)中有37個(gè)關(guān)鍵變量,其中13個(gè)屬于自由變量。模型的收斂時(shí)間為500~1 500 s,大部分情況下在1 200 s之內(nèi)可收斂求得優(yōu)解。不同物性方法組合對(duì)常減壓蒸餾裝置模擬結(jié)果的影響見(jiàn)表1。

    表1 Aspen Plus不同物性計(jì)算方法組合對(duì)常減壓蒸餾裝置模擬結(jié)果的影響

    從表1可見(jiàn),誤差平方和最小的模型是初餾塔、常壓塔部分用BK10,減壓部分用Maxwell-Bonnell方法。(BK10+BK10)方法的誤差是(BK10+MX)方法的1.05倍,(PR+PR)方法的誤差是(BK10+MX)方法的2.96倍??梢?jiàn)用正確選用物性計(jì)算方法對(duì)于常減壓蒸餾裝置的模擬十分重要。從表1還可以看出,初餾塔和常壓塔部分的物性計(jì)算方法按BK10,MX,GS,CS,RKS,SRK,PR的順序,誤差逐漸遞增;減壓部分的物性計(jì)算方法按MX,BK10,CS,GS,RKS,SRK,PR的順序,誤差逐漸遞增;物性計(jì)算方法的使用范圍越寬,如RKS、SRK和PR,準(zhǔn)確度越低;方法的針對(duì)性越強(qiáng),準(zhǔn)確度就越高;越是有物理意義的模型,如PR,不如回歸模型準(zhǔn)確,如MX方法適于低壓下的重組分(尤其適合減壓深拔),因此模擬減壓部分越準(zhǔn);CS和GS法的差別較小,只是在液相的活度系數(shù)、Gibbs自由能和熵計(jì)算方法上有區(qū)別。

    以下分初餾塔、常壓塔和減壓塔三部分,從表1的34組方法中挑選出最重要的方法組合,分析其對(duì)初餾塔、常壓塔、減壓塔模擬結(jié)果的具體影響,結(jié)果分別見(jiàn)表2、表3和表4。表中所列物性計(jì)算方法,如PR MX表示初餾塔和常壓塔用Peng-Robinson方法,減壓塔用Maxwell Bonnell方法。

    由表2可見(jiàn),初餾塔的進(jìn)料壓力誤差較大,初餾塔的塔盤效率在70%左右,各方法對(duì)初餾塔均能較準(zhǔn)確模擬,差別不大。

    由表3可見(jiàn),常壓塔塔頂循環(huán)抽出溫度、一中段抽出溫度、塔頂冷回流量、塔頂石腦油、常一線、二線和三線流量模型預(yù)測(cè)值和實(shí)際有偏差,特別是冷回流和常一線流量。由于采樣時(shí)間不同步等原因,實(shí)際化驗(yàn)分析數(shù)據(jù)可能不具有代表性。如提高常一線終餾點(diǎn),常一線流量模擬值將增加,常二線流量將減小,計(jì)算值和標(biāo)定數(shù)據(jù)將更為接近。上述結(jié)果也表明煉油廠的標(biāo)定數(shù)據(jù)不完全是可信的,在EO建模法下,可以鑒別數(shù)據(jù)的合理性。

    從表4可見(jiàn),MX和BK10模型對(duì)于預(yù)測(cè)減壓渣油、減四線流量較準(zhǔn)確,MX模型更為準(zhǔn)確,這與表1結(jié)果相符。模型預(yù)測(cè)結(jié)果表明,汽化段有2%~5%霧沫夾帶量,這對(duì)于減四線性質(zhì)模擬的準(zhǔn)確性有較大幫助。同時(shí)是否考慮轉(zhuǎn)油線的相分離對(duì)于減壓塔模擬的準(zhǔn)確度也有明顯影響。

    表2 不同物性計(jì)算方法對(duì)初餾塔模擬結(jié)果的影響

    表3 不同物性計(jì)算方法對(duì)常壓塔模擬結(jié)果的影響

    表4 不同物性計(jì)算方法對(duì)減壓塔模擬結(jié)果的影響

    4 不同方法影響模擬結(jié)果的理論解釋

    多元體系氣液兩相平衡要滿足氣相逸度等于液相逸度[12]。液相逸度有兩種計(jì)算方法:一是活度系數(shù)法,Maxwell Bonnell、BK10、Grayson-Streed和Chao-Sea采用該方法;二是液相逸度用逸度系數(shù)計(jì)算,Peng-Robinson(PR)、Soave-Redlich-Kwong(SRK)和Redlich-Kwong-Soave(RKS)采用該方法。分離過(guò)程計(jì)算準(zhǔn)確度主要受相平衡常數(shù)K是否準(zhǔn)確的影響。RKS和SRK的PVT關(guān)系式是相同的,只是液體摩爾體積的計(jì)算方法有所不同。RKS采用API方法,SRK采用Peneloux體積校正因子。

    GS和CS方法中氣相逸度系數(shù)兩者均采用RK狀態(tài)方程,液相活度系數(shù)的計(jì)算都采用Scatchard-Hildebrand正規(guī)溶液模型。

    在減壓塔汽化段溫度365 ℃、壓力5 kPa的條件下,以減壓渣油中平均沸點(diǎn)為579 ℃的虛擬組分為例,應(yīng)用不同方法計(jì)算其K值,計(jì)算結(jié)果見(jiàn)表5。K值與液相活度系數(shù)、純組分飽和蒸氣壓、純組分飽和蒸氣逸度系數(shù)、氣相逸度系數(shù)有關(guān)。但從表5的計(jì)算結(jié)果可見(jiàn),K值的大小主要與飽和蒸氣壓有關(guān),且呈線性關(guān)系,相關(guān)系數(shù)達(dá)0.999 9,而其它系數(shù)幾乎都等于1,影響很小。

    從表5可見(jiàn),RKS和SRK法計(jì)算的K值相同,但模型誤差不同。這是因?yàn)樗钠胶夥匠讨?,除了相平衡方程外,還有能量平衡方程(焓平衡),而液相焓計(jì)算需要用到液相摩爾體積,RKS和SRK的液相摩爾體積計(jì)算方法不同,從而導(dǎo)致最終計(jì)算結(jié)果不同。

    表5 減壓條件下不同物性計(jì)算方法所得相平衡常數(shù)

    常壓塔的計(jì)算類似,計(jì)算結(jié)果表明K值的大小主要與飽和蒸氣壓有關(guān),但并不呈線性關(guān)系,常壓下K值與液相活度系數(shù)、純組分飽和蒸氣壓、純組分飽和蒸氣逸度系數(shù)、氣相逸度系數(shù)均有關(guān)。

    對(duì)比常壓和減壓下的K值,發(fā)現(xiàn)采用不同的物性計(jì)算方法對(duì)減壓部分模擬時(shí),K值相差較大,如PR法的K值是MX的2.74倍;不同的物性計(jì)算方法對(duì)常壓部分影響不大,如PR法的K值是BK10法的1.06倍。這與多套裝置的模擬經(jīng)驗(yàn)完全相符。

    (1)石油窄餾分的飽和蒸氣壓。應(yīng)用API 5A1.19 Maxwell Bonnell關(guān)聯(lián)公式計(jì)算579 ℃的虛擬組分在365 ℃下的飽和蒸氣壓,得0.327 kPa,又從文獻(xiàn)[6]可知,該預(yù)測(cè)值比實(shí)際值高。

    在Aspen Plus中MX方法計(jì)算的579 ℃虛擬組分在365 ℃下的飽和蒸氣壓是0.337 kPa,而實(shí)際值小于0.327 kPa。換言之,與實(shí)際值相比,用MX法預(yù)測(cè)時(shí)K值較大,常壓渣油更易揮發(fā),計(jì)算得到的減壓渣油收率會(huì)偏低,這和實(shí)際經(jīng)驗(yàn)吻合。從表5可見(jiàn),K值按MX,BK10,CS,GS,RKS,SRK,PR依次增大,這和表1的減壓部分模型按MX,BK10,CS,GS,RKS,SRK,PR的順序,誤差逐漸遞增一致。

    與實(shí)際值相比,在Aspen Plus中BK10法預(yù)測(cè)的飽和蒸氣壓數(shù)據(jù)比MX法偏差大。BK10法回歸自大量相平衡測(cè)量數(shù)據(jù),文獻(xiàn)[9]指出BK10法所用原始數(shù)據(jù)中壓力下限是68.95 kPa。由此可見(jiàn),在減壓塔和減壓爐轉(zhuǎn)油線的壓力范圍內(nèi),K值是外推的,因而誤差比MX法大。但是在初餾塔和常壓塔部分,由于BK10法直接根據(jù)實(shí)驗(yàn)值回歸出K值,因此其誤差最小。

    CS,GS,RKS,SKR,PR方法采用不同的立方型狀態(tài)方程模擬飽和蒸氣壓,但在低壓和低溫下,如平均沸點(diǎn)為579 ℃的虛擬組分在365 ℃下預(yù)測(cè)的飽和蒸氣壓數(shù)據(jù)偏差仍較大[13]。

    綜上所述,由于MX和BK10法是基于低壓下實(shí)測(cè)數(shù)據(jù)回歸的蒸氣壓方程或相平衡常數(shù)方程,因此這兩個(gè)方法預(yù)測(cè)的飽和蒸氣壓最準(zhǔn)確。CS,GS,RKS,SRK,PR法是基于立方形狀態(tài)方程預(yù)測(cè)飽和蒸氣壓數(shù)據(jù),而立方型方程本身就較難準(zhǔn)確預(yù)測(cè)飽和蒸氣壓,因此誤差較大。

    (2)石油及常壓渣油溶液的理想性。Hariu等[14]證明一般石油餾分在壓力低于300 kPa條件下基本符合理想體系,石油餾分中某一窄餾分的相平衡常數(shù)可按Raoult定律求取。其液相活度系數(shù)、飽和蒸氣逸度系數(shù)、氣相逸度系數(shù)均可為1。表明對(duì)于處于低壓的常壓渣油溶液,氣相和液相都可以近似作為理想體系和理想溶液處理。

    CS和GS法把液相按正規(guī)溶液處理,計(jì)算出活度系數(shù)大于1,CS,GS,RKS,SRK,PR法用狀態(tài)方程計(jì)算氣相,計(jì)算出氣相逸度系數(shù)小于1,這些都小幅增加了氣液平衡常數(shù)。

    由于GS法在CS法基礎(chǔ)上對(duì)計(jì)算純組分液相逸度系數(shù)的參數(shù)進(jìn)行了更新,使GS法更適合在溫度和壓力提高后相平衡常數(shù)的預(yù)測(cè)。因此在減壓部分,CS法誤差小,而在常壓部分,GS法誤差小。

    從表5可見(jiàn),如果要更準(zhǔn)確地模擬減壓部分,需要進(jìn)一步降低相平衡常數(shù),液相的活度系數(shù)則應(yīng)小于1,常壓渣油可能是較弱的負(fù)偏差溶液。其原因可能是原油中的一些高分子極性化合物,如膠質(zhì)、瀝青質(zhì)、稠環(huán)芳烴等,由于偶極相互作用、電荷轉(zhuǎn)移π-π鍵作用和氫鍵作用,彼此發(fā)生締合,形成“超分子結(jié)構(gòu)”。這些超分子結(jié)構(gòu)外表面的分子帶有過(guò)剩的能量,形成一個(gè)附加的引力場(chǎng),能夠吸引一些分子量較小、芳香度較低的烴類。在原油蒸餾過(guò)程中,這部分烴類,由于受膠核附加引力的作用,在達(dá)到沸點(diǎn)時(shí)難以轉(zhuǎn)入氣相而滯留于重油中[15]。

    5 結(jié) 論

    (1)對(duì)于常減壓蒸餾裝置模擬,在Aspen Plus中推薦采用BK10法模擬常壓部分、Maxwell-Bonnell法模擬減壓部分,模擬結(jié)果與實(shí)際值最接近,且計(jì)算速率快。BK10法的計(jì)算誤差是推薦方法的1.045倍。減壓部分不能應(yīng)用立方型狀態(tài)方程預(yù)測(cè)氣液平衡常數(shù),誤差大,計(jì)算速率慢。

    (2)除正確選擇物性計(jì)算方法外,常減壓蒸餾裝置的模型結(jié)構(gòu)也很重要。減壓部分模擬需要考慮減壓爐轉(zhuǎn)油線的兩相分離以及減壓塔汽化段的霧沫夾帶,否則減四線的誤差較大。

    (3)從理論上解釋了不同物性計(jì)算方法對(duì)常減壓蒸餾裝置模擬結(jié)果的影響。表明在減壓環(huán)境下,常壓渣油不是正規(guī)溶液,常壓渣油可能是弱的負(fù)偏差溶液。

    (4)應(yīng)用EO建模法和理論分析,可以快速確定原油評(píng)價(jià)數(shù)據(jù)、各種工藝參數(shù)、化驗(yàn)分析數(shù)據(jù)的誤差對(duì)模擬結(jié)果的影響,估算塔盤效率、塔盤壓降、塔進(jìn)料壓力等不可測(cè)量或無(wú)測(cè)量的參數(shù),以益于裝置診斷、優(yōu)化。

    參考文獻(xiàn)

    [1] Remesat D.Improving crude vacuum unit performance[J].Petroleum Technology Quarterly,2008(3):107-113

    [2] Pilling M,Roza M,Wong S M.Entrainment issues in vacuum column flash zones[J].Digital Refining,2010(1):1-10

    [3] Golden S,Barletta T,White S.Vacuum unit performance[J].Digital Refining,2012(4):1-6

    [4] Golden S,Barletta T,White S.Myth of high cutpoint in dry vacuum units[J].Digital Refining 2014(2):1-7

    [5] 王健紅,馮樹(shù)波,杜增智,等.化工系統(tǒng)工程理論與實(shí)踐[M].北京:化學(xué)工業(yè)出版社,2009:119-124

    [6] Maxwell J B,Bonnell L S.Deviation and precision of a new pressure correlation for petroleum hydrocarbons[J].Ind Eng Chem,1957,49(7):1187-1196

    [7] Soave G.Equilibrium constants from a modified Redlich-Kwong equation of state[J].Chem Eng Sci,1972,27(6):1197-1203

    [8] Cajander B C,Hipkin H G,Lenoir J M.The prediction of equilibrium ratios—A nomogram of improved accuracy[J].J Chem Eng Data,1960,5(3):251-259

    [9] Chao K C,Seader J D.A general correlation of vapor-liquid equilibria in hydrocarbon mixtures[J].AIChE J,1961,7(4):598-605

    [10] Grayson H G,Streed C W.Vapor-liquid equilibria for high temperature,high pressure hydrogen-hydrocarbon systems[C]Sixth World Petroleum Conference.Paper No.20 PD7.Frankfurt,1963-06

    [11] Peng D Y,Robinson D B.A new two-constant equation of state[J].Industrial and Engineering Chemistry Fundamental,1976,15(1):59-64

    [12] 郭天民.多元汽-液平衡和精餾[M].北京:化學(xué)工業(yè)出版社,1983:14-19

    [13] 汪萍.立方型狀態(tài)方程性能的評(píng)價(jià)研究[D].青島:青島科技大學(xué),2004

    [14] Hariu O H,Sage R C.Crude split figured by computer[J].Hydrocarbon Processing,1969,48(4):143

    [15] 梁文杰.石油化學(xué)[M].東營(yíng):石油大學(xué)出版社,1995:64-65

    猜你喜歡
    減壓塔逸度減壓蒸餾
    煤制油行業(yè)產(chǎn)品提質(zhì)項(xiàng)目創(chuàng)新探索
    山西化工(2023年12期)2024-01-12 01:51:38
    常減壓蒸餾裝置的換熱網(wǎng)絡(luò)綜合設(shè)計(jì)及優(yōu)化研究
    山西化工(2022年2期)2023-01-14 05:20:05
    地幔氧逸度的時(shí)空變化
    減壓塔塔底熱油泵管道設(shè)計(jì)淺析
    云南化工(2020年11期)2021-01-14 00:51:04
    常減壓塔頂污水鐵離子氯離子問(wèn)題分析
    云南化工(2020年11期)2021-01-14 00:51:00
    不同溫度、壓強(qiáng)、氧逸度條件下斜方輝石含水性的實(shí)驗(yàn)研究
    幔源巖漿氧化還原狀態(tài)及對(duì)巖漿礦床成礦的制約*
    減壓蒸餾應(yīng)用于粗苯工段的可行性分析
    山東冶金(2018年6期)2019-01-28 08:14:52
    氣體混合爐中氧逸度控制
    水噴射真空泵水箱降溫方案
    亚洲内射少妇av| 午夜精品一区二区三区免费看| 国产亚洲欧美98| 午夜福利成人在线免费观看| 亚洲国产精品国产精品| 精品久久久久久久久久免费视频| 欧美变态另类bdsm刘玥| 三级经典国产精品| 1000部很黄的大片| 2022亚洲国产成人精品| 搡女人真爽免费视频火全软件| 亚洲美女视频黄频| 亚洲av.av天堂| 精品久久久久久久久亚洲| 久久这里只有精品中国| 久久久久久久久大av| 长腿黑丝高跟| 中文字幕熟女人妻在线| 综合色丁香网| 卡戴珊不雅视频在线播放| 高清毛片免费观看视频网站| 午夜免费激情av| 国产成人aa在线观看| 午夜精品国产一区二区电影 | 国产成人一区二区在线| 91aial.com中文字幕在线观看| 久久久久久大精品| 日韩一区二区三区影片| 欧美+日韩+精品| 久久久午夜欧美精品| 国产一区二区激情短视频| 1024手机看黄色片| 欧美最黄视频在线播放免费| 久久热精品热| 99九九线精品视频在线观看视频| 国产精品久久久久久久久免| 老熟妇乱子伦视频在线观看| 少妇人妻一区二区三区视频| 免费搜索国产男女视频| 中文在线观看免费www的网站| 国产白丝娇喘喷水9色精品| 黄片wwwwww| 国产黄片视频在线免费观看| 九九爱精品视频在线观看| 国产免费男女视频| 亚洲欧洲国产日韩| 亚洲婷婷狠狠爱综合网| 国产av一区在线观看免费| www.av在线官网国产| 午夜福利在线在线| 搡老妇女老女人老熟妇| 国产精品蜜桃在线观看 | 一边亲一边摸免费视频| 最近的中文字幕免费完整| 成人综合一区亚洲| 久久欧美精品欧美久久欧美| 一级二级三级毛片免费看| 韩国av在线不卡| 少妇的逼水好多| 一进一出抽搐gif免费好疼| 国产一区亚洲一区在线观看| 1024手机看黄色片| 97超碰精品成人国产| 91午夜精品亚洲一区二区三区| 又粗又硬又长又爽又黄的视频 | 色综合色国产| 久久久久久伊人网av| 一边摸一边抽搐一进一小说| 黄色欧美视频在线观看| 此物有八面人人有两片| 18禁在线播放成人免费| 校园春色视频在线观看| 高清午夜精品一区二区三区 | 亚洲成人久久爱视频| 啦啦啦观看免费观看视频高清| 精品一区二区三区人妻视频| 国产精品三级大全| 国产日本99.免费观看| av专区在线播放| 大型黄色视频在线免费观看| 久久综合国产亚洲精品| 亚洲一区高清亚洲精品| 三级经典国产精品| 国产精品人妻久久久影院| 69人妻影院| 18禁在线无遮挡免费观看视频| 久久热精品热| av在线观看视频网站免费| 日韩欧美一区二区三区在线观看| 丰满人妻一区二区三区视频av| 国产精品精品国产色婷婷| 18禁裸乳无遮挡免费网站照片| 亚洲图色成人| 日本免费一区二区三区高清不卡| а√天堂www在线а√下载| 午夜福利成人在线免费观看| av福利片在线观看| 亚洲内射少妇av| 久久人人爽人人片av| 欧美+日韩+精品| 在线观看一区二区三区| 性色avwww在线观看| 久久久久久久亚洲中文字幕| 国产精品久久久久久亚洲av鲁大| 我要搜黄色片| av黄色大香蕉| 日日啪夜夜撸| 日韩成人av中文字幕在线观看| 婷婷色av中文字幕| 你懂的网址亚洲精品在线观看 | 丰满人妻一区二区三区视频av| 久久久国产成人免费| 99久久精品一区二区三区| 亚洲真实伦在线观看| 久久综合国产亚洲精品| 久久久a久久爽久久v久久| 成人毛片a级毛片在线播放| 淫秽高清视频在线观看| 亚洲人成网站在线播| 国产高清不卡午夜福利| 久久精品综合一区二区三区| 一卡2卡三卡四卡精品乱码亚洲| av天堂中文字幕网| 黑人高潮一二区| 国产真实乱freesex| 久久亚洲国产成人精品v| 欧美日韩精品成人综合77777| av卡一久久| 国产高清三级在线| 最好的美女福利视频网| 青青草视频在线视频观看| 麻豆精品久久久久久蜜桃| 成人高潮视频无遮挡免费网站| 国产蜜桃级精品一区二区三区| 精品欧美国产一区二区三| 欧美变态另类bdsm刘玥| 国产精品久久久久久精品电影| 国内精品宾馆在线| 免费看光身美女| 国产蜜桃级精品一区二区三区| 久久午夜亚洲精品久久| 亚洲欧美日韩东京热| www.色视频.com| 久久久久九九精品影院| 91精品一卡2卡3卡4卡| 3wmmmm亚洲av在线观看| 国产精品美女特级片免费视频播放器| www.av在线官网国产| 国产精品免费一区二区三区在线| 免费看av在线观看网站| 亚洲一区二区三区色噜噜| 亚洲无线在线观看| 亚洲国产精品国产精品| 2021天堂中文幕一二区在线观| 99久久无色码亚洲精品果冻| 亚洲人成网站在线观看播放| 热99在线观看视频| 性色avwww在线观看| 成人午夜精彩视频在线观看| 亚洲成人中文字幕在线播放| 亚洲欧美日韩高清专用| 国产精品国产高清国产av| 爱豆传媒免费全集在线观看| 欧美3d第一页| 成人毛片a级毛片在线播放| 在线观看美女被高潮喷水网站| 国产久久久一区二区三区| 色哟哟哟哟哟哟| 超碰av人人做人人爽久久| 别揉我奶头 嗯啊视频| 悠悠久久av| 日本av手机在线免费观看| 中文字幕人妻熟人妻熟丝袜美| 12—13女人毛片做爰片一| 亚洲av一区综合| 又爽又黄a免费视频| 日韩人妻高清精品专区| 精品久久久久久久末码| 一级毛片aaaaaa免费看小| 黄片无遮挡物在线观看| 国产精品嫩草影院av在线观看| 中国美白少妇内射xxxbb| 99久久精品热视频| 国产黄色视频一区二区在线观看 | 国产不卡一卡二| 久久九九热精品免费| 欧美日韩国产亚洲二区| 深夜精品福利| 精品日产1卡2卡| 精品熟女少妇av免费看| 校园人妻丝袜中文字幕| 欧美潮喷喷水| 九九久久精品国产亚洲av麻豆| 国产精品美女特级片免费视频播放器| 99视频精品全部免费 在线| 久久久精品94久久精品| 国产在线精品亚洲第一网站| 久久精品国产自在天天线| 在线免费十八禁| 国产精品乱码一区二三区的特点| 国产一区二区亚洲精品在线观看| 国产大屁股一区二区在线视频| 久久久午夜欧美精品| 一区二区三区免费毛片| 成人永久免费在线观看视频| 又粗又硬又长又爽又黄的视频 | 亚洲真实伦在线观看| 日韩欧美精品免费久久| 亚洲最大成人中文| 国产一区二区三区在线臀色熟女| 尤物成人国产欧美一区二区三区| 日本与韩国留学比较| 亚洲欧美清纯卡通| 国产毛片a区久久久久| 床上黄色一级片| 精品日产1卡2卡| 成人毛片a级毛片在线播放| 日韩欧美 国产精品| 国产精品爽爽va在线观看网站| 精品久久久久久成人av| 久久热精品热| 日本欧美国产在线视频| 黄色配什么色好看| 黄色视频,在线免费观看| 日韩欧美国产在线观看| 最新中文字幕久久久久| 国产片特级美女逼逼视频| 日本一本二区三区精品| 日本与韩国留学比较| 日本av手机在线免费观看| 久久精品国产99精品国产亚洲性色| 你懂的网址亚洲精品在线观看 | 老司机福利观看| 我的老师免费观看完整版| а√天堂www在线а√下载| 一级毛片aaaaaa免费看小| 久久久午夜欧美精品| 麻豆国产97在线/欧美| 日韩 亚洲 欧美在线| 在线播放无遮挡| 亚洲精品成人久久久久久| 亚洲国产色片| 小说图片视频综合网站| 99国产精品一区二区蜜桃av| 成人永久免费在线观看视频| 91久久精品国产一区二区三区| 在线播放国产精品三级| 亚洲精品成人久久久久久| 亚洲成人久久爱视频| 小说图片视频综合网站| 国产成人a∨麻豆精品| 成人午夜精彩视频在线观看| 长腿黑丝高跟| 国产成人a区在线观看| 国产一区二区在线观看日韩| 国产69精品久久久久777片| 色哟哟哟哟哟哟| 毛片一级片免费看久久久久| 久久久久久久久久成人| 欧美性猛交╳xxx乱大交人| 国内精品久久久久精免费| 最新中文字幕久久久久| 亚洲欧美精品综合久久99| 欧美一区二区精品小视频在线| 成人无遮挡网站| 熟妇人妻久久中文字幕3abv| 噜噜噜噜噜久久久久久91| 亚洲在线观看片| 久久精品夜夜夜夜夜久久蜜豆| 午夜激情福利司机影院| 国产爱豆传媒在线观看| 亚洲国产精品成人久久小说 | 国产麻豆成人av免费视频| 99久久九九国产精品国产免费| 久久久久久久久久黄片| 六月丁香七月| 欧美性猛交╳xxx乱大交人| 精品无人区乱码1区二区| 亚洲色图av天堂| 深夜精品福利| 啦啦啦啦在线视频资源| 两性午夜刺激爽爽歪歪视频在线观看| 哪里可以看免费的av片| 又爽又黄无遮挡网站| 白带黄色成豆腐渣| 国产精品一区二区三区四区久久| 亚洲人与动物交配视频| 国产精品1区2区在线观看.| 欧美性猛交黑人性爽| 尤物成人国产欧美一区二区三区| 国产极品天堂在线| 性色avwww在线观看| 午夜视频国产福利| 美女国产视频在线观看| 麻豆av噜噜一区二区三区| 插逼视频在线观看| 综合色av麻豆| 一级黄片播放器| 日日摸夜夜添夜夜添av毛片| 国产亚洲5aaaaa淫片| 亚洲欧洲日产国产| 日韩欧美精品免费久久| 日本一本二区三区精品| 国产三级中文精品| 一夜夜www| 国产精品久久电影中文字幕| 午夜福利高清视频| 久久九九热精品免费| 小说图片视频综合网站| 一级毛片我不卡| 色综合色国产| 一个人观看的视频www高清免费观看| 国产成人一区二区在线| 免费观看人在逋| 国产在线男女| 丰满乱子伦码专区| 在线国产一区二区在线| 国产国拍精品亚洲av在线观看| 国产 一区 欧美 日韩| 99久久精品热视频| 国产在视频线在精品| 精品久久久久久久久亚洲| 女人被狂操c到高潮| 亚洲真实伦在线观看| 狂野欧美白嫩少妇大欣赏| 国产成人精品久久久久久| 久久久久九九精品影院| 夫妻性生交免费视频一级片| 最后的刺客免费高清国语| 国产白丝娇喘喷水9色精品| 成人午夜高清在线视频| 黄色一级大片看看| 夜夜爽天天搞| 波野结衣二区三区在线| 精品久久久久久久末码| 日本成人三级电影网站| 一本久久精品| 国产极品精品免费视频能看的| 能在线免费看毛片的网站| 黄色日韩在线| 六月丁香七月| 亚洲人成网站在线播放欧美日韩| 欧美xxxx黑人xx丫x性爽| 国产精品1区2区在线观看.| 免费搜索国产男女视频| 三级男女做爰猛烈吃奶摸视频| 日日摸夜夜添夜夜爱| 国产免费男女视频| 一夜夜www| 国产探花极品一区二区| 丝袜美腿在线中文| 国语自产精品视频在线第100页| 搡老妇女老女人老熟妇| 少妇高潮的动态图| videossex国产| 波野结衣二区三区在线| 高清毛片免费观看视频网站| 国内精品美女久久久久久| 男人和女人高潮做爰伦理| av.在线天堂| 久久国内精品自在自线图片| 悠悠久久av| 亚洲国产精品成人久久小说 | 亚洲欧美清纯卡通| 色哟哟哟哟哟哟| 在现免费观看毛片| 男的添女的下面高潮视频| 能在线免费观看的黄片| 99国产极品粉嫩在线观看| 日韩一区二区三区影片| 少妇人妻一区二区三区视频| 国产91av在线免费观看| 日韩欧美一区二区三区在线观看| 国产高清三级在线| 麻豆国产97在线/欧美| 亚洲丝袜综合中文字幕| 久久久精品94久久精品| 国产伦一二天堂av在线观看| 国产黄色小视频在线观看| 久久鲁丝午夜福利片| 一区二区三区四区激情视频 | 国产老妇伦熟女老妇高清| 国产免费一级a男人的天堂| 日韩亚洲欧美综合| 中文字幕av成人在线电影| 熟女人妻精品中文字幕| 91午夜精品亚洲一区二区三区| 少妇人妻精品综合一区二区 | 麻豆成人午夜福利视频| 久久热精品热| 啦啦啦啦在线视频资源| 久久欧美精品欧美久久欧美| 色视频www国产| 少妇人妻一区二区三区视频| 欧美高清性xxxxhd video| 九九久久精品国产亚洲av麻豆| 精品久久久久久久久久久久久| 国产精品一及| 国产乱人偷精品视频| 色综合亚洲欧美另类图片| 国产又黄又爽又无遮挡在线| 国产av不卡久久| 成人午夜精彩视频在线观看| 免费大片18禁| 亚洲色图av天堂| 日韩欧美精品免费久久| 99久久久亚洲精品蜜臀av| 欧美xxxx性猛交bbbb| 此物有八面人人有两片| 国产精品久久久久久久久免| 精品一区二区免费观看| 国产色婷婷99| 青春草亚洲视频在线观看| 熟妇人妻久久中文字幕3abv| 最近中文字幕高清免费大全6| 亚洲国产精品合色在线| 成人无遮挡网站| 亚洲中文字幕日韩| 国产精品美女特级片免费视频播放器| 亚洲乱码一区二区免费版| 如何舔出高潮| 国产私拍福利视频在线观看| 日韩一区二区视频免费看| 舔av片在线| 麻豆成人av视频| 免费人成视频x8x8入口观看| 成人特级黄色片久久久久久久| 久久久国产成人免费| 欧美又色又爽又黄视频| 午夜精品国产一区二区电影 | 中文字幕免费在线视频6| 丰满的人妻完整版| 成人一区二区视频在线观看| 一级毛片久久久久久久久女| 国产精品1区2区在线观看.| 男人的好看免费观看在线视频| 99久久久亚洲精品蜜臀av| 久久精品国产亚洲网站| 国产不卡一卡二| 中文字幕人妻熟人妻熟丝袜美| 中文精品一卡2卡3卡4更新| 国产在线男女| av视频在线观看入口| 欧美区成人在线视频| 国产精品久久久久久精品电影小说 | 直男gayav资源| 亚洲自偷自拍三级| 人人妻人人看人人澡| 午夜亚洲福利在线播放| 精品国产三级普通话版| 国产91av在线免费观看| 国产高清视频在线观看网站| 18+在线观看网站| av视频在线观看入口| 亚洲人成网站在线播放欧美日韩| 国产精品美女特级片免费视频播放器| 亚洲精品日韩av片在线观看| 国产综合懂色| 欧美性猛交黑人性爽| 国产成人freesex在线| 一边摸一边抽搐一进一小说| 精品人妻偷拍中文字幕| 国产v大片淫在线免费观看| 热99在线观看视频| 成人一区二区视频在线观看| 日韩成人av中文字幕在线观看| av福利片在线观看| 国产成人aa在线观看| 美女大奶头视频| 国产91av在线免费观看| 国产亚洲91精品色在线| 特级一级黄色大片| 亚洲精品国产av成人精品| 久久精品人妻少妇| 中文字幕免费在线视频6| 精品午夜福利在线看| 欧美日本视频| 天天躁日日操中文字幕| 国产爱豆传媒在线观看| 精品日产1卡2卡| 免费不卡的大黄色大毛片视频在线观看 | 日韩亚洲欧美综合| 欧美日韩在线观看h| 久久久久久久久久久丰满| 床上黄色一级片| 两性午夜刺激爽爽歪歪视频在线观看| 男插女下体视频免费在线播放| 日本免费一区二区三区高清不卡| 亚洲av不卡在线观看| 国产一区亚洲一区在线观看| 在线播放无遮挡| 六月丁香七月| 又粗又爽又猛毛片免费看| 亚洲av电影不卡..在线观看| 久久人人爽人人爽人人片va| 免费看光身美女| 国产亚洲av嫩草精品影院| 欧美又色又爽又黄视频| 久久精品国产鲁丝片午夜精品| 国产成人福利小说| 黄色一级大片看看| 亚洲国产精品成人久久小说 | 国产在线精品亚洲第一网站| 美女xxoo啪啪120秒动态图| 国产黄片美女视频| 51国产日韩欧美| 看非洲黑人一级黄片| 成人综合一区亚洲| а√天堂www在线а√下载| 亚洲第一区二区三区不卡| 夜夜夜夜夜久久久久| 国产蜜桃级精品一区二区三区| 国产伦精品一区二区三区四那| av视频在线观看入口| 国产精品一区二区性色av| 国产精品一二三区在线看| 国产伦精品一区二区三区视频9| 欧美最黄视频在线播放免费| 亚洲av熟女| 寂寞人妻少妇视频99o| 欧美丝袜亚洲另类| 久久这里只有精品中国| 小蜜桃在线观看免费完整版高清| 免费搜索国产男女视频| 男女啪啪激烈高潮av片| 亚洲久久久久久中文字幕| 一边亲一边摸免费视频| 色综合色国产| 亚洲人与动物交配视频| 99久久精品热视频| 免费av观看视频| 欧美成人a在线观看| 日韩一本色道免费dvd| 精品久久久久久成人av| 欧美一区二区国产精品久久精品| 人人妻人人澡人人爽人人夜夜 | 黄色视频,在线免费观看| 在线观看免费视频日本深夜| 国产日本99.免费观看| 亚洲av.av天堂| av在线亚洲专区| 国产精华一区二区三区| 成人三级黄色视频| 黄色欧美视频在线观看| 欧美激情在线99| 午夜福利高清视频| 日本在线视频免费播放| 中国国产av一级| 熟女电影av网| 男人狂女人下面高潮的视频| 99热只有精品国产| 亚洲最大成人中文| 亚洲成人久久爱视频| 搡女人真爽免费视频火全软件| 波多野结衣巨乳人妻| 国模一区二区三区四区视频| 亚洲真实伦在线观看| 丝袜喷水一区| 亚洲七黄色美女视频| 一进一出抽搐gif免费好疼| 亚洲久久久久久中文字幕| a级毛片免费高清观看在线播放| 国产精品1区2区在线观看.| 秋霞在线观看毛片| 免费电影在线观看免费观看| 伦精品一区二区三区| 麻豆国产av国片精品| 国内揄拍国产精品人妻在线| 中文字幕av在线有码专区| 日韩欧美精品免费久久| 深夜a级毛片| 一级毛片aaaaaa免费看小| 六月丁香七月| 久久99蜜桃精品久久| 久久久久久久亚洲中文字幕| 岛国在线免费视频观看| 欧美+亚洲+日韩+国产| 亚洲国产精品国产精品| 免费观看人在逋| 91久久精品电影网| 别揉我奶头 嗯啊视频| 一级二级三级毛片免费看| 寂寞人妻少妇视频99o| 亚洲人成网站在线播| 舔av片在线| 日本黄大片高清| 校园春色视频在线观看| 中文字幕熟女人妻在线| 国产午夜福利久久久久久| 中文字幕免费在线视频6| or卡值多少钱| 久久精品91蜜桃| 午夜精品国产一区二区电影 | 日韩一区二区视频免费看| 欧美成人一区二区免费高清观看| 久久99蜜桃精品久久| 国产精品电影一区二区三区| 中国美女看黄片| 亚洲精品日韩av片在线观看| 成人午夜精彩视频在线观看| 一本久久精品| 黄片wwwwww| 亚洲成av人片在线播放无| 在现免费观看毛片| 男女下面进入的视频免费午夜| 日韩在线高清观看一区二区三区| 午夜精品在线福利| 精品国产三级普通话版| 人妻制服诱惑在线中文字幕| 久久人人爽人人片av| 中文精品一卡2卡3卡4更新| av天堂中文字幕网| 国产精品伦人一区二区| 一级毛片电影观看 | 国产午夜精品久久久久久一区二区三区|