沉積相數(shù)據(jù)可視化的研究與應(yīng)用

尚福華，李惠娟

(東北石油大學(xué) 計(jì)算機(jī)與信息技術(shù)學(xué)院，黑龍江 大慶 163318)

沉積相數(shù)據(jù)可視化的研究與應(yīng)用

尚福華，李惠娟

(東北石油大學(xué) 計(jì)算機(jī)與信息技術(shù)學(xué)院，黑龍江 大慶 163318)

為了將各種沉積相數(shù)據(jù)的特征以及油水井的布局疊加展示，以便更加直觀地表現(xiàn)出剩余油層的分布情況，從而為油田的開發(fā)決策者提供決策依據(jù),研究了沉積相數(shù)據(jù)可視化的一些方法，因?yàn)榈刭|(zhì)數(shù)據(jù)可視化技術(shù)是建設(shè)數(shù)字油田的關(guān)鍵步驟。以地質(zhì)數(shù)據(jù)成圖平臺(tái)為依托，結(jié)合常見的數(shù)據(jù)可視化技術(shù)、三角剖分技術(shù)對(duì)沉積相數(shù)據(jù)進(jìn)行可視化的研究與應(yīng)用。主要介紹了數(shù)據(jù)可視化的發(fā)展、實(shí)現(xiàn)的步驟、關(guān)鍵核心技術(shù)以及三角剖分的主要思想、步驟、實(shí)現(xiàn)的方法等。其中可視化映射的方法主要以三角網(wǎng)格剖分為主，也描述了數(shù)據(jù)的預(yù)處理、凸包建立步驟以及三角網(wǎng)優(yōu)化等算法。系統(tǒng)以Java2D繪圖技術(shù)為核心，可視化映射速度快、精準(zhǔn)度高、人機(jī)交互方便快捷，達(dá)到了預(yù)想的實(shí)驗(yàn)?zāi)康摹?/p>

數(shù)據(jù)可視化；沉積相；仿射變換；三角剖分

1 概 述

數(shù)據(jù)可視化(Data Visualization)是利用計(jì)算機(jī)圖形學(xué)和圖像處理技術(shù)，將數(shù)據(jù)轉(zhuǎn)換成圖形或圖像在屏幕上顯示出來，并進(jìn)行交互處理的理論、方法和技術(shù)[1]。數(shù)據(jù)可視化的出現(xiàn)為人們對(duì)數(shù)據(jù)的研究方式產(chǎn)生了極大的影響，相比于大型數(shù)據(jù)集中的數(shù)據(jù)以及冗長(zhǎng)繁多的數(shù)據(jù)離鏈表，現(xiàn)在人們能直觀地觀察數(shù)據(jù)展示出來的各種圖表，從而更加方便用戶去查看數(shù)據(jù)的走向、位置分布、分布的疏密度以及不同位置的不同形態(tài)或者物件等等。數(shù)據(jù)可視化的核心思想就是用圖形元素(點(diǎn)、線、面)來表示數(shù)據(jù)關(guān)系。隨著數(shù)據(jù)倉(cāng)庫(kù)技術(shù)、數(shù)據(jù)挖掘技術(shù)、網(wǎng)絡(luò)技術(shù)以及電子商務(wù)技術(shù)等的發(fā)展，數(shù)據(jù)可視化技術(shù)的使用范圍也越來越廣泛，帶來的效益也越來越大。

所謂的沉積相數(shù)據(jù)可視化，就是以油田地質(zhì)數(shù)據(jù)中的沉積相數(shù)據(jù)為研究對(duì)象，通過可視化技術(shù)將沉積相數(shù)據(jù)以二維圖幅的形式展示出來，實(shí)現(xiàn)由抽象不可見到直觀可見的過程。將各種沉積相數(shù)據(jù)的特征以及油水井的布局疊加展示，更加直觀地表現(xiàn)出剩余油層的分布情況，從而對(duì)下一步油田開發(fā)起到?jīng)Q策支持的作用。

由于計(jì)算機(jī)技術(shù)的快速發(fā)展，使得數(shù)據(jù)可視化的能力越來越強(qiáng)大，也越來越能滿足用戶的需求，尤其是現(xiàn)在的可視化技術(shù)提供了大量的交互功能，為用戶的使用帶了極大的便利，用戶可以根據(jù)自己的喜好對(duì)感興趣的全部或局部數(shù)據(jù)進(jìn)行各種角度的交互處理。因此數(shù)據(jù)可視化的應(yīng)用領(lǐng)域也得到了大幅度的拓寬，除了原先的生物醫(yī)藥、交通、氣象預(yù)報(bào)、海洋領(lǐng)域，又增加了數(shù)據(jù)挖掘、社交、金融、通信以及網(wǎng)絡(luò)等等[2]。文中著重講述了油田地質(zhì)數(shù)據(jù)中的沉積相的可視化過程。

一般用三個(gè)較為重要的步驟來講述數(shù)據(jù)可視化技術(shù)：數(shù)據(jù)預(yù)處理(源數(shù)據(jù)轉(zhuǎn)化到DB(DataBase)中，包括數(shù)據(jù)選樣、數(shù)據(jù)變換等)、可視化映射(數(shù)據(jù)到可視化結(jié)構(gòu)的轉(zhuǎn)換以及可視化結(jié)構(gòu)到視圖的轉(zhuǎn)換，包括坐標(biāo)映射、幾何映射等)、可視化交互(建立人與系統(tǒng)之間的一種交流方式，包括選擇交互、過濾交互、細(xì)節(jié)交互等)?？梢暬牟襟E模型如圖1所示[3-4]。

圖1 可視化的步驟流程

2 數(shù)據(jù)預(yù)處理

沉積相是對(duì)單層砂體規(guī)模、展布情況、接觸連通情況以及儲(chǔ)層成因的詳細(xì)刻畫，是對(duì)垂直方向上的非均勻物質(zhì)演化規(guī)律的有效反映。在該系統(tǒng)中，沉積相數(shù)據(jù)可視化映射主要用圖形學(xué)中最基本的點(diǎn)、線和面構(gòu)成。其中用點(diǎn)來說明不同的井位信息，指的是對(duì)于不同的井別(油井和水井)應(yīng)該按照不同的映射方式進(jìn)行區(qū)別展示，并根據(jù)不同的坐標(biāo)數(shù)據(jù)確定映射位置，并且確定該井所在位置包含的其他有關(guān)沉積相形成的數(shù)據(jù)信息，例如相別等。另外相別也是沉積相數(shù)據(jù)可視化的基本數(shù)據(jù)，根據(jù)沉積相數(shù)據(jù)來形成沉積相帶線，沉積相帶線是沉積相數(shù)據(jù)可視化的重要組成部分，由于沉積相是對(duì)單層沉積物質(zhì)的數(shù)字語言描述，因此比較適合用二維平面的角度來反映當(dāng)前區(qū)域具有的沉積規(guī)律。閉合的相帶線所圍成的面描述了映射范圍中具體地層單元的該相別的空間分布情況，顏色代表的是不同相別之間的區(qū)分，需要用圖例來展示。

系統(tǒng)用到的數(shù)據(jù)是大慶油田中獲取到的真實(shí)數(shù)據(jù)，通過專家對(duì)數(shù)據(jù)進(jìn)行處理，形成一系列的數(shù)據(jù)文件，系統(tǒng)對(duì)這些文件進(jìn)行解編處理，建立專門的數(shù)據(jù)庫(kù)，將坐標(biāo)點(diǎn)以BLOB的形式存儲(chǔ)于數(shù)據(jù)表中，再針對(duì)這些數(shù)據(jù)進(jìn)行坐標(biāo)轉(zhuǎn)換等。從某種程度上來說，采取合適的坐標(biāo)系轉(zhuǎn)換方法能夠簡(jiǎn)化沉積微相的可視化過程。由于該系統(tǒng)主要針對(duì)大慶油田中數(shù)字油田的建立而設(shè)計(jì)，而且大慶地區(qū)的井網(wǎng)排列有一定程度的傾角，這就為井位點(diǎn)的排序、井網(wǎng)的剖分帶來了諸多不便。為了便于可視化映射，將原始井網(wǎng)的區(qū)塊形態(tài)進(jìn)行逆時(shí)針旋轉(zhuǎn)等稍微的轉(zhuǎn)變，變?yōu)檩^為規(guī)則的形態(tài)[5]。

經(jīng)過大量的資料查閱，系統(tǒng)采用仿射變換處理坐標(biāo)變換相關(guān)操作，其中包括井位屏幕坐標(biāo)變換、圖幅的平移和縮放等。

仿射變換的一般形式為：

矩陣形式為：

3 數(shù)據(jù)點(diǎn)剖分映射

二維數(shù)據(jù)進(jìn)行面幾何可視化映射時(shí)主要用到的是網(wǎng)格剖分技術(shù)，由于普通的三角網(wǎng)無法自動(dòng)解決地質(zhì)數(shù)據(jù)二維可視化中的一系列問題，因此在普通的三角網(wǎng)剖分上提出一種基于改進(jìn)約束三角網(wǎng)的地質(zhì)二維數(shù)據(jù)可視化方法[9-11]。

三角網(wǎng)建立的基本思路是：掃描數(shù)據(jù)點(diǎn)進(jìn)行繪制凸包；采用逐點(diǎn)插入法構(gòu)建初始的三角網(wǎng)；將離散點(diǎn)插入到三角網(wǎng)中進(jìn)行三角網(wǎng)的局部修改；為了使剖分的三角形效率和精準(zhǔn)度更高，再對(duì)三角網(wǎng)進(jìn)行局部?jī)?yōu)化；對(duì)凸包進(jìn)行合并處理。

(1)掃描離散點(diǎn)構(gòu)建凸包：凸包在可視化映射中的幾何映射中有廣泛應(yīng)用，常常是其他幾何形態(tài)構(gòu)建進(jìn)行所必要的圖形基礎(chǔ)，也是進(jìn)行可視化映射的必要步驟，并且在很多實(shí)際的科研項(xiàng)目中，難點(diǎn)的解決都需要借助凸包的構(gòu)建作為工具。凸包就是將二維平面上所用點(diǎn)集構(gòu)成的最外層的一個(gè)凸多邊形[12-15]。相當(dāng)于在一個(gè)木板上訂上很多圖釘，用一條橡皮筋將所有的圖釘圈起來，這個(gè)橡皮筋就相當(dāng)于是所有點(diǎn)的凸包。遍歷繪制中遇到的所有點(diǎn)，將位置坐標(biāo)相同的點(diǎn)合并成一個(gè)點(diǎn)，并且找到x(橫坐標(biāo))和y(縱坐標(biāo))方向上的極大值與極小值。以找出來的這四個(gè)點(diǎn)(特殊情況下其中的兩點(diǎn)會(huì)重合，形成的初始凸包為三角形)連成線形成初始的凸包四邊形，用相鄰的兩點(diǎn)按逆時(shí)針組成有向線段，搜索其右側(cè)能與之形成最大面積的三角形的另一個(gè)點(diǎn)，將該點(diǎn)記錄下來，重復(fù)這個(gè)操作直到鏈表中沒有其他的坐標(biāo)點(diǎn)，這個(gè)時(shí)候初始的凸包就建立完成了。三角形面積可以用鏈表中點(diǎn)的坐標(biāo)通過行列式來表示：

圖2 凸包三角剖分

(3)逐步向區(qū)域內(nèi)插入新點(diǎn)：在建立的初始三角網(wǎng)中還存在很多離散的點(diǎn)，如圖3中點(diǎn)G。通過逐步向區(qū)域中插入離散點(diǎn)，可以將點(diǎn)G，H按照一定的規(guī)則插入到三角網(wǎng)中。首先找出含有點(diǎn)G的三角形即

圖3 插入離散點(diǎn)

ΔBFD，然后尋找與ΔBFD相接臨的三角形，看其外接圓是否包含點(diǎn)G，若包含如ΔABF，則將ΔBFD與ΔABF形成一個(gè)區(qū)域，然后刪除掉兩個(gè)三角形的公共邊BF，最后將點(diǎn)G與形成區(qū)域的所有頂點(diǎn)連接起來形成AG、BG、DG、FG四條邊，此時(shí)已完成將離散點(diǎn)G插入到三角網(wǎng)中的全部過程。

(4)優(yōu)化三角網(wǎng)：使用Lawson局部?jī)?yōu)化算法，凡是經(jīng)過此算法進(jìn)行局部?jī)?yōu)化的三角網(wǎng),它在該局部區(qū)域內(nèi)的三角形都會(huì)變成Delaunay三角形。Delaunay三角形有唯一性(以任何一個(gè)點(diǎn)為起點(diǎn)都不會(huì)影響最終的三角網(wǎng)結(jié)構(gòu))、局部性(不影響之前剖分好的三角網(wǎng)結(jié)構(gòu))、最優(yōu)性(形成的三角網(wǎng)是最優(yōu)的)、最規(guī)則性(最小角進(jìn)行排序)、最近性(保證三角網(wǎng)的均勻)、存在凸多邊形的外殼(網(wǎng)格的最外層界限所有點(diǎn)均為凸點(diǎn))、空?qǐng)A特性(任一外接圓中不含有別的點(diǎn))以及對(duì)偶性等。Lawson優(yōu)化算法的主要實(shí)施方法是判斷擁有公共邊的兩個(gè)三角形是否需要重新設(shè)定公共邊，遵循的原則是空外接圓和最小角最大的原則。實(shí)現(xiàn)起來就是判斷一個(gè)三角形的外接圓是否包含另一個(gè)三角形的另一個(gè)頂點(diǎn)，若包含(如圖4中的點(diǎn)C在ΔABD的外接圓內(nèi))，則將BD邊刪除，連接AC，即將對(duì)角線進(jìn)行交換重繪。這種方法易于理解且便于實(shí)施到程序中。

圖4 局部?jī)?yōu)化三角網(wǎng)

(5)合并凸包：根據(jù)兩個(gè)凸包投影的關(guān)系，先找出兩個(gè)凸包的支撐點(diǎn)，根據(jù)凸包的合并原理，利用公共支撐線和支撐點(diǎn)在兩個(gè)凸包之間生成新的三角形，同時(shí)為了保證兩個(gè)凸包合并時(shí)沒有交叉，也要對(duì)原有凸包里面的公共支撐線進(jìn)行新的三角剖分。

4 應(yīng)用實(shí)例

該軟件是針對(duì)實(shí)現(xiàn)數(shù)字油田而設(shè)計(jì)的，用于實(shí)現(xiàn)“同一平臺(tái)的一口井計(jì)劃”。系統(tǒng)首先將沉積相可視化的映射結(jié)果結(jié)合井位坐標(biāo)圖作為底層展示。系統(tǒng)主要以Java語言開發(fā)為核心，使用了先進(jìn)的RIA(RichInternetApplications)豐富互聯(lián)網(wǎng)技術(shù)，前端采用.NET網(wǎng)頁(yè)框架和Applet網(wǎng)頁(yè)嵌入技術(shù)，后端主要采用Java2D繪圖技術(shù)。在繪圖過程中采用Java2DAPI，通過一個(gè)Web瀏覽器使用Java插件軟件，使用JavaWebStart(Java運(yùn)行時(shí)環(huán)境中較新的一門技術(shù))來部署應(yīng)用程序。

該系統(tǒng)通過仿射變換將原有的沉積相數(shù)據(jù)坐標(biāo)點(diǎn)轉(zhuǎn)換成用于剖分三角形的坐標(biāo)點(diǎn)集，通過三角剖分將點(diǎn)集剖分成三角形區(qū)域塊，將相同的沉積相通過曲線平滑形成相帶線。圖5為結(jié)果展示。

圖5 可視化映射展示圖

5 結(jié)束語

文中詳細(xì)研究了數(shù)據(jù)可視化的發(fā)展及應(yīng)用技術(shù)，針對(duì)數(shù)據(jù)可視化映射的研究較為詳細(xì)，將沉積相數(shù)據(jù)進(jìn)行可視化映射的意義非常重大，并且文中算法具有以下特點(diǎn)：

(1)思路簡(jiǎn)單，易于實(shí)現(xiàn)，算法較成熟；

(2)算法的效率較高，將時(shí)間和空間復(fù)雜度進(jìn)行了較好的搭配；

(3)針對(duì)不同的數(shù)據(jù)做了不同的數(shù)據(jù)預(yù)處理，簡(jiǎn)化了后面可映射難度。

另外，該系統(tǒng)的可視化映射不僅速度快并且精準(zhǔn)度高，大大提高了以往地質(zhì)數(shù)據(jù)可視化的效率。系統(tǒng)又對(duì)映射結(jié)果進(jìn)行了交互處理，主要包括選擇交互、過濾交互以及抽象交互等。選擇交互可以顯示數(shù)據(jù)的一部分，對(duì)該數(shù)據(jù)進(jìn)行操作；過濾交互是對(duì)感興趣的部分進(jìn)行約束交互，不滿足的部分不予考慮或被隱藏；縮放交互是對(duì)局部或者整體進(jìn)行縮放，抽象交互及在過濾交互的基礎(chǔ)上對(duì)約束數(shù)據(jù)進(jìn)行詳細(xì)細(xì)節(jié)的抽象操作。建立了人與系統(tǒng)之間的一種交流方式，能更好地幫助用戶理解和揭示數(shù)據(jù)信息。

[1] 韓麗娜.數(shù)據(jù)可視化技術(shù)及其應(yīng)用展望[J].煤礦現(xiàn)代化,2005(6):39-40.

[2] 楊彥波,劉 濱,祁明月.信息可視化研究綜述[J].河北科技大學(xué)學(xué)報(bào),2014,35(1):91-102.

[3] 王媛媛,丁 毅,孫媛媛,等.數(shù)據(jù)可視化技術(shù)的實(shí)現(xiàn)方法研究[J].現(xiàn)代電子技術(shù),2007,30(4):71-74.

[4]LeeMD,ButaviciusMA,ReillyRE.Visualizationsofbinarydata:acomparativeevaluation[J].InternationalJournalofHuman-ComputerStudies,2003,59(3):569-602.

[5] 尚福華,冀 宇,曹茂俊.改進(jìn)約束三角網(wǎng)在沉積相帶圖中的應(yīng)用研究[J].計(jì)算機(jī)技術(shù)與發(fā)展,2014,24(11):242-245.

[6] 宋麥玲,陳云亮.基于仿射變換的圖象定位方法[J].電腦知識(shí)與技術(shù),2007(14):478-480.

[7] 張 偉,任仙怡,張桂林,等.基于對(duì)數(shù)極坐標(biāo)變換和仿射變換的目標(biāo)識(shí)別定位方法[J].中國(guó)圖象圖形學(xué)報(bào),2006,11(9):1255-1259.

[8] 孫 燾,王秀坤,邵 剛,等.二維點(diǎn)模式圖像的仿射變換配準(zhǔn)[J].計(jì)算機(jī)輔助設(shè)計(jì)與圖形學(xué)學(xué)報(bào),2005,17(7):1497-1503.

[9] 徐道柱,劉海硯.大量約束邊條件下Delaunay三角網(wǎng)的快速生成[J].測(cè)繪工程,2007,16(3):6-10.

[10] 鄒永貴,張 濤.改進(jìn)的平面域Delaunay三角網(wǎng)生成算法[J].計(jì)算機(jī)工程與應(yīng)用,2013,49(20):171-174.

[11] 徐 旭,李 源,陳學(xué)工.一種基于插入法的Delaunay三角網(wǎng)生成算法[J].電腦與信息技術(shù),2010,18(4):29-31.

[12] 寧化展,徐炳喜,田茂義,等.帶有邊界條件的Delaunay三角網(wǎng)生成算法的研究與實(shí)現(xiàn)[J].全球定位系統(tǒng),2010,35(4):39-42.

[13] 蔡 強(qiáng),楊 欽,陳其明.地質(zhì)結(jié)構(gòu)重疊域的限定Delaunay三角剖分研究[J].計(jì)算機(jī)輔助設(shè)計(jì)與圖形學(xué)學(xué)報(bào),2004,16(6):766-771.

[14]ChewLP.ConstrainedDelaunaytriangulations[J].Algorith-mica,1989(4):97-108.

[15]LeeDT,SchachterBJ.TwoalgorithmsforconstructingaDelaunaytriangulation[J].InternationalJournalofComputer&InformationSciences,1980,9(3):219-242.

Research and Application of Data Visualization of Sedimentary Facies

SHANG Fu-hua，LI Hui-juan

(School of Computer and Information Technology,Northeast Petroleum University,Daqing 163318,China)

In order to overlay display the characteristics of sedimentary facies data and distribution of oil and water wells，as well as be more convenient and intuitive to show the distribution of remaining oil and helpful for the decision makers of foundation for decision making,some methods of sedimentary facies data visualization are studied because visualization of geological data is a key step in the construction of digital oil field.Based on the geological data mapping platform,data visualization technology and triangulation technology is applied to research the distribution characteristic and data rules of sedimentary facies.In this paper,the development,implementation steps and key core technology of data visualization and main ideas,steps and implementation method of triangulation technology are introduced.The method of visualization mapping is mainly based on triangulation and the methods of the data preprocessing,steps of establishing convex closure and triangular mesh optimization is discussed.The system relied on takes Java2D as the core of rendering technology.Its high speed of visual mapping,high accuracy and convenient man-machine interaction make the system very excellent,which achieves the expected objective.

data visualization;sedimentary facies;affine transformation;triangulation technology

2015-12-13

2016-04-06

時(shí)間：2016-08-01

國(guó)家自然科學(xué)基金資助項(xiàng)目(61170132)；國(guó)家重大專項(xiàng)(2011ZX05020-007)

尚福華(1962-)，男，教授，博士(后)，研究方向?yàn)槿斯ぶ悄堋?shù)據(jù)挖掘、計(jì)算機(jī)理論與方法；李惠娟(1989-)，女，碩士研究生，研究方向?yàn)橛?jì)算機(jī)應(yīng)用技術(shù)。

http://www.cnki.net/kcms/detail/61.1450.TP.20160801.0909.062.html

TP39

A

1673-629X(2016)08-0191-04

10.3969/j.issn.1673-629X.2016.08.041