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

    基于RNN的鋰離子電池SOC估算研究進(jìn)展

    2021-03-07 02:29:42劉雨辰
    電源技術(shù) 2021年2期
    關(guān)鍵詞:優(yōu)化模型

    徐 帥,劉雨辰, 周 飛

    (南京航空航天大學(xué)直升機(jī)傳動(dòng)技術(shù)國(guó)家級(jí)重點(diǎn)實(shí)驗(yàn)室,江蘇南京210016)

    電池荷電狀態(tài)(SOC)是電池管理系統(tǒng)中最重要的部分[1]。精確的SOC 估算可以及時(shí)更新電池的剩余容量,防止發(fā)生過充和過放電現(xiàn)象,并提高電池的充放電效率和延長(zhǎng)其使用壽命。雖然目前SOC 估算方法眾多,但是主流的方法主要包含以下四類:安時(shí)積分法[2]、開路電壓法[3]、模型法[4-5]、機(jī)器學(xué)習(xí)法[6]。深度學(xué)習(xí)算法屬于機(jī)器學(xué)習(xí)法??裳h(huán)神經(jīng)網(wǎng)絡(luò)(RNN)作為深度學(xué)習(xí)算法中的一種,是目前最流行的算法之一。RNN 是一類以序列數(shù)據(jù)為輸入,在序列的演進(jìn)方向進(jìn)行遞歸且所有循環(huán)節(jié)點(diǎn)按鏈?zhǔn)竭B接的遞歸神經(jīng)網(wǎng)絡(luò),但其在訓(xùn)練時(shí),會(huì)出現(xiàn)“梯度消失”和“梯度爆炸”的現(xiàn)象。為了解決這兩個(gè)問題,學(xué)者們開發(fā)了一些進(jìn)階版循環(huán)神經(jīng)網(wǎng)絡(luò),比如長(zhǎng)短期記憶神經(jīng)網(wǎng)絡(luò)(LSTM),門控循環(huán)單元神經(jīng)網(wǎng)絡(luò)(GRU)以及雙向長(zhǎng)短期記憶神經(jīng)網(wǎng)絡(luò)(BiLSTM)。這些由各種門結(jié)構(gòu)組成的人工神經(jīng)網(wǎng)絡(luò),擁有回憶過去信息的能力,這使得它們特別適用于解決較長(zhǎng)時(shí)間序列的問題,例如電池SOC 估算。本文將估算電池SOC 的可循環(huán)神經(jīng)網(wǎng)絡(luò)算法根據(jù)其結(jié)構(gòu)的不同分為三種:簡(jiǎn)單循環(huán)神經(jīng)網(wǎng)絡(luò)算法、深度循環(huán)神經(jīng)網(wǎng)絡(luò)算法、擴(kuò)展循環(huán)神經(jīng)網(wǎng)絡(luò)算法。

    1 簡(jiǎn)單循環(huán)神經(jīng)網(wǎng)絡(luò)算法

    1.1 門結(jié)構(gòu)循環(huán)神經(jīng)網(wǎng)絡(luò)

    1.1.1 長(zhǎng)短期記憶神經(jīng)網(wǎng)絡(luò)

    長(zhǎng)短期記憶神經(jīng)網(wǎng)絡(luò)(LSTM)層結(jié)構(gòu)如圖1 所示。存儲(chǔ)記憶單元Ck來存儲(chǔ)并傳輸先前和當(dāng)前狀態(tài)的信息,以備將來使用;遺忘門fk負(fù)責(zé)控制Ck-1中應(yīng)該傳遞多少先前信息到下一次;輸入門ik負(fù)責(zé)調(diào)節(jié)應(yīng)該將多少過去的信息傳遞到存儲(chǔ)記憶單元Ck中;輸出門ok負(fù)責(zé)控制將多少在存儲(chǔ)記憶單元中計(jì)算出的信息傳遞到輸出hk中。

    圖1 長(zhǎng)短期記憶神經(jīng)網(wǎng)絡(luò)(LSTM)層結(jié)構(gòu)

    Bockrath 等[7]使用LSTM 建立電池測(cè)量信號(hào)(端電壓,負(fù)載電流和電池溫度)與電池SOC 之間非線性關(guān)系,并與擴(kuò)展卡爾曼濾波算法(EKF)進(jìn)行了對(duì)比。相較于EKF,LSTM 不需要建立復(fù)雜的電池模型,而且LSTM 的平均絕對(duì)誤差(MAE)和均方根誤差(RMSE)分別為4.1%和5.0%,也遠(yuǎn)低于EKF 的8.5%和9.5%。然而,這篇文章并沒有采用標(biāo)準(zhǔn)的汽車工況數(shù)據(jù),而是普通的充放電循環(huán)數(shù)據(jù)(CC),環(huán)境溫度也只考慮了常溫。在這種條件下對(duì)比LSTM 和EKF,顯然沒有很大的說服力。Li 等[8]同樣應(yīng)用了LSTM 模型,不過他使用多種汽車工況循環(huán)數(shù)據(jù)集(UDDS,US06,HWFET,LA92)作為測(cè)試文件,而且還考慮了多種環(huán)境溫度(0、10 和25 ℃)。同時(shí)為了進(jìn)一步提高估算精度,作者還對(duì)測(cè)試數(shù)據(jù)進(jìn)行了歸一化處理,將數(shù)據(jù)縮放在-1 到1 之間。不同溫度下的MAE 均小于1.3%的估算結(jié)果證明了LSTM 優(yōu)秀的魯棒性。后來Zhang 等[9]考慮到電池老化對(duì)SOC 估算的影響,增加了電池剩余容量和電池內(nèi)阻作為L(zhǎng)STM 的輸入?yún)?shù),使用不同電池健康狀態(tài)(SOH)條件下的美國(guó)動(dòng)態(tài)應(yīng)力測(cè)試(DST)和中國(guó)電動(dòng)汽車測(cè)試標(biāo)準(zhǔn)(QCT)工況循環(huán)數(shù)據(jù)來測(cè)試LSTM。SOC 估算平均誤差均小于2%的估算結(jié)果證明了在不同的電池SOH 下,LSTM 也適用于估算SOC。以上文獻(xiàn)中的LSTM 模型都是針對(duì)某一種類型的電池,然而現(xiàn)實(shí)生活中,電池的類型肯定不止一種。為了讓LSTM 可以適用于不同類型的電池,Vidal 等[10]提出了基于遷移學(xué)習(xí)的LSTM 模型,首先使用其他類型電池的數(shù)據(jù)對(duì)LSTM 神經(jīng)網(wǎng)絡(luò)進(jìn)行訓(xùn)練,然后使用來自目標(biāo)電池類型的數(shù)據(jù)對(duì)神經(jīng)網(wǎng)絡(luò)進(jìn)行二次訓(xùn)練,最后應(yīng)用幾種不同的實(shí)驗(yàn)來研究其潛力。從實(shí)驗(yàn)估算結(jié)果來看,使用遷移學(xué)習(xí)方法后,LSTM 估算SOC 的準(zhǔn)確度提高了64%,訓(xùn)練時(shí)間從9 h 縮減到了1 h,這證明了遷移學(xué)習(xí)有提高SOC 估算精度并減少神經(jīng)網(wǎng)絡(luò)訓(xùn)練所需數(shù)據(jù)的潛力。這為實(shí)際中應(yīng)用LSTM 去估算不同類型的電池提供了可能。雖然以上的作者們證明了LSTM在各種條件下都可以很好地估算電池SOC,但是估算結(jié)果的好壞主要是由他們所搭建LSTM 模型的超參數(shù)決定的,這些超參數(shù)包括但不局限于模型隱藏層數(shù),模型隱藏層神經(jīng)元數(shù),模型優(yōu)化器,訓(xùn)練批量數(shù),模型輸入時(shí)間步長(zhǎng)。遺憾的是,他們只是介紹了他們選擇的超參數(shù),并沒有過多地解釋這么選的原因。

    Chemali 等[11]提出了一個(gè)隱藏層神經(jīng)元數(shù)為500 的單層LSTM 模型來估算電池SOC,并討論了三種訓(xùn)練批量數(shù)對(duì)訓(xùn)練結(jié)果的影響。相同條件下,當(dāng)訓(xùn)練批量數(shù)從250 增加到500時(shí),平均誤差下降了50%,而繼續(xù)增加到1 000 時(shí),平均誤差只下降了15%。但是在SOC 初始值存在30%的誤差時(shí),訓(xùn)練批量數(shù)為1 000 的神經(jīng)網(wǎng)絡(luò)的收斂速率比500 的神經(jīng)網(wǎng)絡(luò)快了近50 s。雖然這篇文章只考慮了一種超參數(shù)對(duì)LSTM 網(wǎng)絡(luò)估算結(jié)果的影響,但是也從側(cè)面說明了,在用LSTM 對(duì)電池SOC進(jìn)行估算時(shí),超參數(shù)的選擇是個(gè)非常重要的問題。Hong 等[12]綜合考慮了四種超參數(shù),討論了不同的輸入數(shù)據(jù)時(shí)間步長(zhǎng),輸出數(shù)據(jù)時(shí)間步長(zhǎng),和訓(xùn)練樣本總數(shù)以及訓(xùn)練批量數(shù)對(duì)估算結(jié)果的影響。首先固定訓(xùn)練樣本總數(shù)以及訓(xùn)練批量數(shù),對(duì)比了不同的輸入數(shù)據(jù)時(shí)間步長(zhǎng),輸出數(shù)據(jù)時(shí)間步長(zhǎng)的估算誤差。從表1 誤差結(jié)果可以發(fā)現(xiàn),隨著輸入數(shù)據(jù)時(shí)間步長(zhǎng)和輸出數(shù)據(jù)時(shí)間步長(zhǎng)的增大,均方根誤差也隨之增大;然后選擇最優(yōu)的輸入和輸出數(shù)據(jù)時(shí)間步長(zhǎng),分析不同訓(xùn)練樣本總數(shù)對(duì)估算結(jié)果的影響,對(duì)于LSTM 來說,輸入數(shù)據(jù)當(dāng)然是越多越好;最后分析了不同訓(xùn)練批量數(shù)的影響。從圖2 來看,估算誤差并沒有隨著訓(xùn)練批量數(shù)的增加而減少,當(dāng)其值為16 時(shí),誤差最小??上У氖乾F(xiàn)在沒有成熟的算法來選擇循環(huán)神經(jīng)網(wǎng)絡(luò)模型超參數(shù),而它們又決定了模型的估算精度,研究者往往都是根據(jù)經(jīng)驗(yàn)或通過枚舉對(duì)比來確定它們,這不僅耗時(shí),而且還不一定能找到最優(yōu)值。

    表1 不同時(shí)間步長(zhǎng)的均方根誤差[12]

    圖2 基于不同訓(xùn)練批量數(shù)的均方根誤差[12]

    1.1.2 門控循環(huán)單元神經(jīng)網(wǎng)絡(luò)

    圖3 GRU 層結(jié)構(gòu)

    另一種能夠處理長(zhǎng)期依賴關(guān)系的RNN 方法是門控循環(huán)單元神經(jīng)網(wǎng)絡(luò)(GRU),它也使用門結(jié)構(gòu)來學(xué)習(xí),記憶和確定過去和現(xiàn)在的哪些信息將用于生成其輸出。它的層結(jié)構(gòu)如圖3所示,圖中的zk和rk分別表示更新門和重置門。更新門用于控制前一時(shí)刻的狀態(tài)信息hk-1被帶入到當(dāng)前狀態(tài)hk中的程度,更新門的值越大說明前一時(shí)刻的狀態(tài)信息被帶入的越多,重置門控制前一狀態(tài)有多少信息被寫入到當(dāng)前的候選集h?k上,重置門越小,前一狀態(tài)的信息被寫入的越少。與LSTM 相比,GRU 使用單個(gè)門單元同時(shí)控制遺忘量和更新狀態(tài)單元,所以GRU 的參數(shù)相對(duì)要少很多,訓(xùn)練時(shí)間相對(duì)較短,但可以達(dá)到與LSTM 相當(dāng)?shù)男Ч?/p>

    Li 等[13]首次提出使用GRU 來估算電池SOC,并使用了和文獻(xiàn)[11]中相同的電池?cái)?shù)據(jù)集,而且在25 ℃的LA92 工況下的估算誤差MAE 為0.32%。由于本文中使用的GRU 的超參數(shù)與文獻(xiàn)[11]中LSTM 的超參數(shù)不同,所以無法直接對(duì)比兩種模型。不過這篇論文還是證明了GRU 是可以像LSTM 一樣用于SOC 估算的。Yang 等[14]提出了GRU 模型來估算不同環(huán)境溫度下NMC 和LFP 兩種電池SOC。值得一提的是,他在模型中添加了一個(gè)參數(shù)為0.5 的dropout 層,可以在訓(xùn)練的時(shí)候隨機(jī)丟棄50%的信息,從而避免過擬合現(xiàn)象。不過其缺點(diǎn)就是丟棄的信息也有可能是非常重要的,從而導(dǎo)致估算精度不高,可惜的是,文中沒有更加深入地分析dropout 對(duì)估算結(jié)果的影響。和文獻(xiàn)[12]一樣,文中也分析了多種模型超參數(shù)對(duì)估算結(jié)果的影響,比如訓(xùn)練迭代次數(shù),隱藏層數(shù)和隱藏層神經(jīng)元數(shù),并在訓(xùn)練中增加了驗(yàn)證集,使得對(duì)比結(jié)果更有說服力。其分析策略是先給超參數(shù)設(shè)定一個(gè)常見的初始值,然后固定其中兩個(gè)參數(shù),分析另一個(gè)參數(shù)的變化給估算結(jié)果帶來的影響。從誤差結(jié)果來看,在訓(xùn)練前期,隨著迭代次數(shù)的增加,訓(xùn)練集和驗(yàn)證集誤差RMSE 都開始急劇下降,當(dāng)?shù)螖?shù)超過1 500 后,誤差開始逐漸穩(wěn)定,并在2 000 時(shí)達(dá)到最小值。但是隨著迭代次數(shù)進(jìn)一步增加,驗(yàn)證集誤差開始出現(xiàn)上升跡象,表明訓(xùn)練過程開始出現(xiàn)過擬合。綜合考慮測(cè)試性能和與訓(xùn)練迭代次數(shù)呈線性關(guān)系的訓(xùn)練時(shí)間后,作者選擇訓(xùn)練迭代次數(shù)為2 000。至于隱藏層數(shù)和隱藏層神經(jīng)元數(shù),其對(duì)估算結(jié)果的影響也是:當(dāng)他們一定大的時(shí)候,訓(xùn)練集和測(cè)試集誤差會(huì)非常小,但是進(jìn)一步增大都會(huì)導(dǎo)致過擬合現(xiàn)象的發(fā)生,從而增大了驗(yàn)證集誤差,導(dǎo)致模型的泛化能力變差。

    為了提高GRU 估算電池SOC 的精度,Xiao 等[15]使用均方根誤差(MSE)作為GRU 模型的訓(xùn)練損失誤差,并使用集成優(yōu)化器來優(yōu)化它,即在訓(xùn)練迭代次數(shù)小于閾值p 時(shí),用優(yōu)化器Nadam 來優(yōu)化模型,因?yàn)樗哂惺諗克俣瓤斓膬?yōu)點(diǎn);當(dāng)訓(xùn)練迭代次數(shù)超過p 后,改換優(yōu)化器AdaMax 對(duì)模型進(jìn)行微調(diào),從而確定模型最終的訓(xùn)練權(quán)重和偏值。為了進(jìn)一步突出GRU 模型在集成優(yōu)化器下估算SOC 時(shí)的性能,文中還添加了LSTM模型作對(duì)比。測(cè)試數(shù)據(jù)使用的都是LFP 電池在DST,F(xiàn)UDS,US06 三種工況下的循環(huán)數(shù)據(jù)。從表2 誤差結(jié)果來看,同種模型下,相比單個(gè)優(yōu)化器,集成優(yōu)化器可以顯著地提高估算精度,而所花費(fèi)的時(shí)間只是略微增加。另一方面,在同種優(yōu)化器作用下,GRU 無論是在估計(jì)誤差,還是訓(xùn)練時(shí)間,都比LSTM 的小,其中GRU 在集成優(yōu)化器作用下的RMSE,MAE 和訓(xùn)練時(shí)間分別是1.13%、0.84%和2.97 h,與LSTM 在單個(gè)優(yōu)化器下的相應(yīng)值相比分別降低了46.7%、46.84%和17.96%。不過,由于不知道作者使用的LSTM 模型超參數(shù)是否與GRU 一樣,所以這個(gè)對(duì)比結(jié)果僅供參考,并不能證明GRU 一定比LSTM 好。

    表2 基于不同優(yōu)化器下GRU 和LSTM 的估算誤差[15]

    1.2 非門結(jié)構(gòu)循環(huán)神經(jīng)網(wǎng)絡(luò)

    除了門結(jié)構(gòu)的循環(huán)神經(jīng)網(wǎng)絡(luò)外,也有學(xué)者提出使用非門結(jié)構(gòu)的循環(huán)神經(jīng)網(wǎng)絡(luò)來估算電池SOC。非線性自回歸與外來輸入(NARX)神經(jīng)網(wǎng)絡(luò)是RNN 的子類,適用于預(yù)測(cè)非線性和時(shí)間序列問題。NARX 神經(jīng)網(wǎng)絡(luò)使用有限的反饋形成輸出層而不是隱藏層,并將模型k-1 時(shí)刻的輸出轉(zhuǎn)化成k 時(shí)刻的輸入,從而實(shí)現(xiàn)循環(huán)。在學(xué)習(xí)能力,收斂速度,泛化性能和高精度方面,NARXNN 比傳統(tǒng)的RNN 表現(xiàn)更好。Chaoui 等[16]開發(fā)了基于非線性自回歸與外來輸入神經(jīng)網(wǎng)絡(luò)架構(gòu)的動(dòng)態(tài)驅(qū)動(dòng)遞歸網(wǎng)絡(luò)(DDRN)來估計(jì)兩個(gè)鋰離子化學(xué)電池的SOC。雖然與門結(jié)構(gòu)循環(huán)神經(jīng)網(wǎng)絡(luò)相比,使用NARX 會(huì)給DDRN 帶來一些限制,但是這也賦予了DDRN 關(guān)聯(lián)儲(chǔ)存功能,并且減少了模型訓(xùn)練所需的數(shù)據(jù)量。其結(jié)構(gòu)如圖4 所示,文中使用0、10、25 和40 ℃溫度下包含電壓Vk,電流Ik,環(huán)境溫度Tk和來自先前的時(shí)間步的輸出SOCk-1四種變量的數(shù)據(jù)集來訓(xùn)練和測(cè)試DDRN。從估算結(jié)果來看,與“非循環(huán)”神經(jīng)網(wǎng)絡(luò)相比,DRNN的訓(xùn)練時(shí)間減少了1 000 倍,均方根誤差下降了78%。盡管用于測(cè)試模型的數(shù)據(jù)循環(huán)曲線是動(dòng)態(tài)的,但由于在循環(huán)過程中沒有再生制動(dòng)充電脈沖,因此導(dǎo)致了單調(diào)的SOC 輸出曲線。此類測(cè)試文件可能無法很好地應(yīng)用在電動(dòng)汽車領(lǐng)域。因此,這種方法的有效性需要進(jìn)一步研究。

    圖4 DDRN-NARX 結(jié)構(gòu),輸入變量端電壓Vk,電流Ik和環(huán)境溫度Tk都是由傳感器測(cè)量得到,回歸輸入SOCk-1是上一時(shí)刻k-1的輸出[16]

    前面已經(jīng)提到,門結(jié)構(gòu)循環(huán)神經(jīng)網(wǎng)絡(luò)估算SOC 的精度由模型的超參數(shù)所決定,對(duì)于非門結(jié)構(gòu)循環(huán)神經(jīng)網(wǎng)絡(luò)來說,亦是如此。而這些作者都是根據(jù)經(jīng)驗(yàn)或者通過枚舉對(duì)比的方法來確定這些超參數(shù),步驟十分繁瑣且低效。因此Lipu 等[17]提出了基于改進(jìn)的非線性具有基于外源輸入的自回歸神經(jīng)網(wǎng)絡(luò)算法來估算電池SOC,并結(jié)合了線性搜索算法(LSA)來查找模型最優(yōu)超參數(shù):神經(jīng)元數(shù)量,輸入數(shù)據(jù)時(shí)間步長(zhǎng)和輸出數(shù)據(jù)時(shí)間步長(zhǎng)的最佳組合,從而提高模型的估算精度。該模型使用了NCM 電池在0、10 和45 ℃三種溫度下基于FUDS 和US06 兩種工況的數(shù)據(jù)集作為訓(xùn)練集和測(cè)試集,并獲得了均小于0.9%的MAE 和均小于1.3%的RMSE。文中還添加了與同樣由LSA 優(yōu)化的其他SOC 估計(jì)方法的比較,并證明了NARX和LSA 結(jié)合可以獲得更高的精度。盡管使用LSA 和BSA 之類的優(yōu)化算法可以幫助實(shí)現(xiàn)構(gòu)建最佳神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)的自動(dòng)化,但是其他重要的超參數(shù),例如訓(xùn)練批量大小,損失函數(shù)優(yōu)化器的選擇同樣會(huì)對(duì)估算結(jié)果產(chǎn)生重大影響,而且作者沒有給出使用LSA 算法來查找模型超參數(shù)的時(shí)間,如果這種算法需要耗費(fèi)大量的時(shí)間,那么顯然不適合實(shí)際應(yīng)用。

    2 深度循環(huán)神經(jīng)網(wǎng)絡(luò)算法

    2.1 堆疊循環(huán)神經(jīng)網(wǎng)絡(luò)

    與簡(jiǎn)單循環(huán)神經(jīng)網(wǎng)絡(luò)法相比,深度循環(huán)神經(jīng)網(wǎng)絡(luò)法使用多層LSTM,GRU 單元。循環(huán)神經(jīng)網(wǎng)絡(luò)能夠以多種方式由單層加深至多層,其中最常見的策略是使用堆疊的循環(huán)單元。由于在序列演進(jìn)方向已經(jīng)存在復(fù)雜結(jié)構(gòu),因此不同于深度的前饋神經(jīng)網(wǎng)絡(luò),深度循環(huán)神經(jīng)網(wǎng)絡(luò)在輸入和輸出間不會(huì)堆疊太多層次,一個(gè)3 層的深度循環(huán)神經(jīng)網(wǎng)絡(luò)已經(jīng)具有很大規(guī)模。堆疊循環(huán)神經(jīng)網(wǎng)絡(luò)是在全連接的單層循環(huán)神經(jīng)網(wǎng)絡(luò)的基礎(chǔ)上堆疊形成的深度算法。其內(nèi)部循環(huán)單元的狀態(tài)更新使用了其前一層相同時(shí)間步的狀態(tài)和當(dāng)前層前一時(shí)間步的狀態(tài)。

    Khalid 等[18]使用一個(gè)含兩層LSTM 的神經(jīng)網(wǎng)絡(luò)(SLSTM)來估算電池SOC 并分析了兩種優(yōu)化器,RMSprop 和AdaGrad對(duì)估算結(jié)果的影響。從表3 估算誤差來看,在不同的工況下,無論是訓(xùn)練還是測(cè)試,Adagrad 的誤差均小于RMSprop。RMSprop 的RMSE 和MAE 分 別 是2.611% 和1.967%,而AdaGrad 的RMSE 和MAE 分別是1.527%和1.153%,比前者降低了約50%。LSTM 的優(yōu)化器有很多種,比如NAdam,Adadelta,AdaGrad,RMSprop,Adamax,SGD。雖然本文只比較了其中兩種,不過對(duì)比結(jié)果也說明了,在用LSTM 估算電池SOC 時(shí),選擇合適的優(yōu)化器,可以得到更好的估算結(jié)果。

    表3 不同優(yōu)化器下的估算誤差[18]

    為了提高循環(huán)神經(jīng)網(wǎng)絡(luò)估算SOC 的精度,Yang 等[19]使用了一個(gè)堆疊了三層LSTM 的模型來估算SOC,其中每個(gè)LSTM 含有50 個(gè)神經(jīng)元。堆疊的LSTM 層可以重組從先前層學(xué)到的表示形式,并在高度抽象的層次上創(chuàng)建新的表示形式。雖然該模型在三種工況FUDS,DST,US06 下,均獲得了MAE 和RMSE 分別小于2.1%和1.4%的良好估算結(jié)果,但是由于沒有添加單層LSTM 模型做對(duì)比,而且本文只考慮一種環(huán)境溫度,所以無法證明多層LSTM 比單層LSTM 更好。

    2.2 雙向循環(huán)神經(jīng)網(wǎng)絡(luò)

    圖5 雙向循環(huán)神經(jīng)網(wǎng)絡(luò)展開圖

    雙向循環(huán)神經(jīng)網(wǎng)絡(luò)是一種兩層的深度循環(huán)神經(jīng)網(wǎng)絡(luò)。如圖5 所示,其結(jié)構(gòu)由兩個(gè)分別向前和向后的單向循環(huán)層(LSTM,GRU)組成,其中,向前的循環(huán)層使用k-n 時(shí)刻到當(dāng)前k 時(shí)刻的輸入數(shù)據(jù)Ψ,即[Ψk-n→Ψk],而向后的循環(huán)層則使用相反順序的輸入數(shù)據(jù)[Ψk→Ψk-n],其中n 是時(shí)間步長(zhǎng)的總數(shù),y 是雙向循環(huán)層的輸出,而且這兩個(gè)循環(huán)層是同時(shí)進(jìn)行反饋的,并且每個(gè)參數(shù)都是獨(dú)立更新的,但是它們的輸出hk,會(huì)用一個(gè)函數(shù)(求和、求積、求商、求平均)結(jié)合起來。由于這種雙向結(jié)構(gòu)從數(shù)據(jù)的兩個(gè)端捕獲時(shí)間或上下文相關(guān)性,所以通常用于文本翻譯中。Bian 等[20]提出使用堆疊式雙向長(zhǎng)短期記憶神經(jīng)網(wǎng)絡(luò)(SBLSTM)來估算電池SOC,并與使用了同樣的測(cè)試數(shù)據(jù)的LSTM[11]和GRU[13]模型進(jìn)行了對(duì)比。每個(gè)模型都被訓(xùn)練了五次,然后取這五次誤差的平均值來比較。從結(jié)果來看,在基于0、10 和25 ℃三種環(huán)境溫度的US06 和HWET 工況下,SBLSTM 估算精度是最高的,其中25 ℃下的MAE 只有0.46%。雖然這些模型使用了同樣的測(cè)試數(shù)據(jù),但是由于不知道文獻(xiàn)[11]和文獻(xiàn)[13]的作者是否也對(duì)他們的模型進(jìn)行了重復(fù)訓(xùn)練。因此,不同論文中的模型交叉對(duì)比還是很困難的。本文同樣比較了不同的隱藏層數(shù)和隱藏層神經(jīng)元數(shù)對(duì)模型估算結(jié)果的影響,并最終確定最優(yōu)的模型結(jié)構(gòu)擁有兩個(gè)各含有64 個(gè)神經(jīng)元的BLSTM 層,相當(dāng)于四個(gè)單向LSTM 堆疊層,并含有130 000 個(gè)可學(xué)習(xí)參數(shù),即結(jié)構(gòu)中所有權(quán)重和偏差的總和。過多的學(xué)習(xí)參數(shù)不僅會(huì)導(dǎo)致訓(xùn)練時(shí)間成倍的增加,還容易導(dǎo)致訓(xùn)練中出現(xiàn)過擬合現(xiàn)象,這極大限制了這種方法的實(shí)際應(yīng)用。

    3 擴(kuò)展循環(huán)神經(jīng)網(wǎng)絡(luò)算法

    3.1 與非循環(huán)神經(jīng)網(wǎng)絡(luò)結(jié)合

    為了進(jìn)一步提高循環(huán)神經(jīng)網(wǎng)絡(luò)估算電池SOC 的精度,有學(xué)者提出了將循環(huán)神經(jīng)網(wǎng)絡(luò)算法與其他非循環(huán)神經(jīng)網(wǎng)絡(luò)算法結(jié)合的想法,使它們可以相互彌補(bǔ),充分發(fā)揮自己的優(yōu)點(diǎn),從而提高估算精度。Song 等[21]提出了卷積神經(jīng)網(wǎng)絡(luò)(CNN)和LSTM 復(fù)合模型來估算SOC。這個(gè)復(fù)合模型綜合了CNN 和LSTM 模型的優(yōu)點(diǎn),其中CNN 用于從輸入數(shù)據(jù)中提取空間特征,而LSTM 用于提取時(shí)間特征。在環(huán)境溫度0~50 ℃下基于FUDS、DST、US06 三種工況的LFP 電池?cái)?shù)據(jù)被用來測(cè)試模型的性能和針對(duì)未知SOC 初始狀態(tài)的魯棒性。實(shí)驗(yàn)結(jié)果證明了CNN-LSTM 模型可以很好地捕捉電池SOC 和測(cè)量變量(電壓,電流和溫度)之間的非線性關(guān)系,并且比單一的LSTM和CNN 具有更優(yōu)的追蹤性能。在初始SOC 未知的情況下,CNN-LSTM 模型可以很快地收斂到SOC 真實(shí)值,然后保持平滑準(zhǔn)確的估算結(jié)果,最大MAE 和RMSE 分別在1%和2%以下。此外,該復(fù)合方法在不同的環(huán)境溫度下仍可以準(zhǔn)確地估算電池SOC,最大MAE 在1.5%以內(nèi),最大RMSE 在2%以內(nèi)。無獨(dú)有偶,Huang 等[22]提出了CNN 和BGRU 復(fù)合模型。雖然在常溫下,該模型的MAE 和RMSE 均小于2%,但是在0 ℃,MAE 和RMSE 分別是3.04%和4.36%。而且與其他循環(huán)神經(jīng)網(wǎng)絡(luò)一樣,該模型也有需要確定的超參數(shù)。為了確定模型的最優(yōu)結(jié)構(gòu),作者根據(jù)經(jīng)驗(yàn)手動(dòng)設(shè)定了幾組不同神經(jīng)元數(shù)組合并進(jìn)行估算對(duì)比,但是這些超參數(shù)是否最優(yōu),仍值得商榷。

    3.2 與卡爾曼濾波方法結(jié)合

    無跡卡爾曼濾波算法(UKF)是一種常用于去除數(shù)據(jù)流中的噪聲的技術(shù)。與EKF 相比,UKF 在具有相同復(fù)雜度等級(jí)的非線性系統(tǒng)中表現(xiàn)出更好的性能,并且易于實(shí)現(xiàn)。Yang 等[23]提出LSTM 和UKF 結(jié)合的模型(LSTM-UKF)來估算電池SOC,整個(gè)算法的框架如圖6 所示。它包括離線訓(xùn)練和在線測(cè)試兩個(gè)過程:離線訓(xùn)練階段,使用由0~50 ℃,間隔10 ℃的環(huán)境溫度下DST、US06 兩種工況循環(huán)數(shù)據(jù)來訓(xùn)練LSTM,輸入變量有電壓、電流和電池溫度,輸出變量則是SOC;在線測(cè)試階段,UKF 被用來過濾掉訓(xùn)練好的LSTM 的輸出噪聲,從而提高SOC 估算精度,其中狀態(tài)變量為由庫(kù)侖公式得到的SOC,輸出變量為L(zhǎng)STM 估算得到的SOC。由實(shí)驗(yàn)結(jié)果來看,在不同的環(huán)境溫度下,相比單一的LSTM,LSTM-UKF 估算精度至少提高50%,MAE 和RMSE 均在1.1%以內(nèi)。此外,作者還將其與其他機(jī)器學(xué)習(xí)法NN,SVM 和GPR,在常溫下進(jìn)行了對(duì)比,LSTM-UKF 在估算精度上完勝而且估算時(shí)間也只需1.20 s。除去EKF 和UKF 外,還有很多優(yōu)秀的卡爾曼濾波算法的進(jìn)階版,比如容積卡爾曼濾波算法(CKF),相比前兩者,它在高階系統(tǒng)中的表現(xiàn)更加穩(wěn)定和高效。Tian 等[24]提出了LSTM 和ACKF 結(jié)合的復(fù)合模型(LSTM-ACKF)。在這里,ACKF 與文獻(xiàn)[23]中UKF 一樣,也是用來優(yōu)化LSTM 模型的估算結(jié)果。不過與后者不同的,ACKF 具有自適應(yīng)功能,換言之就是,可以在優(yōu)化過程中自動(dòng)更新過程噪聲方差Q 和測(cè)量誤差方差R。首先,作者使用不同環(huán)境溫度下US06 和FUDS兩種工況數(shù)據(jù)來驗(yàn)證LSTM-ACKF 的估算精度并引入LSTM和LSTM-CKF 兩種模型作對(duì)比。結(jié)果顯示,LSTM 的估算誤差RMSE 均大于2.8%,最大RMSE 達(dá)到了3.5%,LSTM-CKF的估算誤差RMSE 均大于2.2%,最大RMSE 為3.2%,而LSTM-ACKF 的RMSE 均小于2.2%,最小RMSE 只有0.9%。然后,在10 ℃的FUDS 工況下,使用不同的Q 和R 初始值來驗(yàn)證LSTM-ACKF 的抗干擾性。結(jié)果顯示,LSTM-CKF 的估算結(jié)果非常依賴Q 和R 的初始值,尤其是Q 很大且R 很小的時(shí)候,此時(shí)的誤差非常大,反觀LSTM-ACKF 的估算誤差幾乎不受Q 和R 的影響。最后,使用不同的SOC 初始值來驗(yàn)證LSTM-ACKF 的收斂性。結(jié)果顯示,即使SOC 初始值誤差達(dá)到40%,該模型也能在短短20 s 內(nèi),收斂到正確值。

    圖6 LSTM-UKF框架[23]

    3.3 與時(shí)間序列模型結(jié)合

    Caliwag 等[25]提出了時(shí)間序列模型(VARMA)和LSTM 結(jié)合的算法來預(yù)測(cè)電動(dòng)摩托車的鋰離子電池電壓和SOC,其中輸入由電動(dòng)機(jī)速度,輸入功率和扭矩以及電池評(píng)估電壓,電流和溫度組合而成。該復(fù)合模型估算電池SOC 的原理是:首先使用VARMA 預(yù)測(cè)電池輸出電流的線性分量,然后用LSTM預(yù)測(cè)電池輸出電流的非線性分量,最后將它們加起來就是預(yù)測(cè)的電池輸出電流,將預(yù)測(cè)的輸出電流代入到安時(shí)積分法中即可獲得預(yù)測(cè)的SOC。作者僅使用韓國(guó)的駕駛循環(huán)CVS-40在0 和25 ℃下測(cè)試了該模型,用于訓(xùn)練模型的數(shù)據(jù)直接來自駕駛摩托車。文中沒有提供電池的有關(guān)信息,LSTM 的結(jié)構(gòu)也沒有介紹,而且尚不清楚在有更大的數(shù)據(jù)集時(shí)是否有必要將VARMA 與LSTM 結(jié)合使用。

    4 討論

    本節(jié)總結(jié)了本文中提到的估算電池SOC 的循環(huán)神經(jīng)網(wǎng)絡(luò)算法的輸入變量、數(shù)據(jù)預(yù)處理、循環(huán)神經(jīng)網(wǎng)絡(luò)(LSTM,GRU,NARX)的層數(shù)及其神經(jīng)元數(shù)、優(yōu)化器、損失函數(shù)、測(cè)試文件、研究的溫度、使用的電池類型以及估算誤差,為讀者未來使用循環(huán)神經(jīng)網(wǎng)絡(luò)方法估算電池SOC 提供理論指導(dǎo)。從表4、表5 和表6 中可以發(fā)現(xiàn),大多數(shù)研究都使用了相同的模型輸入變量,如電壓、電流和電池溫度,只有極少數(shù)使用了平均值或其他的計(jì)算值。訓(xùn)練前進(jìn)行數(shù)據(jù)預(yù)處理的模型大概占了一半,從理論上講,數(shù)據(jù)預(yù)處理可以提高模型估算電池SOC的精度,可惜暫時(shí)沒有論文對(duì)此進(jìn)行討論。另外大多數(shù)模型都只使用了一層循環(huán)神經(jīng)網(wǎng)絡(luò),可見在估算電池SOC 這個(gè)領(lǐng)域,使用多層循環(huán)神經(jīng)網(wǎng)絡(luò)并不一定能帶來更好的結(jié)果,相反還會(huì)增加訓(xùn)練難度和時(shí)間,而且神經(jīng)元數(shù)也不用太多。對(duì)于循環(huán)神經(jīng)網(wǎng)絡(luò)來說,優(yōu)化器的選擇有很多,比如RMSprop、LM、SGD、Adamax、Nadam 等。不過在估算電池SOC 方面,主流的優(yōu)化器還是Adam。至于損失函數(shù)的選擇,只有MAE、RMSE 和MSE 三種,其中MSE 應(yīng)用的最多。對(duì)于測(cè)試數(shù)據(jù)的獲取,絕大部分研究都選擇了標(biāo)準(zhǔn)的汽車工況和可變的環(huán)境溫度,其余的則選擇了簡(jiǎn)單的恒流和脈沖電流循環(huán)或固定溫度。一般來說,模型的測(cè)試數(shù)據(jù)越復(fù)雜,其呈現(xiàn)的估算精度可信度越高。從估算結(jié)果來看,基于RNN 的估算方法非常適合于電池SOC 估算領(lǐng)域。

    表7 對(duì)比了方法之間的優(yōu)缺點(diǎn),并進(jìn)行了詳細(xì)的論述,為未來循環(huán)神經(jīng)網(wǎng)絡(luò)法估算SOC 的發(fā)展提供理論指導(dǎo)。

    5 結(jié)論

    本文介紹了幾種基于循環(huán)神經(jīng)網(wǎng)絡(luò)的鋰離子電池SOC估算方法,并將它們分為簡(jiǎn)單循環(huán)神經(jīng)網(wǎng)絡(luò)算法、深度循環(huán)神經(jīng)網(wǎng)絡(luò)算法以及擴(kuò)展循環(huán)神經(jīng)網(wǎng)絡(luò)算法。簡(jiǎn)單循環(huán)神經(jīng)網(wǎng)絡(luò)算法雖然在某種條件下也能獲得較好的估算精度,不過容易受到模型超參數(shù)的影響,而目前并沒有成熟的算法來確定超參數(shù);深度循環(huán)神經(jīng)網(wǎng)絡(luò)算法通過增加模型中循環(huán)神經(jīng)層數(shù)來提高估算精度,不過模型的學(xué)習(xí)參數(shù)(權(quán)重和偏值)數(shù)量也隨之增多,從而增大了模型過擬合的風(fēng)險(xiǎn)。因此,為了讓循環(huán)神經(jīng)網(wǎng)絡(luò)可以在各種條件下都能很精確地估算電池SOC,除了需要開發(fā)一種成熟的算法來確定循環(huán)神經(jīng)網(wǎng)絡(luò)的超參數(shù)之外,還需要研究將循環(huán)神經(jīng)網(wǎng)絡(luò)法和其他算法結(jié)合起來,取長(zhǎng)補(bǔ)短,實(shí)現(xiàn)鋰離子電池SOC 在不同條件下的快速、準(zhǔn)確的估計(jì)。為此,研究人員仍需要不斷地探索和實(shí)踐。

    表4 簡(jiǎn)單循環(huán)神經(jīng)網(wǎng)絡(luò)模型超參數(shù)比較

    表5 深度循環(huán)神經(jīng)網(wǎng)絡(luò)模型超參數(shù)比較

    表6 擴(kuò)展循環(huán)神經(jīng)網(wǎng)絡(luò)模型超參數(shù)比較

    表7 基于RNN 的鋰離子電池SOC 估算方法優(yōu)缺點(diǎn)

    猜你喜歡
    優(yōu)化模型
    一半模型
    超限高層建筑結(jié)構(gòu)設(shè)計(jì)與優(yōu)化思考
    民用建筑防煙排煙設(shè)計(jì)優(yōu)化探討
    關(guān)于優(yōu)化消防安全告知承諾的一些思考
    一道優(yōu)化題的幾何解法
    由“形”啟“數(shù)”優(yōu)化運(yùn)算——以2021年解析幾何高考題為例
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    3D打印中的模型分割與打包
    FLUKA幾何模型到CAD幾何模型轉(zhuǎn)換方法初步研究
    给我免费播放毛片高清在线观看| 亚洲中文av在线| 嫩草影视91久久| 久久精品91蜜桃| 全区人妻精品视频| 黄色日韩在线| 天天一区二区日本电影三级| xxx96com| 欧美最黄视频在线播放免费| 欧美日本视频| 99国产精品一区二区三区| 亚洲欧美精品综合一区二区三区| 国产精品一区二区精品视频观看| 熟女电影av网| 黑人巨大精品欧美一区二区mp4| 亚洲精品在线美女| 亚洲电影在线观看av| 成人av一区二区三区在线看| 三级国产精品欧美在线观看 | 狠狠狠狠99中文字幕| 观看美女的网站| 国产伦人伦偷精品视频| 日日干狠狠操夜夜爽| 午夜免费激情av| 日韩有码中文字幕| aaaaa片日本免费| 99久久精品国产亚洲精品| 色综合欧美亚洲国产小说| 级片在线观看| 久久热在线av| 色综合婷婷激情| 特级一级黄色大片| 欧美不卡视频在线免费观看| 成人国产综合亚洲| 又紧又爽又黄一区二区| 精品国产超薄肉色丝袜足j| ponron亚洲| 亚洲乱码一区二区免费版| 亚洲自拍偷在线| 一进一出好大好爽视频| 88av欧美| 听说在线观看完整版免费高清| 久久久久久大精品| 深夜精品福利| 97超视频在线观看视频| 亚洲熟妇熟女久久| 国产精品爽爽va在线观看网站| 午夜福利免费观看在线| 麻豆国产97在线/欧美| 日本精品一区二区三区蜜桃| 99热只有精品国产| 丰满人妻一区二区三区视频av | 激情在线观看视频在线高清| 成人高潮视频无遮挡免费网站| 国产免费男女视频| 久久久久国产精品人妻aⅴ院| 国产精品1区2区在线观看.| 亚洲欧美日韩东京热| 一级黄色大片毛片| 国产精品亚洲一级av第二区| 成人特级黄色片久久久久久久| 国产视频一区二区在线看| 亚洲自拍偷在线| 不卡一级毛片| 亚洲真实伦在线观看| 美女大奶头视频| 久久精品夜夜夜夜夜久久蜜豆| 久久精品国产清高在天天线| 国产欧美日韩一区二区三| 久久精品国产亚洲av香蕉五月| 他把我摸到了高潮在线观看| 日韩欧美三级三区| 久久久久久大精品| 欧美zozozo另类| 亚洲av熟女| 亚洲国产精品成人综合色| 美女cb高潮喷水在线观看 | 亚洲人成网站在线播放欧美日韩| 一个人看的www免费观看视频| 身体一侧抽搐| 国产av在哪里看| 美女大奶头视频| 色综合亚洲欧美另类图片| 亚洲精品美女久久av网站| 男女床上黄色一级片免费看| 久久精品影院6| 听说在线观看完整版免费高清| 日韩av在线大香蕉| 最近最新免费中文字幕在线| 久久九九热精品免费| 亚洲五月婷婷丁香| 男插女下体视频免费在线播放| 色精品久久人妻99蜜桃| 99国产精品一区二区蜜桃av| 久久久久久久精品吃奶| 在线观看午夜福利视频| 久久精品91无色码中文字幕| 99在线视频只有这里精品首页| 国产成人aa在线观看| 首页视频小说图片口味搜索| 91麻豆av在线| 法律面前人人平等表现在哪些方面| 国产av一区在线观看免费| 人妻丰满熟妇av一区二区三区| 国产精品电影一区二区三区| 不卡av一区二区三区| 亚洲欧美精品综合一区二区三区| 国产一区二区三区视频了| 日韩精品青青久久久久久| 欧美av亚洲av综合av国产av| 久久天堂一区二区三区四区| 天堂√8在线中文| 国产精品av久久久久免费| 精品久久久久久久末码| 两性夫妻黄色片| 99热6这里只有精品| 精品日产1卡2卡| 脱女人内裤的视频| 香蕉久久夜色| 一区二区三区激情视频| 亚洲av成人不卡在线观看播放网| 黄色女人牲交| 午夜福利免费观看在线| 精品国产乱码久久久久久男人| www.www免费av| 18禁黄网站禁片午夜丰满| 欧美中文综合在线视频| 精品一区二区三区av网在线观看| 亚洲 欧美一区二区三区| 国产亚洲精品av在线| 国产精品一区二区精品视频观看| 国产一区二区三区视频了| 国产精品久久久久久亚洲av鲁大| 午夜福利欧美成人| 国产亚洲欧美在线一区二区| 欧美zozozo另类| 一级作爱视频免费观看| 国产v大片淫在线免费观看| 真实男女啪啪啪动态图| 国产亚洲精品av在线| 欧美不卡视频在线免费观看| 国产高清视频在线播放一区| 日本与韩国留学比较| 色综合亚洲欧美另类图片| 欧美性猛交黑人性爽| www.www免费av| 青草久久国产| 免费在线观看亚洲国产| 欧美日韩亚洲国产一区二区在线观看| 人妻丰满熟妇av一区二区三区| 少妇的逼水好多| 国产亚洲精品av在线| 亚洲av成人精品一区久久| 午夜精品一区二区三区免费看| 国产aⅴ精品一区二区三区波| 中文字幕最新亚洲高清| 国产探花在线观看一区二区| 中文字幕高清在线视频| 免费av不卡在线播放| 一进一出抽搐gif免费好疼| 亚洲成人中文字幕在线播放| 久久亚洲精品不卡| 欧美一区二区国产精品久久精品| 亚洲av熟女| 热99re8久久精品国产| 日韩欧美在线二视频| 久久精品国产综合久久久| 国产精品亚洲一级av第二区| 欧美av亚洲av综合av国产av| 亚洲男人的天堂狠狠| 久久午夜综合久久蜜桃| 9191精品国产免费久久| 国产精品98久久久久久宅男小说| 国产精品综合久久久久久久免费| 中文字幕人成人乱码亚洲影| 99精品欧美一区二区三区四区| 亚洲无线在线观看| 久久性视频一级片| 亚洲第一电影网av| 熟女电影av网| 日韩国内少妇激情av| 很黄的视频免费| 国产欧美日韩精品一区二区| 欧美不卡视频在线免费观看| 草草在线视频免费看| 久久久久久久精品吃奶| 婷婷丁香在线五月| 又黄又粗又硬又大视频| 免费在线观看成人毛片| 99精品在免费线老司机午夜| 黄色 视频免费看| 亚洲国产精品999在线| 国产精品一区二区免费欧美| av片东京热男人的天堂| 欧美激情久久久久久爽电影| 亚洲电影在线观看av| 亚洲 国产 在线| 成年女人永久免费观看视频| 久9热在线精品视频| 校园春色视频在线观看| 欧美高清成人免费视频www| 午夜福利在线观看免费完整高清在 | 午夜影院日韩av| 日韩三级视频一区二区三区| 最近在线观看免费完整版| 天堂√8在线中文| 香蕉av资源在线| www.999成人在线观看| 久9热在线精品视频| 可以在线观看毛片的网站| 夜夜看夜夜爽夜夜摸| 国产一区二区激情短视频| 黄色片一级片一级黄色片| svipshipincom国产片| 国产成人av激情在线播放| 午夜亚洲福利在线播放| 免费看十八禁软件| 欧美中文日本在线观看视频| 精品国内亚洲2022精品成人| 亚洲精品一卡2卡三卡4卡5卡| 嫩草影视91久久| 国产精品九九99| 成年人黄色毛片网站| 国产成人啪精品午夜网站| 免费在线观看视频国产中文字幕亚洲| 日韩免费av在线播放| 91麻豆精品激情在线观看国产| а√天堂www在线а√下载| 国产v大片淫在线免费观看| 少妇熟女aⅴ在线视频| 18禁黄网站禁片免费观看直播| 欧美激情在线99| 精品日产1卡2卡| 亚洲最大成人中文| 日本黄大片高清| 宅男免费午夜| 久久天躁狠狠躁夜夜2o2o| 欧美国产日韩亚洲一区| h日本视频在线播放| 国产不卡一卡二| 俄罗斯特黄特色一大片| 最近在线观看免费完整版| 九色国产91popny在线| 欧美精品啪啪一区二区三区| 99re在线观看精品视频| 国内毛片毛片毛片毛片毛片| 在线免费观看不下载黄p国产 | 在线观看美女被高潮喷水网站 | 欧美日韩一级在线毛片| xxx96com| 久久香蕉精品热| 久久天堂一区二区三区四区| 偷拍熟女少妇极品色| 99riav亚洲国产免费| 欧美一级毛片孕妇| 一区二区三区国产精品乱码| 夜夜看夜夜爽夜夜摸| 国产蜜桃级精品一区二区三区| 亚洲五月天丁香| 午夜成年电影在线免费观看| 亚洲欧美日韩东京热| 久久精品国产亚洲av香蕉五月| 国产一区二区三区在线臀色熟女| 亚洲精品久久国产高清桃花| 亚洲精品456在线播放app | 成人三级做爰电影| 看黄色毛片网站| 人妻夜夜爽99麻豆av| 床上黄色一级片| 别揉我奶头~嗯~啊~动态视频| 怎么达到女性高潮| 国内精品美女久久久久久| 中文字幕人妻丝袜一区二区| www.www免费av| 国产高清videossex| 亚洲成人免费电影在线观看| xxx96com| 国产精品一区二区精品视频观看| 日韩成人在线观看一区二区三区| 国产一区二区三区在线臀色熟女| 午夜激情福利司机影院| 男女那种视频在线观看| 小说图片视频综合网站| 亚洲激情在线av| 丝袜人妻中文字幕| 色尼玛亚洲综合影院| 欧美精品啪啪一区二区三区| 亚洲自拍偷在线| 亚洲性夜色夜夜综合| 两人在一起打扑克的视频| 好看av亚洲va欧美ⅴa在| 99国产极品粉嫩在线观看| 精品国产乱码久久久久久男人| 狂野欧美激情性xxxx| 欧美日韩国产亚洲二区| 男女午夜视频在线观看| 偷拍熟女少妇极品色| 桃色一区二区三区在线观看| 亚洲无线观看免费| 久久国产乱子伦精品免费另类| 精品免费久久久久久久清纯| 精品日产1卡2卡| 国产毛片a区久久久久| 最新中文字幕久久久久 | 别揉我奶头~嗯~啊~动态视频| 久久久久精品国产欧美久久久| 搞女人的毛片| 亚洲乱码一区二区免费版| 少妇裸体淫交视频免费看高清| 在线看三级毛片| 韩国av一区二区三区四区| 国产精品影院久久| 99国产综合亚洲精品| 久久久国产成人免费| 真人一进一出gif抽搐免费| 亚洲国产欧美人成| 国产精品综合久久久久久久免费| 精品国产乱子伦一区二区三区| 极品教师在线免费播放| 国产欧美日韩一区二区三| 国产精品,欧美在线| 午夜福利视频1000在线观看| 久久国产精品人妻蜜桃| 日韩中文字幕欧美一区二区| 长腿黑丝高跟| 国产成人一区二区三区免费视频网站| 99久久99久久久精品蜜桃| 亚洲乱码一区二区免费版| 久久中文字幕人妻熟女| 日本黄色片子视频| 天天一区二区日本电影三级| 一进一出好大好爽视频| 免费观看的影片在线观看| 精品久久久久久久久久免费视频| 夜夜看夜夜爽夜夜摸| 长腿黑丝高跟| 国产成人一区二区三区免费视频网站| 香蕉久久夜色| 亚洲国产看品久久| а√天堂www在线а√下载| 99久久成人亚洲精品观看| 三级国产精品欧美在线观看 | 国产高清三级在线| 国产一区二区三区视频了| 亚洲国产欧洲综合997久久,| 禁无遮挡网站| 脱女人内裤的视频| 亚洲乱码一区二区免费版| 久久亚洲真实| 久久伊人香网站| 亚洲成人久久性| 久久久国产精品麻豆| 国产成人欧美在线观看| 激情在线观看视频在线高清| 亚洲,欧美精品.| 久久这里只有精品19| 搡老熟女国产l中国老女人| 中文字幕人成人乱码亚洲影| 久久亚洲真实| 国产在线精品亚洲第一网站| av天堂在线播放| 99久久无色码亚洲精品果冻| 精品国内亚洲2022精品成人| 国内少妇人妻偷人精品xxx网站 | 免费观看人在逋| 老司机福利观看| 一级黄色大片毛片| 天堂动漫精品| 热99re8久久精品国产| 亚洲色图av天堂| 久久精品影院6| 国产亚洲欧美98| 一二三四在线观看免费中文在| 在线免费观看的www视频| 男人和女人高潮做爰伦理| 欧美最黄视频在线播放免费| 老汉色∧v一级毛片| 亚洲av电影在线进入| 美女扒开内裤让男人捅视频| 国产精品久久久久久久电影 | 国产主播在线观看一区二区| 18美女黄网站色大片免费观看| 久久人妻av系列| 久久久久久九九精品二区国产| 51午夜福利影视在线观看| 亚洲天堂国产精品一区在线| 身体一侧抽搐| 我要搜黄色片| 国产激情久久老熟女| 老熟妇仑乱视频hdxx| 亚洲精品中文字幕一二三四区| 99久久成人亚洲精品观看| 不卡一级毛片| 女生性感内裤真人,穿戴方法视频| 精品一区二区三区四区五区乱码| netflix在线观看网站| 男女做爰动态图高潮gif福利片| 久久九九热精品免费| 亚洲无线观看免费| 国产一区二区在线观看日韩 | 亚洲美女黄片视频| 午夜福利在线观看吧| 国产av麻豆久久久久久久| 国产伦精品一区二区三区四那| 在线观看66精品国产| 国产成年人精品一区二区| 中文字幕精品亚洲无线码一区| 国产精品99久久久久久久久| 91麻豆精品激情在线观看国产| 成人性生交大片免费视频hd| 久久精品国产综合久久久| 国产伦精品一区二区三区视频9 | 亚洲成人精品中文字幕电影| 日本a在线网址| 国内精品久久久久精免费| 色在线成人网| 男女视频在线观看网站免费| 美女被艹到高潮喷水动态| 久久国产乱子伦精品免费另类| 国产成人影院久久av| 国产熟女xx| 免费av毛片视频| 真人一进一出gif抽搐免费| 国产精品亚洲美女久久久| 久久天堂一区二区三区四区| 狂野欧美白嫩少妇大欣赏| 网址你懂的国产日韩在线| 麻豆国产97在线/欧美| 又爽又黄无遮挡网站| 亚洲国产高清在线一区二区三| 欧美日韩乱码在线| 一级毛片女人18水好多| 97人妻精品一区二区三区麻豆| 黄片大片在线免费观看| 亚洲自偷自拍图片 自拍| 久久久久久久久久黄片| 欧美中文日本在线观看视频| 成年女人看的毛片在线观看| 天天添夜夜摸| 久久国产精品人妻蜜桃| 午夜影院日韩av| 一本久久中文字幕| 国产精品乱码一区二三区的特点| 久久精品影院6| 日韩欧美在线二视频| 12—13女人毛片做爰片一| 国产精品久久视频播放| 真人做人爱边吃奶动态| 国产蜜桃级精品一区二区三区| 欧美丝袜亚洲另类 | 欧美绝顶高潮抽搐喷水| 香蕉av资源在线| 日本黄色片子视频| 老汉色∧v一级毛片| 69av精品久久久久久| 国产精品综合久久久久久久免费| 无遮挡黄片免费观看| 99精品久久久久人妻精品| 国内精品一区二区在线观看| 亚洲中文av在线| 久久久久精品国产欧美久久久| 国产av在哪里看| 精品欧美国产一区二区三| 天天一区二区日本电影三级| 久久精品综合一区二区三区| 久久久久精品国产欧美久久久| 午夜福利18| 一进一出抽搐动态| 亚洲精品在线美女| 一个人观看的视频www高清免费观看 | 久久久国产欧美日韩av| 久久午夜亚洲精品久久| 变态另类成人亚洲欧美熟女| 日本a在线网址| 国产亚洲精品一区二区www| 欧美性猛交黑人性爽| 99热精品在线国产| 精品乱码久久久久久99久播| 高清在线国产一区| avwww免费| 亚洲av电影不卡..在线观看| 成人三级黄色视频| 欧美丝袜亚洲另类 | 一个人免费在线观看电影 | 成人性生交大片免费视频hd| 国产成人系列免费观看| 黑人欧美特级aaaaaa片| 精品久久久久久久毛片微露脸| 老司机深夜福利视频在线观看| 欧美黑人巨大hd| 每晚都被弄得嗷嗷叫到高潮| 老司机午夜福利在线观看视频| 国产精品一及| 国产亚洲精品综合一区在线观看| a级毛片在线看网站| 一区二区三区高清视频在线| 91麻豆av在线| 亚洲av成人av| 天天一区二区日本电影三级| 成年人黄色毛片网站| 国产男靠女视频免费网站| svipshipincom国产片| 国产av不卡久久| 欧美绝顶高潮抽搐喷水| 久久性视频一级片| 国产野战对白在线观看| 国产aⅴ精品一区二区三区波| 亚洲无线在线观看| 黄色日韩在线| 三级毛片av免费| avwww免费| 国产精品电影一区二区三区| av视频在线观看入口| 欧美日韩黄片免| 国产精品亚洲av一区麻豆| 免费电影在线观看免费观看| 91九色精品人成在线观看| 麻豆av在线久日| 91字幕亚洲| 国产精品久久电影中文字幕| 变态另类丝袜制服| 曰老女人黄片| 可以在线观看的亚洲视频| 午夜福利欧美成人| 国产激情欧美一区二区| 91久久精品国产一区二区成人 | 国产三级在线视频| 国产 一区 欧美 日韩| 最近视频中文字幕2019在线8| 在线a可以看的网站| 中文字幕精品亚洲无线码一区| 亚洲片人在线观看| 99国产精品99久久久久| 精品久久久久久久末码| 久久欧美精品欧美久久欧美| 久久久久精品国产欧美久久久| 淫秽高清视频在线观看| 小说图片视频综合网站| 日本撒尿小便嘘嘘汇集6| 亚洲激情在线av| 黄色丝袜av网址大全| 一个人免费在线观看的高清视频| 在线a可以看的网站| 免费在线观看亚洲国产| 免费电影在线观看免费观看| 国产伦一二天堂av在线观看| 国产亚洲av嫩草精品影院| 午夜影院日韩av| 国产精品一区二区免费欧美| 搡老妇女老女人老熟妇| 99久久久亚洲精品蜜臀av| 亚洲av电影不卡..在线观看| 成年女人永久免费观看视频| 国产精品98久久久久久宅男小说| x7x7x7水蜜桃| 黄片小视频在线播放| 无限看片的www在线观看| 午夜精品在线福利| 午夜福利欧美成人| 99久久精品一区二区三区| a级毛片在线看网站| 久久久久亚洲av毛片大全| 一本精品99久久精品77| 在线观看美女被高潮喷水网站 | 真实男女啪啪啪动态图| 国产伦人伦偷精品视频| 成人av一区二区三区在线看| www.精华液| 国产一区二区三区在线臀色熟女| 国产亚洲精品久久久久久毛片| 亚洲一区高清亚洲精品| 国产欧美日韩精品亚洲av| 手机成人av网站| 久久这里只有精品19| 国产精品免费一区二区三区在线| 一个人看视频在线观看www免费 | 可以在线观看毛片的网站| 一个人观看的视频www高清免费观看 | АⅤ资源中文在线天堂| 欧美日韩精品网址| aaaaa片日本免费| 88av欧美| 手机成人av网站| 精品一区二区三区四区五区乱码| 国产私拍福利视频在线观看| 欧美绝顶高潮抽搐喷水| 午夜精品在线福利| 搡老岳熟女国产| 亚洲人与动物交配视频| 91麻豆精品激情在线观看国产| 免费在线观看视频国产中文字幕亚洲| 精品熟女少妇八av免费久了| 欧美午夜高清在线| 超碰成人久久| 黄色片一级片一级黄色片| 丰满人妻一区二区三区视频av | 亚洲电影在线观看av| 成人三级做爰电影| 噜噜噜噜噜久久久久久91| 亚洲中文字幕一区二区三区有码在线看 | 日韩精品青青久久久久久| 两个人看的免费小视频| 日日干狠狠操夜夜爽| 神马国产精品三级电影在线观看| 91老司机精品| 久久亚洲真实| 国产三级在线视频| 日日夜夜操网爽| 十八禁人妻一区二区| 国产真人三级小视频在线观看| 人人妻人人看人人澡| 精品欧美国产一区二区三| 久久久成人免费电影| 亚洲va日本ⅴa欧美va伊人久久| 午夜激情福利司机影院|