基于程函方程反演的跨孔地震岩溶探测数值模拟研究

岩溶探测对于岩溶发育地区的工程建设有重大意义，跨孔地震波探测是工程勘察有效的地球物理探测技术。本文旨在利用基于程函方程的初至旅行时层析方法，实现对跨孔地震波岩溶探测的反演，并通过多个工况的数值模拟验证其有效性，探索其适用性和局限性。首先，基于波动方程的时域有限差分对跨孔地震岩溶探测进行正演模拟；然后从模拟的地震记录中提取初至旅行时，最后以基于程函方程的初至层析方法反演得到钻孔间的地层速度模型。本文设置了单一溶洞，单一断层，溶洞断层组合及多溶洞组合4种工况，通过数值试验得到了地震记录、旅行时、反演图像、迭代残差等结果。结果表明基于程函方程的初至旅行时层析反演方法对溶洞位置、大小，断层的方位反演较为准确，对于复杂地质条件，多溶洞和裂隙组合反演较为稳定，相互之间干扰很小，能定性区分对溶洞填充物的速度与地层背景波速的高低。但岩溶异常体中的反演速度较真实速度有偏高，纵向分辨率比横向分辨率高。相比溶洞，该方法对断层较不敏感。考虑到基于程函方程的跨孔层析反演在多个岩溶地质体共存的复杂地质条件下依然保持高效和稳定，该方法可以为实际工程提供重要勘察信息。