A three-dimensional numerical model based on the smoothed particle hydrodynamics (SPH) method is utilized to model such interactions. The effects from air and water are well captured using suitable force components and without explicit inclusion of air and water phases. This reduces the four-phase SPH model into a two-phase model, significantly reducing computational costs and potentially enabling the use of this model for large-scale simulations. This study includes two simulation domain size. The results show that the remained oil slick in the icy domain is completely depend on properties of different crude oils and current velocity. Then we performed a numerical experiment to study the effects of different ice parameters (e.g., ice concentration, floe shape) on the drift of fresh and weathered oil. Ice fields with the size of 750m$\times$1350 m are used in the experiment. The ice concentration varies between 30%, 50% and 70%. Quadrilateral and circular ice floes, in addition to real-shape floes digitize from an aerial image are all used in the experiment. Oil spills of 1257 ${\mathrm{m}}^3$ are released in the open water and later drift into the ice under the influence of steady current of 50 $\frac{\mathrm{cm}}{\mathrm{s}}$. The type of spilled crude oil varies from fresh Troll B, Troll B 250 °C, and Troll B 250 °C and 75%emuls. The behavior of oil spills while crossing through the ice infested area are investigated by measuring the thickness of crude oil, the average and maximum distance travelled by crude oil.

基于平滑粒子流体动力学（SPH）方法的三维数值模型被用来模拟这种相互作用。使用适当的力分量可以很好地捕获空气和水的影响，并且无需明确包含空气和水相。这将四相 SPH 模型简化为两相模型，显着降低了计算成本，并有可能将该模型用于大规模模拟。本研究包括两个模拟域大小。结果表明，冰域中残留的浮油完全取决于不同原油的性质和海流速度。然后我们进行了数值实验，研究不同冰参数（例如冰浓度、浮冰形状）对新鲜和风化石油漂移的影响。冰原面积达750m$\times$实验中使用1350 m。冰浓度在30%、50%和70%之间变化。实验中使用了四边形和圆形的浮冰，以及从航空图像中数字化的真实形状的浮冰。1257 起漏油事件${\mathrm{米}}^3$被释放到开放水域，然后在 50 的稳定电流的影响下漂入冰中$\frac{\mathrm{厘米}}{\mathrm{s}}$。泄漏原油的类型各不相同，有新鲜的 Troll B、Troll B 250 °C、Troll B 250 °C 和 75% 乳化液。通过测量原油厚度、原油移动的平均和最大距离来研究穿过冰覆盖区域时的石油泄漏行为。