<jupyter_start><jupyter_text>Сравним полученные значения со стационарными формулами:<jupyter_code>import math as m
pi = 3.14
Rк = lx/2 # м, радиус контура питания
dP = (p_init-int(prod_bhp[0]))/10*10**6 # Па, перепад давления
k = permx/10**15 # м2, проницаемость (в числовой модели изотропный пласт)
h = 10 # м, толщина пласта
L = 200 # м, длина горизонтального участка в пласте
mu = 0.0005 # Па*с, вязкость
rc = 0.073 # м, радиус скважины
Bo=1.2 # коэффициент объемного расширения
a = lx # м, сторона модели

# Вариант №1
# Батлер
Q4 = (2*3.14*k*h*dP)/(mu*(m.log1p(Rк/rc)-0.75))*86400/Bo
print('Приток по 1 варианту: ' + str(round(Q4)) + ' м3/сут')

# Вариант №2
# Батлер - приток на псевдоустановившемся режиме к протяженной горизонтальной скважины 147 с.
kh = k
kv = permz/10**15
Rm = 0.5*rc*(1+m.sqrt(kh/kv))
hm = h*m.sqrt(kh/kv)
Q5 = (6*kh*L*h*dP)/(mu*(a/2+3*hm/3.14*m.log1p(hm/2/3.14/Rm)))*86400/Bo
print('Приток по 2 варианту: ' + str(round(Q5)) + ' м3/сут')

# Вариант №3
# Батлер - неуст. режим в центре прямоуг зоны с учетов вертикального схождения
S=0
A = lx*ly
prow = A/(L**2) 
print('Условие выбора формулы: ' + str(prow))
if prow < 35:
    P = m.log1p(1+prow)/12
else:
    P = (m.log1p(prow)+0.1515)/(4*3.14)
P_ = hm/(2*3.14*L)*(m.log1p(hm/(2*3.14*Rm*m.sin(3.14/2)))+S)
Q1 = kh*h*dP/mu/(P_+P)*86400
print(h, dP, mu, P)
print('Приток по 3 варианту: ' + str(round(Q1)) + ' м3/сут')<jupyter_output>Приток по 1 варианту: 2557 м3/сут
Приток по 2 варианту: 4841 м3/сут
Условие выбора формулы: 56.25
10 12000000.0 0.0005 0.3343095256715246
Приток по 3 варианту: 9800 м3/сут
{'fraction_non_alphanumeric': 0.17405063291139242, 'fraction_numerical': 0.11139240506329114, 'mean_word_length': 4.305369127516778, 'pattern_counts': {'":': 0, '<': 5, '<?xml version=': 0, '>': 4, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 3, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
{'chain_length': '1', 'dir': 'jupyter-structured-clean-dedup', 'id': '448633', 'license': 'Unlicense', 'n_tokens_mistral': 921, 'n_tokens_neox': 923, 'n_words': 171, 'path': 'generator/horiz_well_res/Сравнение с аналитическим решением.ipynb', 'repo_name': 'unifloc/opm_examples', 'size': '1580'}