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

    廣西南流江流域降雨侵蝕力時(shí)空變化特征*

    2020-08-02 06:34:14陳國清趙銀軍莫德麗謝余初
    廣西科學(xué) 2020年3期
    關(guān)鍵詞:侵蝕性輸沙量雨量

    陳國清,趙銀軍**,莫德麗,梁 珊,謝余初,童 凱

    (1.南寧師范大學(xué),北部灣環(huán)境演變與資源利用教育部重點(diǎn)實(shí)驗(yàn)室/廣西地表過程與智能模擬重點(diǎn)實(shí)驗(yàn)室,廣西南寧 530001;2.南寧師范大學(xué),地理科學(xué)與規(guī)劃學(xué)院,廣西南寧 530001)

    0 引言

    降雨是導(dǎo)致土壤侵蝕的主要因素之一[1-2],降雨侵蝕力(Rainfall Erosivity,簡稱R)是指降雨過程引起地表土壤侵蝕的潛在能力,能客觀表征降雨量大小、降雨強(qiáng)度、降雨動(dòng)能和降雨歷時(shí)等對土壤侵蝕的綜合影響程度[1-3]。估算和分析降雨侵蝕力及其時(shí)空分布格局,對區(qū)域水土保持現(xiàn)狀評價(jià)、預(yù)報(bào)、治理和規(guī)劃具有重要意義[4-5]。目前,定量計(jì)算與表征降雨侵蝕力的經(jīng)典方法是Wischmeier 等[2]提出的基于降雨動(dòng)能(E)和降雨強(qiáng)度(I)乘積的算法,但由于長序列的短時(shí)尺度(秒、分或小時(shí))雨量信息等降雨動(dòng)能數(shù)據(jù)難以獲取,且計(jì)算繁瑣費(fèi)時(shí),使得該方法在一般地區(qū)難以直接應(yīng)用[5-7]。而各氣象站的雨量資料較容易獲取,基于日、月、年尺度雨量資料的降雨侵蝕力簡易算法應(yīng)運(yùn)而生,并得到實(shí)踐與應(yīng)用[6-7]。相對于年或月雨量,日雨量能夠提供更豐富的降雨特征信息,其模型算法估算精度相對較高[7-8],不僅被全國第一次水利普查采用[9],而且適用于我國大部分地區(qū),在各地降雨侵蝕力研究中得到較為廣泛的應(yīng)用,如珠江流域[7]、怒江[10]、四川省[11]、南方紅壤區(qū)[12]、曲靖市[13]等。

    南流江流域是我國西部大開發(fā)三大經(jīng)濟(jì)區(qū)之一的廣西北部灣經(jīng)濟(jì)區(qū)的重要組成區(qū)域,也是廣西典型的紅壤和亞熱帶生態(tài)系統(tǒng)水土保持重要區(qū)域之一。近年來,隨著全球環(huán)境變化以及北部灣經(jīng)濟(jì)區(qū)的不斷發(fā)展、森林轉(zhuǎn)型(林權(quán)改革)等影響,土地覆被變化強(qiáng)烈,流域水土流失防控工作形勢嚴(yán)峻[14-16]。流域內(nèi)降水空間分異大,降雨尤其是臺(tái)風(fēng)帶來的暴雨等是該區(qū)域水土流失的主要外動(dòng)力因素之一[15-17]。為此,本研究利用南流江流域內(nèi)氣象站點(diǎn)1961—2006年間日降雨量資料和下游泥沙觀測數(shù)據(jù),運(yùn)用日降雨量侵蝕力模型估算和分析流域內(nèi)40多年平均降雨侵蝕力及其時(shí)空變化特征與規(guī)律,并討論研究區(qū)侵蝕性降雨標(biāo)準(zhǔn),旨在為廣西沿海紅壤區(qū)土壤侵蝕監(jiān)測與水土流失防治及規(guī)劃工作等提供科學(xué)依據(jù)。

    1 材料與方法

    1.1 研究區(qū)概況

    南流江流域(21°34′—22°52′N,108°51′—110°22′E)地處廣西東南沿海,發(fā)源于廣西北流市大容山,干流長約287 km,流域面積約9 337 km2。流域東部、西部和北部地勢相對較高,中部和南部地勢偏低,上游有玉林盆地,中游是博白盆地和丘陵,下游是合浦沖積平原[17-18]。南流江流域多年平均徑流量為166 m3/s,水力資源較為豐富,是廣西沿海紅壤和磚紅壤區(qū)的最大入海河流。屬于亞熱帶季風(fēng)氣候區(qū),夏季高溫多雨,冬季溫涼少雨,年平均氣溫21—23℃,雨量豐沛,大雨暴雨多發(fā),年平均降雨量達(dá)1 400—1 960 mm,豐水期集中在4—10月[14-16]。主要土壤類型為紅壤、磚紅壤、棕壤、赤紅壤、黃棕壤、潮土、褐土、水稻土等,植被以亞熱帶常綠闊葉林為主[14-16]。

    1.2 數(shù)據(jù)來源與處理

    南流江流域8個(gè)氣象站點(diǎn),1961—2006年逐日降水?dāng)?shù)據(jù)和下游常樂水文站輸沙量觀測數(shù)據(jù),來自于中華人民共和國水文年鑒珠江流域水文資料。流域內(nèi)氣象站點(diǎn)分布均勻,且各站點(diǎn)降雨量觀測數(shù)據(jù)時(shí)間跨度均大于40 a,同時(shí)單個(gè)站點(diǎn)缺失數(shù)據(jù)均未超過10%。在數(shù)據(jù)處理時(shí),由于歷史原因浦北站在1962—1966年的數(shù)據(jù)記錄不完整,因此在計(jì)算流域1962—1966年降雨侵蝕力時(shí)剔除浦北站。

    1.3 方法

    1.3.1 降雨侵蝕力模型

    采用第一次全國水利普查水土保持專項(xiàng)普查使用的日降雨侵蝕力估算方程,該方法在中國南方降水豐沛的地區(qū)表現(xiàn)較好、計(jì)算精度相對較高[7-8,12-13]。具體計(jì)算公式如下:

    (1)

    β=0.8363+18.144/Pd12+24.455/Py12,

    (2)

    α=21.586β-7.1891,

    (3)

    式中:Mi為第i個(gè)半月的降雨侵蝕力,MJ·mm·hm-2·h-1·a-1;k為半月的天數(shù),Dj為半月內(nèi)第j天的侵蝕性日雨量;α、β為模型參數(shù),根據(jù)區(qū)域的降水量確定。本文中,半月時(shí)段的劃分方法是以每月的1—15天作為第一個(gè)半月,余下的天數(shù)作為第二個(gè)半月。Pd12是日雨量大于12 mm的日平均雨量,Py12是日雨量大于12 mm的年平均雨量,根據(jù)南流江流域特征選取我國南方較常用的日降雨量≥20 mm 作為閾值[4,19-20]。

    1.3.2 降雨侵蝕力時(shí)空特征

    南流江流域降雨侵蝕力的時(shí)空變化主要采取趨勢分析、周期分析和克里金(Kriging) 插值法。即利用線性回歸方程分析流域降雨侵蝕力的時(shí)間趨勢變化[21];然后,采用小波分析(Wavelet Analysis)方法計(jì)算分析降雨侵蝕力的時(shí)間序列周期變化特征,具體方法原理與步驟詳見文獻(xiàn)[22]。同時(shí),利用ArcGIS空間分析克里金插值方法對流域降雨量和降雨侵蝕力進(jìn)行空間內(nèi)插,得到連續(xù)分布的流域內(nèi)多年平均降雨侵蝕力空間分布圖。

    2 結(jié)果與分析

    2.1 南流江流域降雨侵蝕力的時(shí)間變化

    2.1.1 降雨侵蝕力年際變化

    首先,1961—2006年間南流江流域平均降雨侵蝕力為13 935.5 MJ·mm·hm-2·h-1·a-1,年降雨侵蝕力呈現(xiàn)波動(dòng)變化、輕微增加的趨勢。其中,南流江流域年降雨侵蝕力在1961—1964年間表現(xiàn)為不明顯的上升趨勢;在1965—1978年間則表現(xiàn)出波動(dòng)下降趨勢;之后在1979—1991年間也呈現(xiàn)波動(dòng)下降趨勢。在1992—2006年間,南流江流域降雨侵蝕力則表現(xiàn)出較大的波動(dòng)變化,振幅高達(dá)7 802 .2 MJ·mm·hm-2·h-1·a-1,年際變化大(圖1a)。其次,南流江流域年降雨侵蝕力(圖1a)與年降雨量(圖1b)變化趨勢基本相似。其中,流域內(nèi)年降雨侵蝕力和年降雨量較低的峰值出現(xiàn)年份相一致,兩者均出現(xiàn)在1962年、1989年;而年降雨侵蝕力和年降雨量最高峰值出現(xiàn)年份略有輕微差異,年降雨侵蝕力較高的峰值出現(xiàn)在1981年和2006年,分別為22 856.0 MJ·mm·hm-2·h-1·a-1和21 737.8 MJ·mm·hm-2·h-1·a-1,年降雨量的最大值出現(xiàn)的年份則是1981年(2 423.6 mm)和2002年(2 320.2 mm)。

    圖1 1961—2006年南流江流域年降雨侵蝕力(a)和年降雨量(b)的變化Fig.1 Changes of annual rainfall erosivity (a) and annual rainfall (b) in Nanliujiang River Basin during 1961—2006

    運(yùn)用小波分析對南流江流域46 a降雨侵蝕力進(jìn)行分析發(fā)現(xiàn)(圖2a),流域年降雨侵蝕力存在3個(gè)時(shí)間尺度的周期變化規(guī)律,第1主周期為23 a,第2主周期為7 a,第3主周期為12 a,其中第3主周期在小波方差曲線的峰值不明顯,降雨侵蝕力變化周期約14 a (圖2b)。

    圖2 小波分析的年降雨侵蝕力動(dòng)態(tài)變化(a)和年降雨侵蝕力第1主周期的更替情況(b)Fig.2 Dynamic change of annual rainfall erosivity by wavelet analysis (a) and replacement of the first main period of annual rainfall erosivity (b)

    2.1.2 降雨侵蝕力季節(jié)變化

    南流江流域降雨侵蝕力在各季節(jié)變化趨勢略有不同,春季呈現(xiàn)不顯著的降低趨勢,夏季、秋季和冬季呈現(xiàn)不明顯的上升趨勢(圖3),上升幅度依次為夏季>冬季>秋季。1961—2006年流域內(nèi)春季的平均降雨侵蝕力為3 092.1 MJ·mm·hm-2·h-1·a-1,占年均降雨侵蝕力的24.1%(圖3a)。夏季降雨侵蝕力變化趨勢與流域年際降雨侵蝕力變化最相似,平均降雨侵蝕力為7 326.1 MJ·mm·hm-2·h-1·a-1,占年均降雨侵蝕力的57.0%(圖3b)。秋季和冬季的平均降雨侵蝕力則分別占年均降雨侵蝕力的15.1% 和3.8%。

    圖3 1961—2006年南流江流域不同季節(jié)降雨侵蝕力的變化曲線Fig.3 Change curve of rainfall erosivity in different seasons in Nanliujiang River Basin during 1961—2006

    2.2 南流江流域年均降雨侵蝕力的空間分布

    在空間上,南流江流域降雨侵蝕力呈現(xiàn)南高北低的格局,沿流域下游到上游逐漸減弱。降雨侵蝕力高值區(qū)(≥17 000 MJ·mm·hm-2·h-1·a-1)主要分布在沿海地區(qū);低值區(qū)(≤12 000 MJ·mm·hm-2·h-1·a-1)則分布在上游內(nèi)陸等區(qū)域。分縣區(qū)統(tǒng)計(jì)發(fā)現(xiàn),欽南區(qū)和合浦縣降雨侵蝕力及年均侵蝕性降雨最高,博白縣和浦北縣次之,北流市、興業(yè)縣和玉州區(qū)最小。從流域降雨侵蝕力等值線與降雨量等值分布的情況上發(fā)現(xiàn),南流江流域年均降雨侵蝕力在10 000 MJ·mm·hm-2·h-1·a-1等值線上,與年均降雨量為1 000 mm等值線相鄰近;在14 000 MJ·mm·hm-2·h-1·a-1等值線上與年均降雨量為1 600 mm等值線相鄰近;在19 000 MJ·mm·hm-2·h-1·a-1等值線上,與年均降雨量為2 500 mm等值線相鄰近(圖4)。同時(shí),降雨侵蝕力隨著降雨量的增大而增大,但降雨侵蝕力的增長速率比降雨量的增長速率高。

    圖4 1961-2006年南流江流域年均降雨侵蝕力(a)和年均降雨量(b)的空間分布圖Fig.4 Spatial distribution of annual rainfall erosivity (a) and annual rainfall (b) in Nanliujiang River Basin during 1961-2006

    從年代尺度上來看(表1),近46年來,流域內(nèi)降雨侵蝕力在不同年代間略有差異,在偏中高值區(qū)占的流域面積呈現(xiàn)上升趨勢。1981—1990年流域內(nèi)降雨侵蝕力數(shù)值最高,達(dá)22 856.0 MJ·mm·hm-2·h-1·a-1,1971—1980年的降雨侵蝕力數(shù)值最低,為16 837.1 MJ·mm·hm-2·h-1·a-1。從年均值上看,2001—2006年的流域年均降雨侵蝕力最大(16 381.3 MJ·mm·hm-2·h-1·a-1),1971—1980年的年均降雨侵蝕力最小(13 132.3 MJ·mm·hm-2·h-1·a-1)。從降雨侵蝕力中高值區(qū)(≥15 000.0 MJ·mm·hm-2·h-1·a-1)所占的面積上來看,1981—1990年的中高值區(qū)的面積最大,約占流域總面積的36.7%;1961—1970年的最少,約占22.0%。在整個(gè)研究期間,中高值區(qū)以上的降雨侵蝕力所在的流域面積比重超過22.1%以上,且呈現(xiàn)緩慢增大的趨勢,表明南流江流域中高值段降雨侵蝕力的區(qū)域可能存在較嚴(yán)重的水土流失發(fā)展態(tài)勢。

    表1 1961—2006年各時(shí)間段的降雨侵蝕力差異比較Table 1 Comparison of rainfall erosivity in different periods during 1961—2006

    流域內(nèi)各雨量觀測站點(diǎn)觀測數(shù)據(jù)顯示(未作圖表),1961—2006年間,玉林站、北流站和靈山站的年降雨侵蝕力呈現(xiàn)輕微下降趨勢,陸川站、浦北站的年降雨侵蝕力變化不大,趨于平穩(wěn)水平狀態(tài)。常樂站的年降雨侵蝕力則呈現(xiàn)明顯的下降趨勢,這可能與南流江上中游暴雨發(fā)生頻數(shù)相對下游較少以及2000年以后退耕還林有關(guān)。博白站年降雨侵蝕力呈現(xiàn)輕微上升趨勢,合浦站年降雨侵蝕力呈現(xiàn)明顯的上升趨勢,這可能是合浦位于南流江下游臨海口易受臺(tái)風(fēng)、暴雨和海潮等降雨事件的影響,且來水量大,尤其是在洪水暴雨頻發(fā)的6—9月。由此可見,地形和強(qiáng)降雨事件等對區(qū)域降雨侵蝕力及其水土流失影響較大。

    同時(shí),由于2007—2009年間個(gè)別月份逐日降雨數(shù)據(jù)的缺失,本文主要分析過去1961—2006年南流江流域年降雨侵蝕力,但通過將1961—2006年與2010-2016年年均降雨侵蝕力對比分析(圖4a和圖5),研究結(jié)果能較客觀地反映流域降雨侵蝕力的變化情況。由圖4a和圖5可知,1961—2006年流域年均降雨侵蝕力在空間分布與2010—2016年的相似,兩者均呈現(xiàn)自上游到下游遞增的趨勢。1961—2006年間流域平均降雨侵蝕力為13 945.5 MJ·mm·hm-2·h-1·a-1,年降雨侵蝕力呈現(xiàn)波動(dòng)變化及輕微增加的趨勢,這與2010—2016年南流江流域平均降雨侵蝕力為15 491.1 MJ·mm·hm-2·h-1·a-1,呈現(xiàn)緩慢上升的趨勢相似??梢?,盡管本研究資料偏舊,但仍較好地反映流域降雨侵蝕力的時(shí)空變化特征,對當(dāng)前南流江流域降雨侵蝕及水土流失防治工作也具有一定的科學(xué)意義。

    圖5 2010—2016年南流江流域年均降雨侵蝕力空間分布圖Fig.5 Spatial distribution of annual rainfall erosivity in Nanliujiang River Basin during 2010—2016

    2.3 南流江流域降雨侵蝕力與下游輸沙量的相關(guān)性

    以南流江流域日降雨量≥10 mm,≥12 mm,≥15 mm,≥20 mm,≥25 mm,≥30 mm,≥35 mm,≥40 mm,≥50 mm分別作為9種不同級(jí)別雨量標(biāo)準(zhǔn),結(jié)合下游常樂水文站1961—2006年觀測的輸沙量,計(jì)算分析不同級(jí)別雨量下降雨侵蝕力與流域輸沙量的相關(guān)性。研究表明,侵蝕性雨量標(biāo)準(zhǔn)為10—20 mm/d 時(shí),日雨量越大,降雨侵蝕力與河流輸沙量的相關(guān)系數(shù)越強(qiáng);當(dāng)侵蝕性雨量標(biāo)準(zhǔn)為20—50 mm時(shí),其相關(guān)系數(shù)隨日雨量的增加而呈下降減弱趨勢。選取日降雨量≥20 mm作為侵蝕性雨量標(biāo)準(zhǔn)時(shí),南流江流域降雨侵蝕力與下游輸沙量的相關(guān)系數(shù)相對較高,R2約為0.773;選取日降雨量≥25 mm時(shí),兩者相關(guān)性系數(shù)次之(約為0.752)。當(dāng)選取日降雨量≥12 mm(全國侵蝕性降雨標(biāo)準(zhǔn))計(jì)算時(shí),流域降雨侵蝕力與下游輸沙量的相關(guān)系數(shù)為0.741(圖6)。其次,從南流江流域降雨侵蝕力(采用日降雨量≥20 mm作為侵蝕性雨量標(biāo)準(zhǔn)計(jì)算)與常樂水文站輸沙量的變化曲線對比來看(圖7),降雨侵蝕力的變化與常樂站輸沙量的變化曲線的趨勢基本保持一致,土壤侵蝕除受降雨量、降水強(qiáng)度、降雨時(shí)間等因素關(guān)系密切外,還受植被覆蓋度變化及坡度坡長等因素影響,降雨直接擊濺侵蝕地面,雨量達(dá)到一定程度后形成地表徑流沖刷,降雨侵蝕力的變化在一定程度上影響著水土流失,造成流域的輸沙量的變化。因此,選取日降雨量≥20 mm/d作為南流江流域侵蝕性降雨標(biāo)準(zhǔn)是比較合理的。

    圖6 不同侵蝕性雨量標(biāo)準(zhǔn)下降雨侵蝕力與河流輸沙量的相關(guān)系數(shù)Fig.6 Correlation coefficient between rainfall erosivity and river sediment transport under different erosive rainfall standards

    圖7 1961—2006年南流江流域降雨侵蝕力與常樂站輸沙量的變化曲線Fig.7 Change curve of rainfall erosivity and sediment discharge of Changle station in Nanliujiang River Basin during 1961—2006

    3 討論

    1)1961—2006年間南流江流域降雨侵蝕力空間分布格局變化不大,總體呈現(xiàn)由下游臨海向上游玉林市和北流市等區(qū)域逐漸減弱的格局,且合浦和常樂站點(diǎn)年降雨侵蝕力變化較大。時(shí)間上,流域年降雨侵蝕力總體呈輕微增大的趨勢,這與其他學(xué)者研究廣西[15-16]或涉及廣西區(qū)域的(如珠江流域[7])降雨侵蝕力的結(jié)果相似。劉斌濤等[4]在對中國大陸1960—2009年降雨侵蝕力的分析中,也指出廣西沿海地區(qū)的降雨侵蝕力呈桂東南-桂南向桂西北-桂西不顯著變化的遞減趨勢,與本文研究結(jié)果相似。其次,流域內(nèi)年降雨侵蝕力與降雨量空間格局相似,但年降雨侵蝕力年際變化比年降雨量更劇烈[23],這可能與我國南方降雨量變化特征及夏秋時(shí)節(jié)強(qiáng)降水事件密切相關(guān)[24]。同時(shí),南流江流域?qū)儆趶V西人工速生桉樹種植集中區(qū),其采伐周期為5—6 a,在桉樹采伐期和種植幼苗期,強(qiáng)降雨事件更容易沖刷地表土壤,不僅使該地區(qū)可能存在較大的水土流失潛在風(fēng)險(xiǎn)[16],而且增大了區(qū)域降雨侵蝕力觀測的不確定性。因此,把握好侵蝕性降雨標(biāo)準(zhǔn)并精準(zhǔn)計(jì)算降雨侵蝕力對區(qū)域土壤侵蝕防治工作具有重要的現(xiàn)實(shí)意義。

    2)目前,我國學(xué)者在應(yīng)用日降雨量方法計(jì)算降雨侵蝕力時(shí),通常采用全國性的12 mm作為侵蝕性降雨標(biāo)準(zhǔn),然而由于我國地形地貌復(fù)雜,不同區(qū)域間土壤、植被、降雨及其降雨歷時(shí)和強(qiáng)度的不同,各不同地區(qū)侵蝕性降雨標(biāo)準(zhǔn)也存在區(qū)域的差異性或者不同程度上的不一致性[25-26]。在我國南方地區(qū),劉斌濤等[4]和王萬中等[20]認(rèn)為與全國通用的侵蝕性降雨標(biāo)準(zhǔn)12 mm/d相比較,日降雨量≥20 mm/d 對南方地區(qū)年降雨侵蝕力影響較大。盧程隆等[19]在福建東南地區(qū)實(shí)驗(yàn)分析也指出侵蝕性降雨約為次降雨量≥20 mm;林春蕾等[27]指出云南玉溪紅壤區(qū)侵蝕性雨量標(biāo)準(zhǔn)約為15—25 mm。黃俊等[28]在廣東五華縣做土壤侵蝕實(shí)驗(yàn)時(shí)指出人工林灌草地侵蝕性雨量標(biāo)準(zhǔn)為19.9—27.8 mm??梢姡谌狈搅髟囼?yàn)小區(qū)和水土保持監(jiān)測數(shù)據(jù)的情況下,選取日降雨量≥20 mm/d 可能更符合廣西南流江流域?qū)嶋H特征[29]。本文通過計(jì)算分析河流輸沙量與流域降雨侵蝕力的相關(guān)性,進(jìn)一步支持了這一判斷。

    4 結(jié)論

    1961—2006年南流江流域平均降雨侵蝕力為13 935.5 MJ·mm·hm-2·h-1·a-1,年降雨侵蝕力變化周期約為14 a,最大值出現(xiàn)在1981年,為22 856.0 MJ·mm·hm-2·h-1·a-1,最小值出現(xiàn)在1962年,為6 104.0 MJ·mm·hm-2·h-1·a-1,各年際和年代間變化差異不大。在各時(shí)間段尺度上,2001—2006年的年均降雨侵蝕力最大(16 381.3 MJ·mm·hm-2·h-1·a-1),1971—1980年的降雨侵蝕力平均值最小(13 132.3 MJ·mm·hm-2·h-1·a-1)。降雨侵蝕力在各季節(jié)上也略有差異,表現(xiàn)為夏季>春季>秋季>冬季。空間分布格局上,降雨侵蝕力與年降雨量空間分布變化趨勢相似,呈現(xiàn)從下游西南沿海向東北內(nèi)陸遞減,降雨侵蝕力最大的區(qū)域是合浦縣。同時(shí),通過討論發(fā)現(xiàn),選取日降雨量≥20 mm/d作為南流江流域侵蝕性降雨標(biāo)準(zhǔn)是合理的。

    猜你喜歡
    侵蝕性輸沙量雨量
    青海省1961—2018 年侵蝕性降水時(shí)空變化特征
    寧夏紅柳溝流域水沙變化及產(chǎn)沙分析
    基于能量的坡面侵蝕性徑流及其水沙傳遞關(guān)系
    2011—2017年蘇北沿海侵蝕性降雨特征研究
    江蘇水利(2020年5期)2020-06-14 09:08:44
    基于小波去噪的稱重雨量數(shù)據(jù)分析
    米非司酮配伍化療治療侵蝕性葡萄胎的療效觀察
    遼河干流遼中站河道沖淤變化趨勢研究
    氣候變化和人類活動(dòng)對祖厲河輸沙量變化的影響分析
    頻率曲線在荊江三口輸沙量計(jì)算中的應(yīng)用
    SL—1 型雨量傳感器故障分析排除和維護(hù)
    西藏科技(2016年5期)2016-09-26 12:16:40
    亚洲五月天丁香| 免费观看在线日韩| 成人综合一区亚洲| 久久草成人影院| 乱系列少妇在线播放| 俺也久久电影网| 一个人免费在线观看电影| 免费不卡的大黄色大毛片视频在线观看 | 国语自产精品视频在线第100页| 淫妇啪啪啪对白视频| 十八禁网站免费在线| 午夜爱爱视频在线播放| 亚洲欧美清纯卡通| 色播亚洲综合网| 中国国产av一级| 激情 狠狠 欧美| 一夜夜www| 免费av观看视频| 成人鲁丝片一二三区免费| 日本爱情动作片www.在线观看 | 少妇熟女aⅴ在线视频| 丝袜喷水一区| aaaaa片日本免费| 国产不卡一卡二| 男人舔奶头视频| 人人妻人人澡欧美一区二区| 久久国产乱子免费精品| 18禁黄网站禁片免费观看直播| 3wmmmm亚洲av在线观看| 欧美一区二区精品小视频在线| 日韩高清综合在线| 国产黄色小视频在线观看| 有码 亚洲区| 天堂动漫精品| 色综合亚洲欧美另类图片| 亚洲成人精品中文字幕电影| 日本成人三级电影网站| 亚洲中文字幕日韩| 在线观看66精品国产| 少妇高潮的动态图| 国产伦一二天堂av在线观看| 搞女人的毛片| 国产成人一区二区在线| 伊人久久精品亚洲午夜| 熟妇人妻久久中文字幕3abv| 黄色一级大片看看| 免费无遮挡裸体视频| 久久鲁丝午夜福利片| 亚洲,欧美,日韩| 成人高潮视频无遮挡免费网站| 老熟妇乱子伦视频在线观看| 欧美不卡视频在线免费观看| 亚洲精品一卡2卡三卡4卡5卡| 一级黄色大片毛片| 婷婷亚洲欧美| 免费av不卡在线播放| 99久久九九国产精品国产免费| 老熟妇乱子伦视频在线观看| 在线观看免费视频日本深夜| 中文亚洲av片在线观看爽| 亚洲中文字幕一区二区三区有码在线看| 插逼视频在线观看| 久久精品人妻少妇| 男女啪啪激烈高潮av片| 午夜福利18| 日日啪夜夜撸| 欧美日韩精品成人综合77777| av在线观看视频网站免费| 村上凉子中文字幕在线| 男女下面进入的视频免费午夜| www.色视频.com| 久久中文看片网| 丰满人妻一区二区三区视频av| 欧美精品国产亚洲| 久久久久久国产a免费观看| 丰满乱子伦码专区| 最近视频中文字幕2019在线8| 搡老妇女老女人老熟妇| 少妇猛男粗大的猛烈进出视频 | 天堂网av新在线| 精品久久久久久久人妻蜜臀av| 亚洲,欧美,日韩| 欧美人与善性xxx| 人人妻人人澡人人爽人人夜夜 | 国产aⅴ精品一区二区三区波| av在线老鸭窝| 中文字幕av在线有码专区| 精品少妇黑人巨大在线播放 | 亚洲精品影视一区二区三区av| 最后的刺客免费高清国语| 熟女人妻精品中文字幕| 精品人妻熟女av久视频| 夜夜夜夜夜久久久久| 一卡2卡三卡四卡精品乱码亚洲| 中文字幕人妻熟人妻熟丝袜美| 成人一区二区视频在线观看| 欧美高清成人免费视频www| 日韩,欧美,国产一区二区三区 | 丝袜美腿在线中文| 亚洲欧美清纯卡通| 久久99热6这里只有精品| 国内久久婷婷六月综合欲色啪| 老师上课跳d突然被开到最大视频| 联通29元200g的流量卡| 国产精品亚洲美女久久久| 一个人看的www免费观看视频| 99国产极品粉嫩在线观看| 亚洲国产精品合色在线| 联通29元200g的流量卡| 免费一级毛片在线播放高清视频| 赤兔流量卡办理| 人妻制服诱惑在线中文字幕| 国产精品99久久久久久久久| 99久久无色码亚洲精品果冻| 日日摸夜夜添夜夜爱| 亚洲av熟女| av在线天堂中文字幕| 最好的美女福利视频网| 日产精品乱码卡一卡2卡三| 久久精品国产亚洲av香蕉五月| 久久精品综合一区二区三区| 狠狠狠狠99中文字幕| 久久久国产成人免费| 晚上一个人看的免费电影| 国产在线精品亚洲第一网站| 亚洲丝袜综合中文字幕| 免费看光身美女| 啦啦啦韩国在线观看视频| 搞女人的毛片| 麻豆一二三区av精品| 国产伦在线观看视频一区| 欧美三级亚洲精品| 欧美日本亚洲视频在线播放| 91麻豆精品激情在线观看国产| 18禁在线播放成人免费| 久久天躁狠狠躁夜夜2o2o| a级毛色黄片| 国产高清三级在线| 免费看光身美女| 免费在线观看影片大全网站| 亚洲乱码一区二区免费版| 日韩制服骚丝袜av| 成人特级av手机在线观看| 18禁在线无遮挡免费观看视频 | 午夜福利在线观看免费完整高清在 | av黄色大香蕉| 国产亚洲精品av在线| 日本成人三级电影网站| av中文乱码字幕在线| 亚洲综合色惰| 少妇的逼好多水| 美女内射精品一级片tv| 午夜激情福利司机影院| 亚洲不卡免费看| 亚洲欧美日韩东京热| 免费看a级黄色片| 国产精品永久免费网站| 亚洲国产精品国产精品| 18+在线观看网站| 国产高清三级在线| a级毛色黄片| 午夜视频国产福利| 三级经典国产精品| 国产高清激情床上av| 狠狠狠狠99中文字幕| 日韩高清综合在线| 观看免费一级毛片| 亚洲国产高清在线一区二区三| av免费在线看不卡| 国产精品久久久久久av不卡| 久久热精品热| 久久精品国产亚洲av天美| 美女被艹到高潮喷水动态| 看片在线看免费视频| 亚洲av第一区精品v没综合| 晚上一个人看的免费电影| 蜜桃久久精品国产亚洲av| 国产大屁股一区二区在线视频| 内射极品少妇av片p| 亚洲精品日韩在线中文字幕 | 国产探花极品一区二区| 国产成人精品久久久久久| av天堂在线播放| 亚洲自偷自拍三级| 日韩av在线大香蕉| 非洲黑人性xxxx精品又粗又长| 国产av一区在线观看免费| 在线播放无遮挡| 免费观看在线日韩| 在线观看一区二区三区| 亚洲无线在线观看| 性插视频无遮挡在线免费观看| 麻豆乱淫一区二区| 人妻丰满熟妇av一区二区三区| 一级毛片电影观看 | 亚洲av不卡在线观看| 国产伦精品一区二区三区四那| 你懂的网址亚洲精品在线观看 | 国产av一区在线观看免费| 亚洲欧美日韩东京热| 日产精品乱码卡一卡2卡三| 色综合色国产| 国产黄片美女视频| 国产综合懂色| 看免费成人av毛片| 日产精品乱码卡一卡2卡三| 亚洲国产欧美人成| videossex国产| 欧美另类亚洲清纯唯美| 日韩在线高清观看一区二区三区| 三级国产精品欧美在线观看| 丰满的人妻完整版| 嫩草影院新地址| 精品无人区乱码1区二区| 欧美一级a爱片免费观看看| 精品福利观看| 日韩精品有码人妻一区| 深爱激情五月婷婷| 在线免费十八禁| av黄色大香蕉| 国产毛片a区久久久久| 国产一区二区三区av在线 | 女同久久另类99精品国产91| 日本色播在线视频| 九九热线精品视视频播放| 久久九九热精品免费| 亚洲不卡免费看| 国产伦在线观看视频一区| 又粗又爽又猛毛片免费看| 亚洲熟妇中文字幕五十中出| 久久热精品热| 午夜影院日韩av| 国产一区二区在线av高清观看| 亚洲欧美清纯卡通| 免费黄网站久久成人精品| 人妻制服诱惑在线中文字幕| 蜜臀久久99精品久久宅男| or卡值多少钱| 一本一本综合久久| 欧美国产日韩亚洲一区| 国语自产精品视频在线第100页| 少妇人妻一区二区三区视频| 亚洲最大成人av| 中文字幕精品亚洲无线码一区| 国产在线男女| 国产爱豆传媒在线观看| 亚洲自偷自拍三级| 99久国产av精品国产电影| 日韩欧美国产在线观看| 精品不卡国产一区二区三区| 日韩中字成人| а√天堂www在线а√下载| 97超级碰碰碰精品色视频在线观看| 午夜视频国产福利| 两个人的视频大全免费| 高清日韩中文字幕在线| av卡一久久| 国产精品久久电影中文字幕| 久久鲁丝午夜福利片| 亚洲精品日韩av片在线观看| 久久精品综合一区二区三区| 国产高清有码在线观看视频| 蜜桃亚洲精品一区二区三区| 久久久久久久久久成人| 国产午夜精品久久久久久一区二区三区 | 亚洲欧美成人综合另类久久久 | 九九在线视频观看精品| 国产蜜桃级精品一区二区三区| 91久久精品国产一区二区成人| 国产亚洲欧美98| 97超碰精品成人国产| 久久久久久九九精品二区国产| 亚州av有码| 国产欧美日韩一区二区精品| 国产色婷婷99| 色哟哟哟哟哟哟| 国产欧美日韩精品亚洲av| 最近在线观看免费完整版| 天堂影院成人在线观看| 夜夜看夜夜爽夜夜摸| 99在线人妻在线中文字幕| 可以在线观看的亚洲视频| 国产在线男女| 日韩一本色道免费dvd| 嫩草影院精品99| 欧美3d第一页| 国产精品99久久久久久久久| 久久精品久久久久久噜噜老黄 | 高清午夜精品一区二区三区 | 国产精品一区二区三区四区久久| 一本久久中文字幕| 十八禁国产超污无遮挡网站| 色哟哟哟哟哟哟| 国产一区二区激情短视频| 人妻少妇偷人精品九色| 成人综合一区亚洲| 亚洲一区二区三区色噜噜| 女生性感内裤真人,穿戴方法视频| 人妻少妇偷人精品九色| 亚洲av中文字字幕乱码综合| 亚洲激情五月婷婷啪啪| 99九九线精品视频在线观看视频| 桃色一区二区三区在线观看| 国产探花极品一区二区| av国产免费在线观看| 少妇猛男粗大的猛烈进出视频 | 搡老妇女老女人老熟妇| 成熟少妇高潮喷水视频| 色综合色国产| 色综合亚洲欧美另类图片| 男人的好看免费观看在线视频| 久久久久久久亚洲中文字幕| 国产精品一区二区性色av| 麻豆国产97在线/欧美| 日韩 亚洲 欧美在线| 91在线观看av| 亚洲美女黄片视频| 亚洲av第一区精品v没综合| 欧美高清性xxxxhd video| 国产精品亚洲一级av第二区| 99久久中文字幕三级久久日本| 在现免费观看毛片| 免费一级毛片在线播放高清视频| 亚洲精品一区av在线观看| 69人妻影院| 亚洲图色成人| 午夜福利在线观看免费完整高清在 | 国产精品1区2区在线观看.| 国产一区亚洲一区在线观看| 亚洲欧美中文字幕日韩二区| av专区在线播放| 蜜臀久久99精品久久宅男| 卡戴珊不雅视频在线播放| 亚洲人成网站高清观看| 91久久精品电影网| 综合色av麻豆| 午夜精品国产一区二区电影 | 亚洲七黄色美女视频| 高清毛片免费看| 三级国产精品欧美在线观看| 亚洲国产欧美人成| 国产精品一二三区在线看| 特级一级黄色大片| 男女做爰动态图高潮gif福利片| 久久天躁狠狠躁夜夜2o2o| 青春草视频在线免费观看| 亚洲美女黄片视频| 亚洲国产精品成人综合色| 久久久久久久亚洲中文字幕| 日韩人妻高清精品专区| 联通29元200g的流量卡| 精品少妇黑人巨大在线播放 | 蜜桃久久精品国产亚洲av| 12—13女人毛片做爰片一| 国产精品福利在线免费观看| 男人舔奶头视频| 69人妻影院| 淫妇啪啪啪对白视频| 你懂的网址亚洲精品在线观看 | 波多野结衣高清作品| 亚洲国产精品成人久久小说 | 国产午夜精品论理片| 国产高清激情床上av| 国产亚洲精品久久久com| 日韩在线高清观看一区二区三区| 亚洲自拍偷在线| 一区二区三区免费毛片| 欧美潮喷喷水| 色在线成人网| 午夜老司机福利剧场| 国产欧美日韩精品一区二区| 亚洲成人精品中文字幕电影| 一区二区三区免费毛片| 国产精品野战在线观看| 97人妻精品一区二区三区麻豆| 卡戴珊不雅视频在线播放| 级片在线观看| 免费在线观看影片大全网站| 噜噜噜噜噜久久久久久91| 99热这里只有精品一区| 两个人的视频大全免费| 成人漫画全彩无遮挡| 亚洲不卡免费看| 人妻夜夜爽99麻豆av| 日本a在线网址| 中文在线观看免费www的网站| 91久久精品电影网| 日本爱情动作片www.在线观看 | 亚洲成人精品中文字幕电影| 97热精品久久久久久| 99久久精品一区二区三区| 久久精品国产自在天天线| 男人和女人高潮做爰伦理| 日日摸夜夜添夜夜添av毛片| 久久精品久久久久久噜噜老黄 | 男女啪啪激烈高潮av片| 日韩av在线大香蕉| 99精品在免费线老司机午夜| 美女大奶头视频| 欧美成人一区二区免费高清观看| 大又大粗又爽又黄少妇毛片口| 搡老岳熟女国产| 男人和女人高潮做爰伦理| 天天一区二区日本电影三级| 亚洲精品色激情综合| 国产69精品久久久久777片| 国产精品亚洲一级av第二区| a级毛片a级免费在线| 如何舔出高潮| av在线蜜桃| 欧美激情在线99| 99热精品在线国产| 国产精品无大码| 日韩亚洲欧美综合| 国产探花极品一区二区| 天堂√8在线中文| 九九爱精品视频在线观看| 寂寞人妻少妇视频99o| 超碰av人人做人人爽久久| 国产大屁股一区二区在线视频| 1000部很黄的大片| 欧美zozozo另类| 久久久久久久久中文| 99热网站在线观看| aaaaa片日本免费| 三级经典国产精品| 亚洲最大成人av| 精品久久久噜噜| 久久这里只有精品中国| 秋霞在线观看毛片| 成人亚洲欧美一区二区av| 级片在线观看| 欧美日韩在线观看h| 成人性生交大片免费视频hd| 欧美不卡视频在线免费观看| 少妇的逼好多水| 亚洲欧美精品综合久久99| 日韩欧美 国产精品| 午夜福利在线观看吧| 日韩一区二区视频免费看| 午夜福利18| 欧美色视频一区免费| 久久久久国产网址| 国产黄a三级三级三级人| 亚洲图色成人| 亚洲自拍偷在线| 简卡轻食公司| 欧美最新免费一区二区三区| 黄色视频,在线免费观看| 久久人妻av系列| 九九在线视频观看精品| 波多野结衣高清作品| 国产乱人偷精品视频| 最好的美女福利视频网| 别揉我奶头~嗯~啊~动态视频| 伦理电影大哥的女人| 97在线视频观看| 又黄又爽又刺激的免费视频.| 看非洲黑人一级黄片| 国产淫片久久久久久久久| 大香蕉久久网| 久久久久久国产a免费观看| 亚洲不卡免费看| www日本黄色视频网| 日日干狠狠操夜夜爽| 日韩亚洲欧美综合| 午夜爱爱视频在线播放| 国产高清不卡午夜福利| 国产精品久久视频播放| av福利片在线观看| 久久6这里有精品| 97超碰精品成人国产| 九色成人免费人妻av| 精品一区二区三区av网在线观看| 国产不卡一卡二| 成年免费大片在线观看| 永久网站在线| 国产视频一区二区在线看| 国产一区二区在线av高清观看| 亚洲美女搞黄在线观看 | 欧美在线一区亚洲| 又爽又黄a免费视频| 亚洲精品亚洲一区二区| 18禁裸乳无遮挡免费网站照片| 亚洲欧美成人精品一区二区| 天堂√8在线中文| 久久久久精品国产欧美久久久| 男女下面进入的视频免费午夜| a级毛色黄片| 不卡一级毛片| 最近手机中文字幕大全| 三级男女做爰猛烈吃奶摸视频| 男女之事视频高清在线观看| 伦精品一区二区三区| 淫秽高清视频在线观看| 国产精品电影一区二区三区| 国产高清视频在线观看网站| 国产av一区在线观看免费| 日韩欧美精品v在线| 亚洲成a人片在线一区二区| 成人毛片a级毛片在线播放| 久久九九热精品免费| 一级黄色大片毛片| 免费观看人在逋| 日日摸夜夜添夜夜添av毛片| 国产精品日韩av在线免费观看| 黄色欧美视频在线观看| 久久这里只有精品中国| 老司机福利观看| 91av网一区二区| 日本黄色片子视频| 国产av在哪里看| 久久精品国产亚洲网站| 91午夜精品亚洲一区二区三区| 精品欧美国产一区二区三| 三级男女做爰猛烈吃奶摸视频| 天天躁日日操中文字幕| 一卡2卡三卡四卡精品乱码亚洲| 欧美+日韩+精品| 男女之事视频高清在线观看| 亚洲欧美日韩东京热| 天天躁日日操中文字幕| 成年女人看的毛片在线观看| 麻豆成人午夜福利视频| 久久久午夜欧美精品| 男人狂女人下面高潮的视频| 久久精品国产自在天天线| 免费黄网站久久成人精品| 国产精品久久久久久精品电影| av黄色大香蕉| 日韩欧美三级三区| 久久久久久久亚洲中文字幕| 五月玫瑰六月丁香| av黄色大香蕉| 久久精品国产清高在天天线| 丝袜美腿在线中文| 99热这里只有精品一区| www日本黄色视频网| 欧美中文日本在线观看视频| 欧美日本亚洲视频在线播放| 久久99热6这里只有精品| 欧美日本亚洲视频在线播放| 美女高潮的动态| 最好的美女福利视频网| 国产精品一二三区在线看| 午夜福利18| 欧美3d第一页| 国产一区二区在线观看日韩| 成人亚洲精品av一区二区| 欧美bdsm另类| 亚洲国产精品久久男人天堂| 久久精品久久久久久噜噜老黄 | 我的老师免费观看完整版| 毛片女人毛片| 两个人的视频大全免费| 亚洲成人久久爱视频| 久久精品人妻少妇| 少妇熟女欧美另类| 亚洲自拍偷在线| 国产精品久久电影中文字幕| 九九久久精品国产亚洲av麻豆| 亚洲第一电影网av| 亚州av有码| 可以在线观看毛片的网站| a级毛片免费高清观看在线播放| 欧美日本视频| 日本免费a在线| 内射极品少妇av片p| 日日摸夜夜添夜夜爱| 亚洲专区国产一区二区| 亚洲va在线va天堂va国产| 我的老师免费观看完整版| 亚洲国产精品成人久久小说 | 日韩一本色道免费dvd| 久久99热6这里只有精品| 精华霜和精华液先用哪个| 亚洲久久久久久中文字幕| 丰满的人妻完整版| 观看免费一级毛片| 国产亚洲91精品色在线| 欧美最黄视频在线播放免费| or卡值多少钱| 亚洲欧美日韩东京热| 亚洲婷婷狠狠爱综合网| 亚洲国产高清在线一区二区三| 99久国产av精品国产电影| 国产爱豆传媒在线观看| 日本黄色片子视频| 亚洲国产精品sss在线观看| 如何舔出高潮| 国产一区二区在线观看日韩| 亚洲美女搞黄在线观看 | 国产91av在线免费观看| 久久久久久伊人网av| 联通29元200g的流量卡| 久久久精品欧美日韩精品| 伊人久久精品亚洲午夜| 中国美女看黄片| 国产精品一区二区性色av| 99热全是精品| 亚洲精品成人久久久久久| 日本在线视频免费播放| 亚洲18禁久久av| 欧美高清成人免费视频www| 精品免费久久久久久久清纯| 国产白丝娇喘喷水9色精品| 一区二区三区免费毛片| 身体一侧抽搐| 国产日本99.免费观看| 在线观看一区二区三区| av在线老鸭窝| 成年女人永久免费观看视频| 两性午夜刺激爽爽歪歪视频在线观看| 国产成人a∨麻豆精品| 免费观看在线日韩|