CN111079085A - 一种考虑提升量的生态流量确定方法 - Google Patents
一种考虑提升量的生态流量确定方法 Download PDFInfo
- Publication number
- CN111079085A CN111079085A CN201911236748.0A CN201911236748A CN111079085A CN 111079085 A CN111079085 A CN 111079085A CN 201911236748 A CN201911236748 A CN 201911236748A CN 111079085 A CN111079085 A CN 111079085A
- Authority
- CN
- China
- Prior art keywords
- ecological
- flow
- water
- river
- annual
- 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
- 238000000034 method Methods 0.000 title claims abstract description 82
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 173
- 239000004576 sand Substances 0.000 claims abstract description 61
- 238000000746 purification Methods 0.000 claims abstract description 17
- 238000012795 verification Methods 0.000 claims abstract description 6
- 230000007613 environmental effect Effects 0.000 claims abstract description 5
- 238000004364 calculation method Methods 0.000 claims description 15
- 238000007865 diluting Methods 0.000 claims description 13
- 230000008859 change Effects 0.000 claims description 10
- 239000003344 environmental pollutant Substances 0.000 claims description 8
- 231100000719 pollutant Toxicity 0.000 claims description 8
- 238000009826 distribution Methods 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 7
- 239000000356 contaminant Substances 0.000 claims description 6
- 230000004941 influx Effects 0.000 claims description 6
- 238000003860 storage Methods 0.000 claims description 6
- 239000010865 sewage Substances 0.000 claims description 5
- 241000251468 Actinopterygii Species 0.000 claims description 4
- 238000004140 cleaning Methods 0.000 claims description 3
- 230000002708 enhancing effect Effects 0.000 claims description 3
- 238000011010 flushing procedure Methods 0.000 claims description 3
- 238000009991 scouring Methods 0.000 claims description 3
- 239000013049 sediment Substances 0.000 claims description 3
- 230000007480 spreading Effects 0.000 claims description 3
- 238000003892 spreading Methods 0.000 claims description 3
- 239000013589 supplement Substances 0.000 claims description 3
- 238000007726 management method Methods 0.000 description 8
- 230000001965 increasing effect Effects 0.000 description 4
- 238000013468 resource allocation Methods 0.000 description 3
- 241000282414 Homo sapiens Species 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000011144 upstream manufacturing Methods 0.000 description 2
- 244000025254 Cannabis sativa Species 0.000 description 1
- 241000196324 Embryophyta Species 0.000 description 1
- 241001465754 Metazoa Species 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009395 breeding Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- 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
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0637—Strategic management or analysis, e.g. setting a goal or target of an organisation; Planning actions based on goals; Analysis or evaluation of effectiveness of goals
-
- E—FIXED CONSTRUCTIONS
- E02—HYDRAULIC ENGINEERING; FOUNDATIONS; SOIL SHIFTING
- E02B—HYDRAULIC ENGINEERING
- E02B1/00—Equipment or apparatus for, or methods of, general hydraulic engineering, e.g. protection of constructions against ice-strains
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C13/00—Surveying specially adapted to open water, e.g. sea, lake, river or canal
- G01C13/002—Measuring the movement of open water
-
- 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/08—Construction
-
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2219/00—Indexing scheme relating to application aspects of data processing equipment or methods
- G06F2219/10—Environmental application, e.g. waste reduction, pollution control, compliance with environmental legislation
-
- 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
- Y02A20/00—Water conservation; Efficient water supply; Efficient water use
- Y02A20/152—Water filtration
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Economics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Operations Research (AREA)
- General Business, Economics & Management (AREA)
- Educational Administration (AREA)
- Marketing (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Software Systems (AREA)
- Evolutionary Biology (AREA)
- Probability & Statistics with Applications (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Entrepreneurship & Innovation (AREA)
- Development Economics (AREA)
- Civil Engineering (AREA)
- Mechanical Engineering (AREA)
- Structural Engineering (AREA)
- Hydrology & Water Resources (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Game Theory and Decision Science (AREA)
Abstract
本发明提供了一种考虑提升量的生态流量确定方法,属于环境工程领域,包括:收集计算生态流量所需的相关数据;确定生态基流;选取生态基流的上限Qmax和下限Qmin,初步确定生态基流的区间[Qmin,Qmax];验证生态基流Q的下限Qmin;计算三个服务对象景观湿地需水B1、输沙需水B2和稀释自净需水B3;对比三个服务对象的需水量,确定提升量S;找出最小值和最大值确定区间的提升量下限Smin和提升量上限Smax;将生态基流Q和提升量S结合起来综合确定生态流量。该方法将复杂问题简单化,同时提出将生态基流作为考核目标,更直接有效;此外,对生态流量的区间化管理也能体现其自然属性和功能属性,同时强调人服务于生态的服务属性。
Description
技术领域
本发明属于环境工程领域,具体涉及一种考虑提升量的生态流量确定方法。
背景技术
河湖本身不存在生态流量问题,因为水生动植物本来就是根据河湖生境而自然生长。但长期以来,由于人类在河湖内外用水需求日益增加,不断占用自然生态系统依赖的河湖水量,基流不够用,从而出现了所谓的生态流量问题。
与生态流量相关的两个概念是生态基流和生态需水,生态基流是为维持河流基本形态和基本生态功能所需的最小流量,功能范围是自然生态系统,侧重点是河道不断流。由于径流时空分布不均,需水总量并不能保证在单位时间生态基流满足河流基本生态功能,因此需要在需水过程上对基流提出要求,对于生态基流的界定,社会各界基本上达成一致。目前水资源短缺、水生态损害等成为常态问题,在此基础上提出用水的优先级是生活、生态、生产,因此生态需水成为广泛讨论的对象,与之联系的就是生态流量,关于这两个名词大家说法不一,有一部分人认为生态需水是将生态系统结构、功能和生态过程维持在一定水平所需要的水量,指一定生态保护目标对应的水生态系统对于水量的需求,而生态流量就是生态需水中的某个流量,也有一部分人认为生态流量是生态基流和敏感期生态流量的组合,对于这几个概念的认识影响了具体的做法,最主要的是生态流量难以确定,给管理造成了混乱。
其中涉及到的技术方法是关于生态基流、生态流量的确定,生态基流是基于算法导向的,确定的方法有四大类,分别是水文学方法、水力学方法、栖息地法和整体法,生态流量的确定要根据具体对象的需水情况来定。
目前关于生态流量的确定争议很多,存在概念不清的情况,没有说清与生态基流、生态需水的关系;由于生态服务对象众多,考虑不同时空分布下的生态需水难以量化;生态流量是与生态需水紧密相关的,而对应的算法太多太乱,一是理解有误,算出的结果不是生态流量,二是算法所需的资料无法获得、参数难以确定,大多数计算方法在考虑众多服务对象之后得到的是一个值或一组值,无法适应时空变化、需水变化、来水变化,且难以作为考核的依据,给管理造成不便。
由于不同时空分布下生态服务对象的多样性,生态需水难以量化,且相关部门现阶段是根据断面的实测流量来考核的,目的也是保证河道不断流,但是所定断面的流量考核标准过高,考核失去意义。
因此,本申请提出一种考虑提升量的生态流量确定方法。
发明内容
为了克服上述现有技术存在的不足,本发明提供了一种考虑提升量的生态流量确定方法。
为了实现上述目的,本发明提供如下技术方案:
一种考虑提升量的生态流量确定方法,包括以下步骤:
步骤1:收集计算生态流量所需的相关数据;
步骤2:从收集的相关数据中确定生态基流,确定方法包括改进蒙大拿法、年内展布法和最小月平均流量法;
步骤3:对所述改进蒙大拿法、年内展布法和最小月平均流量法三种方法分别确定的生态基流M1、M2和M3进行比较分析,选取最大值和最小值分别作为生态基流的上限Qmax和下限Qmin,初步确定生态基流的区间[Qmin,Qmax],将区间的下限Qmin初步确定为生态基流,见式(1)、(2);
Qmin=min{M1,M3,M3} (1)
Qmax=max{M1,M3,M3} (2)
步骤4:对所述步骤3的生态基流的下限Qmin进行验证;
若初步确定的下限Qmin大于规定值,则下限Qmin为生态基流的下限值;否则将规定值作为生态基流的下限值;
步骤5:在所述生态基流的基础上计算河段3个服务对象需水量,3个服务包括景观湿地需水Wl、输沙需水Ws和稀释自净需水Wc;
步骤6:对比分析所述步骤5计算的三个服务对象的需水量,确定提升量S;提升量是一个区间,因此找出最小值和最大值确定区间的提升量下限Smin和提升量上限Smax,为体现年内变化,提升量区间按月份分别确定;
Smin=min{Wl,Ws,Wc} (3)
Smax=max{Wl,Ws,Wc} (4)
步骤7:将所述生态基流和提升量S结合起来综合确定生态流量,将步骤4 验证过后的生态基流作为考核的最低标准,也就是提升量的下限。
优选地,所述步骤1中收集的相关数据包括河段含沙量、输沙量、换水周期、多年平均蓄水量、污染物浓度、目标水质浓度和长系列历史流量数据。
优选地,所述步骤2中,改进蒙大拿法是在蒙大拿法的基础上考虑典型年的选取及分期百分比的确定,典型年的选取根据计算不同季节模数最接近1.0 的年份,分期是将一年分为一般用水期和鱼类产卵育幼期,分别取典型年平均流量的10%和30%作为生态基流;
年内展布法是用多年月均径流量乘以同期均值比得到一年的生态基流过程,同期均值比是根据多年年均径流量与年内各月最小月均径流量的年均值这两个典型的水文特征变量进行确定与量化的,具体见式(5)-(10);
qmin(i)=min(qij) (8)
式中:为第i个月的多年月均径流量,为多年年均径流量,qij为第j年第i个月的月均径流量,qmin(i)为第i个月的多年最小月均径流量,为最小年均径流量,n为统计年数,η为同期均值比,Qi为各月的生态基流;
最小月平均流量法是将最小月平均流量的多年平均值作为河道的生态基流,具体见式(11);
式中:Qij为第i个月第j天的平均流量,Qm为第m个月的生态基流,n为统计年数。
优选地,所述步骤5中景观湿地需水Wl、输沙需水Ws和稀释自净需水Wc的具体计算方法如下:
(1)景观湿地需水:指为维持景观或者湿地生态和环境功能所要消耗和补充的水量,采用换水周期法进行计算,具体见式(12);
式中:T为换水周期,w”为多年平均蓄水量,Q为景观湿地需水流量;
(2)输沙需水:河流输沙需水量是指为了维持河床冲刷与淤积的动态平衡,必须在河道内保持一定的水量,该水量与水流含沙量直接相关,且随着产沙量、含沙量等的年内分配而动态变化,计算公式见式(13)-(15);
w'=ηα*ww (13)
ww=w-ws/γs (14)
η=s进-s出 (15)
式中:W′为输沙需水量,η′为输沙效率,α为指数,Ww为净水量;W为径流量,Ws为输沙量,γs为泥沙容重,通常取2.65t/m3;S进为进入河段的含沙量,S出为流出河段的含沙量;当η′<1时,河段处于冲刷状态,进口站含沙量小于出口站含沙量,此时,α=1;当η′>=1时,河段淤积或冲淤平衡,进口站含沙量大于或等于出口站含沙量,此时,α=0;
(3)稀释自净需水:指在加强污染源防治的前提下,使河流水体达到水功能区划确定的水质目标所要求的最小水量;采用混合模型与衰减模型对河段稀释自净需水量进行模拟;在河段i处,上断面为Ai-1,下断面为Ai,河段长为Li,假设在河段内共有多个排污口及支流的污染物汇入,为便于计算,记做断面B,具体见式(16);
式中:Woi为稀释自净需水量,Ci-1为Ai-1处的污染物浓度,si为B点处汇入的污染物浓度,qi为B点处的汇入流量,Csi为Ai段面处的目标水质浓度,Li为河段长、ui为河段i处的流速,xi为断面B到断面Ai-1的距离,Qi-1为Ai-1处的流量。。
本发明提供的考虑提升量的生态流量确定方法具有以下有益效果:
(1)从实际应用的角度,从概念出发给出了生态流量的确定办法,即生态基流加上提升量,相比于具体服务对象的生态需水难以量化的特点,本发明将复杂问题简单化,同时提出将生态基流作为考核目标,更直接有效;此外,对生态流量的区间化管理也能体现其自然属性和功能属性,同时强调人服务于生态的服务属性;
(2)解决了生态流量难以量化的困难,避免了因生态流量算法众多、服务对象复杂所带来的管理混乱,并且给出的是一个区间,适应性更强,同时考虑到考核的需求,用生态流量作为考核达标与否的判断,更容易被接受,使考核变的简单有效,清晰可管理。
附图说明
图1为本发明实施例1的考虑提升量的生态流量确定方法的流程图。
具体实施方式
为了使本领域技术人员更好的理解本发明的技术方案并能予以实施,下面结合附图和具体实施例对本发明进行详细说明。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。
实施例1
本实施例在辨析生态基流、生态需水、生态流量三者概念的基础上给出生态流量的定义,由定义出发,说明生态流量的确定方法;生态基流前面已作出说明,需要强调的是生态基流最重要的一个特点就是要保证不断流,现在阐述生态需水与生态流量的关系,讲生态需水,首先要搞清楚何为生态对象,如何去服务,之后回过头看,生态的需水量是多少,要用多少的生态流量来保证对应上,故不讲生态需水,生态流量就没有依据。
生态需水可分为河道内和河道外两种,具体的服务对象包括景观、鱼类、水生生物、河道输沙、自净等,功能范围是自然生态系统、社会生态系统,需水条件为需水量、需水过程。不同河道在不同时空尺度下的生态服务对象都有所差别,生态需水不仅与自然状态下的河道变化有关,也受人类社会活动的影响,一些功能性需水例如输沙和自净等,既要求需水过程,也需要从需水量上来保证。
在上述认识的基础上,提出生态流量的概念,第一步在生态基流的基础上,使生态恢复到原生态几乎不可能,只能改善和修复,增加一个提升量,提升量可大可小,使用区间化的思想,给出生态流量的界限;第二步一定要规制化,建立规制库,如:不同河段不同季节用哪种规制,不同服务对象用哪种规制;考虑时空分布下水量如何调配,例如:丰增枯减,在枯水期,有好的手段、技术、较强的调节能力,能在基流基础上满足提升量。
下面结合附图和具体实现方式对本发明进行详细说明。
本发明提供了一种考虑提升量的生态流量确定方法,具体如图1所示,包括以下步骤:
步骤1:收集计算生态流量所需的相关数据,本实施例中收集的相关数据包括长系列历史流量数据、河段含沙量、输沙量、换水周期、多年平均蓄水量污染物浓度、目标水质浓度;
步骤2:从收集的相关数据中确定生态基流,确定方法包括改进蒙大拿法、年内展布法和最小月平均流量法;
步骤3:对改进蒙大拿法、年内展布法和最小月平均流量法三种方法分别确定的生态基流M1、M2和M3进行比较分析,选取最大值和最小值分别作为生态基流的上限Qmax和下限Qmin,初步确定生态基流的区间[Qmin,Qmax],将区间的下限Qmin初步确定为生态基流,见式(1)、(2);
Qmin=min{M1,M3,M3} (1)
Qmax=max{M1,M3,M3} (2)
步骤4:对步骤3的生态基流的下限Qmin进行验证;
根据2005年的《建设项目水资源论证导则(试行)》,北方河流生态基流指标原则上应不小于多年平均流量的10%,枯水时段不应低于同期流量均值的 20%,若初步确定的下限值Qmin大于导则规定,则初步确定的下限值Qmin为生态基流的下限值,若小于,则将导则规定作为生态基流的下限值;对于任一确定的考核断面而言,可以将验证过后所确定的生态基流作为考核的最低标准;
步骤5:在生态基流的基础上计算河段3个服务对象需水量,3个服务包括景观湿地需水Wl、输沙需水Ws和稀释自净需水Wc;
步骤6:对比分析步骤5计算的三个服务对象的需水量,确定提升量S;提升量是一个区间,因此找出最小值和最大值确定区间的提升量下限Smin和提升量上限Smax,为体现年内变化,提升量区间按月份分别确定;
Smin=min{Wl,Ws,Wc} (3)
Smax=max{Wl,Ws,Wc} (4)
例如三个服务对象为景观湿地需水B1、输沙需水B2、稀释自净需水B3,其在1月份的需水量分别为W’1=1.5、Q1=3.2、Wci1=2.7,其中最大值为Smax1=3.2,最小值为Smin1=1.5,那么1月份的提升量区间就是S1=[1.5,3.2],按照这种方法将2-12月的提升量区间全部确定,来体现提升量的年内变化;
步骤7:将生态基流和提升量S结合起来综合确定生态流量,将步骤4验证过后的生态基流作为考核的最低标准,也就是提升量的下限。
区间化管理时,首先要满足河道不断流的需求,即生态基流达到考核的最低标准,在此基础上通过增加提升量对生态进行改善和修复,生态需水水平不同,水资源调配能力不同,提升量可大可小;
生态需水水平的不同,比如:10m3/s的流量可以保证水生生物正常生长, 12m3/s的流量可以保证短时间内正常生长,15m3/s的流量则可以保证长时间内的健康生长,那么提升量首先满足当下的正常生长,再依次提升;
水资源调配能力的不同,比如1,水库在洪水期末没有蓄存够足够的水量,导致用水期水量不足,提升量难以保证;2,水库在用水期开始蓄存了足够的水量,可以将一部分水量用于生态,增加提升量;
一般情况下,生态流量自然越多越好,出现问题时,采取针对性的做法,在此后一段持续期内,尽快恢复和保障水生态;
例如:如果生态基流为8m3/s,提升量的区间为[5,14],最终提升量在5m3/s 的基础上提升,以生态基流8m3/s作为生态流量理论区间的下限,生态流量的理论区间变成[8,14],现在服务对象需水可增加流量的区间是[0,6],如果水资源调配能力对应的是[0,6]的2m3/s,那么就在生态基流8m3/s的基础上提升 2m3/s,那么生态流量就是10m3/s。
进一步地,本实施例中,步骤2中,改进蒙大拿法是在蒙大拿法的基础上考虑典型年的选取及分期百分比的确定,典型年的选取根据计算不同季节模数最接近1.0的年份,具体见式(17),分期是将一年分为一般用水期(11月- 次年3月)和鱼类产卵育幼期(4-10月),分别取典型年平均流量的10%和30%作为生态基流;
式中:Ki为第i年的模数,Kij是第i年第j季的模数,一般用于典型年的选取;
年内展布法是用多年月均径流量乘以同期均值比得到一年的生态基流过程,同期均值比是根据多年年均径流量与年内各月最小月均径流量的年均值这两个典型的水文特征变量进行确定与量化的,具体见式(5)-(10);
qmin(i)=min(qij) (8)
式中:为第i个月的多年月均径流量,为多年年均径流量,qij为第j年第i个月的月均径流量,qmin(i)为第i个月的多年最小月均径流量,为最小年均径流量,n为统计年数,η为同期均值比,Qi为各月的生态基流;
最小月平均流量法是将最小月平均流量的多年平均值作为河道的生态基流,具体见式(11);
式中:Qij为第i个月第j天的平均流量,Qm为第m个月的生态基流,n为统计年数。
进一步地,由于生态服务对象众多,无法细化到一草一木,生态需水难以量化,对于具体河道而言,为了说明不同河段的提升量,本实施例将河段划分为上游河段A1、中游河段A2、下游河段A3三部分;本实施例将服务对象分为景观湿地需水Wl、输沙需水Ws和稀释自净需水Wc,在生态基流的基础上考虑这三类给出提升量,例如:A1位于上游,水质较好,Ws、Wc就少,A2位于中游,两岸工厂较多,排污口就多,Wc就多,A3位于下游,泥沙淤积较为严重,那对应的Ws就多。河段A1-A3的各自提升量可能是不同的,但都可以用区间来表示;
步骤5中景观湿地需水Wl、输沙需水Ws和稀释自净需水Wc的具体计算方法如下:
(1)景观湿地需水:指为维持景观或者湿地生态和环境功能所要消耗和补充的水量,采用换水周期法进行计算,具体见式(12);
式中:T为换水周期,w”为多年平均蓄水量,Q为景观湿地需水流量;
(2)输沙需水:河流输沙需水量是指为了维持河床冲刷与淤积的动态平衡,必须在河道内保持一定的水量,该水量与水流含沙量直接相关,且随着产沙量、含沙量等的年内分配而动态变化,计算公式见式(13)-(15);
w'=ηα*ww (13)
ww=w-ws/γs (14)
η=s进-s出 (15)
式中:W′为输沙需水量,η′为输沙效率,α为指数,Ww为净水量;W为径流量,Ws为输沙量,γs为泥沙容重,通常取2.65t/m3;S进为进入河段的含沙量,S出为流出河段的含沙量;当η′<1时,河段处于冲刷状态,进口站含沙量小于出口站含沙量,此时,α=1;当η′>=1时,河段淤积或冲淤平衡,进口站含沙量大于或等于出口站含沙量,此时,α=0;
(3)稀释自净需水:指在加强污染源防治的前提下,使河流水体达到水功能区划确定的水质目标所要求的最小水量;采用混合模型与衰减模型对河段稀释自净需水量进行模拟;在河段i处,上断面为Ai-1,下断面为Ai,河段长为Li,假设在河段内共有多个排污口及支流的污染物汇入,为便于计算,记做断面B,具体见式(16);
式中:Wct为稀释自净需水量,Ci-1为Ai-1处的污染物浓度,si为B点处汇入的污染物浓度,qi为B点处的汇入流量,Csi为Ai段面处的目标水质浓度,Li为河段长、ui为河段i处的流速,xi为断面B到断面Ai-1的距离,Qi-1为Ai-1处的流量。
本实施例中,生态基流计算还可以采用水力学法、栖息地法或整体法;输沙需水计算还可采用冲淤修正法;景观需水计算还可采用湖泊形态法;生态基流的验证还可采用其他河道演算方法;分期还可采用丰、平、枯的方式;典型年的确定还可以采取频率分析法。
本实施例中提供的生态流量确定方法能够适应流域或区域考核管理的需要,在生态基流基础上给出提升量,提升量是考虑特定河段服务对象不同时空分布下的需水情况给出的需水水平,之后使用区间化的思想给出范围,以生态基流作为考核的标准,就可以保证河道不断流,这样考核是合理有效的,在提升量区间化管理的基础上通过调配决策来改善和修复水生态环境,达到生态管理的需求。
以上所述实施例仅为本发明较佳的具体实施方式,本发明的保护范围不限于此,任何熟悉本领域的技术人员在本发明披露的技术范围内,可显而易见地得到的技术方案的简单变化或等效替换,均属于本发明的保护范围。
Claims (4)
1.一种考虑提升量的生态流量确定方法,其特征在于,包括以下步骤:
步骤1:收集计算生态流量所需的相关数据;
步骤2:从收集的相关数据中确定生态基流,确定方法包括改进蒙大拿法、年内展布法和最小月平均流量法;
步骤3:对所述改进蒙大拿法、年内展布法和最小月平均流量法三种方法分别确定的生态基流M1、M2和M3进行比较分析,选取最大值和最小值分别作为生态基流的上限Qmax和下限Qmin,初步确定生态基流的区间[Qmin,Qmax],将区间的下限Qmin初步确定为生态基流,见式(1)、(2);
Qmin=min{M1,M2,M3} (1)
Qmax=max{M1,M2,M3} (2)
步骤4:对所述步骤3的生态基流的下限Qmin进行验证;
若初步确定的下限Qmin大于规定值,则下限Qmin为生态基流的下限值;否则将规定值作为生态基流的下限值;
步骤5:在所述生态基流的基础上计算河段3个服务对象需水量,3个服务包括景观湿地需水Wl、输沙需水Ws和稀释自净需水Wc;
步骤6:对比分析所述步骤5计算的三个服务对象的需水量,确定提升量S;提升量是一个区间,因此找出最小值和最大值确定区间的提升量下限Smin和提升量上限Smax,为体现年内变化,提升量区间按月份分别确定;
Smin=min{Wl,Ws,Wc} (3)
Smax=max{Wl,Ws,Wc} (4)
步骤7:将所述生态基流和提升量S结合起来综合确定生态流量,将步骤4验证过后的生态基流作为考核的最低标准,也就是提升量的下限。
2.根据权利要求1所述的考虑提升量的生态流量确定方法,其特征在于,所述步骤1中收集的相关数据包括河段含沙量、输沙量、换水周期、多年平均蓄水量、污染物浓度、目标水质浓度和长系列历史流量数据。
3.根据权利要求1所述的考虑提升量的生态流量确定方法,其特征在于,所述步骤2中,改进蒙大拿法是在蒙大拿法的基础上考虑典型年的选取及分期百分比的确定,典型年的选取根据计算不同季节模数最接近1.0的年份,分期是将一年分为一般用水期和鱼类产卵育幼期,分别取典型年平均流量的10%和30%作为生态基流;
年内展布法是用多年月均径流量乘以同期均值比得到一年的生态基流过程,同期均值比是根据多年年均径流量与年内各月最小月均径流量的年均值这两个典型的水文特征变量进行确定与量化的,具体见式(5)-(10);
qmin(i)=min(qij) (8)
式中:为第i个月的多年月均径流量,为多年年均径流量,qij为第j年第i个月的月均径流量,qmin(i)为第i个月的多年最小月均径流量,为最小年均径流量,n为统计年数,η为同期均值比,Qi为各月的生态基流;
最小月平均流量法是将最小月平均流量的多年平均值作为河道的生态基流,具体见式(11);
式中:Qij为第i个月第j天的平均流量,Qm为第m个月的生态基流,n为统计年数。
4.根据权利要求1所述的考虑提升量的生态流量确定方法,其特征在于,所述步骤5中景观湿地需水Wl、输沙需水Ws和稀释自净需水Wc的具体计算方法如下:
(1)景观湿地需水:指为维持景观或者湿地生态和环境功能所要消耗和补充的水量,采用换水周期法进行计算,具体见式(12);
式中:T为换水周期,w”为多年平均蓄水量,Q为景观湿地需水流量;
(2)输沙需水:河流输沙需水量是指为了维持河床冲刷与淤积的动态平衡,必须在河道内保持一定的水量,该水量与水流含沙量直接相关,且随着产沙量、含沙量等的年内分配而动态变化,计算公式见式(13)-(15);
w'=ηα*ww (13)
ww=w-ws/γs (14)
η=s进-s出 (15)
式中:W′为输沙需水量,η′为输沙效率,α为指数,Ww为净水量;W为径流量,Ws为输沙量,γs为泥沙容重,通常取2.65t/m3;S进为进入河段的含沙量,S出为流出河段的含沙量;当η′<1时,河段处于冲刷状态,进口站含沙量小于出口站含沙量,此时,α=1;当η′>=1时,河段淤积或冲淤平衡,进口站含沙量大于或等于出口站含沙量,此时,α=0;
(3)稀释自净需水:指在加强污染源防治的前提下,使河流水体达到水功能区划确定的水质目标所要求的最小水量;采用混合模型与衰减模型对河段稀释自净需水量进行模拟;在河段i处,上断面为Ai-1,下断面为Ai,河段长为Li,假设在河段内共有多个排污口及支流的污染物汇入,为便于计算,记做断面B,具体见式(16);
式中:Wci为稀释自净需水量,Ci-1为Ai-1处的污染物浓度,si为B点处汇入的污染物浓度,qi为B点处的汇入流量,Csi为Ai段面处的目标水质浓度,Li为河段长、ui为河段i处的流速,xi为断面B到断面Ai-1的距离,Qi-1为Ai-1处的流量。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911236748.0A CN111079085B (zh) | 2019-12-05 | 2019-12-05 | 一种考虑提升量的生态流量确定方法 |
PCT/CN2020/129023 WO2021109848A1 (zh) | 2019-12-05 | 2020-11-16 | 一种考虑提升量的生态流量确定方法 |
US17/361,327 US20210326408A1 (en) | 2019-12-05 | 2021-06-29 | Ecological flow determination method for considering lifting amount |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911236748.0A CN111079085B (zh) | 2019-12-05 | 2019-12-05 | 一种考虑提升量的生态流量确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111079085A true CN111079085A (zh) | 2020-04-28 |
CN111079085B CN111079085B (zh) | 2023-03-24 |
Family
ID=70313250
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911236748.0A Active CN111079085B (zh) | 2019-12-05 | 2019-12-05 | 一种考虑提升量的生态流量确定方法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US20210326408A1 (zh) |
CN (1) | CN111079085B (zh) |
WO (1) | WO2021109848A1 (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112184070A (zh) * | 2020-10-28 | 2021-01-05 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种协同生态流量需求的梯级水电站多目标优化调度方法及系统 |
WO2021109848A1 (zh) * | 2019-12-05 | 2021-06-10 | 西安理工大学 | 一种考虑提升量的生态流量确定方法 |
CN113343413A (zh) * | 2021-04-22 | 2021-09-03 | 中国环境科学研究院 | 水环境承载力评价方法、装置、设备和介质 |
CN116151592A (zh) * | 2023-04-17 | 2023-05-23 | 南昌工程学院 | 生态流量的确定和保障方法及系统 |
CN116563025A (zh) * | 2023-05-19 | 2023-08-08 | 生态环境部环境规划院 | 一种基于入河排污口的智能化排污权交易系统 |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111507540B (zh) * | 2020-04-28 | 2022-04-26 | 长江水利委员会长江科学院 | 一种资料贫乏地区长系列泥沙数据综合重构方法 |
CN113792957B (zh) * | 2021-08-02 | 2024-03-05 | 东北农业大学 | 一种基于水生生物对生存环境要求的生态稳定期识别方法 |
CN114438952B (zh) * | 2022-01-24 | 2023-05-16 | 中国长江三峡集团有限公司 | 水库沙峰调度的试验模拟系统与试验模拟方法 |
CN114548799B (zh) * | 2022-02-28 | 2022-12-13 | 中国环境科学研究院 | 一种三象限水生态管理耦合的方法及其应用 |
CN115392725B (zh) * | 2022-08-26 | 2023-08-01 | 中国长江三峡集团有限公司 | 一种生态效益评估方法、装置、存储介质及电子设备 |
CN115423182B (zh) * | 2022-08-31 | 2023-07-11 | 中国长江三峡集团有限公司 | 水电站下泄生态流量评价方法、装置、存储介质及设备 |
CN115795258B (zh) * | 2022-10-27 | 2023-11-28 | 华能伊敏煤电有限责任公司 | 一种定量分析改变取水方式对河流下游污染物浓度影响的方法 |
CN115860421B (zh) * | 2022-12-26 | 2023-08-29 | 长江水利委员会水文局 | 适配目标保证率的生态流量动态计算方法及系统 |
CN116933679B (zh) * | 2023-07-11 | 2024-05-21 | 广东省科学院广州地理研究所 | 流域的水生态安全格局构建方法、装置、设备及存储介质 |
CN117541031B (zh) * | 2024-01-09 | 2024-03-29 | 云南华电金沙江中游水电开发有限公司 | 一种全周期适宜生态流量区间推求方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101865687A (zh) * | 2009-07-31 | 2010-10-20 | 北京师范大学 | 流域生态需水整合分析技术 |
RU2606555C1 (ru) * | 2015-07-27 | 2017-01-10 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Уфимский государственный авиационный технический университет" | Способ определения среднегодовой самоочищающей способности водотока в фарватере |
CN108985576A (zh) * | 2018-06-26 | 2018-12-11 | 华北水利水电大学 | 一种河流生态系统服务功能的价值评估方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008185455A (ja) * | 2007-01-30 | 2008-08-14 | Institute Of Freshwater Biology | 水域環境の評価方法 |
CN101615285A (zh) * | 2009-08-04 | 2009-12-30 | 北京师范大学 | 一种生态环境需水阈值的计算方法 |
CN108876022A (zh) * | 2018-05-31 | 2018-11-23 | 北京东方园林环境股份有限公司 | 一种用于确定天然湿地恢复生态环境需水量的计算方法 |
CN109784527A (zh) * | 2018-12-05 | 2019-05-21 | 北京师范大学 | 一种确定河流生态系统对流量和水质需求的方法 |
CN111079085B (zh) * | 2019-12-05 | 2023-03-24 | 西安理工大学 | 一种考虑提升量的生态流量确定方法 |
-
2019
- 2019-12-05 CN CN201911236748.0A patent/CN111079085B/zh active Active
-
2020
- 2020-11-16 WO PCT/CN2020/129023 patent/WO2021109848A1/zh active Application Filing
-
2021
- 2021-06-29 US US17/361,327 patent/US20210326408A1/en active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101865687A (zh) * | 2009-07-31 | 2010-10-20 | 北京师范大学 | 流域生态需水整合分析技术 |
RU2606555C1 (ru) * | 2015-07-27 | 2017-01-10 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Уфимский государственный авиационный технический университет" | Способ определения среднегодовой самоочищающей способности водотока в фарватере |
CN108985576A (zh) * | 2018-06-26 | 2018-12-11 | 华北水利水电大学 | 一种河流生态系统服务功能的价值评估方法 |
Non-Patent Citations (3)
Title |
---|
侯世文: "基于多种水文学法分析大汶河干流生态基流", 《水文》 * |
卫凤莲等: "延安市河流系统生态环境需水量测评与分析", 《陕西师范大学学报(自然科学版)》 * |
潘扎荣等: "河道基本生态需水的年内展布计算法", 《水利学报》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2021109848A1 (zh) * | 2019-12-05 | 2021-06-10 | 西安理工大学 | 一种考虑提升量的生态流量确定方法 |
CN112184070A (zh) * | 2020-10-28 | 2021-01-05 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种协同生态流量需求的梯级水电站多目标优化调度方法及系统 |
CN112184070B (zh) * | 2020-10-28 | 2021-12-17 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种协同生态流量需求的梯级水电站多目标优化调度方法及系统 |
CN113343413A (zh) * | 2021-04-22 | 2021-09-03 | 中国环境科学研究院 | 水环境承载力评价方法、装置、设备和介质 |
CN113343413B (zh) * | 2021-04-22 | 2023-06-30 | 中国环境科学研究院 | 水环境承载力评价方法、装置、设备和介质 |
CN116151592A (zh) * | 2023-04-17 | 2023-05-23 | 南昌工程学院 | 生态流量的确定和保障方法及系统 |
CN116563025A (zh) * | 2023-05-19 | 2023-08-08 | 生态环境部环境规划院 | 一种基于入河排污口的智能化排污权交易系统 |
Also Published As
Publication number | Publication date |
---|---|
WO2021109848A1 (zh) | 2021-06-10 |
CN111079085B (zh) | 2023-03-24 |
US20210326408A1 (en) | 2021-10-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111079085B (zh) | 一种考虑提升量的生态流量确定方法 | |
CN110189059A (zh) | 一种流域水系统协同承载力评价指标体系构建方法 | |
CN108305006B (zh) | 基于供水保证率约束的水库水资源使用权量初始分配方法 | |
CN107622363B (zh) | 城市水环境质量评估方法 | |
CN108171001A (zh) | 一种评估水利水电工程生态流量下泄效果的方法 | |
CN109598428A (zh) | 一种基于行政单元和水系的污染物削减分配方法 | |
Yang et al. | Development of a long-term, ecologically oriented dam release plan for the Lake Baiyangdian sub-basin, Northern China | |
CN104679993A (zh) | 一种基于二元水循环的纳污能力计算方法 | |
CN115953061A (zh) | 一种水资源短缺度估算方法、装置及计算机可读存储介质 | |
CN113177727B (zh) | 一种河道分级分期旱限水位或流量的计算方法 | |
CN109002946B (zh) | 一种河湖补水的“两库—两站”系统水资源优化调度方法 | |
Deksissa et al. | Control options for river water quality improvement: A case study of TDS and inorganic nitrogen in the Crocodile River (South Africa) | |
O'Keeffe | Environmental flow allocation as a practical aspect of IWRM | |
Abou El Hassan | Improving water quality in the Nile Delta irrigation network by regulating reuse of agricultural drainage water | |
Shelby | Streamflow and recreation | |
CN111597693A (zh) | 一种入河污染负荷最优分配方法及系统 | |
Wang et al. | Study on benefits evaluation of water diversion project: Case study in water transfer from the Yangtze River to Lake Taihu | |
CN114912664A (zh) | 基于水质-水位二元响应关系的适宜生态水位预测方法 | |
CN113449983A (zh) | 一种新型城市湿地规划设计方法 | |
Li et al. | Effect of pollution load reduction on water quality in rural lakes in the shallow hill water network area | |
Chavula | The evaluation of the present and potential water resources management for the Lake Chilwa basin | |
CN112906313A (zh) | 基于系统协同理论的生态流量计算方法 | |
KR20160088030A (ko) | 다중수원 지능적 배분공급 스케줄링 방법 | |
Kamal | Environmental flows Indus river system in Pakistan | |
CN115330259A (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 | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20231218 Address after: 710000 1-1-918, Wanda Plaza, Keji 7th Road, Zhangba street, high tech Zone, Xi'an, Shaanxi Province Patentee after: Xi'an Xinghan Chuangyi Information Technology Co.,Ltd. Address before: Xi'an University of technology, No.5, Jinhua South Road, Beilin District, Xi'an City, Shaanxi Province Patentee before: XI'AN University OF TECHNOLOGY |