主权项 |
1.一种基于遗传算法的Kriging插值的地面沉降空间监控方法,其特征在于:该方法包括如下步骤:(1)采集监测点的坐标(x<sub>i</sub>,y<sub>i</sub>),i=1,2,...,N,以及地面沉降量Z<sub>i</sub>,进行规范化处理;(2)计算监测点之间的距离:<img file="DEST_PATH_FDA00002033352700011.GIF" wi="605" he="88" />式中,i=1,2,...,N;j=1,2,...,N;(3)划分距离组:<img file="DEST_PATH_FDA00002033352700012.GIF" wi="929" he="135" />式中,N<sub>H</sub>表示距离组的个数,N<sub>H</sub>≥4;(4)计算各距离组所对应的变异函数值r*(h′<sub>m</sub>):<img file="DEST_PATH_FDA00002033352700013.GIF" wi="1122" he="135" />式中,N(h′<sub>m</sub>)表示相隔距离矢量为h的所有监测点对的个数,Z<sub>i</sub>表示地面沉降量;(5)选择指数变异函数模型并基于遗传算法求定所述模型中的参数,求得变异函数r(h),指数变异函数模型为:r(h)=C<sub>0</sub>+C(1-e<sup>-h/a</sup>),式中,C<sub>0</sub>为块金值,C为偏基台值,a表示变程。(6)计算预测点P邻域内参与预测的参考点之间及各参考点与预测点P的变异函数值γ(i,j)和<img file="DEST_PATH_FDA00002033352700014.GIF" wi="158" he="67" />(7)计算预测点Kriging插值的预测值:<img file="DEST_PATH_FDA00002033352700015.GIF" wi="522" he="70" />其中,λ表示Kriging权系数,λ=K<sup>-1</sup>M,<img file="DEST_PATH_FDA00002033352700016.GIF" wi="622" he="365" /><img file="DEST_PATH_FDA00002033352700017.GIF" wi="251" he="393" />(8)重复步骤(1)至步骤(7),以网格方式遍历沉降监测结果。 |