本发明涉及水利工程技术领域，尤其涉及一种压力管道输水系统非恒定流模型的优化方法。
背景技术：
长距离输水工程的特点为输水管线长、地形起伏多变和水流条件复杂等，其过渡过程模拟计算的影响因素比较多，除了输水系统自身影响因素外，自身影响因素如摩阻、压力、温度等，计算模型因素包括时间步长选取、采用的计算方法等也至关重要。长距离输水管道瞬变流计算时，时间步长的取值大小是影响水力过渡过程计算结果的重要因素之一，时间步长取值太大，无法满足计算精度要求，方程的解难以达到稳定；取值太小，若输水管线长，导致计算工作量较大。因此合理确定计算管道分段数，既保证计算的精度，又减少计算的工作量，是长距离管道输水泵站工程的水力过渡过程数值模拟必须考虑的一个重要因素。
现有通常调整波速法计算时间步长，调整波速法是通过在一定范围内对波速进行调整，使压力管道输水系统内各管段的分段数是整数且时间步长相等。该法是目前最为常用的压力管道(复杂管系)分段方法。
调整波速法存在的问题如下：
(1)由于调整波速法是依据压力管道输水系统内波传播时间最短的一条管段(通常为管长最短的管段)作为计算基础，在计算时间步长时，不对短管的波速并没有调整，使得压力管道输水系统的管道分段受到限制，有可能使时间步长Δt的取值过小。
(2)波速调整后，水锤波从压力管道输水系统起点传播到终点的总时间T′与未调整波速时的总时间T很难一致，存在一个差值ΔT＝T′-T引起误差。
(3)现有调整波速法没有体现优化思想,分段时只要满足公式a就算满足要求。
尽管调整波速分段法的提出主要是基于水锤波速影响因素多，不能精确得知时间步长，得到的时间步长与真实值间存在误差。
技术实现要素：
本发明的目的在于提供一种压力管道输水系统非恒定流模型的优化方法，从而解决现有技术中存在的前述问题。
为了实现上述目的，本发明所述压力管道输水系统非恒定流模型的优化方法，所述方法包括：
首先，获取压力管道输水系统中的有压管段，对水锤波传播时间最短有压管段的水锤波速进行调整，获取最大时间步长Δtmax，计算获取最大时间步长Δtmax所对应的最短有压管段的最小分段数N0min；
然后，通过最小分段数N0min计算其他任意一个有压管段i的最小波速调整幅度fimin及最小波速调整幅度fimin所对应的最优分段数Ni最优；
最后，通过最小波速调整幅度fimin，调整水锤波速，保证调整前、后水锤波从压力管道输水系统起点传播到压力管道输水系统终点的总时间相同。
优选地，对水锤波传播时间最短有压管段的水锤波速进行调整，获取最大时间步长Δtmax，计算获取最大时间步长Δtmax所对应的最短有压管段最小分段数，具体按照下述步骤实现：
计算压力管道输水系统内水锤波传播时间最短有压管段的管长L0、波速a0和分段数N0，其中，管长L0、波速a0存在关系式(1)：
压力管道输水系统中除水锤波传播时间最短有压管段外任意一个有压管段分段数Ni满足公式(2)，所述Ni为正整数:
其中，i表示压力管道输水系统中除水锤波传播时间最短有压管段外任意一个有压管段的编号，i为大于等于1的正整数，m表示压力管道输水系统中除水锤波传播时间最短有压管段外有压管段的总数量，α表示压力管道输水系统中除水锤波传播时间最短有压管段外其他有压管段的水锤波波速；fi表示编号为i的有压管段的波速调整系数；Li表示编号为i的有压管段的长度；△t表示水锤波传播时间最短管段的关系时间步长；
所述△t的取值范围符合公式(3)：
fmax表示|fi|的最大取值，在计算水锤波传播时间最短有压管段的波速调整的关系时间步长时，|fi|＝0.15；T表示未调整波速时，水锤波从压力管道输水系统起点传播到终点的实际总时间；
由于水锤波传播时间最短有压管段的波速调整不会影响对压力管道输水系统的计算结果，故，在获取△t的取值范围时，设|fi|＝0.15，则fmax＝-0.15或fmax＝0.15，因此，△t的取值范围符合公式(4)：
通过公式(4)，得到最大关系时间步长Δtmax的计算公式(5)；
在公式(5)的基础上，根据公式(6)计算最大的关系时间步长Δtmax所对应的最小分段数N0min为：
优选地，通过最小分段数N0min计算其他任意一个有压管段i的最小波速调整幅度fimin及最小波速调整幅度fimin所对应的最优分段数Ni最优，具体按照下述步骤实现：
S1，采用公式(7)计算任意一个有压管段i的关系时间步长Δtj；
△tj＝L0/(N0mina0(1+f0))，-0.15≤f0≤0.15 (7)；
L0、a0分别表示压力管道输水系统内水锤波传播时间最短有压管段的管长和波速；f0表示压力管系内波传播时间最短管段的波速调整系数；当f0在[-0.15,0.15]范围内以一微小增量j改变时，由公式(7)计算当f0＝f0j时，与f0j对应的Δtj值，所述微小增量j的取值增量为0.001；
S2，根据关系时间步长Δtj，计算压力管道输水系统内编号为i的有压管段的分段数Ni，i＝1,2,…,m，m表示压力管道输水系统中除水锤波传播时间最短有压管段外有压管段的总数量；判断分段数Ni是否为式(8)所示的范围内的正整数；
如果是，进入S3，计算得到编号为i的有压管道的最小波速调整幅度fimin：
如果否，表示f0＝f0j求得的关系时间步长Δtj不符合要求，则返回S1，应用公式(7)计算当f0＝f0(j+1)时，与f0(j+1)所对应的△t(j+1)，继续S2，判断通过f0＝[-0.15,0.15]中的所有取值是否能够计算得到符合式(8)的正整数分段数；如果是，则进入S3；如果否，则返回S1，并将公式(7)中的N0min修改为N0min+1，继续计算，直至计算得到编号为i的管段至少存在一个符合式(8)的正整数分段数为止；
S3，将符合式(8)的编号为i的有压管道分段数按下式(9)进行计算，得到编号为i的有压管道的最小波速调整幅度fimin；
fimin所对应的分段数即为有压管段i的最优分段数Ni最优，Li表示有压管段i的管长，ai表示管段i的波速。
优选地，通过最小波速调整幅度fimin，调整水锤波速，保证调整前、后，水锤波从压力管道输水系统起点传播到压力管道输水系统终点的总时间相同，具体按照下述步骤实现：
水锤波速未调整，水锤波从压力管道输水系统起点传播到压力管道输水系统终点的实际总时间为T，水锤波速调整后，水锤波从压力管道输水系统起点传播到压力管道输水系统终点的总时间为T′，T与T′的关系用公式(10)表示：
T'ε＝T (10)；
其中，系数ε符合公式(11)；
则，最优波速调整系数为f′i：f′i＝(1+fimin)/ε-1 (12)；
最优时间步长为Δt′：△t′＝△tjε (13)；
最后，判断压力管道输水系统中除水锤波传播时间最短有压管段外任意一个有压管段的最优波速调整系数为fi′是否小于0.10，如果是，则该有压管段的最优时间步长为Δt′，如果否，则将水锤波传播时间最短有压管段的分段数加1后，重新计算最大的时间步长Δtmax及最短有压管段的最小分段数N0min，然后进入S2，直至计算得到的所有结果均满足条件为止。
优选地，调整水锤波速调整前、后的总时间相同后，还包括以下步骤：
计算与有压管道相邻的无压渠道的时间步长；设无压管道的时间步长为△t0，管道流的时间步长为△tc，令△t0与△tc符合公式(14)：
△tc＝N△t0，N为整数 (14)；
设明渠和管道流的空间步长分别为△x0与△xc，为保证明渠与管道流的衔接，将明渠在连接断面附近的空间步长△x0设置为△x0/N，然后按库朗条件将明渠在连接断面的△t0缩小1/N×△t0，从而与管道的时间步长△tc相匹配。
本发明是一种通过改进调整波速法压力管道确定时间步长的压力管道输水系统非恒定流模型的优化技术，采用改进的调整波速法进行分段，首先对波传播时间最短管段的波速进行调整，选取尽可能大的时间步长Δt，最小化整个管系波速调整幅度，再通过各管段最优波速调整系数，使波速调整前后水锤波从压力管道起点传播到终点的总时间保持一致。
本发明的有益效果是：本发明所述方法使压力管道输水系统的波速调整幅度与比传统方法相比显著减小；通过ε的调整,使时差ΔT＝0s，保证了波速调整前后水锤波从压力管系起点传播到终点的总时间完全一致。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施方式仅仅用以解释本发明，并不用于限定本发明。
实施例
本实施例所述压力管道输水系统非恒定流模型的优化方法，所述方法包括：
首先，获取压力管道输水系统中的有压管段，对水锤波传播时间最短有压管段的水锤波速进行调整，获取最大时间步长Δtmax，计算获取最大时间步长Δtmax所对应的最短有压管段的最小分段数N0min；
然后，通过最小分段数N0min计算其他任意一个有压管段i的最小波速调整幅度fimin及最小波速调整幅度fimin所对应的最优分段数Ni最优；
最后，通过最小波速调整幅度fimin，调整水锤波速，保证调整前、后水锤波从压力管道输水系统起点传播到压力管道输水系统终点的总时间相同。
更详细的解释说明：
(一)对水锤波传播时间最短有压管段的水锤波速进行调整，获取最大时间步长Δtmax，计算获取最大时间步长Δtmax所对应的最短有压管段最小分段数，具体按照下述步骤实现：
计算压力管道输水系统内水锤波传播时间最短有压管段的管长L0、波速a0和分段数N0，其中，管长L0、波速a0存在关系式(1)：
压力管道输水系统中除水锤波传播时间最短有压管段外任意一个有压管段分段数Ni满足公式(2)，所述Ni为正整数:
其中，i表示压力管道输水系统中除水锤波传播时间最短有压管段外任意一个有压管段的编号，i为大于等于1的正整数，m表示压力管道输水系统中除水锤波传播时间最短有压管段外有压管段的总数量，α表示压力管道输水系统中除水锤波传播时间最短有压管段外其他有压管段的水锤波波速；fi表示编号为i的有压管段的波速调整系数；Li表示编号为i的有压管段的长度；△t表示水锤波传播时间最短管段的关系时间步长；
所述△t的取值范围符合公式(3)：
fmax表示|fi|的最大取值，在计算水锤波传播时间最短有压管段的波速调整的关系时间步长时，|fi|＝0.15；T表示未调整波速时，水锤波从压力管道输水系统起点传播到终点的实际总时间；
由于水锤波传播时间最短有压管段的波速调整不会影响对压力管道输水系统的计算结果，故，在获取△t的取值范围时，设|fi|＝0.15，则fmax＝-0.15或fmax＝0.15，因此，△t的取值范围符合公式(4)：
通过公式(4)，得到最大关系时间步长Δtmax的计算公式(5)；
在公式(5)的基础上，根据公式(6)计算最大的关系时间步长Δtmax所对应的最小分段数N0min为：
(二)通过最小分段数N0min计算其他任意一个有压管段i的最小波速调整幅度fimin及最小波速调整幅度fimin所对应的最优分段数Ni最优，具体按照下述步骤实现：
S1，采用公式(7)计算任意一个有压管段i的关系时间步长Δtj；
△tj＝L0/(N0mina0(1+f0))，-0.15≤f0≤0.15 (7)；
L0、a0分别表示压力管道输水系统内水锤波传播时间最短有压管段的管长和波速；f0表示压力管系内波传播时间最短管段的波速调整系数；当f0在[-0.15,0.15]范围内以一微小增量j改变时，由公式(7)计算当f0＝f0j时，与f0j对应的Δtj值，所述微小增量j的取值增量为0.001；
S2，根据关系时间步长Δtj，计算压力管道输水系统内编号为i的有压管段的分段数Ni，i＝1,2,…,m，m表示压力管道输水系统中除水锤波传播时间最短有压管段外有压管段的总数量；判断分段数Ni是否为式(8)所示的范围内的正整数；
如果是，进入S3，计算得到编号为i的有压管道的最小波速调整幅度fimin：
如果否，表示f0＝f0j求得的关系时间步长Δtj不符合要求，则返回S1，应用公式(7)计算当f0＝f0(j+1)时，与f0(j+1)所对应的△t(j+1)，继续S2，判断通过f0＝[-0.15,0.15]中的所有取值是否能够计算得到符合式(8)的正整数分段数；如果是，则进入S3；如果否，则返回S1，并将公式(7)中的N0min修改为N0min+1，继续计算，直至计算得到编号为i的管段至少存在一个符合式(8)的正整数分段数为止；
S3，将符合式(8)的编号为i的有压管道分段数按下式(9)进行计算，得到编号为i的有压管道的最小波速调整幅度fimin；
fimin所对应的分段数即为有压管段i的最优分段数Ni最优，Li表示有压管段i的管长，ai表示管段i的波速。
(三)通过最小波速调整幅度fimin，调整水锤波速，保证调整前、后，水锤波从压力管道输水系统起点传播到压力管道输水系统终点的总时间相同，具体按照下述步骤实现：
水锤波速未调整，水锤波从压力管道输水系统起点传播到压力管道输水系统终点的实际总时间为T，水锤波速调整后，水锤波从压力管道输水系统起点传播到压力管道输水系统终点的总时间为T′，T与T′的关系用公式(10)表示：
T'ε＝T (10)；
其中，系数ε符合公式(11)；
则，最优波速调整系数为fi′：fi'＝(1+fimin)/ε-1 (12)；
最优时间步长为Δt′：△t′＝△tjε (13)；
最后，判断压力管道输水系统中除水锤波传播时间最短有压管段外任意一个有压管段的最优波速调整系数为fi′是否小于0.10，如果是，则该有压管段的最优时间步长为Δt′，如果否，则将水锤波传播时间最短有压管段的分段数加1后，重新计算最大的时间步长Δtmax及最短有压管段的最小分段数N0min，然后进入S2，直至计算得到的所有结果均满足条件为止。
(四)调整水锤波速调整前、后的总时间相同后，还包括以下步骤：
计算与有压管道相邻的无压渠道的时间步长；设无压管道的时间步长为△t0，管道流的时间步长为△tc，令△t0与△tc符合公式(14)：
△tc＝N△t0，N为整数 (14)；
设明渠和管道流的空间步长分别为△x0与△xc，为保证明渠与管道流的衔接，将明渠在连接断面附近的空间步长△x0设置为△x0/N，然后按库朗条件将明渠在连接断面的△t0缩小1/N×△t0，从而与管道的时间步长△tc相匹配。
通过采用本发明公开的上述技术方案，得到了如下有益的效果：本发明所述方法使压力管道输水系统的波速调整幅度与比传统方法相比显著减小；通过ε的调整,使时差ΔT＝0s，保证了波速调整前后水锤波从压力管系起点传播到终点的总时间完全一致。
以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视本发明的保护范围。