工程巖體三維空間結構建模與關鍵塊識別方法

2023-05-06 18:51:36

專利名稱：：工程巖體三維空間結構建模與關鍵塊識別方法
技術領域：
：本發明屬於巖土工程、礦山開採、工程建築、水利工程、地下結構等巖體工程的實體結構建模和分析領域，尤其涉及到一種基於不連續介質、各向異性巖體、工程擾動過程三維數值模擬分析方法與工程應用。
背景技術：
：工程巖體穩定性影響因素比較複雜，與傳統結構穩定可靠性分析存在較大的差別，它與工程巖體的地質特徵、工程因素、力學機制、破壞機理、擾動過程等各方面的因素緊密相關。在漫長的地質歷史演變過程中，巖石經受了不同大小和不同方向的多期構造應力場的長期作用，加之巖石材料本身的非均質性，從而導致了巖體結構面分布的隨機性，形態的多樣性，空間組合的複雜性，密度的不均勻性，以至於對巖體的綜合性能變得非常複雜，僅僅利用有限幾個統計量如均值、方差、階矩，難以全面描述巖體特性，因此在巖體性能方面需要更先進的分析手段和方法。傳統的分析方法是基於有限元、離散元、邊界元及其耦合等數值方法考慮，模擬過程很難考慮巖體結構的真實特徵和工程的擾動過程，目前，大多採取塊體理論解決上述問題，但是塊體理論主要針對邊坡塊體穩定的分析，研究對象是邊坡，並且缺乏如赤平投影等塊體理論的基本分析功能，因此不是真正意義上的塊體理論計算方法，其它方法如楔形體穩定性分析等，也與塊體理論方法有較大的差別。
發明內容基於現場工程巖體結構特性，建立反映工程巖體真實結構特性三維模型，實現了工程巖體結構和施工過程結合起來模擬真實巖體結構，本發明中同時實現了三種塊體搜索與顯示的方法，即確定性關鍵塊體搜索與顯示、非確定性結構面和隨機塊體的空間模擬、確定性和非確定性結構面模型顯示與分析，同時，通過對現場數據的對比以及現場觀測數據的反饋，可以進一步預測工程巖體破裂與失穩的狀態；工程巖石塊體通常分為兩類關鍵塊體和一般塊體，關鍵塊體就是不用其他塊體阻礙，自己能塌下來、掉下來的塊體；即在工程作用力和自重作用下，由於滑動面上的抗剪強度不足以抵禦滑動力，若不施加工程錨固措施，必將失穩的塊體。工程巖體三維空間結構建模與關鍵塊識別方法步驟如下如圖1所示，—、工程巖體現場結構原始數據獲取，數據包括工程巖體結構參數、工程巖體結構面信息和工程巖體結構狀態，其中工程巖體結構參數包括巷道、隧洞和斷面尺寸；工程巖體結構面信息和工程巖體結構狀態包括斷層參數和節理參數；二、工程巖體原始數據的處理與提取通過巖體結構面的結構特點會產生多種原始數據，同時，在對這些數據進行處理過程中，也會產生大量的中間數據，如何對這些數據進行科學的管理和有效的處理，以充分利用和挖掘其中的有效信息，從而保證成果的準確可靠性，是一個關鍵問題，工程巖體裂隙發育的特點決定了對工程巖體裂隙的研究必須建立在大量的統計分析基礎之上，以實現對研究區具體工程部位的工程巖體裂隙進行優勢方位分析及配套參數提取；在現場採集的原始數據基礎上，利用赤平投影和概率統計方法，整理出1)巖體裂隙方位，以及每一條裂隙所在的硐號、方位、層位、工程部位、風化帶信息，進行裂隙優勢方位分析；2)存儲工程普遍測網法獲得的每條裂隙的基本信息，包括每條裂隙對應的測點號、位置和產狀，進行繪製裂隙平面展布圖的坐標信息；3)存儲工程巖體測點的各條裂隙的產狀及相對測網原點的局部幾何坐標；在對普遍測網法實測數據有效管理的基礎上，根據現場實測數據對各測點的巖體結構參數由計算機進行數據處理和統汁分析。在此基礎上，為本發明提供的初始數據，工程地下結構參數和結構面參數的數據範圍及精度要求1)工程地下結構參數tableseeoriginaldocumentpage52)結構面參數tableseeoriginaldocumentpage5根據前述的裂隙幾何數據對表徵巖體結構指標進行有效的測試，包括結構面的優勢方位、間距、密度、連通率、各類跡長、測網延時質量指標(RQD)值、裂隙係數、三維空間結構模式、裂隙網絡分形維數、損傷張量、滲透張量等十餘個巖體結構參數，並具有對選擇的多個測點進行綜合統計分析和結構面參數概率分布擬合及裂隙網絡模擬計算，包括1)基本分析計算①測點分布及坐標顯示②巖體裂隙網絡顯示及三維結構可視化模型③優勢方位計算④間距計算⑤連通率計算⑥跡長計算⑦裂隙率計算；2)網參數統計①測網延時質量指標(RQD)值統計②測網分形維數統③測網滲透張量統計④測網損傷張量統計；3)參數統計計算；4)結構面網絡模擬功能。三、構建工程巖體模型；1)、建立巖體模型判斷巖體幾何體被分成凸體和凹體兩類，建立模型時加以區別；凸體結構首先，把它分解為多邊形，依次添加多邊形，多邊形的建立是按同一順序即逆時針或順時針順序，依次輸入頂點坐標到文本文件中，把實際工程巖體劃分為凸體遵循一個原則，即巖體被劃分後，任何兩個相鄰的巖體公面必須是重合的，如圖4(a)和4(b)所示，巖體被分成A、B、C、D、E、F六個凸體，其中相鄰的塊間公共面都是完全重合的；凹體結構按照相鄰兩個塊體公共面重合的原則，把凹體分成多個凸體，並把重合面賦上虛擬面的標記值，將被分成的凸體按照凸體的建模規則建模；2)、建立數據結構A數據結構如圖2、圖3(a)、圖3(b)和圖(c)所示，先確定塊體稜線的兩個端點，把放入線段存儲結構中，再把形成的線段放入多邊形存儲結構中，最後把多邊形放入塊體的數據結構中；B建立數據組織，存儲數據，形成數據鏈，採用程式語言中自帶的動態類對數據進行管理；3)、確定結構面信息方法包括確定性結構面、非確定性結構面和耦合結構面確定性結構面現場調查所得的結構面，包括地球物理勘探、數字圖象拍攝等手段直接採集的結構面數據，作為確定性結構面數據，方位、傾向、傾角和跡長，輸入模型中；非確定性結構面通過概率統計的方法來模擬產生的結構面數據，根據現場鄰近工程、專家經驗，採用統計方法獲得結構面參數；耦合結構面模型中輸入工程中確定性結構面，將不能通過鑽孔、照相、物理方法獲取的數據，用概率統計的方法模擬產生結構面數據；四、塊體識別方法把巖體模型劃分為有限個網格，然後在依次加入結構面，對現有的小塊體進行切割，最後在去掉網格，合併小的塊體，從而識別出有結構面切割得到的所有塊體；1)、網格劃分根據區域內結構面的密度，平均半徑進行網格劃分，首先確定區域的範圍，然後設置三組相互垂直的虛擬結構面對巖體模型進行切割，最後記錄下小單元的數據結構，如圖5(a)、5(b)和圖6所示；①切割過程裂隙面對塊體切割可分解為對其表面多邊形的切割，進而分解為裂隙面對直線的切割，實現裂隙與多邊形的切割，裂隙切割多邊形；有三種情況A如果多邊形頂點全在裂隙面的上方，此時把這個多邊形編號後放到上面的新塊體中；B如果多邊形頂點全在裂隙面的下方，此時把這個多邊形編號後放到下面的新塊體中；C如果多邊形頂點分布在裂隙面的兩側，用裂隙平面把這個多邊形切割成兩個新多邊形，編號後分別放入上下兩個塊體中；當塊體表面多邊形屬於情況B時，我們進一步研究裂隙面如何對直線切割，如圖7中，以DG為例，設裂隙面所在的平面方程為Ax+By+Cz+D=0，點D坐標為(Xl，yi，Zl)，點G坐標為(x2，y2，z2)，分別代入方程中結果為Ax一By,Cz,D>0Ax2+By2+Cz2+D<0所以點D在平面上，點G在平面下，求出DG所在直線與裂隙平面的交點K，然後把KD相關數據放入上面塊體相應的多邊形中，KG相關數據放入下面塊體相應的多邊形中；②採用這個裂隙切割其他與其相關的塊體，直到切割完與其相關的所有塊體，然後進入下一個裂隙的切割過程；如圖8所示；2)添加結構面網格換分後，加入實測結和擬合出來的結構面，對現有的小塊體進行再次切割；3)開挖面處理，近似用平面代替曲面，曲面對塊體的切割就是對單元塊體的切割，求出單元體所有稜線與曲面的交點後依次連接各點，形成一個新的面，生成兩個新的塊體；如圖9所示，曲面對塊體的切割就是對單元塊體的切割，求出單元體所有稜線與曲面的交點後依次連接各點，形成一個新的面，從而形成兩個新的塊體，這樣，得到的開挖面就隨著網格數量的增加而無限接近實際的開挖曲面；4)、消除網格，生成塊體首先判斷相鄰兩個小單元面的交集是否與虛擬結構面有交集，如果有交集則合併兩個小塊，沒有交集就不必合併，按照此方法，依次對存在的所有小塊體進行處理，最後記錄下合併後塊體的數據；如圖10、圖ll(a)和圖ll(b)所示；5)、判別塊體的可移移動性首先，找出一個帶有兩個角度參數的試算向量，這兩個參數的初始值為90。和0°，依次試算，如果這個試算向量與塊體每個面的方向向量相乘都為正，則判斷次塊體幾何可移動，進一步利用莫爾庫侖強度理(Mohr-Coulomb)準則對其進行穩定性判斷，如果不穩定，就是關鍵塊體；經過網格劃分，結構面切割塊體，消除虛擬結構面合併塊體的運算後，所有獨立塊體的數據都已經得到，接下利用所得到的塊體數據進行可移動行的判別；關鍵塊體是巖體中的最薄弱環節，圍巖中巖石的冒落或滑動首先由關鍵塊開始，7從而導致相鄰塊體的失穩，並有可能引起連鎖反應，產生由圍巖冒落或者滑動導致的巖體工程失穩，因此，其核心是找出巖體臨空面處的關鍵塊體，以便對它們採取加固措施，保持巖體的穩定，由於將空間各組結構面和臨空面平移，使之通過坐標原點，則空間平面將構成以坐標原點為頂點的一系列稜錐，進行塊體可動性判斷時，首先判斷裂隙錐是否非空，然後加上臨空面形成塊體，判斷塊體錐是否為空，進行穩定性係數計算方法計算時採用莫爾庫侖強度理(Mohr-Coulomb)準則，主動力只考慮塊體自重，分兩種情況塊體可移動滿足三個條件具有出露面，幾何可移動，穩定係數小於規定值；(1)幾何可移動的判別如圖12所示；表示塊體所受驅動力，^表示塊體表面的方向矢量，指向塊體內部，i為多邊形在塊體內的編號，不包括出露面，幾何可動表示為^=&|^[520,^20}非空，;(S。^max(^S)，^)為塊體的運動方向；[ooes](2)穩定性的判別穩定性係數計算方法計算時採用莫爾庫侖強度理(Mohr-Coulomb)準則，主動力只考慮塊體自重，分兩種情況A.當塊體沿單面i滑動時，按下式計算穩定係數《=gcosa,.tan^("式中Q為塊體重量，a,為滑動面傾角，St為滑動面i的面積，&和^分別為滑動面上的內聚力和內摩擦角；B.當塊體沿雙面i和j滑動時，按下式計算穩定係數X二乂tan化+7V;ta，++②式中Q為塊體重量，Ci、P,和Cj、A分別為滑動面i和j上的內聚力及內摩擦角，Si和Sj分別為滑動面i和j的面積，a為滑動面i和j交線稜的傾角，Ni和Nj為作用在二滑動面上的法向力；並不能用解析的方法求得。的，所以考慮到用數值的方法求得，作一個可變的向量，令其指向下方即Z方向小於O，並按照一定的方向旋轉(間隔為1°)，依次試算出；；如圖13為判斷塊體可移動性的流程圖；6)、顯示圖像首先，按照塊體順序進行搜索，得到每個塊體的數據，進而對塊體中面進行搜索，最後對組成面的線段進行數據查詢，利用圖像處理的畫出線段，最後進行渲染；圖形處理處理開放式圖形庫(0penGL)作為本發明圖象接口，作為一個性能優越的圖形應用程式設計界面(API)，適用於廣泛的計算機環境，0penGL已成為目前的三維圖形開發標準，是從事三維圖形開發工作的技術人員通用的開發工具，為此，本模塊選用了該平臺作為圖象處理工具；本模塊中幾何頂點數據包括模型的頂點集、線集、多邊形集，這些數據經過流程圖的上部，包括運算器、逐個頂點操作等；圖像數據包括像素集、影像集、位圖集等，圖像像素數據的處理方式與幾何頂點數據的處理方式是不同的，但都經過光柵化、逐個片元處理直至把最後的光柵數據寫入幀緩衝器，在0penGL中的所有數據包括幾何頂點數據和像素數據都可以被存儲在顯示列表中或者立即可以得到處理，OpenGL中，顯示列表技術是一項重要的技術；OpenGL要求把所有的幾何圖形單元都用頂點來描述，這樣運算器和逐個頂點計算操作都可以針對每個頂點進行計算和操作，然後進行光柵化形成圖形碎片對於像素數據，像素操作結果被存儲在紋理組裝用的內存中，再類似於幾何頂點操作一樣光柵化形成圖形片元；整個流程操作的最後，圖形片元都要進行一系列的逐個片元操作，這樣最後的像素值送入幀緩衝器實現圖形的顯示；如圖14所示。本發明的優點通過該方法實現了工程巖體結構體快速搜索，關鍵塊體的識別，考慮工程實施過程中快速辨別新結構面，並能自動搜索關鍵塊體，統計出可移動塊體數量方法。圖1是本發明識別方法流程圖；圖2是本發明塊體的數據組成結構圖；圖3(a)是本發明多邊形數據結構圖；圖3(b)是本發明塊體數據結構圖；圖3(C)是本發明數據組織形式圖；圖4(a)是本發明巖體實際凹體結構圖；圖4(b)是本發明巖體分解為凸體結構圖；圖5(a)是本發明網格劃分方框流程圖；圖5(b)是本發明劃分網格後等效結構面圖；圖6是本發明網格劃分流程圖；圖7是本發明結構面切割塊體示意圖；圖8是本發明切割塊體流程圖；圖9是本發明開挖面處理圖；圖IO是本發明合併塊體流程圖；圖11(a)是本發明網格消除導致塊體合併示意11(b)是本發明消除網格前後塊體形態圖；圖12是本發明塊體可移動性判別示意圖；圖13是本發明塊體可移動性判別流程圖；圖14是本發明圖像顯示流程圖。具體實施例方式本實施例選取某高速公路隧道作為實例進行說明，具體實現過程如下1、工程概況該隧道軸線方向進口218°，出口21(T，斷面呈近似橢圓形，隧洞淨寬19.24m，高10.39m，開挖寬度為21.242m，軸線處開挖高度15.52m，隧道自樁號K344+670m開始至樁號K345+130m結束，全長520m，隧道縱坡1.6%，隧道位於R=5500m的平曲線上。n、野外地質調查野外地質調查從k345+360公路邊坡開始，向北300m，分兩個測段，測段南部長100m，北部長200m，中間衝溝處間斷約30m，邊坡巖體的巖性以灰巖為主，頁巖僅見於010m處，灰紫色，薄層狀，層面產狀58/45，此測段其餘10100m間皆為灰巖，在北測段146.5150m和165182m間為泥灰巖，層面產狀65/52，其餘皆為灰巖，自北測段200m處向北又為泥灰巖，並漸變為灰巖和泥灰巖交互沉積，根據現場地質調查，可以將該區域斷裂構造分為七組，各斷裂組的主要特徵如下第1組斷裂主要為層間節理和斷層，包括J7、J9、J22、J27等節理和f12、f13斷層，平均產狀66/56，屬壓扭性；第2組斷裂主要為近東西向的節理和斷層，包括J28、J32、J33、J37、J38、J44、J45等節理和f2、f8和f16等斷層，其中，如f8斷層，產狀175/88，壓扭性，斷層面較平直，鈣質泥質充填，破碎帶寬1020cm，延伸大於30m，f8斷層錯斷f7斷層，呈右旋扭動，即斷層上盤巖層向SWW方向錯動，由於被鈣泥質充填覆蓋，見不到原位移檫痕；第3組為北北東向的陡傾斷裂，包括J15、J34等節理和f1Q、fu和f18等斷層，平均產狀120/87，也包括產狀300/87的斷裂，如^。斷層，傾向120°，傾角86°，壓扭性，斷層面較平直，呈左旋扭動，即上盤巖層向NNE方向錯動，斷層面上見檫痕，與水平線呈20°;第4組為近南北向的斷層，包括J4。、J43、J47、J49等節理和f2。斷層，平均產狀280/48，此組斷裂在坡面上呈反傾切割巖體；第5組為北東東向的陡傾斷裂，包括4、f9、f14、f15斷層和Jn節理等，平均產狀148/82，它們多與邊坡平行，如f9斷層產狀140/89，斷層面充填鈣質(方解石)和泥質，表面粗糙，沿走向延伸13m，形成坡面；第6組為北東東向斷裂，傾向北西，對邊坡為反傾，包括J『J『，等節理，平均產狀310/60，為扭性;第7組分布在頁巖中，如衛節理，其產狀250/56，為張性節理；綜上所述，隧道區的七組斷裂構造及其特徵可以概括為表5.1表51隧道斷裂構造分組及主要特徵tableseeoriginaldocumentpage11III、模型的建立與關鍵塊體的搜索工程巖體三維空間結構建模與關鍵塊識別方法步驟如下—、工程巖體現場結構原始數據獲取1)該隧道工程巖體結構參數巷道、隧洞、斷面尺寸巖體模型參數為寬70m，高40m，長450m;2)工程巖體結構面信息，在工程中通過現場鑽孔勘探和隧道沿線結構面測量獲取結構面tableseeoriginaldocumentpage11數據初步分析後如下二、工程巖體原始數據的處理與提取，在現場結構調查的原始數據基礎上，採用蒙特卡羅法對不能測量和讀取的數據進行賦值，對現場的結構面進行統計，確定蒙特卡洛法所需要的參數tableseeoriginaldocumentpage12三、工程巖體模型構建1)巖體模型的建立本模型中建立較為簡單的凸體模型，不需要建立複雜的凹體模型，所以可以利用程序提供的對話框的固定格式輸入巖體尺寸即可；2)結構面信息，實際調查和產生的結構面的密度為0.008/m平均半徑為6.3m;四、塊體識別方法1):網格劃分，網格數因為巖體尺寸簡單，所以採用人工劃分，由於結構面密度較小，所以選取基元的尺寸為0.5mX0.5mX0.5m，對模型劃分1008萬個小基元；2):添加結構面，添加現場調查7組確定性結構面，18組非確定性結構面；3):開挖面尺寸，隧洞開挖面尺寸寬21.242m，高10.39m;4):消除網格，生成塊體，清除第1步人工生成虛擬網格1008萬個，生成塊體；5):判斷塊體的可移移動性，在本模型中，利用幾何參數和力學參數判斷可移動塊體4個；6):圖像顯示顯示關鍵塊體及信息tableseeoriginaldocumentpage12與現場測得的可移動塊體數據對比相吻合。權利要求一種工程巖體三維空間結構建模與關鍵塊識別方法，其特徵在於該識別方法步驟如下一、工程巖體現場結構原始數據獲取，數據包括工程巖體結構參數、工程巖體結構面信息和工程巖體結構狀態，其中工程巖體結構參數包括巷道、隧洞和斷面尺寸；工程巖體結構面信息和工程巖體結構狀態包括斷層參數和節理參數；二、工程巖體原始數據的處理與提取，按現場採集的原始數據，利用赤平投影和概率統計方法，整理出1)巖體裂隙方位，以及每一條裂隙所在的硐號、方位、層位、工程部位、風化帶信息，進行裂隙優勢方位分析；2)存儲工程普遍測網法獲得的每條裂隙的基本信息，包括每條裂隙對應的測點號、位置和產狀，進行繪製裂隙平面展布圖的坐標信息；3)存儲工程巖體測點的各條裂隙的產狀及相對測網原點的局部幾何坐標；三、構建工程巖體模型；1)、建立巖體模型判斷巖體幾何體被分成凸體和凹體兩類，建立模型時加以區別；凸體結構首先，把它分解為多邊形，依次添加多邊形；多邊形的建立是按同一順序即逆時針或順時針順序，依次輸入頂點坐標到文本文件中；凹體結構按照相鄰兩個塊體公共面重合的原則，把凹體分成多個凸體，並把重合面賦上虛擬面的標記值，將被分成的凸體按照凸體的建模規則建模；2)、建立數據結構A數據結構先確定塊體稜線的兩個端點，把放入線段存儲結構中，再把形成的線段放入多邊形存儲結構中，最後把多邊形放入塊體的數據結構中；B建立數據組織，存儲數據，形成數據鏈，採用程式語言中自帶的動態類對數據進行管理；3)、確定結構面信息方法包括確定性結構面、非確定性結構面和耦合結構面確定性結構面現場調查所得的結構面，包括地球物理勘探、數字圖象拍攝等手段直接採集的結構面數據，作為確定性結構面數據，方位、傾向、傾角、跡長等，輸入模型中；非確定性結構面通過概率統計的方法來模擬產生的結構面數據，根據現場鄰近工程、專家經驗，採用統計方法獲得結構面參數；耦合結構面模型中輸入工程中確定性結構面，將不能通過鑽孔、照相、物理方法獲取的數據，用概率統計的方法模擬產生結構面數據；四、對塊體進行識別；塊體識別方法把巖體模型劃分為有限個網格，然後在依次加入結構面，對現有的小塊體進行切割，最後在去掉網格，合併小的塊體，從而識別出有結構面切割得到的所有塊體；1)、網格劃分根據區域內結構面的密度，平均半徑進行網格劃分。首先確定區域的範圍，然後設置三組相互垂直的虛擬結構面對巖體模型進行切割，最後記錄下小單元的數據結構。2)、添加結構面網格換分後，加入實測結和擬合出來的結構面，對現有的小塊體進行再次切割；3)、開挖面處理，近似用平面代替曲面，曲面對塊體的切割就是對單元塊體的切割，求出單元體所有稜線與曲面的交點後依次連接各點，形成一個新的面，生成兩個新的塊體；4)、消除網格，生成塊體首先判斷相鄰兩個小單元面的交集是否與虛擬結構面有交集，如果有交集則合併兩個小塊，沒有交集就不必合併。按照此方法，依次對存在的所有小塊體進行處理，最後記錄下合併後塊體的數據；5)、判別塊體的可移移動性；首先，找出一個帶有兩個角度參數的試算向量，這兩個參數的初始值為90°和0°，依次試算，如果這個試算向量與塊體每個面的方向向量相乘都為正，則判斷次塊體幾何可移動，進一步利用莫爾庫侖強度理(Mohr-Coulomb)準則對其進行穩定性判斷，如果不穩定，就是關鍵塊體；6)、顯示圖像首先，按照塊體順序進行搜索，得到每個塊體的數據。進而對塊體中面進行搜索，最後對組成面的線段進行數據查詢，利用圖像處理的畫出線段，最後進行渲染。2.按權利要求1所述的工程巖體三維空間結構建模與關鍵塊識別方法，其特徵在於所述的步驟四中網格劃分按如下方法進行①切割過程裂隙面對塊體切割可分解為對其表面多邊形的切割，進而分解為裂隙面對直線的切割，實現裂隙與多邊形的切割，裂隙切割多邊形；有三種情況A如果多邊形頂點全在裂隙面的上方，此時把這個多邊形編號後放到上面的新塊體中；B如果多邊形頂點全在裂隙面的下方，此時把這個多邊形編號後放到下面的新塊體中；C如果多邊形頂點分布在裂隙面的兩側，用裂隙平面把這個多邊形切割成兩個新多邊形，編號後分別放入上下兩個塊體中；②採用這個裂隙切割其他與其相關的塊體，直到切割完與其相關的所有塊體，然後進入下一個裂隙的切割過程。全文摘要本發明涉及一種工程巖體三維空間結構建模與關鍵塊識別方法，屬於巖土工程、礦山開採、工程建築、水利工程、地下結構等巖體工程的實體結構建模和分析領域，步驟如下一、工程巖體現場結構原始數據獲取；二、工程巖體原始數據的處理與提取；三、構建工程巖體模型；四、塊體識別方法。本發明的優點通過該方法實現了工程巖體結構體快速搜索，關鍵塊體的識別，考慮工程實施過程中快速辨別新結構面，並能自動搜索關鍵塊體，統計出可移動塊體數量方法。文檔編號G06T17/00GK101739716SQ200910220460公開日2010年6月16日申請日期2009年12月4日優先權日2009年12月4日發明者朱萬成,楊勇,楊天鴻,王述紅申請人:東北大學