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

    全國小流域分布式單位線綜合分析

    2019-10-09 08:53:28翟曉燕劉榮華李照會
    關(guān)鍵詞:匯流洪峰流量坡度

    郭 良 ,翟曉燕 ,劉榮華 ,李照會

    (1.中國水利水電科學(xué)研究院,北京 100038;2.水利部防洪抗旱減災(zāi)工程技術(shù)研究中心,北京 100038)

    1 研究背景

    單位線技術(shù)是應(yīng)用水文學(xué)最有效的工具,方法應(yīng)用簡便、計算效率高[1-2],已廣泛用于我國流域水文預(yù)報、水利工程設(shè)計等[3-4]。然而,流域降雨徑流過程實質(zhì)上是一個復(fù)雜的水文非線性系統(tǒng),傳統(tǒng)單位線法一般將其概化為線性時不變系統(tǒng)[5-6],忽略了凈雨空間分布和下墊面空間異質(zhì)性的影響,此外,單位線匯流參數(shù)的確定大多依賴于實測暴雨洪水資料,限制了其在缺資料地區(qū)中小流域的應(yīng)用。流域地形地貌特征與降雨特性是影響小流域暴雨洪水過程的主要因子[7],Maidment等[8]和郭良等[9-11]基于數(shù)字高程模型探索流域水文響應(yīng)與地形地貌、雨強間的非線性定量關(guān)系,提出了分布式單位線法,認(rèn)為流域內(nèi)各點匯流時間的概率密度分布等價于瞬時單位線,實現(xiàn)了缺資料地區(qū)匯流計算,為缺資料地區(qū)山洪預(yù)報預(yù)警提供了技術(shù)支撐。

    我國中小流域暴雨洪水暴漲暴落、突發(fā)性強、歷時短,洪水預(yù)警期較短,開展分布式單位線綜合分析,有助于實現(xiàn)缺資料地區(qū)山洪的早期快速預(yù)報預(yù)警。近年來,“3S”技術(shù)的發(fā)展促進了流域地形地貌特征的高效準(zhǔn)確獲取,為辨識流域?qū)傩詫π×饔騾R流過程的定量影響提供了數(shù)據(jù)基礎(chǔ)[12]。流域地形地貌屬性一般可分為線性尺度測量指標(biāo)和無因次特性指標(biāo)兩類[13],常用于缺資料地區(qū)流域單位線綜合分析和匯流計算。已有研究中采用的流域?qū)傩灾笜?biāo)多為流域面積、溝道長度、流域坡降等參數(shù),未考慮流域內(nèi)部地形地貌空間分布異質(zhì)性對匯流過程的影響,導(dǎo)致設(shè)計洪水估算不合理、洪水預(yù)報預(yù)警精準(zhǔn)度低等問題[14-16]。因此,需要進一步挖掘能定量描述小流域地形地貌空間差異的屬性參數(shù),構(gòu)建具有匯流成因的單位線綜合公式,提高缺資料地區(qū)中小流域洪水預(yù)報預(yù)警的精準(zhǔn)度和時效性。本研究期望為缺資料地區(qū)中小流域暴雨洪水分析計算提供技術(shù)參考。

    2 全國小流域數(shù)據(jù)集

    2010年以來,中國水利水電科學(xué)研究院結(jié)合全國山洪災(zāi)害防治項目,基于1∶5萬數(shù)字高程模型、2.5 m分辨率數(shù)字正影像圖、1∶5萬數(shù)字線劃圖等高精度國家基礎(chǔ)地理空間信息,構(gòu)建了大范圍小流域精細(xì)劃分和屬性分析技術(shù)體系,將全國劃分為53萬個小流域(平均面積16 km2)和378萬條溝道河流,形成了全國小流域數(shù)據(jù)集[9],覆蓋國土面積868.67萬km2,為全國山洪災(zāi)害防御提供了海量基礎(chǔ)數(shù)據(jù),為不同自然地理格局下暴雨山洪關(guān)鍵因子辨識量化提供依據(jù),實現(xiàn)缺資料地區(qū)洪水預(yù)報預(yù)警。全國小流域數(shù)據(jù)集主要包括53萬個小流域的75項主要基礎(chǔ)屬性和分布式單位線組等。

    2.1 小流域基礎(chǔ)屬性小流域基礎(chǔ)屬性包括反映小流域形狀、坡度、高程、匯流路徑長度等常用屬性,如小流域面積、周長、形狀系數(shù)、平均坡度、最長匯流路徑長度及其比降、溝河長度及平均比降等。通過分析小流域內(nèi)部高程、匯流路徑長度與集水面積間的關(guān)系,提出了小流域不均勻系數(shù)、加權(quán)平均坡度等反映小流域匯流非均質(zhì)性的參數(shù)。

    小流域平均匯流路徑長度為小流域內(nèi)水滴由各網(wǎng)格點匯流至小流域出口的匯流路徑平均值。以小流域平均匯流路徑長度與最長匯流路徑長度一半的比值,反映小流域內(nèi)水流流路分布的非均質(zhì)性,即小流域不均勻系數(shù),計算公式如下:

    式中:L為小流域平均匯流路徑長度;Lm為小流域最長匯流路徑長度;C為小流域不均勻系數(shù)。

    對于小流域坡面而言,以各網(wǎng)格點上游集水面積占所有網(wǎng)格點集水面積之和的比例作為權(quán)重,基于各點坡度計算小流域加權(quán)平均坡度。加權(quán)平均坡度通過小流域內(nèi)各點集水面積的空間分布反映小流域內(nèi)地形分布對水流匯集的影響,平均坡度反映了小流域地形的平均起伏狀態(tài)對水流匯集的影響。

    式中:S′為小流域加權(quán)平均坡度;Sj為小流域坡面網(wǎng)格點j的坡度;ij為坡面網(wǎng)格點j的集水面積;n為小流域坡面網(wǎng)格數(shù)。

    全國小流域不均勻系數(shù)和加權(quán)平均坡度頻數(shù)分布如圖1所示。全國26.4%的小流域內(nèi)各點匯流路徑分布趨于均勻分布,小流域流路質(zhì)心與形心基本重合,小流域形狀呈上下游對稱(|C-1|≤0.05);73.6%的小流域流路分布較不均勻,其中,9.7%的小流域流路質(zhì)心偏下游,水流匯集路徑較短,較易形成山洪(C<0.95),63.9%的小流域流路質(zhì)心偏上游,水流匯集路徑較長,洪水集中偏慢(C>1.05)。以加權(quán)平均坡度與平均坡度之差進一步反映小流域坡度的非均質(zhì)分布特性。全國約半數(shù)小流域坡面上、下游坡度分布較為一致(|S′-S|≤0.05);約34.5%的小流域坡面上游地勢陡峭、下游地勢平緩,小流域坡面?zhèn)绕拭娑喑省鞍肌毙头植迹⊿′-S<-0.05),坡面的侵蝕強度較大,易于形成山洪滑坡地質(zhì)災(zāi)害;其余小流域側(cè)剖面多呈“凸”型分布(S′-S>0.05),坡面有利于排水,但水流沖刷作用也較強,容易造成水土流失問題。

    圖1 全國小流域不均勻系數(shù)及加權(quán)平均坡度頻數(shù)分布

    2.2 分布式單位線應(yīng)用1∶5萬DEM和2.5 m土地利用和植被類型信息數(shù)據(jù),充分考慮流域下墊面空間分布異質(zhì)性,提出基于DEM網(wǎng)格、考慮雨強影響匯流非線性特征的分布式單位線方法。在全國范圍內(nèi)共提取了53萬個小流域15組不同雨強(5~100 mm/h)、不同時段(10~60 min)的分布式單位線組,使用不同地貌類型區(qū)361個流域1萬多場暴雨洪水資料進行檢驗,93.4%的流域雨洪模擬精度較高(即洪峰流量相對誤差不超過20%,峰現(xiàn)時間誤差不超過2 h),并已用于全國不同類型區(qū)中小流域洪水模擬及預(yù)報預(yù)警[9,17-19]。

    坡面流流速計算采用改進的SCS流速公式,綜合考慮雨強、坡面地形對流域匯流非線性的影響。基于D8算法,統(tǒng)計流域各網(wǎng)格(25 m×25 m)的徑流匯集路徑及時間,如公式(3)所示。統(tǒng)計流域內(nèi)各網(wǎng)格的匯流時間場,得到流域匯流時間的概率密度分布和累計S曲線,進一步考慮流域的調(diào)蓄作用,即可得到小流域非線性分布式匯流單位線。

    式中:V為水流速度;K為流速系數(shù),主要反映土地利用特征對流速摩阻影響的經(jīng)驗參數(shù);S為流域內(nèi)某網(wǎng)格沿著水流方向的坡降;i為無因次雨強;Tj為第j個網(wǎng)格的匯流時間;Lk為第k個網(wǎng)格的流路長度;c為系數(shù),c=1或為第j個網(wǎng)格匯流路徑上網(wǎng)格的數(shù)量。

    以降雨量30 mm、時段長10 min的單位線為例(圖2),除去湖泊、沼澤和水庫外,全國小流域單位洪峰模數(shù)處于0.01~0.51 m3/(s·km2)之間,平均值為0.26 m3/(s·km2),四分位數(shù)為0.18 m3/(s·km2)和0.33 m3/(s·km2);全國小流域匯流時間處于0.5~2.0 h之間,平均值為1.09 h,四分位數(shù)為0.83 h和1.33 h。各小流域差異明顯,成果直觀反映了不同地區(qū)小流域下墊面條件對洪水集中度和洪水匯集時間的影響。

    圖2 全國小流域單位線洪峰模數(shù)、匯流時間頻數(shù)分布

    2.3 典型流域在全國選取30個典型中小流域共4610個小流域作為研究區(qū),其基本信息如表1所示。典型流域面積為755~5874 km2,平均流域面積為2275 km2,分別包括44~406個小流域,研究區(qū)跨越三大階梯,在高原、山地、丘陵、盆地和平原均有分布,地形地貌差異較大,平均高程為128~5017 m,相對高差為528~3327 m,平均坡度為0.05~0.65;研究區(qū)分布在19個水文亞區(qū),涉及豐平枯水帶及多種水文情勢特征;水系發(fā)育程度不一,對徑流的調(diào)蓄能力不一,河網(wǎng)密度為0.25~0.41 km/km2,河網(wǎng)頻度為0.02~0.03條/km2,發(fā)育系數(shù)為2.32~21.16,水系不均勻系數(shù)為0.05~2.99。典型流域反映了我國自然地理空間特征的復(fù)雜性及差異性,部分典型流域的地形、水系和小流域分布如圖3所示。

    圖3 部分典型流域地形、水系及小流域分布

    表1 典型流域基本信息統(tǒng)計

    3 典型流域單位線公式優(yōu)選方法

    采用Pearson相關(guān)系數(shù)(r)分析流域?qū)傩耘c單位線特征值(即洪峰流量、匯流時間和峰現(xiàn)時間)之間的線性相關(guān)性,以|r|在0.2~0.5之間為中度相關(guān),以|r|≥0.5為強相關(guān)。對相關(guān)系數(shù)進行顯著性檢驗,原假設(shè)(H0)為:r=0,顯著性檢驗的統(tǒng)計值服從自由度為n-2的T分布。檢驗統(tǒng)計量計算公式為:

    式中:t為統(tǒng)計量;r為相關(guān)系數(shù);n為樣本數(shù);檢驗的顯著性水平α取為0.01,若計算統(tǒng)計量的顯著性水平p<0.01,則拒絕原假設(shè),說明流域?qū)傩耘c單位線特征值間存在顯著的線性相關(guān)性。

    我國中小流域產(chǎn)匯流非線性特質(zhì)十分顯著,采用基于最小二乘法的多元回歸非線性模型識別分布式單位線特征值與小流域?qū)傩蚤g的定量關(guān)系?;貧w模型的一般形式如下:

    式中:X為小流域單位線的特征值,如單位線洪峰流量、匯流時間和峰現(xiàn)時間;M為小流域?qū)傩?,如小流域面積、坡度、匯流路徑長度等;i為所選小流域?qū)傩詡€數(shù);k和n分別為擬合系數(shù)和擬合指數(shù)。

    優(yōu)選最佳擬合公式型的目標(biāo)函數(shù)為均方根誤差(RMSE)最小,采用合格率(ST)和相關(guān)系數(shù)(r)進一步評估擬合公式的優(yōu)劣。RMSE、ST和r的最優(yōu)值分別為0、100%和1,各評估指標(biāo)計算公式如下:

    式中:RMSE為均方根誤差;ST為合格率;r為相關(guān)系數(shù);X為分布式單位線的特征值,如洪峰流量Qm,峰現(xiàn)時間Tp和匯流時間T;N為區(qū)域或流域內(nèi)分布式單位線的總數(shù);n為擬合的單位線特征值合格的分布式單位線數(shù),以單位線特征值的相對誤差在±20%以內(nèi)為合格,即;Xo,i為區(qū)域或流域內(nèi)第i個分布式單位線的特征值;Xs,i為區(qū)域或流域內(nèi)第i個分布式單位線的擬合特征值;Xˉo和Xˉs分別為區(qū)域或流域內(nèi)分布式單位線特征值和擬合特征值的均值。

    4 相關(guān)性分析及結(jié)果評估

    4.1 相關(guān)性分析小流域匯流過程受降雨特性影響較大,不同降雨強度下單位線過程變化顯著,本文主要以雨量為30 mm、時段長為10 min的分布式單位線為例進行單位線綜合分析。結(jié)合小流域匯流特性分析,選取的小流域?qū)傩园ㄐ×饔蛎娣eA、最長匯流路徑長度Lm及其比降J、平均匯流路徑長度L、加權(quán)平均坡度S′、平均坡度S、河段比降Jr等。部分典型流域小流域?qū)傩耘c單位線特征值間的相關(guān)系數(shù)絕對值分布如圖4所示。對于小流域單位線洪峰流量而言,其與所有典型流域小流域面積的正相關(guān)性均較強(0.70≤r≤0.95,p<0.01),與所有典型流域小流域平均/最長匯流路徑長度均存在中度~強度正相關(guān)關(guān)系(0.22≤r≤0.86,p<0.01),與37%~73%的典型流域小流域最長匯流路徑比降、加權(quán)平均比降、平均比降、河段比降存在中度正相關(guān)關(guān)系(0.20≤r<0.50,p<0.01)。對于小流域單位線匯流時間而言,其與所有典型流域小流域平均/最長匯流路徑長度的正相關(guān)性均較強(0.81≤r≤0.98,p<0.01),與所有典型流域小流域最長匯流路徑比降均存在中度負(fù)相關(guān)關(guān)系(-0.81≤r≤-0.23,p<0.01),與17%~25%的典型流域小流域加權(quán)平均比降、平均比降、河段比降存在中度負(fù)相關(guān)關(guān)系(-0.50<r≤-0.20,p<0.01),此外,小流域單位線匯流時間與洪峰流量之間存在中度正相關(guān)關(guān)系(0.20≤r≤0.64,p<0.01)。小流域單位線峰現(xiàn)時間與小流域?qū)傩蚤g的相關(guān)關(guān)系與小流域單位線匯流時間類似。

    圖4 部分典型流域相關(guān)系數(shù)分布

    所選小流域?qū)傩灾饕碚髁诵×饔虻拿娣e、比降、匯流路徑長度等特性,小流域?qū)傩耘c分布式單位線特征值之間存在顯著的相關(guān)關(guān)系,具有匯流成因概念,但不同典型流域的部分屬性與特征值間的相關(guān)性不一,反映了流域下墊面的空間異質(zhì)性對洪水集中度和匯集時間的影響。小流域面積表征了小流域集水區(qū)域的大小,面積越大,形成的洪水量級越大,相應(yīng)的洪峰流量也越大;小流域匯流路徑越長,水流由坡面和主河道匯集至小流域出口的時間越長,此外,小流域匯流路徑長度與小流域面積間存在顯著的正相關(guān)關(guān)系(r=0.82,p<0.01),小流域單位面積匯流路徑越大,小流域?qū)樗恼{(diào)蓄作用越大,洪峰流量的坦化現(xiàn)象越明顯;河段比降、平均坡度、最長匯流路徑比降等特征參數(shù)表征了小流域地形的起伏狀態(tài),地形越陡峭,水流流速越快,對水流的匯集作用越強,在同一時段內(nèi)匯集至小流域出口的水流越多、時間越短。對于小流域內(nèi)部的水流摩阻因素(如糙率)等,暫不考慮,以經(jīng)驗關(guān)系中的綜合系數(shù)來表征,可在后續(xù)研究中進一步細(xì)化該因素的影響。

    綜合已有相關(guān)研究成果[20-21]、小流域匯流特性分析和上述相關(guān)性分析,提出單位線特征值的擬合公式型,如表2、表3所示。

    表2 單位線洪峰流量擬定公式型

    表3 單位線匯流時間和峰現(xiàn)時間擬定公式型

    4.2 洪峰流量擬合評估不同公式型擬合效果不一樣,同一公式型在不同流域擬合效果也有差異。各公式型的評估指標(biāo)如圖5所示??傮w而言,所有公式型的擬合效果均較好,均方根誤差為0.14~0.86 m3/s,合格率為73%~100%,相關(guān)系數(shù)為0.85~0.99。12個典型流域的最優(yōu)公式型為擬定公式型4,即長干河、富爾河、辰清河、岷江、龍河、綿河、五龍河、遂川江、西溪河、布柳河、琴江和新安江,主要分布在橫斷山脈、云貴高原以及太行山沿線等地貌復(fù)雜、地形多變的山區(qū),相應(yīng)的均方根誤差為0.21~0.46 m3/s,合格率為90%~98%,相關(guān)系數(shù)為0.93~0.98;其余18個典型流域的最優(yōu)公式型為擬定公式型7,即貢曲、佳蘆河、小河、湯河、馬灘河、亮江、通道河、秋香江、鴨江、昌曲、俄科河、清水河、尼木瑪曲、九龍溪、益曲、比江、官渡河和伊河,相應(yīng)的均方根誤差為0.14~0.75 m3/s,合格率為86%~100%,相關(guān)系數(shù)為0.93~0.99。單位線洪峰流量與流域面積、最長匯流路徑長度、河段比降、流域(加權(quán))平均坡度關(guān)系密切。小流域源短流急,將已有公式型中的最長匯流路徑比降指標(biāo)細(xì)化為坡面和河段兩部分,精細(xì)地反映了流域地形起伏狀態(tài)對洪水集中度的影響,在一定程度上提高了擬合精度。

    圖5 單位線洪峰流量不同公式型評估指標(biāo)

    相比于已有公式型,最優(yōu)公式型的均方根誤差減少了0.01~0.08 m3/s,合格率基本一致,相關(guān)系數(shù)提高了0~0.03;相比于其他擬定公式型,最優(yōu)公式型的均方根誤差減少了0~0.28 m3/s,合格率提高了0~20%,相關(guān)系數(shù)提高了0~0.09。圖6為部分典型流域最優(yōu)公式型的單位線洪峰流量擬合散點分布,散點均圍繞在1∶1線周圍,且大部分小流域散點落在合格范圍內(nèi),各流域單位線洪峰流量與擬合值吻合較好,構(gòu)建的公式型和所選小流域?qū)傩阅茌^好地反映不同地形地貌類型區(qū)洪水集中度的差異性。

    圖6 部分典型流域單位線洪峰流量最優(yōu)公式型擬合散點圖

    4.3 匯流時間和峰現(xiàn)時間擬合評估單位線匯流時間和峰現(xiàn)時間各公式型的評估指標(biāo)如圖7所示。對于單位線峰現(xiàn)時間而言,所有公式型的均方根誤差為3.15~8.46 min,合格率為61%~100%,相關(guān)系數(shù)為0.82~0.98。所有流域的最優(yōu)擬合公式型均為擬定公式型3,相應(yīng)的均方根誤差為3.15~7.33 min,合格率為65%~100%,相關(guān)系數(shù)為0.83~0.98,單位線峰現(xiàn)時間與流域最長匯流路徑長度及其比降、平均匯流路徑長度關(guān)系密切。相比于其他擬定公式型,所有典型流域最優(yōu)公式型的均方根誤差減少了0~1.91 min,合格率變化了-6%~9%,相關(guān)系數(shù)提高了0~0.03。

    圖7 單位線匯流時間和峰現(xiàn)時間不同公式型評估指標(biāo)

    對于單位線匯流時間而言,所有公式型擬合的均方根誤差為3.16~9.91 min,合格率為79%~100%,相關(guān)系數(shù)為0.81~0.99。所有流域的最優(yōu)擬合公式型為擬定公式型4,相應(yīng)的均方根誤差為3.16~7.21 min,合格率為85%~100%,相關(guān)系數(shù)為0.89~0.99,單位線匯流時間與流域最長匯流路徑長度及其比降、單位線洪峰流量關(guān)系密切。相比于其他擬定公式型,所有典型流域最優(yōu)公式型的均方根誤差減少了0.02~3.29 min,合格率提高了0~16%,相關(guān)系數(shù)提高了0~0.09。圖8為部分典型流域最優(yōu)公式型的單位線匯流時間擬合散點分布,散點均圍繞在1∶1線周圍,且大部分小流域散點落在合格范圍內(nèi),各流域單位線匯流時間與擬合值吻合較好,構(gòu)建的公式型和所選小流域?qū)傩阅茌^好地反映不同地形地貌類型區(qū)的洪水匯集時間的差異性。總體而言,采用反映小流域匯流非均質(zhì)分布特性的流域地形地貌屬性進行分布式單位線綜合分析,具有匯流成因概念,避免了傳統(tǒng)綜合單位線法普遍存在的主觀性,擬合精度高,適于缺資料地區(qū)中小流域暴雨洪水計算[22]。

    圖8 部分典型流域單位線匯流時間最優(yōu)公式型擬合散點圖

    5 結(jié)論

    本文以全國30個不同地形地貌類型區(qū)的中小流域為典型流域,基于全國小流域數(shù)據(jù)集,采用相關(guān)性分析和基于最小二乘法的多元回歸非線性模型進行小流域分布式單位線的綜合分析。主要結(jié)論為:(1)小流域不均勻系數(shù)、加權(quán)平均坡度較好地反映了小流域匯流的非均質(zhì)性,可用于分布式單位線的綜合分析。全國73.6%的小流域流路分布較不均勻,其中9.7%的小流域流路質(zhì)心偏下游,水流匯集路徑較短,較易形成山洪;34.5%的小流域坡面?zhèn)绕拭娉省鞍肌毙头植?,侵蝕強度較大,易于形成山洪滑坡地質(zhì)災(zāi)害,9.2%的小流域側(cè)剖面呈“凸”型分布,坡面有利于排水,但水流沖刷作用也較強,容易造成水土流失。(2)30個典型流域小流域單位線洪峰流量與小流域面積、平均/最長匯流路徑長度存在中度~強度正相關(guān)關(guān)系,與單位面積匯流路徑長度存在負(fù)相關(guān)關(guān)系;單位線匯流時間、峰現(xiàn)時間與小流域平均/最長匯流路徑長度存在較強的正相關(guān)關(guān)系,與最長匯流路徑比降存在中度負(fù)相關(guān)關(guān)系,單位線匯流時間與洪峰流量存在中度正相關(guān)關(guān)系;其余屬性與單位線特征值間的相關(guān)程度不一。(3)12個典型流域小流域單位線洪峰流量擬合的最優(yōu)公式型為即長干河、富爾河、辰清河、岷江、龍河、綿河、五龍河、遂川江、西溪河、布柳河、琴江和新安江,其余18個典型流域的最優(yōu)公式型為即貢曲、佳蘆河、小河、湯河、馬灘河、亮江、通道河、秋香江、鴨江、昌曲、俄科河、清水河、尼木瑪曲、九龍溪、益曲、比江、官渡河和伊河;匯流時間和峰現(xiàn)時間的最優(yōu)公式型分別為和洪峰流量和匯流時間的合格率均在85%以上,相關(guān)系數(shù)在0.90以上,峰現(xiàn)時間的合格率和相關(guān)系數(shù)分別大于65%和0.80。(4)本文主要分析了降雨量為30 mm、時段長為10 min的分布式單位線,今后應(yīng)進一步開展不同雨強、時段長的分布式單位線的綜合分析,繼續(xù)開展小流域內(nèi)部水流摩阻因素對小流域匯流特性的影響量化研究,以進一步提高單位線綜合分析的精度。

    猜你喜歡
    匯流洪峰流量坡度
    關(guān)于公路超高漸變段合成坡度解析與應(yīng)用
    退耕還林工程對渭河洪峰流量的影響
    流域匯流計算方法在風(fēng)峪沙河段的應(yīng)用探討
    佛岡縣潖江流域年洪峰流量P-Ⅲ分布參數(shù)估算
    大南川流域設(shè)計洪峰流量計算分析
    某特小流域設(shè)計洪峰流量計算分析
    基于圖像處理的定位器坡度計算
    電氣化鐵道(2016年4期)2016-04-16 05:59:46
    基于ARM CortexM0的智能光伏匯流采集裝置設(shè)計及應(yīng)用
    一種球載雷達匯流環(huán)設(shè)計
    雷達與對抗(2015年3期)2015-12-09 02:39:00
    坡度在巖石風(fēng)化層解譯中的應(yīng)用
    河北遙感(2015年2期)2015-07-18 11:11:14
    建设人人有责人人尽责人人享有的 | 网址你懂的国产日韩在线| .国产精品久久| 亚洲精品日韩在线中文字幕| 日韩制服骚丝袜av| 最近中文字幕高清免费大全6| 简卡轻食公司| 欧美97在线视频| h日本视频在线播放| 亚洲精品久久久久久婷婷小说| 亚洲精品乱久久久久久| 六月丁香七月| 久久久久性生活片| 国产探花在线观看一区二区| 大话2 男鬼变身卡| 又大又黄又爽视频免费| 色播亚洲综合网| 久99久视频精品免费| 成人性生交大片免费视频hd| 婷婷色综合www| 国产精品日韩av在线免费观看| .国产精品久久| 亚洲精品乱码久久久v下载方式| 日日干狠狠操夜夜爽| 3wmmmm亚洲av在线观看| 色综合色国产| 亚洲国产精品国产精品| 精品久久久精品久久久| 国产精品国产三级国产av玫瑰| 亚洲欧美一区二区三区黑人 | 啦啦啦韩国在线观看视频| 在线观看人妻少妇| 成人欧美大片| 亚洲欧美成人精品一区二区| 黄色配什么色好看| 国产老妇女一区| 伦理电影大哥的女人| 国产伦一二天堂av在线观看| 国产精品综合久久久久久久免费| 国产成人福利小说| 人妻少妇偷人精品九色| 国产伦精品一区二区三区四那| 大陆偷拍与自拍| 亚洲成人中文字幕在线播放| 国产黄片视频在线免费观看| 婷婷色av中文字幕| 一个人看视频在线观看www免费| 日韩一本色道免费dvd| 亚州av有码| 97热精品久久久久久| 99久国产av精品| 久久久精品免费免费高清| 亚洲国产高清在线一区二区三| 26uuu在线亚洲综合色| 美女xxoo啪啪120秒动态图| 国产精品一区www在线观看| 国产 一区 欧美 日韩| 成年版毛片免费区| 精品一区二区免费观看| 一级毛片我不卡| 熟妇人妻不卡中文字幕| 国产精品.久久久| 美女大奶头视频| 国产午夜福利久久久久久| 国产一区二区三区av在线| 欧美另类一区| 国产黄色视频一区二区在线观看| 小蜜桃在线观看免费完整版高清| 国产午夜福利久久久久久| 亚洲综合色惰| 国产亚洲精品久久久com| 日日啪夜夜爽| 亚洲欧美精品自产自拍| 在线免费观看的www视频| 一级毛片我不卡| 亚洲精品一区蜜桃| 日本免费在线观看一区| 亚洲av成人av| 成年人午夜在线观看视频 | 久久国产乱子免费精品| 免费观看无遮挡的男女| 国产在线一区二区三区精| 九草在线视频观看| 久久久成人免费电影| 在线观看免费高清a一片| 晚上一个人看的免费电影| 麻豆成人av视频| freevideosex欧美| 日本爱情动作片www.在线观看| 日本黄色片子视频| 看免费成人av毛片| 菩萨蛮人人尽说江南好唐韦庄| 日韩av在线大香蕉| 中文字幕人妻熟人妻熟丝袜美| h日本视频在线播放| 亚洲av成人精品一二三区| 丝袜美腿在线中文| 亚洲美女搞黄在线观看| 国产精品久久久久久久电影| 国产成人精品婷婷| 午夜激情福利司机影院| 日日摸夜夜添夜夜添av毛片| 美女主播在线视频| 日日干狠狠操夜夜爽| 国产v大片淫在线免费观看| 欧美激情国产日韩精品一区| 热99在线观看视频| 18禁裸乳无遮挡免费网站照片| 伊人久久国产一区二区| 亚洲精品aⅴ在线观看| 国产成人福利小说| 最近的中文字幕免费完整| 欧美变态另类bdsm刘玥| 亚洲乱码一区二区免费版| 舔av片在线| 99热这里只有是精品在线观看| 免费看不卡的av| 婷婷色麻豆天堂久久| or卡值多少钱| 久久久久久久大尺度免费视频| 少妇熟女aⅴ在线视频| 亚洲成人中文字幕在线播放| 91精品伊人久久大香线蕉| 国产单亲对白刺激| 成人欧美大片| 最近中文字幕2019免费版| 国产女主播在线喷水免费视频网站 | 亚洲av免费在线观看| 淫秽高清视频在线观看| 午夜免费激情av| 久久精品夜色国产| 晚上一个人看的免费电影| 在线观看人妻少妇| 国产亚洲精品av在线| 不卡视频在线观看欧美| 成人欧美大片| 免费看av在线观看网站| 禁无遮挡网站| 久久精品久久久久久久性| 99久久人妻综合| 成人午夜精彩视频在线观看| 插逼视频在线观看| 久久久午夜欧美精品| 久久综合国产亚洲精品| 亚洲av电影不卡..在线观看| 亚洲第一区二区三区不卡| 久久久国产一区二区| 非洲黑人性xxxx精品又粗又长| 中国美白少妇内射xxxbb| 国产 一区精品| 亚洲国产精品成人综合色| 亚洲精品一区蜜桃| 亚洲精品自拍成人| 性插视频无遮挡在线免费观看| 中文字幕人妻熟人妻熟丝袜美| 身体一侧抽搐| 久久综合国产亚洲精品| 观看免费一级毛片| 亚洲乱码一区二区免费版| 日日摸夜夜添夜夜爱| 中文字幕人妻熟人妻熟丝袜美| 国产免费一级a男人的天堂| 国国产精品蜜臀av免费| 午夜福利成人在线免费观看| 国产熟女欧美一区二区| 97在线视频观看| 久久精品久久久久久久性| 少妇高潮的动态图| 国产一区亚洲一区在线观看| 欧美3d第一页| 亚洲色图av天堂| 黑人高潮一二区| 日韩 亚洲 欧美在线| 人妻一区二区av| 日韩精品有码人妻一区| 国产精品一区www在线观看| 久久久久免费精品人妻一区二区| 亚洲国产精品专区欧美| 精品一区二区三区视频在线| 色综合色国产| 激情五月婷婷亚洲| 久久久久久久大尺度免费视频| 国产有黄有色有爽视频| 精品久久久久久久久av| 国产精品久久久久久久久免| 国产一级毛片七仙女欲春2| 激情 狠狠 欧美| 亚洲精品中文字幕在线视频 | 国产高清国产精品国产三级 | 精品一区二区三区视频在线| 精品酒店卫生间| 丰满人妻一区二区三区视频av| 水蜜桃什么品种好| 三级男女做爰猛烈吃奶摸视频| 国产av不卡久久| 一级毛片我不卡| 亚洲丝袜综合中文字幕| 久久久久精品久久久久真实原创| 欧美xxxx性猛交bbbb| 天天一区二区日本电影三级| 国产成年人精品一区二区| 一级黄片播放器| 欧美激情久久久久久爽电影| 麻豆国产97在线/欧美| 欧美3d第一页| 日本免费在线观看一区| 国产高清有码在线观看视频| 国产伦精品一区二区三区四那| 五月伊人婷婷丁香| 欧美成人a在线观看| 联通29元200g的流量卡| 成人国产麻豆网| 久久久久精品性色| 欧美最新免费一区二区三区| 91久久精品国产一区二区三区| 亚洲av电影不卡..在线观看| 神马国产精品三级电影在线观看| 亚洲va在线va天堂va国产| 欧美高清成人免费视频www| 婷婷色综合大香蕉| 三级毛片av免费| 人人妻人人澡人人爽人人夜夜 | 天堂影院成人在线观看| 亚洲精品乱码久久久久久按摩| 两个人视频免费观看高清| 一级爰片在线观看| 3wmmmm亚洲av在线观看| 在线免费观看不下载黄p国产| 人妻制服诱惑在线中文字幕| 色综合站精品国产| 日本三级黄在线观看| 久久久久久九九精品二区国产| 成人午夜高清在线视频| kizo精华| 国产男人的电影天堂91| 熟妇人妻不卡中文字幕| 街头女战士在线观看网站| 一区二区三区乱码不卡18| 亚洲精品一二三| 久久久成人免费电影| 一级黄片播放器| 成人亚洲欧美一区二区av| 欧美另类一区| 亚洲在线自拍视频| 久久精品久久久久久噜噜老黄| 好男人在线观看高清免费视频| 久久久久久久久久成人| 久久人人爽人人爽人人片va| 极品教师在线视频| 久久久午夜欧美精品| 午夜老司机福利剧场| 久久亚洲国产成人精品v| 免费av毛片视频| 一级黄片播放器| 精品国产露脸久久av麻豆 | 日韩视频在线欧美| 美女主播在线视频| a级毛色黄片| 国产av码专区亚洲av| 乱人视频在线观看| 久久久久久久亚洲中文字幕| 国内揄拍国产精品人妻在线| 少妇熟女欧美另类| 亚洲第一区二区三区不卡| 免费观看性生交大片5| 可以在线观看毛片的网站| 男女下面进入的视频免费午夜| 亚洲国产精品成人久久小说| 欧美日韩一区二区视频在线观看视频在线 | 国产精品av视频在线免费观看| 国产美女午夜福利| 国产熟女欧美一区二区| 亚洲精品影视一区二区三区av| 波多野结衣巨乳人妻| 精品久久久久久电影网| 如何舔出高潮| 2021天堂中文幕一二区在线观| 日本一二三区视频观看| 久久久久久久久久久免费av| 日本免费a在线| av黄色大香蕉| 久久久久久久久久人人人人人人| av在线观看视频网站免费| 丰满乱子伦码专区| 又黄又爽又刺激的免费视频.| 亚洲自拍偷在线| 精品欧美国产一区二区三| 国内精品一区二区在线观看| 日韩强制内射视频| 又黄又爽又刺激的免费视频.| 黄片wwwwww| 非洲黑人性xxxx精品又粗又长| 免费观看在线日韩| 亚洲婷婷狠狠爱综合网| 女人久久www免费人成看片| 狂野欧美激情性xxxx在线观看| 在线观看人妻少妇| 免费av毛片视频| 91精品国产九色| 午夜日本视频在线| 国产黄色视频一区二区在线观看| 两个人的视频大全免费| 国产亚洲精品久久久com| 伦理电影大哥的女人| 国产精品精品国产色婷婷| 欧美97在线视频| www.色视频.com| 国产一级毛片在线| 亚洲精品乱码久久久久久按摩| 亚洲最大成人手机在线| 国产熟女欧美一区二区| 丰满少妇做爰视频| 国产精品1区2区在线观看.| 日本免费a在线| 亚洲人成网站在线观看播放| 2018国产大陆天天弄谢| 亚洲av成人精品一区久久| 久久国产乱子免费精品| 免费看av在线观看网站| av播播在线观看一区| 亚洲国产精品成人久久小说| 国产免费一级a男人的天堂| 最后的刺客免费高清国语| 国产探花极品一区二区| 偷拍熟女少妇极品色| 久久久精品免费免费高清| 黑人高潮一二区| 久久精品国产鲁丝片午夜精品| 婷婷色综合大香蕉| 91精品国产九色| 只有这里有精品99| 免费在线观看成人毛片| 日韩精品青青久久久久久| 毛片女人毛片| 亚洲欧美一区二区三区黑人 | 波多野结衣巨乳人妻| 久久久午夜欧美精品| 久久久精品94久久精品| 国产亚洲一区二区精品| 免费看光身美女| 精品午夜福利在线看| 免费看av在线观看网站| h日本视频在线播放| 亚洲精品国产成人久久av| 99热这里只有是精品50| 日韩不卡一区二区三区视频在线| 欧美高清性xxxxhd video| 黄片无遮挡物在线观看| 韩国高清视频一区二区三区| 日产精品乱码卡一卡2卡三| 亚洲欧美一区二区三区国产| 中文字幕av成人在线电影| 亚洲高清免费不卡视频| 亚洲经典国产精华液单| 国模一区二区三区四区视频| av福利片在线观看| av在线亚洲专区| 国产高清国产精品国产三级 | 精品亚洲乱码少妇综合久久| 国产精品av视频在线免费观看| 毛片女人毛片| 欧美xxⅹ黑人| 日韩欧美精品免费久久| 夜夜看夜夜爽夜夜摸| 亚洲成人中文字幕在线播放| 午夜精品在线福利| 日韩国内少妇激情av| 免费电影在线观看免费观看| 午夜爱爱视频在线播放| 欧美最新免费一区二区三区| 亚洲精品色激情综合| a级毛色黄片| 最近中文字幕2019免费版| 乱人视频在线观看| 可以在线观看毛片的网站| 午夜激情久久久久久久| 午夜精品一区二区三区免费看| 亚洲一级一片aⅴ在线观看| 伊人久久国产一区二区| 亚洲一级一片aⅴ在线观看| 色尼玛亚洲综合影院| 国产精品av视频在线免费观看| 成年av动漫网址| eeuss影院久久| 欧美日韩亚洲高清精品| 亚洲欧洲国产日韩| 久久精品国产鲁丝片午夜精品| 日本wwww免费看| 亚洲图色成人| 国产伦精品一区二区三区视频9| 午夜精品国产一区二区电影 | 久久久成人免费电影| 亚洲av在线观看美女高潮| 日韩制服骚丝袜av| 丝袜喷水一区| 国产成人精品婷婷| 国产在线男女| 中国国产av一级| 色5月婷婷丁香| 99九九线精品视频在线观看视频| 麻豆成人av视频| 日本免费a在线| 99re6热这里在线精品视频| 一夜夜www| 黄色一级大片看看| 欧美成人一区二区免费高清观看| 欧美最新免费一区二区三区| 国产亚洲av嫩草精品影院| 国产免费福利视频在线观看| 国产老妇伦熟女老妇高清| 91久久精品电影网| 亚洲成人久久爱视频| 国产精品综合久久久久久久免费| 人妻系列 视频| 亚洲欧美精品专区久久| 午夜精品国产一区二区电影 | 亚洲18禁久久av| 国产精品麻豆人妻色哟哟久久 | 高清视频免费观看一区二区 | 亚洲第一区二区三区不卡| 久久久色成人| 亚洲精品日韩在线中文字幕| 色视频www国产| 人人妻人人澡欧美一区二区| 如何舔出高潮| 色综合亚洲欧美另类图片| 老师上课跳d突然被开到最大视频| 免费av毛片视频| 国产高清三级在线| 国产亚洲午夜精品一区二区久久 | 1000部很黄的大片| 99热6这里只有精品| 久久久亚洲精品成人影院| 水蜜桃什么品种好| 久久99蜜桃精品久久| 成人鲁丝片一二三区免费| 精华霜和精华液先用哪个| 亚洲av不卡在线观看| 又大又黄又爽视频免费| videos熟女内射| 国产精品国产三级国产av玫瑰| 成人综合一区亚洲| 男女国产视频网站| 天天一区二区日本电影三级| videossex国产| 永久网站在线| 精品人妻熟女av久视频| 人人妻人人澡欧美一区二区| 精品99又大又爽又粗少妇毛片| 日韩亚洲欧美综合| 免费看日本二区| 成人毛片a级毛片在线播放| 国产精品久久视频播放| av专区在线播放| 日韩大片免费观看网站| 狠狠精品人妻久久久久久综合| 在线观看av片永久免费下载| 禁无遮挡网站| 伦理电影大哥的女人| 乱人视频在线观看| 一级毛片aaaaaa免费看小| 一级毛片 在线播放| 99热这里只有是精品50| 综合色av麻豆| 人妻制服诱惑在线中文字幕| 午夜免费激情av| 22中文网久久字幕| 亚洲精品456在线播放app| 日本wwww免费看| 天堂网av新在线| 国产视频首页在线观看| 少妇高潮的动态图| 日本欧美国产在线视频| 亚洲精品日韩在线中文字幕| 亚洲av不卡在线观看| 国产 亚洲一区二区三区 | 高清毛片免费看| 最近最新中文字幕免费大全7| 日韩电影二区| 日本免费a在线| 老司机影院毛片| 麻豆成人午夜福利视频| 人人妻人人澡欧美一区二区| 色尼玛亚洲综合影院| 国产成人精品婷婷| 高清视频免费观看一区二区 | 国产女主播在线喷水免费视频网站 | 看十八女毛片水多多多| 国产精品嫩草影院av在线观看| 中文字幕av成人在线电影| 亚洲精品456在线播放app| 精品一区二区免费观看| 日韩欧美一区视频在线观看 | 天天一区二区日本电影三级| 国产在视频线精品| 久久久午夜欧美精品| 午夜精品国产一区二区电影 | 搡老乐熟女国产| 人妻制服诱惑在线中文字幕| 99久国产av精品国产电影| 亚洲精品久久久久久婷婷小说| 中文字幕制服av| 国产欧美另类精品又又久久亚洲欧美| 亚洲国产日韩欧美精品在线观看| 精品久久久久久久末码| 91久久精品电影网| 欧美 日韩 精品 国产| 我的女老师完整版在线观看| 精品久久久久久久人妻蜜臀av| 免费观看性生交大片5| 肉色欧美久久久久久久蜜桃 | 国产又色又爽无遮挡免| 69av精品久久久久久| 日本-黄色视频高清免费观看| 非洲黑人性xxxx精品又粗又长| 在线观看一区二区三区| 一级毛片aaaaaa免费看小| 你懂的网址亚洲精品在线观看| 亚洲精品第二区| 热99在线观看视频| 秋霞伦理黄片| 一二三四中文在线观看免费高清| 1000部很黄的大片| 欧美 日韩 精品 国产| av又黄又爽大尺度在线免费看| 国产午夜精品论理片| 午夜老司机福利剧场| 亚洲av中文av极速乱| 欧美极品一区二区三区四区| 少妇的逼好多水| 搞女人的毛片| 日本欧美国产在线视频| 国产 一区精品| 男女那种视频在线观看| 中国国产av一级| 国产老妇女一区| 麻豆乱淫一区二区| 精品久久久久久电影网| 国产又色又爽无遮挡免| 日日干狠狠操夜夜爽| 日本黄大片高清| 成人鲁丝片一二三区免费| 日韩大片免费观看网站| 人妻一区二区av| 婷婷色综合www| 亚洲欧美精品专区久久| 亚洲人成网站高清观看| 18禁在线播放成人免费| 欧美区成人在线视频| 国产精品国产三级国产专区5o| 99热网站在线观看| 欧美高清成人免费视频www| 一区二区三区四区激情视频| 日韩欧美 国产精品| 晚上一个人看的免费电影| 亚洲av成人av| 大香蕉久久网| 最近2019中文字幕mv第一页| 男人舔女人下体高潮全视频| a级一级毛片免费在线观看| 亚洲aⅴ乱码一区二区在线播放| 纵有疾风起免费观看全集完整版 | 2022亚洲国产成人精品| 久久久久久久久久久丰满| 国产高清国产精品国产三级 | 久久鲁丝午夜福利片| 麻豆成人午夜福利视频| 国内揄拍国产精品人妻在线| 国产午夜精品一二区理论片| 午夜老司机福利剧场| 久久99蜜桃精品久久| 99热这里只有是精品50| 91在线精品国自产拍蜜月| 日韩欧美 国产精品| 免费黄网站久久成人精品| 激情五月婷婷亚洲| 人体艺术视频欧美日本| 韩国高清视频一区二区三区| 亚洲18禁久久av| 免费大片黄手机在线观看| 精品人妻熟女av久视频| 国产精品久久久久久精品电影| 久久亚洲国产成人精品v| 免费观看的影片在线观看| 亚洲成人中文字幕在线播放| 可以在线观看毛片的网站| 嫩草影院入口| 国产乱人视频| 久久精品国产鲁丝片午夜精品| 亚洲av电影在线观看一区二区三区 | 成人毛片60女人毛片免费| 午夜精品国产一区二区电影 | 日韩av在线大香蕉| 91午夜精品亚洲一区二区三区| 国产在线一区二区三区精| 99九九线精品视频在线观看视频| 午夜福利在线观看免费完整高清在| 嘟嘟电影网在线观看| 国产亚洲午夜精品一区二区久久 | 亚洲av.av天堂| 国产午夜精品论理片| 亚洲精品日本国产第一区| 亚洲美女视频黄频| 中文资源天堂在线| 久久久色成人| 国产精品av视频在线免费观看| 中文资源天堂在线| 好男人在线观看高清免费视频| 日韩不卡一区二区三区视频在线| 欧美日本视频| 国产精品嫩草影院av在线观看| 精品国产三级普通话版| 日日摸夜夜添夜夜添av毛片| 国产老妇伦熟女老妇高清| 日本wwww免费看|