用Fermat原理作折射波正演算法.docVIP

 物探化探计算技术 GEO PH YS ICA L A N D 第 21 卷 第 4 期 1999 年 11 月 COM PU T IN G T ECHN IQ U E S FO R V o l121 N o. 4 N o v. 1999 GEO CH EM ICA L EX PL O RA T ION 用 F e rm a t 原理作折射波正演算法 戴 云 张建中 钟本善 周熙襄 (成都理工学院) 【摘 要】 本文提出了基于 F e rm a t 原理的折射波正演算法, 它以模型的离散划分和射线路 径搜索为基础, 适用于任意地形和起伏折射界面的速度可横向变化的模型。 本算法快速且有 效, 对作沙漠和山地折射静校的研究有重要的意义。 【关键词】 折射波正演 费马原理 折射静校正 【中图法分类号】 P 631. 4+ 26 0 前言 静校正是沙漠和山地等低速带结构复杂地区地震勘探的难题和关键, 长期以来未能得到 圆满的解决1 。利用大炮折射初至的折射静校正是目前静校正技术中广泛使用且有效的方法。 对于有复杂地表低速带结构的山地, 如何作折射静校正是一世界性的难题, 因此需要对复杂的 山地建模作折射初至正演, 来研究山地静校。 当折射界面起伏不平时, 常用的折射波正演射线 追踪法2 不易实现, 而用 F e rm a t 原理作折射波正演却易实现, 它适用于任意地形和起伏折射 界面的速度可横向变化的模型, 且计算有效。 1 方法原理 由费马原理, 地震波在介质中沿射线传播的时间比沿其它任何路径传播的时间要少; 或者 说波沿所花的时间最少的路径传播。 由此可以确定地震波在已知速度介质中的射线形状和分 布情况。 对于均匀介质, 波的传播路线是直线, 对于非均匀介质, 波的传播路线是曲线3 。 对于水平界面或倾斜界面的层状介质, 可以由公式法根据折射定律导出折射旅行时的正 演计算公式。 对于起伏不平的折射界面, 地震折射初至的正演计算没有现成的计算公式, 因此 本文采用费马原理进行近似计算。 假设一地表起伏、折射界面起伏的如图 1 所示的二层模型。 地面的高程函数为 Z 1 (x ) , 折 射界面的函数为 Z 2 (x )。若模型的水平范围为 (X m in , X m ax ) , 在 x 方向上用m 成 m 21 段, m 条等分线的水平坐标为 X i = X m in + ?X ( i- 1) 条等分线将模型分 收稿日期: 1999- 08- 30 地址: 成都市二仙桥东三路 1 号 ( 邮政编码: 610059) 图 1 用费马原理进行折射初至正演示意图 其中, X 1 = Xm in , X m = X m ax , ?X = (X m ax 2X m in ) ?(m - 1)。 第 i 条等分线与地面的交点坐标为 ( x i , Z 1 (x i ) ) , 与折射界面的交点坐标为 (x i , Z 2 (x i ) ) , 这里近似认为折射界面在相邻等分线间 为直线。模型的第一层是速度为 V 1 的均匀介质; 第二层速度是变化的, 记为 V 2 ( i) , ( i= 1, m - 1) ; V 2 ( i) V 1。 下面讨论折射正演问题。 设炮点 S 和检波点 R 都位于地面上, 且分别位于等分线与地面 相交的第 i 和第 l 处, 起始折射点 A 和出射折射点 B 分别位于等分线与折射面相交的第 j 和 第 k 处, i j k l。当折射面为弯曲界面时, A 和B 位于交点 j 和 k 的假设是近似的, 且最大 走时误差m ax (E t (x ) ) ?x ?V 2 (x )。这样折射波从 S 到A 沿折射面滑行到B 再出射到 R 所用 的旅行时 T : T = tS A + tA B + tB R k - 1 rS A + rB R + 6 rn (1) = n= j V 2 (n ) V 1 其中, rS A = (Z 1 (x i ) - Z 2 (x j ) ) × (Z 1 (x i ) - Z 2 (x j ) ) + (x i - x j ) × (x i - x j ) rn = (Z 2 (x n ) - Z 2 (x n+ 1 ) ) × (Z 2 (x n ) - Z 2 (x n+ 1 ) ) + ? x ×? x (Z 1 (x l ) - Z 2 (x k ) ) × (Z 1 (x l ) - Z 2 (x k ) ) + (x l - x k ) × (x l - x k ) rB R = 利用 F e