国产日韩欧美一区二区三区三州_亚洲少妇熟女av_久久久久亚洲av国产精品_波多野结衣网站一区二区_亚洲欧美色片在线91_国产亚洲精品精品国产优播av_日本一区二区三区波多野结衣 _久久国产av不卡

?

SNESIM和FILTERSIM算法在礦體建模中的比較研究

2021-06-06 09:51:36高路萍劉曉明王濤王雨鄧?yán)?/span>
采礦技術(shù) 2021年3期
關(guān)鍵詞:樣板礦體建模

高路萍,劉曉明,王濤,王雨,鄧?yán)?/p>

(1.武鋼資源集團(tuán)大冶鐵礦有限公司, 湖北 黃石市 435006;2.長(zhǎng)沙迪邁數(shù)碼科技股份有限公司, 湖南 長(zhǎng)沙 410025)

0 引言

現(xiàn)有的礦體建模軟件多是基于人工交互的方式,并遵循以下流程:手動(dòng)繪制地質(zhì)剖面圖;通過人機(jī)交互的方式生成線框模型[1]。對(duì)于復(fù)雜的地質(zhì)模型,步驟二需要花費(fèi)地質(zhì)人員大量的時(shí)間和精力,同時(shí)易出現(xiàn)退化三角形、開口、自相交等問題[2?3]。而自動(dòng)建模方法可以通過隱式函數(shù)插值自動(dòng)重建三維地質(zhì)模型。因此,各類隱式建模方法備受地質(zhì)人員的廣泛關(guān)注,其中油氣藏領(lǐng)域的多點(diǎn)地質(zhì)統(tǒng)計(jì)學(xué)隨機(jī)建模方法是隱式建模研究的熱點(diǎn)方向之一[4]。

GUARDIANO and SRIVASTAVA等人最先提出使用訓(xùn)練圖像的高階統(tǒng)計(jì)量進(jìn)行三維空間模擬[5]。此后經(jīng)過30年的發(fā)展,學(xué)者們提出了多種不同的多點(diǎn)地質(zhì)統(tǒng)計(jì)學(xué)(MPS)模擬算法,總體上可劃分為數(shù)理統(tǒng)計(jì)算法和計(jì)算機(jī)圖形學(xué)算法等[6?8]。STREBELLE[9]提出了SNESIM算法,使用搜索樹來存儲(chǔ)訓(xùn)練圖像的掃描結(jié)果;在此基礎(chǔ)上,MARIETHOZ[10]提出了一種直接采樣DS算法,大大提高了模擬效率。ARPAT等[11]基于計(jì)算機(jī)圖形學(xué)中模式識(shí)別方法提出了SIMPAT算法,但只能應(yīng)用于分類變量模擬。ZHANG等[12?13]提出了FILTERSIM算法,該方法使用多種二維過濾器來分辨空間模式。CHATTERJEE[14]提出了WAVESIM算法,使用小波濾波對(duì)訓(xùn)練模式進(jìn)行分類。MPS模擬核心為3部分:訓(xùn)練圖像、數(shù)據(jù)事件、多點(diǎn)概率,其中訓(xùn)練圖像應(yīng)包含模擬區(qū)域完整的地質(zhì)特征[15]。因此,MPS模擬通常需要構(gòu)建龐大的訓(xùn)練圖像庫,以便數(shù)據(jù)事件掃描提取多點(diǎn)概率[16]。對(duì)于金屬礦床建模,圈定地質(zhì)剖面圖、平面圖是既定的工作流程,這些包含地質(zhì)學(xué)家認(rèn)識(shí)和認(rèn)知的二維圖像是構(gòu)建訓(xùn)練圖像庫最好的數(shù)據(jù)來源[17]。

目前MPS方法大多應(yīng)用于油氣藏領(lǐng)域建模,鮮有學(xué)者將其應(yīng)用于礦體建模。為了研究MPS算法在金屬礦床建模中的適用性,本文使用地質(zhì)統(tǒng)計(jì)學(xué)軟件(SGeMS)中的SNESIM和FILTERSIM算法來模擬一個(gè)合成的礦體模型,最后將這兩種算法的模擬結(jié)果與地質(zhì)學(xué)家人機(jī)交互的建模結(jié)果進(jìn)行比較,以此來判定兩種算法的適用性。

1 多點(diǎn)地質(zhì)統(tǒng)計(jì)學(xué)原理及關(guān)鍵算法

多點(diǎn)地質(zhì)統(tǒng)計(jì)學(xué)方法[18]是模擬空間多個(gè)點(diǎn)(>2)聯(lián)合分布的地質(zhì)統(tǒng)計(jì)學(xué)方法,具體思路為:根據(jù)地勘數(shù)據(jù)構(gòu)建多點(diǎn)搜索樣板掃描訓(xùn)練圖像,獲得相應(yīng)的多點(diǎn)概率密度函數(shù)或相似訓(xùn)練模式,最后用模擬結(jié)果填充建??臻g。

如圖1所示,當(dāng)數(shù)據(jù)樣板中條件數(shù)據(jù)dn與訓(xùn)練圖像中模式相同時(shí),將其記為有效模擬。因此,多點(diǎn)概率計(jì)算公式為:

式中,Nn為訓(xùn)練圖像中樣板數(shù)目;c(dn)為數(shù)據(jù)樣板在訓(xùn)練圖像中有效重復(fù)次數(shù)。

圖1 數(shù)據(jù)樣板、訓(xùn)練圖像對(duì)應(yīng)關(guān)系

當(dāng)模擬空間中存在K種不同物質(zhì)時(shí),空間任意位置的條件概率分布函數(shù)為P{S(u)=sk|dn},其中dn為鉆孔采樣數(shù)據(jù)。

未知節(jié)點(diǎn)為屬性1的數(shù)目有1個(gè),即c1(dn)=1;未知點(diǎn)為屬性2的數(shù)目有2個(gè),即c2(dn)=2。因此,未知點(diǎn)為屬性1的概率可定為1/3,屬性2的概率為2/3。

1.1 SNESIM算法

STREBELLE在ENESIM算法的基礎(chǔ)上提出了一種非迭代的SNESIM算法[19]。為了提高解算效率,STREBELLE提出用數(shù)據(jù)樣板一次性整體掃描訓(xùn)練圖像,然后將結(jié)果保存在搜索樹中,在之后的未知點(diǎn)模擬時(shí)只需要從搜索樹中提取相應(yīng)統(tǒng)計(jì)結(jié)果。REMY等[20]在SGeMS軟件中開發(fā)了SNESIM仿真算法,詳細(xì)算法流程如圖2所示。

VRIES等[21]認(rèn)為搜索樹結(jié)構(gòu)能夠大幅減少M(fèi)PS模擬的時(shí)間,因MPS單次模擬過程中,數(shù)據(jù)樣板尺寸是恒定的,因此僅需特定數(shù)據(jù)樣板在訓(xùn)練圖像中掃描一次即可提取所有有效的訓(xùn)練模式。圖3為訓(xùn)練圖像與搜索樣板示意圖。搜索樣板掃描訓(xùn)練圖像后生成的搜索樹結(jié)構(gòu)如圖4所示。搜索樹的第一個(gè)節(jié)點(diǎn)對(duì)應(yīng)位置u,相應(yīng)的搜索樹的級(jí)數(shù)對(duì)應(yīng)ui(0≤i≤4)。搜索樹的每個(gè)節(jié)點(diǎn)在下一級(jí)中均會(huì)被分為K個(gè)新節(jié)點(diǎn),其中K是模擬空間中可能的類別總數(shù)。搜索樹每個(gè)節(jié)點(diǎn)的值代表當(dāng)前樣板結(jié)構(gòu)在訓(xùn)練圖像中搜索到的有效模式數(shù)目。顯然,搜索樣板模型越復(fù)雜,搜索樹的層級(jí)就越多,對(duì)應(yīng)的模擬結(jié)果也就越準(zhǔn)確。

圖2 SNESIM算法流程

圖3 訓(xùn)練圖像與搜索樣板

圖4 搜索樹結(jié)構(gòu)

1.2 FILTERSIM算法

FILTERSIM算法使用一些線性過濾器對(duì)訓(xùn)練模式進(jìn)行分類,從而實(shí)現(xiàn)降維的目的[22]。模擬期間首先確定距離條件數(shù)據(jù)事件最接近的原型分類,然后選擇與條件數(shù)據(jù)最相似的訓(xùn)練模式,將其粘貼回模擬網(wǎng)格。SNESIM算法將所有訓(xùn)練的重復(fù)結(jié)果保存在搜索樹中,而FILTERSIM算法只保存了每個(gè)訓(xùn)練模式的中心位置,這樣顯著減少了對(duì)內(nèi)存RAM的需求。FILTERSIM算法主要分為3步:濾波分( filter score)計(jì)算、模式分類和模式模擬,詳細(xì)算法流程如圖5所示。

圖5 FILTERSIM算法流程

FILTERSIM算法能夠接受兩種形式的過濾器:默認(rèn)過濾器和用戶自定義過濾器。默認(rèn)情況下,F(xiàn)ILTERSIM算法為X、Y、Z方向提供了3個(gè)過濾器f(均值、梯度和曲率),見式(4)~式(6),其中過濾器結(jié)構(gòu)和搜索樣板完全相同。ni是i方向的模板尺寸(i表示的是X、Y、Z);mi=(ni?1)/2和αi=?mi,…,+mi是過濾器節(jié)點(diǎn)在i方向的偏移量。過濾器可以將訓(xùn)練模式概括為一組分?jǐn)?shù)ST(u),其中u為模板中心節(jié)點(diǎn);pat(u+(αi))為模擬節(jié)點(diǎn)值;f=nx*ny*nz,n為某一方向模板尺寸。

2 建立三維礦體模型

本文選擇將SNESIM和FILTERSIM算法的模擬結(jié)果與人工圈定的礦體模型進(jìn)行比較,以此來評(píng)價(jià)這兩種算法在礦體建模中的適用性。為了便于比較MPS中兩種算法的優(yōu)越性,本文開發(fā)了合成的礦體模型,礦體模型中44個(gè)鉆孔數(shù)據(jù)為MPS模擬硬數(shù)據(jù)。此外,礦山地質(zhì)人員在N?S方向繪制了8個(gè)地質(zhì)平面圖作為訓(xùn)練圖像。圖6為鉆孔數(shù)據(jù)及人工圈定礦體模型圖,鉆孔數(shù)據(jù)可以作為MPS模擬的硬數(shù)據(jù),圖中包含的3種顏色分別代表不同的屬性(藍(lán)色為圍巖,綠色為1號(hào)礦體,紅色為2號(hào)礦體)。為了保證模擬結(jié)果的可比性,SNESIM和FILTERSIM算法參數(shù)基本一致且均是較優(yōu)參數(shù)。

圖6 鉆孔數(shù)據(jù)及人工圈定礦體模型

2.1 SNESIM算法模擬結(jié)果

SGeMS軟件中SNESIM算法計(jì)算參數(shù)為:模擬實(shí)現(xiàn)的數(shù)量為3;初始化偽隨機(jī)數(shù)為211175;分類變量類別為圍巖、1號(hào)礦體、2號(hào)礦體,比例為0.78:0.17:0.05;搜索模板包含90個(gè)節(jié)點(diǎn),最大值為150,中間點(diǎn)為100,最小值為75;伺服系統(tǒng)參數(shù)為0.5;多層網(wǎng)格參數(shù)為3;每次只輸出最終的模擬結(jié)果。

應(yīng)用圖6所示的鉆孔數(shù)據(jù)以及8個(gè)訓(xùn)練圖像進(jìn)行SNESIM算法模擬,模擬網(wǎng)格基礎(chǔ)塊尺寸為1m×1m×1m;模型規(guī)模為75m×250m×100m。圖7為SNESIM算法模擬結(jié)果,其中圖7(a)為三維模擬結(jié)果,圖7(b)~(d)分別為X在0,25,75時(shí)的剖面圖。顯然圖7的三維模型中重現(xiàn)了1號(hào)礦體和2號(hào)礦體,但整個(gè)建模區(qū)域存在較多零散的“噪聲”,并沒有得到地質(zhì)學(xué)家預(yù)期的地質(zhì)模型。

2.2 FILTERSIM算法模擬結(jié)果

SGeMS軟件中FILTERSIM算法計(jì)算參數(shù)為:模擬實(shí)現(xiàn)的數(shù)量為3;初始化偽隨機(jī)數(shù)為211175;搜索模板尺寸為41m×41m×3m;粘貼樣板尺寸為3m×3m×1m;多重網(wǎng)格參數(shù)為3;模式原型的拆分準(zhǔn)則為10×10×10;各類數(shù)據(jù)的權(quán)重指數(shù)為0.8:0.2;每次只輸出最終的模擬結(jié)果;模式分類算法為K-Means,初始聚類數(shù)量為200,之后的聚類數(shù)目為2;相似性判斷依據(jù)為模式分?jǐn)?shù)差最小。

FILTERSIM算法模擬使用了圖6所示的數(shù)據(jù)庫,目的是比較兩種方法的結(jié)果,模擬網(wǎng)格基礎(chǔ)塊尺寸以及模型規(guī)模與SNEISM算法相同。圖8為FILTERSIM算法模擬結(jié)果。顯然圖8的三維模型很好地重現(xiàn)了1號(hào)礦體和2號(hào)礦體,與SNESIM算法相比呈現(xiàn)了更好的連續(xù)性,模擬結(jié)果中“噪音”較少,得到的模擬結(jié)果符合地質(zhì)學(xué)家的預(yù)期。

圖8 FILTERSIM算法模擬結(jié)果

圖9 各數(shù)據(jù)模型中屬性占比

圖9為各數(shù)據(jù)模型中屬性占比,比較了鉆孔數(shù)據(jù)、訓(xùn)練圖像、SNESIM算法模擬結(jié)果以及FILTERSIM算法模擬結(jié)果中各類屬性的占比情況。顯然FILTERSIM算法模擬結(jié)果中各類屬性的占比與鉆孔數(shù)據(jù)、訓(xùn)練圖像中的屬性占比較為吻合,而SNESIM算法模擬結(jié)果則與之存在偏差。因此,各屬性占比統(tǒng)計(jì)結(jié)果與上述視覺對(duì)比結(jié)論一致。

3 結(jié)論

兩種典型的多點(diǎn)地質(zhì)統(tǒng)計(jì)學(xué)算法(SNESIM和FILTERSIM算法)均能夠重現(xiàn)礦體模型復(fù)雜的空間特征,兩種算法模擬過程均可在90 s內(nèi)完成。相較于人機(jī)交互的顯式建模方法,本文算法自動(dòng)化程度更高。FILTERSIM算法的模擬結(jié)果連續(xù)性更好,更具有代表性;而SNESIM算法模擬結(jié)構(gòu)隨機(jī)性更強(qiáng),導(dǎo)致三維模型存在較多“噪音”。多點(diǎn)地質(zhì)統(tǒng)計(jì)學(xué)方法仍需要進(jìn)一步的研究,使其在礦體建模中得到更廣泛的應(yīng)用。訓(xùn)練圖像庫越詳細(xì),模擬結(jié)果就越準(zhǔn)確,而采礦業(yè)往往缺乏準(zhǔn)確且詳細(xì)的地勘信息。目前該方法在礦體建模領(lǐng)域應(yīng)用的主要難點(diǎn)在于構(gòu)建合適、準(zhǔn)確的訓(xùn)練圖像庫。

猜你喜歡
樣板礦體建模
近地表礦體地下組合式連續(xù)開采技術(shù)研究
聯(lián)想等效,拓展建?!浴皫щ娦∏蛟诘刃?chǎng)中做圓周運(yùn)動(dòng)”為例
打造辣椒種植“樣板田”
Chronicle of An Epic War
Beijing Review(2020年25期)2020-06-24 06:21:32
論甲烏拉礦區(qū)斷裂構(gòu)造及控礦作用
打贏脫貧攻堅(jiān)戰(zhàn)的“人大樣板”
基于PSS/E的風(fēng)電場(chǎng)建模與動(dòng)態(tài)分析
電子制作(2018年17期)2018-09-28 01:56:44
不對(duì)稱半橋變換器的建模與仿真
3D Mine不規(guī)則礦體三維地質(zhì)建模和儲(chǔ)量計(jì)算
樣板:不成熟的臺(tái)州
能源(2018年5期)2018-06-15 08:56:00
404 Not Found

404 Not Found


nginx
南投市| 文安县| 泰和县| 上犹县| 广灵县| 吉安市| 临沧市| 葫芦岛市| 米林县| 永济市| 松溪县| 韶关市| 福泉市| 思茅市| 萨嘎县| 宜君县| 清水县| 娱乐| 乌鲁木齐市| 永顺县| 竹溪县| 福清市| 且末县| 汉川市| 襄汾县| 凤凰县| 鸡东县| 元朗区| 三门县| 灵台县| 滦南县| 台州市| 彭州市| 涿鹿县| 高青县| 白河县| 报价| 尖扎县| 丹寨县| 丹棱县| 余干县|