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

    現(xiàn)貨市場環(huán)境下新能源并網(wǎng)接入對市場出清的影響

    2022-01-05 11:02:46魏利屾馮宇昂方家琨艾小猛文勁宇
    上海交通大學(xué)學(xué)報 2021年12期
    關(guān)鍵詞:邊界條件現(xiàn)貨出力

    魏利屾, 馮宇昂, 方家琨, 艾小猛, 文勁宇

    (華中科技大學(xué) 強電磁與新技術(shù)國家重點實驗室, 武漢 430074)

    2015年國務(wù)院發(fā)布的《關(guān)于進一步深化電力體制改革的若干意見》[1]拉開了我國新一輪電力市場改革的序幕.其中,電力現(xiàn)貨市場建設(shè)是中國當(dāng)前電力體制改革的關(guān)鍵環(huán)節(jié)[2].電力現(xiàn)貨市場的建設(shè)會改變現(xiàn)有的計劃調(diào)度模式,計算得到不同時間、不同地點的電價分布情況,還原了電能的商品屬性,促進發(fā)電資源的優(yōu)化配置,同時形成了滿足電網(wǎng)安全約束、可實際物理執(zhí)行的發(fā)電調(diào)度指令.

    與此同時,中國的新能源如風(fēng)電、光伏等的裝機容量在過去十年快速提升,已經(jīng)成為全世界新能源裝機最多的國家[3].據(jù)國家統(tǒng)計局的數(shù)據(jù),2020年末全國發(fā)電裝機容量達到2.2×1012W,其中,并網(wǎng)風(fēng)電與光伏裝機容量分別為2.8×1011W、2.5×1011W[4].隨著“二氧化碳排放力爭于2030年前達到峰值,努力爭取2060年前實現(xiàn)碳中和”重大國家戰(zhàn)略的提出,電力行業(yè)作為碳排放的最主要來源,面臨更加迫切的向清潔低碳能源系統(tǒng)轉(zhuǎn)型的壓力[5].在可預(yù)見的未來,新能源的裝機容量還會持續(xù)提升[6],當(dāng)前以化石能源為主體的電力系統(tǒng)也將逐步轉(zhuǎn)變?yōu)橐孕履茉礊橹黧w的新型電力系統(tǒng)[7].

    為了促進新能源的發(fā)展,中國政府按照不同地區(qū)的風(fēng)光資源稟賦,規(guī)定了不同地區(qū)新能源的標桿上網(wǎng)電價[8],通過在資金、政策上的扶持,在初期新能源技術(shù)成本較高、滲透率較低時,極大地促進我國新能源產(chǎn)業(yè)的快速發(fā)展[9].然而,在新能源的快速發(fā)展過程中,也遇到了許多問題和挑戰(zhàn)[10].首先是隨著新能源滲透率提高帶來的棄風(fēng)、棄光問題[11],這極大地影響了新能源發(fā)電企業(yè)的收益情況,使投資者望而卻步,不利于新能源產(chǎn)業(yè)的長期發(fā)展[12].此外,隨著新能源裝機容量的不斷升高,相關(guān)的補貼資金赤字越來越大[9],并且缺少市場競爭環(huán)節(jié)的反饋,難以反映真實的供需關(guān)系與發(fā)電成本,長期的高收益也無法激勵新能源企業(yè)主動降低成本.因此,急需進一步完善市場化交易機制,同時分析與常規(guī)機組存在諸多差異的新能源機組對于市場出清結(jié)果的影響.

    本文主要對現(xiàn)貨市場環(huán)境下新能源對市場出清結(jié)果的影響進行定量分析評估,可為新能源的發(fā)展與市場機制的設(shè)計提供參考[13].關(guān)于新能源對于現(xiàn)貨市場的影響,現(xiàn)有文獻主要是基于國內(nèi)外的電力市場的實踐經(jīng)驗定性討論高比例新能源下現(xiàn)貨市場的規(guī)則設(shè)計[1, 14-15],而缺乏定量的分析評估結(jié)果,其分析是否適用于未來的新能源發(fā)展與市場運行也存在疑問.

    因此,本文基于安全約束機組組合和安全約束經(jīng)濟調(diào)度模型,對電力現(xiàn)貨市場運行規(guī)則與出清流程進行建模,搭建了電力現(xiàn)貨市場仿真分析框架,實現(xiàn)電力現(xiàn)貨市場的運行模擬,以得到電力系統(tǒng)與現(xiàn)貨市場的定量運行結(jié)果.以某實際省級電網(wǎng)為例,定量分析了現(xiàn)貨市場環(huán)境下新能源對市場出清價格、系統(tǒng)運行情況和機組收益的影響,提前發(fā)現(xiàn)潛在的問題與現(xiàn)象,為未來市場機制的設(shè)計與新能源發(fā)展規(guī)劃提供參考.

    1 現(xiàn)貨市場環(huán)境下新能源參與方式

    2017年8月國家發(fā)改委、能源局發(fā)布了《關(guān)于開展電力現(xiàn)貨市場建設(shè)試點工作的通知》[16],選擇了8個省級地區(qū)作為第一批電力現(xiàn)貨市場改革的試點地區(qū).以中國南方(以廣東為起步)電力現(xiàn)貨市場為例,電力現(xiàn)貨市場的組織流程主要包括邊界條件準備與公布、市場成員報價申報、市場出清計算、市場收益結(jié)算等部分.各試點地區(qū)根據(jù)自身的電力需求、網(wǎng)絡(luò)阻塞、經(jīng)濟水平等情況,在市場組織形式、機組申報方式、電價定價機制等規(guī)則設(shè)計方面存在一些差異,具體差異詳見文獻[17].

    新能源如何參與現(xiàn)貨市場是未來“構(gòu)建以新能源為主體的新型電力系統(tǒng)”背景下必須要面對的一個問題,表1給出了國內(nèi)外部分區(qū)域的新能源參與現(xiàn)貨市場方式[15].目前,部分試點地區(qū),例如廣東電力現(xiàn)貨市場,為落實“可再生能源保障性收購”的政策要求,將新能源作為現(xiàn)貨市場的邊界條件,即新能源出力不參與市場定價.與之相反,則是允許新能源電站報量報價,與常規(guī)的火電機組一道參與電力現(xiàn)貨市場,進行市場出清,如圖1所示.但新能源電站與常規(guī)火電機組在諸多方面存在不同.新能源電站的一次能源來源為風(fēng)光等自然資源,邊際發(fā)電成本極低,這與依靠化石能源發(fā)電的常規(guī)火電機組的成本特性存在極大差異.此外,在進行市場競爭的同時,還要保障新能源的消納水平,減少系統(tǒng)的碳排放.因此,新能源與常規(guī)機組進行同臺競價,會對電力市場形態(tài)造成較大改變,而如何對其影響進行量化分析,乃至針對存在的問題設(shè)計出適應(yīng)于高比例新能源下的電力現(xiàn)貨市場規(guī)則是必須要面對的問題.

    表1 國內(nèi)外新能源參與現(xiàn)貨市場方式Tab.1 Participation patterns of renewable energy in spot electricity market at home and abroad

    圖1 現(xiàn)貨市場環(huán)境下系統(tǒng)結(jié)構(gòu)圖Fig.1 System structure diagram under spot market environment

    2 電力現(xiàn)貨市場仿真分析框架

    2.1 概述

    本文基于實際某省級電網(wǎng)現(xiàn)貨市場規(guī)則,搭建了如圖2所示的現(xiàn)貨市場仿真分析框架,以實現(xiàn)現(xiàn)貨市場環(huán)境下新能源對市場出清結(jié)果的影響分析.

    圖2 電力現(xiàn)貨市場仿真框架Fig.2 Simulation framework of electricity spot market

    首先是形成市場出清的邊界條件,具體包括:負荷預(yù)測曲線、新能源出力曲線、聯(lián)絡(luò)線功率計劃以及電網(wǎng)拓撲結(jié)構(gòu)、機組運行特性參數(shù)和報價等,其中新能源的出力時序曲線可由歷史出力數(shù)據(jù)[18]或氣象數(shù)據(jù)[19]生成.然后,基于給出的邊界條件,由市場出清模型模擬市場出清計算,并形成最后的市場出清結(jié)果,即各機組的中標電量與市場電價情況、系統(tǒng)的新能源消納與碳排放情況等.其中,市場出清模型是電力現(xiàn)貨市場仿真的核心組件,由安全約束機組組合與安全約束經(jīng)濟調(diào)度組成.

    2.2 市場出清數(shù)學(xué)模型

    2.2.1安全約束機組組合 安全約束機組組合是在滿足相關(guān)安全運行約束的前提下,以調(diào)度范圍內(nèi)最小化購電費用為目標,計算每臺機組在每個時間步長內(nèi)(15 min或1 h)的開關(guān)機狀態(tài)與功率輸出.

    (1) 目標函數(shù).

    (1)

    (2) 功率平衡約束.

    (2)

    式中:PG,i(t)、PR,j(t)分別為火電機組與新能源機組的中標電量,等于對應(yīng)的所有分段的中標電量之和;T(t)為省間聯(lián)絡(luò)的凈輸入功率;D(t)為該省在時刻t的負荷功率.該約束保證系統(tǒng)的電力供給與需求相平衡.

    (3) 系統(tǒng)備用約束.

    (3)

    Rdn(t)

    (4)

    (4) 機組出力上下限約束.

    (5)

    (6)

    對火電機組,其出力受到最大技術(shù)出力和最小技術(shù)出力的限制;而對于新能源機組,其出力受到其理論最大出力的限制(與風(fēng)速、光照條件有關(guān)).對于新能源,若市場規(guī)則規(guī)定其出力作為邊界條件參與電力現(xiàn)貨市場,則新能源出力應(yīng)固定為其理論最大出力,即式(6)應(yīng)修改為

    (7)

    (5) 機組爬坡能力約束.

    PG,i(t)-PG,i(t-1)≤

    RU,i+M[1-ui(t-1)]+M[1-ui(t)]

    (8)

    PG,i(t-1)-PG,i(t)≤

    RD,i+M[1-ui(t-1)]+M[1-ui(t)]

    (9)

    (10)

    (11)

    式中:RU,i、RD,i分別為常規(guī)機組i的上爬坡能力與下爬坡能力(一個時間步長內(nèi));M為一個很大的權(quán)重系數(shù).該約束反映了火電機組在連續(xù)時段的出力變化還受到其爬坡能力的限制.

    (6) 最小啟停機時間約束.

    (12)

    (13)

    TU,i(t)=min {MU,i,T-t+1}

    (14)

    TD,i(t)=min {MD,i,T-t+1}

    (15)

    式中:TU,i(t)為常規(guī)機組i在時刻t仍然需要繼續(xù)保持運行的時間;τ為時段索引號;TD,i(t)為常規(guī)機組i在時刻t仍然需要繼續(xù)保持停機的時間;MU,i為常規(guī)機組i的最小開機時間;MD,i為常規(guī)機組i的最小停機時間.因為常規(guī)機組不能頻繁啟停,一則是技術(shù)條件不允許,二則頻繁啟停會導(dǎo)致機組的壽命減少,因此設(shè)置了最小開關(guān)機時間約束,保證機組的啟動與停機都需要一定的時間.

    (7) 潮流約束.

    (16)

    式中:Fl,max為線路(即斷面)l的最大傳輸容量;b為節(jié)點序號;l為線路序號;Gl,b為節(jié)點b對于線路l的功率傳輸分配因子;Pinj,b(t)為節(jié)點b的注入功率,由節(jié)點上的發(fā)電機出力與負荷功率值相減計算而得.實際上,線路傳輸功率應(yīng)該由交流潮流計算得到,但由于交流潮流模型是非線性模型,大大提高了問題的計算復(fù)雜度[22],所以國內(nèi)外電力市場都采用了線性模型的直流潮流法進行建模與計算.

    2.2.2安全約束經(jīng)濟調(diào)度 安全約束經(jīng)濟調(diào)度是指在發(fā)電機啟停狀態(tài)已知的情況下,在滿足系統(tǒng)運行的相關(guān)安全約束的前提下,對發(fā)電機組的出力進行優(yōu)化從而實現(xiàn)系統(tǒng)的購電費用最小化.安全約束經(jīng)濟調(diào)度的目標函數(shù)為

    (17)

    由于機組的啟停計劃是由安全約束機組組合所確定,所以目標函數(shù)(17)相較于目標函數(shù)(1)減少了機組的啟停費用項,其余項不變.對于約束條件,由于機組的啟停計劃已經(jīng)確定,經(jīng)濟調(diào)度模型中不再考慮啟停爬坡約束與最小開關(guān)機時間約束(式(10)~(15)),其余約束不變.

    通過安全約束機組組合后,得到了機組在調(diào)度時段內(nèi)的啟停狀態(tài),然后基于得到的機組啟停計劃,經(jīng)過安全約束經(jīng)濟調(diào)度,就能得到機組的功率輸出情況、全網(wǎng)的新能源消納情況與全系統(tǒng)的發(fā)電費用、碳排放量等信息.

    2.2.3電價計算 根據(jù)邊際價格理論,在t時刻,電價m(t)由下式進行計算:

    (18)

    在t時刻,增加單位負荷,系統(tǒng)的總購電費用的增量,為系統(tǒng)的邊際價格.在網(wǎng)絡(luò)沒有阻塞時,出清電價即為式(2)的拉格朗日乘子;當(dāng)網(wǎng)絡(luò)發(fā)生阻塞時則是式(2)、(16)的拉格朗日乘子的線性組合[23].所建立的安全約束機組組合、安全約束經(jīng)濟調(diào)度以及電價計算均是在Anaconda環(huán)境下使用Python 3.7調(diào)用CPLEX求解器進行求解.

    3 案例介紹

    本文選擇某實際省級電網(wǎng)的電網(wǎng)運行數(shù)據(jù)為研究對象,基于本文搭建的現(xiàn)貨市場仿真分析框架,定量分析現(xiàn)貨市場環(huán)境下新能源對市場出清結(jié)果的影響.當(dāng)前,該省級電網(wǎng)具有燃煤、燃氣機組共200余臺,新能源滲透率達19.52%.預(yù)計到2025年,該省電網(wǎng)將規(guī)劃建設(shè)較多的風(fēng)電機組,此時新能源滲透率將達到25.68%.

    本文分別從春、夏、秋、冬四季選擇典型周進行相關(guān)分析,代替全年的市場運行結(jié)果[24],其負荷曲線形狀如圖3所示,圖中D*為電力負荷標幺值.目前,電力現(xiàn)貨市場的市場化機組主要是由燃煤、燃氣機組組成,核電機組、抽水蓄能機組、省間聯(lián)絡(luò)線等均是按照各自的發(fā)電/輸電計劃作為市場出清的邊界條件.對于新能源電站,考慮兩種參與方式:① 新能源電站出力作為邊界條件參與市場出清;② 新能源電站報量報價參與現(xiàn)貨市場.

    圖3 四季典型周電力負荷曲線Fig.3 Representative weekly demand profiles of different seasons

    燃煤、燃氣機組的分段報價通過機組的成本特性曲線進行生成,而燃煤、燃氣、新能源機組的年化投資成本則使用凈現(xiàn)值法進行計算[19]:

    (19)

    式中:Q為工程單位容量固定投資費用;x為貼現(xiàn)率,一般為7%;y為機組使用年限.

    對于新能源機組,由于其邊際發(fā)電成本極低,趨近于0,而在邊際定價機制下,市場成員會傾向于按照邊際成本進行報價[14],所以假設(shè)新能源機組在現(xiàn)貨市場中報價為0[25].

    4 案例分析

    4.1 對市場出清電價與系統(tǒng)運行情況的影響

    為分析新能源對現(xiàn)貨市場出清價格與系統(tǒng)運行情況的影響,本文對不同新能源滲透率下,考慮上述兩種新能源參與方式,得到全年的市場出清結(jié)果,分別如表2所示.

    表2 不同水平年下的市場出清結(jié)果對比Tab.2 Comparison of the market-clearing results in different years

    從系統(tǒng)運行情況來看,首先,新能源消納率上,當(dāng)新能源作為市場化機組報量報價參與現(xiàn)貨市場后,相較于直接作為邊界條件,新能源消納率有一定幅度的下降,且隨著新能源滲透率的提高,該現(xiàn)象更加明顯(下降幅度從0.08%增加到1.23%).這是因為當(dāng)新能源作為邊界條件進行出清時,是將新能源電站作為優(yōu)先發(fā)電機組,然后在進行其余市場化機組的購電費用最小化,這樣是最大化利用了電力系統(tǒng)的新能源消納能力.然而當(dāng)新能源以報量報價的形式參與市場后,新能源電站會與常規(guī)機組一道在市場里進行同臺競爭,雖然新能源機組由于在市場競價中具有較大優(yōu)勢,但由于新能源的波動性、間歇性等負外部性,優(yōu)先消納新能源反而會導(dǎo)致系統(tǒng)總的運行成本提高.即在現(xiàn)貨市場環(huán)境下,市場決策會通過一定數(shù)量的棄風(fēng)棄光來降低系統(tǒng)總的運行成本.

    從市場出清電價上看,新能源參與現(xiàn)貨市場會較大幅度地降低系統(tǒng)的出清電價,同樣隨著新能源滲透率的提高,該現(xiàn)象會更加明顯(下降幅度從3.44 元/(MW·h)增加到18.82 元/(MW·h)).以2025年冬季典型周為例解釋該現(xiàn)象出現(xiàn)的原因,其電力電量平衡圖如圖4所示,圖中P*為功率標幺值,市場平均出清價格如圖5所示,圖中λ為出清電價.可以看到,當(dāng)允許新能源機組報價后,在部分時刻會出現(xiàn)零電價的情況.這是因為市場的出清電價是由邊際機組決定,而新能源參與市場后,如果出現(xiàn)棄風(fēng)棄光,則新能源機組將成為邊際機組(即新增 1 MW 負荷將優(yōu)先調(diào)用被削減的新能源發(fā)電資源),市場出清價格將由新能源機組的報價所決定,即出現(xiàn)了零電價.而如果將新能源發(fā)電作為邊界條件,系統(tǒng)則會盡可能消納新能源,即使由于系統(tǒng)的物理消納極限受限制而導(dǎo)致出現(xiàn)了少量的棄風(fēng)棄光現(xiàn)象,也因為其不能成為邊際機組而不會出現(xiàn)零電價的現(xiàn)象.并且,隨著新能源滲透率的提高,將新能源作為邊界條件出清對于靈活性的需求會越來越高,這導(dǎo)致燃料成本更高的氣電機組在更多時刻成為邊際機組.因此,當(dāng)新能源滲透率提高時,雖然平均發(fā)電成本下降了8.25 元/(MW·h),但出清電價反而增加了8.37 元/(MW·h).

    圖4 冬季典型周電力電量平衡結(jié)果Fig.4 Power balancing results of a representative week in the winter

    圖5 冬季典型周出清電價Fig.5 Clearing prices of a representative week in the winter

    4.2 對市場化機組的收益影響

    在現(xiàn)貨市場環(huán)境下,新能源不但會對市場出清電價造成影響,還會擠占常規(guī)機組的發(fā)電空間,這些因素勢必會對機組的盈利情況造成巨大影響.在新能源電力系統(tǒng)中,機組參與現(xiàn)貨市場是否有充足的收入覆蓋機組的燃料成本與投資成本,并存在合理的利潤空間,對于現(xiàn)貨市場乃至電力行業(yè)的健康發(fā)展具有重要意義.

    圖6所示為2個水平年下,燃煤機組與燃氣機組的年度收益情況,圖中φ為利潤.對于煤電機組,隨著新能源裝機比例的不斷提高,在新能源參與現(xiàn)貨市場后,其收益均能覆蓋其燃料成本與投資成本,并保證一定收益,這是因為煤電機組的燃料費用較低,在現(xiàn)貨市場中發(fā)電順序較高(僅次于新能源機組).

    圖6 常規(guī)機組收益情況Fig.6 Revenue of conventional units

    對于氣電機組,由于新能源的裝機比例提高,為應(yīng)對新能源發(fā)電的波動性與間歇性,更為靈活、成本較高的氣電機組越來越多地承擔(dān)調(diào)峰任務(wù),導(dǎo)致其出力被擠壓在最小技術(shù)出力處或停機,更多時刻成為非邊際機組,其報價沒有參與市場定價,出清電價反而由成本更低的煤電機組與新能源機組報價決定,進而導(dǎo)致氣電機組的虧損.

    此外,在新能源作為邊界條件出清的前提下,隨著新能源裝機的增加,出現(xiàn)了機組利潤略微增加的現(xiàn)象,這是由于新能源滲透率增加后,由于靈活性緊張,氣電機組成為邊際機組的時刻增加.而與之相反地,在新能源報量報價情況下,由于其對出清電價的下降作用與擠占傳統(tǒng)電源的發(fā)電空間等原因,常規(guī)機組利潤有不同程度的下降.

    對于新能源機組,當(dāng)作為邊界條件出清時,按照政府規(guī)定的上網(wǎng)電價進行支付,其中已經(jīng)考慮了其運行成本與投資成本.但當(dāng)新能源報量報價參與市場時,新能源的收益就由市場定價所決定.新能源機組的收益情況如圖7所示,雖然新能源機組的邊際發(fā)電成本很低,不考慮投資成本時機組收益很高(發(fā)電順序高于燃煤機組),但由于目前新能源機組,特別是海上風(fēng)電機組的技術(shù)成熟度相對常規(guī)機組較低,而投資成本較高,這將導(dǎo)致新能源機組參與現(xiàn)貨市場可能存在虧損的現(xiàn)象.對于陸上風(fēng)電,僅需要0.03元/(kW·h)左右的補貼,就能在現(xiàn)貨市場中實現(xiàn)收支相抵,且考慮到風(fēng)電成本的持續(xù)下降[26],補貼力度可以進一步退坡.但對于海上風(fēng)電,還至少需要0.4~0.5元/(kW·h)的補貼方能在現(xiàn)貨市場實現(xiàn)收支平衡.

    圖7 新能源機組收益情況Fig.7 Revenue of renewable energy units

    5 結(jié)論

    在電力市場化改革的背景下,要構(gòu)建以新能源為主體的新型電力系統(tǒng),盡快實現(xiàn)“碳達峰、碳中和”目標,就需要定量分析現(xiàn)貨市場環(huán)境下新能源對于市場出清結(jié)果的影響.在此背景下,本文基于搭建的現(xiàn)貨市場仿真分析框架,得到如下結(jié)論,以期對未來市場機制設(shè)計與新能源發(fā)展規(guī)劃起到指導(dǎo)作用.

    (1) 相較于直接作為邊界條件,新能源報量報價參與現(xiàn)貨市場后,市場會通過一定量的棄風(fēng)棄光,降低系統(tǒng)總的發(fā)電成本,但新能源報量報價參與市場的情況下,不能最大化發(fā)揮系統(tǒng)的新能源消納作用.

    (2) 新能源報量報價參與現(xiàn)貨市場后,電價波動性更加劇烈,平均電價降低甚至?xí)诓糠謺r刻現(xiàn)貨市場出現(xiàn)零電價的現(xiàn)象.而作為邊界條件出清時,隨著新能源滲透率的提高,會有更多時刻氣電機組成為邊際機組,反而可能導(dǎo)致出清電價略有增加.

    (3) 在現(xiàn)貨市場環(huán)境下,新能源會造成部分市場化機組成本難以回收的問題,尤其是氣電機組,囿于高額的燃料成本,短期運行成本都無法完全覆蓋;而對于新能源機組,特別是海上風(fēng)電機組,仍需要依賴于高額的補貼方能實現(xiàn)在現(xiàn)貨市場的收支相抵.

    猜你喜歡
    邊界條件現(xiàn)貨出力
    美國大豆現(xiàn)貨價格上漲
    豆粕:美豆新高后回落振蕩,現(xiàn)貨臨近春節(jié)成交轉(zhuǎn)淡
    一類帶有Stieltjes積分邊界條件的分數(shù)階微分方程邊值問題正解
    帶有積分邊界條件的奇異攝動邊值問題的漸近解
    浙江電力現(xiàn)貨啟幕
    能源(2017年12期)2018-01-31 01:43:13
    備戰(zhàn)電力現(xiàn)貨
    能源(2017年12期)2018-01-31 01:43:11
    風(fēng)電場有功出力的EEMD特性分析
    要爭做出力出彩的黨員干部
    河南電力(2016年5期)2016-02-06 02:11:35
    風(fēng)電場群出力的匯聚效應(yīng)分析
    電測與儀表(2015年5期)2015-04-09 11:31:12
    帶Robin邊界條件的2維隨機Ginzburg-Landau方程的吸引子
    12—13女人毛片做爰片一| 国产高清videossex| 法律面前人人平等表现在哪些方面| 色尼玛亚洲综合影院| 狂野欧美激情性xxxx| 亚洲精品国产精品久久久不卡| 伊人久久大香线蕉亚洲五| 50天的宝宝边吃奶边哭怎么回事| 免费观看a级毛片全部| 亚洲av第一区精品v没综合| 免费在线观看完整版高清| 亚洲精品在线观看二区| 蜜桃国产av成人99| 亚洲欧美精品综合一区二区三区| 丁香六月欧美| av福利片在线| 国产精品久久久久久人妻精品电影 | 91精品国产国语对白视频| 99精品久久久久人妻精品| 欧美精品一区二区免费开放| 成在线人永久免费视频| 天堂动漫精品| 中文亚洲av片在线观看爽 | 精品国产乱码久久久久久男人| avwww免费| 97在线人人人人妻| 日本黄色视频三级网站网址 | svipshipincom国产片| 日韩一卡2卡3卡4卡2021年| 精品欧美一区二区三区在线| 黄色丝袜av网址大全| 在线观看免费高清a一片| 热re99久久国产66热| 亚洲自偷自拍图片 自拍| 99国产精品一区二区三区| 一级a爱视频在线免费观看| 国产一区二区三区视频了| 亚洲成人免费av在线播放| 18禁国产床啪视频网站| 成年人免费黄色播放视频| 日本av免费视频播放| 久久久久久久久免费视频了| 12—13女人毛片做爰片一| 十八禁网站免费在线| 男人操女人黄网站| 啦啦啦中文免费视频观看日本| 动漫黄色视频在线观看| 国产在线免费精品| 99精品久久久久人妻精品| 99精品在免费线老司机午夜| 亚洲伊人色综图| 女人久久www免费人成看片| 一进一出抽搐动态| 国产精品熟女久久久久浪| 免费黄频网站在线观看国产| 国产成人精品久久二区二区91| 狂野欧美激情性xxxx| 国产精品久久久久成人av| videos熟女内射| 黑丝袜美女国产一区| 女人精品久久久久毛片| 女人爽到高潮嗷嗷叫在线视频| 亚洲中文字幕日韩| 啦啦啦 在线观看视频| 男人操女人黄网站| 天天躁日日躁夜夜躁夜夜| 亚洲精品在线观看二区| cao死你这个sao货| 在线播放国产精品三级| 国产亚洲av高清不卡| 国产不卡一卡二| 天天添夜夜摸| 国产亚洲欧美精品永久| 麻豆乱淫一区二区| 亚洲成人国产一区在线观看| av福利片在线| 两性午夜刺激爽爽歪歪视频在线观看 | 91精品国产国语对白视频| 日韩制服丝袜自拍偷拍| 精品久久久久久久毛片微露脸| 老司机午夜十八禁免费视频| 两性午夜刺激爽爽歪歪视频在线观看 | a级毛片黄视频| 国产淫语在线视频| 欧美日韩视频精品一区| 午夜91福利影院| 国产男靠女视频免费网站| 少妇裸体淫交视频免费看高清 | 国产福利在线免费观看视频| 国产精品一区二区在线不卡| 少妇粗大呻吟视频| 精品久久久久久久毛片微露脸| 99热网站在线观看| 又紧又爽又黄一区二区| 日本黄色日本黄色录像| 夜夜夜夜夜久久久久| 久久精品亚洲熟妇少妇任你| 国产91精品成人一区二区三区 | 乱人伦中国视频| 99精国产麻豆久久婷婷| 伦理电影免费视频| 女人高潮潮喷娇喘18禁视频| 久久国产精品大桥未久av| 亚洲国产av影院在线观看| 亚洲欧美一区二区三区久久| 欧美黑人精品巨大| 后天国语完整版免费观看| 国产1区2区3区精品| 嫩草影视91久久| 国产精品自产拍在线观看55亚洲 | 亚洲一区中文字幕在线| 国产av一区二区精品久久| 丝袜人妻中文字幕| 国产99久久九九免费精品| 男人舔女人的私密视频| 大香蕉久久网| 成年女人毛片免费观看观看9 | 亚洲国产欧美日韩在线播放| 天天添夜夜摸| 国产精品一区二区在线不卡| 国产欧美日韩一区二区三区在线| 亚洲欧美一区二区三区黑人| 丰满迷人的少妇在线观看| 欧美日韩亚洲高清精品| 91av网站免费观看| 国产日韩一区二区三区精品不卡| 国产日韩欧美亚洲二区| 纵有疾风起免费观看全集完整版| 窝窝影院91人妻| 亚洲伊人色综图| 一夜夜www| 国产精品久久电影中文字幕 | 国产精品影院久久| 又大又爽又粗| 亚洲欧美一区二区三区久久| 黄片小视频在线播放| 最黄视频免费看| 亚洲午夜理论影院| 他把我摸到了高潮在线观看 | 色在线成人网| 亚洲一区二区三区欧美精品| 一区二区三区国产精品乱码| 国产不卡一卡二| 黄片大片在线免费观看| 国产成人免费无遮挡视频| 亚洲美女黄片视频| 成年人黄色毛片网站| 啦啦啦视频在线资源免费观看| 97在线人人人人妻| av网站免费在线观看视频| 国产日韩一区二区三区精品不卡| 日韩成人在线观看一区二区三区| 如日韩欧美国产精品一区二区三区| 99国产精品99久久久久| 99re在线观看精品视频| 中文欧美无线码| 王馨瑶露胸无遮挡在线观看| 另类亚洲欧美激情| 亚洲人成电影免费在线| av不卡在线播放| a级片在线免费高清观看视频| 99精品久久久久人妻精品| 一本大道久久a久久精品| av在线播放免费不卡| 亚洲五月婷婷丁香| 两性午夜刺激爽爽歪歪视频在线观看 | 黑人巨大精品欧美一区二区mp4| 99国产精品免费福利视频| 久久这里只有精品19| 色婷婷久久久亚洲欧美| 国产一区二区三区综合在线观看| 国产成人一区二区三区免费视频网站| 国产成人av激情在线播放| 亚洲三区欧美一区| 韩国精品一区二区三区| 欧美精品高潮呻吟av久久| 成年人免费黄色播放视频| 精品福利永久在线观看| 两个人看的免费小视频| 啦啦啦免费观看视频1| 亚洲av国产av综合av卡| 精品国产国语对白av| 亚洲第一av免费看| 成人黄色视频免费在线看| 精品国内亚洲2022精品成人 | 欧美黄色淫秽网站| 成人18禁高潮啪啪吃奶动态图| 精品少妇内射三级| 国产高清激情床上av| 国产av又大| 精品久久久精品久久久| 另类亚洲欧美激情| 老司机午夜福利在线观看视频 | 人成视频在线观看免费观看| 精品国产国语对白av| 国产精品久久久人人做人人爽| 亚洲一卡2卡3卡4卡5卡精品中文| 曰老女人黄片| 亚洲人成77777在线视频| 午夜久久久在线观看| 免费在线观看影片大全网站| 午夜福利在线观看吧| av国产精品久久久久影院| 午夜福利视频在线观看免费| 操出白浆在线播放| 欧美日韩av久久| 国产欧美亚洲国产| 午夜久久久在线观看| 在线亚洲精品国产二区图片欧美| 麻豆av在线久日| 亚洲国产欧美网| 久久久精品免费免费高清| 日韩有码中文字幕| 国产精品亚洲一级av第二区| 欧美日韩亚洲综合一区二区三区_| 久久久精品94久久精品| 免费看a级黄色片| a级毛片在线看网站| 脱女人内裤的视频| 另类亚洲欧美激情| 亚洲第一欧美日韩一区二区三区 | 天天操日日干夜夜撸| 一级毛片女人18水好多| 在线观看一区二区三区激情| 午夜精品久久久久久毛片777| 无限看片的www在线观看| 日韩欧美三级三区| 欧美一级毛片孕妇| 亚洲 国产 在线| cao死你这个sao货| 香蕉久久夜色| 亚洲久久久国产精品| 国产单亲对白刺激| 91字幕亚洲| 天天操日日干夜夜撸| 国产成人免费无遮挡视频| 午夜福利乱码中文字幕| 欧美黄色淫秽网站| 大型黄色视频在线免费观看| 国产成人影院久久av| 香蕉久久夜色| 19禁男女啪啪无遮挡网站| 国产精品 国内视频| 首页视频小说图片口味搜索| 999久久久国产精品视频| 亚洲黑人精品在线| 欧美日韩国产mv在线观看视频| 国产精品国产av在线观看| 久久九九热精品免费| 超碰成人久久| 一个人免费在线观看的高清视频| 丰满迷人的少妇在线观看| 波多野结衣av一区二区av| 69精品国产乱码久久久| 国产视频一区二区在线看| 肉色欧美久久久久久久蜜桃| 美女福利国产在线| 欧美午夜高清在线| 十分钟在线观看高清视频www| 亚洲欧美日韩另类电影网站| 亚洲欧美一区二区三区黑人| 国产精品电影一区二区三区 | 母亲3免费完整高清在线观看| 99精品在免费线老司机午夜| 中文字幕人妻丝袜一区二区| 美女扒开内裤让男人捅视频| 欧美av亚洲av综合av国产av| e午夜精品久久久久久久| 亚洲中文字幕日韩| 亚洲av第一区精品v没综合| 国产成人免费无遮挡视频| 亚洲精品一二三| 一区二区av电影网| 久久久久久人人人人人| 亚洲天堂av无毛| 亚洲一区中文字幕在线| 日本a在线网址| 丁香六月欧美| www.999成人在线观看| svipshipincom国产片| 亚洲成人免费av在线播放| 国产av又大| 久久国产精品影院| 亚洲av欧美aⅴ国产| netflix在线观看网站| 一本一本久久a久久精品综合妖精| 国产精品亚洲av一区麻豆| 每晚都被弄得嗷嗷叫到高潮| 十八禁人妻一区二区| 两性午夜刺激爽爽歪歪视频在线观看 | 99热网站在线观看| 妹子高潮喷水视频| 99国产综合亚洲精品| 久久久久视频综合| 超碰97精品在线观看| 99国产精品免费福利视频| 两人在一起打扑克的视频| 999精品在线视频| 亚洲av国产av综合av卡| 久久国产精品影院| 久久精品成人免费网站| 黄色视频在线播放观看不卡| 久久99一区二区三区| 国产成人精品久久二区二区免费| videosex国产| 啦啦啦在线免费观看视频4| 久久久精品免费免费高清| 亚洲伊人久久精品综合| 中文字幕制服av| 最近最新免费中文字幕在线| 曰老女人黄片| 久久人妻熟女aⅴ| 国产1区2区3区精品| 黄片小视频在线播放| 777米奇影视久久| 亚洲美女黄片视频| 18禁裸乳无遮挡动漫免费视频| 日韩欧美国产一区二区入口| 精品一区二区三区四区五区乱码| 亚洲精品在线观看二区| 亚洲伊人久久精品综合| 日本五十路高清| 成人黄色视频免费在线看| 亚洲欧美日韩另类电影网站| 最近最新免费中文字幕在线| 老司机在亚洲福利影院| 欧美在线一区亚洲| 久久久久网色| 一夜夜www| 成人永久免费在线观看视频 | 免费在线观看视频国产中文字幕亚洲| 俄罗斯特黄特色一大片| 夜夜骑夜夜射夜夜干| videosex国产| 黄色毛片三级朝国网站| 亚洲精品国产精品久久久不卡| 我的亚洲天堂| 99国产精品免费福利视频| 国产一区二区 视频在线| 国产欧美日韩一区二区三| 国产亚洲欧美精品永久| 国产高清激情床上av| 人妻久久中文字幕网| 91成人精品电影| 亚洲人成伊人成综合网2020| 99久久人妻综合| 欧美日韩一级在线毛片| 纯流量卡能插随身wifi吗| 欧美精品一区二区免费开放| 每晚都被弄得嗷嗷叫到高潮| 国产极品粉嫩免费观看在线| av线在线观看网站| av片东京热男人的天堂| 91字幕亚洲| 女警被强在线播放| 成人国产一区最新在线观看| 三级毛片av免费| 成人国产一区最新在线观看| 免费不卡黄色视频| 国产亚洲精品第一综合不卡| 成人国产av品久久久| 国产极品粉嫩免费观看在线| 国产深夜福利视频在线观看| 菩萨蛮人人尽说江南好唐韦庄| 国产深夜福利视频在线观看| 夜夜夜夜夜久久久久| 精品亚洲成a人片在线观看| 精品国产一区二区久久| 男女午夜视频在线观看| 国产成人欧美在线观看 | 免费观看人在逋| 丰满少妇做爰视频| 性色av乱码一区二区三区2| 女人爽到高潮嗷嗷叫在线视频| 国产av精品麻豆| 黑人操中国人逼视频| 国产精品自产拍在线观看55亚洲 | 国产一区有黄有色的免费视频| 国产单亲对白刺激| 久久国产精品男人的天堂亚洲| 中文字幕色久视频| 男女无遮挡免费网站观看| www.精华液| 考比视频在线观看| 人人澡人人妻人| 人人妻人人爽人人添夜夜欢视频| 精品国内亚洲2022精品成人 | 男女下面插进去视频免费观看| 一边摸一边抽搐一进一出视频| 亚洲国产精品一区二区三区在线| 亚洲欧美精品综合一区二区三区| 久久精品亚洲av国产电影网| 在线亚洲精品国产二区图片欧美| 色婷婷av一区二区三区视频| 国产一区有黄有色的免费视频| 婷婷丁香在线五月| 日韩精品免费视频一区二区三区| 黑人猛操日本美女一级片| 国产高清激情床上av| 国产片内射在线| 久久久欧美国产精品| 国产高清videossex| 亚洲,欧美精品.| 精品国产乱子伦一区二区三区| 99国产精品免费福利视频| 国产av国产精品国产| 国产欧美亚洲国产| 美国免费a级毛片| 肉色欧美久久久久久久蜜桃| 欧美精品亚洲一区二区| a在线观看视频网站| 亚洲精品在线美女| 国产成人av教育| 亚洲av成人不卡在线观看播放网| 久久久久视频综合| av视频免费观看在线观看| 国产日韩一区二区三区精品不卡| 国产精品九九99| 黑人操中国人逼视频| 久久人妻熟女aⅴ| 亚洲色图av天堂| 纵有疾风起免费观看全集完整版| 女同久久另类99精品国产91| 欧美乱妇无乱码| 老汉色∧v一级毛片| 18禁黄网站禁片午夜丰满| 老熟女久久久| 一级片'在线观看视频| h视频一区二区三区| 一级毛片电影观看| 欧美黑人欧美精品刺激| 国产成人精品久久二区二区免费| 国产成人精品在线电影| 国产在视频线精品| 日韩熟女老妇一区二区性免费视频| 欧美大码av| 国产99久久九九免费精品| 在线观看一区二区三区激情| 久久精品亚洲av国产电影网| 精品国产国语对白av| 久久精品国产亚洲av香蕉五月 | 欧美日韩视频精品一区| 男人舔女人的私密视频| 日韩大片免费观看网站| av欧美777| 两个人看的免费小视频| 亚洲av电影在线进入| 国产在视频线精品| 免费看十八禁软件| 狠狠精品人妻久久久久久综合| 亚洲免费av在线视频| 在线观看人妻少妇| 夫妻午夜视频| 99re6热这里在线精品视频| 窝窝影院91人妻| 丝袜美腿诱惑在线| 黄色视频在线播放观看不卡| 成人黄色视频免费在线看| 丝瓜视频免费看黄片| 欧美一级毛片孕妇| 亚洲男人天堂网一区| 日日爽夜夜爽网站| 中文字幕另类日韩欧美亚洲嫩草| 水蜜桃什么品种好| 男人操女人黄网站| 在线观看人妻少妇| 搡老乐熟女国产| 久久久水蜜桃国产精品网| 久久精品国产综合久久久| 亚洲欧美一区二区三区黑人| 日本黄色视频三级网站网址 | 久久精品国产综合久久久| 最近最新中文字幕大全免费视频| 在线观看www视频免费| 国产欧美日韩综合在线一区二区| 啦啦啦免费观看视频1| 最近最新中文字幕大全电影3 | 视频在线观看一区二区三区| 男人舔女人的私密视频| 欧美另类亚洲清纯唯美| 欧美日本中文国产一区发布| 一本大道久久a久久精品| 中国美女看黄片| 三上悠亚av全集在线观看| 十分钟在线观看高清视频www| 久久精品国产亚洲av香蕉五月 | 日韩制服丝袜自拍偷拍| 大片电影免费在线观看免费| 91精品三级在线观看| 91九色精品人成在线观看| 一级片免费观看大全| 欧美日韩亚洲高清精品| 狂野欧美激情性xxxx| 欧美老熟妇乱子伦牲交| 久久精品国产综合久久久| 电影成人av| 纯流量卡能插随身wifi吗| 精品欧美一区二区三区在线| 男女床上黄色一级片免费看| 亚洲专区字幕在线| 久久精品亚洲av国产电影网| 亚洲美女黄片视频| 丰满迷人的少妇在线观看| 亚洲av电影在线进入| 欧美精品亚洲一区二区| 国产国语露脸激情在线看| 国产淫语在线视频| 丁香六月欧美| 欧美精品人与动牲交sv欧美| 久久久久久久国产电影| 国产精品免费大片| 亚洲熟妇熟女久久| www日本在线高清视频| 欧美日韩成人在线一区二区| 久久精品国产a三级三级三级| 国产精品欧美亚洲77777| 大香蕉久久网| 精品福利观看| av福利片在线| av欧美777| 黄色毛片三级朝国网站| 夜夜爽天天搞| 国产欧美日韩一区二区三区在线| 美女扒开内裤让男人捅视频| 国产三级黄色录像| 水蜜桃什么品种好| 久久这里只有精品19| 国产免费av片在线观看野外av| 成人精品一区二区免费| 久久精品91无色码中文字幕| 亚洲精品在线观看二区| 国产精品一区二区在线不卡| 久久中文字幕一级| 成人特级黄色片久久久久久久 | 亚洲全国av大片| av网站在线播放免费| www.999成人在线观看| 热99国产精品久久久久久7| 男人操女人黄网站| 精品国产一区二区三区四区第35| 免费在线观看日本一区| 免费观看av网站的网址| 麻豆成人av在线观看| 性少妇av在线| 宅男免费午夜| 一区二区三区乱码不卡18| 免费人妻精品一区二区三区视频| 一个人免费看片子| 久久久久久久久免费视频了| 亚洲精品中文字幕一二三四区 | 成人国语在线视频| 亚洲男人天堂网一区| 老司机深夜福利视频在线观看| 欧美日韩视频精品一区| 久久狼人影院| 热re99久久国产66热| 又大又爽又粗| tube8黄色片| 免费在线观看日本一区| 夜夜爽天天搞| 97在线人人人人妻| 天堂8中文在线网| 国产福利在线免费观看视频| 纯流量卡能插随身wifi吗| 国产精品久久电影中文字幕 | 最近最新免费中文字幕在线| 亚洲伊人色综图| 精品第一国产精品| 狠狠狠狠99中文字幕| 日韩免费av在线播放| 午夜成年电影在线免费观看| 国产精品自产拍在线观看55亚洲 | 91九色精品人成在线观看| 午夜视频精品福利| 黄片大片在线免费观看| 国产成人欧美| 国产一区二区在线观看av| 在线观看人妻少妇| 成人三级做爰电影| av欧美777| 久久久久久久久免费视频了| 大香蕉久久网| 超碰成人久久| 视频在线观看一区二区三区| 老汉色∧v一级毛片| 精品国产一区二区久久| 国产精品麻豆人妻色哟哟久久| svipshipincom国产片| 亚洲久久久国产精品| 色综合婷婷激情| 国产成人精品无人区| 一二三四在线观看免费中文在| 亚洲色图 男人天堂 中文字幕| 黑人巨大精品欧美一区二区蜜桃| 精品久久蜜臀av无| 国产免费av片在线观看野外av| 久久99热这里只频精品6学生| 一个人免费在线观看的高清视频| 免费在线观看视频国产中文字幕亚洲| 国产精品 国内视频| 久久婷婷成人综合色麻豆| 亚洲成a人片在线一区二区| 男女无遮挡免费网站观看| 久久久久久久久免费视频了| 成人亚洲精品一区在线观看| 亚洲va日本ⅴa欧美va伊人久久| 亚洲av国产av综合av卡| 搡老岳熟女国产| 国产免费视频播放在线视频| 亚洲av国产av综合av卡| 日本欧美视频一区| 亚洲午夜理论影院| 久久久国产成人免费| 80岁老熟妇乱子伦牲交| 天天躁日日躁夜夜躁夜夜|