乘潮潮位的理论分布探讨.docVIP

乘潮潮位的（中国海洋大学工程学院, 山东 青岛 266100） 摘要：我国《海港水文规范》在确定港口工程设计潮位时，采用手绘累积频率曲线，再摘取设计潮位值的方法，由于缺乏理论曲线，确定的数值存在任意性。本文将改进最大熵分布应用于海岸港的设计高、低潮位的推算，计算精度高，克服了传统方法推算设计潮位值的不足，有重要理论与工程意义。 关键词：；分布； 在港口工程中，当考虑乘潮作业时（如船舶乘高潮进出港口、乘低潮施工），需要统计确定乘潮潮位的累积频率，其步骤如下： （1）当考虑船舶进出港口时，首先须确定乘潮所需持续时间t。 （1） 式中，Lk表示航道长度，v为航道内允许航速；k为时间富裕系数，可取1.1~1.3。 （2）在潮位过程线上，量取各次潮峰上的历时等于t的潮位值，统计其在不同潮位级的出现次数。 （3）一般采用10厘米为一级，对潮位分级。 （4）由高到低逐级进行累积出现次数的统计。 （5）计算各潮位级的累积频率p，公式如下： （2） 式中，i为各潮位级相应的累积出现次数；n为高潮位总次数。 （6）在方格纸上，以纵坐标表示潮位，以横坐标表示累积频率，把各累积频率值点于相应潮位级的下限处，连绘成高潮位累积频率曲线。 （7）读取累积频率为75%的潮位，即为所需的高潮乘潮潮位。 如系计算低潮乘潮历时t为2小时，累积频率为80%的低潮乘潮施工水位，则应将第2步改为“在潮位过程线上，量取各次潮谷上的历时等于2小时的潮谷值”，同时，在计算低潮乘潮作业水位时，系指低于该潮位所占的累积频率，其余统计步骤与高潮乘潮水位的相同。 上述方法在绘制累积频率曲线时，缺乏理论曲线，的潮位值存在任意性。本文将最大熵分布 1 最大熵分布 1975 年Ulrych 等将最大熵原理应用在海浪的谱分析中[2]；1996 年，吴克俭等将该原理应用于海浪波高的统计分布，指出Rayleigh分布为理想状况下，波高达到最大混乱程度时必然遵循的一种分布，但若用于描述实际波浪则有偏差，更提出Weibull分布为某些约束条件下的最大熵分布的特殊形式[3]；徐福敏于2000年推导出二阶矩意义下的最大熵概率密度分布表达式，并由最大熵法得到了深水充分发展波波高概率分布及风浪槽试验中不同水深时的波高概率分布[4]； Xu Delun推导出非线性海浪波高的最大熵分布，并将其应用于多年一遇海洋要素估计[]；周良明等指出了波高熵与波浪的状态参量及平均波高的关系，将最大熵分布与广泛应用的瑞利分布作了比较并指出最大熵分布具有明显的优越性[]。[7]。董胜等提出了泊松最大熵分布，并将其用于计算重现值波高[]。 在推导最大熵分布的过程中，Xu Delun 等提出适用于海洋环境中的约束条件，并推导得到了适用于描述海洋工程环境条件的最大熵分布密度函数[]。约束条件如下： （） 代入Euler方程，并令，，，可得最大熵分布的概率密度函数为[]： （） 为了拟合的资料，在此引入位置参数a0，由式（）可得改进的最大熵分布 （） 差系数Cv偏态系数CsCv与Cs这2个参数，直至观测点均匀分布在理论分布曲线的两侧。 为了估计改进最大熵分布的参数，本文选取理论分布频率与经验频率的离差平方和最小为，来确定拟合观测值的最佳分布函数。 若设计潮位为随机变量，对任一统计序列，假设数学期望、二阶中心矩及三阶中心矩存在，则它们的表达式如下： （） 将式（）代入式（a），可得 （） 将式（）和（）代入，则 （） 将式（）代入式（）可得 （） 联立式（a）与（b）、式（a）与（c），消去可得 （8） 改进最大熵分布拟合首先由样本求出Cv值，选取不同的k得到Cs=k(Cv，改变k和a0，限定(在某范围内有意义，可以得出对应(的值，进而求出(和(，得到不同的参数组对应的改进最大熵分布密度函数，计算其估计概率值与经验概率值的平均离差平方和Q，以使离差平方和Q最小为标准所求得的参数即为所求的参数，相应的分布函数即为所求分布函数。 2 工程应用 按《海港水文规范》规定，确定高（低），应有1年的潮位资料 选取威海湾1992年的验潮资料，进行了乘高、低潮2小时频率分析，手工绘制求得了乘2小时高潮位值计算结果如表1。 表1 威海湾2小时乘高潮潮位 Table 1 Ridable high tidal level at Weihai Bay 计算方法 累积频率 75% 80% 85% 90% 手工绘制方法/cm 172 167 163 156 最大熵分布法/cm 172 167 161 154 采用最大熵分布，对高潮数据进行拟合，其如图1，累积分布曲线如图，计算得出潮位值此结果与手绘曲线得到的结果cm。 图1 高潮潮