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.

اجرای مدل سازی

Flopy همچنین می تواند برای اجرای مدل مورد استفاده قرار گیرد. شی مدل (mf در این مثال) یک روش متصل دارد که مدل را اجرا می کند. برای این کار، برنامه MODFLOW باید در جایی در مسیر سیستم یا در داخل دایرکتوری کار قرار بگیرد. در این مثال، ما مشخص کرده ایم که نام برنامه اجرایی 'mf2005' است. برای اجرای این مدل زیر را امتحان کنید:

بسته پایه

در گام بعد ما می توانیم جسم FloPy را ایجاد کنیم که MODFLOW Basic Pack را نشان می دهد. جزئیات در کلاس BAS فلاپی در: flopy.modflow.mfbas موجود است. برای این مدل ساده، مقادیر سطح پایه 10 و 0 به ترتیب به ستون اول و آخر مدل (در همه لایه ها) اختصاص داده می شود. کد پایتون برای انجام این کار:

ایجاد مدل MODFLOW

یکی از موارد مفید در مورد ایجاد مدل در پایتون این است که فرآیندی بسیار آسان داشته که با تغییر یک یا دو جزء، به طور کلی شبکه برای مدل شما را تغییر می دهد. بنابراین در این مثال، ما اسکریپت پایتون را طراحی می کنیم تا تعداد لایه ها، ستون ها و ردیف ها به راحتی تغییر کند.

ما می توانیم مدل MODFLOW بسیار ساده ای را تهیه کنیم که یک بسته پایه (BAS)، دیسک ورودی دیسکاوری (DIS)، بسته LPF یا package property، کنترل خروجی (OC) و حل کننده موتور گرادیان (PCG) را ایجاد کند. هر یک از این موارد فایل ورودی خود را دارد، که به صورت خودکار بوسیله FloPy ایجاد می شود.

این مجموعه آموزشی در سایت بیسین نشان خواهد داد که استفاده از FloPy برای توسعه مدل MODFLOW ساده است. توجه داشته باشید که می توانید این اسکریپت آموزش داده شده در پایتون را از اینجا دریافت کنید.

شروع آموزش: مقدمه و نحوه فعال سازی

اگر FloPy به درستی نصب شده باشد، می توان آن را به صورت زیر وارد کرد:

FloPy چیست؟

بسته FloPy شامل مجموعه ای از اسکریپت های پایتون برای اجرای MODFLOW، MT3D، SEAWAT و دیگر برنامه های آب زیرزمینی مرتبط با MODFLOW است. FloPy شما را قادر می سازد تمام این برنامه ها را با اسکریپت های پایتون اجرا کنید. پروژه FloPy در سال 2009 آغاز شده و به یک مجموعه کامل از اسکریپت ها با پایگاه کاربری رو به رشد افزوده است. FloPy3 در دسامبر 2014 با چند پیشرفت عالی که FloPy3 را به مدل سازی پسرو سازگار می کند، منتشر شد. اولین تغییر قابل توجه این است که FloPy3 از همه جا در فهرست صفر استفاده می کند، بدین معنا که تمام لایه ها، ردیف ها، ستون ها و دوره های استرس شروع به صفر شدن می کنند. این تغییر برای سازگاری صورت گرفته بود، چون همه شاخصه های آرایه قبلا صفر بود (همانطور که همه آرایه ها در پایتون) بود. این ممکن است نسبتا به مقدار کمی به کار گرفته شده و استفاده شود، اما امیدوارم از آشفتگی در آتی جلوگیری شود. دومین تقویت قابل توجه در مورد توانایی های این کد، مشخص کردن شرایط مرزی متفاوت و متغیر است

PEST.cloud سایتی است که شما را قادر می سازد تا مدل خود را بر روی ابر مایکروسافت Azure با استفاده از PEST_HP کالیبره کنید. اجرای PEST در ابر آسان است. شما فقط باید مدل خود را آپلود کنید، تعدادی از عوامل محاسبه را انتخاب کنید و استقرار دهید. برای یک مرور کلی از آنچه انتظار می رود از رابط PEST.cloud، اینجا کلیک کنید.

PEST_HP یک نسخه از BEOPEST است که برای استفاده در محیط های بسیار موازی مورد استفاده قرار می گیرد. همچنین دارای تحمل بیشتری برای نقایص مدل است که عملا محاسبه مشتقات مختلفی را در فرآیند واسنجی به خطر می اندازد. دستورالعمل PEST_HP را بخوانید تا چیزهای دیگری را که می بینید بدانید. هنگام استفاده از PEST.cloud شما به طور خودکار از PEST_HP استفاده می نمایید.

مرکز مدل سازی یکپارچه آب های زیرزمینی

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

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

چه ویژگی جدیدی در PEST نسخه 14 است؟ چند رفع اشکال و برخی از بهبود کارایی های جزئی. با این حال، تغییر بزرگ سازگاری با PEST_HP - در نسخه جدید PEST است که طراحی آن برای استفاده در محیط های بسیار موازی مانند حالت ابری بهینه سازی شده است. شما می توانید PEST_HP را از وب سایت PEST.cloud دریافت کنید.

کد PEST بسته نرم افزاری استاندارد صنعت برای ارزیابی پارامترها و تجزیه و تحلیل عدم قطعیت مدل های پیچیده زیست محیطی و دیگر مدل های کامپیوتری است. کارایی این کد  به این گونه است که همراه با طیف وسیعی از برنامه های کاربردی که از آن در ارتباط با مدل های آب زیرزمینی و آب سطحی استفاده می شود؛

PEST - کتاب

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

هدایت هیدرولیکی رسوبات و لایه های زمین شناسی

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

پنجره زنجیره مارکف عمودی

پنجره کمکی در تعریف زنجیره مارکف عمودی. این پنجره از سه بخش اصلی تشکیل شده است:

بخش طرح

بخش زنجیره مارکوف 

بخش صفحات گسترده

هر سه بخش توسعه 1D زنجیره مارکوف در جهت عمودی را میسر می سازد.

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

سه روش برای تولید مدل سازی تصادفی در GMS  با استفاده از MODFLOW در ورژن  2000 یا 2005 وجود دارد. اول، پهنه بندی پارامتر با استفاده از هر یک از روش نمونه گیری تصادفی، نمونه برداری مکعب لاتین، و یا زمینه گاوسی برای تولید نتایج مختلف می باشد. روش دوم استفاده از شبیه سازی شاخص های تولید شده توسط ابزار T-PROGS می باشد. سوم، روش فضایی تهی مونت کارلو (NSMC) است؛ که به تولید مدل های متعدد با مجموعه های مختلف از پارامترها اقدام می کند (در GMS 9.0 و ورژن های بعد از آن).

پس از تولید نتایج شبیه سازی تصادفی، کاربر می تواند این نتایج را با استفاده از جستجوگر پروژه بررسی کرده و ببینید. کاربر همچنین می تواند نتایج را با استفاده از Wizard و یا ابزار تجزیه و تحلیل آماری به منظور بررسی ریسک راه حل های تصادفی اصلاح نماید.

داده های اولیه و توسعه مدل آبخوان

به منظور توسعه و آماده سازی سامانه های پیچیده جریان به داده های با وضوح بالا نیاز است.

در سال 1988، شرکت DuPont با تکمیل تحقیقات آب های زیرزمینی در دوپانت بل در غرب ویرجینیا، تشخیص داد که به منظور اقناع سازمان حفاظت محیط زیست ایالات متحده در حفاظت از منابع و قانون بهبود (RCRA) به مجوز اقدامات اصلاحی و  تسهیلات ویژه ای نیاز است. تمرکز تحقیقات بر تعدادی از واحد های مدیریت مواد زائد جامد با سابقه مورد استفاده، در منطقه کوهستانی 600 هکتاری در مجاورت یک کارخانه بود. این کار توسط شرکت DuPont و بخش گروه بازسازی شرکت (DCRG) و مشاوران URS انجام شد.

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

تغییرات آب و هوایی و افزایش سطح آب دریا می تواند تغییرات قابل توجهی در تولید رواناب ها و جاری شدن سیل در مناظر سواحل کم ارتفاع را سبب شود. یک چالش عمده برای مقامات دولت های محلی و تصمیم گیرندگان در ترجمه این اثرات جهانی بالقوه تغییرات آب و هوایی در انطباق با استراتژی های عملی و مقرون به صرفه در شهرستان ها و در اصل مقیاس های شهری است. از مدل پردازش MODFLOW برای نمایش زیر مقیاس شبکه هیدرولوژی در مناطق شهری برای کمک به رسیدگی به این مسائل استفاده شده است. رهگیری توأم، آب های سطحی، فروافتادگی ها، و ذخیره منطقه غیر اشباع، ارائه گردیده است. ویژگی های اضافی شامل ساختار، مانع ها، جریان مجاری آب بین سلول های مجاور، مرزهای ایستگاه های مشخص، مرزهای جریان بحرانی، شرایط آب سطحی source/sink، و رواناب دو جهته به پردازشگر MODFLOW در روند یابی آب سطحی اعمال گردیده است. برخی از توانایی های مدل پردازش رواناب شبکه (URO) با یک مسأله مصنوعی با استفاده از چهار کاربری اراضی و پوشش متفاوت سلول بررسی شده است. برای بارش از طوفان فرضی استفاده شد و سلول به سلول عمق آب، سطح آب های زیرزمینی، سرعت نفوذ، و نرخ تغذیه آب های زیرزمینی نشان داده شد. نتایج مشخص کرد مدل پردازش URO با موفقیت توانایی تولید متغیرهای زمانی، آب-محتوا وابسته به نفوذ و نشت را به همراه مدل MODFLOW دارا است.

تغییرات در (a) سطح آب های زیر زمینی و (B) ضخامت لایه اشباع در آبخوان High Plains از زمان پیش از توسعه تا سال 1997.

 ( سازمان زمین شناسی ایالات متحده، 1998)

آبخوان High Plains

High Plains یک منطقه 174 مایل مربعی مسطح زمین است که شامل بخش هایی از کلرادو، کانزاس، نبراسکا، نیومکزیکو، اوکلاهما، داکوتای جنوبی، تگزاس، و وایومینگ است. این منطقه توسط بارش متوسط مشخص و گسسته سازی شده است اما به طور کلی دارای یک نرخ تغذیه طبیعی پایین به سمت سیستم آب زیرزمینی است. رسوبات آبرفتی نامستحکم که تشکیل یک سطح آب سفره به نام آبخوان High Plains را می دهد (شامل حد زیادی از آبخوان اوگالالا) در ساختار منطقه است. آب آبیاری پمپ از آبخوان High Plains را به یکی از مناطق مهم کشاورزی کشور تبدیل کرده است.

نمونه هایی واقعی در زمینه چگونگی پاسخ سیستم آبهای زیر زمینی به تغییر در تقاضای پمپاژ

لانگ آیلند از شمال به مناطقی معروف به همان ناحیه آیلند، از شرق و جنوب به اقیانوس اطلس و از غرب به خلیج نیویورک و رودخانه شرق محدود شده است. لانگ آیلند به چهار شهرستان تقسیم شده است. دو شهرستان های غربی، تحت عنوان Kings و Queens، بخشی از شهر نیویورک است.

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

مثال فرضی از چگونگی پاسخ سیستم آبهای زیر زمینی به تغییرات تقاضای پمپاژ

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

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

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

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