The current study examines the methane gas hydrate stability zone (GHSZ) in the Shenhu area in the northern South China Sea (SCS) as an example to calculate the thickness of the GHSZ and reconstruct its evolution since 8.2 Ma. Two mechanisms for typical double BSRs in the Shenhu area are shown, and the relationship between the evolving thickness of the GHSZ and the dynamic accumulation of NGHs at typical stations in the Shenhu area is clarified. The results show that the thickness of the GHSZ varies over time with overall thickening in the Shenhu area. The current thickness of the GHSZ is between 160.98 and 267.94 m. Two mechanisms of double BSRs in the Shenhu area are summarized: the double BSRs pattern based on changes in formation temperature, pressure and other conditions and the double BSRs pattern based on differences in gas source and composition. The formation process and occurrence characteristics of double BSRs and hydrate at site SH-W07-2016 in the Shenhu area are also closely related to the changes in thickness of the GHSZ. In addition, the age when gas source first enters the GHSZ has a considerable influence on the dynamic accumulation process of hydrate. Since the formation of hydrate above the BSR at site SH-W07-2016, the GHSZ has experienced up to two periods of thickening and two periods of thinning at this site. With the changes in the thickness of the GHSZ, up to two stages of hydrate formation and at most two stages of hydrate decomposition have occurred. This paper is of great value for understanding the formation of multiple bottom-simulating reflectors (BSRs) as well as the migration, accumulation and dissipation of natural gas hydrate (NGH) during the dynamic accumulation process.

It is of great significance to study the spatial-temporal variations of the thickness of the gas hydrate stability zone (GHSZ) to understand the decomposition, migration, accumulation and dissipation of gas hydrate, the corresponding relationship between bottom-simulating reflectors (BSRs) and gas hydrate, and the distribution of heterogeneous gas hydrate. We selected the Dongsha region in the South China Sea (SCS) as the research object to calculate the spatial-temporal variation of the GHSZ since 10 Ma, analyzed the main factors affecting the thickness of the GHSZ, discussed the dynamic accumulation processes of gas hydrate, and proposed an accumulation model of gas hydrate in the Dongsha region. The results show that the thicknesses of the GHSZ in the study area were between 0 m and 100 m from 10 to 5.11 Ma, and the relatively higher bottom water temperature (BWT) was the key factor leading to the thinner thickness of the GHSZ during this period. From 5.11-0 Ma, the thickness of the GHSZ gradually increased but showed several fluctuations in thickness due to changes in the geothermal gradient, seawater depth, BWT, and other factors. The decrease in the BWT was the main factor leading to GHSZ thickening from 5.11 to 0 Ma. The thicknesses of the GHSZ are between 110 m and 415 m at present. The present spatial distribution features show the following characteristics. The GHSZ in the deep canyon area is relatively thick, with thicknesses generally between 225 m and 415 m, while the GHSZ in other areas is relatively thin, with thicknesses between 110 m and 225 m. Based on the characteristics of the GHSZ, two hydrate accumulation models are proposed: a double-BSRs model due to thinning of the GHSZ and a multilayer hydrate model due to thickness changes of the GHSZ, with single or multiple BSRs.