CN104156559A - 确定含有重现期修正系数的暴雨强度公式的一种方法 - Google Patents
确定含有重现期修正系数的暴雨强度公式的一种方法 Download PDFInfo
- Publication number
- CN104156559A CN104156559A CN201410330737.XA CN201410330737A CN104156559A CN 104156559 A CN104156559 A CN 104156559A CN 201410330737 A CN201410330737 A CN 201410330737A CN 104156559 A CN104156559 A CN 104156559A
- Authority
- CN
- China
- Prior art keywords
- formula
- parameter
- intensity
- rain
- reoccurrence 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.)
- Granted
Links
Landscapes
- Complex Calculations (AREA)
Abstract
本发明确定含有重现期修正系数的暴雨强度公式的一种方法,对于选配4个参数的暴雨强度公式i=(A1+C1lgP)/(t+b)n,拟合重现期P、降雨历时t、雨强i关系值效果不理想的情况,提出了判断暴雨强度公式中是否需引入重现期修正系数d的方法以及采用高斯-牛顿法或麦夸尔特法确定含有该系数的暴雨强度公式i=[A1+C1lg(P+d)]/(t+b)n中参数的方法。结合两个实施例论述了具体实施步骤。该发明可使求得的暴雨强度公式更能客观地反映当地的暴雨特性、所求暴雨强度公式的参数为最佳拟合参数、使拟合精度显著提高。该公式是城市排水工程设计中确定雨水设计流量的基本依据。
Description
一、技术领域
本发明属于市政工程技术领域。暴雨强度公式是城市排水工程设计中确定雨水设计流量的基本依据之一,该公式的合理性与精度直接影响城市排水工程的规模。
二、背景技术
(一)常用的暴雨强度公式的形式
目前,我国城市常用的暴雨强度公式为[1]:
单一重现期的分公式
各重现期的总公式
式中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为暴雨衰减指数;C1=A1×C。
单一重现期的分公式对降雨历时t,雨强i的关系值(t,i)的拟合精度较高,但只能应用于一个重现期的情况,故目前我国常用式(2)形式的公式。
由于公式(2)为纯非线性回归方程,可采用高斯-牛顿法或麦夸尔特法求解式(2)的拟合参数。
(二)高斯-牛顿法与麦夸尔特法
1.高斯牛顿法
设一般的非线性方程为
y=f(X,θ)+ε
式中f为一般函数;X可以是单个自变量,也可以是r个自变量X=(x1,x2,…,xr);θ为P维参数向量,即θ=(θ1,θ2,…,θp)′;ε为随机误差项,且∈N(0,σ2)。设y和X具有m组观察值(x1j,x2j,…,xrj;yj),j=1~m。求非线性模型参数的“最小二乘”估计,就是求θ的待估计值使残差平方和
为最小。
为方便起见,用记号fj(θ)代替f(Xj,θ),高斯-牛顿法求解满足式(3)Q(θ)为最小的参数递推公式,并写成矩阵形式为[2,3]
θ(k+1)=θ(k)+[J′(θ(k))J(θ(k))]-1J′(θ(k))[(y-f(θ(k))] (4)
式中k为递推次数;y=(y1,y2,…,yj…,ym)′,yj为因变量的第j个观察值;f(θ(k))=[f1(θ(k)),f2(θ(k)),…,fj(θ(k)),…,fm(θ(k))]′,fj(θ(k))为由非线性方程及第k次迭代参数计算的第j个因变量y的估计值,j=1~m。
利用式(4)从参数初值θ(0)开始,一步步递推下去,直到θ(k)收敛稳定,即|θ(k+1)-θ(k)|小于或等于预先指定的小正数δ(例如δ=0.0001),从而得到θ的估计值。
2.麦夸尔特法
麦夸尔特法,放宽了迭代计算时对参数初值的要求,其求解式(3)中暴雨公式参数θ的递推公式,写成矩阵形式为[4,5]
θ(k+1)=θ(k)+[J′(θ(k))J(θ(k))+λI]-1J′(θ(k))[(y-f(θ(k))] (6)
式中:λ为阻尼因子;I为P阶单位阵;其他符号含义与式(4)中相应符号含义相同;偏导数矩阵J(θ(k))的计算式见式(5)。
关于阻尼因子λ的确定方法详见文献[5]。利用式(6)从参数初值θ(0)开始,一步步递推下去,直到θ(k)收敛稳定,即|θ(k+1)-θ(k)|小于或等于预先指定的小正数δ(例如δ=0.0001),从而得到θ的估计值。
(三)暴雨强度公式的推求路径及其拟合误差
推求暴雨强度公式中的参数,有两种不同路径。
路径一:利用理论频率计算所得的雨强选配公式,称其为传统方法。该法涉及两个环节[1,6]。(1)理论频率曲线拟合:首先选配理论频率分布模型,例如,皮尔逊III型分布曲线、耿贝尔分布曲线和指数分布曲线[6],并对各历时相应的实测雨强样本系列进行频率曲线拟合;然后,根据各历时t的理论频率曲线,确定各计算重现期P的设计雨强ip,进而得计算重现期P、历时t、设计雨强ip关系表。此环节存在理论频率曲线与实测雨强样本的拟合误差。(2)暴雨公式参数拟合:根据计算重现期、历时、设计雨强关系表,选配暴雨强度公式。此环节存在由暴雨强度公式拟合由理论频率曲线计算的设计雨强的误差,对于此误差,GB50014-2006《室外排水设计规范》[6](2014版)附录A指出,对于采用年最大值法取样时,“计算重现期在2年~20年时,在一般强度的地方,平均绝对均方差不宜大于0.05mm/min(注:原文为“平均绝对方差不宜大于0.05mm/min”。显然是错误的,均方差才与暴雨强度具有相同量纲)。在较大强度的地方,平均相对均方差不宜大于5%”。
路径二:直接拟合法[4]。文献[4]采用非线性回归方法,根据重现期,历时,实测雨强的关系值(Pj,tj;ij),j=1~m,直接拟合暴雨强度公式参数,该法可实现拟合实测雨强样本的误差平方和为最小,求得的暴雨强度公式参数为整体上的最佳参数,该法不仅减少了理论频率分布模型的选择与拟合环节,从而避免了选择何种理论频率分布模型的纠结及避免推求暴雨强度总公式参数的任意性,而且拟合实测雨强的均方差小于传统方法的结果,提高了拟合精度。
需要指出,实测雨强样本是暴雨强度公式选配与拟合的基础,应与其拟合最佳。直接拟合法的拟合误差则是拟合实测雨强样本的误差,其含义不同于规范[6]中规定的误差。因此,直接拟合法的误差不应该按规范[6]中规定的误差上限来衡量。
三、发明内容
一些地区的计算成果表明,雨力A(即t+b=1时的雨强)与重现期P的对数lgP不是直线关系[1],即A≠A1+C1lgP,此种情况,即使采用非线性最小二乘法或其他最优化方法求解式(2)的拟合参数,所得暴雨强度公式拟合某些重现期的雨强的绝对均方差或平均绝对均方差(即标准差)仍超过0.05mm/min[1,7,8]。有鉴于此,为了更密切的吻合当地多年统计的降雨强度子样点的分布规律,提高拟合精度,不应再拘泥于4参数公式,而应采用式(7)5参数的公式
式中:d为重现期修正系数,简称修正系数;其他符号含义与式(2)中相应符号含义相同。
文献[1]给出了确定式(7)中参数的手工算法(第932-939页)。到目前为止,尚无采用解析法确定式(7)中参数的方法。
为了克服手工算法确定式(7)中参数效率低、所求参数具有一定的任意性、拟合精度不如解析法高的缺陷,本发明提出考虑修正系数d的判别方法及含有该系数的暴雨强度公式(7)中参数的确定方法。
本发明采用的技术方案是:对于4参数暴雨强度公式(2)拟合重现期P、历时t以及由理论频率曲线确定的计算重现期相应的设计雨强ip(或实测雨强i)关系值效果不理想的情况,首先,根据求得的式(2)的参数n,b,利用最小二乘法计算各重现期的雨力A;并点绘A与lgP关系点,若不为直线趋势,即A≠A1+C1lgP,则要考虑修正系数d,选配暴雨强度公式(7);然后,采用非线性回归方法确定包含修正系数d的暴雨强度公式中的5个参数。
本项发明的有益效果是,当A与lgP关系不为直线趋势时,采用含有修正系数d的5参数暴雨强度公式更能客观地反映当地的暴雨特性;采用高斯-牛顿法或麦夸尔特法使所求暴雨强度公式的参数为最佳拟合参数,且计算效率高;拟合精度显著提高。本项发明已利用我国有关城市的暴雨特征值资料进行了验证,效果很好。
四、附图说明
下面结合附图和实施例对本发明进一步说明。
附图1为实例1的A与lgT关系图。
附图2为实例1修正系数初值d(0)=-0.6时A与lg(T-0.6)关系图。
附图3为实例2的A与lgT关系图。
附图4为实例2修正系数初值d(0)=-0.98时A与lg(T-0.98)关系图。
五、具体实施方式
(一)路径一——利用理论频率计算所得的雨强选配公式
实例1:文献[8]中昆明市松华坝雨量站,具有1964~2000年年最大10、20、30、45、60、90、120、180min暴雨系列,将各时段样本系列采用皮尔逊III型分布频率计算后,得到重现期P、降雨历时t、设计暴雨强度ip的72组样本点(Pj,tj;ip,j),j=1~72,见表1。
表1实例1重现期历时设计雨强P-t-ip关系值
利用本发明确定暴雨强度公式,其具体操作如下。
1.根据表1数据,采用文献[2,4]方法,推求4参数暴雨强度公式:
并计算由式(8)拟合表1数据的误差平方和Q=0.2078。
进一步计算由式(8)拟合表1中各重现期相应的不同历时雨强的均方差σ:
式中ip,j为根据计算重现期由理论频率曲线得到的雨强,mm/min,本例即为表1中某一重现期相应的不同历时的雨强;ig,j为根据重现期、降雨历时由暴雨强度公式计算的雨强,mm/min;m1为降雨历时个数,本例m1=8。
计算结果见表2。
表2实例1四参数公式拟合雨强的均方差
可见,4参数暴雨强度公式对于重现期较大或较小的情况拟合效果不理想,拟合均方差超过0.05mm/min,需判断暴雨公式中是否需要考虑修正系数d。
2.A与lgT关系趋势的判别方法
对于4参数暴雨强度公式中的参数n,b一定时,本发明计算某一重现期P相应的雨力A的方法是:
根据某一重现期P的m1组实测关系点(tj,ip,j),j=1~m1,基于式(1),求满足
为最小的A值。即对A求偏导数:
整理得:
因此,利用式(12)即可计算各重现期P相应的雨力A,进而点绘(A,lgP)关系点,若不为直线趋势,则要考虑修正系数d,使A=A1+C1lg(P+d),进而选配5参数暴雨强度公式(7),否则,可通过其他方法,提高拟合精度。
对实例1,根据表1数据、式(8)中参数n=0.7342,b=11.4082,利用式(12)计算各重现期P相应的雨力A,并点绘(A,lgP)关系点,如图1,可见关系不为直线趋势,则要考虑修正系数d。
3.含有重现期修正系数d的暴雨强度公式参数的确定方法
(1)对式(7)中参数A1,C1,d,b,n分别求偏导数
(2)确定暴雨强度公式(7)中参数的初值θ(0)=(A1(0),C1(0),d(0),b(0),n(0))′
将求得的4参数暴雨强度公式中的参数值作为初值A1(0),C1(0),b(0),n(0);根据A与lgT关系点趋势,在d>-Pmin(Pmin为表1样本点中的最小重现期)范围初拟修正系数d(0),使A与lg(P+d(0))近似为直线趋势。本发明实例1中初拟d(0)=-0.6,A=A1+C1(lgP-0.6)近似为直线关系,如图2所示。故本发明实例1参数的初值θ(0)=(9.8857,9.5129,-0.6,11.4082,0.7342)′
(3)由参数的初值θ(0)、式(13)~式(17)以及各个重现期的M组观察值(Pj,tj;ip,j),j=1~M及式(5)计算偏导数矩阵J(θ(0))及矩阵f(θ(0))。
(4)根据式(4)或式(6)计算θ(1)。
(5)再以θ(1)作为参数的初始值θ(0),重复步骤(3),步骤(4),根据给定的δ值(例如,δ=0.0001),经若干次递推迭代,可求得公式(7)中参数的估计值θ=(A1,C1,d,b,n)′。
本发明实例1的计算结果:θ=(A1,C1,d,b,n)′=(12.3965,7.9181,-0.7194,11.4300,0.7347)′,则
由式(18)拟合表1雨强的误差平方和Q=0.0309。进一步由式(18)计算拟合表1中各重现期不同历时雨强的均方差σ,结果见表3。
表3实例1五参数公式拟合雨强的均方差
对实例1,比较4参数、含有修正系数d的5参数暴雨强度公式的拟合表1全部雨强数据的误差平方和Q及拟合各重现期雨强的均方差(表2、表3),可见,考虑修正系数d的5参数暴雨强度公式显著提高了拟合精度。
(二)路径二——直接拟合法
根据暴雨资料选样得到的重现期P、实测暴雨强度i和降雨历时t的m组关系值(Pj,tj;ij),j=1~m,即为直接拟合法确定暴雨强度公式(7)中拟合参数的样本点。下面结合实例2说明采用直接拟合法确定含有修正系数d的暴雨强度公式的方法。
实例2:A城市雨量站为国家基本气象站,利用该站1980~2011年的短历时自记雨量计观测的雨量资料,采用年最大值选样法,分别得到降雨历时5min、10min、15min、20min、30min、45min、60min、90min、120min、180min的10个历时的实测雨强样本,对其进行经验频率及重现期计算后,得重现期历时实测雨强P-t-i关系表,如表4所示,共320组数据(限于篇幅,仅列出部分数据)。
表4实例2重现期历时实测雨强P-t-i关系值
1.根据表4数据,采用文献[2,4]方法,推求4参数暴雨强度公式:
并计算由式(19)拟合表4中全部实测雨强的误差平方和为Q=5.6316。进一步计算由式(19)拟合表4中各重现期不同历时相应雨强的均方差σ,结果见表5。
表5实例2直接法四参数暴雨强度公式的拟合均方差(单位:mm/min)
由表5可见,4参数暴雨强度公式对于重现期较大或较小的情况拟合效果不理想,需判断暴雨公式中是否需要考虑修正系数d。
2.A与lgT关系趋势的判别方法
根据式(19)中b=14.5130,n=0.7921及表4中各重现期33.00a,16.50a,…,1.03a的雨强,利用式(12)计算各重现期的雨力A,并绘制A与lgP关系点如图3所示,关系点不是直线趋势,而应考虑修正系数d,采用5参数暴雨强度公式(7)。
3.确定式(7)中参数的初值θ(0)=(A1(0),C1(0),d(0),b(0),n(0))′
将式(19)中的各参数值分别作为初值A1(0),C1(0),b(0),n(0);在d>-Pmin(Pmin为表4样本点中的最小重现期)范围初拟修正系数d(0)=-0.98,绘制A与lg(P-0.98)关系点如图4所示,关系点近似为直线趋势。故实例2暴雨强度公式参数的初值θ(0)=(18.4451,18.1057,-0.98,14.5130,0.7921)′。
4.利用式(4)或式(6)确定含有修正系数d的5参数暴雨强度公式,具体方法与路径一相应环节的计算方法类似,不再赘述。对实例2,得
由式(20)拟合表4中全部实测雨强的误差平方和为Q=1.6593。进一步计算由式(20)拟合各重现期不同历时实测雨强的均方差及全部重现期的均方差,结果见表6。
强调指出,表6为直接拟合各重现期实测雨强的均方差,不同于传统方法拟合由理论频率曲线计算雨强的误差,故不能用规范[6]中的误差限制值0.05mm/min来衡量。
表6实例2五参数暴雨强度公式的拟合均方差(单位:mm/min)
5.实例2不同公式拟合效果的比较。4参数公式(19)与5参数公式(20)拟合实测雨强样本的误差平方和分别为5.6316、1.6593;比较表5、表6拟合各重现期实测雨强的均方差,在32个重现期中,5参数公式有28个重现期的拟合均方差均显著小于4参数公式的结果,仅有4个重现期的拟合均方差与4参数公式的结果相等;5参数公式拟合全部重现期实测雨强的均方差也显著小于4参数的拟合结果。故5参数公式显著提高拟合精度,特别是能显著降低重现期较小与较大时的拟合均方差。
有关文献
[1]北京市市政工程设计研究总院.给水排水设计手册,第5册,城镇排水,第二版[M].北京:中国建筑工业出版社,2002.932-939.
[2]张子贤.用高斯-牛顿法确定暴雨公式参数[J].河海大学学报,1995,23(5):106-111.
[3]方开泰.实用回归分析[M].北京:科学出版社,1988.168-172.
[4]张子贤,孙光东,孙建印,等.城市暴雨强度公式拟合方法的研究[J].水利学报,2013,44(11):1263-1271.
[5]袁志发,周静芋.多元统计分析[M].北京:科学出版社,2002.99-102.
[6]GB50014-2006.室外排水设计规范(2014版)[s].
[7]樊建军,王峰,陈鹏飞.利用MATLAB推导城市暴雨强度公式[J].中国给水排水,2010,26(11):114-115.
[8]柏绍光,黄英,方绍东,等.城市暴雨总公式与分公式精度分析[J].人民长江,2006,37(8):34-35.
Claims (2)
1.确定含有重现期修正系数的暴雨强度公式的一种方法,根据不同重现期P、各历时t以及由理论频率曲线所得的设计暴雨强度或实测暴雨强度i的关系值(统称为重现期降雨历时暴雨强度关系点),确定各重现期的总公式
式中i为暴雨强度,简称雨强,mm/min;P为重现期,a;t为降雨历时,min。而A1,C,b,n,d是与地方暴雨特性有关的参数,其中A1为重现期为1a的设计降雨的雨力,mm/min;C为雨力变动参数,b为降雨历时修正参数,min;n为暴雨衰减指数;d为重现期修正系数,简称修正系数;而C1=A1×C。
其特征在于,具体操作步骤为:
步骤1,根据不同重现期P、各历时t以及由理论频率曲线所得的设计暴雨强度或根据不同重现期P、各历时t以及实测暴雨强度,推求4参数暴雨强度公式
并计算式(2)拟合样本点的均方差,若拟合均方差较大(例如,拟合由理论频率曲线所得的设计暴雨强度时,某些重现期的绝对均方差或平均绝对均方差超过0.05mm/min),则需判断暴雨强度公式中是否引入修正系数d,其特征在于,具体操作步骤为:
(1)根据某一重现期P的m1组关系点(tj,ij),j=1~m1及式(2)中的参数n,b,利用式(3)计算各重现期的雨力A
(2)点绘A与lgP关系点,若不为直线趋势,即A≠A1+C1lgP,则要引入修正系数d,选配暴雨强度公式(1)。
步骤2,确定含有修正系数d的暴雨强度公式(1)中的参数:
步骤2.1,确定公式(1)中参数的初值θ(0)=(A1(0),C1(0),d(0),b(0),n(0))′;
步骤2.2,根据不同重现期P、各历时t以及由理论频率曲线所得的设计暴雨强度关系值,或根据不同重现期P、各历时t以及实测暴雨强度关系值,利用高斯-牛顿法及其公式(4)
θ(k+1)=θ(k)+[J′(θ(k))J(θ(k))]-1J′(θ(k))[(y-f(θ(k))] (4)
确定含有修正系数d的暴雨强度公式(1)中的参数。
或采用麦夸尔特法及其公式(5)
θ(k+1)=θ(k)+[J′(θ(k))J(θ(k))+λI]-1J′(θ(k))[(y-f(θ(k))] (5)
确定含有修正系数d的暴雨强度公式(1)中的参数。
2.根据权利1中所述确定含有修正系数d的暴雨强度公式参数的一种方法,其特征在于,步骤2.1中,将求得的4参数暴雨强度公式中的参数值作为初值A1(0),C1(0),b(0),n(0);根据A与lgT关系点趋势及重现期降雨历时暴雨强度关系点中的最小重现期Pmin,在d>-Pmin,范围拟定修正系数的初值d(0)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410330737.XA CN104156559B (zh) | 2014-07-08 | 2014-07-08 | 确定含有重现期修正系数的暴雨强度公式的一种方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410330737.XA CN104156559B (zh) | 2014-07-08 | 2014-07-08 | 确定含有重现期修正系数的暴雨强度公式的一种方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104156559A true CN104156559A (zh) | 2014-11-19 |
CN104156559B CN104156559B (zh) | 2016-02-24 |
Family
ID=51882057
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410330737.XA Expired - Fee Related CN104156559B (zh) | 2014-07-08 | 2014-07-08 | 确定含有重现期修正系数的暴雨强度公式的一种方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104156559B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106250669A (zh) * | 2016-07-07 | 2016-12-21 | 河海大学 | 一种干旱重现期计算中确定干旱阈值的方法 |
CN106529140A (zh) * | 2016-10-26 | 2017-03-22 | 沈阳建筑大学 | 一种暴雨强度公式确定方法 |
CN113821759A (zh) * | 2021-09-27 | 2021-12-21 | 河海大学 | 一种基于泰勒级数展开的暴雨强度公式参数率定方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1790346A (zh) * | 2005-12-20 | 2006-06-21 | 沈阳建筑大学 | 一种城市雨水管道的设计方法 |
CN102663256A (zh) * | 2012-04-11 | 2012-09-12 | 江苏建筑职业技术学院 | 确定暴雨衰减指数为函数型的暴雨强度公式的一种方法 |
-
2014
- 2014-07-08 CN CN201410330737.XA patent/CN104156559B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1790346A (zh) * | 2005-12-20 | 2006-06-21 | 沈阳建筑大学 | 一种城市雨水管道的设计方法 |
CN102663256A (zh) * | 2012-04-11 | 2012-09-12 | 江苏建筑职业技术学院 | 确定暴雨衰减指数为函数型的暴雨强度公式的一种方法 |
Non-Patent Citations (2)
Title |
---|
范泽华: "天津市降雨趋势分析及设计", 《中国优秀硕士学位论文全文数据库》, 31 August 2012 (2012-08-31) * |
黄津辉: "天津市设计暴雨方法比较及公式修正", 《天津大学学报(自然科学与工程技术版)》, 30 April 2013 (2013-04-30) * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106250669A (zh) * | 2016-07-07 | 2016-12-21 | 河海大学 | 一种干旱重现期计算中确定干旱阈值的方法 |
CN106250669B (zh) * | 2016-07-07 | 2019-07-12 | 河海大学 | 一种干旱重现期计算中确定干旱阈值的方法 |
CN106529140A (zh) * | 2016-10-26 | 2017-03-22 | 沈阳建筑大学 | 一种暴雨强度公式确定方法 |
CN106529140B (zh) * | 2016-10-26 | 2019-02-19 | 沈阳建筑大学 | 一种暴雨强度公式确定方法 |
CN113821759A (zh) * | 2021-09-27 | 2021-12-21 | 河海大学 | 一种基于泰勒级数展开的暴雨强度公式参数率定方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104156559B (zh) | 2016-02-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Yao et al. | Renewable energy, carbon emission and economic growth: A revised environmental Kuznets Curve perspective | |
Bouchouicha et al. | Estimating the global solar irradiation and optimizing the error estimates under Algerian desert climate | |
Chen et al. | Uncertainty of downscaling method in quantifying the impact of climate change on hydrology | |
Yao et al. | New decomposition models to estimate hourly global solar radiation from the daily value | |
Bakirci | Correlations for estimation of daily global solar radiation with hours of bright sunshine in Turkey | |
Sun et al. | Assessment of surface water resources and evapotranspiration in the Haihe River basin of China using SWAT model | |
CN106600060B (zh) | 基于相似日晴朗系数修正的太阳辐照量预测方法 | |
CN104156559B (zh) | 确定含有重现期修正系数的暴雨强度公式的一种方法 | |
CN102663256B (zh) | 确定暴雨衰减指数为函数型的暴雨强度公式的一种方法 | |
CN102945508A (zh) | 一种基于模型校正的风电功率预测预报系统及方法 | |
CN103902828A (zh) | 城市24小时长历时暴雨强度的确定方法 | |
CN103017668B (zh) | 基于高光谱遥感反射率的海冰厚度计算方法 | |
CN102305792B (zh) | 基于非线性偏最小二乘优化模型的森林碳汇遥感估算方法 | |
CN111984700A (zh) | 一种基于水热耦合平衡原理的月径流预测方法 | |
CN103257000B (zh) | 一种桥梁结构日照作用分析的极值温度预测方法 | |
CN105809271A (zh) | 一种基于组合预测法的生物量模型估计方法 | |
Wei et al. | Nierji reservoir flood forecasting based on a Data-Based Mechanistic methodology | |
CN103530521B (zh) | 基于傅立叶级数和arma模型的日照温度时程模拟方法 | |
CN103473434A (zh) | 城市设计暴雨计算方法 | |
CN104732110B (zh) | 基于混合weibull分布的森林生物量多尺度估测方法 | |
Liu et al. | Study of rainfall prediction model based on GM (1, 1)-Markov chain | |
庄稼成 et al. | Attribution analysis of runoff change based on the abcd model coupled with the snowmelt module in the source region of the Yellow River | |
CN103048056B (zh) | 一种日照温差采集样本概率密度的测定方法 | |
CN104598982A (zh) | 年最大降温负荷的获取方法及系统 | |
CN113762729A (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160224 Termination date: 20180708 |