Quantitative Characterization of Excess Pressure Gradient in the Upper Interval of Es4 Member of Dongying Depression and Its Indicative Signiﬁcance for Oil Migration and Accumulation

: Excess pressure is the main driving force of oil migration in the source-reservoir system of overpressured petroliferous basins. It can reﬂect the change in driving force for oil migration and its inﬂuence on oil accumulation in overpressure transport layers. The drilling stem test (DST) data, well logging data, and seismic velocity data are used to describe the plane distribution of the excess pressures in the Es4s member of the Dongying Depression. Then, the values and directions of the excess pressure gradient, which can indicate oil migration and accumulation, are calculated based on the plane distribution of the excess pressure in the Es4s member of the Dongying Depression. The results suggest that overpressure is widely developed in the Es4s member of the Dongying Depression, and the excess pressure gradually decreases from the center to the edge of each sag, while the variation in the excess pressure gradient is characterized by “low-high-low” in a circular band around the sags. The excess pressure in the sag areas exceeds 15 MPa, but the excess pressure gradient is mainly between 0 and 1 MPa/km. The excess pressure in the northern steep slope zone of the Lijin sag and the northern steep slope zone of the Minfeng sag are less than 15 MPa, while the excess pressure gradient ranges from 1 to 7 MPa/km. The excess pressure in the central anticline belt and the gentle slope belt in the south of the Niuzhuang sag are between 0 and 15 MPa, and the excess pressure gradient is from 0 MPa/km to 2 MPa/km. From geochemical evidence, local oil migration directions indicated by the excess pressure gradient are consistent with those indicated by the ratio parameters of carbazole compounds in crude oil samples, indicating that the direction of the excess pressure gradient can indicate the dominant direction of oil migration driven by excess pressure, and the oil from the Es4s source rock is mainly distributed in the areas with a high excess pressure gradient or the areas with a low excess pressure gradient and low excess pressure (area II).


Introduction
Overpressure is widely developed in petroliferous basins [1][2][3]. The driving force of oil migration in overpressured petroliferous basins has always been important content of petroleum geology research [4], and it is also a difficult point of oil and gas accumulation research [5,6]. Bethke (1991) proposed that overpressure was the main driving force for oil and gas migration in the Illinois Basin [7]. Schegg (1999) believed that overpressure was more important than buoyancy in driving oil and gas migration [8]. Hubbert (1959) where ∅ is the fluid potential (mechanical energy per unit mass of fluid), kg*m 2 /s 2 ; P is the pore pressure of the measuring point, MPa; ρ is the density of the pore fluid, g/cm 3 ; g is the acceleration of gravity, 9.8 m/s 2 ; z is the elevation from the measuring point to the reference plane, m; q is the flow velocity of pore fluid, m/s. It can be seen that the fluid potential mainly consists of three parts: pressure energy, potential energy, and kinetic energy. Due to the slow flow velocity of oil and gas in the formation, the kinetic energy can be ignored, and the compressibility of the pore fluid in the formation is finite, so the density of the pore fluid can be regarded as a constant. That means the fluid potential can be calculated as: According to the second law of thermodynamics, oil and gas spontaneously flow from the high-potential area to the low-potential area. Based on Formula (2), oil and gas spontaneously flow from the areas with high excess pressure to the areas with low excess pressure if they are at the same depth. That means that the change in excess pressure controls oil and gas migration. Early research on the driving force of oil and gas migration mainly used the plane distribution map of the excess pressure to qualitatively describe the excess pressure variations. Feng et al. (2019) proposed the plane pressure decrease gradient in the Bonan sag of the Zhanhua Depression and explained the variation characteristics of the oil and gas migration model under the control of overpressure [10]. Liu et al. (2020) calculated the excess pressure gradient of the third member of the Shahejie Formation in the Bonan sag and analyzed the variation in the excess pressure gradient and its influence on the distribution of oil and gas [11]. However, few studies have noted the directional characteristics of the excess pressure gradient and its indicating significances on oil and gas migration. In addition, it is challenging to choose a reasonable step size when calculating the excess pressure gradient. In this paper, a new method is proposed to calculate the excess pressure gradient. We hope to obtain the direction and value of the excess pressure gradient, and, based on the results of the excess pressure gradient, we expect to explain the control on oil and gas migration and accumulation of the excess pressure gradient.
Previous studies have reported the significance of abnormal overpressure in the upper interval of the Es4 member of the Dongying Depression and the excess pressure range from 5.6 MPa to 30.9 MPa [1,[12][13][14]. However, there are no publications concerning the excess pressure gradient in Dongying Depression at present. In this paper, the new method is used to calculate the excess pressure gradient in upper interval of the Es4 member, Dongying Depression. Furthermore, we expect to explain the relationship between the excess pressure gradient and the oil fields that are from Es4s source rocks.

Geological Setting
The Dongying Depression is located in the East of China and is part of the Bohai Bay Basin. The Dongying Depression is NE-SW trending with an area of about 5700 km 2 . The depression is adjacent to Qingtuozi Uplift in the east, connected to Guanggrao Uplift and Luxi Uplift in the south, bounded by Qingcheng-Binxian Uplift in the west, and terminated at Chenjiazhuang Uplift in the north. Four groups of faults with different trending can be distinguished in the Dongying Depression: East-West (EW), North East-East to (NEE), North West-West (NWW), and North East (NE). The Dongying Depression consists of the Mingfeng Sag, Lijin Sag, Niuzhuang Sag, Boxing Sag, Northern Steep Slope zone, Southern Ramp zone, and Central Anticline Belt ( Figure 1) [15,16]. The basement of the Dongying Sag is Archaean and Paleozoic, and the sedimentary cover is Mesozoic and Cenozoic. The dark mudstones in the lower interval of the third member of the Paleogene Shahejie Formation (Es3x) and the upper interval of the fourth member of the Paleogene Shahejie Formation (Es4s) are the main source rocks in the study area with a total organic content (TOC) of 0.5-3.6% [1,17]. The thickness of the Es4s member in the Dongying Depression is about 344 m. The thickness of mudstones in the Es4s member ranges from 0 m to 172 m with an average of 58.92 m, while thickness of the sandstone reservoirs in the Es4s member ranges from 0 m to 322 m with an average of 273.24 m. In addition, the average ratios of mudstones and sandstone in the Es4s member account for 79.43% and 17.13%, respectively. The muddy source rocks interbedded with sandstone reservoirs are the most important petroleum system [14]. The sandstone lithology is mainly fine sandstone. The basement fractured metamorphic rocks, Mesozoic conglomerates, and Cenozoic sandstones are important reservoirs in the study area, and the mudstones of the First member of the Shahejie Formation of the Paleogene and the Minghuazhen Formation of the Neogene are stable regional cap rocks in the study area ( Figure 2). The exploration results show that there are many types of reservoirs in the Dongying Depression, including the structural reservoir, stratigraphic lithologic reservoir, and structure-stratigraphic lithologic compound reservoir [13,14]. A large amount of drill-stem test data (DST) show that a large-scale overpressure is developed in Es4 and Es3 of the Dongying Depression, and the maximum measured pore pressure coefficient (the pore pressure coefficient is the ratio of pore pressure and hydrostatic pressure) of the Es4s member can reach 1.99.
can be distinguished in the Dongying Depression: East-West (EW), North East-East to (NEE), North West-West (NWW), and North East (NE). The Dongying Depression consists of the Mingfeng Sag, Lijin Sag, Niuzhuang Sag, Boxing Sag, Northern Steep Slope zone, Southern Ramp zone, and Central Anticline Belt ( Figure 1) [15,16]. The basement of the Dongying Sag is Archaean and Paleozoic, and the sedimentary cover is Mesozoic and Cenozoic. The dark mudstones in the lower interval of the third member of the Paleogene Shahejie Formation (Es3x) and the upper interval of the fourth member of the Paleogene Shahejie Formation (Es4s) are the main source rocks in the study area with a total organic content (TOC) of 0.5-3.6% [1,17]. The thickness of the Es4s member in the Dongying Depression is about 344 m. The thickness of mudstones in the Es4s member ranges from 0 m to 172 m with an average of 58.92 m, while thickness of the sandstone reservoirs in the Es4s member ranges from 0 m to 322 m with an average of 273.24 m. In addition, the average ratios of mudstones and sandstone in the Es4s member account for 79.43% and 17.13%, respectively. The muddy source rocks interbedded with sandstone reservoirs are the most important petroleum system [14]. The sandstone lithology is mainly fine sandstone. The basement fractured metamorphic rocks, Mesozoic conglomerates, and Cenozoic sandstones are important reservoirs in the study area, and the mudstones of the First member of the Shahejie Formation of the Paleogene and the Minghuazhen Formation of the Neogene are stable regional cap rocks in the study area ( Figure 2). The exploration results show that there are many types of reservoirs in the Dongying depression, including the structural reservoir, stratigraphic lithologic reservoir, and structure-stratigraphic lithologic compound reservoir [13,14]. A large amount of drill-stem test data (DST) show that a large-scale overpressure is developed in Es4 and Es3 of the Dongying Depression, and the maximum measured pore pressure coefficient (the pore pressure coefficient is the ratio of pore pressure and hydrostatic pressure) of the Es4s member can reach 1.99.  sociated aqueous inclusions is 135.6 °C. On the basis of the burial-thermal maturation history of the well, it can be seen that the time of the period of oil charge was around 4 Ma ( Figure 3). The current vitrinite reflectance (Ro, %) of source rocks in Es4s ranges from 0.65% to 1.30%, which is still in the stage of massive oil generation [1]. The late period of overpressure in the Es4s member was mainly caused by oil generation [1,12,22]. Therefore, the current pressure distribution in the Es4s member of the Dongying Depression can approximately represent the pressure distribution during the late period of oil charge. The Dongying Depression experienced two periods of oil charging. The early oil charging occurred during the deposition period of the Dongying Formation   [18,19], dominated by low-mature oil, and the late oil charging occurred from the late deposition of the Guantao Formation to the present (10-0 Ma), which was the main accumulation period of the Dongying Sag [20,21]. Well L218 is located in the south of the Central Anticline zone in the Dongying Depression ( Figure 1). The sandstone sample from 3656.3 m in the Es4s member was investigated using an integrated fluid inclusion workflow. The results showed that the average homogenization temperature of oil inclusions in the sandstone sample is 99.35 • C, and the average homogenization temperature of associated aqueous inclusions is 135.6 • C. On the basis of the burial-thermal maturation history of the well, it can be seen that the time of the period of oil charge was around 4 Ma ( Figure 3). The current vitrinite reflectance (R o , %) of source rocks in Es4s ranges from 0.65% to 1.30%, which is still in the stage of massive oil generation [1]. The late period of overpressure in the Es4s member was mainly caused by oil generation [1,12,22]. Therefore, the current pressure distribution in the Es4s member of the Dongying Depression can approximately represent the pressure distribution during the late period of oil charge.

Samples and Carbazole Analysis
Oil samples from eight wells in the Es4s member were collected for carbazole analysis to define the main oil migration pathways ( Figure 1 and Table 1). N-hexane was used to precipitate the asphaltene and a chromatographic column was used to fractionate saturates, aromatics, and polar compounds. The chromatographic column consisted of silica gel and alumina with a ratio of 3:2 in volume. The silica gel must be activated by heating at 150 °C for 8 h in an oven, and the alumina would be activated by heating in a muffle furnace for 4 h with 400 °C. The saturated hydrocarbons were separated by sequential elution of n-hexane; the aromatic hydrocarbons were separated by the admixed dichloromethane and n-hexane admixture (ratio of 2:1 in volume); the polar compounds were separated by a chloroform and methanol solution admixture (ratio of 98:2 in volume). The amounts of n-hexane, the admixed dichloromethane and n-hexane admixture, and the chloroform and methanol solution admixture used in the measurement were 50 mL, 50 mL, and 30 mL, respectively. Carbazoles were separated from the resins on a chromatographic column (4 g of silicic acid) and eluted by the admixed n-hexane and toluene solution with a ratio of 1:1 in volume [23]. The 9-phenylcarbazole was added to the carbazoles to quantify the concentrations of carbazoles before performing GC-MS analysis, and the peak-area ratio method was used to calculate the concentrations of selected individual compounds.

Samples and Carbazole Analysis
Oil samples from eight wells in the Es4s member were collected for carbazole analysis to define the main oil migration pathways ( Figure 1 and Table 1). N-hexane was used to precipitate the asphaltene and a chromatographic column was used to fractionate saturates, aromatics, and polar compounds. The chromatographic column consisted of silica gel and alumina with a ratio of 3:2 in volume. The silica gel must be activated by heating at 150 • C for 8 h in an oven, and the alumina would be activated by heating in a muffle furnace for 4 h with 400 • C. The saturated hydrocarbons were separated by sequential elution of n-hexane; the aromatic hydrocarbons were separated by the admixed dichloromethane and n-hexane admixture (ratio of 2:1 in volume); the polar compounds were separated by a chloroform and methanol solution admixture (ratio of 98:2 in volume). The amounts of n-hexane, the admixed dichloromethane and n-hexane admixture, and the chloroform and methanol solution admixture used in the measurement were 50 mL, 50 mL, and 30 mL, respectively. Carbazoles were separated from the resins on a chromatographic column (4 g of silicic acid) and eluted by the admixed n-hexane and toluene solution with a ratio of 1:1 in volume [23]. The 9-phenylcarbazole was added to the carbazoles to quantify the concentrations of carbazoles before performing GC-MS analysis, and the peak-area ratio method was used to calculate the concentrations of selected individual compounds. The carbazoles were analyzed at the Key Laboratory of Tectonics and Petroleum Resources Ministry of Education, China by a gas chromatograph (Agilent 7890A), coupled with a mass spectrometer (GC-MS), equipped with a fused silica capillary column (DB-5MS 60 m × 0.25 mm × 0.25 µm). The carrier gas was helium (purity > 99.995%) and the inject rate of carrier gas was 1.0 mL/min. The experimental temperature during the test was kept at 50 • C for 1 min, and then increased up to 100 • C at a rate of 20 • C/min, up to 315 • C at a rate of 3 • C/min, and eventually remained at 315 • C for 31.8 min. The ion monitoring (SIM) mode was selected for the mass spectrometer at an electron impact of 70 eV during the test.

Homogenization Temperature of Oil and Aqueous Inclusions
The sandstone sample from Well L218 at 3656.3 m in the Es4s member was made into doubly polished sections at a thickness about 100 µm for fluid inclusion analysis. The homogenization temperature (Th) of oil and aqueous inclusions was carried out using a calibrated Linkam TH-600 stage equipped with a Zeiss Axiovert 200 microscope in the Key Laboratory of Tectonics and Petroleum Resources Ministry of Education, China University of Geosciences (Wuhan). The Th was measured at a heating rate of 10 • C/min to 60 • C and then at a heating rate of 5 • C/min. The measured temperature precision for Th was ±1 • C.

Drill-Stem Test Data
In this paper, a total of 592 drill-stem test data from 426 wells of the Es4s member in the Dongying Depression were collected from the Research Institute of Exploration and Development, Sinopec Shengli Oilfield Company. In order to avoid the invalid data influenced by seal failures, the drill-stem test data were carefully checked by combining mud weight log.

Excess Pressure Prediction Using Well Log Data
Excess pressure is the difference between pore pressure and hydrostatic pressure: where ∆P is the excess pressure, MPa; P p is the pore pressure, MPa; P h is the hydrostatic pressure, MPa. Previous studies have suggested that the overpressure generated by different mechanisms should be predicted in different methods [24]. The overpressure in the Dongying Depression is mainly caused by oil generation from the source rocks [1,25], while compaction disequilibrium may have contributed to the overpressure in the north Steep Slope zone of the Dongying Depression [26,27]. Thus, the pore pressure can be calculated by the Eaton Formula [28,29]: where P r is the vertical stress loaded by overlying strata, MPa; ∆t n is the acoustic travel time under the normal compaction trend, µs/m; ∆t i is the observed acoustic travel time from the well log, µs/m; N is an exponent, N = 2.3 in this paper; ρ w is the pore fluid density, g/cm 3 ; h is the depth of investigation, m. Therefore, pore pressure can be calculated by Formulas (4) and (5). The pore pressure calculated by Formula (4) must be corrected by the measured pore pressure before putting it into Formula (3). In this paper, a total of 110 wells were collected to calculate the pore pressure. A total of 158 measured pore pressures in the Es4s member were collected to correct the predicted pore pressure in adjacent mudstones. The plot of measured pore pressure versus predicted pore pressures indicates that the predicted pore pressures calculated by Formula (4) matched well with the measured pore pressure ( Figure 4). It reliably predicted pore pressures calculated by the Eaton Formula with N = 2.3. Furthermore, the fitting curve can be used to correct the predicted pore pressures. The correction formula can be inferred as: where P p is the predicted pore pressure after correction, MPa. Then, the predicted pore pressure after correction was put into Formula (3), and we obtained the excess pressure.
where is the vertical stress loaded by overlying strata, MPa; ∆ is the acoust time under the normal compaction trend, μs/m; ∆ is the observed acoustic tra from the well log, μs/m; is an exponent, = 2.3 in this paper; is the po density, g/cm 3 ; ℎ is the depth of investigation, m.
Therefore, pore pressure can be calculated by Formulas (4) and (5). The pore calculated by Formula (4) must be corrected by the measured pore pressure before it into Formula (3). In this paper, a total of 110 wells were collected to calculate pressure. A total of 158 measured pore pressures in the Es4s member were coll correct the predicted pore pressure in adjacent mudstones. The plot of measur pressure versus predicted pore pressures indicates that the predicted pore press culated by Formula (4) matched well with the measured pore pressure ( Figure 4) bly predicted pore pressures calculated by the Eaton Formula with N = 2.3. Furth the fitting curve can be used to correct the predicted pore pressures. The correc mula can be inferred as: where is the predicted pore pressure after correction, MPa. Then, the predic pressure after correction was put into Formula (3), and we obtained the excess pr

Excess Pressure Gradient
The excess pressure gradient is defined as the variation in the value of exce sure per unit distance in a specified direction, and in this direction, the excess changes the fastest. The excess pressure gradient is a vector with both magnit direction, and it can be calculated as:

Excess Pressure Gradient
The excess pressure gradient is defined as the variation in the value of excess pressure per unit distance in a specified direction, and in this direction, the excess pressure changes the fastest. The excess pressure gradient is a vector with both magnitude and direction, and it can be calculated as: where G p is the excess pressure gradient, MPa/km; dP x , dP y , and dP z are the variations in excess pressure in the directions of X, Y, and Z, respectively, MPa/km; d x , d y , and d z are the variations in distance in the direction of X, Y, and Z, respectively, km. The excess pressure gradient calculation in three-dimensional space is difficult to achieve at present because of the lack of continuous three-dimensional excess pressure data. In this paper, the plane distribution of excess pressure in the Es4s member with the contour map was used to calculate the excess pressure gradient in the two-dimensional plane. The pressure system consists of continuously distributed particles, and each particle has a unique excess pressure value ( Figure 5). As shown in Figure 5, the particle O has eight potential directions for the gradient calculation. Each direction can obtain a value of the excess pressure gradient. Then, the direction with the maximum value is the direction of the excess pressure gradient and the maximum value is the magnitude of the excess pressure gradient, such as the excess pressure gradient between O and e in Figure 5. The excess pressure gradient can be shown in the form of arrows. In this paper, the calculation of the excess pressure gradient in arrows was carried out on the Grid Vector Map in Surfer. The direction and the length of the arrow represent the direction and the magnitude of the excess pressure gradient, respectively. The length of the arrows was digitized and converted into the magnitude of the excess pressure gradient. The magnitude of the excess pressure gradient of each particle can be shown in the form of a contour map.
Energies 2022, 15, x FOR PEER REVIEW data. In this paper, the plane distribution of excess pressure in the Es4s member contour map was used to calculate the excess pressure gradient in the two-dim plane. The pressure system consists of continuously distributed particles, and ea cle has a unique excess pressure value ( Figure 5). As shown in Figure 5, the parti eight potential directions for the gradient calculation. Each direction can obtain a the excess pressure gradient. Then, the direction with the maximum value is the of the excess pressure gradient and the maximum value is the magnitude of th pressure gradient, such as the excess pressure gradient between O and e in Figu excess pressure gradient can be shown in the form of arrows. In this paper, the ca of the excess pressure gradient in arrows was carried out on the Grid Vector Map The direction and the length of the arrow represent the direction and the magn the excess pressure gradient, respectively. The length of the arrows was digit converted into the magnitude of the excess pressure gradient. The magnitude of t pressure gradient of each particle can be shown in the form of a contour map.

Drill Stem Tests
More than 590 measured pore pressure points of the Es4s member sandst 426 exploratory wells were collected to investigate the distribution of overpressu Dongying Depression. The profiles of pore pressures, excess pressure, and press ficients versus depths in the Es4s member are presented in Figure 4. Accordin profiles of pore pressures, excess pressure, and pressure coefficients versus de overpressure occurred from 2100 m. The excess pressure was over 10 MPa when t reached 2600 m and the excess pressure exceeded 20 MPa when the depth reac m ( Figure 6). The maximum excess pressure at the depth of 3574.5 m (Well L98) w MPa, while the pore pressure was 66.67 MPa.

Drill Stem Tests
More than 590 measured pore pressure points of the Es4s member sandstone from 426 exploratory wells were collected to investigate the distribution of overpressure in the Dongying Depression. The profiles of pore pressures, excess pressure, and pressure coefficients versus depths in the Es4s member are presented in Figure 4. According to the profiles of pore pressures, excess pressure, and pressure coefficients versus depths, the overpressure occurred from 2100 m. The excess pressure was over 10 MPa when the depth reached 2600 m and the excess pressure exceeded 20 MPa when the depth reached 3100 m ( Figure 6). The maximum excess pressure at the depth of 3574.5 m (Well L98) was 30.93 MPa, while the pore pressure was 66.67 MPa.

Well-Log Responses to Overpressure
The drill stem test (DST) can directly measure pore pressure in sandstones, while the mudstones cannot be directly measured, because of their low permeability. Therefore, the pore pressure in mudstones is commonly estimated based on well-log data. Well-log data from more than 300 wells were investigated. The acoustic travel time, resistivity log, and density log were used to analyze the differences between the overpressured and normally pressured regimes. Two representative wells, which were located on different structural zones in the Dongying Depression (Figure 1

Well-Log Responses to Overpressure
The drill stem test (DST) can directly measure pore pressure in sandstones, while the mudstones cannot be directly measured, because of their low permeability. Therefore, the pore pressure in mudstones is commonly estimated based on well-log data. Well-log data from more than 300 wells were investigated. The acoustic travel time, resistivity log, and density log were used to analyze the differences between the overpressured and normally pressured regimes. Two representative wells, which were located on different structural zones in the Dongying Depression (Figure 1 According to the well-log data, we suggest that the acoustic travel time is a reliable well-log response to overpressure in the Dongying Depression. The mudstones under the top of overpressure had a higher acoustic travel time than the normally pressured mudstones at a given depth, while the resistivity log and density log had no obvious indication of overpressure in the Dongying Depression (Figure 7). The varying trends of the acoustic travel time versus depths was different between normally pressured mudstones and overpressured mudstones. The acoustic travel time of normally pressured mudstones decreased with depth, while the acoustic travel time of overpressured mudstones increased with depth at the upper interval and then decreased with depth at the lower interval. According to the well-log data, we suggest that the acoustic travel time is a reliable well-log response to overpressure in the Dongying Depression. The mudstones under the top of overpressure had a higher acoustic travel time than the normally pressured mudstones at a given depth, while the resistivity log and density log had no obvious indication of overpressure in the Dongying depression ( Figure 7). The varying trends of the acoustic travel time versus depths was different between normally pressured mudstones and overpressured mudstones. The acoustic travel time of normally pressured mudstones decreased with depth, while the acoustic travel time of overpressured mudstones increased with depth at the upper interval and then decreased with depth at the lower interval.

Distribution of Excess Pressure
The excess pressure calculated by drill stem tests are direct evidence for identifying overpressure. However, there are no drill stem tests in mudstones, because of the low permeability. The relationship between acoustic travel time and pore pressure was used to obtain the pore pressure in mudstones. The vertical distribution of excess pressure showed an increasing trend from the shallow burial regions to the deep burial regions (Figures 6 and 7). According to the plot of drill stem tests versus depths (Figure 6

Distribution of Excess Pressure
The excess pressure calculated by drill stem tests are direct evidence for identifying overpressure. However, there are no drill stem tests in mudstones, because of the low permeability. The relationship between acoustic travel time and pore pressure was used to obtain the pore pressure in mudstones. The vertical distribution of excess pressure showed an increasing trend from the shallow burial regions to the deep burial regions (Figures 6 and 7). According to the plot of drill stem tests versus depths (Figure 6), the excess pressure was less than 3 MPa at depths from 0 m to 2000 m, most of the excess pressure ranged from 5 MPa to 10 MPa at depths from 2000 m to 2600 m, and the excess pressure exceeded 20 MPa at the depth of 3000 m. The excess pressure in mudstones of the Es4s member calculated by the acoustic travel time using the Eaton formula continuously varied with depth ( Figure 7). The excess pressure of mudstones in the Es4s member exceeded 20 MPa at Well T716 and was approximately 15 MPa at Well L105.
The lateral distribution of excess pressure in mudstones of the Es4s member are shown in Figure 8. It can be seen that the overpressure area of the Es4s member was mainly distributed in the Lijin Sag, Minfeng Sag, and Niuzhuang Sag. The excess pressure in mudstones of the Es4s member in the Lijin Sag ranged from 15 MPa to 35 MPa, the excess pressure in mudstones of the Es4s member in the Minfeng Sag and Niuzhuang Sag ranged from 15 MPa to 27 MPa, while the excess pressure in mudstones of the Es4s member in the Boxing Sag and Central Anticline Belt were mainly less than 15 MPa. The area between the excess pressure contour of 0 MPa and the boundary of the Dongying Depression was the normally pressured area without color filling in the depression (Figure 8). mainly distributed in the Lijin Sag, Minfeng Sag, and Niuzhuang Sag. The excess pressure in mudstones of the Es4s member in the Lijin Sag ranged from 15 MPa to 35 MPa, the excess pressure in mudstones of the Es4s member in the Minfeng Sag and Niuzhuang Sag ranged from 15 MPa to 27 MPa, while the excess pressure in mudstones of the Es4s member in the Boxing Sag and Central Anticline Belt were mainly less than 15 MPa. The area between the excess pressure contour of 0 MPa and the boundary of the Dongying Depression was the normally pressured area without color filling in the depression (Figure 8).

Distribution of Excess Pressure Gradient
The lateral distribution of the excess pressure gradient in mudstones of the Es4s member can be shown in the form of arrows or in the form of a contour map (Figure 9). The direction of the arrows is expressed as the excess pressure gradient direction, while the length of the arrows represents the value of the excess pressure gradient (Figure 9a). The results showed that the direction of the pore pressure gradient of the Es4s member in the Dongying Depression diverged from the overpressured area of each sag to the periphery of the Dongying Depression. The excess pressure gradient direction of each point in the depression took the form of a series of continuous lines. The excess pressure gradient in the Es4s member at the Lijin Sag had two main directions and either pointed to the Steep Slope zone converging on Chengjiazhuang Uplift and Bingxian Uplift, or pointed to the Central Anticline Belt. The excess pressure gradient in the Es4s member at the Niuzhuang Sag could be distinguished into two main directions, but most of the excess pressure gradients pointed to the Southern Ramp zone. The excess pressure gradient in the

Distribution of Excess Pressure Gradient
The lateral distribution of the excess pressure gradient in mudstones of the Es4s member can be shown in the form of arrows or in the form of a contour map (Figure 9). The direction of the arrows is expressed as the excess pressure gradient direction, while the length of the arrows represents the value of the excess pressure gradient (Figure 9a). The results showed that the direction of the pore pressure gradient of the Es4s member in the Dongying Depression diverged from the overpressured area of each sag to the periphery of the Dongying Depression. The excess pressure gradient direction of each point in the depression took the form of a series of continuous lines. The excess pressure gradient in the Es4s member at the Lijin Sag had two main directions and either pointed to the Steep Slope zone converging on Chengjiazhuang Uplift and Bingxian Uplift, or pointed to the Central Anticline Belt. The excess pressure gradient in the Es4s member at the Niuzhuang Sag could be distinguished into two main directions, but most of the excess pressure gradients pointed to the Southern Ramp zone. The excess pressure gradient in the Es4s member at the Mingfeng Sag mainly pointed to Chenjiazhuang Uplift and Qingtuozi Uplift. In addition, the length of the arrows at the northern Steep Slope zone was longer than the arrows at other areas in the Dongying Depression. That means that the values of the excess pressure gradient in the northern Steep Slope zone were higher than those in other areas in the Dongying Depression.

The Excess Pressure Gradient Indicates the Direction of Oil Migration
Carbazoles, methylcarbazoles, dimethylcarbazoles, and benzocarbazoles were detected in the eight oil samples from the Dongying Depression (Table 1). The concentration of carbazoles, methylcarbazoles, dimethylcarbazoles, and benzocarbazoles varied from 0.87 μg/g to 5.01 μg/g, 3.30 μg/g to 13.44 μg/g, 4.79 to 23.66 μg/g, and 1.23 to 4.93 μg/g, respectively (Table 2). For dimethylcarbazoles, the ratios of 1,8-/2,7-dimethylcarbazoles The values of the excess pressure gradient presented "low-high-low" from the sag area to the periphery of the Dongying Depression (Figure 9b). The low excess pressure gradient could be identified as two types. Type I was the area with a low excess pressure gradient and high excess pressure: the excess pressure gradient was less than 1 MPa/km, while the excess pressure was more than 15 MPa. This type was mainly distributed in the sag areas (area I in Figure 9b). Type II was the area with a low excess pressure gradient and low excess pressure: the excess pressure gradient was less than 1 MPa/km and the excess pressure was less than 15 MPa. This type was mainly distributed in the periphery of the Dongying Depression (area II in Figure 9b). The high excess pressure gradients were mainly distributed in the faults around the sags and ranged from 1 MPa/km to 7 MPa/km. Carbazoles, methylcarbazoles, dimethylcarbazoles, and benzocarbazoles were detected in the eight oil samples from the Dongying Depression (Table 1). The concentration of carbazoles, methylcarbazoles, dimethylcarbazoles, and benzocarbazoles varied from 0.87 µg/g to 5.01 µg/g, 3.30 µg/g to 13.44 µg/g, 4.79 to 23.66 µg/g, and 1.23 to 4.93 µg/g, respectively ( Table 2). For dimethylcarbazoles, the ratios of 1,8-/2,7-dimethylcarbazoles (1,8-/2,7-DMC) ranged from 2.17 to 7.01 (  Pyrrolic nitrogen compounds within crude oil are widely applied in the research of oil migration, and the change in the carbazole parameters can be used to trace the direction of oil migration [30][31][32]. During the secondary migration of oil, carbazoles with strong polarity are more easily absorbed by reservoir rocks than carbazoles with less polarity. As the secondary migration distance of oil increases, the relative content of nitrogen-shielded isomers increases, while the relative content of nitrogen-exposed isomers decreases. In this study, carbazole parameters such as 1,  Figure 10).

Driven Force of Oil Migration
Oil migration in overpressured petroliferous basins is mainly affected by buoyancy, gravity, excess pressure, and capillary force. The resultant force of buoyancy and gravity is net buoyancy, net buoyancy and excess pressure are the driving force of oil migration, while the capillary force is the resistance of oil migration. The driving force must be great enough to overcome capillary resistance during oil migration. If net buoyancy is the only driving force for oil migration, enough oil column is needed to produce a sufficient net buoyancy to overcome the capillary force. The mercury injection test showed that the maximum connecting pore throat radius of fine sandstone in the Es4s member of the Dongying Depression ranged from 0.551 μm to 7.129 μm. It was calculated that the capillary resistance during oil migration in the presence of oil-water phases in the fine sandstone reservoir was between 1.333 MPa and 6.431 MPa. Taking it as the boundary condition, the critical oil column height in the vertical direction for oil migration in the Es4s formation must be 755 m to 3645 m at least. The dip angle of the Es4s formation was about 5°; therefore, the critical oil column contiguous length for oil migration in the Es4s member must The oil migration direction traced by carbazoles concurs with that indicated by the excess pressure gradient, and the area where the excess pressure gradient arrows converge is often the oil accumulation area, which suggests that the traced oil migration pathways are reliable in the Dongying Depression (Figure 9a). In addition, there are few oil accumulations in the area where the direction of the excess pressure gradient arrows is parallel or divergent, which implies that the direction of the excess pressure gradient can be used to trace oil migration and accumulation.

Driven Force of Oil Migration
Oil migration in overpressured petroliferous basins is mainly affected by buoyancy, gravity, excess pressure, and capillary force. The resultant force of buoyancy and gravity is net buoyancy, net buoyancy and excess pressure are the driving force of oil migration, while the capillary force is the resistance of oil migration. The driving force must be great enough to overcome capillary resistance during oil migration. If net buoyancy is the only driving force for oil migration, enough oil column is needed to produce a sufficient net buoyancy to overcome the capillary force. The mercury injection test showed that the maximum connecting pore throat radius of fine sandstone in the Es4s member of the Dongying Depression ranged from 0.551 µm to 7.129 µm. It was calculated that the capillary resistance during oil migration in the presence of oil-water phases in the fine sandstone reservoir was between 1.333 MPa and 6.431 MPa. Taking it as the boundary condition, the critical oil column height in the vertical direction for oil migration in the Es4s formation must be 755 m to 3645 m at least. The dip angle of the Es4s formation was about 5 • ; therefore, the critical oil column contiguous length for oil migration in the Es4s member must be 8.6 km to 41.7 km, while the distance from the overpressured center of each sag in the Dongying Depression to the sag boundary was about 5-20 km. It can be determined that it is difficult to generate enough net buoyancy to overcome capillary resistance and drive oil migration in the overpressure system of the Es4s member. However, the excess pressure in the overpressure system of the Es4s member exceeded 15 MPa (Figures 6 and 8). It was significantly greater than the capillary resistance of oil migration with the presence of two oil-water phases in the fine sandstone reservoir of the Es4s formation. Therefore, it can be inferred that the excess pressure is the main driving force of oil migration in the overpressured area of the Es4s formation.
The representative profile of the excess pressure gradient and oil reservoirs in the Dongying Depression are shown in Figure 11 (the location of the profile is shown in Figure 9a). It can be seen that the excess pressure in the Es4s member exceeded 20 MPa. The high excess pressure gradient (2 MPa/km-7 MPa/km) mainly distributed between Well T152 and Well T125 around the ShengBei Fault. The crude oil was generated in the sag, and oil laterally migrated under the driving by excess pressure. When the oil reaches the fault, the overpressure hydrocarbon-bearing fluid is released rapidly if the fault is under poor preservation conditions. Furthermore, the excess pressure gradient increases and the excess pressure decreases rapidly. When the excess pressure is insufficient for the oil to overcome the capillary resistance, the oil migration stops under the influence of capillary resistance for the moment. If there is a trap, the oil accumulates and forms a reservoir there, such as the oil reservoirs near the fault and Well T764 ( Figure 11).  (Figures 6 and 8).
It was significantly greater than the capillary resistance of oil migration with the presence of two oil-water phases in the fine sandstone reservoir of the Es4s formation. Therefore, it can be inferred that the excess pressure is the main driving force of oil migration in the overpressured area of the Es4s formation. The representative profile of the excess pressure gradient and oil reservoirs in the Dongying Depression are shown in Figure 11 (the location of the profile is shown in Figure  9a). It can be seen that the excess pressure in the Es4s member exceeded 20 MPa. The high excess pressure gradient (2 MPa/km-7 MPa/km) mainly distributed between Well T152 and Well T125 around the ShengBei Fault. The crude oil was generated in the sag, and oil laterally migrated under the driving by excess pressure. When the oil reaches the fault, the overpressure hydrocarbon-bearing fluid is released rapidly if the fault is under poor preservation conditions. Furthermore, the excess pressure gradient increases and the excess pressure decreases rapidly. When the excess pressure is insufficient for the oil to overcome the capillary resistance, the oil migration stops under the influence of capillary resistance for the moment. If there is a trap, the oil accumulates and forms a reservoir there, such as the oil reservoirs near the fault and Well T764 ( Figure 11).  Figure 9b).
The distribution of oil reservoirs (related to the Es4s source rock) and the excess pressure gradient indicates the close affinity between them ( Figure 12). It can be seen that the oil reservoirs that are sourced from the Es4s source rock are mainly accumulated in highexcess-pressure-gradient areas or a low excess pressure gradient with low-excess-pressure areas (area II), and only a few oil reservoirs are distributed in a low excess pressure gradient with high-excess-pressure areas (area I). The oil reservoirs in area I are mainly lithologic reservoirs in lentoid sandy [12,33]. For the lack of large-scale sandstone reservoirs in Es4s, the oil generated in the Es4s source rock laterally migrates to periphery of the sag. The high excess pressure gradient in the nearby sag means the excess pressure decreasing fast here. That indicates that there exists a good carrier bed for oil migration.  Figure 9b).
The distribution of oil reservoirs (related to the Es4s source rock) and the excess pressure gradient indicates the close affinity between them ( Figure 12). It can be seen that the oil reservoirs that are sourced from the Es4s source rock are mainly accumulated in highexcess-pressure-gradient areas or a low excess pressure gradient with low-excess-pressure areas (area II), and only a few oil reservoirs are distributed in a low excess pressure gradient with high-excess-pressure areas (area I). The oil reservoirs in area I are mainly lithologic reservoirs in lentoid sandy [12,33]. For the lack of large-scale sandstone reservoirs in Es4s, the oil generated in the Es4s source rock laterally migrates to periphery of the sag. The high excess pressure gradient in the nearby sag means the excess pressure decreasing fast here. That indicates that there exists a good carrier bed for oil migration. It can be seen that the faults match well with the high-excess-pressure-gradient areas. Thus, it is inferred that the faults in the high-excess-pressure-gradient areas are the main carrier bed for oil migration (Figures 11 and 12). The trap in the high-excess-pressure-gradient area will be the advantageous area for forming oil reservoirs. If there is no trap, the oil from the Es4s source rock will migrate through faults to the upper formation. When the oil migrates to area II, it will be difficult for the excess pressure to overcome the capillary resistance. The oil will accumulate until the oil column provides enough net buoyancy to overcome the capillary resistance and continue migrating. gradient area will be the advantageous area for forming oil reservoirs. If there is no trap, the oil from the Es4s source rock will migrate through faults to the upper formation. When the oil migrates to area II, it will be difficult for the excess pressure to overcome the capillary resistance. The oil will accumulate until the oil column provides enough net buoyancy to overcome the capillary resistance and continue migrating.

Conclusions
(1) The measured excess pressure in the Es4s member ranges from 5.65 MPa to 30.93 MPa. The results of overpressure prediction show that the overpressure in the Es4s member is mainly distributed in the Lijin Sag, Minfeng Sag, and Niuzhuang Sag, and the excess pressure ranges from 5 MPa to 30 MPa, while the excess pressure in the Boxing Sag is mainly less than 15 MPa. (2) The new method can calculate the values and the directions of the excess pressure gradient. The values of the excess pressure gradient in the Es4s member of the Dongying Depression present as "low-high-low" from the sag area to the periphery of the depression. The directions of the excess pressure gradient are consistent with the directions of oil migration indicated by the ratio of carbazole compounds, which suggests that the directions of the excess pressure gradient can be used to trace oil migration and accumulation.

Conclusions
(1) The measured excess pressure in the Es4s member ranges from 5.65 MPa to 30.93 MPa.
The results of overpressure prediction show that the overpressure in the Es4s member is mainly distributed in the Lijin Sag, Minfeng Sag, and Niuzhuang Sag, and the excess pressure ranges from 5 MPa to 30 MPa, while the excess pressure in the Boxing Sag is mainly less than 15 MPa. (2) The new method can calculate the values and the directions of the excess pressure gradient. The values of the excess pressure gradient in the Es4s member of the Dongying Depression present as "low-high-low" from the sag area to the periphery of the depression. The directions of the excess pressure gradient are consistent with the directions of oil migration indicated by the ratio of carbazole compounds, which suggests that the directions of the excess pressure gradient can be used to trace oil migration and accumulation. (3) The oil migration in overpressured areas is driven by excess pressure. With the excess pressure decreasing, the net buoyancy kicks gradually in the high-excess-pressuregradient areas. Before the net buoyancy overcomes the capillary resistance and drives the oil migration, the oil will continue accumulating at the reservoir until the oil column provides enough net buoyancy to overcome the capillary resistance. (4) The distribution of the excess pressure gradient in the Es4s formation in the study area has a good matching relationship with the distribution of the oil field. The oil from the Es4s source rock is mainly distributed in the areas with a high excess pressure gradient or areas with a low excess pressure gradient and low excess pressure (area II).