• <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
    日韩 欧美 亚洲 中文字幕| tocl精华| 99久国产av精品| 亚洲人成网站高清观看| 小说图片视频综合网站| 九九久久精品国产亚洲av麻豆 | 一个人观看的视频www高清免费观看 | 日韩有码中文字幕| 亚洲人成网站高清观看| 人妻夜夜爽99麻豆av| 亚洲国产看品久久| 91久久精品国产一区二区成人 | 99久久精品国产亚洲精品| 美女高潮的动态| 国内久久婷婷六月综合欲色啪| 脱女人内裤的视频| 亚洲自偷自拍图片 自拍| 全区人妻精品视频| 偷拍熟女少妇极品色| 他把我摸到了高潮在线观看| 丰满的人妻完整版| 婷婷精品国产亚洲av在线| 欧美成狂野欧美在线观看| 国产91精品成人一区二区三区| 成人无遮挡网站| 一区福利在线观看| 欧美乱妇无乱码| 欧美成人性av电影在线观看| 岛国在线免费视频观看| 亚洲欧美一区二区三区黑人| 亚洲国产欧美一区二区综合| 久久精品91无色码中文字幕| 91av网站免费观看| 变态另类丝袜制服| 国产私拍福利视频在线观看| 久久这里只有精品中国| 女生性感内裤真人,穿戴方法视频| 51午夜福利影视在线观看| 毛片女人毛片| 亚洲男人的天堂狠狠| 亚洲中文日韩欧美视频| 国产精品久久久久久久电影 | 激情在线观看视频在线高清| 不卡一级毛片| 日韩成人在线观看一区二区三区| 亚洲片人在线观看| 免费高清视频大片| 天天添夜夜摸| 国产黄色小视频在线观看| 国产成人精品无人区| 欧美一区二区国产精品久久精品| 又爽又黄无遮挡网站| 99国产极品粉嫩在线观看| 男人舔女人的私密视频| 美女午夜性视频免费| 久久久久久久久久黄片| 噜噜噜噜噜久久久久久91| 国产成人精品久久二区二区91| 亚洲五月婷婷丁香| 亚洲人成伊人成综合网2020| 亚洲avbb在线观看| 免费观看精品视频网站| 国产精品99久久久久久久久| 国产亚洲av高清不卡| 伊人久久大香线蕉亚洲五| 色噜噜av男人的天堂激情| АⅤ资源中文在线天堂| 黄频高清免费视频| 久久久国产成人免费| svipshipincom国产片| av黄色大香蕉| 亚洲中文日韩欧美视频| av福利片在线观看| 好男人电影高清在线观看| 亚洲国产看品久久| 麻豆国产97在线/欧美| 亚洲欧美一区二区三区黑人| 1024手机看黄色片| 精品久久久久久成人av| 怎么达到女性高潮| 手机成人av网站| 国产高潮美女av| 99re在线观看精品视频| 天堂√8在线中文| 91av网站免费观看| 真人一进一出gif抽搐免费| 国产69精品久久久久777片 | 亚洲乱码一区二区免费版| 亚洲激情在线av| 最近视频中文字幕2019在线8| 日本 欧美在线| 午夜激情欧美在线| 国产精品永久免费网站| 人人妻人人看人人澡| 一级黄色大片毛片| 国产av在哪里看| 国产91精品成人一区二区三区| 婷婷亚洲欧美| www日本在线高清视频| 国产极品精品免费视频能看的| 2021天堂中文幕一二区在线观| 他把我摸到了高潮在线观看| 国产精品自产拍在线观看55亚洲| 精品免费久久久久久久清纯| 黄片大片在线免费观看| xxxwww97欧美| 啦啦啦观看免费观看视频高清| 亚洲va日本ⅴa欧美va伊人久久| 国内精品美女久久久久久| 亚洲欧美日韩无卡精品| 亚洲18禁久久av| 国产高清有码在线观看视频| 久9热在线精品视频| 成人三级做爰电影| 免费一级毛片在线播放高清视频| 变态另类成人亚洲欧美熟女| av天堂在线播放| 久久中文看片网| 18禁美女被吸乳视频| 亚洲欧美日韩高清专用| 精品久久久久久成人av| 一级毛片高清免费大全| 国产精品1区2区在线观看.| 国产精品影院久久| 亚洲熟妇中文字幕五十中出| 欧美zozozo另类| 国产麻豆成人av免费视频| 亚洲天堂国产精品一区在线| 亚洲自偷自拍图片 自拍| 国产乱人视频| 午夜免费观看网址| 搡老妇女老女人老熟妇| 精品熟女少妇八av免费久了| 日韩欧美国产在线观看| 久久久久国产精品人妻aⅴ院| 18禁观看日本| 男女床上黄色一级片免费看| 亚洲欧美日韩卡通动漫| 99热6这里只有精品| 国产又黄又爽又无遮挡在线| 一个人看视频在线观看www免费 | 99视频精品全部免费 在线 | 亚洲第一欧美日韩一区二区三区| 国产1区2区3区精品| 在线a可以看的网站| 18禁黄网站禁片免费观看直播| 久久久水蜜桃国产精品网| 亚洲精华国产精华精| e午夜精品久久久久久久| 啦啦啦观看免费观看视频高清| 亚洲专区国产一区二区| 日韩成人在线观看一区二区三区| 久久久久久久午夜电影| 桃红色精品国产亚洲av| 熟女少妇亚洲综合色aaa.| 黑人操中国人逼视频| 999久久久精品免费观看国产| 男人舔女人下体高潮全视频| 91av网站免费观看| 精品一区二区三区四区五区乱码| 亚洲人成网站高清观看| 小说图片视频综合网站| av欧美777| 免费看美女性在线毛片视频| 亚洲人成电影免费在线| 亚洲天堂国产精品一区在线| 国产精品国产高清国产av| av中文乱码字幕在线| 日韩三级视频一区二区三区| 给我免费播放毛片高清在线观看| 精品久久久久久成人av| 国产精品乱码一区二三区的特点| 嫩草影视91久久| 波多野结衣高清作品| 男插女下体视频免费在线播放| 亚洲国产精品久久男人天堂| 欧美绝顶高潮抽搐喷水| 国产精品久久久久久亚洲av鲁大| 岛国视频午夜一区免费看| 变态另类成人亚洲欧美熟女| www.自偷自拍.com| 美女高潮喷水抽搐中文字幕| 99久久无色码亚洲精品果冻| 国产97色在线日韩免费| 国产精品野战在线观看| 亚洲电影在线观看av| 国产高清激情床上av| 观看美女的网站| 女警被强在线播放| 欧美一级a爱片免费观看看| 欧美高清成人免费视频www| 精品国产超薄肉色丝袜足j| 在线视频色国产色| 变态另类丝袜制服| 国产成人影院久久av| 日本撒尿小便嘘嘘汇集6| 欧美一级毛片孕妇| 精品国产乱子伦一区二区三区| 婷婷亚洲欧美| 1024香蕉在线观看| 亚洲成人久久爱视频| 别揉我奶头~嗯~啊~动态视频| 999精品在线视频| 午夜免费成人在线视频| 国产精品影院久久| 国产欧美日韩一区二区三| 欧美性猛交╳xxx乱大交人| 天天躁日日操中文字幕| 成年免费大片在线观看| 国产精品1区2区在线观看.| 国产一区二区在线av高清观看| 狂野欧美白嫩少妇大欣赏| 国内精品一区二区在线观看| 国产午夜精品久久久久久| 色视频www国产| 中文亚洲av片在线观看爽| 国产精品自产拍在线观看55亚洲| 天堂网av新在线| 美女大奶头视频| 青草久久国产| 日韩精品中文字幕看吧| 国产亚洲欧美在线一区二区| 亚洲第一电影网av| 欧美又色又爽又黄视频| 长腿黑丝高跟| 亚洲五月天丁香| 在线观看日韩欧美| 亚洲狠狠婷婷综合久久图片| 啪啪无遮挡十八禁网站| 一级作爱视频免费观看| 亚洲精品乱码久久久v下载方式 | 无人区码免费观看不卡| 亚洲男人的天堂狠狠| 欧美日韩瑟瑟在线播放| 五月伊人婷婷丁香| 我要搜黄色片| 麻豆一二三区av精品| 国产又色又爽无遮挡免费看| 黄色片一级片一级黄色片| 国产黄片美女视频| 国产精品精品国产色婷婷| 女人被狂操c到高潮| 国内精品美女久久久久久| 91在线精品国自产拍蜜月 | 亚洲国产高清在线一区二区三| 亚洲九九香蕉| 窝窝影院91人妻| 久久亚洲精品不卡| 亚洲人与动物交配视频| 久久精品91蜜桃| 国产高清有码在线观看视频| 国产熟女xx| 欧美午夜高清在线| 99久久综合精品五月天人人| 男女午夜视频在线观看| 两性夫妻黄色片| 一个人看视频在线观看www免费 | 成年免费大片在线观看| 免费在线观看影片大全网站| 高清毛片免费观看视频网站| 国产精品98久久久久久宅男小说| 亚洲av第一区精品v没综合| 99精品在免费线老司机午夜| 国产私拍福利视频在线观看| 欧洲精品卡2卡3卡4卡5卡区| 啦啦啦免费观看视频1| 久久久国产成人免费| 老司机在亚洲福利影院| 精品国产乱码久久久久久男人| 中文字幕久久专区| 亚洲激情在线av| 99久国产av精品| 91字幕亚洲| 母亲3免费完整高清在线观看| 日日摸夜夜添夜夜添小说| 黄色女人牲交| 欧美中文综合在线视频| 在线视频色国产色| e午夜精品久久久久久久| 老汉色∧v一级毛片| av视频在线观看入口| 看片在线看免费视频| 国产精品久久久久久精品电影| x7x7x7水蜜桃| 1000部很黄的大片| 久久性视频一级片| 欧美高清成人免费视频www| 97人妻精品一区二区三区麻豆| 最近视频中文字幕2019在线8| 成人av一区二区三区在线看| 国产精品爽爽va在线观看网站| 一级作爱视频免费观看| 少妇的逼水好多| 亚洲 国产 在线| 亚洲性夜色夜夜综合| 久9热在线精品视频| 国内精品久久久久精免费| 亚洲电影在线观看av| 亚洲av成人av| 久久九九热精品免费| 日本在线视频免费播放| 日本成人三级电影网站| 全区人妻精品视频| 一二三四在线观看免费中文在| 久久欧美精品欧美久久欧美| 久久精品国产亚洲av香蕉五月| 国产精品自产拍在线观看55亚洲| 亚洲欧美激情综合另类| 免费无遮挡裸体视频| 美女大奶头视频| 成人特级av手机在线观看| 国产真实乱freesex| 成熟少妇高潮喷水视频| 黄色丝袜av网址大全| 一二三四在线观看免费中文在| 免费电影在线观看免费观看| 久久午夜亚洲精品久久| 看黄色毛片网站| av黄色大香蕉| 精品国产乱子伦一区二区三区| 午夜两性在线视频| 丁香欧美五月| а√天堂www在线а√下载| 中出人妻视频一区二区| 少妇熟女aⅴ在线视频| 午夜影院日韩av| 搡老妇女老女人老熟妇| 成熟少妇高潮喷水视频| 国产精品99久久99久久久不卡| av天堂在线播放| 叶爱在线成人免费视频播放| 99久久精品国产亚洲精品| 在线十欧美十亚洲十日本专区| 午夜a级毛片| 丝袜人妻中文字幕| 真人做人爱边吃奶动态| 国产高潮美女av| 视频区欧美日本亚洲| 床上黄色一级片| 久久久久久久久免费视频了| 精品乱码久久久久久99久播| 国产高清三级在线| 国产精品免费一区二区三区在线| 国产高清videossex| 18禁黄网站禁片免费观看直播| av片东京热男人的天堂| 亚洲成人精品中文字幕电影| 免费看日本二区| 精品人妻1区二区| 亚洲人成电影免费在线| 午夜福利高清视频| 国产精品美女特级片免费视频播放器 | 少妇丰满av| 首页视频小说图片口味搜索| 亚洲乱码一区二区免费版| 欧美色欧美亚洲另类二区| 男人舔奶头视频| 精品熟女少妇八av免费久了| 国产乱人视频| aaaaa片日本免费| 又紧又爽又黄一区二区| 一级a爱片免费观看的视频| 亚洲天堂国产精品一区在线| 变态另类成人亚洲欧美熟女| 亚洲无线在线观看| 法律面前人人平等表现在哪些方面| 夜夜爽天天搞| 国产高清videossex| 两个人的视频大全免费| av视频在线观看入口| 看免费av毛片| 国产精品久久久久久精品电影| 日本一二三区视频观看| 窝窝影院91人妻| 亚洲 国产 在线| 国产高清有码在线观看视频| 一级毛片精品| 亚洲av成人不卡在线观看播放网| 国产黄色小视频在线观看| 两人在一起打扑克的视频| 国产视频一区二区在线看| 午夜亚洲福利在线播放| 欧美大码av| 亚洲欧美日韩无卡精品| 成人三级做爰电影| 久久久成人免费电影| 国产成人啪精品午夜网站| 国产伦精品一区二区三区视频9 | 一卡2卡三卡四卡精品乱码亚洲| 成人三级黄色视频| 欧美中文综合在线视频| 国产精品一及| 亚洲在线自拍视频| 久久精品综合一区二区三区| 好看av亚洲va欧美ⅴa在| 舔av片在线| 欧美高清成人免费视频www| 最近视频中文字幕2019在线8| 好看av亚洲va欧美ⅴa在| 免费观看人在逋| 国产精品99久久久久久久久| 亚洲黑人精品在线| 亚洲中文字幕日韩| 亚洲欧美日韩高清专用| 国产精品久久久久久亚洲av鲁大| 国产单亲对白刺激| 男女做爰动态图高潮gif福利片| 成人18禁在线播放| 五月玫瑰六月丁香| 成人欧美大片| 嫩草影院精品99| 国产精品99久久久久久久久| 一个人观看的视频www高清免费观看 | 亚洲欧美日韩无卡精品| 午夜成年电影在线免费观看| 久久婷婷人人爽人人干人人爱| 欧美国产日韩亚洲一区| 中文字幕人妻丝袜一区二区| 99久久精品一区二区三区| 国产精品九九99| 亚洲国产欧美网| 听说在线观看完整版免费高清| 国产成+人综合+亚洲专区| 欧美日韩乱码在线| 亚洲av五月六月丁香网| a在线观看视频网站| 夜夜看夜夜爽夜夜摸| 免费观看的影片在线观看| 亚洲中文字幕一区二区三区有码在线看 | 国内精品一区二区在线观看| 又大又爽又粗| 欧美日韩黄片免| 日韩有码中文字幕| 日韩欧美精品v在线| 精品久久蜜臀av无| 九九久久精品国产亚洲av麻豆 | 男女下面进入的视频免费午夜| 丰满人妻熟妇乱又伦精品不卡| 丰满人妻一区二区三区视频av | 国产精品99久久久久久久久| 嫩草影院入口| 两个人的视频大全免费| 桃色一区二区三区在线观看| 国产又色又爽无遮挡免费看| 黑人欧美特级aaaaaa片| 琪琪午夜伦伦电影理论片6080| 欧美乱码精品一区二区三区| 99久久无色码亚洲精品果冻| 亚洲中文字幕一区二区三区有码在线看 | 国产精品99久久久久久久久| 特级一级黄色大片| 黄色女人牲交| 女人被狂操c到高潮| 国产极品精品免费视频能看的| 欧美性猛交╳xxx乱大交人| 精品国产乱子伦一区二区三区| 男人舔女人下体高潮全视频| 国产亚洲欧美98| 一区二区三区高清视频在线| 18美女黄网站色大片免费观看| 一个人看视频在线观看www免费 | 天堂动漫精品| 久久久国产成人免费| 亚洲电影在线观看av| 国产亚洲精品久久久久久毛片| 欧美av亚洲av综合av国产av| 亚洲五月婷婷丁香| 白带黄色成豆腐渣| 亚洲激情在线av| 女生性感内裤真人,穿戴方法视频| 狂野欧美白嫩少妇大欣赏| 国产精品久久久久久久电影 | 99久久久亚洲精品蜜臀av| 午夜福利在线观看吧| 网址你懂的国产日韩在线| 亚洲无线在线观看| 一级毛片女人18水好多| 99久久精品热视频| 99久久综合精品五月天人人| 男女午夜视频在线观看| 在线观看免费午夜福利视频| av在线蜜桃| 高清毛片免费观看视频网站| 黄色成人免费大全| 99国产精品一区二区蜜桃av| 国产熟女xx| 成人特级av手机在线观看| 久久久精品大字幕| 国产精品久久久久久精品电影| www.精华液| 午夜福利欧美成人| 成人亚洲精品av一区二区| 国产三级在线视频| 日韩有码中文字幕| 一区二区三区激情视频| 亚洲黑人精品在线| 手机成人av网站| 久久久久久久精品吃奶| 国产精品99久久久久久久久| 极品教师在线免费播放| 日日夜夜操网爽| 亚洲精品美女久久久久99蜜臀| 亚洲avbb在线观看| 午夜a级毛片| 国产精品 欧美亚洲| 亚洲美女视频黄频| 成人三级做爰电影| 国产高清videossex| 久久精品aⅴ一区二区三区四区| 国产精品久久电影中文字幕| 亚洲自拍偷在线| 久久精品91无色码中文字幕| 亚洲精品一卡2卡三卡4卡5卡| 国产成人精品无人区| 小蜜桃在线观看免费完整版高清| 色综合站精品国产| 淫妇啪啪啪对白视频| 国产成年人精品一区二区| 国产又黄又爽又无遮挡在线| 丰满人妻一区二区三区视频av | 青草久久国产| 亚洲精品一卡2卡三卡4卡5卡| 亚洲狠狠婷婷综合久久图片| 两个人看的免费小视频| 亚洲美女视频黄频| 很黄的视频免费| 亚洲电影在线观看av| 欧美xxxx黑人xx丫x性爽| 精品福利观看| 搡老妇女老女人老熟妇| 国产免费男女视频| 噜噜噜噜噜久久久久久91| 免费一级毛片在线播放高清视频| 国产一区二区在线av高清观看| 久久国产精品人妻蜜桃| 亚洲av成人不卡在线观看播放网| av片东京热男人的天堂| 国产高清视频在线播放一区| 日本 欧美在线| 国产亚洲精品久久久久久毛片| 免费观看的影片在线观看| 999精品在线视频| 别揉我奶头~嗯~啊~动态视频| 91字幕亚洲| 国产精品美女特级片免费视频播放器 | 丁香欧美五月| 欧美一级毛片孕妇| 在线视频色国产色| 成人亚洲精品av一区二区| 久久婷婷人人爽人人干人人爱| 欧美大码av| 国产精品一区二区三区四区久久| 高潮久久久久久久久久久不卡| 麻豆国产97在线/欧美| 一级毛片高清免费大全| 淫秽高清视频在线观看| 亚洲中文日韩欧美视频| 精品国产乱码久久久久久男人| 欧美日韩瑟瑟在线播放| 夜夜爽天天搞| h日本视频在线播放| a级毛片在线看网站| 一级黄色大片毛片| 亚洲五月婷婷丁香| 午夜福利在线观看免费完整高清在 | 美女cb高潮喷水在线观看 | 精品欧美国产一区二区三| 亚洲一区二区三区色噜噜| a在线观看视频网站| 动漫黄色视频在线观看| 国模一区二区三区四区视频 | 母亲3免费完整高清在线观看| 小说图片视频综合网站| 给我免费播放毛片高清在线观看| 国产爱豆传媒在线观看| aaaaa片日本免费| 亚洲国产精品sss在线观看| 欧美中文综合在线视频| 国产私拍福利视频在线观看| 国内精品久久久久精免费| 国产精品爽爽va在线观看网站| 亚洲自偷自拍图片 自拍| 欧美成人性av电影在线观看| 日本成人三级电影网站| 精品一区二区三区视频在线 | 国产又色又爽无遮挡免费看| 久久久久精品国产欧美久久久| 十八禁网站免费在线| 免费在线观看日本一区| or卡值多少钱| 亚洲午夜精品一区,二区,三区| 久久久久九九精品影院| 三级国产精品欧美在线观看 | 国产精品综合久久久久久久免费| 国产淫片久久久久久久久 | 中亚洲国语对白在线视频| 欧美日韩黄片免| 久久这里只有精品19| 少妇人妻一区二区三区视频| 国产精品 欧美亚洲| 亚洲精品色激情综合| 免费一级毛片在线播放高清视频| 久久中文字幕一级| 国产精品一区二区精品视频观看| 香蕉国产在线看| 此物有八面人人有两片| 欧美成人免费av一区二区三区| 久久久久久久久久黄片| 亚洲午夜理论影院| 色播亚洲综合网| 国内精品美女久久久久久| 在线看三级毛片| 黄色女人牲交| 亚洲av中文字字幕乱码综合| 狠狠狠狠99中文字幕|