• <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ǔ)分用
    亚洲成人免费av在线播放| 精品亚洲成国产av| 手机成人av网站| 校园人妻丝袜中文字幕| 一级片'在线观看视频| 老鸭窝网址在线观看| 成在线人永久免费视频| 亚洲国产日韩一区二区| 1024香蕉在线观看| 丝袜喷水一区| 亚洲国产欧美日韩在线播放| 五月天丁香电影| 国产成人av激情在线播放| 多毛熟女@视频| 精品欧美一区二区三区在线| 精品国产一区二区久久| 别揉我奶头~嗯~啊~动态视频 | 尾随美女入室| 亚洲欧美色中文字幕在线| 女人被躁到高潮嗷嗷叫费观| 老司机在亚洲福利影院| 日韩中文字幕欧美一区二区 | 日韩一区二区三区影片| 好男人视频免费观看在线| 男女高潮啪啪啪动态图| 91国产中文字幕| 中文字幕人妻丝袜一区二区| 欧美黑人精品巨大| 日本91视频免费播放| 丝瓜视频免费看黄片| 日韩 亚洲 欧美在线| 日韩制服骚丝袜av| 麻豆国产av国片精品| 高清黄色对白视频在线免费看| 日韩视频在线欧美| 黄网站色视频无遮挡免费观看| 97人妻天天添夜夜摸| 老司机亚洲免费影院| 中文字幕精品免费在线观看视频| 好男人视频免费观看在线| 丰满迷人的少妇在线观看| 青春草亚洲视频在线观看| 国产一区二区三区av在线| 男女之事视频高清在线观看 | 久久精品人人爽人人爽视色| h视频一区二区三区| 亚洲一区中文字幕在线| 一本综合久久免费| 精品国产国语对白av| 99国产精品99久久久久| 一区二区三区四区激情视频| 欧美久久黑人一区二区| 熟女少妇亚洲综合色aaa.| 国产成人啪精品午夜网站| 日韩人妻精品一区2区三区| 欧美黄色片欧美黄色片| 在线观看免费视频网站a站| 女人高潮潮喷娇喘18禁视频| 婷婷成人精品国产| 亚洲欧美精品综合一区二区三区| 97精品久久久久久久久久精品| 丰满饥渴人妻一区二区三| 久久久久久久大尺度免费视频| 久久久久久免费高清国产稀缺| 免费日韩欧美在线观看| 精品久久久精品久久久| 一本久久精品| 精品卡一卡二卡四卡免费| 午夜福利视频精品| 黄色a级毛片大全视频| 欧美日韩亚洲国产一区二区在线观看 | 欧美精品人与动牲交sv欧美| 人人妻人人澡人人爽人人夜夜| 国产xxxxx性猛交| 日韩精品免费视频一区二区三区| 男女午夜视频在线观看| 日韩视频在线欧美| 久久精品aⅴ一区二区三区四区| 人妻 亚洲 视频| 久久精品aⅴ一区二区三区四区| 考比视频在线观看| 看十八女毛片水多多多| 丝瓜视频免费看黄片| 黑人巨大精品欧美一区二区蜜桃| 一边摸一边做爽爽视频免费| 国产老妇伦熟女老妇高清| 久久精品亚洲熟妇少妇任你| 91九色精品人成在线观看| 国产av国产精品国产| 在线观看国产h片| 国产片内射在线| 国产精品欧美亚洲77777| 大码成人一级视频| 成人国语在线视频| 男人操女人黄网站| 国产精品偷伦视频观看了| 又大又爽又粗| 日韩人妻精品一区2区三区| av网站免费在线观看视频| 国产伦人伦偷精品视频| 亚洲成色77777| 国产一区二区三区综合在线观看| 赤兔流量卡办理| 午夜福利,免费看| 亚洲精品一二三| 亚洲国产精品999| 91麻豆精品激情在线观看国产 | 欧美成狂野欧美在线观看| 超色免费av| 日韩av在线免费看完整版不卡| 美女大奶头黄色视频| 欧美精品人与动牲交sv欧美| 国产精品av久久久久免费| 精品国产一区二区三区久久久樱花| 久久鲁丝午夜福利片| 黄色视频不卡| 十八禁网站网址无遮挡| 99re6热这里在线精品视频| 精品人妻熟女毛片av久久网站| 亚洲成人国产一区在线观看 | 首页视频小说图片口味搜索 | 中文字幕制服av| 999久久久国产精品视频| av线在线观看网站| av天堂在线播放| 精品国产一区二区三区四区第35| 成年动漫av网址| 欧美黑人欧美精品刺激| 一级毛片黄色毛片免费观看视频| 免费观看av网站的网址| 免费看十八禁软件| 精品久久久久久电影网| 中文字幕最新亚洲高清| 看免费av毛片| 亚洲国产欧美网| 在线 av 中文字幕| 国产精品国产三级国产专区5o| 午夜免费成人在线视频| 久久精品熟女亚洲av麻豆精品| 另类亚洲欧美激情| 欧美老熟妇乱子伦牲交| 国产一区二区激情短视频 | 久久久欧美国产精品| 久久亚洲国产成人精品v| 啦啦啦 在线观看视频| 亚洲欧美一区二区三区黑人| 在线观看免费日韩欧美大片| 国产成人欧美在线观看 | 国产高清国产精品国产三级| 91九色精品人成在线观看| 国产视频一区二区在线看| 侵犯人妻中文字幕一二三四区| 丝瓜视频免费看黄片| 另类亚洲欧美激情| av有码第一页| 色精品久久人妻99蜜桃| 成人黄色视频免费在线看| 如日韩欧美国产精品一区二区三区| 777米奇影视久久| 国产高清videossex| 国产成人免费无遮挡视频| 免费日韩欧美在线观看| 妹子高潮喷水视频| 精品人妻1区二区| 亚洲熟女精品中文字幕| 两性夫妻黄色片| 99热网站在线观看| 精品一区在线观看国产| 热re99久久国产66热| 国产精品99久久99久久久不卡| 欧美精品高潮呻吟av久久| av视频免费观看在线观看| 天天添夜夜摸| 丝袜人妻中文字幕| √禁漫天堂资源中文www| 亚洲 欧美一区二区三区| 超碰成人久久| 亚洲av综合色区一区| 在线观看人妻少妇| 国产日韩欧美亚洲二区| 国产精品久久久人人做人人爽| 免费高清在线观看日韩| 嫁个100分男人电影在线观看 | 午夜激情av网站| 欧美成人精品欧美一级黄| 国产xxxxx性猛交| 日韩熟女老妇一区二区性免费视频| 又大又爽又粗| 国产高清videossex| 久久精品亚洲熟妇少妇任你| 午夜激情av网站| 免费在线观看影片大全网站 | 19禁男女啪啪无遮挡网站| 五月开心婷婷网| 久久精品国产亚洲av高清一级| 久久天躁狠狠躁夜夜2o2o | 男女午夜视频在线观看| www.av在线官网国产| 中文欧美无线码| 日日摸夜夜添夜夜爱| 久久久久久免费高清国产稀缺| 热99久久久久精品小说推荐| 久久99一区二区三区| 亚洲av男天堂| 免费一级毛片在线播放高清视频 | 王馨瑶露胸无遮挡在线观看| 一级毛片我不卡| 日韩 亚洲 欧美在线| 亚洲国产成人一精品久久久| 国产一卡二卡三卡精品| 高清黄色对白视频在线免费看| 日韩伦理黄色片| av一本久久久久| 欧美日韩国产mv在线观看视频| 建设人人有责人人尽责人人享有的| 久久精品国产a三级三级三级| www.999成人在线观看| 这个男人来自地球电影免费观看| 国产黄频视频在线观看| 亚洲激情五月婷婷啪啪| 亚洲精品日韩在线中文字幕| 欧美精品一区二区大全| 亚洲av日韩在线播放| 久久天躁狠狠躁夜夜2o2o | 国产一级毛片在线| 日本91视频免费播放| 国产伦人伦偷精品视频| 成年女人毛片免费观看观看9 | 日韩伦理黄色片| 在线观看免费日韩欧美大片| 中文字幕人妻熟女乱码| 亚洲国产欧美在线一区| 高清欧美精品videossex| 韩国精品一区二区三区| 国产精品一区二区在线观看99| 少妇人妻 视频| 国产日韩一区二区三区精品不卡| 丝袜美足系列| 一二三四在线观看免费中文在| 嫩草影视91久久| 一区二区日韩欧美中文字幕| 天天添夜夜摸| 一本—道久久a久久精品蜜桃钙片| 精品第一国产精品| 欧美性长视频在线观看| av天堂在线播放| 亚洲精品美女久久久久99蜜臀 | a级毛片黄视频| 无遮挡黄片免费观看| 99国产精品一区二区蜜桃av | 欧美日韩视频精品一区| xxx大片免费视频| 国产成人精品久久二区二区免费| 亚洲三区欧美一区| 日韩视频在线欧美| 亚洲国产日韩一区二区| 成年人黄色毛片网站| 脱女人内裤的视频| 中文字幕亚洲精品专区| 免费观看人在逋| 欧美精品av麻豆av| 久久人人爽人人片av| 国产av一区二区精品久久| 久久99一区二区三区| 亚洲国产欧美在线一区| 中文乱码字字幕精品一区二区三区| 丝袜在线中文字幕| 91精品三级在线观看| 亚洲熟女精品中文字幕| 久久久久久久久免费视频了| 男人舔女人的私密视频| 国产伦理片在线播放av一区| 久久精品人人爽人人爽视色| 婷婷色麻豆天堂久久| 首页视频小说图片口味搜索 | 国产xxxxx性猛交| 大型av网站在线播放| 日韩大码丰满熟妇| 国产亚洲欧美精品永久| 免费一级毛片在线播放高清视频 | 青春草视频在线免费观看| 亚洲欧美色中文字幕在线| 亚洲久久久国产精品| 亚洲av欧美aⅴ国产| 国产成人91sexporn| 丝袜美足系列| 黄网站色视频无遮挡免费观看| 国产野战对白在线观看| 国产一卡二卡三卡精品| 国产欧美亚洲国产| 高清黄色对白视频在线免费看| 亚洲伊人久久精品综合| av国产精品久久久久影院| 麻豆乱淫一区二区| 18禁观看日本| 亚洲国产精品999| 久久精品成人免费网站| 自拍欧美九色日韩亚洲蝌蚪91| 大片免费播放器 马上看| 性色av一级| 精品视频人人做人人爽| 国产亚洲精品久久久久5区| 亚洲av成人精品一二三区| 免费在线观看日本一区| 国产精品一区二区精品视频观看| 午夜两性在线视频| 美女午夜性视频免费| 99精国产麻豆久久婷婷| 欧美日韩一级在线毛片| 青青草视频在线视频观看| 欧美 亚洲 国产 日韩一| 午夜视频精品福利| 午夜福利,免费看| 中国美女看黄片| av网站在线播放免费| 亚洲自偷自拍图片 自拍| 乱人伦中国视频| 中文字幕人妻丝袜一区二区| 麻豆av在线久日| 久久久精品94久久精品| 亚洲精品av麻豆狂野| 午夜视频精品福利| 欧美黑人欧美精品刺激| 男女边摸边吃奶| 自拍欧美九色日韩亚洲蝌蚪91| 少妇猛男粗大的猛烈进出视频| 最近手机中文字幕大全| 乱人伦中国视频| 男女之事视频高清在线观看 | 好男人视频免费观看在线| 国产欧美日韩一区二区三 | 午夜免费观看性视频| 在线观看免费日韩欧美大片| 脱女人内裤的视频| 精品国产超薄肉色丝袜足j| 另类亚洲欧美激情| 这个男人来自地球电影免费观看| 国产精品一区二区免费欧美 | 国产精品免费视频内射| 国产精品麻豆人妻色哟哟久久| 精品福利永久在线观看| 国产av国产精品国产| 自线自在国产av| 一二三四社区在线视频社区8| 大香蕉久久成人网| av一本久久久久| 久热爱精品视频在线9| 久久久精品免费免费高清| 满18在线观看网站| 国产亚洲欧美精品永久| 亚洲三区欧美一区| 香蕉丝袜av| 啦啦啦在线观看免费高清www| 久久久久国产精品人妻一区二区| 一区二区日韩欧美中文字幕| 久久精品亚洲av国产电影网| 亚洲情色 制服丝袜| 亚洲精品在线美女| 亚洲,欧美精品.| 成人手机av| av在线播放精品| 超碰97精品在线观看| 九草在线视频观看| 如日韩欧美国产精品一区二区三区| 50天的宝宝边吃奶边哭怎么回事| 悠悠久久av| 国产在线视频一区二区| 99香蕉大伊视频| 高清欧美精品videossex| 精品国产一区二区三区四区第35| 丝袜人妻中文字幕| 丁香六月欧美| 亚洲国产看品久久| 1024香蕉在线观看| 日本av免费视频播放| 在线观看国产h片| 男人爽女人下面视频在线观看| 亚洲精品久久成人aⅴ小说| 伊人亚洲综合成人网| 成年女人毛片免费观看观看9 | 久久久精品国产亚洲av高清涩受| 亚洲 欧美一区二区三区| 一区福利在线观看| 欧美成狂野欧美在线观看| 日日摸夜夜添夜夜爱| 无遮挡黄片免费观看| 男女免费视频国产| 一级毛片电影观看| 国产爽快片一区二区三区| 久久亚洲国产成人精品v| 亚洲一卡2卡3卡4卡5卡精品中文| 少妇猛男粗大的猛烈进出视频| 成人国产一区最新在线观看 | 青青草视频在线视频观看| 久久久久久免费高清国产稀缺| 看免费成人av毛片| 免费在线观看黄色视频的| 久久国产亚洲av麻豆专区| 新久久久久国产一级毛片| 久久人妻熟女aⅴ| 日韩av免费高清视频| 黄色 视频免费看| 在线看a的网站| 午夜福利视频精品| 国产91精品成人一区二区三区 | 1024香蕉在线观看| 国产精品秋霞免费鲁丝片| 在线观看www视频免费| 香蕉丝袜av| 亚洲免费av在线视频| 一二三四在线观看免费中文在| av有码第一页| 国产男女内射视频| 国产爽快片一区二区三区| 亚洲国产欧美一区二区综合| 国产成人免费无遮挡视频| 精品人妻在线不人妻| 日韩一卡2卡3卡4卡2021年| 久久人人爽人人片av| 人体艺术视频欧美日本| 制服人妻中文乱码| 免费黄频网站在线观看国产| 日本av免费视频播放| 中文乱码字字幕精品一区二区三区| 国产欧美亚洲国产| 黑人猛操日本美女一级片| 亚洲少妇的诱惑av| 成年人黄色毛片网站| 五月天丁香电影| 国产一区二区激情短视频 | 欧美大码av| 国产精品麻豆人妻色哟哟久久| 日本色播在线视频| 欧美国产精品va在线观看不卡| 十八禁网站网址无遮挡| 亚洲精品国产一区二区精华液| 中文字幕亚洲精品专区| 久久天躁狠狠躁夜夜2o2o | 国产精品av久久久久免费| 9色porny在线观看| 国产真人三级小视频在线观看| 国产国语露脸激情在线看| 国产精品98久久久久久宅男小说| 欧美性猛交╳xxx乱大交人| 制服人妻中文乱码| www日本黄色视频网| 免费在线观看黄色视频的| 国产黄片美女视频| 久久中文字幕人妻熟女| 中文字幕久久专区| 久久99热这里只有精品18| 日本撒尿小便嘘嘘汇集6| 麻豆久久精品国产亚洲av| 国产成人精品无人区| 亚洲第一电影网av| 亚洲国产精品sss在线观看| 日本一区二区免费在线视频| 国产精品免费视频内射| 午夜亚洲福利在线播放| 成人一区二区视频在线观看| av欧美777| 成人国语在线视频| 日本免费一区二区三区高清不卡| 1024视频免费在线观看| 国产伦在线观看视频一区| 国产欧美日韩精品亚洲av| 国产高清videossex| 99re在线观看精品视频| 三级毛片av免费| 欧美亚洲日本最大视频资源| 亚洲一码二码三码区别大吗| 搡老熟女国产l中国老女人| 免费高清视频大片| 国产精品永久免费网站| 久久国产精品男人的天堂亚洲| 亚洲欧美激情综合另类| 久久久久久人人人人人| 日韩大码丰满熟妇| 精品福利观看| 免费观看人在逋| 真人做人爱边吃奶动态| 免费在线观看完整版高清| 亚洲国产精品久久男人天堂| 色在线成人网| 中文字幕精品亚洲无线码一区 | 波多野结衣高清作品| 国产片内射在线| 一a级毛片在线观看| 久久久国产成人精品二区| 18禁观看日本| 国产熟女午夜一区二区三区| 一级黄色大片毛片| 人成视频在线观看免费观看| 不卡av一区二区三区| 99国产综合亚洲精品| 久久久久久久久中文| 亚洲国产欧洲综合997久久, | 亚洲国产看品久久| 少妇的丰满在线观看| 男女午夜视频在线观看| 热99re8久久精品国产| 国产高清视频在线播放一区| 一级作爱视频免费观看| 悠悠久久av| 久久久国产成人免费| 亚洲三区欧美一区| 一区二区三区精品91| 精品不卡国产一区二区三区| 欧美成人性av电影在线观看| 一本综合久久免费| 俄罗斯特黄特色一大片| 中文字幕高清在线视频| 国产亚洲欧美在线一区二区| 精品国产乱子伦一区二区三区| 岛国视频午夜一区免费看| 国产又黄又爽又无遮挡在线| 免费av毛片视频| 校园春色视频在线观看| 国产在线观看jvid| 两人在一起打扑克的视频| av天堂在线播放| 亚洲欧美激情综合另类| 午夜免费成人在线视频| 亚洲人成伊人成综合网2020| 日韩精品青青久久久久久| 18禁黄网站禁片午夜丰满| 亚洲专区中文字幕在线| 久99久视频精品免费| 18禁观看日本| 亚洲成人久久爱视频| 一进一出抽搐gif免费好疼| 好男人在线观看高清免费视频 | 国产精品久久久久久人妻精品电影| 欧美激情久久久久久爽电影| 成人午夜高清在线视频 | 成人午夜高清在线视频 | 国产高清有码在线观看视频 | 高清毛片免费观看视频网站| 两性午夜刺激爽爽歪歪视频在线观看 | 国产精品久久电影中文字幕| 日韩免费av在线播放| 俺也久久电影网| 日韩精品中文字幕看吧| 亚洲人成77777在线视频| 黄色a级毛片大全视频| 男女做爰动态图高潮gif福利片| 国产亚洲欧美98| 日本免费一区二区三区高清不卡| 这个男人来自地球电影免费观看| 超碰成人久久| 欧美日韩中文字幕国产精品一区二区三区| 一本精品99久久精品77| 中文字幕人成人乱码亚洲影| aaaaa片日本免费| 久久久久亚洲av毛片大全| 国产麻豆成人av免费视频| 日韩免费av在线播放| 午夜精品久久久久久毛片777| 午夜两性在线视频| 日韩高清综合在线| 国产精品亚洲一级av第二区| 午夜福利在线观看吧| 一级a爱视频在线免费观看| 成人午夜高清在线视频 | 黄色视频,在线免费观看| 人人妻人人看人人澡| 成人免费观看视频高清| 日日夜夜操网爽| 免费高清视频大片| 脱女人内裤的视频| 亚洲精品粉嫩美女一区| 黄色毛片三级朝国网站| 亚洲熟妇熟女久久| 窝窝影院91人妻| 亚洲中文av在线| 日日摸夜夜添夜夜添小说| av有码第一页| 波多野结衣巨乳人妻| 亚洲国产欧洲综合997久久, | 欧美日本亚洲视频在线播放| 精品欧美国产一区二区三| 在线观看www视频免费| 精品欧美国产一区二区三| 12—13女人毛片做爰片一| 亚洲人成电影免费在线| 伊人久久大香线蕉亚洲五| 久久这里只有精品19| 免费看a级黄色片| 十分钟在线观看高清视频www| 久久国产乱子伦精品免费另类| www国产在线视频色| 97人妻精品一区二区三区麻豆 | netflix在线观看网站| or卡值多少钱| 国产成+人综合+亚洲专区| 久久精品国产亚洲av高清一级| 欧美 亚洲 国产 日韩一| cao死你这个sao货| 亚洲色图av天堂| 18禁裸乳无遮挡免费网站照片 | 丰满的人妻完整版| 国产精品一区二区免费欧美| 99在线人妻在线中文字幕| 99国产精品99久久久久| 国产精品av久久久久免费| 天天一区二区日本电影三级| av视频在线观看入口| 一级黄色大片毛片| 制服诱惑二区| 此物有八面人人有两片| 桃色一区二区三区在线观看| 巨乳人妻的诱惑在线观看| 熟女少妇亚洲综合色aaa.| 精品国产美女av久久久久小说| 91老司机精品|