原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
维普资讯
维普资讯
622 石 油 物 探 第 46卷
Gx( 一蔷 ，r，s)w ，r， 当此矢量与反射面单位法向矢量的标量积大于某
一 给定阈值时，才在此粗网格附近的细网格上进行
P (X，r，s)AN(X，r，s) · 计算。图 2为Kirchhoff积分的阈值及实际计算
(t— Tx(X，，，s))dS (1) 区域，Kirchhoff积分只在稳相点附近的单元 (黄色
区域)上进行。
式中：GN为 3×3二阶张量，其分量对应于地震波
不同的激发及接收方向；s，r和X分别为震源、检
波器和地下界面反射点的位置 ；N表示地震波之
间不同的转换模式 (本文只考虑 PP，PS，SP和 SS
地震波，如果需要考虑多次反射或波型转换，则需
对地震波的传播路径进行编码，然后利用 (1)式分
别进行计算 ])；Ⅱ表示地下某一反射面；R为反射
系数；T和A分别为地震波的旅行时和振幅；W 和
P分别为震源波场的辐射花样和极性 。(1)式的离
散式为
地下反射界面
GV．(r f)一∑R…XWxp～XAV．·
i，J 图2 Kirchhoff积分的阈值及实际计算区域
alI (f— )皿．， (2)
皿 ．… … … 。 一
3 模型试算
式中：Ⅱ 为离散化后反射面上的粗网格单元。(2)
式中的积分可进一步离散化为细网格上的离散相
为验证方法的有效性，我们选择了 —lOooo1TI，
加，这隐含 了在粗网格上如下假设成立 ：
y=lO000m，z=5000m的弹性介质模型进行了数
AR 一 0 △ N 一 0 △P 一 0 △AN一 0
值模拟，模型底部为一弯曲的反射面。模型中的P
(3) 波速度、S波速度和密度均只沿 和z方向变化，
即R ，w ，P 和A 的变化可忽略不计。地震波 方向保持不变。沿 方向共布置了21条测线，在
旅行时可表示为 每条测线上设置了81个炮点，每一炮由81个检波
器接收。根据射线追踪要求，速度均用高斯函数进
△ 一 T ·Ax+去△ ·( T )·Ax
行了平滑。波场计算在一台主频为443MHz的
(4) SUNUltra-10工作站上进行，耗时约 20h。图3
(4)式中二阶项 TN可由旁轴射线追踪求得。 给出了模型和数值模拟的速度参数。
如果将(1)式或(2)式变换到频率域，则可看出上述 利用三分量震源和三分量检波器进行了数值
假设成立，因为在高频近似下旅行时 T的形式为 模拟，检波器位于 一5000m和 z一0m处。图4a
eL或 e-i。显然，在一个粗网格的范围内，R和 和图4b分别为 Z分量震源激发、Z分量检波器和
w 及 P 的变化应 比关于 T 的指数项更为平缓。 X分量检波器接收的地震记录，震源位于 一4950
(2)式及(4)式为数值模拟实现过程的基础 。 I
文档评论（0）