Understanding the role of thin films in porous media is vital to elucidate wettability at the pore level. The type and thickness of films coating pore walls determine reservoir wettability and whether or not reservoir rock can be altered from its initial state of wettability. Pore shape, especially pore wall curvature, is important in determining wetting‐film thicknesses. Yet, pore shape and physics of thin wetting films are generally neglected in flow models in porous rocks. Thin‐film forces incorporated into a collection of star‐shaped capillary tubes model describe the geological development of mixed wettability in reservoir rock. Here, mixed wettability refers to continuous and distinct oil and water‐wetting surfaces coexisting in the porous medium. This model emphasizes the remarkable role of thin films.New pore‐level fluid configurations arise that are quite unexpected. For example, efficient water displacement of oil (low residual oil saturation) characteristic of mixed‐wettability porous media is ascribed to interconnected oil lenses or rivulets which bridge the walls adjacent to pore corners. Predicted residual oil saturations are approximately 35% less in mixed‐wet rock than in completely water‐wet rock. Calculated capillary pressure curves mimic those of mixed‐wet porous media in the primary drainage of water, imbibition of water, and secondary drainage modes. Amott‐Harvey indices range from −0.18 to 0.36 also in good agreement with experimental values (Morrow et al., 1986; Jadhunandan and Morrow, 1991).

Foam in porous media exhibits an unusually high apparent viscosity, making it useful in many industrial processes. The rheology of foam, however, is complex and not well understood. Previous pore-level models of foam are based primarily on studies of bubble flow in circular capillaries. A circular capillary, however, lacks the corners that characterize the geometry of the pores. We study the pressure–velocity relation of bubble flow in polygonal capillaries. A long bubble in a polygonal capillary acts as a leaky piston. The ‘piston’ is reluctant to move because of a large drag exerted by the capillary sidewalls. The liquid in the capillary therefore bypasses the bubble through the leaky corners at a speed an order higher than that of the bubble. Consequently, the pressure work is dissipated predominantly by the motion of the fluid and not by the motion of the bubble. This is opposite to the conclusion based on bubble flow in circular capillaries. The discovery of this new flow regime reconciles two groups of contradictory foam-flow experiments.Part 1 of this work studies the fluid films deposited on capillary walls in the limit Ca → 0 (Ca ≡ μU/σ, where μ is the fluid viscosity, U the bubble velocity, and σ the surface tension). Part 2 (Wong et al. 1995) uses the film profile at the back end to calculate the drag of the bubble. Since the bubble length is arbitrary, the film profile is determined here as a general function of the dimensionless downstream distance x. For 1 [Lt ] x [Lt ] Ca−1, the film profile is frozen with a thickness of order Ca2/3 at the centre and order Ca at the sides. For x ∼ Ca−1, surface tension rearranges the film at the centre into a parabolic shape while the film at the sides thins to order Ca4/3. For x [Gt ] Ca−1, the film is still parabolic, but the height decreases as film fluid leaks through the side constrictions. For x ∼ Ca−5/3, the height of the parabola is order Ca2/3. Finally, for x [Gt ] Ca−5/3, the height decreases as Ca1/4x−1/4.

This work determines the pressure–velocity relation of bubble flow in polygonal capillaries. The liquid pressure drop needed to drive a long bubble at a given velocity U is solved by an integral method. In this method, the pressure drop is shown to balance the drag of the bubble, which is determined by the films at the two ends of the bubble. Using the liquid-film results of Part 1 (Wong, Radke & Morris 1995), we find that the drag scales as Ca2/3 in the limit Ca → 0 (Ca μU/σ, where μ is the liquid viscosity and σ the surface tension). Thus, the pressure drop also scales as Ca2/3. The proportionality constant for six different polygonal capillaries is roughly the same and is about a third that for the circular capillary.The liquid in a polygonal capillary flows by pushing the bubble (plug flow) and by bypassing the bubble through corner channels (corner flow). The resistance to the plug flow comes mainly from the drag of the bubble. Thus, the plug flow obeys the nonlinear pressure–velocity relation of the bubble. Corner flow, however, is chiefly unidirectional because the bubble is long. The ratio of plug to corner flow varies with liquid flow rate Q (made dimensionless by σa2/μ, where a is the radius of the largest inscribed sphere). The two flows are equal at a critical flow rate Qc, whose value depends strongly on capillary geometry and bubble length. For the six polygonal capillaries studied, Qc [Lt ] 10−6. For Qc [Lt ] Q [Lt ] 1, the plug flow dominates, and the gradient in liquid pressure varies with Q2/3. For Q [Lt ] Qc, the corner flow dominates, and the pressure gradient varies linearly with Q. A transition at such low flow rates is unexpected and partly explains the complex rheology of foam flow in porous media.