油氣儲層連井剖面展開的算法及其實(shí)現(xiàn).docx

1、ISSN1001-9081CODENJYIIDUJournalofComputerApplications計(jì)算機(jī)應(yīng)用，2014,34(S2):305-309文章編號：1001-9081(2014)S2-0305-05油氣儲層連井剖面展開的算法及其實(shí)現(xiàn)王家華.，韓家新，翟飛飛，陳雨馨(西安石油大學(xué)計(jì)算機(jī)學(xué)院，西安710065)(通信作者電于郵箱jhwang_l)摘荽：對于各種儲層建模結(jié)果的圖形顯示，開發(fā)在面向?qū)ο蟮能浖h(huán)境下的可祝化子系統(tǒng)已經(jīng)成為當(dāng)務(wù)之急。實(shí)現(xiàn)油氣儲層連井制而圖可視化的過程中，把一個(gè)三維連井刮面展開成二維的制面圖，是關(guān)鍵的一步。給出的油氣儲屈連井制而展開的算法由以下三個(gè)部分姐成：

2、三維制面圖生成、就捋、平移。遵楣該算法，通過OpenInventor有關(guān)的類和函數(shù)，將三雄空間中的連井剖面圖展升成為二維空間中的剖面圖。油田數(shù)據(jù)的實(shí)驗(yàn)結(jié)果表明，這個(gè)算法具有直觀、效率高等優(yōu)點(diǎn)。關(guān)fit詞：連井剖面圖；可視化;OpenInventor;油氣儲層；建模中圖分類號：文獻(xiàn)標(biāo)志瑪:ASpreadingalgorithmandrealizationofwellsectionsofoilandgasreservoirWANGJiahua',HANJiaxin,ZHAIFeifei,CHENYuxin(SchoolofComputerScience,Xi'anShiyouUni

3、venity,Xi'anShaanxi710065.China)Abstract：Developingvisualsubsystemintheobject-orientedenvironment,forgraphicdisplayofreservoirhasbecomeanurgentmatter.Intheprocessofrealizingvisualizationofwellsectionforoilandgasreservoir,spreadinga3Dwellsectiontoobtaina2Dwellsectionisakeystep.Theproposedalgorith

4、mforspreading3Dwellsectionconsistsofthefollowing3parts:3Dwellsectionmaking,rotation,andtranslation.Toobeythisalgorithm,andtouserelatedclassesandfunctionsinOpenInventor,a3Dwellsectionwasspeadtoobtainthe2Dwellsection.Theexperimentalresultswithdataofaoilfielddemonstratethattheproposedalgorithmisintuiti

5、veandhighefficient.Keywords：wellsection;visualization;OpenInventor;oilandgasreservoir,modeling收稿日期：2014-03-21；修回日W：2014-06-08o基金項(xiàng)目：國家科技放大專項(xiàng)(2011ZX05032Q01)；國家自然科學(xué)驀金資助項(xiàng)目(50874091)。作者簡介:王家華(1945-),界，浙江紹興人，教授，主要研究方向:地質(zhì)統(tǒng)計(jì)學(xué)、油氣藏地質(zhì)建模算法及應(yīng)用、計(jì)算機(jī)圖形學(xué)在油田開發(fā)中的應(yīng)用；休家新(1968-),男，湖北宜昌人，教授，主要研究方向：油氣儲房地質(zhì)遂模、地質(zhì)圖件町視化；覆飛飛(19

6、89-),男，江蘇鹽城人，碩士研究生，主要研究分析：面向?qū)ο蠓椒?、地質(zhì)圖件可視化；陳雨馨(I987-),女，陜西渭南人.碩上研究生，主要研究方向：面向?qū)ο蠓椒?、地質(zhì)圖件可視化°。引言三十多年來，隨著計(jì)算機(jī)軟硬件技術(shù)、地質(zhì)統(tǒng)計(jì)學(xué)和儲層沉積學(xué)的飛速發(fā)展，油氣儲層地質(zhì)建模已經(jīng)成為國內(nèi)外油氣田開發(fā)的一項(xiàng)普遍采用的常規(guī)技術(shù)"-，：o隨著儲層建模在國內(nèi)外油氣田開發(fā)中不斷應(yīng)用，國內(nèi)外的儲層建模軟件和可視化系統(tǒng)，也有了長足的進(jìn)步4-51o特別，從油氣勘探開發(fā)的不同角度，國內(nèi)在地質(zhì)剖面圖的可視化方面做了很多工作。其中包括，利用開放圖形庫OpenGL進(jìn)行儲層剖面圖的繪制，利用面向?qū)ο蟮娜S圖形

7、可視化平臺OpenInventor進(jìn)行三維地層的可視化，以及分別基于插值算法、地理信息系統(tǒng)(GeographicInformationSystem,GIS)方法和專家知識的剖面圖可視化研究。在油氣儲層地質(zhì)建模的實(shí)施過程中，直井與直井間在空間中形成的三維連井剖面是一種重要的可視化手段。它能夠同時(shí)在垂直方向和水平方向上顯示地質(zhì)參數(shù)或油氣藏工程參數(shù)。這對于描述建模結(jié)果在三維空間中變化和油氣藏的非均質(zhì)性是特別有利的，同時(shí)對于驗(yàn)證井中數(shù)據(jù)和井間預(yù)測結(jié)果的一致性也具有直觀的作用。連井剖面可以在三維空間中生成,也可以在二維空間中生成。一般地，剖面圖是指把三維空間中生成的連井剖面圖，經(jīng)過展開、拉直后形成的二維

8、空間中所獲得的圖件。剖面圖中包含有兩部分信息。第一部分是在各口井上的測井?dāng)?shù)據(jù)。第二部分則是在各井之間的測最信息，例如地震數(shù)據(jù),也可以是通過某一種方法計(jì)算所得到的信息。但是，由于三維連井剖面是在三維空間中顯示的,且各個(gè)剖面之間按一定的角度構(gòu)成。這樣，難以在一個(gè)平面內(nèi)直觀地觀察完整的各個(gè)剖面的情況，因此將三維連井剖面展開到二維平面中，以便更加清晰地反映各剖面的特點(diǎn)。利用C+語言和工業(yè)標(biāo)準(zhǔn)OpenGL.對一個(gè)陜北油田沉積微相分布的連井削面圖進(jìn)行了可視化，如圖1所示。其中，紅色標(biāo)注的是5口井，下部的黑白圖則表示了所在剖面在井位圖中的位置。橘黃、黃、淺黃和淺藍(lán)分別代表河道、天然堤、決口扇和河漫灘。不難看

9、出，該剖面圖中各井的參數(shù)和其周圍地層的參數(shù)是一致和連續(xù)的。所用的可視化算法很簡單直接，但是其編程卻十分麻煩。在這里，為了把三維空間中的剖面展開成為二維空間中的剖面,需要把剖面中的每一個(gè)點(diǎn)的三維空間位置確定后，再求出相應(yīng)在二維空間中的位置，以及剖面上各點(diǎn)的屬性值。在圖1所示的剖面圖中，井間的沉積微相作為三維網(wǎng)格上取值的一種地質(zhì)參數(shù)，是利用井中數(shù)據(jù)，根據(jù)序貫指示模擬方法求得的。由模擬方法的特點(diǎn)所決定，井中數(shù)據(jù)和井間數(shù)據(jù)是相互匹配，并相互一致的。后來，由于測井和地震數(shù)據(jù)結(jié)合建模結(jié)果的產(chǎn)生,井間各種參數(shù)的預(yù)測成為了一種現(xiàn)實(shí)。相對利用井中數(shù)據(jù)求得的井間的柄性值，這種結(jié)果使得剖面圖的表示井間參數(shù)的可信度有

10、了極大提高，從而其應(yīng)用效果和應(yīng)用范圍有了根本的提高。本文之所以稱之為連井剖面圖，是從油/藏地質(zhì)建模的角度.出發(fā)，強(qiáng)調(diào)對于削面圖中井中數(shù)據(jù)、井間數(shù)據(jù)，以及它們兩者之間的空間相互關(guān)系。2010年以后，國內(nèi)外油氣儲層地質(zhì)建模技術(shù)的篷勃發(fā)展，商品化建模軟件系統(tǒng)的出現(xiàn)，才促進(jìn)了連井剖面圖的制作朝著實(shí)用化的目標(biāo)發(fā)展。OpenInventor是目前國際上應(yīng)用最為廣泛的面向?qū)ο蠛徒换ナ降娜S圖形開發(fā)軟件包。它是在OpenGL的基礎(chǔ)開發(fā)而成的,是一種相對獨(dú)立的窗口系統(tǒng)，可以方便地移植到使用不同窗口系統(tǒng)的硬件平臺上【5。OpenInventor具有跨平臺的能力，因此只要撰寫一份程序代碼，可編譯成在Unix/Lin

11、ux和MicrosoftWindows系統(tǒng)下的可執(zhí)行程序。OpenInventor目前支持的程序開發(fā)語言有C#、C+和Java三種，而且同時(shí)遵循調(diào)用C語言的基本炭范。這意味肴C或C+語言的編程環(huán)境下可以方便地調(diào)用OpenInventor庫中的API函數(shù)。OpenInventor是一個(gè)建立在OpenGL基礎(chǔ)t.的對象庫.開發(fā)人員可以任意使用、修改和擴(kuò)展對象庫。它擁有1300多個(gè)易于使用的類,由一系列的對象模塊組成。利用這些對象模塊,開發(fā)人員花費(fèi)敲小的編程代價(jià)就能開發(fā)出具有強(qiáng)大圖形硬件特性的應(yīng)用程序。本文算法的實(shí)現(xiàn)平臺就是OpenInventor0本文以地震波阻抗數(shù)據(jù)為例,生成連井剖面圖的算法是：

12、首先由三維地震波阻抗數(shù)據(jù)生成三維連井剖面，再通過旋轉(zhuǎn)變換和平移變換將三維連井制面展開到二維平面中，由此完成整個(gè)油氣儲層連井剖面展開的過程。1連井剖面的展開算法本文所述算法以定義在三維網(wǎng)格上的一個(gè)地震波阻抗數(shù)據(jù)體為例C1.1三雄連井剖面的生成生成三維連井剖面時(shí)，首先需要確定剖面所需的各口直井。根據(jù)被剖切的地震數(shù)據(jù)，獲得各口宜井在空間中的坐標(biāo)。確定完各口直井后，就能獲得順序的，每兩口井生成的一個(gè)面，稱之為削切面。各個(gè)剖切面在所在層位的頂界和底界之間的部分，就稱為剖面段。在圖2中,一共有4口井:井0、井1、井2和井3,生成三個(gè)割面段。將剖面段的頂界和底界沿井柱方向投影到一個(gè)平面內(nèi)（此處為YOZ平面）

13、形成的線段，稱為剖切線。在垂直于井柱的YOZ平面內(nèi)，對于直井的連井削面，一個(gè)井柱上的所有點(diǎn)在該平面的投影只是“Z平面中一個(gè)點(diǎn)。直井與直井間的一個(gè)連井剖面投影到VOZ平面就成為一條剖切線。各個(gè)剖面段的位置由坐標(biāo)原點(diǎn)到剖切線的距離和制切線法向址這兩個(gè)鼠確定。下面以井0和井1兩口直井構(gòu)成的一條剖切線為例，求出由這兩口直井構(gòu)成的剖面段的位置“生成的剖切線如圖3所示。其中井0與井1所構(gòu)成的剖切線為點(diǎn)為）和B（z”為）之間的線段。含有箭頭的線段表示的是剖切線到坐標(biāo)原點(diǎn)的距離，記為E根據(jù)井0和井1兩點(diǎn)的坐標(biāo)計(jì)算.這條削切線的直線方程為：k,kt°耳xz.yifx后=0其中4=yi-=藥-Z。根據(jù)上

14、述的剖切線方程求出坐標(biāo)原點(diǎn)到剖切線距離。:r=1為-W虹I+1剖切線的法向量是從原點(diǎn)。（0,0）到垂足M（z.y）的向量,記作"=（z,y）,如圖2中箭頭所示。因?yàn)辄c(diǎn)M在剖切線的方程上，所以有：MA,y=xx"z°xl7+Zo又點(diǎn)（0,0）和點(diǎn)（z.y）所在的直線與點(diǎn)（Zo.ro）和點(diǎn)（電,為）剖切線所在直線垂直正交，所以有：（力-y0）xy+（Z-Zo）xz=0把以上的兩式聯(lián)立，便可以計(jì)算出：=為外/&2-A）_W虹一為'-Aj/fcj+ki/kl'Z*|/&2+上/A這樣就可以求得一般情況下法向址”的值。分析以上兩個(gè)公式后，不難

15、發(fā)現(xiàn)，當(dāng)&=0或者刈=0時(shí)，該公式不成立。因此需要對法向呈的，值根據(jù)站和k2是否等于零分四種情況進(jìn)行確定。1）在用=0且刈=0的情況下，此時(shí)剖切線為一個(gè)點(diǎn)，剖切線不存在，因此法向最也不存在，說明兩口井在軸上重合。2）在如=0,虹，。的情況下：a）如果氣=0,則剖面的法向域?yàn)椤?（0,1）；b）如果乃共0,則剖面的法向量為=（O,y0）03）在知00,=0的情況下：a）如果兒=0,則剖面的法向最為=（1,0）;b）如果七X0,則剖面的法向量為“=（初,0）。4）在4R0,知第0的情況下：）如果七=0,則剖面的法向欣為=（1,-k2/k,）；b）如果七K0,則剖面的法向雖為=（尋J,織：，

16、_為+Zoh/虹kg+虹/虹）。1.2二維連井剖面的生成三維連井剖面展開為二維連井剖面的算法，分為旋轉(zhuǎn)變換和平移變換兩個(gè)步驛。1.2.1族捋變換在三維連井剖面中，前一節(jié)雖然確定了各條剖切線的位置，但剖切線也是有方向的。在進(jìn)行選井剖切時(shí)，將所選的第一口井的索引值置為0,第二口井的索引值為1,第三口井的索引值為2,，第/V口井的索引值為N-1。剖切線的正方向規(guī)定為小索引值的井到大索引值的井所指的方向。在圖3中，這個(gè)正方向即為點(diǎn)X（z0,y0）到點(diǎn)8（勾，為）所指的方向。各剖切線的角度為剖切線正方向與Z軸正方向所成的夾角。如圖4所示的角a即為剖切線的角度。通過反正切歐數(shù)atan2（doublek、,

17、double站）來計(jì)算削切線的角度a,其中4=力-yo，M=玉-初。在本文中規(guī)定始終以K軸正方向?yàn)橛沂中D(zhuǎn)的法向量，所有的剖切線都繞著該法向最旋轉(zhuǎn)a角度。經(jīng)過旋轉(zhuǎn)后，各個(gè)剖切線都與Z軸正方向平行且同向。圖4中的剖切線在經(jīng)過旋轉(zhuǎn)a角度后（如圖4為逆時(shí)針旋轉(zhuǎn)a角度），生成一條新的剖切線°圖5所示的點(diǎn)4（z。,為）到點(diǎn)8*導(dǎo)'）的線段即為旋轉(zhuǎn)后新的削切線。8（玲M）圖5旋轉(zhuǎn)變換O(0,0)Z剖切線經(jīng)過旋轉(zhuǎn)a角度后，其剖切線正方向與Z軸正方向平行且同向，因此剖切線的角度a即為旋轉(zhuǎn)角度。利用反正切函數(shù)atan2（doublek、,doublek2）計(jì)算旋轉(zhuǎn)角度，此函數(shù)的返回值為（-&q

18、uot;,們，避免了一般的反正切函數(shù)atan（）在n/2、-n/2和宣這幾個(gè)特殊點(diǎn)的討論。因此只需要對使用的參數(shù)燈和處進(jìn)行判定，分兩類確定。1）當(dāng)知=0且如=。時(shí),此時(shí)剖切線為一個(gè)點(diǎn)，剖切線不存在,因此剖切線的角度也不存在，不可以進(jìn)行旋轉(zhuǎn)。2）當(dāng)&=0和刈=0不同時(shí)成立時(shí)，此時(shí)反正切函數(shù)一定會返P1個(gè)值，因此只需要按這個(gè)值進(jìn)行旋轉(zhuǎn)即可（需要說明的是,反正切函數(shù)得到的值可能為正也可能為負(fù)或者為零,如果值為正時(shí)，就繞丫軸順時(shí)針旋轉(zhuǎn),如果值為負(fù)時(shí)，就繞Y軸逆時(shí)針旋轉(zhuǎn),如果值為零就不需要旋轉(zhuǎn)）01.2.2平移變換平移變換就是把旋轉(zhuǎn)變換后在空間中的各個(gè)剖面段平移到同一個(gè)平面內(nèi)。規(guī)定將各剖面段都平

19、移至XOZ平面。為了完成這樣的平移，該平面與NOZ平面垂直,于是只需將在YOZ平面中的各個(gè)剖面段對應(yīng)的各剖切線平移至Z軸上。仍以第一條剖切線為例，將旋轉(zhuǎn)變換后的點(diǎn)"與,為）到點(diǎn)伊（',為'）的剖切線平移后得到點(diǎn)0（0,0）到點(diǎn)8*3后術(shù)）的剖切線。如圖6所示。（0.0）b（£Q”）z8（z必）xy/仁奶）圖6平移變換其余的各條剖切線按照井的索引值從小到大的順序依次平移到Z軸上,直至所有的剖切線平移完畢。具體地,若需將點(diǎn)。.力“）和點(diǎn)（Ze，L構(gòu)成的削切線平移到Z軸上，只要計(jì)算出該剖切線的長度"=/（zz-z】了,將其平移到前面一條已經(jīng)完成平移的由點(diǎn)

20、（ZzJz）和（七，以）構(gòu)成的剖切線之后即可。這樣平移結(jié)果的放置就是從Z軸的起點(diǎn)（零點(diǎn)）開始，各口井按索引值從小到大依次排列。通過以上處理就明確了每條剖切線平移時(shí)的次序和位置，但平移的方向卻不能確定。因?yàn)榻?jīng)過旋轉(zhuǎn)后的剖面可能在Z軸的上側(cè)也可能在下側(cè)或者就在Z軸上.因此需要對向上、向下平移或者是不平移進(jìn)行判斷。本文是根據(jù)過原點(diǎn)剖面法向量所在直線與索引值較小井的位置關(guān)系來確定Z坐標(biāo)值符兮，即判斷剖面在Z軸下側(cè)還是非下側(cè)來確定Z坐標(biāo)值的符號。通過Z坐標(biāo)值符號的變化設(shè)置一個(gè)標(biāo)記變量來控制剖切線的平移方向。過原點(diǎn)剖面法向量所在的直線方程為：y-xz+y=0（穌#0H.k2#0）過原點(diǎn)作一個(gè)與剖面平行的Z

21、'軸，然后再作一個(gè)與Z1軸垂直的r軸。其中,0是zy坐標(biāo)系中剖面在r軸的方向上的截距。根據(jù)旋轉(zhuǎn)角度所在象限及,。'與o的大小關(guān)系，分以下八種情況對于z的符號進(jìn)行討論：1）在0<a<出2,為'>0的悄況下，剖面位于軸的下側(cè).此時(shí)2符號為負(fù)；2）在宣/2<a<>0的情況下，剖面位于Z'軸的下側(cè)，此時(shí)z符號為正；3）在0<a<”/2,為'W0的情況下，剖面位于Z'軸的非下側(cè)，此時(shí)z符號為正；4) 在ir/2<a<n,o'W。的情況下，剖面位于軸的非下側(cè).此時(shí)z符號為負(fù)；5) 在-”/2

22、<aWO,%'>0的情況下,剖面位于Z，軸的下側(cè)，此時(shí)z符號為正；6) 在-it<a<-ir/2,y0*>0的情況下，削面位于Z，軸的下側(cè)，此時(shí)z符號為負(fù)；7) 在-”/2<aW0,/。'W0的情況下，剖面位于Z軸的非下側(cè)，此時(shí)z符號為負(fù)；8) 在-ar<a<-竹/2,為'0的情況下，削面位于8軸的非下側(cè)，此時(shí)z符號為正。2利用OpenInventor的實(shí)現(xiàn)以下，描述了利用OpenInventor實(shí)現(xiàn)的連井剖面展開的算法過程。首先,描述了利用研究區(qū)的波阻抗數(shù)據(jù)；然后，展示了在OpenInventor平臺上實(shí)現(xiàn)的連井剖面展開

23、的算法的中間結(jié)果和最終結(jié)果;最后，對于展開的結(jié)果進(jìn)行了分析。所述實(shí)例是利用r一個(gè)油田的三維地震波阻抗數(shù)據(jù)，以及所鉆成的五口井。地震波阻抗數(shù)據(jù)是從0-10小層的頂界，到0-13小層的底界.也就姑0-10+0-11+0-12+03。這些井是CN-42.CN-41,CN-12O,CJS-1-O,CN-5O(圖7)。通過這五口井的那個(gè)制面是由該油田的波阻抗數(shù)據(jù)體生成的。這個(gè)削面利用所切開的剖面，形象地對波阻抗數(shù)據(jù)進(jìn)行了可視化。觀察圖7波阻抗從底面到頂面的變化大致呈現(xiàn)出從大到小的趨勢，這說明在頂部主要呈現(xiàn)為泛濫平原微相，而在其底部則主要呈現(xiàn)為河道、心灘等砂巖微相,油氣主要是集中在底部的砂巖中。地K波阻抗

24、圖7研究區(qū)的波阻抗數(shù)據(jù)2.1三維連井削面的生成由OpenInventor中的SbPlane類來生成三維的連井剖面，調(diào)用該類的構(gòu)造函數(shù)SbPlane()來確定每個(gè)剖面段在空間中的位置。SbPlane()包含兩個(gè)參數(shù)，第一個(gè)參數(shù)是坐標(biāo)原點(diǎn)到剖面段的距離,第二個(gè)參數(shù)是剖面段的法向量。生成三維連井削面的關(guān)鍵代碼如下：SbPlaneplane;Plane=SbPlane(constSbVec3f&normal,floatdistance):三維連井制面顯示的效果如圖8所示，在該圖中顯示的有5口井、4個(gè)剖面段。剖面段排列順序按輸入各口井的順序。第一口井到第二口井為削面段1,剖面段2,，剖面段4的意

25、義以此類推。CN-50圖8三維連井剖面2.2連井副面的旋轉(zhuǎn)由OpenInventor中的SoTransform類實(shí)現(xiàn)各剖面段的旋轉(zhuǎn)。調(diào)用該類的rotation()方法對生成的三維連井剖面圖進(jìn)行旋轉(zhuǎn)。rotation()方法包含兩個(gè)參數(shù).一個(gè)是剖面段旋轉(zhuǎn)所陶繞的旋轉(zhuǎn)法向量，另一個(gè)是旋轉(zhuǎn)的角度。旋轉(zhuǎn)變換的關(guān)鍵代碼如下：transfonnfg_countSlice|->rotation.setValue(SbVec3f(xOfjcountSIice.y0_(g_<*ountSlice,10(g_countSlice),i0g_countSlice);經(jīng)過旋轉(zhuǎn)后的連井剖面如圖9所示。其中，

26、旋轉(zhuǎn)后的各剖面段按圖中標(biāo)定的數(shù)值依次進(jìn)行排列。CJS-1-0CN-50圖9施轉(zhuǎn)后的剖面2.3連井削面的平移展開由OpenInventor中的SoTransIation類實(shí)現(xiàn)各剖面段的平移。調(diào)用該類的translation)方法對旋轉(zhuǎn)后的連井剖面進(jìn)行平移。translation()方法包含一個(gè)參數(shù),是剖面段進(jìn)行平移時(shí)的平移向缺。平移變換的關(guān)鍵代碼如下：transfonng_countSlice->translation,setValue(SbVec3f(xl|g_counlSlice,ylg_countSlice,zlR_counlSlice);經(jīng)過平移后的二維展開剖面如圖10所示。其中，

27、平移后的各剖面段按圖中從左到右順序依次排列，而且垂直的5條線段代表了5口井。CN-42CN-41CN-120CJS-1-0CN-50圖10二維連井削面展開圖2.4實(shí)例分析通過儲層連井剖面圖展開算法，首先獲得如圖8所示的三維空間中的一個(gè)剖面圖。這個(gè)剖面由5口井組成,從左面往右面排列，最后一口井(CN-50)在第4口井的北面。整個(gè)剖面含有一個(gè)明顯的轉(zhuǎn)彎。獲得的如圖10所示的二維剖面圖中,這個(gè)研究區(qū)的上部主要由波阻抗較低的巖性組成，也即由泛濫平原組成,而其下部的波阻抗數(shù)值多為較大，說明由河道和心灘等砂巖組成。這個(gè)認(rèn)識和對圖7的解釋是一致的。不一樣的是，圖7關(guān)于波阻抗數(shù)值的解釋是基于研究區(qū)的側(cè)面的波阻

28、抗數(shù)值作出的，而剖面圖關(guān)于波阻抗的解釋是基于在研究區(qū)內(nèi)部的一個(gè)剖面上的數(shù)值。因此，剖面圖可以顯示三維研究區(qū)內(nèi)部的沉積微相的變化。3結(jié)語本文的連井剖面展開算法實(shí)現(xiàn)了從數(shù)據(jù)到三維圖形，再到二維圖形的轉(zhuǎn)變。它不僅僅實(shí)現(xiàn)了數(shù)-形轉(zhuǎn)換，而且將三維地震數(shù)據(jù)的三維剖面圖轉(zhuǎn)換為一個(gè)二維平面內(nèi)展開的剖面OOpenInventor具有強(qiáng)大的功能，是一種面向?qū)ο蟮目梢暬_發(fā)方法，特別有利于儲層地質(zhì)建?？梢暬到y(tǒng)的開發(fā)和完善。這個(gè)可視化系統(tǒng)的其他子系統(tǒng)的研發(fā)，諸如等值線圖、三維地質(zhì)圖件、柵狀圖等，都可以在OpenInventor平臺上實(shí)現(xiàn)。在本文算法的基礎(chǔ)上，可以再標(biāo)上各井處的測井?dāng)?shù)據(jù)、井名，橫向和縱向的距離標(biāo)尺等技

29、術(shù)數(shù)據(jù)，從而最終生成完整的、實(shí)用性很強(qiáng)的剖面圖。本文算法還適用于利用剖面圖表征孔隙度、滲透率、含油飽和度等物性參數(shù)的空間變化。參考文獻(xiàn)：1 王家華.迎接油氣儲層建模理論、應(yīng)用的大發(fā)展一從2007年國際石油地質(zhì)統(tǒng)計(jì)學(xué)大會談起J.地學(xué)前緣，2008,15(1):16-25.2 王家華.張團(tuán)峰.油弋儲層隨機(jī)建模M.北京：石油工業(yè)出版社,2001.3 王家華.高海余，周葉.克里金地質(zhì)繪圖技術(shù)一計(jì)算機(jī)的模型(上接第273頁)在本次實(shí)驗(yàn)中，經(jīng)過與原版紙張的字體核對，匹配錯(cuò)誤的是146號和117號(如圖5(a)所示)，其匹配度為3.88(匹配度閾值P=3),89號和4號碎片(如圖5(b)所示)，其匹配度為3

30、.031。筆者分析r拼接復(fù)原過程中產(chǎn)生誤差的原因，本文算法只考慮了左右邊緣像素匹配，而沒有考慮拼接后的字體是否成實(shí)存在，雖然這兩個(gè)匹配對完全符合要求且肉眼無法分辨，但這樣的拼接仍然是錯(cuò)誤的。在這種概率極小的情形下，需要人工干預(yù)去修正錯(cuò)誤，但相比整個(gè)拼接復(fù)原過程而言，工作量是很小的。因此，該拼接實(shí)驗(yàn)說明本文的拼接模型具有很大的可行性和操作性。5結(jié)語本文提出了基于碎片的邊緣像素、字體大小、行間距以及段落間距特征的中文文檔碎片自動(dòng)拼接算法。主要介紹了基于匹配度函數(shù)的邊緣像素列橫向拼接方法、基于邊緣像素漸變性的像素列橫向拼接方法以及基于固定行高和行間距特征的碎條縱向拼接方法，并根據(jù)模型編寫了Matla

31、b函數(shù)。基于上述理論，筆者對一實(shí)際中文文檔進(jìn)行了拼接實(shí)驗(yàn),實(shí)驗(yàn)表明本文提出的自動(dòng)拼接算法可行性強(qiáng)，匹配度高，是有效的拼接復(fù)原算法。由于中文與英文文本特征不同，單就本文算法用于拼接英文文檔的匹配度不是很高，但若在本文算法基礎(chǔ)上再結(jié)合另個(gè)針對大小寫英文像素特征的算法，即可得到針對英文文檔拼接的算法C和算法M.北京:石油工業(yè)出版社J999.4 WUR.AprocedurefortheconfigurationofaninflowconlroldevicecompletionusingreservoirmodellingandsimulationintheNorthAmethystPoolC/Proc

32、eedingsoftheSPEReucrvoirCharacterisationandSimulationConferenceandExhibition.S.1.:SocietyofPetruleumEngineers,2011:SPE-147960-MS.5 SENG.SequencestratigraphicmodelingusingOutcropdatain3DspaceC/SPEMiddleEastOilandGasShowandConference.S.i.:SocietyofPetroleumEngineers,2013:SPE-164396-MS.6 潘少偉.楊少春，李胡送.基于

33、開放圖形厚的儲層削面圖的繪制J.中國石油大學(xué)學(xué)報(bào)：自然科學(xué)版，2009.33(1):157-162.7 李華，王明君,傅旭志，等.基于PenInventor的三維地層可視化IJ.物探化探計(jì)算技術(shù),2008,30(4)：345-347.8 劉彩云，陳忠.插值算法在彩色地質(zhì)剖面圖繪制中的應(yīng)用J).石油天然氣學(xué)報(bào)(江漢石油學(xué)院學(xué)報(bào))，2008,30(5):242-243.(9|朱瑩.劉學(xué)宮.陳較忠.基于GIS的地質(zhì)剖面圖自動(dòng)繪制軟件的研究JJ.南京師大學(xué)報(bào)：自然科學(xué)版，2007.30(4):104-108.10 陳錫民.徐文立.基于知識的井間油層連通0動(dòng)對比系統(tǒng)J.清華大學(xué)學(xué)報(bào)：自然科學(xué)版,1999.39(7):89-92.11 王家華.吳少波.王寶清.等.姬源油田耿19-耿20井區(qū)長2=維地質(zhì)建模