الگوريتم تبريد شبيه‌سازي شده‌ انطباق‌پذير و سريع براي مدل‌سازي شبكه شكاف‌ها در مخازن شكاف‌دار طبيعي

An Adaptive and Fast Simulated Annealing Algorithm for Fracture Network Modeling in Naturally Fractured Reservoirs

دانستن جهت و توزیع فضایی شكاف‌ها در مخازن شكاف‌دار برای پیش‌بینی جریان سیالات در این نوع مخازن امری ضروری است. در میان روش‌های موجود برای مدل‌سازی توزیع شكاف‌ها، روش بهینه‌سازی تبرید شبیه‌سازی شده به دلیل توانایی آن در حل مسائل بزرگ و نیز یافتن مقادیر كمینه مطلق از اهمیت خاصی برخوردار است. در این مطالعه، هدف ارائه یك روش تبرید شبیه‌سازی شده انطباق‌پذیر و سریع، برای مدل‌سازی شكاف‌ها می‌باشد. این روش با ارائه‌ یك مدل برای تخمین مقدار اولیه پارامتر شبه دما، پیشنهاد یك مدل انطباق‌پذیر برای طول زنجیر ماركوف و در نهایت پیشنهاد یك مدل انطباق‌پذیر و سریع برای كاهش مقدار پارامتر شبه دما، باعث افزایش سرعت محاسبات الگوریتم تبرید شبیه‌سازی شده می‌شود. همچنین در این مطالعه، یك تابع هدف خاص مورد بررسی قرار داده شده و این نكته نشان داده می‌شود كه چرا الگوی شكاف‌های تولید شده با این روش منجر به دو دسته شكاف عمود برهم می‌شود. در ادامه تابع هدف بهبود داده شده كه بتواند الگوی شكاف‌های مزدوج با هر زاویه دلخواهی را تولید كند.

Knowing the orientation and spatial configuration of fractures in naturally fractured reservoirs is necessary for the prediction of the flow in those reservoirs. Among the fracture network modeling methods, simulated annealing algorithm (SA), due to its capability in solving large problems and finding a global optimum, is well-known. This paper proposes a new adaptive and fast simulated annealing algorithm for modeling naturally fractured reservoirs. This algorithm improves computation performance without degrading solution quality by incorporating a method for the estimation of the initial value of the temperature like parameter, T0, using an adaptive Markov chain length, NT (inner iterations) and suggesting a new fast and adaptive annealing schedule. Moreover, we discuss some aspects of a special objective function and proof this issue that why the minimum of this objective function occurs when the final configuration is a two orthogonal fracture sets with the same number of fractures in each set. Finally, we extend this objective function and modify it to develop an objective function that could generate conjugate fractures with any arbitrary intersection angle.