海洋声学中缓变分层波导的声波传达计算.pdf

摘 要 在声学、电磁学、地震学和一些其他的应用中，有许多在具有分层缓变波 导中大尺度的波传播的问题。如果使用直接方法，譬如有限元法、有限差分法 等，往往会导致求解大规模的线性方程组，这使得使得算法在内存和计算量的 花费非常大，计算效率不高。在本文中，我们针对具有弯曲界面或者弯曲边界 的波导，提出了适合于步进计算一般的坐标变换变换的方法，在坐标变换后的 to 坐标系中，使用基于DtN(DirichletNeumann)映射的正向基本解算子步进 问题和反问题。 在海洋声学中，海洋环境可以看作是一种分层缓变的具有多个弯曲内部 界面和多层介质的声波导。在这样的波导中如果要使用步迸方法大步长的求 解Helmholtz方程，首先需要用适合于步进计算坐标变换将波导中弯曲的界面或 者弯曲边界拉直。本文中给出两种坐标变换的方法：第一种是局部的解析正交 坐标变换法，局部解析坐标变换用一个解析的积分方程联系新旧坐标，使用牛 顿法求解这～积分方程实现新旧坐标之间的转换。坐标变换必然相应的引起方 程的变换，我们使用的方程变换将Helmholtz方程转化成为一个可以用步进方法 求解的形式。对于具有一个平的顶部、平的底部、礼个弯曲的内部界面的二维波 导，我们给出了在变换方程系数在各个介质层一般形式的解析表达式，从而完 成方程变换。然而，局部解析坐标变换方法需要满足在两个相邻界面之间存在 一条水平直线的介质可分性条件。针对不满足可分性条件的情况，我们又提出 了一种使用经典求解常微分方程的Runge-Kutta方法计算坐标变换的数值坐标 变换方法，这一方法不仅具有更广的适用范围，而且不必引入辅助界面，从而更 好的逼近原始问题。 对于变换之后的计算坐标系中的正问题，我们使用正向基本解步进方法来 求解方程。FFOMM基于DtN算子的单向重建，将边界值问题转化为初始值问 题，对于DtN算子和基本解算子满足的方程使用大步长方法离散在传播方向上 的变量，用截断特征值展开方法逼近这些算子。FFOMM对于求解大尺度缓变 波导中波的传播问题非常有效。 在具有分层缓变特性的波导中的波的大尺度传播问题中，反问题有非常重 ii 海洋声学中缓变分层波导的声波传播计算 要的实际意义。一般说来，海洋声学中实际的反问题都很复杂。求解海洋中反 问题一般使用两类方法：一种是基于正向方法(如FFOMM)的迭代法，另一种 为反向方法，反向方法的运算量和计算速度都相比正向方法有很大的优势。基 于FFOMM，我们提出了逆向基本解算子步进方法。在IFOMM中使用波导中存 在的传播模的个数作为正则化参数，用截断的奇异值分解方法求解步进过程中 出现的病态线性方程。数值计算证明，这个方法用于求解Helmholtz方程的反边 界值有效，精确并且对于初始波源传播部分稳定，对于海洋声波中实际的反问 题研究有着一定的意义。 关键词：Helmholtz方程，局部正交变换，DtN重建，步进方法，界面，边界值 反问题，传播模，泄露模，物理坐标系，计算坐标系，正向基本解算子步进方法， 逆向基本解算子步进方法 Abstract andother In migration applications， acoustics，electro-magnetism，seismic thereaxe scalewave directnumericalcom- problems．A manylarge propagation be sincethese like couldveryexpensive putingfinite-element，finite-difference