劉仁浪　蔣翼杰　江美彥　陳靳松　肖婕妤　鄭全林　吳衛(wèi)

摘要： 在國家防止耕地“非糧化”政策下，為合理利用四川地區(qū)紫色土（坡地、荒地等）進(jìn)行白芷種植，本研究采用3因素5水平二次正交回歸旋轉(zhuǎn)組合設(shè)計，探究紫色土中氮肥、磷肥、鉀肥配施對川芷2號白芷品種產(chǎn)量和有效成分含量的影響。結(jié)果表明，紫色土中氮肥、磷肥、鉀肥配施對白芷產(chǎn)量的影響大小為鉀肥＞磷肥＞氮肥，鉀肥和磷肥對產(chǎn)量的影響都達(dá)到顯著或極顯著水平，氮肥對白芷產(chǎn)量沒有顯著影響；氮肥、磷肥、鉀肥配施對歐前胡素含量的影響為鉀肥>氮肥>磷肥，對異歐前胡素含量的影響為磷肥>氮肥>鉀肥，可見鉀肥顯著影響歐前胡素含量，磷肥顯著影響異歐前胡素含量。通過統(tǒng)計頻數(shù)法優(yōu)化出施肥配方，在氮肥施用量為233.73～283.36 kg/hm²、磷肥施用量為1 380.66～1 522.04 kg/hm²、鉀肥施用量為215.14～240.03 kg/hm²時川芷2號能實現(xiàn)高產(chǎn)，在氮肥施用量為240.39～276.58 kg/hm²、磷肥施用量為1 097.69～1 197.45 kg/hm²、鉀肥施用量為163.12～185.91 kg/hm²時可以均衡地提高川芷2號根中歐前胡素和異歐前胡素含量。最終利用灰色關(guān)聯(lián)度和逼近理想解排序法（DTOPSIS法）綜合獲得在紫色土中川芷2號的最優(yōu)施肥方案：氮肥施用量為193.40kg/hm²、磷肥施用量為1 449.26 kg/hm²、鉀肥施用量為239.88 kg/hm²。

關(guān)鍵詞： 紫色土；白芷；配施；產(chǎn)量；歐前胡素；異歐前胡素

中圖分類號： S567.9;S143 文獻(xiàn)標(biāo)識碼： A 文章編號： 1000-4440（2023）01-0208-10

Effects of combined application of nitrogen fertilizer， phosphorus fertilizer and potassium fertilizer on yield and active components of Angelica dahurica in purple soil

LIU Ren-lang¹， JIANG Yi-jie¹， JIANG Mei-yan¹， CHEN Jin-song¹， XIAO Jie-yu¹， ZHENG Quan-lin²，WU Wei¹

（1.College of Agronomy， Sichuan Agricultural University， Wenjiang 611130， China；2.Suining Tiandi Network Chuanbaizhi Industrial Co.， Ltd.， Suining 629000， China）

Abstract： Under the Chinese policy background of preventing non-grain conversion of cultivated land， a quadratic orthogonal regression rotation combination design using three factors and five levels was adopted in this study to make rational use of purple soil （sloping field， wasteland， etc.） in Sichuan area for Angelica dahurica var. formosana planting and to explore the effects of combined application of nitrogen fertilizer， phosphate fertilizer and potassium fertilizer in purple soil on the yield and effective components content of A. dahurica var. formosana variety of Chuanzhi 2. The results showed that， the effect of combined application of nitrogen fertilizer， phosphate fertilizer and potassium fertilizer on A. dahurica var. formosana yield in purple soil was ranked as follows： potassium fertilizer > phosphate fertilizer > nitrogen fertilizer. Both phosphate fertilizer and potassium fertilizer had significant or extremely significant effect on the yield of A. dahurica var. formosana， while nitrogen fertilizer had no significant effect on the yield of A. dahurica var. formosana. The effect of combined application of nitrogen fertilizer， phosphate fertilizer and potassium fertilizer on the imperatorin content was ranked as potassium fertilizer > nitrogen fertilizer > phosphate fertilizer， while the effect on isoimperatorin content was ranked as phosphate fertilizer > nitrogen fertilizer > potassium fertilizer. It can be seen that potassium fertilizer significantly affected imperatorin content， and phosphate fertilizer significantly affected isoimperatorin content. The fertilization formula was optimized by statistical frequency method， and it was found that Chuanzhi 2 could achieve high yield when the application levels of nitrogen fertilizer， phosphate fertilizer and potassium fertilizer were 233.73-283.36 kg/hm²， 1 380.66-1 522.04 kg/hm²and 215.14-240.03 kg/hm²， respectively. The contents of imperatorin and isoimperatorin in roots of Chuanzhi 2 could equably increase when the application levels of nitrogen fertilizer， phosphate fertilizer and potassium fertilizer were 240.39-276.58 kg/hm²， 1 097.69-1 197.45 kg/hm²and 163.12-185.91 kg/hm²， respectively. Finally， the optimal fertilization scheme for Chuanzhi 2 in purple soil was obtained by comprehensively using the grey relational degree and technique for order preference by similarity to ideal solution （DTOPSIS） method， and the application levels of nitrogenous fertilizer， phosphate fertilizer and potassium fertilizer were 193.40 kg/hm²， 1 449.26 kg/hm²and 239.88 kg/hm²， respectively.

Key words： purple soil;Angelica dahurica;combined application;yield;imperatorin;isoimperatorin

中藥白芷[Angelica dahurica （Fisch. ex Hoffm） Benth. et Hook. f.]藥食兩用，其原植物為傘形科當(dāng)歸屬植物白芷（A. dahurica）和杭白芷（A. dahurica var. formosana），具有驅(qū)風(fēng)散寒、清熱排膿、消腫止痛等功效^[1-2]。川白芷（四川）、祁白芷（河北）、禹白芷（河南）、杭白芷（浙江）是目前市場上主要流通的白芷品種^[3]。據(jù)《遂寧縣志》和《遂寧白芷志》記載，川白芷栽培歷史至今已有400～600年，是著名的川產(chǎn)道地藥材^[4-5]，一般種植于涪江兩岸的灰潮土上，但目前隨著城鎮(zhèn)化進(jìn)程的加快，白芷種植區(qū)域從涪江兩岸擴(kuò)張于四川內(nèi)地的紫色土區(qū)域^[6]?！度珖恋乩每傮w規(guī)劃綱要（2006-2020年） 》提出，要確保18億畝耕地紅線，而四川紫色土耕地集中分布于四川盆地丘陵區(qū)和海拔800 m以下的低山區(qū)，占全省耕地面積的68%^[7-8]。本課題組調(diào)研發(fā)現(xiàn)，在國家防止耕地“非糧化”政策下，川白芷種植于紫色土區(qū)域（坡地、荒地等）的面積越來越大，且市場認(rèn)可度和售價都超過用灰潮土種植的白芷。

施肥、產(chǎn)地土壤條件、生長環(huán)境、播種時期等因素會影響白芷的生長發(fā)育及次生代謝^[9-14]。前人對白芷施肥的相關(guān)研究中大多探究施肥措施對白芷抽薹苔率和產(chǎn)量的影響以及不同肥料種類對川白芷產(chǎn)量的影響，并未見關(guān)于紫色土中白芷配方施肥的研究報道^[9-10，15]。賈全全等^[16]在江西省紅壤丘陵區(qū)引種4大白芷品種（川白芷、杭白芷、禹白芷、祁白芷），發(fā)現(xiàn)杭白芷最適宜栽培于江西省紅壤丘陵區(qū)。因此應(yīng)依據(jù)白芷產(chǎn)地土壤和品種的特性制定科學(xué)合理的施肥措施，本研究采用3因素5水平二次正交回歸旋轉(zhuǎn)組合設(shè)計方法，探討紫色土中氮肥、磷肥、鉀肥配施對白芷品種川芷2號產(chǎn)量及品質(zhì)的影響，以期為川芷2號在紫色土上種植時合理施肥提供科學(xué)依據(jù)。

1 材料與方法

1.1 試驗材料及試驗點概況

本試驗于2019-2020年在川白芷道地產(chǎn)區(qū)四川省遂寧市船山區(qū)余健村（30°50′N，105°57′E）進(jìn)行。試驗地為平地，前茬為白芷，土質(zhì)為紫色土，土壤pH為8，有機(jī)質(zhì)含量13.41 g/kg，全氮含量0.96 g/kg，全磷含量0.73 g/kg，全鉀含量17.79 g/kg，堿解氮含量76.03 mg/kg，速效磷含量5.83 mg/kg，速效鉀含量144.10 mg/kg。

供試材料為四川農(nóng)業(yè)大學(xué)吳衛(wèi)教授選育的川芷2號白芷品種，原植物為杭白芷，其豐產(chǎn)性、穩(wěn)定性和適應(yīng)性良好^[17]。

供試肥料：氮肥為尿素（氮含量≥46.40%），磷肥為過磷酸鈣（含P₂O₅≥15%），鉀肥為硫酸鉀（含K₂O≥50%）。

1.2 試驗設(shè)計

采用3因素5水平二次正交回歸旋轉(zhuǎn)組合設(shè)計進(jìn)行田間試驗，設(shè)置氮肥（x₁）、磷肥（x₂）、鉀肥（x₃）3個因素，各5個水平，共23個小區(qū)隨機(jī)排列，每個小區(qū)面積為30 m²（長×寬=5 m×6 m）。二次正交回歸旋轉(zhuǎn)組合設(shè)計試驗因素水平編碼見表1，設(shè)計與試驗方案見表2。

1.3 田間管理

2019年10月采用條播方式將川芷2號播種于遂寧市船山區(qū)余健村，當(dāng)苗高4～7 cm時，進(jìn)行1～2次間苗，并按10 cm×40 cm株行距定苗。氮肥（尿素）作為苗肥（12月）、冬肥（次年2月）、春肥（次年3月）施用，施用比例為2∶5∶3；磷肥（過磷酸鈣）與鉀肥（硫酸鉀）作為底肥（9月）、春肥（次年3月）施用，施用比例為5∶5。其余田間栽培措施，均與當(dāng)?shù)剞r(nóng)戶保持一致。

1.4 測定指標(biāo)與測定方法

1.4.1 農(nóng)藝性狀及產(chǎn)量指標(biāo)的測定 2020年8月10日收獲時采用5點取樣法，每個小區(qū)取10株樣品，將白芷植株從莖基部分離為地上部分和地下部分。用清水洗干凈地下部分（根）后稱鮮質(zhì)量，用直尺測量根長，用游標(biāo)卡尺測量根直徑。將白芷根放置于105 ℃烘箱內(nèi)殺青15～30 min后，45 ℃烘干至恒質(zhì)量后稱其干質(zhì)量。各小區(qū)挖取3壟白芷稱量總的根鮮質(zhì)量，再折算出產(chǎn)量。

1.4.2 有效成分含量的測定 參照《中華人民共和國藥典》^[1]（簡稱《中國藥典》）中的方法測定白芷根部歐前胡素及異歐前胡素的含量。

1.5 數(shù)據(jù)分析與處理

采用Excel 2017、 DPS 7.05對數(shù)據(jù)進(jìn)行初步整理與分析，并建立回歸模型，使用Excel 2017和GraphPad Prism 8制圖。

根折干率=根干質(zhì)量（g）/根鮮質(zhì)量（g）×100%

產(chǎn)量（kg/hm²）=挖取白芷總鮮質(zhì)量（kg）/挖取面積（m²）×10 000

2 結(jié)果與分析

2.1 氮肥、磷肥、鉀肥配施對川芷2號產(chǎn)量的影響

2.1.1 模型的建立與檢驗 使用DPS7.05統(tǒng)計軟件對產(chǎn)量（Y₁）、根鮮質(zhì)量（Y₂）、根干質(zhì)量（Y₃）、折干率（Y₄）、根長（Y₅）、根直徑（Y₆）、歐前胡素含量（Y₇）、異歐前胡素含量（Y₈）的數(shù)據(jù)進(jìn)行回歸分析，得到氮肥（x₁）、磷肥（x₂）、鉀肥（x₃）施肥量的回歸模型。通過F檢驗分別驗證方程各項回歸系數(shù)、總回歸系數(shù)和失擬度，得到方差分析結(jié)果（表3）。

由表3可知，產(chǎn)量、折干率的方程總回歸系數(shù)FR（F1R=6.370，F(xiàn)4R=6.244）均達(dá)到極顯著水平，失擬度F1Lf未達(dá)到顯著水平（F1Lf=2.117），失擬度F4Lf達(dá)到顯著水平（F4Lf=3.272），說明產(chǎn)量的二次回歸方程模型擬合得較好，其可靠程度較高。根鮮質(zhì)量、根干質(zhì)量、根長、根直徑、歐前胡素含量、異歐前胡素含量方程的總回歸系數(shù)FR和失擬度FLf均不顯著，說明其二次回歸方程擬合效果不好，不可以進(jìn)行模型決策。在產(chǎn)量的二次回歸方程中，一次項（x2、x3）、二次項（x22、x32）、交互項（x2x3）均達(dá)到顯著或極顯著水平。因此，剔除在α=0.05顯著水平下不顯著的項后，可得到簡化后的產(chǎn)量（Y1′）二次回歸方程：Y1′=1 390.01+200.08x2+148.56x3-137.77x22-164.42x32+174.27x2x3。

2.1.2 因子主效應(yīng)分析 根據(jù)表3中各項回歸系數(shù)的F值，使用以下貢獻(xiàn)率計算公式得出各因素對因變量的貢獻(xiàn)率：

在紫色土中氮肥、磷肥、鉀肥對川芷2號產(chǎn)量的貢獻(xiàn)率分別為：△氮肥=0.346，△磷肥=1.812，△鉀肥=1.842，其貢獻(xiàn)率大小排序為鉀肥＞磷肥＞氮肥，其中鉀肥和磷肥對產(chǎn)量的影響最大，氮肥對產(chǎn)量的影響最小。

2.1.3 單因素效應(yīng)分析 由表3可知，單因素項中，氮肥對川芷2號產(chǎn)量沒有顯著影響，而磷肥和鉀肥對其產(chǎn)量的影響達(dá)到顯著或極顯著水平，其中磷肥的影響達(dá)到極顯著水平。通過降維法，得出磷肥、鉀肥對產(chǎn)量的單因素效應(yīng)方程，由方程作單因素效應(yīng)圖（圖1）。

由圖1可知，隨著磷肥和鉀肥施用水平的升高，川芷2號產(chǎn)量先增加后降低，符合肥料報酬定律，說明磷肥和鉀肥的施用范圍設(shè)計合理。磷肥和鉀肥在0～1水平時川芷2號的產(chǎn)量達(dá)到最高值，說明在紫色土中，當(dāng)磷肥施用量為1 099.44～1 449.26 kg/hm²、鉀肥施用量為179.91～239.88 kg/hm²時有利于川芷2號獲得高產(chǎn)。

2.1.4 互作效應(yīng)分析 由表3的回歸系數(shù)顯著性分析可知，磷肥-鉀肥互作對紫色土中川芷2號產(chǎn)量的影響達(dá)到了顯著水平，而氮肥-磷肥、氮肥-鉀肥互作對其產(chǎn)量影響不顯著。使用降維法固定氮肥因子使其為0，可得磷肥-鉀肥互作與產(chǎn)量的方程，根據(jù)方程可作出互作效應(yīng)圖（圖2）。

由圖2可知，不同磷肥、鉀肥配比對紫色土中種植的川芷2號產(chǎn)量影響差異較大，當(dāng)磷肥、鉀肥水平處于低磷高鉀和低鉀高磷時都對白芷產(chǎn)量有明顯的抑制作用。隨著磷肥和鉀肥施用水平的增加，產(chǎn)量也隨著增加。磷肥水平為1.68，鉀肥水平為1.34（即磷肥施用量為1 679.14 kg/hm²，鉀肥施用量為260.33 kg/hm²）時川芷2號產(chǎn)量最高。說明磷肥-鉀肥互作有利于川芷2號產(chǎn)量的增加。

2.1.5 川芷2號施肥方案的優(yōu)選 使用統(tǒng)計頻數(shù)法進(jìn)行分析可以得到一個合適的高產(chǎn)范圍。由當(dāng)?shù)赝寥馈夂驐l件可知小區(qū)產(chǎn)量高于 27 347.68 kg/hm²即為高產(chǎn)。因此，以產(chǎn)量方程作為相關(guān)函數(shù)目標(biāo)，采用統(tǒng)計頻數(shù)法優(yōu)選出川芷2號產(chǎn)量高于 27 347.68 kg/hm²的氮肥、磷肥、鉀肥配施方案。由表4可知，在紫色土中種植的川芷2號產(chǎn)量高于27 347.68 kg/hm²的方案有40個，通過計算95%的置信區(qū)間，可優(yōu)化出穩(wěn)定的施肥范圍，即氮肥施用量為233.73～283.36 kg/hm²、磷肥施用量為1 380.66～1 522.04 kg/hm²、鉀肥施用量為215.14～240.03 kg/hm²時，在紫色土中種植的川芷2號產(chǎn)量高于27 347.68 kg/hm²。

2.2 氮肥、磷肥、鉀肥配施對川芷2號有效成分的影響

2.2.1 因子主效應(yīng)分析 白芷有效藥用成分主要為歐前胡素、異歐前胡素等香豆素類化合物^[18]。利用2.1.2中貢獻(xiàn)率公式計算各因素對川芷2號有效成分指標(biāo)的貢獻(xiàn)率。由表5可知，各因素對川芷2號根中的有效成分的影響不同，對歐前胡素含量的影響為鉀肥>氮肥>磷肥，對異歐前胡素含量的影響為磷肥>氮肥>鉀肥。鉀肥對根中歐前胡素含量影響較大，而對根中異歐前胡素含量的影響最小，氮肥處于中間調(diào)控地位，磷肥對根中歐前胡素含量影響最小，而對根中異歐前胡素含量的影響較大。

2.2.2 單因素效應(yīng)分析 由表3可知，單因素項中，鉀肥顯著影響川芷2號根中歐前胡素含量，磷肥顯著影響川芷2號根中異歐前胡素含量，通過降維法得到鉀肥對歐前胡素含量及磷肥對異歐前胡素含量影響的單因素方程，由方程得到單因素效應(yīng)圖（圖3）。

由圖3可知，鉀肥對歐前胡素含量及磷肥對異歐前胡素含量的影響各不相同。由圖3A可知，鉀肥對歐前胡素含量的影響整體呈下降的趨勢，即隨著鉀肥施用水平加大，歐前胡素含量呈下降趨勢，說明施用過量的鉀肥不利于歐前胡素的積累，并對歐前胡素含量的積累有負(fù)面影響，鉀肥施用量與歐前胡素含量呈負(fù)相關(guān)。由圖3B可知，磷肥對異歐前胡素含量的影響呈先下降后上升趨勢，隨著磷肥施用水平的加大，異歐前胡素含量先降低再升高。磷肥施用量在-1.68～1.68水平時異歐前胡素含量先降低后再升高，1.68水平時異歐前胡素含量遠(yuǎn)遠(yuǎn)大于-1.68水平時異歐前胡素的含量，說明高磷水平比低磷水平更有利于異歐前胡素的積累，磷肥施用量與異歐前胡素含量呈正相關(guān)。

2.2.3 提高川芷2號有效成分含量的施肥方案的優(yōu)選 《中國藥典》規(guī)定，白芷根中歐前胡素含量不得少于0.08%^[1]。根據(jù)《中國藥典》的規(guī)定和當(dāng)?shù)厣a(chǎn)實際，將紫色土上種植的川芷2號根中的歐前胡素含量≥0.28%，異歐前胡素含量≥0.15%視為優(yōu)質(zhì)白芷標(biāo)準(zhǔn)。通過統(tǒng)計頻數(shù)法得到歐前胡素含量大于0.28%的方案有55個（表6）；異歐前胡素含量大于0.15%的方案有75個（表7）。

由表6、表7可知，通過計算95%的置信區(qū)間，得出當(dāng)?shù)适┯昧繛?33.73～277.06 kg/hm²，磷肥施用量為1 001.35～1 197.45 kg/hm²，鉀肥施用量為144.83～185.91 kg/hm²時川芷2號根中的歐前胡素含量可達(dá)到0.28%及以上；當(dāng)?shù)适┯昧繛?40.39～276.58 kg/hm²，磷肥施用量為1 097.69～1 317.05 kg/hm²，鉀肥施用量為163.12～196.70 kg/hm²時川芷2號根中的異歐前胡素含量可達(dá)到0.15%及以上。再進(jìn)行交集優(yōu)化后可得出當(dāng)?shù)适┯昧繛?40.39～276.58 kg/hm²，磷肥施用量為1 097.69～1 197.45 kg/hm²，鉀肥施用量為163.12～185.91 kg/hm²時可以均衡地提高紫色土中種植的川芷2號根中歐前胡素和異歐前胡素含量。

2.3 基于灰色關(guān)聯(lián)度和DTOPSIS法綜合評價川芷2號的產(chǎn)量及有效成分含量

為了精準(zhǔn)探究在紫色土中種植川芷2號的施肥配方，本研究綜合了產(chǎn)量、歐前胡素含量、異歐前胡素含量，通過灰色關(guān)聯(lián)度法得到其關(guān)聯(lián)度和權(quán)重（表8）。通過關(guān)聯(lián)度和權(quán)重進(jìn)行DTOPSIS法處理，根據(jù)DTOPSIS法中的C_i（相對接近度）進(jìn)行排序，并計算出待評材料的γ_i（關(guān)聯(lián)度）（表9）。

由表9可知，灰色關(guān)聯(lián)度法中最大關(guān)聯(lián)度差異為45.15%，而DTOPSIS法中最大C_i差異為71.03%，說明DTOPSIS法相比灰色關(guān)聯(lián)度法放大了處理間的差異性，更有利于區(qū)分待評材料間的優(yōu)劣。處理5、處理12和處理21等處理采用灰色關(guān)聯(lián)度和DTOPSIS法得到的綜合評價排序結(jié)果一致，其中綜合表現(xiàn)最優(yōu)的是處理5，即當(dāng)?shù)适┯昧繛?93.40 kg/hm²，磷肥施用量為1 449.26 kg/hm²，鉀肥施用量為239.88 kg/hm²（氮肥∶磷肥∶鉀肥=0.75∶1.81∶1.00）時在紫色土中川芷2號的產(chǎn)量和有效成分含量綜合最優(yōu)。

3 討論

遂寧市位于四川盆地中部丘陵低山地區(qū)，其土壤以紫紅色砂土和泥巖為主^[19]。由于川白芷種植區(qū)域由涪江沖積壩上（灰潮土）逐漸向丘陵低山地區(qū)擴(kuò)張，導(dǎo)致白芷種植在紫色土上的面積越來越大。本課題組在川芷2號的多點試驗中發(fā)現(xiàn)紫色土中種植川芷2號的產(chǎn)量與灰潮土基本一致，但其歐前胡素和異歐前胡素含量略高于灰潮土，可能是由于紫色土母巖疏松，易崩解，所以礦質(zhì)養(yǎng)分含量高，肥力好，有機(jī)質(zhì)含量較低，磷、鉀等含量豐富^[20]。本研究中在紫色土中種植的川芷2號氮肥、磷肥、鉀肥的配比為0.75∶1.81∶1.00，其產(chǎn)量和歐前胡素含量最高。陳郡雯等^[21]研究發(fā)現(xiàn)灰潮土中白芷氮肥、磷肥、鉀肥配比為1.3∶1.9∶1.0時白芷產(chǎn)量和歐前胡素含量最高。這可能是由于紫色土與灰潮土的土壤母質(zhì)、理化性質(zhì)、養(yǎng)分狀況及分布形態(tài)存在差異，從而對白芷產(chǎn)量和品質(zhì)產(chǎn)生不同影響^[22-24]；也可能是不同白芷品種需肥量有所差異。紫色土與灰潮土對川芷2號產(chǎn)量和品質(zhì)的影響因素還需進(jìn)一步探究。

土壤耕作層速效氮含量低，速效磷和速效鉀含量較為豐富是白芷高產(chǎn)的基礎(chǔ)^[11]。從本研究單因素分析中可知，氮肥對紫色土中種植的川芷2號產(chǎn)量影響不顯著，磷肥對產(chǎn)量有極顯著影響，鉀肥對產(chǎn)量有顯著影響，同時磷肥-鉀肥對紫色土中種植的川芷2號產(chǎn)量還有交互作用。翟娟園等^[13]研究發(fā)現(xiàn)，土壤有效磷及速效鉀的含量與白芷根質(zhì)量呈現(xiàn)正相關(guān)，這與本研究的結(jié)果一致。

白芷根中總香豆素主要為歐前胡素和異歐前胡素，歐前胡素和異歐前胡素含量能反映川白芷的品質(zhì)^[25-26]。賈蕾等^[27]研究發(fā)現(xiàn)，磷肥和鉀肥有利于提高白芷根中香豆素的含量。陳郡雯等^[28]研究發(fā)現(xiàn)，川白芷在生長過程中異歐前胡素含量與根部氮、磷、鉀的積累量呈顯著正相關(guān)。在本研究中川芷2號根中的歐前胡素和異歐前胡素含量均大于0.1%，符合《中國藥典》的規(guī)定。本研究中氮肥、磷肥、鉀肥配施對白芷歐前胡素含量的影響為鉀肥>氮肥>磷肥，對異歐前胡素含量的影響為磷肥>氮肥>鉀肥，磷肥顯著影響異歐前胡素含量，鉀肥顯著影響歐前胡素含量，磷肥施用量與異歐前胡素含量呈正相關(guān)，而鉀肥施用量與歐前胡素含量呈負(fù)相關(guān)。陳郡雯等^[21]在氮肥、磷肥、鉀肥配施對川白芷歐前胡素和異歐前胡素含量影響的研究中發(fā)現(xiàn)，鉀肥施用量與歐前胡素含量呈負(fù)相關(guān)，磷肥施用量與異歐前胡素含量呈正相關(guān)，本研究結(jié)果與其一致。因此在四川白芷主產(chǎn)區(qū)紫色土上種植白芷，只要選擇合理的氮肥、磷肥、鉀肥施用量及配比，就可以保障其產(chǎn)量和品質(zhì)，白芷產(chǎn)量最高可達(dá)27 347.68 kg/hm²，歐前胡素和異歐前胡素含量分別可達(dá)0.28%、0.15%。這為新產(chǎn)區(qū)白芷生產(chǎn)提供了科學(xué)施肥指導(dǎo)，為合理利用坡地、荒地等開展白芷生產(chǎn)奠定了基礎(chǔ)。

4 結(jié)論

本試驗通過3因素5水平二次正交回歸旋轉(zhuǎn)組合設(shè)計，建立了以氮肥、磷肥、鉀肥施用量為自變量，在紫色土中種植的白芷產(chǎn)量和有效成分含量為因變量的函數(shù)模型，通過模型分析得出，磷肥和鉀肥施用量對川芷2號產(chǎn)量有顯著影響，磷肥-鉀肥對產(chǎn)量還有交互作用。利用頻數(shù)分析法優(yōu)化出紫色土中種植的川芷2號高產(chǎn)的施肥配方：氮肥施用量為233.73～283.36 kg/hm²、磷肥施用量為1 380.66～1 522.04 kg/hm²、鉀肥施用量為215.142～240.03 kg/hm²。鉀肥施用量對歐前胡素含量有顯著作用，磷肥施用量對異歐前胡素含量有顯著作用，使用頻數(shù)分析法進(jìn)行交集優(yōu)化后可得出，當(dāng)?shù)适┯昧繛?40.39～276.58 kg/hm²，磷肥施用量為1 097.69～1 197.45 kg/hm²，鉀肥施用量為163.12～185.91 kg/hm²時可以均衡地提高紫色土中種植的川芷2號根中歐前胡素和異歐前胡素含量。最后通過灰色關(guān)聯(lián)度和DTOPSIS法綜合評價同時兼顧紫色土中川芷2號的產(chǎn)量和有效成分含量的最優(yōu)施肥方案：氮肥施用量為193.40kg/hm²，磷肥施用量為1 449.26 kg/hm²，鉀肥施用量為239.88 kg/hm²。

參考文獻(xiàn)：

[1] 中華人民共和國國家藥典委員會. 中華人民共和國藥典 [M]. 北京： 中國醫(yī)藥科技出版社， 2020.

[2] WEI W L， ZENG R， GU C M， et al. Angelica sinensis in China-A review of botanical profile， ethnopharmacology， phytochemistry and chemical analysis[J]. Journal of Ethnopharmacology， 2016， 190： 116-141.

[3] 張慶芝，王開疆，刀莉芳. 中藥白芷的品種論述 [J]. 云南中醫(yī)學(xué)院學(xué)報， 2000， 23（2）： 22-24.

[4] 王夢月，賈敏如. 白芷本草考證 [J]. 中藥材， 2004， 27（5）： 382-385.

[5] 遂寧市地方志編篡委員會. 遂寧縣志 [M]. 成都： 巴蜀書社， 1992.

[6] 鄭 利，鄧 聰，馮 亮，等. 遂寧川白芷產(chǎn)業(yè)發(fā)展現(xiàn)狀與分析 [J]. 農(nóng)村經(jīng)濟(jì)與科技， 2020， 31（20）： 181-182.

[7] 施業(yè)家，吳賢靜. 生態(tài)紅線概念規(guī)范化探討 [J]. 中南民族大學(xué)學(xué)報（人文社會科學(xué)版）， 2016， 36（3）： 149-153.

[8] 中國科學(xué)院成都分院土壤研究室. 中國紫色土（上篇）[M]. 北京： 科學(xué)出版社， 1991.

[9] 童 文，孫 佩，楊 曉，等. 施用腐植酸肥對白芷產(chǎn)量和質(zhì)量的影響 [J]. 西南農(nóng)業(yè)學(xué)報， 2011， 24（3）： 1236-1238.

[10]張亞琴，雷飛益，陳 雨，等. 鋅硼鉬配施對川白芷藥材農(nóng)藝性狀與產(chǎn)量的影響 [J]. 植物營養(yǎng)與肥料學(xué)報， 2018， 24（3）： 769-778.

[11]陳興福，丁德蓉，劉歲榮，等. 白芷生長土壤的研究 [J]. 中國中藥雜志， 1994， 19（10）： 591-593.

[12]陳興福，丁德蓉，劉歲榮，等. 白芷生態(tài)環(huán)境和土壤理化特性的研究 [J]. 中草藥， 1996， 27（8）： 489-591.

[13]翟娟園，吳 衛(wèi)，廖 凱，等. 土壤環(huán)境對川白芷產(chǎn)量和品質(zhì)的影響研究 [J]. 中草藥， 2010， 41（6）： 984-988.

[14]陳興福，盧 進(jìn)，丁德蓉. 播種期對白芷早期抽苔影響的研究 [J]. 中國中藥雜志， 1999， 24（4）： 211-212.

[15]丁德蓉，盧 進(jìn)，陳興福，等. 肥料種類對白芷早期抽苔與產(chǎn)量的影響研究 [J]. 中國中藥雜志， 1999， 24（1）： 24-25.

[16]賈全全，黃麗莉，楊春霞，等. 白芷在紅壤丘陵區(qū)栽培適應(yīng)性研究 [J]. 中藥材， 2019， 42（1）： 22-24.

[17]四川省種子站. 2012 年通過審定的農(nóng)作物新品種（七） [J]. 四川農(nóng)業(yè)科技，2013（5）： 12-13.

[18]吳媛媛，蔣桂華，馬逾英，等. 白芷的藥理作用研究進(jìn)展 [J]. 時珍國醫(yī)國藥， 2009， 20（3）： 625-627.

[19]景 可，鄭粉莉. 全國水土流失治理典型 [J]. 水土保持研究， 2004， 11（4）： 34-38.

[20]朱新玉，董志新，況福虹，等. 長期施肥對紫色土農(nóng)田土壤動物群落的影響 [J]. 生態(tài)學(xué)報， 2013， 33（2）： 464-474.

[21]陳郡雯，吳 衛(wèi)，侯 凱，等. N、P和K肥配施對川白芷產(chǎn)量及歐前胡素和異歐前胡素量的影響 [J]. 中草藥， 2011， 42（1）： 153-157.

[22]陳滄桑，孫錫發(fā)，何才富，等. 潮土、紫色土及黃壤上施鉀對小麥產(chǎn)量的影響 [J]. 西南農(nóng)業(yè)學(xué)報， 1999， 12（4）： 48-52.

[23]呂忠貴，胡定金，黃隆斌，等. 湖北省主要油菜種植區(qū)油菜及土壤營養(yǎng)元素含量狀況 [J]. 湖北農(nóng)業(yè)科學(xué)， 1997（3）：28-32.

[24]錢曉華，楊 平，周學(xué)軍，等. 安徽省土壤有效硫現(xiàn)狀及時空分布 [J]. 植物營養(yǎng)與肥料學(xué)報， 2018， 24（5）： 1357-1364.

[25]周愛德，李 強(qiáng)，雷海民. 白芷化學(xué)成分的研究 [J]. 中草藥， 2010， 41（7）： 1081-1083.

[26]王夢月，賈敏如. 白芷的化學(xué)成分研究進(jìn)展 [J]. 中藥材， 2002， 25（6）： 446-449.

[27]賈 蕾，孫淑蘭，李存東. 祁白芷的生長規(guī)律及氮磷鉀肥效應(yīng) [J]. 農(nóng)業(yè)與技術(shù)， 2009， 29（4）： 66-70.

[28]陳郡雯，吳 衛(wèi)，侯 凱，等. 川白芷生長發(fā)育、養(yǎng)分及有效成分的動態(tài)研究 [J]. 中國中藥雜志， 2010， 35（21）： 2812-2817.

（責(zé)任編輯：陳海霞）