采空区流场非达西渗流一种新的迭代算法

对Bachmat非线性渗流方程求解问题，结合采空区问题的力学特点，提出基于变渗透率达西(Darcy)渗流求解非线性渗流的迭代算法.建立了冒落采空区非线性漏风流态的有限元数值模型，通过对复杂边界采空区的漏风渗流的计算，得到与实际流态更接近的流动规律(风压等值线和流函数线)和采空区漏风强度分布(速度场).结果表明，迭代算法是振荡性收敛的，收敛速度快；与Darcy渗流相比，采空区非线性渗流速度场趋向于平缓，在工作面附近30\,m范围内最大风速降低到原来的0.61倍，其余大部分区域速度略有增加.迭代方法满足工程要求.