CN102663256A - 确定暴雨衰减指数为函数型的暴雨强度公式的一种方法 - Google Patents
确定暴雨衰减指数为函数型的暴雨强度公式的一种方法 Download PDFInfo
- Publication number
- CN102663256A CN102663256A CN2012101146942A CN201210114694A CN102663256A CN 102663256 A CN102663256 A CN 102663256A CN 2012101146942 A CN2012101146942 A CN 2012101146942A CN 201210114694 A CN201210114694 A CN 201210114694A CN 102663256 A CN102663256 A CN 102663256A
- Authority
- CN
- China
- Prior art keywords
- formula
- parameter
- lgp
- relation
- rainfall
- 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
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
确定暴雨衰减指数为函数型的暴雨强度公式的一种方法,属于市政工程技术领域。根据重现期P,暴雨强度i,历时t的关系值,首先,采用高斯-牛顿法确定各单一重现期分公式的参数;其次,通过回归计算选配暴雨衰减指数n与重现期P为半对数型或幂函数型相关关系式,并取其相关指数R2较大者作为采用的n与P的关系式,同时,回归计算雨力A与重现期P关系式中的参数,并将上述求得的各个参数值作为各重现期总公式中相应参数的初值,对于参数b,取各单一重现期的平均值作为总公式中b的初值;然后,采用高斯-牛顿法确定暴雨衰减指数n为函数型的总公式的参数,进而得暴雨强度公式。该公式是城市排水工程设计中确定雨水设计流量的基本依据。
Description
一、技术领域
本发明属于市政工程技术领域。暴雨强度公式是城市排水工程设计中确定雨水设计流量的基本依据之一,该公式的合理性与精度直接影响城市排水工程的规模。
二、背景技术
目前,我国城市暴雨强度公式的一般形式为:
单一重现期的分公式
各重现期的总公式
式中i为暴雨强度,简称雨强,mm/min;P为重现期,a;t为降雨历时,min。而A,A1,C,b,n是与地方暴雨特性有关的参数,其中A为雨力,即t+b=1时的雨强,mm/min;A1为重现期为1a的设计降雨的雨力,mm/min;C为雨力变动参数,b为降雨历时修正参数,min;n为暴雨衰减指数。
单一重现期的分公式对降雨历时t,雨强i的关系值(t,i)的拟合精度较高,但只能应用于一个重现期的情况,故目前我国常用式(2)形式的成果,且其参数n,b采用固定值。
文献[1]提出了采用高斯-牛顿法确定式(2)公式中的参数。
高斯牛顿法[2]:
设一般的非线性方程为
y=f(X,θ)+ε
式中f为一般函数;X可以是单个自变量,也可以是r个自变量X=(x1,x2,…,xr);θ为P维参数向量,即θ=(θ1,θ2,…,θp)′;ε为随机误差项,且∈N(0,σ2)。设y和X具有m组观察值(x1i,x2i,…,xri;yi),i=1~m。求非线性模型参数的“最小二乘”估计,就是求θ的待估计值 使残差平方和
为最小。
对于非线性模型,无法直接求参数的“最小二乘”解。高斯-牛顿法将f(X,θ)在参数初值 θ(0)=(θ1(0),θ2(0),…,θp(0))′处展开成只包含一次项的泰勒级数,从而使非线性模型线性化。为方便起见,用记号fi(θ)代替f(xi,θ),可导出满足式(3)Q(θ)为最小的参数递推公式,并写成矩阵形式为[2]
θ(k+1)=θ(k)+[J′(θ(k))J(θ(k))]-1 J′(θ(k))[(y-f(θ(k))] (4)
式中k为递推次数。
y=(y1,y2,…,ym)′,f(θ(k))=[f1(θ(k)),f2(θ(k)),…,fm(θ(k))]。
利用式(4)从参数初值θ(0)开始,一步步递推下去,直到θ(k)收敛稳定,即|θ(k+1)-θ(k)|小于或等于预先指定的小正数δ,从而得到θ的估计值。
对于暴雨强度公式的拟合精度,GB50014-2006《室外排水设计规范》[3](2011版)附录A第A.0.10条指出:“计算抽样误差和暴雨强度公式均方差。宜按绝对均方差计算,也可辅以相对均方差计算。计算重现期在2年~20年时,在一般强度的地方,平均绝对均方差不宜大于0.05mm/min(注:原文为“平均绝对方差不宜大于0.05mm/min”。显然是错误的,均方差才与暴雨强度具有相同量纲)。在较大强度的地方,平均相对均方差不宜大于5%”。
一些地区的计算成果表明,即使采用非线性最小二乘法或其他最优化方法求解式(2)的拟合参数,某些重现期或总公式的平均绝对均方差(即标准差)仍超过0.05mm/min[4,5],特别是参数n随重现期变化较为显著的地区,这种现象尤为突出。文献[4]指出:为了更密切的吻合当地多年统计的降雨强度子样点分布规律,提高精度,总公式的形式也可采用(第931页):
文献[4]还指出:暴雨衰减指数n值精确与否对公式的均方差有很大的影响(第929页)。
因此,从既不使公式过于复杂,又提高拟合精度的角度,可采用暴雨衰减指数n为函数型的暴雨强度公式:
式中C1=A1×C。
文献[4]给出了n与P的关系式为n=n1+n2lgP时确定式(6)中参数的手工算法(第932-939页)。到目前为止,尚无采用解析法确定暴雨衰减指数为函数型的暴雨强度公式,即总公式式(6)的方法。
三、发明内容
为了克服手工算法确定式(6)中暴雨衰减指数的函数形式及式(6)中参数效率低、所求参数具有一定的任意性、拟合精度不如解析法高的缺陷,本发明采用高斯-牛顿法确定暴雨衰减指数为函数型的暴雨强度公式,即总公式式(6)。
本发明采用的技术方案是:根据重现期P,暴雨强度i,历时t的关系值,首先,采用高斯-牛顿法确定各单一重现期分公式(1)的参数A,b,n;其次,通过回归计算选配暴雨衰减指数n与重现期P为半对数型或幂函数型的相关关系式,并取其相关指数R2较大者作为采用的n与P的关系式,同时,回归计算雨力A与重现期P关系式A=A1+C1lgP中的参数,并将上述求得的各个参数值作为总公式中相应参数的初值,对于参数b,取各单一重现期的平均值作为总公式中b的初值;然后,采用高斯-牛顿法确定暴雨衰减指数n为函数型的总公式(6)的参数,进而得暴雨强度公式。
本项发明的有益效果是,计算效率高、拟合精度高,使所求参数为最佳拟合参数,且暴雨衰减指数n的函数形式的确定科学合理。本项发明已利用我国有关城市的暴雨特征值资料进行了验证,效果很好。
四、具体实施方式
已知某城市重现期P、降雨强度i和降雨历时t的m组关系值(Pj,tj;ij),j=1~m,其中各个单一重现期有m1组观察值。
1.确定式(6)中暴雨衰减指数n与重现期P的关系式以及式(6)中各参数的初值。
(1)用高斯-牛顿法确定各单一重现期分公式(1)的参数A,b,n
1)对式(1)参数A,b,n分别求偏导数
2)确定各单一重现期分公式(1)的参数A,b,n的初值。具体方法可参考文献[1],并将求得的三参数的初值表示为:θ(0)=(A(0),b(0),n(0))′。
3)应用式(7)~式(9)、θ(0)以及每一重现期相应的m1组观察值(tj,ij),j=1~m1,计算偏导数矩阵J(θ(0))以及矩阵f(θ(0))。
4)根据式(4)计算θ(1)。
5)再以θ(1)作为参数的初始值θ(0),重复步骤3),步骤4),根据给定的δ值(例如,δ=0.0005),经若干次递推迭代,可求得各单一重现期分公式(1)参数的估计值θ=(A,b,n)′。
(3)根据参数A与重现期P的关系值,确定A=A1+C1lgP关系式中参数A1、C1,并分别将其作为式(6)中参数A1、C1的初值。
(4)将各单一重现期的参数b,取平均值作为式(6)中参数b的初值。
至此,得暴雨衰减指数采用n=n1+n2lgP或 时总公式(6)参数的初值θ(0)=(A1(0),C1(0),b(0),n1(0),n2(0))。
2.用高斯-牛顿法确定总公式(6)的参数。
当n-P关系采用n=n1+n2lgP形式时,总公式式(6)为:
1)对式(10)中的参数A1,C1,b,n1,n2分别求偏导数:
2)由前述确定的总公式参数的初值θ(0)=(A1(0),C1(0),b(0),n1(0),n2(0))′、式(11)~式(15)以及全部重现期的m组观察值(Pj,tj;ij),j=1~m,计算偏导数矩阵J(θ(0))以及矩阵f(θ(0))。
3)根据式(4)计算θ(1)。
4)再以θ(1)作为参数的初始值θ(0),重复步骤2),步骤3),根据给定的δ值(例如,δ=0.0005),经若干次递推迭代,可求得总公式(10)中参数的估计值θ=(A1,C1,b,n1,n2)′。
1)对式(16)中的参数A1,C1,b,n1,n2分别求偏导数:
第2)步~第4)步与推求公式(10)中参数的方法类似,不再赘述。
有关文献
[1]张子贤.用高斯-牛顿法确定暴雨公式参数[J].河海大学学报,1995,23(5):106-111.
[2]方开泰.实用回归分析[M].北京:科学出版社,1988.168-172.
[3]GB50014-2006.室外排水设计规范(2011版)[s].
[4]北京市市政工程设计研究总院.给水排水设计手册,第5册,城镇排水,第二版[M].北京:中国建筑工业出版社,2002.921-975.
[5]樊建军,王峰,陈鹏飞.利用MATLAB推导城市暴雨强度公式[J].中国给水排水,2010,26(11):114-115。
Claims (3)
1.确定暴雨衰减指数为函数型的暴雨强度公式的一种方法,根据各重现期P,暴雨强度i,历时t的关系值,确定各重现期的总公式
其特征是:根据暴雨衰减指数n与重现期P的关系值,选配n与P的关系式,并采用高斯-牛顿法确定暴雨衰减指数为函数型的总公式的参数。
3.根据权利要求1所述的确定暴雨衰减指数为函数型的暴雨强度公式的一种方法,其特征是:采用高斯-牛顿法确定暴雨衰减指数n为函数型的总公式
中的参数,或采用高斯-牛顿法确定暴雨衰减指数n为函数型的总公式
中的参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210114694.2A CN102663256B (zh) | 2012-04-11 | 2012-04-11 | 确定暴雨衰减指数为函数型的暴雨强度公式的一种方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210114694.2A CN102663256B (zh) | 2012-04-11 | 2012-04-11 | 确定暴雨衰减指数为函数型的暴雨强度公式的一种方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102663256A true CN102663256A (zh) | 2012-09-12 |
CN102663256B CN102663256B (zh) | 2014-06-04 |
Family
ID=46772746
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210114694.2A Expired - Fee Related CN102663256B (zh) | 2012-04-11 | 2012-04-11 | 确定暴雨衰减指数为函数型的暴雨强度公式的一种方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102663256B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103902828A (zh) * | 2014-04-02 | 2014-07-02 | 北京工业大学 | 城市24小时长历时暴雨强度的确定方法 |
CN104156559A (zh) * | 2014-07-08 | 2014-11-19 | 江苏建筑职业技术学院 | 确定含有重现期修正系数的暴雨强度公式的一种方法 |
CN111949928A (zh) * | 2020-08-10 | 2020-11-17 | 河海大学 | 一种基于暴雨衰减特性的暴雨强度公式确定方法 |
CN113111295A (zh) * | 2021-04-19 | 2021-07-13 | 黄河勘测规划设计研究院有限公司 | 一种淤地坝设计洪峰时暴雨递减指数的计算方法 |
CN113821759A (zh) * | 2021-09-27 | 2021-12-21 | 河海大学 | 一种基于泰勒级数展开的暴雨强度公式参数率定方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1790345A (zh) * | 2005-12-20 | 2006-06-21 | 沈阳建筑大学 | 城市雨水管道的设计方法 |
-
2012
- 2012-04-11 CN CN201210114694.2A patent/CN102663256B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1790345A (zh) * | 2005-12-20 | 2006-06-21 | 沈阳建筑大学 | 城市雨水管道的设计方法 |
Non-Patent Citations (2)
Title |
---|
卢金锁等: "西安市暴雨强度公式的推求研究", 《中国给水排水》 * |
王敏等: "北京城市暴雨和雨型的研究", 《水文》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103902828A (zh) * | 2014-04-02 | 2014-07-02 | 北京工业大学 | 城市24小时长历时暴雨强度的确定方法 |
CN104156559A (zh) * | 2014-07-08 | 2014-11-19 | 江苏建筑职业技术学院 | 确定含有重现期修正系数的暴雨强度公式的一种方法 |
CN104156559B (zh) * | 2014-07-08 | 2016-02-24 | 江苏建筑职业技术学院 | 确定含有重现期修正系数的暴雨强度公式的一种方法 |
CN111949928A (zh) * | 2020-08-10 | 2020-11-17 | 河海大学 | 一种基于暴雨衰减特性的暴雨强度公式确定方法 |
CN113111295A (zh) * | 2021-04-19 | 2021-07-13 | 黄河勘测规划设计研究院有限公司 | 一种淤地坝设计洪峰时暴雨递减指数的计算方法 |
CN113111295B (zh) * | 2021-04-19 | 2022-06-17 | 黄河勘测规划设计研究院有限公司 | 一种淤地坝设计洪峰时暴雨递减指数的计算方法 |
CN113821759A (zh) * | 2021-09-27 | 2021-12-21 | 河海大学 | 一种基于泰勒级数展开的暴雨强度公式参数率定方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102663256B (zh) | 2014-06-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102663256B (zh) | 确定暴雨衰减指数为函数型的暴雨强度公式的一种方法 | |
Wang et al. | Investigation on the influencing factors of energy consumption and thermal comfort for a passive solar house with water thermal storage wall | |
Radhi et al. | Global warming implications of facade parameters: A life cycle assessment of residential buildings in Bahrain | |
CN106875060B (zh) | 一种基于全局自动优化算法的洪水实时校正方法 | |
Fan et al. | Quantifying the magnitude of the impact of climate change and human activity on runoff decline in Mian River Basin, China | |
CN107085658B (zh) | 一种山洪灾害成灾时间确定方法 | |
CN104915575B (zh) | 一种基于层次物元可拓法的海绵城区生态指数的评价方法 | |
CN104345131B (zh) | 一种稻田氮磷径流流失负荷的田间尺度估算方法 | |
GUO et al. | Prediction of future runoff change based on Budyko hypothesis in Yangtze River basin | |
CN107761528A (zh) | 一种基于实时高温天气和交通量的沥青路面车辙预报方法 | |
CN105046051A (zh) | 一种实际环境冻融作用次数预估方法 | |
CN104156559B (zh) | 确定含有重现期修正系数的暴雨强度公式的一种方法 | |
Cui et al. | Effect of length of the observed dataset on the calibration of a distributed hydrological model | |
CN104732110B (zh) | 基于混合weibull分布的森林生物量多尺度估测方法 | |
Thein | Modelling of short duration rainfall IDF equation for Sagaing region, Myanmar | |
CN103473434A (zh) | 城市设计暴雨计算方法 | |
Vaisi et al. | Towards urban energy density (UED) definition | |
Huo et al. | Effectiveness evaluation of LIDs through SWMM: A case study of typical urban unit in Handan, China | |
CN113762729A (zh) | 一种城市排水系统改造工程综合评价方法 | |
Yoo et al. | Estimation of areal reduction factors using a mixed gamma distribution | |
Luo et al. | Using SWAT to simulate runoff under different land use scenarios in Xiangjiang River Basin | |
Lave et al. | Characterizing local high-frequency solar variability for use in distribution studies | |
Liang | Modeling of Modern Eco-Environment Landscape Planning in Sponge City. | |
CN204188166U (zh) | 一种用于长距离测量可调节觇板 | |
CN113821759B (zh) | 一种基于泰勒级数展开的暴雨强度公式参数率定方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140604 Termination date: 20160411 |
|
CF01 | Termination of patent right due to non-payment of annual fee |