We simulate confined droplets in microchannels by depth-averaged equations solved by a boundary element method. The retarding effect due to film formation is absent in the depth-averaged approach and added by a nonlinear boundary condition. Although deformable the terminal velocity of the droplets streamed in these channels changes only little from the undeformed state. Using the singularity method we develop an analytical model for the droplet velocity at low Reynolds and low capillary number.