Редактирование: Эффективное решение задачи гидроразрыва с использованием модифицированной постановки на примере модели ХГД

Перейти к: навигация, поиск

Внимание! Вы не авторизовались на сайте. Ваш IP-адрес будет публично видимым, если вы будете вносить любые правки. Если вы войдёте или создадите учётную запись, правки вместо этого будут связаны с вашим именем пользователя, а также у вас появятся другие преимущества.

Правка может быть отменена. Пожалуйста, просмотрите сравнение версий, чтобы убедиться, что это именно те изменения, которые вас интересуют, и нажмите «Записать страницу», чтобы изменения вступили в силу.
Текущая версия Ваш текст
Строка 145: Строка 145:
  
 
Было проведено исследование зависимости решения от индекса поведения жидкости <math> n </math>.  
 
Было проведено исследование зависимости решения от индекса поведения жидкости <math> n </math>.  
Величина обратная индексу поведения жидкости входит в уравнения Пуазейля и скорости как показатель степени, поэтому при стремлении <math> n </math> к нулю, можно ожидать снижения точности решения. Погрешность вычислений для разных <math> n </math> представлена на рисунке [[Файл:Accuracy_dif_n.jpg|400px|thumb|right|Точность решения в момент времени  <math> t_{fin} = 100 </math>  в зависимости от индекса поведения жидкости  <math> n </math>]]. Из графика, представленного на рисунке, видно, что погрешность раскрытия в источнике существенно зависит от индекса поведения жидкости, а при его стремлении к нулю погрешность увеличивается в несколько раз. Эту погрешность можно устранить, перестроив алгоритм. Однако, это лежит в стороне от главной темы этой работы. Кроме того, в практике гидроразрыва обычно применяются жидкости с индексом поведения, превышающим <math> 0.5 </math> , т.е. из диапазона, в котором метод обеспечивает достаточную точность.
+
Величина обратная индексу поведения жидкости входит в уравнения Пуазейля и скорости как показатель степени, поэтому при стремлении <math> n </math> к нулю, можно ожидать снижения точности решения. Погрешность вычислений для разных <math> n </math> представлена на рисунке [[Файл:Accuracy_dif_n.jpg|400px|thumb|right|Точность решения в момент времени  <math> t_{fin} = 100 </math>  в зависимости от индекса поведения жидкости  <math> n </math> $]]. Из графика, представленного на рисунке, видно, что погрешность раскрытия в источнике существенно зависит от индекса поведения жидкости, а при его стремлении к нулю погрешность увеличивается в несколько раз. Эту погрешность можно устранить, перестроив алгоритм. Однако, это лежит в стороне от главной темы этой работы. Кроме того, в практике гидроразрыва обычно применяются жидкости с индексом поведения, превышающим <math> 0.5 </math> , т.е. из диапазона, в котором метод обеспечивает достаточную точность.
 
 
==Выводы==
 
 
 
# Установлено, что при использовании конечных разностей и квадратурных формул, учитывающих специфику задачи, задача ХГД может быть решена в глобальных координатах без обращения гиперсингулярного оператора.
 
# Разработанный подход, использующий удвоение размера сетки при достижении фронтом заранее заданной границы, позволяет не пересчитывать матрицу дискретизированного упругого оператора на шагах интегрирования по времени. Новые коэффициенты влияния получаются из элементов матрицы делением на масштабный фактор.
 
# В отличии от большинства предшествующих работ, метод не требует обращения матриц. В нем используются только скалярные произведения. Это благоприятствует объединению его с быстрым методом мультиполей.
 
# Полученная динамическая система можеть быть эффективно решена с помощью методов решения задачи Коши для ОДУ, например, методом Рунге-Кутты, Адамса и др.
 
# Установлено, что в рассматриваемой задаче метод Адамса по схеме предиктор-корректор имеет преимущество, поскольку, будучи более устойчивым, требует меньших вычислительных затрат при сохранении точности получаемого решения.
 
# В случае, когда не привлекаются специальные асимптотические методы, отвечающие близкому к нулю индексу поведения жидкости, точность решения падает. Например, для сетки с <math>M+1=41 </math> узлом погрешность вычислений раскрытия в источнике и полудлины трещины составляет соответственно:
 
#* <math> 7.7\% </math> и <math> 4.8\% </math> для индекса поведения жидкости <math> n=0.3 </math>,
 
#* <math> 10.9\% </math> и <math> 5.7\% </math> для индекса поведения жидкости <math> n=0.1 </math>. Тем не менее в практически важном интервале $ 0.5 < n < 1 $ падение точности не происходит.
 
# Разработанный метод, будучи избавленным от упрощений, доступных только в одномерных задачах, допускает распространение на трехмерную задачу. Его использование для решения задачи ХГД на прямоугольной сетке дало результаты с точностью не меньшей, чем в данной работе.
 
 
 
Направление дальнейших исследований --- распространение метода на трехмерную задачу и применение более эффективных, чем метод Адамса, методов численного интегрирования, например, формулы дифференцирования назад.
 
 
 
==Список литературы==
 
* Adachi J. et al. Computer simulation of hydraulic fractures //International Journal of Rock Mechanics and Mining Sciences. – 2007. – Т. 44. – №. 5. – С. 739-757.
 
 
* Geertsma, J., F. De Klerk. A rapid method of predicting width and extent of hydraulically induced fractures //Journal of Petroleum Technology. – 1969. – Т. 21. – №. 12. – С. 1,571-1,581.
 
*Khristianovic S., Zheltov Y. Formation of vertical fractures by means of highly viscous fluids //Proc. 4th world petroleum congress, Rome. – 1955. – Т. 2. – С. 579-586.
 
* Linkov A. M. Speed equation and its application for solving ill-posed problems of hydraulic fracturing //Doklady Physics. – MAIK Nauka/Interperiodica, 2011. – Т. 56. – №. 8. – С. 436-438.
 
*Linkov A. M. On efficient simulation of hydraulic fracturing in terms of particle velocity //International Journal of Engineering Science. – 2012. – Т. 52. – С. 77-88.
 
*Linkov A. M., Mishuris G. Modified formulation, $ \varepsilon $-regularization and the efficient solution of hydraulic fracture problems //ISRM International Conference for Effective and Sustainable Hydraulic Fracturing. – International Society for Rock Mechanics, 2013.
 
* Linkov A. M. The particle velocity, speed equation and universal asymptotics for the efficient modelling of hydraulic fractures //Journal of Applied Mathematics and Mechanics. – 2015. – Т. 79. – №. 1. – С. 54-63.
 
*Pierce A. P., Siebrits E. A dual multigrid preconditioner for efficient solution of hydraulically driven fracture problem //International Journal of Numerical Methods and Engineering. – 2005. – Т. 65. – С. 1797-1823.
 
*Peirce A., Detournay E. An implicit level set method for modeling hydraulically driven fractures //Computer Methods in Applied Mechanics and Engineering. – 2008. – Т. 197. – №. 33. – С. 2858-2885.
 
*Peirce A. Modeling multi-scale processes in hydraulic fracture propagation using the implicit level set algorithm //Computer Methods in Applied Mechanics and Engineering. – 2015. – Т. 283. – *Rice J. R. Mathematical analysis in the mechanics of fracture //Fracture: an advanced treatise. – 1968. – Т. 2. – С. 191-311.
 
*Settari A., Cleary M. P. Development and testing of a pseudo-three-dimensional model of hydraulic fracture geometry //SPE Production Engineering. – 1986. – Т. 1. – №. 06. – С. 449-466.
 
*Mack M. G., Warpinski N. R. Mechanics of hydraulic fracturing //Reservoir stimulation. – 2000. – С. 6-1.
 
*Алексеенко О. П. и др. Двумерная пошаговая модель распространения трещины гидроразрыва //Вестник НГУ. Серия: Математика, механика, информатика. – 2011. – Т. 11. – №. 3. – С. 36-59.
 
*Линьков А. М. Решение осесимметричной задачи о гидроразрыве для утончающихся жидкостей // ПММ. -2016. - Т. 80. - №. 2. - С. 207-217.
 
*Линьков А. М. Численное решение плоской задачи о гидроразрыве в модифицированной постановке при произвольных начальных условиях //Физико-технические проблемы разработки полезных ископаемых. -2016. -№. 2.
 
*Самарский А. А., Гулин А. В. Численные методы. – Москва: Наука, 1989.
 
*http://vseonefti.ru/upstream/frac.html
 
Вам запрещено изменять защиту статьи. Edit Создать редактором

Обратите внимание, что все добавления и изменения текста статьи рассматриваются как выпущенные на условиях лицензии Public Domain (см. Department of Theoretical and Applied Mechanics:Авторские права). Если вы не хотите, чтобы ваши тексты свободно распространялись и редактировались любым желающим, не помещайте их сюда.
Вы также подтверждаете, что являетесь автором вносимых дополнений или скопировали их из источника, допускающего свободное распространение и изменение своего содержимого.
НЕ РАЗМЕЩАЙТЕ БЕЗ РАЗРЕШЕНИЯ МАТЕРИАЛЫ, ОХРАНЯЕМЫЕ АВТОРСКИМ ПРАВОМ!

To protect the wiki against automated edit spam, we kindly ask you to solve the following CAPTCHA:

Отменить | Справка по редактированию  (в новом окне)