Метод определения фильтрационных характеристик многократно подработанного массива

Побудована дифузійно-дісперсійна модель розподілу забруднюючих речовин у багатократно підробленому масиві гірських порід, яка адекватно описує реальні геомеханічні процеси. Шляхом розв’язання оберненої задачі розсіювання, були одержані коефіцієнти фільтрації багатократно підробленого масиву гірських...

Ausführliche Beschreibung

Gespeichert in:
Bibliographische Detailangaben
Datum:2015
Hauptverfasser: Дрибан, В.А., Дуброва, Н.А.
Format: Artikel
Sprache:Russian
Veröffentlicht: Український науково-дослідницький і проектно-конструкторський інститут гірничої геології, геомеханіки і маркшейдерської справи НАН України 2015
Schriftenreihe:Наукові праці УкрНДМІ НАН України
Online Zugang:http://dspace.nbuv.gov.ua/handle/123456789/99582
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
Назва журналу:Digital Library of Periodicals of National Academy of Sciences of Ukraine
Zitieren:Метод определения фильтрационных характеристик многократно подработанного массива / В.А. Дрибан, Н.А. Дуброва // Наукові праці УкрНДМІ НАН України. — 2015. — № 15. — С. 192-208. — Бібліогр.: 14 назв. — рос.

Institution

Digital Library of Periodicals of National Academy of Sciences of Ukraine
id irk-123456789-99582
record_format dspace
spelling irk-123456789-995822016-05-01T03:02:24Z Метод определения фильтрационных характеристик многократно подработанного массива Дрибан, В.А. Дуброва, Н.А. Побудована дифузійно-дісперсійна модель розподілу забруднюючих речовин у багатократно підробленому масиві гірських порід, яка адекватно описує реальні геомеханічні процеси. Шляхом розв’язання оберненої задачі розсіювання, були одержані коефіцієнти фільтрації багатократно підробленого масиву гірських порід у Центральному районі Донбасу. Diffusion-dispersion model for distribution of pollutants in repeatedly undermined rocks is constructed that satisfactorily describes real geomechanical processes. By solving inverse dispersion problem filtration factors for repeatedly undermined rocks in the Central area of Donbass were obtained. 2015 Article Метод определения фильтрационных характеристик многократно подработанного массива / В.А. Дрибан, Н.А. Дуброва // Наукові праці УкрНДМІ НАН України. — 2015. — № 15. — С. 192-208. — Бібліогр.: 14 назв. — рос. 1996-885X http://dspace.nbuv.gov.ua/handle/123456789/99582 622.02:539.2/.8 ru Наукові праці УкрНДМІ НАН України Український науково-дослідницький і проектно-конструкторський інститут гірничої геології, геомеханіки і маркшейдерської справи НАН України
institution Digital Library of Periodicals of National Academy of Sciences of Ukraine
collection DSpace DC
language Russian
description Побудована дифузійно-дісперсійна модель розподілу забруднюючих речовин у багатократно підробленому масиві гірських порід, яка адекватно описує реальні геомеханічні процеси. Шляхом розв’язання оберненої задачі розсіювання, були одержані коефіцієнти фільтрації багатократно підробленого масиву гірських порід у Центральному районі Донбасу.
format Article
author Дрибан, В.А.
Дуброва, Н.А.
spellingShingle Дрибан, В.А.
Дуброва, Н.А.
Метод определения фильтрационных характеристик многократно подработанного массива
Наукові праці УкрНДМІ НАН України
author_facet Дрибан, В.А.
Дуброва, Н.А.
author_sort Дрибан, В.А.
title Метод определения фильтрационных характеристик многократно подработанного массива
title_short Метод определения фильтрационных характеристик многократно подработанного массива
title_full Метод определения фильтрационных характеристик многократно подработанного массива
title_fullStr Метод определения фильтрационных характеристик многократно подработанного массива
title_full_unstemmed Метод определения фильтрационных характеристик многократно подработанного массива
title_sort метод определения фильтрационных характеристик многократно подработанного массива
publisher Український науково-дослідницький і проектно-конструкторський інститут гірничої геології, геомеханіки і маркшейдерської справи НАН України
publishDate 2015
url http://dspace.nbuv.gov.ua/handle/123456789/99582
citation_txt Метод определения фильтрационных характеристик многократно подработанного массива / В.А. Дрибан, Н.А. Дуброва // Наукові праці УкрНДМІ НАН України. — 2015. — № 15. — С. 192-208. — Бібліогр.: 14 назв. — рос.
series Наукові праці УкрНДМІ НАН України
work_keys_str_mv AT dribanva metodopredeleniâfilʹtracionnyhharakteristikmnogokratnopodrabotannogomassiva
AT dubrovana metodopredeleniâfilʹtracionnyhharakteristikmnogokratnopodrabotannogomassiva
first_indexed 2025-07-07T08:12:27Z
last_indexed 2025-07-07T08:12:27Z
_version_ 1836975074080456704
fulltext Наукові праці УкрНДМІ НАН України, № 15, 2015 Transactions of UkrNDMI NAN Ukraine, № 15, 2015 192 УДК 622.02:539.2/.8 МЕТОД ОПРЕДЕЛЕНИЯ ФИЛЬТРАЦИОННЫХ ХАРАКТЕРИСТИК МНОГОКРАТНО ПОДРАБОТАННОГО МАССИВА Дрибан В. А., Дуброва Н. А. (УкрНИМИ НАНУ, г. Донецк, Украина) Побудована дифузійно-дісперсійна модель розподілу забруд- нюючих речовин у багатократно підробленому масиві гірських порід, яка адекватно описує реальні геомеханічні процеси. Шля- хом розв’язання оберненої задачі розсіювання, були одержані ко- ефіцієнти фільтрації багатократно підробленого масиву гірсь- ких порід у Центральному районі Донбасу. Diffusion-dispersion model for distribution of pollutants in re- peatedly undermined rocks is constructed that satisfactorily describes real geomechanical processes. By solving inverse dispersion problem filtration factors for repeatedly undermined rocks in the Central area of Donbass were obtained. В ранее опубликованных авторами работах [1-3] установле- но, что тектонические нарушения разбивают массив горных по- род на крупные блоки, экранирующие распространение поллю- тантов, либо наоборот, задающие направление пространственно- го распределения загрязнителей. Для выявления пространствен- ных закономерностей характера распределения различных ВЗВ внутри тектонических блоков предложен следующий методиче- ский подход интерпретации исходных данных о распределении ВЗВ и построения безразмерной модели, позволяющей оценить общие характеристики процесса распределения загрязнителей в массиве. Суть данного методического подхода заключается в следующем: полученные статические поэлементные модели были нарезаны горизонтальными и вертикальными сечениями с интер- Наукові праці УкрНДМІ НАН України, № 15, 2015 Transactions of UkrNDMI NAN Ukraine, № 15, 2015 193 валом 100-200 м. Затем, на каждый горизонтальный и вертикаль- ный срез статистической модели были вынесены основные направления координатной сетки. Координатная сетка в каждом конкретном случае ориентирована таким образом (рис. 1), чтобы точка пересечения координатных осей находилась в центре "окружности", характеризующей максимальные концентрации загрязнителя на данной глубине (для горизонтальных сечений) и на данном направлении (для вертикальных). Часть моделей была нарезана в направлении простирания основных тектонических структур, задающих направление простирания ореолов загрязне- ния ВЗВ [4] и формирующих краевые условия распределения. Рис. 1. Горизонтальное и вертикальное сечение модели (за- грязнитель О-ксилол) Использование средств моделирования RockWorks2002 поз- волило визуализировать тело загрязнения в объеме массива гор- ных пород и интерполировать значение концентрации загрязни- теля в любой точке массива независимо от пространственного расположения места отбора проб. Благодаря этому, на опреде- ленных ранее направлениях горизонтальных и вертикальных се- Наукові праці УкрНДМІ НАН України, № 15, 2015 Transactions of UkrNDMI NAN Ukraine, № 15, 2015 194 чений, были выделены численные значения концентраций ВЗВ. Данные значения отбирались через равный интервал, соответ- ствующий шагу изолиний, отображающих концентрации соот- ветствующего загрязнителя (например, для О-ксилола шаг равен 0,05 – см. рис. 1) с пространственной привязкой данной точки от- носительно центра очага максимальных концентраций для кон- кретного среза модели. Отобранные значения концентраций ис- следуемых загрязнителей были нормированы относительно мак- симального значения, встретившегося именно на конкретном сре- зе [4, 5]. В качестве пространственного показателя распределения ВЗВ использован параметр, характеризующий удаленность точки с известной концентрацией загрязнителя от центра очага загряз- нения (от пространственного положения максимального значе- ния) в единицах, приведенных к размеру радиуса очага макси- мально обнаруженной концентрации загрязнителя на данном конкретном срезе. Затем, все расстояния, характеризующие удаленность точки с фиксированной концентрацией загрязнителя от центра очага максимальной его концентрации были приведены к значению Rmax для каждого конкретного среза, то есть нормированы отно- сительно размеров очага максимальных концентраций. Графиче- ское пояснение предложенной методики приведено на рисунке 2. На основании полученных данных были построены графики зависимости нормированных концентраций поллютантов от про- странственной составляющей их обнаружения, приведенной к размеру очага максимальных концентраций для каждого рас- сматриваемого среза модели. Таким образом, мы получили без- размерную модель распределения ВЗВ, позволяющую оценить общие характеристики процесса распространения. Всего было построено около пятидесяти графиков по четы- рем основным загрязнителям (О-ксилол, бутилацетат, стирол, моноэтилфосфат), с числом горизонтальных и вертикальных се- чений от 2-х до 4-х. Наукові праці УкрНДМІ НАН України, № 15, 2015 Transactions of UkrNDMI NAN Ukraine, № 15, 2015 195 а) бутилацетат; б) стирол Рис. 2. Графическое пояснение использованной методики Так для О-ксилола срезы модели выполнены на горизонтах +20, -100, -200 и -400 м в четырех основных направлениях света (С, Ю, З, В); для бутилацетата на горизонтах +120, 0 и -400 м в направлениях С, З, В; для стирола на горизонтах +120, 0 м в направлениях С, З, В; для моноэтилфосфата на горизонтах +200 и +100 м в направлениях С, Ю, З, В. Полученные для каждого за- грязнителя зависимости были разнонаправленно сгруппированы, т.е. рассмотрены совместно как в плоских сечениях (на одном го- ризонте во всех направлениях – рис. 3), так и в вертикальных (в одном направлении на разных горизонтах – рис. 4). В работах [6-9] был разработан и реализован инкременталь- ный подход для временной оценки трансформаций полей напря- жений, деформаций и проницаемости массива горных пород в процессе подработки, который является базисом для оценки из- менения фильтрационных характеристик массива горных пород условиях многократной систематической долговременной подра- ботки массивов. Было установлено, что долгосрочная многократ- ная подработка угольных пластов в условиях крутого падения приводит к образованию в массиве горных пород чередующихся пространственно-сопряженных зон сжатий и растяжений, заме- Наукові праці УкрНДМІ НАН України, № 15, 2015 Transactions of UkrNDMI NAN Ukraine, № 15, 2015 196 щающих друг друга в процессе ведения горных работ и генери- рующих квазирегулярную сеть вертикально ориентированных техногенных коллекторов, которые представляют собой области повышенной проницаемости, возникшие в результате воздей- ствия предельных и запредельных горизонтальных деформаций растяжения, и являются миграционными коридорами, способ- ствующими проникновению поллютантов на более глубокие го- ризонты. Рис. 3. Графики распределения ВЗВ по плоским сечениям Наукові праці УкрНДМІ НАН України, № 15, 2015 Transactions of UkrNDMI NAN Ukraine, № 15, 2015 197 Рис. 4. Графики однонаправленных погоризонтных распре- делений ВЗВ Наукові праці УкрНДМІ НАН України, № 15, 2015 Transactions of UkrNDMI NAN Ukraine, № 15, 2015 198 Данное обстоятельство позволило сформулировать гипоте- зу, о том, что указанные трансформации преобразовали массив горных пород в квазиоднородную по фильтрационным характе- ристикам среду. Для подтверждения адекватности полученной оценки изменения фильтрационных характеристик вмещающих пород и выдвинутой гипотезы, представляется возможным ис- пользовать экспериментальные данные о распределении в массиве вредных загрязняющих веществ (ВЗВ), как наиболее представительные, полные и удовлетворяющие задаче исследования. В качестве примера, иллюстрирующего интенсивность тех- ногенной нагрузки на исследуемом участке приведен рисунок 5 [5], где на разрезе вкрест простирания горных пород совместно визуализированы зоны влияния аппроксимированных контуров горных работ и очаг загрязнения массива ацетоном. Согласно [10], в сравнительно однородных породах, гомо- генных (квазигомогенных) по фильтрационным свойствам, ми- грация описывается в рамках моделей микродисперсии: предпо- лагается, что все механизмы процесса идут на одном микро- уровне, отвечающем репрезентативному элементу гомогенной пористой или трещиноватой среды. Для трещиноватых пород это предполагает: 1) пренебрежимо малое влияние пористости породной мат- рицы; 2) рассмотрение объемов массива, удовлетворяющих усло- вию сплошности среды; 3) синхронность заполнения веществом трещин разного по- рядка в пределах физической точки репрезентативного мини- мального объема среды. В работе [10] отмечено, что для крупноблочных трещинова- тых пород или пород с каналовым механизмом миграции послед- ние два условия выполнимы лишь для больших пространственно- временных масштабов, до достижения которых такие породы должны рассматриваться как гетерогенная среда. Наукові праці УкрНДМІ НАН України, № 15, 2015 Transactions of UkrNDMI NAN Ukraine, № 15, 2015 199 Р ис . 5 . Р аз ре з вк ре ст п ро ст ир ан ия г ор ны х по ро д с ви зу ал из ац ие й зо н вл ия ни я ап пр ок си - м ир ов ан ны х ко нт ур ов г ор ны х ра бо т и оч аг а за гр яз не ни я ац ет он ом Наукові праці УкрНДМІ НАН України, № 15, 2015 Transactions of UkrNDMI NAN Ukraine, № 15, 2015 200 Заметим, что описанные в [10] принципы, в рамках решае- мой нами задачи практически выполнены, что дает нам возмож- ность построения аналитической модели распространения ВЗВ в первом приближении. Таким образом, задачу можно формализо- вать следующим образом. Требуется определить поле распределения концентраций ВЗВ в плоскости при наличии постоянно действующего во вре- мени источника загрязнений интенсивностью 1, занимающего некоторую область . В общем случае диффузионно-дисперсионный перенос опи- сывается уравнением, базирующимся на феноменологическом за- коне Фика [11-14]: 0s c n J W t       , (1) где функция J(c) отвечает массовому потоку (статистически осредненному в пределах минимальных репрезентативных объе- мов пористой или трещиноватой среды); Ws – объемная интенсивность поглощения-выделения, обу- словленная проявлением внутренней гетерогенности пород (двойная пористость, регулярная слоистость и т.п.) или физико- химическими взаимодействиями. Вид функционала J(c) вообще говоря, зависит от мерности процесса. В общем трехмерном случае в декартовой системе ко- ординат описывается формулой:  i ij i i j c J c D x x x               , (2) где Dij – компоненты тензора микродисперсии, зависящие от рассеивающих свойств среды. Высказанные выше соображения о трансформации много- кратно подработанного массива в квазигомогенную среду позво- ляют существенно упростить выражения (1, 2). В безразмерных координатах указанные формулы принимают вид: Наукові праці УкрНДМІ НАН України, № 15, 2015 Transactions of UkrNDMI NAN Ukraine, № 15, 2015 201 2 2 2 2x y            , (3) где  - обобщённый параметр увязывающий время протека- ния процесса, диффузионно-фильтрационные характеристики массива и геометрические параметры области источника. На основании частного решения уравнения диффузии (фильтрации, теплопроводности) поле распределений нормирован- ных концентраций поллютанта можно записать в следующем виде:   2 0 1 4 0 0 , r rt kr t e d d            , (4) где r0 – радиус-вектор данной точки плоскости сечения. Анализ фактических распределений загрязнителей (рис. 1, 2) позволяет, по крайней мере, в первом приближении аппрокси- мировать область максимальных концентраций в виде круга ра- диуса R. В этом случае выражение (4) упрощается и принимает приемлемый для счета и анализа одномерный вид:   2 02 1 1 4 0 0 0 0 , r rt kr t r e drd d             , (5) где r,  - полярные координаты. На рисунке 6 представлены типовые кривые распределения концентраций для различных значений обобщенного параметра в зависимости от относительного (приведенного к радиусу области максимальных концентраций) расстояния. Итак, построена теоретическая модель, описывающая поля распределения поллютантов. Следующим необходимым шагом является верификация указанной модели по экспериментальным данным. Для этого были проанализированы функции невязки ти- повых распределений и фактических данных:   2 ( ) ,i i i f x yf    , (6) Наукові праці УкрНДМІ НАН України, № 15, 2015 Transactions of UkrNDMI NAN Ukraine, № 15, 2015 202 где yfi – фактические относительные концентрации; xi – относительные расстояния до центра пятна загрязнения;  – обобщенный параметр, определяемый по формуле 24ktR  . Рис. 6. Типовые кривые распределения концентраций для различных значений обобщенного параметра в зави- симости от относительного расстояния В связи с тем, что достоверно не известно время начала уте- чек ВЗВ из поверхностных источников загрязнения, расчет про- изводился для временного интервала длительностью десять лет. На рисунке 7 представлена кривая невязки для О-ксилола направление восток-запад гор. минус 120 м. Очевидно, что минимум указанной функции соответствует наилучшему приближению теоретических кривых к эксперимен- тальным данным. Заметим, что высокие корреляционные отно- шения свидетельствуют о хорошем приближении построенной тео- ретической модели протекающим геомеханическим процессам. В соответствии с данными принципами, были произведены перерасчеты и построены планиметрические модели распростра- нения ВЗВ, а также решена обратная задача: по данным имитаци- Наукові праці УкрНДМІ НАН України, № 15, 2015 Transactions of UkrNDMI NAN Ukraine, № 15, 2015 203 онного моделирования найдены параметры диффузионно- дисперсионного процесса распределения поллютантов, описыва- ющие параметры указанного процесса, протекающего в массиве. Рис. 7. Кривая невязки для О-ксилола (направление восток- запад, горизонт минус 120 м) Так как мы рассматриваем процесс квазистационарного со- стояния, при котором утечки носят постоянный во времени дол- госрочный характер, то параметр времени t, не вносит принципи- альных ошибок в полученные результаты, что связано с времен- ной структурой решения уравнения диффузии. Также заметим, что источник загрязнения не является точечным, а представляет собой некую локализованную область. Следует отметить, что перерасчет выполнялся совместно для сопряженных направлений (восток-запад, север-юг), однако в связи со спецификой исходных данных, некоторые направления рассмотрены одиночно. Таким образом, в соответствии с вышеизложенными прин- ципами, были произведены перерасчеты и построены планимет- рические модели погоризонтного распределения для каждого за- грязняющего вещества. Примеры полученных моделей представ- лены на рисунках 8-10. Наукові праці УкрНДМІ НАН України, № 15, 2015 Transactions of UkrNDMI NAN Ukraine, № 15, 2015 204 Рис. 8. Планиметрическая модель распределения моноэтил- фосфата на горизонте +100 м Рис. 9. Планиметрическая модель распределения О-ксилола на горизонте -120 м Наукові праці УкрНДМІ НАН України, № 15, 2015 Transactions of UkrNDMI NAN Ukraine, № 15, 2015 205 Рис. 10. Типовые кривые распределения концентраций О-ксилола на горизонте минус 120 м, направление запад, восток Сводные данные о полученных результатах погоризонтного и поэлементного диффузионно-дисперсионного распределения приведены в таблице 1. Полученные результаты (см. табл. 1) показывают высокие корреляционные соотношения (от 0,68 до 0,9 и выше), что свиде- тельствует о том, что выбранная модель верна. Необходимо отметить, что полученные результаты не толь- ко подтверждают факт того, что массив квазиоднороден, но и весьма проницаем. А поскольку вода является главным «распро- странителем, носителем, переносчиком» поллютантов в массиве, то появляется возможность отождествить диффузионно- дисперсионные коэффициенты с коэффициентами фильтрации и проницаемости. Выводы. 1. Построена диффузионно-дисперсионная модель распре- деления ВЗВ адекватно описывающая реальные геомеханические процессы в многократно подработанном массиве горных пород. 0,0 0,2 0,4 0,6 0,8 1,0 1,2 0 1 2 3 4 5 6 7 8 9 10 О-кс, отметка -120 тренд запад запад восток тренд восток Наукові праці УкрНДМІ НАН України, № 15, 2015 Transactions of UkrNDMI NAN Ukraine, № 15, 2015 206 Таблица 1 Параметры погоризонтного диффузионно-дисперсионного распределения поллютантов В ЗВ , г ор из он т, на пр ав ле ни е Р ад иу с пя тн а с м ак си м ал ьн ой ко нц ен тр ац ие й, м К ол ич ес тв о да нн ы х Зн ач ен ие п ри к от о- ро м д ос ти га ет ся м ин им ум н ев яз ки К оэ ф ф иц ие нт ди ф ф уз ии (ф ил ьт ра ци и) К ор ре ля ци он но е от но ш ен ие С ре дн ек ва др ат ич - но е от кл он ен ие О-кс., З -120 122 14 5 6,132×10-5 0,976 3,916×10-3 О-кс., В -120 122 14 1 2,219×10-5 0,999 1,751×10-4 О-кс., Ю -120 94 8 4 3,027×10-5 0,993 3,719×10-4 О-кс., С -120 94 6 28 1,992×10-4 0,937 1,825×10-3 О-кс., В-З -200 50 27 17 3,473×10-5 0,910 0,013 О-кс., Ю -200 51 8 13 2,821×10-5 0,998 8,868×10-5 О-кс., С -200 51 8 133 2,78×10-4 0,881 5,990×10-3 О-кс., В-З -400 15 12 1550 2,99×10-4 0,820 6,913×10-3 О-кс., Ю -400 88 6 13 8,323×10-5 0,997 9,168×10-5 О-кс., С -400 88 4 94 5,906×10-4 0,982 2,243×10-4 О-кс., В +20 246 10 8 4,278×10-4 0,983 1,436×10-3 О-кс., З +20 246 10 42 2,021×10-3 0,993 6,037×10-4 О-кс., Ю +20 271 5 18 1,107×10-3 0,985 3,124×10-4 О-кс., С +20 271 5 54 3,165×10-3 0,845 2,921×10-3 Бут.-ац., З-В 0 104 16 23 2,031×10-4 0,684 0,01 Бут.-ац., С 0 83 11 7 4,028×10-5 0,998 9,331×10-5 Бут.-ац., ЗВ-400 54 10 80 1,879×10-4 0,834 2,217×10-3 Бут.-ац., ЗВ+120 352 10 6 6,247×10-4 0,897 9,982×10-4 Бут.-ац.,С +120 250 10 14 7,17×10-4 0,928 2,610×10-3 Стирол, З 0 69 10 60 2,265×10-4 0,743 0,019 Стирол, В 0 69 9 13 4,907×10-5 0,729 0,016 Стирол, С 0 87 10 4 2,904×10-5 0,831 0,013 Стирол, В +120 367 7 12 1,318×10-3 0,905 3,697×10-3 Стирол, З +120 367 7 27 2,942×10-3 0,840 6,005×10-3 Стирол, С +120 298 7 6 4,427×10-4 0,975 9,879×10-4 MEP, З +100 174 6 24 5,936×10-4 0,711 5,646×10-3 MEP, В +100 174 6 9 2,381×10-4 0,959 9,161×10-4 MEP, С-Ю +100 70 6 35 1,403×10-4 0,917 5,688×10-4 MEP, З-В +200 111 12 146 1,441×10-3 0,795 5,939×10-3 MEP, С-Ю +200 162 12 56 1,176×10-3 0,685 6,563×10-3 Наукові праці УкрНДМІ НАН України, № 15, 2015 Transactions of UkrNDMI NAN Ukraine, № 15, 2015 207 2. Путем решения обратной задачи рассеяния, были получе- ны коэффициенты фильтрации многократно подработанного мас- сива горных пород. 3. Предложена методика обработки данных результатов фактических наблюдений о распределении ВЗВ в массиве, позво- ляющая привести тело рассеяния ВЗВ к безразмерному виду и по относительно малому объему фактических данных восстанавли- вать фильтрационные параметры массива, то есть решать прямую и обратную задачу распределения поллютантов. СПИСОК ССЫЛОК 1. Дуброва Н. А. Влияние тектонических структур на ореолы рассеяния загрязняющих веществ в подработанном массиве горных пород / Дуброва Н. А., Дьяченко Н. А. // Наукові праці УкрНДМІ: зб. наук. праць. – Донецьк, УкрНДМІ НАН Украї- ни. – 2011. – №9. – С. 485 – 501. 2. Дуброва Н. А. Распространение вредных загрязняющих ве- ществ в подработанном массиве горных пород со сложной те- ктонической нарушенностью (на примере ЦРД) /Н. А. Дуброва // Форум гірників – 2012: матеріали міжнар. конф., м. Дніпропетровськ.– Д.: Національний гірничий уні- верситет,– 2012. – Т. 1.– С. 166-172. 3. Дуброва Н. А. Закономерности распределения вредных загря- зняющих веществ в сложнонарушенном массиве горных по- род (на примере ЦРД) / Н. А. Дуброва // Проблеми гірничої технології: зб. матеріалів регіональної науково-практичної конференції., Красноармійськ. – 2012. – С. 310- 317. 4. Дрибан В. А. Изучение количественных и пространственных параметров распределения вредных загрязняющих веществ в техногенно нарушенном массиве горных пород сложного строения / Дрибан В. А., Дуброва Н. А. // Наукові праці Укр- НДМІ: зб. наук. праць. – Донецьк, УкрНДМІ НАН України. – 2011. – № 11. – С. 306-318. 5. Дуброва Н. А. Оценка изменения фильтрационных свойств деформируемого массива в результате многократной подра- ботки / Н. А. Дуброва // Вісник КрНУ. – Кременчук: КДПУ, 2014. – Вип. 1/2014 (84). – C. 96-105. Наукові праці УкрНДМІ НАН України, № 15, 2015 Transactions of UkrNDMI NAN Ukraine, № 15, 2015 208 6. Дрибан В. А. Особенности деформирования массивов горных пород ЦРД при многократной долгосрочной подработке / Дрибан В. А., Грищенков Н. Н., Дуброва Н. А. // Форум гір- ників – 2013: матеріали міжнар. конф., м. Дніпропетровськ.– Д.: Національний гірничий університет,– 2013. – Т. 2. – С. 226-232. 7. Дрибан В. А. Особенности формирования техногенных кол- лекторов при отработке угольных пластов в условиях ЦРД / Дрибан В. А., Грищенков Н. Н., Ходырев Е. Д., Дуброва Н. А. // Наукові праці УкрНДМІ: зб. наук. праць. – Донецьк, УкрН- ДМІ НАН України. – 2011. – №13 (I). – С. 220-237. 8. Driban V., Dubrova N. Evaluation of changes in rock mass per- meability due to long-time repeated mining // PROGRESSIVE TECHNOLOGIES OF COAL, COALBED METHANE, AND ORES MINING. – CRC Press/Balkema, Netherlands, 2014. – pp. 167-174. 9. Driban, V. Measurement of change of rock mass permeability in consequence of long-time repeated mining / V. Driban, N. Dubrova // Материалы научно-практической конференции посвященной 80-летию стахановского движения. Луганский государственный университет им. Владимира Даля, – Лу- ганск, 2015. – С.72-78. 10. Мироненко В. А. Проблемы гидрогеоэкологии. Т. 1. Теорети- ческое изучение и моделирование геомиграционных процес- сов / В. А. Мироненко, В. Г. Румынин . – М.: Изд. Московско- го государственного горного университета, 1998. – 611 с. 11. Бочевер Ф. М. Защита подземных вод от загрязнения / Боче- вер Ф. М., Лапшин Н. Н, Орадовская А. Е. – М.: Недра, 1979. – 245 с. 12. Лукнер Л. Моделирование миграции подземных вод / Лук- нер Л., Шестаков В. М. – М.: Недра. – 1986. 13. Bear J. Dynamics of fluids in porous media. Amer. Els. N. Y., 1972. – 764 p. 14. Шестаков В. М. Динамика подземных вод / В. М. Шестаков – М.: изд. МГУ, 1979. – 368 с.