基于WKBZ理论耦合简正波_抛物方程理论.pdfVIP

( ) 中　国　科 　学 　　　　 A 　辑 　　　　　　 　 第 31 卷 　第 2 期 SCIENCE IN CHINA ( Series A) 　 2001 年 2 月 基于 W KBZ 理论的耦合简正波抛物方程理论 彭朝晖 　李风华 ( 中国科学院声学研究所声场声信息国家重点实验室 ,北京 100080) 摘要 　　研究水平变化海洋环境中声传播的快速数值预报方法. 在广义相积分 ( ) WKBZ 理论的基础上 ,提出了一种能够快速 、准确地求解简正波本征值的方法 ,并将 ( ) 此方法应用于耦合简正波抛物方程 CMPE 理论. 分析了两种典型声速剖面的简 ( ) 正波本征值 ,并计算了一个楔形海区中的声传播问题 JASA 标准问题 . 计算结果表 明本方法是一种快速而精确的数值方法. 关键词 　　WKBZ 　本征值 　耦合简正波抛物方程 　声传播 对于水平不变和水平变化相对缓慢的海区 ,近年来发展的声场计算模型很多 ,较为著名的 [ 1] ( ) [2 ] [3 ] 有绝热简正波理论 、抛物方程 PE 方法 、射线理论 等等. 在这方面 ,我国学者提出的广 ( ) [4 ] ( ) [5 ] 义相积分 WKBZ 理论 及波束位移射线简正波 BDRM 理论 分别在深海和浅海环境中 ,实 现了声场的快速 、精确预报. 但是 ,对于水平变化较大的海洋环境 ,必须考虑水平变化对声场 的影响. 这时 ,用绝热近似方法来计算声场会产生较大的偏差 ; 目前国际上流行的耦合简正波 [6 ] ( ) 理论 ,尽管计算精度很高 ,但是计算时间很长 ,不能广泛应用 ;抛物方程 PE 方法计算精度 比耦合简正波的稍差 ,只能计算远场 ,且处理向后散射有困难. 文献[7 ]将耦合简正波理论与抛物方程方法结合起来 ,提出了一种全新的声场计算模型 : ( ) 耦合简正波抛物方程 CMPE 方法. CMPE 方法在垂直方向采用本地简正波分析 ,这就克服了 抛物方程方法只能计算远场 ,且在频率较高时垂直网格的划分必须加密 ,使得计算时间成几何 级数增加 ,因而很难用于高频问题的缺点[2 ] . CMPE 方法在水平方向采用抛物方程法求解简 正波系数方程 ,可以克服由于耦合简正波理论的分段水平均匀近似带来的缺陷 ,耦合系数中考 虑了海底倾斜的影响 ,可以加大水平步长 , 同时可以很方便地应用于三维声传播问题. 在用 CMPE 方法计算声传播问题的过程中 ,本地简正波和耦合系数的计算将花去很大部 分的时间. 因此 ,快速而精确地计算本地简正波本征值和本征函数成为提高计算效率的关 键[8 ] . 另外 ,在水平变化声场的分析中 ,尤其在近场的分析中 , 需要计算高号简正波本征值. 本文在 WKBZ 理论的基础上 ,提出一种可以快速而精确地求解本地简正波高号本征值的算法. 然后 ,将改进后的WKBZ 理论应用于