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

    西江中游突發(fā)性水體污染水庫群應(yīng)急調(diào)控初探

    2019-10-17 02:37:54黃代忠方神光張義樂
    人民珠江 2019年9期
    關(guān)鍵詞:大藤峽電站峰值

    黃代忠,方神光,張義樂

    (珠江水利科學(xué)研究院,水利部珠江河口海岸工程技術(shù)研究中心,廣東廣州510611)

    水是生命之源,但隨著工業(yè)和第三產(chǎn)業(yè)的粗獷式發(fā)展,水污染事故正成為危害人民生存健康的重大環(huán)境問題。由于突發(fā)性水污染事故具有發(fā)生、發(fā)展、危害的不確定性,流域性、影響的長期性和應(yīng)急主體不明確等顯著特征[1-3],一旦發(fā)生,不僅對人群安全與生存帶來極大危害,還會對生態(tài)環(huán)境造成極大影響[4]。據(jù)統(tǒng)計從2000—2011年期間就發(fā)生1 176起突發(fā)性水污染事故。如2004年2月沱江支流由于約2 000 t純氨氮排入,造成江水嚴重污染;2005年11月苯類污染物污染松花江重大水污染事故和2005年12月粵北北江流域嚴重鎘超標水污染事故等。這一系列水污染事故的發(fā)生,給人民的身體健康和財產(chǎn)造成了巨大損失,因此突發(fā)性水體污染的應(yīng)對和解決方法已成為當前水環(huán)境研究的熱點問題。

    在修建有梯級水庫的中上游河道,通過梯級水庫的聯(lián)合調(diào)度來應(yīng)對突發(fā)性水體污染事故的方法已經(jīng)在實際情況中得到初步使用并取得了不錯的結(jié)果。如2005年北江鎘污染事件中采用了加大上游水庫排量稀釋水體污染物、利用人工小洪峰加快污染物運移到下游處置區(qū)的應(yīng)急措施以控制污染事故的惡化;2012年廣西龍江發(fā)生的錫污染事件中,通過調(diào)度柳江上游大埔、麻石等水庫水量來進行稀釋,保證下游柳州市的飲水安全。由于突發(fā)性水污染事件頻繁發(fā)生,對于該領(lǐng)域的研究也得到快速發(fā)展。方神光等[5]對西江中游突發(fā)性污水團沿河運動特性進行了分析,提出了其影響變化范圍主要分為隨流剪切擴散階段和稀釋影響衰減階段。辛小康等[6]建立MIKE21[7]模型,以宜昌江段為研究對象,計算選取了三峽水庫對三種排放類型不同的調(diào)水方式對污染物輸移傳播的影響;于磊等[8]對大寧水庫發(fā)生突發(fā)性水污染事故時,綜合考慮風(fēng)、水庫調(diào)度的因素,設(shè)定了5種場景來模擬污染物擴散的規(guī)律。為了解西江中游水庫群對污染物輸移擴散的影響,尋找合適的梯級水庫調(diào)度方案以應(yīng)對突發(fā)性水污染事故,本文建立了西江中游河網(wǎng)水動力及污染物輸移擴散數(shù)學(xué)模型,通過不同調(diào)度方式下污染物擴散結(jié)果的對比分析,探討梯級水庫在突發(fā)性水體污染下的合理調(diào)度模式。

    1 數(shù)學(xué)模型

    1.1 河道梯級水庫水動力及污染物輸移擴散數(shù)學(xué)模型

    水動力及污染物輸移擴散數(shù)學(xué)模型如下:

    (1)

    (2)

    (3)

    式中s——距離,m;Q——流量,m3/ s;Z——水位,m;g——重力加速度,m/s2;t——時間,s;A——過水斷面面積,m2;R——水力半徑,m;n——河道糙率系數(shù);D——擴散系數(shù),m2/s;C——某斷面t時刻的濃度,mg/m3;q——網(wǎng)格單位長度上的源或匯,m3/(s·m);cl——源項或匯項的污染物濃度。

    為保障污染物輸移擴散模型的迎風(fēng)特性,在對式(3)中對流項進行離散和計算過程中,需要判斷每個網(wǎng)格節(jié)點位置的流向,同時河道水系中的污染物輸移擴散主要以對流剪切為主,紊動擴散影響相對微小,因此在計算過程中進行了忽略處理。一維水動力數(shù)學(xué)模型采用Preissmann四點隱式差分格式進行離散,通過追趕法對離散得到的大型對角稀疏方程組進行求解。

    1.2 模型驗證

    通過水庫調(diào)度來改變污染物輸移轉(zhuǎn)化的特性,能在很大程度上減輕污染物的危害程度和危害范圍。但由于突發(fā)性水污染事故具有發(fā)生、發(fā)展等特性,真正通過大型水庫及梯級水庫群調(diào)水的方式來處理水污染的案例和研究很少。為驗證此污染物輸移擴散數(shù)學(xué)模型的準確性和可靠性,本文通過對相同工況下數(shù)值模擬和解析公式的結(jié)果分析,對模型進行驗證。

    設(shè)置一個長10 km,底寬67.5 m,邊坡1∶2.5,底坡為0.15‰的棱柱形梯形河道,河道糙率系數(shù)取0.027,上邊界條件為保持恒定流量2 000 m3/s,下邊界根據(jù)明渠均勻流公式推算到的水深確定水位并保持恒定不變。假定在河道上游某斷面突然泄露了1 t的可溶性難降解污染物,試估算下游5、10 km斷面處污染物濃度的變化過程以及0.5、1 h后河道的污染物濃度分布[9]。此類問題可給出輸移擴散式(3)的解析解如下:

    (4)

    式中u——斷面平均流速,m/s;x——下游距離污染源的距離,m;M0——進入水體污染初始面源強度,g/m2;k——污染物衰減系數(shù),s-1。

    忽略污染物在垂向和橫向的擴散,推算斷面摩阻流速0.11 m/s;初始面源強度934.929 g/m,縱向擴散系數(shù)D=6.01hu*=7.4 m2/s。圖1給出了下游不同控制斷面位置污染物濃度隨時間的變化以及不同時刻污染物濃度沿程分布的模型解和解析解比較,驗證了建立的污染物輸移擴散數(shù)學(xué)模型的準確性和可靠性。

    a) 斷面濃度變化過程線

    b) 不同時間點污染物濃度沿程分布圖1 模擬解與解析解的對比

    1.3 西江中游河網(wǎng)水動力及污染物輸移數(shù)學(xué)模型建立

    西江為典型的樹狀河網(wǎng),主干河流有紅水河、黔江、潯江;柳江、郁江為主要支流。由于西江中游梯級水庫的建設(shè),河道上下游水流失去了水力連續(xù)性,圣維南方程組在此處不能適用,為達到同步計算的目的,并融入水庫的調(diào)度方案,此處引入子河網(wǎng)的概念,將水利樞紐作為河網(wǎng)細分為子河網(wǎng)的特殊節(jié)點,以便計算過程中考慮水利樞紐的調(diào)度作用,從而將各個河段連接成一個計算整體,詳細可參見相關(guān)文獻[10]。

    2 紅水河突發(fā)性水污染事故下的水庫調(diào)度分析

    2.1 突發(fā)性水體污染地點及調(diào)度水庫的選定

    大化瑤族自治縣位于廣西壯族自治區(qū)中部偏西北的紅水河中游,隸屬于河池市。大化鎮(zhèn)是該縣的政治、經(jīng)濟、文化、交通中心,是集建材、冶煉、建筑、加工等產(chǎn)業(yè)為一體的工業(yè)城鎮(zhèn)。大化工業(yè)集中區(qū)是由原來的大化縣城南工業(yè)園區(qū)和巖灘工業(yè)區(qū)合并而成,在園內(nèi)集中著以金屬硅、硅錳合金冶煉,電解錳等為主的企業(yè),給當?shù)氐纳鷳B(tài)環(huán)境帶來巨大的隱患。

    本次計算分析假定在大化瑤族自治縣大化鎮(zhèn)發(fā)生突發(fā)性水污染事故,事故地點上距天峨水文站147 km,下距巖灘水電站2.2 km,具體位置見圖2。假設(shè)污染物從2月1日8時開始進入河道,在10時流量達到50 m3/s的最大流量,一直以該流量維持到18時,然后再20時降為0。污水所含的污染物濃度值初始值設(shè)為1,計算范圍內(nèi)污染物濃度的本底值是0,設(shè)河道內(nèi)的污染物濃度值超過0.000 1時,視為超標。

    選取參與調(diào)度的水庫有巖灘、紅花、大藤峽,這些水庫在枯水期調(diào)度的規(guī)則見表1,在不影響水庫運行的前提條件下調(diào)度一定的水量去稀釋及加快污染物的輸移,讓其危害程度降到最低。

    表1 突發(fā)性污染物事故應(yīng)對方案

    2.2 選取的水文條件及還原

    枯季天然來流小,斷面流速小,污染物停留在河道中時間變長,影響距離增大。因此,選擇了2014年2月1日至3月13日的流量、水位實測資料作為水文邊界條件(圖3)。

    區(qū)間入流采用分段還原,將西江中游河網(wǎng)分為天峨—橋鞏段,橋鞏、柳州—武宣段,武宣、南寧—大湟江口、大湟江口—梧州段,梧州—高要段共計5部分。結(jié)合上下游洪量差值還原區(qū)間入流,區(qū)間入流分別設(shè)于刁江、洛清江、蒙江、桂江和賀江入流口。天峨—橋鞏分段還原時,上、下邊界條件分別由天峨站、橋鞏電站實測資料確定,巖灘、大化、百龍灘、樂灘壩前水庫都采用實測資料值,區(qū)間入流設(shè)在刁江口。橋鞏、柳州—武宣分段還原時,上、下游邊界條件分別由遷江、柳州和武宣站實測資料確定,紅花水電站壩前水位為正常蓄水位,區(qū)間入流設(shè)在洛清江口。武宣、南寧—梧州分段還原時,上、下游邊界條件分別由武宣、南寧和梧州站實測資料確定,西津、貴港、桂平長洲樞紐壩前水位由實測資料確定,區(qū)間入流分別設(shè)在蒙江口、桂江口。梧州—高要分段還原時,上、下游邊界條件根據(jù)梧州和高要站實測資料確定,區(qū)間入流設(shè)在賀江口。各河段還原獲得的水位、流量計算值與實測值在允許誤差范圍內(nèi),其中天峨—橋鞏段還原與驗證結(jié)果見圖3。

    圖2 突發(fā)性水污染發(fā)生地點示意

    a) 上下游邊界條件

    b) 區(qū)間入流過程還原

    c) 橋鞏流量驗證圖3 天峨—橋鞏段還原與驗證

    d) 天峨水位驗證續(xù)圖3 天峨—橋鞏段還原與驗證

    2.3 無應(yīng)對措施下污染物輸移特性

    假設(shè)污染物進入河道后,下游各梯級水庫未采取調(diào)度措施進行應(yīng)對,事故發(fā)生地點下游巖灘電站、大化電站、樂灘電站、橋鞏電站、石龍三江口、桂平三江口、梧州水文站點位置污染物濃度變化曲線見圖4。表2給出了各特征斷面污染物濃度變化特征值。

    圖4 無應(yīng)對措施下各控制站點污染物濃度變化過程線

    典型斷面開始超標時刻開始達標時刻出現(xiàn)峰值時刻峰值濃度/(mg·L-1)超標時長/h巖灘電站2月1日15時2月3日9時2月2日4時0.033 644大化電站2月5日17時2月8號9時2月6日17時0.023 866樂灘電站2月10號10時2月12號18時2月11號12時0.019 159橋鞏電站2月13號9時2月15日14時2月14號10時0.017 262石龍三江口2月17號4時2月19號21時2月18號2時0.011 167武宣水文站2月21日3時2月24日6時2月22日13時0.011 076桂平三江口2月27號5時3月2日1時2月28日11時0.004 970梧州水文站3月5號9時3月8日23時3月7號3時0.006 176

    巖灘電站從2月1日15時開始超標,在2月2日4點達到最大值0.033 6,至2月3日9時消失,總共超標歷時44 h;大化電站壩址污染物濃度從2月5日17時開始超標,在2月6日17時達到最大值0.023 8,至2月8號9時開始達標,總共超標歷時66 h;樂灘電站壩址污染物濃度從2月10號10時開始超標,在2月11號12點達到最大值0.019 1,至2月12號18時消失,總共超標歷時59 h;橋鞏電站壩址處濃度在2月13號9時開始超標,在2月14號10時達到最大值0.017 2,至2月15日14時消失,總共超標歷時62 h;石龍三江口污染物濃度從2月17號4時開始超標,在2月18號2時達到最大值0.011 1,至2月19號21時消失,總共超標歷時67 h;武宣水文站位置從2月21日3時開始超標,在2月22日13點達到最大值0.011,至2月24日6點消失,總共超標歷時76 h;桂平三江口從2月27號5時開始超標,在2月28日11時達到最大值0.004 9,至3月2日1時消失,總共超標歷時70 h;梧州水位站位置污染物濃度從3月5號9時開始超標,在3月7號3時達到最大值0.006 1,至3月8日23時消失,總共超標歷時76 h。

    從圖4中看出各典型斷面位置污染物濃度變化呈拋物線狀,拋物線從上游向下游呈現(xiàn)由尖瘦型向矮胖型過渡,顯示各控制站點污染物濃度都呈現(xiàn)某個時間段逐漸增大,達到最大值后逐漸減少直至消失,即污染團通過了該站點。各斷面污染物濃度的超標時間從上游的巖灘電站44 h增大到下游梧州水文站76 h,峰值濃度從0.033 6降低到0.006 1,在樂灘電站總超標歷時減少是由于刁江支流的匯入加快了污染物的輸移,桂平三江口位置總超標歷時減少是郁江水流的匯入加快了污染物的輸移傳播。

    2.4 水庫調(diào)度方案設(shè)置

    2.4.1梯級水庫基本特性

    事故發(fā)生地下游沿線各梯級水庫情況見表3,由表中可知,事故發(fā)生地下游各梯級水庫興利庫容較大的有巖灘、紅花和大藤峽,由于本次模擬事故發(fā)生時間為枯水期,考慮水庫沖污、興利等綜合因素,本次調(diào)度方案主要對巖灘、紅花和大藤峽3個梯級水庫進行調(diào)度。

    2.4.2水庫應(yīng)對調(diào)度方案設(shè)置

    針對該突發(fā)性水污染事故,設(shè)置了4種水庫調(diào)度方案,見表4。

    a) 由于水污染事故發(fā)生的突然性,現(xiàn)實發(fā)生水污染事件時采取的應(yīng)對手段多采用單庫調(diào)度進行應(yīng)對,因此水庫調(diào)度方案1為當巖灘電站壩址處的污染物濃度從2月1日15時開始超標時,巖灘水庫開始調(diào)度。

    表3 事故地點下游各梯級水庫情況

    b) 水庫調(diào)度方案2在方案1調(diào)度巖灘水庫的基礎(chǔ)上,加入紅花水庫一起調(diào)度,為并聯(lián)式雙庫調(diào)度,考慮到紅花調(diào)電站下泄流量與巖灘流量在石龍三江口匯合以增強調(diào)度效果的需要,紅花水庫的起調(diào)時間定為2月12日9時。

    c) 水庫調(diào)度方案3在方案2同時調(diào)度巖灘和紅花水庫的基礎(chǔ)上,加入大藤峽水庫,大藤峽水庫的起調(diào)時間為當壩址處污染物濃度開始超標時,即從2月22日2時。

    d) 大藤峽水庫興利庫容較大,如有需要可進一步增加下泄流量促進污染物輸移擴散,因此水庫調(diào)度方案4在方案3調(diào)度巖灘、紅花和大藤峽水庫的基礎(chǔ)上,將大藤峽水庫調(diào)度下泄流量從方案3的5 000 m3/s增大為6 000 m3/s。

    表4 水庫調(diào)度方案設(shè)置

    2.5 調(diào)度結(jié)果分析

    2.5.1水庫調(diào)度方案1結(jié)果分析

    水庫調(diào)度方案1實施后的計算結(jié)果見表5。方案1單獨調(diào)度巖灘水庫,選擇以增大下泄流量至4 000 m3/s,直至庫內(nèi)水位降至219 m的最低限制水位,調(diào)度時長共計26 h。圖5、6給出了典型斷面在調(diào)度方案1實施前后污染物濃度峰值削減比例和超標時長減少值(正值表示超標時長減少,負值表示增加)。結(jié)果顯示:巖灘電站壩址處水質(zhì)達標的時間從2月3日9時提前到2月2日的1時,污染物濃度的峰值降低了22.3%,污染物超標時間縮短32 h。武宣水文站處的污染物濃度超標時刻由2月21日3時提前到2月16日9時,達標時刻從2月22日13時提前到2月20日8時,污染物濃度的峰值降低了41.8%,但污染物超標時長增加了20 h;梧州水文站處污染物濃度超標時刻由3月5號9時提前到3月2日3時,污染物濃度達標時間從3月8日23時提前到3月6日12時,污染物濃度峰值降低了49.1%,但污染物超標時長增加了32 h。分析調(diào)度方案1得到以下結(jié)論。

    表5 典型斷面污染物濃度特征值(水庫調(diào)度方案1)

    a) 在突發(fā)性水體污染發(fā)生在巖灘庫區(qū)時,通過調(diào)度巖灘水庫,加大下泄流量,可以顯著降低庫區(qū)及下游典型斷面位置處的污染物濃度峰值,起到稀釋污染物的作用,典型污染物峰值削減幅度總體呈現(xiàn)向下游增大趨勢。

    b) 由于加大巖灘出庫下泄流量,相應(yīng)增加了巖灘壩址下游河道的流速,下泄水流攜帶污染物以更快的速度向下游運動,巖灘水庫下游典型斷面位置污染物濃度超標時刻及達標時刻都顯著提前。

    c) 巖灘水庫調(diào)度,可以顯著削減污染物在巖灘庫區(qū)和大化庫區(qū)滯留的時長,但大化壩址以下河道典型斷面的污染物超標時長呈現(xiàn)顯著增加趨勢。

    圖5 調(diào)度方案1下的典型斷面污染物濃度峰值削減比例

    圖6 調(diào)度方案1下的典型斷面污染物濃度超標時長削減值

    2.5.2其他水庫調(diào)度結(jié)果分析

    表6給出了水庫調(diào)度方案2實施后的主要斷面濃度特征值結(jié)果。由于調(diào)度方案2是在方案1調(diào)度巖灘水庫調(diào)度的基礎(chǔ)上,加入了紅花電站,紅花電站控制下泄流量為2 000 m3/s,于石龍三江口與巖灘下泄流量匯合,起到降低污染物濃度并快速向下游推移的作用。因此,從武宣水文站斷面污染物濃度變化過程來看,污水團抵達該斷面的時間為2月15日18時,較方案1提前了15 h,污染物峰值濃度也較方案1下降了21.9%,但該斷面污染物濃度超標時長則較方案1增加了11 h。污水團抵達梧州水文站斷面的時間同樣較方案1提前了18 h,污染物峰值濃度較方案1下降了19.4%,但污染物超標時長則增加了7 h。

    表7給出了水庫調(diào)度方案3實施后的主要斷面污染物濃度特征值。該方案是在方案2巖灘和紅花電站聯(lián)合調(diào)度的基礎(chǔ)上,加入即將建成的大藤峽水庫,大藤峽水庫調(diào)度開始時間為污水團抵達壩址斷面的時刻,大藤峽水庫調(diào)度期間控制下泄流量5 000 m3/s。在大藤峽水庫參與調(diào)度的情形下,可以看到,污水團抵達大藤峽斷面的時間與方案2一致,但由于大藤峽水庫通過加大下泄流量的方式加入調(diào)度,加快了庫區(qū)水體的流動,因此污水團峰值較方案2提前25 h抵達壩址斷面,污染物峰值濃度與方案2相差不大,但該斷面污染物濃度超標時長較方案2顯著減少了44 h;污水團抵達梧州水文站斷面較方案2提前了26 h,污染物濃度峰值較方案2提前18 h抵達,峰值濃度相差不大,但超標時長較方案2顯著減少17 h,在峰值濃度差別不大的情況下,顯然梧州斷面污染物濃度超標時長減少是由于大藤峽水庫加大下泄流量從而對污染物起到極大稀釋作用形成的。

    表8給出了調(diào)度方案4下的大藤峽壩址斷面和梧州水文站斷面的污染物特征值。該方案與調(diào)度方案3基本一致,只是大藤峽電站調(diào)度期間下泄流量比方案3多1 000 m3/s??傮w來看,該方案實施后的變化趨勢與方案3一致,污水團及濃度峰值更快的抵達大藤峽壩址斷面和梧州水文站斷面,峰值濃度基本一樣,兩斷面污染物濃度超標時長都進一步縮短。

    表6 典型斷面污染物濃度特征值(方案2)

    表7 典型斷面污染物濃度特征值(方案3)

    表8 典型斷面污染物濃度特征值(方案4)

    3 結(jié)論

    本文在假定大化瑤族自治縣大化鎮(zhèn)發(fā)生突發(fā)性水污染事故的前提下,設(shè)置不同的調(diào)度方案,聯(lián)合調(diào)度西江流域上一座或某幾座水庫來應(yīng)對該事故,得到的結(jié)論如下。

    a) 單獨調(diào)度巖灘水庫,加大下泄流量??杉铀傥廴疚锵乱扑俣龋@著降低庫區(qū)及下游典型斷面位置處的污染物濃度峰值,削減污染物在巖灘庫區(qū)和大化庫區(qū)滯留時長;同時使巖灘水庫下游典型斷面位置污染物濃度超標時刻及達標時刻都顯著提前;但大化壩址以下河道典型斷面的污染物超標時長呈現(xiàn)顯著增加趨勢。當增大巖灘下泄流量至4 000 m3/s時,巖灘電站壩前水質(zhì)污染物超標時間縮短32 h,污染物濃度的峰值降低了22.3%;武宣水文站斷面污染物濃度的峰值降低41.8%,超標時長增加了20 h;梧州水文站斷面污染物濃度峰值降低49.1%,污染物超標時長增加32 h。

    b) 同時調(diào)度巖灘水庫和紅花電站??蛇M一步降低石龍三江口下游典型斷面位置處的污染物濃度峰值;促進下游典型斷面位置污染物濃度超標時刻及達標時刻的提前;也會進一步延長下游河道典型斷面的污染物超標時長。當增加紅花電站控制下泄流量為2 000 m3/s時,武宣水文站斷面污染物超標時間較方案1提前15 h,超標時長增加11 h,污染物峰值濃度下降21.9%;梧州水文站斷面污染物超標時間提前18 h,超標時長增加7 h,污染物峰值濃度下降19.4%。

    c) 同時調(diào)度巖灘水庫、紅花電站和大藤峽水庫。大藤峽上游污染物移動未發(fā)生改變,但是大藤峽壩址處污染物峰值到達時間提前,持續(xù)時間降低;大藤峽水庫下游各典型斷面污染物到達時間提前,持續(xù)時間降低,峰值濃度不變。當增加大藤峽水庫控制下泄流量5 000 m3/s時,大藤峽壩前水質(zhì)超標時間較方案2提前25 h,超標時長減少了44 h,污染物峰值濃度基本不變;梧州水文站斷面水質(zhì)超標時間較方案2提前26 h,超標時長減少17 h,峰值濃度相差不大。

    d) 同時調(diào)度巖灘水庫、紅花電站和大藤峽水庫,增加大藤峽水庫下泄流量。大藤峽壩址及下游各典型斷面污染物移動速度進一步加快,但是峰值濃度不變。當大藤峽水庫控制下泄流量6 000 m3/s時,大藤峽壩前水質(zhì)超標時間較方案3提前11 h,超標時長減少了8 h,污染物峰值濃度基本不變;梧州水文站斷面水質(zhì)超標時間較方案3提前22 h,超標時長減少6 h,峰值濃度相差不大。

    e) 對比以上4種調(diào)度方式中污染物運動情況可知,當西江中游發(fā)生突發(fā)性水體污染時:①對于干流沿線的水庫,充分利用各梯級水庫的可利用庫容,根據(jù)污染物移動的具體情況,上游水庫向下游逐級補水;②對于支流上的水庫,通過調(diào)控使支流水庫的調(diào)度流量過程傳遞到三江匯合口與污染物輸移傳播的過程保持同步。通過聯(lián)合調(diào)度西江中游各梯級水庫,可以取得比單庫調(diào)度更好的效果。

    猜你喜歡
    大藤峽電站峰值
    “四單”聯(lián)動打造適齡兒童隊前教育峰值體驗
    少先隊活動(2022年9期)2022-11-23 06:55:52
    基于R/S分析的大藤峽出山徑流灰色預(yù)測
    人民珠江(2022年10期)2022-10-28 07:44:56
    三峽電站再創(chuàng)新高
    大跨度、超高雙向門機快速安裝工藝解析
    低影響開發(fā)(LID)在光伏電站中的應(yīng)用
    大藤峽放歌
    歌海(2016年6期)2017-01-10 01:35:42
    寬占空比峰值電流型準PWM/PFM混合控制
    基于峰值反饋的電流型PFM控制方法
    大藤峽公司已撥付征地移民安置資金11.69億元
    中國水利(2015年13期)2015-02-28 15:14:14
    更正
    很黄的视频免费| 国产无遮挡羞羞视频在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 日本黄色日本黄色录像| 19禁男女啪啪无遮挡网站| 亚洲av电影在线进入| 亚洲av熟女| 又大又爽又粗| 成年女人毛片免费观看观看9 | 嫩草影视91久久| 久久久国产欧美日韩av| 日韩熟女老妇一区二区性免费视频| av超薄肉色丝袜交足视频| 人人妻人人澡人人看| 50天的宝宝边吃奶边哭怎么回事| 亚洲精品一卡2卡三卡4卡5卡| 日韩免费av在线播放| 18禁国产床啪视频网站| 国产有黄有色有爽视频| 精品一区二区三区四区五区乱码| 黄色毛片三级朝国网站| 亚洲五月色婷婷综合| 亚洲伊人色综图| 黄色视频不卡| 高清毛片免费观看视频网站 | 亚洲精品粉嫩美女一区| 国产欧美日韩一区二区精品| 亚洲欧洲精品一区二区精品久久久| 手机成人av网站| 老熟妇仑乱视频hdxx| 久久国产乱子伦精品免费另类| 露出奶头的视频| 国产精品 国内视频| 操出白浆在线播放| 亚洲色图综合在线观看| 亚洲欧美精品综合一区二区三区| 在线观看免费视频网站a站| 久久久国产精品麻豆| 母亲3免费完整高清在线观看| 午夜精品在线福利| 亚洲精品av麻豆狂野| 9热在线视频观看99| 免费在线观看影片大全网站| 女人久久www免费人成看片| 啪啪无遮挡十八禁网站| 国产精品亚洲av一区麻豆| 精品国内亚洲2022精品成人 | 一二三四社区在线视频社区8| 亚洲人成77777在线视频| 国产黄色免费在线视频| 久久久国产一区二区| 国产精品国产高清国产av | 大香蕉久久网| 极品人妻少妇av视频| 黄色丝袜av网址大全| 美国免费a级毛片| 午夜免费成人在线视频| 女人久久www免费人成看片| 这个男人来自地球电影免费观看| 成人黄色视频免费在线看| 黄色a级毛片大全视频| 91精品国产国语对白视频| 国产乱人伦免费视频| 午夜福利欧美成人| 男女之事视频高清在线观看| 99久久精品国产亚洲精品| 免费一级毛片在线播放高清视频 | 国精品久久久久久国模美| 一级黄色大片毛片| 日韩欧美一区二区三区在线观看 | 成人国语在线视频| 无限看片的www在线观看| 国产无遮挡羞羞视频在线观看| avwww免费| 欧美在线一区亚洲| 美女扒开内裤让男人捅视频| 9191精品国产免费久久| 国产真人三级小视频在线观看| 国产亚洲欧美在线一区二区| 麻豆国产av国片精品| 欧美最黄视频在线播放免费 | 国产欧美亚洲国产| 欧美 日韩 精品 国产| 黄色视频不卡| 欧美精品人与动牲交sv欧美| 久久久久精品国产欧美久久久| 超碰成人久久| 国产不卡av网站在线观看| 久久人妻熟女aⅴ| 国产精品美女特级片免费视频播放器 | 国产主播在线观看一区二区| 少妇裸体淫交视频免费看高清 | 69av精品久久久久久| 黄色片一级片一级黄色片| 男男h啪啪无遮挡| 成人特级黄色片久久久久久久| 人人妻人人爽人人添夜夜欢视频| 日本wwww免费看| 黑人操中国人逼视频| 亚洲av成人不卡在线观看播放网| а√天堂www在线а√下载 | 国产又色又爽无遮挡免费看| 久久人妻熟女aⅴ| 美女扒开内裤让男人捅视频| 757午夜福利合集在线观看| 无遮挡黄片免费观看| 亚洲av日韩在线播放| 精品久久久久久电影网| 亚洲国产欧美网| 在线视频色国产色| 狠狠婷婷综合久久久久久88av| 亚洲色图 男人天堂 中文字幕| 91老司机精品| 999久久久精品免费观看国产| 亚洲国产精品sss在线观看 | 日韩三级视频一区二区三区| 超色免费av| 怎么达到女性高潮| 女人被躁到高潮嗷嗷叫费观| 18禁国产床啪视频网站| 日本撒尿小便嘘嘘汇集6| 国产精品电影一区二区三区 | 男人的好看免费观看在线视频 | 国产精品久久电影中文字幕 | 丁香欧美五月| 人成视频在线观看免费观看| 法律面前人人平等表现在哪些方面| 女同久久另类99精品国产91| 久久天堂一区二区三区四区| bbb黄色大片| 怎么达到女性高潮| 日本五十路高清| 美女高潮到喷水免费观看| 欧美乱妇无乱码| 免费一级毛片在线播放高清视频 | 日韩免费高清中文字幕av| 在线观看一区二区三区激情| 一夜夜www| 国产精品欧美亚洲77777| 免费日韩欧美在线观看| 99re在线观看精品视频| 一区在线观看完整版| 久久人妻av系列| 成人av一区二区三区在线看| 精品久久久久久久毛片微露脸| 69精品国产乱码久久久| 亚洲在线自拍视频| 久久国产乱子伦精品免费另类| 99精国产麻豆久久婷婷| 欧美日韩福利视频一区二区| 亚洲自偷自拍图片 自拍| 黑人巨大精品欧美一区二区蜜桃| 18禁美女被吸乳视频| 夜夜夜夜夜久久久久| 亚洲av片天天在线观看| 国产91精品成人一区二区三区| 亚洲五月天丁香| 久久久国产成人精品二区 | 下体分泌物呈黄色| 丁香欧美五月| 亚洲人成电影免费在线| 91九色精品人成在线观看| 日韩 欧美 亚洲 中文字幕| 少妇裸体淫交视频免费看高清 | 精品乱码久久久久久99久播| 欧美成狂野欧美在线观看| a级片在线免费高清观看视频| 在线观看免费日韩欧美大片| 首页视频小说图片口味搜索| 国产一区二区三区视频了| 在线免费观看的www视频| 国产欧美日韩综合在线一区二区| 免费观看精品视频网站| www.自偷自拍.com| 亚洲精品自拍成人| 91老司机精品| 亚洲伊人色综图| 18禁美女被吸乳视频| 777米奇影视久久| www日本在线高清视频| 91麻豆精品激情在线观看国产 | 欧美精品人与动牲交sv欧美| 亚洲av片天天在线观看| 黄色 视频免费看| 久久久水蜜桃国产精品网| 女人精品久久久久毛片| 777久久人妻少妇嫩草av网站| 午夜视频精品福利| 90打野战视频偷拍视频| 日本vs欧美在线观看视频| 亚洲成人手机| 国产精品二区激情视频| 美女 人体艺术 gogo| www日本在线高清视频| 久久狼人影院| 嫁个100分男人电影在线观看| 热re99久久国产66热| 国产精品综合久久久久久久免费 | 午夜91福利影院| 99久久人妻综合| 母亲3免费完整高清在线观看| 国产片内射在线| av免费在线观看网站| 国产午夜精品久久久久久| 日本一区二区免费在线视频| 黑人操中国人逼视频| 国产欧美日韩一区二区三| 18在线观看网站| 在线国产一区二区在线| 午夜福利欧美成人| 制服诱惑二区| 成年版毛片免费区| 高清视频免费观看一区二区| 国产精品久久电影中文字幕 | 纯流量卡能插随身wifi吗| 国产亚洲精品久久久久5区| 精品欧美一区二区三区在线| 黄色 视频免费看| 俄罗斯特黄特色一大片| 免费不卡黄色视频| 日日爽夜夜爽网站| 亚洲国产精品sss在线观看 | svipshipincom国产片| 三上悠亚av全集在线观看| videosex国产| 国产在线精品亚洲第一网站| 日韩制服丝袜自拍偷拍| 欧美激情高清一区二区三区| 精品少妇久久久久久888优播| 亚洲全国av大片| 国产精品一区二区在线观看99| 免费不卡黄色视频| 精品少妇一区二区三区视频日本电影| 91av网站免费观看| 黄片播放在线免费| 精品久久久久久,| 如日韩欧美国产精品一区二区三区| 99精品在免费线老司机午夜| 精品人妻在线不人妻| 在线观看免费日韩欧美大片| 好看av亚洲va欧美ⅴa在| 国产91精品成人一区二区三区| 香蕉国产在线看| 国产精品电影一区二区三区 | cao死你这个sao货| 国产亚洲精品一区二区www | 91老司机精品| 每晚都被弄得嗷嗷叫到高潮| 国产黄色免费在线视频| 丝袜美足系列| 国产高清激情床上av| 99久久国产精品久久久| 欧美日韩精品网址| 亚洲精品美女久久av网站| av网站在线播放免费| 在线观看免费视频日本深夜| 久久精品国产亚洲av高清一级| www.精华液| av线在线观看网站| 夜夜爽天天搞| 天天躁日日躁夜夜躁夜夜| 亚洲精品久久成人aⅴ小说| 欧美精品人与动牲交sv欧美| 免费在线观看日本一区| 亚洲全国av大片| 丰满饥渴人妻一区二区三| www.自偷自拍.com| 国产一区在线观看成人免费| 久久人妻福利社区极品人妻图片| 18禁黄网站禁片午夜丰满| aaaaa片日本免费| 在线十欧美十亚洲十日本专区| 久久久国产欧美日韩av| 香蕉丝袜av| 男人操女人黄网站| 日日爽夜夜爽网站| 国产精品国产高清国产av | 国产成人精品久久二区二区免费| 精品一区二区三区视频在线观看免费 | 黄频高清免费视频| 欧美人与性动交α欧美精品济南到| 色尼玛亚洲综合影院| 欧美乱码精品一区二区三区| 午夜视频精品福利| 国产精品九九99| 后天国语完整版免费观看| 9色porny在线观看| 变态另类成人亚洲欧美熟女 | 久久人人爽av亚洲精品天堂| 涩涩av久久男人的天堂| 久久天堂一区二区三区四区| 日韩欧美在线二视频 | av片东京热男人的天堂| 麻豆乱淫一区二区| 老司机影院毛片| 夫妻午夜视频| 韩国av一区二区三区四区| 成年人免费黄色播放视频| 久久久久国产一级毛片高清牌| 成年版毛片免费区| 80岁老熟妇乱子伦牲交| 国产麻豆69| 欧美精品亚洲一区二区| 母亲3免费完整高清在线观看| 免费一级毛片在线播放高清视频 | svipshipincom国产片| videos熟女内射| 18在线观看网站| 激情在线观看视频在线高清 | 极品人妻少妇av视频| 建设人人有责人人尽责人人享有的| 久久国产精品大桥未久av| 脱女人内裤的视频| 国产免费av片在线观看野外av| 国产精品久久视频播放| 亚洲va日本ⅴa欧美va伊人久久| 高清视频免费观看一区二区| 国产熟女午夜一区二区三区| 欧美亚洲日本最大视频资源| 久久天躁狠狠躁夜夜2o2o| 免费在线观看视频国产中文字幕亚洲| 中文字幕人妻丝袜一区二区| 亚洲情色 制服丝袜| 国产亚洲精品第一综合不卡| 欧美日韩视频精品一区| 丝袜在线中文字幕| 亚洲国产欧美一区二区综合| 国产精品乱码一区二三区的特点 | 女人被躁到高潮嗷嗷叫费观| 欧美在线一区亚洲| 女性被躁到高潮视频| 国产欧美亚洲国产| 50天的宝宝边吃奶边哭怎么回事| 久久久精品免费免费高清| 在线av久久热| 动漫黄色视频在线观看| 一级片免费观看大全| 少妇猛男粗大的猛烈进出视频| 日韩三级视频一区二区三区| 亚洲欧洲精品一区二区精品久久久| 在线观看免费视频网站a站| 国产熟女午夜一区二区三区| 王馨瑶露胸无遮挡在线观看| 久久中文字幕人妻熟女| 久久青草综合色| 99热只有精品国产| www.精华液| 久久精品熟女亚洲av麻豆精品| 一级作爱视频免费观看| 午夜两性在线视频| 久久久久久久久久久久大奶| 曰老女人黄片| 久久中文字幕一级| 黑人欧美特级aaaaaa片| 亚洲av美国av| 欧美丝袜亚洲另类 | 亚洲 国产 在线| 中文字幕另类日韩欧美亚洲嫩草| 国产精品乱码一区二三区的特点 | 欧美亚洲日本最大视频资源| 99香蕉大伊视频| 欧美激情久久久久久爽电影 | 亚洲国产欧美网| 久久天堂一区二区三区四区| 99久久精品国产亚洲精品| 丝瓜视频免费看黄片| 这个男人来自地球电影免费观看| 黄色毛片三级朝国网站| 一进一出好大好爽视频| 欧美丝袜亚洲另类 | 亚洲熟女毛片儿| 一a级毛片在线观看| 国产欧美日韩综合在线一区二区| 国产精品99久久99久久久不卡| 90打野战视频偷拍视频| 亚洲国产欧美网| 亚洲精品一二三| 国产人伦9x9x在线观看| 日韩成人在线观看一区二区三区| 99精品久久久久人妻精品| 国产精品综合久久久久久久免费 | 国产高清视频在线播放一区| 精品福利永久在线观看| 美女福利国产在线| 99精国产麻豆久久婷婷| 真人做人爱边吃奶动态| 女人久久www免费人成看片| 成人国语在线视频| 悠悠久久av| 一级片免费观看大全| 九色亚洲精品在线播放| 免费观看人在逋| 91大片在线观看| 亚洲精品一二三| 久久久国产一区二区| 国产主播在线观看一区二区| 久久草成人影院| 日本欧美视频一区| 在线av久久热| 精品视频人人做人人爽| 777米奇影视久久| 在线观看免费视频网站a站| 午夜老司机福利片| 久久久国产欧美日韩av| 国产97色在线日韩免费| av欧美777| 亚洲一区高清亚洲精品| 精品国产一区二区三区四区第35| 国产亚洲欧美98| 无遮挡黄片免费观看| 女警被强在线播放| 国产日韩欧美亚洲二区| 国产熟女午夜一区二区三区| 精品福利观看| 最近最新中文字幕大全免费视频| 每晚都被弄得嗷嗷叫到高潮| 国产成人欧美在线观看 | 一本综合久久免费| 免费黄频网站在线观看国产| 欧美日韩视频精品一区| 免费一级毛片在线播放高清视频 | 国产在线一区二区三区精| 又大又爽又粗| 搡老乐熟女国产| 亚洲avbb在线观看| 一本综合久久免费| 热99久久久久精品小说推荐| 亚洲第一青青草原| 欧美丝袜亚洲另类 | 日日摸夜夜添夜夜添小说| 国产欧美日韩一区二区三| 一级作爱视频免费观看| 亚洲av成人一区二区三| 丁香欧美五月| 久久久久国产一级毛片高清牌| 窝窝影院91人妻| 天天影视国产精品| 欧美一级毛片孕妇| 亚洲熟女毛片儿| 十八禁网站免费在线| 丝袜在线中文字幕| 亚洲国产精品sss在线观看 | 久久久精品区二区三区| 免费在线观看完整版高清| av福利片在线| 在线观看免费视频日本深夜| 国产男女超爽视频在线观看| 婷婷丁香在线五月| 成年人午夜在线观看视频| 亚洲精品美女久久av网站| 久久久精品国产亚洲av高清涩受| 高清欧美精品videossex| 久热这里只有精品99| 黄片大片在线免费观看| 欧美日韩瑟瑟在线播放| 亚洲第一青青草原| 久99久视频精品免费| 一本大道久久a久久精品| av有码第一页| 熟女少妇亚洲综合色aaa.| 女性被躁到高潮视频| 亚洲九九香蕉| 国产人伦9x9x在线观看| ponron亚洲| 免费女性裸体啪啪无遮挡网站| 国产1区2区3区精品| 国产精品久久久av美女十八| 麻豆国产av国片精品| 黑人巨大精品欧美一区二区mp4| 高清黄色对白视频在线免费看| 一本大道久久a久久精品| 国产不卡一卡二| 操美女的视频在线观看| 色在线成人网| 久久久国产成人免费| videos熟女内射| 精品无人区乱码1区二区| 午夜福利欧美成人| 女人被狂操c到高潮| av欧美777| 亚洲国产中文字幕在线视频| 亚洲精品美女久久久久99蜜臀| 精品第一国产精品| 亚洲片人在线观看| 窝窝影院91人妻| 成年人免费黄色播放视频| 成人永久免费在线观看视频| 久久人妻熟女aⅴ| 一边摸一边做爽爽视频免费| 一区二区三区精品91| 久久精品国产a三级三级三级| 久久久精品国产亚洲av高清涩受| 精品视频人人做人人爽| 亚洲精品国产一区二区精华液| 99热只有精品国产| 夜夜躁狠狠躁天天躁| 国产精品av久久久久免费| 精品久久久精品久久久| 电影成人av| 九色亚洲精品在线播放| 亚洲九九香蕉| 精品少妇久久久久久888优播| 黑人猛操日本美女一级片| 日韩免费av在线播放| 一边摸一边抽搐一进一小说 | 久久精品aⅴ一区二区三区四区| 久久久精品国产亚洲av高清涩受| 欧美精品啪啪一区二区三区| 久久久久久久午夜电影 | 精品人妻1区二区| 久久久久久亚洲精品国产蜜桃av| 老司机靠b影院| 18禁黄网站禁片午夜丰满| 在线观看www视频免费| videosex国产| 色综合婷婷激情| 久久精品aⅴ一区二区三区四区| 国产精品亚洲av一区麻豆| 精品久久久久久,| 在线观看免费日韩欧美大片| 欧美人与性动交α欧美软件| av有码第一页| svipshipincom国产片| 极品少妇高潮喷水抽搐| 久久天躁狠狠躁夜夜2o2o| 欧美日韩av久久| 91成年电影在线观看| 国产aⅴ精品一区二区三区波| 黄色女人牲交| 日本vs欧美在线观看视频| 不卡一级毛片| 波多野结衣一区麻豆| 下体分泌物呈黄色| 国产精品自产拍在线观看55亚洲 | 午夜影院日韩av| 日韩大码丰满熟妇| 咕卡用的链子| a在线观看视频网站| 亚洲精品国产区一区二| 中文字幕色久视频| ponron亚洲| 在线观看免费日韩欧美大片| 中文字幕制服av| www.精华液| 亚洲国产毛片av蜜桃av| 国产成人免费观看mmmm| 深夜精品福利| 亚洲色图综合在线观看| 夜夜爽天天搞| 一级黄色大片毛片| 国产亚洲欧美98| 女人高潮潮喷娇喘18禁视频| 久久热在线av| 亚洲精品美女久久久久99蜜臀| 色婷婷久久久亚洲欧美| 夜夜爽天天搞| 人人妻,人人澡人人爽秒播| 精品一区二区三区四区五区乱码| 国产一卡二卡三卡精品| 亚洲色图综合在线观看| 亚洲五月天丁香| 超碰97精品在线观看| 可以免费在线观看a视频的电影网站| 亚洲精品久久成人aⅴ小说| 精品高清国产在线一区| 国产有黄有色有爽视频| 老汉色∧v一级毛片| 女警被强在线播放| 免费观看精品视频网站| 色94色欧美一区二区| 又紧又爽又黄一区二区| 嫁个100分男人电影在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 欧美日韩亚洲国产一区二区在线观看 | 欧美日韩黄片免| 日韩一卡2卡3卡4卡2021年| 久久狼人影院| 日日夜夜操网爽| av有码第一页| 人妻丰满熟妇av一区二区三区 | 日本撒尿小便嘘嘘汇集6| 成人av一区二区三区在线看| 美女扒开内裤让男人捅视频| 他把我摸到了高潮在线观看| 国产成+人综合+亚洲专区| 大型av网站在线播放| 亚洲全国av大片| 午夜老司机福利片| 精品一区二区三卡| 十分钟在线观看高清视频www| 中文字幕人妻丝袜一区二区| 王馨瑶露胸无遮挡在线观看| 欧美成狂野欧美在线观看| 亚洲 国产 在线| 久久人妻av系列| 久久精品人人爽人人爽视色| 在线观看免费视频日本深夜| 久久久国产成人免费| 少妇被粗大的猛进出69影院| 淫妇啪啪啪对白视频| 欧美在线一区亚洲| 久久香蕉精品热| 成熟少妇高潮喷水视频| 久热这里只有精品99| 日韩免费高清中文字幕av| 久久人妻福利社区极品人妻图片| 国产成人免费无遮挡视频| 国产亚洲精品久久久久5区| 久久香蕉精品热| 久久久久久久久久久久大奶| 妹子高潮喷水视频| cao死你这个sao货| 久久久久久免费高清国产稀缺| 国产成人av教育| 午夜福利视频在线观看免费|