适用于低频率波场传播随机边界条件.pdfVIP

513 SPG/SEG 深圳2011 国际地球物理会议论文集 适用于低频率波场的传播的随机边界条件 沈旭凯*, Robert G.Clapp, 斯坦福大学 随机边界条件（Clapp，2009 ）最早用于逆时偏移，因为它产生时间可逆的，不连贯的散射波场。其优点是源波场 可以先传播从零时刻到最大时刻，再同时与接收波场从最大时刻传播到零时刻。这大大降低存储波场的内存要求。这 种减少内存要求特别适合逆时偏移在GPU （Micikevicius 2008）等特殊硬件上的实现。随机边界散射的不连贯波场在 自相关成像条件后成为随机噪声，对最终的成像结果基本没有影响。 从理论上讲，这种随机的边界条件可以很容易地扩展到在所有不同频段的波场。然而，由于内存和计算需求，我们需 要避免使用过大的边界区域。因此，在低频波场使用随机边界时，使用相同的随机边界会导致边界区域的扩大。但如 果我们能保持随机边界与在逆时偏移中使用的随机边界相同的大小，计算效率会大大提高，特别是在三维情况下，我 们可以避免避免传播的波在边界域的额外的计算。接下来我们通过例子展示通过扰动随机颗粒的尺寸和形状，既可以 保持边界区域相对小，同时仍然有效的散射低频波场。 我们使用含有不同的“颗粒”大小和形状的随机边界的恒定速度场进行波场模拟。此外，当波场达到每个边界区域尽头 时我们使用反射边界条件。在所有的速度场的内部区域速度是相同的常量的值。边界区域中速度或者是恒定的，或者 是含有不同的“颗粒”大小和形状的随机边界。随机边界使用三种类型的“颗粒” （图一）：两个不同边长的立方体颗粒， 和有效的边长一样大的随机形颗粒。我们采用三维波动方程进行数值模拟，但我们只展示速度模型和波场的二维切面. 图一：使用不同的边界区域的速度场：a ）常速度; b)小边长立方体颗粒； c ）大边长立 图二显示了用图一中的速度场和六赫兹峰值频率震源进行数值模拟后3.5 秒时的波场。图二第一行显示了不同速度场 方体颗粒; d）具有相同的有效大边长值的随机形颗粒。 对应的波场的二维切片。对于如此低的频率源，小边长立方体颗粒的边界效果仅略高于常速度边界。只有当颗粒尺寸 变得与峰值波长相似时波场开始出现随机散射。不过使用大边长立方体随机颗粒时，波场中也有一些连贯的低频波长 信号的残余。图二的第二行显示了使用相同的震源和使用十六个不同的随机边界实现所得到的波场的叠加结果。现在 很明显，使用大边长立方体随机颗粒的边界在数值模拟后仍然有较强低频信号残余。由于这些大边长立方体随机颗粒 立方米的形状规则，散布特定的低频波场不能被散射。但大边长随机形颗粒能有效的散射低频信号。 1158 513 SPG/SEG 深圳2011 国际地球物理会议论文集 图二：使用六赫兹震源及不同的随机边界条件的速度场的模拟时域波场快照。a ）常速度; b)小边长 立方体颗粒； c ）大边长立方体颗粒; d）具有相同的有效大边长值的随机形颗粒。第二行显示使用 16 次速度实现的波场快照叠加，对应于边界条件: e）常速度; f)小边长立方体颗粒； g ）大边长立方 体颗粒; h）具有相同的有效大边长值的随机形颗粒。 从这个例子中，可以看出，使用大尺寸随机颗粒的随机的边界条件能有效的散射低频信号。随机形颗粒比规则形状的颗 粒散射效果更好。这样的随机边界条件能有逆时偏移以外的应用，例如波形反演。 1159