Seismic 3D full-horizon tracking is a fundamental and crucial step in sequence analysis and reservoir modeling. Existing automatic full-horizon tracking approaches lack effective methods for representing the stratigraphic sequence relationships in seismic data. However, the inability to represent the stratigraphic sequence relationships fully and accurately makes it challenging to address discontinuous areas affected by faults and unconformities. To address this issue, we develop a knowledge graph representing the stratigraphic sequence relationship, which enables the simultaneous extraction of all horizon surfaces once the stratigraphic distribution of the seismic data is obtained. This method first generates horizon patches and calculates the fault attributes, followed by the construction of an initial knowledge graph that characterizes the overall distribution of horizon patches and faults. The initial knowledge graph comprises nodes and edges. The nodes represent horizon patches, and their attributes cover the geographical location information of the patches and faults. Simultaneously, the edges represent the relationship between the horizon patches, including the stratigraphic sequence relationship, and their attributes illustrate the potential for connecting these patches. Furthermore, we introduce a multilayer knowledge graph based on the point-set topology to fuse the nodes. This allows for the continuous merging of the horizon patches to obtain horizon surfaces across discontinuities with the constraints of fault attributes and stratigraphic sequence relationships in 3D space. Synthetic and field examples demonstrate that our approach can effectively represent stratigraphic sequence relationships and accurately track horizons located in discontinuous areas with faults and unconformities.