The efficient development of shale reservoirs relies on multi-stage hydraulic fracturing technology in horizontal wells involving the difficulty of multi-scale flow for either gas or water phases, which is characterized by the changes in the stress and temperature fields. The original simulation methods are unable to accurately solve the multi-physics field-coupled problems under such complex mechanical conditions or have low stability for two-phase gas-water problems. For this purpose, this study conducted the shale gas production simulation by the meshless finite difference method (GFDM) for the first time to ensure the stability and accuracy of calculations. Bases on the previous multi-sector coupled physical model for shale gas development and proposes a coupled thermo-hydro-mechanical multi-physical field mathematical model considering two-phase gas-water flow. The model is postulated based on different sectors considering various characteristics of gas/water flow, namely desorption attributes, deformation degree and heat transfer characteristics. This method was used to calculate the variation of the multi-physical field and compare the properties of gas-water flow within different reservoir conditions without any non-convergence situation during the calculation. The simulation results indicate that the stress field has a crucial effect on the yield, such that discard of its coupling causes the yield to be overestimated by 12.5%. In contrast, the temperature field has a marginal effect (less than 1%) and can be neglected according to the operating conditions. Furthermore, the stress field significantly characterizes the micro-fracture sector, which is also the major flow-contributor, such that the yield is overestimated by 11% in absence of its coupling, while the effect of stress field in the matrix sector is lower than 1%.

页岩油藏的高效开发依赖于水平井多级水力压裂技术，涉及气水两相多尺度流动的困难，其特点是应力场和温度场的变化。原有的模拟方法无法准确解决如此复杂力学条件下的多物理场耦合问题，或者对气水两相问题的稳定性较差。为此，本研究对页岩气首次采用无网格有限差分法（GFDM）进行生产模拟，确保计算的稳定性和准确性。在以往页岩气开发多部门耦合物理模型的基础上，提出了考虑气水两相流动的热-水-力-多物理场耦合数学模型。该模型是基于考虑气体/水流的各种特征的不同部门假设的，即解吸属性，变形程度和传热特性。该方法用于计算多物理场的变化，比较不同储层条件下的气水流动特性，计算过程中没有出现不收敛的情况。模拟结果表明，应力场对产量有至关重要的影响，丢弃其耦合会导致产量被高估 12.5%。相比之下，温度场具有边际效应（小于 1%），根据运行条件可以忽略不计。此外，应力场显着表征了微裂缝部分，这也是主要的流动贡献者，在没有耦合的情况下，产量被高估了 11%，而基质部分的应力场影响低于1%。