We consider the problem of modeling flow through naturally fractured porous media. In this type of media, various physical phenomena occur on disparate length scales, so it is difficult to properly average their effects. In particular, gravitational forces pose special problems. In this paper we develop a general understanding of how to incorporate gravitational forces into the dual-porosity concept. We accomplish this through the mathematical technique of formal two-scale homogenization. This technique enables us to average the single-porosity, Darcy equations that govern the flow on the finest (fracture thickness) scale. The resulting homogenized equations are of dual-porosity type. We consider three flow situations, the flow of a single component in a single phase, the flow of two fluid components in two completely immiscible phases, and the completely miscible flow of two components.