The production of sand particles associated with the reservoir hydrocarbons becomes one of the most common problems a well of multi cycle steam stimulation may experience during reservoir lifetime. Therefore, it is essential to the investigation of sand particles migration and sand production mechanisms coupled with fluid flow for optimization design of injection and profile control parameters. This paper presents a model based on hydraulic erosion with heat transfer that can analysis release of sand fines under high temperature. The model describes the continuous hydraulic erosion of rock core by the mass conservation equation and particle transport equation, moreover, the sand production and heat transfer are connected through it is coupled with steam injection model. Then, a coupled nonlinear governing equation is derived by using porosity and liquefied sand concentration to evaluate pore enlargement behavior of flow channel and predict sand production behavior under different steam breakthrough channel. Finally, in terms of numerical processing for model, an effectively difference solution method is adopted based on PDE module in COMSOL Multiphysics. The simulation results show that the sand production model can effectively couple heat and sand production, and the relationship between sand production and channeling flow can be studied through seepage process. The degree of steam channeling affects sand production, that is, sand production increases with the increase of the number of steam breakout channels. In addition, it reveals why the sand production can be effectively slowed down by adjusting reasonable injection parameters for multicycle steam stimulation.

与储层碳氢化合物相关的砂粒的生产成为储层寿命期间多循环蒸汽增产井可能遇到的最普遍的问题之一。因此，对于砂粒运移和产砂机理以及流体流动的研究对于注入和剖面控制参数的优化设计至关重要。本文提出了一种基于水力侵蚀和传热的模型，该模型可以分析高温下沙粒的释放。该模型通过质量守恒方程和颗粒输运方程描述了岩心的连续水力侵蚀，并通过与注水模型耦合，将出砂与传热联系起来。然后，利用孔隙度和液化砂浓度推导了非线性控制方程，以评价流道的扩孔行为并预测不同蒸汽穿透通道下的出砂行为。最后，在模型的数值处理方面，采用了COMSOL Multiphysics中基于PDE模块的有效差分求解方法。仿真结果表明，出砂模型可以有效地将热量与出砂耦合，并且可以通过渗流过程研究出砂与窜流的关系。蒸汽窜流的程度会影响出砂量，即，随着蒸汽破裂通道数量的增加，出砂量也会增加。此外，