First, the main difference remains: given a Lawvere theory you can only takes its models in a cartesian $\infty$-category, while you can look at the models of an operads in any monoidal $\infty$-category. So for example if you want to say something like "commutative ring spectra are commutative monoids in the category of spectra", that, as you are using the smash product of spectra which is not the cartesian product, then you really need the operads for commutative monoids and not the corresponding Lawvere theory (or to put it another way, you need to know that the Lawvere theory for $E_\infty$-algebras is "operadic").
The discussion about modules for an $O$-algebras still applies completely unchanged to the $\infty$-world.
Now, the way $\infty$-operads identifies with special monads/Lawvere theory become much cleaner in the $\infty$ setting: The $\infty$-category of $\infty$-operads identifies exactly with the $\infty$-category of "finitary polynomial monads" (and cartesian morphisms between them, and where finitary is interpreted to mean that their arities are finite sets, not finite spaces). The details of this have been worked out by Gepner, Haugseng and Kock in $\infty$-Operads as Analytic Monads.