[發明專利]一種地震和測井時深關系自動轉換的方法

有效

• 申請號: 201210250721.9
• 申請日: 2012-07-19
• 公開（公告）號: CN103576199A
• 公開（公告）日: 2014-02-12
• 發明（設計）人: 賀維勝;張立彬;康海霞
• 申請（專利權）人: 中國石油天然氣集團公司;中國石油集團東方地球物理勘探有限責任公司
• 主分類號: G01V1/32
• 分類號: G01V1/32
• 代理公司: 北京市中實友知識產權代理有限責任公司 11013
• 代理人: 劉天語;任潔

1.一種地震和測井時深關系自動轉換的方法，特點是提供以下具體實施步驟：

1）采集探區地震勘探數據，獲得井A位置的地震子波；在井A周圍選取已知測井資料的井B，提取井B位置地震子波；將井B的地震子波作為A井的地震子波，長度為WL；

2）利用井A聲波時差測井數據按照下式計算第i個采樣點處深度d_i所對應的時間t_i，得到初始的時深關系；

d_i＝d₁+(i-1)×Δd????i＝2，…，n

ti=t1+Σj=2iΔd×DTi]]>

其中：t₁是第一個采樣點d₁對應的時間，由下式得到：t₁=d₁×DT₁；d₁是在第一個采樣點的深度；DT₁是第一個采樣點的聲波時差，單位為秒/米；i是采樣點數；d_i是在第i個采樣點對應的深度；DT_i是第i個采樣點的聲波時差；n是測井數據中總的采樣點數；△d是測井數據的采樣間隔；式中深度、時間和聲波時差的單位均為米制標準單位；

3）根據時深關系把井A測井數據轉換到時間域，由下式計算反射系數R_i：

Ri=ρi+1vi+1-ρiviρi+1vi+1+ρivi]]>

其中：ρ和v分別為速度和密度，單位為米制標準單位；i為采樣點數；

4）由下式計算合成地震數據A_s：

A_s＝R_i*W

其中W是步驟1）提取的井B地震子波，R_i是步驟3）得到的反射系數；

5）計算整體時間移動大小；

取一個時窗長度WIN₁,這個時窗長度小于測井數據的實際長度，并且滿足：Wsh≦WIN₁≦TS₂-TS₁+Wsh；

其中：Wsh是時間值；TS₁和TS₂是測井中根據地質分層得到的目標儲層時間范圍；TS₁-Wsh至TS₂+Wsh是合成地震數據的時間范圍；

然后，

①從合成地震數據A_S中取以TS₁-Wsh+WIN₁/2為中心，長度為WIN₁時窗內的數據，記為As₁¹；把As₁¹以△t為采樣間隔，等于地震數據的采樣間隔，在地震數據的范圍內逐漸上、下移動，用下式計算移到位置處的觀測地震數據Ao₁¹和As₁¹的相關系數c：

c=Cor(Ao11,As11)]]>

記錄下相關系數c取得最大值對應的移動時間T_shift¹；

②將合成地震數據的時窗中點下移△t，重復步驟①，得到對應的T_shift²；

③重復步驟②，直到時窗的中點移動到TS₂+Wsh-WIN₁/2，且時窗在有測井數據的范圍以內，得到各個位置處的T_shiftⁱ；

④計算T_shift的概率密度函數；

⑤最大概率密度函數對應的T_shift^max，確定為需要計算的整體時間移動大小；

⑥將合成地震數據整體移動得到新的時深關系；

6）用下式計算合成地震數據A_s2和對應時窗內的合成地震數據A_o2之間的相關系數c_old：

c_old=Cor(A_o2,A_s2)

式中：A_s2是TS₁-Wsh+T_shift^max至TS₂+Wsh+T_shift^max時間內的合成地震數據；

然后，利用下面的步驟局部調整時深關系：

①從合成地震數據A_s2中隨機地選取一個時間點，把該時間點對應的采樣點上移或者下移一個采樣間隔的距離，得到另一合成地震數據A_s2’，計算新相關系數c_new；把迭代次數k的值增加1；

②若新相關系數c_new大于相關系數c_old，則用另一合成地震數據A_s2’替代合成地震數據A_s2，并且讓相關系數c_old等于新相關系數c_new；跳到步驟④；否則，進入下一步；

③隨機產生一個0到1之間的數；如果這個數大于那么用另一合成地震數據A_s2’替代合成地震數據A_s2，并且讓相關系數c_old等于新相關系數c_new；否則，保留移動前的合成地震數據A_s2；其中e為自然對數的底，w為引入數學變量；

④如果迭代次數k是w遞減迭代次數k_L的整數倍，那么用w的遞減因子l和w的積替代w；

⑤如果從新相關系數c_new小于相關系數的最大值c_m并且k小于迭代次數的最大值k_m，那么返回到步驟①；否則，進入下一步；

⑥完成地震和測井時深關系自動轉換。