Вопросы методологии и новых технологий разработки месторождений природного газа. Часть III (Сборник научных трудов), страница 41

Рассуждая аналогично предыдущему, выполним предельный переход при Ах -» 0, после чего левая часть уравнения (8) примет вид [11]

\pv         \pvvxdF\pwdL,                                                                             (16)

Сл F                           (Ж рL

Используя предельные переходы для произвольной функции координат

90


lim

Д.Г--И) Ax f.


S6OK

Ax

выполним предельные переходы в выражениях, стоящих в правой части уравнения (8).

Для первого выражения имеем

где F- вектор массовой силы, отнесенный к единице массы флюида; dF - элемент поверхности поперечного сечения горизонтального ство­ла.

При выполнении предельного перехода для второго члена уч­тем, что на поверхности F[ имеем п = i, а на поверхности F п = - /; 101 да


1 lm — Г pnds - lim —

&x->o Ax i              Лх-»п Ах


\\ndF- [pidF



д

lim

О


J pwt/tV = —J p/ £//• -+- j pndL


При вычислении предельного перехода для третьего члена бу­дем иметь в виду, что на поверхности F] т„ ~ T\v'a на поверхности F Тп--тх'- тогда


lim — f r^dS = lim


Ax


\rxdF-[zxdF


lim



С учетом полученных предельных переходов закон сохранения количества движения потока флюида в необсаженном горизон­тальном стволе скважины примет вид следующего интегродиффе-ренциального уравнения [11]:


<?

— j pvdF + — jpvvxdF - \pvwdL - JpA/F---------------- / p/ dF -

г1                                                Г                                  i,                               /'                                    я

+\zndL.

L


(17)

91


Уравнение (17) представляет собой векторную форму закона со­хранения количества движения и в гаком виде мало пригодно для практического применения. Более удобно записать его в проекциях на гидравлическую ось х, для чего необходимо вычислить проекции на ось нх" всех членов этого уравнения.

Для левой части уравнения (17) получаем

{a\p~vdF) =


> 

ApvvxdF\  = — \pvxdF +


f


pvx


v    + v \


dF.


(последний  член  учитывает  влияние  возможной  кривизны  ги­дравлической оси).

В дальнейшем будем считать касательную составляющую век­тора скорости на стенке ствола равной нулю, т.е. v - vn • п; учитывая также, что (п)х = Cos(nJ) = Sin(J j), получим

! | р \hvdL    = j | pw vnndL j = - J pw vnSinu , i VlL.

Для правой части уравнения (17) находим UpFdF\  =\pFxdF

VJ x     f

(FK - проекции вектора массовых сил на ось "х");


(


д


х


F

pndL\ =-


F


"XX


.dx.,


lrv


(at


dF.


Здесь было использовано разложение вектора гг по ортам ',/Д, где ^ог- нормальное напряжение, действующее в плоскосги x=const в направлении оси "х";

92


тху  - касательное напряжение, действующее в плоскости x=const в направлении оси "у";

Txz - касательное напряжение, действующее в плоскости x=const в направлении оси "z".

Разложим вектор напряжения х„ по ортам элемента боковой поверхности

(?п)х = rnnCo.s{nJ)+ TntCos{lj)+ TnmCos{mJ\

Здесь Coslmj} = 0, т.к. орт т  лежит в плоскости, параллельной плоскости F, которая, в свою очередь, перпендикулярна орту i С учетом предыдущего, имеем

| JzndL   = }[- TnmSin(I, J)+ rnl( :os(IJ)]dL.

Используя вышеприведенные соотношения, запишем уравнение сохранения количества движения в проекции на гидравлическую ось горизонтального ствола в следующем виде:

дх..

(18)


F


/■■


"X2 I

Перейдем к выводу уравнения сохранения энергии для потока флюида в горизонтальном стволе. Аналогично предыдущему, выпол­ним предельный переход при Дх->0 для левой и правой частей уравнения (10). Левая часть в силу граничных условий (12) при­водится к виду


д


д


\wdL


Предельные переходы для членов правой части уравнения (10)

выполняются аналогично предыдущему; в результате получим

93