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

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

0 前言

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

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

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

1 微震數據分析

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

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

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

BM浮點數格式。但這種浮點數格式不能直接在微機上使用,需要轉換為在微機上通用的IEEE浮點數格式。

在采集微震數據時,常采用三分量檢波器。每一個檢波器分量采集的數據對應一個地震道,即一個檢波器對應三個地震道。因為微震事件能量低頻率高,其采樣間隔比常規(guī)地震勘探更小,通常為0.25 ms。

2 LOD的應用

細節(jié)層次(Level of Detail,LOD)技術根據觀察點位置的變化而選擇不同的精度模型進行成像。該算法首先由Clark于1976年提出,最初是用來簡化采樣密集的多面體網狀物體數據結構。其基本思想是:在不影響視覺效果的前提下,對場景中每一原始物體建立不同逼近精度的幾何模型,在繪制時,根據一定的標準,選擇適當的層次模型來表示物體,從而簡化場景的描述[3]。

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

(1)靜態(tài)LOD。在預處理過程中,產生一個物體的幾個離散的不同細節(jié)層級模型,在繪制時根據特定的標準,選擇合適的細節(jié)層次模型來表示物體。

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

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

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

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

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

3 內存映射

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

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

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

//建立映射

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

4 實例分析

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

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

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

為了準確拾取微震事件的初至時間,可對地震道波形進行局部放大顯示,以方便拾取初至時間(見圖4)。

在圖4中,可以準確地找到微震事件的初至時間(S波的初至時間)。在圖4中標示的初至時間為1 295.00 ms。由于圖4中的S波的能量較強,同時由于程序處理的原因,所以P波的波形在圖4中不明顯。

圖4 拾取初至時間Fig.4 Picking up the first arrival

5 結論

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

[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技術及其算法[J].軟件導刊,2008,7(10):48.

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

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

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

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