MODFLOW 6 آخرین نسخه کد مدل سازی زیرزمینی MODFLOW توسط USGS است. این نسخه شامل برخی از ابزارهای نوآورانه و یک شکل مجدد کامل ساخت مدل است، اما تا به امروز این کد هیچ رابط کاربری گرافیکی (GUI) تجاری / باز نداشته است که باعث بهبود استفاده از کد برای طراحان مبتدی و مبتکران علوم زمین می شود.

ماژول Gempy به عنوان کتاب منبع باز پایتون برای تولید مدل های زمین شناسی ساختاری کامل 3D است. این کتابخانه کامل برای ایجاد مدل های زمین شناختی از interfaces، گسل ها و جهت گیری های لایه است، همچنین توالی لایه های زمین شناسی را برای نشان دادن نفوذ و خطاهای سنگ انجام می دهد.

الگوریتم برای مدل سازی زمین شناسی مبتنی بر واسنجی cokriging جهانی با حمایت از کتابخانه های ریاضی Python تحت عناوین Numpy، PyMC3 و Theano است.

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

به منظور تخمین دقیق تابع هدف، مدلسازی آب های زیرزمینی با MODFLOW و سایر کدها به صورت مدلسازی معکوس تعریف می شود که در آن پارامترهای کلیدی از طریق مقایسه نتایج مدل با داده های مشاهده شده محاسبه می شود. این فرایند مقایسه طولانی است، استفاده از معیارهای پذیرش و تحلیل روند متاثر شرایط مرزی است.

با توجه به گسترش تغییرات شرایط آب و هوایی، تدارکات ذخیره سازی (Tailings (TSF در تعامل با جریان آب های زیرزمینی با رفتار دوگانه از مناطق تخلیه و شارژ است. ارزیابی بلندمدت اثرات زیست محیطی از زباله ها باید تاثیر بالقوه ای را از نفوذ آب شور به رژیم جریان آب زیرزمینی و تأثیر آن روی کیفیت آب در کانال های پایین دست بررسی نماید. مدل سازی عددی یک ابزار قوی برای ارزیابی جریانها و کیفیت در ناحیه تاثیر کلی از سدها در شرایط فعلی و آینده است.

MODPATH یک مدل ردیابی ذرات است که بر روی مدل های جریان آب زیرزمینی مدلسازی شده در MODFLOW اجرا می شود. ردیابی ذرات سه بعدی است و می تواند گزینه های گوناگونی را در جهت مدلسازی رو به جلو / عقب بین سایرین داشته باشد. مدل Muse قبل و بعد از پردازش مدل های عددی در MODFLOW، PHAST و SUTRA است، این نرم افزار می تواند برخی از مدل های مرتبط MODFLOW را به عنوان MODPATH و MT3D (ردیابی ذرات و حمل و نقل جمعی) شبیه سازی کند. این آموزش روش اولیه را برای راه اندازی یک مدل ردیابی ذرات با MODPATH با استفاده از مدل Muse بر یک مدل MODFLOW پایه نشان می دهد.

مدل سازی منطقه ای آب های زیرزمینی یک وظیفه مهم در مدیریت آب استراتژیک است که شامل تمام کاربران، فعالیت ها و اکوسیستم های مربوطه و استفاده پایدار از شرایط فعلی و آینده است. برخی از ملاحظات خاص در مورد مدلسازی منطقه ای با توجه به تقسیم بندی اولیه و فضایی وجود دارد؛ یک مدل منطقه ای در نظر گرفته نشده است که پاسخ منطقه آبخوان را برای منطقه مشخص تعیین کند، در عوض شامل ارزیابی جریان آب های زیرزمینی منطقه و اندازه گیری شارژ، تخلیه و سایر فرآیندهای موجود در تعادل آب است.

SEAWAT یک مدل توسعه یافته توسط USGS برای شبیه سازی چگالی جریان آب زیرزمینی سه بعدی با انتقال شوری و گرما است. این نرم افزار بر اساس MODFLOW-2000 و MT3DMS است و در آخرین نسخه آن می تواند تغییرات ویسکوزیته را شبیه سازی کند و زمان اجرای سریع تر را فراهم کند. SEAWAT در Flopy، کتابخانه پایتون برای ساخت، اجرا و نمایش مدل MODFLOW اجرا می شود. این آموزش دارای جریان کاری کامل برای ایجاد و ارائه مثال پایه ای از ابزار سالین با SEAWAT و Flopy در یک نوت بوک Jupyter است.

نشان دادن مسیر جریان آب زیرزمینی برای درک شرایط واقعی و پیش بینی شده رژیم جریان آب زیرزمینی مفید است. این جهت و بزرگنمایی جهت و  چشم انداز سریع جریان های اصلی آب زیرزمینی و بین خطوط منابع و نقاط تخلیه را نشان می دهد. بردارهای جریان ویژگی های بسیار رایج در مدل سازی آب های زیرزمینی هستند که با استفاده از پشته منبع باز تا زمانی که توسعه کتابخانه Flopy وجود ندارد، دارای یک ابزار ویژه برای ارائه است.

آبخوان ها می توانند متخلخل، شکسته یا کارستی باشد. با توجه به تنظیمات زمین شناسی و فرآیندهای مربوط به تشکیل رسانه متخلخل، این مورد می تواند درجه بالایی ناهمگونی را نشان دهد که بر الگوهای آب های زیرزمینی و انتقال و توزیع آلاینده تأثیر می گذارد.

آلاینده ها برای منابع puntual / aerial با ناهمگونی ها و ناهنجاری های آبخوان روبرو می شوند. شناخت و مفهوم سازی توزیع واحدهای هیدروژئولوژیکی مختلف و خواص آنها در زیر زمین یک چالش برای هیدروژئولوژیست، مدلسازان عددی و متخصصان اصلاح است.

ساخت مدل های هیدروژئولوژیکی می تواند وقت زیادی برای ساخت و ساز، تجسم نتایج و کالیبراسیون را استفاده کند. مهم است که از ابزارهایی استفاده کنید که می توانند این کارها را بهینه سازی کنند و در زمان صرفه جویی کرده و در تجزیه و تحلیل سیستم آن را به کار گیرید.

در این فرصت ما Flopy را برای تکرار یک مدل حمل و نقل 2D از یک پست قبلی استفاده خواهیم کرد  با مجموعه ای متنوع از اسکریپت های پایتون است که می تواند برای اجرای MODFLOW و MT3D، از جمله برنامه های زیرزمینی مرتبط با MODFLOW، به روش ساده و کارآمد مورد استفاده قرار گیرد.

کادر بسته Package MODFLOW نشان دادن گزینه گزینه.

بسته MODFLOW Subidence) SUB) ابزاری برای شبیه سازی فرونشست آبخوان و تخریب زمین است. بسته SUB شامل توانایی شبیه سازی تراکم در هر دو شرایط الاستیک (قابل انعطاف) و غیر قابل الاستیک (غی قابل انعطاف) را داراست. این همچنین شامل توانایی شبیه سازی interbeds است که در آن زهکشی از interbed فوری (بدون تأخیر) یا با تأخیر است. بسته MOD MODEFLOW SUB می تواند یک پروژه جابجایی عمودی باینری ایجاد کند که توسط GMS خوانده شده و به عنوان یک مجموعه داده در موتور حل MODFLOW در Project Explorer نمایش داده شود. شرح مفصلی ورودی SUB در وب سایت USGS موجود است.

برای Visual MODFLOW Flex 5.1 آماده باشید

لذت بردن از نسخه جدید ماژول نرم افزاری حمل و نقل جریان های آب زیرزمینی و حمل و نقل آلودگی. در ژوئن 27 برای ارتقا در دسترس خواهد بود.

به روز رسانی های جدید در 5.1 شامل موارد زیر خواهد بود:

آزمایش های آبخوان (آزمایش های پمپاژ، آزمایش های اسلاگ و آزمون های سطح-ثابت) برای تخمین مقادیر مشخصه مکان برای خواص هیدرولیک آبخوان ها و آکیتاردها انجام می شود. با این حال، در شرایط خاص، اطلاعات مربوط به ویژگی های هیدرولیکی یک مکان ممکن است در صورت لزوم در دسترس نباشد. به عنوان مثال، مطالعات شناسایی یا محاسبات محدوده ممکن است نیاز به ارزش های هیدرولیکی قبل از تحقیق در محل داشته باشد.

در بخش های زیر ارزش های مقادیر موجود ضرایب هیدرولیکی برای هدایت هیدرولیکی افقی و عمودی، ضریب نگهداری، آبدهی ویژه و تخلخل ارائه شده است. اگر این اطلاعات مربوط به مکان مطالعه شما در دسترس نیست یا برای بررسی نتایج آزمایشات آزمایشگاهی به چنین اطلاعاتی نیازمند هستید، به این مقادیر در جداول زیر مراجعه کنید.

فرمولاسیون جریان XT3D

فرمول جدید جریان آب زیرزمینی به طور خاص برای مدل GWF در MODFLOW 6 طراحی شده است. این فرمول جدید جریان، گزینه XT3D نامیده می شود. گزینه XT3D توانایی های MODFLOW را گسترش می دهد و امکان شبیه سازی بی نهایت کامل سه بعدی را در شبکه های منظم یا نامنظم فراهم می کند به طوری که به درستی تانسور هدایت هیدرولیکی کامل سه بعدی در نظر گرفته شود. همچنین می تواند دقت شبیه سازی جریان آب زیرزمینی را در مواردی که شبکه مدل نقض الزامات خاص هندسی را نقض می کند، بهبود دهد. بنابراین، گزینه XT3D جایگزینی برای (Packet Correction Ghost-Node Correction (GNC است که برای MODFLOW-USG طراحی شده است.

این شکل یک شبکه مثلثی (مش) را نشان می دهد که در آن اندازه های سلول های مثلثی در مناطقی با شیب هیدرولیکی نسبتا بزرگ مانند اطراف خط ساحلی دریاچه، چاه های پمپاژ و در امتداد یک جریان کاهش می یابد. این نوع شبکه لایه ای را می توان با استفاده از Discretization by Vertices Package در MODFLOW 6. ارائه کرد.

MODFLOW یک مدل جریان آب زیرزمینی منبع باز است که توسط سازمان تحقیقات زمین شناسی ایالات متحده توزیع شده است. برای بیش از 30 سال، برنامه MODFLOW به طور گسترده ای توسط دانشگاهیان، مشاوران خصوصی و دانشمندان دولتی به طور دقیق، قابل اعتماد و موثر برای شبیه سازی جریان آب های زیرزمینی استفاده شده است. با توجه به رشد روزافزون در تعاملات سطحی و آبهای زیرزمینی، پالایش محلی با شبکه های غوطه ور و غیر ساختاری، جریان آبهای زیرزمینی کارست، حمل و نقل محلول و نفوذ آب شور، منجر به توسعه نسخه های متعدد MODFLOW شده است. اگر چه این نسخه های MODFLOW اغلب بر اساس نسخه اصلی MODFLOW (که قبلا MODFLOW-2005) بود، اغلب ناسازگاری هایی است که استفاده آنها را با سایر نسخه های MODFLOW محدود می کند. در بسیاری از موارد، توسعه این نسخه های MODFLOW جایگزین با توجه به ساختار برنامه پایه، برای شبیه سازی یک مدل جریان جریان آب زیرزمینی با استفاده از یک شبکه منظم MODFLOW که شامل لایه ها، ردیف ها و ستون ها بود، به چالش کشیده شد.

گاها در هنگام مدل سازی و در مرحله واسنجی (نوعی بهینه سازی)، با خطای Floating Overflow برخورد می کنید. علت اصلی این خطا در نرم افزارها، تلفیق ارقام نامعقول پارامتری برای معادلات و در اینجا ماتریس های تفاضل محدود و یا الان محدود است. به این معنی که در توسعه مدل مفهومی، ارقام وارد شده ابتدایی برای پارامترهای آبخوان به شدت ناهمساز و نامرتبط بوده اند. 

پس پردازش نتایج

یک بار دیگر، می توانیم از فایل خروجی باینری MODFLOW، با استفاده از ماژول flopy.utils.binaryfile، سطوح آب را بخوانیم. همراه با Object HeadFile چندین روش وجود دارد که ما در اینجا استفاده خواهیم کرد: get_times لیستی از بارهای موجود در فایل head file را بازگرداند. get_data آرایه سطح سه بعدی را برای زمان مشخص بازگرداند. get_ts آرایه سری زمانی [ntimes، headval] برای سلول مشخص شده را برمیگرداند.

با استفاده از این روش، می توانیم نتایج حاصل از مدل ها، نمونه های سطح آب و هیدروگراف را ایجاد کنیم:

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

Packet Transient Well

حالا ما می توانیم شیء بسته چاه را ایجاد کنیم که از نوع flopy.modflow.ModflowWel است.

تعریف دوره های تنش

برای ایجاد مدل با دوره های تنش های مختلف، ما نیاز به تعریف nper، perlen، nstp و steady داریم. این کار در بلوک زیر انجام می شود به طوری که به ما اجازه می دهد این متغیر را مستقیما به شیء تقسیم بندی منتقل کنیم:

در این مثال، ما مدل آموزش های قبلی را به یک مدل جریان غیرقابل کنترل و مداوم با مرزهای زمانی متفاوت تبدیل خواهیم کرد. به جای استفاده از سطح آب های ثابت برای مرزهای چپ و راست (با تنظیم ibound to -1)، ما از مرزهای سطح آب کلی استفاده خواهیم کرد. ما در این مدل باید شرایط زیر را در نظر داشته باشید:

شرایط اولیه - سطح آب 10.0 در همه جا یکسان است.

دوره 1 (1 روز) - حالت پایدار با مرحله چپ و راست GHB = 10.

دوره 2 (100 روز) - GHB سمت چپ با مرحله = 10، راست GHB با مرحله تعیین شده به 0.