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

    基于智能預測和機理模型的換熱網絡清洗決策

    2022-04-26 09:47:50蔣寧張元毅范偉趙世超徐新杰徐英杰
    化工進展 2022年4期
    關鍵詞:公用結垢換熱器

    蔣寧,張元毅,范偉,趙世超,徐新杰,徐英杰

    (浙江工業(yè)大學機械工程學院,浙江 杭州 310023)

    隨著工業(yè)大數(shù)據(jù)和工業(yè)物聯(lián)網的發(fā)展,對過程工業(yè)建立智能化監(jiān)控平臺可以進一步提高生產安全和生產效率?,F(xiàn)役的石化裝置中布置了較多的傳感器,但僅實現(xiàn)關鍵指標的閾值報警,而不能實現(xiàn)對生產運行性能的預測及優(yōu)化。換熱網絡是過程工業(yè)系統(tǒng)中重要的能量回收單元,提高換熱網絡的運行性能可以有效減少公用工程能耗。目前,針對換熱網絡的研究,主要是通過換熱網絡的優(yōu)化設計和優(yōu)化改造方法,提高換熱網絡的能量回收能力。對于石化裝置的換熱網絡而言,換熱器長時間處在高溫高鹽的運行條件,導致?lián)Q熱器性能衰退,進一步影響換熱網絡的運行性能。

    目前提出的換熱器結垢模型,主要是以Kern提出的污垢增長的微分方程為基礎,綜合考慮其他不同因素對換熱器結垢的沉積率和剝蝕率的影響。Ebert 和Panchal提出的閾值結垢模型認為化學反應可以加快污垢生長,剪切效應可以減緩污垢生長,當化學反應和剪切效應的作用處于平衡時,此時污垢生長緩慢;Knudsen 等和Polley 等對該模型進行了改進,獲得了更好的預測精度。隨著智能算法的發(fā)展,可以進一步挖掘數(shù)據(jù)隱藏的關系,更多地關注運行條件對換熱器性能的影響。樊紹勝基于影響污垢變化的因素有部分關系已知的特點,提出一種基于灰色理論的污垢預測方法。徐志明等采用改進的灰色理論的污垢預測方法,綜合考慮未來的一些擾動因素對污垢預測的影響,通過在建模過程中及時去掉數(shù)據(jù)集中的舊數(shù)據(jù)補充新數(shù)據(jù),實現(xiàn)換熱器性能的中長期預測。Sun 等采用遺傳算法優(yōu)化后的支持向量機模型,優(yōu)化后的訓練精度要明顯優(yōu)于標準的支持向量機模型。神經網絡算法(ANN)具有較強的學習能力和構建非線性的復雜關系的優(yōu)點,被廣泛用于工業(yè)數(shù)據(jù)的建模。Aminian 和Shahhosseini采用前饋網絡和反向傳播算法相結合的方法對神經網絡模型進行訓練,通過試錯法找出隱藏層數(shù)和神經元數(shù)的最佳組合,與實驗數(shù)據(jù)相比,ANN模型的平均相對誤差為26.23%,ANN 模型對原油污垢熱阻變化率的預測精度明顯優(yōu)于閾值結垢模型。Sundar等提出以深度學習為基礎的預測模型,對錯流式換熱器中的污垢形成進行預測分析?;谠撋疃葘W習模型的結果表明,污垢熱阻絕對平均預測誤差低于10kW,系數(shù)超過99%。

    換熱器結垢是一個長期存在的問題,而由結垢所造成的換熱器換熱效果衰退,一方面需要通過更嚴格的操作條件來減緩換熱器結垢,同時需要增加冷熱物流末端的公用工程供應以滿足工藝生產要求。通過對結垢的換熱器制定定期清洗方案,可以重新恢復換熱器的傳熱性能。Sma?li 等提出,污垢可被大致分為兩層,即新鮮層和老化層,老化的程度決定了沉積物的狀態(tài),也決定了污垢層去除的難易程度。針對這兩層,操作員可以選擇不同的清潔方式,清洗方式可分為溶劑清洗和機械清洗。溶劑清洗無法徹底清除沉積物,而在污垢沉積物老化超過一定程度之前,機械清洗是不經濟的?;趦蓪映练e模型優(yōu)化換熱器清洗周期,Ishiyama 等提出沉積模型和經濟模型,也提出了超清洗周期的概念,即一個清洗周期包括數(shù)個子循環(huán),每個子循環(huán)的長度由給出的經濟模型進行評估。Pogiatzis等使用兩種方法來生成混合清潔調度問題的解,第一個基于啟發(fā)式行進算法,選定下一個子周期的清潔方法,最終選擇機械清洗;第二種以每一個子周期長度為變量,選定溶劑清洗的數(shù)量,模擬3種情況,即基本情況、沉積速率變化、老化速率變化,找到每種情況最佳溶劑清洗次數(shù)和一個清洗周期的天數(shù)。但上述方法僅僅局限于對單臺換熱器清洗方案的研究,無法從系統(tǒng)集成的角度分析換熱網絡整體性能的變化,來制定清洗方案。

    目前對換熱網絡的清洗決策主要是建立相應的混合整數(shù)非線性規(guī)劃(MINLP)模型,通過不同算法進行求解以獲得經濟性的清洗方案。Pogiatzis等提出了一個MINLP模型,研究了熱交換器的清洗調度問題;基于此,Georgiadis 等建立新的混合整數(shù)線性規(guī)劃(MILP)模型,該模型沒有考慮有關傳熱或結垢模型的非線性模型,對小型換熱網絡清洗問題,該模型能夠實現(xiàn)全局最優(yōu)性。樊婕等提出以換熱器清洗的最大允許污垢熱阻為優(yōu)化變量,取代表示換熱器是否清洗的二進制變量,將混合整數(shù)非線性規(guī)劃問題(MINLP)轉化成非線性規(guī)劃問題(NLP),降低求解難度。這種求解模型對石化現(xiàn)場換熱網絡的性能監(jiān)測比較少,同時為了獲得最佳的經濟性方案,而忽略了石化現(xiàn)場清洗換熱網絡其他因素的影響。Biyanto等采用帝國競爭算法(ICA)求解,通過該方法求解得到的清洗計劃相比原來節(jié)省了約44.83%的傳熱損失,在44個月內節(jié)省了105 萬美元,占結垢損失費用的22.12%,具有較好的經濟性。Xiao等提出了一種遺傳/模擬退火算法(GA/SA)同時優(yōu)化柔性換熱網絡和清洗計劃的新方法。與以往的MINLP 模型不同,通過初次綜合策略可以有效地消除同時多周期MINLP 模型中的大量的二元變量,進一步優(yōu)化了換熱網絡中換熱器的匹配關系及換熱器傳熱面積和運行周期。Diaby 等采用遺傳算法來求解在結垢和不同老化情況下,原油煉廠預熱機組最佳清洗時間表的問題。與未優(yōu)化清洗和不進行清洗時的原油煉廠預熱機組相比,通過遺傳算法優(yōu)化得到的最佳清洗計劃,可以實現(xiàn)較大的經濟效益。

    與單臺換熱器性能預測模型不同的是,換熱網絡由多臺換熱器組成,不同換熱器性能變化對換熱網絡整體性能變化的影響不同,換熱網絡整體性能變化要比單臺換熱器性能變化更顯著。目前對換熱網絡的清洗決策主要是建立相應的MINLP 模型,通過不同算法進行求解以獲得經濟性的清洗方案。這種求解模型對石化現(xiàn)場換熱網絡的性能監(jiān)測比較少,同時為了獲得最佳的經濟性方案,忽略了石化現(xiàn)場清洗換熱網絡時一些其他因素,包括不同清洗準則的制定、公用工程成本的定價、停用清洗時間和換熱器清洗成本的定價等的影響,這主要影響了最佳清洗的選擇及經濟性結果,而且也沒有考慮清洗時無備用換熱器的影響,在換熱器停用清洗期間換熱網絡需要更多的公用工程來滿足生產工藝要求,需要考慮停用清洗所帶來的損失費用和清洗成本。因此,本文提出了一種基于智能預測和機理模型的換熱網絡清洗決策方法,通過將換熱器智能預測模型與換熱網絡性能模擬模型相結合,獲得換熱網絡的性能變化趨勢,從而從網絡整體性能變化角度來制定清洗方案,實現(xiàn)換熱網絡的高效運行。與原有的方法相比,該方法基于現(xiàn)役運行數(shù)據(jù)建立換熱器性能智能預測模型,同時對換熱網絡整體性能進行監(jiān)測,更能反映現(xiàn)役換熱網絡的性能狀態(tài)。

    1 基于數(shù)據(jù)的換熱器性能智能預測

    原油裝置的換熱器長時間處在高鹽高溫的運行狀態(tài),導致?lián)Q熱器結垢、性能衰退嚴重,傳統(tǒng)的時間序列預測模型和污垢閾值模型對換熱器性能預測差異較大,因此本文基于神經網絡算法建立換熱器性能智能預測模型,以獲得換熱器的性能變化趨勢。

    1.1 原油換熱器閾值結垢模型

    傳統(tǒng)污垢熱阻的計算公式中參數(shù)較多,測量不方便,同時計算誤差較大,對于給定的流體組成和換熱表面狀況,換熱設備換熱面的溫度和流速是影響設備中原油污垢生長的重要因素。因此,針對流速和溫度兩個重要因素,學者提出了有關溫度和流速的原油換熱設備結垢的經驗公式,主要包括以下4種閾值模型。

    Saleh 等提出的模型只考慮了溫度的影響,未考慮流速對污垢形成速率的影響,該模型主要運用在過渡區(qū)域,見式(1)。

    1995 年,Ebert 和Panchal提 出 了 原 油 污 垢 生長的閾值模型,該模型表明污垢的生長不僅取決于化學反應和溫度,還與管壁的剪切應力有關,見式(2)。

    2002 年,Polley 等對原油結垢的閾值模型進行了改進,模型中引入了普朗特數(shù),見式(3)。

    2006 年,Jafari 和Mejidi提 出 了 一 種 新 的 模型,見式(4)。

    表1 的實驗數(shù)據(jù)是某原油換熱器在60h 內的11 組實驗數(shù)據(jù),包括時間、流速、雷諾數(shù)、溫度和污垢熱阻變化率。

    表1 各維度數(shù)據(jù)分布指標表

    該實驗的污垢熱阻變化曲線如圖1所示。通過擬合上述閾值模型的經驗公式,發(fā)現(xiàn)無法通過以往的經驗公式較好地擬合該臺換熱器的閾值結垢曲線。因此,根據(jù)閾值結垢模型的一般方程表達式,擬合得到式(5),擬合公式的為0.07886,為0.0018,絕對平均誤差為0.0005m·K/(kW·h),相對誤差率為77%。從表1中可見,該臺換熱器的流速變化區(qū)間小,導致不同時刻雷諾數(shù)()基本相同,以致傳統(tǒng)閾值模型中的關鍵參數(shù)成為一個固定參數(shù),使得針對該實驗的閾值結垢模型不能反映流速的影響,導致經驗公式在該實驗數(shù)據(jù)上的失效。根據(jù)以往的研究可知,閾值結垢模型反映的是溫度和流速之間的制約和平衡關系,往往實驗中需要較大范圍的溫度和流速的實驗,以探究不同溫度和不同流速之間的關系;而對于較小范圍的溫度和流速的實驗,由于溫度范圍或流速范圍較小,無法發(fā)現(xiàn)兩者之間的內涵關系,傳統(tǒng)理論研究的經驗模型無法適用。

    圖1 經驗模型的擬合結果

    1.2 神經網絡預測模型建立

    隨著工業(yè)數(shù)據(jù)規(guī)模的擴大,支持向量回歸(SVR)、自適應增強算法(Adaboost)等傳統(tǒng)的機器學習對大數(shù)據(jù)模型很難有較好的訓練效果,而神經網絡等深度學習算法的發(fā)展,具有廣泛的適用性,可以充分挖掘數(shù)據(jù)的內在關系,可以對換熱器運行數(shù)據(jù)建立更好的智能預測模型,神經網絡模型的建立步驟如下。

    (1)收集換熱器結垢數(shù)據(jù),該數(shù)據(jù)集由時間、黏度、壓力、流速、氧氣含量、熱負荷、管壁溫度、管腔溫度和污垢熱阻9個維度數(shù)據(jù)組成。

    (2)根據(jù)歸一化公式,對該數(shù)據(jù)集各維度數(shù)據(jù)采用線性函數(shù)歸一化處理,采用歸一化處理的數(shù)據(jù)集,各維度數(shù)據(jù)都分布在[0,1]的范圍內。

    (3)選擇神經網絡神經網絡層數(shù),輸入層為時間、黏度、壓力、流速、氧氣含量、熱負荷、管壁溫度、管腔溫度8 個結點,輸出層節(jié)點數(shù)為1 個(污垢熱阻)。選擇恰當?shù)碾[藏層層數(shù),神經元數(shù)和激活函數(shù)。

    (4)采用不同的優(yōu)化器訓練BP 模型,直到模型達到理想的預測性能,從而建立BP 神經網絡預測模型。采用不同的回歸模型評估指標平均方差(MSE)、均方根誤差(RMSE)和差等反映回歸模型的預測結果的優(yōu)劣。

    本文通過對文獻[25]中的原油換熱器的污垢數(shù)據(jù)采用神經網絡預測模型,經參數(shù)訓練后,得到最佳的參數(shù)組合為神經元數(shù)為100,隱藏層數(shù)為1層,激活函數(shù)為linear,優(yōu)化器為adam方法,訓練樣本為600組,與SVR模型和Adaboost模型訓練的結果曲線如圖2所示。從圖2中可以看出,采用神經網絡(BP)模型和SVR 模型的預測結果明顯優(yōu)于Adaboost 模型。BP 模型和SVR 模型在訓練集和測試集上都逼近原始數(shù)據(jù)曲線,訓練效果較好,而Adaboost模型在訓練集和測試集均無法逼近原始數(shù)據(jù)曲線,導致訓練效果較差。不同算法的污垢預測模型的評價指標結果如圖3所示。從圖3中可以看出,采用BP 模型和SVR 模型在MAE、RMSE 和三個指標上相近,采用Adaboost 模型在MAE、RMSE和三個指標上要明顯劣于神經網絡模型和SVR模型,神經網絡模型在此三個指標上反映出的結果最優(yōu)。因此,小數(shù)據(jù)集上,BP 模型和SVR 模型都具有較好的預測效果,Adaboost模型預測結果較差。隨著工業(yè)數(shù)據(jù)規(guī)模的增加,神經網絡算法可提供更好的優(yōu)化器算法,能夠快速地訓練得到最佳的參數(shù)組合。傳統(tǒng)的機器學習方法如SVR 算法,當數(shù)據(jù)量較大時,需要求解的高階矩陣存儲和計算將耗費大量的機器內存和運行時間,導致復雜度增加,不適合大規(guī)模數(shù)據(jù)的建模。隨著不同類型算法的提出,為不同結垢機理的換熱器運行數(shù)據(jù)建立智能預測模型提供了更多選擇。

    圖2 與傳統(tǒng)機器學習模型結果比較

    圖3 不同預測模型評價指標結果對比

    2 換熱網絡的性能模擬模型

    本節(jié)建立基于數(shù)據(jù)和機理的換熱網絡性能預測方法,實現(xiàn)對現(xiàn)役換熱網絡的性能預測。

    換熱網絡是由數(shù)量較多的換熱器,通過換熱器之間的串聯(lián)和并聯(lián),實現(xiàn)冷熱流體的熱量交換。因此,換熱網絡性能的模擬可以先模擬單個換熱器的性能,再根據(jù)換熱網絡中換熱器之間上下游和串并聯(lián)關系求解獲得。

    在以往的研究中對單個換熱器性能的模擬,通常假設換熱器沒有發(fā)生相變,基于能量守恒定律和換熱器傳熱方程,此時換熱器性能模擬的具體計算見式(6)~式(8)。

    其中,LMTD 為對數(shù)平均溫差,計算公式見式(9)。

    為簡化計算,設參數(shù),表達式見式(10)。

    當≠1時,有式(11)。

    依據(jù)對以及的定義,可將公式轉換為式(12)、式(13)。

    當=1時,公式可以轉換為式(14)~式(16)。

    由此可見,在已知換熱器的換熱面積、冷熱物流通過此換熱器的進口溫度以及熱容流率、傳熱系數(shù)的信息后,可求得該換熱器的或以及值,由此構成換熱網絡性能模擬非線性方程組的、、系數(shù)矩陣,則式(12)、式(13)或式(15)、式(16)這兩個方程僅有兩個未知量(熱物流的出口溫度和冷物流的出口溫度),對于一個由臺換熱器組成的換熱網絡而言,則共有2個方程,通過求解這些方程,可以得到換熱網絡的溫度分布。最后根據(jù)溫度分布與物流目標溫度,可求得換熱網絡公用工程的消耗量。

    3 換熱網絡性能預測方法

    3.1 方法介紹

    隨著工業(yè)物聯(lián)網的發(fā)展,石化公司對現(xiàn)有的石化裝置布置了多種傳感器。雖然石化現(xiàn)場能夠采集到換熱器的溫度、流速、壓力等工況運行數(shù)據(jù),但是對采集數(shù)據(jù)的應用僅停留在簡單的計算和閾值報警,因此失去了讓數(shù)據(jù)促進生產節(jié)能增效的作用。通過對采集的換熱器運行數(shù)據(jù)建立神經網絡預測模型,可以進一步改進對換熱器運行性能的監(jiān)測,實現(xiàn)石化生產的智能化。

    由于換熱網絡中換熱器之間存在上下游和串并聯(lián)的關系,導致?lián)Q熱網絡中換熱器之間相互耦合,從而導致不同換熱器的性能衰退對下游換熱器性能和換熱網絡整體性能的影響不同。單純通過換熱網絡的進出口溫度數(shù)據(jù)建立神經網絡性能預測模型來實現(xiàn)對網絡性能預測,模型的預測精度有限,同時不能反映換熱網絡中不同換熱器的性能衰退對網絡性能的影響,因此將單臺換熱器智能預測模型與換熱網絡性能模擬模型相結合,實現(xiàn)對換熱網絡性能的預測。

    目前石化現(xiàn)場采用的清洗方案,主要是根據(jù)換熱器熱性能衰退到一定程度時,對該臺換熱器進行停用清洗,該清洗方案對換熱器的干預較少,換熱器性能可能長時期處在較差的運行狀態(tài),從而導致?lián)Q熱網絡需要較多額外的公用工程能耗供給以滿足生產工藝要求。因此本文提出基于換熱網絡整體性能來制定清洗方案的方法,當多臺換熱器同時發(fā)生性能衰退時,換熱網絡的整體性能變化要比單臺換熱器性能變化更顯著,因此從換熱網絡整體性能變化的角度制定清洗方案,可以提前對換熱器運行性能進行清洗干預,使得換熱網絡整體性能保持在較好的狀態(tài)。

    3.2 求解流程

    基于智能預測和機理的換熱網絡性能預測方法的求解流程如圖4所示,具體實現(xiàn)步驟如下。

    圖4 本文方法的流程圖

    (1)確定網絡拓撲結構

    根據(jù)現(xiàn)役的換熱網絡結構圖,確定換熱網絡中換熱器的拓撲結構關系,將網絡結構信息編寫為可識別的矩陣信息。

    (2)輸入傳熱系數(shù)矩陣

    通過換熱器運行工況數(shù)據(jù)的采集、數(shù)據(jù)歸一化處理、預測算法的選擇和訓練預測,獲得換熱器運行周期內預測的傳熱系數(shù)值,建立換熱器智能預測模型,以滿足模型預測誤差精度要求。將不同換熱器的傳熱系數(shù)數(shù)據(jù)構建換熱網絡的傳熱系數(shù)矩陣,傳遞給換熱網絡性能模擬模型。

    (3)換熱網絡物流信息和換熱器信息

    收集換熱網絡的物流信息,包括冷熱物流的進出口溫度、熱容流率、質量流量等;收集網絡中換熱器的信息,包括換熱器的面積、傳熱系數(shù)等。

    (4)獲得、、系數(shù)矩陣

    對步驟(1)和步驟(3)得到的換熱器參數(shù)信息,獲得換熱網絡性能模擬非線性方程組的、、系數(shù)矩陣。

    (5)構建性能模擬線性方程組

    對步驟(1)~步驟(4)得到的換熱網絡結構與參數(shù)信息,建立換熱網絡的性能模擬模型,給出該網絡性能模擬的線性方程組。

    (6)采用高斯消去法求解

    對步驟(5)得到的線性方程組,采用高斯消去方法進行求解,獲得所建立的線性方程組的解集。

    (7)計算每天換熱網絡的性能指標

    對步驟(6)求解得到線性方程組的解,根據(jù)步驟(1)換熱網絡拓撲結構信息,確定換熱網絡的溫度分布,從而計算每天公用工程能耗和運行成本的變化情況。

    本文在MATLAB2018b 環(huán)境下對上述方法進行了具體實現(xiàn)。在此基礎上,進一步研究不同清洗方案對換熱網絡的影響,開展換熱網絡的清洗決策。

    4 案例研究與分析

    該案例為一個原油精餾單元的換熱網絡,該換熱網絡包括7條熱物流和3條冷物流,物流數(shù)據(jù)如表2 所示,網絡結構如圖5 所示,圖5 中紅色線條代表熱物流,藍色線條代表冷物流,流股上的數(shù)字是流股的進出口溫度,流股末端的斜體數(shù)字代表物流熱容流率,位于中間換熱的圈內數(shù)字代表中間換熱器,位于流股末端的圈內字母代表公用換熱器。該換熱網絡的換熱器的總傳熱系數(shù)都為0.4kW/(m·℃),為換熱器使用前總傳熱系數(shù)的初始值,換熱網絡公用工程總能耗為166000kW,年度公用工程費用6330000USD/a, 年操作時間8400h/a。

    圖5 現(xiàn)有換熱網絡結構圖

    表2 物流數(shù)據(jù)表

    4.1 換熱網絡性能預測

    根據(jù)該換熱網絡的基本數(shù)據(jù),靈敏度分析的計算結果如圖6所示,該換熱網絡中HE5是最重要的設備,HE5性能衰退所導致的換熱網絡公用工程能耗增加最多;HE1 和HE2 是次重要設備,這兩臺換熱器的性能衰退也會引起換熱網絡公用工程能耗的增加;HE6的性能衰退對換熱網絡的公用能耗增加影響較小,這主要是因為該臺換熱器換熱面積較小,換熱能力較弱,同時處于冷物流起始段第一個換熱器,因此引起的冷物流終端溫度變化較小,進而所增加的公用工程能耗較小。

    圖6 換熱網絡關鍵設備靈敏度分析

    為進一步反應換熱網絡的整體性能隨換熱器性能衰退的變化,考慮關鍵的HE1、HE2和HE5傳熱系數(shù)年度的變化情況,其中HE1 和HE2 的傳熱系數(shù)變化曲線如圖7所示,HE5的傳熱系數(shù)變化曲線如圖8 所示,經過60 天后,HE1 和HE2 的傳熱系數(shù)降為初值的74.63%,HE5 的傳熱系數(shù)降為初值的67.43%;HE1和HE2的傳熱系數(shù)在275天后趨于平緩,為初值的58.69%,HE5的傳熱系數(shù)在200天以后趨于平緩,為初值的58.84%。

    圖7 HE1和HE2換熱器傳熱系數(shù)時間變化曲線

    圖8 HE5換熱器傳熱系數(shù)變化曲線

    經過換熱網絡性能模擬模型計算得到的熱交換網絡公用能耗絕對增加量年度的變化曲線如圖9所示。位于C2 物流上的HE1 和HE2,雖然HE1 和HE2具有相同的傳熱系數(shù)變化趨勢,根據(jù)靈敏度分析可知,相同的性能衰退條件下,HE2所導致的熱交換網絡性能下降更嚴重,以致需要更多的公用工程能耗。位于C1 物流上的HE5,一方面其傳熱性能不斷衰退,另外作為網絡中最關鍵的設備,隨著其性能的下降,導致在換熱網絡整個運行周期內,需要較多的公用工程能耗。若同時考慮HE1、HE2和HE5 性能發(fā)生衰退,因為HE1 和HE2 都位于冷物流C2 上,雖然兩臺換熱器都以相同的性能衰退曲線,但是兩臺換熱器的性能變化對整個網絡的性能影響不同,HE2的性能衰退對網絡整體的影響更大。研究發(fā)現(xiàn),當兩臺換熱器都發(fā)生性能衰退時,HE1 和HE2 兩臺換熱器對熱交換網絡的性能影響并不是HE1 和HE2 兩臺換熱器單獨作用的之和,這主要是因為HE1 和HE2 換熱器在熱交換網絡中存在著上下游耦合的關系。當HE1、HE2和HE5三臺換熱器同時發(fā)生性能衰退時,導致的熱交換網絡的性能衰退所需要的額外公用工程能耗是HE1 和HE2作用于C2物流和HE5作用C3物流兩個部分分別引起的額外增加的公用工程能耗之和,這主要是因為HE5和HE1或HE2之間不存在耦合關系。

    圖9 公用能耗絕對增加量隨時間的變化

    根據(jù)圖9的公用工程能耗絕對增加量隨時間的變化曲線,通過計算可以得到日公用工程能耗絕對增加量隨時間的變化,如圖10 所示。熱交換網絡經運行一年后,當HE1 性能發(fā)生衰退時,日額外公用工程能耗最高為4.2×10kJ/d;當HE2性能發(fā)生衰退時,日額外公用工程能耗最高為7.5×10kJ/d;當HE5 性能發(fā)生衰退時,日額外公用工程能耗最高為1.3×10kJ/d;當HE1 和HE2 同時發(fā)生性能衰退時,日額外公用工程能耗最高為1.2×10kJ/d;當HE1、HE2和HE5都發(fā)生性能衰退時,日額外公用工程能耗最高為2.5×10kJ/d,由此可知當多臺換熱器同時發(fā)生性能衰退時,換熱網絡的整體性能變化要比單臺換熱器性能變化更顯著。

    圖10 日公用能耗絕對增加量隨時間的變化

    根據(jù)圖10 日公用工程能耗絕對增加量隨時間的變化,可計算得到年度總的額外公用工程能耗,如圖11 所示。當HE1 性能發(fā)生衰退時,年度額外公用工程能耗增加量為1.09×10kJ/a;當HE2 性能發(fā)生衰退時,年度額外公用工程能耗2.07×10kJ/a;當HE5 性能發(fā)生衰退時,年度額外公用工程能耗3.87×10kJ/a;當HE1 和HE2 同時發(fā)生性能衰退時,年度額外公用工程能耗3.25×10kJ/a;當HE1、HE2 和HE5 都發(fā)生性能衰退時,年度額外公用工程能耗7.12×10kJ/a,其中年度額外熱公用工程能耗和年度冷公用工程能耗為3.55×10kJ/a,年度熱公用工程年費用增加704365USD,年度冷公用工程費用增加58697USD,年度公用工程能耗總費用增加763062USD,當HE1、HE2和HE5三臺換熱器都發(fā)生性能衰退時,年度額外需要公用工程能耗相對增加高達14.1%,年度額外公用工程能耗費用相對增加12.1%。

    圖11 年度公用能耗絕對增加量和年度公用能耗相對變化率

    4.2 無備用換熱器情況下的清洗計劃比較

    換熱器的性能衰退一方面會影響換熱器的運行安全,另一方面也會導致?lián)Q熱網絡的換熱效果下降嚴重,需要消耗更多的公用工程能耗滿足工藝生產要求。因此需要通過制定合理的清洗方案使得換熱器始終處于較高的換熱性能狀態(tài)運行,以減少公用能耗的增加。根據(jù)文獻[27]得到,單臺換熱器每次清洗的成本為2.5USD,因此可得HE1、HE2和HE5的清洗成本分別為359.7USD/次、363.2USD/次、439.3USD/次。在實際的工業(yè)換熱網絡中,如果換熱器存在備用換熱器,則清洗時只需要切換到備用換熱器,對生產的影響較??;如果換熱器不存在備用換熱器,則清洗時需要將換熱器從網絡中斷開。當換熱網絡中的關鍵換熱器沒有備用換熱器時,對性能衰退嚴重的換熱器進行清洗時,需要將換熱器從換熱網絡中斷開,導致在該設備停用期間換熱網絡需要更多的公用工程來滿足生產工藝要求,將這部分換熱器停用期間的公用工程能耗費用作為清洗的損失費用。

    (1)清洗方案1 當某臺換熱器性能衰退30%時,對換熱器進行清洗。

    如圖12 所示,當HE1、HE2 和HE5 三臺換熱器性能都衰退30%時制定清洗方案1,由于三臺換熱器都沒有備用換熱器,此時清洗時需要將三臺換熱器停用一天進行清洗維修。此時HE1、HE2兩臺換熱器的運行周期為85 天,HE1 和HE2 換熱器每年各需要清洗4 次,HE5 換熱器的運行周期47 天,HE5換熱器每年需要清洗7次,即每年總共需要清洗15 臺次。對圖12 中日額外需要的公用工程能耗進行累加,清洗方案1的年度額外公用工程能耗為3.43×10kJ/a,其中年度額外熱公用工程能耗為1.72×10kJ/a,年度額外冷公用工程能耗為1.72×10kJ/a,年度熱公用工程費用增加340578USD/a,年度冷公用工程費用增加28381USD/a,年度公用工程費用增加368959USD/a。如圖13所示,原始換熱網絡中,HE1回收熱量20800kW,HE2回收熱量25000kW,HE5回收熱量33000kW,停用期間需要將需要停用清洗的換熱器從網絡中移出,因此需要通過流股末端的公用換熱器提供額外的公用工程能耗來滿足生產工藝要求,因此計算可得,HE1 和HE2 停用時的公用工程能耗為7.91×10kJ/d,HE5停用時的公用工程能耗為5.7×10kJ/d。HE1和HE2換熱器每年需要停用4天,HE5換熱器每年需要停用7天,則因停用所需要的額外公用工程費用作為清洗方案的損失費用,清洗方案1 的損失費用為76923USD/a。換熱器清洗成本為5967USD/a,該清洗方案下年度總費用增加451849USD/a。與不清洗的原始網絡額外增加的公用工程能耗和額外增加的公用工程費用相比,年度額外公用工程能耗下降51.7%,年度額外公用工程費用減少51.6%,年度額外總費用減少40.8%。

    圖12 當HE1、HE2和HE5性能衰退30%、停用24h的清洗方案

    圖13 HE1、HE2和HE5停用24h網絡需要的額外公用工程能耗

    (2) 清洗方案2 當換熱器網絡性能衰退10%,綜合開展清洗決策。

    如圖14 所示,當HE1、HE2 和HE5 性能衰退導致?lián)Q熱網絡公用能耗增加10%時,對HE1、HE2和HE5 三臺換熱器同時進行清洗,制定清洗方案2,換熱器停用清洗時間為24h。按照該清洗方案2,此時HE1、HE2 和HE5 三臺換熱器的運行周期為50 天,HE1、HE2 和HE5 換熱器每年各需要清洗6 次,即每年總共需要清洗18 臺次。對圖14 的日額外需要的公用工程能耗進行累加,清洗方案2的年度額外公用工程能耗為2.98×10kJ/a,其中年度額外熱公用工程能耗為1.49×10kJ/a,年度額外冷公用工程能耗為1.49×10kJ/a,年度熱公用工程費用增加295968USD/a,年度冷公用工程費用增加24664USD/a,年度公用工程費用增加320632USD/a。如圖15 所示,HE1、HE2 停用時的公用工程能耗為7.91×10kJ/d,HE5 停用時的公用工程能耗為5.7×10kJ/d。HE1、HE2、HE5 每年都需要停用6天,則因停用所需要的額外公用工程費用作為清洗方案的損失費用,清洗方案的損失費用為87806USD/a。清洗成本為6793USD/a,該清洗方案下年度總費用增加415411USD/a。與不清洗的原始網絡額外增加的公用工程能耗和額外增加的公用工程費用相比,年度額外公用工程能耗下降58%,年度額外公用工程費用減少58%,年度額外總費用減少45.6%。

    圖14 當換熱器網絡性能衰退10%、停用24h的清洗方案

    圖15 HE1、HE2和HE5停用24h網絡需要的額外公用工程能耗

    以上兩種清洗方案分別代表了兩種不同的應用場景下,提高現(xiàn)役換熱器運行性能的清洗方案,換熱器停用清洗時間為24h。從圖16中可以看出,與清洗方案1相比,清洗方案2的年度額外公用工程費用減少13.1%,損失費用減少14.1%,清洗次數(shù)增加3臺次,清洗費用增加8.1%,但年度總費用減少13.8%。這主要是因為基于單臺換熱器性能衰退幅度制定的清洗方案1,換熱器運行周期較長,換熱器長時間處在較差的運行性能,導致?lián)Q熱網絡的額外公用工程能耗較多,額外公用工程費用較高,清洗次數(shù)較少,清洗費用和由于停用清洗所帶來的損失費用較少。由于換熱網絡性能的整體變化要比單臺換熱器性能變化更顯著,監(jiān)測網絡性能變化,不僅可以實現(xiàn)對單臺換熱器性能的監(jiān)測,同時可以反映不同換熱器性能衰退對網絡性能的影響。因此,根據(jù)換熱網絡性能衰退幅度制定清洗方案較傳統(tǒng)的清洗方案具有更好的經濟優(yōu)勢。

    圖16 兩種方案之間的比較

    5 結論

    本文針對原油精餾單元的換熱網絡,從換熱網絡集成的角度,根據(jù)網絡中各臺換熱器性能變化獲得換熱網絡的整體性能變化,提出基于智能預測和機理模型的換熱網絡性能預測方法;在此基礎上開展清洗決策,并開展了相應的案例研究,討論分析網絡性能變化和不同清洗準則下清洗方案的優(yōu)劣,主要得到了以下結論。

    (1)對原油換熱器性能的預測,采用污垢閾值模型的經驗公式,計算過程較為復雜且結果精度較差。通過對換熱器運行數(shù)據(jù)建立神經網絡算法的智能預測模型,可以較好地實現(xiàn)對換熱器性能的預測。

    (2)換熱網絡中換熱器之間存在耦合關系,不同換熱器的性能衰退對換熱網絡的影響不同,多臺換熱器同時發(fā)生性能衰退時,換熱網絡的性能衰退要比單臺換熱器性能衰退更顯著。當HE1、HE2和HE5三臺換熱器都發(fā)生性能衰退時,年度額外需要公用工程能耗相對增加高達14.1%,年度額外公用工程能耗費用相對增加12.1%。

    (3)與傳統(tǒng)基于單臺換熱器性能衰退制定的清洗方案相比,根據(jù)換熱網絡整體性能衰退制定的清洗方案,年度額外公用工程費用減少13.1%,損失費用減少14.1%,年度總費用減少13.8%,而清洗次數(shù)僅增加3臺次。

    —— 換熱面積,m

    —— 熱容流率,kW/℃

    CU—— 冷公用工程

    —— 活化能,kJ/mol

    HE—— 換熱器

    HU—— 熱公用工程

    LMTD—— 對數(shù)平均溫差,℃

    MAE—— 平均絕對誤差

    MSE—— 平均方差

    —— 壓強,Pa

    —— 普朗特數(shù)

    —— 熱負荷,kW

    —— 氣體常數(shù),J/(mol·K)

    R—— R

    —— 污垢熱阻,m·K/kW

    —— 雷諾數(shù)

    RMSE—— 均方根誤差

    —— 冷物流進口溫度,℃

    —— 冷物流出口溫度,℃

    —— 薄膜溫度,K

    —— 熱物流進口溫度,℃

    —— 熱物流出口溫度,℃

    —— 管壁溫度,K

    —— 總傳熱系數(shù),kW/(m·℃)

    —— 時間,s

    —— 平均流速,m/s

    —— 結垢模型中的常數(shù)

    —— 結垢模型中的常數(shù)

    —— 結垢模型中的常數(shù)

    —— 管壁剪應力,Pa

    下角標

    C—— 冷物流

    H—— 熱物流

    w—— 管壁物理量

    f—— 結垢條件下評估的物理量

    猜你喜歡
    公用結垢換熱器
    ASM-600油站換熱器的國產化改進
    能源工程(2021年3期)2021-08-05 07:26:14
    油田采出水結垢影響因素分析
    翅片管式換熱器的傳熱研究進展
    制冷(2019年2期)2019-12-09 08:10:36
    一個公用品牌的養(yǎng)成——橫山羊肉是咋樣“吃香”的
    當代陜西(2019年6期)2019-04-17 05:04:08
    Study on Scaling Law of Oilfield Produced Water with High Ca2+ and Ba2+
    公用電梯自動取消??垦b置初步設計
    電子制作(2017年17期)2017-12-18 06:40:45
    醫(yī)生私車公用撞傷人 醫(yī)院擔責
    310S全焊式板殼型換熱器的焊接
    焊接(2016年10期)2016-02-27 13:05:35
    “私車公用”打錯“方向盤”
    清風(2014年10期)2014-09-08 13:11:04
    大直徑重疊式浮頭式換熱器的制造
    低溫與特氣(2014年4期)2014-03-20 13:36:49
    亚洲图色成人| 久久久亚洲精品成人影院| 国产高清国产精品国产三级 | 精品一品国产午夜福利视频| 夜夜骑夜夜射夜夜干| 久久综合国产亚洲精品| 亚洲色图av天堂| 国产免费视频播放在线视频| 日日啪夜夜撸| 国产免费又黄又爽又色| 大片免费播放器 马上看| 成人毛片60女人毛片免费| 成年av动漫网址| 日日摸夜夜添夜夜爱| 日本av免费视频播放| 精品国产一区二区三区久久久樱花 | av天堂中文字幕网| 毛片女人毛片| 亚洲伊人久久精品综合| 亚洲,欧美,日韩| 午夜激情福利司机影院| 男女国产视频网站| 日本午夜av视频| 伊人久久精品亚洲午夜| 在现免费观看毛片| 精品国产三级普通话版| 国产精品99久久久久久久久| av线在线观看网站| 能在线免费看毛片的网站| 特大巨黑吊av在线直播| 日韩 亚洲 欧美在线| 超碰av人人做人人爽久久| 国产成人免费观看mmmm| 波野结衣二区三区在线| videos熟女内射| 久久国内精品自在自线图片| 日韩av免费高清视频| 日韩人妻高清精品专区| 尾随美女入室| 国产 精品1| 搡老乐熟女国产| www.色视频.com| 久久99蜜桃精品久久| 交换朋友夫妻互换小说| 欧美日本视频| 少妇熟女欧美另类| 黑人猛操日本美女一级片| 亚洲av中文av极速乱| 国产在线视频一区二区| 亚洲欧美中文字幕日韩二区| 国国产精品蜜臀av免费| 日本vs欧美在线观看视频 | 亚洲国产欧美在线一区| 久久久久久久久久人人人人人人| 老女人水多毛片| 在现免费观看毛片| 久久久精品94久久精品| 美女主播在线视频| 国产精品一区二区性色av| 网址你懂的国产日韩在线| 免费人成在线观看视频色| 一区二区三区免费毛片| 少妇人妻 视频| 夜夜看夜夜爽夜夜摸| 老师上课跳d突然被开到最大视频| 麻豆国产97在线/欧美| 日本午夜av视频| 夫妻性生交免费视频一级片| 久久av网站| 超碰97精品在线观看| 婷婷色综合www| 欧美xxxx性猛交bbbb| 成人国产av品久久久| 日韩不卡一区二区三区视频在线| 亚洲久久久国产精品| 亚洲人与动物交配视频| 卡戴珊不雅视频在线播放| 一区在线观看完整版| 99国产精品免费福利视频| 亚洲精品国产av蜜桃| 亚洲欧美中文字幕日韩二区| 欧美xxxx黑人xx丫x性爽| 国产精品女同一区二区软件| 亚洲怡红院男人天堂| 欧美区成人在线视频| 伊人久久国产一区二区| 干丝袜人妻中文字幕| 久久国产亚洲av麻豆专区| 亚洲国产精品专区欧美| 日韩av免费高清视频| 亚洲av综合色区一区| a 毛片基地| 成人无遮挡网站| 亚洲国产欧美在线一区| 欧美精品一区二区大全| 久久久久久久久久成人| 直男gayav资源| 亚洲av免费高清在线观看| 免费观看无遮挡的男女| 午夜免费鲁丝| 深爱激情五月婷婷| 激情五月婷婷亚洲| 日产精品乱码卡一卡2卡三| 国产熟女欧美一区二区| 成人免费观看视频高清| 一个人看的www免费观看视频| 18+在线观看网站| 在线观看人妻少妇| 91精品国产九色| 国产在视频线精品| 精品一区在线观看国产| 日韩大片免费观看网站| 七月丁香在线播放| 国产成人精品婷婷| av免费观看日本| 人体艺术视频欧美日本| 自拍偷自拍亚洲精品老妇| 成人亚洲欧美一区二区av| 成人午夜精彩视频在线观看| 狠狠精品人妻久久久久久综合| 国产亚洲精品久久久com| 美女福利国产在线 | 视频中文字幕在线观看| 色哟哟·www| 日韩中字成人| 欧美xxⅹ黑人| 伊人久久精品亚洲午夜| 成人午夜精彩视频在线观看| 99热这里只有是精品50| 久久久久久久大尺度免费视频| 在线免费观看不下载黄p国产| 精品久久久久久久末码| 国产高清国产精品国产三级 | 日韩成人av中文字幕在线观看| 国产亚洲午夜精品一区二区久久| 高清av免费在线| 少妇的逼好多水| 亚洲av男天堂| 国产精品福利在线免费观看| 免费观看av网站的网址| 精品酒店卫生间| 国产精品国产av在线观看| 又爽又黄a免费视频| 国产美女午夜福利| 夜夜看夜夜爽夜夜摸| h视频一区二区三区| 新久久久久国产一级毛片| 国产 精品1| 少妇猛男粗大的猛烈进出视频| 美女主播在线视频| 全区人妻精品视频| 黄色一级大片看看| 日韩欧美一区视频在线观看 | 亚洲欧洲国产日韩| 日韩一本色道免费dvd| 久久热精品热| 免费高清在线观看视频在线观看| 国产极品天堂在线| 人体艺术视频欧美日本| 在线观看一区二区三区激情| 在线免费十八禁| 最近手机中文字幕大全| 91精品伊人久久大香线蕉| 在线观看美女被高潮喷水网站| 一级毛片 在线播放| 精华霜和精华液先用哪个| 五月伊人婷婷丁香| 日日撸夜夜添| 激情五月婷婷亚洲| 七月丁香在线播放| www.色视频.com| 99热网站在线观看| av国产免费在线观看| 亚洲熟女精品中文字幕| 日韩电影二区| 日本爱情动作片www.在线观看| 天美传媒精品一区二区| av国产久精品久网站免费入址| 国产精品久久久久成人av| 夜夜骑夜夜射夜夜干| 91精品国产国语对白视频| 精品少妇久久久久久888优播| 97在线视频观看| 我的女老师完整版在线观看| 免费不卡的大黄色大毛片视频在线观看| 久久精品国产自在天天线| 亚洲av.av天堂| 深爱激情五月婷婷| 国产精品久久久久久av不卡| 国产成人免费无遮挡视频| 国产在视频线精品| 丝袜脚勾引网站| 美女高潮的动态| 老师上课跳d突然被开到最大视频| 久久精品国产鲁丝片午夜精品| 如何舔出高潮| 亚洲av不卡在线观看| 欧美日韩亚洲高清精品| 国产91av在线免费观看| 人妻少妇偷人精品九色| h日本视频在线播放| 精品熟女少妇av免费看| 26uuu在线亚洲综合色| 狂野欧美激情性xxxx在线观看| 亚洲婷婷狠狠爱综合网| 性色av一级| 少妇的逼水好多| 大陆偷拍与自拍| 国产精品偷伦视频观看了| 欧美少妇被猛烈插入视频| 欧美xxxx性猛交bbbb| 极品少妇高潮喷水抽搐| 国产精品人妻久久久久久| 精品国产露脸久久av麻豆| 日本午夜av视频| 18禁在线播放成人免费| 国产精品久久久久成人av| 九九爱精品视频在线观看| 91狼人影院| 亚洲精品aⅴ在线观看| 少妇人妻 视频| 又粗又硬又长又爽又黄的视频| 日日摸夜夜添夜夜爱| 亚洲综合精品二区| 99热全是精品| 久久毛片免费看一区二区三区| 国产永久视频网站| 最近最新中文字幕免费大全7| 午夜福利影视在线免费观看| 精品亚洲成国产av| 国产高清不卡午夜福利| 久久99蜜桃精品久久| 精品久久久久久久久av| 人人妻人人添人人爽欧美一区卜 | 精品人妻一区二区三区麻豆| 人妻制服诱惑在线中文字幕| av在线观看视频网站免费| 久久精品熟女亚洲av麻豆精品| 亚洲欧美一区二区三区黑人 | 成人一区二区视频在线观看| 国产精品.久久久| 91aial.com中文字幕在线观看| 大片电影免费在线观看免费| 国产精品嫩草影院av在线观看| tube8黄色片| 黄色日韩在线| 99久久精品热视频| 麻豆成人午夜福利视频| 亚洲欧美成人精品一区二区| 国产精品福利在线免费观看| 国产高清有码在线观看视频| 久久精品国产鲁丝片午夜精品| 插逼视频在线观看| 成人国产av品久久久| 成人特级av手机在线观看| 纯流量卡能插随身wifi吗| 亚洲av欧美aⅴ国产| av在线播放精品| 亚洲丝袜综合中文字幕| 18禁裸乳无遮挡动漫免费视频| 超碰97精品在线观看| 高清黄色对白视频在线免费看 | 国产精品蜜桃在线观看| kizo精华| 免费久久久久久久精品成人欧美视频 | 亚洲丝袜综合中文字幕| 色吧在线观看| 国产精品久久久久久精品古装| 国产黄片美女视频| 久久99热这里只频精品6学生| 免费少妇av软件| 精品国产乱码久久久久久小说| 熟女人妻精品中文字幕| 我的老师免费观看完整版| 国产av国产精品国产| 成人综合一区亚洲| 成人国产av品久久久| 亚洲av免费高清在线观看| 欧美激情极品国产一区二区三区 | 国产成人精品福利久久| 少妇裸体淫交视频免费看高清| h视频一区二区三区| 久久久久久伊人网av| 日本黄色片子视频| 婷婷色麻豆天堂久久| 午夜福利影视在线免费观看| 黄色日韩在线| 高清不卡的av网站| 亚洲国产精品999| 高清av免费在线| 国产一级毛片在线| 国产男女超爽视频在线观看| 日本wwww免费看| 看十八女毛片水多多多| 欧美成人午夜免费资源| 国产精品.久久久| 国产成人91sexporn| 热99国产精品久久久久久7| 久久这里有精品视频免费| 久久99蜜桃精品久久| 久久 成人 亚洲| 国产 一区 欧美 日韩| 精品人妻偷拍中文字幕| 精品一区在线观看国产| 男人狂女人下面高潮的视频| 女人久久www免费人成看片| 亚洲精品乱久久久久久| 精品人妻熟女av久视频| 国产成人一区二区在线| 亚洲怡红院男人天堂| 久久久久网色| 欧美激情国产日韩精品一区| 日本wwww免费看| 黄色视频在线播放观看不卡| 伦精品一区二区三区| 欧美bdsm另类| 婷婷色av中文字幕| 欧美精品一区二区免费开放| 一级黄片播放器| 国产免费福利视频在线观看| 2021少妇久久久久久久久久久| 国产高清不卡午夜福利| 精品久久久久久久末码| 国产精品一区二区在线观看99| 18禁在线播放成人免费| 中文资源天堂在线| 亚洲精品国产色婷婷电影| 青春草视频在线免费观看| 一级毛片 在线播放| 成年av动漫网址| 欧美日韩视频精品一区| 亚洲av综合色区一区| 这个男人来自地球电影免费观看 | 黄色怎么调成土黄色| 赤兔流量卡办理| 爱豆传媒免费全集在线观看| 高清在线视频一区二区三区| 色吧在线观看| 久久精品国产亚洲av天美| 精华霜和精华液先用哪个| 日本黄色片子视频| 黄片无遮挡物在线观看| 久热久热在线精品观看| av国产免费在线观看| 中文乱码字字幕精品一区二区三区| 欧美日本视频| 久久99蜜桃精品久久| 一本久久精品| 国产成人一区二区在线| 亚洲综合色惰| 日本欧美国产在线视频| 久久这里有精品视频免费| 欧美精品一区二区大全| 日本av手机在线免费观看| 在线观看人妻少妇| 午夜福利网站1000一区二区三区| 97超视频在线观看视频| 五月天丁香电影| 精品一区二区三卡| freevideosex欧美| av国产精品久久久久影院| 亚洲精品456在线播放app| 国产av一区二区精品久久 | av免费在线看不卡| 日韩 亚洲 欧美在线| 日韩欧美一区视频在线观看 | 国产精品成人在线| 男人和女人高潮做爰伦理| 欧美+日韩+精品| 国产成人91sexporn| av不卡在线播放| av专区在线播放| 亚洲成人中文字幕在线播放| 九九爱精品视频在线观看| 午夜免费观看性视频| 精品国产乱码久久久久久小说| av不卡在线播放| 欧美老熟妇乱子伦牲交| 91aial.com中文字幕在线观看| 久久青草综合色| 日韩电影二区| 精品国产三级普通话版| 在线 av 中文字幕| 亚洲色图av天堂| 久久久精品94久久精品| 女的被弄到高潮叫床怎么办| 高清午夜精品一区二区三区| 噜噜噜噜噜久久久久久91| 成年av动漫网址| 久久久久久九九精品二区国产| 一区二区三区四区激情视频| 欧美一级a爱片免费观看看| 五月伊人婷婷丁香| 高清欧美精品videossex| 国产毛片在线视频| 亚洲av.av天堂| 久久人人爽人人片av| 人体艺术视频欧美日本| 激情 狠狠 欧美| 精品亚洲乱码少妇综合久久| 国产av一区二区精品久久 | 久久99热这里只有精品18| 久久久色成人| 久久毛片免费看一区二区三区| 亚洲一区二区三区欧美精品| 青春草亚洲视频在线观看| 男人狂女人下面高潮的视频| 黄色配什么色好看| 精品国产乱码久久久久久小说| 久久韩国三级中文字幕| 久久久久久久久久久丰满| 久久国产精品大桥未久av | 99热国产这里只有精品6| 最近手机中文字幕大全| 国产精品不卡视频一区二区| 日韩中字成人| 亚洲av免费高清在线观看| 中文资源天堂在线| 七月丁香在线播放| 乱系列少妇在线播放| 18禁裸乳无遮挡免费网站照片| 国产日韩欧美亚洲二区| 久久久亚洲精品成人影院| 啦啦啦中文免费视频观看日本| 99热这里只有精品一区| 男女边摸边吃奶| 欧美xxxx黑人xx丫x性爽| 久久这里有精品视频免费| 亚洲欧洲国产日韩| 欧美丝袜亚洲另类| 欧美日韩一区二区视频在线观看视频在线| 亚洲国产精品一区三区| 精品一品国产午夜福利视频| 一区在线观看完整版| 在线免费观看不下载黄p国产| 国产男女超爽视频在线观看| 中国国产av一级| av黄色大香蕉| 欧美日韩在线观看h| 在现免费观看毛片| 性高湖久久久久久久久免费观看| 国产亚洲最大av| 久久久a久久爽久久v久久| 精华霜和精华液先用哪个| 久久精品国产亚洲网站| 精品人妻偷拍中文字幕| 亚洲欧美清纯卡通| 人妻制服诱惑在线中文字幕| 久久人人爽av亚洲精品天堂 | 亚洲真实伦在线观看| 我要看黄色一级片免费的| 人妻制服诱惑在线中文字幕| 亚洲国产欧美人成| 内地一区二区视频在线| 亚洲精品,欧美精品| 亚洲精华国产精华液的使用体验| 直男gayav资源| 亚洲欧美日韩无卡精品| 下体分泌物呈黄色| 午夜福利在线观看免费完整高清在| 国产黄片视频在线免费观看| 少妇的逼水好多| 免费久久久久久久精品成人欧美视频 | 多毛熟女@视频| 精品亚洲成国产av| 26uuu在线亚洲综合色| 久热久热在线精品观看| 亚洲真实伦在线观看| 最近2019中文字幕mv第一页| 亚洲第一区二区三区不卡| 日韩免费高清中文字幕av| 亚洲av二区三区四区| 午夜免费鲁丝| 国产色婷婷99| 美女国产视频在线观看| 韩国高清视频一区二区三区| 高清黄色对白视频在线免费看 | 最黄视频免费看| 亚洲人与动物交配视频| 国产乱人视频| 热99国产精品久久久久久7| 精品国产露脸久久av麻豆| 欧美zozozo另类| 97超碰精品成人国产| xxx大片免费视频| 国产高清三级在线| 99久久精品热视频| 欧美 日韩 精品 国产| 国产极品天堂在线| 能在线免费看毛片的网站| 另类亚洲欧美激情| 26uuu在线亚洲综合色| 久久精品国产亚洲网站| 午夜福利在线观看免费完整高清在| 免费av中文字幕在线| 欧美xxxx性猛交bbbb| 高清视频免费观看一区二区| 亚洲美女视频黄频| 91久久精品国产一区二区成人| 熟妇人妻不卡中文字幕| 日日啪夜夜撸| 欧美日韩视频精品一区| 久久国产亚洲av麻豆专区| 国产精品人妻久久久久久| 精品久久久久久久久亚洲| 久久久色成人| 免费看不卡的av| 99久久精品热视频| 国产片特级美女逼逼视频| 五月开心婷婷网| 日韩不卡一区二区三区视频在线| 亚洲精品第二区| 精品少妇黑人巨大在线播放| 各种免费的搞黄视频| 亚洲国产最新在线播放| 99re6热这里在线精品视频| 赤兔流量卡办理| 亚洲国产av新网站| 熟女人妻精品中文字幕| 久久av网站| 又粗又硬又长又爽又黄的视频| 少妇熟女欧美另类| videossex国产| 久久毛片免费看一区二区三区| 成人黄色视频免费在线看| 亚洲欧美中文字幕日韩二区| 久久精品国产亚洲网站| 永久网站在线| 免费久久久久久久精品成人欧美视频 | 一区二区三区乱码不卡18| 日本欧美视频一区| 深夜a级毛片| 欧美成人一区二区免费高清观看| 久久久久久久久久成人| 九九久久精品国产亚洲av麻豆| 韩国av在线不卡| 麻豆成人av视频| 多毛熟女@视频| 熟妇人妻不卡中文字幕| 国产黄频视频在线观看| av卡一久久| 久久亚洲国产成人精品v| 欧美日韩视频精品一区| 秋霞在线观看毛片| 在线看a的网站| 有码 亚洲区| 精品亚洲成国产av| 久久99精品国语久久久| 久久久久精品性色| 国产一区二区三区综合在线观看 | 久久久久精品性色| 九草在线视频观看| 美女国产视频在线观看| 国产精品久久久久久久久免| 你懂的网址亚洲精品在线观看| 一区二区三区四区激情视频| 国产成人精品一,二区| 身体一侧抽搐| 久久精品久久精品一区二区三区| 国产精品人妻久久久久久| 久久精品久久久久久噜噜老黄| 性色av一级| 日本免费在线观看一区| 大片免费播放器 马上看| 午夜免费观看性视频| 国产精品人妻久久久影院| 欧美丝袜亚洲另类| 成人毛片a级毛片在线播放| 日本爱情动作片www.在线观看| 97超碰精品成人国产| 久久久久久久久大av| a 毛片基地| 亚洲av欧美aⅴ国产| 精品一区二区三区视频在线| 亚洲欧美日韩东京热| 久久综合国产亚洲精品| 观看免费一级毛片| 菩萨蛮人人尽说江南好唐韦庄| 狂野欧美白嫩少妇大欣赏| 高清黄色对白视频在线免费看 | 婷婷色麻豆天堂久久| 国产 精品1| 亚洲最大成人中文| 色视频在线一区二区三区| 老师上课跳d突然被开到最大视频| www.色视频.com| 九九爱精品视频在线观看| 午夜日本视频在线| 在线免费十八禁| 成人国产麻豆网| 日本黄色片子视频| 日本黄大片高清| 在线 av 中文字幕| 午夜福利在线观看免费完整高清在| 性色avwww在线观看| 久久亚洲国产成人精品v| 搡老乐熟女国产| 免费黄频网站在线观看国产| 亚洲av.av天堂| 亚洲精品中文字幕在线视频 | 狂野欧美白嫩少妇大欣赏| 免费看光身美女| 深爱激情五月婷婷| 美女中出高潮动态图| 国产国拍精品亚洲av在线观看| 亚洲精品456在线播放app| 婷婷色综合大香蕉| 特大巨黑吊av在线直播| 国产精品国产三级专区第一集| 熟女人妻精品中文字幕| 大陆偷拍与自拍| 一本久久精品| 久久鲁丝午夜福利片| 亚洲色图综合在线观看| 80岁老熟妇乱子伦牲交|