The Yinggehai Basin is a typical high temperature and high pressure (HTHP) gas-bearing basin. The pressure coefficient exceeds 2.2 in deeply-buried Miocene reservoirs in the Ledong Slope, a nondiapir zone in the Yinggehai Basin. Determining the overpressure mechanisms and predicting the pore pressure are key issues for natural gas exploration and development in the Ledong Slope. In this paper, overpressure mechanisms were investigated according to the analysis of vertical effective stress-logging responses and geological evaluations, and the pore pressure was predicted using the Bowers method. The loading-unloading crossplots indicated that the overpressure that existed in reservoirs mainly consists of two types: neighbor-source and allo-source overpressure. The neighbor-source overpressure is mainly caused by the pressure transmission from the adjacent mudstone to the reservoir, with a pressure coefficient less than 1.5 ~ 1.6. The high-magnitude overpressure points with pressure coefficients greater than 1.6 show a typical unloading response, indicating elevated sandstone pressures rather than in situ mudstone pressures, which are most likely to be generated by overpressure vertical transfer. The high-magnitude overpressure fluid generated by the high mature ultradeep buried N1s source rock migrated to the shallower reservoirs via hidden faults/microfractures, which led to the vertical transfer of overpressure. Vertically transferred overpressure was generated at 1.5 ~0.2 Ma, which is beneficial for the preservation of overpressure in lenticular sandbodies. The estimated pore pressure by the Bowers method is in good agreement with the measured pressure and provides a meaningful reference for predrilling pressure prediction in nondiapir or diapir zones in the Yinggehai Basin.

1. Introduction

The pore pressure is called abnormally high pressure or overpressure when it is greater than the hydrostatic pressure expected at a given depth [1, 2]. Great risk can be posed due to unexpected overpressure zones during drilling operations, which is harmful to both human life and environment [3, 4]. Significant progress has been made in research on the origin of overpressure and pore pressure prediction in low-permeability formations. The low-permeability mudstones play a critical role in the generation and preservation of abnormal pressure [1, 3, 5]. In situ overpressures in sedimentary basins are commonly attributed to disequilibrium compaction or fluid expansion (gas generation) [13, 5, 6]. Generally, the overpressure in mudstone cannot be sustained for long periods of geological time, with the overpressure dissipating via fluid leakage [79].

The overpressure redistribution phenomenon commonly occurs in permeable aquifers. It is equally important to research the redistribution of overpressure as it is to understand the overpressure mechanism. According to pore fluid sources, overpressure in permeable formations can be divided into three types: self-source, neighbor-source, and allo-source overpressure [10]. Self-source overpressure is caused by the buoyancy of the hydrocarbon column, oil-to-gas cracking, precipitation of dissolved gas [3, 5, 10]. Neighbor-source overpressure refers to the overpressure phenomenon in permeable formations caused by the transmission of pressure from adjacent high-pressure low-permeability formations, and the pressure in the mudstone is equal to that in the sandstone. Commonly used pressure prediction methods including the equivalent depth method and Eaton method can be successfully applied to the prediction of neighbor-source overpressure [11, 12]. Allo-source overpressure means that the pressure results from the overpressure transfer due to the hydrodynamic connection between different overpressure systems, including lateral transfer in inclined laterally continuous sandstones (centroid effect) and vertical transfer via opening faults [7, 13, 14]. In particular, faults can extend to ultradeep high-magnitude overpressure systems, which lead to a very high overpressure in shallow permeable formations [15].

The vertical effective stress-porosity/sonic velocity plots (also known as loading-unloading curves) are used to distinguish between neighbor-source overpressure and allo-source overpressure [3, 7, 16, 17]. Pore pressures relate to disequilibrium compaction commonly parallel to vertical stress; consequently, the vertical effective stress remains constant with burial, thus, the relationship will remain on the loading curve. While the increasing pore pressure causes the vertical effective stress to decrease when formations become overpressured either by fluid expansion or vertical transfer. However, the compaction is mostly an irreversible process, overpressures generated by fluid expansion or vertical transfer that follow a path away from the loading curve [18, 19]. Unfortunately, the loading-unloading analysis cannot distinguish between the fluid expansion and vertical transfer mechanisms [3, 7].

Accurate prediction of allo-source overpressure is another difficult issue during drilling engineering [20, 21]. Vertically transferred overpressure is always generated after normal compaction of formations, which is only associated with a small change in porosity because of the slight elastic contraction of sediment grains [19]. The sonic log response is smaller in sediments where overpressure is generated by vertical transfer. Tingay et al. [2] attempted to predict vertically transferred overpressure using the Eaton method and proposed that a higher Eaton exponent is required to amplify the weak sonic log response. However, the larger Eaton exponent also amplifies any small sonic variation not associated with overpressure, resulting in a noisy pore pressure prediction in hydrostatic and slightly overpressure sediment. The Eaton method seems cannot satisfactorily deal with the problems of vertically transferred overpressure prediction.

As a whole, although allo-source overpressure caused by pressure transfer is not a new observation, the lateral or vertical transfer of overpressure within sedimentary basins has not been widely reported, with a few notable exceptions [7, 15, 17, 2022]. At present, studies on allo-source overpressure are mostly limited to the qualitative description of measured overpressure characteristics and the possible generation process. The identification and accurate prediction of transfer overpressure using measured pressure and logging data are still weaknesses in overpressure research.

The Yinggehai Basin is a typical high temperature and high pressure (HTHP) basin with rich natural gas resources and is as famous as the Gulf of Mexico and North Sea Basin. The overpressure origins in the Gulf of Mexico and North Sea Basin have been thoroughly studied by researchers and recognized that the disequilibrium compaction, gas generation, and chemical compaction are dominated overpressure mechanisms [2326]. However, the temperature and pore pressure in Yinggehai Basin is higher than that in the Gulf of Mexico and North Sea Basin, and the overpressure mechanism is more complicated. The highest temperature in the Yinggehai Basin is approximately 250°C, and the pressure coefficient reaches 2.38. The depth of overpressure is shallow, overpressure increases quickly with depth, and the safety drilling window is very narrow [27]. Some studies have incorporated detailed petrophysical, geochemical, geophysical, and basin modeling into the analysis of the overpressure mechanism and prediction in the Yinggehai Basin [2831]. It has been recognized that the disequilibrium compaction of mudstone caused by high deposition rates since the Miocene is the most common mechanism of overpressure in the Yinggehai Basin [15, 28, 30]. The overpressure contribution of hydrocarbon generation is negligible in the Miocene formations that buried shallower than 4500 m [15]. The development of high-magnitude overpressure in the shallow permeable sandstone is closely related to overpressure transfer by opening faults [15, 28]. In particular, the mud diapir structures in the Yinggehai Basin have greatly changed the distribution pattern of pore pressure [29, 32]. The overpressure is transferred from the deep overpressure system in the diapir belt to the shallow system through faults and fractures, which leads to the development of overpressure in shallow reservoirs [15, 28, 32]. Self-source overpressure caused by disequilibrium compaction is believed to be the most important overpressure mechanism in the nondiapir zone [31]. Previous researchers have tried to improve the existing prediction methods, but due to the complex mechanism of overpressure and the puzzling logging response of overpressure, overpressure prediction in the Yinggehai Basin still has difficulty meeting the needs of safe drilling, which easily causes overflow leakage or even blowout [33, 34].

In recent years, the exploration focus of the Yinggehai Basin has begun to shift to lithologic reservoirs in the Ledong Slope, which shows a large burial depth (more than 4000 m), high temperature (a temperature greater than 200°C, and a temperature gradient reaching 46.8°C/km), and high overpressure (a pressure coefficient close to 2.30) [35]. There have been few detailed studies on the mechanism of overpressure generation in the LD-B Wellblock until now. Considering that the Ledong Slope is located in nondiapir zone, researchers attribute the mechanism of overpressure in LD-B Wellblock to disequilibrium compaction and hydrocarbon generation [31, 36]. However, the equivalent depth method and Eaton method hardly accurately determine the pore pressure [34, 36]. In addition, it is difficult to reasonably explain why the pore pressure among the same sandbodies shows obvious differences in different wells. In fact, the pore pressure of the deep buried reservoirs in the LD-B Wellblock is controlled by sedimentary and burial conditions, sandbody connectivity, and fault activity [35]. There are multiple overpressure mechanisms, and it is difficult to predict and monitor the pore pressure. Therefore, it is worth further understanding the overpressure mechanisms of the LD-B Wellblock.

This paper identifies the intrinsic mudstone pressure and determines whether the reservoir pressures are appropriate or anomalous using loading-unloading curves combined with 1-D basin modeling. The possibility of vertically transferred overpressure is confirmed with geological condition analysis. On this basis, the pressure prediction model parameters are optimized to realize the accurate prediction of the allo-source overpressure. The results of this study are expected to reduce possible accident such as kicks and blowouts during drilling to ensure drilling safety and to provide effective guidance for developing oil and gas resources.

2. Geological Setting

The Yinggehai Basin is a Cenozoic sedimentary basin located in the South China Sea and is diamond-shaped along the northwest-southeast direction, which is controlled by the Honghe strike-slip fault. The northeastern side of the basin is adjacent to the Beibuwan Basin and Hainan Island, and it is connected to the Kuntum Uplift to the west and the Qiongdongnan Basin to the southeast. The Yinggehai Basin is further divided into the Hanoi Depression, Lingao Uplift, Eastern Slope, Central Depression, and Western Slope. Due to dextral strike-slip, rows of left-order enechelon diapir structures with nearly north-south strikes are developed in the Central Depression, namely, the Central Mud Diapir Belt (Figure 1(a)). The LD-B Wellblock is in the southeastern part of the central depression of the Yinggehai Basin (Figure 1(b)), and the average seawater depth is approximately 85 m. The tectonic evolution in the Yinggehai Basin may be divided into two stages: the Paleogene rifting stage and the Neogene-Quaternary postrifting thermal subsiding stage.

The thickness of the Cenozoic sedimentary strata in the Yinggehai Basin is more than 17 km and consists of the Oligocene Yacheng Formation (E3y) and Lingshui Formation (E3l); the Miocene Sanya Formation(N1s); Meishan Formation (N1m) and Huangliu Formation (N1h); the Pliocene Yinggehai Formation (N2y); and the Quaternary Ledong Formation (Qld) (Figure 2). The Yacheng Formation and Lingshui Formation comprise carbonaceous mudstones and coal seams formed in semienclosed shallow sea and delta front.

The Sanya Formation can be divided into two members: the first member (N1s1) is mainly interbedded with sandstone and mudstone, and the second member (N1s2) composed of sandstones intercalated with thin mudstone. The Meishan Formation is composed of argillaceous siltstone and gray mudstone deposited in shallow to semi-deep sea environments, and it also contains submarine fan sandstone. The Huangliu Formation is further divided into two members. The first member (N1h1) interbedded argillaceous siltstone and silty mudstone, and the second member (N1h2) contains gray mudstone and fine sandstone. The Yinggehai Formation is dominated by siltstone and gray mudstone of semideep sea and delta deposition. The Ledong Formation contains shallow marine deposits, and the lithology is mainly gray clay and sandy conglomerates.

The marine mudstone of the N1s and N1m formations are rich in type III kerogen, which is the dominant source rock in the Yinggehai Basin [37]. Multiple sets of reservoirs were developed in the overlying N1h, N2y, and Qld formations. The N2y2 to Qld formations belong to shallow assemblages, the burial depth is relatively shallow, and the pressure coefficient ranges from 1.0 to 1.5. The N1m to N1h formations are deep buried assemblages, the temperature gradient is more than 45°C/km, and the pressure coefficient exceeds 1.8.

The gravity flow sediments and submarine fan sandstones in the N1h and N1m formations are the main gas reservoirs and are overlain by the thick mudstone of the N1y2 in the LD-B Wellblock. The dextral strike-slip tensional faults and fractures provide the migration pathway for natural gas in the area, forming a large HTHP lithological gas reservoir [35].

3. Data and Method

3.1. Data Source

The accuracy of the modular formation dynamics tester (MDT) pressure is affected by many factors, such as the lithology and permeability of the formation, borehole conditions, and mud properties. In most cases, trustless data are shown in the form of seat sealing failure, superpressure, and dry tests. The mobility (the ratio of permeability and viscosity) was chosen as a key factor to determine the validity of the MDT results, and the relative values of mobility corresponding to the effective points, low-permeability points, and dry points were 1 ~ 20md/cp, 0.01~0.5md/cp, and lower than 0.01md/cp, respectively. Finally, total of 269 points, including 63 effective points, 82 low-permeability points, 59 dry points, 60 sealing failure points, and 5 superpressure points, were classified; only the effective points can be used for pore pressure analysis.

The final shut-in pressure recorded by the drilling stem test (DST) is close to the original formation pressure when the shut-in time is sufficiently long, and the permeability of the formation is high enough. In addition to the depressurization caused by gas development, the shut-in pressure is often affected by factors such as inherently low permeability of formations, incomplete pressure recovery, and local mud-cake problems, so the shut-in pressure cannot be used to reflect the pore pressure directly. Therefore, the measured pressure data need to be screened to determine whether DST data are acceptable [39]. DST has rarely been performed in the LD-B Wellblock, and only 5 effective DST points were selected.

Although the mud weight is not an accurate representation of pore pressure, when balanced or underbalanced drilling technology is used, the mud weight can be a proxy for pore pressure calibrated by the measured pressure [39]. The mud weights from more than 20 wells in LD-B Wellblock are selected and converted into equivalent pressures.

By comparing the pore pressure from different sources at similar depths, the MDT pressures are in good agreement with the pressures from DST or mud weight equivalent pressures at corresponding depths (Figure 3). This shows that the measured data can reliably reflect the characteristics of underground pressure after screening.

The carbon isotope of the methane (C1), ethane (C2) in natural gas, and the carbon dioxide (CO2) was measured using SP3400 gas chromatography-combustion-isotope ratio MAT252 mass spectrometer (GC-CIRMS). The carbon isotopic data of the source rocks were derived from Huang et al. [37]. The discrimination plot of different organic gases using δ13C1-δ13C2-δ13C3 of natural gas samples was modified from [40].

As the most widely used method for studying the stages of hydrocarbon accumulation, fluid inclusions also record good geohistorical information, which can effectively determine the fluid properties at the time of accumulation. The homogenization temperature is a minimum trapping temperature of fluid inclusions and represents the minimum precipitation temperature of the host mineral. Homogenization temperature measurement of fluid inclusions was carried by Shui et al. [41], using a calibrated Linkam THMSG-600 heating-freezing stage.

3.2. Loading-Unloading Crossplot

The curve reflecting the sonic velocity changes with the increase in the vertical effective stress during the compaction of mudstone is called the loading curve. The unloading effect refers to the decreasing process of the vertical effective stress, and the unloading curve reflects the relationship between the vertical effective stress and the sonic velocity [18]. The loading and unloading curve of mudstone can be used not only to identify the overpressure mechanisms but also to predict pore pressure.

Disequilibrium compaction impedes the process of mudstone compaction, so the vertical effective stress and sonic velocity remain on the loading curve (Figure 4(a)). The vertical effective stress can be calculated corresponding to the given mudstone sonic velocity based on the loading equation. Overpressure caused by fluid expansion or vertical transfer can significantly reduce sonic velocity and vertical effective stress and plots on the unloading curve (Figure 4(b)). Similarly, when the unloading equation is determined, the vertical effective stress can be determined based on the sonic velocity. Furthermore, it is easy to obtain the pore pressure according to the principle of effective stress, that is, the difference between the vertical stress and vertical effective stress.

The loading curve constructed by the relationship between the sonic velocity and vertical effective stress of mudstone can be expressed as follows [18]: where is the sonic velocity of mudstone, ft/s; is the sonic velocity of the surface or mudline, generally in 5000 ft/s; is the vertical effective stress, MPa; and and are correlation coefficients, which can be determined by the regression of the vertical effective stress versus the sonic velocity from adjacent wells.

Furthermore, the pore pressure related to the loading mechanism can be calculated by the following formula: where is the pore pressure in the loading case, MPa; and is vertical stress, MPa, which is calculated by integrating the density value.

The unloading curve can be described by the following equation [18]: where is the maximum vertical effective stress, MPa; and are consistent with the parameters of the loading curve; and is the elastoplastic coefficient of mudstone.

usually equals the vertical effective stress value corresponding to the maximum sonic velocity (), and is usually equal to the velocity at the start of the velocity reversal.

indicates complete elastic deformation of mudstone, that is, the unloading curve coincides with the original loading curve. represents complete plastic deformation. The value of generally ranges from 3 to 8 with a light local change [11].

The pore pressure can be predicted for the unloading case by the following formula: where is the pore pressure in the unloading case, MPa.

3.3. 1-D Basin Modeling

The 1-D basin modeling is constructed using Schlumberger’s PetroMod v2012 to determine the overpressure mechanisms and recover the thermal history. For 1-D modeling, the informations of each unit, including the names, top depth, and thickness, are obtained from the well reports, and the ages are determined from the known ages of the key layers. The mixed lithology of each layer is determined based on cutting descriptions and (mudstone content, which is calculated from gamma-ray logging). The kerogen type, total organic content (TOC), and hydrogen index (HI) are determined based on published data in [42]. Burnham-III kinetics are selected to model hydrocarbon generation, and EASY% Ro is used to restore the maturity evolution of the source rock. The paleo-water depth (PWD), sediment-water interface temperature (SWIT), paleo-temperature gradient (PTG), and paleo-heat flow (PHF) used in the 1-D modeling refer to [43].

Regarding the reconstruction of thermal and maturity histories, the vitrinite reflectance and temperatures in holes were used for calibration purposes. Based on thermal evolution history, combined with the homogenization temperature of fluid inclusions, the charging time of natural gas is determined. Based on the burial history, temperature, and maturity evolution simulation results, the pore pressure coupled with the disequilibrium compaction and hydrocarbon generation is modeled. The final modeled pressure is constrained by the measured pressure data and mud weights.

4. Overpressure Development within the Reservoir

The screened measured pressure from the MDT and DST and the equivalent pressure are used to analyze the pore pressure characteristics in different reservoirs in the LD-B Wellblock. Overpressure means that the pressure coefficient (ratio of pore pressure and hydrostatic pressure) exceeds 1.2 in this paper.

All the measured pressures in the LD-B Wellblock are overpressures, and the pressure coefficient is between 1.2 and 2.3. The pressure distribution is extremely complicated, the pressure coefficient does not increase with burial depth, and the pressures at different depths and layers are quite different (Figure 5).

Mild-moderate magnitude overpressure developed in N1h1 and N1h2 above 3950 m, the pressure coefficient ranged from 1.29 to 1.60, and the excess pressure was approximately 10.65-23.46 MPa. For depths of 3950-4500 m, high-magnitude overpressure () was observed in the N1h and N1m formations. The pressure coefficient suddenly increased to 2.05-2.27, which is basically close to the fracture pressure gradient, and the excess pressure was between 41.15 and 51.84 MPa. The pressure coefficient of N1s in well B-11 at about 4700 m decreased back to 1.72, and the excess pressure was 32.78 MPa (Figure 5).

The lateral and vertical features of pore pressure in the LD-B Wellblock are analyzed based on the measured pressures and mud weight equivalent pressures. Overall, the overpressure in the LD-B Wellblock began to develop from the N1y formation, and the overpressure top was approximately 2000 m (Figure 6). The pore pressure varies significantly in different wells vertically. The pore pressure in well B-11 shows a gradual increase, and there is no sudden pressure change. The measured pressure is close to the mud weight equivalent pressure. However, the pore pressure distribution in well B-7 is quite different, the pressure coefficient gradually increases to 1.6 in N1h1, and it rapidly increases to about 2.27 in N1h2 within only 200 meters. The sudden increase in pore pressure is clearly reflected in the change in the mud weight (Figure 6).

5. Logging Response to the Presence of Overpressure

5.1. Sonic and Density Response to Overpressure

The comparison of sonic transit time (AC), which reflects transport properties, with density (DEN), which reflects bulk properties, is helpful for identifying the existence and origin of overpressure. The disequilibrium compaction overpressure shows that both AC and DEN deviate from the normal compaction trend (NCT), while fluid expansion is distinguished by more reversals of AC than that of DEN [20].

Well C-1 is taken as an example to discuss the mudstone compaction response to overpressure. The completion depth of well C-1 is about 5026 m, which corresponds to N1m. The measured pressure data are absent in this well, but continuous mud weights are available. The mud weight equivalent pressure shows that the pore pressure is normal above 3000 m. The overpressure gradually increases with depth over 3000 m, and the mud weight reaches 1.8-2.0 in the deep layer (>4000 m). The sonic transit time and density of mudstone in the normal pressure section above 3000 m have a good trend with increasing depth that conforms to the normal compaction trend (Figure 7).

Compared with the normal compaction trend, the sonic transit time of mudstone in the overpressure section deeper than 3000 m shows abnormally high values, but the density shows only a small deviation. Both the sonic transit time and the density deviate from the normal compaction trend for the overpressure intervals with burial depths exceeding 3500 m. At the top of the overpressure, there is an obvious response to the sonic transit time, but the density response is relatively lagging. The ratio of smectite to illite in the mixed-layer illite/smectite in mudstone decreases gradually within this depth [44], which indicates the chemical compaction. Chemical compaction of mudstone increases the density of mudstone and weakens the logging response to low-magnitude overpressure to a certain extent [45]. When the overpressure is much higher at deeper depths, the density begins to show anomalies, but the amplitude is lower than the sonic transit time (Figure 7). In short, compared with the normal pressure section, the deviation of AC and DEN from the normal compaction trend can be used as an index of overpressure in the LD-B Wellblock.

5.2. Vertical Effective Stress-Logging Properties

The loading-unloading curves constructed by the vertical effective stress and logging response can be used to effectively distinguish the overpressure mechanisms in sandstone. If the overpressure points fall on the loading curve, it indicates that the overpressure is mainly generated by disequilibrium compaction of mudstone. If the overpressure points deviate from the loading curve, it can be considered the unloading overpressure mechanism, including aquathermal pressurization, hydrocarbon generation, and overpressure transfer [20, 46].

The rocks will enter the brittle or even plastic stage in the HTHP environment and possibly change the loading and unloading behaviors [47]. According to the previous studies, the critical condition from the brittle to brittle-ductile transition of the Miocene formations is 32 MPa, the corresponding depth is approximately 2500 m, the critical stress from the brittle-ductile transition to plasticity is 123 MPa, and the depth is about 8000 m [48]. Therefore, almost all the Miocene formations are in the brittle-ductile stage, so the difference in the loading and unloading caused by different deformation behaviors may not be significant in the Ledong Slope.

The vertical stress is calculated by integrating density logging, and the vertical effective stress equals the vertical stress minus the pore pressure. The sonic transit time and density are taken as the average value of the mudstones within the interval of the measured pressure data. The normal pressure data of well C-1 above 3000 m and the LD-A Wellblock are used to establish the loading curve, which is used as a datum baseline to identify possible overpressure mechanisms in the LD-B Wellblock.

The loading-unloading crossplots of the LD-B Wellblock suggest that 10 of the 51 overpressure points plot on the loading curves, and the pressure coefficient is about 1.2-1.6 (Figure 8). The results indicate that the overpressures in sandstone were transmitted from the adjacent mudstone, which is mainly formed by disequilibrium compaction. Notably, the points with pressure coefficients close to 1.5-1.6 begin to deviate from the loading curve. A key understanding of this study is that all high-magnitude overpressure (pressure coefficient greater than 1.6) points in the LD-B Wellblock deviate the loading curve. With the increase in the pressure coefficient, the corresponding vertical effective stress decreases significantly, and the logging response falls on the unloading curve, while the points on the unloading curve express little or no changes in density, indicating the typical characteristics of unloading overpressure (Figure 8).

6. Discussions

6.1. Mechanism of Overpressure Generation

The overpressure sections in LD-B Wellblock show both sonic and density logging response anomalies (Figure 7), indicating the contribution of disequilibrium compaction to overpressure. This agrees with the understanding of Luo et al. [15] and Lei et al. [30], and it is also consistent with the results that some overpressure points fall on the loading limb in loading-unloading curves (Figure 8).

Disequilibrium compaction is the most common overpressure mechanism for thick mudstone-dominated successions in continuous rapid burial young basins, which is favored as the mechanism to explain overpressure in a number of basins, including the Yinggehai Basin, Gulf Coast, and North Sea [15, 23, 24]. The geological conditions of the sedimentation in the LD-B Wellblock lead to the occurrence of disequilibrium compaction. The mudstone contents of N1h and N1m in the LD-B Wellblock are approximately 60~80%, and the deposition rates are greater than 200 m/Ma, showing the rapid deposition of mud-rich layers. The mudstone content of the overlying N1y is more than 85%, the deposition rate reaches 1000 m/Ma, and the deposition rate of Qld is also over than 500 m/Ma. Such high sedimentation and burial rate combined with mud-dominated sequences are beneficial to the development of disequilibrium compaction. The numerical simulation results confirm that the pore pressure coefficient can reach approximately 1.5-1.6 due to disequilibrium compaction of mudstone [15, 43], which is consistent with the results that the points with pressure coefficients smaller than 1.6 plot on the loading curve (Figure 8).

The amplitude of the sonic transit time deviating from the normal compaction trend in the deep high overpressure section is much larger than that of the density, and the points with pressure coefficients greater than 1.6 fall on the unloading curve. The evidence shows that disequilibrium compaction alone is not enough to generate the high overpressure observed in reservoirs and indicates the contribution of unloading mechanism. In theory, there are four mechanisms leading to unloading effect: aquathermal pressurization, clay diagenesis, gas generation, and overpressure transfer. This paper analyzes these mechanisms in detail and evaluates the possible efficiency of each mechanism through the analysis of regional geology.

Due to the higher temperature gradient, the temperature reaches 180°C at a depth of 4000 m in the LD-B Wellblock, and aquathermal pressurization is considered to be one of the causes of overpressure [36, 49]. However, quantitative simulation results confirmed that even the permeabilities of rocks were down to , and the aquathermal pressuring was negligible [50]. In fact, this permeability is much less than the permeabilities of real shales () [50, 51]. In addition, a high temperature will reduce the viscosity of pore fluid, thereby promoting the dissipation of pore fluid [50]. Therefore, we conclude that significant aquathermal pressurization does not occur in the LD-B Wellblock.

The diagenetic transformation from smectite to illite has been thought to generate overpressure by releasing bound water or load transfer mechanisms [1, 46, 52]. It has been proven that clay diagenesis can generate very slight overpressure through the fluid expansion process [45]. The load transfer process related to clay diagenesis may generate significant overpressure [46, 52], especially when the temperature is higher than 90°C. And the load transfer has been demonstrated in the Gulf of Mexico and North Sea Basin, with suggestions that overpressure of as much as 14 MPa can be generated through this process [25, 26, 52, 53]. However, the biggest difference between the Yinggehai Basin and Gulf of Mexico and North Sea Basin is that the Gulf of Mexico and North Sea Basin are anomalously smectite-rich system [52, 53]. Small volumes of smectite in mudstones in the LD-B Wellblock have been reported [44], suggesting that the influence of the load transfer effect on overpressure in the LD-B Wellblock is not significant due to the limited transformation of smectite.

Kerogen or oil-to-gas cracking leads to a substantial increase in fluid volume, which may theoretically result in overpressure close to lithostatic pressure [3, 4]. The gas generation has been proved to be an effective pressurization mechanism in Gulf of Mexico and North Sea Basin [1, 23, 24], while it is not in the Yinggehai Basin. The average total organic carbon (TOC) of mudstone in N1y and N1h is less than 0.6%, and the TOC of the N1m and N1s source rocks is about 1.55%. The pore pressure of well C-1 was simulated to quantitatively evaluate the hydrocarbon generation pressurization effect. For Miocene source rocks with burial depths less than 5000 m, liquid hydrocarbon is the dominated product, and gas production is limited. Moreover, oil generation has little effect on overpressure, so the effect of hydrocarbon generation is not an important overpressure mechanism [15]. Obvious gas generation occurs below a depth of 5000 m, and overpressure caused by gas generation in the N1s source rock is significant, the ratio can exceed 50% in ultradeep source rocks [15, 43].

Even considering the possible self- and neighbor-source overpressure, it is not sufficient to explain some extreme overpressures phenomenon in permeability formations. Thus, some other more effective mechanisms, such as overpressure transfer, must be considered. A common concept used to explain the presence of overpressure within sandstones is that of the centroid, in which lateral pressure from deep overpressured formations transfer through an inclined sandbody to shallow intervals, resulting in a higher pressure in the latter than expected [13, 14]. The presence of faults/fractures within the overpressure system has the potential to transfer overpressure into shallower formation via enhanced permeability through fault/fracture planes [7, 15, 17, 54]. The most likely mechanism of high-magnitude overpressure in Miocene reservoirs in the LD-B Wellblock is overpressure vertical transfer along faults/fractures. Previous studies have proposed allo-source overpressure to explain overpressure generation in the Yinggehai Basin. The vertical transfer phenomenon was once thought to occur mainly in the diapir zone of the Yinggehai Basin due to the development of faults, which is easily caused the upward flow of overpressure fluid [28, 29]. Overpressure is considered to be less susceptible to vertical transfer in nondiapir zones because the faults generally do not develop [28, 31]. Researchers have suggested that the overpressure is mainly caused by disequilibrium compaction and hydrocarbon generation considering the LD-B Wellblock belonging to nondiapir zone [36, 49]. However, the evidence in this study shows that allo-source overpressure caused by vertical transfer in the nondiapir zone is still an extremely important overpressure mechanism. In fact, there is sufficient geological evidence to support the viewpoint that overpressure in Miocene reservoirs was affected by the transfer process of ultradeep overpressure fluid along faults/fractures.

6.2. Evidence for Overpressure Generation by Vertical Transfer

The pore pressure in Miocene reservoirs in the LD-B Wellblock shows abruptly transition from moderate to high-magnitude overpressure, and the vertical effective stress-logging crossplots show a typical unloading response. According to the above analysis, the ratio of overpressure caused by hydrocarbon generation is limited. Therefore, the abrupt pressure transition zone and unloading effect can only be explained by the overpressure vertical transfer. We will carry out a detailed geological analysis to discuss the source of deep overpressure fluid, the transfer pathway of allo-source overpressure, and the mode of the overpressure transfer process.

6.2.1. The Fluid Source of Overpressure Transfer

Since the high-magnitude overpressure belongs to allo-source overpressure in reservoirs in the LD-B Wellblock, the source of overpressure coming from deeper away from the Miocene reservoirs. The presence of high-magnitude overpressure in the deep system is prerequisites for this mechanism. The overpressure transfer process in permeable formations is accompanied by fluid migration. Therefore, determining the source of natural gas also indicates clarifying the source of overpressure fluid.

The carbon isotope composition of natural gas is a good indicator to identify its genetic type and source [40]. The δ13C1 of the gas sample from the N1h reservoir in the LD-B Wellblock is about -34.04‰ ~ -29.44‰, and the value of δ13C2 is about -25.88‰ ~ -23.13‰. The δ13C1 and δ13C2 of the gas sample from the N1m gas pool are -29.24‰ ~ -27.40‰ and -23.13‰ ~ -20.26‰, respectively. Both the δ13C2 values of natural gas in N1h and N1m are greater than 28‰, indicating typical characteristics of natural gas generated from coal-bearing source rocks. The δ13C1-δ13C2 diagram of natural gas shows that almost all gas samples from the LD-B Wellblock plot in the coal-type gas area (Figure 9(a)).

The mudstones of E1l and N1s are two sets of source rocks in the Yinggehai Basin and show greatly different carbon isotope compositions. The δ13C2 value of kerogen in the N1s source rock is obviously heavier, mainly ranging from -25‰ to -21‰, with an average of -23.8‰. It is about -30‰ ~ -26‰ with an average of -27.2‰ for kerogen in the E1l source rock [55]. Available data show that the natural gas in the LD-B Wellblock is more closely related to the N1s source rocks (Figure 9(b)).

It has been confirmed that the overpressure fluid mainly came from the N1s source rock, and a 1-D model is established to simulate the overpressure generation of N1s. The pore pressure of the two wells is restored by coupling disequilibrium compaction and hydrocarbon generation. The modeled pore pressure in the N1y and overlying formations is basically consistent with the mud weight equivalent pressure, and the pressure in the permeable formation is similar to that in the adjacent mudstone (Figure 10). The measured pressure in N1h and N1m is larger than that in the simulated result, indicating that the pore pressure is affected by overpressure transfer. The geological conditions of N1s source rocks are benefit to disequilibrium compaction and substantial gas generation potential, which lead to a significant increase in the pore pressure, and the pressure coefficient of N1s reaches about 1.9-2.0. If the overpressure is affected by vertical transfer, the pressure gradient in sandstones equals to the hydrostatic pressure gradient. The estimated pore pressure in N1s according to the measured pressure in shallow reservoirs is basically close to or slightly lower than the modeled result (Figure 10), which further confirms that ultradeep buried source rocks have the potential for overpressure transfer.

6.2.2. The Pathway for Overpressure Transfer

Studies on faults in nondiapir zones are very preliminary for a long time due to the limited resolution of seismic data and unobvious fault throw. It is considered that there is a lack of pathways required to form allo-source overpressures in the Ledong Slope. Therefore, the migration pathway of gas from the ultradeep buried Miocene source rocks to the shallower reservoir is a key factor for the formation of allo-source overpressures.

In recent years, three-dimensional seismic coherence and curvature analysis technologies have been applied to obtain a new understanding in the identification of the high-pressure diapir system and the fracture migration system in the Yinggehai Basin [56]. The seismic profile shows that there are a series of hidden faults/microfractures with little throw in the Miocene formations that are densely distributed in nearly E-W direction (Figure 11(a)), and the faults/microfractures have also been fully verified in the 3D seismic profiles [35, 55, 56]. The hidden faults or microfractures are mainly developed at N1m, extending upward to the bottom of N1h (Figures 11(b) and 11(c)).

Although these hidden faults/microfractures stopped extending upward after 10.5 Ma (the time of the base of N1h), due to the low friction coefficient of existing faults, the faults/microfractures will open when the pore pressure exceeds the failure threshold of incohesive faults [57]. The opening hidden faults/microfractures connect the reservoirs and source rocks, and the deep overpressure natural gas can episodically migrate upward to the reservoirs. The migration and accumulation process of gas will inevitably lead to the vertical transfer of overpressure, resulting in high overpressure in shallower reservoirs.

The overpressure vertical transfer induces the opening faults to close due to the overpressure decreases in the deep reservoirs. The low-fluid diffusion rate of mudstone maintains a higher overpressure than surrounding sandstones [9, 15]. Therefore, fluid in the deep buried system continuously flows from low-permeability mudstones to permeable sandstones during and after fault opening, which leads to an increase in overpressure in sandstones, and eventually facilitates the fault opening for next time. Deeply-buried source rocks generally have enough overpressure fluid to maintain episodic fault opening, which ensures the preservation of high overpressure in shallow formations.

6.2.3. Preservation of Vertically Transferred Overpressure

The overpressuring effects in permeable sandstone are generally instantaneous phenomena, and the overpressure dissipates quickly unless effectively sealed. Generally, the later the faults open, the more effective the sealing, and the longer time for transferred overpressure preservation in sandstone.

The time of source rock maturity in the Ledong Slope is often too late, as is the natural gas charging period. The results of the kinetics of gas generation and carbon isotope fractionation suggested that the gas reservoir in the LD-B Wellblock was formed since 1.2 ~ 0.1 Ma [37, 55]. The homogenization temperatures of aqueous inclusions formed with simultaneous gas-bearing inclusions are mainly distributed between 155~165°C and 170~185°C in the N1h and N1m reservoirs [41]. The formation age of natural gas pools can be obtained by coupling the homogenization temperatures of the inclusions with the burial-thermal history. The result shows that there is one phase of continuous gas charging in the LD-B Wellblock, and the charging time is approximately1.50~0.2 Ma (Figure 12).

The fault opening is current or geologically recent, which means that the transfer overpressure was formed at 1.50 Ma or even later. The channel sandbodies in the LD-B gas reservoir are deeply buried, the reservoir properties are poor, and the sandbodies are surrounded by overpressure mudstone. The large amount of gas in the reservoir further decreases the rock permeability due to the coexistence of two phases. These comprehensive effects lead to the preservation of transfer overpressure. It has been documented that mudstone of several tens of meters can maintain overpressure in sandstones for more than 1 Ma and that such sealing ability is common in basins [51, 54].

6.2.4. Model for the Overpressure Generation in the Ledong Slope

This study recognizes that the overpressure caused by disequilibrium compaction can reach about 1.5-1.6 in N1h and N1m in the Ledong Slope, and the overpressure vertical transfer accounts for most of the observed high-magnitude overpressure. Based on this knowledge, a model of overpressure generation in the LD-B Wellblock can be established and used as the basis for overpressure identification and pressure prediction in the future.

The N1s source rocks have the potential to develop high-magnitude overpressure with a mudstone content greater than 75%, deposition rate exceeding 280 m/Ma, organic matter content greater than 1.55%, maturity higher than 1.3 ~ 2.0%, and the pressure coefficient reaching 1.9-2.0. When the opening faults/microfractures connect source rocks and reservoirs, overpressure is transferred accompanied by the migration of natural gas along the faults/fractures, resulting in high overpressure in the channel sandbodies of the shallow N1h and N1m formations (Figure 13). Faults/microfractures are rarely developed in the Miocene and Quaternary formations, so the overpressure is confined to the fault-affected N1h (Figure 13). The gas-bearing channel sandbodies are isolated in LD-B gas reservoirs, indicating well lateral sealing conditions. The thick and widely distributed mudstone above the gas reservoir is not affected by faults/fractures and exhibits good vertical sealing. Therefore, overpressure transferred by faults/fractures can be preserved for a long time in ultralow permeability and well-sealed reservoirs.

Overpressure sequences are usually considered undesirable targets for oil and gas exploration because of the increased risk of caprock or fault seal failure [58]. Hydraulic fracturing occurs when the overpressure exceeds the fracture pressure of the cap rock, which leads to fluid and overpressure dissipation. Notably, the fracture pressure gradient increases significantly with increasing overpressure (Figure 5), a phenomenon known as pore pressure-stress coupling [59]. The pore pressure-stress coupling in turn increases the ability of the caprock to maintain the maximum overpressure [38, 59].

6.3. Pore Pressure Prediction from Log Data

Understanding the causes of overpressure is the basis for accurate pore pressure prediction; different overpressure-generating mechanisms affect the responses of rock properties to overpressure. The commonly used pressure prediction methods include the equivalent depth, Eaton, and Bowers methods, all of which are carried out based on the sonic transit time. The equivalent depth method predicts overpressure depending on the abnormally high value of porosity, and it is assumed that the overpressure is mainly generated by disequilibrium compaction of mudstone. As an empirical relationship, the Eaton method does not truly indicate the cause of overpressure, although the pressure prediction results can be modified by adjusting the Eaton index. When the loading and unloading curves are established, the Bowers method systematically considers the mechanism of disequilibrium compaction and unloading overpressure [11, 12].

Due to the overpressure in the LD-B Wellblock is a combination of the results of multiple mechanisms, including disequilibrium compaction and overpressure vertical transfer, the Bowers method is used to predict the pore pressure. The sonic transit time is converted into sonic velocity, and the vertical effective stress-sonic velocity relationship is constructed for the LD-B Wellblock (Figure 14). The key parameters for the loading-unloading equations can be determined as follows: , , , , and .

The calculation equations of pore pressure in the loading and unloading cases are and where and are the pore pressure in loading and unloading case, respectively, MPa.

The pore pressure prediction results show that the shallower intervals in well B-1 are mainly normal pressure, the pressure equation under loading conditions is adopted in the prediction, and the predicted results also show normal values. The measured pressure at about 4000 m shows obvious unloading characteristics, and the predicted pressure under unloading conditions is basically consistent with the measured pressure (Figure 15(a)). Moreover, the pore pressure calculated using the equivalent depth method is also superimposed on the plot; for the loading process at shallow depths, the two results are similar. The pressures predicted by the equivalent depth method can approximately reflect the mudstone pressure, but seriously underestimates the reservoir pressure in the target sections (Figure 15(a)). This implies that the transfer overpressure often does not have the corresponding logging anomaly amplitude. The Bowers method can be applied to pore pressure prediction in wells D-1 and E-7, which are nearby the LD-B Wellblock, and more reliable pressure prediction results can also be obtained (Figures 15(b) and 15(c)).

It is worth noting that the pressure predicted by the equivalent depth method reflects continuous pressure changes in mudstone. However, the predicted pressure of sandstone only reflects the pressure in the sandbodies, which are affected by overpressure vertical transfer, rather than the pressure distribution trend of the whole strata. The pressure in any interval should be between the predicted pressure of mudstone and the pressure of sandstone. To determine the pressure in sandstone, the relationship between sandbodies and faults/fractures should be confirmed. If the sandstone does not have the conditions necessary to connect with the deep overpressure systems, the pressure may be basically close to the mudstone pressure. When the sandstone is connected by faults with the higher overpressure sandstone, the pressure should refer to as the prediction pressure by the Bowers method.

7. Conclusions

(1)The disequilibrium compaction of mudstone and fluid expansion mechanisms are insufficient to explain the high-magnitude overpressure in reservoirs in the LD-B Wellblock in the Yinggehai Basin. All the high-magnitude overpressure points deviate from the loading curve and show a typical response of the unloading mechanism. Allo-source overpressure caused by vertical transfer in the nondiapir zone is still an extremely important overpressure mechanism(2)The ultradeep buried Miocene source rocks developed high-amplitude overpressure, which ensured a sufficient source for overpressure transfer fluid. The hidden faults/microfractures connecting source rocks and reservoirs provided a possible pathway for overpressure vertical transfer. The formation time of allo-source overpressure was no earlier than 1.5 Ma, allowed for overpressure preservation in well-sealed lenticular sandbodies(3)The allo-source overpressure related to vertical transfer does not always have enough logging response amplitudes corresponding to the overpressure values. The pore pressure calculated by the Bowers method is in good agreement with the measured pressure, especially the vertically transferred overpressure can be accurately predicted by the Bowers method(4)Recognitions of overpressure mechanisms and pressure prediction provide new insight into high-magnitude overpressure generation and natural gas formation, which is especially meaningful for predrilling pressure prediction in nondiapir zones of the Yinggehai Basin, and it helps to reduce drilling risks

Data Availability

The main data used to support the findings of this study are included within the article, and the others are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no competing interests.


This study has been funded by the Strategic Priority Research Program of the Chinese Academy of Sciences (XDA14010202), National Science and Technology Major Project (2017ZX05008-004), and National Natural Science Foundation of China (U20B2014). We are thankful to the Zhanjiang Branch of the CNOOC for providing the original data. Editors and reviewers are appreciated for their critical comments and suggestions, which significantly improve the quality of this work.