摘要: 根据1976年唐山7.8级地震的前兆资料，结合岩石破裂物理、震源物理、构造物理、地震前兆、地震力源、地壳结构等方面的研究结果，采用含有菱形硬包体的饱和水孔隙介质模型，数值模拟了唐山大地震的孕育过程．首先导出了考虑岩石非弹性体积膨胀和应变软化、孔隙水渗透的固-液两相介质的孕震动力学方程组，提出了用增量形式求解该方程组的有限差分法，然后计算模拟了硬包体内含有断层软弱带时的唐山大震的孕育过程，及其时空变化特征．数值模拟结果表明，唐山大地震的孕育主要经历了弹性积累、早期非弹性体积膨胀、断层第1次加速软化(加速破坏和蠕滑)、天津-宁河一带第2次非弹性膨胀、断层第2次加速软化(临震加速破坏和蠕滑)和发震等过程．这里，与断层的加速破坏相对应，孕震介质的某些部分出现了非弹性膨胀的弹性恢复和二次非弹性膨胀．与前兆观测资料比较分析可以发现，上述模拟结果可以比较好地解释实际观测资料的变化特征．同时，数值模拟还进一步揭示了唐山地震孕育场的时空变化规律和复杂图象，为地震前兆场变化的物理机制的研究奠定了基础．