CN107977505A - 一种前期降水衰退系数k确定的新方法 - Google Patents
一种前期降水衰退系数k确定的新方法 Download PDFInfo
- Publication number
- CN107977505A CN107977505A CN201711213384.5A CN201711213384A CN107977505A CN 107977505 A CN107977505 A CN 107977505A CN 201711213384 A CN201711213384 A CN 201711213384A CN 107977505 A CN107977505 A CN 107977505A
- Authority
- CN
- China
- Prior art keywords
- coefficient
- precipitation
- antecedent
- antecedent precipitation
- scenes
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Testing And Monitoring For Control Systems (AREA)
Abstract
一种前期降水衰退系数k确定的新方法,包括:步骤一:在0.8‑0.98区间,设置不同情景的衰退系数k。步骤二:计算所有k情景的前期降水指数API值,获得不同k情景的前期降水指数API序列。步骤三:计算研究对象如枯水与前期降水指数API的相关系数。步骤四:分析不同k情景的相关系数,如果前期某一天之前的所有相关系数均为零,表明该天之前的前期降水对当前天的研究对象如枯水已不再产生影响,前期降水与当前天的枯水呈零相关,此时对应的衰退系数k即为前期降水的最佳衰退系数。该方法从机理上保证了所确定的衰退系数k的合理性与科学性,克服了凭主观经验选取衰退系数k的不足,计算成果更准确、更科学,有重要的理论意义和实用价值,应用前景广阔。
Description
技术领域
本发明涉及与流域下垫面有关的前期降水影响领域,尤其涉及一种前期降水衰退系数k确定的新方法。
背景技术
降水是水循环过程的重要环节。在降水完成海陆循环的过程中,由于下垫面对降水的截留作用减缓了降水的海陆循环过程,将前期吸收的降水储存一段时间,然后以蒸发、潜水、地下水等形式完成水循环。一方面前期降水对后期的土壤湿度、径流等产生重要影响,另一方面下垫面每次吸收的水量差别各异,导致不同区域的相同降水、或者相同区域的相同降水在不同时期结果不相同。因此,在研究降水对土壤湿度和径流尤其是枯水径流的影响时,必须考虑前期降水的影响。为量化前期降水的影响,M.A.Kohler于20世纪40年代提出前期降水指数(Antecedent precipitation index,API),1951年M.A.Kohler和R.K.Linsley给出如下具体计算公式:
式中,Pt是前期第t天的降水量,T为前期天数,k为衰退系数。
前期降水指数API的值用前期降水量与衰退系数k的乘积的累积和表示,由于前期降水是已知的,API的值取决于衰退系数k,k值越大,表示前期降水衰退越慢、在下垫面中储存量的比例越大,前期降水对后期的影响越大。因此,确定合适的k值对准确计算API和量化前期降水的影响意义重大。
目前国内外的学者多沿用M.A.Kohler和R.K.Linsley给出的k经验取值范围0.85-0.98,以及W.Viessman和G.L.Lewis给出的k取值范围0.8-0.98,之后国内外的学者均采用此经验值,并广泛应用于前期降水对后期的径流、土壤湿度、滑坡、泥石流、极端干旱和洪涝灾害等的影响研究中,但均没有准确计算衰退系数k的可靠方法。
不同的研究区下垫面差异很大,对前期降水的储存和衰退能力差异较大,k值差异较大;相同研究区的不同时期,下垫面的初始状态也受到前期降水的影响,k值差异也很大。而前期降水的影响程度取决于衰退系数k,依照传统方法取经验值具有很大的人为主观性,得到的前期降水指数API值也难以准确量化前期降水对后期水文要素如径流、土壤湿度等的影响,对后期山洪暴发临界雨量、极端干旱或洪涝灾害预报等也会出现极大误差。因此,准确量化k值具有重要的科学意义和实用价值,应用前景广阔。
发明内容
针对上述问题,本发明的目的是提供一种前期降水衰退系数k确定的新方法。它的核心思想在于当前期降水足够远时,其对当前的水文要素如径流、土壤湿度等不再产生影响,两者呈现零相关。本发明首先设定不同情景的衰退系数k方案;根据公式(1)算出所有k情景的前期降水指数API值;然后计算不同k情景的研究对象如枯水序列与API序列的相关系数;最后分析不同k情景的相关系数,如果前期某一天之前的所有相关系数均为零,表明该天之前的前期降水对当前天的研究对象如枯水已不再产生影响,进而确定前期降水的最佳衰退系数k值。
为解决上述问题,本发明采取以下技术方案:
一种前期降水衰退系数k确定的新方法,它适用于不同研究领域、不同时间尺度的前期降水衰退系数k的准确计算。该方法具体步骤如下:
步骤一:在0.8-0.98区间,设置不同情景的衰退系数k。
步骤二:根据前期降水指数API计算公式(1),计算所有k情景的前期降水指数API值,获得不同k情景的前期降水指数API序列。
步骤三:计算研究对象如枯水序列与前期降水指数API序列的相关系数。
步骤四:分析不同k情景的相关系数,如果前期某一天之前的所有相关系数均为零,表明该天之前的前期降水对当前天的研究对象如枯水已不再产生影响,前期降水与当前天的枯水呈零相关,此时对应的衰退系数k即为前期降水的最佳衰退系数。
本发明由于采取以上技术方案,其具有以下优点:
1.基于前期降水与当前天研究对象如枯水的零相关性来确定最佳的衰退系数k,从机理上保证了所确定的衰退系数k的合理性与科学性,克服了凭主观经验选取衰退系数k的不足。
2.该方法简单易操作,条理清楚,计算成果更准确、更科学。
3.该方法具有更好的适用性,既适用于前期降水对当前径流影响的情况,也适用于前期降水对当前土壤湿度影响等的情况,有重要的理论意义和实用价值,应用前景广阔。
附图说明
图1为本发明方法的流程框图。
图2为泾河流域张家山水文站1956-2010年日枯水与流域前期降水指数之间不同衰退系数k情景下的相关系数。
图3为确定的泾河张家山水文站以上流域前期降水的最佳衰退系数k=0.9675。
具体实施方式
下面结合附图和具体实施方式,对本发明作进一步详细说明。
如图1所示,本发明一种前期降水衰退系数k确定的新方法,包括设置0.8-0.98区间不同情景的衰退系数k值、确定不同衰退系数k的前期降水指数API值、计算研究对象与前期降水指数API的相关系数以及获得最佳的衰退系数k四部分。
以水文系统的日枯水和前期1-365天的日降水为例,本发明的具体实施按照以下步骤进行:
步骤一:在0.8-0.98区间,取计算步长为0.0001,设置不同情景的衰退系数k值。
步骤二:根据前期降水指数API计算公式(1),计算所有k情景的前期降水指数API值,获得不同k情景的前期降水指数API序列。
步骤三:计算枯水序列与前期降水指数API序列的相关系数。
步骤四:分析不同k情景的相关系数,如果前期某一天之前的所有相关系数均为零,表明该天之前的前期降水对当前天的枯水已不再产生影响,前期降水与当前天的枯水呈零相关,此时对应的衰退系数k即为前期降水的最佳衰退系数。
实施案例
本发明以泾河流域张家山水文站1956-2010年的日枯水和张家山水文站以上流域的日降水为研究对象,在0.8-0.98区间,取计算步长为0.0001,设置不同情景的衰退系数k值,根据式(1)计算不同衰减系数k情景的前期1-365天的前期降水指数API,确定枯水与前期降水指数API的相关系数,从而根据前期某一天之前的相关系数均为零的特点确定前期降水的最佳衰退系数k值。
其结果,分别见图2、图3。
图2为泾河流域张家山水文站1956-2010年日枯水与流域前期降水指数之间不同衰退系数k情景下的相关系数。
图3为确定的泾河张家山水文站以上流域前期降水的最佳衰退系数k=0.9675。
从上述实例可以看出,本发明提供的一种前期降水衰退系数k确定的新方法,基于前期日降水与当前日枯水相关性为零的特点,从机理上保证了所确定的衰退系数k的合理性与科学性,克服了凭主观经验选取衰退系数k的不足,有重要的理论意义和实用价值,应用前景广阔。
以上所述实施例仅表达了本发明的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对本发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。因此,本发明专利的保护范围应以所附权利要求为准。
Claims (5)
1.一种前期降水衰退系数k确定的新方法,包括以下步骤:
步骤一:在0.8-0.98区间,设置不同情景的衰退系数k;
步骤二:计算所有k情景的前期降水指数API值,获得不同k情景的前期降水指数API序列;
步骤三:计算研究对象如枯水序列与前期降水指数API序列的相关系数;以及
步骤四:分析不同k情景的相关系数,如果前期某一天之前的所有相关系数均为零,表明该天之前的前期降水对当前天的研究对象如枯水已不再产生影响,前期降水与当前天的枯水呈零相关,此时对应的衰退系数k即为前期降水的最佳衰退系数。
2.根据权利要求1所述的一种前期降水衰退系数k确定的新方法,其特征在于步骤一中:在0.8-0.98区间,设置不同情景的衰退系数k。由于衰退系数k对前期降水指数API有很大的影响,因此把衰退系数k的计算步长设置得小一些,设置计算步长为0.0001。
3.根据权利要求2所述的一种前期降水衰退系数k确定的新方法,其特征在于步骤二中:计算所有k情景的前期降水指数API值,获得不同k情景的前期降水指数API序列。前期降水指数API的计算公式为式中,Pt是前期第t天的降水量,T为前期天数,k为衰退系数。比如研究前期日降水对当前日枯水的影响,可以计算出不同k情景下的前期1-365天的前期降水指数序列。
4.根据权利要求3所述的一种前期降水衰退系数k确定的新方法,其特征在于步骤三中:计算研究对象如枯水序列与前期降水指数API序列的相关系数。比如研究前期日降水对当前日枯水的影响,可以计算出不同k情景下的前期1-365天的前期降水指数序列与当前日枯水序列的相关系数。
5.根据权利要求4所述的一种前期降水衰退系数k确定的新方法,其特征在于步骤四中:分析不同k情景的相关系数,如果前期某一天之前的所有相关系数均为零,表明该天之前的前期降水对当前天的研究对象如枯水已不再产生影响,前期降水与当前天的枯水呈零相关,此时对应的衰退系数k即为前期降水的最佳衰退系数。比如研究前期日降水对当前日枯水的影响,可以分析不同衰退系数k情景下的1-365天前期降水指数序列与当前日枯水序列之间的相关系数,当前期某天之前所有的相关系数均为零时,前期降水与当前天的枯水呈零相关,则此时对应的k值即为前期降水的最佳衰退系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711213384.5A CN107977505A (zh) | 2017-11-28 | 2017-11-28 | 一种前期降水衰退系数k确定的新方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711213384.5A CN107977505A (zh) | 2017-11-28 | 2017-11-28 | 一种前期降水衰退系数k确定的新方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107977505A true CN107977505A (zh) | 2018-05-01 |
Family
ID=62012034
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711213384.5A Pending CN107977505A (zh) | 2017-11-28 | 2017-11-28 | 一种前期降水衰退系数k确定的新方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107977505A (zh) |
Citations (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB1356700A (en) * | 1970-12-14 | 1974-06-12 | Texaco Development Corp | Automatically controlling mole ratios |
EP1689383A2 (en) * | 2003-11-19 | 2006-08-16 | Metabasis Therapeutics, Inc. | Novel phosphorus-containing thyromimetics |
US20090174142A1 (en) * | 2008-01-09 | 2009-07-09 | Sullivan Richard J | Methods and apparatus for educational spelling games |
US20100131248A1 (en) * | 2008-11-26 | 2010-05-27 | Microsoft Corporation | Reference model for data-driven analytics |
JP2010257091A (ja) * | 2009-04-23 | 2010-11-11 | Fujitsu Ltd | 危険度判定プログラム、危険度判定装置及び方法 |
US20110308996A1 (en) * | 2010-06-22 | 2011-12-22 | Conocophillips Company | Methodology to determine feed quality (di+ring aromatic content) of fcc and hydrocracking feeds |
US20150206255A1 (en) * | 2011-05-13 | 2015-07-23 | HydroBio, Inc | Method and system to prescribe variable seeding density across a cultivated field using remotely sensed data |
CN105604544A (zh) * | 2014-11-10 | 2016-05-25 | 中国石油化工股份有限公司 | 储层水敏感性室内评价方法 |
US20160160598A1 (en) * | 2014-12-05 | 2016-06-09 | National Oilwell Varco, L.P. | Method of closing a blowout preventer seal based on seal erosion |
CN106021366A (zh) * | 2016-05-10 | 2016-10-12 | 浙江大学 | 一种基于异构信息的api标签推荐方法 |
CN106202163A (zh) * | 2016-06-24 | 2016-12-07 | 中国环境科学研究院 | 通江湖泊生态监测信息管理及预警系统 |
CN106777913A (zh) * | 2016-11-29 | 2017-05-31 | 兰州大学 | 一种近似熵和样本熵共同最优参数m、r确定的新方法 |
CN106874544A (zh) * | 2017-01-05 | 2017-06-20 | 西南石油大学 | 一种页岩储层改造体积的地质表征方法 |
US20170177763A1 (en) * | 2015-12-18 | 2017-06-22 | Triple Ring Technologies, Inc. | Method and apparatus for x-ray ionizing radiation control |
WO2017168425A1 (en) * | 2016-03-31 | 2017-10-05 | Technion Research & Development Foundation Limited | Method for separation of magnesium and calcium ions from saline water, for improving the quality of soft and desalinated waters |
CN107392509A (zh) * | 2017-09-01 | 2017-11-24 | 河海大学 | 一种基于物元分析的河道生态径流过程评价方法 |
-
2017
- 2017-11-28 CN CN201711213384.5A patent/CN107977505A/zh active Pending
Patent Citations (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB1356700A (en) * | 1970-12-14 | 1974-06-12 | Texaco Development Corp | Automatically controlling mole ratios |
EP1689383A2 (en) * | 2003-11-19 | 2006-08-16 | Metabasis Therapeutics, Inc. | Novel phosphorus-containing thyromimetics |
US20090174142A1 (en) * | 2008-01-09 | 2009-07-09 | Sullivan Richard J | Methods and apparatus for educational spelling games |
US20100131248A1 (en) * | 2008-11-26 | 2010-05-27 | Microsoft Corporation | Reference model for data-driven analytics |
JP2010257091A (ja) * | 2009-04-23 | 2010-11-11 | Fujitsu Ltd | 危険度判定プログラム、危険度判定装置及び方法 |
US20110308996A1 (en) * | 2010-06-22 | 2011-12-22 | Conocophillips Company | Methodology to determine feed quality (di+ring aromatic content) of fcc and hydrocracking feeds |
US20150206255A1 (en) * | 2011-05-13 | 2015-07-23 | HydroBio, Inc | Method and system to prescribe variable seeding density across a cultivated field using remotely sensed data |
CN105604544A (zh) * | 2014-11-10 | 2016-05-25 | 中国石油化工股份有限公司 | 储层水敏感性室内评价方法 |
US20160160598A1 (en) * | 2014-12-05 | 2016-06-09 | National Oilwell Varco, L.P. | Method of closing a blowout preventer seal based on seal erosion |
US20170177763A1 (en) * | 2015-12-18 | 2017-06-22 | Triple Ring Technologies, Inc. | Method and apparatus for x-ray ionizing radiation control |
WO2017168425A1 (en) * | 2016-03-31 | 2017-10-05 | Technion Research & Development Foundation Limited | Method for separation of magnesium and calcium ions from saline water, for improving the quality of soft and desalinated waters |
CN106021366A (zh) * | 2016-05-10 | 2016-10-12 | 浙江大学 | 一种基于异构信息的api标签推荐方法 |
CN106202163A (zh) * | 2016-06-24 | 2016-12-07 | 中国环境科学研究院 | 通江湖泊生态监测信息管理及预警系统 |
CN106777913A (zh) * | 2016-11-29 | 2017-05-31 | 兰州大学 | 一种近似熵和样本熵共同最优参数m、r确定的新方法 |
CN106874544A (zh) * | 2017-01-05 | 2017-06-20 | 西南石油大学 | 一种页岩储层改造体积的地质表征方法 |
CN107392509A (zh) * | 2017-09-01 | 2017-11-24 | 河海大学 | 一种基于物元分析的河道生态径流过程评价方法 |
Non-Patent Citations (4)
Title |
---|
AIHUA JIANG ET,: "Utilization of Renewable Energy Extracted from Reservoirs in Air Conditioning System: Case study", 《IEEE XPLORE》 * |
L.DESCROIX ET,: "Evalution of an antecedent precipitation index to model runoff yield in the western sierra madre", 《JOURNAL OF HYDROLOGY》 * |
时兴合 等,: "黄河上游径流变化特征及其影响因素初步分析", 《中国沙漠》 * |
杨春辉 等,: "改进的API 模型在尼尔基洪水预报中的应用", 《东北水利水电》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Perrin et al. | Improvement of a parsimonious model for streamflow simulation | |
Setegn et al. | Modeling of Sediment Yield From Anjeni‐Gauged Watershed, Ethiopia Using SWAT Model 1 | |
Jajarmizadeh et al. | Application of SVM and SWAT models for monthly streamflow prediction, a case study in South of Iran | |
Yue et al. | Unit hydrographs to model quick and slow runoff components of streamflow | |
Li et al. | Correlating check dam sedimentation and rainstorm characteristics on the Loess Plateau, China | |
Deng et al. | Estimation of nonfluctuating reservoir inflow from water level observations using methods based on flow continuity | |
Wang et al. | Improvement of forecast skill for severe weather by merging radar-based extrapolation and storm-scale NWP corrected forecast | |
Lu et al. | Water balance assessment of an ungauged area in Poyang Lake watershed using a spatially distributed runoff coefficient model | |
Chaube et al. | Synthesis of flow series of tributaries in Upper Betwa basin | |
Turhan | A comparative evaluation of the use of artificial neural networks for modeling the rainfall-runoff relationship in water resources management | |
Hui et al. | The Fractal Characteristics of Drainage Networks and Erosion Evolution Stages of Ten Kongduis in the Upper Reaches of the Yellow River, China | |
Zhou et al. | Modified CN method for small watershed infiltration simulation | |
CN106446359A (zh) | 基于流式大数据处理方式的降雨径流预测计算方法 | |
CN107977505A (zh) | 一种前期降水衰退系数k确定的新方法 | |
Hassanzadeh et al. | Performance evaluation of correction coefficients to optimize sediment rating curves on the basis of the Karkheh dam reservoir hydrography, west Iran | |
CN108920571B (zh) | 基于分布式水文模型的室内水文实验模型比例尺确定方法 | |
Bao et al. | Hydraulic model with roughness coefficient updating method based on Kalman filter for channel flood forecast | |
Fathy et al. | Runoff hydrograph modeling for arid regions: case study—Wadi Sudr-Sinai | |
Dalavi et al. | Assessment of empirical methods for runoff estimation in Chaskaman catchment of Western Maharashtra, India | |
Vahidi | Bivariate analysis of river flow and suspended sediment load in Aharchai Basin, Iran | |
Jiao et al. | Conceptualizing antecedent runoff condition using recurrence relation to modify SCS model | |
Adgolign et al. | Assessment of Spatio-temporal occurrence of water resources in Didissa Sub-Basin, West Ethiopia | |
Khidir | Evaluation of the Storms Direct Runoff Prediction Methods used for Goizha-Dabashan Watershed | |
Shekar | Rainfall-Runoff Modelling of a River Basin Using HEC HMS: A Review Study | |
Yang et al. | Applicability evaluation of the SWIM at river basins of the black soil region in Northeast China: A case study of the upper and middle Wuyuer River basin |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20180501 |