施工导流超标洪水风险率计算的最大熵法.pdfVIP

第４１卷 第３期 人　民　长　江 Ｖｏｌ．４１，Ｎｏ．３ ２０１０年 ２月　 　 Ｙａｎｇｔｚｅ　Ｒｉｖｅｒ　 　 Ｆｅｂ．， ２０１０ 　　文章编号：１００１－４１７９（２０１０）０３－００７５－０３ 施工导流超标洪水风险率计算的最大熵法 李 旭 东，郭 红 民 （三峡大学 土木水电学院，湖北 宜昌４４３００２） 摘要：利用最大熵原理推导了超标洪水间隔时间的熵概率密度函数模型，在此基础上给出了施工导流超标洪 水风险率的计算方法。方法可直接根据样本的数字特征值推断洪水类型，而不是先假设典型的概率分布，因 而减少了主观假定，具有更充分的数学和物理意义。利用ＭＡＴＬＡＢ软件对具体工程实例进行了计算，并与传 统假定为威布尔分布的方法进行了对比分析。结果表明：给出的计算方法精度较高，且算法简单、便于应用。 关　键　词：施工导流；超标洪水；风险率；最大熵法 中图法分类号：ＴＶ５５１．１　　　文献标志码：Ａ 　　施工导流系统超标洪水风险率的计算，传统方法 ｉ ｘｆ（ｘ）ｄｘ＝ 　（ｉ＝１，２，…，Ｎ） （３） ∫ μｉ 是首先预估或假定超标洪水间隔时间的分布类型，然 Ｒ 式中，ｆ（ｘ）是连续性随即变量 ｘ的分布密度函数；Ｒ 后通过抽样或根据历史资料对该假定进行拟合优度检 为积分空间； 为第 ｉ阶原点矩，可由计算样本得到； 验，通过后据以计算其风险率。文献［１］即按照上述 μｉ 思路对假设超标洪水间隔时间服从威布尔分布进行了 Ｎ为所用矩的阶数。事实上，随机变量的特性基本上 有关讨论。 可以用它的前４阶矩来描述。因为４阶矩方法已能达 近年来，最大熵理论在结构可靠度分析、岩土工程 到精度要求且其收敛性优于５阶矩法，故本文均采用 反分析及岩石力学参数概率分布等许多方面的应用都 ４阶矩方法。 取得了较好的效果。本文应用最大熵原理，直接根据 通过调整ｆ（ｘ），使熵 Ｈ（ｘ）达到最大值，采用拉 有关样本的数字特征值进行计算，就可以得到精度较 格朗日乘子法来解决此问题。建立如下拉格朗日函 高的概率分布密度函数，进而求出施工导流系统的超 数： 标洪水风险率。 Ｌ＝Ｈ（ｘ）＋（ ＋１）｜ ｆ（ｘ）ｄｘ－１｜＋ λ ０ ∫ Ｒ １　最大熵密度函数 ｎ ｉ ｜ ｘｆ（ｘ）ｄｘ－ ｜　　　 （４） λ μ