Помощь в написании студенческих работ
Антистрессовый сервис

Обоснование эффективных параметров прогнозных моделей геомиграции в неоднородной водовмещающей среде

ДиссертацияПомощь в написанииУзнать стоимостьмоей работы

Ключевым аспектом для надежного моделирования массопереноса является правильное изображение неоднородности (см. главу 1). Однако при моделировании массопереноса на региональном масштабе возникает проблема, связанная с тем, как учесть локальную мелкомасштабную неоднородность. Во-первых, на практике для описания структуры неоднородности необходимо нереально большое количество данных. Во-вторых… Читать ещё >

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

Содержание

  • 1. Анализ физико-математических моделей миграции в неоднородных средах
    • 1. 1. Классическое конвективно-дисперсионное уравнение
    • 1. 2. Признаки и причины аномального переноса
    • 1. 3. Моделирование массопереноса в неоднородных средах
      • 1. 3. 1. Подход, использующий КДУ и эмпирическую макродисперсивность
      • 1. 3. 2. Стохастический подход к моделированию массопереноса в неоднородных средах
      • 1. 3. 3. Подход, использующий осреднение и КДУ на микромасштабе
      • 1. 3. 4. Альтернативные подходы
    • 1. 4. Проблемы параметрического обеспечения долгосрочных моделей миграции компонентов ЖРО
    • 1. 5. Выводы
  • 2. Методы моделирования геологической неоднородности осадочных отложений
    • 2. 1. Классификация методов построения моделей геологической неоднородности
    • 2. 2. Гауссовские геостатистические методы
    • 2. 3. Негауссовские геостатистические методы
    • 2. 3. Детерминированные и стохастические геостатистические методы
    • 2. 4. Выводы
  • 3. Ремасштабирование коэффициента фильтрации
    • 3. 1. Эффективная и эквивалентная проницаемость
    • 3. 2. Методы ремасштабирования
    • 3. 3. Выводы
  • 4. Геолого-гидрогеологические условия района СХК
    • 4. 1. Краткие сведения о СХК
    • 4. 2. Геологическое строение района СХК
    • 4. 3. Гидрогеологические условия района
    • 4. 4. Условия осадконакопления верхнемеловых и эоценовых отложений района СХК
    • 4. 5. Выводы
  • 5. Геостатистический анализ данных
    • 5. 1. Анализ исходных данных
    • 5. 2. Обоснование параметров для геостатистического моделирования
      • 5. 2. 1. Параметры для двухкомпонентной модели литологической неоднородности
      • 5. 2. 2. Параметры для четырехкомпонентной модели литологической неоднородности
  • 5.
  • Выводы
  • 6. Построение моделей геофильтрационной неоднородности и оценка эффективных геофильтрационных параметров
    • 6. 1. Построение трехмерных моделей геофильтрационной неоднородности водовмещающих отложений
    • 6. 2. Определение эффективной проницаемости неоднородных песчано-глинистых толщ
    • 6. 3. Выводы
  • 7. Калибрация модели и моделирование истории закачки ЖРО
    • 7. 1. Калибрация модели геофильтрационной неоднородности
    • 7. 2. Крупномасштабное численное моделирование миграции нейтрального компонента
    • 7. 3. Выводы
  • 8. Оценка эффективных геомиграционных параметров
    • 8. 1. Методика
    • 8. 2. Моделирование
    • 8. 3. Результаты и обсуждение
    • 8. 3. Влияние схематизации неоднородности осадочных отложений на прогноз миграции загрязнения
    • 8. 4. Выводы

Актуальность исследования. Одной из актуальных проблем современной России является защита среды обитания человека от отходов предприятий химической, атомной промышленности и энергетики. Одним из способов длительной изоляции токсичных отходов, включая жидкие радиоактивные отходы, является их захоронение в глубокозалегающие геологические формации. В 60-х годах в России были созданы полигоны глубинного захоронения жидких радиоактивных отходов (ПГЗ ЖРО). Самым крупным из них является полигон Сибирского химического комбината (СХК) (г.

Северск), пласты-коллекторы которого сложены неоднородными песчаноглинистыми породами. В соответствии с действующими в настоящее время нормативами допускается захоронение ЖРО на существующих подземных хранилищах при условии локализации отходов в границах установленных горных отводов [СанПиН 2.6.1.07−03]. Учитывая большие объемы уже захороненных на.

— 1 полигоне СХК к настоящему времени ЖРО (более 43 млн. м общей активностью около 1 млрд. Ки на глубину 270 — 400 м от поверхности земли [Румынии, 2011]) и то, что эксплуатация полигона продолжается, актуальной является оценка безопасности последствий подземного захоронения ЖРО. Для такой оценки необходимы численные модели, с помощью которых выполняются долговременные прогнозы миграции компонентов ЖРО в подземной гидросфере. Ключевым аспектом при моделировании геомиграции контаминантов является учет влияния геофильтационной неоднородности [Шестаков, 2003; Мироненко, Румынии, 1998; Riva et al., 2008]. При создании прогнозных моделей геомиграции в неоднородной водовмещающей среде, возникают следующие проблемы: (1) высокая изменчивость водовмещающих отложений исключает возможность получения полного знания о геофильтрационных свойствах среды, в которой происходит массоперенос- (2) локальная неоднородность не может быть явно представлена в модели, имеющей размеры ячейки в сотни метров. Поэтому возникает необходимость в разработке методик оценки эффективных параметров для прогнозных моделей. В полной мере все вышесказанное относится и к моделям, созданным и разрабатываемым в настоящее время для ПГЗ ЖРО СХК. За последнее десятилетие было разработано несколько численных моделей, которые описывают геологическую среду в районе СХК с различным уровнем детальности, причем одними из главных проблем является отсутствие согласованности между моделями различного масштаба и параметрическое обеспечение геомиграционного моделирования. Поэтому построение трехмерной локальной модели геофильтрационной неоднородности водовмещающих отложений ПГЗ, и определение на ее основе эффективных параметров для прогнозных моделей необходимо для успешного решения задач, связанных с безопасной эксплуатацией полигонов глубинного захоронения жидких отходов.

Целью работы является обоснование трехмерных моделей геофильтрационной неоднородности и разработка методики оценки эффективных параметров для прогноза миграции загрязнения в неоднородной песчано-глинистой водовмещающей среде.

Для достижения поставленной цели решались следующие задачи:

1) Анализ существующих физико-математических моделей переноса загрязнения в подземных водах.

2) Обзор методов моделирования геологической неоднородности и обоснование оптимального метода построения трехмерной модели неоднородных песчано-глинистых отложений.

3) Анализ и обработка данных по ПГЗ ЖРО СХК.

4) Разработка и обоснование трехмерной детальной модели геофильтрационной неоднородности ПГЗ ЖРО СХК и ее калибрация на основе полевых данных.

5) Крупномасштабное численное моделирование миграции нейтрального компонента ЖРО.

6) Разработка методики оценки эффективных параметров прогнозных моделей на основе трехмерных моделей геофильтрационной неоднородности водовмещающих отложений с высоким разрешением.

7) Оценка возможности применимости упрощенных моделей массопереноса.

Объектом исследования является неоднородная водовмещающая среда ПГЗЖРО СХК.

Предметом исследования является геомиграция нейтрального компонента ЖРО в неоднородной водовмещающей среде.

Фактический материал и методы исследования. Работа основана на фактических данных лаборатории геотехнологического мониторинга СХК, полученных в процессе разведки и эксплуатации ПГЗ ЖРО, которые включают данные детального литологического расчленения разрезов 295 скважин, данные по гранулометрическому составу 1117 проб песка из 66 скважин, данные замеров уровней подземных вод в наблюдательных скважинах СХК за 1995 — 1996 г. г., данные по расходам скважин СХК с 1967 по 2003 г. г. Для решения поставленных задач были проанализированы отечественные и зарубежные публикации по вопросам физико-математических моделей массопереноса в неоднородных средах и создания трехмерных (3D) моделей геологической неоднородности. В работе использовались статистические методы обработки данных, метод геостатистического моделирования для построения 3D модели геофильтрационной неоднородности водовмещающих отложений, численный метод для моделирования геофильтрации и геомиграции. Для выполнения и оформления работы использовались следующие стандартные и специализированные пакеты программ: Microsoft Office 2003, Surfer 8.0, Grafer 7.0, STATISTICA, MODFLOW 2000, MT3DMS, T-PROGS, SgemS, Adobe Photoshop 9.0.

Научная новизна исследования:

1) Разработана методика моделирования геофильтрационной неоднородности, на основе которой построена новая 3D модель геофильтрационной неоднородности песчано-глинистых водовмещающих отложений ПГЗ ЖРО СХК пригодная для оценки эффективных параметров и оценки риска быстрой вертикальной миграции загрязнения.

2) Предложена методика оценки эффективных параметров для прогнозных геофильтрационных и геомиграционных моделей.

3) Выявлены основные закономерности миграции загрязнения в неоднородной песчано-глинистой среде на основе созданной детальной 3D модели неоднородности и показана ограниченность применения стандартной модели конвективно-дисперсионного переноса для условий ПГЗ ЖРО СХК.

Личный вклад: обработка и комплексный анализ данных, построение моделей, многовариантное моделирование геофильтации и геомиграции компонента ЖРО.

Практическая значимость исследования:

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

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

3) Эффективные параметры, полученные с помощью ЗБ модели геофильтрационной неоднородности водовмещающих отложений ПГЗ СХК, использованы для выполнения среднесрочных и долгосрочных прогнозов распространения компонентов ЖРО в районе СХК.

По результатам проведенных исследований сформулированы следующие защищаемые положения:

1) Для определения эффективных геофильтрационных и геомиграционных параметров прогнозных моделей неоднородных песчано-глинистых водовмещающих отложений необходимо обоснование и создание ЗЭ модели геофильтрационной неоднородности с высоким разрешением. В качестве основного метода построения такой модели целесообразно использовать метод вероятностей перехода, основанный на цепях Маркова.

2) Геостатистический анализ данных детального литологического расчленения разрезов скважин показал, что водовмещающие отложения района ПГЗ ЖРО СХК существенно анизотропны. Основные литологические разности (песчаная и глинистая) имеют примерно одинаковую долю в разрезе и отношение горизонтального характерного размера неоднородности к вертикальному больше 47, что обусловливает значительную фильтрационную анизотропию водовмещающих отложений: вертикальный эффективный коэффициент фильтрации на 2 — 3 порядка меньше, чем горизонтальный.

3) Модель для долгосрочного прогноза миграции компонентов ЖРО в районе СХК должна строиться не в виде пластовой системы, согласно которой моделируемая толща состоит из водоносных пластов, разделенных водоупорами, а в виде анизотропно-слоистой среды, состоящей из слоев с примерно одинаковыми горизонтальными коэффициентами фильтрации и низкими значениями вертикального коэффициента фильтрации. Подобная модель менее консервативна, чем пластовая, и именно она должна использоваться для выполнения оценок долговременной миграции отходов.

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

Апробация результатов работы. Основные результаты исследования и положения диссертации докладывались и обсуждались на российских и международных конференциях: Первой Сибирской международной конференции молодых ученых по наукам о Земле (Новосибирск, 2002), Международной научной конференции студентов, аспирантов и молодых ученых «Ломоносов-2003» (Москва, 2003), конференции «Ломоносовские чтения"-2003 (Москва, 2003), IV межвузовской молодежной конференции «Школа экологической геологии и рационального недропользования-2003» (Санкт-Петербург, 2003), Международной конференции «Underground Injection Science and Technology» (California, 2003), а также на научно-практическом семинаре «Обоснование безопасности пунктов захоронения РАО на базе современных моделей переноса радионуклидов в геологических формациях и инженерных барьерах» (Москва, 2011).

Публикации. Основные положения работы изложены в 12 публикациях, в том числе в трех статьях в журналах, рекомендованных ВАК.

Структура и объем диссертации

Работа изложена на 194 страницах, состоит из введения, восьми глав и заключения, содержит 46 рисунков и 24 таблицы.

Список использованных источников

включает 183 наименования (из них 57 отечественных и 126 зарубежных).

7.3 Выводы.

Для обоснования применимости модели геофильтрационной неоднородности и с целью уточнения фильтрационных и емкостных параметров была проведена калибрация двухкомпонентной модели фильтрационной неоднородности. Для калибрации использовались данные замеров уровней подземных вод во II, III и IV водоносных горизонтах в наблюдательных скважинах промплощадки СХК за 1995 — 1996 г. г. Методика калибрации локальной модели по данным напоров может быть использована для генерирования начальных параметров прогнозной модели (например, распределения концентраций, полученные на 3D детальной среднемасштабной модели, можно использовать в качестве начальных концентраций для региональной модели). Результаты показали хорошее.

149 воспроизведение изменений напора подземных вод в процессе закачки даже для двухфациальной модели среды, что свидетельствует о соответствии построенной модели геологической среде ПГЗ СХК.

На основе детальной трехмерной высокоразрешимой двухкомпонентной модели водовмещающих отложений промплощадки 18 и, используя данные по расходам нагнетательных скважин, было проведено крупномасштабное численное моделирование миграции нейтрального компонента с 1967 по 2003 г. г. Рассматривалась миграция нейтрального компонента (например, нитрат-иона, присутствующего во всех видах ЖРО СХК) как наиболее жесткий вариант с максимально возможным распространением загрязнения в геологической среде, т. к. из результатов полевых экспериментов и литературных данных известно, что существует механизм быстрой миграции (в малых концентрациях) радионуклидов со скоростью нейтральных компонентов раствора. При моделировании учитывались локальная продольная и поперечная микродисперсивность, а также молекулярная диффузия. Результаты моделирования миграции загрязнения в песчано-глинистых отложениях показали, что ореол загрязнения имеет весьма сложную форму в пространстве в зависимости от распределения (трехмерной топологии) в пласте песчаных и глинистых тел, контролирующих флуктуации поля скоростей и, как следствие, дисперсию нейтрального компонента. Через 36 лет от начала закачки ореол загрязнения с безразмерной концентрацией 0,001 — 1 находится в пределах промплощадки 18, причем загрязнение фиксируется ниже подошвы относительного водоупора Э не смотря на его прерывистое литологическое строение. Главная масса отходов находится внутри песчаных пропластков, несмотря на длительный (36 лет) период миграции. Заполнение глин в результате диффузии и адвекции играет малую роль в транспорте отходов.

8. Оценка эффективных геомиграционных параметров.

Ключевым аспектом для надежного моделирования массопереноса является правильное изображение неоднородности (см. главу 1). Однако при моделировании массопереноса на региональном масштабе возникает проблема, связанная с тем, как учесть локальную мелкомасштабную неоднородность. Во-первых, на практике для описания структуры неоднородности необходимо нереально большое количество данных. Во-вторых, локальная мелкомасштабная неоднородность не может быть явно представлена в модели регионального масштаба (см. раздел 3.3). Однако возможно следующее решение проблемы: с помощью геостатистического моделирования ТР/МС методом создание 3D локальной модели геофильтрационной неоднородности водовмещающих отложений и затем с ее помощью получение эффективных геомиграционных параметров, которые будут использоваться в прогнозной модели массопереноса. Целью главы является определение эффективных геомиграционных параметров для неоднородной водоносной системы.

8.1 Методика.

Метод прослеживания случайного блуждания частиц (в англоязычной литературе — random walk particle tracking method или RWPT) рассматривает транспорт контаминанта как перенос большого числа движущихся в фильтрационном потоке частиц. К преимуществам этого метода относятся отсутствие численной дисперсии и вычислительная эффективность, т.к. время вычисления для моделей с очень большим количеством ячеек и характеризующихся сильной неоднородностью оказывается значительно меньше, чем в традиционном Эйлеровском методе или в TVD схемах [Zhang et al., 2006; Salomon et al., 2006]. Поэтому метод RWPT часто используется для моделирования транспорта в сложных неоднородных средах [LaBolle, Fogg, 2001; Weissmann et al, 2002; Fernandez-Garcia et al, 2005; Engdahl, Weissmann, 2010].

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

1) построение 3D детальной модели геофильтрационной неоднородности водовмещающих отложений с помощью ТР/МС метода (см. раздел 6.1);

2) решение стационарной фильтрационной задачи на полученной 3D модели (см. раздел 6.2);

3) моделирование конвективного переноса нейтрального компонента из неточечного источника на полученной модели с помощью RWPT метода;

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

— получение выходных кривых и оценка с их помощью коэффициента макродисперсии для горизонтального и вертикального переноса.

Значения действительной скорости фильтрации U и коэффициента макродисперсии D могут быть получены из временных моментов экспериментальных выходных кривых [Jury, 1990; Yu et. al., 1999; Fernandez-Garcia et al., 2005c]:

U = -k <8Л>

8.2).

2x где x [м] - расстояние между местом запуска частиц и контрольной плоскостью, ц2 — второй центральный временной момент (дисперсия случайной переменной, которая описывает среднеквадратичное отклонение значений вокруг среднего значения) который может быть определен как [Yu et. al., 1999; Cirpka, 2000]: f0°°(t — т)2 С (х, i) dt т2 ц2(х) =—ъг———-=—— = М2 — М{ (8.3).

0 С (х, t) dt rnQ т0/ где mn — это временной момент n-порядка случайной переменной: рсо тп = I tnC (x, t) dt (8.4).

•^о где C (x, t) — локальная концентрация на расстоянии х, t — время. Мп — нормализованный момент: тп С tn С (х, t) dt.

Мп = —- = ^—(8.5) т0 Cc (x, f) dt и цп — центральный момент: it-M^ctx.m lin =———(8.6).

0 C{x, t) dt.

Расчетная пористость определялась по формуле:

V U где средняя скорость фильтрации (полученная по результатам геофильтрационного neff = - (8.7) моделирования в разделе 6.2) V = кеГТ х / = —, где I, С2 и со — градиент, расход и со площадь поперечного сечения фильтрационного потока.

Продольная дисперсивность аь рассчитывалась по формуле: а, = § (8.8).

Эта методика была опробована на примере неоднородных песчано-глинистых водовмещающих отложений ПГЗ ЖРО СХК.

8.2 Моделирование.

Размеры четырехкомпонентной модели геофильтрацигонной неоднородности приведены в разделе 6.1. Активная пористость задавалась равной 0,2 для всей области моделирования. Использовали четыре варианта моделей, построенные для наиболее вероятных границ горизонтальных размеров гидрофаций Ь (первые четыре варианта из табл. 6.4 раздела 6.2):

1) 2) с Ls =143 м (4 гидрофации) — вариант с минимальными размерамис Ls =243 м (4 гидрофации) — базовый вариант;

3) с Ls =443 м (4 гидрофации) — вариант с максимальными размерами;

4) с Ls=243 м (2 гидрофации) — базовый вариант без промежуточных гидрофаций.

Первые три варианта состоят из четырех гидрофаций (песка, глинистого песка, песчаной глины и глины), последний — из двух гидрофаций (песка и глины).

Для моделирования транспорта и для вычисления траектории и времени движения частиц в конвективном потоке использовалась программа РМРАТН [Chiang, Kinzelbach, 2001]. Входной информацией для РМРАТН являлось распределение скоростей фильтрации во всех узлах модели, полученное в процессе моделирования фильтрации на MODFLOW 2000;2005 [Harbaugh et al., 2000] с помощью многосеточного метода (Geometric Multigrid Solver) и начальные координаты точек. Выходной информацией — координаты х, у, z и время прибытия каждой частицы на противоположную границу (контрольную плоскость).

В каждом из 8 численных экспериментов 4500 частиц помещались на границу с большим значением напора в ячейки модели, имеющие значение коэффициента фильтрации больше 0,002 м/сут, либо на нижнюю грань параллелепипеда (во второй снизу слой для вертикальной миграции), либо на правую грань (во вторую справа колонку для горизонтальной миграции). Очевидно, что при одной и той же модельной сетке и значениях напоров на противоположенных гранях, распределение времен прибытия частиц зависит только от распределения коэффициентов фильтрации внутри моделируемой области. Для дальнейшего расчета брались данные только по тем частицам, которые прошли расстояние равное расстоянию до противоположной грани модели (расстояние в 4275 м — для горизонтального потока и 129 м — для вертикального), и затем делался пересчет позволяющий перейти от времен прибытия для каждой частицы к безразмерным концентрациям.

8.3 Результаты и обсуждение.

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

Показать весь текст

Список литературы

  1. , В. С. Массоперенос в водонасыщенных горных породах / В. С. Алексеев, Г. М. Коммунар, Б. С. Шержуков. Итоги науки и техники. ВИНИТИ. Серия: Гидрогеология, инженерная геология. Том 11. — М.: Информационное издание ВИНИТИ, 1989. — 145 с.
  2. , В. Г. Математическое моделирование радиационно-химического разложения органических примесей жидких радиоактивных отходов при глубинном захоронении / В. Г. Балахонов, А. А. Зубков, В. А. Матюха // Радиохимия. 2001. — Т.43. — № 1. — С. 82−86
  3. , А. В. Вопросы внедрения методов многоточечной статистики /
  4. А. В. Безруков, М. В. Рыкус, А. Р. Давлетова, В. И. Савичев // Научно-технический вестник ОАО «НК «Роснефть». 2010.- № 2. — С. 2−8
  5. , Г. К. Основы теории изменчивости инженерно-геологических свойств пород / Г. К. Бондарик. М.: Недра, 1971. — 272 с.
  6. , Ф. М., Орадовская, А. Е. Гидрогеологическое обоснование защиты подземных вод и водозаборов от загрязнения / Ф. М. Бочевер, А. Е. Орадовская. -М.: Недра, 1972. 128 с.
  7. , А. А., Сурков, В. С. Геологическое строение и перспективы нефтегазоносности юго-восточной части Западно-Сибирской низменности. / А. А. Булынникова, В. С. Сурков М.: Гостоптехиздат, 1962. — 73 с.
  8. , В. А. Ресурсы подземных вод южной части ЗападноСибирской низменности / В. А. Всеволожский. М.: Наука, 1973. — 88 с.
  9. , И. К., Ковалевский, В. С., Язвин, Л. С. и др. Основы гидрогеологии. Гидрогеодинамика / И. К. Гавич, В. С. Ковалевский, Л. С. Язвин и др. -Новосибирск: Наука, 1983. 238 с.
  10. Геология СССР, том XIV, Западная Сибирь (Алтайский край, Кемерово, Новосибирская, Омская, Томская области). Часть 1. Геологическое описание / Под ред.: В. Д. Фомичев, И. Н. Звонарев, Л. Д. Староверов. М.: Недра, 1967. — 664 с.
  11. , И. И. Введение в теорию случайных процессов / И. И. Гихман, А. В. Скороход. М.: Наука, 1977. — 568 с.
  12. , Л. Б. Позднемеловая флора Сибири: автореф. дисс. .д-ра биол. наук 03.00.05 / Л. Б. Головнева. СПб., 2004. — 35 с.
  13. , А. В. Основы региональной палеоклиматологии / А. В. Гольберт. М.: Недра, 1987.-224 с.
  14. , А. В. Палеоландшафты Западной Сибири в юре, мелу и палеогене / А. В. Гольберт, Л. Г. Маркова, И. Д. Полякова. М.: Наука, 1968. -144с.
  15. , В. В. Математическое моделирование глубинного захоронения жидких радиоактивных отходов (на примере Сибирского химического комбината): автореф. дисс. .канд геол.-минер.наук 25.00.07 / В. В. Данилов. Томск, 2010. — 23 с.
  16. , В. В. Геостатистика: теория и практика / В. В. Демьянов, Е. А. Савельева. М.: Наука, 2010. — 327 с.
  17. , Дж. С. Статистический анализ данных в геологии / Дж. С. Дэвис. -М.: Недра, 1990.-319 с.
  18. , О. Использование геостатистики для включения в геологическую модель сейсмических данных / О. Дюбрул. 2е1з1: ЕАвЕ, 2002. — 296 с.
  19. , О. Геостатистика в нефтяной геологии / О. Дюбрул. Москва-Ижевск: Институт компьютерных исследований, НИЦ «Регулярная и хаотическая динамика», 2009. — 256 с.
  20. Западная Сибирь / Геология и полезные ископаемые России. В 6 томах. / Ред. А. Э. Конторович. СПб.: Изд-во ВСЕГЕИ, 2000. Т. 2. — 477с.
  21. , И. М. Ресурсы пресных и маломинерализованных подземных вод южной части Западно-Сибирского артезиансткого бассейна / И. М. Земскова, Ю. К. Смоленцев, М. П. Полканов. М., Недра, 1991. — 262 с.
  22. Г. М. Оценка неоднородности и прогноз нефтеизвлечения по ГИС / Г. М. Золоева М.: Недра. 1995. — 212 с.
  23. , М. М. Прогнозирование эффективных толщин в межскважинном пространстве: методология, тенденции, оценка результатов / М. М. Кавун,
  24. A. В. Степанов, П. В. Ставинский // Геофизика.- 2008, — № 4.- С. 17−21.
  25. , Г. Н., Корчебоков, Н. А., Разин, К. И. Движение подземных вод в неоднородных пластах / Г. Н. Каменский, Н. А. Корчебоков, К. И. Разин. -М.: Объединенное научно-техническое издательство, 1935. 169 с.
  26. , М. Ф., Демьянов, В. В. и др. Основные понятия и элементы геостатистики / М. Ф. Каневский, В. В. Демьянов, Е. А. Савельева, С. Ю. Чернов // Проблемы окружающей среды и природных ресурсов. Обзорная информация. -1999.-Вып. 11.-С. 15−21.
  27. , Н. П., Величкин, В. И., Омельяненко, Б. И., Петров,
  28. B. А., Тарасов, Н. Н. Новые подходы к подземному захоронению высокоактивных отходов в России / Н. П. Лаверов, В. И. Величкин, Б. И. Омельяненко, В. А. Петров, Н. Н. Тарасов // Геоэкология. 2000. — № 1. — С. 3−12.
  29. , Н. П., Омельяненко, Б. И., Величкин, В. И. Геоэкологические аспекты проблемы захоронения радиоактивных отходов / Н. П. Лаверов, Б. И. Омельяненко, В. И. Величкин // Геоэкология. 1994. — № 6. — С. 3−20.
  30. , Л. Моделирование миграции подземных вод / Л. Лукнер, В. М. Шестаков. М.: Наука, 1986. — 208 с.
  31. , В. А. Изучение загрязнения подземных вод в горнодобывающих районах / В. А. Мироненко, Е. В. Мольский, В. Г. Румынии. -Л.: Недра, 1988.-279 с.
  32. , В. А. Проблемы гидрогеоэкологии. Монография в 3-х томах. Том 2. Опытно-миграционные исследования / В. А. Мироненко, В. Г. Румынии. -М.: Издательство Московского государственного горного университета, 19 986. -394 с.
  33. Объяснительная записка к атласу литолого-палеогеографических карт юрского и мелового периодов Западно-Сибирской равнины масштаба 1:5 000 000 // Под. ред. И. И. Нестерова. Тюмень: ЗапСибНИГНИ, 1976. (Труды ЗапСибНИГНИ- Вып. 93). -276 с.
  34. , О. В. Особенности геологического моделирования продуктивных пластов флювиального происхождения // О. В. Пинус, К. В. Пайразян. Геология нефти и газа. — 2008. — № 1. — С. 25−30
  35. , В. М. Фораминиферы биостратиграфия верхнего мела и палеогена в Западной Сибири / В. М. Подобина. Томск: Томский государственный университет, 2009. — 432 с.
  36. , С. П., Бакшевская, В. А. Численный анализ макродисперсии при закачке мигранта в водоносный пласт со слабопроницаемыми включениями /
  37. С. П. Поздняков, В. А. Бакшевская / Проблемы гидрогеологии XXI века. Наука и образование. М.: Изд-во РУДН, 2003. — С. 156−170.
  38. , В. К. Формирование и эксплуатация подземных вод Обь-Томского междуречья / В. К. Попов, Г. М. Рогов, О. Д. Лукашевич. Томск: Изд-во ТАСУ, 2002.- 173 с.
  39. Рац, М. В. Неоднородность горных пород и их физических свойств / М. В. Рац. М.: Наука, 1968. — 108 с.
  40. Рац, М. В. Структурные модели в инженерной геологии / М. В. Рац М.: Недра, 1973.-215 с.
  41. , В. Г. Геомиграционные модели в гидрогеологии / В. Г. Румынии. СПб.: Наука, 2011. — 1157 с.
  42. , А. И. Глубинное захоронение жидких радиоактивных отходов / А. И. Рыбальченко, М. К. Пименов, П. П. Костин. М.: ИздАТ, 1994. -256 с.
  43. , И. В. Стратиграфия и палеогеография верхнего мела — среднего палеогена юго-востока Западной Сибири: автореф. дисс. .канд геол.-минер.наук 25.00.02 / И. В. Рычкова. Томск, 2006. — 22 с.
  44. Стратиграфия СССР. Меловая система. В 2-х полутомах / Под ред. М. М. Москвина. М., Недра, 1987. — Полутом 2. — 326 с.
  45. , И. В. Оценка долгосрочной безопасности захоронения радиоактивных отходов. 1. Палеореконструкция условий формирования подземных вод / И. В. Токарев, А. А. Зубков, В. Г. Румынии // Водные ресурсы. 2009. — Т. 36. -№ 2.-С. 219−227
  46. , Е. В. Геолого-геофизическая модель Северской площади / Е. В. Черняев // Известия Томского политехнического университета. Т.305. Вып. 6. -2002. — С.414−433
  47. , М. И. Статистическая гидрогеодинамика пористых сред / М. И. Швидлер. М.: Недра, 1985. — 288 с.
  48. , В. М. Гидрогеодинамика / В. М. Шестаков. М.: Изд-во МГУ, 1995. — 368 с.
  49. , В. М. Учет геологической неоднородности — ключевая проблема гидрогеодинамики / В. М. Шестаков // Вестник МГУ. Сер. 4. Геология. 2003.-№ 1.-С. 25−28
  50. , JI. С. Достоверность гидрогеологических прогнозов при оценке эксплуатационных запасов подземных вод / Л. С. Язвин. М., 1972. 168 с.
  51. Benson, D. Application of a fractional advection-dispersion equation / D. Benson, S. Wheatcraft, M. Meerschaert // Water Resources Research. American Geophysical Union, 2000. — Vol. 36. — N. 6 — P. 1403−1412.
  52. Berkowitz, В., Scher, H. The role of probabilistic approaches to transport theoryin heterogeneous media / B. Berkowitz, H. Scher // Transport in porous media. 179
  53. Springer Science+Business Media B.V., Formerly Kluwer Academic Publishers B.V., 2001. Vol.42.-P. 241−263.
  54. Burns, E. R. Upscaling facies models to preserve connectivity of designated facies / E. R. Burns, L. R. Bentley, R. Therrien, C. V. Deutsch // Hydrogeology Journal. Springer-Verlag GmbH, 2010. — Vol.18. — P. 1357−1373.
  55. Carle S. F. T-PROGS: Transition Probability Geostatistical Software. Version 2.0. Davis, California: University of California, 1998. — 76 p.
  56. Carle S. F. High-resolution simulation of basin-scale nitrate transport considering aquifer system heterogeneity / S. F. Carle, B. K. Esser, J. E. Moran // Geosphere. Geological Society of America, 2006. — Vol. 2. — P. 195−209.
  57. Carle S. F., Fogg G. E. Modeling spatial variability with one and multidimensional Markov chains / S. F. Carle, G. E. Fogg // Mathematical Geology. -Springer, 1997.-Vol. 29, N7.-P. 891−918.
  58. Carle S. F., Fogg G. E. Transition probability-based indicator geostatistics / S. F. Carle, G. E. Fogg // Mathematical Geology. Springer, 1996. — Vol. 28, N 4. — P. 453 476.
  59. Carrera, J. Inverse problem in hydrogeology / J. Carrera, A. Alcolea, A. Medina, J. Hidalgo, L. Slooten // Hydrogeology Journal. Springer-Verlag GmbH, 2005. -Vol. 13.-P. 206−222.
  60. Chiang, W. H., Kinzelbach, W. 3D-Groundwater modeling with PMWIN. -First Edition. Berlin, Heidelberg, New York: Springer, 2001. 346 p.
  61. Cirpka, O.A., Kitanidis, P.K. Characterization of mixing and dilution in heterogeneous aquifers by means of local temporal moments / O. A. Cirpka, P. K. Kitanidis // Water Resources Research. American Geophysical Union, 2000. -Vol. 36(5). — p. 122−134
  62. Dagan, G. Solute transport in heterogeneous porous formations / G. Dagan // J. Fluid Mech. Cambridge University Press, 1984. — Vol. 145. — P. 151- 177.
  63. Dentz, M., Berkowitz, B. Transport behavior of a passive solute in continuous time random walks and multirate mass transfer / M. Dentz, B. Berkowitz // Water Resources Research. American Geophysical Union, 2003. — Vol. 39, N 5. — P. 11 111 131.
  64. Dentz, M. Time behavior of solute transport in heterogeneous media: transition from anomalous to normal transport / M. Dentz, A. Cortis, H. Sher, B. Berkowitz // Advances in Water Resources. Elsevier, 2004. — Vol. 27. — P. 155−173.
  65. Deutsch C., Journel A. GSLIB: Geostatistical Software Library and Usres Guide. Oxford University Press, New York, 1992. — 340 p.
  66. Doherty J. MODFLOW-ASP Using MODFLOW-2000 with PEST-ASP
  67. Watermark Computing. Australia, 2001. 156 p.
  68. Elfeki, A., Dekking, M. A markov chain model for subsurface characterization: theory and applications / A. Elfeki, M. Dekking // Mathematical Geology. Springer, 2001. — Vol. 33, N 5. P. 569−589.
  69. Engdahl, N., Weissmann, G. Anisotropic transport rates in heterogeneous porous media / N. Engdahl, Weissmann // Water Resour. Res. 2010. V. 46. — W02507
  70. Flach, G. P. Comparison of singl-domain and dual-domain subsurface transport models / G. P. Flach, S. A. Crisman, F J. Molz // Ground Water. Proquest Academic Research Library, 2004. — Vol.42, N 6. — P. 815−828.
  71. Fogg, G. E. Geologically based model of heterogeneous hydraulic conductivity in an alluvial setting / G. E. Fogg, C. D. Noyes, S. F. Carle // Hydrogeology Journal. -Springer-Verlag GmbH, 1998. Vol. 6. — P. 131−143.
  72. Gelhar, L. Stochastic subsurface hydrology / L. Gelhar Prentice-Hall, 1993. -390 p.
  73. Gelhar, L. Stochastic subsurface hydrology from theory to application / L. Gelhar // Water Resources Research. American Geophysical Union, 1986. — Vol. 22, N 9.-P. 135−145.
  74. Gelhar, L. W., Axness, C. L. Three-dimensional stochastic-analysis of macrodispersion in aquifers / L. W. Gelhar, C. L. Axness // Water Resources Research. -American Geophysical Union, 1983. Vol. 19, N 1. — P. 161- 180.
  75. Gelhar L. W. A critical review of data on field-scale dispersion in aquifers / L. W. Gelhar, C. Welty, K. R. Rehfeldt // Water Resources Research. American Geophysical Union, 1992. — Vol. 28. — P. 1955−1974.
  76. Godeke S. Simulation of a reactive tracer experiment using stochastic hydraulic conductivity fields / S. Godeke, H. Geistlinger, A. Fischer, H. Richnow, M. Schirmer // Environmental Geology. Springer-Verlag GmbH, 2008. — Vol. 55. — P. 1255−1261.
  77. Gomez-Hernandez J. J. Simulation of block permeability conditioned upon data measured at a different scale. / J. J. Gomez-Hernandez / ModelCARE 90: Calibration and Reliability in Groundwater Modelling. IAHS Publ. 1990. — N. 195. — P. 407−416.
  78. Gueguen Y. Upscaling: effective medium theory, numerical methods and the fractal dream / Y. Gueguen, M. Le Ravalec, L. Ricard // Pure appl. Geophys. 2006. Vol. 163.-P. 1175−1192.
  79. Haggerty R. On the late-time behavior of tracer test breakthrough curves / R. Haggerty, S. A. McKenna, L. C. Meigs // Water Resources Research. American Geophysical Union, 2000. — Vol. 36, N 12. — P. 3467−3479.
  80. Hill, M.C., Tiedeman, C.R., Effective Groundwater Model Calibration: With Analysis of Data, Sensitivities, Predictions, and Uncertainty: Wiley and Sons, 2007. 464 P
  81. Janza, M. Modelling heterogeneity of Ljubljana polje aquifer using Markov chain and geostatistics // Geologija. Ljubljana, 2009. — V. 52/2. — P. 233−240.
  82. Jones, N. L. Hydrogeologic Unit Flow Characterization Using Transition Probability Geostatistics / N. L. Jones, J. R. Walker, S. F. Carle // Ground Water. 2005. -V. 43 (2). P. 285−289
  83. Journel, A. G., Deutsch C. V. Entropy and spatial disorder / A. G. Journel, C. V. Deutsch // Mathematical Geology. Springer, 1993. — Vol. 25, N 3. — P. 329−355.
  84. Jury, W. A. Transfer functions and solute movement through soil: theory and applications. Boston, Berlin, Birkhauser: Birkhauser Verlag Basel, 1990. — 346 p.
  85. Kerrou, J. Issues in characterizing heterogeneity and connectivity in non-multiGaussian media / J. Kerrou, Ph. Renard, H. Franssen, I. Lunati // Advances in Water Resources. Elsevier Science Publishing Company, Inc., 2008. — Vol. 31. — P. 147−159.
  86. Levy, M., Berkowitz, B. Measurement and analysis of non-Fickian dispersion in heterogeneous porous media / M. Levy, B. Berkowitz // Journal of Contaminant
  87. Hydrology. Elsevier Science Publishing Company, Inc., 2003. — Vol. 64, N 3−4. — P. 203−226.
  88. Mariethoz, G. Truncated plurigaussian simulations to characterize aquifer heterogeneity / G. Mariethoz, P. Renard, F. Cornaton, O. Jaquet // Ground Water. -Proquest Academic Research Library, 2009. Vol. 47, N 1. — P. 13−24.
  89. Neuman, S. P. Universal scaling of hydraulic conductivities and dispersivities in geologic media / S. P. Neuman // Water Resources Research. American Geophysical Union, 1990. Vol. 26, N 8. — P. 1749−1758.
  90. Neuman, S. P., Di Federico, V. On advective dispersion in fractal velocity and permeability fields / S. P. Neuman, V. Di Federico // Water Resources Research. -American Geophysical Union, 2003. Vol. 31, N 6. — P. 1455−1460.
  91. Neuman, S. P., Tartakovsky, D. M. Perspective on theories of non-Fickian transport in heterogeneous media / S. P. Neuman, D. M. Tartakovsky // Advances in Water Resources Elsevier Science Publishing Company, Inc., 2009. — Vol. 32. — P. 670−680.
  92. Neuman, S. P. Stochastic theory of field-scale Fickian dispersion in anisotropic porous media / S. P. Neuman, C. L. Winter, C. M. Newman // Water Resources Research. American Geophysical Union, 1987. — Vol. 23, N 3. P. 453166.
  93. Novak, W. Geostatistical methods for the identification of flow and transport parameters in the subsurface / Institut fur Wasserbau, Universitat Stuttgart. Stuttgart: Inst. fiirWasserbau, 2005. — 150 p.
  94. Proce, C. J. Modeling Multiscale Heterogeneity and Aquifer Interconnectivity / C. J. Proce, R. W. Ritzi, D. F. Dominic, Z. Dai // Ground Water. Proquest Academic Research Library, 2004. — Vol. 42, N 5. — P. 658−670.
  95. Renard, Ph. Stochastic hydrogeology: What professionals really need? // Ground Water. Proquest Academic Research Library, 2007. — Vol. 45, N 5. — P.531−541.
  96. Renard, Ph., de Marsily, G. Calculating equivalent permeability: a review / Ph. Renard, G. de Marsily // Advances in Water Resources. Elsevier Science Publishing Company, Inc., 1997. Vol. 20, N 5−6. — P. 253−278.
  97. Ritzi, R. W. Behavior of indicator variograms and transition probabilities in relation to the variance in lengths of hydrofacies / R. W. Ritzi // Water Resources Research. American Geophysical Union, 2000. — Vol. 36, N 11. — P. 3375−3381.
  98. Ritzi, R. W. Comparing statistical models of physical heterogeneity in buried-valley aguifers / R. W. Ritzi, D. F. Dominic, A. J. Slesers // Water Resources Research. -American Geophysical Union, 2000. Vol. 36, N 11. — P. 3179−3192.
  99. Riva, M. Relative importance of geostatistical and transport models in describing heavily tailed breakthrough curves at tht Lauswiesen site / M. Riva, A. Guadagnini, Fernandez-Garcia et al. // J. of Contaminant Hydrology, v. 101, 2008. P 113.
  100. Rubin, Y. Spatial variability in river sediments and its link with river channel geometry / Y. Rubin, I. A. Lunt, J. S. Bridge // Water Resources Research. American Geophysical Union, 2006. — Vol 42. — W06D16 doi:10.1029/2005WR004853.
  101. Salamon, P. A review and numerical assessment of the random walk particle tracking method / P. Salamon, D. Fernandez-Garcia, J. Comez-Hernandez // J. of Contaminant Hydrology, v/ 87, 2006, 277−305.
  102. Salandin, P., Fiorotto, V. Solute transport in highly heterogeneous aquifers / P. Salandin, V. Fiorotto // Water Resources Research. American Geophysical Union, 1998. Vol. 34, N. 5. — P. 949−961.
  103. Sanchez-Vila, X. A synthesis of approaches to upscaling of hydraulic conductivities / X. Sanchez-Vila, J. P. Girardi, J. Carrera // Water Resources Research. -American Geophysical Union, 1995. Vol. 31, N. 4. — P. 867−882.
  104. Sanchez-Vila, X. Representative hydraulic conductivities in saturated groundwater flow / X. Sanchez-Vila, A. Guadagnini, J. Carrera // Reviews of geophysics. American Geophysical Union, 2006. — Vol. 44. RG 3002. 46 p.
  105. Schulze-Makuch, D. Longitudinal dispersivity data and implication for scaling behavior // Ground Water. Proquest Academic Research Library, 2005.- Vol. 43, N 3. -P. 443−456.
  106. Shestakov, V. M. Development of relationship between specific storage anddepth of sand and clay formation // Environmental Geology. 2002. — V.42,, N 2−3. -P.127−129.
  107. Sivakumar, B. A fractal investigation of solute travel time in a heterogeneous aquifer: transition probability/Markov chain representation / B. Sivakumar, T. Harter, H. Zhang // Ecological modelling. Elsevier, 2005. Vol. 182, N 3−4. — P. 355−370.
  108. Vidstrand, P. Comparison of upscaling methods to estimate hydraulic conductivity // Ground Water. Proquest Academic Research Library, 2001. — Vol. 39, N3.P. 401−407.
  109. Weissmann, G. S., Fogg, G. E. Multi-scale alluvial fan heterogeneity modeled with transitionprobability geostatistics in a sequence stratigraphic framework // Journal of Hydrology. Elsevier Science Publishing Company, Inc., 1999. — Vol. 226. P. 48−65.
  110. Weissmann, G. S. Dispersion of groundwater age in an alluvial aquifer system / G. S. Weissmann, Y. Zhang, E. M. LaBolle, G. E. Fogg // Water Resources Research. -American Geophysical Union, 2002. Vol. 38, N 10. — P. 1198−1206.
  111. Willmann, M. Transport upscaling in heterogeneous aquifers: What physical parameters control memory functions? / M. Willmann, J. Carrera, X. Sanchez-Vila // Water Resources Research. American Geophysical Union, 2008. doi: 10.1029/2007WR006531.
  112. Willmann, M. Transport upscaling in a clay-rich sedimentary formation using multi-porosity models / M. Willmann, W. Kinzelbach, F. Stauffer // Model care 2009. -P.417−420.
  113. Yu, C. A moment method for analyzing breakthrough curves of step inputs / C. Yu, A. W. Warrick, M.H. Conklin // Water Resour. Res. 1999. — V. 11. — p. 35 673 572.
  114. Zhang, Y. Upscaling conductivity and porosity in three-dimensional heterogeneous porous media // Chinese Science Bulletin. 2004. Vol.49, N 22. P. 24 152 423.
  115. Zhang, Y., Benson, D., Baeumer, B. Predicting the tails of breakthrough curves in regional scale alluvial systems // Ground Water. — National Ground Water Association, 2007. — Vol. 45, N 4. — P. 473−484.
  116. Zheng, C., Gorelick, S. M. Analysis of solute transport in flow fields influenced by preferential flowpaths at the decimeter scale // Ground Water. Proquest Academic Research Library, 2003. — Vol. 41, N 2. — P. 142−155.
  117. Отчет по теме «Переоценка эксплуатационных запасов подземных вод Томского месторождения» / ОАО «Томскгеомониторинг», Томск, 2005 г. 436 с.
  118. А. А., Лехов А. В., Поздняков С. П., Рыбальченко А. И., Шестаков В. М. Отчет «Гидрогеологические исследования на объектах захоронения промышленных стоков» / М.: МГУ, 2003. 231 с.
  119. А. А. и др. Отчет «Результаты эксплуатации полигонов подземного захоронения жидких радиоактивных отходов (пл. 18, 18а) и наземных хранилищ жидких радиоактивных отходов в 1996 г.» / Северск, 1998. 179 с.
Заполнить форму текущей работой