[發明專利]一種基于分頻AVO反演的儲層含氣性半定量識別方法

有效

• 申請號: 201310439924.7
• 申請日: 2013-09-24
• 公開（公告）號: CN104459771B
• 公開（公告）日: 2017-08-18
• 發明（設計）人: 肖鵬飛;王世星;呂慧;張亞紅
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司石油物探技術研究院
• 主分類號: G01V1/28
• 分類號: G01V1/28
• 代理公司: 北京思創畢升專利事務所11218
• 代理人: 郭韞

1.一種基于分頻AVO反演的儲層含氣性半定量識別方法，其特征在于：所述方法包括：

(1)利用三維地震數據，分別計算縱波速度衰減梯度和橫波速度衰減梯度；分別繪制縱波速度衰減梯度與生產井產能之間的關系圖和橫波速度衰減梯度與生產井產能之間的關系圖，開展線性擬合，計算兩個衰減梯度與產能之間的相關系數R²；縱波速度衰減梯度與橫波速度衰減梯度之間的權重值x的初始值選擇為1,即x₀＝1；利用地球物理反演技術對公式(1)開展反演，求得縱波速度衰減梯度和橫波速度衰減梯度：

R(θ,f)≈A(θ)ΔVpVp(f0)+(f-f0)×A(θ)ddf(ΔVpVp)(f)+B(θ)ΔVsVs(f0)+(f-f0)B(θ)ddf(ΔVsVs)(f)---(1)]]>

其中，R(θ,f)表示地震反射系數；θ,f分別是地震反射系數的入射角與頻率；A(θ),B(θ)分別是只與入射角相關的系數函數；V_p,V_s,ΔV_p,ΔV_s分別表示地震縱波速度、橫波速度、地震縱波速度變化、橫波速度變化；分別表示是頻率f的函數；

(2)沿各井提取縱波速度衰減梯度和橫波速度衰減梯度，計算各井生產層段反演值的加權和g_i,0，其中i表示井號；所述步驟(2)中是利用公式(5)計算各井生產層段反演值的加權和g_i,0：

gi,0=aix0+bi,ai=ddf(ΔVpVp)(f),bi=ddf(ΔVsVs)(f)---(5),]]>

每改變一次x_k，利用公式(5)得到一個對應的g_i,k；

(3)以井產能Y為縱坐標、以G為橫坐標，進行最小二乘線性擬合，形成如下擬合公式：

Y＝CG+E (2)

其中C、E為常數；G是g_i,0的大寫表示；

(4)計算各井點產能與步驟(3)得到的擬合線之間的距離和d₀；

所述步驟(4)所用的計算公式為：

d0=Σi=1l|Cgi-yi+E|C2+1=Σi=1l|C(aix0+bi)-yi+E|C2+1---(3)]]>

其中d₀為距離和、y_i為第i井的產能，C、E為最小二乘線性擬合常數；每改變一次x_k，利用公式(3)得到一個對應的d_k；

(5)令x_k+1＝x_k+t，其中t為步長，k表示迭代次數；

(6)重復步驟(2)至步驟(5)，直到d_k>d_k-1，轉入步驟(7)；

(7)輸出g_i,k，g_i,k＝a_ix_k+b_i (4)；

根據反演結果數值進行含氣性識別，反演結果數值越高則含氣量越大。