CN106777959A - 人工干扰无水文资料地区河流环境流量分区界定计算方法 - Google Patents
人工干扰无水文资料地区河流环境流量分区界定计算方法 Download PDFInfo
- Publication number
- CN106777959A CN106777959A CN201611139996.XA CN201611139996A CN106777959A CN 106777959 A CN106777959 A CN 106777959A CN 201611139996 A CN201611139996 A CN 201611139996A CN 106777959 A CN106777959 A CN 106777959A
- Authority
- CN
- China
- Prior art keywords
- flow
- unit
- water
- river
- section
- 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.)
- Granted
Links
- 238000000205 computational method Methods 0.000 title claims abstract description 19
- 238000000034 method Methods 0.000 claims abstract description 84
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 76
- 230000007613 environmental effect Effects 0.000 claims abstract description 16
- 238000004458 analytical method Methods 0.000 claims abstract description 15
- IYLGZMTXKJYONK-ACLXAEORSA-N (12s,15r)-15-hydroxy-11,16-dioxo-15,20-dihydrosenecionan-12-yl acetate Chemical compound O1C(=O)[C@](CC)(O)C[C@@H](C)[C@](C)(OC(C)=O)C(=O)OCC2=CCN3[C@H]2[C@H]1CC3 IYLGZMTXKJYONK-ACLXAEORSA-N 0.000 claims abstract description 14
- IYLGZMTXKJYONK-UHFFFAOYSA-N ruwenine Natural products O1C(=O)C(CC)(O)CC(C)C(C)(OC(C)=O)C(=O)OCC2=CCN3C2C1CC3 IYLGZMTXKJYONK-UHFFFAOYSA-N 0.000 claims abstract description 14
- 239000000203 mixture Substances 0.000 claims abstract description 6
- 239000003344 environmental pollutant Substances 0.000 claims description 36
- 231100000719 pollutant Toxicity 0.000 claims description 36
- 238000001704 evaporation Methods 0.000 claims description 23
- 239000003673 groundwater Substances 0.000 claims description 23
- 230000008020 evaporation Effects 0.000 claims description 22
- 238000001556 precipitation Methods 0.000 claims description 16
- 238000013461 design Methods 0.000 claims description 14
- 238000010790 dilution Methods 0.000 claims description 14
- 239000012895 dilution Substances 0.000 claims description 14
- 238000005192 partition Methods 0.000 claims description 13
- 238000012423 maintenance Methods 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 9
- 241000251468 Actinopterygii Species 0.000 claims description 8
- 241000269799 Perca fluviatilis Species 0.000 claims description 6
- 238000004140 cleaning Methods 0.000 claims description 6
- 239000010865 sewage Substances 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 5
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 claims description 4
- 241000269800 Percidae Species 0.000 claims description 4
- 239000002689 soil Substances 0.000 claims description 4
- 238000011144 upstream manufacturing Methods 0.000 claims description 4
- 238000012544 monitoring process Methods 0.000 claims description 3
- 238000011160 research Methods 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 238000012937 correction Methods 0.000 claims description 2
- 230000008595 infiltration Effects 0.000 claims description 2
- 238000001764 infiltration Methods 0.000 claims description 2
- 230000010354 integration Effects 0.000 claims description 2
- 230000035515 penetration Effects 0.000 claims description 2
- 230000009467 reduction Effects 0.000 claims description 2
- 239000013589 supplement Substances 0.000 claims description 2
- 230000008901 benefit Effects 0.000 abstract description 4
- 230000008439 repair process Effects 0.000 abstract description 2
- 230000006870 function Effects 0.000 description 14
- 238000011161 development Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000004888 barrier function Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000004907 flux Effects 0.000 description 2
- 239000009537 qingyi Substances 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- XKMRRTOUMJRJIA-UHFFFAOYSA-N ammonia nh3 Chemical compound N.N XKMRRTOUMJRJIA-UHFFFAOYSA-N 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000008512 biological response Effects 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 241000894007 species Species 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
- 230000004083 survival effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000003911 water pollution Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
-
- 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
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Strategic Management (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明之目的就是提供一种人工干扰无水文资料地区河流环境流量分区界定计算方法,可有效解决人工干扰无水文资料地区河流环境流量的计算问题。以流域水生态功能分区为前提,界定环境流量的构成,然后采用水文比拟法和参数等值线法来推求无水文资料地区河流多年平均径流量及年内分配,以此为基础流量数据,对不同分区环境流量计算,确定河流不同节点的环境流量,本发明方法新颖独特,科学,是一种基于功能需求的环境流量计算方法,可以为受人工干扰程度较大且缺乏水文资料的河流进行水资源优化调配以及水利工程生态调度,保障河流功能正常发挥所需最小环境流量,修复河流生态系统,保护河流的生态环境,经济和社会效益巨大。
Description
技术领域
本发明涉及环保,特别是一种人工干扰无水文资料地区河流环境流量分区界定计算方法。
背景技术
环境流量的保障是维持河流生态系统健康发展的基础,对维持河流水生态系统的结构与功能正常发挥具有重要意义。随着社会经济的发展,人类生存发展与水资源有效供给之间的矛盾日益突出,河湖生态环境用水被挤占,河流环境流量难以满足,同时由于大量闸坝等水利工程的存在,严重改变了河流的天然流量过程,对河流生态系统造成严重影响。为保护和修复河流生态系统,开发适宜的环境流量计算方法十分必要。
目前,国内外环境流量的计算方法较多,国外比较前沿的有利用水资源规划方案“ribasim”的新功能计算高、低脉冲洪水下的环境流量,基于熵的生态水文系统动态相关的环境流量计算模型。国内环境流量计算方法主要有水文学法、水力学法、栖息地模拟法和整体法等,其中水文学法需要的历史资料及现场调查工作较少,多用于宏观规划和水生生物资料缺乏地区使用;与水文学法相比,水力学法考虑了生物栖息地要求,对数据的时间尺度要求不高,但其计算结果受河道断面形状较大;栖息地模拟法将生物响应与水力、水文状况直接联系,可以确定很多物种的最佳流量,但需要收集广泛的数据,成本较高;整体法明确着眼于广泛的环境效果,不限于单一的某种方法,需要科学领域内广泛的专家意见和技术,成本高且过程冗长。
尽管上述方法各有利弊,但多使用于受人工干扰程度小的自然河道环境流量的计算,同时是在具有水文站监测资料前提条件下开展的,侧重计算河流生态系统所需的环境流量。而我国河流现状即面临生态系统受损,又面临水质污染,在上述双重威胁存在的条件下,为保障河流经济社会功能的正常发挥,水库、闸坝等水利工程又大力兴建,河流人工干扰程度大,河流天然的水文情势丧失,维持河流生态系统功能正常发挥的环境流量难以保障,目前常用的计算方法也难以用于此类特殊流域的环境流量的计算。为解决上述问题,十分有必要在现有计算方法上进行改进,为受人工高度干扰及无水文资料地区的河流生态修复以及水利工程的科学调度提供一种适宜的环境流量计算方法。
发明内容
针对上述情况,为克服现有技术之缺陷,本发明之目的就是提供一种人工干扰无水文资料地区河流环境流量分区界定计算方法,可有效解决人工干扰无水文资料地区河流环境流量的计算问题。
本发明解决的技术方案是,以流域水生态功能分区为前提,界定环境流量的构成,然后采用水文比拟法和参数等值线法来推求无水文资料地区河流多年平均径流量及年内分配,以此为基础流量数据,对不同分区环境流量计算,确定河流不同节点的环境流量,包括以下步骤:
(1)确定流域范围,进行水生态功能分区:
进行水生态功能区二级分区,一级分区包括能够反映流域的气候、地势、土壤、植被自然环境要素指标的平均降水量、平均高程、地表粗糙度、土地利用强度;在一级分区基础上,进行二级分区,二级分区包括能够反映河段尺度上的栖息地环境以及人工干扰影响指标的小流域坡度、河道弯曲度、水网密度、水利工程密度和小流域污染负荷,利用GIS系统(地理信息系统软件)采用空间插值法及赋值来确定各参数值,用权重赋值分析方法确定分区;
(2)以水生态功能分区为基础,确定不同分区的环境流量组成及不同类别环境流量的计算方法:
水生态功能分区为城市人工高度干扰区和自然河道区,所述的城市人工高度干扰区包括维持稀释自净、景观娱乐、蒸发和地下水位的流量;所述的自然河道区包括维持水生生物栖息、蒸发和地下水位的流量。
维持稀释自净功能所需流量采用水质水量平衡模型,同时考虑水功能区纳污能力核算中的设计流量综合确定;景观娱乐所需流量用水深和水面面积方法来确定,景观娱乐用水水深为0.7~1.5m;维持水生生物栖息所需最小流量采用Tennant方法和湿周法确定,Tennant法计算时,每年10月~翌年3月用水期的流量采用30年以上的平均月流量的10%作为河道内的维持水生生物栖息最小环境流量,每年4~9月鱼类产卵育幼期采用30年以上的平均月流量的30%作为河道内维持鱼类产卵育幼期的最小环境流量;湿周法采用斜率求法来确定临界点;维持河流地下水位需水流量根据达西定律来确定;蒸发需水流量采用水面面积、水面蒸发深以及降雨量的参数进行计算;
(3)确定参证流域,计算无水文资料流域不同节点的基础流量数据,采用水文比拟法和参数等值线法来推求多年平均径流量及年内分配:
以参证流域近30年的天然流量数据为基础,利用GIS划分流域不同汇水区域,计算闭合汇水区域面积,以此为基础采用水文比拟法计算不同汇水区域节点流量;参数等值线法是依据多年平均径流深等值线图来计算年平均径流量,采用水文比拟法和参数等值线法计算结果平均值作为最终节点流量数据;
(4)建立河段环境流量计算模型及河流环境流量整合模型:
上游河段环境流量为非消耗性流量的最大值与消耗性流量之和,河流下游河段环境流量采用汇流式整合模型计算公式计算,分别是:
城市人工高度干扰区域河段计算公式:
自然河道区域河段计算公式:
其中:和分别为第n河段稀释自净、景观娱乐、水生生物栖息地、维持地下水位、水面蒸发和环境流量;为非消耗性需水流量,即稀释自净和景观娱乐所需流量的最大值;
同一水系下游河段汇流整合模型公式为:
式中:Q为下游河段的环境流量,Qifxh为上游所有河段的消耗性流量,Q11、Q21、Q31为上游河段非消耗性流量,通过上述数学模型计算,从而实现人工干扰无水文资料地区河流环境流量分区界定的计算。
本发明方法新颖独特,科学,是一种基于功能需求的环境流量计算方法,可以为受人工干扰程度较大且缺乏水文资料的河流进行水资源优化调配以及水利工程生态调度,保障河流功能正常发挥所需最小环境流量,修复河流生态系统,保护河流的生态环境,经济和社会效益巨大。
附图说明
图1为本发明的工艺流程图。
图2为本发明同一水系河流下游河段汇流整合模型图。
图3为本发明水质水量平衡构建模型图。
具体实施方式
以下结合具体情况对本发明的具体实施方式作详细说明。
由图1-3所示,本发明在具体实施中,包括以下步骤:
(1)确定流域范围,进行水生态功能分区:
进行水生态功能区二级分区,一级分区包括能够反映流域的气候、地势、土壤、植被自然环境要素指标的平均降水量、平均高程、地表粗糙度、土地利用强度;在一级分区基础上,进行二级分区,二级分区包括能够反映河段尺度上的栖息地环境以及人工干扰影响指标的小流域坡度、河道弯曲度、水网密度、水利工程密度和小流域污染负荷,利用GIS系统(地理信息系统软件)采用空间插值法及赋值来确定各参数值,用权重赋值分析方法确定分区;
(2)以水生态功能分区为基础,确定不同分区的环境流量组成及不同类别环境流量的计算方法:
控制河段的水生态功能分区为城市人工高度干扰区和自然河道区,所述的城市人工高度干扰区包括维持稀释自净、景观娱乐、蒸发和地下水位流量,其中维持稀释自净功能所需流量采用水质水量平衡模型,同时结合水功能区纳污能力核算中的设计流量综合确定,水质水量平衡模型构建时结合河段上游来水流量Q1、单位m3/s和污染物量P1、单位t/a;工业点源排入控制河段的排入量Q2、单位为m3/s和污染物排入量P2、单位为t/a;污水处理厂排入控制河段的排入量Q3、单位m3/s和污染物排入量P3、单位为t/a;为河段直排生活污水流量Q4、单位m3/s和污染物量P4、单位t/a;河段已有生态补水流量Q5、单位m3/s和污染物量P5、单位t/a;河段区间降雨径流的流量Q6、单位m3/s和污染物量P6、单位t/a;河段支流汇入的流量Q7、单位m3/s和污染物量P7、单位t/a;河段下游流量Q8、单位m3/s和污染物量P8、单位t/a;河段蒸发耗水流量Q耗、单位m3/s和污染物量P耗、单位t/a;河段自净削减污染物量Pzj、单位t/a;河段取水/渗漏平衡流量ΔQ、单位m3/s和污染物量ΔP、单位t/a(如图3所示)。
上述自净削减量采用综合衰减系数法来确定,其余采用实际监测数据。
区段污染物的综合削减系数ksj计算公式为:
式中:Q1,C1—分别为河段上游来水的流量、单位m3/s和污染物浓度、单位mg/L;
qi,ci—分别为排污口或支流的流量、单位m3/s和污染物浓度、单位mg/L;
Q8,C8—分别为河段下游的流量、单位m3/s和污染物浓度、单位mg/L;
ksj—污染物综合削减系数;
河流蒸发需水流量采用水面面积、水面蒸发深以及降雨量的参数,蒸发需水量采用如下公式计算:
式中:A为水面面积、单位m2;P为多年平均降雨量,单位mm/a;E为多年平均蒸发量,单位mm/a;
景观娱乐所需流量用水深和水面面积方法来确定,景观娱乐用水水深为0.7~1.5m。维持河流地下水位需水流量根据达西定律来确定,公式如下:
Q河补=10-4·Ks·I·A·L·t 式(6)
Q河补为单侧河道渗漏量,单位:万m3/a;Ks为剖面位置不同岩性的渗透系数,单位m/d;I为垂直于剖面的水力坡降;A为单位长度河道垂直于地下水流向的剖面面积,单位m2/m;L为河道或河段长度,单位m;t为河道或河段过水渗漏时间,单位d;
所述的自然河道区包括维持水生生物栖息、蒸发和地下水位流量,维持水生生物栖息所需最小流量采用Tennant方法和湿周法确定,Tennant法计算时,每年10月~翌年3月用水期的流量采用30年平均月流量的10%作为河道内的维持水生生物栖息最小环境流量,每年4~9月鱼类产卵育幼期采用30年平均月流量的30%作为河道内维持鱼类产卵育幼期的最小环境流量;湿周法采用斜率求法来确定临界点;
(3)采用水文比拟法和参数等值线法计算无水文资料流域的不同节点的基础流量数据:
水文比拟法是将参证流域的径流资料特征参数、分析按要求有选择地移置到设计流域上来的一种方法,参证流域应有较长的实测径流资料,其气候条件、下垫面条件情况与设计流域相近,具体公式如下:
Q设=k1k2Q参 式(7)
其中Q设、Q参分别为设计流域和参证流域的30年平均径流量,m3/s;k1、k2分别为流域面积和年降水量的修正系数;A设、A参分别为设计流域和参证流域的流域面积,km2;P设、P参分别为设计流域和参证流域的30年平均降水量,单位mm;
参数等值线法主要通过水文特征值的参数等值线图来进行计算,以1956~2000年平均径流深值线图来计算年平均径流量;
首先利用GIS系统找出各控制单元的形心,根据每个控制单元形心附近的两条等值线按比例内插计算各形心处的30年平均径流深,再乘以对应控制单元的流域面积,得到各控制单元的30年平均径流量:
Q设=mRA 式(9)
式中Q设为30年平均径流量,单位m3;R为30年平均径流深,单位mm;A为流域面积,单位km2;m为单位换算系数,m为1000;
为保证流量计算结果更为准确,采用两种方法计算结果的平均值作为无水文资料地区不同节点的30年平均径流量;
根据计算所得的30年平均径流量,采用水文比拟法来推求年径流量的年内分配,直接移用参证流域各种设计代表年的月径流量分配比,乘以研究流域的年径流量即得到年径流量的年内分配;
(4)上述流量都为河道内所需流量存在交叉重叠的现象,在进行环境流量整合时,不能简单的加和:
上游河段河流非消耗性需水流量取各个分项的最大值,河流环境流量等于河流非消耗性需水流量与消耗性需水流量之和,其中河流消耗性需水为蒸发及维持地下水位所需流量:
不同分区上游河段流量整合模型分别为:
城市人工高度干扰区:
自然河道区域河流环境流量模型为:
式中:和分别为第n河段稀释自净、景观娱乐、水生生物栖息地、维持地下水位、水面蒸发和环境流量,为非消耗性需水流量,为稀释自净和景观娱乐所需流量的最大值;
根据上述不同分区不同河段计算环境流量结果,对其进行整合得到下游河段(即可代表整条河流)所需环境流量,以汇流式河流水系为例来说明具体整合计算模型:
式中:Q为下游河段的环境流量,Qifxh为上游所有河段的消耗性流量,Q11、Q21、Q31为上游河段非消耗性流量。
本发明经实地应用,效果非常好,具有很好的实地应用价值,有关试验情况如下:
按照上述方法对我国北方无水文资料流域某人工高度干扰河流进行环境流量计算。
该河流属淮河流域沙颍河水系,全长149km,流域面积2362km2。流域四季气候特征为:春季干旱多风沙,夏季炎热多降水,秋季凉爽日照长,冬季寒冷少雨雪,全年平均气温在14.3~14.6℃之间,年极端最高气温为44℃,年极端最低气温为-17.5℃,多年平均降水量为727mm,属于北温带大陆性季风气候,但降水量年内分配不均,雨量集中在夏秋季节,汛期降水量占全年降水量的80%,降水在年内分配不均决定了径流在年内的分配不均,同时由于流域内闸坝较多(共29座),由于闸坝的蓄水拦截严重影响了河道内的径流分配。水面蒸发量为937.5~1000mm,多年平均径流量深为75mm,陆面蒸发量为653mm,根据颍河化行和黄桥1953~2013年的多年天然径流量,计算所得该流域(许昌段)多年平均径流量为1.12亿m3。
本发明以该河流所处流域为研究对象,分为城市人工高度干扰区和自然河道区域,其中人工高度干扰区主要包括稀释自净、景观娱乐、蒸发、维持地下水位所需流量,自然河道区域包括栖息地、蒸发和维持地下水位所需流量。
(1)以其汇入前后黄桥和化行两个水文站点长系列(1956-2013)的水文资料为基础,采用水文比拟法和参数等值线法确定不同节点的流量数据,其中参数等值线法中多年平均径流深参数分别为:清潩河禄马桥段72.95mm,清潩河高村桥段74.09mm,清泥河南外环桥段74.55mm,小泥河大石桥段78.89mm,石梁河三张闸段79.17mm,小洪河地方铁路桥段73.13mm。计算所得不同节点的基础流量数据为:
表1 清潩河(许昌段)流域流量情况表 流量单位:m3/s
(2)城市人工高度干扰区涉及不同类型环境流量计算时的主要参数确定如下:维持稀释自净所需流量中综合衰减系数ksj取值为:化学需氧量(COD)的衰减系数为0.15/d,氨氮为0.18/d;维持景观娱乐需水时涉及河段清潩河禄马桥段水深取值为1~1.5m,清潩河高村桥段水深为1.5~3m,清泥河南外环桥段1.0~2.0m,依据上述主要参数计算所得非消耗需水流量结果见表2。
表2 城市人工高度干扰区域河段非消耗需水流量
按照对应计算方法计算上述河段的蒸发和维持地下水位所需流量,与非消耗流量加和所得三个河段的环境流量。1月至12月流量范围分别为:清潩河长葛段0.903~1.323m3/s、清潩河高村桥段1.041~3.251m3/s、清泥河南外环桥段0.196~0.746m3/s。
(3)自然河道区涉及不同类型环境流量计算时主要参数确定如下:其中湿周法计算时三个河段沿河流方向的高程差与相应的河流长度的比值为0.0004,河道糙率取0.025。
表3 自然河道区河段维持水生生物栖息地所需流量 单位m3/s
按照对应计算方法计算上述河段的蒸发和维持地下水位所需流量,与维持水生生物栖息地所需流量加和所得三个河段的环境流量,1月至12月流量范围分别为:石梁河三张闸河段0.176~0.751m3/s、小洪河地方铁路桥河段0.173~0.998m3/s、小泥河大石桥河段0.294~1.096m3/s。
本发明在方法经对其他河流控制河段进行实验和实测,均取得了与上述实验相同和相挖的结果,表明方法稳定可靠,具有实际的应用价值,可以为受人工干扰程度较大且缺乏水文资料的河流进行水资源优化调配以及水利工程生态调度,保障河流功能正常发挥所需最小环境流量,修复河流生态系统,保护河流生态环境,造福于人类做出创造性的贡献,经济和社会效益巨大。
Claims (2)
1.一种人工干扰无水文资料地区河流环境流量分区界定计算方法,其特征在于,以流域水生态功能分区为前提,界定环境流量的构成,然后采用水文比拟法和参数等值线法来推求无水文资料地区河流多年平均径流量及年内分配,以此为基础流量数据,对不同分区环境流量计算,确定河流不同节点的环境流量,包括以下步骤:
(1)确定流域范围,进行水生态功能分区:
进行水生态功能区二级分区,一级分区包括能够反映流域气候、地势、土壤、植被自然环境要素指标的平均降水量、平均高程、地表粗糙度、土地利用强度;在一级分区基础上,进行二级分区,二级分区包括能够反映河段尺度上的栖息地环境以及人工干扰影响指标的小流域坡度、河道弯曲度、水网密度、水利工程密度和小流域污染负荷,利用GIS系统(地理信息系统软件)采用空间插值法及赋值来确定各参数值,用权重赋值分析方法确定分区;
(2)以水生态功能分区为基础,确定不同分区的环境流量组成及不同类别环境流量的计算方法:
水生态功能分区为城市人工高度干扰区和自然河道区,所述的城市人工高度干扰区包括维持稀释自净、景观娱乐、蒸发和地下水位的流量;所述的自然河道区包括维持水生生物栖息、蒸发和地下水位的流量;
维持稀释自净功能所需流量采用水质水量平衡模型,同时考虑水功能区纳污能力核算中的设计流量综合确定;景观娱乐所需流量用水深和水面面积方法来确定,景观娱乐用水水深为0.7~1.5m;维持水生生物栖息所需最小流量采用Tennant方法和湿周法确定,Tennant法计算时,每年10月~翌年3月用水期的流量采用30年以上的平均月流量的10%作为河道内的维持水生生物栖息最小环境流量,每年4~9月鱼类产卵育幼期采用30年以上的平均月流量的30%作为河道内维持鱼类产卵育幼期的最小环境流量;湿周法采用斜率求法来确定临界点;维持河流地下水位需水流量根据达西定律来确定;蒸发需水流量采用水面面积、水面蒸发深以及降雨量的参数进行计算;
(3)确定参证流域,计算无水文资料流域不同节点的基础流量数据,采用水文比拟法和参数等值线法来推求多年平均径流量及年内分配:
以参证流域近30年的天然流量数据为基础,利用GIS划分流域不同汇水区域,计算闭合汇水区域面积,以此为基础采用水文比拟法计算不同汇水区域节点流量;参数等值线法是依据多年平均径流深等值线图来计算年平均径流量,采用水文比拟法和参数等值线法计算结果平均值作为最终节点流量数据;
(4)建立河段环境流量计算模型及河流环境流量整合模型:
上游河段环境流量为非消耗性流量的最大值与消耗性流量之和,河流下游河段环境流量采用汇流式整合模型计算公式计算,分别是:
城市人工高度干扰区域河段计算公式:
自然河道区域河段计算公式:
其中:和分别为第n河段稀释自净、景观娱乐、水生生物栖息地、维持地下水位、水面蒸发和环境流量;为非消耗性需水流量,即稀释自净和景观娱乐所需流量的最大值;
同一水系下游河段汇流整合模型公式为:
式中:Q为下游河段的环境流量,Qifxh为上游所有河段的消耗性流量,Q11、Q21、Q31为上游河段非消耗性流量,通过上述数学模型计算,从而实现人工干扰无水文资料地区河流环境流量分区界定的计算。
2.根据权利要求1所述的人工干扰无水文资料地区河流环境流量分区界定计算方法,其特征在于,包括以下步骤:
(1)确定流域范围,进行水生态功能分区:
进行水生态功能区二级分区,一级分区包括能够反映流域的气候、地势、土壤、植被自然环境要素指标的平均降水量、平均高程、地表粗糙度、土地利用强度;在一级分区基础上,进行二级分区,二级分区包括能够反映河段尺度上的栖息地环境以及人工干扰影响指标的小流域坡度、河道弯曲度、水网密度、水利工程密度和小流域污染负荷,利用GIS系统采用空间插值法及赋值确定一级分区、二级分区的指标,用权重赋值分析方法确定分区;
(2)以水生态功能分区为基础,确定不同分区的环境流量构成及不同类别环境流量的计算方法:
控制河段的水生态功能分区为城市人工高度干扰区和自然河道区,所述的城市人工高度干扰区包括维持稀释自净、景观娱乐、蒸发和地下水位流量,其中维持稀释自净功能所需流量采用水质水量平衡模型,同时结合水功能区纳污能力核算中的设计流量综合确定,水质水量平衡模型构建时结合河段上游来水流量Q1、单位m3/s和污染物量P1、单位t/a;工业点源排入控制河段的排入量Q2、单位为m3/s和污染物排入量P2、单位为t/a;污水处理厂排入控制河段的排入量Q3、单位m3/s和污染物排入量P3、单位为t/a;为河段直排生活污水流量Q4、单位m3/s和污染物量P4、单位t/a;河段已有生态补水流量Q5、单位m3/s和污染物量P5、单位t/a;河段区间降雨径流的流量Q6、单位m3/s和污染物量P6、单位t/a;河段支流汇入的流量Q7、单位m3/s和污染物量P7、单位t/a;河段下游流量Q8、单位m3/s和污染物量P8、单位t/a;河段蒸发耗水流量Q耗、单位m3/s和污染物量P耗、单位t/a;河段自净削减污染物量Pzj、单位t/a;河段取水/渗漏平衡流量ΔQ、单位m3/s和污染物量ΔP、单位t/a。
上述自净削减量采用综合衰减系数法来确定,其余采用实际监测数据。
区段污染物的综合削减系数ksj计算公式为:
式中:Q1,C1—分别为河段上游来水的流量、单位m3/s和污染物浓度、单位mg/L;
qi,ci—分别为排污口或支流的流量、单位m3/s和污染物浓度、单位mg/L;
Q8,C8—分别为河段下游的流量、单位m3/s和污染物浓度、单位mg/L;
ksj—污染物综合削减系数;
河流蒸发需水流量采用水面面积、水面蒸发深以及降雨量的参数,蒸发需水量采用如下公式计算:
式中:A为水面面积、单位m2;P为多年平均降雨量,单位mm/a;E为多年平均蒸发量,单位mm/a;
景观娱乐所需流量用水深和水面面积方法来确定,景观娱乐用水水深为0.7~1.5m。维持河流地下水位需水流量根据达西定律来确定,公式如下:
Q河补=10-4·Ks·I·A·L·t 式(6)
Q河补为单侧河道渗漏量,单位:万m3/a;Ks为剖面位置不同岩性的渗透系数,单位m/d;I为垂直于剖面的水力坡降;A为单位长度河道垂直于地下水流向的剖面面积,单位m2/m;L为河道或河段长度,单位m;t为河道或河段过水渗漏时间,单位d;
所述的自然河道区包括维持水生生物栖息、蒸发和地下水位流量,维持水生生物栖息所需最小流量采用Tennant方法和湿周法确定,Tennant法计算时,每年10月~翌年3月用水期的流量采用30年平均月流量的10%作为河道内的维持水生生物栖息最小环境流量,每年4~9月鱼类产卵育幼期采用30年平均月流量的30%作为河道内维持鱼类产卵育幼期的最小环境流量;湿周法采用斜率求法来确定临界点;
(3)采用水文比拟法和参数等值线法计算无水文资料流域的不同节点的基础流量数据:
水文比拟法是将参证流域的径流资料特征参数、分析按要求有选择地移置到设计流域上来的一种方法,公式如下:
Q设=k1k2Q参 式(7)
其中Q设、Q参分别为设计流域和参证流域的30年平均径流量,m3/s;k1、k2分别为流域面积和年降水量的修正系数;A设、A参分别为设计流域和参证流域的流域面积,km2;P设、P参分别为设计流域和参证流域的30年平均降水量,单位mm;
参数等值线法通过水文特征值的参数等值线图来进行计算,以30年平均径流深值线图来计算年平均径流量:
首先利用GIS系统找出各控制单元的形心,根据每个控制单元形心附近的两条等值线按比例内插计算各形心处的多年平均径流深,再乘以对应控制单元的流域面积,得到各控制单元的30年平均径流量:
Q设=mRA 式(9)
式中Q设为30年平均径流量,单位m3;R为30年平均径流深,单位mm;A为流域面积,单位km2;m为单位换算系数,m为1000;
为保证流量计算结果更为准确,采用两种方法计算结果的平均值作为无水文资料地区不同节点的30年平均径流量;
根据计算所得的30年平均径流量,采用水文比拟法来推求年径流量的年内分配,直接移用参证流域各种设计代表年的月径流量分配比,乘以研究流域的年径流量即得到年径流量的年内分配;
(4)上述流量都为河道内所需流量存在交叉重叠的现象,在进行环境流量整合时,不能简单的加和:
上游河段河流非消耗性需水流量取各个分项的最大值,河流环境流量等于河流非消耗性需水流量与消耗性需水流量之和,其中河流消耗性需水为蒸发及维持地下水位所需流量:
不同分区上游河段流量根据式(1)计算,自然河道区域河流环境流量根据式(2)计算,根据计算结果,由式(3)进行计算,得到下游河段所需环境流量,从而实现人工干扰无水文资料地区河流环境流量分区界定的计算。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611139996.XA CN106777959B (zh) | 2016-12-12 | 2016-12-12 | 人工干扰无水文资料地区河流环境流量分区界定计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611139996.XA CN106777959B (zh) | 2016-12-12 | 2016-12-12 | 人工干扰无水文资料地区河流环境流量分区界定计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106777959A true CN106777959A (zh) | 2017-05-31 |
CN106777959B CN106777959B (zh) | 2019-02-01 |
Family
ID=58880169
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611139996.XA Active CN106777959B (zh) | 2016-12-12 | 2016-12-12 | 人工干扰无水文资料地区河流环境流量分区界定计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106777959B (zh) |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107341354A (zh) * | 2017-07-11 | 2017-11-10 | 辽宁省环境科学研究院 | 一种计算整个河段生态需水量的方法 |
CN107480884A (zh) * | 2017-08-14 | 2017-12-15 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种多闸坝平原河流冬季生态流量推求方法 |
CN107563647A (zh) * | 2017-09-05 | 2018-01-09 | 郑州大学 | 强干扰区域基于河流主导生态环境功能分区的环境流量界定方法 |
CN108287950A (zh) * | 2017-12-27 | 2018-07-17 | 环境保护部环境规划院 | 基于控制单元水环境质量目标管理的水质模拟方法 |
CN108305022A (zh) * | 2018-04-15 | 2018-07-20 | 李善综 | 一种基于水质约束的中小河流水资源可利用量计算方法 |
CN108664453A (zh) * | 2018-05-15 | 2018-10-16 | 长江水利委员会水文局 | 一种缺乏水文资料河流悬移质输沙量的计算方法 |
CN109325206A (zh) * | 2018-09-10 | 2019-02-12 | 柳创新 | 一种降雨径流模型参数优化方法 |
CN109472102A (zh) * | 2018-11-21 | 2019-03-15 | 北京师范大学 | 一种基于蚂蚁河流域的桥涵构造推算方法 |
CN110163423A (zh) * | 2019-04-24 | 2019-08-23 | 中国科学院新疆生态与地理研究所 | 一种河流生态恢复系统 |
CN110175948A (zh) * | 2019-05-24 | 2019-08-27 | 郑州大学 | 一种基于河流整体健康的生态环境需水阈值量化方法 |
CN110400049A (zh) * | 2019-06-24 | 2019-11-01 | 水利部农村电气化研究所 | 农村水电站的生态流量的核定方法及核定装置 |
CN110598936A (zh) * | 2019-09-17 | 2019-12-20 | 中国水利水电科学研究院 | 一种河流生态需水计算系统及信息传输方法 |
CN110782068A (zh) * | 2019-09-23 | 2020-02-11 | 中国水利水电科学研究院 | 生态流量保障方法、装置、电子设备及计算机存储介质 |
CN111581586A (zh) * | 2020-04-28 | 2020-08-25 | 生态环境部卫星环境应用中心 | 基于套合模型的湖库水质各向异性插值方法和装置 |
CN111651855A (zh) * | 2020-04-23 | 2020-09-11 | 中国地质大学(武汉) | 基于物理模型的河网区水质断面污染贡献率核算方法 |
CN113240171A (zh) * | 2021-05-10 | 2021-08-10 | 福建省水利水电勘测设计研究院 | 河口叶脉状河网干流最优引水流量确定方法 |
WO2021238567A1 (zh) * | 2020-05-27 | 2021-12-02 | 中国长江三峡集团有限公司 | 一种考虑区间支流来水的干流水库生态调度方法 |
CN116681214A (zh) * | 2023-07-26 | 2023-09-01 | 四川省生态环境科学研究院 | 一种受干扰涉水型河流水生态参照点的选择及构建方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101645113A (zh) * | 2009-09-03 | 2010-02-10 | 北京师范大学 | 多目标整合河口生态需水计算方法 |
CN101865687A (zh) * | 2009-07-31 | 2010-10-20 | 北京师范大学 | 流域生态需水整合分析技术 |
CN104794584A (zh) * | 2015-04-30 | 2015-07-22 | 中国科学院南京地理与湖泊研究所 | 一种湖泊型流域水生态功能一二级分区方法 |
CN105022932A (zh) * | 2015-08-14 | 2015-11-04 | 北京师范大学 | 一种基于水生态功能的流域分区方法 |
-
2016
- 2016-12-12 CN CN201611139996.XA patent/CN106777959B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101865687A (zh) * | 2009-07-31 | 2010-10-20 | 北京师范大学 | 流域生态需水整合分析技术 |
CN101645113A (zh) * | 2009-09-03 | 2010-02-10 | 北京师范大学 | 多目标整合河口生态需水计算方法 |
CN104794584A (zh) * | 2015-04-30 | 2015-07-22 | 中国科学院南京地理与湖泊研究所 | 一种湖泊型流域水生态功能一二级分区方法 |
CN105022932A (zh) * | 2015-08-14 | 2015-11-04 | 北京师范大学 | 一种基于水生态功能的流域分区方法 |
Non-Patent Citations (4)
Title |
---|
R.E.THARME: "A global perspective on environmental flow assessment: emerging trends in the development and application of environmental flow methodologies for rivers", 《RIVER RESEARCH AND APPLICATIONS》 * |
刘新华 等: "塔里木河下游生态需水估算", 《中国沙漠》 * |
周振民 等: "河道生态需水量计算方法及应用研究", 《中国农村水利水电》 * |
李兴春: "松嫩平原湿地生态环境需水量研究", 《中国博士学位论文全文数据库 工程科技Ⅰ辑》 * |
Cited By (33)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107341354A (zh) * | 2017-07-11 | 2017-11-10 | 辽宁省环境科学研究院 | 一种计算整个河段生态需水量的方法 |
CN107480884B (zh) * | 2017-08-14 | 2019-10-22 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种多闸坝平原河流冬季生态流量推求方法 |
CN107480884A (zh) * | 2017-08-14 | 2017-12-15 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种多闸坝平原河流冬季生态流量推求方法 |
CN107563647A (zh) * | 2017-09-05 | 2018-01-09 | 郑州大学 | 强干扰区域基于河流主导生态环境功能分区的环境流量界定方法 |
CN107563647B (zh) * | 2017-09-05 | 2020-08-04 | 郑州大学 | 强干扰区域基于河流主导生态环境功能分区的环境流量界定方法 |
CN108287950A (zh) * | 2017-12-27 | 2018-07-17 | 环境保护部环境规划院 | 基于控制单元水环境质量目标管理的水质模拟方法 |
CN108287950B (zh) * | 2017-12-27 | 2021-09-14 | 生态环境部环境规划院 | 基于控制单元水环境质量目标管理的水质模拟方法 |
CN108305022A (zh) * | 2018-04-15 | 2018-07-20 | 李善综 | 一种基于水质约束的中小河流水资源可利用量计算方法 |
CN108305022B (zh) * | 2018-04-15 | 2021-09-21 | 水利部珠江水利委员会技术咨询(广州)有限公司 | 一种基于水质约束的中小河流水资源可利用量计算方法 |
CN108664453A (zh) * | 2018-05-15 | 2018-10-16 | 长江水利委员会水文局 | 一种缺乏水文资料河流悬移质输沙量的计算方法 |
CN108664453B (zh) * | 2018-05-15 | 2021-11-16 | 长江水利委员会水文局 | 一种缺乏水文资料河流悬移质输沙量的计算方法 |
CN109325206B (zh) * | 2018-09-10 | 2023-03-24 | 柳创新 | 一种降雨径流模型参数优化方法 |
CN109325206A (zh) * | 2018-09-10 | 2019-02-12 | 柳创新 | 一种降雨径流模型参数优化方法 |
CN109472102A (zh) * | 2018-11-21 | 2019-03-15 | 北京师范大学 | 一种基于蚂蚁河流域的桥涵构造推算方法 |
CN109472102B (zh) * | 2018-11-21 | 2023-02-28 | 北京师范大学 | 一种基于蚂蚁河流域的桥涵构造推算方法 |
CN110163423A (zh) * | 2019-04-24 | 2019-08-23 | 中国科学院新疆生态与地理研究所 | 一种河流生态恢复系统 |
CN110175948A (zh) * | 2019-05-24 | 2019-08-27 | 郑州大学 | 一种基于河流整体健康的生态环境需水阈值量化方法 |
CN110400049A (zh) * | 2019-06-24 | 2019-11-01 | 水利部农村电气化研究所 | 农村水电站的生态流量的核定方法及核定装置 |
CN110400049B (zh) * | 2019-06-24 | 2023-01-20 | 水利部农村电气化研究所 | 农村水电站的生态流量的核定方法及核定装置 |
CN110598936A (zh) * | 2019-09-17 | 2019-12-20 | 中国水利水电科学研究院 | 一种河流生态需水计算系统及信息传输方法 |
CN110598936B (zh) * | 2019-09-17 | 2022-07-22 | 中国水利水电科学研究院 | 一种河流生态需水计算系统及信息传输方法 |
CN110782068A (zh) * | 2019-09-23 | 2020-02-11 | 中国水利水电科学研究院 | 生态流量保障方法、装置、电子设备及计算机存储介质 |
CN111651855A (zh) * | 2020-04-23 | 2020-09-11 | 中国地质大学(武汉) | 基于物理模型的河网区水质断面污染贡献率核算方法 |
CN111651855B (zh) * | 2020-04-23 | 2022-11-18 | 中国地质大学(武汉) | 基于物理模型的河网区水质断面污染贡献率核算方法 |
CN111581586B (zh) * | 2020-04-28 | 2021-02-05 | 生态环境部卫星环境应用中心 | 基于套合模型的湖库水质各向异性插值方法和装置 |
CN111581586A (zh) * | 2020-04-28 | 2020-08-25 | 生态环境部卫星环境应用中心 | 基于套合模型的湖库水质各向异性插值方法和装置 |
GB2600368A (en) * | 2020-05-27 | 2022-04-27 | China Three Gorges Corp | Method for performing ecological dispatching on trunk stream reservoir by taking inflow water from interregional tributary into consideration |
WO2021238567A1 (zh) * | 2020-05-27 | 2021-12-02 | 中国长江三峡集团有限公司 | 一种考虑区间支流来水的干流水库生态调度方法 |
GB2600368B (en) * | 2020-05-27 | 2022-10-19 | China Three Gorges Corp | Main stream reservoir ecological modulation method considering incoming water from interval tributaries |
CN113240171B (zh) * | 2021-05-10 | 2022-06-03 | 福建省水利水电勘测设计研究院有限公司 | 河口叶脉状河网干流最优引水流量确定方法 |
CN113240171A (zh) * | 2021-05-10 | 2021-08-10 | 福建省水利水电勘测设计研究院 | 河口叶脉状河网干流最优引水流量确定方法 |
CN116681214A (zh) * | 2023-07-26 | 2023-09-01 | 四川省生态环境科学研究院 | 一种受干扰涉水型河流水生态参照点的选择及构建方法 |
CN116681214B (zh) * | 2023-07-26 | 2023-10-31 | 四川省生态环境科学研究院 | 一种受干扰涉水型河流水生态参照点的选择及构建方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106777959B (zh) | 2019-02-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106777959B (zh) | 人工干扰无水文资料地区河流环境流量分区界定计算方法 | |
Xu et al. | The hydrological calibration and validation of a complexly-linked watershed–reservoir model for the Occoquan watershed, Virginia | |
CN107885958A (zh) | 一种平原感潮河网区纳污能力计算方法 | |
Şen | Reservoirs for water supply under climate change impact—a review | |
CN111062125B (zh) | 海绵型综合管廊水文效应评估方法 | |
Baena-Ruiz et al. | Summarizing the impacts of future potential global change scenarios on seawater intrusion at the aquifer scale | |
Gu et al. | Modeling effects of low impact development on road salt transport at watershed scale | |
Gotkowitz et al. | Groundwater flood of a river terrace in southwest Wisconsin, USA | |
Zhang et al. | Quantitative analysis of self-purification capacity of non-point source pollutants in watersheds based on SWAT model | |
Lyu et al. | Towards sustainable water regulation based on a distributed hydrological model for a heavily polluted urban river, northwest China | |
CN109033589B (zh) | 一种基于swmm模型考虑lid净化作用的污染物去除率计算方法 | |
CN110728062A (zh) | 一种基于swmm的农村非点源污染模拟方法 | |
CN112942230B (zh) | 面向水污染控制的城市人工湖水域面积值域计算方法 | |
Zhang et al. | Restoration of urban water environment landscape system and SWMM technology integration in water-scarce cities | |
CN109472102B (zh) | 一种基于蚂蚁河流域的桥涵构造推算方法 | |
Safriani et al. | Evaluation of Krueng Tripa River Capacity in Ujung Krueng Village, Nagan Raya Regency, Indonesia | |
Fang et al. | Zero increase in peak discharge for sustainable development | |
Shi et al. | Research of Sponge City Landscape Planning Based on Landscape Security Pattern | |
Greenberg | Urban Hydrological Modeling of the Malden River Using the Storm Water Management Model (SWMM) | |
Ladson | Using WSUD to Restore Predevelopment Hydrology | |
Berris | Daily flow-routing simulations for the Truckee River, California and Nevada | |
Doulgeris et al. | Use of modern technology for the protection and management of water resources in Strymonas/Struma River basin | |
Li et al. | POLLUTION LOAD AND ECOLOGICAL REPLENISHMENT PLAN OF LIJIANG RIVER, CHINA. | |
Jafari | Modeling of flood and retention in the Sogn allotment garden | |
Lim et al. | Revitalization of an urban stream in terms of water quality and recreational usability: Assessment of options through hydrologic and hydraulic modeling |
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 |