Abstract

The Eocene lacustrine sediments are the primary source rocks in the Huizhou Sag of the Pear River Mouth Basin. This study employs basin modeling for four representative wells and two profiles in the Huizhou Sag to reconstruct the process of generation, expulsion, migration, and accumulation of hydrocarbon fluids. The Eocene source rocks started to generate hydrocarbon at 33.9 Ma and are currently in a mid-mature and postmature stage. Hydrocarbons are mainly expelled from the Eocene Wenchang Fm, and the contribution of the Eocene Enping formation is minor. Under the driving forces of buoyancy and excess pressure, major hydrocarbons sourced from the Eocene source rocks firstly migrated laterally to the adjacent Eocene reservoirs during the postrift stage, then vertically via faults to Oligo-Miocene carrier beds, and finally laterally to the structural highs over a long distance during the Pliocene-Quaternary Neotectonic stage, which is controlled by both structural morphology and heterogeneity of carrier beds. Fault is the most important conduit for hydrocarbon fluid migration during the Neotectonic stage. Reactivation of previous faults and new-formed faults caused by the Dongsha Movement (9.8–4.4 Ma) served as vertical migration pathways after 10.0 Ma, which significantly influenced the timing of hydrocarbon accumulation in the postrift traps.

1. Introduction

As the largest basin in the northern continental margin of the South China Sea (SCS), the Pearl River Mouth Basin (PRMB) has become the exploration focus due to a series of hydrocarbon discoveries. Hydrocarbon exploration in the PRMB began in 1973, and reservoirs were successively found in several hydrocarbon generation sags. According to the existing commercial discoveries, the northern basins located in the continental shelf mainly produce oil, and the southern basins located in the continental slope are gas producing basins. The Huizhou Sag (HZS), located in the north of the PRMB, is a typical rift basin during the Cenozoic. The recent explorations show that about 64% of petroleum reserves in the PRMB were found in intradepression highs on both northern and southern sides of the HZS [1]. In the 1990s, the petroleum exploration is mainly focused on the Oligocene postrift Zhuhai transitional formation (Fm) and overlying strata. Recently, there were more discoveries found in the Eocene syn-rift Wenchang and Enping Fm. Commercial oil flow was obtained in the Enping Fm exampled by well 254, and oil was found in the Wenchang Fm exampled by well 254 and well 81. In 2011, test results of the well 257 revealed oil in the Eocene tight sandstone reservoir with more than 4000 m depth in the HZS, and the oil flow of well 257 is up to 173 m3 per day during production testing stage [2]. Three years later, tight oil was discovered in the LF13 fractured anticlinal structure belt in the Huilu area, and the oil flow of well 141 and well 81 is up to 210 m3 per day during production testing stage [3]. The geological reserves of the new tight sandstone oil discovery are more than  m3 in the well 257, well 144, and well 81 (according to the reserves report of Shenzhen Company of CNOOC China Ltd.). The Eocene Wenchang and Enping Fm, developed in the lacustrine and swamp environment during the midcontinental rift stage, are generally considered to be two main sets of source rocks in the PRMB; the new discoveries not only indicated the huge hydrocarbon potential of the deep Eocene strata but also brought the new question about the petroleum migration of the Eocene source rocks in the HZS.

The driving forces for secondary petroleum migration are buoyancy and groundwater flow, and the restraining force is capillary pressure, which increases with decreasing pore-throat size, increasing interfacial tension and wettability [47]. However, the nature of petroleum migration pathway is still controversial. Some researchers [726] believed that the petroleum fluids migration is mainly driven by structural morphology, and modeling of petroleum fluids migration pathways can be a powerful tool to reduce exploration risk. Others [2731] propose that petroleum fluids tend to pass through a high-permeable, thin sandstone bed than through a low-permeable and thick sandstone bed, controlled by the heterogeneity of the porosity and permeability of the carrier beds.

Because of incomplete geological data and few drilling wells, the mechanism of secondary hydrocarbon fluids migration in the PRMB is still unclear. Long et al. [32] believed that hydrocarbon migration was mainly controlled by the heterogeneity of carrier beds, while Shi (2013) concluded that structural morphology in the HZS was the primary control of petroleum migration pathway based on 3D computer modeling. The migration process and dominant pathway in the HZS are still poorly understood. In this paper, four wells and two regional sections are selected to analyze the process of hydrocarbon generation, fluids migration, and accumulation in the HZS of the PRMB.

2. Geological Setting

2.1. Tectonic Evolution and Sedimentary Character

The PRMB is located in the northern continental margin of the SCS, which covered an area of approximately  km2. It consists of four depression zones, including the Zhu-I Depression, the Zhu-II Depression, the Zhu-III Depression, and the Chaoshan Depression, and three uplift zones including the North Uplift Zone, the Central Uplift Zone, and the Southern Uplift Zone (Figure 1(a)). The Zhu-I Depression has typical double-layer geologic structure of rift overlain by sag and was controlled by a series of major boundary faults [33]. The HZS is located in the central part of the Zhu-I Depression and characterized by half grabens or strongly asymmetric grabens (Figure 1 profile AA′ and profile BB′). The geological evolution of the PRMB is significantly influenced by the Philippine Sea plate, the Eurasian plate, and the Indo-Australian plate [3436]. The HZS mainly has experienced three tectonic events during the Cenozoic: the Paleogene continental rift stage, the Miocene passive continental margin stage (postrift stage), and the Late Miocene-Quaternary Neotectonic stage (Figure 2). During the continental rift stage, the HZS developed the Paleocene Shenhu Fm, the Eocene Wenchang, Enping Fm, and the Oligocene Zhuhai Fm. The Eocene Wenchang Fm, the most important source rock, consists of dark lacustrine mudstone with gray fluvial sandstone and coal interbedded. The Eocene Enping Fm, minor source rock, is mainly composed of lacustrine-swamp mudstone, coal bed, and fluvial sandstone. The Oligocene transitional Zhuhai Fm is composed of light gray delta-littoral siltstone, which is the major reservoir in the HZS. During the passive continental margin stage (postrift stage), the HZS developed the Miocene Zhujiang and Hanjiang Fm. The Miocene formations of Zhujiang and Hanjiang consist mainly of littoral siltstone, neritic mudstone, marine mudstone, and carbonate. The widespread and muddy upper Zhujiang Fm acts as a regional cap rock. The Miocene Yuehai, Pliocene Wanshan, and the Quaternary Fm were deposited in a marine environment during Neotectonic stage. From the Palaeocene to the Neogene, there are three intervals of combinations of oil source rock, reservoir, and cap rock (Figure 2).

Due to regional extensional stresses, the HZS was a lacustrine basin during the sedimentary period of the Eocene Wenchang Fm. The provenance is mainly from the northwest and northeast direction and shoreline-shallow lacustrine and deep lacustrine facies developed in the Huixi half-graben, Huibei half-graben, Huizhong low-uplift, Huinan half-graben, and Huilu low-uplift. Under the regional extensional environment, high quality lacustrine source rocks were widely developed during the sedimentary period of the Wenchang Fm (Figure 3(a)). After that, the ancient lake basin shrank gradually during the sedimentary period of the Eocene Enping Fm, and a large amount of fluvial deposits came from the southwest, northwest, and northeast direction. The shoreline-shallow lacustrine mudstone and swamp coal-bearing strata developed in the Huixi half-graben, Huibei half-graben, and Huinan half-graben during the sedimentary period of the transgression systems tract of the Enping Fm (Figure 3(b)), and the fluvial-delta developed in the most areas during the sedimentary period of the highstand systems tract of the Enping Fm (Figure 3(c)).

2.2. Source Rock
2.2.1. Distribution of Source Rock

The lowstand and highstand systems tracts of the Wenchang Fm developed in the local areas, while the transgression systems tracts of the Wenchang Fm were deposited over a wide area in the HZS (Figure 3(a)). During the extensional rifting stage, shoreline-shallow lacustrine and deep lacustrine mudstones were formed. The thickness of lacustrine mudstone ranges from 200 m to 1000 m, and more than 750 m source rocks developed in the Huixi half-graben. The percentage of deep lacustrine mudstone in the well 132 reaches up to 82.5%, which is much more than that of shoreline-shallow lacustrine mudstone in the other wells (Table 1).

The lowstand systems tract of the Enping Fm developed in the limited area and mainly consists of fluvial-delta facies with lower percentage of source rocks. The transgression and highstand systems tracts of the Enping Fm are widely distributed in the HZS (Figures 3(b) and 3(c)). The transgression systems tract of the Enping Fm mainly consists of shoreline-shallow lacustrine mudstone and swamp carbonaceous mudstone. The highstand systems tract of the Enping Fm mainly consists of shoreline-shallow lacustrine mudstone, swamp carbonaceous mudstone, and fluvial-delta mudstone. The percentage of mudstone in the Enping Fm is generally from 15%~40% (Table 1).

2.2.2. Kerogen Type

The composition and structure of kerogen depend on the origin of the organic matter from which it has evolved, as well as its degree of thermal evolution. In this study, total organic carbon (TOC) content was measured by using a LECO CS-400 analyzer, and Rock-Eval pyrolysis and organic petrology were used to identify the kerogen type of organic matter. Rock-Eval pyrolysis was performed by using a Rock-Eval II instrument. All the samples were heated to 600°C in a helium atmosphere to obtain the primary parameters, such as hydrocarbon generation () at a low temperature peak within a range of 150~300°C, hydrocarbon generation () at a high temperature peak in the range of 420~550°C, carbon dioxide (), and the temperature at maximum rate of hydrocarbons generation () [39]. Because the oxygen index was not measured by Rock-Eval II instrument, the cross-plot of HI versus was used to divide the kerogen type. According the pyrolysis experimental data of some wells in the HZS, the kerogen of the Wenchang source rock is mainly Type-II1 and Type-I, and Enping source rocks are generally composed of Type-II1 and Type-II2 kerogen (Figure 4).

Organic petrography of the Eocene source rocks was also applied to characterize and identify the various organic constituent. Microscopic examination of the Eocene source rocks from 6 wells in the HZS was performed under transmission light. Various types of organic materials were recognized (Figure 5(a)): woody organic matter (W), coaly organic matter (I), exinite (E), cutinite (C), algal organic matter, and amorphous organic matter (AOM). Microscopic examination shows that the kerogens in the Eocene Wenchang source rocks are mainly classified into three types (Figure 5(b)). Type-I is dominated by amorphous organic matter, for example, algae, which was formed in the deep lacustrine environment with high productivity. Type-II is dominated by mixed organic matter rich in hydrogen maceral exinite, and Type-III is dominated by woody vitrinite, which were mainly formed in shore-shallow lacustrine and fluvial environment. Microscopic examination of 6 wells shows that the organic matters of the Enping source rocks were mainly composed of terrestrial phytodetritus, formed in the oxygen-rich environment with high hydrodynamic force (Figure 5(b)). Organic maceral and submaceral characteristics of the Eocene source rocks in the HZS have been studied by means of composite optics of organic petrology by Zhu et al. [40]. The results show that relatively abundant exinite rich in hydrogen results in oil-prone Types-I-II1 of organic matter in deep lacustrine facies, while relatively abundant vitrinite results in gas and condensate-prone Types-III-II2 of organic matter in shore-shallow lacustrine and swamp facies.

All these data indicated that the Eocene Wenchang source rocks are of Type-I/II/III mixed kerogens and the Eocene Enping source rocks are dominated by Type-II2 kerogen with a small amount of Type-III. In the HZS, Type-III kerogen is generally poor to fair source rocks with limited hydrocarbon generating potential. Type-II2 kerogen is fair to good source rock, which contributed a lot to hydrocarbon generation, and Type-II1 and Type-I kerogen are good to excellent source rocks, which made significant contribution to hydrocarbon generation in the HZS (Figure 6).

2.2.3. Maturity

Vitrinite reflectance (Ro) was measured by an oil immersion lens and a Leixa MPV Compact II reflected light microscope fitted with a microphotometer. Measured Ro values in six wells are primarily between 0.5% and 1.2%, indicating that the Eocene source rocks have entered the mature generation stage (Figure 7). of the Eocene source rocks are mainly in range of 430°C to 460°C (Figure 8), which suggests that most source rocks in the HZS are sufficiently mature to have generated oil [41, 42]. Because samples from deeper sags are not available in the study area, these measured Ro data can not reflect the whole maturity of the Eocene source rocks. The theoretical maturity of deeper interval was calculated by using 1D burial and thermal history simulations. According to the previous research [3], the Eocene Wenchang source rock in the HZS is mature within Ro range of 0.8%~1.8%, and the Enping source rock is in range of 0.6% to 1.4%.

2.2.4. Geochemical Characteristics

Table 2 lists the geochemistry characteristics of source rocks with different facies in the HZS. The Wenchang source rocks consist of various facies from swamp to deep lacustrine, with HI ranging from 83.5 mg/g·TOC to 613.1 mg/g·TOC. The TOC content of the Wenchang source rock ranges from 0.5% to 7.7%, with mean value of 2.0%. The highest TOC value of 7.7% is from the deep lacustrine dark mudstones of the well 132 (Figure 9). The Enping source rocks consist of swamp mudstone, swamp carbonaceous mudstone, and fluvial-delta mudstone. The fluvial-delta mudstone belongs to medium quality hydrocarbon source rocks, with TOC value ranging from 0.3% to 11.9%, of which hydrocarbon generation potentials () ranged from 0.2 mg/g to 32.8 mg/g. The swamp-lacustrine mudstone belongs to medium-high quality hydrocarbon source rocks, with TOC value ranging from 0.6% to 6.9%, of which hydrocarbon generation potentials () ranged from 0.6 mg/g to 27.1 mg/g. Because of the limited number and distribution of the source rock samples, the lateral variation of TOC in deeper undrilled part of the HZS was predicted by using the logging-seismic prediction technology. The relationship between measured TOC and well-logging values of sonic transit time, natural gamma, resistivity, density, and neutron was established to obtain the prediction models (Figure 9).

The ratio of pristane (Pr) and phytane (Ph) value is 6.35 in the Wenchang source rocks from well 243 with low content of 4-methyl sterane (Figure 10), which means that the source rocks were formed in a shallow water environment. The ratio of Pr/Ph value is 2.70 in the Wenchang source rocks from well 132 with high content of 4-methyl sterane, which indicates that the source rock is originated from lacustrine with algal input ([43, 44], Fu and Zhu, 2007). The ratio of Pr/Ph values is more than 6 in the Enping source rocks. The 4-methyl sterane content is low and the content of C29 sterane is relatively higher than the levels of C27 sterane and C28 sterane (Figure 10), which indicates that the Enping Fm developed in the environment containing both algal and land plant material ([43, 44], Fu and Zhu 2007).

2.3. Reservoir and Caprock

The Eocene sandstone and andesite of Wenchang Fm reservoir are the main reservoir [46]. Percentage of sandstone in the Wenchang Fm is 4.86% and 37.34%, respectively. Porosities of the Eocene Wenchang Fm from well 254 range from 4.9% to 25.3% (Figure 11(a)). The Enping reservoir is fluvial sandstone, with sandstone content of single well ranging from 24.13% to 71.19%. Porosities of the Enping Fm range from 9.0% to 20.0%, and permeabilities of the Enping Fm range from 0.1 to 28.4 mD (Figure 11(a)). Porosities of the Oligocene Zhuhai Fm range from 4.6% to 15.2% with average value of 9.0%, and permeabilities range from 0.004 mD to 223.0 mD. The sandstones of the Zhujiang Fm show considerable variation in porosity and permeability, with the measured porosities ranging from 5% to over 30% (Figure 11(b)) and the measured permeabilities ranging from 0.1 to over 1000 mD (Figure 11(c)). The average porosity and permeability of the upper Zhujiang Fm are 17.6% and 750 mD, respectively, and the average porosity and permeability of the lower Zhujiang Fm are 18.7% and 587 mD, respectively [47]. Overall, the Wenchang and Zhujiang Fm are strongly heterogeneous.

The Hanjiang, Zhujiang, and Wenchang Fm are the three regional seals in the HZS, which are characterized by wide distribution and thick cumulative thickness (Table 3). The lower Zhujiang Fm together with Enping Fm is the regional seals and Zhuhai Fm acted as the local seal.

3. Modeling Method and Parameters

3.1. Method

In this study, 1D and 2D models were reconstructed using the petroleum modeling software BasinMod 2012. Four representative wells from the Huixi half-graben and the Huilu low-uplift (Figure 1(a)) were selected to perform 1D simulation by BasinMod 1D module in order to better understand the burial and thermal history in the rift and postrift setting and reveal the timing of hydrocarbon generation and expulsion. The calculated and calibrated results of 1D modeling were then input into 2D models [48]. As most oil fields are concentrated in the Huixi half-graben and the Huilu low-uplift, two seismic lines (Figure 1(a)) covering two areas with available measured data were selected to reconstruct the process of petroleum generation, fluid migration, and accumulation by BasinMod 2D module. With the available wells and seismic data, the hydrocarbon generation could be estimated quantitatively, while the process of migration and accumulation has been studied qualitatively to semiquantitatively during the rift and postrift stage. The 3D modeling of fluid flow will not be discussed in this paper due to commercial confidentiality.

The reconstruction of burial history was based on geological information, tectonostratigraphic history of the region, and drilling and seismic data. Exponential compaction correlation method [49] was employed in this paper to reconstruct the burial history of strata in the HZS. Porosity is calculated by the following equation: where is porosity (%), is the initial surface porosity (%), is the compaction factor for each lithology (cm−1), and is depth (m). Parameters for various lithologies under normal pressure are shown in Table 4.

Permeability was calculated by the following equation [50]:where is permeability (mD), is porosity (%), and is the specific surface area of the rock (m2).

Maturity and hydrocarbon generation were simulated by using the LLNL EASY% Ro model [51, 52], which was calibrated by the measured vitrinite reflectance (Ro%) and bottom hole temperature data. Expulsion calculation followed the saturation expulsion method [53] which assumed that the hydrocarbon will begin to be expelled from the source rocks when saturation of hydrocarbons reached a certain threshold related to the rock types and properties [54]. The irreducible water saturation was defined as 1.0.

Hydrocarbon migration was simulated by the transient migration model, which is based on the state equation, the mass conservation law, and the movement equation. The relative permeabilities of gas, oil, and water are calculated from the saturation-permeability relationship, and capillary pressure of gas and oil is calculated from saturation. Two main assumptions of migration in BasinMod are as follows: There is no mass transfer between oil and gas phases during migration in one time step. Migration of gas is much faster than that of oil and water, which occurs before the migration of water and oil. The migration of three-phase (oil, gas, and water) is accomplished by a two-step and two-fluid (gas-liquid and oil-water) migration calculation in BasinMod. The fluid flow equations ([55], Aziz & Settari, 1986, Freeze & Cherry, 1979) during the first migration step are written as follows:

The fluid flow equations (Aziz & Settari, 1986, Freeze & Cherry, 1979) during the second migration step are written as follows:

The whole migration step was calculated with closed boundary ():

And initial condition is as follows:where (gas, liquid, oil, and water phases, resp.), is the definition of thickness in -direction, is elevation (, positive downwards), is -direction transmissibility (), is -direction transmissibility (), is the relative permeability of phase in -direction, is the relative permeability of phase in -direction, is pressure of phase (MPa), is the capillary pressure (MPa), is the density in terms of pressure/distance (), is the density of phase (kg/m3), is the gravitational acceleration (m/s2), is the conversion constant (32.2 /·ft/sec2), is the saturation of phase (%), and is the formation volume factors defined by where is the volume occupied by a fixed mass of component () at the reservoir condition and is the volume occupied by the same component () at standard conditions.

Density of generated oil and gas is calculated by the hydrocarbon kinetics module. During the migration process, the density of fluids is calculated by pressure and temperature as follows:where is the density at standard surface conditions (kg/m3), α is constant (), and β is constant ( MPa−1).

Density of water can be used as surface water density or a calculated water density under calculated temperature and pressure in the BasinMod software. The default value of surface water density is 1.3 gm/cm3. The equation used for calculating water density considering temperature and pressure is written as follows:where is the calculated water density (kg/m3), is the initial water density (kg/m3), (Pa), and is  (K).

The viscosity of fluids for the modeling is obtained from an empirical formula fitted to viscosity versus temperature data as follows (after CRC Handbook of Chem. & Physics, 1986):where is temperature (°C) and is viscosity (Pa·s).

Fluid flow equations together with sediment compaction are written below. The measured pressure data were used in this study as calibration data.where is void ratio, is the initial void ratio of Fraction , is the initial void ratio of Fraction , is frame, or matrix, pressure (Pa), is initial frame, or matrix, pressure (Pa), is the exponential compaction factor, is the linear compaction factor, and Fraction is the portion of lithology which compacts exponentially, that is, all lithologies except sands. Fraction is the portion of lithology which compacts linearly, that is, all sands.

3.2. Modeling Parameters
3.2.1. Formation Tops and Lithologies

The formation tops and lithologies were obtained from drilling results and seismic data. The lithology of each stratum (Table 5) was defined by mixing pure lithologies with specific petrophysical parameters and corrected according to the published work [3, 38, 56]. The functions “Facies Palette,” “Thin Bed,” and “Lenses and Mounds” were used to capture the laterally and vertically lithological variation and the heterogeneity of reservoir (Figure 12).

3.2.2. Tectonic Event

Several tectonic events have been recognized in the HZS (Figure 2). During the Paleogene continental rift stage, the Zhu-Qiong I (49.0 Ma) and Zhu-Qiong II Movement (38.0 Ma) not only formed the structural configuration of the HZS [43, 44, 57] but also developed regional unconformities (Figure 13(f)). In the early Oligocene, the Nanhai Movement (33.9 Ma) shifted the HZS from rifting to steady thermal subsidence (Figure 13(e)). These three movements developed most important unconformities in the HZS, which caused previous lacustrine source rocks to be eroded in the uplift area. The HZS entered the marine sedimentary environment after the Baiyun Movement (23.0 Ma), and fault activity decreased significantly (Figures 13(c) and 13(d)). During the Late Miocene-Quaternary Neotectonic stage, Dongsha Movement (10.0 Ma) led old faults to be reactivated and many new west-northwest-trending faults developed (Figure 13(b)). The most significant erosion events are related to the SCS opening process after the deposition of Eocene Enping Fm.

As regional tectonic events, the most significant unconformities formed at the end of the sedimentary period of the Wenchang and Enping Fm. As available measured Ro data is insufficient, the eroded thickness was estimated by extrapolation method of stratum trend (Figure 14). Erosion of the Wenchang Fm was mainly distributed in four areas, which are the Huixi and Huibei half-graben with erosion thickness of 600~1200 m, Huidong half-graben with erosion thickness of 500~1200 m, and Luxi half-graben with erosion thickness of 200~600 m [58]. Erosion of the Enping Fm was thinner than that of the Wenchang Fm, ranging from 100~400 m [58]. For the structure, only those faults with significant throw were included into the models in order to simplify the geology for saving the computation time.

3.2.3. Kinetic Models

Li et al. [60] suggested that the activation energy of the Wenchang dark mudstone primarily ranges from 46 to 63 kcal/mole, and their hydrocarbon generation potentials have a negatively skewed distribution. Jiang et al. [61] analyzed the chemical kinetic equations of the Enping dark mudstone and carbonaceous mudstone, which showed that dark mudstone has an activation energy of 49.5~64 kcal/mole with a normal distribution of the hydrocarbon generation potential, and carbonaceous mudstone has an activation energy of 46~66 kcal/mole, showing a skewed distribution of the hydrocarbon generation potential. These results of our previous papers were adopted in this paper.

3.2.4. Geothermal Field and Boundary Conditions

Because present water depth is less than 200 m in the HZS, the variation of paleobathymetry has slight influence on the thermal evolution. The paleobathymetry is estimated according to the sedimentary facies characteristics (Table 6) in this paper and published literature [3].

The bottom hole temperature (BHT), surface temperature, and present heat flow were collected from boreholes and corrected according to Yuan et al. [62]. According to the recent researches [6267], the current geothermal gradient of the HZS ranges from 30.0°C/km to 36.0°C/km and increases from north to south.

Present-day basal heat flow is calculated from thermal conductivities of the rock units and geothermal gradients which is determined by measured BHT according the following equations:where is the average conductivity of the sediment column (W/m·°C), is the thickness of each lithology (m), and is the conductivity of each lithology (W/m·°C).where is the basal heat flow at the bottom of the sedimentary pile (mW/m2), is the temperature at the surface of the sediment column (°C), is the temperature at the base of the sediment column (°C), and is the sedimentary thickness from the surface to the base of the stratigraphic column (m).

The present-day basal heat flow of the HZS ranges from 51.2 mW/m2 to 75.33 mW/m2, with average of 62.98 mW/m2 (Table 7). The paleoheat flow was applied in this study based on previous work [6870]. Thermal evolution is calibrated by measured Ro and temperature values. The HZS shares similar thermal history due to the same rifting events. When the calculated basal heat flow value (54.5 mW/m2) from well 131 was used as an input parameter in steady-state heat flow model to calculate the thermal history, the simulated maturity (Ro%) values and temperatures are much higher than the measured Ro values and measured BHT (Figure 15(a)). When the lower basal heat flow value (49.0 mW/m2) was inputted into the same geologic model again, the new “simulated” Ro values fit the “measured” Ro values reasonably well, while the “simulated” temperature was much lower than measured BTH below 3500 meters (Figure 15(b)). After testing the steady-state heat flow method with constant heat flow of 54.5 and 49 mW/m2 and a variable paleoheat flow (Figure 15(c)), the best fit between measured and calculated Ro and temperature was achieved by setting a variable paleoheat flow increasing from background value of 55 mW/m2 [71, 72] to 73.84 mW/m2 during the continental rift stage and then decreasing to the present heat value during the postrift stage (Figure 16). Thus, the steady heat flow method with was a variable paleoheat flow used to calculate geothermal history.

3.2.5. Pressure Field

Based on measured formation pressure, present pressure is characterized by normal pressure within most of the HZS (Figure 17(a)). The pressure of most drilled wells belongs to the range of normal pressure zone, with pressure coefficients ranging from 0.96 to 1.06 (Figure 17(b)). There is no obvious overpressure in the well 132, well 81, and well 241 (Figures 17(c), 17(d), and 17(f)), and formation pressure is equal to the hydrostatic of normal pressure zone. As the buried depth increased, formation pressure curve of well 131 gradually deviates from the hydrostatic pressure, and slight overpressure develops in the lower Zhujiang Fm, Zhuhai Fm, and Enping Fm (Figure 17(e)).

4. Results and Discussion

Two 2D models were constructed based on seismic data and drilling well data. The variation of paleoheat flow calibrated from 1D modeling was adopted in the 2D model to simulate the process of migration and accumulation of hydrocarbon fluids in the HZS [48].

4.1. Maturity History, Petroleum Generation, and Expulsion
4.1.1. Burial and Thermal History

Since the selected wells are situated in the structural high, the mature level of the Eocene source rocks is generally in the early-main nature stage. The well 81 with total depth of 4583.0 m, located in the northern part of the HZS, has drilled into the Eocene Enping Fm. The drilling data were used to establish geologic model for studying the thermal history (Figure 18(a)). The Eocene Enping source rocks are currently in a main mature stage. After the rapid subsidence during the Eocene continental rift stage, the Eocene Enping source rocks entered the early mature stage at 18.5 Ma. With the increase of buried depth, the Eocene Enping Fm further evolved into the main mature stage, with the EASY % Ro ranging from 0.7% to 1.0% during the period from 9.5 Ma to the present. The bottom of the Miocene Zhujiang and Oligocene Zhuhai Fm is currently in an early mature stage.

The well 131 with a total depth of 4804.0 m, located in the central part of the HZS, has drilled into the Eocene Enping Fm. The lithological data obtained from drilling well were used to establish geologic model, from which a better fit between measured and calculated Ro was obtained (Figure 18(b)). The burial and thermal history is quite similar to that of the well 81. The Eocene Enping source rocks are currently in an early-main mature stage, with Ro values ranging from 0.6% to 0.9%. The Eocene Enping Fm evolved into the main mature stage after 9.5 Ma. The upper part of the Eocene Enping source rocks is currently in an early mature stage, and the lower part is in a main mature stage. The lower part of the Miocene Zhujiang and Oligocene Zhuhai Fm is also currently in an early mature stage.

The well 241 with a total depth of 3853.1 m, located in the southern margin of the HZS, has drilled into the Pre-Cenozoic granite. Because this well is located in the uplift area, approximately 1 km of the Eocene Wenchang Fm and whole Enping Fm were eroded due to the Zhu-Qiong II and Nanhai Movement (Figure 18(c)). There is only 165 m thickness of the Eocene Wenchang Fm preserved. The Eocene Wenchang source rocks are currently in an early mature stage, with Ro values ranging from 0.45% to 0.75%, which entered the early mature stage at 12.0 Ma. The lower part of the Miocene Zhujiang and Oligocene Zhuhai Fm is also currently in an early mature stage.

The water depth of well 132 is deeper than that of above three wells, which is located in the Huilu low-uplift. The oldest layer is the Eocene Wenchang Fm with a total drilling depth of 3280.0 m. The thickness of the Eocene Enping and Wenchang Fm is estimated from drilling data. Because of regional tectonic uplift movement, about 400 m of the Eocene Wenchang Fm and 100 m of the Eocene Enping Fm are eroded. It is shown from the simulated results (Figure 18(d)) that the Eocene Wenchang and Enping Fm are currently in a marginally mature stage, with Ro values ranging from 0.5% to 0.6%, and the overlying strata are all in the immature stage (Ro < 0.5%). The Eocene Wenchang and Enping source rocks evolved into the early mature stage at 9.0 Ma.

The Eocene source rocks in the deeper sag are generally in the late-mature stage, with Ro value more than 2.0% (Figures 19(a) and 19(b)). The maturity level in the Huilu low-uplift (Figure 19(c)) is lower than that in the Huixi half-graben (Figure 19(c)). Most of the generated hydrocarbon is from the Eocene Wenchang and Enping source rocks in the deep part of the HZS.

4.1.2. Timing of Hydrocarbon Generation and Expulsion

The hydrocarbon generation process of the Eocene source rocks was simulated using four wells as examples. “Oil” and “gas” in this paper refer to the composition of generated fluids/hydrocarbons. The Eocene Enping source rock of well 81 is mainly in hydrocarbon generation stage, with oil and gas generation amount of 78 mg/g·TOC and 118 mg/g·TOC, respectively (Figure 20(a)). The first stage (approximately 36.0–18.5 Ma) is a minor phase of hydrocarbon generation without any hydrocarbon expulsion. The second stage (from 18.5 Ma to present) is the phase of major hydrocarbon generation and expulsion. The hydrocarbon expulsion began at 18.5 Ma, expelled little in the second stage because of low maturity level of the source rocks. Because of lower maturity level of the Eocene Enping source rocks, the hydrocarbon generation amount is relatively low in well 131, with present oil and gas generation amount of 16.8 mg/g·TOC and 29 mg/g·TOC, respectively (Figure 20(b)). The hydrocarbon expulsion also began at 18.5 Ma but was expelled little. The present hydrocarbon generation amount of well 241 is less than 30 mg/g·TOC, and hydrocarbon expulsion began at 14.0 Ma but was rarely expelled (Figure 20(c)). In general, the hydrocarbon generation and expulsion of the Eocene source rocks on structural high are not high, which is because of the location of structural high and low maturity level.

The situation in the deep part of the HZS is different. The period of large hydrocarbon generation of the Enping source rocks in the Huibei half-graben is generally in the thermal subsidence stage and mainly occurred within the last 20 Ma (Figures 20(e) and 20(f)). Hydrocarbon expulsion of the Enping Fm began at 33.9 Ma but was rarely expelled. The period of large hydrocarbon generation of the Wenchang source rocks in the Huixi half-graben is earlier than that in the Huibei half-graben, which started mainly at 33.9 Ma (Figures 20(g) and 20(h)). The modeling results indicate that hydrocarbon generation potential of the Wenchang Fm in the Huixi half-graben is much higher than that in the Huibei half-graben, with oil generation potential of 70 mg/g·TOC and gas generation potential of around 150 mg/g·TOC. Hydrocarbon expulsion of the Wenchang Fm began at 33.9 Ma and was expelled much more (Figures 20(g) and 20(h)). The oil expelled volume reaches 11 mg/g rock in the center of the Huixi half-graben (Figure 20(h)).

4.2. Petroleum Migration and Accumulation Modeling

In this study, short and long distance migration of the Eocene petroleum have been predicted for the hydrocarbon accumulation in the HZS. Most hydrocarbons generated from the deep Eocene source rocks of the HZS started to generate from 33.9 Ma and shallower source rocks are still generating hydrocarbons. Previous research indicated that most petroleum of oil field originated from the deep lacustrine Wenchang dark mudstone and shoreline-shallow lacustrine Enping dark mudstone and accumulated via lateral migration from the deeply buried source rocks in the HZS [2, 3, 43, 44, 47]. The modeling results indicate that oil expulsion in the HZS is mainly from the Eocene Wenchang source rocks and occurred during postrift stage. The shaly facies of the upper Zhujiang Fm and overlying strata acted as an effective seal rock for the lower Zhujiang and Zhuhai sandy reservoirs, limiting the vertical migration of hydrocarbons from these reservoirs towards higher stratigraphic intervals.

In the eastern HZS, petroleum mainly accumulated not only in structural highs of the shallow Zhujiang and Zhuhai Fm, but also in the deep Wenchang Fm (Figure 21(a)). Petroleum accumulated in the Zhujiang and Zhuhai Fm generally started at 10 Ma, with different oil accumulated volume in different area. In the Huilu low-uplift, oil accumulated volume reached 0.08 m3/m3 at present day (Figures 21(c) and 21(e)). Petroleum accumulated in the southwest of the Luxi half-graben is relatively high (Figures 21(d) and 21(f)), and oil accumulated volume exceeded 0.10 m3/m3 at 5.3 Ma (Figure 21(d)). Petroleum generated from the Eocene Wenchang source rocks began to accumulate in the adjacent Wenchang reservoirs from 16 Ma. The maximum accumulated oil volume in the Luxi half-graben and the Huinan half-graben is both around 0.055 m3/m3 at 5.3 Ma (Figures 21(g) and 21(h)). LF13-2 oil field (Figure 21(i)) was successfully discovered with proven geological reserves of  m3 [38]. Hydrocarbon sourced from the Wenchang Fm migrated from north area via the Zhuhai and Zhujiang carrier and faults to the LF13-2 structural high (Figure 21(j)).

In the Huibei half-graben, hydrocarbons mainly accumulate on the top of the Wenchang and Enping Fm, which have three-stage accumulations (Figure 22(a)). The first stage occurred from 33.9 Ma to 21 Ma with no more than 0.01 m3/m3 oil accumulated. The second stage occurred from 14 Ma to 11 Ma, with more than 0.038 m3/m3 oil accumulated. The third stage started at 5.0 Ma, with accumulated oil of 0.051 m3/m3 (Figure 22(g)). Petroleum currently accumulated in a series of strata, including Wenchang, Enping, Zhuhai, and Zhujiang Fm in the Huixi half-graben (Figure 22(a)). The modeling results indicate that the main period of hydrocarbon accumulation in the Zhujiang Fm is from 5.0 Ma to present day, with oil accumulated volume of 0.045 m3/m3 (Figure 22(c)). The Zhuhai Fm has two-stage petroleum accumulations. First stage is from 16.0 Ma to 11.0 Ma with more than 0.02 m3/m3 oil accumulated, and second stage started from 10.0 Ma with more than 0.02 m3/m3 oil accumulated (Figures 22(d) and 22(e)). Furthermore, hydrocarbons generated from the Eocene source rocks also migrated to the adjacent reservoirs and accumulated in the Wenchang and Enping Fm. Oil accumulated in the Enping Fm started from 15 Ma and temporarily ceased at 11 Ma, and then accumulation started again from 10 Ma to present day (Figure 22(f)). Petroleum accumulated in the Wenchang Fm began at 14 Ma, and accumulated oil volume reached peak of 0.048 0.02 m3/m3 at 10 Ma (Figure 22(h)).

From migration modeling results, the most likely fluid migration pathways are from the center of the depression to the flanks (Figure 22(a)) and structural high (Figure 21(a)), and the main driving forces for secondary petroleum migration are buoyancy and excess pressure caused by hydrocarbon generation of the deep source rocks (Figure 22(b)) and regional cap rocks (Figure 21(b)).

Faults are also important conduits for the Eocene source rocks to expel hydrocarbons towards traps in the structural high. Most major faults acted as migration pathways during the Eocene rift stage and then became sealed during the postrift stage. After that, previous reactivated faults and newly formed faults, resulting from the Dongsha Movement (9.8–4.4 Ma), served as migration pathways for petroleum. The hydrocarbon fluids accumulated in the deep Eocene reservoirs interval and then vertically migrated up to the overlying postrift strata via faults and then laterally migrated to the adjacent reservoirs via the homogeneous carrier beds (Figure 23(c)). The continuous homogeneous Zhuhai and Zhujiang sandstone, combined with structural highs, was considered to act as major pathways for long distance hydrocarbon migration in the HZS (Figure 23(b)).

4.3. Sensitivity Analysis

The crucial parameters were selected for parameter sensitivity analysis of the basin modeling. The sensitivity simulation suggests that paleogeothermal history affects hydrocarbon generation-expulsion histories, total organic carbon content of source rocks affects hydrocarbon expulsion, the permeability of fault zones strongly affects excess pressure in the deep sag, and the reservoir continuity significantly affects hydrocarbon accumulation.

4.3.1. Effect of Paleogeothermal History and TOC Content of Source Rocks on Hydrocarbon Generation-Expulsion Histories

Well 131 was selected to test the sensitivity of hydrocarbon generation-expulsion histories under different heat flow model. Although the amount of the generating and expelling hydrocarbon was the same, the process was slight different (Figure 24(a)). Moreover, the TOC content of source rocks also affects hydrocarbon generation-expulsion process, even though the generated hydrocarbons so far are the same (Figure 24(b)).

4.3.2. Excess Pressure

According to the measured pressure data (Figure 17), overpressure generally does not develop in the HZS. The seal capability of fault zones is considered to be the main controlling factor based on our interacted analysis. A boundary from profile AA′ was selected to analyze the overpressure condition under the different permeabilities of fault zones. Leak fraction is used to present permeability of all faulted cell in BasinMod 2D. No matter how the leak fraction was set, there was not much difference between the excess pressure when the boundary fault was set as permeable fault (Figures 25(a) and 25(b)). However, when the leak fraction was set to 0, which means that the fault was impermeable, the excess pressure in the Huinan half-graben was much higher than that with permeable boundary fault (Figure 25(c)). Hence, permeability of fault zones strongly affects the pressure field in the deep sag.

4.3.3. Reservoir Continuity

A thin continuous sand body was built at the bottom of Zhuhai Fm to test the sensitivity of continuous reservoir. When the lithology was set as pure sandstone, much more oil has migrated through continuous reservoir into the shallower Zhuhai and Zhujiang Fm (Figure 26(a)). However, when the lithology was set as a stratum with laterally variable lithology, most hydrocarbons migrating will be trapped in the deep Wenchang and Enping Fm close to the source rocks for the heterogeneous reservoirs (Figure 26(b)).

5. Conclusion

Basin modeling was conducted for four representative wells and two profiles in the Huizhou Sag (HZS) to reconstruct the process of generation, expulsion, migration, and accumulation hydrocarbon in Cenozoic Period since midcontinental rifting started in Eocene.

The Eocene Wenchang source rocks consist of shoreline-shallow lacustrine and deep lacustrine dark mudstone, and the kerogen is mainly of Type-II1, with minor component of Type-I. The Eocene Enping source rocks consist of swamp mudstone, swamp carbonaceous mudstone, and fluvial-delta mudstone and are generally composed of Type-II1 and Type-II2 kerogen.

The Eocene source rocks deposited in rift stage are currently in a middle-mature and late-mature stage, with % Ro value ranging from 0.7% to 3.0%. The Eocene source rocks developed in structural highs entered early mature stage at 23.0 Ma, while source rocks developed in deep central part of the sag have been in mature stage after 33.9 Ma. The Eocene source rocks started to generate hydrocarbons at 33.9 Ma. Hydrocarbon generation potential of the Wenchang source rock is much higher than that of the Enping source rock. Hydrocarbons were mainly expelled from the Wenchang Fm, and the contribution from the Enping Fm is relatively low.

Both short and long distance hydrocarbon migrations of the Eocene source rocks occurred in the HZS. Short distance hydrocarbon migration occurred in the Eocene reservoirs and is characteristics of multiperiod accumulation, while long distance hydrocarbon migration generally occurred in the Oligocene Zhuhai and Miocene Zhujiang reservoirs, which started to accumulate through 20~50 km homogeneous carrier beds after 10.0 Ma.

The hydrocarbons generated from the Eocene source rocks firstly migrated laterally to the adjacent Eocene reservoirs during the postrift stage, then vertically to the Oligo-Miocene carrier beds, and finally laterally through a long distance to the structural highs during the Neotectonic stage, which is controlled by both structural morphology and heterogeneity of carrier beds.

Fault is the most important conduit for hydrocarbon migration during the Neotectonic stage. Reactivation of previously existing faults and new-formed faults, caused by the Dongsha Movement (9.8–4.4 Ma), served as vertical migration pathways, which significantly influenced the timing of hydrocarbon accumulation in the postrift traps.

According to the sensitivity simulation of different factors, source rocks with high total organic carbon (TOC) content in the deeper sag could expel more hydrocarbons. Zones with impermeable faults are more likely to generate high excess pressure in the deeper sag. The continuous homogeneous reservoir contributes to hydrocarbon migration from the deep source rocks to shallow reservoir in the HZS.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by integrated research project, The Paleogene Sequence Stratigraphy and Conditions of Hydrocarbon Accumulation in Zhu 1 Depression, Pearl River Mouth Basin, sponsored financially by China National Offshore Oil Corporation (2009GYXQ02-06) and National Natural Science Foundation of China (Grant no. 41106108).