王維紅,柯 璇,郭雪豹,張瑩瑩,于 鵬

(1.東北石油大學(xué)地球科學(xué)學(xué)院,黑龍江大慶163318;2.中國(guó)石油化工股份有限公司東北油氣分公司,吉林長(zhǎng)春130062)

基于Qt的Walkaway VSP逆時(shí)偏移軟件開發(fā)研究

王維紅1,柯 璇1,郭雪豹1,張瑩瑩1,于 鵬2

(1.東北石油大學(xué)地球科學(xué)學(xué)院,黑龍江大慶163318;2.中國(guó)石油化工股份有限公司東北油氣分公司,吉林長(zhǎng)春130062)

VSP觀測(cè)系統(tǒng)的檢波器位于井筒內(nèi),地震資料的信噪比高,波場(chǎng)信息豐富,利用基于雙程波的疊前逆時(shí)偏移方法可實(shí)現(xiàn)VSP資料的高精度成像,也為井筒附近復(fù)雜介質(zhì)和目的層流體性質(zhì)的研究提供基礎(chǔ)資料。在Walkaway VSP資料逆時(shí)偏移算法研究的基礎(chǔ)上,在Linux操作系統(tǒng)下,以Qt環(huán)境為開發(fā)平臺(tái),設(shè)計(jì)了一套較為完整的VSP資料逆時(shí)偏移處理流程,初步開發(fā)了Walkaway VSP資料逆時(shí)偏移計(jì)算軟件。該軟件具有CPU和GPU兩套核心算法模塊,可適應(yīng)不同的硬件系統(tǒng),理論模型數(shù)據(jù)測(cè)試結(jié)果表明,該軟件平臺(tái)具有操作性強(qiáng)、移植性好和成像精度高的特點(diǎn)。

VSP;逆時(shí)偏移;Qt;軟件界面

逆時(shí)偏移(RTM)理論自Whitmore等1983年提出之后[1-3],許多地球物理工作者對(duì)其進(jìn)行了深入的研究,但是由于其計(jì)算量巨大,而且在計(jì)算過程中產(chǎn)生大量中間數(shù)據(jù),對(duì)計(jì)算機(jī)的計(jì)算能力和存儲(chǔ)性能都有很高的要求,所以在問世之初未能得到廣泛應(yīng)用。近年來,隨著計(jì)算機(jī)硬件技術(shù)的發(fā)展,逆時(shí)偏移算法也得到很好的發(fā)展和應(yīng)用。逆時(shí)偏移本身基于雙程波方程,同時(shí)由于直接對(duì)波動(dòng)方程進(jìn)行離散,因此能夠準(zhǔn)確模擬地下復(fù)雜構(gòu)造的波場(chǎng)特征,從而大幅度提高地震資料的成像精度。在波場(chǎng)數(shù)值模擬過程中,常用的方法包括有限差分、有限元和偽譜法[4-6],其中有限差分法由于其實(shí)現(xiàn)相對(duì)簡(jiǎn)單而得到普遍應(yīng)用。在邊界條件方面,主要有隨機(jī)邊界和吸收邊界條件等。隨機(jī)邊界的引入主要是為了解決逆時(shí)偏移過程中的波場(chǎng)數(shù)據(jù)存儲(chǔ)問題,但實(shí)際上也會(huì)將隨機(jī)噪聲引入到成像過程中[7];吸收邊界中最常采用的是完全匹配層吸收邊界條件,其吸收效果好,理論上能夠吸收來自任何方向的入射波,但是計(jì)算效率較低,同時(shí)在偏移計(jì)算過程中產(chǎn)生的臨時(shí)波場(chǎng)數(shù)據(jù)量較大,一般的硬件系統(tǒng)限制該方法的有效應(yīng)用。

隨著石油勘探和開發(fā)目標(biāo)的復(fù)雜化,加之從勘探走向開發(fā)是當(dāng)今地震技術(shù)的重要發(fā)展趨勢(shì)之一,井中地震資料的高精度成像處理成為地震數(shù)據(jù)處理的必然要求,VSP地震資料以其精度高和直接測(cè)量目的層的特點(diǎn),得到開發(fā)人員和油藏地球物理工作者的認(rèn)可。針對(duì)VSP地震資料的成像方法,許多學(xué)者進(jìn)行了深入研究。Wyatt[8]提出了用時(shí)間域正演計(jì)算合成VSP地震記錄的方法,為后續(xù)的VSP地震資料成像研究提供數(shù)據(jù)基礎(chǔ);Lee等[9]提出了針對(duì)VSP地震數(shù)據(jù)的處理方法,并指出利用VSP地震數(shù)據(jù)可以有效進(jìn)行子波整形、估算巖石波阻抗;Dillon[10]提出了對(duì)VSP數(shù)據(jù)進(jìn)行Kirchhoff積分偏移計(jì)算,使得VSP地震資料的偏移歸位更為準(zhǔn)確,為資料的應(yīng)用奠定基礎(chǔ)。

前已述及,VSP資料對(duì)井周地層的精細(xì)刻畫及深部構(gòu)造高精度成像具有重要的意義，為進(jìn)一步提高井周微幅構(gòu)造的成像效果,需要研究基于VSP資料的逆時(shí)偏移方法。劉清林等[11]和朱金明等[12-13]提出了VSP資料的逆時(shí)偏移處理方法,但是包括上述作者在內(nèi)的諸多學(xué)者對(duì)VSP資料逆時(shí)偏移的研究都僅限于理論與方法的研究。目前地震資料的逆時(shí)偏移軟件多為國(guó)外產(chǎn)品。業(yè)界針對(duì)VSP地震資料處理的商業(yè)軟件僅局限于速度分析等方面,尚無VSP資料的逆時(shí)偏移算法平臺(tái)。本文在軟件開發(fā)技術(shù)和VSP逆時(shí)偏移方法研究的基礎(chǔ)上,在Linux操作系統(tǒng)環(huán)境下,應(yīng)用Qt開發(fā)了圖形用戶界面,形成了一套基于Qt的Walkaway VSP逆時(shí)偏移處理軟件,從而填補(bǔ)了VSP地震資料逆時(shí)偏移處理軟件的空白,為VSP資料高精度成像方法研究和實(shí)際應(yīng)用奠定了基礎(chǔ)。

1 VSP逆時(shí)偏移原理

不同于傳統(tǒng)的地面地震勘探,VSP觀測(cè)系統(tǒng)是在地面激發(fā),在井筒內(nèi)的一定深度范圍布置檢波器,接收來自地下界面多次覆蓋的地震數(shù)據(jù),這在一定程度上降低了地面噪聲的干擾,井中檢波器記錄的地震反射波無需傳播至地面,旅行路徑更短。因此,VSP地震資料具有更高的信噪比和更小的能量衰減。VSP數(shù)據(jù)偏移是資料有效應(yīng)用的基礎(chǔ),目前常規(guī)VSP偏移處理技術(shù),如VSP-CDP轉(zhuǎn)換、射線法等,難以獲得理想的成像效果,逆時(shí)偏移方法基于雙程波方程,結(jié)合VSP數(shù)據(jù)資料特點(diǎn),可對(duì)井旁復(fù)雜構(gòu)造精確成像。

1.1 VSP逆時(shí)偏移的基本原理

根據(jù)VSP資料和逆時(shí)偏移算法的特點(diǎn),在實(shí)際計(jì)算時(shí)首先進(jìn)行震源波場(chǎng)的正演模擬,再在接收點(diǎn)位置載入檢波點(diǎn)波場(chǎng),并進(jìn)行反傳播模擬,然后將互相關(guān)成像條件應(yīng)用于上述兩種波場(chǎng),進(jìn)行同時(shí)刻成像,最后疊加生成偏移數(shù)據(jù)體。對(duì)Walkaway VSP地震資料而言,正、反傳播波場(chǎng)的計(jì)算公式表述如下：

(1)

(2)

(3)

應(yīng)用有限差分法對(duì)波場(chǎng)進(jìn)行數(shù)值模擬,其時(shí)間二階、空間2M階精度的有限差分格式如下：

(4)

式中：i,j和m分別為波場(chǎng)在x,z和t方向上的離散序號(hào);Δx,Δz和Δt分別為x,z和t方向的采樣間隔;Cl為有限差分系數(shù)。本文應(yīng)用12階的高階差分格式進(jìn)行正演模擬和逆時(shí)偏移計(jì)算,同時(shí)為進(jìn)一步提高波動(dòng)方程數(shù)值解的精度,應(yīng)用了優(yōu)化差分系數(shù)的方法[14]。

1.2 隨機(jī)邊界條件

隨機(jī)邊界通過在速度模型邊界設(shè)置隨機(jī)速度層,從而降低邊界反射的相干性,減小其對(duì)成像結(jié)果的影響。由于該方法僅在速度模型的外圍區(qū)域設(shè)置隨機(jī)速度層,并未對(duì)波動(dòng)方程做出改變,因此不會(huì)造成波場(chǎng)信息的任何損失,對(duì)于震源波場(chǎng)正傳來說,可以通過保存正傳波場(chǎng)最后幾個(gè)時(shí)刻(一般為兩個(gè))的波場(chǎng)信息,再進(jìn)行反傳達(dá)到節(jié)省存儲(chǔ)空間的目的。其在邊界部位的隨機(jī)函數(shù)構(gòu)造如下：

(5)

式中：v(x,z)為構(gòu)造邊界的隨機(jī)速度;v1(x,z)為構(gòu)造邊界的原始速度;r為隨機(jī)數(shù);d為內(nèi)邊界與速度點(diǎn)之間的距離,其保證了速度的隨機(jī)性隨距離的增加而增加[15]。

逆時(shí)計(jì)算時(shí),首先將震源波場(chǎng)正傳至最大時(shí)刻,并存儲(chǔ)最后兩個(gè)時(shí)刻的波場(chǎng)信息,再利用這兩個(gè)時(shí)刻的波場(chǎng)信息逆推之前各個(gè)時(shí)刻的波場(chǎng),這一過程與檢波點(diǎn)波場(chǎng)反傳同時(shí)進(jìn)行,即可同時(shí)得到每一時(shí)刻的震源波場(chǎng)和檢波點(diǎn)波場(chǎng),應(yīng)用一定的成像條件(互相關(guān)成像條件或者震源歸一化成像條件)進(jìn)行成像,直至零時(shí)刻,即可得到逆時(shí)偏移結(jié)果。隨機(jī)邊界方法可以極大地節(jié)省存儲(chǔ)空間,但計(jì)算量需增加一倍,即隨機(jī)邊界是一種以計(jì)算換存儲(chǔ)的計(jì)算方案,其計(jì)算效率的問題可通過后續(xù)給出的GPU加速技術(shù)得到有效解決。

1.3 成像條件

逆時(shí)偏移中常用的成像條件有激發(fā)時(shí)刻成像條件、反褶積成像條件和互相關(guān)類成像條件等,由于互相關(guān)類成像條件實(shí)施簡(jiǎn)易且計(jì)算穩(wěn)定而得到普遍采用,互相關(guān)成像條件可表示為：

(6)

式中：s(x,z,t)代表t時(shí)刻的震源波場(chǎng);r(x,z,t)代表t時(shí)刻的檢波點(diǎn)波場(chǎng);s(x,z,t)與r(x,z,t)的乘積表示在t時(shí)刻對(duì)波場(chǎng)進(jìn)行一次成像運(yùn)算;求和計(jì)算代表最終所成的像是各時(shí)間步長(zhǎng)成像的疊加。由此可知,互相關(guān)成像條件充分利用了雙程波場(chǎng)進(jìn)行成像,是地面和井中地震資料逆時(shí)偏移中廣泛采用的成像條件。

針對(duì)VSP資料的特點(diǎn),為保證目的層的成像效果,提高目的層的成像精度,本文采用的成像條件為震源歸一化互相關(guān)成像條件,其計(jì)算公式如下：

(7)

式中：ε為穩(wěn)定系數(shù),保證成像過程計(jì)算穩(wěn)定。

與(6)式所示的傳統(tǒng)互相關(guān)成像條件相比,震源歸一化互相關(guān)成像在成像數(shù)據(jù)的振幅上保持相對(duì)均衡,并且可在一定程度上補(bǔ)償深層的能量,提高深部的成像效果。

1.4 Laplacian去噪

由逆時(shí)偏移相關(guān)成像條件的原理可知,在震源與檢波點(diǎn)之間的射線路徑上均存在滿足相關(guān)條件的成像點(diǎn),在對(duì)有效構(gòu)造進(jìn)行成像的同時(shí),也會(huì)產(chǎn)生低頻噪聲,可采用Laplacian噪聲壓制方法壓制該噪聲,從而提高疊加數(shù)據(jù)體的信噪比[16-17]。Laplacian算子可表示為：

(8)

式中：kx,ky,kz為成像域波數(shù)在空間方向的分量;kI為成像域的波數(shù)矢量;ω為角頻率;v為介質(zhì)速度;θ為入射角。此處對(duì)Laplacian低頻噪聲壓制方法的實(shí)現(xiàn)步驟不做重點(diǎn)闡述。

2 GPU/CPU并行加速

與地面地震資料的逆時(shí)偏移方法相同,為降低Walkaway VSP逆時(shí)偏移算法的計(jì)算成本,采用基于CUDA平臺(tái)的GPU并行加速算法提高計(jì)算效率。因GPU擁有更多的計(jì)算核心,因此,多線程調(diào)用可實(shí)現(xiàn)大規(guī)模并行計(jì)算,提高Walkaway VSP逆時(shí)偏移方法的計(jì)算效率。

GPU和CPU協(xié)同并行計(jì)算中,將需要GPU加速計(jì)算的部分打包成函數(shù)形式,稱之為kernel函數(shù),可在GPU上直接運(yùn)行加速算法。在調(diào)用kernel函數(shù)前,需劃分計(jì)算量,并將其分配給不同的線程(Thread)計(jì)算,主要以Thread為基本單位。但是,CPU和GPU各具優(yōu)勢(shì),因此,在編程過程中應(yīng)根據(jù)其特性予以考慮,以充分利用兩者優(yōu)點(diǎn)達(dá)到節(jié)省計(jì)算資源和提高計(jì)算速度的目的。

當(dāng)函數(shù)在GPU端計(jì)算時(shí),線程(Thread)為基本單位,若干線程組成塊(Block),若干塊組成柵格(Grid),如圖1所示。Block的大小具有一定限制,主要由顯卡的型號(hào)決定,CUDA有對(duì)應(yīng)變量對(duì)Block和Grid進(jìn)行設(shè)置,并且每個(gè)線程都有編號(hào)與之對(duì)應(yīng)。GPU上具有多種存儲(chǔ)器,其結(jié)構(gòu)如圖2所示。圖中僅顯示常用的存儲(chǔ)器,在調(diào)用kernel函數(shù)時(shí),需先將數(shù)據(jù)傳輸至全局存儲(chǔ)器內(nèi),每個(gè)Block對(duì)應(yīng)一個(gè)共享存儲(chǔ)器,共享存儲(chǔ)器的大小固定,每個(gè)Block內(nèi)的線程均可從共享存儲(chǔ)器中讀取數(shù)據(jù),Block之間的共享存儲(chǔ)器互不通信,其數(shù)據(jù)交互通過全局存儲(chǔ)器進(jìn)行,每個(gè)線程對(duì)應(yīng)一個(gè)寄存器,并且只能從其對(duì)應(yīng)寄存器內(nèi)存取數(shù)據(jù),但寄存器本身很小,在線程進(jìn)行運(yùn)算時(shí),從共享存儲(chǔ)器和寄存器進(jìn)行數(shù)據(jù)讀取較全局存儲(chǔ)器更快,但由于這兩種存儲(chǔ)器本身存在限制,因此需充分考慮實(shí)際計(jì)算的需求。

圖1 GPU編程模式

圖2 GPU存儲(chǔ)器結(jié)構(gòu)

波場(chǎng)數(shù)值模擬是逆時(shí)偏移的核心,也是最耗時(shí)的部分,其本身包含大量的并行運(yùn)算。對(duì)震源波場(chǎng)正傳和檢波點(diǎn)波場(chǎng)反傳進(jìn)行加速計(jì)算,可以明顯提高逆時(shí)偏移的計(jì)算效率。實(shí)際編程時(shí),首先利用CPU在串行計(jì)算方面的優(yōu)勢(shì),從文件中讀取相應(yīng)的數(shù)據(jù),包括震源子波、速度模型等參數(shù),將其傳輸至全局存儲(chǔ)器中,然后對(duì)計(jì)算任務(wù)進(jìn)行優(yōu)化分配,以充分利用計(jì)算資源,調(diào)用kernel函數(shù)進(jìn)行波場(chǎng)正傳的計(jì)算,并采用相應(yīng)的存儲(chǔ)策略。之后進(jìn)行檢波點(diǎn)處波場(chǎng)反傳,在反傳過程中,利用成像條件進(jìn)行成像。最后,將所需要的數(shù)據(jù)從全局存儲(chǔ)器中傳回主機(jī)端進(jìn)行保存。利用GPU并行計(jì)算,在很大程度上降低了逆時(shí)偏移算法的計(jì)算成本。

3 軟件設(shè)計(jì)

3.1 Qt環(huán)境

Qt是基于C++的跨平臺(tái)圖形界面程序開發(fā)軟件,具有跨系統(tǒng)平臺(tái)移植性,可提供“一次編寫,隨處編譯”的開發(fā)框架[18-19]。Qt通過信號(hào)和槽機(jī)制完成任意Qt對(duì)象之間的通信操作,整個(gè)過程由開發(fā)者自行定義,信號(hào)和槽的匹配關(guān)系不受數(shù)量限制,可提高對(duì)象處理的靈活性與多樣性。創(chuàng)建軟件界面主要有兩種方式：一是通過Qt平臺(tái)設(shè)計(jì)工具繪制圖形;二是編寫程序代碼設(shè)計(jì)圖形界面。本文采用后者,這樣可為后續(xù)界面組件的編輯和調(diào)試提供方便,同時(shí)形成的軟件系統(tǒng)也具有良好的可移植性。

3.2 架構(gòu)設(shè)計(jì)

軟件的架構(gòu)設(shè)計(jì)與軟件的執(zhí)行效率密切相關(guān),合理的架構(gòu)設(shè)計(jì)可提高軟件的執(zhí)行效率,大幅縮短地震資料的處理周期。軟件架構(gòu)設(shè)計(jì)思路如圖3所示。對(duì)逆時(shí)偏移算法模塊進(jìn)行預(yù)編譯,生成可執(zhí)行程序,軟件界面向可執(zhí)行程序發(fā)送參數(shù)和命令,反過來可執(zhí)行程序給軟件界面提供功能支持。在界面代碼執(zhí)行過程中,僅進(jìn)行參數(shù)傳遞等操作,也就是說直接調(diào)用前期編譯完成的算法可執(zhí)行程序。這樣做的優(yōu)點(diǎn)主要體現(xiàn)在：形成的軟件平臺(tái)避免了界面代碼與算法代碼的混合編譯,使編譯更容易,計(jì)算效率更高,代碼算法也更易于實(shí)現(xiàn)。另外需要指出的是,按照上述方法設(shè)計(jì)軟件具有很好的可移植性和擴(kuò)展性,只要提供可執(zhí)行軟件和必要的參數(shù)接口,就可以完成相關(guān)計(jì)算程序的界面化移植。

圖3 軟件執(zhí)行路線

3.3 流程設(shè)計(jì)

軟件平臺(tái)的操作流程和技術(shù)流程對(duì)于軟件的執(zhí)行和應(yīng)用效果至關(guān)重要,作業(yè)流程的設(shè)計(jì)將直接影響軟件的執(zhí)行效率。VSP逆時(shí)偏移軟件平臺(tái)的基本流程為工區(qū)建立、數(shù)據(jù)載入、參數(shù)設(shè)置、作業(yè)發(fā)送等4個(gè)關(guān)鍵步驟,根據(jù)軟件系統(tǒng)不同子模塊的結(jié)構(gòu)和VSP逆時(shí)偏移的基本原理,設(shè)計(jì)了基于GPU并行計(jì)算的技術(shù)流程圖(圖4)。

圖4 VSP地震數(shù)據(jù)逆時(shí)偏移處理技術(shù)流程

4 軟件測(cè)試

在VSP逆時(shí)偏移算法的基礎(chǔ)上,基于Qt環(huán)境開發(fā)了垂直地震剖面數(shù)據(jù)逆時(shí)偏移的軟件平臺(tái),為測(cè)試算法的準(zhǔn)確性和軟件平臺(tái)的實(shí)用性,設(shè)計(jì)如圖5 所示的復(fù)雜斷陷盆地模型。該模型包括了高陡構(gòu)造、小斷層和尖滅等復(fù)雜地質(zhì)構(gòu)造,模型在x方向長(zhǎng)度為2700m,z方向深度為2800m。在利用該模型進(jìn)行VSP逆時(shí)偏移軟件測(cè)試前,首先進(jìn)行VSP數(shù)據(jù)的正演模擬,在模型沿x方向0和2700m處共布置兩口VSP井,正演觀測(cè)系統(tǒng)參數(shù)如下：共270炮激發(fā),炮間距10m;井中共280級(jí)檢波器接收,檢波器間距10m;空間網(wǎng)格采用邊長(zhǎng)10m的正方形;時(shí)間采樣間隔1ms;震源子波為主頻25Hz的Ricker子波;記錄長(zhǎng)度2s。圖6為正演模擬的單炮記錄。其中圖6a為模型x=0處VSP井接收到的炮點(diǎn)坐標(biāo)為0的單炮記錄;圖6b 為模型x=2700m處VSP井接收到的炮點(diǎn)坐標(biāo)為1000m的單炮記錄。分析可知,模擬資料的上、下行波波場(chǎng)特征清晰。由于正演模擬所用的模型較為復(fù)雜,導(dǎo)致炮記錄上的時(shí)距曲線特征較為復(fù)雜,該模擬數(shù)據(jù)將用于測(cè)試VSP逆時(shí)偏移算法軟件的成像效果。

將正演模擬得到的270炮記錄作為數(shù)據(jù)輸入,基于圖5所示的速度模型,完成了逆時(shí)偏移試算,經(jīng)過疊加和低頻噪聲壓制后,得到的成像剖面如圖7a 所示。為進(jìn)行成像效果的對(duì)比分析,圖7b給出了基于同樣速度模型的地面地震數(shù)據(jù)的逆時(shí)偏移成像剖面,需要說明的是該數(shù)據(jù)同樣是采取270炮激發(fā),270個(gè)檢波器布置于地面0～2700m處,檢波點(diǎn)間距10m,除檢波點(diǎn)位置不同外,其它的所有參數(shù)都和VSP資料的參數(shù)一致。對(duì)比VSP資料與地面地震資料的成像剖面可知,VSP逆時(shí)偏移算法對(duì)模型底部成像更清晰,主要原因是檢波器在井下,可實(shí)現(xiàn)對(duì)目的層的直接照明和成像,從而使得深層介質(zhì)及陡傾角地層的成像效果比地面地震更清晰。

圖5 復(fù)雜斷陷盆地模型

在成像精度測(cè)試的基礎(chǔ)上,對(duì)軟件平臺(tái)所包括的CPU和GPU并行方式的Walkaway VSP逆時(shí)偏移模塊進(jìn)行了計(jì)算效率對(duì)比。利用前述模擬的270炮VSP地震記錄進(jìn)行了逆時(shí)偏移計(jì)算時(shí)間測(cè)試，測(cè)試結(jié)果如表1所示?；贑UDA平臺(tái)GPU的Walkaway VSP逆時(shí)偏移方法可獲得較高的計(jì)算效率,相比傳統(tǒng)的CPU計(jì)算方式,可加速10倍左右。

圖6 正演模擬的單炮記錄

圖7 VSP資料(a)與地面地震資料(b)逆時(shí)偏移成像效果對(duì)比

表1 CPU和GPU不同計(jì)算方式Walkaway VSP逆時(shí)偏移計(jì)算時(shí)間對(duì)比

計(jì)算方式單炮計(jì)算時(shí)間/s270炮計(jì)算時(shí)間/sCPU12634020GPU123240

5 結(jié)束語

逆時(shí)偏移是當(dāng)前地震資料高精度成像的重要算法,Walkaway VSP地震資料的逆時(shí)偏移可實(shí)現(xiàn)復(fù)雜勘探和開發(fā)研究區(qū)構(gòu)造的高精度成像?；赒t開發(fā)工具集成Walkaway VSP逆時(shí)偏移算法,進(jìn)而形成了軟件系統(tǒng)平臺(tái),為實(shí)際Walkaway VSP資料的高精度成像計(jì)算提供了方便。結(jié)合CPU和GPU兩種計(jì)算架構(gòu),實(shí)現(xiàn)了Walkaway VSP的逆時(shí)偏移計(jì)算,擴(kuò)展了軟件的應(yīng)用范圍,GPU并行加速算法,可在傳統(tǒng)CPU計(jì)算的基礎(chǔ)上,大幅提高計(jì)算效率,縮短成像計(jì)算的周期。復(fù)雜斷陷盆地模型的成像試算結(jié)果表明,所形成的軟件平臺(tái)具有計(jì)算精度高、計(jì)算效率高和實(shí)用性強(qiáng)的優(yōu)點(diǎn)。該軟件開發(fā)的研究實(shí)踐對(duì)同類軟件的設(shè)計(jì)和開發(fā)具有一定的借鑒意義。

[1] Whitmore N D.Iterative depth migration by backward time propagation[J].Expanded Abstracts of 53rdAnnual Internat SEG Mtg,1983,382-385

[2] Baysal E,Kosloff D D,Sherwood J W C.Reverse time migration[J].Geophysics,1983,48(11):1514-1524

[3] McMechan G A.Migration by extrapolation of time-dependent boundary values[J].Geophysical Prospecting,1983,31(3):412-420

[4] 薛東川,王尚旭.波動(dòng)方程有限元疊前逆時(shí)偏移[J].石油地球物理勘探,2008,43(1):17-21 Xue D C,Wang S X.Finite element wave equation prestack reverse time migration[J].Oil Geophysical Prospecting,2008,43(1):17-21

[5] 劉紅偉,李博,劉洪.地震疊前逆時(shí)偏移高階有限差分算法及GPU實(shí)現(xiàn)[J].地球物理學(xué)報(bào),2010,53(7):1725-1733 Liu H W,Li B,Liu H.The algorithm of high order finite difference pre-stack reverse time migration and GPU implementation[J].Chinese Journal of Geophysics(in Chinese),2010,53(7):1725-1733

[6] 郭念民,吳國(guó)忱.基于PML邊界的變網(wǎng)格高階有限差分聲波方程逆時(shí)偏移[J].石油地球物理勘探,2012,47(2):256-265 Guo N M,Wu G C.High-order finite difference method in reverse-time migration with variable grids based on PML boundary condition[J].Oil Geophysical Prospecting,2012,47(2):256-265

[7] 劉紅偉,劉洪,鄒振,等.地震疊前逆時(shí)偏移中的去噪與存儲(chǔ)[J].地球物理學(xué)報(bào),2010,53(9):2171-2180 Liu H W,Liu H,Zou Z,et al.The problems of denoise storage in seismic reverse time migration[J].Chinese Journal of Geophysics(in Chinese),2010,53(9):2171-2180

[8] Wyatt K D.Synthetic vertical seismic profile[J].Geophysics,1981,46(6):880-891

[9] Lee M W,Balch A H.Computer processing of vertical seismic profile data[J].Geophysics,1983,48(3):272-287

[10] Dillon P B.Vertical seismic profile migration using the Kirchhoff integral[J].Geophysics,1988,53(6):786-799

[11] 劉清林,駱毅.VSP資料的疊前波動(dòng)方程逆時(shí)偏移[J].石油地球物理勘探,1989,24(3):281-289 Liu Q L,Luo Y.Prestack reverse time migration of VSP data[J].Oil Geophysical Prospecting,1989,24(3):281-289

[12] 朱金明,嚴(yán)俊華.VSP逆時(shí)偏移處理[J].石油地球物理勘探,1991,26(5):564-570 Zhu J M,Yan J H.VSP reverse-time migration[J].Oil Geophysical Prospecting,1991,26(5):564-570

[13] 朱金明,董敏煜,李承楚.VSP的雙程無反射波動(dòng)方程逆時(shí)偏移[J].石油地球物理勘探,1989,24(3):256-270 Zhu J M,Dong M Y,Li C C.VSP reverse-time migration using two-way nonreflection wave equation[J].Oil Geophysical Prospecting,1989,24(3):256-270

[14] Zhang J H,Yao Z X.Optimized finite-difference operator for broadband seismic wave modeling[J].Geophysics,2012,77(1):A13-A18

[15] Clapp R G.Reverse time migration with random boundaries[J].Expanded Abstracts of 79thAnnual Internat SEG Mtg,2009,2809-2813

[16] 郭念民,馮雪梅,李海山.高階拉普拉斯算子逆時(shí)偏移低頻噪聲去除方法[J].石油物探,2013,52(6):642-649 Guo N M,Feng X M,Li H S.Research on higher order Laplacian operator denoising method in reverse-time migration[J].Geophysical Prospecting for Petroleum,2013,52(6):642-649

[17] 劉詩竹,石穎,郭雪豹,等.VSP數(shù)據(jù)逆時(shí)偏移方法研究[J].地球物理學(xué)進(jìn)展,2014,29(5):2211-2218 Liu S Z,Shi Y,Guo X B,et al.Investigation on reverse-time migration of VSP data[J].Progress in Geophysics(in Chinese),2014,29(5):2211-2218

[18] 鄭阿奇.QT4開發(fā)實(shí)踐[M].北京：電子工業(yè)出版社,2011:30-48 Zheng A Q.QT4 development practices[M].Beijing:Publishing House of Electronics Industry,2011:30-48

[19] 丁仁偉,李振春,仝兆岐.基于Qt的共聚焦點(diǎn)技術(shù)地震資料處理系統(tǒng)研發(fā)[J].石油物探,2009,48(2):199-205 Ding R W,Li Z C,Tong Z Q.Research and development of seismic data processing system based on CFP technology by utilizing Qt graphical user interface library[J].Geophysical Prospecting for Petroleum,2009,48(2):199-205

(編輯：陳 杰)

Qt based software development of reverse time migration for walkaway VSP seismic data

Wang Weihong1,Ke Xuan1,Guo Xuebao1,Zhang Yingying1,Yu Peng2

(1.CollegeofEarthSciences,NortheastPetroleumUniversity,Daqing163318,China;2.Oil&GasExplorationDepartment,NortheastOil&GasBranch,SINOPEC,Changchun130062,China)

The geophones in VSP geometry are located down in the borehole,therefore VSP seismic data has high S/N and with rich wavefield information.Based on two-way acoustic wave raypath,the prestack reverse time migration (RTM) has high precision for imaging,and provides important data for the study of complex subsurface structure near the borehole and fluid properties in targeted formation.According to the algorithm of reverse time migration for walkaway VSP seismic data,the complete RTM processing flow of VSP data was designed and the VSP RTM software was preliminarily developed,which is developed by Qt language as well as in Linux operating system.Two types of core algorithm systems,CPU and GPU can be installed in different hardware.The test of model data demonstrated that the developed software is characterized by easy operating,good transferring and high imaging precision.

VSP,reverse time migration,Qt,software interface

2015-01-23;改回日期：2015-05-21。

王維紅(1975—),男,副教授,主要研究方向?yàn)榈卣鹳Y料數(shù)字處理。

國(guó)家自然科學(xué)基金面上項(xiàng)目(41474118)、國(guó)家高技術(shù)研究發(fā)展計(jì)劃(863計(jì)劃)項(xiàng)目(2012AA061202)聯(lián)合資助。

P631

A

1000-1441(2015)04-0452-07

10.3969/j.issn.1000-1441.2015.04.012