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

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

2、三維制面圖生成、就捋、平移。遵楣該算法，通過OpenInventor有關(guān)的類和函數(shù)，將三雄空間中的連井剖面圖展升成為二維空間中的剖面圖。油田數(shù)據(jù)的實驗結(jié)果表明，這個算法具有直觀、效率高等優(yōu)點。關(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基金項目：國家科技放大專項(2011ZX05032Q01)；國家自然科學(xué)驀金資助項目(50874091)。作者簡介:王家華(1945-),界，浙江紹興人，教授，主要研究方向:地質(zhì)統(tǒng)計學(xué)、油氣藏地質(zhì)建模算法及應(yīng)用、計算機圖形學(xué)在油田開發(fā)中的應(yīng)用；休家新(1968-),男，湖北宜昌人，教授，主要研究方向：油氣儲房地質(zhì)遂模、地質(zhì)圖件町視化；覆飛飛(19

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

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

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

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

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

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

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

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

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

15、發(fā)現(xiàn)，當(dāng)&=0或者刈=0時，該公式不成立。因此需要對法向呈的，值根據(jù)站和k2是否等于零分四種情況進行確定。1）在用=0且刈=0的情況下，此時剖切線為一個點，剖切線不存在，因此法向最也不存在，說明兩口井在軸上重合。2）在如=0,虹，。的情況下：a）如果氣=0,則剖面的法向域為”=（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)變換和平移變換兩個步驛。1.2.1族捋變換在三維連井剖面中，前一節(jié)雖然確定了各條剖切線的位置，但剖切線也是有方向的。在進行選井剖切時，將所選的第一口井的索引值置為0,第二口井的索引值為1,第三口井的索引值為2,，第/V口井的索引值為N-1。剖切線的正方向規(guī)定為小索引值的井到大索引值的井所指的方向。在圖3中，這個正方向即為點X（z0,y0）到點8（勾，為）所指的方向。各剖切線的角度為剖切線正方向與Z軸正方向所成的夾角。如圖4所示的角a即為剖切線的角度。通過反正切歐數(shù)atan2（doublek、,

17、double站）來計算削切線的角度a,其中4=力-yo，M=玉-初。在本文中規(guī)定始終以K軸正方向為右手旋轉(zhuǎn)的法向量，所有的剖切線都繞著該法向最旋轉(zhuǎn)a角度。經(jīng)過旋轉(zhuǎn)后，各個剖切線都與Z軸正方向平行且同向。圖4中的剖切線在經(jīng)過旋轉(zhuǎn)a角度后（如圖4為逆時針旋轉(zhuǎn)a角度），生成一條新的剖切線°圖5所示的點4（z。,為）到點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）計算旋轉(zhuǎn)角度，此函數(shù)的返回值為（-&q

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

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

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

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

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

23、的算法的中間結(jié)果和最終結(jié)果;最后，對于展開的結(jié)果進行了分析。所述實例是利用r一個油田的三維地震波阻抗數(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)。通過這五口井的那個制面是由該油田的波阻抗數(shù)據(jù)體生成的。這個削面利用所切開的剖面，形象地對波阻抗數(shù)據(jù)進行了可視化。觀察圖7波阻抗從底面到頂面的變化大致呈現(xiàn)出從大到小的趨勢，這說明在頂部主要呈現(xiàn)為泛濫平原微相，而在其底部則主要呈現(xiàn)為河道、心灘等砂巖微相,油氣主要是集中在底部的砂巖中。地K波阻抗

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

25、義以此類推。CN-50圖8三維連井剖面2.2連井副面的旋轉(zhuǎn)由OpenInventor中的SoTransform類實現(xiàn)各剖面段的旋轉(zhuǎn)。調(diào)用該類的rotation()方法對生成的三維連井剖面圖進行旋轉(zhuǎn)。rotation()方法包含兩個參數(shù).一個是剖面段旋轉(zhuǎn)所陶繞的旋轉(zhuǎn)法向量，另一個是旋轉(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ù)值依次進行排列。CJS-1-0CN-50圖9施轉(zhuǎn)后的剖面2.3連井削面的平移展開由OpenInventor中的SoTransIation類實現(xiàn)各剖面段的平移。調(diào)用該類的translation)方法對旋轉(zhuǎn)后的連井剖面進行平移。translation()方法包含一個參數(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實例分析通過儲層連井剖面圖展開算法，首先獲得如圖8所示的三維空間中的一個剖面圖。這個剖面由5口井組成,從左面往右面排列，最后一口井(CN-50)在第4口井的北面。整個剖面含有一個明顯的轉(zhuǎn)彎。獲得的如圖10所示的二維剖面圖中,這個研究區(qū)的上部主要由波阻抗較低的巖性組成，也即由泛濫平原組成,而其下部的波阻抗數(shù)值多為較大，說明由河道和心灘等砂巖組成。這個認(rèn)識和對圖7的解釋是一致的。不一樣的是，圖7關(guān)于波阻抗數(shù)值的解釋是基于研究區(qū)的側(cè)面的波阻

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

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

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

31、b函數(shù)?；谏鲜隼碚?，筆者對一實際中文文檔進行了拼接實驗,實驗表明本文提出的自動拼接算法可行性強，匹配度高，是有效的拼接復(fù)原算法。由于中文與英文文本特征不同，單就本文算法用于拼接英文文檔的匹配度不是很高，但若在本文算法基礎(chǔ)上再結(jié)合另個針對大小寫英文像素特征的算法，即可得到針對英文文檔拼接的算法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é)報：自然科學(xué)版，2009.33(1):157-162.7 李華，王明君,傅旭志，等.基于PenInventor的三維地層可視化IJ.物探化探計算技術(shù),2008,30(4)：345-347.8 劉彩云，陳忠.插值算法在彩色地質(zhì)剖面圖繪制中的應(yīng)用J).石油天然氣學(xué)報(江漢石油學(xué)院學(xué)報)，2008,30(5):242-243.(9|朱瑩.劉學(xué)宮.陳較忠.基于GIS的地質(zhì)剖面圖自動繪制軟件的研究JJ.南京師大學(xué)報：自然科學(xué)版，2007.30(4):104-108.10 陳錫民.徐文立.基于知識的井間油層連通0動對比系統(tǒng)J.清華大學(xué)學(xué)報：自然科學(xué)版,1999.39(7):89-92.11 王家華.吳少波.王寶清.等.姬源油田耿19-耿20井區(qū)長2=維地質(zhì)建模