диусом рс, где потенциал равен Фс. Тогда примем, что на расстоянии рк от центра скважины потенциал равен Фк. Для де бита скважины в плоскости £ можно написать формулу Дюпюи
2nh (Фк — Фс)
(11.103)
In (Рк/Рс)
Перейдем снова к плоскости z. При больших значениях у те чение в полосе — будет параллельным оси у. Для этой
оси из (11.101) имеем
р л; shny/o.
Поэтому, согласно рис. 46, можно положить
n L
Из этого выражения соответственно получим lnpK= nLfo—In 2.
При значительных расстояниях по оси у имеем nz^$>a. Тогда можно положить
In рк =55 nLfo.
При незначительных пу/а
я у я у
асГ
е |
— е |
^ я у |
ягс |
|
2 |
а ~ |
а * |
Следовательно,
1прс= 1п(лгс/ст).
Подставляя приведенные значения 1прк н 1прс в формулу (11.104), получим
2лА/1 (рк— Рс) |
2л£А(рк— Рс) |
2akh (рк — рс) |
(11.104) |
(In Рк— 1П Рс) |
|
|
|
|
|
|
r [ L + - r lai k )
По формуле (11.104) можно определить дебит одной скважины из бесконечной цепочки скважин, расположенных в неограни ченном пласте, при условии, что на некотором, достаточно большом расстоянии L от осп х давление равно рЕ, а в сква жинах малого радиуса гс оно составляет рс.
2. Рассмотрим решение одной из основных задач теории теплопроводности, весьма необходимое при расчетах тепловых методов разработки нефтяных месторождений. Пусть имеем полубесконечныи стержень площадью сечения S, полностью
теплоизолированный от окружающей среды. Начальная
6 Ю. И. Желгпю® |
81 |
niiiinimiNUni |
|
|
температура при |
t = 0 |
во |
всем |
|||||||
|
|
стержне |
была равна |
Т0, а |
|
при |
|||||||
|
|
|
|
|
|
на границе |
стержня |
х = 0 |
|||||
|
|
|
|
|
(рис. |
47) |
она |
стала |
равной |
Ти |
|||
|
|
|
|
|
оставаясь |
при |
t— >-оо |
равной |
Т0. |
||||
|
|
|
|
|
Требуется |
определить |
распреде |
||||||
|
|
|
|
|
ление температуры по |
координа |
|||||||
|
|
|
|
|
те л; в различные моменты вре |
||||||||
|
|
|
|
|
мени t. Будем исходить из урав |
||||||||
|
|
|
|
|
нения |
сохранения энергии, |
|
рас |
|||||
|
|
|
|
|
сматривая |
теплоперенос в |
стер |
||||||
|
|
|
|
|
жне только за счет теплопровод |
||||||||
|
|
|
|
|
ности. Для скорости теплопере- |
||||||||
Puc. |
47. |
Схема |
распростране |
носа |
|
vT за счет |
теплопроводнос |
||||||
ния |
температуры |
за счет |
теп |
ти имеем |
следующее |
уравнение: |
|||||||
лопроводности в полубесконеч- |
дх}т |
■ф - ^ = 0. |
|
(11.105) |
|||||||||
ном стержне: |
стержень пло |
|
|||||||||||
1 — полубесконечный |
дх |
|
|
|
|
|
|
|
|
||||
щадью сечения S; |
2 — распределе |
Здесь |
с — удельная |
теплоемко |
|||||||||
ние |
температуры в |
стержне |
в мо |
||||||||||
мент времени t |
|
|
сть |
вещества |
в |
стержне; |
|
р — |
|||||
|
|
|
|
|
|
||||||||
Скорость переноса |
тепла |
плотность вещества. |
|
мож |
|||||||||
vT за |
счет теплопроводности |
||||||||||||
но определить по формуле закона Фурье |
|
|
|
|
|
|
|||||||
= —Я. |
дт |
|
|
|
|
|
|
|
|
(11.106) |
|||
|
|
дх |
|
|
|
|
|
|
|
|
|
|
|
где Хт — коэффициент теплопроводности. |
|
|
|
|
|
|
|||||||
Подставляя (11.106) в (11.105), получим |
|
|
|
|
|
||||||||
д*Т |
дТ |
|
hr |
|
|
|
|
|
|
(11.107) |
|||
* т дх2 — |
dt ’ |
*T~ |
ср • |
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|||||
Уравнение (11.107) есть уравнение теплопроводности при пря молинейном распространении тепла, а входящий в него коэф фициент хт называется к о э ф ф и ц и е н т о м т е м п е р а т у р о п р о в о д н о с т и . В соответствии с условиями задачи
Т = Т0 при х > 0, t = 0; |
t > 0, х ----->■оо, |
(11.108)
Т = Тг при х = 0, t > 0.
Рассмотрим функцию f(x, t), определяемую следующим об разом:
f(x,t) = {T -T,)l(T1- T |
0). |
(11.109) |
||
Тогда |
начальное и граничное |
условия (11.108) запишутся сле |
||
дующим образом: |
|
|
||
/ = 0 |
при я > 0, |
t = 0; |
t y O |
, x ----►сю, |
/ = 1 |
при х = 0, |
t > 0. |
|
(11.110) |
|
|
|||
82
Функция f(x, t), очевидно, также удовлетворяет уравнению теп лопроводности (11.107), как и Т (х, t), т. е.
ау |
(11.111) |
|
'т дх2 а/ ’ |
||
|
Для получения решения рассматриваемой задачи применим преобразование Лапласа, для чего умножим левую и правую части (11.111) на e~st (s — некоторый параметр) и проинтегри руем эти части в пределах от нуля до бесконечности. В резуль тате получим
|
|
(11. 112) |
Будем считать |
преобразованием Лапласа |
функции f(x, t), |
функцию F(x,s), |
причем |
|
оо |
|
|
F{x,s) = ^f(xJ)e~stdt. |
(11.113) |
|
о |
|
|
Учитывая независимость переменных х и t, функцию f{x,t) можно дифференцировать под знаком интеграла. Из (11.113) имеем, помня, что s — некоторый параметр,
|
оо |
|
|
d*F |
а |
i-St dt. |
(11.114) |
dx2 |
дх2 |
|
|
После интегрирования правой части выражения (11.112) полу
чим |
00 |
|
I |
|
|
|-е -* < й = |
—/ (х, t) е_5/-(- s j" / (х, t) e~st dt = sF {x, s). |
(II. 115) |
Первый член в выражении (11.115) равен нулю, так как при верхнем пределе он равен нулю в результате стремления к нулю экспоненты, а при нижнем пределе вследствие того, что f(x, 0) =0 по условию задачи.
После подстановки (11.115) в (II.112) |
|
|
«T^ - s F |
= 0. |
(11.116) |
Решение |
обыкновенного уравнения (11.116) |
имеет вид |
F = С е |
|
(11.117) |
Для установления постоянной интегрирования С выполним граничное условие (11.110). Однако найдем прежде всего, че-
6’ |
83 |
му равно F (О, s): Из граничного условия (11.110)
/40,5) |
e~st dt = \ e~st dt — — . |
(11.118) |
|
s |
|
Тогда |
|
|
|
X |
|
F (х, s) |
е |
(11.119) |
|
Функцию f(x, t ) по ее изображению F(x, s) найдем по табли цам оригиналов функций и их изображений по Лапласу. Имеем
X
2 Via
f(x, 0 = 1 |
2 |
I |
e~z2dz= 1 —erf |
(11. 120) |
|
l/n |
|||||
|
|
|
Получим, наконец, выражение для скорости переноса тепла на границе х=0. Из приведенного решения с учетом (11.106) на ходим
|
дТ |
-ЯТД7\ |
df |
|
|
т IJC=0 = —кт дх дг=0 |
dx Ц=о |
|
|||
|
*2 |
|
|
|
|
|
4 |
ХтДтг |
&Т1 = Т1— Т0 |
(11.121) |
|
= К Ь Тг У |
[х=0 |
||||
|
|||||
|
|
|
|||
Поток тепла qT через сечение стержня площадью 5 при х = 0
_ |
(11. 122) |
тУ nxTt
3.Рассмотрим цриток жидкости (нефти) с постоянным де битом q к точечному стоку, расположенному в однородном бес конечно простирающемся плоском пласте толщиной h при упругом режиме. Сток находится в центре координат, и тече ние к нему в пласте радиальное. В начальный момент време
ни t = Q, |
пластовое |
давление |
постоянно |
и составляет |
рк. |
|||||
При |
|
из точечного стока отбирается |
из пласта нефть с де |
|||||||
битом q= const, а |
пластовое |
давление |
остается равным |
рк |
||||||
только при г— мх>. Требуется определить |
распределение давле |
|||||||||
ния в пласте в любой момент времени. |
|
|
|
|||||||
Уравнение неразрывности массы фильтрующегося в пласте |
||||||||||
вещества имеет в рассматриваемом случае |
следующий вид: |
|
||||||||
a (pi') |
1 |
ри |
| |
д (тр) |
_ |
л |
|
|
(11.123) |
|
дг |
*" |
г |
"* |
dt |
|
'• |
|
|
||
.84
Учитывая закон Дарси и сжимаемость пласта (сжимаемость пород пласта и насыщающей их жидкости) из (11.123) полу чим уравнение упругого режима в следующем виде:
k |
( д2р ,___1_ |
Эр |
R |
др |
|
1(х |
^ дг2 |
г |
дг J |
" |
dt ’ |
(11.124)
где рс и рж — сжимаемость соответственно пород пласта и на сыщающей пласт жидкости. Остальные обозначения такие же, что и принятые выше в формуле закона Дарси. Введем функ цию f(r, t) следующим образом:
2nkh(pK— р)
(11.125)
q\i
и |
подставим |
ее |
в уравнение |
(11.124). В результате получим |
I |
дг2 ' Г |
дг / |
dt • |
(11.126) |
|
||||
Здесь х — пьезопроводность |
пласта. Поскольку сток точечный |
|||
(г— >-0), то для |
него имеем |
следующее граничное условие: |
||
Следовательно, граничное и начальное условия будут
(г |
= —!. |
f(r. |
0) = 0. |
|
|
(11.127) |
|||||
Известно, что рассматриваемое решение задачи зависит от |
|||||||||||
одной |
переменной \,= rl^yit. |
|
В таких случаях считают, что ре |
||||||||
шение |
автомодельное, |
т. е. |
подобное |
самому |
себе. Поэтому |
||||||
/= /(£ ). Имеем |
|
|
|
|
|
|
|
|
|||
df __ |
г, |
г |
|
df |
__ |
г, |
1 |
d2f __ р„ |
1 |
(11.128) |
|
~dt |
2f \/x t |
1 |
~дг |
~ |
' |
Y x t ’ |
dr2 |
x f |
|||
|
|||||||||||
Подставляя эти значения производных в основное уравне |
|||||||||||
ние (11.126), получим |
|
|
|
|
|
|
|
||||
и’ + - ^ - = 0, |
|
« = /'?. |
|
|
|
|
|
(11.129) |
|||
Из |
(11.127) |
имеем следующие условия: |
|
|
|||||||
f = 0 |
при |
6 ----►сю; |
|
|
|
|
|
(11.130) |
|||
|
|
|
|
|
|
|
|
|
|
||
85