• <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ū)納污能力計算及污染物總量控制
    精品人妻偷拍中文字幕| 国产亚洲av片在线观看秒播厂 | 久久久久久国产a免费观看| av在线观看视频网站免费| 91狼人影院| 中文字幕亚洲精品专区| 国产免费视频播放在线视频 | 人妻系列 视频| 一个人观看的视频www高清免费观看| 免费电影在线观看免费观看| 深夜a级毛片| 三级经典国产精品| 国产综合懂色| 久久精品久久久久久久性| 久久久久精品性色| 亚洲国产最新在线播放| 男女国产视频网站| 看非洲黑人一级黄片| 国产有黄有色有爽视频| 国产精品不卡视频一区二区| 国产久久久一区二区三区| a级毛色黄片| 春色校园在线视频观看| 国产免费视频播放在线视频 | 亚洲在线观看片| 国产高清国产精品国产三级 | 亚洲无线观看免费| 观看美女的网站| 国内精品美女久久久久久| 晚上一个人看的免费电影| av国产免费在线观看| 免费观看精品视频网站| 99九九线精品视频在线观看视频| av福利片在线观看| 能在线免费观看的黄片| 欧美bdsm另类| 久久精品熟女亚洲av麻豆精品 | 亚洲欧美一区二区三区国产| 国产亚洲av嫩草精品影院| 最近中文字幕高清免费大全6| 久久6这里有精品| 亚州av有码| 日韩三级伦理在线观看| 日本黄大片高清| 亚洲av日韩在线播放| 三级经典国产精品| 天堂网av新在线| 成人亚洲精品av一区二区| a级毛色黄片| 亚洲国产精品成人久久小说| 欧美日韩国产mv在线观看视频 | 又爽又黄a免费视频| 国产v大片淫在线免费观看| 最近2019中文字幕mv第一页| 日韩亚洲欧美综合| 好男人在线观看高清免费视频| 午夜精品在线福利| 国语对白做爰xxxⅹ性视频网站| 免费观看无遮挡的男女| 一个人免费在线观看电影| 久久精品久久精品一区二区三区| 日日干狠狠操夜夜爽| 亚洲av一区综合| 久久综合国产亚洲精品| 观看美女的网站| 一区二区三区四区激情视频| 亚洲av一区综合| 国产视频首页在线观看| 淫秽高清视频在线观看| 一级片'在线观看视频| 精品午夜福利在线看| 成人国产麻豆网| 看非洲黑人一级黄片| 菩萨蛮人人尽说江南好唐韦庄| 中文欧美无线码| 国产极品天堂在线| 国产三级在线视频| 国产精品一区二区性色av| 网址你懂的国产日韩在线| 亚洲欧美成人精品一区二区| 日韩电影二区| 国产一区有黄有色的免费视频 | 一级爰片在线观看| 亚洲人与动物交配视频| 亚洲欧美中文字幕日韩二区| 69av精品久久久久久| 听说在线观看完整版免费高清| 色吧在线观看| 亚洲伊人久久精品综合| 亚洲怡红院男人天堂| 亚洲国产av新网站| 只有这里有精品99| 22中文网久久字幕| 亚洲欧美精品自产自拍| 成人亚洲欧美一区二区av| 国产午夜精品论理片| 亚洲性久久影院| 三级经典国产精品| 看黄色毛片网站| 亚洲国产日韩欧美精品在线观看| 国产伦一二天堂av在线观看| 免费人成在线观看视频色| 日本一本二区三区精品| 身体一侧抽搐| 亚洲不卡免费看| 干丝袜人妻中文字幕| 亚洲四区av| 51国产日韩欧美| 亚洲精品成人av观看孕妇| 一本一本综合久久| 亚洲婷婷狠狠爱综合网| 久久6这里有精品| 一级二级三级毛片免费看| 欧美精品国产亚洲| 天堂√8在线中文| 精品人妻视频免费看| 哪个播放器可以免费观看大片| 亚洲激情五月婷婷啪啪| 亚洲欧美一区二区三区黑人 | 免费大片18禁| 美女xxoo啪啪120秒动态图| 国产欧美另类精品又又久久亚洲欧美| 国产一区二区在线观看日韩| 亚洲国产精品专区欧美| 纵有疾风起免费观看全集完整版 | 你懂的网址亚洲精品在线观看| 成人综合一区亚洲| 国产精品麻豆人妻色哟哟久久 | 日韩av不卡免费在线播放| 国产免费福利视频在线观看| 人妻夜夜爽99麻豆av| 1000部很黄的大片| 日韩av在线免费看完整版不卡| 免费观看a级毛片全部| 国内揄拍国产精品人妻在线| 天天躁夜夜躁狠狠久久av| 男女啪啪激烈高潮av片| 国产乱人视频| 我的老师免费观看完整版| 人妻夜夜爽99麻豆av| 亚洲av日韩在线播放| 日日啪夜夜撸| 欧美日韩精品成人综合77777| 亚洲av成人精品一二三区| 日韩不卡一区二区三区视频在线| 国产 一区精品| 久久精品夜夜夜夜夜久久蜜豆| 欧美精品国产亚洲| 美女xxoo啪啪120秒动态图| 国精品久久久久久国模美| 插阴视频在线观看视频| 亚洲欧美中文字幕日韩二区| 久久精品久久久久久久性| 麻豆精品久久久久久蜜桃| 亚洲久久久久久中文字幕| videos熟女内射| 国产 一区 欧美 日韩| 亚洲人成网站高清观看| 国产色婷婷99| 少妇的逼水好多| 在线播放无遮挡| 看免费成人av毛片| av.在线天堂| 边亲边吃奶的免费视频| 国产一级毛片七仙女欲春2| 天天躁夜夜躁狠狠久久av| ponron亚洲| 国产高潮美女av| 身体一侧抽搐| 男人舔奶头视频| 国内少妇人妻偷人精品xxx网站| 肉色欧美久久久久久久蜜桃 | 亚洲丝袜综合中文字幕| av在线老鸭窝| 日韩一区二区视频免费看| 久久精品熟女亚洲av麻豆精品 | 成人亚洲精品一区在线观看 | 汤姆久久久久久久影院中文字幕 | 婷婷六月久久综合丁香| 九九久久精品国产亚洲av麻豆| 99久久精品热视频| 久久这里只有精品中国| 又粗又硬又长又爽又黄的视频| 亚洲精品国产成人久久av| 午夜亚洲福利在线播放| 女人十人毛片免费观看3o分钟| 男插女下体视频免费在线播放| 成人性生交大片免费视频hd| 国产麻豆成人av免费视频| 色视频www国产| av在线老鸭窝| 久久久久久久久中文| 成人一区二区视频在线观看| 成人综合一区亚洲| 日日摸夜夜添夜夜添av毛片| 成人性生交大片免费视频hd| 在线播放无遮挡| 成人毛片60女人毛片免费| 亚洲国产精品sss在线观看| 亚洲国产欧美在线一区| 久久久久免费精品人妻一区二区| 九草在线视频观看| 波野结衣二区三区在线| av在线蜜桃| 国产黄频视频在线观看| 美女黄网站色视频| 中文精品一卡2卡3卡4更新| 久久久成人免费电影| 精品久久久久久电影网| 精品久久久久久久久久久久久| 禁无遮挡网站| 亚洲成人一二三区av| 国产激情偷乱视频一区二区| 丝袜喷水一区| 搞女人的毛片| 欧美三级亚洲精品| 亚洲国产精品sss在线观看| 日韩欧美三级三区| 在线观看美女被高潮喷水网站| 女人十人毛片免费观看3o分钟| 国产人妻一区二区三区在| 久久精品国产亚洲av天美| 中国美白少妇内射xxxbb| 黄色日韩在线| 午夜日本视频在线| 亚洲欧美日韩无卡精品| 欧美精品国产亚洲| 亚洲aⅴ乱码一区二区在线播放| 久久久久网色| 中文字幕亚洲精品专区| 国产伦精品一区二区三区四那| 看免费成人av毛片| 亚洲精品乱码久久久久久按摩| 久久精品久久久久久噜噜老黄| 亚洲18禁久久av| 久久6这里有精品| 亚洲精品日本国产第一区| 国产69精品久久久久777片| 日日啪夜夜爽| 国产在线男女| 三级男女做爰猛烈吃奶摸视频| 久久精品人妻少妇| 亚洲精品久久久久久婷婷小说| 亚洲伊人久久精品综合| 国产成人福利小说| 一二三四中文在线观看免费高清| 国产精品美女特级片免费视频播放器| 极品少妇高潮喷水抽搐| 熟妇人妻久久中文字幕3abv| 欧美区成人在线视频| 日韩av不卡免费在线播放| 少妇丰满av| 日韩强制内射视频| 春色校园在线视频观看| 三级国产精品片| 在线观看人妻少妇| 国产成人午夜福利电影在线观看| 色网站视频免费| 汤姆久久久久久久影院中文字幕 | 又大又黄又爽视频免费| 国产精品女同一区二区软件| 免费观看的影片在线观看| 亚洲一区高清亚洲精品| 国产成人精品一,二区| 99re6热这里在线精品视频| 天堂中文最新版在线下载 | 亚洲在线观看片| 久久精品国产亚洲av涩爱| 成人毛片60女人毛片免费| 国内揄拍国产精品人妻在线| 国国产精品蜜臀av免费| 天天躁日日操中文字幕| 国产精品一二三区在线看| 高清在线视频一区二区三区| 秋霞伦理黄片| 三级国产精品欧美在线观看| 最新中文字幕久久久久| 综合色av麻豆| 狂野欧美白嫩少妇大欣赏| 欧美成人午夜免费资源| 国产亚洲精品av在线| 欧美三级亚洲精品| 国产精品av视频在线免费观看| 国产男人的电影天堂91| 色播亚洲综合网| 韩国av在线不卡| 成人美女网站在线观看视频| 亚洲欧美日韩东京热| 人人妻人人澡欧美一区二区| 国产一区二区亚洲精品在线观看| 99久国产av精品国产电影| 晚上一个人看的免费电影| 国产精品人妻久久久久久| 人体艺术视频欧美日本| 精品国产三级普通话版| a级毛片免费高清观看在线播放| 国产精品嫩草影院av在线观看| 中文字幕制服av| 亚洲精品一二三| 亚洲av日韩在线播放| 久久久久久久久久黄片| 亚洲高清免费不卡视频| 99热全是精品| 色5月婷婷丁香| 国产精品一区二区性色av| 99热网站在线观看| 国产高清国产精品国产三级 | 91狼人影院| 日韩av不卡免费在线播放| 国产 亚洲一区二区三区 | 搡老妇女老女人老熟妇| 啦啦啦啦在线视频资源| 免费av观看视频| 亚洲成人av在线免费| 精华霜和精华液先用哪个| 校园人妻丝袜中文字幕| 久久精品国产亚洲网站| 日韩电影二区| 欧美不卡视频在线免费观看| 日韩中字成人| 观看美女的网站| 成人鲁丝片一二三区免费| 午夜福利在线观看吧| 国产v大片淫在线免费观看| 亚洲人成网站在线观看播放| 亚洲av国产av综合av卡| 精品久久国产蜜桃| 久99久视频精品免费| 久久久久久久久久久丰满| 最近中文字幕2019免费版| 久久99热这里只有精品18| 亚洲欧美日韩卡通动漫| 如何舔出高潮| 国内精品一区二区在线观看| 国产精品国产三级专区第一集| 久久久国产一区二区| 亚洲综合精品二区| 嘟嘟电影网在线观看| 国产伦在线观看视频一区| 深夜a级毛片| 校园人妻丝袜中文字幕| 极品教师在线视频| 精品少妇黑人巨大在线播放| 午夜福利高清视频| 精品一区二区免费观看| 色播亚洲综合网| 免费黄频网站在线观看国产| 色网站视频免费| 亚洲国产欧美在线一区| 午夜激情福利司机影院| 欧美人与善性xxx| 看非洲黑人一级黄片| 大片免费播放器 马上看| 极品少妇高潮喷水抽搐| 一级毛片黄色毛片免费观看视频| 日日撸夜夜添| 毛片女人毛片| 国产成人精品一,二区| 亚洲精品亚洲一区二区| 综合色av麻豆| 午夜爱爱视频在线播放| 午夜视频国产福利| 天堂中文最新版在线下载 | 黄色日韩在线| 99久久中文字幕三级久久日本| 午夜福利视频1000在线观看| 亚洲人成网站高清观看| 国产高潮美女av| 深夜a级毛片| av卡一久久| 男女那种视频在线观看| 午夜福利高清视频| 亚洲国产欧美人成| 啦啦啦啦在线视频资源| 免费大片黄手机在线观看| 日本午夜av视频| 午夜福利视频精品| 精品国内亚洲2022精品成人| 一区二区三区高清视频在线| 麻豆国产97在线/欧美| 国产麻豆成人av免费视频| 国产又色又爽无遮挡免| 欧美日韩在线观看h| 久久午夜福利片| 极品少妇高潮喷水抽搐| 麻豆久久精品国产亚洲av| 色综合站精品国产| videossex国产| 中文字幕亚洲精品专区| 成人亚洲精品一区在线观看 | 日本爱情动作片www.在线观看| 欧美日韩亚洲高清精品| 欧美区成人在线视频| 国产黄a三级三级三级人| 男人狂女人下面高潮的视频| 大香蕉久久网| www.色视频.com| 亚洲欧美清纯卡通| 国产乱人偷精品视频| 欧美xxⅹ黑人| 亚洲欧洲日产国产| 三级国产精品片| 久久久久久久久久成人| 久久精品久久久久久噜噜老黄| 亚洲欧美清纯卡通| av专区在线播放| 日韩中字成人| 久久97久久精品| videos熟女内射| 91久久精品电影网| 激情五月婷婷亚洲| 欧美高清性xxxxhd video| 男的添女的下面高潮视频| 深爱激情五月婷婷| 99热这里只有精品一区| 永久网站在线| 乱人视频在线观看| 狠狠精品人妻久久久久久综合| 人妻夜夜爽99麻豆av| 综合色丁香网| 精品亚洲乱码少妇综合久久| 亚洲精品视频女| 免费观看性生交大片5| 亚洲人成网站高清观看| 日本与韩国留学比较| videossex国产| 午夜免费观看性视频| 亚洲欧美日韩无卡精品| 最后的刺客免费高清国语| 久久精品熟女亚洲av麻豆精品 | 舔av片在线| 中文乱码字字幕精品一区二区三区 | 精品久久久久久久久久久久久| 亚洲精品自拍成人| 啦啦啦中文免费视频观看日本| 精品久久久噜噜| 91av网一区二区| 亚洲综合精品二区| 伦精品一区二区三区| 国产精品综合久久久久久久免费| 麻豆国产97在线/欧美| 91精品国产九色| 男的添女的下面高潮视频| 最近中文字幕高清免费大全6| 国产成人精品一,二区| 国产精品福利在线免费观看| 一区二区三区免费毛片| 在线播放无遮挡| 日韩制服骚丝袜av| 丰满少妇做爰视频| 在线播放无遮挡| 免费不卡的大黄色大毛片视频在线观看 | 亚洲精品日韩在线中文字幕| 好男人视频免费观看在线| 欧美激情国产日韩精品一区| 又粗又硬又长又爽又黄的视频| 内地一区二区视频在线| 自拍偷自拍亚洲精品老妇| 亚洲人与动物交配视频| 日本熟妇午夜| 少妇的逼水好多| 美女主播在线视频| 日韩一区二区三区影片| 亚洲第一区二区三区不卡| 国产成人a∨麻豆精品| 中文资源天堂在线| 久久这里只有精品中国| 51国产日韩欧美| 色5月婷婷丁香| 亚洲欧洲日产国产| 日韩欧美国产在线观看| 在线免费观看的www视频| 日韩欧美 国产精品| 亚洲精品视频女| 高清日韩中文字幕在线| 国产亚洲av片在线观看秒播厂 | 日韩av不卡免费在线播放| 久久久色成人| 大话2 男鬼变身卡| 午夜免费激情av| 成人亚洲精品一区在线观看 | 亚洲人成网站在线观看播放| 久久久久久久久久黄片| 国产色爽女视频免费观看| 婷婷六月久久综合丁香| 男人狂女人下面高潮的视频| 男人舔奶头视频| 国产淫片久久久久久久久| 夜夜看夜夜爽夜夜摸| 大片免费播放器 马上看| 97超碰精品成人国产| 成人特级av手机在线观看| 美女被艹到高潮喷水动态| 久久久国产一区二区| 欧美一级a爱片免费观看看| 国产欧美另类精品又又久久亚洲欧美| 欧美高清性xxxxhd video| 亚洲成人中文字幕在线播放| 啦啦啦啦在线视频资源| 午夜福利在线在线| 在线免费观看不下载黄p国产| 午夜福利高清视频| 国产片特级美女逼逼视频| 国产成人精品久久久久久| av国产久精品久网站免费入址| 99久久九九国产精品国产免费| 日日干狠狠操夜夜爽| av在线亚洲专区| 深爱激情五月婷婷| 国产亚洲精品av在线| 激情 狠狠 欧美| 免费看不卡的av| 一级毛片电影观看| 免费少妇av软件| 国产精品一区二区三区四区久久| 十八禁国产超污无遮挡网站| 最近最新中文字幕免费大全7| 日韩欧美国产在线观看| 亚洲综合色惰| 国产 一区 欧美 日韩| 亚洲怡红院男人天堂| 熟妇人妻不卡中文字幕| 亚洲欧洲国产日韩| 亚洲国产精品国产精品| 国产午夜精品久久久久久一区二区三区| 日韩一区二区三区影片| 亚洲伊人久久精品综合| 国产欧美日韩精品一区二区| 久久6这里有精品| 亚洲av日韩在线播放| 亚洲成人一二三区av| 美女黄网站色视频| 亚洲精品一二三| 精品久久久久久久末码| 日韩av在线免费看完整版不卡| 狠狠精品人妻久久久久久综合| 久久人人爽人人爽人人片va| 久久鲁丝午夜福利片| 老司机影院成人| 美女大奶头视频| 亚洲在线观看片| 一夜夜www| 欧美三级亚洲精品| 婷婷色av中文字幕| 亚洲精品色激情综合| 欧美+日韩+精品| 久久国内精品自在自线图片| 亚洲欧美日韩无卡精品| 免费少妇av软件| 嫩草影院新地址| 精品一区二区三区人妻视频| 青春草国产在线视频| 亚洲av男天堂| 久久久久免费精品人妻一区二区| 国产在视频线精品| 欧美日韩在线观看h| 久久久久久久亚洲中文字幕| 在线观看美女被高潮喷水网站| 免费观看av网站的网址| 成人鲁丝片一二三区免费| 亚州av有码| 伊人久久精品亚洲午夜| 国产老妇伦熟女老妇高清| av免费观看日本| 成人一区二区视频在线观看| 91久久精品国产一区二区成人| 国内少妇人妻偷人精品xxx网站| 欧美 日韩 精品 国产| 午夜激情久久久久久久| 大香蕉久久网| 欧美xxxx性猛交bbbb| 久久久成人免费电影| 亚洲国产日韩欧美精品在线观看| 男插女下体视频免费在线播放| 中文欧美无线码| 丝袜喷水一区| 视频中文字幕在线观看| 两个人的视频大全免费| 亚洲久久久久久中文字幕| av一本久久久久| 99视频精品全部免费 在线| 麻豆av噜噜一区二区三区| 爱豆传媒免费全集在线观看| 国产成人freesex在线| 亚洲欧美中文字幕日韩二区| 搡女人真爽免费视频火全软件| 人人妻人人看人人澡| 亚洲综合精品二区| 精品人妻视频免费看| 精品国内亚洲2022精品成人| 午夜福利视频精品| 汤姆久久久久久久影院中文字幕 | 国产午夜精品久久久久久一区二区三区| 日本一本二区三区精品| 国产高清有码在线观看视频| 在线观看av片永久免费下载| 欧美高清性xxxxhd video| 久久久精品欧美日韩精品| 欧美三级亚洲精品| 亚洲av不卡在线观看| 一级a做视频免费观看| av天堂中文字幕网| 啦啦啦啦在线视频资源| 国产国拍精品亚洲av在线观看| 国产成人精品福利久久| 秋霞在线观看毛片| 日韩三级伦理在线观看| 久久久久国产网址| 亚洲人与动物交配视频| 免费电影在线观看免费观看| 亚洲欧美日韩东京热| 午夜激情欧美在线| 久久久久久久久久黄片| 国产精品一区二区三区四区免费观看| 狂野欧美白嫩少妇大欣赏|