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

    農(nóng)村電網(wǎng)故障點(diǎn)定位方法研究

    2017-12-16 07:05:32王亭亭趙玉林
    農(nóng)機(jī)化研究 2017年2期
    關(guān)鍵詞:行波暫態(tài)測(cè)距

    王亭亭,趙玉林

    (東北農(nóng)業(yè)大學(xué),哈爾濱 150030)

    ?

    農(nóng)村電網(wǎng)故障點(diǎn)定位方法研究

    王亭亭,趙玉林

    (東北農(nóng)業(yè)大學(xué),哈爾濱 150030)

    隨著用戶對(duì)供電可靠性要求的提高和農(nóng)村電網(wǎng)的改造升級(jí)的推進(jìn),有些地區(qū)電壓等級(jí)為35kV的變電站實(shí)現(xiàn)了雙電源供電系統(tǒng)。為此,在了解農(nóng)村電網(wǎng)故障特點(diǎn)和現(xiàn)有線路故障測(cè)距方法后,選擇單端行波法進(jìn)行故障定位分析。根據(jù)均勻換位輸電線路的相模變換所具有的性質(zhì),提出改進(jìn)型相模變換矩陣,以達(dá)到解除相間耦合及單一模量表達(dá)故障相的雙重目的,應(yīng)用小波變換對(duì)變換后的故障波形分析提取故障信息,求取故障距離。通過(guò)MatLab進(jìn)行模擬仿真,所得的故障距離相對(duì)誤差值均在1%之內(nèi),具有較高的可靠性,并通過(guò)GUI創(chuàng)建用戶與仿真模型之間的交互界面。

    農(nóng)村電網(wǎng);小波變換;改進(jìn)型相模變換;雙電源供電

    0 引言

    農(nóng)村電網(wǎng)供電區(qū)域大、線路長(zhǎng)、負(fù)荷密度低,因而供電電壓等級(jí)低,35kV電壓等級(jí)線路占農(nóng)村電網(wǎng)線路供電線路的比重較大。由于電壓等級(jí)較低且農(nóng)村電網(wǎng)所處環(huán)境比較惡劣,地理環(huán)境復(fù)雜導(dǎo)致頻繁出現(xiàn)故障,發(fā)生單相接地故障的概率最高,占到線路發(fā)生故障概率的80%[1-3]。由于35kV線路采用中性點(diǎn)經(jīng)消弧線圈接地運(yùn)行方式,當(dāng)系統(tǒng)發(fā)生單相接地故障時(shí),各相間電壓和相位,仍舊維持著一個(gè)電力系統(tǒng)平衡,可以暫時(shí)繼續(xù)運(yùn)行一小段時(shí)間,通常大約為2h左右;但不能長(zhǎng)期運(yùn)行,當(dāng)單相接地故障電容電流增大、長(zhǎng)時(shí)間帶故障運(yùn)行時(shí),就會(huì)使故障擴(kuò)大成兩點(diǎn)接地短路或者多點(diǎn)接地短路。同時(shí),弧光接地還會(huì)引起系統(tǒng)過(guò)電壓,從而損壞設(shè)備,破壞系統(tǒng)安全運(yùn)行,因此必須快速、準(zhǔn)確地找到故障位置并切除故障,恢復(fù)供電[2-3,7]。

    隨著用戶對(duì)電網(wǎng)供電可靠性需求的提高,近幾年很多地區(qū)實(shí)行農(nóng)網(wǎng)改造升級(jí),實(shí)現(xiàn)了35kV變電站雙電源供電。對(duì)于雙電源供電系統(tǒng),現(xiàn)在的研究大多針對(duì)于110kV及以上的高壓輸電線路,因此農(nóng)村電網(wǎng)35kV雙電源供電輸電線路故障的準(zhǔn)確定位成為一個(gè)亟待解決的問(wèn)題。

    當(dāng)前,故障測(cè)距的研究是農(nóng)電領(lǐng)域的熱點(diǎn)問(wèn)題。傳統(tǒng)的阻抗法,其實(shí)現(xiàn)方法簡(jiǎn)單、定位技術(shù)比較成熟,受過(guò)渡電阻的影響很大[2]。文獻(xiàn)[4]中的注入信號(hào)法需在變電所加裝信號(hào)注入設(shè)備,增加了系統(tǒng)的安全隱患,且該方法不能檢測(cè)到瞬時(shí)性故障和間歇性接地故障。行波法從測(cè)量端數(shù)量劃分,又可分為雙端法和單端法,雙端精度較高,但雙端故障測(cè)距需要兩端數(shù)據(jù)采樣同步,很難實(shí)現(xiàn)且需要大量資金,此方法存在一定的局限性[5]。綜合考慮,最終選擇單端行波故障測(cè)距方法,它具有幾乎不受過(guò)渡電阻和系統(tǒng)運(yùn)行方式等因素約束的優(yōu)點(diǎn)。

    在單端行波故障測(cè)距中,主要存下兩方面的的問(wèn)題:一是由于實(shí)際線路中相線與地線之間電磁耦合的存在會(huì)影響測(cè)距精度;二是行波在傳播過(guò)程中存在衰減的問(wèn)題,導(dǎo)致自動(dòng)測(cè)距的可靠性較差。為解決以上兩個(gè)問(wèn)題,本文提出了一種基于改進(jìn)型相模變換和小波變換相結(jié)合的輸電線路故障定位方法。文獻(xiàn)[6]中指出:當(dāng)發(fā)生故障時(shí),電流行波波頭具有較大的幅值且上升速度較快,其定位比故障電壓行波具有更高的靈敏度。因此,在分析時(shí)故障行波時(shí)選擇故障電流暫態(tài)行波。

    1 改進(jìn)型相模變換

    對(duì)于實(shí)際的電力系統(tǒng)來(lái)說(shuō),相線和地線組成了一個(gè)彼此之間存在電磁耦合的具有分布參數(shù)的多導(dǎo)線系統(tǒng)[2-3],當(dāng)輸電線路有一相發(fā)生故障時(shí),非故障相由于電磁耦合的存在也會(huì)產(chǎn)生相對(duì)應(yīng)的故障行波,影響測(cè)距精度。因此,為了提高輸電線路定位的精確度,需要通過(guò)相模變換將耦合去除。

    設(shè)正常運(yùn)行時(shí)測(cè)量端的電流行波采樣值用I表示,S定義為相模變換矩陣,將I的a、b、c三相的相量進(jìn)行相模變換后轉(zhuǎn)化成0、α、β模量的關(guān)系式為

    [Ia,Ib,Ic]T=S[I0,Iα,Iβ]T

    (1)

    相模變換矩陣不唯一且與線路參數(shù)無(wú)關(guān),常用的有凱倫布爾變換、克拉克變換及對(duì)稱分量等[2],但是這些相模變換矩陣沒(méi)有辦法實(shí)現(xiàn)單一模量就能表示故障類型。為此,以凱倫布爾變換為例進(jìn)行分析,并根據(jù)相模變換矩陣的特點(diǎn)推導(dǎo)構(gòu)造出改進(jìn)型相模變換矩陣。

    1.1 凱倫布爾變換

    凱倫布爾變換矩陣為

    (2)

    (3)

    將變換矩陣S-1帶入式(1)中可得

    (4)

    當(dāng)發(fā)生C相故障時(shí),邊界條件為Ia=Ib=0,帶入式(4)中可得α模量的值恒為0,即不能說(shuō)明C相發(fā)生接地故障。當(dāng)發(fā)生B相故障時(shí),邊界條件為Ia=Ic=0,帶入式(4)中可得β模量的值恒為0,即不能說(shuō)明B相發(fā)生接地故障。由此可知:凱倫布爾變換不能夠?qū)崿F(xiàn)根據(jù)單一模量來(lái)判斷故障相。

    經(jīng)過(guò)對(duì)其他的相模變換矩陣進(jìn)行分析后得知:在解決單一模量實(shí)現(xiàn)故障選相的問(wèn)題上,現(xiàn)有相模變換矩陣存在問(wèn)題。本文推導(dǎo)并構(gòu)造改進(jìn)型相模變換矩陣,以達(dá)到單模量表示故障類型的目的。

    1.2 推導(dǎo)并構(gòu)造改進(jìn)型相模變換矩陣

    設(shè)改進(jìn)型相模變換矩陣為,S,則

    (5)

    將S-1帶入式(3)中可得

    (6)

    由相模變換理論可知:0模分量以大地為回路,存在較嚴(yán)重的衰減和參數(shù)隨頻率變換較大的問(wèn)題,其幅值只有在故障處存在對(duì)地放電通道時(shí)才比較大,所以不適用于所有故障類型;而線模分量以導(dǎo)線為回路,基本上與頻率沒(méi)有什么太大的關(guān)系,可用于所有故障類型[3],因此將線模分量作為研究類型。若僅僅通過(guò)α模量或β模量就能夠直接判斷出故障類型,則只需保證下列關(guān)系,即

    (7)

    設(shè)P=a11a22a33-a11a23a32-a12a21a33+a12a23a31+a13a21a32-a13a22a31,可求相模變換矩陣通式S為

    (8)

    根據(jù)式(7)可知:若要達(dá)到改進(jìn)型相模變換的需求,則只需滿足

    (9)

    根據(jù)式(9)可構(gòu)造出改進(jìn)型相模變換矩陣為

    (10)

    (11)

    通過(guò)表1來(lái)顯示不同故障類型時(shí)的α模值和β模值。

    表1 各種故障類型下的電流量模值

    續(xù)表1

    通過(guò)表1中的表達(dá)式可知:運(yùn)用改進(jìn)型相模變換矩陣得到的線模分量的值(α模值或β模值)在各種故障類型的情況下均為非零值且各不相同,即僅僅通過(guò)α模值或β模值就能反應(yīng)判斷出所有的故障類型。

    2 小波變換

    2.1 小波分解與重構(gòu)

    在故障測(cè)距中,故障通常表現(xiàn)在信號(hào)的突變上,這些突變點(diǎn)往往包含著故障的各種信息。在測(cè)量端截取的波形無(wú)論是在何種故障類型下,故障波形中的奇異點(diǎn)包含故障發(fā)生的地點(diǎn)、方向和時(shí)間等信息,這些信息能否準(zhǔn)確地提取出來(lái)是測(cè)距的關(guān)鍵。由于小波變換所特有的特點(diǎn),使得其在對(duì)信號(hào)突變點(diǎn)進(jìn)行捕捉時(shí)具有很大的優(yōu)勢(shì)[2,8-10]。

    本文采用對(duì)故障波形進(jìn)行多尺度分解和重構(gòu)的算法,利用小波變換對(duì)其作多尺度分解,然后對(duì)信號(hào)在高頻下進(jìn)行單支重構(gòu)。單只重構(gòu)信號(hào)突變處也正是高頻信號(hào)出現(xiàn)的時(shí)刻,而模極大值代表的是信號(hào)的幅值和變化的徒度,通過(guò)截取模極大值所對(duì)應(yīng)的時(shí)刻獲取故障波形到達(dá)測(cè)量端的時(shí)間點(diǎn),進(jìn)而求得故障距離[10]。該方法不受其他阻抗不匹配點(diǎn)反射波干擾,不受故障類型的影響,在近區(qū)也不存在無(wú)法識(shí)別反射波的問(wèn)題。

    2.2 小波基函數(shù)的選擇

    不同類型的小波基波形的規(guī)則性和支撐長(zhǎng)度也有較大的差別,因此在選用不同的小波基函數(shù)在不同尺度對(duì)同一個(gè)突變信號(hào)進(jìn)行分解所得到的結(jié)果也會(huì)有很大的差別[2,10]。根據(jù)DbN小波系具有正則性、正交性和緊支集等特點(diǎn),本文選擇Db5小波來(lái)對(duì)故障波形進(jìn)行分析。

    3 故障仿真分析

    MatLab作為成熟的適用于多學(xué)科、多種工作平臺(tái)的大型軟件,以矩陣運(yùn)算為基礎(chǔ)可以實(shí)現(xiàn)工程計(jì)算、算法研究建立模型與仿真、工程繪圖等功能。其中,MatLab中包含的電力系統(tǒng)模塊庫(kù)SIMPOWER,為搭建35kV配電線路系統(tǒng)提供了很好的平臺(tái)[7]。本文應(yīng)用MatLab建立了一個(gè)線路長(zhǎng)度為50km、電壓等級(jí)35kV、工頻50Hz的雙電源供電系統(tǒng)模型。

    MatLab模型主要包括2個(gè)三相電壓源模塊及兩個(gè)配電變壓器模型,2個(gè)消弧線圈模塊(由于消弧線圈是一個(gè)具有鐵芯的可調(diào)電感電路,因此直接用電感代替),1個(gè)電壓電流測(cè)量模塊,2個(gè)三相輸電線路模塊,1個(gè)三相故障發(fā)生器,4個(gè)接地模塊,如圖1所示。

    圖1 35kV雙電源供電配電線路MATLAB模擬仿真圖

    本文研究的故障行波在輸電線路上的傳播速度取v=2.95655397×105km/s。其中,v根據(jù)線路參數(shù)近似得到,不同的線路模型對(duì)應(yīng)的線路參數(shù)也不同,因此波速的取值要根據(jù)實(shí)際情況而定[11]。

    3.1 驗(yàn)證改進(jìn)型相模變換

    以故障發(fā)生距離5km處、經(jīng)300Ω過(guò)渡電阻的B相接地故障為例,故障發(fā)生時(shí)間為0.024s,則故障暫態(tài)波形s1和經(jīng)改進(jìn)型相模變換后的暫態(tài)波形s2分別如圖2和圖3所示。

    在故障區(qū)間截取任一時(shí)刻的相模分量的幅值和對(duì)應(yīng)的線模分量的幅值都滿足2Ib=Iα、-Ib=Iβ。在進(jìn)行大量的仿真后可知:無(wú)論處于何種故障類型,相模分量與線模分量的關(guān)系都滿足表1。

    圖2 暫態(tài)波形s1

    圖3 暫態(tài)波形s2

    3.2 過(guò)渡電阻的大小對(duì)電流波形的影響

    以線路中A相發(fā)生單相接地故障為例,故障距離為35km處。在過(guò)渡電阻Rf分別為3、300、30 000Ω時(shí)的故障電流暫態(tài)仿真波形s如圖4、圖5和圖6所示。

    圖4 當(dāng)過(guò)渡電阻Rf為3Ω時(shí)的波形s

    圖5 當(dāng)過(guò)渡電阻Rf為300Ω時(shí)的波形s

    圖6 當(dāng)過(guò)渡電阻Rf為30000Ω時(shí)的波形s

    從圖4、圖5和圖6中可知:不同過(guò)渡電阻情況下截取到故障時(shí)刻電流暫態(tài)行波幅值的大小是不同的。電流暫態(tài)行波幅值的大小反應(yīng)測(cè)量端可以截取到故障波的難易程度;過(guò)渡電阻越大,測(cè)取到的故障電流沖擊值越小,突變時(shí)刻越不容易被捕捉到,越不利于測(cè)距分析。但是,實(shí)際上對(duì)于35kV線路而言,過(guò)渡電阻通常不會(huì)很大,因此單端行波法還是具有一定的可靠性。

    3.3 模擬結(jié)果分析

    前文提到本文將分析的是經(jīng)改進(jìn)型相模變換后的線模分量,選擇α模值分量。由于單相接地故障發(fā)生概率最大,以A相經(jīng)30Ω過(guò)渡電阻接地故障為例,故障位置為25km處故障波形s的α模分量波形波形如圖7所示。

    對(duì)s暫態(tài)波形的α模分量進(jìn)行小波分解與重構(gòu)得到所需的故障信息,根據(jù)故障信息求取故障距離,表2列出了不同故障距離的情況下所測(cè)取的故障距離和相對(duì)誤差值。

    圖7 s的α模量波形

    故障距離/km測(cè)量距離/m相對(duì)誤差/%55.04680.93671010.07890.78891515.09610.64112020.02470.12372525.02130.08523030.09940.33143534.96590.09744039.95650.10854544.94720.1173

    由表2可知:當(dāng)過(guò)渡電阻為30Ω時(shí),不同故障距離所測(cè)取數(shù)據(jù)相對(duì)誤差的數(shù)值均在1%之內(nèi),具有較高的精度,從而也驗(yàn)證了方法的可靠性。

    限于篇幅,也為了進(jìn)一步直觀地驗(yàn)證單端行波法的可靠性,將通過(guò)圖8和圖9分別展現(xiàn)當(dāng)處于處于同一故障類型A相接地故障而過(guò)渡電阻不同和處于同一過(guò)渡電阻300Ω而故障類型不同時(shí)所得到數(shù)據(jù)的相對(duì)誤差曲線。

    通過(guò)圖8可以直觀地看出:當(dāng)故障類型相同時(shí),過(guò)渡電阻在一定范圍內(nèi)所測(cè)取的故障距離相對(duì)誤差數(shù)值都在1%之內(nèi),有較高的精度。由圖9可知:當(dāng)過(guò)渡電阻一定時(shí)改變故障類型所測(cè)取的故障距離相對(duì)誤差數(shù)值也在1%之內(nèi),進(jìn)一步驗(yàn)證了單端行波法的可靠性。

    圖8 不同過(guò)渡電阻Rf時(shí)的相對(duì)誤差

    圖9 不同故障類型下的相對(duì)誤差

    3.4 GUI界面模擬

    為了方便人們直觀了解到故障波形及故障信息[12],本文將通過(guò)GUI界面來(lái)模擬實(shí)現(xiàn)用戶與仿真模型的交互。當(dāng)發(fā)生A相接地故障、過(guò)渡電阻為30Ω且故障距離為35km處時(shí),小電流接地雙電源供電故障仿真的界面如圖10所示。

    圖10 界面圖

    通過(guò)此界面,可以清楚直觀地觀察到線路是否發(fā)生故障,以及故障時(shí)的電流暫態(tài)行波波形、故障信息和測(cè)距結(jié)果。在改變過(guò)渡電阻、故障距離、故障類型后,都可以通過(guò)此界面便捷地觀察到電流暫態(tài)行波的變化。

    4 結(jié)論

    本文通過(guò)單端行波的方法對(duì)35kV雙電源供電系統(tǒng)輸電線路故障測(cè)距進(jìn)行研究,這種方法可以適用于各種線路和各種母線類型,具有經(jīng)濟(jì)實(shí)用的特點(diǎn)。本文對(duì)相模變換進(jìn)行了改進(jìn),通過(guò)MatLab進(jìn)行仿真模擬,并通過(guò)GUI實(shí)現(xiàn)人機(jī)界面交互。

    1)為了解除線路相線和地線之間電磁耦合,在了解現(xiàn)有的相模變換矩陣不能夠用單一模量表示故障類型的缺點(diǎn)后,結(jié)合相模變換本身具有的特點(diǎn)和實(shí)現(xiàn)單一模量表示故障類型的目的,對(duì)相模變換矩陣進(jìn)行改進(jìn)。

    2)通過(guò)MatLab對(duì)35kV雙電源供電線路進(jìn)行模擬仿真分析。在不同故障距離、不同過(guò)渡電阻(當(dāng)過(guò)渡電阻在一定范圍內(nèi))和不同故障類型下求得的故障距離相對(duì)誤差都在1%之內(nèi),具有較高的精度,證實(shí)了方法的可靠性。

    3)應(yīng)用GUI設(shè)計(jì)出35kV雙電源供電故障仿真的人機(jī)交互界面,方便人們快速地掌握故障信息、分析故障波形及故障點(diǎn),進(jìn)而排除故障恢復(fù)供電。

    隨著用戶對(duì)電網(wǎng)的可靠性要求的提高,原有的供電模型也在發(fā)生變換,因此要隨著模型變化進(jìn)行分析,以使農(nóng)村電網(wǎng)更加安全、可靠地運(yùn)行。

    [1] 郭增為.農(nóng)村電網(wǎng)單相接地故障行波定位方法的研究[D].保定:華北電力大學(xué),2009.

    [2] 徐青山.電力系統(tǒng)故障診斷與故障恢復(fù)[M].北京:中國(guó)電力出版社,2007.

    [3] 楊炎龍.配電網(wǎng)故障定位研究[D].廣州:華南理工大學(xué),2013.

    [4] 杜剛,劉迅,蘇高峰.基于FTU和“S”信號(hào)注入法的配電網(wǎng)接地故障定位技術(shù)的研究[J].電力系統(tǒng)保護(hù)與控制,2010,38(12):73-76.

    [5] 鄭州,呂艷萍.基于小波變換的雙端行波測(cè)距新方法[J].電網(wǎng)技術(shù),2010,34(1):203-207.

    [6] 趙玉林,關(guān)永剛,梁秋艷,等.農(nóng)村電網(wǎng)單相接地故障量的仿真分析[J].農(nóng)機(jī)化研究,2010,32(5):201-203,211.

    [7] 周羽生,彭琢.小電流接地系統(tǒng)單相故障的MatLab仿真[J].電力系統(tǒng)及其自動(dòng)化學(xué)報(bào),2009,21(4):50-52,99.

    [8] 梁睿,靳征.行波時(shí)頻復(fù)合分析的配電網(wǎng)故障定位研究[J].中國(guó)電機(jī)工程學(xué)報(bào), 2013,33(28):130-136.

    [9] 吳樂(lè)鵬,黃純.基于暫態(tài)小波能量的小電流接地故障選線新方法[J].電力自動(dòng)化設(shè)備,2013,5(5):70-75.

    [10] 熊小伏,林金洪.基于小波重構(gòu)的電力電纜故障測(cè)距方法[J].電網(wǎng)技術(shù),2003,27(6):36-38,70.

    [11] 申文,劉亞?wèn)|.輸電線路實(shí)際運(yùn)行狀態(tài)對(duì)行波波速的影響[J].電力系統(tǒng)及其自動(dòng)化學(xué)報(bào),2014, 26(9):12-15,29.

    [12] 李建成,邴銳.基于GUI的數(shù)學(xué)信號(hào)處理課程演示研究系統(tǒng)[J].電子科技,2015,28(4):76-79.

    Study on Rural Power Grid Failure Point Positioning Method

    Wang Tingting, Zhao Yulin

    (College of Electrical and Information of Northeast Agricultural University,Harbin 150030,China)

    With the improvement of user requirements for power supply reliability and advance of the transformation and upgrade of rural power grid, some parts of the voltage level of 35 kV substation have double electricity supply system. Based on the characteristics of rural power grid and the current fault location methods, this paper proposes a novel single-end wave approach to locate the fault occurred in the rural distribution network. According to the nature of averagely various-phase transmission line ,deduces a modified phase-magnitude transformation matrix . The approach eliminates the effects of magnetic coupling and single-magnitude representation. After phase-magnitude transforming of the temporary wave, the wavelet transform is utilized to analyze the magnitude component of the transformed wave. Then the fault information is extracted to locate the fault point. The simulation tests are performed on a 35 kV distribution line with double suppliers in the environment of MATLAB. And the the relative error of fault distance values are within 1%. Then use GUI which is the human and machine interface (HMI) achieve the communication of users and simulation models.

    rural power grid; wavelet transform; modified phase mode transform; double electricity suppliers

    2016-01-04

    黑龍江省教育廳科技成果產(chǎn)業(yè)化前期研發(fā)培育項(xiàng)目(1252CGZH30);黑龍江省教育廳科學(xué)技術(shù)研究項(xiàng)目(12511038)

    王亭亭(1990-),女,哈爾濱人,碩士研究生,(E-mail)yeting1394@163.com。

    趙玉林(1956-),男,黑龍江富錦人,教授,碩士生導(dǎo)師,(E-mail) zyl5631@163.com。

    TM727.1

    A

    1003-188X(2017)02-0226-06

    猜你喜歡
    行波暫態(tài)測(cè)距
    一類非局部擴(kuò)散的SIR模型的行波解
    300Mvar空冷隱極同步調(diào)相機(jī)暫態(tài)特性仿真分析
    類星體的精準(zhǔn)測(cè)距
    科學(xué)(2020年3期)2020-01-06 04:02:51
    電力系統(tǒng)全網(wǎng)一體化暫態(tài)仿真接口技術(shù)
    電子制作(2018年14期)2018-08-21 01:38:28
    Joseph-Egri方程行波解的分岔
    淺談超聲波測(cè)距
    電子制作(2017年7期)2017-06-05 09:36:13
    除氧器暫態(tài)計(jì)算研究
    Kolmogorov-Petrovskii-Piskunov方程和Zhiber-Shabat方程的行波解
    基于PSOC超聲測(cè)距系統(tǒng)設(shè)計(jì)
    (3+1)維Kdv-Zakharov-Kuznetsev方程的亞純行波解
    伦理电影免费视频| 丝袜人妻中文字幕| 国产一区二区在线av高清观看| 每晚都被弄得嗷嗷叫到高潮| 欧美 亚洲 国产 日韩一| 色婷婷久久久亚洲欧美| 老汉色∧v一级毛片| 美女免费视频网站| 亚洲人成77777在线视频| 亚洲精品久久国产高清桃花| 久久久久国产一级毛片高清牌| 可以在线观看毛片的网站| 国产黄片美女视频| 亚洲国产看品久久| 午夜亚洲福利在线播放| 国产欧美日韩一区二区精品| 亚洲人成电影免费在线| 亚洲欧美精品综合久久99| 成人一区二区视频在线观看| 日韩中文字幕欧美一区二区| 国产91精品成人一区二区三区| 中文字幕另类日韩欧美亚洲嫩草| 亚洲成a人片在线一区二区| xxx96com| 亚洲五月色婷婷综合| www日本黄色视频网| 亚洲av电影在线进入| 欧美乱妇无乱码| 久久久久国产一级毛片高清牌| 国产精品久久久久久人妻精品电影| 天天躁夜夜躁狠狠躁躁| 国产精品免费视频内射| av免费在线观看网站| 两个人看的免费小视频| 亚洲国产中文字幕在线视频| 亚洲天堂国产精品一区在线| 国产人伦9x9x在线观看| 亚洲av五月六月丁香网| 两个人免费观看高清视频| 日本免费a在线| 99久久国产精品久久久| 脱女人内裤的视频| 午夜福利高清视频| 欧美日韩黄片免| 久久这里只有精品19| 97超级碰碰碰精品色视频在线观看| 桃红色精品国产亚洲av| 久久精品成人免费网站| 久9热在线精品视频| 欧美色视频一区免费| 在线观看免费午夜福利视频| 国产成人欧美在线观看| 亚洲精品色激情综合| 国产一区二区在线av高清观看| 国产91精品成人一区二区三区| 一进一出好大好爽视频| 一二三四社区在线视频社区8| 少妇裸体淫交视频免费看高清 | 变态另类成人亚洲欧美熟女| 亚洲中文字幕日韩| 欧美国产精品va在线观看不卡| 色播在线永久视频| 精品国内亚洲2022精品成人| 黄色a级毛片大全视频| 国内精品久久久久精免费| 老司机深夜福利视频在线观看| 婷婷精品国产亚洲av在线| 国产黄a三级三级三级人| 制服诱惑二区| 一级毛片高清免费大全| 97碰自拍视频| 亚洲人成77777在线视频| 最近在线观看免费完整版| 中文字幕最新亚洲高清| 亚洲人成网站在线播放欧美日韩| 99精品欧美一区二区三区四区| 高清毛片免费观看视频网站| 高潮久久久久久久久久久不卡| 精品一区二区三区视频在线观看免费| 国产av不卡久久| 真人一进一出gif抽搐免费| 最近最新中文字幕大全电影3 | videosex国产| 黄色片一级片一级黄色片| 久久青草综合色| 国产一级毛片七仙女欲春2 | 香蕉久久夜色| 亚洲精品一区av在线观看| a级毛片a级免费在线| 午夜激情福利司机影院| 搞女人的毛片| 一级a爱片免费观看的视频| 51午夜福利影视在线观看| 亚洲成a人片在线一区二区| 91在线观看av| 亚洲欧美日韩无卡精品| 亚洲一码二码三码区别大吗| 国产av一区在线观看免费| 国产精品98久久久久久宅男小说| 色综合亚洲欧美另类图片| 午夜福利在线观看吧| √禁漫天堂资源中文www| 成年版毛片免费区| 国产成年人精品一区二区| 在线观看免费日韩欧美大片| 成人欧美大片| 亚洲美女黄片视频| 男女视频在线观看网站免费 | 亚洲精品一卡2卡三卡4卡5卡| avwww免费| 国产又黄又爽又无遮挡在线| 在线观看免费视频日本深夜| 黑人巨大精品欧美一区二区mp4| 久9热在线精品视频| www国产在线视频色| 日本一本二区三区精品| 窝窝影院91人妻| 成人一区二区视频在线观看| 人人澡人人妻人| 一级毛片女人18水好多| 亚洲国产精品久久男人天堂| 亚洲精品av麻豆狂野| 黑人操中国人逼视频| 亚洲三区欧美一区| 波多野结衣高清作品| 久久午夜亚洲精品久久| 久久久久久久久中文| 男人舔女人的私密视频| 自线自在国产av| 性欧美人与动物交配| 日本一本二区三区精品| 人人妻人人看人人澡| 久热爱精品视频在线9| 午夜福利欧美成人| 亚洲狠狠婷婷综合久久图片| 国产v大片淫在线免费观看| 91国产中文字幕| 最近在线观看免费完整版| 国产av一区在线观看免费| 中文字幕久久专区| 91成年电影在线观看| 99国产综合亚洲精品| 88av欧美| 熟女电影av网| 男人舔女人的私密视频| 亚洲成av片中文字幕在线观看| av天堂在线播放| 在线播放国产精品三级| 熟女少妇亚洲综合色aaa.| 国产人伦9x9x在线观看| 久久青草综合色| 色播在线永久视频| 亚洲熟女毛片儿| 91字幕亚洲| 亚洲男人天堂网一区| 后天国语完整版免费观看| 成人欧美大片| 国产又爽黄色视频| 久久精品成人免费网站| 一级毛片高清免费大全| www日本黄色视频网| 午夜福利在线观看吧| 淫妇啪啪啪对白视频| 男人舔奶头视频| 极品教师在线免费播放| 可以免费在线观看a视频的电影网站| 色婷婷久久久亚洲欧美| 免费在线观看完整版高清| 亚洲成人久久爱视频| 国产熟女xx| 男女下面进入的视频免费午夜 | 亚洲一区高清亚洲精品| 一二三四在线观看免费中文在| 日韩 欧美 亚洲 中文字幕| www.www免费av| 精品国产乱子伦一区二区三区| 日韩精品青青久久久久久| 国产伦人伦偷精品视频| 欧美国产日韩亚洲一区| 国产精品综合久久久久久久免费| 男女之事视频高清在线观看| 人妻丰满熟妇av一区二区三区| 91字幕亚洲| 男女之事视频高清在线观看| 看免费av毛片| 国内揄拍国产精品人妻在线 | 老熟妇乱子伦视频在线观看| 久久久久国产一级毛片高清牌| 精品久久久久久久久久免费视频| 91在线观看av| 日韩欧美三级三区| 欧美zozozo另类| 成人三级做爰电影| 成人一区二区视频在线观看| 久热爱精品视频在线9| 老司机靠b影院| 美女高潮喷水抽搐中文字幕| 91国产中文字幕| 欧美一级a爱片免费观看看 | netflix在线观看网站| 免费在线观看视频国产中文字幕亚洲| 男人舔奶头视频| 黄色视频不卡| 午夜福利视频1000在线观看| 国产单亲对白刺激| xxx96com| 国产午夜福利久久久久久| 亚洲片人在线观看| 日日干狠狠操夜夜爽| 99久久国产精品久久久| 日韩大尺度精品在线看网址| 国产午夜福利久久久久久| 国内揄拍国产精品人妻在线 | 欧美色视频一区免费| 精品卡一卡二卡四卡免费| 午夜福利欧美成人| 国产区一区二久久| 老汉色∧v一级毛片| 国产黄色小视频在线观看| 国产成人av激情在线播放| 亚洲国产精品sss在线观看| 给我免费播放毛片高清在线观看| 女性被躁到高潮视频| 19禁男女啪啪无遮挡网站| 美女扒开内裤让男人捅视频| 久久久国产欧美日韩av| bbb黄色大片| 老司机在亚洲福利影院| 搡老妇女老女人老熟妇| 成人精品一区二区免费| 国产成年人精品一区二区| 国产单亲对白刺激| 嫩草影视91久久| 99久久无色码亚洲精品果冻| 国产亚洲精品一区二区www| 国产一卡二卡三卡精品| 国产亚洲欧美精品永久| 久久精品91蜜桃| 国产精品久久久av美女十八| 亚洲一区中文字幕在线| 男女午夜视频在线观看| av在线天堂中文字幕| 夜夜夜夜夜久久久久| 国产精品亚洲美女久久久| 久久国产乱子伦精品免费另类| а√天堂www在线а√下载| 精品日产1卡2卡| 白带黄色成豆腐渣| 久久精品国产亚洲av香蕉五月| 午夜福利在线观看吧| 国产在线精品亚洲第一网站| 色哟哟哟哟哟哟| 宅男免费午夜| 视频在线观看一区二区三区| www.熟女人妻精品国产| 一本综合久久免费| 免费看a级黄色片| 久久香蕉国产精品| 午夜精品在线福利| 久久午夜综合久久蜜桃| 一边摸一边做爽爽视频免费| 国产精品免费一区二区三区在线| 色精品久久人妻99蜜桃| 午夜福利成人在线免费观看| 免费在线观看视频国产中文字幕亚洲| 伊人久久大香线蕉亚洲五| 久久久久久大精品| 后天国语完整版免费观看| 免费在线观看日本一区| 妹子高潮喷水视频| 岛国视频午夜一区免费看| 麻豆国产av国片精品| 可以在线观看毛片的网站| 亚洲全国av大片| av视频在线观看入口| 国产精品野战在线观看| 天堂动漫精品| 91大片在线观看| www日本黄色视频网| 99久久国产精品久久久| 国产99白浆流出| av超薄肉色丝袜交足视频| 国产亚洲精品av在线| 亚洲成av片中文字幕在线观看| 法律面前人人平等表现在哪些方面| 男女午夜视频在线观看| 亚洲欧美激情综合另类| 日韩欧美国产一区二区入口| 在线观看免费午夜福利视频| 亚洲成人免费电影在线观看| 丝袜人妻中文字幕| 亚洲色图av天堂| 99国产综合亚洲精品| 成人欧美大片| 日韩三级视频一区二区三区| 欧美最黄视频在线播放免费| 国产精品二区激情视频| 久久伊人香网站| 精品久久久久久,| 国产麻豆成人av免费视频| 在线观看www视频免费| 母亲3免费完整高清在线观看| 国产伦一二天堂av在线观看| 丝袜美腿诱惑在线| 91在线观看av| 极品教师在线免费播放| 久久久久免费精品人妻一区二区 | 变态另类成人亚洲欧美熟女| 日本成人三级电影网站| 中文字幕人妻丝袜一区二区| 女警被强在线播放| 欧美日韩精品网址| 深夜精品福利| 三级毛片av免费| 国产在线精品亚洲第一网站| 亚洲精品国产区一区二| 免费av毛片视频| 国产av又大| 最近最新中文字幕大全电影3 | 亚洲欧美日韩无卡精品| 亚洲国产精品成人综合色| 国产人伦9x9x在线观看| 老熟妇仑乱视频hdxx| 久久久久国产一级毛片高清牌| 国语自产精品视频在线第100页| 黄频高清免费视频| 亚洲va日本ⅴa欧美va伊人久久| 午夜福利免费观看在线| 黄色毛片三级朝国网站| 欧美激情极品国产一区二区三区| 国产成人av激情在线播放| 久久99热这里只有精品18| 可以免费在线观看a视频的电影网站| 欧美+亚洲+日韩+国产| 一级作爱视频免费观看| 动漫黄色视频在线观看| 欧美成狂野欧美在线观看| 亚洲中文字幕日韩| 久久精品成人免费网站| 亚洲三区欧美一区| 一区二区三区精品91| 亚洲精品一卡2卡三卡4卡5卡| 大型av网站在线播放| 无限看片的www在线观看| 欧美性猛交╳xxx乱大交人| 啦啦啦韩国在线观看视频| 午夜免费观看网址| 国产人伦9x9x在线观看| 1024香蕉在线观看| 国产又黄又爽又无遮挡在线| 美女高潮喷水抽搐中文字幕| 脱女人内裤的视频| 在线天堂中文资源库| 日本一本二区三区精品| 99久久久亚洲精品蜜臀av| 国产精品亚洲美女久久久| 亚洲精品中文字幕在线视频| 啦啦啦 在线观看视频| 亚洲免费av在线视频| 国产午夜福利久久久久久| 制服丝袜大香蕉在线| 黑人欧美特级aaaaaa片| 91成人精品电影| 国产成人影院久久av| 美国免费a级毛片| 精品国产美女av久久久久小说| 性欧美人与动物交配| 国产男靠女视频免费网站| 美女 人体艺术 gogo| 亚洲 国产 在线| 操出白浆在线播放| 日本五十路高清| 成年免费大片在线观看| 一级a爱片免费观看的视频| 亚洲第一欧美日韩一区二区三区| 一级a爱片免费观看的视频| 国产精品亚洲av一区麻豆| 国产精品影院久久| 久久精品国产99精品国产亚洲性色| 女人高潮潮喷娇喘18禁视频| 成人国语在线视频| 男女下面进入的视频免费午夜 | 777久久人妻少妇嫩草av网站| а√天堂www在线а√下载| 欧美国产精品va在线观看不卡| 日本a在线网址| www日本在线高清视频| 欧美zozozo另类| 老司机在亚洲福利影院| 中国美女看黄片| 日本成人三级电影网站| 在线观看免费日韩欧美大片| 色哟哟哟哟哟哟| 高清毛片免费观看视频网站| 好看av亚洲va欧美ⅴa在| 一级a爱视频在线免费观看| 久热爱精品视频在线9| 国产精品,欧美在线| 国产单亲对白刺激| 久久婷婷成人综合色麻豆| 最近最新中文字幕大全免费视频| 久久人妻av系列| 午夜福利欧美成人| 成年免费大片在线观看| 成人av一区二区三区在线看| 老鸭窝网址在线观看| 久久青草综合色| 在线观看一区二区三区| www国产在线视频色| 国产男靠女视频免费网站| 色精品久久人妻99蜜桃| 少妇熟女aⅴ在线视频| 88av欧美| 成人国语在线视频| 日本成人三级电影网站| 黑丝袜美女国产一区| 久久午夜综合久久蜜桃| 啦啦啦免费观看视频1| 18禁黄网站禁片午夜丰满| 国产视频一区二区在线看| 国产黄色小视频在线观看| 怎么达到女性高潮| 他把我摸到了高潮在线观看| 一本精品99久久精品77| 免费高清视频大片| 国产亚洲精品综合一区在线观看 | www.www免费av| 国产99白浆流出| 可以在线观看的亚洲视频| 婷婷精品国产亚洲av在线| 少妇的丰满在线观看| 久久国产精品男人的天堂亚洲| 校园春色视频在线观看| 久久中文看片网| 日本一本二区三区精品| 视频区欧美日本亚洲| www.熟女人妻精品国产| 一边摸一边做爽爽视频免费| 婷婷亚洲欧美| 久久精品影院6| 国产久久久一区二区三区| 色哟哟哟哟哟哟| 日本a在线网址| 日日夜夜操网爽| 久久精品aⅴ一区二区三区四区| x7x7x7水蜜桃| 欧美大码av| 亚洲中文字幕日韩| 夜夜爽天天搞| 在线十欧美十亚洲十日本专区| av视频在线观看入口| 国产精品一区二区精品视频观看| 成人手机av| 午夜两性在线视频| 69av精品久久久久久| 啦啦啦观看免费观看视频高清| 日韩 欧美 亚洲 中文字幕| 69av精品久久久久久| 一本综合久久免费| av有码第一页| 看免费av毛片| 国产精品一区二区三区四区久久 | 亚洲熟女毛片儿| 1024视频免费在线观看| 色综合欧美亚洲国产小说| 91在线观看av| 精品国产美女av久久久久小说| 日本一区二区免费在线视频| 又紧又爽又黄一区二区| 国产亚洲精品久久久久5区| 看片在线看免费视频| 亚洲男人天堂网一区| 久久久久精品国产欧美久久久| 欧美午夜高清在线| 91大片在线观看| 亚洲国产欧洲综合997久久, | 亚洲一区二区三区色噜噜| 18美女黄网站色大片免费观看| 色老头精品视频在线观看| av在线天堂中文字幕| 精品国产国语对白av| 天天添夜夜摸| 色综合婷婷激情| 美女高潮到喷水免费观看| 日韩大尺度精品在线看网址| 日韩三级视频一区二区三区| 亚洲av中文字字幕乱码综合 | 一区二区三区国产精品乱码| 美女高潮到喷水免费观看| 男男h啪啪无遮挡| 国产一区二区在线av高清观看| 亚洲第一电影网av| 男女做爰动态图高潮gif福利片| 国产v大片淫在线免费观看| 在线视频色国产色| 97超级碰碰碰精品色视频在线观看| 国产精品影院久久| 在线看三级毛片| 国产野战对白在线观看| 99热这里只有精品一区 | 一边摸一边做爽爽视频免费| 久久久精品国产亚洲av高清涩受| 国产高清激情床上av| 最新美女视频免费是黄的| 亚洲av电影在线进入| 色老头精品视频在线观看| 精品欧美国产一区二区三| 村上凉子中文字幕在线| 国产又黄又爽又无遮挡在线| 极品教师在线免费播放| 熟女电影av网| xxxwww97欧美| 亚洲第一青青草原| 99久久99久久久精品蜜桃| 亚洲国产精品成人综合色| 12—13女人毛片做爰片一| 可以免费在线观看a视频的电影网站| 欧美激情 高清一区二区三区| 看黄色毛片网站| 国产精品 国内视频| 精品电影一区二区在线| 亚洲精品国产一区二区精华液| АⅤ资源中文在线天堂| 在线十欧美十亚洲十日本专区| tocl精华| 99re在线观看精品视频| 日日夜夜操网爽| 久久久久久久精品吃奶| 日本免费一区二区三区高清不卡| 成人亚洲精品av一区二区| 国产欧美日韩精品亚洲av| 超碰成人久久| 久久久国产成人免费| 亚洲国产精品999在线| 91老司机精品| 少妇裸体淫交视频免费看高清 | 亚洲av成人不卡在线观看播放网| 丰满人妻熟妇乱又伦精品不卡| 淫秽高清视频在线观看| 亚洲欧美精品综合一区二区三区| 国产麻豆成人av免费视频| 婷婷六月久久综合丁香| 国产成人精品久久二区二区免费| 一边摸一边抽搐一进一小说| 一进一出抽搐gif免费好疼| 午夜福利成人在线免费观看| 国产三级黄色录像| 亚洲专区国产一区二区| 99热6这里只有精品| 国产精品久久久人人做人人爽| 中文字幕高清在线视频| 丝袜人妻中文字幕| 亚洲自拍偷在线| 黄片播放在线免费| 最近最新中文字幕大全免费视频| ponron亚洲| 法律面前人人平等表现在哪些方面| 色综合站精品国产| av免费在线观看网站| 亚洲人成伊人成综合网2020| 国产91精品成人一区二区三区| 女人爽到高潮嗷嗷叫在线视频| 青草久久国产| 中文字幕人成人乱码亚洲影| 脱女人内裤的视频| 老鸭窝网址在线观看| 久久久久亚洲av毛片大全| 伦理电影免费视频| 日韩欧美免费精品| 悠悠久久av| 一卡2卡三卡四卡精品乱码亚洲| 国产精品日韩av在线免费观看| 变态另类成人亚洲欧美熟女| e午夜精品久久久久久久| 88av欧美| 亚洲午夜精品一区,二区,三区| 18禁黄网站禁片午夜丰满| 一进一出好大好爽视频| 又黄又粗又硬又大视频| 久久香蕉激情| 国产精品久久久人人做人人爽| 国产亚洲精品一区二区www| 后天国语完整版免费观看| 国内揄拍国产精品人妻在线 | 亚洲精品美女久久av网站| 搞女人的毛片| 99riav亚洲国产免费| 999久久久国产精品视频| 国产三级在线视频| 午夜激情福利司机影院| ponron亚洲| 国产区一区二久久| 亚洲国产精品sss在线观看| 非洲黑人性xxxx精品又粗又长| 黑丝袜美女国产一区| 一级毛片高清免费大全| 亚洲精品色激情综合| 久久热在线av| 琪琪午夜伦伦电影理论片6080| 一本综合久久免费| 免费电影在线观看免费观看| 午夜福利在线观看吧| 在线观看66精品国产| 嫁个100分男人电影在线观看| 999久久久精品免费观看国产| 成人手机av| 每晚都被弄得嗷嗷叫到高潮| 国产日本99.免费观看| 国产高清激情床上av| 亚洲精品av麻豆狂野| 免费观看人在逋| 国产精品1区2区在线观看.| 亚洲色图 男人天堂 中文字幕| 男人舔奶头视频| 国产激情久久老熟女|