黑龙江干流上游凌汛演进过程的多维度模拟与分析

黑龙江干流上游凌汛演进过程的多维度模拟与分析一、引言1.1研究背景与意义黑龙江作为我国北方的重要河流，其上游地区受高纬度严寒气候影响，每年11月份左右全线封冻，至次年4-5月份开江，封冻期长达150天以上。在开江期间，由于水力和热力因素的综合作用，凌汛灾害频发。据相关资料记载，建国后黑龙江发生冰坝的年份众多，如1950、1953、1956等年份，其中1971年、1985年、2000年等年份灾情尤为严重，1985年的特大冰坝更是给人民的生命财产造成了巨大损失，大兴安岭和黑河地区的6个县市41个村屯受灾，受灾人口达3.6万人，直接经济损失达2900万元。凌汛灾害的危害是多方面的。冰凌堵塞河道，会对水流产生阻力，进而引起江河水位明显上涨，可能导致堤防溃决，洪水泛滥成灾，冲毁沿岸的房屋、农田、道路等基础设施，威胁人民的生命安全。例如，冰坝形成后，冰块集中堆积在某一河段内，严重阻塞河道的过水面积，影响冰水下泄，使上游水位显著壅高，一旦堤防无法承受这种压力，就会出现决口，洪水会迅速淹没周边地区。凌汛还会对水利工程设施造成严重破坏，如损坏水坝、电站的引水渠、阻塞拦污栅等，影响水利设施的正常运行，进而影响电力供应和水资源的合理利用。凌汛期间，大量冰凌的流动和堆积还可能对桥梁、码头等跨河建筑物造成撞击，影响其结构稳定性，威胁交通运输安全。对黑龙江干流上游凌汛演进过程进行数值模拟及物理模拟研究具有极其重要的意义。通过模拟研究，可以深入了解凌汛的形成机制、发展过程和影响因素，从而提高凌汛预报的准确性和及时性。准确的凌汛预报能够为相关部门提前制定科学合理的防凌减灾措施提供依据，如提前组织人员撤离危险区域、加强堤防巡查和加固、合理调度水利工程设施等，有效减少凌汛灾害造成的损失。模拟研究有助于评估不同防凌措施的效果，为优化防凌方案提供科学参考。例如，通过模拟不同的破冰方案、水库调度方案等对凌汛演进的影响，选择最佳的防凌措施，提高防凌工作的效率和效果。对凌汛演进过程的研究还能为黑龙江流域的水资源管理和开发利用提供科学依据，保障流域的可持续发展，促进经济社会的稳定繁荣。1.2国内外研究现状1.2.1凌汛演进模拟技术发展凌汛演进模拟技术随着科技的进步不断发展，经历了从简单到复杂、从单一到多维度耦合的过程。早期的凌汛模拟主要基于经验公式和简单的水力学原理，对凌汛现象的描述较为粗略。例如，在20世纪中叶以前，研究人员多采用经验性的冰情指标来预测凌汛的发生，这些指标往往仅考虑了部分关键因素，如气温、流量等，对凌汛过程的模拟缺乏系统性和准确性。随着计算机技术的兴起，数值模拟方法逐渐应用于凌汛研究领域。20世纪70-80年代，一维水动力学模型开始被用于模拟河道水流和冰凌运动，这类模型将河道简化为一维线段，通过求解水流连续方程和动量方程，结合冰凌的生消和输移方程，能够初步模拟凌汛期间水位、流量和冰情的变化。但一维模型忽略了河道的横向变化和二维水流特性，对于复杂河道地形和水流条件下的凌汛模拟存在局限性。到了20世纪90年代以后，二维和三维数值模型得到了快速发展。二维模型考虑了河道的横向变化，能够更准确地模拟河道中水流和冰凌的分布情况，对于弯曲河道、分汊河道等复杂地形条件下的凌汛演进模拟具有明显优势。三维模型则进一步考虑了水流在垂向上的变化，能够更真实地反映冰凌在水体中的运动和相互作用，以及冰盖与水流之间的复杂耦合关系。例如，在一些大型江河的凌汛模拟中，三维模型能够准确模拟冰盖下水流的速度分布、温度场变化以及冰凌的堆积和融化过程，为凌汛灾害的预测和防治提供了更可靠的依据。除了数值模拟，物理模拟也是研究凌汛演进的重要手段。物理模拟通过在实验室中构建缩尺模型，模拟实际河道的地形、水流和冰情条件，直观地观察和测量凌汛过程中的各种物理现象。早期的物理模拟主要侧重于观察冰塞、冰坝的形成过程，随着技术的发展，现在的物理模拟能够更精确地控制实验条件，测量各种物理参数，如冰厚、流速、水位等，为数值模拟提供了验证和校准的数据。近年来，随着激光测量技术、粒子图像测速技术（PIV）等先进测量技术的应用，物理模拟的精度和可靠性得到了进一步提高。为了更全面地模拟凌汛演进过程，多维度耦合模拟技术逐渐成为研究热点。这种技术将数值模拟和物理模拟相结合，充分发挥两者的优势，通过数值模拟提供宏观的凌汛演进趋势，利用物理模拟验证和校准数值模型，并获取局部复杂区域的详细物理信息。例如，在一些大型水利工程的凌汛研究中，先利用数值模型进行整体河道的凌汛模拟，确定可能出现冰塞、冰坝的关键区域，然后在物理模型中对这些关键区域进行详细模拟，深入研究冰凌的堆积和融化机制，为工程的防凌设计和调度提供科学依据。随着大数据、人工智能等新兴技术的不断发展，凌汛演进模拟技术也在不断创新，未来有望实现更精准、更智能的凌汛模拟和预测。1.2.2黑龙江凌汛研究成果综述针对黑龙江凌汛，国内外学者开展了多方面的研究并取得了一定成果。在凌汛特点方面，研究表明黑龙江凌汛具有明显的季节性和区域性。其封冻期长，开江期间受多种因素影响，凌汛灾害频发。冰坝是黑龙江凌汛的主要表现形式之一，多发生在黑龙江上游的额尔古纳河入口处至呼玛河入口处约500km的河道上，其中洛古河、连崟（古城岛）等地最为频繁。冰坝出现时间一般在4月末至5月初，其形成后会导致上游水位显著壅高，形成冰坝洪水，对沿岸地区造成严重威胁。在凌汛成因研究上，学者们认为黑龙江凌汛主要受热力、动力和河道边界条件等多种因素的综合影响。热力因素方面，冬季气温的高低和变化影响冰盖的形成和厚度，开江期气温的快速回升是冰盖融化和凌汛发生的重要触发条件。动力因素中，封江水位、秋冬季降水量以及开江期的径流等对冰坝的形成起着关键作用。例如，封江水位高、秋冬季降水量大，会使河槽蓄水量增加，开江时水流动力增大，易形成冰坝。河道边界条件如河道的弯曲程度、宽窄变化以及河床的粗糙度等，会影响冰凌的流动和堆积，导致冰坝的形成。关于黑龙江凌汛灾害，众多研究分析了其对沿岸基础设施、农业生产和居民生活的严重影响。冰坝洪水可能冲毁堤防，淹没农田、房屋和道路，造成巨大的经济损失和人员伤亡。例如1985年黑龙江的特大冰坝，致使大兴安岭和黑河地区的6个县市41个村屯受灾，受灾人口达3.6万人，直接经济损失达2900万元。凌汛还会对水利工程设施、桥梁等造成损坏，影响其正常运行和安全。在凌汛模拟研究方面，虽然针对黑龙江凌汛的数值模拟和物理模拟研究不断推进，但仍存在一些挑战。数值模拟中，如何准确考虑复杂的热力、动力和河道边界条件，提高模型的精度和可靠性是关键问题。物理模拟则面临着如何在实验室条件下准确模拟黑龙江的实际冰情和河道条件，以及如何将实验结果有效应用于实际凌汛预测和防治的问题。尽管已有研究采用多种方法对黑龙江凌汛进行模拟，但由于黑龙江流域的复杂性和特殊性，模拟结果与实际情况仍存在一定偏差，需要进一步改进和完善模拟技术。1.3研究内容与方法1.3.1研究内容概述本研究旨在深入剖析黑龙江干流上游凌汛演进过程，综合运用数值模拟与物理模拟手段，全面揭示凌汛的形成机制、发展规律及其影响因素，为凌汛灾害的有效防治提供科学依据和技术支撑。在数值模拟方面，将构建适用于黑龙江干流上游的水动力学和冰动力学耦合模型。基于研究区域的地形地貌数据，精确刻画河道的几何形状和边界条件，利用高精度的地形测量资料，确保模型对河道地形的准确描述。考虑热力、动力等多种因素对凌汛的影响，引入冰盖生消、冰凌输移等关键过程的数学描述，如建立冰盖厚度与气温、水温之间的关系模型，以及冰凌在水流作用下的运动方程。通过数值模拟，预测不同工况下凌汛的发生时间、发展过程以及水位、流量、冰情等要素的时空变化，分析凌汛的形成机制和影响因素，如研究不同气温变化速率、流量大小以及河道弯曲程度对冰坝形成和凌汛演进的影响。物理模拟则是在实验室中搭建黑龙江干流上游的缩尺模型。依据相似性原理，严格控制模型的几何比尺、糙率比尺等相似条件，确保模型能够真实反映实际河道的水流和冰情特性。利用先进的测量技术，如激光水位测量仪、粒子图像测速仪（PIV）等，精确测量模型中水位、流速、冰厚等物理量的变化，实时监测冰塞、冰坝的形成和发展过程。通过物理模拟，直观展示凌汛现象，验证和校准数值模型，获取数值模拟难以得到的局部复杂区域的物理信息，如冰坝内部的冰块堆积结构和水流特性。将对数值模拟和物理模拟的结果进行对比验证。分析两者在凌汛演进过程模拟中的差异和一致性，评估模型的准确性和可靠性。针对模拟结果与实际情况的偏差，深入分析原因，改进模型参数和模拟方法，提高模拟结果的精度。综合数值模拟和物理模拟的优势，为黑龙江干流上游凌汛灾害的防治提供更科学、准确的决策依据，如通过模拟不同防凌措施下凌汛的演进情况，评估措施的效果，优化防凌方案。1.3.2研究方法介绍本研究采用多种研究方法相结合，以确保研究的全面性、准确性和可靠性。数值模拟方法是重要手段之一，通过建立水动力学和冰动力学耦合模型，对黑龙江干流上游凌汛演进过程进行数值求解。选用合适的数值算法，如有限体积法、有限差分法等，离散控制方程，实现对复杂水流和冰情的模拟。利用专业的数值模拟软件，如MIKEFLOOD、EFDC等，进行模型的搭建、计算和结果分析。通过设置不同的边界条件和初始条件，模拟不同工况下的凌汛过程，分析凌汛的形成和发展机制。物理模拟实验也是关键方法。在实验室中按照一定比例构建黑龙江干流上游的物理模型，模拟实际河道的地形、水流和冰情条件。严格控制实验条件，如水温、气温、流量等，使其与实际情况相似。运用先进的测量设备，对模型中的各种物理量进行实时测量和记录。通过改变实验参数，研究不同因素对凌汛的影响，如调整河道的弯曲程度、粗糙度，观察冰坝的形成和变化规律。物理模拟实验能够直观展示凌汛现象，为数值模拟提供验证和校准数据。为了更好地理解黑龙江干流上游凌汛的历史演变和规律，还将结合历史资料进行分析。收集整理黑龙江流域历年的凌汛观测数据、气象资料、水文数据等，分析凌汛的发生频率、强度、时间和空间分布特征。研究凌汛与气象、水文等因素之间的相关性，总结凌汛的形成条件和影响因素。通过对历史资料的分析，为数值模拟和物理模拟提供背景信息和验证依据，提高模拟结果的可靠性。二、黑龙江干流上游凌汛特性分析2.1流域概况黑龙江流域地理位置独特，其干流上游位于中国东北地区北部，地处北纬49°15′-53°33′，东经121°12′-127°19′之间。该区域是中国与俄罗斯的界河，其源头为额尔古纳河和石勒喀河，两河在漠河市以西的洛古河村附近汇合后始称黑龙江。黑龙江干流上游河段从洛古河村至黑河市，全长约900公里，河道蜿蜒曲折，水流湍急，对周边地区的生态环境和经济发展起着重要作用。从地形地貌来看，黑龙江干流上游流域地势总体西北高、东南低。西北部为大兴安岭山脉，其山势雄伟，海拔多在1000-1500米之间，山脉呈东北-西南走向，是黑龙江流域与蒙古高原的天然分界线。大兴安岭地区森林资源丰富，植被茂密，对涵养水源、保持水土起到了关键作用。东南部则为小兴安岭山脉，地势相对较为和缓，海拔一般在500-1000米左右，山脉呈西北-东南走向。小兴安岭地区同样森林广布，是重要的木材产区。流域内山间多河谷盆地，地势平坦，土壤肥沃，是当地农业和居民生活的集中区域。这些河谷盆地为河流的发育提供了良好的地形条件，使得黑龙江干流上游河道形态复杂，有许多弯曲和分汊的河段。黑龙江干流上游流域属于寒温带大陆性季风气候，冬季漫长而寒冷，夏季短促而温暖。冬季受西伯利亚冷高压控制，气温极低，平均气温在-20℃以下，极端最低气温可达-40℃以下。寒冷的气候使得河流封冻期长，一般从11月上旬开始封冻，直至次年4月下旬或5月上旬才开江，封冻期长达150天以上。在封冻期，冰盖厚度可达1-2米，冰盖表面平整，但在一些水流湍急或河道弯曲的区域，冰盖可能会出现裂缝和破碎现象。夏季受东南季风影响，降水相对集中，多以暴雨形式出现。年降水量一般在500-700毫米之间，其中6-8月降水量占全年的70%以上。夏季气温较为温和，平均气温在18-23℃之间。这种气候条件导致河流的径流量季节变化显著，冬季枯水期流量较小，而夏季丰水期流量较大。在凌汛期间，气温的快速回升和降水的增加会对凌汛的发展产生重要影响。在水文水系方面，黑龙江干流上游水系发达，支流众多。主要支流有额木尔河、呼玛河等，这些支流的汇入丰富了黑龙江的水量。额木尔河发源于大兴安岭伊勒呼里山北麓，全长约469公里，流域面积约1.69万平方公里。其河水清澈，水流较缓，在漠河市附近注入黑龙江。呼玛河是黑龙江上游的较大支流，发源于伊勒呼里山南侧，全长524公里，流域面积2.38万平方公里。呼玛河河道弯曲，河宽水浅，在呼玛县注入黑龙江。这些支流的水文特征与干流有所不同，它们的封冻期、开江期以及径流量变化等都会对黑龙江干流上游的凌汛产生影响。例如，支流的来水可能会改变干流的水流速度和冰情，在凌汛期间，支流与干流的冰块相互作用，可能会导致冰坝的形成和发展。黑龙江干流上游的河流水位、流量等水文要素也呈现出明显的季节性变化。在封冻期，河流水位较低，流量较小；而在开江期，随着气温升高和冰盖融化，河流水位迅速上涨，流量增大，容易引发凌汛灾害。2.2凌汛形成条件2.2.1河道特征黑龙江干流上游河道走向呈现东北-西南向，这种走向使得河流在冬季封冻和春季开河时，上下游的气温和冰情存在明显差异。由于高纬度地区气温较低，封冻时间早，而低纬度地区相对气温较高，封冻时间晚。在春季开河时，低纬度地区的冰先融化，水流携带冰块向下游流动，而此时高纬度地区的河道可能仍处于封冻状态，冰块容易在下游受阻堆积，形成冰坝，进而引发凌汛。河道的弯曲度对凌汛的影响也十分显著。黑龙江干流上游部分河段河道蜿蜒曲折，弯曲系数较大。例如，在洛古河至呼玛河段，河道弯曲系数可达1.5以上。河道弯曲处，水流速度减缓，冰凌容易在此聚集、堆积。当上游来冰量较大时，弯曲处的冰凌逐渐堆积形成冰塞，阻碍水流顺畅通过，导致水位上涨。随着冰塞的不断发展，最终可能形成冰坝，对河道安全构成严重威胁。研究表明，在相同的来冰条件下，河道弯曲度越大，冰塞、冰坝形成的概率越高，凌汛灾害的风险也越大。河道宽窄变化同样是影响凌汛的重要因素。该河段存在许多宽窄不一的区域，如在一些峡谷地段，河道宽度突然变窄，而在河漫滩地区，河道又相对较宽。当冰凌流至河道狭窄处时，过流断面减小，冰凌密度增大，容易相互挤压、碰撞，导致冰凌堆积。例如，在黑河附近的某河段，河道宽度在短距离内从500米骤减至200米，在凌汛期间，此处经常出现冰凌堵塞现象，致使上游水位急剧升高。而在河道较宽的区域，冰凌分散，水流速度相对较慢，也有利于冰凌的堆积和冰坝的形成。河道宽窄变化导致的水流条件改变，使得冰凌的运动和堆积规律变得复杂，增加了凌汛的不确定性和危害性。2.2.2气候因素气温变化是黑龙江干流上游凌汛形成的关键气候因素。在冬季，该地区受西伯利亚冷高压影响，气温急剧下降，当气温降至0℃以下时，河水开始结冰。从11月上旬开始，上游河段逐渐出现初冰，随着气温持续降低，冰盖逐渐增厚，至12月中下旬，大部分河段封冻，冰盖厚度可达1-1.5米。在封冻期，稳定的低温使得冰盖保持相对稳定状态。然而，在春季开江期，气温的变化对凌汛的发生起着决定性作用。当春季气温快速回升时，冰盖开始融化，大量冰块随水流向下游移动。如果气温回升速度过快，冰盖融化迅速，产生的冰块数量多、体积大，下游河道来不及排泄，就容易形成冰坝，引发凌汛。例如，在一些年份，3月下旬气温突然升高，日平均气温上升幅度可达5-8℃，导致冰盖快速融化，大量冰凌顺流而下，在下游狭窄河段堆积，形成冰坝，造成严重的凌汛灾害。降水分布也与凌汛的发生密切相关。黑龙江干流上游地区冬季降水较少，主要以降雪形式存在。冬季降雪量的多少会影响封江时河槽的蓄水量。如果冬季降雪量大，春季气温回升时，积雪融化形成的径流会增加河流水量，为凌汛提供了充足的水源。在凌汛期间，降水的分布也会影响凌汛的发展。如果在开江期，流域内出现降雨，尤其是暴雨，会进一步加大河流水量，使凌汛灾害更加严重。降水还会影响冰盖的稳定性，雨水的冲刷可能导致冰盖破裂，增加冰凌的数量和移动速度，加剧凌汛的危害。2.2.3水文条件流量、流速、水位等水文条件对黑龙江干流上游凌汛演进有着重要作用。流量是水流动力的来源，对冰凌的输移和堆积产生直接影响。在封江前，河流流量的大小决定了河槽的蓄水量。如果封江前流量较大，河槽蓄水量充足，在春季开江时，随着冰盖融化，大量的水和冰块一起向下游流动，增加了凌汛的动力和冰量。例如，在一些丰水年份，封江前流量达到500-800立方米每秒，开江时形成的凌汛规模较大，冰坝更容易形成，对沿岸地区的威胁也更大。流速与流量密切相关，流速的大小影响着冰凌的输移速度和卡冰堵塞的形式。在流速较快的河段，冰凌能够快速向下游移动，不易堆积。而在流速较慢的区域，冰凌容易停滞、聚集，进而形成冰塞、冰坝。在河道弯曲处或水流平缓的河漫滩地区，流速相对较小，冰凌容易在此堆积，导致凌汛的发生。水位是衡量凌汛灾害程度的重要指标。在凌汛期间，水位的变化直接反映了河道的行洪状况和冰凌对水流的阻塞程度。当冰凌堆积形成冰塞、冰坝时，河道过水能力减小，上游水位迅速上涨。水位的升高可能导致堤防漫溢、决口，引发洪水灾害，对沿岸的基础设施、农田和居民生活造成严重影响。在黑龙江干流上游的一些河段，历史上曾出现过因冰坝导致水位上涨超过3-5米的情况，致使大量房屋被淹，农田被冲毁，交通中断。水位的变化还会影响冰盖的稳定性，高水位可能导致冰盖破裂、坍塌，进一步加剧凌汛的危害。2.3凌汛特点2.3.1时空分布规律黑龙江干流上游凌汛在时间上呈现出明显的季节性和周期性。每年11月上旬左右，随着气温下降，河流开始出现初冰，进入流凌期。此时，河面上漂浮着大小不一的冰块，冰块随着水流向下游移动。随着气温进一步降低，11月中下旬至12月上旬，河流逐渐封冻，进入封冻期，冰盖逐渐增厚，冰厚一般可达1-2米。在稳定封冻期，冰盖相对稳定，但在一些特殊天气条件下，如强冷空气过境，可能会导致冰盖破裂、移动。次年4月下旬至5月上旬，气温回升，河流进入开河期，凌汛也随之发生。开河期可分为“文开河”和“武开河”两种情况，“文开河”时，冰盖逐渐融化，冰块缓慢随水流下泄，水位上涨较为平缓，凌汛危害相对较小；而“武开河”时，气温迅速回升，冰盖快速融化，大量冰块随水流迅猛下泄，容易在下游河道狭窄、弯曲处堆积，形成冰坝，引发严重凌汛灾害。从历史数据来看，黑龙江干流上游凌汛具有一定的周期性，约每隔3-5年就会出现一次较为严重的凌汛灾害。例如，在1985年、1990年、1995年等年份都发生了较大规模的凌汛，给沿岸地区带来了严重损失。在空间分布上，黑龙江干流上游凌汛主要集中在一些特定河段。从洛古河至呼玛河段，河道弯曲系数大，河宽变化明显，是凌汛的高发区域。在这些河段，冰凌容易在河道弯曲处和狭窄处堆积，形成冰塞和冰坝。洛古河附近的河道，由于地势起伏，水流速度变化较大，冰凌在此容易聚集，多次引发冰坝堵塞河道的情况。黑河附近的河段，由于河道宽度突然变窄，冰凌受阻堆积，导致水位急剧上升，是凌汛灾害的重点防范区域。此外，支流与干流的交汇处，由于水流条件复杂，冰凌相互作用，也容易出现凌汛现象。例如，呼玛河与黑龙江干流交汇处，在凌汛期间，支流来冰与干流冰块相互碰撞、堆积，增加了冰坝形成的风险。2.3.2冰坝、冰塞特征冰坝和冰塞是黑龙江干流上游凌汛过程中的重要现象，它们的形成机制较为复杂。冰坝的形成通常是由于上游来冰量较大，且下游河道存在阻碍冰凌下泄的因素。在开河期，当上游大量冰块随水流向下游移动时，遇到河道狭窄、弯曲或流速减缓的区域，冰块就会相互挤压、堆积，逐渐形成冰坝。冰坝的规模大小不一，小型冰坝长度可能在几十米到几百米之间，而大型冰坝长度可达数千米。冰坝的高度一般在1-5米左右，最高可达10米以上。1985年黑龙江干流上游形成的大型冰坝，长度超过5千米，高度达到8米，致使上游水位急剧壅高，造成了严重的洪水灾害。冰坝形成后，会对水位和流量产生显著影响。冰坝阻挡了水流的正常下泄，导致上游水位迅速上涨，水位涨幅可达3-10米。由于水位升高，河道过水断面减小，流速增大，进一步加剧了水流对冰坝和河岸的冲击。冰塞的形成主要是由于冰凌在河道中逐渐聚集、堆积，形成了类似堵塞的状态。在封冻期或开河初期，当水流速度较慢时，冰凌容易在河底或岸边堆积，随着时间的推移，堆积的冰凌逐渐增多，形成冰塞。冰塞的长度一般在几百米到数千米之间，厚度在0.5-2米左右。冰塞会导致河道过水能力下降，使上游水位逐渐抬高，影响河流通航和周边地区的用水安全。在冰塞发展过程中，水位会持续上升，当冰塞达到一定规模时，可能会引发冰坝的形成，进一步加剧凌汛灾害。冰塞还会改变河道内的水流结构，使水流流速分布不均，对河道生态环境产生不利影响。例如，冰塞会导致局部水流停滞，水温降低，影响水生生物的生存和繁殖。2.3.3灾害影响凌汛灾害给黑龙江干流上游沿岸地区带来了多方面的严重影响。洪水淹没是凌汛灾害的主要危害之一。当冰坝、冰塞形成后，上游水位急剧上涨，一旦超过堤防的设计标准，洪水就会漫溢，淹没周边的农田、房屋和基础设施。在1985年的特大凌汛中，黑龙江干流上游沿岸的多个村庄被洪水淹没，大量农田被冲毁，农作物绝收，许多居民被迫撤离家园。洪水还会浸泡房屋，导致房屋结构受损，甚至倒塌，严重威胁居民的生命安全。据统计，该次凌汛造成的房屋倒塌数量达到数千间，受灾人口众多。凌汛对堤防破坏也十分严重。冰坝、冰塞形成的巨大压力以及水流的冲击，会使堤防受到强烈的破坏。冰凌的撞击可能导致堤防的堤身出现裂缝、坍塌等情况，降低堤防的防洪能力。在高水位的浸泡下，堤防的基础可能被淘空，进一步加剧堤防的不稳定。一旦堤防决口，洪水将如猛兽般涌出，造成更大范围的灾害。在一些年份的凌汛中，由于堤防被破坏，洪水泛滥，导致交通中断，电力、通信等基础设施受损，给当地的生产生活带来极大的不便。凌汛灾害还会对当地的生态环境造成破坏，影响生物多样性。洪水淹没湿地和野生动物栖息地，可能导致一些珍稀物种的生存受到威胁。2.4历史典型凌汛案例分析1960年4月下旬至5月初，黑龙江上游江段出现了有水文记录以来的特大冰坝凌汛。此次凌汛的冰坝主要分布在额尔古纳河入口处至呼玛河入口处之间的河段。冰坝形成后，致使上游水位急剧壅高，水位涨幅超过5米，大量河水漫溢，淹没了沿岸的大片农田和村庄。从成因上看，当年冬季黑龙江上游地区降雪量较大，封江时河槽蓄水量充足，为凌汛提供了丰富的水源。春季气温回升迅速，冰盖快速融化，大量冰块随水流向下游移动。而此时下游河道尚未完全解冻，冰凌在下游受阻堆积，形成了规模巨大的冰坝。此次凌汛造成了严重的损失，沿岸多个村庄受灾，大量房屋被冲毁，农作物被淹没，交通和通信中断，直接经济损失达数百万元。由于当时的抢险救灾能力有限，受灾群众的生活受到了极大的影响。1985年，黑龙江干流上游再次发生特大冰坝凌汛，这也是黑龙江历史上最为严重的凌汛灾害之一。冰坝主要集中在洛古河至呼玛河段，冰坝长度超过5千米，高度达到8米。此次凌汛导致大兴安岭和黑河地区的6个县市41个村屯受灾，受灾人口达3.6万人。究其原因，1984年冬季降水偏多，封江水位较高，河槽蓄冰量较大。1985年春季，气温异常升高，开河速度加快，大量冰凌顺流而下。在下游河道狭窄、弯曲处，冰凌堆积形成冰坝。这次凌汛造成的直接经济损失达2900万元，大量农田被冲毁，房屋倒塌数千间。许多居民的财产损失殆尽，生活陷入困境。交通、电力、通信等基础设施也遭到严重破坏，对当地的经济发展和社会稳定造成了巨大冲击。三、凌汛演进过程数值模拟3.1数值模拟原理与方法3.1.1水动力模型选择本研究选用一维水动力模型和二维水动力模型相结合的方式来模拟黑龙江干流上游的水流运动。一维水动力模型将河道视为一维线段，主要通过求解圣维南方程组来描述水流的运动规律。圣维南方程组由连续性方程和动量方程组成，连续性方程反映了水流的质量守恒，即单位时间内流入控制体的水量等于流出控制体的水量与控制体内蓄水量的变化之和，其数学表达式为：\frac{\partialA}{\partialt}+\frac{\partialQ}{\partialx}=q其中，A为过水断面面积，t为时间，Q为流量，x为沿河道的距离，q为旁侧入流流量。动量方程则体现了水流的动量守恒，考虑了重力、摩擦力、惯性力等因素对水流的作用，其表达式为：\frac{\partialQ}{\partialt}+\frac{\partial}{\partialx}(\frac{Q^2}{A})+gA\frac{\partialh}{\partialx}+g\frac{Q|Q|}{CA^2R}=0其中，g为重力加速度，h为水位，C为谢才系数，R为水力半径。一维水动力模型在模拟河道整体水流趋势和长距离水流变化方面具有计算效率高、对数据要求相对较低的优点。在模拟黑龙江干流上游较长河段的水流时，能够快速得到水流的总体特征，如流量沿程变化、水位的大致分布等。但它忽略了河道的横向变化和二维水流特性，对于复杂河道地形和水流条件下的局部水流细节模拟不够准确。二维水动力模型考虑了河道的平面二维特性，通过求解二维浅水方程来模拟水流运动。二维浅水方程同样基于质量守恒和动量守恒原理，其连续性方程为：\frac{\partial\zeta}{\partialt}+\frac{\partial(hu)}{\partialx}+\frac{\partial(hv)}{\partialy}=0其中，\zeta为水位，h为水深，u、v分别为x、y方向的流速分量。动量方程在x和y方向分别为：\frac{\partial(hu)}{\partialt}+\frac{\partial(hu^2)}{\partialx}+\frac{\partial(huv)}{\partialy}=-gh\frac{\partial\zeta}{\partialx}-\frac{gu\sqrt{u^{2}+v^{2}}}{C^{2}h}+\frac{\tau_{sx}}{\rho}-\frac{\tau_{bx}}{\rho}\frac{\partial(hv)}{\partialt}+\frac{\partial(huv)}{\partialx}+\frac{\partial(hv^{2})}{\partialy}=-gh\frac{\partial\zeta}{\partialy}-\frac{gv\sqrt{u^{2}+v^{2}}}{C^{2}h}+\frac{\tau_{sy}}{\rho}-\frac{\tau_{by}}{\rho}其中，\rho为水的密度，\tau_{sx}、\tau_{sy}为表面切应力，\tau_{bx}、\tau_{by}为底部切应力。二维水动力模型能够更准确地模拟河道中水流的横向分布、流速变化以及复杂地形条件下的水流特性。在黑龙江干流上游的一些弯曲河道、分汊河道以及冰凌堆积区域，二维模型可以精确地描述水流的运动情况，为冰凌的输移和堆积模拟提供更准确的水流条件。例如，在模拟黑河附近河道的水流时，二维模型能够清晰地展现河道弯曲处的水流速度和流向变化，以及水流在不同区域的分流情况，这对于研究冰凌在该区域的运动和堆积具有重要意义。但二维模型计算量较大，对地形数据和边界条件的精度要求较高。在实际应用中，根据黑龙江干流上游的河道特点和研究需求，将一维和二维水动力模型相结合。对于河道的主要干流和较长河段，先使用一维模型进行整体的水流模拟，确定水流的总体趋势和关键参数。然后，针对河道中地形复杂、冰凌运动和堆积较为关键的区域，如弯曲河段、支流汇入处等，采用二维模型进行精细化模拟，以获取更详细的水流信息。通过这种方式，既提高了模拟效率，又保证了模拟结果的准确性。3.1.2冰情模拟方法冰盖生消过程的模拟是冰情模拟的重要部分。冰盖的生长主要受气温、水温等热力因素的影响。在寒冷的冬季，当气温持续低于0℃时，河水表面热量不断散失，开始结冰并逐渐形成冰盖。冰盖的生长过程可以通过能量平衡方程来描述，考虑冰层与大气、水体之间的热量交换，以及冰层内部的热传导等因素。假设冰盖的厚度为h_i，其生长速率与气温T_a、水温T_w、冰面热交换系数\alpha、冰的密度\rho_i、冰的比热容c_i以及冰的潜热L等参数有关。在稳定状态下，冰盖生长的能量平衡方程可以表示为：\rho_ic_ih_i\frac{\partialT_i}{\partialt}=\alpha(T_a-T_i)+\rho_wc_w(T_w-T_i)当冰盖达到一定厚度后，在气温回升时开始融化。冰盖的融化过程同样受气温和水温的影响，融化速率可以通过冰盖表面的热量吸收和冰层内部的热传导来计算。随着气温升高，冰盖表面吸收的热量增加，冰层逐渐融化变薄。冰塞、冰坝形成与发展的模拟是冰情模拟的关键环节。冰塞是由水内的碎冰在冰盖下堆积形成的，其形成与水流速度、冰凌浓度、河道地形等因素密切相关。在模拟冰塞形成时，考虑冰凌在水流作用下的运动轨迹和相互作用。当水流速度减小，冰凌的输移能力减弱，冰凌就会在河道中某些区域逐渐堆积，形成冰塞。冰塞的发展会导致河道过水断面减小，水流阻力增大，进而引起水位上升。通过建立冰塞的堆积模型，考虑冰凌的粒径分布、堆积密度等因素，模拟冰塞的发展过程。冰坝的形成通常是由于冰塞的进一步发展或者冰凌在河道狭窄、弯曲处的大量堆积。当冰凌堆积到一定程度，形成横跨河道的障碍物，阻碍水流正常通过，就形成了冰坝。冰坝形成后，上游水位会迅速壅高，对河道安全构成严重威胁。在模拟冰坝形成时，除了考虑冰凌的堆积过程外，还需要考虑冰坝的稳定性。冰坝的稳定性受到水流冲击力、冰凌之间的摩擦力以及冰坝与河道边界的相互作用等因素的影响。通过建立冰坝的力学模型，分析冰坝在不同水流和冰凌条件下的稳定性，预测冰坝的溃决风险。在模拟冰坝溃决时，考虑溃决瞬间的水流能量释放和冰凌的快速移动，以及由此引发的下游洪水过程。3.1.3模型耦合技术为了实现水动力模型与冰情模型的有效耦合，采用基于物理过程的耦合方式。在这种耦合方式中，水动力模型和冰情模型相互影响、相互作用。水动力模型为冰情模型提供水流速度、水位等关键信息，这些信息对于冰凌的输移和堆积起着决定性作用。冰情模型则为水动力模型提供冰盖厚度、冰塞和冰坝的位置及规模等信息，这些冰情参数会改变河道的边界条件和水流阻力，进而影响水流的运动。具体实现方式如下：在每个时间步长内，首先运行水动力模型，计算出水流的速度、水位等参数。然后，将这些水动力参数输入到冰情模型中，根据冰情模型的计算方法，模拟冰盖的生消、冰凌的输移和冰塞、冰坝的形成与发展。接着，根据冰情模型的模拟结果，更新河道的边界条件和水流阻力。将更新后的边界条件和水流阻力反馈到水动力模型中，重新计算水流的运动。通过这种反复迭代的方式，实现水动力模型与冰情模型的动态耦合，准确模拟凌汛演进过程。在模型耦合过程中，数据传递的准确性和及时性至关重要。为了确保数据的准确传递，采用统一的数据格式和接口标准。对水动力模型和冰情模型中的变量进行严格的定义和匹配，保证数据在两个模型之间的一致性。利用高效的数据传输算法，减少数据传输过程中的误差和时间延迟。还需要对模型耦合的稳定性进行分析和验证。通过设置不同的初始条件和边界条件，对耦合模型进行多次模拟，检验模拟结果的稳定性和可靠性。根据模拟结果的分析，对模型的参数和耦合方式进行调整和优化，提高耦合模型的模拟精度和稳定性。3.2模型构建与参数设置3.2.1研究区域概化将黑龙江干流上游从洛古河至黑河市的河段作为研究区域，该河段全长约900公里。考虑到研究重点和计算资源的限制，对研究区域进行合理简化。在平面上，将河道视为连续的曲线，忽略一些微小的支流和河汊，重点关注对凌汛演进有重要影响的主河道。对于河道边界，根据实际地形和河岸情况，将其概化为刚性边界，即假设河岸不会因水流和冰凌的作用而发生变形。采用非结构化网格对研究区域进行划分。非结构化网格能够更好地适应复杂的河道地形，在河道弯曲、宽窄变化较大的区域，可以灵活地调整网格的大小和形状，提高模拟的精度。在网格划分过程中，根据河道地形的复杂程度和研究的精度要求，对不同区域采用不同的网格密度。在河道主流区域，网格尺寸相对较大，以减少计算量；而在冰凌容易堆积的区域，如河道弯曲处、狭窄段以及支流汇入处等，网格尺寸则相对较小，以更精确地捕捉冰凌的运动和堆积过程。通过反复调试和验证，确定合适的网格尺寸，使得模拟结果既能满足精度要求，又能保证计算效率。3.2.2数据收集与处理收集研究区域的地形数据，主要通过地形测量和数字高程模型（DEM）获取。地形测量采用先进的测量仪器，如全站仪、GPS测量仪等，对河道的河床、河岸进行详细测量，获取高精度的地形信息。数字高程模型则通过卫星遥感数据或航空摄影测量数据处理得到，其分辨率能够满足模型对地形细节的要求。对收集到的地形数据进行预处理，去除数据中的噪声和异常值，填补缺失的数据，保证地形数据的准确性和完整性。将地形数据转换为模型所需的格式，以便在数值模拟中准确描述河道的地形特征。水文数据的收集包括水位、流量、流速等信息。水位数据通过设置在河道沿线的水位站进行观测，流量数据则通过水文站的流量测验获取。对于流速数据，采用流速仪进行测量。为了获取不同时段的水文数据，收集了多年的历史观测资料，包括凌汛期间和非凌汛期间的数据。对水文数据进行整理和分析，计算不同时段的平均水位、流量和流速等统计参数，为模型的初始条件和边界条件设置提供依据。对水文数据进行质量控制，检查数据的一致性和合理性，剔除明显错误的数据。气象数据的收集主要包括气温、降水、风速、风向等信息。气温数据通过气象站的温度传感器测量得到，降水数据通过雨量计记录。风速和风向数据则通过风速仪和风向标进行观测。收集研究区域及周边地区多个气象站的气象数据，以获取更全面的气象信息。对气象数据进行时空插值处理，将不同气象站的数据统一到模型的时间和空间尺度上。分析气象数据的变化规律，确定不同季节和时段的气象条件，为模型中冰盖生消和凌汛演进的模拟提供气象驱动数据。3.2.3参数率定与验证选择研究区域内有代表性的历史凌汛事件作为参数率定和验证的依据。从历史资料中筛选出多个典型凌汛年份，如1985年、2000年等，这些年份的凌汛过程具有不同的特点和规模，能够全面检验模型的性能。对于每个典型凌汛年份，收集凌汛期间详细的水位、流量、冰情等观测数据。采用试错法和优化算法相结合的方式进行参数率定。首先，根据经验和相关研究成果，初步设定模型中的参数值，如糙率、冰盖热传导系数、冰凌摩擦系数等。然后，运行模型，将模拟结果与观测数据进行对比。如果模拟结果与观测数据存在较大偏差，则通过试错法逐步调整参数值，再次运行模型，直到模拟结果与观测数据达到较好的拟合。为了提高参数率定的效率和准确性，引入优化算法，如遗传算法、粒子群优化算法等。这些算法能够在参数空间中自动搜索最优的参数组合，使得模拟结果与观测数据的误差最小。通过多次迭代计算，确定最优的参数值。在参数率定完成后，利用未参与率定的其他历史凌汛数据对模型进行验证。将验证数据输入到率定好的模型中，运行模型得到模拟结果。采用多种评价指标对模拟结果进行评估，如均方根误差（RMSE）、平均绝对误差（MAE）、相关系数（R）等。均方根误差能够反映模拟值与观测值之间的平均误差程度，其计算公式为：RMSE=\sqrt{\frac{1}{n}\sum_{i=1}^{n}(y_{i}-\hat{y}_{i})^2}其中，n为数据点的数量，y_{i}为观测值，\hat{y}_{i}为模拟值。平均绝对误差则衡量模拟值与观测值之间绝对误差的平均值，计算公式为：MAE=\frac{1}{n}\sum_{i=1}^{n}|y_{i}-\hat{y}_{i}|相关系数用于评估模拟值与观测值之间的线性相关性，取值范围在-1到1之间，越接近1表示相关性越强。通过对评价指标的分析，判断模型的准确性和可靠性。如果模拟结果的评价指标满足预定的精度要求，说明模型能够较好地模拟黑龙江干流上游凌汛演进过程；如果评价指标不理想，则需要进一步分析原因，对模型进行调整和改进，如重新检查数据质量、调整参数值或改进模型结构等。3.3模拟结果与分析3.3.1凌汛演进过程展示通过数值模拟，得到了黑龙江干流上游凌汛演进过程中水位、流量和冰情的动态变化。图1展示了凌汛期间典型断面的水位变化过程。从图中可以看出，在凌汛初期，水位较为稳定，随着冰凌的逐渐堆积，水位开始缓慢上升。当冰坝形成后，水位迅速上涨，达到峰值，此时水位涨幅可达3-5米。随后，随着冰坝的逐渐消融或溃决，水位又逐渐下降。在不同年份和不同河段，水位变化的幅度和速率存在差异，这与冰坝的规模、形成位置以及水流条件等因素有关。图2呈现了流量在凌汛期间的变化情况。在封冻期，流量相对较小且较为稳定。进入凌汛期，随着冰盖融化和上游来水增加，流量逐渐增大。当冰坝形成时，由于河道过水能力减小，流量会出现短暂的下降。而在冰坝溃决后，大量的水和冰凌迅速下泄，流量急剧增大，形成洪峰。洪峰流量的大小与冰坝的规模、上游蓄水量以及河道的泄洪能力等因素密切相关。不同年份的凌汛，流量变化过程也有所不同，一些年份的凌汛可能出现多次洪峰，这是由于冰坝的多次形成和溃决导致的。在冰情变化方面，图3展示了冰凌密度在河道中的分布情况。在凌汛初期，冰凌主要集中在河流的上游和岸边，随着凌汛的发展，冰凌逐渐向下游移动，并在河道狭窄、弯曲处堆积，导致冰凌密度增大。在冰坝形成区域，冰凌密度达到最大值，冰块相互挤压、堆积，形成了紧密的冰体结构。随着凌汛的结束，冰凌逐渐融化，冰凌密度减小。通过对冰凌密度分布的模拟，可以清晰地了解冰凌的输移和堆积规律，为凌汛灾害的防治提供重要依据。3.3.2关键指标分析凌峰流量是衡量凌汛规模和危害程度的重要指标之一。通过数值模拟，对不同年份和不同工况下的凌峰流量进行了计算和分析。结果表明，凌峰流量与多种因素密切相关。封江前的河槽蓄水量对凌峰流量有着重要影响。如果封江前河槽蓄水量大，在春季开江时，随着冰盖融化和水流的作用，更多的水和冰凌会向下游流动，从而导致凌峰流量增大。例如，在一些丰水年份，封江前河槽蓄水量充足，凌峰流量可达到1000-1500立方米每秒。气温回升速度也是影响凌峰流量的关键因素。气温回升速度快，冰盖融化迅速，产生的冰块数量多，会增加凌汛的动力和冰量，进而导致凌峰流量增大。当春季气温在短时间内快速升高时，冰盖大量融化，凌峰流量可能会超过正常年份的2-3倍。最高水位同样是评估凌汛灾害的关键指标。在凌汛期间，冰坝的形成会导致河道水位急剧上升，最高水位的出现往往伴随着冰坝的形成和发展。通过模拟分析发现，冰坝的规模和位置对最高水位有着直接影响。大型冰坝由于其阻挡水流的能力更强，会使上游水位壅高得更多，从而导致最高水位更高。冰坝形成在河道狭窄、弯曲处，由于水流阻力增大，也会使水位升高更为明显。河道的行洪能力也会影响最高水位。如果河道行洪能力不足，无法及时排泄凌汛期间的水量，就会导致水位持续上升，最高水位相应增加。冰坝位置的模拟对于凌汛灾害的防范具有重要意义。通过数值模拟，准确预测了冰坝可能形成的位置。结果显示，冰坝主要出现在河道弯曲系数大、河宽变化明显的河段。洛古河至呼玛河段，由于河道蜿蜒曲折，河宽变化较大，冰凌容易在此堆积形成冰坝。在这些河段，冰凌在水流的作用下，受到河道边界的阻挡，逐渐聚集、堆积，最终形成冰坝。支流与干流的交汇处也是冰坝容易形成的区域。由于支流和干流的水流条件不同，冰凌在交汇处相互碰撞、堆积，增加了冰坝形成的风险。通过对冰坝位置的模拟，可以提前确定凌汛灾害的重点防范区域，为采取有效的防凌措施提供依据。3.3.3敏感性分析为了研究不同参数对模拟结果的影响，进行了敏感性分析。首先考虑了糙率对模拟结果的影响。糙率反映了河道边界对水流的阻力大小。通过改变糙率值，模拟凌汛演进过程，并分析水位、流量等关键指标的变化。结果表明，糙率的变化对水位和流量有较为显著的影响。当糙率增大时，河道水流阻力增大，流速减小，导致水位升高，流量减小。在糙率增加10%的情况下，水位可能会升高0.5-1米，流量则会减小10-20%。这是因为糙率增大使得水流与河道边界的摩擦力增大，水流能量损失增加，从而影响了水流的运动。冰盖热传导系数也是影响模拟结果的重要参数。冰盖热传导系数决定了冰盖与大气、水体之间的热量交换速率。通过调整冰盖热传导系数，观察冰盖生消过程和凌汛演进的变化。当冰盖热传导系数增大时，冰盖与外界的热量交换加快，冰盖融化速度增加。在冰盖热传导系数增大20%的情况下，冰盖融化时间可能会提前5-10天，这会导致凌汛提前发生，且凌汛期间的冰量减少，从而对凌汛的规模和危害程度产生影响。冰盖热传导系数的变化还会影响冰盖下水流的温度和流速分布，进而影响冰凌的输移和堆积。冰凌摩擦系数对模拟结果也有一定的影响。冰凌摩擦系数反映了冰凌之间以及冰凌与河道边界之间的摩擦力大小。当冰凌摩擦系数增大时，冰凌的运动受到更大的阻力，冰凌的输移速度减小，更容易在河道中堆积。在冰凌摩擦系数增大15%的情况下，冰凌在河道弯曲处和狭窄处的堆积量可能会增加15-25%，这会导致冰坝更容易形成，且冰坝的规模可能会增大，从而加剧凌汛灾害的危害程度。通过敏感性分析，可以明确不同参数对凌汛模拟结果的影响程度，为模型参数的优化和调整提供依据，提高模拟结果的准确性和可靠性。四、凌汛演进过程物理模拟4.1物理模拟实验设计4.1.1实验目的与方案物理模拟实验旨在通过在实验室可控条件下，再现黑龙江干流上游凌汛演进过程，直观展示冰塞、冰坝的形成和发展，获取关键物理参数，为数值模拟提供验证和补充，深入揭示凌汛形成机制和影响因素。实验方案以黑龙江干流上游洛古河至呼玛河段为原型，构建缩尺物理模型。考虑到实验场地和设备的限制，选取合适的几何比尺，确保模型能够准确反映原型河道的地形和水流特征。实验分为不同工况进行，分别改变流量、水温、冰量等关键因素，模拟不同条件下的凌汛演进过程。设置高流量、低流量工况，研究流量对冰凌输移和冰坝形成的影响；改变水温，模拟不同气温条件下冰盖的融化和冰凌的生消过程；调整冰量，观察冰凌堆积和冰坝规模的变化。通过多工况实验，全面分析各因素对凌汛的影响规律。4.1.2实验装置与材料实验装置主要包括河道模型、供水系统、造冰系统、测量系统等。河道模型采用有机玻璃制作，根据原型河道的地形数据，精确雕刻出河道的弯曲、宽窄变化等特征，确保模型的几何相似性。供水系统由水箱、水泵、管道等组成，能够稳定地提供不同流量的水流，满足实验对水流条件的要求。造冰系统用于制作模型冰，采用特制的制冰设备，模拟自然条件下冰的形成过程。测量系统配备了激光水位计、粒子图像测速仪（PIV）、冰厚测量仪等先进仪器，能够实时准确地测量水位、流速、冰厚等物理量。实验材料方面，模型冰是关键材料。选用与实际冰块物理性质相似的材料制作模型冰，如采用透明的冰模拟材料，其密度、强度、融化特性等与真实冰相近。在模型冰中添加示踪粒子，以便利用PIV技术测量冰块的运动轨迹和速度。在河道模型表面铺设特定的糙率材料，模拟实际河道的粗糙度，确保模型与原型在水流阻力方面的相似性。4.1.3相似准则确定为保证物理模拟实验的有效性，确定了几何、运动、动力相似准则。几何相似准则要求模型与原型的各部分几何尺寸成一定比例，即长度比尺\lambda_{L}、面积比尺\lambda_{A}和体积比尺\lambda_{V}满足：\lambda_{L}=\frac{L_{m}}{L_{p}}，\lambda_{A}=\lambda_{L}^{2}，\lambda_{V}=\lambda_{L}^{3}其中，L_{m}为模型长度，L_{p}为原型长度。根据实验条件和研究精度要求，确定几何比尺\lambda_{L}=1:100。运动相似准则要求模型与原型的水流运动状态相似，即流速比尺\lambda_{v}、加速度比尺\lambda_{a}和时间比尺\lambda_{t}满足：\lambda_{v}=\frac{v_{m}}{v_{p}}，\lambda_{a}=\frac{a_{m}}{a_{p}}，\lambda_{t}=\frac{t_{m}}{t_{p}}=\frac{L_{m}/v_{m}}{L_{p}/v_{p}}=\frac{\lambda_{L}}{\lambda_{v}}通过弗劳德数相等来实现运动相似，即Fr_{m}=Fr_{p}，其中Fr=\frac{v}{\sqrt{gL}}为弗劳德数，g为重力加速度。由此可得流速比尺\lambda_{v}=\sqrt{\lambda_{L}}，在几何比尺\lambda_{L}=1:100的情况下，流速比尺\lambda_{v}=1:10。动力相似准则要求模型与原型的作用力相似，即重力、摩擦力、惯性力等在模型和原型中的比例相同。通过雷诺数相等来考虑摩擦力相似，即Re_{m}=Re_{p}，其中Re=\frac{vL}{\nu}为雷诺数，\nu为运动粘性系数。由于实验中采用的流体与原型相同，运动粘性系数比尺\lambda_{\nu}=1，结合流速比尺和长度比尺，可确定在满足弗劳德数相似的条件下，雷诺数也能基本满足相似要求。通过这些相似准则的确定和严格控制，保证了物理模拟实验能够准确反映黑龙江干流上游凌汛演进的实际过程。4.2实验过程与数据采集4.2.1实验操作步骤在实验前，首先对实验装置进行全面检查和调试。确保河道模型的表面光滑，无裂缝和凸起，以保证水流和冰凌运动的顺畅。检查供水系统的密封性和稳定性，确保水泵能够正常工作，提供稳定的流量。对造冰系统进行预冷，使其达到制冰所需的温度条件。对测量系统的仪器进行校准，确保测量数据的准确性。实验开始时，启动供水系统，按照预定的流量向河道模型中注水。在注水过程中，通过调节水泵的功率和阀门的开度，精确控制流量的大小。待水位达到一定高度后，开启造冰系统，在河道模型表面制作冰盖。在制作冰盖时，模拟自然条件下冰的生长过程，控制制冰的速度和温度，使冰盖的厚度和强度符合实际情况。在冰盖形成后，利用冰厚测量仪测量冰盖的厚度，确保冰盖厚度均匀。当冰盖达到预定厚度后，调整流量和水温，模拟凌汛期间的水流和热力条件。逐渐增加流量，模拟上游来水的增加，观察冰凌的输移和堆积情况。同时，通过调节水温，模拟气温的变化，观察冰盖的融化和冰凌的生消过程。在实验过程中，保持实验环境的温度稳定，避免外界温度变化对实验结果的影响。在冰凌堆积过程中，密切关注冰塞、冰坝的形成和发展。当发现冰塞、冰坝开始形成时，记录其形成的位置和时间。使用激光水位计实时测量冰塞、冰坝上游和下游的水位变化，利用PIV技术测量水流速度和冰凌的运动轨迹。在冰坝形成后，观察冰坝的稳定性，记录冰坝溃决的时间和过程。实验结束后，停止供水和造冰系统，清理河道模型中的冰块和水。对实验数据进行整理和记录，将测量得到的水位、流速、冰厚等数据进行分类存储，为后续的数据分析和处理做好准备。对实验装置进行维护和保养，检查仪器设备的性能，确保下次实验的正常进行。4.2.2数据采集方法与频率水位数据采用激光水位计进行测量。激光水位计通过发射激光束，测量激光束从发射到接收的时间，根据光速和测量时间计算水位高度。在河道模型中，沿水流方向设置多个测量断面，每个断面布置3-5个测量点，以获取不同位置的水位信息。在实验过程中，每30秒采集一次水位数据，确保能够准确捕捉水位的变化过程。流速数据利用粒子图像测速仪（PIV）进行测量。PIV技术通过在流场中添加示踪粒子，利用激光片光源照射流场，拍摄示踪粒子的图像，通过图像分析计算示踪粒子的位移和速度，从而得到水流的流速分布。在实验中，选择具有代表性的区域进行PIV测量，如河道弯曲处、冰坝形成区域等。每隔5分钟进行一次PIV测量，每次测量持续3-5分钟，以获取不同时刻的流速信息。冰情数据的采集包括冰厚、冰凌密度和冰盖面积等。冰厚使用冰厚测量仪进行测量，冰厚测量仪采用超声波原理，通过发射和接收超声波信号，测量冰盖的厚度。在冰盖表面均匀布置多个测量点，每隔10分钟测量一次冰厚。冰凌密度通过在河道模型中设置采样区域，定期采集冰块样本，测量样本中冰块的体积和质量，计算冰凌密度。冰盖面积则通过图像采集系统拍摄河道模型的照片，利用图像处理软件分析照片，计算冰盖的覆盖面积。在实验过程中，除了上述主要数据的采集外，还同步记录实验过程中的气温、水温等环境参数。气温使用高精度温度计测量，每隔15分钟记录一次。水温通过在河道模型中布置温度传感器进行测量，实时记录水温的变化。通过全面、准确的数据采集，为深入分析黑龙江干流上游凌汛演进过程提供了丰富的数据支持。4.3实验结果分析4.3.1凌汛现象观察与记录在实验过程中，清晰观察到了冰坝的形成过程。随着实验的推进，当水流携带冰凌进入河道弯曲处时，冰凌开始逐渐堆积。由于水流速度在弯曲处减缓，冰凌的输移能力下降，冰块相互挤压、堆叠。起初，冰凌在河道一侧堆积，形成小型的冰堆。随着上游来冰量的增加，冰堆逐渐扩大，并向河道另一侧延伸，最终形成横跨河道的冰坝。冰坝形成后，上游水位迅速上涨，在冰坝上游约1-2米范围内，水位涨幅可达0.2-0.3米。冰坝的厚度在0.1-0.2米之间，其表面较为粗糙，冰块之间相互交错，形成了紧密的结构。冰坝溃决也是实验中重点观察的现象。当冰坝上游水位持续升高，超过冰坝的承受能力时，冰坝开始出现裂缝。裂缝首先出现在冰坝的薄弱部位，如冰块之间的连接处。随着裂缝的逐渐扩大，冰坝的稳定性受到严重影响。最终，冰坝在水流的冲击下突然溃决，大量冰块和水流瞬间向下游涌去。冰坝溃决时，产生了巨大的冲击力，下游水流速度急剧增大，在溃决点下游约5-10米范围内，水流速度可达原来的2-3倍。冰坝溃决后，上游水位迅速下降，水位下降速度在1-2厘米/秒左右。实验中还观察到冰塞的形成过程。在封冻期或开河初期，当水流速度较小时，冰凌在河道底部和岸边逐渐聚集。这些冰凌相互交错，形成了类似于堵塞的冰塞。冰塞的长度在0.5-1米之间，厚度在0.05-0.1米左右。冰塞的存在导致河道过水能力下降，上游水位逐渐抬高。在冰塞发展过程中，水位升高速度较为缓慢，约为0.5-1厘米/小时。随着冰塞的不断发展，当冰凌堆积到一定程度时，可能会引发冰坝的形成。4.3.2实验数据整理与分析对采集到的水位数据进行整理和分析，得到了水位随时间的变化曲线。在冰坝形成前，水位变化较为平稳，呈现出缓慢上升的趋势。当冰坝开始形成时，水位迅速上升，上升速度明显加快。在冰坝形成后的一段时间内，水位保持在较高水平。随着冰坝的逐渐消融或溃决，水位又逐渐下降。通过对不同工况下水位数据的对比分析，发现流量越大，冰坝形成后水位上升的幅度越大。在高流量工况下，冰坝形成后水位涨幅可达0.4-0.5米，而在低流量工况下，水位涨幅仅为0.1-0.2米。流速数据的分析结果表明，在河道的不同位置，流速存在明显差异。在河道主流区域，流速较大，而在河道弯曲处和冰坝形成区域，流速较小。当冰坝形成后，冰坝上游的流速明显减小，而冰坝下游的流速在冰坝溃决瞬间会急剧增大。通过对流速分布的分析，还发现流速的变化与冰情密切相关。当冰凌堆积形成冰塞或冰坝时，会改变河道的过水断面和水流阻力，从而影响流速的分布。冰情数据的分析显示，冰厚在实验过程中呈现出先增加后减小的趋势。在封冻期，冰厚逐渐增加，达到最大值后，随着气温升高和冰盖的融化，冰厚逐渐减小。冰凌密度在冰坝形成区域明显增大，在冰坝形成前，冰凌密度约为0.2-0.3千克/立方米，而在冰坝形成后，冰凌密度可达到0.5-0.8千克/立方米。冰盖面积在实验过程中也会发生变化，随着冰盖的融化和冰凌的输移，冰盖面积逐渐减小。通过对这些实验数据的综合分析，得出了凌汛演进的一些规律。流量和冰量是影响冰坝形成和凌汛规模的重要因素，流量越大、冰量越多，冰坝越容易形成，且规模越大，凌汛灾害的风险也越高。气温的变化对冰盖的生消和冰凌的运动起着关键作用，气温升高会加速冰盖的融化，增加冰凌的数量和移动速度。河道地形对冰凌的堆积和冰坝的形成也有重要影响，河道弯曲处和狭窄处容易导致冰凌堆积，增加冰坝形成的概率。4.3.3与实际情况对比将物理模拟实验结果与历史凌汛实际情况进行对比，发现两者在冰坝形成位置和水位变化趋势上具有一定的相似性。在历史凌汛中，冰坝多形成于河道弯曲、宽窄变化较大的河段，这与实验中冰坝形成的位置相符。实验中冰坝形成后水位迅速上涨，与历史凌汛中冰坝导致水位急剧升高的现象一致。在1985年黑龙江干流上游的特大凌汛中，冰坝形成后上游水位涨幅超过5米，实验中在类似的水流和冰情条件下，冰坝形成后水位涨幅也较为明显。实验结果与实际情况也存在一些差异。在实际凌汛中，由于河道条件更为复杂，受到多种因素的综合影响，冰坝的规模和稳定性可能与实验结果有所不同。实际河道中的水流速度、冰凌的粒径分布等因素可能更加复杂多变，而实验中难以完全模拟这些复杂情况。实际凌汛的影响范围更广，涉及到更长的河段和更多的区域，实验只能模拟部分典型河段的凌汛过程。针对这些差异，进一步分析原因，发现主要是由于实验条件的局限性和实际情况的复杂性所致。为了提高模拟的准确性，需要进一步优化实验方案，考虑更多的实际因素，如河道的糙率变化、冰凌的破碎和融化过程等。五、数值模拟与物理模拟对比验证5.1对比方法与指标为了准确评估数值模拟和物理模拟在黑龙江干流上游凌汛演进过程研究中的效果，采用多种对比方法和指标。在对比方法上，将数值模拟和物理模拟在相同工况下的结果进行直接对比。对于某一特定的流量、冰量和气温条件，分别获取数值模拟和物理模拟得到的水位、流速、冰情等数据，直观地比较两者的差异。在对比指标方面，选取水位、流速、冰情等关键参数作为评估依据。水位是衡量凌汛灾害程度的重要指标，对比数值模拟和物理模拟得到的水位变化过程，包括水位的上升速度、最高水位以及水位变化的时间节点。通过计算两者水位数据的均方根误差（RMSE）和平均绝对误差（MAE）来量化差异。均方根误差能够反映模拟值与观测值之间的平均误差程度，其计算公式为：RMSE=\sqrt{\frac{1}{n}\sum_{i=1}^{n}(y_{i}-\hat{y}_{i})^2}其中，n为数据点的数量，y_{i}为物理模拟得到的水位观测值，\hat{y}_{i}为数值模拟得到的水位模拟值。平均绝对误差则衡量模拟值与观测值之间绝对误差的平均值，计算公式为：MAE=\frac{1}{n}\sum_{i=1}^{n}|y_{i}-\hat{y}_{i}|流速对冰凌的输移和冰坝的形成有着重要影响，对比两者在不同位置的流速分布情况。通过绘制流速分布图，直观地展示数值模拟和物理模拟结果的差异。利用相关系数（R）来评估两者流速数据的相关性，相关系数取值范围在-1到1之间，越接近1表示相关性越强。冰情参数如冰厚、冰凌密度等也是对比的重点。对比数值模拟和物理模拟得到的冰厚变化过程以及冰凌密度在河道中的分布情况。同样采用均方根误差和平均绝对误差来评估冰厚数据的差异，对于冰凌密度，通过计算两者在不同区域的相对误差来衡量差异程度。相对误差的计算公式为：RE=\frac{|y_{i}-\hat{y}_{i}|}{y_{i}}\times100\%其中，RE为相对误差，y_{i}为物理模拟得到的冰凌密度观测值，\hat{y}_{i}为数值模拟得到的冰凌密度模拟值。通过这些对比方法和指标，全面、客观地评估数值模拟和物理模拟的准确性和可靠性，为进一步改进模拟方法和提高模拟精度提供依据。5.2结果对比分析5.2.1水位、流量对比通过对数值模拟和物理模拟结果的分析，发现两者在水位和流量的变化趋势上具有一定的一致性，但在具体数值上存在一些差异。在水位变化方面，图4展示了典型断面在凌汛期间的水位对比情况。从图中可以看出，数值模拟和物理模拟得到的水位变化曲线总体趋势相似，都呈现出在凌汛初期水位缓慢上升，冰坝形成后迅速上升，冰坝溃决或消融后逐渐下降的过程。在冰坝形成后的最高水位上，数值模拟结果略高于物理模拟结果。数值模拟得到的最高水位为4.8米，而物理模拟的最高水位为4.5米。这可能是由于数值模拟中对河道糙率、冰凌堆积等因素的处理与实际情况存在一定偏差，导致水流阻力的模拟不够准确，从而使得水位计算结果偏高。在水位上升和下降的速率上，两者也存在细微差异，物理模拟的水位上升速率相对较慢，这可能与实验中冰坝形成和溃决的过程相对较为缓慢有关。在流量变化方面，图5给出了凌汛期间流量的对比结果。数值模拟和物理模拟的流量变化趋势基本一致，在凌汛初期，随着冰盖融化和上游来水增加，流量逐渐增大。当冰坝形成时，流量出现短暂下降，冰坝溃决后，流量急剧增大形成洪峰。在洪峰流量的数值上，数值模拟结果比物理模拟结果高出约10-15%。数值模拟得到的洪峰流量为1200立方米每秒，而物理模拟的洪峰流量为1050立方米每秒。这可能是因为数值模拟在计算冰凌对水流的阻挡作用时，对冰凌堆积的密实程度和过水能力的估计与实际情况有所不同，导致流量计算结果偏大。数值模拟在处理复杂水流条件下的流量变化时，可能存在一定的模型误差，进一步影响了洪峰流量的模拟精度。5.2.2冰情对比在冰情模拟方面，数值模拟和物理模拟在冰坝、冰塞的发展过程和特征上也存在一定的差异。图6展示了冰坝形成过程中数值模拟和物理模拟的对比情况。从模拟结果可以看出，两者在冰坝形成的位置上基本一致，都出现在河道弯曲系数大、河宽变化明显的河段。在冰坝的规模和形态上存在差异。物理模拟中冰坝的厚度相对较薄，约为0.1-0.2米，而数值模拟得到的冰坝厚度在0.2-0.3米之间。物理模拟中冰坝的长度较短，约为10-15米，数值模拟的冰坝长度可达15-20米。这可能是由于物理模拟在实验过程中，受到模型尺寸和实验条件的限制，难以完全模拟实际河道中冰凌的大规模堆积和相互作用。数值模拟在计算冰凌堆积时，对冰凌之间的摩擦力、碰撞力等因素的考虑不够准确，导致冰坝的规模模拟偏大。在冰塞的发展过程中，数值模拟和物理模拟也表现出不同的特征。图7给出了冰塞厚度随时间的变化对比。物理模拟中冰塞厚度的增长较为缓慢，在实验开始后的2-3小时内，冰塞厚度从0逐渐增加到0.05-0.1米。而数值模拟中冰塞厚度的增长速度相对较快，在相同时间内，冰塞厚度可达到0.1-0.15米。这可能是因为数值模拟在计算冰塞形成时，对冰凌的输移和堆积过程进行了一定的简化，忽略了一些实际因素的影响，如冰凌在堆积过程中的融化和破碎等，导致冰塞厚度的模拟结果偏大。物理模拟中冰塞的分布相对较为均匀，而数值模拟中冰塞在某些区域的分布较为集中，这可能与两者对河道水流条件和冰凌运动的模拟差异有关。5.3模拟结果验证与评估通过对比分析，对数值模拟和物理模拟结果进行验证与评估。在水位模拟方面，虽然数值模拟和物理模拟结果存在一定差异，但两者的变化趋势基本一致，都能较好地反映凌汛期间水位的变化过程。在冰坝形成前，水位的模拟误差相对较小，均方根误差在0.1-0.2米之间，平均绝对误差在0.05-0.1米之间。这表明在凌汛初期，两种模拟方法都能较为准确地预测水位变化。当冰坝形成后，由于数值模拟对冰凌堆积和水流阻力的处理与实际情况存在偏差，导致水位模拟误差增大，均方根误差达到0.3-0.4米，平均绝对误差在0.2-0.3米之间。但总体来说，数值模拟和物理模拟的水位变化趋势与历史观测数据相符，验证了模拟结果在一定程度上的可靠性。在流速模拟方面，数值模拟和物理模拟在河道主流区域的流速模拟较为接近，相关系数达到0.8-0.9之间。这说明两种模拟方法在描述主流区域的水流速度方面具有较高的准确性。在河道弯曲处和冰坝形成区域，由于水流条件复杂，数值模拟和物理模拟的流速差异较大。在河道弯曲处，数值模拟的流速分布相对较为平滑，而物理模拟能够更准确地反映流速的局部变化。这是因为物理模拟能够直观地展示水流在弯曲河道中的复杂流动形态，而数值模拟在处理复杂地形下的水流时存在一定的局限性。在冰坝形成区域，数值模拟对流速的减小趋势模拟较好，但在流速的具体数值上与物理模拟存在一定偏差。通过与实际观测数据对比，发现物理模拟的流速结果更符合实际情况，这表明物理模拟在复杂水流区域的流速模拟具有优势。在冰情模拟方面，数值模拟和物理模拟在冰坝、冰塞的形成位置和发展趋势上基本一致。这说明两种模拟方法都能较好地预测冰情的变化。在冰坝的规模和冰塞的厚度模拟上存在差异。如前所述，数值模拟的冰坝厚度和长度偏大，冰塞厚度增长速度偏快。通过对模拟结果的验证，发现物理模拟的冰情结果与实际情况更为接近。在实际凌汛中，冰坝的厚度和长度受到多种因素的影响，物理模拟能够更真实地模拟这些因素的作用，而数值模拟在模型简化和参数设置上可能存在不足，导致冰情模拟结果存在偏差。综合来看，数值模拟和物理模拟在黑龙江干流上游凌汛演进过程模拟中都具有一定的准确性和可靠性，但也都存在各自的局限性。在实际应用中，可以结合两种模拟方法的优势，相互验证和补充，提高凌汛模拟的精度和可靠性。5.4模拟结果差异原因探讨数值模拟和物理模拟结果存在差异，主要源于以下几方面原因。数值模拟依赖于模型的建立，模型中存在诸多假设和简化。在水动力模型中，对河道地形的简化可能无法完全反映实际河道的复杂性。在构建模型时，可能将一些微小的河道起伏、支流汇入等细节忽略，这会影响水流的模拟精度。冰情模拟中，对冰盖生消、冰凌堆积等过程的描述采用了简化的数学模型，难以精确刻画冰凌之间复杂的相互作用和物理过程。对冰盖融化的模拟，可能只考虑了气温和水温的影响，而忽略了冰盖表面的粗糙度、太阳辐射等因素对融化速率的影响。这些模型假设和简化导致数值模拟与实际情况存在偏差，进而影响模拟结果。物理模拟存在实验条件的局限性。虽然物理模拟遵循相似准则，但在实际操作中，难以完全满足所有相似条件。在模型冰的制作上，尽管选用了与实际冰