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

    土壤水氮遷移轉(zhuǎn)化模型研究進展

    2016-04-07 15:08:14張金萍
    水利科學(xué)與寒區(qū)工程 2016年10期
    關(guān)鍵詞:土壤水溶質(zhì)硝態(tài)

    張金萍,許 靈,倪 偉

    (1.水利部松遼水利委員會規(guī)研中心,吉林 長春 130021;2. 松遼水利水電開發(fā)有限責(zé)任公司,吉林 長春 130021)

    ?

    土壤水氮遷移轉(zhuǎn)化模型研究進展

    張金萍1,許 靈2,倪 偉1

    (1.水利部松遼水利委員會規(guī)研中心,吉林 長春 130021;2. 松遼水利水電開發(fā)有限責(zé)任公司,吉林 長春 130021)

    水資源短缺和農(nóng)業(yè)面源污染問題已成為農(nóng)業(yè)可持續(xù)發(fā)展的關(guān)鍵制約因素,因此提高農(nóng)田水氮的利用效率、減少氮素淋失及其對環(huán)境污染的影響具有十分重要的意義。隨著土壤溶質(zhì)運移理論研究的發(fā)展, 國內(nèi)外學(xué)者對土壤水氮運移、轉(zhuǎn)化機理的研究不斷深入。文章綜述了國內(nèi)外水氮運移、轉(zhuǎn)化理論研究及模擬模型研究的現(xiàn)狀, 并討論了該領(lǐng)域未來的研究趨勢和有待解決的問題。

    運移轉(zhuǎn)化機理;模擬模型;研究現(xiàn)狀

    肥料對糧食作物的增產(chǎn)、穩(wěn)產(chǎn)至關(guān)重要,生育期適時灌水和科學(xué)施肥是提高糧食產(chǎn)量的關(guān)鍵因素。氮素營養(yǎng)既是保障作物增產(chǎn)穩(wěn)產(chǎn)的關(guān)鍵營養(yǎng)元素之一[1],又是相對缺乏的關(guān)鍵元素[2]。土壤中氮素的含量和科學(xué)配比直接關(guān)系農(nóng)作物的健康生長[3]。為增加旱地土壤的氮素供給,氮素肥料在我國農(nóng)業(yè)生產(chǎn)實踐中被廣泛施用,年投入量高達2000萬t以上,居世界首位。2005年,國內(nèi)氮素化肥施用量4500萬t以上,約為1983年的2倍,而相應(yīng)的糧食產(chǎn)量為48 000萬t,僅是1983年的1.2倍。氮素肥料利用率較低,不僅造成農(nóng)業(yè)投入與產(chǎn)出嚴重失衡的浪費現(xiàn)象,更增加了農(nóng)業(yè)面源污染的風(fēng)險。已有研究[4]報道:國內(nèi)每年約有1000萬t氮素肥料轉(zhuǎn)化為污染物,造成農(nóng)業(yè)面源污染負荷,給生態(tài)環(huán)境和農(nóng)業(yè)生產(chǎn)環(huán)境帶來嚴重影響。農(nóng)田中氮素養(yǎng)分損失的途徑主要包括:徑流沖刷損失、氨的揮發(fā)、反硝化脫氮、銨的固定、硝態(tài)氮的淋失損失等。其中,硝態(tài)氮的淋失損失量最大[5],土壤硝態(tài)氮殘留在隨水流遷移到下層土壤進而污染淺層地下水。我國飲用水衛(wèi)生標準明確了水中的硝酸鹽含量的安全標準,每升水中不應(yīng)超過20 mg,雖高于發(fā)達國家標準[6],但部分地區(qū)飲水的水質(zhì)指標仍不達標。硝酸鹽對水體的污染具有隱蔽性,呈逐年加劇的態(tài)勢,治理難度大、代價高昂[7],已引起了社會的廣泛關(guān)注。因此,對于氮素污染問題的研究,已引起科學(xué)界廣泛重視[8],土壤中氮素的遷移轉(zhuǎn)化規(guī)律也已成為國內(nèi)外環(huán)境科學(xué)、生態(tài)科學(xué)和土壤科學(xué)研究的熱點問題[9-10]。

    水是農(nóng)業(yè)生產(chǎn)的命脈,在用水總量中,農(nóng)業(yè)用水所占比重最大。以華北平原為例,該地區(qū)一直是中國最大的糧食生產(chǎn)基地之一,谷物糧食產(chǎn)量占全國總量的21.6%[11],其中冬小麥的產(chǎn)量占全國總產(chǎn)量的71%左右[12],冬小麥全生育期耗水量為450 mm左右,而小麥生長季的降雨不足,一半不超過180 mm,小麥灌溉用水約占該地區(qū)農(nóng)業(yè)用水總量的70%[13],而其中64%的灌溉水來自地下[14-16],由于多年大量開采地下水進行灌溉,導(dǎo)致地下水位每年以1.0~1.5 m的速度下降,地下水位已從20世紀50年代的5 m 降到了現(xiàn)在的30 m以下[15],由此產(chǎn)生的環(huán)境矛盾也日趨尖銳。

    綜上,粗放的灌溉管理及不科學(xué)的施肥方式已成為生態(tài)環(huán)境風(fēng)險的主要因素,制定合理的灌溉制度和科學(xué)的施肥方案、提高灌溉水及氮素利用效率,減少不必要的耗散和損失,對平衡糧食增產(chǎn)與水資源及環(huán)境污染之間的矛盾有極為重要的現(xiàn)實意義。如何確定灌水量、施肥量與氮素淋失量、作物產(chǎn)量等之間的相關(guān)性具有重要意義,然而大田實驗成本高且周期長,不同農(nóng)業(yè)管理措施下的研究過程既費時又耗資,對環(huán)境的影響的評價往往難于達到預(yù)期的效果。因此,研究土壤中氮素的遷移轉(zhuǎn)化規(guī)律及模型,利用模型對不同管理措施進行模擬評估具有十分重要的理論與現(xiàn)實意義。

    1 土壤水氮遷移轉(zhuǎn)化理論研究現(xiàn)狀

    相比國內(nèi),國外對土壤溶質(zhì)運移方面問題的研究開展較早,始于20世紀50年代。由Lapidus開創(chuàng)的對流—彌散方程(CDE)模型,揭開了溶質(zhì)運移研究的序幕。Scheidegg將Lapidus 的方程擴展,進行三維均質(zhì)土壤、穩(wěn)態(tài)流條件下的模擬研究,優(yōu)化了算法,同時考慮了溶質(zhì)運移時水的動力彌散作用。Rizfai在Schidegger研究成果的基礎(chǔ)上,考慮了溶質(zhì)運移時的分子擴散作用,并引入了彌散度的概念,來表征土壤特性對溶質(zhì)運移的影響,使溶質(zhì)運移理論的研究更加深入。Nielson和Biggaz從理論上推導(dǎo)建立了對流—彌散方程,并根據(jù)實驗結(jié)果,對Lapidus、Schidegger和Nielson的模型進行了比較分析,因此實驗和理論進一步表明對流—彌散方程能較好地描述非反映性物質(zhì)在多孔介質(zhì)中的遷移規(guī)律[17]。在溶質(zhì)運移研究歷史上,開創(chuàng)了應(yīng)用數(shù)學(xué)模型來說明和解釋溶質(zhì)運移過程的局面。Biggerh和Nieison提出了考慮土壤的空間變異性問題,開始了對土壤性質(zhì)空間變異性的研究[18]。美國加州大學(xué)Jury提出了“黑箱模型”[19]。國外許多學(xué)者為了量化溶質(zhì)運移規(guī)律,提出了溶質(zhì)運移的定量計算模型,早期的研究多采用確定性機理模型即對流—彌散方程(CDE)來模擬土壤中的溶質(zhì)運移,隨著隨機理論的發(fā)展,隨機過程方程例如TEM模型,傳遞函數(shù)修正模型,此外還有蒙特卡洛方法[20]、譜分析法[21]、矩分析法[22]等開始應(yīng)用于溶質(zhì)運移數(shù)量特征的模擬研究。

    近年來,國內(nèi)土壤物理學(xué)者及農(nóng)學(xué)者開展了一些室內(nèi)、室外的溶質(zhì)運移試驗研究,并結(jié)合模型模擬系統(tǒng)研究了氮素的遷移轉(zhuǎn)化規(guī)律[23]。武曉峰等[24]對冬小麥田間根層氮素遷移轉(zhuǎn)化規(guī)律進行了研究;馮紹元等[25]系統(tǒng)地綜述了氮素遷移轉(zhuǎn)化規(guī)律及其損失對水環(huán)境的影響等。羅陽等[26]對大田中土壤剖面上水分及三氮含量等進行了長期監(jiān)測,并結(jié)合室內(nèi)的土柱模擬實驗,應(yīng)用LEACHN模型進行了模擬計算,研究結(jié)果表明土壤根區(qū)底部的氮素的淋溶損失量與降雨量和灌溉量、施肥量呈正向相關(guān)。張瑜芳等[27]提出了排水條件下氮素遷移轉(zhuǎn)化和流失的簡化計算方法,系統(tǒng)地研究了排水條件下氮肥遷移轉(zhuǎn)化規(guī)律。

    隨著節(jié)水灌溉的推廣應(yīng)用,國內(nèi)學(xué)者對噴灌和滴灌條件下氮素遷移轉(zhuǎn)化進行了系統(tǒng)研究。王康[28]建立了節(jié)水條件下土壤氮素損失和環(huán)境評價概念型模型,用以研究節(jié)水條件下氮素對生態(tài)環(huán)的境影響。馮紹元等[29]對噴灌條件硝態(tài)氮含量與施肥量的相關(guān)性進行了系統(tǒng)研究。李久生等[30]對滴灌點源施肥條件下硝態(tài)氮在距滴頭附件的分布累積規(guī)律進行了系統(tǒng)的研究。

    2 土壤水氮遷移轉(zhuǎn)化模型研究進展

    應(yīng)用數(shù)學(xué)方法模擬農(nóng)田中土壤水氮遷移轉(zhuǎn)化規(guī)律,是定量化研究的必然趨勢。根據(jù)土壤水氮遷移轉(zhuǎn)化模型的求解形式不同,Addisott和Wagenet[31]將其劃分為確定性模型和隨機模型。其中確定性模型又可劃分為確定性機理模型和確定性函數(shù)模型[32]。確定性機理模型即對流—彌散方程(CDE)幾十年來應(yīng)用最為普遍,該模型是根據(jù)達西定律和連續(xù)原理建立的,可以較好地揭示溶質(zhì)在均質(zhì)多孔介質(zhì)中的運移機理及時間、空間對溶質(zhì)運移的影響。確定性函數(shù)模型如活塞流、半解析模型、層狀模型等對溶質(zhì)運移和水流運動做了概化,大多采用物質(zhì)平衡的方法估算土壤水氮的行為[33]。

    隨機模型近20 a來逐漸發(fā)展成熟,并受到重視和廣泛應(yīng)用。其主要特點是輸入變量及輸出結(jié)果均以概率分布的形式表示,而不是單一的數(shù)值。與國外相比,國內(nèi)在隨機模型應(yīng)用研究處于起步較晚,應(yīng)用隨機方法建模和大致可分為兩種方式,即單獨考慮時間變量和空間變量,或同時考慮兩者的隨機模型。目前隨機模型主要有蒙特卡洛模擬、拉格朗日方法、歐拉法和隨機連續(xù)模型,在土壤水力特性參數(shù)及溶質(zhì)運移參數(shù)、水流通量、溶質(zhì)運移速度及彌散系數(shù)等方面應(yīng)用較多[33]。

    現(xiàn)就目前國際上有影響的描述土壤水氮遷移轉(zhuǎn)化的幾個模型做下簡要介紹。

    (1)HYDRUS-1D

    美國鹽土實驗室開發(fā)的HYDRUS-1D,包括3個模塊:水分模塊、溶質(zhì)運移模塊及溫度模塊,廣泛用于模擬飽和—非飽和介質(zhì)中水、熱和溶質(zhì)的運移。水分模塊采用Richards方程計算土壤中水分含量,溶質(zhì)運移模塊和溫度模塊采用對流彌散方程,分別模擬溶質(zhì)運移轉(zhuǎn)化過程和溫度的對流傳導(dǎo)過程。

    局限性:沒有考慮溫度、水分對氮素反應(yīng)速率的影響,且無作物生長模塊,只是對作物吸水做了簡單的處理,因此不能模擬土壤水氮遷移轉(zhuǎn)化對作物生長的影響。

    (2)LEACHM(Leaching Estimation And Chemistry Model)

    美國康奈爾大學(xué)研究開發(fā)的LEACHM模型包括4個模塊:LEACHW模塊、LEACHP模塊、LEACHN模塊和LEACHC模塊,被廣泛應(yīng)用于土壤中水分和溶質(zhì)的運動模擬。其中LEACHW模塊用于模擬水分運動,LEACHP模塊用于模擬殺蟲劑的遷移轉(zhuǎn)化,LEACHN用于模擬土壤中碳、氮和磷的遷移轉(zhuǎn)化,LEACHC用于模擬土壤中鹽分運移。

    局限性: ① 土壤分層厚度必須相等; ②不能模擬地表徑流; ③不能模擬作物對土壤和環(huán)境影響的響應(yīng) ;④不能模擬作物產(chǎn)量。

    (3)RZWQM(Root Zone Water Quality Model)

    由美國農(nóng)業(yè)部大平原系統(tǒng)研究單位在NTRM模型的基礎(chǔ)上修正完善的RZWQM模型包括6個模塊:作物生長模塊、物理運移模塊、化學(xué)運移模塊、養(yǎng)分循環(huán)模塊、殺蟲劑模塊和管理模塊。該模型廣泛應(yīng)用于模擬不同農(nóng)業(yè)管理措施及其對溶質(zhì)運移狀況和作物產(chǎn)量的影響。

    局限性:模型將植物殘茬分為快、慢兩種降解形式,將土壤有機物分成快、中、慢3種速率庫,各個庫之間轉(zhuǎn)換復(fù)雜,參數(shù)較多,模型對這些降解庫的動態(tài)非常敏感,因此模型不易率定。

    (4)NLEAP模型

    主要用于評價由于農(nóng)業(yè)措施的變化而產(chǎn)生的硝酸鹽淋失量對地下水的潛在威脅,用于田間尺度。主要動力學(xué)過程包括:水和硝態(tài)氮的運移過程、作物吸收過程、氮的礦化、硝化與反硝化過程、氨揮發(fā)過程、作物腐殖質(zhì)與有機肥和其他有機物的礦化和固定過程等。

    局限性:①模擬時間步長必須大于等于一天 ;②不能模擬地表徑流產(chǎn)生的氮損失; ③沒有考慮表層有機肥的氨揮發(fā)和反硝化影響。

    此外,還有模擬氮素在土壤中遷移轉(zhuǎn)化機理的DRAINMOD模型、CREAMS 模型,模擬土壤-作物系統(tǒng)中土壤水與氮素動力過程的DAISY模型,用于評估土壤侵蝕對土地生產(chǎn)力、植物競爭、作物產(chǎn)量和水質(zhì)影響的NTRM模型等。

    中國在該領(lǐng)域的研究雖然沒有達到獨立建立機理模型的階段,但是結(jié)合大田實驗和室內(nèi)實驗監(jiān)測,近年來也在模型研究方面取得了較好的研究成果。武曉峰和謝森傳[34]建立了模擬土壤-作物系統(tǒng)中不同形態(tài)氮素遷移轉(zhuǎn)化的數(shù)學(xué)模型,以冬小麥田間實驗觀測為基礎(chǔ),模擬研究了不同狀態(tài)氮素的移轉(zhuǎn)化規(guī)律以及土壤吸附、作物吸收等多種影響因素,利用溶質(zhì)擴散-對流方程模擬了冬小麥生長期間田間水分、氨氮、硝態(tài)氮含量及其分布變化;張瑜芳等[35]建立了稻田中氨態(tài)氮遷移轉(zhuǎn)化數(shù)值模擬,并與土柱試驗觀測資料進行了對比分析。黃元仿等[36]應(yīng)用對流—彌散方程,建立了田間土壤氮素運移的數(shù)值模型,將土壤氮素、水和熱運移綜合考慮,并將土壤氮素轉(zhuǎn)化和作物根系吸收也加入在模型的源匯項中;馮紹元等[37]建立了土壤—水—植物系統(tǒng)中氮素遷移轉(zhuǎn)化與綜合吸收數(shù)學(xué)模型,模型考慮了氮素的礦化、硝化、反硝化和作物吸收等過程,并對轉(zhuǎn)化參數(shù)進行了濕度和溫度的修正,礦化、硝化和反硝化過程采用一級動力學(xué)方程描述,作物吸氮采用米氏競爭方程計算,并將模型運用于玉米、西紅柿和萵苣三種作物生長過程的模擬;沈榮開等[38]建立了一個田間一維飽和—非飽和土壤中NH4+-N和NO3--N在排水條件下的遷移轉(zhuǎn)化模型,模型中綜合考慮了對流、彌散、吸附、降解和硝化作用,對氮素在非飽和土壤中的運移進行了探討,并對模型進行了驗證。

    在模型應(yīng)用研究方面,張思聰?shù)萚39]以長期田間觀測資料為基礎(chǔ),應(yīng)用LEACHM模型,對根區(qū)以下硝態(tài)氮的遷移轉(zhuǎn)化規(guī)律和滲漏損失的關(guān)鍵影響因素進行了研究;曹巧紅等[40]應(yīng)用Hydrus-1D模型模擬了非飽和介質(zhì)中一維水、熱和溶質(zhì)遷移過程,獲得了取得了較好的結(jié)果;任理和馬軍花[41]建立了硝態(tài)氮在土壤中淋失的傳遞函數(shù)模型,并與實測資料對比驗證,結(jié)果表明所建模型比常規(guī)溶質(zhì)的傳遞函數(shù)模型具有更高精度。

    3 結(jié)論與討論

    3.1 區(qū)域水氮運移模型與地理信息系統(tǒng)的結(jié)合

    地理信息系統(tǒng)(GIS)具有強大的數(shù)據(jù)輸入、存儲、檢索、運算、顯示功能,應(yīng)用其地理模型分析方法,可以提供多種空間和動態(tài)的地理信息,為地理研究和管理決策提供技術(shù)支持?,F(xiàn)在GIS被廣泛地應(yīng)用于繪制如作物產(chǎn)量分布圖、種植結(jié)構(gòu)圖、硝態(tài)氮淋洗分布圖等。因此,在研究區(qū)域土壤水氮遷移轉(zhuǎn)化過程時,可以將基于土壤水氮遷移轉(zhuǎn)化機理的點模型與GIS結(jié)合,可以更加直觀地展示模擬運算結(jié)果,便于理解和解釋。模型與GIS結(jié)合方式分為松結(jié)合、緊結(jié)合和完全結(jié)合三種類型[42-43],目前國內(nèi)外大多是采用緊結(jié)合的方式研究土壤水氮遷移轉(zhuǎn)化過程,并向應(yīng)用完全結(jié)合方向轉(zhuǎn)移。

    3.2 多種模型結(jié)合應(yīng)用,加強多因素交互作用研究

    目前的土壤水氮運移、轉(zhuǎn)化模型大多考慮了主要的氮素動力學(xué)過程,包括氮肥施用過程、硝化與反硝化過程、氨揮發(fā)過程、作物吸收過程和硝態(tài)氮淋失過程等。但是,不同的模型各有其優(yōu)缺點,由于模型的局限性或者說模型的適用條件所限,沒有哪一種模型能把氮素遷移轉(zhuǎn)化所有過程都模擬得很好,每一種模型都側(cè)重于一個或者幾個過程,因此根據(jù)研究問題的實際情況,將幾種模型結(jié)合進行研究,可以提供解決問題的思路。

    3.3 存在問題

    (2)非飽和帶氮素運移機理與模型研究將注意力集中于硝態(tài)氮淋失、氮素利用效率與作物產(chǎn)量的關(guān)系方面,氮素對地下水的污染模型研究則集中于含水層。將飽和帶—非飽和帶作為整個系統(tǒng),研究氮素轉(zhuǎn)化與運移規(guī)律并分析對環(huán)境影響的研究相對較少。

    (3)土壤氮素是一個相互聯(lián)系、多因素交互作用的統(tǒng)一體,須要采用系統(tǒng)分析的方法進行研究,而國內(nèi)外對于土壤氮素遷移轉(zhuǎn)化規(guī)律的研究多集中在單個過程,忽視了氮素轉(zhuǎn)化各過程的交互作用,對整個系統(tǒng)進行整體研究時也存在關(guān)注單一因素的情況,今后應(yīng)加強對多因素交互作用的研究。

    [1] 趙平, 孫谷疇, 彭少麟. 植物氮素營養(yǎng)的生理生態(tài)學(xué)研究[J]. 生態(tài)科學(xué), 1998(2):37-42.

    [2] 李秧秧,邵明安.小麥根系對水分和氮肥的生理生態(tài)反應(yīng)[J].植物營養(yǎng)與肥料學(xué)報,2000,6(4):383-388.

    [3] 王夏暉, 劉軍, 王益權(quán). 不同施肥方式下土壤氮素的運移特征研究[J]. 土壤通報, 2002, 33(3):202-206.

    [4] 趙其國. 發(fā)展與創(chuàng)新現(xiàn)代土壤科學(xué)[J]. 土壤學(xué)報, 2003, 40(3): 321-327.

    [5] HUTSON ┼, A. CASS. A retentivity function for use in soil-water simulation models[J]. European Journal of Soil Science, 1987, 38(1):105-113.

    [6] Almasri M N, Kaluarachchi J J. Modeling nitrate contamination of groundwater in agricultural watersheds[J]. Journal of Hydrology, 2007, 343(3-4):211-229.

    [7] 袁新民,王周瓊.硝態(tài)氮的淋洗及其影響因素[J].干旱區(qū)研究,2000,17(4):46-52.

    [8] Bergstr?m L, Johnsson H, Torstensson G. Simulation of soil nitrogen dynamics using the SOILN model[J]. Nutrient Cycling in Agroecosystems, 1991, 27(27):181-188.

    [9] Knighton R E.Simulation of solute transport using a CTMP[J].Water Reasor.Res.,1987,28(10):1917-1925.

    [10] Chang C, Entz T J. Nitrate leaching losses under repeated cattle feed lot manure applications in Southern Alberta[J].Journal of Environmental Quality, 1996,25(1):145-153.

    [11] Yu Q, Saseendran S A, Ma L, et al. Modeling a wheat-maize double cropping system in China using two plant growth mdules in RZWQM[J]. Agricultural Systems, 2006, 89(2-3):457-477.

    [12] Sun H Y, Liu C M, Zhang X Y, et al. Effects of irrigation on water balance, yield and WUE of winter wheat in the North China Plain[J]. Agricultural Water Management, 2006, 85(1-2): 211-218.

    [13] 居輝,蘭霞,李建民,等.不同灌溉制度下冬小麥產(chǎn)量效應(yīng)與耗水特征研究[J]. 中國農(nóng)業(yè)大學(xué)學(xué)報, 2000, 5(5): 23-29.

    [14] 王淑芬,張喜英,裴冬. 不同供水條件對冬小麥根系分布產(chǎn)量及水分利用效率的影響[J]. 農(nóng)業(yè)工程學(xué)報, 2006, 22(2):27-32.

    [15] Li J M, Inanaga S, Li Z H, et al. Optimizing irrigation scheduling for winter wheat in the North Chian Plain[J]. Agricultural Water Management, 2005, 76(1): 8-23.

    [16] Hu C S, Saseendran S A, Green T R, et al. Evaluating nitrogen and water management in a double-cropping system using[J]. Vadose Zone Journal, 2006, 5(1): 493-505.

    [17] 魏新平,王文焰,王全九,等.溶質(zhì)運移理論的研究現(xiàn)狀和發(fā)展趨勢[J].灌溉排水,1998 (4): 58-63.

    [18] Nielsen D R, Van Genuchten M T, Biggar J W. Flow and Solute Transport Processes in the Unsaturated Zone [J]. Water Resources,1986,22(9):89-108.

    [19] Jury W A, Sposito G, White R E. A transfer function model of solute transport through soil [J]. Water Resources Research, 1986, 22(2):243-247.

    [20] Persaud N, Giraldez J V, Chang A C. Monte-Carlo Simulation of Noninteracting Solute Transport in a Spatially Heterogeneous Soil1[J]. Soil Science Society of America Journal, 1985, 49(3):562-568.

    [21] 黃冠華,葉自桐,楊金忠.一維非飽和溶質(zhì)隨機運移模型的譜分析[J].水利學(xué)報,1995(11):1-7.

    [22] Aris R. On the dispersion of linear kine matic waves[J]. Proc.R.Soc..Lon-don.Ser.A.,1958,245(1241):268-277.

    [23] 徐曉榮,李恒輝,陳良.利用15N 研究氮肥對土壤及植物內(nèi)硝酸鹽的影響[J].核農(nóng)學(xué)報,2000,14(5):301-304.

    [24] 武曉峰,謝森傳.冬小麥田間根層中氮素遷移轉(zhuǎn)化規(guī)律研究[J].灌溉排水,1996 (4):10-15.

    [25] 馮紹元,王鳳新,黃冠華.噴灌條件下花生水肥耦合效應(yīng)的田間試驗研究[J].農(nóng)業(yè)工程學(xué)報,1998(4):98-102.

    [26] 羅陽, 張思聰. 灌溉施肥條件下田間氮素在土壤中遷移情況的研究[J]. 水資源保護, 2000(4):7-11.

    [27] 張瑜芳,張蔚榛,沈榮開,等.排水農(nóng)田中氮素轉(zhuǎn)化運移和流失[M]. 武漢:中國地質(zhì)大學(xué)出版社,1997.

    [28] 王康,沈榮開.節(jié)水條件下土壤氮素的環(huán)境影響效應(yīng)研究[J].水科學(xué)進展,2003,14(4):437-441.

    [29] 馮紹元,王鳳新,黃冠華.噴灌條件下花生水肥耦合效應(yīng)的田間試驗研究[J].農(nóng)業(yè)工程學(xué)報,1998(4):98-102.

    [30] Li Jiusheng,Zhang Jianjun,Ren Li. Nitrogen Distributions in Soil Under Fertigation From a Point Source[J]. Transactions of the CSAE,2002,18(5): 61-66.

    [31] Addiscott T M, Wagenet R J. Concepts of solute leaching in soils: a review of modelling approaches[J]. European Journal of Soil Science, 1985, 36(3):411-424.

    [32] Wolf J, Keulen H V. Modeling long-term crop response to fertilizer and soil nitrogen[J].Plant & Soil, 1989, 120(1):11-22.

    [33] 高如泰. 黃淮海平原農(nóng)田土壤水氮行為模擬與管理分析[D].北京: 中國農(nóng)業(yè)大學(xué), 2005.

    [34] 武曉峰,謝森傳. 冬小麥田間根層中氮素遷移轉(zhuǎn)化規(guī)律研究[J].灌溉排水, 1996, 15(4):10-15.

    [35] 張瑜芳,劉培斌.不同滲漏強度條件下淹水稻田中氨態(tài)氮轉(zhuǎn)化和運移的研究[J].水利學(xué)報,1994(6):9-19.

    [36] 黃元仿,李韻珠,陸錦文.田間條件下土壤氮素運移的模擬模型[J].水利學(xué)報,1996(6):9-14.

    [37] 馮紹元,張瑜芳,沈榮開.非飽和土壤中氮素運移與轉(zhuǎn)化試驗及其數(shù)值模擬[J].水利學(xué)報,1996(8):8-15.

    [38] 沈榮開, 楊詩秀, 謝森傳. 土壤水動力學(xué)[M]. 北京: 清華大學(xué)出版社,1988.

    [39] 張思聰,呂賢弼,黃永剛.灌溉施肥條件下氮素在土壤中遷移轉(zhuǎn)化的研究[J].水利水電技術(shù),1999,30(5):6-8.

    [40] 曹巧紅,龔元石.應(yīng)用Hydrus-1D模型模擬分析冬小麥農(nóng)田水分氮素運移特征[J].植物營養(yǎng)與肥料學(xué)報,2003,9(2):139-145.

    [41] 任理,馬軍花. 考慮土壤中硝態(tài)氮轉(zhuǎn)化作用的傳遞函數(shù)模型[J]. 水利學(xué)報. 2001( 5): 38-44.

    [42] Tim U S. Coupling Vadose Zone Models with GIS: Emerging Trends and Potential Bottlenecks[J]. Journal of Environmental Quality, 1996, 25(3):535-544.

    [43] Fedra K. Distributed models and embedded GIS: integration strategies and case studies, GIS and environmental modelling: progress and research issues[M]. USA: GIS World, Inc., Fort Collins, CO, 1996:413-417.

    Progress of research on soil water and nitrogen transport and transformation

    ZHANG Jinping1,XU Ling2,NI Wei1

    (1.SongliaoWaterResourcesCommit,MWR.Changchun130021,China;2.SongliaowaterresourcesandhydropowerdevelopmentCo.,Ltd,Changchun130021,China)

    The shortage of water resource and agricultural non-point source pollution are the limiting factors for sustainable agriculture in the North China Plain (NCP). It is necessary to improve the water use efficiency (WUE), nitrogen use efficiency (NUE) and reduce the nitrogen leaching from the root zone and its effect on environmental pollution. With the development of soil solute transport theory, the study on the water and nitrogen transport and transformation mechanisms is developing in home and abroad. This paper summarizes the current status of research in area of nitrogen transport and transformation mechanisms and simulating models.The development tendency of research on this area was briefly discussed.

    transport and transformation;Simulating Model;Present Status

    張金萍(1980-),女,高級工程師,主要從事水利規(guī)劃方面的工作。

    S153;S154.4

    A

    2096-0506(2016)10-0013-05

    猜你喜歡
    土壤水溶質(zhì)硝態(tài)
    有關(guān)溶質(zhì)質(zhì)量分數(shù)的計算
    滴水成“冰”
    溶質(zhì)質(zhì)量分數(shù)考點突破
    改進的PSO-RBF模型在土壤水入滲參數(shù)非線性預(yù)測中的應(yīng)用研究
    錦州市土壤水動態(tài)過程及影響因素
    灌水定額對土壤水鹽分布及作物產(chǎn)量的影響
    “溶質(zhì)的質(zhì)量分數(shù)”計算歸類解析
    低C/N比污水反硝化過程中亞硝態(tài)氮累積特性研究
    硝態(tài)氮供應(yīng)下植物側(cè)根生長發(fā)育的響應(yīng)機制
    控釋復(fù)合肥對冷季型草坪氨揮發(fā)和硝態(tài)氮淋洗的影響
    免费女性裸体啪啪无遮挡网站| 日韩有码中文字幕| 午夜免费鲁丝| 午夜老司机福利片| 黑丝袜美女国产一区| 乱人伦中国视频| 国产精品久久电影中文字幕 | 亚洲国产中文字幕在线视频| 成人手机av| 久久精品国产亚洲av香蕉五月 | 欧美+亚洲+日韩+国产| 免费在线观看视频国产中文字幕亚洲| 制服诱惑二区| 亚洲五月色婷婷综合| 国产精品免费视频内射| 国产熟女午夜一区二区三区| tocl精华| 国产成人精品在线电影| 久久久久精品国产欧美久久久| 又大又爽又粗| 久久久国产成人免费| 99热网站在线观看| 成年女人毛片免费观看观看9 | 久久久欧美国产精品| 国产成人精品在线电影| 日本欧美视频一区| 精品熟女少妇八av免费久了| 亚洲男人天堂网一区| 国产成人欧美在线观看 | 黄色视频在线播放观看不卡| 中文字幕人妻熟女乱码| 国产精品免费一区二区三区在线 | 亚洲精品美女久久av网站| 老熟妇仑乱视频hdxx| 亚洲色图综合在线观看| 男女之事视频高清在线观看| 涩涩av久久男人的天堂| 两个人看的免费小视频| 亚洲伊人久久精品综合| 欧美日韩成人在线一区二区| 亚洲伊人久久精品综合| √禁漫天堂资源中文www| 国产精品欧美亚洲77777| 国产高清国产精品国产三级| 国产精品一区二区在线不卡| 国产精品美女特级片免费视频播放器 | 国产av精品麻豆| 亚洲欧美日韩另类电影网站| 日本黄色视频三级网站网址 | 韩国精品一区二区三区| av线在线观看网站| 欧美 日韩 精品 国产| 黄色 视频免费看| 少妇精品久久久久久久| 欧美亚洲日本最大视频资源| 欧美乱码精品一区二区三区| 欧美精品一区二区免费开放| 精品国产乱码久久久久久男人| 国产av一区二区精品久久| 久久热在线av| 国产免费av片在线观看野外av| 免费在线观看日本一区| 国产国语露脸激情在线看| 人妻久久中文字幕网| 精品少妇久久久久久888优播| 精品久久久精品久久久| 麻豆乱淫一区二区| 黑人巨大精品欧美一区二区蜜桃| 亚洲国产看品久久| 久久亚洲精品不卡| 中文字幕高清在线视频| 桃花免费在线播放| 夜夜夜夜夜久久久久| 伦理电影免费视频| 啦啦啦视频在线资源免费观看| av电影中文网址| 国产亚洲一区二区精品| 亚洲色图av天堂| 搡老岳熟女国产| 99精国产麻豆久久婷婷| 国产激情久久老熟女| 老司机亚洲免费影院| 免费黄频网站在线观看国产| 精品免费久久久久久久清纯 | 黑人欧美特级aaaaaa片| 手机成人av网站| 在线观看一区二区三区激情| 国产成人精品久久二区二区91| 一级片'在线观看视频| 国产老妇伦熟女老妇高清| 免费少妇av软件| 精品福利永久在线观看| 久久人妻av系列| 亚洲精品国产区一区二| 少妇粗大呻吟视频| 三上悠亚av全集在线观看| 一本久久精品| 99国产极品粉嫩在线观看| 亚洲视频免费观看视频| 国产日韩欧美在线精品| 99国产综合亚洲精品| 99久久国产精品久久久| 一区二区三区乱码不卡18| 99久久人妻综合| 视频区欧美日本亚洲| 亚洲人成电影免费在线| 夜夜骑夜夜射夜夜干| 欧美在线黄色| 久久久久国内视频| 免费一级毛片在线播放高清视频 | 这个男人来自地球电影免费观看| 日韩大片免费观看网站| 日韩中文字幕欧美一区二区| 啦啦啦 在线观看视频| 亚洲熟女精品中文字幕| 大片免费播放器 马上看| 午夜免费成人在线视频| 国产高清激情床上av| 亚洲人成电影观看| 91九色精品人成在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲精华国产精华精| 黄色片一级片一级黄色片| 美女福利国产在线| 丁香六月天网| 久久中文字幕一级| 操出白浆在线播放| 成人精品一区二区免费| 97在线人人人人妻| 91麻豆精品激情在线观看国产 | 国产精品影院久久| 免费在线观看影片大全网站| www.999成人在线观看| 男女床上黄色一级片免费看| 亚洲一码二码三码区别大吗| 美女高潮到喷水免费观看| 国产熟女午夜一区二区三区| 亚洲欧美日韩高清在线视频 | 天堂俺去俺来也www色官网| 女人被躁到高潮嗷嗷叫费观| 欧美激情 高清一区二区三区| 亚洲国产欧美一区二区综合| 国产97色在线日韩免费| 香蕉久久夜色| 丰满饥渴人妻一区二区三| 男女无遮挡免费网站观看| 99在线人妻在线中文字幕 | 国产不卡一卡二| 中文字幕人妻丝袜一区二区| 午夜福利在线免费观看网站| 美女高潮喷水抽搐中文字幕| 99riav亚洲国产免费| 正在播放国产对白刺激| 国产成人av激情在线播放| 一本久久精品| 每晚都被弄得嗷嗷叫到高潮| 他把我摸到了高潮在线观看 | 五月天丁香电影| 亚洲精品在线美女| 在线永久观看黄色视频| 美女高潮喷水抽搐中文字幕| 99热网站在线观看| 久久国产精品影院| 午夜福利在线免费观看网站| 久久热在线av| 国产真人三级小视频在线观看| 中国美女看黄片| 日日摸夜夜添夜夜添小说| 久久精品国产a三级三级三级| 国产免费福利视频在线观看| 黑人欧美特级aaaaaa片| 天天影视国产精品| 亚洲人成电影免费在线| 变态另类成人亚洲欧美熟女 | 黄色视频,在线免费观看| 一区二区av电影网| 777久久人妻少妇嫩草av网站| 黄片小视频在线播放| 黄色片一级片一级黄色片| 国产老妇伦熟女老妇高清| 国产精品二区激情视频| e午夜精品久久久久久久| 久久国产精品影院| 免费观看人在逋| 精品欧美一区二区三区在线| 久久人人97超碰香蕉20202| 亚洲精品中文字幕一二三四区 | 91成年电影在线观看| 一本—道久久a久久精品蜜桃钙片| 搡老熟女国产l中国老女人| 欧美精品高潮呻吟av久久| 汤姆久久久久久久影院中文字幕| 热99国产精品久久久久久7| 日韩欧美一区视频在线观看| 久久久精品国产亚洲av高清涩受| 国产精品久久久久久精品电影小说| 国产不卡av网站在线观看| 国产成人av教育| 日韩免费高清中文字幕av| 精品乱码久久久久久99久播| 欧美另类亚洲清纯唯美| 老熟女久久久| 后天国语完整版免费观看| 免费在线观看日本一区| 久久久欧美国产精品| 日韩一卡2卡3卡4卡2021年| 一边摸一边做爽爽视频免费| 大陆偷拍与自拍| 国产男女内射视频| 777久久人妻少妇嫩草av网站| 美女高潮喷水抽搐中文字幕| 国产成人一区二区三区免费视频网站| 成人特级黄色片久久久久久久 | 国产成人免费无遮挡视频| 在线 av 中文字幕| 一区福利在线观看| 亚洲久久久国产精品| 免费在线观看完整版高清| 男女免费视频国产| 久久久久国产一级毛片高清牌| 热99国产精品久久久久久7| 无遮挡黄片免费观看| 男女之事视频高清在线观看| 日本wwww免费看| 欧美日韩国产mv在线观看视频| 女同久久另类99精品国产91| tube8黄色片| a级片在线免费高清观看视频| 最近最新中文字幕大全免费视频| 日本一区二区免费在线视频| 777米奇影视久久| 亚洲欧美一区二区三区久久| 午夜视频精品福利| 中文字幕人妻丝袜一区二区| 69精品国产乱码久久久| 精品国内亚洲2022精品成人 | 黑人猛操日本美女一级片| 侵犯人妻中文字幕一二三四区| 成年人免费黄色播放视频| 日本撒尿小便嘘嘘汇集6| 亚洲五月色婷婷综合| 久久人妻福利社区极品人妻图片| 国精品久久久久久国模美| 成年动漫av网址| 99热国产这里只有精品6| 手机成人av网站| 精品久久久久久久毛片微露脸| 丰满人妻熟妇乱又伦精品不卡| 新久久久久国产一级毛片| 1024视频免费在线观看| 91九色精品人成在线观看| 中文字幕人妻丝袜一区二区| 国产高清国产精品国产三级| 午夜福利乱码中文字幕| 亚洲免费av在线视频| 亚洲中文字幕日韩| 免费观看av网站的网址| 久久ye,这里只有精品| 欧美精品高潮呻吟av久久| 99热国产这里只有精品6| 又紧又爽又黄一区二区| 夜夜爽天天搞| 国产老妇伦熟女老妇高清| 亚洲av日韩精品久久久久久密| 午夜久久久在线观看| 蜜桃国产av成人99| 亚洲av国产av综合av卡| 国产精品一区二区精品视频观看| 中文字幕人妻丝袜制服| 国产av国产精品国产| 精品视频人人做人人爽| 两人在一起打扑克的视频| 欧美亚洲日本最大视频资源| 黄片大片在线免费观看| 水蜜桃什么品种好| 91av网站免费观看| 真人做人爱边吃奶动态| 黄色视频不卡| 50天的宝宝边吃奶边哭怎么回事| 免费在线观看完整版高清| 免费久久久久久久精品成人欧美视频| 亚洲国产毛片av蜜桃av| 十分钟在线观看高清视频www| 国产日韩欧美视频二区| 视频在线观看一区二区三区| 午夜日韩欧美国产| 亚洲久久久国产精品| www.熟女人妻精品国产| 91成人精品电影| 久9热在线精品视频| 久久国产精品影院| 亚洲精华国产精华精| 免费在线观看视频国产中文字幕亚洲| 欧美日韩亚洲综合一区二区三区_| 午夜成年电影在线免费观看| 免费黄频网站在线观看国产| 国产97色在线日韩免费| 涩涩av久久男人的天堂| 精品国产一区二区久久| 黑人巨大精品欧美一区二区蜜桃| 一进一出抽搐动态| 亚洲精品中文字幕在线视频| 日本av免费视频播放| 侵犯人妻中文字幕一二三四区| 久久99一区二区三区| 91九色精品人成在线观看| 亚洲精品粉嫩美女一区| 日日爽夜夜爽网站| 午夜福利,免费看| 欧美日韩亚洲国产一区二区在线观看 | 麻豆成人av在线观看| 99国产综合亚洲精品| 日日摸夜夜添夜夜添小说| 男女床上黄色一级片免费看| 国产精品一区二区免费欧美| 人人妻人人爽人人添夜夜欢视频| 国产亚洲欧美精品永久| 在线观看一区二区三区激情| av又黄又爽大尺度在线免费看| 亚洲精品在线观看二区| 国产精品美女特级片免费视频播放器 | 夜夜爽天天搞| 热99国产精品久久久久久7| 丝袜喷水一区| 色综合婷婷激情| 50天的宝宝边吃奶边哭怎么回事| 国产精品影院久久| 午夜福利在线观看吧| 亚洲中文日韩欧美视频| 国产精品麻豆人妻色哟哟久久| 国产精品 国内视频| 日日摸夜夜添夜夜添小说| av在线播放免费不卡| 欧美精品一区二区免费开放| 欧美国产精品va在线观看不卡| 国产精品免费一区二区三区在线 | 欧美午夜高清在线| 天天添夜夜摸| 国产精品成人在线| 精品国产国语对白av| 人妻 亚洲 视频| 亚洲五月婷婷丁香| 欧美日韩国产mv在线观看视频| 欧美午夜高清在线| 一区二区av电影网| 九色亚洲精品在线播放| 国产一区二区三区在线臀色熟女 | 欧美日本中文国产一区发布| 成人国产av品久久久| 免费在线观看完整版高清| 黑人猛操日本美女一级片| 亚洲精品久久成人aⅴ小说| 色综合欧美亚洲国产小说| 亚洲国产毛片av蜜桃av| 一级片'在线观看视频| 日韩欧美国产一区二区入口| 欧美亚洲 丝袜 人妻 在线| 黑人欧美特级aaaaaa片| 国产真人三级小视频在线观看| 美女福利国产在线| 男人舔女人的私密视频| 精品免费久久久久久久清纯 | 国产成人啪精品午夜网站| 老熟妇仑乱视频hdxx| av网站免费在线观看视频| 91老司机精品| 99国产极品粉嫩在线观看| 久久这里只有精品19| 黑人巨大精品欧美一区二区蜜桃| 操美女的视频在线观看| 久久久国产欧美日韩av| 最新在线观看一区二区三区| 丝袜喷水一区| 熟女少妇亚洲综合色aaa.| 午夜福利乱码中文字幕| 制服人妻中文乱码| 日本黄色日本黄色录像| 国产av国产精品国产| 1024香蕉在线观看| 91麻豆av在线| 精品国产一区二区三区久久久樱花| 亚洲欧美色中文字幕在线| 国产一区二区三区在线臀色熟女 | 新久久久久国产一级毛片| 黄片大片在线免费观看| 国产精品国产高清国产av | 黄频高清免费视频| 久久热在线av| 国产精品国产av在线观看| 日本wwww免费看| 黄色 视频免费看| 在线观看免费午夜福利视频| 成人av一区二区三区在线看| 精品久久久久久电影网| 天天影视国产精品| 免费av中文字幕在线| h视频一区二区三区| svipshipincom国产片| 热99久久久久精品小说推荐| 久久人人爽av亚洲精品天堂| 亚洲精品粉嫩美女一区| 丰满少妇做爰视频| 日韩欧美免费精品| 欧美精品亚洲一区二区| 母亲3免费完整高清在线观看| 女性被躁到高潮视频| av超薄肉色丝袜交足视频| 国产视频一区二区在线看| 久久性视频一级片| 日韩一区二区三区影片| 最黄视频免费看| 大型黄色视频在线免费观看| 国产精品亚洲av一区麻豆| 91成人精品电影| 777米奇影视久久| 久热这里只有精品99| 成人影院久久| 成年人免费黄色播放视频| 中文字幕高清在线视频| 精品国产乱码久久久久久男人| 757午夜福利合集在线观看| 母亲3免费完整高清在线观看| 亚洲欧美日韩另类电影网站| 久9热在线精品视频| 女人精品久久久久毛片| 黄色成人免费大全| 18禁美女被吸乳视频| 国产免费视频播放在线视频| 女人高潮潮喷娇喘18禁视频| av欧美777| 三级毛片av免费| 美女视频免费永久观看网站| 国产伦人伦偷精品视频| 在线观看免费高清a一片| 免费看a级黄色片| 中国美女看黄片| 手机成人av网站| 国产亚洲精品久久久久5区| 国产成人免费无遮挡视频| 亚洲男人天堂网一区| 国产精品欧美亚洲77777| 国产真人三级小视频在线观看| 久久久久久亚洲精品国产蜜桃av| 天堂俺去俺来也www色官网| a级毛片在线看网站| 色婷婷久久久亚洲欧美| 久久中文字幕人妻熟女| 免费在线观看日本一区| 搡老岳熟女国产| 一本—道久久a久久精品蜜桃钙片| 亚洲伊人色综图| 最新在线观看一区二区三区| 欧美变态另类bdsm刘玥| 国产av国产精品国产| 99精国产麻豆久久婷婷| 国产精品98久久久久久宅男小说| 这个男人来自地球电影免费观看| 狠狠精品人妻久久久久久综合| 久久这里只有精品19| 啦啦啦免费观看视频1| 欧美大码av| 国产在线免费精品| 99国产精品一区二区三区| 免费看十八禁软件| 精品国产亚洲在线| 色婷婷av一区二区三区视频| 人妻 亚洲 视频| 亚洲精品久久午夜乱码| 国产av精品麻豆| 777久久人妻少妇嫩草av网站| 如日韩欧美国产精品一区二区三区| av线在线观看网站| 亚洲欧洲日产国产| 99re6热这里在线精品视频| 国产免费视频播放在线视频| 最近最新中文字幕大全免费视频| 天堂俺去俺来也www色官网| 国产成人一区二区三区免费视频网站| 男女无遮挡免费网站观看| 91老司机精品| 久久人妻福利社区极品人妻图片| 亚洲精品成人av观看孕妇| 久久久国产成人免费| 丁香六月天网| 国产精品一区二区在线不卡| 国产成人精品在线电影| av免费在线观看网站| 午夜精品久久久久久毛片777| 999久久久国产精品视频| 久久久久国产一级毛片高清牌| 亚洲视频免费观看视频| 亚洲成人手机| 中亚洲国语对白在线视频| 国产一区二区三区视频了| 无人区码免费观看不卡 | 国产精品秋霞免费鲁丝片| 国产淫语在线视频| av国产精品久久久久影院| 黄频高清免费视频| 精品国产超薄肉色丝袜足j| 日本一区二区免费在线视频| 成年动漫av网址| av国产精品久久久久影院| 黄频高清免费视频| 另类精品久久| 一区二区av电影网| 考比视频在线观看| 丝袜美足系列| 视频在线观看一区二区三区| 老熟女久久久| 黄色丝袜av网址大全| 亚洲av日韩精品久久久久久密| 美女国产高潮福利片在线看| 在线观看免费日韩欧美大片| 国产欧美日韩一区二区三区在线| 亚洲一区二区三区欧美精品| 老司机在亚洲福利影院| 国产欧美亚洲国产| 91av网站免费观看| 色播在线永久视频| 久久99一区二区三区| 岛国在线观看网站| 高清黄色对白视频在线免费看| 欧美精品啪啪一区二区三区| 国产成人精品无人区| 黄色视频在线播放观看不卡| 一边摸一边抽搐一进一小说 | 国产精品1区2区在线观看. | 午夜福利在线免费观看网站| 一本综合久久免费| 男女免费视频国产| 欧美精品亚洲一区二区| 欧美日韩视频精品一区| 国产精品影院久久| 日韩大片免费观看网站| 黑人巨大精品欧美一区二区蜜桃| 可以免费在线观看a视频的电影网站| 国产又爽黄色视频| 亚洲综合色网址| av视频免费观看在线观看| 成人国语在线视频| 美女国产高潮福利片在线看| 超碰97精品在线观看| 欧美激情高清一区二区三区| 亚洲专区字幕在线| 国产亚洲精品久久久久5区| 99九九在线精品视频| 国产高清激情床上av| 午夜久久久在线观看| 久久狼人影院| 国产97色在线日韩免费| 日韩 欧美 亚洲 中文字幕| www.熟女人妻精品国产| 欧美黑人欧美精品刺激| 我的亚洲天堂| 日韩熟女老妇一区二区性免费视频| 精品国产国语对白av| 免费观看av网站的网址| av电影中文网址| 人成视频在线观看免费观看| 日韩中文字幕欧美一区二区| 久久香蕉激情| 国产男靠女视频免费网站| 欧美精品一区二区大全| 欧美黄色淫秽网站| 狠狠精品人妻久久久久久综合| 三级毛片av免费| av天堂在线播放| 日韩大片免费观看网站| avwww免费| 大码成人一级视频| 日韩欧美一区二区三区在线观看 | 欧美日韩成人在线一区二区| 每晚都被弄得嗷嗷叫到高潮| e午夜精品久久久久久久| 成年动漫av网址| 久久人人爽av亚洲精品天堂| 久久久国产精品麻豆| 超碰成人久久| 一边摸一边抽搐一进一小说 | 国产成人系列免费观看| 一级毛片电影观看| 亚洲视频免费观看视频| 脱女人内裤的视频| 90打野战视频偷拍视频| 高清欧美精品videossex| 国产伦人伦偷精品视频| 777久久人妻少妇嫩草av网站| 欧美乱妇无乱码| 一级黄色大片毛片| 超色免费av| 国产亚洲av高清不卡| 欧美日韩视频精品一区| 国产日韩欧美亚洲二区| 黄色丝袜av网址大全| 99在线人妻在线中文字幕 | 又紧又爽又黄一区二区| 欧美日韩一级在线毛片| 女人高潮潮喷娇喘18禁视频| 国产亚洲欧美精品永久| 人妻一区二区av| xxxhd国产人妻xxx| 国内毛片毛片毛片毛片毛片| 中文字幕精品免费在线观看视频| 亚洲av国产av综合av卡| 亚洲av成人不卡在线观看播放网| 亚洲免费av在线视频| avwww免费| 精品久久久精品久久久| 69精品国产乱码久久久| 亚洲视频免费观看视频| 色婷婷久久久亚洲欧美| 日韩一卡2卡3卡4卡2021年| 国产在线视频一区二区| 亚洲五月色婷婷综合|