Three-dimensional geological modeling is essential for simulating natural gas hydrate (NGH) productivity and formulating development strategies. Current approaches primarily concentrate on the single-phase modeling of either hydrate or free gas layers. However, an increasing number of instances suggest that the three-phase coexistence zone, which includes hydrate, gas, and water, is common and has become a focal point of international research, as this type of reservoir may present the most viable opportunities for exploitation. At present, there exists a significant gap in the research regarding modeling techniques for such reservoirs. This study undertakes a comprehensive modeling investigation of the three-phase zone reservoir situated in the sand layer of the Qiongdongnan Basin. By employing deterministic complex geological modeling techniques and integrating existing seismic and logging data, we have developed a three-phase coexistence zone model that precisely characterizes the interactions between geological structures and utilizes them as auxiliary constraints. This approach effectively mitigates the potential impact of complex geological conditions on model accuracy. Through a comprehensive analysis of 105 seismic profiles, we enhanced the model’s accuracy, resulting in the creation of a three-phase coexistence zone model comprising 350,000 grids. A comparison between the modeling results and well data indicates a relatively small error margin, offering valuable insights for future development efforts. Furthermore, this method serves as a reference for modeling hydrates in marine environments characterized by three-phase coexistence on a global scale.

Based on comprehensive interpretation of three-dimensional seismic data and quantitative analysis of basin-boundary fault activity in the Nanpu Sag, this study employs subsidence history backstripping and equilibrium profile techniques to reconstruct the structural evolution of the main profile. The results indicate that the Cenozoic evolution of the Nanpu Sag can be divided into a syn-rift stage and a post-rift stage, with the syn-rift stage further subdivided into Rift I and Rift II episodes. During Rift I, tectonic activity was primarily controlled by the NE- and NEE-trending Xinanzhuang Fault, Shabei Fault, and No. 2 Fault Zone, which formed under a NW–SE extensional stress regime and governed the development of NE- or NEE-trending faults and associated sedimentary subsidence centers. In Rift II, tectonic activity was dominated by a southward-curved normal fault system, composed of the Xinanzhuang, Gaoliu, and Baigezhuang faults, as well as the Shabei Fault, reflecting the influence of a near N–S ex-tensional stress field. The progressive southward migration of the Sag’s subsidence center over time—from the Linque sub-sag in the third section of the Shahe Formation to the Caofeidian sub-sag in the Dongying Formation—and noting, coupled with the pronounced left-lateral strike-slip characteristics of the Baigezhuang Fault and No. 4 Fault, and regional tectonic evolution analysis of the Bohai Bay Basin, support the proposal that a strike-slip extension mechanism—characterized by lateral strike-slip and forward extension—constitutes the fundamental developmental model of the Nanpu Sag. This study deepens the understanding of the tectonic evolution of the Nanpu Sag and provides new insights in-to the dynamic mechanisms governing the formation of similar Sags in the Bohai Bay Basin.