一种利用旋回分析法恢复不整合剥蚀厚度的方法与流程

技术特征：

1.一种利用旋回分析法恢复不整合剥蚀厚度的方法，其特征在于，包括以下步骤：

沉积速率曲线绘制：选取研究区，结合研究区的地质资料，通过去压实校正计算压实系数和压实前地层的原始厚度，根据地层年代表确定每个地层沉积的起始时间和结束时间，计算各沉积期的平均地层沉积速率，绘制沉积速率曲线；

沉积速率曲线的频谱分析：根据沉积速率曲线的旋回规律确定缺失曲线段；根据缺失曲线段前后沉积速率曲线的周期值，推算缺失曲线段对应的周期、缺失曲线段波峰和波谷对应的横坐标；对沉积速率曲线进行分频处理得到对应低频子波和高频子波，根据缺失曲线段前后对应的低频子波和高频子波的波峰和波谷的纵坐标,推算缺失曲线段波峰和波谷对应的纵坐标；

计算不整合剥蚀厚度：根据频谱分析结果恢复完整的沉积速率曲线，计算不整合剥蚀厚度，所述不整合剥蚀厚度的计算公式如下：

其中，H_e为不整合剥蚀厚度，单位：m；F(t)为完整沉积速率曲线函数；H_oc为在沉积时间段内对F(t)积分得到的地层原始厚度；H_r为现今残余厚度，单位：m；C为总压实系数；t₁为剥蚀地层沉积开始时间，单位：Ma；t₂为剥蚀地层沉积结束时间，单位：Ma。

2.根据权利要求1所述的利用旋回分析法恢复不整合剥蚀厚度的方法，其特征在于，在沉积速率曲线绘制步骤中，所述压实前地层的原始厚度的计算公式为：

H₀＝H_s/C_s+H_m/C_m

其中，H₀为地层压实前的原始厚度，单位：m；H_s为砂岩厚度，单位：m；C_s为砂岩压实系数，利用去压实校正方法计算获得；H_m为泥岩厚度，单位：m；C_m为泥岩压实系数，利用去压实校正方法计算获得；

在计算不整合剥蚀厚度时，所述总压实系数的计算公式为：

其中，C为总压实系数；H_s为砂岩厚度，单位：m；C_s为砂岩压实系数，利用去压实校正方法计算获得；H_m为泥岩厚度，单位：m；C_m为泥岩压实系数，利用去压实校正方法计算获得。

3.根据权利要求1所述的利用旋回分析法恢复不整合剥蚀厚度的方法，其特征在于：所述地质资料包括地层分层数据、地层年代表，以及录井资料或孔隙度资料。

4.根据权利要求3所述的利用旋回分析法恢复不整合剥蚀厚度的方法，其特征在于：所述去压实校正方法为经验公式法或孔隙度-岩性方法；当所述地质资料包括录井资料时，所述去压实校正方法采用经验公式法；当所述地质资料包括孔隙度资料时，所述去压实校正方法采用孔隙度-岩性方法。

5.根据权利要求1所述的利用旋回分析法恢复不整合剥蚀厚度的方法，其特征在于：在计算不整合剥蚀厚度步骤后，还包括利用地层分层数据或现有恢复不整合剥蚀厚度的方法验证所得不整合剥蚀厚度准确性的步骤。

6.根据权利要求5所述的利用旋回分析法恢复不整合剥蚀厚度的方法，其特征在于，利用地层分层数据进行准确性验证的具体步骤为：利用恢复的完整沉积速率曲线，通过对未发生剥蚀的地层对应的沉积速率曲线段进行积分得到未剥蚀地层的原始厚度，通过压实计算得到未剥蚀地层的压实后地层厚度，与地层分层数据中的现今地层厚度进行对比，验证准确性。

7.根据权利要求5所述的利用旋回分析法恢复不整合剥蚀厚度的方法，其特征在于：所述现有恢复不整合剥蚀厚度的方法选自泥岩声波时差法、镜质体反射率法、磷灰石裂变径迹法、物质平衡法、包裹体均一温度法、沉积速率法、砂岩孔隙度法、波动分析法、地震趋势分析法、邻层厚度比值法和参考层厚度变化率法中的任一种。