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

说明书

本发明提供基于米兰科维奇旋回的古水深定量恢复方法，包括选取标定井和岩样，获取岩样岩相敏感测井数据，并进行微常量元素分析化验计算单点古水深；对岩相敏感测井数据降噪，将降噪后的测井数据分别进行时间域和深度域频谱分析及滤波分析，确认降噪测井数据是否符合米兰科维奇天文信号，然后通过演化谐波功率谱分析，并通过蒙特卡洛模拟计算演化平均光谱失配，获得降噪测井数据的深度域谱图；基于米兰科维奇天文旋回理论，追踪零假设显著性水平最小值，并将从中获取的沉积速率与单点古水深等深度对应，拟合沉积速率‑古水深公式，校验公式后计算标定井的全序列古水深。该方法更直观地表征不同构造‑沉积单元在不同地质时期的沉积环境及演化序列。

技术领域

本发明涉及地球化学和地球物理中古水深恢复技术领域，具体涉及一种基于米兰科维奇旋回的古水深定量恢复方法和存储介质。

背景技术

现有技术中，在文献(杨万芹,朱德顺,银燕,丁桔红,2015.古水深的地球化学恢复方法及在层序地层划分中的应用.地质论评,61(S1):756-757.)中，根据Fe/Co、Fe/Mn和Th/U三个元素比值综合，在结合地区平均古水深进行古水深曲线恢复。此方法主要建立对古水深变化反映灵敏的地球化学指标，进行多个指标综合，进行古水深的量化恢复，并应用恢复的古水深曲线。但是此方法无法对无取心或无微量元素测试的井段进行水深的恢复。

在另一篇专利文件(申请公布号CN 106019401 A)“一种古水深的定量恢复方法以及装置”中，确定古湖盆的研究层位中滩坝发育区域；确定所述滩坝发育区域中至少一口典型井的单点古水深；确定所述滩坝发育区域中的至少一口所述典型井中的Th含量和U含量，并根据所述Th含量与所述U含量计算Th与U比值；根据所有典型井的单点古水深以及单点古水深所在位置所对应的Th与U比值，拟合单点古水深与Th/U值的关系方程，恢复所述典型井的定量古水深；根据单点古水深与Th/U值的关系方程，恢复所述古湖盆中非滩坝发育区域其它井的定量古水深。通过该方法能够获取古湖盆同一沉积时期水体深度垂向和平面变化的连续定量古水深数据。但是此方法限于滩坝沉积体系特殊性，无法对无滩坝发育区域建立古水深与Th/U值的关系方程；此外，表征古水深的关键参数选取及恢复路径相对独立。

发明内容

本发明解决的一个主要问题是现有古水深量化恢复方法无法更直观地恢复沉积时期的全序列的水体深度，从而难以建立沉积时期古地貌和古环境的问题。

根据本发明的一个方面，本发明提供一种基于米兰科维奇旋回的古水深定量恢复方法，包括：

选取标定井；

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

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

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

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

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

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

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

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