我国煤矿水文地质条件复杂,导水通道(陷落柱、采空区、断层和裂隙等)种类繁多,导致煤矿水害事故频发。其中,作为煤矿中分布最为广泛的一种隐蔽致灾源,陷落柱具有隐蔽性强、突发性高以及突水量大等特点,对于煤矿生产以及人民的生命财产安全危害尤为突出,所以准确查明陷落柱的位置及其导水性,对于水害防治和保障煤矿的安全生产意义重大。目前,国内外陷落柱探测技术主要以地震勘探和电磁法为主。地震勘探尤其是三维地震勘探对地质构造有较高的探测精度,可以准确识别陷落柱的空间位置及形态,但该方法无法实现陷落柱富水性探测;而电磁法对低阻异常反映灵敏,可以有效探测陷落柱的富水性,但是该方法分辨率较低,无法准确圈定陷落柱的边界。也就是说,单独利用地震勘探或者电磁法中的任何一种方法都无法实现对陷落柱空间位置和富水性的精细探测,因此有必要通过有效的参数耦合方法来进行联合反演,克服单一方法的局限性,减少反演的多解性,实现陷落柱的精细定位及解释。论文首先采用交错网格有限差分法,对一阶速度-应力弹性波方程进行求解。以华北煤田主要地层结构为基础,建立陷落柱模型,实现了三维地震波场正演。研究和分析了不同直径和不同埋深陷落柱模型的地震波场响应特征,研究表明:①随着陷落柱直径减小、埋深增加,陷落柱顶底界面反射波振幅变弱;当陷落柱直径较小时,陷落柱的多次反射纵波在地震记录上呈现“烟圈”状。②对于本文研究的陷落柱模型,当埋深250m、直径与埋深比小于8%时,三维地震勘探方法仅能有效识别其顶界面,几乎无法识别出其底界面。其次采用有限单元法对麦克斯韦方程组进行求解,编程实现了2.5维和三维可控源音频大地电磁法的正演。分别与层状模型的解析解以及地下含低阻异常体模型的积分方程法计算结果进行对比,验证了本文程序的正确性和精度。对煤系地层中含有陷落柱模型进行电磁场数值模拟,着重分析了不同直径、不同埋深以及不同电阻率陷落柱模型的可控源音频大地电磁响应特征。研究表明,陷落柱的电磁场异常响应随其直径减小、埋深增加而迅速减弱可控源音频大地电磁法的分辨能力也随着陷落柱的直径减小和埋深增加而逐渐降低。研究了地震与可控源音频大地电磁的联合反演算法,阐述了人工蜂群算法的基本原理,包括人工蜂群算法的基本行为、算法流程以及主要参数等。针对该方法中的随机性过大和易陷入局部最优等问题,在人工蜂群算法的随机性、收益度公式以及新食物源的产生策略等方面进行改进,提出了改进的人工蜂群算法。利用标准测试函数,对改进后的人工蜂群算法的优化结果与粒子群算法、标准型人工蜂群算法进行了对比,评估了其优化能力。为了解决不同的地球物理方法中参数耦合问题,将交叉梯度函数加入到联合反演的目标函数中,并对交叉梯度函数的定义、离散和性质方面进行了研究。针对地震解释剖面中缺少地震波速或波阻抗等量化信息,导致无法计算交叉梯度的问题,利用聚类的思想,对传统交叉梯度法进行改进,提出了修正的交叉梯度法,并给出了算法流程。通过对二维速度和电阻率模型进行传统交叉梯度法与修正的交叉梯度法试算对比,验证了该方法的正确性和可行性。在地震波场正演和可控源音频大地电磁法正演的基础上,分别开展了陷落柱模型的可控源音频大地电磁法单独反演、基于精确地震解释的联合反演、基于“不精确”地震解释的联合反演以及基于“不精确”地震解释的多聚类联合反演,并对反演结果进行对比分析。结果表明,联合反演结果明显优于可控源音频大地电磁法单独反演,能够有效提高陷落柱的定位精度。最后,结合山东某煤矿疑似陷落柱的探测实例,分别对比了可控源音频大地电磁法二维反演、可控源音频大地电磁法2.5维反演和联合反演结果。结果表明,联合反演能够使陷落柱位置处的低阻异常与围岩呈现出明显的高低阻差异,提高了陷落柱的横向定位精度,并且在陷落柱形态方面能提供更多的细节信息,说明本文提出的联合反演算法具有较好的适用性。