基于单粒分离-稳定投种的小麦宽幅匀播技术与装置研究

王龙宝 · 工学博士（农业机械化工程 · 中国农业大学） · 评审提交版

一、总体评价

本文面向现行小麦宽幅播种由外槽轮排种器配合开沟器撒播后种沟内种子分布不均、断条漏播重播突出的瓶颈，提出基于单粒分离与稳定投种的宽幅匀播方法，研制了双气室小麦精量排种器与气流附壁式投种装置。作者建立小麦种子多面体离散元模型并以休止角标定接触参数，对排种与投种各过程作理论受力分析，搭建CFD-DEM耦合仿真，并完成排种台架、投种台架与整机田间四级试验，给出各环节最优参数组合（排种最优单粒率80.7%、投种最优粒距合格率76.7%、分布变异系数33.1%、田间一类苗占比86.7%），核心装置有授权发明专利与实用新型专利支撑。论文章节结构完整，理论、仿真、台架、田间链条清晰，工作量饱满，工艺路线可行，图表与中英文摘要较规范，具一定工程参考价值。论文也存在若干不宜忽略的问题：理论播种粒距5cm这一设计基准与式（2-1）回算结果不符并传至下游结构参数，标定响应面最优落约束边界且预测误差低于全部实测点，三章寻优最优解均落设计空间边界、验证为重测已测点，投种达标声明混用了规范无阈值的自定指标，田间核心对照四因素同时变化致归因不可分离且全程无统计检验。上述问题分布于设计基准、标定可靠性、寻优方法、达标声明与核心验证多个环节，多属数值核对、口径界定与试验设计完善，可通过修改与补充说明解决，答辩前须较多修改完善。

二、主要问题与修改建议

问题1

式（2-1）推算的理论播种粒距5cm与同式回算结果不符，并传播为全文设计基准。第2章用式（2-1）推算理论播种粒距，正文写最小理论播种粒距约为5cm，第3章式（3-4）据此定出每列型孔数Z等于36，第4章投种台架又以保持理论播种粒距5cm为前提。但把论文自给的同一组输入（亩播量9kg、千粒重43.7g、行距15cm）代回同一印刷式（2-1），得到的是8.63cm，与5cm相差约73%。由于这个5cm是后续结构参数链的起点，而它与同一公式直接回算的8.63cm自相矛盾，型孔数与投种粒距判据链的起点随之存疑，若取8.63cm则Z约为21。建议作者：1. 核对式（2-1）中系数4的来历，并明确理论播种粒距的定义标准，是每行单粒距还是苗带内当量粒距。2. 统一式（2-1）的计算结果、式（3-4）的Z等于36与第4章5cm基准三者口径，若为当量粒距则写出其与单行粒距的换算关系。3. 复核式（2-1）中行距应取单行行距15cm还是苗带间距，使量纲标准一致。

问题2

式（2-20）休止角标定的响应面最优落约束下界，预测误差低于全部实测点。第2章以休止角相对误差最小为目标，在式（2-20）回归模型上求得最优接触参数，称其预测误差为1.96%、验证仿真2.3%。但该最优组合中摩擦参数落在式（2-21）约束区间下边界、不是模型内部驻点，且式（2-20）对该参数只有一次项，极小被模型结构推到下界；预测的1.96%又低于全部15个BBD实测点的最小实测误差2.4%，属外推到数据之下。由于这组参数是第3、4章全部CFD-DEM仿真的接触参数输入，最优落边界与预测低于实测使第2章响应面优化得到可靠最优参数的论证强度有待加强。需说明，用印刷式（2-20）回代得1.19%、与论文1.96%不一致，或源于求解时使用未截断的全精度系数，若作者给出全精度系数即可解释，故此项不直接证伪标定参数本身。建议作者：1. 说明最优参数落约束下界是否因区间设定偏窄，复核最陡爬坡的高低水平是否覆盖更优区，必要时外扩区间重做寻优。2. 给出式（2-20）的全精度回归系数或公开优化复算过程，使预测1.96%可复现。3. 对预测误差低于全部实测点给出解释，或以模型内部驻点解替代当前边界解。

问题3

第2、3、4章响应面最优解均落设计空间边界，验证试验为重测已测点。三章都用BBD响应面求多目标最优并各做一次验证试验，但三章最优解分别落在设计空间的边界顶点或棱边上而非内部驻点，其中第3章最优编码恰为表3-4第11组、第4章最优编码恰为表4-7第14组，所谓验证试验实为对已测边界点的同点重测。由于响应面优化的意义在于在所测试验点之外找到更优工作点、再用一个不等于任何已测行的新工况独立校验，三章最优落边界且验证非独立，使响应面优化获得最优工作参数这一贯穿三章的方法学结论缺少独立验证支撑；若真实最优确在边界，则提示三章因素水平区间可能均未覆盖最优区。三章回归式回代均自复现，偏差不超过0.01，故问题在寻优落边界与验证非独立，不是算术错误。建议作者：1. 对三章分别说明最优落因素水平边界的原因，复核各因素区间是否需要外扩。2. 若保留响应面优化结论，各补一组不等于任何已测BBD行的独立工况做验证试验。3. 在结论中如实表述最优参数的适用区间，避免把边界点当作内部最优陈述。

问题4

第2章接触参数标定以10次取极值定因素水平、休止角仅给单一均值，输出却报到5位有效数字。第2章测定杨氏模量、泊松比、堆积密度、各类摩擦与碰撞恢复系数时，均以重复10次取极值的上下限直接作为PB试验高低水平（如杨氏模量30.51至90.63MPa，区间跨约3倍），休止角作为唯一寻优目标只给均值32.9度、未给标准差，而标定输出的动摩擦系数却给到5位有效数字（如0.11694）。由于用极值定区间会被极端样本支配、得到的不是均值附近常规波动，而休止角无标准差则相对误差2.3%相对实测分散度是大是小无从判断，加之标定结果是第3、4章仿真的接触参数输入，输入端取极值、输出端报5位数使参数标定可靠这一结论的论证完整性有所欠缺。这是论证交代方式的问题，并不直接否定标定数值本身。建议作者：1. 以均值加减标准差或置信区间报告本征与接触参数，并据此设定PB试验高低水平，替代取极值。2. 补报休止角的标准差或变异系数，给出相对误差2.3%的可比分散度基线。3. 将标定输出的有效数字位数与测量精度匹配，动摩擦系数保留2至3位。

问题5

第4章投种各项指标满足NY/T1143-2006规范这一声明，混用了规范无阈值的自定指标。第4章投种台架对标NY/T1143-2006，表4-5列出该规范理论粒距小于等于10cm档的四项阈值，其中第四项合格粒距变异系数小于等于40%，论文p101自注该项需剔除重播与漏播粒距。但论文用于评价的分布变异系数（含全部粒距、二维，最优值33.1%）与规范第四项不是同一个量，小结却称投种各项指标满足规范要求。由于满足规范成立的前提是所报指标与规范指标同量同口径，合格率76.7%、重播14.2%、漏播9.2%三项确可对标该档且达标，但33.1%是论文自定的分布变异系数、规范对其并无阈值，把它计入满足规范使达标声明的覆盖范围被高估，而该声明又上行到摘要与结论的稳定投种可行。建议作者：1. 把达标声明限定为合格率、重播率、漏播率三项满足NY/T1143-2006理论粒距小于等于10cm档要求。2. 对自定的分布变异系数33.1%明确说明其无规范阈值，仅作种沟内分布均匀性的自评指标，不计入满足规范。3. 若要以分布变异系数佐证均匀性，补充该指标的对照基线或目标值。

问题6

第5章田间匀播与撒播对照四因素同时变化，核心主张归因不可分离。第5章以宽幅匀播与宽幅撒播对比作业质量，结论称匀播种沟内分布均匀性相对较高。但两组之间开沟器类型、排种器类型、有无投种装置、驱动方式四者同时不同（匀播为芯铧开沟器加双气室精量排种器加气流附壁投种装置加电机驱动，撒播为圆盘或翼铲开沟器加外槽轮排种器加镇压轮驱动加软管撒播），没有任何一个只切换排种或投种系统、其余固定的解耦对照单元。由于把均匀性改善归因于本文核心的排种与投种系统，需要把开沟器、驱动方式等其它变量固定、只切换排种或投种系统，四变量同时改变使该改善也可能来自芯铧开沟器或恒速电机驱动，而第5章是全文核心主张的整机验证环节，归因不可分离会削弱该主张在整机层面的实证支撑。建议作者：1. 补充解耦对照，至少固定开沟器与驱动方式、仅切换排种与投种系统，单独考察其对均匀性的贡献。2. 或把第5章结论表述限定为整机方案层面的对比，不表述为排种或投种系统使均匀性提高。3. 说明撒播侧镇压轮驱动与匀播侧电机驱动的差异对对照公平性的影响。

问题7

第5章田间播深合格率在0.9m/s时为73.3%、未达行业标准，小结只呈现有利指标。第5章宽幅匀播播深合格率在0.6、0.75、0.9m/s三个速度下为83.3%、81.7%、73.3%，随速度递减，论文p121自述播深合格率是评价核心指标、p123自认0.9m/s时低于NY/T996-2006阈值80%，且圆盘式撒播在播深上整体优于匀播。但小结只称匀播播幅与播种均匀性变异系数相对较低，未把播深合格率随速度递减、0.9m/s不达标、反劣于撒播这三点写入对比结论，满足标准也未限定速度范围。由于核心质量指标在高速档不达标却未在结论呈现，第5章提高播种质量的总主张与播深这一关键指标的实测水平之间出现张力，读者会以为整体可行。建议作者：1. 在小结如实写出播深合格率随速度递减、0.9m/s未达标、圆盘撒播在播深上更优的事实。2. 把满足标准限定到达标的0.6至0.75m/s速度段。3. 如保留0.9m/s高速档，分析播深合格率下降的原因与改进方向。

问题8

第5章田间试验全程无方差分析、显著性检验与重复离散度。第5章田间播种质量各指标每组重复3次、植株生长每模式取60株，给出一类苗占比86.7%对48.3%对53.3%、地上部干物质0.43对0.36对0.32g以及株高、分蘖、叶龄等组间差异并直接作结论，但全章无方差分析、无显著性检验、无置信区间，也未给3次重复及60株样本的离散度，而第3、4章台架试验是配有方差分析的。由于组间差异是否超出试验误差需要方差分析、显著性检验或置信区间支撑，田间作为核心验证环节却在统计推断上弱于台架，读者无法判断一类苗占比、干物质等差异是真实差异还是落在试验波动之内，这使苗齐苗壮、分布均匀性相对较高等田间结论的证据强度低于台架。建议作者：1. 对田间各指标的组间差异补做方差分析或显著性检验，给出置信区间或误差棒。2. 报告3次重复及60株样本的离散度。3. 使田间结论的统计严谨性与台架一致。

问题9

第3、4章排种与投种台架的最优性能偏低且无部件级对照基线。排种台架最优单粒率80.7%、漏吸5.8%、重吸13.5%，投种台架最优分布变异系数33.1%、粒距合格率76.7%，均为各自最优参数组合下的结果。但两台架都没有与外槽轮撒播或其它精量排种器、投种方式的部件级对照，80.7%只与自身回归预测81.29%比较，33.1%只与自身回归预测32.64%比较。由于精量、稳定、提高均匀性都是相对性判断、成立需要一个可比基线，两台架的最优指标只与自身回归预测比较属自比而非他比，加之整机环节又有问题6所述四因素混淆，排种与投种环节的相对性结论无法单独建立，核心机制的有效性主要靠自比支撑。建议作者：1. 在排种与投种台架补与外槽轮撒播或同类精量排种器、投种方式的对照试验，给出可比基线。2. 或明确说明80.7%、33.1%相对何种基线属提高。3. 在结论中限定精量、稳定的相对范围，避免无基线的绝对表述。

问题10

第4章投种仿真所示均匀性与台架实测离散度量级悬殊，外推未说明。第4章CFD-DEM仿真（单管50粒）给出投种速度变异系数仅3.1至3.9%，同一装置的台架最优分布变异系数却为33.1%，二者量级相差约10倍；第3章另有隔板平板扰流致8、9、10号型孔总压突变、漏吸增大这一来自仿真的机理结论，但台架以整体单粒率、漏吸率、重吸率为指标，未对该几列型孔做分列统计印证。由于仿真投种速度变异系数与台架分布变异系数度量的不是同一个量（前者不含触轮弹跳、机具速度耦合与种沟着床，后者含重漏、二维），10倍差距本身不构成数值矛盾，但CFD-DEM仿真是第3、4章机制论证的主要手段，仿真所示均匀性与特定型孔漏吸都未被台架同口径指标隔离印证，若不补充说明，仿真结论向台架的外推效度尚不充分。此项指向外推说明的补充，并不否定仿真本身的正确性。建议作者：1. 以与台架同口径的指标（如仿真同时输出种沟粒距分布变异系数）重新比对，使仿真与台架在同一量上对接。2. 或补对8、9、10号特定型孔列的分列漏吸统计，使仿真机理结论获得台架印证。3. 若口径无法统一，明确说明仿真仅作机理定性、不直接预测台架的分布离散度。

三、次要问题

以下为校对、勘误、取值交代与完善性问题，集中说明，建议作者全文通校核对。校对与勘误方面，第3章p64变量名标注串行，分别将种层高度设置为2、4、6、8cm并将种层高度设置为中间水平27.5r/min，其中被赋转速值的应为排种转速，式（3-48）寻优约束负3kPa小于等于X1小于等于负5kPa数值上下限方向异常（负3大于负5），应为负5至负3或注明负压绝对值口径，表3-2引用固相参数如3.3.21章节所述疑为3.3.2.1编号脱点，同一量两名两处（休止角与堆积角并存、气流附壁式投种装置与气力式小麦投种装置并用），建议全文各取一名统一。取值依据与试验交代方面，型孔直径2mm落经验式范围2.41至2.48mm的下界、取2而非中值的依据未给，投种通孔11mm非式（4-1）取整值（取整应为9mm、当量直径取谐均与系数1.8无引文），建议补两处取值依据，BBD编码中心与单因素中间水平两处不一致（影响中心点重复即纯误差的代表性），投种台架工作气压固定取最低值2kPa、依据偏笼统而仿真示气压2至8kPa投种速度差约50%，投种最优作业速度0.6m/s即最低水平、合格率随速度单调降、与宽幅高效作业取向的权衡未讨论，投种台架共测120个粒距未说明测区与重复构成，越冬期匀播株高18.9cm为三组最小、解释为苗齐苗壮但株高较小是否为正面指标未给农学依据，建议逐项补充依据或调整表述。文献方面，蒋杉（2025）一器四行气吸方案、丁启朔（2024）与王淞（2024）二维分布变异系数三处近新引文因检索限制未独立核到原文（其存在性合理，不判作者有误），建议作者补全可核出处，并就创新点2相对既有一器四行方案的差异点作一句辨析。

四、评审结论

较大修改后答辩。

附：评审表（凡科博士）