Phase Behavior Identification and Formation Mechanisms of the BZ19-6 Condensate Gas Reservoir in the Deep Bozhong Sag, Bohai Bay Basin, Eastern China

Significant developments have been observed in recent years, in the field of deep part exploration in the Bozhong Sag, Bohai Bay Basin in eastern China. The BZ19-6 large condensate gas field, the largest gas field in the Bohai Bay Basin, was discovered for the first time in a typical oil-type basin. The proven oil and gas geological reserves in the deeply buried hills of the Archean metamorphic rocks, amount to approximately 
 
 3
 ×
 
 
 10
 
 
 8
 
 
 
 tons of oil equivalent. However, the phase behavior and genetic mechanisms of hydrocarbon fluids are still unclear. In this study, the phase diagram identification method and various empirical statistical methods, such as the block diagram method, 
 
 
 
 φ
 
 
 1
 
 
 
 parameter method, rank number method, and 
 
 Z
 
 -factor method were implemented to comprehensively identify the phase behavior types of the BZ19-6 condensate gas reservoir. The genetic mechanism of the BZ19-6 condensate gas reservoir was investigated in detail through analyses of physical properties of the fluid at high temperatures and pressures, organic geochemical characteristics of condensate oil and gas, and regional tectonic background. Consequently, this study is shown as follows: (1) The BZ19-6 condensate gas reservoir is a part of a secondary condensate gas reservoir with oil rings, formed synthetically since the Neogene period due to multiple factors, such as retrograde evaporation from deep burial and high temperature, inorganic CO2 charging from the deep mantle, and late natural gas invasion. (2) The hydrocarbon accumulation process of the BZ19-6 condensate gas reservoir is as follows: First, a large amount of oil is accumulated at the end of the lower Minghuazhen deposition (5 Ma BP); second, a large amount of natural gas is generated in the deep-source kitchen and mantle-derived inorganic CO2 charged into the early formed oil reservoirs at the end of the upper Minghuazhen deposition (2 Ma BP). As a result, the content of gaseous hydrocarbons in the hydrocarbon system of the reservoir increased, which led to large amounts of liquid hydrocarbons dissolved in gaseous hydrocarbons and significantly reduced the critical temperature of the hydrocarbon system. Therefore, existing secondary condensate gas reservoirs are formed when the critical temperature is lower than the formation temperature and it enters the critical condensate temperature range.


Introduction
A series of deep of strata and deep water and unconventional areas represent the global development trends of oil and gas exploration [1][2][3]. Deep oil and gas are the most suitable replacement fields for hydrocarbon exploration in China and the world [2,3]. The Bohai Bay Basin is recognized as one of the two super petroliferous basins in China and is a typical oil-type basin [4,5]. In recent years, a strategic breakthrough has been achieved pertaining to the deep gas exploration in the Bohai Bay Basin. The proven oil and gas geological reserves of the first discovered BZ19-6 large condensate gas field in the deeply buried hills of the Archean metamorphic rocks, amounted to an approximate equivalent of 3 × 10 8 tons of oil (according to http://www.china.com.cn, 2020-1- 10 14 : 16). Here, the proven gas geological reserve and the proven condensate oil geological reserve amounted to more than 1 × 10 8 m 3 each, revealing immense exploration potential of deep gas [4,6]. The existence of large deep gas reservoirs in the Bozhong Sag can be attributed to certain favorable conditions. First, the Bozhong Sag is the most likely source kitchen to form large gas fields, with a high hydrocarbon potential and thermal maturity [5,7]. Second, the Bozhong Sag exhibits favorable petroleum geological conditions with several sets of reservoir-cap assemblages [8,9]. In addition, the regional thick mudstone of the Dongying Formation is favorable for deep hydrocarbon accumulation [4-6, 8, 9].
It was found that the BZ19-6 deep buried hill gas reservoir is a condensate gas reservoir. A condensate gas reservoir is a special hydrocarbon deposit located between oil and gas reservoirs and is more valuable compared with conventional oil and gas reservoirs. Condensate gas exhibits the characteristics of complex genesis, complex thermodynamics, and multiple phase behaviors [10][11][12][13], of which phase behavior and formation mechanism are of great significance in the evaluation of deep oil and gas resources. The study on the phase behavior of condensate gas reservoirs began in the 1930s. A large number of gas-liquid state experiments have expanded our understanding of condensate oil and gas and have promoted further research and development of condensate gas reservoirs [14,15]. A majority of the condensate oil and gas, which is the product of organic matter formed in the higher maturity stage (R o > 1:3%), higher temperature, and higher pressure environment, are initially generated in the late hypogenic stage. In addition, condensate gas is directly generated by the resin component of type III kerogen in the immature stage (0:4% < R o < 0:6%) [16][17][18]. Secondly, the chemical composition and temperature-pressure system of hydrocarbon systems are the key factors controlling the hydrocarbon phase behavior [17][18][19][20][21][22]; multiple sets of source rocks and multiple stages of hydrocarbon migration along with secondary alteration after hydrocarbon accumulation lead to diversity in the hydrocarbon phase behavior and complexity in its genesis [23]. The condensate gas reservoirs can be classified as primary and secondary, among which the secondary type can be subdivided into crude oil cracking, migration fractionation, and gas invasion and enrichment types.
The BZ19-6 large condensate gas reservoir possesses the following characteristics: humic-sapropelic and kerogencracking gases are the origin gas, the gas field has a large gas thickness and high gas column, and the accumulation process began with oil, followed by gas [8]. The hydrocarbon accumulation process of "early oil and late gas" was characterized by near-source early oil and late gas continuous charging, in which the crude oil and the natural gas was charged at 12-5.1 Ma BP and 5.1-0 Ma BP, respectively [24]. Phase-controlled miscibility is the main formation mechanism of condensate gas in the BZ19-6 structure; when the local temperature and pressure rise, the early crude oil dissolves in the late charging natural gas and forms the present condensate gas reservoir [25,26]. However, phasecontrolled miscibility cannot completely explain the genesis and formation processes of the BZ19-6 condensate gas reser-voirs. Moreover, the deep oil and gas phase behavior affects the migration mode and accumulation mechanism and influences the direction of deep oil and gas exploration, restricting the evaluation of oil and gas resources in deep condensate gas reservoirs. Based on the composition of the BZ19-6 condensate gas and the oil test results, the phase behavior types of the condensate gas reservoir were comprehensively identified using various methods. The formation mechanism of the condensate gas reservoir and its hydrocarbon accumulation process were systematically analyzed based on previous studies of sedimentary and tectonic evolution characteristics and key hydrocarbon accumulation elements in this area. The results of our study can potentially guide future resource evaluation of deep condensate gas reservoirs in the Bozhong Sag and development of deep hydrocarbon accumulation in the geological theory.

Petroleum Geological Background
The BZ19-6 condensate gas field is an anticlinal uplift located on the deep structural ridge in the southwest of the Bozhong Sag, adjacent to the Chengbei low uplift in the west, the Bonan low uplift in the southeast, the Shaleitian uplift in the north, and the main subsag in the northeast (Figure 1). Drilling revealed that the BZ19-6 and CFD18 structures in the north are both direct outcropping areas of the Archean. The strata in the BZ19-6 structure from new to old includes the Neogene Minghuazhen Formation (Fm.) (N 2 m), Guantao Fm. (N 1 g), Paleogene Dongying Fm. (E 3 d), Shahejie Fm. (E 2 s), Kongdian Fm. (E 1-2 k), and the Archaeozoic Fm. (Ar). The N 2 m is a mixed deposition of mudstone and sandstone, whereas the N 1 g primarily comprises sandstone; both are Neogene reservoirs. Both E 3 d and E 2 s primarily comprise dark mudstone and a small amount of sandstone, which are the major source rocks. E 1-2 k is primarily glutenite and reservoirs. The Ar is primarily a metamorphic rock buried hill, which is the main reservoir. From top to bottom, the strata can be divided into three sets of hydrocarbon accumulation assemblages: upper Neogene reservoir assemblage, middle Paleogene reservoir assemblage, and lower buried hill reservoir assemblage [5].
According to the structural characteristics, the BZ19-6 low buried hill is divided into south and north blocks in the plane, and strike-slip and extensional faults cut each other, forming several complex fault-block traps with anticline and broken noses. Vertically, it includes the single-layer structures of the Archean buried hill and the double-layer structures in the glutenite of the Kongdian Formation, with the highest point of the trap being between 3825 and 5475 m [8].
The metamorphic rocks of the Archean low buried hill and glutenite of the Kongdian Fm. overlying the low buried hill are both the primary gas-bearing strata of the BZ19-6 condensate gas reservoir. The lithology of the Archean reservoirs is primarily monzonitic gneiss, plagioclase gneiss, and mixed gneiss. The reservoir is dominated by fractures, along with the presence of dissolution pores. The porosity ranges from 2.7% to 12.8% and the permeability ranges from 0.01 to 11.81 mD [4]. Since 2016, several wells have been drilled 2 Geofluids into the BZ19-6 structural belt and have exhibited good productivity. The highest daily gas production in the Kongdian Fm. and Archean region is >30 × 10 4 m 3 (600 m 3 /d oil equivalent), the daily oil production is >300 m 3 , and the gas-oil ratio (GOR) is 1023-1539 m 3 /m 3 ( Table 1). The average temperature of the Archean buried hill is 150.8°C, and the geothermal gradient is 3.3°C/100 m. The average pressure is 47.46 MPa, and the pressure coefficient is 1.15-1.26; this implies that the Archean buried hill is a weak overpressure condensate gas field.

Methods and Results of Fluid Phase Identification
The accurate identification of the condensate gas reservoir phase behavior type is related to the estimation of oil and gas reserves, optimization of the development mode, and evaluation of the production performance. The identification methods include the phase diagram and empirical statistical discrimination methods, among which the pressure-volume-temperature (PVT) experimental phase diagram discrimination method is the most reliable; however, the effort and cost of field sampling and experimental analysis associated with this method are relatively high. The empirical statistical discriminant method requires less data and is simple and easy to operate; however, it also has certain limitations. Therefore, a comprehensive analysis using both the methods should be considered. Several PVT experiments have been implemented for verification of existing experimental results; based on these PVT experiments, the phase pattern discriminant method along with several empirical statistical discriminant methods were used to comprehensively identify the phase behavior types of the BZ19-6 condensate gas reservoir. Moreover, PVTsim20 software was used to simulate the change characteristics of the hydrocarbon phase behavior with an increase in the CH 4 and CO 2 content. The PVTsim20 software was used during numerical simulation; the P-R state equation, introducing Peneloux volume correction integrated with the thermodynamic theories and gas-liquid phase equilibrium equations, was selected. Furthermore, the parameters including dew-point pressure, dew-point pressure, constant composition expansion (P-V relation), and phase envelope were determined by phase equilibrium calculation of the hydrocarbon system. Meanwhile, the results of the theoretical calculation were verified with the measured results to confirm the accuracy of the simulation results.

Phase Diagram Identification
3.1.1. Fluid Properties of the Gas Reservoir. The methane content in the natural gas of the BZ19-6 condensate gas reservoir ranged from 70.85% to 78.25%, with an average of 75.63%. The nitrogen content was low, ranging from 0.12% to    3 Geofluids 0.15%. The CO 2 content ranged from 9.35% to 15.76%, with an average of 11.24% ( Table 2). The average surface condensate density was 0.798 g/cm 3 (at 20°C), the viscosity ranged from 1.326 to 2.438 mPa·s (at 50°C), with an average of 1.796 mPa·s, and the average contents of wax, sulfur, resin, and asphaltene were 13.7%, 0.016%, 0.836%, and 0.274%, respectively. The Pr/Ph ratio is 1.29-1.51, indicating a sedimentary environment of weak reduction and weak oxidation.
Formation fluids include gaseous hydrocarbons, liquid hydrocarbons, and nonhydrocarbons. The components of liquid hydrocarbons are particularly complex, usually considered as a whole (C 7 + ) in PVT simulation [27,28], and the critical temperature and critical pressure of the subsur-face multicomponent fluids can be further calculated [29,30]. The previous PVT experimental results of the formation fluid indicate that the BZ19-6 structural fluid is a single-gas phase underground, with C 1 +N 2 content ranging from 61.19% to 71.62%, C 2-6 +CO 2 content ranging from 21.09% to 30.90%, and the C 7+ content ranging from 6.45% to 8.06%. It is a condensate gas reservoir in the triangle diagram of the fluid type of the oil and gas reservoir ( Figure 2). The critical pressure of the formation fluid is low (29.53-46.22 MPa), with an average of 39.00 MPa. The critical temperature is low (between -22.20°C and 64.70°C), with an average of 33.75°C. However, the dew-point pressure was high (45.50-48.28 MPa) and the dew-point pressure difference was small  The P-T phase diagram is shown in Figure 3. The formation temperature isoline was in the critical condensate temperature range, and the dew-point pressure was lower than the formation pressure; the hydrocarbon fluid existed in a single-gas phase. The difference between the formation pressure and dewpoint pressure was 1.15 MPa in well BZ19-6-2Sa; the critical pressure was 46.22 MPa, critical temperature was 64.70°C, critical condensate pressure was 48.19 MPa, and critical condensate temperature was 393.40°C. In the constant volume depletion test, the maximum reverse condensate pressure was 27.62 MPa, maximum amount of condensate was 26.79%, production GOR was 1381 (m 3 /m 3 ), and condensate content in the ground was 568.41 g/m 3 (Figure 3(a)). The difference between the formation pressure and the dew-point pressure was 2.92 MPa in well BZ19-6-5; the critical pressure was 43.08 MPa, critical temperature was -9.90°C, critical condensate pressure was 48.19 MPa, and critical condensate temperature was 364.30°C. In the constant volume depletion test, the maximum reverse condensate pressure was determined to be 19.00 MPa; the maximum amount of condensate was 19.74%, production GOR was 1421 (m 3 /m 3 ), and condensate content on the ground was 594.31 g/m 3 ( Figure 3(b)). The BZ19-6 condensate gas reservoir exhibits a small ground surface pressure difference and high condensate oil content.

Empirical Statistical Identification.
Several empirical statistical discrimination methods of phase behavior types have been proposed using statistics based on physical property data, such as fluid composition, average molecular weight, density, and gas-oil ratio, which reflect the phase behavior types of oil and gas reservoirs [31]. In this study, five reliable empirical discriminant methods such as the block diagram method, average relative molecular weight, φ 1 parameter method, rank number method, and Z-factor method were used to comprehensively identify the phase behaviors of hydrocarbon fluids in the BZ19-6 structure.

Four-Parameter Block Diagram
Method. The parameters and criteria used in the four-parameter block diagram method are presented in Ref. [30]. The parameters of the Kongdian Fm. and Archean buried hill is primarily distributed in the lower-left corner of the box of the oil ring gas Volatile oil Oil condensate gas C 2 ~ C 6 + CO 2 /mol%   : Four-parameter block diagram and the types of condensate gas reservoir in the BZ19-6 structure. The four parameters are C 2+ , C 2 /C 3 , 100C 2 /(C 3 +C 4 ), and 100C 2+ /C 1 , respectively. Types of condensate gas reservoir are both in boxes of gas reservoir with oil rings or condensate gas cap reservoir, and condensate gas reservoir without oil rings. 6 Geofluids reservoir or condensate gas cap reservoir and the upper-right corner of the box of the condensate gas reservoir without oil rings, indicating that the reservoir fluid type belongs to the gas reservoir with oil rings (Figure 4).

Identification of Formation Fluid Density and Average
Relative Molecular Mass. The density of the formation fluid is directly measured by sampling; it can also be obtained by an empirical formula, in case measured data are unavailable. The average relative molecular mass and discriminant criteria are shown in "Technical requirements for determining phase behaviors of gas condensate reservoirs" (SY/T 6101-2012) [32]. The average molecular mass of the reservoir fluid in the Kongdian Fm. and the Archean buried hill in the BZ19-6 structure is 32.25-39.97; the density of the formation fluid is 0.417-0.468 g/cm 3 ( Figure 5), which is between those of the condensate gas reservoir and volatile oil reservoir. However, the reservoir is speculated to be a condensate gas reservoir with oil rings.

φ 1 Parameter, Rank Number, and Z-Factor Method.
Considering the C 5+ content criterion, C 1 /C 5+ ratio, φ 1 parameter, rank number, and Z-factor methods, the fluid phase behavior of the Kongdian Fm. and the Archean buried hill reservoir in the BZ19-6 structure was analyzed, and the oil ring characteristics were determined. Various empirical statistical discrimination methods and standards were presented in Ref. 32. The results are shown in Table 3. The results of all the methods show that BZ19-6 is a condensate gas reservoir with oil rings.

Heavy oil reservoir
Black oil reservoir Condensate gas reservoir

Properties of Underground
Fluid. The density of the underground fluid ranged from 0.4172 to 0.4678 g/cm 3 , with an average of 0.44 g/cm 3 , and an average viscosity of 6.05 ×10 -5 mPa·s; the formation temperature and pressure are higher than the critical temperature and pressure of the hydrocarbon fluid, which may be a supercritical fluid. This fluid exists above the critical temperature and critical pressure and exhibits unique physical properties. Its density and solubility are close to those of liquid and are sensitive to temperature and pressure changes. Its diffusion coefficient is close to that of gas and has an extraordinary ability to extract organic matter. There is no obvious gas and liquid interface, and the event phase effect is not obvious or does not exist [33]. These special properties are conducive to effective hydrocarbon migration and accumulation.

Discussion on Formation Mechanism of Condensate Gas Reservoirs
Condensate gas reservoirs of different genetic types in hydrocarbon-rich sags significantly differ in terms of fluid properties, phase characteristics, geochemical characteristics, and distribution rules, which can be used to distinguish condensate gas reservoirs of different genetic types. The study of natural gas and condensate oil in the condensate oil-gas sys-tem is crucial to understand the genetic mechanisms of condensate gas reservoirs. Some studies have shown that natural gas in the BZ19-6 Archean group buried hills and the overlying Paleogene Kongdian Fm. is an oil-type gas with the characteristics of different mature natural gas mixtures. There are three suites of hydrocarbon source rocks for deep oil and gas, among which the infilling from the third member of the Shahejie Fm. is the primary source. The formation of deep condensate gas reservoirs is characterized by oil-gas charging and oiltype gas formation in different periods, which have experienced early reservoir formation and destruction (from 24 Ma to 5 Ma BP), late reservoir adjustment, and condensate gas reservoir conversion (from 5 Ma BP) [8]. Retrograde evaporation is a notable characteristic of condensate gas reservoirs. As the temperature and pressure increases with burial depth, the liquid hydrocarbons in the oil and gas reservoirs gradually dissolve and evaporate, forming a condensate gas reservoir of inverse evaporation origin, with oil rings or condensate gas cap.
The burial history of BZ19-6-1 shows that the burial depth from the lower member of the Minghuazhen Fm. (5 Ma BP) increased substantially by approximately 1100 m, during the major accumulation period. The geothermal gradient in this area is approximately 3.40°C/100 m, thereby increasing the formation temperature by approximately 37.40°C during the hydrocarbon accumulation period.
The rise in the formation temperature and pressure leads to the gradual dissolution and evaporation of liquid hydrocarbon in the reservoir containing the gaseous hydrocarbon, that is, retrograde evaporation. The BZ19-6 oil and gas reservoirs were identified as a condensate gas reservoir with oil rings, which is consistent with the phase characteristics of the retrograde evaporation of genetic condensate gas reservoir. This implies that the BZ19-6 condensate gas reservoir experienced retrograde evaporation. The formation process is presented in Figure 6. The original formation temperature was less than the critical temperature; the hydrocarbon fluid existed in a single-liquid phase, and the formation temperature exceeded the critical temperature with increased burial depth. The above indicators support that the isothermal depressurization line moves downward, and the hydrocarbon fluid changes from the liquid phase to the condensate gas phase, forming a retrograde evaporation genetic condensate gas reservoir.

Gas
Invasion. Significant increase of the burial depth of the BZ19-6 structure since the Neogene resulted in the rapid increase in the maturity of the source rocks in the Bozhong Sag with the generation of a large amount of gas. The deep natural gas was charged into the early oil reservoirs along the unconformity and deep faults, forming gas invasion condensate gas reservoirs. The evidence of gas invasion in the BZ19-6 structure is as follows.   [34,35]. The R o values calculated using the two formulas ranged from 1.61% to 1.85% (Table 4). The values of δ 13 C 2 -δ 13 C 1 are also an effective index for measuring the maturity of natural gas, which decreases with maturity [36]. The values of δ 13 C 2 -δ 13 C 1 of the gas in the BZ19-6 structure lie between 11.50‰ and 12.89‰, and the corresponding R o is approximately 1.5%, indicating a primary gas stage. This indicates that the natural gas of the BZ19-6 structure is the product of the deep-source rock in the Bozhong Sag formed during the high-maturity stage.
The heavy hydrocarbon gas of the BZ19-6 structure is poor in light aromatics (benzene and toluene) and rich in cycloalkanes (methylcyclopentane, cyclohexane, and methylcyclohexane) ( Table 5). This feature may be related to evaporative fractionation in the process of gas invasion [37], which belongs to the product of the initial stage of evaporative fractionation according to Thompson's cross-plot ( Figure 7). However, Thompson's cross-plot is ineffective in validating the formation of a condensate gas by evaporation fractionation. However, regardless of the influence of organic matter type and maturity, the presence of an unusually high or abnormally low content of low-carbon aromatic hydrocarbons in the crude oil or condensate gas, could potentially indicate whether the crude oil or condensate gas experienced a secondary alteration and whether evaporation fractionation was the reason [38]. The late charging of the natural gas transforms the early oil reservoirs into condensate gas reservoirs and increases the wax content of condensate oil and paraffinicity of the heavy hydrocarbon gas, ensuring that the oil quality is higher than that of the normal condensate [13]. The average condensate density of the BZ19-6 structure is 0.798 g/cm 3 (at 20°C), which is higher than the normal condensate density (approximately 0.74 g/cm 3 ); the wax content and the paraffinicity is high, with an average of    (Tables 2 and 5), respectively. This observation confirms that the BZ19-6 reservoirs have experienced gas invasion.
(2) Characteristics of Fluid Inclusions and Asphalt in Reservoirs. Oil inclusions in the BZ19-6 condensate gas reservoir are widely developed with GOI values of up to 80%. Observation under a microscope revealed that the low-mature light oil inclusions primarily exhibited a yellow-green fluorescence color (Figure 8(a)), and the highly mature gas-rich inclusions in the deep reservoirs (E 1-2 k and Pt) exhibited a blue-white color (Figure 8(b)). The abundance of gas inclusions in the deep layer is higher than that in the shallow layer, which is also an evidence of gas invasion.
Dark brown solid asphalt can be observed in the fractures and pores of the upper reservoir of the BZ19-6 structure (Figures 8(c) and 8(d)) and is also visible under the microscope (Figures 8 (e1), 8(f1), and 8(g)); it displays a green color under ultraviolet (UV) excitation fluorescence (Figures 8(e2) and 8(f2)), which demonstrates the characteristics of oil asphalt. The equivalent vitrinite reflectance of the asphalt is only 0.9%, which indicates a gas invasion origin for the asphalt [8].
The late gas invasion leads to an increase in the light component (mainly CH 4 ) and the phase behavior of the oil and gas changes accordingly. PVTsim20 software was used to reconstruct the fluid components of the BZ19-6 reservoir and analyze the influence of a single factor on the fluid phase by controlling the variables. The changes in the oil and gas phase behavior during the gas invasion were simulated by increasing the CH 4 content in the hydrocarbon system, and the results are shown in Figure 9. With an increase in the CH 4 content (gas invasion intensity), the critical temperature   Critical point of C 1 increased by 20% Phase diagram of C 1 increased by 20% Phase diagram of C 1 increased by 10% Critical point of C 1 increased by 30% Critical point of C 1 increased by 10% Formation temperature and pressure Original critical point Phase diagram of C 1 increased by 30%  Crust-derived inorganic CO 2 Figure 10: CO 2 genesis discrimination chart in the BZ19-6 structure (after [39]). According to 3 He/ 4 He and δ 13 C CO2 data and genetic partition, CO 2 in the BZ19-6 structure belongs to mantle-derived gas.

Geofluids
hydrocarbon fluid changed from the liquid phase to the condensate gas phase, forming a gas invasion condensate gas reservoir.
The phase behavior of the condensate gas reservoirs is closely related to the gas invasion intensity. When the gas invasion intensity is high, the reservoir would be a high GOR pure condensate gas reservoir without the oil ring; when the gas invasion intensity is low, it would be a low GOR condensate gas reservoir with oil rings. Therefore, the intensity of the gas invasion affects the phase behavior of  Figure 11: Change process of phase diagram with the increase of CO 2 content in the BZ19-6 gas field. If CO 2 content is increased by 10%, 20%, 30%, 35%, and 45%, the critical temperature is decreased from 506.10°C to 448.10°C, 321.30°C, 115.7°C, 20.20°C, and 5.80°C, respectively. Therefore, the original oil reservoir easily transformed into the condensate gas reservoir.

Retrograde evaporation
Deep CO 2 and gas charging Figure 12: Schematic diagram of oil and gas phase behavior genesis from BZ19-6 condensate gas reservoir, it may be a secondary condensate gas reservoir, formed by retrograde evaporation, inorganic CO 2 charging from the deep mantle, and late natural gas invasion in order.
oil and gas. The gas invasion intensity of the BZ19-6 condensate gas reservoir is relatively weak, implying that it is a condensate gas reservoir with oil rings with an average GOR of 1269 m 3 /m 3 .

Mantle-Derived
Inorganic CO 2 Charging. The average CO 2 content of the BZ19-6 condensate gas reservoir was 11.24%, and the carbon isotope value of CO 2 ranged from -3.6‰ to -7.0‰, implying that CO 2 is an inorganic gas.
The δ 13 C CO2 and 3 He/ 4 He (Ra) values confirmed that CO 2 is an inorganic gas derived from the volcanic mantle ( Figure 10), primarily from CO 2 -rich mantle degassing [36]. The phase behavior of the hydrocarbon fluid is affected by the increase in the CO 2 content [40]. PVTsim20 software was used to simulate the change characteristics of the hydrocarbon phase behavior with increasing CO 2 content [41], as shown in Figure 10. As the CO 2 content increased, the critical temperature of hydrocarbons decreased, the critical pressure and saturation pressure increased, the critical point shifted from right to left, and the phase diagram changed to light characteristics (shrinking to the left). When the CO 2 content increased by 30%, the critical temperature decreased from 506.10°C to 115.70°C, and the original oil reservoir transformed into a condensate gas reservoir ( Figure 11). The increase in the CO 2 content is conducive to the formation of the condensate gas reservoir. The simulation results are consistent with the mercury-free visual DBR-PVT experiment by Hou et al. [42]. The technology of CO 2 injection to enhance oil recovery in the later stage of condensate gas reservoir development is a practical application based on this theory. It can increase natural gas production by at least 1.39 times compared with the water flooding mode [43,44].
Considering with the above factors, the BZ19-6 condensate gas reservoir may be a secondary condensate gas reservoir formed by retrograde evaporation, inorganic CO 2 charging from the mantle, and late natural gas invasion. The changes in the phase behavior of the hydrocarbon are shown in Figure 12. Retrograde evaporation was caused by the burial depth and temperature increase since the Neogene. Concurrently, a large amount of natural gas and mantlederived inorganic CO 2 was charged into the early formed oil reservoirs at the end of the upper Minghuazhen deposition (2 Ma BP). As a result, the content of gaseous hydrocarbons (C 1 +CO 2 ) increased, which provided a material basis for the liquid hydrocarbon dissolving in the gaseous hydrocarbon in the reservoir, and significantly reduced the critical temperature of the hydrocarbon system. When the critical temperature was below the formation temperature, it shifted into the critical condensate temperature range and formed the BZ19-6 condensate gas reservoir.

Differences in Composition and Genesis of Natural Gas.
The BZ21-22 structure is located to the east of the BZ19-6 structure ( Figure 1). In recent years, numerous natural gas reserves have been discovered in the buried hills of Ordovician carbonate rocks in the BZ21-22 structure. The average contents of C 1 , C 2-6 , C 7+ , and CO 2 in the BZ19-6 structure were 66.80%, 14.65%, 7.31%, and 24%, respectively, according to the PVT experimental results. Moreover, the natural gas of BZ19-6 was primarily a product of high-maturity source rock (R o : 1.61-1.85%). Component analysis of the natural gas of the BZ21-22 structure revealed that the content .5 % 9 C 2 /C 3 12 15 Figure 13: Genesis of the natural gas between in BZ19-6 and BZ21-22 structures (after [47]). The results suggest that the gas genesis in BZ21-22 structure belongs to NSO (nitrogen, sulfur, and oxygen organic compounds) macromolecular compounds cracked gas, where the gas in BZ19-6 structure belongs to Kerogen cracked gas. 14 Geofluids of CH 4 was 59.6%, C 2-6 was 5.34%, with a minute amount of C 7+ , and CO 2 content as high as 34.6%. Moreover, the maturity of natural gas of BZ21-22 was substantially higher than that of BZ19-6, and the R o was distributed in the range of 1.61-2.73% (Table 4). The carbon isotope value of CO 2 was -3.2‰, which was much higher than that of the BZ19-6 structure; however, it was of the same origin as the mantle-derived inorganic gas (Figure 9). The high content of CO 2 is related to several deep and large basement faults in the BZ21-22 structure. During the Neotectonic movement, the activity of deep and large faults in the basement was enhanced, converting them into dominant channels for the migration of mantle-derived inorganic CO 2 . The difference in the activity of deep and large faults may result in a difference in the content of mantle-derived inorganic CO 2 in the BZ19-6 and BZ21-22 structures [45,46]. The formation temperature of the BZ21-22 structure was as high as 173°C, which was higher than that at the same depth in the BZ19-6 structure. This may be due to the thermal anomaly caused by the strong activity of the deep mantle-derived inorganic CO 2 , which promoted the maturity of the source rocks higher than that of the BZ19-6 structure. The R o value of the Shahejie source rocks of the BZ21-22 structure was 0.1-1.6% during the main gas generation stage and that of the source rocks with burial depths >5000 m was >2.0%; the R o of the Shahejie source rocks of the BZ19-6 structure is generally <0.9%. High temperatures (more than 160°C) promoted the cracking of crude oil into gas. The gas genesis of the BZ21-22 structure revealed that it is a nitrogen-sulfur-oxygen (NSO)-cracked gas (Figure 13), which is different from that of the BZ19-6 structure.

Differences in Genetic
Mechanisms. Based on the tectonic location, natural gas composition, and genesis, the origin of pure gas reservoirs of BZ21-22 can be summarized. The mantle-derived inorganic CO 2 fluid activity along the basement deep fault caused the hydrocarbon system to be rich in CO 2 and substantially reduced its critical temperature.
This further induced thermal anomalies that consequently increased the maturity of the source rock during the main gas generation stage. The geothermal temperature rise caused the NSO cracking of crude oil into gas. In addition, the BZ21-22 structure was located closer to the Bozhong Sag than the BZ19-6 structure. The natural gas generated in the late stage of the hydrocarbon source rock preferentially charged into the BZ21-22 structure, and the gas invasion intensity became greater than that of the BZ19-6 structure. Therefore, BZ21-22 pure gas reservoirs were formed by the combined action of deep mantle-derived inorganic CO 2 fluid activity and late gas invasion.   (e2)). There are at least two stages for the inclusions of the BZ19-6 structure: the early stage comprises low-maturity light oil inclusions with green and yellow-green fluorescence, and the late stage includes high-maturity gas-rich inclusions with blue-white fluorescence. Light brown hydrocarbon-bearing brine inclusions and dark gray gas inclusions (Figures 14(h) and 14(i)) were observed in the deep Archean quartz grain microcracks, which contain significant late natural gas charging characteristics.

4.3.2.
Charging Period of Hydrocarbon. Hydrocarbon charge and accumulation in the Shijiutuo uplift in the northern Bozhong Sag rapidly occurred in the late stage, with a charging time not earlier than 5.1 Ma BP [48]. With this method, the oil and gas charging period of the BZ19-6 reservoir could be restored by the homogenization temperature and thermal evolution history of the hydrocarbon-containing brine inclusions. The reservoirs of the Kongdian Fm. and the Archean buried hill experienced at least two stages of oil and gas charging. The homogenization temperature of low-maturity light oil is 120-125°C, which corresponds to the end of the lower Minghuazhen Fm. (5 Ma BP). The homogenization temperature of the rich gas inclusion is 130-135°C, which corresponds to the end of the upper Minghuazhen Fm. (2 Ma BP). The main accumulation period was 5 Ma BP, with evident late rapid accumulation characteristics ( Figure 15).

Hydrocarbon Accumulation
Process. The oil-source correlation indicates that the condensate oil and gas mainly originate from the high-maturity source rock of the third section of the Shahejie Fm. (E 2 s 3 ) [9]. The BZ19-6 condensate gas reservoir is a regional overpressure mudstone enrichment model that has undergone a continuous natural gas charging process after the oil accumulates in the trap [4]. The excess pressure in mudstone is widespread, which is an important factor affecting hydrocarbon accumulation in the Bohai Bay Basin [49]. Based on the study of the BZ19-6 condensate gas reservoir genetic mechanism, the period of oil and gas charging, and comparison with the gas reservoir in the BZ21-22 structure, the hydrocarbon accumulation process of the gas reservoir of the BZ19-6 structure was analyzed. At the end of the Paleogene (24 Ma BP), the source rocks of the Bozhong Sag entered the mature stage. Hydrocarbon fluid was transported to the BZ19-6 buried hill along the fault and unconformity and accumulated into small-scale reservoirs. Due to the shallow depth (<2000 m) and the low diagenetic degree of the mudstone caprock of the Dongying Fm., oil and gas could escape easily (Figure 16(a)). At the end of the lower Minghuazhen deposition (5 Ma BP), the source rocks of the Bozhong Sag entered maturity and high-maturity stage, a large amount of hydrocarbon was generated and expelled, and a large amount of crude oil began to accumulate in the BZ19-6 structure. Some crude oil migrated upward along the fault and accumulated in the shallow Guantao Fm., forming the BZ19-4 reservoirs. This was the main accumulation period since 5 Ma (Figure 16(b)).
At the end of the upper Minghuazhen deposition (2 Ma BP), the source rocks of the Bozhong Sag entered a highly mature stage, with the generation of a large quantity of gas. Owing to the invasion by natural gas, combined with the mantle-derived inorganic CO 2 in the early reservoir, the gaseous hydrocarbon content increased and the critical temperature decreased significantly [8]. Under high-temperature and high-pressure conditions, the soluble components in the crude oil dissolved into natural gas and the heavy components precipitated to form asphalt. As a result, the early oil reservoir was transformed into a current condensate gas reservoir (Figure 16(c)).

Conclusions
(1) The BZ19-6 condensate gas reservoir is an ultra-high condensate gas reservoir. The average concentrations of C 1 , C 2-6 , and C 7+ in the hydrocarbon fluid were 66.8%, 14.66%, and 7.3%, respectively. The CO 2 content of the mantle-derived inorganic origin was as high as 11.24%. The average content of ground condensate was 615.96 g/m 3 , which indicates a gas reservoir with extremely high condensate content (2) Based on the analysis using different methods, the hydrocarbon fluid phase type of the BZ19-6 structure was revealed to be a condensate gas reservoir with oil rings. Its genesis resulted from the retrograde evaporation caused by burial depth, temperature increase since the Neogene, inorganic CO 2 charging from the deep mantle, and late natural gas invasion (3) The hydrocarbon accumulation process of the BZ19-6 condensate gas reservoir can be summarized as follows: a large amount of oil accumulated at the end of the lower Minghuazhen deposition (5 Ma BP), followed by a large amount of natural gas and mantle-derived inorganic CO 2 charged into the early formed oil reservoirs at the end of the upper Minghuazhen deposition (2 Ma BP). As a result, the content of gaseous hydrocarbons (C 1 +CO 2 ) increased, which led to large amounts of liquid hydrocarbons being dissolved in gaseous hydrocarbons, and significantly reduced the critical temperature of the hydrocarbon system. When the critical temperature is lower than the formation temperature and enters the critical condensate temperature range, secondary condensate gas reservoirs are formed

Data Availability
The data used to support the findings of this study are included within the article.

Conflicts of Interest
The authors declare that there are no conflicts of interest.