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

    混雜免疫多目標(biāo)優(yōu)化算法及對動態(tài)經(jīng)濟(jì)環(huán)境調(diào)度問題優(yōu)化

    2023-10-18 05:06:59唐湘黔錢淑渠武慧虹
    關(guān)鍵詞:多目標(biāo)優(yōu)化免疫系統(tǒng)

    唐湘黔 錢淑渠 武慧虹

    摘 要:動態(tài)經(jīng)濟(jì)環(huán)境調(diào)度(DEED)問題是電力系統(tǒng)調(diào)度中一類含大規(guī)模約束的高維多目標(biāo)優(yōu)化問題,傳統(tǒng)的進(jìn)化算法易于陷入局部最優(yōu),使得所獲的Pareto前沿分布性和收斂性差。為了充分挖掘免疫系統(tǒng)的克隆選擇原理,提出一種混雜免疫多目標(biāo)優(yōu)化算法(HIMOA)。該算法以傳統(tǒng)進(jìn)化算法為基本框架,面對高維決策變量優(yōu)化易于陷入局部最優(yōu)的缺陷,改進(jìn)外部存檔更新機(jī)制以保存歷代優(yōu)秀的多樣性個體,采用克隆、高斯突變策略強(qiáng)化局部開采能力,有效地迫使算法跳出停滯搜索狀態(tài)。為應(yīng)對大規(guī)模約束,提出逐步微調(diào)機(jī)組出力策略,提高進(jìn)化群體的可行性。數(shù)值仿真實(shí)驗(yàn)以10機(jī)系統(tǒng)為測試算例,將HIMOA與著名的六種算法MODE、NSGA-Ⅱ、IMOEA/D-CH、ADEA、MOHDE_SAT、MONNDE進(jìn)行比較分析,結(jié)果表明,HIMOA能為DEED問題的10機(jī)系統(tǒng)提供較好的Pareto解,所獲的Pareto前沿收斂性和分布性優(yōu)越于其他算法,各評價指標(biāo)的箱型圖表明HIMOA具有優(yōu)越于其他算法的統(tǒng)計(jì)特征。

    關(guān)鍵詞:動態(tài)經(jīng)濟(jì)環(huán)境調(diào)度; 免疫系統(tǒng); 克隆選擇原理; 大規(guī)模約束; 多目標(biāo)優(yōu)化

    中圖分類號:TP304.7?? 文獻(xiàn)標(biāo)志碼:A

    文章編號:1001-3695(2023)09-024-2720-09

    doi:10.19734/j.issn.1001-3695.2023.02.0037

    Hybrid immune multi-objective optimization algorithm fordynamic emission economic dispatch

    Tang Xiangqian1,2, Qian Shuqu2, Wu Huihong2

    (1.College of Mathematics & Statistics, Guizhou University, Guiyang 550025, China; 2.School of Mathematics & Computer Sciences, Anshun University, Anshun Guizhou 561000, China)

    Abstract:Dynamic economic emission dispatch(DEED) problem in power system is a kind of high-dimensional multi-objective optimization problem with large-scale constraints. The traditional evolutionary algorithm is easy to fall into local optimization, and the distribution and convergence of the obtained Pareto frontier are poor. This paper fully explored the clonal selection principle of immune system, and proposed a hybrid immune multi-objective optimization algorithm(HIMOA) . The proposed algorithm took the traditional evolutionary algorithm as the basic framework. Since the deficiency of falling local optimization for the high-dimensional decision variables optimization, it improved the external archive update mechanism to preserve the excellent individuals from previous generations, and adopted the cloning and Gaussian mutation strategy to strengthen the local exploitation ability, forcing effectively the HIMOA to jump out of the stagnant search state. In order to cope with the large-scale constraints,it designed the fine-turning output power step by step to improve the feasibility of the evolutionary population. In the numerical simulation experiment, taking a 10-unit system as a test example, and compared HIMOA with the famous algorithms MODE, NSGA-Ⅱ, IMOEA/D-CH, ADEA, MOHDE_SAT, MONNDE. The results show that HIMOA can provide a better Pareto solution to the 10-unit system of the DEED problem, and the convergence and distribution of the Pareto frontier obtained are better than other algorithms. The box diagram of each evaluation metric shows that HIMOA has better statistical characteristics than other algorithms.

    Key words:dynamic economic emission dispatch; immune system; clonal selection principle; large-scale constraints; multi-objective optimization

    0 引言

    動態(tài)經(jīng)濟(jì)調(diào)度(dynamic economic dispatch,DED)是電力系統(tǒng)調(diào)度中一類復(fù)雜的單目標(biāo)優(yōu)化問題,該問題的目標(biāo)是確定調(diào)度周期內(nèi)各機(jī)組的最優(yōu)輸出功率[1,2],以滿足負(fù)載需求及線路功率傳輸損失。由于DED考慮了發(fā)電機(jī)組的爬坡率約束,使得其非常適合實(shí)際電力調(diào)度場景,這對電力公司的機(jī)組調(diào)度具有較大的指導(dǎo)價值。然而,DED未考慮發(fā)電機(jī)組在燃煤過程中的污染排放,這與企業(yè)和政府越來越關(guān)注的人類健康環(huán)境問題背道而馳。因此,近年來許多研究工作將污染排放納入動態(tài)經(jīng)濟(jì)調(diào)度模型的構(gòu)建,從而衍生出一種更為復(fù)雜的動態(tài)經(jīng)濟(jì)環(huán)境調(diào)度(dynamic economic emission dispatch,DEED)模型。該模型同時優(yōu)化發(fā)電成本和污染排放,并考慮機(jī)組爬坡率、負(fù)荷需求等約束[3,4],使得該問題屬一類含大規(guī)模約束的高維多目標(biāo)優(yōu)化問題。鑒于該問題的復(fù)雜性和實(shí)用性,開發(fā)高效的求解算法以獲得滿意的調(diào)度方案顯得尤為重要,這對實(shí)際工程問題的求解和實(shí)現(xiàn)我國“雙碳”目標(biāo)具有重要的理論和現(xiàn)實(shí)意義。

    求解多目標(biāo)DEED模型的傳統(tǒng)優(yōu)化方法是通過權(quán)重法將其轉(zhuǎn)換為單目標(biāo)問題,采用拉格朗日乘子[5]、梯度[6]、二次規(guī)劃[7]和動態(tài)規(guī)劃[8]等數(shù)學(xué)規(guī)劃方法獲取最優(yōu)調(diào)度方案。這些方法充分利用梯度信息搜索起始點(diǎn)附近的最優(yōu)解。然而,由于機(jī)組頻繁啟停產(chǎn)生損耗劇增的閥點(diǎn)效應(yīng),考慮閥點(diǎn)效應(yīng)成本的DEED模型具有多個局部最優(yōu)解且目標(biāo)函數(shù)非凸和非光滑,而數(shù)學(xué)優(yōu)化方法對選取的起始點(diǎn)極為敏感,往往難以有效處理這類含多個局部最優(yōu)解的復(fù)雜問題優(yōu)化。所以,許多智能優(yōu)化方法被開發(fā)來解決DEED問題,如進(jìn)化算法(evolutionary algorithm,EA)[9]、粒子群優(yōu)化(particle swarm optimization,PSO)[10,11]、差分進(jìn)化(differential evolutionary,DE)[12,13]、和諧搜索算法[14]、教學(xué)優(yōu)化[15]、生物地理優(yōu)化[16]等。在文獻(xiàn)[9]中,基于非支配排序遺傳算法Ⅱ(NSGA-Ⅱ),提出一種二次多項(xiàng)式求根法獲得滿足功率平衡約束的機(jī)組功率分配策略。然而,該方法獲得的第N個發(fā)電機(jī)的功率水平仍可能超出機(jī)組的功率上下界。文獻(xiàn)[12]針對DEED問題開發(fā)了一種多目標(biāo)差分進(jìn)化(multi-objective differential evolutionary,MODE)算法,采用上述約束處理機(jī)制,以10機(jī)系統(tǒng)測試表明,MODE所獲的Pareto前沿(pareto frontier,PF)優(yōu)于NSGA-Ⅱ。文獻(xiàn)[17]提出了一種模糊自適應(yīng)修正θ粒子群優(yōu)化算法。文獻(xiàn)[18]設(shè)計(jì)了一種改進(jìn)的θ多目標(biāo)教學(xué)優(yōu)化算法,利用模糊聚類技術(shù)對知識庫的規(guī)模進(jìn)行處理,從決策者的角度獲得可行的調(diào)度方案。文獻(xiàn)[2]開發(fā)了一種非優(yōu)勢DE算法,采用比例動態(tài)調(diào)整變量法調(diào)節(jié)機(jī)組功率范圍的差異以處理功率平衡約束。文獻(xiàn)[19]提出一種基于精英克隆的局部搜索DE算法,為提高算法開采和探索能力,增設(shè)精英群的克隆和突變機(jī)制,有效增強(qiáng)算法的全局搜索能力。文獻(xiàn)[20]提出了一種自適應(yīng)多目標(biāo)差分進(jìn)化算法,通過設(shè)計(jì)自適應(yīng)的current to best/1交叉算子,有效地提高了傳統(tǒng)進(jìn)化算法的探索與開采能力。

    雖然智能優(yōu)化方法能獲得DEED問題的有效解,但仍存在一些亟待解決的問題,如所獲Pareto解的收斂性、分布性和多樣性難于協(xié)調(diào)統(tǒng)一等問題。正如文獻(xiàn)[21]所述,混合策略可能是解決復(fù)雜的DEED問題的更有效方法。文獻(xiàn)[22]提出一種基于克隆選擇的PSO算法,以傳統(tǒng)的PSO算法為基本框架,嵌入克隆選擇算子,從而加速算法的收斂性。本文為克服傳統(tǒng)EA難以處理多個局部最優(yōu)解以及收斂性差等問題,提出一種混雜免疫多目標(biāo)優(yōu)化算法(hybrid immune multiobjective optimization algorithm,HIMOA)解決帶有閥點(diǎn)效應(yīng)的DEED問題。處理該問題的難點(diǎn)在于決策空間高維性和約束的大規(guī)模性,致使已有算法直接求解易于陷入局部搜索或無法獲取可行的非支配解,故解決該問題的策略是提高算法的搜索和開采能力,以獲取更多的可行解,而克隆選擇機(jī)制具有局部搜索和開采能力,能利用克隆機(jī)制產(chǎn)生大量的有潛力的解,并通過高斯變量擾動后,提高開采復(fù)雜區(qū)域的能力,同時采用修復(fù)機(jī)制對不可行的個體進(jìn)行修復(fù),從而提升種群中個體的可行性。故在HIMOA中,為了提升算法在高維決策空間的收斂速度,充分保存歷代多樣性的非支配解,選取存檔中優(yōu)秀個體進(jìn)行克隆,采用改進(jìn)的高斯突變策略對克隆體突變,實(shí)現(xiàn)收斂與開采的平衡。同時改進(jìn)Deb[23]的可行性規(guī)則對存檔中的個體進(jìn)行更新,確保存檔中非支配個體的多樣性和分布均勻性。另外,為了應(yīng)對功率平衡等式約束,提出一種逐步修復(fù)策略微調(diào)部分機(jī)組的出力,以提高群體中個體的可行性。數(shù)值實(shí)驗(yàn)以DEED問題的10機(jī)系統(tǒng)測試HIMOA的性能。結(jié)果表明,與同類算法相比,HIMOA獲得PF的分布均勻性和延展性更具有優(yōu)勢,且近似于真實(shí)PF的程度更高。

    本文的主要貢獻(xiàn)是:a)構(gòu)建了求解具有閥點(diǎn)效應(yīng)的DEED問題的混雜多目標(biāo)免疫優(yōu)化算法框架;b)研究了HIMOA在處理大規(guī)模約束時的運(yùn)行機(jī)制;c)設(shè)計(jì)了功率平衡方程的修復(fù)策略,以獲得更多可行的解;d)以考慮閥點(diǎn)效應(yīng)成本的10機(jī)系統(tǒng)為算例驗(yàn)證了HIMOA的求解能力,并與同類算法進(jìn)行了比較,驗(yàn)證了本文算法的優(yōu)越性。

    1 DEED問題模型

    電力系統(tǒng)DEED模型屬一類含大規(guī)模約束的多目標(biāo)優(yōu)化問題,其目標(biāo)是在滿足大量的等式和不等式約束下,使發(fā)電成本和污染排放盡可能小。

    然而,由HIMOA、MODE和NSGA-Ⅱ獲得的最好折中值是不可比較的,但由HIMOA的污染排放2.9774×105相對于MODE的污染排放2.9642×105增加了1.32×103個單位,增加率為0.44%小于1,而HIMOA的成本2.5188×106相對于MODE的成本2.5463×106減少了2.75×104個單位,降低率為1.09%大于1。類似地,對于NSGA-Ⅱ,HIMOA的污染排放2.9774×105相對于NSGA-Ⅱ的污染排放2.9588×105增加了1.86×103個單位,增加率為0.63%,而HIMOA的成本2.5188×106相對于NSGA-Ⅱ的成本2.5414×106減少了2.26×104個單位,降低率為0.90%,HMIOA所獲的折中解以增大少量的污染排放換取更低的發(fā)電成本。

    為了可視化各算法25次獨(dú)立執(zhí)行所獲的最好折中解在目標(biāo)空間的分布情況,圖2給出了基于模糊決策方法的各算法在25次運(yùn)行中得到的最好折中解分布,其中approximative PF為近似的Pareto前沿(PF)。由圖2(d)可以看出,HIMOA得到的25個最好折中解非常接近于近似的PF中心位置,這表明HIMOA每次執(zhí)行能獲得非常好的收斂性和分布性。然而,由圖2(a)和(c)可以看出,MODE和NSGA-Ⅱ得到的最好折中解偏向于成本目標(biāo)方向,特別是MODE的25個最好折中解遠(yuǎn)離近似的PF中心位置,且接近近似的PF程度時遠(yuǎn)時近,這表明MODE處理大規(guī)模約束DEED問題的穩(wěn)定性和收斂性非常差。觀察圖2(b)發(fā)現(xiàn),IMOEA/D-CH獲得的25個最好折中解能位于近似的PF中部,然而其分布空間區(qū)域非常大且遠(yuǎn)離近似PF,這表明IMOEA/D-CH處理該問題時在某些次獨(dú)立運(yùn)行中能獲得較好的收斂性,但一致收斂性差、穩(wěn)定性弱。由上述比較分析表明,新提出的HIMOA在處理該問題時能表現(xiàn)出較穩(wěn)定的搜索能力和較強(qiáng)的收斂性能。

    為了觀察各算法所獲的PF效果,選取25次獨(dú)立執(zhí)行獲得的PF中IGD指標(biāo)值最小的一次PF進(jìn)行比較。圖3展示了各算法所獲的PF分布。

    由圖3(d)可以明顯看出,HIMOA得到的PF在目標(biāo)空間中呈現(xiàn)出均勻的分布性,且分布范圍非常廣,所獲的PF接近近似PF的程度非常高,這表明HIMOA所獲的Pareto解多樣性和收斂性非常好。由圖3(a)和(c)可以看出,MODE和NSGA-Ⅱ所獲的PF分布范圍非常狹窄,這是由于該算法所獲的最后Pareto解收斂于某一區(qū)域,陷入了局部最優(yōu),非支配解的多樣性非常差,表明這兩種算法難以處理這類復(fù)雜的優(yōu)化問題,特別是MODE在接近近似PF方面表現(xiàn)非常差,這表明其對復(fù)雜的解空間開采能力弱,進(jìn)化種群發(fā)生早熟現(xiàn)象。觀察圖3(b)發(fā)現(xiàn),IMOEA/D-CH所獲的PF分布性具有一定的均勻性,但PF的分布范圍相對于HIMOA顯得短窄,且在成本目標(biāo)方向收斂性優(yōu)于污染排放目標(biāo)方向,這表明IMOEA/D-CH不能很好地平衡兩個目標(biāo)收斂性。由以上分析表明,HIMOA呈現(xiàn)較強(qiáng)的搜索和開采能力,能獲得分布均勻且范圍廣和收斂性優(yōu)越的PF,NSGA-Ⅱ和MODE的性能較差,僅能獲得局部的PF,且收斂性弱。

    為了從不同角度分析各算法25次獨(dú)立執(zhí)行所獲性能指標(biāo)HV、IGD、CR和SP的統(tǒng)計(jì)性能,圖4給出這四種評價指標(biāo)的箱型圖。

    圖4(a)的HV箱型圖顯示,與MODE、IMOEA/D-CH和NSGA-Ⅱ相比,HIMOA實(shí)現(xiàn)超體積大的概率高,而MODE最差,這表明HIMOA所獲的PF接近真實(shí)PF的概率高,而MODE所獲的25次PF或遠(yuǎn)離真實(shí)PF,或PF的范圍狹窄。對于CR指標(biāo),圖4(b)顯示HIMOA獲得最好的CR統(tǒng)計(jì)特征,而NSGA-Ⅱ最差,這表明HMIOA所獲的PF覆蓋率大的概率高,而NSGA-Ⅱ所獲得的25次PF或狹窄,或僅能得到部分非支配解。圖4(c)箱型圖表明,HIMOA獲得IGD值較小的概率高,MODE獲得大的IGD值概率高,這表明HIMOA所獲得的PF均勻性和收斂性優(yōu)于其他算法的概率大,其他算法中MODE獲效果差的概率高。對于SP評價指標(biāo),圖4(d)表明獲SP最大的概率為IMOEA/D-CH,其次為HIMOA,這表明對于SP指標(biāo)IMOEA/D-CH比HIMOA優(yōu)越,這是由于IMOEA/D-CH所獲的PF廣度較HIMOA的概率高。

    為了給出每臺發(fā)電機(jī)組在HIMOA優(yōu)化下的功率輸出,選取HIMOA最好折中解時10臺機(jī)組的各時段輸出功率(表2),燃料成本為2.5188E+6(美元),污染排放為2.9774E+5(磅)。表2的最后列出了考慮傳輸損失的各時間段的傳輸損失(loss)。由表2可以看出,對于考慮傳輸損耗的10機(jī)系統(tǒng)的多目標(biāo)DEED優(yōu)化,各時段各機(jī)組的出力都得到了控制,各機(jī)組的爬坡率也得到了適當(dāng)?shù)臐M足。

    4.4 策略分析及討論

    為了分析算法HIMOA中改進(jìn)策略的優(yōu)越性,本節(jié)以10機(jī)系統(tǒng)DEED問題作為測試實(shí)例,分析不同的約束閾值、不同的變異策略、不同的存檔更新策略、有/無修復(fù)策略對算法收斂性能的影響。實(shí)驗(yàn)中HIMOA獨(dú)立執(zhí)行25次循環(huán),每次循環(huán)的最大迭代數(shù)為2 000,計(jì)算每代所獲PF與理論P(yáng)F的IGD平均值,形成IGD隨迭代數(shù)的變化曲線,分析改進(jìn)的策略對算法收斂性和所獲PF上點(diǎn)的分布性的影響。

    圖5給出了約束閾值ε取0、0.05、0.1、0.15和0.2五種情況時HIMOA對10機(jī)系統(tǒng)優(yōu)化所獲的IGD隨迭代數(shù)的變化曲線。IGD是刻畫算法所獲PF收斂真實(shí)PF的程度,由圖5的IGD變化曲線獲知,ε取0時收斂性能最差,但ε在0.05~0.2較為相似,具體為:在ε取0.1時表現(xiàn)初始收斂速度快,隨后陷入局部搜索;ε取0.15和0.2時初始收斂速度稍微慢,但最終與ε取0.05表現(xiàn)相當(dāng)。圖6分析了對克隆體采用三種不同的變異策略獲得的IGD隨迭代數(shù)的變化曲線。由圖6可知,非均勻突變和多項(xiàng)式突變策略在初始階段表現(xiàn)出快速的收斂,但非均勻突變到500代左右陷入了局部搜索,而多項(xiàng)式突變到600代左右收斂速度減慢,隨后下降速度緩慢并一直保持平衡趨勢。高斯突變初始下降速度具有一定的波動,并表現(xiàn)弱于其他兩種變異策略,在1 000代左右所獲的IGD值開始低于其他兩種變異策略。這是由于本文采取的高斯突變表達(dá)式為uij=wij+1/g·N(0,1),其中第2項(xiàng)中1/g使得變異量隨迭代數(shù)增大而變小,動態(tài)調(diào)整該系數(shù)使得算法具有精準(zhǔn)的探索能力,N(0,1)為均值為0,方差為1的高斯變量。

    圖7給出了NSGA-Ⅱ中原始存檔更新策略所獲的IGD變化曲線與本文改進(jìn)的存檔更新策略下的IGD變化曲線。由變化曲線明顯可知,改進(jìn)的存檔更新策略獲得了更好的收斂速度,這表明改進(jìn)的存檔更新策略能有效提高算法對10機(jī)系統(tǒng)DEED問題的優(yōu)化效果。這是由于改進(jìn)的存檔更新策略既能加大算法對可能空間的探索與開采能力,又能提升存檔中非支配解的多樣性,確保存檔中非支配解的均勻分布。

    圖8給出了25次獨(dú)立執(zhí)行所獲PF分布。對比圖8(a)和(b)明顯可以看出,改進(jìn)的更新策略25次所獲PF均能保持很好的收斂性和分布性,而原始的更新策略僅能獲得局部的PF,每次執(zhí)行所獲PF分布性和收斂性不統(tǒng)一。

    圖9分析了修復(fù)策略對算法HIMOA搜索性能的影響。其中圖9(a)比較了有/無修復(fù)策略情況下進(jìn)化種群中個體的非可行性比率隨迭代數(shù)的變化曲線。由圖9(a)可知,無修復(fù)策略時,種群中個體的非可行性比率隨迭代數(shù)增加呈現(xiàn)下降趨勢,但運(yùn)行到2 000代時仍未達(dá)到0,這表明種群中還存在不可行的個體。而含修復(fù)策略情況下,在算法運(yùn)行到1 000代左右時,個體的非可行性比率已達(dá)到0,這充分表明修復(fù)算子對提升個體的可行性比率有重大影響。圖9(b)為有/無修復(fù)策略所獲IGD隨迭代數(shù)的變化曲線。由圖9可知,含修復(fù)策略情況下算法所獲的IGD變化曲線明顯下降速度快,這表明算法的收斂性能好。

    5 結(jié)束語

    動態(tài)環(huán)境經(jīng)濟(jì)調(diào)度優(yōu)化是計(jì)算機(jī)、管理學(xué)交叉學(xué)科中一類極具挑戰(zhàn)性的研究課題,該問題具有多個局部最優(yōu)點(diǎn),且目標(biāo)函數(shù)非光滑和非凸,傳統(tǒng)的進(jìn)化算法難以尋找全局最優(yōu)解。本文提出了一種混合免疫多目標(biāo)優(yōu)化算法(HIMOA),并以10機(jī)系統(tǒng)的DEED問題進(jìn)行測試。在算法設(shè)計(jì)中,采用修復(fù)方案處理功率平衡等式約束,提高了進(jìn)化種群的可行性。為了增強(qiáng)傳統(tǒng)EA的收斂能力,設(shè)計(jì)外部存檔集保存非支配個體并提出新的裁剪策略更新外部存檔,選擇存檔中優(yōu)秀個體進(jìn)行克隆,然后采用高斯策略對克隆體進(jìn)行突變。將該算法作為一種搜索方法處理10機(jī)系統(tǒng)DEED問題,實(shí)驗(yàn)結(jié)果表明,相對于最新文獻(xiàn)報道的方法,本文算法能夠獲得更好的性能,評價指標(biāo)統(tǒng)計(jì)結(jié)果表明本文算法具有更好的統(tǒng)計(jì)性能。

    與此同時,本文的工作還有待進(jìn)一步的研究,以便應(yīng)用于更復(fù)雜的DEED模型求解,如將所開發(fā)的框架引入模擬退火、蟻群、魚群等其他搜索機(jī)制,應(yīng)用于含其他能源的混雜DEED問題。另外,本文只在HIMOA中選取了一組固定的算法參數(shù),而沒有對本文算法進(jìn)行參數(shù)調(diào)節(jié)或自適應(yīng)設(shè)計(jì),這將為開發(fā)自適應(yīng)算法提供后期研究思路。

    參考文獻(xiàn):

    [1]Elattar E E. A hybrid genetic algorithm and bacterial foraging approach for dynamic economic dispatch problem [J]. International Journal of Electrical Power & Energy Systems, 2015,69:18-26.

    [2]Qiao Baihao, Liu Jing, Hao Xingxing. A multi-objective differential evolution algorithm and a constraint handling mechanism based on variables proportion for dynamic economic emission dispatch problems[J]. Applied Soft Computing, 2021,108: 107419.

    [3]Jebaraj L, Venkatesan C, Soubache I, et al. Application of differential evolution algorithm in static and dynamic economic or emission dispatch problem: a review[J]. Renewable & Sustainable Energy Reviews, 2021,77: 1206-1220.

    [4]Li Lingling, Liu Zhifeng, Tseng M L, et al. Improved tunicate swarm algorithm: solving the dynamic economic emission dispatch problems[J]. Applied Soft Computing, 2021,108: 107504.

    [5]Yang Lingjian, Fraga E S, Papageorgiou L G. Mathematical programming formulations for nonsmooth and non-convex electricity dispatch problems[J]. Electric Power Systems Research, 2013,95: 302-308.

    [6]Wood A J, Wollenberg B. Power generation operation and control[J]. Fuel & Energy Abstracts, 1996,37(3): 90-93.

    [7]Coelho L S, Mariani V C. Combining of chaotic differential evolution and quadratic programming for economic dispatch optimization with valve-point effect[J]. IEEE Trans on Power Systems, 2021,21(2): 989-996.

    [8]Ross D W, Kim S. Dynamic economic dispatch of generation[J]. IEEE Trans on Power Apparatus & Systems, 2007,99(6): 2060-2068.

    [9]Basu M. Dynamic economic emission dispatch using nondominated sorting genetic algorithm-Ⅱ[J]. International Journal of Electrical Power & Energy Systems, 2008,30(2): 140-149.

    [10]Zhang Yong , Gong Dunwei, Geng Na, et al. Hybrid bare-bones PSO for dynamic economic dispatch with valve-point effects[J]. Applied Soft Computing Journal, 2014,18:248-260.

    [11]Mason K, Duggan J, Howley E. Multi-objective dynamic economic emission dispatch using particle swarm optimisation variants[J]. Neurocomputing, 2017,270: 188-197.

    [12]Basu M. Multi-objective differential evolution for dynamic economic emission dispatch[J]. International Journal of Emerging Electric Power Systems, 2014,15(2): 141-150.

    [13]Zhang Huifeng, Yue Dong, Xie Xiangpeng, et al. Multi-elite guide hybrid differential evolution with simulated annealing technique for dynamic economic emission dispatch[J]. Applied Soft Computing, 2015,34: 312-323.

    [14]Fakhim-Babaei A, Nazari-Heris M, Mohammadi-Ivatloo B, et al. Economic dispatch of large-scale integrated heat and power systems by application of a novel harmony search approach[J]. Hybrid Energy System Models, 2021,52(3): 25-49.

    [15]Yang Zhile, Li Kang, Niu Qun, et al. A self-learning TLBO based dynamic economic/environmental dispatch considering multiple plu-gin electric vehicle loads[J]. Journal of Modern Power Systems & Clean Energy, 2014,2(4): 298-307.

    [16]Ma Haiping, Yang Zhile, You Pengcheng, et al. Multi-objective biogeography-based optimization for dynamic economic emission load dispatch considering plug-in electric vehicles charging[J]. Energy, 2017,135: 102-111.

    [17]Bahmanifirouzi B, Farjah E, Niknam T. Multi-objective stochastic dynamic economic emission dispatch enhancement by fuzzy adaptive modified theta particle swarm optimization[J]. Journal of Renewable & Sustainable Energy, 2012,4(2): 013102-6.

    [18]Niknam T, Golestaneh F, Sadeghi M S. θ-multiobjective teaching learning-based optimization for dynamic economic emission dispatch[J]. IEEE Systems Journal, 2012,6(2): 341-352.

    [19]武慧虹, 錢淑渠, 劉衍民, 等. 精英克隆局部搜索的多目標(biāo)動態(tài)環(huán)境經(jīng)濟(jì)調(diào)度差分進(jìn)化算法[J]. 山東大學(xué)學(xué)報: 工學(xué)版, 2021,51(1):11-23. (Wu Huihong, Qian Shuqu, Liu Yanmin, et al. Multiobjective dynamic economic emission dispatch differential evolution algorithm based on elites cloning local search[J]. Journal of Shandong University:Engineering Science, 2021,51(1): 11-23.)

    [20]武慧虹, 林妤, 曾茜, 等. 自適應(yīng)差分進(jìn)化算法及對動態(tài)環(huán)境經(jīng)濟(jì)調(diào)度問題應(yīng)用[J]. 計(jì)算機(jī)應(yīng)用研究, 2021,38(5):1443-1448,1454. (Wu Huihong, Lin Yu, Zeng Xi, et al. Adaptive diffe-rent evolutionary algorithm and its application of dynamic emission economic dispatch problem[J]. Application Research of Compu-ters, 2021,38(5): 1443-1448,1454.)

    [21]Qian Shuqu, Wu Huihong, Xu Guofeng. An improved particle swarm optimization with clone selection principle for dynamic economic emission dispatch[J]. Soft Computing, 2020,24(20): 15249-15271.

    [22]Elaiw A M, Xia X, Shehata A M. Hybrid DE-SQP and hybrid PSO-SQP methods for solving dynamic economic emission dispatch problem with valve-point effects[J]. Electric Power Systems Research, 2013,103: 192-200.

    [23]Deb K. An efficient constraint handling method for genetic algorithms[J]. Computer Methods in Applied Mechanics & Engineering, 2000,186(24): 311-338.

    [24]De Castro L N, Von Zuben F J. Learning and optimization using the clonal selection principle[J]. IEEE Trans on Evolutionary Computation, 2002,6(3): 239-251.

    [25]Panigrahi B K, Yadav S R, Agrawal S, et al. A clonal algorithm to solve economic load dispatch[J]. Electric Power Systems Research, 2007,77(10): 1381-1389.

    [26]El-Sharkh M Y. Clonal selection algorithm for power generators maintenance scheduling[J]. International Journal of Electrical Power & Energy Systems, 2014,57(5): 73-78.

    [27]Aragon V S, Esquivel S C, Coello C C. An immune algorithm with power redistribution for solving economic dispatch problems[J]. Information Sciences, 2015,295: 609-632.

    [28]Deb K, Pratap A, Agarwal S, et al. A fast and elitist multiobjective genetic algorithm: NSGA-Ⅱ[J]. IEEE Trans on Evolutionary Computation, 2002,6(2): 182- 197.

    [29]Zhu Yongsheng, Qiao Baihao, Dong Yan, et al. Multiobjective dynamic economic emission dispatch using evolutionary algorithm based on decomposition[J]. IEEE Trans on Electrical and Electronic Engineering, 2019,14(9): 1323-1333.

    [30]Kukkonen S, Deb K. Improved pruning of non-dominated solutions based on crowding distance for bi-objective optimization problems[C]//Proc of IEEE International Conference on Evolutionary Computation. Piscataway, NJ: IEEE Press, 2006: 1179-1186.

    [31]Qian Shuqu, Ye Yongqiang, Jiang Bin, et al. Constrained multi-objective optimization algorithm based on immune system model[J]. IEEE Trans on Cybernetics, 2016,46(9): 2056-2069.

    [32]Mason K, Duggan J, Howley E. A multi-objective neural network trained with differential evolution for dynamic economic emission dispatch[J]. International Journal of Electrical Power & Energy Systems, 2018,100: 201-221.

    [33]Basu M. Particle swarm optimization based goal-attainment method for dynamic economic emission dispatch[J]. Electric Power Components & Systems, 2006,34(9): 1015-1025.

    [34]Wang Weiwei, Dai Shijuan, Zhao Wanzhong, et al. Multi-objective optimization of hexahedral pyramid crash box using MOEA/D-DAE algorithm[J]. Applied Soft Computing, 2022,118: 108481.

    [35]Zitzler E, Thiele L, Laumanns M, et al. Performance assessment of multiobjective optimizers: an analysis and review[J]. IEEE Trans on Evolutionary Computation, 2003,7(2): 117-132.

    [36]Basu M. Dynamic economic emission dispatch using evolutionary programming and fuzzy satisfying method[J]. International Journal of Emerging Electric Power Systems, 2007,8(4): 14-19.

    收稿日期:2023-02-02;修回日期:2023-03-27? 基金項(xiàng)目:國家自然科學(xué)基金資助項(xiàng)目(62241301,61762001);貴州省教育廳創(chuàng)新群體重大研究資助項(xiàng)目(黔教合KY字[2019]069);貴州省教育廳青年科技人才成長資助項(xiàng)目(黔教合KY字[2020]131號);安順學(xué)院研究生創(chuàng)新專項(xiàng)資助項(xiàng)目(asxysrt(202223)號)

    作者簡介:唐湘黔(1994-),男,貴州遵義人,碩士研究生,主要研究方向?yàn)橹悄軆?yōu)化算法及應(yīng)用;錢淑渠(1978-),男(通信作者),安徽樅陽人,教授,碩導(dǎo),博士,主要研究方向?yàn)橛?jì)算智能、系統(tǒng)建模與控制(shuquqian@163.com);武慧虹(1980-),女,山西太原人,教授,碩士,主要研究方向?yàn)橹悄軆?yōu)化算法、群與圖.

    猜你喜歡
    多目標(biāo)優(yōu)化免疫系統(tǒng)
    讓免疫系統(tǒng)適當(dāng)“放松”
    身體的保護(hù)傘——免疫系統(tǒng)
    一類具有抗原性的腫瘤-免疫系統(tǒng)的定性分析
    保護(hù)好你自己的免疫系統(tǒng)
    Staying healthy
    神經(jīng)-免疫系統(tǒng)相互作用的分子機(jī)制
    改進(jìn)的多目標(biāo)啟發(fā)式粒子群算法及其在桁架結(jié)構(gòu)設(shè)計(jì)中的應(yīng)用
    群體多目標(biāo)優(yōu)化問題的權(quán)序α度聯(lián)合有效解
    云計(jì)算中虛擬機(jī)放置多目標(biāo)優(yōu)化
    狼群算法的研究
    看免费成人av毛片| 少妇猛男粗大的猛烈进出视频 | 日日夜夜操网爽| 18禁裸乳无遮挡免费网站照片| 午夜a级毛片| 三级男女做爰猛烈吃奶摸视频| 国产精品电影一区二区三区| 免费看日本二区| 亚洲四区av| 久久亚洲真实| 欧美激情在线99| 床上黄色一级片| 亚洲精品一区av在线观看| aaaaa片日本免费| 久久午夜亚洲精品久久| 久久精品国产亚洲网站| av视频在线观看入口| 国产伦精品一区二区三区视频9| 亚洲人成网站在线播放欧美日韩| 欧美一级a爱片免费观看看| 在线免费十八禁| 舔av片在线| 欧美日韩精品成人综合77777| av福利片在线观看| 成人美女网站在线观看视频| 午夜爱爱视频在线播放| 黄色视频,在线免费观看| 一a级毛片在线观看| 国产视频一区二区在线看| 亚洲av中文字字幕乱码综合| 91精品国产九色| 国模一区二区三区四区视频| 国产精品一区二区三区四区久久| 日本免费一区二区三区高清不卡| 色5月婷婷丁香| 亚洲不卡免费看| 亚洲成人久久爱视频| 成人特级av手机在线观看| 色在线成人网| 久久精品人妻少妇| 亚洲精品日韩av片在线观看| 久久亚洲精品不卡| 久久婷婷人人爽人人干人人爱| 在线天堂最新版资源| 日韩一区二区视频免费看| 啦啦啦啦在线视频资源| 一个人观看的视频www高清免费观看| 99热网站在线观看| 欧美一区二区精品小视频在线| 亚洲av二区三区四区| 欧美日本视频| 99热精品在线国产| xxxwww97欧美| 日本 欧美在线| 哪里可以看免费的av片| 国产午夜精品久久久久久一区二区三区 | 久久欧美精品欧美久久欧美| 亚洲va日本ⅴa欧美va伊人久久| 天堂影院成人在线观看| 国产精品一区www在线观看 | 日日摸夜夜添夜夜添小说| 成人二区视频| 欧美成人一区二区免费高清观看| 日本精品一区二区三区蜜桃| 人人妻,人人澡人人爽秒播| 一级黄色大片毛片| 午夜视频国产福利| 蜜桃亚洲精品一区二区三区| 好男人在线观看高清免费视频| 国产精品女同一区二区软件 | 成年人黄色毛片网站| 国内毛片毛片毛片毛片毛片| 99久久中文字幕三级久久日本| 日日啪夜夜撸| 网址你懂的国产日韩在线| 国产激情偷乱视频一区二区| 午夜福利成人在线免费观看| 久久天躁狠狠躁夜夜2o2o| 九色成人免费人妻av| 国产真实伦视频高清在线观看 | 岛国在线免费视频观看| 亚洲国产欧美人成| www.色视频.com| 国产av在哪里看| 国产亚洲精品久久久久久毛片| 特级一级黄色大片| 国产成人一区二区在线| 国产大屁股一区二区在线视频| 欧美丝袜亚洲另类 | 在线a可以看的网站| 美女被艹到高潮喷水动态| 亚洲人成网站高清观看| 国产久久久一区二区三区| 成年女人毛片免费观看观看9| 国产美女午夜福利| 99久久精品国产国产毛片| 亚洲五月天丁香| 18禁在线播放成人免费| 亚洲五月天丁香| 能在线免费观看的黄片| 久久99热6这里只有精品| bbb黄色大片| 色尼玛亚洲综合影院| 亚洲最大成人中文| 国产白丝娇喘喷水9色精品| 国产乱人伦免费视频| 欧美激情国产日韩精品一区| 欧美+日韩+精品| 人人妻人人看人人澡| 大又大粗又爽又黄少妇毛片口| 亚洲图色成人| 国产视频内射| 久久久色成人| 亚洲五月天丁香| 国产精品自产拍在线观看55亚洲| 一进一出好大好爽视频| 免费看a级黄色片| 亚洲第一区二区三区不卡| 中国美白少妇内射xxxbb| 级片在线观看| 啦啦啦啦在线视频资源| 国产淫片久久久久久久久| 欧美成人性av电影在线观看| 长腿黑丝高跟| eeuss影院久久| 日本免费一区二区三区高清不卡| 国产精品98久久久久久宅男小说| 亚洲中文日韩欧美视频| 国产又黄又爽又无遮挡在线| 欧美丝袜亚洲另类 | 亚洲va日本ⅴa欧美va伊人久久| 午夜日韩欧美国产| 啦啦啦韩国在线观看视频| 天堂影院成人在线观看| 日本在线视频免费播放| 午夜福利成人在线免费观看| 亚洲精品456在线播放app | 日韩精品中文字幕看吧| 免费在线观看日本一区| 欧美成人免费av一区二区三区| 亚洲四区av| 一区二区三区四区激情视频 | 九色国产91popny在线| 亚洲精品乱码久久久v下载方式| 国产精品伦人一区二区| 一级黄片播放器| 别揉我奶头~嗯~啊~动态视频| 亚洲精品久久国产高清桃花| 99国产精品一区二区蜜桃av| 欧美zozozo另类| 免费观看人在逋| 国产精品爽爽va在线观看网站| 欧美精品国产亚洲| 天堂动漫精品| 亚洲av.av天堂| 精品久久久久久久久亚洲 | 日韩欧美精品免费久久| 男人狂女人下面高潮的视频| 免费搜索国产男女视频| 亚洲狠狠婷婷综合久久图片| 美女高潮喷水抽搐中文字幕| 国产高清激情床上av| 亚洲无线在线观看| 两个人的视频大全免费| 在线观看66精品国产| 国产成人一区二区在线| 精品久久久久久久人妻蜜臀av| 国产综合懂色| 欧美极品一区二区三区四区| 午夜激情福利司机影院| 麻豆国产av国片精品| 天天躁日日操中文字幕| 亚洲av第一区精品v没综合| 男女边吃奶边做爰视频| 免费电影在线观看免费观看| 国内揄拍国产精品人妻在线| 精品久久久久久,| 精品人妻偷拍中文字幕| 12—13女人毛片做爰片一| 全区人妻精品视频| 欧美不卡视频在线免费观看| 国产伦一二天堂av在线观看| 中文字幕高清在线视频| 精品国内亚洲2022精品成人| 国产又黄又爽又无遮挡在线| 免费黄网站久久成人精品| 亚洲人成伊人成综合网2020| 最近在线观看免费完整版| 观看美女的网站| 九九热线精品视视频播放| АⅤ资源中文在线天堂| 十八禁网站免费在线| 美女高潮喷水抽搐中文字幕| 丰满人妻一区二区三区视频av| 亚洲综合色惰| 日日撸夜夜添| 国产蜜桃级精品一区二区三区| 精华霜和精华液先用哪个| 亚洲成人中文字幕在线播放| 无遮挡黄片免费观看| 无遮挡黄片免费观看| 午夜福利18| 2021天堂中文幕一二区在线观| 亚洲精品一区av在线观看| 亚洲精品一区av在线观看| 亚洲乱码一区二区免费版| 成人三级黄色视频| 精品一区二区免费观看| 久久久久久伊人网av| 国产精品免费一区二区三区在线| 天堂动漫精品| 麻豆国产97在线/欧美| 精品久久国产蜜桃| 偷拍熟女少妇极品色| 少妇猛男粗大的猛烈进出视频 | 亚洲av中文av极速乱 | 精品一区二区三区视频在线| 国内精品久久久久久久电影| 欧美3d第一页| 中出人妻视频一区二区| 欧美极品一区二区三区四区| 精品人妻一区二区三区麻豆 | 99九九线精品视频在线观看视频| 嫩草影院精品99| 日韩欧美免费精品| 免费大片18禁| 最近中文字幕高清免费大全6 | 1024手机看黄色片| 男人和女人高潮做爰伦理| 国产淫片久久久久久久久| 成人av在线播放网站| 在线免费十八禁| 热99re8久久精品国产| 麻豆av噜噜一区二区三区| 啦啦啦啦在线视频资源| 精品人妻熟女av久视频| 麻豆成人午夜福利视频| 国产精华一区二区三区| 午夜精品在线福利| 日韩亚洲欧美综合| 久久久色成人| 久久99热6这里只有精品| 欧美成人免费av一区二区三区| 亚洲人成伊人成综合网2020| 熟女人妻精品中文字幕| 51国产日韩欧美| 国内精品久久久久久久电影| 国产单亲对白刺激| 国内久久婷婷六月综合欲色啪| 欧美国产日韩亚洲一区| 国产 一区 欧美 日韩| 久久久久久久亚洲中文字幕| 可以在线观看毛片的网站| 日韩欧美在线乱码| 久久久久国产精品人妻aⅴ院| 99久久精品一区二区三区| 真人做人爱边吃奶动态| 国产精品一区二区三区四区免费观看 | 日韩高清综合在线| 国产精品美女特级片免费视频播放器| 性欧美人与动物交配| 国产aⅴ精品一区二区三区波| 看黄色毛片网站| 少妇高潮的动态图| 欧美日韩黄片免| 神马国产精品三级电影在线观看| 一a级毛片在线观看| 美女高潮喷水抽搐中文字幕| 久久久久久大精品| 最近中文字幕高清免费大全6 | 美女cb高潮喷水在线观看| 日本爱情动作片www.在线观看 | 男女边吃奶边做爰视频| 美女xxoo啪啪120秒动态图| 嫩草影院入口| 别揉我奶头~嗯~啊~动态视频| 久久久国产成人精品二区| 亚洲男人的天堂狠狠| 男女那种视频在线观看| 免费av观看视频| 亚洲经典国产精华液单| 亚洲中文字幕日韩| 亚洲精品在线观看二区| 少妇被粗大猛烈的视频| 婷婷六月久久综合丁香| 丰满人妻一区二区三区视频av| 亚洲av日韩精品久久久久久密| 少妇裸体淫交视频免费看高清| 啪啪无遮挡十八禁网站| 久久国内精品自在自线图片| aaaaa片日本免费| 欧美国产日韩亚洲一区| av专区在线播放| 欧美激情国产日韩精品一区| 联通29元200g的流量卡| 天堂av国产一区二区熟女人妻| 直男gayav资源| 99热这里只有是精品50| 伦精品一区二区三区| 国产精品久久久久久久电影| 日韩精品有码人妻一区| 久久久久久久午夜电影| 真人一进一出gif抽搐免费| 免费黄网站久久成人精品| 国产精品不卡视频一区二区| 成年女人看的毛片在线观看| 午夜a级毛片| a级毛片免费高清观看在线播放| 久久国内精品自在自线图片| 精品午夜福利视频在线观看一区| a级一级毛片免费在线观看| 国产高清不卡午夜福利| 国产真实伦视频高清在线观看 | 成人二区视频| 淫妇啪啪啪对白视频| 亚洲欧美日韩东京热| 亚洲第一区二区三区不卡| 一级a爱片免费观看的视频| 波多野结衣高清作品| 制服丝袜大香蕉在线| 中国美女看黄片| 国内精品久久久久久久电影| 又黄又爽又刺激的免费视频.| 国产午夜精品论理片| 国产熟女欧美一区二区| 国产av在哪里看| 88av欧美| 免费人成视频x8x8入口观看| 国产精品美女特级片免费视频播放器| 欧美国产日韩亚洲一区| 亚洲无线观看免费| 99久久久亚洲精品蜜臀av| 免费电影在线观看免费观看| 欧美成人a在线观看| 国产v大片淫在线免费观看| 色综合婷婷激情| 一区二区三区高清视频在线| 国产私拍福利视频在线观看| 亚洲aⅴ乱码一区二区在线播放| xxxwww97欧美| 真人做人爱边吃奶动态| 动漫黄色视频在线观看| 亚洲美女视频黄频| 日本免费a在线| 人妻少妇偷人精品九色| 色尼玛亚洲综合影院| 亚洲无线观看免费| 老熟妇仑乱视频hdxx| 美女高潮喷水抽搐中文字幕| 欧美激情在线99| 精品欧美国产一区二区三| 桃色一区二区三区在线观看| 又爽又黄a免费视频| 日本熟妇午夜| 日本五十路高清| 久久久久久久午夜电影| 精品欧美国产一区二区三| 观看免费一级毛片| 99精品在免费线老司机午夜| 尾随美女入室| 午夜免费成人在线视频| 亚洲无线在线观看| 亚洲av日韩精品久久久久久密| 最近最新免费中文字幕在线| 久久久国产成人精品二区| 可以在线观看毛片的网站| 999久久久精品免费观看国产| 国产高清激情床上av| 欧美成人性av电影在线观看| 亚洲一区高清亚洲精品| 亚洲中文日韩欧美视频| 国内少妇人妻偷人精品xxx网站| 亚洲自偷自拍三级| 一区二区三区免费毛片| 99热这里只有精品一区| 又黄又爽又免费观看的视频| 午夜福利欧美成人| 日韩av在线大香蕉| 午夜福利成人在线免费观看| 国产精品人妻久久久久久| 日本-黄色视频高清免费观看| 久久精品国产亚洲av香蕉五月| 一个人免费在线观看电影| 很黄的视频免费| 国产精品人妻久久久久久| 日韩人妻高清精品专区| 长腿黑丝高跟| 在线免费十八禁| 赤兔流量卡办理| 日本五十路高清| 美女xxoo啪啪120秒动态图| 国产亚洲欧美98| 成人永久免费在线观看视频| 日韩欧美国产一区二区入口| 国产国拍精品亚洲av在线观看| 久久国内精品自在自线图片| 中文字幕人妻熟人妻熟丝袜美| 国产精品国产高清国产av| 精品久久久久久久久亚洲 | 天堂网av新在线| 精品一区二区三区视频在线观看免费| 免费不卡的大黄色大毛片视频在线观看 | av视频在线观看入口| 成熟少妇高潮喷水视频| 国产毛片a区久久久久| 黄片wwwwww| 亚洲五月天丁香| 欧美另类亚洲清纯唯美| 精品福利观看| 亚洲美女视频黄频| 91在线精品国自产拍蜜月| 女人被狂操c到高潮| 亚洲,欧美,日韩| 又粗又爽又猛毛片免费看| 国产av不卡久久| 久久午夜亚洲精品久久| 精品久久久噜噜| 99久久久亚洲精品蜜臀av| 国产国拍精品亚洲av在线观看| 乱系列少妇在线播放| 亚洲成av人片在线播放无| 日本色播在线视频| 精品免费久久久久久久清纯| 亚洲av不卡在线观看| 天天一区二区日本电影三级| 男女下面进入的视频免费午夜| 亚洲黑人精品在线| 别揉我奶头 嗯啊视频| 免费av不卡在线播放| 丰满的人妻完整版| 日韩亚洲欧美综合| 天堂av国产一区二区熟女人妻| 国产三级中文精品| 美女xxoo啪啪120秒动态图| 简卡轻食公司| 日韩大尺度精品在线看网址| 免费在线观看影片大全网站| 亚洲熟妇熟女久久| 美女 人体艺术 gogo| 在线播放无遮挡| 深夜a级毛片| 欧美xxxx黑人xx丫x性爽| 乱码一卡2卡4卡精品| 国产精品一区二区性色av| 欧美一区二区国产精品久久精品| 国产日本99.免费观看| 久久精品国产亚洲av涩爱 | 国产精品久久视频播放| 中文在线观看免费www的网站| 欧美区成人在线视频| 亚洲自偷自拍三级| 亚洲成人久久爱视频| 色5月婷婷丁香| 色噜噜av男人的天堂激情| 天天躁日日操中文字幕| 又黄又爽又刺激的免费视频.| 免费观看在线日韩| 色噜噜av男人的天堂激情| 欧美成人一区二区免费高清观看| 成人欧美大片| 中文亚洲av片在线观看爽| 免费高清视频大片| xxxwww97欧美| 欧美最新免费一区二区三区| 国产av不卡久久| 啦啦啦啦在线视频资源| 国产精品98久久久久久宅男小说| 天堂动漫精品| 免费在线观看成人毛片| 伦精品一区二区三区| 国产 一区精品| 禁无遮挡网站| 麻豆一二三区av精品| 一级a爱片免费观看的视频| 18禁裸乳无遮挡免费网站照片| 99精品在免费线老司机午夜| 亚洲精品在线观看二区| 天天躁日日操中文字幕| 亚洲欧美日韩高清在线视频| 亚洲精品在线观看二区| 色播亚洲综合网| 免费观看在线日韩| 黄色视频,在线免费观看| 亚洲国产精品久久男人天堂| 国产极品精品免费视频能看的| 国产精品自产拍在线观看55亚洲| 特大巨黑吊av在线直播| 永久网站在线| 美女xxoo啪啪120秒动态图| 国产精品久久久久久久久免| 久久精品91蜜桃| av在线天堂中文字幕| 亚洲图色成人| 国产女主播在线喷水免费视频网站 | 18禁在线播放成人免费| 级片在线观看| 婷婷精品国产亚洲av在线| 久久午夜福利片| 国产激情偷乱视频一区二区| 黄片wwwwww| 亚洲国产色片| 亚洲色图av天堂| 亚洲av成人av| 亚洲av不卡在线观看| 日本黄大片高清| 欧美日韩乱码在线| 又爽又黄a免费视频| 男女啪啪激烈高潮av片| 色精品久久人妻99蜜桃| 午夜免费男女啪啪视频观看 | 国产乱人视频| 男女视频在线观看网站免费| 在线播放无遮挡| 哪里可以看免费的av片| 十八禁国产超污无遮挡网站| 91麻豆精品激情在线观看国产| 男人狂女人下面高潮的视频| 中文字幕av成人在线电影| 亚洲五月天丁香| 麻豆一二三区av精品| 男女边吃奶边做爰视频| 欧美一级a爱片免费观看看| 最新中文字幕久久久久| 麻豆av噜噜一区二区三区| 成人av一区二区三区在线看| 观看美女的网站| 欧美潮喷喷水| 亚洲熟妇中文字幕五十中出| 精品久久久久久成人av| 日韩精品青青久久久久久| 禁无遮挡网站| 免费无遮挡裸体视频| 在线观看美女被高潮喷水网站| 亚洲在线观看片| 婷婷色综合大香蕉| 国产麻豆成人av免费视频| av天堂中文字幕网| 亚洲欧美日韩卡通动漫| 日韩欧美在线乱码| 国产精品一及| 精品99又大又爽又粗少妇毛片 | 色精品久久人妻99蜜桃| 国产精品电影一区二区三区| 国产精品福利在线免费观看| 色播亚洲综合网| 18+在线观看网站| 国产高清视频在线观看网站| 免费黄网站久久成人精品| 国产三级在线视频| 狂野欧美白嫩少妇大欣赏| 国产精品久久久久久久电影| 一个人观看的视频www高清免费观看| 亚洲成av人片在线播放无| 亚洲av成人av| av黄色大香蕉| 男女下面进入的视频免费午夜| 女人被狂操c到高潮| 91午夜精品亚洲一区二区三区 | 少妇熟女aⅴ在线视频| 亚洲性夜色夜夜综合| 亚洲人与动物交配视频| 国产女主播在线喷水免费视频网站 | 成年人黄色毛片网站| 欧美最新免费一区二区三区| 天美传媒精品一区二区| 波多野结衣巨乳人妻| 色播亚洲综合网| 成人永久免费在线观看视频| 琪琪午夜伦伦电影理论片6080| 国产在线精品亚洲第一网站| 一区福利在线观看| 免费看av在线观看网站| 最好的美女福利视频网| 国产欧美日韩一区二区精品| 在线观看免费视频日本深夜| 亚洲七黄色美女视频| 91狼人影院| 丰满乱子伦码专区| 免费人成在线观看视频色| 国产男人的电影天堂91| 午夜久久久久精精品| 久久这里只有精品中国| 极品教师在线视频| 亚洲成人精品中文字幕电影| 非洲黑人性xxxx精品又粗又长| 22中文网久久字幕| 国产欧美日韩一区二区精品| 一区二区三区激情视频| 亚洲欧美清纯卡通| 日本黄大片高清| 午夜老司机福利剧场| 久久久久久久久久成人| 国产精品自产拍在线观看55亚洲| 成人精品一区二区免费| 我的女老师完整版在线观看| 直男gayav资源| 久久久成人免费电影| 免费黄网站久久成人精品| 久久精品久久久久久噜噜老黄 | 国产真实乱freesex| 麻豆久久精品国产亚洲av| 99久久精品国产国产毛片| 国产亚洲精品久久久久久毛片| 久久久久久久精品吃奶| 婷婷亚洲欧美| 亚洲精品久久国产高清桃花| 三级男女做爰猛烈吃奶摸视频| 性插视频无遮挡在线免费观看| 又黄又爽又免费观看的视频| 亚洲黑人精品在线| h日本视频在线播放| 免费搜索国产男女视频| 精品福利观看| 成人无遮挡网站| 亚洲第一区二区三区不卡|