## 3 Analytical Methods  
It is challenging to identify all faults based upon well information ascribed to the large spacing and insignificant density of wells in the research area. As the exploration technology advances, the technique of interpreting the location and scale of fault development through three-dimensional seismic data is becoming mature step by step. Considering the condition of low well density in the explored area, the identification of faults is grounded in 3D seismic data, and calibrated by well data, and the faults are identified, in line with the reflection's continuity. As suggested by our research findings, the faults in the seismic data often display the following characteristics: (1) the seismic wave reflection axis shows a misalignment; (2) the number of reflection axes suddenly increases, decreases, or disappears, and the interval of wave groups suddenly changes; (3) the shape of the reflection wave co-axis changes, and the reflection is messy or there is a blank zone; and (4) the standard reflection wave co-axis is bifurcated, merged, distorted, and with phenomena such as the strong phase transition, and the appearance of anomalous waves [51]. Grounded in the characteristics of the fault in the seismic section, this paper identifies faults in multiple directional sections. To be specific, the seismic attribute (ant attribute) is employed to analyze and identify trivial faults. Last but not least, the fault closure is completed and the precise interpretation of the faults in the investigated area is realized. The Zhuhai Formation of the Baiyun Sag was deposited during the transition period, and its tectonic evolution is featured by the superposition of multiple rotations and the composite of multiple genesis mechanisms. Faults of dissimilar scales, diverse orientation [52], and distinct natures are criscrossed and stacked up and down in the depression, which makes the Zhuhai Formation form a relatively complex fault control sand coupling system [53]. In an effort to clarify the type of sandstone distribution controlled by faults in the explored area, this paper summarizes the combination styles of syn-sedimentary faults in the research area on the basis of the interpretation of the faults, based upon the formation status of the faults [54], the fault-cutting-stacking relationship, the influence of the faults on the palaeomorphology, etc., and establishes three kinds of fault combinations that have distinct influences on the spatial distribution of sand bodies [55; 56] and have dissimilar impacts on the spatial distribution of sands. The extraction of seismic attributes, followed by their conversion into geological information [57], and, ultimately, the prediction of sandstone distribution guided by sedimentological principles, constitutes the paramount approach in seismic sedimentological analysis (i.e., seismic attribute analysis) [44; 58; 59]. A geological interpretation of seismic properties is the key to seismic properties analysis. It is the most extensively employed method in seismic sedimentology to throw light upon the distribution of sand bodies by analyzing the plane characteristics of seismic properties [60]. Ascribable to the deficiency of drilling in the investigated area, we choose to employ the measured data of the sand ratio to select the seismic attributes. Nevertheless, seismic attributes can be chosen by utilizing the sandstone thickness as a percentage of the recordeddata. Moreover, grounded in the link between seismic attributes and sandstone thickness, seismic attributes can be favored based on the link between seismic attributes and sandstone thickness as a percentage at statistical well locations. The intersection analysis of well-point seismic attributes and sandstone thickness as a percentage was performed within the statistical unit of the explored region. Afterwards, the correlation coefficients between the seismic attributes and sandstone thickness as a percentage of each well were counted by implementing this approach.