Summary
We present a simple 1-D and one-lithology model of integrated hydrocarbon generation and vertical migration. Although the model is very simple, it captures some of the most basic features of the vertical migration process, like expulsion from source rocks, and vertical migration in carrier beds above and below the source rock. The dynamics in the model are caused by sedimentation, which is done at a constant rate. Owing to sedimentation, the source rock is brought down to the depths where kerogen breakdown starts. This two-phase model is based on fractional flow theory, which makes it possible to avoid a pressure equation. The solution of the saturation equation is characterized by a dimensionless number, the gravity number, and the ‘size’ of the capillary forces. When the gravity number is less than 1, buoyancy is less important, while for gravity numbers greater than 1, buoyancy is dominating. We discuss capillary effects, and estimate a Peclet number for the migration process, which shows that capillary forces are less important in this one-lithology model. In the special case of no compaction, negligible capillarity and no generation of pore space due to the kerogen breakdown, we can approximate a solution to the saturation equation by considering Riemann problems. We find that above the source region, we get at least a shock and possibly a rarefaction, too. Below the source region, we have at least a rarefaction wave, and possibly a shock. We present explicit expressions for the shock position, and the characteristics propagating saturations. These solutions enable us to decide if, and when, the shock will reach the top of the basin. If we know the total oil production, we can read off from the flux function curve an approximation of the saturation propagated upwards. Although these solutions apply to the case of no compaction, negligible capillarity and no generation of pore space, they are useful approximations if any of these restrictions is only ‘weakly violated.
[1]
J. E. Smith,et al.
The dynamics of shale compaction and evolution of pore-fluid pressures
,
1971
.
[2]
D. F. Oltz,et al.
Long-Range Petroleum Migration in the Illinois Basin
,
1991
.
[3]
T. Gallouet,et al.
Hydrocarbons Migration in Basin Modelling: is the Combined Use of Finite Element and Control Volume Possible?
,
1992
.
[4]
J. Norbury,et al.
Two-stage, Two-level Finite Difference Schemes for Non-linear Parabolic Equations
,
1982
.
[5]
Alan K. Burnham,et al.
Comparison of methods for measuring kerogen pyrolysis rates and fitting kinetic parameters
,
1987
.
[6]
K. Aziz,et al.
Petroleum Reservoir Simulation
,
1979
.
[7]
M. Hubbert,et al.
Entrapment of Petroleum Under Hydrodynamic Conditions
,
1953
.
[8]
G. Korvin,et al.
Shale compaction and statistical physics
,
1984
.
[9]
Thomas F. Corbet,et al.
Linear and nonlinear solutions for one‐dimensional compaction flow in sedimentary basins
,
1988
.
[10]
C. Marlé.
Multiphase Flow in Porous Media
,
1981
.
[11]
W. A. England,et al.
The movement and entrapment of petroleum fluids in the subsurface
,
1987,
Journal of the Geological Society.
[12]
R. R. Berg,et al.
Capillary Pressures in Stratigraphic Traps
,
1975
.
[13]
T. M. Quigley,et al.
Principles of Geochemical Prospect Appraisal
,
1988
.
[14]
D. W. Peaceman.
Fundamentals of numerical reservoir simulation
,
1977
.
[15]
P. Ungerer,et al.
Basin Evaluation by Integrated Two-Dimensional Modeling of Heat Transfer, Fluid Flow, Hydrocarbon Generation, and Migration
,
1990
.
[16]
I. Lerche,et al.
The hydrodynamic shape of oil stringers
,
1989
.
[17]
J. Bear.
Dynamics of Fluids in Porous Media
,
1975
.
[18]
L. F. Athy.
Density, Porosity, and Compaction of Sedimentary Rocks
,
1930
.
[19]
K. Terzaghi.
Theoretical Soil Mechanics
,
1943
.