融合InSAR和GA的开采沉陷预计参数反演模型研究

针对基于常规合成孔径雷达差分干涉测量（differential interferometry synthetic aperture radar，D-InSAR）技术成果无法反演全部概率积分参数问题，开展了融合D-InSAR和遗传算法（genetic algorithm，GA）的概率积分参数反演方法研究。①依据采动区D-InSAR视线（line of sight，LOS）方向变形为下沉和南北、东西方向水平移动沿LOS方向投影的关系，基于GA理论，首次构建了融合D-InSAR和GA的概率积分法预计全参数反演模型，并编制了求参软件。模拟实验结果表明，q、tanβ、b、θ的相对误差不超过6.7%，拐点偏移距S相对误差不超过20%（参数敏感度低，对求参整体效果影响小）；求参拟合下沉误差为-5.90~6.10 mm，拟合中误差约为±2.20 mm，整体求参效果较好。②利用融合D-InSAR和GA的概率积分法预计全参数反演模型，在一定的假设条件下，对2012年3月11日9310工作面开采沉陷的动态概率积分参数进行了求取，得到q=0.172，b=0.13，tanβ=2.08，θ=88°，S₁=-10 m，S₂=9 m，S₃=-85 m，S₄=30 m，拟合误差约为-35.00~45.00 mm，拟合中误差为±15.16 mm，并认为S₃=-85 m主要为9310工作面开采导致邻近9312采空区边界的悬臂梁或砌体梁失稳"活化"所致。

Parameter Inversion Model for Mining Subsidence Prediction Based on Fusion of InSAR and GA

Considering that conventional differential interferometry synthetic aperture radar (D-InSAR) techniques are impossible to calculate all the parameters of probability integral method, this paper studies the parameters inversion of probability integral method based on fusion of D-InSAR and genetic algorithm(GA). ①In mining area, the deformation along line of sight(LOS) direction by D-InSAR is the projection of sinking, north-south and east-west horizontal movement along LOS direction. Based On GA theory, the probability integral full parameters inversion model is firstly established based on D-InSAR and GA. Moreover, the parameter calculation software is compiled. The simulation experimental results show that the relative errors of q, tanβ, b, θ are not more than 6.7%, the relative error of the inflection point offset S is not more than 20%. The fitting error is between -5.90-6.10 mm, and the fitting mean square error is within±2.20 mm. The effect of the whole parameter calculation is better. ②Taking advantage of the probability integral full parameters inversion model based on fusion of D-InSAR and GA, under certain assumptions, we obtain the dynamic predicting parameters of mining subsidence in 9310 working face on Mar. 11, 2012. The results are as follows:q=0.172, b=0.13, tanβ=2.08, θ=88°, S₁=-10 m, S₂=9 m, S₃=-85 m, S₄=30 m, the fitting error is -35.00-45.00 mm, and the fitting mean square error is within ±15.16 mm. The mining of 9310 working face causes the instability of cantilever beam or masonry beam in the boundary of the adjacent 9312 goaf, thus making S₃ be -85 m. The research results have significant reference value to monitoring mining deformation by D-InSAR.