井间地震层析成像方法的研究及应用.pdfVIP

井间地震层析成像方法研究及其应用 裴正林¨ 余钦范∞ 石油大学，北京，100083；2)中国地质大学，北京，100083 内謇提要 目前．井间地震层析成像技术已发展到宴用化阶段。然而，我国层析成像技术 在地质灾害方面的应用属刚雕起步，尚无成热的方祛技术。本文以实际工程的应用研究为倒来展 示我们几年来在发展这一技术及应用方面的进展和研究水平。 美t词 地震层析成像左岿窀f摩成席 0 币 、 j t ～ 井间地震波、声波的层析成像技术．能够提供被探测目的物的精细结构和岩性变化的直 观图像，因此，得到国内外地质学家和地球物理学家的重视。目前，井间地震层析成像技术已 发展到实用化阶段。然而，我国层析成像技术在地质灾害方面的应用属刚刚起步，尚无成熟 的方法技术。 现行的井间地震层析成像方法对初始模型依桢性强，仅适合于简单模型和低对比度(速度 差小)异常体成像，图像分辨率偏低。而地质灾害体如岩溶、陷落拄、裂隙、裂缝、断层破碎带、软 弱夹层等井间地震层析成像属于复杂模型或大扰动异常体的非线性成像，显然，对于地质灾 害体，采用现行的井间地震层析成像方法无法获得好的成像效果。引起图像分辨率偏低、成像 效果不好的的主要原因是：①有限观测角；②介质的非均匀性和各向异性；③观测数据含有 误差；④成像方法本身的缺陷(如线性化、正则化和弱散射近似等)。因此，井间地震CT仍需 从反演理论及方法、反演算法，野外施工与数据采集等方面进行研究、改进，进而很好地解决上 述几方面的问题，才能提高地震CT技术在地质灾害勘查中的应用效果。本文以实际工程的应 用研究为倒来展示我们几年来在发展这一技术及应用方面的进展和研究水平。 1井间地震层析成像方法 地震层析成像问题的非线性是与介质的不均匀性联系在一起的。地震层析成像反演的困 难主要是目标函数中存在许多局部极小．因而影响了迭代反演达到全局极小。为了有效的解 决多极值日标函数的优化问题，即减小非线性反问题的复杂性，我们将小波多尺度分解的思想 引入其中，提出了小波多尺度井间地震层析成像方法(裴正林等，1999a)。该方法能明显改善成 像效果，提高了图像分辨率，而且有较强的抗噪能力，适合于大扰动非均匀模型和复杂模型成 像。 小波多尺度井间地震层析成像方法是一种多尺度逐次逼近的线性化迭代方法，本质上仍 属于确定性反演方法，因此，在各尺度上的迭代反演仍存在搜索陷入局部极小的可能性，自然 无法彻底地解决井间地震层析成像的非线性问题。与确定性反演方法相比，随机擐优化方法无 注；本文为“九五”国家科技项目(编号95—913—06—02—03)的部分研究成果 731 须对非线性目标函数做任何近似和求导数。这类最优化算法仅要求模型的高精度正演计算和 目标函数(或误差函数)值计算。因此，随机最优化方法是一种完全的非线性反演方法，其结果 一般优于确定性算法的结果，即更能真实反映实际模型。 井闯地震层析成像是高维的非线性反问题，其可行解空间巨大，因此，从计算成本而言，一 拟自然界进化中优胜劣汰的优化过程，原则上能以较大的概率找到全局的最优解，具有鲁棒性 强，全局收敛性好等优点，在模型搜索时具有继承性和并行性，因此优于SA法，而简单的遗传 进化算法存在收敛速度慢和编码问题。为此，本文提出了一种适合于井间地震层析成像的遗传 进化反演方法(裴正林等，2000aj。为了提高反演速度，我们将小波多尺度分析(MRA)思想引 人到随机最优化过程中，即将多尺度反演和进化算法反演的优点结台起来，形成多尺度逼近遗 传进化算法(裴正林，2000b)。 2并间地震层析成像的分辨率与有效观测系统 我们利用数值模拟试验，讨论了井间地震弯曲射线走时成像分辨能力、在不降低分辨率条 件下最佳观测系统等问题(裴正林，2000b)。 2．1井问地震屡析成像的分辨率 井间地震层析成像反演算法的解的物理意义是绐出地质结构和岩性分布，因此，对于解的 可靠性和分辨率的研究就显得非常重要