论文部分内容阅读
随着社会的不断发展,人们日益注重与个人生命财产安全息息相关的工程建设安全、生态环境安全,以及地质灾害防护等问题。该类问题的解决离不开浅层地球物理勘探方法。地震CT技术与直流电阻率法作为两个重要的浅层地球物理勘探手段,被广泛的应用于工程勘察等领域。该两种方法在常规浅部地层区域内勘探效果明显,但当地下存在低速地质体时,受射线稀疏影响,地震成像效果不佳。当地下存在高阻体或中间薄层时,电场线分布较少,直流电阻率勘探精度不足。为了弥补该两种方法的劣势,前人在矩形网格剖分的基础上,实现了地震初至走时数据与直流电阻率数据的交叉梯度联合反演。该方法在一定程度上实现了不同物性方法之间的优势互补,提高了二者的反演精度,但矩形网格剖分难以模拟起伏地形和复杂模型。为此,本文开发了一种基于非规则三角网格剖分的井间地震CT与直流电阻率的交叉梯度联合反演算法。该算法与前人算法的主要不同之处在于:(一)采用了非规则三角网格剖分;(二)井间地震初至波走时射线追踪方法采用的是基于非规则三角网格剖分的井间地震初至波走时线性插值法;(三)在相邻三角单元物性参数分布线性假设条件下,实现了交叉梯度函数的非差分形式离散。此外,开展了井间地震CT与直流电阻率法的单独反演,以便于分析联合反演效果。为提高非结构化三角网格井间地震CT初至时间计算精度,正演使用本文开发的新算法,反演则是采用阻尼最小二乘反演方法;直流电阻率法正演使用有限单元法,反演则采用基于最小梯度支撑的直流电阻率正则化聚焦反演算法。通过对理论模型的计算及对比,得出了以下三点结论:(一)基于非规则三角网格剖分的井间地震初至波走时线性插值法,能够对复杂的初至波走时模型进行数值模拟;(二)以相邻三角单元物性参数线性假设为基础的交叉梯度函数求解方法,是可行的;(三)基于非规三角网格剖分的井间地震CT技术与直流电阻率法交叉梯度联合反演,使得不同物性模型之间相互约束,但获得的反演结果仍有待进一步提高。