張亞琴，雷飛益，陳 雨，李思佳，竇明明，馬留輝，石 峰，陳興福

（四川農(nóng)業(yè)大學(xué)農(nóng)學(xué)院/農(nóng)業(yè)部西南作物生理生態(tài)與耕作重點(diǎn)實(shí)驗(yàn)室，成都 611130）

白芷 (Angelica dahurica) 為傘形科當(dāng)歸屬植物[1]，其根藥食兩用，常用于治療感冒、偏頭痛、皮膚炎癥等疾病，也用作食物香料[2]。四川省遂寧市是白芷的道地產(chǎn)區(qū)，種植白芷的土壤為紫色土[3–4]。研究表明其有效態(tài)養(yǎng)分中缺硼和較缺鋅、鉬[5–6]，是遂寧白芷增產(chǎn)的限制因子之一。長(zhǎng)期栽培白芷，產(chǎn)地土壤日趨貧瘠，導(dǎo)致近年來(lái)白芷體積普遍瘦小，商品性狀不佳[7]，優(yōu)質(zhì)的川白芷資源常處于供不應(yīng)求的狀態(tài)[8]。為了滿足白芷生長(zhǎng)過(guò)程中對(duì)微量元素的需求，提高藥材川白芷的產(chǎn)量和商品性，增加藥農(nóng)的經(jīng)濟(jì)效益，在生產(chǎn)中補(bǔ)充施用微量元素肥料十分必要。

微量元素是植物生長(zhǎng)代謝及養(yǎng)分吸收運(yùn)轉(zhuǎn)相關(guān)酶的重要組成成分，在植物生長(zhǎng)發(fā)育過(guò)程中扮演著重要角色[9–11]，使用微量元素肥料來(lái)調(diào)節(jié)植物生長(zhǎng)在經(jīng)濟(jì)作物生產(chǎn)實(shí)踐中已較為廣泛[12–13]。中藥材作為一種特殊商品，始終遵循著“看貨論級(jí)，分檔議價(jià)”的指導(dǎo)思想，即藥材的農(nóng)藝性狀可以區(qū)分藥材質(zhì)量，從而決定了藥材的商品等級(jí)，是判別藥材價(jià)格檔次的重要指標(biāo)[14]。研究證明，施用微量元素肥料可以有效地調(diào)控藥材性狀，同時(shí)提高產(chǎn)量。孟杰等[15]研究表明噴施適宜配比的鋅、硼、鉬肥料有利于北柴胡 (Bupleurum chinense) 干物質(zhì)的向地下部分的轉(zhuǎn)移，可顯著提高北柴胡的產(chǎn)量。李瑤等[16]研究了化肥與微肥配施對(duì)川澤瀉 (Alisma orientalis) 的影響，發(fā)現(xiàn)化肥、微肥配施有利于增加澤瀉的產(chǎn)量，并優(yōu)化了最優(yōu)施肥量配方。羅意等[17]在附子 (Aconitum carmichaeli Debx) 上葉面噴施微肥，發(fā)現(xiàn)適宜濃度的鐵、鋅、硼、錳肥能提高附子的干物質(zhì)積累及產(chǎn)量，增加一等品的比例。綜合前人研究和當(dāng)前川白芷產(chǎn)區(qū)土壤基本情況，本研究選擇了鋅、硼和鉬三個(gè)因素研究其四川遂寧白芷的農(nóng)藝性狀和產(chǎn)量的影響，以期為該地區(qū)白芷生產(chǎn)上微肥的施用提供理論參考和技術(shù)支撐。

1 材料與方法

1.1 試驗(yàn)地概況

試驗(yàn)于2014年至2015年在四川省遂寧市大英縣隆盛鎮(zhèn)天丹中藥材科技有限責(zé)任公司試驗(yàn)基地(105°21′E，30°36′N(xiāo)) 進(jìn)行，該地區(qū)日平均溫度為17.5℃，年平均降雨量為925 mm，無(wú)霜期為298 d。試驗(yàn)地土壤pH 7.3，全氮 (N) 1.47 g/kg，全磷 (P)0.30 g/kg，全鉀 (K) 15.30 g/kg，堿解氮 (N) 14.60 mg/kg，有效磷 (P) 5.15 mg/kg，速效鉀 (K) 112.13 mg/kg，有效鋅 (Zn) 9.47 mg/kg，有效硼 (B) 0.22 mg/kg，有效鉬 (Mo) 0.38 mg/kg。

1.2 試驗(yàn)設(shè)計(jì)

供試材料經(jīng)四川農(nóng)業(yè)大學(xué)陳興福教授鑒定為傘形科當(dāng)歸屬植物白芷[Angelica dahurica (Fisch. ex Hoffm.) Benth.et Hook. f]。試驗(yàn)采用三因素二次回歸正交旋轉(zhuǎn)組合設(shè)計(jì)，設(shè)置鋅肥 (ZnSO4·7H2O)、硼肥(H3BO3) 和鉬肥[(NH4)2Mo4O13]三個(gè)因素，各5個(gè)水平，試驗(yàn)的23個(gè)小區(qū)隨機(jī)排列，另外設(shè)置1小區(qū)噴施等量的清水為對(duì)照。供試微肥均為分析純，含量均 ≥ 99.0%。試驗(yàn)設(shè)定的各因素水平編碼見(jiàn)表1，各小區(qū)微肥用量見(jiàn)表2。將表2中微肥量按小區(qū)面積折算用量并均分成3份，每份混合配制成500 mL溶液。分別在2015年4月12日﹑2015年4月27日和2015年5月12日清晨用微型噴霧器均勻噴施于川白芷葉面，噴施每小區(qū)時(shí)用塑料薄膜遮擋周邊小區(qū)。

表 1 因素水平編碼表(kg/hm2)Table 1 Factor level coding table

表 2 試驗(yàn)設(shè)計(jì)與實(shí)施方案Table 2 Experimental design and implementation plan

1.3 田間管理

白芷于2014年10月播種，行距0.5 m，株距0.27 m，小區(qū)面積為21 m2(5 m × 4.2 m)，每穴定苗3株。2015年4月12日起，每15天葉面噴施試驗(yàn)設(shè)計(jì)微肥量共3次，其余農(nóng)業(yè)栽培措施與當(dāng)?shù)厮庌r(nóng)習(xí)慣保持一致。

1.4 測(cè)定項(xiàng)目與方法

1.4.1 藥材性狀的測(cè)定 試驗(yàn)材料于2015年7月11日田間取樣，每小區(qū)采用五點(diǎn)法隨機(jī)取5株完整白芷植株，用清水洗凈，除去少量須根，晾干表面水分。將地上部分和地下部分從莖基部分開(kāi)，測(cè)量根長(zhǎng)、根直徑、根體積和根重。所有樣品置于105℃烘箱殺青15 min，45℃烘干至恒重后稱(chēng)干重。

1.4.2 產(chǎn)量測(cè)定 于2015年7月11日，白芷地上部分枯黃時(shí)取樣，將各小區(qū)內(nèi)川白芷全部取出，去掉地上部分，全部洗凈后稱(chēng)鮮重，45℃烘干至恒重后稱(chēng)干重。

1.4.3 數(shù)據(jù)分析與處理 使用Excel 2013軟件進(jìn)行數(shù)據(jù)的初步換算和處理，使用DPS 7.05數(shù)據(jù)處理系統(tǒng)對(duì)數(shù)據(jù)進(jìn)行回歸模型的建立和分析，最后用Excel 2013軟件制圖。根據(jù)川白芷生產(chǎn)要求，參考孟杰等[18]對(duì)藥材綜合品質(zhì)評(píng)價(jià)方法，結(jié)合灰色關(guān)聯(lián)度和DTOPSIS法對(duì)川白芷外觀性狀，包括根體積、根密度、根折干率、根長(zhǎng)、根直徑和產(chǎn)量6個(gè)指標(biāo)進(jìn)行權(quán)重計(jì)算和綜合評(píng)價(jià)。

根折干率 (%) = 平均根干重/平均根鮮重 × 100

根密度 (g/cm3) = 平均單株鮮重/平均單株體積

川白芷產(chǎn)量 (kg/hm2) = 平均單株根重 × 每公頃株數(shù)

2 結(jié)果與分析

2.1 模型的建立與檢驗(yàn)

使用DPS 7.05統(tǒng)計(jì)軟件對(duì)所得數(shù)據(jù)進(jìn)行回歸分析，分別得到鮮產(chǎn)量 (Y1)、干產(chǎn)量 (Y2)、根折干率(Y3)、根體積 (Y4)、根密度 (Y5)、根長(zhǎng) (Y6)、根直徑(Y7) 與鋅 (X1)、硼 (X2)、鉬 (X3) 的回歸模型。用F檢驗(yàn)法分別檢驗(yàn)方程各項(xiàng)回歸系數(shù)、方程總回歸系數(shù)和失擬度，得到方差分析表 (表3)。

由表3可知，在顯著水平α = 0.10的水平下通過(guò)方差分析發(fā)現(xiàn)方程Y1、Y2、Y4、Y5的總回歸系數(shù)F1均達(dá)到了顯著水平，方程Y1、Y2、Y4的失擬度未達(dá)到顯著水平，這說(shuō)明方程Y1、Y2、Y4的模型成立，具有較好的預(yù)測(cè)性，可以進(jìn)行模型決策。剔除顯著水平α ＜ 0.10不顯著的項(xiàng)后，得到優(yōu)化的回歸方程 Y1′、Y2′和 Y4′。

表 3 方差分析表Table 3 Analysis of variance

2.2 微肥配施對(duì)川白芷外觀性狀的影響

2.2.1 因子主效應(yīng)分析 利用表3中各項(xiàng)回歸系數(shù)的F值計(jì)算各因素對(duì)因變量的貢獻(xiàn)率，可分析回歸方程中各因素的重要性。使用貢獻(xiàn)率計(jì)算公式計(jì)算，分別得到各因素對(duì)川白芷外觀性狀的貢獻(xiàn)率 (表4)。

由表4可知，各因素對(duì)川白芷外觀性狀的影響不同。對(duì)根折干率、根體積、根長(zhǎng)和根直徑的影響表現(xiàn)為：鉬 ＞ 鋅 ＞ 硼。可以看出，硼肥對(duì)外觀性狀的影響最小，鉬肥對(duì)外觀品質(zhì)影響較大，鋅肥處于中間調(diào)控的地位。對(duì)根密度的影響表現(xiàn)為：硼 ＞鋅 ＞ 鉬。

2.2.2 單因素效應(yīng)分析 使用降維法，固定其他因子為零水平，可分析鋅、硼、鉬對(duì)川白芷外觀性狀的影響。由表3可知，單因素項(xiàng)中，鋅分別對(duì)根體積和根直徑有顯著影響，硼分別對(duì)根密度和根長(zhǎng)有顯著影響，鉬分別對(duì)根體積和根密度有顯著影響。由此分別計(jì)算出鋅、硼、鉬肥對(duì)根體積、根密度、根長(zhǎng)和根直徑的單因素效應(yīng)方程，由方程作單因素圖(圖 1)。

表 4 單因素貢獻(xiàn)率Table 4 Single factor contribution rate

由圖1A可看出，在本研究地區(qū)的地力水平下，鉬肥對(duì)川白芷的根體積的影響存在一個(gè)最佳濃度值X3= 0.90，即鉬肥濃度為0.096 kg/hm2時(shí)，其根體積最大，說(shuō)明噴施適宜濃度的鉬肥有利于根體積的增加。由圖1B可看出隨著鉬肥水平的增加，根密度先降低后增加，與鉬肥對(duì)根體積的影響恰好相反，說(shuō)明鉬肥在增加川白芷根體積的同時(shí)，根密度在降低。由圖1A可知，隨著鋅肥濃度的升高，根體積逐漸降低。由圖1B、C可看出隨著硼肥水平的增加，根密度逐漸降低，根長(zhǎng)在縮短，這說(shuō)明硼肥抑制了川白芷根密度的增長(zhǎng)和根的伸長(zhǎng)，濃度越高，抑制效果越明顯。由圖1D可看出，鋅肥對(duì)川白芷根直徑存在最佳濃度值X1= – 1.18，即鋅肥濃度為0.043 kg/hm2時(shí)，其根最粗。這說(shuō)明適宜濃度的鋅肥有利于川白芷根的增粗，但對(duì)其體積的增加有抑制作用。2.2.3 互作效應(yīng)分析 由表3可知，鋅、鉬互作對(duì)根體積有顯著影響，硼、鉬互作對(duì)根密度有較顯著影響。使用降維法固定單因子為零，可得到鋅、鉬互作對(duì)根體積和硼、鉬互作對(duì)根密度的方程，分別作互作效應(yīng)圖 (圖2、圖3)。

由圖2可知，在鉬肥水平為1.68，鋅肥水平為–1.68時(shí)，川白芷根體積最高。當(dāng)鉬肥處于1水平以上時(shí)隨著鋅肥水平的增加，根體積越來(lái)越低。當(dāng)鋅肥處于1以上水平時(shí)，隨著鉬肥水平的增加，根體積越來(lái)越低。在鋅肥和鉬肥均處于1.68水平時(shí)體積最小。這說(shuō)明鋅、鉬互作對(duì)川白芷根體積增長(zhǎng)有抑制作用。

由圖3可知，當(dāng)硼肥和鉬肥處于 –1.68水平時(shí)，川白芷根密度最大。當(dāng)硼肥保持一定水平時(shí)，鉬肥水平越高，川白芷根密度先降低后升高，在鉬肥為零水平左右時(shí)達(dá)到最低。這說(shuō)明鉬肥與硼肥互作不利于川白芷根密度的增加，尤其是在鉬肥處于零水平左右時(shí)抑制作用最為明顯。

2.3 微肥配施對(duì)川白芷產(chǎn)量的影響

2.3.1 因子主效應(yīng)分析 使用2.2.1中貢獻(xiàn)率法計(jì)算得到各因素對(duì)產(chǎn)量的貢獻(xiàn)率為Δ鋅=1.713，Δ硼=0.725，Δ鉬=1.985。按從大到小排序?yàn)殂f肥 ＞ 鋅肥 ＞ 硼肥，說(shuō)明鉬肥對(duì)產(chǎn)量的影響最大，鋅肥和硼肥次之。2.3.2 單因素效應(yīng)分析 由表3分析可知，鋅硼鉬配施對(duì)川白芷的折干率無(wú)影響，且川白芷的鮮產(chǎn)量和干產(chǎn)量回歸方程趨勢(shì)一致，由此本文僅對(duì)干產(chǎn)量進(jìn)行分析，以下簡(jiǎn)稱(chēng)產(chǎn)量。使用降維法對(duì)產(chǎn)量的回歸方程進(jìn)一步分析，可分別得到鋅、硼、鉬肥對(duì)產(chǎn)量的單因素效應(yīng)方程，分別由方程作單因素圖 (圖4)。

圖 1 單因素效應(yīng)Fig. 1 Single factor effect

圖 2 鋅、鉬互作對(duì)根體積的影響Fig. 2 Effect of Zn-Mo interaction on root volume

圖 3 硼、鉬互作對(duì)根密度的影響Fig. 3 Effect of B-Mo interaction on root density

由圖4可知，隨著鋅肥水平的升高，產(chǎn)量逐漸降低，施用鋅肥對(duì)產(chǎn)量有抑制作用；隨著鉬肥水平的增加，產(chǎn)量先增加后降低，鉬肥在0.5～1.5水平川白芷產(chǎn)量達(dá)到最高，這說(shuō)明0.139～0.203 kg/hm2的鉬肥水平可使產(chǎn)量達(dá)到較高水平。

圖 4 產(chǎn)量單因素效應(yīng)Fig. 4 Single factor effect of yield

2.3.3 互作效應(yīng)分析 由表3的回歸系數(shù)顯著性分析可知，鋅和硼、鋅和鉬對(duì)川白芷產(chǎn)量的影響達(dá)到了顯著水平，硼和鉬對(duì)產(chǎn)量的影響不顯著。通過(guò)降維法對(duì)互作效應(yīng)進(jìn)一步分析，可得到鋅和硼互作、鋅和鉬互作與川白芷產(chǎn)量的關(guān)系方程，由此可根據(jù)方程作出互作效應(yīng)的曲面圖 (圖5A和圖5B)。

由圖5A可知在高水平的鋅肥和硼肥水平下川白芷產(chǎn)量最低，而高水平的硼肥和低水平的鋅肥條件下川白芷產(chǎn)量達(dá)到最高，即1.682水平的硼肥和0水平的鋅肥條件下產(chǎn)量可達(dá)到極值11626 kg/hm2。隨著硼肥水平的提高，鋅肥水平越低，產(chǎn)量越高。這說(shuō)明高水平的鋅肥和硼肥配合不利于川白芷產(chǎn)量的增加。由圖5B可知鋅肥和鉬肥的互作效應(yīng)與鋅肥和硼肥類(lèi)似，在1.682水平的硼肥和0水平的鋅肥水平下產(chǎn)量最高可達(dá)12345 kg/hm2。這說(shuō)明在互作關(guān)系中，鋅肥對(duì)硼、鉬肥料的肥效產(chǎn)生了抑制作用。

2.3.4 川白芷產(chǎn)量?jī)?yōu)化配方分析 通過(guò)已建立的產(chǎn)量的回歸模型，可計(jì)算出各因素的最適施用配方的理論值。考慮到實(shí)際生產(chǎn)時(shí)存在土壤、氣候環(huán)境的差異影響，通過(guò)使用統(tǒng)計(jì)頻數(shù)法進(jìn)行分析可以得到一個(gè)合適的高產(chǎn)范圍。按照當(dāng)?shù)赝寥?、氣候條件可知小區(qū)產(chǎn)量高于8250 kg/hm2即為高產(chǎn)，在此使用統(tǒng)計(jì)頻數(shù)法計(jì)算出產(chǎn)量高于8250 kg/hm2的最佳施肥量范圍 (表 5)。

圖 5 鋅硼、鋅鉬互作效應(yīng)分析Fig. 5 Interactions analysis of Zn-B, Zn-Mo

由表5可知，當(dāng)鋅肥在低水平編碼 – 1.682～ – 1時(shí)頻率最高，即鋅肥施用濃度在0～0.058 kg/hm2間，川白芷的產(chǎn)量呈現(xiàn)較高狀態(tài)；硼肥頻率分布較平均，這說(shuō)明施用硼肥濃度對(duì)產(chǎn)量的影響并不大；鉬肥在中高水平編碼1～1.682間頻率較大，即當(dāng)鉬肥濃度在0.171～0.214 kg/hm2時(shí)川白芷產(chǎn)量較高。綜合以上的頻數(shù)分析，可以優(yōu)化最佳農(nóng)藝措施為鋅肥0.0305～0.0960 kg/hm2，硼肥1.0270～1.7928 kg/hm2，鉬肥 0.1354～0.1774 kg/hm2。

2.4 基于灰色關(guān)聯(lián)度和DTOPSIS法綜合評(píng)價(jià)川白芷外觀性狀和產(chǎn)量

為了用川白芷的各項(xiàng)指標(biāo)對(duì)川白芷整體做出一個(gè)客觀評(píng)價(jià)，本文綜合了產(chǎn)量、根折干率、根體積、根密度、根長(zhǎng)、根直徑，使用灰色關(guān)聯(lián)度法計(jì)算出各參評(píng)材料的權(quán)重值和灰關(guān)聯(lián)度 (表6)。使用灰關(guān)聯(lián)權(quán)重值進(jìn)行了DTOPSIS法處理，得到DTOPSIS法的相對(duì)接近度值Ci及其排序 (表7)，并計(jì)算出待評(píng)材料與理想材料的關(guān)聯(lián)度 (表7)。

由表7可知，采用灰色關(guān)聯(lián)度法和DTOPSIS法的綜合評(píng)價(jià)結(jié)果幾乎一致。兩種方法所得綜合表現(xiàn)最優(yōu)的處理5和處理7，即當(dāng)使用鋅肥0.06 kg/hm2、硼肥2.28 kg/hm2、鉬肥0.17 kg/hm2和鋅肥0.06 kg/hm2、硼肥0.58 kg/hm2、鉬肥0.17 kg/hm2時(shí)川白芷的經(jīng)濟(jì)性狀和產(chǎn)量綜合最優(yōu)，這與頻數(shù)分析法得到的最優(yōu)方案相同。這兩種配方僅硼肥施用量不同，這印證了單因素分析過(guò)程中，硼對(duì)產(chǎn)量無(wú)影響的結(jié)果。由表7可知灰色關(guān)聯(lián)度法所得關(guān)聯(lián)度最大差異為25.39%，DTOPSIS法所得Ci值最大差異可達(dá)64.04%，這說(shuō)明DTOPSIS法能放大樣品間的差異性，更便于區(qū)分待評(píng)材料間表現(xiàn)的優(yōu)劣。

3 討論

微量元素在植物生長(zhǎng)代謝過(guò)程中起到重要的作用，缺乏微量元素會(huì)影響植物根系營(yíng)養(yǎng)供給和生理活動(dòng)，過(guò)量時(shí)植物會(huì)受到毒害，而每種植物對(duì)微量元素的耐受性和需求量不同，不同地區(qū)的土壤微量元素含量也存在較大差異[19]，施用適合的微量元素用量才能有效調(diào)控植物的生長(zhǎng)。本研究的三個(gè)因素鋅、硼、鉬對(duì)植物生長(zhǎng)發(fā)育影響較大[20–22]，已在眾多植物中得到印證[23–25]。本研究結(jié)果表明，各因素對(duì)產(chǎn)量、根體積等的貢獻(xiàn)率排序?yàn)殂f肥 ＞ 鋅肥 ＞ 硼肥，對(duì)根密度的貢獻(xiàn)率排序?yàn)殇\肥 ＞ 鉬肥 ＞ 硼肥。單因素分析結(jié)果表明，鉬肥能降低川白芷根的密度、增加體積和增產(chǎn)，在施用鉬肥0.1 kg/hm2左右濃度時(shí)體積最大，密度最小，在0.17 kg/hm2濃度左右，川白芷產(chǎn)量最高。這與鉬肥對(duì)甘草[26]和萵筍[27]的影響相同。這可能與鉬肥促進(jìn)了植物光合作用，有利于植物干物質(zhì)積累有關(guān)[28]，鉬肥對(duì)增加根類(lèi)藥材的性狀和產(chǎn)量有較好的提升效果。鋅肥能夠增加川白芷根的長(zhǎng)度，在濃度為0.043 kg/hm2時(shí)根長(zhǎng)最長(zhǎng)，但是施用鋅肥會(huì)抑制根體積和產(chǎn)量的增長(zhǎng)，鋅和鉬互作對(duì)根體積和產(chǎn)量有抑制效應(yīng)，鋅和硼互作對(duì)產(chǎn)量也產(chǎn)生了抑制效應(yīng)。這與柴胡[15]、澤瀉[29]不同，可能是由于白芷對(duì)鋅肥的耐受性較低，使鋅肥對(duì)白芷的根系干物質(zhì)積累產(chǎn)生了抑制，也影響了鉬肥對(duì)根體積和產(chǎn)量的促進(jìn)效果。硼對(duì)川白芷的根密度和根直徑都有抑制作用，對(duì)產(chǎn)量影響較小，硼和鉬互作會(huì)影響根密度的增加。這可能與硼肥移動(dòng)性較弱，對(duì)根部干物質(zhì)積累的控制力度不夠?qū)е拢疸f互作影響根密度可能與硼影響細(xì)胞伸展有關(guān)，影響了鉬肥對(duì)根密度的調(diào)控效果[30]。已有研究證明硼肥對(duì)川白芷有效成分含量影響較大，能有效增加川白芷中香豆素成分的含量[31]，這說(shuō)明硼肥對(duì)有效成分含量的提高效果更佳，其原因有待進(jìn)一步研究探索。藥材市場(chǎng)上通常按白芷的粗細(xì)、重量、粉性和香氣分3等，一等品的價(jià)格最貴而三等品的藥效成分更高[32]。以本試驗(yàn)結(jié)果看，鉬肥既能增加藥材性狀，提高白芷的商品等級(jí)和經(jīng)濟(jì)效益，也能有效增產(chǎn)；鋅肥能影響根長(zhǎng)，也能提高白芷等級(jí)，從而增加經(jīng)濟(jì)效益。

表 5 產(chǎn)量 ≥ 8250 kg/hm2的頻率分布及農(nóng)藝措施Table 5 The frequency distribution of yield ≥ 8250 kg/hm2 and agronomic measures

表 6 參評(píng)材料的關(guān)聯(lián)度γi和權(quán)重值ωkTable 6 Correlative degree weight value of the evaluation materials

表 7 待評(píng)價(jià)材料的加權(quán)關(guān)聯(lián)度ri與Ci排序值比較Table 7 Comparison of the weighted correlation degree ri and Ci sorting value of the material

鋅硼鉬的調(diào)控效果不同，但是三者協(xié)同作用時(shí)可使川白芷達(dá)到高產(chǎn)的目標(biāo)。竇明明等[29]研究微肥和化肥在澤瀉上的應(yīng)用，結(jié)果表明適宜配比的鋅、硼、鉬肥料和化肥配施可顯著提高澤瀉的產(chǎn)量。羅意等[17]研究發(fā)現(xiàn)適當(dāng)濃度鐵、錳、鋅、硼的微肥配比能顯著提高附子產(chǎn)量，并優(yōu)化出4∶1∶1∶2的附子最優(yōu)微肥施用配方。本研究結(jié)果表明，微肥鋅、硼、鉬配施可調(diào)控川白芷根密度、根體積、根長(zhǎng)和根粗，增加川白芷產(chǎn)量，并優(yōu)化出微肥配施方案：在噴施鋅肥0.0305～0.0960 kg/hm2、硼肥1.0270～1.7928 kg/hm2、鉬肥0.1354～0.1774 kg/hm2時(shí)，川白芷能夠高產(chǎn)。本文還使用了灰色關(guān)聯(lián)度和DTOPSIS法綜合分析比較了本試驗(yàn)條件下所有處理的優(yōu)劣，結(jié)果表明，使用鋅肥0.06 kg/hm2、硼肥2.28 kg/hm2、鉬肥0.17 kg/hm2和鋅肥0.06 kg/hm2、硼肥0.58 kg/hm2、鉬肥0.17 kg/hm2時(shí)川白芷的藥材性狀和產(chǎn)量綜合最優(yōu)，這與頻數(shù)分析法的結(jié)果一致，也說(shuō)明了硼肥對(duì)產(chǎn)量的影響很小。但因素之間是通過(guò)什么途徑對(duì)川白芷產(chǎn)量產(chǎn)生影響，怎樣協(xié)同互作還需進(jìn)一步研究。

4 結(jié)論

通過(guò)正交旋轉(zhuǎn)試驗(yàn)設(shè)計(jì)，本研究建立了以鋅、硼、鉬為自變量，川白芷藥材性狀和產(chǎn)量為因變量的函數(shù)模型。通過(guò)對(duì)模型分析得出，鋅、鉬對(duì)川白芷藥材性狀和產(chǎn)量有顯著影響，且因素間存在互作效應(yīng)。三因素對(duì)產(chǎn)量和根體積等的影響排序?yàn)殂f ＞鋅 ＞ 硼，對(duì)根密度的影響排序?yàn)殇\ ＞ 鉬 ＞ 硼。

鋅、硼、鉬三種微肥配施能調(diào)節(jié)川白芷藥材性狀，增加川白芷的產(chǎn)量，但存在一個(gè)合適的濃度配比，施用濃度過(guò)高或過(guò)低都會(huì)影響其產(chǎn)量和藥材農(nóng)藝性狀。本研究使用頻數(shù)分析法和灰色關(guān)聯(lián)度和DTOPSIS法綜合優(yōu)化出川白芷施肥配方為鋅肥0.0305～0.0960 kg/hm2、硼肥 1.0270～1.7928 kg/hm2、鉬肥0.1354～0.1774 kg/hm2。川白芷生產(chǎn)上可根據(jù)土壤、產(chǎn)量、經(jīng)濟(jì)效益預(yù)期等，依據(jù)此配方調(diào)整鋅、硼、鉬肥施用量，以提高川白芷產(chǎn)量。

[1]中國(guó)科學(xué)院中國(guó)植物志編輯委員會(huì). 中國(guó)植物志[M]. 北京: 科學(xué)出版社, 1992.Editorial Committee of Flora of China. Flora Reipublicae popularis sinicae[M]. Beijing: Science Press, 1992.

[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]郭志剛. 遂寧市天福鎮(zhèn)土壤質(zhì)量地球化學(xué)評(píng)價(jià)[D]. 成都: 成都理工大學(xué)碩士學(xué)位論文, 2009.Guo Z G. Geochemical evaluation of soil quality in Tianfu town of Suining city [D]. Chengdu: MS Thesis of Chengdu University of Technology, 2009.

[4]楊遠(yuǎn)祥, 朱雪梅, 邵繼榮, 等. 遂寧組母質(zhì)封禁后的土壤侵蝕規(guī)律分析[J]. 水土保持通報(bào), 2007, 27(1): 24–28.Yang Y X, Zhu X M, Shao J R, et al. Analysis of soil erosion on suining group parent material after closing management[J]. Bulletin of Soil and Water Conservation, 2007, 27(1): 24–28.

[5]易德超. 川西南土壤微量元素特征[J]. 中國(guó)土壤與肥料, 2001, (6):24–26.Yi D C. The characteristics of soil microelements in southwest of Sichuan[J]. Soil and Fertilizer Sciences in China, 2001, (6): 24–26.

[6]許偉, 劉應(yīng)平, 彭培好, 陳文德. 四川省簡(jiǎn)陽(yáng)市三岔鎮(zhèn)紫色土微量元素全量與有效性評(píng)價(jià)[J]. 物探化探計(jì)算技術(shù), 2011, 33(3): 335–339.Xu W, Liu Y P, Peng P H, Chen W D. Evaluation on total and available content of micro-elements in perple soil of Sancha town in Jianyang city of Sichuan province[J]. Computing Techniques for Geophysical and Geochemical Exploration, 2011, 33(3): 335–339.

[7]桂鏡生. 中藥商品學(xué)[M]. 昆明: 云南大學(xué)出版社, 2015.Gui J S. Traditional Chinese medicine commodity science[M].Kunming: Yunnan University Press, 2015.

[8]馬逾英, 郭丁丁, 蔣桂華, 等. 白芷種質(zhì)資源的調(diào)查報(bào)告[J]. 華西藥學(xué)雜志, 2009, 24(5): 457–460.Ma Y Y, Guo D D, Jiang G H, et al. Investigation report on the germplasm resources of Angelica dahurica[J]. West China Journal of Pharmaceutical Sciences, 2009, 24(5): 457–460.

[9]黃益宗. 植物對(duì)硼素不足的反應(yīng)及其成因探討[J]. 土壤與環(huán)境,2002, 11(4): 434–438.Huang Y Z. Response of plant to low boron and factors of boron deficiency[J]. Soil and Environmental Sciences, 2002, 11(4):434–438.

[10]何忠俊, 華珞, 白玲玉, 等. 土壤-植物系統(tǒng)中氮鋅交互作用研究進(jìn)展[J]. 土壤與環(huán)境, 2001, 10(2): 133–137.He Z J, Hua L, Bai L Y, et al. Research progress of interaction between nitrogen and zinc in plant-soil system[J]. Soil and Environmental Sciences, 2001, 10(2): 133–137.

[11]何忠俊, 梁社往, 曾波, 等. 鉬對(duì)滇重樓生長(zhǎng)、養(yǎng)分和總皂甙含量的影響[J]. 植物營(yíng)養(yǎng)與肥料學(xué)報(bào), 2011, 17(6): 1481–1486.He Z J, Liang S W, Zeng B, et al. Effects of Mo fertilization on growth, nutrient and total saponin contents of Pairs polyphylla var.yunnanensis[J]. Plant Nutrition and Fertilizer Science, 2011, 17(6):1481–1486.

[12]Chatterjee R, Bandyopadhyay S. Effect of boron, molybdenum and biofertilizers on growth and yield of cowpea (Vigna unguiculata L.Walp.) in acid soil of eastern Himalayan region[J]. Journal of the Saudi Society of Agricultural Sciences, 2017, 16(4): 332–336.

[13]Liu L, Xiao W, Ji M L, et al. Effects of molybdenum on nutrition,quality, and flavour compounds of strawberry (Fragaria×ananassa Duch. cv. Akihime) fruit[J]. Journal of Integrative Agriculture, 2017,16(7): 1502–1512.

[14]秦雪梅, 孔增科, 張麗增, 等. 中藥材“辨狀論質(zhì)”解讀及商品規(guī)格標(biāo)準(zhǔn)研究思路[J]. 中草藥, 2012, 43(11): 2093–2098.Qing X M, Kong Z K, Zhang L Z, et al. Interpretation of “assessing quality by distinguishing TCMM features” and research ideas about commodity specifications[J]. Chinese Traditional and Herbal Drugs,2012, 43(11): 2093–2098.

[15]孟杰, 陳興福, 楊文鈺, 等. 鋅硼鉬配施對(duì)北柴胡產(chǎn)量及柴胡皂苷a,d含量的影響[J]. 中國(guó)中藥雜志, 2014, 39(22): 4297–4303.Meng J, Chen X F, Yang W Y, et al. Effect of combined application of zinc, boron and molybdenum on yield and saikosaponin a,saikosaponin d contents of Bupleurum chinense[J]. China Journal of Chinese Materia Medica, 2014, 39(22): 4297–4303.

[16]李瑤, 陳興福, 彭世明, 等. 二次正交旋轉(zhuǎn)優(yōu)化澤瀉施肥[J]. 中藥材,2015, 38(4): 664–668.Li Y, Chen X F, Peng S M, et al. Quadratic orthogonal rotation combination design on Alisma orientalis of fertilization[J]. Journal of Chinese Medicinal Materials, 2015, 38(4): 664–668.

[17]羅意, 陳興福, 劉莎, 等. 微肥配施對(duì)附子產(chǎn)量和質(zhì)量的調(diào)控效應(yīng)[J]. 中國(guó)中藥雜志, 2011, 36(2): 102–106.Luo Y, Chen X F, Liu S, et al. Effects of micronutrient fertilizer application on yield and quality of Aconitum carmichaeli[J]. China Journal of Chinese Materia Medica, 2011, 36(2): 102–106.

[18]孟杰, 陳興福, 楊文鈺, 等. 基于灰色關(guān)聯(lián)度分析和DTOPSIS法綜合評(píng)價(jià)青川柴胡資源[J]. 中國(guó)中藥雜志, 2014, 39(3): 433–437.Meng J, Chen X F, Yang W Y, et al. Evaluation of Bupleuri Radix resources in Qingchuan based on DTOPSIS and grey related degree[J]. China Journal of Chinese Materia Medica, 2014, 39(3):433–437.

[19]謝佰承, 張春霞, 薛緒掌. 土壤中微量元素的環(huán)境化學(xué)特性[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào), 2007, 26(S1): 132–135.Xie B C, Zhang C X, Xue X Z. Characteristics of environmental chemistry for trace elements in soil[J]. Journal of Agro-Environment Science, 2007, 26(S1): 132–135.

[20]Camacho-Crist6bal J J, Rexach J, González-Fontess A. Boron in plants: deficiency and toxicity[J]. Journal of Integrative Plant Biology, 2008, 50(10): 1247–1255.

[21]Broadley M R, White P J, Hammond J P, et al. Tansley review: zinc in plants[J]. New Phytologist, 2007, 173(4): 677–702.

[22]Mendel R R, H?nsch R. Molybdoenzymes and molybdenum cofactor in plants[J]. Journal of Experimental Botany, 2002, 53(375):1689–1698.

[23]王利紅, 徐芳森, 王運(yùn)華. 硼鉬鋅配合對(duì)甘藍(lán)型油菜產(chǎn)量和品質(zhì)的影響[J]. 植物營(yíng)養(yǎng)與肥料學(xué)報(bào), 2007, 13(2): 318–323.Wang L H, Xu F S, Wang Y H. Effects of B, Mo, Zn interaction on the seed yield and quality of Brassica napus[J]. Plant Nutrition and Fertilizer Science, 2007, 13(2): 318–323.

[24]梁新華, 張風(fēng)俠, 王俊, 何文壽. 鉬、硼、錳和鋅對(duì)人工種植烏拉爾甘草品質(zhì)的影響[J]. 植物營(yíng)養(yǎng)與肥料學(xué)報(bào), 2011, 17(6): 1487–1494.Liang X H, Zhang F X, Wang J, He W S. Effects of molybdenum,boron, manganese and zinc on quality of cultured Glycyrrhiza uralensis Fisch[J]. Plant Nutrition and Fertilizer Science, 2011, 17(6):1487–1494.

[25]湯璐, 林江輝, 閆廣軒, 等. 銅、鋅、硒對(duì)藥用菊花主要有效成分和花中硒含量的影響[J]. 植物營(yíng)養(yǎng)與肥料學(xué)報(bào), 2009, 15(6):1475–1480.Tang L, Lin J H, Yan G X, et al. Effects of Cu, Zn and Se on contents of total flavonlid, chlorogenic acid and Se in the flower of Chrysanthemum morifolium Ramat[J]. Plant Nutrition and Fertilizer Science, 2009, 15(6): 1475–1480.

[26]謝建軍, 張佳杰, 姜俊基, 等. 施用微肥對(duì)甘草生長(zhǎng)動(dòng)態(tài)和產(chǎn)量的影響[J]. 廣東農(nóng)業(yè)科學(xué), 2012, 39(9): 62–66.Xie J J, Zhang J J, Jiang J J, et al. Influence of microfertilizers on yield and growth and development of Glycyrrhiza uralensis[J].Guangdong Agricultural Sciences, 2012, 39(9): 62–66.

[27]楊保川, 王昌全, 項(xiàng)虹艷, 盧益武. 葉面噴施鉬錳對(duì)萵筍生長(zhǎng)發(fā)育和硝酸鹽含量的影響[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào), 2006, 25(s1): 96–99.Yang B C, Wang C Q, Xiang H Y, Lu Y W. Effects of foliar spraying molybdenum, manganese on growth and nitrate content of lettuce[J].Journal of Agro-Environment Science, 2006, 25(s1): 96–99.

[28]張紀(jì)利, 李余湘, 羅紅香, 等. 施鉬對(duì)煙草葉綠素含量、光合速率、產(chǎn)量及品質(zhì)的影響[J]. 中國(guó)煙草科學(xué), 2011, 32(2): 24–28.Zhang J L, Li Y X, Luo H X, et al. Effects of molybdenum application on chlorophyll contents, photosynthetic rates, yield and quality of tobacco leaves[J]. Chinese Tobacco Science, 2011, 32(2):24–28.

[29]竇明明, 石峰, 馬留輝, 等. 微肥配施對(duì)澤瀉產(chǎn)量的影響[J]. 中藥材,2017, 40(1): 7–11.Dou M M, Shi F, Ma L H, et al. Effects of applying micro-element fertilizer on the yield of Alisma orientalis[J]. Journal of Chinese Medicinal Materials, 2017, 40(1): 7–11.

[30]夏金嬋, 何奕騉. 植物對(duì)硼元素的吸收轉(zhuǎn)運(yùn)機(jī)制[J]. 中國(guó)生物化學(xué)與分子生物學(xué)報(bào), 2009, 25(8): 702–707.Xia J C, He Y K. Mechanisms of boron transport in plants[J].Chinese Journal of Biochemistry and Molecular Biology, 2009,25(8): 702–707.

[31]張亞琴, 楊正明, 石峰, 等. 葉面噴施微肥對(duì)川白芷主要有效成分含量的影響[J]. 應(yīng)用生態(tài)學(xué)報(bào), 2017, 28(11): 3505–3514.Zhang Y Q, Yang Z M, Shi F, et al. Effects of micro－fertilizers foliar spray on the content of main effective components of Angelica dahurica [J]. Chinese Journal of Applied Ecology, 2017, 28(11):3505–3514.

[32]鐘世紅. 川白芷規(guī)范化種植質(zhì)量評(píng)價(jià)的研究[D]. 成都: 成都中醫(yī)藥大學(xué)碩士學(xué)位論文, 2005.Zhong S H. Study on the quality evaluation of the standardized cultivation of Angelica dahurica [D]. Chengdu: MS Thesis of Chengdu University of Traditional Chinese Medicine, 2005.