[发明专利]基于米兰科维奇旋回的古水深定量恢复方法和存储介质

权利要求书

1.一种基于米兰科维奇旋回的古水深定量恢复方法，其特征在于，包括：

选取标定井；

基于所述标定井选取岩样，并获取所述岩样的岩相敏感测井数据；

对所述岩样进行微常量元素分析化验，并利用微常量元素丰度计算获得单点古水深；

对所述岩相敏感测井数据进行降噪，获得降噪测井数据；

将所述降噪测井数据分别进行时间域和深度域频谱分析和滤波分析，确认所述降噪测井数据是否符合米兰科维奇天文信号；

所述降噪测井数据符合所述米兰科维奇天文信号后，通过演化谐波功率谱分析，并通过蒙特卡洛模拟计算演化平均光谱失配，获得所述降噪测井数据的深度域谱图；

基于米兰科维奇天文旋回理论，追踪零假设显著性水平最小值，获取所述深度域谱图中的沉积速率；

将获取的所述沉积速率与所述单点古水深等深度对应，建立并拟合沉积速率-古水深公式；

对所述沉积速率-古水深公式进行校验，获得校验后的沉积速率-古水深公式；

使用校验后的所述沉积速率-古水深公式计算所述标定井的全序列古水深，恢复目标区域的古水深时空分布。

2.如权利要求1所述的一种基于米兰科维奇旋回的古水深定量恢复方法，其特征在于，所述对所述沉积速率-古水深公式进行校验，获得校验后的沉积速率-古水深公式还包括：

使用微常量元素丰度计算获得的所述单点古水深与利用米兰科维奇天文旋回公式计算的古水深进行比对，验证所述沉积速率-古水深公式的准确性。

3.如权利要求1所述的一种基于米兰科维奇旋回的古水深定量恢复方法，其特征在于，所述对所述岩相敏感测井数据进行降噪还包括：

利用dmey函数将所述岩相敏感测井数据分解为9层，去除分层后的所述岩相敏感测井数据中的背景a9和最高频d1；

并将去除背景a9和最高频d1后所述岩相敏感测井数据保存为含有深度和表头的深度域数据以及不含所述深度和所述表头的时间域数据。

4.如权利要求3所述的一种基于米兰科维奇旋回的古水深定量恢复方法，其特征在于，将所述降噪测井数据分别进行时间域和深度域频谱分析以及滤波分析，确认所述降噪测井数据是否符合米兰科维奇天文信号还包括：

导入时间域数据，进行频谱分析验证特征峰频率是否满足米氏旋回周期的反比；

导入深度域数据，使用“astrochron”软件包中的mtm代码和bandpass代码验证所述深度域数据是否符合米兰科维奇天文周期。

5.如权利要求1所述的一种基于米兰科维奇旋回的古水深定量恢复方法，其特征在于，所述深度域谱图还包括沉积速率、天文周期和零假设显著性水平信息。

6.如权利要求1所述的一种基于米兰科维奇旋回的古水深定量恢复方法，其特征在于，所述微常量元素为Co元素。

7.如权利要求1所述的一种基于米兰科维奇旋回的古水深定量恢复方法，其特征在于，所述使用校验后的所述沉积速率-古水深公式计算所述标定井的全序列古水深包括：

利用校验后的所述沉积速率-古水深公式计算所有标定井的全序列古水深，进而形成井网外推全目标区域的古水深水深分布；

基于所述全目标区域的古水深水深分布，定量恢复全目标区域的各时期的古水深分布。

8.如权利要求1所述的一种基于米兰科维奇旋回的古水深定量恢复方法，其特征在于，所述建立并拟合沉积速率-古水深公式包括：

将微常量元素丰度计算获得的所述单点古水深与对应深度的所述沉积速率整合到一起，利用线性回归软件拟合所述单点古水深与对应深度的所述沉积速率间的函数关系。

9.一种存储介质，其特征在于，所述存储介质为计算机可读存储介质，所述计算机可读存储介质上存储有如权利要求1-8任一一项所述的一种基于米兰科维奇旋回的古水深定量恢复方法。