为了尽可能避免煤自燃导致的资源浪费和煤矿火灾事故对人民生命和财产带来的危害,进行煤自燃规律的研究,及时预报预测煤自燃的发生是十分必要的,其中使用计算机进行数值模拟是研究煤自燃规律及影响因素的重要方法和手段之一。通过对煤自然发火实验的数值模拟,可以揭示煤田火区的形成机理、发展进程、高温点空间位置转移规律,研究煤田火区自燃、熄灭和复燃特性与规律及其影响因素,为制定防灭火措施提供理论依据。　　本文结合我校煤自然发火实验台的实际情况,根据多孔介质渗流力学和传热传质学理论,通过分析煤自燃过程中气体浓度场、漏风流场、温度场的多场耦合作用,研究煤自燃过程的多场耦合模型,并将耦合模型的控制方程归结为非稳态的对流扩散方程。针对非稳态二维对流扩散方程,构造高精度紧致差分格式,并结合Richardson外推法进行求解;将特征线法与有限体积法结合,构造特征有限体积法,用差分格式进行求解。结合实验和现场实测参数,确定边界条件和初始条件,分别用两种方法对宁夏枣泉矿煤的自然发火实验进行数值模拟,分析煤自燃过程中高温点的发生、发展及变化过程,模拟计算的结果与相同实验条件下的实验结果相符,从而印证了模拟方法的正确性,为替代实验台实验提供了依据。