Abstract. On 20 April 2013, Lushan experienced a magnitude 7.0 earthquake. In seismic assessments, borehole strain meters, recognized for their remarkable sensitivity and inherent reliability in tracking crustal deformation, are extensively employed. However, traditional data processing methods encounter challenges when handling massive datasets. This study proposes using a graph wavenet graph neural network to analyze borehole strain data from multiple stations near the earthquake epicenter and establishes a node graph structure using data from four stations near the Lushan epicenter, covering years 2010–2013. After excluding the potential effects of pressure, temperature, and rainfall, we statistically analyzed the pre-earthquake anomalies. Focusing on the Guza, Xiaomiao, and Luzhou stations, which are the closest to the epicenter, the fitting results revealed two accelerations of anomalous accumulation before the earthquake. Approximately four months before the earthquake event, one acceleration suggests the pre-release of energy from a weak fault section. Conversely, the acceleration a few days before the earthquake indicated a strong fault section reaching an unstable state with accumulating strain. We tentatively infer that these two anomalous cumulative accelerations may be related to the preparation phase for a large earthquake. This study highlights the considerable potential of graph neural networks in con-ducting multi-station studies of pre-earthquake anomalies.