[發明專利]基于EEMD的高分辨率層序地層劃分方法

在審

• 申請號: 201710156428.9
• 申請日: 2017-03-16
• 公開（公告）號: CN106707341A
• 公開（公告）日: 2017-05-24
• 發明（設計）人: 楊巍;楊應;朱仕軍;張桓
• 申請（專利權）人: 西南石油大學
• 主分類號: G01V1/30
• 分類號: G01V1/30;G06F17/14
• 代理公司: 成都方圓聿聯專利代理事務所(普通合伙)51241
• 代理人: 周佳

1.基于EEMD的高分辨率層序地層劃分方法，其特征在于，包括以下步驟：

1)、輸入測井數據，進行EEMD經驗模態分解，得到IMF分量，將待處理測井信號中加入等長度不等幅高斯白噪聲，對復合信號EMD分解，重復操作k次，得到IMF分量c_ik和余項r_ik，白噪聲幅值遵循如下規律：

a＝ek或ln e+0.5a ln k＝0

式中，e為原信號標準差，即信號與EMD重構結果的偏離，a為白噪聲幅值；

2)、對IMF分量整體平均：

${\mathrm{IMF}}_{ik}=\frac{1}{n}\underset{k=1}{\overset{n}{\Σ}}{\mathrm{IMF}}_{ik}$

$r_i=\frac{1}{n}\underset{k=1}{\overset{n}{\Σ}}{r}_{ik}$

式中，n為EMD分解的次數；

3)、對輸出的IMF分量進行信噪比SNR計算，并去除掉低信噪比的IMF第一分量；對從IMF2分量開始，分別對個IMF進行希爾伯特變換，并得到HHT時頻譜；

4)、地質層位與各個IMF得到HHT時頻譜進行匹配，并標定出地質層位對應的時頻譜突變點；

5)、在進行深度匹配后的HHT時頻譜的基礎上，根據工區精度需要，進行各級別的層序劃分，劃分準則為：HHT時頻譜橫坐標為頻率，縱坐標同測井數據一致為深度；在HHT時頻譜上，隨著深度變化，頻率升高即為海/湖平面降低；隨著深度變化，頻率降低即為海/湖平面升高。