[发明专利]一种实现地震波场数值模拟方法有效
申请号: | 201310726137.0 | 申请日: | 2013-12-25 |
公开(公告)号: | CN103699798B | 公开(公告)日: | 2017-01-25 |
发明(设计)人: | 秦臻;宋建勇;马晓宇;宋雪娟 | 申请(专利权)人: | 中国石油天然气股份有限公司 |
主分类号: | G06F17/50 | 分类号: | G06F17/50 |
代理公司: | 北京三友知识产权代理有限公司11127 | 代理人: | 贾磊 |
地址: | 100007 *** | 国省代码: | 北京;11 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | 本发明涉及地球物理勘探技术领域,特别是关于地震波场数值模拟技术领域,具体的讲是一种地震波场数值模拟的方法,包括在波数域使用低阶有理分式有限差分算子逼近超高阶显式有限差分算子;计算所述低阶有理分式有限差分算子的系数;将有理分式有限差分算子进行分解,利用傅里叶变换将低阶有理分式有限差分算子从波数域变换到空间域;在空间域使用有理分式有限差分法分别计算不同方向上的空间二阶偏导数项;在时间方向上进行迭代实现地震波场的数值模拟。通过本发明实施例的方法,使用低阶有理分式有限差分法来实现地震波数值模拟,提高了地震波场数值模拟的精度和计算效率。 | ||
搜索关键词: | 一种 实现 地震波 数值 模拟 方法 | ||
【主权项】:
一种实现地震波场数值模拟方法,其特征在于包括,步骤101,在波数域使用低阶有理分式有限差分算子逼近超高阶显式有限差分算子;其中,使用如下公式逼近超高阶显式有限差分算子:2Σm=0Mbmcos(mkhΔh)≈Σl=1Lcl[2-2cos(lkhΔh)]Πl=1L[dl+2cos(khΔh)],khΔh∈[-π,π];]]>其中,方程左端项为波数域的具备2M阶差分精度的超高阶显式差分算子,右端项的分子项为波数域的2L阶差分精度显式差分算子的求和,在空间域通过显式差分法计算;方程右端项的分母项为波数域的二阶隐式差分算子的连乘,在空间域通过隐式差分法求解L个三对角线性方程组,cl和dl为待求的低阶有理分式有限差分的系数,kh为波数域的自变量,Δh为空间采样间隔,bm为超高阶显式差分算子的系数,其中L、M为正整数;步骤102,计算所述低阶有理分式有限差分算子的系数;步骤103,将有理分式有限差分算子进行分解,利用傅里叶变换将低阶有理分式有限差分算子从波数域变换到空间域;步骤104,在空间域使用有理分式有限差分法分别计算不同方向上的空间二阶偏导数项;步骤105,在时间方向上进行迭代实现地震波场的数值模拟。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中国石油天然气股份有限公司,未经中国石油天然气股份有限公司许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/patent/201310726137.0/,转载请声明来源钻瓜专利网。