Numerical models of petroleum migration via buoyancy-driven porosity waves in viscously deformable sediments

The apparent ability of petroleum to migrate rapidly through low permeability sediments in sedimentary basins has led to many questions about the manner of its transport. One possible explanation is suggested from observations of the compaction of viscously deformable porous media. These systems have been found in some cases to give rise to regions of locally elevated liquid fraction, in the form of fluid‐filled porosity waves that can ascend at rates much greater than that of the background flow. Previous research on the phenomenon has focused on its implications for magma transport, but recognition of the fact that the compaction of viscous porous media can take place in sedimentary basins has suggested the possibility that porosity waves could similarly be important for hydrocarbon transport. The purpose of the present study was to test this hypothesis by quantifying the transport that would occur as the result of porosity waves initiated during the conversion of kerogen to petroleum. A one‐dimensional numerical model was constructed solving equations for the mechanics of viscous compaction and for the kinetics of reactions describing the formation of petroleum from kerogen. The results showed that porosity waves would develop readily in viscously deformable regions of sedimentary basins, but would not necessarily provide enhanced transport over that of the background flow regime. In order for the waves to achieve this enhanced transport, they must develop high amplitudes, i.e. high porosities relative to the background porous medium. To achieve the high wave amplitudes, the background porosity must be very low in absolute terms. In addition, high kerogen contents are needed in the source layer, and the source layer needs to be buried rapidly to a high temperature region of the oil window. Considerable uncertainty exists as to the value of the matrix shear viscosity of sediments in basins. However, the wave volumetric transport capacity was not found to be significantly altered as a result of variations in the value of this parameter. The physical form of the waves was strongly altered by the matrix shear viscosity, with higher values leading to lower amplitudes and generation frequency, but higher wavelengths. Thus the waves become less recognizable physically at higher values of the matrix shear viscosity. As the waves ascend to higher stratigraphic levels, where the porosity is higher, they gradually lose their physical definition and become absorbed into the background.