Effect of nitrogen fertigation management on soil nitrogen fractions and N2O emissions from ratoon sugarcane fields
-
摘要:目的
研究滴灌追氮管理对宿根蔗田土壤氮组分和N2O排放的影响，揭示影响土壤N2O通量的土壤因子。
方法以二代宿根蔗Saccharum officinarum为研究对象，在移动防雨棚内进行2个滴灌灌水量[田间持水量的70%~80%(W0.8)和田间持水量的80%~90%(W0.9)]及3种滴灌追氮比例(等氮量250 kg·hm−2，其中，N0为用作追肥的氮肥全部施用到土壤中，N5为50%土施追氮、50%用滴灌系统施用，N7为30%土施追氮、70%滴灌追氮)的田间试验。在甘蔗生长的各个时期测定蔗田土壤N2O通量、pH和氮组分含量，并分析土壤N2O通量与土壤pH和氮组分含量的关系。
结果土壤N2O通量在施用氮肥和灌水后2 d较高，其中，分蘖后期和成熟期W0.9N5处理的土壤N2O通量显著低于其他处理。W0.9条件下，分蘖后期N5处理的土壤N2O累积排放量分别比N0和N7低47.3%和11.8%，伸长后期N5处理的土壤N2O累积排放量比N7低21.5%。相同滴灌追氮比例下，土壤硝态氮含量表现为W0.9>W0.8，随着灌水量的增加，土壤硝态氮含量有所增加。在伸长初期和成熟期，W0.8 N5处理的土壤铵态氮含量比W0.8 N0高56.4%和71.8%、比W0.8 N7高68.5%和160.3%。在分蘖后期，相同滴灌灌水量下，土壤微生物量氮含量表现为N5>N7>N0，W0.8和W0.9两种灌水量下，N5处理的土壤微生物量氮含量分别较N0处理高120.0%和100%。土壤N2O通量与铵态氮含量之间呈正相关关系(r=0.313)，与硝态氮含量之间呈负相关关系(r= −0.391)。
结论W0.9N5处理可以降低土壤N2O排放，且土壤铵态氮和硝态氮含量影响土壤N2O通量，即土壤铵态氮含量越高，土壤N2O通量越高，而土壤硝态氮含量越高，土壤N2O通量却越低。
Abstract:ObjectiveTo study the effects of drip irrigation nitrogen (N) chasing management on soil N fraction and N2O emission in host cane fields, and reveal the soil factors affecting soil N2O fluxes.
MethodTwo water drip irrigation levels of W0.8 (70%−80% of field water holding capacity) and W0.9 (80%−90% of field water holding capacity) and three drip irrigation N chasing ratios (250 kg·hm−2 of equal N) were carried out in a mobile rain-proof shelter with the second generation host cane (Saccharum officinarum) as the research object. N0 was all the N fertilizer applied to the soil (referred to as soil-applied N); N5 was 50% soil-applied N and 50% applied by drip irrigation system (referred to as drip-applied N); N7 was 30% soil-applied N and 70% drip-applied N in the field trials. Soil N2O flux, pH and nitrogen fraction content were measured at each period of sugarcane growth, and the relationship between soil N2O flux and soil pH/nitrogen fraction content was analyzed.
ResultSoil N2O flux was higher after two days of applying N fertilizer and water irrigation, and the soil N2O flux of W0.9N5 was significantly lower than other treatments at the late tillering stage and maturity stage. Under W0.9, the cumulative soil N2O emissions in N5 at the late tillering stage were 47.3% and 11.8% lower than those of N0 and N7, and 21.5% lower than those of N7 at the late elongation stage. Soil nitrate N content showed W0.9 > W0.8 for the same drip irrigation N rate, and increased with the increase of water irrigation. At the early elongation and maturity stage, the soil ammonium N content of N5 under W0.8 was 56.4% and 71.8% higher than those of N0, was 68.5% and 160.3% higher than those of N7. At the late tillering stage, soil microbial biomass nitrogen showed N5 > N7 > N0 under the same water drip irrigation level, and soil microbial biomass nitrogen of N5 was 120.0% and 100% higher than N0 under both water irrigation rates, respectively. There was a positive correlation between soil N 2O flux and ammonium N content (r=0.313) and a negative correlation with nitrate N (r=−0.391).
ConclusionW0.9N5 treatment can reduce soil N2O emissions, and soil ammonium N content and nitrate N content affect soil N2O flux, i.e., higher soil ammonium N content is associated with higher soil N2O flux, while higher soil nitrate N content is associated with lower soil N2O flux.
-
Keywords:
- Ratoon sugarcane /
- Fertigation /
- Sugarcane field soil /
- Inorganic nitrogen /
- N2O emission
-
三七Panax notoginseng是我国名贵的中药材，市场需求量大[1-2]。根据农艺要求，三七育苗播种行株距均为50 mm，播深为10 mm，属小行株距、浅播深精密播种[3]。目前人工点播为主的模式存在用工成本高、作业效率低、作业质量不均匀等问题，机械化播种能有效提高播种质量进而提高效率、降低成本[4]；但三七种子三轴尺寸集中度低、似圆度低、硬度低、表面粗糙度和湿度大、流动性不好，与其他大田大行株距作物种子相比，机械化播种难度大。开沟、排种、覆土各部分系统协调才能实现三七育苗播种机械化。导种管将种子导入种沟，并约束种子流；由于三七种子的特性，导种管对投种过程及种子落点位置影响较大[5-6]。因此结合三七播种的农艺要求和种子特性研究三七种子离开排种器后的投种轨迹，对设计合理的导种管进而提升播种精度和播种效率具有现实意义。
本研究前期采用机械重力式排种原理，设计了窝眼轮式排种结构，经试验，排种器能排出均匀的种子流[7-8]；因导种管设计不科学，种子的纵向分布均匀性与生产要求差距大，本文将这种差距界定为播种精确性，此差距越大播种精确性越低。精确性是衡量种粒纵向分布均匀性和播种机性能的重要指标，受排种、投种和种床构建等多个环节影响[9-10]。种粒与管壁的碰撞和触土时的弹跳所引起的不规则运动易造成种粒的纵向分布不均匀[11-13]，影响播种精确性，合理的导种管能减少种粒的碰撞和弹跳[14-15]。因此，设计出合理的导种管是提高播种精确性的重要手段[16]。
多年来，国内外学者在导种管方面展开了较多的研究。在低速作业时，改变导种管外形能改变种子的下落轨迹，减轻种子的弹跳，保证种子纵向分布的均匀性[17-19]；在高速作业时，可利用二次投种的方法将种子运移到种床，通过带式导种或V型凹槽拨轮能准确地将种子输送到种沟，大幅提高播种的纵向分布的均匀性[20-23]。此类方法设计的结构较为复杂，不适于三七的精密播种。在导种管内分配均匀气流将种子运输的气力集排式分种器，可以实现宽幅均匀播种[24-25]，在导种管内增加波纹和保证导种管长度一致可增加种子行内稳定性[26]。在投种过程中，种子在导种管内的弹跳与碰撞情况复杂，很难通过理论计算做全面分析，高速摄像技术广泛应用于投种过程分析，运用高速摄像技术可以获得种子速度随位移变化规律、不同投种参数下的投种轨迹[27-29]。上述研究分别从部件结构和运动规律等方面展开，但这些研究主要针对油菜、大豆、玉米，小麦等大田大行株距作物，研究结果不适用于小行株距、浅播深精密播种的三七，因此，有必要开展有针对性的研究。
本研究拟应用力学分析确定影响播种精确性的因素，通过台架试验研究精密播种的最佳工作参数组合，运用高速摄像与图像处理技术研究种子离开排种器后的运动轨迹，以拟合曲线为基础，设计出一种新的导种管，并通过土槽试验验证此设计的合理性。
1. 槽式育苗播种方式及播种机结构
1.1 三七槽式育苗播种方式
云南广泛采用的育苗槽规格为1 440 mm宽，据此设计播种机的播幅为1 400 mm，根据50 mm行距的播种农艺要求，播幅内等距排列28个开沟器和28根导种管，机械化播种时，播种机在槽肩上行走，槽内播种，三七育苗槽及播种示意图如图1所示。
![]()
图 1 三七育苗槽及播种示意图1：槽肩，2：三七种子，3：土槽基质，4：育苗槽Figure 1. Structure sketch of Panax notoginseng seedling trough and seeding1: Slotted shoulder, 2: Panax notoginseng seeds, 3: Soil trough matrix, 4: Nursery trough1.2 三七育苗播种机整机结构
三七育苗播种机整机结构如图2a所示，包括种箱、排种器、导种装置、开沟装置、限位装置、控制柜、传动装置、机架和升降螺杆行走装置等。排种器由窝眼滚筒、V型导种槽、圆弧形护种板等组成，导种装置安装在排种器下部，由28根导种管和投种器组成。如图2b所示，开沟装置由28只仿形开沟轮组成，安装在播种机前部，开沟轮厚20 mm，每两只开沟轮之间的中心距离为50 mm，每两根导种管之间的中心距离也为50 mm，其中导种管中心线与开沟轮中心线处于一条直线上。
![]()
图 2 三七育苗播种机结构示意图1：种箱，2：窝眼排种滚筒，3：导种装置，4：开沟器，5：镇压弹簧，6：限位拉伸手柄，7：限位轮，8：被动地轮，9：电动机，10：控制柜，11：传动装置，12：机架，13：升降螺杆，14：主动轮，15：万向轮Figure 2. Structure sketch of Panax notoginseng seedling planter1: Seed box, 2: Hole roller seed metering device, 3: Seed guide device, 4: Ditching device, 5: Crack spring, 6: Limit stretch handle, 7: Limit wheel, 8: Passive ground wheel, 9: Electric motor, 10: Control cabinet, 11: Transmission, 12: Frame, 13: Lifting screw, 14: Driving wheel, 15: Universal wheel1.3 三七育苗播种机工作原理
播种机播种过程中，电机带动主动轮使播种机前进，安装在机架底部的限位轮保证播种机沿直线行走，主动轮轴通过传动系统使窝眼滚筒转动，种子由窝眼滚筒排种器携带，当种子到达投种点时因重力作用离开排种器，经导种管输送到开沟器开出的种沟内。
2. 投种过程种子力学分析
种子的运动轨迹与种子下落过程中的受力有关，有导种管和无导种管时受力有所不同[30]，无导种管时种子的运动只受重力作用，有导种管时还会受导种管管壁的作用。本研究对2种情况下的种子运动进行力学分析，探究影响种子运动轨迹的因素，并通过对比最终确定影响播种精确性的主要因素。
2.1 有导种管时种子的运动分析
理想情况下种子离开排种器后就与导种管管壁接触并沿管壁下滑，然后离开管壁做类斜抛运动，直至落地。为便于研究三七种子的运动规律，以投种点为坐标原点、以播种机前进速度的反方向为x0轴正方向、以竖直向下为y0轴正方向，建立如图3所示的坐标系O0。
![]()
图 3 种子与导种管管壁接触过程的受力及运动分析图建立一个以刚进入导种管种子中心点为坐标原点、以播种机前进速度的反方向为x0轴正方向、以竖直向下为y0轴正方向的坐标系O0；G为种子重力，N；N为导种管直线部分对种子的支持力，N；f为种子与导种管之间的摩擦力，N；β为导种管与水平线的夹角，(°)；H为种子离开排种器时离沟底的高度，m；V0为种子脱离排种器后初始运动速度，m/s；θ为投种角度，(°)；R1为排种器上型孔分布半径，m；n1为排种器转速，r/min；V2为育苗播种机的前进速度，m/s；V1为种子离开导种管时的初速度，m/sFigure 3. Force and motion analysis graph of seeds which contact with seed guide tube wallEstablish a coordinate system O0 with the center of the seed just entering the seed guide tube as the coordinate origin, the reverse direction of the seeder’s forward speed as the positive direction of x0 axis, the vertical downward direction as positive direction of y0 axis; G is seed gravity, N; N is the support force of the straight part of the seed guide tube to the seed, N; f is the friction between seed and seed guide tube, N; β is the angle between the guide tube and the horizontal line, (°); H is the height of the seed from the seed bottom of the ditch when it leaves the seed metering device, m; V0 is the initial velocity of the seed after it is removed from the seed metering device, m/s; θ is the seed release angle, (°); R1 is the upper hole distribution radius on the seed metering device, m; n1 is the rotation speed of the seed metering device, r/min; V1 is the initial velocity of the seed when it leaves the seed guide tube, m/s; V2 is the forward speed of seedling planter, m/s种子刚进入导种管时，种子的速度与种子在排种器投种口处的速度(V0)相同，此时V0=2πR1n1。在种子沿导种管下滑的过程中，种子在重力(G)、支持力(N)、摩擦力(f)的作用下做匀加速运动，此时的加速度(a)为
$$ a = g\sin \beta - {\mu _1}g\cos \beta \text{,} $$ (1) 式中，μ1为种子与导种管间的摩擦系数。
对式(1)进行2次积分得：
$$ V = at + {V_0} \text{,} $$ (2) $$ {S_1} = \frac{{\text{1}}}{{\text{2}}}a{t^{\text{2}}} + {V_0}t\text{,} $$ (3) 式中，V为种子运动的瞬时速度，m/s；S1为种子运动的位置到坐标原点的距离，m；t表示种子运动时间，s。
种子沿管壁下滑时，可分解为沿x0轴和沿y0轴2个方向的运动，设种子沿导种管内壁滑动的时间为t1，在0≤t<t1时，种子的运动轨迹为
$$ \left\{ \begin{array}{l} x(t) = \left[ {\dfrac{1}{2}\left( {g\sin \beta - {\mu _1}g\cos \beta } \right){t^2} + {V_0}t} \right]\cos \beta \\ y(t) = \left[ {\dfrac{1}{2}\left( {g\sin \beta - {\mu _1}g\cos \beta } \right){t^2} + {V_0}t} \right]\sin \beta \end{array} \right.\text{,} $$ (4) 三七种子沿导种管位移的终点位置满足关系式
$$ y{|}_{t{ = t}_1} = x{|}_{t{ = t}_1}\tan \beta \text{,} $$ (5) 此时，种子的行程S1与导种管长度尺寸L相等，有S1=L。
由式(4)和式(5)可知，种子在导种管内的运动与导种管与水平线的夹角、种子离开排种器时的初速度和运动时间有关。
种子离开导种管后以Vt的初速度做类斜抛运动，此时种子运动微分方程为
$$ \left\{ {\begin{array}{*{20}{l}} {m\dfrac{{{\rm{d}}{x^\prime}_0^2}}{{{\rm{d}}{t^2}}} = 0}\\ {m\dfrac{{{\rm{d}}{y^\prime}_0^2}}{{{\rm{d}}{t^2}}} = mg} \end{array}} \right. \text{,} $$ (6) 式中，
$x_0^\prime $ 、$y_0^\prime $ 分别指种子在离开导种管后在x轴、y轴方向移动的距离。设种子落地的时间是t2，在t1≤t<t2时，种子在坐标系中的运动轨迹为$$ \left\{ \begin{array}{l} x(t) = L\cos \beta + {V_t}t\cos \beta \\ y(t) = L\sin \beta {\text{ + }}\dfrac{1}{2}g{t^2} + {V_t}t\sin \beta \end{array} \right.\text{,} $$ (7) 可以求得
$$ {t_{_2}} = \frac{{\sqrt {2g(H - L\sin \beta ) + {{({V_t}\sin \beta )}^2}} - {V_t}\sin \beta }}{a} \text{,} $$ (8) 整理式(7)并把式(8)带入其中，求得种子着地时的位置满足关系式
$$ \begin{array}{c} y{|_{t = {t_2}}} = \\ \dfrac{1}{2}g{\left[ {\dfrac{{x{|_{t = {t_2}}} - L\cos \beta }}{{\cos \beta \sqrt {V_0^2 + 2\left( {g\sin \beta - {\mu _1}g\cos \beta } \right)L} }}} \right]^2} + x{|_{t = {t_2}}}\tan \beta {\rm{,}} \end{array} $$ (9) 根据能量守恒定律，种子在接触土壤时的相对速度(VH，m/s)为
$$ {V_{H}} = \sqrt {V_0^2 + 2gH - 2\mu_1 gL\cos \beta } \text{。} $$ (10) 由式(7)、(9)和(10)可看出，种子的运动受种子与导种管间的摩擦系数、导种管的倾斜角度、导种管长度、投种高度和种子离开排种器时的初速度等因素的影响。其中摩擦系数、导种管的倾斜角度和导种管长度及结构有关，说明不同结构参数的导种管使种子产生不同形状的运动轨迹。在导种管结构参数形式及安装位置确定的前提下，影响种子运动轨迹的因素包括种子的运动时间、投种高度、种子离开排种器时的初速度。
2.2 无导种管时种子的运动分析
对于无导种管的情况，建立一个以投种点为坐标原点、以播种机前进速度的反方向为x1轴正方向、以竖直向下为y1轴正方向的坐标系O1，种子运动的各参数如图4所示。
![]()
图 4 无导种装置时种子运动结构简图建立一个以刚脱离排种器种子中心点为坐标原点、以播种机前进速度的反方向为x1轴正方向、以竖直向下为y1轴正方向的坐标系O1；H'为种子离开排种器时离沟底的高度，m；V'0为种子脱离排种器后初始运动速度，m/s；θ'为投种角度，(°)；R'1为排种器上型孔分布半径，m；n'1为排种器转速，r/min；V'2为育苗播种机的前进速度，m/s；V'1为种子离开导种管时的速度，m/s；V'1x为V'1的水平分速度；V'1y为V'1的竖直分速度Figure 4. Structural sketch of seed motion without seed guiding deviceEstablish a coordinate system O1 with the center of the seed just separated from seed metering device as the coordinate origin, the reverse direction of the seeder’s forward speed as the positive direction of x1 axis; the vertical downward direction as positive direction of y1 axis; H' is the height of the seed from the bottom of the ditch when it leaves the seed metering device, m; V'0 is the initial movement speed of seeds after they are separated from the seed metering device, m/s; θ' is the seed release angle, (°); R'1 is the distribution radius of the upper hole on the seed metering device, m; n'1 is the rotation speed of the seed metering device, r/min; V'1 is the speed of seeds after leaving the seed guide tube, m/s; V'1x is the horizontal speed of V'1; V'1y is the vertical partial velocity of V'1; V'2 is the forward speed of seedling planter, m/s种子在下落过程中的运动可分解为竖直和水平2个运动[19]。
$$ \left\{\begin{array}{l} V_{1 x}^{\prime}=V_{0}^{\prime} \cos \theta^{\prime} \\ V_{1 y}^{\prime}=V_{0}^{\prime} \sin \theta^{\prime}+g t \end{array}\right.\text{,} $$ (11) 式中，
$V_{1x}^{^\prime } $ 为种子相对排种器在x1轴方向上的速度，m/s；$V_{1y}^{^\prime } $ 为种子相对排种器在y1轴方向上的速度，m/s；t为种子运动时间，s。由式(11)可知，种子离开排种器后的运动速度与种子离开排种器时相对排种器的运动速度有关，同时也受投种角度的影响。
在运动过程中种子只受重力的作用，设种子落地的时间是t3，在0≤t<t3,时，种子的运动轨迹为
$$ \left\{\begin{array}{l} x^{\prime}(t)=V_{0}^{\prime} t \cos \theta^{\prime} \\ y^{\prime}(t)=V_{0}^{\prime} t \sin \theta^{\prime}+\dfrac{1}{2} g t^{2} \end{array}\right. \text{,} $$ (12) 在种子下落过程中，种子在水平方向做匀速直线运动，在竖直方向做自由落体运动。
种子第一次触土时的位置满足关系式
$$ y^{\prime}{|_{t = {t_3}}} = x^{\prime}{|_{t = {t_3}}}\tan \theta^{\prime} + \frac{1}{2}g{\left( {\frac{{x^{\prime}{|_{t = {t_3}}}}}{{{V^{\prime}_0}\cos \theta^{\prime} }}} \right)^2}\text{,} $$ (13) 此时种子的相对速度为
$$ V_{{H^\prime }}^\prime = \sqrt {{V^{\prime}}_0^2 + 2g{H^{\prime} }} \text{,} $$ (14) 种子触土弹跳与种子速度大小有关，由式(14)可知，种子离开排种器时的初速度和投种高度影响种子入土时的相对速度进而影响种子的弹跳。通过式(12)和(13)可知，种子的运动轨迹受种子离开排种器时的初速度、投种角度和种子离开排种器后的运动时间的影响。其中种子离开排种器时的初速度与排种器的圆周速度有关，而运动时间与投种高度有关。因此，在无导种管的情况下影响种子运动轨迹的因素包括排种器的圆周速度、投种高度和投种角度。
对比式(7)、(10)、(12)和(14)可知导种管结构影响种子的运动轨迹；排种器的圆周速度、投种高度约束种子着地时的相对速度；种子离开排种器时的初速度、种子运动时间和投种角度影响种子运动；导种管的曲线形式不同，种子的运动轨迹不同。综合以上内容确定在投种过程中影响种子运动轨迹的主要因素为排种器的圆周速度、投种高度和投种角度。
3. 最优工作参数试验
求出三七育苗播种机的最优工作参数是进行导种管曲线设计的前提，后续将在最优工作参数下通过高速摄像分析种子离开排种器后的下落轨迹，进而研究导种管的曲线。为研究种子离开排种器后的运动规律和运动轨迹，搭建了播种精确性试验台，并进行试验研究。
3.1 试验条件
试验在云南农业大学土槽实验室进行。土槽按照三七机械化播种要求制作，宽1 440 mm，土槽内基质与三七播种时基质一致。设计的三七播种精确性试验台如图5所示，投种角度通过不同角度护种板组合实现调节。
![]()
图 5 播种精确性试验台1：种箱，2：排种器，3：投种器，4：摄像机支架，5：机架，6：电控柜，7：传动链条，8：电动机，9：地轮Figure 5. Seeding accuracy test bench1: Seed box, 2: Seed metering device, 3: Seed pitching device, 4: Camera bracket, 5: Frame, 6: Control cabinet, 7: Transmission chain, 8: Electric motor, 9: Ground wheel试验所用种子为经分级处理后的云南文山三七种子，种粒直径范围为5.5~7.0 mm，千粒质量为119.4 g，平均含水率为60%。
3.2 试验因素与指标
对于播种机而言，当排种滚筒直径确定时，由于排种器圆周速度不便测量，且与机器前进速度为固定传动比，故排种器圆周速度可以由机器前进速度来代替。所以选取播种机前进速度、投种高度、投种角度为此次试验的试验因素。
种子纵向(种子沿机具前进方向)分布主要以播种粒距合格指数来评价，文献[31]以理论播种点与实际播种点距离偏差的标准差作为评价种子纵向分布均匀性的新指标。考虑三七播种的实际情况，选取理论播种点与实际播种点距离偏差的标准差为试验指标。由于三七种子的特性，种子与导种管内壁产生较多碰撞引起触土时弹跳，使种子最终位置变化，影响种子纵向分布；因此增加弹跳率作为一个试验指标来探究播种机的播种精确性，各指标的计算公式如下[9, 31]：
$$ C = \sqrt {\frac{1}{n}\sum\limits_{i = 1}^{n} {{{ {{X_i}} }^2}} } \text{,} $$ (15) $$ T = \frac{{{n_t}}}{{N'}} \times 100{\text{%}} \text{,} $$ (16) 式中，C为理论播种点与实际播种点距离偏差的标准差(以下简称标准差)；Xi为第i个种子的理论播种点与实际播种点的距离偏差，mm；n为实际测量种子数；T为弹跳率，%；nt为发生弹跳的种粒数；N'为观察到的排种总数。
弹跳率的测定方法：试验时，用安放在支架上的摄像机对种子下落过程进行拍摄，试验结束后通过对种子下落视频的慢速播放来观察种子是否弹跳。若弹跳记此种子状态为1，否则为0，最后通过计算发生弹跳的种粒数与观察的总的种粒数的比例，得出此时的弹跳率。
试验包括单因素试验和二次回归正交旋转组合试验，为降低试验误差，每组试验重复3次，取平均值。
3.3 单因素试验
根据三七育苗播种机设计要求，播种机的播种速度应在6~10 m/min，设定试验时播种机的前进速度为8 m/min；按照文献[32]中的试验结果，投种高度在20~40 cm时播种效果较好，设定投种高度为20 cm。选取投种角度15°、25°、35°、45°、55°和65°共6个水平，开展投种角度对播种精确性影响试验，试验结果见图6。
![]()
图 6 投种角度与播种精确性关系曲线Figure 6. The curve of the relation between seed release angle and seeding accuracy由图6可知，投种角度在15°~25°和45°~65°时，弹跳率与标准差呈现出相同的变化趋势。分析原因是三七种子在触土时弹跳的增加使重播、漏播增多，使得种子的纵向分布均匀性变差，标准差变低；相反，种子弹跳减少，种子纵向分布的均匀性变好。弹跳率在投种角度为35°时达到了最低；在投种角度为55°时，标准差达到最大，并和弹跳率一起出现拐点。综合可以得出使标准差小、弹跳率低的投种角度的范围为25°~45°，此时整体播种精确性较好。
3.4 二次回归正交旋转组合试验
3.4.1 试验设计
以播种机前进速度、投种高度、投种角度为试验因素，以标准差和弹跳率为试验指标在土槽上开展三因素五水平正交旋转组合试验。根据本研究前期工作及单因素试验确定投种高度为20~40 cm，投种角度为25°~45°，播种机前进速度为6~10 m/min，因素编码如表1所示，试验方案及结果如表2所示，Z1、Z2、Z3为因素编码值。
表 1 播种精确性试验因素水平表Table 1. Seeding accuracy test factor level table水平
level投种高度(Z1)/cm
Seed release height播种机前进速度(Z2)/(m·min−1)
Seeder forward
speed投种角度(Z3)/(°)
Seed release
angle1.682 40 10.0 45 1 36 9.2 41 0 30 8.0 35 −1 24 6.8 29 −1.682 20 6.0 25 Δj 6 1.2 6 表 2 播种精确性试验方案与结果Table 2. Seeding accuracy test scheme and result序
号
No.投种高度 (Z1)
Seed release height播种机前进速度 (Z2)
Seeder forward
speed投种角度 (Z3)
Seed release
angle标准差(P1)/mm
Standard
deviation弹跳率(P2)/%
Bounce
rate1 −1 −1 −1 82.34 79.27 2 1 −1 −1 82.79 81.67 3 −1 1 −1 96.48 80.06 4 1 1 −1 101.52 82.53 5 −1 −1 1 51.08 67.98 6 1 −1 1 56.32 79.84 7 −1 1 1 56.80 75.95 8 1 1 1 79.59 73.25 9 −1.682 0 0 57.24 75.81 10 1.682 0 0 62.87 80.43 11 0 −1.682 0 60.04 78.79 12 0 1.682 0 104.95 83.70 13 0 0 −1.682 60.17 79.86 14 0 0 1.682 63.07 74.74 15 0 0 0 54.43 77.07 16 0 0 0 54.95 81.07 17 0 0 0 51.66 78.80 18 0 0 0 75.10 82.80 19 0 0 0 61.60 79.33 20 0 0 0 61.05 80.93 21 0 0 0 60.23 82.13 22 0 0 0 59.80 83.65 23 0 0 0 75.80 82.53 3.4.2 回归分析
利用Design-Expert 10.0.3对结果进行多元回归拟合和方差分析，得到标准差(P1)和弹跳率(P2)对试验因素实际值的二次多项式回归方程，并对其显著性进行检验。
1)标准差(P1)的显著性分析：由表3可知，Z2、
$ {Z}_{\text{2}}^{\text{2}} $ 对模型有极显著的影响，Z3对模型有显著的影响，其余因素对模型影响不显著。由F值分析可知，影响标准差因素主次顺序为Z2>Z3>Z1。将P>0.1的因素剔除[33]，并再次进行方差分析，得到各因素对标准差的二次回归方程如式(17)所示：表 3 播种精确性试验方差分析1）Table 3. Variance analysis of seeding accuracy test来源
Source标准差
Standard deviation弹跳率
Bounce rate平方和
Sum of square均方
Mean squareF P 平方和
Sum of square均方
Mean squareF P 模型
Model3 914.67/
(3 600.15)434.96/
(1 200.05)3.27/
(11.14)0.026 3**/
(0.000 2**)224.35/
(161.81)24.93/
(53.94)3.89/
(7.03)0.013 5*/
(0.002 3**)Z1 135.32 135.32 1.02 0.331 8 34.80/
(34.80)34.80/
(34.80)5.44/
(4.54)0.036 5*/
(0.046 5*)Z2 1382.15/
(1 382.15)1 382.15/
(1 382.15)10.38/
(12.84)0.006 7**/
(0.002**)9.33 9.33 1.46 0.248 9 Z3 959.35/
(959.35)959.35/
(959.35)7.20/
(8.91)0.018 8*/
(0.007 6**)90.32/
(90.32)90.32/
(90.32)14.11/
(11.77)0.002 4**/
(0.002 8**)Z1Z2 61.27 61.27 0.46 0.509 5 26.25 26.25 4.10 0.063 9 Z1Z3 63.51 63.51 0.48 0.502 0 2.30 2.30 0.36 0.559 2 Z2Z3 1.88 1.88 0.014 0.907 2 0.01 0.01 0.00142 0.970 5 $ {Z}_{\text{1}}^{\text{2}} $ 15.36 15.36 0.12 0.739 5 24.40 24.40 3.81 0.072 8 $ {Z}_{\text{2}}^{\text{2}} $ 1 263.42/
(1 258.64)1 263.42/
(1 258.64)9.49/
(11.69)0.008 8**/
(0.002 9**)0.29 0.29 0.045 0.835 7 ${Z}_{\text{3}}^{\text{2}} $ 37.52 37.52 0.28 0.604 5 37.15/
(36.70)37.15/
(36.70)5.80/
(4.78)0.031 6*/
(0.041 4*)残差
Residual1 731.32/
(2 045.85)133.18/
(107.68)83.23/
(145.76)6.40/
(7.67)失拟项
Lack of fit1 147.58/
(1 462.11)229.52/
(132.92)3.15/
(1.82)0.072 9*/
(0.202 1)46.31/
(108.85)9.26/
(9.90)2.01/
(2.14)0.181 8/
(0.144 1)误差
Pure error583.74/
(583.74)72.97/
(72.97)36.91/
(36.91)4.61/
(4.61)总和
Total5645.99/
(5645.99)307.57/
(307.57)1)Z1：投种高度，Z2：播种机前进速度，Z3：投种角度；“*”表示在P≤0.05水平差异显著，“**”表示在P≤0.01水平差异显著；“/”下“( )”内数字为剔除不显著项后重新进行方差分析的结果
1) Z1: Seed release height, Z2: Seeder forward speed, Z3: Seed release angle;“*” indicates significant difference at P≤0.05 level, “**” indicates significant difference at P≤0.01 level; Data in “( )” under “/” are variance analysis result after rejecting unsignificant items$$ {P_1} = 447.39 - 92.23{Z_2} - 1.41{Z_3} + 6.29Z_2^2\text{。} $$ (17) 对式(17)进行失拟性检验，如表3所示，其中P>0.1，即失拟不显著，说明回归方程拟合较好。
2)弹跳率(P2)的显著性分析：由表3可知，Z3对模型有极显著的影响，Z1和
$ {Z}_{\text{3}}^{\text{2}} $ 对模型有显著的影响，其余因素对模型影响不显著。由F值分析可知，影响弹跳率因素主次顺序为Z3>Z1>Z2。将P>0.05的不显著因素剔除，并再次进行方差分析，得到各因素对弹跳率的二次回归方程如式(18)所示：$$ {P_2} = 34.56 + 0.27{Z_1} + 2.58{Z_3} - 0.043Z_3^2 。$$ (18) 对式(18)进行失拟性检验，如表3所示，其中P>0.1，失拟不显著，说明回归方程拟合较好。
3.5 最佳参数优化
为找出既满足种子纵向分布均匀，又符合播种机播种精确性要求的参数组合，将标准差最小、弹跳率最小作为评价指标，结合因素边界条件建立参数化数学模型：
$$ \left\{ \begin{array}{l}\mathrm{min}{P}_{1}\left({Z}_{1}、{Z}_{2}、{Z}_{\text{3}}\right)\\ \mathrm{min}{P}_{2}\left({Z}_{1}、{Z}_{2}、{Z}_{\text{3}}\right)\\ \text{s}.\text{t}.\left\{ \begin{array}{l}20\;\text{cm}\leqslant {Z}_{\text{1}}\leqslant 40\;\text{cm}\\ 6\;\text{m}/\text{min}\leqslant {Z}_{2}\leqslant 10\;\text{m}/\text{min}\\ \text{25}^\circ \leqslant {Z}_{\text{3}}\leqslant \text{45}^\circ \end{array} \right.\end{array} \right. \text{。} $$ (19) 利用Design-Expert 10.0.3软件对模型进行优化求解，得到最优参数组合为投种高度20 cm、前进速度7.781 m/min、投种角度42°；此时标准差为51.553 mm、弹跳率为72.31%。为方便安装调试，对参数进行圆整，设置投种高度为20 cm、前进速度为7.8 m/min、投种角度为42°；此时标准差为51.66 mm、弹跳率为72.31%。
4. 投种轨迹曲线研究
“3.5”研究结果显示，无导种管情况下播种的弹跳率和标准差不理想。为设计新型导种管，在“3.5”最优参数下，通过高速摄像技术和图像处理技术对种子的运动轨迹进行研究，以获得导种管最优曲线。
4.1 试验装置与材料
试验地点为云南农业大学试验基地，试验材料和装置包括经分级后直径范围为5.5~7.0 mm的文山三七种子、三七播种精确性试验台、千眼狼5F04M型高速摄像仪、智云CR110手持云台、电脑、白色纸板、光源等。试验过程中仪器安装及拍摄如图7所示。
![]()
图 7 图像试验仪器安装图1：高速摄像仪；2：手持云台；建立一个以排种器中心点为坐标原点、以水平正方向为x2轴正方向、以竖直向下为y2轴正方向的世界坐标系O2；$H^{\prime \prime} $为种子离开排种器时离沟底的高度，m；$V_0^{\prime \prime} $为种子脱离排种器后初始运动速度，m/s；$\theta ^{\prime \prime} $为投种角度，(°)；$R_1^{\prime \prime} $为排种器上型孔分布半径，m；$n_1^{\prime \prime} $为排种器转速，r/minFigure 7. Installation picture of image testing instrument1: High speed camera; 2: Handle stabilizer; Establish a world coordinate system O2 with the center of the seed metering device as the coordinate origin, the horizontal positive direction as x2 axis positive direction, the vertical downward direction as positive direction of y2 axis; $ H^{\prime \prime}$ is the height of the seed from the bottom of the ditch when it leaves the seed metering device, m; $V_0^{\prime \prime} $ is the initial movement speed of seeds after they are separated from the seed metering device, m/s; $\theta ^{\prime \prime} $ is the seed release angle, (°); R"1 is the distribution radius of the upper hole on the seed metering device, m; $n_1^{\prime \prime} $ is the rotation speed of the seed metering device, r/min4.2 试验方法
由于包衣的三七种子为暗红色，表面有白色粉末状的药粉，为方便观察运动轨迹和提取位置的坐标，拍摄的背景选择白色纸板。设置高速摄像仪的拍摄速度为240 帧/s，像素分辨率为1 280×720，保证镜头水平，适当对光线进行补偿，使其能清晰地观测到三七种子。投种角度调至42°，调整试验台减速电机，使播种机前进速度稳定在7.80 m/min，启动并设置恰当的录制时间，待试验结束后保存视频文件至电脑端。
4.3 距离标定
在摄像头视野内采集边长为28 mm × 28 mm的正方形网格图像。计算正方形边长的物理长度与其对应的像素长度的比值[34]，记K=B/A，其中A为正方形网格边长对应的像素长度的平均值；B=28 mm，为标定正方形网格的实际边长。
根据试验结果，图像中水平距离换算标准为1像素对应0.210 93 mm；竖直距离换算标准为1像素对应0.211 11 mm。
4.4 特征提取
将种子下落的视频以时间顺序转换为JPG格式的帧图像序列，其尺寸为1 280像素×720像素。通过对图片进行筛选并编号，共得319粒完整的三七种子下落图像。图像处理过程如图8所示。
4.4.1 图像预处理
对34号种子进行分析，图8a和8b为种子下落图像，图像中种子有一定的拖影，通过拖影可看出种子的轨迹。由图8可以看出前一帧图像(图8a)和后一帧图像(图8b)之间有一定的差距，可运用背景差分法进行目标提取，背景差分法是通过前图像与背景图像的差分来检测目标[35]，其原理可表示为
$$ {d_k}\left( {x,y} \right) = \left| {{f_k}\left( {x,y} \right){\text{ }} - {b_k}\left( {x,y} \right){\text{ }}} \right| \text{,} $$ (20) 式中，fk(x, y)为前一帧，bk(x, y)为背景帧，dk(x, y)是前一帧与背景帧的差分结果。
运用差分法，可得到三七种子提取轨迹图。试验发现，在HSV空间下V分量能更易于将下落种子的轨迹与背景分离，效果如图8c所示。由于显示效果不好，为方便观察，对图8c做简单处理，把图片色温由6 500降到1 500。
![]()
图 8 种子下落轨迹序列图像处理过程Figure 8. Sequential image processing of seed falling trajectory sequence4.4.2 图像分割和轮廓提取
预处理后的图像更易于分割，本研究运用Otsu法进行自适应阈值二值化对图像进行分割。Otsu法是基于统计特性来计算目标和背景的最大类间方差，并找出图像的分割阈值的一种方法[36]。由于预处理后图像会存在一些噪声，阈值分割会错误地将噪声分割出来，运用形态学方法可以对二值化后的图像进行修正[37]。形态学的基本原理包括腐蚀、膨胀运算等。腐蚀是一种消除边界点，使边界向内部收缩的过程，可以用来消除小且无意义的物体；而膨胀是将与物体接触的所有背景点合并到该物体中，使边界向外部扩张的过程，可以用来填补物体中的空洞[37]。本研究将二值化后的V分量图先通过连通域面积阈值过滤，然后再进行腐蚀和膨胀运算，得到如图8d所示的缺失的种子轨迹图，最后将2幅图像进行异或运算得到近似完整的种子轨迹图(图8e)。为更方便地提取出种子的运动轨迹，对异或后的图像进行形态学骨架提取，并取反得到轨迹线的大致形状，如图8f所示。
4.4.3 曲线拟合
首先在每粒种子的图像中拾取3~5个点以方便对其轨迹进行二次曲线拟合，如图9a所示。然后根据像素坐标系与世界坐标系之间的不同关系[37]，对坐标系进行转换，得到三七种子在世界坐标系中实际的运动情况。最后按照此方法处理109颗种子图像并对拟合出的曲线进行汇总，得到图9b的整体拟合曲线。由图9b可知，种子的运动轨迹呈抛物线的形式。
由于曲线过于密集，为方便找出规律，对种子拾取的点进行分析，其散点坐标如图9c所示。在寻找曲线过程中，为了保证数据的可靠性，对数据进行数理统计，发现80%种子集中在区域一，另外的20%集中在区域二，如图9c所示。为达到更明显的拟合效果，选择区域一，对其边界进行拟合，并以拟合好的2条边界曲线横坐标中点为标准点，拟合出曲线，拟合效果曲线如图9d所示，拟合得到曲线方程为式(21)，拟合方程的确定系数R2=0.9997。
$$ {y_3} = 0.063\;58x_3^2 + 5.735{x_3} + 217.2 \text{。} $$ (21) 依据前期试验基础设定播种机前进速度为7.8 m/min，排种器转速为13.7 r/min，投种角度为42°；则种子离开排种器时的初速度为0.138 m/s。由式(21)计算出三七种子的理论下落曲线图9e。由于种子的运动轨迹更加接近二次曲线，在利用Matlab软件进行曲线拟合过程中，通过Cure Fitting Tool工具箱，利用多项式逼近(Polynomial approximation)中的二次多项式(Quadratic polynomial)的形式对曲线进行拟合。首先，结合2条曲线进行描点，然后，通过所描绘的点进行曲线拟合，得到最终的曲线方程式，如式(22)，拟合后确定系数为0.999 5，说明拟合度较高，据此曲线形状设计导种管的轮廓曲线。
$$ {y_{4}} = 0.061 \, 94x_4^2 + 5.709{x_{4}} + 217.2 \text{。} $$ (22) 5. 基于EDEM仿真的导种管设计与试验
导种管的研究包括导种管曲线及截面的研究，“4.4.3”已研究出导种管轮廓曲线，现对导种管截面尺寸进行研究。三七种子在离开排种器落入种沟前，由于种子质量不同、在型孔内姿态不同等因素，可能有横向(在水平面内垂直机具前进方向)位移，种子的横向位移分布决定着导种管截面尺寸。对种子横向位移采用EDEM 软件进行研究，根据“3.5”结果，设置投种角度为42°，投种高度为20 cm，播种机前进速度为7.8 m/min；将种子设置为多球面聚合颗粒(图10a)，根据三七种子三轴尺寸[38]，将种子长、宽、高分别设置为6.35、5.67、5.24 mm。排种滚筒轴向不同窝眼的排种质量相互独立，为缩减仿真所需时间，取排种滚筒轴向1/5段进行仿真。通过三维建模软件建立简化后的排种器三维模型并导入EDEM中，仿真模型如图10b所示。根据文献[39]所研究的种子与ABS塑料的接触参数，在EDEM进行相应参数设置，仿真结束后对排种器下方种子流切片，分析其横向偏移量。
将种子无横向偏移下落时的横向坐标作为种子下落理论横向坐标，将仿真获得的种子在投种口下方20 cm处时横向坐标与理论横向坐标做差运算以求得种子横向偏移量；仿真结束后，种子的横向偏移量分布如图11所示。
从图11中可以看出，93.3%的种子横向偏移量为0~13 mm，故将排种器出口截面的横向尺寸取为26 mm。根据图9c可知区域一种子流下方的宽度约为30 mm，故将出口截面尺寸设置为26 mm×30 mm。导种管采用向后(播种机前进方向为前方)倾斜逐渐收缩的矩形端面有利于形成均匀粒距[40]；“3.5”试验得出的最佳投种角度为42°，投种高度为20 cm，结合导种管安装空间的综合考虑，将导种管入口截面尺寸设置为出口截面纵向尺寸的2倍即26 mm × 60 mm。为了避免卡种，将导种管截面的横向边界设计成圆弧；导种管中间部分根据公式(22)和图9f曲线进行过渡，末端对曲线进行延长并优化。
根据曲线方程以及截面尺寸，在Solidworks中建立导种管的三维模型(图12a)，运用3D打印技术对导种管进行加工。为减少种子在导种管内碰撞引起的触土弹跳，导种管材质采用韧性较好的高性能软尼龙材料，加工效果如图12b所示，安装到试验台上之后效果如图12c所示。
![]()
图 12 拟合曲线型导种管R2为导种管圆弧半径，mm；H2为拟合曲线设计高度，mm；H3为过渡放样曲线高度，mmFigure 12. Fitting curve-shaped seed guide tubeR2 is the arc radius of seed guide tube, mm; H2 is the design height of fitted curve, mm; H3 is the height of transition lofting curve, mm导种管圆弧半径R2的大小与护种板圆弧半径的大小一致，保证了导种管的安装角度；拟合曲线设计高度H2为170 mm，过渡放样曲线高度H3设计为15 mm。将导种管安装到样机上进行土槽试验。试验时，采用和“3.3”“3.4”一致的土槽条件，先对种床进行精整，用本研究前期开发的压轮仿形开沟器[41]开沟，建立良好的种床条件，试验用种子为云南文山三七种子，设置投种高度为20 cm、播种机前进速度为7.8 m/min、投种角度为42°，验证试验重复3次，开沟及播种效果如图13所示。3次验证试验结果平均值如表4所示。
表 4 新型导种管播种验证试验结果Table 4. Verification experiment results of seeding with new seed guide tube试验编号
Number of test标准差/mm
Standard deviation弹跳率/%
Bounce rate1 23.10 41.80 2 23.90 45.30 3 33.70 48.50 均值 Mean 26.90 45.20 根据验证试验统计结果，标准差平均值由无导种管时的51.66 mm降至26.90 mm，弹跳率平均值由72.31%降至45.20%，标准差和弹跳率明显减少；证明在有导种管约束的情况下播种，从导种管排出的种子流更加稳定有效，播出的种子纵向分布更加均匀，粒距均匀性增加，播种精确性提高。
6. 结论与讨论
1)通过分别建立有无导种管情况下种子运动的力学模型，找出影响种子离开排种器后运动的共同因素为排种器的圆周速度、投种高度和投种角度。
2)对投种角度进行单因素试验，得到整体播种精确性较好的范围为25°~45°。以播种机前进速度、投种高度、投种角度为试验因素，以标准差和弹跳率为试验指标进行三因素五水平二次正交旋转组合试验，通过Design-Expert 10.0.3软件优化得出无导种管条件下排种器最优投种参数组合为投种高度20 cm、投种角度42°、播种机前进速度7.8 m/min，此条件下标准差为51.66 mm，弹跳率为72.31%。
3)在最优投种参数组合条件下，基于高速摄像和图像处理技术对种子离开排种器后的轨迹曲线进行分析，得出种子的轨迹曲线方程为
$ {y_4} = 0.061\;94x_4^2 + 5.709{x_4} + 217.2 $ ，并通过种子轨迹分布范围得出种子纵向偏移量为30 mm；通过EDEM软件仿真得到种子横向偏移量为26 mm；为了提高播种精确性，最终得到出口和入口截面尺寸分别为26 mm×30 mm和26 mm×60 mm，据此设计试制出导种管。4)根据得到的导种管尺寸数据，运用3D打印技术试制出导种管，并进行验证试验。经验证试验得出，种子粒距纵向分布均匀性增加，标准差为26.90 mm，弹跳率为45.20%，新型导种管满足了三七育苗播种农艺要求，为播种机导种管的田间播种应用提供了参考依据。
从本研究结果看，加装导种管之后的标准差与弹跳率明显降低，虽然弹跳率降低幅度小于标准差，但种子纵向分布均匀性增加,说明导种管的约束作用较好，使得在弹跳率不是很理想的情况下，播种精确性也得到保障。深入分析弹跳率还不够理想的原因主要有3个：第1，三七种子质量不均匀导致三七种子在导种管内的运动轨迹有偏差；第2，在导种管曲线优化设计时，选用80%的种子轨迹作为约束，导致其余20%的种子可能会与导种管壁碰撞；第3，本研究导种管材料为能吸收种子碰撞能量的高性能软尼龙材料，若导种管采用更软质材料则吸收碰撞效果会更佳，后续将对导种管材料展开研究。
-
![]()
图 1 不同滴灌灌水量(W)和滴灌追氮比例(N)下土壤N2O通量的变化
Figure 1. Changes of cumulative N2O emission in soil under different water drip irrigation amounts (W) and N fertigation proportions (N)
![]()
图 2 不同滴灌灌水量(W)和滴灌追氮比例(N)下土壤N2O累积排放量的变化
同一生育期柱子上方的不同小写字母表示处理间差异显著 (P<0.05，Duncan’s法)
Figure 2. Changes of cumulative N2O emission in soil under different water drip irrigation amounts (W) and N fertigation proportions (N)
Different lowercase letters above the columns at the same growth stage show significant difference among treatments (P<0.05, Duncan’s method)
表 1 滴灌灌水量(W)和滴灌追氮比例(N)对土壤N2O累积排放量影响的方差分析
Table 1 ANOVA of the effects of different drip irrigation amounts (W) and N fertigation proportions (N) on cumulative soil N2O emission
参数 Parameter P 分蘖后期 Late tillering stage 伸长初期 Early elongation stage 伸长后期 Late elongation stage 成熟期 Ripening stage W 0.641 0.004 0.800 0.092 N 0.001 0.040 0.000 0.118 W×N 0.000 0.053 0.000 0.000 
下载: 导出CSV
表 2 滴灌灌水量(W)和滴灌追氮比例(N)对土壤硝态氮含量的影响及方差分析
Table 2 Effects and ANOVA of water drip irrigation amounts (W) and N fertigation proportions (N) on soil nitrate nitrogen content
灌水量 Irrigation amount 追氮比例 N fertigation proportion w(硝态氮)1)/(mg· kg−1) Nitrate nitrogen content 分蘖后期 Late tillering stage 伸长初期 Early elongation stage 伸长后期 Late elongation stage 成熟期 Ripening stage W0.8 N0 7.2±1.6d 8.2±4.1c 18.2±0.8c 39.5±3.5b N5 19.0±1.2bc 13.1±3.1c 43.9±7.1ab 46.7±6.9ab N7 15.0±0.4c 37.9±3.3a 34.9±2.1b 43.7±3.4b W0.9 N0 16.9±2.1c 21.8±4.3b 55.0±2. 6a 50.4±15.6ab N5 26.8±1.3a 36.2±3.9a 33.4±4.2b 49.2±2.1ab N7 22.4±0.9b 40.7±5.2a 40.2±5.5b 71.2±6.3a P W <0.01 <0.01 <0.01 <0.05 N <0.01 <0.01 0.872 0.270 W×N 0.657 <0.01 <0.01 0.284 1) 表中数据为平均值±标准误差；同列数据后的不同小写字母表示处理间差异显著 (P<0.05，Duncan’s法) 1)Values in table are means ± standard errors; Different lowercase letters in the same column indicate significant differences among treatments (P<0.05, Duncan’s test)
下载: 导出CSV
表 3 滴灌灌水量(W)和滴灌追氮比例(N)对土壤铵态氮含量的影响及方差分析
Table 3 Effects and ANOVA of water drip irrigation amounts (W) and N fertigation proportions (N) on soil ammonium nitrogen content
灌水量 Irrigation amount 追氮比例 N fertigation proportion w(铵态氮)1)/(mg·kg−1) Ammonium nitrogen content 分蘖后期 Late tillering stage 伸长初期 Early elongation stage 伸长后期 Late elongation stage 成熟期 Ripening stage W0.8 N0 20.6±1.8c 14.0±0.4b 7.0±0.30c 10.3±1.6b N5 18.1±3.7c 21.9±0.6a 13.4±1.5b 17.7±3.5a N7 22.0±0.8c 13.0±0.3c 19.3±3.8a 6.8±1.9b W0.9 N0 26.5±0.6b 13.5±0.3c 6.9±0.7c 6.9±0.6b N5 32.2±3.1a 9.1±0.3c 7.7±0.8c 5.1±0.9b N7 21.2±0.4c 13.0±0.6c 8.4±1.0bc 5.8±0.8b P W <0.01 <0.01 <0.05 <0.01 N 0.089 <0.01 <0.05 <0.05 W×N <0.01 <0.01 <0.05 <0.05 1) 表中数据为平均值±标准误差；同列数据后的不同小写字母表示处理间差异显著 (P<0.05，Duncan’s法) 1)Values in table are means ± standard errors; Different lowercase letters in the same column indicate significant differences among treatments (P<0.05, Duncan’s test)
下载: 导出CSV
表 4 滴灌灌水量(W)和滴灌追氮比例(N)对土壤微生物量氮的影响及方差分析
Table 4 Effects and ANOVA of water drip irrigation amounts (W) and N fertigation proportions (N) on soil microbial biomass nitrogen
灌水量 Irrigation amount 追氮比例 N fertigation proportion w(微生物量氮)1)/(mg·kg−1) Microbial biomass nitrogen content 分蘖后期 Late tillering stage 伸长初期 Early elongation stage 伸长后期 Late elongation stage 成熟期 Ripening stage W0.8 N0 21.9±0.3e 9.0±0.4d 29.2±1.4ab 83.3±1.9a N5 48.2±1.7a 31.4±2.5b 21.8±2.5b 54.3±5.2b N7 32.5±1.2d 55.4±2.4a 37.0±4.9a 30.6±4.0cd W0.9 N0 21.9±1.9e 17.0±0.7c 30.9±3.9ab 46.0±0.2bc N5 43.8±1.4b 14.4±3.3c 26.4±3.6ab 30.4±9.5cd N7 37.9±1.1c 13.3±2.7c 37.1±1.5a 27.6±2.9d P W 0.77 <0.01 0.44 <0.01 N <0.01 <0.01 <0.05 <0.01 W×N <0.05 <0.01 0.788 <0.05 1) 表中数据为平均值±标准误差；同列数据后的不同小写字母表示处理间差异显著 (P<0.05，Duncan’s法) 1)Values in table are means ± standard errors; Different lowercase letters in the same column indicate significant differences among treatments (P<0.05, Duncan’s test)
下载: 导出CSV
表 5 不同滴灌灌水量(W)和滴灌施氮比例(N)下的土壤pH
Table 5 Soil pH under different water drip irrigation amounts (W) and N fertigation proportions (N)
灌水量 Irrigation amount 追氮比例 N fertigation proportion 分蘖后期 Late tillering stage 伸长初期 Early elongation stage 伸长后期 Late elongation stage 成熟期 Ripening stage W0.8 N0 6.7 6.4 6.3 6.7 N5 6.9 6.9 6.7 7.0 N7 6.7 7.1 7.0 7.1 W0.9 N0 6.5 6.5 6.6 6.9 N5 6.8 6.0 6.6 7.2 N7 7.2 7.1 6.7 6.9
下载: 导出CSV
-
[1] PACHAURI K, MEYER A. Climate change 2014: Synthesis report[J]. Environmental Policy Collection, 2014, 27(2): 408.
[2] CARLSON K M, GERBER J S, MUELLER N D, et al. Greenhouse gas emissions intensity of global croplands[J]. Nature Climate Change, 2017, 7(1): 63-68. doi: 10.1038/nclimate3158
[3] 雷豪杰, 李贵春, 柯华东, 等. 滴灌施肥对两种典型作物系统土壤N2O排放的影响及其调控差异[J]. 中国农业科学, 2021, 54(4): 768-779. doi: 10.3864/j.issn.0578-1752.2021.04.009 [4] 李若楠, 武雪萍, 张彦才, 等. 滴灌氮肥用量对设施菜地硝态氮含量及环境质量的影响[J]. 植物营养与肥料学报, 2015, 21(6): 1642-1651. doi: 10.11674/zwyf.2015.0632 [5] LI Z, ZHANG R, WANG X, et al. Effects of plastic film mulching with drip irrigation on N2O and CH4 emissions from cotton fields in arid land[J]. The Journal of Agricultural Science, 2014, 152(4): 534-542. doi: 10.1017/S0021859613000701
[6] 李平, 魏玮, 郎漫. 不同水分对半干旱地区砂壤土温室气体排放的短期影响[J]. 农业环境科学学报, 2021, 40(5): 1124-1132. doi: 10.11654/jaes.2021-1377 [7] 陈津赛, 孙玮皓, 王广帅, 等. 不同施氮量对麦田土壤水稳性团聚体和N2O排放的影响[J]. 应用生态学报, 2021, 32(11): 3961-3968. [8] 王敬, 程谊, 蔡祖聪, 等. 长期施肥对农田土壤氮素关键转化过程的影响[J]. 土壤学报, 2016, 53(2): 292-304. [9] SHAABAN M, WU Y, KHALID M S, et al. Reduction in soil N2O emissions by pH manipulation and enhanced nosZ gene transcription under different water regimes[J]. Environmental Pollution, 2018, 235: 625-631.
[10] BAGGS E M, SMALES C L, BATEMAN E J. Changing pH shifts the microbial source as well as the magnitude of N2O emission from soil[J]. Biology and Fertility of Soils, 2010, 46(8): 793-805. doi: 10.1007/s00374-010-0484-6
[11] BASALIRWA D, SUDO S, WACAL C, et al. Assessment of crop residue and palm shell biochar incorporation on greenhouse gas emissions during the fallow and crop growing seasons of broccoli (Brassica oleracea var. italica)[J]. Soil and Tillage Research, 2020, 196: 104435. doi: 10.1016/j.still.2019.104435.
[12] SHAO Y, CHEN J, WANG L, et al. Effects of fermented organic fertilizer application on soil N2O emission under the vegetable rotation in polyhouse[J]. Environmental Research, 2021, 200: 111491. doi: 10.1016/j.envres.2021.111491.
[13] 赵岩. 土壤氨氮、亚硝酸盐氮、硝酸盐氮的测定试料保存实验[J]. 干旱环境监测, 2014, 28(4): 182-186. doi: 10.3969/j.issn.1007-1504.2014.04.009 [14] 李振高, 骆永明, 滕应. 土壤与环境微生物研究法[M]. 北京: 科学出版社, 2008. [15] 赵苗苗, 张文忠, 裴瑶, 等. 农田温室气体N2O排放研究进展[J]. 作物杂志, 2013(4): 25-31. [16] 梁东丽, 同延安, OVE E, 等. 灌溉和降水对旱地土壤N2O气态损失的影响[J]. 植物营养与肥料学报, 2002, 8(3): 298-302. doi: 10.3321/j.issn:1008-505X.2002.03.008 [17] GIGUERE A T, TAYLOR A E, SUWA Y, et al. Uncoupling of ammonia oxidation from nitrite oxidation: Impact upon nitrous oxide production in non-cropped Oregon soils[J]. Soil Biology and Biochemistry, 2017, 104: 30-38. doi: 10.1016/j.soilbio.2016.10.011
[18] 冯洁, 张克强, 陈思, 等. 土壤N2O吸收和消耗机制及研究进展[J]. 农业环境科学学报, 2014, 33(11): 2084-2089. doi: 10.11654/jaes.2014.11.002 [19] STEVENS R J, LAUGHLIN R J, MALONE J P. Soil pH affects the processes reducing nitrate to nitrous oxide and di-nitrogen[J]. Soil Biology and Biochemistry, 1998, 30(8/9): 1119-1126.
[20] 张忠学, 陈鹏, 聂堂哲, 等. 不同水氮调控模式对稻田土壤氮素分布与有效性的影响[J]. 农业机械学报, 2018, 49(11): 210-219. doi: 10.6041/j.issn.1000-1298.2018.11.025 [21] 武际, 郭熙盛, 王允青, 等. 氮肥基追比例对烟农19小麦氮素吸收利用及产量和品质的影响[J]. 麦类作物学报, 2008, 28(6): 1021-1027. [22] 陈静, 王迎春, 李虎, 等. 滴灌施肥对免耕冬小麦水分利用及产量的影响[J]. 中国农业科学, 2014, 47(10): 1966-1975. doi: 10.3864/j.issn.0578-1752.2014.10.010 [23] 党建友, 裴雪霞, 张定一, 等. 灌溉方法与施氮对土壤水分、硝态氮和小麦生长发育的调控效应[J]. 应用生态学报, 2019, 30(4): 1161-1169. doi: 10.13287/j.1001-9332.201904.034 [24] 尤昆明, 刘士余, 聂小飞, 等. 基于增强回归树的红壤旱坡花生地气态氮损失影响因素分析[J]. 生态与农村环境学报, 2021, 37(5): 651-660. doi: 10.19741/j.issn.1673-4831.2020.0527 [25] LI Y, ZHENG Q, YANG R, et al. Evaluating microbial role in reducing N2O emission by dual isotopocule mapping following substitution of inorganic fertilizer for organic fertilizer[J]. Journal of Cleaner Production, 2021, 326: 129442. doi: 10.1016/j.jclepro.2021.129442.
[26] 赵国胜, 李伏生, 农梦玲. 等氮滴灌对宿根蔗产量及土壤氧化亚氮排放的影响[J]. 广西植物, 2022, 42(3): 413-421. [27] 吕金东, 张丽媛, 虞娜, 等. 水氮耦合对设施土壤N2O和NO排放的影响[J]. 农业环境科学学报, 2021, 40(6): 1366-1376. doi: 10.11654/jaes.2020-1217 [28] 叶丹丹, 谢立勇, 郭李萍, 等. 华北平原典型农田CO2和N2O排放通量及其与土壤养分动态和施肥的关系[J]. 中国土壤与肥料, 2011(3): 15-20. doi: 10.3969/j.issn.1673-6257.2011.03.003 -
期刊类型引用(4)
1. 王慕川，衣淑娟，李衣菲，张紫恒. 指夹式排种器亚高速导种管设计与试验. 农机化研究. 2025(08): 140-148 . 
百度学术
2. 曹钦洲 ，张兆国 ，王法安 ，解开婷 . 三七生产机械化技术与装备研发现状. 中国农机装备. 2023(07): 58-65 . 百度学术
3. 廖宜涛，孙迈，廖庆喜，王传奇，郑娟，武安阳. 油菜精量播种导种环节种子流有序性分析. 农业工程学报. 2023(23): 23-35 . 百度学术
4. 陈玉龙 ，韩杰 ，兰玉彬 ，张猛 ，金亚琛 ，张正 ，王文君 . 精密排种器组合式导种管设计与试验. 农业工程学报. 2022(24): 14-24 . 百度学术
其他类型引用(14)
计量
- 文章访问数: 500
- HTML全文浏览量: 2
- PDF下载量: 324
- 被引次数: 18
