CC BY
62
В.В. Данилов, А.Д. Истомин, М.Д. Носков
МНОГОУРОВНЕВАЯ ЦИФРОВАЯ МОДЕЛЬ ОСАДОЧНОЙ ТОЛЩИ РАЙОНА РАСПОЛОЖЕНИЯ СИБИРСКОГО ХИМИЧЕСКОГО КОМБИНАТА
Работа поддержана грантами РФФИ № 06-01-00073-а, № 06-07-96907-р_офи.
Представлена многоуровневая иерархическая цифровая модель района полигона глубинного захоронения жидких радиоактивных отходов Сибирского химического комбината, которая согласованно описывает различные участки геологической среды с разным уровнем детальности. Представленная модель состоит из пространственно вложенных моделей трех рангов. В работе описаны структуры, принципы построения и схемы согласования моделей различных рангов. Эта модель может использоваться для анализа литолого-стратиграфической структуры района полигона глубинного захоронения, выполнения прогнозных расчетов изменения состояния пласта-коллектора и распространения ореола загрязнения.
Ключевые слова: многоуровневая иерархическая цифровая модель; жидкие радиоактивные отходы; глубинное захоронение.
В процессе работы предприятий атомного энергопромышленного комплекса образуется значительное количество жидких радиоактивных отходов (ЖРО). Одним из способов обращения с ЖРО является глубинное захоронение - закачка их в стратифицированные осадочные геологические формации без намерения последующего извлечения и переработки [1]. Этот метод обращения с ЖРО применяется на Сибирском химическом комбинате (СХК) - одном из крупнейших предприятий атомной отрасли России. Глубинное захоронение осуществляется через сеть нагнетательных скважин на специально оборудованном полигоне. Геологическое строение района размещения полигона глубинного захоронения является основным фактором, обеспечивающим экологически безопасную изоляцию ЖРО от области активной жизнедеятельности человека. В результате смешения ЖРО с пластовыми водами и взаимодействия с вмещающими породами формируется сложная геотехнологическая система - пласт-коллектор. Для обеспечения безопасности эксплуатации полигона глубинного захоронения необходим постоянный контроль состояния пласта-коллектора и прогнозирование процессов, протекающих в геологической среде [1]. Для решения этих задач целесообразно использовать современные геоинформационные технологии и методы математического моделирования. В последние годы был разработан ряд цифровых моделей, которые имели различные пространственные масштабы, и с разным уровнем детальности описывали геологическую среду района полигона глубинного захоронения ЖРО СХК. Например, в работах [2, 3] в ка-
честве объектов моделирования рассматривались непосредственно прифильтровые зоны нагнетательных скважин. Другие модели имели более мелкий масштаб и описывали пласты-коллекторы, а также примыкающие к ним буферные горизонты в границах полигона глубинного захоронения [1, 4, 5]. Мелкомасштабные цифровые модели рассматривали всю толщу осадочного чехла в районе расположения полигона глубинного захоронения ЖРО СХК [1, 6, 7]. Модели разных масштабов [1-7] использовали различные исходные данные, основывались на разных методах описания геологической среды и не согласовывались друг с другом. В связи с необходимостью комплексного обеспечения безопасности эксплуатации полигона глубинного захоронения ЖРО СХК, была разработана представленная в настоящей работе иерархическая многоуровневая модель, которая согласованно описывает различные участки геологической среды с разным уровнем детальности.
Район полигона глубинного захоронения ЖРО СХК расположен в краевой части Западно-Сибирской платформы, вблизи ее сочленения с Саяно-Алтайской складчатой областью [1]. Геологический разрез района представлен терригенными, преимущественно континентальными, слабосцементированными отложениями мезо-кайнозойского платформенного чехла (мелового, палеогенового, неогенового и четвертичного возраста), которые несогласно залегают на метоморфизирован-ных, сцементированных, собранных в складки породах фундамента, представленных глинистыми сланцами (рис. 1). Возраст кровли фундамента датируется верхним палеозоем.
Рис. 1. Геологический разрез района полигона глубинного захоронения ЖРО СХК
В толще осадочного чехла прослеживается неравномерное чередование песчаных и глинистых пород, что предопределяет этажное развитие водоносных и водоупорных горизонтов, характеризующихся фильтрационной неоднородностью в плане и разрезе. Наряду с общепринятой системой стратификации гидрогеологического разреза для района расположения полигона глубинного захоронения ЖРО СХК была разработана условная система, которая используется в целом ряде работ, посвященных вопросам поиска источников водоснабжения и захоронения ЖРО в этом районе [1-7]. В соответствии с данной системой, по совокупности геологических, фильтрационных и емкостных признаков, в разрезе мезо-кайнозойских отложений выделяется семь литологических горизонтов преимущественно песчаного состава и восемь относительно водоупорных глинистых слоев. Песчаные горизонты нумеруются снизу вверх римскими цифрами: I, II, III, IV, 1Уа, V, VI, а водоупорные слои обозначаются буквами латинского алфавита: А1, А2, В, С, Б, Е, Б, в. Захоронение ЖРО осуществляется в эксплуатационные горизонты II и III осадочного чехла, которые отделены от вышележащих водоносных слоев и палеозойского фундамента пачкой из глинистых водоупоров (А1, А2, В, Б, Е) и буферных горизонтов (I, IV, ^а) (рис. 1).
Для решения различных задач по обеспечению безопасности и оценки последствий эксплуатации полигона глубинного захоронения ЖРО СХК необходимо иметь согласованные модели геологической среды разных пространственных масштабов от прифильтровых областей отдельных нагнетательных скважин до всего района расположения полигона. В соответствии с этим представляемая в работе многоуровневая модель включает в себя несколько цифровых моделей различных масштабов (рангов). Ранг модели возрастает с укрупнением ее масштаба и повышением уровня детальности описания геологической среды. Анализ результатов применения цифровых моделей в работах [1-7] показал, что для контроля состояния пласта-коллектора и прогнозирования процессов, протекающих в геологической среде, необходимо использование моделей трех рангов:
- района расположения полигона (модель I ранга);
- эксплуатационных и буферных горизонтов и разделяющих их водоупорных пластов в пределах полигона (модель II ранга);
- прискважинной области (модель III ранга).
Каждая из моделей состоит из совокупности субго-
ризонтальных слоев, что обусловлено стратифициро-
Иерархическая модель включает в себя несколько моделей III ранга (прискважинных моделей), для обозначения порядкового номера которых используется индекс д. В каждой из моделей III ранга число состав-
ванной структурой моделируемой геологической толщи. Слой представляет собой тело, сверху и снизу ограниченное поверхностями, которые имеют не более одной точки пересечения с любой вертикальной линией. При этом каждая из поверхностей раздела слоев толщи одновременно является подошвой перекрывающего слоя и кровлей для слоя подстилающего, исключение составляют поверхности, соответствующие кровле и подошве всей области моделирования.
Прослои модели большего ранга выделяются из слоев модели меньшего ранга путем их расщепления. Выделение слоев осуществляется в соответствии со стратифицированным строением и технолого-гидрогеоло-гическими особенностями геологической среды. В зависимости от ранга модели, каждому слою соответствуют элементы стратифицированной толщи различного масштаба. В модели I ранга в качестве слоя рассматриваются водоносные или водоупорные горизонты, выделяемые в соответствии с системой технологической стратификации СХК. Слоям модели II ранга соответствуют преимущественно песчаные или глинистые прослои стратифицированной литологической толщи пластов-коллекторов, примыкающих к ним буферных горизонтов и водоупоров, которые их разделяют. В моделях III ранга слои описывают пропластки отдельных литологических разностей в осадочных толщах, вскрываемых фильтрами нагнетательных скважин.
Для обозначения слоев используются комбинированные индексы, которые однозначно определяют соотношения вложенности между слоями моделей различных рангов. Нумерация слоев всех моделей осуществляется сверху вниз. Слоям модели I ранга присваивается индекс 5, значение которого изменяется от 1 до 15. Соотношение между системой стратификации, принятой на СХК, и индексацией слоев в модели I ранга представлено в табл. 1. Например, модельному слою с 5 = 1 соответствует неоген-четвертичный водоносный комплекс (VI водоносный горизонт), а слою модели с 5 = 15 - первичные коры выветривания палеозойского фундамента (водоупорный пласт А1).
Прослои модели II ранга детально описывают структурные особенности эксплуатационных и буферных горизонтов, а также разделяющих их водоупоров (слои модели I ранга 5 = 7-13) в пределах территории полигона. Общее число слоев модели II ранга составляет 159 (табл. 1). Для их идентификации используется комбинированный индекс 5,5 ', в котором &•' соответствует порядковому номеру прослоя модели II ранга в слое 5 модели I ранга.
ляющих ее слоев уникально и определяется литологоструктурными особенностями геологической среды и строением фильтра конкретной нагнетательной скважины. Эти слои идентифицируются с помощью комбини-
Т а б л и ц а 1
Соотношение между различными системами индексации слоев
Система стратификации, принятая на СХК VI О V Б № Е IV Б III С II В I А2 А1
Индекс 5 слоя в модели I ранга 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
Количество слоев в модели II ранга 33 41 31 19 14 11 10 - -
рованного индекса 5,5 ,5 ,д, в котором 5,5 наследуются от индексов моделей меньших рангов (5 - от слоя модели I ранга, 5' - от модели II ранга), а 5" обозначает порядковый номер прослоя в слое 5' модели II ранга. Значение индекса д соответствует номеру модели III ранга.
Плановая дискретизация слоев модели каждого ранга осуществляется с помощью равномерной квадратной сетки. Положение узла сетки модели I ранга определяется индексом у. Сопряжение сеток моделей различных рангов достигается путем разбиения ячейки сетки модели меньшего ранга на п2 ячеек сетки моделей большего ранга (п - целое число, равное соотношению длин сторон ячеек моделей меньшего и большего рангов). Индексы, обозначающие узлы сеток моделей II и III рангов, являются комбинированными. Для модели II ранга - 1,1’,],]’, а для модели III ранга - 1,1',1",},}’,}". В обоих случаях последняя величина в индексе обозначает положение узла в сетке, а предыдущие части индекса наследуются от ячейки модели меньшего ранга.
Геометрическое строение стратифицированной толщи описывается заданием абсолютных отметок для всех поверхностей, соответствующих кровлям модельных слоев, и абсолютных отметок подошвы самого нижнего слоя, которая, одновременно, является подошвой всей моделируемой геологической толщи. Для обозначения кровли каждого из слоев модели применяется индексация, аналогичная индикации, используемой для обозначения этих слоев.
Положение каждой поверхности раздела 5 (для модели II ранга - 5,5', для модели III ранга - 5, 5,5') определяется значениями абсолютной отметки к] в узле /, ] (в модели II ранга - к,^/"5' , III ранга - к,,. ■ ш5'* ). Мощность слоя однозначно определяется значениями соседних поверхностей раздела:
- для модели I ранга
£ і з і 3
т — п — п
'пї,і 'Ч, і пї,і
- для модели II ранга
т’1. ., = КІ. .,— П^+1,
;,;,.М,і,і М,і,і
- для модели III ранга
— „з,з ,5 7 5,5 ,5 7 5,5 ,5 +1
-і -н • -і -и — П- -і -и • • і -м — П- -і -и • -і -и
М,", і , і, і м',Л і, і ', і '' г,г',;", і, і ',.''
(1)
(2)
(3)
Расчет значений сеточных функций абсолютных отметок поверхностей раздела слоев осуществляется путем выполнения площадной интерполяции данных, полученных на основе анализа геологических колонок скважин, вскрывающих исследуемую литологическую толщу. Интерполяция проводится с помощью геостати-стического метода Кригинг [8].
Геофильтрационные и емкостные свойства геологической среды описываются набором двумерных сеточных функций для каждого слоя. Значения этих функций в узлах сетки рассчитываются путем площадной интерполяции аналогично значениям абсолютных отметок. Диапазон изменения параметров геологической среды, задаваемых в модели, приводится в табл. 2.
Т а б л и ц а 2
Параметры геологической среды
Наименование параметра Диапазон изменения
Горизонтальный коэффициент фильтрации (кх, к), м/сут 0,02 - 15
Вертикальный коэффициент фильтрации (к), м/сут 1,5 • 10-5 - 1,0 • 10-2
Полная пористость (Л), % 30 - 40
Активная пористость (Л*), % 5 - 15
Упругоемкость (п), 1/м 1,0 • 10-6 - 1,0 • 10-3
Пространственная вложенность слоев и ячеек моделей обусловливает согласование всех параметров при переходе от модели меньшего ранга к модели большего ранга. Для примера рассмотрим согласование параметров для моделей I и II рангов. Из условия сохранения объема среды, при переходе от модели меньшего ранга к модели большего ранга, следуют соотношения согласования мощности:
полной пористости:
X т-;![,,
N. — ;,. ,3
™ 25
п т.,
;,3
активной пористости:
X тй
3 ДТ И«5,5
из'
N — ‘,з ,5 ;,3
п2 т*, ‘,3
упругой емкости:
(4)
(5)
(6)
Соотношения согласования латеральных и вертикального коэффициентов фильтрации вытекают из условий равенства гидропроводности ячеек в моделях большего и меньшего рангов. Коэффициенты фильтрации в модели I ранга связаны с коэффициентами фильтрации модели II ранга следующими соотношениями:
(
(Кх). — т. X X
X
ч. т;..,(Кх)
(К),/ — тГ—1X X
1
ХЛ ,,■( Кг )5
(8)
(9)
т
к) і,—-т. X X
т
55
*;,;',і',і''
п ;,з V 5і (Kz )
5,5
;,;',і',і ' у
. (10)
Модель I ранга (модель района расположения полигона) описывает всю геологическую толщу осадочного чехла Западно-Сибирской плиты в районе размещения полигона глубинного захоронения ЖРО СХК (рис. 2).
Рис. 2. Область (а) и слои (б) модели I ранга. Сплошной линией показаны границы моделируемой области модели I ранга, пунктирной - границы модели II ранга
В гидрогеологическом плане этот район является краевой, относительно обособленной частью ЗападноСибирского артезианского бассейна. Структура фильтрационного потока является главным фактором, обусловливающим выбор боковых границ модели I ранга. Направление потока фильтрации во всех водоносных горизонтах района западно-юго-западное. Областью питания всех водоносных горизонтов служит локальная водораздельная возвышенность, а областью разгрузки - долина р. Томь [9]. Поэтому в качестве восточной границы модели была выбрана линия местного водораздела, а западно-юго-западная граница проводилась по руслу р. Томь. С юга и севера модель ограничивалась, соответственно, руслами рр. Самуська и Б. Киргизка (рис. 2, а), которые являются наиболее крупными правобережными притоками р. Томь в данном районе.
Описываемая моделью I ранга территория вписана в прямоугольную область размером 25*25 км, которая покрыта сеткой с квадратными ячейками со стороной 100 м. Площадь территории, рассматриваемой в модели I ранга, составляет около 400 км2. Ячейки, которые
оказываются за границами моделируемой области, считаются неактивными и не рассматриваются в дальнейшем. Общее количество активных ячеек в каждом слое модели составляет 39569.
Модель II ранга описывает стратифицированную толщу осадочных пород на территории полигона глубинного захоронения ЖРО (рис. 3). Боковые границы модели II ранга заданы таким образом, чтобы обеспечить нахождение фильтрата отходов в пределах моделируемой области в течение всего проектного времени эксплуатации полигона. Общая площадь территории модели II ранга на порядок меньше площади модели I ранга и составляет 40 км2. В плане моделируемый массив имеет форму прямоугольника, длинная сторона которого ориентирована в субмеридиональном направлении (рис. 3, а). Сетки поверхностей раздела слоев модели II ранга формируются путем вычленения из области модели I ранга прямоугольного блока, состоящего из 50 * 80 ячеек (рис. 2). После этого каждая ячейка делится на 16 частей (п = 4). Таким образом, общее количество ячеек в одном слое модели II ранга составляет 200 * 320.
Рис. 3. Область (а) и слои (б) модели II ранга. Сплошной линией показана граница прогнозируемой области распространения фильтрата отходов в пластах-коллекторах на момент окончания эксплуатации полигона. Темным выделены водоупорные слои, светлым - водоносные горизонты. Для удобства восприятия на рисунке приводится только каждая пятая из 160 поверхностей раздела слоев
Модели III ранга описывают непосредственно при-фильтровые зоны отдельных нагнетательных скважин, используемых для захоронения нетехнологических отходов, и области прифильтровых зон групп скважин для захоронения технологических отходов. Таким образом, для каждой нагнетательной скважины нетехнологических отходов создается отдельная модель III ранга. Для участка захоронения технологических отходов формируются десять прискважиных моделей, опи-
сывающих прифильтровые зоны всех нагнетательных скважин, когда-либо эксплуатировавшихся на этом участке (рис. 4). При переходе к моделям III ранга каждая ячейка модели II ранга разбивается на 625 ячеек модели III ранга (п = 25). В результате площадь каждой из ячеек модели III ранга составляет 1 м2. Столь высокий уровень детализации обусловлен необходимостью описания отдельных элементов литологической толщи в непосредственной близости от фильтра скважины.
Рис. 4. Области прискважинных моделей на участке для захоронения технологических отходов.
На сетке модели II ранга темными прямоугольниками показаны области каждой из моделей III ранга. Точки соответствуют отдельным нагнетательным скважинам, кругами показаны области их прифильтровых зон
Иерархическая модель реализована в виде структурированной совокупности текстовых файлов. Для моделей всех рангов созданы отдельные директории с соответствующими названиями. В этих директориях располагаются папки, содержащие файлы, послойно описывающие распределение одного из параметров геологической среды. Имя папки соответствует названию этого параметра. Имена файлов содержат индексы описываемых слоев в соответствии с иерархической структурой данных. Кроме того, в каждой папке находится текстовый файл, описывающий ее внутреннюю структуру и внешние связи с другими файлами и папками модели.
Для визуализации структурного и статистического анализа моделей любого ранга могут использоваться широко распространенные пакеты программ, такие как картографическая система Surfer (версия 6.0 и выше) [10] или пакет программ для трехмерного конструирования AutoCAD (версия 2000І и выше). Также могут применяться специализированные программные средства для отображения геологических структур: 3D Groundwater Explorer, Stanford Geophysical Modeling Software [11]. Результаты визуализации могут быть экспортированы в другие пакеты программ для обработки графических данных растрового и векторного форматов. Выполнение геофильтрационных и геоми-грационных расчетов на основе построенных моделей возможно с использованием специализированных программных средств: Feflow [12], Processing Modflow
[13], Groundwater Modeling System [14]. Прогнозные расчеты состояния геологической среды с применением моделей всех рангов в составе единой иерархической модели возможны с использованием оригинального программного комплекса, разрабатываемого авторами в настоящий момент.
Разработанная иерархическая модель геологической среды полигона глубинного захоронения ЖРО СХК может быть использована для анализа литологостратиграфической структуры объектов различного пространственного масштаба с различным уровнем детальности. На основе данной модели возможно проведение расчетов физико-химических процессов, протекающих при взаимодействии ЖРО с пластовыми водами и вмещающими породами в прифильтровых зонах нагнетательных скважин. Кроме того, эта модель может использоваться для прогнозирования изменения состояния пласта-коллектора и распространения ореола загрязнения как в пределах территории полигона глубинного захоронения ЖРО, так и в целом для всего района. Предложенный метод построения иерархической модели может применяться для создания цифровых многоуровневых моделей геологической среды природно-антропогенных систем, возникающих при захоронении радиоактивных и токсичных отходов в осадочные геологические формации, эксплуатации артезианских водозаборов площадного и линейного типов, разработке осадочных месторождений полезных ископаемых геотехнологическими методами и т.д.
ЛИТЕРАТУРА
1. Рыбальченко А.И., ПименовМ.К., Костин П.П. и др. Глубинное захоронение жидких радиоактивных отходов. М.: ИздАТ, 1994. 256 с.
2. Зубков А.А., Макарова О.В., Данилов В.В. и др. Техногенные геохимические процессы в песчаных пластах-коллекторах при захоронении
жидких радиоактивных отходов // Геоэкология. 2002. № 2. С. 133-144.
3. Истомин А.Д., Носков М.Д., Кеслер А.Г. и др. Моделирование распределения радионуклидов в пласте-коллекторе при глубинном захороне-
нии кислых жидких радиоактивных отходов // Радиохимия. 2007. Т. 49, № 2. С. 182-187.
4. Pozdniakov S.P., Bakshevskay V.A., Zubkov A.A. et al. Modeling of Waste Injection in Heterogeneous Sandy Clay Formations // Underground Injec-
tion Science and Technology / Еd. by C.F. Tsang and J. Apps. Elsevier, 2005. Р. 203-218.
5. Зубков А.А., Данилов В.В. Носков М.Д., Истомин А.Д. Прогнозное моделирование распространения фильтрата ЖРО в пластах-коллекторах
полигона глубинного захоронения ЖРО СХК // Геоинформатика. 2007. № 4. С. 36-43.
6. Alexandrova L.N., Glinsky M.L., Danilov V.V. et al. Prediction of Radioactive Liquid Waste Migration at the Siberian Chemical Combine Under-
ground injection Site Based on Mathematical Modeling // Abs. of Sec. International symposium «Underground injection science and technology». Berkeley, 2003. 77 p.
7. Слинко А.М., Тер-Саакян Н.Ф., Зильберштейн Б.М. и др. Оценка эксплуатационных запасов подземных вод для хозяйственно-питьевого
водоснабжения г. Северска Томской области методом математического моделирования. Отчет ФГУГП «Гидроспецгеология». М., 2004. 187 с.
8.МатеронЖ. Основы прикладной геостатистики. М.: Мир, 1968. 408 с.
9. Всеволожский В.А. Ресурсы подземных вод южной части Западно-Сибирской низменности. М.: Наука, 1973. 88 с.
10. Surfer. URL: http://www.goldensoftware.com
11. Stanford Geophysical Modeling Software (S-GeMS). URL: http://sgems.sourceforge.net/
12. Finite Element subsurface FLOW system (Feflow). URL: http://www. scintificsoftwaregroup.com
13. Processing Modflow Pro. URL: http://www.webtech360.com
14. Groundwater Modeling System (GMS) for the PC. URL: http://www.ems-i.com/gms/downloads.html Статья представлена научной редакцией «Науки о Земле» 10 сентября 2009 г.
