999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

基于多源數據的寧夏賀蘭縣現代化生態灌區健康評價

2021-03-17 07:20:12楊振峰田軍倉馮克鵬
農業工程學報 2021年24期
關鍵詞:現代化生態水平

楊振峰,田軍倉,2,3※,馮克鵬,2,3,龔 雪

基于多源數據的寧夏賀蘭縣現代化生態灌區健康評價

楊振峰1,田軍倉1,2,3※,馮克鵬1,2,3,龔 雪1

(1. 寧夏大學土木與水利工程學院,銀川 750021;2. 寧夏節水灌溉與水資源調控工程技術研究中心,銀川 750021;3. 旱區現代農業水資源高效利用教育部工程研究中心,銀川 750021)

現代化生態灌區是灌區發展方向,因此,灌區健康水平評價具有重要意義。該研究根據現代化生態灌區的概念、內涵要求,應用知識圖譜可視化工具VOSviewer對國內外灌區綜合評價研究的指標體系進行系統梳理,篩選共現頻率高的關鍵指標,在結合實際、廣泛調研的基礎上,構建了包括生態環境、輸配水現代化、田間工程現代化、生產及管理現代化、生產效益5類一級指標,共計18個二級指標的寧夏賀蘭縣現代化生態灌區綜合評價指標體系,并以鄉鎮為尺度劃分評價單元,結合遙感數據、調研統計、田間試驗數據實現指標定量化,利用層次-模糊綜合評價方法,對寧夏賀蘭縣現代化生態灌區健康水平進行綜合評價。結果表明:1)對于權重較大的二級指標信息化管理水平(權重為0.100)、主要排水溝水質(權重為0.091)、灌溉保證能力(權重為0.090),各評價單元健康水平較好且單元之間差異性小;2)各評價單元綜合評價結果顯示,南梁臺子及京星農牧場健康水平等級值最高(0.81),最低為洪廣鎮(0.69),所有評價單元均處于二至三級健康水平之間;3)賀蘭縣灌區整體處于一級健康水平的測度值(0.325)最大,表明灌區整體處于一級健康水平,但一至四級健康水平的測度值差別較小,存在健康水平惡化的風險較大,可從發展節水灌溉措施、提高規模化生產程度、完善管理制度、增大水域面積等方面優化灌區結構。該研究可為寧夏乃至西北現代化生態灌區建設提供參考。

遙感;生態;灌區;指標體系;健康評價;寧夏賀蘭縣

0 引 言

灌區是糧食生產的基本保障,是現代化農業發展的基礎,對社會經濟發展起重要作用[1]。中國社會經濟的快速發展,為現代化生態灌區的建設提供了有利條件[2],建設現代化、生態型灌區成為灌區發展的主要方向。Rapport等[3]對生態系統健康的定義是:隨著時間推移,外界條件改變時,生態系統自身結構及功能依然能夠保持穩定,并具有可恢復的能力。生態灌區是生態系統健康概念的延伸,其目的是維持灌區系統各要素之間相對平衡,使整體發展具有可持續性。在生態灌區及其評價研究領域,茆智院士是國內最早提倡建立節水、防污生態型灌區的學者,他提出節水灌溉措施可以提高單位面積產量;田間節水+小型濕地+生態溝的工程形式能夠有效抑制灌區污染,防止灌區生態環境惡化[4]。除嚴峻的生態環境狀況之外,生產力水平及現代化程度低也是制約灌區發展的重要因素,建設現代化生態灌區是提升灌區生產能力的關鍵。國內有關現代化生態灌區的研究起步較晚,目前關于現代化生態灌區的嚴格定義尚不統一。綜合來說,實現“自然-社會-經濟”可持續和諧發展、達到水利-農業-農村協調發展、生產力水平高、水資源配置合理、生物多樣性高是現代化生態灌區的核心內容,也是現代化生態灌區健康的內涵要求[4-7]。依據現代化生態灌區內涵要求,針對灌區實際生產中突出的生態環境惡化、生產力水平低等問題,進行灌區健康評價有助于掌握灌區運行狀況,及時優化灌區建設方案,從而實現灌區可持續發展。國內外學者依據特定區域的實際問題做了一些研究。Balasubramani等[8]在GIS支持下,評估了安迪帕蒂流域(印度)土地和水資源的潛力、利用水平及土地對農業的適宜性,制定了優化后的管理措施。Oad等[9]針對美國西部人口用水、生態環境用水及農業用水之間的突出矛盾,利用決策支持系統為奧格蘭德保護區提供最優的供配水制度。Sun等[10]建立了灌區農業用水管理評價指標體系,包括技術指標、工程指標、管理指標、環境指標和經濟指標,通過改進的層次分析法確定各指標的權重,采用灰色關聯法和模糊綜合評價法評估了黃淮海流域汾河、石津和人民勝利渠灌區農業用水管理等級。

在綜合評價方法方面,模糊評價法[11-12]、灰色關聯法[13-15]、數據包絡法及神經網絡方法[16-17]等廣泛應用于灌區綜合評價研究。其中,模糊評價法是利用模糊數學對受到多種因素制約的事物或對象做出總體評價,適合各種非確定性問題的解決。在指標體系方面,文獻[18-20]探討了衡量灌區水利設施現代化程度的評價指標,也有學者提出了寧夏節水型灌區建設及運行狀況的評價標準[21-22]。通過層次分析(Analytic Hierarchy Process,AHP)確定指標權重的方法具有靈活、易于實現及多維客觀決策的優勢,使復雜問題層次化、系統化,被廣泛應用于綜合評價研究[10,23]。方延旭等[12]基于層次結合模糊的綜合評價模型對北京市北野廠灌區進行綜合評價,結果顯示該模型能較好地評價灌區生態系統健康狀況。

然而,灌區現代化程度及生態狀況處于動態變化中,及時更新灌區健康水平信息對于制定管理計劃、調整管理措施具有重要意義。以往有關灌區綜合評價的研究以統計資料為主,且以灌區整體為評價單元,忽略了灌區內各行政單元之間的差異性。遙感技術具有時效性高、空間特征明顯的特點,成為灌區管理及快速決策的首選工具。但利用遙感數據產品構建灌區現代化程度及生態環境狀況評價的具體指標有待深入。基于此,本文應用VOSviewer可視化工具分析灌區評價相關文獻常用的評價指標,結合現代化生態灌區概念、內涵及實際調研構建評價指標體系;通過遙感技術定量化指標體系,依據灌區建設規劃要求制定健康水平等級;劃分以鄉鎮為單位的評價單元,應用層次-模糊評價模型評價灌區現代化程度及生態健康水平。旨在探討利用遙感數據融合統計、田間試驗等數據定量化評價指標,在灌區綜合評價中的應用效果,以期為寧夏現代化生態灌區建設提供參考。

1 研究區概況及評價體系構建

1.1 研究區概況

賀蘭縣地處銀川平原河西灌區中部(東經105°57′~106°36′、北緯38°27′~38°52′(圖1),屬中溫帶大陸性氣候,干旱少雨,蒸發強烈,光照時間長,晝夜溫差大。東臨黃河,西倚賀蘭山。轄區總面積1 201.7 km2。地勢西高東低,自東向西3種地貌形態,即黃河沖積平原、山前洪積平原及賀蘭山山脈。瀕臨黃河,獨特的地理及氣候條件使賀蘭縣農業主要依靠灌溉系統,因此,賀蘭縣有著悠久的灌溉歷史。西干渠、唐徠渠、漢延渠、惠農區四大干渠自西向東流經賀蘭縣,總長為112.3 km,灌溉控制面積可達3 900 hm2。支斗渠縱橫交錯,干渠直開口支斗渠總計189條。賀蘭縣是寧夏回族自治區現代化生態灌區建設的示范區域,評估現代化灌區建設成效,將有助于灌區管理政策的制定。

1.2 評價指標體系構建

為明確灌區綜合評價研究的重點問題,本文利用VOSviewer文獻可視化工具進行文獻知識圖譜共現分析。在Web of Science核心數據庫以“Irrigation District + Evaluation”為主題進行文獻檢索,檢索時間是2020年12月20號,時間范圍是1958—2020年,共得到448條記錄。圖2a是檢索結果關鍵詞知識圖譜,從圖中可以看出,“Irrigation”“management”“quality”“impact”“groundwater quality”“salinity”“water”“suitability”“system”“agriculture”等關鍵詞占比較大,且彼此之間的關聯性較好,說明灌區評價關注灌溉、水質量、環境質量及系統性管理等問題。在CNKI數據庫以“灌區+評價”為主題進行檢索,共得到643條檢索記錄,對其進行可視化分析,如圖2b所示。從圖中可以看出,灌區評價大致可以分為灌區綜合評價和灌區污染評價兩大模塊,灌區綜合評價包含的范圍更廣泛一些,其中“水資源”“節水改造”“層次分析法”“指標體系”“灌溉”“水質”“地下水”等關鍵詞是重要組成部分,灌區污染模塊主要包括“重金屬污染”“污染灌區”“土壤重金屬”“健康風險”等關鍵詞。通過文獻關鍵詞知識圖譜分析,本研究從灌區管理、灌溉、鹽漬化水平、節水、水質等方面構建指標體系。

依據VOSviewer文獻關鍵詞可視化分析結果,結合現代化生態灌區概念[4-7,24]、灌區評價研究綜述[25]及寧夏現代化生態灌區指標考核標準[26-27],本研究針對灌區健康評價(目標層)所取指標體系如表1所示。

表1 現代化生態灌區健康評價指標體系及表征意義

2 數據處理與方法

2.1 指標數據來源及處理

2.1.1 生態環境

1)土壤肥力

土壤肥力包括有機質、全氮、有效磷、速效鉀含量4項內容,構建評價準則如表2所示,采用包長征等[28]關于賀蘭縣土壤質地調查的樣點平均值數據。該數據采用隨機等距、分層的采樣方法,共布置樣點4 587個,數據具有良好的可靠性。

2)鹽漬化水平

研究證明,基于遙感數據進行土壤鹽漬化水平評價具有可靠性[29-31],本研究參考陳實等[32]基于鹽度指數(Salinity Index,SI)、歸一化植被指數(Normalized Difference Vegetation Index,NDVI)構建土壤鹽分狀態模型評價寧夏賀蘭縣土壤鹽漬化水平,該方法是基于植被生長狀況與土壤含鹽量具有高度的相關性,NDVI與土壤鹽分呈現負相關趨勢,鹽度指數能較好地區分土壤鹽漬化程度。計算如下:

式中1、2分別為藍光與紅光波段反射率。

式中3、4分別代表近紅光和紅光波段反射率。

表2 土壤肥力等級

通過計算NDVI和SI指數,構成二維特征空間,從而建立土壤鹽漬化的遙感監測模型(Salinization Remote Sensing Index,SRSI),計算公式如下:

選用Landsat8遙感影像數據(中國遙感數據網http://ids.ceode.ac.cn),空間分辨率30 m,數據獲取時間為2020年7月19日、2020年8月13日,此時間避免了灌溉高峰期對土壤鹽漬化評價帶來的誤差。

3)水域面積確定

利用遙感影像提取水體面積、幾何形態、水質狀況和水體生態環境等信息取得廣泛應用[33]。水體指數法通過光譜知識構建分類模型進行水體信息的提取,具有簡單可行的特點,國內外學者構建了大量的水體指數,可參考文獻[34]。依據各水體指數特點結合賀蘭縣實際情況,本文選擇歸一化水體指數(Normalized Difference Water Index,NDWI)來提取水體信息。采用Landsat8遙感影像數據(中國遙感數據網http://ids.ceode.ac.cn),數據獲取時間為2020年8月13日,此時間水稻覆蓋度較大,從而降低了水田對提取結果的影響,經過輻射校正、大氣校正之后,計算NDWI指數提取水體信息。水體信息主要為賀蘭縣小型湖泊面積,利用無人機影像目視解譯面積對計算結果進行校正。無人機影像來源于2018年課題組對賀蘭縣小型湖泊的航測數據。

4)主要排水溝水質

主要排水溝水質根據銀川市生態環境監測中心對轄區內黃河、湖泊濕地、排水溝水質監測的年平均數據。主要排水溝水質是灌區生態環境污染程度的重要組成部分,影響排水溝水質的原因除了不同種植作物的施肥量,田間用藥等不同之外,還與排水溝結構有關,對于跨度在多個評價單元之間的排水溝,排水溝之間的疊加作用也無法消除。因此,為了量化同一排水溝在各評價單元水質情況,利用進入和出評價單元時的水質情況變化顯示評價單元排水溝受污染及去污染能力(表3)。

表3 排水溝水質健康水平等級準則

5)植被干旱指數確定

干旱給農業生產帶來嚴重威脅,建立快速響應的干旱評價是現代化生態灌區健康評價的重要組成部分。國內外學者對于干旱指數已開展廣泛的研究,大致可以分為4類:氣象干旱指數、水文干旱指數、農業干旱指數以及基于遙感技術的植被指數[35]。本文選擇基于遙感技術的植被供水指數(Vegetation Supply Water Index,VSWI)作為干旱程度評價的指數,VSWI是由NDVI與地表溫度(Land Surface Temperature,LST)的比值計算得到。數據來源于NASA網站(http://reverb.echo.nasa.gov/reverb/)2020年5—8月份的MODIS13A3植被指數數據集,該數據集包括16 d合成的NDVI數據,LST數據采用8 d合成MODIS11A2數據集,空間分辨率為500 m。VSWI值越大代表土壤含水率越高,植物受水分脅迫程度越小。

2.1.2 輸配水現代化

1)骨干渠砌護率、排水溝整治率及量水設備配比數據均來自于2019年1月—2021年3月實際調研結果。

2)灌溉保證能力計算過程參考文獻[36]。遙感數據利用2020年5—8月的MODIS16A2數據集,包括蒸散發、潛在蒸散發、潛熱通量、潛在熱通量數據,空間分辨率為1 km,經過裁剪之后得到研究區的月合成數據集。5 —8月為灌溉高峰期,也是蒸散發最大的時間段。同時,利用同一時期研究區氣象數據(銀川氣象站(站點號53614))計算參考蒸散量(ET0)。為了量化主要種植結構信息,實際采樣時間是2020年9月10—12日,共采集樣本85個,均勻布滿整個研究區。采集樣本時,單個樣方大小為80 m×80 m,用RTK記錄樣方點的位置信息,坐標系統采用WGS84坐標系。

3)灌溉水利用系數計算采用“灌區首尾”比較的綜合方法,計算公式如下:

2.1.3 田間工程現代化

高效節水灌溉面積占比(高效節水灌溉面積/總面積)、高標準農田整治率(整治農田面積/總面積)數據來自于2019年1月—2021年3月實際調研及統計數據。

2.1.4 生產及管理現代化

機械化使用率、規模化生產率、水利人才配比(水利專業人數/667 hm2)、管理制度完善度、信息化管理水平數據均來自于2019年1月—2021年3月實際調研及統計數據,其中,管理制度完善度及信息化管理水平由專家評審。

2.1.5 生產效益

水分生產效率(Water Production Efficiency,WPE)采用課題組2017年1月—2020年12月在賀蘭縣針對主要作物(糧食作物:水稻、玉米、小麥;蔬菜、葡萄及枸杞)進行的大田試驗觀測數據,其計算如下:

WPE=/(5)

式中WPE為水分生產效率,kg/m3;為產量,kg;為灌溉定額,m3。

農民人均收入數據來自于縣統計局。

2.2 評價單元劃分

綜合分析調研、訪談、監測數據資料發現,土地質量、農田水利設施、農作物種植種類、特色農業產業等在鄉鎮尺度上表現出差異,土地流轉的規模化種植大多集中于鄉鎮及村級尺度,高標準農田建設及農田水利建設項目單元在鄉鎮尺度上。結合現代化生態灌區“自然-社會-經濟”的綜合效應,本研究評價單元的最優選取為在鄉鎮級尺度,共建立習崗鎮、金貴鎮、立崗鎮、洪廣鎮、常信鄉、京星農牧場、暖泉農場及南梁臺子8個評價單元。

2.3 指標權重計算

現代化生態灌區建設是按照人的意愿對灌區進行升級改造,灌區健康評價即對灌區各項指標達標情況進行分析。因此,選擇AHP法計算指標權重[37-38]。通過專家對各指標的評價打分,計算得到指標權重如表4所示。

2.4 綜合評價模型構建

表4 指標權重

3 結果與分析

3.1 單元評價結果

3.1.1 生態環境評價

土壤肥力方面:除地勢較高的洪廣鎮、暖泉農場及南梁臺子之外,其余地區有機質(D111)均屬于I級肥力水平,南梁臺子土壤有機質含量水平最低;各評價單元全氮含量(D112)均較低,洪廣鎮、暖泉農場及南梁臺子屬于IV級肥力水平,其余為III級肥力水平;各評價單元有效磷含量(D113)具有較高水平,除南梁臺子屬于III級肥力水平之外,其余評價單元均為II級肥力水平;除習崗鎮、金貴鎮及立崗鎮速效鉀含量(D114)屬于II級肥力水平之外,其余評價單元均低于II級肥力水平。依據表2土壤肥力C11構成要素,得出各評價單元土壤肥力情況,除南梁臺子土壤肥力水平低于其他評價單元之外,其余各評價單元土壤肥力較為一致,但均低于二級健康水平,因此,改善土壤肥力狀況、防止營養元素流失依然是各評價單元需要關注的問題,尤其是對于南梁臺子,較低的土壤肥力制約生產力水平。

利用LandSat 8遙感影像,計算鹽漬化水平模型SRSI,結果如圖3所示。陳實等[32]對北疆鹽漬化的研究中,認為SRSI在0.40~0.05之間為中度鹽漬,在>0.50~0.70之間為重度鹽漬,>0.70為鹽土。從圖3可以看出,洪廣鎮SRSI普遍較高,且集中在賀蘭山前的洪水沖積區,這與實際調查相符,證明構建的SRSI鹽度模型能夠反映灌區內鹽漬化情況。同時,圖3也顯示習崗鎮SRSI較大,多集中在0.70~1.00之間,習崗鎮位于賀蘭縣城區,轄區內多為建筑物,表明鹽度指數SRSI在裸土與建筑物的區分上存在一定的交集。立崗鎮SRSI多集中在0.35~0.55之間,常信鄉SRSI多集中在0.45~0.70之間。總體來說,各評價單元SRSI整體較高,鹽漬化狀況均低于三級健康水平,有效控制土壤鹽漬化程度是所有評價單元都面臨的生態環境問題。

通過計算水體指數NDWI,削弱植被、建筑等信息的影響,增強水體信息,進而提取水域面積(C13)。從圖 4提取的水體信息可以看出,常信鄉包含最多的水體信息。立崗鎮、洪廣鎮、南梁臺子水體信息較少,這主要與地勢高低有關。總體來說,各評價單元水域面積占比均較小,水域面積均低于三級健康水平。現代化生態灌區提倡“山、水、田、林、路、湖”交錯的布局形式,水系與湖泊是生物多樣性的保障,這也是水域面積指標權重較大的重要依據,但水域面積受限于地理等因素。近年來,各評價單元通過“稻漁空間”生態種植、農村宜居環境改善、修建人工湖泊等多種方式增加水域面積,改善生態環境。

水質是衡量灌區污染的重要手段,這使得排水溝水質指標權重高于其他指標權重。各評價單元主要排水溝水質(C14)健康水平如表5。京星農牧場主要排水溝水質為一級健康水平,這與京星農牧場距離黃河較近,排水溝易于排泄,且排水溝在評價單元之間的累積疊加程度較小有關。習崗鎮、金貴鎮、暖泉農場及南梁臺子為二級健康水平。暖泉農場和南梁臺子由于地勢較高,使得排水溝排水較易,另外,從種植結構上看,玉米為主要的種植作物,而種植玉米的田間管理帶來的污染較小(如施肥、施藥等)。常信鄉和洪廣鎮主要排水溝水質為四級健康水平,常信鄉和洪廣鎮是面積最大的評價單元,也是排水溝較密集的評價單元,從所處的地理位置及排水溝結構特點上來看,評價單元之間的疊加作用比其他單元明顯,種植面積大,單元內排水溝長度較長,滯排明顯使得水質惡化加劇,水質情況較其他評價單元差。

表5 各評價單元指標健康水平等級

注:構造健康等級評分向量,一級健康為1,二級健康為0.8,三級健康為0.6,四級健康為0.4;介于兩級之間均按高一級計算。

Note: The health level scale vectors are constructed for levels 1, 2, 3 and 4, corresponding to 1, 0.8, 0.6 and 0.4 respectively, and any level between any two health levels is calculated as the higher level.

圖5為VSWI指數(C15)空間分布圖,從圖中可以看出,洪廣鎮、習崗鎮VSWI值較小,表明所受的水分脅迫較大,而這兩個評價單元作物覆蓋面積小,因此,干旱對于洪廣鎮及習崗鎮農業生產的影響相對較小。計算各評價單元VSWI指數的平均值可知,京星農牧場、常信鄉及南梁臺子3個評價單元的VSWI值最高,表明干旱程度低,京星農牧場緊鄰黃河,地勢較低,補水容易,干旱程度最低;常信鄉水稻種植面積大,而水稻灌溉以大水漫灌為主,使得田間蓄水量大,干旱程度較低。總體來說,各評價單元干旱程度較小,表明灌區灌溉用水保證程度高。

3.1.2 輸配水現代化評價

所有評價單元干、支渠砌護率為100%,而在主要排水溝整治水平上表現出差異。因此,指標C21主要受排水溝整治情況影響。常信鄉主要排水溝整治情況較好,處于一級健康水平,其余評價單元存在部分排水溝未整治的情況,處于三、四級健康水平(如表5)。各評價單元在灌溉高峰期(5—8月)的灌溉保證能力空間分布如圖6所示。

從圖中可以看出,除洪廣鎮,其余評價單元有效灌溉量等于或大于灌溉需水量,表明各評價單元灌溉保證能力較好,一定程度上也反映了賀蘭縣具有節水灌溉的潛力和空間。灌溉水利用系數(C23)主要受各評價單元節水灌溉措施(如噴灌、滴灌)種植面積的影響,且與種植結構密切相關。各評價單元種植結構統計及田間試驗結果顯示,洪廣鎮、暖泉農場及南梁臺子灌溉水利用系數處于一級健康水平,其余評價單元為三級健康水平。各評價單元內各支渠口、用戶取水口均配有量水設備,各評價單元渠道輸水計量設施現代化程度高,各評價單元量水設備配比(C24)差異性主要體現在使用地下水灌溉的用戶,灌溉用水無法有效計量,這也反映了各評價單元應加強地下水的開采管理。

3.1.3 田間工程現代化水平評價

各評價單元高效節水面積占比(C31)存在較大差異,其中,京星農牧場高效節水實施面積占比最大,處于一級健康水平,洪廣鎮及暖泉農場處于二級健康水平,其余評價單元均處于四級健康水平。所有評價單元農田整治率(C32)均處于一級健康水平,表明農田整治建設情況較好。

3.1.4 生產及管理現代化評價

生產及管理現代化指標評價結果如表5所示。對于機械化使用率(C41)和管理制度完善度(C45),各評價單元均處于一級和二級健康水平,體現了各評價單元灌區現代化建設的同步性。規模化生產率(C42)主要受土地流轉面積的影響,因此,各評價單元具有較大的差異性,立崗鎮土地流轉面積比例最高,規模化生產水平處于一級健康水平。灌區信息化管理水平(C43)主要體現在農業用水的信息化管理水平,而其管理程度與種植結構及灌溉方式緊密相關,通過實際調研及統計計算,各評價單元信息化管理水平處于一至三級健康水平。水利人才配比(C44)主要受評價單元行政管理制度的影響,除習崗鎮、京星農牧場及南梁臺子處于一級健康水平外,其余評價單元均處于四級健康水平,這顯示了各評價單元水利人才配比情況受行政管理的政策及措施影響。

3.1.5 生產效益

田間試驗及調研結果顯示,洪廣鎮、暖泉農場、京星農牧場及南梁臺子4個評價單元的水分生產效率(C51)處于一級健康水平,其余評價單元處于二、三級健康水平,主要原因是受地勢影響,洪廣鎮、暖泉農場、京星農牧場及南梁臺子位于揚黃灌溉區域,水分生產效率相對較高,這也表明其余評價單元具有很大的節水潛力。居民人均收入(C52)各評價單元差異較小,整體處于二、三級健康水平。

3.2 綜合評價

3.2.1 各單元健康水平綜合評價

為了定量化各評價單元健康水平差異,以健康水平等級評分與具體指標權重的乘積累加和作為總體健康水平程度,計算結果如圖7所示。所有評價單元均高于三級健康水平,京星農牧場及南梁臺子總體健康水平最高(0.81),洪廣鎮綜合健康水平最低(0.69)其余評價單元在二、三級建康水平之間(0.71~0.76)。

3.2.2 賀蘭縣代現化生態灌區綜合評價結果

將指標權重集合與隸屬度矩陣代入模糊綜合評價模型中,得到健康等級的測度值。賀蘭縣現代化生態灌區屬于一級健康的測度值最大,為0.325,即灌區整體建設情況趨于一級健康水平,但各評價單元在自然等因素影響下,依然存在差異。而隸屬于二級健康的測度值(0.272)僅次于一級健康的測度值,且相差較小,說明灌區處于一二級健康水平的趨勢較明顯。一至四級健康水平的測度值依次減小(三級和四級分別為0.212和0.191),各級之間的測度值相差不大,說明灌區健康水平存在明顯的下降風險,依然要加強灌區建設的一些短板,維護建設成果。

4 討 論

關于評價指標的選取及評價指標值的確定是灌區綜合評價的難點,評價指標太多、指標值太高或太低,不利于指導實際生產,評價指標少又不能全面反映灌區現狀。因此,本研究在確定灌區評價指標時,采用了文獻分析熱點問題、實際調研、參考已有研究及寧夏現代化生態灌區建設的有關標準,優化指標體系的同時具有簡單易行的特點。盡管如此,本文構建的評價指標依然最優適用于評價灌區短期運行狀況,這主要是由于灌區自身條件決定適用于本身的評價指標體系及指標值,當灌區本身結構發生變化時,評價指標體系及指標值也應做出改變,因此,灌區綜合評價處于動態變化中。如史海濱等[41]研究表明,沈烏灌區節水改造工程對土壤鹽漬化時空分布產生影響。進行灌區健康水平評價及時發現這種動態變化,利用遙感技術是最有效的方法,本文探討了將遙感影像融入灌區綜合評價之中,這種方法不僅具有較高的準確性,還具有空間可視化優勢,有利于明確評價單元之間的差異。另外,本研究將評價指標值劃分為一至四級健康水平,一級健康水平值為灌區長遠期規劃設計值,二、三級健康水平值為短期規劃設計值,四級健康水平值表示現狀值低于短期規劃設計值,這一劃分標準參考《全國現代灌溉發展規劃》及有關寧夏引黃灌區對具體指標規劃年設計值。例如,杜軍等[26]認為到規劃年(2025年)灌溉水利用系數≥0.57,本文確定灌溉水利用系數一級健康水平為0.60。

應用構建的評價體系和評價方法,對寧夏賀蘭縣現代化生態灌區進行評價,結果顯示賀蘭縣整體處于一級健康水平,主要原因是權重較大的指標如生態環境B1(排水溝水質C14)、輸配水現代化B2(灌溉保證能力C22)、生產及管理現代化B4(灌區信息化管理水平C43)在灌區整體尺度上具有較好的健康水平,從而使得整體處于一級健康水平。但對于各評價單元,各項指標的健康水平存在差異(表5)。同時,也表明了將灌區視為一個整體給評價結果帶來的誤差,即權重較大指標的健康水平很大程度上決定灌區整體健康水平。基于這一考慮,本研究將灌區劃分為行政管理上相對獨立的若干評價單元,各評價單元健康水平的衡量相互獨立,使得評價結果能真實地反映灌區現狀,易于決策者制定優化措施。

5 結 論

本研究在總結前人關于灌區評價研究的基礎上,構建了現代化生態灌區健康水平評價體系,對寧夏賀蘭縣現代化生態灌區建設健康水平進行評價,得到如下結論:

1)在生態環境方面,受地理條件等因素影響,各評價單元表現出較大的差異,特別是對于權重較大的水域面積、主要排水溝水質指標,而干旱程度指標計算結果顯示賀蘭縣有一定的節水空間,因此,大力發展節水灌溉措施依然是提高水分生產效益的有效途徑;在田間工程現代化、生產及管理現代化、水利設施現代化方面,各評價單元表現出較小差異。

2)各評價單元綜合評價結果處于三級健康水平以上,健康等級值最高為京星農牧場及南梁臺子(0.81),最低為洪廣鎮(0.69),其余評價單元健康等級在0.71~0.76之間,這表明賀蘭縣各鄉鎮現代化生態灌區建設保持一致性,建設具有同步性。對于具體指標,各評價單元也存在差異,基于這些差異可以為各評價單元現代化及生態建設提供指導。

3)按照隸屬度最大原則,模糊綜合評價結果顯示賀蘭縣灌區整體處于一級健康水平(測度值為0.325),但依然存在健康水平惡化的風險,這主要是因為一至四級健康水平的測度值相差較小,應加強并保持灌區現代化及生態建設取得的成果。

[1] 翟浩輝. 加大灌區改造力度保障國家糧食安全[J]. 求是,2004(6):31-33.

[2] 呂純波. 關于現代生態灌區發展方向的思考[J]. 水利科學與寒區工程,2018,1(6):125-130.

Lyu Chunbo. Thoughts on the development direction of modern ecological irrigation area[J]. Hydro Science and Cold Zone Engineering, 2018, 1(6): 125-130. (in Chinese with English abstract)

[3] Rapport D J, Costanza R, McMichael A J. Assessing ecosystem health[J]. Trends in Ecology and Evolution, 1998, 13(10): 397-402.

[4] 茆智. 構建節水防污型生態灌區[J]. 中國水利,2009(19):28.

Mao Zhi. Construction of water-saving and pollution-proof ecological irrigation area[J]. China Water Resources, 2009(19): 28. (in Chinese with English abstract)

[5] 顧斌杰,王超,王沛芳. 生態型灌區理念及構建措施初探[J]. 中國農村水利水電,2005(12):7-9.

Gu Binjie, Wang Chao, Wang Peifang. Study on concept of ecological irrigation district and it’s construction measures[J]. China Rural Water and Hydropower, 2005(12): 7-9. (in Chinese with English abstract)

[6] 張展羽,孔莉莉,張國華. 我國生態灌區建設淺議[J]. 中國水利學會2006學術年會論文集[C]. 北京:中國水利水電出版社,2006.

[7] 楊培嶺,李云開,曾向輝,等. 生態灌區建設的理論技術及支撐技術體系研究[J]. 中國水利,2009(14):32-36.

Yang Peiling, Li Yunkai, Zeng Xianghui, et al. Ecological irrigation district: Studies on theory basis and supporting technical system[J]. China Water Resources, 2009(14): 32-36. (in Chinese with English abstract)

[8] Balasubramani K. Physical resources assessment in a semi-arid watershed: An integrated methodology for sustainable land use planning[J]. ISPRS Journal of Photogrammetry and Remote Sensing, 2018, 142:358-379.

[9] Oad R, Garcia L, Kinzli K D, et al. Decision support systems for efficient irrigation in the middle rio grande valley[J]. Journal of Irrigation and Drainage Engineering, 2009, 135(2): 177-185.

[10] Sun H, Wang S, Hao X. An improved analytic hierarchy process method for the evaluation of agricultural water management in irrigation districts of north China[J]. Agricultural Water Management, 2017, 179: 324-337.

[11] 趙振國. 生態友好灌區水資源聯合調度研究[D]. 新鄉:中國農業科學院農田灌溉研究所,2012.

Zhao Zhenguo. Research of Water Resources Management in Irrigation Aeras by Taking Into Account the Eco-Friendly[D]. Xinxiang: Farmland Irrigation Research Institute, 2012. (in Chinese with English abstract)

[12] 方延旭,楊培嶺,宋素蘭,等. 灌區生態系統健康二級模糊綜合評價模型及其應用[J]. 農業工程學報,2011,27(11):199-205.

Fang Yanxu, Yang Peiling, Song Sulan, et al. Secondary fuzzy comprehensive evaluation model for ecosystem health of irrigation district and its application[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2011, 27(11): 199-205. (in Chinese with English abstract)

[13] 王錦輝,費良軍,謝芳,等. 灰色關聯理論-TOPSIS法的大型灌區運行狀況綜合評價[J]. 排灌機械工程學報,2015,33(11):985-990.

Wang Jinhui, Fei Liangjun, Xie Fang, et al. Comprehensive evaluation of operational status of large-scale irrigation district based on grey relation theory and TOPSIS method[J]. Journal of Drainage and Irrigation Machinery Engineering, 2015, 33(11): 985-990. (in Chinese with English abstract)

[14] 門寶輝,梁川,趙燮京. 井灌區節水技術模式綜合評價模型[J]. 節水灌溉,2003(5):42-45.

Men Baohui, Liang Chuan, Zhao Xiejing. Synthetically evaluating model for the mode of water saving techniques in well irrigation district[J]. Water Saving Irrigation, 2003(5): 42-45. (in Chinese with English abstract)

[15] 費良軍,王錦輝,王光社,等. 基于灰色關聯理論-密切值法的大型灌區運行狀況綜合評價[J]. 干旱地區農業研究,2016,34(1):242-246.

Fei Liangjun, Wang Jinhui, Wang Guangshe, et al. Comprehensive evaluation on the operation condition of large irrigation district based on the grey relation theory-osculating value method[J]. Agricultural Research in the Arid Areas, 2016, 34(1): 242-246. (in Chinese with English abstract)

[16] 高峰,雷聲隆,龐鴻賓. 節水灌溉工程模糊神經網絡綜合評價模型研究[J]. 農業工程學報,2003,19(4):84-87.

Gao Feng, Lei Shenglong, Pang Hongbin. Model of fuzzy nervous network integrated assessment of water-saving irrigation projects[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2003, 19(4): 84-87. (in Chinese with English abstract)

[17] 王征,束銳,王天祎,等. 基于數據包絡分析的灌區農業生產綜合效率評價與分析[J]. 灌溉排水學報,2017,36(11):107-114.

Wang Zheng, Shu Rui, Wang Tianyi, et al. Efficiency of agriculture production in irrigation districts of china: assessed using data envelopment analysis method[J]. Journal of Irrigation and Drainage, 2017, 36(11): 107-114. (in Chinese with English abstract)

[18] 彭濤,高旺盛,隋鵬. 農田生態系統健康評價指標體系的探討[J]. 中國農業大學學報,2004,9(1):21-25.

Peng Tao, Gao Wangsheng, Sui Peng. Discussion on indicator system of farmland ecosystem health assessment[J]. Journal of China Agricultural University, 2004, 9(1): 21-25. (in Chinese with English abstract)

[19] 范習超,呂謀超,王曉森,等. 生態灌區的概念、屬性與指標體系研究[J]. 中國水利,2018(11):51-53.

Fan Xichao, Lyu Mouchao, Wang Xiaoseng, et al. Study on concept, attributes and indicator system of ecological irrigation districts[J]. China Water Resources, 2018(11): 51-53. (in Chinese with English abstract)

[20] 穆建新,呂振豫,許迪,等. 農田水利現代化評價指標體系及評價方法研究[J]. 中國農村水利水電,2016(8):33-40.

Mu Jianxin, Lyu Zhenyu, Xu Di, et al. Study on assessment indictor system and evaluation method of irrigation and drainage modernization[J]. China Rural Water and Hydropower, 2016(8): 33-40. (in Chinese with English abstract)

[21] 朱潔,劉學軍,陸立國. 寧夏節水型灌區考核評價標準研究[J]. 中國農村水利水電,2019(1):13-15,21.

Zhu Jie, Liu Xuejun, Lu Liguo. Research on evaluation standard of water-saving irrigation district in Ningxia[J]. China Rural Water and Hydropower, 2019(1): 13-15, 21. (in Chinese with English abstract)

[22] 劉瑾亮,劉學軍. 寧夏節水型灌區評價標準研究[J]. 水利技術監督,2018(4):23-26.

Liu Jingliang, Liu Xuejun. Study on evaluation standard of water-saving irrigation district in Ningxia[J]. Technical Supervision in Water Resources, 2018(4): 23-26. (in Chinese with English abstract)

[23] Karamouz M , Kerachian R , Zahraie B, et al. Monitoring and evaluation scheme using the multiple criteria decision making technique: Application to irrigation projects[J]. Journal of Irrigation and Drainage Engineering, 2002, 128(6): 341-350.

[24] 王超,王沛芳,侯俊,等. 生態節水型灌區建設的主要內容與關鍵技術[J]. 水資源保護,2015,31(6):1-7.

Wang Chao, Wang Peifang, Hou Jun, et al. Main contents and key technologies of constructing ecological water-saving irrigation district[J]. Water Resources Protection, 2015, 31(6): 1-7. (in Chinese with English abstract)

[25] 易銀蓮,崔寧博,張青雯,等. 灌區運行水平綜合評價研究進展[J]. 中國農村水利水電,2018(8):45-54.

Yi Yinlian, Cui Ningbo, Zhang Qingwen, et al. Research on the comprehensive evaluation model of irrigation operational level[J]. China Rural Water and Hydropower, 2018(8): 45-54. (in Chinese with English abstract)

[26] 杜軍,哈岸英. 寧夏引黃現代化生態灌區頂層設計與建設[J]. 中國水利,2018(7):40-43.

Du Jun, Ha Anying. Top-level design and construction of modernized ecological-friendly irrigation districts by diverting water from Yellow River in Ningxia[J]. China Water Resources, 2018(7): 40-43. (in Chinese with English abstract)

[27] 楊自健. 寧夏吳忠市利通區現代化生態灌區建設實踐探索[J]. 中國水利,2018(23):56-58.

Yang Zijian. Practices of building modernized ecological-friendly irrigation districts in Litong district of Wuzhong city in Ningxia[J]. China Water Resources, 2018(23): 56-58. (in Chinese with English abstract)

[28] 包長征,李廣成,馬廣福,等. 寧夏賀蘭縣耕地地力評價與測土配方施肥[M]. 銀川:陽光出版社,2015.

[29] 吉別克·哈力克巴義,丁建麗,謝力扎提·哈布爾. 干旱區土地鹽漬化時空動態變化及影響因子分析[J]. 水土保持研究,2011,18(4):112-116,278.

Jibieke·Halikebayi, Ding Jianli, Xielizat·Habul. Analysis on tempo-spatial dynamic change and driving factors on soil salinification in arid area[J]. Research of Soil and Water Conservation, 2011, 18(4): 112-116, 278. (in Chinese with English abstract)

[30] 張添佑,王玲,王輝,等. 瑪納斯河流域鹽漬化灌區生態環境遙感監測研究[J]. 生態學報,2017,37(9):3009-3018.

Zhang Tianyou, Wang Ling, Wang Hui, et al. Assessment of soil salinization ecological environment change in the Manas river basin using remote sensing technology[J]. Acta Ecologica Sinica, 2017, 37(9): 3009-3018. (in Chinese with English abstract)

[31] 關元秀,劉高煥,劉慶生,等. 黃河三角洲鹽堿地遙感調查研究[J]. 遙感學報,2001(1):46-52,86.

Guan Yuanxiu, Liu Gaohuan, Liu Qingsheng, et al. The study of salt-affected soils in the Yellow river delta based on remote sensing[J]. National Remote Sensing Bulletin, 2001(1): 46-52, 86. (in Chinese with English abstract)

[32] 陳實,徐斌,金云翔,等. 北疆農區土壤鹽漬化遙感監測及其時空特征分析[J]. 地理科學,2015,35(12):1607-1615.

Chen Shi, Xu Bin, Jin Yunxiang, et al. Remote sensing monitoring and spatial-temporal characteristics analysis of soil salinization in agricultural area of northern Xinjiang[J]. Scientia Geographica Sinica, 2015, 35(12): 1607-1615. (in Chinese with English abstract)

[33] 李丹,吳保生,陳博偉,等. 基于衛星遙感的水體信息提取研究進展與展望[J]. 清華大學學報:自然科學版,2020,60(2):147-161.

Li Dan, Wu Baosheng, Chen Bowei, et al. Review of water body information extraction based on satellite remote sensing[J]. Journal of Tsinghua University: Science and Technology, 2020, 60(2): 147-161. (in Chinese with English abstract)

[34] 童慶禧. 中國典型地物波譜及其特征分析[M]. 北京:科學出版社,1990.

[35] 吳志勇,程丹丹,何海,等. 綜合干旱指數研究進展[J]. 水資源保護,2021,37(1):36-45.

Wu Zhiyong, Cheng Dandan, He Hai, et al. Research progress of composite drought index[J]. Water Resources Protection, 2021, 37(1): 36-45. (in Chinese with English abstract)

[36] 黃健熙,李荔,張超,等. 基于遙感蒸散發數據的耕地灌溉保證能力評價方法[J]. 農業工程學報,2015,31(5):100-106.

Huang Jianxi, Li Li, Zhang Chao, et al. Evaluation of cultivated land irrigation guarantee capability based on remote sensing evapotranspiration data[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(5): 100-106. (in Chinese with English abstract)

[37] 孔令健,章啟兵. 基于層次分析法的清流河健康綜合評價[J]. 安徽農學通報,2018,24(19):105-107,113.

Kong Lingjian, Zhang Qibing. Comprehensive health evaluation of Qingliu river based on analytic hierarchy process[J]. Anhui Agricultural Science Bulletin, 2018, 24(19): 105-107, 113. (in Chinese with English abstract)

[38] 蔣光昱,王忠靜,索瀅. 西北典型節水灌溉技術綜合性能的層次分析與模糊綜合評價[J]. 清華大學學報:自然科學版,2019,59(12):981-989.

Jiang Guangyu, Wang Zhongjing, Suoying. Hierarchical analysis and fuzzy evaluation of comprehensive performance of typical water-saving irrigation techniques in Northwest China[J]. Journal of Tsinghua University: Science and Technology, 2019, 59(12): 981-989. (in Chinese with English abstract)

[39] 王振寶,徐海濤,趙天力,等. 水利現代化評價指標體系及評價方法研究[J].中國水利水電科學研究院學報,2013,11(3):236-240.

Wang Zhenbao, Xu Haitao, Zhao Tianli, et al. Evaluation index system and method of water modernization[J]. Journal of China Institute of Water Resources and Hydropower Research, 2013, 13(3): 236-240. (in Chinese with English abstract)

[40] 徐超,高芮,王明田,等. 苗期高溫對草莓果實營養品質影響的模糊綜合評價及模型建立[J]. 中國農業氣象,2020,41(12):785-793.

Xu Chao, Gao Rui, Wang Mingtian, et al. Fuzzy comprehensive evaluation and model establishment of the effect of high temperature in the seedling stage on the nutritional quality of strawberry fruit[J]. Chinese Journal of Agrometeorology, 2020, 41(12): 785-793. (in Chinese with English abstract)

[41] 史海濱,吳迪,閆建文,等. 鹽漬化灌區節水改造后土壤鹽分時空變化規律研究[J]. 農業機械學報,2020,51(2):318-331.

Shi Haibin, Wu Di, Yan Jianwen, et al. Spatial-temporal variation of soil salinity after water saving transformation in salinized irrigation district[J]. Transactions of the Chinese Society for Agricultural Machinery, 2020, 51(2): 318-331. (in Chinese with English abstract)

Health evaluation of the modern ecological irrigation district using multisource data in Helan County, Ningxia, China

Yang Zhenfeng1, Tian Juncang1,2,3※, Feng Kepeng1,2,3, Gong Xue1

(1.,,750021,;2.-,750021,.3.,,750021,)

An ecological irrigation district has emerged as the modern management with healthy ecosystems in precision agriculture. A health assessment of irrigation districts has been one of the most important operations for the decision-making on agricultural water management and pollution. In this study, a novel health evaluation was made on the ecological irrigation district using multisource data. A study area was selected as a part of the Yellow River basin, located at Helan County in Ningxia, China. A knowledge mapping was also applied to systematically sort out the index systems of comprehensive evaluation on the irrigation districts using a VOSviewer platform. The key indicators were then screened with a high co-occurrence frequency. An indicator system was constructed for the comprehensive evaluation of the ecological irrigation districts. Five categories of the primary indicators were divided in the indicator system, including the ecological environment (B1), the modernization of water transmission and distribution (B2), and the modernization of field engineering (B3), the modernization of production and management (B4), and production efficiency (B5). There was also 18 second-grade indicators in total. All indicators were classified into one to four levels of health. With the support of GIS, these indicators were quantified by combining the remote sensing data (LandSat8, MODIS) with research statistics and field trial data. The weight of each indicator was confirmed through Analytic Hierarchy Process (AHP). The evaluation data was divided by the scale of townships, where each administrative division was used as an evaluation unit. Finally, a fuzzy evaluation model was used to implement the comprehensive evaluation on the health level of the ecological irrigation district in the study area. The results show that there were large differences in the indicator B1 among the evaluation units, especially from the indicator of watershed area (C13) and the indicator of water quality in the main drainage ditches (C14), indicating a higher weight. At the same time, the degree of aridity (vegetation index VSWI) (C16) indicated a great potential to save water in the study area. Therefore, it was still encouraged to develop water-saving irrigation. More importantly, there were small differences in the indicators of B2, B3, B4, and B5 among the evaluation units, indicating the consistent upgrading of irrigation districts. Consequently, the highest value of health level (0.81) was found in Nanliangtaizi and Jingxing Farm, whereas, the lowest value was in Hongguang Town (0.69), where all evaluation units were in the range of health level 3 to 2. In addition, 0.325 of the irrigation districts were presented at the first level of health, whereas, 0.272 of the irrigation districts were at the second level, indicating an excellent operating condition. This study can provide a strong reference to efficiently construct the ecological irrigation districts in Northwest China.

remote sensing; ecology; irrigation districts; indicator system; health evaluation, Helan County of Ningxia

2021-07-02

2021-12-09

寧夏回族自治區重點研發計劃重大項目(2018BBF02022);寧夏高等學校一流學科建設項目(NXYLXK2017A03;NXYLXK2021A03))

楊振峰,博士,研究方向為節水灌溉理論與技術。Email:yzf112082@163.com

田軍倉,教授,博士生導師,研究方向為節水灌溉與水資源高效利用。Email:slxtjc@163.com

10.11975/j.issn.1002-6819.2021.24.034

S274

A

1002-6819(2021)-24-0307-10

楊振峰,田軍倉,馮克鵬,等. 基于多源數據的寧夏賀蘭縣現代化生態灌區健康評價[J]. 農業工程學報,2021,37(24):307-316. doi:10.11975/j.issn.1002-6819.2021.24.034 http://www.tcsae.org

Yang Zhenfeng, Tian Juncang, Feng Kepeng, et al. Health evaluation of the modern ecological irrigation district using multisource data in Helan County, Ningxia, China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2021, 37(24): 307-316. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2021.24.034 http://www.tcsae.org

猜你喜歡
現代化生態水平
張水平作品
邊疆治理現代化
“生態養生”娛晚年
保健醫苑(2021年7期)2021-08-13 08:48:02
住進呆萌生態房
學生天地(2020年36期)2020-06-09 03:12:30
生態之旅
加強上下聯動 提升人大履職水平
人大建設(2019年12期)2019-05-21 02:55:32
Hair Highway--Studio Swine 的現代化詮釋
工業設計(2016年5期)2016-05-04 04:00:35
現代化
生態
領導文萃(2015年4期)2015-02-28 09:19:05
國內某1 700 mm熱軋廠現代化改造
上海金屬(2014年4期)2014-12-15 10:40:43
主站蜘蛛池模板: 色综合五月| 亚洲熟女中文字幕男人总站| 国产成人h在线观看网站站| 亚洲欧洲日产国码无码av喷潮| 99免费视频观看| 国产一区二区福利| 全部免费毛片免费播放| 亚洲精品欧美日韩在线| 欧美 亚洲 日韩 国产| 国产精品乱偷免费视频| 爱爱影院18禁免费| 中文一级毛片| 成人在线不卡视频| 国产精品欧美激情| 国产成人午夜福利免费无码r| 日韩无码视频播放| av在线5g无码天天| 免费人成网站在线观看欧美| 成人国产精品网站在线看| 好紧太爽了视频免费无码| 91丝袜乱伦| 一本大道视频精品人妻| 国产精品区网红主播在线观看| 久久久成年黄色视频| 国产精品黄色片| 国产成人高清精品免费软件 | 欧美成人一级| 国产91小视频在线观看| 99尹人香蕉国产免费天天拍| 成人福利在线观看| 欧美在线三级| 国内精品九九久久久精品| 国产乱子精品一区二区在线观看| 久久精品这里只有国产中文精品| 国产精品区视频中文字幕| 在线欧美国产| 97影院午夜在线观看视频| 久草视频精品| 青青草国产精品久久久久| 国产免费怡红院视频| 日韩区欧美区| 国产在线观看成人91| 99久久精品久久久久久婷婷| 久久免费观看视频| 在线欧美日韩国产| 国产乱人伦AV在线A| 伊人久久久久久久久久| 国产福利微拍精品一区二区| 小说区 亚洲 自拍 另类| 青青草91视频| 亚洲色欲色欲www网| 国产美女自慰在线观看| 国产精品xxx| 亚洲人成影院在线观看| 99re这里只有国产中文精品国产精品 | 成人综合网址| 无码AV高清毛片中国一级毛片 | 色哟哟精品无码网站在线播放视频| 国产一级视频久久| 99热这里只有成人精品国产| 亚洲天堂网在线播放| 玖玖精品视频在线观看| 国产麻豆另类AV| 中文无码影院| 国产亚洲精品yxsp| 国产亚洲精品97AA片在线播放| 99精品伊人久久久大香线蕉| 狠狠色综合网| 国产你懂得| 中文字幕亚洲电影| 国产精品午夜福利麻豆| 一本大道香蕉中文日本不卡高清二区| 99re在线视频观看| 欧美精品成人| 亚洲欧美人成电影在线观看| 日本五区在线不卡精品| 中文精品久久久久国产网址| 不卡午夜视频| 性视频久久| 国产一区在线观看无码| 国产精品对白刺激| 国产黑丝视频在线观看|