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

    原油調(diào)度優(yōu)化中的分段線性松弛研究

    2016-04-12 02:09:05周智菊郭錦標(biāo)
    石油煉制與化工 2016年8期
    關(guān)鍵詞:廠區(qū)分段線性

    周 祥,周智菊,呂 寧,郭錦標(biāo)

    (中國石化石油化工科學(xué)研究院,北京 100083)

    原油調(diào)度優(yōu)化中的分段線性松弛研究

    周 祥,周智菊,呂 寧,郭錦標(biāo)

    (中國石化石油化工科學(xué)研究院,北京 100083)

    針對(duì)混合原油性質(zhì)計(jì)算導(dǎo)致原油調(diào)度優(yōu)化模型中出現(xiàn)非線性約束并嚴(yán)重影響模型求解性能的情況,采用分段線性松弛法將非線性約束轉(zhuǎn)化為線性形式,使模型可應(yīng)對(duì)原油多次混合的復(fù)雜情況。以實(shí)際煉油廠的原油調(diào)度優(yōu)化為例,對(duì)混合原油硫含量計(jì)算作線性轉(zhuǎn)化,并將分段線性松弛法與兩種常見方法的應(yīng)用效果進(jìn)行了對(duì)比。計(jì)算結(jié)果表明:分段線性松弛法可使模型的求解性能達(dá)到工業(yè)應(yīng)用要求;當(dāng)加工油種的硫質(zhì)量分?jǐn)?shù)為0.50%~3.00%時(shí),分段線性松弛法更有利于獲得優(yōu)化的調(diào)度方案;當(dāng)加工油種硫質(zhì)量分?jǐn)?shù)的差異由0.50%增加至2.50%時(shí),分段線性松弛法對(duì)廠區(qū)罐內(nèi)混合原油硫含量的計(jì)算值均接近實(shí)際值,相對(duì)偏差低于5%,可確保模型的準(zhǔn)確度滿足工業(yè)應(yīng)用要求。

    原油調(diào)度 優(yōu)化模型 分段線性松弛 硫含量

    國內(nèi)煉油廠加工的原油品種較多、性質(zhì)差異較大,多油種混煉是操作常態(tài),并時(shí)常摻煉劣質(zhì)原油和機(jī)會(huì)原油以提高效益,因此原油調(diào)合越來越受到煉油廠重視,成為其加工過程的第一道工序。原油調(diào)度優(yōu)化可提供合理的原油調(diào)合操作方案,保障蒸餾裝置及后續(xù)裝置的平穩(wěn)運(yùn)行并緩解庫存壓力[1]。利用數(shù)學(xué)模型實(shí)現(xiàn)原油調(diào)度優(yōu)化具有可行性[2-7],然而實(shí)際原油儲(chǔ)運(yùn)過程比較復(fù)雜,期間一般發(fā)生多次原油混合,使混合原油的性質(zhì)計(jì)算呈非線性,模型中存在此類約束時(shí)會(huì)嚴(yán)重影響其求解性能[3-5]。

    求解含非線性約束的原油調(diào)度優(yōu)化模型多采用迭代法[6-7],此方法的求解效率取決于迭代初值的設(shè)置是否合理,而調(diào)整初值主要依靠調(diào)度人員的經(jīng)驗(yàn),自動(dòng)化程度較低。約束規(guī)劃建模方法可允許模型中出現(xiàn)較多的非線性約束[5],其優(yōu)勢(shì)在于快速獲得可行的調(diào)度方案,模型中一旦設(shè)置優(yōu)化目標(biāo)則求解效率同樣較低。本課題采用分段線性松弛法將模型中混合原油性質(zhì)的計(jì)算轉(zhuǎn)化為線性形式,并與兩種常見的處理方法進(jìn)行比較,考察線性松弛對(duì)模型求解性能及準(zhǔn)確度的影響。

    1 非線性約束轉(zhuǎn)化

    1.1 非線性約束

    沿海煉油廠的原油儲(chǔ)運(yùn)流程較長,原油調(diào)合操作復(fù)雜,其原油調(diào)度優(yōu)化問題具有代表性。圖1為沿海煉油廠原油調(diào)度過程示意。由圖1可知,沿海煉油廠的原油調(diào)度需考慮油輪到港卸油、碼頭罐收油、原油經(jīng)長輸管線付至廠區(qū)罐、廠區(qū)罐向蒸餾裝置供料的全過程。一般而言,各品種的原油在碼頭罐內(nèi)單儲(chǔ)以便于收付計(jì)量,多個(gè)碼頭罐付油到長輸管線時(shí)發(fā)生原油的第一次混合,廠區(qū)罐接收多段由長輸管線送來的混合原油時(shí)發(fā)生第二次混合,國內(nèi)煉油廠多在蒸餾裝置進(jìn)料環(huán)節(jié)進(jìn)行劣質(zhì)原油和機(jī)會(huì)原油的摻煉,即存在多個(gè)廠區(qū)罐同時(shí)向一套蒸餾裝置供料的情況,這將導(dǎo)致第三次原油混合。

    圖1 沿海煉油廠原油調(diào)度過程示意

    原油調(diào)度優(yōu)化的主要目的是確保廠區(qū)罐所存混合原油的性質(zhì)滿足蒸餾裝置進(jìn)料要求,并盡可能減少進(jìn)料性質(zhì)波動(dòng)。Lee等[3]將硫等影響混合原油性質(zhì)的因素定義為“關(guān)鍵組分”,直接對(duì)這些“關(guān)鍵組分”在廠區(qū)罐內(nèi)的存量進(jìn)行物料衡算,通過限制廠區(qū)罐向蒸餾裝置供應(yīng)“關(guān)鍵組分”的數(shù)量來滿足蒸餾裝置的進(jìn)料性質(zhì)要求,如式(1)和式(2)所示。

    WT,S=WT,S,0+WT,S,I-WT,S,O

    (1)

    (2)

    式(1)和式(2)均呈線性,而實(shí)際上廠區(qū)罐供給蒸餾裝置的是混合原油,其中所含“關(guān)鍵組分”的數(shù)量受廠區(qū)罐內(nèi)混合原油性質(zhì)的約束,如式(3)所示。

    WT,S,O=WT,O·ST

    (3)

    式中:WT,O為廠區(qū)罐付出混合原油的質(zhì)量,t;ST為廠區(qū)罐內(nèi)混合原油的硫質(zhì)量分?jǐn)?shù),%。

    廠區(qū)罐內(nèi)原油二次混合的情況在沿海煉油廠的儲(chǔ)運(yùn)過程中普遍存在,廠區(qū)罐從長輸管線接收“關(guān)鍵組分”的數(shù)量可表示為式(4)。

    WT,S,I=∑(WP,O·SP)

    (4)

    式中:WP,O為長輸管線付出混合原油的質(zhì)量,t;SP為長輸管線付出混合原油的硫質(zhì)量分?jǐn)?shù),%。

    式(3)和式(4)均呈非線性,導(dǎo)致原油調(diào)度優(yōu)化模型中出現(xiàn)非線性約束,而劣質(zhì)原油和機(jī)會(huì)原油的摻煉將使模型中的非線性約束進(jìn)一步增加。

    1.2 分段線性松弛法

    “關(guān)鍵組分”法為保障模型的求解效率,回避了式(3)和式(4)所示的非線性約束[3],然而求解模型所獲得的調(diào)度方案中,廠區(qū)罐內(nèi)的混合原油性質(zhì)及蒸餾裝置的進(jìn)料性質(zhì)與實(shí)際情況往往存在較大的偏差,嚴(yán)重時(shí)會(huì)導(dǎo)致方案不具備可操作性。

    由式(3)和式(4)可知,模型中的非線性約束均具有雙線性形式(bilinear)。有研究者提出將此類雙線性約束轉(zhuǎn)化為一系列線性約束,這種處理方法會(huì)使模型中的約束更松弛,因此松弛模型的最優(yōu)解可能不同于原模型(松弛模型的最優(yōu)解通常更加優(yōu)化)。然而,松弛模型不包含非線性約束,其求解效率大大高于原模型,若能減小松弛模型與原模型的差異,則其優(yōu)化解的準(zhǔn)確性亦有保證[8-9]。以前的研究曾嘗試?yán)肅onvex and Concave Envelopes(簡稱CCE)法實(shí)現(xiàn)非線性約束的線性轉(zhuǎn)化[10],但此方法對(duì)雙線性約束的松弛不夠緊致,且國內(nèi)煉油廠所加工原油的性質(zhì)差異較大,因此對(duì)混合原油性質(zhì)的計(jì)算也存在較大的偏差。

    為實(shí)現(xiàn)更為緊致的線性松弛,本研究對(duì)雙線性約束中變量的取值區(qū)間進(jìn)行劃分,并據(jù)此對(duì)變量乘積作分段松弛。以式(3)所示雙線性約束為例,將廠區(qū)罐內(nèi)混合原油硫含量的取值區(qū)間均分為N段,各段端點(diǎn)值如式(5)所示。

    (5)

    (6)

    (7)

    從廠區(qū)罐向蒸餾裝置的供料量與蒸餾裝置的加工量相關(guān)聯(lián),在實(shí)際操作中蒸餾裝置的加工量一般比較穩(wěn)定,因此對(duì)供料量的取值區(qū)間進(jìn)行分段意義不大。定義變量ΔWT,O,n,與硫含量各子區(qū)間相對(duì)應(yīng),當(dāng)BT,S,n= 1時(shí)ΔWT,O,n為供料量相對(duì)其下限值的偏差,以式(8)和式(9)約束。

    (8)

    (9)

    對(duì)式(3)進(jìn)行線性松弛,如式(10)~式(13)所示。

    (10)

    (11)

    (12)

    (13)

    由式(10)~式(13)可知,當(dāng)N= 1時(shí),此分段線性松弛法與CCE法是等效的,而不同之處在于,此方法為使線性松弛更緊致而在模型中引入了新的01變量,增加01變量也可能會(huì)降低模型的求解效率。目前,混合整數(shù)線性規(guī)劃模型求解算法和軟件工具的開發(fā)比非線性規(guī)劃模型更為完善,對(duì)于大規(guī)模工業(yè)應(yīng)用級(jí)的模型,其優(yōu)勢(shì)更加明顯,因此01變量對(duì)模型求解性能的影響弱于非線性約束,而調(diào)整N的取值也可平衡松弛緊致程度和01變量數(shù)。

    2 計(jì)算結(jié)果與討論

    2.1 案例介紹

    以某沿海煉油廠的原油調(diào)度優(yōu)化為案例,考察線性松弛對(duì)模型求解性能及準(zhǔn)確度的影響。該煉油廠以含硫和高硫原油為主要加工油種,儲(chǔ)運(yùn)系統(tǒng)包括1個(gè)油輪??坎次?、10個(gè)碼頭罐、1條長輸管線、8個(gè)廠區(qū)罐,原油在廠區(qū)罐內(nèi)進(jìn)行二次混合,有時(shí)在蒸餾裝置進(jìn)料環(huán)節(jié)摻煉劣質(zhì)原油和機(jī)會(huì)原油。

    以180 h為一個(gè)調(diào)度周期,采用異步連續(xù)時(shí)間模式建立該煉油廠的原油調(diào)度優(yōu)化模型,綜合考慮油輪滯船期、蒸餾裝置進(jìn)料性質(zhì)波動(dòng)及儲(chǔ)罐收付切換次數(shù)等指標(biāo),并以其加權(quán)和最小為優(yōu)化目標(biāo)[10]。蒸餾裝置的進(jìn)料性質(zhì)要求重點(diǎn)考慮硫含量指標(biāo),并選擇不同的初始狀態(tài)使調(diào)度周期內(nèi)加工油種的硫含量位于不同的區(qū)間。為比較模型的求解性能及準(zhǔn)確度,分別應(yīng)用“關(guān)鍵組分”法、CCE法與上文所述分段線性松弛法對(duì)模型進(jìn)行轉(zhuǎn)化,在分段線性松弛法中對(duì)混合原油硫含量劃分10個(gè)子區(qū)間(N= 10)。

    2.2 模型規(guī)模及求解性能

    模型的求解環(huán)境為:32個(gè)主頻為2.8 GHz的CPU核心、128 G內(nèi)存、64位XPRESS 7.5混合整數(shù)線性規(guī)劃求解器。表1為3種方法所得模型的規(guī)模(均為預(yù)求解后),其中“關(guān)鍵組分”法忽略了非線性約束,因此所得模型規(guī)模最小;CCE法并未使模型中增加新的變量,而是將非線性約束松弛化使模型中的約束增加;分段線性松弛法最為復(fù)雜,使模型規(guī)模明顯擴(kuò)大,尤其是01變量增加了60%以上。

    表1 3種方法所得模型的規(guī)模

    當(dāng)加工油種的硫含量位于不同區(qū)間時(shí),3種方法所得模型的求解情況如表2所示。工業(yè)應(yīng)用對(duì)模型的求解效率有較高的要求,因此將模型求解時(shí)間上限設(shè)為1 h,若求解時(shí)間達(dá)到上限時(shí)仍未獲得最優(yōu)解,則記錄當(dāng)前優(yōu)化解,若在求解時(shí)間上限之前獲得最優(yōu)解,則記錄計(jì)算時(shí)間和最優(yōu)解。

    表2 3種方法所得模型的求解情況

    由表2可知:3種方法所得模型均可在1 h內(nèi)獲得優(yōu)化解,求解性能均滿足工業(yè)應(yīng)用要求;當(dāng)加工油種的硫質(zhì)量分?jǐn)?shù)差異較小時(shí)(≤1.00%),“關(guān)鍵組分”法和CCE法所得模型的求解效率更高,較易確定最優(yōu)解,然而當(dāng)加工油種的硫質(zhì)量分?jǐn)?shù)位于2.00%~3.00%時(shí),CCE法已確定最優(yōu)方案的優(yōu)化目標(biāo)值不低于130,而“關(guān)鍵組分”法則已將最優(yōu)目標(biāo)值降至130以下,可見后者獲得的優(yōu)化調(diào)度方案未必可靠;當(dāng)加工油種的硫質(zhì)量分?jǐn)?shù)差異較大時(shí)(>1.00%),“關(guān)鍵組分”法和CCE法所得模型的求解效率下降,達(dá)到求解時(shí)間上限時(shí)優(yōu)化目標(biāo)值相對(duì)較差;分段線性松弛法所得模型也可確定最優(yōu)解,同時(shí)其求解效率對(duì)加工油種硫含量的差異不敏感,而且當(dāng)差異較大時(shí),應(yīng)用此方法在相同計(jì)算時(shí)間內(nèi)獲得的優(yōu)化解質(zhì)量最好,表明其具有更廣泛的適用性。

    此外,雖然分段線性松弛使模型中的01變量大幅增加,但模型的求解性能仍然較好,其原因在于01變量的取值是通過分支定界算法確定的,分支定界過程在求解環(huán)境內(nèi)以并行計(jì)算的方式完成,從而緩和了01變量規(guī)模對(duì)模型求解效率的影響。

    2.3 模型準(zhǔn)確度

    準(zhǔn)確計(jì)算混合原油的性質(zhì)需采用非線性約束,而線性松弛法將式(3)和式(4)中的非線性等式約束轉(zhuǎn)化為線性不等式約束,因此可能導(dǎo)致混合原油性質(zhì)的計(jì)算值與實(shí)際值之間存在偏差。為考察模型的準(zhǔn)確度,對(duì)應(yīng)用3種方法獲得的優(yōu)化調(diào)度方案進(jìn)行驗(yàn)算,依據(jù)方案中各加工油種的儲(chǔ)運(yùn)過程確定廠區(qū)罐內(nèi)混合原油硫含量的實(shí)際值,將計(jì)算值與實(shí)際值進(jìn)行比較,最大相對(duì)偏差(|計(jì)算值-實(shí)際值|/實(shí)際值×100%)如表3所示。

    表3 3種方法所得模型計(jì)算廠區(qū)罐混合原油硫含量的準(zhǔn)確度

    由表3可知:“關(guān)鍵組分”法所得模型的準(zhǔn)確度較差,當(dāng)加工油種的硫質(zhì)量分?jǐn)?shù)差異較小(0.50%~1.00 %)或較大(0.50%~3.00%)時(shí),其計(jì)算偏差可能達(dá)50%以上,這是忽略混合原油性質(zhì)計(jì)算造成的影響,應(yīng)用此方法獲得的調(diào)度方案不具備可操作性;CCE法所得模型的準(zhǔn)確度高于“關(guān)鍵組分”法,然而其計(jì)算偏差同樣受限于加工油種硫含量的差異;與CCE法相比,分段線性松弛法對(duì)混合原油性質(zhì)計(jì)算進(jìn)行了更為緊致的轉(zhuǎn)化,因此所得模型的準(zhǔn)確度在3種方法中是最高的,并且當(dāng)加工油種的硫含量位于不同區(qū)間時(shí),計(jì)算偏差均在5%以下,可滿足工業(yè)應(yīng)用要求,獲得的調(diào)度方案具備可操作性。

    3 結(jié) 論

    采用分段線性松弛法將混合原油性質(zhì)計(jì)算轉(zhuǎn)化為線性形式,避免原油調(diào)度優(yōu)化模型中出現(xiàn)非線性約束,適用于煉油廠摻煉劣質(zhì)原油和機(jī)會(huì)原油的復(fù)雜情況。以實(shí)際煉油廠的原油調(diào)度優(yōu)化為例,對(duì)混合原油硫含量計(jì)算作線性轉(zhuǎn)化,將分段線性松弛法與兩種常見方法的應(yīng)用效果進(jìn)行了對(duì)比,計(jì)算結(jié)果表明:應(yīng)用分段線性松弛法可使模型的求解性能達(dá)到工業(yè)應(yīng)用要求;當(dāng)加工油種的硫質(zhì)量分?jǐn)?shù)為0.50%~3.00%時(shí),分段線性松弛法更有利于獲得優(yōu)化的調(diào)度方案,具有更廣泛的適用性;當(dāng)加工油種硫質(zhì)量分?jǐn)?shù)的差異由0.50%增至2.50%時(shí),分段線性松弛法始終能將廠區(qū)罐內(nèi)混合原油硫含量的計(jì)算偏差控制在5%以下,確保模型的準(zhǔn)確度滿足工業(yè)應(yīng)用要求,所獲得的調(diào)度方案具備可操作性。

    [1] 郭錦標(biāo),楊明詩.化工生產(chǎn)計(jì)劃與調(diào)度的優(yōu)化[M].北京:化學(xué)工業(yè)出版社,2006:163-165

    [2] Shah N.Mathematical programming techniques for crude oil scheduling[J].Computers & Chemical Engineering,1996,20(S):1227-1232

    [3] Lee H,Pinto J M,Grossmann I E,et al.Mix-integer linear programming model for refinery short-term scheduling of crude oil unloading with inventory management[J].Industrial & Engineering Chemistry Research,1996,35(5):1630-1641

    [4] Jia Z,Ierapetritou M G,Kelly J D.Refinery short-term scheduling using continuous time formulation:Crude oil operations[J].Industrial & Engineering Chemistry Research,2003,42(13):3085-3097

    [5] 周祥,郭錦標(biāo),周涵,等.原油混輸調(diào)度的模型法研究[J].石油煉制與化工,2008,39(1):49-52

    [6] Li Wenkai,Hui Chiwai,Hua Ben,et al.Scheduling crude oil unloading,storage,and processing[J].Industrial & Engineering Chemistry Research.2002,41(26):6723-6734

    [7] Reddy P C P,Karimi I A,Srinivasan R.Novel solution approach for optimizing crude oil operations[J].AIChE Journal,2004,50(6):1177-1197

    [8] Meyer C A,F(xiàn)loudas C A.Global optimization of a combinatorially complex generalized pooling problem[J].AIChE Journal,2006,52(3):1027-1037

    [9] Karuppiah R,Grossmann I E.Global optimization for the synthesis of integrated water systems in chemical processes[J].Computers & Chemical Engineering.2006,30(4):650-673

    [10]周智菊,周祥,郭錦標(biāo),等.基于異步時(shí)間段的連續(xù)時(shí)間原油混輸調(diào)度模型[J].石油煉制與化工,2015,46(6):96-101

    STUDY ON PIECEWISE LINEAR RELAXATION FOR CRUDE SCHEDULING OPTIMIZATION

    Zhou Xiang,Zhou Zhiju,Lü Ning,Guo Jinbiao

    (SINOPECResearchInstituteofPetroleumProcessing,Beijing100083)

    The efficiency of solving crude scheduling optimization models is seriously affected by nonlinear constraints,which are produced by the property calculations of mixed crudes. In order to empower the model to cope with multiple crude pooling situations,a piecewise linear relaxation approach was employed to convert the nonlinear constraints into linear ones. Crude scheduling optimization of a real refinery was used as an example and the sulfur content calculations of mixed crudes were converted into linear forms. The effect of the piecewise linear relaxation approach was compared with another two frequently-used approaches. The results indicated that by use of the piecewise linear relaxation approach,the efficiency of solving the converted model met industrial application demand;while the sulfur contents of the processed crudes lied between 0.50% and 3.00%,the piecewise linear relaxation approach was more conducive to acquire optimal scheduling schemes;as the sulfur content difference of the processed crudes increased from 0.50% to 2.50%,the calculated sulfur contents of mixed crudes in charging tanks remained close to the actual value with the piecewise linear relaxation approach,the relative deviations were under 5% and therefore the industrial demand for model accuracy was satisfied.

    crude scheduling;optimization model;piecewise linear relaxation;sulfur content

    2016-02-15;修改稿收到日期:2016-04-19。

    周祥,博士,主要研究方向?yàn)闊捰瓦^程模型及煉油企業(yè)調(diào)度優(yōu)化。

    周祥,E-mail:zhoux.ripp@sinopec.com。

    猜你喜歡
    廠區(qū)分段線性
    哈軸集團(tuán)南北廠區(qū)
    哈爾濱軸承(2023年4期)2024-01-11 11:55:38
    漸近線性Klein-Gordon-Maxwell系統(tǒng)正解的存在性
    一類連續(xù)和不連續(xù)分段線性系統(tǒng)的周期解研究
    線性回歸方程的求解與應(yīng)用
    分段計(jì)算時(shí)間
    二階線性微分方程的解法
    承載厚重記憶的蓮花山老廠區(qū)
    西江月(2017年4期)2017-11-22 07:24:24
    3米2分段大力士“大”在哪兒?
    太空探索(2016年9期)2016-07-12 10:00:04
    小型水電站廠區(qū)用電設(shè)計(jì)
    關(guān)于年齡分段的描述
    如何舔出高潮| 亚洲最大成人手机在线| 久久99热6这里只有精品| 丝袜喷水一区| 国产av麻豆久久久久久久| 黄色配什么色好看| 亚洲av中文av极速乱| 色综合亚洲欧美另类图片| 3wmmmm亚洲av在线观看| 超碰av人人做人人爽久久| 成人漫画全彩无遮挡| 国产精品久久视频播放| 亚洲成人久久爱视频| 97热精品久久久久久| 亚洲专区国产一区二区| 国产成人freesex在线 | 国产成人一区二区在线| 91久久精品电影网| 哪里可以看免费的av片| 亚洲欧美中文字幕日韩二区| 美女被艹到高潮喷水动态| 中文亚洲av片在线观看爽| 国产不卡一卡二| 久久综合国产亚洲精品| 麻豆精品久久久久久蜜桃| 免费一级毛片在线播放高清视频| 国产免费一级a男人的天堂| 国产私拍福利视频在线观看| 国产精华一区二区三区| 男女那种视频在线观看| 色哟哟哟哟哟哟| 91午夜精品亚洲一区二区三区| 亚洲精品日韩在线中文字幕 | 在线a可以看的网站| 久久久久九九精品影院| 久久久久国内视频| 丝袜美腿在线中文| 色播亚洲综合网| 欧美成人一区二区免费高清观看| 亚洲精品国产av成人精品 | 真人做人爱边吃奶动态| 久久亚洲国产成人精品v| 亚洲色图av天堂| 12—13女人毛片做爰片一| 亚洲自拍偷在线| 亚洲综合色惰| 波多野结衣高清作品| 一个人免费在线观看电影| 丝袜美腿在线中文| 九九爱精品视频在线观看| 国产成人freesex在线 | 亚洲精品粉嫩美女一区| 精品国内亚洲2022精品成人| 欧美成人精品欧美一级黄| 蜜桃久久精品国产亚洲av| 老师上课跳d突然被开到最大视频| 啦啦啦韩国在线观看视频| 久久久久久久久久成人| 国产高清不卡午夜福利| 91久久精品电影网| 亚洲成人精品中文字幕电影| 久久久成人免费电影| 成年女人永久免费观看视频| 亚洲高清免费不卡视频| 啦啦啦啦在线视频资源| 欧美激情久久久久久爽电影| 3wmmmm亚洲av在线观看| 欧美在线一区亚洲| 真实男女啪啪啪动态图| 22中文网久久字幕| 久久综合国产亚洲精品| 又爽又黄无遮挡网站| 国产成人一区二区在线| 波野结衣二区三区在线| 国产精品免费一区二区三区在线| 亚洲最大成人手机在线| 免费人成在线观看视频色| 精品熟女少妇av免费看| 亚洲av免费在线观看| 久久九九热精品免费| 亚洲欧美日韩高清在线视频| 成人av一区二区三区在线看| 久久精品综合一区二区三区| 久久精品国产99精品国产亚洲性色| 伦精品一区二区三区| 国产精品福利在线免费观看| 搡老妇女老女人老熟妇| 97超级碰碰碰精品色视频在线观看| 亚洲人成网站在线观看播放| 亚洲精品在线观看二区| av视频在线观看入口| 精品99又大又爽又粗少妇毛片| 国产精品国产高清国产av| 久久精品夜夜夜夜夜久久蜜豆| 天天躁日日操中文字幕| 国产精品一区二区免费欧美| 日韩人妻高清精品专区| 日本爱情动作片www.在线观看 | 99久久久亚洲精品蜜臀av| 亚洲av成人av| 男人和女人高潮做爰伦理| 午夜激情欧美在线| 欧美成人精品欧美一级黄| 99热这里只有是精品在线观看| 在线免费十八禁| 老司机影院成人| 在线看三级毛片| 欧美另类亚洲清纯唯美| 国产精品电影一区二区三区| 白带黄色成豆腐渣| 亚洲精品亚洲一区二区| 国产精品国产三级国产av玫瑰| 深夜精品福利| 日韩人妻高清精品专区| 国产精品久久久久久精品电影| 一级毛片电影观看 | 天美传媒精品一区二区| 国产美女午夜福利| 亚洲在线自拍视频| 亚洲国产日韩欧美精品在线观看| 黄色一级大片看看| 国产成人精品久久久久久| 少妇被粗大猛烈的视频| 亚洲精品粉嫩美女一区| 亚洲婷婷狠狠爱综合网| 亚洲av成人av| 亚洲欧美清纯卡通| 午夜久久久久精精品| 内射极品少妇av片p| 日韩精品青青久久久久久| or卡值多少钱| 麻豆成人午夜福利视频| 久久精品国产99精品国产亚洲性色| 久久亚洲国产成人精品v| 寂寞人妻少妇视频99o| 一进一出好大好爽视频| 一个人免费在线观看电影| 1024手机看黄色片| 少妇的逼水好多| 国产女主播在线喷水免费视频网站 | 亚洲国产精品合色在线| 99热这里只有精品一区| 亚洲国产精品国产精品| 成人永久免费在线观看视频| 亚洲成人久久性| 久久人人精品亚洲av| 久久99热这里只有精品18| 成年女人永久免费观看视频| 一a级毛片在线观看| 国产三级中文精品| 成人毛片a级毛片在线播放| 国产免费一级a男人的天堂| 国产精品久久久久久精品电影| 97超视频在线观看视频| 久久精品国产清高在天天线| 18禁黄网站禁片免费观看直播| 国产av不卡久久| 我的女老师完整版在线观看| 精品午夜福利在线看| 卡戴珊不雅视频在线播放| 久久草成人影院| 可以在线观看毛片的网站| 欧美国产日韩亚洲一区| 乱系列少妇在线播放| 18禁在线无遮挡免费观看视频 | 大香蕉久久网| 两个人的视频大全免费| 亚洲精品日韩av片在线观看| 精品午夜福利在线看| 1024手机看黄色片| 国产片特级美女逼逼视频| 亚洲欧美日韩高清在线视频| 成人亚洲欧美一区二区av| 亚洲欧美日韩东京热| av在线蜜桃| 国产av不卡久久| 一个人免费在线观看电影| 国产69精品久久久久777片| 欧美色视频一区免费| 91久久精品国产一区二区成人| 久久99热6这里只有精品| 麻豆国产av国片精品| 级片在线观看| 国产在视频线在精品| 国产精品一区二区三区四区免费观看 | 免费高清视频大片| 麻豆久久精品国产亚洲av| 一进一出抽搐gif免费好疼| 老熟妇仑乱视频hdxx| 精品乱码久久久久久99久播| 久久久精品欧美日韩精品| 日日摸夜夜添夜夜爱| 精品免费久久久久久久清纯| 黄色视频,在线免费观看| 国内精品宾馆在线| 中文字幕免费在线视频6| 搡老熟女国产l中国老女人| av天堂中文字幕网| 久久中文看片网| 欧美性猛交╳xxx乱大交人| 国产视频内射| 国产免费男女视频| 无遮挡黄片免费观看| 成人漫画全彩无遮挡| 丝袜美腿在线中文| 精品久久久久久久末码| 国产91av在线免费观看| 欧美成人精品欧美一级黄| 日韩大尺度精品在线看网址| 亚洲电影在线观看av| 国产男靠女视频免费网站| 久久久久久久久大av| 97热精品久久久久久| 69人妻影院| 99在线视频只有这里精品首页| 色视频www国产| 欧美中文日本在线观看视频| 欧美日本视频| 又黄又爽又免费观看的视频| av天堂在线播放| 亚洲一级一片aⅴ在线观看| 91久久精品国产一区二区成人| 国产单亲对白刺激| 成人高潮视频无遮挡免费网站| 不卡一级毛片| 亚洲av五月六月丁香网| 中文亚洲av片在线观看爽| 久久久久久九九精品二区国产| 国产一区二区激情短视频| 91久久精品国产一区二区三区| 国产淫片久久久久久久久| 午夜福利视频1000在线观看| 观看免费一级毛片| 午夜日韩欧美国产| 久久精品国产亚洲av涩爱 | 日韩欧美 国产精品| 自拍偷自拍亚洲精品老妇| 久久久久久伊人网av| 丰满的人妻完整版| 禁无遮挡网站| 青春草视频在线免费观看| 欧美日韩一区二区视频在线观看视频在线 | 91在线精品国自产拍蜜月| 最近中文字幕高清免费大全6| 一区二区三区四区激情视频 | 给我免费播放毛片高清在线观看| 亚洲av成人精品一区久久| 成人av一区二区三区在线看| 听说在线观看完整版免费高清| 国产美女午夜福利| 成人精品一区二区免费| 1000部很黄的大片| 天堂动漫精品| 亚洲人与动物交配视频| 亚洲久久久久久中文字幕| 日本成人三级电影网站| 中文字幕熟女人妻在线| 精品久久久噜噜| 久久精品国产自在天天线| 51国产日韩欧美| av福利片在线观看| 精品免费久久久久久久清纯| 精品久久久久久久久久免费视频| 午夜a级毛片| 美女 人体艺术 gogo| 在线a可以看的网站| 亚洲av一区综合| 国产精品一区二区免费欧美| 男人舔女人下体高潮全视频| 国产淫片久久久久久久久| 白带黄色成豆腐渣| 久久久久久国产a免费观看| 久久久久久久久久久丰满| 人人妻人人澡欧美一区二区| av女优亚洲男人天堂| 免费无遮挡裸体视频| a级毛片免费高清观看在线播放| 一卡2卡三卡四卡精品乱码亚洲| 身体一侧抽搐| 久久6这里有精品| 大又大粗又爽又黄少妇毛片口| 久久精品夜夜夜夜夜久久蜜豆| 一级黄片播放器| 国产一级毛片七仙女欲春2| 99国产精品一区二区蜜桃av| 搞女人的毛片| a级毛片免费高清观看在线播放| 久久人妻av系列| 一级a爱片免费观看的视频| 亚洲国产色片| 欧美最新免费一区二区三区| 69av精品久久久久久| 日韩在线高清观看一区二区三区| 欧美精品国产亚洲| 男女那种视频在线观看| 亚洲综合色惰| 别揉我奶头 嗯啊视频| 久久人人精品亚洲av| 成人av在线播放网站| 人人妻人人澡人人爽人人夜夜 | 级片在线观看| 91av网一区二区| 国产美女午夜福利| 舔av片在线| 最近最新中文字幕大全电影3| 国产 一区精品| 可以在线观看的亚洲视频| 亚洲专区国产一区二区| 日日干狠狠操夜夜爽| 国产精品人妻久久久久久| 十八禁网站免费在线| 久久久久久久久久久丰满| 成人三级黄色视频| 桃色一区二区三区在线观看| 国产三级中文精品| 老司机影院成人| 国语自产精品视频在线第100页| 国产一区二区亚洲精品在线观看| 国产大屁股一区二区在线视频| 精品人妻视频免费看| 国产人妻一区二区三区在| 中文亚洲av片在线观看爽| 免费人成视频x8x8入口观看| 不卡视频在线观看欧美| 欧美日韩乱码在线| 两性午夜刺激爽爽歪歪视频在线观看| 久久久精品欧美日韩精品| 精品久久久噜噜| 亚洲av.av天堂| 色av中文字幕| 黄片wwwwww| 美女高潮的动态| ponron亚洲| 欧美xxxx黑人xx丫x性爽| 伊人久久精品亚洲午夜| 国产精品永久免费网站| 国产高清激情床上av| 免费看av在线观看网站| 此物有八面人人有两片| 可以在线观看的亚洲视频| 精品一区二区三区视频在线观看免费| 日本一二三区视频观看| av.在线天堂| 成人无遮挡网站| 成人高潮视频无遮挡免费网站| 嫩草影院精品99| 成人欧美大片| 国产精品日韩av在线免费观看| 99热这里只有是精品50| 亚洲18禁久久av| 老女人水多毛片| 精品国产三级普通话版| 在线播放国产精品三级| 久久99热这里只有精品18| 女的被弄到高潮叫床怎么办| 18禁黄网站禁片免费观看直播| 欧美色欧美亚洲另类二区| 在线观看av片永久免费下载| 亚洲人与动物交配视频| 亚洲精品亚洲一区二区| 三级毛片av免费| 色哟哟哟哟哟哟| 国产亚洲av嫩草精品影院| 中文亚洲av片在线观看爽| 国产成人a∨麻豆精品| 一个人看视频在线观看www免费| 18禁在线播放成人免费| 国产午夜精品论理片| a级毛片免费高清观看在线播放| 99精品在免费线老司机午夜| 久久人人爽人人爽人人片va| 一级毛片aaaaaa免费看小| 丰满人妻一区二区三区视频av| 亚洲图色成人| 男女下面进入的视频免费午夜| 91精品国产九色| 少妇人妻一区二区三区视频| 性欧美人与动物交配| 国产三级在线视频| 97超级碰碰碰精品色视频在线观看| 精品人妻偷拍中文字幕| 中文字幕精品亚洲无线码一区| 12—13女人毛片做爰片一| 日本欧美国产在线视频| 两个人的视频大全免费| 99riav亚洲国产免费| 白带黄色成豆腐渣| 免费搜索国产男女视频| 欧美日韩综合久久久久久| 天天躁日日操中文字幕| 国产白丝娇喘喷水9色精品| 亚洲aⅴ乱码一区二区在线播放| 国产精品99久久久久久久久| 成人漫画全彩无遮挡| av在线亚洲专区| 少妇裸体淫交视频免费看高清| 日日摸夜夜添夜夜添av毛片| 免费大片18禁| 又黄又爽又免费观看的视频| 亚洲国产精品成人综合色| 久久鲁丝午夜福利片| 91在线精品国自产拍蜜月| 国产探花极品一区二区| 亚洲国产精品成人综合色| 久久久久国产精品人妻aⅴ院| 成人午夜高清在线视频| 白带黄色成豆腐渣| 亚洲国产精品成人综合色| 老熟妇乱子伦视频在线观看| 日韩av在线大香蕉| 欧美不卡视频在线免费观看| 午夜福利在线观看吧| 三级男女做爰猛烈吃奶摸视频| 六月丁香七月| 你懂的网址亚洲精品在线观看 | 综合色av麻豆| 伦精品一区二区三区| 男女之事视频高清在线观看| 久久久久国产网址| 99久久无色码亚洲精品果冻| 欧美最新免费一区二区三区| 日本一二三区视频观看| 一进一出抽搐gif免费好疼| av女优亚洲男人天堂| 国产精品一区二区性色av| 别揉我奶头~嗯~啊~动态视频| 三级国产精品欧美在线观看| 51国产日韩欧美| 久久精品人妻少妇| 国产精品一二三区在线看| 日日啪夜夜撸| 精品午夜福利视频在线观看一区| 久久这里只有精品中国| 国产高清三级在线| 男女边吃奶边做爰视频| 亚洲成av人片在线播放无| 免费av不卡在线播放| 美女内射精品一级片tv| 你懂的网址亚洲精品在线观看 | 在线观看午夜福利视频| 99久久久亚洲精品蜜臀av| 中文字幕精品亚洲无线码一区| 无遮挡黄片免费观看| 欧美高清性xxxxhd video| 99热网站在线观看| a级一级毛片免费在线观看| 美女xxoo啪啪120秒动态图| 国产毛片a区久久久久| 插阴视频在线观看视频| 老熟妇仑乱视频hdxx| 国产亚洲欧美98| 国内揄拍国产精品人妻在线| 久久久精品大字幕| 俄罗斯特黄特色一大片| 亚洲色图av天堂| 一级毛片我不卡| 一区二区三区高清视频在线| 欧美中文日本在线观看视频| 日日撸夜夜添| 中文字幕av在线有码专区| 最近中文字幕高清免费大全6| 一级黄色大片毛片| 日韩 亚洲 欧美在线| 日韩欧美精品v在线| 色尼玛亚洲综合影院| 国产欧美日韩精品亚洲av| 欧美日本视频| 午夜福利在线观看吧| 悠悠久久av| 国产亚洲欧美98| 热99在线观看视频| 99精品在免费线老司机午夜| 男女之事视频高清在线观看| 你懂的网址亚洲精品在线观看 | 免费观看的影片在线观看| 亚洲精品亚洲一区二区| 成年av动漫网址| 十八禁网站免费在线| 尾随美女入室| 久久久久国内视频| 亚洲av中文av极速乱| 成熟少妇高潮喷水视频| 天堂动漫精品| 欧美激情国产日韩精品一区| 久久久久九九精品影院| av天堂在线播放| 欧洲精品卡2卡3卡4卡5卡区| 精品久久久久久久久av| 久久久久国产网址| 在现免费观看毛片| 欧美日韩在线观看h| 久久久久久久久久久丰满| 床上黄色一级片| 国产探花极品一区二区| 国产亚洲精品久久久com| 99国产极品粉嫩在线观看| 亚洲久久久久久中文字幕| 精品午夜福利在线看| av在线亚洲专区| 亚洲精品久久国产高清桃花| 亚洲精品日韩在线中文字幕 | av天堂中文字幕网| 亚洲人与动物交配视频| eeuss影院久久| 精品久久久噜噜| 日韩欧美免费精品| 日本黄色视频三级网站网址| 成年免费大片在线观看| 久久中文看片网| 中文资源天堂在线| 久久久久国内视频| av在线观看视频网站免费| 亚洲国产精品国产精品| 露出奶头的视频| 免费av不卡在线播放| 国产黄色视频一区二区在线观看 | 自拍偷自拍亚洲精品老妇| 在线播放国产精品三级| 露出奶头的视频| .国产精品久久| 51国产日韩欧美| 亚洲美女视频黄频| 久久久久久久久中文| 亚洲高清免费不卡视频| 欧美+亚洲+日韩+国产| 国产久久久一区二区三区| 成年女人永久免费观看视频| 看十八女毛片水多多多| 最近2019中文字幕mv第一页| 国产伦精品一区二区三区视频9| 天堂av国产一区二区熟女人妻| 日韩av在线大香蕉| 日韩一本色道免费dvd| 高清毛片免费看| 欧美一区二区国产精品久久精品| 国语自产精品视频在线第100页| 午夜爱爱视频在线播放| 国产高清激情床上av| 成人精品一区二区免费| 欧美xxxx黑人xx丫x性爽| 波多野结衣高清无吗| 精品久久久久久久久久久久久| 噜噜噜噜噜久久久久久91| 在线a可以看的网站| 国产日本99.免费观看| 色吧在线观看| 内射极品少妇av片p| 日韩亚洲欧美综合| 少妇丰满av| 日韩一区二区视频免费看| 国产精品久久久久久av不卡| 此物有八面人人有两片| 午夜久久久久精精品| 国产精品女同一区二区软件| 亚洲欧美精品自产自拍| 色av中文字幕| 欧美性猛交黑人性爽| 亚洲成人av在线免费| 麻豆精品久久久久久蜜桃| 久久久久久久久久黄片| 亚洲欧美日韩东京热| 简卡轻食公司| 日韩欧美一区二区三区在线观看| 亚洲国产日韩欧美精品在线观看| 国产亚洲av嫩草精品影院| 精华霜和精华液先用哪个| 欧美潮喷喷水| 欧美成人一区二区免费高清观看| 日本与韩国留学比较| 精品久久久噜噜| 丝袜美腿在线中文| 午夜福利在线在线| 亚洲av成人精品一区久久| 国产高清有码在线观看视频| 日韩人妻高清精品专区| 级片在线观看| 综合色丁香网| 国产精品综合久久久久久久免费| 国产精品永久免费网站| 一级毛片久久久久久久久女| 1000部很黄的大片| 亚洲图色成人| 超碰av人人做人人爽久久| 悠悠久久av| 欧美国产日韩亚洲一区| 免费大片18禁| 欧美一区二区精品小视频在线| 亚洲国产精品国产精品| 少妇裸体淫交视频免费看高清| 女同久久另类99精品国产91| 乱系列少妇在线播放| 午夜福利18| 精品人妻一区二区三区麻豆 | 男女视频在线观看网站免费| 亚洲国产精品合色在线| 在现免费观看毛片| 午夜福利18| av在线老鸭窝| 久久综合国产亚洲精品| 亚洲中文字幕日韩| 国产精品一二三区在线看| www.色视频.com| 午夜免费男女啪啪视频观看 | 三级经典国产精品| 乱系列少妇在线播放| 看非洲黑人一级黄片| 亚洲在线自拍视频| 亚洲性夜色夜夜综合| 亚洲专区国产一区二区| 久99久视频精品免费| www日本黄色视频网| 国产精品久久电影中文字幕| 国产国拍精品亚洲av在线观看| 特级一级黄色大片| 禁无遮挡网站|