本發(fā)明屬于井中地震勘探領(lǐng)域，具體涉及一種多次波反射點(diǎn)位置計(jì)算方法、系統(tǒng)、存儲(chǔ)介質(zhì)及設(shè)備。

背景技術(shù)：

1、在地表的一階多次波以往經(jīng)常作為干擾波不被利用，在處理中需要被消除掉。近年來(lái)，隨著地震干涉成像技術(shù)研究的發(fā)展，多次波的利用成為一種可能。vsp憑借在井中接收的特點(diǎn)，具有高信噪比、高分辨率的特點(diǎn)，而它在地表的一階多次波又由于射線路徑特殊，照明范圍大的優(yōu)點(diǎn)，在地震勘探領(lǐng)域中有著特殊的應(yīng)用，比如通過(guò)vsp直達(dá)波和地表一階多次波干涉，可以將多次波記錄重構(gòu)成類(lèi)地面地震記錄進(jìn)行成像，進(jìn)而擴(kuò)大vsp的成像范圍，對(duì)井旁特殊構(gòu)造進(jìn)行精細(xì)成像。

2、對(duì)vsp多次波進(jìn)行利用，提前是采集到有利于目標(biāo)體成像的多次波。因此采集前，一般要進(jìn)行多次波采集設(shè)計(jì)論證，而多次波反射點(diǎn)位置的計(jì)算是采集設(shè)計(jì)中的基礎(chǔ)步驟，涉及到后續(xù)覆蓋次數(shù)、成像范圍和拉伸率等參數(shù)的計(jì)算。現(xiàn)有的技術(shù)沒(méi)有專(zhuān)門(mén)針對(duì)多次波采集設(shè)計(jì)的反射點(diǎn)計(jì)算方法，若要計(jì)算多次波反射點(diǎn)，一般通過(guò)通用的射線追蹤來(lái)計(jì)算多次波反射點(diǎn)位置。雖然射線追蹤是一種較高精度的計(jì)算方法，但是它是通過(guò)不斷地嘗試調(diào)整射線的出射方向來(lái)逼近反射點(diǎn)位置，因此計(jì)算非常的大，計(jì)算效率較低(尤其在大型模型或較多炮檢點(diǎn)情況下)，不利于快速計(jì)算多種采集方案，并對(duì)比優(yōu)選方案。

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

1、本發(fā)明的目的在于解決上述現(xiàn)有技術(shù)中存在的難題，提供一種多次波反射點(diǎn)位置計(jì)算方法、系統(tǒng)、存儲(chǔ)介質(zhì)及設(shè)備，既有快速的計(jì)算效率，又有和射線追蹤相當(dāng)?shù)挠?jì)算精度，解決了現(xiàn)有技術(shù)沒(méi)有專(zhuān)門(mén)針對(duì)vsp多次波采集設(shè)計(jì)中反射點(diǎn)位置計(jì)算的高精度、高效率的計(jì)算方法的問(wèn)題。

2、本發(fā)明的目的之一是提供一種多次波反射點(diǎn)位置計(jì)算方法。

3、本發(fā)明的目的之二是提供一種多次波反射點(diǎn)位置計(jì)算系統(tǒng)。

4、本發(fā)明的目的之三是提供一種計(jì)算機(jī)可讀存儲(chǔ)介質(zhì)。

5、本發(fā)明的目的之四是提供一種計(jì)算機(jī)設(shè)備。

6、本發(fā)明是通過(guò)以下技術(shù)方案實(shí)現(xiàn)的：

7、本發(fā)明的第一個(gè)方面，提供一種多次波反射點(diǎn)位置計(jì)算方法，包括以下步驟：

8、步驟1，輸入模型參數(shù)和觀測(cè)系統(tǒng)參數(shù)；

9、步驟2，計(jì)算地表激發(fā)點(diǎn)到目標(biāo)地層深度處的最短時(shí)間；

10、步驟3，計(jì)算地表激發(fā)點(diǎn)到目標(biāo)地層深度處的均方根速度；

11、步驟4，計(jì)算地表激發(fā)點(diǎn)到井中接收點(diǎn)深度處的最短時(shí)間；

12、步驟5，計(jì)算地表激發(fā)點(diǎn)到井中接收點(diǎn)深度處的均方根速度；

13、步驟6，計(jì)算目標(biāo)地層上的多次波反射點(diǎn)位置。

14、本發(fā)明的進(jìn)一步改進(jìn)在于：

15、所述步驟1中輸入模型參數(shù)和觀測(cè)系統(tǒng)參數(shù)，具體包括：

16、輸入模型參數(shù)(hi,vi)，i＝1,2……n，n為模型層數(shù)，hi表示層深度，vi表示層速度；

17、輸入地表激發(fā)點(diǎn)位置xs，井中接收點(diǎn)位置xg；

18、輸入目標(biāo)地層深度hr，地表激發(fā)點(diǎn)深度hs和井中接收點(diǎn)深度hg。

19、本發(fā)明的進(jìn)一步改進(jìn)在于：

20、所述步驟2中計(jì)算地表激發(fā)點(diǎn)到目標(biāo)地層深度處的最短時(shí)間，具體操作包括：

21、步驟201，從地表激發(fā)點(diǎn)深度hs所在層開(kāi)始到目的層深度hr所在層，計(jì)算每層的厚度ssi,其中i＝1,2……n，n為地震波走過(guò)的層數(shù)；

22、步驟202，逐層計(jì)算地震波在每層的最短走時(shí)tsi＝ssi/vi，其中i＝1,2……n，n為地震波走過(guò)的層數(shù)，vi為地震波在第i層的層速度；

23、步驟203，利用公式(1)計(jì)算從地表激發(fā)點(diǎn)到目的層深度處的最短走時(shí)t0sr；

24、

25、本發(fā)明的進(jìn)一步改進(jìn)在于：

26、所述步驟3中計(jì)算地表激發(fā)點(diǎn)到目標(biāo)地層深度處的均方根速度，具體為：

27、利用公式(2)計(jì)算從地表激發(fā)點(diǎn)到反射點(diǎn)深度處的均方根速度vsr：

28、

29、本發(fā)明的進(jìn)一步改進(jìn)在于：

30、所述步驟4中計(jì)算地表激發(fā)點(diǎn)到井中接收點(diǎn)深度處的最短時(shí)間，具體操作包括：

31、步驟401，從地表激發(fā)點(diǎn)深度hs所在層開(kāi)始到井中接收點(diǎn)深度hg所在層，計(jì)算每層的厚度sgi,其中i＝1,2……n，n為地震波走過(guò)的層數(shù)；

32、步驟402，逐層計(jì)算地震波在每層的最短走時(shí)tgi＝sgi/vi，其中i＝1,2……n，n為地震波走過(guò)的層數(shù)，vi為地震波在第i層的層速度；

33、步驟403，利用公式(3)計(jì)算地表激發(fā)點(diǎn)到井中接收點(diǎn)深度處的最短t0sg：

34、

35、本發(fā)明的進(jìn)一步改進(jìn)在于：

36、所述步驟5中計(jì)算地表激發(fā)點(diǎn)到井中接收點(diǎn)深度處的均方根速度，具體為：

37、利用公式(4)計(jì)算從地表激發(fā)點(diǎn)到井中接收點(diǎn)深度處的均方根速度vsg：

38、

39、本發(fā)明的進(jìn)一步改進(jìn)在于：

40、所述步驟6中計(jì)算目標(biāo)地層上的多次波反射點(diǎn)位置，具體為：

41、采用公式(5)計(jì)算計(jì)算多次波反射點(diǎn)位置：

42、

43、本發(fā)明的第二個(gè)方面，提供一種多次波反射點(diǎn)位置計(jì)算系統(tǒng)，包括：

44、輸入模塊，用于輸入模型參數(shù)和觀測(cè)系統(tǒng)參數(shù)；

45、第一計(jì)算模塊，用于計(jì)算地表激發(fā)點(diǎn)到目標(biāo)地層深度處的最短時(shí)間；

46、第二計(jì)算模塊，用于計(jì)算地表激發(fā)點(diǎn)到目標(biāo)地層深度處的均方根速度；

47、第三計(jì)算模塊，用于計(jì)算地表激發(fā)點(diǎn)到井中接收點(diǎn)深度處的最短時(shí)間；

48、第四計(jì)算模塊，用于計(jì)算地表激發(fā)點(diǎn)到井中接收點(diǎn)深度處的均方根速度；

49、第五計(jì)算模塊，用于計(jì)算目標(biāo)地層上的多次波反射點(diǎn)位置。

50、本發(fā)明的第三個(gè)方面，提供一種計(jì)算機(jī)可讀存儲(chǔ)介質(zhì)，所述計(jì)算機(jī)可讀存儲(chǔ)介質(zhì)存儲(chǔ)有計(jì)算機(jī)可執(zhí)行的至少一個(gè)程序，所述至少一個(gè)程序被所述計(jì)算機(jī)執(zhí)行時(shí)使所述計(jì)算機(jī)執(zhí)行如上述所述的多次波反射點(diǎn)位置計(jì)算方法中的步驟。

51、本發(fā)明的第四個(gè)方面，提供一種計(jì)算機(jī)設(shè)備，包括存儲(chǔ)器和處理器，所述存儲(chǔ)器存儲(chǔ)有計(jì)算機(jī)程序，所述計(jì)算機(jī)程序被所述處理器執(zhí)行時(shí)，使得所述處理器執(zhí)行如上述所述的多次波反射點(diǎn)位置計(jì)算方法的步驟。

52、與現(xiàn)有技術(shù)相比，本發(fā)明的有益效果是：

53、通過(guò)實(shí)際測(cè)試表明，本發(fā)明的方法既有快速的計(jì)算效率，又有和射線追蹤相當(dāng)?shù)挠?jì)算精度。

技術(shù)特征：

1.一種多次波反射點(diǎn)位置計(jì)算方法，其特征在于，包括以下步驟：

2.根據(jù)權(quán)利要求1所述的多次波反射點(diǎn)位置計(jì)算方法，其特征在于，所述步驟1中輸入模型參數(shù)和觀測(cè)系統(tǒng)參數(shù)，具體包括：

3.根據(jù)權(quán)利要求2所述的多次波反射點(diǎn)位置計(jì)算方法，其特征在于，所述步驟2中計(jì)算地表激發(fā)點(diǎn)到目標(biāo)地層深度處的最短時(shí)間，具體操作包括：

4.根據(jù)權(quán)利要求3所述的多次波反射點(diǎn)位置計(jì)算方法，其特征在于，所述步驟3中計(jì)算地表激發(fā)點(diǎn)到目標(biāo)地層深度處的均方根速度，具體為：

5.根據(jù)權(quán)利要求4所述的多次波反射點(diǎn)位置計(jì)算方法，其特征在于，所述步驟4中計(jì)算地表激發(fā)點(diǎn)到井中接收點(diǎn)深度處的最短時(shí)間，具體操作包括：

6.根據(jù)權(quán)利要求5所述的多次波反射點(diǎn)位置計(jì)算方法，其特征在于，所述步驟5中計(jì)算地表激發(fā)點(diǎn)到井中接收點(diǎn)深度處的均方根速度，具體為：

7.根據(jù)權(quán)利要求6所述的多次波反射點(diǎn)位置計(jì)算方法，其特征在于，所述步驟6中計(jì)算目標(biāo)地層上的多次波反射點(diǎn)位置，具體為：

8.一種多次波反射點(diǎn)位置計(jì)算系統(tǒng)，其特征在于，包括：

9.一種計(jì)算機(jī)可讀存儲(chǔ)介質(zhì)，其特征在于，所述計(jì)算機(jī)可讀存儲(chǔ)介質(zhì)存儲(chǔ)有計(jì)算機(jī)可執(zhí)行的至少一個(gè)程序，所述至少一個(gè)程序被所述計(jì)算機(jī)執(zhí)行時(shí)使所述計(jì)算機(jī)執(zhí)行如權(quán)利要求1-7任一項(xiàng)所述的多次波反射點(diǎn)位置計(jì)算方法中的步驟。

10.一種計(jì)算機(jī)設(shè)備，其特征在于，包括存儲(chǔ)器和處理器，所述存儲(chǔ)器存儲(chǔ)有計(jì)算機(jī)程序，所述計(jì)算機(jī)程序被所述處理器執(zhí)行時(shí)，使得所述處理器執(zhí)行如權(quán)利要求1-7任一項(xiàng)所述的多次波反射點(diǎn)位置計(jì)算方法的步驟。

技術(shù)總結(jié)
本發(fā)明提供了一種多次波反射點(diǎn)位置計(jì)算方法、系統(tǒng)、存儲(chǔ)介質(zhì)及設(shè)備，屬于井中地震勘探領(lǐng)域。所述計(jì)算方法包括以下步驟：步驟1，輸入模型參數(shù)和觀測(cè)系統(tǒng)參數(shù)；步驟2，計(jì)算地表激發(fā)點(diǎn)到目標(biāo)地層深度處的最短時(shí)間；步驟3，計(jì)算地表激發(fā)點(diǎn)到目標(biāo)地層深度處的均方根速度；步驟4，計(jì)算地表激發(fā)點(diǎn)到井中接收點(diǎn)深度處的最短時(shí)間；步驟5，計(jì)算地表激發(fā)點(diǎn)到井中接收點(diǎn)深度處的均方根速度；步驟6，計(jì)算目標(biāo)地層上的多次波反射點(diǎn)位置。本發(fā)明既有快速的計(jì)算效率，又有和射線追蹤相當(dāng)?shù)挠?jì)算精度。

技術(shù)研發(fā)人員：陳占國(guó),陳浩,張潔,楊心超,楊小慧
受保護(hù)的技術(shù)使用者：中國(guó)石油化工股份有限公司
技術(shù)研發(fā)日：
技術(shù)公布日：2025/5/8