CC BY
60
Мезо-, нано-, биомеханика и механика природных процессов Вестник Нижегородского университета им. Н.И. Лобачевского, 2011, № 4 (2), с. 505-507
УДК 539.3+622.83
НЕЛИНЕЙНЫЕ ПРОЦЕССЫ ЭВОЛЮЦИИ ГЕОМЕХАНИЧЕСКИХ ПОЛЕЙ ПРИРОДНЫХ И ТЕХНОГЕННЫХ ОБЪЕКТОВ
© 2011 г. Л.А. Назарова2, И.Н. Ельцов1, Л.А. Назаров1, М.И. Эповх
'Институт нефтегазовой геологии и геофизики СО РАН, Новосибирск 2Институт горного дела СО РАН, Новосибирск
larisa@misd.nsc.ru
Поступила в редакцию 15.06.2011
Методами механики деформируемого твердого тела описаны протекающие в породном массиве при ведении горных работ процессы необратимого деформирования и разрушения геосред. Предложен метод синтеза определяющих уравнений нарушения сплошности на основе величины фрактальной размерности его берегов. Построена геомеханическая модель процесса глубокого бурения, с использованием которой установлена зависимость размеров зон разрушения в прискважинной области от соотношения компонент природного поля напряжений. Разработан и апробирован на реальном объекте статистический подход, устанавливающий количественную связь между параметрами техногенной сейсмичности и поля напряжений в массиве горных пород при отработке месторождений полезных ископаемых. Сформулированы и решены обратные коэффициентные и граничные задачи геомеханики, позволяющие дать оценку времени жизни конструктивных элементов технологии выемки месторождений и фокальных параметров готовящегося сейсмического события по данным мониторинга состояния геосреды.
Ключевые слова: породный массив, напряжение, деформация, упругопластическая модель с дилатан-сией, эволюция геомеханических полей, скважина, обратная задача, техногенная сейсмичность.
Диагностика состояния массива горных пород по информации об эволюции параметров физических полей, вызванной природными и техногенными факторами, геомеханическое обоснование технологий отработки месторождений полезных ископаемых, оценка устойчивости подземных объектов различного масштабного уровня — вот далеко не полный перечень проблем, связанных с необходимостью описывать процессы необратимого деформирования и разрушения геосред.
Синтез уравнений состояния нарушений сплошности
Процессы необратимого деформирования структурированного массива горных пород локализуются, как правило, на межблочных границах, поэтому необходимо знать соотношения, описывающие закономерности взаимодействия берегов разломных нарушений. Разработаны методика генерации поверхностей с заданной фрактальной размерностью / и механизм взаимодействия двух полупространств, границы которых — случайные поверхности с одинаковой фрактальной размерностью. На основе статистической обработки результатов численных экспериментов установлена зависимость нор-
мального напряжения S(p) = Knp на нарушении от конвергенции берегов р, где нормальная жесткость Kn пропорциональна величине (f — 2)-0'5 и прочности на одноосное сжатие вмещающих пород. Предложен способ определения величины f по длине береговой линии.
Геомеханические аспекты
бурения глубоких скважин
Многие технологии горных работ связаны с созданием в породном массиве скважин различного назначения. При глубоком бурении для разведки и добычи углеводородов могут иметь место нежелательные явления такие, как отклонение скважины от заданной траектории, разрушение прискважинной зоны и изменение ее проницаемости. На основе упругопластической дилатантной модели горных пород [1] описан процесс формирования и эволюции зон необратимых деформаций D в окрестности глубокой скважины, проходящей пачку тонких слоев с пониженными прочностными свойствами — типичной ситуации при бурении скважин на нефть и газ. Оказалось:
— радиальные размеры пропорциональны величине (1 — q)2 (где q = Ch fcv — коэффициент бокового отпора, Oh и ov — горизонтальные и
вертикальные напряжения в нетронутом массиве);
— разрушение менее прочного слоя может наступить до пересечения его скважиной, что нужно учитывать при инверсии каротажных данных;
— перед забоем возникает область растягивающих радиальных напряжений, амплитуда которых пропорциональна разности давления бурового раствора и ch и может достигать предельных значений, что объясняет наблюдаемое в реальных условиях увеличение скорости бурения с глубиной [2];
— проницаемость зоны D может измениться до 15% в зависимости от глубины и прочностных свойств горных пород.
Связь техногенной сейсмичности и вариации параметров напряженного состояния при обработке месторождений полезных ископаемых
Предложен метод поэтапного решения краевых задач на основе иерархии объемных геомеханических моделей. Граничные условия на первом — глобальном уровне формулируются на основе косвенной (сейсмотектонической, геодезической) информации о полях напряжений в литосфере. На втором (региональном) и третьем (локальном) уровне для этой цели используются результаты расчетов с предыдущего иерархического уровня, уточняемые по данным измерений in situ параметров геомеханических полей. Реализация подхода выполнена с использованием метода конечных элементов (сферическая и декартова системы координат, контакт-элементы для описания межблочного взаимодействия) для объектов «Центральной Азия и ее обрамление», «Алтае-Саянская складчатая область», «район рудных месторождений Гор -ной Шории» и «Таштагольское железорудное месторождение». Для последнего объекта построена детальная геомеханическая модель, с использованием которой описан процесс эволюции напряженно-деформированного состояния при его отработке в 1978—2009 гг. Обоснован подход, позволяющий на основе статистического анализа характеристик пространственно-временного распределения очагов индуцированных горными работами динамических событий установить их количественную связь с параметрами напряженного состояния. Предложенный подход, апробированный с использованием базы данных сейсмических событий Таштагольского месторождения в 1989—
2009 гг., позволяет на основе планов горных работ и форвардных расчетов полей напряжений дать прогнозную количественную оценку уровня техногенной сейсмичности и локализации в пространстве очагов динамических явлений при отработке месторождений полезных ископаемых.
Обратные задачи геомеханики
При моделировании геомеханических полей в природных объектах различного масштабного уровня всегда присутствует элемент неопределенности, связанный, например, с недостатком данных о свойствах среды или действующих в породном массиве напряжениях. В таких случаях целесообразно сформулировать и решить обратную задачу соответствующего типа. В качестве входной информации могут выступать как прямые (результаты измерений деформаций, напряжений, смещений), так и косвенные (компоненты магнитного поля, сейсмотектонические деформации и их траектории) данные.
Граничные обратные задачи. Создан метод реконструкции краевых условий по информации о полях смещений или деформаций во внутренних точках области, заключающийся в построении дискретной функции Грина для соответствующей краевой задачи и последующей минимизации целевой функции. На этой основе предложен способ определения параметров эквивалентного точечного источника (типа двойной силы с моментом), моделирующего очаг готовящегося динамического события по изменению смещений на поверхности Земли, что позволяет дать априорную оценку координат гипоцентра и магнитуды предстоящего сейсмического события по GPS данным.
Коэффициентные обратные задачи. Предложен метод определения параметров уравнений состояния, описывающих реологические процессы деформирования и разрушения горных пород, базирующийся на решении обратных задач по данным измерений in situ изменения размеров подземных полостей в процессе отработки месторождения. Оказалось, что соответствующая целевая функция (среднее квадратичное отклонение модельных и измеренных величин конвергенции кровли и почвы очистного пространства) имеет несколько расположенных на одной прямой локальных минимумов, что позволило создать эффективный алгоритм решения обратных задач, комбинирующий методы градиентного спуска и прямого
поиска. На этой основе с использованием моделей ползучести и критерия накопления повреждений разработана методика оценки времени жизни целиков (типичных конструктивных элементов камерно-столбовой системы отработки месторождений полезных ископаемых) и последовательности их разрушения, что позволило обосновать новую технологию добычи медных руд Джезказганского месторождения.
Работа выполнена при поддержке РФФИ (гранты № 10-05-00736, 10-05-00835), Программы Президиума РАН №16.8 и Интеграционного Проекта СО РАН №60.
Список литературы
1. Николаевский В.Н. Геомеханика и флюидоди-намика. М.: Недра, 1996.
2. Калинин А.Г Бурение нефтяных и газовых скважин. М.: ЦентрЛитНефтеГаз, 2008.
NON-LINEAR PROCESSES OF GEOMECHANICAL FIELDS EVOLUTION IN NATURAL AND TECHNOLOGICAL OBJECTS
L.A. Nazarova, I.N. Eltsov, L.A. Nazarov, M.I. Epov
The paper uses methods of the deformable solid mechanics to describe mining-induced irreversible deformation and fracture processes in a rock mass, and proposes to construct constitutive relations of a discontinuity based on the value of fractal dimensionality of edges of the discontinuity. A deep-hole drilling model is constructed and used to relate the dimensions of ruptured zones in the well bore area and the ratio of natural stress components. The authors developed a statistic approach to determining the quantitative relation between parameters of the mining-induced seismic activity and stress field evolution in rocks, and have approved this approach in a real-life object. The formulated and solved inverse coefficients and boundary value problems allow estimating life time of design elements, which extraction technology includes, and the focal parameters of an incipient seismic event by the geomedium monitoring data.
Keywords: rock mass, stress, strain, elastoplastic model, dilatation, evolution of geo-mechanical fields, borehole, inverse problem, mining-induced seismic activity.
