四川省微信网站建设推广,扶贫网站建设方案,美容网站设计,wordpress always主题铁路轨道不平顺数据分析与预测
1.引言
铁路轨道作为铁行车的基础设施#xff0c;是铁路线路的重要组成部分。随着经济和交通运输业的发展#xff0c;我国的铁路运输正朝着高速和重载方向迅速发展#xff0c;与此同时#xff0c;轨道结构承受来自列车荷载、运行速度的冲击…铁路轨道不平顺数据分析与预测
1.引言
铁路轨道作为铁行车的基础设施是铁路线路的重要组成部分。随着经济和交通运输业的发展我国的铁路运输正朝着高速和重载方向迅速发展与此同时轨道结构承受来自列车荷载、运行速度的冲击和列车的振动等各方面的作用力不断增大不仅加速了铁路轨道设备的损坏由此产生的轨道不平顺问题会严重影响车辆行乘客的舒适度以及设备的使用寿命等存在非常严重的安全隐患。 在铁路运营过程中轨道在列车不稳定荷载的反复作用下容易发生一定的几何形变这其中有些是垂直向和横向的动态弹性变形有些则是永久变形而这些现象都统称为轨道不平顺 轨道的不平顺会引起列车的振动和轮轨作用力不断变化不仅影响列车的舒适度降低轨道与列车部件的使用寿命甚至可能导致列车脱轨严重威胁行车安全。其中轨道的高低不平顺容易激起列车的垂向振动导致轮轨间作用力大幅减小甚至造成轮轨分离使得列车存在悬浮和脱轨的危险。
本项目需要对高速铁路实际运营过程中测量的动检数据进行分析得到滑动轨道质量指数TQI。进行高低不平顺劣化预测根据预测结果对模型的预测效果进行评价。 项目目标 1.分析所给轨道动检数据的滑动轨道质量指数并结合相应规范给出轨道质量验收报告。 2.根据计算得到的高低不平顺标准差构建时间序列预测分析模型预测未来一年的高低不平顺劣化趋势并采用均方误差和均方根误差对预测模型的结果进行评价。
2.数据预处理
使用python中的pandas库函数读入数据并对数据进行初步的处理。 对于2016年1月至2021年12月共计12个月的月度数据读取至dataframe结构如图2.1所示一共有里程(m)轨距(mm)超高(mm)左轨向(mm)右轨向(mm)左高低(mm)右高低(mm)三角坑(mm)水平(mm)共计9个属性。 画出数据分布的箱线图使用箱线图判断是否存在在可信范围之外的异常数据。箱形图Box-plot又称箱线图是一种表征一组数据分散情况的统计图主要用于体现原始数据分布的特征还可以进行多组数据分布特征的比较。箱形图包括最小值min下四分位数Q1中位数Xm上四分位数Q3和最大值max如图 3-2 所示绿色矩形框的上下两端边的位置分别对应数据的上、下四分位数Q3、Q1矩形框内部的蓝色线段为中位线对应中位数Xm。上四分位数与下四分位数之间的距离是数据分布的一种简单度量它给出被数据的中间一半所覆盖的范围。该距离称为四分位数极差IQR。
因此最大值max定义为 Q31.5IQR最小值min定义为 Q3-1.5IQR。基于此视大于最大值max或小于最小值min的数据为离群点。
图2.2 数据箱线图
由图2.2可以观察到数据的轨距水平左高低等属性存在较多的离群点即异常数据这里采用最大最小值异常值替换的方式进行异常值处理。处理之后的箱线图如图2.3所示所有数据都在可信范围之内。
3.滑动TQI的计算
轨道质量指数是评价单元区段整体不平顺的重要指标。目前我国釆用的计算方法主要是通过计算200米单元区段内轨距左右高低左右轨向水平超高以及三角坑等共7个单项的不平顺幅值标准差之和来获得的能综合反映该区段的整体不平顺状态以及轨道恶化程度。通过对所有轨道单元区段的TQI数值排序可以确定需要重点维修的区段所以TQI是我国铁路工务部门指导线路维修最重要的指标。
根据工务部门多年的管理经验表明对轨道区段TQI进行长期的监测和分析有利于现场维修人员及时掌握轨道状态的改善或者恶化情况以及未来的发展趋势从而更有效率的编制维修计划。
TQI 以 200 m 长度轨道区段作为计量单元对单元区段内的轨道几何进行统计用标准差来表示单 项轨道几何不平顺状态而 TQI 则为一个单元区段内左高低、右高低、左轨向、右轨向、轨距、水平和 三角坑等七个单项几何不平顺标准差之和计算方法如下列公式 由上述TQI的定义可以得到滑动TQI的计算方法。以1为滑动TQI的移动步长以800为滑动TQI的移动窗长。以2016年1月的测量数据为例可以得到滑动TQI结果图见图3.1
画出2020年1月的TQI沿里程变化图与2016年1月滑动TQI沿里程变化图作比较 由图3.2可以得出从2016年至2020年轨道的区段的状态发生了明显恶化且在部分区段恶化非常显著。
4.高低不平顺劣化预测
在列车运行的过程中因轨道不平顺而产生的列车振动加速度对货车以及客车的运行均有影响主要的影响货车的安全性客车的舒适性等。由于近年来铁路大面积提高运行效率朝着高速和重载方向的发展对轨道结构提出了更高的要求除了应该具备必要的强度外还要保持良好的几何形位。为此国内外铁路部门及研究机构都加强了对轨道几何变形的检测。
轨道几何不平顺的变化特征是轨道上大量样本检测点的几何不平顺随时间及车俩载荷重复作用下表现出共有或者相似的特征反映了轨道几何不平顺的一般规律。本章采用SARIMA模型进行轨道高低不平顺劣化预测。
我们以K44074段的左高低为例进行研究首先画出K44074段的左高低标准差在2016年1月至2020年12月的变化趋势如图所示 我们选取2016年1月至2019年12月的数据作为训练集2020年全年的数据作为验证集使用SARIMA模型进行建立预测模型。 得到ACF与PACF图像如图所示 由图4.1可知最佳的SARIMA模型为SARIMA (0, 1, 1)。 使用SARIMA(0, 1, 1)模型对2020年轨道高低不平顺标准差进行预测可以得到结果如所示 参考代码
135. TQI[]
136. for i in range(0,len(data1)-800):
137. temp_datadata.loc[data.index[i:i800]]
138. a1np.array(temp_data[三角坑(mm)])
139. a2np.array(temp_data[右轨向(mm)])
140. a3np.array(temp_data[右高低(mm)])
141. a4np.array(temp_data[左轨向(mm)])
142. a5np.array(temp_data[左高低(mm)])
143. a6np.array(temp_data[轨距(mm)])
144. a7np.array(temp_data[水平(mm)])
145. rea1.std()a2.std()a3.std()a4.std()a5.std()a6.std()a7.std()
146. TQI.append(re)
147. import matplotlib.pyplot as plt
148. import numpy as np
149. import matplotlib
150. matplotlib.rc(font,familyKaiTi)
151.
152. xnp.arange(0,4800,0.25)
153.
154. plt.title(2016年1月滑动TQI沿里程变化图)
155. plt.xlabel(里程(m))
156. plt.ylabel(TQI)
157.
158. plt.plot(x,TQI)
159.
160. # 2020年1月滑动TQI沿里程变化的序列数据
161. TQI_2[]
162. for i in range(len(data1)*48,len(data1)*49-800):
163. temp_datadata.loc[data.index[i:i800]]
164. # python中对于矩阵的计算尽量不要在DataFrame中效率非常低。应该使用numpy
165. a1np.array(temp_data[三角坑(mm)])
166. a2np.array(temp_data[右轨向(mm)])
167. a3np.array(temp_data[右高低(mm)])
168. a4np.array(temp_data[左轨向(mm)])
169. a5np.array(temp_data[左高低(mm)])
170. a6np.array(temp_data[轨距(mm)])
171. a7np.array(temp_data[水平(mm)])
172. rea1.std()a2.std()a3.std()a4.std()a5.std()a6.std()a7.std()
173. TQI_2.append(re)
174.
175. xnp.arange(0,4800,0.25)
176.
177.
178.
179. plt.plot(x,TQI,label2016年1月)
180. plt.plot(x,TQI_2,label2020年1月)
181.
182. plt.title(滑动TQI沿里程变化图)
183. plt.xlabel(里程(m))
184. plt.ylabel(TQI)
185. plt.legend()
186.
187. input_data[]
188. for i in range(72):
189. temp_datadata.loc[data.index[len(data1)*i16000:len(data1)*i20000]]
190. tempnp.array(temp_data[左高低(mm)])
191. input_data.append(temp.std())
192.
193. from matplotlib.pyplot import MultipleLocator
194. #从pyplot导入MultipleLocator类这个类用于设置刻度间隔
195.
196. fig plt.figure(figsize(15,5)) # 设置画布大小
197. plt.ylim(0.1,1) # y轴显示范围
198. plt.plot(input_data[0:60],o-)
199.
200. from statsmodels.graphics.tsaplots import plot_acf, plot_pacf # 画图定阶
201.
202. def draw_acf(data):
203. # 利用ACF判断模型阶数
204. plot_acf(data)
205. plt.title(序列自相关图(ACF))
206. plt.show()
207.
208. def draw_pacf(data):
209. # 利用PACF判断模型阶数
210. plot_pacf(data)
211. plt.title(序列偏自相关图(PACF))
212. plt.show()
213.
214. def draw_acf_pacf(data):
215. f plt.figure(facecolorwhite)
216. # 构建第一个图
217. ax1 f.add_subplot(211)
218. # 把x轴的刻度间隔设置为1并存在变量里
219. x_major_locator MultipleLocator(1)
220. plot_acf(data, axax1)
221. # 构建第二个图
222. ax2 f.add_subplot(212)
223. plot_pacf(data, axax2)
224. plt.subplots_adjust(hspace0.5)
225. # 把x轴的主刻度设置为1的倍数
226. ax1.xaxis.set_major_locator(x_major_locator)
227. ax2.xaxis.set_major_locator(x_major_locator)
228. plt.show()
229.
230. draw_acf_pacf(np.array(train_data))
全部代码见https://download.csdn.net/download/weixin_44026026/87657283
