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

    面向石化企業(yè)間物流集成計(jì)劃優(yōu)化的模型及應(yīng)用

    2016-11-18 03:09:03張鵬飛王子豪榮岡馮毅萍
    化工學(xué)報(bào) 2016年11期
    關(guān)鍵詞:石腦油企業(yè)級(jí)原料

    張鵬飛,王子豪,榮岡,馮毅萍

    ?

    面向石化企業(yè)間物流集成計(jì)劃優(yōu)化的模型及應(yīng)用

    張鵬飛,王子豪,榮岡,馮毅萍

    (浙江大學(xué)智能系統(tǒng)與控制研究所,工業(yè)控制技術(shù)國家重點(diǎn)實(shí)驗(yàn)室,浙江杭州 310027)

    在石油化工行業(yè),上下游企業(yè)的聯(lián)合生產(chǎn)是供應(yīng)鏈優(yōu)化的主要途徑和重要手段。但由于石化企業(yè)的生產(chǎn)工藝過程復(fù)雜且設(shè)計(jì)相互獨(dú)立,企業(yè)間的協(xié)同生產(chǎn)往往僅存在于大型石化集團(tuán)統(tǒng)一設(shè)計(jì)、建造的煉化聯(lián)合生產(chǎn)系統(tǒng)中。本文提出了適用于獨(dú)立運(yùn)營(yíng)的石化企業(yè)集成的前提條件和建模方法,以獨(dú)立生產(chǎn)的煉油廠和烯烴廠為研究對(duì)象,基于其生產(chǎn)物流計(jì)劃模型,建立統(tǒng)一的聯(lián)合生產(chǎn)計(jì)劃優(yōu)化模型。該集成優(yōu)化模型為混合整數(shù)線性規(guī)劃模型,以生產(chǎn)效益最大化作為優(yōu)化求解目標(biāo),并分析了市場(chǎng)需求平穩(wěn)、萎縮和旺盛3種不同的案例場(chǎng)景。實(shí)例計(jì)算表明,同獨(dú)立生產(chǎn)相比較,石化企業(yè)的集成優(yōu)化和聯(lián)合生產(chǎn)能夠降低物料存儲(chǔ)成本和總成本,節(jié)省優(yōu)化求解所需的計(jì)算資源,并促進(jìn)經(jīng)濟(jì)效益的提升。

    石化企業(yè);集成;優(yōu)化;企業(yè)級(jí)優(yōu)化;混合整數(shù)規(guī)劃模型;生產(chǎn)效益

    引 言

    供應(yīng)鏈管理(supply chain management,SCM)對(duì)企業(yè)運(yùn)營(yíng)十分重要,其具體內(nèi)容包括生產(chǎn)過程的原料管理、物料運(yùn)輸管理,以及產(chǎn)品分配和運(yùn)輸?shù)纫幌盗猩a(chǎn)優(yōu)化命題[1]。圖1展示了供應(yīng)鏈管理的基本組成要素[2]。圖中,feedstocks代表系統(tǒng)的原料輸入部分,包括各種供應(yīng)鏈體系中所涉及的原始物料成分;basics代表供應(yīng)鏈中相互影響的上下游生產(chǎn)企業(yè),performance代表不同生產(chǎn)企業(yè)的生產(chǎn)過程與生產(chǎn)流水線,一條完整的供應(yīng)鏈包括一個(gè)及以上的企業(yè);供應(yīng)鏈的末端直接面向顧客,處理市場(chǎng)的訂單需求并調(diào)整產(chǎn)品策略,以及協(xié)調(diào)運(yùn)輸問題等。供應(yīng)鏈管理以滿足消費(fèi)者需求,提升企業(yè)效益并降低生產(chǎn)成本作為目標(biāo)[3-5]。

    而供應(yīng)鏈管理在具體的企業(yè)生產(chǎn)和執(zhí)行層面則是以企業(yè)作為主體,企業(yè)生產(chǎn)系統(tǒng)管理(process system engineering,PSE)是供應(yīng)鏈條的基本模塊。一個(gè)基本的生產(chǎn)系統(tǒng)管理模塊包括原料的輸入(raw materials),生產(chǎn)過程(performance),以及產(chǎn)品輸出(products)3個(gè)基本組分,其組成結(jié)構(gòu)如圖2所示[6]。

    近年來,隨著企業(yè)生產(chǎn)過程的復(fù)雜化,產(chǎn)品的多元化,以及上下游企業(yè)共享物料、設(shè)備和信息的需求增長(zhǎng),基于供應(yīng)鏈的生產(chǎn)系統(tǒng)的集成優(yōu)化成為企業(yè)提升生產(chǎn)效益重要途徑之一[1,7]。企業(yè)級(jí)優(yōu)化(enterprise-wide optimization,EWO)是企業(yè)集成生產(chǎn)優(yōu)化的重要形式之一,它將不同的生產(chǎn)系統(tǒng)管理模塊進(jìn)行集成,以提升企業(yè)的整體效益,綜合考慮集成企業(yè)的物料運(yùn)輸、生產(chǎn)方案、市場(chǎng)需求、庫存管理等問題[2,6]。Grossmann[8]在2005年就將企業(yè)級(jí)優(yōu)化思想引入流程工業(yè)領(lǐng)域。Wassick[2]以某一大型的石化集團(tuán)為例,將不同的生產(chǎn)子系統(tǒng)定義為“信封”,不同的“信封”可以組合成集成生產(chǎn)系統(tǒng),并對(duì)集成生產(chǎn)系統(tǒng)進(jìn)行優(yōu)化求解。Guyonnet等[9]將原油調(diào)度與產(chǎn)品問題進(jìn)行了集成分析。Quaglia等[10]指出企業(yè)級(jí)優(yōu)化的難點(diǎn)之一是集成系統(tǒng)的搭建與設(shè)計(jì),并以一個(gè)綜合的石化集團(tuán)為對(duì)象建立集成生產(chǎn)模型,綜合考慮集成系統(tǒng)的原料、產(chǎn)品和市場(chǎng)需求等要素,最后采用混合整數(shù)非線性規(guī)劃方法進(jìn)行建模求解。Shah等[11]研究了用于求解大規(guī)模企業(yè)級(jí)優(yōu)化模型的增廣拉格朗日分解求解策略。Quaglia等[12]分析了不確定條件下的企業(yè)級(jí)優(yōu)化命題,并提出了基于分解法的求解方法。Al-Qahtani等[13]綜合考慮了煉油石化系統(tǒng)中多個(gè)生產(chǎn)過程的集成生產(chǎn)過程,并通過協(xié)調(diào)上游生產(chǎn)和下游需求,提供了一套同步解析集成系統(tǒng)的方法策略。Wang等[14]以乙烯廠原料卸載調(diào)度和生產(chǎn)過程為研究對(duì)象,并提出了一種基于同步調(diào)度的優(yōu)化求解策略。Pinto等[15]提出了一種應(yīng)用于石化企業(yè)供應(yīng)鏈的企業(yè)級(jí)優(yōu)化框架,并給出了相應(yīng)的混合整數(shù)非線性數(shù)學(xué)模型。Shah等[16]以基于連續(xù)時(shí)間的集成優(yōu)化模型為對(duì)象,闡述了在企業(yè)級(jí)優(yōu)化命題中運(yùn)行時(shí)間、產(chǎn)量、質(zhì)量等多種生產(chǎn)要素之間的平衡。2014年,Grossmann[17]總結(jié)了近年來諸多學(xué)者在企業(yè)級(jí)優(yōu)化研究領(lǐng)域的研究與成果,評(píng)述了企業(yè)級(jí)優(yōu)化研究領(lǐng)域所涉及的求解算法,并指出了企業(yè)級(jí)優(yōu)化所面臨的重大挑戰(zhàn)在于研究對(duì)象的拓展、算法的改進(jìn)、不確定性干擾等。

    近年來諸多學(xué)者在企業(yè)級(jí)優(yōu)化的研究領(lǐng)域取得了大量成果,但由于不同的生產(chǎn)企業(yè)之間設(shè)計(jì)、生產(chǎn)和運(yùn)營(yíng)都相對(duì)獨(dú)立,這使得企業(yè)之間的系統(tǒng)集成面臨較大難度和挑戰(zhàn),現(xiàn)有的集成生產(chǎn)也往往局限于單一企業(yè)或者大規(guī)模生產(chǎn)系統(tǒng)的內(nèi)部[17-18],最終導(dǎo)致供應(yīng)鏈上下游企業(yè)在資源共享、信息共享的障礙以及企業(yè)生產(chǎn)效益提升受限。因此,本文將原本獨(dú)立的生產(chǎn)企業(yè)進(jìn)行整合,建立統(tǒng)一的集成生產(chǎn)系統(tǒng),而這一研究?jī)?nèi)容也是企業(yè)級(jí)優(yōu)化突破現(xiàn)有集成生產(chǎn)難點(diǎn)的重點(diǎn)方向之一[8]。

    本文以獨(dú)立設(shè)計(jì)建造、生產(chǎn)和經(jīng)營(yíng)的石化生產(chǎn)企業(yè)作為研究對(duì)象,以企業(yè)級(jí)優(yōu)化作為集成生產(chǎn)的方法,提出適用于獨(dú)立石化企業(yè)的集成條件、集成建模方法以及優(yōu)化模型,并根據(jù)集成對(duì)象的物流耦合關(guān)系重新定義產(chǎn)品、中間物料、原料,運(yùn)輸管道等信息集合及其元素組成。最后,以煉油廠和烯烴廠的集成生產(chǎn)為例建立MILP模型[19],所選用案例全面考慮了在市場(chǎng)需求平穩(wěn)、萎縮和旺盛3種市場(chǎng)情形下聯(lián)合生產(chǎn)集成優(yōu)化的效果。

    1 集成前提假設(shè)

    大多數(shù)獨(dú)立石化企業(yè)的生產(chǎn)系統(tǒng)往往是地理隔絕的,并且在生產(chǎn)時(shí)間尺度以及加工方案設(shè)計(jì)標(biāo)準(zhǔn)上不盡相同,從而導(dǎo)致不同企業(yè)之間的企業(yè)級(jí)優(yōu)化困難較大或者經(jīng)濟(jì)效益提升不明顯[15]。

    為了系統(tǒng)地評(píng)估獨(dú)立企業(yè)之間是否適合集成生產(chǎn),并篩選具有集成生產(chǎn)意義的生產(chǎn)系統(tǒng),本文提出了企業(yè)集成生產(chǎn)應(yīng)盡量滿足的前提條件。本文所研究的企業(yè)級(jí)優(yōu)化和集成生產(chǎn)也建立在以下前提條件之上。

    1.1 地理位置相近

    企業(yè)間物流運(yùn)輸成本將隨著相對(duì)距離的增加快速增長(zhǎng),較大的相對(duì)距離也對(duì)應(yīng)著較高的庫存成本和較慢的響應(yīng)。

    仿真計(jì)算表明,一般地,在廠際運(yùn)輸成本不超過普通廠內(nèi)運(yùn)輸成本的10倍時(shí),聯(lián)合生產(chǎn)集成優(yōu)化仍然能夠提升整體生產(chǎn)效益。

    1.2 物料與設(shè)備相同或相似

    相同的物料或加工設(shè)備,能在集成過程中通過物料、設(shè)備的共享,實(shí)現(xiàn)物料、資源的優(yōu)化配置,從而增強(qiáng)生產(chǎn)系統(tǒng)的魯棒性和穩(wěn)定性,最終促進(jìn)經(jīng)濟(jì)效益的總體提升。

    1.3 生產(chǎn)時(shí)間尺度相近或可調(diào)相近

    獨(dú)立生產(chǎn)系統(tǒng)的生產(chǎn)時(shí)間尺度相差過大,會(huì)導(dǎo)致企業(yè)級(jí)優(yōu)化過程中物料、設(shè)備共享受限,使得生產(chǎn)進(jìn)度不一致。因此,用于集成的獨(dú)立企業(yè)需要具有相近或者經(jīng)調(diào)整后相近的生產(chǎn)時(shí)間尺度[18]。

    1.4 生產(chǎn)利益統(tǒng)一

    企業(yè)級(jí)優(yōu)化構(gòu)建的集成生產(chǎn)系統(tǒng)在滿足一定條件時(shí),能夠?qū)崿F(xiàn)生產(chǎn)經(jīng)濟(jì)效益總和的提升,但這不代表任一原本獨(dú)立的生產(chǎn)系統(tǒng)經(jīng)濟(jì)效益均有所提升。因此,聯(lián)合生產(chǎn)的企業(yè)需要具有統(tǒng)一的利益或者合理的利益分配方式,才能保證各生產(chǎn)子系統(tǒng)能夠服務(wù)于集成系統(tǒng)的生產(chǎn)效益最大化目標(biāo)。

    基于以上前提,獨(dú)立生產(chǎn)的企業(yè)也能夠較好實(shí)現(xiàn)企業(yè)級(jí)優(yōu)化,并在集成生產(chǎn)中實(shí)現(xiàn)生產(chǎn)信息、資源的優(yōu)化配置。

    一般地,現(xiàn)有的石化企業(yè)建造往往具有集群效應(yīng),常見的石化企業(yè)生產(chǎn)園區(qū)或石化集團(tuán)下設(shè)的各類生產(chǎn)子廠均能夠很好地滿足以上假設(shè),從而保證了本文的假設(shè)具有實(shí)際物理意義和工程意義。本文案例中涉及的某真實(shí)石化集團(tuán)的生產(chǎn)系統(tǒng)均能夠滿足以上條件。

    2 企業(yè)級(jí)優(yōu)化方法與系統(tǒng)建模

    本節(jié)介紹獨(dú)立生產(chǎn)系統(tǒng)的集成方法以及集成系統(tǒng)的計(jì)劃優(yōu)化模型。

    2.1 集成生產(chǎn)系統(tǒng)構(gòu)建

    現(xiàn)有的研究中,較少提到集成系統(tǒng)的搭建方法和具體步驟,研究者往往通過分析集成生產(chǎn)系統(tǒng)的工藝要求進(jìn)行集成系統(tǒng)的搭建,而這些建模方法往往不通用也無法移植[2,14,20]。本文提出了一種集成系統(tǒng)的構(gòu)建方法,通過分析待集成企業(yè)之間的關(guān)鍵物料并搭建管道建立聯(lián)系,對(duì)同種生產(chǎn)裝置、存儲(chǔ)罐區(qū)的使用進(jìn)行調(diào)整,并重新定義原料、產(chǎn)品等生產(chǎn)信息集合3個(gè)步驟,能夠在不改變或者較少改變?cè)邢到y(tǒng)結(jié)構(gòu)的基礎(chǔ)上,快速構(gòu)建集成生產(chǎn)系統(tǒng)拓?fù)浣Y(jié)構(gòu)。這種集成系統(tǒng)構(gòu)建的方法具有簡(jiǎn)便、通用的優(yōu)點(diǎn)。具體的操作步驟如下。

    2.1.1 上下游企業(yè)的銜接 對(duì)于供應(yīng)鏈上下游的獨(dú)立企業(yè)而言,上游企業(yè)的產(chǎn)品或者半成品通??捎米飨掠纹髽I(yè)的原料。根據(jù)集成生產(chǎn)系統(tǒng)的物流耦合關(guān)系,將上游企業(yè)的終產(chǎn)品或半成品作為下游企業(yè)的生產(chǎn)原料;通過添加物流管道和物料緩存罐的方式,合理進(jìn)行上下游企業(yè)生產(chǎn)系統(tǒng)的拼接,組成集成的拓?fù)浣Y(jié)構(gòu)。

    2.1.2 相同物料與裝置的處理 在通過輸送管道將同種物料進(jìn)行關(guān)聯(lián)之后,便可以將相同的物料進(jìn)行整合,通過管道控制流向,使同種物料統(tǒng)一管理和分配。同時(shí),將相同的生產(chǎn)裝置進(jìn)行整合,通過裝置共用構(gòu)建多條加工線,或?qū)⒉糠滞N裝置用作備用裝置以增強(qiáng)系統(tǒng)的魯棒性。

    本文案例涉及的煉油廠、烯烴廠組成的集成生產(chǎn)系統(tǒng)中,石腦油、輕石腦油是煉油生產(chǎn)系統(tǒng)的重要產(chǎn)品,同時(shí)也是烯烴生產(chǎn)系統(tǒng)的重要原料。因此,集成系統(tǒng)中煉油廠的石腦油、輕石腦油產(chǎn)品可以作為烯烴廠的原料,即通過增加對(duì)應(yīng)的石腦油、輕石腦油輸送管線便可以將兩個(gè)生產(chǎn)系統(tǒng)進(jìn)行銜接。同時(shí),煉油廠、烯烴廠各自擁有的石腦油緩存罐可以實(shí)現(xiàn)共用,進(jìn)而提升石腦油緩存罐的利用率,并增加物料調(diào)度的靈活性。同理,輕石腦油的緩存罐也可共用。

    2.1.3 重新定義原料、產(chǎn)品等集合 集成系統(tǒng)的輸入輸出關(guān)系、物料的流動(dòng)關(guān)系在集成過程中發(fā)生了一定變化。所以,集成系統(tǒng)需要重新定義生產(chǎn)系統(tǒng)中原料、中間物料、產(chǎn)品以及運(yùn)輸管線等生產(chǎn)信息集合及其元素組成。

    本文案例涉及的煉油廠、烯烴廠集成生產(chǎn)系統(tǒng)中,石腦油、輕石腦油是煉油廠的產(chǎn)品,而在烯烴廠中則作為原料;在新建集成生產(chǎn)系統(tǒng)中,石腦油和輕石腦油則作為中間物料直接輸送給乙烯裝置進(jìn)行加工。上述生產(chǎn)流程的改變,要求重新定義系統(tǒng)原料、中間物料以及產(chǎn)品集合,具體地,在煉油廠系統(tǒng)中,石腦油、輕石腦油為產(chǎn)品集合的元素;烯烴生成系統(tǒng)中,石腦油、輕石腦油是原料集合的元素;而新建的集成系統(tǒng)中,則需要將石腦油、輕石腦油調(diào)整為中間物料集合的元素。

    通過以上主要步驟的操作,便可得到集成系統(tǒng)的拓?fù)浣Y(jié)構(gòu)及其基本生產(chǎn)信息。該集成系統(tǒng)保留了原有獨(dú)立生產(chǎn)企業(yè)的主要生產(chǎn)裝置和拓?fù)浣Y(jié)構(gòu),并能實(shí)現(xiàn)聯(lián)合生產(chǎn)過程中的物料、裝置、生產(chǎn)信息共享和統(tǒng)一調(diào)度。

    2.2 計(jì)劃優(yōu)化模型

    基于石化企業(yè)的生產(chǎn)拓?fù)浣Y(jié)構(gòu),本文建立了的基于MILP的計(jì)劃優(yōu)化模型。文中對(duì)生產(chǎn)拓?fù)浣Y(jié)構(gòu)進(jìn)行如下抽象:原料進(jìn)廠站點(diǎn)、產(chǎn)品出廠站點(diǎn)以及生產(chǎn)裝置為邏輯點(diǎn),不同的邏輯點(diǎn)之間由弧線連接代表著實(shí)際生產(chǎn)系統(tǒng)中的運(yùn)輸管道[20];同時(shí)對(duì)系統(tǒng)中的物料進(jìn)行以下劃分。

    原料集合:如原油、MTBE,BOC氧氣等,作為系統(tǒng)的輸入,供加工使用,被定義為系統(tǒng)的原料集合元素,其屬性包括名稱、市場(chǎng)價(jià)格等。

    中間物料集合:減壓瓦斯油、裂解碳四(C4)等在系統(tǒng)中產(chǎn)生,隨后被后級(jí)加工單元消耗,不作為系統(tǒng)的輸入和輸出物料,被劃分為中間物料。中間物料組分不考慮其物料價(jià)格。

    產(chǎn)品集合:汽油、柴油以及乙烯等,作為系統(tǒng)的輸出產(chǎn)品,直接向市場(chǎng)銷售,被定義為系統(tǒng)產(chǎn)品,其屬性包括名稱、市場(chǎng)價(jià)格等。

    2.2.1 優(yōu)化目標(biāo) 該計(jì)劃優(yōu)化模型以生產(chǎn)系統(tǒng)的生產(chǎn)效益最大化作為系統(tǒng)優(yōu)化目標(biāo),并同時(shí)兼顧生產(chǎn)、庫存成本的控制。

    式(1)中第1項(xiàng)為集成生產(chǎn)系統(tǒng)產(chǎn)品銷售的毛利潤(rùn),第2項(xiàng)為原料購買費(fèi)用,第3項(xiàng)為原料和產(chǎn)品的運(yùn)輸費(fèi)用,第4項(xiàng)為所有加工裝置產(chǎn)生的加工費(fèi)用,第5項(xiàng)與第6項(xiàng)為產(chǎn)品罐以及產(chǎn)品銷售站點(diǎn)的庫存費(fèi)用,第7項(xiàng)為市場(chǎng)產(chǎn)品需求未滿足的懲罰項(xiàng)。

    2.2.2 模型約束條件 本節(jié)給出集成生產(chǎn)系統(tǒng)計(jì)劃優(yōu)化模型中各種生產(chǎn)約束條件。

    (3)

    (4)

    (6)

    式(2)定義每一種輸入集成生產(chǎn)系統(tǒng)的原料購買量;式(3)限定了原料購買量的取值范圍。式(4)表示產(chǎn)品產(chǎn)量不能超過市場(chǎng)的實(shí)際需求量。式(5)與式(6)指出當(dāng)產(chǎn)品供給不能滿足預(yù)設(shè)需求時(shí),供給差額等于預(yù)設(shè)的產(chǎn)品需求量減去實(shí)際產(chǎn)品輸出量,且供給差額為非負(fù)值。

    (8)

    (10)

    (11)

    (12)

    式(7)與式(8)定義生產(chǎn)系統(tǒng)物料輸入輸出的流量同與之相連的側(cè)線流量一致。式(9)與式(10)定義生產(chǎn)系統(tǒng)內(nèi)部輸送管線的流量同與之相連的側(cè)線流量相等。式(11)與式(12)分別限定了生產(chǎn)系統(tǒng)輸入輸出管道和內(nèi)部輸送管道流量值的范圍。

    (13)

    式(13)表示在系統(tǒng)管道分流點(diǎn)處,物料的進(jìn)料量等于各條管線該物料出量的總和。式(14)表示在系統(tǒng)管道匯流點(diǎn)處,同種物料各管線進(jìn)料總和等于該物料的總出料量。

    (16)

    (17)

    (18)

    式(15)~式(18)對(duì)生產(chǎn)系統(tǒng)中加工裝置的生產(chǎn)方案進(jìn)行了描述與限制。式(15)限定了任意裝置的特定加工方案處理的物料量的上下限。式(16)表示任意加工裝置的加工總量等于該生產(chǎn)裝置各種加工方案加工量之和。式(17)表示任意加工裝置的加工總量還等于對(duì)應(yīng)生產(chǎn)周期內(nèi)該裝置進(jìn)料量的總和。式(18)描述加工裝置的產(chǎn)率模型,定義在不同產(chǎn)率模型下加工裝置各產(chǎn)出物料的比例。

    計(jì)劃優(yōu)化模型中站點(diǎn)與罐區(qū)具有物料存儲(chǔ)功能,為了區(qū)分站點(diǎn)與罐區(qū)的初始物料量,將生產(chǎn)的第一個(gè)周期定義為初始生產(chǎn)周期0。

    (19)

    (21)

    (23)

    (24)

    式(19)與式(20)定義在初始生產(chǎn)周期0時(shí),站點(diǎn)存儲(chǔ)的物料量等于站點(diǎn)初始存儲(chǔ)量加上站點(diǎn)進(jìn)側(cè)線流量再減去站點(diǎn)出側(cè)線流量;其余生產(chǎn)周期內(nèi),站點(diǎn)存儲(chǔ)的物料量等于上一生產(chǎn)周期末物料存儲(chǔ)量加上本生產(chǎn)周期進(jìn)側(cè)線流量再減去本生產(chǎn)周期出側(cè)線流量。式(21)限定了站點(diǎn)存儲(chǔ)量的上下限。

    同理,式(22)與式(23)定義在初始生產(chǎn)周期0時(shí),罐區(qū)存儲(chǔ)的物料量等于罐區(qū)初始存儲(chǔ)量加上罐區(qū)進(jìn)側(cè)線流量再減去罐區(qū)出側(cè)線流量;其余生產(chǎn)周期內(nèi),罐區(qū)存儲(chǔ)的物料量等于上一生產(chǎn)周期末物料存儲(chǔ)量加上本生產(chǎn)周期進(jìn)側(cè)線流量再減去本生產(chǎn)周期出側(cè)線流量。式(24)限定了罐區(qū)存儲(chǔ)量的上下限。

    3 應(yīng)用案例

    3.1 案例描述

    本文以某石化煉油廠、烯烴廠作為研究對(duì)象,建立具有代表性的石化企業(yè)生產(chǎn)拓?fù)浣Y(jié)構(gòu)。

    煉油生產(chǎn)系統(tǒng)的生產(chǎn)拓?fù)浣Y(jié)構(gòu)如圖3所示,其中包含常減壓裝置,延遲焦化裝置,催化裂化裝置等主要生產(chǎn)裝置;并以原油、甲基叔丁基醚(methyl-butyl ether,MTBE)作為原料,以柴油、90#汽油、93#汽油以及煤油作為產(chǎn)品。

    烯烴生產(chǎn)的生產(chǎn)拓?fù)浣Y(jié)構(gòu)如圖4所示,其中包含乙烯裝置、乙二醇裝置、丁二烯裝置等主要生產(chǎn)裝置;并以石腦油、丙烯、輕石腦油、裂化尾油、減壓瓦斯油(vacuum gas oil, VGO)和叔丁氧羰基(-butyl oxy carbonyl,BOC)氧氣作為原料,以乙烯、乙二醇、二乙二醇、環(huán)氧乙烷、丁二烯等作為產(chǎn)品。

    由于石腦油是煉油廠最為重要的產(chǎn)品并同時(shí)是烯烴廠最主要的原料之一,本文根據(jù)2.1節(jié)中的集成系統(tǒng)構(gòu)造方法,將石腦油和輕石腦油作為二廠銜接的關(guān)鍵物料,將煉油廠和烯烴廠的生產(chǎn)系統(tǒng)進(jìn)行集成,構(gòu)建如圖5所示的集成生產(chǎn)系統(tǒng)。集成系統(tǒng)中包含一次加工裝置:常減壓裝置兩套;二次加工裝置:催化重整裝置,延遲焦化裝置,煤油加氫裝置,柴油加氫裝置,催化裂化裝置兩套,乙烯裝置兩套,乙二醇裝置,丁二烯裝置兩套,汽油加氫裝置兩套;調(diào)和裝置:柴油調(diào)和裝置以及汽油調(diào)和裝置;系統(tǒng)中總計(jì)12種17套不同的生產(chǎn)加工裝置。同時(shí),集成系統(tǒng)還包含原油緩存罐、柴油混合罐、FCC進(jìn)料緩存罐、柴油產(chǎn)品罐、乙烯原料緩存罐等共計(jì)35個(gè)存儲(chǔ)罐。系統(tǒng)的產(chǎn)品包括:90#汽油、93#汽油、柴油、煤油、乙烯、乙二醇、二乙二醇、三乙二醇、環(huán)氧乙烷、丁二烯等。

    本案例中的生產(chǎn)時(shí)間周期為半年總計(jì)6個(gè)月,每個(gè)月以30天計(jì)算。同時(shí)每一天分為兩個(gè)操作批次,即每月具有60個(gè)操作批次。生產(chǎn)周期信息如表1所示。

    表1 生產(chǎn)周期

    參考市場(chǎng)價(jià)格,對(duì)生產(chǎn)系統(tǒng)中每一種物料的價(jià)格進(jìn)行設(shè)定。如表2所示,其中各物料的價(jià)格高低體現(xiàn)了其市場(chǎng)價(jià)值的相對(duì)大小。

    表2 物料價(jià)格

    在實(shí)際生產(chǎn)中,廠際運(yùn)輸?shù)某杀鞠鄬?duì)廠內(nèi)運(yùn)輸增長(zhǎng)較大。本案例中,在滿足第1節(jié)中距離較近的前提條件下,設(shè)定廠際運(yùn)輸?shù)某杀緸閺S內(nèi)運(yùn)輸成本的10倍。

    3.2 求解策略

    本文以生產(chǎn)經(jīng)濟(jì)效益作為主要指標(biāo)對(duì)比集成系統(tǒng)和獨(dú)立系統(tǒng)的生產(chǎn)情況,所建立的計(jì)劃優(yōu)化模型由市場(chǎng)需求(訂單)驅(qū)動(dòng)。

    計(jì)算過程中,首先設(shè)定集成生產(chǎn)系統(tǒng)的市場(chǎng)需求,計(jì)算集成生產(chǎn)系統(tǒng)的生產(chǎn)效益最優(yōu)解。然后,基于集成生產(chǎn)的生產(chǎn)信息,設(shè)定獨(dú)立運(yùn)行的煉油廠、烯烴廠二者獨(dú)立運(yùn)行時(shí)的市場(chǎng)需求,保證優(yōu)化求解得到的生產(chǎn)效益同集成系統(tǒng)具有可對(duì)比性。最后再計(jì)算二廠各自的生產(chǎn)效益最優(yōu)解,將獨(dú)立生產(chǎn)和集成生產(chǎn)的效益進(jìn)行對(duì)比分析。圖6展示了求解分析策略和思路。

    3.3 模型求解與分析

    根據(jù)2.2節(jié)所建立的計(jì)劃優(yōu)化MILP模型,分別對(duì)煉油廠、烯烴廠以及集成生產(chǎn)系統(tǒng)進(jìn)行優(yōu)化求解。案例求解采用IBM ILOG Optimization Studio 12.4作為優(yōu)化求解平臺(tái),所調(diào)用的求解器為:CPLEX 12.4.0。求解平臺(tái)所依托的PC配置為:Intel Xeon E5-2403 CPU,主頻1.8GHz,內(nèi)存8G,操作系統(tǒng)Windows Sever 2008。

    如表3所示,集成生產(chǎn)系統(tǒng)的優(yōu)化模型總共包含6408個(gè)約束條件,整數(shù)變量為186個(gè),連續(xù)變量為2875個(gè);煉油生產(chǎn)系統(tǒng)的優(yōu)化模型包括總計(jì)2916條約束條件,整數(shù)變量102個(gè),連續(xù)變量1297個(gè);烯烴生產(chǎn)系統(tǒng)的優(yōu)化模型包含總計(jì)3198條約束條件,整數(shù)變量54個(gè),連續(xù)變量1411個(gè)。

    表3 優(yōu)化模型參數(shù)對(duì)比

    下文考慮了在市場(chǎng)需求平穩(wěn)、萎縮以及旺盛3種市場(chǎng)情況下,集成生產(chǎn)、獨(dú)立生產(chǎn)對(duì)于企業(yè)生產(chǎn)效益的影響。

    3.3.1 案例一:市場(chǎng)需求平穩(wěn) 當(dāng)市場(chǎng)需求平穩(wěn)時(shí),石化生產(chǎn)系統(tǒng)處于平穩(wěn)生產(chǎn)的狀態(tài)。生產(chǎn)系統(tǒng)的產(chǎn)能處于正常水平,產(chǎn)品產(chǎn)量能夠滿足大部分的市場(chǎng)需求。

    生產(chǎn)系統(tǒng)初始條件的設(shè)定原則為:產(chǎn)品需求參考市場(chǎng)需求平穩(wěn)時(shí)的產(chǎn)品產(chǎn)量設(shè)定,具有存儲(chǔ)功能的罐區(qū)和站點(diǎn)初始庫存略高于安全值的底線并擁有較大的緩存區(qū)間。

    分別求解集成系統(tǒng),煉油生產(chǎn)系統(tǒng)和烯烴生產(chǎn)系統(tǒng),得到的優(yōu)化求解結(jié)果如表4所示。

    表4 案例一優(yōu)化結(jié)果對(duì)比

    如表4所示,系統(tǒng)集成能夠帶來的效益提升較小。獨(dú)立系統(tǒng)生產(chǎn)的利潤(rùn)分別為14042542元和9786211元,集成系統(tǒng)協(xié)同生產(chǎn)的利潤(rùn)為25464324元,生產(chǎn)效益總體提升約6.8%。獨(dú)立生產(chǎn)的總收入與總成本均高于集成生產(chǎn)系統(tǒng),但集成生產(chǎn)系統(tǒng)的最終利潤(rùn)卻更高。

    當(dāng)市場(chǎng)需求穩(wěn)定時(shí),系統(tǒng)生產(chǎn)平穩(wěn),互為上下游企業(yè)的煉油廠和乙烯廠能夠較好地激發(fā)產(chǎn)能,獲得不錯(cuò)的生產(chǎn)效益。在獨(dú)立生產(chǎn)模式下,獨(dú)立生產(chǎn)的系統(tǒng)均以自身效益最大化作為生產(chǎn)目標(biāo),銷售更多產(chǎn)品,也消耗更多的原料,導(dǎo)致原料成本增加;同時(shí),獨(dú)立生產(chǎn)的煉油廠和乙烯廠匹配市場(chǎng)需求的能力較弱,部分產(chǎn)品庫存增加,導(dǎo)致庫存成本增長(zhǎng),而另一部分產(chǎn)品不能滿足市場(chǎng)需求,導(dǎo)致需求未滿足懲罰值增加。而另一方面,集成生產(chǎn)能夠通過調(diào)整產(chǎn)品方案組合,更好地滿足市場(chǎng)需求,控制庫存成本;同時(shí)集成系統(tǒng)內(nèi)部物料可以作為自身的原料,從而節(jié)省了原料的購買、運(yùn)輸成本;最終集成生產(chǎn)能夠?qū)崿F(xiàn)生產(chǎn)效益的提升。

    總體上,集成生產(chǎn)由于較好地控制了庫存成本和原料購買、運(yùn)輸成本,實(shí)現(xiàn)了生產(chǎn)方案的優(yōu)化配置,進(jìn)而獲得了更好的效益。

    3.3.2 案例二:市場(chǎng)需求萎縮 當(dāng)市場(chǎng)經(jīng)濟(jì)蕭條時(shí),各種化工產(chǎn)品的市場(chǎng)需求也進(jìn)入蕭條期。此時(shí)市場(chǎng)極度萎縮,產(chǎn)品的市場(chǎng)需求量大大低于平穩(wěn)時(shí)期甚至銳減為零,生產(chǎn)系統(tǒng)的實(shí)際產(chǎn)能下降。

    生產(chǎn)系統(tǒng)初始條件的設(shè)定原則為:產(chǎn)品需求減少甚至降低為零,設(shè)定的產(chǎn)品需求低于市場(chǎng)平穩(wěn)時(shí)期;緩存罐的初始儲(chǔ)量高于正常水平甚至接近于存儲(chǔ)安全上限值。

    分別求解集成系統(tǒng),煉油生產(chǎn)系統(tǒng)和烯烴生產(chǎn)系統(tǒng),得到的優(yōu)化求解結(jié)果如表5所示。

    表5 案例二優(yōu)化結(jié)果對(duì)比

    如表5所示,獨(dú)立生產(chǎn)系統(tǒng)的獨(dú)立生產(chǎn)得到的利益分別是:6324364和8357961元,而集成生產(chǎn)系統(tǒng)的生產(chǎn)效益為18803812元,生產(chǎn)效益相較于獨(dú)立生產(chǎn)提升了28.1%。

    當(dāng)市場(chǎng)萎縮時(shí),生產(chǎn)的產(chǎn)品被迫滯銷,系統(tǒng)產(chǎn)能受到抑制,產(chǎn)品積壓導(dǎo)致存儲(chǔ)罐區(qū)和站點(diǎn)初始庫存增大,獨(dú)立生產(chǎn)和集成生產(chǎn)的生產(chǎn)效益相較于市場(chǎng)平穩(wěn)時(shí)期都急劇減少。

    集成生產(chǎn)通過調(diào)整生產(chǎn)方案和產(chǎn)品組合,能相對(duì)更好地滿足市場(chǎng)產(chǎn)品需求并出售產(chǎn)品,從而在一定程度控制庫存成本;而另一方面,獨(dú)立生產(chǎn)造成了極端的兩極現(xiàn)象,即部分產(chǎn)品脫銷而部分產(chǎn)品滯銷,造成了高庫存成本和高懲罰值。因此,在市場(chǎng)萎縮的情況下,集成生產(chǎn)能夠更好地控制庫存,調(diào)整產(chǎn)品方案,從而相對(duì)于獨(dú)立生產(chǎn)有著較高的效益提升;這也說明了在市場(chǎng)萎縮時(shí),集成生產(chǎn)能夠極大限度地減少企業(yè)損失,幫助企業(yè)渡過難關(guān)。

    3.3.3 案例三:市場(chǎng)需求旺盛 當(dāng)市場(chǎng)需求旺盛時(shí),各種產(chǎn)品的市場(chǎng)需求普遍高于市場(chǎng)平穩(wěn)時(shí)期。生產(chǎn)系統(tǒng)的產(chǎn)能被釋放,能夠生產(chǎn)和銷售更多的產(chǎn)品。

    (1)情況一:產(chǎn)品價(jià)格不變或同倍率上浮

    正常情況下,供應(yīng)鏈后級(jí)的產(chǎn)品具有更高的成本和附加價(jià)值,其價(jià)格不會(huì)低于前級(jí)產(chǎn)品。在市場(chǎng)需求旺盛時(shí),產(chǎn)品的價(jià)格會(huì)同倍率上浮或受政策控制維持不變。

    生產(chǎn)系統(tǒng)初始條件的設(shè)定原則為:產(chǎn)品需求增加且高于平時(shí)的需求水平,設(shè)定的產(chǎn)品需求高于市場(chǎng)平穩(wěn)時(shí)的產(chǎn)品需求值,產(chǎn)品的價(jià)格同市場(chǎng)平穩(wěn)時(shí)持平;緩存罐的初始儲(chǔ)量低于正常水平甚至接近于安全下限值。

    分別求解集成系統(tǒng),煉油生產(chǎn)系統(tǒng)和烯烴生產(chǎn)系統(tǒng),得到的優(yōu)化求解結(jié)果如表6所示。

    表6 情況一優(yōu)化結(jié)果對(duì)比

    如表6所示,獨(dú)立生產(chǎn)系統(tǒng)的獨(dú)立生產(chǎn)得到的利益分別是:19572897和14299543元;而集成生產(chǎn)系統(tǒng)的生產(chǎn)效益為34461395元,相較于獨(dú)立生產(chǎn)的總和提升1.7%。

    當(dāng)市場(chǎng)需求旺盛時(shí),系統(tǒng)的產(chǎn)能得到最大的釋放,生產(chǎn)裝置滿負(fù)荷運(yùn)行,保證產(chǎn)品的收益總和最大化,此時(shí)產(chǎn)品產(chǎn)量的主要限制是系統(tǒng)本身的物理特質(zhì),如生產(chǎn)系統(tǒng)的設(shè)計(jì)產(chǎn)能上限值等。相較于獨(dú)立生產(chǎn),集成生產(chǎn)的原料成本下降幅度高于銷售收入的減少幅度,因此集成生產(chǎn)由于提高了物料利用率和使用效率,實(shí)現(xiàn)了集成生產(chǎn)的效益提升。

    在市場(chǎng)需求旺盛的情況下,煉油廠作為上游企業(yè),其產(chǎn)品可以銷售給除烯烴廠以外的其他企業(yè)或者個(gè)人,以實(shí)現(xiàn)更大的經(jīng)濟(jì)效益。此時(shí),集成生產(chǎn)需建立在基于1.4節(jié)提出“利益統(tǒng)一”的前提上,只有保證煉油廠和烯烴廠具有統(tǒng)一的利益或合理的利益分配方式時(shí),才能保證煉油廠首先滿足烯烴廠的原料供給,以實(shí)現(xiàn)集成生產(chǎn)。從最終結(jié)果來看,集成生產(chǎn)的生產(chǎn)效益高于獨(dú)立生產(chǎn)效益總和的,即集成生產(chǎn)依舊能夠保障各生產(chǎn)子系統(tǒng)的利益并實(shí)現(xiàn)效益提升。

    (2)情況二:價(jià)格異常波動(dòng)

    極少數(shù)極端情況下,部分前級(jí)產(chǎn)品的價(jià)格或者價(jià)格上浮幅度會(huì)反常提升,從而較大程度影響供應(yīng)鏈下游企業(yè)的生產(chǎn)情況。為了對(duì)此極端情況進(jìn)行說明,在情況一的基本設(shè)定上,設(shè)定煉油廠石腦油的價(jià)格上浮50%,其他條件保持不變。相應(yīng)的求解結(jié)果如表7所示。

    表7 情況二優(yōu)化結(jié)果對(duì)比

    如表7所示,獨(dú)立生產(chǎn)系統(tǒng)的獨(dú)立生產(chǎn)得到的利益分別是:38544911和532784元;而集成生產(chǎn)系統(tǒng)的生產(chǎn)效益為38591135元,生產(chǎn)效益相較于獨(dú)立生產(chǎn)的總和下降1.2%。

    反常情況下,當(dāng)上游企業(yè)的產(chǎn)品價(jià)格持續(xù)上漲到臨界點(diǎn)時(shí)(本案例中為石腦油價(jià)格上浮50%左右),下游企業(yè)會(huì)因?yàn)樯a(chǎn)成本的提升而導(dǎo)致整體效益的下降。本案例中,由于石腦油價(jià)格反常上升,導(dǎo)致烯烴廠的原料成本急劇上升,而銷售的收入沒有增加,從而導(dǎo)致利潤(rùn)減少。此時(shí),下游企業(yè)可能會(huì)通過選擇停車的方式應(yīng)對(duì)價(jià)格反常波動(dòng)。但總體上,此類極端情況是十分罕見的。

    3.4 小結(jié)

    在市場(chǎng)需求平穩(wěn)、萎縮和旺盛3種情況下,獨(dú)立生產(chǎn)和集成生產(chǎn)系統(tǒng)生產(chǎn)效益和優(yōu)化參數(shù)的對(duì)比分析表明:集成生產(chǎn)能夠更好地調(diào)整產(chǎn)品生產(chǎn)方案的組合以適應(yīng)產(chǎn)品需求的波動(dòng),更好地控制產(chǎn)品庫存量降低庫存成本,實(shí)現(xiàn)物料、生產(chǎn)裝置的優(yōu)化配置,從而實(shí)現(xiàn)相較于獨(dú)立生產(chǎn)系統(tǒng)的生產(chǎn)效益提升。特別地,在市場(chǎng)萎縮的情況下,集成生產(chǎn)系統(tǒng)能夠極大程度地減少企業(yè)損失,幫助企業(yè)控制庫存成本并進(jìn)行生產(chǎn)方案的調(diào)整,相較于獨(dú)立生產(chǎn)更有利于幫助企業(yè)渡過經(jīng)濟(jì)蕭條期;在市場(chǎng)平穩(wěn)和旺盛的情況下,集成生產(chǎn)也能夠通過裝置共用,控制庫存和調(diào)整生產(chǎn)方案等方式促進(jìn)生產(chǎn)總效益的提升。

    4 結(jié) 論

    本文討論了石化系統(tǒng)企業(yè)級(jí)優(yōu)化的命題,提出了獨(dú)立生產(chǎn)的石化企業(yè)集成生產(chǎn)的前提條件,保證了獨(dú)立企業(yè)集成生產(chǎn)的有效性;同時(shí),提出了生產(chǎn)系統(tǒng)的集成方法,并建立集成生產(chǎn)系統(tǒng)的計(jì)劃優(yōu)化數(shù)學(xué)模型。文章以某石化煉油廠、烯烴廠作為案例對(duì)象,研究了企業(yè)獨(dú)立生產(chǎn)和集成生產(chǎn)在市場(chǎng)需求平穩(wěn)、需求萎縮和需求旺盛3種典型情況下的生產(chǎn)效益情況。案例實(shí)踐表明,基于企業(yè)級(jí)優(yōu)化的集成生產(chǎn)在滿足一定前提條件時(shí),具有集成優(yōu)化和聯(lián)合生產(chǎn)的商業(yè)價(jià)值;集成優(yōu)化能夠?qū)崿F(xiàn)生產(chǎn)裝置、緩存罐區(qū)站點(diǎn)、生產(chǎn)物料的共用和優(yōu)化配置,控制原料成本與庫存成本,促進(jìn)更優(yōu)的生產(chǎn)方案配置,最終實(shí)現(xiàn)生產(chǎn)效益的提升。

    符 號(hào) 說 明

    C——產(chǎn)品輸出站點(diǎn)集合 CDP——產(chǎn)品需求未滿足懲罰費(fèi)用,元 CHO——庫存管理費(fèi)用,元 COP——裝置加工費(fèi)用,元 CPRO——產(chǎn)品銷售費(fèi)用,元 CPUR——原料購買費(fèi)用,元 CTS——物料運(yùn)輸費(fèi)用,元 DEMcm,t——物料cm在生產(chǎn)周期t的需求量,噸 DPCcm——產(chǎn)品cm在未滿足需求時(shí)單位物料量的懲罰費(fèi)用 DPQcm,t——產(chǎn)品cm在時(shí)間t內(nèi)的未滿足量,噸 DQcm,t——產(chǎn)品cm在時(shí)間t內(nèi)的需求量,噸 HCu,cm, t——裝置u對(duì)物料cm在時(shí)間t內(nèi)的庫存費(fèi)用,元 INum,cm——加工方案um中對(duì)應(yīng)物料cm的進(jìn)料比例 MOu——裝置u的所有加工方案集合 MOIu——裝置u加工方案信息集合 MOQu,um,t——裝置u處于方案um時(shí)在時(shí)間t內(nèi)的加工量,噸 MQout,in,mv,cm,t——從出側(cè)線out到進(jìn)側(cè)線in的移動(dòng)mv在時(shí)間t內(nèi)的物料cm的流量,噸 MT——所有物料集合 MTc——所有產(chǎn)品物料集合 MTR——所有原料物料集合 MTs——所有站點(diǎn)s物料集合 MTss——存儲(chǔ)站點(diǎn)SS物料集合 MTtk——罐區(qū)TK物料集合 MV——生產(chǎn)系統(tǒng)內(nèi)部移動(dòng)運(yùn)輸集合 OBJ——優(yōu)化目標(biāo),即生產(chǎn)經(jīng)濟(jì)效益,元 OCu——加工裝置u的加工單價(jià),元 Pcm——物料cm的單價(jià),元 R——原料供應(yīng)點(diǎn)(供應(yīng)商)集合 RQcm,t——原料rm在時(shí)間t內(nèi)的采購量,噸 RUNu,um,t——裝置u在加工方案m下于時(shí)間t內(nèi)的加工次數(shù),次 S——進(jìn)出廠站點(diǎn)集合 SM——分匯流點(diǎn)集合 SMm——匯流點(diǎn)集合 SMs——分流點(diǎn)集合 SQu,cm,t——存儲(chǔ)站點(diǎn)u在時(shí)間t內(nèi)物料cm的存儲(chǔ)量,噸 SS——存儲(chǔ)站點(diǎn)集合 ST——側(cè)線集合 STi——進(jìn)側(cè)線集合 STo——出側(cè)線集合 STu——生產(chǎn)裝置u的側(cè)線集合 STRIQno,in,cm,t——節(jié)點(diǎn)no(裝置,站點(diǎn)或者分匯流點(diǎn))的進(jìn)側(cè)線st在時(shí)間t內(nèi),進(jìn)側(cè)線上的物料cm的流量,噸 STROQno,st,cm,t——節(jié)點(diǎn)no(裝置,站點(diǎn)或者分匯流點(diǎn))的出側(cè)線st在時(shí)間t內(nèi),出側(cè)線上物料cm的流量,噸 STRQno,st,cm,t——節(jié)點(diǎn)no(裝置,站點(diǎn)或者分匯流點(diǎn))的側(cè)線st上在時(shí)間t內(nèi),側(cè)線上物料cm的流量,噸 T——生產(chǎn)周期時(shí)間集合 TCts——運(yùn)輸移動(dòng)ts的運(yùn)送單價(jià),元 TK——存儲(chǔ)罐區(qū)集合 TQu,cm,t——存儲(chǔ)罐區(qū)u在時(shí)間t內(nèi)物料cm的存儲(chǔ)量,噸 TR——生產(chǎn)系統(tǒng)邊界移動(dòng)運(yùn)輸集合 TRQout,in,ts,cm,t——從出側(cè)線out到進(jìn)側(cè)線in的移動(dòng)運(yùn)輸ts在時(shí)間t內(nèi)的物料cm的流量,噸 U——加工裝置集合 UQu,t——裝置u在時(shí)間t內(nèi)的加工量,噸 YIEum,cm——加工方案um中對(duì)應(yīng)物料cm的產(chǎn)率 上角標(biāo) ini——初始值 lb——下限 up——上限 下角標(biāo) cm——物料 in——進(jìn)側(cè)線 mv——移動(dòng) no——節(jié)點(diǎn) out ——出側(cè)線 rm——原料 s——站點(diǎn) st——側(cè)線 t——時(shí)間周期 tk——存儲(chǔ)罐 ts——運(yùn)輸 u——裝置 um——裝置加工方案

    References

    [1] MARAVELIAS C T, SUNG C. Integration of production planning and scheduling: overview, challenges and opportunities[J]. Computers & Chemical Engineering, 2009, 33(12): 1919-1930.

    [2] WASSICK J M. Enterprise-wide optimization in an integrated chemical complex[J]. Computers & Chemical Engineering, 2009, 33(12): 1950-1963.

    [3] ISLAM D M Z, MEIER J F, ADITJANDRA P T,. Logistics and supply chain management[J]. Research in Transportation Economics, 2012, 41(1): 1-14.

    [4] BEAMON B M. Supply chain design and analysis: models and methods[J]. International Journal of Production Economics, 1998, 55(3): 281-294.

    [5] MIN H, ZHOU G. Supply chain modeling: past, present and future[J]. Computers & Industrial Engineering, 2002, 43(1): 231-249.

    [6] BOK J K, GROSSMANN I E, PARK S. Supply chain optimization in continuous flexible process networks[J]. Industrial & Engineering Chemistry Research, 2000, 39(5): 1279-1290.

    [7] ERENGü? ? S, SIMPSON N C, VAKHARIA A J. Integrated production/distribution planning in supply chains: an invited review[J]. European Journal of Operational Research, 1999, 115(2): 219-236.

    [8] GROSSMANN I. Enterprise-wide optimization: a new frontier in process systems engineering[J]. AIChE Journal, 2005, 51(7): 1846-1857.

    [9] GUYONNET P, GRANT F H, BAGAJEWICZ M J. Integrated model for refinery planning, oil procuring, and product distribution[J]. Industrial & Engineering Chemistry Research, 2008, 48(1): 463-482.

    [10] QUAGLIA A, SARUP B, SIN G,. Integrated business and engineering framework for synthesis and design of enterprise-wide processing networks[J]. Computers & Chemical Engineering, 2012, 38(10): 213-223.

    [11] SHAH N K, IERAPETRITOU M G. Integrated production planning and scheduling optimization of multi-site, multi-product process industry[J]. Computers & Chemical Engineering, 2011, 29(2): 214-226.

    [12] QUAGLIA A, SARUP B, GüRKAN S,. A systematic framework for enterprise-wide optimization: synthesis and design of processing networks under uncertainty[J]. Computers & Chemical Engineering, 2013, 59(10): 47-62.

    [13] AL-QAHTANI K, ELKAMEL A. Multisite refinery and petrochemical network design: optimal integration and coordination[J]. Industrial & Engineering Chemistry Research, 2008, 48(2): 814-826.

    [14] WANG Z, FENG Y, RONG G. Synchronized scheduling approach of ethylene plant production and naphtha oil inventory management[J]. Industrial & Engineering Chemistry Research, 2014, 53(15): 6477-6499.

    [15] NEIRO S M S, PINTO J M. A general modeling framework for the operational planning of petroleum supply chains[J]. Computers & Chemical Engineering, 2004, 28(6): 871-896.

    [16] SHAH N K, IERAPETRITOU M G. Short-term scheduling of a large-scale oil-refinery operations: incorporating logistics details[J]. AIChE Journal, 2011, 57(6): 1570-1584.

    [17] GROSSMANN I E. Challenges in the application of mathematical programming in the enterprise-wide optimization of process industries[J]. Theoretical Foundations of Chemical Engineering, 2014, 48(5): 555-573.

    [18] VARMA V A, REKLAITIS G V, BLAU G E,. Enterprise-wide modeling & optimization-an overview of emerging research challenges and opportunities[J]. Computers & Chemical Engineering, 2007, 31(5/6): 692-711.

    [19] FLOUDAS C A, LIN X. Mixed integer linear programming in process scheduling: modeling, algorithms, and applications[J]. Annals of Operations Research, 2005, 139(1): 131-162.

    [20] GUBITOSO F, PINTO J M. A planning model for the optimal production of a real-world ethylene plant[J]. Chemical Engineering & Processing Process Intensification, 2007, 46(11):1141-1150.

    Optimized model and application of integrated logistics planning for petrochemical enterprises

    ZHANG Pengfei, WANG Zihao, RONG Gang, FENG Yiping

    (State Key Laboratory of Industrial Control Technology, Institute of Cyber-System and Control, Zhejiang University, Hangzhou 310027, Zhejiang, China)

    The synergetic production among upstream and downstream enterprises is a critical step toward supply chain optimization across petrochemical industry. However, because of complexity and independent design of petrochemical process, such a synergetic operation between facilities is available only in a few conglomerate corporations with coordinated manufacturing systems. An optimized model for integrated logistic planning in independently operating petrochemical companies was proposed on a basis of current supply chain model of independent refinery and olefin plants. In this mixed-integer linear programming (MILP) model, maximum operation profit was targeted as optimization condition and analyzed in three scenarios corresponding to normal, depressed and active marketing demands. Results of case study showed that, compared to independent operations, the integrated optimization and united operation among petrochemical companies could reduce material storage and overall cost, improve resource allocation, and increase operation profit.

    petrochemical industry; integration; optimization; enterprise-wide optimization; MILP; operation profit

    2016-03-24.

    Prof.RONG Gang, grong@zju.edu.cn

    10.11949/j.issn.0438-1157.20160338

    TQ 021.8

    A

    0438—1157(2016)11—4678—11

    張鵬飛(1992—),男,碩士研究生。

    國家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃項(xiàng)目(2012CB720500)。

    2016-03-24收到初稿,2016-07-29收到修改稿。

    聯(lián)系人:榮岡。

    supported by the National Basic Research Program of China (2012CB720500).

    猜你喜歡
    石腦油企業(yè)級(jí)原料
    企業(yè)級(jí)BOM數(shù)據(jù)管理概要
    石腦油輸送管道腐蝕減薄分析
    造血原料缺乏引起的貧血
    嚴(yán)把原料采購關(guān),才是對(duì)養(yǎng)殖負(fù)責(zé)
    基于慕課網(wǎng)的“企業(yè)級(jí)應(yīng)用開發(fā)”課堂教學(xué)改革探索
    電子制作(2017年10期)2017-04-18 07:23:03
    企業(yè)級(jí)信息系統(tǒng)應(yīng)用級(jí)災(zāi)備建設(shè)與應(yīng)用
    廣西電力(2016年4期)2016-07-10 10:23:38
    烘焙原料簡(jiǎn)易“識(shí)”
    美食堂(2015年5期)2015-05-30 10:48:04
    城市軌道交通企業(yè)級(jí)BIM應(yīng)用策劃研究
    以鐵泥為原料合成Fe2O3 并制備LiFePO4/C
    石腦油的分儲(chǔ)分用
    久久99热6这里只有精品| 你懂的网址亚洲精品在线观看 | 91麻豆精品激情在线观看国产| 亚洲中文字幕日韩| 人妻制服诱惑在线中文字幕| 日本一二三区视频观看| 久久人人精品亚洲av| 精品久久久久久久久久免费视频| 国产成人a∨麻豆精品| 亚洲av中文av极速乱| 精品久久久噜噜| 国产亚洲av嫩草精品影院| 99热只有精品国产| 欧美区成人在线视频| 欧美另类亚洲清纯唯美| a级毛片a级免费在线| 此物有八面人人有两片| 欧美人与善性xxx| 亚洲av电影不卡..在线观看| 国模一区二区三区四区视频| 女生性感内裤真人,穿戴方法视频| 国产黄片美女视频| 免费看光身美女| 亚洲人成网站在线播| 精品欧美国产一区二区三| 国产成人91sexporn| 悠悠久久av| 亚洲最大成人中文| 啦啦啦观看免费观看视频高清| 美女内射精品一级片tv| 久久精品影院6| 日本精品一区二区三区蜜桃| 国内揄拍国产精品人妻在线| 精华霜和精华液先用哪个| 午夜影院日韩av| 亚洲真实伦在线观看| 国产精品一区二区三区四区免费观看 | 女的被弄到高潮叫床怎么办| 男人和女人高潮做爰伦理| 国产亚洲精品久久久com| 国产淫片久久久久久久久| 国产视频内射| 国产精品久久视频播放| 亚洲av成人av| 人人妻,人人澡人人爽秒播| 一进一出抽搐gif免费好疼| 我的老师免费观看完整版| 精品久久久久久久人妻蜜臀av| 桃色一区二区三区在线观看| 国产精品久久视频播放| 国产精品爽爽va在线观看网站| 99久久无色码亚洲精品果冻| 日韩精品有码人妻一区| 国产午夜福利久久久久久| 亚洲av中文字字幕乱码综合| 日本在线视频免费播放| 日本免费a在线| 亚洲内射少妇av| 日日撸夜夜添| 日韩一区二区视频免费看| 国产中年淑女户外野战色| 九色成人免费人妻av| 婷婷六月久久综合丁香| 悠悠久久av| 亚洲成人av在线免费| 国内精品久久久久精免费| 久久草成人影院| 精品人妻熟女av久视频| 春色校园在线视频观看| 六月丁香七月| 亚洲五月天丁香| 久久天躁狠狠躁夜夜2o2o| 18禁在线无遮挡免费观看视频 | 一本精品99久久精品77| 中文字幕精品亚洲无线码一区| 夜夜爽天天搞| 国产精品久久久久久精品电影| 久久久成人免费电影| 国产精品久久久久久精品电影| 男女边吃奶边做爰视频| 一夜夜www| 一级黄片播放器| 免费看av在线观看网站| av专区在线播放| 国产毛片a区久久久久| 免费av观看视频| 看十八女毛片水多多多| 国产精品人妻久久久久久| 亚洲aⅴ乱码一区二区在线播放| 美女cb高潮喷水在线观看| 亚洲av不卡在线观看| 啦啦啦韩国在线观看视频| 免费av不卡在线播放| 在线观看av片永久免费下载| 国产高清激情床上av| 露出奶头的视频| 亚洲国产精品成人综合色| 国产精华一区二区三区| 一本久久中文字幕| 亚洲av免费高清在线观看| 干丝袜人妻中文字幕| 别揉我奶头~嗯~啊~动态视频| 12—13女人毛片做爰片一| 婷婷精品国产亚洲av在线| 少妇人妻一区二区三区视频| 男女做爰动态图高潮gif福利片| 久久精品国产清高在天天线| 免费电影在线观看免费观看| 国产精品免费一区二区三区在线| 久久人人爽人人爽人人片va| 亚洲18禁久久av| 99热只有精品国产| av免费在线看不卡| 18禁在线播放成人免费| 国产精品精品国产色婷婷| 欧美激情在线99| 欧美丝袜亚洲另类| 免费人成在线观看视频色| 在线天堂最新版资源| 国产午夜精品久久久久久一区二区三区 | 我的老师免费观看完整版| 少妇被粗大猛烈的视频| 久久人人爽人人爽人人片va| 99久久无色码亚洲精品果冻| 日韩亚洲欧美综合| aaaaa片日本免费| 亚洲图色成人| 欧美成人a在线观看| 久久精品国产自在天天线| 亚洲精品粉嫩美女一区| 卡戴珊不雅视频在线播放| 禁无遮挡网站| 色视频www国产| 成人综合一区亚洲| 国产精品1区2区在线观看.| 午夜福利18| 久久久精品94久久精品| 亚洲精品影视一区二区三区av| 给我免费播放毛片高清在线观看| 精品久久久久久成人av| 日日摸夜夜添夜夜添小说| 人人妻人人澡人人爽人人夜夜 | 亚洲精华国产精华液的使用体验 | 亚洲最大成人中文| 全区人妻精品视频| 亚洲色图av天堂| 高清午夜精品一区二区三区 | 国产一区二区三区av在线 | 国产一区二区在线观看日韩| 久久久久久国产a免费观看| 在线观看午夜福利视频| 老司机午夜福利在线观看视频| 久久久久久久久中文| 1000部很黄的大片| 菩萨蛮人人尽说江南好唐韦庄 | a级一级毛片免费在线观看| 黄色日韩在线| 国产真实伦视频高清在线观看| 久久精品国产清高在天天线| 午夜福利视频1000在线观看| 看片在线看免费视频| 午夜影院日韩av| 六月丁香七月| 赤兔流量卡办理| 色噜噜av男人的天堂激情| 18禁在线无遮挡免费观看视频 | 男女之事视频高清在线观看| 亚洲人成网站高清观看| 男女啪啪激烈高潮av片| 日日摸夜夜添夜夜爱| 亚洲在线观看片| 久久精品91蜜桃| 男女边吃奶边做爰视频| 精品久久久久久久久久久久久| 国产精品综合久久久久久久免费| 久久久精品欧美日韩精品| 成人二区视频| 成年女人看的毛片在线观看| 成年女人永久免费观看视频| 亚洲一区高清亚洲精品| 国产精品久久久久久久电影| 日本a在线网址| 亚洲图色成人| 99国产极品粉嫩在线观看| 亚洲性夜色夜夜综合| 欧美国产日韩亚洲一区| 一进一出抽搐gif免费好疼| 搡老妇女老女人老熟妇| 成人一区二区视频在线观看| 97人妻精品一区二区三区麻豆| 在线观看午夜福利视频| 伊人久久精品亚洲午夜| 国产一区二区激情短视频| 国产精品久久久久久久久免| 亚洲一区二区三区色噜噜| 国产精品精品国产色婷婷| 国产乱人偷精品视频| 精品熟女少妇av免费看| 狠狠狠狠99中文字幕| 国产大屁股一区二区在线视频| 久久久久久九九精品二区国产| 岛国在线免费视频观看| 国产亚洲精品综合一区在线观看| 菩萨蛮人人尽说江南好唐韦庄 | 我要看日韩黄色一级片| 中文字幕久久专区| 3wmmmm亚洲av在线观看| 欧美色欧美亚洲另类二区| 国产黄片美女视频| 国产精品日韩av在线免费观看| 成年av动漫网址| 国产精品一区二区性色av| 69人妻影院| 久久久久久久亚洲中文字幕| 欧美bdsm另类| 国产中年淑女户外野战色| 国产高清不卡午夜福利| 乱人视频在线观看| 国产精品一区二区免费欧美| 搡女人真爽免费视频火全软件 | 变态另类成人亚洲欧美熟女| 18禁在线播放成人免费| 日日摸夜夜添夜夜添小说| 日韩精品有码人妻一区| 三级毛片av免费| h日本视频在线播放| 综合色av麻豆| 亚洲在线自拍视频| 久久亚洲精品不卡| 久久精品国产亚洲网站| 国产精品一区二区三区四区免费观看 | 99久久精品一区二区三区| 一a级毛片在线观看| 国产欧美日韩一区二区精品| 九九在线视频观看精品| 一级毛片久久久久久久久女| 波多野结衣巨乳人妻| 美女免费视频网站| 欧美色欧美亚洲另类二区| av.在线天堂| 精品久久国产蜜桃| 国产精品一区二区免费欧美| 国产蜜桃级精品一区二区三区| 久久久久性生活片| 麻豆久久精品国产亚洲av| 狠狠狠狠99中文字幕| 日韩一区二区视频免费看| 黄片wwwwww| 成人性生交大片免费视频hd| h日本视频在线播放| 中国国产av一级| 麻豆av噜噜一区二区三区| 欧美区成人在线视频| 观看美女的网站| 国产精品电影一区二区三区| av免费在线看不卡| 久久精品国产亚洲av天美| 亚洲av中文av极速乱| 亚洲欧美清纯卡通| 国语自产精品视频在线第100页| 少妇人妻一区二区三区视频| 亚洲av.av天堂| 听说在线观看完整版免费高清| 国产麻豆成人av免费视频| 男女那种视频在线观看| 国产精品乱码一区二三区的特点| 亚洲自拍偷在线| 亚洲七黄色美女视频| 日韩欧美国产在线观看| 国产精品1区2区在线观看.| 亚洲av美国av| 亚洲精品久久国产高清桃花| 老女人水多毛片| 一区二区三区高清视频在线| av在线天堂中文字幕| 亚洲精品456在线播放app| 亚洲av二区三区四区| 色综合亚洲欧美另类图片| 在现免费观看毛片| 亚洲在线自拍视频| 人人妻,人人澡人人爽秒播| 99热只有精品国产| 午夜激情福利司机影院| 国产伦一二天堂av在线观看| 日韩精品有码人妻一区| 春色校园在线视频观看| 欧美日本视频| 男人舔女人下体高潮全视频| 精品久久久久久久末码| 国内精品久久久久精免费| 精华霜和精华液先用哪个| 最后的刺客免费高清国语| 国产激情偷乱视频一区二区| 男人舔女人下体高潮全视频| 精品日产1卡2卡| 日产精品乱码卡一卡2卡三| 国产欧美日韩一区二区精品| 成年女人看的毛片在线观看| 精品久久久久久久末码| 亚洲av二区三区四区| 午夜亚洲福利在线播放| 日韩欧美在线乱码| 淫秽高清视频在线观看| 在线天堂最新版资源| 日本色播在线视频| 亚洲精品456在线播放app| 欧美+亚洲+日韩+国产| 亚洲欧美清纯卡通| 在线天堂最新版资源| 少妇猛男粗大的猛烈进出视频 | 久久久久久九九精品二区国产| 插阴视频在线观看视频| 内射极品少妇av片p| 色哟哟·www| 精品午夜福利视频在线观看一区| 欧美日韩综合久久久久久| 久久久久久大精品| 午夜免费男女啪啪视频观看 | 国产成人freesex在线 | 国产熟女欧美一区二区| 日日干狠狠操夜夜爽| 日日摸夜夜添夜夜添av毛片| 久久久久久伊人网av| 精品久久久久久成人av| 97人妻精品一区二区三区麻豆| 在线看三级毛片| 久久热精品热| 三级毛片av免费| 悠悠久久av| 99久久精品国产国产毛片| av天堂中文字幕网| 狠狠狠狠99中文字幕| 国产精品亚洲美女久久久| 99在线人妻在线中文字幕| 日本三级黄在线观看| 3wmmmm亚洲av在线观看| 亚洲人成网站在线播| 搡老岳熟女国产| 在线免费观看不下载黄p国产| 亚洲国产高清在线一区二区三| videossex国产| 看免费成人av毛片| 中国国产av一级| 看黄色毛片网站| 成年女人永久免费观看视频| 赤兔流量卡办理| 亚洲精品粉嫩美女一区| a级一级毛片免费在线观看| 哪里可以看免费的av片| 99久久精品一区二区三区| 天堂√8在线中文| 亚洲欧美精品综合久久99| 婷婷精品国产亚洲av在线| av女优亚洲男人天堂| 在线观看免费视频日本深夜| 精品欧美国产一区二区三| 免费观看精品视频网站| 亚洲综合色惰| 亚洲国产精品久久男人天堂| 久久久久久伊人网av| 精品久久久久久久久久久久久| 久久精品国产亚洲av香蕉五月| 亚洲av不卡在线观看| 无遮挡黄片免费观看| 精品无人区乱码1区二区| 亚洲精品亚洲一区二区| 国产淫片久久久久久久久| 国产亚洲精品久久久久久毛片| 国模一区二区三区四区视频| 亚洲电影在线观看av| 男人的好看免费观看在线视频| 免费av不卡在线播放| 亚洲最大成人av| 青春草视频在线免费观看| 久久精品影院6| 99热精品在线国产| 久久久久久大精品| 亚洲,欧美,日韩| 欧美日韩乱码在线| 天天一区二区日本电影三级| 一个人观看的视频www高清免费观看| 白带黄色成豆腐渣| 欧美在线一区亚洲| 免费无遮挡裸体视频| 一个人免费在线观看电影| 久久精品国产亚洲av涩爱 | 日本黄色片子视频| 日韩中字成人| 99久久成人亚洲精品观看| avwww免费| 精品午夜福利视频在线观看一区| 男女做爰动态图高潮gif福利片| 欧美国产日韩亚洲一区| 99热只有精品国产| 美女被艹到高潮喷水动态| 成年女人永久免费观看视频| 尾随美女入室| 搡老妇女老女人老熟妇| 婷婷精品国产亚洲av| 色播亚洲综合网| av福利片在线观看| 久久久色成人| 在线观看av片永久免费下载| 国产精品一区www在线观看| 嫩草影视91久久| 国产亚洲精品久久久com| 精品一区二区三区人妻视频| 能在线免费观看的黄片| 亚洲精品一区av在线观看| 国产单亲对白刺激| 女同久久另类99精品国产91| 国产 一区精品| 18禁在线无遮挡免费观看视频 | 国产黄片美女视频| 九九热线精品视视频播放| 国产久久久一区二区三区| 男女啪啪激烈高潮av片| 18禁裸乳无遮挡免费网站照片| 中文资源天堂在线| 91麻豆精品激情在线观看国产| 欧美中文日本在线观看视频| 又粗又爽又猛毛片免费看| 99久久精品国产国产毛片| 白带黄色成豆腐渣| 美女xxoo啪啪120秒动态图| 啦啦啦韩国在线观看视频| 久久久久久久久大av| 成熟少妇高潮喷水视频| 又粗又爽又猛毛片免费看| 亚洲av免费高清在线观看| 黄色日韩在线| 最近中文字幕高清免费大全6| 日韩欧美在线乱码| 久久精品影院6| 国产在线精品亚洲第一网站| 久久99热这里只有精品18| 老司机福利观看| 2021天堂中文幕一二区在线观| 日本欧美国产在线视频| 国产精品99久久久久久久久| 久久久久久大精品| 国产精品电影一区二区三区| 免费电影在线观看免费观看| 蜜桃久久精品国产亚洲av| 美女大奶头视频| a级毛色黄片| 亚洲激情五月婷婷啪啪| 免费人成视频x8x8入口观看| 国产一区亚洲一区在线观看| 亚洲熟妇熟女久久| 夜夜夜夜夜久久久久| 国产三级中文精品| 两个人视频免费观看高清| 精品少妇黑人巨大在线播放 | 免费av观看视频| 国产精品,欧美在线| 免费黄网站久久成人精品| 成人三级黄色视频| 国产一区二区激情短视频| 免费搜索国产男女视频| h日本视频在线播放| 中国美女看黄片| 婷婷精品国产亚洲av| 亚洲婷婷狠狠爱综合网| 久久久久久久久久成人| 校园人妻丝袜中文字幕| 中国美白少妇内射xxxbb| or卡值多少钱| 久久亚洲国产成人精品v| 嫩草影院精品99| 国产又黄又爽又无遮挡在线| 亚洲乱码一区二区免费版| 最近手机中文字幕大全| 又黄又爽又免费观看的视频| 久久久久久久午夜电影| 人妻丰满熟妇av一区二区三区| 亚洲熟妇中文字幕五十中出| 老女人水多毛片| 色吧在线观看| 极品教师在线视频| 亚洲精品成人久久久久久| 午夜激情欧美在线| 久久精品91蜜桃| 久久久久九九精品影院| 日韩精品中文字幕看吧| 国产黄色小视频在线观看| 色5月婷婷丁香| 欧美性猛交╳xxx乱大交人| 亚洲丝袜综合中文字幕| 久久中文看片网| 亚洲人成网站在线播| 亚洲精品日韩av片在线观看| 狂野欧美白嫩少妇大欣赏| 亚洲成人av在线免费| 一个人观看的视频www高清免费观看| 熟妇人妻久久中文字幕3abv| 97在线视频观看| 黄色欧美视频在线观看| 色综合亚洲欧美另类图片| 内地一区二区视频在线| 五月伊人婷婷丁香| 男人狂女人下面高潮的视频| 亚洲精品亚洲一区二区| 成人av在线播放网站| 女同久久另类99精品国产91| 一个人观看的视频www高清免费观看| 又黄又爽又免费观看的视频| 麻豆国产av国片精品| 少妇人妻一区二区三区视频| 亚洲人成网站高清观看| 欧美日韩乱码在线| 午夜福利在线观看免费完整高清在 | 在线观看美女被高潮喷水网站| 亚洲av二区三区四区| 乱系列少妇在线播放| 少妇被粗大猛烈的视频| 成人一区二区视频在线观看| 一区二区三区四区激情视频 | 日日撸夜夜添| 亚洲国产日韩欧美精品在线观看| 丝袜美腿在线中文| 国产免费一级a男人的天堂| 麻豆av噜噜一区二区三区| 国产私拍福利视频在线观看| 搡老熟女国产l中国老女人| 亚洲婷婷狠狠爱综合网| 亚洲精品日韩在线中文字幕 | 国产精品亚洲美女久久久| 最后的刺客免费高清国语| 亚洲精品日韩在线中文字幕 | av福利片在线观看| 成人性生交大片免费视频hd| avwww免费| 国产av麻豆久久久久久久| 最近中文字幕高清免费大全6| 国产一级毛片七仙女欲春2| 久久精品夜色国产| www.色视频.com| 亚洲欧美成人精品一区二区| 国产成人影院久久av| 自拍偷自拍亚洲精品老妇| 国产 一区 欧美 日韩| 久久亚洲精品不卡| 亚洲中文日韩欧美视频| 97热精品久久久久久| 久久久午夜欧美精品| 亚洲图色成人| 我要看日韩黄色一级片| 亚洲电影在线观看av| 国产单亲对白刺激| 久久人人爽人人爽人人片va| 婷婷亚洲欧美| 51国产日韩欧美| 国产精品久久久久久精品电影| 亚洲精华国产精华液的使用体验 | 久久精品国产亚洲av天美| 欧美性猛交黑人性爽| 日本成人三级电影网站| 麻豆av噜噜一区二区三区| 亚洲精品色激情综合| 伊人久久精品亚洲午夜| 俺也久久电影网| 日韩欧美精品v在线| 日韩大尺度精品在线看网址| 欧美在线一区亚洲| 97碰自拍视频| 老女人水多毛片| 国产黄色小视频在线观看| 国产免费一级a男人的天堂| 欧美另类亚洲清纯唯美| 亚洲精品国产成人久久av| 国产亚洲精品久久久久久毛片| 亚洲一区高清亚洲精品| 天堂√8在线中文| 欧美激情国产日韩精品一区| 六月丁香七月| 亚洲中文字幕一区二区三区有码在线看| 国内精品宾馆在线| 在线观看免费视频日本深夜| 禁无遮挡网站| 人妻夜夜爽99麻豆av| 韩国av在线不卡| 国产精品无大码| 又爽又黄无遮挡网站| 最好的美女福利视频网| 久久草成人影院| 日韩精品有码人妻一区| 久久精品国产99精品国产亚洲性色| 国产日本99.免费观看| 婷婷精品国产亚洲av| 免费人成视频x8x8入口观看| 在线a可以看的网站| 97热精品久久久久久| 不卡一级毛片| 午夜激情欧美在线| 国产一区二区三区av在线 | 国产精品久久久久久久久免| 看非洲黑人一级黄片| 九九爱精品视频在线观看| 国产高清激情床上av| 99热精品在线国产| 少妇裸体淫交视频免费看高清| 在线观看66精品国产| 午夜福利视频1000在线观看| 日日啪夜夜撸| 中文字幕免费在线视频6| 亚洲国产欧洲综合997久久,| 国产免费男女视频| 国产成人影院久久av| 此物有八面人人有两片| 午夜精品国产一区二区电影 | 欧美激情久久久久久爽电影| 中国美女看黄片| 神马国产精品三级电影在线观看| 女的被弄到高潮叫床怎么办|