The ocean bottom dual-sensor (OBD) observation system places water detectors (pressure sensors) and land detectors on the irregular seabed interface to receive acoustic pressure fields and elastic wave fields, respectively. These data are usually applied to effectively remove ghost waves and ringing interference in the seawater layer, rather than jointly being used for migration. We propose an ocean bottom dual-sensor viscoacoustic and separated-viscoelastic coupled Q-compensated least-squares reverse time migration (OBD-VSV-QLSRTM) in curvilinear coordinates. In this method, to solve the Q-compensated elastic inverse problem, we jointly use the recorded acoustic pressure and elastic multicomponent seismic data from the OBD observation system. A new misfit function in a sense of least-squares inversion is constructed based on viscoacoustic and separated viscoelastic wavefields. Acoustic and separated viscoelastic equations are utilized to describe the acoustic wavefield in the acoustic medium, and Q-compensated P- and S-wavefields in underlying viscoelastic medium, respectively. On the acoustic-viscoelastic interface, We implement a stable governing equation to ensure a continuous and steady transmission of Q-compensated stresses in the subseabed and pressure in the seawater. To address the irregular seabed problem, we mesh the acoustic-viscoelastic models onto curvilinear grids and apply coordinate transformations to map the models and equations to a new curvilinear coordinate system. We derive viscoacoustic and separated-viscoelastic coupled adjoint and demigration operators in the curvilinear coordinates, enabling linear acoustic-viscoelastic modeling and stable receiver-side back-propagated wavefield continuation. Numerical examples conducted on two typical acoustic-viscoelastic models demonstrate that our proposed OBD-VSV-QLSRTM in curvilinear coordinates can produce highly accurate P- and S-wave images efficiently.
