After some investigation I find a viable way.
Insterestingly, I post reply here to answere myself.
In ADS, we can use 1-port Z-parameter component(Z-box) in eqation based linear components. The Z parameters can be read from external file.
There is a parameter set in the Z-box called Ci. This can be used to model mutual coupling. As for how to represent mutual coupling in a complex invironment, I can only think out the way of series connecting sereval Z-box, each of which have different Ci to model the multiple coupling.
I still hope someone can offer better ways!