Simulation method for the management of a pipeline network comprising nodes (S 1, S 2, C 1), said simulation method being implemented using a processing unit and comprising: - defining one or more operating states ST p of the pipeline network characterizing operating conditions of said pipeline network, where p is an integer comprised between 1 and P, P being the number of operating states, characterized by the following steps: - determining, for each operating state ST p, a three-dimensional dynamic matrix DM STp whose each coefficient DM STp ( i,j,t k) corresponds to a pressure variation value from an initial pressure value at a j-th node at a k-th time step t k, following a variation of a predetermined reference value (ΔFR), at a i-th node, of a flow rate value when the operating conditions of the pipeline network correspond to said operating state ST p, where i and j are integers comprised between 1 and N, N being the number of nodes, and k is an integer comprised between 1 and T, T being a number of time steps corresponding to a predetermined time horizon, - estimating, for a given node at a given moment, a pressure value on the basis of an operating schedule of the pipeline network, said operating schedule providing information regarding variations of the flow rate value for each node and evolution of the operating conditions of said pipeline network until said given moment, said estimation using the one or more operating states and the associated three-dimensional dynamic matrices and - comparing the estimated pressure value at the given node at the given moment with a predetermined pressure threshold (P lim) and, if the estimated pressure value is greater than or equal to said predetermined pressure threshold, modifying the operating schedule of the pipeline network until the estimated pressure value is lower than the predetermined pressure threshold. - defining one or more operating states ST p of the pipeline network characterizing operating conditions of said pipeline network, where p is an integer comprised between 1 and P, P being the number of operating states, characterized by the following steps: - determining, for each operating state ST p, a three-dimensional dynamic matrix DM STp whose each coefficient DM STp ( i,j,t k) corresponds to a pressure variation value from an initial pressure value at a j-th node at a k-th time step t k, following a variation of a predetermined reference value (ΔFR), at a i-th node, of a flow rate value when the operating conditions of the pipeline network correspond to said operating state ST p, where i and j are integers comprised between 1 and N, N being the number of nodes, and k is an integer comprised between 1 and T, T being a number of time steps corresponding to a predetermined time horizon, - estimating, for a given node at a given moment, a pressure value on the basis of an operating schedule of the pipeline network, said operating schedule providing information regarding variations of the flow rate value for each node and evolution of the operating conditions of said pipeline network until said given moment, said estimation using the one or more operating states and the associated three-dimensional dynamic matrices and - comparing the estimated pressure value at the given node at the given moment with a predetermined pressure threshold (P lim) and, if the estimated pressure value is greater than or equal to said predetermined pressure threshold, modifying the operating schedule of the pipeline network until the estimated pressure value is lower than the predetermined pressure threshold.