瞬变电磁法S-反演和二维拟地震偏移成像研究.pdfVIP

瞬变电磁法S一反演与二维拟地震 偏移成像研究 龙凡 韩天成 (沈阳军区司令部工程科研设计院 沈阳 110013) 摘 要 瞬变电磁法是国防工程地质探测的新方法，由于不受复杂接地条件限制，加之其穿透能 力强、勘探深度大、体移I效应小、电性分辨率高而被广泛应用于勘探领域，现已成为军事工程地 球物理勘探的一种主要方法。但是．由于时域电磁理论的复杂性，一方面．高维的正、反演定量 解释方法目前正处于试验研究之中，另一方面，现有的解释软件又复杂繁琐、费时，很难适应野 外大批量数据实时处理成像的需要[1’“。我们利用分析地震的方法来分析瞬变电磁法的数据结果， 应用S一反演法和二维拟地震偏移成像原理编制了瞬变电磁反演成像软件，取得了较好的效果。 关键词 瞬变电磁法S一反演法 偏移成像 1 方法原理 1．1 S一反演法 S一反演法是一种基于薄层模型的TEM快速成像技术。水平薄板模型是瞬变电磁场正演 计算中唯一能用初等函数解析表示的地电模型，随着时间的推移，瞬变电磁场向地层深处传 播，因此，可将每个瞬问观测到的电磁信号等效为某一“浮动”导电薄层产生的场，从而直 接把观测的垂直磁场分量时问导数矩阵转换为电导率一深度值，计算方程如下： aHz(，。， do+t／ZoS 一3Pe·q Ot ‰一—丌,uo—r3S百干丌瓦干玩孺丁P 式中，，Pe为发射电矩，q为接收线圈的有效面积。采用恰当的方法解反问题，可由薄层 模型在每时刻所确定的两参数(深度cz和纵向电导5)得到电导随深度变化曲线。而纵向电导 又可由电导率随深度的积分公式表达： dp 5(d)一j嗽)如 (2) 因此，通过对计算曲线沿深度微分，便可得到随深度变化的电导率。若对每个记录点均 作纵向电导反演，对于连续的观测时间．即可形成每条剖面连续的电导率一深度图像。 1．2 电磁拟波动方程偏移成像 I．2．1定解问题方程 理论研究表明，满足扩散方程的低频电磁场与满足波动方程的波场间存在如下积分变换 关系： 一55— 型《掣一—士ihⅥ)_，辜．『，㈡、h∽．(1，／ ㈩‘川 —_厂一一瓦■鬲儿-¨7’_…一卜～、厂“7卜“‘7 通过选择合适的数值算法及参数．可?陌瞬变电磁场随时问的衰减川1线变换为类似．-j二地震 波的波传：潘曲线U(hq)．其传播的慢度为＼／n，d。 与瞬变电磁场变换波场相对应的定解问题方程为： ∥己∥叫“盯嘉L7、邓”嘉‘， ㈩ 弓f|一 aI}‘ U’(r，：，0)一0 (，5lr．一(，≯(_．0．吖) U’Ir。一0一一、厅下≯一cV， 直接利用地震勘探中有限差分偏移算法．可对这种拟波动场进行偏移成像处理．以确定 地电断面的几何参数。 1．2．2二维有限差分逆时偏移算法 为推导方便起见，在这里给出声波方程的一般形式： (’2(拍2)·[等+警]+厂c丁，：，f)一等LaI，r。 az‘J ar‘ cs， 数。 对比(4)与(5)即可得到以上各函数的具体形式。 采用有限差分格式来模拟波的逆向传播过程，求得波动方程的数值解。 对于函数U(．z)．它的Taylor展开式可以写成： 这样，就可以得到：