[發明專利]點約束貝葉斯稀疏脈沖反演方法

無效

• 申請號: 201210417846.6
• 申請日: 2012-10-26
• 公開（公告）號: CN103792571A
• 公開（公告）日: 2014-05-14
• 發明（設計）人: 徐希坤;孫喜新;王志高;穆星;楊培杰;王長江;劉惠民;王永詩;羅紅梅;畢俊鳳;鄧玉珍;劉華夏;王東曄;方旭慶;石砥石
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司勝利油田分公司地質科學研究院
• 主分類號: G01V1/28
• 分類號: G01V1/28;G01V1/40
• 代理公司: 濟南日新專利代理事務所 37224
• 代理人: 謝省法

1.點約束貝葉斯稀疏脈沖反演方法，其特征在于，該點約束貝葉斯稀疏脈沖反演方法包括:

步驟1，提取地震子波；

步驟2，建立初始地質模型；

步驟3，建立反演目標函數；

步驟4，對該反演目標函數各項求導數，得到非線性反演方程；

步驟5，設置反演參數；

步驟6，用迭代法求解該非線性反演方程；以及

步驟7，輸出儲層高分辨率反演結果。

2.根據權利要求1所述的點約束貝葉斯稀疏脈沖反演方法，其特征在于，在步驟1中，采用確定性子波提取方法，利用聲波測井和密度測井資料，計算出反射系數序列，并結合井旁地震道由褶積模型求出該地震子波。

3.根據權利要求1所述的點約束貝葉斯稀疏脈沖反演方法，其特征在于，在步驟1中，采用統計性子波提取方法，假設地震子波是時不變的，地下的反射系數是高斯白噪聲，并對子波的相位譜給出一定的假設，以通過地震道自身來估計該地震子波。

4.根據權利要求1所述的點約束貝葉斯稀疏脈沖反演方法，其特征在于，在步驟2中，在該初始地質模型建立之前，分析斷層、層位接觸關系、地層沉積模式，再以測井資料為基礎，在已有解釋層位的基礎上建立反演用的該初始地質模型，并在該初始地質模型生成后，檢查該初始地質模型生成時產生的層位斷層的空間合理性，以確保后續反演結果的可靠性。

5.根據權利要求1所述的點約束貝葉斯稀疏脈沖反演方法，其特征在于，步驟3包括：

進行基于貝葉斯理論的參數估計；

通過柯西約束體現反演過程的隨機性和稀疏程度；

通過波阻抗點約束來提高反演結果的準確性和穩定性；以及

以貝葉斯理論為基礎，結合反射系數柯西約束和地層模型點約束建立該反演目標函數。

6.根據權利要求5所述的點約束貝葉斯稀疏脈沖反演方法，其特征在于，在進行基于貝葉斯理論的參數估計的步驟中，以貝葉斯理論為基礎，通過貝葉斯公式將先驗信息和樣本信息相結合，得到參數的后驗信息，再進行基于貝葉斯理論的參數估計。

7.根據權利要求6所述的點約束貝葉斯稀疏脈沖反演方法，其特征在于，在進行基于貝葉斯理論的參數估計的步驟中，假定要估計的參數為θ，利用已有的信息來設定該估計的參數θ的先驗分布p(θ)，再將得到的該樣本信息融入該先驗信息來改進先驗分布，假定與該估計的參數θ有關的觀測樣本x，用p(x|θ)表示其條件密度函數，得到該估計的參數θ的后驗分布的公式為：?

其中，p(θ|x)為后驗概率密度，p(θ)為先驗概率密度，p(x|θ)為似然函數，由該后驗分布對該估計的參數θ作出推斷，進行期望、方差等估計。?

8.根據權利要求6所述的點約束貝葉斯稀疏脈沖反演方法，其特征在于，在通過柯西約束體現反演過程的隨機性和稀疏程度的步驟中，以零值為主，有少數非零值的稀疏解柯西公式為：?

其中，m_i表示待反演的參數，表示參數的均值，σ_m表示參數的方差。?

9.根據權利要求6所述的點約束貝葉斯稀疏脈沖反演方法，其特?征在于，在以貝葉斯理論為基礎，結合反射系數柯西約束和地層模型點約束建立該反演目標函數中，得到的該反演目標函數為：?

其中，m為待反演的參數，G為地震子波矩陣，d為實際地震道，R為稀疏約束項，C為單位矩陣，ξ為地層的模型，μ稀疏參數，α穩定性參數。?