[發(fā)明專利]一種時間域地震層速度反演方法

有效

• 申請?zhí)枺?/td>: 201210246809.3
• 申請日: 2012-07-17
• 公開（公告）號: CN103543466A
• 公開（公告）日: 2014-01-29
• 發(fā)明（設(shè)計）人: 吳鵬;孫開峰
• 申請（專利權(quán)）人: 中國石油化工股份有限公司;中國石油化工股份有限公司石油物探技術(shù)研究院
• 主分類號: G01V1/28
• 分類號: G01V1/28
• 代理公司: 北京思創(chuàng)畢升專利事務(wù)所 11218
• 代理人: 郭韞

1.一種時間域地震層速度反演方法，其特征在于：所述方法首先通過地震測井獲取初始瞬時速度，然后建立瞬時速度與均方根速度之間的關(guān)系，將所述初始瞬時速度轉(zhuǎn)換成均方根速度，再用實際拾取的均方根速度進(jìn)行約束構(gòu)建目標(biāo)函數(shù)，然后求解目標(biāo)函數(shù)得到更新速度模型所需的變量，通過多次迭代后，當(dāng)目標(biāo)函數(shù)的值小于閥值ε時，得到最終的層速度。

2.根據(jù)權(quán)利要求1所述的時間域地震層速度反演方法，其特征在于：所述方法包括以下步驟：

第一步：獲取初始瞬時速度：

將測井速度或平均瞬時速度作為初始瞬時速度；

第二步：構(gòu)建速度更新方程如下：

V_n(t)＝V_n-1+ΔV*e^-kt?????????????(1)

(1)式中，V_n(t)是待求解的層速度，ΔV是速度差，n代表第n次迭代，k為速度變化梯度，k的初始值設(shè)為零，*表示相乘，e為自然指數(shù)；

第三步：求取均方根速度，計算公式為：

Vrmscal(ΔV,k,t)=1t·∫0tVn2(τ)dτ=1t·∫0t(Vn2+2VnΔV·e-2kτ+ΔV2·e-2kt)dτ---(2)]]>

(2)式中，n的初始值為1，當(dāng)n＝1時，設(shè)ΔV＝0，k＝0，則V₁＝V₀，V₀為第一步得到的初始瞬時速度；

第四步：構(gòu)建目標(biāo)函數(shù)如下：

Φ(ΔV,k)=Σi=1I[Vrmscal(ΔV,k,i)-Vidata]2→Min---(3)]]>

(3)式中，I為速度的時間采樣點數(shù)，為初始拾取的均方根速度；通過求解所述目標(biāo)函數(shù)求得變量ΔV和k；

第五步，將第四步求得的ΔV和k代入公式(1)，得到新的層速度，同時判斷目標(biāo)函數(shù)值是否小于閥值ε，如果是，則所述新的層速度即為最終的層速度數(shù)據(jù)體，然后轉(zhuǎn)入第六步，否則，使n＝n+1，然后返回第二步；

第六步，輸出所述層速度數(shù)據(jù)體。