薄怀志，宋炳忠，李川建，王猛，卢国宏，陈宗成，李云伟

(1.山东省地质矿产勘查开发局第三水文地质工程地质大队，山东 兖州 272100；2.山东省地质矿产勘查开发局第二地质大队，山东 兖州 272100；3.山东省地矿局采煤塌陷地综合治理工程技术研究中心，山东 兖州 272100；4.武汉大学 测绘学院， 湖北 武汉 430000)

0 引言

目前概率积分法沉陷预计在“三下采煤”、矿区土地复垦、采煤塌陷地治理规划设计、矿区建筑物稳定性评估等工作中应用广泛[1-4]。近年来，对任意多边形的沉陷预计研究较多，吴侃[5]首先将非矩形工作面划分为若干小矩形工作面，分别进行预计，然后将各个小矩形预计结果叠加的方法；许冬等[6]将任意多边形工作面按照角点划分为若干梯形，分别预计后叠加(减)处理实现沉陷预计；徐孟强等[7]提出了利用三角剖分根据角点将多边形工作面划分为若干三角形、分别预计后对各三角形预计结果叠加(减)的方法。可以看出，对于任意多边形开采区域划分多从角点往内进行剖分，容易造成小的剖分区域重复计算或遗漏，难以真正实现任意多变形开采区域沉陷预计。概率积分法预计参数反演常用模矢法和遗传算法。模矢法是传统的直接寻优算法的典型代表，但该算法对初值可靠性要求较高，初值选取不当，求参结果易陷入局部极值[8-9]。遗传算法是典型的智能算法，但算法求解有一定的不稳定性，即使求参数据完全相同，每次求参结果都会有一定的差异，且算法执行效率低、后期收敛慢[10-11]。……