專利名稱：含逆斷層等值線繪圖方法
技術(shù)領(lǐng)域：
本發(fā)明涉及地震資料采集、處理與解釋領(lǐng)域，更具體地講，涉及一種含逆斷層等值線繪圖方法。
背景技術(shù)：
等值線圖是在石油勘探開發(fā)、采礦、地質(zhì)、地球物理、地球化學(xué)和氣象等領(lǐng)域等工程和技術(shù)領(lǐng)域內(nèi)應(yīng)用極廣的一種圖形，是眾多領(lǐng)域成果表示的重要圖件之一。所謂繪制等值線圖就是在二維平面上，把大量離散的、又具有一定規(guī)律的幾何量值或物理量值，用數(shù)學(xué)的方法插值并將具有相同量值的點轉(zhuǎn)換成圖形的過程。等值線圖是最主要也是最常用的一種圖件，國內(nèi)外許多學(xué)者對此進(jìn)行了大量的研究，但目前各種軟件獲得的等值線圖通常難以令人滿意。這一方面是由于已知信息量有限， 且分布不均；另一方面由于研究區(qū)常常存在許多的斷層特別是逆斷層，使得構(gòu)造非常復(fù)雜。 這都使得計算機(jī)自動繪制等值線圖變得很困難。目前我國石油地球物理勘探中使用的等值線繪制軟件大都是國外進(jìn)口的軟件，但當(dāng)研究區(qū)斷層較復(fù)雜時特別是存在逆斷層時，常常繪制不出令地質(zhì)學(xué)家滿意的等值線圖，這是因為傳統(tǒng)的等值線繪制只能繪制不含逆斷層的等值線，即等值線在空間上不重疊；如果要繪制重疊的等值線，只有先繪制其中上盤部分， 再繪制下盤部分，然后用人工重疊。這給工程上繪制等值線圖非常不便。

發(fā)明內(nèi)容
本發(fā)明提供了一種基于空間三角網(wǎng)格的含逆斷層等值線繪圖方法，所述方法采用空間網(wǎng)格化算法，以逆斷層的范圍和層次關(guān)系作為約束條件，對層位數(shù)據(jù)進(jìn)行三角網(wǎng)格化； 對空間網(wǎng)格化數(shù)據(jù)進(jìn)行追蹤，實現(xiàn)含逆斷層的等值線繪圖。本發(fā)明基于空間三角網(wǎng)格的含逆斷層等值線繪圖方法，利用斷層多邊形和數(shù)據(jù)的關(guān)聯(lián)關(guān)系，以及斷層多邊形的層次關(guān)系，通過分點集的三角形連接，實現(xiàn)含逆斷層的構(gòu)造圖繪制。此種方法繪制的構(gòu)造圖效果好、精度高，能夠?qū)⒛鏀鄬拥南卤P信息繪制出來，有效地解決了含逆斷層的構(gòu)造圖繪制問題。此種方法能夠廣泛地應(yīng)用于石油地震勘探、地質(zhì)地礦和建筑工程地況調(diào)查。根據(jù)本發(fā)明的一方面，提供了一種含逆斷層等值線繪圖方法，所述方法包括下述步驟劃分?jǐn)鄬佣噙呅蝺?nèi)外圍的數(shù)據(jù)；采用空間網(wǎng)格化算法對數(shù)據(jù)進(jìn)行三角網(wǎng)格化以得到最終的三角網(wǎng)；對最終的三角網(wǎng)進(jìn)行等值線追蹤操作。根據(jù)本發(fā)明的另一方面，所述劃分?jǐn)鄬佣噙呅蝺?nèi)外圍的數(shù)據(jù)的步驟可包括獲取包括上下盤線段和斷層與斷層之間的位置關(guān)系的斷層數(shù)據(jù)；通過上下盤線段擬合出斷層多邊形邊界數(shù)據(jù)，并通過層位數(shù)據(jù)和斷層數(shù)據(jù)之間的關(guān)系計算出斷層多邊形的ζ值，其中，ζ 值表示時間值或深度值；將各個斷層多邊形進(jìn)行多邊形與多邊形之間的裁剪，獲得多邊形與多邊形的并集或者交集作為結(jié)果斷層多邊形，并記錄各個結(jié)果斷層多邊形之間的位置關(guān)系；根據(jù)各個結(jié)果斷層多邊形之間的位置關(guān)系，確定各個結(jié)果斷層多邊形中每個點在ζ值上的順序關(guān)系。根據(jù)本發(fā)明的另一方面，所述劃分多邊形內(nèi)外圍的數(shù)據(jù)的步驟還可包括獲取層位數(shù)據(jù)，并對層位數(shù)據(jù)段進(jìn)行遍歷，其中，在對層位數(shù)據(jù)段進(jìn)行遍歷的過程中，判斷當(dāng)前層位數(shù)據(jù)段是否與斷層相關(guān)聯(lián)；如果當(dāng)前層位數(shù)據(jù)段與斷層關(guān)聯(lián)，則獲取有封閉關(guān)系的相應(yīng)的結(jié)果斷層多邊形的數(shù)據(jù)，并對所述相應(yīng)的結(jié)果斷層多邊形進(jìn)行遍歷；否則，將當(dāng)前層位數(shù)據(jù)段的層位數(shù)據(jù)放入最外層的多邊形的結(jié)果數(shù)據(jù)中。根據(jù)本發(fā)明的另一方面，在對所述相應(yīng)的結(jié)果斷層多邊形進(jìn)行遍歷的過程中，可使用線段和多邊形之間的裁剪算法獲取層位數(shù)據(jù)段在相應(yīng)的結(jié)果斷層多邊形內(nèi)部的數(shù)據(jù)， 并根據(jù)這個結(jié)果斷層多邊形與其他的結(jié)果斷層多邊形的ζ值，獲取層位數(shù)據(jù)段的層次關(guān)系。根據(jù)本發(fā)明的另一方面，如果對所述相應(yīng)的結(jié)果斷層多邊形的遍歷結(jié)束，則獲取裁剪剩余的層位數(shù)據(jù)段的數(shù)據(jù)，并將層位數(shù)據(jù)段的數(shù)據(jù)直接放入最外層的多邊形的結(jié)果數(shù)據(jù)中。根據(jù)本發(fā)明的另一方面，所述采用空間網(wǎng)格化算法對數(shù)據(jù)進(jìn)行三角網(wǎng)格化以得到最終的三角網(wǎng)的步驟可包括對斷層外的外圍數(shù)據(jù)進(jìn)行構(gòu)網(wǎng)；刪除斷層內(nèi)的三角形，并對斷層內(nèi)的數(shù)據(jù)單獨進(jìn)行構(gòu)網(wǎng)；將對斷層外的外圍數(shù)據(jù)進(jìn)行構(gòu)網(wǎng)所得的三角網(wǎng)與對斷層內(nèi)的數(shù)據(jù)進(jìn)行構(gòu)網(wǎng)所得的三角網(wǎng)合并，獲得最終的三角網(wǎng)。根據(jù)本發(fā)明的另一方面，可使用狄洛尼三角剖分算法對斷層內(nèi)的數(shù)據(jù)和斷層外的外圍數(shù)據(jù)進(jìn)行三角形構(gòu)網(wǎng)。根據(jù)本發(fā)明的另一方面，所述對最終的三角網(wǎng)進(jìn)行等值線追蹤操作的步驟可包括如果給定等值線的高程值介于三角形的邊對應(yīng)的兩個端點的高程值之間，則對該邊進(jìn)行標(biāo)記，并以該邊上相應(yīng)高程值的點作為等值線的起始點；在該邊所在的三角形的其他邊中尋找相應(yīng)高程值的點作為出口 ；重復(fù)在包含出口所在的邊的相鄰的三角形的其他邊中尋找該相應(yīng)高程值的點作為新的出口的操作，直至返回到起始點。


圖1是示出根據(jù)本發(fā)明示例性實施例的含逆斷層等值線繪圖方法的流程圖。圖2是根據(jù)本發(fā)明示例性實施例的劃分?jǐn)?shù)據(jù)的數(shù)據(jù)結(jié)構(gòu)的示例。圖3是示出根據(jù)本發(fā)明示例性實施例的劃分?jǐn)?shù)據(jù)的方法的流程圖。圖4至圖6是示出對數(shù)據(jù)進(jìn)行三角形構(gòu)網(wǎng)的方法的示意圖。圖7是根據(jù)本發(fā)明示例性實施例的等值線追蹤方法的示意圖。
具體實施例方式下面將結(jié)合附圖來詳細(xì)描述本發(fā)明的示例性實施例。圖1是根據(jù)本發(fā)明示例性實施例的含逆斷層等值線繪圖方法的流程圖。根據(jù)本發(fā)明示例性實施例的含逆斷層等值線繪圖方法的基本構(gòu)思是采用空間網(wǎng)格化算法，以逆斷層的范圍和層次關(guān)系作為約束條件，對層位數(shù)據(jù)進(jìn)行三角網(wǎng)格化；對空間網(wǎng)格化數(shù)據(jù)進(jìn)行追蹤，實現(xiàn)含逆斷層的等值線繪圖。主要技術(shù)內(nèi)容包括含逆斷層的空間三角形網(wǎng)格化算法和基于三角網(wǎng)格的等值線追蹤算法。
為了描述簡便，下面所描述的斷層包含逆斷層。參照圖1，根據(jù)本發(fā)明的含逆斷層等值線繪圖方法包括下面所述的步驟在步驟 S10，劃分?jǐn)鄬佣噙呅蝺?nèi)外圍的數(shù)據(jù)。在步驟S20，采用空間網(wǎng)格化算法對數(shù)據(jù)進(jìn)行三角網(wǎng)格化以得到最終的三角網(wǎng)。在步驟S30，對最終的三角網(wǎng)進(jìn)行等值線追蹤操作。下面將參照圖 2至圖7來詳細(xì)描述根據(jù)本發(fā)明示例性實施例的含逆斷層的等值線繪圖方法的流程。這里，如何劃分?jǐn)鄬佣噙呅蝺?nèi)外圍的數(shù)據(jù)就是本算法的核心部分。圖2是根據(jù)C++ 編程語言的劃分?jǐn)?shù)據(jù)的數(shù)據(jù)結(jié)構(gòu)的示例；當(dāng)然，可以以各種語言來定義本發(fā)明的數(shù)據(jù)結(jié)構(gòu)。在圖2中，SpacePoint結(jié)構(gòu)為基本的數(shù)據(jù)結(jié)構(gòu)。其中x，y，ζ表示一個點的三維大地坐標(biāo)；index表示在所有的結(jié)果數(shù)據(jù)中的索引位置；faultageFlag表示這個點是否為斷層多邊形上面的點，0表示這個點不在多邊形邊界上，1表示這個點在多邊形邊界上； IayerFlag表示點的層次關(guān)系，如果為最上面一層的點，取值為0，第二層為1，依次類推。BrokenLine和Points是由一系列的點構(gòu)成的線段。BlockPoints包含了一個多邊形中的點結(jié)構(gòu)，其中，points包含了多邊形中的所有點(包括多邊形邊界上的點和層位數(shù)據(jù)的點)；externalEdge表示多邊形的外邊界。 internalEdges表示多邊形內(nèi)部包含的多邊形的邊界，也叫內(nèi)邊界，可能同時存在多個內(nèi)邊界。InputData是最終劃分成功的數(shù)據(jù)，其中，allPoints包含了所有的點，這些點中的index在這些數(shù)據(jù)里面是唯一的；externalPoints表示斷層多邊形外部的數(shù)據(jù)； internalPoints表示斷層多邊形內(nèi)部的數(shù)據(jù)。InputData是一種嵌套結(jié)構(gòu)，因此通過遞歸遍歷劃分出的數(shù)據(jù)，可以生成子網(wǎng)。圖3是示出根據(jù)本發(fā)明示例性實施例的劃分?jǐn)鄬佣噙呅蝺?nèi)外圍的數(shù)據(jù)的方法的流程圖，下面將結(jié)合圖2和圖3來詳細(xì)介紹根據(jù)本發(fā)明示例性實施例的劃分?jǐn)?shù)據(jù)的方法。在步驟S310，首先獲取所有的斷層數(shù)據(jù)，包括斷層與斷層之間的位置關(guān)系和上下盤線段。在步驟S320，通過上下盤線段擬合出斷層多邊形邊界數(shù)據(jù)，并通過層位數(shù)據(jù)和斷層數(shù)據(jù)之間的關(guān)系計算出斷層多邊形的ζ值(也即斷層多邊形的時間或者深度值)。在步驟S330，將各個斷層多邊形進(jìn)行多邊形與多邊形之間的裁剪，獲得多邊形與多邊形的并集或者交集作為結(jié)果斷層多邊形，并記錄各個結(jié)果斷層多邊形之間的位置關(guān)系。在步驟S340，根據(jù)已經(jīng)獲取的各個結(jié)果斷層多邊形之間的位置關(guān)系，確定各個結(jié)果斷層多邊形中每個點在ζ值方向上的順序關(guān)系，并將結(jié)果斷層多邊形的數(shù)據(jù)加入到具有InputData數(shù)據(jù)結(jié)構(gòu)的結(jié)果數(shù)據(jù)中。其中每個裁剪之后的結(jié)果斷層多邊形都將有一個 BrokenLine實例；同時，每個結(jié)果斷層多邊形在hputData. internalPoints中都將占有一個實例。在步驟S350，獲取所要劃分的層位數(shù)據(jù)，并對層位數(shù)據(jù)段進(jìn)行遍歷。在對層位數(shù)據(jù)段進(jìn)行遍歷的過程中，在步驟S360判斷遍歷層位數(shù)據(jù)段是否結(jié)束。如果結(jié)束，則在步驟S370結(jié)束數(shù)據(jù)的劃分；否則，在步驟S380中進(jìn)一步判斷當(dāng)前層位數(shù)據(jù)段是否與斷層相關(guān)聯(lián)(即，判斷此層位數(shù)據(jù)段的封閉關(guān)系)。如果當(dāng)前層位數(shù)據(jù)段與斷層沒有關(guān)聯(lián)(即，不存在封閉關(guān)系)，則在步驟S390將當(dāng)
6前層位數(shù)據(jù)段的層位數(shù)據(jù)作為最外層的多邊形的結(jié)果數(shù)據(jù)放入斷層多邊形外部的數(shù)據(jù)，即 InputData數(shù)據(jù)結(jié)構(gòu)中的externalPoints中；同時，也將其放入allPoints中，并返回到步驟S360以對下一層位數(shù)據(jù)段重復(fù)上述操作。如果當(dāng)前層位數(shù)據(jù)段與斷層關(guān)聯(lián)(即，該層位數(shù)據(jù)段存在封閉關(guān)系)，則在步驟 S3100中獲取與當(dāng)前層位數(shù)據(jù)段有封閉關(guān)系的各個結(jié)果斷層多邊形的數(shù)據(jù)，并遍歷所述各個結(jié)果斷層多邊形。在對所述各個結(jié)果斷層多邊形進(jìn)行遍歷的過程中，在步驟S3110中判斷對與當(dāng)前層位數(shù)據(jù)段具有封閉關(guān)系的各個結(jié)果斷層多邊形的遍歷是否結(jié)束。如果沒有結(jié)束，則在步驟S3120中使用線段和多邊形之間的裁剪算法獲取當(dāng)前層位數(shù)據(jù)段在當(dāng)前遍歷到的結(jié)果斷層多邊形內(nèi)部的數(shù)據(jù)(交集)。然后在步驟S3130中，獲取當(dāng)前遍歷到的結(jié)果斷層多邊形在具有InputData數(shù)據(jù)結(jié)構(gòu)的結(jié)果數(shù)據(jù)中的internalPoints中的索引位置，并根據(jù)這個結(jié)果斷層多邊形與其他的結(jié)果斷層多邊形的ζ值，獲取這個層位數(shù)據(jù)段的層次關(guān)系，并將獲取的層次關(guān)系寫入IayerFlag信息中，將數(shù)據(jù)放入索引位置對應(yīng)的斷層多邊形的結(jié)果數(shù)據(jù)中的 interalPoints以及allPoints中，然后返回步驟S3110對下一個斷層多邊形進(jìn)行操作。如果遍歷結(jié)果斷層多邊形結(jié)束，則執(zhí)行步驟S3140，獲取裁剪剩余的層位數(shù)據(jù)段的數(shù)據(jù)，并轉(zhuǎn)向執(zhí)行步驟S390將層位數(shù)據(jù)段的數(shù)據(jù)直接放入最外層的多邊形的結(jié)果數(shù)據(jù)中。在根據(jù)圖2和圖3所示的數(shù)據(jù)劃分的基礎(chǔ)上，含逆斷層的數(shù)據(jù)被分為多個點集的數(shù)據(jù)，每個點集的數(shù)據(jù)在平面上不重合。根據(jù)本發(fā)明示例性實施例的空間三角網(wǎng)格化算法采用了分塊的思想，具體步驟如下所示首先，如圖4所示將斷層內(nèi)的數(shù)據(jù)挖空，先對斷層外的外圍數(shù)據(jù)進(jìn)行構(gòu)網(wǎng)；這里，由于外圍數(shù)據(jù)一般來說數(shù)據(jù)量比較大，因此，在外圍數(shù)據(jù)連網(wǎng)上面，同樣可以采用分塊算法來提高速度。然后如圖5所示，刪除斷層內(nèi)的三角形，并且如圖6所示，對于斷層內(nèi)的數(shù)據(jù)進(jìn)行單獨構(gòu)網(wǎng)。最后將斷層外的外圍數(shù)據(jù)進(jìn)行構(gòu)網(wǎng)所得的三角網(wǎng)與對斷層內(nèi)數(shù)據(jù)進(jìn)行構(gòu)網(wǎng)所得的三角網(wǎng)合并，就是最終的三角網(wǎng)，其中，將圖6中所示的上層的數(shù)據(jù)網(wǎng)格的上盤與外圍數(shù)據(jù)相連，下層的數(shù)據(jù)網(wǎng)格的下盤與外圍數(shù)據(jù)相連，最終構(gòu)成一個完整的三角網(wǎng)。在對斷層內(nèi)的數(shù)據(jù)以及斷層外的外圍數(shù)據(jù)進(jìn)行構(gòu)網(wǎng)的過程中，首先，根據(jù)斷層多邊形與測線的相交關(guān)系，將層位數(shù)據(jù)分堆處理；其次，對數(shù)據(jù)進(jìn)行遍歷，根據(jù)數(shù)據(jù)的大小劃分區(qū)間，在每一個區(qū)間使用狄洛尼三角剖分算法進(jìn)行三角形連接(構(gòu)網(wǎng))；再次，遍歷所有的分區(qū)，將分區(qū)之間的空隙采用三角形連接。在完成三角形連網(wǎng)之后，執(zhí)行圖1所示的步驟S30的基于三角網(wǎng)的等值線追蹤操作。下面將結(jié)合圖7來描述等值線追蹤操作的具體步驟。在對三角網(wǎng)的邊進(jìn)行遍歷的過程中，執(zhí)行下述所述的操作(1)若給定等值線的高程值介于三角形的邊對應(yīng)的兩個端點的高程值之間，則對該邊進(jìn)行標(biāo)記并以該條邊上相應(yīng)高程值的點作為等值線的起始點(圖7的(a)) ； (2)在該邊所在的三角形的其他邊中尋找相應(yīng)高程值的點作為出口，并記錄具有出口的邊所在的相鄰的三角形(圖7的(b)) ； (3) 在該相鄰的三角形的其他邊中尋找該相應(yīng)高程值點作為新的出口(圖7的(C))，這樣直至最后返回起始點，則這條等值線追蹤結(jié)束(圖7的(d))。如果在相鄰三角形中沒有找到相應(yīng)高程值的點，則該條等值線追蹤結(jié)束。
根據(jù)本發(fā)明的示例性實施例，利用斷層多邊形和數(shù)據(jù)的關(guān)聯(lián)關(guān)系，以及斷層多邊形的層次關(guān)系，通過分點集的三角形連接，實現(xiàn)含逆斷層的構(gòu)造圖繪制。此種方法繪制的構(gòu)造圖效果好、精度高，能夠?qū)⒛鏀鄬拥南卤P信息繪制出來，有效地解決了含逆斷層的構(gòu)造圖繪制問題。雖然已經(jīng)參照本發(fā)明的示例性實施例具體描述和顯示了本發(fā)明，但是本領(lǐng)域的普通技術(shù)人員應(yīng)該理解，在不脫離由權(quán)利要求限定的本發(fā)明的精神和范圍的情況下，可以對其進(jìn)行形式和細(xì)節(jié)的各種改變。
權(quán)利要求
1.一種含逆斷層等值線繪圖方法，所述方法包括下述步驟 劃分?jǐn)鄬佣噙呅蝺?nèi)外圍的數(shù)據(jù)；采用空間網(wǎng)格化算法對數(shù)據(jù)進(jìn)行三角網(wǎng)格化以得到最終的三角網(wǎng)； 對最終的三角網(wǎng)進(jìn)行等值線追蹤操作。
2.如權(quán)利要求1所述的方法，其中，所述劃分?jǐn)鄬佣噙呅蝺?nèi)外圍的數(shù)據(jù)的步驟包括 獲取包括斷層與斷層之間的位置關(guān)系和上下盤線段的斷層數(shù)據(jù)；通過上下盤線段擬合出斷層多邊形邊界數(shù)據(jù)，并通過層位數(shù)據(jù)和斷層數(shù)據(jù)之間的關(guān)系計算出斷層多邊形的ζ值，其中，ζ值表示時間值或深度值；將各個斷層多邊形進(jìn)行多邊形與多邊形之間的裁剪，獲得多邊形與多邊形的并集或者交集作為結(jié)果斷層多邊形，并記錄各個結(jié)果斷層多邊形之間的位置關(guān)系；根據(jù)各個結(jié)果斷層多邊形之間的位置關(guān)系，確定各個結(jié)果斷層多邊形中每個點在ζ值上的順序關(guān)系。
3.如權(quán)利要求2所述的方法，其中，所述劃分多邊形內(nèi)外圍的數(shù)據(jù)的步驟還包括 獲取層位數(shù)據(jù)，并對層位數(shù)據(jù)段進(jìn)行遍歷，其中，在對層位數(shù)據(jù)段進(jìn)行遍歷的過程中，判斷當(dāng)前層位數(shù)據(jù)段是否與斷層相關(guān)聯(lián)；如果當(dāng)前層位數(shù)據(jù)段與斷層關(guān)聯(lián)，則獲取有封閉關(guān)系的相應(yīng)的結(jié)果斷層多邊形的數(shù)據(jù)，并對所述相應(yīng)的結(jié)果斷層多邊形進(jìn)行遍歷；如果當(dāng)前層位數(shù)據(jù)段與斷層沒有關(guān)聯(lián)，將當(dāng)前層位數(shù)據(jù)段的層位數(shù)據(jù)放入最外層的多邊形的結(jié)果數(shù)據(jù)中。
4.如權(quán)利要求3所述的方法，其中，在對所述相應(yīng)的結(jié)果斷層多邊形進(jìn)行遍歷的過程中，使用線段和多邊形之間的裁剪算法獲取層位數(shù)據(jù)段在相應(yīng)的結(jié)果斷層多邊形內(nèi)部的數(shù)據(jù)，并根據(jù)這個結(jié)果斷層多邊形與其他的結(jié)果斷層多邊形的ζ值，獲取層位數(shù)據(jù)段的層次關(guān)系。
5.如權(quán)利要求3所述的方法，其中，如果對所述相應(yīng)的結(jié)果斷層多邊形的遍歷結(jié)束，則獲取裁剪剩余的層位數(shù)據(jù)段的數(shù)據(jù)，并將層位數(shù)據(jù)段的數(shù)據(jù)直接放入最外層的多邊形的結(jié)果數(shù)據(jù)中。
6.如權(quán)利要求1所述的方法，其中，所述采用空間網(wǎng)格化算法對數(shù)據(jù)進(jìn)行三角網(wǎng)格化以得到最終的三角網(wǎng)的步驟包括對斷層外的外圍數(shù)據(jù)進(jìn)行構(gòu)網(wǎng)；刪除斷層內(nèi)的三角形，并對斷層內(nèi)的數(shù)據(jù)單獨進(jìn)行構(gòu)網(wǎng)；將對斷層外的外圍數(shù)據(jù)進(jìn)行構(gòu)網(wǎng)所得的三角網(wǎng)與對斷層內(nèi)的數(shù)據(jù)進(jìn)行構(gòu)網(wǎng)所得的三角網(wǎng)合并，獲得最終的三角網(wǎng)。
7.如權(quán)利要求6所述的方法，其中，使用狄洛尼三角剖分算法對斷層內(nèi)的數(shù)據(jù)和斷層外的外圍數(shù)據(jù)進(jìn)行三角形構(gòu)網(wǎng)。
8.如權(quán)利要求1所述的方法，其中，所述對最終的三角網(wǎng)進(jìn)行等值線追蹤操作的步驟包括如果給定等值線的高程值介于三角形的邊對應(yīng)的兩個端點的高程值之間，則對該邊進(jìn)行標(biāo)記，并以該邊上相應(yīng)高程值的點作為等值線的起始點；在該邊所在的三角形的其他邊中尋找相應(yīng)高程值的點作為出口 ； 重復(fù)在包含出口所在的邊的相鄰的三角形的其他邊中尋找該相應(yīng)高程值的點作為新的出口的操作，直至返回到起始點。
全文摘要
提供了一種含逆斷層等值線繪圖方法。所述方法包括下述步驟劃分?jǐn)鄬佣噙呅蝺?nèi)外圍的數(shù)據(jù)；采用空間網(wǎng)格化算法對數(shù)據(jù)進(jìn)行三角網(wǎng)格化以得到最終的三角網(wǎng)；對最終的三角網(wǎng)進(jìn)行等值線追蹤操作。此種方法繪制的構(gòu)造圖效果好、精度高，能夠?qū)⒛鏀鄬拥南卤P信息繪制出來，有效地解決了含逆斷層的構(gòu)造圖繪制問題。
文檔編號G01V1/30GK102385067SQ201110220729
公開日2012年3月21日 申請日期2011年8月3日 優(yōu)先權(quán)日2011年8月3日
發(fā)明者吳秋波, 姚興苗, 巫盛洪, 張洞君, 鄒文, 陳小二, 陶正喜, 黃東山 申請人:中國石油集團(tuán)川慶鉆探工程有限公司地球物理勘探公司