[发明专利]一种山区地表情况下的三角网格逆时偏移方法

摘要

本发明涉及一种山区地表情况下的三角网格逆时偏移方法，其特征在于模型剖分的过程为将偏移涉及的地质区域利用三角网格进行剖分，三角网格之间无重叠区域和空白区域；剖分完成后，模型变为大量三角网格的集合体，每一个结点可能是多个三角形的共同顶点；具有共同顶点的多个三角网格构成一个计算单元，该计算单元为一个多边形；对计算单元中的每一个三角形取重心，将各三角形的重心连接起来形成一个新的多边形，该新的多边形为计算单元的对偶单元。

主权项

一种山区地表情况下的三角网格逆时偏移方法，其特征在于：模型剖分的过程为：将偏移涉及的地质区域利用三角网格进行剖分，三角网格之间无重叠区域和空白区域；剖分完成后，模型变为大量三角网格的集合体，每一个结点可能是多个三角形的共同顶点；具有共同顶点的多个三角网格构成一个计算单元，该计算单元为一个多边形；对计算单元中的每一个三角形取重心，将各三角形的重心连接起来形成一个新的多边形，该新的多边形为计算单元的对偶单元；在该三角网格化的速度模型中，波场延拓公式表示为：u0(j-1)=2u0(j)-u0(j+1)+(Δt)2Σi=1NAivi2Σi=1Nsihi[ui(j)-u0(j)]]]>式中，hi表示计算单元中各三角形的共同顶点与计算单元多边形的第i个顶点的连线长度；si表示对偶单元多边形的第i个边的边长；N表示计算单元中三角网格的数量，Ai是包含公共顶点p0的对偶单元中被三角形网格的边所分割成的第i个小的多边形的面积；式中，其余符号的含义如下：ui(j)表示第i个网格点处第j个时刻的波场；同理，u0(j‑1)，u0(j)，u0(j+1)分别表示第0个网格点处第j‑1、j、j+1时刻的波场；vi为与单元的公共顶点p0相邻的第i个网格点处的介质速度；Δt表示时间间隔。