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

    松花江流域干旱突變特征及其遙相關(guān)因子分析

    2024-01-11 13:56:00青,劉龍,郝
    人民長(zhǎng)江 2023年12期
    關(guān)鍵詞:松花江流域共振小波

    曹 青,劉 康 龍,郝 振 純

    (1.水利部水文氣象災(zāi)害機(jī)理與預(yù)警重點(diǎn)實(shí)驗(yàn)室/協(xié)同創(chuàng)新中心,江蘇 南京 210044; 2.南京信息工程大學(xué) 水文與水資源工程學(xué)院,江蘇 南京 210044; 3.河海大學(xué) 水文水資源與水利工程科學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室,江蘇 南京 210098)

    0 引 言

    干旱是災(zāi)害性最大的常發(fā)自然災(zāi)害之一,對(duì)農(nóng)業(yè)、水資源供給等人類(lèi)活動(dòng)構(gòu)成了巨大威脅[1]。近年來(lái),全球氣象學(xué)者提出一系列氣象干旱評(píng)價(jià)指標(biāo)來(lái)研究不同地區(qū)的干旱演化特征,如標(biāo)準(zhǔn)化降水蒸散指數(shù)(SPEI)[2-3]、標(biāo)準(zhǔn)化降水指數(shù)(SPI)[4-5]、帕爾默干旱強(qiáng)度指數(shù)[6]等,并與農(nóng)業(yè)、環(huán)境、水利水電工程等學(xué)科交叉結(jié)合以研究干旱對(duì)不同學(xué)科現(xiàn)象的影響特征[7-8]。李華偉[9]通過(guò)Copula函數(shù)開(kāi)發(fā)出的氣象-農(nóng)業(yè)模型模擬得出,新疆谷物產(chǎn)量相應(yīng)較敏感的干旱時(shí)間尺度較長(zhǎng),而棉花所響應(yīng)尺度較短;Cuartas等[10]基于多種干旱指數(shù)評(píng)估了巴西的水資源儲(chǔ)量情況,發(fā)現(xiàn)SPEI可以用于評(píng)估測(cè)站數(shù)量較少的流域干旱事件,且對(duì)巴西干旱形勢(shì)下水力水資源管理提出警醒與理論依據(jù)。與此同時(shí),研究學(xué)者嘗試尋求地區(qū)干旱的年或年代際等多種時(shí)間尺度演變趨勢(shì)的遙相關(guān)因素和動(dòng)力機(jī)制[11-12],劉仲藜等[13]、張劍明等[14]研究發(fā)現(xiàn),洞庭湖流域以及湖南大部分區(qū)域干旱事件皆與西太平洋副熱帶高壓的形勢(shì)有著顯著的相關(guān)影響。更進(jìn)一步地,部分學(xué)者通過(guò)構(gòu)建數(shù)學(xué)模型并對(duì)其不斷修正以期對(duì)未來(lái)干旱現(xiàn)象做出準(zhǔn)確預(yù)警和處理[15-17]。Dikshit等[18]采用標(biāo)準(zhǔn)降水指數(shù)(SPI)和SHAP預(yù)測(cè)算法闡述了氣候變量在月尺度和年尺度不同范圍內(nèi)的重要性,以此證明了在預(yù)測(cè)模型中添加氣候變量指數(shù)作為預(yù)測(cè)因子的必要性。Malik等[19]通過(guò)有效干旱指數(shù)(EDI)研究得出GWO與SHO等元啟發(fā)式算法對(duì)氣象干旱預(yù)測(cè)具有更好的適用性和優(yōu)越性,且混合SVR-GWO模型更優(yōu)于SVR-SHO模型。

    松花江流域是中國(guó)十大流域之一,作為中國(guó)東北地區(qū)重要的農(nóng)業(yè)產(chǎn)區(qū)和重工業(yè)基地,其耕地面積大、土壤肥沃等優(yōu)勢(shì)使得東北地區(qū)農(nóng)業(yè)發(fā)展十分發(fā)達(dá),因此研究松花江流域的干旱特征對(duì)農(nóng)業(yè)災(zāi)害防治預(yù)警有著重大意義[20-21]。趙蘭蘭等[22]通過(guò)VIC模型模擬土壤含水率距平百分率指數(shù),證明了基于該指數(shù)的干旱等級(jí)劃分標(biāo)準(zhǔn)具有較好的適用性。Faiz等[23]研究發(fā)現(xiàn)基于參考蒸散量的干旱指數(shù)能夠捕捉到更多的干旱事件,且多個(gè)干旱指數(shù)在RCP8.5下皆顯示出最高的干旱發(fā)生率。通過(guò)采用WEAP模型,預(yù)測(cè)氣候變化條件下未來(lái)河流徑流可能會(huì)減少24%,且輕、中度干旱發(fā)生概率增加[24]。Song 等[25]探究了松花江干旱事件與極端降水的演變特征,得出松花江流域春季與秋季干旱發(fā)生率較高。

    綜上所述,松花江流域已有的氣象干旱研究主要集中在干旱指數(shù)的時(shí)空演變和多模態(tài)特征研究,以及與農(nóng)業(yè)、植被、環(huán)境等其他領(lǐng)域的相互影響和干旱頻率強(qiáng)度風(fēng)險(xiǎn)分析等方面,而在干旱現(xiàn)象年代際突變點(diǎn)識(shí)別以及突變前后變化差異等方面仍存在不足。在20世紀(jì)80~90年代之后,由于降水減少與全球變暖的因素,中國(guó)北方呈現(xiàn)“暖干化”趨勢(shì)[26],松花江流域作為中國(guó)重要的商品糧基地,厘清干旱指數(shù)是否存在突變及其突變前后與遙相關(guān)因子的響應(yīng)關(guān)系,有利于加深對(duì)松花江流域旱澇災(zāi)害風(fēng)險(xiǎn)趨勢(shì)的理解。因而本研究結(jié)合多種非參數(shù)檢驗(yàn)方法識(shí)別標(biāo)準(zhǔn)化降水蒸散指數(shù)時(shí)間序列突變點(diǎn),并運(yùn)用交叉小波變換與小波相干分析,探究突變前后標(biāo)準(zhǔn)化降水蒸散指數(shù)與遙相關(guān)因子相互作用關(guān)系的差異,具象化流域干旱突變特性的空間分布模式。研究成果可為松花江流域抗旱減災(zāi)、生態(tài)水資源管理、保障國(guó)家糧食安全等工作提供技術(shù)支撐。

    1 研究區(qū)概況與數(shù)據(jù)來(lái)源

    松花江是黑龍江在中國(guó)境內(nèi)最大支流,為中國(guó)七大河流之一,總長(zhǎng)度約達(dá)1 927 km,流域總面積達(dá)到55.68萬(wàn)km2,橫跨內(nèi)蒙古、吉林、黑龍江3個(gè)省份;流域多年平均氣溫為3~5 ℃,年內(nèi)日平均溫度最高月份為7月,可達(dá)25 ℃,1月最低,約為-20 ℃。多年平均降水量約為500 mm。本研究所用氣溫與降雨數(shù)據(jù)為中國(guó)氣象數(shù)據(jù)網(wǎng)中的中國(guó)地面氣候資料日值數(shù)據(jù)集V3.0;北極濤動(dòng)Arctic Oscillation(AO)指數(shù)與太平洋年代際振蕩Pacific Decadal Oscillation(PDO)數(shù)據(jù)來(lái)自NOAA網(wǎng)站(http:∥www.cpc.noaa.gov);El Nio-Southern Oscillation(ENSO)指數(shù)的計(jì)算采用全球海表溫度再分析數(shù)據(jù),來(lái)自英國(guó)氣象局哈德利中心觀測(cè)數(shù)據(jù)集(HadiSST1)。松花江流域地理信息和氣象站點(diǎn)如圖1所示,選取松花江流域52個(gè)氣象站點(diǎn)1960~2020年共61 a逐日降水與氣溫?cái)?shù)據(jù),氣象站點(diǎn)數(shù)量較多且在流域范圍內(nèi)分散均勻,且本研究選取的時(shí)間序列跨度較長(zhǎng),確保了研究分析與結(jié)論的可靠性。

    圖1 松花江流域氣象站分布Fig.1 Meteorological station distribution in Songhua River basin

    2 研究方法

    2.1 遙相關(guān)因子

    遙相關(guān)研究對(duì)于流域氣象旱澇要素變化與擾動(dòng)機(jī)制具有重大意義,研究顯示中國(guó)東北地區(qū)干旱變化的第一二主周期與厄爾尼諾-南方濤動(dòng)事件具有一致的振蕩準(zhǔn)周期[27-28]。北極濤動(dòng)呈負(fù)位相分布、太平洋年代際振蕩指數(shù)較高時(shí),松花江流域西部易雨澇[29-30]。本文選取北極濤動(dòng)、太平洋年代際濤動(dòng)、厄爾尼諾-南方濤動(dòng)作為遙相關(guān)研究因子。

    2.2 標(biāo)準(zhǔn)化降水蒸散指數(shù)

    SPEI是在標(biāo)準(zhǔn)化降水指數(shù)(SPI)的基礎(chǔ)上提出的標(biāo)準(zhǔn)化蒸散指數(shù),具備多時(shí)間尺度特征,計(jì)算簡(jiǎn)便的同時(shí)還引入了影響干旱的潛在蒸發(fā)項(xiàng),能夠更好地刻畫(huà)中國(guó)區(qū)域的干旱事件[31]。沈國(guó)強(qiáng)等[32]、梁豐等[33]、徐一丹等[34]還通過(guò)干旱災(zāi)情數(shù)據(jù)、土壤水分檢測(cè)值與標(biāo)準(zhǔn)化降水蒸散指數(shù)的顯著關(guān)聯(lián)性等多個(gè)角度論證了標(biāo)準(zhǔn)化降水蒸散指數(shù)在東北地區(qū)干旱定量化研究的適用性。

    標(biāo)準(zhǔn)化蒸散指數(shù)首先采用Thornthwaite方法計(jì)算逐月潛在蒸散發(fā)量,將逐月降水量與潛在蒸散量的差值序列進(jìn)行正態(tài)化處理,構(gòu)建成不同時(shí)間尺度Log-logistic概率分布的累計(jì)函數(shù),對(duì)時(shí)間序列進(jìn)行擬合并標(biāo)準(zhǔn)化,得到不同時(shí)間尺度的SPEI指數(shù)序列。本文取12個(gè)月尺度的SPEI指數(shù)(SPEI-12),SPEI干旱等級(jí)分類(lèi)如表1所列。

    表1 SPEI-12干旱等級(jí)分類(lèi)Tab.1 SPEI-12 drought classification

    2.3 突變檢驗(yàn)方法

    (1) Yamamoto法。該方法是從信噪比的角度來(lái)分析突變問(wèn)題,檢驗(yàn)兩子序列均值的差異是否顯著,被廣泛應(yīng)用于識(shí)別氣象要素年代際突變點(diǎn)[35]。

    (2) Lepage法。該方法為無(wú)分布雙樣本的非參數(shù)檢驗(yàn)方法,將序列中的兩個(gè)子序列視為兩個(gè)獨(dú)立主體,若統(tǒng)計(jì)檢驗(yàn)顯示兩個(gè)子序列有顯著差異,則認(rèn)定在兩個(gè)子序列的劃分基準(zhǔn)線點(diǎn)時(shí)刻發(fā)生了突變[36]。

    2.4 交叉小波變換和小波相干變換

    交叉小波變換(XWT)結(jié)合了小波變換與交叉譜分析,可以識(shí)別出兩個(gè)時(shí)間序列在不同時(shí)頻域相互作用的顯著性,并且能識(shí)別出序列間的相關(guān)一致性及在時(shí)頻空間中的相位關(guān)系。設(shè)WX(s),WY(s)為兩個(gè)時(shí)間序列X與Y的交叉小波變換,則定義其交叉小波功率譜為[37]

    (1)

    式中:左項(xiàng)絕對(duì)值即對(duì)應(yīng)交叉功率譜密度,絕對(duì)值越大,表明兩者高能量區(qū)相關(guān)越顯著。

    (2)

    式中:σX,σY分別為時(shí)間序列X,Y的標(biāo)準(zhǔn)差;v為Morlet小波變換自由度,取2,當(dāng)式中左項(xiàng)超過(guò)紅色噪音功率譜的95%的置信限上界,則認(rèn)為通過(guò)了顯著性水平α=0.05的紅色噪音標(biāo)準(zhǔn)譜的檢驗(yàn)[37];Zv(P)是與概率P有關(guān)的置信度,在顯著性水平α=0.05下,Z2(95%)=3.999。

    小波相干變換(WTC)彌補(bǔ)了交叉小波在低能量區(qū)識(shí)別序列間相關(guān)關(guān)系的不足。定義序列X與Y的小波相干譜為

    (3)

    本文使用MATLAB中的小波工具箱(http:∥grinsted.github.io/wavelet-coherence/)來(lái)進(jìn)行小波變換相關(guān)計(jì)算,將處理后的數(shù)據(jù)導(dǎo)入工具箱,從而分析得出交叉小波與小波相干譜圖。

    3 結(jié)果分析

    3.1 干旱指數(shù)突變的時(shí)空分布特征

    結(jié)合松花江流域52個(gè)氣象站點(diǎn)1960~2020年共61 a逐日降水?dāng)?shù)據(jù),計(jì)算逐年SPEI-12指數(shù),使用Yamamoto法、Lepage法識(shí)別并檢驗(yàn)指數(shù)年序列突變點(diǎn),線性擬合突變前后變化趨勢(shì)直線,并對(duì)該年序列取10 a滑動(dòng)平均,結(jié)果如表2、圖2所示。

    表2 SPEI-12指數(shù)突變識(shí)別與檢驗(yàn)Tab.2 Identification and test of the mutation of the SPEI-12

    圖2 SPEI-12與遙相關(guān)因子指數(shù)變化Fig.2 Variation of SPEI-12 and teleconnection factor index

    表2表明,兩種檢驗(yàn)方法都顯示SPEI-12指數(shù)序列最為顯著突變點(diǎn)為1984年,本研究取1984年為松花江流域SPEI-12指數(shù)的突變節(jié)點(diǎn)。

    由圖2可知,SPEI-12指數(shù)突變前后趨勢(shì)變化明顯,突變前擬合斜率約為-0.017,突變后擬合斜率約為-0.046,突變后下降趨勢(shì)增強(qiáng)。

    利用空間插值方法處理所取突變點(diǎn)前后的SPEI-12指數(shù),繪制出地理空間分布圖,如圖3所示。由圖3(a)、(b)可知松花江流域SPEI-12突變前空間上北部與南部少部分區(qū)域數(shù)值在正常范圍之內(nèi),其他大部分區(qū)域數(shù)值處在輕澇級(jí)別;突變后轉(zhuǎn)變?yōu)槌诤邮心喜颗c齊齊哈爾市西北部小部分區(qū)域有輕旱現(xiàn)象外,其他絕大部分區(qū)域數(shù)值處在正常范圍內(nèi),無(wú)明顯旱澇現(xiàn)象。圖3(c)則顯示出松花江流域整體SPEI-12指數(shù)趨于減小,齊齊哈爾市以西北呼倫貝爾以東南地區(qū)、黑河市周邊地區(qū)和吉林北部地區(qū)變化最為顯著。

    圖3 SPEI-12突變前后空間分布Fig.3 Spatial distribution of SPEI-12 before and after mutation

    3.2 突變前后干旱指數(shù)遙相關(guān)因子分析

    3.2.1交叉小波與小波相干譜分析

    本文對(duì)松花江流域SPEI-12指數(shù)突變前后序列與3個(gè)遙相關(guān)因子指數(shù)對(duì)應(yīng)年際變化序列進(jìn)行交叉小波分析及小波相干分析,探究松花江流域降水集中指數(shù)突變前后與遙相關(guān)因子的共振周期、相位關(guān)系與顯著時(shí)段等遙相關(guān)特征的區(qū)別與變化,繪制突變前后交叉小波能量譜與小波相干譜圖,如圖4、圖5所示。

    注:(a)(c)(e)為突變前,(b)(d)(f)為突變后。圖4 PDO、AO、ENSO與SPEI-12交叉小波功率譜Fig.4 Cross wavelet power spectrum between SPEI-12 and PDO,AO,ENSO

    注:(a)(c)(e)為突變前,(b)(d)(f)為突變后。圖5 PDO、AO、ENSO與SPEI-12小波相干譜Fig.5 Wavelet coherent spectrum between SPEI-12 and PDO,AO,ENSO

    由圖4~5可知:松花江流域SPEI-12指數(shù)在高能量區(qū)突變前與PDO無(wú)明顯共振周期(見(jiàn)圖4(a)),突變后1993~1997年存在位相角為90°、尺度為3~4 a的共振周期,表示在此期間高能量區(qū)SPEI-12變化滯后于PDO變化1/4個(gè)共振周期(見(jiàn)圖4(b));在低能量區(qū),突變前與PDO無(wú)明顯共振周期(見(jiàn)圖5(a)),突變后存在位相角為90°、尺度為2~3 a的共振周期(見(jiàn)圖5(b)),即低能量區(qū)SPEI-12變化滯后于PDO變化1/4個(gè)共振周期。

    SPEI-12與AO在高能量區(qū)突變前1966~1971年存在位相角0°、尺度2~4 a的共振周期(見(jiàn)圖4(c)),即SPEI-12與AO負(fù)相關(guān)變化,突變后無(wú)明顯共振周期(見(jiàn)圖4(d));在低能量區(qū)突變前存在位相角約270°、尺度2~5a的共振周期(見(jiàn)圖5(c)),即SPEI-12變化提前于AO變化1/4個(gè)周期,突變后無(wú)顯著共振周期(見(jiàn)圖5(d))。

    SPEI-12與ENSO高能量區(qū)突變前1966~1975年顯示有位相角180°、尺度為1.5~4.0 a的共振周期(見(jiàn)圖4(e)),即兩者同步變化。突變后在1994~2000年存在位相角0°、尺度為2.5~4.0 a的共振周期(見(jiàn)圖4(f)),即兩者負(fù)相關(guān)變化;低能量區(qū)1962~1972年存在尺度為2.5~5.0 a的同步變化共振周期(見(jiàn)圖5(e)),突變后無(wú)顯著共振周期(見(jiàn)圖5(f))。

    3.2.2時(shí)滯相關(guān)分析

    交叉小波與小波相干分析顯示,突變后PDO和ENSO與SPEI-12共振現(xiàn)象更加顯著,因此對(duì)PDO、ENSO與SPEI-12指數(shù)序列分別進(jìn)行突變前后的同期及滯后1~12個(gè)月的Pearson相關(guān)系數(shù)計(jì)算,結(jié)果如表3、表4所列。

    表3 PDO與SPEI-12指數(shù)突變前后時(shí)滯相關(guān)分析Tab.3 Simultaneous and lag correlation analysis between PDO and SPEI-12 index before and after mutation

    表4 ENSO與SPEI-12指數(shù)突變前后時(shí)滯相關(guān)分析Tab.4 Simultaneous and lag correlation analysis between ENSO and SPEI-12 index before and after mutation

    由表3~4可知:SPEI-12突變前與PDO同期與滯后1~3個(gè)月具有逐漸減弱的顯著正相關(guān)性,滯后9~12個(gè)月轉(zhuǎn)變?yōu)轱@著的負(fù)相關(guān)性;突變后同期與滯后1~12個(gè)月具有極顯著的正相關(guān)性,相關(guān)系數(shù)隨滯后尺度增大逐漸減小。SPEI-12突變前與ENSO同期與滯后1~10個(gè)月具有顯著的正相關(guān)性,其中滯后4~5個(gè)月相關(guān)性更為顯著;突變后僅滯后6~8個(gè)月具有較顯著的正相關(guān)性。

    3.3 討 論

    綜合上述,松花江流域SPEI-12在1984年左右突變后,在時(shí)間、空間上皆整體呈現(xiàn)減小趨勢(shì),齊齊哈爾以西北呼倫貝爾以東南地區(qū)、黑河市周邊地區(qū)和吉林北部地區(qū)變化最為顯著。這與吳燕鋒等[38]在分析松花江水文干旱與氣象干旱關(guān)系中氣象干旱在1982年后趨于減輕的演變結(jié)論相符。曾鵬[39]也在其對(duì)中國(guó)未來(lái)干旱趨勢(shì)模擬研究中以RCP4.5與RCP8.5兩種模擬情景模式模擬出松花江流域的干燥趨勢(shì)顯著,SPEI指數(shù)整體呈現(xiàn)負(fù)值,與本文的研究結(jié)果一致。在PDO、AO、ENSO與SPEI-12的交叉小波與小波相干分析中,松花江流域突變前后SPEI-12指數(shù)與遙相關(guān)因子存在共振周期,其中PDO和ENSO與突變后SPEI-12共振相互作用更為顯著,表明遙相關(guān)因子可能是導(dǎo)致松花江流域SPEI-12指數(shù)降低的重要因素之一。

    而在進(jìn)一步的相關(guān)分析中,PDO和ENSO與SPEI-12突變前后相關(guān)性質(zhì)的變化與小波分析部分結(jié)論相符,即PDO和ENSO是導(dǎo)致松花江流域SPEI-12指數(shù)降低的重要的遙相關(guān)影響因素。ENSO表現(xiàn)為赤道中東部太平洋海面溫度異常增高或降低,暖相與冷相皆會(huì)引起沃克環(huán)流的變化,Ren等[40]也在研究中證明了松花江流域SPEI指數(shù)與不同ENSO指數(shù)皆存在顯著的相關(guān)關(guān)系;PDO體現(xiàn)為北太平洋西部和中部海溫的異常。兩者皆起源于太平洋與大氣的相互作用,影響海洋大氣水平-垂直相互作用,給太平洋沿岸大陸乃至全球氣候帶來(lái)熱量及水汽輸送的異常變化。但除了PDO與ENSO的影響外,松花江流域干旱的成因機(jī)制還與季風(fēng)、人類(lèi)活動(dòng)和地形地質(zhì)等有關(guān)。并且由于數(shù)據(jù)的限制,本文選擇使用了Thornthwaite方法計(jì)算潛在蒸散發(fā),該方法只考慮了溫度的影響,忽略了其他如風(fēng)速、輻射等氣候因子的影響,因此松花江旱澇狀況的評(píng)估與遙相關(guān)因子之間的耦合機(jī)理仍需進(jìn)一步研究。

    4 結(jié) 論

    本文結(jié)合松花江流域52個(gè)氣象站點(diǎn)1960~2020年共61a逐日降水與氣溫?cái)?shù)據(jù),計(jì)算12個(gè)月尺度的標(biāo)準(zhǔn)化降水蒸散指數(shù),識(shí)別其時(shí)間序列突變點(diǎn),并運(yùn)用交叉小波變換與小波相干分析,得出以下結(jié)論:

    (1) 松花江流域SPEI-12指數(shù)突變節(jié)點(diǎn)約為1984年,空間上呈現(xiàn)中部及中部偏東地區(qū)突變前輕澇、突變后輕旱的變化特征,流域整體SPEI-12指數(shù)趨于減小,中部及中部偏東地區(qū)變化最為明顯。

    (2) 交叉小波功率譜與小波相干譜顯示,SPEI-12與PDO、AO、ENSO存在共振周期,其中PDO和ENSO突變前后差異更明顯,對(duì)突變后SPEI-12共振相互作用更加顯著。

    (3) 同期與滯后相關(guān)性分析顯示,突變前后SPEI-12與ENSO相關(guān)時(shí)滯尺度有較大變化。SPEI-12與PDO隨時(shí)滯尺度變化突變前相關(guān)性質(zhì)具有較大差異,突變后轉(zhuǎn)變?yōu)闃O顯著的正相關(guān)性,表明PDO與ENSO是導(dǎo)致松花江流域SPEI-12指數(shù)突變后數(shù)值減小的重要的遙相關(guān)影響因素。

    猜你喜歡
    松花江流域共振小波
    自然災(zāi)害、農(nóng)戶(hù)傳統(tǒng)風(fēng)險(xiǎn)管理手段與天氣指數(shù)保險(xiǎn)需求
    構(gòu)造Daubechies小波的一些注記
    基于MATLAB的小波降噪研究
    電子制作(2019年13期)2020-01-14 03:15:32
    吉林省松花江流域產(chǎn)業(yè)生態(tài)系統(tǒng)適應(yīng)性分析
    安然 與時(shí)代同頻共振
    選硬人打硬仗——紫陽(yáng)縣黨建與脫貧同頻共振
    基于一維水動(dòng)力模型的洪水頂托影響分析
    基于改進(jìn)的G-SVS LMS 與冗余提升小波的滾動(dòng)軸承故障診斷
    CTA 中紡院+ 化纖聯(lián)盟 強(qiáng)強(qiáng)聯(lián)合 科技共振
    農(nóng)業(yè)面源氨氮污染減排體系的構(gòu)建及減排控制對(duì)策
    村上凉子中文字幕在线| 色噜噜av男人的天堂激情| 亚洲天堂国产精品一区在线| 性插视频无遮挡在线免费观看| 91久久精品电影网| 啦啦啦韩国在线观看视频| 久久午夜福利片| 日韩人妻高清精品专区| 观看免费一级毛片| 中亚洲国语对白在线视频| 少妇的逼好多水| 国产一区二区三区av在线 | 久久欧美精品欧美久久欧美| 国产色婷婷99| 午夜精品在线福利| 丰满乱子伦码专区| 国产免费一级a男人的天堂| 国产色爽女视频免费观看| 可以在线观看的亚洲视频| 亚洲美女视频黄频| 欧美三级亚洲精品| 淫妇啪啪啪对白视频| 美女 人体艺术 gogo| 少妇的逼好多水| 少妇被粗大猛烈的视频| av天堂中文字幕网| 成人精品一区二区免费| 男人狂女人下面高潮的视频| 毛片一级片免费看久久久久 | 久久人人爽人人爽人人片va| 在线观看av片永久免费下载| 九色国产91popny在线| 桃红色精品国产亚洲av| 亚洲狠狠婷婷综合久久图片| 国产免费av片在线观看野外av| 国产白丝娇喘喷水9色精品| 久久久久久久久久久丰满 | 久久精品国产99精品国产亚洲性色| 夜夜爽天天搞| 久久久久久国产a免费观看| 中文字幕av成人在线电影| 免费av观看视频| 黄色丝袜av网址大全| 乱系列少妇在线播放| 乱人视频在线观看| 精品人妻1区二区| 亚洲av中文字字幕乱码综合| 亚洲欧美激情综合另类| 俄罗斯特黄特色一大片| 国产老妇女一区| 国产高清三级在线| 国产美女午夜福利| 69av精品久久久久久| 我的女老师完整版在线观看| 91久久精品电影网| 日韩,欧美,国产一区二区三区 | 欧美成人性av电影在线观看| 亚洲美女视频黄频| 国内少妇人妻偷人精品xxx网站| 内射极品少妇av片p| 男女那种视频在线观看| 免费在线观看成人毛片| 久久草成人影院| 干丝袜人妻中文字幕| 丝袜美腿在线中文| 美女高潮喷水抽搐中文字幕| 99久久久亚洲精品蜜臀av| 欧美绝顶高潮抽搐喷水| 亚洲一区高清亚洲精品| 国模一区二区三区四区视频| 国产黄a三级三级三级人| 久久人人精品亚洲av| 小说图片视频综合网站| 国产精品一区二区免费欧美| 国产人妻一区二区三区在| 18禁在线播放成人免费| 91麻豆精品激情在线观看国产| 一进一出抽搐gif免费好疼| 久久久久免费精品人妻一区二区| 国产v大片淫在线免费观看| 99在线视频只有这里精品首页| 麻豆国产97在线/欧美| 2021天堂中文幕一二区在线观| 亚洲经典国产精华液单| 欧美3d第一页| 国产91精品成人一区二区三区| 免费在线观看日本一区| 久久久久久大精品| 真人做人爱边吃奶动态| 听说在线观看完整版免费高清| 欧美绝顶高潮抽搐喷水| 欧美成人性av电影在线观看| .国产精品久久| 久久久久九九精品影院| 欧美一区二区精品小视频在线| 又爽又黄无遮挡网站| 久久久久久久久大av| 欧美一级a爱片免费观看看| 天堂√8在线中文| 亚洲天堂国产精品一区在线| 国产精华一区二区三区| 国产一区二区在线av高清观看| 日韩中字成人| 亚洲色图av天堂| 波多野结衣巨乳人妻| 午夜福利成人在线免费观看| 久久午夜亚洲精品久久| 搡老妇女老女人老熟妇| 国产精品一区二区三区四区久久| 美女大奶头视频| 国产三级在线视频| 成年女人看的毛片在线观看| 中文字幕av在线有码专区| 特大巨黑吊av在线直播| 日本欧美国产在线视频| 国产高清有码在线观看视频| 久久人人爽人人爽人人片va| 91av网一区二区| 搡女人真爽免费视频火全软件 | 亚洲国产色片| 琪琪午夜伦伦电影理论片6080| eeuss影院久久| 2021天堂中文幕一二区在线观| 国产三级中文精品| 亚洲av二区三区四区| 亚洲在线自拍视频| 少妇熟女aⅴ在线视频| 少妇被粗大猛烈的视频| 搞女人的毛片| 精品久久久久久成人av| 禁无遮挡网站| 午夜影院日韩av| 白带黄色成豆腐渣| 国产成人a区在线观看| 少妇人妻精品综合一区二区 | 久久午夜亚洲精品久久| 国产主播在线观看一区二区| 免费av不卡在线播放| 国产探花极品一区二区| 免费在线观看日本一区| 欧美成人一区二区免费高清观看| 婷婷六月久久综合丁香| av.在线天堂| 色综合婷婷激情| 国产麻豆成人av免费视频| 俄罗斯特黄特色一大片| 午夜a级毛片| av福利片在线观看| 国产色婷婷99| bbb黄色大片| 亚洲精华国产精华精| 美女黄网站色视频| 91久久精品电影网| 免费人成在线观看视频色| 麻豆精品久久久久久蜜桃| 午夜亚洲福利在线播放| 成人国产麻豆网| 熟妇人妻久久中文字幕3abv| 亚洲久久久久久中文字幕| 一本久久中文字幕| 欧美精品国产亚洲| 精品一区二区三区人妻视频| 欧美人与善性xxx| 日韩欧美一区二区三区在线观看| 1000部很黄的大片| 久久99热6这里只有精品| 欧美成人一区二区免费高清观看| 婷婷丁香在线五月| 琪琪午夜伦伦电影理论片6080| 无遮挡黄片免费观看| 精品一区二区三区视频在线观看免费| 小蜜桃在线观看免费完整版高清| 中文字幕免费在线视频6| 在线观看午夜福利视频| 美女黄网站色视频| 亚洲av中文av极速乱 | 成人精品一区二区免费| 18禁黄网站禁片午夜丰满| 淫秽高清视频在线观看| 久久这里只有精品中国| 久久精品91蜜桃| 国产激情偷乱视频一区二区| 亚洲人与动物交配视频| 午夜亚洲福利在线播放| 国产一区二区在线av高清观看| 97碰自拍视频| 看免费成人av毛片| 又黄又爽又刺激的免费视频.| aaaaa片日本免费| 一区二区三区四区激情视频 | 女同久久另类99精品国产91| 中出人妻视频一区二区| 久久精品影院6| 成人欧美大片| 成年女人永久免费观看视频| 国产精品永久免费网站| 亚洲精品成人久久久久久| 蜜桃亚洲精品一区二区三区| 舔av片在线| 看片在线看免费视频| 免费看美女性在线毛片视频| 免费av不卡在线播放| 很黄的视频免费| 午夜a级毛片| 亚洲最大成人中文| 麻豆一二三区av精品| 亚洲国产精品久久男人天堂| 久久久久国内视频| 99国产精品一区二区蜜桃av| 搡老岳熟女国产| 伊人久久精品亚洲午夜| 中文字幕人妻熟人妻熟丝袜美| 男人的好看免费观看在线视频| 嫁个100分男人电影在线观看| 麻豆成人午夜福利视频| 亚洲性夜色夜夜综合| 亚洲精品色激情综合| 97碰自拍视频| 日本熟妇午夜| 大又大粗又爽又黄少妇毛片口| 成熟少妇高潮喷水视频| 18禁黄网站禁片午夜丰满| 午夜福利视频1000在线观看| 91久久精品国产一区二区三区| 日本与韩国留学比较| 欧美不卡视频在线免费观看| 老女人水多毛片| 精品午夜福利视频在线观看一区| 久久久久九九精品影院| 日韩欧美免费精品| 国产精品爽爽va在线观看网站| 搡老岳熟女国产| 久久精品人妻少妇| 少妇高潮的动态图| av在线老鸭窝| 久久草成人影院| 性色avwww在线观看| 亚洲av日韩精品久久久久久密| 天天一区二区日本电影三级| 能在线免费观看的黄片| avwww免费| 国产亚洲精品av在线| 免费观看的影片在线观看| 国产国拍精品亚洲av在线观看| av福利片在线观看| 久久99热6这里只有精品| 舔av片在线| 成人毛片a级毛片在线播放| 一本精品99久久精品77| 成年人黄色毛片网站| 久久久久久伊人网av| 久久久久久久久久成人| 欧美xxxx性猛交bbbb| 国产精品免费一区二区三区在线| 三级男女做爰猛烈吃奶摸视频| 亚洲av成人精品一区久久| 日本a在线网址| 搡老岳熟女国产| 人妻夜夜爽99麻豆av| 国产爱豆传媒在线观看| 精品一区二区三区视频在线观看免费| 免费av不卡在线播放| 在线a可以看的网站| 精品久久久久久久人妻蜜臀av| 国产成人a区在线观看| 亚洲一区高清亚洲精品| 久久久国产成人免费| 美女xxoo啪啪120秒动态图| 久久这里只有精品中国| 三级国产精品欧美在线观看| 国产高清不卡午夜福利| 内射极品少妇av片p| 国产激情偷乱视频一区二区| 亚洲国产色片| 日本撒尿小便嘘嘘汇集6| 午夜老司机福利剧场| 国产黄色小视频在线观看| 91在线精品国自产拍蜜月| 国产精品一区www在线观看 | 欧美三级亚洲精品| 国产av在哪里看| 国产精品免费一区二区三区在线| 久久国产乱子免费精品| 国产伦一二天堂av在线观看| 97超级碰碰碰精品色视频在线观看| 少妇人妻精品综合一区二区 | 精品久久久久久久久久久久久| 国产亚洲av嫩草精品影院| 亚洲精品在线观看二区| 美女黄网站色视频| 97超视频在线观看视频| 床上黄色一级片| 露出奶头的视频| 亚洲四区av| 亚洲av电影不卡..在线观看| 联通29元200g的流量卡| 免费av不卡在线播放| 国产乱人视频| 男人舔女人下体高潮全视频| 久久这里只有精品中国| 99在线人妻在线中文字幕| 成人午夜高清在线视频| 熟女人妻精品中文字幕| 国产成人aa在线观看| 美女高潮的动态| 特大巨黑吊av在线直播| 精品午夜福利视频在线观看一区| 亚洲精品456在线播放app | 国产激情偷乱视频一区二区| 午夜视频国产福利| 亚洲最大成人中文| 欧美又色又爽又黄视频| 亚洲av美国av| 久99久视频精品免费| 国产色爽女视频免费观看| av专区在线播放| 国产大屁股一区二区在线视频| 99热精品在线国产| 亚洲无线观看免费| 成年女人看的毛片在线观看| 婷婷色综合大香蕉| 精品乱码久久久久久99久播| 国产一区二区在线av高清观看| 久久人人爽人人爽人人片va| 国产高清有码在线观看视频| 内地一区二区视频在线| 久久久久久九九精品二区国产| 亚洲精品影视一区二区三区av| av专区在线播放| av视频在线观看入口| 免费人成视频x8x8入口观看| 99riav亚洲国产免费| 中文字幕精品亚洲无线码一区| 亚洲精品乱码久久久v下载方式| av在线老鸭窝| 成人永久免费在线观看视频| 国产乱人伦免费视频| 成人av一区二区三区在线看| 欧美黑人欧美精品刺激| 男女之事视频高清在线观看| 精品无人区乱码1区二区| 精品人妻一区二区三区麻豆 | 欧美日韩乱码在线| 精品免费久久久久久久清纯| 欧美xxxx黑人xx丫x性爽| 成人性生交大片免费视频hd| 男人狂女人下面高潮的视频| 国产精品久久久久久久久免| 亚洲图色成人| 日韩欧美在线乱码| 亚洲中文字幕日韩| 亚洲美女黄片视频| 国产精品免费一区二区三区在线| 狂野欧美激情性xxxx在线观看| 国产精品久久视频播放| 亚洲人与动物交配视频| 在现免费观看毛片| 欧美在线一区亚洲| 校园人妻丝袜中文字幕| 在线观看美女被高潮喷水网站| 性色avwww在线观看| 美女高潮的动态| 赤兔流量卡办理| 久久久久久久久久成人| 国产极品精品免费视频能看的| 最好的美女福利视频网| 精品一区二区免费观看| 国产一区二区在线av高清观看| 久久久精品欧美日韩精品| 国产精品爽爽va在线观看网站| 亚洲av熟女| 51国产日韩欧美| 亚洲精品乱码久久久v下载方式| 日本熟妇午夜| 又粗又爽又猛毛片免费看| 99热网站在线观看| 狂野欧美激情性xxxx在线观看| 女人十人毛片免费观看3o分钟| 无人区码免费观看不卡| 日本一本二区三区精品| 国产高清不卡午夜福利| 成人三级黄色视频| 麻豆国产97在线/欧美| 国内精品宾馆在线| 一级黄片播放器| 十八禁国产超污无遮挡网站| 午夜久久久久精精品| 国产亚洲精品久久久久久毛片| 婷婷色综合大香蕉| 国产欧美日韩精品亚洲av| 18禁在线播放成人免费| 国产亚洲av嫩草精品影院| 亚洲avbb在线观看| 在线播放无遮挡| 久久久久久九九精品二区国产| 国产伦精品一区二区三区视频9| 国产黄色小视频在线观看| 在现免费观看毛片| 国产美女午夜福利| 中国美白少妇内射xxxbb| 国产爱豆传媒在线观看| 日韩中字成人| 国产伦人伦偷精品视频| 最近最新免费中文字幕在线| 午夜激情欧美在线| 午夜福利成人在线免费观看| 一个人看视频在线观看www免费| 成人av一区二区三区在线看| 我的女老师完整版在线观看| 亚洲欧美清纯卡通| 精品99又大又爽又粗少妇毛片 | 亚洲专区中文字幕在线| 日日撸夜夜添| 国内毛片毛片毛片毛片毛片| 久久久久性生活片| 99久久成人亚洲精品观看| 婷婷亚洲欧美| 精品久久久久久久久亚洲 | 亚洲成a人片在线一区二区| 国产精品一区二区免费欧美| 日韩中文字幕欧美一区二区| 精品人妻1区二区| 麻豆国产97在线/欧美| 好男人在线观看高清免费视频| 成熟少妇高潮喷水视频| 日日摸夜夜添夜夜添av毛片 | 亚洲av日韩精品久久久久久密| 国产精品久久视频播放| 真实男女啪啪啪动态图| ponron亚洲| 88av欧美| 精品久久久久久久人妻蜜臀av| 欧美三级亚洲精品| 国产蜜桃级精品一区二区三区| 久久午夜福利片| 精品久久久久久成人av| 中文字幕av在线有码专区| 亚洲第一区二区三区不卡| 久久久久久久久久黄片| 午夜爱爱视频在线播放| 亚洲三级黄色毛片| 成年女人看的毛片在线观看| 美女免费视频网站| 国产一区二区亚洲精品在线观看| 给我免费播放毛片高清在线观看| 亚洲不卡免费看| 网址你懂的国产日韩在线| 色在线成人网| 日日干狠狠操夜夜爽| 精品人妻熟女av久视频| 听说在线观看完整版免费高清| 亚洲av.av天堂| 日本 欧美在线| 九九爱精品视频在线观看| 99久久九九国产精品国产免费| 免费不卡的大黄色大毛片视频在线观看 | 久久久久久九九精品二区国产| 97超视频在线观看视频| 亚洲精品久久国产高清桃花| 午夜福利成人在线免费观看| 亚洲av二区三区四区| 国产高清激情床上av| 国产精品久久久久久av不卡| 99久国产av精品| 中文资源天堂在线| 少妇丰满av| 女人十人毛片免费观看3o分钟| 91麻豆av在线| 亚洲精品一卡2卡三卡4卡5卡| 床上黄色一级片| 日本 欧美在线| 欧美xxxx性猛交bbbb| 国产精品99久久久久久久久| av在线观看视频网站免费| 久久精品影院6| 免费看美女性在线毛片视频| 内地一区二区视频在线| 窝窝影院91人妻| 波野结衣二区三区在线| 日韩精品有码人妻一区| 日韩av在线大香蕉| 黄色欧美视频在线观看| 欧美日韩黄片免| 高清日韩中文字幕在线| 在线观看午夜福利视频| 午夜福利成人在线免费观看| 成人毛片a级毛片在线播放| 88av欧美| 最近在线观看免费完整版| 亚洲人成网站高清观看| 三级国产精品欧美在线观看| АⅤ资源中文在线天堂| 99热6这里只有精品| 国产精品野战在线观看| 日韩人妻高清精品专区| 久久久精品欧美日韩精品| 性欧美人与动物交配| 韩国av一区二区三区四区| 亚洲欧美日韩无卡精品| 亚洲精品粉嫩美女一区| 麻豆国产av国片精品| 国产精品电影一区二区三区| 欧美精品啪啪一区二区三区| 又粗又爽又猛毛片免费看| 天堂√8在线中文| 可以在线观看的亚洲视频| 亚洲成av人片在线播放无| 我的老师免费观看完整版| 亚洲无线在线观看| x7x7x7水蜜桃| 国产精品乱码一区二三区的特点| 亚洲真实伦在线观看| 狂野欧美激情性xxxx在线观看| 久久精品国产亚洲网站| 91狼人影院| 男女下面进入的视频免费午夜| 春色校园在线视频观看| 日韩一区二区视频免费看| 婷婷精品国产亚洲av在线| 超碰av人人做人人爽久久| 在现免费观看毛片| 日韩在线高清观看一区二区三区 | 亚洲图色成人| 又紧又爽又黄一区二区| 亚洲av免费在线观看| 午夜激情福利司机影院| 两性午夜刺激爽爽歪歪视频在线观看| 国产一区二区三区av在线 | 51国产日韩欧美| 国产精品国产高清国产av| 欧美日韩中文字幕国产精品一区二区三区| 午夜激情欧美在线| 中文资源天堂在线| 18禁黄网站禁片免费观看直播| 岛国在线免费视频观看| 久久久久久久亚洲中文字幕| 亚洲精品日韩av片在线观看| 91久久精品国产一区二区成人| 欧美国产日韩亚洲一区| 永久网站在线| 国产精品精品国产色婷婷| 国产精品久久久久久久久免| 日韩精品青青久久久久久| 日本a在线网址| 精品久久久久久久久久久久久| 中出人妻视频一区二区| 老熟妇乱子伦视频在线观看| 我的女老师完整版在线观看| 51国产日韩欧美| 搡老熟女国产l中国老女人| 午夜a级毛片| 国产单亲对白刺激| 神马国产精品三级电影在线观看| 真人一进一出gif抽搐免费| 婷婷精品国产亚洲av| 搡老妇女老女人老熟妇| 在线国产一区二区在线| 国产一区二区在线av高清观看| 久久久久久久亚洲中文字幕| 久久久久久伊人网av| 久久精品国产亚洲av天美| 国产不卡一卡二| 免费看av在线观看网站| 免费看a级黄色片| 中文字幕精品亚洲无线码一区| 午夜福利高清视频| 亚洲欧美日韩高清专用| 亚洲专区国产一区二区| 九色成人免费人妻av| 精品99又大又爽又粗少妇毛片 | a级毛片免费高清观看在线播放| 成年女人毛片免费观看观看9| 国产熟女欧美一区二区| 91久久精品国产一区二区成人| 啦啦啦观看免费观看视频高清| 五月伊人婷婷丁香| 热99re8久久精品国产| 身体一侧抽搐| 亚洲精品久久国产高清桃花| 中国美白少妇内射xxxbb| 搡女人真爽免费视频火全软件 | 男女啪啪激烈高潮av片| 亚洲欧美日韩卡通动漫| 赤兔流量卡办理| 国产精品日韩av在线免费观看| 国产久久久一区二区三区| 自拍偷自拍亚洲精品老妇| 欧美成人a在线观看| 少妇的逼水好多| 国产伦精品一区二区三区视频9| 别揉我奶头~嗯~啊~动态视频| 亚洲人成网站在线播放欧美日韩| 国产精品女同一区二区软件 | 又黄又爽又免费观看的视频| АⅤ资源中文在线天堂| 狂野欧美白嫩少妇大欣赏| 黄色一级大片看看| 床上黄色一级片| 91麻豆av在线| 成人国产一区最新在线观看| 欧美成人a在线观看| 人人妻人人看人人澡| 他把我摸到了高潮在线观看| 久久午夜亚洲精品久久| 女同久久另类99精品国产91| 听说在线观看完整版免费高清| 男女那种视频在线观看| 国产一区二区在线观看日韩| 97热精品久久久久久| 麻豆成人av在线观看| 免费大片18禁| 真人做人爱边吃奶动态| 亚洲熟妇熟女久久| 色av中文字幕| 亚洲真实伦在线观看|