CC BY
12
♦
ИЗВЕСТИЯ
ТОМСКОГО ОРДЕНА ТРУДОВОГО КРАСНОГО ЗНАМЕНИ ПОЛИТЕХНИЧЕСКОГО
ИНСТИТУТА им. С. М. КИРОВА
Том 151 1966
ИССЛЕДОВАНИЕ ПЛОТНОСТИ ГОРНЫХ ПОРОД ИТАТСКОГО БУРОУГОЛЬНОГО МЕСТОРОЖДЕНИЯ РАДИАЦИОННЫМ
МЕТОДОМ
В. Е. ОЛЬХОВАТЕНКО, В. К. КУЛЕШОВ, В. В. УЧАЙКИН (Представлена научным семинаром кафедры гидрогеологии и инженерной геологии)
В последние годы радиометрические методы все более широко применяются при инженерно-геологических исследованиях. Этими методами сравнительно быстро и с достаточной точностью можно определять плотность и влажность пород.
Методы и аппаратура, разработанные в настоящее время, позволяют вести исследования как в полевых, так и в лабораторных условиях [1,2,3, 5, 6, 7].
В настоящей работе освещаются результаты исследований плотности осадочных горных пород Итатского буроугольного месторождения методом ослабления гамма-излучения.
Итатское месторождение является одним из крупнейших в Канско-Ачинском бассейне и разработку его предполагается вести открытым способом с применением мощных роторных экскаваторов. При разведке месторождения кафедрой гидрогеологии и инженерной геологии Томского политехнического института проводятся детальные инженерно-геологические исследования. Большой объем в этих исследованиях занимает определение объемного веса пород, который является расчетной характеристикой при расчетах устойчивости бортов проектируемых карьеров. Существующие стандартные методы определения объемного веса пород (метод гидростатического взвешивания, режущих колец, парафиниро-вания) довольно громоздки и требуют значительных затрат времени, поэтому в последнее время в литературе появились работы, посвященные определению плотности пород и почв путем ¡просвечивания их пучком гамма-квантов [2, 3, 5].
При просвечивании коллимированным моноэнергетическим гамма-излучением исследуемых пород интенсивность последнего уменьшается. Ослабление интенсивности гамма-излучения в результате взаимодействия гамма-квантов с породами обусловлено частичным выходом гам-ма-квантов из пучка вследствие актов комптоновского рассеяния, фотоэлектрического поглощения и образования пар.
Закон ослабления интенсивности коллимированного гамма-излучения при прохождении через породу описывается уравнением
/ = /0е-'1 дин*, (1)
где / и /0 — интенсивность излучения, регистрируемого детектором при наличии и отсутствии просвечиваемого образца породы соответственно; (-''лин — линейный коэффициент ослабления интенсивности пучка гамма-излучения данной энергии для данной породы см~"1; х — толщина просвечиваемого слоя вещества, см.
20 Известия ТПИ, т. 151.
305
Величина ^лин характеризует относительное уменьшение интенсивности коллимированного излучения слоем просвечиваемого вещества толщиной в 1 см.
Последняя пропорциональна плотности вещества и зависит от энергии излучения, атомного номера и атомного веса вещества. В нашем случае в качестве источника излучения использовался изотоп Сэ137 (662 Мэв) 0,5 мг-экв Ра. Толщина исследуемых образцов была постоянной и равнялась 20 мм.
Линейный .коэффициент ослабления связан с массовым (объемным) коэффициентом ослабления соотношением
Улин
Г-обм = -, где (2)
Уобм — объемный коэффициент ослабления, см2/г; р —плотность исследуемого образца, г/см3. Последняя определяется, исходя из приведенных выше зависимостей по формуле
1п /п — 1п/ , .> р =---, г/см-'. (3)
Усбм ' X
Условные обозначения те же, что и в приведенных выше формулах. Для исследования плотности пород в нашем случае была использована установка, блок-схема которой приведена на рис. 1. Поток гамма-
Рис. 1. Блок-схема установки для определения плотности пород в образцах радиационным методом.
1 — свинцовый контейнер; 2 — источник излучения Сэ137; 3 — исследуемый образец; 4 — коллиматор; 5 — детектор излучения; б — фотоэлектродный умножитель ФЭУ-52; 7 — свинцовая защита; 8 — катодный повторитель, собранный по схеме Уайта; 9 — интегральный дискриминатор БД-5; 10 — пересчетное устройство ПП-8; 11—стабилизированный выпрямитель ВС-22; 12 — блок питания БН-40; 13 — блок питания БН-41.
квантов от источника 2, проходя через образец 3 и коллиматор 4, попадает на детектор излучения, представляющий из себя сцинтилляционный счетчик, состоящий из сцинтилляционного кристалла 4 NaJ (Те) 70X^0 и фотоэлектронного умножителя 6 ФЭУ-52 жалюзного типа. Импульсы блока ФЭУ, возникающие в результате сцинтилляций от гамма-квантов в кристалле и переобразованные в импульсы напряжения, передаются через катодный повторитель 8, собранный по схеме Уайта [9] на лампе 6 НЗП, на вход интегрального дискриминатора 9 типа БД-5 (10).
Блок интегрального дискриминатора включает в себя инвертор, пре-дусилитель и дискриминирующий каскад. Усилитель, построенный на триод-пентоде по схеме «двойки», обеспечивает усиление ~20. Дискриминатор, выполненный на каскаде одновибратора с обеими проводящими лампами и диодом, срывающим автоколебания в отсутствие внешних сигналов, обладает чувствительностью ~ 100 мв и выдает выходной
стандартный импульс. Чувствительность дискриминатора по входу меняется аттеньюатором, введенным перед инвертором. Высокая стабильность обеспечивается отрицательной обратной связью по сигналу в предусилителе, стабилизацией по постоянному току режима ламп од-новибратора и питанием напряжением канала ламп диода в дискриминаторе. Амплитуда выходного стандартного импульса с дискриминатора достаточна для того, чтобы запустить пересчетное устройство ПП-8.
Питание ФЭУ производилось от высоковольтного стабилизированного выпрямителя 11 типа ВС-22. Питание катодного повторителя и интегрального дискриминатора осуществлялось от блоков питания 12, 13 БН-40 и БН-41. Сцинтилляционный счетчик вместе с катодным повторителем выполнен в виде единого выносного блока, помещенного в свинцовую защиту. Активность источника определялась, исходя из безопасности работы с ним при смене образцов. Размеры коллиматора определяются просвечиваемой поверхностью, возникновением краевых эффектов в кристалле и, возможно, большей скоростью счета с целью уменьшения погрешности измерений. Последняя складывается из погрешностей измерительной аппаратуры и погрешностей, обусловленных дискретностью характера радиоактивного излучения. В нашем случае за счет использования высокостабильных источников питания и выбора оптимальной работы катодного повторителя погрешность измерения аппаратуры не превышала 1,5%.
С точки зрения снижения статистической погрешности измерения важным является выбор метода ^измерения.
В случае интегрального метода измерения погрешность можно оценить следующим образом _
где ш — флуктуации напряжения и на выходе интегрирующей цепи с параметрами Я и С относительно среднего значения и выходного напряжения;
£) (1]) — дисперсия напряжения.
Случайную велинину и можно представить в виде произведения случайных величин ^ и
и=Ъ 1^2, (5)
где ^ принимает значение 1, если фотон зарегистрировался устройством, и О—в противном случае, а £2 — подчиняется распределению Р(и) (Ш, описывающему амплитудный спектр импульсов, идущих с фотоумножителя. Тогда, в соответствии с теорией вероятностей
и = м#х)м(12)% (6)
а 7)(£7) = Я(У0(У + И(У]2ЩУ + [М{12)]*0&), (7)
где М (?) — математическое ожидание величины
Аппроксимируем фотопик в амплитудном распределении (рис. 2) гауссианом с дисперсией а2 и средним значением £/0, а комптоновское распределение—равномерным от 0 до ¿У2 таким_образом, чтобы площадь под всем распределением^ была равна еЛ^0, где в — эффективность регистрации фотонов, а А/0 —среднее число фотонов, падающих на кристалл в единицу времени. Отношение площади под фотопиком к площади под всем распределением—есть фотовклад. Обозначим его через Р. Наличие порога дискриминации ¿7, приведет к тому, что эти величины изменятся. Обозначив новые величины через г* и Я*} можно записать
и.
(1 -р)
р* = р
п ,
1 —— (1 — Р)
и,
(8) (9)
Поскольку число импульсов, зарегистрированных в единицу вре мени, подчиняется распределению Пуассона, можно записать
Тогда
£>00 = 6*^0.
В (£/) = £*Л'0 О (?2) + (г'Ыоу О (У + [М (У]2 (з*УУ0) = 1
(е*и0У
8* М
уИ(^) - [/И(с2)]2
(10) (П) (12)
(13)
Математические ожидания величины ?2 и ее квадрата найдем следующим образом
Ж(У= иР(У)йи =
1 — Я*
и2-их
ши+ ир{и)йи
и!
(1 + ¿7,)
Р* £/,
05
Ж (с?) - J и*Р{Ц)<Ш
о
Отсюда
(1_Я*)2 / и\ — и\
1 "
Ц-
е'М,.
(1 -/>')(£/,+ £/,)
3 \ ¿7, - £/,
(1 - ¿7?)
¿7, - ¿Л 3
+ Р*и о
(14)
+ (/>*)2(=г + ^б)-
(15)
+ (/> )-' (а2 + ¿УЗ)
(16)
Подставляя выражения (12), (14) и (16) в формулу (4), можно получить относительные флуктуации при произвольном пороге дискриминации. Рассмотрим случай, когда дискриминация осуществляется вблизи фотопика (рис. 2). Тогда иу = иг и
и
1
2 РС
иг
г* Я* ли г* Р* 1V,
(17)
Здесь величина а2 связана с энергетическим разрешением кристалла а соотношением
а-
8 1п 2
При сравнительно невысокой скорости счета целесообразно применять нормализацию импульсов. В этом случае и выражение (17) принимает вид
/-1-(19)
У 2RCг*P*N0
При работе в счетном режиме (с использованием пересчетного прибора) относительные флуктуации в числе фотонов, зарегистрированных за I сек, определяются формулой
Ы! = 1 _ . (20)
Из сравнения выражений (17), (19) и (20) видно, что оптимальным методом измерения оказывается счетный, так как основным требованием в нашем случае является высокая точность результатов.
В выражении (20) Л^ определится по формуле
N, = 3.7-10
где
А—активность источника в кюри, у — число квантов на один распад; ¿ — геометрический фактор ослабления; _ ослабление пучка гамма-квантов материалом исследуемого образца толщиной х; В — фактор накопления.
Геометрический фактор ослабления излучения в конусе, ограниченном краем коллиматора, подсчитывается по формуле
О-0.5
1 —
Р2
где
г—радиус окна коллиматора;
Р — расстояние от источника до кристалла детектора.
Погрешность, обусловленная дискретностью характера радиоактивного излучения, не превышала 0,3%, а общая погрешность составляла не более 1,8%, что достигалось выбором оптимальных размеров сцин-тиллятора. С точки зрения эффективности размеры его должны быть как можно больше. Однако увеличение размеров кристалла ведет к удорожанию установки. В нашем случае весьма удачным явилось сочетание ФЭУ-52 и кристалла Ыа J (Те) 70X50. Диаметр фотокатода умножителя полностью совпадает с диаметром кристалла. В то же время фотоумножитель имеет эффективную систему светосбора, а эффективность регистрации квантов сцинтилляционным кристаллом равна 0,75 (13). Высокая эффективность сцинтиллятора дает возможность производить отсечку (дискриминацию) шумовых импульсов ФЭУ и импульсов, обусловленных рассеянием гамма-квантов. На рис. 2 представлено спектральное распределение амплитуд импульсов, поступающих на дискриминатор. Как видно из рис. 2, площадь под кривой, характеризующей распределение шумовых импульсов и импульсов, обусловленных рассеянием гамма-квантов, довольно велика. Эти импульсы не несут в себе полезной информации и являются помехой, увеличивающей погрешность измерений. Установление порога дискриминации в области провала позволяет свести погрешности до минимума.
Определение плотности производилось на песчаниках, алевролитах и аргиллитах, пользующихся наибольшим распространением на территории Итатского буроугольного месторождения. Данные породы имеют близкий минералогический состав, несколько отличаются по гранулометрическому составу и в большей степени структурным и текстурным
ЕЕЗ / Эз
Рис. 2. Спектральное распределение амплитуд импульсов напряжений на входе интегрального дискриминатора.
1 — спектр импульсов при отсутствии образца породы; 2 — спектр импульсов при наличии исследуемого образца; 3—область спектра, состоящая из шумов ФЭУ и рассеянного в коллиматоре излучения; 4 — область, характеризующая комп-тоновское рассеяние в кристалле и исследуемом образце; 5 — область, характеризующая фотопоглощение.
особенностям. Вещественный состав незначительно влияет .на результаты определения плотности пород радиационным методом.
В результате проведенных исследований получена зависимость интенсивности регистрируемого излучения от плотности пород, представленная на рис. 3. Из рис. 3 видно, что эта зависимость является прямолинейной в интервале плотности от 1,70 до 2,36 г/см3 и данной точности измерений. Здесь же показана расчетная и экспериментальная прямые, которые почти полностью совпадают. В связи с этим при определении плотности пород, близких по своему составу и состоянию к исследованным образцам, можно пользоваться как расчетной, так и эксперимент тальной прямой.
Обработка полученных материалов методами математической статистики показала, что частные значения объемного веса и интенсивности излучения подчиняются нормальному закону распределения и имеют тесную корреляционную связь. Коэффициент корреляции оказался 0,83. Наряду с определением плотности отмеченных пород радиационным методом, плотность их по всем пробам определялась еще обычным стан-
й
CZH / CD ? ПП3
Рис. 3. Зависимость скорости счета от плотности пород. I — экспериментальная прямая; II — расчетная прямая; 1—песчаники; 2 — аргиллиты; 3 — алевролиты.
дартным способом. Сопоставления результатов показали, что расхождение не более +0.05 г/см3. Это находится в пределах принятой точности определения данной характеристики.
ЛИТЕРАТУРА
1. А. Е. Б а б и н е ц, С. Т. 3 в о л ь с к и й. Исследование плотности и влажности грунтов методом радиоактивных излучений. Труды института геологических наук АН УССР. Вып. 6. Киев, 1961.
2. М. П. Беликов, В. А. Емельянов, В. Е. Нестеров. Применение радиоактивных изотопов в гидротехническом строительстве. Госстройиздат, 1961.
3. Е. К- В а р в а р и н а, Е. М. Ф и л л и п о в. Определение плотности горных пород в образцах по ослаблению гамма-излучения. Бюллетень научно-технической информации № 1 (51), 1964.
4. Н. В. Г л а з о в. Применение радиоактивных изотопов в инженерных изысканиях. Госатомиздат, 1962.
5. В. А. Емельянов. Гамма-лучи и нейтроны в полевых почвенно-мелиоратнв-ных исследованиях. Госатомиздат, 1962.
6. В. А. Емельянов. Гамма-гамма метод измерения плотности почв и грунтов. Гидротехника и мелиорация № 1, 1965.
7. Изотопы измеряют влажность. Информационный бюллетень ВДНХ № 1, 1965.
8. Э. Б. М и л е в с к и й. Радиационный контроль и измерение изделий. Машгиз, 1963.
9. А. А. Санин. Электронные приборы ядерной физики, 1964.
10. Техническое описание ЖШ 2, 803.01« ТО, 1964.
11. Техническое описание ЖШ 2. 087.164 ТО, 1964. Т2. Техническое описание ЖШ 2. 087.165 ТО, .1964.
13. Anger Н. О. und Davis D. Н. Gamma-Ray Detection Efficiency und Jage Resolution in Sodium Jodide, RSI, 35, № 6, '1-964.
