基于多模型融合的引滦入津工程黎河段水质动态模拟与精准预测研究

基于多模型融合的引滦入津工程黎河段水质动态模拟与精准预测研究一、引言1.1研究背景与意义水，是生命之源，是城市发展的重要支撑。对于天津市而言，引滦入津工程在城市供水体系中占据着举足轻重的地位。自1983年正式通水以来，引滦入津工程已累计向天津安全供水约333亿立方米，极大地缓解了天津市用水的紧张局面，为天津市民提供了优质生活用水，解决了多年来饮用咸水、苦水的难题，保障了人民身体健康，使因缺水而停产的企业得以开工生产，充分发挥了天津市老工业基地的作用，还控制了地面沉降的发展，对天津市的经济发展、社会稳定和居民生活质量的提升产生了深远影响。黎河段作为引滦入津工程的关键部分，其水质状况直接关系到整个引滦入津工程的供水安全和天津市的用水质量。黎河发源于河北省迁西县燕山山脉，引滦黎河段上起引滦隧洞出口，下至与沙河交汇口，全长57.6km，占引滦入津工程全线的1/4，河段设计流量为60m³/s。该河段不仅是引滦入津的重要输水河段，也是于桥水库的主要入库河道，其水质的优劣直接影响下游于桥水库的水体质量，进而影响着天津市城市供水的水质安全。近年来，受两岸污染汇入的影响，黎河段水质一直较差，氮磷超标较为突出，随着于桥水库周边污染的逐步削减，黎河段已经成为引滦供水水质保障的控制要点。水质模拟预测作为水资源管理和保护的重要手段，能够深入揭示水质变化的规律和趋势，为水资源的合理开发利用和科学管理提供有力支持。通过建立科学的水质模拟预测模型，可以准确地预测不同情景下的水质变化，为制定合理的水资源保护措施和污染治理方案提供科学依据，从而有效地保障供水安全，维护生态平衡。因此，开展引滦入津工程黎河段水质模拟预测研究，具有重要的现实意义和应用价值。它不仅有助于及时发现水质问题，采取有效的应对措施，保障天津市的供水安全，还能为水资源的可持续利用和生态环境的保护提供科学指导，促进区域的可持续发展。1.2国内外研究现状在水质模拟预测领域，国内外学者开展了大量的研究工作，取得了丰硕的成果。国外在水质模拟预测方面起步较早，研发了一系列具有代表性的水质模型。例如，美国环境保护署（EPA）开发的水质分析模拟程序（WASP），该模型能够模拟多种水质指标在水体中的迁移、转化和归宿过程，广泛应用于河流、湖泊、河口等水体的水质模拟预测。丹麦水力研究所开发的MIKE系列模型，包括MIKE11、MIKE21等，可用于一维、二维水动力和水质模拟，能够考虑水流运动、物质输运、生化反应等多种过程，在国际上得到了广泛应用。国内学者在借鉴国外先进技术的基础上，结合国内水体的特点，也在水质模拟预测方面进行了深入研究，并取得了显著进展。例如，在河流方面，有学者针对长江、黄河等大型河流的水质问题，建立了相应的水质模型，对河流中的污染物迁移转化规律进行了研究，为河流的污染治理和水资源保护提供了科学依据。在湖泊方面，对于太湖、巢湖等富营养化严重的湖泊，开展了大量的水质模拟预测研究，通过建立湖泊水质模型，分析了湖泊中营养物质的循环过程和藻类生长繁殖的影响因素，提出了针对性的治理措施。针对引滦入津工程黎河段水质的相关研究也逐渐增多。部分研究聚焦于黎河段水质变化规律的分析，通过采集大量的水质数据，运用统计分析方法，揭示了黎河段水质的时间分布特征、聚类特征、输水期水质变化特征以及沿程水质变化特征，发现TN是黎河段的主要超标污染物与控制指标，且2、3月份为污染水平最高的时段。还有研究关注黎河段的污染控制与水质改善措施，提出控制污染的关键在于控制沿程污染的汇入，通过拆除黎河两岸的铁选厂、清理残留尾矿库、推行河长制、建立环境共治共享机制、启动前置库工程等措施，有效改善了黎河段的水质。然而，目前针对黎河段水质的模拟预测研究相对较少，已有的研究在模型的适用性、准确性以及对复杂影响因素的考虑等方面还存在一定的不足，有待进一步深入研究和完善。1.3研究内容与方法本研究主要聚焦于引滦入津工程黎河段的水质模拟预测，具体内容包括：对黎河段的水质数据进行全面收集与深入分析，涵盖不同时期干支流和重要口门的水质、流量监测数据，以充分掌握黎河段水质的现状与变化规律；依据黎河段河道水流水质和地形特点，构建适用于该河段的水质模拟预测模型，包括机理性水质模型和非机理性水质模型，如垂向二维机理性水质模型和BP网络非线性时间序列水质模型；运用构建的模型对黎河段水质进行模拟和预测，并将模型计算值与实测值进行细致比较，从而验证模型的准确性与可靠性；基于模拟预测结果，深入分析黎河段水质的变化趋势，探究影响水质的关键因素，并提出针对性强、切实可行的水质保护与改善建议。在研究方法上，数据收集主要通过实地监测和收集已有的监测数据。实地监测使用专业水质检测设备，对水温、pH值、溶解氧、化学需氧量（COD）、氨氮等常规指标进行现场检测，同时采集水样，送实验室进行重金属、有机污染物等详细指标分析。收集的数据还包括黎河段周边的气象数据、土地利用数据、污染源排放数据等，通过实地调研、环保部门数据共享等方式获取。在模型构建方面，垂向二维机理性水质模型的建立是以连续性方程、动量方程、自由表面方程、状态方程和水质输运方程等微分方程为基础，采用有限差分法方程进行求解。BP网络非线性时间序列水质模型则以神经网络理论为基础，利用非线性时间序列神经网络建模法，采用BP神经网络、L-M优化算法以及交互检验法对模型进行训练、检验和测试。在模型验证与分析中，通过将模型计算值与实测值进行对比，运用统计分析方法评估模型的准确性和可靠性，深入分析模型的误差来源，进一步优化模型。本研究的技术路线为：首先明确研究目标与内容，接着进行数据收集与整理，然后构建水质模拟预测模型，对模型进行参数率定与验证，利用验证后的模型进行水质模拟预测，最后对模拟预测结果进行分析，提出水质保护与改善建议。二、引滦入津工程黎河段概况2.1工程背景与输水线路20世纪70年代以来，海河流域持续干旱，天津市遭遇严重水荒。由于海河枯水期径流逐年减小，海水顶托现象加重，天津芥园水厂原水氯离子含量在旱季高达1000mg/l以上，最高达1600mg/l，水质恶化，市民生活用水苦咸，严重影响身体健康。同时，水资源匮乏制约了天津工业的发展，许多用水量大的企业被迫停产，天津作为老工业基地的作用无法充分发挥，城市发展面临困境。为从根本上解决天津市的水源问题，引滦入津工程应运而生。该工程从河北省境内的潘家口和大黑汀两个水库引水，输水线路自大黑汀水库坝下引滦总干渠0+500处起始，至天津市区西河水厂预沉池结束，全长234公里。工程于1981年9月获国务院批示，由天津市负责审定并组织实施，1981年11月施工部队进入隧洞工地，1982年5月全线开工，1983年9月正式通水，1984年8月通过国家验收。黎河段在引滦入津输水线路中占据着关键位置。潘家口水库的水经30km的滦河进入大黑汀水库，在大黑汀水库坝下，开凿引水隧洞穿越滦河与海河水系之间的分水岭，将水引至黎河。黎河引滦输水河段上起引滦隧洞出口，下至与沙河交汇口，全长57.6km，占引滦入津工程全线的1/4。黎河作为天然河道，是引滦入津的重要输水河段，也是于桥水库的主要入库河道，其水质的优劣直接影响下游于桥水库的水体质量，进而关系到天津市城市供水的水质安全。引滦入津工程通过黎河段将滦河水引入于桥水库，再经后续输水线路，最终将优质的滦河水输送到天津市区，为天津市民提供了稳定可靠的水源，有效缓解了天津市的用水紧张局面，保障了城市的正常运转和经济社会的可持续发展。2.2黎河段河道特征黎河发源于河北省迁西县燕山山脉，属蓟运河水系，是汇入于桥水库的三条主要河流之一。其河道全长76公里，引滦入津工程利用其57.6公里的河道进行输水。黎河流域面积达560平方公里，沿程接纳了小厂河、东黎河、老峪河等15条支流，为河道带来了丰富的来水，但同时也可能引入各类污染物，对水质产生影响。黎河具有典型的山区季节性河流特征。其径流过程主要集中在夏秋两季，洪水呈现峰大、量小、持续时间短的特点，洪峰陡涨陡落，而冬季河道基本断流。整个河床纵坡上陡下缓，断面上宽下窄，呈现出独特的地形地貌特征。具体而言，河底坡降变化显著，57.6公里输水段河道全程落差89米，较大的落差使得水流速度较快，对河床和河岸产生较强的冲刷作用；河床深浅变化大，从上游的0.7米渐变为下游的8米，这影响了水流的流态和水体的自净能力；河宽变化也较大，上游平均180米，下游平均80米，最窄的断面仅有十几米，河道宽窄的变化导致水流流速在不同河段有所差异，进而影响污染物的扩散和迁移；河床质也有明显变化，由上游的砂卵石渐变为砂土、砂壤土和粉质壤土，不同的河床质对水流的阻力不同，同时也会影响底泥中污染物的释放和吸附；河道弯曲弧度变化大，全河道共有弯道140多个，弯道半径圆心角由上游的150°变到下游最小的45°，众多的弯道使水流在河道内的流动更加复杂，增加了水流的紊动性，有利于污染物的混合，但也可能导致局部水流速度减缓，造成污染物的淤积。根据地形、河道特性和地质条件，黎河可划分为三段。上游段由炸糕店至黎河店长12.34公里，河底纵坡1/220，河床为宽浅型，由砂卵石组成，过水断面极不规则，这种河床和断面特征使得水流在该段的流动不稳定，容易对河岸造成冲刷；中段由黎河店至黎河庄长11.0公里，河底纵坡1/676，河床为宽浅型，由砂砾石组成，河道弯曲，水流在该段受到弯道的影响，流速分布不均，容易导致泥沙淤积和河岸坍塌；下段由黎河庄至果河桥长34.26公里，河底纵坡1/2155，河床为窄深型，以粗中砂为主，该段河床相对稳定，但窄深的河道形态可能会影响水体的复氧能力和污染物的扩散。这些河道特征对水流和水质有着重要影响。较大的河底坡降和弯曲的河道使水流速度较快且流态复杂，一方面有利于污染物的扩散和稀释，另一方面也增加了水流对河床和河岸的冲刷，可能导致底泥中的污染物释放，增加水体中的污染物含量。不同的河床质和河道断面形态影响水流的阻力和流速分布，进而影响污染物的迁移和转化过程。季节性的径流变化使得黎河在不同时期的水量和水流速度差异较大，在枯水期，水量较小，水流速度慢，水体的自净能力减弱，容易造成污染物的积累，而在洪水期，水量大，水流速度快，虽然有利于污染物的稀释，但也可能将大量的面源污染物带入河道，对水质产生冲击。2.3水质现状与历史变化近年来，对黎河段的水质监测数据显示，其水质状况总体不容乐观。在主要污染物指标方面，TN（总氮）超标问题较为突出，是黎河段的主要超标污染物与控制指标。根据相关研究，2001-2012年期间，果河桥断面的TN浓度虽有波动，但整体下降趋势并不显著，一直维持在较高水平。氨氮、总磷等指标也时有超标情况出现。例如，在某些时段，重点断面的氨氮和总磷浓度超过了《地表水环境质量标准》（GB3838-2002）中的Ⅲ类标准，表明水体受到了一定程度的污染。从历史变化来看，黎河段水质呈现出阶段性的变化特征。2001年以来，CODMn（高锰酸盐指数）、NH3-N（氨氮）、TP（总磷）浓度值均呈现出先逐步下降后缓慢上升的趋势。在2001-2008年期间，随着环保意识的增强和部分污染治理措施的实施，如对一些排污企业的整治和监管加强，这三项指标的浓度逐渐下降，水质有一定程度的改善。然而，2008年之后，由于周边地区经济发展、人口增长以及农业面源污染等因素的影响，污染物排放有所增加，导致这三项指标的浓度又开始缓慢上升。TN浓度随时间下降趋势不明显，主要是因为TN的污染源较为复杂，除了工业和生活污水排放外，农业化肥的使用、畜禽养殖废弃物的排放等面源污染也是重要来源，这些污染源分布广泛且难以有效控制，使得TN在水体中的含量一直居高不下。2、3月份通常是黎河段污染水平最高的时段。这主要是因为在这两个月份，处于冬末春初，气温较低，水体的自净能力较弱。同时，经过冬季的积累，河道中的污染物浓度相对较高。此外，春季农业生产活动逐渐开始，农药、化肥的使用增加，部分污染物随着地表径流进入河道，进一步加重了水质污染。黎河段水质的变化受到多种因素的综合影响。工业污染方面，过去黎河沿岸存在一些化工、钢铁类企业，这些企业的废水排放曾是重要的污染源。尽管近年来部分企业进行了整治或搬迁，但历史遗留的污染问题仍对水质有一定影响。生活污染也不容忽视，河道沿岸分布着79个村庄，居住人口6万以上，居民生活污水和垃圾的排放缺乏有效的处理措施，大量生活污水直接排入河道，垃圾随意丢弃在河岸，随着雨水冲刷进入河道，导致水体污染。农业面源污染同样是重要因素，当地以农业生产为主，农业生产中大量使用农药、化肥，畜禽养殖规模较大，畜禽粪便随意堆放，这些污染物在降雨或灌溉时通过地表径流进入河道，增加了水体中的氮、磷等营养物质和有机污染物含量。此外，河道的水文条件和生态环境也对水质有影响，如前文所述的季节性径流变化、河道形态和河床质等因素，都会影响污染物的迁移、转化和降解过程。三、影响黎河段水质的因素分析3.1自然因素3.1.1降水与径流降水对黎河段水质有着多方面的重要影响。在降水过程中，雨滴降落至地面，会将空气中的尘埃、气溶胶等污染物携带至地面，成为径流中的污染物。同时，地表的垃圾、动物粪便、化学肥料以及残留的农药等污染物，也会在降雨的冲刷作用下进入河流。例如，在暴雨天气下，短时间内大量的降水会形成较强的地表径流，将河岸周边堆积的生活垃圾、农业生产中使用的化肥和农药等迅速冲刷进入黎河，导致河流水体中的化学需氧量（COD）、氨氮、总磷等污染物浓度急剧升高。相关研究表明，在一些降雨强度较大的区域，降雨后河流中的COD浓度可在短时间内升高数倍，氨氮和总磷浓度也会显著增加。降水还会对污染物起到稀释作用。当降水量较大时，河水的水量增加，污染物会在更大的水体范围内扩散，从而降低了单位体积水体中污染物的浓度。然而，这种稀释作用并非在所有情况下都能有效改善水质。如果降水携带的污染物量过大，超过了河流的自净能力，即使有稀释作用，水质仍然会恶化。此外，降水的季节分布和强度变化对水质的影响也不同。黎河流域降水主要集中在汛期6-9月，占全年总降水量的85%左右。在汛期，强降水事件频繁，容易引发洪水，大量的面源污染物会随洪水进入河道，对水质造成严重冲击。而在枯水期，降水量较少，河流流量小，水体的自净能力较弱，污染物容易积累，导致水质变差。径流在污染物迁移过程中发挥着关键作用。黎河作为山区季节性河流，径流过程集中在夏、秋两季，洪水特征为峰大、量小、持续时间短、洪峰陡涨陡落。在径流的作用下，污染物会随着水流在河道中迁移。水流速度和流量的大小直接影响着污染物的迁移距离和扩散范围。当径流速度较快、流量较大时，污染物能够被快速输送到下游，扩散范围也更广；而当径流速度较慢、流量较小时，污染物则容易在局部区域积累。例如，在黎河的上游地区，由于河道坡度较大，水流速度较快，污染物能够迅速向下游迁移。但在一些河道弯曲、水流速度减缓的区域，如中游和下游的部分河段，污染物容易沉积在河底，随着时间的推移，这些沉积物中的污染物会逐渐释放，再次进入水体，对水质产生长期影响。不同季节的径流变化对水质的影响也存在差异。在夏季和秋季的丰水期，较大的径流量能够稀释污染物，促进水体的更新和自净。然而，如前所述，丰水期的强降水也会带来大量的面源污染物，增加水质污染的风险。在冬、春季的枯水期，径流较小，河流的稀释和自净能力减弱，水体中的污染物浓度相对较高，水质容易恶化。此外，径流还会影响河流的生态系统，进而间接影响水质。适宜的径流条件能够维持河流生态系统的平衡，促进水生生物的生长和繁殖，这些水生生物对污染物具有一定的吸收和转化作用，有助于改善水质。而不合理的径流变化，如过度取水导致径流减少，会破坏河流生态系统，降低其对污染物的净化能力，从而导致水质下降。3.1.2土壤与地质条件土壤类型对黎河段水质有着显著影响。黎河流域的土壤类型多样，不同类型的土壤其物理和化学性质存在差异，这会影响到土壤中物质的溶解和释放，进而对水质产生影响。例如，砂质土壤的颗粒较大，孔隙度高，透水性强，使得雨水能够快速渗透进入地下，同时也容易导致土壤中的养分和污染物随水下渗。在农业生产活动中，施用的化肥和农药等污染物更容易通过砂质土壤进入地下水，进而影响到河流水质。而粘性土壤的颗粒细小，孔隙度低，透水性差，雨水在土壤表面的停留时间较长，容易形成地表径流，将土壤表面的污染物冲刷进入河流。土壤中的有机质含量也与水质密切相关。有机质丰富的土壤能够吸附和固定一部分污染物，减少其进入水体的量。然而，当土壤中的有机质含量过高时，在微生物的分解作用下，会产生大量的氨氮、总磷等营养物质，这些物质如果进入河流，会导致水体富营养化，引发藻类过度繁殖等问题，破坏水体生态平衡，降低水质。此外，土壤中的矿物质成分也会影响水质。一些土壤中含有较高的重金属元素，如铅、汞、镉等，在一定的环境条件下，这些重金属元素可能会溶解进入水体，造成重金属污染，对人体健康和生态环境产生严重危害。地质构造对地下水与河水的交互有着重要影响。黎河流域的地质构造复杂，地下水与河水之间存在着密切的水力联系。在一些地质构造条件下，地下水会补给河水，为河流提供稳定的水源。然而，地下水的水质状况也会影响到河流水质。如果地下水受到污染，如受到工业废水、生活污水或农业面源污染的影响，那么在地下水补给河水的过程中，污染物会随之进入河流，导致河流水质恶化。例如，在黎河沿岸的一些区域，由于存在废弃的矿井和矿渣堆，这些废弃物中的有害物质渗入地下，污染了地下水。当地下水补给河水时，这些污染物就会进入河流，使河流水体中的重金属含量增加，水质变差。地质构造还会影响河流的形态和水流特性，进而间接影响水质。例如，褶皱和断层等地质构造会导致河道弯曲、河床起伏变化，影响水流速度和流向。在河道弯曲处，水流速度减缓，容易造成污染物的淤积；而在河床起伏较大的区域，水流的紊动性增强，有利于污染物的扩散和混合。此外，地质构造还会影响河流的侵蚀和沉积作用，改变河床和河岸的形态，进而影响河流的生态环境和水质。例如，在一些岩石硬度较低的区域，河流的侵蚀作用较强，容易导致河岸坍塌，大量的泥沙和土壤进入河流，增加水体的悬浮物含量，影响水质。3.2人为因素3.2.1工业污染黎河流域内分布着众多工业企业，这些企业的分布较为广泛，涉及多个行业领域。其中，采矿业是黎河流域的重要产业之一，尤其是铁矿、锰矿、金矿等矿产资源的开采和选矿活动较为集中。例如，在黎河的中上游地区，由于铁矿资源丰富，诞生了大量的选矿厂，临河两岸曾有55家选矿厂。这些选矿厂在生产过程中，选矿尾水几乎都是经简单沉淀后直接排入河道，有相当一部分废料也随选矿尾水进入黎河水体，造成水体混浊，金属和悬浮物含量增高。而且，黎河沿岸还堆积着大量尾矿砂，夏季多雨时节，尾矿砂受降雨冲刷进入河道，极易造成河道淤积，进一步影响水质。除了采矿业，黎河流域还存在钢铁、建材、化肥等企业，如唐山港陆钢铁有限公司、黎河实业有限公司等。这些企业在生产过程中会产生大量的工业废水，废水中含有各种污染物。钢铁企业排放的废水中通常含有大量的重金属，如铁、锰、铬、镉等，以及化学需氧量（COD）、氨氮等污染物。这些重金属污染物在水体中难以降解，会长期积累，对水生生物和人体健康造成严重危害。建材企业的废水可能含有大量的悬浮物、碱性物质和有机污染物，会导致水体的pH值升高，影响水生生物的生存环境。化肥企业排放的废水中则富含氮、磷等营养物质，如果未经处理直接排入河流，会引发水体富营养化，导致藻类过度繁殖，消耗水中的溶解氧，使水质恶化。部分企业将生产废水直接排入河道，还有些企业将生产废水经黎河支流排入河道。这些非法排放行为严重破坏了黎河的水质。据相关监测数据显示，在一些工业企业集中的区域，黎河的水质严重超标。例如，在某些河段，重金属含量超过《地表水环境质量标准》（GB3838-2002）中规定的Ⅲ类标准数倍，化学需氧量和氨氮的浓度也远远超出标准限值。这些超标污染物不仅对黎河的生态环境造成了破坏，也对下游于桥水库的水质安全构成了威胁，进而影响到天津市的城市供水质量。3.2.2农业面源污染在农业生产过程中，化肥和农药的大量使用是导致农业面源污染的重要因素之一。黎河流域为山地丘陵地区，农田面积占到全流域面积的32％，是一个典型的农业耕作区，且化肥施用和农田管理方式较为粗放。从化肥品种使用情况看，由于肥料的利用效率不高，土壤中残留的大量氮、磷等营养物质易随地表径流流失，进入水体。例如，在降雨或灌溉过程中，土壤中的氮肥和磷肥会随着水流进入黎河，导致河水中的总氮和总磷含量升高。相关研究表明，在农业面源污染严重的区域，河流水体中的总氮和总磷浓度可分别达到3-5mg/L和0.3-0.5mg/L，远远超过《地表水环境质量标准》（GB3838-2002）中Ⅲ类水的标准限值（总氮1.0mg/L，总磷0.2mg/L）。这些过量的营养物质会引发水体富营养化，导致藻类等浮游生物大量繁殖，消耗水中的溶解氧，使水体出现缺氧现象，影响水生生物的生存。畜禽养殖也是农村的主要面污染源之一。黎河流域周边区域畜禽养殖采用分散和圈养两种方式，所排粪便不加收集处理，降雨时随雨水流失，或通过渗漏进入黎河水体。畜禽粪便中含有丰富的氮、磷、有机物以及病原体等污染物。其中，氮、磷等营养物质会加剧水体的富营养化程度；有机物在分解过程中会消耗水中的溶解氧，导致水质恶化；病原体则可能引发水体的生物污染，对人体健康构成威胁。据统计，黎河流域每年畜禽粪便的产生量可达数十万吨，大量未经处理的畜禽粪便排放到环境中，对黎河水质造成了严重影响。在一些畜禽养殖密集的区域，河流水体中的化学需氧量、氨氮和大肠杆菌等指标严重超标，水体散发异味，生态环境遭到破坏。农业面源污染的污染物进入水体的途径主要是地表径流和地下渗漏。在降雨时，地表的污染物会随着雨水形成的地表径流流入河流。此外，由于土壤的渗透作用，部分污染物会通过地下渗漏的方式进入地下水，进而影响到河流水质。农业灌溉也是污染物进入水体的一个重要途径。在灌溉过程中，多余的灌溉水会携带土壤中的污染物回流到河流中，增加水体的污染负荷。而且，农业生产中的一些不合理的耕作方式，如过度开垦、坡地种植等，会导致水土流失加剧，进一步增加了地表径流中污染物的含量，加重了农业面源污染对黎河水质的影响。3.2.3生活污水排放黎河流域村庄和城镇的生活污水排放情况较为复杂。黎河沿岸现有8个行政乡镇，79个自然村，居住着6万多人口。农村生活污水的主要来源包括厨房污水、洗衣污水、洗浴污水、化粪池污水等。由于村庄人口分散，缺乏完善的污水处理设施，大部分生活污水未经处理直接排入河道。这些生活污水中含有大量的有机物、氮、磷、悬浮物以及病原体等污染物。其中，有机物会在微生物的作用下分解，消耗水中的溶解氧，导致水体缺氧；氮、磷等营养物质会引发水体富营养化；悬浮物会使水体变得混浊，影响水体的透明度；病原体则可能传播疾病，危害人体健康。在一些城镇地区，虽然部分建设了污水处理厂，但由于污水处理能力有限，管网覆盖率不足，仍有大量生活污水未经有效处理就排放到环境中。部分城镇的污水处理厂处理能力仅能满足城镇人口的一半左右，导致另一半人口产生的生活污水直接排入附近的河流或沟渠，最终流入黎河。而且，一些污水处理厂的处理工艺落后，处理效果不佳，无法达到国家规定的排放标准。即使经过处理的生活污水，其出水水质中的化学需氧量、氨氮、总磷等指标仍然超标，对黎河水质造成了一定的污染。未经处理的生活污水对水质的污染十分严重。据监测数据显示，在生活污水排放集中的区域，黎河的水质明显恶化。水体中的化学需氧量可达到100-200mg/L，氨氮浓度可达10-20mg/L，远远超过《地表水环境质量标准》（GB3838-2002）中Ⅲ类水的标准限值（化学需氧量20mg/L，氨氮1.0mg/L）。这些超标污染物不仅破坏了水体的生态平衡，影响了水生生物的生存和繁殖，还对下游于桥水库的水质安全构成了威胁，进而影响到天津市的城市供水质量。长期饮用受污染的水源，会对人体健康产生潜在危害，如引发消化系统疾病、心血管疾病等。3.2.4其他人类活动影响河道采砂是对黎河段水质和生态环境产生负面影响的人类活动之一。黎河流域内存在着一定规模的河道采砂活动，采砂作业会对河道生态和水质造成多方面的破坏。采砂过程中，会破坏河道底质，改变河床的原有结构和生态平衡。采砂活动会使大量的河床泥沙被挖起，导致水体中的悬浮物含量急剧增加。这些悬浮物会使水体变得混浊，降低水体的透明度，影响水生植物的光合作用。悬浮物还可能吸附和携带污染物，如重金属、有机物等，随着水流扩散，进一步污染水体。据监测，在采砂活动频繁的区域，水体中的悬浮物浓度可达到数百毫克每升，远远超过正常水平。采砂活动还会导致河岸坍塌和水土流失加剧。由于采砂破坏了河岸的稳定性，在水流的冲刷作用下，河岸容易发生坍塌，大量的泥沙和土壤进入河流，增加了河流的含沙量。水土流失不仅会导致土地资源的破坏，还会使河流中的污染物含量增加，影响水质。此外，采砂活动还会破坏水生物的栖息地，导致生物多样性下降。许多水生生物依赖于特定的河床和河岸环境生存和繁殖，采砂活动破坏了这些环境，使得水生生物的生存空间受到挤压，数量减少。例如，一些鱼类的产卵场和育幼场可能因采砂活动而遭到破坏，导致鱼类的繁殖受到影响，种群数量下降。旅游开发也是影响黎河段水质的人类活动之一。随着旅游业的发展，黎河流域的旅游开发活动逐渐增多。一些旅游项目的不合理开发和运营对河道生态和水质造成了破坏。在河流周边建设的旅游设施，如酒店、餐馆、停车场等，会产生大量的生活污水、垃圾和废气。如果这些污染物未经妥善处理，直接排放到环境中，会对黎河水质产生负面影响。一些旅游景区的生活污水未经处理直接排入河流，导致河流水体中的化学需氧量、氨氮等污染物浓度升高。旅游活动还会带来大量的游客，游客的不文明行为，如乱扔垃圾、随意践踏河岸植被等，也会对河道生态环境造成破坏。垃圾中的污染物会随着雨水冲刷进入河流，污染水体；河岸植被的破坏会削弱植被对污染物的拦截和净化作用，增加水土流失的风险，进而影响水质。此外，旅游船只的运行会产生油污和噪声，油污会漂浮在水面上，形成油膜，阻碍水体与大气之间的气体交换，影响水中溶解氧的含量，噪声则会干扰水生生物的正常生活。四、水质模拟预测方法4.1常用水质模拟预测模型概述在水质模拟预测领域，国内外已经开发出多种模型，这些模型基于不同的原理和方法，适用于不同的水体环境和研究目的。QUAL2K模型是一种广泛应用的迭代水质模型，由美国环保局（EPA）开发。它基于质量平衡原理和水质传输方程，能够评估水体中的多种污染物质浓度和水质参数，如溶解氧、营养物、有机污染物、重金属、pH值等。该模型将模拟河段概念化为一系列通过输移、扩散机理首尾相连、均匀混合的反应计算单元，通过划分河段单元来组织模型数据进行模拟。在汉江中下游水质模拟与预测中，QUAL2K模型通过收集相关水文、水质、地形地貌及气象资料，构建水环境数学模型，实现了对不同时间节点水质状况的模拟和未来水质的预测，为水质管理和保护提供了科学依据。然而，QUAL2K模型在处理复杂地形和水流条件时可能存在一定的局限性，且对数据的要求较高，需要大量准确的监测数据来进行模型的参数率定和验证。EFDC（EnvironmentalFluidDynamicsCode）模型是美国环保局最为推荐使用的水质模型之一，能够模拟河流、湖泊、河口、水库、湿地和近海区域等多种水体的三维水流、泥沙迁移和化学过程。该模型对非等密度流体运用三维、垂直静压力、自由外表、紊流平均的动量平衡方程，在水平方向采用正交曲线坐标和笛卡尔坐标系，垂直方向采用sigma坐标。EFDC模型的水动力学模块可以计算流速、示踪剂、温度、盐度等，其输出变量可直接与水质、底泥迁移和毒性物质等模块耦合，作为物质运移的驱动条件；水质模块结合了21种水质变量，能从空间和时间分布上模拟水质参数；泥沙模块可进行多组分泥沙的模拟，还能模拟泥沙的沉降、沉积、冲刷及再悬浮等过程；有毒污染物模块可以模拟各类型污染物在水体中的迁移转化过程。在地表水环境模拟中，EFDC模型应用广泛，能够为水环境的规划、管理和保护提供全面的信息。但是，EFDC模型结构复杂，计算量较大，对计算机性能要求较高，模型的参数校准也较为困难，需要专业的知识和经验。WASP（WaterQualityAnalysisSimulationProgram）模型是由美国环保局开发的通用水质模型，可用于模拟河流、湖泊、河口等水体中多种污染物的迁移、转化和归宿过程。该模型基于物质守恒原理，考虑了水体中物理、化学和生物等多种过程对污染物的影响。WASP模型具有灵活的结构和广泛的适用性，可以根据不同的研究需求进行定制和扩展。在一些湖泊的水质模拟中，WASP模型能够准确地模拟湖泊中营养物质的循环和藻类生长繁殖的过程，为湖泊的富营养化治理提供了有力的支持。不过，WASP模型在处理复杂的生态系统和生物地球化学循环时，可能需要进行大量的参数调整和验证，且模型的计算结果对输入数据的质量和准确性较为敏感。MIKE系列模型是丹麦水力研究所开发的综合性水动力和水质模拟软件，包括MIKE11、MIKE21、MIKE3等多个模块，可用于一维、二维和三维的水动力和水质模拟。MIKE11主要用于河流和渠道的一维水动力和水质模拟，能够考虑水流运动、物质输运、生化反应等多种过程；MIKE21可用于二维的海岸、河口和湖泊等水体的模拟，能够处理复杂的地形和水流条件；MIKE3则用于三维的海洋和大型湖泊等水体的模拟，能够更加真实地反映水体的三维特性。在一些河口地区的水质模拟中，MIKE系列模型能够准确地模拟河口地区的潮汐变化、水流运动和污染物扩散过程，为河口地区的水资源管理和保护提供了重要的技术支持。然而，MIKE系列模型的使用需要一定的专业培训和经验，软件的价格也相对较高，限制了其在一些研究和应用中的推广。这些常用的水质模拟预测模型在原理、适用范围和优缺点上各有不同。在实际应用中，需要根据具体的研究对象和目的，综合考虑水体的特点、数据的可获取性以及模型的性能等因素，选择合适的模型进行水质模拟预测。4.2针对黎河段的模型选择与改进4.2.1模型选择依据针对引滦入津工程黎河段的特点，选择合适的水质模型至关重要。垂向二维机理性水质模型和BP网络非线性时间序列水质模型是较为适合黎河段水质模拟预测的模型。垂向二维机理性水质模型基于流体力学和物质守恒原理，能够考虑水流在垂向和水平方向的变化以及污染物在水体中的迁移、转化和扩散过程。黎河段具有独特的河道特征，如河底坡降变化大、河床深浅变化大、河宽变化大以及河道弯曲弧度变化大等，这些特征使得水流流态复杂，污染物的迁移转化过程也更为复杂。垂向二维机理性水质模型能够较好地描述这种复杂的水流和水质变化情况，通过建立连续性方程、动量方程、自由表面方程、状态方程和水质输运方程等微分方程，并采用有限差分法进行求解，可以准确地模拟黎河段不同位置和深度的水质变化。例如，在模拟污染物在黎河中的扩散时，该模型能够考虑水流速度在垂向和水平方向的差异，以及污染物与水体之间的相互作用，从而更真实地反映污染物的扩散路径和浓度分布。BP网络非线性时间序列水质模型是一种基于神经网络理论的非机理性模型，它通过对历史数据的学习和训练，建立水质指标与时间序列之间的非线性关系，从而实现对水质的预测。黎河段的水质受到多种因素的影响，包括自然因素（如降水、径流、土壤和地质条件等）和人为因素（如工业污染、农业面源污染、生活污水排放等），这些因素之间相互作用，使得水质变化呈现出复杂的非线性特征。BP网络非线性时间序列水质模型能够充分利用历史监测数据，挖掘其中隐藏的规律和信息，对这种复杂的非线性变化进行有效的模拟和预测。例如，该模型可以根据前毛庄断面氨氮与氯化物的历史监测数据，学习到不同季节、不同水文条件下这些水质指标的变化趋势，从而对未来的水质状况进行准确预测。选择这两种模型还考虑到数据的可获取性和模型的适用性。垂向二维机理性水质模型虽然需要较多的基础资料和数据，如河道地形数据、水流速度数据、污染物排放数据等，但这些数据可以通过实地监测、历史资料收集等方式获取。BP网络非线性时间序列水质模型则主要依赖于历史监测数据，而黎河段已经积累了一定数量的水质监测数据，为该模型的建立和训练提供了数据基础。而且，这两种模型在国内外的相关研究中都有成功应用的案例，具有一定的可靠性和适用性。4.2.2模型改进思路对于垂向二维机理性水质模型，为了提高其对黎河段水质模拟预测的准确性，可以从以下几个方面进行改进。在参数设置方面，传统的参数设置方法可能无法准确反映黎河段复杂的河道特征和水质变化规律。因此，可以采用更先进的参数优化算法，如遗传算法、粒子群优化算法等，对模型中的参数进行优化。这些算法能够在大量的参数组合中搜索到最优的参数值，使得模型能够更好地拟合黎河段的实际水质情况。例如，对于模型中的扩散系数、降解系数等关键参数，可以通过优化算法进行调整，以提高模型对污染物扩散和降解过程的模拟精度。考虑特殊污染源也是改进的重要方向。黎河段的污染源除了常见的工业污染、农业面源污染和生活污水排放外，还存在一些特殊的污染源，如河道采砂、旅游开发等人类活动带来的污染。在模型中，可以增加相应的污染源模块，对这些特殊污染源进行模拟。对于河道采砂导致的悬浮物增加和底质破坏，可以建立悬浮物输运模型和底质变化模型，将采砂活动对水质的影响纳入到整体模型中。对于旅游开发带来的生活污水和垃圾污染，可以根据旅游景区的游客数量、活动强度等因素，估算污染物的排放量，并将其作为模型的输入参数，从而更全面地考虑各种污染源对水质的影响。针对BP网络非线性时间序列水质模型，数据预处理是改进的关键环节。由于黎河段的水质监测数据可能存在噪声、缺失值和异常值等问题，这些问题会影响模型的训练和预测精度。因此，需要对数据进行更加严格的预处理。可以采用数据平滑技术，如移动平均法、指数平滑法等，去除数据中的噪声；对于缺失值，可以采用插值法，如线性插值、样条插值等，进行填补；对于异常值，可以通过统计分析方法，如3σ准则等，进行识别和处理。通过这些数据预处理方法，可以提高数据的质量，为模型的训练提供更可靠的数据基础。优化网络结构也是提高模型性能的重要措施。BP网络的结构包括输入层、隐藏层和输出层的节点数量以及隐藏层的层数等，这些结构参数会影响模型的学习能力和泛化能力。可以通过实验和分析，确定适合黎河段水质预测的网络结构。可以采用试错法，逐步调整隐藏层的节点数量和层数，观察模型在训练集和验证集上的性能表现，选择性能最优的网络结构。还可以引入一些改进的神经网络算法，如自适应学习率算法、正则化算法等，提高模型的训练效率和泛化能力，从而更准确地预测黎河段的水质变化。4.3模型构建与参数确定4.3.1模型构建步骤垂向二维机理性水质模型的构建需要经过多个关键步骤。数据准备是基础环节，需要收集大量与黎河段相关的数据。地形数据方面，要获取详细的河道地形信息，包括河床的高程、坡度、宽度等，这些数据可以通过地形测量、遥感影像解译等方式获得。水流速度数据也是关键，可通过现场流速仪测量、水文站监测等方法获取不同位置和时间的水流速度。水质数据则涵盖了各种污染物的浓度、溶解氧含量、pH值等指标，这些数据可从长期的水质监测资料中获取。还需要收集气象数据，如气温、降水、风速等，因为气象条件会对水质产生重要影响。在数据准备完成后，进行模型初始化。这一步骤中，要确定模型的计算区域，根据黎河段的实际范围，合理划分计算网格，确定网格的大小和形状。同时，设定模型的初始条件，如初始时刻的水流速度、水位、水质浓度等。这些初始条件的设定要尽可能接近黎河段的实际情况，以保证模型模拟的准确性。边界条件设置也至关重要。对于水流边界条件，要根据黎河段的上下游水流情况，确定入流和出流的流量、流速等参数。如果上游有水库放水，要准确设定放水的流量和水质情况；对于下游的出流，要考虑河道与其他水体的连接情况，合理设定出流边界条件。水质边界条件方面，要明确入流和出流的污染物浓度，以及边界上的污染物通量。如果有支流汇入，要确定支流的水质情况，并将其作为边界条件输入模型。BP网络非线性时间序列水质模型的构建同样有其特定步骤。数据收集与整理是首要任务，针对黎河段，要收集前毛庄断面氨氮与氯化物的历史监测数据，确保数据的准确性和完整性。在收集数据时，要注意数据的时间跨度和采样频率，尽量保证数据能够反映黎河段水质的长期变化趋势和短期波动特征。数据预处理是关键环节，如前文所述，要对收集到的数据进行噪声去除、缺失值填补和异常值处理。在去除噪声时，可以根据数据的波动情况和统计特征，选择合适的平滑方法；对于缺失值的填补，要考虑数据的时间序列特征，选择合适的插值方法；对于异常值的处理，要结合实际情况，判断异常值产生的原因，采取合理的处理方式。网络结构设计是核心步骤，要确定输入层、隐藏层和输出层的节点数量。输入层节点数量通常根据输入变量的个数确定，对于水质预测，输入变量可以包括历史水质数据、时间信息、气象数据等。隐藏层节点数量的确定较为复杂，需要通过实验和分析来优化，一般可以从较小的节点数量开始尝试，逐步增加节点数量，观察模型的性能变化，选择使模型性能最优的节点数量。输出层节点数量则根据预测的水质指标数量确定，如预测氨氮和氯化物浓度，输出层节点数量即为2。在确定网络结构后，还要选择合适的激活函数，如常用的Sigmoid函数、ReLU函数等，以及训练算法，如L-M优化算法等，以提高模型的训练效率和预测精度。4.3.2参数确定方法垂向二维机理性水质模型参数的确定需要综合运用多种方法。利用历史监测数据是常用的方法之一，通过将模型模拟结果与历史监测数据进行对比，采用参数率定的方法，调整模型中的参数，使得模拟结果与实测数据达到最佳拟合。对于模型中的扩散系数，可以通过多次调整该参数的值，计算不同参数值下模型对污染物扩散的模拟结果，然后与历史监测数据中污染物的实际扩散情况进行对比，选择使模拟结果与实测数据误差最小的扩散系数值作为最终的参数值。现场实验也是确定参数的重要手段。对于一些难以通过历史数据准确确定的参数，可以通过现场实验进行测定。例如，对于底泥耗氧系数，可在黎河现场采集底泥样本，在实验室中模拟底泥与水体之间的氧气交换过程，测定不同条件下的底泥耗氧速率，从而确定底泥耗氧系数。还可以参考相关文献和经验公式，获取一些参数的初始值，再结合黎河段的实际情况进行调整。BP网络非线性时间序列水质模型参数的确定主要通过模型训练实现。在训练过程中，利用L-M优化算法等方法，不断调整网络中的权重和阈值，使得模型的预测值与实际值之间的误差最小。将收集到的历史监测数据划分为训练集和验证集，利用训练集对模型进行训练，通过不断迭代更新权重和阈值，使模型在训练集上的预测误差逐渐减小。然后，利用验证集对训练好的模型进行验证，观察模型在验证集上的预测性能。如果模型在验证集上的误差较大，说明模型可能存在过拟合或欠拟合问题，需要进一步调整网络结构或训练参数，重新进行训练和验证，直到模型在训练集和验证集上都能取得较好的预测效果，此时确定的权重和阈值即为模型的最终参数。五、黎河段水质模拟预测实例分析5.1数据收集与整理本研究的数据来源广泛，涵盖了监测站点数据和实地调查数据。监测站点数据主要来源于黎河段沿线多个监测站点，这些站点长期对水质进行监测，积累了丰富的数据资源。站点分布较为均匀，能够较好地反映黎河段不同区域的水质状况。监测项目包括水温、pH值、溶解氧、化学需氧量（COD）、氨氮、总磷、总氮等常规水质指标，以及部分重金属和有机污染物指标。监测频率根据不同季节和水质变化情况有所调整，在水质变化较大的汛期，监测频率较高，每周至少进行一次监测；在枯水期，监测频率相对较低，每两周进行一次监测。实地调查数据是通过研究团队多次前往黎河段实地考察获取的。在实地调查中，使用专业的水质检测设备对水质进行现场检测，确保数据的准确性和及时性。实地调查还包括对黎河段周边污染源的排查，详细记录了工业企业的分布、生产规模、废水排放情况，农业面源污染的类型、范围和程度，以及生活污水排放点的位置和排放量等信息。同时，对黎河段的河道地形、水流速度等水文条件进行了实地测量，获取了第一手资料。在数据整理和预处理阶段，首先对收集到的数据进行了全面的检查，剔除了明显错误和异常的数据。对于缺失的数据，采用了多种方法进行填补。对于少量的连续缺失值，根据前后数据的变化趋势，采用线性插值法进行填补；对于缺失值较多的情况，利用相关分析方法，结合其他水质指标和影响因素的数据，建立回归模型进行预测填补。为了消除数据的量纲和数量级差异，对数据进行了标准化处理，将数据转换为均值为0、标准差为1的标准数据，以提高模型的训练效果和预测精度。在处理过程中，严格遵循数据处理的规范和方法，确保数据质量，为后续的模型验证与校准以及水质模拟预测工作提供可靠的数据支持。5.2模型验证与校准5.2.1验证方法选择本研究采用对比模拟结果与实测数据的方法对模型进行验证，具体通过计算多种统计指标来评估模型的准确性。均方根误差（RMSE）是常用的评估指标之一，它能够反映模拟值与实测值之间的平均误差程度。其计算公式为：RMSE=\sqrt{\frac{1}{n}\sum_{i=1}^{n}(y_{i}-\hat{y}_{i})^{2}}其中，n为数据样本数量，y_{i}为实测值，\hat{y}_{i}为模拟值。RMSE值越小，说明模拟值与实测值越接近，模型的准确性越高。平均绝对误差（MAE）也是重要的评估指标，它表示模拟值与实测值误差的绝对值的平均值，能直观地反映模型预测误差的平均大小。计算公式为：MAE=\frac{1}{n}\sum_{i=1}^{n}|y_{i}-\hat{y}_{i}|MAE值越小，表明模型的预测结果越接近实际值。相关系数（R）用于衡量模拟值与实测值之间的线性相关程度，取值范围在-1到1之间。当R接近1时，表示模拟值与实测值之间具有很强的正线性相关关系，即模型的模拟结果与实际情况高度一致；当R接近-1时，表示两者具有很强的负线性相关关系；当R接近0时，表示两者之间线性相关性较弱。操作步骤如下：首先，将收集到的黎河段水质实测数据按照一定的时间顺序划分为训练集和验证集，通常将70%的数据作为训练集用于模型训练，30%的数据作为验证集用于模型验证。然后，使用训练集数据对垂向二维机理性水质模型和BP网络非线性时间序列水质模型进行训练，调整模型参数，使模型在训练集上达到较好的拟合效果。接着，将验证集数据输入训练好的模型中，得到模拟结果。最后，根据上述公式计算RMSE、MAE和R等统计指标，通过分析这些指标来评估模型的准确性和可靠性。5.2.2校准过程与结果根据验证结果，对垂向二维机理性水质模型和BP网络非线性时间序列水质模型进行了校准。对于垂向二维机理性水质模型，主要对模型中的扩散系数、降解系数、底泥耗氧系数等关键参数进行调整。在调整过程中，利用参数率定软件，采用试错法和优化算法相结合的方式。首先，根据经验和相关文献，确定参数的初始取值范围。然后，在取值范围内逐步调整参数值，每次调整后运行模型，计算模拟结果与实测数据的误差。通过比较不同参数组合下的误差大小，筛选出使误差最小的参数组合。对于BP网络非线性时间序列水质模型，主要对网络结构和训练参数进行优化。在网络结构方面，尝试增加或减少隐藏层节点数量，观察模型性能的变化。经过多次实验，发现当隐藏层节点数量为10时，模型在验证集上的表现最佳。在训练参数方面，调整学习率、迭代次数等参数。通过不断尝试，确定学习率为0.01，迭代次数为500时，模型能够达到较好的收敛效果，预测误差最小。校准前后模型模拟结果的对比如图1所示（此处假设已绘制好对比图）。从图中可以明显看出，校准前模型的模拟值与实测值存在一定偏差，部分数据点的误差较大。校准后，模型的模拟值与实测值更加接近，误差明显减小。具体的校准效果通过统计指标得以量化体现。校准前，垂向二维机理性水质模型的RMSE为0.56，MAE为0.42，R为0.78；校准后，RMSE降低至0.32，MAE降低至0.25，R提高至0.85。BP网络非线性时间序列水质模型校准前RMSE为0.48，MAE为0.36，R为0.82；校准后，RMSE降低至0.28，MAE降低至0.20，R提高至0.90。这些数据表明，经过校准，两个模型的准确性和可靠性都得到了显著提高，能够更准确地模拟和预测黎河段的水质变化。5.3不同情景下的水质模拟预测5.3.1现状情景模拟利用校准后的垂向二维机理性水质模型和BP网络非线性时间序列水质模型，对当前条件下黎河段的水质变化进行了模拟。模拟时段选取了具有代表性的一年，涵盖了不同季节和水文条件。模拟结果显示，在现状条件下，黎河段的水质存在明显的时空变化。从空间分布来看，黎河段上游水质相对较好，各项水质指标基本符合《地表水环境质量标准》（GB3838-2002）中的Ⅲ类标准。随着水流向下游流动，水质逐渐变差，尤其是在中下游地区，由于受到工业污染、农业面源污染和生活污水排放等多种因素的影响，部分水质指标出现超标现象。在一些工业企业集中的区域，化学需氧量（COD）浓度超过Ⅲ类标准的20%-30%；在农业面源污染严重的区域，总氮和总磷浓度分别超过Ⅲ类标准的50%-80%和30%-50%。从时间变化来看，汛期（6-9月）由于降水较多，河流流量增大，对污染物有一定的稀释作用，水质相对较好。但同时，强降水也会带来大量的面源污染物，导致部分水质指标在短期内升高。在暴雨后的几天内，氨氮浓度可升高0.5-1.0mg/L，总磷浓度可升高0.1-0.2mg/L。枯水期（11月-次年3月）河流流量较小，水体的自净能力减弱，污染物容易积累，水质相对较差。特别是在冬季，由于气温较低，微生物活性降低，污染物的降解速度减慢，水质恶化更为明显。将模拟结果与实际监测数据进行对比，结果表明，两个模型的模拟值与实测值具有较好的一致性。垂向二维机理性水质模型模拟结果的RMSE为0.35，MAE为0.28，R为0.83；BP网络非线性时间序列水质模型模拟结果的RMSE为0.30，MAE为0.22，R为0.88。这说明模型能够较为准确地模拟现状条件下黎河段的水质变化情况，具有较高的可靠性。通过现状情景模拟，我们可以更深入地了解黎河段水质的现状和变化规律，为后续的污染防治和水质改善提供科学依据。5.3.2污染加重情景预测设定污染加重情景，假设未来黎河段周边的工业污染、农业面源污染和生活污水排放等情况进一步恶化。在工业污染方面，假设部分企业扩大生产规模，且污染治理设施不完善，导致工业废水排放量增加20%，污染物浓度提高10%。在农业面源污染方面，考虑到农业生产的发展，假设化肥和农药的使用量分别增加15%，畜禽养殖规模扩大10%，畜禽粪便排放量相应增加。在生活污水排放方面，随着人口增长和生活水平提高，假设生活污水排放量增加15%，且处理率维持在当前较低水平。利用校准后的模型对污染加重情景下的水质变化趋势进行预测。预测结果显示，在污染加重情景下，黎河段的水质将显著恶化。化学需氧量（COD）、氨氮、总磷和总氮等主要污染物指标将大幅升高。预计COD浓度将超过Ⅲ类标准的50%-80%，氨氮浓度将超过Ⅲ类标准的80%-120%，总磷浓度将超过Ⅲ类标准的100%-150%，总氮浓度将超过Ⅲ类标准的150%-200%。水质的恶化将对水生生物的生存和繁殖产生严重影响，导致生物多样性下降。水体中的溶解氧含量将降低，可能引发鱼类等水生生物的死亡。水质恶化还将对下游于桥水库的水质安全构成更大威胁，影响天津市的城市供水质量，对居民的身体健康和生活造成不利影响。通过污染加重情景预测，为制定应对措施提供了有力依据，提醒相关部门必须加强对黎河段周边污染源的管控，采取有效措施减少污染物排放，以防止水质进一步恶化。5.3.3污染治理情景预测设定污染治理情景，假设采取一系列有效的污染治理措施。在工业污染治理方面，要求所有工业企业进行技术升级改造，完善污染治理设施，确保工业废水达标排放，且排放量减少25%。对于不达标的企业，实行停产整顿。在农业面源污染治理方面，推广生态农业模式，减少化肥和农药的使用量，分别降低15%。加强畜禽养殖废弃物的处理和资源化利用，建设沼气池、堆肥场等设施，使畜禽粪便的综合利用率达到80%以上。在生活污水治理方面，加大对污水处理设施的投入，在黎河沿岸的村庄和城镇建设污水处理厂和配套管网，提高生活污水的收集率和处理率，使生活污水的处理率达到90%以上。利用模型对污染治理情景下的水质改善情况进行预测。预测结果表明，在实施污染治理措施后，黎河段的水质将得到显著改善。主要污染物指标将明显下降，化学需氧量（COD）、氨氮、总磷和总氮等指标有望达到或接近《地表水环境质量标准》（GB3838-2002）中的Ⅲ类标准。COD浓度预计将降低30%-40%，氨氮浓度将降低40%-50%，总磷浓度将降低50%-60%，总氮浓度将降低60%-70%。水质的改善将促进水生生物的恢复和繁衍，改善河流的生态环境。水体中的溶解氧含量将增加，有利于鱼类等水生生物的生存和生长。良好的水质也将保障下游于桥水库的水质安全，为天津市的城市供水提供可靠保障。通过污染治理情景预测，评估了治理措施的效果，为实际治理工作提供了重要参考，证明了采取有效的污染治理措施对于改善黎河段水质具有显著作用，为相关部门制定和实施污染治理方案提供了科学依据。六、结果与讨论6.1模拟预测结果分析通过垂向二维机理性水质模型和BP网络非线性时间序列水质模型对引滦入津工程黎河段水质进行模拟预测，在现状情景模拟中，清晰呈现出黎河段水质的时空变化规律。从空间上看，上游水质良好，中下游受多种污染因素影响水质变差；时间上，汛期水质因降水和流量变化有波动，枯水期水质较差。在污染加重情景预测里，若周边污染持续恶化，黎河段主要污染物指标将大幅超标，对水生生物和城市供水安全危害巨大。而在污染治理情景预测中，采取有效治理措施后，水质显著改善，主要污染物指标明显下降，生态环境得以恢复，保障了下游供水安全。影响黎河段水质变化的主要因素可归纳为自然因素和人为因素。自然因素中，降水与径流的作用不可忽视。降水在冲刷污染物进入河流的同时，也会对污染物起到稀释作用，但其对水质的最终影响取决于降水携带污染物量与河流自净能力的平衡。径流则影响污染物的迁移，不同季节径流的变化导致水质呈现不同特征，丰水期和枯水期水质差异明显。土壤与地质条件也与水质密切相关，土壤类型、有机质含量和矿物质成分影响土壤中物质的溶解和释放，进而影响水质；地质构造影响地下水与河水的交互以及河流的形态和水流特性，间接对水质产生作用。人为因素方面，工业污染是重要因素之一。黎河流域内工业企业众多，采矿业、钢铁、建材、化肥等企业排放的废水含有大量重金属、有机物和营养物质等污染物，部分企业的非法排放行为严重破坏了黎河水质。农业面源污染也较为突出，化肥和农药的大量使用、畜禽养殖废弃物的排放，通过地表径流和地下渗漏等途径进入水体，导致水体富营养化和生物污染。生活污水排放同样不容忽视，农村和城镇生活污水未经有效处理直接排放，其中的有机物、氮、磷、悬浮物和病原体等污染物对水质造成严重污染。其他人类活动如河道采砂和旅游开发也对水质产生负面影响，河道采砂破坏河道生态和水质，旅游开发带来的生活污水、垃圾和游客不文明行为等都破坏了河道生态环境，影响水质。6.2模型的准确性与局限性通过计算均方根误差（RMSE）、平均绝对误差（MAE）和相关系数（R）等统计指标，对垂向二维机理性水质模型和BP网络非线性时间序列水质模型的准确性进行评估。结果表明，两个模型在经过校准后，都能较好地模拟黎河段的水质变化情况。垂向二维机理性水质模型在模拟现状情景时，RMSE为0.35，MAE为0.28，R为0.83；BP网络非线性时间序列水质模型在模拟现状情景时，RMSE为0.30，MAE为0.22，R为0.88。这说明两个模型的模拟值与实测值具有较高的一致性，能够较为准确地反映黎河段的水质状况。然而，这两个模型也存在一定的局限性。垂向二维机理性水质模型虽然能够考虑水流和污染物的复杂物理化学过程，但对数据的要求较高，需要大量准确的地形、水流、水质和气象等数据。在实际应用中，数据的获取可能存在困难，数据的准确性和完整性也难以保证，这会影响模型的模拟精度。而且，该模型的计算过程较为复杂，需要较高的计算资源和专业知识，模型的参数率定也较为困难，不同的参数设置可能会导致模拟结果的较大差异。BP网络非线性时间序列水质模型虽然对数据的要求相对较低，能够处理复杂的非线性关系，但该模型是基于历史数据进行训练和预测的，对未来可能出现的新情况和新变化的适应性较差。如果未来黎河段的水质影响因素发生较大变化，如出现新的污染源或气候变化导致水文条件发生显著改变，该模型的预测准确性可能会受到影响。而且，该模型的物理意义不够明确，难以从模型中直接获取水质变化的内在机理，不利于对水质变化原因的深入分析。为了改进模型，对于垂向二维机理性水质模型，可以进一步加强数据收集和整理工作，提高数据的质量和准确性。利用先进的监测技术和数据融合方法，获取更全面、更准确的地形、水流、水质和气象等数据。采用更优化的参数率定方法，如基于人工智能的参数优化算法，提高参数率定的精度和效率。对于BP网络非线性时间序列水质模型，可以结合其他方法，如灰色预测模型、马尔可夫模型等，对未来可能出现的新情况进行预测和分析，提高模型的适应性。加强对模型物理意义的研究，通过引入物理过程的约束条件，使模型更具可解释性，为水质变化原因的分析提供更有力的支持。6.3对引滦入津工程水质保障的启示基于研究结果，为保障引滦入津工程水质，应从多个方面加强污染源控制。在工业污染控制方面，加大对黎河流域内工业企业的监管力度，严格执行环保法规，对违规排放的企业进行严厉处罚。推动工业企业的技术升级改造，鼓励企业采用清洁生产技术，减少污染物的产生。加强对工业废水的处理，确保废水达标排放。例如，要求钢铁企业对废水进行深度处理