关于相关辨识技术在谱激电中的应用研究.docxVIP

相关辨识技术在谱激电中的应用研究2014年度工作进展报告中国地质大学（北京）2014年12月扉页相关辨识技术在谱激电中的应用研究2014年度工作进展报告 项目负责人：李梅 报告编写：李梅 罗维斌中国地质大学（北京）2014年12月第一章重要研究进展1.1 课题的研究任务通过理论研究、数值模拟和实验验证等手段，研究不同介质的电场响应（冲激响应和频率响应）特征，总结相关辨识结果与地下介质电性分布的对应关系以及由辨识结果得到多种地电参数的方法，形成一个体系化的基于谱激电法的观测方案和数据处理流程，为方法的实用化提供理论和实践依据。1.2 研究技术路线课题提出了拟解决的关键科学问题：（1）不同地电系统的响应特征问题；（2）辨识结果与地电参数的关系问题。针对这两个问题，编制程序计算水平电流源层状地电模型的电场和垂直磁场的频率域响应特征及时间域阶跃响应特征。通过改变模型参数正演计算电场和磁场响应特征，分析研究电磁场响应与地层地电特征的对应关系。1.3 重要研究进展1编制了程序，计算水平电流源层状地电模型的频率域电场和磁场响应；2利用频率域电场和磁场值计算出全区视电阻率，利用全区视电阻率分析各场分量对目标层的分辨力；3 改变Cole-Cole模型参数，计算频率域电场和磁场响应，研究了激电效应对频率域电场和磁场响应的影响；4 计算了层状模型的阶跃响应；5 建立了基于观测电磁场时间序列的地电系统相关辨识算法；通过模型计算，基本认识如下：1、轴向装置电场响应对地层的分辨力优于赤道装置电场响应；2、赤道装置垂直磁场Hz得到的全区视电阻率比Ex电场得到的全区视电阻率更接近地层真电阻率；3、多道接收及道间变化率对地层有较高分辨力；4、阶跃响应晚期渐近值可指示地层电阻率变化；5、激发极化效应可引起明显的电磁响应。1.4存在问题及解决办法目前的研究仅限于一维层状介质，对于二度体电磁响应及对目标体的响应灵敏度分析还需要大量的计算工作。目前，先完善一维层状模型的响应分析计算。研究成果如下：第二章层状介质的电磁场分量频率响应特征水平电流源层状介质的电磁场响应已有精确的解析表达式（汤井田，2005；A.A. Kaufman，2001），其中涉及的Bessel积分，也称Hankel变换大多以数字滤波算法为主（Kerry Key，2012）。本课题中，采用了基于连分式逼近的数值积分算法（Walter L.Anderson,1989），计算稳定，效率和计算精度较高。2.1 水平电流源层状介质电磁场公式推导如图2-1所示，N层水平层状介质中第n层的电导率和层厚度分别记为和。水平电偶极子（接地导线）位于层状介质表面，偶极矩为P=IdL（I为谐变电流）。选取公共坐标原点位于偶极子中心的柱坐标系和直角坐标系，使x轴指向偶极矩方向（即y=0的方向），z轴垂直向下，求层状介质表面的电磁场分布。图2-1 层状地电模型略去繁冗的数学推导，此处直接写出极坐标系统中地表电磁场各分量的表达式为：（2.1a）（2.1b）（2.1c）（2.1d）（2.1e）式中：；，为自由空间导磁率。称为空间频率，它具有距离倒数的量纲。特别地，当N=1时，可得到均匀半空间表面电磁场各分量的表达式：(2.2a)(2.2b)(2.2c)(2.2d) (2.2e)以上各式中，为自由空间的导磁率。为谐变电流的圆频率。、和、分别为以为宗量的第一和第二类虚宗量贝寒尔函数，下标“0”或“1”表示贝寒尔函数的阶数，k为波数。必须指出，只有当观察点从上半空间趋近地表时，垂直分量Ez才不等于零。而当观察点从下半空间趋近地表时，Ez等于零。2.2 水平电流源层层状介质电磁场频率响应计算方法令：，式（2.1）可写为：（2.3a）（2.3b）（2.3c）（2.3d）（2.3e）式中(2.4a) (2.4b) (2.4c) (2.4d) (2.4e) (2.4f)(2.4g)(2.4h)(2.4i)(2.4j)(2.4k) (2.4l)在计算层状介质表面的电场分量和磁场垂直分量时，减去了均匀半空间的场，这是为了保证在时积分核快速趋于零，从而保证积分收敛，从而节省计算时间。写出(2.1c)，(2.1d)，(2.1e)中核函数为：(2.5a)(2.5b)(2.5c)比较发现上述三式有类似的特点，对其取极限，则有： (2.6a)(2.6b) (2.6c)这三个极限的收敛速度是很快的。据此，将(2.1)各式改写为(2.4)。当采用直角坐标时， (2.7a)