亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        PMSG風(fēng)力發(fā)電系統(tǒng)轉(zhuǎn)速估計(jì)算法的研究

        2016-10-14 15:14:03楊克立王克軍
        關(guān)鍵詞:風(fēng)速

        程 輝,楊克立,王克軍,李 娜

        ?

        PMSG風(fēng)力發(fā)電系統(tǒng)轉(zhuǎn)速估計(jì)算法的研究

        程 輝1,楊克立2,王克軍1,李 娜1

        (1.河南工程學(xué)院電氣信息工程學(xué)院,河南 鄭州 451192; 2.中原工學(xué)院工業(yè)訓(xùn)練中心,河南 鄭州 450007)

        提出了一種無(wú)傳感器直驅(qū)式永磁同步風(fēng)力發(fā)電機(jī)控制系統(tǒng)的轉(zhuǎn)子速度估算法。該算法可由簡(jiǎn)單的定子磁通方程推導(dǎo),且只需檢測(cè)定子電壓和電流。定子磁通估算通過(guò)基于定子電壓模型的可編程低通濾波器而設(shè)計(jì),建立了基于遞歸最小二乘法的電機(jī)轉(zhuǎn)速辨識(shí)模型,并將實(shí)時(shí)辨識(shí)的電機(jī)轉(zhuǎn)速用于轉(zhuǎn)子位置觀測(cè)和轉(zhuǎn)子磁場(chǎng)定向矢量控制。仿真結(jié)果表明,實(shí)現(xiàn)了最大功率跟蹤控制,驗(yàn)證了控制算法的有效性和可行性。

        永磁同步發(fā)電機(jī);轉(zhuǎn)速估計(jì);風(fēng)力發(fā)電系統(tǒng);無(wú)傳感器;最大功率跟蹤控制

        0 引言

        風(fēng)力發(fā)電系統(tǒng)主要有恒速恒頻風(fēng)力發(fā)電系統(tǒng)和變速恒頻風(fēng)力發(fā)電系統(tǒng)兩大類,后者使用永磁同步發(fā)電機(jī),可在很寬的風(fēng)速變化范圍內(nèi)實(shí)現(xiàn)最大風(fēng)能捕獲,具有高功率密度,且無(wú)需齒輪箱和滑環(huán)及勵(lì)磁裝置[1],結(jié)構(gòu)簡(jiǎn)單,可靠性高,是未來(lái)風(fēng)力發(fā)電技術(shù)發(fā)展的主要趨勢(shì)之一。

        風(fēng)能的不確定性以及風(fēng)力機(jī)自身特性會(huì)使風(fēng)力發(fā)電機(jī)輸出功率隨風(fēng)速變化而波動(dòng)。為了提高風(fēng)能的利用效率,最大程度捕獲風(fēng)能,需使風(fēng)力機(jī)在任意風(fēng)速時(shí)都能運(yùn)行在最佳轉(zhuǎn)速,為此需要進(jìn)行高精度的風(fēng)速檢測(cè)和轉(zhuǎn)速檢測(cè),但也降低了系統(tǒng)運(yùn)行的可靠性,另外,風(fēng)力機(jī)周?chē)鷼饬魇苋~片擾動(dòng)較大,難以準(zhǔn)確測(cè)量當(dāng)前有效風(fēng)速[2]。與此同時(shí),定子電流和定子電壓的準(zhǔn)確測(cè)量可用于轉(zhuǎn)速估算,從而無(wú)需進(jìn)行轉(zhuǎn)速測(cè)量。目前,國(guó)內(nèi)外學(xué)者對(duì)直驅(qū)式PMSG無(wú)速度傳感器矢量控制進(jìn)行了大量的研究。對(duì)永磁同步電機(jī)位置和轉(zhuǎn)速估算較為成熟的方法有高頻信號(hào)注入法和基于反電勢(shì)估算法。前者利用電機(jī)的凸極效應(yīng)具有較強(qiáng)的魯棒性,特別適合于零速啟動(dòng)和低速運(yùn)行,但是只適用于具有凸極效應(yīng)的電機(jī),并且注入的高頻信號(hào)會(huì)使電機(jī)損耗、轉(zhuǎn)矩波動(dòng)及噪聲變大,電機(jī)的動(dòng)態(tài)性能變差;后者是從測(cè)量電機(jī)的發(fā)電動(dòng)勢(shì)入手,算法實(shí)現(xiàn)較為簡(jiǎn)單,但由于低速下電機(jī)的反電動(dòng)勢(shì)值較小,且易受系統(tǒng)測(cè)量誤差的影響,效果不夠理想[3-4]。此外,文獻(xiàn)[5]提出了引入軟件鎖相環(huán)的控制策略,通過(guò)鎖相環(huán)得到風(fēng)力機(jī)的轉(zhuǎn)速,但由于發(fā)電機(jī)輸出電壓信號(hào)中夾雜有部分諧波分量等影響相位鎖定的精度。文獻(xiàn)[6-7]提出擴(kuò)展卡爾曼濾波算法,但由于計(jì)算量很大需要利用高精度的數(shù)字信號(hào)處理器來(lái)實(shí)現(xiàn)實(shí)時(shí)控制,工程上不宜實(shí)現(xiàn)。除了采用基于數(shù)學(xué)模型的方法以外,還可以利用基于人工智能理論進(jìn)行轉(zhuǎn)子速度的估算[8-10]。

        本文針對(duì)無(wú)速度傳感器的PMSG變速風(fēng)能發(fā)電系統(tǒng)提出了一種轉(zhuǎn)速估計(jì)矢量控制算法,通過(guò)簡(jiǎn)單的定子磁通模型以進(jìn)行轉(zhuǎn)速估算,無(wú)需測(cè)量裝置就可以得到風(fēng)力機(jī)轉(zhuǎn)速,通過(guò)控制直驅(qū)永磁風(fēng)力發(fā)電機(jī)的電磁功率來(lái)間接控制風(fēng)力發(fā)電機(jī)轉(zhuǎn)速以追蹤風(fēng)力機(jī)最佳功率曲線實(shí)現(xiàn)對(duì)最大風(fēng)能的捕獲。文中,轉(zhuǎn)速采用遞歸最小二乘法的參數(shù)估計(jì)算法進(jìn)行實(shí)時(shí)在線辨識(shí),并將實(shí)時(shí)辨識(shí)的轉(zhuǎn)速應(yīng)用于PMSG轉(zhuǎn)子位置估計(jì)和轉(zhuǎn)子磁場(chǎng)定向矢量控制。定子磁通估算通過(guò)基于定子電壓模型的可編程低通濾波器PLPF設(shè)計(jì),以提高定子磁通的估算性能。

        1 風(fēng)力機(jī)特性

        由空氣動(dòng)力學(xué)可知,風(fēng)力機(jī)從風(fēng)能中捕獲的機(jī)械功率為[11]

        (2)

        風(fēng)力機(jī)的功率-轉(zhuǎn)速特性曲線如圖1所示。

        圖1 風(fēng)力機(jī)功率-轉(zhuǎn)速特性曲線

        2 無(wú)傳感器PMSG 變速風(fēng)能發(fā)電控制系統(tǒng)

        在兩相靜止坐標(biāo)系下,表面式PMSM的數(shù)學(xué)模型[14]可以表示為

        (4)

        定子磁通(5)

        對(duì)式(5)求微分并聯(lián)立式(4)可得:

        (7)

        由式(7)可知,在定子電流很容易測(cè)量獲得的基礎(chǔ)上,如果定子磁通已知的話轉(zhuǎn)子轉(zhuǎn)速就可以估算出來(lái)。

        2.1 定子磁通估計(jì)

        2.2 轉(zhuǎn)子轉(zhuǎn)速估計(jì)

        重新整理式(7)可得,

        本文所辨識(shí)的參數(shù)是時(shí)變的電機(jī)轉(zhuǎn)子轉(zhuǎn)速,根據(jù)極值定理可以推導(dǎo)出引入遺忘因子的最小二乘法的數(shù)學(xué)表達(dá)式為

        (10)

        (12)

        固定值矩陣的跡

        (14)

        簡(jiǎn)化后的恒增益遞歸最小二乘估算法表達(dá)式為

        (16)

        2.3 機(jī)側(cè)變流器控制策略

        根據(jù)轉(zhuǎn)子磁場(chǎng)定向可得到同步旋轉(zhuǎn)坐標(biāo)系下PMSG的數(shù)學(xué)模型為

        (19)

        機(jī)側(cè)變流器的d軸變量和q軸變量間存在耦合,給控制器的設(shè)計(jì)帶來(lái)困難,為此采用前饋解耦控制策略。

        當(dāng)PMSM穩(wěn)定運(yùn)行時(shí),且忽略定子電阻壓降的情況下,電壓方程可以簡(jiǎn)化為

        (21)

        機(jī)側(cè)變流器的矢量控制如圖2所示。采用矢量雙閉環(huán)控制技術(shù)。外環(huán)為轉(zhuǎn)速環(huán),其速度參考值通過(guò)最大功率點(diǎn)跟蹤算法給出,調(diào)節(jié)發(fā)電機(jī)的電磁轉(zhuǎn)矩,使得實(shí)際轉(zhuǎn)速跟蹤此參考轉(zhuǎn)速即可捕獲最大風(fēng)能。參考轉(zhuǎn)速與采用估計(jì)算法估算的轉(zhuǎn)速進(jìn)行比較,將其偏差值送入PI調(diào)節(jié)器即可輸出有功電流的參考值;內(nèi)環(huán)分別實(shí)現(xiàn)d、q電流的閉環(huán)控制,誤差信號(hào)經(jīng)PI調(diào)節(jié)后與各自的解耦補(bǔ)償項(xiàng)相加得到機(jī)側(cè)變流器參考電壓,以此進(jìn)行SVPWM調(diào)制。按照單位功率因數(shù)控制設(shè)定d軸參考電流。

        圖2 機(jī)側(cè)變流器控制系統(tǒng)

        2.4 網(wǎng)側(cè)變流器控制策略

        在同步旋轉(zhuǎn)dq坐標(biāo)系下,網(wǎng)側(cè)變換器的數(shù)學(xué)模型為

        (23)

        式中:和分別是網(wǎng)側(cè)進(jìn)線電抗器的電感值和電阻值;和分別是網(wǎng)側(cè)變流器輸入三相交流電壓的d、q軸分量;和分別為電網(wǎng)側(cè)三相交流電壓的d、q軸分量。同樣,網(wǎng)側(cè)變流器的d軸變量和q軸變量之間也存在耦合,需要采用前饋解耦控制策略??刂瓶驁D如圖3所示。

        圖3 網(wǎng)側(cè)變流器控制系統(tǒng)

        3 仿真結(jié)果

        為了驗(yàn)證文章所提出的控制策略的可行性,利用Matlab/Simulink進(jìn)行了仿真。建立了無(wú)風(fēng)速無(wú)轉(zhuǎn)子轉(zhuǎn)速傳感器的PMSG直驅(qū)變速風(fēng)力發(fā)電系統(tǒng),仿真模型如圖4所示。圖中,機(jī)側(cè)以獲得額定磁通,網(wǎng)側(cè)以獲取單位功率因數(shù)。PMSG仿真參數(shù):直流母線電壓400 V,電樞繞組交直軸電感為4 mH,每相繞組電阻為0.9,電機(jī)極對(duì)數(shù)為18,電機(jī)額定電壓為200 V,額定電流11.6 A,額定功率2 kW,額定轉(zhuǎn)速200 r/min。

        圖4 PMSG風(fēng)能發(fā)電系統(tǒng)

        仿真中用到的風(fēng)速曲線有兩種。一種曲線如圖5(a)所示,采用理想風(fēng),最高風(fēng)速為9 m/s,最低風(fēng)速為7 m/s,采用無(wú)傳感器轉(zhuǎn)子估計(jì)算法控制系統(tǒng)的轉(zhuǎn)子速度、功率因數(shù)及有功功率輸出分別如圖5(b)(c)(d)所示。為方便比較起見(jiàn),圖6(a)給出了對(duì)應(yīng)的轉(zhuǎn)子速度及估計(jì)的轉(zhuǎn)子速度仿真曲線,圖中點(diǎn)劃線表示為估計(jì)的轉(zhuǎn)子速度,實(shí)線為實(shí)際轉(zhuǎn)子速度;圖6(b)為兩者速度誤差曲線。

        圖5 理想風(fēng)速下無(wú)傳感器PMSG風(fēng)能發(fā)電系統(tǒng)仿真結(jié)果

        圖6 轉(zhuǎn)速和估計(jì)轉(zhuǎn)速波形

        從圖中可以看出,估計(jì)風(fēng)速和實(shí)際風(fēng)速誤差的最大值為0.64 rad/s,控制的效果非常理想。

        如圖7(a)所示為自然風(fēng)速曲線,采用無(wú)傳感器轉(zhuǎn)子估計(jì)算法控制系統(tǒng)的轉(zhuǎn)子速度、功率因數(shù)及有功功率輸出分別如圖7(b)(c)(d)所示。為方便比較起見(jiàn),圖8(a)給出了對(duì)應(yīng)的轉(zhuǎn)子速度及估計(jì)的轉(zhuǎn)子速度仿真曲線,圖中點(diǎn)劃線表示為估計(jì)的轉(zhuǎn)子速度,實(shí)線為實(shí)際轉(zhuǎn)子速度;圖8(b)為兩者速度誤差曲線。

        圖7 自然風(fēng)速下無(wú)傳感器PMSG風(fēng)能發(fā)電系統(tǒng)仿真結(jié)果

        圖8 轉(zhuǎn)速和估計(jì)轉(zhuǎn)速波形

        從圖8可以看出,通過(guò)定子磁通估算轉(zhuǎn)子速度跟蹤風(fēng)力轉(zhuǎn)速效果良好,估計(jì)風(fēng)速和實(shí)際風(fēng)速誤差的最大值為0.83 rad/s。結(jié)果表明,本文提出的轉(zhuǎn)子速度估計(jì)算法準(zhǔn)確度較高,控制效果較好。

        4 結(jié)論

        本文研究了PMSG風(fēng)力發(fā)電系統(tǒng)的轉(zhuǎn)子速度估算法,無(wú)需測(cè)量風(fēng)速和風(fēng)力機(jī)轉(zhuǎn)速就可以實(shí)現(xiàn)風(fēng)能的最大功率跟蹤。該算法通過(guò)簡(jiǎn)單的定子磁通方程推導(dǎo),定子磁通估算通過(guò)基于定子電壓模型的可編程低通濾波器設(shè)計(jì),利用遺忘因子的遞推最小二乘算法對(duì)轉(zhuǎn)速進(jìn)行實(shí)時(shí)辨識(shí)。本文建立了系統(tǒng)的數(shù)學(xué)模型,詳細(xì)進(jìn)行了系統(tǒng)的設(shè)計(jì),并進(jìn)行了仿真驗(yàn)證,驗(yàn)證了算法的正確性和可行性。

        [1] 陳家偉, 陳杰, 龔春英. 變速風(fēng)力發(fā)電系統(tǒng)瞬態(tài)載荷分析及其優(yōu)化設(shè)計(jì)途徑[J]. 電工技術(shù)學(xué)報(bào), 2015, 30(4): 233-240.

        CHEN Jiawei, CHEN Jie, GONG Chunying. Way to optimally design transient load of variable-speed wind energy conversion system[J]. Transactions of China Electrotechnical Society, 2015, 30(4): 233-240.

        [2] 馬衛(wèi)東. 風(fēng)力發(fā)電機(jī)組最大功率追蹤[J]. 高壓電器, 2012, 48(7): 57-65.

        MA Weidong. Maximum power point tracking for wind turbine[J]. High Voltage Apparatus, 2012, 48(7): 57-65.

        [3] 樊繼永. 基于EKF的電動(dòng)車(chē)用無(wú)刷直流電機(jī)無(wú)傳感器運(yùn)行[J]. 微特電機(jī), 2011, 39(6): 9-11, 14.

        FAN Jiyong. The operation of sensorless BLDCM used in E-bike based on EKF[J]. Small & Special Electrical Machines, 2011, 39(6): 9-11, 14.

        [4] YANG Junyou, ZHAO Jingming, QIN Chengzhi, et al. Model reference adaptive sensorless control for surface permanent magnet synchronous machine under severe starting conditions[C] // Power Electronics and Motion Control Conference, 2004: 1018-1022.

        [5] 張安安, 郭紅鼎, 于兵, 等. 一種風(fēng)電機(jī)組在低風(fēng)速區(qū)間的功率控制方法[J]. 電力系統(tǒng)保護(hù)與控制, 2015, 43(23): 84-88.

        ZHANG Anan, GUO Hongding, YU Bing, et al. A power control method of wind turbines in low wind speed area[J]. Power System Protection and Control, 2015, 43(23): 84-88.

        [6] CHEN Z, TOMITA M, DOKI S, et al. An extended electromotive motive force model for sensorless control of interior permanent-magnet synchronous motors[J]. IEEE Transactions on Industrial Electronics, 2003, 50(2): 288-295.

        [7] 修春波, 任曉, 李艷晴, 等. 基于卡爾曼濾波的風(fēng)速序列短期預(yù)測(cè)方法[J]. 電工技術(shù)學(xué)報(bào), 2014, 29(2): 254-259.

        XIU Chunbo, REN Xiao, LI Yanqing, et al. Short-term prediction method of wind speed series based on Kalman filtering fusion[J]. Transactions of China Electrotechnical Society, 2014, 29(2): 254-259.

        [8] 張宇輝, 陳峰, 李慧敏, 等. 基于小波變換和矩陣束算法的同步電機(jī)參數(shù)辨識(shí)[J]. 電力系統(tǒng)保護(hù)與控制, 2012, 40(9): 87-92.

        ZHANG Yuhui, CHEN Feng, LI Huimin, et al. Parameter identification of synchronous machine based on wavelet transform and matrix pencil algorithm[J]. Power System Protection and Control, 2012, 40(9): 87-92.

        [9] 崔雙喜, 王維慶, 張強(qiáng). 風(fēng)力發(fā)電機(jī)組獨(dú)立變槳魯棒自適應(yīng)槳距角跟蹤控制[J]. 電力系統(tǒng)保護(hù)與控制, 2015, 43(6): 52-57.

        CUI Shuangxi, WANG Weiqing, ZHANG Qiang. Robust adaptive blade pitch angle tracking control for wind turbines[J]. Power System Protecction and Control, 2015, 43(6): 52-57.

        [10]孟令斌, 朱鳳龍. 混合電源及功率預(yù)測(cè)系統(tǒng)在風(fēng)電并網(wǎng)中的應(yīng)用[J]. 電力系統(tǒng)保護(hù)與控制, 2015, 43(13): 79-85.

        MENG Lingbin, ZHU Fenglong. Applications of the mixing power and power forecasting system in wind power[J]. Power System Protecction and Control, 2015, 43(13): 79-85.

        [11]張子泳, 胡志堅(jiān), 李勇匯, 等. 大型雙饋風(fēng)力發(fā)電寫(xiě)小信號(hào)動(dòng)態(tài)建模及附加阻尼控制器設(shè)計(jì)[J]. 電力系統(tǒng)保護(hù)與控制, 2011, 39(18): 127-130.

        ZHANG Ziyong, HU Zhijian, LI Yonghui, et al. Small signal dynamic modelling and additionaller damping controller designing for large wind generation system based on DFIG[J]. Power System Protecction and Control, 2011, 39(18): 127-130.

        [12]何玉林, 吳德俊, 杜靜, 等. 風(fēng)力發(fā)電機(jī)組的自適應(yīng)轉(zhuǎn)矩控制及載荷優(yōu)化[J]. 電網(wǎng)與清潔能源, 2012, 28(3): 79-84.

        HE Yulin, WU Dejun, DU Jing, et al. Adaptive torque control and load optimization of wind turbines[J]. Power System and Clean Energy, 2012, 28(3): 79-84.

        [13] 年珩, 李嘉文, 萬(wàn)中奇. 基于參數(shù)在線辨識(shí)的永磁風(fēng)力發(fā)電機(jī)無(wú)位置傳感器控制技術(shù)[J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2012, 32(12): 146-153.

        NIAN Heng, LI Jiawen, WAN Zhongqi. Sensorless control technique of PMSG for wind power application based on on-line parameter identification[J]. Proceedings of the CSEE, 2012, 32(12): 146-153.

        [14]秦斌, 姜學(xué)想, 朱萬(wàn)力, 等. 基于無(wú)風(fēng)速傳感器的永磁直驅(qū)風(fēng)力發(fā)電系統(tǒng)的直接轉(zhuǎn)矩最大功率跟蹤控制[J]. 電工技術(shù)學(xué)報(bào), 2014, 29(1): 132-137.

        QIN Bin, JIANG Xuexiang, ZHU Wanli, et al. Maximun power point tracking of PMSG wind turbine system based on the direct torque control without wind speed sensors[J]. Transactions of China Electrotechnical Society, 2014, 29(1): 132-137.

        [15] 趙海森, 杜中蘭, 劉曉芳, 等. 基于遞推最小二乘法與模型參考自適應(yīng)法的鼠籠式異步電機(jī)轉(zhuǎn)子電阻在線辨識(shí)方法[J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2014, 34(30): 5386-5394.

        ZHAO Haisen, DU Zhonglan, LIU Xiaofang, et al. An on-line identification method for rotor resistance of squirrel cage induction motors based on recursive least square method and model reference adaptive system[J]. Proceedings of the CSEE, 2014, 34(30): 5386-5394.

        [16]滕予非, 行武, 張宏圖, 等.風(fēng)力發(fā)電系統(tǒng)短路故障特征分析及對(duì)保護(hù)的影響[J]. 電力系統(tǒng)保護(hù)與控制, 2015, 43(19): 29-36.

        TENG Yufei, XING Wu, ZHANG Hongtu, et al.Analysis of characteristics of short circuit fault of wind power system and the impact on the protection[J]. Power System Protection and Control, 2015, 43(19): 29-36.

        (編輯 張愛(ài)琴)

        Study on a rotor speed estimation algorithm of PMSG wind power system

        CHENG Hui1, YANG Keli2, WANG Kejun1, LI Na1

        (1. College of Electrical and Information Engineering, Henan Institute of Engineering, Zhengzhou 451192, China; 2. Training Center for Industry, Zhongyuan University of Technology, Zhengzhou 450007, China)

        A rotor speed estimation algorithm in a direct vector controlled permanent magnet synchronous generator wind energy conversion system is proposed. The method is based on a simple equation obtainer from stator voltage and current. Stator flux information required by the speed estimator is obtained using the stator voltage equation by implementing a programmable low pass filter. The constant gain recursive least squares algorithm is adopted. Then, the rotor position obtained based on the proposed method are used in the rotor-flux oriented vector control of PMSG. The simulation result shows that the control strategy realizes maximum power point tracking control, which verifies the validity and feasibility of the theoretical model and control strategy. This work is supported by National Natural Science Foundation of China (No. 61405054).

        permanent magnet synchronous generator; speed estimation; wind energy conversion systems; sensorless control; maximum power point tracking control

        10.7667/PSPC150739

        國(guó)家自然科學(xué)青年基金項(xiàng)目(61405054)

        2015-08-05;

        2015-12-20

        程 輝(1979-),女,碩士,講師,研究方向?yàn)殡娏﹄娮优c電力傳動(dòng)、風(fēng)力發(fā)電。E-mail: ch.cxy@163.com

        猜你喜歡
        風(fēng)速
        邯鄲市近46年風(fēng)向風(fēng)速特征分析
        基于Kmeans-VMD-LSTM的短期風(fēng)速預(yù)測(cè)
        基于最優(yōu)TS評(píng)分和頻率匹配的江蘇近海風(fēng)速訂正
        基于時(shí)間相關(guān)性的風(fēng)速威布爾分布優(yōu)化方法
        陜西黃土高原地區(qū)日極大風(fēng)速的統(tǒng)計(jì)推算方法
        陜西氣象(2020年2期)2020-06-08 00:54:38
        基于GARCH的短時(shí)風(fēng)速預(yù)測(cè)方法
        快速評(píng)估風(fēng)電場(chǎng)50年一遇最大風(fēng)速的算法
        風(fēng)能(2016年11期)2016-03-04 05:24:00
        考慮風(fēng)切和塔影效應(yīng)的風(fēng)力機(jī)風(fēng)速模型
        GE在中國(guó)發(fā)布2.3-116低風(fēng)速智能風(fēng)機(jī)
        考慮風(fēng)速分布與日非平穩(wěn)性的風(fēng)速數(shù)據(jù)預(yù)處理方法研究
        亚洲综合小综合中文字幕| 久久无码人妻一区二区三区午夜 | 日韩精品视频在线观看无| 永久天堂网av手机版| 国产自偷自偷免费一区| 国内无遮码无码| 麻豆夫妻在线视频观看| 免费在线观看av不卡网站| 精品三级av无码一区| 亚洲欲色欲香天天综合网| 人妻少妇喷水意淫诱惑| 一区二区三区午夜视频在线| 99爱在线精品免费观看| 美女胸又www又黄的网站| 久久久国产精品粉嫩av| 不卡一区二区三区国产| 久久天天躁狠狠躁夜夜av| 91spa国产无码| 久草91这里只有精品| 91九色人妻精品一区二区三区| 无码人妻久久一区二区三区app| 中文字幕第1页中文字幕在| 亚洲精品一区二区三区蜜臀| 最新中文字幕人妻少妇| 99久久人妻精品免费二区| 色爱无码A V 综合区| 国产一区二区三区成人av | 国产精品无码一本二本三本色| 亚洲欧美日韩综合在线观看| 一区二区三区精品偷拍| 77777亚洲午夜久久多喷| 男受被做哭激烈娇喘gv视频| 成人永久福利在线观看不卡| 亚洲黄色精品在线播放| 内射人妻视频国内| 手机看片福利日韩| 蜜桃av一区二区三区| 丝袜人妻一区二区三区| 成人做爰高潮尖叫声免费观看| 无码av永久免费大全| 国产精品国产三级久久|