MINISTERE DES MINISTERIE VAN AFFAIRES ECONOMIQUES ECONOMISCHE ZAKEN

ADMINISTRATION DE LA BESTUUR QUALITE ET DE LA SECURITE KWALITEIT EN VEILIGHEID

GEOLOGICAL SURVEY OF PROFESSIONAL PAPER 1999/1- N.289

BURIAL HISTORY AND COALIFICATION MODELLING OF WESTPHALIAN STRATA IN THE EASTERN CAMPINE BASIN (NORTHERN BELGIUM)

by S. Helsen and V. Langenaeker

T

Q

-50001+----r----t----+----+---T------T----1 -350 -300 -250 -200 -150 -100 -50 0 Time (myr) MINISTERE DES MINISTERIE VAN AFFAIRES ECONOMIQUES ECONOMISCHEZAKEN

ADMINISTRATION DE LA BESTUUR QUALITE ET DE LA SECURITE KW ALITEIT EN VEILIGHEID SERVICE GEOLOGIQUE DE BELGIQUE BELGISCHE GEOLOGISCHE DIENST

GEOLOGICAL SURVEY OF BELGIUM* PROFESSIONAL PAPER 199911, N. 289

BURIAL HISTORY AND COALIFICATION MODELLING OF WESTPHALIAN STRATA IN THE EASTERN CAMPINE BASIN (NORTHERN BELGIUM)

by

S. Helsen 1 and V. Langenaeker2

1999

1 Historische Geologie, Instituut voor Aardwetenschappen, Redingenstr. 16, B-3000 Leuven present address: ECOREM, Wayenborgstraat 21, B-2800 Mechelen 2 Geologica n.v., Tervuursesteenweg 200, B-3060 Bertem present address: D.H.V., B-3600

Comite editorial : L. Dejonghe, P. Laga, R. Paepe Redactieraad: L. Dejonghe, P. Laga, R. Paepe Service Geologique de Belgique Belgische Geologische Dienst Rue Jenner, 13 - 1000 Bruxelles Jennerstraat 13, 1000 Brussel

* «The Geological Survey of Belgium cannot be held responsible for the accuracy of the contents, the opinions given and the statements made in the articles published in this series~ the responsability resting with the authors >). Toute !'information economique en claire et net sur Aile economische informatie rechtstreeks op http : 1/mineco.fgov.be http : 1/mineco.fgov.be

Editeur responsable: Roland PAEPE Verantwoordelijke uitgever: Roland PAEPE rue Jenner 13 Jennerstraat 13 1000 Bruxelles 1000 Brussel Depot legal: D 1999/0880/3 Wettelijk depot: D 1 999/0880/3 Table of contents

ABSTRACT ...... 3

KEYWORDS ...... 3

1. AIM OF THE STUDY ...... 4

2. GEOGRAPHIC AND GEOLOGIC SETTING ...... 4

3. MODELLING ...... 4 3.1. SOFTWARE ...... 4 3.2. INPUT DATA ...... 5 3.2.1. Stratigraphic data ...... 5 3.2.2. Decompaction data ...... 5 3.2.3. Porosity and thermal conductivity data ...... 5 3.2.4. Heat flow models ...... 6 3.2.5. Surface temperature model ...... 6

4. RESULTS ...... 6 4.1. BURIAL HISTORY ...... 6 4.2. COALIFICATION ...... 7 4.2.1. Carboniferous ...... 7 4.2.1.1. Wells KB69, KB 79, KB8l, KB86 and KBI85 ...... 7 4.2.1.2. Wells KB l 15 and KB !63 ...... 8 4.2. I .3. Wells K.Bl74, KBI86 and KBLommel (coalification minimum north ofZolder) ...... 8 4.2. I .4. Wells KBI68, KBI72, K20I and KBI98 ...... 8 4.2.2. Post-Carboniferous ...... 8 4.3. HYDROCARBON POTENTIAL...... 9

5. CONCLUSIONS ...... 9

6. ACKNOWLEDGEMENTS ...... 10

7. REFERENCES ...... 10

TABLES ...... 13

FIGURES ...... 23 2 ABSTRACT

One dimensional basin modelling techniques are applied in the eastern Campine Basin in Belgium to reconstruct the burial history and thermal maturation of -bearing strata. Separated by two phases of regional uplift and subse­ quent erosion, subsidence occurred during the Carboniferous, the Late Permian through Jurassic and the Late Creta­ ceous to Cenozoic. A model with a first phase of important burial during the Westphalian/Stephanian and a second during the Triassic through Middle Jurassic is proposed. For the western and southern wells it appears that maximum coalification was attained immediately after the first stage of burial. For the other wells, especially for those in the northeastern part of the coalfield, maximum thermal maturation may have evolved from the second phase of heating in the Late Jurassic ahd Early Cretaceous. A best fit between calculated and measured coalification data is obtained by a combination of high heat flows and important packages of Late Carboniferous and Early to Middle Jurassic strata. 2 Heat flows of 102-lll mW/m , assumed for Westphalian and Stephanian times, were succeeded by moderate values of 80-60 m W/m2 in the eastern and 70-60 mW/m 2 in the western part of the Cam pine coalfield. Estimated thicknesses for now eroded Westphalian D/Stephanian strata are in the range of 1250-1600 m, while thicknesses of Rhaetian through Dogger rocks are varying between 1250-1980 m. Pre-Zechstein erosion of Carboniferous rocks was more important in the south than in the north with maximum values of 2480 m. As a result of Late Jurassic to Early Cretaceous erosion related to the Cimmerian uplift, 2320-2740 m of Carboniferous through Jurassic sediments were removed in the central and southeastern part of the study area. Coalification data indicate thermal potential for gas, which was mainly generated during Carboniferous burial in the south, and during both Carboniferous and Mesozoic burial in the central and northern parts.

KEY WORDS: Burial and thermal history, coalification, Carboniferous, Westphalian, Campine Basin (Belgium)

3 1. AIM OF THE STUDY

A modelling study was carried out for the Belgian Geological Survey to examine the burial history and the evolution of the coalification of post-Namurian coal-bearing strata in the Campine coal basin. The aim of this study is to reveal the timing of maximum coalification and the related hydrocarbon generation. In the same area a model study for subsurface methane production in a coal bed methane test well was already evaluated by Fermont et a!. (1994 ). For this study I 3 well locations were selected in the coalfield in the eastern Cam pine Basin, where sufficient stratigraphic and coalification data are available from exploration wells (Fig. 1). For better geographic coverage a 14th, pseudo-well was chosen in an area where data are scarce. Labeled KBLommel, this well is situated 0.5 km west of the centre of the namesake town. The stratigraphic succession is based, with sufficient confidence, on extrapolation from nearby wells. Software used in this modelling was developped by one of us (V.L.), based on the models of Sweeney & Burnham (1990) and Sclater & Christie (I 980).

2. GEOGRAPHIC AND GEOLOGIC SETTING

The Cam pine Basin covers a considerable area of the Belgian provinces of and Limburg, as well as adjacent areas in the southern . The stratigraphy of the gently northward dipping Devonian and Carboniferous rocks is only known from deep exploration wells and collieries in the mining district. Upper Namurian and Westphalian strata consist of a cyclic deposition of shales, siltstones and sandstone with coal. To date, Stephanian rocks have not been recognised in the area. However, this does not exclude their presence in earlier times. In the Campine Basin deformation of the Palaeozoic strata during the Asturian phase of the Variscan orogenic cycle or during the Early Permian Saalian phase resulted in a series of tilted blocks (Bouckaert & Dusar, 1987). A considerable pile of sediment is assumed to have been eroded before the deposition of Late Permian (Zechstein) and Triassic fluviatile and near­ shore sediments (Tys, 1980). Early to Middle Jurassic (Lias-Dogger) rocks were deposited in marine environments, created by a rapid subsidence of the Central Netherlands Basin and the Rur Valley Graben (Geluk et al. I 994; Zijerveld et al. I 992). The northwest-southeast striking Jurassic shoreline probably crossed northeastern Belgium (Ziegler, 1990). A second period of regional erosion occurred as a result of the Late Jurassic (Maim) to Early Cretaceous Cimmerian uplift, during which important packages of Mesozoic and Late Palaeozoic strata were removed. Renewed sedimentation started during a Late Cretaceous highstand. However, inversion movements along major faults (Rossa, 1987) resulted in a temporarily period of non-deposition and erosion in the Rur Valley Graben. At present, Quaternary fluvial and aeolian deposits and Tertiary to uppermost Cretaceous marine deposits cover the older substratum.

3. MODELLING

3.1. SOFTWARE

The modelling software was developed by V.L. while at the Institute of Earth Sciences of the K. U.Leuven. Programs are written in TurboPascal, and integrate equations for several procedures during the modelling of burial and coalification. Decompaction, formed during burial, is based on the exponential relationship between porosity and depth. The equa­ tions of Sclater & Christie ( 1980) are used for detrital sedimentary rocks. The decompaction algorithm is described by Sclater & Christie (1980). The relationship between porosity and thermal conductivity of rocks is described by Lewis & Rose (1970). For sand and sandstone a constant thermal conductivity of 8 meal/ems EC was assumed. For other lithologies a combination of matrix thermal conductivity and porosity was assumed as follows: 5.5 meal/ems EC for shale, 7.75 meal/ems EC for silt, 6 meal/em sEC for chalk and I .5 meal/em sEC for peat-coal series. Based on the thermal conductivity and the heat flow, burial tempertures are calculated. Coalification, expressed as VRm%, is related to the temperature by a model described in Sweeney & Burnham (1990). This model is based on a number of parallel first order kinetical reactions which describe the thermal alteration of kerogen in source rocks.

4 The software has been calibrated on published examples of input data and subsequently calculated maturity values. The present model is less elaborated than other existing commercial software used by oil companies and some geo­ logical surveys (e.g., GENEX, BasinMod), but has sufficient possibilities regarding the present study. As the model takes into account decompaction, porosity, variations in thermal conductivity, surface temperature and heat flow the thermal model is based on the state-of-the-art thermal maturation theories by Sweeney & Burnham (1990):

3.2. INPUT DATA

3.2.1. Stratigraphic data

Tables 3-16 show for each of the studied sites the lithology, thickness and stratigraphy of the observed post-Namurian succession, as well as the inferred depositional and erosional phases. For this purpose, rocks have been grouped in a number of litho-types, i.e., clay-shale (1), sand-sandstone (2), chalk (3), and silt-siltstone (4). Thicknesses preceded by >-= and >zero= codes for the lithology refer to eroded packages. Geochronologic data are from Harland et al. ( 1990) or Lippolt et al. (1984) for the Westphalian. Stratigraphic terms are explained in Table 2. Thicknesses of Westphalian A-C, Permian and Triassic strata are from the various well descriptions, summarised in Langenaeker (in prep.), or have been obtained through extrapolation from nearby boreholes (Table I). Thicknesses of eroded sequences are obtained through modelling processes. Thicknesses for the Westphalian A series range between 900 m in the eastern mining district and 1050 m in the western mining district. Westphalian B strata show an increasing trend in thickness from the west to the east, i.e., from 520 m in KB 186 to 685 m in KB 115. Westphalian C rocks are approximately 700 m thick in the central part but measure up to 1000 m in the southeastern corner (KB81). A number of deep wells and seismic data indicate the presence of Westphalian D rocks in central Limburg. Seismic sections show that a thickness of at least 500 m is locally present without all being traversed by boreholes. Permian strata attain 40 m maximum, and are considered to have been deposited in the northeastern part of the study area only. Buntsandstein sediments are assumed to have covered the entire eastern Campine-Brabant areas but were probably thicking towards the northeast. Cumulative thicknesses indicate values up to 580 m. In KB 198 a thickness of 700 m is assumed. Muschelkalk and Keuper deposits are each known to be up to 86 m thick at KB99, 5 km northeast of KB 168. In our simulation, these strata vary between 50 m in the southern part and !00 m in KB I 98. Although their extension throughout the Rur Valley Graben is evident from seismic investigations, Rhaetian and Liassic (Hettangian to Pliensbachian) strata are drilled at KB99 and at KB 198 only, where they are at least 480 m thick (Demyttenaere & Laga, 1988). The total preserved thickness from seismic evidence may approach 800 min the graben. Stratigraphic information on the Cretaceous and Tertiary cover is from the isopach maps of Legrand (1968), Demyttenaere ( 1988) and Demyttenaere & Lag a ( 1988). Maximum thicknesses of the Cretaceous succession (Santonian through Maastrichtian) are recorded in the southeast, i.e., 224 m in KB 168. Due to inversion tectonics in the Rur Valley Graben, Cretaceous deposits in KB 198 are reduced to 27 m of Maastrichtian rocks. Cenozoic strata thicken from the south to the north, ranging between 292 min KB86 and I 257 m in KB 198.

3.2.2. Decompaction data

The decompaction of rocks is based on the litho-type chosen, and on the exponential relationships between porosity and depth (Sclater & Christie, 1980). Compaction or decompaction is only excecuted for the first burial phase after sedimentation. During uplift and renewed burial variations in porosity are dominated by diagenetic reactions, rather than by mechanical reduction of porosity. To a lesser extend, this is also the case during the first burial phase. As a result, the calculation of porosity, and thus of decompaction, is in all cases only approximative. By applying decompaction corrections the burial history of strata can be calculated more precisely.

3.2.3. Porosity and thermal conductivity data

Porosity data, as calculated during the decompaction process, are translated into thermal conductivity data using the equations of Lewis & Rose ( 1970) and values for matrix thermal conductivity (section 3.1.).

5 3.2.4. Heat flow models

Actual heat flow in Northeast-Belgium approximates 60 mW/m2 (Hurtig et al. 1992). However, during extension of the Cam pine Basin in Carboniferous times, stretching of the crust may have been accompanied by higher heat flows (Leeder, 1987). In this study an average heat flow of 102 mW/m2 is considered between 315-300 myr, with exception of the southerly locations of KB69, KB81 and KB86, where higher values of 103-111 mW/m2 have been taken into account. All these figures fall well within the range of heat flow values in actual active (syn-rift) back-arc basins (Allen & Allen, 1990). After a period of basin extension heat flow probably decreased during the Permian. At that time, heat flow is supposed to have been higher within and nearby the Rur Valley Graben than in the western part of the study area. This may be due to the initiation of the rifting in the graben. As a result, heat flow may have dropped to 80 mW/m2 at 270 myr and 74 mW/m2 at 260 myr in the eastern part of the study area, and to a value of70 mW/m2 at 270 myr in the western part. In the early Tertiary, at 50 myr, heat flow values of respectively 70 and 65 mW/m2 are considered in the eastern, respectively western part, before reaching the present value (Figs 4-17). Heat flow simulations in the Cam pine Basin by Van Keer et al. (in prep.) are different from the above values. Beside 2 2 an average figure of63 mW/m , high values of84 mW/m for Westphalian D and Stephanian (KB172 and KB 174), 74 mW/m2 during Dogger (KB172) and 96 mW/m2 for Oligocene to Pliocene (KB172 and KB174) \\!ere considered. In the German Ruhr basin a Carboniferous heat flow of83 mW/m2 is supposed to be succeeded by a value of 61 mW/m2 by Littke et al. (1994).

3.2.5. Surface temperature model

Using palaeogeographic maps of western Europe, according to Dercourt et al. (1993) and Ziegler (1989) and to con·esponding actual temperatures for these latitudes, palaeo-surface temperatures for the Cam pine Basin were estimated for the geological record between 315-0 myr. As a result, values of 25EC, refering to tropical climate conditions, are maintained until 300 myr. Due to the progressive shifting of Gondwana towards higher latitudes, the surface temperature is assumed to have reached its maximum at 27EC between 280-270 myr, after which it dropped gradually into the lower twenties. Until Middle Cretaceous times (100 myr) the surface temperature is supposed to have been 20EC or more. In Early Tertiary times it has been assumed to drop from 19EC (at 60 myr) to 16EC (at 35 myr). An average of 1OEC is considered from 9 to 7 myr, after which the surface temperature drops to a minimum of IEC during the Pleistocene at I myr, before reaching the actual 9EC (Fig. 2).

4. RESULTS

4.1. BURIAL HISTORY

Assumed thicknesses of now mostly eroded Westphalian D/Stephanian deposits used in this modelling are in the range of 1250-1570 m (Figs 3-17). These figures are much larger than the 600 m suggested by Bless et al. ( 1977) for the eastern part of the basin. However, Van Keer et al. (in prep.) suggested cumulative thicknesses for the Westphalian D and Stephanian strata of 1900 m, respectively 2000 m for the KB 172 and KB 174 wells. Estimates for the eroded Westphalian series as a result from pre-Zechstein unloading show considerable differences within the Cam pine Basin. From our modelling, it appears that the greatest thicknesses were stripped off in the south, i.e., near the sites of KB 185, KB79, KB86 and KB8l, where more than 2300 m were removed by erosion. This confirms earlier theses, e.g. Grosjean ( 1936). Towards the north, eroded packages are considered to be smaller, averaging around 1200-1300 min the central part of the study area. This trend can be explained by the more important Late Palaeozoic uplift in the southern area, as compared to the north. Considerable stratigraphic differences in Pre-Zechstein erosion levels between relatively close areas may be explained by block faulting and tilting during the Carboniferous (Tys, I 980). This is further corroborated by up to I 000 m of differential removal of Westphalian strata on both flanks of the Donderslag Fault, deviding the Campine coalbasin in an eastern and western block (Dusar & Langenaeker, I 992). In the Rur Valley Graben (KB 198) an eroded thickness of 800 m is assumed. Some of these figures are in good agreement with the results of basin modelling in the German Ruhr Basin. According to Littke et al. ( 1994 ), 2200-3500 m of Carboniferous overburden in this basin was eroded during the Stephanian and Rotliegend.

6 According to Ziegler ( 1990), sedimentation in the southern part of the study area occurred until the Early Dogger (Aalenian), but continued until the Late Dogger (Bathonian) in the central part (KB 174, KB 172), and until the Early Maim (Oxfordian) in the very north of the study area (KB 198). This resulted in huge packages of Late Triassic and Jurassic sediments, varying between 1250 m (KB86) and 1980 m thickness (KB 198) as calculated during the modelling (see also Fig. 3). These figures are in the same range as the 1800 m of Jurassic rocks considered for the KB 172 well by Van Keer et al. (in prep.), but they do not agree with their assumption of 400 m of Jurassic strata (Lias only) at the KB 174 site. Van Keer et al. (in prep.) assigned the important variation in stratigraphical thickness to differences in geological history of the various tectonic bocks. Alternatively, Dusar et al. (1987) suggested that the high coalification data in the KB 201-KB 172 area may be related to a basic intrusion, dated post-Carboniferous to pre-Tertiary, and associated with a gravimetric anomaly (Petrofina, 1963). Similar anomalies and coalification highs nearby Krefeld and Erkelenz, Germany, are probably produced by Permo-Carboniferous, respectively Permian intrusions (Buntebarth et al. 1982; Drozdzewski & Wrede, 1994, among others). Possibly, the Erkelenz anomaly may be an east-bound continuation of the anomaly in Limburg, actually separated by the Rur Valley Graben. In our model, however, high heat flow up to 120 m W/m2 and more, related to a possible intrusion, could not produce the actual high thermal maturation and coalification gradient observed in KB 172. Numerous simulation runs, with a timing for the thermal event shifting between the latest Carboniferous and the Late Cretaceous, combined with varying sediment loading, did not result in any satisfying match between calculated and measured coalification data. Post-Cimmerian erosion reaches maximum values of 2320-2740 min the central and southeastern parts of the studied area. Towards the south and southwest, where the Triassic and Jurassic cover was more reduced, eroded thicknesses are in the range of 2000-2100 m. In the subsiding Rur Valley Graben thick Triassic and Jurassic sequences were less affected by erosion. This can be documented with data from KB 198, where the lower part of the Liassic and probably the complete Triassic record are preserved, and probably >only= some 1340 m of Jurassic rocks have been eroded. The great pile of Westphalian through Jurassic sediments considered to be eroded in the Cam pine Basin is in agreement with the overburden of at least 3000 m on the Ordovician basement of the southern Brabant Massif, estimated to be removed no earlier than the Middle Jurassic (Van den haute & Vercoutere, 1989; Vercoutere & Van den haute, 1993). The oldest preserved Cretaceous rocks are of Santonian age. Together with Campanian and Maastrichtian strata they account for 250 m maximum. Although deposition in the Tertiary was far from continuous (major hiatuses occur in the Eocene and Neogene), erosion during this period was likely less important. The thickest sequences occur in the Rur Valley Graben. Quarternary strata account only for some tens of m maximum. The burial history for some Westphalian stratigraphic markers, as well as for the base of the Cretaceous and Tertiary of the studied boreholes, are depicted in Figs 4a-17a. By using these geohistory plots, the selected heat flow models (Figs 4b-l7b) and the surface temperature model (Fig. 2) the evolution of the burial temperature is calculated (Figs. 4c-17c). The latter is needed to calculate the thermal maturation of organic matter.

4.2. COALIFICATION

4.2.1. Carboniferous

Figs 18-31 show the evolution in thermal maturation of Westphalian strata in all of the 14 studied wells. Coalification data are from the basal strata of the lower Westphalian A (Sarnsbank Marine Band), upper Westphalian A (Wasserfall Marine Band), lower Westphalian B (Quaregnon Marine Band), upper Westphalian B (Eisden Marine Band), Westphalian C (Mam·age Marine Band) and Westphalian D units (Neeroeteren Sandstone). For each of the studied wells, actual (maximum) measured and calculated coalification values for the preserved geological record, including Westphalian rocks, are shown in Tables 2-15, based on Langenaeker (1992). Minor differences between measured and calculated Carboniferous VRm% values may be explained by variations in porosity and thermal conductivity of rocks. In general, however, these differences fall well within the standard deviation of the measured values (Langenaeker, in prep.).

4.2.1.1. Wells KB69, KB 79, KB81, KB86 and KB185

In the southern part of the study area significant thermal maturation was the result of Late Carboniferous burial. VRm% values probably increased after maximum burial at 300 myr, until a maximum in coalification was reached

7 some 290 myr ago. With values of more than 2.4 VRm% for lower Westphalian A strata, thermal maturation was likely most important in KB69 and KB81. Apparently, renewed sedimentary burial resulted in only little, if any, additional thermal maturation for the lower Westphalian. The latter is less than 0.003 VRm% for any of the Westphalian strata in KB86, KB79, KB69 and KB81, but approximates 0.008 VRm% for the base of the Westphalian B (Quaregnon M.B.) in KB 185 (Figs I 8-22; Tables 3-7).

4.2.1.2. Wells KB115 and KB163

Like in the southernmost wells, thermal maturation of Westphalian rocks in the KB I I 5 and KB I 63 wells is mainly due to some 15 myr of Late Carboniferous burial. Apparently, the resulting coalification rapidly increased until 295 myr. However, maximum coalification was only attained after a second period of heating by sedimentary burial during the Permian through Middle Jurassic. Situated at approximately 150 myr, this is actually some 20 myr after the time of maximum burial. This additional thermal maturation can be documented by figures for basal Westphalian A rocks (Sarnsbank M.B.), which are in the range of 0.18 VRm% in KB 163, and 0.22 VRm% in KB 115, resulting in a final coalification for the base of the lower Westphalian A of 1.902 VRm% (KB163) and of 2.088 VRm% (KB l 15). For the Quaregnon Marine Band the additional coalification is even higher, i.e., 0.290 VRm% in KB 163 and 0.320 VRm% in KB 115 (Figs 23-24; Tables 8-9).

4.2.1.3. Wells KB174, KB186 and KBLommel (coalification minimum north ofZolder)

In the western part of the study area Late Carboniferous sedimentary burial resulted in important coalification before 295 myr. Additional thermal maturation of Westphalian rocks, produced by heating during Permian and Mesozoic overburden, attained its maximum during Late Jurassic (Maim) times, and is relatively small. Probably, maximum coalification took place at I 40 myr for the KB I 74 site, shifting to younger ages of I 35 myr and 130 myr for the more northerly KB I 86 and KBLommel wells. The overall modest additional coalification is small in KB I 74, with values of less than 0.050 VRm% for the Sarnsbank and Quaregnon Marine Bands, but it gradually increases to the north with values of 0.220 VRm% for Quaregnon in KB I 86 and 0.670 VRm% for the same stratigraphic level at KBLommel. This trend can be illustrated with the actual VRm values for the Sarnsbank and Quaregnon Marine Bands, which vary from 1.653% and 0.872% in KB 174, to 1.852% and 1.088% in KB I 86 and 2.255% and 1.447% in KBLommel (Figs 25-27; Tables 10- I 2).

4.2.1.4. Wells KB168, KB172, K201 and KB198

In the northeastern part of the study area, the evolution of thermal maturation in time is markedly different from that of the other sites. Being the result of Carboniferous overburden, the coalification attained at 295 myr is of the same range as in the south (KB1 15 and KB163) and the west (KBI74, KBI86 and KBLommel). In wells KB172, KB169 and KB201 Permian and Mesozoic, mainly Jurassic sedimentation evolved in an important increase of coalification, culminating between 145-135 myr. The additional coalification can be documented by the increase in VRm of approximately 1.5 %for lower Westphalian A strata (Sarnsbank M.B.). In KB 198 maximum coalification occurred in the Early Cretaceous at 115 myr, some 29 myr after the time of maximum burial. Here, differences in coalification are as large as 1.7 VRm% for the base of the lower Westphalian A, although smaller differences of0.9 VRm% are noticed for the basal strata of the Westphalian D Neeroeteren Sandstone. Resulting VRm values are high: for the Sarnsbank Marine Band they are 2.516% (KB168), 2.648% (KB172), 2.889% (KB201) and up to 3.567% (KB198). Apparently, in none of the I4 studied wells heating from Cretaceous and Cenozoic burial has affected the already attained coalification level, even not in KB 198, where these strata are more than 1280 m thick (Figs 28-31; Tables 13-16).

4.2.2. Post-Carboniferous

Permian VRm% values are scarce in the Campine Basin. Among the studied wells, only KB 172 and KB20 I yielded vitrinite reflectance values in Zechstein sediments. In both wells three populations of vitrinite reflectance data are determined. The lowest reflectance values, assigned to an indigenous population, are of 0.44% in KB 172 and 0.47%

8 in KB20 I. The other, possibly reworked populations have values of 0.84% and 1.39% in KB 172 and of 0. 76% and 0.85% in KB20 1 (Langenaeker, in prep.; Somers, unpubl. documents, Liege). Being substantially lower, the reflectance data of the indigenous populations do not agree with the values of 0.872% and 0.961%, calculated for the KB 172, respectively KB201 well. However, the small number of measured particles (6 in KBl 72 and 9 in KB20 I) indicates that these measured data are less reliable. A better correspondence is observed between calculated VRm data and reflectance values of vitrinites from the reworked populations. Furthermore, the oxidised nature of the Permian and Triassic rocks possibly hampers the exact identification of the type of maceral, which may be an additional explanation for the discrepancy between measured and modelled values.lt is clear that additional measurements on more numerous vitrinite populations may contribute to a better knowledge of the thermal maturation of Permian deposits in the Campine Basin (Petrotina, 1993). In addition, due to the relatively small number of analysed grains and their low degree of thermal metamorphism (0.46-0.30% ), coalification measurements from Liassic, Upper Cretaceous and Paleocene strata in KB 172, KB20 I and KB 198 (Langenaeker, in prep.; Petrofina, 1993) are less reliable than those from Upper Carboniferous deposits (Tables 14-16 ).

4.3. HYDROCARBON POTENTIAL

Regarding the hydrocarbon potential, the Cam pine coal basin needs to be differentiated in zones with only one impor­ tant phase of hydrocarbon generation at the end of the Carboniferous and zones with an important secondary hydrocarbon generation. In the former areas possible hydrocarbon plays were destroyed by subsequent uplift and erosion. Coalbed methane recovery testing in well KB206, approximately 6 km southeast of KB 174, showed that in this area are undersaturated (Fermont et al. 1994; Wenselaers et al. 1996). In the latter areas possible hydrocarbon plays may exist within Buntsandstein reservoirs.

5. CONCLUSIONS

Based on the above explained calculations the following conclusions may be drawn for the burial history of post­ Namurian strata in the eastern Campine Basin.

I. Important sedimentary burial and thermal maturation took place during the latest Westphalian (WtD) or Stephanian. Presumed cumulative thicknesses for these strata range between 1570 and 1250 m. The maximum cumulative thickness of the Westphalian preserved in the Campine coal basin approximates 3 km. For the well sites in the southern part of the study area maximum coalification was attained shortly after this phase of burial.

2. Asturian and/or Saalian deformation and subsequent erosion resulted in important unloading. In the northeastern part of the coalfield at least 800 m were eroded. In the south, where uplift was more important, some 2300-2400 m were removed. 3. A second period of important subsidence and sedimentation which ended in Middle Jurassic (Dogger) times resulted in a maximum burial and additional coalification for most of the sites in the central and northern parts of the area. Permian through Upper Triassic (Keuper) deposits are 250m thick in the southwestern part and 900 m thick in the northeast. Thicknesses for Rhaetian through Jurassic rocks vary between 1250 min the southern part and 1980 m in the Rur graben. 4. The study area shows thermal potential for gas generation from Westphalian rocks. Possibly in the central, but certainly in the northeastern areas hydrocarbon gasses were probably generated from the Carboniferous source rocks during both periods of burial. 5. Post-Cimmeric erosion was most important in the central and southeastern part of the area with values up to 2740 m. Likely, the smallest erosion, probably some 1340 m, occurred in the Rur Valley Graben. 6. Late Cretaceous through Cenozoic burial has not affected the coalification resulting from earlier periods of heating.

9 6. ACKNOWLEDGEMENTS

Dr. M. Dusar (Geological Survey of Belgium, Brussels) is thanked for his suggestions and constructive comments. Financial s.upport by the Geological Survey of Belgium is greatly acknowledged (project NAT/95-2.2.a). Dr. W. Fermont (NITG, Heerlen) made clear that discussion on heat flow models is far from closed. However, our pragmatic approach may provide an impetus for a more thorough international investigation of the crucial problem of explaining the actual coalification values in the Northwestern European Carboniferous.

7. REFERENCES

Allen, P.A. & Allen, J.R. 1990. Basin Analysis. Principles and applications. Blackwell, Oxford, London, 451 pp. Bless, M.J.M., Bouckaert, J., Calver, M.A., Graulich, J.M. & Pap roth E. 1977. Paleogeography of Upper Westphalian deposits in NW Europe with reference to the Westphalian C north of the mobile Variscan belt. Meded. Rijks Geol. Dienst, N.S. 28 (5), 101-127. Bouckaert, J. & Dusar, M. 1987. Arguments geophysiques pour une tectonique cassante en Campine (Belgique), active au Pah~ozoique superieur et reactivee de puis le Jurassique superieur. Ann. Soc. Geol. Nord 151, 201-208. Buntebarth, G., Michel, W. & Teichmiiller, R. 1982. Das permokarbonische Intrusiv von Krefeld und seine Einwirkung auf die Karbon-Kohlen am linken Niederrhein. Fortschr. Geol. Rheinld. u. Westf. 30, 31-45. Demyttenaere, R. 1988. De post-Paleozo'ische geologische geschiedenis van Noord-Belgie. Unpub. Ph.D. thesis K.U.Leuven, 175 pp. Demyttenaere, R. & Laga, P. 1988. Breuken- en isohypsenkaarten van het Belgische gedeelte van de Roerdal Slenk. Bel g. Geol. Surv. Prof. Paper 234, 20 pp. Dercourt, J. Ricou, L.E. & Vrielynck, B. (eds) 1993. Atlas Tethys Palaeoenvironmental Maps. Gauthier-Villars, Paris, 307 pp., 14 maps, I pl. Drozdzewski, G. & Wrede, V. 1994. Faltung und Bruchtektonik. Analyse der Tektonik im Subvariscikum. Fortschr. Geol. Rhld. u. Westf. 38, 159. Dusar, M., Bless, M.J.M., Borremans, G., Burger, K., De Loose, J., Fairon-Demaret, M., Felder, P.J., Gullentops, F., Lie Sun Fan, Muchez, Ph., Paproth, E., Pierart, P., Rossa, H.G., Smolderen, A., Somers, Y., Steurbaut, E., Streel, M., Viaene, W., Witte, H. & Wouters, L. 1987. De steenkoolverkenningsboring Gruitrode-Ophovenderheide (Boring 172 van het Kempens Bekken), kaartblad Opoeteren, 63E224. Belg. Geol. Surv. Prof. Paper 230, 235 pp. Dusar, M. & Langenaeker, V.l992. De oostrand van het Massief van Brabant, met beschrijving van de geologische verkenningsboring te Martenslinde. Belg. Geol. Surv. Prof. Paper 255, 22 pp. Fermont, W., Jegers, L. & Dusar, M. 1994. A model study of coalbed methane generation in well Peer (KB206) Belgium. Rijks Geol. Dienst Intern. Paper. Geluk, M.C., Duin, E.J.T., Dusar, M., Rijkers, R.H.B., Van den Berg, M.W. & van Rooijen, P. 1994. Stratigraphy and tectonics of the Roer Valley Graben. Geol. Mijnbouw 73, 129-141. Grosjean, A. 1936. Les traits essentiels du terrain houiller de Ia Campine. Ann. Mines de Belgique 37, 263-273. Hurtig, E., Cermak, V., Haenel, R. & Zui, V. (eds) 1992. Geothermal Atlas of Europe. Intern. Assoc. Seismology and Physics of the Earth=s Interior, Intern. Heat Flow Comm., Central Inst. for Physics of the Earth. GeoForschungsZentrum Potsdam, pub!. I, 156 pp. Harland, W.B., Armstrong, R.L., Cox, A.V., Craig, L.E., Smith, A.G. & Smith, D.G. 1990. A Geologic Time Scale, 1989. Cambridge Univ. Press, Cambridge, 263 pp. Langenaeker, V. 1992. Coalification maps for the Westphalian of the Campine Coal Basin. Belg. Geol. Surv. Prof. Paper 256, 34 pp. Langenaeker, V. 1998. The Campine basin. Stratigraphy, structural geology, coalification and hydrocarbon potential for the Devonian to Jurassic. Ph.D. thesis K.U.Leuven, 213 pp. Leeder, M.R. 1987. Tectonic and paleographic models for Lower Carboniferous Europe. In: J. Miller, A.E. Adams & V.P. Wright (eds) European Dinantian environments. Wiley, New York 1-20. Legrand, R. 1968. Le Massif du Brabant. Mem. explic. cartes geol. et min. Belgique 9, 148 pp., 5 pl. Lewis, C.R. & Rose, S.C. 1970. A theory relating high temperatures and overpressures. J. Petr. Techno!. 22, ll-16. Lippolt, H.J., Hess, J.C. & Burger, K. 1984. Isotopische Alter von pyroklastischen Sanidinen aus Kaolin­ Kohlentonsteinen als Korrelationsmarken fiir das mitteleuropaische Oberkarbon. Fortschr. Geol. Rheinld. u. Westf. 32,119-150.

I 0 Littke, R., BUker, C., Li.ickge, A., Sachsenhofer, R.F.& Welte, D.H. 1994. A new evaluation of palaeo-heat flows and eroded thicknesses for the Carboniferous Ruhr basin, western Germany. J. Coal Geol. 26, 155-183. Petrofina, I963. Reinterpretation de Ia campagne sismique 1953 a I 956 et complements 1963 (Heibaart-Bree-Kessenich) pour Ia compte de Ia S.C.R.E.M. Unpubl. document. Petrotina, 1993. Campine basin: thermal maturity assessment of Permian through Tertiary sediments. Spot horizon data compendium. Fin a Research, Exploration and Production, Geolab Report EGGIEB I 34. Pillement, Ph. 1982. La houillification des formations Carboniferes en Belgique. Sa place dans l=histoire sedimentaire et structurale varisque. Unpubl. Ph.D. thesis Universite de Droit, d=Economie et des Sciences d=Aix-Marseille, 220 pp. Rossa, H.G. 1987. Upper Carboniferous and Tertiary inversion tectonics of the Western part of the Rhenish-Westphalian coal district (FRG) and in the Cam pine area (N-Belgium). Ann. Soc. Geol. Bel g. l 09, 367-410. Sclater, J.G. & Christie, P.A.F. 1980. Continental stretching: an explanation of the post Mid-Cretaceous subsidence of the central basin. J. Geophys. Res. 85, 3711-3739. Sweeney, J.J. & Burnham, A.K. 1990. Evaluation of a simple model of vitrinite reflectance based on chemical kinet­ ics. Am. Assoc. Petr. Geol. Bull. 74, 1559-1570. Tys, E. 1980. De geologische structuur van het steenkoolterrein ten noorden van het ontginningsgebied der Kempense mijnen. Belg. Geol. Surv. Prof. Paper 179, 44 pp. Van den haute, P. & Vercoutere, C. 1989. Apatite fission-track evidence for a Mesozoic uplift of the Brabant Massif: preliminary results. Ann. Soc. Geol. Belg. 112, 443-452. Vercoutere, C. & Van den haute, P. 1993. Post-Palaeozoic cooling and uplift of the Brabant Massif as revealed by apatite fission track analysis. Geol. Mag. 130, 639-646. Van Keer, I., Ondrak, R., Muchez, Ph., Bayer, U., Dusar, M. & Viaene, W. in prep. Contrasting burial history and thermal evolution of Westphalian coal-bearing strata in the Cam pine Basin (NE-Belgium): the importance of post­ Paleozoic tectonism and heat flow. Geol. Mijnbouw. Wenselaers, P., Dusar, M. & Van , P.C.H. 1996. Steenkoollaag methaangaswinning in het Kempisch Kolenbekken. AHet proefproject te Peer@. Min. VI. Gem., Afd. Nat. Rijkdom. Energie, Brussel, 67 pp. Ziegler, P.A. 1989. Evolution of Laurussia. A study in Late Palaeozoic Plate Tectonics. Kluwer, Dordrecht, 102 pp. Ziegler, P.A. 1990. Geological atlas of Western and Central Europe. Shell Int. Petr. Mij. and Geol. Soc., London, 2"d ed., 239 pp. Zijerveld, L., Stephenson, R., Cloetingh, S., Duin, E. & Vandenberg, M.W. 1992. Subsidence analysis and modelling of the Roer Valley Graben (SE Netherlands). Tectonophysics 208, 159-171.

I I [ 2 TABLES

13 14 Table 1. Overview of the studied wells.

Nr BGS Archives X-Lambert Y-Lambert Year T.D. (m)

KB69 78W203 229013 187040 1908 1050 KB79 62W205 217403 192704 1906 1133 KB81 64WI76 243766 188775 1910 1026 KB86 77WI74 217105 187372 1920 1912 KBI85 46E278 209390 202578 1986 1480 KBII5 63E200 235060 191412 1943 1730 KB163 63E220 235956 191512 1983 1253 KB174 47E196 220085 199406 1985 1500 KB186 47W264 213939 206366 1986 1504 KB168 63E223 240545 194697 1983 1247 KB172 63E224 234022 196268 1984 1599 KB198 49W226 247659 207752 1987 1773 KB201 48WI91 232836 203184 1989 1340 KBLOMMEL 214500 214000

Table 2. Key stratigraphic units used in Tables 3-16.

WfA Wesphalian A Bathon Bathonian wrn Westphalian B Oxford Oxfordian wrc Westphalian C Sant Santonian WID Westphalian D Camp Campanian St Stephanian Maastr Maastrichtian Zech Zechstein Paleoc Paleocene Bunt Buntsandstein Eoc Eocene Musch Muschelkalk Olig Oligocene Keup Keuper Mioc Miocene Rhaet Rhaetian Plioc Pliocene Hctt Hettangian Quart Quarternary Aalen Aalenian I lower Bajoc Bajocian m middle u upper

15 Table 3. Burial history of well KB69 showing lithologies, thicknesses, ages and thermal maturation data of the different stratigraphic units. Key to lithologies: 0: erosion, 1: shale, 2: sand(stone), 3: chalk, 4: silt(stone). Begin (b)and end (e) ages are according to Harland et al. (1990) and Lippolt et al. (1984). Measured (m) VRm(%) values are from Langenaeker (in prep.). Calculated (c) VRm(%) values correspond to the base of the corresponding stratigraphical units.

Lith. Thickness Age-b Age-e Stratigraphy Vrm%-m VRm%-c

4 450.0 315.0 314.0 lWfA 2.357 4 500.0 314.0 313.0 uWfA 1.68 1.744 4 300.0 313.0 312.0 1WfB 1.26 1.198 4 300.0 312.0 311.0 uWfB 4 260.0. 311.0 310.0 1WfC 4 260.0 310.0 309.0 mWfC 4 260.0 309.0 308.0 uWfC 2 550.0 308.0 307.0 IWfD 2 540.0 307.0 306.0 mWfD 2 300.0 306.0 305.0 uWfD 2 100.0 305.0 303.0 I St 2 80.0 303.0 300.0 u St 0 -2114.0 300.0 240.0 2 440.0 240.0 235.0 Bunt I 70.0 235.0 220.0 Musch I 70.0 220.0 210.0 Keup I 50.0 210.0 208.0 Rhaet I 1400.0 208.0 176.0 Hett/Aalen 0 -2740.0 176.0 84.0 4 60.0 84.0 83.0 Sant/1 Camp 0.313 3 80.0 83.0 74.0 u Camp/! Maastr 0.305 3 60.0 74.0 65.0 u Maastr 0.298 4 112.0 65.0 54.0 Paleoc/1 Eoc 0.293 4 112.0 37.0 25.0 u Eoc/Oiig 0.269 2 75.0 19.0 0.0 I Mioc/Quart 0.000

Table 4. Burial history of well KB79 showing lithologies, thicknesses, ages and thermal maturation data of the different stratigraphic units. Key to lithologies: 0: erosion, 1: shale, 2: sand(stone), 3: chalk, 4: silt(stone). Begin (b) and end (e) ages are according to Harland et al. (1990) and Lippolt et al. (1984). Measured (m) and extrapolated (e) VRm(%) values are from Langenaeker (in prep.). Calculated (c) VRm(%) values correspond to the base of the corresponding stratigraphical units.

Lith. Thickness Age-b Age-e Stratigraphy Vrm%-m VRm%-c

4 500.0 315.0 314.0 1WfA 1.739 4 500.0 314.0 313.0 uWfA 1.19 1.237 4 270.0 313.0 312.0 1WfB 0.83(e) 4 270.0 312.0 311.0 uWfB 4 220.0 311.0 310.0 IWfC 4 230.0 310.0 309.0 mWfC 4 250.0 309.0 308.0 uWfC 2 550.0 308.0 307.0 IWfD 2 500.0 307.0 306.0 mWfD 2 300.0 306.0 305.0 uWfD 2 100.0 305.0 300.0 St 0 -2300.0 300.0 240.0 2 150.0 240.0 235.0 Bunt 1 50.0 235.0 220.0 Musch I 50.0 220.0 210.0 Keup 1 1350.0 210.0 176.0 Rhaet/ Aalen 0 -2000.0 176.0 84.0 4 60.0 84.0 83.0 Santi! Camp 0.322 3 100.0 83.0 74.0 u Camp/! Maastr 0.312 3 50.0 74.0 65.0 u Maastr 0.302 4 157.0 65.0 54.0 Paleoc/1 Eoc 0.298 4 150.0 37.0 25.0 u Eoc/Olig 0.270 2 25.0 19.0 16.0 IMioc 0.243 2 50.0 9.0 6.0 u Mioc 0.233 2 10.0 2.0 0.0 Quart 0.000

16 Table 5. Burial history of well KB81 showing lithologies, thicknesses, ages and thermal maturation data of the different stratigraphic units. Key to lithologies: 0: erosion, 1: shale, 2: sand(stone), 3: chalk, 4: silt(stone). Begin (b) and end (e) ages are according to Harland et al. (1990) and Lippolt et al. (1984). Measured (m) and extrapolated (e) VRm(%) values are from Langenaeker (in prep.). Calculated (c) VRm(%) values correspond to the base of the corresponding stratigraphical units.

Lith. Thickness Age-b Age-e Stratigraphy Vrm%-m VRm%-c

4 450.0 315.0 314.0 IWfA 2.468 4 450.0 314.0 313.0 uWfA 1.829 4 260.0 313.0 312.0 IWfB 1.32 1.315 4 250 .. 0 312.0 311.0 uWfB 1.03(e) 4 335.0 311.0 310.0 IWfC 4 335.0 310.0 309.0 mWfC 4 335.0 309.0 308.0 uWfC 2 550.0 308.0 307.0 IWfD 2 550.0 307.0 306.0 mWfD 2 250.0 306.0 305.0 uWfD 2 100.0 305.0 303.0 I St 2 100.0 303.0 300.0 u St 0 -2482.0 300.0 240.0 2 480.0 240.0 235.0 Bunt I 75.0 235.0 220.0 Musch 1 75.0 220.0 210.0 Keup I 1650.0 210.0 156.0 Rhaet/Bathon 0 -2580.0 156.0 84.0 4 75.0 84.0 83.0 Sant/1 Camp 0.311 3 50.0 83.0 74.0 u Camp/1 Maastr 0.302 3 95.0 74.0 65.0 u Maastr 0.298 4 85.0 65.0 54.0 Paleoc 0.285 4 112.0 37.0 25.0 u Eoc/Olig 0.264 2 40.0 19.0 16.0 IMioc 0.242 2 10.0 2.0 0.0 Quart 0.000

Table 6. Burial history of well KB86 showing lithologies, thicknesses, ages and thermal maturation data of the different stratigraphic units. Key to lithologies: 0: erosion, 1: shale, 2: sand(stone), 3: chalk, 4: silt(stone). Begin (b) and end (e) ages are according to Harland eta!. (1990) and Uppolt et al. (1984). Measured (m)VRm(%) values are from Langenaeker (in prep.). Calculated (c) VRm(%) values correspond to the base of the corresponding stratigraphical units.

Lith. Thickness Age-b Age-e Stratigraphy Vrm%-m VRm%-c

4 480.0 315.0 314.0 IWfA 1.82 1.835 4 500.0 314.0 313.0 uWfA 1.326 4 270.0 313.0 312.0 IWfB 4 270.0 312.0 311.0 uWfB 4 220.0 311.0 310.0 IWfC 4 230.0 310.0 309.0 mWfC 4 250.0 309.0 308.0 uWfC 2 500.0 308.0 307.0 I WID 2 500.0 307.0 306.0 mWfD 2 300.0 306.0 305.0 uWfD 2 100.0 305.0 303.0 I St 2 100.0 303.0 300.0 u St 0 -2340.0 300.0 240.0 2 150.0 240.0 235.0 Bunt I 50.0 235.0 220.0 Musch I 50.0 220.0 210.0 Keup I 50.0 210.0 208.0 Rhaet I 1200.0 208.0 176.0 Hett/Aalen 0 -2150.0 176.0 84.0 4 60.0 84.0 83.0 Sant/1 Camp 0.315 3 110.0 83.0 74.0 u Camp/1 Maastr 0.306 3 30.0 74.0 65.0 u Maastr 0.297 4 120.0 65.0 54.0 Paleoc/1 Eoc 0.293 4 152.0 37.0 25.0 u Eoc/Olig 0.270 2 20.0 19.0 0.0 I Mioc/Quart 0.000

l 7 Table 7. Burial history of well KB185 showing lithologies, thicknesses, ages and thermal maturation data of the different stratigraphic units. Key to lithologies: 0: erosion, 1: shale, 2: sand(stone), 3: chalk, 4: silt(stone). Begin (b)and end (e) ages are according to Harland et al. (1990) and Lippolt et al. (1984). Measured (m) VRm(%) values are from Langenaeker (in prep.). (Values) are measured above the base of a stratigraphic unit. Calculated (c) VRm(%) values correspond to the base of the corresponding stratigraphical units.

Lith. Thickness Age-b Age-e Stratigraphy Vrm%-m VRm%-c

4 500.0 315.0 314.0 !WfA 1.822 4 500.0 314.0 313.0 uWfA ( 1.08) 1.302 4 286.0 313.0 312.0 IWfB 0.94 0.901 4 250.,0 312.0 311.0 uWfB 0.74 0.764 4 220.0 311.0 310.0 lWfC 4 230.0 310.0 309.0 mWfC 4 250.0 309.0 308.0 uWfC 2 550.0 308.0 307.0 IWfD 2 500.0 307.0 306.0 mWfD 2 300.0 306.0 305.0 uWfD 2 150.0 305.0 300.0 St 0 -2396.0 300.0 240.0 2 370.0 240.0 235.0 Bunt I 60.0 235.0 220.0 Musch I 60.0 220.0 210.0 Keup I 1550.0 210.0 156.0 Rhaet/Bathon 0 -2090.0 156.0 84.0 4 65.0 84.0 83.0 Sant/1 Camp 0.334 3 133.0 83.0 74.0 u Camp/1 Maastr 0.322 3 50.0 74.0 65.0 u Maastr 0.305 4 174.0 65.0 53.0 Paleoc/1 Eoc 0.301 4 165.0 37.0 25.0 u Eoc/Oiig 0.273 2 35.0 19.0 16.0 I Mioc 0.244 2 125.0 9.0 5.0 u Mioc/1 Plioc 0.234 2 10.0 3.0 0.0 u Plioc/Quart 0.000

Table 8. Burial history of well KB115 showing lithologies, thicknesses, ages and thermal maturation data of the different stratigraphic units. Key to lithologies: 0: erosion, 1: shale, 2: sand(stone), 3: chalk, 4: silt(stone). Begin (b) and end (e) ages are according to Harland et al. (1990) and Lippolt et al. (1984). Measured (m) and extrapolated (e) VRm(%) values are from Langenaeker (in prep.). (Values) are measured above the base of a stratigraphic unit. Calculated (c) VRm(%) values correspond to the base of the corresponding stratigraphical units.

Lith. Thickness Age-b Age-e Stratigraphy Vrm%-m VRm%-c

4 450.0 315.0 314.0 lWfA 2.088 4 450.0 314.0 313.0 uWfA ( 1.54) 1.653 4 395.0 313.0 312.0 IWfB 1.36 1.31 I 4 290.0 312.0 311.0 uWfB 1.10 1.070 4 250.0 311.0 310.0 IWfC 0.855 0.898 4 250.0 310.0 309.0 mWfC 0.67(e) 4 250.0 309.0 308.0 uWfC 2 500.0 308.0 307.0 IWfD 2 500.0 307.0 306.0 mWfD 2 350.0 306.0 300.0 u WfD/St 0 -1852.0 300.0 250.0 4 25.0 250.0 245.0 u Zech 2 490.0 245.0 235.0 Bunt I 70.0 235.0 220.0 Musch I 70.0 220.0 210.0 Keup I 1680.0 210.0 170.0 Rhaet/Bajoc 0 -2335.0 170.0 84.0 4 65.0 84.0 83.0 Sant/1 Camp 0.319 3 60.0 83.0 72.0 u Camp/1 Maastr 0.308 3 75.0 72.0 65.0 u Maastr 0.302 4 132.0 65.0 54.0 Paleoc/1 Eoc 0.296 4 133.0 37.0 25.0 u Eoc/Oiig 0.269 2 65.0 19.0 16.0 I Mioc 0.246 2 10.0 9.0 6.0 u Mioc 0.232 2 27.0 2.0 0.0 Quart 0.000

I 8 Table 9. Burial history of well KB163 showing lithologies, thicknesses, ages and thermal maturation data of the different stratigraphic units. Key to lithologies: 0: erosion, 1: shale, 2: sand(stone), 3: chalk, 4: silt(stone). Begin (b) and end (e) ages are according to Harland et al. (1990) and Lippolt et al. (1984). Measured (m) and extrapolated (e) VRm(%) values are from Langenaeker (in prep.). (Values) are measured above the base of a stratigraphic unit. Calculated (c) VRm(%) values correspond to the base of the corresponding stratigraphical units.

Lith. Thickness Age-b Age-e Stratigraphy Vrm%-m VRm%-c

4 450.0 315.0 314.0 IWfA 1.902 4 450.0 314.0 313.0 uWfA 1.495 4 300.0 313.0 312.0 IWfB ( 1.13) 1.185 4 275.0 312.0 311.0 uWfB 1.02 1.006 4 250.0 31 1.0 310.0 IWfC 0.85 0.866 4 250.0 310.0 309.0 mWfC 0.69(e) 4 250.0 309.0 308.0 uWfC 2 500.0 308.0 307.0 I WID 2 500.0 307.0 306.0 mWfD 2 300.0 306.0 300.0 u WfD/St 0 -1801.0 300.0 250.0 4 20.0 250.0 245.0 u Zech 2 460.0 245.0 235.0 Bunt I 60.0 235.0 220.0 Musch l 60.0 220.0 210.0 Keup l 1650.0 210.0 170.0 Rhaet/Bajoc 0 -2250.0 170.0 84.0 4 66.0 84.0 83.0 Sant/1 Camp 0.319 3 58.0 83.0 72.0 u Campll Maastr 0.308 3 80.0 72.0 65.0 u Maastr 0.302 4 124.0 65.0 53.0 Palcocll Eoc 0.295 4 133.0 37.0 25.0 u Eoc/Oiig 0.269 2 63.0 19.0 16.0 IMioc 0.244 2 22.0 9.0 5.0 u Mioc!l Plioc 0.232 2 14.0 3.0 0.0 u Plioc/Quart 0.000

Table 10. Burial history of well KB174 showing lithologies, thicknesses, ages and thermal maturation data of the different stratigraphic units. Key to lithologies: 0: erosion, 1: shale, 2: sand(stone), 3: chalk, 4: silt(stone). Begin (b) and end (e) ages are according to Harland et al. (1990) and Lippolt et al. (1984). Measured (m) VRm(%) values are from Langenaeker (in prep.). (Values) are measured above the base of a stratigraphic unit. Calculated (c) VRm(%) values correspond to the base of the corresponding stratigraphical units.

Lith. Thickness Age-b Age-e Stratigraphy V1·m%-m VRm%-c

4 500.0 315.0 314.0 lWfA 1.653 4 450.0 314.0 313.0 uWfA (0.98) 1.179 4 287.0 313.0 312.0 lWfB 0.88 0.872 4 250.0 312.0 311.0 uWfB 0.74 0.758 4 220.0 311.0 310.0 IWfC 4 220.0 310.0 309.0 mWfC 4 220.0 309.0 308.0 uWfC 2 500.0 308.0 307.0 !WID 2 500.0 307.0 306.0 mWfD 2 350.0 306.0 300.0 u WfD/St 0 -2111.0 300.0 250.0 4 27.0 250.0 245.0 u Zech 0.709 2 430.0 245.0 235.0 Bunt 0.702 I 80.0 235.0 220.0 Musch I 80.0 220.0 210.0 Keup 1 1550.0 210.0 156.0 Rhaet/Bathon 0 -2049.0 156.0 84.0 2 71.0 84.0 78.0 Sant/1 Camp 0.324 3 89.0 78.0 70.0 u Camp/! Maastr 0.319 3 47.0 70.0 65.0 u Maastr 0.308 4 188.0 65.0 50.0 Paleoc/1 Eoc 0.303 4 193.0 37.0 25.0 u Eoc/Olig 0.277 2 58.0 19.0 16.0 I Mioc 0.246 2 19.0 9.0 6.0 uMioc 0.233 2 64.0 5.0 0.0 Plioc/Quart 0.000

19 Table 11. Burial history of well KB186 showing lithologies, thicknesses, ages and thermal maturation data of the different stratigraphic units. Key to lithologies: 0: erosion, 1: shale, 2: sand(stone), 3: chalk, 4: silt(stone). Begin (b) and end (e) ages are according to Harland et al. (1990) and Lippolt et al. (1984). Measured (m) VRm(%) values are from Langenaeker (in prep.). (Values) are measured above the base of a stratigraphic unit. Calculated (c) VRm(%) values correspond to the base of the corresponding stratigraphical units.

Lith. Thickness Age-b Age-e Stratigraphy Vrm%-m VRm%-c

4 500.0 315.0 314.0 IWfA 1.852 4 500.0 314.0 313.0 uWfA 1.408 4 260.0 313.0 312.0 IWfB 1.088 4 260.0 312.0 311.0 uWfB 0.96 0.933 4 200.0 311.0 310.0 !WfC 0.82 0.832 4 250.0 310.0 309.0 mWfC (0.73) 0.766 4 250.0 309.0 308.0 uWfC 0.708 2 550.0 308.0 307.0 IWfD 2 500.0 307.0 306.0 mWfD 2 400.0 306.0 300.0 u WfD/St 0 -1695.0 300.0 250.0 4 20.0 250.0 245.0 uZech 0.706 2 400.0 245.0 235.0 Bunt I 70.0 235.0 220.0 Musch I 70.0 220.0 210.0 Keup I 1580.0 210.0 156.0 Rhaet/Bathon 0 -2120.0 156.0 84.0 4 67.0 84.0 83.0 Santi! Camp 0.340 3 102.0 83.0 74.0 u Camp/1 Maastr 0.329 3 47.0 74.0 65.0 u Maastr 0.314 4 240.0 65.0 54.0 Paleoc 0.308 4 163.0 37.0 25.0 u Eoc/Olig 0.273 2 45.0 19.0 16.0 !Mioc 0.246 2 159.0 9.0 5.0 u Mioc/1 Plioc 0.235 2 10.0 3.0 0.0 u Plioc/Quart 0.000

Table 12. Burial history of well KBLommel showing lithologies, thicknesses, ages and thermal maturation data of the different stratigraphic units. Key to lithologies: 0: erosion, 1: shale, 2: sand(stone), 3: chalk, 4: silt(stone). Begin (b) and end (e) ages are according to Harland eta!. (1990) and Lippolt eta!. (1984). Measured (m) VRm(%) values are from Langenaeker (in prep.). Calculated (c) VRm(%) values correspond to the base of the corresponding stratigraphical units.

Lith. Thickness Age-b Age-e Stratigraphy VRm%-c

4 500.0 315.0 314.0 IWfA 2.255 4 500.0 314.0 313.0 uWfA 1.820 4 260.0 313.0 312.0 IWfB 1.447 4 260.0 312.0 311.0 uWfB 1.282 4 200.0 311.0 310.0 lWfC 1.135 4 250.0 310.0 309.0 mWfC 1.019 4 250.0 309.0 308.0 uWfC 0.888 2 550.0 308.0 307.0 lWfD 0.795 2 500.0 307.0 306.0 mWfD 2 400.0 306.0 300.0 u WfD/St 0 -1350.0 300.0 250.0 4 20.0 250.0 245.0 u Zech 0.776 2 400.0 245.0 235.0 Bunt 0.771 I 90.0 235.0 220.0 Musch I 90.0 220.0 210.0 Keup 1 1800.0 210.0 144.0 Rhaet/Oxford 0 -2280.0 144.0 84.0 2 80.0 84.0 78.0 Santll Camp 0.337 3 40.0 78.0 70.0 u Camp/! Maastr 0.331 3 55.0 70.0 65.0 u Maastr 0.325 4 300.0 65.0 50.0 Paleoc/l Eoc 0.317 4 140.0 37.0 25.0 u Eoc/Oiig 0.272 2 80.0 19.0 16.0 IMioc 0.249 2 180.0 9.0 6.0 uMioc 0.237 2 50.0 5.0 0.0 Plioc/Quart 0.000

20 Table 13. Burial history of well KB168 showing lithologies, thicknesses, ages and thermal maturation data of the different stratigraphic units. Key to lithologies: 0: erosion, 1: shale, 2: sand(stone), 3: chalk, 4: silt(stone). Begin (b) and end (e) ages are according to Harland et al. (1990) and Lippolt et al. (1984). Measured (m) VRm(%) values are from Langenaeker (in prep.). Calculated (c) VRm(%) values correspond to the base of the corresponding stratigraphical units.

Lith. Thickness Age-b Age-e Stratigraphy Vrm%-m VRm%-c

4 450.0 315.0 314.0 IWfA 2.516 4 450.0 314.0 313.0 uWfA 2.086 4 300.0 313.0 312.0 IWfB 1.683 4 300.0 312.0 311.0 uWfB 1.443 4 250.,0 311.0 310.0 IWfC 1.231 4 250.0 310.0 309.0 mWfC 0.925 1.085 4 250.0 309.0 308.0 uWfC 0.924 2 500.0 308.0 307.0 I WID 0.77 0.814 2 500.0 307.0 306.0 mWfD 2 250.0 306.0 300.0 u WfD/St 0 -1168.0 300.0 250.0 4 20.0 250.0 245.0 uZech 2 460.0 245.0 235.0 Bunt I 70.0 235.0 220.0 Musch I 70.0 220.0 210.0 Keup I 1700.0 210.0 170.0 Rhaet/Bajoc 0 -2320.0 170.0 84.0 4 85.0 83.0 78.0 Sant/1 Camp 0.327 3 80.0 78.0 70.0 u Camp/1 Maastr 0.311 3 59.0 70.0 67.0 u Maastr 0.302 4 134.0 67.0 54.0 Paleoc 0.297 4 151.0 37.0 25.0 u Eoc/Olig 0.272 2 85.0 19.0 16.0 I Mioc 0.246 2 41.0 9.0 5.0 u Mioc/1 Plioc 0.233 2 18.0 3.0 0.0 u Plioc/Quart 0.000

Table 14. Burial history of well KB172 showing lithologies, thicknesses, ages and thermal maturation data of the different stratigraphic units. Key to lithologies: 0: erosion, 1: shale, 2: sand(stone), 3: chalk, 4: silt(stone). Begin (b) and end (e) ages are according to Harland et al. (1990) and Lippolt et al. (1984). Measured (m) VRm(%) values are from Langenaeker (in prep.). Calculated (c) VRm(%) values correspond to the base of the corresponding stratigraphical units.

Lith. Thickness Age-b Age-e Stratigraphy Vrm%-m VRm%-c

4 450.0 315.0 314.0 IWfA 2.648 4 450.0 314.0 313.0 uWfA 2.207 4 270.0 313.0 312.0 IWfB 1.790 4 260.0 312.0 311.0 uWfB 1.532 4 250.0 311.0 310.0 IWfC 1.32 1.323 4 240.0 310.0 309.0 mWfC 1.172 4 240.0 309.0 308.0 uWfC 1.039 2 500.0 308.0 307.0 I WID 0.93 0.903 2 450.0 307.0 306.0 mWfD 2 400.0 306.0 300.0 u WfD/St 0 -1225.0 300.0 250.0 4 22.0 250.0 245.0 u Zech 0.44 0.872 2 540.0 245.0 235.0 Bunt 0.864 I 85.0 235.0 220.0 Musch I 85.0 220.0 210.0 Keup I 1800.0 210.0 156.0 Rhaet/Bathon 0 -2381.0 156.0 84.0 4 78.0 84.0 83.0 Sant/1 Camp 0.39 0.329 3 67.0 83.0 74.0 u Camp/1 Maastr 0.314 3 77.0 74.0 65.0 u Maastr 0.305 3 155.0 65.0 54.0 Paleoc 0.298 4 142.0 37.0 25.0 u Eoc/Olig 0.271 2 106.0 19.0 16.0 IMioc 0.249 2 46.0 9.0 5.0 u Mioc/1 Plioc 0.233 2 10.0 3.0 0.0 u Plioc/Quart 0.000

2 1 Table 15. Burial history of well KB201 showing lithologies, thicknesses, ages and thermal maturation data of the different stratigraphic units. Key to lithologies: 0: erosion, 1: shale, 2: sand(stone), 3: chalk, 4: silt(stone). Begin (b)and end (e) ages are according to Harland et al. (1990) and Lippolt et al. (1984). Measured (m) and extrapolated (e) VRm(%) values are from Langenaeker (in prep.). Calculated VRm(%) (c) values correspond to the base of the corresponding stratigraphical units.

Lith. Thickness Age-b Age-e Stratigraphy Vrm%-m VRm%-c

4 450.0 315.0 314.0 1WfA 2.889 4 450.0 314.0 313.0 uWfA 2.421 4 300.0 313.0 312.0 !WfB 1.993 4 300.0 312.0 311.0 uWfB 1.732 4 260.0 311.0 310.0 1WfC 1.477 4 260.0 310.0 309.0 mWfC 1.294 4 260.0 309.0 308.0 uWfC 1.128 2 550.0 308.0 307.0 !WID 1.05(e) 2 400.0 307.0 306.0 mWID 2 400.0 306.0 300.0 u WID/St 0 -1360.0 300.0 250.0 4 39.0 250.0 245.0 u Zech 0.47 0.961 2 600.0 245.0 235.0 Bunt 0.939 1 85.0 235.0 220.0 Musch 1 85.0 220.0 210.0 Keup 1 1900.0 210.0 156.0 Rhaet/Bathon 0 -2216.0 156.0 84.0 4 90.0 84.0 83.0 Sant/1 Camp 0.337 3 75.0 83.0 74.0 u Camp/1 Maastr 0.319 3 70.0 74.0 65.0 u Maastr 0.308 4 150.0 65.0 54.0 Paleoc/1 Eoc 0.301 4 177.0 37.0 25.0 u Eoc/01ig 0.278 2 120.0 19.0 16.0 1Mioc 0.250 2 75.0 9.0 5.0 u Mioc/1 Plioc 0.234 2 49.0 3.0 0.0 u Plioc/Quart 0.000

Table 16. Burial history of well KB198 showing lithologies, thicknesses, ages and thermal maturation data of the different stratigraphic units. Key to lithologies: 0: erosion, 1: shale, 2: sand(stone), 3: chalk, 4: silt(stone). Begin (b) and end (e) ages are according to Harland et al. (1990) and Lippolt et al. (1984). Measured (m) .VRm(%) values are from Langenaeker (in prep.). (Values) are measured above the base of a stratigraphic unit. Calculated (c) VRm(%) values correspond to the base of the corresponding stratigraphical units.

Lith. Thickness Age-b Age-e Stratigraphy Vrm%-m VRm%-c

4 450.0 315.0 314.0 1WfA 3.567 4 450.0 314.0 313.0 uWfA 3.120 4 330.0 313.0 312.0 IWfB 2.628 4 330.0 312.0 311.0 uWfB 2.285 4 250.0 311.0 310.0 IWfC 1.958 4 250.0 310.0 309.0 mWfC 1.736 4 250.0 309.0 308.0 uWfC 1.508 2 550.0 308.0 307.0 !WID 1.324 2 500.0 307.0 306.0 mWID 2 200.0 306.0 305.0 uWID 2 50.0 305.0 303.0 l St 2 50.0 303.0 300.0 u St 0 -800.0 300.0 250.0 4 40.0 250.0 245.0 uZech 1.117 2 700.0 245.0 235.0 Bunt 1.095 I 100.0 235.0 220.0 Musch 0.862 I 100.0 220.0 210.0 Keup 0.821 I 1980.0 210.0 144.0 Rhaet/Oxford (0.46) 0.782 0 -1280.0 144.0 74.0 3 27.0 74.0 65.0 Maastr. 0.30 0.333 4 150.0 65.0 54.0 Paleoc/1 Eoc 0.325 4 333.0 37.0 25.0 u Eoc/Olig 0.305 2 369.0 19.0 16.0 lMioc 0.259 2 169.0 9.0 5.0 u Mioc/1 Plioc 0.237 2 236.0 3.0 0.0 u Plioc/Quart 0.000

22 FIGURES

23 N @

4-KB

KB186 Bree KB185 KB201 + • + • Mo.o.seik KB 174 + KB172 Beringen • KBI68 + KB79 KB 115 · + ++ KBI63 KBSI +KB86 + KB69 + Genk ~------~------~--~· • Mo.o.sMechelen

\ \ \ " 11' ' "'-..._ • ossrf'...... _ ....___,.,,...... _~ \ Brussels " ' -- ...______

N @ co

0 10 20 kl""l 0 100 KM

Fig. I. Location of the studied wells. Conesponding localities for the \veils are Winterslag-Genk (KB69), Zolder (KB79, KB86), Eisden (KB81), Opglabbeek (KBI15), Opoeteren (KBI63, KB168). Gruitrode (KBI72), Bechtel (KDI74), Oostham (KBI85), Lommel (KB186, KBLom.), Molenbeersel (KB 198) and Bree (KB20 I). SURFACE TEMPERATURE MODEL

30.------~------~------~------~------~------~~----~

25 ...... , ...... ! ·······················r···························r···························-,...... , ...... 6 ; ; ; : i i ~ "••••••••u••h•••••••o•o•OOoo~ooo•oO•h••••oooouo••••••oo• B

0+------r------;------+------,_I I I I ______~------~------~ I I -350 -300 -250 -200 -150 -1 00 -50 0 Time (myr)

Fig. 2. Surf.1ce temperature model for the Campine Basin, showing the evolution of temperature in time. N

N ® ® 10 20 kM 0 100 kM

Fig. 3. Isopach map for Penn ian through Jurassic strata in the east em Campine Basin. KB69 Burial history

-5000·~------~------r------;------~------+------~------~ a -350 -300 -250 -200 -150 -100 -50 0 Time (myr)

KB69 Evolution heat flow

160.------~------~------~------~------~------~------~

1 40 ···························+··························f .. ······ .. ··················[······ .. ····················t···························f····························[····························

E 120 ...... i...... [ ...... f ...... j...... f····························f···························· I1:: L :j:~ --1 : :_-:t : : :: j-_ ~~ ::J:: ~: --~

: r : :r --1 : :::r : r: : ~ 1 : 20 ...... ) ...... ) ...... ) ...... t...... ) ...... ) ...... o+------r------r------;------+------r------;------11 l ~ ~ l ~ b -350 -300 -250 -200 -150 -100 -50 0 Time (myr)

KB69 Burial temperature evolution

2501~------,------~------,------~------~~------~------,

200 ·······················

6 a :s.(!} ci E 1 oo ~ s 50 w Q ___...... ;...-~-..-...... ::~K 01+-----~~------4------~------4------~--~--L-~------~ T -350 -300 -250 -200 -150 -1 00 -50 0 c Time (myr)

Fig. -L Bmial history (a), evolution heat flow (b) and bmial temperature (c) in the deep well KB69. Key to stratigraphic markers: S: Samsbank M.B., W: Wasser£1ll M.B., Q: Quaregnon M.B., K: Base Cretaceous, T: Base Teitiaiy. KB79 Burial history

o•~----~~------~------~------~------~-c---c~------T -500 ...... ······r························-r··························l""" ...... ~·························· K ~ ~~~ :.::. !: i i :: :::. ri i I:- I : : ~v g :o: ----•••-•-•·••---•······ t·-··-----•••• :t :: r rl : J :• -4ooo ...... L ...... i...... \...... \...... \·-·····················:····f··························· -4500 ···························i·········· .. ················f···························l······· .. ···················f···························l····························f··························· ; ; ; ! ~ ~ -5000+------~.------~------~~------r------r------+------~ a -350 -300 -250 -200 -1 50 -1 00 -50 0 Time (myr)

KB79 Evolution heat flow

-350 -300 -250 -200 -150 -100 -50 0 b Time (myr)

KB79 Burial temperature evolution

250•,------~------~------~------~------~~------~------~

200

c. ;::c 1 00] ..... I-

50 I

o! I -350 -300 -250 -200 -150 -100 Time (myr) c

Fig. 5. Blllial histor~y (a), evolution heat flow (b) and bmial temperature (c) in the deep \Veil KB79. Kev to stratigraphic markers: S: Samsbank M.B., W: Wasserfall M.B., K: Base Cretaceous, T: Base Tertiary~ KB81 Burial history 0~--!f~----~----~------r-----~~~~~==~T -500 ··················· ...... , ...... K -1000 ...... ····-··························· Q w s

-5000-~------r------r------r------~------r------r------~ -350 -3oo -250 -200 -150 -1 oo -so o a Time (myr)

KB81 Evolution heat flow

160.------~------~------~------~------~------~------,

1 40 ····························t .. ·························t····························(························-t-··························t····························j···························· ~ :: •.••.••••..•.••••••• j j :: :j:::=: :=!:::= :_j : = :I :: : f 80 ···························+···· ··+ ···1--- -f --1--+ ················· 60 ····························(··························(··························j····························j····························t····························(··················~ 40 ...... ) ...... ~ ...... ) ...... ~ ...... ) ...... ) ...... 20 ...... L ...... [ ...... J...... l...... l...... l ...... l 1 0+.------+------T------+------T------T------T------~~ ~ ~ ~ -350 -300 -250 -200 -150 -100 -50 0 b Time (myr)

KB81 Burial temperature evolution

250.------~------~------~------~------.~------,------.

200 ......

0 0, :s.

Fig. 6. B uri a! history (a), evolution heat flow (b) and burial temperature (c) in the deep well KBSI. Key to stratigraphic markers: S: Samsbank M.B., W: Wasserfall M.B., Q: Quaregnon M.B., K: Base Cretaceous, T: Base Tertiary. KB86 Burial history

T K w s

,.... c. OJ 0

-4500 ...... -r························--r·························r·························r························r························-r··························

-50001 +------+------r------r------~------.------.------~ -350 -300 -:250 -:200 -1 50 -1 00 -50 0 a Time (myr)

KB86 Evolution heat flow

160.------~------~------~------~------~------.

E : ::r : I t: r: : l T I : l1 OOJ ...... ;. ···························t····························l····························l····························-f····························l···························· ; 80 ...... ;...... ; ...... 1····························1-··························+·························+···························

i :: :.. ·::.:::·::::.·::· .. ·:.::.~:.·:::··:.::·::::·.::··:.:··r::::.:·::::::::·::::·::::r::.::::.::::::::::::::.:::L::·:·:::::::::::::.·::::I::····.::·:··· .. ::.::::·::r::::.:.:.:.::.:.:: ...... ~~ · , r L r I -350 -300 -:250 -:200 -150 -1 00 -50 0 b Time (myr)

KB86 Burial temperture evolution 250 I . , , , , .

...... s I' : w ~ : o, " . ~-? -350 -300 -:250 -:200 -150 -100 -50 0 Time (myr) c

Fig. 7. Bwial hist01y (a), c\·olution heat flow (b) and burial temperature (c) in the deep well KB86. Key to stratigraphic markers: S: Samsbank M.B., TV: Wasser£1ll M.B., K: Base Cretaceous. T: Base Te1tiary. KB185 Burial history

T E.K Q g s c. -2500 C,) 0 -3000 -3500 -4000 -4500

-5000' -350 -300 -250 -200 -150 -100 -50 0 a Time (myr)

KB185 Evolution heat flow

Q ; j 1 • I -3SO -300 -250 -200 -1 50 -1 oo -so 0 b Time {myr)

KB185 Burial temperature evolution

250.------~------,------,------,------,,------,------~ i :: ~r I I: I , E co 1-- s Q LK T 0 I , -350 -300 -250 -200 Time (myr) c

Fig. 8. Bmial history (a), evolution heat flow (b) and burial temperature (c) in the deep well KB 185. Ke_y to stratigraphic markers: S: Samsbank M.B., Q: Quaregnon M.B., E: Eisden M.B., K: Base Cretaceous. T: Base Teniary. KB115 Burial history 0 T -500 K -1000 ••••••u•••••••••••• M -1500 Q g -2000 ~ s c. -2500 "-' 0 -3000 -3500 -4000 -4500 -5000 a -350 -300 -250 -200 -150 -100 -50 Time (myr)

KB115 Evolution heat flow

160~------,------,------~------~------~------~------~

140 ············· ••)•• ························-:····························'····························;·····························[··························· :···························· r 1 ,_ _I r:3 t ' _ g

-35o -3oo -2so -2oo -1so -1 oo -5o o Time (myr) b

KB115 Burial temperature evolution

2501~------~------~------~~------~------~------~------~

200

s Q M K T

Time (myr) c fig. 9. Burial histmy (a), evolution heat flow (b) and bmial temperature (c) in the deep \Veil KB 115. Key to stratigraphic markers: S: Samsbank M.B., Q: Quaregnon M.B., M: Maurage M.B .. K: Base Cretaceous, T: Base Te1tiary. KB163 Burial history 0 T -500 K -1000 ...... M -1500 Q g -2000 s c.~ -2500 <:.> 0 -3000

-4500 ··················································

-5000+------~------r------r------r------r------.------, -350 -300 -250 -200 -150 -1 00 -50 0 Time (myr) a

KB163 Evolution heat flow

160.------~~------~------~------~------~------~------140 : : ' . ' '

~ , 201••••••••••••••••••••••••••••:•······; r F r : r•••••••••·····•

!'~rI J. :J: . l . ::1:: . I: . :: . : :r · r , r r 1- -350 -300 -250 -200 -150 -100 -50 0 b Time (myr)

KB163 Burial temperature evolution

2501.------~------,------~------~------~~------~------,

200

6 cs 150 :s<:.> ~ E 1 oo ...... Q I- ] I -350 -300 -250 -200 -150 Time (myr) c

Fig. !0. Burial history (a), evolution heat flow (b) and burial temperature (c) in the deep well KB 163. Key to stratigraphic markers: S: Samsbank M.B., Q: Quaregnon M.B., 1vf: Maurage M.B .. /\: Base Cretaceous, T: Base Tettiary. KB174 Burial history 0 T -500 ...... K -100 E -1500 Q g -2000 ,... s c. -2500 C) 0 -3000 -3500 -4000 -4500 -5000 a -350 -300 -250 -200 Time (myr)

KB174 Evolution heat flow

160~------~------~------~------~

E ~ :: :::·::·::·::·: __ ::_·:: ___ :__ ;-·_:::::::::·:_:.o:::_:-·:_:r::··_:_-_::::·::::_:·::-::~_·::::::.::_:·:::.::::::::::~:.:::::::··::::::::::·:::::.::.:::::::::::::::::::::::·:r::::::::::::::::::::::::::

r••·••• I: :l: ::J••:t :=:••••• !~ ,:: 60 •••••••••••••••••••·•••••··················------<------T·------T·------r------r·------r------":1 ; r r ::: r i : -350 -300 -250 -200 -1 50 -1 00 -50 0 Time (myr) b

KB174 Burial temperature evolution

2501,------~------~------~------~------~~------~------~

200 ...... • ...... j··························+··························r·················· , ...... ··················:····· ......

6 ~ 3

-350 -300 -250 -200 -150 -1 00 -50 0 Time (myr) c

Fig. ll. Bmial history (a), evolution heat flow (b) and burial temperature (c) in the deep \\'dl KB 174. Key to stratigraphic markers: S: Samsbank M.B., Q: Quaregnon M.B., E: Eisdcn M.B. K: Base Cretaceous, T: Base Tertiary. KB186 Burial History ·

T K M Q c.~ ~ s 0 -3000 -3500 -4000 -4500 -50007·~------~------~------~~------~------+------~------~ -350 -300 -250 -200 -150 -100 -50 0 a Time (myr)

KB186 Evolution heat flow

b -350 -300 -250 -200 -150 -1 00 -50 0 Time (myr)

KB186 Burial temperature evolution

250~------;------,------,------,,------~------~------~

200

0 ~ sc. ~ ']············ s Q M K J T -350 -250 -200 -150 -100 -50 0 Time (myr) c

Fig. 12. Burial history (a), evolution heat flow (b) and burial temperature (c) in the deep \veil KB 186. Key to stratigraphic markers: S: Samsbank M.B., Q: Quaregnon M.B., lv/: Maurage M.B., K: Base Cretaceous, T: )3ase Tertiary. KBLOMMEL Burial history 0 -500 T K -1000 N -1500 -2000 g Q .r:::c. -2500 C) 0 -3000 s -3500 -4000 -4500

-5000' -350 -300 -2so -200 -150 -1 00 Time (myr) a

KBLOMMEL Evolution heat flow 160

140

E 120 c- -...."' 100 $: .§. 80 :;: g 60 t5 C) I 40

20

0 -350 -300 -250 -200 0 Time (myr) b

KBLOMMEL Burial temperature evolution

250•.------~------~------,------,------,------,------~

200

G 0 150 E ci. E 1 oo C) I- s Q 50 ··············· N K T 0. -350 -300 -250 -200 -150 -1 00 -50 0 Time (myr) c

Fig. 13. Burial histOI}' (a), evolution heat flow (b) and burial temperature (c) in the deep well KBLommel. Key to stratigraphic markers: S: Samsbank M.B., Q: Quaregnon M.B., N: Base Neeroeteren Sst., K: Base Cretaceous, T: Base Tertiary. KBI68 Burial history

T K N

~ -2000 Q Ci. -2500 ~ -3000 ...... ···'··········· ···············•········ s

::+ r _! : 1 : +: r r -5000-350 -300' -250' -200 -1 50 -1 00 -50 0 Time {myr) a

KBI68 Evolution heat flow

E c- "' .s~ g3 c:; ~ I

0' ' ' -350 -300 -250 -200 -1 50 -1 00 -50 0 Time (myr) b

KBI68 Burial temperature evolution

2501,------~------~------~------,------,,------~------~

s Q

N.K T

-300 -250 -200 -150 -100 -50 0 Time (myr) c

Fig. 14. Buiial history (a), evolution heat flow (b) and buiial temperature (c) in the deep well KB 168. Key to stratigraphic markers: S: Samsbank. M.B., Q: Quaregnon M.B., N: Base Neeroeteren Sst.. K: Base Cretaceous, T: Base Tertiary. KB172 Burial history

Or~----~~------~----~~-----,------:--c-~~------1 -500 ...... ···-················· T K N

Q .r:::.

c.dl 0 s

-4000 -4500 - 500? ~5-o------3~oo------245_o______2oro------1~5-o------1~o~o~------_5To~------~o 3 a Time (myr)

KB172 Evolution heat flow

1601,------~------~------,------,------,------,------~

E 120 ci- "' 100 ~ .s 80 3: ,g 60 d C) I 40

20

or+------~------+------~.------r------T------~------~ -350 -300 -250 -200 -150 -1 00 -50 0 b Time (myr)

KB172 Burial temperature evolution 250

200

0 ....: Cl 150 dl 2 ci E 100 C) 1- s Q 50 ...... N K T 0 -350 -300 -250 -200 -1 50 -1 00 -50 0 Time (myr) c

Fig. 15. Burial histo1y (a), evolution heat flow (b) and burial temperature (c) in the deep well KB 172. Key to stratigraphic markers: S: Samsbank M.B., Q: Quaregnon M.B., N: Base Neeroeteren Sst.. K: Base Cretaceous, T: Base Tertiary. KB201 Burial history

T K

M Q

s

-50001+------T------r------r------~------r------r------~ -350 -300 -250 -200 -1 50 -50 0 a Time (myr)

KB201 Evolution heat flow

1601.------~------~------~------~------~------~------,

J : : : : : : :+ I I i T j :r: -350 -300 -250 -200 -150 -1 00 -50 0 Time (myr) b

KB201 Burial temperature evolution

250.------~------~------~------~------~------~------,

200 ··················································································

,., s Q M g ,:r1r························;·········· ...... K T 0' ' . . -350 -300 -250 -200 -150 -1 00 -50 0 Time (myr) c

Fig. 16. Burial history (a), evolution heat flow (b) and burial temperature (c) in the deep well KB20 I. Key to stratigraphic markers: S: Samsbank M.B., Q: Quaregnon M.B., M: Maurage M.B., K: Base Cretaceous, T: Base Tertiary. KB198 Burial history

0~--~\T------,-----;-----~------~--~~~--~

-1000 ···················· .. ·-r·········

~ -3000 ...... T...... ··r ...... N

Q ::: ••••••••••••••••••• 1:••••••••••• r -~········ s -6000' -350 -300 -250 -200 -150 -100 -50 0 Time (myr) a

KB198 Evolution heat flow

160•.------~------~------~------~------~------~------,

140 ·················. ········~············· ···············;····························t····························j·····························i·········· ·················:····························

-350 -300 -250 -200 -150 -100 -50 0 Time (myr) b

KB198 Burial temperature evolution 250

200

6 150 s a<:.> 2. Q c. E 100 N <:.> f-

50 K,T

o. -350 -300 -250 -200 -1 50 -1 00 -50 0 Time (myr) c

Fig. 17. Burial histOJy (a), evolution heat flow (b) and burial temperature (c) in the deep well KBI98. Key to stratigraphic markers: S: Samsbank M.B., Q: Quaregnon M.B., N: Base Neeroeteren SsL K: Base Cretaceous, T: Base Tertiary. KB69 Evolution VRm% Sarnsbank M.B.

3,..,.----~-----,-----..,..------,----:----...,----,o.4 ! 0.35

0.3

0.25 ~ 'E~ .§ a: 02. 0~ E > a: 0.15 >

-350 -300 -250 -200 -150 -1 00 -50 0 Time (myr) a

KB69 Evolution VRm% Wasserfall M.B. :...----..,.------,-----,------,----;-----;---T0.4

;.:::::::::.::.:.:.::.:::::::.:::;=:.:.:.::.:::.:.:.::.:::::::.::::;=::::::.::.:.:::::.=+.::.:.....:------+------r-----ro.35

0.3

0.25 5o 6'< .§ E a: 0.2 ~ > E

-350 -300 -250 -200 -150 -100 -50 0 Time (myr) b

KB69 Evolution VRm% Quaregnon M.S. 2 0.4 ~ ;:~ ~: ~~~ ~: :I : := :r:= : :-l- := : =:=1:=:-: =:T: :::==::r::=- ~= ::: ~ OE 0.2 ~ a: a: > 0.15 >

0.1

0 . . ' -3'so -3oo -2'so -200 -1 50 -1 oo -50 o Time (myr) c

Fig. 18. Evolution ofVRm(%) (continuous line) and VRm(%)/myr (dashed line) in the deep well KB69 (a: Samsbank M.B., b: Wasserfall M.B., c: Quaregnon M.B.) KB79 Evolution VRm% Sarnsbank M.B.

2'~------~------~------~------~------,------,------T 0.4

1.8 0.35 1.6 0.3 1.4 0.25 1.2 ~ :oR ....._E 'E 0.2 a: *-E > 0.8 a: 0.15 > 0.6 0.1

0.05

ol+-----~~------+------+------~.------.r------~.------+0 -350 -300 -250 -200 -1 50 -1 00 -50 0 a Time (myr)

KB79 Evolution VRm% Wasserfall M.B.

0.2

0.15

0.1

0.05

b

Fig. 19. Evolution ofVRm(%) (continuous line) and VRm(%)/myr (dashed line) in the deep well KB79 (a: Samsbank M.B., b: Wasserfall M.B.) KB81 Evolution VRm% Sarnsbank M.B. 3 ..------~------~------~------,------,------,------~ 0.4

0.35 2.5 0.3 2 0.25 S;, E 0.2 *---E a: 0.15 1 ·················· .. > 0.1

0.05

o+-----L-~------~------+------~~----~~------~------to -350 -300 -250 -200 -150 -1 00 -50 0 Time (myr) a

KB81 Evolution VRm% Wasserfa!! M.B.

2-r------;------,------;------,-----,------;------~0.4

1.8 0.35 1.6 0.3 1.4 0.25 S;, 1.2 ~ *& 1 ...... 02. 0~ E > 0.8 0.15 g; 0.6 ...... 0.1

0.05

-350 -300 -250 -200 -150 -1 00 -50 0 Time (myr) b

KB81 Evolution VRm% Quaregnon M.B.

2~.------~-----,------,------,------;-----~------~ 0.4

0.35

0.3 1.4 ..._ 0.25 >- 1.2 E E 0.2 -- *a: *-E > 0.8 a: 0.15 > 0.6 ··················· 0.1 0.4 ···················

0.2 ··················· 0.05 o+---_J~~----~------~-----+------~----~~-----T 0 -350 -300 -250 -200 -1 50 -1 00 -50 0 Time (myr) c

Fig. 20. Evolution ofVRm(%) (continuous line) and VRm(%)/myr (dashed line) in the deep well KB81 (a: Samsbank M.B., b: Wasserfall M.B., c: Quaregnon M.B.) KB86 Evolution VRm% Sarnsbank M.B. 2 0.4

1.8 0.35 1.6 0.3 1.4 0.25 1.2 >. cf. E ,o...... E 0.2 O' a: E > 0.8 a: 0.15 > 0.6 0.1 0.4 0.2 0.05

-350 -300 -250 -200 -1 50 -1 00 -50 0 Time (myr) a

KB86 Evolution VRm% Wasserfall M.B.

2.------~------~------~------~------~------~------.0.4

:.:1.4 •••••••••·•••••••••:·••••·•••••••••················'·················•••••·••F•••···••••••••••••••••:••••••··••••••••·····•·••i•·•••••·••••·••••••·•...... f······:··············""""""""""""!"""••••··· .. ········· ...... ,...... ,...... !...... •••··········••••••········ , ...... ::, 0.25 >. ~ 1.2 E ,o-.. 'E 0.2 O' a: E > a: 0.15 >

0.1

0.05

-350 -300 -250 -200 -1 50 -1 00 -50 0 Time (myr) b

Fig. 21. Evolution ofVRm(%) (continuous line) and VRm(%)/my:r (dashed line) in the deep well KB86 (a: Samsbank M.B., b: Wasserfall M.B.) KB185 Evolution VRm% Sarnsbank M.B.

2·~------~------~------~------,------,------,------r 0.4

1.8 0.35

0.3

0.25 5-...... E 0.2 0~ E a: 0.15 >

0.1

0.05 0.2 o+-----L_~~----4------4-----~~------~------~------+ 0 -350 -300 -250 -200 -150 -1 00 -50 0 Time (myr) a

KB185 Evolution VRm% Quaregnon M.B.

2.------~------~------~------~------~------~------r 0.4

0.35 : :: .::::::::::::.::::::::.:::r:::.::::::::::::::::::::::r::::::::::::::::::::::::::r:::::::::::::.::::::::::::r::::::::::::::::::::::::::r:::::::::::::::::::::::::r:::::::::::::::::::::::::: 0.3

0.25 5-. ,o o'- E -.0 E 0.2 a~--- a: E > 0.8 a: 0.15 > 0.6 0.1 0.4 0.2 0.05 0 0 -350 -300 -250 -200 -150 -100 -50 0 Time (myr) b

KB185 Evolution VRm% Eisden M.B. 0.4

0.35

0.3

0.25 5-...... E 0.2 0~ E a: 0.15 >

0.1

0.05

0 -300 -250 -200 -150 -1 00 -50 0 Time (myr) c fig. 22. Evolution ofVRm(%) (continuous line) and VRtn(%)/myr (dashed line) in the deep weii KB 18S (a: Samsbank M.B., b: Quaregnon M.B., c: Eisden M.B.) KBii5 Evolution VRm% Sarnsbank M.B.

3~------~------~------~------~------~------,------,-0.4

0.35 2.5 0.3 2 0.25 ~ E 0.2 of2.-- E a: 0.15 >

0.1 0.5 0.05

o+-----J_-r~------r------,_--~~~~----~------T------+0 -350 -300 -250 -200 -150 -1 00 -50 0 Time (myr) a

KBi i5 Evolution VRm% Quaregnon M.B. 2 0.4

1.8 0.35 1.6 0.3 1.4 0.25 1.2 1 ~ (j~ r:i 1:1: E E ~i 0.2 of!.-- a: E > ...... ·---~1- .··;- ···-··. --~ ...... ················ ·-········f••·······················-~---··· ...... : ...... a: 0.8 0.15 > 0.6 ...... 1\ ... i...... L ...... ! ...... ~ ...... 1...... l" ...... 0.1 0.4 ·············:·:····;················-··········:···············-···········:···························:··················-········:···························:··························· 0.05 0.2o+---·_ ... _... _... ,LJ_\TZ~ .. -~ ... _... _... _.... _ ... _ .. _... T.,._ ... _... _ ... _... _... _... _... _.. ·~·T_ ... _ ...~~-~- .. ~~;~....- ...~T_ ..._ ..._ ... _ ... _ ... _... _... _... ~.. j_· .. -_ ..._ ... _ ... _ ... _ ... _ ... _ ... ~ .. 1_ ... _... _...... _ _ ... _ ... _ ... _ ... ~ .. 0 -350 -300 -250 -200 -150 -1 00 -50 0 Time (myr) b

KBii5 Evolution VRm% Maurage M.B.

2,------~------~------~------~------~------~------. 0.4

0.35

0.3

0.25 ~ I ! l I -...E I 0~ i : :I ~ I I ~ 0.2 E a: 0.8 ...... '''"\'"'"''"''''"''""''"'~""''""'''"""'"'""\"'''"""''' '"""""l"""'"""'""""""'\'"''"'""''""""""·1··· ...... 0.15 >

0.1

0.05 :r lJ r :1~;_1: r: r -350 -300 -250 -200 -150 -1 00 -50 0 ° Time (myr) c

Fig. 23. Evolution ofVRm(%) (continuous line) and VRm(%)/myr (dashed line) in the deep well KB 115 (a: Samsbank M.B., b: Quaregnon M.B., c: Maurage M.B.) KB163 Evolution VRm% Sarnsbank M.B. 2 0.4 : : : : 1.8 0.35 1.6 0.3 1.4 0.25 ,o 1.2 >. o'- g ,o E 0.2 o'- a: E > 0.8 a: 0.15 > 0.6 0.1 0.4 0.2 0.05 o,4------L--~------4------r--~~~/--~·'--~------~------+------+o -350 -300 -250 -200 -150 -100 -50 0 Time (myr) a

KB163 Evolution VRm% Quaregnon M.B. 2 1.8 1.6 1.4

,o 1.2 o'- E a: 0.2 > 0.8 0.15 0.6 0.1 0.4 0.2 0.05

b

KBI63 Evolution VRm% Maurage M.B. 2 i ~ 0.4 1.8 0.35 1.6 1.4 0.3 .... 1.2 0.25 >- 0;,e E I I : g a: I ;: ~i ~I :~ I :0.2 ~ > : : : ! ! : E 0.8 ···················· ·····~···························~·························--f··············· ··········~···························~··························-~·-························· a: 0.15 > 0.6 0.4 0.1 0.2 r - ~--~~~ 1: i 0.05 0+-----~~~------~------~~~~--~------~------+------~lJ 0 -350 -300 -250 -200 -150 -100 -50 0 Time (myr) c

Fig. 24. Evolution ofVRm(%) (continuous line) and VRm(%)/myr (dashed line) in the deep well KB 163 (a: Samsbank M.B., b: Quaregnon M.B., c: Maurage M.B.) KBI74 Evolution VRm% Sarnsbank M.B. 0.4 2 1.8 i = 0.35 ~ ··················-~·········· !- : : ~ ] :I -- + .. ·············ur············u············r ... -- t + 1.6 0.3 1.4 ········t·························t:::::::::::::::::::::::::r:::::::::::::::::::::::::: 0.25 >. 1.2 ~ cf 0.2 cf E E a: a: > 0.8 0.15 >

0.6 0.1 0.4 0.05 0.2 0 ~3,~5-0____ JL __ 3~0~0------~25-0------~2~00~--~--~1t50~------71t00~-----:_5fro~------t0 Time (myr) a

KB174 Evolution VRm% Quaregnon M.B. 2 0.4

1.8 0.35 1.6 0.3 1.4 0.25 >. ,o 1.2 E C' ..... E 0.2 ??- a: E > a: 0.8 0.15 > 0.6 0.1 0.4 0.2 0.05

01+-----~-?------r------~--~--~~------r------r------r 0 -350 -300 -250 -200 -150 -1 00 -50 0 Time (myr) b

KB174 Evolution VRm% Eisden M.B.

2~------~------~------~------~------~~------,------, 0.4 1.8 0.35 1.6 0.3 1.4 0.25 ,o 1.2 >. ci' ~ E 0.2 oQ a: E > 0.8 ················-· a: 0.15 > 0.6 0.1 0.4 0.2 0.05 o. -350 ···H················••••••••:·••••••••••••••··•••••••r••••••••••••••·••••·•••~•••••••••••••••••••••·•l••••••••••••••••••••••••••:••••••••••••·•••·•··,-300 -250 -200 -1 50 -100 -50 0 Time (myr) c

Fig. 25. Evolution ofVRm(%) (continuous line) and VRm(%)/myr (dashed line) in the deep well KB 174 (a: Samsbank M.B., b: Quaregnon M.B., c: Eisden M.B.) KB186 Evolution VRM% Sarnsbank M.B.

2~------~------~------~------~------~~------~------~0.4 1.8 ···························i············-··-·-·········1 ...... 1...... --L--i···························i······--·------·--·········i···--·····---·----··--·---· 1 .6 ...... ·!--··--···----·----·--··----i------·--·i ...... l------!------l------0.35 1 .4 ...... :. ---~ ...... \ ...... \ ...... \...... \...... \...... 0.3 0.25 ;;:. E 0~ E 0.2 ?f.--- a: E > a: 0.15 >

0.1

0.05 0 : " : : ____...... : : 0 -350 -300 -250 -200 -150 -1 00 -50 0 Time (myr) a

KB186 Evolution VRM% Quaregnon M.B. 0.4

0.35

0.3

0.25 '-> ?f! 1 .2 E E 0.2 o'<--- a: E > 0.8 a: 0.15 > 0.6 0.1 0.4 0.2 0.05 o+-----~~------r------,_--~~~------~------~------~0 ·350 -300 -250 -200 -150 -100 -50 0 Time (myr) b

KB186 Evolution VRM% Maurage M.B.

2.------~------~------~------~------~------~------~ 0.4

0.35

0.3

0.25 ;;:. { 1 i I : II :: E I 0.2 ?f! ~ ~ ~ --- E a: 0.15 > g; ~.: ••····••••·•••••••••l···l••••••••••••••••·••·••••·•:•••••••••··•••••·•·•• !•••••·•·•·••••••·••·•••··•i•··•••···•·••••·•••••·••L•.••••····•••••••••l:•••••·•••••······••• 0.1 0.05 :.: ····················f\, .•.....•...•••....••.... !·························!···· 0 =:~~~!······················1·························•1••······················· 0 -350 -300 -250 -200 -150 -100 -50 0 Time (myr) c

Fig. 26. Evolution of VRm(%) (continuous line) and VRm(% )/myr (dashed line) in the deep well KB I8j (a: Samsbank M.B., b: Quaregnon M.B., c: Maurage M.B.) KBLOMMEL Evolution VRm% Sarnsbank M.B.

3.------~------~------~------~------~------~------r0.4

2 . Q$

~ ·: ••••••••••••••••••rr•••••••••••••• r: r r : : : : :::. ~ E 1.5 ~ ~ E > ~ 0.15 >

0.5

04.----~L-~------~------~--~~-;~------~------~~------~o -350 -300 -250 -200 -1 50 -100 -50 0 Time (myr) a

KBLOMMEL Evolution VRm% Quaregnon M.B.

2·~------~------~------~------~------~~------~------r0.4 ~ :: ::::::::::::.·::·:::::::::r:::::.:::::::::·:::::::::.::::::::::::::::.:::::::::r::::::::::::.:.:·:·:::::r.:::::::::::::::::::::::r::::::::::::::::::::::::r::::::::::::::::::::::::: :::5 1.4 l-0.25 >. .o 1.2 o' E ,o E o'-­ ~ E > 0.8 a: 0.15 > 0.6

0.4 0.2

01+------L~~------r------~--~---T~------r------~------+ -350 -300 -250 -200 -1 50 -1 00 -50 0 Time (myr) b

KBLOMMEL Evolution VRm% Base Neeroeteren Sst.

2.------~------~------,------,------;------;------r0-4

1.8 0.35 1.6 0.3 1.4

-!< 1.2 'E >~

Time (myr) c

Fig. 27. Evolution of VRm(%) (continuous line) and VRm(%)/myr (dashed line) m the deep well KBLommel (a: Samsbank M.B., b: Quaregnon M.B., c: Base Neeroeteren Sst.) KB168 Evolution VRm% Sarnsbank M.B. 3~------,------,------=------~------~------~------~ 0.4

2.5 0.35 0.3 2 0.25 s:. ?1< .._E fi 1.5 0.2 ?1< E > cc 0.15 >

0.1 0.5 0.05

?3·~.5~0----~-;~------~------~~--~~------~------~------~o -300 -250 -200 -150 -1 00 -50 0 Time (myr) a

KB168 Evolution VRm% Quaregnon M.B.

2,------~------~~------~------~------~------~------~0.4

1·8 ···························j···························:···························:···························j···························:···························:··························· 0.35 1.6 0.3 1.4 1.2 0.25 s:. ?1< E '· E 0.2 o« tr. E > cc 0.15 >

0.1

0.05

' I I I 0 -350 -300 -250 -200 -150 -1 00 -50 0 Time (myr) b

KB168 Evolution VRm% Base Neeroeteren Sst. 2 0.4

1.8 0.35 1.6 0.3 1.4 0.25 s:. ,a 1.2 a' .._E ,o E 0.2 o' tr. E > tr. 0.8 0.15 > 0.6 0.1 0.4 0.2 0.05 0 0 c

Fig. 28. E\·olution of VRm(%) (continuous line) and VRm(% )/myr (dashed line) in the deep well KB 16: (a: Samsbank M.B., b: Quaregnon M.B., c: Base Neeroeteren Sst.) KB172 Evolution VRm% Sarnsbank M.B. 3 0.4 0.35 2.5 0.3 2 0.25 >- E cf...... _ 0.2 0~ E 1.5 E a: a: > 0.15 >

0.1

0.5 0.05

o+---_J~~-----+------~~--~~-----j~-----t~-----t 0 -350 -300 -250 -200 -150 -100 -50 0 Time (myr) a

KB172 Evolution VRm% Quaregnon M.B. 0.4

0.35

0.3

0.25 ,o >- o' ...... _E ,_o E 0.2 O' a: E > 0.8 a: 0.15 > 0.6 0.1 0.4 0.2 0.05

O+------L~~~----+------~.~------~.~------+------~------~0 -350 -300 -250 -200 -150 -1 00 -50 0 Time (myr) b

KB172 Evolution VRm% Base Neeroeteren Sst.

2.------~------~------~------~------~------~------r 0.4

0.35 :: - i 1 i: I 1-:I _0.3 0.25 1 >- ...... _E ~ ·; ·:··.:::::::::::·::::J:r::::·:::.:::::::.::::.:::l:::::::::::::::::::.:::::::l:::::::::::::::::::::::::::l.::::::::::::::::::::::::::r:::::::::::::.::::::::::::j::::::::::::::::::::::::::: 0.2 cP. E > ...... §.. .\ ...... \ ...... : ...... L...... i...... L...... a: :.: ...... , ...... L...... !...... !...... ~ ...... !...... 0.15 > 0.1 0.4 0.05 0.2 1~:~::~~~~ I : l : 0+-----~~~------+-----~~~~~~~------~------~------+r 0 -350 -300 -250 -200 -1 50 -1 00 -50 0 Time (myr) c

fig. 29. Evolution ofVRm(%) (continuous line) and VRm(%)/myr (dashed line) in the deep well KB 172 (a: Samsbank M.B., b: Quaregnon M.B., c: Base Neeroeteren Sst.) KB201 Evolution VRm% Sarnsbank M.B. 3 0.4

0.35 2.5 0.3 2 0.25 >. ?f. -..E E 1.5 0.2 oP. a: E > a: ,. 0.15 > 0.1 0.5 0.05 o+-----JL· -~+h~----~'------~'~~~~----~-~~----~=~----~~=------to -350 -300 -250 -200 -150 -100 -50 0 Time (myr) a

KB201 Evolution VRm% Quaregnon M.B. 3 0.4

0.35 2.5 0.3 2 l: r i r r ' 0.25 ~ ,o O' E ,o E 1.5 0.2 O'-- a: E > a: 0.15 >

0.1 0.5 0.05

0 ?35o -3oo -2so -2oo -1 5o -1 oo -so 0 Time (myr) b

KB201 Evolution VRm% Maurage M.B. 3 0.4

0.35 2.5 0.3 2 0.25 ~ {jC::. J J : I : .§ II :: ,o E 1.5 I 0.2 o' a: E > a: : : 0.15 >

0.1 0.5 ·······~,j······ ·······•••••••••• ·••••••·•••••••••••••••••••E-~:~J••••••••••••••••••••••••••t•••••••••••••••••••••••••••.••••••••••••••••••••••••••• 0.05 o+-----~~------~------~~--~------7~----~~----~ 0 -350 -300 -250 -200 -150 -1 00 -50 0 Time (myr) c

Fig. 30. Evolution ofVRm(%) (continuous line) and VRm(%)/myr (dashed line) in the deep well KB20 I (a: Samsbank M.B., b: Quaregnon M.B., c: Maurage M.B.) KB198 Evolution VRm% Sarnsbank M.B. 4

3.5

3

2.5

0~ E 0.2 a: 2 > 1.5 0.15

0.1

0.5 0.05

Time (myr) a

KB198 Evolution VRm% Quaregnon M.B. 4 0.4

3.5 0.35 3 i [ L ! I i : ! 0.3 2.5 ...... ~ ...... ! ...... ) ...... ( ...... \...... l" ...... l ...... 0.25 >. 0~ E ~ 2 0.2 oR.-- 0: ~ > 0: 1.5 0.15 >

0.1

0.5 0.05

0 0 -350 -300 -250 -200 -150 -100 -50 0 Time (myr) b

KB198 Evolution VRm% Base Neeroeteren Sst. 4 0.4

3.5 0.35

3 0.3

2.5 0.25 >. cR ..._E E 2 ...... 0.2 ?fi a: E > a: 1.5 0.15 >

0.1 1:J ···· ..·.·· .. ·•.•.•.•.•.•.•.•_l:.1:l_1~:, ...... =_,. ~·=.... =.... =.... =.... =. = ...._ ...._ ...._ ...... !...... :...... :...... 0. _ , t·T- _ ::r:::~:;::::F<- _I __ · _ _ 0.05 0 -350 -300 -250 -200 -150 -1 00 -50 0 Time (myr) c

Fig. 31. Evolution ofVRm(%) (continuous line) and VRm(%)/myr (dashed line) in the deep well KB 198 (a: Samsbank M.B., b: Quaregnon M.B., c: Base Neeroeteren Sst.)