地震偏移方程的差分求解新方法.pdf

科赫 科l技l论I坛 地震偏移方程的差分求解新方法 刘铁民1朱学东1王威z (1黑龙江八一农垦大学，黑龙江大庆1633ll2、哈尔滨华夏计算机职业技术学院，黑龙江哈尔滨150025) 摘要：构造了二维45。上行波方程的一种差分格式，该差分格式是隐格式但却是半隐式求解，利用此差分格式求偏移问题数值解时，既保证 了算法的稳定性，又使得运算量大大减少。 关键词：隐格式；稳定性；单程渡方程 Ab翻玎Ict：Thi8 di矗．e他nceschen他of45 it is喜dvedwitlI method． paperconst珈ctB卸implicit deg”弛one—waveequation，锄d sellli—i瑚dUcit Jinle TK亭djfrerenceschemeb鹊蛐exceU朗ts￡abiljlyand calculati加whendojJlgmi鲫i∞． KeywOrds：implicitdi矗bIence∞heme；st幽ility；m伊wavee口Ilati仰 引言 零。因而系数矩阵是正定矩阵，因此求此方程组 ，羲筹譬嘲+刃“)+芸△：以“+￡衅∽一dI川)(5) 在地震勘探的偏移问题中，需对上行波的 的解法很多，可采用直接法中的追赶法或迭代 高阶逼近方程进行数值求解．在实际应用中， 法中的雅可比迭代法和高斯一赛德尔迭代法．当 对于地层倾角不是很大的地区，利用450上行 嚣，妻铷鼍一舡粕㈤ 未知量的个数特别大时，为了保证计算的精度 波方程逼近式偏移处理就能满足要求，文献[1】 ∥+删=碟+塞衅魄一粕) (7)和计算的速度，可采用共轭梯度法。逐层逆推求 给出了一种二维45。上行波方程的逼近公式， 将式(7)代入式(2)得 方程组的解即可求得最终偏移问题的解。由于 为了对该逼近式进行数值求解，本文构造了逼 (1一A△!x面H+瞄1)=【l+(％+反一％鸠X，磊+础)+该差分格式是隐格式保证了算法的无条件稳定 近方程的一个隐式差分格式，再通过一个简单 成△：面墨 (8)性，而具体求解时它却是半隐式求解，大大减 的函数变换将其变为一种半隐式求解的差分格 式中-％：羼祟，店：魄+％一嘞 少了运算量。 式，这样就既能保证算法的稳定性，又使得算 ％2以至石’店2‘以十％一％J 3结论 法的实现变得比较容易．通过对资料的处理， 再由式(6)得 以上构造了450上行波方程的一种隐式差 说明了该方法的有效性．本方法很容易推广到 分格式。文中在算法设计中，应用隐式差分格式 锚=，耋筹△：(面*+面)+筹衅(麻”+，‰)+舍衅‰ 三维问题的处理。 以达到好的稳定性，而在求解时利用简单变换 l二维450上行波方程逼近式 进一步整理得 使差分格式变为半隐式求解，通过对实际资料 在浮动坐标下，由文献…中所导出的二维 的偏移处理验证了所构造的差分格式的有效 45上行波方程，结合定解条件可得偏移求解的 ㈥… 性。 混合问题是 等△：(po十2+此+1)+等△：p：+t蓉凳雹篡攀抵 参考文献