[發明專利]一種確定煤層氣縱波和橫波速度的方法

有效

• 申請號: 201310449931.5
• 申請日: 2013-09-27
• 公開（公告）號: CN104516015B
• 公開（公告）日: 2017-03-15
• 發明（設計）人: 馬光凱;孫鵬遠;錢忠平;任曉喬;黃燕
• 申請（專利權）人: 中國石油天然氣集團公司;中國石油集團東方地球物理勘探有限責任公司
• 主分類號: G01V1/28
• 分類號: G01V1/28
• 代理公司: 北京康信知識產權代理有限責任公司11240
• 代理人: 吳貴明,張永明

1.一種確定煤層氣縱波和橫波速度的方法，特點是通過以下步驟實現：

1）采集地震勘探和測井數據，處理后得到縱波速度x和橫波速度y；

2）按照下式分別對縱波速度x和橫波速度y進行處理，使它們的平均值分別為零：

xj=xj-(Σj=1Lxj)/Lyj=yj-(Σj=1Lyj)/L---(1)]]>

式中：L是組數，j是縱波速度x和橫波速度y的組序列號，x_j和y_j分別是縱波速度x和橫波速度y的第j個元素；

3）按照下式對L組縱波速度x和橫波速度y進行歸一化處理：

xj=xj/max(x)yj=yj/max(y)(j=1,2,···,L)---(2)]]>

式中，L是組數；j是縱波速度x和橫波速度y的組序列號；x_j和y_j分別是是縱波速度x和橫波速度y的第j個元素；max(x)是求取縱波速度x的最大值；max(y)是求取橫波速度y的最大值；

4）按下式分別計算縱波速度x和橫波速度y最佳變換的初始值φ⁰(x)和θ⁰(y)：

φ0(xj)=0θ0(yj)=yj(j=1,2,···,L)---(3)]]>

式中，L是縱波速度x和橫波速度y的組數；j是縱波速度x和橫波速度y的組序列號；x_j和y_j分別是縱波速度x和橫波速度y的第j個元素；

5）按照下式求取縱波速度x和橫波速度y的最佳變換φ^k(x_j)和θ^k(y_j)：

φk(xj)=Sx(θk-1(yj))θk(yj)=Sy[φk(xj)]/||Sy[φk(xj)]||(j=1,2,···,L)---(4)]]>

式中，k是迭代次數；S是指對數據作平滑，平滑方法采用樣條函數法或超級平滑方法；x_j和y_j分別是縱波速度x和橫波速度y的第j個元素；φ^k(x_j)和θ^k(y_j)分別是第k次迭代，縱波速度x和橫波速度y最佳變換φ^k(x)和θ^k(y)的第j個元素；

6）按下式計算迭代誤差：

e2[θk(yj),φk(xj)]=1LΣj=1L[θk(yj)-φk(xj)]2---(5)]]>

式中：k是迭代次數；

7）重復步驟5）和6），直到相鄰兩次迭代誤差滿足下式：

|e²[θ^k(y_j),φ^k(x_j)]-e²[θ^k-1(y_j),φ^k-1(x_j)]|≤δ????????（6）

式中，k是迭代次數；δ是誤差門檻值，取值范圍是0.001到0.01；當相鄰兩次迭代誤差滿足公式（6）時，第k次迭代得到θ(y)稱為橫波速度y的最佳變換函數，得到的φ(x)稱為縱波速度x的最佳變換函數；φ^k(x_j)和θ^k(y_j)分別是第k次迭代，縱波速度x和橫波速度y最佳變換φ^k(x)和θ^k(y)的第j個元素；

8）將步驟7）得到的最佳變換函數θ(y)和φ(x)的函數值，利用線性回歸方法進行回歸，得到θ(y)和φ(x)的函數關系式：

θ(y)=f(φ(x))?????????????????????????????????（7）

式中，f是由線性回歸得方法到的θ(y)和φ(x)的函數關系式；

9）根據函數關系式θ(y)，φ(x)和f，計算得到其他井或地區的縱波速度x或者橫波速度y。

2.根據權利要求1的方法，特點是步驟9）的計算過程是：首先把其他井或地區的縱波速度x或橫波速度y代入函數關系式φ(x)或θ(y)，得到縱波速度x的最佳變換或橫波速度y的最佳變換，再根據公式（7）求取橫波速度y的最佳變換或縱波速度x的最佳變換，最后根據函數關系式θ(y)反求橫波速度y或者根據函數關系式φ(x)反求縱波速度x。