无线电波层析成像技术及在矿井坑透中的应用.pdfVIP

第26卷第5期 煤 炭 学 报 Vd．26No．5 CHINACoALSI。CIETY oct． 2001 2001年10月 J吼瓜NAI，0F 文章编号：0253—9993(2001)05—0468—05 无线电波层析成像技术及在矿井坑透中的应用 宁书年1，张绍红1，杨峰1，郭继如2 (1．中国矿业大学北京校区，北京100083；2．开滦矿务局，河北唐山063018) 摘要：通过对无线电波层析成像中的模型构造和反演算法等方面的研究，对模型数据和实测数 据进行了处理．研究结果表明，本文给出的反演算法是有效的．在实测数据处理中，若背景干扰 因素不大，研究的反演算法可以取得比较好的效果．这表明研究算法具有一定的实用价值． 关键词：无线电波；层析成像；反演算法；模型；数据处理 中图分类号：P631．325文献标识码：A 层析成像技术近几年来已逐渐用到煤矿开采、深部隐伏矿床勘探等领域中．在矿井物探中，无线电波 透视法广泛应用于煤矿井下工作面内的地质构造探测，其资料解释方法也从单一的手工方式向地球物理层 析成像技术方向发展． 国外在20世纪70年代对电磁波在地下吸收介质和煤层中的传播进行了大规模理论研究，把层析成像 技术引人到地下电波法资料处理解释当中；80年代又进一步对不连续煤层中的传播理论展开了研究【10J． 国内也进行了地下电波层析成像方法的研究，其算法主要是基于直射线的m汀，SIRT等算法，并提出了 抗干扰的相对层析成像技术【6J． 目前电波成像的主要困难和问题是：(1)井下工作面干扰因素多，观测数据误差大；(2)投影不完 全，每个发射点的观测数据一般只有20个左右；(3)初始场强确定困难．本文根据这些实际情况，首先 对电波层析成像的模型构造进行了研究，之后给出反演当中的一些约束条件，最后采用掷反演算法 对理论模型和实测数据进行了处理． 1电波成像中的模型构造 模型构造就是规定反演的解类和解的表达形式，因此模型的构造与解的适定性有直接的关系．无线电 波层析成像是根据实际观测数据来反演煤层内部的物性参数，一般反演解是一个无限维的函数空间，而实 际测量的数据总是有限的，以有限的观测资料来确定无限函数空间的解，显然求得的解不唯一． 为了保证解的适定性及真解不被排除在外，必须充分利用已知的有关解的各种先验知识或其它约束条 件来构造模型空间． 在以往的电波层析成像当中，对于模型的建立大部分是把所研究的区域分成若干个常系数块，以每块 中系数扰动值作为反演参数．这种模型虽然考虑了横向不均匀性，但由于在块与块之间引入了衰减系数函 数的不连续性，要反演精细的衰减系数结构，块体就必须分得很小，这样就会增加反演参数，同时为了保 证有射线穿过小块体，就必须加密观测，这种模型对于井下工作面进行的探测是不实际的【l~5J．为了克 服上述模型构造的缺陷，本文选取双三次样条函数作为模型函数，即把待反演的区域划分成许多矩形网 格，网格节点的系数为待反演的参数，区域内任一点的衰减系数由包围该点的16个节点的衰减值进行插 值得到[6|．双三次样条函数由B样条组合而成，解的性质由此样条的阶数决定，这样构成的模型函数为 收稿日期：200卜04一lO 第5期 宁书年等：无线电波层析成像技术及在矿井坑透中的应用 卢(z，z)=∑∑口i怠(z，z)， 式中，口fi为模型系数；毫沩一组B样条；优和，z为模型网格z和z方向上的节点数． 这里卢(i，歹)=卢(i△z，歹△z) 的宽度，它们的大小按观测资料的多少、质量以及图像的分辨率的大小来选取．这样就可以用网格节点处 的卢(z，2)来描述介质模型． 由于实际电波射线并不一定总是经过成像网格节点，故需要用插值方法求p(z，z)及偏导数 格节点处的衰减系数解析关系为 卫 卢(z，z)=∑Kf展， (1)