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

    靜電放電算法-投影尋蹤融合模型及其在水質(zhì)綜合評價中的應(yīng)用

    2019-02-12 14:10:00李祥蓉崔東文
    水資源與水工程學(xué)報 2019年6期
    關(guān)鍵詞:盤龍電子設(shè)備功能區(qū)

    李祥蓉, 崔東文

    (1.云南省水文水資源局文山分局 云南 文山 663000; 2.云南省文山州水務(wù)局, 云南 文山 663000)

    1 研究背景

    科學(xué)、客觀評價水功能區(qū)水質(zhì)狀況對于區(qū)域河長制考核評價、落實最嚴(yán)格水資源管理制度具有重要意義。開展水質(zhì)綜合評價方法研究一直是水資源評價領(lǐng)域的熱點和難點,目前灰色關(guān)聯(lián)法[1]、TOPSIS法[2]、神經(jīng)網(wǎng)絡(luò)法[3]、模糊法[4]、支持向量機(jī)法[5]、集對分析法[6]等方法已在水質(zhì)綜合評價中得到應(yīng)用。投影尋蹤(projection pursuit,PP)技術(shù)是將高維數(shù)據(jù)投影到低維空間,并在低維空間進(jìn)行數(shù)據(jù)分析來獲取最優(yōu)空間投影向量,目前已在水質(zhì)綜合評價[7-9]中得到廣泛應(yīng)用。研究表明,PP在實際應(yīng)用過程中,科學(xué)確定其最優(yōu)投影向量是客觀評價的基礎(chǔ)和關(guān)鍵。傳統(tǒng)遺傳算法(GA)、粒子群優(yōu)化(PSO)算法等普遍存在收斂速度慢、易陷入局部最優(yōu)等問題,很難搜索得到真正意義上的最優(yōu)投影向量[10]。目前一些新型群體智能算法嘗試用于投影尋蹤最優(yōu)空間投影向量的優(yōu)化,并取得了較好的優(yōu)化效果,如足球聯(lián)賽競爭(SLC)算法[11]、雞群優(yōu)化(CSO)算法[12]、蛾群算法(MSA)[13]、風(fēng)力驅(qū)動優(yōu)化(WDO)[14]算法等。

    為能科學(xué)、客觀評價水功能區(qū)水質(zhì)狀況,進(jìn)一步拓展新型群體智能算法與投影尋蹤融合模型在水質(zhì)綜合評價中的應(yīng)用范疇,本文內(nèi)容安排如下:(1)介紹一種新型群體智能算法——靜電放電算法(Electrostatic Discharge Algorithm,ESDA),選取6個標(biāo)準(zhǔn)測試函數(shù)對ESDA進(jìn)行仿真測試,并與當(dāng)前具有較好尋優(yōu)能力的新型群體智能算法——教學(xué)優(yōu)化(TLBO)算法、灰狼優(yōu)化(GWO)算法、鯨魚優(yōu)化算法(WOA)的仿真結(jié)果進(jìn)行比較。(2)以文山州盤龍河干流4個斷面2015-2017年共30個月水功能區(qū)水質(zhì)綜合評價分析為例,采用實測水質(zhì)數(shù)據(jù)和所選評價指標(biāo)在《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》(GB 3838-2002)中Ⅰ~Ⅴ類分級閾值構(gòu)建水質(zhì)綜合評價投影目標(biāo)函數(shù),利用ESDA求解該目標(biāo)函數(shù)最佳投影方向,構(gòu)建ESDA-PP模型對實例水功能區(qū)水質(zhì)進(jìn)行綜合評價及分析。

    2 靜電放電算法及仿真驗證

    2.1 靜電放電算法

    靜電放電算法(Electrostatic Discharge Algorithm,ESDA)是Houssem[15]于2019年受靜電放電事件啟發(fā)提出的一種新型元啟發(fā)式優(yōu)化算法。ESDA假設(shè)電子設(shè)備適應(yīng)度與搜索空間位置有關(guān),利用電子設(shè)備個體間直接或間接靜電放電現(xiàn)象,通過具有最低適應(yīng)度值電子設(shè)備向高適應(yīng)度值電子設(shè)備移動來改變其位置以獲得最佳空間位置,即問題最優(yōu)解[15]。與傳統(tǒng)群智能算法相比,ESDA設(shè)置參數(shù)少、收斂速度快、尋優(yōu)精度高,是一種具有較好競爭力的群體智能算法。

    ESDA數(shù)學(xué)描述簡述如下[15]:

    (1)初始化階段。在搜索空間中隨機(jī)生成群體大小n,即電子設(shè)備數(shù)量。每個電子設(shè)備由類似設(shè)計變量的不同組件組成。電子設(shè)備適應(yīng)度與搜索空間中設(shè)備的位置有關(guān),適應(yīng)度值越大,設(shè)備的免疫力就越強(qiáng)。每個設(shè)備都有一個計數(shù)器,用于計算該設(shè)備作為損害設(shè)備的次數(shù)。

    在搜索空間中,利用公式(1)隨機(jī)生成n個電子設(shè)備初始解:

    xi=LB+rand(UB-LB)

    (1)

    式中:xi為電子設(shè)備群體中第i個電子設(shè)備的空間位置(i=1,2,…,n);UB、LB分別為搜索空間的上、下界。

    (2)在每次迭代過程中,從電子設(shè)備群體中隨機(jī)選擇3個電子設(shè)備的適應(yīng)度值并按降序排序(即最佳適應(yīng)度值個體排序為第1)。生成隨機(jī)數(shù)r1,如果r1>0.5,則靜電放電現(xiàn)象僅在兩個電子設(shè)備間發(fā)生;否則,涉及所有3個電子設(shè)備。每當(dāng)設(shè)備受損害或電子設(shè)備受到靜電放電時,該設(shè)備計數(shù)器將增加1。

    ①若r1>0.5,假設(shè)具有最低適應(yīng)度值的電子設(shè)備2朝著具有最高適應(yīng)度值的電子設(shè)備1移動。當(dāng)電子設(shè)備2移動到電子設(shè)備1附近時,假設(shè)在兩個個體之間發(fā)生靜電放電現(xiàn)象,且電子設(shè)備2受到損害,這相當(dāng)于直接靜電放電事件。電子設(shè)備2新位置更新公式如下:

    x2new=x2+2β1(x1-x2)

    (2)

    式中:x2new為電子設(shè)備2新空間位置;x1、x2分別為電子設(shè)備1、2當(dāng)前空間位置;β1為服從正態(tài)分布的隨機(jī)數(shù),其均值μ= 0.7,標(biāo)準(zhǔn)差σ= 0.2。

    ②若r1≤0.5,靜電放電涉及所有隨機(jī)選擇的3個電子設(shè)備,假設(shè)第3個電子設(shè)備正朝著另外兩個電子設(shè)備移動。當(dāng)電子設(shè)備3移動到電子設(shè)備1和電子設(shè)備2附近時,認(rèn)為發(fā)生靜電放電現(xiàn)象,且電子設(shè)備3受到損害,這相當(dāng)于間接靜電放電事件。電子設(shè)備3新位置更新公式如下:

    x3new=x3+2β2(x1-x3)+2β3(x2-x3)

    (3)

    式中:x3new為電子設(shè)備3新空間位置;x1、x2、x3分別為電子設(shè)備1、2、3當(dāng)前空間位置;β2和β3均為服從正態(tài)分布的隨機(jī)數(shù),其均值μ= 0.7,標(biāo)準(zhǔn)差σ= 0.2。

    (3)假設(shè)電子設(shè)備受靜電放電超過3次,則認(rèn)為該電子設(shè)備已被損壞,必須更新,即在搜索空間中隨機(jī)生成新電子設(shè)備;否則,如果電子設(shè)備受靜電放電事件影響小于或等于3次,則產(chǎn)生隨機(jī)數(shù)r2,如果r2<0.2,則假定該電子設(shè)備的一部分已損壞并且必須改變,否則電子設(shè)備組件是安全的。

    (4)將新增電子設(shè)備補(bǔ)充到原電子設(shè)備群體中并保存,再按照電子設(shè)備個體適應(yīng)度值的大小進(jìn)行降序排序,選擇前n個電子設(shè)備作為下一次迭代群體總量。

    2.2 仿真驗證

    為驗證ESDA尋優(yōu)能力,利用ESDA對文獻(xiàn)[16]中Sphere、Griewank、Quadric、Rastrigin、Schwefel 2.22和Ackley 6個典型測試函數(shù)進(jìn)行仿真驗證,并與當(dāng)前具有較好尋優(yōu)能力的新型群體智能算法TLBO、GWO和WOA的仿真結(jié)果進(jìn)行比較。其中函數(shù)Sphere、Quadric、Schwefel 2.22為單峰函數(shù),用于測試算法的尋優(yōu)精度;函數(shù)Griewank、Rastrigin、Ackley為多峰函數(shù),用于測試算法的全局搜索能力。4種算法對6個基準(zhǔn)測試函數(shù)重復(fù)進(jìn)行20次尋優(yōu)計算,并從平均值、標(biāo)準(zhǔn)差2個方面進(jìn)行評估,見表1。

    驗證實驗參數(shù)設(shè)置如下:ESDA、TLBO、GWO和WOA 4種算法最大迭代次數(shù)T=100,群體規(guī)模N=100。其中,TLBO算法參數(shù)TF為1~10之間隨機(jī)整數(shù);WOA對數(shù)螺旋形狀常數(shù)b=2;其他參數(shù)采用各算法默認(rèn)值。

    表1 6個典型測試函數(shù)優(yōu)化結(jié)果對比

    從表1可以看出,對于Sphere、Griewank、Quadric、Rastrigin、Schwefel 2.22函數(shù),ESDA 20次尋優(yōu)均獲得了理論最優(yōu)值0,尋優(yōu)效果遠(yuǎn)優(yōu)于TLBO、GWO和WOA 3種算法;對于多峰函數(shù)Ackley,ESDA 20次尋優(yōu)均獲得相對理論最優(yōu)值8.88×10-16,尋優(yōu)精度高于TLBO、GWO、WOA算法6個數(shù)量級以上??梢?,無論是單峰函數(shù)還是多峰函數(shù),ESDA均表現(xiàn)出較為理想的尋優(yōu)效果,其具有出色的搜索平衡能力、極值尋優(yōu)能力和跳出局部極值能力。

    3 ESDA-PP評價模型

    3.1 投影尋蹤(PP)技術(shù)

    (1)利用公式(4)、(5)分別對正向、負(fù)向指標(biāo)進(jìn)行處理。

    x(i,j)=(x(i,j)-xmin(j))/(xmax(j)-xmin(j))

    (4)

    x(i,j)=(xmax(j)-x(i,j))/(xmax(j)-xmin(j))

    (5)

    式中:x(i,j)為指標(biāo)歸一化序列;xmax(j)、xmin(j)分別為第j個指標(biāo)值上、下限值。

    (2)構(gòu)造投影值z(i)指標(biāo)函數(shù):

    (6)

    式中:a為單位長度向量。

    (3)將搜尋最優(yōu)投影向量問題轉(zhuǎn)化為單目標(biāo)非線性最優(yōu)求解問題,即:

    (7)

    式中:Sz為投影值z(i)的標(biāo)準(zhǔn)差;Dz為投影值z(i)的局部密度;Sz、Dz表達(dá)式參見文獻(xiàn)[12];Q(a)為待優(yōu)化目標(biāo)函數(shù)。

    3.2 ESDA-PP模型評價實現(xiàn)步驟

    ESDA-PP水質(zhì)綜合評價模型的實現(xiàn)步驟歸納如下:

    從資產(chǎn)流轉(zhuǎn)順暢的角度來看,產(chǎn)權(quán)的完整程度主要體現(xiàn)在集體資產(chǎn)是否能夠自由進(jìn)出市場。與此同時,這也能夠更有效地實現(xiàn)要素的優(yōu)化與組合。當(dāng)集體資產(chǎn)順暢流轉(zhuǎn)時,集體資產(chǎn)的保值與增值效果才能達(dá)到最好。隨著產(chǎn)權(quán)的逐步明確,農(nóng)村經(jīng)濟(jì)產(chǎn)權(quán)的改革也得到進(jìn)一步完善,構(gòu)建與之相匹配的經(jīng)濟(jì)產(chǎn)權(quán)交易市場,相關(guān)的交易規(guī)則與程序也能夠逐步得到完善,最終實現(xiàn)公平交易[1]。

    Step1選取水質(zhì)綜合評價指標(biāo)、序列長度和各評價指標(biāo)Ⅰ~Ⅴ類閾值,對樣本進(jìn)行一致性處理,利用處理后的樣本構(gòu)建水質(zhì)評價待優(yōu)化目標(biāo)函數(shù)Q(a),利用ESDA對Q(a)投影向量進(jìn)行優(yōu)化求解。

    Step2設(shè)置ESDA群體規(guī)模N、最大迭代次數(shù) 和算法終止條件;利用公式(1)生成電子設(shè)備初始空間位置,令當(dāng)前迭代次數(shù)t=0。

    Step3隨機(jī)選擇3個電子設(shè)備,計算待優(yōu)化問題目標(biāo)函數(shù)適應(yīng)度值并按降序進(jìn)行排序。

    Step4 在(0,1)之間生成隨機(jī)數(shù)r1,如果r1>0.5,則靜電放電現(xiàn)象僅在兩個電子設(shè)備間發(fā)生,利用公式(2)更新電子設(shè)備2位置;若r1≤0.5,涉及所有3個電子設(shè)備,利用公式(3)更新電子設(shè)備3位置。

    Step5 檢查電子設(shè)備空間位置是否超出搜索范圍,若電子設(shè)備空間位置位于搜索空間之外,則利用公式(1)將其強(qiáng)制帶回搜索空間。

    Step6 依據(jù)電子設(shè)備受靜電放電次數(shù)更新電子設(shè)備群體,并將新增電子設(shè)備補(bǔ)充到原電子設(shè)備種群中并保存。按電子設(shè)備個體適應(yīng)度值的大小進(jìn)降序排序,選擇前N個電子設(shè)備作為下一次迭代群體規(guī)模。

    Step7 令t=t+1,判斷算法是否滿足終止條件,若是,輸出最優(yōu)解xbest,尋優(yōu)過程結(jié)束;否則重復(fù)Step3~Step7。

    Step8輸出ESDA最優(yōu)解xbest,即最佳投影向量a。利用a計算各評價對象投影值z及各水質(zhì)類別分級投影值z′,利用水質(zhì)類別分級投影值z′對實例水功能區(qū)水質(zhì)進(jìn)行綜合評價及分析。

    4 應(yīng)用實例

    4.1 研究區(qū)概況與數(shù)據(jù)來源

    盤龍河屬紅河流域,瀘江水系,發(fā)源于紅河州蒙自縣三道灣,流經(jīng)硯山、文山等6縣市,于天保船頭注入越南。盤龍河中越國界以上流域面積6 497 km2,河長247 km,落差1 803 m,平均坡降7.24‰,多年平均水資源量26.7×108m3。依據(jù)《文山州地表水功能區(qū)劃(2015修訂)》,盤龍河干流由上至下依次劃分為盤龍河硯山-文山農(nóng)業(yè)工業(yè)用水區(qū),水質(zhì)代表斷面天生橋2020和2030年(規(guī)劃年)的水質(zhì)保護(hù)目標(biāo)為Ⅱ類;盤龍河文山飲用水源區(qū)的水質(zhì)代表斷面為灰土寨,規(guī)劃年水質(zhì)保護(hù)目標(biāo)為Ⅱ類;盤龍河文山景觀農(nóng)業(yè)用水區(qū)的水質(zhì)代表斷面為花橋,規(guī)劃年水質(zhì)保護(hù)目標(biāo)為Ⅲ類;盤龍河文山-麻栗坡保留區(qū)的水質(zhì)代表斷面為天保大橋,規(guī)劃年水質(zhì)保護(hù)目標(biāo)為Ⅲ類。開展盤龍河干流水功能區(qū)水質(zhì)綜合評價對于區(qū)域河長制考核評價、落實最嚴(yán)格水資源管理制度具有重要意義。

    本文選取對水功能區(qū)水質(zhì)影響較大的DO、NH3—N、CODMn、BOD5和TP作為水質(zhì)綜合評價因子,通過構(gòu)建ESDA-PP模型對盤龍河干流天生橋、灰土寨、花橋、天保大橋4個斷面2015-2017年水功能區(qū)水質(zhì)進(jìn)行綜合評價,其中2015年為雙月監(jiān)測,2016、2017年為逐月監(jiān)測。限于篇幅,水質(zhì)監(jiān)測數(shù)據(jù)從略。

    4.2 模型求解

    依據(jù)ESDA-PP水質(zhì)綜合評價模型實現(xiàn)步驟,利用天生橋、灰土寨、花橋、天保大橋4個斷面2015-2017年共30個月斷面水質(zhì)實測數(shù)據(jù)和所選5個評價指標(biāo)在《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》(GB3838-2002)中Ⅰ~Ⅴ類分級閾值構(gòu)建水質(zhì)綜合評價投影目標(biāo)函數(shù)Q(a),利用ESDA求解Q(a)最佳投影方向量a。ESDA參數(shù)設(shè)置同上。經(jīng)求解,20次優(yōu)化目標(biāo)Q(a)值均為18 934.447,標(biāo)準(zhǔn)差均為0;20次優(yōu)化投影向量均相同,即a=(-0.0754 0.5530 0.4440 0.5530 0.4307)。該實例優(yōu)化結(jié)果再次驗證了ESDA具有較好的尋優(yōu)精度和極值尋優(yōu)能力。

    4.3 結(jié)果與分析

    (1)構(gòu)造水質(zhì)評價等級標(biāo)準(zhǔn)。利用上述優(yōu)化獲得的最佳投影向量a和經(jīng)一致性處理后的Ⅰ~Ⅴ類分級閾值,基于公式(6)計算得到Ⅰ~Ⅴ類分級標(biāo)準(zhǔn)投影值z′,z′即為水質(zhì)綜合評價分級標(biāo)準(zhǔn),即:Ⅰ類<0.1801、Ⅱ類 ∈[0.1801,0.3303)、Ⅲ類 ∈[0.3303,0.5795)、Ⅳ類 ∈[0.5795,0.9411)、Ⅴ類 ∈[0.9411,1.4574)和劣Ⅴ類≥1.4574。

    (2)綜合評價。同理,依據(jù)上述優(yōu)化獲得的最佳投影向量a和經(jīng)一致性處理后的天生橋、灰土寨、花橋、天保大橋4個斷面2015-2017年共30個月斷面水質(zhì)實測數(shù)據(jù),基于公式(6)計算得到各斷面各月份綜合投影值z,利用水質(zhì)綜合評價分級標(biāo)準(zhǔn)z′進(jìn)行評價與分析,見表2及圖1~2。

    表2 2015-2017年盤龍河4個斷面水質(zhì)綜合投影值及評價結(jié)果

    圖12015-2017年30個月各斷面水質(zhì)評價結(jié)果 圖22015-2017年30個月各斷面枯、豐水期水質(zhì)評價結(jié)果

    從表2及圖1、2可以得出以下結(jié)論:

    (1)ESDA-PP模型對盤龍河4個斷面評價結(jié)果與單因子評價結(jié)果存在差異,總體評價結(jié)果要優(yōu)于單因子評價結(jié)果Ⅰ~Ⅱ類。從評價數(shù)據(jù)來看,影響盤龍河水質(zhì)主要因子為NH3—N,單因子評價法僅依據(jù)NH3—N單項指標(biāo)評價各斷面水質(zhì)類別,顯然不夠科學(xué)、客觀。

    (2)依據(jù)ESDA-PP模型評價結(jié)果,并按頻次法大于等于80%即為達(dá)標(biāo)來看,天生橋、灰土寨、花橋和天保大橋斷面在研究期的30個月內(nèi)水功能區(qū)水質(zhì)達(dá)標(biāo)率分別為93.3%、83.3%、83.3%和100%,能滿足或基本滿足水功能區(qū)水質(zhì)保護(hù)目標(biāo)要求。其中天生橋斷面水質(zhì)處于Ⅰ~Ⅲ類之間,水質(zhì)保護(hù)目標(biāo)Ⅱ類,30個月水功能區(qū)水質(zhì)達(dá)標(biāo)率為93.3%,水功能區(qū)水質(zhì)能滿足Ⅱ類水質(zhì)保護(hù)目標(biāo)要求;灰土寨斷面水質(zhì)處于Ⅰ~Ⅲ類之間,水質(zhì)保護(hù)目標(biāo)Ⅱ類,30個月水功能區(qū)水質(zhì)達(dá)標(biāo)率為83.3%,水功能區(qū)水質(zhì)基本能滿足Ⅱ類水質(zhì)保護(hù)目標(biāo)要求;花橋斷面總體水質(zhì)處于Ⅱ~Ⅴ類,水質(zhì)保護(hù)目標(biāo)Ⅲ類,30個月水功能區(qū)水質(zhì)達(dá)標(biāo)率為83.3%,水功能區(qū)水質(zhì)基本能滿足Ⅱ類水質(zhì)保護(hù)目標(biāo)要求;天保大橋斷面水質(zhì)處于Ⅰ~Ⅱ類,水質(zhì)保護(hù)目標(biāo)Ⅲ類, 30個月水功能區(qū)水質(zhì)達(dá)標(biāo)率為100%,水功能區(qū)水質(zhì)完全能滿足Ⅱ類水質(zhì)保護(hù)目標(biāo)要求。若按單因子評價法進(jìn)行水功能區(qū)水質(zhì)達(dá)標(biāo)評價,則天生橋、灰土寨、花橋、天保大橋30個月水功能區(qū)水質(zhì)達(dá)標(biāo)率分別為46.7%、56.7%、33.3%、100%,除天保大橋外,其余3個斷面均不能滿足水功能區(qū)水質(zhì)保護(hù)目標(biāo)要求。

    (3)天生橋、灰土寨、花橋斷面枯水期水質(zhì)總體上要優(yōu)于豐水期,說明農(nóng)田和城市面源污染是影響盤龍河水質(zhì)的主導(dǎo)因素。天保大橋位于盤龍河下游,汛期來水較大,豐水期水質(zhì)總體上要優(yōu)于枯水期,說明盤龍河下游段水質(zhì)總體尚好的主要原因是巨大的水量和流量,也間接說明點源污染是影響盤龍河下游段水質(zhì)的主要因素。

    5 結(jié) 論

    本文介紹一種新型群體智能算法——靜電放電算法(ESDA),通過對ESDA進(jìn)行仿真驗證和構(gòu)建ESDA-PP評價模型對盤龍河4個斷面2015-2017年30個月水功能區(qū)水質(zhì)進(jìn)行實例綜合評價,得出以下結(jié)論:

    (1)選取6個標(biāo)準(zhǔn)測試函數(shù)對ESDA進(jìn)行仿真驗證,并與TLBO、GWO、WOA仿真結(jié)果進(jìn)行對比,結(jié)果表明ESDA尋優(yōu)精度遠(yuǎn)優(yōu)于TLBO、GWO、WOA,具有出色的尋優(yōu)精度和全局搜索能力,是一種具有較好競爭力的群體智能算法。

    (2)ESDA-PP模型對實例評價結(jié)果顯示,盤龍河天生橋、灰土寨、花橋和天保大橋斷面2015-2017年30個月水功能區(qū)水質(zhì)達(dá)標(biāo)率分別為93.3%、83.3%、83.3%和100%,能滿足或基本滿足水功能區(qū)水質(zhì)保護(hù)目標(biāo)要求。天生橋、灰土寨、花橋斷面枯水期水質(zhì)總體上要優(yōu)于豐水期,農(nóng)田和城市面源污染是影響盤龍河水質(zhì)的決定因素;天保大橋豐水期水質(zhì)總體上要優(yōu)于枯水期,點源污染是制約盤龍河下游段水質(zhì)的主要因素。

    猜你喜歡
    盤龍電子設(shè)備功能區(qū)
    有機(jī)硅灌封材料對計控電子設(shè)備的保護(hù)應(yīng)用
    山東冶金(2022年4期)2022-09-14 09:00:20
    某型機(jī)載敏感電子設(shè)備電磁干擾抑制技術(shù)研究
    云南省昆明市盤龍區(qū)盤龍小學(xué)五(5)螢火蟲中隊
    少先隊活動(2021年4期)2021-07-23 01:46:36
    武漢盤龍廣場銀興影城
    燒結(jié)混合機(jī)各功能區(qū)的研究與優(yōu)化
    山東冶金(2019年3期)2019-07-10 00:53:54
    盤龍糕
    幼兒100(2018年29期)2018-10-30 01:38:00
    衡水市不同功能區(qū)土壤營養(yǎng)元素的比較
    探析現(xiàn)代醫(yī)用電子設(shè)備的維修
    電子測試(2017年15期)2017-12-18 07:19:55
    小白兔采蘑菇
    小布老虎(2017年2期)2017-07-31 17:14:37
    水功能區(qū)納污能力計算及污染物總量控制
    国产精品无大码| 在线观看www视频免费| 婷婷色麻豆天堂久久| 国产成人系列免费观看| 秋霞在线观看毛片| 亚洲四区av| 美女福利国产在线| 欧美日韩一区二区视频在线观看视频在线| 精品少妇一区二区三区视频日本电影 | 欧美精品av麻豆av| 91成人精品电影| 亚洲欧洲日产国产| 视频区图区小说| 国产在线视频一区二区| 亚洲少妇的诱惑av| 欧美 日韩 精品 国产| 777久久人妻少妇嫩草av网站| 成人影院久久| 18禁国产床啪视频网站| 妹子高潮喷水视频| 久久ye,这里只有精品| 亚洲国产av影院在线观看| 亚洲国产精品一区二区三区在线| 熟妇人妻不卡中文字幕| 一级毛片我不卡| 国产视频首页在线观看| 午夜av观看不卡| 国产欧美亚洲国产| 少妇被粗大的猛进出69影院| 亚洲三区欧美一区| 色婷婷久久久亚洲欧美| 亚洲国产欧美日韩在线播放| 80岁老熟妇乱子伦牲交| 男女午夜视频在线观看| 丝袜在线中文字幕| 狂野欧美激情性bbbbbb| 黄色毛片三级朝国网站| 2018国产大陆天天弄谢| 日本vs欧美在线观看视频| 亚洲av中文av极速乱| 国产精品av久久久久免费| 1024香蕉在线观看| 纯流量卡能插随身wifi吗| 久久久精品区二区三区| 色婷婷久久久亚洲欧美| 夜夜骑夜夜射夜夜干| 成人影院久久| 老司机深夜福利视频在线观看 | 国产精品 欧美亚洲| 极品人妻少妇av视频| 免费不卡黄色视频| 日韩一卡2卡3卡4卡2021年| 久久久久国产精品人妻一区二区| 岛国毛片在线播放| 日韩精品免费视频一区二区三区| 日日摸夜夜添夜夜爱| 欧美另类一区| 免费观看a级毛片全部| 欧美日韩成人在线一区二区| 99国产综合亚洲精品| 亚洲av日韩精品久久久久久密 | 国产av码专区亚洲av| 欧美97在线视频| 国产在线一区二区三区精| 纯流量卡能插随身wifi吗| 在线亚洲精品国产二区图片欧美| tube8黄色片| 国产成人午夜福利电影在线观看| 色94色欧美一区二区| 久久99一区二区三区| 国产极品粉嫩免费观看在线| 久久久久人妻精品一区果冻| 国产一区二区 视频在线| 超碰成人久久| 久久精品国产a三级三级三级| 男人添女人高潮全过程视频| 精品国产一区二区三区久久久樱花| 亚洲av日韩在线播放| 欧美日韩福利视频一区二区| 国产精品无大码| 成人手机av| 亚洲天堂av无毛| 国产熟女午夜一区二区三区| 久久精品国产综合久久久| 人人妻,人人澡人人爽秒播 | 18禁观看日本| 亚洲精品一区蜜桃| 亚洲av在线观看美女高潮| 欧美精品亚洲一区二区| 久久热在线av| 国产精品一区二区在线不卡| 伦理电影免费视频| √禁漫天堂资源中文www| 欧美中文综合在线视频| 国语对白做爰xxxⅹ性视频网站| 日本欧美国产在线视频| 极品少妇高潮喷水抽搐| 国产精品一区二区在线观看99| 男女之事视频高清在线观看 | 亚洲免费av在线视频| 国产精品三级大全| 午夜91福利影院| 啦啦啦啦在线视频资源| 午夜福利一区二区在线看| 99国产综合亚洲精品| 大香蕉久久成人网| 极品人妻少妇av视频| 天堂8中文在线网| 国语对白做爰xxxⅹ性视频网站| 精品国产超薄肉色丝袜足j| 亚洲专区中文字幕在线 | 免费观看人在逋| 久久性视频一级片| 国产亚洲午夜精品一区二区久久| 91精品国产国语对白视频| 秋霞在线观看毛片| 久热这里只有精品99| 亚洲av日韩精品久久久久久密 | 亚洲欧美日韩另类电影网站| 自线自在国产av| 亚洲成人国产一区在线观看 | 热99国产精品久久久久久7| 久久精品久久久久久久性| 国产爽快片一区二区三区| 久久精品aⅴ一区二区三区四区| kizo精华| 亚洲第一区二区三区不卡| 日本午夜av视频| √禁漫天堂资源中文www| 亚洲综合精品二区| 久久久久网色| 天天躁狠狠躁夜夜躁狠狠躁| 欧美变态另类bdsm刘玥| 天堂中文最新版在线下载| 日本vs欧美在线观看视频| 久久青草综合色| 国产精品一区二区精品视频观看| 男女午夜视频在线观看| 啦啦啦 在线观看视频| 在线天堂中文资源库| 满18在线观看网站| 极品人妻少妇av视频| 久久人人爽av亚洲精品天堂| 99九九在线精品视频| 亚洲国产日韩一区二区| 永久免费av网站大全| 一级,二级,三级黄色视频| av网站免费在线观看视频| 国产av精品麻豆| 午夜福利网站1000一区二区三区| 黄片无遮挡物在线观看| 久久午夜综合久久蜜桃| 一级a爱视频在线免费观看| 在线观看一区二区三区激情| 亚洲av日韩在线播放| 国产国语露脸激情在线看| 观看美女的网站| 亚洲欧美中文字幕日韩二区| 国产又爽黄色视频| 国产精品.久久久| 国产精品成人在线| 国产人伦9x9x在线观看| 国产精品一区二区精品视频观看| 日本91视频免费播放| 19禁男女啪啪无遮挡网站| 欧美精品人与动牲交sv欧美| 视频区图区小说| 最近中文字幕高清免费大全6| 日韩免费高清中文字幕av| 精品国产超薄肉色丝袜足j| 日韩电影二区| 中文字幕人妻丝袜制服| 别揉我奶头~嗯~啊~动态视频 | av又黄又爽大尺度在线免费看| 老汉色av国产亚洲站长工具| 成人毛片60女人毛片免费| 男女高潮啪啪啪动态图| 高清av免费在线| 久久久久视频综合| 日韩制服骚丝袜av| bbb黄色大片| 欧美日韩一级在线毛片| 黄片无遮挡物在线观看| 大香蕉久久成人网| 超色免费av| 丝袜人妻中文字幕| 久久天躁狠狠躁夜夜2o2o | 激情视频va一区二区三区| 韩国av在线不卡| 午夜激情av网站| 亚洲欧美精品综合一区二区三区| 国产精品久久久久成人av| 一本久久精品| 夜夜骑夜夜射夜夜干| 女性被躁到高潮视频| 国产精品熟女久久久久浪| 99久国产av精品国产电影| 一边摸一边抽搐一进一出视频| 国产一区亚洲一区在线观看| av片东京热男人的天堂| 色吧在线观看| 高清av免费在线| 国产人伦9x9x在线观看| 我的亚洲天堂| 一区二区三区四区激情视频| 丁香六月欧美| 欧美乱码精品一区二区三区| 精品久久久精品久久久| 亚洲自偷自拍图片 自拍| 在线 av 中文字幕| 丝袜脚勾引网站| 日韩中文字幕欧美一区二区 | 国产 精品1| 侵犯人妻中文字幕一二三四区| 久久精品亚洲熟妇少妇任你| 最近2019中文字幕mv第一页| 国产日韩欧美视频二区| 亚洲精品中文字幕在线视频| 久久久久久久大尺度免费视频| 国产亚洲一区二区精品| 婷婷色av中文字幕| 三上悠亚av全集在线观看| 久久精品国产a三级三级三级| 日本91视频免费播放| 精品一品国产午夜福利视频| 少妇被粗大的猛进出69影院| 亚洲成av片中文字幕在线观看| 99国产精品免费福利视频| 岛国毛片在线播放| e午夜精品久久久久久久| 色婷婷av一区二区三区视频| 久久久久视频综合| av国产精品久久久久影院| 青草久久国产| 99精品久久久久人妻精品| e午夜精品久久久久久久| 一边亲一边摸免费视频| 国产片内射在线| 国产又色又爽无遮挡免| 精品视频人人做人人爽| 人成视频在线观看免费观看| 亚洲欧洲精品一区二区精品久久久 | 七月丁香在线播放| 一本—道久久a久久精品蜜桃钙片| 美女福利国产在线| 日韩电影二区| 天天躁夜夜躁狠狠久久av| 在线观看国产h片| 韩国精品一区二区三区| 精品酒店卫生间| 欧美日韩福利视频一区二区| 啦啦啦视频在线资源免费观看| 亚洲精品国产区一区二| 久久精品国产亚洲av涩爱| 波多野结衣一区麻豆| 久久精品久久精品一区二区三区| 国产精品久久久久久久久免| 免费女性裸体啪啪无遮挡网站| 久久久久精品国产欧美久久久 | 亚洲美女黄色视频免费看| 精品国产国语对白av| 国产成人午夜福利电影在线观看| 男人操女人黄网站| 国产精品熟女久久久久浪| 欧美人与善性xxx| 在线天堂最新版资源| 欧美日韩福利视频一区二区| 老司机在亚洲福利影院| 深夜精品福利| 国产国语露脸激情在线看| netflix在线观看网站| 十八禁高潮呻吟视频| 天天躁夜夜躁狠狠久久av| 卡戴珊不雅视频在线播放| 亚洲美女搞黄在线观看| 国产一区二区三区综合在线观看| 一级片免费观看大全| 午夜免费鲁丝| 国产在线一区二区三区精| 亚洲国产精品一区二区三区在线| 亚洲av电影在线观看一区二区三区| 久久人人爽人人片av| 日本黄色日本黄色录像| 日韩中文字幕视频在线看片| 国产精品一区二区精品视频观看| 美女中出高潮动态图| 午夜av观看不卡| 国精品久久久久久国模美| 精品亚洲成国产av| 欧美精品一区二区免费开放| 一区二区三区乱码不卡18| 午夜免费男女啪啪视频观看| 大码成人一级视频| 日韩一区二区三区影片| 免费观看人在逋| √禁漫天堂资源中文www| 久久久精品区二区三区| 国产av码专区亚洲av| 一级毛片电影观看| 亚洲少妇的诱惑av| 香蕉国产在线看| 亚洲国产精品成人久久小说| 欧美激情极品国产一区二区三区| 少妇人妻 视频| 在现免费观看毛片| 欧美激情 高清一区二区三区| 在线观看免费日韩欧美大片| 日日摸夜夜添夜夜爱| 国产精品麻豆人妻色哟哟久久| 男女免费视频国产| 色视频在线一区二区三区| 欧美变态另类bdsm刘玥| 多毛熟女@视频| 日韩视频在线欧美| 成年女人毛片免费观看观看9 | 乱人伦中国视频| 如日韩欧美国产精品一区二区三区| 美女国产高潮福利片在线看| 欧美xxⅹ黑人| www日本在线高清视频| 国产一区二区激情短视频 | 女性生殖器流出的白浆| 操美女的视频在线观看| 看免费av毛片| 欧美成人精品欧美一级黄| 免费观看av网站的网址| 欧美精品av麻豆av| 国产深夜福利视频在线观看| 亚洲欧美精品自产自拍| 91精品国产国语对白视频| 国产精品久久久人人做人人爽| 操出白浆在线播放| 国产 精品1| 天天影视国产精品| av网站免费在线观看视频| 欧美日韩亚洲国产一区二区在线观看 | 国产精品久久久av美女十八| 亚洲欧美精品自产自拍| 国产男人的电影天堂91| 超碰成人久久| 波多野结衣av一区二区av| 视频在线观看一区二区三区| 亚洲欧美一区二区三区久久| 高清视频免费观看一区二区| 中文字幕色久视频| 夜夜骑夜夜射夜夜干| 日本91视频免费播放| 精品国产乱码久久久久久男人| 少妇被粗大猛烈的视频| 国精品久久久久久国模美| 国产精品蜜桃在线观看| 欧美黄色片欧美黄色片| 国产精品偷伦视频观看了| 欧美久久黑人一区二区| 天堂俺去俺来也www色官网| 亚洲第一青青草原| 十八禁人妻一区二区| 国产亚洲最大av| 久热爱精品视频在线9| 精品人妻熟女毛片av久久网站| bbb黄色大片| 王馨瑶露胸无遮挡在线观看| 自线自在国产av| 国产免费现黄频在线看| bbb黄色大片| 亚洲国产最新在线播放| 99久久人妻综合| 欧美激情高清一区二区三区 | 99久久人妻综合| 最近中文字幕2019免费版| 国产精品亚洲av一区麻豆 | 午夜福利在线免费观看网站| 色婷婷av一区二区三区视频| 欧美人与性动交α欧美软件| 免费少妇av软件| 久久久亚洲精品成人影院| 午夜久久久在线观看| 欧美精品亚洲一区二区| 1024香蕉在线观看| 午夜av观看不卡| 精品少妇一区二区三区视频日本电影 | 欧美日韩综合久久久久久| 9色porny在线观看| 在线观看一区二区三区激情| 免费在线观看黄色视频的| 夜夜骑夜夜射夜夜干| 在线观看一区二区三区激情| 伊人亚洲综合成人网| 人人澡人人妻人| 亚洲图色成人| 黄频高清免费视频| www.熟女人妻精品国产| 黄网站色视频无遮挡免费观看| kizo精华| 两性夫妻黄色片| 国产99久久九九免费精品| 亚洲av成人不卡在线观看播放网 | 人妻人人澡人人爽人人| 日韩 欧美 亚洲 中文字幕| 日韩av不卡免费在线播放| 国产亚洲午夜精品一区二区久久| 国产片特级美女逼逼视频| 日韩一区二区视频免费看| 18在线观看网站| 啦啦啦在线免费观看视频4| 国产日韩欧美亚洲二区| 亚洲第一区二区三区不卡| 国产精品欧美亚洲77777| 一边摸一边做爽爽视频免费| 在线天堂中文资源库| 十分钟在线观看高清视频www| 亚洲自偷自拍图片 自拍| 成人午夜精彩视频在线观看| 两性夫妻黄色片| 中文字幕高清在线视频| 巨乳人妻的诱惑在线观看| 亚洲人成电影观看| 91精品国产国语对白视频| 国产色婷婷99| 亚洲av福利一区| 伊人久久国产一区二区| 久久99精品国语久久久| 国产精品久久久久久精品电影小说| 亚洲av电影在线进入| 精品第一国产精品| 青春草视频在线免费观看| 亚洲人成电影观看| www日本在线高清视频| 欧美日韩亚洲综合一区二区三区_| 成人影院久久| 午夜免费观看性视频| 国产97色在线日韩免费| 女人被躁到高潮嗷嗷叫费观| 婷婷色麻豆天堂久久| 老司机亚洲免费影院| 黑丝袜美女国产一区| 97人妻天天添夜夜摸| 日韩 欧美 亚洲 中文字幕| 欧美日韩成人在线一区二区| 亚洲色图综合在线观看| 日韩制服丝袜自拍偷拍| 肉色欧美久久久久久久蜜桃| 一级毛片电影观看| 黄色 视频免费看| 丝袜美腿诱惑在线| 国产精品欧美亚洲77777| 亚洲精品中文字幕在线视频| 亚洲人成77777在线视频| 2018国产大陆天天弄谢| 久久精品久久久久久久性| 亚洲国产日韩一区二区| 最新在线观看一区二区三区 | 一区在线观看完整版| 国语对白做爰xxxⅹ性视频网站| av在线老鸭窝| 日韩av在线免费看完整版不卡| 国产精品一区二区在线不卡| 女人被躁到高潮嗷嗷叫费观| 中文字幕人妻丝袜制服| 久久午夜综合久久蜜桃| www.自偷自拍.com| 纵有疾风起免费观看全集完整版| 欧美另类一区| 亚洲专区中文字幕在线 | 亚洲精品av麻豆狂野| 免费少妇av软件| 两性夫妻黄色片| 咕卡用的链子| 亚洲视频免费观看视频| 操出白浆在线播放| 一边摸一边抽搐一进一出视频| 午夜免费男女啪啪视频观看| 日韩伦理黄色片| 9热在线视频观看99| 老汉色av国产亚洲站长工具| 亚洲成国产人片在线观看| 国产乱来视频区| 精品国产一区二区三区久久久樱花| 大片电影免费在线观看免费| 999精品在线视频| 51午夜福利影视在线观看| 91老司机精品| 国产精品久久久av美女十八| 中文欧美无线码| 我要看黄色一级片免费的| 宅男免费午夜| 色吧在线观看| 啦啦啦啦在线视频资源| 精品国产露脸久久av麻豆| 国产精品一区二区在线不卡| 女人高潮潮喷娇喘18禁视频| 亚洲国产欧美一区二区综合| 在线观看免费高清a一片| 国产淫语在线视频| bbb黄色大片| 免费黄色在线免费观看| 久久久久久久久久久免费av| 久久青草综合色| 精品一品国产午夜福利视频| 久久久久精品久久久久真实原创| 午夜福利影视在线免费观看| 久久av网站| 亚洲av福利一区| 久久久久久人人人人人| 午夜福利网站1000一区二区三区| 丰满乱子伦码专区| 高清视频免费观看一区二区| 久久性视频一级片| av不卡在线播放| www.av在线官网国产| 女人精品久久久久毛片| 免费看av在线观看网站| 国产女主播在线喷水免费视频网站| 久久久久人妻精品一区果冻| av在线app专区| 亚洲熟女精品中文字幕| 亚洲精品日本国产第一区| 久久鲁丝午夜福利片| 多毛熟女@视频| 国产野战对白在线观看| 成人黄色视频免费在线看| 国产片特级美女逼逼视频| 看十八女毛片水多多多| 9色porny在线观看| 久久精品熟女亚洲av麻豆精品| 99香蕉大伊视频| 欧美成人午夜精品| 丝袜美腿诱惑在线| 久久久久久久久免费视频了| 成年动漫av网址| 成人国产av品久久久| 国产老妇伦熟女老妇高清| 日韩熟女老妇一区二区性免费视频| av网站在线播放免费| 高清在线视频一区二区三区| 日本色播在线视频| 97精品久久久久久久久久精品| av女优亚洲男人天堂| 久久久久久久久久久免费av| 亚洲av电影在线观看一区二区三区| videos熟女内射| 看非洲黑人一级黄片| 建设人人有责人人尽责人人享有的| 男女无遮挡免费网站观看| 色婷婷av一区二区三区视频| 飞空精品影院首页| av在线观看视频网站免费| 亚洲国产精品成人久久小说| 精品国产一区二区三区久久久樱花| 黑人欧美特级aaaaaa片| 人妻 亚洲 视频| 各种免费的搞黄视频| 免费黄网站久久成人精品| 欧美日韩成人在线一区二区| 97人妻天天添夜夜摸| 亚洲欧美精品自产自拍| 天天躁日日躁夜夜躁夜夜| videosex国产| 高清不卡的av网站| 欧美激情 高清一区二区三区| 日韩中文字幕视频在线看片| 十八禁高潮呻吟视频| 日韩精品有码人妻一区| 亚洲精品国产区一区二| 午夜日本视频在线| 97人妻天天添夜夜摸| 国产成人精品无人区| 国产精品久久久久成人av| 99国产综合亚洲精品| 中文字幕高清在线视频| 岛国毛片在线播放| 久久99热这里只频精品6学生| 午夜福利一区二区在线看| 欧美日韩亚洲国产一区二区在线观看 | 国产日韩欧美视频二区| 一级,二级,三级黄色视频| 成年人午夜在线观看视频| 欧美日韩一区二区视频在线观看视频在线| 亚洲精品久久午夜乱码| 国产精品二区激情视频| 中文字幕最新亚洲高清| 亚洲人成网站在线观看播放| 国产精品免费大片| 在线观看一区二区三区激情| 国产精品亚洲av一区麻豆 | 一区二区三区精品91| 日韩中文字幕欧美一区二区 | 99国产综合亚洲精品| 咕卡用的链子| 国产高清不卡午夜福利| 精品卡一卡二卡四卡免费| 啦啦啦 在线观看视频| 一区二区三区四区激情视频| 国产在线视频一区二区| 最近中文字幕高清免费大全6| 尾随美女入室| 国产一区二区 视频在线| 亚洲男人天堂网一区| 色视频在线一区二区三区| 精品少妇黑人巨大在线播放| 午夜免费鲁丝| 九草在线视频观看| 国产成人a∨麻豆精品| 国产精品.久久久| 99久久99久久久精品蜜桃| 黑丝袜美女国产一区| 亚洲精品国产av成人精品| 久久人妻熟女aⅴ| 男人爽女人下面视频在线观看| 你懂的网址亚洲精品在线观看| 精品少妇内射三级| 国产av码专区亚洲av| 我要看黄色一级片免费的| 男女免费视频国产| 亚洲欧美精品综合一区二区三区|