книги из ГПНТБ / Телков, А. П. Подземная гидрогазодинамика
.pdfворяла последнему граничному условию в IX (2). В конечном сче те получается решение для распределения потенциала в пласте, свободное от указанных ограничений:
Ф(о . р', 1, h) = Ф0 |
- 4 |
f (р,р',Е,А); |
|
IX (8 ) |
||||
|
|
|
П0 |
|
|
|
|
|
|
|
|
|
у |
z |
i b |
|
1Х(9) |
а П0 ■ |
;/ - |
л г |
’ |
|
Ао |
й = |
тh0- |
|
где F — некоторая функция, |
выражаемая |
в рядах |
и интегралах |
|||||
через функции |
Бесселя |
[16]. |
|
|
имеем Ф == |
Фс. Тогда |
||
При г — гс |
(второе |
граничное условие) |
||||||
из IX (8 ) следует формула |
IX (5), |
где |
|
|
|
|||
|
Е |
= |
■F(o, h) |
|
|
IX (1 0 ) |
||
3.Точное решение задачи о потенциале точечного стока
(источника) в однородно-анизотропном круговом |
пласте. Приток |
||
к несовершенной |
скважине |
|
|
Как известно, в области, содержащей источники или стоки, |
|||
потенциал удовлетворяет уравнению Пуассона [16] |
|
||
Д<р = ф (*, у, |
г): |
|
IХ (1 1 ) |
Здесь А— оператор Лапласа; |
<р— потенциал; |
<1> (х, у, |
г)— |
плотность стоков как функция координат. |
|
г = |
|
Рассмотрим приток к точечному стоку с координатами |
|||
=0 , 2 = т), расположенному в круговом осесимметричном одно
родно-анизотропном |
пласте |
(рис. 43). |
|
|
|
|
||
В этом случае уравнение IX(1) будет иметь вид: |
|
|
||||||
<Э3ф , |
1 |
д ср |
1 д-ю |
Ь (г) „ / |
. |
2 |
kr |
|
+ -L £» + |
й |
- |
|
*’ = |
f |
,Х(12> |
||
дг2 |
г |
дг |
|
|||||
Р и с. 43. Приток к точечному стоку и линии стоков
91
где |
q —-мощность |
точечного |
стока; |
б (х) — функция Дирака. |
|||||
|
Будем считать |
кровлю ( z= |
0) и подошву (z = |
ft0) |
непроницае |
||||
мыми, т. е. |
|
|
|
|
|
|
|
|
|
|
—■— 0 |
при |
z = 0 |
и |
z — h.Q |
I Х(13) |
|||
На |
внешнем контуре |
(г = |
R0) |
примем |
|
|
|
||
|
|
<р = |
сро |
при |
г — R0 |
|
1Х(14) |
||
|
Математическое |
решение задачи |
IX (12)—IX |
(14) |
выполнено |
||||
Сткляниным Ю. И. с помощью метода интегральных преобразова ний [16].
Точное решение для потенциала точечного стока имеет вид:
~ ch [т * (I_ А)1 ch т - ^ 7о 0 ?о. к)
ср (р, 6, Л) |
= ср0 ------q^- |
- |
2 |
L Р_____________ J |
Р_____________________ |
IX (15) |
||||||
|
|
14 sh |Xj/p /? (н-i) |
|
|||||||||
£ < ft |
|
|
|
*Ro <=1 |
|
|
|
|||||
cp(p, £, ft) |
=cp0 |
Д 7 1 |
|
chtjT |
(1~ £)] ch7 ' /t/o(^ |
^ |
IX(16) |
|||||
l > h |
|
|
|
- R 0 i=i |
|
l^i sh|Aj/p/? (|лг) |
|
|
||||
Здесь |
[1/ — положительный |
корень |
уравнения |
J 0 (p-i) = |
0; |
|||||||
J0 (x) — функция |
Бесселя |
первого рода нулевого порядка; |
||||||||||
J x (х) — функция |
Бесселя |
первого рода первого порядка; |
||||||||||
безразмерные |
координаты и |
параметры: |
|
|
|
|||||||
|
Р |
|
%h0 |
|
|
|
|
|
|
|
|
IХ( 17) |
|
|
|
|
|
|
|
|
|
|
|
||
Найдем потенциал для линии стоков (рис. 43). В формулах IX (15) |
||||||||||||
и IX (16) |
примем |
|
|
|
|
|
|
|
|
|||
q = q (ft) = const |
|
при |
0 |
< |
ft < |
1; |
ft = -£- |
IX (18) |
||||
|
|
|
|
|
|
|
|
|
|
|
«0 |
|
Очевидно, |
что |
потенциал для линии стоков будет |
|
|||||||||
|
|
|
|
|
ф = ftol ? (Р,l,h)dh, |
|
1Х(19) |
|||||
|
|
|
|
|
|
|
о |
|
|
|
|
|
где ср |
определяется по формулам |
IX (15) и IX (16). Интегрируя |
||||||||||
IX (19), |
после некоторых преобразований получим окончательно: |
|||||||||||
|
|
|
|
|
In R0 + |
«о |
c h ^ e s h !^ -(i_ /l) /0( ^ |
|) |
||||
Ф — Ф0 |
= |
± |
|
2 2 |
— i — s-J!-----5----------- |
|
||||||
I < h |
|
|
2nhho |
|
|
i= l |
|
И? sh (Aj/p /? (fif) |
IX (20) |
|||
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
92
Рассматривая приток к элементарному кольцу, интегрируя полученное уравнение вдоль кольца — стоков и усредняя потен циал вдоль стенки цилиндра — стоков, т. е. вдоль вскрытой части пласта, получим формулу для дебита:
^ |
2тс А0 (Ф0 — Фс) |
_ 2к h0(Фб — Фс) |
IX (2 2 ) |
|
|
|
|
|
|
где |
|
|
|
IX (23) |
|
’ |
h 1=1 |
^ 3 sh tL /?(N) |
|
|
|
|||
93
|
|
о |
Ro |
|
|
|
_ | |
|
|
|
q |
r |
1 |
|
К*, г |
^ |
1—< |
||
1 |
||||
|
|
* |
Кg
1
=с:
С......
7
Р и с . 45. Схема притока к несовершенной скважине в двухслойном однородноанизотропном пласте
Функция IX (23) рассчитана на ЭВМ, графическое изображе' ние которой представлено на рис. 44.
Задача о притоке газа специально не ставилась, однако приведенное здесь решение может быть использовано и для прито ка газа к несовершенной скважине при установившемся изотерми ческом процессе по линейному закону фильтрации.
Переходя от потенциала к давлению в формуле IX(22) и при нимая вместо давления р функцию Лейбензона Р, а вместо объем ного расхода Q— весовой расход G = уamQ. а затем делая обрат ный переход, получим
Q = |
. р* ~ р* |
IX (22') |
ЦРат |
Во |
|
4. Потенциал несовершенной скважины в двухслойном однородно-анизотропном круговом пласте
Решение данной задачи (рис. 45) выполнено Ю. И. Сткляниным. Решения для распределения потенциала (Фх— потенциал для верхнего пласта, Ф2— потенциал для нижнего пласта) выра жаются в сложной интегральной форме. Для скважины радиуса гс потенциал на скважине определится интегралом
ь
ф с = I q (tj) <р (гс, г, т]) d г}, |
IX (24) |
О |
|
94
где т), q и ®— ордината, мощность |
и потенциал точечного стока |
соответственно. |
можно добиться постоянства |
Соответствующим подбором q (Д |
|
потенциала на забое вскрытой части |
пласта, т. е. Фс = constant |
при 0 < Z < Ь. |
|
Удобнее за потенциал на стенке скважины принять среднее значение его, вычисленное по формуле:
IX (25)
о
где Oj (g, р, К) определяется сложным выражением [16], которое из-за громоздкости здесь не приводится.
5. Потенциал точечного стока (источника), горизонтальной дрены и несовершенной галереи в полосообразном однородно
анизотропном пласте
Вопросы притока пластовых жидкостей к горизонтальным скважинам, дренам и трещинам рассматривались рядом авторов в различной постановке. Эти вопросы становятся все более ак туальными, поскольку современное развитие техники бурения наклонных и горизонтальных скважин делает вполне реальным их широкое практическое разрешение в ближайшее время.
|
Рассмотрим |
приток к точечному |
стоку |
с координатами (tj, |
||||
/ х ) , |
расположенному несимметрично |
в |
полосообразном горизон |
|||||
тальном |
однородно-анизотропном |
пласте с |
подошвенной водой |
|||||
(рис. 46). |
исходное |
принимается уравнение |
|
|||||
|
За |
|
||||||
|
|
|
|
|
|
|
|
IX (26) |
где |
б (z— tj) |
и |
б (х— 1г) — функции |
Дирака. |
||||
|
Введем потенциал |
|
|
|
|
|||
|
|
|
|
|
|
|
|
IX (27) |
и характеристику |
анизотропности |
пласта |
|
|||||
|
|
|
|
|
|
|
|
IX (28) |
|
Тогда |
уравнение IX (26) примет |
вид |
|
|
|||
|
|
|
|
|
|
|
|
IX (29) |
Кровля и подошва считаются непроницаемыми, т. е.
— = 0 |
при |
z = 0, |
z = h0 |
IX (30) |
дг |
|
|
|
|
95
Р и с. 46. Схема притока к точечному стоку, горизонтальной дрене, несовер шенной щели в полосообразном пласте
На контуре питания для простоты принимается
9 = 0 |
при |
х = 0, |
х = / |
IX(31) |
Уравнение IX (29) с граничными условиями IX (30) и IX (31) можно решить методом интегральных преобразований, применяя последовательно косинус - и синус - преобразование Фурье с конечными пределами и формулы обращения. Такое решение дано в работе [16]:
«, ch m л 2L (й0 — г) ch /я то) —
?! = |
— Щ-Ъ ----------------------—т----------- — |
IX (32) |
||
г>т] |
11 m~i т sh т ъ Л„ ~ I sin т я fi sin т к |
— |
|
|
|
0 / ' |
I |
I |
|
ооch т п ~ (/г0— T | ) c h m i z i
92 = - f f X |
-----------Ц ------------ |
p - i |
-----р |
IX (33) |
г<7: |
m=1 т sh т тс n0 — /sin т и is i n |
т я — |
|
|
Формулы IX (32) и IX (33) дают распределение потенциала, вызванного точечными стоком и источником, в элементе анизо тропного полосообразного пласта. Их можно использовать при экспериментировании на щелевом лотке и для горизонтальной дрены. Можно получить потенциал и для несовершенной щели в том же пласте. При этом щель будем рассматривать как линию стоков, расположенных вдоль прямой х = от г = 0 до г = Ь, при постоянной мощности q (рис. 46). Тогда потенциал линии сто ков определится интегралом
Ф = [ 9 (х >г,п) d у} |
IX (34) |
о |
|
96
Подставляя IX (32) и IX (33) в IX (34), интегрируя последнее уравнение и делая некоторые преобразования, получим
|
|
Ф ,= |
I ср^х, |
z , rt) d |
7] = |
|
|
||
|
|
zyb |
о |
|
|
|
|
|
|
|
|
|
* тих ,, |
|
. , |
т их , |
|
||
|
= _ |
ch -у— (h0— г) sh —j~ |
b |
|
|||||
|
_ , 0 , m г:у. , , |
mu e* |
. m и x |
|
|||||
|
|
1m2 s h ----- |
A0/s |
n -----1 s i n --------- |
|
||||
|
|
|
|
/ |
0 |
|
; |
l |
|
|
|
, ffln |
|
, min |
,, |
,. . |
m - l. . |
mu. x |
|
|
|
ch------ z sh |
----- |
(h,— b) ■sin----- i- s i n ------- |
|||||
Ф» = |
-^ -2 |
|
l |
|
|
______ l |
/ _ |
||
|
rrfi sh |
—/г„ |
|
|
|||||
z < i |
n |
m = l |
|
|
|
||||
IX (35)
+
( — |
2l? " 7 F “ x |
ПР Н |
|
0 < |
X < |
l 1 |
|
+ |
I |
(l — x) |
при |
/, < |
x < |
IX(36) |
|
— 2qj- |
/ |
||||||
Е с л и b = h0, |
то |
из формулы |
IX (36) |
получается выражение |
|||
для потенциала |
совершенней |
щели в полосообразней залежи: |
|||||
( — |
|
|
ПРИ 0 < х < |
||||
Ф ,= |
|
|
|
|
|
|
1X (37) |
V— 2q *2 (/ — х) |
|
при [х < Х<(/ |
|||||
Если за горизонтальною дрену принять горизонтальную сква жину диаметром d, произвольно расположенную в пласте, то потенциал такой дрены определится интегралом
ь
Ф = ) С£г (х, г, т,) d У\ |
1Х(38), |
е—d |
|
Ряды IX (32) и IX (33) сходятся медленно, однако после улуч шения сходимости [16] для них можно получить приближенные выражения с достаточной степенью точности, более удобные для вычислений.
6. Анализ распределения потенциала в однородно-анизотропном: круговом пласте
Обычно аналитические решения о распределении потенциала скоростей фильтрации в пласте получаются в виде бесконечных рядов или интегралов, что связано с трудоемкими численными расчетами. В настоящее время, благодаря использованию элек тронной вычислительной техники, стало возможным проводить сложные расчеты и анализы в широких диапазонах параметров.
Впервые детальный анализ распределения потенциала на поверхности забоя несовершенной по степени вскрытия скважины
4 Заказ 612 |
97 |
|
в однородном пласте с осевой симметрией, а также в пласте по всей мощности относительно несовершенной скважины, дан М. Мас-
кетом [3].
Представляет теоретический и практический интерес (напри мер, при решении некоторых статических и динамических задач теории конусообразования) показать более наглядно характер распределения потенциала в однородно-анизотропном ограничен ном пласте, вскрытом несовершенной скважиной (рис. 43). Исполь зуем решение IX (21), которое в безразмерных параметрах прини мает следующий вид:
|
|
|
= |
^ |
( 1 , |
Ро, |
Л , |
|
|
|
IX (39) |
|
|
|
|
ch 1--- - (ij sh J L iHh(\4 -FA |
|
||||||
F ( £, Ро, h, |
р - |
Е |
|
Po__________ P |
' |
R p J |
IX (40) |
||||
|
|
A sh X i- i U ih) |
|
||||||||
|
|
|
i—1 |
|
|
|
|||||
Расчет функции IX (40) производился |
на ЭВМ для параметра |
||||||||||
^ = 10 и относительного вскрытия |
h = |
0,3 |
при |
фиксированных |
|||||||
значениях |
r/R0■ Результаты |
представлены |
на рис. 47. Степень |
||||||||
отклонения |
от |
плоскорадиального |
потока |
для |
каждого значе |
||||||
ния r/Rо можно |
показать из |
соотношения |
|
|
|
||||||
|
|
5% = |
- ^ ° ’5 |
/ |
Ь1 |
• 10096 |
|
I X(41) |
|||
|
|
|
Г£=0,° |
|
|
|
|
|
|
||
Р и с . 47. Распределение потенциала скоростей фильтрации в круговом пласте (к определению зоны пространственного движения)
9 8
Такое отклонение составляет: 17,2; 11,2; 1,18 и 0,25% для отношения радиусов 0,01; 0,05; 0,10 и 0,20 соответственно.
Таким образом, |
для зоны r0 ~ rJR 0 > 0,10 практически |
можно считать поток |
плоскорадиальным, поскольку на расстоянии |
г0 —0,10 отклонение составляет величину порядка одного процента.
На рис. 47 показаны линии равных потенциалсв, где линия г0 — = 0,10 представляет ссбой почти прямую, параллельную оси абс цисс.
Чтобы нагляднее представить себе размеры зон пространст
венного и |
плоскорадиального течения, |
рассмотрим |
один |
пример. |
||||||||
Пусть |
R 0 = |
200 м, h0 = |
20 м и х ^ |
1 |
(пласт сднорсднс-изотроп- |
|||||||
ный), |
что соответствует |
параметру с0 = |
10. Для |
|
вскрытия |
h = |
||||||
; ^ 0,3, |
принимая га |
0,10 как |
границу |
между |
указанными |
зо |
||||||
нами, |
находим т-- |
0,10 |
• R 0 |
0,10 • |
200 |
20 |
м. |
Как |
видим, |
|||
зона пространственного притока оказалась равной одной мощности пласта. С увеличением относительного вскрытия, как это явствует из физических соображений, эта зона будет уменьшаться, прибли жаясь к скважине, и в предельном случае вскрытия пласта (сква жина совершенная) во всей области дренажа будет иметь место стро го плоскорадиальный приток. Что касается анизотропных пластов, то, очевидно, зона пространственного движения будет меньше по
сравнению с однородно-изотропными пластами и |
в общем случае |
г0 = / (р, Щ |
|
Проведенный анализ подтверждает выводы |
М. Маскета и |
И. А. Чарного о размерах зоны пространственного движения при притоке к несовершенной скважине по степени вскрытия пласта. Маскет показывает, что в однородно-изотропном пласте ПрИ от'
носительном вскрытии |
h = 0,5 |
плоскорадиальность |
потока на |
рушается лишь при г < |
2h0. И. |
А. Чарный предлагает принимать |
|
зону пространственного потока |
радиуса г = (1 -ь 1,5) |
h0. |
|
Произведенная количественная оценка позволяет принять за критерий, характеризующий приток к несовершенной скважине, параметр р. При 10 и ft > 0,3 (что выполняется в большинст ве практических случаев) зона пространственного притока с высо
кой |
степенью точности может быть принята равной r0 = h0. При: |
р > |
Ю этот радиус будет несколько больше и для практических: |
расчетов может быть принят равным удвоенной мощности пласта. Следует заметить, что выбор зоны пространственного движе ния при притоке к несовершенным скважинам имеет важное зна чение во многих задачах подземной гидродинамики. В дальнейшем будет широко использован прием решения многих задач, основан
ный на схеме условного разделения притока на две зоны.
4*
99
7. Приток неньютоновсксй жидкости к несоверненной скважине
Нефти некоторых месторождений СССР (Башкирии, Азер' байджана, Узбекистана) содержат большое количество пара" фино-смолоасфальтеновых компонентов. Наличие таких компонен" тов может обусловливать проявление структурно-механических свойств нефти (динамическое напряжение сдвига и структурная вязкость) при ее фильтрации в пористой среде. Исследования по казали, что чем больше нефть содержит смолоасфальтеновых ком понентов и чем ниже температура, тем большее значение приобре тают ее структурно-механические свойства.
Известно, что для движения обычных (ньютоновских) жидко стей в пористой среде, когда абсолютная вязкость жидкости ос тается постоянной, справедлив закон Дарси. Движение вязко пластичных (неньютоноЕСких) жидкостей характеризуется дина мическим напряжением сдвига и структурной вязкостью. В этом случае фильтрация жидкости описывается так называемым обобщен ным законом Дарси или законом фильтрации с преденльым гра диентом давления сдвига, достаточно полно исследован м в ра ботах А. X. Мирзаджанзаде и др., Султанова Б. И., В. М. Ентова. Обобщенный закон Дарси можно записать по аналогии, вводя коэффициент структурной вязкости Дх'0) как функцию динамиче ского напряжения сдвига (т0) [17]
и |
— grad |
— р + р grad |
k |
I * |
|
Л Со) |
IX (42) |
||||
|
|
Ti Со) |
|
! |
|
При некоторых условиях, |
когда величина второго слагаемого |
||||
в уравнении |
IX (42) |
составляет незначительную долю* от всей |
|||
суммы и ей можно пренебречь без существенного ущерба для
практических расчетов, обобщенное уравнение |
Дарси IX (42) |
|||
может быть |
записано приближенно |
в виде |
|
|
|
и х — grad —— р = |
— grad Ф |
IX (43) |
|
|
|
лСо) |
|
|
Отсюда следует, что |
если известно решение для |
распределения |
||
потенциала |
скорости |
фильтрации |
Ф ньютоновской жидкости |
|
при притоке к совершенной или несовершенной скважине, то оно может быть использовано для приближенных расчетов и для слу
чая фильтрации вязко-пластичной |
жидкости |
путем введения |
||
т) (т0) вместо коэффициента абсолютной |
вязкости |
р. ньютоновской |
||
* Оценка величины grad |
—- — .может |
быть |
произведена по формуле с [17] |
|
|
Wo) |
|
|
|
grad |
Г[1+ к/и (г)] |
|
х = А т |
|
т,(хо) |
|
ц |
|
|
100
