• <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新网站| 久热久热在线精品观看| 亚洲精品自拍成人| 丝袜喷水一区| 别揉我奶头 嗯啊视频| 丰满乱子伦码专区| 亚洲国产高清在线一区二区三| 亚洲,一卡二卡三卡| 91aial.com中文字幕在线观看| 久久久欧美国产精品| 国语对白做爰xxxⅹ性视频网站| 亚洲精品日韩在线中文字幕| 亚洲国产精品999| 亚洲人成网站高清观看| 在线观看三级黄色| 久久久久久久亚洲中文字幕| 亚洲精品久久午夜乱码| 国产精品无大码| 特级一级黄色大片| 亚洲四区av| 日韩制服骚丝袜av| 成年女人看的毛片在线观看| 免费少妇av软件| 欧美成人精品欧美一级黄| 伦理电影大哥的女人| 欧美一级a爱片免费观看看| 蜜桃亚洲精品一区二区三区| 色婷婷久久久亚洲欧美| 中文精品一卡2卡3卡4更新| 国产淫语在线视频| 成人无遮挡网站| 街头女战士在线观看网站| 精品久久久久久久人妻蜜臀av| 内地一区二区视频在线| 嫩草影院入口| 中文精品一卡2卡3卡4更新| av国产免费在线观看| 日韩av不卡免费在线播放| 夫妻性生交免费视频一级片| kizo精华| 日本一本二区三区精品| 国产爽快片一区二区三区| 国产淫片久久久久久久久| 欧美日韩国产mv在线观看视频 | 熟妇人妻不卡中文字幕| 国产美女午夜福利| av国产久精品久网站免费入址| 久久精品久久精品一区二区三区| 一级二级三级毛片免费看| 久久影院123| 亚洲国产av新网站| 国产黄色免费在线视频| 亚洲av二区三区四区| 亚洲av中文字字幕乱码综合| 免费黄频网站在线观看国产| 国产精品一二三区在线看| 男插女下体视频免费在线播放| 高清欧美精品videossex| 一级毛片我不卡| 国产亚洲av嫩草精品影院| 大片免费播放器 马上看| av在线播放精品| 国国产精品蜜臀av免费| 亚洲精品国产色婷婷电影| 嫩草影院精品99| 亚洲电影在线观看av| 国产综合懂色| 在线免费观看不下载黄p国产| 色哟哟·www| 久久久亚洲精品成人影院| 国产精品麻豆人妻色哟哟久久| 建设人人有责人人尽责人人享有的 | 日韩成人av中文字幕在线观看| 久久ye,这里只有精品| 国产精品三级大全| 日本猛色少妇xxxxx猛交久久| 2018国产大陆天天弄谢| 国产成人精品婷婷| 日韩三级伦理在线观看| 22中文网久久字幕| 久久这里有精品视频免费| 久久99热6这里只有精品| 日本av手机在线免费观看| 亚洲伊人久久精品综合| 网址你懂的国产日韩在线| 国产精品.久久久| 国产一区二区三区综合在线观看 | 亚洲内射少妇av| 亚洲自偷自拍三级| 大片电影免费在线观看免费| 亚洲av福利一区| 99久久中文字幕三级久久日本| 亚洲欧美成人综合另类久久久| 久久国产乱子免费精品| 天天躁夜夜躁狠狠久久av| 纵有疾风起免费观看全集完整版| 熟女电影av网| 肉色欧美久久久久久久蜜桃 | 午夜免费鲁丝| 国模一区二区三区四区视频| 国产伦理片在线播放av一区| 别揉我奶头 嗯啊视频| 真实男女啪啪啪动态图| 免费少妇av软件| 秋霞伦理黄片| 熟女av电影| 99久久人妻综合| 亚洲国产精品国产精品| av黄色大香蕉| 亚洲在线观看片| 免费av观看视频| 亚洲国产日韩一区二区| 欧美日韩一区二区视频在线观看视频在线 | av卡一久久| 亚洲自拍偷在线| 国产成人精品福利久久| 日韩欧美精品免费久久| 又粗又硬又长又爽又黄的视频| 成人国产麻豆网| 永久免费av网站大全| 亚洲人成网站在线观看播放| 久久久欧美国产精品| 国产精品人妻久久久久久| 免费看不卡的av| 99re6热这里在线精品视频| 国产精品一区二区三区四区免费观看| 波野结衣二区三区在线| 插逼视频在线观看| 边亲边吃奶的免费视频| 日韩人妻高清精品专区| 亚洲av二区三区四区| 国产毛片在线视频| 中文字幕久久专区| 99热这里只有是精品50| 中文在线观看免费www的网站| 99久久人妻综合| 国产亚洲av片在线观看秒播厂| 乱码一卡2卡4卡精品| 成人美女网站在线观看视频| 日韩电影二区| 人妻 亚洲 视频| 18+在线观看网站| 亚洲无线观看免费| 天天一区二区日本电影三级| 日韩在线高清观看一区二区三区| 搞女人的毛片| 国产成人精品福利久久| 亚洲经典国产精华液单| 老司机影院毛片| 欧美日韩一区二区视频在线观看视频在线 | 80岁老熟妇乱子伦牲交| 久久久久国产网址| 熟女电影av网| 一区二区三区精品91| 精品国产乱码久久久久久小说| 中文欧美无线码| 青春草视频在线免费观看| 精品久久久噜噜| 免费看光身美女| 肉色欧美久久久久久久蜜桃 | 大香蕉97超碰在线| 国产黄色视频一区二区在线观看| 国产一级毛片在线| 尤物成人国产欧美一区二区三区| 国产精品麻豆人妻色哟哟久久| 亚洲精华国产精华液的使用体验| 亚洲图色成人| 亚洲久久久久久中文字幕| 毛片女人毛片| 69人妻影院| 欧美激情在线99| 久久久久久久午夜电影| 亚洲精品国产av成人精品| 亚洲美女搞黄在线观看| 人妻一区二区av| 免费黄频网站在线观看国产| 特大巨黑吊av在线直播| 亚洲在线观看片| 亚洲丝袜综合中文字幕| 日韩欧美一区视频在线观看 | 舔av片在线| 亚洲国产最新在线播放| 国产美女午夜福利| 国产91av在线免费观看| 亚洲三级黄色毛片| 尾随美女入室| 国产精品一区二区三区四区免费观看| 日日摸夜夜添夜夜添av毛片| 成人欧美大片| 亚洲国产精品999| 久久久久久久久久人人人人人人| av天堂中文字幕网| 麻豆久久精品国产亚洲av| 久久久久久九九精品二区国产| av在线观看视频网站免费| 国产精品秋霞免费鲁丝片| xxx大片免费视频| 搡老乐熟女国产| 69av精品久久久久久| 热99国产精品久久久久久7| 欧美亚洲 丝袜 人妻 在线| 日韩一本色道免费dvd| 最新中文字幕久久久久| 国产国拍精品亚洲av在线观看| 综合色丁香网| 最近手机中文字幕大全| 久久精品久久久久久久性| 日日啪夜夜爽| 不卡视频在线观看欧美| 久久精品人妻少妇| 国产欧美另类精品又又久久亚洲欧美| 亚洲欧美日韩卡通动漫| 狂野欧美白嫩少妇大欣赏| 在线观看人妻少妇| 久热这里只有精品99| 丝袜喷水一区| 91精品国产九色| 别揉我奶头 嗯啊视频| 国产亚洲av片在线观看秒播厂| 伦精品一区二区三区| 边亲边吃奶的免费视频| 又大又黄又爽视频免费| 亚洲av免费在线观看| 久久久精品欧美日韩精品| 日本欧美国产在线视频| 女人被狂操c到高潮| 秋霞伦理黄片| 欧美区成人在线视频| 亚洲成人精品中文字幕电影| 亚洲欧洲日产国产| 久久午夜福利片| 久久久久久久大尺度免费视频| 赤兔流量卡办理| 我的老师免费观看完整版| 永久免费av网站大全| 国产免费一级a男人的天堂| 老司机影院毛片| 久久久久久久久久人人人人人人| 国产精品.久久久| 色综合色国产| 国产女主播在线喷水免费视频网站| 制服丝袜香蕉在线| 插逼视频在线观看| 美女主播在线视频| 有码 亚洲区| 亚洲国产精品成人综合色| 亚洲精品第二区| 国产日韩欧美在线精品| 国产成人精品福利久久| 汤姆久久久久久久影院中文字幕| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 久久97久久精品| 午夜老司机福利剧场| 少妇 在线观看| 婷婷色av中文字幕| 极品少妇高潮喷水抽搐| 美女被艹到高潮喷水动态| 日韩国内少妇激情av| 精品国产一区二区三区久久久樱花 | 91精品伊人久久大香线蕉| 国产 精品1| 亚洲精品一区蜜桃| 国产精品99久久久久久久久| 欧美xxⅹ黑人| 精品国产三级普通话版| 久久ye,这里只有精品| 高清欧美精品videossex| 女的被弄到高潮叫床怎么办| 国产又色又爽无遮挡免| 亚洲欧美一区二区三区黑人 | 亚洲精品,欧美精品| www.色视频.com| 青青草视频在线视频观看| 狂野欧美激情性bbbbbb| 免费少妇av软件| 看十八女毛片水多多多| 日韩制服骚丝袜av| 69av精品久久久久久| 亚洲人成网站在线播| 成人欧美大片| 免费观看av网站的网址| 国产欧美另类精品又又久久亚洲欧美| 久久6这里有精品| 韩国av在线不卡| 中文字幕久久专区| 国产老妇女一区| 看非洲黑人一级黄片| av在线播放精品| 女人被狂操c到高潮| 欧美日本视频| 身体一侧抽搐| 爱豆传媒免费全集在线观看| videossex国产| 狠狠精品人妻久久久久久综合| 欧美日本视频| 亚洲精品日韩av片在线观看| 日韩强制内射视频| 秋霞在线观看毛片| 熟女av电影| 日日啪夜夜撸| 夫妻午夜视频| 男人和女人高潮做爰伦理| 国产欧美日韩一区二区三区在线 | 亚洲精品国产av成人精品| 日韩亚洲欧美综合| 蜜桃久久精品国产亚洲av| 观看免费一级毛片| 在线a可以看的网站| 欧美区成人在线视频| 日韩人妻高清精品专区| 国产一区二区亚洲精品在线观看| 精品久久久精品久久久| 成年免费大片在线观看| 可以在线观看毛片的网站| 国产免费福利视频在线观看| av免费在线看不卡| 成人无遮挡网站| 日本一本二区三区精品| 欧美日韩一区二区视频在线观看视频在线 | 在线观看一区二区三区| 97在线视频观看| 男人爽女人下面视频在线观看| 欧美成人精品欧美一级黄| 国产成人一区二区在线| 2022亚洲国产成人精品| 欧美国产精品一级二级三级 | 国产免费福利视频在线观看| 成人国产av品久久久| 免费大片18禁| 狂野欧美白嫩少妇大欣赏| 午夜免费男女啪啪视频观看| 只有这里有精品99| 视频中文字幕在线观看| 亚洲久久久久久中文字幕| 老司机影院毛片| 亚洲综合色惰| 禁无遮挡网站| 久久久a久久爽久久v久久| 国产成人freesex在线| 亚洲精品视频女| 一边亲一边摸免费视频| 久久久久网色| 久久久久久久久久久丰满| 免费大片黄手机在线观看| freevideosex欧美| 一级二级三级毛片免费看| 欧美成人精品欧美一级黄| 免费av不卡在线播放| 黑人高潮一二区| 五月天丁香电影| 赤兔流量卡办理| 国产一区二区三区av在线| 女人十人毛片免费观看3o分钟| 成人亚洲精品av一区二区| 夫妻性生交免费视频一级片| 国产精品一及| 一本一本综合久久| 小蜜桃在线观看免费完整版高清| av福利片在线观看| av在线观看视频网站免费| 精品午夜福利在线看| 欧美zozozo另类| 最近2019中文字幕mv第一页| 精品人妻一区二区三区麻豆| 欧美老熟妇乱子伦牲交| 亚洲激情五月婷婷啪啪| 少妇丰满av| 午夜福利视频精品| 天堂俺去俺来也www色官网| 国产探花极品一区二区| 欧美bdsm另类| 亚洲国产精品成人综合色| 久久精品国产a三级三级三级| 在线天堂最新版资源| 亚洲精品色激情综合| 国产亚洲av嫩草精品影院| 在线观看一区二区三区激情| 午夜福利高清视频| 综合色av麻豆| 男人狂女人下面高潮的视频| 一个人看的www免费观看视频| 七月丁香在线播放| 春色校园在线视频观看| 国产爽快片一区二区三区| 国产av码专区亚洲av| 一边亲一边摸免费视频| 丝袜脚勾引网站| 大又大粗又爽又黄少妇毛片口| 国产高清不卡午夜福利| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 69人妻影院| 国产欧美日韩一区二区三区在线 | 一级毛片黄色毛片免费观看视频| 亚洲国产欧美在线一区| 国产精品一区www在线观看| 日韩在线高清观看一区二区三区| 欧美日韩在线观看h| 大片免费播放器 马上看| 国产一区有黄有色的免费视频| 国产女主播在线喷水免费视频网站| 成人特级av手机在线观看| 久久人人爽人人片av| h日本视频在线播放| 国产一区亚洲一区在线观看| 丝袜脚勾引网站| 色5月婷婷丁香| 国产精品av视频在线免费观看| 亚洲精品日本国产第一区| 一个人看视频在线观看www免费| 草草在线视频免费看| 五月天丁香电影| 国产成人午夜福利电影在线观看| 国产黄色免费在线视频| 精品久久久久久久末码| 欧美一级a爱片免费观看看| 神马国产精品三级电影在线观看| 舔av片在线| 国产精品久久久久久精品古装| 亚洲天堂国产精品一区在线| 国产毛片a区久久久久| 黑人高潮一二区| 久久精品久久久久久久性| 色网站视频免费| 国产欧美另类精品又又久久亚洲欧美| 亚洲欧美日韩另类电影网站 | 99热6这里只有精品| 欧美激情久久久久久爽电影| 久久人人爽人人爽人人片va| 99热这里只有是精品在线观看| 国产高清不卡午夜福利| 国产黄a三级三级三级人| 特大巨黑吊av在线直播| 欧美xxⅹ黑人| 日日撸夜夜添| 日韩亚洲欧美综合| 成人午夜精彩视频在线观看| 九九在线视频观看精品| av免费在线看不卡| 少妇熟女欧美另类| 中国美白少妇内射xxxbb| 中文欧美无线码| 国产黄片视频在线免费观看| 亚洲av不卡在线观看| 精品酒店卫生间| 搡女人真爽免费视频火全软件| av在线播放精品| 99热这里只有精品一区| 熟女电影av网| 777米奇影视久久| 美女视频免费永久观看网站| 欧美日韩一区二区视频在线观看视频在线 | 久久精品国产自在天天线| 国产老妇伦熟女老妇高清| 国产精品人妻久久久久久| 九九久久精品国产亚洲av麻豆| 99热这里只有是精品在线观看| 免费av不卡在线播放| 丰满乱子伦码专区| 老师上课跳d突然被开到最大视频| 69人妻影院| 秋霞在线观看毛片| .国产精品久久| 搡老乐熟女国产| 国产日韩欧美亚洲二区| 亚洲精品,欧美精品| 午夜福利视频1000在线观看| 国产伦在线观看视频一区| 免费不卡的大黄色大毛片视频在线观看| 久久精品夜色国产| 三级男女做爰猛烈吃奶摸视频| 欧美性感艳星| 黄色欧美视频在线观看| 搡女人真爽免费视频火全软件| 一级黄片播放器| 欧美xxⅹ黑人| 亚洲国产精品成人综合色| 精品久久久久久久人妻蜜臀av| 色吧在线观看| 国产一级毛片在线| 日韩一本色道免费dvd| 美女cb高潮喷水在线观看| 成人高潮视频无遮挡免费网站| 欧美少妇被猛烈插入视频| 老女人水多毛片| 男女下面进入的视频免费午夜| 久久人人爽人人片av| 国产成人精品久久久久久| 国产免费一区二区三区四区乱码| 内射极品少妇av片p| 国产男女内射视频| 丝瓜视频免费看黄片| 2022亚洲国产成人精品| 日本欧美国产在线视频| 成年版毛片免费区| 久久综合国产亚洲精品| 国产探花在线观看一区二区| 亚洲综合色惰| 特大巨黑吊av在线直播| 午夜日本视频在线| 男插女下体视频免费在线播放| 人妻制服诱惑在线中文字幕| 我要看日韩黄色一级片| 青春草视频在线免费观看| 七月丁香在线播放| 亚洲久久久久久中文字幕| 亚洲自拍偷在线| 亚洲人与动物交配视频| 三级国产精品欧美在线观看| 国模一区二区三区四区视频| 三级经典国产精品| 久久久久久久精品精品| 18+在线观看网站| 亚洲欧洲日产国产| 日韩视频在线欧美| 国产精品一及| 久久精品国产鲁丝片午夜精品| 国产伦精品一区二区三区视频9| 免费黄网站久久成人精品| 九九久久精品国产亚洲av麻豆| 熟妇人妻不卡中文字幕| 最新中文字幕久久久久| 丝袜脚勾引网站| 亚洲色图av天堂| 国产免费视频播放在线视频| av专区在线播放| 99久久精品热视频| 国产在线一区二区三区精| 亚洲精品自拍成人| 精品久久久精品久久久| .国产精品久久| 欧美性猛交╳xxx乱大交人| 最近最新中文字幕大全电影3| 日韩精品有码人妻一区| 国产精品99久久99久久久不卡 | 国产免费一区二区三区四区乱码| 国产一区二区三区综合在线观看 | 免费av不卡在线播放| 久久99热这里只频精品6学生| kizo精华| 国产精品爽爽va在线观看网站| 成人一区二区视频在线观看| 精品午夜福利在线看| 亚洲天堂国产精品一区在线| 久久久久久久久大av| 亚洲天堂国产精品一区在线| 国产精品一及| 国产成人精品一,二区| 一级毛片 在线播放| av黄色大香蕉| 一区二区av电影网| 亚洲美女搞黄在线观看| 亚洲精品成人久久久久久| 天堂俺去俺来也www色官网| 一区二区三区精品91| 麻豆乱淫一区二区| 男的添女的下面高潮视频| 3wmmmm亚洲av在线观看| 国产视频内射| 日韩av在线免费看完整版不卡| 国产淫片久久久久久久久| 三级国产精品欧美在线观看| 一本一本综合久久| 日本黄大片高清| 五月天丁香电影| 国产免费福利视频在线观看| 亚洲精品国产av成人精品| 国产av码专区亚洲av| 国产在视频线精品| 在线天堂最新版资源| 男插女下体视频免费在线播放| 国产男女超爽视频在线观看| 亚洲国产精品成人综合色| 免费看光身美女| 精品少妇久久久久久888优播| 欧美激情国产日韩精品一区| 欧美日韩视频高清一区二区三区二| 熟女av电影| 午夜亚洲福利在线播放| 欧美激情国产日韩精品一区| 黄色怎么调成土黄色| 人人妻人人爽人人添夜夜欢视频 | 亚洲精品成人av观看孕妇| 亚洲av国产av综合av卡| 国产乱来视频区| 国产精品成人在线| 婷婷色av中文字幕| 免费观看a级毛片全部| 欧美成人精品欧美一级黄| 黄片wwwwww| 直男gayav资源| 爱豆传媒免费全集在线观看| 中文字幕久久专区| 亚洲精品国产成人久久av| 亚洲精品一二三| 久久人人爽av亚洲精品天堂 | 秋霞伦理黄片| 综合色av麻豆| 麻豆成人av视频| 日本爱情动作片www.在线观看| 天天一区二区日本电影三级| 亚洲,一卡二卡三卡| 色综合色国产| 欧美成人一区二区免费高清观看| 国产男女超爽视频在线观看| 亚洲va在线va天堂va国产| 国产午夜福利久久久久久| 亚洲最大成人手机在线| 纵有疾风起免费观看全集完整版| 国产美女午夜福利| 免费高清在线观看视频在线观看| 国产精品一二三区在线看| 少妇人妻一区二区三区视频| 嘟嘟电影网在线观看| 最新中文字幕久久久久| 精品亚洲乱码少妇综合久久|