[發(fā)明專利]一種礦體三維品位模型局部更新方法

在審

• 申請?zhí)枺?/td>: 202010540865.2
• 申請日: 2020-06-15
• 公開（公告）號: CN111739165A
• 公開（公告）日: 2020-10-02
• 發(fā)明（設(shè)計）人: 王海龍;徐連生;郭甲騰;柳小波;王江梅
• 申請（專利權(quán)）人: 鞍鋼集團礦業(yè)有限公司
• 主分類號: G06T17/10
• 分類號: G06T17/10;G06F17/18;G06F16/23
• 代理公司: 鞍山貝爾專利代理有限公司 21223
• 代理人: 王常謙

本發(fā)明涉及一種礦體三維品位模型局部更新方法，包括以下步驟：S1：初始礦體三維品位模型計算，S2、地統(tǒng)計學(xué)計算，S3、礦體三維品位模型估值，S4、更新范圍判斷以及S5、礦體三維品位模型更新。本發(fā)明的優(yōu)點是：在地質(zhì)知識的驅(qū)動下，基于地統(tǒng)計學(xué)方法和搜索橢球?qū)崿F(xiàn)礦體三維品位模型局部更新，能夠有效解決礦體三維品位模型更新困難問題。

技術(shù)領(lǐng)域

本發(fā)明涉及三維礦體建模技術(shù)領(lǐng)域，尤其涉及一種礦體三維品位模型局部更新方法。

背景技術(shù)

在地質(zhì)礦產(chǎn)勘探及其信息化過程中，三維礦體建模是一個重要的研究方向，在勘探開發(fā)和資源管理中具有重要的應(yīng)用意義。目前在地理空間中礦體三維實體建模方法日趨成熟，但精細礦體三維品位模型的構(gòu)建依然是一個亟待解決的難題。

礦體三維品位模型大多是通過對鉆孔或剖面的分層界面數(shù)據(jù)進行插值或拓撲推理得到的，缺少地質(zhì)知識的融入，屬于純數(shù)學(xué)的空間格架建模方式，易造成所構(gòu)建礦體三維品位模型與實際的地質(zhì)內(nèi)容脫節(jié)。這樣的礦體三維品位模型更新困難，難以支持地礦開采和生產(chǎn)需求。

目前，大多是通過多種軟件聯(lián)合操作進行礦體三維品位模型的修正以實現(xiàn)品位模型更新，如Surpac、3DMine等。此種方法難以滿足礦山增量數(shù)據(jù)(日常生產(chǎn)樣品)導(dǎo)致品位模型局部更新的現(xiàn)勢性要求。而且，手動交互式的礦體三維品位模型更新勢必不能實時地反應(yīng)出礦體的形態(tài)變化，不能將礦體三維品位模型更新與礦山生產(chǎn)實踐相結(jié)合。除此之外，現(xiàn)有的方法都是通過人工指定礦體三維品位模型更新范圍，無法依據(jù)變化數(shù)據(jù)的各向異性合理確定礦體三維品位模型更新范圍。因此，現(xiàn)有方法嚴重影響礦山生產(chǎn)的進度、資源回收利用率和經(jīng)濟收益。

發(fā)明內(nèi)容

本發(fā)明針對現(xiàn)有礦體三維品位模型更新方法人工交互量巨大、耗時、效率低等問題，提供了一種礦體三維品位模型局部更新方法，在地質(zhì)知識的驅(qū)動下，基于地統(tǒng)計學(xué)方法和搜索橢球?qū)崿F(xiàn)礦體三維品位模型局部更新，能夠有效解決礦體三維品位模型更新困難問題。

本發(fā)明提供一種礦體三維品位模型局部更新方法，包括以下步驟：

步驟1、初始礦體三維品位模型計算：從鉆孔數(shù)據(jù)中，依據(jù)品位信息、礦體圈定信息及勘探線信息實現(xiàn)礦體三維品位模型初始計算；

步驟2、地統(tǒng)計學(xué)計算：依據(jù)鉆孔數(shù)據(jù)及勘探線信息，建立地統(tǒng)計學(xué)所要用的概率密度函數(shù)；

步驟3、品位估值：通過步驟1獲得的初始礦體三維品位模型及步驟2獲得的概率密度函數(shù)，獲得礦體三維品位模型估值；

步驟4、更新范圍的判斷：根據(jù)新增鉆孔數(shù)據(jù)、減少鉆孔數(shù)據(jù)及變化的鉆孔數(shù)據(jù)，生成用于確定礦體三維品位模型更新范圍的搜索橢球，并通過搜索橢球進行礦體三維品位模型更新范圍的判斷；

步驟5、礦體三維品位模型更新：根據(jù)步驟4的礦體三維品位模型更新范圍的判斷，進行更新范圍內(nèi)的礦體三維品位模型更新，并檢驗結(jié)果。

作為本發(fā)明的進一步優(yōu)化，所述步驟2中，依據(jù)鉆孔數(shù)據(jù)及勘探線信息，建立地統(tǒng)計學(xué)所要用的概率密度函數(shù)是通過多點地質(zhì)統(tǒng)計學(xué)方法獲取的，

在多點地質(zhì)統(tǒng)計學(xué)中，包含n個數(shù)據(jù)點的數(shù)據(jù)事件出現(xiàn)的概率為：

其中，U_a為n個數(shù)據(jù)點，為多點中心指示協(xié)方差。

作為本發(fā)明的進一步優(yōu)化，所述步驟1-步驟5均使用移動四面體作為礦體三維品位模型的基本單元，并將四面體重心作為該單元的坐標，即該重心的屬性值即為該四面體單元的屬性值。

作為本發(fā)明的進一步優(yōu)化，所述步驟4中，生成的用于確定礦體三維品位模型更新范圍的搜索橢球，是通過半變異函數(shù)進行確定的，