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

    蘭州新區(qū)SCS徑流模型改進(jìn)與應(yīng)用研究

    2021-03-28 11:49:17張永明蒲秉華王鵬全
    人民黃河 2021年3期
    關(guān)鍵詞:產(chǎn)流蘭州片區(qū)

    張永明 蒲秉華 王鵬全

    摘 要:近年來(lái)全國(guó)氣候異常,由強(qiáng)降雨引發(fā)的洪澇災(zāi)害已成為城市經(jīng)濟(jì)快速發(fā)展面臨的嚴(yán)峻水問(wèn)題之一。因此,以無(wú)水文資料地區(qū)——蘭州新區(qū)為研究區(qū),開(kāi)展降雨徑流過(guò)程模擬研究是十分必要的。首先根據(jù)黃河干支流區(qū)綜合概化時(shí)程分配雨型表、《甘肅省暴雨洪水圖集》,利用地方經(jīng)驗(yàn)公式計(jì)算得到研究區(qū)不同頻率的產(chǎn)流量;其次借鑒相關(guān)研究成果,考慮地形地貌特征、下墊面情況,結(jié)合SCS徑流模型及DEM-GIS技術(shù),劃分流域片區(qū),以流域片區(qū)中各典型水文斷面為研究單元,構(gòu)建適合蘭州新區(qū)產(chǎn)流計(jì)算的模型參數(shù),并運(yùn)用LH-OAT方法對(duì)模型參數(shù)進(jìn)行敏感度分析;最后通過(guò)對(duì)比地方經(jīng)驗(yàn)公式法和SCS徑流模型的產(chǎn)流成果,經(jīng)科學(xué)分析后認(rèn)為,模型模擬精度較為理想,符合蘭州新區(qū)實(shí)際降雨-徑流過(guò)程,推薦研究區(qū)采用SCS徑流模型計(jì)算成果。

    關(guān)鍵詞:DEM-GIS技術(shù);SCS-CN模型改進(jìn);LH-OAT;無(wú)水文資料地區(qū);蘭州新區(qū)

    Abstract:In recent years, the national climate has been abnormal and floods caused by heavy rainfall have become one of the severe water problems facing the rapid economic development of cities. Therefore, it is very necessary to carry out rainfall-runoff process simulation study with the area without hydrological data, i.e. Lanzhou New District as the research area. Firstly, according to the comprehensive generalized time-history distribution rain pattern table of the main and tributary areas of the Yellow River and the “Atlas of Storm and Flood in Gansu Province”, the local empirical formula was used to calculate the flow rate of the study area at different frequencies; secondly, the relevant research results were used for reference and considering the topographic and geomorphic characteristics and underlying surface conditions, combined with SCS runoff model and DEM-GIS technology the basin area was divided, and the typical hydrological sections in the basin area were used as the research unit, the model parameters suitable for the calculation of runoff generation in Lanzhou New Area were constructed and the LH-OAT method was used to conduct sensitivity analysis on model parameters; finally, by comparing the local empirical formula method and the runoff results of the SCS runoff model, scientific analysis concluded that the model simulation accuracy is relatively ideal, which is consistent with the actual rainfall-runoff process in Lanzhou New Area. Therefore, it is recommended to use the SCS runoff model calculation results in the study area.

    Key words: DEM-GIS technology; SCS-CN model improvement; LH-OAT; area without hydrological data; Lanzhou New District

    “自然-社會(huì)”二元因素對(duì)水文循環(huán)規(guī)律造成了很大影響,尤其近幾年,隨著暴雨強(qiáng)度不斷增大,不僅導(dǎo)致城市部分城區(qū)雨水集聚,而且對(duì)郊區(qū)鄉(xiāng)鎮(zhèn)的群眾生命財(cái)產(chǎn)造成了不同程度的破壞,因此由強(qiáng)降雨引發(fā)的洪澇災(zāi)害對(duì)城市的影響也成為現(xiàn)今重要的研究課題之一。蘭州新區(qū)徑流的主要來(lái)源為流域上游的天然砂溝,洪水主要由暴雨造成,具有峰高、量集中、來(lái)勢(shì)迅猛、歷時(shí)短、沖刷破壞性強(qiáng)的特點(diǎn)[1],不僅對(duì)城市居民造成了較大的影響,而且對(duì)城市基礎(chǔ)設(shè)施構(gòu)成了一定的威脅,因此蘭州新區(qū)面臨的防洪問(wèn)題隨著時(shí)間的推移將會(huì)越發(fā)嚴(yán)重。目前蘭州新區(qū)沒(méi)有主要河道的歷史洪水資料,雨量觀測(cè)資料也十分缺乏,屬于無(wú)水文資料地區(qū)。

    SCS模型是美國(guó)農(nóng)業(yè)部水土保持局開(kāi)發(fā)的小流域設(shè)計(jì)洪水模型,計(jì)算過(guò)程簡(jiǎn)單,所需參數(shù)較少,資料易于獲取,尤其適用于資料缺乏地區(qū),并且有效地考慮了流域下墊面的特點(diǎn),是一種較好的計(jì)算小流域降雨徑流量的方法[2]。在我國(guó),許多專(zhuān)家和學(xué)者將此模型應(yīng)用于小流域工程規(guī)劃、流域水土保持、城市徑流、無(wú)資料流域的徑流模擬等研究中[3]。筆者通過(guò)改進(jìn)SCS模型,提出適合蘭州新區(qū)的模型參數(shù),以期為今后蘭州新區(qū)產(chǎn)匯流計(jì)算提供一定技術(shù)指導(dǎo)。

    1 研究區(qū)概況及資料

    1.1 研究區(qū)概況

    蘭州新區(qū)是全國(guó)第五個(gè)、西北地區(qū)第一個(gè)國(guó)家級(jí)新區(qū),位于隴西黃土高原的西北部,屬青藏高原、內(nèi)蒙古高原和黃土高原的交匯地,南北長(zhǎng)約73 km,東西寬約32 km,規(guī)劃總面積1 744 km2,地理坐標(biāo)為北緯36°17′15″—36°43′29″、東經(jīng)103°29′22″—103°49′46″[4]。蘭州新區(qū)深居內(nèi)陸,屬于大陸性冷溫帶半干旱氣候區(qū)。區(qū)內(nèi)無(wú)常年性天然地表徑流分布,主要溝谷有水阜河、咸水河、堿溝以及龔巴川,主要的水利設(shè)施有引大入秦工程、皋蘭縣西岔電力提灌工程西干渠和水庫(kù)。周邊低山及黃土丘陵植被很差,下墊面消納、蓄滯雨水能力較弱,降雨極易產(chǎn)生洪水。此外,蘭州新區(qū)所在流域并沒(méi)有主要河道的歷史洪水資料,雨量觀測(cè)資料缺乏,氣象信息監(jiān)測(cè)與自動(dòng)傳輸設(shè)施缺位,無(wú)法為暴雨洪水的實(shí)時(shí)監(jiān)測(cè)和有效預(yù)警提供服務(wù),再加上蘭州新區(qū)規(guī)劃范圍內(nèi)現(xiàn)狀防洪系統(tǒng)較為薄弱,若遭遇強(qiáng)暴雨量級(jí)的洪水,將面臨極大的洪水災(zāi)害。

    1.2 實(shí)測(cè)水文資料

    流域內(nèi)無(wú)降雨實(shí)測(cè)資料時(shí),可采用周邊氣象站點(diǎn)實(shí)測(cè)資料進(jìn)行插補(bǔ)。蘭州新區(qū)無(wú)任何實(shí)測(cè)徑流資料,產(chǎn)流計(jì)算依據(jù)附近永登、皋蘭、蘭州和中川機(jī)場(chǎng)等氣象站點(diǎn)暴雨資料以及《甘肅省暴雨洪水特性研究》(2003年)、《甘肅省暴雨洪水圖集》(1988年)、《甘肅省洪水調(diào)查》(第三冊(cè),黃河流域)、《甘肅省水文圖集》等資料。

    依據(jù)經(jīng)驗(yàn)頻率和暴雨統(tǒng)計(jì)參數(shù),按P-Ⅲ型理論頻率曲線分別對(duì)永登、皋蘭、蘭州、中川機(jī)場(chǎng)等4個(gè)雨量站最大24 h降雨系列進(jìn)行頻率分析計(jì)算,適線均值不做調(diào)整,偏態(tài)系數(shù)CS=3.5CV,求得各站年最大24 h暴雨量,見(jiàn)表1。

    2 基于DEM-GIS的水系分析

    運(yùn)用ArcGIS先后對(duì)研究區(qū)DEM數(shù)據(jù)進(jìn)行填洼、流向分析、匯流累計(jì)量分析、河網(wǎng)提取與分級(jí)、流域盆地分割與確定,最終得到天然情況下以流域盆地和矢量河網(wǎng)作為背景的集水流域計(jì)算數(shù)據(jù),并根據(jù)研究區(qū)土地利用規(guī)劃圖、蘭州新區(qū)水利建設(shè)規(guī)劃圖等資料,結(jié)合本次集水盆地分析,將蘭州新區(qū)劃分為西片區(qū)、北片區(qū)和盆地中心的城區(qū)(以不同底色加以區(qū)分),如圖1所示。因?yàn)樘m州新區(qū)徑流主要來(lái)源為流域上游的天然砂溝,洪水主要由暴雨造成,所以選取蘭州新區(qū)上游流域?yàn)檠芯繉?duì)象,即出山口斷面和正路鄉(xiāng)斷面以北區(qū)域,并提取了各片區(qū)集水面積等流域參數(shù),成果見(jiàn)表2。

    3 經(jīng)驗(yàn)公式產(chǎn)流計(jì)算

    利用2004版《甘肅省暴雨特性研究》中點(diǎn)暴雨參數(shù),結(jié)合流域降雨時(shí)—面—深關(guān)系、流域形狀系數(shù)曲線,將不同歷時(shí)的設(shè)計(jì)點(diǎn)暴雨量經(jīng)點(diǎn)面折減系數(shù)及形狀系數(shù)修正后,得到暴雨點(diǎn)面指數(shù),計(jì)算得到不同頻率的設(shè)計(jì)面暴雨量,見(jiàn)表3。由于永登、皋蘭、蘭州、中川機(jī)場(chǎng)4個(gè)雨量站年最大24 h雨量綜合適線結(jié)果與通過(guò)《甘肅省暴雨特性研究》等值線圖查算流域重心處24 h雨量差異較小,同時(shí)考慮到《甘肅省暴雨特性研究》中等值線位置以及暴雨參數(shù)的確定是精準(zhǔn)的、科學(xué)的、全面的,因此認(rèn)為暴雨設(shè)計(jì)是合理的。

    根據(jù)黃河干支流區(qū)綜合概化時(shí)程分配雨型表,結(jié)合《甘肅省暴雨洪水圖集》中各分區(qū)產(chǎn)流期平均入滲率,對(duì)北片區(qū)、西片區(qū)進(jìn)行產(chǎn)流計(jì)算,結(jié)果見(jiàn)表4(下文將地方經(jīng)驗(yàn)公式法得到的計(jì)算成果簡(jiǎn)稱(chēng)為對(duì)比產(chǎn)流量)。

    式中:P為降雨總量,mm;R為徑流深,mm;Ia為初損量,即徑流開(kāi)始前的損失量,主要包括蒸發(fā)、植被截留、地表填洼蓄水等,mm;S為流域當(dāng)前最大可滯留量,mm,它是后損量的上限;λ為初損系數(shù),是一個(gè)無(wú)量綱區(qū)域參數(shù),λ>0。

    由于S的變化范圍很大,不便于取值,因此引入無(wú)因次參數(shù)CN來(lái)推求S[6]。

    CN是反映降雨前流域特征的綜合參數(shù),反映流域前期土壤濕潤(rùn)程度、坡度、植被、土壤類(lèi)型和土地利用狀況的綜合特性,可以較好地反映下墊面條件對(duì)產(chǎn)匯流過(guò)程的影響[7]。

    實(shí)際計(jì)算過(guò)程中,降雨量P為己知量,因此CN值的確定是建立SCS產(chǎn)流模型的關(guān)鍵。《美國(guó)國(guó)家工程手冊(cè)》在第四章中詳細(xì)列出了CN值查算表,按照土壤前期濕潤(rùn)程度、土壤類(lèi)型和土地利用方式,可查算CN值[8]。

    5 模型的應(yīng)用

    5.1 參數(shù)λ、CN的敏感性分析

    為了定性了解參數(shù)λ、CN在蘭州新區(qū)西、北片區(qū)不同頻率下對(duì)產(chǎn)流模擬的影響,利用LH-OAT全局敏感度分析方法[9]對(duì)SCS模型進(jìn)行參數(shù)敏感性分析,分析成果如圖2所示。

    分析對(duì)比表明:在相同頻率不同水文控制斷面下,SCS模型的兩個(gè)參數(shù)中CN是最敏感參數(shù),且隨著頻率的變化,兩參數(shù)的敏感性始終未發(fā)生改變,即CN>λ;CN、λ值的敏感性隨頻率的不斷增大呈現(xiàn)遞增趨勢(shì);CN值的敏感性變化速率遠(yuǎn)大于λ值的,同時(shí)隨著頻率的不斷變大,參數(shù)對(duì)SCS模型產(chǎn)流結(jié)果影響變大。SCS模型在蘭州新區(qū)運(yùn)用過(guò)程中都是對(duì)CN值較λ值敏感,即CN值對(duì)SCS模型產(chǎn)流值影響較大。

    5.2 CN值

    在水文學(xué)中,CN值是一個(gè)無(wú)量綱參數(shù)[10]。SCS模型中參數(shù)CN使得最后的輸出結(jié)果對(duì)該參數(shù)高度敏感,CN值變化±10%就會(huì)引起徑流量計(jì)算結(jié)果分別變化55%、-45%[11],因此CN的取值對(duì)于準(zhǔn)確計(jì)算蘭州新區(qū)地表徑流量具有重要影響。利用MODIS衛(wèi)星遙感監(jiān)測(cè)土地利用類(lèi)型,提取了2013年蘭州新區(qū)的土地利用類(lèi)型,研究區(qū)主要土地利用類(lèi)型有灌木、草地、基本農(nóng)田、城鎮(zhèn)用地,見(jiàn)表5。根據(jù)Landsat 8衛(wèi)星觀測(cè)波段數(shù)據(jù),結(jié)合NDVI公式計(jì)算植被覆蓋程度,利用2016年7月31日Landsat 8影像,計(jì)算出了研究區(qū)各類(lèi)植被覆蓋面積占比,見(jiàn)表5。研究區(qū)土壤類(lèi)型主要為灰鈣土(669.21 km2)、栗鈣土(133.16 km2)[12],根據(jù)SCS模型土壤分類(lèi)定義表和土壤前期濕潤(rùn)程度等級(jí)劃分標(biāo)準(zhǔn)表,選定研究區(qū)土壤為B類(lèi)土壤,并根據(jù)蘭州新區(qū)降雨稀少、蒸發(fā)強(qiáng)烈、風(fēng)大沙多等情況,選定研究區(qū)濕潤(rùn)程度為Ⅱ類(lèi)。

    根據(jù)DEM-GIS技術(shù)提取的片區(qū)面積,考慮輸出成果對(duì)CN值的高度敏感性和參數(shù)本身的關(guān)鍵性,結(jié)合按不同植被覆蓋度查表得到的灌木、草地占研究片區(qū)面積的大小,運(yùn)用面積加權(quán)平均法,計(jì)算得到能反映蘭州新區(qū)不同土地利用類(lèi)型在各流域片區(qū)的平均CN值,分別為69.2、79.7。因?yàn)樘m州新區(qū)土壤類(lèi)型絕大多數(shù)屬灰鈣土,僅有少量栗鈣土等其他土壤,且裸地面積占流域面積比例相對(duì)較小,不會(huì)對(duì)最后結(jié)果造成很大影響,所以裸地仍采用查表所得CN值,即86。調(diào)查蘭州新區(qū)主要農(nóng)作物耕種情況發(fā)現(xiàn),作物主要分為直行種植、等高種植、等高耕種,查詢CN值查算表并求平均值后得到作物CN值,即76.3。這樣計(jì)算滿足了片區(qū)流域在空間上的差異,同時(shí)在一定程度上克服了降雨和下墊面條件在空間上的不均勻性對(duì)計(jì)算結(jié)果造成的影響。

    為了解西片區(qū)、北片區(qū)整體下墊面情況,預(yù)測(cè)蘭州新區(qū)今后徑流變化趨勢(shì),揭示蘭州新區(qū)土地利用變化條件下集水流域的降雨—徑流關(guān)系,通過(guò)面積加權(quán)平均法,計(jì)算得到了西片區(qū)、北片區(qū)綜合CN值,即77.3、77.9,見(jiàn)表6。

    5.3 初損系數(shù)λ

    初損量Ia受土地利用、耕作方式、灌溉條件、枝葉截留、下滲、填洼等因素的影響,它與土壤最大可滯留量S成正比關(guān)系,美國(guó)農(nóng)業(yè)部水土保持局在分析了大量試驗(yàn)結(jié)果的基礎(chǔ)上,提出了二者最合適的比值系數(shù)0.2[13],即初損系數(shù)λ為0.2。鄧景成等[14]研究分析了位于黃土區(qū)楊青川流域的λ值并提出:草地最適宜的λ值為0.13、裸地最適宜的λ值為0.03。首先,參數(shù)λ=0.2主要適用于美國(guó),導(dǎo)致SCS模型在研究區(qū)的應(yīng)用存在偏差和局限性;其次,我國(guó)地域遼闊,土壤土質(zhì)復(fù)雜多樣,特別在干旱區(qū)及半干旱區(qū)尤為明顯,相同土質(zhì)、不同地區(qū)的λ也可能存在較大差異,即λ=0.13、λ=0.03在研究區(qū)的直接使用會(huì)對(duì)模擬結(jié)果產(chǎn)生一定影響。因此,為使模型計(jì)算結(jié)果更加科學(xué)、可信、合理,必須對(duì)λ作進(jìn)一步分析計(jì)算。

    由于研究區(qū)植被主要為低矮稀疏的灌木草地,因此λ同CN的取值方法類(lèi)似,通過(guò)加權(quán)平均法計(jì)算得到:西片區(qū)灌木草地面積與裸地面積的比值為7.18%,λ為0.123 6;北片區(qū)灌木草地面積與裸地面積的比值為2.46%,λ為0.129 8。

    6 結(jié)果與分析

    采用暴雨歷時(shí)為24 h的降雨量,模擬的蘭州新區(qū)各片區(qū)產(chǎn)流量(徑流深)見(jiàn)表7。

    基于單因素方差分析方法[15],逐一分析模擬產(chǎn)流量和對(duì)比產(chǎn)流量在相同斷面、不同頻率下總體平均值差異的顯著性,方差分析結(jié)果見(jiàn)表8。在置信度95%情況下,正路鄉(xiāng)斷面、出山口斷面ANOVA檢驗(yàn)、Bonferroni檢驗(yàn)、Homogeneity檢驗(yàn)的Prob值均大于0.05,認(rèn)為模擬產(chǎn)出值減去對(duì)比值為0,充分說(shuō)明各斷面SCS模型模擬結(jié)果值和根據(jù)主雨峰外雨量同頻率放大所得到的結(jié)果值差異性不顯著,即SCS模型得到的產(chǎn)流結(jié)果基本符合蘭州新區(qū)實(shí)際降雨徑流情況。

    由于研究區(qū)為無(wú)資料地區(qū),缺乏長(zhǎng)系列水文資料佐證最后結(jié)果,因此為更加了解兩種產(chǎn)流計(jì)算結(jié)果的優(yōu)劣性,提高單因素方差分析結(jié)果的可信度,增強(qiáng)SCS模型輸出值的科學(xué)性,將模擬產(chǎn)流量、對(duì)比產(chǎn)流量(用流量表示)分別與降雨量擬定回歸方程,結(jié)果如圖3所示。分析可得:正路鄉(xiāng)斷面、出山口斷面降雨量—流量回歸方程的斜率和相關(guān)系數(shù)的平方R2均是模擬產(chǎn)流量?jī)?yōu)于對(duì)比產(chǎn)流量且相差較小,且模擬產(chǎn)流量小于對(duì)比產(chǎn)流量,可以認(rèn)為模擬產(chǎn)流值精度高于對(duì)比產(chǎn)流值精度,同時(shí)也間接證明了單因素方差分析法所得成果的準(zhǔn)確性。因此,SCS模型得到的產(chǎn)流量具有合理性、可靠性,模擬精度也較為理想。

    7 結(jié) 論

    運(yùn)用DEM-GIS技術(shù),根據(jù)蘭州新區(qū)土地利用類(lèi)型、植被覆蓋程度,將SCS模型應(yīng)用于蘭州新區(qū)進(jìn)行徑流模擬,得出的主要結(jié)論如下。

    (1)利用LH-OAT全局敏感度分析方法,以蘭州新區(qū)西、北兩片區(qū)為背景,對(duì)SCS模型中參數(shù)CN、λ進(jìn)行敏感性分析,發(fā)現(xiàn)CN敏感度遠(yuǎn)大于λ,且CN的敏感度變化速率也遠(yuǎn)大于λ的。

    (2)綜合考慮土地利用類(lèi)型、植被覆蓋程度,結(jié)合面積加權(quán)平均法推薦了適合蘭州新區(qū)不同土地利用類(lèi)型的CN值以及各流域片區(qū)的綜合CN值。西片區(qū)綜合CN值為77.9、λ值為0.123 6,北片區(qū)綜合CN值為77.3、λ值為0.129 8。同時(shí)推薦了正路鄉(xiāng)典型水文控制斷面,頻率為0.33%、1%、2%、5%、10%、20%、30%、40%、50%的產(chǎn)流量(流量)分別為30.5、28.4、17.1、11.0、7.1、3.7、2.1、1.3、0.7 m3/s;出山口典型水文控制斷面,頻率為0.33%、1%、2%、5%、10%、20%、30%、40%、50%的產(chǎn)流量(流量)分別為42.6、35.4、25.3、17.5、12.5、6.9、4.6、3.0、2.1 m3/s。

    (3)根據(jù)改進(jìn)后SCS模型產(chǎn)流成果,運(yùn)用單因素方差分析法及擬定線性回歸方程分析法,與通過(guò)主雨峰外雨量同頻率放大徑流量進(jìn)行分析對(duì)比,得出:改進(jìn)后的SCS模型能較好地反映蘭州新區(qū)實(shí)際降雨徑流過(guò)程,模型模擬精度較為理想。

    參考文獻(xiàn):

    [1] 劉卓,陳敏茹.濕陷性黃土地區(qū)泥石流溝道生態(tài)治理措施[J].甘肅水利水電技術(shù),2020,56(9):58-62.

    [2] 晉華,孫西歡,李仰斌.SCS模型在嵐河流域的應(yīng)用研究[J].太原理工大學(xué)學(xué)報(bào),2003,34(6):735-736,752.

    [3] 王愛(ài)娟,張平倉(cāng),丁文峰.應(yīng)用SCS模型計(jì)算秦巴山區(qū)小流域降雨徑流[J].人民長(zhǎng)江,2008,39(15):49-50,77,112.

    [4] 及志松.蘭州新區(qū)成為第五個(gè)國(guó)家級(jí)新區(qū)意義重大[J].國(guó)際融資,2012(10):33-35.

    [5] 李軍,劉昌明,王中根,等.現(xiàn)行普適降水入滲產(chǎn)流模型的比較研究:SCS與LCM[J].地理學(xué)報(bào),2014,69(7):926-932.

    [6] 張鵬,張雅莉,王剛,等.基于SCS模型的甘肅南部小流域徑流模擬研究[J].人民長(zhǎng)江,2013,44(增刊2):43-46.

    [7] 王瑾杰,丁建麗,張成,等.基于SCS模型的新疆博爾塔拉河流域徑流模擬[J].農(nóng)業(yè)工程學(xué)報(bào),2016,32(7):129-135.

    [8] 劉家福,蔣衛(wèi)國(guó),占文鳳,等.SCS模型及其研究進(jìn)展[J].水土保持研究,2010,17(2):120-124.

    [9] JUNG Y W, OH D S, KIM M, et al. Calibration of LEACHN Model Using LH-OAT Sensitivity Analysis[J].Nutrient Cycling in Agroecosystems,2010,87(2):261-275.

    [10] 王紅雷,王秀茹,王希,等.采用SCS-CN水文模型和GIS確定雨水集蓄工程的位置[J].農(nóng)業(yè)工程學(xué)報(bào),2012,28(22):108-114.

    [11] MISHRA S K, SINGH V P.Validity and Extension of the SCS-CN Method for Computing Infiltration and Rainfall-Excess Rates[J].Hydrological Processes,2004,18(17):3323-3345.

    [12] 石斌.基于GIS與RS的蘭州新區(qū)生態(tài)安全評(píng)價(jià)[D].蘭州:甘肅農(nóng)業(yè)大學(xué),2012:25-27.

    [13] 羅利芳,張科利,符素華.徑流曲線數(shù)法在黃土高原地表徑流量計(jì)算中的應(yīng)用[J].水土保持通報(bào),2002,22(3):58-61,68.

    [14] 鄧景成,高鵬,穆興民,等.模擬降雨條件下黃土區(qū)SCS模型的參數(shù)率定[J].水土保持研究,2018,25(5):205-210.

    [15] 戴金輝,袁靖.單因素方差分析與多元線性回歸分析檢驗(yàn)方法的比較[J].統(tǒng)計(jì)與決策,2016(9):23-26.

    【責(zé)任編輯 張 帥】

    猜你喜歡
    產(chǎn)流蘭州片區(qū)
    產(chǎn)流及其研究進(jìn)展
    不同坡面單元人工降雨產(chǎn)流試驗(yàn)與分析
    我的蘭州夢(mèng)
    黃河之聲(2021年8期)2021-07-23 03:34:32
    蘭州石化推進(jìn)改革正當(dāng)時(shí)
    片區(qū)綜合開(kāi)發(fā)的投資模式探討
    蘭州瑣記
    展望臨港新片區(qū)金融改革與創(chuàng)新
    青島新片區(qū)打造人才聚集高地
    商周刊(2019年19期)2019-10-17 03:23:36
    北京山區(qū)側(cè)柏林地坡面初始產(chǎn)流時(shí)間影響因素
    白塔寺片區(qū)地圖
    世界建筑(2019年7期)2019-08-07 07:57:52
    亚洲色图综合在线观看| 亚洲欧洲日产国产| 精品少妇久久久久久888优播| 亚洲精品久久午夜乱码| 美女视频免费永久观看网站| 国产精品久久久人人做人人爽| 国产免费现黄频在线看| 在线观看免费午夜福利视频| 亚洲精品av麻豆狂野| 动漫黄色视频在线观看| 日韩欧美一区二区三区在线观看 | 亚洲av成人不卡在线观看播放网| 久久精品人人爽人人爽视色| 精品久久久久久电影网| 丝袜美足系列| 欧美日韩亚洲高清精品| 日韩欧美一区视频在线观看| 成在线人永久免费视频| 久久久久久亚洲精品国产蜜桃av| 两性午夜刺激爽爽歪歪视频在线观看 | 天天影视国产精品| 成人亚洲精品一区在线观看| 女警被强在线播放| 丰满少妇做爰视频| 一区二区三区国产精品乱码| 不卡一级毛片| 精品人妻在线不人妻| 精品少妇黑人巨大在线播放| 99香蕉大伊视频| 青青草视频在线视频观看| 久久国产精品影院| 亚洲av美国av| 欧美成狂野欧美在线观看| 丰满迷人的少妇在线观看| 侵犯人妻中文字幕一二三四区| 国产不卡一卡二| 亚洲欧美一区二区三区黑人| √禁漫天堂资源中文www| 757午夜福利合集在线观看| 日韩成人在线观看一区二区三区| 69av精品久久久久久 | 黄片小视频在线播放| 一级,二级,三级黄色视频| 2018国产大陆天天弄谢| 在线 av 中文字幕| 国产精品久久久人人做人人爽| 亚洲国产欧美一区二区综合| 国产aⅴ精品一区二区三区波| av免费在线观看网站| av欧美777| 久久99一区二区三区| 亚洲第一欧美日韩一区二区三区 | av有码第一页| 国产成人一区二区三区免费视频网站| 香蕉久久夜色| av欧美777| 国产精品久久久久久人妻精品电影 | 国产精品99久久99久久久不卡| 日韩视频在线欧美| 久久精品91无色码中文字幕| 欧美久久黑人一区二区| 女人被躁到高潮嗷嗷叫费观| 亚洲av电影在线进入| 亚洲人成伊人成综合网2020| 欧美人与性动交α欧美软件| 午夜福利一区二区在线看| 一级片免费观看大全| 国产精品久久电影中文字幕 | 飞空精品影院首页| 天天添夜夜摸| av网站在线播放免费| 久久久久网色| 少妇被粗大的猛进出69影院| 亚洲av国产av综合av卡| 精品少妇久久久久久888优播| 悠悠久久av| 精品人妻熟女毛片av久久网站| 制服诱惑二区| 王馨瑶露胸无遮挡在线观看| 美女午夜性视频免费| 国产福利在线免费观看视频| 最新美女视频免费是黄的| 国产日韩欧美视频二区| 黄色视频不卡| 免费不卡黄色视频| 不卡av一区二区三区| 国产日韩欧美在线精品| 国产黄色免费在线视频| 啦啦啦在线免费观看视频4| 999久久久国产精品视频| 午夜福利在线免费观看网站| 菩萨蛮人人尽说江南好唐韦庄| 久久久久久亚洲精品国产蜜桃av| 色婷婷久久久亚洲欧美| 国产精品免费一区二区三区在线 | 色婷婷av一区二区三区视频| 久久精品亚洲熟妇少妇任你| 亚洲第一青青草原| 大香蕉久久成人网| 我的亚洲天堂| 精品人妻1区二区| 国产男女超爽视频在线观看| 久久久久国产一级毛片高清牌| 女性被躁到高潮视频| 午夜福利乱码中文字幕| 国产黄频视频在线观看| 9191精品国产免费久久| 精品一品国产午夜福利视频| tocl精华| 国产在视频线精品| 91麻豆精品激情在线观看国产 | 亚洲av欧美aⅴ国产| 色尼玛亚洲综合影院| 女同久久另类99精品国产91| 国产精品.久久久| 最新在线观看一区二区三区| 亚洲精品成人av观看孕妇| 一本一本久久a久久精品综合妖精| 亚洲天堂av无毛| 午夜福利视频精品| 一区福利在线观看| 最黄视频免费看| 午夜日韩欧美国产| 久热这里只有精品99| 久久精品亚洲av国产电影网| 热re99久久国产66热| 无遮挡黄片免费观看| 天天操日日干夜夜撸| 国产aⅴ精品一区二区三区波| 免费在线观看影片大全网站| 亚洲七黄色美女视频| svipshipincom国产片| 久久久久久免费高清国产稀缺| 在线天堂中文资源库| 亚洲伊人久久精品综合| 国产精品国产高清国产av | a在线观看视频网站| 精品国内亚洲2022精品成人 | 亚洲精品久久成人aⅴ小说| 十八禁人妻一区二区| 99精品在免费线老司机午夜| 啦啦啦免费观看视频1| 亚洲av国产av综合av卡| 黑人操中国人逼视频| 99re在线观看精品视频| 美女主播在线视频| 午夜福利影视在线免费观看| 高清视频免费观看一区二区| 精品久久蜜臀av无| 悠悠久久av| 久久人妻福利社区极品人妻图片| 久久国产精品人妻蜜桃| 亚洲欧美色中文字幕在线| 久久天躁狠狠躁夜夜2o2o| 91老司机精品| 国产日韩欧美在线精品| 男人操女人黄网站| 欧美激情久久久久久爽电影 | 69精品国产乱码久久久| 色视频在线一区二区三区| 悠悠久久av| 免费观看av网站的网址| 人人妻人人添人人爽欧美一区卜| 亚洲第一av免费看| av福利片在线| 国产精品久久久久久精品电影小说| 欧美黄色片欧美黄色片| 精品久久蜜臀av无| 午夜精品国产一区二区电影| 91av网站免费观看| 精品国产一区二区久久| 精品少妇久久久久久888优播| 12—13女人毛片做爰片一| 国产精品香港三级国产av潘金莲| 黑丝袜美女国产一区| 99久久国产精品久久久| 正在播放国产对白刺激| 9热在线视频观看99| av线在线观看网站| www.999成人在线观看| 亚洲五月婷婷丁香| 成人18禁在线播放| 99九九在线精品视频| 韩国精品一区二区三区| 日韩熟女老妇一区二区性免费视频| 欧美日韩视频精品一区| 国产av国产精品国产| 大型黄色视频在线免费观看| 十八禁网站网址无遮挡| 99香蕉大伊视频| 日日摸夜夜添夜夜添小说| 中文字幕制服av| 国产在线免费精品| 一区在线观看完整版| 天堂8中文在线网| 中文亚洲av片在线观看爽 | 免费在线观看完整版高清| 18在线观看网站| 亚洲情色 制服丝袜| 久热爱精品视频在线9| 在线天堂中文资源库| 亚洲国产欧美在线一区| 好男人电影高清在线观看| 大型av网站在线播放| 人人妻人人添人人爽欧美一区卜| 午夜福利乱码中文字幕| 欧美黄色淫秽网站| 久久久久久久大尺度免费视频| 日本a在线网址| 国产成人av激情在线播放| 国产亚洲精品第一综合不卡| 免费人妻精品一区二区三区视频| 色老头精品视频在线观看| 十八禁网站免费在线| 成人国产av品久久久| 国内毛片毛片毛片毛片毛片| 日本av手机在线免费观看| 亚洲第一av免费看| 日本撒尿小便嘘嘘汇集6| 中文亚洲av片在线观看爽 | 久久精品国产99精品国产亚洲性色 | 18禁观看日本| 最新的欧美精品一区二区| 亚洲欧美一区二区三区久久| 久久久久精品人妻al黑| av天堂久久9| videos熟女内射| xxxhd国产人妻xxx| 一区二区日韩欧美中文字幕| 极品人妻少妇av视频| 精品少妇一区二区三区视频日本电影| 性色av乱码一区二区三区2| 我的亚洲天堂| 91麻豆精品激情在线观看国产 | 国产国语露脸激情在线看| 水蜜桃什么品种好| 18禁美女被吸乳视频| 黄片播放在线免费| 午夜福利一区二区在线看| 日本撒尿小便嘘嘘汇集6| 成人特级黄色片久久久久久久 | 热re99久久国产66热| 另类亚洲欧美激情| 777米奇影视久久| 女人被躁到高潮嗷嗷叫费观| 国产黄频视频在线观看| 手机成人av网站| 欧美黑人欧美精品刺激| av一本久久久久| 涩涩av久久男人的天堂| 国产成人av教育| 国产视频一区二区在线看| 日韩欧美免费精品| 王馨瑶露胸无遮挡在线观看| 高清视频免费观看一区二区| 搡老乐熟女国产| tube8黄色片| 一区在线观看完整版| 黄色视频在线播放观看不卡| 51午夜福利影视在线观看| 久久热在线av| 欧美变态另类bdsm刘玥| 人妻久久中文字幕网| 自线自在国产av| 这个男人来自地球电影免费观看| 色老头精品视频在线观看| 90打野战视频偷拍视频| 国产99久久九九免费精品| 人妻 亚洲 视频| 日本撒尿小便嘘嘘汇集6| 蜜桃国产av成人99| 黄片大片在线免费观看| 又紧又爽又黄一区二区| 欧美在线一区亚洲| 热99国产精品久久久久久7| 国产三级黄色录像| 日本av手机在线免费观看| 久热这里只有精品99| 国产区一区二久久| av线在线观看网站| 99久久国产精品久久久| 午夜精品国产一区二区电影| 麻豆国产av国片精品| 在线观看66精品国产| 中文字幕人妻丝袜制服| 亚洲av成人一区二区三| 亚洲中文av在线| 最新美女视频免费是黄的| 99re在线观看精品视频| 久久久久国产一级毛片高清牌| 伊人久久大香线蕉亚洲五| 高清毛片免费观看视频网站 | 国产成人啪精品午夜网站| 国产亚洲精品第一综合不卡| 日韩有码中文字幕| 每晚都被弄得嗷嗷叫到高潮| tocl精华| 亚洲精品中文字幕在线视频| 国产黄频视频在线观看| 亚洲中文字幕日韩| a在线观看视频网站| 手机成人av网站| 免费女性裸体啪啪无遮挡网站| 久久国产亚洲av麻豆专区| 99精品欧美一区二区三区四区| 黑人猛操日本美女一级片| 日韩成人在线观看一区二区三区| 国产成人免费观看mmmm| 免费看a级黄色片| 日韩有码中文字幕| 18禁观看日本| bbb黄色大片| 国产一区二区三区视频了| 亚洲av国产av综合av卡| 精品少妇一区二区三区视频日本电影| 1024视频免费在线观看| 麻豆成人av在线观看| 欧美精品一区二区大全| www日本在线高清视频| 亚洲人成77777在线视频| 啦啦啦中文免费视频观看日本| 亚洲精品一二三| 色尼玛亚洲综合影院| 精品视频人人做人人爽| 美女扒开内裤让男人捅视频| 国产又色又爽无遮挡免费看| 国产精品免费大片| 久久久精品免费免费高清| av网站在线播放免费| 国产一区有黄有色的免费视频| 性少妇av在线| 精品人妻1区二区| 无人区码免费观看不卡 | 日韩中文字幕视频在线看片| 久久国产精品男人的天堂亚洲| 夜夜夜夜夜久久久久| 91字幕亚洲| 欧美人与性动交α欧美精品济南到| 日韩一区二区三区影片| 桃红色精品国产亚洲av| 国产亚洲一区二区精品| 久久精品亚洲av国产电影网| 老汉色∧v一级毛片| a在线观看视频网站| 麻豆成人av在线观看| 中亚洲国语对白在线视频| 9热在线视频观看99| 满18在线观看网站| 97人妻天天添夜夜摸| 国产成人精品无人区| 欧美精品啪啪一区二区三区| 亚洲一码二码三码区别大吗| 亚洲色图av天堂| 天天操日日干夜夜撸| 久久久久网色| 色精品久久人妻99蜜桃| 久久国产精品影院| 久久精品亚洲av国产电影网| 丝袜在线中文字幕| 视频区图区小说| 999久久久国产精品视频| 最黄视频免费看| 中文字幕制服av| 欧美日韩福利视频一区二区| 久久婷婷成人综合色麻豆| tocl精华| 91字幕亚洲| 黄色片一级片一级黄色片| 三上悠亚av全集在线观看| 欧美人与性动交α欧美精品济南到| 国产成人欧美| 无人区码免费观看不卡 | 国产亚洲午夜精品一区二区久久| 少妇被粗大的猛进出69影院| 两个人看的免费小视频| 99re6热这里在线精品视频| 黄片播放在线免费| 国产午夜精品久久久久久| 中文字幕av电影在线播放| av有码第一页| 一区二区三区乱码不卡18| 精品国产一区二区三区久久久樱花| 99九九在线精品视频| 极品教师在线免费播放| 老司机亚洲免费影院| 亚洲精品中文字幕一二三四区 | 一边摸一边做爽爽视频免费| 久久精品成人免费网站| 不卡一级毛片| 无人区码免费观看不卡 | 一级毛片电影观看| 在线永久观看黄色视频| 久久国产精品人妻蜜桃| 岛国毛片在线播放| 91九色精品人成在线观看| 日韩欧美免费精品| 欧美日韩av久久| 纵有疾风起免费观看全集完整版| 一级a爱视频在线免费观看| 中国美女看黄片| 建设人人有责人人尽责人人享有的| 超色免费av| 久久久久久免费高清国产稀缺| 国产av又大| 精品国产亚洲在线| 美女主播在线视频| 新久久久久国产一级毛片| 91精品三级在线观看| 黄色怎么调成土黄色| 欧美日本中文国产一区发布| 精品少妇内射三级| 性色av乱码一区二区三区2| 一边摸一边抽搐一进一出视频| 欧美精品人与动牲交sv欧美| 午夜两性在线视频| 精品熟女少妇八av免费久了| 涩涩av久久男人的天堂| 美女扒开内裤让男人捅视频| 日韩欧美三级三区| 黄色毛片三级朝国网站| 久久中文字幕一级| 水蜜桃什么品种好| 在线观看免费日韩欧美大片| 老司机靠b影院| 99久久国产精品久久久| 国产精品麻豆人妻色哟哟久久| 亚洲国产欧美在线一区| 国产一卡二卡三卡精品| 久久青草综合色| 国产精品一区二区在线观看99| 又黄又粗又硬又大视频| 午夜福利在线观看吧| 黄网站色视频无遮挡免费观看| 国产av又大| 自线自在国产av| 欧美黄色淫秽网站| 国产真人三级小视频在线观看| 熟女少妇亚洲综合色aaa.| 蜜桃在线观看..| 国产99久久九九免费精品| 老司机影院毛片| www日本在线高清视频| 中国美女看黄片| 婷婷成人精品国产| 日本撒尿小便嘘嘘汇集6| 日本欧美视频一区| 国产高清国产精品国产三级| 亚洲中文av在线| 又黄又粗又硬又大视频| 亚洲av国产av综合av卡| 日韩免费av在线播放| 亚洲伊人久久精品综合| 在线观看免费视频网站a站| 国产一区二区 视频在线| 国产精品.久久久| 中文字幕精品免费在线观看视频| 一级,二级,三级黄色视频| 69av精品久久久久久 | 久久午夜亚洲精品久久| 久久精品国产a三级三级三级| 亚洲第一欧美日韩一区二区三区 | 亚洲国产欧美在线一区| 精品乱码久久久久久99久播| 91精品国产国语对白视频| 女性生殖器流出的白浆| 欧美激情久久久久久爽电影 | 99久久人妻综合| 自拍欧美九色日韩亚洲蝌蚪91| 国产不卡av网站在线观看| xxxhd国产人妻xxx| 肉色欧美久久久久久久蜜桃| 日韩一区二区三区影片| 精品国产一区二区久久| 午夜福利乱码中文字幕| 啦啦啦中文免费视频观看日本| 亚洲九九香蕉| 黄色成人免费大全| 操出白浆在线播放| 日韩 欧美 亚洲 中文字幕| 久久影院123| 亚洲成av片中文字幕在线观看| 亚洲成人手机| 电影成人av| 每晚都被弄得嗷嗷叫到高潮| 少妇的丰满在线观看| 国产不卡av网站在线观看| 免费一级毛片在线播放高清视频 | 丰满饥渴人妻一区二区三| 一进一出抽搐动态| 国产精品香港三级国产av潘金莲| 怎么达到女性高潮| 肉色欧美久久久久久久蜜桃| 亚洲中文字幕日韩| 极品人妻少妇av视频| 久久精品aⅴ一区二区三区四区| 怎么达到女性高潮| 一夜夜www| 成人特级黄色片久久久久久久 | 悠悠久久av| 久久精品国产99精品国产亚洲性色 | 亚洲国产精品一区二区三区在线| 亚洲中文字幕日韩| 操美女的视频在线观看| 少妇被粗大的猛进出69影院| 国产成人精品久久二区二区免费| a级毛片在线看网站| 成人av一区二区三区在线看| 日韩中文字幕视频在线看片| 天堂动漫精品| 不卡av一区二区三区| 国产精品成人在线| 美女大奶头视频| 黑人欧美特级aaaaaa片| 欧美性猛交╳xxx乱大交人| 亚洲欧美精品综合一区二区三区| 黄色片一级片一级黄色片| avwww免费| 免费观看精品视频网站| 亚洲av第一区精品v没综合| 国产极品精品免费视频能看的| 亚洲欧美激情综合另类| 91麻豆精品激情在线观看国产| 亚洲国产欧美网| 母亲3免费完整高清在线观看| 国内少妇人妻偷人精品xxx网站 | 99在线人妻在线中文字幕| 不卡av一区二区三区| 久久午夜亚洲精品久久| 18禁黄网站禁片免费观看直播| 亚洲在线观看片| 国产黄色小视频在线观看| 精品国产亚洲在线| 又爽又黄无遮挡网站| 99精品欧美一区二区三区四区| 熟妇人妻久久中文字幕3abv| 亚洲色图 男人天堂 中文字幕| 91av网一区二区| 怎么达到女性高潮| 99国产精品一区二区蜜桃av| 一进一出抽搐gif免费好疼| 欧美另类亚洲清纯唯美| 国产黄a三级三级三级人| 欧美日韩瑟瑟在线播放| 成人av在线播放网站| 男人舔女人下体高潮全视频| 丁香六月欧美| 日本 欧美在线| 久久欧美精品欧美久久欧美| 天堂影院成人在线观看| 亚洲九九香蕉| 亚洲中文字幕一区二区三区有码在线看 | 久久国产精品影院| 很黄的视频免费| 一进一出好大好爽视频| 欧美日韩国产亚洲二区| tocl精华| 一区福利在线观看| 在线国产一区二区在线| 又大又爽又粗| 国产三级在线视频| 99热这里只有是精品50| 成人午夜高清在线视频| 亚洲 国产 在线| 国产人伦9x9x在线观看| 免费观看人在逋| 老汉色∧v一级毛片| 校园春色视频在线观看| 90打野战视频偷拍视频| 中文亚洲av片在线观看爽| 久久精品国产99精品国产亚洲性色| 亚洲精品美女久久久久99蜜臀| 亚洲欧美日韩东京热| 欧美日韩一级在线毛片| 精品国产三级普通话版| 国产亚洲精品av在线| 国产精华一区二区三区| 国产一区在线观看成人免费| 精品久久久久久久人妻蜜臀av| 18美女黄网站色大片免费观看| 男女午夜视频在线观看| 久久久精品大字幕| 啪啪无遮挡十八禁网站| 99国产极品粉嫩在线观看| 国产aⅴ精品一区二区三区波| 日本免费a在线| 欧美一区二区精品小视频在线| 在线播放国产精品三级| 国产精品女同一区二区软件 | 美女高潮的动态| АⅤ资源中文在线天堂| 网址你懂的国产日韩在线| 免费在线观看成人毛片| 男女床上黄色一级片免费看| 白带黄色成豆腐渣| av片东京热男人的天堂| 长腿黑丝高跟| 国产精品1区2区在线观看.| 桃色一区二区三区在线观看| 亚洲av成人av| 免费av不卡在线播放| 国产精品亚洲av一区麻豆| 天堂√8在线中文| 三级男女做爰猛烈吃奶摸视频| 国产69精品久久久久777片 | 亚洲欧美日韩东京热| 三级毛片av免费| 国产成年人精品一区二区| 国产乱人伦免费视频| 欧美黄色淫秽网站| 亚洲中文av在线| 欧美三级亚洲精品| 一卡2卡三卡四卡精品乱码亚洲| 成在线人永久免费视频| 12—13女人毛片做爰片一| av黄色大香蕉|