کد پایتون مدل سازی مدل جریان گذرا آبخوان آزاد

۲۷ آبان ۹۶ ، ۱۹:۰۳ نويسنده: بهزاد سرهادی

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

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

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

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

دوره 3 (100 روز) - پمپاژ چاه در مرکز مدل با نرخ = 10-، چپ و راست GHB = 10. و 0.

ما با دستورات مدل انتخاب شده از آموزش قبلی شروع خواهیم کرد.

توجه داشته باشید که می توانید این اسکریپت آموزش پایتون را از اینجا دریافت کنید.

آغاز محاسبات

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

import flopy

ایجاد مدل MODFLOW

اندازه مدل، حلقه شبکه و خصوصیات را تعریف کنید.

اطلاعات مدل را اختصاص دهید:

# Model domain and grid definition
Lx = 1000.
Ly = 1000.
ztop = 0.
zbot = -50.
nlay = 1
nrow = 10
ncol = 10
delr = Lx / ncol
delc = Ly / nrow
delv = (ztop - zbot) / nlay
botm = np.linspace(ztop, zbot, nlay + 1)
hk = 1.
vka = 1.
sy = 0.1
ss = 1.e-4
laytyp = 1

# Variables for the BAS package
# Note that changes from the previous tutorial!
ibound = np.ones((nlay, nrow, ncol), dtype=np.int32)
strt = 10. * np.ones((nlay, nrow, ncol), dtype=np.float32)