[发明专利]复杂构造区地震勘探数据采集观测系统评价方法及装置

权利要求书

1.一种复杂构造区地震勘探数据采集观测系统评价方法，其特征在于，包括：

根据历史地震数据和历史解释数据建立地质模型，利用网格间距离散化所述地质模型，得到网格地质模型；

建立多套地震勘探数据采集观测系统；

建立全波形反演模型的初始模型；

利用多套地震勘探数据采集观测系统中的每一套地震勘探数据采集观测系统执行如下步骤：

利用其中一套地震勘探数据采集观测系统对所述网格地质模型在时间空间域做有限差分正演模拟，得到模拟地震观测数据，所述模拟地震观测数据中包括正演参数；

执行如下迭代步骤：

利用同一套地震勘探数据采集观测系统和所述正演参数对所述初始模型做正演模拟，确定每个时刻炮点正向传播的地震波场和初始模型模拟数据；

计算模拟地震观测数据和初始模型模拟数据之间的观测数据残差；

利用所述初始模型，以观测数据残差的逆时数据作为检波点处的源进行正演模拟，确定不同时刻检波点残差的逆时传播波场；

利用所述地震波场和所述逆时传播波场计算速度更新量；

利用所述速度更新量更新所述初始模型；

使用更新后的初始模型继续执行迭代步骤，直至速度更新量满足预设条件时迭代结束，得到反演速度模型；

计算所述反演速度模型与所述网格地质模型每个网格的速度逼近程度和所述速度逼近程度的方差；

将利用多套地震勘探数据采集观测系统得到多个方差进行比较，获得最小方差，所述最小方差对应的地震勘探数据采集观测系统为满足勘探要求的地震勘探数据采集观测系统。

2.如权利要求1所述的复杂构造区地震勘探数据采集观测系统评价方法，其特征在于，采用如下方式建立全波形反演模型的初始模型：

利用预设的平滑半径平滑网格地质模型，获得所述全波形反演模型的初始模型。

3.如权利要求1所述的复杂构造区地震勘探数据采集观测系统评价方法，其特征在于，所述正演参数包括震源子波频率和记录时间。

4.如权利要求1所述的复杂构造区地震勘探数据采集观测系统评价方法，其特征在于，当地震勘探数据采集观测系统为二维地震勘探数据采集观测系统时，利用所述地震波场和所述逆时传播波场按照如下公式计算速度更新量：

dV(x,z)＝∑_tUf(s,t,x,z)*Ub(s,t,x,z)；

其中，dV(x,z)表示速度更新量；Uf(s,t,x,z)表示每个时刻炮点正向传播的地震波场，Ub(s,t,x,z)表示不同时刻检波点残差的逆时传播波场；s表示炮点编号；x，z表示二维空间位置；t表示记录时刻，t的取值范围为0到T，T表示记录时间长度；

当地震勘探数据采集观测系统为三维地震勘探数据采集观测系统时，利用所述地震波场和所述逆时传播波场按照如下公式计算速度更新量：

dV(x,y,z)＝∑_tUf(s,t,x,y,z)*Ub(s,t,x,y,z)；

其中，dV(x,y,z)表示速度更新量；Uf(s,t,x,y,z)表示每个时刻炮点正向传播的地震波场，Ub(s,t,x,y,z)表示不同时刻检波点残差的逆时传播波场；s表示炮点编号；x，y，z表示三维空间位置；t表示记录时刻，t的取值范围为0到T，T表示记录时间长度。

5.如权利要求4所述的复杂构造区地震勘探数据采集观测系统评价方法，当地震勘探数据采集观测系统为二维地震勘探数据采集观测系统时，利用所述速度更新量按照如下公式更新所述初始模型：

V₀′(x,z)＝V₀(x,z)+dV(x,z)；

其中，V₀′(x,z)表示更新后的初始模型；V₀(x,z)表示初始模型；

当地震勘探数据采集观测系统为三维地震勘探数据采集观测系统时，利用所述速度更新量按照如下公式更新所述初始模型：

V₀′(x,y,z)＝V₀(x,y,z)+dV(x,y,z)；

其中，V₀′(x,y,z)表示更新后的初始模型；V₀(x,y,z)表示初始模型。