联合概率分布
文明上网手抄报-《十二月花名歌》
2023年2月22日发(作者：高二周记)随机波浪联合分布概率模型
郄禄文;张立丹;蒋学炼
【摘要】针对传统统计方法难以准确描述波浪联合分布特征的问题,结合波浪物理
意义及相应统计参数,提出了可描述波浪要素联合分布的概率模型.首先给出了波高、
周期和波向的边缘分布及联合分布理论概率模型,然后利用实测数据进行了验证.结
果表明,本模型可有效重现波浪条件.
【期刊名称】《河北大学学报（自然科学版）》
【年(卷),期】2014(034)003
【总页数】6页(P238-243)
【关键词】波高;周期;波向;联合分布
【作者】郄禄文;张立丹;蒋学炼
【作者单位】河北大学建筑工程学院,河北保定071002;河北大学建筑工程学院,河
北保定071002;天津城建大学土木工程系,天津300384
【正文语种】中文
【中图分类】O353.2
波浪力是防波堤最重要的外部荷载.文献[1]通过规定"波列累积频率"来反映不同类
型建筑物上波浪力的作用性质.从表1可看出，因为波列中个别大波很可能造成直
立结构的破坏，且修复困难，因此直立结构采用较小累积频率，而波列中个别大波
并不能决定斜坡结构的毁坏，且较易修复，故可采用较大累积频率.
文献[1]中也提出了波浪长期分布推演方法，即基于波浪观测或气象图推算得到的
波浪数据分波向统计其波高或周期频率，再利用Pearson-III型、对数正态型、极
值I型、Weibull型等理论累积频率分布对经验累积频率曲线拟合外延.
但上述方法获得的波高、周期为离散值，无法清晰表征出两者之间内含的物理关系
(在一定水深处，波高、周期之间存在波陡的约束)，且无法满足风险评估工作中对
波浪要素连续型联合概率分布的要求.因此，本文基于文献[2]和[3]的理论对该问题
进行探讨，以提供便于概率积分求解的波高、周期、波向联合概率分布函数.
表1设计波高的累积频率标准Tab.1Cumulativefrequencystandardvalueof
designwaveheight类型部位项目波高累积频率F/%直墙/墩柱上部结构/墙
身/墩柱/桩基强度/稳定1基床/护底块石稳定5胸墙/堤顶方块强度/稳定1斜坡护
面块石/护面块体稳定13/5护底块石稳定13
进行斜坡式建筑物的护面块石、护面块体稳定性计算时，平均波高与水深之比时，
F取取5%.
1波浪气候综合模式
依照NMIMET模式[2]中，首先分波向对波浪实测数据进行统计分析，获得给定波
向下的波高与周期的联合概率分布F(H,T)=F(T|H)F(H)，然后加入波向要素以获得
三者的联合分布F(H,T,D)=F(H,T|D)F(D).由于实际海况的随机三维性，实时观测的
波向与相应的波高周期数据往往不匹配，具有很大的波动性，考虑到极值统计针对
大浪，其主波向通常与对应风向偏差很小，因此本文采用分风向代替分波向进行统
计分析.需要说明的是，本文提到的波浪周期均为平均跨零周期Tz，以下简称周期
(T).
1.1指定波向的波高极值分布
文献[1]建议利用Pearson-III型、对数正态型、极值I型、Weibull型等理论概型
拟合经验累积频率曲线以获得波高长期分布.但文献[4]认为真实海况为风浪涌浪并
存的混合浪，上述方法存在较大误差，并建议可借助风速-波高参数化模式推导出
波高长期分布，如下.
1)风速-波高联合分布密度
f(H,Wr)=f(H|Wr)f(Wr)，
(1)
其中，H为波高，Wr为特征风速，f(Wr)为特征风速的概率密度，f(H|Wr)为指定
特征风速Wr下波高H的条件概率密度.
公式(1)的建立基于如下前提：
①指定特征风速Wr对应的波高Hr为合成浪(风浪+涌浪)，即
(2)
其中，HY代表Wr=0时的平均涌浪波高，a,n为风成浪的波高参数[5].
②给定风速下的波高标准差
(3)
其中，b,c,d为根据实测数据拟合出的参数[5].
③给定风速下的波高条件分布密度为伽马分布，即
(4)
式中，
由式(4)得到f(H|Wr)后，根据概率学理论，通过积分即可求得波高的边缘分布
F(H)=∬f(H|Wr)f(Wr)dHdWr.
(5)
2)波高长期概率分布
由公式(5)获得波高边缘分布F(H)后，可引入理论概型进行外延拟合.由于理论概型
对波高长期分布影响很大，故应尝试多种概型，择优选用.工程中常用概型如下：
①对数正态分布
(6)
其中，μh和σh分别是lnH的平均值和标准差.
②Weibull分布
(7)
其中，A为位置参数；B为比例尺参数；C为形状参数.
③极值-I型分布
F(H)=exp(-exp(-(H-A)/B))，
(8)
其中，A为位置参数；B为比例尺参数.
④极值-II型分布
(9)
其中，A为位置参数；B为比例尺参数；C为形状参数.
⑤伽马分布(Pearson-III型曲线)
(10)
其中，A为位置参数；B为比例尺参数；C为形状参数.
⑥广义伽马分布
(11)
其中，A为位置参数；B为比例尺参数；C为形状参数，D为指数参数.
外延拟合过程中，适线法、矩法、概率权重矩法、权函数法及模糊数学法等常用参
数估计方法存在计算过程复杂、需人为调整等问题.文献[6]所提出的多变量优化数
值方法更为简便，拟合结果与传统方法基本吻合.
1.2波高与周期的联合分布
迄今为止，已有众多学者对波高-周期长期联合分布进行了研究.其中文献[7]基于波
高-周期散布资料所提出的5参数对数正态模式较为常用
(12)
其中μh，μt，σh，σt，ρht分别为h=lnH和t=lnT的均值、标准差、相关系数，
可通过计算实测样本(H,T)的对数值(lnH,lnT)样本矩得到.
上述模式虽然获得了较好拟合结果[7]，但其要求波高及周期均符合对数正态边缘
分布却使其使用受到了限制.所以，此处仍以风速-波高参数化模式来描述其联合分
布
f(H,T)=f(T|H)f(H).
(13)
与公式(12)相比，公式(13)的优势在于不限定波高的边缘分布概型f(H).
此外，波周期的条件分布概型[4]
(14)
式中，μt|H为在给定波高H时对数周期t=lnT的条件均值；σt|H为给定波高H
时t的条件标准差
其中，σT|H和μT|H为给定波高H下的周期条件均值、标准差，其计算式为
其中，μH,μT,σH,σT,ρHT为H，T的总体平均值、总体标准差和总体相关系数.综
合上述公式，基于观测数据获得μH,μT,σH,σT,ρHT5个统计参数，结合波高边缘
分布F(H)导出波高周期联合分布F(H，T).其中，μH,σH可由波高边缘分布F(H)直
接导出，μT,σT,ρHT则可根据文献[4]的回归公式可得
(波浪不受遮蔽的海域)
(波浪受遮蔽的沿海海域)
ρHT=0.429+0.0957σγ,
其中，μt,σt分别为对数周期t=lnT的均值、标准差.波浪受遮蔽的非沿海海域采用
上述2种情况的加权值.
1.3波高-周期-波向联合分布
上述2节给出了指定风向的波高边缘分布F(H)和周期条件分布F(T|H).基于“同时
刻主/强波向与主/强风向近似一致”的认识，假定指定波向D与指定风向d条件
下的波高-周期的联合概率分布F(H，T)一致，即F(H,T|D,d)≈F(H,T|d)，
(15)
基于概率理论，波高-周期-波向的联合分布概型为F(H,T,D)≈F(H,T|d)F(d).
(16)
2工程实例分析
为验证上述概型的准确性，采用美国国家数据浮标中心(NDBC)46001测站(位于
北纬56°17'44"、西经148°10'19"，水深为4206m，资料跨度为1972—2004
年)的实测数据进行日极值系列分析.由于篇幅限制，此处仅列出N向分析结果，其
余方向具有类似规律.
利用传统的直接统计方法以及本文提出的多变量函数最小值的数值方法，对
46001测站33年波浪(北向)的平均周期和有效波高列于表2和表3,其频数图见图
1和图2.
表246001测站N向波高-周期联合分布(直接统计)Tab.2Jointdistributionof
N-directionwaveaboutwaveheightandperiodatNo.46001station
(Directstatisticalmethod)有效波高H/m平均周期
/s4.815.766.727.678.639.5810.5411.4912.4513.401.72.16.35.73.60.60.20.40.0
0.00.02.80.45.910.67.63.81.90.60.00.00.43.80.00.46.38.95.11.31.10.20.00.04.8
0.00.01.55.33.82.70.80.00.00.05.90.00.00.01.53.41.30.80.60.20.06.90.00.00.00
.00.40.60.80.20.20.07.90.00.00.00.00.20.60.60.20.00.08.90.00.00.00.00.00.40.
00.00.00.010.00.00.00.00.00.00.00.00.00.00.011.00.00.00.00.00.00.00.00.00.0
0.2
表346001测站N向波高-周期联合分布(本文数值方法)Tab.3Jointdistribution
ofN-directionwaveaboutwaveheightandperiodatNo.46001
station(Numericalmethod)有效波高H/m平均周期
/s4.815.766.727.678.639.5810.5411.4912.4513.401.72.26.16.63.91.50.40.10.0
0.00.02.80.85.19.98.64.31.50.40.10.00.03.80.11.66.08.76.22.70.80.20.00.04.80
.00.21.84.75.33.21.20.30.10.05.90.00.00.31.42.72.41.20.40.10.06.90.00.00.00.
20.81.20.90.40.10.07.90.00.00.00.00.10.30.40.20.10.08.90.00.00.00.00.00.10.1
0.10.10.010.00.00.00.00.00.00.00.00.00.00.011.00.00.00.00.00.00.00.00.00.00.
0
图1散布频数图(直接统计)图2散布频数图(本文方法)Fig.1Spread
frequencychart(Directstatisticalmethod)Fig.2Spreadfrequencychart
(Numericalmethod)
3结论
本文引入波浪气候综合模式(NMIMET)[2]对波浪观测资料进行统计分析，引进随
机波浪多变量联合分布函数,提出一种描述波高、周期、波向联合分布的概型，便
于防波堤风险分析与管理过程中的概率积分求解.通过美国国家数据浮标中心
46001测站波浪观测数据分析表明，采用求取多变量函数最小值的数值方法进行
求解，与传统方法波浪年极值数据拟合结果基本吻合，可有效重现波浪条件，但计
算过程更为简便.
参考文献:
[1]中交第一航务工程勘察设计院有限公司.海港水文规范(JTS145-2-2013)[M].北
京:人民交通出版社,2013.
[2]HOGBENN,DACUNHAN,mentofanewglobal
capabilityforwaveclimatesynthesis[Z].SanFrancisco,ProcOCEANS’83,
1983.
[3]方钟圣.《西北太平洋波浪统计集》与其他图册资料的比较[J].中国海洋平
台,1997,12(1):29-32.
risonofwavedatafrom“wavestatisticsfor
NorthwestPacificOceanAreas”andotherSources[J].Chinaoffshore
Platform,1997,12(1)：29-32.
[4]方钟圣,金承仪,缪泉明.西北太平洋波浪统计集[M].北京:国防工业出版社,1996.
FANGZhongsheng,JINChengyi，atisticsfor
NorthwestPacificOceanAreas[M].Beijing:NationalDefenceIndustrial
Press,1996
[5]方钟圣,金承仪,缪泉明.中国海与西北太平洋波浪长期统计的导算方法[J].中国造
船,1994(4):21-35.
FANGZhongsheng,JINChengyi，vationofwave
statisticsforChinaseasandnorthwestPacificocean[J].Shipbuildingof
China,1994,4:21-35.
[6]LAGARIASJC,REEDSJA,WRIGHTMH,genceproperties
oftheNelder-Meadsimplexmethodinlowdimensions[J].SIAMJournalof
Optimization,1998,9(1):112-147.
[7]atisticsforthedesignofshipsandoceanstructures
[J].TransSNAME,1978,86:47-76.
[8]方钟圣,戴顺孙,金承仪.海洋特征波高和周期的长期联合分布及其应用[J].海洋学
报,1989,11(5):535-543.
FANGZhongsheng,DAIShunsun,-termjointdistribution
andapplicationoftheoceanwaveheightcharacteristicsandwave
period[J].ActaOceanologicaSincia,1989,11(5):535-543.