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

    基于一致性模型的梯次利用鋰離子電池組能量利用率估計(jì)方法

    2021-05-28 00:11:12李章溢
    電工技術(shù)學(xué)報(bào) 2021年10期
    關(guān)鍵詞:內(nèi)阻電池組倍率

    趙 偉 閔 婕 李章溢 孫 瑞 姜 研

    基于一致性模型的梯次利用鋰離子電池組能量利用率估計(jì)方法

    趙 偉1閔 婕2李章溢2孫 瑞2姜 研3,4,5

    (1. 廣東電網(wǎng)有限責(zé)任公司電力科學(xué)研究院 廣州 510080 2. 深圳市欣旺達(dá)綜合能源服務(wù)有限公司 深圳 518000 3. 欣旺達(dá)電子股份有限公司 深圳 518108 4. 華南理工大學(xué)電力學(xué)院 廣州 510641 5. 國家能源主動(dòng)配電網(wǎng)技術(shù)研發(fā)中心(北京交通大學(xué)) 北京 100044)

    鋰離子電池的不一致性導(dǎo)致電池組的容量和壽命等指標(biāo)顯著小于電池單體,且該問題在梯次利用的應(yīng)用場景下更為嚴(yán)重。電池組老化后,電池組內(nèi)電池不一致性參數(shù)之間并非相互獨(dú)立,而是相互耦合,因此造成電池組一致性建模困難,降低了電池組的能量利用率的估計(jì)精度。針對(duì)梯次利用電池篩選成組和成組后的狀態(tài)評(píng)估,該文提出一種基于Copula理論的電池組一致性建模方法,同時(shí)實(shí)現(xiàn)了電池參數(shù)分布特性和參數(shù)間相關(guān)性的準(zhǔn)確刻畫?;陔姵亟M一致性模型,提出一種電池組能量利用率估計(jì)方法,并進(jìn)行實(shí)驗(yàn)驗(yàn)證,估計(jì)誤差在1%以內(nèi)。最后分析電池組能量利用率的影響因素,以指導(dǎo)電池組的優(yōu)化成組和使用。

    鋰離子電池組 梯次利用 一致性 能量利用率

    0 引言

    由于鋰離子電池具有能量密度高、效率高和使用壽命長等優(yōu)點(diǎn),近年來已被廣泛地應(yīng)用于電動(dòng)汽車和儲(chǔ)能系統(tǒng)當(dāng)中[1-2]。隨著電動(dòng)汽車數(shù)量的增長,未來將有大量的鋰離子電池從電動(dòng)汽車上退役,退役電池的出路逐漸成為全行業(yè)乃至全社會(huì)關(guān)注的重點(diǎn)。鋰離子電池從電動(dòng)汽車上退役后,仍然具有70%~80%的容量。若能將退役電池應(yīng)用于梯次利用應(yīng)用場景之中,不僅可以延長電池的使用年限,發(fā)揮出電池的剩余價(jià)值,增加鋰離子電池全生命周期的經(jīng)濟(jì)效益;也可以大大減緩電池回收階段的壓力,帶來巨大的社會(huì)效益[3-4]。因此,近年來國家相繼出臺(tái)了多項(xiàng)保障政策,以推進(jìn)動(dòng)力電池梯次利用的發(fā)展,梯次利用動(dòng)力電池規(guī)模化工程應(yīng)用于儲(chǔ)能系統(tǒng)中成為技術(shù)趨勢之一[5]。

    由于鋰離子電池單體電壓和功率等級(jí)的限制,在實(shí)際使用中需要將其大量串并聯(lián)成組以滿足功率和能量的需求。電池單體出廠時(shí)一致性往往較好,但在電池組實(shí)際使用中,由于電池組溫度場分布不均勻、電池單體在電池組中所處的位置不同等因素[6],導(dǎo)致電池組老化過程中電池單體的一致性參數(shù)(如容量、內(nèi)阻和荷電狀態(tài)(State of Charge, SOC)等)快速發(fā)散,使得電池在梯次利用場景下,其一致性不及新電池,而電池的不一致直接影響電池組的使用效率,降低電池組的經(jīng)濟(jì)效益[7]。因此,亟需研究梯次利用電池組的能量利用率計(jì)算方法,以指導(dǎo)梯次利用電池使用前的篩選成組和成組后的電池組狀態(tài)評(píng)估。

    電池組能量利用率的主要影響因素是電池的一致性。在實(shí)際應(yīng)用中,測量每一個(gè)電池單體的參數(shù)將耗費(fèi)大量人力物力,使得梯次利用電池不再具備成本上的優(yōu)勢。通過適當(dāng)?shù)某闃臃椒y量少部分電池樣本,可以獲得電池參數(shù)的統(tǒng)計(jì)特性,因此學(xué)者多研究基于概率統(tǒng)計(jì)的電池一致性建模。

    當(dāng)前,針對(duì)電池組一致性的相關(guān)研究主要集中于兩方面:一方面研究電池單體參數(shù)的統(tǒng)計(jì)特性。S. F. Schuster等[8]通過將老化的電池組和新三元電池組進(jìn)行對(duì)比,發(fā)現(xiàn)電池容量和內(nèi)阻在新電池狀態(tài)下呈現(xiàn)正態(tài)分布,但在電池組老化過程中,正態(tài)分布逐漸轉(zhuǎn)變?yōu)閃eibull分布。王震坡等[9]對(duì)不一致性的發(fā)展規(guī)律進(jìn)行了研究,提出電壓的不一致性分布遵循正態(tài)分布并給出了正態(tài)分布參數(shù)的計(jì)算方法。魏五星等[10]分析了110節(jié)串聯(lián)電池電壓的一致性,認(rèn)為開路電壓和極化電壓都服從正態(tài)分布。程功 等[11]對(duì)兩輛電動(dòng)公交車電壓的數(shù)據(jù)進(jìn)行了長期跟蹤,發(fā)現(xiàn)電池的開路電壓(Open Circuit Voltage, OCV)服從正態(tài)分布,但隨著車輛的使用,正態(tài)分布的標(biāo)準(zhǔn)差不斷增大。在電池組老化的過程中,電池參數(shù)不僅呈現(xiàn)一定的統(tǒng)計(jì)規(guī)律,參數(shù)間也同時(shí)呈現(xiàn)出一定的相關(guān)性,因此另一方面研究關(guān)注于電池參數(shù)間的耦合機(jī)制。Zheng Yuejiu等[12]研究了電池容量衰退不一致與電池SOC不一致之間的關(guān)系,研究結(jié)果表明,鋰離子在負(fù)極上的損失是電池間產(chǎn)生SOC差異的主要因素。Yang Fangfang等[13]研究了庫侖效率與電池容量衰退之間的關(guān)系,庫倫效率較低的電池不僅具有更快的容量衰退速率,也同時(shí)加大了電池間SOC的不一致。但是,當(dāng)涉及電池一致性多參數(shù)聯(lián)合建模時(shí),現(xiàn)有的一致性刻畫的方法中,多認(rèn)為電池參數(shù)間相互獨(dú)立[14]或服從聯(lián)合分布[15],并沒有同時(shí)考慮參數(shù)的統(tǒng)計(jì)特性和參數(shù)間的相關(guān)性。因此本文提出了一種基于Copula理論的電池一致性建模方法,基于一致性模型實(shí)現(xiàn)了電池組能量利用率的仿真計(jì)算,通過實(shí)驗(yàn)驗(yàn)證了方法的有效性并進(jìn)一步分析了能量利用率的影響因素。

    1 基于Copula理論的一致性模型

    當(dāng)前,電池組的不一致主要表現(xiàn)為電池單體容量、內(nèi)阻和SOC的不一致。每個(gè)參數(shù)都有其特定的統(tǒng)計(jì)特征,與此同時(shí),參數(shù)之間又存在著一定的相依性。因此,本文基于Copula理論,實(shí)現(xiàn)電池參數(shù)分布特性及參數(shù)間相關(guān)性的準(zhǔn)確刻畫。

    1.1 Copula定義

    因此,Copula可被用于刻畫隨機(jī)變量之間的相依性。同樣地,如果已知隨機(jī)變量的邊緣分布函數(shù)和對(duì)應(yīng)的Copula,則可以構(gòu)造出一系列與隨機(jī)變量具有相同分布特性和變量間相依性的隨機(jī)數(shù)。

    1.2 Copula參數(shù)估計(jì)方法

    由式(3)可知,Copula理論在應(yīng)用過程中的未知參數(shù)為隨機(jī)變量的邊緣分布函數(shù)的相關(guān)參數(shù)和Copula模型參數(shù)。本文中采用基于非參數(shù)核密度估計(jì)的兩階段極大似然估計(jì)實(shí)現(xiàn)Copula參數(shù)估計(jì),即首先利用核密度估計(jì)法獲得隨機(jī)變量的累積分布函數(shù),再通過極大似然估計(jì)法估計(jì)Copula模型參數(shù)。

    其中

    對(duì)隨機(jī)變量的邊緣分布函數(shù)進(jìn)行核密度估計(jì)后,開展對(duì)Copula模型參數(shù)的極大似然估計(jì)。本文以三變量為例,其Copula函數(shù)及其概率密度函數(shù)可分別表示為

    根據(jù)隨機(jī)變量邊緣分布函數(shù)的估計(jì)方法,其聯(lián)合分布函數(shù)可表示為

    對(duì)應(yīng)的概率密度函數(shù)為

    樣本的似然函數(shù)為

    對(duì)似然函數(shù)求對(duì)數(shù)可得

    盡管根據(jù)Sklar定理,任意一多元分布函數(shù)都有其對(duì)應(yīng)的特定的Copula函數(shù),但受限于相關(guān)數(shù)學(xué)理論,目前僅提出有限個(gè)類型的Copula函數(shù),每一個(gè)類型的Copula函數(shù)適用于描述一種特定類型的變量間的相關(guān)關(guān)系。由于具備在高維空間的泛化能力,本文選取高斯Copula函數(shù)來描述變量間的相依結(jié)構(gòu)。其分布函數(shù)和概率密度函數(shù)的形式分別為

    1.3 電池組一致性模型抽樣方法

    在給定電池一致性參數(shù)的數(shù)據(jù)后,即可根據(jù)1.2節(jié)中介紹的方法確定最優(yōu)的Copula函數(shù)。因此,便可以模擬生成一系列隨機(jī)數(shù)來建立電池組一致性多元分布模型,其具體步驟如下:

    (1)基于觀測值,采用核密度估計(jì)方法,確定各隨機(jī)變量的累積分布函數(shù)(隨機(jī)變量聯(lián)合分布中的邊緣分布函數(shù))。

    (2)將輸入隨機(jī)變量(-域)通過累積分布函數(shù)變換成Copula域內(nèi)([0, 1]上均勻分布)的隨變量。

    (3)確定Copula函數(shù)的參數(shù)。

    (4)根據(jù)Copula函數(shù)生成Copula域內(nèi)的隨機(jī)數(shù)。

    (5)根據(jù)累積分布函數(shù)的反函數(shù),將Copula域內(nèi)的隨機(jī)數(shù)變換為對(duì)應(yīng)的電池一致性參數(shù)的數(shù)據(jù)。

    經(jīng)過變換后,新產(chǎn)生的隨機(jī)數(shù)將與觀測值具有相同的邊緣分布函數(shù)和隨機(jī)變量間的相依性特性。

    2 能量利用率估計(jì)方法

    電池組的能量利用率為電池組的放電可用能量與電池組額定能量的比值,即

    式中,為電池組的額定能量,即電池組理論上可儲(chǔ)存的能量最大值;為電池組的放電可用能量,即在某一電流倍率下,電池組放電SOC從100%到0的過程中的放電能量。在放電過程中,若任意一電池單體的電壓達(dá)到放電截止電壓,即認(rèn)為電池組SOC=0。本文以如圖1所示的鋰離子電池Thevenin等效電路模型為例,介紹電池組能量利用率的計(jì)算方法。

    圖1中,ocv為電池的開路電壓;U和p分別為歐姆內(nèi)阻o和極化內(nèi)阻P上的壓降;為放電電流。設(shè)為電池歐姆內(nèi)阻和極化內(nèi)阻之和。若電池極化電壓的動(dòng)態(tài)特性可忽略不計(jì),則電池的端電壓可表示為

    當(dāng)電池組對(duì)外放出能量時(shí),一部分能量在內(nèi)阻上消耗。在一個(gè)由個(gè)電池串聯(lián)組成的電池組中,設(shè)Q、R和SOCstart,分別為電池組中第個(gè)電池單體的容量、內(nèi)阻和在放電初始時(shí)刻的SOC。則電池組的放電可用能量可進(jìn)一步推導(dǎo)為

    式中,ocv,i()為電池OCV在放電過程中隨放電時(shí)間的變化;D為電池組的放電時(shí)間。通常來說,可近似認(rèn)為所有電池都具有相同的OCV-SOC曲線,即

    由于電池組內(nèi)的電池串聯(lián)連接,電池組的放電時(shí)間由所有電池中具有最小放電電量的電池決定。設(shè)pack為電池組的容量,DSOC和SOCend,i分別為第個(gè)電池在以電流放電到截止電壓下限cutoff時(shí)SOC的變化范圍和放電截止時(shí)的SOC,電池組的放電時(shí)間可表示為

    在實(shí)際應(yīng)用中,可認(rèn)為SOC已知,而放電截止時(shí)刻的SOC為

    電池組放電可用能量的定義,考慮了電池組放電過程中電流在內(nèi)阻上的損耗。而電池組的額定容量定義為電池單體可儲(chǔ)存的能量之和,反映了電池組儲(chǔ)存能量的能力,因此電池組的額定能量[11]可寫為

    式中,av為電池SOC在從100%到0變化過程中開路電壓OCV的平均值。因此,基于電池組放電可用能量和電池組額定能量的定義,電池組能量利用率為

    需要說明的是,第1節(jié)中所述的電池一致性多元分布模型為關(guān)于電池一致性參數(shù)的概率模型,而式(24)中的EUE仍為一隨機(jī)變量,且該隨機(jī)變量具有復(fù)雜的概率密度函數(shù),因此,電池組能量利用率的數(shù)值結(jié)果不能通過數(shù)學(xué)推導(dǎo)來解決。為了解決這一問題,本文使用蒙特卡羅模擬的方式來獲得能量利用率的數(shù)值結(jié)果。蒙特卡洛的基本原理思想是通過多次重復(fù)模擬的方式來獲得模型輸出的統(tǒng)計(jì)特征。對(duì)于每次模擬中的模型輸入,可根據(jù)參數(shù)的概率分布通過隨機(jī)抽樣的方式獲得。本文以電池組能量利用率模擬仿真結(jié)果的期望值作為能量利用率的估計(jì)值,在次模擬后,能量利用率的估計(jì)值可表示為

    3 驗(yàn)證與討論

    3.1 電池測試

    為了驗(yàn)證本文提出的電池組一致性建模方法的先進(jìn)性以及能量利用率估計(jì)方法的準(zhǔn)確性,本文針對(duì)磷酸鐵鋰電池組開展相關(guān)測試進(jìn)行實(shí)驗(yàn)驗(yàn)證。該電池組在一電動(dòng)乘用車上運(yùn)行超過3年的時(shí)間,電池組由95塊電池單體串聯(lián)組成,電池單體的額定容量為60A·h,在超過3年的使用中,電池組未經(jīng)過均衡維護(hù)。因此,電池間存在著容量、內(nèi)阻以及SOC的差異。整個(gè)電池測試分為電池組倍率性能測試、電池單體容量測試和電池單體內(nèi)阻測試三個(gè)部分:

    (1)電池組倍率性能測試

    在電池組倍率性能測試中,電池組以0.1進(jìn)行恒流充電,并分別以0.1、0.3、0.7和1.2進(jìn)行恒流放電。在充電與放電過程之間,電池組靜置2h。對(duì)于放電倍率大于0.1的放電測試,在電池組靜置過后,增加一個(gè)以0.1放電并靜置2h的過程,使得電池組在下一次充電之前處于電量完全放空的狀態(tài)。整個(gè)電池組倍率性能測試過程中,電池組的充/放電截止條件分別為任意一電池單體達(dá)到充/放電截止電壓(分別為3.65V和2.5V)。

    (2)電池單體容量測試

    在完成電池組倍率性能測試后,首先對(duì)電池組以0.1進(jìn)行恒流充電,截止條件為任一電池單體電壓達(dá)到充電截止電壓。然后將電池組拆解成95個(gè)電池單體,進(jìn)行電池單體的容量測試。在電池單體容量測試中,將電池以0.1進(jìn)行恒流充電與恒流放電,充放電的截止電壓分別為3.65V和2.5V,充電與放電過程之間,電池靜置2h。將電池單體容量測試中第二次的充電容量記為電池單體的容量。

    (3)電池單體內(nèi)阻測試

    本文采用復(fù)合脈沖電流法(Hybrid Pulse Power Characteristic, HPPC)對(duì)電池單體的內(nèi)阻進(jìn)行測試。首先將電池單體的SOC調(diào)至50%,靜置2h后,以1對(duì)電池進(jìn)行恒流放電10s,靜置40s后,再以0.5對(duì)電池進(jìn)行恒流充電10s。整個(gè)測試過程中電池的電流與電壓響應(yīng)如圖2所示,電池的內(nèi)阻可由歐姆定律計(jì)算得出。為了簡化計(jì)算,在本文中以10s的放電內(nèi)阻作為電池的內(nèi)阻,其具體計(jì)算方法為

    式中,Rd為放電內(nèi)阻;Id為放電電流;V0和V1分別為放電0s和10s時(shí)刻電池的端電壓。

    完成電池組倍率性能測試和電池單體的容量測試之后,可計(jì)算得出電池單體在電池組中SOC的差異。設(shè)電池單體容量測試中第一次循環(huán)過程中的充電容量為c,電池單體的容量為,則電池單體在電池組滿電狀態(tài)下的SOC為

    3.2 一致性建模與能量利用率估計(jì)精度驗(yàn)證

    圖3為3.1節(jié)所述的電池測試中電池容量、內(nèi)阻和電池組充滿電時(shí)刻SOC的矩陣散點(diǎn)圖。從圖3中參數(shù)的直方圖中可以看出,電池的容量近似正態(tài)分布,而電池內(nèi)阻和SOC兩個(gè)參數(shù)呈現(xiàn)出明顯的偏態(tài)分布。與此同時(shí),電池容量與內(nèi)阻、內(nèi)阻和SOC之間沒有展現(xiàn)出明顯的相關(guān)性,而容量和SOC之間顯示出較強(qiáng)的正相關(guān)特性。

    圖3 電池容量、內(nèi)阻和SOC的矩陣散點(diǎn)圖

    以電池的測試數(shù)據(jù)作為觀測值,根據(jù)1.2節(jié)和1.3節(jié)所述的電池組一致性建模方法,在采樣點(diǎn)數(shù)量為1 000的情況下建立的電池組一致性參數(shù)的三維散點(diǎn)圖如圖4所示。在所提出的一致性模型中,高斯Copula函數(shù)的選擇保證了所建立的模型與觀測值之間具有相同的相依性關(guān)系,而1.3節(jié)中所述的基于觀測值累積分布函數(shù)的逆函數(shù)的變換方法保證了模型中每個(gè)隨機(jī)變量都和觀測值具有相同的統(tǒng)計(jì)特性。從圖4中可以清楚地看到,SOC和容量之間存在明顯的正相關(guān)關(guān)系,表明所提出的電池一致性建模方法可同時(shí)準(zhǔn)確地刻畫參數(shù)自身的統(tǒng)計(jì)特性和參數(shù)之間的相關(guān)性。

    圖4 電池組一致性模型

    根據(jù)電池組倍率性能測試和電池單體容量測試的結(jié)果,可得到電池組在不同倍率下的放電能量,進(jìn)而獲得電池當(dāng)前狀態(tài)下能量利用率的真實(shí)值。基于電池組一致性多元分布模型,根據(jù)第2節(jié)中所述的方法對(duì)電池組能量利用率進(jìn)行估計(jì),在仿真計(jì)算中,將蒙特卡洛模擬次數(shù)設(shè)置為5 000。電池組能量利用率估計(jì)結(jié)果見表1,在4個(gè)電流倍率下,電池組能量利用率的估計(jì)結(jié)果都顯示出較高的估計(jì)精度,估計(jì)誤差在0.5%以內(nèi)。

    表1 不同放電倍率下電池組EUE估計(jì)結(jié)果

    Tab.1 EUE estimation results at different discharge current rates (%)

    3.3 能量利用率影響因素分析

    電池組的能量利用率受多種因素影響。第一,在基于電池組一致性模型的能量利用率估計(jì)中,采用了蒙特卡洛模擬的方法。蒙特卡洛模擬的次數(shù)在一定程度上影響了電池組EUE估計(jì)的精度。第二,電池組的能量利用率受電池組成組方式和使用方式的影響。當(dāng)未來大量鋰離子電池從電動(dòng)汽車上退運(yùn)后,若采用基于一致性模型的電池組能量利用率評(píng)估方法,從概率統(tǒng)計(jì)的角度去評(píng)估電池組在不同放電電流、不同串聯(lián)數(shù)下電池組的性能,將有助于為退役鋰離子電池尋找到最優(yōu)的成組及使用方式。第三,電池參數(shù)的一致性也在一定程度上影響了電池組能量利用率?;谝恢滦阅P偷碾姵亟M能量利用率評(píng)估方法可定量地分析各參數(shù)的一致性對(duì)電池組能量利用率的敏感性,進(jìn)而為確定電池組老化的主導(dǎo)因素提供理論支撐。為此,本節(jié)中將以根據(jù)3.1節(jié)中電池測試結(jié)果建立的電池組一致性模型為例,定量地分析蒙特卡洛模擬次數(shù)、電池串聯(lián)數(shù)、放電倍率和電池參數(shù)一致性對(duì)電池組能量利用率的影響。

    3.3.1 蒙特卡洛模擬次數(shù)的影響

    基于電池組一致性模型,對(duì)電池組的能量利用率進(jìn)行了5次仿真計(jì)算。在每次仿真的過程中,蒙特卡洛模擬次數(shù)設(shè)置為10 000。圖5顯示了在5次仿真計(jì)算中,電池組能量利用率的估計(jì)結(jié)果隨蒙特卡洛模擬次數(shù)的演變規(guī)律。在前幾百次的模擬過程中電池組能量利用率估計(jì)的結(jié)果波動(dòng)較大,當(dāng)模擬次數(shù)達(dá)到2 500后,5次仿真的能量利用率估計(jì)結(jié)果相差不到0.1%,模擬進(jìn)行到5 000次后,5次仿真的結(jié)果幾乎相同。因此,在后續(xù)的仿真過程中,為了節(jié)省仿真時(shí)間,統(tǒng)一將蒙特卡洛模擬次數(shù)設(shè)定為3 000。

    圖5 電池組EUE計(jì)算結(jié)果隨蒙特卡洛模擬次數(shù)的演變規(guī)律

    3.3.2 電池串聯(lián)數(shù)的影響

    圖6顯示了電池串聯(lián)數(shù)從10增加到150的過程中,電池組的能量利用率在不同放電倍率下的變化情況。在不同的放電倍率下,能量利用率都表現(xiàn)出相似的趨勢。在電池串聯(lián)數(shù)從10增加到30的過程中,電池組能量利用率下降的趨勢明顯,各個(gè)倍率下的能量利用率都下降了5%~6%。當(dāng)電池串聯(lián)數(shù)繼續(xù)增大后,能量利用率下降的速率逐漸變緩,能量利用率漸漸趨于穩(wěn)定。這意味著從概率統(tǒng)計(jì)的角度看,梯次利用電池大規(guī)模串聯(lián)成組后,電池組的能量利用率會(huì)最終趨近于一個(gè)穩(wěn)定值,能量利用率不會(huì)因串聯(lián)電池?cái)?shù)量的繼續(xù)增加而降低。

    圖6 在不同放電電流倍率下電池組EUE隨電池串聯(lián)數(shù)的演變規(guī)律

    3.3.3 放電倍率的影響

    圖7顯示了放電電流倍率從0.1逐漸增加到1.3的過程中,電池組能量利用率在不同電池串聯(lián)數(shù)下的變化情況。隨著放電倍率的增加,電池組能量利用率呈現(xiàn)出近似線性下降的規(guī)律,并沒有能量利用率迅速下降的拐點(diǎn)出現(xiàn)。當(dāng)放電電流倍率從0.1增加到1.3時(shí),各電池串聯(lián)數(shù)下電池組的能量利用率均下降了8%~9%。

    圖7 在不同電池串聯(lián)數(shù)下電池組EUE隨放電電流倍率的演變規(guī)律

    綜合串聯(lián)電池?cái)?shù)和放電倍率的影響分析,可建立電池組能量利用率和關(guān)于電池串聯(lián)數(shù)和放電電流倍率的對(duì)應(yīng)關(guān)系?;谠搶?duì)應(yīng)關(guān)系,可以在電池成組過程中確定最佳的電池串聯(lián)數(shù),并在后續(xù)的使用過程中優(yōu)化放電電流,以達(dá)到電池組能量最大化利用的目的。

    3.3.4 電池參數(shù)一致性的影響

    電池組中電池各個(gè)參數(shù)的一致性同樣影響了電池組的能量利用率。電池參數(shù)分布特性不同,對(duì)電池組能量利用率的影響程度也不同。為此,本文根據(jù)3.1節(jié)的測試結(jié)果,定量地分析容量、內(nèi)阻和SOC三個(gè)參數(shù)的一致性對(duì)電池組能量利用率的影響。在仿真分析中,設(shè)定了四種情景。在情景1中,電池組在不同倍率下的能量利用率為表1所示的實(shí)際測試結(jié)果。而在情景2~4中,則分別消除了電池SOC、容量和內(nèi)阻的不一致(情景2中所有電池SOC均為100%,情景3中所有電池的容量均取電池單體容量測試結(jié)果中的最大值,情景4中所有電池的內(nèi)阻均取電池單體內(nèi)阻測試結(jié)果中的最小值)。四種情景下電池組在不同放電電流倍率下的仿真結(jié)果見表2。從電池組能量利用率的仿真結(jié)果可以看出,針對(duì)本文所測試的電池組,消除電池間SOC的差異可使得電池組能量利用率增加約25%;消除容量的不一致可使電池組能量利用率增加約5%;而消除電池間內(nèi)阻的不一致幾乎不會(huì)提高電池組的能量利用率。電池內(nèi)阻一致性對(duì)電池組能量利用率影響不明顯的主要原因是,本文中所測試電池的內(nèi)阻一致性較好,且仿真計(jì)算中電流倍率較小。由此可見,SOC的差異是本文所測試電池組低能量利用率的主導(dǎo)因素。

    表2 四種情景下電池組能量利用率在不同放電倍率下的仿真結(jié)果

    Tab.2 The battery pack EUE estimation simulation results of the four situations under different discharge current rates (%)

    4 結(jié)論

    針對(duì)梯次利用鋰離子電池組一致性較差進(jìn)而導(dǎo)致能量利用率降低這一問題,本文提出了一種基于Copula理論的電池組一致性建模方法,基于該一致性模型,提出了電池組能量利用率的估計(jì)方法。

    首先,針對(duì)一致性參數(shù)的特點(diǎn),采用Copula函數(shù)來描述變量間的相關(guān)性,介紹了Copula參數(shù)估計(jì)方法和電池組一致性模型抽樣方法。其次,基于電池組一致性模型和蒙特卡洛模擬的思想,提出了電池組能量利用率計(jì)算方法。根據(jù)電池組和電池單體的相關(guān)測試結(jié)果,對(duì)所提出的方法進(jìn)行了精度驗(yàn)證,結(jié)果表明估計(jì)誤差在1%以內(nèi)。最后,定量地分析了蒙特卡洛模擬次數(shù)、電池串聯(lián)數(shù)、放電電流倍率和電池參數(shù)一致性對(duì)電池組能量利用率的影響。結(jié)果表明,當(dāng)蒙特卡洛模擬次數(shù)達(dá)到2 500后,隨機(jī)模擬帶來的能量利用率估計(jì)誤差可降低至0.1%以內(nèi);隨著電池串聯(lián)數(shù)的增多,電池組能量利用率逐漸降低,但下降的速率逐漸減慢,最終能量利用率趨于穩(wěn)定;隨著放電電流倍率的增加,電池組能量利用率呈現(xiàn)出近似線性下降的規(guī)律,并沒有能量利用率迅速下降的拐點(diǎn)出現(xiàn)。針對(duì)文中所測試的電池組,SOC的不一致是電池組低能量利用率的主要影響因素,消除SOC和容量的不一致,電池組能量利用率可分別提高約25%和5%,而消除電池內(nèi)阻的不一致幾乎不會(huì)提高電池組的能量利用率。

    本文提出的基于一致性模型的電池組能量利用率估計(jì)方法和相關(guān)影響因素分析將為梯次利用鋰離子電池組的優(yōu)化成組與使用提供有效指導(dǎo)。

    [1] 焦東升, 王海云, 朱潔, 等. 基于離散Fréchet距離的電動(dòng)汽車電池健康狀態(tài)診斷方法[J]. 電力系統(tǒng)保護(hù)與控制, 2016, 44(12): 68-74.

    Jiao Dongsheng, Wang Haiyun, Zhu Jie, et al. EV battery SOH diagnosis method based on discrete Fréchet distance[J]. Power System Protection and Control, 2016, 44(12): 68-74.

    [2] 劉偉龍, 王麗芳, 王立業(yè). 基于電動(dòng)汽車工況識(shí)別預(yù)測的鋰離子電池SOE估計(jì)[J]. 電工技術(shù)學(xué)報(bào), 2018, 33(1): 17-25.

    Liu Weilong, Wang Lifang, Wang Liye. Estimation of state-of-energy for electric vehicles based on the identification and prediction of driving condition[J]. Transactions of China Electrotechnical Society, 2018, 33(1): 17-25.

    [3] 楊若岑, 冬雷, 廖曉鐘, 等. 電池剩余容量估算方法綜述[J]. 電氣技術(shù), 2019, 20(10): 1-5.

    Yang Ruocen, Dong Lei, Liao Xiaozhong, et al. A review on battery remaining capacity estimation[J]. Electrical Engineering, 2019, 20(10): 1-5.

    [4] 顏湘武, 鄧浩然, 郭琪, 等. 基于自適應(yīng)無跡卡爾曼濾波的動(dòng)力電池健康狀態(tài)檢測及梯次利用研究[J]. 電工技術(shù)學(xué)報(bào), 2019, 34(18): 3937-3948.

    Yan Xiangwu, Deng Haoran, Guo Qi, et al. Study on the state of health detection of power batteries based on adaptive unscented Kalman filters and the battery echelon utilization[J]. Transactions of China Electro- technical Society, 2019, 34(18): 3937-3948.

    [5] 白愷, 李娜, 范茂松, 等. 大容量梯次利用電池儲(chǔ)能系統(tǒng)工程技術(shù)路線研究[J]. 華北電力技術(shù), 2017(3): 39-45.

    Bai Kai, Li Na, Fan Maosong, et al. Research on the technical roadmap for engineering application of large-scale echelon use battery energy storage system[J]. North China Electric Poser, 2017(3): 39-45.

    [6] Klein M, Tong Shijie, Park J W. In-plane nonuniform temperature effects on the performance of a large- format lithium-ion pouch cell[J]. Applied Energy, 2016, 165: 639-647.

    [7] 郭永芳, 黃凱, 李志剛. 基于短時(shí)擱置端電壓壓降的快速鋰離子電池健康狀態(tài)預(yù)測[J]. 電工技術(shù)學(xué)報(bào), 2019, 34(19): 3968-3978.

    Guo Yongfang, Huang Kai, Li Zhigang. Fast state of health prediction of lithium-ion battery based on terminal voltage drop during rest for short time[J]. Transactions of China Electrotechnical Society, 2019, 34(19): 3968-3978.

    [8] Schuster S F, Brand M J, Berg P, et al. Lithium-ion cell-to-cell variation during battery electric vehicle operation[J]. Journal of Power Sources, 2015, 297: 242-251.

    [9] 王震坡, 孫逢春, 張承寧. 電動(dòng)汽車動(dòng)力蓄電池組不一致性統(tǒng)計(jì)分析[J]. 電源技術(shù), 2003(5): 438-441.

    Wang Zhenpo, Sun Fengchun, Zhang Chengning. Study on inconsistency of electric vehicle battery pack[J]. Chinese Journal of Power Sources, 2003(5): 438-441.

    [10] 魏五星. 磷酸鐵鋰動(dòng)力電池組性能測試與分析[D]. 武漢: 武漢理工大學(xué), 2010.

    [11] 程功. 電池組一致性的統(tǒng)計(jì)特性與變化規(guī)律研究[D]. 北京: 北京交通大學(xué), 2017.

    [12] Zheng Yuejiu, Ouyang Minggao, Lu Languang, et al. Understanding aging mechanisms in lithium-ion battery packs: from cell capacity loss to pack capacity evolution[J]. Journal of Power Sources, 2015, 278: 287-295.

    [13] Yang Fangfang, Wang Dong, Zhao Yang, et al. A study of the relationship between coulombic effici- ency and capacity degradation of commercial lithium- ion batteries[J]. Energy, 2018, 145: 486-495.

    [14] Diao Weiping, Xue Nan, Bhattacharjee V, et al. Active battery cell equalization based on residual available energy maximization[J]. Applied Energy, 2018, 210: 690-698.

    [15] Zhou Long, Zheng Yuejiu, Ouyang Minggao, et al. A study on parameter variation effects on battery packs for electric vehicles[J]. Journal of Power Sources, 2017, 364: 242-252.

    [16] Xi Zhimin, Jing Rong, Wang Pingfeng, et al. A copula-based sampling method for data-driven pro- gnostics[J]. Reliability Engineering & System Safety, 2014, 132: 72-82.

    Energy Utilization Efficiency Estimation Method for Second-Use Lithium-Ion Battery Packs Based on a Battery Consistency Model

    12223,4,5

    (1. Electric Power Research Institute of Guangdong Power Grid Co. Ltd Guangzhou 510080 China 2. Sunwoda Energy Solution Co. Ltd Shenzhen 518000 China 3. Sunwoda Electronic Co. Ltd Shenzhen 518108 China 4. School of Electric Power South China University of Technology Guangzhou 510641 China 5. National Active Distribution Network Technology Research Center Beijing Jiaotong University Beijing 100044 China)

    The inconsistency of lithium-ion battery cells has caused the lifetime and capacity of the battery pack to be much smaller than that of the battery cells, and this problem is more serious in the second-use applications. Along with the aging of the battery pack, the battery consistency parameters are mutual dependent, which causes battery pack consistency modeling difficult and state of health estimation inaccurate. Aiming at the grouping and the state evaluation stages of second-use batteries, this paper proposes a Copula-based battery pack consistency modeling method. The proposed method can describe the statistical characteristics of the battery consistency parameters and the dependence structure between parameters. Based on the consistency model, a battery pack energy utilization efficiency (EUE) estimation method is put forward and verified by experiments. The estimation error is within 1%. Finally, the influencing factors of EUE is analyzed to guide the optimal grouping and operation of second-use lithium-ion battery packs.

    Lithium-ion battery pack, second use, consistency, energy utilization efficiency

    TM911

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

    國家重點(diǎn)研發(fā)計(jì)劃資助項(xiàng)目(2018YFB0905304)。

    2020-03-10

    2020-07-10

    趙 偉 男,1981年生,博士,教授級(jí)高工,研究方向?yàn)橄冗M(jìn)電力儲(chǔ)能集成及控制技術(shù)。E-mail: zh_solar@126.com

    姜 研 男,1991年生,博士,研究方向?yàn)殇囯x子電池管理與成組應(yīng)用技術(shù)。E-mail: jiangyan3@sunwoda.com (通信作者)

    (編輯 崔文靜)

    猜你喜歡
    內(nèi)阻電池組倍率
    大型桅桿起重機(jī)起升變倍率方法及其應(yīng)用
    Effect of Xuebijing injection on hematopoietic homeostasis of LPS induced sepsis in mice
    FANUC0iD系統(tǒng)速度倍率PMC控制方法
    2017年7月原電池及原電池組產(chǎn)量同比增長2.53%
    “測定電池的電動(dòng)勢和內(nèi)阻”復(fù)習(xí)課之八問
    一種智能加工系統(tǒng)中的機(jī)床倍率控制方法
    伏安法測電源電動(dòng)勢和內(nèi)阻的測量值與真實(shí)值
    超級(jí)電容器內(nèi)阻測試方法研究
    拉伸倍率對(duì)BOPP薄膜性能的影響
    中國塑料(2015年6期)2015-11-13 03:03:05
    基于LTC6802的電池組均衡電路設(shè)計(jì)
    亚洲内射少妇av| 国产成人a区在线观看| 免费大片黄手机在线观看| 亚洲av中文av极速乱| 亚州av有码| 国产黄片视频在线免费观看| 国产精品一二三区在线看| 久久久久久久亚洲中文字幕| 欧美xxⅹ黑人| 一级毛片 在线播放| 国产免费视频播放在线视频| 欧美 日韩 精品 国产| 国产欧美日韩精品一区二区| 国产午夜福利久久久久久| 亚洲国产精品专区欧美| 天美传媒精品一区二区| 中文字幕制服av| 成人午夜精彩视频在线观看| 亚洲激情五月婷婷啪啪| 大片电影免费在线观看免费| 下体分泌物呈黄色| 日韩欧美一区视频在线观看 | 精品少妇黑人巨大在线播放| 欧美日韩在线观看h| 欧美性猛交╳xxx乱大交人| 国产 精品1| 国产成人a区在线观看| 美女脱内裤让男人舔精品视频| 美女主播在线视频| 性插视频无遮挡在线免费观看| 嘟嘟电影网在线观看| 欧美高清成人免费视频www| 1000部很黄的大片| 亚洲国产欧美人成| 少妇 在线观看| av专区在线播放| 制服丝袜香蕉在线| 欧美另类一区| 秋霞在线观看毛片| 亚洲,欧美,日韩| 欧美亚洲 丝袜 人妻 在线| 国产一区二区亚洲精品在线观看| 久久精品国产亚洲网站| 在线精品无人区一区二区三 | 深爱激情五月婷婷| 性色av一级| 亚洲精品一二三| 久久久午夜欧美精品| 久久久久久久午夜电影| 欧美丝袜亚洲另类| 在线精品无人区一区二区三 | 嫩草影院新地址| 美女内射精品一级片tv| 丝瓜视频免费看黄片| 伦理电影大哥的女人| 久久精品国产自在天天线| av国产精品久久久久影院| 干丝袜人妻中文字幕| 亚洲无线观看免费| 看非洲黑人一级黄片| 免费观看性生交大片5| 香蕉精品网在线| 一个人看视频在线观看www免费| 亚洲精品,欧美精品| 久久ye,这里只有精品| 狂野欧美白嫩少妇大欣赏| av在线老鸭窝| 久久久久久九九精品二区国产| 欧美精品国产亚洲| 亚洲最大成人手机在线| 成人欧美大片| 婷婷色麻豆天堂久久| 国产精品爽爽va在线观看网站| 狂野欧美白嫩少妇大欣赏| 看黄色毛片网站| 亚洲一区二区三区欧美精品 | 亚洲成色77777| 亚洲欧洲国产日韩| 免费观看的影片在线观看| 国产乱人偷精品视频| 老司机影院毛片| 成人国产av品久久久| 国产精品久久久久久久电影| 国产一区二区亚洲精品在线观看| 91精品国产九色| 少妇高潮的动态图| 亚洲自偷自拍三级| 国产精品一区www在线观看| 秋霞伦理黄片| 亚洲天堂av无毛| 女人十人毛片免费观看3o分钟| 欧美高清性xxxxhd video| 成人亚洲精品一区在线观看 | 一级a做视频免费观看| 亚洲国产色片| 麻豆成人av视频| 欧美日韩视频精品一区| av国产久精品久网站免费入址| 最近最新中文字幕大全电影3| 国产av码专区亚洲av| 久久久成人免费电影| 日韩欧美精品免费久久| 一级爰片在线观看| 亚洲av欧美aⅴ国产| 国产精品三级大全| 91精品国产九色| 人妻制服诱惑在线中文字幕| 国产欧美另类精品又又久久亚洲欧美| 99热这里只有是精品在线观看| 国产精品99久久久久久久久| 国产精品一及| 少妇 在线观看| 国产成人一区二区在线| 婷婷色综合大香蕉| 春色校园在线视频观看| av在线老鸭窝| 国产午夜精品久久久久久一区二区三区| 精品少妇久久久久久888优播| 蜜桃久久精品国产亚洲av| 国产成人精品福利久久| 亚洲内射少妇av| 久久精品国产亚洲网站| 欧美日韩综合久久久久久| 大片电影免费在线观看免费| 一级av片app| 少妇高潮的动态图| 大码成人一级视频| 蜜桃亚洲精品一区二区三区| 看十八女毛片水多多多| 少妇的逼好多水| 久久久久国产精品人妻一区二区| 女人久久www免费人成看片| 国产精品久久久久久精品电影| 高清欧美精品videossex| 秋霞伦理黄片| 国产精品av视频在线免费观看| 18禁裸乳无遮挡免费网站照片| 亚洲在线观看片| 日韩欧美 国产精品| 激情 狠狠 欧美| 在线亚洲精品国产二区图片欧美 | 亚洲国产色片| 亚洲在久久综合| 在线观看一区二区三区| 亚洲精品日本国产第一区| 免费av观看视频| 哪个播放器可以免费观看大片| 一个人观看的视频www高清免费观看| 亚洲久久久久久中文字幕| 精品国产露脸久久av麻豆| 国产精品人妻久久久影院| 99久国产av精品国产电影| 免费少妇av软件| 久久热精品热| 午夜亚洲福利在线播放| 亚洲国产精品国产精品| 中文字幕av成人在线电影| 狂野欧美激情性xxxx在线观看| 亚洲精品成人av观看孕妇| 亚洲,欧美,日韩| 99久久精品热视频| 丝袜脚勾引网站| 国产午夜精品一二区理论片| 久久99蜜桃精品久久| 亚洲欧美日韩卡通动漫| 男女那种视频在线观看| 午夜日本视频在线| 一二三四中文在线观看免费高清| 夜夜爽夜夜爽视频| 中文欧美无线码| 国产黄片视频在线免费观看| 亚洲av在线观看美女高潮| 亚洲电影在线观看av| 各种免费的搞黄视频| 精品一区二区三区视频在线| 国产一区二区三区综合在线观看 | 欧美精品一区二区大全| 一级毛片久久久久久久久女| videos熟女内射| 22中文网久久字幕| 91在线精品国自产拍蜜月| 欧美成人午夜免费资源| 国产乱人视频| 午夜免费鲁丝| 国产日韩欧美亚洲二区| 亚洲av中文av极速乱| 免费大片18禁| 国产色爽女视频免费观看| av网站免费在线观看视频| 久久ye,这里只有精品| 如何舔出高潮| 在线看a的网站| 一级毛片电影观看| 狠狠精品人妻久久久久久综合| av天堂中文字幕网| 在线观看国产h片| 香蕉精品网在线| 人体艺术视频欧美日本| av在线天堂中文字幕| 最近的中文字幕免费完整| 亚洲,欧美,日韩| 亚洲怡红院男人天堂| 丝袜脚勾引网站| 亚洲在久久综合| 日日撸夜夜添| 久久精品久久久久久久性| 少妇高潮的动态图| 男女边吃奶边做爰视频| 国产爽快片一区二区三区| 欧美精品人与动牲交sv欧美| 亚洲va在线va天堂va国产| 国产精品国产av在线观看| 国产伦精品一区二区三区视频9| 麻豆成人午夜福利视频| av福利片在线观看| 最近中文字幕高清免费大全6| tube8黄色片| 亚洲精品,欧美精品| 一本久久精品| 久久久久久久久久久免费av| 在线观看美女被高潮喷水网站| 新久久久久国产一级毛片| 亚洲精品自拍成人| 国内精品宾馆在线| 久久鲁丝午夜福利片| 2018国产大陆天天弄谢| 国产在线一区二区三区精| 男女边摸边吃奶| 欧美成人a在线观看| 国产精品国产av在线观看| 亚洲欧美日韩另类电影网站 | 99久久精品国产国产毛片| 新久久久久国产一级毛片| 在线天堂最新版资源| 欧美精品国产亚洲| 亚洲最大成人手机在线| 久久久a久久爽久久v久久| www.av在线官网国产| 成人毛片a级毛片在线播放| 在线观看三级黄色| 国产精品嫩草影院av在线观看| 真实男女啪啪啪动态图| 身体一侧抽搐| 在线观看三级黄色| 国产精品99久久99久久久不卡 | 午夜精品一区二区三区免费看| 日本-黄色视频高清免费观看| 久久久国产一区二区| 亚洲精品第二区| 少妇的逼水好多| 亚洲激情五月婷婷啪啪| 国产老妇伦熟女老妇高清| 国产精品一区二区三区四区免费观看| 国产精品精品国产色婷婷| 在线观看人妻少妇| 国产伦精品一区二区三区视频9| a级毛色黄片| 寂寞人妻少妇视频99o| 久久97久久精品| 国产欧美日韩精品一区二区| 亚洲国产欧美在线一区| 亚洲国产精品国产精品| 免费播放大片免费观看视频在线观看| 最近的中文字幕免费完整| 国产av码专区亚洲av| 国产黄片视频在线免费观看| 在线播放无遮挡| 在线观看三级黄色| 免费黄频网站在线观看国产| 亚洲伊人久久精品综合| 中文字幕亚洲精品专区| 免费高清在线观看视频在线观看| 久久久久精品性色| 麻豆久久精品国产亚洲av| 国产成人免费无遮挡视频| 少妇的逼水好多| 各种免费的搞黄视频| 狠狠精品人妻久久久久久综合| 2022亚洲国产成人精品| 黄色视频在线播放观看不卡| 亚洲成人av在线免费| 国产黄a三级三级三级人| 久久人人爽人人片av| 18+在线观看网站| 精品人妻视频免费看| av女优亚洲男人天堂| 三级男女做爰猛烈吃奶摸视频| 欧美日韩视频高清一区二区三区二| 亚洲三级黄色毛片| 美女国产视频在线观看| 久久精品国产亚洲av天美| 国产一区二区三区av在线| 日韩 亚洲 欧美在线| 国产av码专区亚洲av| 一二三四中文在线观看免费高清| 久久鲁丝午夜福利片| 久久久精品欧美日韩精品| 欧美一级a爱片免费观看看| 99久久精品热视频| 别揉我奶头 嗯啊视频| 深夜a级毛片| 水蜜桃什么品种好| 国产永久视频网站| 麻豆成人午夜福利视频| 永久免费av网站大全| 国产成人a∨麻豆精品| 亚洲在久久综合| 国产视频首页在线观看| 欧美97在线视频| 亚洲av.av天堂| 极品教师在线视频| 日本色播在线视频| 91精品国产九色| 一区二区三区精品91| 在线 av 中文字幕| 在线看a的网站| 亚洲va在线va天堂va国产| 久久久久久久久久久免费av| 少妇猛男粗大的猛烈进出视频 | 国语对白做爰xxxⅹ性视频网站| 天天一区二区日本电影三级| 啦啦啦中文免费视频观看日本| 午夜免费男女啪啪视频观看| 激情五月婷婷亚洲| 另类亚洲欧美激情| 久久久精品欧美日韩精品| a级一级毛片免费在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 天天躁夜夜躁狠狠久久av| 国产精品熟女久久久久浪| 在线天堂最新版资源| 18禁在线无遮挡免费观看视频| 不卡视频在线观看欧美| 天美传媒精品一区二区| 亚洲最大成人中文| 一二三四中文在线观看免费高清| 五月开心婷婷网| 日本爱情动作片www.在线观看| 尾随美女入室| av网站免费在线观看视频| 嫩草影院入口| 亚洲成色77777| 少妇猛男粗大的猛烈进出视频 | 99久国产av精品国产电影| 最近中文字幕2019免费版| 精品久久久久久久末码| 欧美区成人在线视频| 18禁在线无遮挡免费观看视频| 精品酒店卫生间| 亚洲色图综合在线观看| 色综合色国产| 久久午夜福利片| 国产精品一区二区三区四区免费观看| 91精品伊人久久大香线蕉| 人人妻人人看人人澡| 免费播放大片免费观看视频在线观看| 久久久久精品性色| 国产伦精品一区二区三区四那| 日本色播在线视频| 久久久久精品性色| 国产欧美亚洲国产| 国产黄a三级三级三级人| 狠狠精品人妻久久久久久综合| 91久久精品电影网| 日本爱情动作片www.在线观看| 黄色欧美视频在线观看| 一区二区三区精品91| 一区二区三区乱码不卡18| 久久久a久久爽久久v久久| 久久99精品国语久久久| videossex国产| 尾随美女入室| av免费在线看不卡| 欧美3d第一页| 一级毛片黄色毛片免费观看视频| 18禁在线播放成人免费| 国产色爽女视频免费观看| 国产真实伦视频高清在线观看| 大香蕉97超碰在线| 内射极品少妇av片p| 在线观看三级黄色| 在线亚洲精品国产二区图片欧美 | 亚洲自偷自拍三级| 国产爽快片一区二区三区| 久久久午夜欧美精品| 亚洲精品成人av观看孕妇| 亚洲精品亚洲一区二区| 成年免费大片在线观看| 黄片wwwwww| 久久久久久久国产电影| 国产成人91sexporn| 精品人妻一区二区三区麻豆| xxx大片免费视频| 国产探花在线观看一区二区| 一级毛片电影观看| 欧美人与善性xxx| 综合色丁香网| 国产精品一区二区三区四区免费观看| 成人无遮挡网站| 精品国产露脸久久av麻豆| 禁无遮挡网站| 久久久久久九九精品二区国产| 最近的中文字幕免费完整| freevideosex欧美| 国产一级毛片在线| 日韩精品有码人妻一区| 97在线人人人人妻| 国产亚洲5aaaaa淫片| 男人和女人高潮做爰伦理| 99久久九九国产精品国产免费| 我要看日韩黄色一级片| 欧美+日韩+精品| 国精品久久久久久国模美| h日本视频在线播放| av国产久精品久网站免费入址| 国产有黄有色有爽视频| 国产淫片久久久久久久久| 两个人的视频大全免费| 人妻系列 视频| 国产精品成人在线| 欧美成人午夜免费资源| 国产精品人妻久久久影院| av专区在线播放| 听说在线观看完整版免费高清| 人妻一区二区av| 久久97久久精品| 久久久久久久久久人人人人人人| 国产亚洲精品久久久com| 九九在线视频观看精品| 国产免费视频播放在线视频| 特级一级黄色大片| 2018国产大陆天天弄谢| 99久久人妻综合| 日本午夜av视频| 日韩一区二区三区影片| 永久网站在线| 午夜激情久久久久久久| 高清在线视频一区二区三区| 精华霜和精华液先用哪个| 国产中年淑女户外野战色| 国产精品一区二区性色av| 能在线免费看毛片的网站| 欧美激情久久久久久爽电影| 舔av片在线| 久久精品熟女亚洲av麻豆精品| 美女国产视频在线观看| 蜜桃久久精品国产亚洲av| av免费观看日本| 丝袜脚勾引网站| 欧美亚洲 丝袜 人妻 在线| 欧美 日韩 精品 国产| 亚洲国产精品专区欧美| 亚洲丝袜综合中文字幕| 偷拍熟女少妇极品色| 国产精品.久久久| 最近手机中文字幕大全| 22中文网久久字幕| 如何舔出高潮| 亚州av有码| 看黄色毛片网站| av在线播放精品| 美女脱内裤让男人舔精品视频| 在线免费十八禁| 国产精品偷伦视频观看了| 免费黄网站久久成人精品| 国产 一区 欧美 日韩| 男女那种视频在线观看| av播播在线观看一区| 日本一二三区视频观看| 91狼人影院| 嫩草影院精品99| 日韩欧美一区视频在线观看 | 亚洲国产最新在线播放| 亚洲精华国产精华液的使用体验| 一级毛片电影观看| 国产伦精品一区二区三区四那| a级毛片免费高清观看在线播放| 青春草视频在线免费观看| 欧美bdsm另类| 国产久久久一区二区三区| 五月伊人婷婷丁香| 国产精品久久久久久久久免| 精品人妻偷拍中文字幕| av女优亚洲男人天堂| av在线蜜桃| 国产免费一区二区三区四区乱码| 天天一区二区日本电影三级| 欧美日韩视频高清一区二区三区二| 欧美老熟妇乱子伦牲交| 久久精品综合一区二区三区| 精品一区二区三卡| 亚州av有码| 亚洲在线观看片| 在线观看三级黄色| 可以在线观看毛片的网站| 搡老乐熟女国产| 十八禁网站网址无遮挡 | 亚洲精品自拍成人| 亚洲精品第二区| 日产精品乱码卡一卡2卡三| 狂野欧美激情性bbbbbb| 全区人妻精品视频| 一区二区三区免费毛片| 日韩国内少妇激情av| 久久久久久久精品精品| 人人妻人人爽人人添夜夜欢视频 | 色视频在线一区二区三区| 亚洲国产精品999| 久久人人爽av亚洲精品天堂 | 欧美激情久久久久久爽电影| 日本爱情动作片www.在线观看| 亚洲在久久综合| 综合色丁香网| 特级一级黄色大片| 少妇裸体淫交视频免费看高清| 亚洲精华国产精华液的使用体验| 欧美国产精品一级二级三级 | 插阴视频在线观看视频| 国产精品.久久久| 亚洲精品影视一区二区三区av| 久久久久久久久久人人人人人人| 18禁在线无遮挡免费观看视频| 老司机影院毛片| 日韩免费高清中文字幕av| 全区人妻精品视频| 亚洲四区av| 免费av不卡在线播放| 性色avwww在线观看| 亚洲欧美精品自产自拍| 在线免费十八禁| 国产永久视频网站| 一级黄片播放器| 在线观看av片永久免费下载| 国产精品嫩草影院av在线观看| 两个人的视频大全免费| 男女边吃奶边做爰视频| 亚洲丝袜综合中文字幕| 国产精品爽爽va在线观看网站| 中文资源天堂在线| 久久人人爽av亚洲精品天堂 | 成人免费观看视频高清| 亚洲va在线va天堂va国产| 只有这里有精品99| 又黄又爽又刺激的免费视频.| 18禁在线播放成人免费| 国产午夜福利久久久久久| 免费看日本二区| 国产成人免费观看mmmm| 久久久久久久亚洲中文字幕| 日本猛色少妇xxxxx猛交久久| 成人午夜精彩视频在线观看| 在线观看国产h片| 十八禁网站网址无遮挡 | 国产白丝娇喘喷水9色精品| 免费观看的影片在线观看| 国产一区二区三区av在线| 又黄又爽又刺激的免费视频.| 午夜亚洲福利在线播放| 大话2 男鬼变身卡| 狠狠精品人妻久久久久久综合| 久久久久久国产a免费观看| 国产精品一区二区在线观看99| 简卡轻食公司| 亚洲欧美一区二区三区黑人 | 美女脱内裤让男人舔精品视频| 日本一二三区视频观看| 亚洲欧美一区二区三区黑人 | 欧美最新免费一区二区三区| 春色校园在线视频观看| 国内少妇人妻偷人精品xxx网站| 久久精品久久久久久久性| 少妇丰满av| 国产高清三级在线| 国产国拍精品亚洲av在线观看| 国产精品国产av在线观看| 嫩草影院新地址| 三级男女做爰猛烈吃奶摸视频| 高清欧美精品videossex| 午夜福利高清视频| 麻豆乱淫一区二区| 美女主播在线视频| 涩涩av久久男人的天堂| 大码成人一级视频| 赤兔流量卡办理| 亚洲av一区综合| 国产探花极品一区二区| 亚洲aⅴ乱码一区二区在线播放| 色网站视频免费| 成年人午夜在线观看视频| av网站免费在线观看视频| 天堂俺去俺来也www色官网| 欧美bdsm另类| 中文在线观看免费www的网站| 有码 亚洲区| 亚洲无线观看免费| 国产免费视频播放在线视频| 欧美三级亚洲精品| 99热6这里只有精品| 精品国产乱码久久久久久小说| 国产综合懂色| 综合色av麻豆| 国产探花极品一区二区| 一本久久精品| 婷婷色麻豆天堂久久| 亚洲精品国产av蜜桃| 亚洲一级一片aⅴ在线观看| 国产在线男女| 少妇裸体淫交视频免费看高清| 国产免费一级a男人的天堂| 亚洲av在线观看美女高潮| 偷拍熟女少妇极品色| 男的添女的下面高潮视频| 乱码一卡2卡4卡精品| 黄色怎么调成土黄色| 成人毛片a级毛片在线播放|