海岸与河口潮流泥沙模拟技术规程71

'关于发布《海岸与河口潮流泥沙模拟技术规程》的通知交水发[1998]761号各省、自治区、直辖市交通厅（局、委、办），部属及双重领导企事业单位：由我部组织交通部天津水运工程科学研究所等单位制定的《海岸与河口潮流泥沙模拟技术规程》，业经审查，现批准为推荐性标准，编号为JTJ/T233—98，自1999年5月1日起施行。本规程的管理和出版组织工作由部水运司负责，具体解释工作由交通部天津水运工程科学研究所负责。中华人民共和国交通部一九九八年十二月十一日制定说明（条文说明） 本规程根据交通部交基发[1994]873号文制定。主编单位为交通部天津水运工程科学研究所，参加单位为南京水利科学研究院、天津港湾工程研究所、河海大学。本规程编写过程中吸取了国内外有关研究经验，广泛地征求了国内有关单位和专家的意见，于1996年10月完成了送审稿。本规程人员分工如下：1.总则蔡嘉熙曹祖德2.术语、符号吴以喜卢汉才3.一般规定吴以喜杨华4.潮流定床模型试验吴以喜吴明阳5.潮流泥沙模型试验罗肇森杨树森6.波浪泥沙模型试验刘子琪刘顺宽7.河口入海船闸模型试验周华兴8.平面二维潮流泥沙数值模拟徐宏明李蓓9.三维潮流泥沙数值模拟张征李孟国10.波浪泥沙数值模拟张东生本规程于1998年3月31日通过部审，1998年12月11日颁布，1999年5月1日实施。 前言 随着社会广义市场经济建设的发展，海岸与河口水运工程建设及其模拟研究工作也得到了蓬勃的民展。统一、规范海岸与河口潮流泥沙模拟研究工作，确保水运工程建设质量已成为当前迫切需要解决的问题。为此，交通部组织了本规程的编制。本规程包括物理模型和数值模拟两部分，主要内容有海岸与河口潮流定床模型试验、潮流泥沙模型试验、波浪泥沙模型试验、河口入海船闸模型试验、平面二维潮流泥沙数值模拟、三维潮流泥沙数值模拟及波浪泥沙数值模拟等。物理模型和数值模拟各有特点，在使用中应该发挥所长，有机结合。本规程共分10章64节304条，6个附录，并附有条文说明。71 本规程由交通部天津水运工程科学研究所负责解释。海岸与河口潮流泥沙模拟技术随着科技民展和新技术的应用，将会得到不断民展，请各单位在使用中，注意总结经验，将意见告诉解释单位。本规程如进行局部修订，其修订的内容将在《水运工程标准与造价管理信息》上刊登。1总则1.0.1为指导与规范海岸与河口潮流泥沙试验方法和模拟技术，提高研究成果的可靠性，更好地为水运工程建设提供科学依据，制定本规程。1.0.2本规程适用于海岸与河口水运工程建设项目的潮流泥沙模拟，也可适用于海岸与河口其它工程的潮流泥沙模拟研究。1.0.3海岸与河口潮流泥沙模拟研究是集理论与工程实践经验为一体的专业技术性很强的工作，应由专业研究单位承担。1.0.4海岸与河口潮流泥沙模拟研究工作除应符合本规程外，尚应符合国家现行有关标准的规定。条文说明1.0.1海岸与河口地区受到海洋潮汐、波浪、河口径流及泥沙运移等多种自然因素的作用和影响，且水流、泥沙的运动及其内在规律复杂，直到目前还不能全面、充分的认识，直接影响到港口、航道、船闸等水运工程的安全使用与维护。而依靠现场的观测则需费大量人力、物力、财力和时间，且工作条件恶劣，难以达到预期目的，物理模型和数值模拟是研究水流、泥沙运动以及水运工程对水流和海床影响的重要手段。通过模拟技术的应用，则能较好地得以解决。为能更好地运用模拟技术，规范和统一试验技术，提高模拟精度，制定本规程是十分必要的。1.0.2海岸与河口其它工程是指海岸与河口的围垦工程、闸坝工程、取排水工程等一切与海岸、河口岸滩演变、地形冲淤变化有关的工程。2.1术语2.1.1海岸以海洋动力作用为主的海岸地带。2.1.2河口海洋动力与河流动力相互作用的河口区。2.1.3模拟技术模仿自然界物质运动和变化规律的技术，包括物理模型和数值模拟两种方法。2.1.4物理模型71 将研究对象按满足一定相似条件或相似准则，缩制而成的模型，又称实体模型。2.1.5正态模型长、宽、高都按同一几何比尺缩小的物理模型。2.1.6变态模型水平比尺大于垂直比尺的物理模型。2.1.7模型变率模型水平比尺与垂直比尺之比值。2.1.8数值模拟针对研究对象和需要研究问题的数学方程式，按给定的定解条件进行数值求解的方法，又称数学模型。2.1.9入海船闸船舶进出港口和入海河口的过船建筑物。2.1.10潮流定床模型水流为清水，模型床面在潮汐水流作用下不发生变形的物理模型。2.1.11潮流泥沙模型潮流中挟带悬沙和底沙颗粒、床面铺有适当厚度的模型沙，在潮流作用下可冲可淤的物理模型。按模型沙运动型态的不同可分以下两种：以底沙运动为主的底沙动床模型；以悬沙运动为主的悬沙动床冲淤模型。2.1.12悬沙淤积定床模型悬沙动床冲淤模型的一种特例，模型床面做成定床，模型水流挟带悬沙，用以研究潮汐水流挟带悬沙对床面淤积的影响。2.1.13沿岸输沙海岸地区在破波带内由破碎波的作用而形成沿海岸方向的泥沙运动。2.1.14一线模型在海岸横剖面上，某一点水深与该点至岸线的距离成线性关系，并采用一维方法计算沿岸输沙作用下岸滩变化的模型。2.1.15n线模型在海岸横剖面上，某一点水深与该点至岸线的距离n次方成正比，并采用二维方法计算沿岸输沙和横向输沙作用下岸滩变化的模型。2.1.16波浪泥沙模型模拟海岸泥沙运动、床面由模型沙组成，按一定相似比尺复演波浪泥沙运动和地形冲淤变化的模型。2.1.17水工模型观察研究水工建筑物有关水力学和工程设计问题的物理模型。2.1.18几何相似模型与原型的线性度量之间保持固定的比例关系。2.1.19重力相似模型与原型内惯性力与重力在相应方向上的分量比值相等。2.1.20阻力相似模型与原型内阻力与重力在相应方向上的分量比值相等。2.1.21动力相似两个几何相似体系的液流中，相应质点所受性质相同的作用力保持着固定的比例关系。即力场的几何相似。71 2.1.22运动相似两个几何相似体系的液流中，相应质点迹线几何相似，而且质点流过相应线段所需的时间保持固定的比例关系。即速度场的几何相似。2.1.23床面冲淤变形相似模型床面在一定相应时段内与原型床面冲淤变化相似。2.1.24模型沙物理模型中，按相似条件模拟原型悬沙、底沙和模型床面可动的固体颗粒材料，它有天然沙和轻质沙两种。2.1.25边界条件边界上水流、波浪和泥沙控制的条件。2.1.26起始条件模型试验或数值模拟开始时所采取的水流、波浪和泥沙运动的起始状态。2.1.27加糙调整物理模型床面糙率以满足模型水流阻力相似条件的措施。2.1.28密排加糙将加糙的物质（沙、砾或卵石）彼此紧密排列，平铺在模型床面上。2.1.29有间距加糙将加糙的砾或卵石按有间距排列，均匀粘牢在模型床面上。2.1.30水中加糙选择模型沙做泥沙试验时，模型沙的糙率不足，在模型水面或水中增加拉线并挂阻力体，以增加水流阻力的办法。2.1.31桩点法制造物理模型时，根据规则的网格节点高程来塑造模型地形的制模方法。2.1.32断面法制造物理模型时，根据断面地形变化，塑造地形的制模方法。2.1.33等高线法制造物理模型时，根据地形图等高线位置和高程来塑造模型地形的制模方法。2.1.34生潮系统物理模型中，用来复演原型潮汐运动的控制系统。2.1.35验证试验和验证计算物理模型和数值模拟中，为检验和校正模型与原型相似程度的试验和计算。2.1.36方案试验和方案计算物理模型和数值模拟通过验证试验和计算后，根据试验要求研究工程区域流场变化、含沙量分布、地形冲淤变化及其影响等的试验和计算，用以优选工程方案。2.1.37输水系统流量系数船闸中，与输水系统进出口、廊道的几何参数、阀门型式和开启度影响有关的综合系数，它反映水流的畅通程度。2.1.38船闸水力特性曲线船闸输水过程中的各水力因素变化与时间关系的曲线。2.1.39气幕法河口海船闸闸门外侧的底部设置管道，利用压缩空气从小孔排出形成一道向上流动的水气混合屏幕，以破坏盐水楔异重流侵入闸室的方法。2.1.40置换法71 利用咸、淡水密度不同和自然分层的原理，船闸分别设置两套咸、淡水输水系统，海水由闸室底板上支孔进出，淡水从闸墙上部孔口出入，以减少淡水损失的方法。2.1.41潜没咸水厢法利用海水密度大，与淡水分层的原理，将盛满海水的咸水厢潜没在闸室淡水中，根据船舶进出闸室的状况升降咸水厢，从而减少闸室咸淡水体的交换，达到防咸的方法。2.1.42垂向坐标变换一种坐标变换技术，通过变换，将直角坐标系（x，y，z）中三维方程转换成新坐标系（x，y，σ）中的方程，其中x，y坐标不变，仅在垂向坐标中，将自由面和不规则海底转换成新坐标系中平面表面和底层，水体深度变成了单位1。这种变换可增加浅水地区在垂向上的分辨率。2.2符号l——平面长度L——波长h——垂直方向高度，水深，小时H——波高T——波浪周期u——平面坐标中x向流速分量v——平面坐标中y向流速分量——流速矢量V——垂线或断面平均流速V0——泥沙起动流速Vf——泥沙扬动流速ub——浅水波底部轨迹速度Vl——沿岸流流速——底部切应力矢量τ——底部切应力τd——泥沙不淤临界切应力，临界淤积切应力τc——泥沙起动临界切应力，临界冲刷切应力C——波速c——谢才系数α——流向角，波向角，波向线与等深线夹角δ——波坦，f——波浪阻力系数，摩擦系数，柯氏参量fb——底摩擦系数Kθ——波浪绕射系数g——重力加速度d——泥沙粒径γ——水的重度γs——泥沙颗粒的重度71 γo——泥沙的干重度ρ——水体密度ρs——泥沙颗粒密度ω——泥沙静水沉降速度ωf——泥沙有效沉降速度ν——流体的运动粘滞系数S——单位水体含沙重量，含盐度S*——水流挟沙力n——糙率QT——沿岸输沙量Q——流量，输沙率W——潮量，重量G——底沙输沙率G*——底沙输沙能力Re——雷诺数，F——力M——力矩，冲刷系数t——时间to——波浪动力持续作用造成床面冲淤变形的时间t1——潮流、水流持续作用的时间t2——由悬沙造成床面冲淤变形的持续时间t3——由底沙造成床面冲淤变形的持续时间各物理量凡带下角标“ｐ”表示原型的值，带下角标“m”表示模型的值λ——比例常数或相似比尺，当有脚标时，则表示为该下角标物理量的原型与模型比值Pi——i量级波要素出现的频率3.1模拟技术及其应用 3.1.1模拟技术应根据海岸与河口不同工程的性质、研究阶段及当地的水文泥沙地形等自然条件确定。3.1.2在海岸与河口水运工程规划、可行性研究阶段宜采用数值模拟，根据需要情况，也可采用物理模型；初步设计及施工图设计阶段宜采用物理模型；对于重要工程宜同时采用两种模拟方法。3.1.3以潮流作用为主的物理模型可分为潮流定床模型和潮流泥沙模型两大类。3.1.4研究海岸与河口各种工程措施对工程区域及其附近水域的潮位、流场、流态的影响，优化工程方案时，宜采用潮流定床模型。3.1.5当研究要求对工程区域及其附近水域的泥沙淤积和冲刷数量、分布规律进行预报时，宜采用潮流泥沙模型。3.1.6潮流泥沙模型动床部分根据工程影响范围的不同，可采用局部动床潮流泥沙模型。3.1.7潮流泥沙模型根据泥沙运动形态和床面冲淤特点，可分为下列三种情况：（1）当工程区域以悬沙淤积为主时，如天然港湾、港池、岸边挖入式港池、航道开挖等的试验研究，宜采用悬沙淤积定床模型；71 （2）当工程区域泥沙以悬沙形式为主而造成冲淤变化时，宜做悬沙冲淤动床模型；（3）当工程区域泥沙以底沙形式为主而造成冲淤变化时，宜采用底沙动床模型。3.1.8研究以波浪作用为主的沙质海岸的岸滩演变及工程泥沙问题，宜采用波浪泥沙模型。3.1.9当采用潮流和波浪共同作用的泥沙模型时，可参照本规程潮流泥沙模型和波浪泥沙模型的相关规定执行。3.1.10入海河口修建闸坝，需开展船闸水力学和防咸措施等的研究时，宜采用船闸模型。3.1.11海岸与河口潮流泥沙数值模拟根据具体情况可采用平面二维和三维不同数值模拟计算，有时也可采用一维数值模拟计算。3.1.12较大范围的潮流、泥沙冲淤计算和为物理模型及局部范围工程的三维数值模拟提供边界条件，宜用平面二维数值模拟；局部范围工程区域内的潮流、泥沙空间场及床面冲淤计算宜用三维数值模拟。3.1.13沙质海岸由波浪作用引起的岸滩演变和工程泥沙冲淤计算宜用波浪泥沙数值模拟。当需要确定因突堤、丁坝引起岸线变化时，可采用一线模型；当需要确定岸滩的变化时，可采用n线模型。 3.2试验研究大纲 3.2.1模型试验前应根据试验任务的要求，编制试验研究大纲。其内容应包括项目概况、试验研究目的和要求、技术路线、主要研究人员、必备资料、试验研究方法、试验设备及量测仪器、试验研究进度计划、预期目标和研究成果等内容。试验研究过程中，应根据具体情况对试验研究大纲予以不断完善。重要工程的试验研究大纲应增加国内、外研究水平的内容。 3.3仪器设备及记录3.3.1模型试验设备及量测仪器，可根据试验内容和要求参照附录A选择使用。3.3.2试验使用的量测仪器、仪表应通过率定，其技术指标应满足试验精度和稳定性的要求。3.3.3试验记录应根据不同试验的要求和内容制定不同的表格，并按归档要求认真填写。记录值发生错误应划掉重写，不应涂改，原始记录不应重抄。记录应及时整理、计算、校核和装订成册。重要工程的试验应有照相或录像记录。条文说明71 3.3.1～3.3.2针对不同的试验研究应有相应的基本设备和仪器，在我国这些仪器和设备批量小，一般为某一单位生产或自制，并随研究的深入，要求它不断地改进，为不限制其发展，提出了试验仪器、仪表应经率定要求，其技术指标应满足试验精度和稳定性的要求，但又能确保试验成果的质量。3.4报告编写 3.4.1试验报告的格式应符合下列规定：（1）封面应有报告的全称，排在上方居中，单位名称和日期排在下方居中；（2）扉页内容应包括单位主管院（所）长、单位技术负责人、项目负责人、主要参加人、报告编写人等；（3）报告目次；（4）摘要；（5）正文；（6）参考文献和资料：如在正文中引用文献资料，应将文献资料名称、作者、单位、日期按序列于正文后。3.4.2报告正文应包括下列基本内容：（1）前言，说明研究工作的目的、依据、起止时间、工作内容、主要工作方法、所采用的基本资料及其它需要说明的问题；（2）概述，研究地区的自然条件，包括水文、气象、泥沙、地质地貌及工程概况等；（3）模型设计或模拟方法，当采用物理模型时，模型设计应包括设计原则、相似条件、比尺、糙率和模型沙的选择、平面布置及试验仪器设备等；当采用数值模拟时，模拟方法应包括依据的基本方程、数值方法、离散格式、边界条件、参数、模拟依据资料等；（4）模型制作，包括图纸资料、边界和安装精度等；（5）验证资料及验证试验，论证模型设计或数值模拟的合理性及相似性；（6）方案试验或计算及其成果分析；（7）结语。3.4.3报告编写应符合下列规定：（1）报告层次分明、文理清晰、语句精炼通顺，用字规范，标点符号正确；（2）使用法定计量单位；（3）术语应符合有关规定，尚无统一规定的应给出定义；（4）插图应有图名且列于图下方居中，插图中应少用文字说明，宜用阿拉伯数字标以图注号、图注的文字应列于图名的正下方；（5）表格应有表名且列于表上方居中，表格编号应列于右起空一格，表格上部和下部应用实线闭合；（6）结语或结论观点应明确，对存在问题应予明确指出，并提出进一步解决问题的方法或建议。 3.5技术资料归档 3.5.1模拟技术资料是研究工作的重要组成部分，是查询工程试验研究状况及总结提高的重要依据，工作完成后应将技术资料及时归档。3.5.271 归档资料应按项目合同或任务书、依据资料、模型设计、验证试验或验证计算、方案试验或方案计算、终结试验或终结计算、报告成果审批表、项目技术评审意见书等层次进行整理。3.5.3归档资料应符合国家科技档案归档的有关规定。 4.1基本资料4.1.1水深地形图资料应满足下列要求：（1）制模用图范围应比研究区域大，以能满足模型试验要求为准；（2）测图的时间宜尽量与水文测验时间同步，如两者间隔时间较长，应分析论证该间隔时段内水深地形的变化及对水流的影响；（3）测图比尺宜采用1/5000～1/25000，其中试验区测图图比视具体情况宜采用1/5000～1/10000，测图应换算成统一的基准面和坐标系；（4）分析研究用图，包括试验研究工程区域在内的历次测图和海图。4.1.2水文泥沙资料按下列要求确定。4.1.2.1全潮同步水文、泥沙测验资料应满足下列要求：（1）模型试验范围内水文测点垂线数量、垂线布置应根据试验区范围大小和验证试验要求确定，测点垂线不宜少于5条；（2）在设置水文测点垂线处进行枯、洪（或风）季大、中、小潮全潮同步水文泥沙测验；（3）模型试验范围内，设置潮位站不应少于2个，在洪（或风）、枯季全潮水文泥沙测验期间进行半个月以上的潮位连续观测；（4）水文测点垂线和潮位观测内容、方法和要求可按试验研究要求确定。4.1.2.2有径流汇入的模型应包括河道洪水期和枯水期在内的水文泥沙资料，以及相应于全潮同步水文测验期的径流量、输沙量资料。4.1.2.3在河口或海域有水流分汊时，模型试验范围内应有断面流量及流量分配资料。并应有大潮或中潮期的涨、落潮流路和相应的潮位观测资料。4.1.3应有工程设计有关资料及供分析研究用的泥沙资料。条文说明71 4.1.1～4.1.3为潮流定床模型试验所应具备的基本资料及其应具有的基本要求。其中水文测点垂线数量、测点垂线布置对于地形复杂的海岸与河口或重要的工程而言为下限要求，一般在河道（或海域）的上、下游断面和工程区域内应布置有测点垂线，以能反映河道（或海域）潮流场的变化。4.2模型设计4.2.1潮流定床模型应满足几何相似、重力相似和阻力相似条件。相应比尺应分别按下列公式计算：平面比尺（4.2.1-1）垂直比尺（4.2.1-2）流速比尺（4.2.1-3）糙率比尺（4.2.1-4）水流时间比尺（4.2.1-5）流量比尺（4.2.1-6）潮量比尺（4.2.1-7）4.2.2模型设计应符合下列限制条件：（1）模型水流雷诺数（Re）m应大于1000，水流应处于阻力平方区；（2）模型水流应避免表面张力的影响，试验浅滩段最小水深应大于2.0cm；（3）模型糙率不宜大于0.030。4.2.3模型范围应按下列原则确定：（1）在单边界生潮模型中，模型应有过渡段和试验段，在双边界生潮模型中，试验段两侧均应有过渡段，过渡段长度应保证试验段水流相似；71 （2）模型试验段范围应根据试验目的要求和现场潮流具体情况确定，当试验段有工程建筑物时，海岸范围的宽度和长度宜大于3倍建筑物凸出部分长度。4.2.4模型比尺的确定应满足下列要求：（1）海岸与河口潮流模型宜为变态模型，模型变率可取3～10；（2）根据模型范围、试验目的和要求、试验厅（场地）大小和布置，确定模型平面比尺，模型平面比尺宜在1000以内；（3）根据模型相似准则和限制条件、仪器测验精度等，结合电网供电能力、生潮设备生潮能力、径流设备供水能力、水库（蓄水池）蓄水量及模型运转费用等因素，进行综合比较，确定垂直比尺。4.2.5模型糙率及加糙方法按下列条件确定：（1）根据现场水文测验资料及水深图，计算原型床面糙率，由糙率比尺确定模型糙率；（2）模型加糙采用一定粒径的沙、砾、卵石或其它材料粘于模型上，加糙方法宜采用密排加糙和有间距加糙等方式。4.2.6模型潮汐控制的方式及使用范围按下列要求确定：（1）在整体潮汐模型中下边界用水位控制，上边界宜用扭曲水道模拟潮区界段的长度和容积；（2）在试验段距离较长时，上、下边界可采用水位控制，这种形式有潮波反射影响，模型需谨慎调节与操作；（3）在试验段较短、局部工程需放大时，宜采用下边界水位控制，上边界流量控制；（4）在河口外试验工程区域有旋转流时，宜采用下边界流量控制，上边界用流量或水位控制；（5）多口门河道、开敞海域及多边界的控制，可参照上述四种方法进行。4.2.7模型平面布置应满足下列要求：（1）模型边界条件应与天然潮流情况相吻合，应根据工程要求、现场潮流方向、边界情况和模型试验场地、试验设备等具体情况采用合适的生潮控制方式进行单边、双边或多边生潮。（2）当模型范围仅为河口一段或靠岸一侧，涨、落潮为往复流时，应采用双边界生潮控制。71 （3）当模型试验范围处于海岸或河口一部分且流态复杂时，应扩大模型范围或模型进行三面边界流量调节模拟。（4）在模型生潮进出口处和河口边界处应满足模型进出口段水流平顺、潮位变化连续的要求。模型生潮方向应与天然涨、落潮潮流方向基本一致，若不一致可采用人工方法进行调整，需要调整的角度不宜大于15°。（5）在有径流汇入的河口和海岸，在模型中应制有汇入河口一段距离的水深地形，余下河段可采用扭曲河段模拟纳潮量，并在扭曲河段末端施放径流，也可采用流量控制方法控制其纳潮量和径流量。（6）当模型中有径流下泄或余流较强时，在模型边界外应配置水量平衡调配管路和控制系统，及时将多余水量调出到模型外水库（蓄水池），保证模型试验用水量的正常循环。条文说明4.2.2本条为模型设计限制的条件。（1）天然潮流流动形态为紊流，其判别流态的标准是雷诺数（Re）m应大于1000，而要求模型与原型流态基本保持一致，则应保证模型水流处于阻力平方区内；（2）正态模型糙率总是小于原型糙率，为原型糙率np的分之一，变态模型的糙率，它随模型平面比尺λ1和变率η变化而变化，实践经验证明，当0.014＜nm＜0.03时，模型加糙容易实现，因此模型平面比尺λ1和变率η的选取受到一定的限制。4.2.3本条为确定模型范围的原则。（1）在模型水流进（出）口与试验段中间应加有过渡段，使模型水流经过渡段调整到试验段时，能与天然相似，以保证模型试验段水流运动相似。过渡段的长度应根据模型选择的生潮方式、模型潮流进（出）口段宽度及扩展程度、模型生潮方向与潮流间夹角大小等多种因素综合确定。在开敞式海域模型，边界人为处理的较多，为便于模型边界情况的调整过渡段应较长。若过渡段过短，将会导致试验成果的不可信或错误的结论。（2）模型试验段的范围与模型试验方案的范围，设计工程对水流影响的程度、模型流速大小、工程附近的自然条件等因素有关，本条提出的要求是按试验影响范围最小的原则确定的，使用时应按具体情况确定。71 4.2.4本条为确定模型几何比尺的原则。（1）海岸与河口潮流模型要求试验范围都很大，而水深都较小，由于试验场地、供水、供电能力限制，一般要求水平比尺都很大，如采用正态模型，则模型水深就很小，难以保证模型水流流态的相似，因此模型一般做成变态。模型变率应随研究问题的性质、研究范围大小、模型加糙的可能性等因素而变化。根据国内外经验，定床模型一般取3～10，其中局部模型可取较小值，整体模型可取较大值。范围较大的河口、海域模型变率可超过10。在条件容许的情况下，模型变率应尽可能小一些。（2）模型平面比尺应根据试验要求和场地的可能性，并兼顾试验的基本设备如供水、供电能力及人力、物力、财力的消耗来确定。实践表明，局部模型、动床泥沙模型宜取小于600，整体模型宜在1000以内，范围较大的河口、海域模型可超出1000。（3）为保证模型处于阻力平方区内，必须使模型水深不致过小。根据国内外经验，模型垂直比尺宜在100左右，此时模型潮位、制模地形、流速测量精度均能与原型精度基本吻合。若水深很浅，则垂直比尺相应可取小值。4.2.5本条款为模型糙率和加糙的方法。模型水泥浆粉面压光、拉毛或模型床面铺上模型沙，此时模型糙率约在0.014左右；当采用密排或有间距加糙时，模型糙率最大只能达0.03；若模型糙率需再增大则应采用其它方法加糙，这除有一定的难度外并对局部水流有较大的破坏作用，模型设计时应予慎重考虑。4.2.6在模型设计中，应根据模型边界条件、模型场地、纳潮池（水库）水量、生潮固定设备及其布置等条件选用合适的潮汐控制方式。4.3模型制作的准备工作 4.3.1应按模型设计将模型范围、生潮设备、供水循环系统、管路、阀门位置等绘制模型平面布置图。4.3.2应根据模型设计及模型平面布置，编制模型制作施工组织计划。4.3.3应整理与拼接水深图，布设平面控制导线和水深控制断面，摘取模型断面的数据。4.3.4平面控制导线的布设宜采用下列方法：（1）弯曲、分汊河道模型宜采用三角形导线网或平行导线控制；（2）顺直河口和海岸宜采用直导线控制；（3）在主导线两侧应各布置一条辅导线，三根导线应能够覆盖控制模型内全部位置。4.3.5模型地形控制视地形变化复杂程度，可分别采用断面法、桩点法和等高线法，或者几种方法混合使用。选择时可按下列规定进行：（1）当地形变化剧烈、滩槽交错、坡度较大时应采用断面法，断面布设宜垂直于主、辅导线，断面间距宜采用0.5～1.0m；地形复杂地方可适当加密；（2）当地形比较规则、起伏变化不大、岸线比较平顺时，可采用桩点法，间距宜采用1.0m；71 （3）在断面（或桩点）间特殊地形可采用等高线法控制。4.3.6量取导线长度及导线点至各断面距离，其长度应与坐标网格计算值一致，若不一致，应对量取距离按坐标计算值进行修正。4.3.7从导线起始点开始，依次对断面（桩点）编号，分清主导线左右，按图摘取断面距主导线点和辅导线点距离、断面上各点在主导线左右的距离和相应的水深，并作好相应的记录。当断面（或桩点）处无实测水深时，可用相邻等深线或实测点水深内插求得。4.3.8模型断面板可采用三层板或镀锌板等材料制作。断面上应标明断面号、水深点及距离、主辅导线位置，断面裁剪应与断面线平齐，允许偏差为0.5mm。4.3.9当断面处于方案挖深或局部动床处时，断面应做成双层，其厚度由工程方案或冲刷深度要求确定。 4.4模型制作 4.4.1模型场地应进行地基处理，不均匀沉陷不应超过制模精度；模型场地的设置应考虑风、雨等自然条件对试验的影响。4.4.2应在场地上进行导线放样，固定主辅导线点，并进行测量误差计算，使之闭合；在主辅导线间进行断面或桩点放样。4.4.3根据模型长度和范围，在模型范围外便于观测、不易碰撞且不沉陷的地点，设置2～3个固定水准点。试验场地地坪的起始基面应根据模型最大水深或开挖水深加上5～10cm富裕量确定。4.4.4按模型规划断面宽度砌筑模型边墙。边墙应高于模型最高潮位10～15cm。4.4.5断面或桩点应按放样断面线（桩点）位置架设固定。用精密水准仪测量和控制断面（桩点）高程。4.4.6模型内间隙宜采用易密实的砂、煤灰等填料充填并充分压实，按断面（桩点）高程预留水泥沙浆粉面及加糙厚度。4.4.7模型水泥沙浆粉面宜分二次进行，第一次为刮制粗模，第二次为粉面。施工时应使水泥沙浆面压紧，注意施工缝搭接和与边墙的连接，避免漏水、渗水。4.4.8在模型刮制粗模后应对断面或桩点位置和高程进行复核，并做好微地形。制模完成后应对模型断面或桩点平面位置和高程进行校核，并有完整的记录。4.4.9制模断面和桩点高程允许偏差为±1.0mm，平面位置允许偏差为±1.0cm。4.4.10应按模型糙率设计的要求加糙。 4.5模型试验设备4.5.1根据试验场地和试验要求，应选择合适的生潮系统，配置相应的生潮设备、潮水箱或水库。生潮设备的生潮能力可按下式估算：Qm＞（Vmax）m×（hmax）m×Bm＋Q0（4.5.1-1）潮水箱或水库的贮水量可按下式估算：W＞Bm×lm×[（hmax）m－（hmin）m]＋W0（4.5.1-2）式中W——潮水箱或水库的贮水量（m3）；Qm——模型中流量或鼓风机风量（m3/s）；Q0——尾门的泄水量或潮水箱阀门出风量的富裕量（m3/s）；71 W0——潮水箱或水库与供水、回水系统容积的富裕量（m3）；（Vmax）m——模型中最大流速（m/s）；Bm——模型过水断面宽度（m）；lm——模型长度（m）；（hmax）m、（hmin）m——模型中最大水深、最小水深。4.5.2在靠近模型潮水箱或尾门处应设置生潮控制站，模型生潮系统应采用计算机自动控制。4.5.3应根据原型潮位、潮流站观测位置在模型中相应设置潮位仪、流速仪或流速流向仪。4.5.4应根据模型径流量大小和场地具体布置情况，选择不同的供水系统管道和流量控制方式。管道供水流量应大于模型生潮流量，并按供水流量选择供水系统。4.5.5通过模型进入潮水箱生潮池的径流水量，应用水泵或排水控制阀门返回管道（渠道）和水库。4.5.6仪器设备安装完毕，应先单独检查、调试和运转，正常后联合运转。4.6模型验证试验及精度控制4.6.1模型生潮控制站应有给定的边界潮汐水位过程或流量过程，当缺乏此资料时，可用邻近站位资料推算或用数值模拟计算，验证时允许稍微修改边界的控制水位或控制流量过程。4.6.2模型潮汐过程应按水流时间比尺控制。4.6.3应根据原体观测资料进行模型验证试验，内容有潮位、流速、流向、流路和局部流态。4.6.4验证试验应按下列顺序进行，大、中、小潮型的潮位和主流向验证，流速、流向验证；必要的流路验证以及局部流态合理性分析。4.6.5验证试验必须重复进行2～3次，取其平均值作为成果。该成果均以图表表示。4.6.6模型验证偏差换算成原型数值应满足下列精度要求：（1）潮位，高低潮时间的相位允许偏差为±0.5h，最高最低潮位值允许偏差为±10cm；（2）流速，憩流时间和最大流速出现的时间允许偏差为±0.5h，流速过程线的形态基本一致，涨、落潮段平均流速允许偏差为±10%；（3）流向，往复流时测点主流流向允许偏差为±10°，平均流向允许偏差为±10°；旋转流时测点流向允许偏差为±15°；（4）流路，应趋向一致；71 （5）潮量，断面潮量允许偏差为±10%。4.6.7当模型验证试验个别测站流速、潮位结果超出允许偏差时，应对比现场实测资料，分析产生偏差的原因，并应采取相应的措施。条文说明现场全潮水文同步测验是多条测量垂线站位同时作业，人员多，操作精度不一，而自然条件复杂（涨、落潮变化及风、浪影响），使某些资料不尽如人意，在潮位观测上，站位间距较大，岸边情况复杂，水准测量难度较大，且水位观测接尺也很困难，精度上不免受到影响。因此，在验证试验前应首先对各站位垂线流速、流向和各站潮位进行整理，然后比较它们之间存在的关系和差别，分析资料的合理性，对于明显不合理的资料应予剔除。基于上述实际情况，模型相似精度控制指标订得相对比较宽松，其中流速项仅指出“憩流时间和最大流速出现的时间允许偏差为±0.5h，流速过程线形态基本一致，涨、落潮潮段平均流速允许偏差为±10%的要求，在实际控制中，各项指标可能做得吻合更好。在验证试验中，试验结果除上述的原因会产生一定偏差外，还会因水流与局部地形不对应产生偏差，此时应对局部地形进行检查校正，以求模型与原型的相似程度更好。4.7模型水流特性试验4.7.1水流特性试验的测点布置应满足下列基本要求：（1）根据不同工程目的和要求在不同地点布设潮位、测点垂线流速站，观测该处在自然条件下代表潮型时的潮位、流速、流向及流态变化过程；（2）应在拟建港池、航道、导堤等工程区域及上、下游断面布设测点垂线流速站，其数量应能反映工程前、后流场变化和工程影响范围；（3）在工程区上、下游不同地点设置潮位站，该站可与验证试验潮位站吻合。4.7.2水流特性试验宜选用实测大、中潮潮型作为代表潮型，观测内容可包括潮位、流速、流向、流路、流态等。4.7.3水流特性试验必须重复进行2～3次，取其平均值作为试验结果，每二次观测高、低潮位偏差小于原型0.10m，潮段平均流速允许误差为±5%。条文说明71 4.7.1在模型验证试验达到与原型相似的基础上，为取得工程区域及其影响范围内流场、潮位情况，并为工程方案实施后流场、潮位变化提供对比的依据，应进行模型水流特性的观测试验。要求潮位、垂线流速站位布局合理，能反映工程前后流场变化和影响范围。4.7.2大潮潮型潮差大，相对水流动力强，工程前后水流变化大，规律性较为明显，因此观测试验宜选用实测大潮潮型作为代表潮型。对于潮差大的地区或重要工程可增加实测中潮潮型的观测试验；在河口段的重要工程可采用洪、枯季实测大、中潮潮型作试验的代表潮型，以反映不同季节不同潮型的潮流场变化。4.7.3为避免水流脉动的影响和观测的误差，提高试验结果的可靠程度，每组试验应重复观测2次，当2次测验结果误差大于10%或超出观测误差范围时，应进行第三次重复测验，并分析产生的原因，取其接近的2次结果平均值作为测验结果。4.8方案试验 4.8.1应按工程设计方案依次进行方案试验，并应根据试验结果进行修改方案试验。4.8.2方案试验观测内容和方法应与本规程第4.7节规定相同。流速站的数量可随方案的位置不同适当增减。4.8.3通过方案试验，应给出工程前后各方案下列结果，并分别用图和表表示：（1）各潮位站潮位变化，包括最高、最低潮位的变化量，潮汐相位的变化；（2）各测点垂线流速站流速、流向的变化，包括各方案涨、落潮潮段平均流速值和变化量，流向变化量，涨、落潮最大流速及变化量；（3）流路流态的变化，包括主要方案工程附近流路或流态的变化；（4）有断面潮量（流量）时，应增加各方案潮量（或流量）的变化。 4.9试验成果分析4.9.1根据本规程第4.7节和第4.8节试验得到的结果，并结合现场潮流泥沙水力特性和地形变化的分析，论证各方案的优劣。通过综合分析，推荐优化方案。条文说明4.9.1潮流定床模型试验，只能得到各工程方案情况下潮流场、潮量及流态的变化情况，而我们需要通过该模型试验了解与确定工程后泥沙运动、淤积和冲刷、以及对工程的影响情况，此时，只能根据其结果和当地泥沙有关资料，地形演变与水文特征，结合工程实践经验做综合的分析判断，确定优化方案和工程影响。因此，对于潮流定床模型试验研究而言，成果分析是非常重要的一步，且可得到事半功倍的效果。5.1基本资料71 5.1.1潮流泥沙模型所需要的制模地形、水流等资料应符合本规程第4.1节的规定。5.1.2泥沙资料应包括下列内容：（1）陆域一年以上的各月径流量及输沙量资料；（2）研究区域每个水边界至少有一条测点垂线与试验区同步的洪、枯季大、中潮潮流、含沙量变化过程；（3）研究区域内应有与水文测验同步的3～5条垂线测点分层含沙量；（4）有条件时，应搜集工程区不同风况下水体含沙量分布资料；（5）研究区域内洪、枯季底沙及悬沙粒径及级配曲线。5.1.3研究区域内针对研究的对象（如港池、航道）应有两次水下地形图资料，其中一次是初始的，另一次是经过一个水文年或一段较明显变化的水文期（如洪季）的。地形图图比为1∶2000～1∶5000；当研究对象范围较大时，地形图图比可放宽至1∶10000；当地形变化较明显时，可采用固定断面测量资料。5.1.4河口地区不同季节、不同潮型的盐度资料。5.1.5推算的风暴潮资料。条文说明5.1.1海岸与河口天然泥沙运动复杂，有悬沙运动、底沙运动或两者兼有，底沙又包括推移质与悬移质中的床沙质，底沙运动是海岸与河口在潮流、波浪作用下的通常运动形式。研究的问题如港口布置、航道淤积、河口及拦门沙整治等都会遇到不同的泥沙问题，故应根据泥沙运动的性质与需解决问题的实际情况，搜集与测量不同的资料。5.2模型设计5.2.1模型设计前，应按下列要求做好准备工作：（1）收集与分析研究区域已有资料，掌握研究区域动力条件和地形冲淤变化规律；（2）分析研究区域泥沙是以悬沙、底沙运动为主或两者兼有的运动形态类型，根据工程任务的要求，确定进行潮流泥沙模型的类别。71 5.2.2潮流泥沙模型设计除按本规程第4.2节的规定执行外，尚应满足下列泥沙运动的相似条件：（1）泥沙起动相似：（5.2.2-1）（2）泥沙扬动相似：（5.2.2-2）（3）泥沙沉降相似：（5.2.2-3）（4）悬沙挟沙能力相似：（5.2.2-4）（5）底沙单宽输沙率相似：（5.2.2-5）（6）悬沙床面变形相似：（5.2.2-6）（7）底沙床面变形相似：（5.2.2-7）5.2.3悬沙冲淤动床模型泥沙运动相似应满足本节第5.2.2条（1）～（4）及（6）项相似要求，各项比尺应按式（5.2.2-1）～（5.2.2-4）及式（5.2.2-6）计算确定。71 5.2.4悬沙淤积定床模型泥沙运动相似应满足本节第5.2.2条（3）、（4）、（6）项相似要求，各项比尺应按式（5.2.2-3）、（5.2.2-4）及式（5.2.2-6）计算确定。5.2.5底沙动床模型泥沙运动相似应满足本节第5.2.2条（1）、（5）、（7）项相似要求，各项比尺应按式（5.2.2-1）、（5.2.2-5）及式（5.2.2-7）计算确定。5.2.6应根据模型比尺和试验要求选择模型沙。模型沙的密度、重度、粒径大小和沉降速度、起动流速等性能指标，应通过预备试验测定。5.2.7模型沙材料的力学性能和几何特性应保持稳定。条文说明5.2.1～5.2.7（1）潮流泥沙模型设计与试验是运用泥沙运动理论，合理分析原型床面演变规律与泥沙运动状态，正确选择模型沙以及良好的运行操作的综合过程，模型可先按已掌握的资料、知识和必要的辅助试验进行初步设计，在验证过程中允许有修改，完成了验证即完成最后的设计。（2）模型设计宜尽可能使（5.2.2-1）～（5.2.2-7）条件都能满足。例如：定床悬沙淤积模型要满足5.2.2条（3）、（4）、（6）款相似要求，但当试验段较长时，模型沙亦应满足5.2.2（1）、（2）两款的相似条件，以使模型中沉降的泥沙为潮流掀扬带动进入淤积区；又如：底沙动床模型应满足5.2.2（1）、（5）、（7）款的相似要求，但底沙亦有临底悬扬跳跃运动的悬沙，为使这部分泥沙运动相似，则需满足5.2.2（3）款泥沙的沉降相似。国内外有些学者提出悬沙模型尚需满足悬浮相似的条件，即λw=λu*（u*—水流摩阻流速），要满足沉降相似及悬浮条件相似，模型需做正态方可，这是不可能的。国内所进行的河口潮流泥沙模型，未按悬浮条件设计，都取得验证成功，故此处以沉降相似考虑。有异重流发生的盲肠河段的港口、航道的泥沙淤积问题，满足水流运动及悬沙运动的相似条件，即能满足异重流相似的要求。5.3模型制作 5.3.1模型制作除按本规程第4.3节和第4.4节的规定执行外，尚应满足本节的要求。5.3.2模型在考虑的动床范围内宜做成两层，第一层为固定底层，第二层为可动的面层。固定底层断面一般为矩形槽、梯形槽或根据节省模型沙需要做成的其它槽形。固定底层应充分考虑模型可能出现的最大冲刷深度，其高程宜低于最大冲刷深度5～10cm。5.3.3模型采用镀锌板或铁丝与角钢制成的断面板搁置在预制的底槽上，然后充填模型沙浸水密实后，刮制成活动床面。71 5.3.4模型动床的范围应比工程需要的范围有所扩大，定床与动床之间有过渡地带，该处可铺模型沙，用小石子盖面，由密到疏，造成定床到动床的过渡段。5.3.5模型的深槽部分应预留充、排水孔，与模型外的进水管和集水槽连接，充、排水孔应按泥沙的级配设置过滤层。5.3.6模型加糙宜按本规程第4.2.5条规定执行，动床部分的床面糙率不足时可在水面或水中加铅丝等物加大糙率，其值由预备试验确定和模型中调整。5.3.7底沙动床模型在模型尾端应设置集沙槽，悬沙冲淤动床模型在模型回水槽的首端应设置沉沙池。 5.4模型试验设备 5.4.1模型试验设备除应按本规程第4.5节的规定执行外，尚应满足本节的要求。5.4.2模型试验段的两侧宜安置水平导轨。5.4.3根据模型沙的粒径大小、运动形态等因素选用测沙仪、地形仪以及合适的加沙方式和加沙设备等。 5.5模型验证试验及精度控制5.5.1模型试验前，应先做成定床，按本规程第4.6节至第4.9节规定进行模型验证、方案试验和成果分析，根据定床优化的方案进行潮流泥沙模型试验。5.5.2模型应验证试验段的泥沙运动与床面冲淤变化的相似性。5.5.3模型试验时进水与退水应不破坏地形，当使用尾门控制潮汐时，活动的尾门应缓慢升高或降低；当使用潮水箱控制时，应在管路上安装防止模型内涌波产生的排气阀门。5.5.4模型应根据原型泥沙运动情况，涨潮流时在试验段的下游断面加沙，落潮流时在试验段的上游断面加沙；如以单向输沙为主，可在上游或下游断面加沙；如滩地产沙，可在该区增设加沙断面。5.5.5模型加沙可采用下列方式：（1）悬沙可采用潮水箱混合与水库混合方式，或采用加沙泵、加沙水箱连接喷淋头方式；（2）底沙用播沙机或人工加沙。5.5.6模型加沙量的控制应满足下列要求：（1）悬沙淤积、悬沙冲淤动床模型应根据天然的含沙量过程线，按含沙量及时间比尺估算模型进口加入的沙量，实际加沙量要经过床面冲淤验证后确定；（2）根据地形验证，调整含沙量或输沙量比尺和冲淤时间比尺；71 （3）底沙动床模型应根据实测资料或经验公式计算，预估底沙输沙量，其数量、过程和相应比尺由地形验证后确定。5.5.7模型验证试验应随时测定模型中含沙量及相应的有关参数。5.5.8模型的潮汐控制应符合下列规定：（1）潮流泥沙模型应按水流时间比尺控制潮汐水流过程，冲淤时间应根据泥沙模型类型确定，其中悬沙淤积和悬沙冲淤动床模型按悬沙床面冲淤变形时间比尺控制潮汐作用时间，底沙模型按底沙床面冲淤变形时间比尺控制潮汐作用时间；（2）模型潮型应按试验期确定大、中、小潮潮型与个数，并根据潮位记录资料选择大、中、小潮的混合潮型；（3）依据冲淤验证的经验，可使用大潮、大、中潮或大、小潮的组合，以缩短模型冲淤的试验时间，挖入式港池的试验宜选中潮。5.5.9定床悬沙淤积模型应在试验区域内，分区段测量泥沙淤积厚度，或分区段将淤积物取出，测量其体积或烘干称重，确定淤积量。5.5.10模型的冲淤部位应与原型基本符合，总冲淤量的允许偏差为±20%，单纯淤积的模型允许偏差为±15%。5.5.11当一次试验获得模型的冲淤地形和冲淤量与原型基本相似之后，模型应进行重复试验，两次试验的偏差不得大于10%。条文说明5.5.1本条说明潮流泥沙模型复杂，操作运行较难，试验周期较长，根据国内外经验，一般都先在定床模型上优选方案。定床模型上优选方案是以测量流速变化为依据的，泥沙运动与流速、流态关系密切，流速、流态较好的方案即是较优的方案。5.5.6～5.5.8分别说明模型的加沙量和潮汐控制，这两者互为因果。如果模型潮型采用大、中、小潮的组合，则按模型沙特性，先按冲淤时间比尺或71 确定模型冲淤时间，以调整加沙量使模型冲淤变化的地形与天然相似。此种方式模型试验时间较长。根据若干潮汐河口泥沙模型试验的经验，小潮造床的作用微弱，可以忽略。选用大潮或者大潮与中潮的组合也可达到冲淤地形变化的相似，又可节省时间，故模型可选用大、中潮的混合潮型做试验。当模型冲刷过多时，可增加小潮或选用大、小潮的组合，但此时要估算一个大潮、小潮相当于多少个平均潮的作用，以便决定模型冲淤变形的试验时间。决定冲淤时间后，计算试验的加沙量。每次试验结束后，测量模型地形，并与原型实测数值比较，如果符合第5.5.10条要求，则模型动床验证完成；如果不符合，重新调整模型的加沙量，直至使模型冲淤地形变化与原型的冲淤地形变化相似为止。挖入式港池及引排水口的淤积试验，可只考虑工程布置口门附近含沙量过程线及泥沙级配分布的相似。5.5.10目前国内悬沙淤积模型，总淤积量的误差一般都已达到15%范围内；有冲有淤的泥沙冲淤模型较难，总冲淤量的允许偏差可放宽至±20%，但应掌握工程近区相似精度较高的要求。5.6方案试验及精度控制 5.6.1方案试验应按验证试验确定的水、沙条件进行，通过方案比选，选取较优方案。径流影响较大的河口，确定较优方案后，再按不同的水沙组合条件进行试验，以确定该方案在不同水沙条件下的变化。5.6.2潮汐河口的不同水沙条件应包括中水年、枯水年及丰水年三种，重要的工程还应包括特大洪水年和风暴潮的影响。5.6.3需要时可在定床上进行铺沙试验，定性比较各方案的优劣。5.6.4方案试验时，模型加沙量控制应符合下列规定：（1）一个潮流期断面的底沙输沙量允许偏差为±20%，全试验期总输沙量允许偏差为±5%；（2）测点在一个涨落潮周期内悬沙含沙量的允许偏差为±30%，总输沙量的允许偏差为±5%。5.6.5在模型试验过程中，应随时记录潮位、流量、沙量控制的状态，观察水流与泥沙运动受工程建筑物影响的情况。5.6.6方案试验的测量方法和精度控制应与验证试验相同，选择的优化方案应做重复试验。 5.6方案试验及精度控制 5.6.1方案试验应按验证试验确定的水、沙条件进行，通过方案比选，选取较优方案。径流影响较大的河口，确定较优方案后，再按不同的水沙组合条件进行试验，以确定该方案在不同水沙条件下的变化。5.6.2潮汐河口的不同水沙条件应包括中水年、枯水年及丰水年三种，重要的工程还应包括特大洪水年和风暴潮的影响。5.6.3需要时可在定床上进行铺沙试验，定性比较各方案的优劣。5.6.4方案试验时，模型加沙量控制应符合下列规定：（1）一个潮流期断面的底沙输沙量允许偏差为±20%，全试验期总输沙量允许偏差为±5%；（2）测点在一个涨落潮周期内悬沙含沙量的允许偏差为±30%，总输沙量的允许偏差为±5%。5.6.5在模型试验过程中，应随时记录潮位、流量、沙量控制的状态，观察水流与泥沙运动受工程建筑物影响的情况。5.6.6方案试验的测量方法和精度控制应与验证试验相同，选择的优化方案应做重复试验。71  5.7试验成果分析5.7.1在整理分析潮流泥沙模型试验成果时，若有不合理数据，应进行重复试验。5.7.2潮流泥沙模型应分析泥沙的冲淤变化、各部位冲淤量、平均冲淤厚度和最大冲淤厚度，比较各工程方案的优劣。5.7.3优化的工程方案应指明工程布置特点，比较工程前、工程初期及工程后期引起地形变化和流速、流向及流态的变化。5.7.4应从水流、泥沙运动及地形冲淤变化等方面进行综合分析，提出推荐方案。条文说明5.7.4对潮流泥沙模型试验结果进行综合分析，是试验负责人运用泥沙运动理论，已有经验对试验成果进行科学分析，既有理论性，又有实践性，应认真考虑，反复修改，力求准确并简明扼要，在此基础上提出推荐方案与有关建议。6.1基本资料 6.1.1应具有不同季节（如夏季与冬季，大浪期与小浪期）研究海域水深图和典型岸滩剖面图，测图图比宜采用1∶2000～1∶10000；宜具有海岸演变的历史资料。6.1.2应具有年平均高潮位、年平均潮位、年平均低潮位或全年潮位观测资料。6.1.3有测波站时，应具有当地海域一年以上波浪观测资料，包括波高、波周期、波向及观测点的水深；当缺乏测波站常年实测波浪资料时，应具有当地海域连续3个月实测波浪资料及相应时期一年以上的风资料。6.1.4应具有海岸泥沙颗粒分布及级配资料；有条件时，应具有现场实测沿岸输沙率资料。6.1.5有条件时，宜具有沿岸流流速、流向资料和破波带范围的资料。6.1.6应有工程设计图纸资料。 6.2模型设计6.2.1波浪泥沙模型相似律包括：几何相似、波浪运动相似、波浪对岸滩作用相似、泥沙运动相似和地形冲淤变化相似。模型各项相似比尺分别由下列各式确定：（1）几何相似：71 （6.2.1-1）（6.2.1-2）（2）波浪运动相似：波浪折射、波浪陡度、波浪传播速度相似：λL＝λh＝λH（6.2.1-3）（6.2.1-4）波浪绕射、波浪反射相似：λL＝λl＝λh（6.2.1-5）（3）波浪对岸滩作用相似：浅水波底部轨迹速度、沿岸流速度相似：（6.2.1-6）（4）泥沙运动相似：泥沙沉降速度相似：（6.2.1-7）泥沙颗粒加速度相似：（6.2.1-8）泥沙起动相似：（6.2.1-9）71 泥沙起动波高和起动水深，宜通过波浪槽试验确定，也可通过一些半理论、半经验公式初步确定。沿岸输沙量相似，沿岸输沙量比尺可通过试验或通过由下式确定：（6.2.1-10）（5）床面冲淤变形相似：冲淤时间比尺：（6.2.1-11）6.2.2波浪作用下的海岸泥沙运动是复杂的动力过程，波浪及泥沙运动的各自相似条件很难同时满足，应根据现场资料进行全面认真的分析，针对需要解决的问题，抓住主要相似条件，合理地选取模型沙和模型比尺。6.2.3试验波要素及试验水位按下列方法确定。6.2.3.1代表波向可采用下式计算：（6.2.3-1）式中Hi、Ti、Pi、αi——第i量级波浪所对应的有效波高、周期及出现频率与波向。6.2.3.2代表波高可采用对海岸冲淤演变中起作用的大于某量级的有效波高进行统计，按下式计算：（6.2.3-2）6.2.3.3代表波周期可选用与代表波高相对应的周期。6.2.3.4试验水位可按下列条件选取：（1）潮差不大的海岸，取平均潮位；71 （2）潮差较大的海岸取平均高潮位、平均潮位和平均低潮位三种潮位，以其所占的历时比例，依次进行试验。6.2.3.5在有明显的两个方向输沙时，应分别确定该两个方向的代表波向、波高及其它相应的波要素。6.2.4模型几何比尺应按下列条件确定：（1）模型平面比尺应根据模拟的范围、试验目的和要求、以及试验场地的大小确定，平面比尺不宜大于300；（2）模型垂直比尺应满足模型波高大于2.0cm和模型波周期大于0.5s要求；（3）模型变率根据不同的泥沙冲淤情况，宜取1.5～5.0；（4）宜在波浪槽中进行岸滩平衡剖面验证试验，调整确定模型垂直比尺。6.2.5可选用天然沙或轻质沙做模型沙，但应满足泥沙淤积或岸滩冲淤变形相似要求。6.2.6模型平面布置及边界条件应满足下列要求：（1）依据现场实测资料及所需研究的工程问题，确定模型试验的范围。在模型输沙方向首端、尾端应有过渡段，过渡段长度应大于5倍波长；（2）根据波浪作用下泥沙起动水深，确定模型动床范围。动床模型沙厚度应大于最大冲刷坑深度；（3）按试验采用的代表波向布置造波机，造波机与工程建筑物的距离不宜小于6倍波长；（4）波浪水池尾部及两侧应设消波装置，消波装置与建筑物的距离应大于2倍波长，有条件时模型四周宜设置回流槽；（5）集沙槽（盒）应设置在正式试验段的尾部，走向基本垂直于岸线，并与波向线基本平行；集沙槽（盒）在满足拦截全部沿岸输沙的前提下，形状应尽量窄浅；（6）应绘制模型总平面布置图，图中包括模型的范围、工程布置、有关仪器设备的位置及模型导线等。条文说明6.2.2模型沙选择时应注意下列条件的区别。71 （1）在研究破波带以内沙质海岸的岸滩演变时，可在保证岸滩平衡剖面相似、冲淤部位相似、沉降速度相似、起动水深相似及沿岸输沙量相似等主要相似条件和限制条件下，对模型沙进行选取；（2）研究近岸区以外的港区、深水航道冲淤，尤其当波浪较小、泥沙颗粒较细、水流作用因素较强地区的泥沙运动，宜选用轻质沙作模型沙；（3）在研究实际工程问题中，应抓住主要矛盾，宜选用一种模型沙。6.3模型制作 6.3.1波浪泥沙模型制作应满足本规程第4.3节、第4.4节和第5.3节有关规定的要求。6.3.2在变态模型中，对反射影响较大的建筑物，如防波堤、码头等，宜按正态模型制作；也可在变态的模型建筑物表面增加糙度或增加孔隙率的消能措施。建筑物模型可用金属、木材、塑料、有机玻璃、水泥等材料制作。建筑物模型应做防水处理，浸水后不变形，尺寸允许偏差为±1.0mm。 6.4模型试验设备 6.4.1进行波浪泥沙模型试验所用的主要设备为生波机、加沙机、水泵，测量仪器主要有波高仪、流速流向仪、水位仪等，应按模型设计要求进行布置、安装和调试。6.4.2用于量测的各种仪器传感器不宜破坏波形及流场。 6.5模型验证及冲淤时间比尺的确定6.5.1海域波况的验证应符合下列规定：（1）原始波控制点应设在深水区（h/L≥0.5），或设置在海底坡度不影响波向变化的地点；（2）破波带处沿岸线方向应布置三个及以上测波点，测量破碎波波高；在工程建筑物附近，如口门、航道、堤后波影区或港内应布置测波点；（3）控制点处原始波要素应取试验水位条件下，未放置工程建筑物模型时三次测量的平均值，三次测量的允许偏差为±5%；（4）各控制点的波高平均值与模型试验要求的波高值间的允许偏差为±5%。6.5.2应进行破波带位置的验证，破波带位置应与原型基本一致。6.5.3当进行沿岸流流速的验证时，模型与原型允许偏差为±10%。6.5.4各测点波碎波波高和破波角的验证值与原型允许偏差为±10%。71 6.5.5应在模型中进行岸滩剖面验证试验，进一步调整确定模型波浪要素比尺或模型沙，使之满足岸滩形态相似。6.5.6在岸滩形态基本相似的前提下，通过调整加沙量和加沙分布，在上游加沙点和下游集沙点之间达到岸滩形态基本稳定，且加沙量接近于集沙量的情况下，测定模型中的沿岸输沙量。6.5.7应根据输沙量比尺，计算确定冲淤时间比尺。条文说明6.5.1～6.5.6模型验证应注意下列问题：（1）当试验波型为规则波时，采集的波数不能小于10个。不规则波时，采集的波数应多于100个，并给出实测波要素的特征值，一般以H1/3大波作为考核统计波高。数据采集的时间间隔Δt一般取（1/10～1/20）（为平均周期）。（2）破碎波高和破波角的模拟，是岸滩剖面形状和沿岸输沙量、输沙分布模拟的技术关键。在模型验证试验中除了注意模拟破碎波高、破波角外，还应注意破碎水深及破碎点位置的相似。（3）在动床模型中，用初选的代表波要素进行验证试验，目的是在岸滩相似的前提下，进一步调整、确定波要素比尺及模型沙。试验中应注意模型沙加沙量和加沙分布调整同步进行。6.6方案试验6.6.1方案试验应分别测取下列各项资料：（1）各方案在代表波作用下，不同年限的冲淤地形；（2）港池、航道、堤坝等关键工程部位不同冲淤年限对应的冲淤量、冲淤强度；（3）各方案风暴冲淤变化地形，以及各关键部位的骤淤（冲）量和骤淤（冲）强度。6.6.2应通过资料分析，检验各方案的合理性和技术上的可行性，确定试验方案并进行试验。6.6.3各方案试验应重复两次或两次以上，每次试验结果的允许偏差为±20%。71 条文说明6.6.1～6.6.4方案试验中应注意下列各点：（1）完全复演验证试验中的原始控制波浪要素，以及输沙率率定的加沙量和加沙分布；（2）在试验开始，先用较小的波浪作用，使模型沙密实后，再逐渐加大到所要求的波浪。经加沙且状态稳定后，再计算冲淤时间，测取试验大纲所规定的各项资料。为再现试验过程，主要组次应摄影、录像；（3）为消除偶然性和随机性进行的重复试验，地形应恢复到原始状态，在起始条件相一致的条件下进行比较。6.7试验成果分析 6.7.1应通过对所有试验资料的系统分析，得出各方案相应成果。包括不同年限冲淤地形图，港池、航道、堤坝等关键工程部位不同年限对应的冲淤量、冲淤强度以及冲淤变化趋势曲线图，风暴骤淤（冲）量和骤淤（冲）强度。6.7.2应指明各工程方案的特点，从地形冲淤变化和波浪场分布等方面对各方案进行综合比较，提出优化方案。 7.1基本资料 7.1.1船闸模型制作资料应包括船闸平面布置图、纵剖面图、横断面图、细部大样图、上、下游引航道地形图。7.1.2水位资料应包括潮位和内河水位变化，双向水位差组合，河侧与海侧最高、正常、最低的通航水位。7.1.3含盐度资料应包括海侧和河侧水质中的含盐度及其变化。7.1.4船舶资料应包括船舶吨位、船舶型线图、船舶排列组合及型式。7.1.5船闸输水系统设计的水力计算及分析资料。 7.2模型设计7.2.1模型为正态模型，比尺不宜大于25。7.2.2模型应满足几何相似、水流运动相似、动力相似和盐度相似。各项比尺应由下列各式确定：平面比尺（7.2.2-1）71 垂直比尺（7.2.2-2）流速比尺（7.2.2-3）流量比尺（7.2.2-4）时间比尺（7.2.2-5）力的比尺（7.2.2-6）力矩比尺（7.2.2-7）重量比尺（7.2.2-8）糙率比尺（7.2.2-9）盐度比尺λs＝1（7.2.2-10）7.2.3模型比尺除应满足相似条件外，应根据试验要求、试验场地、供水能力、试验时间等条件确定。7.2.4模型设计应符合下列规定：（1）模型水流处于阻力平方区；（2）水流异重流雷诺数应大于1×103；（3）廊道出水孔段水流雷诺数应大于5×104。7.2.5船闸模型应包括上、下游引航道，其长度应大于最大设计船舶（队）长度的3倍。7.2.6船闸模型上、下游平水槽折线型长度，应由平水槽溢流厚度和船闸灌（泄）水过程中最大瞬时流量、按薄壁堰流公式计算确定。7.2.7船闸模型总体布置图应包括船闸、量水系统和仪器操纵系统。71 条文说明7.2.4为模型设计的限制条件，其一，因原型水流为充分紊流，雷诺数的绝对值很大，而模型的雷诺数较小，为了达到模型水流运动与原型相似，即要求模型与原型两者的能线比降与水面比降的比值相等，应控制模型水流在阻力平方区内；其二，廊道出水孔段支孔出流，因受模型缩尺的影响，要求雷诺数应大于5×104，这样可忽略水体粘滞性对出水孔口阻力系数的影响。7.2.5为入海船闸上下游引航道的长度，主要考虑需观测异重流盐水楔入侵闸室时的运动形态，即当盐水楔从上闸首运动到下闸首，再从下闸首返回到上闸首时，它的运动距离相当于盐水楔运动到端部后，尚未返回到闸首时的长度。7.3模型制作与安装 7.3.1应绘制船闸输水系统结构物模型详图、测点布置图，并应提出模型加工及安装要求。7.3.2模型可用塑料、玻璃、有机玻璃、混凝土、木材等材料制作。7.3.3船闸输水廊道关键部位，宜采用便于观察的有机玻璃或玻璃等透明材料制作。7.3.4船闸结构部件加工尺度允许偏差为±0.5mm；船闸模型高程允许偏差为±0.5mm；船闸轴线平面位置允许偏差为±0.1°；引航道地形高程允许偏差为±1.0mm。7.3.5船闸模型应避免漏水。7.3.6船模外形轮廓制作应符合型线要求。 7.4模型设备和量测仪器7.4.1船闸模型试验用水由循环供水系统供给。当试验水质为咸水时，泵房抽水设备应有耐腐蚀性能，平水塔的平水槽应用钢筋混凝土结构，管路系统应选择聚氯乙烯管等。7.4.2上、下游引航道端部应安装能升降的四个支点结构组成的折线型平水槽。7.4.3船闸咸淡水输水系统的阀门、溢流堰闸门、横拉闸门的启闭设备可采用蜗轮蜗杆式、伞齿轮式、链条式传动，由步进电机带动、驱动电源控制。7.4.4防咸船闸的咸淡水交界面宜采用盐度跟踪系统与微机数据采集处理系统控制。7.4.5船闸模型试验采用的常规量测仪器，可按本规程附录A选用。7.4.6船闸模型试验可采用弹簧拉力仪或全环电阻式拉力仪等设备测量缆绳拉力。71 条文说明7.4.2平水槽是船闸水工模型试验的专用设备，功能是保持船闸灌（泄）水时，上、下游水位不会有大的升降。平水槽应为四个支点，实践运用经验说明，采用单点或二点式，不能保持水槽水平，用四点式能保持最大的溢流长度和最佳的平水效果。7.5模型试验内容与方法7.5.1当船闸为单一水质时的模型试验内容与要求，应按现行行业标准《船闸设计规范》（JTJ262）的有关规定执行，试验方法按现行行业标准《内河航道与港口模拟技术规程》（JTJ/T232）的有关规定执行。7.5.2当船闸为两种水质时，应观测船舶过闸时咸水入侵运动形态、掺混情况、沿闸室水深方向的盐度分布和变化。7.5.3船舶（队）进出闸室时，应观察船舶水平运动和螺旋桨搅拌水流引起咸淡水交界面的不稳定和水质掺混。7.5.4应选择船舶（队）在闸室最不利的停泊位置，观测不同船舶吨位、过闸时间、水质盐度、缆绳的拉力。对置换法防咸船闸还应观测非对称停泊时的缆绳拉力。每测次要重复三次。7.5.5船闸与防咸船闸的基本参数，如输水系统阻力系数、廊道转弯处的压力、廊道出口处消能段的流速分布、置换法双层底板上的剩余压力等，应由船闸单一水质的恒定流试验确定。7.5.6船闸模型试验前准备工作应包括调节海侧与河侧中的含盐浓度，校核海侧、河侧水位高程，校核船舶（队）的吃水深度和排水量，校核阀门启闭速度和率定缆绳拉力仪等。7.5.7观测防咸船闸船舶（队）的缆绳拉力时，仪器应平置，测杆应垂直，船舶不得水平移动。7.5.8置换法防咸船闸试验应观测各种水质盐度、船舶大小和不同置换时间时的闸室盐度变化。7.5.9潜没咸水厢法防咸船闸试验应观测咸水厢不同的升降速度、船舶进出咸水厢和闸室的盐度等变化。7.5.10气幕法防咸船闸试验应观测不同供气量、开闸时间、船舶吨位和闸室盐度等变化。71 条文说明7.5.7原型中船舶（队）在闸室停泊是没有位移的，因此在船闸模型试验时也不应有位移，故不宜用弹簧拉力仪测力，它若有水平移动，则影响船侧咸淡水交界的高差，使船只受力值偏小。7.6方案试验7.6.1根据设计要求和基本资料进行原设计方案试验。对模型试验成果进行分析和论证，并应符合现行行业标准《船闸设计规范》（JTJ262）有关规定。当试验成果不满足设计规范有关规定时，应进行模型修改和方案优化试验。7.6.2当船舶（队）在防咸船闸中不能安全停泊时，应分析受力性质，提出改善措施。7.6.3防咸船闸应在优化输水系统设计和方案试验的基础上，进行防咸效果的方案试验。条文说明7.6.1～7.6.2本条在进行原设计方案试验时，若一项不满足设计规范要求，这里指的是输水时间、廊道转弯内侧曲面的压力、引航道流速与涌浪高度、船舶停泊条件等，应进行方案修改试验。7.6.3防咸船闸的方案试验，是指船闸输水系统设计中对加剧水质掺混，影响防咸节淡效果的细部结构及运转方式，应进行方案修改和优化。7.7试验成果分析7.7.1试验成果应包括下列内容：（1）船闸闸室灌（泄）水时的水力特性曲线，包括闸室水位、闸室流量、闸室能量、闸室比能、阀门过水面积、输水系统的流量系数及流量增率与灌（泄）水历时的关系曲线；（2）船闸闸室灌（泄）水过程中，相应阀门各开度时的门槽与转弯半径内侧曲面压力分布，及置换法防咸船闸剩余压力值；（3）船闸闸室灌（泄）水过程中，船舶（队）相应于不同水质盐度和不同停泊位置时的缆绳拉力过程线；（4）气幕法防咸船闸的不同水质盐度、供气量、船舶（队）排水量与咸水入侵量的关系曲线；71 （5）置换法防咸船闸的不同水质盐度、置换时间及船舶（队）排水量与咸水入侵量和淡水损失量的关系曲线；（6）潜没咸水厢法防咸船闸的不同水质盐度、咸水厢升降速度及船舶（队）排水量与咸水入侵量和淡水损失量的关系曲线。7.7.2成果分析应符合下列规定：（1）用船闸运转安全技术指标衡量输水系统布置、灌（泄）水过程中的各项水力特性；（2）在防咸船闸中，船舶（队）在闸室的停泊标准应按现行行业标准《船闸设计规范》（JTJ262）的规定执行；（3）论述船舶（队）过闸时，两种水质交替变化的物理现象，影响掺混的因素；（4）提出防咸船闸的防咸率与淡水损失率；（5）经综合分析，提出推荐方案。条文说明7.7.2防咸船闸试验成果分析应注意下列的问题：（1）防咸船闸在运转过程中，两种水质交替出现，要注意发生掺混的物理现象，一是交界面掺混的控制因素，密度佛氏数Fr＝1；二是掺混层厚度与交界面升降速度有关，宜控制在0.5～1.0cm/s；三是异重流孔口交界面上吸，它与孔口单宽流量有关；四是船舶底部包水现象，它与交界面升降速度有关。（2）防咸船闸除应满足运转安全技术指标和船舶停泊标准外，还应满足防咸率和淡水损失率的指标要求。由于防咸方法不同，指标是有差异的。根据国内外有关防咸船闸的试验研究，防咸船闸的防咸率和淡水损失率分别是：置换法为95%和15%～30%；潜没咸水厢法为93%～95%和5%～7%；气幕法各为50%；集咸坑法为95%和140%～150%。8.1基本资料 8.1.1基本资料除满足本规程第5.1节的规定外，尚应补充工程海区以外计算域的水深图或海图、底质粒径分布图和其它水文泥沙资料。 8.2基本方程71 8.2.1潮流运动可按下列方程控制：（1）连续方程：（8.2.1-1）（2）x向动量方程：（8.2.1-2）（3）y向动量方程：（8.2.1-3）式中ζ——相对某一基面的水位（m）；h——相对某一基面的水深（m）；Nx——x向水流紊动粘性系数（m2/s）；Ny——y向水流紊动粘性系数（m2/s）；f——科氏参量；fb——底部摩阻系数。8.2.2悬沙输移扩散可按下列方程控制：（8.2.2）式中Dx——x向悬沙紊动扩散系数（m2/s）；Dy——y向悬沙紊动扩散系数（m2/s）；Fs——源汇函数（kg/（m2·s））。71 8.2.3床面冲淤变化可按下列方程控制：（8.2.3）式中△h——冲淤厚度（m）；qx——x向底沙单宽输沙率（kg/（m·s））；qy——y向底沙单宽输沙率（kg/（m·s））；γ0——底沙干容重（kg/m3）。8.2.4根据计算方法的不同，亦可选用等价的基本方程。条文说明8.2.1～8.2.4包括潮流运动方程组，悬沙输移扩散方程和床面变化方程。（1）在式（8.2.1-2）、（8.2.1-3）和（8.2.2）中，水流紊动粘性系数Nx、Ny和悬沙紊动扩散系数Dx、Dy的值时空变化不大，均可取为常数，因此扩散项中的Nx、Ny、Dx、Dy可移到偏微分号外。（2）在式（8.2.3）中，同时考虑了底沙和悬沙的影响，但潮流作用下底沙输沙目前研究得很不够，至今尚无有效的计算公式，一般都借用单向水流的研究成果，如Meyer-Peter公式、Englund公式、窦国仁公式和武汉水利电力大学公式进行计算。由于泥质海岸与河口的潮流底沙输沙量所占比重很小，可忽略不计，因此，在式（8.2.3）中，令qx＝0，qy＝0，可满足工程要求。8.3计算模式8.3.1依据计算区地形特征、工程方案等具体情况，可采用ADI法、三角元法、破开算子法、贴体坐标变换法和体积元法等任意一种计算模式进行计算。8.3.2ADI法的计算模式可按附录B执行。8.3.3三角元法的计算模式可按附录C执行。条文说明71 8.3.1计算模式种类很多，按差分网格形状分有三角形、正方形、矩形、四边形、曲线坐标网格及各种形状网格组合等；按计算方法分，有显式、全隐式、半隐半显式；按模拟格式分有ADI法、三角元法、破开算子法、贴体坐标变换法和体积元法等。各种模式均有其实用条件和优缺点，在采用其中任一模式时，应对其稳定性、收敛性、精度和实用效果进行论证和分析。8.3.2ADI法兼有显、隐式优点，是目前我国工程界中最先开展并应用较广泛的一种模式，据不完全统计，我国采用该法进行过海岸与河口工程潮流数值模拟的单位和人员占半数以上，有较广泛的应用基础。8.3.3三角元法原理简单、解法直观、编程容易、节省内存、布设随意，易拟合固边界，有较好的稳定性、收敛性和精度，在海岸与河口二维潮流数值模拟中已得到了广泛应用。8.4计算域的确定及网格剖分 8.4.1计算域的确定应符合下列规定：（1）计算域应足够大，能反映工程海区整体流场特性，并使计算域开边界处的水文要素不受域内工程方案的影响；（2）开边界宜选在流场比较均匀的断面。8.4.2根据计算模式的要求进行计算域网格剖分，网格剖分应符合下列规定：（1）网格剖分后，网格结点水深应能反映水下地形特征和工程前后水深变化；（2）应能概化岸线边界和各工程方案的固边界；（3）网格的疏密应根据计算域内不同部位的工程要求和计算要求确定。 8.5初始条件和边界条件8.5.1初始条件按下列公式确定：ζ（x，y，t）|t＝0＝ζ0（x，y）（8.5.1-1）u（x，y，t）|t＝0＝u0（x，y）（8.5.1-2）v（x，y，t）|t＝0＝v0（x，y）（8.5.1-3）s（x，y，t）|t＝0＝s0（x，y）（8.5.1-4）式中ζ0、u0、v0、s0——分别为ζ、u、v、s初始条件下的已知值。8.5.2固边界应按下列方法确定：（1）法向流速为零：71 （8.5.2-1）（2）法向泥沙通量为零：（8.5.2-2）8.5.3水边界应按下列方法确定：（1）潮流用已知潮位或流速控制：ζ（x，y，t）|Γ＝ζ*（x，y，t）（8.5.3-1）或（8.5.3-2）（2）悬沙按入流和出流情况分别控制：①入流时：s（x，y，t）|Γ＝s*（x，y，t）（8.5.3-3）②出流时：（8.5.3-4）式中Γ——水边界；ζ*——已知潮位（m）；——已知流速（m/s）；s*——已知含沙量（kg/m3）；un——法向流速（m/s）。8.5.4当计算域内存在大面积潮间浅滩时，宜采用动边界技术处理露滩问题。71 条文说明8.5.1初始场的确定与计算潮段的起始时间有关，为提高收敛速度，计算潮段的起始时间宜选在憩流时刻，这时，潮位初始场取起始时间的潮位，流速初始场取零。含沙量初始场依现场含沙量分布情况确定。8.5.2～8.5.4边界条件是数值模拟成败的关键因素之一，因此在确定边界条件时，应对计算域水沙情况和边界特性有充分了解，从而正确确定边界条件。边界条件分固边界和水边界两类。水边界通常应根据已知资料来确定，当水边界处缺乏实测资料时，可依据周围海区已知资料初步选定边界条件，并与计算域内验证点上实测资料进行对比，反复调整初选边界，直至域内验证点上计算结果与实测结果相吻合为止。当计算域内存在大面积潮间浅滩时，处理好边滩处的边界条件是提高模拟质量的重要因素。处理方法可分为两类，一类是变动计算域，另一类是固定计算域。变动计算域法是通过不断地判别，使边界随水边线变动而变动，如水深判别法等。固定计算域法是使边界位置固定而通过边界一定水量的进出以达到相似，如窄缝法等。8.6基本参数的确定8.6.1水流紊动粘性系数Nx和Ny宜由试验确定，也可通过验证计算确定，其值可取50～500m2/s。8.6.2悬沙紊动扩散系数Dx和Dy可取与相应的水流紊动粘性系数Nx、Ny相同数值。8.6.3源汇函数Fs可按下列方法确定：（1）采用挟沙力法按下式计算：Fs＝αω（s*－s）（8.6.3-1）式中Fs——源汇函数；α——泥沙沉降机率；s*——挟沙力（kg/m3），可采用经验公式或试验确定。（2）采用切应力法可通过切应力大小比较确定：71 底部切应力按下式计算：（8.6.3-2）当τ≤τd时，水中泥沙处于落淤状态，则：（8.6.3-3）当τd＜τ＜τe时，床面处于不冲不淤状态，则：Fs＝0（8.6.3-4）当τ≥τe时，床面泥沙起动，则：（8.6.3-5）式中τd——临界淤积切应力（N/m2），宜由试验确定，也可通过验证计算确定；τe——临界冲刷切应力（N/m2），宜由试验确定，也可通过验证计算确定；M——冲刷系数（kg/（m2·s）），宜由试验确定，也可通过验证计算确定。（3）源汇函数也可采用其它等价形式表示，但需经试验验证。8.6.4泥沙静水沉降速度宜通过试验或经验公式确定。条文说明8.6.2在淤泥质海岸与河口，悬沙颗粒很细，含沙量不大，潮流运动紊动粘性系数Nx、Ny与悬沙紊动扩散系数Dx、Dy数值相接近，可采用同一值。8.6.3源汇函数Fs的确定，是悬沙输移扩散计算的关键，由泥沙沉降和起动机理所决定，有多种形式表示和多种方法确定，本规程所列的两种方法是目前应用较广的方法。8.7验证计算及精度控制71  8.7.1验证计算应通过参数和边界条件的调整，使模拟计算结果达到与实测结果基本相符的要求。8.7.2验证计算应包括下列内容：（1）潮位过程线；（2）流速、流向过程线；（3）流路；（4）含沙量过程线；（5）床面冲淤。8.7.3潮位、流速、流向的验证计算精度应符合本规程第4.6.6条的规定，含沙量验证要求变化趋势一致，潮段平均含沙量允许偏差为±30%。8.7.4在验证计算中应对计算与实测的流场与含沙量场合理性进行分析。8.7.5床面冲淤验证计算中，计算与实测平均冲淤厚度偏差应小于30%，并应满足冲淤部位和冲趋势相似的要求。 8.8方案计算 8.8.1应调整工程区边界，进行方案计算。方案计算中有关参数的选取及边界控制条件均应与验证计算的相同。8.8.2通过方案计算，应得出各方案的潮位、潮差、流速、流向、含沙量及床面冲淤变化等结果，并应以相应图表表示。 8.9成果分析 8.9.1应根据工程前后潮位、流场、床面冲淤强度和分布等方面的变化，论证各工程方案的优劣。8.9.2通过对方案的综合分析，应提出推荐方案。8.9.3重要工程的数值模拟应有可视化动态显示。 9.1基本资料 9.1.1除满足本规程第8.1节的要求外，应增加不同水层的流速、流向和含沙量资料。 9.2基本方程9.2.1潮流运动可按下列方程控制：（1）连续方程：（9.2.1-1）（2）x向动量方程：71 （9.2.1-2）（3）y向动量方程：（9.2.1-3）（4）z向动量方程：（9.2.1-4）式中x、y、z——原点置于某一基面，z轴垂直向上的直角坐标系坐标；u、v、w——空间流速矢量沿x、y、z轴的速度分量（m/s）；P——水压力（kg/m2）；Nx、Ny、Nz——分别为潮流沿x、y、z向的紊动粘性系数（m2/s）。9.2.2悬沙输移扩散可按下列方程控制：（9.2.2）式中s——水体含沙量（kg/m3）；Dx、Dy、Dz——分别为泥沙沿x、y、z向紊动扩散系数（m2/s）；ωf——泥沙有效沉速（m/s），ωf＝ω－w。9.2.3床面冲淤变化可按下列方程控制：（9.2.3）71 式中Δh——床面冲淤厚度（m）；qx——沿x向的底沙单宽输沙率（kg/（m·s））；qy——沿y向的底沙单宽输沙率（kg/（m·s））；sb——临底处水体含沙量（kg/m3）；ωfb——临底处泥沙有效沉速（m/s）；Dzb——临底处泥沙垂向紊动扩散系数（m2/s）。条文说明9.2.1～9.2.3给出了在空间直角坐标系中的三维潮流泥沙的基本方程，方程形式为常用形式，亦可写成其它的等价形式。在海岸与河口潮流中，垂向流速w相对于水平方向的流速u、v一般很小，其时间和空间的微分量级要比水平方向流速的微分量级小得多，可以忽略其影响，式（9.2.1-4）可近似写为：（9.2.1-4A）即假定潮流运动服从静水压强分布，这一假定在三维潮流数值模拟中被广泛采用。泥质海岸与河口的潮流底沙输沙量所占比重很小，通常可以忽略不计，式（9.2.3）可写为：（9.2.3A）9.3计算模式9.3.1依据计算域地形特征、工程方案要求等情况，可采用垂向坐标变换法、流速分解法、分层二维法、过程分裂法、边值模型法、破开算子法、谱方法、解析法等任一计算模式进行计算。9.3.2垂向坐标变换法的计算模式可按附录D执行。71 条文说明三维潮流泥沙数值模拟的垂向坐标变换法的计算模式是目前使用最为广泛的一种模式，使用该模式可按附录D执行。在附录D的垂向坐标变换法中采用了潮流运动服从静水压强分布和潮流底沙输沙量可忽略不计的假定。9.4计算域的确定及网格剖分 9.4.1计算域的确定原则应符合本规程第8.4.1条的规定。9.4.2根据计算域地形特征、工程方案计算要求等进行垂向分层和平面网格剖分，平面可采用矩形网格或三角形网格，剖分原则应符合本规程第8.4.2条的规定。垂向分层可根据试验要求、水深及工程性质确定，可分5～10层。 9.5初始条件和边界条件9.5.1初始条件可按下列各式确定：ζ（x，y，t）|t＝0＝ζ0（x，y）（9.5.1-1）u（x，y，z，t）|t＝0＝u0（x，y，z）（9.5.1-2）v（x，y，z，t）|t＝0＝v0（x，y，z）（9.5.1-3）w（x，y，z，t）|t＝0＝w0（x，y，z）（9.5.1-4）s（x，y，z，t）|t＝0＝s0（x，y，z）（9.5.1-5）式中ζ0、u0、v0、w0、s0——分别为ζ、u、v、w、s初始条件下的已知值。9.5.2固边界的边界条件可按下列方法确定：（1）法向流速为零：（9.5.2-1）（2）法向泥沙通量为零：（9.5.2-2）9.5.3水边界的边界条件可按下列方法确定：（1）潮流运动用已知潮位或分层流速（分量）过程控制：ζ（x，y，t）|Γ＝ζ*（x，y，t）（9.5.3-1）或（9.5.3-2）（2）分层含沙量按入流、出流两种情况控制：①入流时：s（x，y，t）|Γ＝s*（x，y，z，t）（9.5.3-3）②出流时：（9.5.3-4）式中Γ——水边界；ζ——相对于某一基面（一般指坐标系原点所在基面）的潮位；ζ*——ζ的已知值；71 u*、v*——u、v的已知值；s*——s的已知值；un——边界法向流速。9.5.4水面边界条件可按下列各式确定：（9.5.4-1）（9.5.4-2）（9.5.4-3）（9.5.4-4）9.5.5床面边界条件可按下列各式确定：（9.5.5-1）（9.5.5-2）（9.5.5-3）（9.5.5-4）式中τx、τy——分别为沿x、y向的分量；h——相对于某一基面（一般指坐标系原点所在基面）的水深。9.5.6当计算域内存在大面积潮间浅滩时，宜采用动边界技术处理露滩问题。9.6基本参数的确定9.6.1Nx、Ny、Dx、Dy宜采用试验或经验公式确定，可取50～500m2/s。9.6.2Nz、Dz宜采用试验或经验公式确定，Dz可取Nz的值。9.6.3τ、τx、τy可按下列各式确定：（9.6.3-1）（9.6.3-2）71 （9.6.3-3）式中、——分别为u、v的垂线平均值；fb——底摩擦系数，fb＝g/c2，。9.6.4τd，τe，M与底质密实度、底质粒径等因素有关，可由试验确定，当缺乏试验资料时也可根据经验选取，最后经验证计算确定。9.6.5有效沉速可由试验及计算确定。条文说明9.6.1潮流水平紊动粘性系数Nx、Ny和泥沙紊动扩散系数Dx、Dy一般来说是时空变量，但其时空变化量很小，可以忽略，通常近似取为常值。9.6.2潮流垂向紊动粘性系数Nz和泥沙垂向紊动扩散系数Dz是三维潮流泥沙基本方程中十分重要的量，Nz可采用Munk型公式、VanRees型公式、混合长理论型公式等经验公式确定，在泥质海岸河口，悬沙颗粒很细，含沙量不大，Dz与Nz量值相近，可取同一数值。9.7验证计算及精度控制 9.7.1验证计算应通过参数和边界条件的调整，使模拟计算结果与实测结果基本相符。9.7.2验证计算应包括下列内容：（1）潮位过程线；（2）分层流速、流向、含沙量过程线；（3）垂线平均流速、流向、含沙量过程线；（4）流速、流向、含沙量的垂向分布；（5）底床冲淤变化。9.7.3潮位、分层流速、流向、垂线平均流速、流向的验证计算精度与本规程第4.6.6条的要求相同；含沙量变化及分布趋势应与实测一致。9.7.4验证计算应对流速场和含沙量场进行合理性分析。9.7.5底床冲淤变化验证中，计算与实测平均冲淤厚度偏差应小于30%，并应使床面冲淤部位及冲淤趋势相似。 9.8方案计算 9.8.1方案计算的内容和要求除满足本规程第8.8节的要求外，还应给出流速、流向、含沙量的空间分布，比较方案的优劣。71  9.9成果分析 9.9.1成果分析内容和要求除满足本规程第8.9节的要求外，尚应分析工程方案实施后的流速、流向、含沙量的垂向变化及其对工程的影响。 10.1基本资料10.1.1应满足本规程第6.1节的要求。其中水下地形图的图比可采用1∶10000～1∶50000。10.1.2应有计算域内突堤、丁坝或丁坝群平面布置图。10.1.3应有计算域内的主要入海河流的输沙量；人工补沙或人工采沙的位置和数量。条文说明10.1.1本条中包括进行波浪泥沙数值模拟时的波浪资料。波浪泥沙数值模拟对波浪资料输入的顺序很敏感，采用定时长期连续波浪观测资料可使数值模拟的结果更近于实际。因此，应尽量采用“长期连续波浪观测资料”。沿岸输沙和岸滩演变计算的时间跨度一般较长，故本条文中要求连续波浪序列的长度一般不宜少于一年。10.2基本方程10.2.1一线模型可按下列基本方程控制：（10.2.1）式中y——岸线位置坐标（m）；x——沿岸线方向的距离坐标（m）；D——计算剖面高度（m），一般取波浪作用下泥沙起动水深处至波浪爬高的最高点间的垂直距离；Q——通过某一断面的沿岸输沙率（m3/s）；71 q——计算域内河流入海或人工采（补）沙引起的进出计算域的单宽输沙率（m3/（s·m））。10.2.2n线模型可按下列基本方程控制：（10.2.2）式中qx——破波带内任一点处的单宽沿岸输沙率（m3/（s·m））；qy——破波带内任一点处的单宽横向输沙率（m3/（s·m））。10.2.3输沙率的计算应符合下列规定。10.2.3.1破波带内总的沿岸输沙率可按下式计算：（10.2.3-1）（10.2.3-2）式中Q——总的沿岸输沙率（m3/s）；P——泥沙孔隙率，可取0.4；Hb——破碎波高（m）；Cb——破碎波速（m/s）；nb——波能传递速度系数，；Lb——破碎波长（m）；hb——破碎水深（m）；αb——波浪破碎时波峰线与岸线的夹角（°），αb恒小于90°；71 αx——波浪破碎时波峰线与x轴的夹角（°）；tanβ——滩面坡度；E——波能，；K1、K2——输沙率系数，宜根据实测沿岸输沙率确定；在缺少实测沿岸输沙率资料时，可取K1＝0.385。当波浪绕射显著时，K2＝1.0K1，当波浪绕射不显著时，可取K2＝0。其余符号意义见第2.2节。注：在以上公式中波高均采用有效波高，有效波高可取平均波高的1.60倍，周期采用平均周期。10.2.3.2n线模型x向的沿岸输沙率和y向的横向输沙率可分别按下列公式计算：（10.2.3-3）（10.2.3-4）式中yb——波浪破碎点的坐标（m）；α——，是底坡；yj-1——第j－1条等深线的y坐标（m）；yj——第j条等深线的y坐标（m）；Qxi——第i断面破波带内总的沿岸输沙率，按式（10.2.3-1）计算；Wi,j——平衡剖面上的第j－1条等深线与第j条等深线之间的距离（m）；C——有速度量纲的系数，可取为3×10－6（m/s）。71 条文说明10.2.1一线模型假定海滩剖面形状不变，在整个计算剖面内滩面坡度不变，海岸的演变表现为岸滩的前进与后退，只存在沿岸方向的沿岸输沙，不考虑垂直于岸线方向的横向输沙。一维沙量平衡方程就是一线模型的基本方程。方程中的输沙率是通过垂直于海岸线的某一断面的总的沿岸输沙率。10.2.2n线模型的基本方程是平面二维沙量平衡方程。方程中的输沙率qx和qy分别是沿岸方向和垂直于海岸线方向的单宽输沙率。10.2.3本条给出输沙率计算公式。（1）常用的一线模型沿岸输沙率公式可分为两类，即流速法和波能流法。波能流法的基本假定是当波浪传入岸边引起破波带内泥沙运动时，总的沿岸输沙与波浪所具有的能量的沿岸分量成正比。行业标准《海洋水文规范》（JTJ213—98）“沙质海岸沿岸输沙率的估算方法”推荐波能流法和流速法。由于流速法计算繁杂，个别参数需根据图确定，不便于海岸演变数值模拟使用。因此，本条规定采用波能流法计算沿岸输沙率。波能流法计算沿岸输沙率的简单公式为：此式将近岸破碎波的波能流的沿岸分量与沿岸输沙率联系起来，适用于计算开敞海岸斜向入射波产生的总的沿岸输沙，为工程界广泛采用。条文10.2.3中式（10.2.3-1）右端的第一项即是。当海岸上修建建筑物后，波隐区波高的沿岸梯度变得显著起来，此时，既要考虑斜向入射波作用，又要同时考虑破碎波波高沿岸不等作用引起的沿岸输沙。Ozasa和Brampton于1980年首先提出了同时考虑斜向入射波和波高沿岸不等的沿岸输沙公式：条文10.2.3中的式（10.2.3-1）即采用了此式。引入71 项，是用来描述建筑物的绕射作用导致的破碎波的沿岸变化对沿岸输沙率的影响。事实上，即使是没有建筑物的开敞海岸，破碎波波高的沿岸不等现象也是存在的，其对沿岸输沙率的影响也是存在的。但是，如果不存在波浪绕射，波高沿岸梯度的变化是随机的，且很小，式中第二项与第一项相比可忽略不计，因此条文中规定：“当波浪绕射不显著时，可取K2＝0”。从理论上讲，破波角沿岸不等对沿岸输沙也有影响，但研究表明，与另外两项比，其影响可以忽略不计，故在条文中未出现。关于沿岸输沙率系数K1、Komar和Inman，根据他们莹光示踪沙实验建议，当用均方根波高时，取K1＝0.77，当采用有效波高时，系数K1应乘以1/2。美国海岸工程研究中心（CoastalEngi－neeringResearchCentre）CERC公式为Q＝7500Pls（yard3/year），相当于输沙率系数K1＝0.87。K2与K1的比值可用来描述建筑物引起的波浪绕射产生的反向水流（Counter-Current）与斜向入射波在建筑物附近海滩上产生的平均沿岸流之间的相对强度。这种平衡取决于沿岸流的横向分布，而影响沿岸流分布的因素众多，如引起波浪绕射的建筑物的离岸距离，建筑物的走向，入射波波高、波周期以及入射角，岸线的平面形态，岸滩坡度等。由于波浪计算和岸线变化模型中有不少近似的和假设，加之现有计算沿岸输沙率公式的不严格，将输沙率系数K1与K2处理为因地而变的参数比较适宜，即通过岸线模型的率定计算来确定K1和K2值。条文中表述为“宜根据实测沿岸输沙率确定”。（2）在n线模型中，不仅需要知道破波带内的总沿岸输沙，而且要知道沿岸输沙在横断面上的分布。关于后者的研究成果不多见。Fulford采用下式表示输沙的断面分布：式中，B是待定常系数，n是待定正整数。上式突出优点是容易积分，并可以选择适当的B，以使积分值等于某一个确定的数值。众多研究表明，取n为3较为合适。泥沙输运还同时发生在静水面以上的滩面上，为能够反映这一特征，以更一般的形式代替上式式中，yb是破碎点的距离；a是常数，a的引进是为了考虑静水位以上波浪增水引起的泥沙运动；c是反映分布曲线宽度的待定常数。下面来确定该方程式中的系数。理论上，泥沙运动是发生在y从0至∞的范围内。即，泥沙运动不仅仅发生在破碎带以内，破碎带以外也会有泥沙的运动。71 令即则可以求出系数a等于破波水深除以底坡，即系数c需根据泥沙输运分布曲线和该方程式采用非线性回归确定。根据实测输沙的沿岸分布，可以推求得c＝1.25，因此，破波带内任一点y处的相对输沙率为破波带内横断面上从y1到y2的相对输沙率为式中Qx[y1，y2]是无维量，乘以整个破波带内的总的沿岸输沙Q即可得到自y1至y2的实际输沙率。条文10.2.3中式（10.2.3-3）即源于此。横向输沙率的计算基于所谓的平衡剖面的概念。平衡剖面是一个理想的岸滩剖面，当岸滩的实际剖面与平衡剖面不一致时，将产生横向输沙。横向输沙的方向取决于实际岸滩剖面的坡度，当实际岸滩剖面比平衡剖面陡时，即产生离岸方向的横向泥沙运动，反之，将产生向岸方向的横向泥沙运动。横向输沙的量值取决于实际岸滩剖面与平衡剖面的差值。因此，两条等深线j－1与j之间的横向输沙为条文式（10.2.3-4）所示。10.3计算模式 10.3.171 根据计算域的地形特征、工程方案等具体情况可选用显式差分或隐式差分模式（见附录E）。 10.4网格剖分10.4.1一线模型应先选定基线。基线宜与x轴重合。然后以Δx为间距对基线剖分，岸线位置的计算点与输沙率计算点相间布置。10.4.2n线模型可采用非正交网格系统。x方向的网格结点位置可固定不变，y方向的结点宜布置在等深线上，等深线随时间变动时，结点的y坐标即反映出水下地形的变化。10.4.3网格点间距应由下列因素确定：（1）模拟岸线的长度；（2）岸线形态的复杂程度；（3）模拟精度的要求。条文说明10.4.3形态比较简单的岸线，网格间距可取100～200m；岸线形态比较复杂或有海岸建筑物，一般取25～50m；也可以在不同的区域采用不同的计算网格。10.5边界条件10.5.1海岸演变数值模拟宜采用计算域两端断面的输沙率作为边界条件。10.5.2边界断面的输沙率可按下列三种情况确定：（1）计算域一端断面有岬头、较长的防波堤、突堤、堤岸或水下深槽，边界外的沿岸输沙在边界处完全受阻不能进入研究区域，即：Qb＝0（10.5.1）（2）边界处海岸线位置不随时间变化，海滩在长时段中处于平衡状态，沿岸输沙率梯度为零，即：Qb＝Qb±1（10.5.2）（3）边界处有丁坝等建筑物，丁坝等拦截部分沿岸输沙，即：71 （10.5.3）式中Qb——有丁坝时的沿岸输沙率（m3/s）；Q0——没有丁坝时的沿岸输沙率（m3/s）；hg——丁坝坝头处水深（m）；hs——波浪作用下泥沙起动水深（m）。10.5.3n线模型边界上整个断面的输沙率宜按本节第10.5.2条规定确定，边界断面上不同网格输沙率可按本规程式（10.2.3-3）确定。条文说明10.5.2说明波浪沿岸输沙率计算的三种边界情况。对于（3）中的情况，边界处有丁坝之类的建筑物，丁坝拦截了部分沿岸输沙，在丁坝坝头处有泥沙绕运，一般可假定破波区内输沙率为矩形分布，绕运输沙量与丁坝坝头以外破波区的宽度成正比。10.6计算波要素及水位的确定10.6.1波浪要素可按下列资料确定：（1）采用不少于一年的长期连续波浪序列。可用现场连续波浪观测资料，也可用根据风资料由波浪预（后）报模式推算的连续波浪资料；（2）采用代表波，代表波的波向、波高和波周期按本规程第6.2.3条要求确定。10.6.2计算水位可按下列条件选取：（1）潮差不大的海岸，取平均潮位；（2）潮差较大的海岸，取平均高潮位、平均潮位和平均低潮位三种潮位，以其所占的历时比例依次进行计算。 条文说明71 10.6.1～10.6.2规定了用于计算波浪沿岸输沙的资料。其中根据风资料用风浪预（后）报模式推算波浪序列时，须用同期的短期实测波浪资料进行验证。计算代表波波高时，可剔除对海岸冲淤演变起作用不大的波高。10.7验证计算 10.7.1应通过参数和边界条件的调整，使计算结果与实测结果基本一致。10.7.2验证计算应包括下列内容：（1）当有沿岸输沙率资料时应进行波浪沿岸输沙率计算；（2）岸滩变形计算。10.7.3波浪沿岸输沙验证计算的结果与实测结果的允许偏差为±20%；岸滩变形验证计算的冲淤部位和冲淤变化总趋势应与实测结果基本一致。 10.8方案计算 10.8.1方案计算中采用的参数应与验证计算的参数一致。10.8.2通过各种方案计算，应提出工程后沿岸输沙沿程分布、岸线演变、岸滩及局部地形变化的趋势，并应绘制有关的图表。 10.9成果分析 10.9.1应分析和比较工程前后岸滩演变情况，工程区域的局部变化及变化原因。10.9.2通过各种方案比较，应提出推荐优化方案。 附录A模型试验设备和量测仪器A.1模型试验设备A.1.1模型试验设备主要有生潮系统、供水循环系统、生波系统、浑水加沙系统、微机系统。辅助设备有流速率定槽、波浪槽、波浪水池、水流水槽、泥沙颗粒分析设备。A.1.2生潮系统应满足下列要求。A.1.2.1生潮系统的型式主要有：（1）气压式潮水箱生潮系统，它是利用密封倒置的箱体存水，由箱顶盖或侧面进气管和出气管控制箱内压力变化，使水压出（涨潮）或流进（落潮）产生模型水位变化；由罗茨鼓风机或叶氏鼓风机、电机、电源开关柜、进风管道、阀门、潮水箱、伺服电机、排水管道、蝶阀与生潮控制仪等设备组成；（2）流量式生潮系统，它通过尾门的提升或下降和阀门的启闭使进入模型的水量变化达到涨落潮的目的。它包括尾门式生潮系统（升降、横拉、翻板型式）、管道式生潮系统和可逆泵式生潮系统等。尾门式生潮系统由水库、水泵、电机、电源开关柜、管道、阀门、尾门、伺服电机、传动杆（齿轮）和生潮控制仪等设备组成。管道式生潮系统由水库、水泵、电机、电源开关柜、进水和回水管道、阀门及启闭传动机构、生潮控制仪等设备组成。71 可逆泵式生潮系统，由可逆水泵或多台小水泵抽水和排水，使模型内水量产生变化，达到涨落潮的目的。它由水库、水泵、电源开关柜、管道或明渠和生潮控制仪等设备组成。A.1.2.2生潮系统的选择按下列要求进行：（1）潮水箱或水库能够利用的水量应大于模型最高潮位和最低潮位变化的水量及在渠道中循环的水量；（2）生潮系统的生潮能力（气压式生潮系统为罗茨鼓风机或叶氏鼓风机的供风量、流量式生潮系统为水泵供水流量、可逆泵式生潮系统为可逆泵供水流量）应大于模型中涨（落）潮水流最大变化速度的要求；（3）气压式生潮系统为试验大厅的固定生潮系统，一般用于生潮口门较宽、潮量较大的边界，其它型式生潮系统一般用于生潮口门较窄的边界，并可随模型生潮边界的变化而灵活布置；当模型为双边生潮时，可根据具体情况选择合适的生潮系统组合使用；（4）采用流量式生潮系统在模型进口段应采取消能措施；（5）采用尾门式生潮系统生潮，尾门宜采用非淹没出流，采用管道式生潮系统生潮，回水管道宜保持非淹没出流；（6）生潮系统的型式选用应考虑供电网路的负荷能力。A.1.3供水循环系统应符合下列规定。A.1.3.1供水循环系统是水流槽、率定水槽和潮流模型流量供给的主要设备，它由水泵、泵房、平水塔、水库（蓄水池）、输水管道、阀门、回水渠道（槽）等组成。A.1.3.2水泵是提取试验用水的供水设备，一般选用离心式水泵，水泵的流量应大于试验用水最大流量，可由一台或多台水泵同时供给。A.1.3.3水库（蓄水池）是贮存试验用水的设备，一般呈圆形或矩形，水库（蓄水池）的最高水位宜低于试验地坪0.5～1.0m，最低水位应保证水泵吸水的浸没深度，贮水量应满足试验用水量（包括循环水量）的要求。A.1.3.4平水塔是保证管道水头不变，管道流量稳定的供水设备。塔高6～8m，容积50～70m3，塔中溢水槽宽度和深度为15cm，溢水槽长度按每0.1m3/s的流量设20m长度计算。当模型需单独供水循环时，可设高2～3m的平水槽代替。A.1.3.5回水渠道（槽）是保证试验用水自流回到水库（蓄水池）的设施，宜置于地面以下，上铺混凝土盖板。回水渠道（槽）布置应保证回水的通畅，且在适当部位加设集水井，使沉积的泥沙和污物易于清理。当模型为泥沙浑水试验时，在回水渠道（槽）前应设置沉沙池（槽）。A.1.4生波系统应符合下列规定。A.1.4.1生波系统是产生波浪的主要设备，由造波板、驱动电机和模拟信号发生器或微机与软件控制组成，按生波板运动方式可分为摇板式（造深水波）、推板式、活塞式（造浅水波）和摆式（造各种波型）；按驱动方式可分为电动液压伺服油缸系统，直流伺服马达与滚珠丝杠系统，变速马达与曲柄驱动系统。A.1.4.2生波机的推板运动振幅与产生的波高可由下式确定：摇板式（A.1.4-1）活塞式（A.1.4-2）式中a——水面处推波板振幅；71 K——波数，K＝2π/L；D——水深。A.1.4.3推板式生波机所需最大功率由下式确定：（A.1.4-3）式中Wmax——生波机最大功率；B——水池或水槽宽度；η——机械效率，可取0.7。A.1.5浑水加沙系统应符合下列规定。A.1.5.1浑水加沙系统是专门用来做模型挟沙水流试验的基本设备。模型浑水加沙有全循环浑水系统和高含沙量水流配制两种形式。A.1.5.2全循环浑水系统是借助潮水箱进气管道气量或水库借助搅拌器搅动模型沙，使水体含沙量在进入模型时达到要求。它不需另外设备。高含沙量水流配制加沙系统是清水进入模型后由高含沙量水掺混，达到水体含沙量的方法，可由泥沙搅拌筒、搅拌电机、进水管、输送浑水管道、喷淋管（头）和控制阀门组成，模型外应设置泥沙沉淀池（槽）。A.1.6微机系统采集、处理试验数据和控制是提高模型试验精度，加快试验周期、实行自动控制的主要设备。它通过潮位、流速、含沙量、压力、应变等传感器进行模/数（A/D）转换输入，经处理后将成果用显示、打印、绘图等方式输出或存储。在生潮、生波、含沙量控制等自动控制系统中，将数字经数/模（D/A）转换送给执行机构产生控制动作，微机系统由微机软件、传感器、输入设备、执行机构、输出设备等组成。A.1.7流速率定水槽是流速量测仪器率定校正的基本设备，水槽分上、下两级，上级装复式量水堰进行流量计量，下级装标准管嘴，由循环供水系统供给流量，在标准管嘴孔口中测定流速量测仪器参数，孔口流速可由堰上水头关系曲线或计算得到，水槽全长约15～20m。A.1.8波浪水槽是研究泥沙（或模型沙）在波浪作用下起动流速，岸滩平衡剖面和防波结构物试验的基本设备，由生波机、消浪器（坡）以及给排水设备组成。槽宽宜60～100cm，槽高应满足模型比尺1∶10的要求，槽长一般应大于10倍波长。水槽两端应设有1∶5～1∶10的碎石斜坡，或设置装有消浪充填物的金属笼消波器，在水槽的一端设置生波机，并装有滤波器。A.1.9波浪水池是波浪泥沙试验的场地，大多为长方形，长宽比可取1.5∶1～2.0∶1，水池可为混凝土结构，底部及池壁应作防水处理，水池尾部及两侧应设消波装置，模型四周应设置回水槽。A.1.10水流水槽是试验研究水流特性和泥沙运动基本规律的主要设备。型式有直槽和环形槽。A.1.11泥沙颗粒分析设备，用于测定泥沙样品的含沙量、各类泥沙样品的颗粒级配和泥沙物理特性。设备主要有烘箱、电热板、感量为0.001g及0.0001g的天平、感量为1g的台秤、振动机、光电颗粒分析仪、标准筛、250cm3与500cm3的比重瓶、温度计、玻璃器皿、玻璃漏斗、干燥器、毛刷、盛沙杯、高压冲洗球等。A.2量测仪器71 A.2.1水位测量：测针、跟踪式自记水位计。A.2.2流量测量：量水堰（矩形堰、三角堰）、电磁流量计、蜗轮流量计。A.2.3流速测量：光电式旋转流速仪、激光流速仪、电磁流速仪、氢泡、摄影等。A.2.4流向测量：水流中流向可用跟踪式流速流向仪；表面流向可采用纸屑、干木屑、泡沫塑料小球测定；河底流向可采用湿木屑、高锰酸钾颗粒等方法测定。A.2.5含沙测量：光电测沙仪、γ-射线测沙仪。A.2.6地形测量：测针、超声波地形仪、电阻式地形仪。A.2.7波高测量：电容式波高仪、电阻式波高仪。A.2.8力的测量：缆绳应力仪、测力计。A.2.9压力测量：测压管，压力传感器。A.2.10盐度测量：电极式盐度计、感应式盐度计、离子选择电化学式盐度计。附录BADI法计算模式B.1差分网格B.1.1基本方程离散时，流速、水位、水深在网格中采用交错排列。B.2基本方程的离散及数值求解B.2.1水流方程的离散及数值求解应按下列要求进行：（1）将一个时间步长Δt分为两个半步长，在每半个时间步长内，依下述求解过程计算潮位及x向、y向流速；（2）在（即前半个时间步长）内，将水流方程离散为下列方程：（B.2.1-1）（B.2.1-2）（B.2.1-3）式中BZ1=171 先将式（B.2.1-1）与式（B.2.1-2）联立，采用追赶法按x方向逐行求解un+1/2和ζn+1/2，然后利用式（B.2.1-3）逐点求出vn+1/2；71 （3）在内（即后半个时间步长）内，将水流方程离散为下列方程：（B.2.1-4）（B.2.1-5）（B.2.1-6）式中BZ2=171 先将式（B.0.1-4）与式（B.2.1-6）联立，采用追赶法按y方向逐行求解vn+1和ζn+1，然后利用式（B.2.1-5）逐点求出un+1；B.2.2悬沙输移扩散方程的离散及数值求解应按下列要求进行。（1）将一个时间步长Δt分为两个半步长，在每半时间步长内，依下述求解过程计算含沙量；（2）在内（即前半个时间步长）内，将泥沙输移扩散方程离散为以下方程：（B.2.2-1）式中71 cx和cy分别为x向和y向上风差分格式与中心差分格式的权系数。对式（B.2.2-1），采用追赶法按x方向逐行求解。（3）在内（即后半个时间步长）内，将悬沙输移扩散方程离散为以下方程：（B.2.2-1）式中对式（B.2.2-2），采用追赶法按y方向逐列求解sn+1。B.2.3地形变化方程的离散及数值求解应按下列要求进行：（1）在nΔt→（n＋1）Δt内，将地形变化方程离散为以下方程：（B.2.3-1）（2）对式（B.2.3-1）显示求解Δhn+1。附录C三角元法计算模式C.1偏导数的求法C.1.1一阶偏导数的求法按下列过程进行：（1）将计算域剖分成若干个任意三角形，对其中的任意一个三角形，都可将水文变量写成下述一阶近似形式：F＝a＋bx＋cy（C.1.1-1）式中F——代表ζ、u、v、s中的任意一个变量；71 a、b、c——为系数。如果已知三角形三个顶点的坐标（x1，y1）、（x2，y2）、（x3，y3）及水文参量的值F1、F2、F3，则可唯一求得式（C.1.1-1）中的三个系数。（2）在任意一个三角形中，水文变量F的一阶偏导数为：（C.1.1-2）（C.1.1-3）（3）在平面域中，如果某一点同为m个三角形的顶点，则该点的一阶偏导数可取这m个三角形在该点的一阶偏导数的加权平均：（C.1.1-4）（C.1.1-5）式中Sei——代表第i个三角形的面积。C.1.2在一阶偏导数求出后，将一阶偏导数作为变量，依与求解一阶偏导数相同的过程，求出水文参量的二阶偏导数。C.2基本方程的离散及数值求解C.2.1水流方程的离散及数值求解应按下列要求进行：（1）在nΔt→（n＋1）Δt内，将水流方程离散为下列方程：（C.2.1-1）（C.2.1-2）（C.2.1-3）（2）依式（C.2.2-1）～式（C.2.1-3），逐点求出ζn＋1、un＋1和vn＋1。C.2.2悬沙输移扩散方程的离散及数值求解应按下列要求进行：（1）在nΔt→（n＋1）Δt内，将悬沙输移扩散方程离散为以下方程：71 （C.2.2-1）（2）依式（C.2.2-1），逐点求出sn＋1。C.2.3地形变化方程的离散及数值求解应按下列要求进行：（1）在nΔt→（n＋1）Δt内，将地形变化方程离散为以下方程：（C.2.3-1）（2）依式（C.2.3-1），逐点求出Δhn＋1。C.3时间步长的确定C.3.1潮流计算的时间步长可按下式确定：（C.3.1-1）式中a——系数，可取1～3；ΔLmin——三角形网格中最小边长（m）；Hmax——计算域中最大水深（m）。C.3.2悬沙输移扩散及床面冲淤计算的时间步长可采用1～10倍的潮流计算时间步长。附录D三维潮流泥沙垂向坐标变换法的计算模式D.1潮流泥沙基本方程D.1.1潮流运动可按下列方程控制：（1）连续方程：（D.1.1-1）（2）x向动量方程：（D.1.1-2）（3）y向动量方程：（D.1.1-3）（4）σ向和z向的速度分量方程：（D.1.1-4）71 （D.1.1-5）式中H——实际水深，H＝h＋ζ；σ——σ坐标，；W——σ向的流速分量。D.1.2悬沙输移扩散可按下列方程控制：（D.1.2）D.1.3床面冲淤可按下列方程控制：（D.1.3）D.2初始条件和边界条件D.2.1初始条件可按下列各式确定：ζ（x，y，t）|t＝0＝ζ0（x，y）（D.2.1-1）u（x，y，σ，t）|t＝0＝u0（x，y，σ）（D.2.1-2）v（x，y，σ，t）|t＝0＝v0（x，y，σ）（D.2.1-3）w（x，y，σ，t）|t＝0＝w0（x，y，σ）（D.2.1-4）s（x，y，σ，t）|t＝0＝s0（x，y，σ）（D.2.1-5）式中ζ0、u0、v0、w0、s0——分别为初始条件下的ζ、u、v、w、s的已知值。D.2.2固边界和水边界的边界条件分别按本规程第9.5.2条和第9.5.3条的规定确定。D.2.3水面边界条件可按下列各式确定：（D.2.3-1）（D.2.3-2）W＝0（D.2.3-3）（D.2.3-4）D.2.4床面边界条件可按下列各式确定：（D.2.4-1）（D.2.4-2）W＝0（D.2.4-3）71 （D.2.4-4）D.3潮流泥沙基本方程的离散及数值求解D.3.1潮流泥沙基本方程可采用有限差分方法离散及数值求解，时间偏导数采用前差，水平方向空间偏导数采用显式差分，垂向空间偏导数采用隐式中心差分。D.3.2采用有限差分方法离散方程式（D.1.1-1）～（D.1.1-4）、（D.1.2）、（D.1.3）可得到下列求解ζ、u、v、W、w、s、Δh的差分方程：（1）潮位ζ：（D.3.2-1）（D.3.2-2）（D.3.2-3）式中i，j——分别为空间网格结点在x、y方向的位置坐标；n——时间层。（2）流速分量u、v、W、w：（D.3.2-4）（D.3.2-5）（D.3.2-6）（D.3.2-7）式中71 i、j、k——分别为空间网格结点在x、y、σ方向的位置坐标；Δσ——垂向空间步长，一般情况下Δσ取0.2或0.1，即垂向分5层或10层。（3）悬沙含沙量S：（D.3.2-8）式中（4）床面冲淤厚度Δh：（D.3.2-9）D.3.3差分方程的求解过程可按下列方法进行：（1）用式（D.3.2-1）显式计算出，用式（D.3.2-4）和式（D.3.2-5）分别隐式计算出和，用式（D.3.2-2）显式计算出，用式（D.3.2-3）计算出，用式（D.3.2-6）和式（D.3.2-7）分别显式计算出和，至此即完成从nΔt到（n＋1）Δt一个时间步长的潮流计算；（2）用式（D.3.2-8）隐式计算出，即可完成从nΔt到（n＋1）Δt一个时间步长的悬沙输移扩散计算；（3）用式（D.3.2-9）显式计算出，即可完成从nΔt到（n＋1）Δt71 一个时间步长的床面冲淤计算。D.3.4水平方向空间偏导数的离散可按下列方法进行：（1）采用矩形网格时：（D.3.4-1）（D.3.4-2）（D.3.4-3）（D.3.4-4）式中F——代表u，v，ζ，W，w，s中的任一量；Δx——x向的空间步长；Δy——y向的空间步长。（2）采用三角形网格时参见附录C。D.3.5时间步长可按下列方法确定：（1）潮流计算的时间步长：（D.3.5）式中ΔLmin——计算域内网格结点间的最短距离；Hmax——计算域内最大水深；α——系数，可取0.5～2.0之间。（2）悬沙输移扩散及床面冲淤计算的时间步长可采用1～10倍的潮流计算的时间步长。附录E波浪泥沙数值模拟离散格式E.0.1一线模型离散格式可按下列方法进行：（1）显式差分格式：（E.0.1-1）式中Δt——岸滩演变计算的时间步长；下标i代表第i个断面；上标n代表第n时步。显式差分格式的允许最大时间步长按下式确定：（E.0.1-2）式中（2）隐式差分格式：71 （E.0.1-3）E.0.2n线模型离散格式可按下列方法进行：（E.0.2-1）式中式（E.0.2-1）可写成下列形式：（E.0.2-2）式中F=Δt/(2ΔxΔh)在式（E.0.2-2）中左端（n＋1）Δt时刻的y值共5个，是待求的未知量，它们前面的系数值和方程的右端项均仅与（n－1）Δt和nΔt时刻的量值有关，故均为已知量。对所有离散结点写出各自的离散方程，构成一个五对角线性代数方程组，辅以上、下游边界条件和海滩的初始条件，即可数值求解。附录F本规程用词用语说明 F.0.1为便于在执行本规程条文时区别对待，对要求严格程度不同的用词说明如下：（1）表示很严格，非这样做不可的：正面词采用“必须”；71 反面词采用“严禁”。（2）表示严格，在正常情况下均应这样做的：正面词采用“应”；反面词采用“不应”或“不得”。（3）对表示允许稍有选择，在条件许可时首先应这样做的：正面词采用“宜”或“可”；反面词采用“不宜”。F.0.2条文中指定应按其它的有关标准、规范执行时，写法为“应符合………的规定”或“应按……执行。” 附加说明 本规程主编单位、参加单位和主要起草人名单 主编单位交通部天津水运工程科学研究所参加单位南京水利科学研究院天津港湾工程研究所河海大学主要起草人吴以喜曹祖德卢汉才罗肇森（以下按姓氏笔画为序）刘子琪刘顺宽李蓓李孟国杨华杨树森吴明阳张征张东生周华兴徐宏明蔡嘉熙71'