• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于有限元法的同走廊交直流線路混合電場的計算研究

    2017-03-09 03:14:03王建林李博文
    關(guān)鍵詞:交直流場強(qiáng)走廊

    王建林,靳 斌,李博文,李 興,熊 明

    (西華大學(xué) 電氣與電子信息學(xué)院,四川 成都 610039)

    基于有限元法的同走廊交直流線路混合電場的計算研究

    王建林,靳 斌,李博文,李 興,熊 明

    (西華大學(xué) 電氣與電子信息學(xué)院,四川 成都 610039)

    特高壓直流線路與超高壓交流線路同走廊架設(shè)將會有效地減少架設(shè)輸電線路的土地成本,但卻帶來了新的電場計算問題.為了研究二者的混合電場,首先使用有限元分析軟件Comsol Multiphysics計算了±800kV特高壓直流線路下方的合成電場及500kV超高壓交流輸電線路下方的工頻電場,進(jìn)而計算了混合電場,最后重點(diǎn)探討了導(dǎo)線布局對混合電場的影響.結(jié)果表明,基于Comsol Multiphysics的混合電場計算是合理可靠的;合理的接近距離和導(dǎo)線最小對地高度能分別有效地減少混合區(qū)域場強(qiáng)和混合電場最大值.

    特高壓直流線路;Comsol;同走廊;混合電場;有限元

    近年來,為滿足持續(xù)增長的用電負(fù)荷需求,需要建設(shè)更多的輸電線路.與此同時,原有的土地資源變得越來越緊缺,尤其是在經(jīng)濟(jì)發(fā)達(dá)地區(qū),新輸電走廊的獲取變得越來越困難,輸電工程中勢必會出現(xiàn)交直流線路同走廊或同塔架設(shè)的情況,雖然這樣可以節(jié)約土地資源并提高單位走廊的輸送容量,但也帶來了新的電磁環(huán)境問題[1].因此對其混合電場的計算就變得很有必要了.

    近年來,中外學(xué)者基于ANSYS、Ansoft等有限元分析軟件對輸電線路的電場分布作了不少的研究.但使用Comsol Multi- physics來研究輸電線路電場分布的還不多.文獻(xiàn)[2]基于Comsol研究了特高壓直流輸電線路的電場分布,對使用Comsol求解輸電線路的電場分布有一定的借鑒意義.本文中,筆者首先基于Comsol分別計算了特高壓直流輸電線路和超高壓交流輸電線路的電場,并在此基礎(chǔ)上計算了二者的混合電場,最后分析了不同的交直流線路接近距離和導(dǎo)線最小對地高度對混合電場分布的影響規(guī)律,以期為交直流線路同走廊建設(shè)和運(yùn)行提供參考.

    1 直流線路電場的計算

    1.1 幾何模型

    以同塔雙回同相序交流輸電線及一條雙極直流輸電線路同走廊并行為例(圖1),直流導(dǎo)線型號為6×LGJ-720/50,分裂間距為0.45m,子導(dǎo)線半徑為1.68cm,極間距22m,高度為21m.交流導(dǎo)線型號為4×LGJ-400/35,分裂間距為0.45m,弧垂最低點(diǎn)對地高度為15m.交、直流輸電線路接近距離D為70m,為簡化計算不考慮避雷線的影響.

    圖1 同走廊交、直流輸電線路的幾何位置示意圖

    1.2 直流線路電場的數(shù)學(xué)控制方程及邊界條件

    高壓直流輸電線路正常運(yùn)行時,導(dǎo)線表面電場強(qiáng)度高于電暈起始電場強(qiáng)度會產(chǎn)生電暈[3],由電暈產(chǎn)生的帶電離子所形成空間電場將會與輸電線路上的電荷產(chǎn)生的標(biāo)稱電場疊加并形成合成電場.目前,計算分析直流輸電線路的合成場強(qiáng)時主要使用的方法有:半經(jīng)驗公式法,基于Deustch假設(shè)的解析法和二維數(shù)值計算方法[4].本文采取在Comsol上建立直流輸電線路模型[5-6],并結(jié)合經(jīng)驗公式的方法來計算直流輸電線路的合成場強(qiáng).其控制方程為

    (1)

    (2)

    (3)

    (4)

    (5)

    (6)

    式中:φ為標(biāo)量電位;ρ+、ρ-分別為正、負(fù)空間電荷密度;E為合成電場的電場強(qiáng)度;J+和J-分別為正、負(fù)離子流密度;K+、K-分別為正、負(fù)離子遷移率;R為離子復(fù)合系數(shù);ε0為空氣介電常數(shù);e為電子電量.

    邊界條件:

    導(dǎo)線表面電壓為其運(yùn)行電壓φ±=±U;

    地面上的電位φ0=0.

    1.3 網(wǎng)格剖分

    在分析分裂導(dǎo)線的離子流場時通常將將分裂導(dǎo)線用等效半徑導(dǎo)線代替,等效公式為[7]

    (7)

    式中:req為等效半徑;Rf為分裂間距;n為分裂數(shù);r為子導(dǎo)線半徑.

    以本文的直流輸電線為例,采用公式(7)計算可知req為0.38m.

    考慮到精度要求和盡量減少系統(tǒng)運(yùn)算量,選擇2D場進(jìn)行處理,使用靜電穩(wěn)態(tài)求解器對模型進(jìn)行求解,所設(shè)置的求解域為半徑為50m的半圓形域(圖2).

    為了提高計算的精度,筆者截取求解域邊緣附近5m寬的半圓環(huán)狀帶作無限元處理,這樣可將被人工截斷的物理域作1 000倍放大,即實際物理寬度變?yōu)?×103m;有限元中,在進(jìn)行網(wǎng)格剖分時,理論上網(wǎng)格數(shù)目越多計算精度越高,但是過密的網(wǎng)格剖分將會占用大量的計算內(nèi)存并降低計算速度,甚至于導(dǎo)致內(nèi)存溢出[8].綜合考慮,本文設(shè)置自由剖分三角形網(wǎng)格的尺寸大小為“細(xì)化”,最大單元尺寸6.7m,最小單元尺寸0.03m, 從較小的單元到較大單元的最大單元生長率為1.3.

    圖2 網(wǎng)格剖分圖

    1.4 計算結(jié)果分析

    1.4.1 標(biāo)稱場強(qiáng)計算

    從圖3可以看出,兩極導(dǎo)線的標(biāo)稱場強(qiáng)以線路中心對稱分布,最大值出現(xiàn)在輸電線路外側(cè)為24.2kV/m,最小值出現(xiàn)在兩極線路中心為0.圖3中,每極導(dǎo)線最中間的空白部分為輸電線路的內(nèi)徑,其電場值為0,通過高斯通量定理進(jìn)行計算也可得到相同結(jié)果.

    圖3 等效直流輸電線路標(biāo)稱電場分布

    求解直流線路的電場時,一般將電場評價參考點(diǎn)設(shè)在地面上0m或者1.5m處,這兩處的電場值差別不大.考慮到混合電場的求解,在下文的合成場強(qiáng)計算中,筆者將距地1.5m處作為計算的參考點(diǎn).

    1.4.2 合成場強(qiáng)計算

    對于使用半經(jīng)驗法求解直流線路的合成場強(qiáng),法國電力科學(xué)院和美國EPRI等科研機(jī)構(gòu)都作過相應(yīng)的研究并得到了不同形式的求解公式.其中,最具代表性的是美國EPRI提出的方法.基于大量模擬試驗,EPRI認(rèn)為直流線路下的電場有兩種極限情況:一是沒有電暈時,僅有導(dǎo)線上電荷決定的靜電場(標(biāo)稱電場);二是飽和電暈時僅由空間電荷決定的電場.計算直流線路下的空間電場時,首先計算出上述兩種極限情況的電場分布,在此基礎(chǔ)上再計算未飽和電暈放電時的合成場強(qiáng).并提出了計算合成電場的經(jīng)驗公式[4,8-9]:

    (8)

    式中:E、E0、ED分別為實際線路合成電場、標(biāo)稱電場、飽和電暈電場;Ke為常數(shù),查相關(guān)曲線可知;gmax、g0分別為導(dǎo)線表面最大場強(qiáng)和起暈場強(qiáng),gmax可通過仿真計算得出,g0查表可知.

    圖4 直流輸電線路合成場強(qiáng)分布

    圖4為距地1.5m處的標(biāo)稱場強(qiáng)和合成場強(qiáng)的分布圖.從圖中可以看出,合成場強(qiáng)的最大值是標(biāo)稱場強(qiáng)的2倍以上,其最大值約為27.2kV/m.可見,空間電荷的存在對于合成場強(qiáng)的增大有明顯的作用.

    2 交流線路電場的計算

    2.1 計算分析

    目前計算交流線路工頻電場的方法主要有模擬電荷法和有限元法等.雖然交流輸電線路的工頻電場是交變電場,但由于工頻頻率很低,且本文研究范圍遠(yuǎn)小于200m,在討論工頻電場時可忽略工頻磁場的影響而將其視作準(zhǔn)靜電場處理[9-10],這就為使用ComsolMultiphysics求解交流工頻電場提供了可行性.將輸電線路視為理想的三相對稱,根據(jù)圖1的幾何布局圖,在ComsolMultiphysics上建立同塔雙回同相序的交流輸電線路模型.與上文求解直流線路電場類似,對分裂導(dǎo)線進(jìn)行等效,網(wǎng)格剖分選為“細(xì)化”,求解域為半徑50m的半圓形域,不同的是,選取的求解器為靜電瞬態(tài)求解器.在施加激勵源時,通過調(diào)用函數(shù)的方式實現(xiàn),以回路一A相輸電線路為例:

    (9)

    相較于A相輸電線路,B相輸電線路和C相線路的激勵源幅值不變,僅在相位上依次滯后120°.

    2.2 計算結(jié)果

    圖5為輸電線路電場強(qiáng)度幅值的分布圖,從圖中可看出來,電場強(qiáng)度兩側(cè)呈對稱分布,最大值出現(xiàn)在導(dǎo)線外側(cè)為514kV/m,離導(dǎo)線越遠(yuǎn)電場強(qiáng)度越小,當(dāng)距地只有1.5m左右時,中心處的電場最大值大概為后文所提的5.7kV/m.圖6為回路一的A相輸電線路的電場分布,從圖中可看出,由于此時不存在其他導(dǎo)線的屏蔽作用,電場最大值為532kV/m,距離A相輸電線路越遠(yuǎn),電場強(qiáng)度下降幅度越大.

    圖5 等效交流輸電線路電場強(qiáng)度幅值分布

    圖6 回路一的A相輸電線周圍電場分布

    圖7為對仿真結(jié)果處理后所作距地1.5m處的工頻電場分布圖.從圖中可以看出,隨著距線路中心距離的增加,電場強(qiáng)度先增加后減小,最大電場強(qiáng)度出現(xiàn)在線下距中心10m處左右,其最大值約為6.4kV/m,距離超過40m后,電場強(qiáng)度緩慢變化并接近0.

    圖7 離地1.5m處電場強(qiáng)度

    3 混合電場的計算

    3.1 計算方法

    文獻(xiàn)[11-12]指出,在求解混合電場時,由于三相交流輸電線路下方的空間各點(diǎn)電場是旋轉(zhuǎn)的矢量,且直流線路電場可以視作恒定值,計算時,選擇將交流線路工頻電場和直流線路合成場強(qiáng)直接疊加得到混合電場的方法在工程中是可行的.則交、直流線路同走廊架設(shè)時,對空間中任意一點(diǎn)的混合電場強(qiáng)度有:

    Eh=Ed+Ea

    (10)

    式中:Eh為混合電場強(qiáng)度,Ed為直流線路電場,Ea為交流線路電場,其中Ed、Ea上文已經(jīng)給出.

    需要注意的是,求解混合電場時,直流線路電場強(qiáng)度取其絕對值.

    3.2 計算結(jié)果

    以接近距離D的中心作為作圖參考點(diǎn)得到圖8.從圖中可以看出,由于場的疊加作用,相較于直流線路單獨(dú)作用時,同走廊區(qū)域中靠近交流線路側(cè)的電場有明顯的增大,但是混合電場的最大值主要還是由直流線路決定的,最大值27.9kV/m出現(xiàn)在距走廊中心約20m處,且混合電場滿足我國電力行業(yè)規(guī)定的限值.

    圖8 同走廊交、直流輸電線路混合電場分布圖

    3.3 結(jié)果驗證

    目前,對于交直流同走廊的實測數(shù)據(jù)較少,文獻(xiàn)[13]在試驗室架設(shè)兩條架空線來模擬輸電線路,給其中一條導(dǎo)線通以+50kV直流電壓,另一條則通以有效值為36kV的交流電壓,交直流導(dǎo)線距地高度h均為1m,交直流導(dǎo)線間距D為1.8m,直徑d均為3.08mm,將兩條導(dǎo)線的幾何中心作為參考點(diǎn).從表1的計算結(jié)果可以看出,本文的混合電場計算值和試驗值基本一致.在交、直流線路共走廊的地方由于忽略了其相互影響,計算值比試驗值偏大,以距中心0.6m內(nèi)的數(shù)據(jù)為例,誤差最大的計算值出現(xiàn)在0.3m處,其誤差為8.9%,屬于工程設(shè)計中允許的范圍,通過理論分析,筆者也得到了相似結(jié)論.可見本文基于有限元的交、直流輸電線路同走廊時混合電場的計算是合理可靠的,對工程設(shè)計和預(yù)測具有一定參考價值.

    表1 混合電場計算值和試驗值比較

    水平距離/m試驗值/kV·m-1本文方法/kV·m-1文獻(xiàn)[10]/kV·m-1-1.213.813.713.6-0.916.316.516.0-0.615.816.015.6-0.315.316.415.8016.818.217.20.318.019.618.50.619.120.018.30.919.519.819.11.216.516.515.9

    4 導(dǎo)線布局對混合電場的影響

    對于單獨(dú)架設(shè)的直流輸電線路或交流輸電線路,導(dǎo)線對地高度、導(dǎo)線布置方式,分裂子導(dǎo)線間距和半徑,大氣條件等都會影響其線路下方的電場強(qiáng)度.同走廊架設(shè)的交、直流線路輸電線路混合電場的影響因素很多.筆者采用上文所提方法重點(diǎn)研究了改變交、直流輸電線路的接近距離和導(dǎo)線對地高度對混合電場的影響規(guī)律.

    4.1 不同接近距離對混合電場的影響

    在上文的基礎(chǔ)上,改變交、直流輸電線的接近距離,分別計算當(dāng)接近距離D為50m和90m時的混合電場.計算結(jié)果如圖9所示.從圖中可以看出,增加或是減少二者的接近距離,場強(qiáng)最大值基本沒有變化.增加兩者的接近距離,混合區(qū)域的場強(qiáng)有一定程度減少;減少兩者的接近距離,混合區(qū)域內(nèi)的場強(qiáng)增加,但是其場強(qiáng)并未超出最大值.交、直流線路共用走廊時, 由于線路內(nèi)側(cè)不必考慮居民區(qū),此時考慮走廊寬度時的原則是:在該公用走廊內(nèi)最大地面合成場強(qiáng)未超過迭加電場的最大允許值, 走廊邊緣處場強(qiáng)的允許值與交直流線路單獨(dú)存在時相同[14].因此在滿足絕緣要求的情況下,適當(dāng)?shù)販p少接近距離,既滿足了工程設(shè)計的要求,又能減少走廊寬度,從而有效地提高土地的利用率.對于同走廊架設(shè)的特高壓直流線路和超高壓交流直流線路建議接近距離為50m~70m.

    圖9 不同接近距離時的混合電場分布

    4.2 不同對地高度對混合電場的影響

    分別改變交流輸電線路和直流輸電線路的對地高度,分析其對混合電場的影響規(guī)律.接近距離D=70m,計算結(jié)果如圖10所示,從圖中可以看出,保持直流線路的最小對地高度h=21m不變,當(dāng)交流線路的最小對地高度從15m增加到18m時,混合電場的最大值基本無變化,混合區(qū)域內(nèi)的場強(qiáng)有小幅度的減少;保持交流線路的最小對地高度h=15m不變,當(dāng)直流線路的高度從21m增加到24m時,混合區(qū)域內(nèi)的場強(qiáng)基本無變化,混合電場的最大值從27.9kV/m變?yōu)?5.5kV/m,減幅達(dá)8.6%.可見,適當(dāng)?shù)靥Ц咧绷骶€路線路的最小對地高度可以明顯的減少混合電場的最大值.考慮到輸電線路實際運(yùn)行時的電場強(qiáng)度控制原則:交流工頻電場不超過7kV/m,直流電場不超過30kV/m.對于同走廊架壓直流線路和超高壓交流直流線路,建議交流導(dǎo)線最小對地高度大于12m,直流導(dǎo)線最小對地高度大于21m.

    圖10 導(dǎo)線對地高度不同時的混合電場分布

    5 結(jié)束語

    本文基于ComsolMultiphysics分別計算了同走廊架設(shè)的直流線路和交流線路的電場分布,在此基礎(chǔ)上,使用疊加法計算了二者的混合電場.通過與相關(guān)試驗數(shù)據(jù)的對比,進(jìn)一步說明了計算的準(zhǔn)確性和可靠性;重點(diǎn)討論了交、直流線路的接近距離和導(dǎo)線最小對地高度對混合電場分布的影響.結(jié)果表明:合理的接近距離能限制混合區(qū)域場強(qiáng),合理的導(dǎo)線最小對地高能減少混合電場最大值.

    [1] 趙永生, 張文亮. 同走廊高壓交、直流輸電線路混合電場分析[J]. 高電壓技術(shù), 2014,40(3): 923-929.

    [2] 趙丹, 張彼德, 李明坤,等. 特高壓直流輸電線路合成電場計算方法及其影響因素分析[J].水電能源科學(xué), 2015,33(5):182-185.

    [3] 黎小林, 張波, 王琦,等. 特高壓直流線路的場分布影響因素分析[J]. 高壓電器, 2006,42(6): 407-409.

    [4] 劉振亞. 特高壓直流輸電工程電磁環(huán)境[M]. 北京:中國電力出版社,2009.

    [5]張玉寶.基于COMSOLMultiphysics的MEMS建模及應(yīng)用[M].北京:冶金工業(yè)出版社,2007.

    [6]ROGERW.Pryor.MultiphysicsModeling:UsingCOMSOL5andMATLAB(Engineering).ElectromagneticForcesonParallelCurrent-CarryingWires[M/CD].MercuryLearning&Information:Har/DVDedition,2015.

    [7] 周愷.高壓直流輸電線路電磁環(huán)境的計算研究[D].武漢: 華中科技大學(xué), 2007.

    [8] 李海峰, 吳冀川, 劉建波,等. 有限元網(wǎng)格剖分與網(wǎng)格質(zhì)量判定指標(biāo)[J]. 中國機(jī)械工程, 2012, 23(3): 368-377.

    [9] 國家環(huán)境保護(hù)總局.HJ/T24-1998 500kV超高壓送變電工程電磁輻射環(huán)境影響評價技術(shù)規(guī)范[S].

    [10] 彭湃, 程漢湘, 陳杏燦,等. 基于有限元的超高壓輸電線路電場的數(shù)值分析與測量[J]. 電測與儀表, 2016, 53(3): 56-61.

    [11] 吳桂芳, 袁春峰, 陸家榆,等.特高壓直流線路與交流線路同走廊時混合電磁環(huán)境的計算[J].電網(wǎng)技術(shù),2010, 34(2): 14-19.

    [12]楊勇, 陸家榆, 鞠勇. 交流線路與±800kV直流線路同走廊時的地面混合電場研究[J]. 電網(wǎng)技術(shù), 2009, 33(15): 54-59.

    [13] 魯成棟, 肖登明, 秦松林. 交直流并行輸電線路地面混合電場的計算及試驗研究[J]. 高壓電器, 2014, 50(1): 81-86,91.

    [14] 陸國慶, 何宏明, 張軍. 交直流輸電線路相鄰架設(shè)或共用走廊的探討[J]. 高電壓技術(shù), 1997, 23(4): 68-70.

    (編輯:姚佳良)

    Calculation of hybrid electric field of AC/DC transmission lines erected in a common corridor based on the finite element method

    WANG Jian-lin, JIN Bin, LI Bo-wen, LI Xing, XIONG Ming

    (School of Electrical and Electronic Information,Xihua University,Chengdu 610039,China)

    The cost of land will be effectively reduced when UHVDC and EHVAC transmission lines are erected in a common corridor, but it brings a new calculation problem of the electric field. In order to study the hybrid electric field, Comsol Multiphysics, a kind of finite element analysis software was used to calculate total electric field below ±800kV UHVDC transmission lines and power frequency electric field below EHVAC transmission lines. Then the hybrid electric field was calculated,and the influence of lines layout on hybrid electric field is discussed lastly. The results show that the hybrid electric field calculation based on Comsol Multiphysics is reasonable and reliable.The reasonable approaching distance and minimum clearance to ground can effectively reduce field strength of mixed region and maximum value of hybrid electric field, respectively.

    UHVDC transmission lines; Comsol;common corridor; hybrid electric field; finite element method

    2016-07-28

    王建林,男,1157304821@qq.com;通信作者:靳斌,男,jb123456@163.com

    1672-6197(2017)03-0046-05

    TM 711

    A

    猜你喜歡
    交直流場強(qiáng)走廊
    神奇的走廊
    求解勻強(qiáng)電場場強(qiáng)的兩種方法
    交直流混合微電網(wǎng)多時間尺度協(xié)同控制
    能源工程(2021年1期)2021-04-13 02:06:12
    場強(qiáng)與電勢辨析及應(yīng)用
    走廊上的時光
    文苑(2020年4期)2020-05-30 12:35:20
    基于K-means聚類的車-地?zé)o線通信場強(qiáng)研究
    走廊上
    LTE-R場強(qiáng)測試系統(tǒng)的實現(xiàn)
    基于兩次采樣的交直流電功率迭代算法
    電測與儀表(2016年1期)2016-04-12 00:35:08
    配網(wǎng)20kV變電站交直流一體化電源的應(yīng)用與研究
    9191精品国产免费久久| 精品电影一区二区在线| 最近在线观看免费完整版| 最好的美女福利视频网| av福利片在线观看| 国产淫片久久久久久久久 | 亚洲中文字幕日韩| 亚洲色图 男人天堂 中文字幕| 日韩欧美精品v在线| 99热这里只有精品一区 | 亚洲成人中文字幕在线播放| 亚洲avbb在线观看| АⅤ资源中文在线天堂| 国产精品乱码一区二三区的特点| 欧美zozozo另类| 两个人的视频大全免费| 国产v大片淫在线免费观看| 久久精品91无色码中文字幕| 日韩欧美一区二区三区在线观看| 91老司机精品| 国产精品精品国产色婷婷| 免费无遮挡裸体视频| 又大又爽又粗| 亚洲成人久久性| 桃色一区二区三区在线观看| 色综合亚洲欧美另类图片| 丁香六月欧美| 99国产精品一区二区蜜桃av| 午夜福利高清视频| 亚洲精品乱码久久久v下载方式 | 国产爱豆传媒在线观看| 欧美国产日韩亚洲一区| 日韩三级视频一区二区三区| 国产精品国产高清国产av| 99热6这里只有精品| 少妇丰满av| 日日夜夜操网爽| 国产成人精品无人区| 婷婷精品国产亚洲av| 亚洲成a人片在线一区二区| 亚洲av免费在线观看| 亚洲午夜理论影院| 69av精品久久久久久| ponron亚洲| 真实男女啪啪啪动态图| 午夜亚洲福利在线播放| 国产一级毛片七仙女欲春2| 色尼玛亚洲综合影院| 日本 欧美在线| av在线蜜桃| 高清在线国产一区| 亚洲国产精品成人综合色| 悠悠久久av| 欧洲精品卡2卡3卡4卡5卡区| 国产一区二区三区在线臀色熟女| 亚洲欧美精品综合一区二区三区| 午夜福利视频1000在线观看| 男女午夜视频在线观看| 男女那种视频在线观看| 亚洲aⅴ乱码一区二区在线播放| 美女 人体艺术 gogo| 亚洲国产看品久久| 可以在线观看毛片的网站| 性色avwww在线观看| 天天躁狠狠躁夜夜躁狠狠躁| netflix在线观看网站| 久久久国产成人免费| 久久久成人免费电影| 亚洲欧美日韩东京热| 成人亚洲精品av一区二区| 亚洲国产欧美一区二区综合| 日本五十路高清| 久久久久久久久中文| 可以在线观看的亚洲视频| 亚洲精品中文字幕一二三四区| 午夜免费观看网址| 夜夜看夜夜爽夜夜摸| 中文字幕久久专区| 国产亚洲av嫩草精品影院| 在线免费观看的www视频| 亚洲国产高清在线一区二区三| 国内精品久久久久精免费| 国产亚洲精品久久久com| 免费电影在线观看免费观看| 国产乱人伦免费视频| 欧美国产日韩亚洲一区| 69av精品久久久久久| 久久婷婷人人爽人人干人人爱| 国内精品一区二区在线观看| 性欧美人与动物交配| 午夜亚洲福利在线播放| 热99在线观看视频| 午夜影院日韩av| 精品一区二区三区视频在线观看免费| 黑人欧美特级aaaaaa片| 国产精品久久视频播放| www.999成人在线观看| 成年女人看的毛片在线观看| 中文字幕久久专区| a在线观看视频网站| 九色成人免费人妻av| 悠悠久久av| 搞女人的毛片| 熟女少妇亚洲综合色aaa.| 日韩有码中文字幕| 亚洲av熟女| 五月伊人婷婷丁香| 亚洲成人免费电影在线观看| 国产成人av激情在线播放| 久久九九热精品免费| 亚洲国产欧洲综合997久久,| 身体一侧抽搐| 脱女人内裤的视频| 在线观看午夜福利视频| 色综合亚洲欧美另类图片| 亚洲中文字幕一区二区三区有码在线看 | 久久久久国内视频| 亚洲精品粉嫩美女一区| 国产精品久久久久久亚洲av鲁大| 亚洲色图av天堂| 99热这里只有是精品50| 精品一区二区三区四区五区乱码| 日本免费a在线| 国产伦精品一区二区三区视频9 | 成年女人毛片免费观看观看9| 成人特级av手机在线观看| 男人舔奶头视频| 亚洲欧美日韩东京热| 久久久久久久久免费视频了| 女人被狂操c到高潮| 亚洲欧美日韩高清专用| 淫秽高清视频在线观看| 中文字幕人成人乱码亚洲影| 亚洲电影在线观看av| 亚洲自偷自拍图片 自拍| 成年女人毛片免费观看观看9| 十八禁人妻一区二区| 岛国视频午夜一区免费看| 成人av一区二区三区在线看| 九色国产91popny在线| 国产三级黄色录像| 此物有八面人人有两片| 岛国在线观看网站| 黄色丝袜av网址大全| 欧美激情在线99| 九九久久精品国产亚洲av麻豆 | 国产高清视频在线播放一区| 免费观看精品视频网站| 最新在线观看一区二区三区| 欧美激情在线99| 国产精品一区二区三区四区久久| av在线蜜桃| 一级毛片高清免费大全| 可以在线观看的亚洲视频| 嫩草影院精品99| 男人舔奶头视频| 在线观看66精品国产| 久久久久久久久中文| 99久久综合精品五月天人人| 欧美乱妇无乱码| 91av网站免费观看| 国产精品一及| 久久欧美精品欧美久久欧美| 免费无遮挡裸体视频| 可以在线观看毛片的网站| 午夜成年电影在线免费观看| 午夜激情欧美在线| 一本久久中文字幕| 亚洲无线在线观看| 国产精品av视频在线免费观看| 国产高清视频在线播放一区| 又黄又爽又免费观看的视频| 啦啦啦观看免费观看视频高清| 夜夜夜夜夜久久久久| e午夜精品久久久久久久| 这个男人来自地球电影免费观看| 黄片小视频在线播放| 日韩欧美国产一区二区入口| 一级黄色大片毛片| 国产成人精品久久二区二区免费| 伦理电影免费视频| 精品无人区乱码1区二区| 精品久久久久久久末码| 国产一区二区激情短视频| 国产欧美日韩一区二区三| 波多野结衣高清无吗| 欧美日韩黄片免| 欧美乱色亚洲激情| 亚洲人与动物交配视频| 日本在线视频免费播放| 国产精品久久久久久久电影 | 国产亚洲av高清不卡| 日本一本二区三区精品| 亚洲人与动物交配视频| 午夜精品在线福利| 在线观看午夜福利视频| 51午夜福利影视在线观看| 91麻豆精品激情在线观看国产| 女人高潮潮喷娇喘18禁视频| av天堂中文字幕网| 欧美绝顶高潮抽搐喷水| 国产激情偷乱视频一区二区| 伦理电影免费视频| 国产69精品久久久久777片 | 超碰成人久久| 制服人妻中文乱码| 婷婷丁香在线五月| 国产精品亚洲av一区麻豆| 免费观看的影片在线观看| 欧美黑人巨大hd| 琪琪午夜伦伦电影理论片6080| 国产伦精品一区二区三区视频9 | 国产精品一区二区免费欧美| 欧美成人免费av一区二区三区| 亚洲欧美激情综合另类| 麻豆av在线久日| 免费看美女性在线毛片视频| 久久99热这里只有精品18| 动漫黄色视频在线观看| 亚洲自偷自拍图片 自拍| 香蕉久久夜色| 国产伦精品一区二区三区四那| 床上黄色一级片| 国产野战对白在线观看| 国产精品电影一区二区三区| 精品久久蜜臀av无| 色精品久久人妻99蜜桃| 午夜免费观看网址| 欧美日韩精品网址| 欧美乱色亚洲激情| 久久久久久久久中文| 日本一本二区三区精品| 国产高清videossex| 成年女人毛片免费观看观看9| 操出白浆在线播放| 毛片女人毛片| 母亲3免费完整高清在线观看| 男女之事视频高清在线观看| 最近视频中文字幕2019在线8| 美女高潮喷水抽搐中文字幕| 久久婷婷人人爽人人干人人爱| av国产免费在线观看| 女人高潮潮喷娇喘18禁视频| 国产乱人视频| 色哟哟哟哟哟哟| 亚洲一区二区三区色噜噜| 久久精品91蜜桃| 在线观看免费视频日本深夜| 国产精品国产高清国产av| 国产伦一二天堂av在线观看| 亚洲av片天天在线观看| 国产欧美日韩精品一区二区| 国产综合懂色| 国产aⅴ精品一区二区三区波| 亚洲精品乱码久久久v下载方式 | 精品国产乱码久久久久久男人| 国产一区二区激情短视频| 久久精品影院6| 国产黄a三级三级三级人| АⅤ资源中文在线天堂| 制服丝袜大香蕉在线| 啪啪无遮挡十八禁网站| 国产黄色小视频在线观看| 1024香蕉在线观看| 国内揄拍国产精品人妻在线| 国产精品久久久av美女十八| 波多野结衣高清无吗| 我的老师免费观看完整版| 欧美日韩综合久久久久久 | 18美女黄网站色大片免费观看| 国产黄色小视频在线观看| 亚洲狠狠婷婷综合久久图片| a级毛片a级免费在线| 久久久久久国产a免费观看| 亚洲精品粉嫩美女一区| av视频在线观看入口| 国产真实乱freesex| 中文字幕最新亚洲高清| 精品一区二区三区av网在线观看| 99久久精品一区二区三区| 一级黄色大片毛片| 黑人操中国人逼视频| 好看av亚洲va欧美ⅴa在| 色综合婷婷激情| 亚洲av片天天在线观看| 欧美乱色亚洲激情| 久久精品91蜜桃| 亚洲成人久久性| 999精品在线视频| 日本免费一区二区三区高清不卡| 欧美一区二区精品小视频在线| 国产精品综合久久久久久久免费| 欧美一区二区精品小视频在线| 亚洲 国产 在线| 人人妻,人人澡人人爽秒播| 婷婷精品国产亚洲av在线| 午夜福利欧美成人| 最新中文字幕久久久久 | 免费av毛片视频| 国产伦精品一区二区三区视频9 | 久久人妻av系列| x7x7x7水蜜桃| 性色avwww在线观看| 久久精品91无色码中文字幕| 国模一区二区三区四区视频 | 日日夜夜操网爽| 最近最新中文字幕大全电影3| 精品久久久久久久久久免费视频| 国产亚洲欧美98| 最新在线观看一区二区三区| 亚洲人与动物交配视频| 成人av在线播放网站| 1024手机看黄色片| 两性夫妻黄色片| 亚洲熟妇中文字幕五十中出| 搡老岳熟女国产| 色噜噜av男人的天堂激情| 久久天堂一区二区三区四区| 色综合亚洲欧美另类图片| 黄片小视频在线播放| 亚洲,欧美精品.| 黄频高清免费视频| 桃红色精品国产亚洲av| 久久精品91无色码中文字幕| 夜夜爽天天搞| 亚洲精华国产精华精| а√天堂www在线а√下载| 亚洲av熟女| 最近最新中文字幕大全免费视频| 亚洲七黄色美女视频| 国模一区二区三区四区视频 | 亚洲欧美日韩高清在线视频| 一级毛片精品| 国产激情偷乱视频一区二区| 国产激情欧美一区二区| 99久久精品一区二区三区| 午夜精品一区二区三区免费看| 精品久久久久久久久久久久久| 国产成人av激情在线播放| 亚洲av中文字字幕乱码综合| 国产伦精品一区二区三区四那| 看黄色毛片网站| 欧美乱码精品一区二区三区| 国产熟女xx| 亚洲最大成人中文| 全区人妻精品视频| 久久精品91蜜桃| 久久久久久人人人人人| 国产不卡一卡二| 1024香蕉在线观看| 久久久国产精品麻豆| 极品教师在线免费播放| 久久久久久久精品吃奶| 国产亚洲av嫩草精品影院| 午夜精品久久久久久毛片777| 国产成人啪精品午夜网站| 精品久久久久久久久久免费视频| 免费无遮挡裸体视频| 午夜两性在线视频| 99久久久亚洲精品蜜臀av| 97人妻精品一区二区三区麻豆| 久久亚洲精品不卡| 国产精品一及| 精品久久久久久久毛片微露脸| 日本免费a在线| 精品久久久久久久末码| 欧美日韩福利视频一区二区| 首页视频小说图片口味搜索| 亚洲欧洲精品一区二区精品久久久| 啦啦啦观看免费观看视频高清| 久久精品国产综合久久久| 黄色日韩在线| 久久久久久九九精品二区国产| av国产免费在线观看| 宅男免费午夜| 最近视频中文字幕2019在线8| 亚洲五月天丁香| bbb黄色大片| 特大巨黑吊av在线直播| 精品国内亚洲2022精品成人| 国产伦在线观看视频一区| 亚洲欧美精品综合一区二区三区| 午夜福利高清视频| 男女下面进入的视频免费午夜| 男女床上黄色一级片免费看| 噜噜噜噜噜久久久久久91| 色综合亚洲欧美另类图片| 亚洲18禁久久av| 国产精品一区二区三区四区久久| 欧美一区二区国产精品久久精品| 国产黄色小视频在线观看| 免费观看的影片在线观看| 国产美女午夜福利| 蜜桃久久精品国产亚洲av| 真实男女啪啪啪动态图| 国产黄片美女视频| 亚洲国产欧美人成| 熟女人妻精品中文字幕| 日本三级黄在线观看| 免费av不卡在线播放| 久久99热这里只有精品18| 香蕉国产在线看| 香蕉国产在线看| 桃色一区二区三区在线观看| svipshipincom国产片| netflix在线观看网站| 香蕉久久夜色| 国产亚洲精品久久久久久毛片| 一区二区三区高清视频在线| 亚洲国产精品999在线| 九九久久精品国产亚洲av麻豆 | 中文字幕熟女人妻在线| a级毛片a级免费在线| 老鸭窝网址在线观看| 亚洲欧美日韩高清专用| 99精品久久久久人妻精品| 色视频www国产| 中文字幕最新亚洲高清| 亚洲中文日韩欧美视频| 国语自产精品视频在线第100页| 成年女人毛片免费观看观看9| 亚洲真实伦在线观看| 亚洲国产中文字幕在线视频| 亚洲国产精品999在线| 熟女人妻精品中文字幕| 精品一区二区三区视频在线 | 亚洲av免费在线观看| 国产av不卡久久| 色精品久久人妻99蜜桃| 日本黄色片子视频| 久久久精品大字幕| 黄色日韩在线| 亚洲欧美日韩无卡精品| 长腿黑丝高跟| 亚洲精华国产精华精| 老司机午夜十八禁免费视频| 真实男女啪啪啪动态图| 亚洲va日本ⅴa欧美va伊人久久| 精品国产亚洲在线| av国产免费在线观看| 波多野结衣高清作品| 成人午夜高清在线视频| 精品久久久久久久末码| 午夜激情福利司机影院| 制服人妻中文乱码| 日韩精品中文字幕看吧| 1000部很黄的大片| 国产精品久久久av美女十八| 国产成年人精品一区二区| 亚洲精品中文字幕一二三四区| 久久亚洲真实| 国内揄拍国产精品人妻在线| 国产高清有码在线观看视频| 三级毛片av免费| 怎么达到女性高潮| 精品不卡国产一区二区三区| 久久精品国产亚洲av香蕉五月| 免费看美女性在线毛片视频| 一进一出抽搐动态| 窝窝影院91人妻| 色哟哟哟哟哟哟| 日韩欧美三级三区| 久久久久久久精品吃奶| www国产在线视频色| 在线观看舔阴道视频| 久久欧美精品欧美久久欧美| 亚洲中文av在线| 欧美成狂野欧美在线观看| 丰满人妻熟妇乱又伦精品不卡| 国产伦在线观看视频一区| 国产精品综合久久久久久久免费| 丰满人妻一区二区三区视频av | 久久精品人妻少妇| av在线天堂中文字幕| 一卡2卡三卡四卡精品乱码亚洲| 激情在线观看视频在线高清| 久久精品夜夜夜夜夜久久蜜豆| 丝袜人妻中文字幕| 精品国产乱子伦一区二区三区| 每晚都被弄得嗷嗷叫到高潮| 午夜亚洲福利在线播放| 国内久久婷婷六月综合欲色啪| 又黄又粗又硬又大视频| 女人高潮潮喷娇喘18禁视频| 国产高清视频在线观看网站| 性色avwww在线观看| 亚洲男人的天堂狠狠| 人人妻人人看人人澡| 国产三级中文精品| 亚洲精品国产精品久久久不卡| 禁无遮挡网站| 国产av一区在线观看免费| 手机成人av网站| 村上凉子中文字幕在线| 亚洲无线在线观看| svipshipincom国产片| 悠悠久久av| 欧美午夜高清在线| 嫩草影院入口| 可以在线观看的亚洲视频| 人妻久久中文字幕网| 一级毛片高清免费大全| 两个人的视频大全免费| 精品乱码久久久久久99久播| 欧美zozozo另类| 午夜成年电影在线免费观看| 亚洲av五月六月丁香网| 啦啦啦免费观看视频1| 日本a在线网址| 美女高潮的动态| 欧美绝顶高潮抽搐喷水| 午夜免费成人在线视频| 亚洲av日韩精品久久久久久密| 热99re8久久精品国产| 国产亚洲精品综合一区在线观看| 日韩人妻高清精品专区| 亚洲激情在线av| 欧美xxxx黑人xx丫x性爽| 亚洲男人的天堂狠狠| 国语自产精品视频在线第100页| a级毛片a级免费在线| 99热6这里只有精品| 在线观看舔阴道视频| 亚洲美女黄片视频| 国产精品女同一区二区软件 | 99热只有精品国产| 国模一区二区三区四区视频 | 男人舔女人下体高潮全视频| a级毛片在线看网站| 欧美中文日本在线观看视频| 在线十欧美十亚洲十日本专区| 国产野战对白在线观看| 亚洲黑人精品在线| 欧洲精品卡2卡3卡4卡5卡区| 一二三四在线观看免费中文在| 久久久水蜜桃国产精品网| 99精品久久久久人妻精品| 看片在线看免费视频| 91字幕亚洲| 舔av片在线| 久久久国产成人精品二区| 成年版毛片免费区| 丰满人妻一区二区三区视频av | 欧美在线黄色| 啪啪无遮挡十八禁网站| 午夜精品一区二区三区免费看| av在线蜜桃| 听说在线观看完整版免费高清| 国产单亲对白刺激| 一级黄色大片毛片| 一二三四社区在线视频社区8| 久久久久性生活片| 欧美激情久久久久久爽电影| 他把我摸到了高潮在线观看| 麻豆国产av国片精品| 中文字幕高清在线视频| 午夜福利18| 亚洲精华国产精华精| 国产 一区 欧美 日韩| 精品电影一区二区在线| 叶爱在线成人免费视频播放| 亚洲第一电影网av| 999精品在线视频| 国产精品精品国产色婷婷| 成人三级做爰电影| 免费av不卡在线播放| 国产成人福利小说| 国产视频内射| 国产精品国产高清国产av| 久久中文看片网| 18禁黄网站禁片午夜丰满| 免费在线观看视频国产中文字幕亚洲| 国产私拍福利视频在线观看| 在线观看免费午夜福利视频| 亚洲国产日韩欧美精品在线观看 | 欧美精品啪啪一区二区三区| 婷婷六月久久综合丁香| 国产精华一区二区三区| 少妇的逼水好多| 国产成人系列免费观看| 一本久久中文字幕| 搡老熟女国产l中国老女人| 国产精品av视频在线免费观看| 精品一区二区三区视频在线 | 波多野结衣巨乳人妻| 91av网一区二区| 国产真实乱freesex| 又黄又爽又免费观看的视频| 丰满人妻一区二区三区视频av | avwww免费| 亚洲欧美日韩高清专用| 搞女人的毛片| 99热只有精品国产| 97超级碰碰碰精品色视频在线观看| 曰老女人黄片| 男女床上黄色一级片免费看| www国产在线视频色| 伊人久久大香线蕉亚洲五| 操出白浆在线播放| 午夜两性在线视频| 日韩三级视频一区二区三区| av片东京热男人的天堂| 免费av不卡在线播放| 亚洲国产欧美人成| 亚洲精品乱码久久久v下载方式 | 波多野结衣巨乳人妻| 美女大奶头视频| 午夜精品久久久久久毛片777| 亚洲午夜精品一区,二区,三区| 亚洲性夜色夜夜综合| 亚洲欧美日韩高清在线视频| 国内揄拍国产精品人妻在线| 在线看三级毛片| 欧美日韩黄片免| 日本 欧美在线|