曹 渊 王文科 王铁良 刘 峰 王英杰

1 (长安大学 西安 710054)

2 (西北核技术研究所 西安 710024)

3 (西安通信学院 西安 710106)

核素迁移的复杂机制和地介质的高度非均质性，使核素迁移计算参数较多并有较大不确定性。参数灵敏度分析可有效识别系统的敏感参数，为参数校准和不确定性分析奠定基础[1,2]。目前常用的参数灵敏度分析方法可分为局部灵敏度分析法和全局灵敏度分析法[3–8]。前者用于检验单个参数变化对数值模拟结果的影响，分析时仅改变某一待分析参数的值，其他参数保持不变。该方法易于操作，但忽略了参数间相互作用对模型输出结果的影响，故局限性较大。全局灵敏度分析允许多个参数同时变化，考虑参数间相互作用的影响，分析结果更接近实际，主要有多元回归法、Morris法、傅里叶幅度灵敏度检验法、基于方差分析的Sobol法、偏秩相关分析法[9–11]等。国外在数值模拟研究中，特别重视参数灵敏度分析研究，如美国地质调查局开发的JUPITER[12]和美国 Sandia国家实验室开发的Dakota(专门的模型分析平台库)中均包含参数灵敏度分析的模型和程序。

本文建立了基于拉丁超立方抽样和 Spearman偏秩相关分析方法的核素迁移参数灵敏度分析模型，定量计算了3H在地下水中迁移的参数灵敏度。

1 参数灵敏度分析方法

1.1 参数不确定性的量化

布类型、取值范围和均值。

1.2 参数空间的拉丁超立方抽样(LHS)

分析不确定性参数的灵敏度，需进行参数抽样。对于多维参数空间抽样，LHS方法获得的数据点的代表性远好于常规抽样方法，从而大大节省样本量。……