A computer-implemented method for obtaining reconstructed seismic data for determining a subsurface feature, includes: determining an initial training velocity model, training a machine learning model based on first training seismic data and second training seismic data generated from the training velocity model, the first training seismic data corresponding to one or more first frequencies, the second training seismic data corresponding to one or more second frequencies lower than the one or more first frequencies, obtaining, based on measured seismic data and the machine learning model, reconstructed seismic data corresponding to the one or more second frequencies, generating a velocity model based on the measured seismic data, the reconstructed seismic data, and a full waveform inversion (FWI), and when the generated velocity model does not satisfy a preset condition, updating the training velocity model based on the generated velocity model, to obtain updated reconstructed seismic data for determining a subsurface feature.