Определение и верификация параметров расчетной модели соляных пород с учетом разупрочнения и ползучести

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

Для прогнозирования развития деформационных процессов во времени в массиве горных пород в окрестности подземного сооружения активно используются численные методы. Достаточно много специализированных моделей и методов решения задач механики горных пород и массивов, и подземного строительства реализованы в программных комплексах Abaqus, DIANA, FLAC, Midas, PLAXIS, ZSoil и др.

В нормативном документе СП 91.13330.2011 «Подземные горные выработки» приведены инженерные методики расчета долговременных деформаций и устойчивости подземных выработок, при этом никаких сведений о базовых моделях, лежащих в основе этих методик, не приводится. Важно отметить, что СП 91.13330.2011 в ряде случаев не позволяет произвести расчеты на устойчивость выработок глубокого залегания в соляных и калийных рудниках, а использование численных методов при проектировании – обязательно. Другие геотехнические своды правил рекомендуют использовать модели упруго-идеально-пластической среды Мора – Кулона и Хоека – Брауна, не приводя при этом конкретных рекомендаций по учету реологических эффектов, ограничиваясь лишь требованиями об их учете. На практике введение предположения об идеально пластическом поведении в ряде задач может привести к некорректным выводам как в части взаимодействия вмещающего массива с крепью [5], так и в части оценки размеров зон нарушения выбранных критериев или учета ползучести, особенно когда образование пластических зон больших размеров неизбежно. В таких случаях учет разупрочнения позволяет избежать указанной проблемы и корректнее проводить расчет деформаций ползучести за счет другого характера распределения напряжений в массиве. Кроме того, появляется возможность выполнения более точных расчетов на исследование случаев проявления прогрессирующего разрушения.

Постановка проблемы. Модели, построенные на основе интегрального оператора с Абелевым ядром [4], а также модели Мора – Кулона и Хоека – Брауна, отражающие идеально пластическое поведение, с эмпирической зависимостью типа Нортона [11] и подобными степенными зависимостями [2, 3] получили широкое распространение при изучении реономных эффектов. Вместе с тем наиболее корректно использовать описанные модели при определенных горно-геологических и горнотехнических условиях (малые глубины, заполнение выработанного пространства искусственными целиками и др.), либо при выполнении предварительного анализа. Но при проектировании сооружений на больших глубинах, создании сложных и крупномасштабных подземных геотехнических комплексов, инженерных подземных объектов повышенной опасности и высокой ответственности, зачастую результаты бывают недостаточно адекватными и корректными. Так, например, «сценарий» процесса разрушения массива горных пород в окрестности подземного сооружения, согласно модельным расчетам, может лишь отдаленно соответствовать реальным условиям [19, 26], когда процесс ползучести вмещающего массива описывается некорректно и отличается от данных мониторинга в натурных условиях [11, 23].

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

Методология. Для более корректного учета поведения соляных пород в сложных условиях нагружения предлагается использовать комбинированную вязко-упруго-пластическую модель, состоящую из модифицированной модели Бюргерса и модели Мора – Кулона, где основные параметры (сцепление \( c(\bar {\epsilon}_{pl})\), угол внутреннего трения \( \phi(\bar {\epsilon}_{pl})\), угол дилатансии \( \psi(\bar {\epsilon}_{pl})\) и предельные растягивающие напряжения \( \sigma_t(\bar {\epsilon}_{pl})\)) зависят от эффективных пластических относительных деформаций \( (\bar{\epsilon}_{pl})\) (рис.1). При этом полные относительные деформации \(\epsilon\) состоят из обратимой (упругой \( {\epsilon}_{el}\) и вязкоупругой \( {\epsilon}_{ve}\)) и необратимой (вязкой \( {\epsilon}_{visc}\) и пластической \( {\epsilon}_{pl}\)) частей, а изменение объемных относительных деформаций \( {\epsilon}_{v}\) происходит за счет упругой \( {\epsilon}_{v. el}\) и пластической \( {\epsilon}_{v. pl}\) частей:

Многими авторами предлагались подходы к построению модификаций модели Мора – Кулона с учетом деградации коэффициента сцепления и возрастания коэффициента трения. Основные предшествующие подходы заключались в учете линейного роста величины коэффициента трения, линейного уменьшения сцепления [13] и обеих величин [12], использования зависимостей экспоненциального характера [21, 25] или более сложных зависимостей с определением влияния напряженного состояния и введением логарифмических функций для угла дилатансии [18, 20]. При этом, например используются нелинейные зависимости для «мобилизованных» коэффициента сцепления c' и угла внутреннего трения φ', предложенные Вермеером и де Борстом [25], а также угла дилатансии ψ' согласно зависимости, предложенной Ченом с соавторами [18]:

где c, sin(φ') – максимальные значения сцепления и синуса угла внутреннего трения; \( {\epsilon}_{c}\), \( {\epsilon}_{f}\) – константы, влияющие на характер изменения кривых c' и sin(φ') соответственно; a, b, d, f – численные коэффициенты.

Обсуждение. На основании статистической обработки результатов лабораторных испытаний были получены модули упругости и спада, предел прочности при различных величинах обжатия образца в стабилометре σ_h, значения относительных продольных деформаций \( {\epsilon}_{3}\) для характерных точек диаграмм. Полученные на их основании значения модифицированной модели Мора – Кулона позволили воспроизвести усредненные условно-мгновенные испытания в стабилометре (кривые деформирования при разных величинах обжатия) (рис.2).

Для этапа вторичной ползучести, по данным лабораторных испытаний, выполнялось аналитическое определение параметров ядер ползучести для степенной и синус-гиперболической зависимостей и их сравнение. Установлено, что степенная зависимость не включает в доверительный диапазон напряжений значения выше 18,8 МПа, в то время как функция гиперболического синуса удовлетворительно «вписалась» в лабораторные испытания при всех значениях интенсивности напряжений (рис.3, 4). Аналогичная зависимость используется в варианте модификации модели Минкли [7, 10]). Окончательный вид используемой функции:

где \( {\dot\epsilon}_{visc}\) – скорость вязких относительных деформаций; a, b, c – коэффициенты; \( \sigma_i\) – интенсивность напряжений, \( \sigma_i=\sqrt{3j_2}\); \( \sigma_i^{'}\) – опорное давление.

Для верификации предложенной модели с использованием данных полевых наблюдений все наблюдательные станции делились на группы по схожим горно-геологическим и горнотехническим условиям и для каждой были проведены расчеты (результаты по одной из групп приведены на рис.5). Первоначальные расчеты с эмпирическим учетом влияния забоя (рис.5, линия «2D модель»), в которых на контуре выработки использовался понижающий множитель к компонентам начального поля напряжений согласно зависимости Булычева, показали завышенные показатели по скорости реализации этапа первичной ползучести. При этом скорости интервала установившейся ползучести имели удовлетворительное соответствие, которое легко объясняется качественным описанием вязкого элемента через функцию гиперболического синуса, что позволило достоверно описать широкий диапазон интенсивностей напряжений. Для исключения значительного влияния эмпирического метода учета забоя были решены задачи с проходкой выработки во времени в пространственной постановке (рис.5, линия «с учетом влияния забоя»). Результаты показали, что подобранные параметры модели для разных групп позволяют достоверно рассчитать величину первичных деформаций ползучести, что очень важно при проектировании конструкций обделок и крепи. Данный показатель позволяет использовать в проектах расчетно-обоснованные величины минимального и максимального времени возведения постоянного крепления после проходки выработки.

Для решения проблемы чрезмерно быстрой реализации первичных деформаций ползучести при модельных расчетах была выполнена оптимизация нелинейного вязкоупругого элемента путем итерационного уменьшения одного из коэффициентов. Диапазон уменьшения для разных групп – от 3 до 5 раз (рис.5, линия «после оптимизации»).

Заключение. Верификация рассмотренной модифицированной модели Бюргерса показала, что она является более эффективной, чем модели Мора – Кулона или Хоека – Брауна в комбинации с эмпирическими степенными законами ползучести типа Нортона. Преимущества заключаются в следующем:

• модель корректна для более широкого диапазона напряжений, что актуально для больших глубин;

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

Ограничения, требующие учет при анализе решений:

• пластический элемент не зависит непосредственно от скорости деформирования, в отличие от вязко-пластических моделей типа Деринга – Киля [17, 27], получивших широкое применение на Верхнекамском месторождении калийно-магниевых солей [1];

• модифицированная модель не учитывает влияние влажности на ползучесть, в отличие от модели Гампеля – Шульца, которая используется для соляных пород [14] и учитывает соответствующую зависимость Гюнша – Шульца;

• не учитывает восстановления нарушенных горных пород, например, вследствие перекристаллизации [15, 22];

• модель может быть расширена введением экспоненциальной зависимости Аррениуса от температуры в функции вязкости при необходимости учета сопряженного температурного воздействия [24];

• по сравнению с оригинальной вязко-упруго-пластической моделью Минкли [6] и ее обновленными формулировками [10, 16], используемая функция поверхности разрушения является линейной и не имеет границы максимальных сдвиговых напряжений σmax. Данное ограничение может быть критично в случае фиксации низких величин σmax (незначительный рост прочности при увеличении обжатия) при испытаниях некоторых соляных пород с низкой прочностью, таких как карналлит.

Предложенная модель, в сравнении с другими распространенными моделями реологического поведения соляных пород, соответствует модели Минкли [8, 9], сохраняя ее основные преимущества и недостатки. Следует подчеркнуть, что некоторые проблемы, касающиеся учета эффекта дилатансии в модели, могут быть исправлены, например, путем введения дополнительных уточненных нелинейных зависимостей угла дилатансии от напряжений в используемую зависимость от эффективных пластических деформаций [18].

В статье приведен процесс определения, оптимизации и верификации параметров предложенной расчетной модели поведения соляных пород с учетом разупрочнения и ползучести. Показано, что предложенную модель корректно использовать для широкого спектра геотехнических задач. Она комбинирует модифицированные модели: Бюргерса с зависимостью коэффициентов вязкости от напряжений через функцию гиперболического синуса и Мора – Кулона с нелинейными зависимостями коэффициента сцепления, углами внутреннего трения и дилатансии от эффективных пластических деформаций.

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

Литература