保水劑A對橡膠樹專用肥養(yǎng)分溶出的影響

鄭國亮等

摘 要 在橡膠樹專用肥中添加不同比例的保水劑，研究其對橡膠樹專用肥養(yǎng)分釋放的影響，并利用經(jīng)驗(yàn)方程對肥料養(yǎng)分釋放曲線進(jìn)行擬合，以期為開發(fā)新型橡膠專用肥提供依據(jù)。結(jié)果表明：緩釋效果隨著保水劑A用量的增加而增強(qiáng)，尋找到較好的模型進(jìn)行擬合研究，運(yùn)用雙常數(shù)方程和一級動力學(xué)方程對溶出速率曲線和累積溶出率曲線分別進(jìn)行擬合，相關(guān)系數(shù)R2分別為0.990 8～1.000 0和0.992 3～0.999 9，并得出保水劑A用量與緩釋效果的關(guān)系，提出以模型參數(shù)的變化規(guī)律、溶出抑制率的顯著性來確定緩/控釋材料的最佳用量，保水劑A的用量在5.0%～6.0%之間較為適宜。

關(guān)鍵詞 保水劑；橡膠樹專用肥；緩釋效果；模型擬合；溶出抑制率

中圖分類號 S794.1 文獻(xiàn)標(biāo)識碼 A

Abstract Water retention agent with different proportions added in the rubber fertilizer was employed to investigate its influence on the nutrient stripping of specialized rubber fertilizer. An empirical equation was modeled for the nutrient release curve in order to provide a basis for the development of new type rubber fertilizers. Results showed that the slow-release effect increased with more dosage of water retention agent A； Dissolution rate and cumulative dissolution rate curve were modeled by the double constant curve equation and first order kinetics equation， and the correlation coefficient R2 was 0.990 8-1.000 0， and 0.992 3-0.999 9， respectively， and the relationship between the dosage of water retention agent A and slow-release effect was determined. It was put forward that the optimal dosage of slow/controlled release materials was determined by the change rule of model parameters and the significant of dissolution inhibition rate and the amount of water retention agent was more appropriate between 5.0% and 6.0%.

Key words Water retention agent；Specialized rubber fertilizer；Slow-release effect；Model；Dissolution inhibition rate

doi 10.3969/j.issn.1000-2561.2014.09.030

我國橡膠種植園區(qū)主要位于熱帶雨林、季雨林地區(qū)，該地區(qū)雨熱同期，降雨量大且集中，5～10月為雨季，降水占全年的80%～85%，正處于橡膠樹一年中的生長旺盛期，也是化肥大量施用期，因而橡膠樹生長需要的肥料在施用后易遭淋失[1]。通過研發(fā)合適的緩/控釋肥料，減少肥料養(yǎng)分淋失，提高肥料利用率具有重要的意義[2-3]。緩/控釋肥料主要是通過包膜限制、基質(zhì)限制、化學(xué)結(jié)合、化學(xué)抑制以及物理與化學(xué)方法結(jié)合進(jìn)行限制和抑制起緩釋作用[4-8]。保水劑作為緩釋材料與肥料復(fù)合的工藝有物理吸附型、混合造粒、包膜型、化學(xué)結(jié)合型4種，其中物理吸附型和混合造粒是通過基質(zhì)限制起緩釋作用；包膜型和化學(xué)結(jié)合型分別是通過包膜限制和化學(xué)結(jié)合起緩釋作用[9-11]。保水劑緩/控釋肥料除了對肥料養(yǎng)分有緩釋作用外，還有保水、改土、可被生物降解的優(yōu)點(diǎn)[12-13]。物理吸附型、混合造粒復(fù)合的加工工藝比包膜型、化學(xué)結(jié)合的復(fù)合工藝簡便易行，再加工成本較低，而且混合造粒比物理吸附在運(yùn)輸和使用方便。但是非水溶性保水劑與肥料造粒的效果很不理想，只有水溶性保水劑與肥料有良好的造粒效果[9]。因此，筆者選擇1種適合的水溶性保水劑進(jìn)行試驗(yàn)，以不同的比值與橡膠樹專用肥混合，研究混合后的橡膠樹專用肥養(yǎng)分的靜水培養(yǎng)溶出規(guī)律，旨在為研制適宜在橡膠生產(chǎn)上推廣應(yīng)用的緩/控釋肥料提供初步參考。

1 材料與方法

1.1 材料

供試肥料：橡膠樹專用肥，用尿素（化學(xué)純，廣州化學(xué)試劑廠生產(chǎn)）、氯化鉀（化學(xué)純，廣州化學(xué)試劑廠生產(chǎn)）和過磷酸鈣（中化化肥有限公司海南分公司生產(chǎn)）按照橡膠樹專用肥的配比進(jìn)行配制，配比為N ∶ P2O5 ∶ K2O=15 ∶ 9 ∶ 6。

供試材料：保水劑A，一種粉末狀高分子可溶性保水材料，本身分子的結(jié)構(gòu)可使無機(jī)離子和養(yǎng)分暫時被吸水溶脹的保水劑分子包裹起來，也可通過范德華力、靜電引力、螯合等機(jī)制增加對養(yǎng)分的吸附作用。

1.2 方法

將保水劑與橡膠樹專用肥（5 g）均勻混合配制供試肥料，0.0%、0.5%、1.0%、2.0%、3.0%、4.0%、5.0%、6.0%，分別用CK、A1、A2、A3、A4、A5、A6、A7表示。試驗(yàn)設(shè)10個重復(fù)。每個處理的每個重復(fù)分別用3層醫(yī)用紗布進(jìn)行包裹，壓緊，用細(xì)尼龍繩扎口制成肥料包，并把肥料包分別放在裝有100 mL去離子水的燒杯里在室溫條件下進(jìn)行靜水培養(yǎng)。在培養(yǎng)兩個小時后，分別將肥料包取出，更換裝有100 mL去離子水的燒杯進(jìn)行靜水培養(yǎng)，接下來分別在浸泡第4、6、8、10、12、14、16 h后更換裝有100 mL去離子水的燒杯中進(jìn)行靜水培養(yǎng)，測定每一次燒杯里溶液中水溶性K和全N含量。K和N分別用火焰光度計和NPK流動分析儀進(jìn)行測定。

1.3 數(shù)據(jù)分析

1.3.1 計算公式 溶出率D[t1，t2]=（處理在區(qū)間[t1，t2]內(nèi)溶出養(yǎng)分的量/肥料中該養(yǎng)分的含量）×100；溶出速率DVt=D[（t-1），（t+1）]/2；累積溶出率SDt=D2+D4+……+Dt-2+Dt；溶出抑制率DIt=[SDt（對照）-SDt（處理）]×100/SDt（對照）；式中t、D、DVt、SDt和DIt的單位分別為h、%、%/h、%和%。

1.3.2 數(shù)據(jù)分析處理 利用Microsoft office Excel 2010軟件計算各處理溶出速率、累積溶出率、溶出相對抑制率和繪制圖表，利用JMP 9.0軟件的非線性回歸分析模塊進(jìn)行模型擬合和Tukey HSD測驗(yàn)對溶出抑制率進(jìn)行方差分析。

2 結(jié)果與分析

2.1 保水劑A用量對橡膠樹專用肥N素溶出特征的影響

2.1.1 對N素溶出速率的影響 保水劑A用量對橡膠樹專用肥N素溶出速率的影響，如圖1中A所示：在第1小時時，各處理的溶出速率由大到小的順序?yàn)镃K>A1>A2>A3>A4>A5>A6>A7，其中處理A1、A2和CK的相差不明顯。在第3小時時，各處理及空白對照的溶出速率由大到小的順序?yàn)锳6>A5>A4>A7>A3>A2>A1>CK，其中處理A3、A4、A5、A6、A7的速率明顯大于處理A1、A2和CK。在第5小時后，處理A1、A2和CK的溶出速率幾乎為0，而其它處理的溶出速率由大到小的順序?yàn)锳7>A6>A5>A4>A3，并且速率隨時間的增加而逐漸減小。

2.1.2 N素溶出速率曲線的擬合 以雙常數(shù)方程模型Ln（Vt）=a+bLn（t）對各處理的N素溶出速率曲線進(jìn)行擬合，結(jié)果如表1所示，參數(shù)b由大到小的順序?yàn)锳7>A6>A5>A4>A3>A2>A1>CK，緩釋性能參數(shù)b隨著保水劑A用量的增加而增大，說明對橡膠樹專用肥釋放抑制的效果隨著材料添加量的增加而增加。N素溶出速率曲線擬合模型的材料緩釋性能參數(shù)b與材料添加量x之間的數(shù)量關(guān)系用二次函數(shù)擬合的方程為b=-0.012 3x2+0.167 5x+0.216 5（相關(guān)系數(shù)R2=0.985 5）。該二次函數(shù)的頂點(diǎn)坐標(biāo)為（6.808 9、0.786 6），說明隨著保水劑A添加量的增加，對橡膠樹專用肥養(yǎng)分釋放的抑制效果的增加在逐漸變小，并預(yù)計在添加量約為6.8%時達(dá)到比較好的抑制效果，材料添加量繼續(xù)增加后效果增加將不明顯。

2.1.3 對N素累積溶出率的影響 保水劑A用量對橡膠樹專用肥N素溶出速率的影響，如圖1中B所示：處理A1、A2 的N素累積溶出率曲線與CK相差不明顯，而處理A3、A4、A5、A6、A7的N素累積溶出率曲線的斜率均明顯小于CK，隨著用量繼續(xù)增加，曲線的斜率逐漸變小，說明保水劑A對橡膠樹專用肥N素溶出的抑制效果隨著用量增加而增加。各處理及對照CK的累積溶出率達(dá)到80%的先后順序?yàn)镃K>A1>A2>A3>A4>A5>A6>A7，到達(dá)80%時的時間段分別為第2、2、2、4、5、7、8和10小時，并且處理A1、A2和CK的累積溶出率到達(dá)80%時的時間點(diǎn)極為相近，說明保水劑A的用量達(dá)到2.0%后才對橡膠樹專用肥N素的溶出有抑制效果，并且抑制效果隨著用量的增加而增加。

2.1.4 N素累積溶出率曲線的擬合 以一級動力學(xué)方程N(yùn)t=N0[1-exp（-kt）]對各處理的N素累積溶出率曲線進(jìn)行擬合，結(jié)果如表1所示，根據(jù)速率常數(shù)k的大小，各處理的順序?yàn)镃K>A1>A2>A3>A4>A5>A6>A7，說明肥料養(yǎng)分釋放速率隨著材料添加量的增加而減小，緩釋效果逐漸增強(qiáng)。N素累積溶出率曲線擬合模型的速率常數(shù)k與材料添加量x之間的數(shù)量關(guān)系用二次函數(shù)擬合的方程為k=0.042 6x2-0.455 4x+1.453 7（相關(guān)系數(shù)R2=0.987 3）。該二次函數(shù)的頂點(diǎn)坐標(biāo)為（5.345 1，0.236 6），說明隨著保水劑A添加量的增加，對橡膠樹專用肥養(yǎng)分中N素釋放速率抑制效果的增加量在逐漸變小，頂點(diǎn)顯示添加量在5.3%時達(dá)到較好的抑制效果，材料添加量繼續(xù)增加后效果增加不明顯。

2.1.5 對N素溶出抑制率的影響 保水劑A用量對橡膠樹專用肥N素養(yǎng)分的溶出抑制情況如表2所示：各處理的N素溶出抑制率隨著培養(yǎng)時間的延長而減小，且各處理對N素溶出抑制效果的大體順序?yàn)锳7>A6>A5>A4>A3>CK>A2>A1。處理A3、A4、A5、A6、A7的N素溶出抑制率分別在2～6、2～8、2～12、2～14、2～16 h時，和對照處理CK差異顯著；其中處理A3、A4、A5、A6、A7分別在2～4、2～8、2～12、2～12、2～16 h時，和對照處理CK差異極顯著。處理A7的N素溶出抑制率和其他處理均差異顯著。說明處理A7對N素溶出的抑制效果最好，處理A6次之。

2.2 保水劑A用量對橡膠樹專用肥K素溶出特征的影響

2.2.1 對K素溶出速率的影響 保水劑A用量對橡膠樹專用肥K素溶出速率的影響，如圖2中A所示：在第1小時時，各處理及空白對照的溶出速率由大到小的順序?yàn)镃K>A1>A2>A3>A4>A5>A6>A7，其中處理A1、A2和CK、處理A6和A7相差不明顯。在第3小時時，各處理及空白對照的溶出速率由大到小的順序?yàn)锳6>A4>A5>A7>A3>A2>A1>CK，其中處理A3、A4、A5、A6、A7的速率明顯大于處理A1、A2和CK，且處理A4、A5、A6、A7不明顯。在第5小時后，各處理及空白對照的溶出速率由大到小的順序?yàn)锳7>A6>A5>A4>A3>A2>A1>CK，并且速率隨時間的增加逐漸減小。

2.2.2 K素溶出速率曲線的擬合 以雙常數(shù)方程模型Ln（Vt）=a+bLn（t）對各處理的K素溶出速率曲線進(jìn)行擬合，結(jié)果如表3所示，按參數(shù)b由大到小的順序?yàn)锳7>A6>A5>A4>A3>A2>A1>CK，說明保水劑A對橡膠樹專用肥養(yǎng)分緩釋效果隨著材料添加量的增加而增加。K素溶出速率曲線擬合模型的材料性能參數(shù)b與材料添加量x的之間的數(shù)量關(guān)系用二次函數(shù)擬合的方程為，b=-0.008 1x2+0.113 9x+0.300 8（相關(guān)系數(shù)R2=0.966 0）。該二次函數(shù)的頂點(diǎn)坐標(biāo)為（7.030 9，0.701 2），說明隨著保水劑A添加量的增加，對肥料釋放的抑制效果的回報率在逐漸變小，并在添加量將在7.0%時達(dá)到較好的抑制效果，材料添加量繼續(xù)增加后效果增加將不明顯。

2.2.3 對K素累積溶出率的影響 保水劑A用量對橡膠樹專用肥K素累積溶出率的影響，如圖2中B所示：處理A1、A2的 K素累積溶出率曲線與CK相差不明顯，而處理A3、A4、A5、A6、A7的K素累積溶出率曲線的斜率均明顯小于CK處理，隨著用量繼續(xù)增加，曲線的斜率逐漸變小，說明保水劑A對橡膠樹專用肥K素溶出的抑制作用隨著用量增加而增加。各處理及對照CK的累積溶出率達(dá)到80%的先后順序?yàn)镃K>A1> A2>A3>A4>A5>A6>A7，到達(dá)80%時的時間段分別為第2、2、2、4、5、6、7和7小時，并且處理A1、A2和CK的累積溶出率到達(dá)80%時的時間點(diǎn)極為相近，處理A6和A7相近，說明保水劑A在添加量達(dá)到2%后才對橡膠樹專用肥K素溶出有抑制效果，并且效果隨著用量的增加而增加，當(dāng)用量增加到5%以后效果增加不明顯。

2.2.4 K素累積溶出率曲線的擬合 以一級動力學(xué)方程N(yùn)t=N0[1-exp（-kt）]對各處理的K素累積溶出率曲線進(jìn)行擬合，結(jié)果如表3所示，根據(jù)速率常數(shù)k的大小，各處理的順序?yàn)镃K>A1>A2>A3>A4>A5>A6>A7，說明肥料養(yǎng)分釋放速率隨著材料添加量的增加而減小。K素累積溶出率曲線擬合模型的速率常數(shù)k與材料添加量x之間的數(shù)量關(guān)系用二次函數(shù)擬合的方程為k=0.022 6x2-0.271 7x+1.152（相關(guān)系數(shù)R2=0.970 2）。該二次函數(shù)的頂點(diǎn)坐標(biāo)為（6.011 1，0.335 4），說明隨著保水劑A添加量的增加，對橡膠樹專用肥養(yǎng)分中K素釋放速率抑制效果的增加量在逐漸變小，頂點(diǎn)顯示添加量在6.0%時達(dá)到較好的抑制效果，材料添加量繼續(xù)增加后效果增加將不明顯。

2.2.5 對K素溶出抑制率的影響 保水劑A用量對橡膠樹專用肥K素養(yǎng)分的溶出抑制情況如表4所示：各處理的K素溶出抑制率隨著培養(yǎng)時間的延長而減小，且各處理對K素溶出的抑制效果大體順序?yàn)锳7>A6>A5>A4>A3>A2>A1>CK。處理A2、A3、A4、A5、A6、A7的K素溶出抑制率分別在8～16、2～16、2～16、2～16、2～16、2～16小時時，和對照處理CK差異顯著；其中處理A3、A4、A5、A6、A7分別在2～10、2～16、2～16、2～16、2～16小時時，和對照處理CK差異極顯著。處理A4、A5、A6的K素溶出抑制率分別在12～16、8～16、2～16小時時，和處理A7差異不顯著。說明處理A6、A7對K素溶出的抑制效果最好，處理A5次之。

3 討論與結(jié)論

3.1 保水劑A的緩釋效應(yīng)

隨著保水劑A用量的增加，對橡膠樹專用肥中的N、K養(yǎng)分的緩釋效果顯著增強(qiáng)，而且保水劑A用量在1.0%～2.0%間，處理間的緩釋效果有一個飛躍式的增強(qiáng)過程。這說明保水劑A主要通過吸漲對肥料進(jìn)行包裹，限制肥料養(yǎng)分的溶出，達(dá)到緩釋的目的：用量小于1.0%時，肥料不能被完全包裹，肥料與水溶液直接接觸導(dǎo)致養(yǎng)分迅速溶出，達(dá)不到緩釋效果；用量在2.0%以上時，保水劑A吸漲后能將肥料包裹完全，肥料不與水溶液直接接觸，溶出受到限制，釋放放緩。學(xué)者張富倉等[12]指出保水劑本身分子的結(jié)構(gòu)特性，使無機(jī)離子和養(yǎng)分可以暫時被吸水溶脹的保水劑分子包裹起來，也可通過范德華力、靜電引力、離子吸附、離子交換、螯合等機(jī)制增加對養(yǎng)分的吸附作用，抑制流失達(dá)到緩釋效果。說明保水劑膠體對肥料包裹起到一定的緩釋效果。

各處理N素和K素的溶出率隨著保水劑A用量增加而逐漸減小，說明保水劑A對肥料養(yǎng)分的緩釋效果，隨著用量的增加而增強(qiáng)；在增加同樣的用量情況下，保水劑A對N素緩釋效果的增幅比對K素的緩釋效果增幅大。用量低于1.0%時，保水劑A對K素的緩釋效果比N素的強(qiáng)，但是用量超過2.0%時，保水劑A對N素的緩釋效果比K素的強(qiáng)。說明材料用量超過2.0%后，保水劑A通過吸漲體積膨大將肥料完全包裹，此時主要抑制肥料養(yǎng)分溶出的可能是分子間的作用力，而不是離子吸附作用。因?yàn)樵谟兴那闆r下KCl中的K素是以K+的形式存在，分子比以酰胺態(tài)存在的N素的分子小，在膠體中遷移快，所以保水劑A對K素的緩釋效果比N素弱。但也可能是由于保水劑中的基團(tuán)與尿素分子鏈接減緩尿素的水解速率，減緩N素的溶出，隨著保水劑A用量的增加，對N素緩釋效果增強(qiáng)。王旭東等[14]也指出大分子中的羧基、酚羥基、甲氧基等基團(tuán)可以和尿素分子進(jìn)行鏈接，形成共聚體，從而減少尿素分子的移動性，增加其穩(wěn)定性，減緩尿素的水解速率。

3.2 保水劑A用量與其緩釋性能的關(guān)系

保水劑A對于K素的緩釋效果，隨著材料用量的增加，材料用量在0.0%～1.0%間緩釋效果增強(qiáng)速度緩慢；材料用量在1.0%～2.0%緩釋效果增強(qiáng)速度迅速；材料用量在2.0%～5.0%間緩釋效果增強(qiáng)速度逐漸減緩慢；材料用量大于5.0%時，緩釋效果不繼續(xù)增強(qiáng)。對N素的緩釋效果變化趨勢也與K素的相似，但是N素的緩釋效果潛力比K素的大。說明保水劑A對橡膠樹專用肥的緩釋效果隨著材料用量的增加呈“S型”的增加的趨勢。材料的用量與緩釋效果的關(guān)系趨勢與學(xué)者郝世雄[15]用聚乙烯醇水溶液制備緩/控釋包膜尿素的結(jié)果類似。

3.3 模型參數(shù)與緩釋性能的關(guān)系

試驗(yàn)中，一級反應(yīng)動力學(xué)方程模型Nt=N0[1-exp（-kt）]對養(yǎng)分累積溶出率擬合較好，擬合度在0.989 9～0.999 9間；雙常數(shù)方程模型Ln（Vt）=a+bLn（t）（為方程模型Vt=a·tb的衍生方程Vt=ea·tb的變形方程）對養(yǎng)分溶出速率擬合較好，擬合度在0.990 8～1.000 0間。兩個方程模型的參數(shù)k和參數(shù)b分別為溶出速率常數(shù)和材料性能參數(shù)，參數(shù)k的值越大養(yǎng)分釋放越快，參數(shù)b越大緩釋效果愈好[16-18]。試驗(yàn)中，參數(shù)k隨著保水劑A用量的增加逐漸減小，而參數(shù)b隨著用量增加逐漸增大，都說明隨著保水劑A用量的增加緩釋效果逐漸增強(qiáng)。其中參數(shù)b和k有極好的相關(guān)性：bN=0.959 5exp（-0.958kN），R2=0.999 4；bK=0.953 9exp（-0.968kK），R2=0.998 9。說明兩個模型的反映材料對養(yǎng)分的緩釋能力是相關(guān)的，兩個參數(shù)的隨著材料用量變化關(guān)系均能反映材料用量對提高緩釋性能的回報情況。參數(shù)b、k與材料添加量之間以二次函數(shù)擬合的相關(guān)性達(dá)到0.966 0以上，說明在緩釋回報率隨著用量的增加逐漸減少，當(dāng)用量增加到一定值后緩釋效果增強(qiáng)不明顯。

從本試驗(yàn)的結(jié)果看來，在考慮成本和緩釋效果等情況下，添加量在5.0%～6.0%之間較為適宜。

參考文獻(xiàn)

[1] 潘長兵. 橡膠園磚紅壤中氮、 磷、 鉀——經(jīng)濟(jì)管理學(xué)院地表徑流流失特征初步研究[D]. ?？冢?海南大學(xué)， 2012.

[2] 牛永生， 明大增， 李志祥， 等. 緩/控釋肥料的研究進(jìn)展[J]. 化學(xué)工業(yè)與工程技術(shù)， 2012， 33（6）： 36-39.

[3] 李 楊. 保水劑與肥料及土壤的互作機(jī)理研究[D]. 北京： 北京林業(yè)大學(xué)， 2012.

[4] 王 亮， 秦玉波， 于閣杰， 等. 新型緩控釋肥的研究現(xiàn)狀及展望[J]. 吉林農(nóng)業(yè)科學(xué)， 2008， 33（4）： 38-42.

[5] 谷佳林， 曹 兵， 李亞星， 等. 緩控釋氮素肥料的研究現(xiàn)狀與展望[J]. 土壤通報， 2008， 39（2）： 431-434.

[6] 劉 英， 熊海蓉， 李 霞， 等. 緩/控釋肥料的研究現(xiàn)狀及發(fā)展趨勢[J]. 化肥設(shè)計， 2013， 50（6）： 54-57.

[7] 王玉倩. 我國緩控釋肥行業(yè)現(xiàn)狀分析[J]. 化學(xué)工業(yè)， 2013， 6： 34-36.

[8] 史和平. 緩釋肥料肥效研究[J]. 現(xiàn)代農(nóng)業(yè)科技， 2014， 04： 214-216.

[9] 何緒生， 廖宗文， 黃培釗， 等. 保水緩/控釋肥料的研究進(jìn)展[J]. 農(nóng)業(yè)工程學(xué)報， 2006， 22（5）： 184-190.

[10] 黃幫裕. 保水緩釋肥料研究進(jìn)展[J]. 價值工程， 2013， 32（23）： 303-305.

[11] 莫 云， 周福亞， 尹應(yīng)武， 等. 高分子材料在保水緩/控釋肥料中的應(yīng)用進(jìn)展[J]. 中國農(nóng)業(yè)信息， 2013， 19： 85-86.

[12]張富倉， 李繼成， 雷 艷， 等. 保水劑對土壤保水持肥特性的影響研究[J]. 應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報， 2010， 18（1）： 120-128.

[13]黃 震， 黃占斌， 李文穎， 等. 不同保水劑對土壤水分和氮素保持的比較研究[J]. 中國生態(tài)農(nóng)業(yè)學(xué)報， 2010， 18（2）： 245-249.

[14] 王旭東， 關(guān)文玲， 楊雪芹， 等. 添加不同材料對肥料中氮磷在土壤中淋出特性的影響[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報， 2005， 24（3）： 544-547.

[15] 郝世雄. 包膜型緩釋/控釋尿素的研制[D]. 成都：四川大學(xué)， 2004.

[16] 王小利， 周建斌， 段建軍. 包膜控釋肥料氮素釋放動力學(xué)研究[J]. 楊凌：西北農(nóng)林科技大學(xué)學(xué)報（自然科學(xué)版）， 2003， 31（5）： 35-38， 42.

[17] 韓艷玉. 不同種類包膜肥料氮素溶出特性及對環(huán)境影響效應(yīng)的研究[D]. 沈陽： 沈陽農(nóng)業(yè)大學(xué)， 2011.

[18] 陳金佩. 緩/控釋肥料用粘土/聚丙烯酸（鉀）的制備及性能研究[D]. 合肥： 合肥工業(yè)大學(xué)， 2012.