CN106951674B - 一种无历史资料小流域洪水重现期推求方法 - Google Patents

一种无历史资料小流域洪水重现期推求方法 Download PDF

Info

Publication number
CN106951674B
CN106951674B CN201710071487.6A CN201710071487A CN106951674B CN 106951674 B CN106951674 B CN 106951674B CN 201710071487 A CN201710071487 A CN 201710071487A CN 106951674 B CN106951674 B CN 106951674B
Authority
CN
China
Prior art keywords
flood
design
years
section
return period
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201710071487.6A
Other languages
English (en)
Other versions
CN106951674A (zh
Inventor
田宏伟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shaanxi Hydrology And Water Resources Survey Bureau
Original Assignee
Shaanxi Hydrology And Water Resources Survey Bureau
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shaanxi Hydrology And Water Resources Survey Bureau filed Critical Shaanxi Hydrology And Water Resources Survey Bureau
Priority to CN201710071487.6A priority Critical patent/CN106951674B/zh
Publication of CN106951674A publication Critical patent/CN106951674A/zh
Application granted granted Critical
Publication of CN106951674B publication Critical patent/CN106951674B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
    • G16Z99/00—Subject matter not provided for in other main groups of this subclass
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
    • Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明涉及一种无历史资料小流域通过洪水频率曲线推求洪水重现期的方法，包括在工作地图确定设计流域边界，在设计流域选取一个或多个设计断面，量算各设计断面面积、河长比降等流域特征值；使用雨量等值线、变差系数等值线计算流域各设计断面设计暴雨；采用推理公式法根据设计断面的暴雨推求设计洪水，计算设计断面的洪水；将各设计断面频率计算成果点绘在海森频率纸上，得到无资料地区小流域各设计断面的洪水频率曲线；根据已知洪水查洪水频率曲线、确定洪水重现期。该方法为下一步计算机编程，便捷的确定设计洪水、确定洪水重现期提供了基础。

Description

一种无历史资料小流域洪水重现期推求方法
技术领域
本发明涉及水利计算技术领域，涉及小流域发生的洪水重现期推求方法，特别涉及一种无资料地区小流域洪水重现期推求方法。
背景技术
小流域是指流域面积小于300km2的流域，现有的无资料地区小流域洪水重现期确定，一般采用试算法，采用水文比拟法或设计暴雨推求断面设计洪水，与小流域洪水比较、多次试算，确定重现期。
水文比拟法，选取同流域或临近降雨产流条件相似、流域面积较为接近的设有水文站流域，将该水文站作为参证站，推求参证站洪水频率曲线，在参证站频率曲线上，查取频率，采用水文比拟法推算到计算流域同频洪水，与小流域发生洪水比较、经多次试算，确定洪水重现期。
设计暴雨推求设计洪水，依据省或当地水文手册雨量等值线成果图、变差系数等值线图推求设计暴雨，采用推理公式法由设计暴雨推求设计洪水，与小流域发生洪水比较、经多次试算，确定洪水重现期。
水文比拟法适用于计算流域与参证站流域面积较为接近、产汇流条件相似的情况，一般来说，部分小流域，难以找到满足计算精度的参证站。因为绝大部分水文站流域面积都在500km2以上，而小流域的流域面积在300km2以下，很多流域面积甚至小于50km2，由于流域面积相差过大，通过水文比拟法计算的设计洪水误差较大，难以反映小流域洪水重现期真实情况。
现有方法采用试算法计算，通过推求设计断面的设计洪水与发生洪水比较，确定发生洪水的重现期，要多次计算洪水频率才能得到与发生洪水相同的洪水，试算法计算较为复杂，步骤多，计算周期长，给防汛抢险指挥带来不便。
现有方法通过推求参证站的洪水频率推求设计站的设计洪水，确定设计站洪水重现期，不是直接制作小流域设计频率曲线。随着计算机软件编程的发展，需要将频率曲线编程，放入计算机软件，根据发生洪水流量，迅速自动查线，得到洪水频率、确定洪水重现期，因此，过去现有的经验方法难以满足新技术发展的需要。
发明内容
为克服上述现有技术存在的缺陷或不足，本发明的目的在于，提供一种无历史资料小流域洪水重现期推求方法。
为了实现上述任务，本发明采用的技术解决方案是：
一种无历史资料小流域洪水重现期推求方法，其特征在于，包括下列步骤：
1)在工作地图确定设计流域边界，在设计流域选取一个或多个设计断面，量算各设计断面面积、河长比降等流域特征值；
2)使用雨量等值线、变差系数等值线计算流域各设计断面设计暴雨，即依据省或当地水文手册成果的年最大1小时、3小时、6小时、12小时、24小时雨量等值线成果图和变差系数等值线图，推求计算各设计断面的2年、3年、5年、10年、15年、20年、33年、50年、100年、150年、200年、……一遇的暴雨；
3)采用推理公式法根据设计断面的暴雨推求设计洪水，计算设计断面2年、3年、5年、10年、15年、20年、33年、50年、100年、150年、200年、……一遇的洪水；
4)将各设计断面频率计算成果点绘在海森频率纸上，得到无资料地区小流域各设计断面的洪水频率曲线；
5)根据已知发生洪水流量，查洪水频率曲线，得到洪水重现期。
本发明的无历史资料小流域洪水重现期推求方法与现有技术相比，带来的有益效果是：
在无资料地区小流域，附近没有降雨、产流条件相似、面积接近的参证水文站情况下，通过设计暴雨推求设计洪水，点绘到海森频率纸，可以直接得到频率曲线，通过查频率曲线确定发生洪水的重现期，而不是传统通过试算确定洪水重现期，因此，本发明的无历史资料小流域洪水重现期推求方法更加直接、便捷。
设计暴雨推求设计洪水过程中使用的雨量等值线成果图、变差系数等值线图，推理公式法计算公式采用的参数，通过收集资料、洪水调查、分析计算、地区综合得到，更能反映流域实际情况，也充分考虑与临近流域的关系，设计洪水成果更加合理，确定洪水重现更加科学准确。
能够直接通过计算、点绘得到无资料地区小流域频率曲线，为下一步计算机编程，便捷的确定洪水重现期提供了基础。尤其是在防洪抢险时需要快速确定某一小流域断面发生洪水的重现期，可以同于预先做好的频率曲线或计算机程序迅速得到，如果按照传统的经验方法计算推求洪水频率、确定重现期，需要2-5天工时，而且计算出的成果还需要审查、修改，总共需要5-7天。采用本发明的无历史资料小流域洪水频率曲线推求方法可以提前制作频率曲线、提前完成审查、修改，预先制图或编程，确定设计洪水、重现期只需不到1分钟时间，大大提高了工作效率，在防汛抢险的紧急关头为缩短决策时间、提高工作效率提供了保障。
充分利用的地区降雨等值线、变差系数、推理公式计算参数等成果，提高了计算精度和工作效率。
附图说明
图1是本发明的一种无历史资料小流域洪水重现期推求方法流程图；
图2为两河皇冠镇区域流域图；
图3为推理公式法计算图，该计算图为两河下游断面，十年一遇(p＝10％)洪水Qm～τ、Qm～t曲线。
图4为两河上游、下游断面PⅢ频率曲线计算成果图。
下面结合附图和实施例对本发明作进一步详细说明。
具体实施方式
按照本发明的技术方案，本实施例给出一种无历史资料两河流域某场洪水重现期推求方法，包括下列步骤：
1)在工作地图确定设计流域边界，在设计流域选取一个或多个设计断面，量算各设计断面面积、河长比降等流域特征值；
2)使用雨量等值线、变差系数等值线计算流域各设计断面设计暴雨，即依据省或当地水文手册成果的年最大1小时、3小时、6小时、12小时、24小时雨量等值线成果图和变差系数等值线图，推求计算各设计断面的2年、3年、5年、10年、15年、20年、33年、50年、100年、150年、200年、……一遇的暴雨；
3)采用推理公式法根据设计断面的暴雨推求设计洪水，计算设计断面2年、3年、5年、10年、15年、20年、33年、50年、100年、150年、200年、……一遇的洪水；
4)将各设计断面频率计算成果点绘在海森频率纸上，得到无资料地区小流域各设计断面的洪水频率曲线；
5)根据已知发生洪水流量，查洪水频率曲线，得到洪水重现期。
本实例给出的两河流域位于陕西省安康市宁陕县，属汶水河上游一级支流，在宁陕县黄冠镇有支流朝阳沟汇入，汇入口以下流域面积205.8km2，河流长度40.5km，河流比降41.4‰，本实施例计算选择汇入口下游作为下游断面，在汇入口上游1km出选择上游断面进行频率曲线推算，分别计算出上游断面、下游断面的频率曲线，经测量某场洪水在该河段下游段面洪峰流量为850m3/s，推求该场洪水的重现期。上、下游断面水文特征值见表1，流域概况见图2所示。
设计断面选取应能够满足重现期查算需要，结合工作需要，在防汛重点位置选取设计断面，当有支流加入流域面积发生明显变化应在支流下游选取设计断面，为提高重现期计算精度、便于查算，也可加密选择设计断面。
表1：两河上游、下游断面水文特征值统计表
断面名称 流域面积(km<sup>2</sup>) 河长(km) 比降(‰)
上游断面 205.8 40.5 41.4
下游断面 129.8 39.5 43.5
暴雨是形成洪水的主要因素，在假定设计暴雨与设计洪水的频率相应和暴雨全面积笼罩的前提下，分析产流和汇流过程，进而推求设计洪峰流量。
1、设计暴雨面雨量计算
①设计暴雨历时
不同集水面积采用不同设计暴雨历时，根据《安康市实用水文手册》，流域面积在300km2以下,采用设计暴雨历时为12h。
②设计点暴雨
根据《安康市实用水文手册》查得各河流设计断面1h、3h、6h、12h、24h的点暴雨量均值Ht和变差系数Cvt，取Cs＝3.5Cv，从P—Ⅲ型曲线模比系数表上查得不同频率相应的KP值，按公式Htp＝KpHt，计算各频率不同历时的点暴雨值HtP。
③设计面暴雨量
各历时设计面雨量由相应点暴雨量乘点面系数而得，点面系数αt的计算公式为：
αt＝1/(1+atF)bt
式中，at、bt表示参数，按《安康市实用水文手册》取值，见表2。
F表示流域面积(km2)。
表2：参数at、bt取值表
设计面暴雨计算公式为：HtF＝αt×Ht
式中，HtF表示某历时t的设计面暴雨量(mm)；Ht表示某历时t的设计点暴雨量(mm)。
设计点暴雨计算成果见表3和表4。
表3：设计暴雨计算表(上断面)
表4：设计暴雨计算表(下断面)
④设计面暴雨量的时程分配
本实施例计算涉及河流，流域暴雨分区属陕南Ⅰ分区，设计面雨量时程按《安康市实用水文手册》分配表确定,时程分配见表5，上断面、下断面降雨过程见表5、表6。
表5：降雨过程成果(上断面)
单位：mm
表6：降雨过程成果(下断面)
单位：mm
2、产流过程计算
根据《安康市实用水文手册》，本实施例计算管线的河流穿越属Ⅰ分区，Im＝100mm，该区域前期影响雨量Pa＝(2/3)Im，确定前期影响雨量为67mm，扣除初损33.5mm，求得产流过程。
将产流过程各时段的产流量累加起来，即为产流总量，乘以“潜流百分数”可得潜流总量，除以相应历时即得平均潜流量，再从单项产流中扣除平均潜流量，即得地面径流过程。
3、汇流计算
本实施例评价的河流面积小于300km2，适合用推理公式法进行汇流计算，推求设计洪水。
①计算Qm～τ曲线
计算流域汇流参数分区属于陕南Ⅰ分区，有关参数按相应经验公式确定根据《安康市实用水文手册》该流域参数按下式计算：
θ＝L/(J1/3F1/4)
将各参数带入后计算，根据计算的θ值确定m的公式为：
m＝0.128×θ0.5927
式中：L表示沿河槽汇流方向至分水岭的河长，km；
F表示集水面积，km2；
J表示与L相应的河道平均比降；
按公式τ＝0.278×L/(mJ1/3Qm 1/4)计算Qm～τ曲线。
②计算Q～t曲线
按公式Q＝0.278F×h/t计算Qm～t曲线。
将Qm～τ及Qm～t两曲线绘于同一坐标系中，见图2，该图为推理公式法计算图，该图为两河下游断面，十年一遇(p＝10％)洪水Qm～τ、Qm～t曲线。
两曲线交点的纵横坐标数值即为所求的洪峰流量和汇流历时，计算成果见表7、表8。
表7：设计断面汇流推理公式法计算成果表(上断面)
序号 频率 τ(h) Q<sub>m</sub>(m<sup>3</sup>/s) 序号 频率 τ(h) Q<sub>m</sub>(m<sup>3</sup>/s)
1 50 12.77 33.5 8 1 5.92 721
2 30 9.89 90 9 0.5 5.63 872
3 20 9.34 119 10 0.3 5.5 951
4 10 7.79 236 11 0.2 5.38 1044
5 5 7.01 363 12 0.1 5.22 1220
6 3 6.54 421 13 0.05 5.06 1365
7 2 6.33 555 14 0.01 4.74 1775
表8：设计断面汇流推理公式法计算成果表(下断面)
序号 频率 τ(h) Q<sub>m</sub>(m<sup>3</sup>/s) 序号 频率 τ(h) Q<sub>m</sub>(m<sup>3</sup>/s)
1 50 12.6 47 8 1 5.76 1110
2 30 9.91 133 9 0.5 5.52 1340
3 20 9.12 179 10 0.3 5.38 1495
4 10 7.64 357 11 0.2 5.29 1592
5 5 6.89 549 12 0.1 5.08 1845
6 3 6.43 723 13 0.05 4.9 2110
7 2 6.13 862 14 0.01 4.59 2742
4、制作上、下断面频率曲线
将计算成果点会在海森频率纸上，得到上断面、下断面频率曲线，见图4。
根据本实施例的设计过程中的经验，点绘应注意以下技术环节：
①针对洪水频率计算，海森频率纸的纵坐标为流量值，单位m3/s。
②应根据计算成果最大值，即10000年(p＝0.01％)的流量，确定海森频率纸纵坐标最大值。
③海森频率纸纵坐标流量均匀排列，应合理布置，每格数据长度宜选择1、2、2.5、5的整数倍，以方便将来查图计算。
④海森频率纸横坐标为频率(P)，单位％，横坐标不均匀分布，不同频率位置由坐标纸给定，不需要重新调整。
⑤当计算同一小流域2个以上断面洪水频率时，宜将各断面频率曲线绘制在同一副海森频率纸上，便于对比定线，也方便今后查算。
⑥同一小流域2个以上断面频率曲线时，应注意上下游关系，一般来说下游断面的设计洪水，应大于同一频率上游断面设计洪水。
⑦同一小流域2个以上断面频率曲线成果，应为一组PⅢ频率曲线，下游断面频率曲线位于上游断面频率曲线上方。
⑧实现过程中,应反复对比，合理分析，综合确定频率曲线。
5、根据断面发生洪水流量，查洪水频率曲线，得到洪水重现期。
经测量某次洪水下游断面洪峰流量为850m3/s，查图4，850m3/s流量对应洪水频率为p＝1.8％，洪水重现期为：
1/p＝1/0.018＝56。
则该场洪水重现期为56年。
重现期查算应注意，查频率曲线得到洪水频率，洪水频率取倒数、并取整得到洪水重现期；当已知流量断面与设计断面重合或断面面积差小于3％时，直接使用已知流量查算洪水重现期；当已知洪水断面与设计断面不重合或断面面积差大于3％时，采用水文比拟法计算最近设计断面流量，查算洪水重现期。

Claims (1)

1.一种无历史资料小流域洪水重现期推求方法，其特征在于，包括下列步骤：
1)在工作地图确定设计流域边界，在设计流域选取一个或多个设计断面，量算各设计断面面积、河长比降流域特征值；
2)使用雨量等值线、变差系数等值线计算流域各设计断面设计暴雨，即依据省或当地水文手册成果的年最大1小时、3小时、6小时、12小时、24小时雨量等值线成果图和变差系数等值线图，推求计算各设计断面的2年、3年、5年、10年、15年、20年、33年、50年、100年、150年、200年、……一遇的暴雨；
3)采用推理公式法根据设计断面的暴雨推求设计洪水，计算设计断面2年、3年、5年、10年、15年、20年、33年、50年、100年、150年、200年、……一遇的洪水；
4)将各设计断面频率计算成果点绘在海森频率纸上，得到无资料地区小流域各设计断面的洪水频率曲线；
5)根据已知发生洪水流量，查洪水频率曲线，得到洪水重现期；
6)已知流量断面与设计断面重合或断面面积差小于3％时，直接使用已知流量查算洪水重现期；当已知洪水断面与设计断面不重合或断面面积差大于3％时，采用水文比拟法计算最近设计断面流量，查算洪水重现期。
CN201710071487.6A 2017-02-09 2017-02-09 一种无历史资料小流域洪水重现期推求方法 Active CN106951674B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710071487.6A CN106951674B (zh) 2017-02-09 2017-02-09 一种无历史资料小流域洪水重现期推求方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710071487.6A CN106951674B (zh) 2017-02-09 2017-02-09 一种无历史资料小流域洪水重现期推求方法

Publications (2)

Publication Number Publication Date
CN106951674A CN106951674A (zh) 2017-07-14
CN106951674B true CN106951674B (zh) 2019-07-02

Family

ID=59466144

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710071487.6A Active CN106951674B (zh) 2017-02-09 2017-02-09 一种无历史资料小流域洪水重现期推求方法

Country Status (1)

Country Link
CN (1) CN106951674B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108509695A (zh) * 2018-03-13 2018-09-07 广州地理研究所 一种中小流域暴雨洪水过程线形状推求方法
CN108875119B (zh) * 2018-04-17 2022-05-27 长江水利委员会水文局 基于分期设计暴雨推求分期设计洪水的计算方法
CN111209525B (zh) * 2020-04-24 2020-09-01 中国水利水电科学研究院 一种基于优化推理公式的洪峰流量计算方法及系统
CN112052602B (zh) * 2020-09-18 2022-06-07 福建省水土保持工作站 一种应用于山区小流域设计暴雨量时程分配计算方法
CN112231907A (zh) * 2020-10-15 2021-01-15 浪潮云信息技术股份公司 一种流域设计洪水时空分布推算方法
CN112765852B (zh) * 2021-01-19 2021-11-02 黄河勘测规划设计研究院有限公司 一种用于特小流域的可能最大暴雨计算方法
CN113095597B (zh) * 2021-05-07 2024-03-12 黄河勘测规划设计研究院有限公司 一种适用于无资料小流域的洪水过程设计方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103276686A (zh) * 2013-05-24 2013-09-04 武汉大学 一种梯级水库下游设计洪水的推求方法
CN104615907A (zh) * 2015-03-11 2015-05-13 武汉大学 基于多变量条件最可能组合推求设计洪水过程线的方法
CN105740969A (zh) * 2016-01-21 2016-07-06 水利部交通运输部国家能源局南京水利科学研究院 一种基于数据驱动的小流域洪水实时预报方法
CN105975672A (zh) * 2016-04-29 2016-09-28 河海大学 一种小流域设计洪水计算方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103276686A (zh) * 2013-05-24 2013-09-04 武汉大学 一种梯级水库下游设计洪水的推求方法
CN104615907A (zh) * 2015-03-11 2015-05-13 武汉大学 基于多变量条件最可能组合推求设计洪水过程线的方法
CN105740969A (zh) * 2016-01-21 2016-07-06 水利部交通运输部国家能源局南京水利科学研究院 一种基于数据驱动的小流域洪水实时预报方法
CN105975672A (zh) * 2016-04-29 2016-09-28 河海大学 一种小流域设计洪水计算方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
"基于暴雨推流原理的洪水预报应用系统研发";胡海娇;《中国优秀硕士学位论文全文数据库 信息科技辑》;20150915(第9期);正文第2.2节，第3.2.3节，第5.1节
"青岛市暴雨统计参数特征与规律研究分析";盛茂刚等;《治淮》;20130415;正文第26页中栏第4段

Also Published As

Publication number Publication date
CN106951674A (zh) 2017-07-14

Similar Documents

Publication Publication Date Title
CN106951674B (zh) 一种无历史资料小流域洪水重现期推求方法
CN105714729B (zh) 一种自适应气候变化的水库多变量设计洪水推求方法
CN105678067B (zh) 城市河涌防洪排涝动态预报控制方法及系统
CN105930931A (zh) 一种电力工程造价的管理方法
CN103276686A (zh) 一种梯级水库下游设计洪水的推求方法
CN103106625A (zh) 水库、闸泵群联合抑咸调度方法
CN109740285A (zh) 一种基于三级汇流结构的分布式侵蚀产沙模拟方法
CN106875048A (zh) 基于典型经验模式的应急洪水预报方法
CN106971013A (zh) 山洪灾害小流域设计暴雨洪水计算方法
CN107085658B (zh) 一种山洪灾害成灾时间确定方法
CN106287233B (zh) 一种管网漏失预警方法及系统
CN107401140B (zh) 支流尾闾河段造床流量的确定方法
CN104182634B (zh) 一种梯级水电站联合运行水位控制断面优选方法
CN113095597B (zh) 一种适用于无资料小流域的洪水过程设计方法
Shouyu et al. Fuzzy Optimization Neural Network Approach for Ice Forecast in the Inner Mongolia Reach of the Yellow River/Approche d'Optimisation Floue de Réseau de Neurones pour la Prévision de la Glace Dans le Tronçon de Mongolie Intérieure du Fleuve Jaune
CN106845116A (zh) 一种河系的洪水预报系统
CN103942325A (zh) 一种融合气候分区思想的海陆气候事件关联规则挖掘方法
CN115907536B (zh) 一种基于区间谬误率的水库特征曲线校正方法
CN108491664A (zh) 一种分布式水文预报模型聚合系统
CN105354416B (zh) 一种基于代表电站的流域降雨径流电量宏观预测方法
CN114548609B (zh) 一种面向生态发电供水的水资源协同调配方法
CN103810401B (zh) 一种人类活动影响分离的二维径流还原方法
CN102912760B (zh) 水库下游河道最低通航水位生成方法
CN106871876A (zh) 面向通江湖泊生态需水的水情监测与分析系统及其方法
CN103823936B (zh) 一种确定等流时线的方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant