本發(fā)明涉及三維地質(zhì)建模領(lǐng)域，特別是涉及一種礦集區(qū)尺度區(qū)域三維地質(zhì)建模方法及系統(tǒng)。
背景技術(shù)：
：進(jìn)入21世紀(jì)以來，國際上礦產(chǎn)品需求旺盛、價格不斷攀升、勘查投資連創(chuàng)新高，礦業(yè)快速發(fā)展。在國內(nèi)，中國經(jīng)濟(jì)的持續(xù)、快速、健康發(fā)展，拉動能源和礦物原材料需求的快速增長，令世人矚目。中國礦業(yè)和地質(zhì)工作總體上與全球發(fā)展同步。礦業(yè)和地質(zhì)工作走出低谷，地勘經(jīng)濟(jì)快速發(fā)展，礦業(yè)行業(yè)利潤增長迅速；地質(zhì)找礦深受重視，勘查投人主體多元化，勘查經(jīng)費(fèi)快速增長；礦產(chǎn)勘查新發(fā)現(xiàn)不斷，儲量有新的增長，國內(nèi)礦產(chǎn)生產(chǎn)和供應(yīng)能力增強(qiáng)；地質(zhì)工作融入地方，與經(jīng)濟(jì)社會的聯(lián)系更加密切，服務(wù)領(lǐng)域不斷拓展。但是，與此同時，隨著我國工業(yè)化、城鎮(zhèn)化進(jìn)程的加快和經(jīng)濟(jì)社會的不斷發(fā)展，礦產(chǎn)資源匱乏的態(tài)勢在我國愈發(fā)顯現(xiàn)。在45種主要礦產(chǎn)中有11種國民經(jīng)濟(jì)支柱性礦產(chǎn)出現(xiàn)嚴(yán)重短缺趨勢，其中石油、鐵礦石、銅、鉀的對外依存度已超過50％。如不加強(qiáng)地質(zhì)勘查和加快轉(zhuǎn)變經(jīng)濟(jì)發(fā)展方式，礦產(chǎn)資源對經(jīng)濟(jì)發(fā)展的支撐力和承載力將面臨很大挑戰(zhàn)，成為制約經(jīng)濟(jì)可持續(xù)發(fā)展的關(guān)鍵因素。目前我國許多大中型礦產(chǎn)由于地質(zhì)儲量枯竭成為資源危機(jī)性礦山，但這些礦山許多是由于地質(zhì)勘查工作投入不夠，在其深部及其外圍資源潛力不清所致。為解決這些礦產(chǎn)資源危機(jī)問題，我國立專項(xiàng)進(jìn)行危機(jī)礦山外圍及深部礦產(chǎn)資源勘查工作。目前此項(xiàng)工作正在全國幾十個危機(jī)礦山全面實(shí)施。危機(jī)礦山外圍及深部礦產(chǎn)資源找礦的突破關(guān)鍵在于隱伏礦床三維立體定位預(yù)測的實(shí)現(xiàn)，成功的隱伏礦床三維立體定位預(yù)測，既依賴于新的理論和勘探技術(shù)方法，又依賴于使用先進(jìn)計算機(jī)及信息技術(shù)，對危機(jī)礦山地質(zhì)勘查工作獲取的多方法信息資料，進(jìn)行集成、綜合。危機(jī)礦產(chǎn)地質(zhì)勘探資料數(shù)據(jù)包括了從礦山勘探發(fā)現(xiàn)和開采后補(bǔ)充勘探等幾十年的資料，同時新的危機(jī)礦山找礦工作也會產(chǎn)生大量的基礎(chǔ)信息，這些信息資料涵蓋了地質(zhì)、地球物理、地球化學(xué)、遙感和深部勘探等多專題的數(shù)據(jù)；同時，危機(jī)礦山地質(zhì)勘探工作要求在三維立體空間進(jìn)行，即從空中、地表、地下進(jìn)行探測，形成全空間的立體三維數(shù)據(jù)。隨著危機(jī)礦山找礦工作的深入，對于已經(jīng)獲取的三維空間多元地質(zhì)勘查數(shù)據(jù)的有效管理，和在此基礎(chǔ)上的儲量估算處理，無疑將對危機(jī)礦山外圍及深部礦產(chǎn)資源找礦工作起到重要的推動作用。面對嚴(yán)峻的現(xiàn)實(shí)，提高礦產(chǎn)資源對經(jīng)濟(jì)社會發(fā)展的保障能力，探索創(chuàng)新地質(zhì)找礦新機(jī)制、實(shí)現(xiàn)地質(zhì)找礦新突破迫在眉睫。全國能源和重要礦產(chǎn)資源潛力分析結(jié)果表明，我國成礦地質(zhì)條件有利，礦產(chǎn)勘查程度較低，總體資源探明程度不足三分之一，還有很大的資源潛力。在具體找礦工作中，隨著地表礦、淺部礦及易識別礦的日益減少，找礦難度日益增大，找礦效果日益降低，進(jìn)入21世紀(jì)80年代以來，重點(diǎn)找尋深部隱伏礦、新類型礦及新領(lǐng)域礦已經(jīng)成為世界各國所關(guān)注的找礦方向，其中使用三維立體模型開展深部隱伏礦體(即，第二找礦空間)的尋找，已經(jīng)成為許多國家和地區(qū)找礦的主要對象，因此大比例尺(通常大于1:50000)的成礦預(yù)測的作用更加突出，它已經(jīng)成為礦床勘查工作的一個重要組(趙鵬大等，1992)?！暗诙傻V找礦空間”問題的提出，更強(qiáng)調(diào)的是在地殼深部有別于淺表部位的地質(zhì)構(gòu)造環(huán)境中尋找類似或有別于淺表部位已知礦床的新礦種或新類型礦床，因此，現(xiàn)代礦產(chǎn)勘查必須更新找礦觀念，依托多學(xué)科、新理論、新技術(shù)與新方法進(jìn)行聯(lián)合攻關(guān)。2011年6月和2012年9月，中國地質(zhì)調(diào)查局分別主辦了的“深部礦產(chǎn)遠(yuǎn)景調(diào)查研討會”和“深部礦產(chǎn)遠(yuǎn)景調(diào)查培訓(xùn)班”，其目的是貫徹落實(shí)國土資源部找礦新機(jī)制，實(shí)現(xiàn)地質(zhì)找礦戰(zhàn)略行動計劃目標(biāo)，促進(jìn)地質(zhì)調(diào)查從傳統(tǒng)向現(xiàn)代、從單一向綜合、從二維向三維轉(zhuǎn)變，推動新理論、新技術(shù)、新方法子在實(shí)際工作中的應(yīng)用。礦區(qū)隱伏礦體的資源預(yù)測要從二維平面預(yù)測向三維空間立體預(yù)測轉(zhuǎn)化，所涉及的三維技術(shù)方法，如：三維找礦信息的獲取、三維空間的可視化技術(shù)、基于三維空間的資源評價方法等，都是當(dāng)前學(xué)科前沿的重要課題，對它的研發(fā)和試驗(yàn)應(yīng)用具有重要的理論和實(shí)際意義。然而，當(dāng)前，針對礦集區(qū)尺度區(qū)域的地質(zhì)模型處于空白狀態(tài)，而且，由于礦集區(qū)尺度區(qū)域的面積遠(yuǎn)遠(yuǎn)比礦區(qū)的面積要大，這使得應(yīng)用于礦區(qū)資源勘探的建模方法不再適用于礦集區(qū)尺度區(qū)域，例如，在礦區(qū)我們需要鉆出多個鉆孔以輔助地質(zhì)模型的建立，而在礦集區(qū)尺度區(qū)域，由于礦集區(qū)尺度區(qū)域的面積很大，鉆孔的數(shù)量相較于礦區(qū)要多很多，成本也會急劇增加，可見，傳統(tǒng)的建模方法來對礦集區(qū)尺度區(qū)域地質(zhì)進(jìn)行建模并不理想和現(xiàn)實(shí)。技術(shù)實(shí)現(xiàn)要素：本發(fā)明的目的是提供一種三維地質(zhì)建模方法及系統(tǒng)，能夠在保障低成本的基礎(chǔ)上，對礦集區(qū)尺度區(qū)域的地質(zhì)進(jìn)行準(zhǔn)確的建模。為實(shí)現(xiàn)上述目的，本發(fā)明提供了如下方案：一種礦集區(qū)尺度區(qū)域三維地質(zhì)建模方法，所述方法包括：獲取地質(zhì)調(diào)查數(shù)據(jù)，所述地質(zhì)調(diào)查數(shù)據(jù)包括地層區(qū)特征數(shù)據(jù)、構(gòu)造線特征數(shù)據(jù)、產(chǎn)狀點(diǎn)特征數(shù)據(jù)、鉆孔特征數(shù)據(jù)、勘探線特征數(shù)據(jù)、槽探特征和坑道特征數(shù)據(jù)；定義地層與地層之間的時空關(guān)系、斷層之間的時空關(guān)系以及地層與斷層之間的時空關(guān)系；采用協(xié)同克里格算法對所述地質(zhì)調(diào)查數(shù)據(jù)和所述時空關(guān)系進(jìn)行隱式空間插值；根據(jù)插值結(jié)果構(gòu)建虛擬地質(zhì)剖面；對所述虛擬地質(zhì)剖面進(jìn)行編輯；根據(jù)所述虛擬地質(zhì)剖面確定礦集區(qū)三維地質(zhì)模型。可選的，所述根據(jù)插值結(jié)果構(gòu)建虛擬地質(zhì)剖面，具體包括：在構(gòu)建所述地質(zhì)剖面的過程中，根據(jù)鉆孔控制的地層巖性段和產(chǎn)狀信息，對生成的地質(zhì)體界線進(jìn)行修正；根據(jù)鉆孔中地層界限、巖性界線及礦體邊界對構(gòu)建的所述地質(zhì)剖面進(jìn)行優(yōu)化?？蛇x的，所述對所述虛擬地質(zhì)剖面進(jìn)行編輯，具體包括：對所述地質(zhì)剖面上的地質(zhì)體界線進(jìn)行檢查和編輯，使所述地質(zhì)體界線符合實(shí)際地質(zhì)情況。可選的，所述根據(jù)所述虛擬地質(zhì)剖面確定礦集區(qū)三維地質(zhì)模型，具體包括：獲取所述地質(zhì)剖面上的地質(zhì)體界線；以構(gòu)造分區(qū)為基礎(chǔ)，從所述地質(zhì)體界線中篩選出地質(zhì)體的邊界；將所述地質(zhì)體的邊界連接在一起，形成曲面。本發(fā)明還提供了一種礦集區(qū)尺度區(qū)域三維地質(zhì)建模系統(tǒng)，所述系統(tǒng)包括：地質(zhì)調(diào)查數(shù)據(jù)獲取單元，用于獲取地質(zhì)調(diào)查數(shù)據(jù)，所述地質(zhì)調(diào)查數(shù)據(jù)包括地層區(qū)特征數(shù)據(jù)、構(gòu)造線特征數(shù)據(jù)、產(chǎn)狀點(diǎn)特征數(shù)據(jù)、鉆孔特征數(shù)據(jù)、勘探線特征數(shù)據(jù)、槽探特征和坑道特征數(shù)據(jù)；時空關(guān)系定義單元，用于定義地層與地層之間的時空關(guān)系、斷層之間的時空關(guān)系以及地層與斷層之間的時空關(guān)系；插值單元，用于采用協(xié)同克里格算法對所述地質(zhì)調(diào)查數(shù)據(jù)和所述時空關(guān)系進(jìn)行隱式空間插值；虛擬地質(zhì)剖面構(gòu)建單元，用于根據(jù)插值結(jié)果構(gòu)建虛擬地質(zhì)剖面；地質(zhì)剖面編輯單元，用于對所述虛擬地質(zhì)剖面進(jìn)行編輯；三維地質(zhì)模型確定單元，用于根據(jù)所述虛擬地質(zhì)剖面確定礦集區(qū)三維地質(zhì)模型?？蛇x的，所述虛擬地質(zhì)剖面構(gòu)建單元，具體包括：修正子單元，用于在構(gòu)建所述地質(zhì)剖面的過程中，根據(jù)鉆孔控制的地層巖性段和產(chǎn)狀信息，對生成的地質(zhì)體界線進(jìn)行修正；優(yōu)化子單元，用于根據(jù)鉆孔中地層界限、巖性界線及礦體邊界對構(gòu)建的所述地質(zhì)剖面進(jìn)行優(yōu)化?？蛇x的，所述地質(zhì)剖面編輯單元，具體包括：地質(zhì)剖面編輯子單元，用于對所述地質(zhì)剖面上的地質(zhì)體界線進(jìn)行檢查和編輯，使所述地質(zhì)體界線符合實(shí)際地質(zhì)情況?？蛇x的，所述三維地質(zhì)模型確定單元，具體包括：地質(zhì)體界線獲取子單元，用于獲取所述地質(zhì)剖面上的地質(zhì)體界線；地質(zhì)體邊界篩選子單元，用于以構(gòu)造分區(qū)為基礎(chǔ)，從所述地質(zhì)體界線中篩選出地質(zhì)體的邊界；曲面生成子單元，用于將所述地質(zhì)體的邊界連接在一起，形成曲面。根據(jù)本發(fā)明提供的具體實(shí)施例，本發(fā)明公開了以下技術(shù)效果：本發(fā)明通過獲取礦集區(qū)尺度區(qū)域的地層區(qū)特征數(shù)據(jù)、構(gòu)造線特征數(shù)據(jù)、產(chǎn)狀點(diǎn)特征數(shù)據(jù)、鉆孔特征數(shù)據(jù)、勘探線特征數(shù)據(jù)、槽探特征和坑道特征數(shù)據(jù)，定義地層與地層之間的時空關(guān)系、斷層之間的時空關(guān)系以及地層與斷層之間的時空關(guān)系，并采用協(xié)同克里格算法對所述地質(zhì)調(diào)查數(shù)據(jù)和所述時空關(guān)系進(jìn)行隱式空間插值，根據(jù)插值結(jié)果構(gòu)建地質(zhì)剖面圖，以及對得到的剖面圖進(jìn)行編輯和修正，最終得到了礦集區(qū)尺度區(qū)域的三維地質(zhì)模型。實(shí)現(xiàn)了基于較少深部勘探工程控制點(diǎn)與地質(zhì)認(rèn)識約束的三維地質(zhì)體界面自動擬合，為大深度的找礦勘查工作提供有效的技術(shù)支撐，促進(jìn)地質(zhì)工作由二維向三維轉(zhuǎn)換。附圖說明為了更清楚地說明本發(fā)明實(shí)施例或現(xiàn)有技術(shù)中的技術(shù)方案，下面將對實(shí)施例中所需要使用的附圖作簡單地介紹，顯而易見地，下面描述中的附圖僅僅是本發(fā)明的一些實(shí)施例，對于本領(lǐng)域普通技術(shù)人員來講，在不付出創(chuàng)造性勞動性的前提下，還可以根據(jù)這些附圖獲得其他的附圖。圖1為本發(fā)明實(shí)施例礦集區(qū)尺度區(qū)域三維地質(zhì)建模方法流程圖；圖2為本發(fā)明實(shí)施例地質(zhì)填圖原始數(shù)據(jù)圖；圖3為本發(fā)明實(shí)施例1:50,000礦集區(qū)地質(zhì)填圖數(shù)據(jù)圖；圖4為本發(fā)明實(shí)施例地質(zhì)對象時空關(guān)系設(shè)置圖；圖5為本發(fā)明實(shí)施例剖面編輯圖；圖6為本發(fā)明實(shí)施例三維地質(zhì)模型圖；圖7為本發(fā)明實(shí)施例區(qū)域構(gòu)造分區(qū)圖；圖8為本發(fā)明實(shí)施例拓?fù)渲亟▓D；圖9為本發(fā)明實(shí)施例礦集區(qū)尺度區(qū)域三維地質(zhì)建模系統(tǒng)結(jié)構(gòu)圖。具體實(shí)施方式下面將結(jié)合本發(fā)明實(shí)施例中的附圖，對本發(fā)明實(shí)施例中的技術(shù)方案進(jìn)行清楚、完整地描述，顯然，所描述的實(shí)施例僅僅是本發(fā)明一部分實(shí)施例，而不是全部的實(shí)施例。基于本發(fā)明中的實(shí)施例，本領(lǐng)域普通技術(shù)人員在沒有做出創(chuàng)造性勞動前提下所獲得的所有其他實(shí)施例，都屬于本發(fā)明保護(hù)的范圍。本發(fā)明的目的是提供一種三維地質(zhì)建模方法及系統(tǒng)，能夠在保障低成本的基礎(chǔ)上，對礦集區(qū)尺度區(qū)域的地質(zhì)進(jìn)行準(zhǔn)確的建模。一種礦集區(qū)比例尺區(qū)域三維地質(zhì)建模方法，屬于地學(xué)空間建模領(lǐng)域。該地學(xué)空間建模的核心內(nèi)容是地質(zhì)屬性庫-時空關(guān)系模型-幾何算法-交互式編輯。該方法模型與裝置主要用于由1:50,000二維地質(zhì)圖出發(fā)，構(gòu)建礦集區(qū)尺度的三維地質(zhì)幾何模型。該方法主要內(nèi)容是，從地質(zhì)空間的固有特征分析入手，解決基于較少深部勘探工程控制點(diǎn)與地質(zhì)認(rèn)識約束的三維地質(zhì)體界面自動擬合關(guān)鍵技術(shù)問題，為大深度的找礦勘查工作提供有效的技術(shù)支撐，促進(jìn)地質(zhì)工作由二維向三維轉(zhuǎn)換?？傮w技術(shù)路線圖如圖1所示。首先，基于1:50,000二維地質(zhì)圖、少量鉆孔數(shù)據(jù)、實(shí)測剖面創(chuàng)建地質(zhì)屬性數(shù)據(jù)庫；其次，定義地層與地層之間、構(gòu)造之間以及地層與構(gòu)造之間的時空關(guān)系；第三，應(yīng)用協(xié)同克里格矩陣建立空間勢場；第四，進(jìn)行人機(jī)交互編輯，定義虛擬剖面并自動生成剖面圖；第五步，曲面連接與實(shí)體組裝。研究通用的、可更新的地質(zhì)體建模技術(shù)一直是地學(xué)空間建模中的重點(diǎn)和難點(diǎn)問題。這里“通用的”是指能夠滿足絕大多數(shù)地質(zhì)工作者需求的地質(zhì)體模擬方法。在這些需求中，當(dāng)然應(yīng)該包括那些缺少鉆孔和地球物理數(shù)據(jù)的用戶需求；可更新則是指地質(zhì)體的建模結(jié)果能夠隨著用戶需求的變化，快速升級，易于修改甚至是動態(tài)重構(gòu)。大比例尺空間地質(zhì)建模主要是面向礦集區(qū)的，因此，其所需的體建模方法不僅需要滿足已開采或有較多勘探工程的地質(zhì)體建模的需要，同時還要滿足地質(zhì)工作程度較低的礦集區(qū)的地質(zhì)體重構(gòu)的建模需求。與前者相比，后者具有研究區(qū)面積大、工程控制數(shù)據(jù)少、建模結(jié)果更加依賴于地質(zhì)專家的地質(zhì)解釋等特點(diǎn)。上述特性，決定了常用的基于勘探工程的地質(zhì)體三維建模方法，如輪廓線拼接、移動立方格等，并不能完全滿足礦集區(qū)尺度三維地質(zhì)體模擬的實(shí)際需要。為使本發(fā)明的上述目的、特征和優(yōu)點(diǎn)能夠更加明顯易懂，下面結(jié)合附圖和具體實(shí)施方式對本發(fā)明作進(jìn)一步詳細(xì)的說明。圖1為本發(fā)明實(shí)施例礦集區(qū)尺度區(qū)域三維地質(zhì)建模方法流程圖，如圖1所示，礦集區(qū)尺度區(qū)域三維地質(zhì)建模的具體步驟為：步驟101：獲取地質(zhì)調(diào)查數(shù)據(jù)，所述地質(zhì)調(diào)查數(shù)據(jù)包括地層區(qū)特征數(shù)據(jù)、構(gòu)造線特征數(shù)據(jù)、產(chǎn)狀點(diǎn)特征數(shù)據(jù)、鉆孔特征數(shù)據(jù)、勘探線特征數(shù)據(jù)、槽探特征和坑道特征數(shù)據(jù)；步驟102：定義地層與地層之間的時空關(guān)系、斷層之間的時空關(guān)系以及地層與斷層之間的時空關(guān)系；步驟103：采用協(xié)同克里格算法對所述地質(zhì)調(diào)查數(shù)據(jù)和所述時空關(guān)系進(jìn)行隱式空間插值；步驟104：根據(jù)插值結(jié)果構(gòu)建虛擬地質(zhì)剖面；步驟105：對所述虛擬地質(zhì)剖面進(jìn)行編輯；步驟106：根據(jù)所述虛擬地質(zhì)剖面確定礦集區(qū)三維地質(zhì)模型。其中，步驟101中，礦集區(qū)尺度地質(zhì)調(diào)查數(shù)據(jù)具體包括區(qū)域地質(zhì)圖(如1:5萬地形地質(zhì)圖)、區(qū)域地質(zhì)研究程度圖、構(gòu)造地質(zhì)圖、構(gòu)造分區(qū)圖、構(gòu)造綱要圖、區(qū)域礦產(chǎn)圖、巖性-巖相分布圖、巖相古地理圖、第四紀(jì)地質(zhì)圖等。其中，區(qū)域地質(zhì)圖可分為成果地質(zhì)圖和地質(zhì)填圖基本信息，成果地質(zhì)圖又包括平面地質(zhì)圖、地質(zhì)剖面圖和綜合柱狀圖。平面地質(zhì)圖提供了建模相關(guān)的地層、巖性、產(chǎn)狀、斷層等地質(zhì)要素以及地形、地貌特征、主要居民點(diǎn)等地理信息。地質(zhì)剖面圖描述了巖性、褶皺、斷層等地質(zhì)要素的空間展布形態(tài)以及產(chǎn)狀、地貌特征等。綜合柱狀圖主要描述平面圖區(qū)域內(nèi)的地層層序、厚度、巖性變化及接觸關(guān)系等。地質(zhì)填圖基本信息是指地質(zhì)填圖野外工作獲取的成果數(shù)據(jù)，通常以表格形式保存，如圖2所示(圖2為本發(fā)明實(shí)施例地質(zhì)填圖原始數(shù)據(jù)圖)。三維建?；緮?shù)據(jù)主要包括地質(zhì)調(diào)查三維數(shù)據(jù)(如產(chǎn)狀等)、地質(zhì)勘探數(shù)據(jù)(如鉆孔、槽探等)、地球物理解譯與三維反演數(shù)據(jù)。地質(zhì)調(diào)查三維數(shù)據(jù)需要包含地質(zhì)剖面、區(qū)域地層產(chǎn)狀、構(gòu)造產(chǎn)狀、構(gòu)造性質(zhì)、侵入巖體產(chǎn)狀、巖性變化等信息。地質(zhì)勘探工作使用鉆探、坑道、槽探等手段揭示近地表、地下地質(zhì)巖石分層、構(gòu)造、礦產(chǎn)等分布、產(chǎn)狀、質(zhì)量、數(shù)量等。地質(zhì)勘探產(chǎn)生的信息是通常所說的“硬數(shù)據(jù)”，是最寶貴的勘探資料。由于鉆探、坑道、槽探等手段都是在三維空間進(jìn)行的，數(shù)據(jù)具有真正三維特性。地球物理解譯與三維反演數(shù)據(jù)主要包括物探成果剖面圖，如等ab/2視電阻率剖面圖，物性-地質(zhì)剖面圖，平洞聲波測試剖面圖等，為隱伏地區(qū)三維地質(zhì)建模提供必要的信息支持。三維建?；緮?shù)據(jù)主要表現(xiàn)為圖層形式，如地層區(qū)文件、構(gòu)造線文件、產(chǎn)狀點(diǎn)文件、鉆孔文件以及由平面等值線數(shù)據(jù)生成的dem、tin等。其中，地層區(qū)文件的主要屬性為區(qū)編號、地層、巖性、顏色等；區(qū)域構(gòu)造線文件的主要屬性包括線編號、類型、傾向、傾角、級別、最大延伸深度等；產(chǎn)狀點(diǎn)文件的主要屬性包括傾向、傾角、左巖性、右?guī)r性等(表1)；鉆孔文件包含了鉆孔位置表、鉆孔形態(tài)表、采樣信息表中的所有屬性信息，來控制鉆孔的位置、形態(tài)及見礦鉆孔的礦體品位信息(表2、3、4)；探槽文件描述了導(dǎo)線位置分布、巖礦信息、層位、產(chǎn)狀、接觸關(guān)系等信息(表5)；坑道文件描述了坑道方位、長度、地質(zhì)情況等信息(表6)。表1地質(zhì)調(diào)查三維數(shù)據(jù)格式表2鉆孔位置表表3鉆孔形態(tài)表表4采樣信息表表5探槽記錄表表6坑道記錄表字段名稱數(shù)據(jù)類型長度小數(shù)位數(shù)字段狀態(tài)空值狀態(tài)坑道號文本型320必須可有空值編錄日期日期型80必須可有空值總長度雙精度型327必須可有空值比例尺文本型320必須可有空值制表人文本型200必須可有空值地質(zhì)描述可變文本500必須可有空值備注可變文本500可選可有空值表7坑道記錄表步驟102時空關(guān)系主要指定義構(gòu)造、地層邊界等的關(guān)系數(shù)據(jù)。根據(jù)地質(zhì)學(xué)原理，礦床的賦存是一個長期的地質(zhì)演化過程。在很多研究區(qū)內(nèi)，歷史上存在重大歷史事件和多期次的強(qiáng)烈構(gòu)造運(yùn)動等，導(dǎo)致構(gòu)造、地層等地質(zhì)對象之間存在復(fù)雜關(guān)系。這些關(guān)系的厘定需要通過礦床學(xué)、年代學(xué)等多種科學(xué)技術(shù)手段的組合。在模擬地層等曲面之前，定義構(gòu)造與構(gòu)造之間、地層與地層之間、以及構(gòu)造與地層之間的空間約束條件，如圖3所示(圖3為本發(fā)明實(shí)施例1:50,000礦集區(qū)地質(zhì)填圖數(shù)據(jù)圖)。步驟103基于協(xié)同克里格算法的快速隱式空間插值方法。應(yīng)用隱式空間插值算法，整合產(chǎn)狀、位置以及時空關(guān)系等地質(zhì)屬性數(shù)據(jù)。該方法模型的輸入數(shù)據(jù)集主要包括①.點(diǎn)集p＝{pj∈li，j＝1，2，…，n；i＝1，2，...，m}。其中，pk為點(diǎn)數(shù)據(jù)，li為地質(zhì)邊界線；②.產(chǎn)狀集合v＝{vk(xk，yk，zk)，k＝1，2，…，s}。其中vk為某地質(zhì)邊界的產(chǎn)狀；時空關(guān)系為集合r，如圖4所示(圖4為本發(fā)明實(shí)施例地質(zhì)對象時空關(guān)系設(shè)置圖)。協(xié)同克里格算法：傳統(tǒng)的地質(zhì)統(tǒng)計學(xué)方法如普通克立格法、泛克立格法等，只限于對單變量的研究與估算，并未能充分考慮空間尺度對區(qū)域化變量的影響。為解決此問題，matherongeorges在上世紀(jì)60年代提出了多元地質(zhì)統(tǒng)計學(xué)，其核心方法就是協(xié)同克立格。在此之后，goovaerts(1998)對其進(jìn)行了改進(jìn)。該方法思想是以協(xié)同區(qū)域化理論為基礎(chǔ)，以協(xié)變異函數(shù)為基本工具，研究定義于同一空間域內(nèi)，既有統(tǒng)計相關(guān)又具有空間位置相關(guān)的多元信息的空間結(jié)構(gòu)。在地質(zhì)找礦應(yīng)用中，一個具體的實(shí)例是，在研究多金屬礦床時，常遇到多個金屬元素伴生的現(xiàn)象，如銀，銅，鉛，鋅等(候景儒和黃競先，1990)。在數(shù)學(xué)上，我們可以將每一個化學(xué)元素看作為一個變量，則當(dāng)某一個變量的取樣不足時，單獨(dú)對該變量進(jìn)行普通克立格計算往往達(dá)不到要求的精度。但如果其他變量有較充足的取樣量時，研究該變量與其他變量間的空間關(guān)系，借助其它元素的樣品信息，應(yīng)用協(xié)同克立格法估算該元素的品位與空間關(guān)系。顯而易見，該方法相比于普通克立格法能明顯改進(jìn)估計精度。設(shè)有k個變量構(gòu)成協(xié)同區(qū)域化變量集合{zk(x)，k＝1，2，…，k}。假設(shè)它們是二階平穩(wěn)的，即一階矩和二階矩存在且平穩(wěn)：設(shè)k0為k＝1，2，…，k中某一特定值，即設(shè)為我們要估計的主變量，要估計的是中心點(diǎn)在x0處的承載上的平均值樣品有效數(shù)據(jù)集是而是確定在承載上的平均值，即要得到待估樣本點(diǎn)的無偏最優(yōu)估計值，必須滿足最小方差無偏估計的要求：第一，無偏條件，即第二，估計方差最小性，即滿足二階平穩(wěn)的協(xié)同區(qū)域化變量必須每一個區(qū)域化變量zk(x)的期望存在且平穩(wěn)，其互協(xié)方差函數(shù)也存在且平穩(wěn)即交叉變異函數(shù)的定義為：交叉協(xié)方差函數(shù)為ckk(h)＝e[zk(x+h)·zk(x)]公式5要求權(quán)重系數(shù)滿足協(xié)同克立格方程組，即：方程組未知數(shù)為個權(quán)系數(shù)及k個拉格朗日參數(shù)μk，共有個未知數(shù)，同時也有個線性方程構(gòu)成這個方程組。算法的矩陣形式為：kλ＝m公式7其中：基于上述協(xié)同克里格算法，由二維地質(zhì)圖出發(fā)，構(gòu)建區(qū)域大比例尺三維地質(zhì)體，需要回答三個問題。即建模方法，如何在稀疏數(shù)據(jù)場解譯并推斷出合理的地質(zhì)曲面；如何解決整合主要構(gòu)造曲面的問題；如何檢驗(yàn)復(fù)雜地質(zhì)體模型質(zhì)量。本方法模型包含以下三個前提假設(shè)。即①.地質(zhì)曲面是界定兩個目標(biāo)地質(zhì)體的幾何對象；②.產(chǎn)狀數(shù)據(jù)決定了地層等地質(zhì)曲面的空間展布；③.在未定義約束關(guān)系的情況下，地質(zhì)曲面的延展是無限的。應(yīng)用該模型所需要的輸入數(shù)據(jù)包括①.地質(zhì)邊界的出露(至少一部分)是已知的；②.每個地質(zhì)邊界至少采集一個產(chǎn)狀點(diǎn)。由該方法模型所考慮的前提假設(shè)、輸入數(shù)據(jù)以及適用條件來看，可以解決如下地質(zhì)曲面的重構(gòu)問題：①.沉積層系的模擬，即不同年代的地層是近似平行的；②.侵入巖體的模擬，即存在已知出露以及地質(zhì)剖面作為約束條件；③.變質(zhì)巖體的模擬，即存在已知邊界線及其產(chǎn)狀作為約束條件。該方法模型的輸入數(shù)據(jù)集主要包括①.點(diǎn)集p＝{pj∈li，j＝1，2，…，n；i＝1，2，...，m}。其中，pk為點(diǎn)數(shù)據(jù)，li為地質(zhì)邊界線；②.產(chǎn)狀集合v＝{vk(xk，yk，zk)，k＝1，2，…，s}。其中vk為某地質(zhì)邊界的產(chǎn)狀。步驟104構(gòu)建虛擬(圖切)地質(zhì)剖面。要獲得地質(zhì)體曲面，首先要創(chuàng)建虛擬地質(zhì)剖面。虛擬地質(zhì)剖面中地層界線的模擬是通過對地層區(qū)產(chǎn)狀的綜合分析得到的結(jié)果。在創(chuàng)建圖切剖面的過程中，可以加入鉆孔進(jìn)行綜合計算，在自動生成剖面的過程中，計算機(jī)會根據(jù)鉆孔控制的地層巖性段和產(chǎn)狀信息，對生成的地質(zhì)體界線進(jìn)行修正，以便于更符合實(shí)際地質(zhì)情況。同時，鉆孔的選擇主要考慮孔深和巖性，地層方面需打穿含礦層，選擇盡量深的鉆孔；在巖性方面，需盡量保持巖性完整，即鉆孔盡量包括圍巖到礦，再到圍巖多個界線。在定義勘探線的時候，就需要把參與模擬的鉆孔依次關(guān)聯(lián)到離鉆孔最近的剖面上，圖切剖面自動擬合算法會根據(jù)鉆孔中地層、巖性界線及礦體邊界對模擬結(jié)果進(jìn)行優(yōu)化處理。步驟105依次對剖面上的地質(zhì)體界線進(jìn)行檢查和編輯，使其符合實(shí)際地質(zhì)情況。重點(diǎn)檢查地質(zhì)體界線是否相交，地質(zhì)體界線的形態(tài)是否符合地質(zhì)規(guī)律，地質(zhì)體界線與地表的交點(diǎn)與實(shí)際出露情況是否一致，地質(zhì)體界線是否缺失，構(gòu)造線的錯動距離是否合理。針對檢查結(jié)果，對剖面上的地質(zhì)體界線進(jìn)行編輯，依據(jù)已有三維地質(zhì)概念模型中層序及地層厚度，將缺失的地質(zhì)體界線補(bǔ)充完整，如圖5所示(圖5為本發(fā)明實(shí)施例剖面編輯圖)。步驟106在獲取地質(zhì)體界線的基礎(chǔ)上，以構(gòu)造分區(qū)為基礎(chǔ)，逐一篩選地質(zhì)體的邊界，將地質(zhì)體的邊界線連接成曲面。對于地層來說，首先篩選地層的上下界面，分別連接成兩個曲面，進(jìn)而生成其他側(cè)面。曲面的生成是基于三角網(wǎng)的插值實(shí)現(xiàn)的，只對曲面整體進(jìn)行編輯和修改，在某些特殊情況無法解決具體問題，如曲面中某三角形位置不對，或出現(xiàn)了三角形交叉，三角形重疊等現(xiàn)象。這就需要對于基礎(chǔ)三角形的一系列編輯功能，如刪除三角形，添加三角形，修改三角形頂點(diǎn)坐標(biāo)，修改三角形名稱，修改三角形顏色等。生成的完整曲面還需再次進(jìn)行拓?fù)錂z查，查找在編輯過程中產(chǎn)生的懸掛點(diǎn)、懸掛邊、相交線、未封閉線等問題。根據(jù)計算機(jī)生成的拓?fù)錂z查結(jié)果，進(jìn)行三維空間編輯，刪除懸掛點(diǎn)及懸掛邊，將相交線、未封閉線等通過打斷線、連接線、結(jié)點(diǎn)平差等具體功能解決其拓?fù)鋯栴}。封閉曲面通過計算機(jī)拓?fù)錂z查后，進(jìn)行拓?fù)渲亟?，可形成三維地質(zhì)體。生成的三維地質(zhì)體應(yīng)根據(jù)需要修改其顏色，使其與地質(zhì)圖中顏色近似，便于三維空間分析。使用ascii碼格式文件保存礦集區(qū)三維地質(zhì)模型，如圖6所示(圖6為本發(fā)明實(shí)施例三維地質(zhì)模型圖)。作為本發(fā)明的又一個實(shí)施例，首先，數(shù)據(jù)收集與整理，花垣-漁塘鉛鋅礦自20世紀(jì)以來，進(jìn)行了多次勘查工作，取得了海量的勘查數(shù)據(jù)，如鉆孔數(shù)據(jù)、探槽數(shù)據(jù)、坑道工程數(shù)據(jù)、物化探數(shù)據(jù)、相關(guān)勘查報告等。研究區(qū)通過1：5萬的區(qū)調(diào)和礦調(diào)工作，積累了花垣縣幅和麻栗場幅等基礎(chǔ)地質(zhì)資料，為重點(diǎn)區(qū)域的區(qū)域三維地質(zhì)建模提供了數(shù)據(jù)基礎(chǔ)。通過對相關(guān)資料進(jìn)行收集、整理發(fā)現(xiàn)，隨著數(shù)字化建設(shè)的不斷深入，近年來的相關(guān)資料和圖件基本采用電子文檔和計算機(jī)制圖，而早期的勘查資料大多為紙質(zhì)資料。電子版資料需要按照研究區(qū)三維空間數(shù)據(jù)庫的格式要求進(jìn)行格式轉(zhuǎn)換，或重新錄入，經(jīng)檢查通過后進(jìn)入數(shù)據(jù)庫。而紙質(zhì)文檔和圖件則首先要進(jìn)行數(shù)字化，信息提取與重新錄入，經(jīng)格式檢查后進(jìn)入數(shù)據(jù)庫。通?；A(chǔ)數(shù)據(jù)主要來源于地質(zhì)、物探、化探、遙感、重砂等方面。研究區(qū)收到的勘查數(shù)據(jù)中收集圖件具體包括花垣幅1:5萬地形地質(zhì)圖、麻栗場幅1:5萬地形地質(zhì)圖、花垣-鳳凰整裝勘查區(qū)地形地質(zhì)圖(1:20萬)、湖南省鉛鋅礦成礦規(guī)律圖(1:500000)、湖南省鉛鋅礦成礦要素圖(1:500000)、湖南省鉛鋅礦預(yù)測成果圖(1:500000)、湖南省龍山-鳳凰預(yù)測工作區(qū)李梅式層控?zé)嵋盒豌U鋅礦成礦要素圖(1:100000)、湖南省龍山-鳳凰預(yù)測工作區(qū)李梅式層控?zé)嵋盒豌U鋅礦預(yù)測要素圖(1:100000)、湖南省龍山-鳳凰預(yù)測工作區(qū)李梅式層控?zé)嵋盒豌U鋅礦預(yù)測成果圖(1:100000)、湖南省花垣縣漁塘寨礦區(qū)鉛鋅礦成礦要素圖(1:10000)、湖南省花垣縣漁塘寨礦區(qū)鉛鋅礦預(yù)測要素圖(1:10000)、花垣-鳳凰整裝勘查區(qū)最小預(yù)測區(qū)分布圖(1:200000)、花垣-鳳凰整裝勘查區(qū)化探異常圖(1:200000)、花垣-鳳凰整裝勘查區(qū)物探異常圖(1:200000)、花垣-鳳凰整裝勘查區(qū)遙感解譯成果圖(1:200000)等、白巖礦區(qū)鉛鋅礦勘探線剖面圖(1:5000)、白巖礦區(qū)鉛鋅礦鉆孔柱狀圖、清水塘礦區(qū)白巖礦區(qū)鉛鋅礦勘探線剖面圖(1:5000)、清水塘礦區(qū)鉛鋅礦鉆孔柱狀圖、三塘礦區(qū)白巖礦區(qū)鉛鋅礦勘探線剖面圖(1:5000)、三塘礦區(qū)鉛鋅礦鉆孔柱狀圖、豬圈坪礦區(qū)鉛鋅礦鉆孔柱狀圖、大腦坡礦區(qū)鉛鋅礦勘探線剖面圖、楊家寨礦區(qū)鉛鋅礦勘探線剖面圖、野外實(shí)測剖面圖、野外素描圖等。湖南湘西北鉛鋅礦位于揚(yáng)子陸塊東南緣與江南地軸過渡部位的八面山褶皺帶內(nèi)(陳明輝等,2011；羅衛(wèi)等,2009；孫玉嫻等,1985)。在經(jīng)歷了武陵、加里東、海西、印支、燕山、喜馬拉雅期的構(gòu)造運(yùn)動之后，在區(qū)域上形成了總體以北東向褶皺變形和深大斷裂為主的構(gòu)造格架(楊紹祥，勞可通,2007a；張攀華,1984)。深大斷裂主要包括花垣—張家界斷裂、吉首—古丈斷裂、麻栗場斷裂，呈北北東—北東—北東東—近東西方向的帚狀弧形展布形態(tài)，北東側(cè)朝著保靖、張家界方向收斂，南西側(cè)朝著貴州境內(nèi)發(fā)散(付勝云等,2006；楊紹祥和勞可通,2007b)。研究區(qū)位于湘西北帚狀弧形構(gòu)造帶的外側(cè)(北西側(cè))，是一個在揚(yáng)子旋回基礎(chǔ)上發(fā)展起來的長期沉降單元(羅衛(wèi),2009)。區(qū)內(nèi)主要地層分布有元古界震旦系、古生界寒武系、奧陶系、志留系、泥盆系。其中下古生界尤為發(fā)育(楊紹祥和勞可通,2007b)，地層之間沒有明顯的角度不整合接觸現(xiàn)象。區(qū)內(nèi)主要含礦層為寒武系下統(tǒng)清虛洞組中段(∈2q2)的厚層灰?guī)r，平均厚度約為200m。清虛洞組巖石類型以灰?guī)r為主，自下而上由灰?guī)r向白云巖過渡(周云等,2011)?？刂蒲芯繀^(qū)鉛鋅礦的巖相為早寒武世晚期碳酸鹽巖臺地邊緣淺灘-藻礁相(付勝云,2011)，藻礁分布情況與鉛鋅礦化帶的分布基本吻合。藻礁越發(fā)育，藻礁灰?guī)r越厚，則礦化越有利。礦體以淺黃色閃鋅礦為主，呈斑塊狀、浸染狀充填于藻礁灰?guī)r的空隙中，似層狀多層產(chǎn)出，脈石礦物為方解石、重晶石、螢石等(劉文均等,1994)?？偨Y(jié)并提出的地質(zhì)屬性庫-時空關(guān)系-幾何算法-交互式編輯模型出發(fā)，開展與mvt鉛鋅礦床成礦預(yù)測相關(guān)的三維區(qū)域地質(zhì)建模驗(yàn)證工作。(1)區(qū)域鉆孔模型在區(qū)域地質(zhì)建模驗(yàn)證工作中，仍然會包含少量三維鉆孔模型。與傳統(tǒng)方法相似，基于描述鉆孔的三個數(shù)據(jù)表，鉆孔位置表、鉆孔形態(tài)表和采樣信息表。它們能夠精確地描述了鉆孔上的空間采樣點(diǎn)在三維空間中的形態(tài)和位置關(guān)系。通過對三個數(shù)據(jù)表的有效設(shè)置，就可以建立三維鉆孔模型。(2)構(gòu)建地表地質(zhì)體表面模型將平面地質(zhì)圖融合到地形等高線生成的數(shù)字高程模型(dem)中，形成全地表地質(zhì)體表面，即可表達(dá)出露地質(zhì)體地表分布情況。同時，也便于驗(yàn)證建模結(jié)果與原始地質(zhì)信息的一致性。首先，將矢量地質(zhì)圖中地層、巖體等圖層整理為一個區(qū)圖層，進(jìn)行必要的編輯，如根據(jù)具體需要刪除河流等要素，確保每個區(qū)都是完整封閉的。接下來利用地形等高線生成數(shù)字高程模型，并將第一步中整理好的區(qū)圖層融合到dem中，生成全地表地質(zhì)體表面。(3)地質(zhì)體界線三維模擬①.構(gòu)造分區(qū)首先統(tǒng)計研究區(qū)內(nèi)構(gòu)造的性質(zhì)、產(chǎn)狀，及區(qū)域內(nèi)地層的產(chǎn)狀信息，依據(jù)三維地質(zhì)格架模型，理解整個區(qū)域內(nèi)的構(gòu)造條件及地層演化特征。然后根據(jù)研究區(qū)內(nèi)主要構(gòu)造的性質(zhì)、產(chǎn)狀及對地層、成礦的影響，以茶垌-花垣-張家界大斷裂為邊界，將產(chǎn)狀相近的地質(zhì)體構(gòu)造分區(qū)，將整個研究區(qū)分為兩個大的區(qū)塊，一個為大斷裂北西側(cè)，總體傾向?yàn)楸北蔽?北西向；另一個區(qū)塊為大斷裂南東側(cè)，總體傾向?yàn)槟蠔|-南東東向。本項(xiàng)目分別對兩個區(qū)塊進(jìn)行地質(zhì)體曲面的三維模擬，如圖7所示(圖7為本發(fā)明實(shí)施例區(qū)域構(gòu)造分區(qū)圖)。②.虛擬剖面的模擬要獲得地質(zhì)體曲面，首先要創(chuàng)建圖切剖面。圖切剖面中地層界線的模擬是通過對地層區(qū)產(chǎn)狀的綜合分析得到的結(jié)果。在創(chuàng)建圖切剖面的過程中，可以加入鉆孔進(jìn)行綜合計算，在自動生成剖面的過程中，計算機(jī)會根據(jù)鉆孔控制的地層巖性段和產(chǎn)狀信息，對生成的地質(zhì)體界線進(jìn)行修正，以便于更符合實(shí)際地質(zhì)情況。研究區(qū)主干勘探線方向?yàn)楸蔽?南東方向，近似于圖幅的對角線方向?？碧骄€剖面按主干剖面方向設(shè)置，共設(shè)置了與主干剖面方向平行的54條剖面。在李梅礦區(qū)共418個鉆孔，投影到1：5萬的花垣幅之后，鉆孔分布過于集中，因此，考慮到建模需要，在已知鉆孔中選擇了8個鉆孔參與圖切剖面中地質(zhì)界線模擬。鉆孔的選擇主要考慮孔深和巖性，地層方面需打穿含礦層，選擇盡量深的鉆孔；在巖性方面，需盡量保持巖性完整，即鉆孔盡量包括圍巖到礦，再到圍巖多個界線。在定義勘探線的時候，就需要把參與模擬的鉆孔依次關(guān)聯(lián)到離鉆孔最近的剖面上，圖切剖面自動擬合算法會根據(jù)鉆孔中地層、巖性界線及礦體邊界對模擬結(jié)果進(jìn)行優(yōu)化處理。③.人機(jī)交互編輯依次對剖面上的地質(zhì)體界線進(jìn)行檢查和編輯，使其符合實(shí)際地質(zhì)情況。重點(diǎn)檢查地質(zhì)體界線是否相交，地質(zhì)體界線的形態(tài)是否符合地質(zhì)規(guī)律，地質(zhì)體界線與地表的交點(diǎn)與實(shí)際出露情況是否一致，地質(zhì)體界線是否缺失，構(gòu)造線的錯動距離是否合理。針對檢查結(jié)果，對剖面上的地質(zhì)體界線進(jìn)行編輯，依據(jù)三維地質(zhì)格架模型中層序及地層厚度，將缺失的地質(zhì)體界線補(bǔ)充完整。結(jié)合地球物理解譯成果對剖面上的地質(zhì)體界線進(jìn)行修正，特別是對于隱伏地區(qū)的地質(zhì)體界線的推斷時，要充分考慮地球物理解譯的推斷結(jié)果，對于巖性界面及巖體邊界進(jìn)行合理修正。(4)曲面連接在獲取地質(zhì)體界線的基礎(chǔ)上，以構(gòu)造分區(qū)為基礎(chǔ)，逐一篩選地質(zhì)體的邊界，將地質(zhì)體的邊界線連接成曲面。對于地層來說，首先篩選地層的上下界面，分別連接成兩個曲面，進(jìn)而生成其他側(cè)面。曲面的生成是基于三角網(wǎng)的插值實(shí)現(xiàn)的，只對曲面整體進(jìn)行編輯和修改，在某些特殊情況無法解決具體問題，如曲面中某三角形位置不對，或出現(xiàn)了三角形交叉，三角形重疊等現(xiàn)象。這就需要對于基礎(chǔ)三角形的一系列編輯功能，如刪除三角形，添加三角形，修改三角形頂點(diǎn)坐標(biāo)，修改三角形名稱，修改三角形顏色等。(5)曲面成體生成的完整曲面還需再次進(jìn)行拓?fù)錂z查，查找在編輯過程中產(chǎn)生的懸掛點(diǎn)、懸掛邊、相交線、未封閉線等問題。根據(jù)計算機(jī)生成的拓?fù)錂z查結(jié)果，進(jìn)行三維空間編輯，刪除懸掛點(diǎn)及懸掛邊，將相交線、未封閉線等通過打斷線、連接線、結(jié)點(diǎn)平差等具體功能解決其拓?fù)鋯栴}。封閉曲面通過計算機(jī)拓?fù)錂z查后，進(jìn)行拓?fù)渲亟ǎ尚纬扇S地質(zhì)體。生成的三維地質(zhì)體應(yīng)根據(jù)需要修改其顏色，使其與地質(zhì)圖中顏色近似，便于三維空間分析，如圖8所示(圖8為本發(fā)明實(shí)施例拓?fù)渲亟▓D)。本發(fā)明提供的方法基于少量工程控制點(diǎn)、地形地質(zhì)圖、實(shí)測剖面、產(chǎn)狀、地質(zhì)體界面露頭線以及構(gòu)造等地質(zhì)認(rèn)識的地質(zhì)體界面自動模擬。由于采用隱式空間插值作為三維地質(zhì)建模算法的核心機(jī)，因此本發(fā)明具有建模速度快速、動態(tài)可更新、易于擴(kuò)展性等特性以及在避免大量重復(fù)勞動等，相較于目前流行的單純的人機(jī)交互建模技術(shù)，具有明顯優(yōu)勢。本發(fā)明還提供了一種礦集區(qū)尺度區(qū)域三維地質(zhì)建模系統(tǒng)，圖9為本發(fā)明實(shí)施例礦集區(qū)尺度區(qū)域三維地質(zhì)建模系統(tǒng)結(jié)構(gòu)圖，如圖9所示，建模系統(tǒng)包括：地質(zhì)調(diào)查數(shù)據(jù)獲取單元901，用于獲取地質(zhì)調(diào)查數(shù)據(jù)，所述地質(zhì)調(diào)查數(shù)據(jù)包括地層區(qū)特征數(shù)據(jù)、構(gòu)造線特征數(shù)據(jù)、產(chǎn)狀點(diǎn)特征數(shù)據(jù)、鉆孔特征數(shù)據(jù)、勘探線特征數(shù)據(jù)、槽探特征和坑道特征數(shù)據(jù)；時空關(guān)系定義單元902，用于定義地層與地層之間的時空關(guān)系、斷層之間的時空關(guān)系以及地層與斷層之間的時空關(guān)系；插值單元903，用于采用協(xié)同克里格算法對所述地質(zhì)調(diào)查數(shù)據(jù)和所述時空關(guān)系進(jìn)行隱式空間插值；虛擬地質(zhì)剖面構(gòu)建單元904，用于根據(jù)插值結(jié)果構(gòu)建虛擬地質(zhì)剖面；地質(zhì)剖面編輯單元905，用于對所述虛擬地質(zhì)剖面進(jìn)行編輯；三維地質(zhì)模型確定單元906，用于根據(jù)所述虛擬地質(zhì)剖面確定礦集區(qū)三維地質(zhì)模型。其中，所述虛擬地質(zhì)剖面構(gòu)建單元904，具體包括：修正子單元，用于在構(gòu)建所述地質(zhì)剖面的過程中，根據(jù)鉆孔控制的地層巖性段和產(chǎn)狀信息，對生成的地質(zhì)體界線進(jìn)行修正；優(yōu)化子單元，用于根據(jù)鉆孔中地層界限、巖性界線及礦體邊界對構(gòu)建的所述地質(zhì)剖面進(jìn)行優(yōu)化。所述地質(zhì)剖面編輯單元905，具體包括：地質(zhì)剖面編輯子單元，用于對所述地質(zhì)剖面上的地質(zhì)體界線進(jìn)行檢查和編輯，使所述地質(zhì)體界線符合實(shí)際地質(zhì)情況。所述三維地質(zhì)模型確定單元906，具體包括：地質(zhì)體界線獲取子單元，用于獲取所述地質(zhì)剖面上的地質(zhì)體界線；地質(zhì)體邊界篩選子單元，用于以構(gòu)造分區(qū)為基礎(chǔ)，從所述地質(zhì)體界線中篩選出地質(zhì)體的邊界；曲面生成子單元，用于將所述地質(zhì)體的邊界連接在一起，形成曲面。本發(fā)明提供的系統(tǒng)能夠基于少量工程控制點(diǎn)、地形地質(zhì)圖、實(shí)測剖面、產(chǎn)狀、地質(zhì)體界面露頭線以及構(gòu)造等地質(zhì)認(rèn)識的地質(zhì)體界面自動模擬。由于采用隱式空間插值作為三維地質(zhì)建模算法的核心機(jī)，因此本發(fā)明具有建模速度快速、動態(tài)可更新、易于擴(kuò)展性等特性以及在避免大量重復(fù)勞動等，相較于目前流行的單純的人機(jī)交互建模技術(shù)，具有明顯優(yōu)勢。本說明書中各個實(shí)施例采用遞進(jìn)的方式描述，每個實(shí)施例重點(diǎn)說明的都是與其他實(shí)施例的不同之處，各個實(shí)施例之間相同相似部分互相參見即可。對于實(shí)施例公開的系統(tǒng)而言，由于其與實(shí)施例公開的方法相對應(yīng)，所以描述的比較簡單，相關(guān)之處參見方法部分說明即可。本文中應(yīng)用了具體個例對本發(fā)明的原理及實(shí)施方式進(jìn)行了闡述，以上實(shí)施例的說明只是用于幫助理解本發(fā)明的方法及其核心思想；同時，對于本領(lǐng)域的一般技術(shù)人員，依據(jù)本發(fā)明的思想，在具體實(shí)施方式及應(yīng)用范圍上均會有改變之處。綜上所述，本說明書內(nèi)容不應(yīng)理解為對本發(fā)明的限制。當(dāng)前第1頁12