[發明專利]處理地震剖面圖像的方法和裝置

有效

• 申請號: 201610192345.0
• 申請日: 2016-03-30
• 公開（公告）號: CN107290786B
• 公開（公告）日: 2019-01-08
• 發明（設計）人: 劉旭躍;周巍;張兵
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司石油物探技術研究院
• 主分類號: G01V1/36
• 分類號: G01V1/36;G06T5/00
• 代理公司: 北京思創畢升專利事務所 11218
• 代理人: 孫向民

1.一種處理地震剖面圖像的方法，該方法包括：

基于下式對地震剖面圖像進行TV去噪處理：

其中，所述地震剖面圖像共有N條線，每條線具有N個cdp點，表示第n次迭代得到的第i條線中第j個cdp點的像素值，表示原始地震資料中第i條線中第j個cdp點的像素值，Δt表示時間步長，λ為預設參數，邊界條件滿足

利用相似性學習法對經TV去噪處理后的地震剖面圖像進行插值，得到處理后的地震剖面圖像；

其中，所述相似性學習法為逐層相似性學習法，所述逐層相似性學習法包括：

第一步，建立數學模型

圍繞未知像素點選擇鄰域窗口B(ε,p)，其中，ε為該鄰域窗口的半徑，p為該鄰域窗口中的像素點的數目；

針對鄰域窗口B(ε,p)中的每一個像素點形成對應于未知像素點的訓練窗口,其中針對鄰域窗口B(ε,p)中存在的未知像素點，定義對應的訓練窗口中包含的像素點與鄰域窗口在方向上一一對應，并從訓練窗口的各個子模式拷貝學習以得到邊緣導向屬性；

基于下面的公式得到權值β：

其中，β_{_tmp}是1*p的行向量，X∈B(ε,p)，X為p*1的列向量，Y為該未知像素點的像素值，Q為n*1階向量，Q中的第k個元素q_k為：

其中C_k是未知像素點到已知像素點的距離，C_k為p*1的列向量，Δ是用來避免分母為0的預設參數；

第二步，填補未知點

采用逐層遞進的分組行進法進行填補，逐層遞進的分組行進法是找到滿足到達的時間的一組數據B：

其中，α是權重因子，表示曲線到達點的時刻，F表示在到達點處的速度，T_m表示第m個的時間，T_L表示第L組的時間，ΔT_L表示第L組的時間步長，B_L表示待填補層；

該組數據B對應于修補邊界層，對此數據同時進行更新，每層數據的更新只受上層數據或已知數據的影響，和本層已求出的數值無關。

2.根據權利要求1所述的方法，其中，基于下式得到公式1中的

其中，L(*)為拉普拉斯算子，

3.根據權利要求1所述的方法，其中，公式1中的λ是基于下列方法確定的：

設其中，M表示總的迭代次數，z表示當前的圖像函數，K表示降質函數并取歸一化的值，表示[K^TK+λ_tmpL(x)]^-1K^T，x表示圖像上任一像素點的像素值，trace(*)表示求矩陣的跡，||*||₂表示求范數的組合數；

取是f(λ)帽，指統計參數的估計量。