Natural Gas Migration Pathways and Their Influence on Gas Hydrate Enrichment in the Qiongdongnan Basin, South China Sea

2D and 3D seismic data and basin simulation were used to investigate the gas hydrate distribution and natural gas migration pathways in the Qiongdongnan Basin (QDNB). Hydrate-related amplitude anomalies and extensive bottom simulating re ﬂ ectors (BSRs) were mapped within the uppermost part. Based on their seismic re ﬂ ection characteristics, the three main types of natural gas migration pathways and their distributions in the QDNB were identi ﬁ ed through high-resolution seismic data. Basin modeling was carried out to document the migration e ﬃ ciency of di ﬀ erent migration pathways and their e ﬀ ects on hydrate enrichment. The basin modeling results show the following: (1) Diapirs, fault structures, and fractures constitute the three main types of natural gas migration pathways that transport the thermogenic gas from the deep to shallow layers in the QDNB. (2) The three migration pathways impact hydrate enrichment in di ﬀ erent ways. Diapirs and faults contribute signi ﬁ cantly to hydrate enrichment due to their higher migration e ﬃ ciency. In comparison, the migration e ﬃ ciency of the fracture systems is lower, with minimal bene ﬁ t to hydrate enrichment. (3) The natural gas hydrate in the QDNB is mainly distributed along the diapirs and deep faults and generally scattered around the fracture system. These conclusions indicate that the migration pathways in the QDNB are regionally distributed and are closely related to hydrate accumulation.


Introduction
Gas hydrates, also known as methane hydrates, are ice-like crystalline solids in which gas (mostly methane) is trapped within a water cage under high-pressure and lowtemperature conditions [1]. Since 1999, natural gas hydrate surveys have been conducted on the northern slope of the South China Sea (SCS) in China. Research has been carried out on the macroscopic geological conditions, such as the influence of tectonic activity, sedimentary environment, temperature and pressure conditions on the gas saturation, burial depth, and distribution of hydrate reservoirs [2][3][4][5]. In some studies, geological and geophysical data have been used to study the evolution of reservoir-forming processes [6,7].
The migration pathway is the key controlling factor for petroleum accumulation in conventional reservoirs [8,9], which can be identified through the reflection characteristics on the seismic profile. Migration is crucial for the formation and accumulation of natural gas hydrate, since most of the gas source of hydrates comes from deep thermogenic gas. Globally discovered hydrate enrichment areas have confirmed that the natural gas migration pathways are composed of faults and diapirs [10][11][12]. Mud diapirs and gas chimney structures developed in the Nankai trough and the Yuling basin of Korea, directly control the distribution of natural gas hydrate, whereas fault systems control hydrate formation in the Gulf of Mexico, USA [12][13][14][15]. In the SCS, some studies found that hydrate enrichment is related to diapiric structures, faults, and unconformities [4,10,16,17]. Furthermore, fluid migration pathways such as faults and fractures, submarine slumps, mud diapirs, and gas chimneys are relatively common in the north SCS [7].  proposed that the fluid migration pathways developed in the QDNB include gas chimneys, diapir structures, and faults (including synsedimentary faults) [18]. There are many kinds of migration pathways in QDNB, but the diapir structure and fault systems mainly controlled the fluid migration in the basin.
Different types of migration pathways have different effects on hydrate saturation and accumulation.  found that in the Shenhu area, SCS, the deep fluids flowed rapidly upward along a high-permeability fracture channel, forming a high-saturation natural gas hydrate in the gas hydrate stability zone (GHSZ) [19]. Other studies have shown that the migration pathways, GHSZ, gas source condition, and reservoir-forming characteristics are the theoretical basis for hydrate research [7,17,20]. These studies provide a good theoretical basis for gas migration and hydrate accumulation systems. In recent years, basin modeling has been widely used in hydrate research, and 2D and 3D models have been used successfully to reproduce hydrate distributions inferred from seismic data [21][22][23][24]. For example,  reconstructed the gas generation, migration, and accumulation process by restoring the burial history and thermal evolution history of the basin and analyzed the petroleum generation potential of organic matter in sediments and its contribution to gas hydrate accumulation. In addition, a 2D-restored structural model was used to simulate gas hydrate properties over time, focusing on the presence and extent of the GHSZ [25]. Kroeger et al. (2019) investigated the mechanisms for gas hydrate formation and resulting hydrate distribution using a 3D model of the southern Hikurangi Margin, New Zealand. Therefore, it is feasible to use basin modeling technique to simulate the process of deep fluid migration and gas hydrate formation.
The QDNB has been proved to have enormous gas hydrate resource potential, but it is still in the early stage of gas hydrate exploration. Previous studies on the gas migration pathway are only limited to the reflection characteristics description [17,20,[26][27][28]. Still, the influence of different migration pathways on the dynamic process of gas migration and hydrate enrichment has not been considered, which limits the understanding of the timing, mechanisms, and abundance of gas hydrate formation. In this study, the seismic interpretation and dynamic basin modeling have been integrated to study the seismic reflection characteristics and migration efficiency of different types of natural gas  2 Geofluids migration mechanisms and their influence on hydrate accumulation in the deep-water area of the QDNB, so as to provide a useful guideline for the evaluation of high-saturation hydrate accumulation area in the basin.

Geological Setting
The QDNB is located in the northwestern continental margin of the SCS, covering an approximate area of 3:8 × 10 4 k m 2 [19]. Several uplifts and depressions are distributed across the basin. The basin consists of four primary structural units, including the northern depression zone, northern uplift zone, central depression zone, and southern uplift zone, which can be further divided into more than 20 tertiary structural units ( Figure 1) [29]. The QDNB experienced rift (Eocene-Oligocene) and depression stages (after the early Miocene), with a typical "two-layer" structural pattern ( Figure 2) [30,32]. During the rifting period, with intense tectonic activities, it generated a lot of faults and mud diapirs. On a basinal scale, many magmatic diapirs occurred in the western central depression and the southern uplift in the deep-water area of the basin (Ledong sag and Lingshui sag, and Songnan Low Salient areas) [17,33,34]. In contrast, in the eastern part of the basin, a large number of faults were developed in the Baodao and Changchang sags, without diapirism. Basement faults are widely distributed in the basin, with generally three groups trending nearly EW, NE, and NW [35][36][37]. Moreover, the structural layer of depression represents the deposition in the thermal subsidence period. During the postrifting stage, the basin subsided with minor additional faulting [38].

Geofluids
Oligocene Yacheng and Lingshui Formations, the Miocene Sanya, Meishan, and Huangliu Formations, the Pliocene Yinggehai Formation, and the Quaternary Ledong Formation, bounded by the Sequence Boundary T80, T70, T60, T50, T40, T30, and T20 [26] (Figure 2). Previous studies revealed that the sedimentary environment in the Paleogene changed from the Eocene lacustrine to the Oligocene marine [26,40] (Figure 2). Since the Miocene, moderately deep to deep-water marine sediments have been deposited, and after the Pliocene, a clear shelf-slope system has been gradually differentiated [17,36]. Under these environmental conditions, three sets of source rocks are developed in the basin: (1) the Eocene lacustrine high-quality source rocks [41,42]; (2) the Oligocene coastal plain, offshore, and shallow marine facies source rocks. The marine and continental transitional and shallow marine source rocks in the Yacheng Formation are mainly gas-generating, and the marine mudstones in the Lingshui Formation are secondary source rocks [17,43,44]; (3) the Miocene moderately deep to deep-water marine source rocks are also secondary source rocks but are in the immature to low maturity stage with respect to petroleum generation [45,46].

Methods and Data
The dataset for this study includes 800 km 2 of recently acquired 3D seismic data and 5700 km of 2D seismic data within the QDNB. In the study, the software PetroMod™ v.2016 2D was used to simulate the study area's petroleum migration and hydrate formation process [47].
3.1. Model Construction. The seismic profiles through well A, well B, and well C are used to define the architecture of the model (Figure 1), which, respectively, represent different migration pathways: gas chimney (line 1), fault structure (line 2), and faults and fracture system (line 3) (the positions of the seismic profiles are shown in Figure 1). The stratigraphic framework is based on the interpretation of CNOOC for petroleum exploration and gas hydrate exploration for shallower sediments [48,49]. Then, the faults, gas chimneys, diapirs, and mass-flow deposits were delineated in the stratigraphic framework to complete the construction of the basin model. The major gas chimneys and Meishan formations are represented by high permeable sediments, and faults are considered open faults in the model. The T2 and T3 models were chosen as the petroleum generation kinetic models corresponding to kerogen types II and III, and the petroleum saturation method was selected as the expulsion model [48].

Geological Framework and Lithologic
Parameters. The lithology model for shallow sediments (<300 m) is defined by the core descriptions and gamma-ray (GR) logs of well A, well B, and well C. The lithologies for deeper layers are based on the research results on regional sedimentary facies. In the process of lithology assignment of a given model layer, several different lithologies were mixed according to the statistical composition. Lines 1, 2, and 3 show some differences in lithologic assemblages. The depositional environment and sedimentary facies changed progressively in response to the paleotopography. Different depositional facies must be identified and assigned with other parameters. The proportion of sandstone for each formation in line 1 is listed in Table 1, derived from the research results of Kong (2010).

Erosion.
The erosion parameters of the study mainly refer to the research data of Kong (2010). Referring to the previous estimation of erosion in a single well [50], combined with the geological reality of the QDNB, the thickness eroded during the unconformities T60, T40, and T70 was calculated to be 200-300 m, 100-300 m, and 100-200 m, respectively. The specific amount depends on the location and characteristics of each section.

Paleo-Heat
Flow, Paleo-Temperature, and Paleo-Water Depth Parameters. Paleo-heat flows, paleo-temperatures, and paleo-water depths in this study were taken from the   Table 2). The basin modeling software was used to calculate paleo-temperatures, and the simulated present-day temperature gradient was adjusted using the measured temperature gradient. The water depth also varied significantly in different areas and sedimentary periods. The paleo-water depths were extracted according to the sedimentary facies and sea-level change curves during different periods (Table 2) [50].

Organic Geochemical Parameters of Source Rocks.
The organic geochemical parameters of source rocks mainly refer to the research data of Kong (2010). The main types of kerogen in the study are III and II [48], and the geochemical parameters of three sets (the Miocene, the Oligocene, and the Eocene) of source rocks are determined by reference to Kong (2010) ( Table 3).

Types and Characteristics of the Natural Gas Migration
Pathways. Based on the detailed description and interpretation of seismic profiles in the key target areas, drilling and logging results of the QDNB, three main gas migration pathways, including magmatic diapirs and gas chimneys, fault structures, and fracture systems were identified in the basin. Gas chimneys act as conduits for fluid transport and storage [53][54][55], which shows abnormal reflections on the seismic profile due to the vertical migration of natural gas or other fluids. As shown in Figure 3, there is a large-scale gas chimney structure in profile line 1, and the lower part of the gas chimney is a magmatic diapir (developed in Songnan Low Uplift). The gas chimney in this study area showed a pull-down, narrow, vertical column with fuzzy reflection or chaotic reflection interior, and the reflection events at the flanks of a gas chimney are also interrupted or terminated across the border of a gas chimney (Figure 3). What is more, gas chimneys often end at the BSRs. The occurrence of the BSRs with high amplitude and poor continuity is the most important and effective indicator of layered gas hydrate [5,[56][57][58][59]. As we can see in Figure 3, BSRs with strong amplitude, wide distribution, and reversed polarity with seabed are developed at the top of the gas chimney ( Figure 3). In addition, the top of the gas chimney is usually blocked by the caprock, and the free gas widely accumulated under the caprock and around the gas chimney, as well as the deposition of shallow mass transport deposits (MTDs), leading to the lateral extension of the strong reflectors and the bright spots at the periphery of the chimney [17,[59][60][61]. In the seismic profile of line 1, the periphery of the chimney is surrounded by bright spots and strong amplitude reflection, and above the BSRs at the top of the gas chimney is the GHSZ (Figure 3).
Faulting is widely developed in the QDNB, and there are mainly interpreted based on high-resolution seismic profiles by recognizing abrupt changes in the continuity of seismic reflectors [61]. The seismic profile shows that the deep faults cut through the thick formations, connecting the deep and shallow stratum, and in the upper of the large faults, some BSRs and bright spots were observed ( Figure 4). As shown in Figure 4, the BSRs in the upper faults showed a relatively high reflection magnitude and continuity (Figure 4), which often indicates the presence of hydrates. Due to the extensive accumulation of free gas and the deposition of MTDs, many bright spots with strong reflections can be observed at the top of faults. Besides deep faults, natural gas also accumulates and diffuses at the end of some shallow faults, forming bright spot reflections in seismic profiles ( Figure 4). In sedimentary strata, fractures are produced by strata deformation, forming complex networks in the Earth's crust and serving as critical pathways for fluid migration in subsurface geological media [62][63][64]. The fracture system is hard to identify on seismic profile due to its small scale and chaotic characteristics. As is shown in profile line 3, some gas fuzzy zone related to the fracture systems can be observed above the T40 interface, with bright spots and BSRs distributed at the top zones ( Figure 5).

Numerical Modeling Results.
In the numerical simulation model, the density of the arrow was used to represent the amount of oil and gas migration to compare the migration efficiency of different natural gas migration pathways. The numerical modeling results show that the gas was generated from the deep coal-measure and lacustrine source rock from Yacheng and Lingshui formation, migrated along with the diapir/sediment interface and then vertically upwards, forming a gas chimney to the shallow sediments above (Figure 6(a)). The gas chimney migration pathway is characterized by continuous vertical gas distribution, sufficient gas supply, and high migration efficiency. Furthermore, large amounts of natural gas migrated upward along large faults, then accumulated and diffused at the upper end of the faults in profile line 1 (Figure 6(a)) and 2.5-7.5 km distance in profile line 2 ( Figure 6(b)). But there are few faults in the lower part of well B, and only a small amount of oil and gas are migrated to the shallow layer ( Figure 6(b)). The result of the simulation shows that the characteristics of deep fault are high efficiency, long distance, and relatively large-volume petroleum migration and     Figure 1). 8 Geofluids that present the characteristics of "wide but sparse" and are discontinuous in the vertical direction ( Figure 6(c)). Compared with gas chimneys and deep faults, fracture systems are likely to have lower migration efficiency.  (Figure 7(a)). In the lower gas  9 Geofluids hydrate layer, no effective AC logs were obtained, but the RT curve shows obvious enhancement (Figure 7(a)), indicating that there are a lot of hydrates in this layer. However, no hydrate layers were drilled in well B and well C, and the raw log data of the RT and AC curves have no distinct indications of their trends. The AC curves showed a normal increasing trend with the depth, while the change of the resistivity curve was not significant with the mean value of 1.5 Ω·m and 1.2 Ω·m, respectively (Figures 7(b) and 7(c)).

Distribution of the Natural Gas Migration Pathways.
The seismic interpretation results show that the distribution of different migration pathways in the basin has certain regional variations. The diapir structures and gas chimneys were well developed in the western deep-water area of the basin, while in the eastern part of the basin, especially in Baodao and Changchang sags, faults are widely developed (Figure 8).
Gas chimneys are suggested to act as a major pathway for fluid migration from the intermediate gas reservoir to shallow sediment above the BSRs [65]. From the magmatic intrusions identified in seismic profiles, there are differences in scale and intrusion period of these diapirs across the area. The intensity of magmatic diapirism is greater in the west of the basin than in the east. A large volume of magma was intruded along deep faults, leading to the formation of low uplift and a series of secondary fault steps in the central depression zone (Figure 9). Magmatic diapirs and gas chimneys are the most widespread type of fluid migration structures in the low salient areas in the southwest of the Central depression [66]. Magmatic diapirs and gas chimneys are relatively concentrated in the deep-water area. Most of them distribute in the Ledong and Lingshui sags and Songnan and Lingnan low salient zones [11,17,67]. It is seen from the seismic profile that there is more faulting in the east than in the west, especially below surface boundary T60. Since the Neogene, the faults usually terminated at the T40 interface and connected the deep strata below T60, with high fault density and over a large area ( Figure 10). In addition, there are a series of normal faults developed in the Baodao and Changchang depressions and the Lingnan Low Uplift area, and the southern uplift area also developed some normal faults in the Ledong Formation [37,68,69]. It is believed that the Neogene deep normal faults are the main migration pathways in the Baodao and Changchang Sags in the eastern area of the basin. Fractures, however, are commonly developed in the whole basin and banding distribution or associated with diapirs and faults [20,70].

The Source of Natural Gas Forming
Hydrates. The thermogenic gas generated from the deep source rocks and the microbiological gas generated from the shallow organic matter can provide natural gas for hydrate accumulation. Previous studies revealed that the natural gas composition for hydrate was biogenic gas, thermogenic gas, and mixed biogenic-thermogenic gas [11,28,41,71,72]. Methane and C 2+ hydrocarbon were dominant in the hydrate-related gas component. According to the C 1 /ðC 2 + C 3 Þ value to distinguish the origin of natural gas, the biogenic gas is methane-dominated with a very high dry coefficient, and C 1 /ðC 2 + C 3 Þ > 1000, while the thermogenic gas usually has a large amount of C 2+ hydrocarbons and C 1 /ðC 2 + C 3 Þ < 100 [73][74][75]. In addition, it is accepted that the values of the CH 4 less than −55‰ indicate microbial origin, while the values more than −50‰ generally represent their thermogenic origin [11,28,73]. Accordingly, in this study area's well A, hydrate-decomposed gas has C 1 /ðC 2 + C 3 Þ of 20 to 200, the isotope of methane of -70‰ to -50‰, and δ 13 C 2 of -32‰ to -24‰, indicating a mixed biogenicthermogenic gas source and mainly coal-type gas  10 Geofluids ( Figure 11). Previous studies showed that the BD 13-1 field, the BD 19-2 field, and the YC 13-1 field all have a thermogenic gas origin, sourcing from deeply buried coaly type gas pools and they preliminarily concluded that the C 2+ hydrocarbon gases (thermogenic gas) within the hydraterelated gas from site GMGS5-W08 predominantly originated from deeply buried coaly type gas pools (Figure 11), while the methane has a mixed thermogenic and microbial origin [28]. In addition, Xu et al. (2021) analyzed some wells (Well-b and Well-c) in the study area and found that the hydrate decomposition gas samples were mainly within the thermogenic gas and mixed thermogenic and biogenic gas range [18] (Figure 11). Therefore, well A's gas source analysis results in the study area are completely consistent with previous research. Furthermore, researches show that in the QDNB, the biogenic gas in the hydrates was derived from the Marine source rocks deposited in the Huangliu formation and their overlying strata, and the thermogenic gas may have come from the coaly type gas pools or source rocks developed in the deep Yacheng and Lingtou formations in the Paleogene [11,30,41,75,76]. Therefore, hydraterelated gas in the shallow strata has a close relationship with the deep petroleum reservoirs in the QDNB [11].

The Efficiency of Natural Gas
Migration. The migration pathways identified here, namely gas chimneys above diapirs, faults, and fracture systems, have different migration efficiencies in the basin [4].

Gas Chimney.
Gas chimneys are widely developed in the basinal area, suggesting a viable migration pathway for gas-bearing fluids. In the seismic profile, there are many bright spots on the top of the gas chimney with strong amplitude and continuous BSR distribution with a wide range ( Figure 3). Besides, according to the numerical simulation results, many gases migrated from the deep to the shallow layer along with the gas chimney ( Figure 6(a)). The GHSZ in the basin models correlates well with the BSR locations in the seismic profile of line 1 (Figure 3). Furthermore, three sets of gas hydrate-bearing sections were identified in

Geofluids
well A with approximately 40 m thickness. All of this shows that chimneys effectively connect the deep natural gas source to shallow traps or submarine GHSZ [17,77,78]. The gas migrates upward through the gas chimney to the GHSZ, where it forms gas hydrates. The important natural gas pathway motivated the migration of deep thermogenic gas upward to the shallow layers. However, the distribution of diapirs in the Neogene is limited with their intrusion to shallow layers, and their effects on large-scale upward migration of oil and gas are limited, so small faults and fractures developed in the upper of the diapirs are used to migrate the gas further into the shallows.

Fault
Structure. The formation and evolution of faults play important roles in the petroleum migration process [79,80], as they can connect source rock and trap as "a bridge and link" [58,81,82]. Since the Baodao-Changchang area with relatively developed faults has not been drilled, we selected the profile with relatively developed faults in Lingnan low salient to research (the profile of line 2). In the seismic profile, there are many bright spots on the top and around the faults, especially in the area where there are many faults in the NW of the profile, and the distribution of bright spots and BSRs is larger (Figure 4). The gas migration channels shown on profile line 2 correspond with the bright spots in the seismic profile, suggesting that the fault structures are also an important natural gas migration path-way, and the deep faults control the distribution of bright spots in shallow layers. In addition, the simulation results show that oil and gas can migrate efficiently and over long distances through faults ( Figure 6(b)). The drilling results of Well B showed that no hydrate was encountered, which is presumed to be due to fewer faults development at the drilling location, making it more difficult for gas to migrate to the GHSZ. The long-term activity of deep faults cut down to deep, mature source rocks facilitates oil and gas movement to shallow reservoirs over geological time. In some cases, oil and gas can migrate to the surface (e.g., the seafloor). Therefore, the combination of deep and shallow faults is essential for forming integrated migration channels to the surface or near-surface. The fault structures are characterized by high efficient, long-distance transportation of fluids.

Fracture System. Fractures and microcracks have been
proposed as effective pathways for natural gas and oil in the QDNB [20,83,84]. In addition, overpressure occurred in the basin [31,85]. When overpressure exists in the lower strata, numerous microfractures have the potential to form in the upper strata. So, the episodic opening activities of the fractures would not only constitute the main channel for energy release in an overpressured system but also form vertical fluid migration pathways [31,86]. On the one hand, a large number of primary joints and contraction joints were developed during the process of magma condensation of intrusive   Figure 1).  Figure 11: Cross-plot of δ 13 C 1 and δ 13 C 2 values for the gas samples collected from the QDNB (the δ 13 C values of natural gas from BD13-1, BD19-2 gas fields are cited from [28], the δ 13 C values of natural gas from YC13-1 is cited from [41], the δ 13 C values of natural gas from wellb, well-c are cited from [18]). 13 Geofluids during the postrift thermal subsidence period in the study area [88,90]. In addition, studies indicate that the intense tectonic activity controls the nature, occurrence, and scale of fractures that develop in the buried hill [59,91]. Due to the tectonic deformation, a fracture network with highangle shear fractures and tensile shear fractures almost parallel to the edge of the Magmatic diapir was developed near the contact surface between the diapir and surrounding rocks [92,93].
In this study, the numerical simulation results show great correspondence with the bright spots on the seismic profile, indicating that the fracture system played a role in upward gas migration from depth. Gas hydrate samples have not been found in well C, indicating that fractures poorly affect oil and gas migration. Many basement faults terminate below the T60 surface, which can migrate the gas and oil to intermediate and shallow layers. So, the fracture system above surface T40 is usually regarded as a relay pathway to continue upward migration to the GHSZe and formation of the gas hydrates. The fracture system has a short transportation distance and is usually combined with other migration pathways like fault structures and magmatic diapirs to be regionally effective.

The Influence of Different Migration Pathways on
Hydrate Enrichment. We suggest here that the distribution of natural gas hydrate is closely related to deep structures, especially faulting and diapirs, along with fracture systems and gas chimneys at shallower depths. Deep faults correlate with the distribution of shallow bright spots and hydrate accumulation areas, but the influence of different migration pathways on hydrate enrichment differs [4]. Based on this, three hydrate enrichment models controlled by gas  Figure 12: Hydrate enrichment models controlled by different migration pathways: (a) gas chimney; (b) fault structure; and (c) fracture system. 14 Geofluids chimneys, fault structures, and fracture systems are proposed in this paper ( Figure 12). It has been documented that abundant gas has been generated in the QDNB, with both microbial and thermogenic origins [57,92]. The presence of hydrate in the QDNB is closely related to the diapiric structures, deep faults, and fracture systems. Gas chimneys often connect the deep and shallow strata, providing pathways for the deep thermogenic gas to migrate to the shallow GHSZ, and correlate with highabundance and high-saturation hydrate reservoirs [4,94]. Some large faults extend from depth to the shallow strata, producing favorable vertical migration channels and high migrating ability within the basin. Therefore, the deepsourced thermogenic gas migrated vertically through the faults to the shallow area, forming high-abundance hydrate resources [78,94]. Furthermore, fractures make up the favorable conduits and channels for continuous gas supply to the shallow section [95]. Although the migration efficiency of the fracture systems is lower than gas chimneys and fault structures, it usually develops along with the shallow small faults, leading to a relatively short distance of oil and gas migration or merely in situ accumulation. Therefore, the distribution of gas hydrate controlled by the fracture systems is scattered, and the large-scale, contiguous hydrate deposits are hardly formed by this mechanism alone.
Currently, drilling is mainly concentrated in the areas where diapirs and gas chimneys developed, and the results show that magmatic diapirs and gas chimneys control the distribution of hydrates. Baodao and Changchang sags are where fault development has not been drilled. However, the study shows that the faults in this area also impact hydrate enrichment; therefore, this area may be a favorable area for future hydrate exploration.

Conclusions
This study proposes three main gas migration pathways in the QDNB, SCS, namely gas chimneys, fault structures, and fracture systems. The gas chimneys appear as a vertical feature on seismic, often thicker at the base and narrowing towards the top, with fuzzy bands or wiped-out reflectors. Bright spots are generally observed on top of the gas chimneys. Fault structures show obvious stratigraphic breaks on the seismic profile, with blank reflections and bright points seen at the end of, and adjacent to, the fault. The fracture system is in the form of dendrites and flowers on the seismic profile, it is characterized by steep angles, small fault distance, and blank reflections, and it is not easy to identify.
The diverse influences of the three migration pathways on natural gas migration and hydrate enrichment are as follows: (1) The gas chimney migration pathway is characterized by continuous vertical gas distribution, sufficient gas supply, and high migration efficiency, and it makes a significant contribution to hydrate enrichment. (2) The characteristics of deep faults are high efficiency and long distance, and it also has a significant influence on hydrate enrichment. (3) Compared with the other two migration pathways mentioned above, the fracture systems have low migration efficiency and make less contribution to hydrate formation, and hence, it is not conducive to the formation of widely distributed and high-saturation gas hydrate.
The natural gas hydrates in the QDNB are mainly correlated with diapir structures, gas chimneys, and deep faults and scattered around the fracture systems. The gas chimney is the main gas hydrates migration pathway in the Ledong -Lingshui sags and Lingnan -Songnan low salient in the central depression zone, while in the Baodao and Changchang sags, the deep faults play an important role in natural gas migrating and hydrates enrichment.

Data Availability
The basin numerical simulation data used to support the findings of this study are included in the article. Previously researched data were used to support this study. These previous studies (and datasets) are cited at relevant places within the text as references [47,[50][51][52].

Conflicts of Interest
The authors declare that there is no conflict of interest regarding the publication of this paper.