Page tree

Versions Compared

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.
Show If
grouparax
Panel
bgColorpapayawhip
titleARAX

@wikipedia


Consider a differential volume 

LaTeX Math Inline
body\delta V
 with surface area 
LaTeX Math Inline
body\delta \Sigma
 and mass 
LaTeX Math Inline
body\delta M
.

Пусть в этом объеме за время
LaTeX Math Inline
bodydt
изменится масса сплошной среды за счет внешних источников (например, за счет притока флюида из ствола скважины, который вскрыл объем
LaTeX Math Inline
body\delta V
).

Изменение массы
LaTeX Math Inline
bodyd \, \delta M
 сплошной среды в этом объеме за время
LaTeX Math Inline
bodydt
приводит к изменению плотности среды
LaTeX Math Inline
bodyd \rho_V \, \delta V
за то же время и притоку/оттоку флюида через поверхность объема
LaTeX Math Inline
body\delta V
за то же время:

LaTeX Math Block
anchor1
alignmentleft
\frac{d \, \delta M}{dt} = \int_{\delta V} \frac{d \rho_{\delta V}}{dt} dV +  \oint_{\delta \Sigma} {\bf j}_m \, dS

где

LaTeX Math Inline
body{\bf j}_m = \frac{dm}{dt \, dS} {\bf e}_{\bf u} = \rho \, {\bf u}
, где
LaTeX Math Inline
body\rho
– плотность флюида,
LaTeX Math Inline
body\bf u
– вектор скорости потока флюида.

В случае порового коллектора насыщенного флюидом, масса сплошной среды в объеме
LaTeX Math Inline
body\delta V
можно представить как сумму массы скелета пород
LaTeX Math Inline
body \delta m_{rock}
и массы флюида в порах
LaTeX Math Inline
body \delta m
:
LaTeX Math Inline
body\delta M = \delta m_{rock} + \delta m
.

В процессе фильтрации изменением плотности материала пород можно пренебречь по сравнению с изменениями плотности флюида и, следовательно, масса скелета пород в данном объеме не меняется

LaTeX Math Inline
body\frac{d \, \delta m_{rock}}{ dt} =0
т. е.
LaTeX Math Inline
body\frac{d \, \delta M}{ dt} = \frac{d \, \delta m}{ dt}
.

При этом плотность среды

LaTeX Math Inline
body\rho_{\delta V}
выражается через плотность флюида как
LaTeX Math Inline
body\frac{ d\rho_{\delta V} }{dt}= \frac{ d (\rho \phi) }{dt}
.

Тогда закон сохранения массы среды запишется в виде:

LaTeX Math Block
anchor1
alignmentleft
\frac{d \, \delta m}{dt} = \int_{\delta V} \frac{d (\rho \phi)}{dt} dV +  \oint_{\delta \Sigma} \rho {\bf u} \, dS

Выражая поверхностный интеграл через объемный

LaTeX Math Block
anchor1
alignmentleft
\oint_{\delta \Sigma} \rho {\bf u} \, dS = \int_{\delta V} \nabla \cdot (\rho {\bf u}) dV

получим закон сохранения массы в виде:

LaTeX Math Block
anchor1
alignmentleft
\frac{d \, \delta m}{dt} = \int_{\delta V} \bigg( \frac{d (\rho \phi)}{dt} + \nabla \cdot (\rho {\bf u}) \bigg) dV 

Обозначим изменение массы флюида в выделенном объеме пласта за счет внешних источников через скорость

LaTeX Math Inline
bodyq_m
притока массы
LaTeX Math Inline
body\frac{ d \, \delta m}{dt} = \int_{\delta V} q_m dV
и перепишем закон сохранения массы как


LaTeX Math Block
anchor1
alignmentleft
 \frac{d \, \delta m}{dt} =  \int_{\delta V} q_m dV = \int_{\delta V} \bigg( \frac{d (\rho \phi)}{dt} + \nabla \cdot (\rho {\bf u}) \bigg) dV 

и в силу произвольности выбора объема

LaTeX Math Inline
body\delta V
сплошной среды можно переписать закон сохранения массы в дифференциальном виде как

LaTeX Math Block
anchor123
alignmentleft
\frac{\partial ( \rho \phi)}{\partial t} + \nabla \cdot ({\rho \bf u}) = q_m

и называется уравнением непрерывности сплошной среды.