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

    農(nóng)田灌排系統(tǒng)雨洪調蓄模擬及分析

    2024-12-31 00:00:00龔來紅才碩趙樹君王瓊
    安徽農(nóng)業(yè)科學 2024年24期

    摘要為探究農(nóng)田、溝渠、塘堰、小型水庫等組成的農(nóng)田灌溉排水系統(tǒng)(簡稱“田-溝-塘-庫”系統(tǒng))對雨洪調蓄的影響,以江西省灌溉試驗中心站研究基地的農(nóng)田灌溉排水系統(tǒng)為研究對象,運用環(huán)境流體動力學模型(environmentalfluiddynamicscode,EFDC)構建典型小型灌溉排水系統(tǒng)的水量平衡模型,評估不同降雨頻率下農(nóng)田、塘堰和水庫采用不同水管理模式對排水口洪峰流量、排水量的影響。結果表明:系統(tǒng)排水溝在正常水管理狀態(tài)下,通過蓄積一定的雨洪資源可以削減排水口洪峰流量29.55%,減少排水量31.08%;與淹水灌溉模式相比,間歇灌溉模式系統(tǒng)排水量減少12.69%;當塘堰控制水深為0.6m時,系統(tǒng)排水量最小。模擬結果表明,農(nóng)田灌排系統(tǒng)可以削減徑流峰值、減少排水量,具有良好的蓄雨調洪能力。

    關鍵詞灌溉排水;環(huán)境流體動力學模型;雨洪調蓄;情景分析

    中圖分類號S27"文獻標識碼A

    文章編號0517-6611(2024)24-0177-07

    doi:10.3969/j.issn.0517-6611.2024.24.038

    開放科學(資源服務)標識碼(OSID):

    SimulationandAnalysisofRainandFloodRegulationandStorageinFarmlandIrrigationandDrainageSystem

    GONGLai-hong1,CAIShuo2,ZHAOShu-jun3etal

    (1.JiangxiAuthorityofWaterConservancyProjectoftheGanfuPlain,Nanchang,Jiangxi330096;2.JiangxiCentralStationofIrrigationExperiment,Nanchang,Jiangxi330201;3.ChinaThreeGorgesUniversity,Yichang,Hubei443002)

    AbstractInordertoexploretheinfluencesoffarmlandirrigationanddrainagesystemcomposedoffarmland,ditch,pondweirandsmallreservoir(referredtoas“farmland-ditch-pond-reservoir”system)onrainwaterregulationandstorage,wetookfarmlandirrigationanddrainagesystemoftheresearchbaseinJiangxiCentralStationofIrrigationExperimentastheresearchobject,andusedtheenvironmentalfluiddynamicscode(EFDC)modeltoconstructwaterbalancemodelofthetypicalsmallirrigationanddrainagesystem,andevaluatetheinfluencesofdifferentwatermanagementmodeloffarmlandandpondweir,andreservoironoutletpeakdischargeandthe"drainagevolumeofthesystemunderdifferentrainfallfrequencies.Theresultsshowedthatthesystemcouldreducethepeakoutletrunoffby29.55%,andreducethedrainagevolumeby31.08%underthenormalwatermanagementstatus.Thedrainagevolumeoftheintermittentirrigationmodewasreducedby12.69%comparedwiththatoffloodedirrigationmode.Thedrainagevolumeofthesystemwastheminimumwhenthecontrolledwaterdepthofthepondweirwas0.6m.Thesimulationresultsshowedthatthefarmlandirrigationanddrainagesystemcouldreducethepeakrunoff"anddrainagevolume,andithadagoodrainwaterstorageandfloodregulationcapacity.

    KeywordsIrrigationanddrainage;EFDC;Rainandfloodregulationandstorage;Scenarioanalysis

    稻田、排水溝、塘堰和水庫共同組成了我國南方灌區(qū)典型的農(nóng)田灌溉排水系統(tǒng),在保障糧食生產(chǎn)上發(fā)揮了重要的作用。稻田在一定程度上具備蓄雨調洪的功能。Smiley等[1根據(jù)稻田對雨洪資源的調節(jié)作用,提出了稻田的防洪蓄雨功能。Matsuno等[2對稻田功能價值進行了評價,發(fā)現(xiàn)稻田的防洪蓄雨功能得到了廣泛研究。黃璜3通過調查研究發(fā)現(xiàn)湖南省內稻田具備“隱形水庫”的功能。王傳娟等4通過構建稻田水量平衡模型,研究發(fā)現(xiàn)稻田能有效蓄存雨水減緩澇災,且稻田不同灌溉模式中淺濕調控灌溉模式具有更好的節(jié)水、蓄雨效果。塘堰系統(tǒng)也有一定的納洪調蓄能力,鄭祖金等[5在研究塘堰攔蓄雨洪資源的能力時,提出了攔蓄系數(shù),并將其用于衡量塘堰對地表徑流的調蓄能力。Ferrati等[6通過構建水平衡模型,模擬分析了Esteros del Ibera濕地的水位動態(tài)變化過程。李玲君等[7改進了WRSIS(灌溉-排水-濕地綜合管理系統(tǒng))中的濕地系統(tǒng)、灌排系統(tǒng)、田間水位控制系統(tǒng),并將其應用在我國南方水稻灌區(qū),減少了灌溉水量。董斌等[8、魏小華等9在廣西青獅潭灌區(qū)和湖北漳河灌區(qū)的田間試驗表明,農(nóng)田灌排系統(tǒng)削減了灌溉水量,攔截了地表徑流,此外它對水質有凈化作用。此外,通過合理控制灌區(qū)灌排系統(tǒng)的灌排水量,可以重復利用農(nóng)田排水,減輕下游防洪壓力[10,減少農(nóng)田氮磷污染物排放11-12。

    目前,大部分學者[4,13-14以稻田、排水溝、塘堰濕地、水庫及灌溉排水系統(tǒng)的單一個體為研究對象來研究水平衡過程,而關于“田-溝-塘-庫”系統(tǒng)聯(lián)合調度水量的報道較少。對于灌溉排水系統(tǒng),大多采用水文模型或數(shù)值模型,從流域或灌區(qū)角度評價節(jié)水灌溉措施對系統(tǒng)出口水量、水質的影響,而關于水庫管理和濕地水管理對水量影響的研究較少。環(huán)境流體動力學模型(environmental fluid dynamics code,EFDC)是一種兼?zhèn)渌畡恿退|模擬計算的綜合性數(shù)學模型[15-16,主要應用于湖泊、河口、濕地、排水溝的水量、水質、泥沙等過程的水動力模擬17-19,比如Okeechobe湖[20、Tasman and Golden海灣[21、St.Lucie河口[22等,在水庫、濕地和排水溝的水動力、水質和泥沙模擬中也有較好的模擬效果23-26,可作為灌區(qū)農(nóng)田灌排系統(tǒng)納洪減污效應綜合評價研究的方法[27

    筆者選取江西省灌溉試驗中心站內已建成的“田-溝-塘-庫”農(nóng)田灌溉排水系統(tǒng)為研究對象,基于實地調研、實測數(shù)據(jù)資料建立EFDC水量平衡模型,模擬不同情景模式下的水力動力學過程,以系統(tǒng)總排水口的實測和模擬徑流過程評估“田-溝-塘-庫”系統(tǒng)的蓄雨調洪能力,以期為江西乃至南方地區(qū)的雨洪調蓄提供理論參考和實踐指導。

    1研究區(qū)概況與方法

    1.1研究區(qū)概況

    江西省灌溉試驗中心站位于鄱陽湖流域贛撫平原灌區(qū)內,灌區(qū)內農(nóng)業(yè)生產(chǎn)模式主要為雙季水稻,灌區(qū)流域經(jīng)緯度為115°49′~116°46′E,28°24′~29°46′N,平均海拔約20m,南北走向170km,最大寬度約70km。試驗區(qū)土壤多為紅壤性水稻黏土,黏土多以高嶺石-石英-蒙脫石為礦物組成成分,鹽基飽和度較低(10.00%~25.00%)。土壤呈酸性,pH4.50~6.50。

    試驗以江西省灌溉試驗中心站內已建成的“田-溝-塘-庫”農(nóng)田灌溉排水系統(tǒng)(圖1)為研究對象。系統(tǒng)中,農(nóng)田區(qū)域地勢呈北高南低,主農(nóng)田排水溝位于稻田東側,區(qū)域內田塊按水稻不同處理分別管理。該研究涉及面積8400m2的試驗稻田4塊,長度110m的排水溝1條,面積460m2的塘堰3個。稻田排水進入排水溝后輸運至塘堰,最終匯入面積680m2的水庫進行儲水調蓄。

    1.2數(shù)據(jù)來源

    以水稻全生育期稻田水量變化過程為研究對象,通過試驗和模型相結合的方法開展研究,所需資料來源于江西省灌溉試驗中心站長系列灌溉試驗和現(xiàn)場實地調研。其中,降雨資料來源于研究基地氣象站1978—2019年的氣象資料。50%年單日最大降雨量為111.70mm,出現(xiàn)在2008年5月28日,并以該次降雨為典型降雨條件開展模擬。

    EFDC模型水動力模塊需要搜集地形資料、流量資料、氣象要素資料、水工建筑物資料、水管理措施等資料;輸出項包括三維流速、水深、流量等隨時間變化的變量。

    稻田的水平衡過程由試驗場內的需水量測坑測得,降雨量數(shù)據(jù)來自試驗場內的標準氣象場,蒸發(fā)蒸騰量、排水量、灌溉用水量由測坑實測。稻田排水處、排水溝、塘堰和水庫的進水、出水位置均設矩形薄壁堰,用于水量計量;在薄壁堰處設置HOBO-U20-001-04型號水位計,記錄數(shù)據(jù)間隔5min。塘堰和水庫通過管道連接,水庫的出水口為管道,管道上安裝水表,逐日計量塘堰至水庫進水量和水庫排水量。

    1.3模型原理

    EFDC模型水動力過程的計算采用Sigma坐標下的動量和質量守恒原理及方程,主控方程采用外模與內模分裂的方法進行數(shù)值求解,網(wǎng)格內和網(wǎng)格間可以在水量平衡過程計算的同時,對動力學過程進行計算,能夠滿足農(nóng)田灌排系統(tǒng)最小單元及日尺度上的徑流過程計算。喬斌[28、王曉玲等29利用EFDC模型對研究區(qū)內的農(nóng)田、排水溝、塘堰進行了概化,研究了農(nóng)田排水溝和塘堰降雨徑流過程及其氮磷凈化效果。該研究將小尺度的灌溉排水系統(tǒng)概化為一個有自由水面的大水箱,將稻田、塘堰、水庫概化為小水箱,小水箱內又劃分成多處網(wǎng)格,可開展水量平衡計算,系統(tǒng)內各單元的水動力計算過程適用于排水溝、系統(tǒng)出口徑流模擬。該研究選取EFDC模擬農(nóng)田灌排系統(tǒng)對雨洪調蓄的影響。

    研究區(qū)采用正交笛卡爾網(wǎng)格劃分,二維地形數(shù)據(jù)、高程數(shù)據(jù)由全站儀測得。稻田區(qū)域面積8400m2,長170.00m,寬68.00m,矩形網(wǎng)格尺寸大小為5.70m×5.70m;塘堰總面積1380m2,共3個塘堰,矩形網(wǎng)格尺寸0.60m×0.60m;水庫面積680m2,矩形網(wǎng)格尺寸0.60m×0.60m;“田-溝-塘-庫”系統(tǒng)共劃分7320個網(wǎng)格,網(wǎng)格劃分結果如圖2所示。網(wǎng)格劃分后,將高程數(shù)據(jù)導入EFDC模型,網(wǎng)格區(qū)域變?yōu)槿S網(wǎng)格計算模型。

    1.4模型的構建和率定

    1.4.1模型的構建。模型的構建涉及邊界條件、初始條件、模型運行參數(shù)等。邊界條件:模型的邊界條件需要確定模型流量邊界、建筑物邊界、氣象要素等邊界條件。模型的初始條件包括水位初始條件、底部糙度等初始條件;模型的其他參數(shù)包括水動力計算的參數(shù)、時間步長參數(shù)等。

    1.4.1.1

    邊界條件。①流量邊界包括入流和出流邊界,模擬時以稻田的逐日灌溉流量、逐日排水流量的時間序列為入、出流邊界;②水工建筑物邊界,包含稻田、塘堰、水庫等設置的田埂和稻田-排水溝以及排水溝-塘堰所設置的薄壁堰;③開、閉邊界,包含模型區(qū)域與外界發(fā)生流量交換的部分,該研究中閉邊界為稻田、排水溝、塘堰、水庫等,其設置目的是避免與外界發(fā)生水體交換;④三維網(wǎng)格計算模型的上邊界為自由水面邊界,下邊界為固壁邊界。

    1.4.1.2初始條件。①地形條件。包括稻田、排水溝、塘堰、水庫和地面底層標高,分別為20.19、20.19、19.29、18.59和20.85m。稻田隨生育期的不同有不同的蓄水深度,通過與排水溝連接的薄壁堰堰口高程調整;稻田、排水溝、塘堰、水庫等單元均設置田埂,且田埂高度30cm,高程20.49m。②底部糙度。糙度值與水動力的計算過程有關,其選取對于紊流計算結果的影響較大,根據(jù)經(jīng)驗取為0.01。③水位初始條件。初始時,稻田水深為泡田結束時的水層深度(30mm),3塊塘堰的水深分別為200、400、600mm,水庫為排干狀態(tài),排水溝中無水。

    1.4.1.3

    模型的其他運行參數(shù)。①干濕節(jié)點參數(shù),用于判斷是否參與漫灘分析計算。若超過濕網(wǎng)格高程,參與水動力過程計算;若低于該值,為干網(wǎng)格,不參與計算;②時間步長。當時間步長較長時,模型網(wǎng)格發(fā)生水體溢出;若時間步長較短,則需要較長的計算時間。取0.03s為計算時間步長。

    1.4.2 模型率定。對1978—2019年江西省雙季稻生長期內的長系列降雨資料進行排頻分析,確定2008年為50%降雨頻率水平年。以2008年作為50%典型年進行模擬,對實測值最齊全的2018年開展模型率定。參考EFDC模型參數(shù)敏感性分析相關文獻[30,發(fā)現(xiàn)風速、底部糙度為對水動力過程模擬的不確定性貢獻率較大的2個參數(shù)。選用上述2個參數(shù)進行手動調參,平均風速取灌溉試驗站氣象站實測的逐日風速、風向,底部糙度取0.01。基于實測資料,開展稻田水量平衡與塘堰水量平衡模型驗證。

    1.4.2.1稻田水量平衡。

    稻田水量平衡影響因素主要包含降雨量、騰發(fā)量、滲漏量、灌水量及各單元排水量。滲漏量包含深層滲漏和側向滲漏,該研究中側向滲漏較小,不予考慮。灌水量通過實測數(shù)據(jù)獲得,而排水量通過設定不同條件模擬獲得。水量平衡原理如圖3所示。

    利用EFDC模型進行稻田水量平衡模擬,根據(jù)水稻不同生育期內水分控制方式的差異,在稻田內設置不同的田埂高度,以改變稻田雨后上限,模擬獲得系統(tǒng)排水量。在模型矩形薄壁堰堰口高程設置過程中,無法實現(xiàn)設置堰口高程隨時間序列的變化,必須手動輸入水稻不同生育期內控制堰口高程。故在模型模擬過程中,不同生育期改變堰口高程的同時,利用上一生育期末稻田模擬水深作為下一生育期稻田初始水位,根據(jù)不同生育期稻田水分控制方式模擬水量變化過程。

    該研究利用水位作為驗證稻田水量平衡模型的重要特征。在排水量、滲漏量以及騰發(fā)量等去水量的作用下,田面水層會下降。在降雨量、灌水量等來水量的作用下,稻田水位又會逐漸上升。模型在設置初始條件的過程中,不能直接反映稻田的滲漏水量,故將滲漏量概化為實測的騰發(fā)量,經(jīng)驗證準確有效。對2008年水稻全生育期內不同灌溉模式下稻田實測水位進行模型驗證,其結果如圖4所示。利用SPSS軟件分析2組數(shù)據(jù)的相關性,結果見表1。相對誤差計算公式如下:

    RE=Pt-QtQt×100%

    式中,RE為模型模擬相對誤差,Pt為模擬值,Qt為實測值。

    通過對2008年水稻全生育期內不同灌溉模式下水位實測數(shù)據(jù)和模型模擬計算結果進行對比發(fā)現(xiàn),間歇灌溉模式下2018年水稻全生育期內實測水位與模擬值相對誤差僅0.30%,相關系數(shù)為0.96;淹灌模式下實測水位與模擬值相對誤差為0.19%,相關系數(shù)為0.98。由此可見,模型模擬水位變化趨勢與實測數(shù)據(jù)相一致,模擬平均相對誤差為0.25%,相關系數(shù)均大于0.95,擬合程度較高,說明EFDC模型可用于模擬稻田的水量平衡過程。在模型模擬過程中,前期誤差范圍較大。究其原因,這是因為模型采用冷啟動方式,很難避免帶來起始誤差。在后期模型運行過程中主要受邊界條件的控制,誤差較小。

    1.4.2.2塘堰水量平衡。塘堰水量平衡與稻田有一定的差異,其主要包含大氣降雨量、滲漏量、蒸發(fā)量、承接排水溝來水的入流量和匯入水庫的出流量,形式和結構與稻田相比更為簡單。塘堰降雨量、入流量等來水量以蒸發(fā)、滲漏和排水的形式離開。以2008年實測水位數(shù)據(jù)作為驗證塘堰水量平衡的表征指標,驗證結果如圖5所示。實測水位與模擬值十分接近,平均相對誤差僅4.4%,擬合程度較好,說明EFDC模型可用于模擬塘堰水量變化過程。

    2結果與分析

    基于所構建的“田-溝-塘-庫”EFDC水動力模型,對灌區(qū)“田-溝-塘-庫”農(nóng)田灌排體系開展水動力過程模擬研究。模擬情景模式包括系統(tǒng)不同水位工況、不同灌溉模式和不同塘堰控制水深條件下系統(tǒng)的水量變化過程,研究其變化規(guī)律。

    2.1不同水管理模式對灌排系統(tǒng)納洪效應的影響

    2.1.1不同水管理模式試驗情景設置。

    水稻生育期50%降雨頻率水平年(2008年)日最大降雨出現(xiàn)在5月28日,日降雨量為111.7mm。以該次降雨為平均水平年內日最大降雨,設置3種情景模式,分別模擬稻田、塘堰和水庫在正常水管理、蓄滿狀態(tài)、全部排干等水管理模式采用不同控制水深對排水過程的影響,研究不同水管理模式對系統(tǒng)出口排水過程的影響。由表2可知,當系統(tǒng)處于正常水管理模式時,稻田、塘堰、水庫水深分別為0.03、0.60和1.00m;當系統(tǒng)處于全部蓄滿狀態(tài)時,稻田、塘堰、水庫水深分別為0.08、0.90和1.60m。

    2.1.2不同水管理模式下灌排系統(tǒng)對雨洪的響應。

    由表3可知,在系統(tǒng)正常蓄水狀態(tài)下,在一日模擬降雨過程中,稻田水位上升至雨后上限值為80mm,即水位20.27m,通過薄壁堰的堰口調整,排水溝排水量和排水口洪峰流量分別為495.64m3和0.0062m3/s,較蓄滿狀態(tài)系統(tǒng)下減少排水量31.08%,排水口洪峰流量削減29.55%,同時水庫排水量為0m3,說明系統(tǒng)在正常水位狀態(tài)下同樣具備蓄積一日雨洪資源的能力。當系統(tǒng)處于全部蓄滿狀態(tài)時,系統(tǒng)均在較短時間內達到穩(wěn)定的出流狀態(tài),排水溝、塘堰、水庫排水口洪峰流量分別為0.0088、0.0115和0.0138m3/s,排水量分別為719.19、911.37和1093.46m3。系統(tǒng)全干狀態(tài)下的一日排水過程中,排水溝排水量為384.54m3,相較于全滿狀態(tài)排水量減少334.65m3,而塘堰、水庫均無排水量,說明系統(tǒng)在全干狀態(tài)下完全具備蓄積一日雨洪資源的能力。

    2.2不同塘堰控制水深對灌排系統(tǒng)納洪效應的影響

    2.2.1不同塘堰控制水深試驗情景設置。

    在模擬過程中,通過改變塘堰出水口處堰高,實現(xiàn)對塘堰水深的控制。根據(jù)塘堰最大蓄水深度0.90m,正??刂扑?.60m,設定以0.20m為梯度的塘堰不同控制水深,具體情景設置見表4。根據(jù)不同情景模式設置,研究系統(tǒng)排水量在不同塘堰控制水深下的變化規(guī)律。為方便排水量和洪峰流量的計算,在EFDC模型水庫排水口處設置底部高程1.0m、面積2160m2的虛擬蓄水池(簡稱蓄水池),以便于利用蓄水池水深來計算系統(tǒng)的排水量。

    2.2.2模擬結果分析。

    圖6給出了不同控制水深的情景模式下模擬塘堰、水庫、蓄水池的水位變化。在情景1模式下,當塘堰控制水深為0m時,早稻生育期前期在少量降雨和蒸發(fā)蒸騰的作用下,整個系統(tǒng)水位變化不大;當集中發(fā)生大的降雨而由稻田產(chǎn)生地表徑流進入塘堰時,迅速通過排水管道匯入水庫。當水庫水位達到蓄滿狀態(tài)時,體系中多余的水量排入蓄水池蓄積。由此可見,在系統(tǒng)聯(lián)動蓄水、排水的作用下,可以有效滯后峰值排水流量。

    在不同控制水深情景模擬下,根據(jù)系統(tǒng)水位變化過程計算出系統(tǒng)排水量。在塘堰不同控制水深情景模擬下,系統(tǒng)排水量為3731.43~4135.16m3。其中,系統(tǒng)排水量最大的是控制水深0m的塘堰,達4135.16m3;其次是控制水深0.80m的塘堰,系統(tǒng)排水量為3837.85m3;排水量較小的是正??刂扑?.60m與控制水深0.20m時的塘堰,排水量分別為3731.43和3747.71m3。因此,在水管理上可以通過合理控制塘堰水深,減少系統(tǒng)排水量,從而達到蓄水納洪的目的。

    2.3不同灌溉模式對灌排系統(tǒng)納洪效應的影響

    2.3.1不同灌溉模式試驗情景設置。

    通過研究塘堰不同控制水深對“田-溝-塘-庫”系統(tǒng)納洪效應的影響,進一步開展不同灌溉模式下水稻全生育期內系統(tǒng)的水量模擬試驗,研究其在間歇灌溉和淹水灌溉2種灌溉模式下系統(tǒng)的排水量變化。2種灌溉模式下控制水深見表5。

    2.3.2模擬結果分析。

    在整個生育期,由于降雨量變化較大,利用模型模擬流量值來推算排水量誤差較大,故在模型系統(tǒng)出水口處建蓄水池(長54.00m、寬40.00m)。通過蓄水池水深變化推算系統(tǒng)排水量,模擬整個生育期不同灌溉模式下系統(tǒng)水位、蓄水池水位的變化,如圖7、8所示。

    稻田不同灌溉模式下系統(tǒng)排水量計算結果如表6所示。由表6可知,2種灌溉模式下均有大量排水,間歇灌溉和淹水灌溉模式水庫排水量分別為3617.20和4142.75m3,在間歇灌溉模式下排水量較淹水灌溉模式減少12.69%。分析2種模式下系統(tǒng)水位的變化可知,在淹灌條件下水庫較間歇灌溉更快達到蓄滿狀態(tài),減小了后期蓄水庫容,增大了系統(tǒng)排水量,從而降低了系統(tǒng)蓄雨納洪的能力。以上結果說明間歇灌溉條件下系統(tǒng)具備更好的納洪效應。

    3結論

    基于構建的“田-溝-塘-庫”水動力模型,研究在不同情景模式下系統(tǒng)的納洪效應,通過EFDC模型模擬不同塘堰控制水深(包括正常水管理狀態(tài)、蓄滿狀態(tài)和全干狀態(tài))下的排水量變化,分析其對系統(tǒng)納洪效應的影響以及不同灌溉模式下系統(tǒng)排水量的變化。其主要結論如下:

    (1)“田-溝-塘-庫”系統(tǒng)具備一定的雨洪資源調蓄能力。系統(tǒng)排水溝在正常水管理狀態(tài)下,可通過蓄積一定的雨洪資源削減系統(tǒng)洪峰流量29.55%,減少排水量30.18%,能夠有效降低灌區(qū)洪澇災害。

    (2)塘堰不同控制水深對系統(tǒng)排水量有一定的影響。在正??刂扑顬?.60m時,系統(tǒng)排水量最?。?731.43m3),可以為農(nóng)田水管理提供一定的依據(jù)。

    (3)在間歇灌溉條件下,相較于淹灌模式,“田-溝-塘-庫”系統(tǒng)可減少排水量525.55m3,減少幅度達12.69%,表現(xiàn)出較好的納洪效應。

    參考文獻

    [1]

    SMILEYPC,ALLREDBJ.Differencesinaquaticcommunitiesbetweenwetlandscreatedbyanagriculturalwaterrecyclingsystem[J].Wetlandsecologyandmanagement,2011,19(6):495-505.

    [2]MATSUNOY,NAKAMURAK,MASUMOTOT,etal.ProspectsformultifunctionalityofpaddyricecultivationinJapanandothercountriesinmonsoonAsia[J].Paddywaterenvironment,2006,4(4):189-197.

    [3]黃璜.湖南境內隱形水庫與水庫的集雨功能[J].湖南農(nóng)業(yè)大學學報(自然科學版),1997,23(6):499-503.

    [4]王傳娟,王少麗,陳皓銳,等.稻田水量調控模擬計算及分析[J].中國農(nóng)村水利水電,2016(8):137-143.

    [5]鄭祖金,崔遠來,董斌,等.灌區(qū)塘堰攔蓄地表徑流能力的研究[J].中國農(nóng)村水利水電,2005(1):39-40.

    [6]FERRATIR,CANZIANIGA.AnanalysisofwaterleveldynamicsinEsterosdelIberawetland[J].Ecologicalmodelling,2005,186(1):17-27.

    [7]李玲君,李琳娜,邵孝侯,等.濕地-水塘-地下灌排改進系統(tǒng)的效應[J].水資源保護,2012,28(3):75-78,87.

    [8]董斌,茆智,李新建,等.灌溉-排水-濕地綜合管理系統(tǒng)的引進和改造應用[J].中國農(nóng)村水利水電,2009(11):9-12,15.

    [9]魏小華,王修貴,單真瑩,等.改進的農(nóng)田灌溉排水系統(tǒng)[J].灌溉排水學報,2007,26(S1):6-8.

    [10]劉鳳麗,熊玉江,范樂,等.水稻灌區(qū)降雨徑流特征及其影響因素研究[J].中國農(nóng)村水利水電,2022(5):95-100.

    [11]華克驥,何軍,張宇航,等.不同灌溉和施肥方式對稻田土壤氮、磷遷移轉化的影響[J].灌溉排水學報,2022,41(7):35-43.

    [12]伯彥萍,張春雷,孫一迪,等.不同灌溉方式下稻田氮磷排放特征[J].節(jié)水灌溉,2022(9):93-99,107.

    [13]奚歌,劉紹民,賈立.黃河三角洲濕地蒸散量與典型植被的生態(tài)需水量[J].生態(tài)學報,2008(11):5356-5369.

    [14]王帥.巴南區(qū)雙河口鎮(zhèn)采空區(qū)稻田蓄水現(xiàn)狀及建議[J].南方農(nóng)業(yè),2015,9(25):46-47,67.

    [15]HAMRICKJM.Athree-dimensionalenvironmentalfluiddynamicscomputercode:Theoreticalandcomputationalaspects[R].Gloucester,Massachusetts:VirginiaInstituteofMarineScience,theCollegeofWilliamandMary,1992:1-60.

    [16]陸耀烽,丁志斌,黎煒,等.近岸海域水質模型研究現(xiàn)狀及展望[J].海洋科學,2020,44(2):161-170.

    [17]陳異暉.基于EFDC模型的滇池水質模擬[J].云南環(huán)境科學,2005,24(4):28-30,46.

    [18]謝銳,吳德安,嚴以新,等.EFDC模型在長江口及相鄰海域三維水流模擬中的開發(fā)應用[J].水動力學研究與進展(A輯),2010,25(2):165-174.

    [19]王翠,孫英蘭,張學慶.基于EFDC模型的膠州灣三維潮流數(shù)值模擬[J].中國海洋大學學報(自然科學版),2008,38(5):833-840.

    [20]JINKR,JIZG.Calibrationandverificationofaspectralwind:WavemodelforLakeOkeechobee[J].Oceanengineering,2001,28(5):571-584.

    [21]TUCKEYBJ,GIBBSMT,KNIGHTBR,etal.TidalcirculationinTasmanandgoldenbays:Implicationsforriverplumebehaviour[J].NewZealandjournalofmarineandfreshwaterresearch,2006,40(2):305-324.

    [22]JIZG,HUGD,SHENJ,etal.Three-dimensionalmodelingofhydrodynamicprocessesintheSt.LucieEstuary[J].Estuarine,coastalandshelfscience,2007,73(1/2):188-200.

    [23]徐婉珍.基于EFDC的潘家口水庫富營養(yǎng)化模擬[J].安徽農(nóng)業(yè)科學,2015,43(34):115-118,290.

    [24]王敏,姜利兵,黃海真,等.大型水庫水質數(shù)值模擬及應用[J].環(huán)境影響評價,2021,43(1):47-51.

    [25]王佳兒,崔遠來,張博超,等.基于EFDC模型的表面流人工濕地水力性能和凈化效果影響因素分析[J].農(nóng)業(yè)工程學報,2024,40(7):172-181.

    [26]WUDA,LIX,WANGYG.StabilityofthestratificationofwaterbodiesintheNorthPassageoftheYangtzeRiverEstuarybasedontheEFDCmodel[J].JournalofOceanUniversityofChina,2024,23(1):23-32.

    [27]時紅,才碩,孫占學,等.田-溝-塘系統(tǒng)納洪減污效應研究進展[J].節(jié)水灌溉,2024(1):16-22.

    [28]喬斌.農(nóng)田生態(tài)溝渠對稻田降雨徑流氮磷的去除實驗與模擬研究[D].天津:天津大學,2015.

    [29]王曉玲,張福超,李松敏,等.穩(wěn)定塘內稻田降雨徑流中氮磷輸移的三維數(shù)值模擬[J].天津大學學報(自然科學與工程技術版),2019,52(6):649-660.

    [30]李一平,唐春燕,余鐘波,等.大型淺水湖泊水動力模型不確定性和敏感性分析[J].水科學進展,2012,23(2):271-277.

    亚洲,欧美,日韩| 中文字幕熟女人妻在线| 国产真实乱freesex| 中文字幕免费在线视频6| 91在线观看av| 午夜影院日韩av| 欧美一区二区国产精品久久精品| 成人av在线播放网站| 国产色爽女视频免费观看| 757午夜福利合集在线观看| 神马国产精品三级电影在线观看| 亚洲性夜色夜夜综合| 久久久国产成人免费| 精品久久久久久久久久久久久| 在线十欧美十亚洲十日本专区| 国产男靠女视频免费网站| 国产亚洲欧美在线一区二区| 最近在线观看免费完整版| 日韩欧美国产在线观看| 欧美三级亚洲精品| 天堂动漫精品| 亚洲欧美日韩高清专用| 成人欧美大片| 此物有八面人人有两片| 18禁在线播放成人免费| 国产v大片淫在线免费观看| 亚洲无线观看免费| 少妇人妻一区二区三区视频| 免费人成视频x8x8入口观看| 别揉我奶头 嗯啊视频| 国产毛片a区久久久久| 国产午夜精品久久久久久一区二区三区 | 十八禁国产超污无遮挡网站| 观看免费一级毛片| 一级黄色大片毛片| 精品乱码久久久久久99久播| 国产欧美日韩一区二区精品| 18+在线观看网站| 久久热精品热| 听说在线观看完整版免费高清| 可以在线观看的亚洲视频| 国内精品久久久久精免费| 亚洲成人免费电影在线观看| 天堂√8在线中文| 亚洲五月婷婷丁香| 成人一区二区视频在线观看| 亚洲一区二区三区色噜噜| 亚洲av免费在线观看| 嫁个100分男人电影在线观看| av黄色大香蕉| 动漫黄色视频在线观看| 日本五十路高清| 国内揄拍国产精品人妻在线| 免费观看人在逋| 色播亚洲综合网| 国产成人啪精品午夜网站| 又黄又爽又免费观看的视频| 美女大奶头视频| 天美传媒精品一区二区| 蜜桃久久精品国产亚洲av| 国内精品一区二区在线观看| 亚洲激情在线av| 国模一区二区三区四区视频| 高潮久久久久久久久久久不卡| 动漫黄色视频在线观看| 女生性感内裤真人,穿戴方法视频| 十八禁人妻一区二区| 中文在线观看免费www的网站| 在线a可以看的网站| 九九久久精品国产亚洲av麻豆| 国内毛片毛片毛片毛片毛片| 亚洲欧美精品综合久久99| 天天一区二区日本电影三级| 色综合亚洲欧美另类图片| 97热精品久久久久久| 亚洲第一电影网av| 午夜福利在线观看免费完整高清在 | 亚洲精品456在线播放app | 男人舔女人下体高潮全视频| 日本一二三区视频观看| 丁香六月欧美| 亚洲欧美日韩无卡精品| 亚洲精华国产精华精| 黄色丝袜av网址大全| 小说图片视频综合网站| 国产精品三级大全| 夜夜夜夜夜久久久久| 欧美日韩国产亚洲二区| 久久这里只有精品中国| www.www免费av| 久久国产精品人妻蜜桃| 国产精品综合久久久久久久免费| 欧美性猛交╳xxx乱大交人| 国产亚洲欧美98| 久久久久久久久久成人| 久久久成人免费电影| 99riav亚洲国产免费| 十八禁国产超污无遮挡网站| 日韩欧美 国产精品| 精品久久久久久久久av| 99热这里只有是精品在线观看 | 亚洲欧美清纯卡通| 午夜日韩欧美国产| 一区福利在线观看| 男人和女人高潮做爰伦理| 免费人成在线观看视频色| 久久热精品热| av视频在线观看入口| 久久精品国产亚洲av天美| 全区人妻精品视频| 国产精品爽爽va在线观看网站| 国产日本99.免费观看| 久久精品国产亚洲av涩爱 | 亚洲在线自拍视频| 国产野战对白在线观看| 成人无遮挡网站| 真实男女啪啪啪动态图| 国产精品影院久久| 亚洲黑人精品在线| 国产色婷婷99| 亚洲精品在线观看二区| 亚洲精华国产精华精| 亚洲人成网站高清观看| 日本与韩国留学比较| 91在线精品国自产拍蜜月| 亚洲精品乱码久久久v下载方式| 在线a可以看的网站| 波多野结衣高清作品| 国产麻豆成人av免费视频| 亚洲国产日韩欧美精品在线观看| 国产高清激情床上av| 亚洲五月婷婷丁香| 亚洲成av人片在线播放无| 午夜福利18| 69人妻影院| 久久久久精品国产欧美久久久| 三级国产精品欧美在线观看| 欧美日韩乱码在线| 午夜久久久久精精品| 精品一区二区三区人妻视频| 搡老妇女老女人老熟妇| 99热这里只有是精品50| 日韩精品青青久久久久久| 最近中文字幕高清免费大全6 | 国产激情偷乱视频一区二区| 日韩中字成人| 久久久国产成人免费| 韩国av一区二区三区四区| 少妇丰满av| 国产精华一区二区三区| 免费av观看视频| 精华霜和精华液先用哪个| 久久午夜福利片| 脱女人内裤的视频| 色噜噜av男人的天堂激情| 男插女下体视频免费在线播放| 丁香六月欧美| 国产一区二区在线观看日韩| 久久九九热精品免费| 99久久99久久久精品蜜桃| avwww免费| 亚洲成av人片免费观看| 国产精品不卡视频一区二区 | 国产色婷婷99| av国产免费在线观看| 午夜老司机福利剧场| 国产av麻豆久久久久久久| 久久国产乱子伦精品免费另类| 夜夜爽天天搞| 亚洲国产精品合色在线| 丰满人妻一区二区三区视频av| bbb黄色大片| 久久久久国产精品人妻aⅴ院| 免费电影在线观看免费观看| 99热6这里只有精品| 757午夜福利合集在线观看| 欧美黄色片欧美黄色片| 国产v大片淫在线免费观看| 成人特级黄色片久久久久久久| 午夜福利免费观看在线| 午夜久久久久精精品| 亚洲精品亚洲一区二区| 亚洲国产精品成人综合色| 最近最新中文字幕大全电影3| 男人舔奶头视频| 性插视频无遮挡在线免费观看| 国产欧美日韩精品一区二区| 久久精品国产亚洲av香蕉五月| 国产淫片久久久久久久久 | 一个人免费在线观看的高清视频| 级片在线观看| 亚洲成人中文字幕在线播放| 简卡轻食公司| 少妇高潮的动态图| 舔av片在线| 国产精品永久免费网站| 国产精品久久久久久亚洲av鲁大| 免费观看的影片在线观看| 亚洲午夜理论影院| 美女高潮喷水抽搐中文字幕| 人妻久久中文字幕网| 国产免费一级a男人的天堂| 一个人免费在线观看电影| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产精品嫩草影院av在线观看 | 欧美三级亚洲精品| 色播亚洲综合网| 可以在线观看毛片的网站| www日本黄色视频网| 欧美成狂野欧美在线观看| 国产免费男女视频| 有码 亚洲区| 亚洲avbb在线观看| 村上凉子中文字幕在线| 日本一二三区视频观看| 欧美不卡视频在线免费观看| 我要看日韩黄色一级片| 高潮久久久久久久久久久不卡| 18禁黄网站禁片免费观看直播| 给我免费播放毛片高清在线观看| 国产精品98久久久久久宅男小说| 男女之事视频高清在线观看| 亚洲中文日韩欧美视频| 精品午夜福利在线看| 久久草成人影院| 国内精品美女久久久久久| 天美传媒精品一区二区| 中文字幕av在线有码专区| 黄色女人牲交| 亚洲精品456在线播放app | 久久久久久国产a免费观看| 日韩欧美在线二视频| 亚洲精品亚洲一区二区| 麻豆国产av国片精品| 欧美日韩瑟瑟在线播放| 噜噜噜噜噜久久久久久91| 欧美日韩亚洲国产一区二区在线观看| 婷婷精品国产亚洲av在线| 一边摸一边抽搐一进一小说| 亚洲国产精品合色在线| 国产精品女同一区二区软件 | 99久久久亚洲精品蜜臀av| 久久6这里有精品| 国产在线男女| 五月玫瑰六月丁香| 日韩欧美国产在线观看| 精品国产三级普通话版| 91麻豆av在线| 久久人人精品亚洲av| 亚洲欧美清纯卡通| 久久国产精品影院| 欧美zozozo另类| 18+在线观看网站| 欧美黄色片欧美黄色片| 精品久久久久久成人av| 美女大奶头视频| 日韩有码中文字幕| 亚洲中文字幕一区二区三区有码在线看| 亚洲欧美日韩卡通动漫| 搞女人的毛片| 俄罗斯特黄特色一大片| 色综合婷婷激情| 亚洲av免费在线观看| 欧美日韩亚洲国产一区二区在线观看| 国产真实伦视频高清在线观看 | 国产伦精品一区二区三区四那| 久久99热这里只有精品18| 欧美色欧美亚洲另类二区| 国内精品美女久久久久久| 夜夜爽天天搞| 少妇高潮的动态图| 91在线精品国自产拍蜜月| 赤兔流量卡办理| 免费电影在线观看免费观看| 亚洲av不卡在线观看| 国产乱人视频| 亚洲人成电影免费在线| 一区二区三区激情视频| 精品久久久久久久末码| 黄色一级大片看看| 性色av乱码一区二区三区2| 麻豆一二三区av精品| 香蕉av资源在线| 欧美黄色淫秽网站| av在线观看视频网站免费| 最近在线观看免费完整版| 久久国产精品影院| 日本 欧美在线| 欧美激情国产日韩精品一区| 国内毛片毛片毛片毛片毛片| 欧美日本视频| 色综合站精品国产| av专区在线播放| 国产精品免费一区二区三区在线| 欧美日韩亚洲国产一区二区在线观看| 搡老岳熟女国产| 久久久久性生活片| 久久欧美精品欧美久久欧美| 国产麻豆成人av免费视频| 午夜福利在线观看免费完整高清在 | 国产精品不卡视频一区二区 | 欧美日韩瑟瑟在线播放| 偷拍熟女少妇极品色| 99精品在免费线老司机午夜| 国产av不卡久久| 级片在线观看| 亚洲av.av天堂| 亚洲美女搞黄在线观看 | 97碰自拍视频| 熟女人妻精品中文字幕| 老熟妇仑乱视频hdxx| 十八禁人妻一区二区| 成年版毛片免费区| 毛片一级片免费看久久久久 | 国产一区二区三区视频了| 最近最新免费中文字幕在线| 一个人看的www免费观看视频| 亚洲av第一区精品v没综合| 99久久99久久久精品蜜桃| 国产一区二区在线av高清观看| 国产又黄又爽又无遮挡在线| 丰满人妻一区二区三区视频av| 亚洲精品乱码久久久v下载方式| 五月伊人婷婷丁香| 成人亚洲精品av一区二区| 人人妻,人人澡人人爽秒播| 99久久无色码亚洲精品果冻| 国产成人啪精品午夜网站| 亚洲人与动物交配视频| 国产精品一区二区免费欧美| 尤物成人国产欧美一区二区三区| 国产精品久久久久久亚洲av鲁大| 一夜夜www| 丁香六月欧美| 国产蜜桃级精品一区二区三区| 在线免费观看不下载黄p国产 | 国产伦一二天堂av在线观看| 免费高清视频大片| 久久精品国产亚洲av天美| 亚洲自拍偷在线| 亚洲人成电影免费在线| 欧美日韩综合久久久久久 | 久久亚洲精品不卡| 国产成人aa在线观看| 在线十欧美十亚洲十日本专区| 国产亚洲欧美在线一区二区| 又爽又黄a免费视频| 亚洲精品一区av在线观看| 禁无遮挡网站| 免费无遮挡裸体视频| 久久精品影院6| 欧洲精品卡2卡3卡4卡5卡区| 成人毛片a级毛片在线播放| 51午夜福利影视在线观看| 国产探花在线观看一区二区| 亚洲片人在线观看| 久久久久性生活片| 国产精品久久久久久久电影| 国产黄a三级三级三级人| 成熟少妇高潮喷水视频| 成人国产一区最新在线观看| 日本一二三区视频观看| 欧美绝顶高潮抽搐喷水| 日本精品一区二区三区蜜桃| 在线免费观看不下载黄p国产 | 91午夜精品亚洲一区二区三区 | 欧美黑人巨大hd| 欧美+日韩+精品| 亚洲精品粉嫩美女一区| 成年女人毛片免费观看观看9| 99久久九九国产精品国产免费| 久久精品影院6| 成年免费大片在线观看| 成人国产综合亚洲| 亚洲激情在线av| 日韩亚洲欧美综合| 国产亚洲精品久久久久久毛片| 嫩草影院入口| .国产精品久久| 日韩av在线大香蕉| 99久久九九国产精品国产免费| 岛国在线免费视频观看| 亚洲最大成人中文| 日日干狠狠操夜夜爽| 一本综合久久免费| 免费观看人在逋| 在线观看舔阴道视频| 色综合站精品国产| 波多野结衣巨乳人妻| 99在线人妻在线中文字幕| 国产在线精品亚洲第一网站| 男女床上黄色一级片免费看| 国产成人a区在线观看| 免费电影在线观看免费观看| 国产大屁股一区二区在线视频| 国产伦精品一区二区三区视频9| 精品日产1卡2卡| 成人三级黄色视频| 国产v大片淫在线免费观看| 久久人人精品亚洲av| 亚洲一区高清亚洲精品| 亚洲国产精品sss在线观看| 国产亚洲欧美98| av中文乱码字幕在线| 国产中年淑女户外野战色| 国内少妇人妻偷人精品xxx网站| 好男人在线观看高清免费视频| 日韩中文字幕欧美一区二区| 久久性视频一级片| 色哟哟哟哟哟哟| 欧美成狂野欧美在线观看| 三级国产精品欧美在线观看| 午夜精品久久久久久毛片777| 色哟哟哟哟哟哟| 午夜精品在线福利| 一进一出抽搐gif免费好疼| 国产欧美日韩一区二区精品| 亚洲一区二区三区色噜噜| 国产亚洲精品久久久com| 国产老妇女一区| 一边摸一边抽搐一进一小说| 欧美三级亚洲精品| 1024手机看黄色片| 午夜福利高清视频| 97超级碰碰碰精品色视频在线观看| 亚洲真实伦在线观看| 一级黄片播放器| 婷婷精品国产亚洲av在线| aaaaa片日本免费| 日本五十路高清| 久久99热6这里只有精品| 最好的美女福利视频网| 蜜桃亚洲精品一区二区三区| eeuss影院久久| 欧美极品一区二区三区四区| 国产伦精品一区二区三区视频9| av黄色大香蕉| 五月玫瑰六月丁香| 午夜免费男女啪啪视频观看 | 国产视频一区二区在线看| 亚洲精品久久国产高清桃花| 亚洲电影在线观看av| 日日夜夜操网爽| 亚洲精品色激情综合| 淫妇啪啪啪对白视频| 十八禁网站免费在线| 久久久久久久亚洲中文字幕 | 脱女人内裤的视频| 欧美日韩中文字幕国产精品一区二区三区| 国产亚洲精品综合一区在线观看| 欧美性猛交╳xxx乱大交人| 99国产精品一区二区三区| 一区二区三区四区激情视频 | 久久99热这里只有精品18| www.熟女人妻精品国产| 亚洲精品色激情综合| 免费在线观看日本一区| 色综合站精品国产| 欧美日韩国产亚洲二区| 女人十人毛片免费观看3o分钟| 亚洲精品一卡2卡三卡4卡5卡| 97热精品久久久久久| 成人特级黄色片久久久久久久| 成年人黄色毛片网站| 身体一侧抽搐| 国产美女午夜福利| 久久人人精品亚洲av| 亚洲第一欧美日韩一区二区三区| 久久亚洲精品不卡| 九色国产91popny在线| 亚州av有码| 午夜视频国产福利| 丝袜美腿在线中文| 免费在线观看影片大全网站| 女人被狂操c到高潮| www.色视频.com| 美女被艹到高潮喷水动态| 观看美女的网站| 国产老妇女一区| bbb黄色大片| 3wmmmm亚洲av在线观看| 久久久精品大字幕| 欧美日韩乱码在线| 久久精品91蜜桃| 一级作爱视频免费观看| 18禁黄网站禁片免费观看直播| 日韩欧美国产在线观看| or卡值多少钱| 波野结衣二区三区在线| 国产精品美女特级片免费视频播放器| 亚洲成人中文字幕在线播放| 淫妇啪啪啪对白视频| 99国产综合亚洲精品| 在线a可以看的网站| 一进一出抽搐动态| 久久精品久久久久久噜噜老黄 | 亚洲精品亚洲一区二区| 欧洲精品卡2卡3卡4卡5卡区| 午夜福利在线观看免费完整高清在 | 国产黄片美女视频| 国产真实乱freesex| 久久国产乱子免费精品| 国产成人a区在线观看| 欧美一区二区亚洲| 日韩中文字幕欧美一区二区| 亚洲成a人片在线一区二区| 国产主播在线观看一区二区| 欧美潮喷喷水| 欧美激情在线99| 悠悠久久av| 国产v大片淫在线免费观看| 亚洲自偷自拍三级| 老司机深夜福利视频在线观看| 18美女黄网站色大片免费观看| 好男人电影高清在线观看| 特大巨黑吊av在线直播| 99在线人妻在线中文字幕| 伊人久久精品亚洲午夜| 毛片女人毛片| 97人妻精品一区二区三区麻豆| 亚洲av成人精品一区久久| 免费观看人在逋| 亚洲 欧美 日韩 在线 免费| 又爽又黄无遮挡网站| 精品乱码久久久久久99久播| 亚洲一区二区三区不卡视频| 美女cb高潮喷水在线观看| 亚洲美女黄片视频| 亚洲在线自拍视频| 51午夜福利影视在线观看| 精品一区二区三区人妻视频| 色综合欧美亚洲国产小说| 欧美黑人欧美精品刺激| 内射极品少妇av片p| 精品欧美国产一区二区三| 亚洲国产精品999在线| 久久久精品大字幕| 一区二区三区四区激情视频 | 精品久久久久久久末码| 神马国产精品三级电影在线观看| 黄色丝袜av网址大全| 舔av片在线| 亚洲人成网站在线播放欧美日韩| 最新在线观看一区二区三区| 日本精品一区二区三区蜜桃| 少妇熟女aⅴ在线视频| 精品一区二区三区视频在线| 在线观看一区二区三区| 成人精品一区二区免费| 一边摸一边抽搐一进一小说| av在线观看视频网站免费| 久久人妻av系列| 深爱激情五月婷婷| 身体一侧抽搐| 亚洲人成网站在线播| 九九热线精品视视频播放| 在线观看66精品国产| aaaaa片日本免费| 日韩欧美三级三区| 国产免费男女视频| 别揉我奶头~嗯~啊~动态视频| 午夜a级毛片| 午夜日韩欧美国产| 日本一二三区视频观看| 18禁黄网站禁片免费观看直播| 欧美最黄视频在线播放免费| 老鸭窝网址在线观看| 婷婷精品国产亚洲av在线| 国产色爽女视频免费观看| 69人妻影院| 亚洲专区中文字幕在线| 日韩欧美国产在线观看| 一个人免费在线观看电影| 国产高潮美女av| 欧美国产日韩亚洲一区| 深爱激情五月婷婷| 免费在线观看影片大全网站| 内地一区二区视频在线| 乱人视频在线观看| 身体一侧抽搐| 日韩欧美在线二视频| 国产精品精品国产色婷婷| 日本黄大片高清| 国产精品久久电影中文字幕| 国产精品自产拍在线观看55亚洲| 国产极品精品免费视频能看的| 欧美日韩亚洲国产一区二区在线观看| 亚洲无线在线观看| 精品国产亚洲在线| 免费看日本二区| 久久久久久国产a免费观看| 欧美区成人在线视频| 成年女人毛片免费观看观看9| 精品久久国产蜜桃| 亚洲人成网站在线播放欧美日韩| 蜜桃亚洲精品一区二区三区| 亚洲国产欧美人成| 成人特级黄色片久久久久久久| 亚洲av不卡在线观看| 亚洲av电影在线进入| 精品国产亚洲在线| 精品久久久久久成人av| 欧美性猛交╳xxx乱大交人| 我的老师免费观看完整版| 丰满人妻熟妇乱又伦精品不卡| 欧美一区二区亚洲| 国产精品综合久久久久久久免费| 级片在线观看| 国产高清有码在线观看视频| 亚洲av成人不卡在线观看播放网| 亚洲av成人精品一区久久| 久久亚洲真实| 亚洲欧美日韩无卡精品| 麻豆一二三区av精品| 毛片女人毛片|