本發(fā)明涉及儲層預(yù)測研究方法，屬于油田儲層預(yù)測領(lǐng)域，具體涉及一種井震層位逐級匹配的速度場建立方法、介質(zhì)及設(shè)備。

背景技術(shù)：

1、油田開發(fā)后期，剩余油主要富集在注采不完善區(qū)，如斷層邊部、局部構(gòu)造高點等。斷層邊部布井、水平井等挖潛方式已成為常見的挖潛方式，這對地震構(gòu)造精度要求越來越高，需要一種與之適應(yīng)的高精度速度場建立方法。

2、目前開發(fā)地震比較常用的三維速度場建立方式是在地震層位約束下利用井點時深關(guān)系插值建立速度場。該方法存在以下兩點不足：一是用于約束趨勢的地震層位和與之對應(yīng)的井點地質(zhì)分層未完全匹配(僅部分井分層能夠與地震層位匹配上)，這是因為合成地震記錄的質(zhì)量是由合成地震記錄與井旁地震道的波形匹配程度所決定的，而不是依靠井分層與地震層位之間的匹配；二是約束趨勢的地震層位較少，一般為油層組級的或者是個別在地震剖面由明顯反射特征的地質(zhì)界面，該方法無法得到沉積單元級地震層位，也就無法做到沉積單元級層位匹配與約束。

3、因此，如何提高速度場的精度成為了本領(lǐng)域亟需解決的問題。

技術(shù)實現(xiàn)思路

1、針對上述問題，本發(fā)明提供了一種井震層位逐級匹配的速度場建立方法、介質(zhì)及設(shè)備，該方法、介質(zhì)及設(shè)備基于多井一致性分析，可有效提高地震構(gòu)造解釋精度，保障了剩余油挖潛效果。

2、為實現(xiàn)上述技術(shù)目的，本發(fā)明采用的技術(shù)方案如下：

3、根據(jù)本發(fā)明的一方面，提供一種井震層位逐級匹配的速度場建立方法，該方法包括下列步驟：

4、步驟1：準備目標井位的測井數(shù)據(jù)、地震數(shù)據(jù)和目的層頂、底面地震層位數(shù)據(jù)h頂和h底；

5、步驟2：通過制作地震合成記錄得到每口井的初始時深曲線td1；

6、步驟3：在目的層頂、底面地震層位數(shù)據(jù)h頂和h底趨勢約束下，利用井點初始時深曲線td1插值建立第一速度場v1；

7、步驟4：利用第一速度場v1將目的層頂、底分層top深頂、top深底由深度域轉(zhuǎn)換到時間域top時頂、top時底；

8、步驟5：在目的層頂、底地震層位趨勢約束下，通過時間域的頂、底分層top時頂、top時底插值得到新的頂、底面地震層位h頂1和h底1，對h頂1和h底1分別進行一致性分析和異常值處理得到最終的頂、底面地震層位h頂和h底；

9、步驟6：利用深度域頂、底面分層top深頂、top深底與步驟5得到的最終的頂、底面地震層位h頂和h底對應(yīng)匹配得到時深對t-d，應(yīng)用時深對t-d對初始時深曲線td1進行校正得到每口井校正后的時深曲線td2；

10、步驟7：在最終的頂、底面地震層位h頂和h底趨勢約束下，利用校正后的時深曲線td2插值建立第二速度場v2；

11、步驟8：利用第二速度場v2將深度域目的層內(nèi)沉積單元級分層top深1、top深2......top深n轉(zhuǎn)換到時間域top時1、top時2......top時n；

12、步驟9：在目的層最終的頂、底面地震層位h頂和h底趨勢約束下，通過時間域的目的層內(nèi)沉積單元級分層top時1、top時2......top時n插值得到初始目的層內(nèi)各沉積單元頂面地震層位h1、h2......hn，對h1、h2......hn分別進行一致性分析和異常值處理得到最終的目的層內(nèi)各沉積單元頂面地震層位hd1、hd2......hdn；

13、步驟10：利用深度域目的層內(nèi)沉積單元級分層top深1、top深2......top深n、目的層底面分層與步驟9得到的最終的目的層內(nèi)各沉積單元頂面地震層位hd1、hd2......hdn、h底對應(yīng)匹配得到每口井的最終時深曲線td3；

14、步驟11：在最終的目的層內(nèi)各沉積單元頂面地震層位hd1、hd2......hdn、h底趨勢約束下，利用最終時深曲線td3插值建立第三速度場v3。

15、根據(jù)本發(fā)明的一個實施例，步驟1中的數(shù)據(jù)通過常規(guī)測井解釋和地震構(gòu)造解釋得到。

16、根據(jù)本發(fā)明的一個實施例，步驟2包括：由聲波和密度測井曲線計算得到反射系數(shù)，將反射系數(shù)與提取的地震子波進行褶積得到初始合成地震記錄，初始合成地震記錄與井旁地震道匹配調(diào)整，得到最終合成地震記錄，最終合成地震記錄上每一個深度點都對應(yīng)一個時間值，從而得到每口井的初始時深曲線td1。

17、根據(jù)本發(fā)明的一個實施例，步驟3包括：在時間域按照等比例將目的層分成n個小層，每個小層的厚度為d，則：

18、d＝(h底-h頂)/n

19、其中，h底和h頂是地震構(gòu)造解釋得到的目的層頂、底面地震層位，n表示目的層所分的小層個數(shù)；

20、在h頂和h底之間得到n-1個地震層位，通過計算得出每個小層的地震層位hi：

21、hi＝h頂+i*d

22、其中，i表示是第幾個小層，1≤i≤n-1；

23、基于初始時深曲線td1將時間域每口井與地震層位hi相交點的時深對提取出來進行井間插值得到小層hi的平面速度vi，當n足夠大時，hi的平面速度vi就是hi層的速度，從而得到目的層位的第一速度場v1。

24、根據(jù)本發(fā)明的一個實施例，步驟4包括：根據(jù)h＝vt/2得到t＝2h/v，其中t是地震雙層旅行時間，h是深度，v是速度，通過第一速度場v1利用t＝2h/v公式將目的層頂、底分層top深頂、top深底由深度域轉(zhuǎn)換到時間域top時頂、top時底。

25、根據(jù)本發(fā)明的一個實施例，步驟5包括：以間隔預(yù)定距離的兩井為界，通過插值得到新的頂、底面地震層位h頂1和h底1。

26、根據(jù)本發(fā)明的一個實施例，其中當兩井間距大于預(yù)定距離時，井間分層插值趨勢以地震層位趨勢為準，當兩井間距小于預(yù)定距離時，井間趨勢以兩井時間域分層插值趨勢為準。

27、根據(jù)本發(fā)明的一個實施例，步驟5中對h頂1和h底1分別進行一致性分析和異常值處理包括：當個別井時間域的分層與其附近的井差異較大時在插出的新的地震層位上出現(xiàn)異常凸起或凹陷在該井位置與地震層位的整體趨勢不一致的情況下，通過對該井位置進行異常值剔除或平滑處理得到最終的頂、底面地震層位h頂和h底。

28、根據(jù)本發(fā)明的一個實施例，步驟6包括：用深度域井點處油層組頂面分層與最終的頂面地震層位h頂進行匹配，用油層組底面分層與最終的底面地震層位h底進行匹配得到頂、底界面時深對t-d，應(yīng)用時深對t-d對初始時深曲線td1進行校正得到每口井校正后的時深曲線td2。

29、根據(jù)本發(fā)明的一個實施例，步驟7包括：在時間域按照等比例將目的層位分成n’個小層，每個小層的厚度為d’，則：

30、d’＝(h底-h頂)/n’

31、其中，h底和h頂是步驟5得到的最終的頂、底面地震層位，n’表示所分的小層個數(shù)；

32、在h頂和h底之間得到n’-1個地震層位，通過計算得出每個小層的地震層位hi：

33、hi＝h頂+i’*d’

34、其中，i’表示是第幾個小層，1≤i’≤n’-1；

35、基于校正后的時深曲線td2將時間域每口井與地震層位hi相交點的時深對提取出來進行井間插值得到小層hi的平面速度vi’，當n’足夠大時，hi的平面速度vi’就是hi’層的速度，從而得到目的層位的第二速度場v2。

36、根據(jù)本發(fā)明的一個實施例，步驟8包括：根據(jù)h＝vt/2得到t＝2h/v，其中t是地震雙層旅行時間，h是深度，v是速度，通過第二速度場v2利用t＝2h/v公式將各個沉積單元分層頂、底分層由深度域轉(zhuǎn)換到時間域。

37、根據(jù)本發(fā)明的一個實施例，步驟9包括：以間隔第二預(yù)定距離的兩井為界，通過對時間域的目的層內(nèi)沉積單元級分層top時1、top時2......top時n插值得到初始目的層內(nèi)各沉積單元頂面地震層位h1、h2......hn。

38、根據(jù)本發(fā)明的一個實施例，其中當兩井間距大于第二預(yù)定距離時，井間分層插值趨勢以地震層位趨勢為準，當兩井間距小于第二預(yù)定距離時，井間趨勢以兩井時間域分層插值趨勢為準。

39、根據(jù)本發(fā)明的一個實施例，步驟9中對h1、h2......hn分別進行一致性分析和異常值處理包括：當個別井時間域的分層與其附近的井差異較大時在插出的新的地震層位上出現(xiàn)異常凸起或凹陷在該井位置與地震層位的整體趨勢不一致的情況下，通過對該井位置進行異常值剔除或平滑處理得到最終的目的層內(nèi)各沉積單元頂面地震層位hd1、hd2......hdn。

40、根據(jù)本發(fā)明的一個實施例，步驟10包括：用深度域井點處目的層內(nèi)沉積單元級分層top深1與沉積單元頂面地震層位hd1匹配、分層top深2與地震層位hd2匹配......分層top深n與地震層位hdn匹配、目的層底面分層與h底匹配，得到n+1組時深對，采用這n+1組時深對來對校正后的時深曲線td2進行校正得到每口井的最終時深曲線td3。

41、根據(jù)本發(fā)明的一個實施例，步驟11包括：將每個沉積單元級層段都按照步驟3的方法應(yīng)用最終時深曲線td3進行各個層的速度求取，當所有沉積單元級層段的層的速度都求取完，得到最終要建立的第三速度場v3。

42、根據(jù)本發(fā)明的一個實施例，在步驟10中油層組頂與第一個沉積單元的頂是同一分層。

43、根據(jù)本發(fā)明的一個實施例，步驟3和步驟5中根據(jù)目的層厚度和所要建立的速度的精度要求確定小層個數(shù)。

44、根據(jù)本發(fā)明的另一方面，提供一種計算機可讀存儲介質(zhì)，所述計算機可讀存儲介質(zhì)存儲有計算機程序，其特征在于，所述計算機程序被處理器執(zhí)行時執(zhí)行如上所述方法的步驟。

45、根據(jù)本發(fā)明的又一方面，提供一種計算機設(shè)備，包括：

46、至少一個處理器；以及

47、存儲器，所述存儲器存儲有可在所述處理器上運行的計算機程序，其特征在于，所述處理器執(zhí)行所述程序時執(zhí)行如上所述方法的步驟。

48、由于采用以上技術(shù)方案，本發(fā)明提供的方法、介質(zhì)及設(shè)備與現(xiàn)有技術(shù)相比具有以下有益效果中的至少一項：

49、(1)通過對井點地質(zhì)分層與地震層位進行逐級匹配分析與一致性校正得到井震層位一致性較好的時深曲線，從而得到研究區(qū)目的層的高精度速度場；

50、(2)在速度場求取過程中得到的沉積單元級的地震層位可用于地震屬性分析、地震反演的格架模型建立等；

51、(3)最終得到的高精度速度場實現(xiàn)地震數(shù)據(jù)體、地震層位等的時深轉(zhuǎn)換，特別是地震反演數(shù)據(jù)體的時深轉(zhuǎn)換，可用于指導(dǎo)措施方案的編制。