三维大地电磁场的图像性质的研究.pdfVIP

三维大地电磁场的图像性质研究 张大海1～L．Szarka2 (1中国科学院上海天文台上海200030；2匈牙利大地测量与地球物理研究所匈牙利) 1．研究意义 本文主要是对复杂地下介质情况下3维大地电磁场的数值模拟结果进行了相关性分析和图 像性质研究，研究结果对于多维大地电磁问题的数据处理、反演的目标函数选择以及资料解释都 有很重要的意义。 2．方法 在已有的旋转不变量中．例如平方和不变量ssq、主元差不变量z1等多种形式，我们选择了模 值不变量1引作为例子来说明这个问题，其他形式的旋转不变量具有相似的情况。在数值模拟过 程中，我们应用了2种3维大地电磁正演模拟程序，一个是俄罗斯Spichak教授的程序，另一个是 and 美国MackieMadden的程序，这2个程序都使用了有限差分法作为原理。在数值模拟过程里， 我们不仅使用了计算机合成的模型(近地表u形薄板异常体)，而且我们也使用了野外的多个测 井数据形成的盆地模型(位于匈牙利中部的／3ekes盆地)；同时我们又在野外实际模型的基础添加 了一个电导基底，增加了模型的复杂程度。 针对基于不同大地电磁阻抗张量旋转不变量定义的视电阻率响应和相位响应的图像特性进 行了分析．同时结合不同响应值和地下复杂介质中的电性分界面的相关性研究，我们分别使用Re 2000)，其中det代表模值，Re代表实部，Im代表虚部，z代表阻抗张量矩阵。 3．分析和结果 我们发现对于所有模型结果，从相关曲线图的角度来讲由Re_det(z)定义的响应值与地下介 质电性分界面的相关性与负相位的相关性曲线最为接近，几乎在所有频段上接近重合，detRe(z) 的效果次之；同时我们引用了Baskour的定义进行比较，结果显示他的定义也获得了较好的结果； 在模型Bekes盆地+电导基底中，随着基底的深度加深，从相关曲线上可以看出基底对浅部异常 体的影响在减弱，同时相位曲线趋于0的频率值在降低；对于近地表高、低阻u型薄板模型，我们 同样观察到相关系数曲线在高频段存在一个极值的突变过程，而在图像性质分析中在相同频段存 在着一个图像强烈对比现象，我们称这个为”关键孔”现象；针对以上现象我们集中在关键7L附近 频段加密了采集数据，对得到的数据同时进行相关性分析和图像研究；研究结果显示，无论高阻异 常体还是低阻异常体当相关系数在一个很窄的周期内(t0．02s)由一个极值跳变到另一个极值的 时候，从图像上我们就可以看到视电阻率和相位值发生一个灰度显示的一个逆转，我们可以发现 在这2个极值中间的某一点的相关系数接近0；当我们对比高阻异常体和电导异常体的相关系数 曲线，正如我们期望的那样，从时间角度看电导异常体曲线比电阻异常体曲线要想后推一段距离， 也就是说它们需要更大的时间间隔；对于另外一个模型(U型薄板异常体+阶梯状导电基底)，相 关曲线图同时也显示了实部基的旋转不变量定义的响应较为符合负相位的变化趋势，而且对电性 的变化具有一个较敏感的特性，从全频段视电阻率和相位的图像中也可以看出Re—det(Z)基的视 电阻率和相位图像较为真实的反映了地下的电性分布情况。 总之，通过大量合成模型和野外真实模型的数值模拟试验，我们认为通过对基于不同阻抗张 ·579· 量旋转不变量定义的大地电磁响应和地下电性地质结构的相关性分析以及图像分析．相关性曲线 在～定程度上与地下电性结构的分布是有其本质的联系的．甚至在近地表的3维异常体和3维深 部构造情况下对真实的地下电性分布的认识是有意义的t实部基定义的视电阻率较其他定义更接 近相位相关系数的曲线变化，结合相关系数和图像分析，我们可以确定”关键孔”的位置，从而为分 辨地下异常体几何形状提供更丰富的信息；慎重的选择目标函数，提高反演质量。当然本研究也 存在着一些问题。希望能够深入研究”关键孔”现象形成的物理机制，以及对噪声的影响等问题。 本研究由中科院知识创新I程重要方向性项目(No．KJCX2一sw—T1)和匈牙利自然科学研 究基金(No．TS40848)资助。 参考文献 [1]IkHskour九’F．，GeophysicalProspecting，42．1，1994． 1．．eta1．AclaG舯d [2]Szarka Geoph,Hung．35，1，2000．