手算的话不准确。如下:

上图300s关闭,水压升高3.2m。

接下来用aft impulse试试。
先建立模型,具体教程见后面的知乎。在建立模型前,需要用excel表格计算出这种管路的水头损失和对应的流量。反算出海曾威廉系数。上述3.6km,管径3m,水损6m左右,流量15m3/s,海曾威廉系数是98

各个节点无非是水库(水池),因为配水池是在阀前,所以要用调压井来代替,因为水库的水位是不变的,只有调压井可以模拟水位变化;然后是关阀规律,这是瞬态,其实一般是按直线关闭。管道的话,就是直径和波速啦。不详细列出了,看教程。
分析时长先按500s。50s直线关阀,管道高程20m,静水位32.24m,压力12.24m,水锤水位42m,压力升高10m,超过50%。

上图,50s关阀,配水井水位达到41m高程
在界面右侧scenario——base scenario下建立一个children,命名为300s,设置阀的关闭时间到300s,如下。修改分析时长为1000s。300s直线关阀,管道高程20m,静水位32.24m,压力12.24m,水锤水位40m,压力升高8m,超过50%。
十分注意:阀门直线关闭起码有3种:蝶阀阀板角度从90°(或0°)关闭到0°的过程与时间是直线关系;阀门的过流面积变化随时间是直线关系;阀门的流量系数随时间是直线关系等。因为阀门关闭特性对水锤影响极大,而这个直线关闭规律必须取得阀门制造企业的曲线。本次简化处理是不正确的。会另外开辟新的教程来探讨。
另一个话题:阀门的流量系数Cv,总的概念是阀前后的水头差和流量的关系。我们可以根据假定的水头损失,反求Cv。而在本例中,当阀门全开时,整个系统的流量基本与表格计算的流量一致,则说明这个阀门的Cv差不多了。也可以通过查看系统纵剖面水头来看看阀门的局部水损是否超出常规。


300s直线关阀,配水井水位40m高程
效果不明显,继续延长关阀时间到1000s,调整分析时长为3000s,水面线也就降到38m。
即1000s直线关阀,管道高程20m,静水位32.24m,压力12.24m,水锤水位38m,压力升高6m,接近50%。

再来回答第二个问题,没有这个配水井情况下,1000s关阀,上游管道末端压力如下:

上图管道的压力波动在一个水锤相8s(波从末端出发,又回到末端的时间)里,从39m变到25m。
有配水井情况下上游管道末端压力如下图,管道压力稍微小一点。相长一点。