Effects of applying compound fertilizer on yield, nutrient absorption and oligosaccharide accumulation of Morinda officinalis
-
摘要:目的
研究复合肥施用量对巴戟天产量、养分吸收和主要有效成分寡糖累积量的影响，为巴戟天的施肥管理提供理论依据。
方法以巴戟天为研究对象，通过盆栽试验，设置每株施用复合肥0、2、4、6、8和10 g (CK、T1、T2、T3、T4和T5)共6个处理，测定巴戟天在不同施肥处理下的产量、养分含量和寡糖含量，分析各指标间的相关性，并运用主成分分析方法进行综合评价。
结果施用复合肥可以显著提高巴戟天产量和N、P、K含量，并有利于提高寡糖的累积量。其中，单株施肥量为6 g处理的巴戟天产量和寡糖总累积量最高，与CK相比，分别提升114.47%和75.36%。施肥在一定程度上降低了巴戟天的寡糖含量，与CK相比，降幅最高达35.65%。巴戟天对复合肥中养分的利用效率表现为N > K > P，其中T3处理N、P、K累积量均高于其他处理。影响巴戟天寡糖累积量的最关键因素是产量，巴戟天产量与全株N/P/K含量、4种寡糖累积量呈显著或极显著正相关，但耐斯糖、1F−果呋喃糖基耐斯糖含量与产量、全株N/P/K含量呈显著或极显著负相关。主成分分析的综合得分由高到低为T3 > T5 > T4 > T2 > T1 > CK，T3处理最佳。
结论建议巴戟天的种植以N肥为主，适量补充K和P肥，适当施肥可显著提高巴戟天药材产量和寡糖累积量。
Abstract:ObjectiveTo study the effects of compound fertilizer application amount on yield, nutrient absorption and accumulation of main effective component oligosaccharides of Morinda officinalis, and provide a theoretical basis for the fertilization management of M. officinalis.
MethodUsing M. officinalis as the research object, through the pot experiment, a total of six treatments were set up to apply 0, 2, 4, 6, 8 and 10 g compound fertilizer (CK, T1, T2, T3, T4 and T5) per plant. The yield, nutrient content and oligosaccharide content of M. officinalis under different fertilization treatments were determined. Then we analyzed the correlation between the indicators, and comprehensively evaluated them by principal component analysis.
ResultThe application of compound fertilizer significantly improved the yield of M. officinalis, increased the N, P and K contents, and contributed to the increase in the accumulation of oligosaccharides. The highest yield and total oligosaccharide accumulation were achieved in the treatment of 6 g compound fertilizer per plant (T3), compared with no fertilizer (CK), with increases of 114.47% and 75.36% respectively. Fertilizer application reduced the oligosaccharide content of M. officinalis to some extent, with a maximum decrease of 35.65% compared to CK. The utilization efficiency of nutrient in compound fertilizer by M. officinalis showed N > K > P, and the accumulations of N, P and K in T3 treatment were higher than those of other treatments. The yield was the most critical factor affecting the accumulation of oligosaccharides in M. officinalis. The yield of M. officinalis was significantly or highly significantly positively correlated with the whole plant N, P, K contents and the accumulation of four oligosaccharides, but the contents of nystose and 1F-fructofuranosyl nystose were significantly or highly significantly negatively correlated with the yield and the whole plant N, P, K contents. The comprehensive score of principal component analysis from high to low was T3>T5>T4>T2>T1>CK, with T3 treatment being the best.
ConclusionIt is suggested that the cultivation of M. officinalis should focus on N fertilizer, and supplement K and P fertilizer appropriately. Proper fertilization can significantly increase the yield of M. officinalis and the accumulation of oligosaccharides.
-
Keywords:
- Morinda officinalis /
- Compound fertilizer /
- Yield /
- Oligosaccharide /
- Nutrient absorption
-
浆果，是一种外表层较为脆弱的肉质果。截至2020年6月，全球浆果的产量约为9932 t，包含小番茄、葡萄、蓝莓、猕猴桃等。随着智能农业装备的发展，为了提升生产效率，越来越多的智能机械被应用于浆果采摘、分拣等生产环节。智能机械末端在夹持果实过程中对浆果造成的损伤程度，将直接影响后续的加工、储存等作业质量，以致影响产品的鲜销价值；因此，智能机械末端的夹持机构是决定浆果智能化生产加工品质的关键部件[1-3]。近年来，国内外学者对浆果的夹持机构进行了大量研究，如柔顺梁[4]、多关节[5]、四叶片软体[6]、三指气动软体[7-9]等夹持机构(图1)。上述四类新型夹持机构结合材料科学，利用柔软材料自身具备的大变形和多自由度特性，增强了末端夹持机构的灵活性和顺应性，这正是末端夹持机构设计的主要发展趋势[10-11]。然而，由于浆果夹持工艺特点，采用单一材料难以同时满足承载能力和变形能力对末端夹持机构刚度的差异化要求。现有夹持机构往往由不同材料的多组件装配而成，不仅动力学建模设计复杂，后期组件的加工、装配等繁杂操作也对机构的动作精度有极大影响[12-14]。因此，如何在保证柔性夹具稳固夹持的前提下，获得易加工、易驱动的轻质柔性夹持装置，实现低损甚至无损完成浆果自动夹持，是有效推进农业机器人在浆果生产领域推广应用所面临的关键问题。
柔性机构是一种少构件、无装配设计的驱动结构，可通过机构自身的弹性形变传递运动至输出端产生驱动响应[15-16]。结合以增材制造技术[17-19]为代表的先进制造技术，可快速成型集驱动连接和柔性夹具为一体的柔性夹持机构。同时，为适应不同种类和不同外形浆果，需调整夹持机构不同区域的刚度来提高其顺应性，以达到稳固、无损地夹持采摘的目的。已有研究通过引入形状记忆聚合物[20]、电活性聚合物智能材料[13]等方法来调整结构不同区域的刚度，但上述2类方法需对结构进一步分段或分层设计，加大了设计和加工难度。而在柔性夹持机构中引入梯度点阵结构[21-23]，可以根据设计需求在设计域内实现刚度的梯度变化。本研究前期提出一种同时优化宏观结构拓扑和介观单胞分布的多层级点阵结构优化方法，实现了具有刚度梯度变化的连续体结构优化设计[23]。鉴于梯度点阵结构在结构变刚度设计中体现出的优异性能，将其引入到以上的柔性夹持机构设计中，通过优化机构整体拓扑构型和内部不同刚度单胞的分布，即可在一体化结构中同时实现高刚度的驱动连接和低刚度的无损夹持。这不仅大大缩短了结构优化设计的时间，并且优化后的多层级结构还具备改善结构整体刚度、可靠性高和局部稳定性好等多方面的优点[24-25]。
综上，本文将针对典型浆果——小番茄，引入多层级拓扑优化理论，提出基于梯度点阵结构的变刚度柔性夹持机构设计方法，分析小番茄夹持过程对夹持机构的性能需求，选取特定的点阵单胞构型并建立其特征尺寸与等效刚度的调控模型，构建变刚度柔性夹持机构的优化设计模型，获得可实现浆果无损夹持的变刚度柔性夹持机构的设计。同时，通过Polyjet增材制造技术加工优化设计获得柔性夹持机构样件，进行茂名‘千禧’小番茄夹持试验，验证设计方法的可行性和柔性夹持机构的有效性。
1. 浆果夹持机构的性能需求
小番茄作为典型的浆果，受损后具备明显的浆果损伤特性，能直观反映夹持机构的夹持效果。以下将通过分析小番茄对夹持机构的性能需求，来体现浆果采摘对末端夹持机构的性能需求。
1.1 夹持机构的可夹持范围
测量随机挑选的100颗茂名‘千禧’小番茄，统计得到该品种小番茄的外形规格为：径向直径平均值为31 mm，样本方差为3.92；横向直径平均值为30 mm，样本方差为2.19，即果实的横径向直径比约为1∶1。夹持机构可从径向和横向2个方向夹持小番茄(图2)：若夹持机构从径向夹持小番茄，夹持口可张开最大距离需大于31 mm；若从横向夹持，夹持口可张开最大距离需大于30 mm。
![]()
图 2 不同加载方式下小番茄的压缩载荷与压缩率的关系CL、CR分别为横向、径向压缩率，RL2、RR2分别为横向、径向决定系数Figure 2. The relationship between compression load and compression ratio of small tomatoes under different loading modesCL and CR are lateral and radial compression ratios respectively, RL2 and RR2 are lateral and radial determination coefficients respectively1.2 夹持机构的材料要求
夹持机构材料首先应具备一定刚度且轻质，既能保证夹持时的可承载度，也可降低夹持机构的自身质量，提高工作的灵活度和准确度。本研究采用新型环保TPU材料作为夹持结构的基体材料，其具有轻质、高强度、高弹性以及良好的耐磨损性。另外，浆果的外果皮薄，且其中果皮和内果皮肉质多汁，不易区分。夹持机构最低损伤地夹持浆果时，应保证果实软接触，即在接触区域设置特殊的缓冲材料，降低浆果的损伤程度。根据常见缓冲材料的弹性模量[26]，指端与果实的缓冲层弹性模量应设置为2.5~11.5 MPa的点阵单胞构型。
1.3 夹持机构的驱动力控制
机械进行夹持作业时，小番茄与夹持指端发生接触部位的皮下组织易变形和积聚液体，对小番茄的品质造成影响，该现象称为小番茄受到机械损伤[27]。在夹持机构驱动力控制下，影响小番茄机械损伤程度的因素主要是夹持载荷、加载位置和压缩率[28]。果实压缩率为压缩变形量与其压缩方向上果实初始直径的比值。通过以下2个试验确定以上驱动力控制下3个因素对小番茄的实际影响。
1) 选择若干个未受损、符合试验条件的小番茄，通过SC-500N数显推拉力计试验机(深测智能有限公司，量程：0~500 N，精度：0.1 N)分别采取横向和径向2种加载方式压缩果实至表皮破裂，测量并处理数据得到2种加载位置下不同压缩载荷的小番茄压缩率线性拟合直线(图2)。由图2可知，从横向加载压缩小番茄对其压缩率的影响比采取径向加载更为敏感；夹持机构的设计将侧重考虑横向加载时驱动力对小番茄夹持品质的影响。
2) 选择若干个未受损、符合试验条件的小番茄，通过深测智能有限公司数显推拉力计试验机，采取横向加载方式，压缩载荷Fp于小番茄的赤道面从0开始压缩小番茄，每间隔0.3 N为1组，每组重复3次“加载—卸载”操作，至其表皮出现裂纹时停止，期间测量并记录每个被压缩果实的初始横向直径和压缩变形量。所有果实完成压缩操作后，均将其置于同一环境(室温26 ℃)的货架进行储藏，定期观察并记录被测小番茄(选择其中1个未被压缩的小番茄为对照)的货架期(t)[29]。被测小番茄的货架期为小番茄完成压缩后置于同一环境货架的储藏时间，若表皮出现裂纹，小番茄货架期为0；若表皮没有裂纹，小番茄货架期从加载试验后开始计算至出现皱痕的那一刻结束，试验小番茄的货架期均记录于表1。最终根据公式(1)(2)计算得到各个小番茄的横向压缩率(CL)和机械损伤度(η)。
表 1 试验结果Table 1. Test result压缩载
荷/N
Compression
load横向压缩率/%
Lateral
compression
ratio货架
期/d
Shelf
life机械损伤度/%
Degree of
mechanical
damage0 0 8.0 0 0.3 0.22 8.0 0 0.6 1.36 8.0 0 0.9 1.91 8.0 0 1.2 2.37 7.5 6.25 1.5 3.39 7.5 6.25 1.8 3.52 7.0 12.50 2.1 4.28 7.0 12.50 2.4 4.55 7.0 12.50 2.7 5.01 6.0 25.00 3.0 6.15 5.5 31.25 3.3 6.30 5.0 37.50 3.6 6.17 3.0 62.50 3.9 7.26 1.0 87.50 4.2 7.82 0 100.00 $$ C_{\rm{L}}=横向压缩变形量/初始横向直径\times 100{\text{%}}， $$ (1) $$ \eta = \dfrac{{t - {t_1}}}{t} \times 100{\text{%}}， $$ (2) 式中，t为对照组果实的货架期；t1为被压缩果实的货架期。
由表1可知，横向加载下，当小番茄受到压缩载荷为0.9 N时，其横向压缩率为1.91%，货架期与未被压缩的样本相同，此时机械损伤度为0；当小番茄受到压缩载荷超过4.2 N后，表皮出现了裂纹，其横向压缩率为7.82%，此时机械损伤程度为100.00%，该值又称小番茄的组织损伤阈值。因此，小番茄末端夹持机构的设计若要实现无损夹持，需控制小番茄末端夹持机构驱动力所引起的小番茄压缩率不超过1.91%。
针对以上对小番茄末端夹持性能需求的分析可知：小番茄柔性夹持机构的设计，将从横向位置加载夹持小番茄，且其夹持口可张开最大距离需大于30 mm；将选用新型环保TPU材料加工机构，设计指端弹性模量为2.5~11.5 MPa的中密度缓冲结构，以制作出满足夹持需求的轻质高品质夹持机构；需控制机构的输入驱动力范围，以保证横向夹持输出载荷对果实的压缩率为0~1.91%，实现无损夹持。
2. 基于梯度点阵结构的浆果变刚度柔性夹持机构的优化设计
2.1 浆果变刚度柔性夹持机构的优化模型
针对典型的浆果夹持工况，浆果变刚度柔性夹持机构优化模型的设计方法如图3所示。柔性机构采用单输入−单输出模型[30]，设计域左端为输入端，输出位移位于右边界，设计域上下部分固定。鉴于该模型具有对称性，取下半部分优化区域进行设计。在设计域内考虑几何非线性大变形情况[31-32]，从输入端施加驱动载荷(Fin)，通过机构的弹性变形传递运动，在输出端产生与驱动方向垂直的虚拟激励位移(Uout)。夹持接触区域设置弹性模量为Eb的中密度缓冲层，以实现指端与果实的软接触。输入、输出弹簧刚度分别设为kin、kout，以模拟输入驱动端和夹持果实端的刚度。点阵单胞构型采用O形构型，其特征尺寸(xi)为孔洞壁厚在水平方向上的宽度，通过数值均匀化方法[33](公式3)建立特征尺寸与单胞等效刚度的调控模型(图4)。
![]()
图 3 浆果变刚度柔性夹持机构设计图Fin：输入端驱动载荷；Uout：输出端位移；kin、kout：输入、输出端弹簧刚度；Eb：中密度缓冲层的弹性模量Figure 3. The design diagram of berry compliant clamping mechanism with variable stiffnessFin: Drive load at input end; Uout: Displacement at output end; kin, kout: Spring stiffness at input and output ends; Eb: Elastic modulus of medium density buffer layer$$ \begin{split} &\displaystyle\int_{{\Omega ^{{\rm{MA}}}}} {\varepsilon \left( u \right)} :D\left( {{x_i}} \right):\varepsilon \left( \nu \right){\rm{d}}{\Omega ^{{\rm{MA}}}} = \displaystyle\int_{{\Gamma ^{{\rm{MA}}}}} F_{\rm{t}} \nu {\rm{d}}{\Gamma ^{{\rm{MA}}}} +\\ &\quad \quad \displaystyle\int_{{\Gamma ^{{\rm{MA}}}}} {F_{\rm{b}} \nu {\rm{d}}{\Omega ^{{\rm{MA}}}}} ,\nu \in {U_{{\rm{ad}}}} ， \end{split} $$ (3) 式中，ε是应变场，u是宏观下产生的应变，υ是动态空间Uad内虚位移场，Ft与Fb分别是作用于机构的面力与体力，ΩMA为宏观设计域，D(xi)是每个单元的弹性张量。
![]()
图 4 随特征尺寸变化的O形单胞等效刚度变化图D11、D12、D33均为以特征尺寸xi为变量的多项式Figure 4. Change chart of equivalent stiffness of O-shape unit cell structure with characteristic sizeD11, D12 and D33 are polynomials with characteristic size xi as variable设机构输出端位移(Uout)与夹持区域宽度(L)的比值为其输出端的变形率。以机构输出端变形率的最大化为目标，其优化列式如下：
$$ \left\{ \begin{array}{l}find:{x}_{i}\left(i=1,2,\cdots ,N\right)\\ \underset{{x}_{i}}{\mathrm{max}}:\dfrac{{U}_{{\rm{out}}}}{L}=\dfrac{{U}^{{\rm{T}}}{{\boldsymbol{l}}}_{{\bf{out}}}}{L}\\ {\rm{s.t.}}:G\left({x}_{i}\right)={ {\displaystyle\int }_{{\Omega }^{{\rm{MA}}}}V\left({x}_{i}\right){\rm{d}}{\Omega }^{{\rm{MA}}}-{V}_{\mathrm{max}}\leqslant 0}\\ KU={F}_{{\rm{in}}}\\ 0 < {x}_{\mathrm{min}} < {x}_{i}\leqslant 1 \end{array} \right.，$$ (4) 式中，xi为介观单元的特征尺寸，max为函数值最大化，xmin为最小介观单元的特征尺寸，N为设计域内单元总数，lout为输出自由度处值为1的单位矢量，G(xi)为设计域内材料使用量，Vmax为材料使用总量上限值，V(xi)为第i个单元的材料使用量，K为总刚度矩阵，U为设计域位移。
2.2 小番茄变刚度柔性夹持机构的优化
通过Matlab代码构建二维的基于点阵结构的小番茄变刚度柔性夹持机构优化设计。优化模型中具体参数设定为：TPU材料的弹性模量E=29 MPa、泊松比=0.4(各向同性的假设材料)，设计域大小为60 mm × 120 mm，单胞数目为15 × 30，夹持区域为20 mm × 40 mm，缓冲层区域宽度为8 mm、弹性模量Eb=10.15 MPa，输入端驱动载荷Fin=1 N，设计域内材料使用总量上限值Vmax=30%，输入弹簧刚度kin=0.1 N/mm、输出弹簧刚度kout=2.0 N/mm。迭代过程采用优化准则法[34]，其目标函数的迭代过程如图5A所示。结合图5B可知，当迭代次数为20~27时，得到的二维拓扑构型机构的输出端压缩率最大，则该拓扑构型为小番茄柔性夹持机构的最优二维拓扑构型。值得注意的是：夹持机构的设计域和优化参数可根据实际场景需求定义，再通过优化算法求得当前参数下的优化构型，但指端接触区域的刚度是根据果实特征提前设置，因此，即使优化得出的夹持机构构型会随设计域和参数发生改变，依然可以保证其无损夹持的特性。
3. 浆果变刚度柔性夹持机构性能试验
3.1 小番茄柔性夹持机构试验过程
采用Polyjet增材制造技术设计小番茄柔性夹持机构。将“2.2”得到的最优二维拓扑构型(图5)对称操作得到完整的小番茄柔性夹持机构的二维拓扑构型(图6A)。整体尺寸为120 mm × 120 mm，单胞大小为4 mm × 4 mm，夹持区域为40 mm × 40 mm，缓冲层区域宽度为8 mm。柔性夹持机构采用柔性TPU材料进行打印，考虑到样品厚度对结构在变形过程中的屈曲影响，选择沿Z轴打印厚度为15 mm，使其与整体结构尺寸相匹配，最终获得质量约45 g的小番茄柔性夹持机构样件(图6B)。
![]()
图 6 小番茄变刚度柔性夹持机构Figure 6. Compliant clamping mechanism with variable stiffness of small tomato如图7所示，将柔性夹持机构样件固定于数显推拉力计试验机上进行试验。本试验采用不锈钢圆柱形压头施加驱动载荷，机构样件上方的中心区域作为驱动载荷输入端(图7A)。当柔性夹持机构夹持小番茄时，一对方向相反的法向力施加于指端，即从小番茄的赤道面的2个接触点(图7B)进行夹取。将型号为DF9-16、量程为20~500 g、精度为0.01 N的柔性压力传感器安装于柔性夹持机构的指端(图7C)，用于测试指端压力。
选择20个未受损、符合试验条件的小番茄进行试验。试验前，测量每个小番茄的横向直径(L0)。试验期间，测量夹持机构不同的驱动力载荷以及对应夹持口两端的距离(L1)，此时小番茄的压缩率为(L0−L1)/L0。以夹持机构夹持住果实，且机构受到轻微晃动果实不掉落的情况，记为成功夹持果实。夹持机构的驱动力载荷继续增大至小番茄的压缩率为1.91%，此时对应指端柔性压力传感器的数值为1.2 N。试验后，计算并记录该夹持过程中小番茄压缩率为0~1.91%时机构对应的输出端驱动力载荷。
3.2 试验结果与分析
所有被测小番茄的压缩率为0~1.91%时，对应夹持机构在输入端施加的驱动载荷为0~14.56 N(图8)。当小番茄的压缩率为0.90%时，输入端施加驱动载荷为11.00 N，夹持机构可成功夹持果实；当小番茄的压缩率为0.90%~1.91%，对应输入端施加驱动载荷范围为11.00~14.56 N，夹持机构可稳固无损夹持果实。因此，根据小番茄变刚度柔性夹持机构设计，为达到作业时稳固且几乎无损夹持小番茄，在设置夹持机构的驱动装置时，驱动力应当大于11.00 N，但不能超过14.56 N。由于小番茄尺寸差异较小，不同尺寸样本下成功夹持与无损夹持对应的压缩率具有重叠区域，可根据该重叠区域对应的驱动力确定夹持机构的驱动力控制区间。若夹持对象的尺寸差异较大，将出现成功夹持与无损夹持对应的压缩率无单一的重叠区域；此时，可根据重叠区域的分布情况选取若干驱动力控制区间及其对应的果实尺寸区间，但需要借助其他辅助装置对果实的尺寸进行预判。
![]()
图 8 柔性机构样件夹持小番茄时输出端压缩率与输入端驱动载荷的关系Figure 8. Relationship between compression ratio at output end and driving load at input end of the compliant clamping mechanism sample while picking up small tomatoes4. 结论
本文以夹持外表层较为脆弱的浆果类果实为研究背景，在柔性机构设计中引入梯度点阵结构，通过构建柔性夹持机构优化模型，获得可无损夹持浆果类果实的变刚度柔性夹持机构设计。以茂名‘千禧’小番茄为例，通过加工获得质量约45 g的小番茄变刚度柔性夹持机构样件，并进行试验得出：该夹持机构样件在输入端驱动载荷为11.00~14.56 N时，小番茄的压缩率为0.90%~1.91%，对应的机械损伤度为0。若将样件应用于智能机械末端以实现浆果的智能化高效无损夹持，在设置夹持机构的驱动装置时，驱动力应当大于11.00 N，但不能超过14.56 N。
本文提出了基于梯度点阵结构的浆果变刚度柔性夹持机构设计方法，可得到轻质、易驱动、顺应性强且几乎无损夹持的浆果夹持机构，在一定程度上促进了浆果自动化采摘、分拣等环节中夹持装备设计的发展。
-
![]()
图 1 复合肥不同施用量对巴戟天不同部位N、P、K含量的影响
各小图中，相同部位柱子上方的不同小写字母表示处理间差异显著(P<0.05，Duncan’s法)
Figure 1. Effects of different application amounts of compound fertilizer on contents of N, P and K in different parts of Morinda officinalis
In each figure, different lowercase letters on the columns of the same part indicate significant differences among different treatments (P<0.05, Duncan’s method)
![]()
图 2 复合肥不同施用量对巴戟天各部位N、P、K单株累积量的影响
各小图中，不同小写字母表示处理间地上部、根部养分累积量差异显著，不同大写字母表示处理间全株养分累积量差异显著(P<0.05，Duncan’s法)
Figure 2. Effects of different application amounts of compound fertilizer on the accumulation of N, P and K per plant in different parts of Morinda officinalis
In each figure, different lowercase letters on the columns indicate significant differences in aboveground and root nutrient accumulation among treatments, different capital letters indicate significant differences in whole plant nutrient accumulation among treatments (P < 0.05, Duncan’s method)
![]()
图 3 复合肥不同施用量对巴戟天4种寡糖含量的影响
相同寡糖柱子上方的不同小写字母表示处理间差异显著(P<0.05，Duncan’s法)
Figure 3. Effects of different application amounts of compound fertilizer on contents of four oligosaccharides in Morinda officinalis
Different lowercase letters on the columns of the same oligosaccharide indicate significant differences among treatments (P<0.05, Duncan’s method)
![]()
图 4 复合肥不同施用量对巴戟天4种寡糖单株累积量的影响
相同寡糖柱子上方的不同小写字母表示处理间差异显著(P<0.05，Duncan’s法)
Figure 4. Effects of different application amounts of compound fertilizer on accumulation of four oligosaccharides in Morinda officinalis
Different lowercase letters on the columns of the same oligosaccharide indicate significant differences among treatments (P<0.05, Duncan’s method)
![]()
图 5 巴戟天不同指标的相关性分析
“*”和“**”分别表示在P<0.05和P<0.01水平显著相关(Pearson法)
Figure 5. Correlation analyses of different indicators for Morinda officinalis
“*” and “**” indicate significant correlations at P < 0.05 and P < 0.01 levels, respectively (Pearson method)
![]()
图 6 复合肥不同施用量的主成分分析
PC1：第一主成分，PC2：第二主成分
Figure 6. Principal component analyses of different application amounts of compound fertilizer
PC1: The first principal component, PC2: The second principal component
表 1 复合肥不同施用量对巴戟天生长指标的影响1)
Table 1 Effects of different application amounts of compound fertilizer on growth indexes of Morinda officinalis
处理
Treatment地径/mm
Ground diameter总生物量/g
Total biomass根部产量/g
Root yieldCK 5.19±0.35a 6.12±0.31d 0.76±0.04c T1 5.25±0.14a 8.00±0.14c 0.89±0.05c T2 6.32±0.54a 9.65±0.61bc 1.11±0.08b T3 6.07±0.31a 13.96±0.78a 1.63±0.04a T4 6.22±0.44a 10.88±0.24b 1.22±0.08b T5 5.72±0.25a 10.40±0.82b 1.26±0.05b 1)同列数据后的不同小写字母表示处理间差异显著(P<0.05，Duncan’s法)
1) Different lowercase letters in the same column indicate significant differences among different treatments (P<0.05, Duncan’s method)
下载: 导出CSV
表 2 复合肥不同施用量的肥料利用效率1)
Table 2 Fertilizer utilization efficiency of different application amounts of compound fertilizer
% 处理 Treatment N P K T1 36.97±1.37a 2.69±0.31ab 9.16±0.53b T2 27.14±4.41b 2.31±0.07b 8.04±2.05b T3 41.35±2.37a 3.23±0.20a 16.62±2.63a T4 20.97±1.29bc 2.11±0.17b 7.27±1.00b T5 18.38±1.20c 1.18±0.10c 6.11±0.56b 1)同列数据后的不同小写字母表示处理间差异显著(P<0.05，Duncan’s法)
1) Different lowercase letters in the same column indicate significant differences among treatments (P<0.05, Duncan’s method)
下载: 导出CSV
表 3 主成分分析的载荷系数
Table 3 Loading coefficients for principal component analysis
评价指标
Evaluation
index第一主成分
The first
principal
component第二主成分
The second
principal
component总生物量 Total biomass 0.952 −0.083 产量 Yield 0.959 −0.154 全株N含量
Whole plant N content0.843 0.353 全株P含量
Whole plant P content0.708 0.475 全株K含量
Whole plant K content0.961 −0.141 蔗糖积累量
Sucrose accumulation0.946 −0.189 1−蔗果三糖积累量
1-Kestose accumulation0.935 −0.232 耐斯糖积累量
Nystose accumulation0.755 −0.372 1F−果呋喃糖基耐斯糖累积量
1F-Fructofuranosyl nystose accumulation0.618 0.655 特征值 Eigenvalue 6.682 1.057 方差贡献率/%
Contribution of variance74.246 11.748 累计贡献率/%
Cumulative contribution rate74.246 85.994
下载: 导出CSV
表 4 复合肥不同施用量的主成分分析综合得分
Table 4 Comprehensive score of principal component analysis of different application amounts of compound fertilizer
处理
Treatment主成分分值
Principal component value综合得分(F)
Comprehensive
score排序
RankingF1 F2 CK −3.92 −0.49 −2.96 6 T1 −2.28 0.27 −1.66 5 T2 −0.01 −0.86 −0.11 4 T3 3.48 −0.52 2.53 1 T4 1.07 1.17 0.93 3 T5 1.65 0.42 1.27 2
下载: 导出CSV
-
[1] 国家药典委员会. 中华人民共和国药典: 1部[M]. 北京: 中国医药科技出版社, 2020: 83. [2] YANG Z, HU J, ZHAO M. Isolation and quantitative determination of inulin-type oligosaccharides in roots of Morinda officinalis[J]. Carbohydrate Polymers, 2011, 83(4): 1997-2004. doi: 10.1016/j.carbpol.2010.11.006
[3] CHEN D, LI N, LIN L, et al. Confocal mirco-Raman spectroscopic analysis of the antioxidant protection mechanism of the oligosaccharides extracted from Morinda officinalis on human sperm DNA[J]. Journal of Ethnopharmacology, 2014, 153(1): 119-124. doi: 10.1016/j.jep.2014.01.021
[4] YANG X, CHEN D, CHEN T, et al. Oligosaccharides from Morinda officinalis slow the progress of aging mice by regulating the key microbiota-metabolite pairs[J]. Evidence-Based Complementary and Alternative Medicine, 2019, 2019: 9306834. doi: 10.1155/2019/9306834.
[5] ZHU J, PENG Q, XU Y, et al. Morinda officinalis oligosaccharides ameliorate depressive-like behaviors in poststroke rats through upregulating GLUT3 to improve synaptic activity[J]. FASEB Journal, 2020, 34(10): 13376-13395. doi: 10.1096/fj.201902546RR
[6] 葛阳, 康传志, 万修福, 等. 生产中氮肥施用及其对中药材产量和质量的影响[J]. 中国中药杂志, 2021, 46(8): 1883-1892. [7] 刘哲, 钟海蓉, 威则日沙, 等. 氮、磷、钾配施对华重楼根茎的产量和有效成分含量的影响[J]. 中草药, 2019, 50(24): 6103-6113. [8] 邹庆军, 汪涛, 郭巧生, 等. N、P、K调控对野菊药用部位性状及活性成分含量的影响[J]. 中草药, 2019, 50(2): 535-540. doi: 10.7501/j.issn.0253-2670.2019.02.036 [9] 苏文华, 张光飞, 王泽明, 等. 氮、磷和钾肥对灯盏花生长和有效成分积累的影响[J]. 中草药, 2009, 40(12): 1963-1966. [10] 苏现明, 王洪庆, 陈若芸, 等. 巴戟天属植物化学成分及药理活性研究进展[J]. 中药材, 2017, 40(4): 986-991. [11] LIANG J, LIANG J, HAO H, et al. The extracts of Morinda officinalis and its hairy roots attenuate dextran sodium sulfate-induced chronic ulcerative colitis in mice by regulating inflammation and lymphocyte apoptosis[J]. Frontiers in Immunology, 2017, 8: 905. doi: 10.3389/fimmu.2017.00905.
[12] 代文豪. 发酵法炮制巴戟天工艺优化及其活性成分变化研究[D]. 广州: 华南农业大学, 2018. [13] 李倩. 巴戟天炮制过程中化学成分变化规律的研究[D]. 广州: 广州中医药大学, 2015. [14] 朱露林, 沈燚, 易晓薇, 等. 不同产地及种植年限巴戟天中6种寡糖类成分的定量分析[J]. 中草药, 2021, 52(14): 4348-4352. [15] 冼丽铧, 吴道铭, 隆曼迪, 等. 土壤养分、种植年限和种植方式对巴戟天寡糖含量的影响[J]. 华南农业大学学报, 2021, 42(3): 75-85. [16] 邵玲, 李美映, 黎学荣, 等. 肇庆巴戟天两种种植模式的植物学性状及药效成分研究[J]. 中药材, 2019, 42(11): 2480-2485. [17] 隆曼迪. 种植环境及栽培管理措施对巴戟天生长和品质的影响[D]. 广州: 华南农业大学, 2019. [18] 吴礼树. 土壤肥料学[M]. 北京: 中国农业出版社, 2011. [19] 邵镪钎, 李丹, 蒋攀, 等. 氮、磷、钾配施对川明参产量和品质的影响[J]. 中草药, 2018, 49(16): 3926-3932. doi: 10.7501/j.issn.0253-2670.2018.16.029 [20] 梁琴, 陈兴福, 李瑶, 等. 化肥与有机肥配施对川芎产量的影响[J]. 中药材, 2015, 38(10): 2015-2020. [21] 陈晓玉, 贺超, 闫滨滨, 等. N、P、K配施对丹参产量和质量的影响[J]. 中草药, 2019, 50(3): 722-730. [22] 范巧佳, 刘灵, 郑顺林, 等. 春季施氮时期和数量对川芎产量、生物碱和阿魏酸的影响[J]. 四川农业大学学报, 2013, 31(2): 136-139. [23] 贾袭伟, 晋小军, 李慧, 等. 磷酸二氢钾施肥量对款冬生长及产量、质量的影响[J]. 中成药, 2022, 44(1): 314-317. [24] DIMKPA C O, FUGICE J, SINGH U, et al. Development of fertilizers for enhanced nitrogen use efficiency: Trends and perspectives[J]. Science of the Total Environment, 2020, 731: 139113. doi: 10.1016/j.scitotenv.2020.139113.
[25] 洪丕征, 刘世荣, 王晖, 等. 正常和环割条件下不同形态氮素添加对红椎幼苗光合特性的影响[J]. 生态学报, 2016, 36(14): 4485-4495. [26] 潘超美, 詹若挺, 丁平, 等. 不同有机肥对巴戟天生长及矿质养分吸收的影响[J]. 中药材, 2002, 25(10): 699-701. doi: 10.3321/j.issn:1001-4454.2002.10.002 [27] 漆小雪, 韦霄, 陈宗游, 等. 黄花蒿干物质的积累及青蒿素与N、P、K量的动态变化研究[J]. 中草药, 2011, 42(12): 2541-2544. [28] 董娟娥, 张康健, 梁宗锁. 植物次生代谢与调控[M]. 杨凌: 西北农林科技大学出版社, 2009. [29] 冯嘉仪, 谢姗宴, 吴道铭, 等. 氮磷钾配施对银杏果实和外种皮产量及品质的影响[J]. 生态学杂志, 2021, 40(6): 1650-1659. [30] 胡佳栋, 毛歌, 张志伟, 等. 不同施肥处理对党参产量和次生代谢物含量的影响研究[J]. 中国中药杂志, 2017, 42(15): 2946-2953. [31] 卢挺, 杨全, 唐晓敏, 等. 氮磷钾配比施肥对广金钱草产量及质量的影响[J]. 广西植物, 2014, 34(3): 426-430. [32] 吴向维. 基于巴戟天寡糖的喘可治注射剂质量评价研究[D]. 广州: 广州中医药大学, 2014. [33] 邓秋林, 杨正明, 陈雨, 等. 氮磷钾配施对瓦布贝母产量及总生物碱质量分数的影响[J]. 西北农业学报, 2019, 28(7): 1138-1146. [34] 鲁泽刚, 朱永全, 卢迎春, 等. 氮、磷、钾施用对灯盏花产量和主要提取物收获量的影响[J]. 植物科学学报, 2019, 37(1): 55-62.
计量
- 文章访问数: 0
- HTML全文浏览量: 0
- PDF下载量: 0
