复杂地区3D近地表速度建模方法.pdfVIP

大会交流论文 面，所以，计算时我们将三维速度体设计为沿Y方向速度不变)，图2是随机抽出一炮的比较结果．可 见，两种方法的计算结果基本上没有差别；尽管只是用二维剖面计算结果进行比较，但也在一定程度上 说明了“三维变步长扰动线追踪法”计算的初至波走时是正确的。 图l三维介质射线追踪示意图 o+分别为变步长扰动法炮点左、右侧T0值 。+分别为波动方程有限差分法炮点左、右侧T0值 图2变步长扰动法射线追踪结果与波动方程有限差分法正演结果比较 (2)速度反演 根据初至波走时来建立复杂地区3I)近地表速度模型，我们使用了最小二乘法。根据最小二乘法原 理，对于已知的每一条射线，我们可以近似地用一线性方程表示： 24 ——一一一一叁!兰垫垒苎 、“(1) 7 AI：。艺a0矿T,AV， 其中：M为网格节点个翥。7疗巧 综合辨舂E籍黪熬线方程(缓翔共有N条)，l§么，我g】就可敬缮爨～线性方糕组： allAv,+Q12△K+-··+aimA％=△五 a2l△K+疗22AV：+…+a2占f△％=A是 码l△K+吼2△吒+t··+群3埘△‰=矗正(2) l口州△巧+area巧牛…+a!vMA％=△巧 蓄轰南蠛影巢纛儿扣1’2，…∥。 疆论计算走时之差，AV为阏格节点速度修改缀。(3)武是我们求解反演阍题的方糕组。镲蟹矩阵A一 般在正演计算时的射线追踪过程中求得。以往槛建立偏肆矩阵时，仅考虑射线所经过的节点遵度矿，对走 露瓣影响，这里，我霞j不但考虑y，对走时的影响，疑对还考虑致邻近节焱速度对惫霹△置豹影蛹。 通常方程(3)是一大型的、稀疏的超定方程组，截这里，我们采用求解线性最小二乘闷题的豪斯 犄尔德(Housellolder)变换法来解方程蛆(3)。 由于钢至渡_鸯于绫轨迹与横鍪|参数露关，莲戴，秘至波走靖爱滚是一嵩囊菲线性溺题。在计算过程孛 采用了多次迭代的办法，反复修改模型，直至爝终得到满意的结果为止。 3。疆澎壤叠l试验及效暴分橱 为了检验本方法的正确性及其计算程序的可靠性，我们设计了一地面起伏、纵、横向速度都变化的三 绻速凄模爨搭。逢舞起佚形态如整3，(其中(a)建蘑等态线，(b)是三维立终瑟)，模型体瓣搔划分始 前述图1(a)所示。 试验珏寸“实测”初至波惫时资料娥通过对设计好的三维理论模型，用“三维变步长扰动射线追踪法” 正演}壤得，然后根据走骶信息送行反演骨簿。为了舞壹鬣戆流鹾反演效莱，鼯疆论接登、裙戆稹裂 和反演结果都在x、Y及水平方向上用切片图的形式来加以说明。 y 黔 扭豫l澎等裹缱自霪鳢鸯韧纛经过处， 岛)鼗髂避彭爨 图3理论模型地形图 复杂地区3D近地表速度建模方法 作者： 李满树 作者单位： 中石化石油勘探开发研究院物探研究所(江苏南京) 引用本文格式：李满树 复杂地区3D近地表速度建模方法[会议论文] 2003