CC BY
80
ВЕСТНИК ПНИПУ
2014 Механика № 4
DOI: 10.15593/perm.mech/2014.4.03 УДК 51-72:539.4.013:539.411.4
А.В. Димаки1,2, Е.В. Шилько1,2, С.В. Астафуров1,2, С.Ю. Коростелев1, С.Г. Псахье1,2
1Институт физики прочности и материаловедения СО РАН, Томск, Россия Томский государственный университет, Томск, Россия
СВЯЗАННАЯ МОДЕЛЬ ФЛЮИДОНАСЫЩЕННЫХ ПОРИСТЫХ МАТЕРИАЛОВ, ОСНОВАННАЯ НА СОВМЕЩЕНИИ ДИСКРЕТНОГО И КОНТИНУАЛЬНОГО ПОДХОДОВ
Предложена численная модель флюидонасыщенных трещиновато-пористых хрупких материалов, основанная на гибридном подходе, совмещающем метод частиц и метод сеток. В рамках данной модели вмещающее пористое тело описывается на основе метода дискретных элементов. На ансамбле дискретных элементов моделируются процессы деформации пористого твердого тела и фильтрации однофазного флюида в связанной сети «микропор» (пор, каналов и других несплошностей, содержащихся в объеме дискретных элементов). Предложены соотношения, связывающие напряжения и деформации дискретного элемента, изменение объема его порового пространства и поровое давление флюида в «микропорах». Вычисление массоперено-са жидкости между «микропорами» и «макропорами» (которые рассматриваются как области между пространственно разнесенными и невзаимодействующими дискретными элементами) осуществляется на более мелкой сетке, «вмороженной» в лабораторную систему координат. С использованием развитой связанной модели проведено изучение механического отклика образцов упруго-хрупкого материала, поровое пространство которых насыщено водой, при одноосном сжатии. Показано, что эффективная прочность таких образцов определяется не только прочностными свойствами «сухого» (ненаполненного жидкостью) материала и давлением жидкости в порах, но в значительной степени - соотношением геометрических размеров образцов, скоростью деформирования и характеристиками пористости материала. По результатам моделирования предложена обобщающая зависимость прочности проницаемого водонасыщенного хрупкого материала при одноосном сжатии от приведенного диаметра фильтрационных каналов, представляющего собой отношение характерного диаметра фильтрационных каналов к квадратному корню скорости деформирования образца. Приведенные результаты демонстрируют возможности построенной модели для изучения нестационарных процессов деформирования и разрушения флюидонасыщенных материалов в условиях динамического нагружения.
Ключевые слова: трещиновато-пористая среда, жидкость, численное моделирование, метод дискретных элементов, метод конечных разностей, поровое давление, упруго-хрупкий материал, разрушение, прочность.
A.V. Dimaki12, E.V. Shilko12, S.V. Astafurov12, S.Yu. Korostelev1, S.G. Psakhie12
Institute of Strength Physics and Materials Science SB RAS, Tomsk, Russian Federation Tomsk State University, Tomsk, Russian Federation
COUPLED MODEL OF FLUID-SATURATED POROUS MATERIALS BASED ON A COMBINATION OF DISCRETE AND CONTINUUM APPROACHES
The numerical model of fluid-saturated porous brittle materials is proposed. The model is based on a hybrid approach combining the particle-based numerical method with finite difference method. In the framework of the model an enclosing porous solid is described with the discrete element method. An ensemble of discrete elements is used to model the processes of deformation of a porous solid and filtration of single-phase fluid in interconnected network of "micropores" (which are the pores, channels and other discontinuities enclosed in the volume of discrete elements). Relations between stress and strain of a discrete element, the change in volume of its pore space and fluid pore pressure in the "micropores" are proposed. Fluid mass transfer between the "micropores" and "macropores" (which are considered as the areas between spatially separated and non-interacting discrete elements) is calculated on the finer grid freezed in the laboratory system of coordinates. The developed coupled model was applied to study the mechanical response of samples of elastic-brittle material with water-saturated pore space under uniaxial compression. It is shown that the strength of fluid-saturated samples is determined by not only the strength properties of "dry" (unfilled) material and the fluid pore pressure, but largely by sample aspect ratio, rate of deformation and the characteristics of porosity of the material. Analysis of simulation results allowed authors to suggest a generalizing dependence of the uniaxial compressive strength of water-saturated permeable brittle material on the reduced diameter of filtration channels, which is the ratio of the characteristic diameter of the filtration channels to the square root of the sample strain rate. Presented results demonstrate the ability of the developed model to study non-stationary processes of deformation and fracture of fluid-saturated materials under dynamic loading.
Keywords: fractured porous medium, fluid, numerical modeling, discrete element method, finite difference method, pore pressure, elastic-brittle material, fracture, strength.
Введение
Многие природные и техногенные материалы и среды, такие как проницаемые горные породы (в том числе угольные и нефтеносные пласты) [1, 2], костная ткань [3], фильтрующие материалы, эндопроте-зы и т.д., представляют собой пористые или трещиновато-пористые флюидонасыщенные среды. Механический отклик проницаемых флюидонасыщенных материалов характеризуется рядом особенностей, которые отличают его от твердофазных композиционных материалов, в которых фазы также имеют весьма различные упругие и реологические характеристики. Эти особенности связаны со способностью жидкой фазы перераспределяться в трещинно-поровом пространстве вме-
щающего материала (макроскопический модуль сдвига жидкости традиционно полагается равным нулю). В результате такого перераспределения могут иметь место выравнивание среднего напряжения в различных областях материала, а также его значительные колебания, связанные с дилатансией и заполнением новых несплошностей, а также с обменом жидкостью с окружающей средой (в случае проницаемых границ).
Общим теоретическим подходом к изучению проницаемых флюидонасыщенных материалов является использование связанных моделей, учитывающих три важнейших аспекта функционирования таких систем: 1) взаимосвязь деформации твердофазного скелета с изменением объема и просветности трещинно-порового пространства; 2) взаимосвязь порового давления и напряжений во вмещающем объеме твердого тела; 3) процессы перераспределения флюида в трещинно-поровом пространстве. Наиболее известными представителями этого подхода являются аналитические макроскопические модели пороупру-гости, теоретические основы которых были заложены в работах Био [4-6]. Дальнейшее их развитие связано с учетом многомасштабности реальных материалов, процессов накопления повреждений, дилатансии и их влияния на упругие свойства каркаса и поровое давление флюида [7-10]. Возможности указанных и аналогичных макроскопических связанных моделей, как правило, ограничены проведением оценок напряженно-деформированного состояния проницаемых материалов (в первую очередь геологических) в стационарных условиях. При этом «за кадром» остается важный класс проблем, связанных с нестационарными и переходными процессами при динамическом нагружении, динамическом изменении граничных условий или разрушении. Основным инструментом теоретического решения подобных задач является численное моделирование с привлечением различных связанных моделей на основе концепции пороупругости.
В настоящее время численное моделирование флюидонасыщен-ных проницаемых материалов осуществляется преимущественно с использованием методов механики сплошных сред (методов конечных элементов и конечных разностей) [11-22]. Это связано в первую очередь с возможностью реализации в рамках этих методов различных реологических (в том числе связанных и многоуровневых) моделей твердых тел и изучения механического отклика материалов на различ-
ных масштабных уровнях. Несмотря на очевидные и известные достоинства подхода сплошной среды, в его рамках затруднено изучение процессов во флюидонасыщенных материалах и средах, сопровождающихся множественным разрушением твердого каркаса и контактным взаимодействием образующихся поверхностей.
Указанные ограничения при моделировании разрушения не присущи вычислительным методам, относящимся к дискретной концепции представления среды [23-28]. Известным представителем данного класса методов является метод дискретных элементов (МДЭ), в рамках которого моделируемый материал представляется ансамблем взаимодействующих частиц конечного размера [29-34]. Основные его преимущества связаны со способностью дискретных элементов изменять окружение, что является принципиально важным при моделировании сложных явлений контактного взаимодействия, трещинообразования и разделения тел на фрагменты и т.д.
Несмотря на это, область применения традиционных реализаций МДЭ для численного решения задач механики флюидонасыщенных материалов, как правило, ограничивается микроскопическим масштабным уровнем, на котором поры, трещины и каналы учитываются явным образом. При этом широко используются комбинированные схемы, в которых МДЭ применяется для моделирования механического отклика твердофазного каркаса, а метод решеточных уравнений Больцмана, или метод сеток, - для моделирования течения слабосжи-маемых жидкостей [35-39]. В то же время возможности использования МДЭ для моделирования отклика проницаемых материалов на более высоких масштабных уровнях являются весьма ограниченными ввиду недостаточного развития математического формализма данного метода для адекватного описания многомасштабных систем. В частности, одним из широко распространенных недостатков существующих реализаций МДЭ является применение парного приближения при описании взаимодействия дискретных элементов, что затрудняет моделирование сред со сложной реологией, в том числе упруго-пластических.
Авторами предложен подход к описанию механического отклика многомасштабных флюидонасыщенных упруго-пластических сред в рамках МДЭ, основанный на использовании многочастичного приближения при построении соотношений для сил межэлементного взаимодействия [40, 41].
Ранее на примере изотропных материалов было показано, что использование предложенных соотношений делает возможным реализацию в рамках МДЭ различных моделей и критериев упругости, пластичности и прочности (в том числе макроскопических моделей пластичности горных пород). Дальнейшее развитие предложенного формализма позволило расширить область применения МДЭ на широкий спектр пространственных и структурных масштабов, а также предложить комбинированный численный метод для изучения проницаемых сред на мезо- и макроскопическом масштабных уровнях. Данный метод получил название гибридного метода клеточных автоматов (ГКА). В методе ГКА механический отклик вмещающего твердого тела описывается методом дискретных элементов, а процессы фильтрации и диффузии флюида в трещинно-поровом пространстве (учитываемом неявно) - методом конечных разностей. Отметим, что ранее в рамках метода ГКА была развита связанная модель газоносных материалов [42-44].
В настоящей работе предложена и реализована в рамках метода ГКА связанная макроскопическая модель трещиновато-пористого хрупкого материала, насыщенного жидкостью. Данная модель учитывает взаимовлияние процессов деформирования твердофазного скелета и перераспределения жидкости в поровом пространстве. Возможности метода продемонстрированы на примере параметрического исследования зависимостей прочности упруго-хрупких флюидонасыщенных материалов от размеров пор, геометрии образца и скорости деформирования.
1. Формализм метода гибридных клеточных автоматов
Основу метода ГКА составляет декомпозиция решаемой задачи на две подзадачи: 1) описание механического поведения вмещающего твердого тела (каркаса) и 2) описание переноса жидкости в фильтрационном объеме, представленном системой взаимосвязанных каналов, пор, трещин и т.д. Для описания механического поведения вмещающего тела используется метод подвижных клеточных автоматов (ПКА), развиваемый авторами и являющийся одной из реализаций метода дискретных элементов [40, 41]. В рамках метода ПКА размеры и форма дискретных элементов могут изменяться в процессе деформирования. В зависимости от особенностей структуры рассматриваемого про-
ницаемого материала и моделируемого масштабного уровня размеры дискретных элементов могут значительно превышать линейные размеры несплошностей, быть сопоставимыми с ними или меньшими. Влияние «микропор» твердого каркаса (т.е. пор, каналов и других несплош-ностей, характерный размер которых меньше размера частицы) на механические свойства и отклик дискретного элемента учитывается неявно. Также в рамках метода ПКА решается задача фильтрационного переноса жидкости в связанной сети «микропор» вмещающего твердого тела. В то же время вычисление массопереноса жидкости между «микропорами» и «макропорами» (которые рассматриваются как области между пространственно разнесенными и невзаимодействующими частицами) осуществляется на более мелкой сетке, «вмороженной» в лабораторную систему координат (рис. 1, а). Эта же сетка используется и для расчета объемов «макропор».
Модель механического отклика флюидонасыщенного твердого каркаса
Метод ПКА, использующийся для моделирования деформирования и разрушения твердофазного каркаса, относится к группе методов однородно деформируемых дискретных элементов (в английской терминологии - simply deformable discrete element method) [30, 45-47], в рамках которых используется приближение однородного распределения напряжений и деформаций в объеме элемента. Для моделирования отклика трещиновато-пористых хрупких материалов в рамках метода ПКА реализована модель пластичности горных пород с неассо-циированным законом течения и критерием Мизеса-Шлейхера достижения предела упругости (модель Николаевского) [48, 49]. Выбор данной модели связан с тем, что она адекватно описывает отклик широкого класса хрупких материалов (геоматериалов, керамик) на различных масштабных уровнях с учетом вкладов «нижележащих» структурных масштабов. Особенностью модели Николаевского является постулируемая линейная связь между скоростями объемной и сдвиговой составляющих пластической деформации с коэффициентом пропорциональности Л, называемым коэффициентом дилатансии.
Адаптация модели Николаевского к методу ПКА осуществлена с использованием алгоритма Уилкинса [50]. В рамках этого алгоритма решение упруго-пластической задачи сводится к решению упругой за-
дачи в приращениях и последующей корректировке потенциальных сил взаимодействия частиц с соблюдением необходимых требований модели Николаевского, предъявляемых к величинам локального давления и девиатора напряжений [49].
Конечно-разностная сетка
а б
Рис. 1. Слои, образованные дискретными элементами (подвижными клеточными автоматами) и конечно-разностной сеткой (а); ячейки сетки на границе твердого
каркаса и макропоры (б)
Решение упругой задачи на ансамбле частиц сводится к определению приращений центральной и тангенциальной компонент силы упругого отклика подвижного клеточного автомата г на воздействие со стороны соседнего автомата у. Соответствующие выражения записываются на основе обобщенного закона Гука в гипоупругой форме [41]:
Ас.(.) = АКсе)7£, = 20, Ае.(.) +
>(.) >0) / . > '0)
(
1-
20,
Л
К>
Да"
Ат
{.)
:АК (7/ Бу = 20 Ах
> У
(1)
где символ А означает приращение величины соответствующего параметра за шаг по времени Аt численной схемы интегрирования; с,, и т,, -удельные значения центральной (КУ) и тангенциальной (Ку ) составляющих потенциальной силы реакции дискретного элемента г на воздействие со стороны соседа у; Б. - площадь поверхности контакта пары; 0 и К - модули сдвига и всестороннего сжатия материала элемента >'; Ае>( j) и Ау>(.) - приращения нормальной и сдвиговой деформации элемента г в паре г—с"еап - среднее напряжение в объеме эле-
• /-ч —mean
мента г. Среднее напряжение аг вычисляется с использованием соотношения Лява, связывающего компоненты тензора усредненных напряжений (огаР) в объеме дискретного элемента i с силами, действующими на поверхность данного объема [40, 41].
Напряженное состояние проницаемого твердого тела, содержащего связанную систему пор, каналов и трещин, является сложным и обусловлено не только величиной удельного объема несплошностей, но и особенностями их геометрии и пространственного распределения [51]. При отсутствии выраженной геометрии и ориентации (текстуры) трещинно-порового пространства адекватным является приближение, в соответствии с которым на мезо- и макроскопическом масштабных уровнях учитывается только вклад давления внутрипорового флюида в величину гидростатического напряжения в скелете (гидростатическое растяжение). В этом приближении влияние флюида, содержащегося в «микропорах», на механический отклик дискретного элемента может быть учтено посредством следующей модификации соотношения (1) для силы центрального взаимодействия:
( Дpfluid ^
Кг J
С 2G ^
Да,. = 2G Дг.,. - + 1 Да", (2)
1{]) '(1) / V г t(]) Y Y
К К J
где Р/1""1 - вклад порового давления флюида в «микропорах» в величину среднего напряжения в объеме дискретного элемента г. Отметим, что выражение (2) аналогично записи закона Гука в линейной модели пороупругости [6].
Величина Р/1""1 линейно связана со средним поровым давлением
флюида рроге в «микропорах» дискретного элемента г.
РАШ = а11,1роге, (3)
где аг - безразмерный коэффициент (0< аг <1), величина которого связана с характеристиками «микропор» (пористостью, топологией поро-вого пространства, размерами пор и т. д.).
После решения упругой задачи для элемента г на текущем временном шаге проверяется условие достижения предельного состояния (критерий Мизеса-Шлейхера), учитывающее вклад порового давления флюида.
ф ( + ЪрГ )Д/3 = у, (4)
где Уг - предел упругости материала элемента 1 при сдвиге; (3, - коэффициент внутреннего трения; о®9 - интенсивность усредненных напряжений в объеме элемента г; Ъг - безразмерный коэффициент. Ввиду конечной величины шага интегрирования уравнений движения Дt фактически проверяется условие Фг > Уг. При его выполнении осуществляется приведение тензора усредненных напряжений огаР в объеме
дискретного элемента г к предельной поверхности. Как показано в [41], такое приведение реализуется путем корректировки удельных значений сил центрального и тангенциального взаимодействия элемента г с соседями у.
I _ -=.твап\ \ 4 . (-=твап лт \
0'(У)=()-0>- )) + ( - ),
т'(.)=т.1.М
,(.) '0) 1
(5)
где ((.),Т',.)) - приведенные значения удельных сил реакции автомата; Мг = 1 -((/ар )( (Ф, -у )/(К,АД + 3Ог)) - коэффициент приведения девиаторов напряжений; N. = КЛ (Ф, - У )/(КгЛД. + ^) -
поправка к величине среднего напряжения, вычисленного в упругом приближении; Л, - коэффициент дилатансии материала автомата '.
Под влиянием напряжений в твердом каркасе изменяется объем вмещающего твердого тела и, соответственно, объем связанной системы несплошностей (именуемой также открытой пористостью или фильтрационным объемом). Удельный объем связанной системы «микропор» характеризуется величиной «микропористости» ф. В рамках предложенной модели связь между открытой микропористостью дискретного элемента и его объемной деформацией устанавливается следующим образом.
ф'=ф0 (1 + сг) + , (6)
где ф0 и ф - начальное (при нулевом давлении в каркасе) и текущее значения открытой микропористости; с, - безразмерный коэффициент
пропорциональности (с>>1); и О.рш - упругая и пластическая со-
ставляющие объемной деформации элемента, определяемые следующим образом:
Здесь егаа - диагональные компоненты тензора усредненных деформаций в объеме дискретного элемента г [40, 41].
Уравнения (2)-(7) составляют математическую формулировку модели взаимодействия твердого каркаса и жидкого флюида в «микро-поровом» пространстве, реализованной в рамках метода гибридных клеточных автоматов. Важным вопросом при использовании этой модели является определение значений безразмерных коэффициентов аг, Ь> и с> для изучаемого материала. Одним из вариантов такого определения для коэффициентов аг и с> является использование линейной модели пороупругости Био [6]. В этом случае коэффициент аг соответствует экспериментально определяемой константе а скорости изменения удельного объема жидкости при деформировании вмещающего твердого тела, а коэффициент с> - константе а для частично заполненных жидкостью пор либо произведению (В/Ки/К/г), где В - коэффициент Скемптона, Ки и К^ - модули всестороннего сжатия флюидонасыщен-ного материала и флюида. Следует отметить, что такое определение коэффициентов аг и с> может рассматриваться как оценка сверху. Оценкой снизу является приближение аг = фг и с> =1, которое становится адекватным в случае малых значений «микропористости» и «микро-пор» наноскопического масштаба [52]. Величина коэффициента Ьг традиционно полагается равной 1 [53, 54]. Однако необходимо отметить, что такое приближение справедливо для материалов, проницаемость и конфигурация порового пространства которых обеспечивают возможность равномерного распределения гидростатического давления в твердофазном каркасе (здесь также действует условие химически нейтральной жидкости) [55, 56]. Неявно здесь используется предположение о том, что образование новых повреждений (с которыми и связана макроскопическая квазипластичность хрупких материалов) привязано непосредственно к сети существующих микропор. В то же время
Qelast _з (с"еап + р—иЫ ) )К
г >
(7)
особенности порового пространства и микроскопической структуры твердофазного каркаса могут определять более сложную связь между поровым давлением и накоплением повреждений, в которой величина Ьг меньше 1 и является функцией параметров пористости и соотношения приложенного и порового давлений [52, 53, 56, 57-59]. В качестве нижней границы параметра Ьг может рассматриваться величина пористости ф> [53].
Развитый подход к построению многочастичного взаимодействия дает возможность использовать многопараметрические критерии разрушения (Кулона-Мора, Друккера-Прагера и т.д.) в качестве критерия разрыва химической связи (потери сцепления) в парах взаимодействующих дискретных элементов [41]. В представленных далее результатах численного моделирования применяется модифицированный критерий разрушения Друккера-Прагера, учитывающий вклад локального порового давления флюида аналогично (4):
с0р _ 0,5 (Я +1) а +1,5 (Я, -1) (а"еап + ЬРроге) _с<, (8)
где Я _ сс/оt - отношение прочностей сцепления пары элементов при
одноосном сжатии (ос) и растяжении (о^, коэффициент Ь соответствует тому же в соотношении (4). Детальное описание расчета критерия прочности для пар связанных дискретных элементов приведено в работах [40, 41]. Отметим, что, хотя критерии наступления неупругой деформации (4) и разрушения (8) имеют одинаковую структурную форму, смысл коэффициентов при слагаемом стеап в них различный. Так, величина коэффициента в в (4) определяет степень влияния давления на зарождение и развитие сети микроскопических повреждений и трещин в объеме материала. В то же время коэффициент 1,5(Я-1) определяет влияние давления на условие локализованного объединения отдельных микроповреждений в единый мезо- или макроразрыв. Сказанное свидетельствует о том, что значения этих двух коэффициентов не обязательно являются согласованными, поскольку роль локального давления в эволюции повреждения, как правило, существенно возрастает с ростом пространственного масштаба повреждения.
Уравнение состояния и модель фильтрационного переноса жидкости
В рамках развитой модели переноса жидкости используются следующие приближения: 1) жидкость может занимать объем поры как полностью, так и частично; 2) жидкость считается сжимаемой; 3) адсорбция жидкости на стенках пор, капиллярные эффекты и эффект адсорбционного понижения прочности проницаемого твердого тела не учитываются; 4) распределение «микропор» (то есть связанных между собой несплошностей в объеме дискретных элементов) по размерам не учитывается. Отметим, что в рамках последнего приближения трещин-но-поровое пространство дискретного элемента описывается двумя параметрами: величиной открытой «микропористости» ф и характерным диаметром фильтрационного канала ё^. Важно отметить, что величина ёс}, в общем случае равна не среднему поперечному размеру несплошностей, но размеру наименьших сужений в поровом пространстве, определяющих скорость фильтрации жидкости через проницаемый объем дискретного элемента. При корректном определении параметра ёс}, описанное упрощение при описании трещинно-порового пространства является адекватным для широкого класса проницаемых материалов, включая геоматериалы [60].
В развитой модели сжимаемая жидкость описывается следующим уравнением состояния [61]:
р(Р) _Р0 (1 + (р-Р0)/К1), (9)
где р и Р - текущие значения плотности и давления жидкости; р0 и Р0 -
значения плотности и давления жидкости при атмосферных условиях; К-1 - модуль всестороннего сжатия жидкости. Важной особенностью жидкости является ее способность частично занимать объем предоставленного порового пространства (в то время как газ занимает предоставленный объем целиком). Когда жидкость заполняет объем поры лишь частично, ее давление на твердый каркас может полагаться равным нулю (фактически давление определяется в этом случае весом жидкости, которым в рамках данной модели пренебрегаем). Соответственно, если отношение массы жидкости в поре к объему занимаемой поры меньше р0, ее давление полагается равным значению Р0.
Соотношение (9) определяет связь между плотностью и давлением жидкости как в «макропорах» (учитываемых явным образом как области между поверхностями несвязанных дискретных элементов), так и «микропорах» в объеме дискретных элементов. В последнем случае предполагается, что проницаемость и конфигурация «микропоро-вого» пространства дискретного элемента обеспечивают возможность равномерного распределения гидростатического давления в его твердофазном каркасе. Это позволяет оперировать для каждого дискретного элемента г усредненными по его объему значениями плотности р, и порового давления Ргроге содержащегося флюида (данное приближение отвечает формализму метода однородно деформируемого дискретного элемента).
При моделировании фильтрации используется предположение об определяющей роли градиента давления жидкости как движущей силы этого процесса. В пренебрежении гравитационными эффектами уравнение фильтрации (переноса плотности жидкости) в «микропоровом» пространстве твердого тела может быть записано следующим образом [61].
Ф^ = К Л
к У п
(10)
где п - вязкость жидкости; к - коэффициент проницаемости твердого каркаса, который может быть определен следующим образом [62].
к = Ф<Й. (11)
Уравнение (10) численно решается на сетке дискретных элементов методом Эйлера. Отметим, что с учетом используемых приближений массоперенос между узлами сетки, в которых р < р0, отсутствует.
Используемые в предложенной модели приближения подразумевают пренебрежение массопереносом жидкости между дискретными элементами за счет капиллярных сил. Для большинства горных пород это является справедливым в случае «микропор», характерный размер которых ёс}, не превышает первых микрометров. При больших значениях «капиллярный» вклад в фильтрацию может учитываться отдельным соотношением.
Моделирование массопереноса между микро-и макропорами
При описании твердого каркаса как ансамбля частиц возникает задача определения объемов макроскопических пор и полостей в твердом каркасе, а также вычисление протяженности границ раздела между твердым каркасом, содержащим микропоры, и макропорами. Как отмечалось выше, под макропорой (полостью) в данной модели понимается пора, размер которой больше размера подвижного клеточного автомата (дискретного элемента твердого каркаса). Для определения расположения и размеров макропор система дискретных элементов проецируется на конечно-разностную сетку (см. рис. 1, б). Далее, для расчета переноса жидкости на границе «твердофазный скелет-макропора» определяются узлы сетки, принадлежащие подвижным клеточным автоматам и граничащие с узлами, принадлежащими макропорам. Количество таких узлов для каждого подвижного автомата определяет протяженность границы его контакта с макропорой (либо с несколькими макропорами). Расчет переноса жидкости между макропорами и граничащими с ними подвижными клеточными автоматами, содержащими в своем объеме микропоры, осуществляется на основе уравнения (10).
Для описания перераспределения жидкости в объеме одной макропоры или между соединяющимися макропорами в работе используется приближение равного давления. В рамках данного приближения полагается, что в замкнутом макрообъеме на каждом шаге интегрирования плотность и давление жидкости распределены однородно. Такое упрощение является адекватным при низких скоростях процессов и в случае частично заполненных жидкостью макрообъемов. Отметим, что в более общем случае моделирование процессов массопереноса жидкости в системе соединяющихся макропор может осуществляться, например, с использованием метода решеточных уравнений Больцмана.
Верификация построенной модели проводилась на основе следующих тестов: 1) просачивание жидкости через тонкий слой и 2) выдавливание жидкости из образца при медленном нагружении (без разрушения). Результаты, получаемые в численных расчетах, отличались от соответствующих аналитических оценок для профилей давления, времени установления стационарного режима фильтрации и остаточного количества жидкости в образце не более чем на 1 %.
2. Изучение прочности пористых флюидонасыщенных образцов упруго-хрупкого материала при одноосном сжатии
Развитая модель была применена для моделирования одноосного сжатия пористых образцов, фильтрационный объем которых заполнен водой. Образцы были закреплены между неподвижной матрицей (снизу) и пуансоном (сверху), который двигался вниз, сжимая образец с постоянной скоростью Уу. Направление сжатия совпадало с вертикальной осью образца (рис. 2, а).
Г
Рис. 2. Схема нагружения моделируемого образца (а), диаграмма одноосного сжатия модельного упруго-хрупкого материала образца (б) (точкой на рис. (б) отмечено
значение ос)
Задача решалась в двумерной постановке в приближении плосконапряженного состояния. Структура образца полагалась однородной, макроскопические поры и включения отсутствовали. Свойства модельного материала образца соответствовали свойствам песчаника, диаграмма одноосного сжатия модельного материала приведена на рис. 2, б. В проведенных численных экспериментах были использованы следующие значения параметров модельного материала:
К = 37,5 ГПа, О = 5,77 ГПа, р = 2000 кг/м3, ос = 70 МПа, X = 7, а = ф = 0,1,
с = 1, К = 2,2 ГПа, значение параметра Ь варьировалось в пределах от Ь = ф = 0,1 до Ь = 1. Высота образца Н = 0,1 м, ширина образца Ж = 0,05 м.
о-к-—-—;—-—t—>—;—>—,—>—,
0,000 0,001 0,002 0,003 0,004 0,005 Деформация б
Матрица (Vy=0) а
Результаты моделирования выявили существенное влияние поро-вого давления жидкости на механический отклик хрупких пористых образцов, при этом имеющее более сложный характер, чем влияние давления газа в порах [44]. Помимо порового давления, на прочность образцов оказывают влияние такие факторы, как скорость нагружения, характерный диаметр фильтрационного канала (величина, совместно со значением открытой пористости характеризующая проницаемость материала) и геометрические размеры образца.
В условиях одноосного сжатия давление жидкости в порах твердого тела определяется конкуренцией двух процессов: 1) деформированием твердого каркаса, сопровождающимся уменьшением объема порового пространства в соответствии с (6)-(7), и ростом порового давления (и, следовательно, усилением влияния жидкости на величину напряжений в каркасе); 2) фильтрацией жидкости во внешнее пространство через боковые поверхности образца, приводящей к уменьшению порового давления и влияния жидкости на напряженное состояние каркаса.
Рис. 3. Зависимости прочности образцов модельного упруго-хрупкого материала при их одноосном сжатии от среднего диаметра фильтрационного канала при различных давлениях жидкости. Переменная Рш в легенде обозначает начальное давление жидкости в порах твердого каркаса. Скорость нагружения во всех случаях составляла
Уу = 0,1 м/с
Баланс этих двух процессов связан, в частности, с величиной коэффициента проницаемости материала, определяемого характерным диаметром фильтрационных каналов . Так, на рис. 3 приведены полученные путем компьютерного моделирования зависимости прочно-
сти образцов модельного материала от величины dh для различных
начальных значений порового давления жидкости Pinit. Можно видеть, что приведенные зависимости характеризуются наличием трех участков. На участке I (при достаточно низкой проницаемости образцов) отток жидкости из образца не компенсирует увеличение ее плотности в процессе деформирования пор. Вследствие этого по мере деформирования давление содержащейся в фильтрационном объеме жидкости постоянно возрастает, что приводит к снижению эффективной прочности образца (4). При больших значениях dch (участок III) обеспечивается скорость оттока жидкости из образца, достаточная для падения порового давления до атмосферного. Влияние жидкости на напряженное состояние таких образцов практически отсутствует, а их прочность стремится к прочности «сухого» образца. На переходном участке II скорость истечения жидкости из образца сопоставима со скоростью увеличения ее плотности.
Очевидно, что прочность образцов определяется, в том числе, распределением давления жидкости в поровом пространстве по сечению образца, которое, в свою очередь, зависит от отношения ширины образца W к его высоте H. Зависимость прочности образцов от dch при различных соотношениях W/H наглядно иллюстрируется результатами, приведенными на рис. 4, а. Как видно из данного рисунка, с увеличением ширины образца зависимость прочности от среднего диаметра фильтрационного канала становится более пологой. Иными словами, с увеличением соотношения W/H влияние внутрипоровой жидкости на величину прочности образца усиливается. Данный эффект объясняется следующим образом. С увеличением ширины образца доля материала, примыкающего к боковой поверхности и, соответственно, более подверженного оттоку жидкости во внешнее пространство, снижается, как видно из распределений порового давления в горизонтальном сечении образца (рис. 4, б). Соответственно, удельное количество и давление жидкости, остающееся в порах к началу разрушения, возрастает с увеличением соотношения W/H. Как следует из рис. 4, начиная с некоторого критического (достаточно большого) значения отношения W/H влияние проницаемости на прочность образцов становится несущественным, а значение прочности стремится к нижнему пределу (т. е. к прочности насыщенного жидкостью образца в отсутствие фильтрационного перераспределения жидкости при dch ^ 0 ).
Координата X,. а б
Рис. 4. Зависимости прочности образцов модельного упруго-хрупкого материала при их одноосном сжатии от среднего диаметра фильтрационного канала при различных значениях ширины W образца (а) и распределения порового давления в горизонтальном сечении образцов различной ширины к началу разрушения (б). Во всех случаях
Vy = 0,1 м/с, Pmlt = Po
Обобщенное выражение для прочности образца упруго-хрупкого флюидонасыщенного материала при одноосном сжатии
Как было отмечено выше, прочность флюидонасыщенных образцов при их одноосном сжатии определяется конкуренцией процессов механического деформирования образца под действием приложенной внешней нагрузки и оттока жидкости из порового пространства во внешнюю среду. Для выявления общих закономерностей такой конкуренции проведено детальное изучение влияния скорости нагружения на прочность образцов, насыщенных жидкостью. Моделировалось одноосное сжатие насыщенных жидкостью образцов, характеризующихся различными значениями dch, с различными скоростями нагружения
Vy. Начальное поровое давление жидкости Pimit полагалось равным атмосферному давлению (поровое пространство образцов полностью заполнено жидкостью). Рассматривались два приближения при учете влияния порового давления на прочность твердофазного каркаса.
1. Микропоры распределены однородно, размер их много меньше характерного размера формирующихся в процессе разрушения повреждений. В таком случае при относительно небольших значениях пористости можно полагать, что возникновение повреждений жестко не привязано к местам расположения микропор и в большей степени определяется наличием в твердом каркасе других дефектов большего (в сравнении с микропорами) масштаба. При этом влияние поровой
жидкости на прочность каркаса определяется величиной пористости и учитывается посредством следующего определения коэффициента Ь в (8): Ь = ф.
2. Микропоры распределены неоднородно, повреждения в материале формируются путем соединения трещиной нескольких микро-пор. В этом случае влияние поровой жидкости на прочность каркаса определяется непосредственно величиной порового давления: Ь = 1.
На рис. 5 приведены зависимости прочности образцов от характерного диаметра фильтрационного канала при различных скоростях деформации. Можно видеть, что значение параметра Ь, определяющего вклад порового давления в критерии пластичности (4) и разрушения (8), оказывает значительное влияние на прочность флюидонасыщен-ных образцов. Так, в случае использования приближения однородного распределения микроскопических пор (Ь = ф = 0,1, рис. 5, а) в твердом каркасе вклад порового давления оказывается достаточно слабым (максимальное снижение прочности образцов не превышает 3 %). В рамках второго приближения (Ь = 1, рис. 5, б) вклад порового давления в прочность образцов возрастает на порядок величины, а величина предельной деформации снижается. При этом имеет место «выполажи-вание» зависимостей прочности от диаметра фильтрационного канала, обусловленное более сильным влиянием остаточной жидкости в поро-вом пространстве на прочность образцов.
Рис. 5. Обобщенная зависимость прочности водонасыщенных образцов от диаметра фильтрационного канала при различных скоростях нагружения: а - приближение однородного распределения микроскопических пор (Ь = ф = 0,1); б - приближение относительно крупных микропор, развитие которых определяет формирование макротрещин (Ь = 1). Начальное давление жидкости в порах твердого каркаса Рш = Р0
Анализ полученных зависимостей прочности образцов от эффективного диаметра фильтрационных каналов при различных скоростях нагружения показал возможность сведения их в единую кривую зависимости прочности от приведенного диаметра фильтрационного канала
(12)
где е = ¥у / Н - скорость деформации; е0у - масштабный множитель,
имеющий размерность скорости деформации. Кривые, показанные на рис. 5, построены в указанных переменных.
Зависимости прочности от приведенного диаметра (см. рис. 5) могут быть аппроксимированы аналитическим выражением следующего вида:
1 + exP (( - d0)/
_min _ _0
°c (dch , tyy ) = °0---C - ° I-T~~—\ , (13)
...../ё yy/ e °y
где g0 - прочность образца при одноосном сжатии в отсутствие жидкости в поровом пространстве; от™ - прочность водонасыщенного образца в условиях «запирания» жидкости в порах (т.е. при отсутствии фильтрационного движения жидкости в твердом каркасе); d0 и dc -
параметры аппроксимирующей зависимости, имеющие размерность расстояния и характеризующие соответственно положение и ширину переходного участка зависимости прочности от характерного диаметра фильтрационных каналов. Величина omin может быть оценена на основе соотношений (6)-(8) следующим образом:
от1П
Г1 +1,5 (Х-l)bcf. (14)
K
J
Как следует из выражения (14), на прочность флюидонасыщен-ных образцов оказывают влияние не только прочностные свойства твердого каркаса, но и свойства жидкости, находящейся в поровом пространстве, а также характер пористости и удельный вклад порового давления в напряженное состояние каркаса. Для использованных в на-
стоящих расчетах параметров материала твердого каркаса и жидкости оценка о™" составляет 63,92 МПа для Ь = 0,1 и 44 МПа для Ь = 1, что хорошо согласуется с результатами моделирования для ёсЬ ^0
(см. рис. 5). Заметим, что оценка (14) дает немного заниженные значения прочности в связи с тем, что в ней не учитывается влияние порово-го давления жидкости на изменение объема пор.
В целом сложный характер взаимосвязей между параметрами, характеризующими механический отклик твердого каркаса, физико-механические свойства жидкости и динамику ее фильтрационного перераспределения в системе пор, определяет нелинейную зависимость прочности образцов от совокупности данных параметров и обусловливает необходимость применения численных методов для изучения отклика пористых флюидонасыщенных образцов и материалов.
Заключение
В работе рассмотрен подход к моделированию пористого материала, насыщенного жидкостью. В рамках данного подхода, названного методом гибридных клеточных автоматов, рассматриваемая среда моделируется как суперпозиция двух взаимосвязанных слоев, один из которых представлен ансамблем частиц (подвижных клеточных автоматов), а другой - конечно-разностной сеткой. Развитый подход позволяет явным образом учитывать как деформацию и разрушение твердого каркаса, так и перераспределение жидкости в поровом пространстве и в макроскопических полостях.
Развитый в работе формализм применен для параметрического исследования зависимостей прочности упруго-хрупких флюидонасы-щенных материалов от размеров пор, геометрии образца и скорости деформирования. Показано, что прочность упруго-хрупких образцов, насыщенных водой, при одноосном сжатии существенным образом зависит не только от свойств материала и давления жидкости в поровом пространстве, но и от размера пор, физико-механических свойств жидкости и геометрических размеров образца. Построена аналитическая зависимость прочности флюидонасыщенного образца при одноосном сжатии от указанных параметров, обобщающая результаты численного моделирования.
Приведенные результаты демонстрируют возможности построенной модели, а также показывают важность применения численного моделирования для изучения механических свойств (включая прочность) флюидонасыщенных материалов при динамическом нагружении.
Работа выполнена в рамках Программы фундаментальных научных исследований государственных академий наук на 2013-2020 годы.
Библиографический список
1. Zhejun P., Connell L.D. A theoretical model for gas adsorption-induced coal swelling // Int. J. of Coal Geol. - 2007. - Vol. 69. - No. 4. -Р.243-252.
2. Taylor D. Fracture and repair of bone: a multiscale problem // J. Mater. Sci. - 2007. - Vol. 42. - Р. 8911-8918.
3. Jing L., Stephansson O. Fundamentals of discrete element method for rock engineering: theory and applications. - London: Elsevier, 2007. -562 p.
4. Biot M.A. General theory of three-dimensional consolidation // J. of Appl. Phys. - 1941. - Vol. 12. - Р. 155-164.
5. Biot M.A. The elastic coefficients of the theory of consolidation // J. appl. Mech. - 1957. - Vol. 24. - Р. 594-601.
6. Detournay E., Cheng A.H.-D. Fundamentals of poroelasticity. Chapter 5 in Comprehensive Rock Engineering: Principles, Practice and Projects. Vol. II. Analysis and Design Method / ed. C. Fairhurst. - Perga-mon Press, 1993. - Р. 113-171.
7. Hamiel Y., Lyakhovsky V., Agnon A. Coupled evolution of damage and porosity in poroelastic media: theory and applications to deformation of porous rocks // Geophys. J. Int. - 2004. - Vol. 156. - Р. 701-713.
8. Lyakhovsky V., Hamiel Y. Damage Evolution and Fluid Flow in Poroelastic Rock // Izvestiya, Physics of the Solid Earth. - 2007. - Vol. 43. -No. 1. - Р. 13-23.
9. Мейрманов А.М. Метод двухмасштабной сходимости Нгует-сенга в задачах фильтрации и сейсмоакустики в упругих пористых средах // Сибирский математический журнал. - 2007. - Т. 48, № 3. -С. 645-667.
10. Бочаров О.Б., Рудяк В.Я., Серяков А.В. Простейшие модели деформирования пороупругой среды, насыщенной флюидами // Физи-
ко-технические проблемы разработки полезных ископаемых. - 2014. -№ 2. - С. 54-68.
11. Coupled fluid flow and geomechanical deformation modeling / S.E. Minkoff, C.M. Stone, S. Bryant, M. Peszynska, M.F. Wheeler // J. of Petroleum Sci. and Eng. - 2003. - Vol. 38. - Р. 37-56.
12. Micromechanically Based Poroelastic Modeling of Fluid Flow in Haversian Bone / C.C. Swan, R.S. Lakes, R.A. Brand, K.J. Stewart // J. of Biomech. Eng. - 2003. - Vol. 125. - Р. 25-37.
13. Silbernagel M.M. Modeling Coupled Fluid Flow and Geomechani-cal and Geophysical Phenomena Within a Finite Element Framework. -Golden, CO: Colorado School of Mines, 2007. - 218 р.
14. White J.A., Borja R.I. Stabilized low-order finite elements for coupled solid-deformation/fluid-diffusion and their application to fault zone transients // Comput. Methods Appl. Mech. Eng. - 2008. - Vol. 197. -Р. 4353-4366.
15. Jha B., Juanes R. Coupled multiphase flow and poromechanics: A computational model of pore pressure effects on fault slip and earthquake triggering // Water Resour. Res. - 2014. - Vol. 50. - Iss. 5. - Р. 3776-3808.
16. Turner D.Z., Nakshatrala K.B., Martinez M.J. Framework for Coupling Flow and Deformation of a Porous Solid // Int. J. Geomech. -2014. - Р. 04014076. (in press).
17. Masson Y.J., Pride Y.J., Nihei K.T. Finite difference modeling of Biot's poroelastic equations at seismic frequencies // J. Geophys. Res. -2006. - Vol. 111. - B10305.
18. Finite element modelling of the effective elastic properties of partially saturated rocks / D. Makarynska, B. Gurevich, R. Ciz, C.H. Arns, M.A. Knack-stedt // Computers & Geosciences. - 2008 - Vol. 34. - Р. 647-657.
19. Наседкина А. А., Наседкин А.В., Иоване Ж. Моделирование нестационарного воздействия на многослойный пороупругий пласт с нелинейными геомеханическими свойствами // Физико-технические проблемы разработки полезных ископаемых. - 2009. - № 4. - С. 23-32.
20. Наседкина А.А. Моделирование нестационарных процессов фильтрации в пороупругих средах с физическими нелинейностями // Вестник Нижегород. ун-та им. Н.И. Лобачевского. - 2011. - № 4-3. -С.1006-1008.
21. Доброскок А. А., Линьков А.М. Моделирование течения, напряженного состояния и сейсмических событий в породах при сбросе
давления в трещине гидроразрыва // Физико-технические проблемы разработки полезных ископаемых. - 2011. - № 1. - С. 12-22.
22. Численное моделирование динамики составного пороупруго-го тела / Л.А.Игумнов, С.Ю. Литвинчук, Д.В. Тарлаковский, Н.А. Локтева // Проблемы прочности и пластичности. - 2013. - № 75, Т. 2. -С. 130-136.
23. Rieth M. Nano-engineering in science and technology: An Introduction to the World of Nano-Design. - Singapore: World Scientific, 2003. - 164 p.
24. Cundall P.A., Strack O.D.L. A discrete numerical model for granular assemblies // Geotechnique. - 1979. - Vol. 29. - Р. 47-65.
25. Mustoe G.G.W. A generalized formulation of the discrete element method // Eng. Comp. - 1992. - Vol. 9. - Р. 181-190.
26. Shi G.H. Discontinuous deformation analysis - a new numerical model for statics and dynamics of block systems // Eng. Comp. - 1992. -Vol. 9. - Iss. 2. - Р. 157-168.
27. Munjiza A.A., Knight E.E., Rougier E. Computational mechanics of discontinua. - Chichester: Wiley, 2012. - 276 p.
28. Lisjak A., Grasseli G. A review of discrete modeling techniques for fracturing processes in discontinuous rock masses // Int. J. Rock. Mech. Min. Sci. - 2014 - Vol. 6. - Р. 301-314.
29. Munjiza A. The combined finite-discrete element method. -Chichester: Wiley, 2004. - 352 p.
30. Bicanic N. Discrete element methods // Encyclopedia of computational mechanics / eds. Stein E., Borst R., Hughes T.J.R. Vol. 1: Fundamentals. - Chichester: Wiley, 2004. - P. 311-371.
31. Jing L., Stephansson O. Fundamentals of discrete element method for rock engineering: theory and applications. - Elsevier, 2007. - 562 p.
32. Williams J.R., Hocking G., Mustoe G.G.W. The Theoretical basis of the discrete element method // Numerical methods of engineering, theory and applications / ed. A.A. Balkema. - Rotterdam: NUMETA, 1985.
33. Potyondy D.O., Cundall P.A. A bonded-particle model for rock // Int. J. Rock. Mech. Min. Sci. - 2004. - Vol. 41. - P. 1329-1364.
34. Хан Г.Н. О несимметричном режиме разрушения массива горных пород в окрестности полости // Физ. мезомех. - 2008. - Т. 11, № 1. - С. 109-114.
35. Zhao G.F., Khalili N. A Lattice Spring Model for Coupled Fluid Flow and Deformation Problems in Geomechanics // Rock Mech. and Rock Eng. - 2012. - Vol. 45. - P. 781-799.
36. Cook B.K., Noble D.R. A direct simulation method for particle-fluid systems // Eng. Comp. - 2011. - Vol. 21. - No. 2/3/4. - P. 151-168.
37. Sakaguchi H., Muhlhaus H.-B. Hybrid Modelling of Coupled Pore Fluid-solid Deformation Problems // Pure appl. geophys. - 2000. -Vol. 157. - P. 1889-1904.
38. Han Y., Cundall P.A. Lattice Boltzmann modeling of pore-scale fluid flow through idealized porous media // Int. J. Numer. Meth. Fluids. -2011. - Vol. 67. - P. 1720-1734.
39. Han Y., Cundall P A. LBM-DEM modeling of fluid-solid interaction in porous media // Int. J. Numer. Analyth. Meth. Geomech. - 2013. -Vol 37. - Iss. 10. - P. 1391-1407.
40. Development of a formalism of movable cellular automaton method for numerical modeling of fracture of heterogeneous elastic-plastic materials / S. Psakhie, E. Shilko, A. Smolin [et al.] // Fracture and Structural Integrity. - 2013. - Vol. 24. - P. 59-91.
41. A mathematical model of particle-particle interaction for discrete element based modeling of deformation and fracture of heterogeneous elastic-plastic materials / S.G. Psakhie, E.V. Shilko, A.S. Grigoriev, S.V. Asta-furov, A.V. Dimaki, A.Yu. Smolin // Engineering Fracture Mechanics, 2014 (in press).
42. Hybrid Cellular Automata Metod. Application to Research on Mechanical Response of Contrast Media / S. Zavsek, A.V. Dimaki, A.I. Dmitriev, E.V. Shilko, J. Pezdic, S.G. Psakhie // Phys. Mesomech. - 2013. - Vol. 1. -P. 42-51.
43. Approach to simulation of deformation and fracture of hierarchically organized heterogeneous media, including contrast media / S.G. Psakhie, E.V. Shilko, A.Yu. Smolin [et. al.] // Phys. Mesomech. - 2011. - Vol. 14. -No. 5-6. - P. 224-248.
44. Развитие формализма метода частиц для моделирования отклика флюидонасыщенных пористых геологических материалов / А.В. Димаки, Е.В. Шилько, С.В. Астафуров, С.Г. Псахье // Известия ТПУ. - 2014. - Т. 324, № 1. - С. 102-111.
45. Cundall P.A. Formulation of a three-dimensional distinct element model - Part I: A scheme to detect and represent contacts in a system com-
posed of many polyhedral blocks // Int. J. Rock Mech, Min. Sci. Geomech. Abstr. - 1988. - Vol. 25. - Iss. 3. - P. 107-116.
46. Hwang J.-Y., Ohnishi Y., Wu J. Numerical analysis of discontinuous rock masses using three-dimensional discontinuous deformation analysis (3D DDA) // Geotech. Eng. - 2004. - Vol. 8. - Iss. 5. - P. 491-496.
47. Hahn M., Wallmersperger T., Kroplin B.-H. Discrete element representation of discontinua: proof of concept and determination of material parameters // Comp. Mat. Sci. - 2010. - Vol. 50. - P. 391-402.
48. Garagash I.A., Nikolaevskiy V.N. Non-associated laws of plastic flow and localization of deformation // Adv. Mech. - 1989. - Vol. 12. -Iss. 1. - P. 131-183.
49. Stefanov Yu.P. Deformation localization and fracture in geomate-rials. Numerical simulation // Phys. Mesomech. - 2002. - Vol. 5-6. - P. 67-77.
50. Wilkins M.L. Computer simulation of dynamic phenomena. -Heidelberg: Springer-Verlag, 1999. - 246 p.
51. Kushch V.I., Shmegera S.V., Sevostianov I. SIF statistics in micro cracked solid: effect of crack density, orientation and clustering // Int. J. Eng. Sci. - 2009. - Vol. 47. - P. 192-208.
52. Permeability of Wilcox shale and its effective pressure law / O. Kwon, A.K. Kronenberg, A.F. Gangi, B. Johnson // J. Geophys. Res. -2001. - Vol. 106. - No. b9. - P. 19339-19353.
53. Paterson M.S., Wong T.F. Experimental Rock Deformation. The Brittle Field. - Berlin-Heidelberg: Springer-Verlag, 2005. - 347 p.
54. Yamaji A. An Introduction to Tectonophysics: Theoretical Aspects of Structural Geology. - Tokyo: TERRAPUB, 2007. - 378 p.
55. Experimental Deformation of Sedimentary Rocks Under Confining Pressure: Pore Pressure Tests / J. Hangin, R.V.Hager, M. Friedman, J.N. Feather // AAPG Bulletin. - 1963. - Vol. 47. - Iss. 5. - P. 717-755.
56. Ставрогин А.Н., Тарасов Б.Г. Экспериментальная физика и механика горных пород. - СПб.: Наука, 2001. - 343 с.
57. Robin P.-Y.F. Note on Effective Pressure // J. Geophys. Res. -1973. - Vol. 78. - No. 14. - P. 2434-2437.
58. Gangi A.F., Carlson R.L. An asperity-deformation model for effective pressure // Tectonophysics. - 1996. - Vol. 256. - P. 241-251.
59. Boitnott G.N., Scholz C.H. Direct Measurement of the Effective Pressure Law: Deformation of Joints Subject to Pore and Confining Pressure // J. Geophys. Res. - 1990. - Vol. 95. - No. B12. - P. 19279-I9298.
60. Хавкин А.Я. Наноявления и нанотехнологии в добыче нефти и газа. - М.; Ижевск: Регулярная и хаотическая динамика, 2010. - 692 с.
61. Басниев К.С., Кочина И.Н., Максимов В.М. Подземная гидромеханика. - М.: Недра, 1993. - 416 с.
62. Лейбензон Л. С. Движение природных жидкостей и газов в пористой среде. - М.; Л.: Гос. изд-во технико-теорет. лит., 1947. - 244 с.
References
1. Zhejun P., Connell L.D. A theoretical model for gas adsorption-induced coal swelling. Int. J. of Coal Geol., 2007, vol. 69, no. 4, pp. 243252.
2. Taylor D. Fracture and repair of bone: a multiscale problem. J. Mater. Sci, 2007, vol. 42, pp. 8911-8918.
3. Jing L., Stephansson O. Fundamentals of discrete element method for rock engineering: theory and applications. London: Elsevier, 2007. 562 p.
4. Biot M.A. General theory of three-dimensional consolidation. J. of Appl. Phys, 1941, vol. 12, pp. 155-164.
5. Biot M.A. The elastic coefficients of the theory of consolidation. J. appl. Mech, 1957, vol. 24, pp. 594-601.
6. Detournay E., Cheng A.H.-D. Fundamentals of poroelasticity. Chapter 5 in Comprehensive Rock Engineering: Principles, Practice and Projects, Vol. II, Analysis and Design Method, ed. C. Fairhurst. Pergamon Press, 1993. P.113-171.
7. Hamiel Y., Lyakhovsky V., Agnon A. Coupled evolution of damage and porosity in poroelastic media: theory and applications to deformation of porous rocks. Geophys. J. Int., 2004, vol. 156, pp. 701-713.
8. Lyakhovsky V., Hamiel Y. Damage Evolution and Fluid Flow in Poroelastic Rock. Izvestiya, Physics of the Solid Earth, 2007, vol. 43, no. 1, pp. 13-23.
9. Meirmanov A.M. Metod dvukhmasshtabnoi skhodimosti Nguet-senga v zadachakh fil'tratsii i seismoakustiki v uprugikh poristykh sredakh. [The Nguetseng method of two-scale convergence in the problems of filtration and seismoacoustics in elastic porous media]. Sibirskii matematicheskii zhurnal, 2007, vol. 48, no. 3, pp. 645-667.
10. Bocharov O.B., Rudiak V.Ia., Seriakov A.V. Prosteishie modeli de-formirovaniia porouprugoi sredy, nasyshchennoi fliuidami [Simple models of
deformation of poroelastic medium saturated with fluids]. Fiziko-tekhnicheskie problemy razrabotkipoleznykh iskopaemykh, 2014, no. 2, pp. 54-68.
11. Minkoff S.E., Stone C.M., Bryant S., Peszynska M., Wheeler M.F. Coupled fluid flow and geomechanical deformation modeling. J. of Petroleum Sci. and Eng., 2003, vol. 38, pp. 37-56.
12. Swan C.C., Lakes R.S., Brand R.A., Stewart K.J. Micromechani-cally Based Poroelastic Modeling of Fluid Flow in Haversian Bone. J. of Biomech. Eng., 2003, vol. 125, pp. 25-37.
13. Silbernagel M.M. Modeling Coupled Fluid Flow and Geomechanical and Geophysical Phenomena Within a Finite Element Framework. Golden, CO: Colorado School of Mines, 2007. 218 p.
14. White J.A., Borja R.I. Stabilized low-order finite elements for coupled solid-deformation/fluid-diffusion and their application to fault zone transients. Comput. Methods Appl. Mech. Eng., 2008, vol. 197, pp. 4353-4366.
15. Jha B., Juanes R. Coupled multiphase flow and poromechanics: A computational model of pore pressure effects on fault slip and earthquake triggering. Water Resour. Res., 2014, vol. 50, iss. 5, pp. 3776-3808.
16. Turner D.Z., Nakshatrala K.B., Martinez M.J. Framework for Coupling Flow and Deformation of a Porous Solid. Int. J. Geomech., 2014, p. 04014076. (in press).
17. Masson Y.J., Pride Y.J., Nihei K.T. Finite difference modeling of Biot's poroelastic equations at seismic frequencies. J. Geophys. Res, 2006, vol. 111, B10305.
18. Makarynska D., Gurevich B., Ciz R., Arns C.H., Knackstedt M.A. Finite element modelling of the effective elastic properties of partially saturated rocks. Computers & Geosciences, 2008, vol. 34, pp. 647-657.
19. Nasedkina A.A., Nasedkin A.V., Iovane G. Modeling and finite element analysis of the nonstationary action on a multi-layer poroelastic seam with nonlinear geomechanical properties. J. of Mining Sci, 2009, vol. 45, no. 4, pp. 324-333.
20. Nasedkina A.A. Modelirovanie nestatsionarnykh protsessov fil-tratsii v porouprugikh sredakh s fizicheskimi nelineinostiami [Simulation of non-stationary filtration processes in poroelastic media with physical non-linearity]. Vestnik Nizhegorodskogo universiteta im. N.I. Lobachevskogo, 2011, no. 4-3, pp. 1006-1008.
21. Dobroskok A.A., Linkov A.M. Modeling of fluid flow, stress state and seismicity induced in rock by an instant pressure drop in a hydrofracture. J. of MiningSci, 2011, vol. 47, no. 1, pp. 10-19.
22. Igumnov L.A., Litvinchuk S.Iu., Tarlakovskii D.V., Lokteva N.A. Chislennoe modelirovanie dinamiki sostavnogo porouprugogo tela [Numerical simulation of dynamics of a composite poroelastic body]. Problemy prochnosti iplastichnosti, 2013, vol. 75, no. 2, pp. 130-136.
23. Rieth M. Nano-engineering in science and technology: An Introduction to the World of Nano-Design. Singapore: World Scientific, 2003. 164 p.
24. Cundall P.A., Strack O.D.L. A discrete numerical model for granular assemblies. Geotechnique, 1979, vol. 29, pp. 47-65.
25. Mustoe G.G.W. A generalized formulation of the discrete element method. Eng. Comp, 1992, vol. 9, pp. 181-190.
26. Shi G.H. Discontinuous deformation analysis - a new numerical model for statics and dynamics of block systems. Eng. Comp, 1992, vol. 9, iss. 2, pp. 157-168.
27. Munjiza A.A., Knight E.E., Rougier E. Computational mechanics of discontinua. Chichester: Wiley, 2012. 276 p.
28. Lisjak A., Grasseli G. A review of discrete modeling techniques for fracturing processes in discontinuous rock masses. Int. J. Rock. Mech. Min. Sci, 2014, vol. 6, pp. 301-314.
29. Munjiza A. The combined finite-discrete element method. Chichester: Wiley, 2004. 352 p.
30. Bicanic N. Discrete element methods. Encyclopedia of computational mechanics. Eds. E. Stein, R. Borst, T.J.R. Hughes. Vol. 1: Fundamentals. Chichester: Wiley, 2004. P. 311-371.
31. Jing L, Stephansson O. Fundamentals of discrete element method for rock engineering: theory and applications. Elsevier, 2007. 562 p.
32. Williams J.R., Hocking G, Mustoe G.G.W. The Theoretical basis of the discrete element method. Ed. A.A. Balkema. Numerical methods of engineering, theory and applications. Rotterdam: NUMETA, 1985.
33. Potyondy D.O., Cundall P.A. A bonded-particle model for rock. Int. J. Rock. Mech. Min. Sci, 2004, vol. 41, pp. 1329-1364.
34. Khan G.N. O nesimmetrichnom rezhime razrusheniia massiva gornykh porod v okrestnosti polosti [Asymmetric failure of the rock mass near a hollow]. Fiz. Mesomekh., 2008, vol. 11, no. 1, pp. 109-114.
35. Zhao G.F., Khalili N. A Lattice Spring Model for Coupled Fluid Flow and Deformation Problems in Geomechanics. Rock Mech. and Rock Eng, 2012, vol. 45, pp. 781-799.
36. Cook B.K., Noble D.R. A direct simulation method for particle-fluid systems. Eng. Comp., 2011, vol. 21, no. 2/3/4, pp. 151-168.
37. Sakaguchi H., Muhlhaus H.-B. Hybrid Modelling of Coupled Pore Fluid-solid Deformation Problems. Pure appl. geophys., 2000, vol. 157, pp. 1889-1904.
38. Han Y., Cundall P.A. Lattice Boltzmann modeling of pore-scale fluid flow through idealized porous media. Int. J. Numer. Meth. Fluids, 2011, vol. 67, pp. 1720-1734.
39. Han Y., Cundall P.A. LBM-DEM modeling of fluid-solid interaction in porous media. Int. J. Numer. Analyth. Meth. Geomech., 2013, vol 37, iss. 10, pp. 1391-1407.
40. Psakhie S., Shilko E., Smolin A. [et al.] Development of a formalism of movable cellular automaton method for numerical modeling of fracture of heterogeneous elastic-plastic materials. Fracture and Structural Integrity, 2013, vol. 24, pp. 59-91.
41. Psakhie S.G., Shilko E.V., Grigoriev A.S., Astafurov S.V., Dima-ki A.V., Smolin A.Yu. A mathematical model of particle-particle interaction for discrete element based modeling of deformation and fracture of heterogeneous elastic-plastic materials. Engineering Fracture Mechanics, 2014 (in press).
42. Zavsek S., Dimaki A.V., Dmitriev A.I., Shilko E.V., Pezdic J., Psakhie S.G. Hybrid Cellular Automata Metod. Application to Research on Mechanical Response of Contrast Media. Phys. Mesomech, 2013, vol. 1, pp. 42-51.
43. Psakhie S.G., Shilko E.V., Smolin A.Yu. [et al.] Approach to simulation of deformation and fracture of hierarchically organized heterogeneous media, including contrast media. Phys. Mesomech, 2011, vol. 14, no. 5-6, pp. 224-248.
44. Dimaki A.V., Shilko E.V., Astafurov S.V., Psakhie S.G. Razvitie formalizma metoda chastits dlia modelirovaniia otklika fliuidonasyshchen-nykh poristykh geologicheskikh materialov [Development of particle method formalism for numerical modeling of the response of fluid_saturated porous geological materials]. Izvestiya Tomskogo politekhnicheskogo uni-versiteta, 2014, vol 324, no. 1, pp. 102-111.
A.B. fíuMaKU, E.B. ^rnbKO, C.B. Acma$ypoe, C.W. Kopocmeßee, C.F. ncaxbe
45. Cundall PA. Formulation of a three-dimensional distinct element model. Part I: A scheme to detect and represent contacts in a system composed of many polyhedral blocks. Int. J. Rock Mech, Min. Sci. Geomech. Abstr, 1988, vol. 25, iss. 3, pp. 107-116.
46. Hwang J.-Y., Ohnishi Y., Wu J. Numerical analysis of discontinuous rock masses using three-dimensional discontinuous deformation analysis (3D DDA). Geotech. Eng, 2004, vol. 8, iss. 5, pp. 491-496.
47. Hahn M., Wallmersperger T., Kroplin B.-H. Discrete element representation of discontinua: proof of concept and determination of material parameters. Comp. Mat. Sci, 2010, vol. 50, pp. 391-402.
48. Garagash I.A., Nikolaevskiy V.N. Non-associated laws of plastic flow and localization of deformation. Adv. Mech., 1989, vol. 12, iss. 1, pp. 131-183.
49. Stefanov Yu.P. Deformation localization and fracture in geomate-rials. Numerical simulation. Phys. Mesomech, 2002, vol. 5-6, pp. 67-77.
50. Wilkins M.L. Computer simulation of dynamic phenomena. Heidelberg: Springer-Verlag, 1999. 246 p.
51. Kushch V.I., Shmegera S.V., Sevostianov I. SIF statistics in micro cracked solid: effect of crack density, orientation and clustering. Int. J. Eng. Sci, 2009, vol. 47, pp. 192-208.
52. Kwon O., Kronenberg A.K., Gangi A.F., Johnson B. Permeability of Wilcox shale and its effective pressure law. J. Geophys. Res, 2001, vol. 106, no. b9, pp 19339-19353.
53. Paterson M.S., T.F. Wong. Experimental Rock Deformation. The Brittle Field. Berlin-Heidelberg: Springer-Verlag, 2005. 347 p.
54. Yamaji A. An Introduction to Tectonophysics: Theoretical Aspects of Structural Geology. Tokyo: TERRAPUB, 2007. 378 p.
55. Hangin J., Hager R.V., Friedman M., Feather J.N. Experimental Deformation of Sedimentary Rocks Under Confining Pressure: Pore Pressure Tests. AAPG Bulletin, 1963, vol. 47, iss. 5, pp. 717-755.
56. Stavrogin A.N., Tarasov B.G. Eksperimentalnaya fizika i mekha-nika gornykh porod [Experimental physics and mechanics of rocks]. Saint Petersburg: Nauka, 2001. 343 p.
57. Robin P.-Y.F. Note on Effective Pressure. J. Geophys. Res, 1973, vol. 78, no. 14, pp. 2434-2437.
58. Gangi A.F., Carlson R.L. An asperity-deformation model for effective pressure. Tectonophysics, 1996, vol. 256, pp. 241-251.
59. Boitnott G.N., Scholz C.H. Direct Measurement of the Effective Pressure Law: Deformation of Joints Subject to Pore and Confining Pressure. J. Geophys. Res, 1990, vol. 95, no. B12, pp. 19279-I9298.
60. Khavkin A.Ia. Nanoiavleniia i nanotekhnologii v dobyche nefti i gaza [Nano-effects and nanotechnology in oil and gas mining]. Moscow; Izhevsk: Reguliarnaya i khaoticheskaia dinamika, 2010. 692 p.
61. Basniev K.S., Kochina I.N., Maksimov V.M. Podzemnaya gidro-mekhanika [Underground fluid mechanics]. Moscow: Nedra, 1993. 416 p.
62. Leibenzon L.S. Dvizhenie prirodnykh zhidkostei i gazov v poristoi srede [A motion of natural liquids and gases in a porous medium]. Moscow-Leningrad: Gosudarstvennoe izdatelstvo tekhniko-teoreticheskoi literatury, 1947. 244 p.
Об авторах
Димаки Андрей Викторович (Томск, Россия) - кандидат технических наук, научный сотрудник лаборатории компьютерного конструирования материалов Института физики прочности и материаловедения СО РАН, старший научный сотрудник Томского государственного университета (634021, г. Томск, пр. Академический 2/4, e-mail: dav@ispms.tsc.ru).
Шилько Евгений Викторович (Томск, Россия) - доктор физико-математических наук, ведущий научный сотрудник лаборатории компьютерного конструирования материалов Института физики прочности и материаловедения СО РАН, профессор Томского государственного университета (634021, г. Томск, пр. Академический 2/4, e-mail: shilko@ispms.tsc.ru).
Астафуров Сергей Владимирович (Томск, Россия) - кандидат физико-математических наук, научный сотрудник лаборатории компьютерного конструирования материалов Института физики прочности и материаловедения СО РАН, старший научный сотрудник Томского государственного университета (634021, г. Томск, пр. Академический 2/4, e-mail: svastafurov@gmail.com).
Коростелев Сергей Юрьевич (Томск, Россия) - кандидат физико-математических наук, научный сотрудник лаборатории компьютерного конструирования материалов Института физики прочности и материаловедения СО РАН (634021, г. Томск, пр. Академический 2/4, e-mail: sergeyk@ispms.tsc.ru).
Псахье Сергей Григорьевич (Томск, Россия) - член-корреспондент РАН, директор Института физики прочности и материаловедения СО РАН, профессор Томского государственного университета (634021, г. Томск, пр. Академический 2/4, e-mail: svastafurov@gmail.com).
About the authors
Andrey V. Dimaki (Tomsk, Russian Federation) - Ph.D. in Technical Sciences, Researcher, Laboratory of computer-aided design of materials, Institute of Strength Physics and Materials Science Siberian Branch of Russian Academy of Science, Senior researcher of Tomsk State University (2/4, Akademicheskii av., 634021, Tomsk, Russian Federation, e-mail: dav@ispms.tsc.ru).
Evgeny V. Shilko (Tomsk, Russian Federation) Doctor of Physical and Mathematical Sciences, Prof., Principal researcher, Laboratory of computer-aided design of materials, Institute of Strength Physics and Materials Science Siberian Branch of Russian Academy of Science, Professor of Tomsk State University (2/4 Akademicheskii av., 634021, Tomsk, e-mail: shilko@ispms.tsc.ru).
Sergei V. Astafurov (Tomsk, Russian Federation) - Ph.D. in Physical and Mathematical Sciences, Researcher, Laboratory of computer-aided design of materials, Institute of Strength Physics and Materials Science Siberian Branch of Russian Academy of Science, Senior researcher of Tomsk State University (2/4, Akademicheskii av., 634021, Tomsk, Russian Federation, e-mail: svastafurov@gmail.com).
Sergei Yu. Korostelev (Tomsk, Russian Federation) - PhD in Physical and Mathematical Sciences, Senior Researcher, Laboratory of Computer-Aided Design of Materials, Institute of Strength Physics and Materials Science, Siberian Branch of the Russian Academy of Sciences (2/4, Akademicheskii av., 634021, Tomsk, Russian Federation, e-mail: sergeyk@ispms.tsc.ru).
Sergei G. Psakhie (Tomsk, Russian Federation) - Corresponding Member of the Russian Academy of Sciences, Director of Institute of Strength Physics and Materials Science, Siberian Branch of the Russian Academy of Sciences, Professor of Tomsk State University (2/4, Akademicheskii av., 634021, Tomsk, Russian Federation, e-mail: sp@ispms.tsc.ru).
Получено 28.08.2014
Просьба ссылаться на эту статью в русскоязычных источниках следующим образом:
Димаки А.В., Шилько Е.В., Астафуров С.В., Коростелев С.Ю., Псахье С.Г. Связанная модель флюидонасыщенных пористых материалов, основанная на совмещении дискретного и континуального подходов // Вестник Пермского национального исследовательского политехнического университета. Механика. -2014. - № 4. - С. 68-101. DOI: 10.15593/perm.mech/2014.4.03
Please cite this article in English as:
Dimaki A.V., Shilko E.V., Astafurov S.V., Korostelev S.Yu., Psakhie S.G. Coupled model of fluid-saturated porous materials based on a combination of discrete and continuum approaches. PNRPUMechanics Bulletin. 2014. No. 4. Pp. 68-101. DOI: 10.15593/perm.mech/2014.4.03
