表1 各处理施肥量
Tab.1 The amount of fertilizer used in each treatment kg/hm2
处理Treatment施肥量 Amount of fertilizerNP2O5K2ONPK90120100NP901200NK900100PK0120100CK000
芸豆(Phaseolus vulgaris L.)是我国一种主要的杂粮作物，生育期短，耐旱耐瘠，抗逆性强，籽粒营养丰富，在全国高寒冷凉区广泛种植，其中，山西省是我国芸豆主产地之一。红芸豆属于食用豆中普通菜豆矮生直立型的一个品种，因其色泽艳丽，营养价值和经济价值较高，在国内外市场受到欢迎[1]，近年来种植面积逐年扩大。
我国对芸豆的研究起步较晚，始于1978年，而且科研力量薄弱，育成品种较少，栽培手段落后，田间管理粗放，造成芸豆单产水平较低，极大地影响了芸豆的规模化生产[1-3]。
近年来对红芸豆的研究主要集中在施肥、种植密度、光合和水分方面。已有研究表明，施肥对红芸豆产量及产量构成因素具有促进作用[4-12]，施用氮、磷肥可以提高不同品种芸豆叶绿素含量、氨基酸含量及组分[13-16]，不同的种植密度对红芸豆产量和品质也形成一定的影响[16-18]。这些研究主要集中在施肥及氮肥对芸豆产量和农艺性状的影响上，而关于施肥对红芸豆的养分吸收、干物质累积、分配特征及对产量构成因子的影响以及之间关系的研究未见报道。但在大豆上有一些关于施肥对大豆品质、干物质积累、源库调节效应等方面的研究[20-26]，杂豆高效生产关键技术研究与示范课题组在试验设计上参考了这些试验的思路。
根据红芸豆产区生产特点，课题试验方案以土壤测定结果结合农民习惯施肥确定施肥量，设置氮、磷、钾配施及缺素试验，系统研究不同氮磷钾配施对红芸豆养分吸收、干物质积累及产量相关因子的影响，以期为红芸豆高效生产及适量施肥提供理论依据。
本试验于2014-2015年在山西省农业科学院旱农农业研究中心的河村试验基地(N38°05′，E112°90′)进行。试验地位于山西省中部，海拔1 270 m，无霜期140 d左右，年均降水量450 mm，年均蒸发量1 995 mm，年均气温6～7 ℃，≥10 ℃活动积温2 600 ℃，为典型的半干旱区。试验地土壤属于褐土性土，土壤pH值8.02，有机质含量为16.50 g/kg，全N为1.09 g/kg，碱解N为55.12 mg/kg，有效P为14.63 mg/kg，速效K为161.09 mg/kg。前茬作物为玉米。
供试作物为红芸豆，属于普通菜豆(Phaseolus vulgaris Linn)矮生直立型品种，供试品种为英国红。施用的氮肥为尿素，N含量为46%；磷肥为粒状过磷酸钙，P2O5含量为12%；钾肥为硫酸钾，K2O含量为45%。
本试验设5个处理，分别为NPK、NP、NK、PK、CK，各处理具体施肥量如表1所示。3次重复，随机区组排列，小区面积28 m2(7 m×4 m)。0.5 m等行距种植，穴距0.2 m，每穴2株。试验在自然降水条件下进行，无补充灌溉。2015年试验设计、田间布置与2014年相同。
表1 各处理施肥量
Tab.1 The amount of fertilizer used in each treatment kg/hm2
处理Treatment施肥量 Amount of fertilizerNP2O5K2ONPK90120100NP901200NK900100PK0120100CK000
1.4.1 测定日期 在红芸豆的生长发育过程中对其进行采样，生育期的划分及取样、测量日期如表2所示。
表2 红芸豆的生育期划分及取样、测量日期
Tab.2 Red kidney bean growth stage and the sampling and measurement date 月-日
年份Year播种Seeding出苗Period of emergenceV2R1R3R6R8201405-1805-2806-1407-0107-1908-0808-26201505-1505-2506-1107-0207-1908-1108-26
注：V2.苗期，第2片未完全展开的3出复叶出现在第4节上；R1.初花期，植株上第1朵花开放；R3.盛花期，植株上有1个荚达到最大长度；R6.结荚期，植株上50%的豆荚中豆粒完全发育；R8.成熟期，植株上80%的豆荚达到了成熟颜色。表5-12同。
Notes: V2.Seedling stage, the trifoliolate leaf of second not fully expand appears in fourth section; R1.Early flowering stage, the first flower on the plant opens; R3.Full bloom stage, there is 1 pod on the plant to reach the maximum length; R6.Pod bearing stage, 50% of the pods on the plant were fully developed; R8.Maturity stage, 80% of the pod on the plant reached its mature color. The same as Tab.5-12.
1.4.2 红芸豆植株的干物质积累、分配测定及计算方法 分别在红芸豆V2、R1、R3、R6和R8共5个时期选有代表性的植株4～6株，用自来水小心冲洗干净，然后迅速将采集的红芸豆植株按照器官分为茎、叶、荚皮、籽粒等，将分好的各部分器官装入纸袋，在105 ℃烘箱中杀青30 min，于80 ℃恒温条件下烘至恒质量，称取各器官的干质量，并计算干物质量，进一步推算各时期各部位的干物质积累和分配转移状况。
1.4.3 田间测定方法 挑选长势均匀连续的10株红芸豆，各处理分别挂牌标记。在红芸豆成熟后，即将收获时，对各处理小区内挂牌的10株红芸豆进行考株，分别考察单株有效荚数、单株粒数、百粒质量等。各小区单收单打，经去杂和晾晒计产。
采用Excel 2003、SPSS 14.0等进行数据处理和统计分析。
从表3可以看出，与CK相比，各施肥处理的产量均有提高，并与对照间差异显著；其中，2014年缺氮、缺磷处理的产量与全施肥处理差异显著，缺钾处理的产量与全施肥处理差异不显著，从产量构成因子看，差异主要体现在百粒质量上，但缺磷处理与全施肥处理之间的百粒质量差异不显著，而单株有效荚数和每荚粒数差异不显著；2015年3个缺素处理均与全施肥处理的产量差异显著，从产量构成因子看，其差异表现与2014年基本吻合；2014年各处理产量比CK增产11.48%～25.15%，其中氮、磷、钾配施比缺施氮、磷、钾肥处理分别增产12.27%，4.20%，3.34%；2015年各处理产量比CK增产21.14%～37.67%。其中氮、磷、钾配施比缺施氮、磷、钾肥处理分别增产13.65%，10.55%，5.39%，由此可见，说明养分是红芸豆产量提高的主要限制因子。从产量构成因子来看，各施肥处理对产量的影响主要表现在对百粒质量的影响上，除2014年缺氮处理与对照差异不显著外，其他施肥处理的百粒质量与CK均差异显著；而对单株有效荚数和每荚粒数则影响较小，除2014年缺钾处理与对照差异显著外，其他施肥处理均与对照差异不显著(表3)。
经过对2年产量及构成因子进行相关性分析可知，2年的趋势高度一致，对产量影响的重要性依次为百粒质量>每荚粒数>单株有效荚数，但只有百粒质量与产量的相关性达到显著水平(P<0.05)，2014，2015年的相关系数分别为0.88和0.90(表4)。
表3 红芸豆各处理的产量及构成因子
Tab.3 Red kidney bean yield and constituent factor
年份Year处理Treatment单株有效荚数Pod number on single plant每荚粒数The number of seeds per pod百粒质量/g100-grain weight 平均产量/(kg/hm2)Average yield比对照增产/%Gain in yield than contrast2014NPK11.03±0.27a3.11±0.16ab52.20±2.46a2 538.09±66.28a25.15NP9.33±0.07a3.65±0.11a47.22±2.34bc2 455.95±71.15ab21.10NK10.10±1.00a2.99±0.23ab51.87±1.66ab2 435.71±29.03b20.11PK9.30±0.15a2.98±0.20ab45.28±0.63c2 260.71±43.06c11.48CK9.51±1.14a2.70±0.37b42.74±2.02c2 027.97±83.85d0.002015NPK8.50±0.69a3.70±0.12a53.78±1.95a2 419.05±45.42a37.67NP7.90±0.15ab3.93±0.15a49.28±1.53c2 295.24±12.60b30.62NK7.30±0.38ab3.80±0.10a52.37±0.60ab2 188.10±41.72c24.53PK6.90±0.60b3.77±0.19a49.92±1.48bc2 128.57±32.99c21.14CK7.80±0.52ab3.63±0.23a44.12±0.63d1 757.14±68.14d0.00
表4 产量及构成因子相关性分析
Tab.4 Correlation analysis of yield and constituent factors
年份Year相关系数Correlation coefficient单株有效荚数Pod number onsingle plant每荚粒数The number of seeds per pod百粒质量100-grain weight平均产量Average yield2014单株有效荚数-0.090.790.58每荚粒数-0.090.290.67百粒质量0.790.290.88*平均产量0.580.670.88*2015单株有效荚数-0.170.170.36每荚粒数-0.170.310.54百粒质量0.170.310.90*平均产量 0.360.540.90*
注：*.P<0.05；**.P<0.01。表11-12同。
Note:*.P<0.05；**.P<0.01. The same as Tab.11-12。
2.2.1 红芸豆各处理的养分含量 从红芸豆植株样本的分析测试结果来看，各处理氮含量变化各不相同，在V2-R1阶段，与CK相比，除NP处理的含量略有降低外，其他4个处理均升高；在R3-R8阶段，除NPK处理是一直降低外，其他4个处理均为先降低后升高；氮含量总的趋势呈高-低-高的变化。各处理磷含量随着生育期的推进，除NK在R6和CK在R8略有升高外，其他4个处理均呈现降低趋势。钾含量在V2-R1阶段，除NPK处理略有升高外，其他处理均在降低，在R3-R8阶段的含量与前期相比则降低，但2个阶段内差异不显著(表5)。
2.2.2 红芸豆各器官的养分含量 从表6可以看出，根系含氮量在不同生育阶段差异显著，呈U型趋势，从V2到R6逐渐降低，R8又有回升；茎含氮量呈倒S型趋势，从V2到R1先升高，然后逐渐降低，R8又有回升，但V2和R1这2个时期差异不显著；叶片含氮量的趋势与根系一样，呈U型趋势，从V2到R6逐渐降低，R8又有回升。这3个器官氮含量在R8的回升可能与氮在成熟期的回流有关。豆荚和荚皮含氮量呈降低趋势，籽粒含氮量呈升高趋势，这是干物质转移的表现。
从植株各器官含磷量数据来看，根系含磷量除在V2与R1和R6有显著差异外，R1-R8阶段差异并不显著，整个生育期变化比较小；茎含磷量在V2-R6阶段逐渐降低，R8略有回升，但与前期差异不显著；叶片含磷量在V2和R1阶段比较高，R3-R8则较低，且这2个生育区间差异显著，但区间内则差异不显著；豆荚、荚皮和籽粒的含磷量均呈降低趋势。
从植株各器官含钾量数据来看，根系含钾量呈明显S型变化；茎含钾量呈抛物线变化，但R3-R8含量明显降低；叶片含钾量呈显著下降趋势，R8虽略有回升，但与R6间差异不显著；豆荚含钾量呈降低趋势；荚皮含钾量无显著差异，籽粒则降低，差异显著。
从整个植株数据来看，各器官中根、叶片、豆荚及籽粒中氮>钾>磷，茎和荚皮中钾>氮>磷；植株各器官在不同生育期的氮、磷、钾含量变化较大，从整个生育期各器官的氮、磷、钾含量来看，各器官中含氮量依次为豆荚>叶>籽粒>茎>根>荚皮，含磷量依次为豆荚>籽粒>叶>茎>根>荚皮，含钾量依次为豆荚>荚皮>茎>籽粒>叶>根(表6)。
表5 各处理不同生育期的氮磷钾含量
Tab.5 NPK content in different growth stages of each treatment g/kg
养分Nutrient处理TreatmentV2R1R3R6R8NNPK29.75±0.58a32.34±0.98a25.22±0.49b16.29±2.07c15.62±1.32cNP29.38±0.53a27.50±0.93a21.26±0.43b13.93±0.43d17.16±1.23cNK29.86±0.16a30.85±1.68a20.92±1.21b14.79±1.29c15.74±0.89cPK29.00±0.14a29.11±1.45a19.70±0.94b12.09±0.48d15.99±0.36cCK29.43±0.57a29.63±0.65a21.53±0.99b11.70±0.34d16.36±0.83cPNPK5.90±0.17a5.88±0.27a5.42±0.26ab4.76±0.10b3.58±0.45cNP5.86±0.18a5.34±0.27ab4.67±0.25bc4.18±0.45c3.88±0.12cNK6.00±0.52a5.55±0.22ab4.48±0.24bc4.53±0.17bc3.50±0.48cPK6.36±0.41a5.84±0.14a4.50±0.02b3.72±0.21b4.32±0.27bCK5.55±0.04a5.38±0.16ab4.66±0.27bc3.65±0.27d4.13±0.38cdKNPK22.17±2.91a24.83±2.44a15.08±0.02b13.85±1.25b11.97±1.32bNP26.01±2.19a23.25±2.88a13.76±0.73b11.61±1.06b12.36±0.43bNK23.48±1.71a23.40±3.57a14.37±1.31b12.84±1.84b11.41±0.78bPK24.69±0.88a20.86±1.27b15.96±1.19c12.78±1.28c13.54±0.58cCK28.60±3.87a24.96±1.39a14.75±0.24b11.07±0.46b12.26±0.95b
表6 各器官不同生育期的氮磷钾含量
Tab.6 NPK content in different growth stages of each organs g/kg
养分Nutrient器官Plant organV2R1R3R6R8N根20.23±0.21a15.17±0.12b13.52±0.24c8.24±0.27e10.71±0.69d茎25.50±0.91a26.15±1.07a17.81±1.80b7.74±0.39d10.59±0.47c叶42.73±0.52a41.88±1.08a26.32±0.88b19.16±1.74c24.47±0.78b豆荚36.34±2.41a29.25±1.56b荚皮9.14±0.99a4.94±0.21b籽粒24.52±0.67a30.17±3.43aP根4.73±0.66a3.03±0.13b4.38±0.61ab3.03±0.13b3.82±0.21ab茎5.37±0.28a5.11±0.15ab3.91±0.19bc3.15±0.38c3.33±0.77c叶7.70±0.10a6.77±0.36a4.26±0.49b4.83±0.31b4.38±0.22b豆荚7.48±0.48a6.42±0.41b荚皮2.56±0.17a1.76±0.04a籽粒6.13±0.27a6.11±0.59aK根21.63±1.91a15.92±1.67b8.16±0.65c15.92±1.67b6.91±0.33c茎33.14±2.58a33.87±5.96a16.09±1.83b9.74±0.64b8.34±0.49b叶20.21±0.15a16.69±0.92b11.42±0.53c9.07±0.95d9.19±0.51d豆荚27.36±0.36a23.47±0.50b荚皮22.30±1.26a22.32±3.70a籽粒15.97±0.64a14.79±0.37b
2.3.1 红芸豆各器官氮、磷、钾养分积累和分配特征 从整个生育期看，红芸豆植株总的氮、磷、钾积累量均呈增加趋势，从V2-R3是营养器官的快速积累期，此阶段叶片和茎秆是养分的分配中心；到R3后，生殖器官逐渐成为养分的分配中心，到R8时，籽粒的氮、磷、钾积累量所占比例分别达到植株总积累量的71.73%，64.94%，46.07%；养分积累量为氮>钾>磷。
从氮的积累量来看，根、茎、叶中氮积累量在生育前期呈快速升高趋势，生育后期开始下降，但根和茎在成熟期均有升高；而生殖器官中氮积累量在生育后期主要转向籽粒；成熟期各器官氮累积量依次为籽粒>茎>叶>荚皮>根，其所占总量的比例分别为71.73%，10.55%，10.09%，6.23%，1.40%。
从磷的积累量来看，所有器官中磷积累量的趋势和氮积累量的趋势完全相同，只是叶片和荚皮的磷积累量前后排序颠倒；成熟期各器官磷积累量依次为籽粒>茎>荚皮>叶>根，其所占总量的比例分别为64.94%，14.83%，9.92%，8.07%，2.24%。
从钾的积累量来看，根在V2-R6呈抛物线趋势，到V8有所上升，茎和叶中钾积累量均呈抛物线趋势，只是根系的峰值在R1，茎秆和叶片的峰值在R3；而生殖器官中荚皮的钾积累量和氮、磷趋势不同，呈现增加的趋势；成熟期各器官钾累积量依次为籽粒>荚皮>茎>叶>根，其所占总量的比例分别为46.07%，36.89%，10.89%，4.96%，1.19%(表7)。
表7 各器官不同生育期氮磷钾积累量及占比
Tab.7 NPK accumulation and percentage of total in different growth stages of each organs
养分Nutrient积累量/占总量比例Accumulation amount/Percentage of total器官Plant organV2R1R3R6R8N积累量根0.72±0.02d1.84±0.15bc2.69±0.38a1.53±0.15c2.25±0.13ab/(kg/hm2)茎1.43±0.09d13.13±1.77bc24.80±3.58a11.61±1.94c16.93±1.77b叶7.54±0.42c36.51±4.00a36.47±5.74a25.76±5.92b16.18±5.28bc豆荚1.20±0.16b34.01±1.73a荚皮17.75±1.68a10.00±0.51a籽粒66.74±0.91b115.11±4.10a占总量根7.433.492.751.241.40比例/%茎14.7224.9225.329.4110.55叶77.8569.3037.2320.8710.09豆荚2.2934.70荚皮14.396.23籽粒54.0971.73P积累量根0.17±0.00c0.37±0.03b0.87±0.12a0.77±0.07a0.80±0.05a/(kg/hm2)茎0.30±0.02c2.57±0.35b5.44±0.78a4.72±0.79a5.32±0.56a叶1.36±0.08b5.89±0.65a5.92±0.93a6.49±1.49a2.90±0.95b豆荚0.25±0.03b7.46±0.38a荚皮4.97±0.47a3.56±0.18a籽粒16.69±0.23b23.31±0.83a占总量比例/%根9.134.044.422.302.24茎16.4228.2827.6414.0414.83叶74.4564.9330.0419.298.07豆荚2.7537.90荚皮14.789.92籽粒49.5964.94K积累量根0.77±0.03c1.93±0.15a1.62±0.23ab0.94±0.09c1.45±0.08b/(kg/hm2)茎1.85±0.11c17.01±2.29b22.41±3.24a14.61±2.44b13.33±1.39b叶3.56±0.20b14.55±1.60a15.83±2.49a12.19±2.80a6.08±1.98b豆荚0.90±0.12b27.29±1.39a荚皮43.31±4.09a45.18±2.29a籽粒43.47±0.59b56.43±2.01a占总量根12.495.622.410.821.19比例/%茎29.9449.4533.3712.7510.89叶57.5742.3023.5710.654.96豆荚2.6340.65荚皮37.8236.89籽粒37.9646.07
2.3.2 红芸豆植株各生育阶段氮、磷、钾养分吸收、积累和分配特征 随着生育期的推进，红芸豆植株氮、磷、钾的积累量在逐渐累积增加，到收获时达到最大，此时氮、磷、钾的积累量分别达到161.15，38.27，126.7 kg/hm2，其积累量之间的比例为4.21∶1.00∶3.31；从红芸豆整个生育期各阶段的氮、磷、钾积累量来看，其积累量呈为抛物线形式，这是由于生育前期植株较小，其对养分的吸收量也相对较小，随着植株生长，对养分的需求也在增加，中期进入养分需要旺盛期，后期进入生殖生长阶段后，养分主要向籽粒转移；但氮、磷、钾三者的积累量的顶点不同，氮出现在R1-R3，而磷、钾出现在R3-R6(表8)。
表8 各生育阶段氮磷钾积累量、占比及吸收速率
Tab.8 NPK accumulation and percentage of total and absorption rate in different growth stages
养分Nutrient积累量/占总量比例Accumulation amount/Percentage of totalV2V2-R1R1-R3R3-R6R6-R8N积累量/(kg/hm2) 8.92±0.06d36.64±0.59c56.74±1.52a48.88±2.53b9.97±0.41d占总量比例/%5.54±0.14d22.76±0.34c35.22±0.16a30.29±0.64b6.18±0.13d吸收速率/(kg/(hm2·d)) 0.52±0.00d2.15±0.03c3.15±0.09a2.44±0.13b0.56±0.02dP积累量/(kg/hm2)1.80±0.05d6.86±0.24c12.40±0.31b14.82±0.63a2.39±0.17d占总量比例/% 4.70±0.17d17.92±0.20c32.43±1.08b38.70±0.80a6.25±0.41d吸收速率/(kg/(hm2·d)) 0.10±0.00c0.40±0.01b0.69±0.02a0.74±0.03a0.13±0.01cK积累量/(kg/hm2)7.56±0.44d28.76±2.63c38.61±1.36b44.18±1.03a7.59±0.49d占总量比例/% 5.97±0.25d22.62±1.29c30.50±0.65b34.94±1.04a5.98±0.22d吸收速率/(kg/(hm2·d)) 0.45±0.03c1.69±0.16b2.15±0.08a2.21±0.05a0.42±0.03c
从红芸豆整个生育期单株干物质积累量可以看出，V2-R6干物质累积量差异显著，R6-R8则无显著差异，呈生育前期快速积累，生育后期逐渐平缓的变化趋势；干物质积累速率呈抛物线形状，R3-R6的20 d达到0.952 g/d，是全生育期的干物质积累高峰。
从红芸豆各器官干物质积累量来看，根系和叶片在生育前期快速增长，积累速率到R1达到高峰，在V2-R1的17 d中分别为0.027，0.218 g/d，后期趋向缓慢；茎秆的趋势和根系、叶片相同，只是积累速率到R3达到高峰，在R1-R3的18 d中达到为0.263 g/d；从R3以后，生殖器官的干物质积累成为主导，豆荚在R1-R3的18 d中达到0.336 g/d，荚皮、籽粒的最高积累速率在R3-R6的20 d中达到0.520，0.728 g/d；由此可见，生育前期营养器官是积累分配中心，在R8干物质占总量比例达到72.22%，生育后期生殖器官是积累分配中心，在R8干物质占总量比例达到70.26%，这与红芸豆养分吸收是一致的(表9-10)。
表9 各器官在不同生育期单株的干物质积累量及占比
Tab.9 Dry matter accumulation amount and percentage of total in different growth stages of each organs
干物质积累量/占总量比例Dry matter accumulation amount/Percentage of total器官Plant organV2R1R3R6R8干物质积累量/g根0.19±0.01c0.65±0.05b1.06±0.15a0.99±0.09a1.13±0.07aDry matter 茎0.30±0.02c2.69±0.36b7.43±1.07a8.02±1.33a8.58±0.88aaccumulation 叶0.95±0.05c4.66±0.51b7.38±1.15a7.19±1.66a3.53±1.16bamount豆荚0.18±0.02b6.22±0.30a荚皮10.39±0.95a10.85±0.56a籽粒14.55±0.23b20.43±0.73a植株1.44±0.07d8.18±0.95c22.10±2.60b41.14±3.79a44.53±2.56a占总量比例/%根13.437.944.752.412.53Percentage茎20.8332.8733.1819.4919.28of total叶65.7456.9934.2917.497.94豆荚2.2027.78荚皮25.2624.37籽粒35.3645.89植100.00100.00100.00100.00100.00
表10 红芸豆单株干物质积累量及积累速率
Tab.10 Dry matter accumulation amount and accumulation rate in plants
器官Plant organ干物质积累量/累积速率Dry matter accumulation amount/Accumulation rateV2V2-R1R1-R3R3-R6R6-R8根干物质积累量/g 0.190.460.41-0.070.14Root累积速率/(g/d) 0.0110.0270.023-0.0040.008茎干物质积累量/g 0.302.394.740.590.56Stem累积速率/(g/d)0.0180.1410.2630.0300.031叶干物质积累量/g 0.953.712.72-0.19-3.66Leaf累积速率/(g/d) 0.0560.2180.151-0.010-0.203豆荚干物质积累量/g 0.186.04Pod累积速率/(g/d) 0.0110.336荚皮干物质积累量/g 10.390.46Pod shell累积速率/(g/d) 0.5200.026籽粒干物质积累量/g 14.555.88Grain累积速率/(g/d) 0.7280.327植株干物质积累量/g 1.446.7413.9219.043.39Plant累积速率/(g/d) 0.0850.3960.7730.9520.188
2.5.1 各生育期养分积累量与干物质积累量的关系 通过对红芸豆各生育期的养分积累量与干物质积累量之间关系进行相关性分析可以看出，随着生育期的推进，氮积累量与干物质积累量之间的相关性呈逐渐增加趋势，到R6和R8时其相关性极显著；磷积累量与干物质积累量之间的相关性呈倒S型，到R6时其相关性极显著；钾积累量与干物质积累量之间的相关性与磷含量相似，到R6时其相关性达到最大，但没有达到显著水平。由此可见，施肥对红芸豆植株的干物质积累量起到至关重要的作用(表11)。
表11 各生育期养分积累量与干物质积累量的相关性分析
Tab.11 Correlation analysis of nutrient accumulation amount and dry matter accumulation
amount at different growth stages
相关系数Correlation coefficientV2R1R3R6R8氮积累量 N accumulation amount0.260.360.680.95**0.96**磷积累量 P accumulation amount0.560.380.690.93**0.28钾积累量 K accumulation amount0.360.070.440.790.19
2.5.2 各生育期干物质积累量与产量及构成因子的关系 通过对红芸豆各生育期干物质积累量与产量及构成因子之间关系进行相关性分析可以看出，单株有效荚数和生育前期的相关性较小，和生育后期的相关性较大，但均没有达到显著水平；每荚粒数与初花期的干物质积累量有极显著的相关性，说明此时的干物质积累量会极大地影响到红芸豆每荚粒数的多少；百粒质量的多少与结荚期的干物质积累量呈极显著的相关性，此生育阶段的干物质积累量决定了百粒质量的多少；产量和干物质积累量在盛花期达到显著水平，到结荚期达到极显著水平。可以看出，生殖生长阶段的干物质积累量是影响产量结果的关键因素(表12)。
表12 各生育期干物质积累量与产量及构成因子的相关性分析
Tab.12 Correlation analysis of dry matter accumulation amount and yield and constituent factor at different growth stages
相关系数Correlation coefficientV2R1R3R6R8单株有效荚数Pod number on single plant0.240.010.550.710.74每荚粒数The number of seeds per pod-0.630.95**0.760.54-0.67百粒质量100-grain weight0.000.390.660.94**0.41平均产量Average yield-0.170.690.86*0.98**0.06
氮、磷、钾配施对红芸豆具有显著的增产作用，产量构成因子中百粒质量对产量的形成影响最大；氮、磷、钾配施比缺施氮、磷、钾肥处理2014年增产分别为12.27%，4.20%，3.34%，2015年增产分别为13.65%，10.55%，5.39%，影响产量的养分限制因子氮>磷>钾,与畅建武等[9]的试验结果相同，而韩彦龙等[5]的试验结果是氮>钾>磷。
各器官中含氮量依次为豆荚>叶>籽粒>茎>根>荚皮，含磷量依次为豆荚>籽粒>叶>茎>根>荚皮，含钾量依次为豆荚>荚皮>茎>籽粒>叶>根；根、茎、叶这3个器官氮含量在R8的回升可能与氮在成熟期的回流有关。
红芸豆植株在整个生育期氮、磷、钾积累量均呈增加趋势，积累量分别达到161.15，38.27，126.7 kg/hm2，比例为4.21∶1.00∶3.31；氮积累量的顶点出现在初花到盛花期，而磷、钾出现在盛花到结荚期；生育前期叶片和茎秆是养分的分配中心；生育后期籽粒是养分的分配中心，氮、磷、钾积累量占总量的比例分别为71.73%，64.94%，46.07%；养分积累量为氮>钾>磷。
红芸豆干物质积累趋势与养分吸收趋势基本一致，生育前期营养器官是干物质积累分配中心，其干物质占总量比例达到72.22%，生育后期生殖器官是干物质积累分配中心，其干物质占总量比例达到70.26%；整个生育期干物质积累速率呈抛物线形状，盛花到结荚期达到积累高峰，单株积累量为0.952 g/d；各器官积累速率的高峰不同，根系和叶片在初花期，茎秆和豆荚在盛花期，荚皮和籽粒在结荚期。韩彦龙等[5]的也有相似结论，只是在各器官的积累速率上有所不同。
红芸豆干物质积累量与养分积累量具有相关性，氮积累量与干物质积累量之间的相关系数在结荚期到成熟期分别为0.95和0.96，且呈极显著水平(P<0.01)；磷积累量与干物质积累量之间的相关系数在结荚期为0.93(P<0.01)，且呈极显著水平；钾积累量与干物质积累量之间的相关性没有达到显著水平；红芸豆干物质积累量与产量构成因子具有相关性，每荚粒数与干物质积累量的相关系数在初花期为0.95，且达到极显著水平(P<0.01)，百粒质量与干物质积累量在结荚期的相关系数为0.94，且呈极显著水平(P<0.01)，产量与干物质积累量的相关系数在盛花期为0.86，达到显著水平(P<0.05)，在结荚期为0.98，达到极显著水平(P<0.01)，生殖生长阶段的干物质积累量是影响产量结果的关键因素。
氮磷钾配施是一个大田肥料试验，由于降水、温度和农田小气候的不可控性，受到诸多环境因素影响，使得试验结果进行纵向和横向对比时重现性差。因此，仍需要进行多年定点的连续试验才有可能得出比较客观的结论。
[1] 林汝法. 中国小杂粮[M]. 北京: 中国农业科学技术出版社, 2002.
[2] 龙静宜. 食用豆类作物[M]. 北京: 科学出版社, 1989.
[3] 郑卓杰. 中国食用豆类学[M]. 北京: 中国农业出版社, 1997.
[4] 宋谨同, 赵宏伟. 氮肥用量对芸豆淀粉含量及产量的影响[J]. 作物杂志, 2010(5): 76-78.
[5] 韩彦龙, 晋凡生, 郑普山, 等. 红芸豆养分限制因子及养分吸收、积累和分配特征研究[J]. 中国生态农业学报, 2016,24(7):902-909.
[6] 乔秀平. 不同施氮水平对芸豆产量的影响[J]. 山西农业科学, 2014, 42(7): 694-696,703.
[7] 宋谨同, 赵宏伟, 杨 亮, 等. 氮肥用量对芸豆氮肥利用率和产量影响的研究[J]. 农业现代化研究, 2013, 34(6): 749-753.
[8] 曾玲玲, 崔秀辉, 李清泉, 等. 氮磷钾配施对芸豆产量的效应研究[J]. 黑龙江农业科学, 2013(2): 39-43.
[9] 畅建武, 郝晓鹏, 王 燕, 等. 红芸豆氮磷钾肥效试验研究[J]. 中国农学通报, 2015, 31(15): 108-113.
[10] 高运青, 徐东旭, 尚启兵, 等. 华北高寒区施肥量对芸豆产量和经济效益的影响[J]. 河北农业科学, 2012, 16(8): 28-30.
[11] 韩彦龙, 晋凡生, 李晓平, 等. 干旱冷凉区红芸豆施肥效应与养分吸收研究[J]. 中国农学通报, 2016,32(15):73-78.
[12] 李 洁, 晋凡生, 李晓平, 等. 施肥和覆膜对旱地垄沟种植红芸豆耗水特征及产量效应的研究[J]. 中国农学通报, 2018,34(13):19-23.
[13] 孟宪欣. 磷肥对芸豆几个品质性状的影响[J]. 中国农学通报, 2010, 26(24): 183-187.
[14] 杨 亮, 赵宏伟, 宋谨同. 氮肥用量对芸豆叶绿素含量和子粒营养品质影响的研究[J]. 作物杂志, 2013(1): 81-87.
[15] 晋凡生，韩彦龙.氮磷钾配比对红芸豆光合特性及产量相关性状的影响[J]. 山西农业科学，2017，45(8):1295-1230.
[16] 华劲松, 戴红燕, 夏明忠. 不同光照强度对芸豆光合特性及产量性状的影响[J]. 西北农业学报, 2009, 18(2): 136-140.
[17] 杨广东, 张亚芝, 魏淑红. 种植密度对粒用芸豆的产量及商品品质的影响[J]. 黑龙江农业科学, 2012(12): 24-26.
[18] 张晓艳. 种植密度对芸豆产量和品质的影响研究[D]. 太谷: 山西农业大学, 2004.
[19] 程丽娟, 杨锦忠, 杜天庆, 等. 种植密度对粒用芸豆品质的影响研究[J]. 山西农业大学学报：自然科学版, 2009，29(2): 129-131,142.
[20] 闫春娟, 韩晓增, 王树起, 等. 钾对大豆干物质积累、产量及品质的影响[J]. 大豆科学, 2008, 27(1): 113-117.
[21] 徐本生, 籍玉尘, 杨建堂. 夏大豆的干物质积累和氮磷钾吸收分配动态的研究[J]. 大豆科学, 1989, 8(1): 47-53.
[22] 王景伟, 金喜军, 杜文言, 等. 干旱胁迫对芸豆籽粒干物质积累的影响及动态模型的建立[J]. 干旱地区农业研究, 2014, 32(2): 147-150.
[23] 李春杰, 王建国, 许艳丽, 等. 钾对大豆产量及品质的影响[J]. 农业系统科学与综合研究, 2005, 21(2): 154-160.
[24] 刘建国, 帕尼古丽, 董志新, 等. 单位叶面积负荷量对大豆源库调节效应的研究[J]. 石河子大学学报：自然科学版, 2003, 7(4): 260-263.
[25] 黄中文, 王 伟, 徐新娟, 等. 大豆动态株高及其生长速率与产量的相关分析[J]. 河南科技学院学报, 2010, 38(2): 16-19.
[26] 许海涛, 王友华, 许 波, 等. 氮磷钾优化施肥对高蛋白大豆籽粒及蛋白质产量的影响[J]. 土壤通报, 2008, 39(1): 195-196.