王 鵬,桂志先,謝宋雷,王朝鋒,張建萍

(1.油氣資源與勘探技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室(長江大學(xué)) 長江大學(xué) 地球物理與石油資源學(xué)院,湖北荊州 434023;2.中原油田地球物理測井公司,河南濮陽 457001)

0 前言

壓裂自從上個(gè)世紀(jì)七十年代在美國堪薩斯州西部Hugoton氣田首次使用以來,已成為一種常見的提高油氣采收率的油藏工程方法[1]。壓裂作業(yè)會誘生微震,可對微震信號進(jìn)行采集,處理和解釋,還可對壓裂作業(yè)進(jìn)行監(jiān)測,并對壓裂效果進(jìn)行評價(jià)。壓裂作業(yè)的時(shí)間可長達(dá)幾個(gè)小時(shí),甚至幾十個(gè)小時(shí)。由于微震事件能量弱,信號頻率高,易被地層吸收,所以這些特點(diǎn)決定了微震數(shù)據(jù)采集時(shí)間較長,采樣頻率高,最終將導(dǎo)致得到的微震數(shù)據(jù)相當(dāng)龐大。

利用井下多個(gè)三分量檢波器接收微震信號,并在信號中拾取出微震事件的初至?xí)r間。對微震事件做極化分析,就可以確定微震事件的傳播方向。利用測井和地震資料,可建立速度模型。結(jié)合以上數(shù)據(jù),利用多種反演方法(例如:縱、橫波時(shí)差震源定位法)計(jì)算出微震震源的空間位置,從而可以畫出裂縫的空間圖像。從裂縫的空間圖像可以直觀地看出裂縫的發(fā)育方向和范圍,再與壓裂后試井結(jié)果相比對,就能對壓裂效果進(jìn)行評價(jià),進(jìn)而改進(jìn)壓裂作業(yè)。

在壓裂作業(yè)的各個(gè)階段,微震事件的發(fā)生率有明顯的差異。相對于采集得到的數(shù)據(jù),實(shí)際反演需要的數(shù)據(jù)只占很少的一部份。這樣,就需要用快捷的方法來讀取和顯示采集的數(shù)據(jù),并從中找到微震事件。作者在本文中,以VC6.0作為編程工具,并利用LOD和內(nèi)存映射來實(shí)現(xiàn)快速讀取和顯示微震數(shù)據(jù)。

1 微震數(shù)據(jù)分析

經(jīng)采集得到的微震數(shù)據(jù)解編后,使用地震勘探中廣泛采用的SEGY數(shù)據(jù)格式來存儲。SEGY數(shù)據(jù)格式由卷頭和道信息二部份組成,其中道信息又分為道頭和道數(shù)據(jù)二部份(見圖1)。

圖1 SEGY文件格式Fig.1 The structure of SEGY

卷頭包括3 200字節(jié)的字符信息和400字節(jié)的二進(jìn)制信息。400字節(jié)的二進(jìn)制信息內(nèi)容,將會影響到整個(gè)SEGY文件,這些信息對文件的數(shù)據(jù)處理至關(guān)重要,特別是采樣間隔,道長度和編碼格式[2]。240字節(jié)的道頭,包含了單道的基本信息和道間的差異信息,道數(shù)據(jù)包含單道的所有采樣點(diǎn)數(shù)據(jù)信息。對于多數(shù)SEGY數(shù)據(jù)文件,其道數(shù)據(jù)采用I

BM浮點(diǎn)數(shù)格式。但這種浮點(diǎn)數(shù)格式不能直接在微機(jī)上使用,需要轉(zhuǎn)換為在微機(jī)上通用的IEEE浮點(diǎn)數(shù)格式。

在采集微震數(shù)據(jù)時(shí),常采用三分量檢波器。每一個(gè)檢波器分量采集的數(shù)據(jù)對應(yīng)一個(gè)地震道,即一個(gè)檢波器對應(yīng)三個(gè)地震道。因?yàn)槲⒄鹗录芰康皖l率高,其采樣間隔比常規(guī)地震勘探更小,通常為0.25 ms。

2 LOD的應(yīng)用

細(xì)節(jié)層次(Level of Detail,LOD)技術(shù)根據(jù)觀察點(diǎn)位置的變化而選擇不同的精度模型進(jìn)行成像。該算法首先由Clark于1976年提出,最初是用來簡化采樣密集的多面體網(wǎng)狀物體數(shù)據(jù)結(jié)構(gòu)。其基本思想是:在不影響視覺效果的前提下,對場景中每一原始物體建立不同逼近精度的幾何模型,在繪制時(shí),根據(jù)一定的標(biāo)準(zhǔn),選擇適當(dāng)?shù)膶哟文Ｐ蛠肀硎疚矬w,從而簡化場景的描述[3]。

LOD模型的實(shí)現(xiàn)方式:

(1)靜態(tài)LOD。在預(yù)處理過程中,產(chǎn)生一個(gè)物體的幾個(gè)離散的不同細(xì)節(jié)層級模型,在繪制時(shí)根據(jù)特定的標(biāo)準(zhǔn),選擇合適的細(xì)節(jié)層次模型來表示物體。

(2)動態(tài)LOD。在動態(tài)LOD算法中生成一個(gè)數(shù)據(jù)結(jié)構(gòu)。在繪制時(shí),可以從這個(gè)數(shù)據(jù)結(jié)構(gòu)中抽取所需要的細(xì)節(jié)層次模型,從而得到大量不同的分辨率,并且是可以連續(xù)變化的[4]。

微震波型有體波(P波和S波二種)和導(dǎo)波二類。在三分量檢波器記錄上,每個(gè)分量上P波和S波成對出現(xiàn),并且三個(gè)分量上的P波和S波的波至?xí)r間分別相同[5]。在每個(gè)分量上,P波后面緊跟一個(gè)S波,S波比P波的波幅更大(見圖2)。在噪音干擾不強(qiáng)的情況下,微震事件在地震記錄中是比較明顯的。同時(shí)根據(jù)大量的文獻(xiàn)分析,單個(gè)微震的持續(xù)時(shí)間從十幾毫秒到三百毫秒不等,這與其能量強(qiáng)弱直接相關(guān),也與檢波器到震源距離有關(guān)。

圖2 微震波形Fig.2 The wavefor m ofmicroseisms

以持續(xù)時(shí)間為50 ms的微震事件為例,微震事件對應(yīng)的采樣點(diǎn)個(gè)數(shù)為200。這是因?yàn)槲⒄鹗录腟波的振幅比P波振幅大,所以在地震記錄中比較明顯。同時(shí),微震事件對應(yīng)的采樣點(diǎn)多,所以利用LOD技術(shù)來顯示微震數(shù)據(jù)是可行的。依據(jù)靜態(tài)LOD模型,將可分為三個(gè)顯示等級:①每隔十個(gè)樣點(diǎn)進(jìn)行顯示。這時(shí)能夠顯示能量較強(qiáng),延續(xù)時(shí)間較長的微震事件。這些微震事件由于能量較強(qiáng),所以能較為準(zhǔn)確地提取出其初至?xí)r間;②每隔五個(gè)樣點(diǎn)進(jìn)行顯示,能進(jìn)一步發(fā)現(xiàn)其它能量較弱的微震事件。③顯示所有的樣點(diǎn)。

由于壓裂作業(yè)時(shí)間較長,采集到的微震數(shù)據(jù)往往會被依次保存到多個(gè)文件中。使用LOD就可減少讀取的數(shù)據(jù)量,加快讀取數(shù)據(jù)的速度,并且快速找到有效的微震事件。

3 內(nèi)存映射

處理大數(shù)據(jù)量SEGY文件時(shí),往往會出現(xiàn)內(nèi)存不足的問題。作者本文中,利用文件內(nèi)存映射的方法來解決這個(gè)問題。

內(nèi)存映射是W indows的一種內(nèi)存管理方法,與虛擬內(nèi)存類似[6]。建立文件內(nèi)存映射的過程大致可分為:①利用CreateFile()函數(shù)創(chuàng)建或打開一個(gè)文件內(nèi)核對象;②利用CreateFileMapping()函數(shù)創(chuàng)建一個(gè)文件映射內(nèi)核對象;③利用MapViewOfFile()函數(shù),將文件映射對象的全部或部分映射到進(jìn)程的地址空間中;④利用UnmapViewOfFile()函數(shù),撤銷文件數(shù)據(jù)的映射,并調(diào)用CloseHandle()函數(shù)關(guān)閉前面創(chuàng)建的文件映射對象和文件對象。部份代碼如下:

//創(chuàng)建文件內(nèi)存映射內(nèi)核對象

//建立映射

同時(shí),當(dāng)窗口由于任何原因需要重繪時(shí),總是先用背景色將顯示區(qū)清除,然后才調(diào)用相應(yīng)的繪圖函數(shù),而背景色往往與繪圖內(nèi)容反差很大,在短時(shí)間內(nèi),背景色與顯示圖形交替出現(xiàn),形成屏閃現(xiàn)象。為了解決這一問題,應(yīng)采用雙緩沖的方法。具體實(shí)現(xiàn)方法是:先在內(nèi)存中將圖像繪制完成,再將繪制好的圖像一次性拷貝到屏幕中。

4 實(shí)例分析

實(shí)例數(shù)據(jù)來源于油田壓裂監(jiān)測數(shù)據(jù)。監(jiān)測方式采用單井連續(xù)監(jiān)測,井下共有七個(gè)三分量檢波器。在壓裂作業(yè)開始前10 min,進(jìn)行壓裂監(jiān)測,以記錄背景噪音。壓裂監(jiān)測持續(xù)約2 h。利用LOD技術(shù)來顯示微震數(shù)據(jù),在圖像中可以清晰識別出微震事件(見圖3)。

圖3 微震事件的識別Fig.3 The recognition ofmicroseis ms event

圖3是對同一個(gè)微震事件,使用三種不同的顯示方式。圖3(a)采用常規(guī)的顯示方式,圖3(b)和圖3(c)采用了LOD技術(shù)來顯示。在圖3(a)中,雖然微震事件在圖像中很清晰,但繪制所需的數(shù)據(jù)較多,繪制的時(shí)間比較長。在繪制圖3(b)和圖3(c)時(shí),對監(jiān)測數(shù)據(jù)進(jìn)行了適當(dāng)?shù)某槿『秃喕?在繪制圖3(b)時(shí)抽取了1/5的數(shù)據(jù)量,而繪制圖3(c)時(shí)抽取了1/10的數(shù)據(jù)量。在圖3(b)和圖3(c)中,都能識別出微震事件。

為了準(zhǔn)確拾取微震事件的初至?xí)r間,可對地震道波形進(jìn)行局部放大顯示,以方便拾取初至?xí)r間(見圖4)。

在圖4中,可以準(zhǔn)確地找到微震事件的初至?xí)r間(S波的初至?xí)r間)。在圖4中標(biāo)示的初至?xí)r間為1 295.00 ms。由于圖4中的S波的能量較強(qiáng),同時(shí)由于程序處理的原因,所以P波的波形在圖4中不明顯。

圖4 拾取初至?xí)r間Fig.4 Picking up the first arrival

5 結(jié)論

微震事件的初至拾取是進(jìn)行壓裂監(jiān)測的基礎(chǔ)之一。同時(shí)在接收到的微震信號中,存在眾多的微震事件。這些微震事件的能量大小,持續(xù)時(shí)間都存在明顯的差異。這樣就需要對監(jiān)測數(shù)據(jù)進(jìn)行快速讀取,識別出其中能量較強(qiáng)的微震事件,并對這些微震事件的初至?xí)r間進(jìn)行準(zhǔn)確的拾取。通過以上分析和實(shí)例驗(yàn)證,在微震數(shù)據(jù)的讀取和顯示中,利用LOD技術(shù)和內(nèi)存映射方法,可以加快讀取速度和改進(jìn)顯示效果。同時(shí),由于微震事件的自身特點(diǎn)也會限制LOD技術(shù)的使用,如微震事件的持續(xù)時(shí)間存在明顯差異等。當(dāng)采集的微震數(shù)據(jù)信噪比不高時(shí),利用LOD技術(shù)的作用就不明顯,這時(shí)需要對采集的數(shù)據(jù)進(jìn)行必要的處理。

[1] L I YINGPING,C H CHENG,TOKSOZ M N.seismic monitoring of the growth of a hydraulic fracture zone at Fenton Hill[J],New Mexico.Geophysics,1998,63(1):120.

[2] MICHAEL W.Norris,Alank.Faichaney(Eds)SEG YrevI data Exchange for mat[Z].Tulsa,Oklahoma Society of Exploration Geophysicists,2002.

[3] 邵艷紅.LOD技術(shù)及其算法[J].軟件導(dǎo)刊,2008,7(10):48.

[4] 曾維,韓占校,朱學(xué)芳.LOD算法在3D地表模擬中的應(yīng)用研究[J].系統(tǒng)仿真學(xué)報(bào),2009,21(1):292.

[5] 梁兵,朱廣生.油氣田勘探開發(fā)中的微震監(jiān)測方法[M].北京:石油工業(yè)出版社,2004.

[6] 李曉軍.采用內(nèi)存映射文件設(shè)計(jì)超大文件數(shù)據(jù)編輯器[J].電腦編程技巧與維護(hù),2008,(9):12.

[7] 雷兵,馬在田.三維地震數(shù)據(jù)體繪制[J].物探化探計(jì)算技術(shù),2004,26(3):206.