1. 引言

长期以来，烘干称重法仍然是水流泥沙含量测定应用最广泛的方法。这种方法虽然精度高，但测量周期长、劳动强度大、检测过程繁琐。同时在水位较高流速较快时难以保证测验人员的安全。另一方面采样获得的是瞬时泥沙的浓度，在时间上无连续性，且受取样次数限制，不能反映泥沙实时变化过程，也难以捕捉含沙量峰值数据[1]。随着国家最严格的水资源管理制度的强制实行，利用在线自动化设备实现悬移质泥沙的在线实时连续监测是今后悬移质泥沙测量的必然发展方向，这给我们提供了进一步研究的机遇和挑战。

鉴于传统泥沙观测方法存在的问题，相关学者开展了大量观测泥沙新技术的研究，这些方法在特定的实验室条件或实验流域的应用中都取得了一定的成功[2] [3]。许多研究表明：基于红外光的悬移质泥沙在线监测技术是一种高效、安全和可靠的泥沙测验方法[4]-[6]，但在不同条件的水域中，仍受到悬移质泥沙的颜色、粒径及水体颜色和气泡等影响[7]。因此，在应用光学法仪器测验时需根据不同条件水域，将仪器示值和人工测沙数据进行比测分析率定后，才能把仪器示值用于悬移质含沙量计算中。本文以宜昌水文站为研究对象，论证TES-91泥沙在线监测系统在宜昌水文站的适用性，开展了泥沙监测仪精度、稳定性、可靠性比测试验，为该系统在长江流域的推广应用提供一定的经验支持。

2. TES-91泥沙在线监测系统

2.1. 系统组成

TES-91泥沙在线监测系统由泥沙监测仪、遥测终端设备(Remote Terminal Unit, RTU)、通讯网络、数据中心等组成。泥沙监测仪基于红外技术测定水体中的悬移质泥沙含量，服务器通过无线路由器采集监测水体的泥沙数据，通过服务器将相应数据传输至网络进行存储和后续分析处理，用户可在泥沙在线监测管理软件上进行远程管理和浏览。

2.2. 测量原理

系统运用了组合红外吸收散射光线法的原理。系统核心是一个红外光学传感器，可监测散射角为90˚~135˚之间的红外光散射信号。根据不同泥沙的反射率差异，利用传感器软件的逆投影成像技术，在不受色度影响下测定悬移质泥沙含量，直接通过数理模型转换并输出含沙量数据，实现在线连续水体中悬移质泥沙含量的监测。

2.3. 主要参数

系统主要技术参数见表1。

表1. 泥沙在线监测红外光学传感器主要技术参数

3. 测验河段概况

3.1. 测站基本情况

宜昌水文站位于长江中游干流上段(湖北省宜昌市滨江公园)，控制流域面积1,005,501 km²，占全流域面积55%，是国家基本水文站、国家重点水文站、三峡梯级水库总出库控制站。宜昌站测验项目有水位、水温、降水量、流量、悬移质输沙率及颗粒分析、沙质推移质、卵石推移质、床沙，长期以来为基本水文资料收集、长江中下游防汛抗旱、水资源管理和葛洲坝及三峡水利枢纽工程建设和运行服务。三峡水库蓄水运用后，宜昌站服务对象发生变化，增强了河道航运、水生态监测、水利枢纽工程运行调度服务要求，为水生态、水环境、水资源保护和利用提供监测资料。

枝城水文站位于长江中游干流上段(湖北省宜都市枝城镇滨江路)，集水面积1,024,100 km²，为国家基本水文站，是长江荆江入口重要控制站，三峡大坝泥沙出库验证站。监测项目包括降水量、水位、流量、悬沙、悬移质颗分、沙质推移质、卵石推移质、床沙等。该站为荆江河段及洞庭湖防汛抗旱、水资源合理调配、河道整治提供水文资料。

3.2. 河流特征

宜昌站测验河段长约3 km，较顺直，近几十年来河段河势未有大变化。断面上游约2.2 km处有一大弯，右岸为弯道凹岸，坡度较大。测验断面形态呈“U”型，左岸为砾卵石夹沙河床，低水时约有100 m的边滩裸露，高程42 m以上为混凝土护岸；右岸系基岩，地形陡峭。断面河床左边(起点距210~300 m)为沙滩，坡度较缓；河床中泓(起点距约300~700 m)为砾卵石夹沙质河床，近右岸(起点距约700 m以右)为礁板乱石组成。

枝城站测验河段在两弯道之间的顺直过渡段上，顺直段长度约3 km，略显上窄下宽状。断面上游700 m处有石矶，高水影响主泓摆动，河槽中高水位河宽1200~1400 m，属于宽浅型河流。测流断面河床组成为沙质和礁岩，冲淤变化不大，河床较为稳定。上游72 km有葛洲坝水利枢纽，上游20 km右岸有清江入汇，入汇口沿清江上溯10 km有高坝洲水利枢纽，测站下游2 km有枝城长江大桥，下游82 km右岸有1953年建成的荆江分洪区进水闸北闸，区间有松滋口、太平口分流入洞庭湖。枝城站径流量来自上游长江干流和清江。

图1. 垂线含沙量横向分布

图2. 垂线各测点含沙量分布

3.3. 测站断面含沙量横向分布

根据2017~2021年宜昌站选点法测验资料，点绘不同水位级下含沙量横向分布(图1)，各级水位级下垂线含沙量横向变化形态较为均匀，横向分布较为相似，一般主泓含沙量略大于其他区域含沙量，靠近左岸区域含沙量与断面平均含沙量具有较好的相关性，且波动性较为稳定；右岸含沙量最小；随着水位或流量变化，含沙量横向分布差异略有不同，当断面平均含沙量小于0.010 kg/m³时，横向分布差异逐渐缩小。

根据2017~2021年宜昌站选点法测验资料，统计近5年选点法13条取样垂线水深相对位置0.0、0.2、0.6、0.8、1.0的测点含沙量平均值，点绘选点法取样垂线沿垂向含沙量分布(图2)，各垂线水深相对位置测点含沙量具有从水面到河底逐渐增加的特点，但各测点含沙量差距较小，垂向分布较为均匀。

随着长江上游梯级水库修建和水土流失治理，三峡水库蓄水后上游来沙逐渐减小[8]，上游补给满足不了现有悬移质和床沙交换的需要，水流中悬浮的泥沙含量常处于不饱和状态，导致沿程河道床沙被清水冲刷至下游，河底上方悬浮的较大粒径泥沙含量逐渐减小，直至整个垂向含沙量变化减小。

4. TES-91现场比测试验

4.1. 仪器安装

系统于2020年11月6日安装在宜昌水文趸船内侧中部水下1.5 m处，在天然水体中进行测量。

4.2. 含沙量率定方法

1) 仪器稳定性率定分析：在监测仪探头旁同水深处，采用横式采样器独立取样，利用烘干法处理、计算测点含沙量，开展传感器处传统法测点含沙量与断面平均含沙量的比测，寻找相关关系，分析仪器监测数据的稳定性。

2) 建立仪器探头同位置的传统法测点含沙量与同步施测的断面平均含沙量相关关系，分析仪器安装位置处含沙量与断面含沙量之间是否有稳定关系，仪器安装位置是否合理，该处含沙量是否具有代表性。

3) 利用泥沙在线监测数据与断面平均含沙量建立相关关系，分析关系变化规律，判断在线监测数据过程推算断面平均含沙量的代表性及精度。

4.3. 比测率定资料选取

比测时间为2020年12月01日~2021年12月31日，人工测验含沙量样本160份。实测泥沙含沙量范围为0.002~0.235 kg/m³，本次比测试验含沙量级配均匀，代表性强。

系统收集的全年资料完整，未出现仪器故障、数据中断现象，通过数据采集面板实时进行含沙量采集。仪器每30 s自动采集一次。每5分钟将采集到的10组数据自动进行一次平滑处理，并存储1次[9] [10]。可通过登录系统，查询或导出已收集到的数据。

5. 比测数据分析

5.1. 稳定性分析

按照设计要求，在监测仪探头旁同深度处同步取样(采用1000 ml手提式横式采样器单独取样，利用烘干法处理和计算测点含沙量)，根据不同沙量级共取样46次，并利用与取样同时的在线监测数据建立相关关系，其相关性为0.99。建立传感器处传统法测点含沙量与断面平均含沙量相关关系，传统法测含沙量与断面平均含沙量关系值的相关度为0.9951，接近1.00，相关性较好，仪器安装位置较合理。

5.2. 相关性分析

当仪器感应器同位置测点含沙量与断面平均含沙量相关性很好情况下，对在线测沙仪器示值与实测断面含沙量建立相关关系(图3)，结果显示两者相关性显著。将在线测沙仪示值与实测断面含沙量建立的关系线节点录入整编程序中，按一元三点插值法推求在线断面平均含沙量。

图3. 2021年宜昌站TES-91自动监测含沙量与断面平均含沙量相关关系

图4. 2021年6~8月TES-91自动监测过程

对上述关系线断面含沙量大于0.020 kg/m³的测次进行检验，样本容量N = 71，三项检验均合格，随机不确定度为17.8%，系统误差为0.8%，满足《水文资料整编规范》中随机不确定度与系统误差分别不超过±18%和±2%的要求。因此，瞬时断面平均含沙量可由上述关系线和泥沙在线测沙仪器测量值推算而得。

5.3. 与下游站自动监测含沙量过程线比较

为检验监测数据的合理性，将宜昌站与下游已投产的枝城水文站TES-91自动监测含沙量过程进行比较。枝城站于2019年4月在荆江河段开展了比测试验，现已成功投产用于泥沙在线监测报汛[4]。由于2021年宜昌站较大的含沙量集中在5~8月(图4)，绘制该时段自动监测数据绘制成过程线图。其中7月16日、8月24~26日为人工配沙试验数据，7月28日为漂浮物遮挡造成数据失真。

剔除掉人工配沙试验和漂浮物遮挡影响的监测数据后，与枝城水文站TES-91自动监测数据进行比较。含沙量变化过程基本一致，说明泥沙在线监测系统测定的宜昌站含沙量变化过程合理性好，仪器用于宜昌站悬移质含沙量监测的性能稳定。

5.4. 仪器示值与实测含沙量洪水泥沙过程对比

将自动测沙监测数据通过仪器示值与实测断面含沙量建立的关系线进行换算，得到瞬时断面含沙量过程，并与宜昌站实测含沙量过程进行比较后可以看出，自动监测仪换算后的断面平均含沙量与实测断面含沙量关系较好，过程一致，吻合程度高，但最大沙峰峰值相差较大。

为进一步分析最大含沙量相差较大的原因，绘制自动监测含沙量过程线及实测断面含沙量测点关系图(图5)，河流断面含沙量浓度的变化过程可能受降雨的短期影响而产生短暂、剧烈的变化[11]，同一月份的最大含沙量能

图5. 2021年宜昌站TES-91自动监测含沙量过程线及实测含沙量测点

达到平日的数倍。这种短暂剧烈的变化过程表明人工测验含沙量的测次密集布置在流量发生较大变化的时候，具有较好的代表性，但由于人工测验的时间一般较为固定，难以捕捉到含沙量变化过程的转折点[12] [13]，比如2021年最大沙峰出现在7月21日04点，且沙峰消退速度特别快，最大沙峰峰顶处附近未布置实测断面含沙量测点。

5.5. 仪器示值与含沙量整编成果分析

根据2021年1月1日~12月31日在线测沙数据和已建立的关系线，计算断面平均含沙量，并利用该过程进行时段含沙量资料整编，生成在线逐日平均含沙量和逐日平均输沙率，并与传统人工测验含沙量取样方法计算的含沙量整编成果对照分析(图6和图7)，可以看出，两种逐日平均含沙量过程线较为吻合，变化趋势基本一致，但最大沙峰峰后段存在一定的过程变化偏差，与人工测验含沙量测次疏密有关。

图6. 自动监测逐日平均含沙量与实测逐日平均含沙量过程线

图7. 自动监测逐日平均输沙率与实测逐日平均输沙率过程线

将两种方法得到的含沙量各月特征值及误差进行统计，6~7月最大含沙量相对误差较大，分别为33.3%、101.3%，原因为最大沙峰出现在7月21日04点，且沙峰消退速度特别快，因安全考虑，人工测验含沙量未布置测次(断面含沙量测次布置时间在7月21日7:00~8:00)而在线监测基本实测到最大沙峰过程，造成两种方法含沙量沙峰数据出现较大误差；加上6月含沙量整体较小，绝对差为0.004 kg/m³，导致相对误差较大。统计2021年的年输沙量，实测值和自动监测值分别为1110 t和1100 t，两者相对误差为−0.9%，误差较小。综上，2021年宜昌站TES-91自动监测逐日平均含沙量与实测逐日平均含沙量过程线较为吻合，变化趋势基本一致，在线监测系统能真实反映含沙量变化过程。

5.6. 人工配沙试验

因2021年宜昌站断面含沙量较小，考虑到未来可能出现较大含沙量，在保证沙型情况下，人工增加单位体积内泥沙含量，检验TES-91测沙仪在大含沙量情况下的监测数据可靠性和代表性。

将宜昌站送水质泥沙室的悬移质泥沙样品烘干称重后进行配样，配样区间为0.010~3.00 kg/m³，按照沙量级均匀递增，配置96个沙样。

将测沙仪感应器放置在人工配沙的已知含沙量水样中进行测量，待数据稳定后，读取面板5组瞬时数据，取平均值作为不同沙量级下监测数据，最后将测得的96组数据建立关系线，并与今年已收集到传统法断面平均含沙量和在线监测数据对应关系，以及枝城站投产方案中传统法断面平均含沙量和在线监测数据对应关系共同进行比较(图8)。

图8. 三组关系对比图

如图，三组关系线(宜昌站实测断沙与TES-91关系线、枝城站实测断沙与TES-91关系线、人工配沙与TES-91关系线)较一致，其中宜昌站实测断沙与TES-91断沙关系与人工配样建立的关系线基本重合。人工配样建立的(含沙量范围为0.010~3.00kg/m³)关系中，当含沙量小于0.8 kg/m³时，关系线呈60˚角；当含沙量大于0.8 kg/m³时，关系线呈45˚角，当含沙量在0.4~0.8 kg/m³区间，关系线呈曲线，相关性较差。

对TES-91监测示值与人工配沙试验含沙量关系进行检验，三项检验均合格；随机不确定度为12.8%，系统误差为−1.0%，满足《水文资料整编规范》相关要求。

6. 结论

1) 经比测分析TES-91在线测系统监测值与同水深下传统人工测验含沙量取样得到的测点含沙量相关性，发现两者相关系数为0.99，表明TES-91在线测沙仪示值与实测断面含沙量相关性好。

2) 为检验监测数据的合理性，将宜昌站与下游已投产的枝城水文站TES-91自动监测含沙量过程进行比较，发现两站的含沙量变化过程基本一致，说明在线监测系统测定的宜昌站含沙量变化过程合理性好，仪器用于宜昌站悬移质含沙量监测的性能稳定。

3) 自动监测数据与宜昌断面平均含沙量关系在0.002~0.235 kg/m³范围内良好，通过三项检验和误差分析，误差满足《水文资料整编规范》相关要求，2021年宜昌站自动监测逐日平均含沙量与实测逐日平均含沙量过程线较为吻合，变化趋势基本一致，监测系统能真实反映含沙量变化过程。

4) 因2021年宜昌断面含沙量较小，采用人工配沙试验，比测试验的实测含沙量范围有待于进一步扩大，为实现长江流域的悬移质泥沙实时在线监测的推广应用奠定基础。

NOTES

作者简介：张释今(1990-)，女，工程师，主要从事水文测验、水环境监测和水资源研究，Email: 970734026@qq.com

^*通讯作者Email: 1211410268@qq.com

参考文献