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

計及鐵心損耗和磁滯效應改進的變壓器模型

2016-12-27 05:58:30李慧奇李曉孟李金忠張書琦閆彥剛
電工技術學報 2016年21期
關鍵詞:變壓器模型

李慧奇 李曉孟 李金忠 張書琦 閆彥剛

(1.新能源電力系統國家重點實驗室(華北電力大學) 保定 071003 2.中國電力科學研究院 北京 100192)

?

計及鐵心損耗和磁滯效應改進的變壓器模型

李慧奇1李曉孟1李金忠2張書琦2閆彥剛1

(1.新能源電力系統國家重點實驗室(華北電力大學) 保定 071003 2.中國電力科學研究院 北京 100192)

分析電力系統暫態過程或鐵磁諧振等問題時,建立變壓器精確的等效非線性模型非常關鍵。基于變壓器的空載實驗,利用非線性電阻和考慮磁滯效應的非線性電感的并聯來等效變壓器鐵心勵磁支路,并給出了變壓器的等效電路模型和勵磁支路參數特性的計算方法。在工頻條件下,分別建立了模型變壓器和單相三柱變壓器JMB-500的等效模型,計算結果與測量結果吻合很好。最后,建立了模型變壓器在不同頻率下的等效模型,該文的計算結果與測量結果也同樣吻合,證明了該文方法的有效性和準確性。

損耗 變壓器 鐵心 磁滯 模型

0 引言

變壓器是電力系統中重要的電氣設備之一,針對變壓器模型,國內外學者進行了大量的研究工作,采用不同的方法對變壓器進行了建模。文獻[1]中,考慮到變壓器鐵心的非線性主要受飽和、磁滯和渦流的影響,指出磁鏈-電流曲線是鐵心飽和非線性模型的重要一部分。文獻[2]用非線性電阻與非線性電感并聯來模擬變壓器,并用數學函數來表示鐵心勵磁特性。文獻[3]用基于損耗函數磁滯模型來描述非線性電感的磁滯回線,但缺少物理機理,并且解析式計算中需要已知鐵心的幾何以及物理參數等,這些參數不易得到。文獻[4,5]基于變壓器飽和非線性特性進行了建模,但沒有考慮到磁滯的影響。文獻[6-8]提出了變壓器時域模型和頻域模型,計算過程中需要通過實驗測得變壓器的幾何以及物理參數,在很多情況下,這些參數測量不易實現。變壓器勵磁支路參數的準確性對于變壓器研究至關重要,特別是在分析變壓器鐵磁諧振、電磁暫態過程、直流偏磁及勵磁涌流等問題[9-11]時,需建立變壓器的準確模型。文獻[12]對330 kV系統中切除空載變壓器的過電壓進行了建模研究,利用所建變壓器非線性模型所得的過電壓幅值為396.519 kV,較變壓器線性模型過電壓幅值582.036 kV有了很大改進,但所建變壓器模型并未考慮到磁滯的影響。文獻[13]對變壓器電磁暫態仿真進行了建模研究,所建變壓器模型同樣未考慮到磁滯的影響,所計算勵磁電流與PSCAD仿真結果仍有一定差距。上述所提出的變壓器模型,以韋-安曲線描述鐵心非線性磁化特性[14],其勵磁支路參數模型的確定并未體現出磁滯對變壓器的影響,這與變壓器的實際工作機理不完全一致。

本文基于變壓器空載實驗,根據變壓器鐵心磁化機理和電路原理,對變壓器勵磁支路等效電路模型中非線性電阻和非線性電感參數進行改進,提出了一種考慮鐵磁材料磁滯特性和鐵心非線性的變壓器等效電路模型。所建的變壓器模型,其勵磁支路的非線性電阻和非線性電感能夠準確反映出變壓器勵磁過程。依據本文方法,在兩種不同鐵磁材料的變壓器中進行了實驗測量和仿真計算,模型的計算結果和實驗結果吻合很好。最后分析了不同頻率下變壓器對應的等效模型,本文模型的計算結果與測量結果吻合很好,證明了本文方法的有效性和準確性。

1 改進變壓器新模型

1.1 變壓器空載的損耗特性

圖1為單相變壓器的等效電路,圖中用非線性電阻RN及非線性電感LN替代變壓器勵磁阻抗。由于變壓器漏阻抗相對于勵磁阻抗很小,本文忽略[15]。

圖1 單相變壓器的等效電路Fig.1 The equivalent circuit of single-phase transformer

對變壓器進行空載實驗,測量變壓器有效值U、勵磁電流有效值I和空載損耗P,記錄變壓器一次電流i波形和一次電壓u波形。變壓器在正弦激勵情況下,有

(1)

式中,u為電壓瞬時值;U為電壓有效值;ω為電壓角頻率。

非線性電阻RN參數特性可以用一條過原點的曲線u-iR描述。將u-iR曲線上原點到(uk,iRk)共k+1個點分段線性化為k段,有

(2)

式中,j=1,2,3,…,k。

由電阻功率定義,結合式(2)有

(3)

式中,ωtRj= arccos(uk-j/uk),j=1,2,3,…,k-1。

由空載變壓器鐵心損耗可知,在正弦交變磁場下鐵心的損耗由磁滯損耗、渦流損耗和附加損耗[16-19]三部分組成。式(3)在計算變壓器勵磁支路參數時,將變壓器的空載損耗完全歸結于一個非線性電阻的損耗,這與空載時變壓器鐵心損耗機理不一致。本文根據變壓器鐵心磁化機理和電路原理,進一步改進了變壓器模型。

1.2 改進變壓器模型參數的計算

基于1.1節中求得的u-iR曲線和變壓器的一次電流i,由基爾霍夫電流定律可求出電感支路瞬時電流

iL=i-iR

(4)

非線性電感的磁鏈為

(5)

由式(4)和式(5)所求得的非線性電感參數特性Ψ-iL曲線是一條回線,而不是一條簡單的單值曲線,該回線考慮了非線性電感的磁滯效應。則此時非線性電感LN對應損耗PL為

(6)

根據鐵磁材料的損耗機理,變壓器空載損耗應為變壓器等效模型中非線性電阻的損耗與非線性電感的損耗(這部分損耗由磁滯效應引起)之和,故對非線性電阻所消耗的損耗做如下修正

(7)

重復式(2)~(7),直至求得的非線性電阻的損耗滿足誤差,即△P=PRn+1-PRn<ε則停止循環,得到非線性電阻和非線性電感的最終參數。

2 實驗驗證和應用

為了驗證本文提出的變壓器模型的有效性,對一臺電壓比為1∶1的變壓器模型進行了建模和驗證。變壓器的一次側、二次側的匝數均為700匝,鐵心為電力變壓器常用硅鋼30Q120,鐵心的截面積為30 mm3.6 mm。電源為WF1973信號源與功率放大器NF4510構成的標準正弦電壓源,變壓器的空載損耗、電壓和勵磁電流同時由精密功率分析儀WT3000測得。通過空載實驗獲得的變壓器的空載損耗、一次電壓和勵磁電流的有效值見表1。

表1 變壓器空載實驗數據Tab.1 Transformer no-load test data

由實驗數據和本文提出的改進建模方法,得到的變壓器等效模型的非線性電阻RN的參數u-iR曲線如圖2所示。非線性電感LN的參數Ψ-iL曲線如圖3所示,圖中的每一條Ψ-iL回線對應圖2中u-iR曲線上的一個點。為便于觀察,圖3僅給出圖2中7個點對應的回線,對應的電壓值分別為3.533 V、14.191 V、19.438 V、22.992 V、26.527 V、29.368 V和30.067 V。圖4為勵磁電壓值為30.067 V時對應的一條非線性電感的韋-安回線。從圖4中可明顯地看出,改進模型的非線性電感不再是單調曲線而是考慮了磁滯影響的回線,該回線反映了變壓器在空載時的磁滯效應。

圖2 變壓器模型的非線性電阻參數Fig.2 The parameter of non-linear resistor of transformer’s model

圖3 變壓器模型非線性電感的參數Fig.3 The parameter of non-linear inductance of transformer’s model

圖4 變壓器模型的磁鏈電流特性曲線Fig.4 Flux-current characteristic curve of transformer’s model

為驗證所建變壓器模型的正確性,圖5給出了電壓為30.067 V時實驗測得的勵磁電流、本文算法及未考慮磁滯效應所得到的勵磁電流的對比。可以看出,未考慮磁滯效應所得到的勵磁電流為每個峰值左右對稱的曲線,這與實驗結果不吻合。

圖5 變壓器勵磁電流計算和測量結果的對比Fig.5 Comparison between calculated and measured results of excited currents of transformer

然而要建立一個完整的,能夠反映空載、負載、過載等多種工況下的變壓器模型,可基于一定數量的實驗數據結合神經網絡算法(從本文的工作經驗,數據以15組為宜)得到任意電壓下變壓器等效模型的參數。例如計算表1之外的任一電壓(如22.27 V)的變壓器參數。圖6為通過神經網絡算法結合本文模型得到的勵磁電流和實際勵磁電流對比結果。從圖5、圖6中可看出:基于本文方法得到的變壓器模型的勵磁支路參數特性更符合變壓器空載損耗機理。從計算結果與測量結果的對比可明顯看出,本文方法與實驗結果吻合很好,在計算精度上比以往發表的結果有了很大的提高。

圖6 變壓器勵磁電流計算和測量結果的對比Fig.6 Comparison between calculated and measured results of excited currents of transformer

為進一步驗證方法的有效性,對一臺單相三柱式變壓器JMB-500進行了實驗和建模。變壓器JMB-500等效模型的非線性電阻RN和非線性電感LN如圖7和圖8所示。圖9給出了電壓在23 V時本文方法的結果、測量結果及與未考慮磁滯效應的勵磁電流的對比。結果表明,本文所建的模型考慮了鐵心磁滯效應的影響,得到的變壓器勵磁支路參數能更準確地反映變壓器的電磁過程。模型計算結果和實驗結果吻合很好,驗證了方法的準確性。

圖7 JMB-500變壓器的非線性電阻參數Fig.7 The parameter of non-linear resistor of JMB-500

圖8 JMB-500非線性電感的參數Fig.8 The parameter of non-linear inductance of JMB-500

圖9 JMB-500勵磁電流計算和測量結果的對比Fig.9 Comparison between calculated and measured results of excited currents of JMB-500

3 不同頻率下變壓器模型

近年來,隨著電力電子技術的廣泛應用,電力系統諧波污染[20-23]日益嚴重。不同頻率正弦激勵下變壓器的模型是研究諧波對變壓器運行特性影響的基礎。在不同頻率下變壓器的建模有多種方法[24,25],應用本文改進的變壓器建模方法,在不同頻率下對模型變壓器進行了空載實驗、建模和驗證。

圖10給出了頻率分別為50 Hz、100 Hz、150 Hz以及200 Hz時的激磁支路非線性電阻RN參數特性(由于實驗設備功率有限,功率放大器NF4510的最大輸出功率為1 kV·A,實驗的最高頻率為200 Hz)。圖11僅給出了頻率為150 Hz對應的電感的Ψ-iL回線。圖12給出了在150 Hz條件下電壓為90.201 V時勵磁電流的測量結果、本文算法及未考慮磁滯算法所得勵磁電流的對比結果。圖13給出了頻率分別為50 Hz、100 Hz、150 Hz以及200 Hz,相同磁感應強度B時所對應的非線性電感參數。可看出,不同頻率下,變壓器損耗特性和勵磁特性均不相同,在不同頻率下不能僅用工頻下的變壓器參數特性來代替非工頻下的情況,需要建立不同頻率下的變壓器模型。應用本文方法所建變壓器模型能夠準確地計算不同頻率下變壓器勵磁情況,是研究含有高次諧波工況下變壓器模型的基礎。限于篇幅,關于本文中所建變壓器模型求解及復雜諧波下變壓器建模等問題不再展開。

圖10 不同頻率下非線性電阻參數Fig.10 Nonlinear resistance parameters at different frequencies

圖11 變壓器模型的非線性電感參數Fig.11 The parameter of non-linear inductance of transformer’s model

圖12 變壓器勵磁電流計算和測量結果的對比Fig.12 Comparison between calculated and measured results of excited currents of transformer

圖13 不同頻率下非線性電感參數Fig.13 Nonlinear inductance parameters at different frequencies

4 結論

變壓器是電力系統中重要的電氣設備之一,分析變壓器鐵磁諧振、電磁暫態過程、直流偏磁以及勵磁涌流等問題需要建立準確的變壓器模型。本文利用非線性電阻與考慮磁滯效應的非線性電感的并聯來建立變壓器模型,并給出了變壓器勵磁支路非線性電阻和非線性電感參數的計算方法。在工頻條件下,針對一臺模型變壓器和一臺實際的單相三柱變壓器,建立了變壓器的等效模型,模型計算結果與測量結果吻合很好,證明了本文方法的有效性和準確性。最后,針對模型變壓器,建立了不同頻率下的變壓器模型,所建模型的計算結果與測量結果吻合很好。本文工作對于研究變壓器鐵磁諧振、電磁暫態過程、直流偏磁和勵磁涌流等問題具有重要意義,也是研究含高次諧波工況下變壓器模型的基礎。

[1] Swift G W.Power transformer core behavior under transient conditions[J].IEEE Transactions on Power Apparatus and Systems,1971,PAS-90(5):2206-2210.

[2] Santesmases J G,Ayala J,Cachero A H.Analytical approximation of dynamic hysteresis loop and its application to a series ferroresonant circuit[J].Proceedings of the Institution of Electrical Engineers,1970,117(1):234-240.

[3] Chua L O,Stromsmoe K A.Lumped-circuit models for nonlinear inductors exhibiting hysteresis loops[J].IEEE Transactions on Circuit Theory,1970,17(4):564-574.

[4] Neves W L A,Dommet H W.On modeling iron core nonlinearities[J].IEEE Transactions on Power Systems,1993,8(5):417-425.

[5] Neves W L A,Dommet H W.Saturation curves of delta-connected transformers from measurements[J].IEEE Transactions on Power Delivery,1995,10(3):1432-1437.

[6] Greene J D,Gross C A.Nonlinear modeling of transformers[J].IEEE Transactions on Industry Applications,1988,24(3):434-438.

[7] Pedra J,Corcoles F,Sainz L,et al.Hormonic nonlinear transformer modeling[J].IEEE Transactions on Power Delivery,2004,19(2):884-890.

[8] Pedra J,Sainz L,Córcoles F,et al.Pspice computer model of a nonlinear three-legged transformer[J].IEEE Transactions on Power Delivery,2004,19(1):200-207.

[9] 白保東,趙曉璇,陳德志,等.基于J-A模型對直流偏磁條件下變壓器勵磁電流的模擬及實驗研究[J].電工技術學報,2013,28(2):162-166. Bai Baodong,Zhao Xiaoxuan,Chen Dezhi,et al.Simulation and experiment research on transformer excitation current under DC masnetic bias based on J-A model[J].Transactions of China Electrotechnical Society,2013,28(2):162-166.

[10]張炳達,黃杰.基于查分雙正交小波熵的變壓器勵磁涌流識別算法[J].電力系統保護與控制,2014,42(18):9-13. Zhang Bingda,Huang Jie.Identification of transformer inrush current based on difference biorthogonal wavelet entropy[J].Power System Protection and Control,2014,42(18):9-13.

[11]司馬文霞,鄭哲人,楊慶,等.用參數不匹配混沌系統的脈沖同步方法抑制鐵磁諧振過電壓[J].電工技術學報,2012,27(6):218-225. Sima Wenxia,Zheng Zheren,Yang Qing,et al.Suppressing the ferroresonance overvoltage by impulsive synchronization of chaotic systems with parameter mismatched[J].Transactions of China Electrotechnical Society,2012,27(6):218-225.

[12]阮全榮,王亮,施圍.一種計及損耗的變壓器非線性模型[J].西安交通大學學報,2006,40(10):1115-1119. Ruan Quanrong,Wang Liang,Shi Wei.A kind of transformer’s nonlinear model considering loss[J].Journal of Xi’an Jiaotong University,2006,40(10): 1115-1119.

[13]林濟鏗,王超,呂曉艷,等.采用統一磁路及空載試驗的變壓器飽和模型[J].高電壓技術,2011,37(2):422-428. Lin Jikeng,Wang Chao,Lü Xiaoyan,et al.Saturation model of transformer using unified magnetic equivalent circuit model and no-load experiment[J].High Voltage Engineering,2011,37(2):422-428.

[14]李謙,陳猷清,王曉瑜,等.變壓器和電磁式電壓互感器鐵心勵磁特性轉換方法的探討[J].變壓器,1993(12):10-13. Li Qian,Chen Youqing,Wang Xiaoyu,et al.Transformers and electromagnetic voltage transformer explore the core excitation characteristic conversion method[J].Transformer,1993(12):10-13.

[15]黃晶晶,徐習東,曾平,等.電力變壓器鐵心的非線性磁化特性[J].高電壓技術,2010,36(10):2582-2587. Huang Jingjing,Xu Xidong,Zeng Ping,et al.Iron core nonlinearities of power transformer[J].High Voltage Engineering,2010,36(10):2582-2587.

[16]郭滿生,梅桂華,張喜樂,等.直流偏磁條件下單相三柱電力變壓器的損耗計算[J].電工技術學報,2010,25(7):67-71. Guo Mansheng,Mei Guihua,Zhang Xile,et al.Calculation of losses in single phase three limb transformer under DC-biasing[J].Transactions of China Electrotechnical Society,2010,25(7):67-71.

[17]曹林,何金良,張波.直流偏磁狀態下電力變壓器鐵心動態磁滯損耗模型及驗證[J].中國電機工程學報,2008,28(24):141-146. Cao Lin,He Jinliang,Zhang Bo.Dynamic hysteresis loss model of power transformer under DC current biasing and its verification[J].Proceedings of the CSEE,2008,28(24):141-146.

[18]Bertotti G,Fiorillo F,Mazzetti P,et al.Satistical models of losses in soft magnttic materials[J].Journal of Magnetism and Magnetic Materials,1984,46(1):68-80.

[19]Bertotti G.General properties of power losses in soft ferromagnetic materials[J].IEEE Transactions on Magnetics,1988,24(1):621-630.

[20]李貞,李慶民,李長云,等.直流偏磁條件下變壓器的諧波畸變特征[J].電力系統保護與控制,2010,38(24):52-55. Li Zhen,Li Qingmin,Li Changyun,et al.Harmonic dictortion feature of AC transformers caused by DC bias[J].Power System Protection and Control,2010,38(24):52-55.

[21]黃純,朱智軍,曹一家,等.一種電網諧波與簡諧波分析新方法[J].電工技術學報,2013,28(9):32-39. Huang Chun,Zhu Zhijun,Cao Yijia,et al.A novel power system harmonic and interharmonic analysis method[J].Transactions of China Electrotechnical Society,2013,28(9):32-39.

[22]謝鵬,王善銘,王祥珩,等.三繞組諧波屏蔽變壓器的研究[J].電工電能新技術,2005,24(2):42-43,80. Xie Peng,Wang Shanming,Wang Xiangheng,et al.Study of harmonic suppression scheme based on three-winding transformer[J].Advanced Technology of Electrical Engineering and Energy,2005,24(2):42-43,80.

[23]李勇,羅隆福,Christian REHTANZ,等.新型自耦感應諧波屏蔽電力變壓器的通用模型及諧波模型[J].電力系統自動化,2011,35(11):43-49. Li Yong,Luo Longfu,Christian REHTANZ,et al.General model and harmonic model for a kind of new auto-inductive and harmonic-suppressing transformer[J].Automation of Electric Power Systems,2011,35(11):43-49.

[24]方少麟,李建華,黃瑩,等.變壓器諧波模型的研究與評價[J].電力系統及其自動化學報,2013,25(2):103-108. Fang Shaolin,Li Jianhua,Huang Ying,et al.Research and evaluation of the transformer harmonic model[J].Proceedings of the CSU-EPSA,2013,25(2):103-108.

[25]Densem T J,Bodger P S,Arrillaga J.Three phase transmission system modelling for harmonic penetration studies[J].IEEE Transactions on Power Apparatus and Systems,1984,103(2):310-317.

An Improved Transformer Model Considering of the Losses and Hysteresis of the Core

Li Huiqi1Li Xiaomeng1Li Jinzhong2Zhang Shuqi2Yan Yangang1

(1.The State Key Laboratory of Alternate Electrical Power System With Renewable Energy Sources North China Electric Power University Baoding 071003 China 2.China Electric Power Research Institute Beijing 100192 China)

In order to analyze the power system transient process or ferromagnetic resonance problems,a precise equivalent nonlinear model of a transformer is important.Based on no-load transformer test,an equivalent circuit model of the transformer and a calculation method of magnetizing branch parameter are proposed,which uses nonlinear resistors and parallel connection between inductance considering nonlinear and hysteresis.The equivalent circuit of a model transformer magnetizing branch circuit model is implemented through no-load test measurement in the paper.The calculated results are found to be in good agreement with the measured data.Simultaneously,this method is also used to establish a circuit model of an origin single-phase three-limb transformer,JMB500.Again,and the computed results agree well with the measured data.Finally,the equivalent circuits of transformer model are obtained at different frequencies and the comparison between the computed data and the measured data also shows good agreement,which proves the validity and accuracy of the modeling method proposed in this paper.

Losses,transformer,core,hysteresis,model

國家自然科學基金資助項目(51577066,51307057)。

2015-10-10 改稿日期2016-04-30

TM401

李慧奇 男,1970年生,副教授,博士,研究方向為電磁場數值計算和電磁兼容。

E-mail:huiqili@263.net(通信作者)

李曉孟 男,1991年生,研究生,研究方向為電磁場數值計算。

E-mail:1194253042@qq.com

猜你喜歡
變壓器模型
一半模型
理想變壓器的“三個不變”與“三個變”
重要模型『一線三等角』
重尾非線性自回歸模型自加權M-估計的漸近分布
開關電源中高頻變壓器的設計
3D打印中的模型分割與打包
一種不停電更換變壓器的帶電作業法
變壓器免維護吸濕器的開發與應用
FLUKA幾何模型到CAD幾何模型轉換方法初步研究
基于RFE-SA-SVM的變壓器故障診斷
主站蜘蛛池模板: 亚洲精品桃花岛av在线| 一级香蕉视频在线观看| 成年A级毛片| 亚洲国产日韩欧美在线| 亚洲国产精品日韩专区AV| 国产资源站| 永久免费av网站可以直接看的| 四虎永久在线精品影院| 国产欧美日韩精品第二区| 91综合色区亚洲熟妇p| 99精品热视频这里只有精品7| 99久久精品国产综合婷婷| 欧美天堂久久| 免费av一区二区三区在线| 国产精品55夜色66夜色| 亚洲第一香蕉视频| 亚洲狠狠婷婷综合久久久久| 无码区日韩专区免费系列| 国产精品视频久| 欧美在线视频a| 女人18毛片一级毛片在线| 国产激情国语对白普通话| 亚洲丝袜第一页| 国产女人在线观看| 福利姬国产精品一区在线| 国产亚洲精品91| 香蕉久久国产精品免| 国产极品美女在线| 强奷白丝美女在线观看| 国产靠逼视频| 欧美区国产区| 日韩欧美网址| 欧美劲爆第一页| 直接黄91麻豆网站| 天天操天天噜| 亚洲天堂色色人体| 亚洲最新在线| a在线亚洲男人的天堂试看| 毛片大全免费观看| 粗大猛烈进出高潮视频无码| 日本人妻一区二区三区不卡影院| 欧美一级高清视频在线播放| 日韩一二三区视频精品| 视频一本大道香蕉久在线播放| 欧美三级视频网站| 亚洲区欧美区| 久久综合国产乱子免费| 97se亚洲| 日韩在线第三页| 亚洲91在线精品| 日韩欧美国产中文| 久夜色精品国产噜噜| 5388国产亚洲欧美在线观看| 99精品视频播放| 日本久久网站| 国内黄色精品| 亚洲综合九九| 99视频国产精品| 真实国产乱子伦视频| 搞黄网站免费观看| 午夜影院a级片| 找国产毛片看| 亚洲AV人人澡人人双人| 在线看片国产| 国产精品私拍99pans大尺度| 免费A级毛片无码无遮挡| 亚洲视频免| 国产丰满大乳无码免费播放| 天堂网亚洲系列亚洲系列| 国产丝袜91| 国产成人精品综合| 日本精品中文字幕在线不卡| 91在线免费公开视频| 亚洲日韩国产精品综合在线观看| 中文字幕调教一区二区视频| 国产精品久久国产精麻豆99网站| 欧美综合成人| 最新无码专区超级碰碰碰| 亚洲日本中文字幕乱码中文| 亚洲日本一本dvd高清| 久久久久久高潮白浆| 国产视频入口|