Los métodos de campos potenciales, basados en la explotación de campos gravitacionales y magnéticos, son uno de los métodos más importantes para recuperar información fundamental sobre la estructura de la corteza terrestre a escalas globales, regionales y locales. El cuello de botella para este tipo de métodos geofísicos a menudo está representado por el desarrollo de técnicas para explotar completamente los datos disponibles. De hecho, cada técnica diferente puede observar el efecto de una sola propiedad del subsuelo y cuando queremos estimar esta propiedad a partir del campo observado (el llamado problema inverso), surgen varios problemas como la no unicidad y la inestabilidad. Una posible solución a estos problemas consiste en invertir conjuntamente, de manera consistente, diferentes campos observados, incorporando posiblemente todas las restricciones geológicas disponibles. En el trabajo actual, presentamos un innovador algoritmo bayesiano destinado a realizar una inversión conjunta completa en 3D de campos gravitacionales y magnéticos restringidos por información geológica cualitativa. El algoritmo se prueba en un escenario de caso real, a saber, un estudio local para estimar un modelo completo en 3D del complejo carbonatítico de Oka. Este complejo es un plutón compuesto en Quebec (Canadá), importante para las operaciones mineras relacionadas con materias primas críticas como el niobio y otros elementos de tierras raras. Este ejemplo muestra la fiabilidad del algoritmo de inversión desarrollado y da pistas sobre el papel fundamental que los campos potenciales pueden desempeñar en las actividades mineras.