The work described in this thesis is primarily concerned with the construction and study of laboratory scale models for the process of air injection into liquid-saturated grain packing. Experiments, both in two-dimensional (2D) and three-dimensional (3D) setups, were carried out using water-saturated packings of glass beads and/or packings of crashed fused silica glass grains saturated with a glycerin-water solution. High resolution digital images of the invasion patterns were recorded and analyzed. During air injection into a vertically-placed 2D glass bead packing saturated with water, three stages were identified, a tree-like pattern, a fluidized pattern, and a migrating single-channel pattern. The expansion of the tree-like pattern behaves in a diffusion-like manner as the air branches advance upward randomly, and finally reach a more or less constant width. The starting position of the fluidized pattern was quantitatively estimated via balancing the pressure forces between the effective stress due to the weight of the grains and the pressure resistance on the displaced fluid combined with the capillary pressure. Four dynamic regimes were distinguished: regime (i) where the fluidization stops somewhere between the top of the packing and the injection orifice, a transition regime (ii), regime (iii) where the fluidization reaches the injection orifice, and regime (iv) where the deformation of the packing appears as soon as the air is injected. A critical injection rate