streamlit 快速搭建芝加哥雨型应用_暴雨强度公式编制-CSDN博客

本文链接：https://blog.csdn.net/weixin_42984235/article/details/139756063

目的

通过streamlit快速搭建芝加哥雨型应用，该应用实现以下功能：

内置重庆各区县的暴雨强度计算公式，
支持自定义城市的暴雨公式参数，
支持雨型生成，输出图片和文字
提供降雨历时、汇水面积、径流系数输入，支持径流量计算
在线部署，网页使用。

原型

用 obsidian excalidraw 简单绘制原型。

主要原理

《关于发布重庆市暴雨强度修订公式及设计暴雨雨型的通知》（渝建〔2017〕443 号）
芝加哥合成暴雨过程线的公式推导 (vivifree.com)

主要功能实现

雨强计算

def intensity(A, B, C, N, t, P):
    """
    雨强计算。
    参数:
    - a,b,c,n: 参数。
    - p (float): 设计重现期（单位：年）。
    - t (np.ndarray): 分钟数组
    返回:
    - its: 雨强（单位：mm/min）。
    """
    a = A * 0.4 * (1 + C * math.log10(P))
    its = a * ((1 - N) * t + B) / np.power(t + B, N + 1)
    return its
    
q = (A * (1 + C * math.log10(P))) / ((duration_minutes + B) ** N)

雨强分布计算

def rainCalc_single_period(A, B, C, N, T: int, p: float, peak_ratio: float):
    """
    计算单一时段内的降雨强度分布。
    参数:
    - T (int): 降雨持续时间（单位：分钟）。
    - p (float): 设计重现期（单位：年）。
    - peak_ratio (float): 雨强峰值所在时间占总降雨历时的比例。
    返回:
    - np.ndarray: 随时间变化的降雨强度数组（单位：mm/min）。
    内部参数:
    - t (np.ndarray): 分钟数组
    - peak_time (float): 峰值时间
    """
    # ...函数实现代码...
    t = np.arange(0, T)
    peak_time = T * peak_ratio
    itAr = np.zeros(len(t))
    # 计算雨强
    for i in range(len(t)):
        if t[i] < peak_time:
            itAr[i] = intensity(A, B, C, N, (peak_time - t[i]) / peak_ratio, p) / 60
        else:
            itAr[i] = intensity(A, B, C, N, (t[i] - peak_time) / peak_ratio, p) / 60
    return itAr