It is assumed that a pipeline is laid through a vertical fault fracture zone,and is excited by seismic ground motion modelled as stationary stochastic process.For horizontal incidence of waves,the cross-PSD (Power Spectral Density) function is developed using wave propagation theory,while for vertical incidence of waves the cross-PSD function is composed by auto-PSD model,coherence model and site response model.As the seismic input,the cross-PSD function is used to calculate the the axial and lateral seismic responses of underground pipeline through the fracture zone.The results show that the incident directions of seismic waves,width and soil property of the fracture zone have great influence on underground pipeline.It is suggested that the flexible joints with appropriate stiffness should be added into the pipeline near the interfaces between the fracture zone and the surrounded media. It is assumed that a pipeline is laid through a vertical fault fracture zone, and is excited by a seismic ground motion modelled as a stationary stochastic process. For horizontal incidence of waves, the cross-PSD (Power Spectral Density) function is developed using wave propagation theory while for vertical incidence of waves the cross-PSD function is composed by auto-PSD model, coherence model and site response model. As the seismic input, the cross-PSD function is used to calculate the the axial and lateral seismic responses of underground pipeline through the fracture zone. The results show that the incident directions of seismic waves, width and soil property of the fracture zone have great influence on underground pipeline .It is suggested that the flexible joints with appropriate stiffness should be added into the pipeline near the interfaces between the fracture zone and the surrounding media.