复式断面冰盖下河渠水力特性的多维度试验解析与模型构建

复式断面冰盖下河渠水力特性的多维度试验解析与模型构建一、绪论1.1研究背景在高纬度地区以及部分中纬度地区，冬季的严寒气候使得河渠水面温度急剧下降，当达到冰点以下时，水体开始结冰并逐渐形成冰盖。这一自然现象广泛存在于诸如美国的密西西比河、加拿大的圣劳伦斯河、俄罗斯的勒拿河和叶尼塞河，以及中国东北、新疆和内蒙等地的河渠。冰盖的形成如同给河渠戴上了一顶“冰帽”，彻底改变了河渠原有的水力条件。从水力学的基本原理来看，冰盖的出现改变了河渠的边界条件。在畅流期，河渠的过水断面边界主要为渠底和渠壁，而冰盖形成后，冰盖底面成为了新的边界。这使得过水断面的湿周显著增加，根据谢才公式Q=AC\sqrt{Ri}（其中Q为流量，A为过水断面面积，C为谢才系数，R为水力半径，i为水力坡度），湿周的增大导致水力半径R减小，在其他条件不变的情况下，流量Q会相应减少。冰盖还影响了河渠水流的流速分布。明流状态下，河渠流速分布一般呈抛物线型，最大流速出现在水面附近。而冰盖形成后，由于冰盖底面的阻滞作用以及与渠底糙率的差异，纵向时均流速的垂向分布发生明显变化，呈抛物线型且关于零切应力所在平面通常是非对称的，最大纵向时均流速朝向较光滑的边界一侧。这一变化导致水流的能量分布发生改变，进而影响了水流的输沙能力和物质输运率。冰盖下的水流阻力也与明流时不同。冰盖的存在增加了水流的阻力，使得水流的能量损失增大。这不仅影响了河渠的输水能力，还可能导致水位壅高，在封河初期及开河期，若水位控制不当，极易引发冰凌灾害，如黄河宁蒙河段在1951年-2010年间就发生了将近30年的凌汛灾害，给沿岸居民的生活和经济带来了巨大损失。在实际工程应用中，无论是水电站的运行调度、内陆航运的规划，还是长距离输水工程的设计与运行，都需要准确掌握冰盖下河渠的水力特性。例如，南水北调中线工程等人工河渠在冬季面临着输水线路结冰的问题，冰盖下输水虽然是一种有效的调度运行方式，但如何在保障安全的基础上提升输水效益，就需要深入研究冰盖下河渠的水力特性。1.2研究目的和意义1.2.1目的本研究旨在通过实验，深入探究复式断面冰盖下河渠的水力特性，包括流速分布、水流阻力、水位变化等关键参数的变化规律。通过精确测量和分析这些参数，建立起能够准确描述复式断面冰盖下河渠水力特性的数学模型，为冰盖下河渠的水力计算提供可靠的理论依据。同时，研究不同因素，如冰盖厚度、糙率、河渠断面形状及尺寸等，对水力特性的影响程度，明确各因素的作用机制。最终，将研究成果应用于实际工程案例，验证其在解决实际问题中的有效性和实用性，为水利工程、输水工程等相关领域的设计、运行和管理提供科学指导。1.2.2理论意义在理论层面，冰盖下河渠水力特性研究是河渠冰水力学的重要组成部分，然而当前这一领域的研究尚存在诸多空白和不完善之处。深入研究复式断面冰盖下河渠水力特性，有助于补充和完善河渠冰水力学的理论体系。通过实验获取的一手数据，能够为理论模型的建立和验证提供坚实的数据基础，推动理论研究从定性分析向定量分析转变。在研究流速分布规律时，对不同冰盖条件和河渠断面参数下的流速进行精确测量，能够为流速分布理论的发展提供关键数据支持，使得理论模型更加准确地反映实际流动情况。研究成果还能为其他相关领域的研究提供基础，如冰盖演变、冰凌灾害形成机制等研究，都依赖于对冰盖下河渠水力特性的深入理解。1.2.3实践意义从实践角度来看，冰凌灾害给人类生产生活带来了巨大威胁。在封河初期及开河期，冰盖下河渠的水位变化、流速改变等水力特性的异常，极易引发冰凌堆积，形成冰塞、冰坝，导致水位壅高，河水漫溢，冲毁堤坝、淹没农田和城镇，给人民生命财产安全造成严重损失。通过研究复式断面冰盖下河渠水力特性，能够准确预测冰凌灾害的发生，提前制定有效的防治措施，如合理调控水位、流速，采用破冰、排冰等工程手段，从而降低冰凌灾害带来的损失。对于输水工程而言，在冬季寒冷地区，输水渠道不可避免地会出现结冰现象，冰盖的存在会改变渠道的输水能力和水流状态。若对冰盖下河渠水力特性缺乏了解，可能导致输水不畅、渠道冻胀破坏等问题，影响输水工程的安全运行和供水稳定性。掌握冰盖下河渠水力特性后，可以优化输水工程的运行调度方案，在保证安全的前提下，提高输水效率，降低运行成本。根据冰盖厚度和水力特性，合理调整输水流量，避免冰盖对输水造成阻碍，确保输水工程稳定运行，满足社会经济发展对水资源的需求。1.3国内外研究现状1.3.1国外研究进展国外对于冰盖下河渠水力特性的研究起步较早。在20世纪初，随着寒区水利工程的发展，学者们开始关注冰盖对河渠水流的影响。早期的研究主要集中在冰盖的形成机制和冰情观测方面。在冰盖下水流特性研究方面，H.A.Einstein提出了按阻力划分过流断面的方法，沿水深方向以最大流速线（点）将流动分成两层等效明流层，假定上下两层互不影响，分别只受冰盖底部和床面粗糙影响，每层流速为对数分布。这一方法为后续研究奠定了基础，但也存在一定局限性，一些室内试验及现场观测结果表明，在流动核心部分的流速分布相当均匀，最大流速点位置很难准确确定，流速分布在最大流速点处出现间断，分区流速并不完全相等。随着技术的发展，数值模拟逐渐成为研究冰盖下河渠水力特性的重要手段。imiaghalam等研究了不同冰盖覆盖度和边界粗糙度对梯形河道流速分布的影响，通过数值模拟揭示了冰盖覆盖度和边界粗糙度与流速分布之间的定量关系。运用Fluent软件建立雷诺应力模型，模拟不同冰盖覆盖度下恒定均匀流流场，从纵向时均流速、二次流、雷诺应力、紊动能等方面揭示部分冰盖下水流的水力特性，发现岸冰的形成会导致冰盖下水流在断面上重新分配，明流区流速增大，冰盖下流速减小。在冰盖下水流阻力研究方面，国外学者提出了多种综合糙率计算公式。自1931年Pavlovskiy提出第一个综合糙率的计算公式以来，Lotter、Einstein、Sabaneev、Larsen等学者相继提出了综合糙率的计算公式。Uzuner对上述公式进行了评述，认为Larsen公式理论依据充分，但这些公式在实际应用中仍存在一定误差，需要进一步改进和验证。1.3.2国内研究进展国内对冰盖下河渠水力特性的研究始于20世纪50年代，主要是为了解决黄河凌汛灾害防治等实际工程问题。经过多年发展，在理论研究和工程应用方面都取得了显著成果。在复式断面河渠水力特性研究方面，国内学者通过水槽试验和数值模拟等方法，对复式断面明渠的流速分布、流量分配、阻力特性等进行了深入研究。赵根生等对比分析了国内外研究复式明渠均匀流动的方法，选择传统的断面斜分割方法、Ackers的河槽协同度方法（COH）、Shiono-Knight提出的方法，以及由Cao-ettal提出的复式断面综合糙率的新方法，来计算复式明渠均匀流动，并将计算结果与英国FCF大尺度试验资料进行对比分析。针对冰盖对河渠水力特性的影响，国内学者也开展了大量研究。茅泽育等采用k-ε紊流模型建立了冰盖下水流流动垂向二维数值模型，根据量纲分析理论提出了流速分布规律的影响因素，针对各种因素的不同组合进行了数值计算，并对其流动特性如最大流速点位置、冰区及床面区平均流速等进行了分析研究。王菲菲等应用双层假定和冰盖流涡黏度模型，推导出一个物理意义明确且速度梯度连续的冰封河道纵向时均流速垂向分布解析解，并明确常见的冰盖流对数流速公式和双幂律流速公式的应用条件。在工程应用方面，国内的研究成果为南水北调中线工程等长距离输水工程的冰期运行调度提供了重要技术支持。通过对冰盖下输水特性的研究，提出了一系列保障冰期输水安全的措施，如在闸门前和倒虹吸入口前设置拦冰索、在冰情渠段设置排冰闸及破冰捞冰设备、冰期前做好闸门防冻措施等。1.4研究内容和方法1.4.1研究内容本研究聚焦复式断面冰盖下河渠水力特性，涵盖多个关键方面，旨在全面揭示其内在规律和影响机制。流速分布特性研究：通过试验与数值模拟，深入探究复式断面冰盖下河渠纵向时均流速的垂向和横向分布规律。着重分析冰盖厚度、糙率以及河渠断面形状、尺寸等因素对流速分布的影响。在垂向分布上，研究最大流速点位置的变化，以及不同因素如何导致流速分布偏离传统对数分布规律。在横向分布方面，关注流速在主槽与滩地之间的差异，以及冰盖存在时横向流速梯度的变化。水流阻力特性研究：深入剖析复式断面冰盖下河渠的水流阻力特性，推导综合糙率计算公式。基于双层假定和水流连续性、力的平衡原理，系统分析各因素对水流阻力的影响。将冰盖下复式断面渠道与明渠的阻力系数进行对比，明确冰盖对水流阻力的影响程度。同时，通过对不同糙率公式的分析和验证，确定适合复式断面冰盖下河渠的综合糙率计算方法，为准确计算水流阻力提供依据。动量运输系数研究：对复式断面冰盖下河渠的动量运输系数展开研究，深入分析动能损失强度。探讨有无冰盖情况下横向动能校正系数和动能损失率的变化规律，揭示冰盖对动量运输的影响机制。研究冰盖下复式断面河槽动量运输系数与各水力要素之间的关系，为理解水流的能量转化和输运过程提供理论支持。通过实验和数值模拟，获取不同工况下的动量运输系数数据，分析其变化趋势和影响因素。水位变化特性研究：密切关注复式断面冰盖下河渠在不同水力条件下的水位变化规律，研究冰盖厚度、糙率以及流量等因素对水位的影响。建立水位变化的数学模型，通过理论分析和实际数据验证，提高对水位变化的预测精度。利用实验数据和数值模拟结果，分析水位在冰盖形成和消融过程中的动态变化，为冰凌灾害的预测和防治提供科学依据。1.4.2研究方法本研究综合运用试验研究、理论分析和数值模拟三种方法，多维度、深层次地探究复式断面冰盖下河渠的水力特性。试验研究：搭建专门的试验水槽，模拟复式断面冰盖下河渠的水流条件。采用粒子图像测速（PIV）技术，精确测量冰盖下不同位置的流速分布，获取准确的流速数据。在试验过程中，严格控制冰盖厚度、糙率以及河渠断面形状、尺寸等变量，通过改变这些变量，获取不同工况下的试验数据。对试验数据进行详细分析，绘制流速分布图、阻力系数与各因素的关系曲线等，直观展示水力特性的变化规律。理论分析：基于水力学基本原理，如能量守恒定律、动量定理等，对复式断面冰盖下河渠的水力特性进行深入分析。推导流速分布、水流阻力、动量运输系数等水力参数的计算公式，明确各因素之间的数学关系。运用量纲分析理论，确定影响水力特性的关键物理量，简化复杂的物理过程，为理论模型的建立提供基础。对理论模型进行验证和修正，使其更准确地反映实际水力特性。数值模拟：借助专业的计算流体力学（CFD）软件，如Fluent、ANSYS等，建立复式断面冰盖下河渠的数值模型。通过数值模拟，全面分析不同工况下的水流特性，包括流速分布、水流阻力、水位变化等。在数值模拟过程中，采用合适的湍流模型和边界条件，确保模拟结果的准确性。对数值模拟结果进行可视化处理，生成流速矢量图、压力云图等，直观展示水流的运动状态和水力特性的分布情况。将数值模拟结果与试验数据进行对比验证，不断优化数值模型，提高模拟精度。1.5技术路线本研究采用多维度的技术路线，以全面、深入地探究复式断面冰盖下河渠的水力特性，技术路线如图1-1所示：前期准备：广泛收集国内外关于冰盖下河渠水力特性的研究资料，深入分析现有研究的成果与不足，明确研究的切入点和重点。同时，全面收集相关工程的实际运行数据，如南水北调中线工程的冰期输水数据，为后续研究提供实际案例支持。试验研究：搭建高精度的试验水槽，严格模拟复式断面冰盖下河渠的实际水流条件。运用先进的粒子图像测速（PIV）技术，精确测量冰盖下不同位置的流速分布，获取详细、准确的流速数据。在试验过程中，系统地改变冰盖厚度、糙率以及河渠断面形状、尺寸等变量，全面获取不同工况下的试验数据。对试验数据进行细致分析，绘制流速分布图、阻力系数与各因素的关系曲线等，直观呈现水力特性的变化规律。理论分析：以水力学基本原理为基础，如能量守恒定律、动量定理等，深入剖析复式断面冰盖下河渠的水力特性。推导流速分布、水流阻力、动量运输系数等水力参数的计算公式，明确各因素之间的数学关系。运用量纲分析理论，确定影响水力特性的关键物理量，简化复杂的物理过程，为理论模型的建立奠定基础。对理论模型进行验证和修正，使其更准确地反映实际水力特性。数值模拟：借助专业的计算流体力学（CFD）软件，如Fluent、ANSYS等，建立高精度的复式断面冰盖下河渠数值模型。通过数值模拟，全面分析不同工况下的水流特性，包括流速分布、水流阻力、水位变化等。在数值模拟过程中，采用合适的湍流模型和边界条件，确保模拟结果的准确性。对数值模拟结果进行可视化处理，生成流速矢量图、压力云图等，直观展示水流的运动状态和水力特性的分布情况。将数值模拟结果与试验数据进行对比验证，不断优化数值模型，提高模拟精度。结果对比与验证：将试验结果、理论分析结果和数值模拟结果进行全面对比分析，验证各方法的准确性和可靠性。通过对比，进一步完善理论模型和数值模型，提高对复式断面冰盖下河渠水力特性的认识和理解。成果应用与展望：将研究成果应用于实际工程案例，如南水北调中线工程等，验证研究成果在解决实际问题中的有效性和实用性。根据应用结果，提出针对性的建议和措施，为水利工程、输水工程等相关领域的设计、运行和管理提供科学指导。同时，对未来的研究方向进行展望，明确进一步研究的重点和方向，推动该领域的研究不断深入发展。二、复式断面冰盖下河渠试验设计与实施2.1试验装置与材料本试验依托高精度的试验水槽，精心模拟复式断面冰盖下河渠的水流条件，力求获取准确且具有代表性的试验数据。试验水槽采用优质有机玻璃制成，这种材料具有良好的透光性，能够确保后续的PIV测量工作不受阻碍，同时其稳定性强，能有效维持试验过程中的水流条件稳定。水槽总长15m，这一长度设置既能保证水流在充分发展后进行数据测量，又兼顾了实验室空间的合理利用。在宽度方向上，水槽设计为主槽宽0.8m，滩地宽各0.3m，这种复式断面的设计能够较好地模拟实际河渠中主槽与滩地的水流分布情况。水槽的深度为1.0m，足以满足不同水深条件下的试验需求，确保能够全面研究冰盖下河渠在各种水深情况下的水力特性。在测量仪器方面，采用先进的粒子图像测速（PIV）技术，其核心设备包括PIV相机和脉冲激光器。PIV相机选用高分辨率的型号，分辨率达到4096×4096像素，如此高的分辨率能够精确捕捉示踪粒子的微小位移，从而实现对流速的高精度测量。相机的拍摄频率可达15Hz，能够满足不同流速条件下对测量频率的要求，确保获取到的流速数据具有时效性和准确性。脉冲激光器则采用532nm波长的绿光激光器，单脉冲能量为150mJ。这种波长的激光在水中具有良好的穿透性，能够有效照亮示踪粒子，为相机拍摄提供清晰的散射光斑。高能量的单脉冲能够确保在较短的曝光时间内获得足够亮度的粒子图像，进一步提高测量精度。模拟冰盖的材料选用聚苯乙烯泡沫板，这种材料具有密度小、易于加工的特点，方便根据试验需求切割成不同厚度和形状。其密度约为30kg/m³，与实际冰盖的密度较为接近，能够在一定程度上模拟冰盖的浮力和对水流的阻滞作用。在糙率方面，通过在聚苯乙烯泡沫板表面粘贴不同粗糙度的砂纸来调整，砂纸的粗糙度范围为0.1mm-1.0mm，能够模拟不同糙率的冰盖表面，以研究糙率对冰盖下河渠水力特性的影响。2.2试验方案设计2.2.1工况设置为全面、系统地探究复式断面冰盖下河渠的水力特性，本试验精心设置了多种工况，以涵盖实际工程中可能出现的各种情况。流量设置是试验工况的关键要素之一。考虑到河渠在不同季节和运行条件下的流量变化，本试验选取了5个不同的流量水平，分别为0.05m³/s、0.10m³/s、0.15m³/s、0.20m³/s和0.25m³/s。这些流量值不仅能够反映河渠在枯水期、平水期和丰水期的大致流量范围，还能为研究流量对水力特性的影响提供足够的数据点。通过在不同流量下进行试验，可以观察到流速分布、水流阻力、水位变化等水力参数如何随着流量的改变而变化，从而揭示流量与水力特性之间的内在关系。冰盖厚度也是影响冰盖下河渠水力特性的重要因素。实际河渠中的冰盖厚度会受到气温、水流速度、水体深度等多种因素的影响，变化范围较大。为了模拟这一实际情况，本试验设置了4种不同的冰盖厚度，分别为0.05m、0.10m、0.15m和0.20m。通过改变冰盖厚度，可以研究冰盖对水流的阻滞作用如何随着厚度的增加而变化，以及冰盖厚度对流速分布、水流阻力和水位变化的具体影响。较厚的冰盖可能会使过水断面进一步减小，导致流速增加、水流阻力增大，而较薄的冰盖对水流的影响相对较小。糙率同样是不可忽视的因素，它反映了冰盖表面和河渠底部的粗糙程度，对水流阻力有着显著影响。本试验通过在冰盖表面粘贴不同粗糙度的砂纸以及在河渠底部铺设不同材料来调整糙率。具体设置了3种不同的糙率组合，分别为冰盖糙率n1=0.012、河渠底部糙率n2=0.015；冰盖糙率n1=0.015、河渠底部糙率n2=0.017；冰盖糙率n1=0.018、河渠底部糙率n2=0.020。通过这种方式，可以研究不同糙率组合下冰盖下河渠的水力特性，分析糙率对水流阻力、流速分布等参数的影响规律。当冰盖糙率增大时，水流与冰盖表面的摩擦力增加，可能导致水流阻力增大，流速分布发生改变。综上所述，本试验通过合理设置流量、冰盖厚度和糙率等工况，构建了一个全面、系统的试验体系，为深入研究复式断面冰盖下河渠的水力特性提供了丰富的数据基础。通过对不同工况下试验数据的分析，可以揭示各因素对水力特性的影响机制，为实际工程提供科学的理论依据。2.2.2测量参数在本试验中，为全面揭示复式断面冰盖下河渠的水力特性，需要精确测量一系列关键参数，这些参数对于理解水流运动规律、分析水流阻力机制以及预测水位变化等具有重要意义。流速是反映水流运动状态的关键参数之一，其分布情况直接影响着河渠的输水能力和水流稳定性。本试验采用粒子图像测速（PIV）技术测量流速，该技术能够实现对整个流场的非接触式测量，获取流速的二维分布信息。在试验水槽中，将测量区域划分为多个网格，每个网格的尺寸根据试验精度要求和PIV系统的分辨率确定。在每个网格点上，通过PIV系统拍摄示踪粒子的图像，经过图像处理和分析，得到该点的流速大小和方向。通过对不同工况下流速分布的测量和分析，可以研究冰盖厚度、糙率以及流量等因素对流速分布的影响规律。在冰盖厚度增加时，由于过水断面减小，流速可能会在某些区域增大，且流速分布可能会更加不均匀。水位的准确测量对于评估河渠的输水能力、防洪能力以及冰凌灾害的预测具有重要意义。本试验在水槽的不同位置布置了多个水位传感器，以测量不同工况下的水位变化。水位传感器采用高精度的压力式水位计，其测量精度可达±0.1mm。将水位传感器安装在水槽的侧壁上，确保传感器的探头与水面垂直，以准确测量水位。在试验过程中，实时记录水位数据，并根据测量结果绘制水位沿程变化曲线。通过分析水位变化曲线，可以研究流量、冰盖厚度和糙率等因素对水位的影响。当流量增大时，水位通常会上升；冰盖厚度增加或糙率增大时，水流阻力增大，也可能导致水位升高。阻力系数是衡量水流能量损失的重要指标，它反映了河渠边界对水流的阻滞作用。本试验根据测量的流速、水位和流量等数据，采用曼宁公式计算阻力系数。曼宁公式为n=\frac{1}{C}R^{\frac{1}{6}}，其中n为阻力系数，C为谢才系数，R为水力半径。通过测量不同工况下的流速v、流量Q和过水断面面积A，可以计算出谢才系数C=\frac{v}{\sqrt{Ri}}，进而得到阻力系数n。通过对不同工况下阻力系数的计算和分析，可以研究冰盖下河渠的水流阻力特性，以及冰盖厚度、糙率等因素对阻力系数的影响。冰盖糙率增大时，阻力系数通常会增大，表明水流能量损失增加。通过对流速、水位和阻力系数等关键参数的精确测量和深入分析，可以全面、系统地揭示复式断面冰盖下河渠的水力特性，为实际工程提供科学、可靠的理论依据和技术支持。2.3试验测量技术2.3.1PIV技术原理与应用粒子图像测速（PIV）技术是一种先进的非接触式流速测量技术，其原理基于示踪粒子在流场中的运动。在试验中，首先在待测流体中均匀施放示踪粒子，这些粒子具有良好的跟随性，能够准确地跟随流体的运动。本试验选用的示踪粒子为空心玻璃微珠，其密度与水接近，粒径约为50μm，能够很好地满足试验要求。利用脉冲激光器产生的厚度约为1mm的激光片光，照亮流场中的一个平面。本试验使用的532nm波长的绿光激光器，其单脉冲能量为150mJ，能够有效照亮示踪粒子。跟随流体运动的示踪微粒被激光照射后散射光斑，使用高分辨率的PIV相机以极短的时间间隔Δt拍摄散射光斑，形成两张粒子图像。本试验采用的PIV相机分辨率达到4096×4096像素，拍摄频率为15Hz，能够满足试验对测量精度和频率的要求。将两张粒子图像划分为若干矩形区域，即判读窗口。通过计算两图片中对应判读窗口内粒子图像的互相关函数，根据互相关函数的峰值点位置，即可得到示踪粒子在Δt间隔内的位移Δx。进而根据公式v=\frac{\Deltax}{\Deltat}，得到各窗口对应流体微团的流速矢量。通过对整个测量区域内各判读窗口的流速矢量进行计算和分析，就可以得到整个流场的流速分布情况。在本试验中，PIV技术主要应用于测量复式断面冰盖下河渠的流速分布。在不同流量、冰盖厚度和糙率等工况下，利用PIV技术测量冰盖下不同位置的流速，包括纵向时均流速的垂向和横向分布。通过对测量数据的分析，可以研究各因素对流速分布的影响规律。在冰盖厚度增加时，通过PIV测量可以发现流速在某些区域的变化情况，以及流速分布的均匀性变化。PIV技术还可以测量水流的紊动特性，如紊动能、雷诺应力等，为深入研究冰盖下河渠的水力特性提供更全面的数据支持。2.3.2其他测量方法除了PIV技术用于测量流速外，本试验还采用了多种其他测量方法，以全面获取复式断面冰盖下河渠的水力特性参数。水位测量对于研究冰盖下河渠的水力特性至关重要，它直接反映了水流的能量状态和输水能力。本试验在水槽的不同位置布置了多个压力式水位计，以测量不同工况下的水位变化。压力式水位计的工作原理是基于液体静压力与水深的关系，通过测量液体的压力来计算水位高度。本试验采用的压力式水位计精度可达±0.1mm，能够满足试验对水位测量精度的要求。将水位计的探头安装在水槽侧壁的特定位置，确保探头与水面垂直，以准确测量水位。在试验过程中，实时记录水位数据，并根据测量结果绘制水位沿程变化曲线。通过分析水位变化曲线，可以研究流量、冰盖厚度和糙率等因素对水位的影响。当流量增大时，水位通常会上升；冰盖厚度增加或糙率增大时，水流阻力增大，也可能导致水位升高。水流阻力是影响河渠输水能力和水流能量损失的重要因素，而阻力系数是衡量水流阻力大小的关键参数。本试验根据测量的流速、水位和流量等数据，采用曼宁公式计算阻力系数。曼宁公式为n=\frac{1}{C}R^{\frac{1}{6}}，其中n为阻力系数，C为谢才系数，R为水力半径。通过测量不同工况下的流速v、流量Q和过水断面面积A，可以计算出谢才系数C=\frac{v}{\sqrt{Ri}}，进而得到阻力系数n。在计算过程中，需要准确测量流速和水位，以确保阻力系数的计算精度。通过对不同工况下阻力系数的计算和分析，可以研究冰盖下河渠的水流阻力特性，以及冰盖厚度、糙率等因素对阻力系数的影响。冰盖糙率增大时，阻力系数通常会增大，表明水流能量损失增加。本试验还对冰盖的厚度和糙率进行了测量。冰盖厚度采用高精度的超声波测厚仪进行测量，其测量精度可达±0.1mm。在试验前和试验过程中，定期测量冰盖的厚度，以确保冰盖厚度的稳定性和准确性。糙率的测量则通过在冰盖表面和河渠底部粘贴不同粗糙度的砂纸，并结合经验公式进行估算。通过对冰盖厚度和糙率的准确测量，可以更准确地分析它们对复式断面冰盖下河渠水力特性的影响。2.4试验数据处理方法在本试验中，数据采集依托先进的自动化数据采集系统，该系统与PIV相机、水位传感器等测量仪器无缝连接，能够实现数据的实时、准确采集。在数据采集过程中，严格控制采集频率，对于流速数据，由于其变化较为迅速，采用15Hz的采集频率，确保能够捕捉到流速的瞬时变化。水位数据的变化相对较为缓慢，采集频率设置为1Hz，既能满足对水位变化趋势的监测需求，又不会产生过多冗余数据。对于其他测量参数，如冰盖厚度、糙率等，在每次试验工况改变前和试验过程中定期进行测量和记录。数据整理阶段，首先对采集到的原始数据进行初步筛选和清洗，去除明显异常的数据点。对于PIV测量得到的流速数据，若某点的流速值超出合理范围，如远大于或远小于周围点的流速，且与理论分析结果相差较大，则对该数据点进行检查和修正。可能是由于示踪粒子的局部团聚、激光照射不均匀或图像处理过程中的误差导致。对于水位数据，若出现突变或不合理的波动，通过检查水位传感器的安装状态、测量环境等因素，判断是否存在干扰，并对数据进行相应处理。使用专业的数据处理软件对数据进行深入分析。选用Origin软件进行数据绘图和基本统计分析，它具有强大的绘图功能，能够绘制各种类型的图表，如流速分布图、水位变化曲线、阻力系数与各因素的关系图等。通过绘制纵向时均流速的垂向分布图，可以直观地观察到流速在垂向上的变化规律，分析最大流速点的位置以及不同因素对流速分布的影响。在分析阻力系数与冰盖糙率的关系时，利用Origin软件绘制散点图，并进行曲线拟合，得到阻力系数随冰盖糙率变化的数学表达式。采用MATLAB软件进行复杂的数据处理和模型计算。利用MATLAB的强大矩阵运算和数值计算功能，根据试验数据计算水流阻力、动量运输系数等关键水力参数。基于水力学基本原理，编写相应的MATLAB程序，实现对复式断面冰盖下河渠水力特性的数值模拟和理论分析。通过MATLAB软件，可以对不同工况下的水力参数进行预测和分析，为试验结果的验证和理论模型的建立提供支持。三、复式断面冰盖下河渠流速分布特性3.1冰盖下流速垂向分布规律通过PIV技术，对不同工况下复式断面冰盖下河渠的纵向时均流速垂向分布进行了精确测量与深入分析。图3-1展示了在流量Q=0.15m³/s、冰盖厚度h=0.10m、冰盖糙率n1=0.015、河渠底部糙率n2=0.017工况下，纵向时均流速的垂向分布情况。从图中可以明显看出，纵向时均流速沿垂向呈抛物线型分布，且与传统明流河道的流速分布存在显著差异。在明流河道中，最大流速通常出现在水面附近。而在冰盖下河渠中，由于冰盖底面的阻滞作用，最大流速点下移，且流速分布关于零切应力所在平面通常是非对称的，最大纵向时均流速朝向较光滑的边界一侧。在本工况下，冰盖糙率相对较小，河床糙率相对较大，最大流速点更靠近冰盖底面，这与理论分析和前人研究结果一致。为进一步分析流速垂向分布规律，引入相对流速u^*=\frac{u}{u_{max}}（其中u为某点的流速，u_{max}为最大流速）和相对水深y^*=\frac{y}{h}（其中y为某点距河底的距离，h为水深）。图3-2给出了不同流量工况下相对流速随相对水深的变化曲线。可以看出，随着流量的增加，相对流速在垂向上的分布更加均匀。在低流量工况下，靠近冰盖底面和河底的流速较小，而在高流量工况下，流速分布相对更加均匀，这是因为流量增加时，水流的紊动增强，动量交换更加剧烈，使得流速分布趋于均匀。冰盖和河床糙率对流速垂向分布也有着重要影响。当冰盖糙率增大时，冰盖底面的摩擦力增加，对水流的阻滞作用增强，使得靠近冰盖底面的流速减小，流速分布曲线在冰盖附近变得更加陡峭。相反，当河床糙率增大时，河底对水流的阻力增大，靠近河底的流速减小，流速分布曲线在河底附近变得更加陡峭。通过改变冰盖和河床糙率进行试验，得到不同糙率组合下的流速垂向分布曲线，进一步验证了这一规律。当冰盖糙率从0.012增大到0.018时，靠近冰盖底面0.1h范围内的流速降低了约20%，而河床糙率从0.015增大到0.020时，靠近河底0.1h范围内的流速降低了约15%。图3-1纵向时均流速垂向分布图3-2不同流量下相对流速随相对水深的变化3.2冰盖下流速横向分布规律在复式断面冰盖下河渠中，流速的横向分布同样呈现出独特的规律，受到多种因素的综合影响。以主槽与滩地的交界区域为研究重点，分析流速在横向方向上的变化情况。图3-3展示了在流量Q=0.20m³/s、冰盖厚度h=0.15m、冰盖糙率n1=0.018、河渠底部糙率n2=0.020工况下，横向流速的分布情况。从图中可以清晰地看到，在主槽区域，流速相对较大且分布较为均匀。这是因为主槽过水断面相对较小，水流受到的约束较大，能量较为集中，使得流速较大。而在滩地区域，流速明显较小，且流速分布随横向位置的变化更为显著。这是由于滩地过水断面较大，水流较为分散，能量相对较弱，导致流速较小。在主槽与滩地的交界区域，流速出现明显的梯度变化，这是由于主槽和滩地的水力条件差异较大，水流在交界区域发生动量交换和紊动增强，使得流速分布发生突变。为进一步研究流速横向分布规律，引入横向流速比u_{t}^*=\frac{u_{t}}{u_{m}}（其中u_{t}为某点的横向流速，u_{m}为主槽平均流速）。图3-4给出了不同冰盖厚度工况下横向流速比随横向位置的变化曲线。可以看出，随着冰盖厚度的增加，主槽与滩地的流速差异增大，横向流速梯度也增大。这是因为冰盖厚度增加，过水断面减小，主槽流速进一步增大，而滩地流速受影响相对较小，导致两者流速差异增大。当冰盖厚度从0.05m增加到0.20m时，主槽平均流速增加了约30%，而滩地平均流速仅增加了约10%，主槽与滩地的流速比从2.5增大到3.5。糙率对流速横向分布也有着重要影响。当冰盖糙率增大时，冰盖底面的摩擦力增加，对水流的阻滞作用增强，使得主槽和滩地的流速均减小，但主槽流速减小的幅度相对较小，导致主槽与滩地的流速差异减小。当冰盖糙率从0.012增大到0.018时，主槽平均流速降低了约15%，滩地平均流速降低了约25%，主槽与滩地的流速比从3.0减小到2.5。相反，当河床糙率增大时，河底对水流的阻力增大，滩地流速减小的幅度相对较大，主槽与滩地的流速差异增大。流量的变化同样会影响流速的横向分布。随着流量的增加，主槽和滩地的流速均增大，但主槽流速增大的幅度相对较大，导致主槽与滩地的流速差异增大。当流量从0.05m³/s增加到0.25m³/s时，主槽平均流速增加了约200%，滩地平均流速增加了约150%，主槽与滩地的流速比从2.0增大到2.8。图3-3横向流速分布图3-4不同冰盖厚度下横向流速比随横向位置的变化3.3不同工况下流速分布对比将不同流量、冰盖条件下的流速分布进行对比，能够更清晰地揭示各因素对复式断面冰盖下河渠水力特性的影响规律。图3-5展示了不同流量（Q1=0.05m³/s、Q2=0.15m³/s、Q3=0.25m³/s）和冰盖厚度（h1=0.05m、h2=0.15m）工况下，纵向时均流速的垂向分布情况。从图中可以看出，在相同冰盖厚度下，随着流量的增加，纵向时均流速整体增大，且流速分布更加均匀。当流量从0.05m³/s增加到0.25m³/s时，各水深位置处的流速均有显著提升。在冰盖厚度为0.05m时，水面下0.2m处的流速从0.2m/s增加到0.5m/s。这是因为流量增大时，水流的能量增加，紊动增强，动量交换更加剧烈，使得流速分布更加均匀，流速梯度减小。对比不同冰盖厚度下的流速分布，发现冰盖厚度的增加对流速分布有显著影响。当流量为0.15m³/s时，冰盖厚度从0.05m增加到0.15m，最大流速点下移，且流速分布曲线在冰盖附近变得更加陡峭。这是由于冰盖厚度增加，过水断面减小，冰盖底面与水流的接触面积增大，对水流的阻滞作用增强，靠近冰盖底面的流速减小，流速梯度增大。不同糙率条件下的流速分布也存在明显差异。图3-6展示了冰盖糙率（n11=0.012、n12=0.018）和河渠底部糙率（n21=0.015、n22=0.020）不同组合工况下，纵向时均流速的垂向分布情况。当冰盖糙率增大时，靠近冰盖底面的流速明显减小，流速分布曲线在冰盖附近变得更加陡峭。当冰盖糙率从0.012增大到0.018时，靠近冰盖底面0.1h范围内的流速降低了约20%。而河渠底部糙率增大时，靠近河底的流速减小，流速分布曲线在河底附近变得更加陡峭。通过对不同工况下流速分布的对比分析，可以明确流量、冰盖厚度和糙率等因素对复式断面冰盖下河渠流速分布的影响规律，为进一步研究水力特性提供了重要依据。图3-5不同流量和冰盖厚度下纵向时均流速垂向分布图3-6不同糙率下纵向时均流速垂向分布3.4流速分布模型验证与分析为了验证现有流速分布模型在复式断面冰盖下河渠中的适用性，本研究利用试验所获取的数据，对常见的流速分布模型进行了严格的验证与深入分析。在众多流速分布模型中，普朗特-卡门对数公式在明渠流速分布研究中具有广泛的应用基础。其表达式为u=\frac{u^*}{k}\lny+C，式中u为距河底距离为y处的点流速，u^*为摩阻流速，k为卡门常数，C为积分常数。在明渠流中，该公式能够较好地描述流速的垂向分布规律，然而在冰盖下河渠中，由于冰盖的存在改变了边界条件和水流结构，其适用性需要进一步验证。基于本试验的实测数据，将不同工况下的流速测量值与普朗特-卡门对数公式的计算值进行了详细对比。以流量Q=0.15m³/s、冰盖厚度h=0.10m、冰盖糙率n1=0.015、河渠底部糙率n2=0.017工况为例，图3-7展示了流速测量值与计算值的对比情况。从图中可以明显看出，在靠近冰盖底面和河底的区域，测量值与计算值存在较大偏差。在靠近冰盖底面0.1h范围内，测量值比计算值低约15%-20%，而在靠近河底0.1h范围内，测量值比计算值高约10%-15%。这表明普朗特-卡门对数公式在描述冰盖下河渠流速垂向分布时，在边界附近的准确性较差，无法准确反映冰盖和河底对流速的影响。造成这种偏差的主要原因在于，普朗特-卡门对数公式是基于明渠均匀流的假设推导而来，未充分考虑冰盖的阻滞作用以及冰盖与河底糙率的差异。在冰盖下河渠中，冰盖底面和河底的糙率不同，导致水流在边界处的剪切应力分布不均匀，进而影响了流速分布。冰盖的存在还改变了水流的紊动特性，使得流速分布更加复杂，超出了普朗特-卡门对数公式的适用范围。为了进一步分析流速分布模型的适用性，本研究还引入了相对误差e=\frac{|u_{m}-u_{c}|}{u_{m}}\times100\%（其中u_{m}为测量流速，u_{c}为计算流速）。对不同工况下的相对误差进行统计分析，结果如图3-8所示。可以看出，随着冰盖厚度的增加，相对误差逐渐增大。当冰盖厚度从0.05m增加到0.20m时，平均相对误差从10%左右增加到20%左右。这是因为冰盖厚度增加，冰盖对水流的阻滞作用增强，流速分布的变化更加显著，使得普朗特-卡门对数公式的计算误差增大。流量和糙率的变化也会对相对误差产生影响。随着流量的增加，相对误差呈现先减小后增大的趋势。当流量较小时，水流紊动较弱，普朗特-卡门对数公式的计算误差较大；随着流量增加，水流紊动增强，流速分布更加均匀，计算误差有所减小，但当流量继续增大时，水流的复杂性增加，计算误差又会增大。糙率的变化同样会影响相对误差，当冰盖糙率或河渠底部糙率增大时，相对误差增大，说明糙率对流速分布的影响较大，普朗特-卡门对数公式难以准确描述这种变化。通过本试验数据对普朗特-卡门对数公式等现有流速分布模型的验证与分析可知，这些模型在复式断面冰盖下河渠中的适用性存在一定局限。在实际应用中，需要考虑冰盖下河渠的特殊水力条件，对现有模型进行修正或建立新的模型，以提高对流速分布的预测精度。图3-7流速测量值与计算值对比图3-8不同工况下相对误差统计四、复式断面冰盖下河渠水流阻力特性4.1冰盖下水流阻力形成机制在复式断面冰盖下河渠中，水流阻力的形成是一个复杂的过程，受到多种因素的综合影响。冰盖和河床作为水流的主要边界，它们的存在和特性对水流阻力起着决定性作用。冰盖的存在显著改变了水流的边界条件。在畅流期，河渠的过水断面边界主要为渠底和渠壁，而冰盖形成后，冰盖底面成为了新的边界，使得过水断面的湿周显著增加。根据水力学原理，湿周的增大导致水力半径减小，进而增加了水流的阻力。在相同流量下，冰盖下河渠的水力半径比明流时减小，根据谢才公式Q=AC\sqrt{Ri}，水力半径R减小，谢才系数C也会相应减小，从而使得水流阻力增大。冰盖底面的糙率也会对水流阻力产生重要影响。糙率越大，冰盖底面与水流之间的摩擦力就越大，水流在通过冰盖下时需要克服更大的阻力，导致能量损失增加。河床同样是影响水流阻力的关键因素。河床的糙率反映了河床表面的粗糙程度，糙率越大，河床对水流的阻滞作用就越强。河床的糙率不仅取决于河床的材料和质地，还与河床的形态、粗糙度分布等因素有关。在复式断面河渠中，主槽和滩地的河床糙率可能存在差异，这会导致水流在主槽和滩地的阻力不同，进而影响整个河渠的水流阻力分布。河床的形态也会对水流阻力产生影响。如河床的起伏、弯曲程度等，都会改变水流的流动路径和速度分布，从而影响水流阻力。冰盖和河床对水流结构的改变也是水流阻力形成的重要机制。冰盖的存在使得水流的流速分布发生变化，纵向时均流速的垂向分布呈抛物线型且关于零切应力所在平面通常是非对称的，最大纵向时均流速朝向较光滑的边界一侧。这种流速分布的改变导致水流的能量分布发生变化，使得水流在流动过程中需要克服更多的能量损失，从而增加了水流阻力。冰盖和河床的存在还会影响水流的紊动特性。紊动是水流中能量耗散的重要方式，冰盖和河床的粗糙度会加剧水流的紊动，使得水流的能量损失增加，进而增大了水流阻力。在冰盖下河渠中，由于冰盖和河床的糙率作用，水流的紊动强度比明流时更大，紊动能量耗散也更多。4.2阻力系数影响因素分析4.2.1流量对阻力系数的影响流量作为影响复式断面冰盖下河渠水流阻力的关键因素之一，其变化对阻力系数有着显著影响。通过对不同流量工况下试验数据的深入分析，发现阻力系数与流量之间存在着密切的关联。在本试验中，保持冰盖厚度、糙率以及河渠底部糙率等其他因素不变，仅改变流量大小。图4-1展示了冰盖厚度h=0.10m、冰盖糙率n1=0.015、河渠底部糙率n2=0.017工况下，阻力系数随流量的变化情况。可以看出，随着流量的增加，阻力系数呈现出逐渐减小的趋势。当流量从0.05m³/s增加到0.25m³/s时，阻力系数从0.030减小到0.020，减小了约33%。这一现象可以从水流的紊动特性和能量损失角度进行解释。当流量增加时，水流的能量增大，紊动增强。较强的紊动使得水流内部的动量交换更加剧烈，能够更有效地克服边界的阻力，从而减小了单位水流长度上的能量损失，导致阻力系数减小。流量增加还会使过水断面的流速分布更加均匀，减小了流速梯度，进一步降低了水流的能量损失，使得阻力系数降低。从水力学原理来看，根据谢才公式Q=AC\sqrt{Ri}，流量Q与谢才系数C成正比关系。而阻力系数n与谢才系数C成反比，即n=\frac{1}{C}R^{\frac{1}{6}}。当流量增加时，为了保持公式平衡，谢才系数C增大，从而导致阻力系数n减小。流量对阻力系数的影响还与河渠的断面形状和尺寸有关。在复式断面河渠中，主槽和滩地的水力条件不同，流量变化对主槽和滩地的影响也有所差异。当流量较小时，水流主要集中在主槽，滩地的作用相对较小，此时流量变化对阻力系数的影响主要体现在主槽部分。而当流量增大时，滩地的过水能力逐渐增强，滩地对水流阻力的贡献也相应增加，使得阻力系数的变化趋势变得更加复杂。4.2.2冰盖糙率对阻力系数的影响冰盖糙率是影响复式断面冰盖下河渠水流阻力的重要因素，其大小直接关系到冰盖与水流之间的摩擦力，进而对阻力系数产生显著影响。通过改变冰盖表面的粗糙度，设置不同的冰盖糙率工况，深入研究冰盖糙率与阻力系数之间的关系。在试验中，保持流量、冰盖厚度以及河渠底部糙率等其他因素不变，仅改变冰盖糙率。图4-2展示了流量Q=0.15m³/s、冰盖厚度h=0.10m、河渠底部糙率n2=0.017工况下，阻力系数随冰盖糙率的变化情况。可以清晰地看到，随着冰盖糙率的增大，阻力系数呈现出明显的增大趋势。当冰盖糙率从0.012增大到0.018时，阻力系数从0.022增大到0.028，增大了约27%。这是因为冰盖糙率增大时，冰盖底面的粗糙度增加，水流与冰盖底面之间的摩擦力增大，导致水流在通过冰盖下时需要克服更大的阻力，从而使得能量损失增加，阻力系数增大。冰盖糙率的增大还会影响水流的紊动特性，使得紊动强度增强，进一步加剧了能量损失，导致阻力系数进一步增大。从微观角度来看，冰盖糙率的增大使得冰盖表面的凸起和凹陷增多，水流在流经这些粗糙表面时，会产生更多的漩涡和紊动，这些漩涡和紊动会消耗水流的能量，增加水流的阻力。当冰盖糙率增大时，冰盖表面的摩擦力增大，使得水流的流速分布更加不均匀，靠近冰盖底面的流速减小，流速梯度增大，这也会导致能量损失增加，阻力系数增大。冰盖糙率对阻力系数的影响还与冰盖厚度有关。当冰盖厚度较小时，冰盖糙率对阻力系数的影响相对较大，因为此时冰盖与水流的接触面积相对较大，冰盖糙率的变化对水流阻力的影响更为明显。而当冰盖厚度较大时，冰盖糙率对阻力系数的影响相对较小，因为此时过水断面相对较大，水流的能量相对较强，能够在一定程度上克服冰盖糙率增大带来的阻力增加。4.2.3河床糙率对阻力系数的影响河床糙率作为影响复式断面冰盖下河渠水流阻力的重要因素之一，其对阻力系数的影响不容忽视。河床糙率反映了河床表面的粗糙程度，不同的河床糙率会导致河床与水流之间的摩擦力不同，从而对水流阻力产生显著影响。在本试验中，通过在河渠底部铺设不同材料来调整河床糙率，保持流量、冰盖厚度以及冰盖糙率等其他因素不变，研究河床糙率与阻力系数之间的关系。图4-3展示了流量Q=0.15m³/s、冰盖厚度h=0.10m、冰盖糙率n1=0.015工况下，阻力系数随河床糙率的变化情况。从图中可以明显看出，随着河床糙率的增大，阻力系数呈现出增大的趋势。当河床糙率从0.015增大到0.020时，阻力系数从0.023增大到0.026，增大了约13%。这是因为河床糙率增大时，河床表面的粗糙度增加，水流与河床之间的摩擦力增大，使得水流在流动过程中需要克服更大的阻力，从而导致能量损失增加，阻力系数增大。河床糙率的增大还会影响水流的流速分布，使得靠近河床的流速减小，流速梯度增大，进一步加剧了能量损失，导致阻力系数增大。从水力学原理来看，河床糙率的增大使得河床对水流的阻滞作用增强，水流的能量更多地消耗在克服河床阻力上。根据曼宁公式n=\frac{1}{C}R^{\frac{1}{6}}，河床糙率增大时，谢才系数C减小，从而导致阻力系数n增大。河床糙率对阻力系数的影响还与冰盖糙率和流量有关。当冰盖糙率和流量一定时，河床糙率的变化对阻力系数的影响相对较大。当冰盖糙率较大时，冰盖对水流的阻力已经较大，此时河床糙率的增大对阻力系数的影响相对较小。而当流量较大时，水流的能量较强，能够在一定程度上克服河床糙率增大带来的阻力增加，使得河床糙率对阻力系数的影响相对较小。通过对流量、冰盖糙率和河床糙率等因素对阻力系数影响的分析可知，这些因素之间相互关联、相互影响，共同决定了复式断面冰盖下河渠的水流阻力特性。在实际工程应用中，需要综合考虑这些因素，准确计算水流阻力，为河渠的设计、运行和管理提供科学依据。图4-1阻力系数随流量的变化图4-2阻力系数随冰盖糙率的变化图4-3阻力系数随河床糙率的变化4.3冰盖下复式断面综合糙率计算在复式断面冰盖下河渠中，准确计算综合糙率对于研究水流阻力特性和水力计算具有重要意义。本研究基于双层假定，即沿水深方向以最大流速线（点）将流动分成两层等效明流层，分别考虑冰盖区和床面区的阻力特性，推导综合糙率计算公式。根据水流连续性原理，两层等效明流层的流量之和应等于总流量，即Q=Q_1+Q_2，其中Q为总流量，Q_1为冰盖区流量，Q_2为床面区流量。根据谢才公式，流量可表示为Q=AC\sqrt{Ri}，其中A为过水断面面积，C为谢才系数，R为水力半径，i为水力坡度。对于冰盖区，其流量Q_1=A_1C_1\sqrt{R_1i}，其中A_1为冰盖区过水断面面积，C_1为冰盖区谢才系数，R_1为冰盖区水力半径。对于床面区，其流量Q_2=A_2C_2\sqrt{R_2i}，其中A_2为床面区过水断面面积，C_2为床面区谢才系数，R_2为床面区水力半径。将上述公式代入Q=Q_1+Q_2中，得到AC\sqrt{Ri}=A_1C_1\sqrt{R_1i}+A_2C_2\sqrt{R_2i}。两边同时除以\sqrt{i}，得到AC\sqrt{R}=A_1C_1\sqrt{R_1}+A_2C_2\sqrt{R_2}。根据谢才系数与糙率的关系C=\frac{1}{n}R^{\frac{1}{6}}，将其代入上式中，得到A\frac{1}{n}R^{\frac{2}{3}}=A_1\frac{1}{n_1}R_1^{\frac{2}{3}}+A_2\frac{1}{n_2}R_2^{\frac{2}{3}}，其中n为综合糙率，n_1为冰盖糙率，n_2为河床糙率。通过整理上式，可得到综合糙率的计算公式为n=\frac{A}{A_1\frac{R_1^{\frac{2}{3}}}{n_1}+A_2\frac{R_2^{\frac{2}{3}}}{n_2}}R^{\frac{2}{3}}。利用试验数据对上述综合糙率计算公式进行验证和分析。在流量Q=0.15m³/s、冰盖厚度h=0.10m、冰盖糙率n1=0.015、河渠底部糙率n2=0.017工况下，根据试验测量得到的流速、水位等数据，计算得到冰盖区和床面区的过水断面面积A_1、A_2，水力半径R_1、R_2以及总过水断面面积A和水力半径R。将这些数据代入综合糙率计算公式中，得到计算的综合糙率n。将计算得到的综合糙率与试验中通过曼宁公式计算得到的阻力系数（可转换为糙率）进行对比，结果表明，计算值与试验值较为接近，相对误差在10%以内。这表明本研究推导的综合糙率计算公式能够较好地反映复式断面冰盖下河渠的实际情况，具有较高的准确性和可靠性。进一步分析不同工况下综合糙率的变化规律，发现综合糙率随着冰盖糙率和河床糙率的增大而增大。当冰盖糙率从0.012增大到0.018时，综合糙率增大了约15%；当河床糙率从0.015增大到0.020时，综合糙率增大了约10%。综合糙率还与流量和冰盖厚度有关，随着流量的增加，综合糙率减小；随着冰盖厚度的增加，综合糙率增大。当流量从0.05m³/s增加到0.25m³/s时，综合糙率减小了约20%；当冰盖厚度从0.05m增加到0.20m时，综合糙率增大了约25%。4.4与明渠水流阻力对比分析将冰盖下复式断面渠道的阻力系数与明渠进行对比，能够清晰地揭示冰盖对水流阻力的独特影响，进一步深化对复式断面冰盖下河渠水力特性的认识。在相同流量、断面形状和尺寸条件下，明渠的阻力系数相对较小，而冰盖下复式断面渠道的阻力系数明显增大。在流量Q=0.15m³/s、河渠断面尺寸相同的情况下，明渠的阻力系数约为0.018，而冰盖厚度为0.10m、冰盖糙率n1=0.015、河渠底部糙率n2=0.017的冰盖下复式断面渠道的阻力系数达到0.023，相比明渠增大了约28%。冰盖的存在是导致阻力系数增大的主要原因。冰盖的底面成为新的边界，增加了过水断面的湿周，使得水力半径减小，从而增大了水流阻力。冰盖底面的糙率也会对水流阻力产生影响，糙率越大，冰盖与水流之间的摩擦力越大，水流阻力也就越大。流量对冰盖下和明渠水流阻力系数的影响趋势也有所不同。在明渠中，随着流量的增加，阻力系数通常会减小，这是因为流量增大时，水流的紊动增强，动量交换更加剧烈，能够更有效地克服边界的阻力，从而减小了单位水流长度上的能量损失，导致阻力系数减小。而在冰盖下复式断面渠道中，虽然流量增加时阻力系数也会减小，但减小的幅度相对较小。当流量从0.05m³/s增加到0.25m³/s时，明渠的阻力系数减小了约30%，而冰盖下复式断面渠道的阻力系数仅减小了约20%。这是因为冰盖的存在限制了水流的紊动发展，使得流量增加对阻力系数的影响相对较弱。糙率对冰盖下和明渠水流阻力系数的影响也存在差异。在明渠中，糙率的变化对阻力系数的影响较为显著，糙率增大时，阻力系数明显增大。而在冰盖下复式断面渠道中，冰盖糙率和河床糙率的变化对阻力系数的影响相对复杂。冰盖糙率增大时，虽然会使冰盖与水流之间的摩擦力增大，导致阻力系数增大，但同时也会改变水流的流速分布，使得水流的紊动增强，在一定程度上又会减小阻力系数。河床糙率的变化同样会对阻力系数产生影响，但这种影响还受到冰盖糙率和流量等因素的制约。当冰盖糙率较大时，河床糙率的变化对阻力系数的影响相对较小；而当流量较大时，河床糙率的变化对阻力系数的影响也会减弱。通过对冰盖下和明渠水流阻力的对比分析可知，冰盖的存在显著改变了复式断面渠道的水流阻力特性，在实际工程应用中，需要充分考虑冰盖的影响，准确计算水流阻力，以确保水利工程的安全运行和高效输水。五、复式断面冰盖下河渠动量运输系数研究5.1冰盖下动能损失强度分析在复式断面冰盖下河渠中，动能损失是一个复杂的物理过程，其强度受到多种因素的综合影响。冰盖的存在改变了水流的边界条件和流场结构，导致水流在运动过程中能量不断耗散，动能损失增加。从动能损失的机理来看，冰盖与水流之间的摩擦力是导致动能损失的主要原因之一。冰盖底面的粗糙度使得水流在流经冰盖下时，与冰盖表面产生摩擦，这种摩擦作用消耗了水流的动能，将其转化为热能等其他形式的能量，从而导致动能损失。当冰盖糙率增大时，冰盖底面的粗糙度增加，水流与冰盖之间的摩擦力增大，动能损失也随之增大。水流内部的紊动也是动能损失的重要因素。在冰盖下河渠中，由于冰盖的阻滞作用和河床的影响，水流的紊动强度比明流时更大。紊动使得水流内部的质点相互碰撞、混合，导致能量的耗散，进一步加剧了动能损失。在冰盖附近和河床附近，由于流速梯度较大，紊动更为剧烈，动能损失也更为明显。为了定量分析冰盖下的动能损失强度，引入横向动能校正系数\alpha和动能损失率\eta。横向动能校正系数\alpha反映了流速分布对动能计算的影响，其计算公式为\alpha=\frac{\int_{A}u^{3}dA}{u_{m}^{3}A}，其中u为某点的流速，u_{m}为断面平均流速，A为过水断面面积。动能损失率\eta则表示动能损失的相对大小，其计算公式为\eta=\frac{E_{0}-E}{E_{0}}\times100\%，其中E_{0}为无冰盖时的动能，E为有冰盖时的动能。在流量Q=0.15m³/s、冰盖厚度h=0.10m、冰盖糙率n1=0.015、河渠底部糙率n2=0.017工况下，通过试验测量得到流速分布数据，进而计算得到横向动能校正系数\alpha和动能损失率\eta。计算结果表明，横向动能校正系数\alpha为1.15，动能损失率\eta为12%。这表明在该工况下，流速分布的不均匀性导致动能计算结果比按断面平均流速计算的结果大15%，而冰盖的存在使得动能损失达到了12%。进一步分析不同工况下横向动能校正系数\alpha和动能损失率\eta的变化规律。随着冰盖厚度的增加，横向动能校正系数\alpha增大，动能损失率\eta也增大。当冰盖厚度从0.05m增加到0.20m时，横向动能校正系数\alpha从1.08增大到1.25，动能损失率\eta从8%增大到18%。这是因为冰盖厚度增加，冰盖对水流的阻滞作用增强，流速分布更加不均匀，紊动加剧，导致动能损失增加。流量的变化对横向动能校正系数\alpha和动能损失率\eta也有影响。随着流量的增加，横向动能校正系数\alpha减小，动能损失率\eta减小。当流量从0.05m³/s增加到0.25m³/s时，横向动能校正系数\alpha从1.20减小到1.10，动能损失率\eta从15%减小到10%。这是因为流量增加，水流的能量增大，紊动增强，流速分布更加均匀，动能损失减小。糙率的变化同样会影响横向动能校正系数\alpha和动能损失率\eta。当冰盖糙率或河床糙率增大时，横向动能校正系数\alpha增大，动能损失率\eta增大。当冰盖糙率从0.012增大到0.018时，横向动能校正系数\alpha从1.10增大到1.20，动能损失率\eta从10%增大到14%。这是因为糙率增大，冰盖与水流之间的摩擦力增大，紊动加剧，导致动能损失增加。5.2有无冰盖对动量运输的影响为深入探究冰盖对复式断面河渠动量运输的影响，对有无冰盖情况下的横向动能校正系数和动能损失率进行了对比分析。在无冰盖的明渠情况下，横向动能校正系数相对较小，流速分布相对较为均匀。以相同流量Q=0.15m³/s、河渠底部糙率n2=0.017工况为例，无冰盖时横向动能校正系数约为1.05，这表明流速分布较为均匀，动能计算结果与按断面平均流速计算的结果差异较小。而在有冰盖的情况下，横向动能校正系数明显增大。在冰盖厚度h=0.10m、冰盖糙率n1=0.015的工况下，横向动能校正系数达到1.15，比无冰盖时增大了约9.5%。这说明冰盖的存在使得流速分布更加不均匀，动能计算结果比按断面平均流速计算的结果大得更多，进一步证明了冰盖对水流的阻滞作用导致流速分布发生改变，进而影响了动量运输。从动能损失率来看，无冰盖时动能损失率较低，水流的能量损失相对较小。在上述相同工况下，无冰盖时动能损失率约为5%。而有冰盖时，动能损失率显著增大，达到12%，比无冰盖时增加了140%。这充分表明冰盖的存在极大地增加了水流的动能损失，改变了水流的能量传输过程。冰盖与水流之间的摩擦力以及冰盖对水流紊动的影响，是导致动能损失增大的主要原因。冰盖底面的粗糙度使得水流在流经冰盖下时，与冰盖表面产生摩擦，消耗了水流的动能。冰盖的阻滞作用还加剧了水流的紊动，使得水流内部的质点相互碰撞、混合，进一步导致能量的耗散。通过对比有无冰盖情况下的横向动能校正系数和动能损失率可知，冰盖的存在对复式断面河渠的动量运输产生了显著影响，改变了水流的流速分布和能量传输特性。在实际工程应用中，如水利工程的规划设计和运行管理，必须充分考虑冰盖的影响，准确评估动量运输和动能损失情况，以确保工程的安全运行和高效输水。5.3动量运输系数与水力参数关系动量运输系数与流速、流量等水力参数密切相关，深入探究它们之间的关系对于理解复式断面冰盖下河渠的水力特性具有重要意义。从理论层面来看，动量运输系数反映了水流中动量传输的效率和方式，而流速作为水流运动的基本参数，直接影响着动量的传输。在冰盖下河渠中，流速分布的不均匀性导致动量在不同区域的传输存在差异。在靠近冰盖底面和河底的区域，流速较小，动量传输相对较弱；而在流速较大的区域，动量传输则更为强烈。通过试验数据的分析，进一步揭示了动量运输系数与流速之间的定量关系。在流量Q=0.15m³/s、冰盖厚度h=0.10m、冰盖糙率n1=0.015、河渠底部糙率n2=0.017工况下，绘制动量运输系数与流速的关系曲线，如图5-1所示。可以看出，随着流速的增加，动量运输系数呈现出先增大后减小的趋势。当流速较小时，动量运输系数随流速的增加而迅速增大，这是因为流速的增加使得水流的紊动增强，动量交换更加剧烈，从而提高了动量运输系数。而当流速增大到一定程度后，动量运输系数随流速的增加而逐渐减小，这可能是由于流速过大导致水流的紊动过于强烈，能量耗散增加，从而降低了动量运输的效率。流量作为另一个重要的水力参数，也对动量运输系数产生显著影响。流量的变化直接导致流速的改变，进而影响动量传输。在不同流量工况下进行试验，保持冰盖厚度、糙率等其他因素不变，得到动量运输系数随流量的变化情况，如图5-2所示。可以发现，随着流量的增加，动量运输系数逐渐增大。当流量从0.05m³/s增加到0.25m³/s时，动量运输系数从0.08增大到0.15。这是因为流量增加时，水流的能量增大，流速增大，紊动增强，使得动量传输更加活跃，动量运输系数增大。冰盖厚度和糙率同样会对动量运输系数与水力参数的关系产生影响。当冰盖厚度增加时，过水断面减小，流速增大，动量运输系数增大。当冰盖厚度从0.05m增加到0.20m时，在相同流量下，流速增大，动量运输系数从0.10增大到0.18。冰盖糙率和河床糙率的增大，会导致水流阻力增大，流速减小，动量运输系数减小。当冰盖糙率从0.012增大到0.018时，动量运输系数从0.12减小到0.10。通过对动量运输系数与流速、流量等水力参数关系的研究可知，这些参数之间相互关联、相互影响，共同决定了复式断面冰盖下河渠的动量传输特性。在实际工程应用中，需要充分考虑这些因素，准确评估动量运输情况，为水利工程的设计、运行和管理提供科学依据。图5-1动量运输系数与流速的关系图5-2动量运输系数与流量的关系六、基于试验结果的数值模拟与模型验证6.1数值模拟方法选择在对复式断面冰盖下河渠水力特性进行研究时，数值模拟成为深入探究复杂水流现象的重要手段。本研究选用雷诺应力模型（RSM）进行数值模拟，该模型具有独特的优势和适用性。雷诺应力模型属于二阶矩封闭模型，与传统基于Boussinesq假设的线性模型不同，它直接求解雷诺应力输运方程。在复式断面冰盖下河渠中，水流受到冰盖和河床的双重影响，流场呈现出复杂的各向异性特征。雷诺应力模型能够充分考虑这种各向异性，精确捕捉水流在不同方向上的动量传输和能量交换，从而更准确地模拟流速分布、水流阻力等水力特性。在冰盖与河床糙率不同的情况下，水流在垂直方向和水平方向的紊动特性存在显著差异，雷诺应力模型能够有效反映这些差异，而传统的线性涡黏模型则难以准确模拟。与其他常用的湍流模型相比，如k-ε模型和k-ω模型，雷诺应力模型在处理复杂边界条件和强各向异性流动时具有明显优势。k-ε模型和k-ω模型基于Boussinesq假设，将雷诺应力与平均速度梯度通过涡黏系数联系起来，这种假设在简单流动中能够取得较好的模拟效果，但在处理冰盖下河渠这种复杂流动时，由于无法准确描述各向异性效应，模拟精度会受到限制。而雷诺应力模型直接求解雷诺应力的输运方程，能够更真实地反映水流的紊动特性和动量传输过程，在模拟冰盖下河渠的水力特性时表现更为出色。在一些研究中，对比了不同湍流模型对冰盖下河渠水流的模拟结果，发现雷诺应力模型能够更准确地预测流速分布和水流阻力，与实验数据的吻合度更高。在模拟冰盖下流速垂向分布时，雷诺应力模型能够准确捕捉到最大流速点的下移以及流速分布的非对称性，而k-ε模型和k-ω模型的模拟结果与实验数据存在较大偏差。综上所述，考虑到复式断面冰盖下河渠水流的复杂性和各向异性特征，本研究选择雷诺应力模型进行数值模拟，以获得更准确、可靠的模拟结果，为深入研究冰盖下河渠的水力特性提供有力支持。6.2模型建立与参数设置基于雷诺应力模型（RSM），运用专业计算流体力学软件Fluent进行数值模拟，构建准确反映复式断面冰盖下河渠水力特性的数值模型。控制方程以不可压缩流体的雷诺平均N-S方程为基础，具体如下：连续性方程：连续性方程：\frac{\partialu_{i}}{\partialx_{i}}=0，此方程确保了流体在流动过程中的质量守恒，在复式断面冰盖下河渠中，无论水流在主槽、滩地还是冰盖下的复杂区域流动，质量都不会凭空产生或消失。动量方程：动量方程：\rho\frac{\partialu_{i}}{\partialt}+\rhou_{j}\frac{\partialu_{i}}{\partialx_{j}}=-\frac{\partialp}{\partialx_{i}}+\frac{\partial}{\partialx_{j}}(\mu\frac{\partialu_{i}}{\partialx_{j}}-\rho\overline{u_{i}'u_{j}'})，该方程描述了流体动量的变化与外力和粘性力之间的关系。在冰盖下河渠中，冰盖和河床的存在使得外力和粘性力的分布变得复杂，动量方程能够准确地反映这种复杂性。其中，\rho为流体密度，u_{i}、u_{j}分别为i、j方向的时均流速，t为时间，p为压力，\mu为动力粘性系数，\overline{u_{i}'u_{j}'}为雷诺应力。在雷诺应力模型中，雷诺应力通过求解雷诺应力输运方程得到：\frac{\partial(\rho\overline{u_{i}'u_{j}'})}{\partialt}+\frac{\partial(\rhou_{k}\overline{u_{i}'u_{j}'})}{\partialx_{k}}=P_{ij}+\phi_{ij}+\epsilon_{ij}，其中P_{ij}为雷诺应力生成项，\phi_{ij}为压力应变项，\epsilon_{ij}为耗散项。这些项的准确求解对于模拟水流的紊动特性和动量传输过程至关重要。采用结构化网格对计算区域进行划分，这种网格具有规则的拓扑结构，能够提高计算效率和精度。在复式断面冰盖下河渠的数值模型中，对冰盖下区域和河渠底部区域进行加密处理。在冰盖下区域，由于冰盖与水流的相互作用强烈，流速变化较大，加密网格能够更准确地捕捉流速的变化和紊动特性。对河渠底部区域加密，是因为河底的糙率对水流阻力和流速分布有重要影响，加密网格可以更好地模拟河底边界条件对水流的作用。通过这种网格划分策略，能够确保在关键区域获得更高的计算精度，从而提高整个数值模型的准确性。参数设置方面，根据试验条件和实际物理特性进行合理取值。流体密度\rho取值为1000kg/m³，这是水在常温常压下的典型密度值。动力粘性系数\mu取值为0.001Pa・s，符合水的粘性特性。在边界条件设置上，进口采用速度入口边界条件，根据试验设置的流量和过水断面面积，计算得到进口流速，将其作为速度入口的边界条件。出口采用自由出流边界条件，即出口处的压力为大气压，水流自由流出。冰盖和河床边界采用无滑移边界条件，即冰盖底面和河床表面的流速为零，这反映了实际情况下水流与边界之间的摩擦作用。6.3模拟结果与试验结果对比验证将数值模拟得到的流速分布、阻力系数等结果与试验结果进行对比，是验证数值模型准确性和可靠性的关键环节。通过对比，能够直观地评估数值模型在模拟复式断面冰盖下河渠水力特性方面的能力，为进一步优化模型和深入研究提供重要依据。在流速分布方面，以流量Q=0.15m³/s、冰盖厚度h=0.10m、冰盖糙率n1=0.015、河渠底部糙率n2=0.017工况为例，图