1、估算儲(chǔ)量誤差的方法探討勘探階段所獲得的各級(jí)儲(chǔ)量，其誤差范圍實(shí)際上是不可知的，無(wú)論用現(xiàn)行的一些什么方法，都 不能肓接冋答礦量精度。而探采對(duì)比已為時(shí)過(guò)晚，不能指導(dǎo)同一礦床的勘探。例如，蘇聯(lián)克里活洛 格恢礦多年探采實(shí)踐證明，規(guī)范定的網(wǎng)度a密，放稀2 -4倍尚能滿(mǎn)足礦量粕度要求。雖然通過(guò)反 復(fù)研究已得出這個(gè)結(jié)論，但在勘探中只能對(duì)法定網(wǎng)度稍許放寬而不能徹底改變，造成很人的浪費(fèi)。本文試圖論證在勘探過(guò)程中，直至交報(bào)告時(shí)，用一定的置信概率估算礦量誤差范圍的途徑。 以便有效地少用工程量，兼顧勘探工程的必要性和充分性，適時(shí)終止勘探。一、統(tǒng)計(jì)單元和定量工程首先要確定估算課差的對(duì)象（統(tǒng)計(jì)單元）是總量還是分量（分層、塊

2、段、類(lèi)型、品級(jí)等），是 單體還是復(fù)合體（多層復(fù)合，分支復(fù)合或復(fù)合透鏡體等）。總量誤差尤其是工業(yè)礦量的總量誤差， 是影響礦產(chǎn)評(píng)價(jià)的關(guān)鍵因素，是一定要計(jì)算的。首采區(qū)、儲(chǔ)量集屮區(qū)也應(yīng)當(dāng)計(jì)算，但不宜劃分過(guò)細(xì)。 分的過(guò)細(xì)，樣點(diǎn)減少，誤差范圍相應(yīng)增人，使問(wèn)題復(fù)朵化。因此，要根據(jù)地質(zhì)規(guī)律和采選工藝基本 要求，適宜地劃分統(tǒng)計(jì)單元。同時(shí)，更要明確統(tǒng)計(jì)變量確實(shí)來(lái)自同一統(tǒng)計(jì)對(duì)象。勘探工程一般具有定性、定量、定位和定價(jià)的多重目的性。或以一項(xiàng)為重，或兼而有之。般多以定量為重。當(dāng)達(dá)到一定的礦量精度要求，根據(jù)需要還可以布置定位工程（解決構(gòu)造錯(cuò)斷， 邊界座標(biāo)等），或?yàn)槿〈髽雍万?yàn)證而布置工程。定位、定性工程同時(shí)也提供了定量數(shù)據(jù)

3、。圈邊工程 對(duì)邊緣塊段儲(chǔ)量相對(duì)課差影響人，但對(duì)總礦呈影響有限，也是以定位為主的工程。落空的圈邊工程 當(dāng)然不能參與誤差計(jì)算。按丄業(yè)指標(biāo)確定的邊界工程視其數(shù)量和所占的儲(chǔ)量份額決定其是否參與計(jì) 算。在一個(gè)統(tǒng)計(jì)單元內(nèi)部的落空工程應(yīng)參與計(jì)算，它反映了礦體的連續(xù)性，包含在厚度變化的計(jì)算 之屮。但不可取零，以最小數(shù)代之。謀差性質(zhì)和影響謀差的主要參數(shù)謀差分兩人類(lèi)。一類(lèi)是測(cè)量性謀差，如工程測(cè)量、厚度換算、化驗(yàn)質(zhì)蜃等。這類(lèi)謀差主要通過(guò) 嚴(yán)格技術(shù)管理來(lái)解決。另一類(lèi)是代表性謀差，或稱(chēng)外推謀差、抽樣謀差。也可以說(shuō)以各變蜃平均值（期望值）代替總體的真值時(shí)，可能有多大的誤差。本文討論僅限于后一類(lèi)誤差。從理論上講，總體誤差是

4、不可知的，抽樣誤差只代表樣木的誤差。勘探過(guò)程則是逐步抽樣 的過(guò)程。在各個(gè)階段，每投入一項(xiàng)或一批工程，都有可以用一定的置信概率去估算礦量的課差范用。 隨著工程數(shù)目（樣本）增大，誤差必然逐漸減小。直至達(dá)到相應(yīng)勘探階段（或級(jí)別）的允許誤差范圍內(nèi)，即可滿(mǎn)足該階段勘探的定量粹度要求。計(jì)算礦量公式的諸參數(shù)屮，影響誤差的最主要的因索是厚度和品位。為簡(jiǎn)便起見(jiàn)，本文試圖以平均厚度誤差代表鐵礦石量誤差，將規(guī)范允許誤差降低20%,即c1級(jí)±36%, b級(jí)±20%。二、誤差計(jì)算和最少工程量k=it變量平均值相對(duì)誤差公式是實(shí)，9 " （1）式中px為相對(duì)誤差vx為變化系數(shù)，n為樣本（樣數(shù)）

5、，t為一定置信概率的概率系數(shù)（查表）。 公式表明，課差的人小與變化系數(shù)（表征變量的離散性）成正比，與樣數(shù)的平方根成反比，這與類(lèi) 型、網(wǎng)度、abc的基本原則是一致的。最關(guān)鍵的問(wèn)題是厚度不是隨機(jī)變量，而是區(qū)域性變量，套用隨機(jī)變量公式去計(jì)算區(qū)域性變量的變化系數(shù)與實(shí)際偏差太人（另有專(zhuān)文研究）。據(jù)b、m卡扎克的研究，按普通方法計(jì)算的厚度變化系數(shù)vm比探采資料確定的加高04倍，厚度誤差高0. 2-4倍。這樣高倍數(shù)的誤差，毫無(wú)近似可言,也沒(méi)有實(shí)用價(jià)值。因此，対vm必須進(jìn)行區(qū)域化修正,式中為vq區(qū)域化變化系數(shù)，vj為隨機(jī)變化系數(shù)，r為變量與空間的相關(guān)系數(shù)（由最大值 向各個(gè)方向變量與距離的相關(guān)系數(shù)，再求其中平均

6、值用于式中）o按筆者的研究，經(jīng)區(qū)域化修正的 厚度變化系數(shù)，約減少一半左右（天湖鐵礦），比較合乎實(shí)際，有使用價(jià)值。根據(jù)林徳柏格萊維定理,經(jīng)驗(yàn)，n>30就算充分大（大數(shù)定律）。再加大n對(duì)變量平均值精度提高的幅度逐漸減小，不是必要的。對(duì)于區(qū)域性變量厚度來(lái)說(shuō)，可以用更小的樣本即能達(dá)到同樣的粹度，例如對(duì)任何的兒何形體, 用很小的n （n條線）即能獲得相當(dāng)精確的平均值。因此，用區(qū)域化修正的vx,求最小工程數(shù)n,可以用（1）式的反轉(zhuǎn)公式,(3)它指出了一個(gè)統(tǒng)計(jì)單元（礦體或塊段）最少有多少個(gè)丄程即可達(dá)到精度要求。總之，利用區(qū)域化厚度變化系數(shù)，從初探階段開(kāi)始，隨吋對(duì)統(tǒng)計(jì)單元作定量誤差計(jì)算，不 必等到工程很

7、密時(shí)再作網(wǎng)度驗(yàn)證對(duì)比。當(dāng)估算課弟在允許范陽(yáng)內(nèi)，則定量工程是充分的，如綜合各 方而的研究，或部分地段工程控制空檔過(guò)大，還有網(wǎng)度過(guò)稀之慮，可選擇有代表性或變化最人的重 點(diǎn)塊段，插入個(gè)別驗(yàn)證工程。驗(yàn)證結(jié)果如一再落入允許范用內(nèi)，應(yīng)終止再用工程，結(jié)束該階段的勘 探。三、天湖恢礦勘探實(shí)例新鳴天湖鐵礦為埋深200m以下的盲體礦。60年磁法掃而發(fā)現(xiàn)了大界常，按界常形態(tài)良好, 附近地表有小鐵礦，判定為磁鐵礦界常。64年底界常中心的zk4在300米深入穿過(guò)4m富礦，接著 又在500m深處,打到40米貧礦和14米富礦體,證實(shí)了 1:2000磁探的推算。為億噸以上的大型鐵 礦（磁探結(jié)論是長(zhǎng)3km,厚40米，無(wú)限延深）

8、，從而結(jié)束了定性評(píng)階段。此吋對(duì)礦床評(píng)價(jià)的關(guān)鍵問(wèn)題是確定工業(yè)儲(chǔ)量的埋藏深度和質(zhì)量。有關(guān)方而經(jīng)初步計(jì)算認(rèn)為，如果六百米以上有四t萬(wàn)噸礦量即可開(kāi)發(fā)。因此決定勘探。當(dāng)時(shí)筆者提出用30個(gè)鉆孔約二力米 進(jìn)尺一完成這個(gè)使命，作了總體設(shè)計(jì)。以400200x200100米不等距網(wǎng)度勘探工業(yè)儲(chǔ)量，深部遠(yuǎn) 景控制到900米，總長(zhǎng)控制3600米。經(jīng)兩年多勘探，連定性階段共用了 32個(gè)鉆孔21000米進(jìn)尺。其中項(xiàng)端圈邊孔7個(gè)，兩側(cè)狹縮地段4個(gè)孔，1300米長(zhǎng)的主要礦段只用了 18個(gè)定量鉆孔。68年經(jīng)鞍山設(shè)計(jì)院和酒鋼等單位共 同議定，勘探程度是合理的，可以終止勘探，提供設(shè)計(jì)。與勘探前相比，總礦量變化不大，工業(yè)礦量比預(yù)計(jì)減

9、少四分之一，而且下延了 100米，降 低了天湖鐵礦的“身價(jià)”,未能及時(shí)轉(zhuǎn)入開(kāi)發(fā)。后來(lái)在七六年和八五年，為審批報(bào)告和補(bǔ)勘，圍繞勘探程度（礦量誤差）展開(kāi)了兩次爭(zhēng)論。 筆者始終認(rèn)為，天湖鐵礦的地質(zhì)規(guī)律已經(jīng)查明，嚴(yán)格受地層控制，是較穩(wěn)定的層狀一透鏡狀礦體。 下層富礦為層狀一1層狀，上層貧礦為似層狀一較規(guī)則的透鏡狀，厚薄分布與磁界常形態(tài)一致，主 礦段內(nèi)無(wú)落空鉆孔，連續(xù)穩(wěn)定。與克里沃羅格同類(lèi)型鐵礦相比，網(wǎng)度不算稀。個(gè)別鉆孔偏斜人且方 位測(cè)不準(zhǔn)，主要影響定位和度。其厚度值測(cè)斜數(shù)據(jù)和軸心夾角兩種方法計(jì)算結(jié)果一致，不至于影響 定量謀差。又據(jù)梅山鐵礦經(jīng)驗(yàn)，磁探預(yù)算和勘探儲(chǔ)量相差17%,而天湖鐵礦僅差13%,因此總

10、儲(chǔ)量和 工業(yè)礦量是對(duì)信的。最后又以抽樣誤差方法計(jì)算了塊段、工業(yè)礦量和總儲(chǔ)量的誤差范圍，證明其誤 差在允許范圍內(nèi)。至于個(gè)別地段空擋較大，厚度變化也較人，建議用少量鉆孔驗(yàn)證，如果用了 3至5個(gè)鉆孔 驗(yàn)證，其療度值都落在允許誤差范圍內(nèi)，則可以證明達(dá)到精度耍求1佃不必大動(dòng)干戈。其實(shí)在68年就 作過(guò)個(gè)別驗(yàn)證。如0號(hào)和3號(hào)剖面6個(gè)鉆孔之間空擋較人，厚度變化也大，在屮間插入zk73驗(yàn)證（鞍 山院也提過(guò)），其結(jié)果示意如圖。圖中點(diǎn)間數(shù)字為間距（米），橫線上下分別為上層礦（m2）和下 層礦51）厚度（米）如圖按周用6個(gè)鉆孔厚度計(jì)算如下卜層(ml)上層（m2）算術(shù)平均厚6. 68m1& 72m隨機(jī)變化系數(shù)v

11、mj51.6%98. 2%空間相關(guān)系數(shù)r0. 8840. 94區(qū)域化變化系數(shù)vmq24. 0%33. 6%誤差限（概率0. 95）am± 1. 6m±6. 16m驗(yàn)證孔實(shí)際厚度6. 67m13. 50m計(jì)算結(jié)果表明，驗(yàn)證孔ml與期望值相等，m2落在允許范圍內(nèi)。證明定量工程是充分的。 該地段平均厚度的相對(duì)誤差：ml為23.9%, m2為32. 9%。如把zk73加入成7個(gè)孔計(jì)算,誤差還會(huì)降 低，當(dāng)在c1級(jí)要求范圍內(nèi)。上表應(yīng)用區(qū)域化厚度變化系數(shù)計(jì)算不同置信度的誤差，得到較滿(mǎn)意的結(jié)果。工業(yè)塊段的總 誤差為30. 4%（a =0.95）或25. 4%（a二0. 90）。上層誤差較大

12、,f層誤差多小于20% （總誤差1& 1%和 15.3%）,都在允許誤差范圍內(nèi)。反過(guò)來(lái)以給定允許誤差限，求得最少工程數(shù)與實(shí)際孔數(shù)基本相符，這里就不專(zhuān)門(mén)列表論證了。四、討論和推論1. 85-88年間，犬湖鐵礦準(zhǔn)備開(kāi)發(fā)，乂進(jìn)行了詳勘。用了三年多時(shí)間，投資八百多萬(wàn)元, 打了 3. 7萬(wàn)米鉆探。五十多個(gè)鉆孔主耍用在工業(yè)塊段上。與二十年前的成果相比，無(wú)論勘探程度、 勘探質(zhì)量、綜合研究等各方面無(wú)疑是大人地捉高了。但就礦量精度而言，工業(yè)儲(chǔ)量（b+c1級(jí)）只減 少了 10%,上層品位提高了 10%,全礦總儲(chǔ)量減少了 13%左右（注：三個(gè)%是當(dāng)年聽(tīng)說(shuō)的），增減幅度 是很小的。比本文計(jì)算的誤差限小得多。證

13、明當(dāng)年對(duì)天湖鐵礦勘探的指導(dǎo)思想是對(duì)路的，所獲儲(chǔ)量 是可信的，應(yīng)用抽樣誤差理論估算礦量誤差是可行的。因此，在勘探中后期，可采用本文使用的方 法，逐步插入驗(yàn)證，以一定的置信概率估算謀差，直至每一項(xiàng)或一小批工程所得結(jié)果都落在允許范 圍內(nèi),即可保證礦量精度。2. 品位變化一般隨機(jī)性較大，單工程按單樣求得的品位變化系數(shù)，只反映該工程品位沿 厚度方向的變化性質(zhì)（離散性）。它既取決于客觀的變化性，也取決于人為的分段長(zhǎng)度。例如某一 工程穿過(guò)10米礦體，分別按lm、2m、5m甚至10m取樣，其總品位（或平均值）是一致的（變化則 不相同）。它直接代表該取樣線（厚度點(diǎn)）的品位。如冇誤差也是測(cè)量性誤差。只冇當(dāng)它與其他工 程求塊段或礦體平均品位時(shí)，才產(chǎn)生代表性謀差。因此，計(jì)算金屬量的誤差應(yīng)以工程品位為基點(diǎn)（與 厚度計(jì)算為等階量綱）。應(yīng)計(jì)算其空間相關(guān)系數(shù)，進(jìn)行區(qū)域化修正。3. 木文選用的公式早已有之，用于勘探求儲(chǔ)量課差者并不多見(jiàn)，主要是沒(méi)有解決區(qū)域化 修正問(wèn)題，而門(mén)，這種數(shù)理統(tǒng)計(jì)方法沒(méi)有考慮礦體規(guī)模、形態(tài)、邊界特征、工程分布特點(diǎn)、參數(shù)平 均的方法及