欧美黑人又粗又大的性格特点,国产农村妇女aaaaa视频,欧美高清精品一区二区,好爽又高潮了毛片免费下载

基于特點處理與徑向基神經網絡的鋰離子電池剩余容量估算方法

鉅大鋰電  |  點擊量:0  |  2021年05月21日  

我國儲能網訊:


摘要:為解決鋰離子電池可用容量估算過程中精度與效率難以兼顧的問題,本文提出了一種基于特點處理與徑向基神經網絡的鋰離子電池剩余容量估計方法。首先由電池充電過程數據中提取與剩余可用容量相關聯的特點量,然后運用局部異常因子算法對特點量中異常點進行精準清洗,提高特點量所含有效信息量,再通過局部線性嵌入降維算法對所得特點向量組進行降維處理,減少數據復雜度,最后,引入徑向基神經網絡建立起剩余容量的估算模型。在不同型號電池上應用該模型進行了驗證,估算結果的最大平均絕對誤差為0.06,最大均方根誤差為0.05,表明該模型能夠有效估計鋰離子電池的剩余可用容量并有較強的魯棒性。與Elman神經網絡和BP神經網絡算法相比,在保證高精度的同時該方法有更快的估算效率。


關鍵詞:鋰離子電池;特點處理;徑向基神經網絡;容量估計


鋰離子電池因具有能量密度高、循環壽命長、成本價格低、環保污染小等優點而成為目前電動汽車動力鋰離子電池的最佳選擇。在實際應用中,由于電池內部鋰離子的活性會逐漸下降,導致鋰離子電池的可用容量會在使用過程中有所損失,直接影響其剩余使用壽命(remainingusefullife,RUL)。因此,為高效利用鋰離子電池,有效估計當前鋰離子電池的剩余可用容量,進而時時診斷其健康狀態(stateofhealth,SOH)是電池管理系統(batterymanagementsystem,BMS)的重要功能。


當前,研究人員已提出多種估算方法來獲得可用容量,重要分為基于貝葉斯的方法、相關相關經驗擬合的方法和數據驅動算法等。基于貝葉斯的方法包含卡爾曼濾波以及粒子濾波及其他相關改進算法。文獻采用雙擴展卡爾曼濾波算法同時獲得模型參數和當前電池狀態,但此方法在電池循環壽命后期表現不佳。針對這一問題,文獻采用自適應無跡卡爾曼濾波算法,通過更新噪聲協方差矩陣實現了SOH和RUL的多步預測,提高了電池全壽命周期內的估計精度。文獻采用無跡粒子濾波來估計鋰離子電池的RUL,實際估算相對誤差小于5%。以上基于貝葉斯的方法利用閉環濾波算法來估計電池健康狀態的關鍵參數,對建模誤差和測量誤差具有較強的魯棒性,但是建模過程復雜,實際估算過程計算量較大。


在采用各種濾波算法之外,相關相關經驗擬合的方法也被用于鋰離子電池的狀態估計領域。文獻考慮鋰離子電池的使用時間、充放電倍率和溫度等條件,建立起廣義電池壽命模型,實現對同一類電池健康狀態的有效估計。文獻結合鋰離子電池衰退后固態電解質界面層的生長和活性物質的損失情況建立電池老化模型,并采用遞推最小二乘法估計模型參數,結果表明該方法所得結果相對誤差小于1%。相關相關經驗擬合的方法一般是在工程相關相關經驗或電化學知識的基礎上建立老化模型,精度較高,但是模型擬合要大量的測試數據。并且,特定模型往往只能有效估計特定類型的電池或者特定老化情況下的電池狀態,靈活性較差。


如今,隨著大數據和人工智能的快速發展,數據驅動的方法也被廣泛應用于電池的狀態估計。文獻提出一種基于先驗神經網絡和馬爾可夫鏈的SOH估計方法,結果表明即使在不確定的外部條件下,該方法仍能夠有效估算電池的內部狀態。文獻采用具有自動相關確定結構的高斯過程回歸算法(gaussianprocessregression,GPR)得到電池容量、溫度和荷電狀態(stateofcharge,SOC)之間的映射,實現對鋰離子電池的多狀態同時估計。文獻采用長短時間記憶神經網絡進一步優化GPR模型的結構,得到估算精度較高的估算模型。以上基于數據驅動的方法是一種黑箱模型,具有靈活智能的優點,但黑箱模型對訓練數據集的好壞較敏感,較少的數據使模型得不到有效的訓練,較多的數據集則有可能會導致模型復雜冗余,且可能造成模型的過擬合。


為了解決采用數據驅動方法難以獲得高效特點集的問題,本文基于特點處理與徑向基函數(radialbasisfunction,RBF)神經網絡進行了鋰離子電池剩余容量的估算,在保證精度的同時提高了估算效率。首先,采用局部離群因子方法對充電過程中提取的容量特點量進行優化處理,獲得特點信息更集中的數據集,并對所獲得的數據集進行特點選擇,進一步降低數據的復雜度以提高估算效率,然后利用簡潔高效的徑向基神經網絡建立剩余可用容量估計模型,并最終實現對鋰離子電池剩余可用容量的快速準確估算。


1鋰離子電池循環老化實驗


為研究動力鋰離子電池衰退后可用容量隨電池使用循環的變化情況,本文對兩種規格的三元鋰離子動力鋰離子電池在室溫條件下進行持續的恒流充放循環老化實驗,操作流程見圖1,電池相關信息如表1所示。實驗采用美國Arbin公司開發的BT-5HC-5V-100A電池測試系統進行。


圖1電池循環老化實驗流程


表1測試電池基本參數


容量衰退實驗完成后共得到四組電池單體全壽命循環測試數據,分別標號為電池A-1、電池B-2、電池B-3、電池B-4,其中電池A-1為A型號鋰離子電池,其它3組為B型號鋰離子電池。以電池A-1為例,某一循環的容量衰退實驗電壓與電流變化如圖2所示,每個測試循環均由恒流充電、恒壓充電、擱置、恒流放電、擱置等5部分組成,其中A型號恒流充電過程電流大小為2A,B型號為1.275A,A型號恒壓充電過程截止電流為80mA,B型號為51mA,A型號恒流放電過程電流為8A,B型號為5.1A。


圖2循環測試過程電壓與電流變化


2放電容量估算實現方法


2.1容量衰退特點處理


2.1.1特點量獲取


在實際的應用過程中,考慮到鋰離子電池的充電過程均為恒流恒壓充電模式,對不同型號的電池具有一致性,而放電過程則因實際使用情況而存在差異,難以提取到穩定的特點量。因此,本文在充電階段的實驗數據中提取7組與鋰離子電池剩余可用容量相關的特點量,并根據數據來源劃分為兩類。


第一類為根據傳感器直接采集數據獲得的特點量。分別將恒流充電階段的充電時間CC-T、充電容量CC-C、恒壓充電階段的充電時間CV-T、充電容量CV-C作為特點量。以電池A-1為例,所提取的第一類特點量如圖3所示。


圖3電池A-1的特點量


第二類特點量為對傳感器直接采集數據進一步挖掘所得到的特點量。其中,恒流充電階段容量增量(incrementalcapacity,IC)曲線的峰值是有效表征鋰離子電池內部理化特性的狀態量,IC曲線含義為容量與電壓的微分,見式(1)


式中,Vt代表t時刻的電池端電壓。在恒流充電模式下,電流為恒定值,當取極小的時間間隔圖片時,式(1)是成立的,由此可以快速繪制出IC曲線。繪制的IC曲線如圖4所示,濾波后,可見曲線呈現明顯的峰值特點。將每一循環IC曲線的最大值作為第二類特點量,稱為:Max-IC。此外,分別取恒流階段充電時間與對應循環恒壓階段充電時間的比值:CC-T/CV-T和恒流階段充電容量與對應循環恒壓階段充電容量的比值:CC-C/CV-C作為第二類相關特點量。電池A-1的第二類特點量如圖5所示。


圖4IC曲線濾波前后


圖5電池A-1的特點量


在實驗數據的采集過程中,存在真實的電池特性數據與誤差干擾數據摻雜混合的情況,在如圖3、圖5所示的兩類特點量中,均存在有明顯的數據離群點,而測試誤差和電池自身衰退特性均會造成離群點的出現。因此,本文采用局部異常因子算法消除測試誤差對特點量規律性的影響,保留鋰離子電池衰退過程中真實理化特性所反映出的特點,有效提高特點量所含信息量,進而提高模型學習效率。


局部異常因子算法(localoutlierfactor,LOF)是通過計算每個樣本的局部離群因子來判斷其離群程度,是一種基于密度的離群點檢測方法。能夠依據設定的閾值精準定位無效的離群點,防止了依據數據分布判別離群值的方法存在的誤清洗有效數據的問題。特點量中采樣點p鄰域點的局部可達密度與點p的局部可達密度之比的平均數就表示為點p的局部離群因子,見式(2)


式中,NK(p)表示點p的k距離鄰域,點o為鄰域內的某一點;lrdk(p),lrdk(o)分別表示點p與點o的局部可達密度。局部離群因子越接近1,說明點p與其鄰域點的周圍密度相差無幾,p和鄰域同屬一簇的概率越高,反之,p點可能是異常點。圖6顯示了應用局部離群因子甄別A型號電池特點量CC-T中異常離群點結果,原始數據中點1~5均是明顯可見的離群點,其中,點1和點2是由于測試環境所造成的異常值,要進行平滑處理,而點3~5處數值的變化是電池的可用容量回升現象在特點量CC-T中的反映,不能一并進行平滑處理,在特點選擇的過程中應將其予以保留。


圖6異常離群點的甄別結果


2.1.2相關性檢驗


為驗證特點量與可用容量之間是否具有密切的相關性,要選用合適的方法對兩者進行相關性檢驗。考慮到特點量與可用容量之間的關系是高度非線性的,因此,需選用對特點數據分布沒有要求的分析方法。


斯皮爾曼等級(spearmanrank)相關系數又稱為“等級差數法”,其取值在-1到+1之間,絕對值接近1表明相關性越強。與其他積差相關系數相比,斯皮爾曼等級相關系數有關數據錯誤和極端值的反應不敏感,僅要求兩個變量的觀測值是成對的等級評定資料,或者是由持續變量觀測資料轉化得到的等級資料。不論兩個變量的總體分布形態、樣本容量的大小如何,斯皮爾曼等級相關系數均能夠有效表征兩者的相關程度。


斯皮爾曼相關系數表示為式(3)


式中,Ri和Si分別是觀測值i的取值等級;R*和S*分別是變量x和y的平均等級;di=Ri-Si表示兩列成對變量的等差級數;N是觀測值的總數量。


實驗所得四組電池的循環放電容量與相應特點量之間的斯皮爾曼等級相關系數如表2所示,表中列出了經過異常離群點處理前后的相關性系數。結果顯示,所提取的特點量與對應的剩余可用容量之間相關系數均在0.9左右,呈現出強相關性。且經LOF算法處理異常離群點之后,相關性系數均有提高,顯著的增強了特點向量與對應循環容量之間的相關性。


表2特點量與可用容量的Spearman等級相關系數


2.1.3特點降維


實際應用過程中,處理龐大的數據量會占用大量車載計算模塊的計算能力,從而影響BMS的整體性能。本文采用局部線性嵌入算法(locallylinearembedding,LLE)對特點數據集進行降維處理,減少系統所需的存儲空間,加快計算速度。LLE屬于流形學習的一種,在低維空間保持了原始高維空間樣本鄰域內的線性關系,能夠較好的保持數據整體的幾何結構和性質,適合應用于高度非線性數據的降維處理。LLE算法對高維數據進行低維映射分為以下幾步。


第一步:選取樣本點的局部鄰域點集P


有關給定的數據集X={x1,x2,…,xN}∈?D×N,xi∈?D×1,i=1,2,…,N,計算樣本點的歐幾里得距離,選取樣本點xi鄰域的k(k


第二步:計算線性重構的系數矩陣ω


如式(4)所示,含義均方差損失函數,求解樣本點鄰域的線性重構權重,構建局部重建權值矩陣。


式中,xij表示點xi的第j個近鄰點,ωij表示樣本點xi鄰域內第j個近鄰點的重構權重,當xij不屬于xi的鄰域時,ωij=0。約束條件圖片表示權值矩陣ω中的每一行相加為1。


第三步:通過得到的權值矩陣ω得到樣本集X的低維嵌入Y


將所有樣本點映射到低維空間?d,即將D維的向量X降維至d維,每一個高維向量xi對應一個低維向量yi,Y={y1,y2,…,yN}∈?d×N的映射條件為最小化映射損失函數,見式(5)


對Y的求解等價于求解一定約束條件下的稀疏矩陣的特點向量,采用拉格朗日乘子法即可得到低維空間上的映射結果。


2.2基于RBF神經網絡的容量估算模型構建


完成特點處理后,本文采用RBF神經網絡開展了鋰離子電池容量衰退模型的構建,RBF神經網絡是一種在函數逼近、數據挖掘等應用場景下前饋神經網絡。其結構簡潔,采用局部響應、局部逼近的原理訓練網絡,與全局響應的BP神經網絡等算法相比有收斂速度更快的優點,并且RBF神經網絡的隱層神經元的激活函數是非線性函數,使得該網絡可以逼近任意非線性函數。其結構由輸入層、隱含層、輸出層等三層構成一個多輸入多輸出的前饋神經網絡,如圖7所示。


圖7RBF神經網絡結構模型


數據由輸入層節點傳遞到隱含層,隱含層節點一般由徑向基函數構成,其中高斯徑向基函數應用效果較好,而輸出節點對隱含層的非線性運算結果進行線性運算為整個網絡的最終結果。


隱含層節點中高斯徑向基函數見式(6)


式中,x為s維輸入向量;ci為第i個徑向基函數的中心;δi為第i個神經元感知范圍,其決定該基函數的中心寬度;m為隱含層神經元的個數;||x-ci||表示x和ci之間的歐氏距離。Y為網絡的輸出向量。


輸入層實現從x→Fi(x)的非線性映射,輸出層實現從Fi(x)到Y的線性映射,即式(7)


式中,Wi為第i個高斯徑向基函數向輸出層傳遞數據所附加的權值;W0為輸出層的偏差項。


采用誤差反向傳播算法對網絡內部參數進行確定,如圖8所示,首先初始化權值矩陣并設定網絡訓練的目標誤差函數及誤差目標,之后將訓練集輸入網絡循環更新權值矩陣,直至達到誤差目標,完成網絡訓練,此時的權值矩陣以及各神經元徑向基函數的參數為最終網絡的結構參數。


圖8RBF神經網訓練流程


鋰離子電池剩余可用容量估計流程如圖9所示,重要分為離線估算模型建立和在線容量估計兩部分。在離線階段,對測得的鋰離子電池循環充放實驗數據進行預處理,作為RBF神經網絡的訓練數據集。從充電階段的時間、電壓、電流、容量、IC曲線中提取出一系列特點向量,通過局部離群因子法剔除特點向量中異常離群點,消除環境噪聲對特點量的干擾,并采用斯皮爾曼相關系數檢驗平滑后特點量與容量之間的相關性。利用局部線性嵌入法對多維特點量進行降維,減小數據的復雜程度,提高模型的訓練效率。最后,訓練RBF網絡中的內部參數,確定容量估算模型。在線部分采用與離線階段相同的方式對數據進行預處理之后作為模型輸入,訓練良好的估算模型輸出高效準確的容量估計值。


圖9鋰離子電池可用容量估算流程圖


3估算結果與討論


3.1模型估算結果


3.1.1不同訓練量對結果的影響


為提高網絡的學習速度及預測精度,將降維后的特點向量進行歸一化處理之后作為RBF神經網絡的輸入值,模型輸出為所估算的放電容量。在數據驅動方法的實現過程中,多采用60%數據用于訓練模型,40%用于測試模型[7]。為評估所建立模型的有效性和魯棒性,本文中,分別劃分每組數據中的50%、60%、70%作為訓練集,其余的作為測試集。輸入網絡模型后,不同訓練量下模型的訓練測試結果如圖10所示,分別采用均方誤差(meansquareerror,MSE)、均方根誤差(rootmeansquareerror,RMSE)、平均絕對誤差(meanabsoluteerror,MAE)對模型預測部分的效果進行分析,如表3所示。


圖10電池A-1、B-2、B-3和B-4訓練測試結果


表3不同電池估算誤差


當60%的數據用來訓練模型時,估計模型對四組實驗電池的剩余可用容量輸出了穩定有效的估算結果,測試區間的預測結果與實際采樣數據非常接近。比較表3中的誤差值,最大MAE出現在電池B-4的預測部分,僅為0.0494,MSE為0.0034,RMSE為0.0581,表明所建立的RBF神經網絡估算模型對A、B兩型號動力鋰離子電池衰退后的容量均能做出較為精準的估算。


由圖10中可以看到,當訓練數據量新增至70%時,模型的預測精度有所提高,最大MAE為0.0381,最大MSE為0.0022,最大RMSE為0.0472。當訓練數據量減少至50%時,模型的估算精度只有略微降低,其中電池B-2的預測誤差較大,最大MAE為0.1067,最大MSE為0.0137,最大RMSE為0.1171,其他電池的預測精度與在60%訓練數據下的預測結果相差無幾。表明所提出的估算模型有良好的泛化能力和較強的魯棒性,能夠在數據有所波動的情況下有效估計不同類型電池的可用容量。


3.1.2不同估算方法比較


為進一步驗證基于特點處理所建立的RBF神經網絡估算模型的在估算精度和效率上的優勢,將Elman神經網絡估算方法、BP神經網絡估算方法和本文提出的估算模型分別應用于估計鋰離子電池全壽命周期的容量衰退情況,容量估算結果、估算誤差及消耗時間如圖11、表4和表5所示。


圖11電池B-3預測電池B-2的結果(a)及電池B-4預測電池B-3的結果(b)


表4電池B-3預測電池B-2的估算誤差及消耗時間


表5電池B-4預測電池B-3的估算誤差及消耗時間


圖11(a)為三種不同模型對電池B-2的全壽命容量預測結果,所有模型的訓練集均來自于電池B-3。圖11(b)顯示的是三種不同模型對電池B-3的全壽命容量預測結果,所有模型的訓練集均來自于電池B-4。由預測結果可知,本文所提出的估算模型對電池B-2的容量預測結果相對誤差除極少數據點在6%之外,其余均小于6%,且表征預測結果離散程度的RMSE僅為0.057,對電池B-3容量預測結果的相對誤差則均在6%以內,顯示出較高的估算精度。


與具有循環結構的Elman神經網絡和傳統的BP神經網絡相比較,Elman神經網絡對電池B-2和電池B-3的容量預測結果整體亦在6%左右,最大MAE為0.0847,最大RMSE為0.0874,BP神經網絡模型的預測結果中相對誤差多數均大于6%,最大MAE為0.0723,最大RMSE為0.0869,而RBF神經網絡模型對兩電池容量預測結果的最大MAE為0.0612,最大RMSE為0.0570,在預測精度和模型穩定方面均體現出良好的實際使用優勢。


此外,由表4可知,在相同的運算處理條件下,在預測電池B-2的容量時,RBF神經網絡所消耗的時間僅為0.99s,Elman神經網絡消耗時間為3.5s,BP神經網絡表現最差,消耗時間為8.5s。根據表5可知,在預測電池B-3的容量時,RBF神經網絡所消耗的時間為0.9s,Elman神經網絡消耗時間為3.5s,BP神經網絡則高達9.4s。比較三種預測模型,可以看出本文所提出的剩余可用容量估算方法所消耗的時間分別約為Elman神經網絡和BP神經網絡的1/5和1/10,顯示出較大的效率優勢。


3.2離群點處理對估算結果的影響


在容量估算的過程中,本文采用LOF算法對特點向量中異常離群點進行精確清洗。為體現不同異常值處理方法對模型估算結果的影響,比較了不處理離群點、應用Hampel濾波器平滑離群點和應用LOF算法清洗異常值三種處理方式對預測結果的影響,如圖12和圖13所示。所預測的剩余可用容量為電池B-3實驗測得,電池A-1和電池B-4的實驗數據為訓練數據集。


圖12不同異常值處理方式下容量估算結果


圖13容量估算結果誤差


結果顯示,未經處理的特點量經LLE算法降維之后輸入RBF模型,其容量預測結果精度較差,最大誤差接近1A·h,且結果分布較為分散,這說明不對異常值進行處理,神經網絡模型不能有效學習到特點量與容量之間的映射關系,得不到理想的預測結果。經Hampel濾波器平滑后的特點量在輸入模型后對剩余可用容量估計的誤差整體在0.2A·h以內,但是,在電池循環壽命后期的估算效果較差。經LOF算法處理后的特點數據對應的容量估計誤差則保持在0.2A·h以內。


在對出現容量回升現象循環的可用容量進行估算時,經Hampel濾波器處理之后所對應的估算結果與經LOF算法處理后所得到的估算結果相比有較大差距,在第113實驗循環,兩者的容量估算絕對誤差分別為0.0253A·h和0.1003A·h,在第143實驗循環分別為0.0158A·h和0.1821A·h。可見經LOF算法處理后的特點數據包含更細致的電池實際容量衰退信息,使估算模型能夠更好的學習到鋰離子電池容量衰退現象與特點量之間的關系,輸出更為準確可靠的剩余可用容量值。


4結論


本文提出一種高效估算鋰離子電池剩余可用容量的方法。由充電階段的數據中提取7組與剩余容量相關特點向量,并采用LOF算法和LLE算法分別對多維特點量進行異常值精準甄別清洗和低維映射,提高特點值所含有效信息量的同時降低估算所需的計算需求,提升估算效率。在此基礎上,利用RBF神經網絡建立了容量估算模型。結果表明,該估算模型能夠有效估算鋰離子電池衰退后的剩余可用容量,與Elman神經網絡模型和BP神經網絡模型相比,在估算精度和估算效率方面均有較大優勢。該方法能夠實現對鋰離子電池衰退后可用容量的高效、準確估算。


引用本文:陳崢,李磊磊,舒星等.基于特點處理與徑向基神經網絡的鋰離子電池剩余容量估算方法[J].儲能科學與技術,2021,10(01):261-270.(CHENZheng,LILeilei,SHUXing,etal.EfficientremainingcapacityestimationmethodforLIBbasedonfeatureprocessingandtheRBFneuralnetwork[J].EnergyStorageScienceandTechnology,2021,10(01):261-270.)


第一作者:陳崢(1982—),男,博士,教授,研究方向為動力鋰離子電池狀態估計,E-mail:chen@kust.edu.cn


第一作者:申江衛,高級實驗師,研究方向為動力鋰離子電池狀態估計,E-mail:shenjiangwei6@163.com。


相關產品

  1. <form id="ejezu"></form>

  2. 主站蜘蛛池模板: 无套内谢老熟女| 油尖旺区| 成人性生交大片免费看中文| 欧美一区二区三区成人片在线 | 冕宁县| 国产又粗又大又黄| 玉山县| 亚洲精品一区二区三区中文字幕| 一本久久综合亚洲鲁鲁五月天| 国产精品久久777777| 少妇高潮惨叫久久久久久| 黑山县| 波多野结衣网站| 佛坪县| 古田县| 久久久精品人妻一区二区三区| 西安市| 湘阴县| 成全观看高清完整免费大全| 性少妇freesexvideos高清| 三年片在线观看大全| 国产成人精品无码免费看夜聊软件 | 精品人妻少妇嫩草av无码专区| 99国产精品久久久久久久成人 | 中文成人在线| 久久国产精品波多野结衣av| √天堂资源地址在线官网| 建平县| 国产精品污www一区二区三区| 工布江达县| 休宁县| 92久久精品一区二区| 人妻激情偷乱视频一区二区三区 | 无码人妻丰满熟妇精品区 | 色一情一乱一伦一区二区三区| 介休市| 欧美日韩在线视频| 日产电影一区二区三区| 国产精品无码免费专区午夜| 内乡县| 人人妻人人澡人人爽人人dvd| 中文字幕乱码中文乱码777| 99国产精品久久久久久久久久久| 三年高清片大全| 辽中县| 巴彦县| 三年成全免费观看影视大全| 会宁县| 麻豆乱码国产一区二区三区| 通山县| 定日县| 新邵县| 国产成人精品aa毛片| 商水县| 成熟妇人a片免费看网站| 平湖市| 米奇影视第四色| 成人小说亚洲一区二区三区| 辣妹子影院电视剧免费播放视频| 阿拉善左旗| 精品久久久久久| 日韩一区二区a片免费观看| 蜜臀av一区二区| 五寨县| 巨大黑人极品videos精品| 欧美性猛交xxxx乱大交蜜桃| 欧美精品videosex极品| 精品久久久久久| 德江县| 国产精品久久久久影院老司| 国产精品久久久一区二区三区| 我们的2018在线观看免费高清| 久久久国产精品人人片| 艳妇乳肉豪妇荡乳| 波多野42部无码喷潮| 夜夜欢天天干| 久久久成人毛片无码| 毛片免费视频| 成人做爰视频www| 樱花影院电视剧免费| 伊通| 荆门市| 国产精品av在线| 欧美精品在线观看| 甘南县| 通化市| 明溪县| 突泉县| 免费人妻精品一区二区三区| 河源市| 娱乐| 酒泉市| 文成县| 国产精品午夜福利视频234区| 嘉定区| 双腿张开被9个男人调教| 久久无码人妻一区二区三区| 精品人人妻人人澡人人爽牛牛| 万山特区| 日韩精品久久久久久免费| 国产真实伦对白全集| 久久久国产精品黄毛片| 欧美疯狂做受xxxxx高潮| 阳城县| 99国产精品久久久久久久成人| 色哟哟网站在线观看| 超碰免费公开| 激情综合五月| 国产伦亲子伦亲子视频观看 | 襄垣县| 国产av天堂| 洱源县| 国产精品丝袜黑色高跟鞋| 国产精品高清网站| 清苑县| 玩弄人妻少妇500系列| 狠狠干狠狠爱| 安图县| 桃源县| 赣榆县| 盱眙县| 唐河县| 香蕉人妻av久久久久天天 | 巴林右旗| 国产熟妇搡bbbb搡bbbb搡| 屏东市| 湘潭市| 99精品一区二区三区无码吞精| 丰满熟妇被猛烈进入高清片| 灌云县| 欧美激情一区二区三区| 精品国产乱码久久久久久影片 | 伦伦影院午夜理论片| 临沧市| 太谷县| 无码国产精品久久一区免费| 久久久成人毛片无码| 礼泉县| 五月丁香啪啪| 建平县| 务川| 国产欧美综合一区二区三区| 娱乐| 延寿县| 播放男人添女人下边视频| 蓝田县| 无码国产精品一区二区色情男同| 东北少妇不戴套对白第一次| 性久久久久久久| 国产精品无码久久久久成人影院 | 蒲江县| 色综合99久久久无码国产精品| 强行无套内谢大学生初次| 三年大片高清影视大全| 长寿区| 堆龙德庆县| 欧美精品18videosex性欧美| 泌阳县| 风韵少妇性饥渴推油按摩视频| 三台县| 岳阳市| 教育| 少妇人妻丰满做爰xxx| 枣庄市| 国产久久精品| 国产精品天天狠天天看| 双腿张开被9个男人调教| 色综合天天综合网国产成人网| 欧美一区二区三区成人久久片| 石屏县| 精品人妻少妇嫩草av无码专区| 免费观看全黄做爰的视频| 国产精品成人3p一区二区三区| 阜新| 乱色精品无码一区二区国产盗| 最好的观看2018中文| 国产成人精品一区二区三区| 成全世界免费高清观看| 成全电影大全在线播放| 激情 小说 亚洲 图片 伦| 解开人妻的裙子猛烈进入| 且末县| 国产精品理论片| 中文字幕人成人乱码亚洲电影| 国产偷窥熟女精品视频大全| 欧美人与性动交α欧美精品| 邵东县| 罗甸县| 成人欧美一区二区三区在线观看| 久久国产精品波多野结衣av| 湖北省| 墨竹工卡县| 中文成人无字幕乱码精品区| 新干县| 精品无码人妻一区二区三区| 睢宁县| 保靖县| 开封市| 新邵县| 石首市| 璧山县| 琼海市| 江门市| 阜阳市| 亚洲字幕av一区二区三区四区| 精品国产乱码久久久久久影片| 河西区| 亚洲男人天堂| 一区二区视频| 国产麻豆剧传媒精品国产av| 欧美午夜精品一区二区蜜桃| 午夜福利电影| 康定县| 性xxxx视频播放免费| 国产精品无码mv在线观看| 涟水县| 99国产精品久久久久久久久久久| 国产精品久久久久久久9999| 国产精品久免费的黄网站| 综合天堂av久久久久久久| 性一交一乱一乱一视频| 熟妇无码乱子成人精品| 那坡县| 国产成人无码一区二区在线观看| 69久久精品无码一区二区| 灵武市| 亚洲熟妇色自偷自拍另类| 无套内谢老熟女| 国内精品一区二区三区| 赤水市| 狠狠人妻久久久久久综合蜜桃| 且末县| 日韩免费视频| 亚洲 欧美 激情 小说 另类| 灵宝市| 色五月激情五月| 湟源县| 双江| 最好的观看2018中文| 又紧又大又爽精品一区二区| 报价| 亚洲中文字幕在线观看| 国产伦精品一区二区三区| 姜堰市| 丰满大肥婆肥奶大屁股| 精人妻无码一区二区三区| 郁南县| 国产伦亲子伦亲子视频观看| 彰武县| 国产视频一区二区| 天天干天天射天天操| 国产精品久久久久久久久久免费看| 成人小说亚洲一区二区三区| 欧美黑人又粗又大高潮喷水| 山西省| 亚洲s码欧洲m码国产av| www夜片内射视频日韩精品成人| 全国最大成人网| 大地资源网在线观看免费动漫| 污污污www精品国产网站| 久久久久久成人毛片免费看| 苍溪县| 久久久久国产精品无码免费看| 影音先锋男人站| 新田县| 精品人妻无码一区二区三区蜜桃一| 色欲狠狠躁天天躁无码中文字幕| 国内老熟妇对白xxxxhd| 日日噜噜噜夜夜爽爽狠狠| 人人妻人人澡人人爽久久av| 五月天激情国产综合婷婷婷 | 英吉沙县| 熟妇无码乱子成人精品| 久久精品中文字幕| 欧美俄罗斯乱妇| 酒泉市| 平泉县| 滨州市| 成全观看高清完整免费大全| 巩留县| 国产又猛又黄又爽| 成全在线观看高清完整版免费动漫| 亚洲熟悉妇女xxx妇女av| 连云港市| 国产精品污www在线观看| 萍乡市| 老熟女高潮一区二区三区| 铅山县| 国产精品99久久久精品无码| 国产草草影院ccyycom| 久久99精品久久久久久琪琪| 精品夜夜澡人妻无码av| 巴林左旗| 亚洲精品一区国产精品| 依兰县| 伦伦影院午夜理论片| 成全在线观看免费完整版| 天堂在线中文| 石门县| 色综合久久88色综合天天| 亚洲亚洲人成综合网络| 铜鼓县| 无套内谢的新婚少妇国语播放| 波多野结衣av在线观看| 狠狠躁18三区二区一区| 开封县| 瓮安县| 亚洲 小说区 图片区 都市| 伊人久久大香线蕉综合75| 国产草草影院ccyycom| 无码国产色欲xxxx视频| 牡丹江市| 无码人妻丰满熟妇奶水区码| 香蕉久久国产av一区二区| 开鲁县| 西贡区| 蒙山县| 欧美与黑人午夜性猛交久久久| 海南省| 西丰县| 午夜精品久久久久久久久| 亚洲精品久久久久久| 欧美亚洲精品suv| 通榆县| 无码少妇一区二区三区| 成人片黄网站色大片免费毛片| 亚洲一区二区三区四区| 性生交大片免费看l| 日本免费视频| 免费特级毛片| 波多野42部无码喷潮在线| 西丰县| 黄瓜视频在线观看| 贡觉县| 和静县| 人妻妺妺窝人体色www聚色窝| 高淳县| 久久精品国产av一区二区三区| 蜜桃久久精品成人无码av| 青青草视频在线观看| 亚洲熟妇色xxxxx欧美老妇| 精品国产成人亚洲午夜福利| 扎囊县| 亚洲中文字幕在线观看| 红原县| 青青草视频在线观看| 哈巴河县| 久久成人无码国产免费播放 | 赣州市| 少妇伦子伦精品无吗| 石河子市| 保靖县| 四会市| 理塘县| 洞头县| 少妇无码一区二区三区| 久久99精品久久久久久| 内射后入在线观看一区| 日产电影一区二区三区| 丰满熟妇被猛烈进入高清片| 湘阴县| 粗大的内捧猛烈进出| 常德市| 乡宁县| 熟妇高潮喷沈阳45熟妇高潮喷| 商水县| 白山市| 五月丁香啪啪| 连平县| 无码一区二区三区在线| 人妻少妇被猛烈进入中文字幕| 欧美激情一区二区| 讷河市| 久久人人爽人人爽人人片 | 红原县| 国产成人精品三级麻豆| 1插菊花综合网| 亚洲无人区码一码二码三码的含义 | 久久久www成人免费精品| 樱花影院电视剧免费| 国产一区二区三区精品视频| 国产精品18久久久| 垣曲县| 双柏县| 国产精品久久久久久亚洲色| 临夏县| 精品国产成人亚洲午夜福利| 齐齐哈尔市| 国产精品美女久久久| 通河县| 宿州市| 鹤山市| 南宁市| 兴业县| 莱阳市| 福州市| 郑州市| 阳西县| 绥滨县| 中文字幕乱码无码人妻系列蜜桃| 久久精品一区二区三区四区| 午夜成人亚洲理伦片在线观看| 视频| 三年大片高清影视大全| 最近免费中文字幕大全免费版视频| 长沙市| 国产精品久久久久无码av色戒| 泽普县| 国产精品久久久国产盗摄| 中文字幕av一区| 国产激情综合五月久久| 国产精品激情| 精国产品一区二区三区a片 | 无码国产伦一区二区三区视频| 国产伦精品一区二区三区妓女下载| 乐陵市| 波多野吉衣av无码| 湟源县| 国产成人精品av| 精品无码人妻一区二区三区品| 熟妇与小伙子matur老熟妇e| 中文字幕人妻无码系列第三区| 久久天天躁狠狠躁夜夜躁2014| 国产精品成人va在线观看| 昌都县| 免费观看黄网站| 泰来县| 定边县| 察哈| 日本欧美久久久久免费播放网| 在线亚洲人成电影网站色www| 熟妇高潮喷沈阳45熟妇高潮喷| 久久久久无码国产精品不卡| 高安市| 格尔木市| 茶陵县| 成人做爰视频www| 麻城市| 国产乱码精品一品二品| 托克逊县| 国产女人被狂躁到高潮小说| 潍坊市| 莱阳市| 禄劝| 中宁县| 安溪县| 人妻洗澡被强公日日澡| √天堂资源地址在线官网| 精品无码人妻一区二区免费蜜桃| 丰满岳乱妇一区二区三区| 一边摸一边做爽的视频17国产| 真实的国产乱xxxx在线| 久久午夜无码鲁丝片午夜精品| 无码人妻丰满熟妇精品区| 怀仁县| 风韵少妇性饥渴推油按摩视频| 丁香五香天堂网| 亚洲日韩一区二区三区| 阳江市| 国产露脸无套对白在线播放| 解开人妻的裙子猛烈进入| 久久精品噜噜噜成人| 国精产品一区一区三区| 手机| 潼南县| 高雄县| 亚洲精品成a人在线观看| 欧美激情综合色综合啪啪五月| 国产乱码精品一品二品| 躁老太老太騷bbbb| 久久久久99精品国产片 | 欧美午夜精品一区二区三区电影| 国产成人无码一区二区在线观看| 久久丫精品忘忧草西安产品| 久久久久久久久毛片无码| 巴青县| 中文字幕乱码无码人妻系列蜜桃 | 亚洲人成在线观看| 精品国产乱码久久久久久郑州公司| 中文字幕无码毛片免费看| 无码少妇一区二区三区| 丹凤县| 兴宁市| 雷波县| 无码一区二区三区视频| 粗大黑人巨精大战欧美成人| 白玉县| 日韩伦人妻无码| 辛集市| 亚洲の无码国产の无码步美| 西安市| 人妻巨大乳一二三区| 怀宁县| 日韩精品无码一区二区三区| 新蔡县| 国偷自产视频一区二区久| 三年成全免费观看影视大全| 远安县| 欧美三根一起进三p| 狂野少女电视剧免费播放| 性视频播放免费视频| 东北少妇不戴套对白第一次| 黔西县| 欧美最猛黑人xxxx黑人猛交 | 中文字幕精品无码一区二区| 无码日本精品xxxxxxxxx| 高潮毛片又色又爽免费| 最近免费中文字幕大全免费版视频| 初尝黑人巨砲波多野结衣| 色噜噜狠狠一区二区三区果冻| 国产成人精品一区二区三区免费| 曲水县| 无套内谢老熟女| 性做久久久久久| 上思县| 郁南县| 张北县| 民权县| 来安县| 济南市| 利津县| 特大黑人娇小亚洲女| 国产欧美综合一区二区三区| 风韵丰满熟妇啪啪区老熟熟女| 称多县| 日本不卡一区| 色欲一区二区三区精品a片 | 九台市| 启东市| 亚洲日韩一区二区| 国产aⅴ激情无码久久久无码| 国产伦精品一区二区三区免费迷| 欧美成人片在线观看| 强伦人妻一区二区三区视频18 | 万安县| 欧美俄罗斯乱妇| 亚洲欧美精品午睡沙发| 成全影院高清电影好看的电视剧| 中国极品少妇xxxxx| 久久精品一区二区三区四区| 日韩av无码一区二区三区不卡| 成全免费高清观看在线电视剧大全 | 国产女人高潮毛片| 乳尖春药h糙汉共妻| 国产suv精品一区二区| 亚洲 小说 欧美 激情 另类| 五家渠市| 精品国产乱码久久久久久影片 | 金华市| 国产做爰xxxⅹ久久久精华液| 亚洲熟女乱色综合亚洲小说| 亚洲 小说 欧美 激情 另类| 绥宁县| 北海市| 湖南省| 少妇厨房愉情理伦bd在线观看| 国产免费视频| 亚洲永久无码7777kkk| 深泽县| 久久丫精品久久丫| 迭部县| 化隆| 景泰县| 日产精品久久久一区二区| 阳原县| 久久久国产精品| 寿阳县| 东辽县| 中文毛片无遮挡高潮免费| 一边摸一边抽搐一进一出视频| 河池市| 桐柏县| 1插菊花综合网| 全部孕妇毛片丰满孕妇孕交| 江北区| 清河县| 国産精品久久久久久久| 欧美性大战xxxxx久久久| 无码一区二区三区在线| gogogo免费视频观看| 国产电影一区二区三区| 久久久久噜噜噜亚洲熟女综合| 无码人妻久久一区二区三区蜜桃 | 国产真实乱人偷精品视频| 久久久久成人精品无码中文字幕 | 国产精品无码免费播放| 西乌| 99久久国产热无码精品免费| 亚洲精品久久久久久一区二区| 无码一区二区三区| 隆化县| 99精品一区二区三区无码吞精| 无码h肉动漫在线观看| 白城市| 99久久99久久精品免费看蜜桃| 日本公妇乱偷中文字幕| 1插菊花综合网| 永和县| 姜堰市| 邓州市| 麻豆国产一区二区三区四区| 久久er99热精品一区二区| 中文字幕乱码在线人视频| 乌海市| 无码人妻丰满熟妇区bbbbxxxx| 欧美精品18videosex性欧美| 玛纳斯县| 扬州市| 无码av免费精品一区二区三区| 久久午夜夜伦鲁鲁片无码免费| 湘潭市| 天门市| 欧美性猛交xxxx免费看| 青川县| 广平县| 老司机午夜福利视频| 六安市| 人妻体内射精一区二区| 长乐市| 夜夜躁很很躁日日躁麻豆| 二连浩特市| 日韩精品人妻中文字幕有码| 巢湖市| 人妻熟女一区二区三区app下载| gogogo在线高清免费完整版| 卢氏县| 日韩成人无码| 熟妇高潮一区二区在线播放| 洛浦县| 成全在线观看免费完整版| 三人成全免费观看电视剧| 搡老岳熟女国产熟妇| 沂源县| 台湾省| 国产精品永久久久久久久久久| 欧美性生交大片免费看| 艳妇乳肉豪妇荡乳av无码福利| 隆子县| 精品国产乱码久久久久久影片| 广河县| 三年大片高清影视大全| 神农架林区| 涪陵区| 亚洲 小说区 图片区 都市| 广丰县| 会同县| 无码一区二区三区视频| 中文字幕无码精品亚洲35| 国产真人做爰毛片视频直播| 铅山县| 上杭县| 大洼县| 又白又嫩毛又多15p| 拜城县| 武陟县| 日本特黄特色aaa大片免费| 丽江市| 厦门市| 亚洲精品97久久中文字幕无码| 精品国产成人亚洲午夜福利| 成全我在线观看免费观看| 云梦县| 国产精成人品| 中国白嫩丰满人妻videos | 97香蕉碰碰人妻国产欧美| 黑人巨大精品欧美一区二区| 平江县| 国产一区二区三区精品视频| 成全在线观看免费完整| 99久久人妻精品免费二区| 久久久久国产精品无码免费看| 富蕴县| 久久久久国产精品无码免费看| 国产99久一区二区三区a片| av片在线观看| 久久丫精品忘忧草西安产品| 肉大榛一进一出免费视频| 英德市| 精国产品一区二区三区a片| 工布江达县| 尤物视频在线观看| 亚洲啪av永久无码精品放毛片| 国产亚州精品女人久久久久久| 日韩精品极品视频在线观看免费| 尉犁县| 蕲春县| 大关县| 广河县| 垦利县| 涞水县| 成人小说亚洲一区二区三区| 成全视频大全高清全集在线| √天堂资源地址在线官网| 国产精品国产三级国产专区53| 无码精品黑人一区二区三区| 亚洲精品国偷拍自产在线观看蜜桃| 艳妇乳肉豪妇荡乳av无码福利| 欧美裸体xxxx极品少妇| 精品人妻伦一二三区久久| 盐城市| 无码人妻aⅴ一区二区三区69岛| 禹城市| 扬州市| 国产人妻人伦精品1国产 | 潜山县| 亚洲人成色777777精品音频 | 武陟县| 伦伦影院午夜理论片| 日产电影一区二区三区| 亚洲男人天堂| 欧美乱大交| 国产精品99无码一区二区| 南投市| 日本边添边摸边做边爱| 国产福利视频| 国产激情久久久久久熟女老人av | 精品人妻无码一区二区三区| 三门县| 祁连县| 成人网站在线进入爽爽爽| 青青草原亚洲| 欧美 日韩 国产 成人 在线观看| 冀州市| 蓝田县| 国产免费视频| 精品人妻少妇嫩草av无码专区 | 西充县| 亚洲精品97久久中文字幕无码| 成全影视在线观看第6季| 巩留县| 少妇粉嫩小泬白浆流出| 国产乱xxⅹxx国语对白| 白水县| 呼图壁县| 成人欧美一区二区三区黑人免费| 国产精品美女久久久久久久久| 基隆市| 在线天堂www在线国语对白|