[发明专利]一种基于沉积速率估算的三角洲沉积微相定量划分方法

权利要求书

1.一种基于沉积速率估算的三角洲沉积微相定量划分方法，其特征在于，具体过程为：

步骤S1：采集地层测井数据序列，地层测井数据序列包括自然伽马GR和自然电位SP数据；

步骤S2：数据预处理，将步骤S1的地层测井数据序列进行重采样、剔除异常值、去均值、去趋势化以及预白化处理，以消除地层测井数据中的各种环境“噪声”；

步骤S3：基于eCOCO的地层沉积速率估算，该步骤包括相关系数评价、天文驱动检验和有贡献的天文周期参数个数统计三个环节，其中，相关系数评价即计算待测试沉积速率下的地层时间序列与天文目标序列的相关系数；天文驱动检验即通过零假设检验某沉积速率下是否存在天文周期信号；统计对相关系数有贡献的天文周期参数的个数用于避免仅有少数天文周期参数得到高相关系数的情况，待测试的沉积速率在一定取值范围内按特定步长取值，同时满足相关系数高、零假设p值低、有贡献的天文周期参数数量多三个条件的沉积速率即为最优沉积速率；待测试沉积速率的取值范围受Nyquist频率和Reyleigh频率的约束，其中，时间域内周期图的Nyquist频率f_nyq和Reyleigh频率f_ray分别定义为：

f_nyq＝sr/(2*Δd)

f_ray＝sr/(M*Δd)

式中，sr为沉积速率，单位m/ka，Δd为采样间距，单位m，M为地层测井数据序列的总样本数量，根据f_nyq大于等于最高天文频率确定sr最小值，根据f_ray小于等于最低天文频率确定sr最大值；

待测试沉积速率下的地层时间序列功率谱与天文目标序列功率谱之间的相关系数由皮尔逊积矩相关系数ρ获得：

式中，T和D分别为天文目标序列和地层时间序列的功率谱，N为T或D的观测值数量，μ_T、μ_D分别为T和D的均值，σ_T、σ_D分别为T和D的标准差，功率谱通过周期图法获得，利用一阶自回归模型对地层转换后的时间序列进行红噪估计，若功率谱幅值小于红噪均值则赋值为0；

天文驱动检验的具体过程为：建立零假设H₀，即假设地层时间序列中不存在天文周期信号且所有超过一阶自回归模型红噪背景值的n个峰值f1，f2，...，fn均是偶然产生的，在每个待测试的沉积速率下，利用蒙特卡洛模拟法生成随机数据功率谱和目标功率谱之间相关系数ρ值的分布，建立零分布，从而估计该沉积速率下随机产生相应相关系数的概率p，如果检验的p值＜0.01，认为假设检验是显著的，拒绝零假设H₀，即在该沉积速率下存在天文周期信号的置信度为1-p；

步骤S4：生成沉积速率曲线，基于滑动窗口技术，eCOCO跟踪地层不同深度的最优沉积速率，获得地层剖面深度域的沉积速率曲线；

步骤S5：沉积微相划分，以步骤S4获得的沉积速率曲线为基础，寻找沉积速率显著变化的界线，结合沉积背景及岩性组合地质信息完成三角洲沉积环境下微相的划分。

2.根据权利要求1所述基于沉积速率估算的三角洲沉积微相定量划分方法，其特征在于，所述S3中待测试沉积速率取值范围参照已知地层剖面两个以上层位的地质年龄而估算的平均沉积速率；对于缺少年龄约束的剖面，待测试沉积速率取值范围参照相似沉积系统的已知沉积速率。

3.根据权利要求2所述基于沉积速率估算的三角洲沉积微相定量划分方法，其特征在于，所述S3中天文周期参数根据地层的大致年龄确定。

4.根据权利要求3所述基于沉积速率估算的三角洲沉积微相定量划分方法，其特征在于，所述S4中滑动窗口技术中滑动窗口的长度大于1个地层长偏心率周期对应的沉积厚度，步距为采样间距的整数倍。