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

    基于目標級聯(lián)分析法的輸電網(wǎng)結(jié)構(gòu)優(yōu)化

    2023-02-27 07:03:48張玉敏吉興全楊子震
    電力自動化設(shè)備 2023年2期
    關(guān)鍵詞:輸電網(wǎng)下層約束

    張 旋,張玉敏,吉興全,楊子震,劉 健,楊 建

    (1.山東科技大學 電氣與自動化工程學院,山東 青島 266590;2.國網(wǎng)青島供電公司,山東 青島 266002)

    0 引言

    在能源廣域互聯(lián)的背景下,電網(wǎng)互聯(lián)有利于電能遠距離傳輸,提高能源利用效率,進而為大容量、高效能機組的優(yōu)化管理提供支持。然而,我國的能源與負荷呈逆向分布,源荷不平衡矛盾較為突出,而各類間歇性可再生能源發(fā)電占比日益增加,源網(wǎng)矛盾逐漸突出[1],輸電網(wǎng)運行過程中時常出現(xiàn)網(wǎng)絡(luò)阻塞等問題[2],使多區(qū)域互聯(lián)電網(wǎng)間的調(diào)控能力面臨嚴峻考驗,因此增強電力系統(tǒng)區(qū)域間的互聯(lián)和非同調(diào)效應具有重要意義。

    輸電網(wǎng)中的非同調(diào)現(xiàn)象,在不額外增加輸電網(wǎng)建設(shè)投資的前提下,可以根據(jù)電網(wǎng)的實際工作條件,通過開斷線路和電磁環(huán)網(wǎng)開環(huán)運行進行輸電網(wǎng)結(jié)構(gòu)優(yōu)化,可在一定程度上消除網(wǎng)絡(luò)阻塞,從而緩解了源、網(wǎng)、荷之間的沖突,提高了系統(tǒng)的靈活性和經(jīng)濟性?;诖?,輸電網(wǎng)結(jié)構(gòu)優(yōu)化(optimal transmission switching,OTS)成為一種通過改變電力系統(tǒng)中部分線路的開閉狀態(tài)以改變潮流分布的有效方法[3]。目前已有部分學者對輸電網(wǎng)的結(jié)構(gòu)優(yōu)化進行了研究。文獻[4]提出了考慮結(jié)構(gòu)優(yōu)化的含可再生能源的OTS模型,通過OTS降低了系統(tǒng)運行成本及故障后網(wǎng)絡(luò)阻塞造成的損失成本。文獻[5]結(jié)合可再生能源,構(gòu)建了計及網(wǎng)絡(luò)結(jié)構(gòu)優(yōu)化和儲能配置的兩階段隨機優(yōu)化模型,提高了電力系統(tǒng)的經(jīng)濟性。文獻[6]提出了基于廣義短路比靈敏度分析的OTS方法,其有助于協(xié)調(diào)系統(tǒng)強度與短路電流水平之間的矛盾。文獻[7]在矩估計理論的基礎(chǔ)上構(gòu)建OTS模型,使得特高壓接入下的電力系統(tǒng)安全性有所提高。由于考慮OTS的經(jīng)濟調(diào)度模型是典型的混合整數(shù)規(guī)劃問題,模型中引入了大量的離散變量,使得模型求解難度增加,上述研究雖然驗證了OTS能夠有效提升系統(tǒng)的經(jīng)濟性與可靠性,但沒有對提升模型的求解速度進行研究。

    為此,諸多學者圍繞模型的求解算法展開了深入研究。文獻[8]通過使用交替方向乘子法(alternating direction method of multipliers,ADMM)實現(xiàn)交流系統(tǒng)的分布式最優(yōu)潮流(optimal power flow,OPF)計算。由于ADMM源自增廣拉格朗日乘子法和鄰近點算法,它并沒有消除一階算法固有的缺點[9]。文獻[10]采用分布式內(nèi)點法對電網(wǎng)OPF進行分散式求解。為實現(xiàn)各個區(qū)域的同步迭代,文獻[11]采用同步交替方向乘子法解決多子系統(tǒng)的并行協(xié)調(diào)優(yōu)化問題,進而確定整體最優(yōu)解。上述優(yōu)化算法能有效提升模型的求解速度,但未考慮到電網(wǎng)規(guī)模龐大、受分區(qū)管理、信息不共享等因素影響,難以建立起統(tǒng)一的互聯(lián)電網(wǎng)模型。

    目標級聯(lián)分析(analytical target cascading,ATC)法可用于加速難以用集中方式解決大規(guī)模優(yōu)化問題的處理過程,并用于管理具有多個獨立控制實體的系統(tǒng)。文獻[12]利用ATC法尋找日前調(diào)度的最優(yōu)發(fā)電計劃,在每次迭代中,只有邊界母線的電壓和相角需要在主問題和子問題之間共享。因此,每次迭代的數(shù)據(jù)交換量低且通信簡單。文獻[13]利用ATC法搭建了雙層調(diào)度框架,提出了一種面向主動配電網(wǎng)和虛擬微電網(wǎng)的線性化交互式調(diào)度模型,解決了由于上下層之間存在交互變量而導致的耦合問題。上述研究為本文開展輸電網(wǎng)結(jié)構(gòu)優(yōu)化提供了良好的借鑒作用。

    為提升互聯(lián)電網(wǎng)整體運行水平和模型的計算效率,本文基于輸電系統(tǒng)中的非同調(diào)現(xiàn)象,提出了一種基于ATC法的計及可再生能源的OTS模型。針對區(qū)域互聯(lián)的輸電網(wǎng)結(jié)構(gòu)模型,通過母線撕裂法對互聯(lián)電網(wǎng)進行切分,在切分后的各分區(qū)內(nèi)進行獨立的經(jīng)濟調(diào)度,利用ATC法進行解耦和求解,有效降低了求解難度?;谒憷治鲵炞C了所提模型和方法的有效性。

    1 ATC法

    ATC法首先將上層系統(tǒng)中的耦合變量傳遞給子系統(tǒng),然后對各個子系統(tǒng)分別進行求解,從而獲得全局 最 優(yōu) 解[14]。ATC法 與ADMM和 輔 助 問 題 原 理(auxiliary problem principle,APP)法類似[15],APP、ADMM應用了對偶概念并引入了一組罰函數(shù),然后將原優(yōu)化問題分解成若干個子問題;ATC法將整個系統(tǒng)分解成若干個子系統(tǒng),然后應用約束松弛的概念對各問題進行并行求解[16]。本文選取區(qū)域間聯(lián)絡(luò)線交換功率和電壓相角作為共享變量,通過設(shè)置目標變量和響應變量對相鄰層間的共享變量進行建模,形成各獨立區(qū)域相關(guān)的目標函數(shù)和約束條件,其中目標變量為上層(系統(tǒng)層)傳送至下層(子系統(tǒng)層)的共享變量,響應變量為下層傳送至上層的共享變量,具體如式(1)所示。

    在式(1)所示的目標函數(shù)中,第2、3項是由子系統(tǒng)響應和耦合變量表示的偏差,子系統(tǒng)可表示為式(2)。

    式中:gsub,i(·)和hsub,i(·)分別為子系統(tǒng)i的不等式約束和等式約束集合;xsub,i為子系統(tǒng)i的設(shè)計變量。

    2 基于ATC法的輸電網(wǎng)結(jié)構(gòu)優(yōu)化模型

    考慮到可再生能源與電網(wǎng)結(jié)構(gòu)優(yōu)化的經(jīng)濟調(diào)度是一個復雜的非線性優(yōu)化問題,基于母線撕裂法和ATC法構(gòu)建OTS數(shù)學模型。為不失一般性,本文通過母線撕裂法將系統(tǒng)分成區(qū)域A—C這3個部分,分別建立相應的OTS模型。本節(jié)以區(qū)域A為例,構(gòu)建基于ATC法的輸電網(wǎng)多區(qū)域經(jīng)濟調(diào)度模型,同理可得區(qū)域B、C的數(shù)學模型。

    2.1 目標函數(shù)

    目標函數(shù)為系統(tǒng)常規(guī)機組的發(fā)電成本最小化,如式(3)所示。

    2.2 約束條件

    1)發(fā)電機出力約束。

    2)爬坡約束。

    3)支路功率約束。4)電壓相角約束。

    5)節(jié)點功率平衡約束。

    6)區(qū)域耦合約束。

    式中:ΔA為區(qū)域A的相鄰區(qū)域集合;ΓA,B為區(qū)域A和B的聯(lián)絡(luò)支路集合;(i,j)∈ΓA,B表示區(qū)域A內(nèi)節(jié)點i和區(qū)域B內(nèi)節(jié)點j分別為聯(lián)絡(luò)支路的首、末節(jié)點;Pij和θij分別為聯(lián)絡(luò)支路的功率和相角差。

    7)支路開斷數(shù)量約束。

    式中:JA為區(qū)域A中最大允許的線路開斷數(shù)目。

    8)風電和光伏出力約束。

    2.3 模型簡化

    以圖1所示的上、下兩層系統(tǒng)結(jié)構(gòu)為例,對本文所提模型進行分解。圖中:上層系統(tǒng)為區(qū)域A,下層系統(tǒng)為區(qū)域B和區(qū)域C;t為從上層系統(tǒng)發(fā)送到下層系統(tǒng)的共享變量組成的向量,將其作為目標向量;r為下層系統(tǒng)向上層系統(tǒng)傳輸?shù)墓蚕碜兞拷M成的向量,將其作為響應變量。由圖1(a)可見,為追求系統(tǒng)總運行成本的最小化,將上下層系統(tǒng)決策量進行迭代優(yōu)化,上層系統(tǒng)向下分流,下層系統(tǒng)向上反饋,從而實現(xiàn)互聯(lián)電網(wǎng)的整體優(yōu)化。為便于描述ATC算法的求解過程,本文給出模型的矩陣表達式為:

    圖1 上、下層系統(tǒng)的分解Fig.1 Decomposition of upper- and lower-layer system

    式中:gA、hA分別為上層系統(tǒng)中區(qū)域A需滿足的不等式和等式約束;x為上層系統(tǒng)中除聯(lián)絡(luò)支路以外的決策變量,x、t1、t2、…、tn為上層系統(tǒng)的區(qū)域變量;、分別為第n個下層系統(tǒng)(例如圖1(a)區(qū)域B)的不等式、等式約束條件;y為第n個下層系統(tǒng)的變量需滿足的本地約束,y、r1、r2、…、rn、t′1、t′2、…、t′n為第n個下層系統(tǒng)區(qū)域變量;c表示上層系統(tǒng)區(qū)域變量與第n個下層系統(tǒng)區(qū)域變量的耦合約束。

    3 模型轉(zhuǎn)化與求解

    3.1 模型解耦

    為便于形成各區(qū)域相關(guān)的目標函數(shù)及約束條件,設(shè)置2組變量集合。根據(jù)ATC解算思想,將約束式(19)分別在上層和下層系統(tǒng)的模型中進行求解。這里利用罰函數(shù)ζ來松弛耦合約束,如式(20)所示。

    式中:⊙表示Hadamard積;λ和μ分別為罰函數(shù)ζ的一次項和二次項的乘子向量。

    上下層系統(tǒng)的優(yōu)化模型可表示為:

    通過式(20)將上、下層系統(tǒng)間的耦合約束進行松弛處理后,上層系統(tǒng)和下層系統(tǒng)只剩下區(qū)域變量中需要滿足的本地約束和本地決策變量,由此上、下兩層系統(tǒng)實現(xiàn)了解耦。

    3.2 解耦的上層系統(tǒng)與下層系統(tǒng)優(yōu)化模型

    在基于ATC法的互聯(lián)電網(wǎng)結(jié)構(gòu)優(yōu)化調(diào)度中,上層系統(tǒng)的目標函數(shù)為:

    式中:λn、μn分別為拉格朗日一次項、二次項的乘子,通過不斷更新拉格朗日乘子使目標函數(shù)滿足收斂條件,以至趨于最優(yōu);為rn優(yōu)化后的新變量。

    上層的區(qū)域A在求解自身模型時,將虛擬負荷tn優(yōu)化后的值以參數(shù)的形式傳遞給下層系統(tǒng)。上層系統(tǒng)在更新時,需要將虛擬負荷聯(lián)動優(yōu)化,并對每個區(qū)域的虛擬發(fā)電機進行優(yōu)化。

    由式(20)可見,與上層系統(tǒng)相同,在下層系統(tǒng)進行獨立優(yōu)化時,需要將虛擬發(fā)電機rn和虛擬負荷tn結(jié)合進行優(yōu)化,通過引入罰函數(shù)對耦合約束進行松弛,并加入到下層系統(tǒng)的目標函數(shù)中。以區(qū)域B為例,下層系統(tǒng)的目標函數(shù)表示為:

    通過不斷更新乘子使目標函數(shù)滿足收斂條件,以趨于最優(yōu),從而上、下層系統(tǒng)可獨立進行更新、求解。

    3.3 收斂判據(jù)與乘子更新原則

    優(yōu)化調(diào)度算法的收斂判據(jù)為:式(25)表示在第k次迭代中,下層系統(tǒng)虛擬發(fā)電機rnk與上層系統(tǒng)虛擬負荷tnk作為耦合變量,其差值需要滿足所要求的精度ε1;式(26)表示分解后的上層系統(tǒng)和下層系統(tǒng)的總體效益在相鄰2次迭代中需要滿足精度ε2。若不能同時滿足式(25)和式(26),則依據(jù)式(27)更新乘子。

    式中:β的作用是為了加快收斂速度,一般取值為2≤β≤3;λ和μ的初值一般選取較小的常數(shù),本文取值為λ=μ=1。

    3.4 基于ATC法的OTS優(yōu)化求解流程

    基于ATC法的OTS優(yōu)化求解流程如附錄A圖A1所示,具體步驟如下。

    1)設(shè)定ATC法的最大迭代次數(shù)為k′,并置當前迭代次數(shù)k=1,輸入常規(guī)機組參數(shù),設(shè)置各個優(yōu)化變量以及罰函數(shù)乘子等初值。

    2)求解下層系統(tǒng)。根據(jù)下層系統(tǒng)模型,即式(24)和式(18)對每個區(qū)域進行并行優(yōu)化,并將求解得到的功率傳遞給上層系統(tǒng)。

    3)求解上層系統(tǒng)。根據(jù)上層系統(tǒng)模型,即式(23)和式(17)進行優(yōu)化,并將求解得到的功率,傳遞給下層系統(tǒng)。

    4)判斷是否同時滿足式(25)和式(26),若同時滿足,則終止迭代過程并輸出最優(yōu)調(diào)度結(jié)果,否則根據(jù)式(27)更新乘子,置迭代次數(shù)k=k+1并返回步驟2)。如果k>k′,則說明算法不收斂,終止計算。

    4 算例分析

    選擇改進的IEEE 14節(jié)點系統(tǒng)及IEEE 118節(jié)點測試系統(tǒng)對模型進行仿真分析。計算機配置為intel i7-6300處理器、8 GB內(nèi)存,采用GAMS軟件進行編程,并選用CPLEX工具對模型進行求解。

    4.1 IEEE 14節(jié)點系統(tǒng)

    IEEE 14節(jié)點系統(tǒng)包含5臺發(fā)電機、20條線路,機組和線路參數(shù)見附錄A表A1—A3,風電場參數(shù)和光伏發(fā)電參數(shù)分別見附錄A表A4和表A5??紤]到光伏與風力發(fā)電具有較大的間歇性與隨機性,故選取其出力上、下限中的隨機變量進行測試。通過母線撕裂法將該系統(tǒng)分解成三區(qū)域系統(tǒng),如圖2所示。設(shè)置區(qū)域A和區(qū)域B的開斷輸電線數(shù)為1。

    圖2 三區(qū)域互聯(lián)系統(tǒng)Fig.2 Three regional interconnection system

    4.1.1 ATC法有效性分析

    為了驗證所采用ATC法的求解優(yōu)勢,將ATC、APP法所得結(jié)果進行對比,如表1與圖3所示。

    表1 不同算法的求解結(jié)果Table 1 Results of different algorithms

    圖3 算法結(jié)果對比Fig.3 Comparison of algorithm results

    由表1和圖3可知:在運行成本方面,采用APP法求解的運行成本為 $ 18 212.23,而ATC法求得的運行成本為 $ 18 394.46,由于初始參數(shù)的設(shè)置與網(wǎng)絡(luò)拓撲的改變,使得ATC法的運行成本略高于APP法;在收斂性能上,ATC法的迭代次數(shù)和迭代時間分別為34次和7.71 s,APP法在迭代次數(shù)和收斂時間上的表現(xiàn)劣于ATC法,說明了本文所采用的ATC法收斂性能較好。APP法的收斂性能取決于其參數(shù)和電網(wǎng)分區(qū),由于OTS問題的復雜性,很難得到充分而有效的解;而ATC法可擴展性強,對系統(tǒng)的OTS優(yōu)化問題具有較好的適用性。

    4.1.2 輸電網(wǎng)結(jié)構(gòu)優(yōu)化及可再生能源接入影響分析

    為了證明OTS的效果和可再生能源發(fā)電對系統(tǒng)運行的影響,設(shè)置如下4種方案并進行對比分析:方案1,不考慮OTS及可再生能源發(fā)電;方案2,不考慮OTS,考慮可再生能源發(fā)電;方案3,考慮OTS,不考慮可再生能源發(fā)電;方案4,同時考慮OTS及可再生能源發(fā)電。

    1)不考慮輸電網(wǎng)結(jié)構(gòu)優(yōu)化的影響分析。

    當不考慮OTS時,為了驗證可再生能源接入對輸電網(wǎng)運行的影響,對比分析了方案1、2下各區(qū)域運行成本和機組出力,結(jié)果如表2所示。

    表2 方案1和方案2的運行結(jié)果對比Table 2 Operating result comparison between Scheme 1 and Scheme 2

    由表2分析可知:在方案2中,由于考慮了可再生能源發(fā)電的接入,極大地釋放了機組GA1的發(fā)電能力,使得經(jīng)濟性較好的機組GA1出力增加,經(jīng)濟性最差的機組GA2出力減少;采用方案1時,各區(qū)域成本分別為 $ 8 476.62、$ 7 363.38及 $ 4 363.69,而采用方案2時,各區(qū)域運行成本與采用方案1時相比分別降低了19.97 %、1.39 %和-0.13 %,其中區(qū)域A成本降低幅度較大,而區(qū)域B、C的成本變化幅度較小,這說明下層區(qū)域B、C能夠憑借自身的發(fā)電資源實現(xiàn)分布自治,從而極大地減少了上層區(qū)域A的功率救濟,使得區(qū)域A的經(jīng)濟性得到明顯提高。

    2)輸電網(wǎng)結(jié)構(gòu)優(yōu)化影響分析。

    為了證明可再生能源接入時考慮OTS的效果,方案2和方案4的機組出力和運行成本如表3所示,輸電線路運行情況如表4所示。

    表4 考慮OTS的輸電線路運行情況Table 4 Transmission line operation condition considering OTS

    由表3分析可知:與不考慮OTS的方案2相比,采用方案4時,機組GA1出力增加了43.27 %,經(jīng)濟性最差的機組GA2出力減少了35.27 %,這說明考慮OTS可實現(xiàn)源出力模式與網(wǎng)架結(jié)構(gòu)的配合,使得發(fā)電資源合理分配;在經(jīng)濟性上,采用方案4時,各區(qū)域運行成本分別為 $ 5 814.34、$ 6 086.17和 $ 4 261.58,與采用方案2時相比分別下降了14.29 %、16.18 % 和2.47 %,由此說明將電網(wǎng)結(jié)構(gòu)視為“動態(tài)可變”的經(jīng)濟調(diào)度方式能夠達到緩解網(wǎng)絡(luò)阻塞的目的,提升了系統(tǒng)的運行經(jīng)濟性。

    表3 方案2和方案4的運行結(jié)果對比Table 3 Operating result comparison between Scheme 2 and Scheme 4

    由表4分析可知:當考慮OTS時,系統(tǒng)可通過負荷模式對輸電線路狀態(tài)進行優(yōu)化調(diào)節(jié),調(diào)節(jié)過程中區(qū)域A內(nèi)支路1-101處于停運狀態(tài),從而降低了傳輸阻塞程度,輸電線路1-2的傳輸功率得以大幅度增加,進一步表明了風、光的接入及開斷輸電線路的措施同時應用于電網(wǎng)調(diào)度中,可有效提升系統(tǒng)運行的經(jīng)濟性,更大程度地實現(xiàn)源網(wǎng)協(xié)同調(diào)度;方案3和方案4的開斷線路分別為2-101、202-2和1-101、202-2,這說明在考慮OTS后,輸電線路的運行狀態(tài)可以依據(jù)相應負荷模式進行適當調(diào)整,可有效提高電網(wǎng)運行的靈活性。

    4.2 IEEE 118節(jié)點系統(tǒng)

    IEEE 118節(jié)點測試系統(tǒng)結(jié)構(gòu)圖如附錄B圖B1所示,系統(tǒng)發(fā)電機組特性數(shù)據(jù)、輸電線路數(shù)據(jù)以及系統(tǒng)負荷數(shù)據(jù)見文獻[18]。

    4.2.1 算法性能對比

    為證明ATC法在大規(guī)模系統(tǒng)中的適用性,在IEEE 118節(jié)點測試系統(tǒng)中比較了ATC法和APP法,結(jié)果如圖4所示。

    圖4 不同收斂精度的算法優(yōu)化結(jié)果Fig.4 Optimization results of algorithms with different convergence accuracy

    由圖4可知:隨著收斂精度的提高,ATC法趨于最優(yōu)的速度最快,展現(xiàn)了良好的收斂性能;收斂精度的提高導致迭代次數(shù)也隨之增加;相比于APP法,在任一精度要求下,ATC法始終保持良好的收斂性能和優(yōu)化結(jié)果,這說明ATC法對于大規(guī)模系統(tǒng)仍具有較好的適用性。

    4.2.2 優(yōu)化結(jié)果對比

    針對IEEE 118節(jié)點系統(tǒng),分別采用上述4種方案對模型進行測試分析,所得的結(jié)果如表5所示,具體分析如下。

    表5 4種方案的運行成本對比Table 5 Comparison of operating costs of four cases

    方案4的總成本相較于方案1—3分別下降了18.99 %、17.65 %、2.50 %;三區(qū)域運行成本分別為$ 47 311.34、$ 63 384.13、$ 52 117.87,與方案1—3相比都有不同程度的降低。這說明在保證系統(tǒng)安全運行的前提下,通過電網(wǎng)結(jié)構(gòu)調(diào)度,提升了可再生能源的接納水平,使系統(tǒng)運行的經(jīng)濟性得以提高。方案4中區(qū)域A與區(qū)域B的運行成本較方案1—3的下降程度高于區(qū)域C,這是由于區(qū)域A、B阻塞較嚴重,故節(jié)省成本占比較為明顯。綜上所述,將風電、光伏系統(tǒng)及開斷輸電線路同時應用于電網(wǎng)中,實現(xiàn)系統(tǒng)區(qū)域間的協(xié)同合作,能夠降低發(fā)電成本,提高能源利用效率,從而提升系統(tǒng)整體運行的靈活性和經(jīng)濟性。

    5 結(jié)論

    本文對多區(qū)域電力系統(tǒng)調(diào)度優(yōu)化問題進行研究,提出了一種基于ATC法的輸電網(wǎng)結(jié)構(gòu)優(yōu)化模型。算例分析所得結(jié)論如下。

    1)所提模型能夠有效提高互聯(lián)電網(wǎng)整體運行的經(jīng)濟性,使電網(wǎng)靈活運行的能力得以充分發(fā)揮,高效利用全網(wǎng)資源,實現(xiàn)互聯(lián)電網(wǎng)整體的統(tǒng)一優(yōu)化。

    2)所提模型增加了表示輸電線路是否運行的離散變量,使其能夠合理決策出系統(tǒng)所需的輸電網(wǎng)架構(gòu)。將電網(wǎng)結(jié)構(gòu)優(yōu)化、可再生能源與ATC相結(jié)合,極大地釋放了網(wǎng)架結(jié)構(gòu)的傳輸能力,提高了模型和方法的通用性。

    3)本文所采用的ATC法收斂速度快,應用于多區(qū)域電網(wǎng)分布式求解時具有明顯的效率優(yōu)勢,對較大規(guī)模的系統(tǒng)也具有較好的適應性。

    后續(xù)研究可考慮如何提高模型的計算效率與計算精度。

    附錄見本刊網(wǎng)絡(luò)版(http://www.epae.cn)。

    猜你喜歡
    輸電網(wǎng)下層約束
    “碳中和”約束下的路徑選擇
    約束離散KP方程族的完全Virasoro對稱
    一類多個下層的雙層規(guī)劃問題
    積雪
    陜西橫山羅圪臺村元代壁畫墓發(fā)掘簡報
    考古與文物(2016年5期)2016-12-21 06:28:48
    計及多重不確定因素的輸電網(wǎng)隨機潮流計算
    含光伏電站的輸電網(wǎng)不對稱故障分析方法
    電測與儀表(2016年4期)2016-04-11 09:44:34
    基于差分和聲搜索算法的輸電網(wǎng)差異化規(guī)劃
    適當放手能讓孩子更好地自我約束
    人生十六七(2015年6期)2015-02-28 13:08:38
    有借有還
    777米奇影视久久| 欧美97在线视频| 欧美日本中文国产一区发布| 国产熟女午夜一区二区三区 | 精品视频人人做人人爽| 久久精品人人爽人人爽视色| 九九久久精品国产亚洲av麻豆| 多毛熟女@视频| 美女cb高潮喷水在线观看| 中文字幕免费在线视频6| 久久av网站| 亚洲激情五月婷婷啪啪| 欧美另类一区| 一区二区三区四区激情视频| 欧美xxxx性猛交bbbb| 国产日韩欧美视频二区| 99久久人妻综合| 亚洲精品美女久久av网站| 一区二区三区精品91| 黑人猛操日本美女一级片| a级毛片黄视频| 亚洲av成人精品一区久久| videosex国产| 我的女老师完整版在线观看| 在现免费观看毛片| 日韩一区二区视频免费看| 午夜福利在线观看免费完整高清在| 亚洲国产精品专区欧美| 亚洲国产最新在线播放| 九九爱精品视频在线观看| 国产在线免费精品| 亚洲欧洲国产日韩| 国模一区二区三区四区视频| 国产一区有黄有色的免费视频| 欧美日韩综合久久久久久| av免费观看日本| 国产成人av激情在线播放 | 制服丝袜香蕉在线| 午夜久久久在线观看| 国产精品嫩草影院av在线观看| 日本91视频免费播放| 亚洲人成77777在线视频| 日韩伦理黄色片| 国产成人av激情在线播放 | 天天躁夜夜躁狠狠久久av| 777米奇影视久久| 五月天丁香电影| 又黄又爽又刺激的免费视频.| 色94色欧美一区二区| 99热国产这里只有精品6| 精品少妇久久久久久888优播| 蜜桃国产av成人99| 久久精品国产亚洲av涩爱| 国产黄频视频在线观看| 日日爽夜夜爽网站| 国产精品蜜桃在线观看| 精品国产一区二区三区久久久樱花| 欧美日韩精品成人综合77777| 校园人妻丝袜中文字幕| 99久久人妻综合| 大香蕉久久网| 成年美女黄网站色视频大全免费 | av电影中文网址| 啦啦啦啦在线视频资源| 国产精品一区www在线观看| 涩涩av久久男人的天堂| 人妻 亚洲 视频| 大片免费播放器 马上看| 97精品久久久久久久久久精品| 欧美97在线视频| 久久精品国产亚洲av涩爱| 亚洲精品久久久久久婷婷小说| 男女边吃奶边做爰视频| 国模一区二区三区四区视频| 欧美xxxx性猛交bbbb| 涩涩av久久男人的天堂| 亚洲,欧美,日韩| 国产成人精品在线电影| av.在线天堂| 久久久亚洲精品成人影院| 亚洲激情五月婷婷啪啪| 国产国语露脸激情在线看| 女人精品久久久久毛片| 欧美另类一区| 嫩草影院入口| 国产视频内射| 嘟嘟电影网在线观看| av黄色大香蕉| 久久精品熟女亚洲av麻豆精品| 另类亚洲欧美激情| 久久久久久久久久人人人人人人| 久久免费观看电影| 草草在线视频免费看| 日韩av免费高清视频| 美女内射精品一级片tv| 狂野欧美激情性bbbbbb| 欧美成人午夜免费资源| 18禁观看日本| 国产精品人妻久久久久久| 亚洲精品aⅴ在线观看| 免费久久久久久久精品成人欧美视频 | 亚洲精品,欧美精品| 丰满迷人的少妇在线观看| 久久精品久久久久久久性| 伊人久久精品亚洲午夜| 亚洲国产精品专区欧美| 午夜激情福利司机影院| 亚洲久久久国产精品| 五月伊人婷婷丁香| 能在线免费看毛片的网站| 又黄又爽又刺激的免费视频.| 成人毛片60女人毛片免费| 午夜久久久在线观看| a级毛色黄片| 国产探花极品一区二区| 色婷婷av一区二区三区视频| 黄色视频在线播放观看不卡| 成年人免费黄色播放视频| 精品一区二区免费观看| 九色亚洲精品在线播放| 伊人久久精品亚洲午夜| 大码成人一级视频| 麻豆精品久久久久久蜜桃| 亚洲在久久综合| 国产伦精品一区二区三区视频9| 亚洲欧美清纯卡通| kizo精华| 又大又黄又爽视频免费| 国产精品一区www在线观看| a级片在线免费高清观看视频| 日本色播在线视频| 国模一区二区三区四区视频| 亚洲精品国产av蜜桃| 日韩,欧美,国产一区二区三区| 视频在线观看一区二区三区| 国国产精品蜜臀av免费| 久久久久久人妻| 99国产综合亚洲精品| av又黄又爽大尺度在线免费看| 水蜜桃什么品种好| 国产免费现黄频在线看| 免费av中文字幕在线| 国产无遮挡羞羞视频在线观看| 亚洲中文av在线| 97超视频在线观看视频| 国产69精品久久久久777片| 女人久久www免费人成看片| 我的女老师完整版在线观看| 青春草国产在线视频| 岛国毛片在线播放| 免费久久久久久久精品成人欧美视频 | av电影中文网址| 一级毛片aaaaaa免费看小| 国产黄频视频在线观看| 我要看黄色一级片免费的| 一区二区日韩欧美中文字幕 | 菩萨蛮人人尽说江南好唐韦庄| 美女cb高潮喷水在线观看| 国产午夜精品一二区理论片| 国产一区二区在线观看av| 一本大道久久a久久精品| 日韩欧美一区视频在线观看| 99热国产这里只有精品6| 国内精品宾馆在线| 不卡视频在线观看欧美| 又粗又硬又长又爽又黄的视频| 少妇熟女欧美另类| 你懂的网址亚洲精品在线观看| 精品久久国产蜜桃| 99久久人妻综合| 成人毛片a级毛片在线播放| 久久久精品94久久精品| 亚洲少妇的诱惑av| 久久免费观看电影| 日韩欧美精品免费久久| 伦精品一区二区三区| 九九久久精品国产亚洲av麻豆| 丝袜脚勾引网站| 欧美激情国产日韩精品一区| 日韩av免费高清视频| 嘟嘟电影网在线观看| 肉色欧美久久久久久久蜜桃| 久久青草综合色| av在线老鸭窝| 国产熟女午夜一区二区三区 | 日本与韩国留学比较| 亚洲婷婷狠狠爱综合网| 日韩一区二区视频免费看| 最黄视频免费看| av不卡在线播放| 国产毛片在线视频| 中文精品一卡2卡3卡4更新| 亚洲av综合色区一区| 插逼视频在线观看| 久久97久久精品| 欧美97在线视频| 哪个播放器可以免费观看大片| 免费观看的影片在线观看| 亚洲精华国产精华液的使用体验| 国产一区亚洲一区在线观看| 亚洲精品乱码久久久v下载方式| 国语对白做爰xxxⅹ性视频网站| 亚洲三级黄色毛片| 精品久久久久久电影网| 国产精品国产三级国产专区5o| 女人久久www免费人成看片| 国产精品免费大片| 99久久人妻综合| 亚洲精品视频女| 王馨瑶露胸无遮挡在线观看| 在线精品无人区一区二区三| 九色成人免费人妻av| 欧美日韩亚洲高清精品| 国产永久视频网站| 性色avwww在线观看| 欧美老熟妇乱子伦牲交| 色婷婷av一区二区三区视频| 91精品国产九色| 纵有疾风起免费观看全集完整版| 精品卡一卡二卡四卡免费| 欧美变态另类bdsm刘玥| 亚洲精品乱久久久久久| 亚洲不卡免费看| 亚洲精品自拍成人| 日韩亚洲欧美综合| 国产成人精品福利久久| 亚洲av电影在线观看一区二区三区| 精品亚洲成a人片在线观看| 最新的欧美精品一区二区| 最后的刺客免费高清国语| a级片在线免费高清观看视频| 最近最新中文字幕免费大全7| 欧美丝袜亚洲另类| 日韩av不卡免费在线播放| 热re99久久国产66热| 国产不卡av网站在线观看| 日本vs欧美在线观看视频| a级毛色黄片| 欧美xxⅹ黑人| videossex国产| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲久久久国产精品| 精品一区二区免费观看| 水蜜桃什么品种好| 欧美精品高潮呻吟av久久| 男人爽女人下面视频在线观看| 另类精品久久| kizo精华| 哪个播放器可以免费观看大片| 亚洲情色 制服丝袜| 搡老乐熟女国产| 国产精品蜜桃在线观看| 国产有黄有色有爽视频| 男女国产视频网站| 欧美亚洲日本最大视频资源| 免费高清在线观看日韩| www.色视频.com| 色婷婷久久久亚洲欧美| 久热这里只有精品99| 视频在线观看一区二区三区| 天美传媒精品一区二区| 女的被弄到高潮叫床怎么办| 伊人久久精品亚洲午夜| 日韩制服骚丝袜av| 久久国内精品自在自线图片| 久久久欧美国产精品| 国产成人精品无人区| 午夜激情久久久久久久| a 毛片基地| 精品视频人人做人人爽| av网站免费在线观看视频| 国产日韩一区二区三区精品不卡 | 又黄又爽又刺激的免费视频.| 寂寞人妻少妇视频99o| 国产精品欧美亚洲77777| 亚洲精品成人av观看孕妇| 国产成人freesex在线| 成人国产麻豆网| 日韩免费高清中文字幕av| 欧美一级a爱片免费观看看| 国产黄色免费在线视频| 成人漫画全彩无遮挡| 晚上一个人看的免费电影| 久久精品人人爽人人爽视色| 国产精品国产三级国产av玫瑰| 欧美 亚洲 国产 日韩一| 国产成人精品一,二区| 日韩精品有码人妻一区| 在线亚洲精品国产二区图片欧美 | 在线亚洲精品国产二区图片欧美 | 十八禁高潮呻吟视频| 亚洲色图 男人天堂 中文字幕 | 18禁观看日本| 看非洲黑人一级黄片| 日韩一区二区三区影片| 日韩av不卡免费在线播放| 在现免费观看毛片| 亚洲欧美一区二区三区黑人 | xxxhd国产人妻xxx| 伦精品一区二区三区| 激情五月婷婷亚洲| 80岁老熟妇乱子伦牲交| 免费久久久久久久精品成人欧美视频 | 久久精品国产亚洲网站| 亚洲一区二区三区欧美精品| 最近2019中文字幕mv第一页| 91久久精品电影网| 曰老女人黄片| 22中文网久久字幕| 熟女人妻精品中文字幕| 成年美女黄网站色视频大全免费 | 日本黄色日本黄色录像| 国产一区二区三区av在线| 街头女战士在线观看网站| 亚洲高清免费不卡视频| 纯流量卡能插随身wifi吗| 久久久久国产网址| 国产欧美日韩综合在线一区二区| 免费av不卡在线播放| 精品一区在线观看国产| 精品人妻一区二区三区麻豆| 免费人妻精品一区二区三区视频| 免费日韩欧美在线观看| 精品酒店卫生间| 亚洲一级一片aⅴ在线观看| 啦啦啦啦在线视频资源| 国产女主播在线喷水免费视频网站| 一二三四中文在线观看免费高清| 热99久久久久精品小说推荐| 日韩中文字幕视频在线看片| 欧美日韩av久久| 日韩,欧美,国产一区二区三区| 五月玫瑰六月丁香| 国产69精品久久久久777片| 九九在线视频观看精品| 国产精品一区二区在线不卡| 亚洲不卡免费看| 亚洲精品自拍成人| 久久午夜福利片| 国产不卡av网站在线观看| www.色视频.com| 国产精品久久久久久精品电影小说| 伦理电影大哥的女人| 久久综合国产亚洲精品| 国产免费又黄又爽又色| 国产熟女午夜一区二区三区 | 一级毛片 在线播放| 最近手机中文字幕大全| 午夜免费男女啪啪视频观看| 亚洲怡红院男人天堂| 免费黄色在线免费观看| 久久国产精品大桥未久av| 男女边吃奶边做爰视频| 国产深夜福利视频在线观看| 视频在线观看一区二区三区| 精品国产国语对白av| 国产又色又爽无遮挡免| 国产国语露脸激情在线看| 日韩 亚洲 欧美在线| 免费看不卡的av| 极品人妻少妇av视频| 在线免费观看不下载黄p国产| 午夜免费男女啪啪视频观看| 少妇人妻久久综合中文| 国产精品久久久久久久久免| 中文天堂在线官网| 91精品国产国语对白视频| 久久久欧美国产精品| 精品久久蜜臀av无| 91精品国产九色| 91午夜精品亚洲一区二区三区| 日韩免费高清中文字幕av| 亚洲国产精品999| 成人国产av品久久久| 亚洲欧洲日产国产| 亚洲不卡免费看| av播播在线观看一区| 最后的刺客免费高清国语| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 建设人人有责人人尽责人人享有的| 国产有黄有色有爽视频| 日韩成人av中文字幕在线观看| 精品一区二区三卡| 亚洲国产av新网站| 十八禁高潮呻吟视频| 99热国产这里只有精品6| 久久99热6这里只有精品| 男女免费视频国产| 国产精品免费大片| 丰满迷人的少妇在线观看| 久久精品国产亚洲av天美| 久久免费观看电影| 精品久久久精品久久久| 人人妻人人爽人人添夜夜欢视频| xxx大片免费视频| 久久久久久久大尺度免费视频| av网站免费在线观看视频| 王馨瑶露胸无遮挡在线观看| 欧美激情 高清一区二区三区| 国产亚洲欧美精品永久| 国产深夜福利视频在线观看| 日本色播在线视频| 99热6这里只有精品| 日本黄大片高清| 国产免费一区二区三区四区乱码| 国产成人一区二区在线| 国产男人的电影天堂91| 久久久久久久久久久丰满| 日本av免费视频播放| 精品人妻熟女av久视频| 免费高清在线观看日韩| 在线观看美女被高潮喷水网站| 日韩一区二区视频免费看| 中国三级夫妇交换| 国产爽快片一区二区三区| 91国产中文字幕| 国产片特级美女逼逼视频| 99精国产麻豆久久婷婷| 成年人午夜在线观看视频| 美女大奶头黄色视频| 少妇 在线观看| 免费高清在线观看视频在线观看| 欧美日韩亚洲高清精品| 我的女老师完整版在线观看| 久久这里有精品视频免费| 一级a做视频免费观看| 亚洲精品一二三| 97在线视频观看| 国产色爽女视频免费观看| 桃花免费在线播放| av专区在线播放| 久久久久视频综合| 亚洲av在线观看美女高潮| 日本vs欧美在线观看视频| 在线观看人妻少妇| 日韩制服骚丝袜av| 国产爽快片一区二区三区| 成人亚洲欧美一区二区av| 老司机影院毛片| 一级片'在线观看视频| 精品99又大又爽又粗少妇毛片| 亚洲av国产av综合av卡| 99九九在线精品视频| 亚洲久久久国产精品| 国产亚洲最大av| 国产伦理片在线播放av一区| 中国三级夫妇交换| 精品人妻在线不人妻| 欧美性感艳星| 亚洲人成77777在线视频| 日本欧美国产在线视频| 蜜臀久久99精品久久宅男| 亚洲精品中文字幕在线视频| 亚洲av.av天堂| 亚洲丝袜综合中文字幕| 亚洲国产精品一区三区| 免费看光身美女| 大片免费播放器 马上看| 亚洲欧美一区二区三区黑人 | 日韩强制内射视频| 国产一区二区三区综合在线观看 | 欧美精品高潮呻吟av久久| √禁漫天堂资源中文www| 国产毛片在线视频| 99九九在线精品视频| 纵有疾风起免费观看全集完整版| 国产无遮挡羞羞视频在线观看| 亚洲怡红院男人天堂| 国产色婷婷99| 欧美日韩视频高清一区二区三区二| videosex国产| 下体分泌物呈黄色| 亚洲色图 男人天堂 中文字幕 | 国产av一区二区精品久久| 9色porny在线观看| 夫妻性生交免费视频一级片| 亚洲国产精品成人久久小说| 国产一区二区三区av在线| 国产乱人偷精品视频| 国产成人精品婷婷| 欧美激情国产日韩精品一区| 日韩免费高清中文字幕av| 亚洲精品视频女| 亚洲精品国产av蜜桃| 精品久久蜜臀av无| 国产成人一区二区在线| 精品久久久精品久久久| 婷婷色麻豆天堂久久| 久久久久久人妻| 亚洲综合色网址| 伦精品一区二区三区| 日韩欧美精品免费久久| 婷婷色麻豆天堂久久| 自线自在国产av| 国产亚洲欧美精品永久| 国产在线视频一区二区| 2021少妇久久久久久久久久久| 蜜臀久久99精品久久宅男| 久久精品国产亚洲av天美| 有码 亚洲区| 伊人久久精品亚洲午夜| 大香蕉97超碰在线| 这个男人来自地球电影免费观看 | 国国产精品蜜臀av免费| 亚洲中文av在线| 亚洲国产精品国产精品| 一级毛片电影观看| 亚洲av二区三区四区| 国产精品一区www在线观看| 五月伊人婷婷丁香| 男人添女人高潮全过程视频| 最近中文字幕高清免费大全6| 狠狠婷婷综合久久久久久88av| 97精品久久久久久久久久精品| 国产极品粉嫩免费观看在线 | 成人国产av品久久久| 搡老乐熟女国产| 色94色欧美一区二区| 午夜激情av网站| 日韩伦理黄色片| 免费人成在线观看视频色| 国产成人av激情在线播放 | 啦啦啦视频在线资源免费观看| 日本免费在线观看一区| 丰满饥渴人妻一区二区三| 国产欧美日韩综合在线一区二区| 在线观看免费视频网站a站| 日韩成人av中文字幕在线观看| 肉色欧美久久久久久久蜜桃| 另类精品久久| 777米奇影视久久| 伊人久久精品亚洲午夜| 精品一区二区三区视频在线| 99热国产这里只有精品6| 在线观看免费视频网站a站| 黄色毛片三级朝国网站| a级毛片在线看网站| 内地一区二区视频在线| 男女国产视频网站| 大陆偷拍与自拍| 99热这里只有是精品在线观看| 18禁动态无遮挡网站| 亚洲综合色网址| 人妻制服诱惑在线中文字幕| 丝袜脚勾引网站| 免费播放大片免费观看视频在线观看| 一区二区日韩欧美中文字幕 | 国产成人aa在线观看| 亚洲美女视频黄频| 蜜桃在线观看..| av电影中文网址| 成人手机av| 少妇精品久久久久久久| 色视频在线一区二区三区| 男女免费视频国产| 在线亚洲精品国产二区图片欧美 | 大香蕉久久网| 久久精品夜色国产| 亚洲av男天堂| 91精品三级在线观看| 中文乱码字字幕精品一区二区三区| 最黄视频免费看| 久久久精品免费免费高清| 毛片一级片免费看久久久久| 亚洲精品视频女| 91成人精品电影| 丰满迷人的少妇在线观看| 欧美三级亚洲精品| 曰老女人黄片| 超色免费av| 在线亚洲精品国产二区图片欧美 | 午夜日本视频在线| 亚洲精品第二区| 亚洲精品,欧美精品| 嘟嘟电影网在线观看| 久久久久精品久久久久真实原创| 久久精品国产a三级三级三级| 亚洲一级一片aⅴ在线观看| 一级毛片aaaaaa免费看小| 成人18禁高潮啪啪吃奶动态图 | 亚洲欧美成人综合另类久久久| 91成人精品电影| 黑人猛操日本美女一级片| 在线精品无人区一区二区三| 精品少妇内射三级| 亚洲欧美精品自产自拍| 狂野欧美白嫩少妇大欣赏| 国产精品无大码| 内地一区二区视频在线| 大码成人一级视频| 国产不卡av网站在线观看| 一二三四中文在线观看免费高清| 夜夜看夜夜爽夜夜摸| 久久久久国产精品人妻一区二区| 一区在线观看完整版| 国产日韩欧美亚洲二区| 久久人人爽人人爽人人片va| 亚洲精品一二三| 嫩草影院入口| 亚洲av日韩在线播放| 亚洲精品第二区| 黑人高潮一二区| 看十八女毛片水多多多| 国产乱来视频区| 国产亚洲欧美精品永久| 亚洲人成77777在线视频| 超碰97精品在线观看| 男女边吃奶边做爰视频| 高清午夜精品一区二区三区| 嘟嘟电影网在线观看| 男女高潮啪啪啪动态图| 免费观看在线日韩| 久久久久久人妻| 黄片播放在线免费| 美女大奶头黄色视频|