一种基于波动方程的正演模拟方法及装置制造方法及图纸

本发明专利技术提供一种基于波动方程的正演模拟方法及装置，所述方法包括：根据声波测井数据，获得纵波速度场；根据背景速度模型对应的格林函数和地震子波频谱，获得频率域的背景波场；其中，所述背景速度模型是预设的，与所述纵波速度场对应；根据所述纵波速度场、所述背景速度模型、所述背景波场和全波场的线性方程组，获得频率域的声波全波场；其中，所述全波场的线性方程组是预设的。所述装置用于执行上述方法。本发明专利技术实施例提供的基于波动方程的正演模拟方法及装置，提高了正演模拟的效率和准确性。

【技术实现步骤摘要】
一种基于波动方程的正演模拟方法及装置
本专利技术涉及地质勘探
，具体涉及一种基于波动方程的正演模拟方法及装置。
技术介绍
现今，全波形反演在地下介质参数估计中起到了重要的作用，在反演过程中通常需要进行多次迭代，在迭代过程中需要通过正演模拟计算波场，正演模拟计算波场的效率和精度是是制约全波形反演应用的重要因素。现有技术中，正演模拟计算波场通常使用波动方程进行数值计算，常用的数值求解方法包括有限差分法、有限元法、伪谱法以及边界积分方程法等。例如，将短波逼近(WentzelKramers—Brillouin—Jeffery，简称WKBJ)方法近似引入背景介质格林函数的求解中，利用诺伊曼(Neumann)级数迭代方法求解一维声波方程的近似解，但该方法忽略了地下介质和波场之间严格的非线性特征；在求解声波方程散射问题时采用了一种介于Neumann级数迭代法和共轭梯度迭代方法之间的数值求解方法，但计算过程比较复杂。因此，如何提出一种基于波动方程的正演模拟方法，能够提高正演模拟的效率和准确性成为本领域需要解决的重要课题。专利技术内本文档来自技高网...

【技术保护点】
1.一种基于波动方程的正演模拟方法，其特征在于，包括：/n根据声波测井数据，获得纵波速度场；/n根据背景速度模型对应的格林函数和地震子波频谱，获得频率域的背景波场；其中，所述背景速度模型是预设的，与所述纵波速度场对应；/n根据所述纵波速度场、所述背景速度模型、所述背景波场和全波场的线性方程组，获得频率域的声波全波场；其中，所述全波场的线性方程组是预设的。/n

【技术特征摘要】
1.一种基于波动方程的正演模拟方法，其特征在于，包括：
根据声波测井数据，获得纵波速度场；
根据背景速度模型对应的格林函数和地震子波频谱，获得频率域的背景波场；其中，所述背景速度模型是预设的，与所述纵波速度场对应；
根据所述纵波速度场、所述背景速度模型、所述背景波场和全波场的线性方程组，获得频率域的声波全波场；其中，所述全波场的线性方程组是预设的。


2.根据权利要求1所述的方法，其特征在于，所述根据所述纵波速度场、所述背景速度模型、所述背景波场和全波场的线性方程组，获得频率域的声波全波场包括：
将所述纵波速度场、所述背景速度模型和所述背景波场输入到所述全波场的线性方程组，获得待求解线性方程组；
利用克雷洛夫子空间数值迭代法对所述待求解线性方程组进行数值求解，获得频率域的声波全波场。


3.根据权利要求2所述的方法，其特征在于，所述利用克雷洛夫子空间数值迭代法对所述待求解线性方程组进行数值求解，获得频率域的声波全波场包括：
引入对角矩阵预条件算子对所述待求解线性方程组进行数值求解。


4.根据权利要求1至3任一项所述的方法，其特征在于，所述背景速度模型为匀速背景速度模型或者速度均匀变化的背景速度模型。


5.一种基于波动方程的正演模拟装置，其特征在于，包括：
第一获得单元，用于根据声波测井数据，获得纵波速度场；

【专利技术属性】
技术研发人员：孙夕平，李景叶，李凌高，郭康康，于永才，陈小宏，刘国昌，马继涛，
申请(专利权)人：中国石油天然气股份有限公司，
类型：发明
国别省市：北京;11