[發明專利]儲層孔隙度、含水飽和度和泥質含量參數同時反演新方法

有效

• 申請號: 201610191787.3
• 申請日: 2016-03-30
• 公開（公告）號: CN107290782B
• 公開（公告）日: 2019-07-12
• 發明（設計）人: 楊培杰;董立生;邵卓娜;李逢英;石峪;顏世翠;徐仁;亓雪靜;鄭文召;初春光;鄒文勇;范云
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司勝利油田分公司勘探開發研究院
• 主分類號: G01V1/30
• 分類號: G01V1/30
• 代理公司: 濟南日新專利代理事務所 37224
• 代理人: 崔曉艷

1.儲層孔隙度、含水飽和度和泥質含量參數同時反演新方法，其特征在于，該儲層孔隙度、含水飽和度和泥質含量參數同時反演新方法包括：

步驟1，采集巖石彈性參數、測井數據和巖心數據信息；

步驟2，將巖石彈性參數與儲層物性參數相結合，建立巖石物理模型；

步驟3，建立儲層孔隙度、含水飽和度和泥質含量參數三參數反演的目標函數，并求解目標函數；

步驟4，通過求解目標函數，輸出儲層孔隙度、含水飽和度和泥質含量參數反演結果；

在步驟1中，通過疊前反演獲得巖石彈性參數，疊前反演是通過分析疊前道集來研究地下巖性的地震分析方法，借助于Zoeppritz方程或其近似式，最終反演出巖石彈性參數，反演出的巖石彈性參數包括巖石的體積模量、剪切模量和密度信息；

在步驟2中，首先，將巖石彈性參數與儲層物性參數相結合，建立確定性的巖石物理模型；其次，針對不同儲集類型的儲層，以確定性巖石物理模型為基礎，通過統計信息，建立統計的巖石物理模型，統計巖石物理模型用來描述測井上沒有體現的所有可能的儲層條件實現；所建立的巖石物理模型如公式(1)所示：

[K,μ,ρ]＝f_RPM(φ,s_w,C)+ε (1)

其中，f_RPM表示確定性巖石物理模型，K,μ,ρ表示巖石彈性參數，包括體積模量、剪切模量、密度，φ,s_w,C表示儲層物性參數，包括儲層孔隙度、含水飽和度、泥質含量，ε表示實際觀測數據的誤差；

在步驟3中，根據貝葉斯算法建立物性參數提取的目標函數，目標函數為儲層物性參數在已知彈性參數條件下的最大后驗概率分布，表示為：

[φ,sw,C]＝argMax{P([φ,sw,C]_j|[K,μ,ρ])} (3)

其中，φ,s_w,C分別為儲層孔隙度、含水飽和度、泥質含量，P()表示概率密度，argMax表示求最大值,K表示體積模量，μ表示剪切模量，ρ表示密度，j表示各個儲層物性參數的類別數；

P([φ,sw,C]_j|[K,μ,ρ])為儲層物性參數的后驗分布，該公式的意義是，已知彈性參數，分別計算儲層孔隙度、含水飽和度和泥質含量分屬于各個儲層物性參數類的后驗概率分布；當后驗概率取最大值時，該儲層物性參數即為最終反演結果；

在步驟3中，通過遺傳算法和蒙特卡羅方法相結合求解目標函數，首先生成初始化群體，并計算每個個體的適應度函數，如果結果滿意，則遺傳算法結束，如果結果不滿意，更新初始化的群體，通過蒙特卡羅進行仿真，并再次計算每個個體的適應度函數，更新群體，如此循環，直到計算結果滿意為止。

2.根據權利要求1所述的儲層孔隙度、含水飽和度和泥質含量參數同時反演新方法，其特征在于，在步驟1中，測井數據包括儲層孔隙度、含水飽和度和泥質含量信息；巖心數據提供測井數據的輔助信息，包括巖石類型、巖石孔隙度。