同时研究两个交互项用什么模型？

请问各位老师，我论文中使用了三个核心解释变量x1，x2，x3，基准回归方程是Y=x1+x2+x3+ε。现在想研究x1和x2、x1和x3之间的交互效应，是用Y=x1+x2+x1#x2+ε和Y=x1+x3+x1#x3+ε来分别回归，还是像下面这样放在同一个方程中？Y=x1+x2+x3+x1#x2+x1#x3+ε
还有就是这样做会不会被质疑为什么不放入x2#x3以及x1#x2#x3的交互项？

楼主，问题解决了吗？