5.9. Погрешности подсчета запасов и методы их оценки

Опыт эксплуатации месторождений показывает, что между данными разведки и эксплуатации всегда имеются расхождения в значениях параметров оруденения (контуров рудных тел, количества и качества руды и пр.). Если они незначительные, то их относят к категории погрешностей подсчета запасов. Если же расхождения велики, то принято говорить о неподтверждении запасов или качества минерального сырья. Теоретически расхождения неизбежны, так как разведка месторождений осуществляется по прерывистой сети наблюдений. За пределами разведочных выработок параметры оруденения находят путем интерполяции и экстраполяции, и тем самым они являются приближенным отражением действительности. С другой стороны, при эксплуатации имеют место потери части руды и попадание вмещающих пород в руду (разубоживание или засорение), которые оцениваются приближенно и также вносят свой вклад в расхождения.

Подсчет запасов не только должен содержать достоверные данные по количеству и качеству руды, но и возможно более точно показывать пространственное положение, морфологию, строение и условия залегания рудных тел, расположение и характер тектонических деформаций, изменчивость мощности рудных тел и качества руд в пространстве. Искажение или излишняя схематизация условий залегания, пространственного положения и строения рудных тел нередко приводят к тому, что при добыче рудные тела оказываются не там, где их предполагали, мощность и качество руд отличаются от ожидаемых и т.п. Все это может служить причиной неоправданных затрат, а иногда и убытков предприятия, нарушения ритма его работы, невыполнение планов добычи, увеличения потерь и разубоживания руды при добыче и переработке.

Различают две группы погрешностей подсчета запасов: технические, возникающие в процессе измерений исходных параметров, и методические  в процессе обобщения геологической и разведочной информации.

Технические погрешности включают погрешности измерения координат точек входа и выхода разведочных выработок из рудного тела, мощности рудных тел, элементов их залегания, качества руд и других параметров. Технические погрешности делят на случайные и систематические.

Случайные погрешности неизбежны при любых измерениях, и нужно стремиться к уменьшению их до разумных пределов. Иногда они лимитированы нормативными документами, например допустимые случайные погрешности химических анализов (см. раздел 3.11). В отличие от случайных систематические погрешности считаются недопустимыми, и нужно принимать меры к их предупреждению и устранению, хотя, к сожалению, не всегда удается обнаружить их в процессе разведки.

Для предупреждения и устранения технических погрешностей необходимо знать возможные причины их появления. Наиболее надежно параметры оруденения определяют при разведке месторождений горными выработками. Если рудное тело имеет четкие границы, то их можно измерить мерной лентой или рулеткой с точностью 1-2 см. Элементы залегания измеряют горным компасом с точностью 3-5. Если границы рудного тела определяют по данным опробования, то точность их измерения зависит от длины секционных проб и может составлять 0,3-0,5 м и более.

Бурение  наиболее распространенный способ разведки месторождений  по точности определения параметров оруденения уступает горным выработкам. При колонковом бурении контакты рудного тела, его мощность, строение и состав руд устанавливаются по керну скважин. На достоверность определения параметров влияют искривление скважин, избирательное истирание и неполный выход керна. Искривление скважин начинает сказываться с отметок 100-150 м и усиливается с глубиной. Для контроля искривления скважин производятся инклинометрические замеры через 20-50 м. Зенитные углы измеряются с точностью 15-30, азимутальные  с точностью 1‑3. При малом зенитном угле (обычно менее 3) азимутальные углы измерить невозможно. Несмотря на инклинометрию, положение ствола скважины в пространстве все равно определяют с погрешностью, которая зависит от характера неоднородности горных пород, пересекаемых скважиной, и растет с глубиной. Практика показывает, что на глубине 500 м ошибка в определении положения ствола скважины в плане составляет 1-2 м, а на глубине 1000 м возрастает до 10 м и более. Для уменьшения погрешности инклинометрии применяется направленное бурение скважин и сгущение до 5-10 м сети инклинометрических замеров в участках сильного искривления скважин.

Глубину забоя скважины измеряют по длине бурового снаряда (колонны штанг и колонковой трубы). За счет деформации штанг возникает погрешность в определении глубины скважины и, как следствие, в определении местоположения рудных тел, которое на глубине порядка 1000 м может достигать 0,5-2 м. Для оценки этой погрешности производят систематические контрольные замеры глубин скважин через 100-200 м углубки.

На достоверность разведки влияет неполный выход и избирательное истирание керна. Для большинства месторождений установлен минимальный выход керна 70 % на один рейс бурения. Для руд, сложенных весьма мягкими или хрупкими минералами (молибденитом, киноварью и пр.), минимальный выход керна принимается 80‑90 %. Если фактический выход керна меньше допустимого, то данные опробования по керну считаются недостоверными. Но и при допустимом выходе керна 70-90 % границы рудных тел и мощности рудных тел определяются с погрешностью 0,1-0,5 м. При избирательном истирании керна возможно занижение содержания в руде мягких и хрупких минералов и связанных с ним полезных компонентов, в результате чего искажаются и состав руды и границы рудных тел, определяемые по данным опробования. Для повышения выхода керна применяют двойные колонковые трубы, укороченные рейсы бурения и другие технологические приемы.

В скважинах ударно-канатного бурения и при низком выходе керна при колонковом бурении границы рудного тела и качество полезного ископаемого определяют по шламу. При этом необходимо считаться с тем, что имеет место запаздывание подъема и неполное извлечение шлама, особенно тяжелых минералов в нем (например, самородного золота), в результате чего границы рудного тела смещаются вниз, иногда уменьшается мощность рудного тела. Для предупреждения погрешностей принимают меры к более полному извлечению шлама после каждого рейса бурения. Если стенки скважин не закреплены обсадными трубами, то возможна потеря части шлама в трещинах вмещающих пород и, наоборот, попадание безрудных вмещающих пород в шлам. Для контроля и предупреждения возможных погрешностей объем вынутого шлама сопоставляют с теоретическим объемом скважины по каждому рейсу бурения.

Данные бурения во многих случаях могут быть проконтролированы и исправлены геофизическими методами  каротажем скважин. Каротаж позволяет уточнить контакты, мощность и строение рудного тела, а иногда определить и качество руд (см. раздел 3.9). Например, на месторождениях магнетитовых руд данные опробования по керну контролируются каротажем магнитной восприимчивости. Если данные по керну и по каротажу различаются более чем на 0,5-1 м, осуществляется корректировка интервалов химического опробования керна по каротажным диаграммам, а при низком выходе керна содержание железа принимается по данным каротажа.

К техническим относятся погрешности опробования, подробно рассмотренные в главе 3, ошибки геологической документации, которые могут возникнуть из-за недостаточной квалификации геологического персонала: неправильное определение типов горных пород и руд, неточное установление контактов между ними, пропуски важных текстурных и структурных признаков, характерных минералов, органических остатков, признаков зон разрывных нарушений, неправильное определение элементов залегания горных пород и тектонических контактов. Во избежание ошибок необходимо поручать геологическую документацию квалифицированному геологическому персоналу, проводить систематическое его обучение, иметь эталонные коллекции горных пород и руд и периодически контролировать качество геологической документации.

Методические погрешности можно разделить на геологические погрешности и погрешности аналогии.

Геологические погрешности играют наиболее существенную роль. Они вызваны ошибками в расшифровке геологического строения месторождения, в определении количества и условий залегания рудных тел, в неправильной увязке их между разведочными выработками. Причинами ошибок могут быть недостаточно густая сеть разведочных выработок, низкое качество геологической документации, субъективизм в понимании генезиса и факторов контроля оруденения, недостаточный учет пострудных факторов преобразования месторождений (тектонических деформаций, размыва, выщелачивания) и т.п.

В рудных зонах, состоящих из множества рудных тел, формальное оконтуривание зон по крайним рудным пересечениям может привести к занижению объемов и запасов руды (рис.38). В этих условиях контуры подсчетных блоков следует проводить не по данным опробования, а по геологическим границам рудных зон.

контуры подсчетных блоков следует проводить не по данным опробования, а по геологическим границам рудных зон.

Наиболее сложные проблемы в оконтуривании рудных тел возникают на месторождениях с весьма неравномерным прерывистым оруденением. На таких месторождениях по мере сгущения разведочной сети в ранее выделенном контуре рудного тела появляются безрудные или некондиционные участки  "окна", которые пропускают при редкой сети. Чем гуще разведочная сеть, тем больше выявляется таких "окон", в результате уменьшаются запасы руды, а иногда и снижается ее качество. Эта погреш

ность носит систематический характер и, по-видимому, неизбежна на месторождениях с прерывистым оруденением, что показано в ряде работ САИГИМСа, например на сурьмяно-ртутных месторождениях Средней Азии.

Возможна и обратная ситуация, когда внутри рудного тела, сложенного рядовыми рудами, имеются небольшие участки высокосортных руд. При редкой разведочной сети они выявляются не полностью, при сгущении сети доля высокосортных руд растет, что установлено, например, на Соколовском железорудном месторождении.

Погрешности аналогии обусловлены способом распространения данных по разведочным выработкам на пространство между ними, включая построение контуров рудных тел и подсчетных блоков и обоснование наиболее вероятных значений параметров, лежащих в основе подсчета запасов. Погрешность аналогии зависит от изменчивости оруденения, от достигнутой плотности разведочной сети и от выбранного метода подсчета запасов. Чем сильнее изменчивость оруденения и чем реже расположены разведочные выработки, тем, естественно, больше будет погрешность аналогии.

В каждом методе подсчета запасов существуют определенные правила оконтуривания рудных тел, определения их объема и распространения данных по разведочным выработкам на эти объемы. На проекциях и разрезах сложные по форме рудные тела изображают упрощенными геометрическими фигурами, чаще всего многоугольниками из прямолинейных отрезков. Рудные тела заменяют простыми геометрическими фигурами: призмами, клиньями и пирамидами. Перечисленные аппроксимации реальных контуров и объемов рудных тел являются неизбежными допущениями, обусловленными дискретностью разведочной сети, и служат одной из причин появления погрешностей аналогии.

Вторая группа погрешностей аналогии обусловлена способом распространения данных внутри геометрических тел. Так, в методах геологических блоков и параллельных сечений средний состав руд, подсчитанный по разведочным выработкам, распространяется на объем подсчетного блока, т.е. состав руды в блоке принимается условно- постоянным. В методе ближайшего района данные по разведочной выработке распространяются на пространство вблизи выработки  на ближайший район, в пределах которого мощность рудного тела и состав руд также предполагаются постоянными.

Более точным отражением действительности являются интерполяционные методы распространения данных по разведочным выработкам, которые реализуются в различных сеточных методах подсчета запасов. В этих методах предполагается, что значения параметров оруденения меняются в пространстве по тому или иному закону, что является некоторым приближением к действительности, хотя и не совпадает с ней. Практика показывает, что наиболее объективное представление о пространственном поведении параметров оруденения в рудных телах получается в сеточных методах подсчета запасов с использованием геостатистической модели и модели обратных расстояний.

Оценка погрешности подсчета запасов в процессе разведки является довольно сложной операцией. Существуют лишь некоторые частные приемы ее решения.

В существующих нормативных документах рекомендуется осуществлять подсчет запасов несколькими методами. Удовлетворительная сходимость результатов нескольких методов служит некоторой гарантией от крупных ошибок, хотя и не дает представления о вероятной погрешности подсчета запасов. Естественно, что при подсчете запасов нельзя устранить систематические погрешности, имеющиеся в исходных данных, они должны быть выявлены и исправлены ранее.

Один из приемов оценки ожидаемой случайной погрешности запасов разработан для метода геологических блоков. По этому методу запасы руды в блоке (см. раздел 5.6)

,

а запасы компонента в руде

.

Предполагая, что площадь блока S, средняя мощность m_ср, средняя плотность _ср и среднее содержание компонента C_ср являются независимыми случайными величинами, ожидаемую относительную среднеквадратичную случайную погрешность запасов руды можно выразить формулой

,

а запасы компонентов в руде

где _С  относительная случайная погрешность определения среднего содержания руды. Эти формулы получены на основе теоремы сложения дисперсий независимых случайных величин.

Абсолютная погрешность площади _S, согласно исследованиям С.Н.Куличихина, равна одной четвертой части так называемой приконтурной полосы: _S = S_прик/4. Отсюда относительная погрешность площади _S = S_прик /4S. Приконтурная полоса  полоса неопределенности положения границы рудного тела. Она занимает площадь между минимально и максимально возможными границами рудного тела. Приемы определения приконтурной полосы зависят от способа оконтуривания рудного тела. Например, на рис.39 контур рудного тела построен путем интерполяции на половину расстояния между рудными и безрудными скважинами. Минимальный контур проходит по заведомо рудным скважинам (по опорным точкам), а максимальный  по заведомо безрудным. Если безрудных скважин недостаточно для оконтуривания, их можно мысленно добавить с учетом принятой плотности разведочной сети.

(17)

Если эти параметры являются зависимыми (автокоррелированными) случайными величинами, то применяется модель стационарной случайной функции, и относительная погрешность

где r(h)  коэффициенты автокорреляции мощности или содержания. Возможно применение и других математических моделей изменчивости.

Погрешность средней плотности может быть найдена аналогичными приемами. Обычно измерения плотности рассматриваются как независимые случайные величины, тогда

где n  количество определений плотности.

Если плотность определяется по зависимости от состава руды, то погрешность

где r  коэффициент корреляции (для линейных зависимостей) или корреляционное отношение (для нелинейных зависимостей) между плотностью и составом руды.

При расчете относительных погрешностей можно учесть технические погрешности измерений. Например, к относительной погрешности среднего содержания _C, определенной по формуле (17), можно добавить случайную погрешность опробования или химического анализа:

Возможны случаи, когда погрешность определения запасов трудно или невозможно оценить, например, когда запасы категории С₂ оконтурены путем неограниченной экстраполяции или подвески на глубину.

Пример. На рис.39 приведен план опробования рудного тела по скважинам по сети 100  100 м. Контур рудного тела получен путем интерполяции на половину расстояния между рудными и безрудными скважинами. Средняя плотность руды  = 3,5 т/м³, погрешность ее определения 3 %. Требуется подсчитать запасы руды и металла и определить погрешность подсчета запасов.

Площадь рудного тела, подсчитанная путем разбиения на простые геометрические фигуры, S = 305 тыс.м², средняя мощность m_ср= 2,5 м, среднее содержание металла C_ср = 1,08 %. Отсюда получаем объем рудного тела V = 762 тыс.т³, запасы руды P = 2,67 млн т, запасы металла q = 28,8 тыс.т.

Площадь приконтурной полосы равна 210 тыс.м², следовательно, абсолютная погрешность определения площади 52,5 тыс.м², а относительная погрешность 17,2 %.

Коэффициент автокорреляции мощности в широтном направлении r(h) = -0,186 при погрешности _r = 0,20, в меридиональном направлении r(h) = -0,154 при погрешности _r = 0,21, следовательно, мощность является независимой случайной величиной и к ней нужно применить статистическую модель. Имеем коэффициент вариации V_m = 29,5 %, число разведочных выработок n = 31, относительную погрешность средней мощности _m = 5,3 %.

Содержание металла в руде является зависимой случайной величиной, в чем можно убедиться, если рассчитать автокорреляцию содержаний по широтному направлению r(h) = 0,393 (_r = 0,18) и по меридиональному направлению r(h) = 0,400 (_r = 0,18). Близость коэффициентов автокорреляции по обоим направлениям свидетельствует об изотропности изменчивости. Для расчетов примем среднее значение r(h) = 0,396. На основе стационарной случайной функции найдем коэффициент вариации содержания V_С = 22,4 % и относительную погрешность _С = 3,7 %.

Теперь имеются все данные для расчета относительной погрешности запасов руды и металла:

Главное значение в данном примере имеет погрешность измерения площади рудного тела. Ее можно уменьшить, сгустив сеть разведочных выработок по периферии рудного тела.