美章網(wǎng) 資料文庫(kù) 變壓器動(dòng)態(tài)電感參數(shù)辨識(shí)方法范文

變壓器動(dòng)態(tài)電感參數(shù)辨識(shí)方法范文

本站小編為你精心準(zhǔn)備了變壓器動(dòng)態(tài)電感參數(shù)辨識(shí)方法參考范文,愿這些范文能點(diǎn)燃您思維的火花,激發(fā)您的寫(xiě)作靈感。歡迎深入閱讀并收藏。

變壓器動(dòng)態(tài)電感參數(shù)辨識(shí)方法

《中國(guó)電機(jī)工程學(xué)報(bào)》2015年第S1期

摘要:

電力變壓器在交直流混合模式下運(yùn)行時(shí)的電感和電流變化與傳統(tǒng)情況具有明顯區(qū)別,該文結(jié)合電磁耦合機(jī)理提出一種動(dòng)態(tài)電感參數(shù)辨識(shí)方法。建立磁場(chǎng)-電路耦合模型,以時(shí)域電流和動(dòng)態(tài)電感為關(guān)鍵耦合參數(shù),采用棱邊有限元法求解磁場(chǎng)模型,利用局部線性化的方法獲取系統(tǒng)能量,基于能量原理計(jì)算動(dòng)態(tài)電感。采用四階龍格庫(kù)塔法計(jì)算時(shí)域電流,通過(guò)非線性磁場(chǎng)和微分電路迭代實(shí)現(xiàn)磁場(chǎng)與電路的耦合。研究無(wú)直流和交直流混合條件下變壓器的電磁特性,結(jié)合勵(lì)磁飽和狀態(tài),分析動(dòng)態(tài)電感波動(dòng)情況并總結(jié)其變化規(guī)律。通過(guò)仿真計(jì)算與實(shí)驗(yàn)數(shù)據(jù)對(duì)比,驗(yàn)證該文方法的有效性和所得結(jié)果的正確性。

關(guān)鍵詞:

變壓器;動(dòng)態(tài)電感;電磁耦合;能量增量;交直流混合

高壓直流輸電單極大地回路運(yùn)行[1-2]和地球磁暴條件下交流電網(wǎng)的接地變壓器-輸電線路-大地回路中存在(準(zhǔn))直流電流[3-4],同時(shí)大量非線性電力電子元件(如相控交流負(fù)載、整流器等)的投入和運(yùn)行均可能產(chǎn)生直流分量[5],從而在交流電網(wǎng)中形成了交直流混雜的特殊環(huán)境,造成變壓器等電磁設(shè)備直流偏磁,由此導(dǎo)致變壓器勵(lì)磁飽和、局部過(guò)熱、振動(dòng)噪聲等一系列問(wèn)題[6-8]。另一方面,變壓器直流偏磁時(shí)參數(shù)發(fā)生畸變,使得傳統(tǒng)差動(dòng)、諧波制動(dòng)等保護(hù)措施不再適用,以致出現(xiàn)大范圍停電事故[9-10],嚴(yán)重威脅了設(shè)備乃至電網(wǎng)的穩(wěn)定運(yùn)行。目前,變壓器乃至電網(wǎng)遭受直流擾動(dòng)危害的嚴(yán)重性已引起人們廣泛關(guān)注,但是關(guān)于變壓器直流偏磁的異常故障和相關(guān)保護(hù)未做深入研究。因此,結(jié)合變壓器直流偏磁的特點(diǎn)研究關(guān)鍵電磁參數(shù)受直流擾動(dòng)的影響對(duì)設(shè)備及系統(tǒng)的安全運(yùn)行和保護(hù)設(shè)計(jì)都具有十分重要的意義。

國(guó)內(nèi)外學(xué)者針對(duì)變壓器直流偏磁問(wèn)題主要開(kāi)展了電磁特性的模擬分析與實(shí)驗(yàn)研究。通過(guò)能量守恒原理對(duì)變壓器外端電路與內(nèi)部磁路間的能量轉(zhuǎn)換過(guò)程進(jìn)行描述和分析能有效克服傳統(tǒng)方法的建模復(fù)雜性,從而提升計(jì)算效率[11-13]。文獻(xiàn)[14]依據(jù)能量守恒原理,提出一種變壓器電感參數(shù)測(cè)量的新方法,能有效對(duì)變壓器內(nèi)部磁路進(jìn)行模擬、具有一定準(zhǔn)確性。文獻(xiàn)[15]利用能量轉(zhuǎn)換手段分析利用變壓器電感參數(shù)設(shè)計(jì)保護(hù)判據(jù),具體實(shí)驗(yàn)驗(yàn)證了思路的可行性,但忽略了勵(lì)磁遭受直流擾動(dòng)因數(shù)的影響。上述研究均認(rèn)為變壓器電感近似為常數(shù),而在交直流混合條件下,變壓器勵(lì)磁飽和加劇,電流畸變,漏磁增大[16-20],因此,對(duì)應(yīng)的電感參數(shù)也可能出現(xiàn)較大變化,但少有文獻(xiàn)進(jìn)行深入的研究。

本文基于變壓器電磁耦合機(jī)理研究其直流偏磁問(wèn)題,結(jié)合能量擾動(dòng)原理建立計(jì)及直流偏磁的間接場(chǎng)路耦合模型,提出一種直流偏磁電感參數(shù)辨識(shí)方法,模擬分析變壓器直流偏磁電磁特性,利用時(shí)域電流和動(dòng)態(tài)電感等關(guān)鍵參數(shù)表征變壓器直流偏磁飽和狀態(tài),研究不同條件下變壓器動(dòng)態(tài)電感與勵(lì)磁、諧波及內(nèi)部漏磁在直流擾動(dòng)下的波動(dòng)情況及相互關(guān)系,總結(jié)其變化規(guī)律,并通過(guò)仿真與實(shí)驗(yàn)對(duì)比,對(duì)本文方法進(jìn)行驗(yàn)證。

1基于電磁耦合的動(dòng)態(tài)電感辨識(shí)原理

1.1電路模型變壓器基本電磁耦合原理[21]如圖1所示。圖1中,u1為交流激勵(lì)源,Φ1和Φ2分別為鐵芯的主磁通和漏磁通,L和M分別為變壓器的自感和互感。

1.2磁場(chǎng)模型采用基于矢量磁位A的棱邊有限元法[22-27],不考慮磁滯效應(yīng),建立磁場(chǎng)模型。

1.3能量擾動(dòng)原理依據(jù)能量轉(zhuǎn)換原理,由時(shí)變系統(tǒng)參數(shù)求取動(dòng)態(tài)電感參數(shù)。若線圈勵(lì)磁電流改變量為Δip(0≤Δ≤1)時(shí),將電源總能量與動(dòng)態(tài)電感和電流關(guān)聯(lián)。利用繞組電流i與動(dòng)態(tài)電感LD,即可實(shí)現(xiàn)電磁耦合,其計(jì)算原理如圖3所示。

2仿真計(jì)算

針對(duì)變壓器直流偏磁問(wèn)題,采用時(shí)域場(chǎng)路耦合模型求解動(dòng)態(tài)電感參數(shù),其中電路模型通過(guò)編寫(xiě)四階龍格庫(kù)塔法程序進(jìn)行求解,磁場(chǎng)模型通過(guò)ANSYS軟件實(shí)現(xiàn)。變壓器的額定參數(shù)見(jiàn)表1,鐵心硅鋼片型號(hào)為DW360-50。計(jì)算變壓器空載運(yùn)行時(shí)的耦合參數(shù),不同電壓情況下的結(jié)果如圖4所示。圖4中i1表示原邊電流,與勵(lì)磁特性有關(guān),當(dāng)i1的數(shù)值較大時(shí),表明變壓器的勵(lì)磁狀態(tài)趨于飽和;當(dāng)i1數(shù)值較小時(shí),則說(shuō)明變壓器勵(lì)磁狀態(tài)不飽和。另一方面,計(jì)算過(guò)程受相關(guān)參數(shù)影響存在振蕩現(xiàn)象,如圖4(a)中的區(qū)域I和II所示。分析不同電壓下的變壓器勵(lì)磁特性,結(jié)果表明,隨著電壓升高,變壓器的勵(lì)磁狀態(tài)逐漸由不飽和狀態(tài)轉(zhuǎn)向飽和狀態(tài)。圖4(b)中原邊電感L1的波峰對(duì)應(yīng)變壓器不飽和勵(lì)磁狀態(tài),原邊電感L1的波谷對(duì)應(yīng)變壓器的飽和勵(lì)磁狀態(tài)。在直流擾動(dòng)狀態(tài)下分析耦合參數(shù)的變化特性。

變壓器原邊空載電流最大值為Im,直流源UDC接入時(shí)產(chǎn)生的直流電流大小為IDC。存在直流時(shí),變壓器鐵心勵(lì)磁受直流影響出現(xiàn)偏置,i1波形發(fā)生畸變,對(duì)應(yīng)的L1也發(fā)生明顯變化。在直流偏置的正半周,勵(lì)磁飽和程度加劇,鐵磁材料磁導(dǎo)率μ下降,L1減小;在直流偏置的負(fù)半周,勵(lì)磁飽和程度減弱,μ上升,L1增大。并隨著直流水平升高,ie波形畸變嚴(yán)重,L1正負(fù)半周的不對(duì)稱程度加劇。通過(guò)仿真分析,可以得出以下結(jié)論。1)變壓器的耦合參數(shù)在交變電磁耦合的過(guò)程中是時(shí)變的,仿真表明動(dòng)態(tài)電感參數(shù)并不是恒定的常數(shù),其波動(dòng)情況取決于變壓器的運(yùn)行狀態(tài)和勵(lì)磁飽和特性。2)在交直流混合條件下,變壓器勵(lì)磁飽和程度加深,電流和電感在每周期的波形不再對(duì)稱,并隨著直流水平升高,這種不對(duì)稱程度加劇。

3實(shí)驗(yàn)驗(yàn)證

實(shí)驗(yàn)變壓器型號(hào)為BK300,如圖6所示,參數(shù)見(jiàn)表1,實(shí)驗(yàn)接線原理如圖7所示。在不同交直流混合條件下進(jìn)行實(shí)驗(yàn)量測(cè),電流結(jié)果見(jiàn)圖8。由圖5和圖8可知,空載運(yùn)行直流偏磁時(shí),電流的計(jì)算結(jié)果與實(shí)驗(yàn)測(cè)量基本相同;兩者存在的誤差可能由磁滯所導(dǎo)致,在直流偏磁情況下兩者誤差更小。設(shè)計(jì)低通濾波模塊降低實(shí)驗(yàn)過(guò)程中存在的高次諧波分量影響,以一次側(cè)交流電壓有效值UAC=50V時(shí)的實(shí)驗(yàn)為例,結(jié)果如圖9所示。利用時(shí)域差分模塊測(cè)量記錄變壓器動(dòng)態(tài)電感,結(jié)果如圖10所示。無(wú)直流時(shí)不同電壓下的動(dòng)態(tài)電感結(jié)果如圖11所示。當(dāng)交流電壓有效值UAC=50V時(shí),變壓器勵(lì)磁處于不飽和區(qū),i1的波形近似呈正弦波,L1波動(dòng)較小。隨著UAC升高,勵(lì)磁飽和程度加深,i1的波形趨于尖頂波,L1波動(dòng)逐漸增大,具體數(shù)據(jù)及波形見(jiàn)附錄A。

交直流混合實(shí)驗(yàn)原理同上,UAC=220V時(shí)的電感結(jié)果見(jiàn)圖12。當(dāng)UAC較低時(shí),變壓器勵(lì)磁不飽和,電流波動(dòng)較小,結(jié)果存在一定誤差。隨著UAC升高,變壓器勵(lì)磁趨于飽和,電流波動(dòng)增大,誤差逐漸減小,電感參數(shù)在直流擾動(dòng)的影響下出現(xiàn)半周不對(duì)稱的問(wèn)題。結(jié)果表明,動(dòng)態(tài)電感受直流擾動(dòng)時(shí)的變化情況與仿真分析所得結(jié)論基本一致,驗(yàn)證了本文方法的正確性。表2中I1為空載電流有效值。通過(guò)對(duì)比變壓器在不同條件下耦合參數(shù)的數(shù)值,表明電流和電感受直流影響發(fā)生顯著變化。交直流混合時(shí),勵(lì)磁飽和程度加深,I1增加,L1AV減小。因此,可以根據(jù)不同直流水平下的電感參數(shù)變化情況制定合理有效地保護(hù)策略。

4結(jié)論

利用變壓器電磁耦合模型計(jì)算交直流混合條件下的動(dòng)態(tài)電感參數(shù),通過(guò)研究其變化特性,得出以下結(jié)論。1)變壓器電磁耦合模型能夠有效模擬交直流混合運(yùn)行時(shí)的電磁特性,該過(guò)程中的動(dòng)態(tài)電感具有時(shí)變性,其波動(dòng)情況與變壓器的勵(lì)磁飽和狀態(tài)對(duì)應(yīng)。無(wú)直流時(shí)隨著交流電壓升高,變壓器勵(lì)磁飽和加劇,電感均值先增大后減小;受直流擾動(dòng)影響,電流與電感發(fā)生畸變。2)仿真計(jì)算與實(shí)驗(yàn)結(jié)果對(duì)比有效驗(yàn)證了本文方法的正確性,為進(jìn)一步開(kāi)展變壓器等電磁設(shè)備在交直流混合條件下的電磁特性分析和相關(guān)保護(hù)研究奠定基礎(chǔ)。

參考文獻(xiàn)

[1]蒯狄正,萬(wàn)達(dá),鄒云.直流輸電地中電流對(duì)電網(wǎng)設(shè)備影響的分析與處理[J].電力系統(tǒng)自動(dòng)化,2005,29(2):81-82.KuaiDizhang,WanDa,ZouYun.AnalysisandhandlingoftheimpactofgeomagneticallyinducedcurrentuponelectricnetworkequipmentinDCtransmission[J].AutomationofElectricPowerSystem,2005,29(2):81-82(inChinese).

[2]李曉萍,文習(xí)山.三相五柱變壓器直流偏磁計(jì)算研究[J].中國(guó)電機(jī)工程學(xué)報(bào),2010,30(1):127-131.LiXiaoping,WenXishan.DCbiascomputationstudyonthree-phasefivelimbstransformer[J].ProceedingsoftheCESS,2010,30(1):127-131(inChinese).

[3]趙小軍,李琳,程志光,等.應(yīng)用諧波平衡有限元法的變壓器直流偏磁現(xiàn)象分析[J].中國(guó)電機(jī)工程學(xué)報(bào),2010,30(21):103-108.ZhaoXiaojun,LiLin,ChengZhiguang,etal.AnalysisoftheDCbiasphenomenonintransformersbasedonharmonic-balancedfiniteelementmethod[J].ProceedingsoftheCSEE,2010,30(21):103-108(inChinese).

[4]LiuCM,LiuLG,PirjolaR.Geomagneticallyinducedcurrentsinthehigh-voltagepowergridinChina[J].IEEETransactionsonPowerDelivery,2010,24(4):2368-2374.

[5]劉連光,劉春明,張冰.磁暴對(duì)我國(guó)特高壓電網(wǎng)的影響[J].地球物理學(xué)報(bào),2008,51(4):976-981.LiuLianguang,LiuChunming,Zhangbing.EffectsofgeomagneticstormonUHVpowergridsinChina[J].ChineseJournalofGeophysics,2008,51(4):976-981(inChinese).

[6]陳國(guó)棟,宋晉峰,張亮,等.動(dòng)態(tài)電壓恢復(fù)器注入變壓器的直流偏磁抑制策略[J].中國(guó)電機(jī)工程學(xué)報(bào),2014,34(28):4983-4989.ChenGuodong,SongJinfeng,ZhangLiang,etal.ControlstrategyforeliminatingDCmagneticfluxintheinjectiontransformerofadynamicvoltagerestorer[J].ProceedingsoftheCSEE,2014,34(28):4983-4989(inChinese).

[7]鄧祥力,熊小伏,高亮,等.基于參數(shù)辨識(shí)的變壓器繞組變形在線監(jiān)測(cè)方法[J].中國(guó)電機(jī)工程學(xué)報(bào),2014,34(28):4950-4958.DengXiangli,XiongXiaofu,GaoLiang,etal.Methodofon-linemonitoringoftransformerwindingdeformationbasedonparameteridentification[J].ProceedingsoftheCSEE,2014,34(28):4950-4958(inChinese).

[8]劉曲,李立涅,鄭健超.考慮海洋影響的直流輸電單極大地運(yùn)行時(shí)變壓器中性點(diǎn)直流電流研究[J].電網(wǎng)技術(shù),2007,31(2):57-60.LiuQu,LILicheng,ZhengJianchao.StudyonDCcurrentthroughtransformerneutralscausedbygroundreturnoperationmodeofHVDCsystemwithseainfluenceconsidered[J].Powersystemtechnology,2007,31(2):57-60(inChinese).

[9]焦在濱,馬濤,屈亞軍,等.基于勵(lì)磁電感參數(shù)識(shí)別的快速變壓器保護(hù)[J].中國(guó)電機(jī)工程學(xué)報(bào),2014,54(10):158-163.JiaoZaibin,MaTao,QuYajun,etal.Anovelexcitationinductance-basedtransformerprotectionscheme[J].ProceedingsoftheCSEE,2014,54(10):1658-1663(inChinese).

[10]鄭濤,陳佩璐,劉連光,等.計(jì)及直流偏磁的電流互感器特性對(duì)差動(dòng)保護(hù)的影響[J].電力系統(tǒng)自動(dòng)化,2012,36(20):89-93.ZhengTao,ChenPeilu,LiuLianguang,etal.TransferringcharacteristicsofcurrenttransformeraffectedbyDCmagneticviasanditsimpactondifferentialprotection[J].AutomationofElectricPowerSystems,2012,36(20):89-93(inChinese).

[11]ZhaoXJ,LiL,ChengZG.Characteristicanalysisofthesquarelaminatedcoreunderdc-biasedmagnetizationbythefixed-pointharmonic-balancedFEM[J].IEEETransactionsonMagnetics,2012,48(2):747-750.

[12]史海濤,王建元,潘超,等.變壓器電磁耦合模型中動(dòng)態(tài)電感的計(jì)算方法[J].電工電能新技術(shù),2015,34(8):69-74.ShiHaitao,WangJianyuan,CaiGuowei,etal.Dynamicinductancecomputationinelectromagneticcoupledmodeloftransformer[J].AdvancedTechnologyofElectricalEngineeringandEnergy,2015,34(8):69-74(inChinese).

[13]潘超,何雯婷,郝文波,等.一種變壓器時(shí)域場(chǎng)路耦合計(jì)算的自適應(yīng)算法[J].電工技術(shù)學(xué)報(bào),2014,29(12):113-119.PanChao,HeWenting,HaoWenbo,etal.Aself-adaptivealgorithmoftime-domainmagneticfieldandelectriccircuitcoupledcomputationfortransformer[J].TransactionsofChinaElectrotechnicalSociety,2014,29(12):113-119(inChinese).

[14]HoSL,NiuSX.Apower-balancedtime-steppingfiniteelementmethodfortransientmagneticfieldcomputation[J].IEEETransactionsonMagnetics,2012,48(2):291-294.

[15]王雪,王增平.基于有限元法的變壓器電感參數(shù)計(jì)算方法研究[J].電力系統(tǒng)保護(hù)與控制,2009,37(24):11-14.WangXue,WangZengping.Aninvestigationontransformerinductancecalculationbasedonfiniteelementmethod[J].PowerSystemProtectionandControl,2009,37(24):11-14(inChinese).

[16]ZhangHaijun,YangBin.Dynamicdeformationanalysisofpowertransformerwindingsinshort-circuitfaultbyFEM[J].IEEETransactionsonAppliedSuperconductivity,2014,28(3):502-504.[17]潘超,王澤忠,楊敬瑀,等.變壓器直流偏磁瞬態(tài)場(chǎng)路耦合計(jì)算的穩(wěn)定性分析[J].電工技術(shù)學(xué)報(bào),2012,27(12):226-239.PanChao,WangZezhong,YangJingyu,etal.Stabilityanalysisbasedontransientmagnetic-circuitcoupledmethodforDC-biasedtransformer[J].TransactionsofChinaElectrotechnicalSociety,2012,27(12):226-239(inChinese).

[18]索南加樂(lè),焦在濱,康小寧,等.Y/Δ接線變壓器漏感參數(shù)的識(shí)別方法[J].中國(guó)電機(jī)工程學(xué)報(bào),2008,28(13):84-90.SuonanJ,JiaoZaibin,KangXiaoning,etal.AlgorithmtoidentifyleakageinductancesofpowertransformerwithY-deltaconnection[J].ProceedingsoftheCSEE2008,28(13):84-90(inChinese).

[19]王增平,馬靜.基于等效瞬時(shí)漏感與回路方程的變壓器保護(hù)原理[J].中國(guó)電機(jī)工程學(xué)報(bào),2007,27(19):39-44.WangZengping,MaJing.Anovelprincipleoftransformerprotectionbasedonequivalentinstantaneousleakageinductanceandloopequation[J].ProceedingsoftheCSEE2007,27(19):39-44(inChinese).

[20]林湘寧,劉沛,高艷.基于數(shù)學(xué)形態(tài)的電流互感器飽和識(shí)別判據(jù)[J].中國(guó)電機(jī)工程學(xué)報(bào),2005,25(5):44-48.LinXiangning,LiuPei,GaoYan.Anovelmethodtoidentifythesaturationofthecurrenttransformerusingmathematicalmorphology[J].ProceedingsoftheCSEE2005,25(5):44-48(inChinese).

[21]潘超,王澤忠,李海龍,等.基于瞬態(tài)場(chǎng)路耦合模型的變壓器直流偏磁計(jì)算[J].電工技術(shù)學(xué)報(bào),2013,28(5):174-181.PanChao,WangZezhong,LiHailong,etal.DC-biascalculationforsinglephasetransformersbasedontransientfield-circuitcoupledmodel[J].TransactionsofChinaElectrotechnicalSociety2013,28(5):174-181(inChinese).

[22]梁振光,唐任遠(yuǎn).采用場(chǎng)路耦合的三維有限元法分析變壓器突發(fā)短路過(guò)程[J].中國(guó)電機(jī)工程學(xué)報(bào),2003,23(3):137-140.LiangZhenguang,TangRenyuan.Faultsimulationpowertransformersusing3Dfiniteelementmodelcoupledtoelectriccircuitequations[J].ProceedingsoftheCSEE,2003,23(3):137-140(inChinese).

[23]謝若冰,陳橋夫,康崇皓,等.基于組合式場(chǎng)路耦合法的多繞組變壓器建模與阻抗參數(shù)設(shè)計(jì)[J].中國(guó)電機(jī)工程學(xué)報(bào),2009,29(9):104-111.XieRuobing,ChenQiaofu,KangChonghao,etal.Modelingandimpedanceparameterdesignformulti-windingtransformerbasedoncombinedfield-circuitcoupledmethod[J].ProceedingsoftheCSEE,2009,29(9):104-111(inChinese).

[24]ChishanY.DetectionandcorrectionofsaturatedcurrenttransformermeasurementusingdecayingDCcomponents[J].IEEETransactionsonPowerDelivery,2010,25(3):1340-1347.

[25]楊浩,曾鑫,羅建.基于約束最小二乘法的變壓器三相漏電感辨識(shí)算法[J].電力系統(tǒng)自動(dòng)化,2009,30(13):68-72.YangHao,ZengXin,LuoJian.Athree-phasetransformerleakageinductancesparameteridentificationalgorithmbasedontheconstrainedleastsquaresmethod[J].AutomationofElectricPowerSystem,2009,30(13):68-72(inChinese).

[26]索南加樂(lè),張健康,張軍民,等.交直流混聯(lián)系統(tǒng)對(duì)變壓器保護(hù)性能的影響及解決措施[J].電力系統(tǒng)自動(dòng)化,2010,34(3):101-106.SuonanJiale,ZhangJiankang,ZhangJunmin,etal.TheadverseimpactsontransformerprotectionperformancebyhybridAC-DCtransmissiongridsystemandcounter-measures[J].AutomationofElectricPowerSystems,2010,34(3):101-106(inChinese).

[27]MauricioVFerreiraDL,JeanVL,etal.Three-phasetransformermodelingusingavectorhysteresismodelandincludingtheeddycurrentandtheanomalouslosses[J].IEEETransactionsonMagnetics,2010,46(8):3201-3204.

作者:潘超 孟濤 王夢(mèng)純 蔡國(guó)偉 單位:東北電力大學(xué)電氣工程學(xué)院

主站蜘蛛池模板: 久久精品视频16| 午夜视频在线观看一区二区| 99精品国产在热久久| 无码日韩精品一区二区免费| 亚洲午夜久久久精品影院| 男人肌肌桶女肌肌网站| 成人自拍小视频| 亚洲av无码不卡在线播放| 老师小sao货水好多真紧h视频| 国产精品2019| 99久久精品费精品国产| 成人综合激情另类小说| 久久精品国产亚洲av日韩| 欧美换爱交换乱理伦片免费| 国产卡一卡二卡3卡乱码免费| 2021国内精品久久久久精免费| 女人18毛片水真多免费播放| 久久久久777777人人人视频| 欧美jlzz18性欧美| 亚洲欧美日韩一区在线观看| 青娱乐精品视频| 国产精品久久久小说| 两性色午夜视频免费播放| 日韩欧美中文字幕在线视频| 亚洲最大无码中文字幕| 男人天堂官方网站| 午夜爽爽爽男女免费观看hd| 超清中文乱码字幕在线观看| 国产特级毛片aaaaaaa高清 | 男女一进一出无遮挡黄| 啊灬啊灬啊灬快灬高潮少妇| 青青青国产精品一区二区| 国产欧美日韩另类一区乌克兰| 2021国产麻豆剧传媒官网| 国语高清精品一区二区三区| xxxxx做受大片在线观看免费 | 天天爽夜夜爽夜夜爽| 一级片黄色免费| 手机国产乱子伦精品视频| 久久亚洲精品成人777大小说| 暖暖在线日本免费中文|