水文频率分析以水文数据为基础。在实际水文计算中,由于观测条件限制,一些测站数据较短,难以满足水文频率分析计算条件,不能为工程规划设计与管理提供合理的数据支撑。建立多维联合分布模型,合理利用邻近长序列测站的特征信息,提高水文频率分析精确度,是不完全水文序列频率计算方法的一种新途径。多变量分布模型在水文领域中的应用已较为成熟。大量应用研究表明,传统多变量分布模型以变量间线性关系为基础,而Copula函数可刻画变量间线性或者非线性关系,优势明显。本文在吸收国内外不完全水文频率分析方法的基础上,应用Copula函数理论,探究基于Copula函数的不完全降水序列频率计算方法。以陕西省关中地区年降水序列为研究对象,选取短序列测站为设计站,邻近长序列测站为参证站。首先,采用P-Ⅲ分布、Gamma分布和Gumbel分布对设计站及其参证站等时段长度序列进行单变量降水频率分析,采用矩法、最大熵法、极大似然法和概率权重矩法4种参数估计方法计算3种线型分布参数值。用2种拟合优度评价指标对每一种线型下各参数估计值进行拟合优度检验,获得设计站及其参证站等时段长度序列联合分布模型的边缘分布。通过拟合检验,评判所选单变量概率分布模型的拟合效果。然后,计算设计站及其参证站等时段长度序列秩相关系数,在此基础上,采用相关性指标法和极大似然法计算Copula函数参数q的值,并选用均方根误差法(RMSE)、赤池信息准则(AIC)和贝叶斯信息准则(BIC)作为拟合优度评价标准,选出各站拟合度较好的Copula函数。通过分析经验频率与理论频率联合拟合图效果,并结合A-D检验,判定所选Copula函数的优劣。最后,计算出基于Copula函数的不完全年降水序列参数估计值,并计算出设计站在该方法下年降水量设计值,并与单变量降水频率分析下年降水量设计值进行对比,探究基于Copula函数的不完全降水频率分析,以期为工程水文提供理论依据。研究取得以下主要结论:(1)采用P-Ⅲ分布、Gamma分布和Gumbel分布对设计站及其参证站等时段长度序列进行单变量年降水频率分析,拟合研究区年降水量值。拟合检验结果表明:P-Ⅲ分布、Gamma分布和Gumbel分布对关中地区年降水序列拟合效果良好,可以用来进行关中地区年降水量频率分析;P-Ⅲ分布采用最大熵法估计参数值,更能符合年降水量取值要求。(2)用Gumbel-Hougaard(G-H)copula函数、Frank copula函数、Clayton copula函数,构建设计站及其参证站等时段长度序列多变量概率分布模型。二维联合分布拟合检验结果表明:蓝田站与长安站、临潼站与长安站、临潼站与咸阳站、旬邑站与彬县站、岐山站与宝鸡站、凤翔站与宝鸡站、凤县站与太白站选用Frank copula函数,临潼站与渭南站、华阴站与潼关站选用Clayton copula函数,其余各组联合分布选用G-H copula函数,均能较好地反映实测点据的联合分布,可用来对设计站及其参证站等时段长度序列年降水量进行多变量频率分析。(3)在确立单变量分布模型和Copula联合分布模型基础上,计算出基于Copula函数的不完全降水序列参数估计值,并采用拟合标准差法(SEF)和最大相对误差绝对值法(MARD)两种拟合度评价方法进行基于Copula函数的不完全降水频率计算结果评价。Copula函数模拟结果和拟合度评价结果表明:基于Copula函数的不完全降水序列频率计算方法合理可行,可为短序列测站降水频率分析提供新的计算途径。(4)根据基于Copula函数的不完全降水序列参数估计值,绘制设计站在新分布参数下的年降水量频率曲线,并推算出不同频率下的年降水量设计值。结果表明,新频率曲线对实测数据的拟合效果较好,尤其在频率曲线的中低部拟合效果好。(5)以Frank Copula为例,采用蒙特卡洛试验法探究基于Copula函数的不完全降水序列频率计算方法统计性能。偏差和标准误差值计算结果表明:基于Copula函数的不完全降水序列频率计算方法不偏性与单变量降水频率分析基本相同,有效性优于单变量降水频率分析方法。