專利名稱：一種基于變形場(chǎng)的地層面計(jì)算機(jī)圖形生成方法
技術(shù)領(lǐng)域：
本發(fā)明涉及一種地層面模型的生成方法，特別是涉及計(jì)算機(jī)領(lǐng)域的一種基于變形場(chǎng)的地層面計(jì)算機(jī)圖形模型生成方法。
背景技術(shù)：
野外勘察和量測(cè)獲得的地質(zhì)信息是一些離散的數(shù)據(jù)，面對(duì)這些地質(zhì)信息，地質(zhì)工程師必然提出這樣的問題，如何運(yùn)用這些資料來推斷其在研究區(qū)內(nèi)的發(fā)育分布規(guī)律。即使能夠預(yù)測(cè)地質(zhì)信息在區(qū)域的分布情況，面對(duì)大量的輸出數(shù)據(jù)，地質(zhì)人員仍然會(huì)難于想象和分析。圖形、圖象是最直觀的數(shù)據(jù)解釋，且地質(zhì)人員又習(xí)慣于用圖件來反映地質(zhì)信息，自然會(huì)希望利用計(jì)算機(jī)來顯示結(jié)構(gòu)面等在地質(zhì)體中的分布情況。
圖1所示為一正逆斷層示意圖，圖2所示為一逆斷層俯視圖，圖3所示為一斷層相交時(shí)的示意圖，圖4所示為存在兩個(gè)相交斷層的地層形態(tài)圖，圖5所示為一斷層示意圖。目前，生成地層主要存在以下困難斷層數(shù)量較多，斷層之間關(guān)系復(fù)雜，存在逆斷層；斷層有相交的情況；一個(gè)地質(zhì)模型可能包括很多斷層；存在逆斷層時(shí)，根據(jù)地層被撕裂的情形，將斷層分為兩種類型正斷層和逆斷層。在地層撕裂處以垂直方向看去，若被撕裂的地層兩邊沒有重疊，這種斷層稱為正斷層，否則稱為逆斷層，如圖1所示。在正斷層情況下，建立地層的難度相對(duì)較小，在逆斷層的情況下生成的地層有重疊的部分，如圖2所示，很多軟件在處理過程中采用了正斷層和逆斷層分開處理的方式，正斷層較少的認(rèn)為控制，逆斷層加入較多的人為控制因素。當(dāng)斷層有相交的情況發(fā)生時(shí)，生成地層也有困難，兩個(gè)或兩個(gè)以上的斷層相交，相交處附近的地層受兩個(gè)斷層的影響，兩個(gè)斷面的交線處的地層要完全和斷層走向貼和，如圖3所示，其中，30為斷面，31為兩斷面的相交處；斷面相交時(shí)生成的地層如圖4所示，其中40表示兩斷層相交處地層也在兩斷面上，41表示地層在斷層處的斷裂。圖5所示為某油田數(shù)據(jù)生成的斷層示意圖，說明了一個(gè)地質(zhì)模型可能包括很多斷層。在進(jìn)一步根據(jù)地層求解地層相關(guān)的物性參數(shù)，小層和地層剖面圖的過程中也會(huì)存在地層生成困難，生成地層只是建模的第一步工作，在地層生成后還有一系列建模工作要做，比如油田部門要計(jì)算儲(chǔ)油情況，要生成各地層的縱向剖切圖等工作，這些工作如果要通過手工或者計(jì)算交線的方法可能會(huì)造成不準(zhǔn)確和耗費(fèi)時(shí)間。
發(fā)明內(nèi)容
本發(fā)明所要解決的技術(shù)問題在于提供一種基于變形場(chǎng)的地層面計(jì)算機(jī)圖形生成方法，用以解決現(xiàn)有技術(shù)中用戶大量通過手動(dòng)操作來生成地層面，對(duì)斷層較多情況無法自動(dòng)生成地層面，從而導(dǎo)致精確度不高和效率較低的問題。
為了實(shí)現(xiàn)上述目的，本發(fā)明提供了一種基于變形場(chǎng)的地層面計(jì)算機(jī)圖形生成方法，用以通過地震數(shù)據(jù)和測(cè)井解釋結(jié)果數(shù)據(jù)進(jìn)行計(jì)算機(jī)處理而自動(dòng)生成地層面模型，其特征在于，包括步驟一，根據(jù)一地層斷面的空間位置建立一斷面樹；步驟二，根據(jù)所述斷面樹劃分地層離散點(diǎn)；步驟三，由所述地層離散點(diǎn)建立一變形場(chǎng)；步驟四，建立連續(xù)的地層曲面；及步驟五，通過所述變形場(chǎng)的作用力將所述地層曲面變形為實(shí)際地層面。
上述的基于變形場(chǎng)的地層面計(jì)算機(jī)圖形生成方法，其中，所述方法進(jìn)一步包括步驟A，根據(jù)地層斷面空間位置建立的斷面樹決定斷面斷裂的先后順序和各斷面之間的影響關(guān)系；步驟B，由所述斷面樹所表示的空間位置關(guān)系將地層的離散點(diǎn)分布到所述斷面樹的各個(gè)葉子節(jié)點(diǎn)；步驟C，通過所述斷面樹和斷面樹中分布的地層離散點(diǎn)逐個(gè)建立斷面變形場(chǎng)，并通過該變形場(chǎng)的逆作用將全部該離散點(diǎn)恢復(fù)為一連續(xù)地層面的離散點(diǎn)；步驟D，由全部恢復(fù)的離散點(diǎn)來建立一連續(xù)的地層曲面；及步驟E，所述變形場(chǎng)按所述先后順序作用于所述地層曲面，將該地層曲面撕裂變形為實(shí)際地層面。
上述的基于變形場(chǎng)的地層面計(jì)算機(jī)圖形生成方法，其中，所述步驟A中還包括一建立地層斷面的步驟，該地層斷面由所輸入數(shù)據(jù)中斷面的離散點(diǎn)集合采用計(jì)算機(jī)自由曲面生成算法來建立。
上述的基于變形場(chǎng)的地層面計(jì)算機(jī)圖形生成方法，其中，所述輸入數(shù)據(jù)包括地層面的離散點(diǎn)集合、斷面的離散點(diǎn)集合。
上述的基于變形場(chǎng)的地層面計(jì)算機(jī)圖形生成方法，其中，所述步驟C中還包括一為建立斷面變形場(chǎng)步驟，包括步驟C1，根據(jù)斷面兩側(cè)的地層離散點(diǎn)插值出兩片小的地層曲面；步驟C2，由所述兩個(gè)地層曲面分別求出兩條交線；及步驟C3，通過所述兩條交線的Z值差，插值出斷層面上Z方向的變形量。
上述的基于變形場(chǎng)的地層面計(jì)算機(jī)圖形生成方法，其中，所述先后順序?yàn)閿鄬用娣e大的斷面先斷裂，斷層面積小的斷面后斷裂。
上述的基于變形場(chǎng)的地層面計(jì)算機(jī)圖形生成方法，其中，所述步驟C、E中，所述變形場(chǎng)的建立和離散點(diǎn)的恢復(fù)過程按所述先后順序的逆序進(jìn)行；所述實(shí)際地層面的生成過程按所述順序進(jìn)行。
上述的基于變形場(chǎng)的地層面計(jì)算機(jī)圖形生成方法，其中，所述斷面變形場(chǎng)是根據(jù)斷層兩側(cè)離散點(diǎn)的落差和斷面的形態(tài)來建立的。
上述的基于變形場(chǎng)的地層面計(jì)算機(jī)圖形生成方法，其中，所述步驟E中，所述實(shí)際地層面為一由計(jì)算機(jī)三角網(wǎng)格或自由曲面表示的地層面。
上述的基于變形場(chǎng)的地層面計(jì)算機(jī)圖形生成方法，其中，所述地層面為一不連續(xù)曲面。
上述的基于變形場(chǎng)的地層面計(jì)算機(jī)圖形生成方法，其中，所述斷面樹為一二叉樹。
上述的基于變形場(chǎng)的地層面計(jì)算機(jī)圖形生成方法，其中，所述二叉樹的節(jié)點(diǎn)包括節(jié)點(diǎn)主斷面和子斷面列表。
上述的基于變形場(chǎng)的地層面計(jì)算機(jī)圖形生成方法，其中，所述節(jié)點(diǎn)的左子節(jié)點(diǎn)及該左子節(jié)點(diǎn)的子節(jié)點(diǎn)在所述節(jié)點(diǎn)主斷面的左側(cè)；所述節(jié)點(diǎn)的右子節(jié)點(diǎn)及該右子節(jié)點(diǎn)的子節(jié)點(diǎn)在所述節(jié)點(diǎn)主斷面的右側(cè)。
以下結(jié)合附圖和具體實(shí)施例對(duì)本發(fā)明進(jìn)行詳細(xì)描述，但不作為對(duì)本發(fā)明的限定。
圖1為一正逆斷層示意圖；圖2為一逆斷層俯視圖；圖3為一斷層相交時(shí)的示意圖；圖4為存在兩個(gè)相交斷層的地層形態(tài)圖；
圖5為一斷層示意圖；圖6為現(xiàn)有的生成三維地層的流程圖；圖7為基于斷面區(qū)塊方式的自動(dòng)生成地層的流程圖；圖8為本發(fā)明的地層面的離散點(diǎn)與斷面的示意圖；圖9為本發(fā)明的地層面生成的總流程圖；圖10為一地層面的斷層俯視示意圖；圖11為一斷面關(guān)系樹；圖12為本發(fā)明的一斷面示意圖；圖13～圖20為本發(fā)明實(shí)施例的一離散點(diǎn)生成地層的過程；及圖21～圖36為本發(fā)明實(shí)施例的一地層生成具體過程示意圖。
其中，附圖標(biāo)記30 斷面31 兩斷面的相交處40 兩斷層相交處地層也在兩斷面上41 地層在斷層處的斷裂70、71、72 斷面73 地層的離散點(diǎn)
具體實(shí)施方式
圖6所示為現(xiàn)有的生成三維地層的流程圖。其在建立三維地層面模型時(shí)需要大量的手工編輯操作工作，在通過地層計(jì)算地層的物性參數(shù)和儲(chǔ)層模型也需要大量的手動(dòng)編輯計(jì)算，使用本方法建立三維地質(zhì)模型—包括地層，地層的物性參數(shù)和儲(chǔ)層。
對(duì)于所述流程圖采用的方法，地質(zhì)工作人員只要負(fù)責(zé)提供原始的建模離散點(diǎn)就可以建立斷裂的地層模型并進(jìn)行相關(guān)的進(jìn)一步建模。具體步驟包括步驟600，首先建立斷面模型；步驟610，根據(jù)斷面建立連續(xù)地層面模型；步驟620，根據(jù)地層和斷層的交線編輯斷裂地層；及步驟630，判斷地層是否符合要求，若符合要求，則結(jié)束過程，若不符合要求，則執(zhí)行步驟620。
地層面又稱層面，地層面是地層與地層的分界面，一般來說同一地層內(nèi)的密度是不變的，但地層與地層之間通常具有不同的地質(zhì)屬性，這是因?yàn)樗鼈兊牡刭|(zhì)構(gòu)造不同；一般來說地層面就是地震射線反射面。一個(gè)地層面可以唯一的代表其下面緊鄰的地層，因此，在地質(zhì)處理系統(tǒng)中通常將對(duì)地層的研究轉(zhuǎn)化為對(duì)地層面的研究。
每個(gè)地層都具有不同的地質(zhì)屬性，如地層的密度、地震波在地層中的傳播速度等；這些屬性是地質(zhì)解釋中不可缺少的數(shù)據(jù)來源。通常，某個(gè)屬性在地層內(nèi)是處處可變的，但為了簡(jiǎn)化模型，一般認(rèn)為同一個(gè)地層內(nèi)的屬性值是常數(shù)。
斷層面又稱斷面，是因受擠壓或拉伸等作用過度，使地層發(fā)生斷裂而形成的面。
地質(zhì)面包括地層面與斷層面，地質(zhì)面與地質(zhì)面之間可能相交。斷層面是由于拉伸等作用過度使地層發(fā)生扭曲、變形直至出現(xiàn)裂口而形成的面，因此斷層面是與部分地層面相交的；地下地質(zhì)構(gòu)造可能相當(dāng)復(fù)雜。斷層描述地層的斷裂情況，若斷層與某個(gè)地層A相交，則表明在相交處A會(huì)出現(xiàn)裂口。斷層與斷層也可以相交。
圖7所示為本發(fā)明的地層面的離散點(diǎn)與斷面的示意圖，圖8所示為本發(fā)明的地層面生成的總流程圖。地層在某一年代因地殼運(yùn)動(dòng)擠壓或拉伸產(chǎn)生斷裂，在地殼上的力不僅使裂口處地層發(fā)生形變也帶動(dòng)整個(gè)地層發(fā)生拉伸或擠壓形變，經(jīng)過地層產(chǎn)生形變，經(jīng)過多次形變地層達(dá)到目前的形態(tài)，地質(zhì)工作者采集到地層上的離散點(diǎn)數(shù)據(jù)。本發(fā)明的工作是根據(jù)離散點(diǎn)利用計(jì)算機(jī)來建立地層曲面。直接根據(jù)地層離散點(diǎn)來生成斷裂的地層要考慮很多斷層的信息，問題復(fù)雜。為了降低地層生成的難度，本方法將地層是逐漸斷裂的情況作為已知條件來看待，地層的每次發(fā)生斷裂形變就是受到力的作用發(fā)生形變，可以用一個(gè)變形力場(chǎng)來表示。為每次斷裂建立一個(gè)變形場(chǎng)，將變形場(chǎng)按順序逐個(gè)作用于一個(gè)連續(xù)的地層面就可以建立一個(gè)斷裂的地層面。每個(gè)斷層描述地層的一個(gè)裂口，所以方法為每個(gè)斷層建立一個(gè)變形場(chǎng)，地層的離散點(diǎn)通過變形場(chǎng)的逆作用恢復(fù)為連續(xù)地層面的離散點(diǎn)，根據(jù)連續(xù)面的離散點(diǎn)生成曲連續(xù)的地層曲面，連續(xù)的地層曲面就可以通過變形場(chǎng)作用斷裂為不連續(xù)的曲面，本地層生成方法的優(yōu)點(diǎn)是無須求曲面之間交線，可以處理裂口及斷面形態(tài)復(fù)雜的情況，數(shù)值穩(wěn)定性更高。有利于進(jìn)一步的三維模型處理及計(jì)算，如生成剖面圖，小層等。本方法的流程包括方法輸入數(shù)據(jù)和計(jì)算結(jié)果輸出。
生成地層面計(jì)算機(jī)圖形包括的輸入數(shù)據(jù)主要有地震數(shù)據(jù)和測(cè)井解釋結(jié)果兩個(gè)部分地震解釋結(jié)果數(shù)據(jù)包括地層離散點(diǎn)數(shù)據(jù)、地層上的斷層線數(shù)據(jù)、斷面離散點(diǎn)數(shù)據(jù)。這些數(shù)據(jù)可以是地震解釋軟件生成的文件，也可以是從構(gòu)造圖掃描后采集得到的數(shù)據(jù)。測(cè)井解釋結(jié)果數(shù)據(jù)包括層位信息表、井分層信息表、斷面表、井?dāng)帱c(diǎn)表、井斜信息表，系統(tǒng)也可直接從用戶的數(shù)據(jù)庫中讀取上述信息。
在本發(fā)明中，輸入數(shù)據(jù)包括地層面的離散點(diǎn)集合ScatterSet，描述地層裂口的斷面集合FaultSet，斷層面很多且互相相交截?cái)?，關(guān)系復(fù)雜，如圖7所示，其中，70、71、72皆為斷面，73為地層的離散點(diǎn)。
計(jì)算結(jié)果是以任何形式表示的地層面及相關(guān)的儲(chǔ)層面；1、自由曲面表示法(參見，S.Auerbach，H.Schaeben，Computer aided geometric design of geologicsurfaces and bodies[J]，Journal of Math.Geology，Vol.22，No.8，1990.8；或者，張劍秋，張福炎，地震層位信息三維可視化[J]，石油地球物理勘探，Vol.33，No.1，1998.2，P117～124)；2、三角網(wǎng)格(TIN)曲面表示法(參見，楊欽，限定Delaunay三角剖分，北京航空航天大學(xué)博士論文，2001.8；或者，牛文杰，三維數(shù)據(jù)場(chǎng)可視化克里金建模及其算法的理論和應(yīng)用研究，北京航空航天大學(xué)博士論文，2002.9)；3、PEBI網(wǎng)格曲面表示法(參見，Z.E.Heinemann，C.W.Brand，Griddingtechniques in reservoir simulation[C]，presented at First and Second InternationalForum on Reservoir Simulation，1989，P339～426；或者，謝海兵，馬遠(yuǎn)樂，恒冠仁，郭尚平，非結(jié)構(gòu)網(wǎng)格油藏?cái)?shù)值模擬方法研究[J]，石油學(xué)報(bào)，Vol.22，No.1，2001.1，P63～67)。在本發(fā)明中，計(jì)算結(jié)果輸出是由三角網(wǎng)格表示的地層面，即不連續(xù)曲面LayerFace。該結(jié)果也為以自由曲面來表示的地層面。該地層面在斷層處出現(xiàn)裂口并且裂口邊緣貼在斷面上，地層上的已知離散點(diǎn)在生成的地層面上。
根據(jù)以上的思路，地層生成方法的內(nèi)容如下根據(jù)斷面的空間位置來建立斷面樹；建立好的斷面樹是一棵二叉樹，通過這個(gè)樹決定地層裂口斷裂的先后順序及之間的相互影響順序。通過這棵樹可以查找空間的拓?fù)潢P(guān)系；根據(jù)斷面樹表示的位置關(guān)系將地層的離散點(diǎn)分布到斷面樹的葉子節(jié)點(diǎn)上；根據(jù)斷面樹和斷面樹中分布好的地層離散點(diǎn)來逐個(gè)建立斷面變形場(chǎng)。到根節(jié)點(diǎn)時(shí)斷層全部恢復(fù)，恢復(fù)好的地層離散點(diǎn)可以創(chuàng)建一個(gè)連續(xù)地層面；用恢復(fù)好的離散點(diǎn)建立連續(xù)的地層面；通過變形場(chǎng)的作用力將連續(xù)的地層面變形為撕裂的地層面。
圖8所示的該總流程圖包括步驟800，首先根據(jù)斷面建立斷面樹；步驟810，根據(jù)斷面樹劃分地層離散點(diǎn)；步驟820，根據(jù)地層離散點(diǎn)建立變形場(chǎng)；步驟830，建立連續(xù)的地層曲面；及步驟840，根據(jù)變形場(chǎng)將連續(xù)的地層面變形為實(shí)際地層面。
斷面樹的建立。作用于地層面的變形場(chǎng)可以通過裂口兩側(cè)離散點(diǎn)的落差和斷層的角度模擬出來。在斷層兩側(cè)地層離散點(diǎn)有明顯的落差對(duì)于只有一個(gè)的裂口地層的變形場(chǎng)建立起來相對(duì)容易，只要考慮單個(gè)斷面的形狀和傾斜角度。但當(dāng)裂口的數(shù)目增加時(shí)，為每個(gè)描述裂口的斷面建立一個(gè)變形場(chǎng)時(shí)要考慮其他裂口相關(guān)，每個(gè)裂口的斷裂時(shí)間不同，導(dǎo)致變形場(chǎng)的形態(tài)，影響范圍也不同，如何為多個(gè)互相有干擾的斷面建立變形場(chǎng)是本算法解決問題的基礎(chǔ)。
裂口在不同年代產(chǎn)生的，早發(fā)生的地層裂口的形態(tài)會(huì)在下一次斷裂中發(fā)生進(jìn)一步的形變，比較理想的狀況是，得知地層斷裂的順序，逐個(gè)按年代逆序恢復(fù)裂口。本算法設(shè)計(jì)為斷面建立變形場(chǎng)，最終恢復(fù)到地層沒有斷裂的情況，生成連續(xù)曲面，再根據(jù)時(shí)代的順序斷裂連續(xù)曲面生成不連續(xù)的曲面了。生成連續(xù)曲面算法目前已經(jīng)有較成熟的方法。在實(shí)際情況中無法得知斷裂的準(zhǔn)確順序，所以本算法為斷面的斷裂順序制定了如下規(guī)則(1)斷面兩側(cè)有較大的落差要后斷裂；(2)分布在斷面兩側(cè)的斷面之間互相沒有影響，恢復(fù)順序可以不考慮。該利用規(guī)則可使得斷面之間的影響減小。
為此，建立一個(gè)二叉樹來排列斷面斷裂的順序和各個(gè)斷面之間的影響關(guān)系。每個(gè)二叉樹節(jié)點(diǎn)由兩部分信息組成主變形場(chǎng)斷面和子斷面列表。斷面樹的父節(jié)點(diǎn)先于子斷面斷裂，影響子斷面的斷裂形態(tài)。斷面的面積決定斷面的斷裂順序，斷面地理位置信息決定斷面的影響關(guān)系。每個(gè)斷面作為主變形場(chǎng)斷面并不唯一出現(xiàn)，斷面的變形場(chǎng)由一個(gè)或多個(gè)樹節(jié)點(diǎn)決定。
本發(fā)明在設(shè)計(jì)過程中考慮到地層是分年代斷裂的，所以根據(jù)斷面的斷裂順序建立變形場(chǎng)是最接近實(shí)際地理情況下，但是在實(shí)際的數(shù)據(jù)中斷面的斷裂年代順序難以確定，所以本發(fā)明根據(jù)斷面的自身特點(diǎn)設(shè)計(jì)了斷面的斷裂順序——按斷層面積大小斷裂，面積大的認(rèn)為斷裂先發(fā)生，面積小的則后發(fā)生，因?yàn)橥ǔＧ闆r下，斷面描述的裂口通常具有斷面面積較大的地層的落差也比較大，對(duì)其他斷面的影響也比較大，在恢復(fù)斷面的過程中放在其他斷面之后恢復(fù)能減少地層數(shù)據(jù)的波動(dòng)，減小對(duì)變形場(chǎng)建立原則的依賴性更尊重原始數(shù)據(jù)。
因?yàn)橄葦嗔训牧芽谝彩芎髷嗔训男巫兊挠绊?，在某些情況下，因?yàn)閿嗔演^小影響的也較小或因?yàn)榭缌芽跀嗔阉杂绊懸草^小，所以在建立變形場(chǎng)時(shí)要考慮斷面之間的影響關(guān)系。
在變形過程中，變形場(chǎng)的建立和恢復(fù)過程是按斷面面積逆序，斷裂生成地層的順序是斷裂的正序，為了簡(jiǎn)化生成變形場(chǎng)時(shí)許多斷面之間的關(guān)系，為此本發(fā)明建立一棵斷面二叉樹來描述斷面之間的關(guān)系和變形的先后順序，劃分地層空間。斷面樹節(jié)點(diǎn)的結(jié)構(gòu)如下表1-1
每個(gè)斷層面將建?？臻g分為兩個(gè)部分，一部分為稱為左側(cè)空間，另一部分稱為右側(cè)空間，該空間可任意確定，只要保證在整個(gè)計(jì)算過程中，每個(gè)斷面的左側(cè)和右側(cè)唯一確定即可。樹節(jié)點(diǎn)的主斷面MainFault是當(dāng)前斷面樹節(jié)點(diǎn)要描述的斷面，節(jié)點(diǎn)主斷面MainFault將存在的子斷面節(jié)點(diǎn)集合FaultSet劃歸到兩個(gè)子集合，一個(gè)集合的所有斷面都和斷面的左側(cè)空間相干涉，另一個(gè)斷面都在斷面的右側(cè)空間相干涉。
斷面樹描述的斷面關(guān)系如下(1)，斷面樹的父節(jié)點(diǎn)要晚于子節(jié)點(diǎn)建立變形場(chǎng)，子節(jié)點(diǎn)建立變形場(chǎng)時(shí)要受到所有父節(jié)點(diǎn)的影響。(2)，斷面樹的不同分支互相之間不受影響，斷面樹節(jié)點(diǎn)受所有父節(jié)點(diǎn)的影響，此父節(jié)點(diǎn)包括節(jié)點(diǎn)的父節(jié)點(diǎn)及其父節(jié)點(diǎn)的所有父節(jié)點(diǎn)。(3)，斷面的主斷面節(jié)點(diǎn)可能是同一斷面，但同一斷面作為主斷面的樹節(jié)點(diǎn)之間不存在直接或間接互為父節(jié)點(diǎn)的情況。
斷面樹的建立是一個(gè)遞歸的過程，通過逐級(jí)判定斷面的關(guān)系來生成斷面樹。每個(gè)樹的節(jié)點(diǎn)主斷面在是空間劃分?jǐn)嗝?，其左子?jié)點(diǎn)及其子節(jié)點(diǎn)的節(jié)點(diǎn)主斷面在當(dāng)前節(jié)點(diǎn)主斷面的左側(cè)，其右子節(jié)點(diǎn)在當(dāng)前節(jié)點(diǎn)主斷面的右側(cè)。
對(duì)輸入斷面集Faultset的斷面按面積由大到小進(jìn)行排序sort(Faultset)，得到有序斷面集合Faultset’＝{f1，f2，….fn}，對(duì)于這個(gè)集合建立一個(gè)樹節(jié)點(diǎn)FaultTreeNode，F(xiàn)aultTreeNode(MainFault)＝f1，F(xiàn)aultTreeNode(FaultSet)＝{f2，…fn}；根據(jù)GetLeftFaultSet(fault，S)計(jì)算出該子斷面集合中所有與斷面左側(cè)(及其延長(zhǎng)線上)的空間區(qū)域相干予的斷面集合FaultSetLeft，并根據(jù)這個(gè)斷面集合建立子節(jié)點(diǎn)，樹節(jié)點(diǎn)的右子節(jié)點(diǎn)也同理，直到樹節(jié)點(diǎn)FaultTreeNode的FaultSet集合為空，根據(jù)此法則遞歸建立一個(gè)斷面樹描述斷面的地理位置。每個(gè)斷面作為Mainfault一次或多次。此方法的程序示意描述如下。
<pre listing-type="program-listing">FaultTreeNode *creatFaultTree(FaultSet){FaultTreeNode *pRetTreeNode=createATreeNodebyFaultSet(FaultSet)；//根據(jù)斷面集建立一個(gè)斷面樹節(jié)點(diǎn)FaultSetLeft=GetLeftFaultSet(pRetTreeNode-＞Fault，F(xiàn)aultSet)；FaultSetRight=GetLeftFaultSet(pRetTreeNode-＞Fault，F(xiàn)aultSet)；If(IsNotEmpt(FaultSetLeft))LeftTreeNode=creatFaultTree(FaultSet-pRetTreeNode->Fault，)if IsNotEmpt(FaultSetLeft))RightTreenode=creatFaultTree(FaultSet-_pRetTreeNode->Fault，)Return pRetTreeNode}</pre>圖9所示為一地層面的斷層俯視示意圖，圖10所示為一斷面關(guān)系樹。根據(jù)如上方法建立如下關(guān)系樹這樣建立的斷面樹有如下的邏輯關(guān)系，節(jié)點(diǎn)的左子節(jié)點(diǎn)及其子節(jié)點(diǎn)干預(yù)父斷面節(jié)點(diǎn)的空間位置左側(cè)，右子節(jié)點(diǎn)及其子節(jié)點(diǎn)干預(yù)父斷面空間位置右側(cè)，本算法以斷面樹為基礎(chǔ)建立變形場(chǎng)。將ScatterSet中的點(diǎn)按空間位置分配到建好的樹的葉子節(jié)點(diǎn)。即空間位置在斷面左側(cè)的點(diǎn)分配到斷面左側(cè)的子節(jié)點(diǎn)，空間位置在斷面右側(cè)的點(diǎn)分配到斷面右側(cè)的子節(jié)點(diǎn)，直至分配到葉節(jié)點(diǎn)。具體計(jì)算干預(yù)主斷面區(qū)域的方法如上所述。如圖9為斷面俯視示意圖；圖9中所示的斷面位置關(guān)系根據(jù)斷面樹的生成方法生成斷面樹如圖所示。根據(jù)建立好的斷面樹，參見圖10，將地層上的離散點(diǎn)分配到斷面樹的葉子節(jié)點(diǎn)上，因?yàn)殡x散點(diǎn)的影響范圍很小，所以只要根據(jù)空間位置判定即可，空間位置在斷面左側(cè)的點(diǎn)分配到樹節(jié)點(diǎn)主斷面左側(cè)的子節(jié)點(diǎn)，空間位置在樹節(jié)點(diǎn)主斷面右側(cè)的點(diǎn)分配到斷面右側(cè)的子節(jié)點(diǎn)，逐級(jí)按直至分配到葉節(jié)點(diǎn)。如圖9為點(diǎn)p1，p2，p3，p4，p5在斷面樹葉子節(jié)點(diǎn)的位置。
下面給出一地層斷面的斷面二叉樹生成的一實(shí)施例。參見圖5所示，其中，地層信息包括地層離散點(diǎn)的數(shù)量為2337個(gè)，斷層的數(shù)目為15，其斷面是由TIN曲面表示；由該斷面生成的斷面二叉樹為斷面f13的左子樹是f17，斷面f13的右子樹是f17，斷面f17的左子樹是f6，斷面f17的右子樹是f6，斷面f6的左子樹是f9，斷面f6的右子樹是f8，斷面f9的左子樹是f7，斷面f9的右子樹是f7，斷面f7的左子樹是f4，斷面f7的右子樹是f4，斷面f4的左子樹是f18，斷面f4的右子樹是f18，斷面f7的左子樹是f4，斷面f7的右子樹是f4，斷面f4的左子樹是f18，斷面f4的右子樹是f18，斷面f8的左子樹是f10，斷面f8的右子樹是f10，斷面f10的左子樹是f11，斷面f10的右子樹是f11，斷面f11的左子樹是f7，斷面f11的右子樹是f7，斷面f7的左子樹是f4，斷面f7的右子樹是f4，斷面f4的左子樹是f18，斷面f4的右子樹是f18，斷面f18的右子樹是f19，斷面f4的左子樹是f18，斷面f4的右子樹是f18，斷面f7的左子樹是f4，斷面f7的右子樹是f4，斷面f11的左子樹是f7，斷面f11的右子樹是f7，斷面f7的左子樹是f4，斷面f7的右子樹是f4，斷面f4的左子樹是f18，斷面f4的右子樹是f18，斷面f4的左子樹是f18，斷面f4的右子樹是f18，斷面f7的左子樹是f4，斷面f7的右子樹是f4，斷面f10的左子樹是f7，斷面f10的右子樹是f7，斷面f7的左子樹是f4，斷面f7的右子樹是f4，斷面f4的左子樹是f18，斷面f4的右子樹是f18，斷面f4的左子樹是f18，斷面f4的右子樹是f18，斷面f7的左子樹是f4，斷面f7的右子樹是f4，斷面f4的左子樹是f18，斷面f4的右子樹是f18，斷面f4的左子樹是f18，斷面f4的右子樹是f18，斷面f6的左子樹是f9，斷面f6的右子樹是f8，斷面f9的左子樹是f7，斷面f9的右子樹是f7，斷面f7的左子樹是f4，斷面f7的右子樹是f4，斷面f4的左子樹是f18，斷面f4的右子樹是f18，斷面f4的左子樹是f18，斷面f4的右子樹是f18，斷面f7的左子樹是f4，斷面f7的右子樹是f4，斷面f4的左子樹是f18，斷面f4的右子樹是f18，斷面f4的左子樹是f18，斷面f4的右子樹是f5，斷面f18的右子樹是f22，斷面f22的左子樹是f23，斷面f5的右子樹是f18，斷面f18的右子樹是f21，斷面f21的右子樹是f22，斷面f8的左子樹是f10，斷面f8的右子樹是f10，斷面f10的左子樹是f11，斷面f10的右子樹是f11，斷面f11的左子樹是f7，斷面f11的右子樹是f7，斷面f7的左子樹是f4，斷面f7的右子樹是f4，斷面f4的左子樹是f18，斷面f4的右子樹是f18，斷面f4的左子樹是f18，斷面f4的右子樹是f18，斷面f7的左子樹是f4，斷面f7的右子樹是f4，斷面f11的左子樹是f7，斷面f11的右子樹是f7，斷面f7的左子樹是f4，斷面f7的右子樹是f4，斷面f4的左子樹是f18，斷面f4的右子樹是f18，斷面f4的左子樹是f18，斷面f4的右子樹是f18，斷面f7的左子樹是f4，斷面f7的右子樹是f4，斷面f10的左子樹是f7，斷面f10的右子樹是f7，斷面f7的左子樹是f4，斷面f7的右子樹是f4，斷面f4的左子樹是f18，斷面f4的右子樹是f18，斷面f4的左子樹是f18，斷面f4的右子樹是f18，斷面f7的左子樹是f4，斷面f7的右子樹是f4，斷面f4的左子樹是f18，斷面f4的右子樹是f18，斷面f4的左子樹是f18，斷面f4的右子樹是f5，斷面f5的右子樹是n8，斷面f17的左子樹是f7，斷面f17的右子樹是f7，斷面f7的左子樹是f4，斷面f7的右子樹是f4，斷面f4的左子樹是f18，斷面f4的右子樹是f18，斷面f4的左子樹是f18，斷面f4的右子樹是f18，斷面f7的左子樹是f4，斷面f7的右子樹是f4，斷面f4的左子樹是f18，斷面f4的右子樹是f18，斷面f18的右子樹是f21。
圖11所示為本發(fā)明的一斷面示意圖。下面介紹基于變形場(chǎng)的地層生成。本發(fā)明為每個(gè)斷面建立變形場(chǎng)，在場(chǎng)的作用下地層的離散點(diǎn)依據(jù)某一原則發(fā)生彈性形變，通過變形場(chǎng)將這一側(cè)的地層離散點(diǎn)在斷層兩側(cè)的落差恢復(fù)，模擬為地層產(chǎn)生形變出現(xiàn)此斷層之前的狀態(tài)?？紤]到地層在斷層處的形變是受到已斷裂的裂口影響，所以變形場(chǎng)變形依據(jù)是當(dāng)前斷面和地層已經(jīng)斷裂的裂口斷面。所以采用斷面樹來決定當(dāng)前斷面受哪部分?jǐn)嗝嬗绊?，即要盡量減少恢復(fù)斷層后離散點(diǎn)與實(shí)際離散點(diǎn)的偏差，采用的策略是a、當(dāng)前斷面決定變形的Z值差；b、本斷面及其斷面樹中所有作為其父節(jié)點(diǎn)主斷面的形態(tài)來決定的。
為每個(gè)斷面建立變形場(chǎng)時(shí)，建?？臻g的地質(zhì)元素通過變形場(chǎng)的逐級(jí)作用完成存在斷面時(shí)的建模工作。變形場(chǎng)對(duì)空間斷面的作用有兩個(gè)方向。a，恢復(fù)變形場(chǎng)變形場(chǎng)沿?cái)嗝娴幕謴?fù)方向，在這個(gè)方向上，通過變形場(chǎng)的作用，斷面兩側(cè)的離散點(diǎn)恢復(fù)到斷面不存在的狀態(tài)。b，斷裂變形場(chǎng)變形場(chǎng)沿?cái)嗝鏀嗔训姆较?，在這個(gè)方向上，通過變形場(chǎng)的作用，斷面兩側(cè)的幾何元素由連續(xù)的狀態(tài)斷裂到離散的狀態(tài)。
變形場(chǎng)按斷面逐個(gè)建立變形場(chǎng)，按斷面的面積從小到大逐個(gè)恢復(fù)，當(dāng)前待恢復(fù)斷面一定是斷面樹的葉子節(jié)點(diǎn)的父節(jié)點(diǎn)，斷面變形場(chǎng)建立在斷面的正負(fù)任意一側(cè)，通過恢復(fù)變形場(chǎng)將這一側(cè)的地層離散點(diǎn)變形為未斷裂的狀態(tài)。每建立一個(gè)斷面恢復(fù)變形場(chǎng)就將斷面兩側(cè)的地層離散點(diǎn)恢復(fù)為不存在斷面的地面情況如圖11所示，此時(shí)該斷層面在地質(zhì)模型中就視為不存在了，斷面樹中該斷面節(jié)點(diǎn)也就可以刪除了，該斷面中所有的地層離散都作為葉子節(jié)點(diǎn)掛載在作為其父節(jié)點(diǎn)的一個(gè)葉子子節(jié)點(diǎn)。斷面將地層的離散點(diǎn)所有的斷面恢復(fù)后，生成連續(xù)的曲面，再由斷裂變形場(chǎng)連續(xù)的曲面變形為斷裂的地層面。
本發(fā)明設(shè)計(jì)的變形場(chǎng)滿足如下的特性1.同一斷面的恢復(fù)變形場(chǎng)和斷裂變形場(chǎng)互為逆作用；2.變形場(chǎng)作用于任意一點(diǎn)時(shí)，點(diǎn)不會(huì)跨斷面；3.在斷面上的點(diǎn)沿?cái)嗝孀冃危患?.變形場(chǎng)的范圍以為恢復(fù)的斷面為界，不跨過任意未恢復(fù)的斷面。
每個(gè)變形場(chǎng)力的作用大小和方向由該變形場(chǎng)變形斷面兩側(cè)的地層離散點(diǎn)落差和由變形斷面及其斷面樹上以該斷面作為節(jié)點(diǎn)主斷面的所有節(jié)點(diǎn)的父節(jié)點(diǎn)集合共同作用的方向，該變形場(chǎng)的區(qū)域即為所有父節(jié)點(diǎn)斷面圍成的區(qū)域。
對(duì)于任意一個(gè)斷面建立變形場(chǎng)的過程如下1.根據(jù)斷層兩側(cè)的地層離散點(diǎn)插值出兩片小的地層面；2.根據(jù)兩個(gè)曲面分別求出兩條交線；及3.根據(jù)兩條交線的Z值差，插值出斷層面上Z方向的變形量。
變形量記錄在斷層的每個(gè)網(wǎng)格上，斷面上的變形量是完全變形量，變形場(chǎng)中的其他空間位置的變形量是依據(jù)待計(jì)算點(diǎn)到斷面的投影點(diǎn)的變形量為基礎(chǔ)來計(jì)算的。
將變形場(chǎng)是通過該斷面作為主斷面的樹節(jié)點(diǎn)的左右兩側(cè)離散點(diǎn)的Z值差和本斷面及其斷面樹的中所有以該斷面作為MainFault的所有父節(jié)點(diǎn)主斷面的形態(tài)來決定的。斷面兩側(cè)的離散點(diǎn)Z值差決定變形場(chǎng)的Z值變形量，斷面及其父斷面的形態(tài)決定離散點(diǎn)在x，y方向的變形量。變形場(chǎng)建立好后，通過變形場(chǎng)的作用將地層面的離散點(diǎn)恢復(fù)到在斷面處恢復(fù)為地層撕裂之前的狀態(tài)，變形后的點(diǎn)合并斷面另一側(cè)的點(diǎn)合在一起作為其父節(jié)點(diǎn)建立變形場(chǎng)的采樣點(diǎn)，當(dāng)斷面的離散點(diǎn)恢復(fù)到未斷裂以前時(shí)，對(duì)于為恢復(fù)的斷面來說可以認(rèn)為恢復(fù)的斷面已經(jīng)不存在了，所以所有正在恢復(fù)的斷面都是樹的子節(jié)點(diǎn)。
圖12～圖19為本發(fā)明實(shí)施例的一離散點(diǎn)生成地層的過程。具體為一變形場(chǎng)變形實(shí)際過程地層離散點(diǎn)數(shù)目為2337，斷面數(shù)目為3。變形場(chǎng)的Z方向變形量決定于斷面兩側(cè)的變形量Z和空間點(diǎn)到變形斷面的距離Dis，并認(rèn)為當(dāng)空間點(diǎn)距離斷面的距離大于R時(shí)變形量為0；任意空間點(diǎn)的x，y偏移取決于所有父斷面和當(dāng)前斷面的在Z提升后的x，y偏移值。本方法采用的Z方向的變形量計(jì)算公式如下。其中，Z方向的變形規(guī)則ΔZ′=ΔZ×cos(π2×Disi÷R)0≤Dis≤R]]>(式1-0)X，Y方向的變形公式為 (式1-1)0≤Dis≤R地層通過離散點(diǎn)生成的過程如圖12～圖19所示。地層的恢復(fù)和斷裂的過程并不與實(shí)際斷裂的過程一致，可以存在多種斷面樹的建立過程和變形場(chǎng)的建立方法，只要在變形過程中保證了變形場(chǎng)建立的合理性及正逆變形的準(zhǔn)確性，就可以實(shí)現(xiàn)本方法。其中，圖12所示為初始地層離散點(diǎn)數(shù)據(jù)；圖13所示為恢復(fù)了一個(gè)小裂口后的地層面；圖14所示為恢復(fù)另一個(gè)小裂口后的地層面；圖15所示為完全恢復(fù)了離散點(diǎn)的數(shù)據(jù)后的地層面；圖16所示為根據(jù)變形的離散點(diǎn)生成一個(gè)完整的地層面；圖17所示為撕裂一個(gè)斷面后的地層面；圖18所示為撕裂兩個(gè)斷面后的地層面；及圖19所示為計(jì)算完成后的地層面。
綜上所述，給出了如下實(shí)施例。圖20～圖36所示為本發(fā)明實(shí)施例的一地層生成具體過程示意圖。該圖采用“利津地區(qū)”的數(shù)據(jù)，選取了三個(gè)斷層，該斷層互相相交，并互相裁減，產(chǎn)生斷面的形態(tài)如圖20所示；其斷層的俯視圖見圖21；根據(jù)該斷面的空間位置建立一棵斷面樹，如圖22所示為建立好的斷面樹；將地層的離散點(diǎn)分布到各個(gè)葉子節(jié)點(diǎn)中各個(gè)節(jié)點(diǎn)的離散點(diǎn)，圖23所示為葉子節(jié)點(diǎn)1的空間點(diǎn)的位置示意圖，圖24所示為葉子節(jié)點(diǎn)2的空間點(diǎn)的位置示意圖，圖25所示為葉子節(jié)點(diǎn)3的空間點(diǎn)的位置示意圖，圖26所示為葉子節(jié)點(diǎn)4的空間點(diǎn)的位置示意圖，葉子節(jié)點(diǎn)5不存在空間點(diǎn)；地層離散點(diǎn)分配好后為各個(gè)斷面建立變形場(chǎng)，每個(gè)斷面兩側(cè)地層的離散點(diǎn)變形量，圖27所示給出了斷面F3的變形；圖28所示給出了斷面F2的變形場(chǎng)；圖29所示給出了斷面F1的變形場(chǎng)；根據(jù)變形場(chǎng)的變形，地層離散點(diǎn)統(tǒng)一為未斷裂的地層面，見圖30；其斷裂地層面生成的過程見下圖，圖31、圖32顯示了斷裂斷面F1；圖33、圖34顯示了斷裂斷面F2；圖35、圖36顯示了斷裂斷面F3。
本發(fā)明通過斷面樹來建立變形場(chǎng)的好處是不需要考慮對(duì)當(dāng)前斷面無影響的地層采樣點(diǎn)和斷面，能表達(dá)出斷面之間的邏輯關(guān)系，簡(jiǎn)單易行。雖然地層的恢復(fù)和斷裂的過程并不與實(shí)際斷裂的過程一致，但是在變形過程中保證了變形場(chǎng)建立的合理性及逆變形的準(zhǔn)確性，所以生成的曲面與離散點(diǎn)偏差小，擬合曲面精度高，可實(shí)現(xiàn)過點(diǎn)插值，在連續(xù)處的光滑程度依插值方法而定。本發(fā)明解決了在斷層較多、情況復(fù)雜的地質(zhì)環(huán)境的情況下三維地質(zhì)建模得問題，能夠自動(dòng)生成地層面，有效提高了精確度和效率。
當(dāng)然，本發(fā)明還可有其他多種實(shí)施例，在不背離本發(fā)明精神及其實(shí)質(zhì)的情況下，熟悉本領(lǐng)域的技術(shù)人員當(dāng)可根據(jù)本發(fā)明作出各種相應(yīng)的改變和變形，但這些相應(yīng)的改變和變形都應(yīng)屬于本發(fā)明所附的權(quán)利要求的保護(hù)范圍。
權(quán)利要求
1.一種基于變形場(chǎng)的地層面計(jì)算機(jī)圖形生成方法，用以通過地震數(shù)據(jù)和測(cè)井解釋結(jié)果數(shù)據(jù)進(jìn)行計(jì)算機(jī)處理而自動(dòng)生成地層面圖形，其特征在于，包括步驟一，根據(jù)一地層斷面的空間位置建立一斷面樹；步驟二，根據(jù)所述斷面樹劃分地層離散點(diǎn)；步驟三，由所述地層離散點(diǎn)建立一變形場(chǎng)；步驟四，建立連續(xù)的地層曲面；及步驟五，通過所述變形場(chǎng)的作用力將所述地層曲面變形為實(shí)際地層面。
2.根據(jù)權(quán)利要求1所述的基于變形場(chǎng)的地層面計(jì)算機(jī)圖形生成方法，其特征在于，所述方法進(jìn)一步包括步驟A，根據(jù)地層斷面空間位置建立的斷面樹決定斷面斷裂的先后順序和各斷面之間的影響關(guān)系；步驟B，由所述斷面樹所表示的空間位置關(guān)系將地層的離散點(diǎn)分布到所述斷面樹的各個(gè)葉子節(jié)點(diǎn)；步驟C，通過所述斷面樹和斷面樹中分布的地層離散點(diǎn)逐個(gè)建立斷面變形場(chǎng)，并通過該變形場(chǎng)的逆作用將全部該離散點(diǎn)恢復(fù)為一連續(xù)地層面的離散點(diǎn)；步驟D，由全部恢復(fù)的離散點(diǎn)來建立一連續(xù)的地層曲面；及步驟E，所述變形場(chǎng)按所述先后順序作用于所述地層曲面，將該地層曲面撕裂變形為實(shí)際地層面。
3.根據(jù)權(quán)利要求2所述的基于變形場(chǎng)的地層面計(jì)算機(jī)圖形生成方法，其特征在于，所述步驟A中還包括一建立地層斷面的步驟，該地層斷面由所輸入數(shù)據(jù)中斷面的離散點(diǎn)集合采用計(jì)算機(jī)自由曲面生成算法來建立。
4.根據(jù)權(quán)利要求3所述的基于變形場(chǎng)的地層面計(jì)算機(jī)圖形生成方法，其特征在于，所述輸入數(shù)據(jù)包括地層面的離散點(diǎn)集合、斷面的離散點(diǎn)集合。
5.根據(jù)權(quán)利要求2所述的基于變形場(chǎng)的地層面計(jì)算機(jī)圖形生成方法，其特征在于，所述步驟C中還包括一為建立斷面變形場(chǎng)步驟，包括步驟C1，根據(jù)斷面兩側(cè)的地層離散點(diǎn)插值出兩片小的地層曲面；步驟C2，由所述兩個(gè)地層曲面分別求出兩條交線；及步驟C3，通過所述兩條交線的Z值差，插值出斷層面上Z方向的變形量。
6.根據(jù)權(quán)利要求2所述的基于變形場(chǎng)的地層面計(jì)算機(jī)圖形生成方法，其特征在于，所述先后順序?yàn)閿鄬用娣e大的斷面先斷裂，斷層面積小的斷面后斷裂。
7.根據(jù)權(quán)利要求6所述的基于變形場(chǎng)的地層面計(jì)算機(jī)圖形生成方法，其特征在于，所述步驟C、E中，所述變形場(chǎng)的建立和離散點(diǎn)的恢復(fù)過程按所述先后順序的逆序進(jìn)行；所述實(shí)際地層面的生成過程按所述順序進(jìn)行。
8.根據(jù)權(quán)利要求2所述的基于變形場(chǎng)的地層面計(jì)算機(jī)圖形生成方法，其特征在于，所述斷面變形場(chǎng)是根據(jù)斷層兩側(cè)離散點(diǎn)的落差和斷面的形態(tài)來建立的。
9.根據(jù)權(quán)利要求2所述的基于變形場(chǎng)的地層面計(jì)算機(jī)圖形生成方法，其特征在于，所述步驟E中，所述實(shí)際地層面為一由計(jì)算機(jī)三角網(wǎng)格或自由曲面表示的地層面。
10.根據(jù)權(quán)利要求9述的基于變形場(chǎng)的地層面計(jì)算機(jī)圖形生成方法，其特征在于，所述地層面為一不連續(xù)曲面。
11.根據(jù)權(quán)利要求2所述的基于變形場(chǎng)的地層面計(jì)算機(jī)圖形生成方法，其特征在于，所述斷面樹為一二叉樹。
12.根據(jù)權(quán)利要求11所述的基于變形場(chǎng)的地層面計(jì)算機(jī)圖形生成方法，其特征在于，所述二叉樹的節(jié)點(diǎn)包括節(jié)點(diǎn)主斷面和子斷面列表。
13.根據(jù)權(quán)利要求12所述的基于變形場(chǎng)的地層面計(jì)算機(jī)圖形生成方法，其特征在于，所述節(jié)點(diǎn)的左子節(jié)點(diǎn)及該左子節(jié)點(diǎn)的子節(jié)點(diǎn)在所述節(jié)點(diǎn)主斷面的左側(cè)；所述節(jié)點(diǎn)的右子節(jié)點(diǎn)及該右子節(jié)點(diǎn)的子節(jié)點(diǎn)在所述節(jié)點(diǎn)主斷面的右側(cè)。
全文摘要
本發(fā)明公開了一種基于變形場(chǎng)的地層面計(jì)算機(jī)圖形生成方法，用以通過地震數(shù)據(jù)和測(cè)井解釋結(jié)果數(shù)據(jù)進(jìn)行計(jì)算機(jī)處理而自動(dòng)生成地層面模型，其特征在于，包括步驟一，根據(jù)一地層斷面的空間位置建立一斷面樹；步驟二，根據(jù)所述斷面樹劃分地層離散點(diǎn)；步驟三，由所述地層離散點(diǎn)建立一變形場(chǎng)；步驟四，建立連續(xù)的地層曲面；及步驟五，通過所述變形場(chǎng)的作用力將所述地層曲面變形為實(shí)際地層面。本發(fā)明可處理斷層較多、情況復(fù)雜的地質(zhì)環(huán)境的情況，能夠自動(dòng)生成地層面，有效提高生成地層面的精確度和效率。
文檔編號(hào)G01V1/28GK1825139SQ20051001133
公開日2006年8月30日 申請(qǐng)日期2005年2月21日 優(yōu)先權(quán)日2005年2月21日
發(fā)明者楊欽, 洪錦玲, 張俊安, 劉瑞剛, 金宇林, 董艷秋 申請(qǐng)人:北京航空航天大學(xué)