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

شکست هیدرولیکی روشی نوین و پرکاربرد جهت استحصال ذخایر و منابع انرژی موجود در اعماق زمین است. با توجه به افزایش مصرف انرژی از یک سو و اتمام ذخایر انرژی از سوی دیگر، در آینده‌ای نه چندان دور استفاده از این روش به عنوان امری ضروری مبدل خواهد شد. پارامترهای مهم زیادی در این روش  ﺗﺃثیرگذار هستند که از جمله‌ی آن‌ها می‌توان به فشار و نحوه‌ی اعمال آن به منظور ایجاد شکستگی و پیشرفت درزه درون لایه‌های سنگی‌ اشاره نمود. پارامتر مهم دیگر در این روش، اندرکنش میان شکستگی‌های طبیعی از پیش موجود با زوایای مختلف و شکست هیدرولیکی ایجاد شده، می‌باشد. با توجه به هزینه‌های بالای این فرآیند، هدف از انجام این پژوهش رسیدن به بهینه ‌ترین حالت ممکن برای بیشترین میزان پیشرفت شکست هیدرولیکی و در عین حال استفاده از کمترین مقدار فشار شکست است که با استفاده از مدل‌سازی عددی به وسیله‌ی نرم افزار المان مجزاء pfc از شرکت آیتسکا به صورت دو بعدی بر روی نمونه سنگ‌های پوچئون گرانیت با رفتار شکننده و مقاومت بالا انجام می‌شود. برای ساخت نمونه به روش المان مجزا در نرم‌افزار pfc از دیسک‌های دایره‌ای در دو بعد استفاده می‌شود که این ذرات به وسیله پیوندهایی با یکدیگر در تماس هستند و برای اعمال بارگذاری بر روی نمونه‌ دیوارهایی وجود دارند که با این دیسک‌ها در تعامل می‌باشند. به منظور ایجاد پیوند میان ذرات از مدل باند مسطح و برای ایجاد شکستگی‌های طبیعی از پیش موجود از روش شبکه شکستگی استفاده می‌شود. با توجه به نتایج به دست آمده مشخص گردید که با افزایش نرخ اعمال بارگذاری، نمونه در زمان کمتری به فشار شکست رسیده و شکسته می شود ولی مقدار فشار شکست افزایش می یابد. هم‌چنین با افزایش زاویه شکست طبیعی نسبت به افق (در جهت ساعتگرد)، نمونه در فشار شکست کمتری شکسته می شود. در نهایت هم با افزایش فاصله شکست طبیعی از مرکز نمونه، تاثیر حضور درزه در نمونه کاهش می‌یابد و فشار شکست به حالت بدون درزه نزدیک می شود. 

numerical modeling of hydraulic fracture considering the pre-existing joints using pfc software

hydraulic fracturing is a new and widely used method for extracting reserves and energy resources in the depths of the earth. in the near future, due to increase in energy consumption on the one hand and depletion of energy reserves on the other, using of this method will become a necessity. one of the most important and effective parameters in this process is the pressure and how it is applied in order to create fractures and fracture progression in rock layers. another important parameter is the interaction between pre-existing natural fractures with different angles and hydraulic fracture.  due to the high costs of this process, the purpose of this study is to achieve the optimal state for the maximum progress of hydraulic fracture and the lowest amount of breakdown pressure at the same time. numerical modeling was performed in two dimensions by particle flow code ( pfc ) software from itasca company on samples of pocheon granite rocks with brittle behavior using distinct element method. pfc software uses circular disks in two dimensions to construct and make the sample using distinct element method. these particles are in contact with each other through bonds. in this program there are walls that interact with these disks to apply load on the sample. the flat joint model is used in order to create contacts between particles, and discrete fracture network is uded to to create pre-existing natural fractures for interacting with hydraulic fracture.given that in the distinct element method (pfc software) our sample consists of a large number of particles in two dimensions that the general characteristics of the sample are formed based on the interaction between these particles, so we need parameters as input data to our software, existing disks and the link between them, to finally obtain the specifications of the same laboratory sample after modeling. these specifications and input data are referred to as micro parameters, and the final specifications, which are the same as our mechanical parameters in the laboratory, are referred to as macro parameters.to find the micro parameters of the sample we use the trial and error method. here, our modeling is based on brazilian and uniaxial compression experiments and … performed by zhuang et al. laboratory investigations.due to the limitations of pfc software for fluid flow modeling and limitations for using of cfd relationships, pressure equals to fluid pressure can be used as a new solution. in this way, by modeling a number of walls that form a complete circle with overlapping each other, and by considering the servo control mechanism, we move them in the opposite direction of their normal vectors and off-center, creating a comprehensive pressure which is actually same as the fluid pressure. from the obtained results, it was found that with increasing the loading rate, the sample reaches the breakdown pressure and breaks in less time, but the amount of breakdown pressure increases. also, by increasing the natural fracture angle relative to the horizon (clockwise), the specimen breaks at a lower breakdown pressure. finally, by increasing the natural fracture distance from the center of the sample, the effect of the presence of the joint in the sample decreases and the breakdown pressure approaches the seamless state.