This study applies the Load/Unload Response Ratio(LURR)method to identify pre-seismic gravity anomalies using high-precision superconducting gravity observations from the Yilan and Nianzishan stations as response variables. the results show that detectable LURR anomalies appeared before the studied earthquakes when the preprocessed gravity data were obtained by subtracting only the theoretical solid Earth tide. the occurrence of these anomalies was not substantially affected by whether the linear instrumental drift was removed. This suggests that actual tidal and barometric perturbations may act as external loading factors related to earthquake triggering. At the Yilan station,the background LURR value was approximately 1.0. During the 1-3 months preceding the Tangwang M_L4.0 earthquake,the superconducting-gravity-derived LURR anomaly showed distinct peak-like variations. the inferred preparatory process of the Tangwang earthquake may have involved regional tectonic deformation that generated fractures in highly stressed rocks,allowing deep fluids to migrate into the fractured zones. As stress accumulation approached a critical state,weak tidal stress or barometric-pressure perturbations may have promoted rock damage and instability,ultimately contributing to earthquake occurrence. At the Nianzishan station,the background LURR value was also approximately 1.0. Preseismic LURR anomalies associated with the Tailai M_L3.6 earthquake similarly exhibited peak-like variations. the preparation process of the Tailai earthquake may have been related to tectonic compression of near-surface and subsurface media. Near the epicentral area,particularly along the margins of high-compression zones,the rock medium may have become unstable and fractured,thereby promoting earthquake occurrence. Overall,this retrospective case study indicates that the superconducting-gravity-based LURR method can extract potential preseismic anomalies associated with small-to-moderate earthquakes in Heilongjiang Province. the results suggest that this method has practical value for retrospective earthquake-case verification and may provide a useful observational constraint for regional earthquake anomaly tracking.