[發明專利]基于米蘭科維奇旋回的古水深定量恢復方法和存儲介質

在審

• 申請號: 202110947608.5
• 申請日: 2021-08-18
• 公開（公告）號: CN113687440A
• 公開（公告）日: 2021-11-23
• 發明（設計）人: 劉強虎;朱紅濤;李志垚;賈悅銳;鄭小波;鮑怡晨
• 申請（專利權）人: 中國地質大學（武漢）
• 主分類號: G01V9/00
• 分類號: G01V9/00
• 代理公司: 武漢知產時代知識產權代理有限公司 42238
• 代理人: 萬文廣

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任一一項所述的一種基于米蘭科維奇旋回的古水深定量恢復方法。