ارزیابی روش‌های تکراری زیرفضای کریلف برای محاسبه جریان سیال در شبکه شکستگی‌های مجزای سه‌بعدی

محاسبه جریان سیال در محیط‌های سنگی درزه‌دار یکی از موضوعات بسیار مهم در مباحث تراوایی در مهندسی سنگ است. توده‌سنگ شامل شبکه درهم‌تنیده‌ای از ماده‌سنگ و شکستگی‌ها است. تعداد و الگوی اتصال‌پذیری شکستگی‌ها دو عامل کلیدی کنترل‌کننده جریان سیال در توده‌سنگ هستند. یکی از روش‌هایی که می‌تواند ساختار هندسی توده‌سنگ را به دقیق‌ترین شکل ممکن نشان دهد، روش شبکه شکستگی‌های مجزا (dfn) است. با توجه به‌این‌که ناهمگنی و ناهمسانگردی توده‌های سنگی در جهات مختلف می‌تواند در محاسبه جریان سیال بسیار اثرگذار باشد، روش‌ شبکه شکستگی‌های مجزای سه‌بعدی بیشتر در کانون توجه قرار دارد. محاسبه‌ی عددی جریان در توده‌سنگ مستلزم حل معادلات پرتعداد و پیچیده‌ای است که توسط روش‌های عددی نظیر روش المان محدود تولید می‌شوند. حل دستگاه‌های معادلات مذکور معمولاً ساده نیست و نیازمند روش‌های خاصی است که یکی از بهترین ‌آنها، روش‌های تکراری زیرفضای کریلف است. در این تحقیق پس از اعتبارسنجی محاسبات روش‌های مختلف زیرفضای کریلف با یک روش مستقیم و مدلسازی 3dec، ارزیابی کارایی این روش‌ها مورد بررسی قرارگرفته و کارآمدترین آنها از طریق یک سری تحلیل‌های حساسیت بر روی پارامتر‌های دقت و سرعت انجام محاسبات، تعیین شده است. با توجه به تحلیل‌های انجام شده در این تحقیق، روش‌های cg، cr و iom جزو سریع‌ترین و دقیق‌ترین روش‌های زیرفضای کریلف انتخاب شده‌اند. نتایج این پژوهش می‌تواند راهنمای خوبی برای پژوهشگرانی باشد که بر روی حل جریان سیال در توده‌سنگ متمرکز هستند.

Evaluation of Krylov Subspaces Iterative Methods for Calculating Fluid Flow in ThreeDimensional Discrete Fracture Networks

SummaryComputation of fluid flow in fractured rocks is very important. The rockmass is consisted of intact rock and fractures. Number and connectivity pattern of the fractures are two key factors controlling the fluid flow in the rockmasses. One of the most accurate methods to model geometrical structure of the rockmasses is discrete fracture network (DFN). Anisotropy and heterogeneity of the rock masses often affects the computations of the flow, therefore, threedimensional DFN has been more desirable in literatures. Numerical calculation of the fluid flow requires solving a large system of equations which are generated by discretization schemes. Solving these systems are not usually straightforward and it needs more special and complex methods to converge the result. One of the best methods in this regard are Krylov subspaces methods. Evaluation of different Krylov subspaces methods which have been validated in comparison with a direct method and 3DEC modeling, has been considered in this research and the most optimized methods have been determined using a series of sensitivity analyses. Therefore, CG, CR and IOM have been characterized as the most accurate and fastest Krylov subspaces methods. The provided results in this research can be a sufficient guideline for the researchers who want to study the fluid flow in fractured rocks. IntroductionIn this research, a numerical model is developed to calculate flow field in fractured rocks. The rockmass is modeled geometrically from crude mapping data and an optimized algorithm is chosen to triangulate the geometrical framework of the model. The flow filed is discretize using a numerical scheme and the generated system of equations are solved using different iterative methods. The application of different iterative solving methods has been validated in comparison with direct one and a series of sensitivity analysis is performed to determine the most optimal iterative methods. Methodology and Approaches3DDFN forms the geometrical framework of the present geometrical model as one of the most accurate and used methods to simulate fractured rocks. The algorithm provided by Erhel et. al. is the foundation of meshing process and FEM method is used to discretize the geometrical structure. The large system of equation generated by FEM is solved using different Krylov subspaces methods and the results are validated by LQ factorization. The sensitivity analysis is performed on two key parameters, precise and CPU runtime of the model to determine the most optimal method. Results and ConclusionsThe results of sensitivity analysis show that the methods CG, CR and IOM are the most optimal methods of Krylov subspaces which can be well compatible to calculation of fluid flow in 3DDFN models.