[发明专利]一种利用三维地震恢复相对微古地貌的工业化流程

权利要求书

1.一种利用三维地震恢复相对微古地貌的工业化流程，其特征在于，包括以下步骤：

步骤(1)：在三维地震资料上，利用井震结合的地震资料解释方法，分别获得目的地层顶面深度的数据网格T_top和底面深度的数据网格T_bot；

步骤(2)：目的地层底面与顶面相减得到目的地层的厚度数据网格T_h；

其中，步骤(2)所述的计算公式为T_h＝T_bot-T_top；

步骤(3)：选取研究区任意一口井A，在地层厚度网格T_h上读取A井处地震解释得到的地层厚度h_t，利用A井的测井曲线数据获得A井处的地层厚度h_l，两者相减得到A井点处地震解释结果的误差D_A；

其中，步骤(3)所述的计算公式为D_A＝h_l-h_t；

步骤(4)：利用步骤(3)的方法得到研究区每一口井点处的地层厚度误差，并对这些数据点进行插值和外推计算，得到覆盖研究区的误差网格T_d；

步骤(5)：将步骤(2)的T_h与步骤(4)的T_d相加，得到井震校正后的地层厚度数据网格T_hc；

其中，步骤(5)所述的计算公式为T_hc＝T_h+T_d；

步骤(6)：利用三维地震属性分析计算目的层顶面的地震倾角数据网格T_dat，对该网格上每个点求取余弦值，获得目的层顶面的倾角余弦值网格T_cost；

步骤(7)：利用三维地震属性分析计算目的层底面的地震倾角数据网格T_dab，对该网格上每个点求取余弦值，获得目的层底面的倾角余弦值网格T_cosb；

步骤(8)：将步骤(6)和步骤(7)分别求取的目的层顶面和底面地层倾角余弦值网格进行算数平均，获得研究区目的层的平均地层倾角余弦值网格T_cosave；

其中，步骤(8)所述的计算公式为T_cosave＝(T_cost+T_cosb)/2；

步骤(9)：将步骤(5)获得的井震校正后的地层厚度数据网格与步骤(8)获得的目的层的平均地层倾角余弦值网格相乘，得到目的地层真厚度数据网格T_hr；

其中，步骤(9)所述的计算公式为T_hr＝T_hc×T_cosave；

步骤(10)：利用步骤(9)得到的目的地层真厚度数据网格T_hr制作等值线图，作为目的层沉积时期的相对古地貌形态。