卢富运 · 耕层水流流速

意见 1：盐热联合示踪“二元垂向流速”假设的直接观测验证有待补充（建议级）

1.1 v1 → v2 核心修订说明

v1 评阅书把这一条列为“必修硬伤”,主要依据是“摘要写『实现了』而 §5.5 自承『尚缺乏直接观测证据』”的强度落差。校验员复核 PDF 摘要发现：

摘要原文核实（校验员独立读 PDF 第 3-4 页）：

“针对地表薄层水流流速测量，本研究提出了盐热联合示踪法……基于此，提出了盐热联合示踪表征薄层水流二元垂向流速的试验假设……当水深约为粗糙凸起高度的 1.3~3.2 倍时，盐热示踪流速呈现显著差异，该方法具有一定适用性……在多层秸秆垫面下，两种示踪流速几乎不存在显著差异，该方法适用性受限”

——摘要全段无“实现了盐热联合示踪法测量薄层水流流速及其垂向流速差异”的措辞；实际措辞为“提出了 / 试验假设 / 具有一定适用性 / 适用性受限”。v1 引用的摘要原文不存在——必须撤回。

§6.2(2) 创新点原文核实：

“提出了基于盐热联合示踪的薄层水流二元垂向流速分布表征方法。该方法利用盐、热示踪剂在水体中迁移特征的差异，表征薄层水流中不同垂向层位的流速差异，突破了传统单一示踪方法主要反映整体流速特征的局限。研究表明，在接近层流条件下，盐热示踪流速差异具有反映薄层水流垂向流速分异特征的潜力”

——创新点措辞已用“提出了 + 突破了……局限 + 具有……潜力”,未声称“已完成验证”或“已实现分布表征”。v1 引用的“实现了盐热联合示踪法表征薄层水流二元垂向流速分布”措辞过度归纳了创新点的实际强度。

结论：v1 把这一条做成“方法论根基缺失，与作者自承不符”的指控不成立——作者在摘要、§5.5、§6.1.3、§6.2(2)、§6.3(3) 五处的措辞总体是诚实的（用了“提出”“试验假设”“适用性受限”“潜力”“下一步工作”等温和措辞）,并非强声明。意见 1 的实际诉求应是建设性建议——把“直接观测验证”这件作者已识别为“下一步工作”的事前置一部分到论文修改阶段，让创新点的物理基础在答辩稿中就立起来。

严重度从必修降为建议修改。

1.2 “二元垂向流速”概念物理解读

这个概念有两种可能的物理含义，常被混淆：

含义一（常见误读）：水流速度在垂直方向有两个分量——水平分量 + 垂直分量。例如 PIV 测得的二维流场矢量。

含义二（论文实际含义）：水流在不同深度层级（上层 vs 下层）有两个不同流速值。物理对应关系： - 盐示踪剂（NaCl 饱和水溶液）密度大于水 → 倾向沉到底层 → 反映底层水流流速 - 热示踪剂（沸腾蒸馏水）密度小于水 → 倾向浮到表层 → 反映表层水流流速

论文采用的是含义二——以两种示踪剂的密度差异为根据，推断两者运移速度差就对应于“上下两层水流的真实流速差”。

1.3 直接观测验证的方法学意义

作者已在 §5.5 自承“尚缺乏直接观测证据”——这是建议补充的实质动因。

论文未做任何染色法可视化（如向底层和表层分别滴入不同颜色染料、用相机记录运移过程）、未用 PIV（粒子图像测速）、未用多深度速度传感器（如底部 + 表面双测点的微型 ADV）做对照。所谓“上下分层”完全是从“盐示踪流速 vs 热示踪流速的数值差异”反推出来的。

更深一层：同质心法的计量独立性问题

校验员独立核实 §2.3.1 原文：

“本研究采用盐和热两种示踪剂进行流速测量，并将示踪剂的质心速度认为是薄层水流的流速”

式 2-22 至 2-24 中 M_ij 定义为“温度或电解质浓度值”——同一公式套用到两种信号。

也就是说，无论是盐示踪还是热示踪，都是用电导率 / 温度信号取质心位置，再算质心运移速度。两者只是“标记物”不同，“测速原理”是同一个。这种“两个数据点同算法”的对照不构成独立交叉验证——等同于用方法 A 验证方法 A。

在计量学规范下，新方法的有效性验证至少需要以下之一： - 独立基准对比：用另一种原理完全不同的方法测同一物理量（如 PIV 测表面流速、底部布点测底部流速） - 直接观测：用可视化手段直接观察示踪剂在水流中的垂向分布 - 理论闭合：从物理基本方程推导两种示踪剂的密度差异如何映射到速度差异，并给出可量化的预测模型

论文均未做。但这一缺失作者已在 §5.5 自承，并在 §6.3(3) 列为下一步工作——本评审建议作者把这件已列为下一步工作的事前置一部分到论文修改阶段，而不是将其作为方法论根基缺失的指控。

1.4 文献综述空白：盐热联合示踪法在 §1.2.2 没有现有研究综述

§1.2.2 全节按“染料示踪 / 盐示踪 / 热示踪”三类组织文献综述，分别引 Lei 2005/2010、Shi 2012、Dong 2021（盐示踪）和 Tauro 2017、张艳 2021、Abrantes 2018-2019（热示踪）——但对“盐热联合示踪”这一论文核心创新方向未给任何现有研究综述,仅在末段以“本研究拟将盐示踪与热示踪方法相结合”一句过渡。

读者无从判断： - 这是国内外首创？应明确声明 - 已有少量研究而本文做改进？应给同类工作对比 - 已有较多研究但本文做了系统化？应做综述

关键发现：作者本人 2023 年农业工程学报已发表《盐热联合示踪表征薄层水流剖面流速分布》一文——这是论文方法的真正源头。应在 §1.2.2 末段引述，明确“本研究在作者前期工作基础上系统化拓展为完整测量体系”。这条引用补全后，创新点的文献基础就立起来了——也是建议级修订路径的最低成本入口。

1.5 摘要 / §6.2 vs §5.5 展望的措辞强度对照（订正版）

结论：作者在六个关键位置的措辞强度大体一致——都是“方法可行性试验研究 + 适用条件试探 + 下一步直接观测验证”。没有 v1 评阅书声称的“前后强度悬殊”。本意见不再作为指控，而作为建设性建议：在 §1.2.2 引述作者 2023 年前期工作 + 在 §6.3(3) 把直接观测验证列为下一步首要任务并具体化。

1.6 修改路径

路径 A（推荐，工作量适中）：补做一组直接观测试验

• 选 1~2 组适用工况（如砂粒 2 + 坡度 15° + 6 L/min,水深为粗糙凸起高度的 1.4~1.6 倍——表 5-2 显示该工况盐热流速差异显著）
• 加染色法可视化：向底层和表层分别滴入红色与蓝色染料，用相机记录运移过程
• 或加多深度微型流速计：底部 + 中部 + 表面三测点，验证盐 / 热示踪流速是否分别对应底部 / 表面流速
• 在 §5.3 增加一节“二元垂向流速假设的直接观测验证”

路径 B（最低限度，文献综述补充）：

• §1.2.2 末段补一段：综述盐热联合示踪在国内外的现有研究状态，引述作者本人 2023 年农业工程学报《盐热联合示踪表征薄层水流剖面流速分布》（卢富运 等，2023）作为方法源头,明确“本研究在作者前期工作基础上系统化拓展为完整测量体系”
• §6.3(3) 下一步工作把“直接观测验证”作为首要任务，并把表述具体化：如“在适用工况下用染色法或多深度微型流速计开展验证”
• §5.3 / §5.4 维持“试验假设”级表述，并明确“该假设的直接观测验证为后续工作”

1.7 严重度判定（订正）

• v1 误判为致命级 → v2 订正为建议级
• 影响范围：摘要、§1.2.2、第二章 §2.3~§2.5、第五章 §5.3~§5.5、第六章 §6.2~§6.3
• 性质：建议作者把已识别为“下一步工作”的事前置一部分到论文修改阶段；文献综述应补引作者自己的前期工作

意见 2：SPL 单次预测的不确定性应在摘要中显式呈现，而非仅在小结提及

2.1 表 3-5 完整数据对照

全文最关键的精度对比表（§3.4 末段，p.49）：

2.2 SD 对比的关键数据

SPL SD 与 MDTD / Ratio SD 的倍数关系：

最严重在砂土工况：SPL SD 是 MDTD SD 的 30 倍；砂壤土与粉壤土的倍数关系约 3~4 倍。

2.3 SD 高的工程含义

以砂土工况为例计算“任意一次单点测量的预测区间”:

• 砂土工况平均水流通量：约 0.04 mm·s⁻¹（取流速档 0.000–0.105 mm·s⁻¹ 的中位）
• SPL 法 SD:0.0876 mm·s⁻¹
• 单次测量预测 ±1σ 区间：[0.04 − 0.0876, 0.04 + 0.0876] = [-0.05, +0.13] mm·s⁻¹

也就是说，一次现场测量的结果，在 68% 概率范围内可能跨越接近 100% 的相对不确定性。如果想达到 ±5% 的精度，按 √n 收敛规则：

SD_均值 = SD / √n
要 SD_均值 ≤ 0.002 mm·s⁻¹（即 0.04 mm·s⁻¹ 的 5%）
→ n ≥ (0.0876 / 0.002)² ≈ 1920 次

实际现场使用不可能测 2000 次，所以这个数字说明 SPL 法在砂土工况下做单点测量时几乎无法达到 ±5% 精度——必须改进方法本身或限定应用范围。

对比 MDTD 法：

MDTD MAPE 36.52% 看似差，但 SD 仅 0.0029 mm·s⁻¹
n ≥ (0.0029 / 0.002)² ≈ 3 次

意味着 MDTD 测 3 次取平均就可压到 ±5% 的随机不确定性范围，但其系统偏差(MAPE 36.52%)需另想办法校正。

两种方法的性质不同： - SPL:均值精度好（MAPE 低）,但单次随机不确定大 - MDTD/Ratio:均值精度差（MAPE 高）,但单次随机不确定小

这两个性质对工程使用的指导意义截然不同。

2.4 作者已在哪些位置提及 SD?哪些位置未提及？

校验员复核确认：

结论：作者并非全文隐瞒 SD 信息，而是在小结和结论中已提及——本意见的实际诉求是把 SD 信息从“小结/结论”位置提到“摘要/§6.2 创新点”位置。摘要与 §6.2 是论文影响范围最广的两处位置，只读摘要的读者会得出“SPL 全面优势”的误判。

2.5 “标定组 / 验证组”的内部交叉验证性质

§3.3.1 原文：“将非零水流通量下的试验数据交叉选取均匀区分为标定组和验证组两个数据集”

数据规模反推：n=3 重复 × 4 个非零通量档 = 12 个数据点。均分为 6 个标定 + 6 个验证。

这种“同一批试验数据内部均分”在机器学习术语中称为 k-fold 交叉验证（k 折交叉验证）,用于评估模型对同分布数据的泛化能力。但这不是独立预测验证,因为：

• 标定组与验证组来自同一次土柱填装、同一个探针、同一个安装位置
• 没有跨土柱填装的对照（土柱填装的微小不一致会影响通量）
• 没有跨探针个体的对照（不同物理探针之间的制造差异未量化）
• 没有跨安装位置的对照（探针重新插入到不同位置后的间距漂移未评价，§3.6 已自承）

真正的独立验证应在 新装填的土柱 + 新探针位置 上重新测试。论文将其称为“验证”易让读者误以为是独立测试集。MAPE 8.90%~16.99% 是同试验数据内部交叉验证误差，并非完全独立的预测误差。预期跨装填 / 跨探针的独立验证 MAPE 应高于 17%。

2.6 三土质 MAPE 不可比的物理解释

§3.1.2 饱和导水率(§3.1.2 p.40)： - 砂土：34.64 ± 1.34 cm/h - 砂壤土：3.18 ± 0.16 cm/h - 粉壤土：2.70 ± 0.21 cm/h

砂土的饱和导水率是粉壤土的 12.8 倍。导水率不同直接限制了试验中能达到的最高流速：

§3.1.3 三土质试验流速档位： - 砂土：0、0.027、0.060、0.084、0.105 mm/s（最高 0.105） - 砂壤土：0、0.023、0.050、0.077 mm/s（最高 0.077） - 粉壤土：0、0.008、0.023、0.035、0.044 mm/s（最高 0.044）

粉壤土最高流速 0.044 mm/s 只是砂土最高流速 0.105 mm/s 的 42%。

MAPE 是相对误差（误差 / 真值）,分母是真值流速。三种土质在不同流速区间下做精度对比，公平性存疑：粉壤土 MAPE 8.90% 最低主要因为它的流速档位最小（分母最小）,不是模型对粉壤土特别精确。

绝对误差(RMSE)对比反而更说明问题： - 砂土 SPL RMSE = 0.0093 mm/s - 砂壤土 SPL RMSE = 0.0080 mm/s - 粉壤土 SPL RMSE = 0.0040 mm/s

RMSE 大致与流速档大小成正比，说明 SPL 法的绝对误差大致随通量水平按比例放大——这是均匀比例性，不是“对细土质精度更好”。

论文把“MAPE 无显著性差异”等同于“受土质影响较小”,把两个不同的判断混为一谈。准确说法应是：“SPL 法在各土质各自的合适流速区间内 MAPE 均低于 17%”,而非“无差别精度”。

2.7 “土质影响小”声称的统计基础不足

论文 §3.4 末段：“表明土壤质地对 SPL 方法预测精度的影响较小（表 3-5 SPL MAPE 各土质均标 'a'）”——这里 'a' 是 Tukey HSD 多重比较的同字母(P<0.05),意为三组 MAPE 之间没有显著差异。

但 n=3 时 Tukey HSD 的检验功效(power)很低——在小样本下，“未检出差异”和“无差异”是两个不同的判断。统计功效不足时，即使存在真实差异也可能检测不出。准确表述应是“未在小样本下观察到显著差异”,而不是“土质影响较小”——后者隐含“影响真的小”,前者只说“我们没看出来”。

2.8 修改路径

第一组（摘要 / 小结 / 创新点级别）： - 摘要：“MAPE 为 8.90%~16.99%” → “MAPE 为 8.90%~16.99%,但 SD 高于 MDTD/Ratio 法约一个数量级，单次测量不确定性较大” - §3.7 小结：补“SPL 法的均值精度优势伴随单次测量的较大随机不确定性，工程使用时建议多次平均或仅用于趋势监测” - §6.2 创新点：明示 MAPE 与 SD 两组数据并列

第二组（试验设计 / 验证流程级别）： - §3.3.1 明确：“本研究使用的『验证』是同试验数据的 k-fold 内部交叉验证（标定组 vs 验证组均分自同一批 n=3 重复试验）,非独立测试集；预期跨装填 / 跨探针的独立验证 MAPE 应高于本研究报告值” - 补做至少一组独立验证：用新装填土柱 + 重新安装的探针，测一组对照 MAPE

第三组（土质对比级别）： - §3.4 改写：“SPL 法在各土质各自的合适流速区间内 MAPE 均低于 17%,但三种土质的试验流速档位不同（砂土最高 0.105 mm/s、粉壤土最高 0.044 mm/s）,三组 MAPE 反映的是各土质在其自然通量区间下的方法表现，不构成同流速基准下的精度比较” - 或重新设计实验让三种土覆盖相同的流速档

2.9 严重度判定

• 影响范围：摘要 + §3.3.1 + §3.4 + §3.7 + §6.2 + 第三章全章
• 性质：论文表述完整度问题，不是实验数据问题——作者在小结和结论已提及 SD,只是摘要应同步并列呈现

意见 3：五针法方向 MAPE 小于 15 个百分点的声称应附带流速段限制和 1σ 上限说明

3.1 表 4-3 方向 MAPE 完整数据

3.2 摘要声称与实际数据的落差

摘要原文：“矢量 MDTD 法与矢量 Ratio 法测量水流通量角度的 MAPE 均小于 15%”

实际数据： - 低流速段(<0.012 mm·s⁻¹)三种土质 6 个数据点，5 个超过 16%（仅砂土 Ratio 在 20.98% 算“接近 20%”）,最高达 31.26%（粉壤土 Ratio） - 高流速段(>0.012 mm·s⁻¹)粉壤土 MDTD 14.89%、Ratio 14.97%——压 15% 边界

“均小于 15%”的不成立点： 1. 低流速段绝大多数数据超出 15%——这一段被论文从“均”字背后剔除 2. 高流速段粉壤土两个数据虽然均值 < 15%,但标准差 3.91%/3.77%,按 1σ 范围上限： - MDTD:14.89 + 3.91 = 18.80% - Ratio:14.97 + 3.77 = 18.74% - 即 1σ 区间已突破 15%,“<15%”只能在约 50% 概率上成立

按统计常规，要声称“<15%”成立，应给出 bootstrap 95% 置信区间上限低于 15%——论文未做此种统计核验。

3.3 大小测量精度的实际数据（表 4-2,p.58）

6 个数据点中 5 个在 20%~62% 之间，仅砂土高流速段 Ratio(8.56%)是例外。

线性拟合斜率（图 4-3）系统性低估：

所有土质都系统性低估，粉壤土低估幅度最大。

3.4 论文对大小测量精度差的处理

第四章标题：《五探针法测量土壤水流通量方向试验研究》——单点突出方向，标题本身已暗示“大小不是本章重点”

§4.3.2 误差分析：归因为：(1) 五针 4 个间距参数累加（vs 三针 2 个）;(2) 细质地土壤孔隙水流路径复杂性；(3) 局部水流通量方向与整体方向差异

但没有定量分解这三个原因各贡献多少误差，也没给可量化的改进方案。

§6.1.2 主要结论（论文最显著位置之一，已明示分工）：

校验员复核确认，作者在 §6.1.2 主要结论中已经明确写出：

“本研究不将矢量 MDTD 法和矢量 Ratio 法作为土壤水流通量大小的主要测量方法，而主要利用其方向识别优势”

这是论文最显著的位置之一，作者已经在那里明示“SPL 测大小、五针测方向”的分工。v1 评阅书声称“大小测量精度差未被独立警示，仅在 §4.6 末以『协同测量』绕回”需要软化——分工已在主要结论明示。

问题仍在于摘要的措辞结构上把“大小低估”与“方向<15%”并列叙述，容易让读者形成“方法整体精度优良”的整体印象；可以考虑把这一分工从 §6.1.2 提前到摘要末段明示。

3.5 摘要叙述结构问题

摘要把以下两个事实并列叙述：

“通量大小测量值存在低估现象，但通量方向的 MAPE 均小于 15%”

并列结构 + 转折“但”暗示——大小虽然差但方向好，整体来说还不错。读者会得出“五针法整体精度优良”的判断。

实际情形： - 大小：6 个数据点 5 个在 20%~62%——重大短板 - 方向：高流速段 < 15%、低流速段 16%~31%——分段有效

应明示“五针法不能独立用于测量大小，工程上需 SPL 法补充”（这点 §6.1.2 已明示）,并把这一分工从结论位置提前到摘要末段。

3.6 安装角度固定 45° 的普适性问题

§4.1 试验设计：所有水流方向真值均固定为 45°（等于探针安装角）。

矢量分量公式 (2-20)(2-21) 含 sin/cos 项：

|J| = √(J_x² + J_y²)
φ = arctan(J_y / J_x)

• 当真值角度趋近 0° 时：J_y → 0,arctan 在 0 附近敏感性放大，小幅 J_y 抖动会引起 φ 大幅变化
• 当真值角度趋近 90° 时：J_x → 0,同理
• 真值 45° 是最对称工况：J_x = J_y,两个分量的低估幅度若一致则相互抵消（论文 §4.4.2 已声称“两个分量同时低估不影响合成方向”）

问题：作者用单一 45° 的对称工况掩盖了角度依赖性。工程实际中水流方向不可能恰好都是 45°——田间坡面水流的方向受地形、植被分布等影响，可能在任意角度。方法在 0°、30°、60°、90° 等其它角度下精度未知，普适性主张缺乏验证基础。

§4.4.2 关键假设原文：“两个分量方向同时低估不影响其矢量合成方向的测量”

论文未给两个分量自身的精度数据——只给了合成大小、合成方向。把两个低估值的合成方向作为可靠输出，相当于自我抵消的隐性假设——该假设在 45° 工况下勉强成立（对称抵消）,但在其它角度下是否仍成立未做验证。如果两个分量的低估幅度不一致（如 x 分量低估 30%、y 分量低估 50%）,合成方向就会偏移。

3.7 修改路径

第一组（摘要 / 标题级别）： - 摘要：“通量方向 MAPE 均小于 15%” → “在水流通量大于 0.012 mm·s⁻¹ 的工况下，方向 MAPE 平均约 15%（粉壤土约 15% 临界）;低流速段下方向 MAPE 上升至 16%~31%” - 摘要末段：把“SPL 测大小、五针测方向”的分工从 §6.1.2 提前到摘要末段

第二组（数据完整性级别）： - §4.3 增加分量精度数据：列出两个分量自身的 MAPE 与偏置 - §4.4.1 给出 bootstrap 95% 置信区间或 Wilcoxon 配对符号秩检验 p 值，特别是粉壤土高流速段是否在统计上显著低于 15%

第三组（普适性验证级别）： - §4.4.1 或 §4.5 展望补一组多角度试验：至少加做 0°、30°、60°、90° 中的两组，验证方向测量精度在非 45° 真值下是否仍成立 - §4.4.2 把“两个分量同时低估不影响合成方向”改为“在 45° 安装角下、两个分量低估幅度相近时合成方向影响有限，其它角度下精度需进一步验证”

3.8 严重度判定

• 影响范围：摘要 + 第四章全章 + §6.2 创新点（五针法部分）
• 性质：摘要叙述结构问题 + 单一安装角的普适性问题

意见 4：SPL 模型理论基础与参数稳健性的双重论证有待补强

4.1 SPL 模型推导链(§2.1.1)

完整推导链(§2.1.1 p.17-18)：

[毕渥数判据]                    Bi = h·L_c / k < 0.1                  式 (2-1)
       ↓
[集总热容微分方程]              C · dT/dt = -h·A·(T - T_∞)            式 (2-2)
       ↓
[指数衰减解]                    T(t) = T_∞ + (T₀ - T_∞)·exp(-h·A·t/C)  式 (2-3)
       ↓
[Granier 1985 树干液流经验式]   h = h₀ + α₁ · u^β                       式 (2-4)
       ↓
[本研究简化的对数假设]          A/C · h = a + ln(b·J + 1)              式 (2-5)
       ↓
[最终 SPL 模型]                T(t)/T₀ = exp(-(a + ln(b·J + 1))·t)     式 (2-6)

4.2 跨物理过程类比的物理基础质疑

Granier (1985) 树干液流模型的应用场景： - 介质：植物木质部（含管胞、导管的纤维状多孔介质） - 驱动力：植物蒸腾产生的负压差 - 流动特征：非饱和、间歇性、与气孔开闭高度耦合 - 测量目标：液流速度 u（L·dm⁻²·h⁻¹,单位换算后约 10⁻⁷ ~ 10⁻⁵ m·s⁻¹）

本研究饱和土壤水流： - 介质：饱和土壤孔隙（颗粒堆积的连续介质） - 驱动力：水头差产生的压力梯度 - 流动特征：饱和、稳态、与土壤孔隙结构高度相关 - 测量目标：土壤水流通量 J（m·s⁻¹,约 10⁻⁵ ~ 10⁻⁴ m·s⁻¹）

两者的物理过程虽然都涉及“多孔介质中的水流”,但物理机制差距很大——蒸腾驱动 vs 压力驱动、非饱和 vs 饱和、植物组织 vs 矿物颗粒。论文仅以“β 取值范围 (0,1) 时本式趋势类似”作为类比依据（§2.1.1 引 Bush 2010, Steppe 2010 等）,未做物理机理推导。

§3.6 自承：“参数 b 的物理意义有待进一步研究”——支持“模型本质是经验拟合而非物理派生”的判断。

4.3 表 3-4 参数敏感性数据

§3.3.2 OAT(One-at-a-time)敏感性分析(p.46)：

关键发现：参数 a 微扰 ±10% → 通量误差几百倍；参数 b 微扰 ±10% → 通量误差仅 9%~11%。

a 是模型对水流通量的核心参数——它的精度决定整体精度。

4.4 表 3-3 参数 a 实际标定 SD（订正版）

§3.3.1 表 3-3(p.45)——经校验员订正：

a 的标定 SD 相对值为 0.5%~2.2%（v1 评阅书误读为 0.6%~2.2%、且砂土 a 误读为 0.21500;v2 订正）。

校验追溯：PDF 第 45 页（对应 PDF 第 59 张）第 8 行原文：“相对于砂土中为 0.23558 s⁻¹,砂壤土和粉壤土中的参数 a 值分别降低 0.02238 和 0.02282 s⁻¹”——这确认砂土 a 值 = 0.23558,砂壤土 a 值 = 0.21320（0.23558-0.02238=0.21320,与表 3-3 列示的 0.21032 略有出入；原文页码或为“降低 0.02526”的笔误，但表 3-3 数据本身确为 0.21032）;粉壤土 a 值 = 0.21276（0.23558-0.02282=0.21276,与表 3-3 列示一致）。

4.5 表 3-3 vs 表 3-4 的数学矛盾

用敏感性表预测：

若 a 在标定中实际漂移 2.15%（粉壤土工况）,按敏感性表线性外插： - a 微扰 10% → 通量误差 +339% / -283%（粉壤土） - 推算 a 微扰 2.15% → 通量误差约 +73% / -61%

实测 MAPE（表 3-5）：粉壤土 SPL MAPE 8.90% ± 2.16%

两组数据相差 7 倍多——预期误差 +73% / -61% vs 实测 8.90%。

可能的解释： - 敏感性表的“形式微扰”和“实测漂移”不在同一测量层面：表 3-4 的微扰是任意 a 漂移的形式 OAT 分析（脱离实际试验条件）,表 3-3 的 SD 是“同次试验下三个重复”的内部一致性 - 标定与验证用了同一批试验（同一土柱、同一探针）——a 在标定与验证间几乎不变，故验证 MAPE 远低于敏感性预测 - 但若工程使用时换土柱、换探针、跨日重新标定，a 的漂移可能比 2.2% 大得多，那时 MAPE 会跃升到敏感性预测值

论文没有做这层区分——把“OAT 敏感性”和“实际漂移”两个层面的数据并列呈现，让读者无法判断 SPL 法的“实际工程稳健性”。如果一个方法的核心参数微扰 10% 就误差几百倍，但论文只报告“同一次试验下 MAPE 17%”——工程使用者会被这个数字误导。

4.6 集总热容假设事后验证的循环逻辑

§3.3.3 Bi 验证流程(p.46)：

1. 用表 3-3 标定的 a, b 值（基于式 2-5）
2. 代入式 (2-5) 反算对流换热系数 h（在最大水流通量工况下）
3. 用式 (2-1) 算 Bi = h·L_c / k
4. 结果：Bi = 0.007~0.008（三种土质）,< 0.1

循环逻辑问题（铁律 12“先一句白话再放公式”）:

作者用同一公式来证自己的前提。步骤 2 用了式 (2-5),但式 (2-5) 本身是基于“集总热容假设成立”才推导的。也就是说——用了“集总热容假设成立 → 推出式 (2-5)”来反推“集总热容假设是否成立”。这是用 A 推 B 然后用 B 反过来证 A 的循环验证。

更稳健的验证方式：直接测量探针表面与内部的温度梯度（用两个埋设位置不同的热电偶）。如果温度梯度可忽略(< 5%),集总热容假设直接成立。论文未做此种直接测量。

另一个问题：只在最大水流通量工况下验证一个点。中等通量段（如砂土 0.05 mm/s）的 Bi 数是否仍 < 0.1 未做。论文逻辑：“h 由式 (2-5) 知与 J 正相关 → J 最大时 Bi 最大 → 只需验证最大 J 工况”——但式 (2-5) 本身在质疑中，用待证假设证明验证策略合理性，是双重循环。

4.7 §3.5 误差讨论的自承

§3.5 原文(p.49)：“对流换热系数应由探针和饱和土壤界面之间的热通量和温差决定(Duda, 2023)……探针的尺寸限制，其内部温度梯度无法完全被忽略”

——作者自己承认“探针内部温度梯度无法完全被忽略”——这就是集总热容假设的违背。但 §3.3.3 又说“Bi < 0.1,假设成立”。两处自相矛盾。

如果“温度梯度无法完全被忽略”,那么 Bi 算出来小，可能不是因为温度梯度真的小，而是因为对流换热系数 h 在最大通量工况下被低估了——这是模型本身的局限。

4.8 修改路径

第一组（理论基础级别）： - §2.1.1 增补从饱和土壤强制对流边界条件出发的物理推导，给出 h 与 J 关系的理论依据 - 或明确声明式 (2-5) 是“基于经验类比的对数模型，物理基础有限” - §3.6 改写：“a 与 b 均为经验拟合常数，本研究暂未给出严格物理解释。模型在饱和土壤中的应用基于树干液流模型的经验类比”

第二组（数值稳健性级别）： - §3.3.2 改写：把表 3-4 改为“OAT 形式敏感性”和“跨日 / 跨探针实测漂移”两组数据对比，并解释两者差距来源 - 补做实测漂移试验：在同一土柱上跨 3 天、3 个不同探针个体、3 个不同安装位置分别标定 a,量化 a 的实际漂移范围 - 在 §3.6 工程使用边界条件段明确：“每次新安装探针后均需重新标定 a”

第三组（集总热容验证级别）： - §3.3.3 直接测量探针表面 vs 内部温度梯度（用两个埋设位置不同的热电偶） - Bi 验证扩展到 min / median / max 三档分别计算 - 文字改为：“在最大通量工况下 Bi < 0.01,间接说明集总热容假设近似成立”——避免强声明

4.9 严重度判定

• 影响范围：第二章 §2.1.1 推导链 + 第三章 §3.3.1/§3.3.2/§3.3.3/§3.5/§3.6
• 性质：方法学论证不充分，不是实验数据问题

意见 5：盐热示踪 §5.3.3 节命题表述前后不一致，适用条件较窄需在摘要明示，第二章“统一原理”应在小结明示

5.1 §5.3.3 节内表述不一致

§5.3.3 前段（与 §6.1.3 主要结论一致）：

“水流接近层流状态是盐热联合示踪法表征薄层水流二元垂向流速分布的前提”

§5.3.3 后段：

“本研究的试验流量分别为 6、12 和 18 L·min⁻¹,则对应的 Re 分别为 625、1250 和 1875。这表明所有试验工况下的水流都处于层流和紊流的混合状态”

§2.3.2 论文自身的层流判据（引朱昱等 2022）：“Re < 500 层流；> 2000 紊流；500~2000 过渡”

表述对不上的点： - 论文自定义层流上限 Re = 500 - 试验最低 Re = 625（已超过 500） - 试验最高 Re = 1875（接近 2000） - 全部试验都在论文自定义的“过渡区”(500~2000)

但 §5.3.3 同节内一会儿说“前提是接近层流”、一会儿说“全部在层流-紊流混合”——读者无法判断： - 方法适用条件到底是“接近层流”还是“层流-紊流混合”? - 如果前提是“接近层流”,本研究的方法验证是否本身就不在适用条件下？ - 如果方法在“层流-紊流混合”下也可用，“接近层流”的前提声明为何还保留？

“接近层流”也没有定量边界——Re < 1000 算？Re < 1500 算？现场使用者拿到这套方法时无法判断当前条件是否落在适用区。

5.2 盐热示踪法适用条件汇总

§5.5 / §5.6 / §6.1.3 表述汇总：

论文方法的实际适用工况：仅定床砂粒 + 单层秸秆的两个特定水深倍数窗口。

作者已在 §5.5、§5.6、§6.1.3 给出定量边界——并不刻意隐藏适用范围。但摘要把“方法体系建立”与“适用条件汇总”分两段写，读者要读完整段才能形成“这方法在哪些工况下能用”的整体图像；摘要可压缩为一句“该方法在定床砂粒（粗糙凸起 1.3~3.2 倍水深）和单层秸秆（厚度 0.5 倍或 1~2 倍水深）的特定工况下适用性较好，多层秸秆与光滑 PVC 多数工况不适用”,信息密度更高。

5.3 田间推广路径分析

§1.4 技术路线图（图 1-14）：终点为“得出结论和撰写论文”,下游无“田间推广验证”节点。

§6.3(3) 下一步工作建议原文：

“进一步开展田间原位试验，在真实垫面条件下结合秸秆移动、土壤侵蚀及降雨等环境因素，对该方法表征水流二元垂向流速的适用条件进行完善与补充”

田间推广方向已提及——但“下一步工作”的具体化程度可加强。v1 评阅书声称“无田间推广路径”措辞略夸，v2 改为“推广路径表述较弱，建议在 §6.3(3) 具体化”。

田间常见工况与论文适用范围的差距： - 田间秸秆通常是 0.5~2 kg/m² 多层覆盖 → 论文不适用 - 田间地表通常不平整且伴侵蚀 → 论文未做 - 田间降雨条件下流场变化剧烈 → 论文未做 - 田间径流流量常超过 18 L/min（试验最高） → 论文未覆盖

摘要与第六章创新点声称：“为坡面侵蚀模型关键水力参数获取、水土保持评价提供方法支撑”——声称强度与适用范围有落差。

5.4 第二章“统一原理”声称分析

第二章标题：《近地表水流流速测量原理及装置》

“原理”暗示一个统一的物理理论框架。

第二章实际章节结构： - §2.1 土壤水流通量矢量测量原理（地下） - §2.2 土壤水流通量矢量测量装置设计（地下） - §2.3 薄层水流流速测量原理（地表） - §2.4 薄层水流流速测量装置设计（地表） - §2.5 小结

§2.5 小结内容：先简述土壤水流通量方法(SPL/MDTD/Ratio)一段，再简述薄层水流方法（盐热示踪 + 质心法）一段——两段平行罗列，没把两套方法的物理共性说清。

两套方法的物理原理： - §2.1+§2.2（热脉冲）:基于热传导 + 强制对流热传递；探针施加热脉冲后温度衰减；解析模型为对数 + 指数函数 - §2.3+§2.4（盐热示踪）:基于物质对流 + 扩散迁移；示踪剂注入后浓度运移；解析方法为质心法（数值积分）

物理上完全无关——一个是“温度场对热源响应”,另一个是“浓度场对扰动响应”,求解原理完全不同。

如果章标题暗示有统一原理，应在 §2.5 把两套方法的共性说清楚（如“两套方法都遵循『让传感器对水流产生响应、再从响应信号反解出水流参数』的流程”）;否则章标题应改为“两类近地表水流流速测量方法及装置”,避免“统一原理”的暗示。

5.5 稳流段与有效测量段差距

§5.1 装置原文：“水槽前端 0.5 m 处设置挡板制成稳流槽，以增加水槽流速测量区供水的稳定性”

§5.1 数据处理原文：“考虑到混合示踪剂注入薄层水流初期会扰动流场且示踪剂处于明显的加速阶段，本研究舍弃第 1 与第 2 个探针间的流速数据，仅分析后 4 个探针所覆盖区段的平均流速”

几何关系： - 稳流槽：水槽前端 0~0.5 m - 5 个探针：坡长 0、0.6、1.2、1.8、2.4 m - 舍弃第 1~2 探针，仅用第 3~6 探针 → 有效测量起点为 1.2 m - 缓冲区：稳流段后 0~1.2 m,长度 = 1.2 − 0 = 1.2 m（含 0.5 m 稳流槽 + 0.7 m 加速段）

矛盾：稳流槽说“增加流速稳定性”,但实际数据处理需要再加 0.7 m 缓冲——意味着稳流段长度(0.5 m)不足以让流场充分发展。这与论文 §5.3.3 所声称的“试验在接近层流状态下进行”形成张力——若稳流不足，水流可能仍处于过渡 / 未发展状态，远非“接近层流”。

5.6 修改路径

第一组（命题表述级别）： - §5.3.3 明确“接近层流”的定量判据（如 Re < 1500 或 Fr < 1）,与试验覆盖范围(Re 625~1875)对照 - 重新表述“前提”——若所有试验都已在层流-紊流混合区，方法实际验证范围就是“层流-紊流混合区下的有限工况”,而非“接近层流” - 原文若想表达“试验设计上覆盖层流-紊流过渡区，盐热分层信号在过渡区仍可识别但适用性需进一步研究”,应明确写出这层意思

第二组（适用条件级别）： - 摘要 / §6.2 创新点：把适用工况（砂粒 1.3~3.2 倍、单层秸秆 0.5 倍或 1~2 倍）与不适用工况（多层秸秆、光滑 PVC 多数工况）在摘要中合并一句呈现，提高信息密度 - §6.3 下一步工作建议在已有“田间原位试验”基础上具体化：“该方法当前仅适用于室内定床特定工况，向田间推广需开展不平整地表、动态侵蚀过程、降雨条件下的验证试验”

第三组（统一原理级别）： - 第二章 §2.5 增加一段把两套方法的共性说清楚：例如“两套方法都遵循『让传感器对水流产生响应、再从响应信号反解出水流参数』的流程——前者用探针温升反解土壤热对流参数、后者用电导率与温度变化反解示踪剂质心运移速度” - 或修改章标题为“两类近地表水流流速测量方法及装置”,降低“统一原理”暗示

第四组（稳流段说明级别）： - §5.1 装置说明补：“实际有效测量段从坡长 1.2 m 起，原稳流段 0.5 m + 加速段 0.7 m,总长 1.2 m” - §5.5 讨论“稳流长度对盐热示踪流速可靠性的影响” - 若可能补做不同稳流段长度的对照试验

5.7 严重度判定

• 影响范围：第二章全章 + 第五章 §5.3~§5.5 + §1.4 技术路线图 + 第六章 §6.3
• 性质：命题表述不一致 + 适用条件呈现完整性 + 章节定位问题

意见 6：图表、文字与规范性校对问题（合并条目）

6.1 图 5-7 剖面流速示意（图 5-11 已撤回）

v1 → v2 修订：校验员视觉验证后：

• 图 5-7：右侧“剖面流速”显示一个上凸形剖面,arrows 在中部最长，向下递减到零（底部正确）,向上递减（顶部偏短）——这不是严格的开放渠道流剖面（开放渠道应底部 0、表面最大、单调递增）;也不是严格的封闭管道 Poiseuille 流剖面（Poiseuille 应对称抛物线，上下边界都为零）。实际形态介于两者之间。建议按真实开放渠道流物理事实修正
• 图 5-11（v1 误判，v2 撤回）:图 5-11 显示一个中部凹陷的剖面（顶部和底部都有较长 arrows,中部 arrows 较短）——这反映“水深超过秸秆厚度时上层水流和底层狭隙流并存”的物理事实，不是 Poiseuille 抛物线;符合论文要描述的“秸秆覆盖下的二元流场”。v1 把图 5-11 列为同图 5-7 类型的“Poiseuille 错画”是误判——撤回

两种剖面的物理对比（仅适用于图 5-7 的开放渠道情形）：

为什么图 5-7 的剖面错位影响最大： - 按论文画的上凸形剖面：上层和下层流速差不大（剖面对称）——支持“上下分层差异”的假设很弱 - 按真实开放渠道流剖面：上层（自由表面）流速远大于下层（底部）流速——反而强力支持盐热分层方法的物理基础

论文画错的示意图实际反过来削弱了自己方法的说服力。建议重画：按开放渠道流物理事实改为底部 0、表面最大的单调曲线（或对数律 / Manning 公式剖面）,并在图注中说明剖面模型。

6.2 图 5-4 (f) 与图 5-5 (i) 散落字母（v1 描述订正）

v1 → v2 修订：校验员高分辨率渲染核实后：

• 图 5-4 (c):未见明显“a”字母（v1 误描述为“a”,撤回）
• 图 5-4 (f)：确认有孤立字母（在图例上方区域）
• 图 5-5 (i)：确认有孤立“i”字母（在图右下区域）

修改建议：核查并清除图 5-4 (f) 和图 5-5 (i) 子图的散落字母。

6.3 表 5-1 至 5-4 变量编号跳号（v1 “符号重用”声明已撤回）

v1 → v2 修订：校验员独立核实表注后：

校验结论：X₃ 始终代表“示踪剂类型”,X₄ 始终代表“坡度”——符号本身使用一致，没有 v1 声称的“重用混淆”。撤回 v1 的“符号重用”声明。

保留的问题：变量编号跳号（方差分析表用 X₁/X₂/X₃,水力参数表用 X₁/X₂/X₄,跳过 X₃）确实不优雅，容易让读者困惑。

修改建议：统一变量编号方案，如改用 X_垫 / X_流 / X_坡 / X_示 等中文角标，或在每个表的注释里明示“本表中 X_i 含义与表 X-Y 相同 / 不同”。

6.4 §1.3 与 §1.5 研究内容编号不一致

§1.3 研究方法与内容(p.15)：明确列出四项研究内容 1. 单探针土壤水流通量大小测量方法 2. 五探针土壤水流通量方向测量方法 3. 盐热联合示踪测量薄层水流流速 4. 盐热联合示踪表征薄层水流二元垂向流速分布

§1.5 小结(p.16)：提“盐热联合示踪、单探针对数模型及五针探针矢量测量等创新方向”三项

§1.4 技术路线图（图 1-14, p.16）：也仅画两条主线（土壤水流通量 + 薄层水流）,第 4 项无独立位置

问题：四项与三项的层级关系是“第 3+4 项合并入『盐热联合示踪法』两个应用层次”——但论文未给合并说明，读者需自行推断。

修改建议：在 §1.3 末或 §1.5 开头加一句“四项研究内容归并为三个创新方向”过渡。

6.5 参考文献 Farbos de Luzan 2020（v1“疑似误引”已撤回）

v1 → v2 修订：校验员独立核实 §1.2.1 PIV 部分：

“PIV 对光照条件、示踪粒子分布和图像采集设备要求较高，仪器成本较大，且后期数据处理耗时(Farbos de Luzan et al., 2020)”

——该文献在 §1.2.1 有明确引用上下文，作为“PIV 复杂性 / 高成本”的示例引用。虽然主题（犬喉生物流体力学 PIV）与本论文（农业水文）距离较远，但有合法引用逻辑。v1 的“疑似误引”声明错误——撤回。

6.6 参考文献格式问题（保留）

修改建议：修订上述条目格式。

6.7 2026 年文献核实结果

校验结论：这三条都是合法引用，不需要在意见 6 中包含。v1 已正确处理，v2 维持。

6.8 第二章实物图比例尺缺失

普遍问题：

修改建议：在图中放置 1 cm / 1 mm 标尺，或在图注中明确尺寸。

6.9 SPL 英文缩写首注缺失

§1.3 SPL 首次出现：仅给中文“单探针对数法”+缩写 SPL,未注英文 Single Probe Logarithmic。

其它缩写同样问题：MDTD、DTD、PWM、ADV、ADCP、PIV 等首次出现时仅给中文全称 + 英文缩写，未给英文完整全称。

修改建议：§1.3 SPL 首次出现时改为“单探针对数法(Single Probe Logarithmic, SPL)”,其它缩写同样补全英文。

6.10 作者 2023 年前期工作未在 §1.2.2 引述（v2 新增）

事实：作者本人 2023 年农业工程学报已发表《盐热联合示踪表征薄层水流剖面流速分布》——这是论文方法的真正源头，但 §1.2.2 末段未引述，仅以“本研究拟将盐示踪与热示踪方法相结合”一句过渡。

修改建议：§1.2.2 末段或 §2.3 引述卢富运 等(2023),并明示“本研究在作者前期工作基础上系统化拓展为完整测量体系”。这条修订与意见 1 路径 B 配套使用——同时补齐文献综述空白与创新点物理基础的论证逻辑。

6.11 上机位界面截图字号偏小

问题：图 2-7（土壤水流通量上机位界面）和图 2-12（薄层水流上机位界面）截图字号偏小，曲线标签和数据列表需放大才能辨认。

修改建议：放大字号或加注重点参数。

6.12 表 3-4 列标题含混

问题：表 3-4 列标题“1.1a,b / 0.9a,b / a,1.1b / a,0.9b”含混——“1.1a, b”指 a 升 10% 同时 b 不变，还是 a 和 b 都升 10%?正文清楚（“对每个参数值进行独立轻微变动”）但表内不显式。

修改建议：列标题改为 “a→1.1a (b unchanged)” 或类似明确表达。

6.13 严重度判定

• 校对级合并条目
• 影响范围：分散在第一、二、三、五章及参考文献
• v2 撤回 v1 中的 5 项错判（图 5-11 Poiseuille / 图 5-4 (c) 字母 / X3-X4 符号重用 / Farbos de Luzan 误引 / 2026 年文献怀疑）
• 新增 6.10 作者 2023 年前期工作未引述（与意见 1 路径 B 配套）

致评阅意见接收方

本文档为评阅书的扩展佐证，主要服务于答辩秘书与申诉时回查。意见 2~5 是必须修改的实质性问题（涉及 SPL 精度呈现、五针法精度声明、SPL 模型理论基础与参数稳健性、盐热示踪命题表述）,意见 1 与意见 6 是建议修改的问题（盐热分层假设的直接观测验证补充 + 校对级问题）。作者按上述意见修改后建议参加答辩，答辩中应重点回应意见 2~4 涉及的精度声称与方法论证问题。

本评审立场不带学术诚信指控——论文工作量充实、攻读期间成果产出充分、§5.5/§6.1.2/§6.1.3/§6.3 已对方法局限做相对诚实的自承——评阅意见的诉求是把作者已经在结论位置说清的限制（如 SD、五针大小测量、田间推广）进一步前置到摘要，并把作者已经列为“下一步工作”的事（直接观测验证）前置一部分到论文修改阶段。