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

風電變槳電機參數對控制系統影響分析及改進辨識方法

2019-08-05 07:20:49李輝謝翔杰劉行中姚然柴兆森夏桂森
電機與控制學報 2019年7期
關鍵詞:風力發電

李輝 謝翔杰 劉行中 姚然 柴兆森 夏桂森

摘?要:針對風電機組變槳電機參數變化和負載不確定等帶來變槳控制系統魯棒性差的問題,分析變槳電機參數對控制系統的影響并提出參數辨識改進方法。首先,基于變槳電機控制系統中的復合磁鏈估算模型,針對不同參數的變化對磁鏈估算影響不同的問題,利用頻率響應函數法(FRF)進行定量分析,獲取影響較大的參數。其次,基于模型參考自適應辨識方法(MRAS),將辨識初始值作為補償項,以克服電機轉速和負載轉矩變化對辨識結果的影響,建立轉子時間常數Tr的在線辨識改進模型,提出考慮參數辨識的風電機組變槳控制策略。最后,對風電機組變槳電機及其控制系統進行動態仿真和樣機實驗。結果表明,Tr對變槳電機控制系統的影響較為明顯,考慮在線辨識后的改進模型及控制策略具有更好的魯棒性。

關鍵詞:風力發電;變槳電機;矢量控制;參數辨識;模型參考自適應

中圖分類號:TM 34

文獻標志碼:A

文章編號:1007-449X(2019)07-0009-10

Abstract:Uncertainties of parameters and load in pitch motor may lead to low robustness of pitch control system.Therefore,influence of motor parameters on control system was analyzed and an improved parameter identification method was presented.Firstly,based on the complex flux estimator in pitch control system,impact of motor parameters on the flux estimator were analyzed using frequency response functions (FRF),which could derive the key parameter.Secondly,based on the model reference adaptive system (MRAS),the initial value of identification was regarded as the compensation term to weaken the effect of variable motor speeds and load torques on identification result,then the improved online identification of the rotor time constant Tr and pitch control strategy considering parameter identification in wind turbine system were presented.Finally,simulation and prototype test of pitch motor control system were carried out.The result validates that Tr has a prominent impact on the pitch control system,and the improved identification model and control strategy considering online identification would be more robust.

Keywords:wind power generation; pitch motor; vector control; parameter identification; model reference adaptive system

0?引?言

近年來,變槳距風力發電機組由于能夠有效利用風能并減少風力沖擊,提高系統運行的可靠性,已成為目前的主流機型[1-3]。電動變槳系統可采用直流、永磁同步、感應電機等實現驅動,其中直流電機由于換向器、電刷等結構帶來的低可靠性限制了其應用[4];永磁同步電機具有體積小、效率高、調速性能好等優點,但其成本高,且高溫振動條件下永磁體存在退磁風險,可靠性較差;而三相感應電機因具有成本低、結構簡單可靠、轉矩大等優點[5],目前仍然廣泛應用于風電機組的變槳驅動中[6]。然而,風湍流、風剪切、塔影效應、偏航偏差等因素會使風電機組的葉輪產生不均衡載荷[7],且變槳系統在輪轂中隨著葉片的旋轉而工作[8],勢必會導致變槳電機負載轉矩的不確定性,此外,變槳電機在運行過程中由于內部溫度變化和磁路飽和等因素的影響,參數也存在不確定性。轉矩和參數的不確定性將直接影響風電機組變槳控制的魯棒性。因此,研究變槳電機參數辨識,尋求適應變槳電機參數變化和負載轉矩不確定性的變槳控制策略對提高風電機組運行可靠性和穩定性具有重要的現實意義。

針對系統參數變化以及風速隨機性等導致風電變槳控制系統魯棒性差的問題,國內外學者提出了較多的控制策略[9-11]。文獻[10]基于滑模控制理論,提出一種降低風電機組載荷的多目標變槳距控制策略,可有效地抑制機械振動并對風電機組動力學模型的參數變化具有較好的魯棒性。文獻[11]提出了一種基于模糊邏輯的變槳控制策略并與PID控制方法進行對比,結果表明模糊邏輯控制可以適應變槳系統的非線性,使得風電機組的輸出功率保持平穩。然而,從現有文獻對變槳距風電機組控制策略的研究來看,很少考慮到變槳電機參數變化對其控制性能的影響。針對感應電機轉子磁場定向矢量控制中參數變化的問題,已有眾多的辨識方法,其中模型參考自適應控制(model reference adaptive system,MRAS)因其易于實現、穩定可靠而得到了廣泛地應用[12-15]。文獻[12]證明了基于無功功率的轉子電阻與轉速MRAS在線辨識模型的穩定性,并建立了小信號模型分析定、轉子電阻對磁場定向的敏感性,但是未涉及如轉子時間常數等其他參數,也沒有考慮負載轉矩波動對辨識準確性的影響。文獻[14]利用Popov超穩定性理論,分析了基于無功功率的MRAS轉子時間常數在線辨識方法的穩定性,并研究了該辨識方法對所涉及電機參數的敏感性,但沒有考慮到電機轉速和負載轉矩變化對辨識準確性的影響。因此,現有文獻涉及的參數辨識方法是否適用于變槳電機的運行工況,還有待進一步研究,此外,也缺乏針對變槳電機參數辨識與變槳控制相結合的魯棒控制策略的研究。

基于此,考慮風電變槳電機參數變化和負載不確定性的影響,本文提出一種適應變槳電機運行工況的參數辨識改進方法及相應的變槳控制策略。首先,利用頻率響應函數法定量分析變槳電機不同參數的變化對控制系統中復合磁鏈估算的影響,獲取影響較為明顯的轉子時間常數Tr。其次,基于模型參考自適應方法,考慮電機轉速和負載轉矩變化的情況,以辨識初始值作為補償項,建立Tr在線辨識的改進模型并提出考慮參數辨識的風電機組變槳控制策略。最后,對變槳電機控制系統進行動態仿真和樣機實驗,驗證改進模型和控制策略的有效性。

1?變槳電機參數變化對磁鏈估算的影響

1.1?變槳電機控制系統簡介

以感應電機為例的風電變槳控制系統如圖1所示,其中控制環包括位置環、速度環、磁鏈環和電流環,脈沖占空比通過空間矢量脈寬調制(space vector pulse width modulation,SVPWM)計算得到,圖中符號的下標ref代表各物理量的參考值,d、q分別表示dq坐標系下的分量,α、β為αβ坐標系下的分量。符號γ表示槳距角,ωr為轉子的電角速度,|ψr |、θ分別為轉子磁鏈的幅值與相位,is、us分別為定子電流、電壓,Ta、Tb、Tc分別為逆變器abc相上橋臂開關器件的占空比,Vdc為直流母線電壓。

由圖1可知,磁鏈估算模塊輸出的|ψr|用于反饋跟蹤,θ用于坐標變換,因此磁鏈估算的準確性對矢量控制尤為重要。本文采用文獻[16]提出的磁鏈估算方法,利用電流模型與電壓模型所估算的磁鏈之差生成補償項,以保證在寬調速范圍內轉子磁鏈的估算精度。圖1中磁鏈估算模塊的原理框圖如圖2所示,其中上標i和v分別表示電流模型和電壓模型中的物理量,|ψs|為定子磁鏈的幅值,ucom_αs、ucom_βs分別為生成的定子補償電壓在αβ坐標系下的分量,Rs為定子電阻。

由式(5)可知,Lm僅對磁鏈的幅值產生影響,且不隨ωs變化,僅與Lm變化的比值成正比;而Tr對幅值和相位均有影響,并與ωs有關。

頻率響應函數可以直觀地表征系統在給定頻率下的穩態輸出與輸入的幅值和相位關系,適合于分析參數變化對磁鏈估算的影響,文獻[17]也利用此方法對磁鏈觀測器的參數靈敏度做出了分析,但僅僅考慮了參數變為1.5倍之后的影響。為了定量分析單一參數變化與磁鏈估算的關系,獲取影響較大的參數,將Tr和Lm的實際值分別變為估算值的2~5倍,則實際磁鏈的幅值、相位與觀測值之差關于參數的變化倍數、轉差率s的關系如圖3所示,其中轉差率s=ωs/ωe。

由圖3可知,在參數變化相同倍數的情況下,當s>0.1時,Tr對幅值的影響與Lm相當,s在0.01附近時Tr對相位的影響最大。由于在達到給定槳距角之前,變槳電機會在高速狀態下運行,s一般在0.02左右;而達到槳距角之后,為了平衡負載轉矩,s為1。由此可見,就變槳電機運行工況而言,轉差率一般在0.02~1的范圍內,即高速低速狀態均會出現,且大部分為低速運行的工況,因此Tr對電流模型所觀測出的磁鏈影響較為明顯。

1.3?參數變化對電壓模型磁鏈估算的影響

由式(12)、式(13)可知,Tr、Lm、Rs、Ls和Lr等參數對電壓模型所觀測出的磁鏈幅值與相位均有影響。但定轉子漏感相對Lm而言較小,可忽略其對Ls和Lr的影響,因此在考慮磁路飽和對電機電感的影響時,可認為Ls、Lr和Lm的變化基本相等[18],僅分析Lm對磁鏈觀測的影響即可。基于此,當實際的Tr、Lm和Rs分別變為估算值的2~5倍時,實際磁鏈的幅值、相位與觀測值之差如圖4所示。

對比圖4(a)、圖4(b)可知,Tr與Rs對幅值的影響比Lm更大,而在圖4(a)中,隨著轉差率s的增大,即低速情況下Tr的影響更為顯著;就相位而言,由圖4(c)可知,僅在s=1附近的低速情況下,Lm的影響比Tr更為顯著,對比圖4(c)、圖4(d)可以看出,低速時,當參數變化超過4倍,Rs的影響最大,而在s=0.1左右,Tr對相位的影響較Lm和Rs更大。

在圖2所示的變槳電機磁鏈估算模型中,選擇適當的PI參數可以使電流模型在低速時占主導,電壓模型則主要在高速時起作用。因此,綜合考慮變槳電機全速范圍內的運行工況,對變槳電機磁鏈估算影響較為明顯的參數主要是Tr。本文以Tr的在線辨識模型及控制策略開展研究。

2?轉子時間常數在線辨識的改進模型

2.1?考慮初始值作為補償項的辨識改進模型

利用Popov超穩定性定律推導自適應律,與轉速估計[19]類似,以式(7)所示電壓模型為參考模型,式(2)所示電流模型為可調模型,上標*表示標幺值,~表示可調模型,取比例加積分的自適應律為

在常規的MRAS模型中,辨識初始值一般取為額定值。由于轉子磁場定向實現了解耦控制[5],隨著變槳電機負載轉矩的變化,參考模型計算的磁鏈會隨之改變,轉速的變化也會使電流模型的計算產生偏差,由此可知,若初始值仍保持不變,則會導致由式(14)辨識出的Tr變化。因此,在變槳電機轉速與負載轉矩不確定的情況下,僅將T~r(0)當作常數可能使辨識結果存在較大誤差。為此,本文考慮將T~r(0)作為補償項。

在兩相同步旋轉的dq軸系下有[5]

由于實際的控制系統均為離散系統,利用k-1時刻的數據計算k時刻的Tr,式(15)可以表示為

2.2?不同轉速和負載轉矩下的仿真對比驗證

為了驗證Tr辨識改進模型的有效性,本文對變槳電機在不同轉速、不同負載轉矩下的辨識效果進行仿真,并將結果與未改進的方法進行比較,仿真中變槳電機參數如表1所示,上一節分析參數變化對磁鏈估算的影響也是基于表1中的參數。

仿真中分別給定轉速為0.97pu、0.5pu和0.2pu,負載轉矩為額定值36 N·m,t=1 s時MRAS辨識模型開始作用,變槳電機Tr初始值為0.155 s,t=3 s時變為0.310 s,改進前后的辨識結果對比如圖6所示。

對比圖6(a)~圖6 (c)可知,高速時改進前后的結果差別不大,但隨著轉速的降低,改進后的辨識精度有所提升。改進前后辨識結果的相對誤差如表2所示。

由表2中的誤差對比可知,改進后的辨識模型能夠更好地適應電機參數和轉速的變化。

考慮變槳電機的負載具有不確定性,在轉速為0.2pu的情況下模擬負載轉矩的變化,t=2 s時轉矩突增為50 N·m,4 s時變為額定值36 N·m,改進前后的辨識結果如圖7所示。

由圖7可知,改進辨識模型可較好地抑制Tr隨負載變化的波動,轉矩為50 N·m時改進前后的辨識誤差對比如表3所示。

由表3可知,轉矩變化后,改進方法仍然能夠準確辨識出實際的Tr。因此,改進的辨識模型可以更好地適應變槳電機負載的不確定性。

2.3?樣機實驗驗證

圖8所示為本文所用實驗平臺,其中控制算法采用TI公司的TMS320F28335 DSP芯片實現,PWM開關頻率為10 kHz,電流環、速度環的控制周期分別為100 μs和1 ms。以直流電機作為負載,可通過調節其勵磁以及電樞電流來改變轉矩。由于實驗條件限制,采用的感應電機為小功率電機,具體參數見表4,其中轉子時間常數的額定值約為0.107 s。需要說明的是,前面所提的參數辨識改進方法并非針對特定的感應電機而言,因此,采用小功率電機并不影響參數辨識及控制的效果,實驗結果也可與仿真相對應,進而驗證理論的正確性。

由于變槳電機的Tr對控制系統的影響主要體現在低速段。因此在實驗中給定轉速0.2pu,提取1.6 s時間內的實驗數據,其中0~0.8 s時負載為額定值,0.8~1.6 s時負載變為1.5倍額定值,得出加入改進、未改進辨識模型以及未加辨識3種情況下的效果如圖9所示。

由圖9(a)可知,未加辨識環節時,轉速存在較大的振蕩,且轉矩增大后速度無法控制在0.2pu;加入未改進的辨識模型,負載轉矩增大后轉速控制效果仍不理想;改進的辨識模型可以使轉速控制較好地適應負載轉矩變化的情況。加入改進、未改進辨識模型以及未加辨識3種情況下轉速控制的相對誤差如表5所示,由此可知改進的辨識模型可以適應負載轉矩的變化,使轉速控制達到較好的效果。

圖9(b)表明在感應電機的運行過程中,溫度升高、磁路飽和等現象會使Tr相比額定值0.107 s發生一定的變化,改進的辨識模型能更準確地辨識出Tr的變化,使得整個變槳系統在受到參數和負載轉矩變化后仍然能達到良好的控制效果。

3?考慮參數辨識的變槳控制策略

3.1?考慮轉子時間常數辨識的變槳控制策略

將圖5所示的改進辨識模型加入到圖1的變槳電機控制系統中,構成如圖10所示的考慮參數辨識后的風電變槳電機控制策略。其中辨識模塊的輸入為變槳電機定子電壓uαs、uβs,定子電流iαs、iβs、ids、iqs,同步轉速ωe以及轉子轉速ωr;輸出的Tr作為磁鏈估算模塊的輸入,實時修正其中的參數Tr,使得估算出的磁鏈更加精確,有利于提高整個變槳系統的魯棒性。

3.2?電機參數變化對變槳性能影響仿真

考慮到變槳電機負載轉矩的不確定性,仿真初始在額定負載轉矩下工作,4 s時突增為50 N·m,仿真時間設定為8 s。若不加入Tr辨識環節,在系統參數均為初始值和變槳電機Tr改為0.310 s 2種情況下,仿真得到變槳系統位置環與速度環的控制效果對比如圖11所示。

從圖11可以看出,將Tr改為初始值2倍后,電機在低速時帶動50 N·m的負載轉矩會出現轉速振蕩,導致位置環的控制出現一定波動與偏差,驗證了感應電機轉子磁場定向矢量控制系統對Tr的變化具有敏感性。加入未改進Tr辨識環節,由圖6、圖7的仿真結果可知,當Tr和負載轉矩變化后,辨識結果仍然不準確,其控制效果與圖11類似。而加入改進Tr辨識環節后,控制和辨識效果如圖12所示。

由圖12可知,變槳電機的Tr改變后,改進的辨識環節能夠在電機轉速及負載變化時精確得出Tr的實際值,以此作為磁鏈估算模塊的輸入,實時地修正輸出磁鏈的幅值和相角,可以使風電變槳控制系統對電機Tr變化的敏感度降低,魯棒性增強。

4?結?論

針對風電變槳電機參數變化和負載不確定性導致控制系統魯棒性差的問題,提出了一種變槳電機參數辨識改進模型及其變槳控制策略。通過仿真和樣機實驗等驗證了本文所提方法的有效性。主要結論有:

1)對磁鏈估算模型的理論分析與變槳電機控制策略的仿真分析表明,影響變槳系統性能的主要參數是變槳電機的轉子時間常數Tr;

2)無論是變槳電機參數、轉速還是負載轉矩的變化,本文提出的Tr在線辨識改進模型都能夠較準確地跟蹤實際參數,使得電機的轉速控制更有效。

3)考慮參數辨識后的變槳電機控制策略能夠在一定程度上抑制變槳電機在低速時轉速的振蕩,提高整個變槳系統控制的魯棒性。

參 考 文 獻:

[1]?李輝,楊超,趙斌,等.風電機組電動變槳系統建模及運行特性評估與測試[J].電力系統自動化,2013,37(11):20.

LI Hui,YANG Chao,ZHAO Bin,et al.Modeling and operation performance evaluation and test of electric pitch system in wind turbine generator system[J].Automation of Electric Power System,2013,37(11):20.

[2]?史運濤,侯彥嬌,孫德輝,等.風能轉換系統隨機建模與H∞容錯控制[J].電機與控制學報,2015,19(3):100.

SHI Yuntao,HOU Yanjiao,SUN Dehui,et al.Stochastic modelling and H∞ fault tolerance control of WECS[J].Electric Machines and Control,2015,19(3):100.

[3]?楊俊友,王海鑫,井艷軍,等.并網型風電機組模擬控制策略[J].電機與控制學報,2016,20(3):43.

YANG Junyou,WANG Haixin,JING Yanjun,et al.Imitation control strategy of gridconnected wind power system[J].Electric Machines and Control,2016,20(3):43.

[4]?周俊宇,鄧文浪.永磁同步電機風力發電變槳控制系統設計[J].中南林業科技大學學報,2011,31(9):165.

ZHOU Junyu,DENG Wenlang.Design of pitch controlled system of wind turbines based on PMSM[J].Journal of Central South University of Forest & Technology,2011,31(9):165.

[5]?高景德,王祥珩,李發海.交流電機及其系統的分析[M].北京: 清華大學出版社,2005.

[6]?李輝,楊超,趙斌,等.考慮變槳驅動電機特性的風電機組運行性能仿真[J].電力系統保護與控制,2011,39(18):26.

LI Hui,YANG Chao,ZHAO Bin,et al.Simulation on the operational performances of wind turbine generator system considering the variable pitch drive motor characteristics[J].Power System Protection and Control,2011,39(18):26.

[7]?BURTON T.風能技術[M].武鑫,譯.北京:科學出版社,2014.

[8]?李輝,楊超,李學偉,等.風機電動變槳系統狀態特征參量挖掘及異常識別[J].中國電機工程學報,2014,34(12):1922.

LI Hui,YANG Chao,LI Xuewei,et al.Conditions characteristic parameters mining and outlier identification for electric pitch system of wind turbine[J].Proceedings of the CSEE,2014,34(12):1922.

[9]?HWAS A,KATEBI R.Wind turbine control using PI pitch angle controller[J].IFAC Proceedings Volumes,2012,45(3):241.

[10]?肖帥,楊耕,耿華.抑制載荷的大型風電機組滑模變槳距控制[J].電工技術學報,2013,28(7):145.

XIAO Shuai,YANG Geng,GENG Hua.Slidingmode pitch control strategy for large wind turbines to reduce loads[J].Transactions of China Electrotechnical Society,2013,28(7):145.

[11]?VAN T L,NGUYEN T H,LEE D C.Advanced pitch angle control based on fuzzy logic for variablespeed wind turbine systems[J].IEEE Transactions on Energy Conversion,2015,30(2):578.

[12]?MAITI S,CHAKRABORTY C,HORI Y,et al.Model reference adaptive controllerbased rotor resistance and speed estimation techniques for vector controlled induction motor drive utilizing reactive power[J].IEEE Transactions on Industrial Electronics,2008,55(2):594.

[13]?徐占國,邵誠,馮冬菊.基于模型參考自適應的感應電機勵磁互感在線辨識新方法[J].中國電機工程學報,2010,30(3):71.

XU Zhanguo,SHAO Cheng,FENG Dongju.Online identification of induction motor mutual inductance based on model reference adaptive system[J].Proceedings of the CSEE,2010,30(3):71.

[14]?王高林,楊榮峰,張家皖,等.一種感應電機轉子時間常數MRAS的在線辨識方法[J].電工技術學報,2012,27(4):48.

WANG Gaolin,YANG Rongfeng,ZHANG Jiawan,et al.Rotor time constant online estimation of induction motors based on MRAS[J].Transactions of China Electrotechnical Society,2012,27(4):48.

[15]?CAO Pengpeng,ZHANG Xing,YANG Shuying.A unifiedmodelbased analysis of MRAS for online rotor time constant estimation in an induction motor drive[J].IEEE Transactions on Industrial Electronics,2017,64(6):4361.

[16]?LASCU C,BOLDEA I,BLAABJERG F.A modified direct torque control for induction motor sensorless drive [J].IEEE Transactions on Industry Applications,2000,36(1):124.

[17]?JANSEN P L,LORENZ R D.A physically insightful approach to the design and accuracy assessment of flux observers for field oriented induction machine drives [J].IEEE Transactions on Industry Applications,2002,30(1):103.

[18]?徐小增,李葉松,秦憶.矢量控制感應電機參數變化的影響研究[J].華中科技大學學報(自然科學版),2003,31(7):44.

XU Xiaozeng,LI Yesong,QIN Yi.Effect of parameters variety on vectorcontrolled induction motor [J].Journal of Huazhong University of Science and Technology (Nature Science Edition),2003,31(7):44.

[19]?LEI Wang,LI Yinghui,LEI Xiao.Sliding mode variablestructure MRAS speed identification for induction motor direct torque control system[C]//2009 ISECS International Colloquium on Computing,Communication,Control,and Management,August 8-9,2009,Sanya,China.2009: 9-62.

(編輯:邱赫男)

猜你喜歡
風力發電
關于“棄風”現象研究的文獻綜述
東方教育(2016年13期)2017-01-12 22:44:43
農用風力發電機葉片設計
風力發電系統故障診斷技術淺析
東汽FD77雙饋式變速恒頻風力發電機的控制策略
電池儲能技術在風力發電系統中的應用
中國市場(2016年32期)2016-12-06 11:06:23
風力發電發展現狀及前景淺析
價值工程(2016年31期)2016-12-03 00:07:43
垂直軸風力發電機葉片建模與研究農用
呼倫貝爾風電場對草原生態環境潛在影響初探
科技視界(2016年16期)2016-06-29 17:02:02
我國風力發電現狀與展望
商(2016年12期)2016-05-09 10:38:04
國外新能源發電技術開發及應用的方式對我國的啟迪
科技視界(2016年8期)2016-04-05 08:56:22
主站蜘蛛池模板: 一级黄色网站在线免费看| 91外围女在线观看| 欧美精品成人一区二区在线观看| V一区无码内射国产| 国内精品免费| 国产成人久视频免费| 亚洲婷婷丁香| 日本手机在线视频| 亚洲人成日本在线观看| 91麻豆精品国产高清在线| 国产小视频a在线观看| 中文成人在线| 免费高清自慰一区二区三区| 国产精品爆乳99久久| 亚洲综合片| 91精品视频在线播放| 国产高清免费午夜在线视频| 亚洲精选高清无码| 国产成在线观看免费视频| 国产三级毛片| 欧美三级视频网站| 亚洲综合久久成人AV| 亚洲第一区精品日韩在线播放| 精品少妇人妻av无码久久| 国产美女丝袜高潮| 久久国语对白| 日本午夜精品一本在线观看 | 精品无码一区二区三区电影| 国内精品自在自线视频香蕉| 日本不卡在线| 孕妇高潮太爽了在线观看免费| 久久久久久久蜜桃| 中文成人在线| 国产www网站| 操操操综合网| 亚洲一区二区三区在线视频| 亚洲AV无码久久精品色欲| 国产成人精品一区二区不卡| 久久久久免费精品国产| 一级毛片在线免费视频| 国产真实乱子伦精品视手机观看 | 国产成在线观看免费视频 | 五月天福利视频| 日韩视频免费| 久久99热这里只有精品免费看| 黄色网址免费在线| 欧美精品H在线播放| 国产又爽又黄无遮挡免费观看| 国产成人亚洲无吗淙合青草| 国内a级毛片| 成人日韩精品| 天天综合天天综合| 亚洲中文字幕在线一区播放| 久操中文在线| 一级一级一片免费| 亚洲AV人人澡人人双人| 久久综合一个色综合网| 国产一级二级三级毛片| 亚洲国产成人精品无码区性色| a亚洲天堂| 成·人免费午夜无码视频在线观看| 国产91熟女高潮一区二区| 亚洲男人天堂久久| 欧美亚洲欧美| 亚洲成人免费看| 黄色成年视频| 午夜欧美在线| 3344在线观看无码| 成人午夜视频免费看欧美| 国内精自视频品线一二区| 日韩高清成人| 国产精彩视频在线观看| 色悠久久综合| 亚欧美国产综合| 欧美成人手机在线视频| 高清亚洲欧美在线看| 欧美国产成人在线| 久久黄色毛片| 亚洲成在线观看| 久久这里只有精品免费| 91色在线视频| 青青草久久伊人|