Размещено на http://www.allbest.ru/

Размещено на http://www.allbest.ru/

Исследование процесса деформирования массива каменной соли, содержащего подземное нефтегазохранилище

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

Приведем результаты испытания прочностных и деформативных свойств образцов солей Солигорского месторождения при одноосном сжатии [1]. Ползучесть призматических образцов размером см исследовалась на гидравлических и пружинных прессах при циклически возрастающей и длительно действующей нагрузках. Установлено, что мгновенное (t=0) деформирование соляных образцов линейно, а при t>0 развивается процесс нелинейной ползучести каменной соли с преобладанием необратимых деформаций.

Подобные процессы могут быть описаны уравнением состояния вида:

, (1)

причем уровень нелинейности деформирования зависит от степени нагруженности образца: для нагрузок, не превышающих 0,4 ( - напряжение разрушения при одноосном сжатии), получены практические линейные изохронны, т.е. в этом случае можно воспользоваться уравнениями линейной теории вязкоупругости [2], успешно применяемыми для описания ползучести горных пород. Если нагрузки не удовлетворяют упомянутому ограничению, то нелинейный член в уравнении (1) необходимо сохранить. Более того, при значительных нагрузках в уравнении (1) можно опустить второе слагаемое без ущерба точности кривых ползучести.

Лабораторные испытания трубчатых образцов различных соляных пород в условиях плоской деформации [3] позволили обобщить уравнение (1) на случай сложного напряженного состояния и для больших нагрузок представить уравнение механического состояния солей в виде уравнения, представляющего собой уравнения нелинейной вязкоупругости

, (2)

ядро которых

, (3)

или уравнения теории старения, если

, (4)

где - интенсивность напряжений, D - параметры ползучести.

Использование теории старения наиболее предпочтительно. В этом случае второй член в уравнениях (2) дает необратимую деформацию ползучести, что вполне согласуется с экспериментами [3]. В результате обработки лабораторных данных и натурных наблюдений за конвергенцией горизонтальных протяженных выработок круглого поперечного сечения определены параметры ползучести ядер (3), (4).

Согласно данным [3], соотношения линейной вязкоупругости (когда Р(t)=0) применимы при расчете мало заглубленных подземных сооружений, возводимых в соляных толщах. В этом случае задача упрощается в силу независимости поля напряжений от времени.

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

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

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

, (5)

где [D] матрица упругих констант.

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

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

. (6)

нелинейный напряжение соль

В этом случае путем подбора начальных напряжений, сводимых к вектору начальных узловых сил, искомые значения напряжений определяются последовательным решением ряда задач линейной теории упругости. Итерационный процесс осуществляется следующим образом. Расчетный временной интервал вновь разбивается на необходимое число точек. При t=0 решается краевая задача линейной теории упругости. По распределению упругих деформаций и согласно выражению (6) строятся начальное поле напряжений и соответствующий ему суммарный вектор начальной узловой нагрузки , затем устанавливаются поправки в результат первого упругого расчета. Для следующего временного шага процедура повторяется вновь, но с откорректированным полем напряжений. Итерационный процесс прекращается, как только в двух смежных временных шагах напряжения станут достаточно близкими.

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

Обратимся к вопросу сходимости метода упругих решений в задачах нелинейной теории ползучести. Математическое обоснование уравнений механического состояния нелинейной теории вязкоупругости предложено в работе [5]. Вводя гильбертовы пространства , функций Е(), , под которыми понимаются соответственно тензора деформаций и напряжений , определенные на отрезке времени [0, t], и задавая нормы и в виде

,

где , положительные функции памяти, авторы показывают, что любые аналитические в окрестности нуля операторы F, Q, отображающие соответственно на и на =F(E), E=Q(), можно представить в виде рядов:

, (7)

где

,

.

Точность аппроксимации операторов F, Q зависит от числа членов, сохраняемых в рядах (7) при замене бесконечных сумм конечными. Авторами получены условия, при которых операторы F и Q взаимно обратны, причем обращение осуществляется методом сжатых отображений.

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

(8)

и введем оператор Р

, (9)

где 0, , интенсивность деформаций, а ядро

.

Предположим, что , элементы пространства С [0,] с нормой . Известно [6], что операторы типа (9), определенные на С [0,], дифференцируемы в этом пространстве по Фреше. В частности, производная по Фреше от Р определяется соотношением:

,

,

при этом суммирование по i не производится.

Выделим в С [0,] шар . Если в R имеет место неравенство , то для любых

. (10)

В силу условия (10) при L<1 оператор Р является сжимающим в R, т.е. его неподвижную точку можно найти методом последовательных приближений, принимая за начальное нулевое приближение, удовлетворяющее неравенству [5]

. (11)

Определим условия, при которых L<1:

.

Таким образом, , если

. (12)

Для того чтобы удовлетворить (11), примем

, (13)

предполагая, что условия (12), (13) выполнены, и используя (9), в результате второго приближения получим уравнение:

, (14)

где интенсивность напряжений.

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

В работе [7] исследована сходимость метода упругих решений для уравнений механического состояния вида

, (15)

где ядра разбиваются на сингулярную и регулярную составляющие:

.

При этом некоторая вектор-функция u называется обобщенным решением краевой задачи для области V с границей S:

, (16)

, (17)

если u удовлетворяет интегральному тождеству:

,

где .

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

(18)

и ищется обобщенное решение в пространстве Н, которое получается замыканием по норме (18) множества дважды непрерывно дифференцируемых векторов-функций, удовлетворяющих (16). При этом метод упругих решений краевой задачи (16), (19) для уравнений механического состояния (15) сходится к единственному решению, если уравнение однозначно разрешимо в виде , а функция таковы, что для любого :

, (19)

,

причем , а начальное приближение выбрано из условия , где М константа, ограничивающая .

Уравнения (8) вытекают из (15), если в последних положить и во втором интеграле вместо ввести функцию .

Приведенная теорема о сходимости метода упругих решений распространяется на (8), а следовательно, и на уравнения (14). Если , то для удовлетворения неравенства (19) достаточно положить . Тогда метод упругих решений для уравнений механического состояния (8), (14) сходится, если и начальное приближение удовлетворяют неравенствам:

(20)

В промежутке времени [ ч] компоненты девиаторов напряжений и деформаций, соответствующие параметрам упругости и ползучести каменной соли, обеспечивают выполнение условий (12), (13) и (20). За начальное приближение принимается упругое решение.

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

1. Ержанов, Ж.С. Об оценке устойчивости формы осесимметричной полости в соляном массиве / Ж.С. Ержанов, Г.А. Аршинов, Э.И. Бергман // Известия АН КазССР. Сер. Физ.-мат. - 1974. - №5.

2. Ержанов, Ж.С. Теория ползучести горных пород и ее приложения / Ж.С. Ержанов. - Алма-Ата: Наука, 1964.

3. Ержанов, Ж.С. Ползучесть соляных пород / Ж.С. Ержанов, Э.И. Бергман. - Алма-Ата: Наука, 1977.

4. Зенкевич, О.С. Метод конечных элементов в технике / О.С. Зенкевич. - М.: Мир, 1975.

5. Илюшин, А.А. Основы математической теории термовязкоупругости / А.А. Илюшин, Б.Е. Победря. - М.: Наука, 1970.

6. Люстерник, Л.А. Элементы функционального анализа / Л.А. Люстерник, В.И. Соболев. - М.: Наука, 1965.

7. Победря, Б.Е. О сходимости метода упругих решений в нелинейной вязкоупругости / Б.Е. Победря // ДАН АН СССР. - 1970. - Т. 195. - №2.

Размещено на Allbest.ru