[發(fā)明專利]儲(chǔ)層裂縫顯示方法及裝置

有效

• 申請?zhí)枺?/td>: 201610887204.0
• 申請日: 2016-10-11
• 公開（公告）號: CN106569271B
• 公開（公告）日: 2018-09-04
• 發(fā)明（設(shè)計(jì)）人: 謝春輝;雍學(xué)善;楊午陽;楊慶;周春雷;王恩利;鄢高韓;杜炳毅;李琳
• 申請（專利權(quán)）人: 中國石油天然氣股份有限公司
• 主分類號: G01V1/30
• 分類號: G01V1/30;G01V1/34
• 代理公司: 北京三友知識產(chǎn)權(quán)代理有限公司 11127
• 代理人: 黨曉林

1.一種儲(chǔ)層裂縫顯示方法，其特征在于，所述儲(chǔ)層裂縫顯示方法包括以下步驟：

獲取目標(biāo)儲(chǔ)層的具有三維參數(shù)的地震資料；

基于所述地震資料獲得所述目標(biāo)儲(chǔ)層的裂縫發(fā)育特征，所述裂縫發(fā)育特征至少包括裂縫密度矩陣、裂縫發(fā)育方向矩陣和裂縫縱橫比矩陣；

基于所述裂縫密度矩陣、所述裂縫發(fā)育方向矩陣和所述裂縫縱橫比矩陣得到裂縫延伸長度和張開度，其具體的計(jì)算過程如下：

ρ_max＝max(P_ni×nx)

ρ_o＝P(i,j)

r_s＝r_lα

φ_o＝Φ_ni×nx(i,j)

其中，r_l表示裂縫延伸長度，r_s表示張開度，ρ_max表示裂縫密度的最大值點(diǎn)，r_max表示根據(jù)地質(zhì)調(diào)查和測井裂縫研究得到裂縫密度最大值點(diǎn)的裂縫延伸長度，α表示裂縫縱橫比矩陣中的縱橫比，ρ_O表示裂縫密度矩陣中某一點(diǎn)的裂縫密度，ρ_O中O點(diǎn)對應(yīng)的橢圓盤長軸方向與裂縫發(fā)育方向相同，P_ni×nx表示多個(gè)點(diǎn)的裂縫密度，ni為目標(biāo)儲(chǔ)層工區(qū)內(nèi)inline線的數(shù)目，nx為目標(biāo)儲(chǔ)層工區(qū)內(nèi)xline線的數(shù)目，i，j為點(diǎn)O對應(yīng)的線道號，i＝1,2,...,ni，j＝1,2,...,nx，Φ_ni×nx(i,j)表示目標(biāo)儲(chǔ)層工區(qū)內(nèi)不同的點(diǎn)的橢圓盤長軸方向，φ_o表示裂縫發(fā)育方向矩陣中的裂縫發(fā)育方向；

基于所述裂縫延伸長度和所述張開度得到裂縫延伸權(quán)重矩陣，所述裂縫延伸權(quán)重矩陣內(nèi)的非零值呈橢圓盤形狀，該橢圓盤的尺寸參數(shù)由裂縫延伸長度和張開度確定，所述裂縫延伸權(quán)重矩陣中包括各點(diǎn)的權(quán)值，權(quán)值的具體計(jì)算過程如下：

γ_M＝d_M/r_M

ω_M＝f(γ_M)

f(γ_M＝0)＝1

f(γ_M≥1)＝0

其中，ω_M表示權(quán)值，權(quán)值是該權(quán)值位置的點(diǎn)與橢圓盤中心O之間的距離d_M與橢圓盤該方向極徑r_M比值γ_M的函數(shù)，橢圓盤內(nèi)各點(diǎn)權(quán)值從中心等于1遞減至邊界處等于0，權(quán)重等值線為橢圓；

通過所述裂縫密度矩陣中各點(diǎn)的裂縫密度與所述裂縫延伸權(quán)重矩陣得到新的裂縫密度。

2.根據(jù)權(quán)利要求1所述的儲(chǔ)層裂縫顯示方法，其特征在于，在所述基于所述地震資料獲得所述目標(biāo)儲(chǔ)層的裂縫發(fā)育特征的步驟中，具體為對所述地震資料進(jìn)行各向異性分析進(jìn)而獲得所述目標(biāo)儲(chǔ)層的裂縫發(fā)育特征。

3.根據(jù)權(quán)利要求1所述的儲(chǔ)層裂縫顯示方法，其特征在于，所述裂縫密度矩陣中至少包括各點(diǎn)的裂縫密度。

4.根據(jù)權(quán)利要求1所述的儲(chǔ)層裂縫顯示方法，其特征在于，在所述通過所述裂縫密度矩陣中各點(diǎn)的裂縫密度與所述裂縫延伸權(quán)重矩陣得到新的裂縫密度的步驟中，將所述裂縫密度矩陣中各點(diǎn)的裂縫密度與所述裂縫延伸權(quán)重矩陣的乘積進(jìn)行疊加以對儲(chǔ)層內(nèi)每點(diǎn)的裂縫進(jìn)行延伸擴(kuò)展進(jìn)而生成新的裂縫密度，其具體計(jì)算公式如下：

或

或

其中，W_O表示裂縫延伸權(quán)重矩陣，ρ_O表示裂縫密度矩陣中各點(diǎn)的裂縫密度，P_ni×nx表示裂縫密度，ni為工區(qū)內(nèi)inline線的數(shù)目，nx為工區(qū)內(nèi)xline線的數(shù)目。