二维不均匀介质中厄米高斯波束.pdfVIP

中国地球物理201I ·489· (15)信息技术与地球物理 二维不均匀介质中的厄米高斯波束 in2D ofwavefields inhomogeneous Computation media卅ermiteGaussianbeams 石秀林+孙建国 ShiXiulin‘Sun Jianguo 吉林大学地球探测科学与技术学院长春130026 自高斯波束方法由Popov和Cerveny引入勘探地震领域以来，凭借其计算稳定且无奇异点的优良特性， 被越来越广泛地应用于该领域的多个方向【l】。高斯波束偏移已被视为一种非常有效的偏移成像技术，并正 致力于深度偏移和真振幅偏移的研究中[21；Cerveny已研究和推导了单高斯波束在不均匀各向异性层状介 质中的性质，并准备进一步讨论高斯波束在该介质中的积分求和问题；传统高斯波束法的推导是基于拉格 朗日机制，受限于中心射线坐标这种局部化坐标，一种基于全局笛卡尔坐标的欧拉形式高斯波束方法已经 提出【3l，既避免了常规高斯波束法计算中存在大量坐标变换的问题，又对复合源点问题有更好的适应性， 其更多的优点仍在进一步的探索中：此外，高斯波束法还被应用于地形波这种大气波的研究。 理论上，高斯波束只是高频波动方程的一个特解，而厄米高斯波束才是该方程的一组完备基础解系。 利用厄米高斯波束进行地震波场的计算，可在保留高斯波束法优良特性的前提下，提供更加丰富的振幅剖 面选择，虽然目前还没有足够的实验证据表明这种丰富振幅剖面选择的优越性，但至少可为高斯波束研究 提供一种新的理论探索。 普通高斯波束法的基本思想是将整个波场近似成是由无数个基本波的波场叠加而成，而每个基本波的波 场能量主要沿中心射线传播，其振幅剖面为钟形曲线，即极大值在中心射线上，垂直方向距离中心射线越远， 幅值成指数衰减。因为振幅幅值的这种分布，符合高斯分布，所以这种基本波被称为高斯波束。理论上，由 于高斯波束是高频波动方程的解，所以波场可近似为高斯波束场的叠加，而通过推导，厄米高斯波束同样是 高频波动方程的解，因此，经证明，波场同样可近似成无数个厄米高斯波束场的叠加。由于高阶厄米高斯含 有厄米多项式，而零阶厄米高斯波束就变为普通高斯波束，所以这样的波束被称为厄米高斯波束。 将普通高斯波束法拓展到厄米高斯波束法，需要解决两个问题，分别是权函数的确定和波束参数的选择。 权函数的物理意义是基本波的初始能量分布，它只跟射线坐标有关，即随中心射线不同而变换，当中心射线 固定时，它沿中心射线是一个常量。普通高斯波束法为了确定权函数，采用的是稳相法，近似求解广义积分， 然后将其与传统射线法等经典方法相比较，求得理论公式。但由于厄米高斯波束包含了厄米多项式，使得稳 相法已经不能得到普遍适用，需要另外寻找新的办法，这个问题最终可由复变函数的一些理论来解决。 波束参数选择的问题，是普通高斯波束法中的经典问题，其物理实质是波束参数对波束物理形态的控制。 普通高斯波束方法中，由于振幅剖面上的形态为钟形曲线，是固定的，所以只需考虑随s增加有效长度上变化 的影响，通常由初始有效长度厶来确定。本文通过研究发现，将高斯波束拓展到厄米高斯波束，由于厄米多 项式的引入，可得到的振幅剖面，其能量分布多样，具体能量分布的形态将更多由线性叠加的常数项来决定。 同时，厄米高斯波束也存在有效长度L的概念，高斯波束法中有关波束参数的讨论，仍将适用于厄米高斯波束法。 本文对厄米高斯波束进行均匀介质和梯度介质的数值模拟，发现如果波束参数选择适当，可得到非常 稳定的计算结果，其整体模拟误差，可控制在非常小的范围内，完全保留了普通高斯波束法的优良特性。 由于目前模型仍然比较简单，无法验证丰富振幅剖面选择的优劣性，但随着数值模拟模型的改进，相信厄 米高斯波束会得到更加深步的研究与应用。 参考文献 L withGaussian