999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

基于中值觀測的自適應變步長MPPT算法研究

2014-11-14 23:36:57徐彭濤趙巧娥張強張麗娜
現代電子技術 2014年22期

徐彭濤+趙巧娥+張強+張麗娜

摘 要: 提出一種應用于光伏發電系統中的基于中值觀測的自適應變步長MPPT方法。該MPPT方法以中值法為觀測基礎,具有自適應且變步長的優勢,跟蹤過程不需要人為干涉與尋找初值和優值,在具有高效率的同時具有很高的適用性。在Matlab/Simulink平臺驗證了此MPPT控制算法,其模型避免用到除法器,有利于加快運算速度與減小電路規模。

關鍵詞: 光伏發電; 中值判斷; 自適應變步長; Matlab/Simulink

中圖分類號: TN911?34 文獻標識碼: A 文章編號: 1004?373X(2014)22?0150?03

Study on adaptive variable?step MPPT algorithm based on mid?value observation

XU Peng?tao1, ZHAO Qiao?e1, ZHANG Qiang2, ZHANG Li?na1

(1. Shanxi University, Taiyuan 030013, China; 2. Metrology Center, Shanxi Electric Power Company, TaiYuan 030001, China)

Abstracts: An adaptive variable?step MPPT (maximum power point tracking) algorithm based on mid?value observation is proposed, which is applied to the photovoltaic power generation system. It bases on mid?value observation and owns the advantages of adaptive variable?step. It does not need any man?made interference or define parameters during tracking. Its high efficiency and applicability are outstanding. This MTTP algorithm is verified in Matlab/Simulink platform. There is no divider in the model, so it is beneficial to accelerate tracking and simplify circuit.

Keywords: photovoltaic power generation; mid?value observation; self?adaption variable?step; Matlab/Simulink

利用太陽能光伏發電是新能源的重要發展方向之一,光伏發電系統由太陽能電池板、控制器、蓄電池、逆變器四個部分組成??刂破骺刂普麄€系統的工作狀態,為了改善太陽能電池的轉換效率,需要其對輸出的最大功率實施跟蹤。太陽能電池作為非線性元件,其輸出易受光照強度、溫度等的影響,最大功率點隨著光照強度和環境溫度的變化而變化。如果光伏陣列不能實時保持工作在最大功率點上,將導致轉換效率降低[1]。為了提高系統的轉換效率,需要使光伏陣列的輸出阻抗和負載阻抗實時自動調節達到匹配狀態,這一技術叫做最大功率點跟蹤(Maximum Power Point Tracking,MPPT)。

本文以最大功率點跟蹤技術為研究方向,提出一種基于中值觀測的自適應變步長MPPT方法。該方法以中值法為觀測與判斷基礎,跟蹤過程是變步長的且步長自適應不需要人為尋找優值。在之后 Matlab/Simulink驗證中可以發現,相比于其他的變步長算法,此方法減少了乘法器、特別是避免用到除法器,加快了運算速度,有效減小了實現電路規模。

1 光伏電池的輸出特性

光伏電池本質上為光電子器件,其利用光生伏特效應將光能轉化為電能。根據光生伏特器件的基本原理,光伏電池相當于一個恒流源J1和一個二極管D1并聯。其等效電路如圖1所示,考慮包含電池內部溫度(結溫)T與光照強度S對最終輸出的影響,根據半導體器件與電路原理得到如下關系:

[Ipv=Isc·S1 000+ID·expq·Vpv+IpvRskT- Vpv+IpvRsRsh]

式中:Isc是光伏電池短路電流[2];S是光照強度;T是電池溫度(單位為開爾文);ID是二極管本征反向電流;k是波爾茲曼常數;q是電子電量;Ipv是光伏電池陣列的輸出電流;Vpv是光伏電池陣列的輸出電壓;Rsh是分流電阻,表征光伏電池內阻;Rs是串聯電阻,等效于光伏電池內部陣列之間及外部到負載傳輸路徑上的傳輸損耗。

利用Matlab仿真來闡明光照強度和溫度對輸出工作點的影響,不同光照和溫度下的輸出電壓、輸出電流和輸出功率之間的關系如圖2和圖3所示。

由圖2可見,隨著光照強度的增加,光伏陣列的輸出電流隨之增加從而導致最大輸出功率點也相應變化。而圖3說明了隨著環境溫度的變化,最大輸出功率點也相應變化。

圖1 光伏電池等效電路

圖2 不同光照下的I?V關系

圖3 不同溫度下的I?V關系

2 MPPT電壓控制

由上述仿真得到的I?V和P?V關系可以看出,光伏電池陣列的最佳工作點是隨著光照強度和溫度隨時變化的,需要改變電壓Vpv使得輸出功率Ppv達到最大。MPPT電壓控制常用的電路包括Buck(降壓)、Boost(升壓)、Buck?Boost(降壓?升壓)等類型[3]。文中所用為Boost升壓電路作為光伏陣列的電壓控制系統。Boost升壓電路如圖4示。

圖4 Boost升壓電路的電路圖

設開關管S通斷周期為T,占空比為D。當S導通時,電感電壓為Vpv,導通時間為DT,;當S斷開時,此時光伏電池和電感L同時給C充電,充電電壓為(Vpv-Vdc),關斷時間T·(1-D)。電感L在一個周期內的電流平衡、電容C穩定Vdc,有關系式如下:

[VpvL·D·T+Vpv-VdcL·1-D·T=0Vpv=Vdc·1-D]

可見,通過改變占空比D即可改變光伏陣列的電壓Vpv,從而使光伏陣列輸出最大的功率,完成最大輸出功率點跟蹤過程。

3 基于中值觀測的自適應變步MPPT方法

為了提高系統的轉換效率,光伏陣列需要實時工作在最大輸出工作點,即最大功率點跟蹤MPPT。傳統的MPPT控制方法有擾動觀察法和電導增量法。擾動觀測法結構簡單,需要人工定義的參數少,但是存在振蕩問題和跟蹤精度速度問題。雖然出現了一些改進的變步長算法[4?5],但是存在著步長改變速度需要人為優化設定,變步長存在速度與精度矛盾問題。電導增量法的優點是精度高,響應速度較快,其缺點是需要手動選擇合適的步長,且運用了除法器,大大影響了跟蹤速度[6]。另外出現了一些在新思想基礎上發展出的MPPT算法如模糊控制[7]、粒子群[8]、SVM 預測模型[9]等算法,亦存在著實現困難,人工設定參數復雜的問題。本文以中值法為觀測與判斷基礎,提出了一種在跟蹤最大功率點過程中變步長且步長自適應的MPPT方法。其步長不僅能夠隨著工作點離最大功率輸出點的遠近而改變,而且變步長過程不必需要人為干涉與尋找初值和優值,在具有高效率的同時具有很高的適用性。此算法的流程圖見圖5。

(1) 在步驟A1中,采樣Vpv1(i),Vpv2(i),Vpv3(i)。以D為單位步進,K為步進方向,此時序的步進變化后的電壓Vpv3(i)=Vpv1(i)+2nK(i)D,并取中值Vpv2(i)=[Vpv1(i)+Vpv3(i)[]2]。根據P?V關系采樣相應的Ppv1(i),Ppv2(i),Ppv3(i)。

(2) 步驟A2,A3,A4判斷中值Ppv2(i)與Ppv1(i),Ppv3(i)的關系Ppv1(i)

(3) 步驟A2,B3,C4判斷中值Ppv2(i)與Ppv1(i),Ppv3(i)的關系Ppv1(i)>Ppv2(i)>Ppv3(i),此時最大功率輸出工作點在P?V曲線的右側,定義下次步進方向為-1即向左。

(4) 步驟A2,A3,B4和步驟A2,B3,B4判斷中值Ppv2(i)與Ppv1(i),Ppv3(i)的關系為Ppv1(i)Ppv3(i),此時最大功率輸出工作點在此次步進跨國的步長內,定義下次步進方向為K=-K即反向。

(5) 步驟A5判斷在前幾個時序中步進方向是否一致。若K(i-3)=K(i-2)=K(i-1)=K(i),則如步驟B6步長加倍n=n+1;如果不一致,在步驟A6中判斷K(i)與上一時序的K(i-1)是否相等,如果不相等n=n-1。

可見,算法明顯以中值Vpv2(i)為觀測與判斷基礎,步進2n在每個時序中根據相應情況自動取合適的改變量,達到了變步長和自適應的效果。

圖5 算法流程圖

4 實驗驗證與結果分析

為了驗證算法的可行性,本文通過Matlab/Simulink工具對流程算法進行系統搭建與仿真。仿真系統的模型分別見圖6、圖7。仿真結果見圖8。

圖6 Boost電路仿真模型

仿真模型左側是光伏陣列在一定光強S和溫度T條件下的P?V對應關系。此MPPT算法模型中,明顯的避免用到除法器,有利于加快運算速度與減小電路規模。為了驗證算法的可行性,分別設定下列幾個條件:

(1) 設定在7.8 s時光照強度從1 000 W/m2下降到600 W/m2;

(2) 設定在7.8 s時光照強度從600 W/m2上升到1 000 W/m2;

(3) 設定在7.8 s時溫度從1 000 K下降到300 K;

(4) 設定在7.8 s時溫度從300 K上升到1 000 K。

圖7 MPPT仿真模型

圖8 Boost升壓電路占空比Simulink仿真結果

Simulink仿真結果如圖9所示。

圖9 Simulink仿真結果

可見在光照強度或者溫度出現突然變化時,無論是正向還是負向應用本文的MPPT算法都能快速準確的響應環境變化,且變化過程無人為干涉步長自適應地變化,驗證了此算法的可行性和有效性。

5 結 論

最大功率點跟蹤技術能有效提高系統的轉換效率,針對MPPT研究中存在的擾動法步長固定、某些改進的變步長參數需要人為優化、電導增量法存在除法運算等問題提出了一種基于中值觀測的自適應變步長MPPT方法。此MPPT方法以中值法為觀測與判斷基礎,最大的特點是跟蹤過程為變步長且步長自適應不需要人為尋找優值。相比于其他的變步長算法,此方法減少了乘法器、特別是避免用到除法器,加快了運算速度,有效減小了實現電路規模。

參考文獻

[1] 艾欣,韓曉男,孫英云.光伏發電并網及其相關技術發展現狀與展望[J].現代電力,2013,30(1):1?7.

[2] 茆美琴,余世杰,蘇建徽.帶有MPPT功能的光伏陣列Matlab通用仿真模型[J].系統仿真學報,2005,17(5):1248?1251.

[3] 王兆按,黃俊.電力電子技術[M].北京:機械工業出版社,2009.

[4] 沈實疊,姚維.基于變步長的光伏系統MPPT 算法研究[J].輕工機械,2013,31(5):14?17.

[5] 路曉,秦立軍.自適應擾動觀察法在光伏MPPT中的應用與仿真[J].現代電力,2011,28(1):80?84.

[6] 李鋒,劉春生,徐曉淳.光伏發電系統改進復合算法[J].計算機與現代化,2013,8(3):15?19.

[7] 王云飛,尹忠東,申燕飛.基于模糊控制的MPPT在光伏并網發電系統中的應用[J].現代電力,2011,28(4):54?58.

[8] 王平,周沖.基于粒子群算法的光伏電池MPPT控制策略[J].可再生能源,2013,31(5):26?28.

[9] 陳文穎,林永君,楊春來,等.基于SVM預測模型的光伏發電系統MPPT研究[J].太陽能學報,2013,34(2):67?69.

圖4 Boost升壓電路的電路圖

設開關管S通斷周期為T,占空比為D。當S導通時,電感電壓為Vpv,導通時間為DT,;當S斷開時,此時光伏電池和電感L同時給C充電,充電電壓為(Vpv-Vdc),關斷時間T·(1-D)。電感L在一個周期內的電流平衡、電容C穩定Vdc,有關系式如下:

[VpvL·D·T+Vpv-VdcL·1-D·T=0Vpv=Vdc·1-D]

可見,通過改變占空比D即可改變光伏陣列的電壓Vpv,從而使光伏陣列輸出最大的功率,完成最大輸出功率點跟蹤過程。

3 基于中值觀測的自適應變步MPPT方法

為了提高系統的轉換效率,光伏陣列需要實時工作在最大輸出工作點,即最大功率點跟蹤MPPT。傳統的MPPT控制方法有擾動觀察法和電導增量法。擾動觀測法結構簡單,需要人工定義的參數少,但是存在振蕩問題和跟蹤精度速度問題。雖然出現了一些改進的變步長算法[4?5],但是存在著步長改變速度需要人為優化設定,變步長存在速度與精度矛盾問題。電導增量法的優點是精度高,響應速度較快,其缺點是需要手動選擇合適的步長,且運用了除法器,大大影響了跟蹤速度[6]。另外出現了一些在新思想基礎上發展出的MPPT算法如模糊控制[7]、粒子群[8]、SVM 預測模型[9]等算法,亦存在著實現困難,人工設定參數復雜的問題。本文以中值法為觀測與判斷基礎,提出了一種在跟蹤最大功率點過程中變步長且步長自適應的MPPT方法。其步長不僅能夠隨著工作點離最大功率輸出點的遠近而改變,而且變步長過程不必需要人為干涉與尋找初值和優值,在具有高效率的同時具有很高的適用性。此算法的流程圖見圖5。

(1) 在步驟A1中,采樣Vpv1(i),Vpv2(i),Vpv3(i)。以D為單位步進,K為步進方向,此時序的步進變化后的電壓Vpv3(i)=Vpv1(i)+2nK(i)D,并取中值Vpv2(i)=[Vpv1(i)+Vpv3(i)[]2]。根據P?V關系采樣相應的Ppv1(i),Ppv2(i),Ppv3(i)。

(2) 步驟A2,A3,A4判斷中值Ppv2(i)與Ppv1(i),Ppv3(i)的關系Ppv1(i)

(3) 步驟A2,B3,C4判斷中值Ppv2(i)與Ppv1(i),Ppv3(i)的關系Ppv1(i)>Ppv2(i)>Ppv3(i),此時最大功率輸出工作點在P?V曲線的右側,定義下次步進方向為-1即向左。

(4) 步驟A2,A3,B4和步驟A2,B3,B4判斷中值Ppv2(i)與Ppv1(i),Ppv3(i)的關系為Ppv1(i)Ppv3(i),此時最大功率輸出工作點在此次步進跨國的步長內,定義下次步進方向為K=-K即反向。

(5) 步驟A5判斷在前幾個時序中步進方向是否一致。若K(i-3)=K(i-2)=K(i-1)=K(i),則如步驟B6步長加倍n=n+1;如果不一致,在步驟A6中判斷K(i)與上一時序的K(i-1)是否相等,如果不相等n=n-1。

可見,算法明顯以中值Vpv2(i)為觀測與判斷基礎,步進2n在每個時序中根據相應情況自動取合適的改變量,達到了變步長和自適應的效果。

圖5 算法流程圖

4 實驗驗證與結果分析

為了驗證算法的可行性,本文通過Matlab/Simulink工具對流程算法進行系統搭建與仿真。仿真系統的模型分別見圖6、圖7。仿真結果見圖8。

圖6 Boost電路仿真模型

仿真模型左側是光伏陣列在一定光強S和溫度T條件下的P?V對應關系。此MPPT算法模型中,明顯的避免用到除法器,有利于加快運算速度與減小電路規模。為了驗證算法的可行性,分別設定下列幾個條件:

(1) 設定在7.8 s時光照強度從1 000 W/m2下降到600 W/m2;

(2) 設定在7.8 s時光照強度從600 W/m2上升到1 000 W/m2;

(3) 設定在7.8 s時溫度從1 000 K下降到300 K;

(4) 設定在7.8 s時溫度從300 K上升到1 000 K。

圖7 MPPT仿真模型

圖8 Boost升壓電路占空比Simulink仿真結果

Simulink仿真結果如圖9所示。

圖9 Simulink仿真結果

可見在光照強度或者溫度出現突然變化時,無論是正向還是負向應用本文的MPPT算法都能快速準確的響應環境變化,且變化過程無人為干涉步長自適應地變化,驗證了此算法的可行性和有效性。

5 結 論

最大功率點跟蹤技術能有效提高系統的轉換效率,針對MPPT研究中存在的擾動法步長固定、某些改進的變步長參數需要人為優化、電導增量法存在除法運算等問題提出了一種基于中值觀測的自適應變步長MPPT方法。此MPPT方法以中值法為觀測與判斷基礎,最大的特點是跟蹤過程為變步長且步長自適應不需要人為尋找優值。相比于其他的變步長算法,此方法減少了乘法器、特別是避免用到除法器,加快了運算速度,有效減小了實現電路規模。

參考文獻

[1] 艾欣,韓曉男,孫英云.光伏發電并網及其相關技術發展現狀與展望[J].現代電力,2013,30(1):1?7.

[2] 茆美琴,余世杰,蘇建徽.帶有MPPT功能的光伏陣列Matlab通用仿真模型[J].系統仿真學報,2005,17(5):1248?1251.

[3] 王兆按,黃俊.電力電子技術[M].北京:機械工業出版社,2009.

[4] 沈實疊,姚維.基于變步長的光伏系統MPPT 算法研究[J].輕工機械,2013,31(5):14?17.

[5] 路曉,秦立軍.自適應擾動觀察法在光伏MPPT中的應用與仿真[J].現代電力,2011,28(1):80?84.

[6] 李鋒,劉春生,徐曉淳.光伏發電系統改進復合算法[J].計算機與現代化,2013,8(3):15?19.

[7] 王云飛,尹忠東,申燕飛.基于模糊控制的MPPT在光伏并網發電系統中的應用[J].現代電力,2011,28(4):54?58.

[8] 王平,周沖.基于粒子群算法的光伏電池MPPT控制策略[J].可再生能源,2013,31(5):26?28.

[9] 陳文穎,林永君,楊春來,等.基于SVM預測模型的光伏發電系統MPPT研究[J].太陽能學報,2013,34(2):67?69.

圖4 Boost升壓電路的電路圖

設開關管S通斷周期為T,占空比為D。當S導通時,電感電壓為Vpv,導通時間為DT,;當S斷開時,此時光伏電池和電感L同時給C充電,充電電壓為(Vpv-Vdc),關斷時間T·(1-D)。電感L在一個周期內的電流平衡、電容C穩定Vdc,有關系式如下:

[VpvL·D·T+Vpv-VdcL·1-D·T=0Vpv=Vdc·1-D]

可見,通過改變占空比D即可改變光伏陣列的電壓Vpv,從而使光伏陣列輸出最大的功率,完成最大輸出功率點跟蹤過程。

3 基于中值觀測的自適應變步MPPT方法

為了提高系統的轉換效率,光伏陣列需要實時工作在最大輸出工作點,即最大功率點跟蹤MPPT。傳統的MPPT控制方法有擾動觀察法和電導增量法。擾動觀測法結構簡單,需要人工定義的參數少,但是存在振蕩問題和跟蹤精度速度問題。雖然出現了一些改進的變步長算法[4?5],但是存在著步長改變速度需要人為優化設定,變步長存在速度與精度矛盾問題。電導增量法的優點是精度高,響應速度較快,其缺點是需要手動選擇合適的步長,且運用了除法器,大大影響了跟蹤速度[6]。另外出現了一些在新思想基礎上發展出的MPPT算法如模糊控制[7]、粒子群[8]、SVM 預測模型[9]等算法,亦存在著實現困難,人工設定參數復雜的問題。本文以中值法為觀測與判斷基礎,提出了一種在跟蹤最大功率點過程中變步長且步長自適應的MPPT方法。其步長不僅能夠隨著工作點離最大功率輸出點的遠近而改變,而且變步長過程不必需要人為干涉與尋找初值和優值,在具有高效率的同時具有很高的適用性。此算法的流程圖見圖5。

(1) 在步驟A1中,采樣Vpv1(i),Vpv2(i),Vpv3(i)。以D為單位步進,K為步進方向,此時序的步進變化后的電壓Vpv3(i)=Vpv1(i)+2nK(i)D,并取中值Vpv2(i)=[Vpv1(i)+Vpv3(i)[]2]。根據P?V關系采樣相應的Ppv1(i),Ppv2(i),Ppv3(i)。

(2) 步驟A2,A3,A4判斷中值Ppv2(i)與Ppv1(i),Ppv3(i)的關系Ppv1(i)

(3) 步驟A2,B3,C4判斷中值Ppv2(i)與Ppv1(i),Ppv3(i)的關系Ppv1(i)>Ppv2(i)>Ppv3(i),此時最大功率輸出工作點在P?V曲線的右側,定義下次步進方向為-1即向左。

(4) 步驟A2,A3,B4和步驟A2,B3,B4判斷中值Ppv2(i)與Ppv1(i),Ppv3(i)的關系為Ppv1(i)Ppv3(i),此時最大功率輸出工作點在此次步進跨國的步長內,定義下次步進方向為K=-K即反向。

(5) 步驟A5判斷在前幾個時序中步進方向是否一致。若K(i-3)=K(i-2)=K(i-1)=K(i),則如步驟B6步長加倍n=n+1;如果不一致,在步驟A6中判斷K(i)與上一時序的K(i-1)是否相等,如果不相等n=n-1。

可見,算法明顯以中值Vpv2(i)為觀測與判斷基礎,步進2n在每個時序中根據相應情況自動取合適的改變量,達到了變步長和自適應的效果。

圖5 算法流程圖

4 實驗驗證與結果分析

為了驗證算法的可行性,本文通過Matlab/Simulink工具對流程算法進行系統搭建與仿真。仿真系統的模型分別見圖6、圖7。仿真結果見圖8。

圖6 Boost電路仿真模型

仿真模型左側是光伏陣列在一定光強S和溫度T條件下的P?V對應關系。此MPPT算法模型中,明顯的避免用到除法器,有利于加快運算速度與減小電路規模。為了驗證算法的可行性,分別設定下列幾個條件:

(1) 設定在7.8 s時光照強度從1 000 W/m2下降到600 W/m2;

(2) 設定在7.8 s時光照強度從600 W/m2上升到1 000 W/m2;

(3) 設定在7.8 s時溫度從1 000 K下降到300 K;

(4) 設定在7.8 s時溫度從300 K上升到1 000 K。

圖7 MPPT仿真模型

圖8 Boost升壓電路占空比Simulink仿真結果

Simulink仿真結果如圖9所示。

圖9 Simulink仿真結果

可見在光照強度或者溫度出現突然變化時,無論是正向還是負向應用本文的MPPT算法都能快速準確的響應環境變化,且變化過程無人為干涉步長自適應地變化,驗證了此算法的可行性和有效性。

5 結 論

最大功率點跟蹤技術能有效提高系統的轉換效率,針對MPPT研究中存在的擾動法步長固定、某些改進的變步長參數需要人為優化、電導增量法存在除法運算等問題提出了一種基于中值觀測的自適應變步長MPPT方法。此MPPT方法以中值法為觀測與判斷基礎,最大的特點是跟蹤過程為變步長且步長自適應不需要人為尋找優值。相比于其他的變步長算法,此方法減少了乘法器、特別是避免用到除法器,加快了運算速度,有效減小了實現電路規模。

參考文獻

[1] 艾欣,韓曉男,孫英云.光伏發電并網及其相關技術發展現狀與展望[J].現代電力,2013,30(1):1?7.

[2] 茆美琴,余世杰,蘇建徽.帶有MPPT功能的光伏陣列Matlab通用仿真模型[J].系統仿真學報,2005,17(5):1248?1251.

[3] 王兆按,黃俊.電力電子技術[M].北京:機械工業出版社,2009.

[4] 沈實疊,姚維.基于變步長的光伏系統MPPT 算法研究[J].輕工機械,2013,31(5):14?17.

[5] 路曉,秦立軍.自適應擾動觀察法在光伏MPPT中的應用與仿真[J].現代電力,2011,28(1):80?84.

[6] 李鋒,劉春生,徐曉淳.光伏發電系統改進復合算法[J].計算機與現代化,2013,8(3):15?19.

[7] 王云飛,尹忠東,申燕飛.基于模糊控制的MPPT在光伏并網發電系統中的應用[J].現代電力,2011,28(4):54?58.

[8] 王平,周沖.基于粒子群算法的光伏電池MPPT控制策略[J].可再生能源,2013,31(5):26?28.

[9] 陳文穎,林永君,楊春來,等.基于SVM預測模型的光伏發電系統MPPT研究[J].太陽能學報,2013,34(2):67?69.

主站蜘蛛池模板: 亚洲精品国偷自产在线91正片| 国产精品网曝门免费视频| 欧美日韩精品一区二区视频| 成人综合网址| 一级成人a毛片免费播放| 天天躁日日躁狠狠躁中文字幕| 国产成人精品午夜视频'| 精品久久久久久中文字幕女| 伊人福利视频| 欧美第一页在线| 扒开粉嫩的小缝隙喷白浆视频| 无码一区18禁| 热99精品视频| 四虎影视国产精品| 四虎国产在线观看| 欧洲熟妇精品视频| 亚洲无码91视频| 国产黄色视频综合| 欧美另类精品一区二区三区| 激情综合激情| 日韩欧美中文字幕在线精品| 美女潮喷出白浆在线观看视频| 内射人妻无套中出无码| 美女扒开下面流白浆在线试听| 亚洲男人的天堂久久精品| 欧美日本二区| 激情国产精品一区| 国产色婷婷| 毛片视频网| 精品福利网| 亚洲成a人片| 免费又爽又刺激高潮网址| 国产va欧美va在线观看| 99热这里只有精品国产99| 国产精品漂亮美女在线观看| 亚洲天堂2014| 91蜜芽尤物福利在线观看| 亚洲成人精品| 久久香蕉国产线| 嫩草在线视频| 国产小视频网站| 国产精品亚欧美一区二区| 亚洲成人免费看| 国产亚洲现在一区二区中文| 色婷婷啪啪| 91青青视频| 国产精品久久久久久影院| av无码久久精品| 操国产美女| 国产精品自在在线午夜| 亚洲另类国产欧美一区二区| 国产精品亚欧美一区二区三区 | 岛国精品一区免费视频在线观看 | 动漫精品啪啪一区二区三区| 99久久精品免费看国产电影| 亚洲香蕉在线| 天天视频在线91频| 国产精品久久久精品三级| 青青青视频免费一区二区| 欧美亚洲中文精品三区| 久久96热在精品国产高清| 国产白浆一区二区三区视频在线| 国产成本人片免费a∨短片| 成人亚洲国产| 日a本亚洲中文在线观看| 男人天堂伊人网| 亚洲不卡无码av中文字幕| 亚洲全网成人资源在线观看| 国产精品制服| 首页亚洲国产丝袜长腿综合| 一区二区欧美日韩高清免费| 青草91视频免费观看| 国产一级裸网站| 亚洲成年网站在线观看| 依依成人精品无v国产| 国产成人福利在线视老湿机| 国产在线八区| 日韩黄色大片免费看| 制服丝袜无码每日更新| 成人中文字幕在线| 亚洲狼网站狼狼鲁亚洲下载| 动漫精品啪啪一区二区三区|