Method and system for subsurface imaging using multi-physics joint migration inversion and geophysical constraints

A system, comprising:

• a reservoir simulator comprising a computer processor, wherein the reservoir simulator is configured to perform a method comprising:
• obtaining seismic data for a geological region of interest;
• obtaining first geophysical data for the geological region of interest;
• determining a first plurality of pressure wavefields using a slowness model and a reflectivity model;
• determining a plurality of slowness gradients for the slowness model based on the first geophysical data, the first plurality of pressure wavefields, the seismic data, and a first geophysical constraint,
• wherein the first geophysical constraint corresponds to an objective function that determines a misfit between the first geophysical data and the slowness model;
• updating the slowness model using the plurality of slowness gradients to produce an updated slowness model;
• and
• generating, based on the updated slowness model and the reflectivity model, a subsurface image of the geological region of interest.