• <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
    欧美激情久久久久久爽电影| 国内精品一区二区在线观看| 窝窝影院91人妻| 中文字幕免费在线视频6| 亚洲av免费高清在线观看| 欧美成人性av电影在线观看| 成人午夜高清在线视频| 性插视频无遮挡在线免费观看| 欧美一区二区亚洲| 欧美成人性av电影在线观看| 日本 欧美在线| 日本免费a在线| 亚洲内射少妇av| 免费高清视频大片| 女生性感内裤真人,穿戴方法视频| 深夜精品福利| a级毛片a级免费在线| 午夜福利免费观看在线| 亚洲五月天丁香| 特大巨黑吊av在线直播| 国产精品国产高清国产av| 欧美日本亚洲视频在线播放| 国产一区二区三区在线臀色熟女| 99热这里只有是精品在线观看 | 国产一区二区激情短视频| 国产亚洲精品久久久com| 日韩欧美在线二视频| 国内精品一区二区在线观看| 久久婷婷人人爽人人干人人爱| 日日干狠狠操夜夜爽| xxxwww97欧美| 国产探花极品一区二区| 97人妻精品一区二区三区麻豆| 国产精品爽爽va在线观看网站| 99在线视频只有这里精品首页| 一本精品99久久精品77| 1000部很黄的大片| 久久精品影院6| 青草久久国产| 麻豆一二三区av精品| 白带黄色成豆腐渣| 国产精品久久久久久亚洲av鲁大| 亚洲人成网站在线播放欧美日韩| 午夜激情欧美在线| 亚洲欧美日韩高清专用| 毛片女人毛片| 亚洲欧美清纯卡通| 男插女下体视频免费在线播放| 毛片一级片免费看久久久久 | 亚洲av免费在线观看| 一个人观看的视频www高清免费观看| 久久性视频一级片| 能在线免费观看的黄片| 欧美色欧美亚洲另类二区| 久久久久久九九精品二区国产| 亚洲真实伦在线观看| 国产三级在线视频| 最新中文字幕久久久久| 狠狠狠狠99中文字幕| 一个人观看的视频www高清免费观看| www.www免费av| 一本一本综合久久| 男人狂女人下面高潮的视频| 又黄又爽又刺激的免费视频.| 免费观看精品视频网站| 少妇熟女aⅴ在线视频| 特级一级黄色大片| 精品一区二区三区视频在线观看免费| 男人舔奶头视频| 国产午夜精品久久久久久一区二区三区 | 久久亚洲精品不卡| 亚洲国产欧美人成| 国产精品久久视频播放| 亚洲人成伊人成综合网2020| 中文在线观看免费www的网站| 高潮久久久久久久久久久不卡| 国产在线男女| 精品福利观看| 在线观看av片永久免费下载| 婷婷六月久久综合丁香| 欧美性感艳星| 午夜精品在线福利| 成熟少妇高潮喷水视频| 成人欧美大片| av福利片在线观看| 热99在线观看视频| 亚洲美女视频黄频| 亚洲av成人精品一区久久| 脱女人内裤的视频| 亚洲va日本ⅴa欧美va伊人久久| 麻豆久久精品国产亚洲av| 国产成年人精品一区二区| 禁无遮挡网站| 韩国av一区二区三区四区| 亚洲精品一区av在线观看| 制服丝袜大香蕉在线| 性欧美人与动物交配| 日韩国内少妇激情av| 午夜福利欧美成人| 毛片女人毛片| 国产精品人妻久久久久久| 日韩欧美国产一区二区入口| 欧美激情在线99| 欧美极品一区二区三区四区| 欧美日韩综合久久久久久 | 搡女人真爽免费视频火全软件 | 91麻豆精品激情在线观看国产| 国产精品久久视频播放| 丁香六月欧美| 自拍偷自拍亚洲精品老妇| 偷拍熟女少妇极品色| 性欧美人与动物交配| 黄色配什么色好看| 国产成年人精品一区二区| 简卡轻食公司| 黄色一级大片看看| 久久人人爽人人爽人人片va | 日本黄大片高清| 亚洲内射少妇av| 很黄的视频免费| 国产亚洲精品av在线| 午夜福利18| а√天堂www在线а√下载| 内地一区二区视频在线| 国产老妇女一区| 亚州av有码| 可以在线观看的亚洲视频| 看黄色毛片网站| 亚洲真实伦在线观看| 亚洲av第一区精品v没综合| 久久精品国产亚洲av天美| 熟女电影av网| 嫩草影院精品99| 国产成人欧美在线观看| 一区二区三区激情视频| 亚洲内射少妇av| 国产视频内射| 99久久九九国产精品国产免费| 久久久成人免费电影| 成年人黄色毛片网站| 午夜视频国产福利| 精华霜和精华液先用哪个| 欧美黑人巨大hd| 国产综合懂色| 村上凉子中文字幕在线| 最近最新免费中文字幕在线| 观看免费一级毛片| 欧美一区二区精品小视频在线| 亚洲av电影在线进入| 两人在一起打扑克的视频| 免费看光身美女| 国产免费av片在线观看野外av| 男人和女人高潮做爰伦理| 欧美潮喷喷水| 少妇丰满av| 一本久久中文字幕| 可以在线观看毛片的网站| 亚洲精品一区av在线观看| 亚洲国产精品成人综合色| 国产在线精品亚洲第一网站| 欧美日本亚洲视频在线播放| 毛片女人毛片| 99久久99久久久精品蜜桃| 少妇丰满av| 高清毛片免费观看视频网站| 亚洲成av人片在线播放无| 少妇被粗大猛烈的视频| 熟女人妻精品中文字幕| 欧美最新免费一区二区三区 | 别揉我奶头 嗯啊视频| 日本黄色视频三级网站网址| 欧美性猛交╳xxx乱大交人| 亚洲av不卡在线观看| 亚洲电影在线观看av| 丰满乱子伦码专区| 又爽又黄a免费视频| 在线观看av片永久免费下载| 欧美xxxx黑人xx丫x性爽| 欧美潮喷喷水| 久久久成人免费电影| 熟女人妻精品中文字幕| 午夜日韩欧美国产| 热99在线观看视频| 99精品久久久久人妻精品| 国产伦一二天堂av在线观看| 亚洲精品色激情综合| 一进一出抽搐动态| 亚洲专区国产一区二区| 免费观看人在逋| 又爽又黄a免费视频| 久久精品人妻少妇| 18美女黄网站色大片免费观看| 精品国内亚洲2022精品成人| 久久精品久久久久久噜噜老黄 | 99久久精品一区二区三区| 精华霜和精华液先用哪个| 国产亚洲精品久久久com| 韩国av一区二区三区四区| 久久久成人免费电影| 51国产日韩欧美| 成人高潮视频无遮挡免费网站| www.色视频.com| 热99re8久久精品国产| 免费无遮挡裸体视频| АⅤ资源中文在线天堂| 亚洲精品456在线播放app | 国产精品,欧美在线| 精品无人区乱码1区二区| 亚洲经典国产精华液单 | 又黄又爽又免费观看的视频| 综合色av麻豆| 国产 一区 欧美 日韩| 男人舔奶头视频| 国产精品久久视频播放| 国产探花极品一区二区| 久久久久久久午夜电影| 欧美高清性xxxxhd video| 久久热精品热| 免费搜索国产男女视频| 9191精品国产免费久久| www.色视频.com| 久久久久性生活片| 日韩欧美 国产精品| 欧美在线一区亚洲| 9191精品国产免费久久| 精品乱码久久久久久99久播| 中文字幕免费在线视频6| 国产精品综合久久久久久久免费| а√天堂www在线а√下载| 亚洲三级黄色毛片| 丰满的人妻完整版| 特大巨黑吊av在线直播| 日本五十路高清| 97热精品久久久久久| 中文资源天堂在线| 少妇的逼好多水| 国产精品久久电影中文字幕| 无遮挡黄片免费观看| 久久久色成人| 国产亚洲精品久久久久久毛片| 免费大片18禁| 人人妻人人看人人澡| 国产91精品成人一区二区三区| 国产日本99.免费观看| 在线观看免费视频日本深夜| 欧美zozozo另类| 久久久久久久久中文| 91麻豆精品激情在线观看国产| 最新在线观看一区二区三区| 亚洲专区中文字幕在线| 欧美xxxx黑人xx丫x性爽| 国产熟女xx| 国产精品,欧美在线| 亚洲精品456在线播放app | 国产精品一区二区性色av| 日本精品一区二区三区蜜桃| 精品人妻熟女av久视频| 好看av亚洲va欧美ⅴa在| av中文乱码字幕在线| 18禁在线播放成人免费| 国产v大片淫在线免费观看| 国产精品乱码一区二三区的特点| 美女大奶头视频| 成人午夜高清在线视频| 九色国产91popny在线| 欧美三级亚洲精品| 中文字幕熟女人妻在线| 精品99又大又爽又粗少妇毛片 | av在线蜜桃| 国产一区二区激情短视频| 精品无人区乱码1区二区| 91九色精品人成在线观看| 亚洲成人精品中文字幕电影| 最近最新免费中文字幕在线| av在线老鸭窝| 自拍偷自拍亚洲精品老妇| 国产伦人伦偷精品视频| 午夜久久久久精精品| 国内精品美女久久久久久| 成年女人永久免费观看视频| 中文字幕免费在线视频6| 国产av在哪里看| 久久久色成人| 国产极品精品免费视频能看的| 亚洲aⅴ乱码一区二区在线播放| 亚洲成av人片在线播放无| 天堂√8在线中文| 日本在线视频免费播放| 3wmmmm亚洲av在线观看| 久久精品国产99精品国产亚洲性色| 每晚都被弄得嗷嗷叫到高潮| 国产亚洲av嫩草精品影院| 成人国产综合亚洲| 国内毛片毛片毛片毛片毛片| 好看av亚洲va欧美ⅴa在| 精品久久国产蜜桃| 可以在线观看毛片的网站| 久久中文看片网| 国产精品久久电影中文字幕| 午夜a级毛片| 欧美中文日本在线观看视频| 精品久久久久久,| 中文字幕精品亚洲无线码一区| 身体一侧抽搐| 日韩亚洲欧美综合| 精品国产亚洲在线| 午夜福利视频1000在线观看| 日本黄色视频三级网站网址| 免费观看的影片在线观看| 一本综合久久免费| 中文字幕熟女人妻在线| netflix在线观看网站| 麻豆国产97在线/欧美| 亚洲人成网站高清观看| 国产成年人精品一区二区| 精品人妻一区二区三区麻豆 | 又紧又爽又黄一区二区| 欧美成人一区二区免费高清观看| 亚洲经典国产精华液单 | 国产精品一区二区性色av| a在线观看视频网站| 成人国产综合亚洲| 夜夜爽天天搞| av天堂中文字幕网| 日韩av在线大香蕉| 熟妇人妻久久中文字幕3abv| 国产精品女同一区二区软件 | 午夜免费成人在线视频| 国产成人av教育| 成人av一区二区三区在线看| 日韩有码中文字幕| 中文字幕av成人在线电影| 久9热在线精品视频| 九九在线视频观看精品| 男人的好看免费观看在线视频| 少妇高潮的动态图| 男女视频在线观看网站免费| 国产伦在线观看视频一区| 国产黄a三级三级三级人| 亚洲在线观看片| 丰满人妻一区二区三区视频av| 在线播放国产精品三级| 亚洲激情在线av| 在线播放国产精品三级| 国产精品,欧美在线| 精品久久久久久久久久久久久| 国产主播在线观看一区二区| 黄色丝袜av网址大全| 少妇丰满av| 人妻久久中文字幕网| 亚洲自拍偷在线| 麻豆av噜噜一区二区三区| av中文乱码字幕在线| 亚洲欧美激情综合另类| 亚洲成人精品中文字幕电影| 亚洲欧美精品综合久久99| 久久久久久久亚洲中文字幕 | 又爽又黄无遮挡网站| 欧美区成人在线视频| 99久久无色码亚洲精品果冻| 人妻制服诱惑在线中文字幕| 麻豆国产97在线/欧美| 黄色丝袜av网址大全| 国产免费av片在线观看野外av| 亚洲不卡免费看| 婷婷六月久久综合丁香| 国产综合懂色| www.999成人在线观看| 久久久精品大字幕| 亚洲熟妇中文字幕五十中出| 国产精品野战在线观看| 久久久精品欧美日韩精品| 9191精品国产免费久久| 午夜亚洲福利在线播放| 黄色配什么色好看| 国产精品久久视频播放| 精品国产亚洲在线| 国产人妻一区二区三区在| 国产高清视频在线播放一区| 可以在线观看的亚洲视频| 亚洲专区中文字幕在线| 97超级碰碰碰精品色视频在线观看| 成人国产一区最新在线观看| 中文字幕精品亚洲无线码一区| 国产精品久久久久久亚洲av鲁大| 性欧美人与动物交配| 亚洲精品日韩av片在线观看| 欧美另类亚洲清纯唯美| 香蕉av资源在线| 九色国产91popny在线| 亚洲自偷自拍三级| 天堂av国产一区二区熟女人妻| 亚洲av日韩精品久久久久久密| 国产在线男女| 欧美一区二区亚洲| 99久久九九国产精品国产免费| 啪啪无遮挡十八禁网站| 欧美不卡视频在线免费观看| 欧美xxxx黑人xx丫x性爽| 国产亚洲精品久久久久久毛片| 桃红色精品国产亚洲av| 一进一出抽搐gif免费好疼| av女优亚洲男人天堂| 中文字幕av成人在线电影| 精品无人区乱码1区二区| 亚洲精品一区av在线观看| 女人被狂操c到高潮| 亚洲成人中文字幕在线播放| 直男gayav资源| 亚洲综合色惰| 国内精品久久久久精免费| 少妇的逼好多水| 小蜜桃在线观看免费完整版高清| 少妇裸体淫交视频免费看高清| 亚洲欧美日韩东京热| 毛片一级片免费看久久久久 | 国产精品影院久久| 99久国产av精品| 亚洲片人在线观看| 中文在线观看免费www的网站| 亚洲自偷自拍三级| 看黄色毛片网站| 免费大片18禁| 波野结衣二区三区在线| 欧美乱妇无乱码| 国产精品综合久久久久久久免费| 国产真实乱freesex| 亚洲专区国产一区二区| 观看美女的网站| 成人欧美大片| 国产精品av视频在线免费观看| 一进一出抽搐动态| 久久国产精品影院| 少妇裸体淫交视频免费看高清| 午夜日韩欧美国产| 国产精品一区二区三区四区久久| 婷婷丁香在线五月| 不卡一级毛片| 天天躁日日操中文字幕| 9191精品国产免费久久| 久久人人爽人人爽人人片va | 在线a可以看的网站| 51午夜福利影视在线观看| 小蜜桃在线观看免费完整版高清| www.999成人在线观看| 夜夜躁狠狠躁天天躁| 在线免费观看的www视频| 精品久久久久久成人av| 免费看a级黄色片| 国产美女午夜福利| 国产一区二区三区视频了| 一边摸一边抽搐一进一小说| 天天一区二区日本电影三级| 亚洲狠狠婷婷综合久久图片| 欧美潮喷喷水| 亚洲熟妇中文字幕五十中出| 午夜久久久久精精品| 国产主播在线观看一区二区| 国产v大片淫在线免费观看| 日韩欧美精品v在线| 美女 人体艺术 gogo| 国产精品一区二区性色av| 亚洲精品影视一区二区三区av| 成人永久免费在线观看视频| 亚洲精品粉嫩美女一区| 十八禁国产超污无遮挡网站| 91麻豆精品激情在线观看国产| 99久国产av精品| 2021天堂中文幕一二区在线观| 欧美成人一区二区免费高清观看| 中文在线观看免费www的网站| 精品乱码久久久久久99久播| 99热这里只有精品一区| 少妇人妻一区二区三区视频| 亚洲欧美激情综合另类| 精品人妻1区二区| 最新在线观看一区二区三区| 亚洲国产精品999在线| 首页视频小说图片口味搜索| 自拍偷自拍亚洲精品老妇| 熟女人妻精品中文字幕| 国产午夜福利久久久久久| 麻豆国产av国片精品| 国产高清视频在线播放一区| 一边摸一边抽搐一进一小说| 一级黄色大片毛片| 成人无遮挡网站| 亚洲专区国产一区二区| www.www免费av| 在线天堂最新版资源| 一本一本综合久久| 国产爱豆传媒在线观看| 18禁黄网站禁片午夜丰满| 成人国产一区最新在线观看| 久99久视频精品免费| 成人国产一区最新在线观看| 欧美xxxx性猛交bbbb| 久久热精品热| 国产美女午夜福利| 大型黄色视频在线免费观看| 男人的好看免费观看在线视频| 亚洲电影在线观看av| 久久久色成人| 又紧又爽又黄一区二区| 久久久久久大精品| 十八禁国产超污无遮挡网站| 亚洲av美国av| 脱女人内裤的视频| 成年女人永久免费观看视频| 人人妻人人看人人澡| 午夜福利免费观看在线| 搡老熟女国产l中国老女人| 欧美日韩乱码在线| 欧美精品啪啪一区二区三区| 亚洲国产精品成人综合色| 波野结衣二区三区在线| 欧美又色又爽又黄视频| 久99久视频精品免费| 97碰自拍视频| 男人狂女人下面高潮的视频| 一a级毛片在线观看| 天堂影院成人在线观看| 97超视频在线观看视频| 国产成人a区在线观看| 深夜精品福利| 久久久精品大字幕| 男人舔女人下体高潮全视频| 亚洲第一电影网av| 色吧在线观看| 色尼玛亚洲综合影院| 日韩欧美一区二区三区在线观看| 国产又黄又爽又无遮挡在线| 一二三四社区在线视频社区8| 成人美女网站在线观看视频| 亚洲av美国av| 亚洲三级黄色毛片| 少妇的逼好多水| 美女被艹到高潮喷水动态| 国产精品一区二区三区四区久久| av在线老鸭窝| 亚洲国产欧洲综合997久久,| 啦啦啦韩国在线观看视频| 国产精品三级大全| 少妇的逼水好多| 欧美性感艳星| 少妇熟女aⅴ在线视频| 啪啪无遮挡十八禁网站| 国产老妇女一区| 欧美xxxx黑人xx丫x性爽| 久久国产精品影院| 高清毛片免费观看视频网站| 国产视频一区二区在线看| 中国美女看黄片| 欧美三级亚洲精品| 亚洲成人久久爱视频| 嫩草影视91久久| 欧美中文日本在线观看视频| av国产免费在线观看| 国产精品永久免费网站| 香蕉av资源在线| 国产大屁股一区二区在线视频| 免费电影在线观看免费观看| 十八禁国产超污无遮挡网站| 日本精品一区二区三区蜜桃| 丁香六月欧美| 欧美一区二区亚洲| 又黄又爽又免费观看的视频| 中国美女看黄片| 亚洲国产高清在线一区二区三| 久久久国产成人免费| 此物有八面人人有两片| 成年免费大片在线观看| 中文字幕人成人乱码亚洲影| 国产精品永久免费网站| www.999成人在线观看| 久久久久久久久中文| 性色av乱码一区二区三区2| 欧美激情久久久久久爽电影| 性色avwww在线观看| 一进一出抽搐gif免费好疼| 日本撒尿小便嘘嘘汇集6| 亚洲精品成人久久久久久| 久久久久精品国产欧美久久久| 长腿黑丝高跟| 搡老熟女国产l中国老女人| 欧美在线一区亚洲| 亚洲中文字幕日韩| 51国产日韩欧美| 国产真实伦视频高清在线观看 | 久久久色成人| 欧美又色又爽又黄视频| 精品99又大又爽又粗少妇毛片 | 99热这里只有是精品在线观看 | 搞女人的毛片| 熟女人妻精品中文字幕| 国产精品一及| 人妻制服诱惑在线中文字幕| 国产成人av教育| 99热只有精品国产| 国产精品一区二区三区四区免费观看 | 亚洲精品色激情综合| 黄色配什么色好看| 很黄的视频免费| 亚洲国产欧洲综合997久久,| h日本视频在线播放| 国产69精品久久久久777片| 变态另类丝袜制服| 在线观看免费视频日本深夜| 如何舔出高潮| 国产一区二区三区在线臀色熟女| 免费无遮挡裸体视频| 久久久色成人| 亚洲精品一区av在线观看| 国产久久久一区二区三区| 国产成人福利小说|