基于SPAC系統(tǒng)機(jī)理模型的火龍果灌溉管理研究

張 麗，張廷強(qiáng)

（廣西水利電力勘測(cè)設(shè)計(jì)研究院有限責(zé)任公司，南寧 530023）

0 引言

火龍果（Hylocereus undulates Britt）是一種熱帶、亞熱帶水果，又稱紅龍果、青龍果、仙蜜果等，屬仙人掌科（Cactaceae）量天尺屬（Hylocereus undatus）和蛇鞭柱屬（Seleniereus Mejalantous）植物[1，2]?；瘕埞a(chǎn)于中、南美洲等地區(qū)，二十世紀(jì)八九十年代最先傳入我國(guó)臺(tái)灣省，在20 世紀(jì)90 年代末開始引入福建等地進(jìn)行栽培。經(jīng)品種馴化后，受地理和氣候條件等的限制，目前只在廣西、廣東、貴州、云南、海南、臺(tái)灣等南方地區(qū)有較大規(guī)模種植，北方地區(qū)在增溫設(shè)施中也有所栽培[2，3]。火龍果喜疏松、富含有機(jī)質(zhì)、排水良好、pH 值范圍為6.0～7.5 的土壤環(huán)境，最適空氣溫度為25～35℃?；瘕埞粌H清甜爽口，且營(yíng)養(yǎng)價(jià)值高，同時(shí)壽命長(zhǎng)、繁殖快、見效快，具有較好的觀賞及經(jīng)濟(jì)價(jià)值，受到種植者和消費(fèi)者的喜愛。

火龍果雖然具有一定的耐旱能力，但是為了正常的生長(zhǎng)發(fā)育，仍需要充足的水分供應(yīng)及適宜的排水措施，特別是在開花結(jié)果期間[4]，否則易產(chǎn)生裂果[5]。目前火龍果種植過(guò)程中灌溉管理方式粗放，對(duì)火龍果需水規(guī)律針對(duì)性不強(qiáng)，甚至出現(xiàn)多余的灌水導(dǎo)致火龍果根部氧氣受限。本研究旨在通過(guò)模擬火龍果生長(zhǎng)過(guò)程，討論不同灌溉情形下火龍果的水分脅迫情況，并結(jié)合實(shí)際管理?xiàng)l件，制定合理的灌溉管理方案。通過(guò)本次研究，一方面可以合理高效地利用水資源，另一方面可以提高種植園灌溉技術(shù)水平和火龍果產(chǎn)量。

1 研究區(qū)及數(shù)據(jù)

1.1 研究區(qū)介紹

研究區(qū)位于廣西壯族自治區(qū)南寧市隆安縣丁當(dāng)鎮(zhèn)金福農(nóng)業(yè)火龍果基地，金福農(nóng)業(yè)火龍果基地面積達(dá)5000多畝，是廣西連片種植面積最大的紅心火龍果種植、種苗培育基地之一。項(xiàng)目區(qū)所在隆安縣年平均氣溫為21.8℃，年平均降水量為1 301.0 mm，年平均日照時(shí)數(shù)為1 528.4 h，年平均蒸發(fā)量1 632.8 mm，年平均相對(duì)濕度為80%。在研究區(qū)內(nèi)選擇代表性觀測(cè)田塊，布設(shè)土壤水分觀測(cè)設(shè)備和氣象觀測(cè)站（見圖1）。代表性觀測(cè)田塊在2020年4月1日至2020年10月15日（火龍果生長(zhǎng)耗水量大的時(shí)段）期間未進(jìn)行灌溉，以制造干旱條件，分析火龍果的水分脅迫規(guī)律。

圖1 研究區(qū)及土壤水分觀測(cè)點(diǎn)位置圖

1.2 數(shù)據(jù)

1.2.1 氣象數(shù)據(jù)

利用研究區(qū)布設(shè)的氣象觀測(cè)站，監(jiān)測(cè)太陽(yáng)輻射、空氣溫度（離地高度2 m）、降雨、濕度、風(fēng)速等氣象數(shù)據(jù)。氣象數(shù)據(jù)主要用于驅(qū)動(dòng)SPAC系統(tǒng)機(jī)理模型，圖2 展示了降雨隨時(shí)間的變化情況。由圖2 可知，在2020 年7 月初至2020 年9 月中旬，有兩個(gè)時(shí)段降雨量少，火龍果生長(zhǎng)可能受到脅迫。

圖2 研究時(shí)段內(nèi)降雨數(shù)據(jù)

1.2.2 土壤水分?jǐn)?shù)據(jù)

如圖1所示，在代表性田塊內(nèi)安裝有10個(gè)土壤水分監(jiān)測(cè)點(diǎn)，每個(gè)點(diǎn)監(jiān)測(cè)5、15 cm兩個(gè)深度（火龍果根系深度約為30 cm）內(nèi)的土壤含水量，為了準(zhǔn)確反映代表性田塊內(nèi)的土壤含水量情況，對(duì)10個(gè)設(shè)備監(jiān)測(cè)的結(jié)果進(jìn)行了平均，減小觀測(cè)數(shù)據(jù)的不確定性。

2 SPAC系統(tǒng)機(jī)理模型

采用土壤-水分-大氣-作物（SWAP）模型模擬土壤-植被-大氣連續(xù)體（SPAC）系統(tǒng)內(nèi)的土壤水分變化和火龍果根系吸水及水分脅迫過(guò)程。下面就土壤水分運(yùn)動(dòng)及水分脅迫模擬和基于SWAP模型的灌溉策略進(jìn)行介紹。

2.1 土壤水分運(yùn)動(dòng)

采用基于達(dá)西定律和質(zhì)量守恒定律推導(dǎo)的垂向一維Richards 方程描述農(nóng)田飽和-非飽和土壤水分運(yùn)動(dòng)：

式中：K（h）為水力傳導(dǎo)度，cm/d；h為壓力水頭，cm；z為垂向位置（以向上為正方向），cm；C（h）為土壤容水度，cm-1；t為時(shí)間，d；Sr（h）為植被根系吸水速率，cm3/d。

土壤本構(gòu)關(guān)系采用Mulem-van Genuchten 模型描述：

2.2 植被水分脅迫

作物生長(zhǎng)需要合適的土壤水分狀態(tài)，若土壤過(guò)濕即土壤含水量過(guò)高將導(dǎo)致作物根部缺氧而無(wú)法吸收土壤中的水分導(dǎo)致水分脅迫。若土壤過(guò)干及土壤含水量過(guò)低，作物根部將因無(wú)水可以吸收而出現(xiàn)水分脅迫。Feddes et al.（1978）根據(jù)作物的這一特征建立了基于土壤負(fù)壓水頭的根系水分脅迫模型，水分脅迫指數(shù)與土壤負(fù)壓的函數(shù)關(guān)系見圖3。圖3中水分脅迫指數(shù)越大代表作物受到的水分脅迫越弱，當(dāng)指數(shù)等于1 時(shí)，作物未遭受脅迫，當(dāng)指數(shù)等于0時(shí)，作物遭受最大水分脅迫，無(wú)生物量合成。

圖3 Feddes et al.（1978）水分脅迫模型

2.3 灌溉策略

SWAP 模型中共有6 種灌溉策略，其啟動(dòng)灌溉的標(biāo)準(zhǔn)分別為：①水分脅迫下限（式5）；②可用水量在總可用水量比值下限（式6）；③實(shí)際土壤含水量占飽和土壤含水量的比值下限（式7）；④絕對(duì)土壤含水量下限，指根系區(qū)的土壤含水量下限（式8）；⑤土壤含水量與田間持水量最大差值（式9）；⑥指定灌溉最小時(shí)間間隔和灌溉啟動(dòng)所對(duì)應(yīng)的最低土壤含水量與田間持水量的差值（式10）。在這6種灌溉策略中，灌溉水量的上限一般為田間持水量，也可以指定單次的灌溉水量。由圖3 可以看出，如果發(fā)生水分脅迫，作物生長(zhǎng)將受到影響，因此在灌溉時(shí)應(yīng)尋求水分脅迫與灌溉水量之間的平衡，做到水分脅迫較小且灌溉水量和灌溉時(shí)間適宜，不能一味地保證水分脅迫較小而無(wú)限增加灌溉水量和灌水頻率，否則不僅造成水資源的浪費(fèi)而且管理困難。本研究采用第一種灌溉策略并考慮實(shí)際管理?xiàng)l件設(shè)計(jì)合理的灌溉方式，結(jié)合第①種策略和圖3 可以看出，通過(guò)水分脅迫可以反映蒸騰量，進(jìn)而可表征產(chǎn)量的損失。SWAP 模型中6 種灌溉策略，其啟動(dòng)灌溉的標(biāo)準(zhǔn)分別為：

式中：Pa，t，P0，t分別為實(shí)際蒸騰量和潛在蒸發(fā)量，mm；WSl為水分脅迫閾值；Rθa，t，Rθw，Rθf(wàn)分別為根系區(qū)實(shí)際含水量、凋萎系數(shù)對(duì)應(yīng)的根系區(qū)土壤含水量和田間持水率對(duì)應(yīng)的根系區(qū)土壤含水量，mm；Rw為可用水比值閾值；Rθs為根系區(qū)土壤飽和含水量，mm；Rs為飽和比閾值；Rθ0根系區(qū)土壤含水量閾值，mm；Rf可用水余量閾值，mm；n灌溉次數(shù)順序序列號(hào)，t0最小間隔時(shí)間，d。

3 灌溉策略分析

3.1 模型參數(shù)校正

為了能準(zhǔn)確模擬土壤水分運(yùn)動(dòng)和火龍果水分脅迫信息，首先基于田間的根系觀測(cè)數(shù)據(jù)構(gòu)建根系模型（根系深度30 cm，須根系，根系分布函數(shù)滿足三角形分布），然后利用根系區(qū)土壤含水量觀測(cè)數(shù)據(jù)對(duì)SWAP 模型參數(shù)（式4、圖3）進(jìn)行校準(zhǔn)，優(yōu)化目標(biāo)是在有觀測(cè)數(shù)據(jù)時(shí)刻，模擬值與觀測(cè)值之間的均方根誤差最小。另外，利用室內(nèi)實(shí)驗(yàn)測(cè)定了火龍果試驗(yàn)區(qū)土壤的田間持水量數(shù)據(jù)，為制定合理的灌溉制度提供基礎(chǔ)數(shù)據(jù)。圖4展示了參數(shù)校準(zhǔn)后試驗(yàn)區(qū)觀測(cè)數(shù)據(jù)與模擬數(shù)據(jù)變化圖，由圖4可知，校準(zhǔn)后的模型能良好地反映土壤含水量的變化情況。

圖4 模型土壤含水量模擬情況

3.2 自然降雨情境下火龍果水分脅迫分析

圖5展示了自然降雨情況下降雨和水分脅迫指數(shù)之間的關(guān)系。由圖5可知，在7月中旬之前，降雨較為充沛，因此作物未受到明顯干旱影響。但是7月中旬至9月中旬時(shí)間段內(nèi)，雖然有降雨，但是研究區(qū)的季風(fēng)性氣候?qū)е陆涤暝跁r(shí)間上不均勻，部分時(shí)段無(wú)降雨，干旱較為嚴(yán)重，從圖5 可以看出，此時(shí)段內(nèi)若無(wú)灌溉，火龍果將遭受明顯的水分脅迫，其生長(zhǎng)將受到嚴(yán)重影響。

圖5 自然降雨情況下火龍果水分脅迫情況

3.3 不同水分脅迫下限閾值的灌溉對(duì)比分析

由圖5可知，在7月中旬至9月中旬時(shí)段內(nèi)應(yīng)進(jìn)行灌溉以保證作物良好生長(zhǎng)，采用第①種灌溉策略即設(shè)定水分脅迫指數(shù)下限來(lái)控制灌溉，下限分別為0.95、0.90、0.85、0.80、0.75、0.70、0.65、0.60、0.55、0.50，灌水上限均為田間持水量，不同水分脅迫指數(shù)下限閾值情況下水分脅迫指數(shù)變化情況見圖6。不同水分脅迫指數(shù)下限閾值和水分脅迫指數(shù)均值及灌溉水量的關(guān)系見圖7。由圖6和圖7可知，考慮到研究區(qū)的降雨情況，并不是水分脅迫指數(shù)下限越高越好，越高可能會(huì)因?yàn)楣喔葧r(shí)間后的降雨導(dǎo)致作物根系因缺氧而產(chǎn)生水分脅迫進(jìn)而影響到作物生長(zhǎng)。由圖7 可知，將水分脅迫指數(shù)下限設(shè)定在0.85～0.90之間最為合適，可同時(shí)兼顧較高的水分脅迫指數(shù)和較低的灌溉水量。

圖6 不同水分脅迫指數(shù)下限閾值情況下水分脅迫指數(shù)變化情況

圖7 不同水分脅迫指數(shù)下限閾值和水分脅迫指數(shù)均值及灌溉水量關(guān)系圖

3.4 考慮實(shí)際管理?xiàng)l件的灌溉對(duì)比分析

前面結(jié)果僅考慮了氣象并未考慮實(shí)際管理情況，在實(shí)際管理中，因研究區(qū)火龍果灌溉采用輪灌方式，且單次灌水一般在30 mm 以內(nèi)，兩次臨近灌溉時(shí)間的間隔不低于5 d?？紤]到實(shí)際管理情況，模擬了不同水分脅迫指數(shù)下限閾值情況下的水分脅迫指數(shù)變化情況（見圖8），不同水分脅迫指數(shù)下限閾值和水分脅迫指數(shù)均值及灌溉水量的關(guān)系見圖9。由圖8 和圖9 可知，考慮實(shí)際情況時(shí)，最優(yōu)的水分脅迫指數(shù)下限為0.95～0.96，可以實(shí)現(xiàn)較小的灌溉水量和較大的水分脅迫指數(shù)均值。相比于不考慮實(shí)際管理情況的結(jié)果而言，雖然灌溉次數(shù)（6 次>3次）增加了，但是灌溉水量（180 mm<191 mm）和水分脅迫指數(shù)均值（0.97>0.96）得到了更優(yōu)的結(jié)果，原因是研究區(qū)的季風(fēng)性氣候?qū)е陆涤晔构喔葲Q策制定變得更為復(fù)雜，從另一方面也可看出，基于現(xiàn)有的管理方式，結(jié)合建立的模型可以實(shí)現(xiàn)研究區(qū)火龍果的優(yōu)化灌溉管理。

圖8 實(shí)際管理?xiàng)l件下不同水分脅迫指數(shù)下限閾值情況下水分脅迫指數(shù)變化情況

圖9 實(shí)際管理?xiàng)l件下不同水分脅迫指數(shù)下限閾值和水分脅迫指數(shù)均值及灌溉水量關(guān)系圖

4 結(jié)論

火龍果作為重要的經(jīng)濟(jì)作物，合理的灌溉管理對(duì)于有效利用水資源和提高火龍果產(chǎn)量具有重要意義。本研究基于SPAC 系統(tǒng)機(jī)理模型，模擬了不同灌溉情形下火龍果的水分脅迫情況，并得到了對(duì)應(yīng)的灌溉時(shí)間和灌溉水量。結(jié)果顯示：無(wú)灌溉僅依賴自然降雨情形下，火龍果在蒸發(fā)顯著的夏秋季節(jié)遭受干旱脅迫明顯，說(shuō)明灌溉對(duì)于保證火龍果正常生長(zhǎng)和高產(chǎn)具有重要作用；結(jié)合實(shí)際管理?xiàng)l件，不同水分脅迫指數(shù)下限閾值的模擬結(jié)果表明單次灌水量在30 mm以內(nèi)，且臨近灌溉時(shí)間間隔不低于5 d時(shí)，最優(yōu)的水分脅迫指數(shù)下限閾值為0.95，灌水量在180 mm 左右。基于現(xiàn)有的管理方式，利用SAPC 系統(tǒng)機(jī)理模型可以實(shí)現(xiàn)火龍果的優(yōu)化灌溉管理。此外，從模擬的水分脅迫結(jié)果可以看出，除干旱外，在降雨集中且水量較大的時(shí)刻，火龍果根系因缺氧而受到脅迫，對(duì)產(chǎn)量也產(chǎn)生了一定的負(fù)面影響，應(yīng)注意做好排水措施。SPAC系統(tǒng)機(jī)理模型可實(shí)時(shí)反映脅迫情況，指導(dǎo)農(nóng)田水分（灌溉與排水）管理。