\phi \cdot \partial_\tau \Psi -   
 \nabla \cdot \left( k \cdot   \vec \nabla \Psi   \right) 
  = 0

-\frac{k}{\mu} \int_{\Sigma} \, \frac{\partial p}{ d n} \, d \sigma = q(t)

where

	reservoir pressure		time
	effective porosity		position vector
	total compressibility		position vector of the -th source
	formation permeability to a given fluid		gradient operator
	dynamic viscosity of a given fluid		Pseudo-Pressure
	fluid compressibility factor		Pseudo-Time
			sandface flowrates to/from the well-reservoir contact

See also