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

    基于功率-溫度自適應控制的多堆質子交換膜電解制氫系統(tǒng)效率優(yōu)化

    2024-04-09 05:21:12韓鵬飛徐瀟源
    電工技術學報 2024年7期
    關鍵詞:適應控制電解槽制氫

    韓鵬飛 徐瀟源 王 晗 嚴 正

    基于功率-溫度自適應控制的多堆質子交換膜電解制氫系統(tǒng)效率優(yōu)化

    韓鵬飛1,2徐瀟源1,2王 晗1,2嚴 正1,2

    (1. 電力傳輸與功率變換控制教育部重點實驗室(上海交通大學) 上海 200240 2. 上海非碳基能源轉換與利用研究院 上海 200240)

    隨著可再生能源發(fā)電及電解制氫技術的發(fā)展,基于質子交換膜(PEM)技術的可再生能源制氫將成為氫能供應的主要手段之一。新能源發(fā)電具有隨機性和波動性,現有多堆PEM電解制氫系統(tǒng)功率分配策略難以在功率波動場景下保證系統(tǒng)高效運行。對此,該文首先建立了包含主要能耗環(huán)節(jié)的制氫系統(tǒng)整體效率模型,探究了制氫效率與電流密度、溫度的特性關系。然后,根據該效率特性提出了基于功率-溫度自適應控制的多堆PEM電解制氫系統(tǒng)效率優(yōu)化策略,包含離線優(yōu)化與在線控制兩部分:離線優(yōu)化部分根據制氫系統(tǒng)效率模型預先制定不同電解總功率下的電解槽功率及溫度設定方案;在線控制部分根據離線優(yōu)化方案、制氫功率進行電解槽功率和溫度調節(jié)。最后,采用實際風電功率進行算例分析,結果表明,相比于傳統(tǒng)的功率平均分配策略和鏈式分配策略,該文提出的效率優(yōu)化策略將產氫量分別提高了6.4%和5.7%。

    質子交換膜電解槽 電解制氫 新能源發(fā)電 效率優(yōu)化 自適應控制

    0 引言

    隨著現代社會對清潔能源的重視,氫能在能源、交通等領域不斷普及[1-2]。預計到2050年,氫能在全球終端能源消費量中的占比可達到25%[3]??稍偕茉窗l(fā)電制氫已成為實現“雙碳”目標、推動能源供給可持續(xù)發(fā)展的重要途徑。

    電解水制氫技術主要有堿性電解、質子交換膜(Proton Exchange Membrane, PEM)電解及高溫固體氧化物電解技術。近年來,PEM電解技術逐步成熟,呈現出電流密度高、響應速度快的優(yōu)勢,適應功率波動的制氫場景[4-5]。多堆PEM電解制氫成為有效消納富余可再生能源的重要手段,其運行方式從配合新能源消納的被動負載逐步轉換為與電力系統(tǒng)雙向互動的柔性負載。在電力系統(tǒng)運行經濟性方面,文獻[6-8]在區(qū)域電網經濟調度研究中通過優(yōu)化電解制氫功率計劃,提升了可再生能源消納水平。在電力系統(tǒng)動態(tài)性能控制方面,通過控制電解制氫系統(tǒng)中的電力電子換流器,可以在電力系統(tǒng)動態(tài)過程中提供電壓[9]/頻率[10]支撐,以及平抑功率波動[11-12]。上述研究使電解制氫系統(tǒng)運行或控制更加靈活,但缺少對電解制氫系統(tǒng)自身能耗特性的關注。

    對于大型PEM電解制氫系統(tǒng),電能成本約占全生命周期成本的70%[13],提升能耗效率對于降低制氫成本至關重要。文獻[14]詳細地分析了電解制氫系統(tǒng)中的電解、干燥、壓縮等環(huán)節(jié)的能耗,得出不同條件下的經濟運行點,為提升電解制氫效率提供了理論指導。文獻[15]建立電解槽動態(tài)模型,探究了不同壓強下電解槽電壓特性及能耗特性。文獻[16]探究了電解槽壓強、PEM膜厚度對制氫全過程效率的影響,得到不同電解電流密度下的最優(yōu)壓強設定值。文獻[17]探究了溫度和壓強變化對電解效率的影響,并通過溫度與壓強優(yōu)化得到最佳運行效率。文獻[18]在風氫系統(tǒng)優(yōu)化配置中考慮電解效率,提出了風氫系統(tǒng)容量配置方法。文獻[19]分析了溫度對電解率的影響,并對電解槽溫度進行了自適應溫度控制,使電解槽運行于最佳溫度點。上述研究著眼于電解槽本體效率特性,但未將制氫系統(tǒng)中輔助設備的能耗特性考慮在內,也未考慮多電解槽之間的功率分配問題。

    在多堆電解槽功率分配方面,文獻[20-21]針對堿性電解槽動態(tài)響應慢的問題,提出適用于可再生能源制氫的堿性電解槽定期輪值策略,從而兼顧電解槽壽命和制氫效率。文獻[22]充分利用PEM電解槽良好的變載能力,提出根據可再生能源出力逐級啟動電解槽的鏈式分配策略,降低系統(tǒng)的控制復雜度。文獻[23]提出多堆PEM電解槽統(tǒng)一控制結構,多組電解槽并列于同一條直流母線,無需獨立的功率控制器,并依靠電解槽相同的外特性實現制氫功率的平均分配。上述多堆電解槽控制策略具有操作簡便、可實施性強的優(yōu)勢,適用于制氫功率相對恒定的場景。然而其設置方式相對主觀,在可再生能源出力頻繁波動時,難以保證系統(tǒng)的高效運行。

    綜上所述,現有針對多堆PEM電解槽制氫效率提升的研究還不夠完善,需進一步探究電解制氫系統(tǒng)的整體效率特性,提出適用于可再生能源制氫場景的多堆電解槽功率分配方法。對此,本文對電解制氫系統(tǒng)主要能耗環(huán)節(jié)進行了詳細分析,建立了包含換流器、電解槽、多級壓縮機的整體效率模型,并基于該模型從優(yōu)化和控制兩個層面提出了基于功率-溫度自適應控制的多堆PEM電解系統(tǒng)效率優(yōu)化方法。該方法首先通過離線優(yōu)化預先制定不同電解功率下的電解槽功率及溫度設定方案;其次,在線控制換流器占空比及循環(huán)水閥門開度,實現電解槽功率和溫度調節(jié),從而使電解槽在不同的制氫負載下最大化制氫效率;最后,通過實時數字仿真平臺RT-LAB與熱力瞬態(tài)仿真工具TRNSYS進行了聯(lián)合仿真驗證,并與傳統(tǒng)功率分配策略進行對比,驗證了本文所提效率優(yōu)化方法的有效性。

    1 電解制氫系統(tǒng)及其效率模型

    PEM電解制氫的主要生產環(huán)節(jié)[14,24]如圖1所示。在陰極生成的氫氣經過氣液分離、干燥、壓縮后進行存儲或管道運輸。對于大規(guī)模PEM電解制氫系統(tǒng),電解過程中產生大量熱量,為防止PEM過熱老化損壞,需要將水進行循環(huán)流動并采用換熱器散熱。在整個制氫流程中,主要耗能環(huán)節(jié)有電解槽制氫、整流器耗能及壓縮耗能,干燥、循環(huán)冷卻等輔助環(huán)節(jié)耗能占比通常低于3%[25-26]。因此,本文分別對電解、壓縮、換流器三個主要能耗環(huán)節(jié)進行效率分析,以獲得表征效率特性的數學模型。

    圖1 PEM電解水制氫系統(tǒng)主要生產環(huán)節(jié)[14,24]

    1.1 電解槽及其效率模型

    大型電解槽通常由多個電解小室串聯(lián)而成。單個電解小室的電解電壓可以分解為能斯特電壓(Nernst voltage)、活化過電勢(activation overpotential)、歐姆過電勢(Ohmic overpotential)三個部分[17,27],即

    其中

    式中,cell和cell分別為單個電解小室的電壓和溫度;N為能斯特電壓;act為活化過電勢;res為歐姆過電勢;an為陽極電勢;cat為陰極電勢;0 rev為可逆電壓;為理想氣體常數;為法拉第常數;為轉移電子摩爾數;為電荷轉移系數;為電流密度;0為交換電流密度;an和cat分別為陽極和陰極壓強;0為參考壓強;m、m、m分別為質子交換膜厚度、膨脹系數、電解膜電導率;m H2O為質子交換膜中的水分活度;0為電解液電阻。

    電解制氫本身為吸熱反應。在標準狀況下,當電解電壓達到理論分解電壓1.23 V時,電解反應開始發(fā)生,并需要從外部吸收熱量48.6 kJ/mol。然而在該電壓下,電解反應非常緩慢。工業(yè)生產中為了提升制氫速率,通常會提升電解電流,將不可避免地增加歐姆過電勢并釋放出焦耳熱。當電壓達到熱中性電壓1.48 V時,焦耳熱與電解吸收熱量動態(tài)平衡,電解反應表現出既不發(fā)熱也不吸熱。此后繼續(xù)提升電流,電解速率將加快,需要采取散熱措施。

    電解效率存在不同的定義方式。若將電解效率定義為反應所生成氫氣的低熱值與消耗電量的比值,則由上述分析可知,電解效率模型可以采用兩階段效率函數表示:在電解電壓低于熱中性電壓1.48 V時,電解消耗的電能完全轉化為氫氣中的化學能和增加的水蒸氣中的焓值,并且需要從外部吸收熱量,這一階段電解效率可用熱力學效率計算,如式(5)所示;當電壓高于熱中性電壓1.48 V時,電能傳輸產生的熱量高于電解吸收的熱能,需要采用散熱措施以保證電解溫度恒定,該階段電解效率可用電壓效率計算,如式(7)所示。

    另一方面,由于氫氣分子直徑很小,反應生成的部分氫氣會透過質子交換膜到陽極發(fā)生逆反應,降低氫氣產量??紤]這部分損失后的效率稱為法拉第效率F H2,可表示為

    綜上所述,電解槽效率ELE可以表示為

    由于電解槽電壓受到電流密度和電解溫度的雙重影響,可基于電解效率模型,得到考慮產氫效率和法拉第效率的電解效率曲面,如圖2所示。由圖2可知,熱力學效率隨著電流密度的增大而增大,而電壓效率隨著電流密度的增大而降低。在電流密度為1~2 A/cm2的范圍內,熱力學效率曲面ηp H2,1和電壓效率曲面ηp H2,2相交。因此,電解效率在電流密度較小的范圍(約1 A/cm2以內)時較低,并隨著電流密度的增大而快速增大,最終電解效率隨著電流密度的增加而降低。此外,電解效率也受到電解溫度的影響。在電流密度很低時受溫度影響較小;而當電流密度較大時,隨著溫度的增加而提高。

    1.2 壓縮功耗及其效率模型

    氫氣壓縮是電解制氫系統(tǒng)中不可忽視的能耗環(huán)節(jié)。氫氣的儲運壓強通常維持在數MPa至數十MPa范圍內,其與電解槽壓強之間存在很大的壓強差。大規(guī)模電解制氫系統(tǒng)一般采用多級壓縮、級間冷卻的方式降低壓縮能耗[28],如圖3a所示。單級壓縮并冷卻后,氫氣恢復至初始溫度,隨后進入下一級壓縮。圖3b中陰影部分即為通過多級壓縮方式降低的功耗。假設電解槽內部壓強恒定,則單級壓縮功耗及壓縮級數由式(11)和式(12)給出。采用多級壓縮時,當每一級壓縮比相同時所消耗的壓縮能量最小[29-30]。因此,可以得出多級壓縮的效率模型,如式(13)所示。

    圖3 多級壓縮示意圖及熱力學過程

    式中,c為單級壓縮作用在氫氣上的有效功率;為氫氣的絕熱指數;in為氫氣壓縮前的溫度;為壓縮因子;c為單級壓縮比;c為壓縮級數;s為儲氫壓強;c為壓縮機工作效率;CH為氫氣多級壓縮過程的整體效率。

    1.3 換流器效率模型

    電解制氫負荷屬于典型低電壓大功率直流負載,正常運行中換流器必不可少[31-32]。以電解系統(tǒng)中廣泛使用的Buck降壓換流器為例分析DC-DC降壓變流器的運行效率。Buck電路功率損耗主要包含通態(tài)損耗、開通損耗和關斷損耗[33],即

    式中,con為換流器通態(tài)損耗;on為開通損耗;off為關斷損耗;為占空比;o為負載電流;T為MOSFET導通電阻;L為電感直流電阻;c為電容器等效串聯(lián)電阻;s為開關頻率;i和i分別為輸入端電壓和電流;on和off分別為MOSFET開通時間和關斷時間;con為降壓換流器效率。

    1.4 系統(tǒng)總效率

    基于上述分析,包含主要功率損耗環(huán)節(jié)的整體效率模型total為

    為了進一步探究與電解槽效率特性相適應的自適應控制策略,以額定功率為2 MW的電解制氫系統(tǒng)為例,繪制電解效率關于電流密度和電解溫度的曲面,如圖4所示。電解槽運行及輔助設備參數分別見表1和表2[17,34]。表中,1 bar=0.1 MPa。

    表1 電解槽運行參數

    Tab.1 Operating parameters of electrolyzers

    表2 輔助設備參數

    Tab.2 Auxiliary equipment parameters

    由圖4可知,電解效率并非隨著電流密度、電解溫度單調變化。當電解溫度固定時,電解效率隨著電流密度的增加先增大后降低,最大效率點出現在1.5~2 A/cm2范圍內,數值約為0.6~0.7。當電流密度固定時,在電流密度小于1.5 A/cm2時,電解效率隨著溫度的上升而下降;當電流密度超過1.5 A/cm2后,電解效率隨著溫度的升高而升高。在實際制氫系統(tǒng)中,通常會通過提升電解溫度來提升電解效率,但為了防止電解膜高溫老化,電解溫度應當保持在360 K以下[35]。

    由于電解溫度影響,電解效率曲面在電解效率total-電流密度平面上的投影是一條寬度隨電流密度變化而變化的區(qū)間,如圖5陰影部分所示。區(qū)間的寬度表示溫度對電解效率影響的程度。陰影區(qū)間寬度越大,表示該電流密度下溫度對電解效率影響越大。該區(qū)間上限表示在電解溫度可靈活調節(jié)情況下的電解槽最大效率曲線ηbest。為方便對比,同時給出了電解系統(tǒng)常用的固定溫度360 K下的制氫效率曲線η360。將ηbest通過圖4所示的電解系統(tǒng)效率曲面可對應得到最佳電解溫度best??梢钥闯?,最佳電解效率并非總是在最高電解溫度下實現。當電流密度較低時,最佳電解溫度相對較低,這是由于電流密度較低時,溫度升高將大大降低法拉第效率。當電流密度超過1.5 A/cm2時,最佳電解溫度隨著電流密度增大而增大。

    圖5 電解效率曲面投影及效率峰值對應的溫度曲線

    基于上述分析,對于大規(guī)模可再生能源電解制氫系統(tǒng),當可再生能源功率波動時,合理調整電解槽功率分配以及電解溫度,可以提升電解效率。電解功率可以通過調節(jié)DC-DC換流器實現,電解槽運行溫度可以通過循環(huán)水流速進行調節(jié)。

    2 制氫功率-溫度自適應控制策略

    現有制氫功率分配策略通?;诠潭▓鼍霸O定,在可再生能源出力頻繁波動的場景下,難以保證系統(tǒng)處于高效的狀態(tài)。對此,需探索一種能夠根據可再生能源實際出力自適應調節(jié)電解槽狀態(tài)的控制策略。

    圖6給出了電解功率-溫度可控的可再生能源制氫的系統(tǒng)結構,其中多組電解槽通過DC-DC換流器并聯(lián)在直流母線上,每組電解槽通過獨立的換流器進行功率控制,并通過冷卻功率可控的換熱器進行溫度調節(jié)。電解槽功率參考值由自適應控制器確定,溫度參考值可根據最大效率曲面獲得??紤]到電解系統(tǒng)效率模型計算相對復雜,在線實時優(yōu)化難以在有效時間內完成,因此將功率自適應分配策略分為離線優(yōu)化和在線控制兩個部分:離線優(yōu)化實現功率-溫度參考值生成;在線控制實現基于功率-溫度參考值的自適應調節(jié)控制。以下進行詳細論述。

    圖6 電解功率-溫度可控的多電解槽制氫系統(tǒng)結構

    2.1 離線優(yōu)化

    電解制氫系統(tǒng)結構通常不會隨時改變,因此可預先計算出制氫總功率從最小值到最大值變化時的電解槽功率和溫度設定方案,形成離線數據表。在實時運行中只需根據實際制氫總功率查表得到電解槽功率和溫度參考值。PEM電解槽功率可以在s級時間范圍內完成調節(jié),而電解槽溫度變化則相對較慢,需要min級時間完成調節(jié)[26]。因此,離線優(yōu)化環(huán)節(jié)分別以最小化制氫功率損耗、最小化溫度調節(jié)量為目標進行兩階段優(yōu)化。在第一階段實現最大化制氫效率,如式(19)所示。

    由于第一階段通常會存在多種電解槽功率分配或啟停方案,為盡快達到高效制氫狀態(tài),第二階段以最小化溫度調整為目標確定電解槽制氫功率,如式(20)所示。

    上述兩階段優(yōu)化模型普遍適用于多堆PEM電解制氫系統(tǒng)。實際的制氫系統(tǒng)中通常由多臺規(guī)格相同的電解槽構成,這種情況下制氫功率分配變得相對簡單??勺C明,對于一個由臺電解槽構成的電解系統(tǒng),如果已開啟(=2,…,)臺電解槽,將電解功率平均分配給臺電解槽時的效率大于其他分配方案,證明過程見附錄,則不同電解功率下最大化制氫效率問題轉變?yōu)殡娊獠坶_啟臺數問題。

    P-1,表示開啟-1臺電解槽與臺電解槽的啟停臨界功率(即當制氫功率為P-1,時,開啟-1臺電解槽與開啟臺電解槽效率相同),則只需在離線優(yōu)化階段計算出各個臨界啟停功率P-1,j,便可以根據實際制氫功率所屬范圍確定電解槽開啟臺數。臨界啟停功率可以通過式(21)獲得。

    圖7給出了由n臺相同規(guī)格電解槽構成的制氫系統(tǒng)在不同制氫功率下的分配方案。在該分配方案下,當制氫總功率PELE滿足Pj-1,j<PELE≤Pj, j+1時,每臺電解槽的功率為PELE/j。該方案可能導致不同電解槽的運行時間不均衡。為避免其使用壽命快速降低,可定期對相同規(guī)格的電解槽進行輪換[20]。

    2.2 在線控制

    在線控制階段,基于離線控制的參考指標對電解功率和溫度進行實時控制。電解功率采用雙閉環(huán)PI控制策略,如圖8下部所示,其中外環(huán)為功率環(huán),內環(huán)為電流環(huán)。電解槽溫度通過調節(jié)換熱器閥門開度進行控制,在滿足電解槽散熱約束的情況下,散熱功率由換流器低溫側液體流量進行調節(jié),如圖8上部所示。忽略電解槽與環(huán)境之間的被動散熱,則電解槽-換熱器溫度模型[36]為

    式中,e為由于焦耳效應引起的電加熱功率;ex為換熱器帶走的熱功率;ELE為電解槽溫度;ELE為電解槽內電解液的總質量;w為水的比熱容;q為閥門流量系數,取0.82;為閥門的流通面積;為閥門開度;為閥門兩側壓強差;pipe為電解槽補水溫度;pipe為電解液循環(huán)的質量流量。由式(22)可知,電解槽溫度屬于一階慣性環(huán)節(jié),因此可以采用PI控制器進行溫度控制。

    圖8 電解槽功率-溫度控制圖

    2.3 基于功率-溫度自適應控制的效率優(yōu)化策略

    經過上述分析,得到由離線優(yōu)化和在線控制兩部分構成的電解制氫系統(tǒng)功率-溫度控制策略,如圖9所示。離線優(yōu)化階段,根據系統(tǒng)結構、設備參數構建系統(tǒng)效率模型,基于該模型獲得電解功率從最小值到最大值的電解槽功率分配以及溫度設定方案;在線控制階段,基于系統(tǒng)當前制氫總功率,結合離線優(yōu)化階段預先制定的功率-溫度設定方案,獲取對應電解槽功率及溫度,然后通過Buck電路和閥門控制器生成PWM調制信號及閥門開度信號。

    圖9 基于功率-溫度自適應控制的效率優(yōu)化策略

    3 仿真對比驗證

    3.1 不同制氫功率情況下功率分配情況

    為驗證本文所提出的基于功率-溫度自適應控制的效率優(yōu)化策略的正確性,采用實時聯(lián)合仿真平臺進行驗證。功率-溫度自適應控制策略基于RT-LAB OP5600實時仿真平臺搭建,多電解槽制氫系統(tǒng)基于熱力瞬態(tài)實驗平臺TRNSYS構建。通過RT-LAB內置的socket接口與TRNSYS下的type155與Matlab COM模塊實現通信。實時仿真系統(tǒng)設置如圖10所示。為方便描述,將平均分配策略、鏈式分配策略、基于功率自適應控制的效率優(yōu)化策略、基于功率-溫度自適應控制的效率優(yōu)化策略分別用策略1~策略4表示;前三種控制策略采用固定電解溫度360 K,策略4的溫度根據圖5中最佳電解溫度曲線best設定。電解功率設置為從0開始逐漸增加至滿載,得到的分配結果及效率曲線如圖11和圖12所示。

    圖10 實時仿真系統(tǒng)設置

    圖11 不同策略下電解槽功率分配曲線

    圖12 不同策略的制氫效率曲線

    圖11表示當制氫總功率從0升至最大負載時,4臺電解槽(ELE1~ELE4)在不同策略下的功率分配情況。由圖11可知,采用平均分配策略(策略1)的制氫系統(tǒng),在制氫總功率較小時無法啟動,降低了低功率場景下新能源的利用率;鏈式分配策略(策略2)可根據制氫功率逐級啟動電解槽,對低功率利用效率高;采用功率自適應控制策略(策略3和4),在總制氫功率為0~4 MW范圍內,逐級開啟電解槽,使制氫效率曲線較大范圍地維持在最高效的區(qū)間。與溫度恒定的功率自適應控制策略(策略3)相比,功率-溫度自適應控制策略(策略4)使制氫效率最高點向左偏移,較為提前啟動下一級電解槽。

    上述分析可以從圖12的制氫總效率曲線上體現:采用平均分配策略時,電解槽在1 MW制氫功率以內不會啟動,因此低功率范圍利用率低,但是高功率范圍利用率相對較高;鏈式分配策略的效率最高峰覆蓋范圍較小,在系統(tǒng)沒有滿載的時候,效率低于其他三種策略;采用功率自適應分配策略時,控制器可以根據制氫總功率的變化自適應地分配電解槽運行功率,從而使制氫系統(tǒng)更大范圍地運行于最佳效率區(qū)間。與平均分配策略相比,功率自適應控制策略提升了對低功率的利用效率;與鏈式分配策略相比,功率自適應控制策略可以使整個制氫功率范圍的效率提升。功率-溫度自適應控制策略由于同時優(yōu)化功率和溫度,效率在0~5 MW范圍內相對于僅考慮功率優(yōu)化的策略3提升約2%;隨著制氫功率的增加,效率優(yōu)勢逐漸降低,制氫功率達到7 MW后,優(yōu)化后的溫度上升至360 K,策略3和策略4效率曲線重合。

    綜上所述,鏈式分配策略適用于制氫功率較低的場景;平均分配策略適用于制氫功率較大的場景;基于功率自適應控制的效率優(yōu)化策略兼具前兩種策略的優(yōu)勢,從低負載區(qū)間到高負載區(qū)間均保持較高的制氫效率;而基于功率-溫度自適應控制的效率優(yōu)化策略通過功率和溫度的雙重控制,可在0%~70%負載區(qū)間內進一步提升制氫效率。

    3.2 不同場景下的電解槽策略對比

    本文所提效率優(yōu)化策略適用于風電、光伏等可再生能源制氫場景,以風電制氫為例進行有效性驗證。制氫系統(tǒng)由額定容量為4×2 MW、過載率為120%的PEM電解槽構成,功率由裝機容量為10 MW的風電場提供?;陲L電場一年實測數據,采用Kantorovich距離[37]將歷史數據削減為低、中、高三種典型風電出力場景,對應概率分別為0.6、0.2、0.2。為便于描述,將三種風功率場景分別表示為場景1~場景3,典型制氫功率曲線如圖13所示。

    圖13 典型制氫功率曲線

    圖14和圖15分別給出了場景1下四種制氫策略的功率分配方案及制氫效率曲線。從功率分配方案可以看出:策略1由于采用平均分配策略,在制氫總功率較低時段,例如,600 min時,每臺電解槽分配的功率不滿足最小運行功率,系統(tǒng)停止制氫,降低了新能源利用率;策略2采用鏈式開啟策略,整體功率始終沒有達到開啟第4臺電解槽的階段;策略3和4采用了最大化制氫效率策略,在高功率時段啟用了四臺電解槽。在四種控制策略中,后三種策略都可以根據風功率變化調整電解槽的啟停。由圖15可以看出,采用功率-溫度自適應控制策略后,除了在1 100 min時總效率較低外,基本上可全過程保持最高的電解效率;功率自適應控制策略與功率平均分配策略效率曲線在大部分時段重合,在總制氫功率很低的時段,平均分配策略導致電解槽停止運行,而功率自適應分配策略還能使得電解槽保持較高的運行效率;鏈式分配策略根據制氫功率的大小調整電解槽啟停,因此可以利用較低的制氫功率,減少了停機的情況。但鏈式分配策略沒有進行效率優(yōu)化,在風功率出力較大的時段(0~600 min)內,運行效率低于其他三種分配策略。

    圖14 場景1下四種策略電解槽功率分配情況

    圖15 場景1下四種策略制氫效率

    圖16和圖17分別給出了場景2下四種制氫策略的功率分配方案及制氫效率曲線。在場景2下,由于功率未下降至最低范圍,采用策略1的四臺電解槽未退出運行,但每臺電解槽分配到的制氫功率很低;相比之下,策略2采用了鏈式分配策略,因此只有很短時段開啟3號和4號電解槽;策略3和策略4的電解槽功率分配基本相同。在功率較高的1 100~1 200 min范圍內,由于控制策略1、3、4的控制信號相同(每臺電解槽的功率和溫度指標相同),因此其效率也相同??梢钥闯?,采用功率-溫度自適應控制策略的電解槽,其效率在制氫高功率時段與其他三種策略相同,而在低、中制氫功率范圍內,制氫效率明顯提升。

    圖16 場景2下四種策略電解槽功率分配情況

    圖17 場景2下四種策略制氫效率

    圖18和圖19分別給出了場景3下四種制氫策略的功率分配方案及制氫效率曲線。由于制氫功率全時段都很高,控制策略1、3、4的功率分配結果相同;而鏈式分配策略難以使電解槽功率相同,因此在400~500 min功率高峰范圍內啟動了4號電解槽。在電解效率方面,控制策略1和3的功率和溫度指標相同,因此兩條效率曲線重合;由于控制策略4采用了溫度自適應控制,部分時段的溫度指標與策略1和3不同,效率也略高于其他兩種控制策略。

    圖18 場景3下四種策略電解槽功率分配情況

    圖19 場景3下四種策略制氫效率

    表3給出了三種典型場景下四種控制策略的制氫總量。相比于傳統(tǒng)的平均分配策略,功率-溫度自適應控制在三種場景下可以分別提升制氫量約9.8%、1.6%、0.9%;相對于鏈式分配策略,可分別提升6.2%、5.7%、4.1%??紤]到三個場景的概率分別為0.6、0.2、0.2,功率-溫度自適應控制策略相對于平均分配策略增加制氫量6.4%,相對于鏈式分配策略增加制氫量5.7%。

    表3 不同場景/策略下的制氫總量

    Tab.3 Total hydrogen production under different scenarios and strategies(單位:kg)

    4 結論

    本文建立了包含主要能耗環(huán)節(jié)的電解制氫系統(tǒng)效率模型,分析了電解制氫系統(tǒng)最優(yōu)效率與電解功率、溫度的關系,提出了基于離線優(yōu)化與在線控制的多堆PEM電解制氫系統(tǒng)效率優(yōu)化策略。通過算例分析,得到以下結論:

    1)制氫效率受到電解功率和溫度的雙重影響,升高電解溫度并不一定會提高制氫效率,其最優(yōu)取值需在制氫效率曲面上依據電解功率對應確定。

    2)相對于平均分配策略和鏈式分配策略,本文提出的功率-溫度自適應控制策略既可以充分利用低功率高效制氫,又可以在較大功率范圍內保持高效率,因此適用于制氫功率波動的場景。

    3)在實際風電制氫場景下,相對于傳統(tǒng)的平均控制策略和鏈式分配策略,本文提出的功率-溫度自適應控制策略分別提升產氫量6.4%和5.7%。

    由不同場景下的制氫效率曲線可知,優(yōu)化后的系統(tǒng)效率維持在70%左右,而大量的能量在電解、壓縮等環(huán)節(jié)以熱能形式耗散。因此,下一步計劃開展電解制氫系統(tǒng)與區(qū)域供熱系統(tǒng)耦合運行的相關研究,實現氫/熱聯(lián)供,提升可再生能源利用效率。

    以兩臺電解槽為例證明,其功率分配示意圖如附圖1所示。若兩臺電解槽分配的功率不相同,設其分別為1和2,對應的制氫量分別為(1)和(2),設(4)為兩臺電解槽的平均產氫量,即

    附圖1 兩臺電解槽功率分配示意圖

    App.Fig.1 Power distribution of 2 electrolytic cells

    令3為平均分配的功率,則有

    即對應的產氫量為2(3)。根據制氫量-功率曲線,有

    因此,兩臺電解槽功率平均分配的產氫量大于不平均分配的產氫量。類似地,如果啟動臺電解槽,則平均分配的產氫量大于其他分配策略。

    [1] 王雨晴, 王文詩, 徐心竹, 等. 面向低碳交通的含新能源汽車共享站電-氫微能源網區(qū)間-隨機混合規(guī)劃方法[J]. 電工技術學報, 2023, 38(23): 6373-6390. Wang Yuqing, Wang Wenshi, Xu Xinzhu, et al. Hybrid interval/stochastic planning method for new energy vehicle sharing station-based electro-hydrogen micro-energy system for low-carbon transportation[J]. Transactions of China Electrotechnical Society, 2023, 38(23): 6373-6390.

    [2] 吳孟雪, 房方. 計及風光不確定性的電-熱-氫綜合能源系統(tǒng)分布魯棒優(yōu)化[J]. 電工技術學報, 2023, 38(13): 3473-3485. Wu Mengxue, Fang Fang. Distributionally robust optimization of electricity-heat-hydrogen integrated energy system with wind and solar uncertainties[J]. Transactions of China Electrotechnical Society, 2023, 38(13): 3473-3485.

    [3] 劉瑋, 萬燕鳴, 熊亞林, 等. 碳中和目標下電解水制氫關鍵技術及價格平準化分析[J]. 電工技術學報, 2022, 37(11): 2888-2896. Liu Wei, Wan Yanming, Xiong Yalin, et al. Key technology of water electrolysis and levelized cost of hydrogen analysis under carbon neutral vision[J]. Transactions of China Electrotechnical Society, 2022, 37(11): 2888-2896.

    [4] Han Bo, Steen S M, Mo Jingke, et al. Electrochemical performance modeling of a proton exchange membrane electrolyzer cell for hydrogen energy[J]. International Journal of Hydrogen Energy, 2015, 40(22): 7006-7016.

    [5] 李軍舟, 趙晉斌, 曾志偉, 等. 具有動態(tài)調節(jié)特性的光伏制氫雙陣列直接耦合系統(tǒng)優(yōu)化策略[J]. 電網技術, 2022, 46(5): 1712-1721. Li Junzhou, Zhao Jinbin, Zeng Zhiwei, et al. Optimization strategy of photovoltaic hydrogen production dual array direct coupling system with dynamic regulation characteristics[J]. Power System Technology, 2022, 46(5): 1712-1721.

    [6] 熊宇峰, 司楊, 鄭天文, 等. 基于主從博弈的工業(yè)園區(qū)綜合能源系統(tǒng)氫儲能優(yōu)化配置[J]. 電工技術學報, 2021, 36(3): 507-516. Xiong Yufeng, Si Yang, Zheng Tianwen, et al. Optimal configuration of hydrogen storage in industrial park integrated energy system based on stackelberg game[J]. Transactions of China Electro-technical Society, 2021, 36(3): 507-516.

    [7] 魯明芳, 李咸善, 李飛, 等. 季節(jié)性氫儲能-混氫燃氣輪機系統(tǒng)兩階段隨機規(guī)劃[J]. 中國電機工程學報, 2023, 43(18): 6978-6992. Lu Mingfang, Li Xianshan, Li Fei, et al. Two-stage stochastic programming of seasonal hydrogen energy storage and mixed hydrogen-fueled gas turbine system[J]. Proceedings of the CSEE, 2023, 43(18): 6978-6992.

    [8] Huang Chunjun, Zong Yi, You Shi, et al. Economic model predictive control for multi-energy system considering hydrogen-thermal-electric dynamics and waste heat recovery of MW-level alkaline electro-lyzer[J]. Energy Conversion and Management, 2022, 265: 115697.

    [9] Nguyen T T, Kim H M. Cluster-based predictive PCC voltage control of large-scale offshore wind farm[J]. IEEE Access, 2020, 9: 4630-4641.

    [10] Kumar S S, Nirmal Mukundan C M, Jayaprakash P. Modified LMS control for a grid interactive PV–fuel cell–electrolyzer hybrid system with power dispatch to the grid[J]. IEEE Transactions on Industry Applications, 2022, 58(6): 7907-7918.

    [11] Abdelghany M B, Faisal Shehzad M, Liuzza D, et al. Modeling and optimal control of a hydrogen storage system for wind farm output power smoothing[C]// 2020 59th IEEE Conference on Decision and Control (CDC), Jeju, Korea (South), 2020: 49-54.

    [12] Huang Chunjun, Zong Yi, You Shi, et al. Cooperative control of wind-hydrogen-SMES hybrid systems for fault-ride-through improvement and power smoothing[J]. IEEE Transactions on Applied Superconductivity, 2021, 31(8): 1-7.

    [13] Guerra O J, Eichman J, Kurtz J, et al. Cost competitiveness of electrolytic hydrogen[J]. Joule, 2019, 3(10): 2425-2443.

    [14] Tjarks G, Gibelhaus A, Lanzerath F, et al. Energetically-optimal PEM electrolyzer pressure in power-to-gas plants[J]. Applied Energy, 2018, 218: 192-198.

    [15] Correa G, Marocco P, Mu?oz P, et al. Pressurized PEM water electrolysis: dynamic modelling focusing on the cathode side[J]. International Journal of Hydrogen Energy, 2022, 47(7): 4315-4327.

    [16] Scheepers F, St?hler M, St?hler A, et al. Improving the efficiency of PEM electrolyzers through membrane-specific pressure optimization[J]. Energies, 2020, 13(3): 612.

    [17] Scheepers F, St?hler M, St?hler A, et al. Temperature optimization for improving polymer electrolyte membrane-water electrolysis system efficiency[J]. Applied Energy, 2021, 283: 116270.

    [18] 鄧智宏, 江岳文. 考慮制氫效率特性的風氫系統(tǒng)容量優(yōu)化[J]. 可再生能源, 2020, 38(2): 259-266. Deng Zhihong, Jiang Yuewen. Optimal sizing of a wind-hydrogen system under consideration of the efficiency characteristics of electrolysers[J]. Renewable Energy Resources, 2020, 38(2): 259-266.

    [19] Keller R, Rauls E, Hehemann M, et al. An adaptive model-based feedforward temperature control of a 100 kW PEM electrolyzer[J]. Control Engineering Practice, 2022, 120: 104992.

    [20] 沈小軍, 聶聰穎, 呂洪. 計及電熱特性的離網型風電制氫堿性電解槽陣列優(yōu)化控制策略[J]. 電工技術學報, 2021, 36(3): 463-472. Shen Xiaojun, Nie Congying, Lü Hong. Coordination control strategy of wind power-hydrogen alkaline electrolyzer bank considering electrothermal charac-teristics[J]. Transactions of China Electrotechnical Society, 2021, 36(3): 463-472.

    [21] 唐司航. 基于風電功率最值預測的電解槽陣列輪換控制策略研究[D]. 石家莊: 河北科技大學, 2022. Tang Sihang. Research on control strategy of alkaline electrolyzer array with rotation mode based on best value forecast for wind power[D]. Shijiazhuang: Hebei University of Science and Technology, 2022.

    [22] 盧昕宇, 杜幫華, 趙波, 等. 基于鏈式分配策略的風氫耦合系統(tǒng)設計與控制[J]. 太陽能學報, 2022, 43(6): 405-413. Lu Xinyu, Du Banghua, Zhao Bo, et al. Design and control of wind-hydrogen coupled system based on chain distribution strategy[J]. Acta Energiae Solaris Sinica, 2022, 43(6): 405-413.

    [23] Luxa A, J?res N, Yá?ez C, et al. Multilinear modeling and simulation of a multi-stack PEM electrolyzer with degradation for control concept comparison[C]// Proceedings of the 12th International Conference on Simulation and Modeling Methodologies, Technologies and Applications, Lisbon, Portugal, 2022: 52-62.

    [24] Barei? K, de la Rua C, M?ckl M, et al. Life cycle assessment of hydrogen from proton exchange membrane water electrolysis in future energy systems[J]. Applied Energy, 2019, 237: 862-872.

    [25] Sakas G, Ibá?ez-Rioja A, Ruuskanen V, et al. Dynamic energy and mass balance model for an industrial alkaline water electrolyzer plant process[J]. International Journal of Hydrogen Energy, 2022, 47(7): 4328-4345.

    [26] Tiktak W J. Heat management of PEM electrolysis[D]. Delft: Delft University of Technology, 2019.

    [27] Zheng Yi, You Shi, Bindner H W, et al. Optimal day-ahead dispatch of an alkaline electrolyser system concerning thermal–electric properties and state-transitional dynamics[J]. Applied Energy, 2022, 307: 118091.

    [28] 游雙矯. 中石油為雄安新區(qū)供氫方式優(yōu)化的研究[D]. 北京: 中國石油大學(北京), 2020. Yon Shuangjiao. Research on optimization of hydrogen supplies for Xiongan New Area by CNPC[D]. Beijing: China University of Petroleum, 2020.

    [29] 徐衛(wèi)君, 張偉, 胡宇濤, 等. 先進絕熱壓縮空氣儲能多能流優(yōu)化調度模型[J]. 電工技術學報, 2022, 37(23): 5944-5955. Xu Weijun, Zhang Wei, Hu Yutao, et al. Multi energy flow optimal scheduling model of advanced adiabatic compressed air energy storage[J]. Transactions of China Electrotechnical Society, 2022, 37(23): 5944-5955.

    [30] Skorek-Osikowska A, Bartela ?, Katla D, et al. Thermodynamic assessment of the novel concept of the energy storage system using compressed carbon dioxide, methanation and hydrogen generator[J]. Fuel, 2021, 304: 120764.

    [31] 張理, 韓民曉, 范溢文. 多相堆疊交錯并聯(lián)制氫變換器控制策略與特性分析[J]. 電工技術學報, 2023, 38(2): 485-495. Zhang Li, Han Minxiao, Fan Yiwen. Control strategy and characteristic analysis of multi-phase stacked interleaved buck converter for hydrogen production[J]. Transactions of China Electrotechnical Society, 2023, 38(2): 485-495.

    [32] 任洲洋, 王皓, 李文沅, 等. 基于氫能設備多狀態(tài)模型的電氫區(qū)域綜合能源系統(tǒng)可靠性評估[J]. 電工技術學報, 2023, 38(24): 6744-6759. Ren Zhouyang, Wang Hao, Li Wenyuan, et al. Reliability evaluation of electricity-hydrogen regional integrated energy systems based on the multi-state models of hydrogen energy equipment[J]. Transactions of China Electrotechnical Society, 2023, 38(24): 6744-6759.

    [33] Ayachit A, Kazimierczuk M K. Averaged small-signal model of PWM DC-DC converters in CCM including switching power loss[J]. IEEE Transactions on Circuits and Systems II: Express Briefs, 2019, 66(2): 262-266.

    [34] Talebian I, Alavi P, Marzang V, et al. Analysis, design, and investigation of a soft-switched buck converter with high efficiency[J]. IEEE Transactions on Power Electronics, 2022, 37(6): 6899-6912.

    [35] Yigit T, Selamet O F. Mathematical modeling and dynamic Simulink simulation of high-pressure PEM electrolyzer system[J]. International Journal of Hydrogen Energy, 2016, 41(32): 13901-13914.

    [36] Rizwan M, Alstad V, J?schke J. Design considerations for industrial water electrolyzer plants[J]. International Journal of Hydrogen Energy, 2021, 46(75): 37120-37136.

    [37] 宋天昊, 李柯江, 韓肖清, 等. 儲能系統(tǒng)參與多應用場景的協(xié)同運行策略[J]. 電力系統(tǒng)自動化, 2021, 45(19): 43-51. Song Tianhao, Li Kejiang, Han Xiaoqing, et al. Coordinated operation strategy of energy storage system participating in multiple application scenarios[J]. Automation of Electric Power Systems, 2021, 45(19): 43-51.

    Operational Efficiency Enhancement of Multi-Stack Proton Exchange Membrane Electrolyzer Systems with Power-Temperature Adaptive Control

    Han Pengfei1,2Xu Xiaoyuan1,2Wang Han1,2Yan Zheng1,2

    (1. Key Laboratory of Control of Power Transmission and Conversion Shanghai Jiao Tong University Shanghai 200240 China 2. Shanghai Non-Carbon Energy Conversion and Utilization Institute Shanghai 200240 China)

    With the growing focus of human society on low-carbon energy, the proportion of hydrogen energy in power systems is increasing and the hydrogen production system from water electrolysis is evolving towards a larger capacity and multiple stacks. Water electrolysis based on the proton exchange membrane (PEM) technology using renewable energy is a promising means of promoting the consumption of renewable energy and building a low-carbon power system. However, the randomness and volatility of renewable energy pose challenges to the efficient operation of a multi-stack PEM electrolytic hydrogen production system under fluctuating power scenarios. To enhance the efficiency of hydrogen production systems with multi-stack PEM electrolyzers, this research proposes an efficiency optimization strategy with power-temperature adaptive control. The research work is described as follows:

    Firstly, an overall efficiency model of the hydrogen production system is established. The relationship between optimal system efficiency and the electrolytic power and temperature is analyzed. It reveals that the optimal electrolytic efficiency might not be achieved at the upper limits of electrolytic temperature. When the current density is low, the optimal electrolytic temperature is relatively low due to the significant reduction in Faraday efficiency caused by temperature rise. However, when the current density exceeds 1.5 A/cm2, the optimal electrolysis temperature increases with the increase of current density. Therefore, to improve the electrolytic efficiency when the power of renewable energy fluctuates, it is necessary to properly adjust the electrolytic power and temperature of the electrolyzers.

    Secondly, offline optimization and online control are employed to enhance electrolytic efficiency. The offline optimization involves a two-stage optimization model, where the first stage optimizes the startup and shutdown status, electrolysis power, and temperature settings of different electrolyzers to maximize hydrogen production efficiency. In the second stage, a specific power-temperature implementation scheme is determined for the electrolyzers, to achieve rapid temperature adjustment when there are multiple solutions in the first stage. In the online control section, the dynamic models of the power and temperature control of the electrolyzers are established, respectively. Based on the current total power of hydrogen production and the power-temperature setting scheme determined in the offline optimization stage, the PWM modulation signal and the valve opening signal are generated by the buck and the valve controllers, respectively.

    Finally, to verify the effectiveness of the proposed efficiency enhancement strategy, a comparison is made between the proposed method and the existing average allocation strategy and chain allocation strategy. The analysis shows that the efficiency of hydrogen production is affected by both electrolytic power and temperature. Increasing the electrolysis temperature does not necessarily improve the hydrogen production efficiency, and its optimal value should be determined based on the electrolysis power. Compared with the average allocation strategy and the chain allocation strategy, the power-temperature adaptive control strategy proposed in this research can make full use of low renewable energy efficiently, while maintaining high efficiency over a large power range, making it suitable for hydrogen production with power fluctuations. In actual scenarios of hydrogen production from wind power, the proposed power-temperature adaptive control strategy can increase hydrogen production by 6.4% and 5.7%, respectively, compared with the average allocation strategy and chain allocation strategy.

    Proton exchange membrane elecrolyzers, power to hydrogen, renewable energy power generation, efficiency optimization, adaptive control

    10.19595/j.cnki.1000-6753.tces.230113

    TQ116.2; TM73

    國家自然科學基金資助項目(52077136,U2166201)。

    2023-02-01

    2023-03-13

    韓鵬飛 男,1993年生,博士研究生,研究方向為可再生能源發(fā)電制氫系統(tǒng)優(yōu)化與控制等。E-mail:han_pf@sjtu.edu.cn

    徐瀟源 男,1989年生,副教授,研究方向為電力系統(tǒng)不確定性分析、電力系統(tǒng)優(yōu)化運行等。E-mail:xuxiaoyuan@sjtu.edu.cn(通信作者)

    (編輯 李 冰)

    猜你喜歡
    適應控制電解槽制氫
    堿性電解槽成本最低
    采用自適應控制的STATCOM-SEIG系統(tǒng)Matlab仿真
    電子測試(2018年15期)2018-09-26 06:01:04
    電解槽焊接施工中的質量控制
    江西建材(2018年1期)2018-04-04 05:26:02
    考慮執(zhí)行器飽和的改進無模型自適應控制
    自動化學報(2016年8期)2016-04-16 03:38:51
    制氫工藝技術比較
    基于 L1自適應控制的無人機橫側向控制
    高活性Al-LiBH4-Bi鋁基復合制氫材料
    電源技術(2015年11期)2015-08-22 08:50:26
    Rucklidge混沌系統(tǒng)的自適應控制
    電解制氫設備開發(fā)入選“863”
    低溫與特氣(2014年4期)2014-03-20 13:36:50
    零極距電解槽改造的幾點思考
    中國氯堿(2014年12期)2014-02-28 01:05:11
    麻豆乱淫一区二区| 人妻 亚洲 视频| 国产又色又爽无遮挡免| 日韩av在线免费看完整版不卡| 久久这里只有精品19| 啦啦啦中文免费视频观看日本| 曰老女人黄片| tube8黄色片| 91麻豆精品激情在线观看国产 | 国产色视频综合| www.精华液| av不卡在线播放| 亚洲国产精品国产精品| 色网站视频免费| 99香蕉大伊视频| 黄色怎么调成土黄色| 狂野欧美激情性bbbbbb| 午夜福利一区二区在线看| 国产成人精品久久二区二区91| 操美女的视频在线观看| 午夜福利影视在线免费观看| www.自偷自拍.com| 欧美激情极品国产一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| 少妇猛男粗大的猛烈进出视频| 国产精品久久久久久精品古装| 午夜福利一区二区在线看| 欧美在线一区亚洲| 国产片内射在线| 久久久国产精品麻豆| 成人手机av| 欧美日韩亚洲国产一区二区在线观看 | 久久国产精品人妻蜜桃| 蜜桃国产av成人99| 亚洲成色77777| 精品熟女少妇八av免费久了| 亚洲中文日韩欧美视频| 精品人妻一区二区三区麻豆| 久久人妻福利社区极品人妻图片 | 国产淫语在线视频| 天天躁夜夜躁狠狠久久av| 亚洲精品国产一区二区精华液| 色婷婷久久久亚洲欧美| 涩涩av久久男人的天堂| 赤兔流量卡办理| 天天添夜夜摸| 久久精品亚洲av国产电影网| 色精品久久人妻99蜜桃| 亚洲美女黄色视频免费看| 免费不卡黄色视频| 久久ye,这里只有精品| 一区二区三区四区激情视频| 老司机亚洲免费影院| 丰满饥渴人妻一区二区三| 丝瓜视频免费看黄片| 精品国产一区二区三区四区第35| 亚洲男人天堂网一区| 中文字幕人妻丝袜制服| 亚洲人成电影观看| 亚洲五月婷婷丁香| 亚洲av综合色区一区| 侵犯人妻中文字幕一二三四区| 黄色 视频免费看| 久久国产精品男人的天堂亚洲| 精品国产一区二区三区四区第35| 在线观看免费视频网站a站| 亚洲 国产 在线| 丝袜在线中文字幕| 丝袜喷水一区| 91麻豆精品激情在线观看国产 | 手机成人av网站| 国产深夜福利视频在线观看| 亚洲国产看品久久| 美女高潮到喷水免费观看| 视频在线观看一区二区三区| 国产极品粉嫩免费观看在线| 男女国产视频网站| 国产精品三级大全| 国产亚洲欧美精品永久| 国产亚洲精品第一综合不卡| 国产成人系列免费观看| 亚洲精品日韩在线中文字幕| 精品第一国产精品| 少妇粗大呻吟视频| 五月开心婷婷网| 免费高清在线观看视频在线观看| 日韩一卡2卡3卡4卡2021年| 午夜精品国产一区二区电影| 九色亚洲精品在线播放| 亚洲欧美日韩另类电影网站| 丰满饥渴人妻一区二区三| 99久久人妻综合| 1024香蕉在线观看| 99国产综合亚洲精品| 久久久亚洲精品成人影院| 日韩一卡2卡3卡4卡2021年| 久久热在线av| 午夜福利影视在线免费观看| 久久久久久免费高清国产稀缺| 国产xxxxx性猛交| 成人免费观看视频高清| 亚洲av国产av综合av卡| 女人精品久久久久毛片| 精品一区二区三区av网在线观看 | 亚洲伊人色综图| 国产精品av久久久久免费| 国产男人的电影天堂91| 狂野欧美激情性bbbbbb| 欧美日本中文国产一区发布| 国产成人91sexporn| 老司机深夜福利视频在线观看 | 国产精品久久久久久人妻精品电影 | www.av在线官网国产| 美女主播在线视频| 欧美人与性动交α欧美软件| h视频一区二区三区| 亚洲人成电影观看| 看免费av毛片| 免费少妇av软件| 国产成人精品久久二区二区91| 一级毛片电影观看| 成人国产一区最新在线观看 | 又紧又爽又黄一区二区| 午夜福利,免费看| 一区在线观看完整版| 亚洲av片天天在线观看| 一区二区日韩欧美中文字幕| 国产熟女午夜一区二区三区| 又黄又粗又硬又大视频| 亚洲,欧美,日韩| 视频区图区小说| 中文字幕精品免费在线观看视频| 亚洲人成网站在线观看播放| 菩萨蛮人人尽说江南好唐韦庄| 日韩免费高清中文字幕av| 日韩 欧美 亚洲 中文字幕| 又大又爽又粗| 精品久久久精品久久久| 亚洲,欧美精品.| 亚洲综合色网址| 欧美精品亚洲一区二区| 日韩欧美一区视频在线观看| 少妇人妻 视频| 老司机亚洲免费影院| 热99久久久久精品小说推荐| 操美女的视频在线观看| 久久久精品区二区三区| 日本黄色日本黄色录像| 国产男女内射视频| 超碰97精品在线观看| 国产成人精品在线电影| 爱豆传媒免费全集在线观看| 999精品在线视频| 国产欧美亚洲国产| 精品欧美一区二区三区在线| 蜜桃国产av成人99| cao死你这个sao货| 中文字幕av电影在线播放| 丁香六月天网| 只有这里有精品99| 免费一级毛片在线播放高清视频 | 韩国高清视频一区二区三区| 亚洲九九香蕉| 久久久久久久精品精品| 亚洲国产看品久久| 午夜老司机福利片| 国产亚洲av高清不卡| 女人久久www免费人成看片| 国产精品人妻久久久影院| 水蜜桃什么品种好| 亚洲中文av在线| 国产男人的电影天堂91| 蜜桃国产av成人99| 伊人亚洲综合成人网| 脱女人内裤的视频| 国产av一区二区精品久久| 99精品久久久久人妻精品| 久久热在线av| 久久99热这里只频精品6学生| 国产成人免费观看mmmm| 精品一区二区三区av网在线观看 | 欧美性长视频在线观看| 精品少妇一区二区三区视频日本电影| 少妇粗大呻吟视频| 精品视频人人做人人爽| 一本一本久久a久久精品综合妖精| 国产精品av久久久久免费| 国产精品欧美亚洲77777| 在线亚洲精品国产二区图片欧美| 日韩人妻精品一区2区三区| 啦啦啦在线免费观看视频4| 午夜影院在线不卡| 青草久久国产| 在线 av 中文字幕| 岛国毛片在线播放| av一本久久久久| 狠狠精品人妻久久久久久综合| 十八禁网站网址无遮挡| 亚洲成色77777| 秋霞在线观看毛片| 悠悠久久av| 亚洲国产精品999| 亚洲欧美成人综合另类久久久| 一边摸一边抽搐一进一出视频| 日本黄色日本黄色录像| 国产高清videossex| 国产爽快片一区二区三区| 国产女主播在线喷水免费视频网站| 91老司机精品| 女人高潮潮喷娇喘18禁视频| 久9热在线精品视频| 欧美精品啪啪一区二区三区 | 久久天堂一区二区三区四区| 国产97色在线日韩免费| 妹子高潮喷水视频| 你懂的网址亚洲精品在线观看| 国产一区二区 视频在线| 在线观看人妻少妇| 99热国产这里只有精品6| 精品福利观看| 亚洲九九香蕉| 一区二区三区四区激情视频| 热99久久久久精品小说推荐| 中文精品一卡2卡3卡4更新| 超碰97精品在线观看| 欧美xxⅹ黑人| 国产99久久九九免费精品| 男人添女人高潮全过程视频| e午夜精品久久久久久久| 天堂俺去俺来也www色官网| 肉色欧美久久久久久久蜜桃| 青青草视频在线视频观看| 欧美精品高潮呻吟av久久| 亚洲天堂av无毛| 亚洲一区中文字幕在线| 精品国产乱码久久久久久男人| 国产精品一区二区免费欧美 | 最近中文字幕2019免费版| 成人18禁高潮啪啪吃奶动态图| 亚洲国产精品国产精品| 精品国产乱码久久久久久小说| 国产精品免费视频内射| 99国产综合亚洲精品| 天堂8中文在线网| 国产免费又黄又爽又色| 操美女的视频在线观看| 久久人人97超碰香蕉20202| 日韩av在线免费看完整版不卡| 91九色精品人成在线观看| 亚洲第一青青草原| svipshipincom国产片| 涩涩av久久男人的天堂| 日韩伦理黄色片| 日韩 亚洲 欧美在线| 老司机影院毛片| 国产精品一区二区免费欧美 | 久久人妻福利社区极品人妻图片 | 免费观看人在逋| 高潮久久久久久久久久久不卡| 日本欧美视频一区| 国产成人一区二区三区免费视频网站 | 国产一区有黄有色的免费视频| 久久九九热精品免费| 久久久久网色| 美女脱内裤让男人舔精品视频| 久久午夜综合久久蜜桃| 80岁老熟妇乱子伦牲交| 国产91精品成人一区二区三区 | h视频一区二区三区| 如日韩欧美国产精品一区二区三区| 亚洲午夜精品一区,二区,三区| 亚洲成色77777| 青青草视频在线视频观看| 久久精品熟女亚洲av麻豆精品| 日韩制服骚丝袜av| 狠狠精品人妻久久久久久综合| 18在线观看网站| 欧美精品啪啪一区二区三区 | 大型av网站在线播放| 精品高清国产在线一区| 亚洲,欧美,日韩| 亚洲欧洲日产国产| 久久狼人影院| 中文字幕色久视频| 91字幕亚洲| 久久99精品国语久久久| 国产精品一区二区精品视频观看| 在线看a的网站| 色播在线永久视频| videosex国产| 欧美日韩视频高清一区二区三区二| h视频一区二区三区| 国产成人免费观看mmmm| 午夜免费成人在线视频| 少妇人妻 视频| 老汉色∧v一级毛片| 日韩人妻精品一区2区三区| 免费在线观看视频国产中文字幕亚洲 | 欧美日韩亚洲综合一区二区三区_| 亚洲色图综合在线观看| 免费在线观看完整版高清| av网站在线播放免费| 国产野战对白在线观看| 99热网站在线观看| 精品人妻一区二区三区麻豆| 老司机午夜十八禁免费视频| 亚洲av在线观看美女高潮| 欧美人与性动交α欧美精品济南到| 麻豆乱淫一区二区| 亚洲精品久久久久久婷婷小说| 手机成人av网站| 黄频高清免费视频| 亚洲av电影在线观看一区二区三区| 亚洲国产成人一精品久久久| 另类精品久久| 91字幕亚洲| 黄片小视频在线播放| 国产欧美亚洲国产| 婷婷色综合大香蕉| 校园人妻丝袜中文字幕| 亚洲av综合色区一区| 国产在线视频一区二区| 国产成人91sexporn| 欧美 亚洲 国产 日韩一| 美女中出高潮动态图| 男女边摸边吃奶| 久久ye,这里只有精品| 婷婷色av中文字幕| 波野结衣二区三区在线| 亚洲av男天堂| 欧美精品啪啪一区二区三区 | 99精品久久久久人妻精品| 少妇猛男粗大的猛烈进出视频| 精品国产一区二区三区四区第35| 日韩av不卡免费在线播放| 欧美大码av| 亚洲精品国产av蜜桃| 欧美日韩亚洲综合一区二区三区_| 考比视频在线观看| videos熟女内射| 大话2 男鬼变身卡| 欧美另类一区| 热re99久久精品国产66热6| 国产精品九九99| 两人在一起打扑克的视频| 老司机深夜福利视频在线观看 | 热re99久久国产66热| 97在线人人人人妻| 一级片免费观看大全| 黄色怎么调成土黄色| 韩国精品一区二区三区| 最近中文字幕2019免费版| 午夜久久久在线观看| 午夜两性在线视频| 国产精品99久久99久久久不卡| 国产又色又爽无遮挡免| 别揉我奶头~嗯~啊~动态视频 | 午夜福利影视在线免费观看| 久久精品人人爽人人爽视色| 精品第一国产精品| 宅男免费午夜| 亚洲精品国产色婷婷电影| 另类亚洲欧美激情| 又粗又硬又长又爽又黄的视频| 亚洲av在线观看美女高潮| 日本色播在线视频| 亚洲三区欧美一区| 日韩一本色道免费dvd| 欧美 亚洲 国产 日韩一| 波野结衣二区三区在线| 午夜免费鲁丝| 亚洲精品一二三| 啦啦啦在线观看免费高清www| 女人高潮潮喷娇喘18禁视频| 国产欧美日韩综合在线一区二区| 老汉色∧v一级毛片| 日韩大码丰满熟妇| 少妇 在线观看| 亚洲av欧美aⅴ国产| 午夜福利影视在线免费观看| 90打野战视频偷拍视频| 日日夜夜操网爽| 高清视频免费观看一区二区| av天堂久久9| 一区二区日韩欧美中文字幕| 精品熟女少妇八av免费久了| 国产欧美日韩一区二区三区在线| 亚洲一卡2卡3卡4卡5卡精品中文| 老司机靠b影院| 这个男人来自地球电影免费观看| 五月天丁香电影| 最近中文字幕2019免费版| 亚洲,一卡二卡三卡| 看免费成人av毛片| 咕卡用的链子| 大型av网站在线播放| 欧美日韩av久久| 亚洲一区二区三区欧美精品| 亚洲精品中文字幕在线视频| 一区二区av电影网| 国产亚洲精品第一综合不卡| 欧美日韩亚洲国产一区二区在线观看 | 亚洲五月色婷婷综合| 亚洲人成77777在线视频| 色婷婷av一区二区三区视频| 妹子高潮喷水视频| 麻豆av在线久日| 亚洲精品国产av成人精品| 91麻豆精品激情在线观看国产 | 你懂的网址亚洲精品在线观看| 精品一区二区三区av网在线观看 | 亚洲人成网站在线观看播放| 久久人人爽av亚洲精品天堂| 久久久久久免费高清国产稀缺| 色网站视频免费| 国产高清国产精品国产三级| 国产日韩欧美亚洲二区| 亚洲激情五月婷婷啪啪| 国产视频首页在线观看| 午夜视频精品福利| 伊人亚洲综合成人网| 脱女人内裤的视频| 精品亚洲成a人片在线观看| 99国产精品免费福利视频| 最新在线观看一区二区三区 | 国产精品 欧美亚洲| 欧美 日韩 精品 国产| 男女免费视频国产| 免费av中文字幕在线| 日韩制服丝袜自拍偷拍| 大型av网站在线播放| 啦啦啦在线免费观看视频4| 十八禁高潮呻吟视频| 精品久久久久久久毛片微露脸 | 免费在线观看视频国产中文字幕亚洲 | 日韩一本色道免费dvd| 9色porny在线观看| 亚洲专区中文字幕在线| 欧美xxⅹ黑人| 国产成人啪精品午夜网站| 老司机靠b影院| 大话2 男鬼变身卡| 亚洲情色 制服丝袜| 人人妻,人人澡人人爽秒播 | e午夜精品久久久久久久| 久热爱精品视频在线9| 精品久久久久久久毛片微露脸 | 精品免费久久久久久久清纯 | 久久人人爽av亚洲精品天堂| 超碰成人久久| 精品亚洲成国产av| 久久人妻福利社区极品人妻图片 | 久久综合国产亚洲精品| 亚洲成人手机| 久久久久视频综合| 日韩中文字幕欧美一区二区 | 久久九九热精品免费| 婷婷色综合www| 国产一区二区三区综合在线观看| 亚洲精品国产一区二区精华液| 午夜视频精品福利| 中文字幕另类日韩欧美亚洲嫩草| 人人澡人人妻人| av天堂在线播放| 亚洲精品自拍成人| 一级毛片我不卡| 亚洲第一青青草原| 亚洲伊人久久精品综合| 久久久久国产精品人妻一区二区| 久久精品人人爽人人爽视色| 国产成人啪精品午夜网站| 人人妻,人人澡人人爽秒播 | 亚洲伊人久久精品综合| 一个人免费看片子| 大片电影免费在线观看免费| 午夜福利一区二区在线看| av又黄又爽大尺度在线免费看| av在线老鸭窝| 桃花免费在线播放| 久久狼人影院| 国产精品一区二区在线不卡| 亚洲精品美女久久av网站| 午夜福利视频精品| 99精国产麻豆久久婷婷| 五月天丁香电影| 久久免费观看电影| 成年人午夜在线观看视频| 黄频高清免费视频| 国语对白做爰xxxⅹ性视频网站| 岛国毛片在线播放| 精品久久久精品久久久| 满18在线观看网站| 99热全是精品| 国产av精品麻豆| 天堂8中文在线网| 国产精品麻豆人妻色哟哟久久| 国产免费福利视频在线观看| 啦啦啦在线免费观看视频4| 国产野战对白在线观看| 美女主播在线视频| 一级黄片播放器| 深夜精品福利| av网站在线播放免费| 妹子高潮喷水视频| 午夜福利免费观看在线| 久久久精品94久久精品| 久久久久久久国产电影| 十八禁人妻一区二区| 亚洲欧美一区二区三区国产| 只有这里有精品99| 亚洲专区中文字幕在线| 色婷婷av一区二区三区视频| 一个人免费看片子| 国产精品三级大全| 亚洲,欧美,日韩| 脱女人内裤的视频| 日韩熟女老妇一区二区性免费视频| 欧美日韩国产mv在线观看视频| 黄频高清免费视频| a 毛片基地| 又大又爽又粗| 免费在线观看影片大全网站 | 亚洲伊人久久精品综合| 在线观看人妻少妇| 9色porny在线观看| 黄频高清免费视频| 美女视频免费永久观看网站| 各种免费的搞黄视频| 欧美黑人欧美精品刺激| 国产在视频线精品| 国产高清videossex| 亚洲 国产 在线| 国产男女内射视频| 欧美激情极品国产一区二区三区| 亚洲国产精品成人久久小说| 成年人黄色毛片网站| 无遮挡黄片免费观看| 亚洲情色 制服丝袜| 国产亚洲精品久久久久5区| 另类精品久久| 精品欧美一区二区三区在线| 男女床上黄色一级片免费看| 不卡av一区二区三区| 亚洲天堂av无毛| 日韩 欧美 亚洲 中文字幕| 亚洲,一卡二卡三卡| 男的添女的下面高潮视频| 国产爽快片一区二区三区| 夫妻午夜视频| 久久国产精品大桥未久av| 国产精品国产三级国产专区5o| 七月丁香在线播放| 尾随美女入室| 欧美性长视频在线观看| 中文字幕人妻熟女乱码| 精品人妻在线不人妻| 亚洲国产欧美网| 欧美日韩成人在线一区二区| 天天躁狠狠躁夜夜躁狠狠躁| 久久99热这里只频精品6学生| 亚洲一码二码三码区别大吗| 天天躁夜夜躁狠狠躁躁| 伦理电影免费视频| 香蕉丝袜av| 国产无遮挡羞羞视频在线观看| 精品欧美一区二区三区在线| 岛国毛片在线播放| 久久国产精品影院| 首页视频小说图片口味搜索 | 国产精品一区二区免费欧美 | 国产男女内射视频| 久久女婷五月综合色啪小说| 天堂俺去俺来也www色官网| 日韩中文字幕视频在线看片| 黄色视频不卡| 侵犯人妻中文字幕一二三四区| 女人精品久久久久毛片| 国产又色又爽无遮挡免| 一本大道久久a久久精品| 欧美人与性动交α欧美精品济南到| avwww免费| 国产免费一区二区三区四区乱码| 91国产中文字幕| 美女中出高潮动态图| 99久久综合免费| 欧美日韩成人在线一区二区| 日本欧美视频一区| 中文字幕高清在线视频| 国产精品偷伦视频观看了| 欧美97在线视频| 99久久综合免费| 久久九九热精品免费| 国产精品香港三级国产av潘金莲 | 日韩人妻精品一区2区三区| 极品少妇高潮喷水抽搐| 两个人看的免费小视频| 天天操日日干夜夜撸| 97在线人人人人妻| 亚洲精品久久午夜乱码| 日韩av不卡免费在线播放| 精品亚洲成a人片在线观看| 看免费av毛片| 亚洲一区二区三区欧美精品| 国产精品麻豆人妻色哟哟久久| 国产成人免费无遮挡视频| 免费高清在线观看视频在线观看| 久久国产精品男人的天堂亚洲| 精品少妇内射三级| 日本av免费视频播放| 成人午夜精彩视频在线观看| 日韩免费高清中文字幕av| 欧美日韩一级在线毛片| 亚洲成人国产一区在线观看 | 69精品国产乱码久久久| h视频一区二区三区|