感潮闸控口门引排水量计算研究.pdf

第4 期2023年8 月DOI:10.19364/j.1674-9405.2023.04.011水利信息化Water Resources Informatization感潮闸控口门引排水量计算研究NO.4Aug.2023毛兴华，陈澄，杨晓斌?（1.上海市水文总站，上海2 0 0 2 32；2.上海市浦东新区水文水资源管理事务中心，上海2 0 0 12 9)摘要：针对上海市平原感潮河网地区闸控河流复杂的水流过程，为获得完整、准确、可靠的引排水量数据，通过对大治河东闸站和三甲港闸站现有监测分析获得的水位、水量、闸门运行工况等数据进行多元相关分析，得到引排水量与水位、引排时间、水闸过水断面面积的相关关系，并进行回归计算。回归计算结果与实测结果的误差分析结果表明：这种相关关系是比较可靠的，相对误差和标准偏差等均在可接受的范围内。说明用这种方法插补计算缺测期间的引排水量具有一定的可行性，可弥补因各种原因导致缺测带来的水量成果缺失的不足。该方法具有一定的普适性，可用于其他闸控河流水量计算分析，也可避免水位流量关系不显著、不稳定带来的困扰。关键词：感潮河流；闸控口门；多元线性回归；水量计算中图分类号：TV66；P332.40引言上海市濒江临海，属于典型的平原感潮河网地区，嘉定、宝山、浦东、奉贤、金山等沿江沿海各区有若干河流通过水闸与长江口、杭州湾等发生水量交换，是上海市出境水的主要通道。掌握各河流的出入水量情况，对研究上海市出境水资源量，实现全市水资源监控全覆盖，深化最严格水资源管理具有重要意义。在水资源分析评价中，获得河流某一断面的水量通常有以下2 种方式：1）通过转子式流速仪、ADCP等监测流速并计算得到水量；2）假设水位与流量存在较好的相关关系，获取监测水位，通过水位流量关系查询获得实时流量，再计算某一时段的水量。方式1可靠性相对较高，但需要开展监测，成本较大；方式2 可以有效降低成本，可对因仪器故障、外部环境干扰等导致缺测的时段进行插补，获得完整的水量数据，但获得稳定的水位流量关系难度较大。近年来，上海市在闸控河流水资源监测方面开展了大量工作，取得明显成效，在建立水位流量关系方面做了许多尝试，基本结论是感潮闸控河道不满足水位流量的单值关系，反映了感潮河流的复杂性。本研究在梳理上海市入海河流及水量监测现状的基础收稿日期：2 0 2 2-12-14基金项目：上海市水务局（上海市海洋局）科研项目（沪海科2 0 2 3-0 9)作者简介：毛兴华（19 7 3一），男，甘肃民乐人，硕士，高级工程师，主要从事河口海岸水文分析、水资源管理、水文信息化等工作。E-mail:文献标识码：A文章编号：16 7 4-9 4 0 5(2 0 2 3)0 4-0 0 59-0 5上，选取代表性测站历年水位、流量、水量资料，分析确定引排水量与上下游水位和水闸工况等影响因素之间的定量关系，为优化低流速、小流量、感潮河流水资源监测，深化入海水量分析评价提供可行的技术思路。1入海河流水资源监测及影响因素分析1.1监测概况上海市大陆片区共有2 4 条主要的入海河流，分别为黄浦江，嘉定区浏河、蒲华塘，宝山区练祁河、新川沙河、顾泾、马路河，浦东新区外环运河、随塘河、赵家沟、张家浜、川杨河、大治河、赤风港、芦潮引河、人民塘随塘河、芦潮港，奉贤区中港、金汇港、航塘港、南门港、南竹港，以及金山区龙泉港和张泾河。目前，黄浦江，嘉定区蒲华塘，宝山区练祁河，浦东新区随塘河、张家浜、川杨河、大治河、赤风港和芦潮引河，奉贤区中港、南门港、金汇港和南竹港等13条河流均开展了水量监测。水资源监测以ADCP在线监测为主，通过监测数据质量控制和定期比测率定获得流量水量成果，在一定程度上保证了成果的准确可靠性。但部分站点误差相对较大，难以得到满足整编要求的监测成果，主要原因在于：1）受潮汐潮流影响，往复流特征明显，60小潮期和涨落潮转换期流速较小，但流态复杂，走航式ADCP在极低流速下的测验精度难以保证，导致难以准确测流；2）大部分河流受水闸启闭的控制，不同工况对流量监测和分析影响较大；3）大多数河流规模较小，水位落差不明显，人为影响因素较多，对流量率定和水位流量关系的建立造成了较大困扰，目前还没有较成熟的成果。1.2影响因素分析影响感潮水闸过闸流量的主要因素包括闸门上下游水位差、水闸开度、开闸孔数等。感潮河网地区在涨落潮转换阶段，有可能存在涨潮落潮流和落潮涨潮流等情形，对流量的影响更加复杂。多年来，感潮闸控河流的流量与影响因素的相关关系一直是水资源监测研究的难点，大多数研究都聚焦于断面流量与上下游水位的关系上，以找出两者的定量关系为目标。仲兆林！开展了常州沿江感潮河道水文站水位流量关系综合定线分析，通过一潮推流法结合图解法开展综合定线分析，取得了较好效果；顾正华等 2 应用人工神经网络理论建立了感潮水闸流量的非线性计算模型，可计算出一次调水的水闸过水总量；万晓凌等 3 以水力学公式为基础，分析水位差流量关系，为沿江水闸的引水量统计提供了比较精确且简易可行的方法；周全等 4-5 利用实测水闸流量数据，分析了在线监测、水文水动力模型和基于水力学多元回归3种方法获取流量水量数据的优缺点；常黎明等 6 探讨了特定流量条件下平板闸门开启高度与上游水位之间的关系，并将计算结果与实际运行测量结果进行对比；高晨晨等7 以上海市龙泉港出海闸为例，基于回归计算方法分析了淹没式堰流水位流量关系，定线精度可达到一类精度标准，并认为可以在受潮汐影响的出海闸水位流量关系率定中推广。21代表站引排水量计算分析2.1基本思路针对闸控口门的流量和水量，通常做法是直接寻求流量和主要影响因素的相关关系，其中最主要的是流量与上下游水位差之间的关系。但诸多尝试表明，在感潮河网地区，这种相关关系是不稳定且不可靠的，流量与水位差的相关系数较低。比较成功的方法是一潮推流法，即在一个涨潮或落潮周期内，在水闸运行工况稳定（开启孔数和开度高度不变）的情况下，流量与上下游水位有相对较好的关系，但该方法普适性较低，难以推广应用。水利信息化从另一个角度讲，水资源分析评价关注的重点实际上是一段时间内河流的水量变化，是流量在时段内的积分。寻求水量与影响因素的相关关系，能够抵消时段内偶然性因素的干扰带来的误差，使相关关系更加稳定、可靠。因此，可以考虑通过闸控口门水文测站引排水期间所测水量、闸内外水位、水闸开度和开启时间等资料，参照水闸流量理论公式，通过回归计算方法，分析水量与各要素的相关关系，得出水量计算的经验公式。同时通过复相关系数和标准偏差判定回归关系的准确性和可靠性，最终确定经验公式计算的可行性。2.2计算方法根据SL247一19 9 9 水文资料整编规范，平底闸半开（淹没）或全开（自由孔流）状态下流量计算公式如下：Q=M,BeVZ,式中：为过闸引水流量；M2为淹没式孔流流量系数；B为闸孔总宽或开启净宽；e为闸门开启高度；Z为上下游水位差。可以看出，水闸流量大小与过水断面面积和上下游水位差有关。将式（1）应用到引排水量计算中，需要考虑引排水时间的长短。同时还要注意，对于感潮河流，闸外是广阔的河口或海洋，水位受潮汐的直接影响，并不会因为引水或者排水而变化。对于闸内水位，由于水域范围相对较小且封闭，不受潮汐的影响，引排水会导致水位的升高或降低。因此在计算过闸水量时，用闸内外水位差是不合适的，用引排水期末和期初闸内的水位差作为计算因子更符合实际情况。引排水量计算公式如下：W=kTasbZ,式中：W为引排水期间的过闸水量；k为水量计算系数；T为引排水时长；S为水闸过水断面面积；Z为引排水期末和期初闸内水位差；a,b，c 为各参数对应的经验指数。为便于Excel 回归计算，对式（(2）取对数后得到：lgW=lgk+algT+blgS+clgZ。式（3）为多元线性关系式，可通过多元回归分析得到k，a，b，c 的回归值和复相关系数，并计算回归水量。本研究以上海市浦东新区2 0 19 年川杨河的三甲港闸（闸内）站、2 0 2 0 年大治河的大治河东闸（闸内）站的监测成果为代表，分析相关关系。这2 处测站都通2023(4)(1)(2)(3)第4 期过ADCP监测获得断面流量，再计算出引排期间的引排水量。2.3计算过程与结果分析2.3.1河流基本情况2.3.1.1川杨河川杨河位于上海市浦东新区南部，全长为2 8.7 km，西通黄浦江，东入长江口，是浦东新区大控制骨干河流和重要的通航河流。川杨河河口宽为54 6 4 m，河底宽为2 0 30 m，河底高程为-1.6 m，边坡比为1：3,常水位在2.52.8 m之间。川杨河西端建有杨思水利枢纽，东端建有三甲港节制闸。后者主要功能为挡潮和引排水，闸孔径为3孔10 m，设计水位闸内最高为3.0 m，最低为2.0 m，闸外最高为6.0 m,最低为-0.5m。2.3.1.2大治河大治河横贯闵行、浦东两区，西起黄浦江，东入东海，是浦东新区主要的排水河流，也是上海市最大的人工河。大治河全长为39.54 km，河面宽为10 2 m，河底宽为6 4 m，河底高程为-2 m，具有航运、灌溉、纳潮泄洪、排涝调蓄、引清调水等功能。历史最高水位为 3.6 5 m。大治河建有大治河西闸和东闸2 座水闸枢纽工程；其中：大治河东闸为6 孔，宽为10 m；北侧孔为通航孔；南侧另设鱼游道1孔，孔径为2 m，孔总净宽为 6 2 m。2.3.2计算过程2019年，三甲港闸站共测得30 0 次引排水流量水量数据，其中排水13次，引水2 8 7 次。2 0 2 0 年，大治河东闸站共测得16 9 次数据，均为排水数据。首先，对数据进行初步的合理性分析，如引水出现内河水位下降，或排水出现内河水位升高，均为不合理样本，应予剔除。根据每一次闸门的开启孔数和高度，计算水闸过水断面面积，并计算引排水期间的水位差和引排水历时，完成数据的初步处理。其次，对引排水量、闸内水位差、断面面积、引排水历时取对数，并按式（3）进行引排水量与其他3个要素的多元线性回归分析，求出k，a,b，c 的回归值和复相关系数，再根据式（2）计算水量，并与原实测水量进行比较。2.3.3结果分析2.3.3.1回归分析经计算，2 0 19 年三甲港闸引水水量W，与引水历时T、断面面积S和闸内水位差Z的相关关系如下：毛兴华等：感潮闸控口门引排水量计算研究方法基本可行。500一2 0%包络线400F叫上/鲁回3002001000图1三甲港闸站回归计算水量与实测水量关系图1200一2 0%包络线10008006004002000图2大治河东闸站回归计算水量与实测水量关系图2.3.3.2误差分析为判别计算引排水量的准确性，分别对单次监测结果，月、年水量的相对误差和标准偏差等进行分析，结果如表1所示。可以看出，个别单次引排水量回归计算结果的相对误差较大，如三甲港闸站的单次相对误差最大为7 9.0%，大治河东闸站的单次相对误差最大为6 2.8%。2 个站单次回归结果的标准偏差分别为15.9%和13.1%，满足SL2471999水文资料整编规范三类精度水文站的精度要求 8 16-18。三甲港闸站和大治河东闸站各月引排水量计算结果如表2 所示。从月度引排水量来看，由于正负61。(4)2020年，大治河东闸排水水量W与排水历时T、断面面积S和闸内水位差Z的相关关系如下：W,-83.31 0990 02362 0.097三甲港闸站、大治河东闸站2 处水闸引排水量与历时、过水断面面积和水位差的相关系数R分别为0.934和0.9 51。实测水量和计算水量的关系分别如图1和2 所示。经计算得出，2 个站分别有8 1.7%和87.0%的计算水量在2 0%误差范围内 8 16-18，说明该50100150200250300350400实测水量/万m3100200300400500600实测水量/万m3(5)70080090062表1三甲港闸站和大治河东闸站引排水量回归计算误差分析%误差计算三甲港闸站大治河东闸站最大值79.0单次相对误差平均值最大值月水量相对误差平均值年水量相对误差单次结果标准偏差月度水量20%误差包络线上线率误差的相互抵消，月度水量的相对误差大幅减小，均在2 5%以内，标准偏差分别为11.9%和7.18%。回归计算得到的年度引排水量的相对误差更小，仅为月份实测水量/万m计算水量/万m1月3075.732月3601.383月4 137.774月3 896.805月9.490.796月6604.987月7216.958月7764.949月6.847.7510月5781.2011月5005.3412月1585.06总计65 008.68本方法具有一定的普适性，可用于其他闸控河流水量的计算分析，同时避免了传统水位流量关系不显著、不稳定带来的困扰。本研究以水位差、断面面积、引排水历时3个因素为自变量，探究了与水量之间的复相关关系。目前初步分析得到的结论是：3个因素都比较敏感，去掉任何1个因素，复相关系数都会降低。每一个因素在复相关关系中的贡献率和敏感性如何，还需要进一步深入分析。参考文献：1仲兆林.常州沿江感潮河道水文站水位流量关系综合定线分析 .广东水利水电，2 0 0 9（11)：33-34,4 1.2顾正华，李荣.感潮水闸流量神经网络计算模型研究 .水利信息化-1.43%和0.9 5%。说明用本研究的方法计算1个时段内（月、年)的引排水量，准确性是能够满足要求的。62.83结语1.410.1924.519.60.380.50-1.430.9515.913.111.97.1881.787.0表2 三甲港闸站和大治河东闸站各月引排水量计算结果三甲港闸站相对误差/%3829.0424.503992.3510.904 539.339.704.471.3514.708726.83-8.056 498.95-1.616 860.60-4.947152.87-7.886.507.47-4.975425.80-6.154.757.07-4.961320.29-16.7064 081.96-1.433万晓凌，陆小明，周毅，等.感潮水闸引水量计算方法研究J.长江科学院院报，,2 0 13,30(4)：17-2 0.4周全，杨笛音，孟秀.平原感潮河网地区边界水闸流量计算方法 .水电能源科学,2 0 14,32（5)：19-2 2.5周全，吴伟华，杨晓斌，等.上海浦东水闸过闸流量和水量率定方法研究 .水利水文自动化，2 0 10（1)：6 6-6 8.6常黎明，胡亚洲,左小益.特定流量下闸门开度与上游水位关系及计算方法 .人民长江，2 0 19,50（增刊2)：17 6-17 8.7高晨晨，谢先坤.基于回归分析的龙泉港出海闸淹没式堰流水位流量关系率定.上海水务，2 0 2 1（2)：2 8-32.8中华人民共和国水利部.水文资料整编规范：SL247一1999S.北京：中国水利水电出版社，19 9 9.2023(4)本研究通过寻求实测引排水量与水位、断面面积、引排水历时等要素的相关关系，得出水量计算的经验公式；对于因各种因素导致监测缺失的时段，只要水位、水闸运行记录完整，就可以利用这种方法进行插补计算，得到相对完整、准确的引排水量，准确把握水资源变化。在水资源分析评价中，一般以月、年水资源变化情况为重点。用本研究提供的方法计算引排水量，尽管某个具体的引排水量计算过程可能存在较大误差，但是汇总到月和年时间尺度，误差相对较小，能够满足水资源分析计算要求。实测水量/万m计算水量/万m35 193.514773.393 611.913563.522.834.272 617.893 028.293211.622521.783015.9111 675.4811993.7616534.6817317.596227.996028.836 492.606482.935500.375336.124 641.734 661.841 801.191 725.8370 063.8170 729.23海洋工程,2 0 0 7,2 5(3)：10 9-114.大治河东闸站相对误差/%-8.09-1.34-7.636.0519.602.734.73-3.20-0.15-2.990.43-4.180.95第4 期毛兴华等：感潮闸控口门引排水量计算研究63Research on water diversion or drainage calculation of sluice controlled gate intidal areasMAO Xinghua,CHEN Cheng,YANG Xiaobin?(1.Shanghai Hydrological Administration,Shanghai 200232,China;2.Shanghai Pudong New Area Hydrology andWater Resources Management Affairs Center,Shanghai 200129,China)Abstract:In order to obtain complete,accurate,and reliable water diversion or drainage data for the complex flowprocess of sluice controlled rivers in the plain tidal river network of Shanghai,this paper conducts a multivariatecorrelation analysis on water level,water volume,gate operating conditions from monitoring and analysis of theDazhi River East Gate Station and Sanjiagang Gate Station,leading to the correlation between water volume and thewater level,the diversion or drainage time,and the cross section area of the sluice.Thus,regression calculation canbe carried out.The error analysis results of regression calculation and measured results indicate that the correlation isrelatively reliable,and the relative error and standard deviation are acceptable.By applying this method to interpolateand calculate the water diversion or drainage amount during measurement absence,it is feasible to compensate for thelack of water yield caused by measurement absence due to various reasons.It is universal for this method be appliedin calculation and analysis of water flow in other sluice controlled rivers,solving the problems of inconspicuous andunstable relationship between water level and discharge.Key Words:tidal rivers;sluice controlled gate;multiple linear regression;water volume calculation(责任编辑：来冰华）(上接第34 页)Application of digital twin Jialu River in flood control during the projectconstruction periodLIU Zihan,LI Yanfeng,LYU Pengju,ZHAO Qianqian(Henan Provincial Hydrology and Water Resources Center,Zhengzhou 450003,China)Abstract:In order to meet the requirements of refined flood control management during the construction period ofJialu River Comprehensive Treatment Project,as well as ensure the project safety while benefiting from flood controlof the river basin,digital twin technology is applied to the accurate risk analysis,warning,and flood control planningduring the projects construction period.The baseplate of digital twin Jialu River is based on self-built collectionof Level-3 data and shared Level-1 and Level-2 data.Digital twin Jialu River is supported by model platform andknowledge platform as the algorithm,and underpinned by informatization infrastructure as the computing power,to build a watershed flood control application with“four pre”functions,such as forecast,early warning,preview andplan.By collecting basic data,establishing digital scenarios and constructing models in the Jialu River Basin,threetypes of flood scenarios are set with rainfal levels exceeding 100-year,approximately 50-year,and approximately 20-year occurrences.Six types of risks during the project construction period are identified,including dike overflow risk,submergence risk,erosion risk,flood backflow risk,bridge water-blocking risk,and out-embankment waterlogging risk,and flood control measures are formulated.The digital twin Jialu River ensures the project safety during flood season in2022,summarizes the types and quantitative evaluation methods of flood risk during the construction period of watercomprehensive treatment project.Thus,it forms a complete solution for terrain surveying,flood simulation,risk warning,and contingency measures formulation under complex environmental conditions of ongoing construction projects.Key words:digital twin;water conservancy project;construction period;flood control;Jialu River(责任编辑：来冰华）