The gas-liquid-solid three-phase mixed flow is the most general in multiphase mixed transportation. It is significant to exactly solve the coupling hydraulic transient problems of this type of multiphase mixed flow in pipelines. Presently, the method of characteristics is widely used to solve classical hydraulic transient problems. However, when it is used to solve coupling hydraulic transient problems, excessive interpolation errors may be introduced into the results due to unavoidable multiwave interpolated calculations. To deal with the problem, a finite difference scheme based on the StegerWarming flux vector splitting is proposed. A flux vector splitting scheme is established for the coupling hydraulic transient model of gas-liquid-solid three-phase mixed flow in the pipelines. The flux subvectors are then discretized by the Lax-Wendroff central difference scheme and the Warming-Beam upwind difference scheme with second-order precision in both time and space. Under the Rankine-Hugoniot conditions and the corresponding boundary conditions, an effective solution to those points located at the boundaries is developed, which can avoid the problem beyond the calculation region directly induced by the second-order discrete technique. Numerical and experimental verifications indicate that the proposed scheme has several desirable advantages including high calculation precision, excellent shock wave capture capability without false numerical oscillation, low sensitivity to the Courant number, and good stability. The gas-liquid-solid three-phase mixed flow is the most general in multiphase mixed flow in pipelines. The present general formula of multiphase mixed flows in pipelines. Presently, the method of characteristics is widely used To solve classical hydraulic transient problems. However, when it is used to solve coupling hydraulic transient problems, excessive interpolation errors may be introduced into the results due to unavoidable multiwave interpolated calculations. To deal with the problem, a finite difference scheme based on the Steger Warming flux vector splitting is proposed. A flux vector splitting scheme is established for the coupling hydraulic transient model of gas-liquid-solid three-phase mixed flow in the pipelines. The flux subvectors are then discretized by the Lax-Wendroff central difference scheme and the Warming-Beam upwind difference scheme with second-order precision in both time and space. Under the Rankine-Hug oniot conditions and the corresponding boundary conditions, an effective solution to those points located at the boundary is developed by the second-order discrete technique. several desirable advantages including high calculation precision, excellent shock wave capture capability without false numerical oscillation, low sensitivity to the Courant number, and good stability.