1.一种煤矿井下二维矿震波速反演降维方法,其特征在于,包括:使用直线将目标区域等间隔划分成m行n列的正方形网格,其中m和n为奇数,其中待反演区域为PT,模型空间MT的维数为MTnum=m×n;
获得目标区域的待反演向量v=(v11,v12,…,v1n,v21,v22,…,vij,…,vmn),所述目标区域的每个正方形的介质均匀,所述目标区域的波速为vij,i∈[1,m],j∈[1,n];
根据所述待反演向量获得所述目标区域的每个正方形的波速上限向量和波速下限向量,所述波速上限向量和所述波速下限向量分别为:获得降维之后的待反演变量mod_angvector=(modvector,angvector),其中所述待反演区域PT的频域为PF,所述频域PF的维数为MFnum=MTnum=m×n,所述待反演变量mod_angvector的维数为MFCnum,根据modvector和angvector形成的频域为PFC,频域PF的中心点坐标为 频域PF的最远点坐标为获得所述待反演变量mod_angvector的上限向量为 其中modvector和angvector的每个元素上限构成的向量分别为:获得所述待反演变量mod_angvector的下限向量为 其中modvector和angvector的每个元素下限构成的向量分别为:当||tobs-ttheo||2最小时,根据所述上限向量 和所述下限向量获得待反演变量mod_angvector,其中tobs为测得到时差向量,ttheo为理论到时差向量。
2.根据权利要求1所述的煤矿井下二维矿震波速反演降维方法,其特征在于,所述根据所述待反演向量获得所述目标区域的每个正方形的波速上限向量和波速下限向量的步骤之后,所述获得降维之后的待反演变量的步骤之前,包括:对所述待反演向量、所述波速上限向量、所述波速下限向量进行归一化处理,获得归一化之后的待反演向量、波速下限向量、波速上限向量分别为:v=(v11,v12,…,v1n,v21,v22,…,vij,…,vmn)根据归一化之后的待反演向量、波速上限向量、波速下限向量获得归一化系数为:其中,min(vmin)为向量vmin中最小的元素,max(vmax)为向量vmax中最大的元素,min(vmin)为向量vmin中最小的元素,max(vmax)为向量vmax中最大的元素。
3.根据权利要求2所述的煤矿井下二维矿震波速反演降维方法,其特征在于,还包括:根据所述归一化系数获得如下公式:
vij=knorm·vij (2)
根据上述公式(2)获得波速归一化公式为:
v=knorm·v (3)。
4.根据权利要求1所述的煤矿井下二维矿震波速反演降维方法,其特征在于,所述频域PF之中拟保留的变量数目为MFCnum=4D2+4D+1。
5.根据权利要求1所述的煤矿井下二维矿震波速反演降维方法,其特征在于,所述理论到时差向量为:ttheo=RAY_TRACE(FFT_2D-1(MAP_TO_VECTOR-1(mod_angvector))) (4)其中,FFT_2D-1为根据所述待反演区域PT获得所述频域PF的逆过程,MAP_TO_VECTOR-1为根据所述频域PFC获得待反演变量mod_angvector的逆过程,RAY_TRACE为根据射线追踪获得理论到时差向量ttheo。