≯拟地震法电磁波场有限元法一了正演.pdfVIP

第柏眷晕6柑 十由稿；l事谯拳挂 Vo1．23 6 1992年 12月 I．CENT．一SOU 。INST．MIN．METALL D ．1992 ≯一 了 拟 地 震 法 电磁 波场有 限元法 正 演 模 型 塑壁董／鲍光 淑 ． (地 质 系 ) 摘 要 在 电法勘探正演模拟 中，利用有限元法提 出7电磁波在二维 电阻率介质模 型 中的传播特征 ．我们将 电磁波波动方程与声波波动方程熹 比，运用适 当的近似 ，取 等效方程 中相同的主要项，利用适当的边界条件，则电磁波场 问题可化为类似于地震 波场 问题，从而可用地震资料处理技术中的任何方法去完成．我们连用的是波动方程 有限元方法．计算可在频域也可在时域进行去模拟=维 电阻率分布的响应即电场视 剖 面．实例表 明：正滴结果所得 到的二维 电阻率成像剖面与实际的二维 电阻率分布 关键词有限元法，电扬视剖詈面}地震学 )硝么予． J饧匀 在电法勘探正演模拟中，有限元法是用来求解稳定场的变分问题的数值方法．用这种方 法可求解复杂构造中稳定电流场的分布．Cog~on (1971)首次根据电场总能量最小原理， 导出了电位的变分方程，并求 出了稳定场的分布．而利用有限元法来仿地震波场的波动特征 至夸来见报道 ．作者基手伽辽金法建立了电磁波有限元方程，构建了与地震剖面相类似的正 演 电磁波视剖面．可以看出，利用有限元涪不但可 以对=维 电阻率模型进行模拟和成像，而 且可对更为复杂的 (亦即更切台实际的)模型模拟和成像，更重要的是为电磁波场新的处理 技柬应用迭加、偏移等提供了可能伍和前提． 电磁渡拟地震珐解释在国内 不论是在理论 上还是实践上还刚刚起步．一方面缺乏新理 论指导，另一方面在实陆资料的应用中还束迈出步子，为此，我们从新的理论——波动方程 理论为指导来研究电磁波拟地震波波动方程的正演模拟．在下面的内容中，我们首先讨论声 波波动方程与电磁波波动方程的相似性，并通过适当的近似，利用伽辽金法推出电磁波有限 元方程．最后用实例给出相应的电场视剖面，证实刹用有限元法进行EM数据成像想法的正 确佳． 1 基本理论 在额率域中，电磁波场 中的电场 与磁场 的关系由Maxwell方程得如下方程， 1991年 6月 2 日收到初稿 ’国家 自拣科学基叠赞助课艇 十钉曲堙童伍举|t 1口口2年1±月 l∞H =一(1／)V ×￡ (口+ic．oe)E=V ×H ( ：[1／((三=川]VxH) ‘ 其中 为导磁率， 为介质的电导率，8为介电常数．方程 (1)在二维傅里叶域 (∞， )中 E极化部分的结果为。 未[：]=[。0， 一 ][：] cz， (2)武 中 [oGco)一+。]-I．H 极化部分亦可类似表示． 在二维无源地震波场中，声压强 与粒子速度向量 = (， )在频率域中的关系由 下式给出 (Clearboutt”，1976)： f∞ =一(1／p)v (3) 、 icop= ～K V · 其中，∞