地球物理学报  2018, Vol. 61 Issue (6): 2548-2559   PDF    
大地电磁资料精细处理和二维反演解释技术研究(六)——交错模型的大地电磁二维反演
郭春玲, 陈小斌     
中国地震局地质研究所地震动力学国家重点实验室, 北京 100029
摘要:交错模型指的是由几组走向不同的线性构造(近似二维)在空间(纵向、横向)上组合而成的模型.对于电性结构,交错模型表现为纵向或横向上的构造存在电性主轴方位的变化.从构造维性的角度来看,交错模型是一种由二维模型组合而成的特殊的三维模型,直接对其进行大地电磁二维反演,不易获得可靠的反演结果.本文针对交错模型的特点,提出分频段-分区段反演方案.该方案首先需要借助于阻抗张量成像技术,在频率域确定组成交错模型的各线性构造的电性主轴,然后,针对不同的频段、区段,选择对应电性主轴的数据进行反演,通过初始模型的构建将不同电性主轴方位的反演结果对接起来.本文通过三维理论模型的研究,系统展示了分频段-分区段反演的全过程,归纳得到:在分频段反演时先做低频段反演,在分区段反演时先做二维性更不显著一侧的反演.最后,本文将这一技术用于郯庐断裂带中南段一条实测剖面的反演中,其结果与常规二维反演结果相比较,深部的信息更为丰富,且与其他已有地质、地球物理结果的可对比性更好,表明在构造复杂地区,大地电磁分频段-分区段二维反演具有较高的模型分辨率和可靠性.
关键词: 交错模型      大地电磁      二维反演      电性主轴      郯庐断裂带     
Refined processing and two-dimensional inversion of Magnetotelluric (MT) data(Ⅵ): Two-dimensional magnetotelluric inversion based on the staggered model
GUO ChunLing, CHEN XiaoBin     
State Key Laboratory of Earthquake Dynamics, Institute of Geology, China Earthquake Administration, Beijing 100029, China
Abstract: The staggered model refers to a combination of several groups of linear structures (approximately two-dimensional) in space (vertical and horizontal directions). For electrical structures, the staggered model is characterized by variations of the geo-electrical strike in longitudinal or transverse directions. From the viewpoint of structural dimension, the staggered model is a special three-dimensional model made up of two-dimensional models, which is not easy to obtain reliable inversion results by two-dimensional Magnetotelluric (MT) inversion directly. Considering such characteristics, a sub-band and sub-segment inversion scheme is proposed in this work, in which the geo-electrical strikes of the linear structures are determined in the frequency domain by means of impedance tensor imaging. Then, the data of the corresponding strikes are selected for different frequency bands and segments, and the inversion results of different strikes are connected to each other through the construction of the initial models. In this paper, the whole process of sub-band and sub-segment inversion is demonstrated by a theoretical three-dimensional model. The results show that the inversion for low frequency band data should be performed at first in sub-band inversion, and in sub-segment inversion the segments with weak two-dimensionality are inverted in advance. This technique is applied to the inversion of a measured profile in the middle and southern sections of the Tanlu fault zone. The result is compared with the conventional two-dimensional inversion results, which shows that the new method permits to reveal more abundant information at depth, consistent with existing geological and geophysical results from previous work. It is shown that the two-dimensional MT inversion of sub-band and sub-segment has higher resolution and reliability in areas with complicated structures.
Key words: Staggerd model    Magnetotelluric    2D inversion    Geo-electrical strike    Tanlu fault    
0 引言

地球在漫长的构造演化过程中,会经历多期次的应力作用,导致岩石圈的纵向和横向结构具有不同的几何特征.张国伟等(1995, 1997)研究发现,秦岭造山带岩石圈深部结构以近南北向为主,而浅部则以近东西向为主.在北非阿特拉斯山脉,也发现地壳和上地幔构造走向的不同(Kiyan et al., 2010; Anahnah et al., 2011Ledo et al., 2011).在本研究系列(五)中,我们利用大地电磁阻抗张量成像技术发现了合肥盆地的深部结构中存在多期次的立交桥组合构造(陈小斌和郭春玲,2017).此外,横向上也经常存在不同走向的构造之间的交叉接触关系,一条剖面会同时穿越多个具有不同走向的构造体系.更多的情况下,存在纵向立交桥交错模型和横向交错接触模型之间的组合.为论述简单,将此类构造模型(纵向立交桥交错模型、横向交错接触模型及它们的组合)统称为“交错模型”.

交错模型在本质上属于三维模型,不满足二维模型所假设的整条剖面具有相同的构造走向的要求,因而导致大地电磁二维反演中电性主轴方位选择的困难.北非阿特拉斯山脉的MT探测导致Jones课题组和Anahnah课题组之间出现了激烈的争论(Anahnah et al., 2012Jones et al., 2012),因为二者在反演中选择了不同的电性主轴,反演结果出现了巨大的差异.但双方最终并没有达成共识——双方均未为阿特拉斯山脉壳幔电性结构这一交错模型的二维反演提出合理的解决方案.

然而,交错模型在空间上可以分解为一系列具有不同走向方位的二维线性构造,这意味着交错模型不是一般意义上的没有明确构造走向的三维模型.能否利用交错模型的这一特征,寻找一种比固定主轴方式更为合适的大地电磁二维反演技术呢?

在传统的二维反演中,遇到横向交错模型时,会采用分段反演,每一段分别有对应的主轴方位,各段之间会有些测点的重叠.然而,这种反演方式会造成各段之间对接的困难.即在重叠段,反演得到的电性结构多数情况下会有所差异.此外,这种传统的分区段反演,各区段主轴方位的确定不是来自大地电磁的信息,而主要依据地表地质调研所给出的构造分布.对于纵向交错模型,迄今并没有好的办法.

本文针对上述问题,利用大地电磁二维正演(陈小斌等,2000)、三维正演(Mackie et al., 1993)、阻抗张量统计成像(陈小斌等,2014)及二维非线性共轭梯度法反演(Rodi and Mackie, 2001)等技术,开展理论交错模型和实测数据的反演对比研究,力求为解决交错模型的大地电磁二维反演提供合理的新技术.

1 理论交错模型的二维反演 1.1 模型一:纵向交错模型(立交桥模型)的二维反演 1.1.1 模型与正演计算

本文研究过程中设计了多个理论模型来研究交错模型的大地电磁二维反演问题.这里选择了其中的一个进行讨论.该模型中,浅部为走向呈45°的高阻体(电阻率10000 Ωm),深部为走向0°的低阻异常体(电阻率为1 Ωm),如图 1所示.很显然,这是纵向交错模型中比较简单的、典型的三维立交桥模型.我们对该三维模型及其目标剖面所对应的二维模型分别进行了三维正演和二维正演计算,获得了相应的正演响应数据,然后,采用MT-pioneer软件系统(陈小斌等,2004)进行后续的阻抗张量成像分析和二维反演的对比研究.

图 1 (a) 纵向交错模型(立交桥模型); (b)该模型的电性结构的平面投影(图中虚线L为进行反演对比分析的目标剖面位置); (c)测线L剖面所对应的二维模型 对(a)所示的三维模型,利用Mackie等(1993)研发的程序进行了三维正演计算,计算网格为86×86×76,频段范围为0.0001~1000 Hz,共40个频点;对(c)所示二维模型采用陈小斌等(2000)研发的程序进行了二维正演计算,计算网格为86×76,频率参数与三维正演完全相同. Fig. 1 (a)The longitudinal staggered model (overpass model); (b)the plan view of the electrical structure of the model (the dash line L is the target profile for inversion contrast analysis); (c)The corresponding two- dimensional model For the three-dimensional model shownin(a), three-dimensional forward modeling was carried out using the code by Mackie et al. (1993) with a grid of 86×86×76 and the frequency range of 0.0001~1000 Hz, and 40 frequency points in total. For the two-dimensional model in (c), two- dimensional forward modeling was carried out using the program by Chen et al. (2000) with a grid of 86×76 and the same frequency parameters as in 3D forward modeling.
1.1.2 阻抗张量成像分析

阻抗张量成像技术(陈小斌等,2014)为交错模型的二维反演提供了频率和空间上的完整的电性主轴分布信息,从而为探寻合理的交错模型二维反演方案提供了技术基础.图 2给出了立交桥模型中测线L上所有测点和频点的电性主轴的统计成像结果.由图可知,测区存在两个主要的电性主轴方位,分别为45°和0°,其中45°主轴方位分布在高频段(频率>0.3 Hz)和剖面中心位置附近,0°主轴方位分布在低频段(频率 < 0.3 Hz)的整个L剖面上.因此,利用阻抗张量成像技术,我们确定了L剖面上深浅不同走向构造所对应的频段:高频段(0.3~1000 Hz)对应于浅部45°走向的构造,而低频段(0.0001~0.3 Hz)对应于0°走向的构造.利用这些信息,我们就有可能采取比传统更合理的二维反演方案.

图 2 立交桥模型电性主轴统计成像结果(蔡军涛-陈小斌-赵国泽法,即CCZ方法,蔡军涛等,2010) (a)统计玫瑰图;(b)频率分布云图(陈小斌等,2014);(c)测点分布云图(陈小斌等,2014).图中色标表示经二维性加权后的统计频次(下同). Fig. 2 Statistical images of geo-electrical strikes for the overpass model using the CCZ method (Cai et al., 2010) (a) Rose diagram; (b) Frequency cloud diagram; (c) Site cloud diagram. In the cloud diagrams, the color scale represents the statistical frequency with dimensional weighting (the same below).
1.1.3 二维反演结果比较

(1) 常规二维反演

对于一条实测的大地电磁剖面,常规的二维反演一般要求全剖面处于同一电性主轴方位下.对于图 1所示的模型而言,存在两个构造主轴方位:0°和45°,那么在常规的二维反演中,可能存在以下几种处理方式:整条剖面都采用0°走向反演、45°走向反演、二者的平均值22.5°走向反演.为了能够全面地对比常规反演与可能的新反演技术,本文分别对45°、0°、22.5°走向三种情形进行了二维反演,反演中初始模型均为100 Ωm的均匀半空间,反演结果见图 3abc.由于TE极化模式对于二维性的要求要远远强于TM模式(蔡军涛等,2010),TE模式数据参与反演,会使得反演结果带有更大的不确定性,影响对比分析,故以下反演中我们均只采用了TM极化模式进行二维反演.

图 3 立交桥模型二维反演结果对比 (a)、(b)和(c)依次为取主轴方位为45°、0°、22.5°的反演结果;(d)为先做高频数据反演(0.3~1000 Hz,主轴方位45°),以其结果为初始模型再做低频数据反演(0.0001~0.3 Hz,主轴方位0°)所得到的结果;(e)为先做低频数据反演(0.0001~0.3 Hz,主轴方位0°),以其结果为初始模型再做高频数据反演(0.3~1000 Hz,主轴方位45°)所得到的结果;(f)为L剖面所对应的二维模型的二维反演结果,可作为比较标准.采用NLCG二维反演算法程序(Rodi and Mackie, 2001),反演网格为91×173.反演最原始的初始模型均为100 Ωm的均匀半空间模型.每幅图上,‘tau’为正则化因子,‘floor Err(%)=TM 2’表示采用TM极化模式数据,反演门限误差为2%,‘Iter’表示迭代次数,‘RMS’为拟合均方差(下同). Fig. 3 Comparison of two-dimension inversion results of the overpass model (a), (b) and (c) are for geo-electric strikes 45°, 0° and 22.5°, respectively; (d) The inversion result for low-frequency band data (0.0001~0.3 Hz, strike azimuth 45°) with the initial model constructed by the inversion result of high-frequency band data (0.3~1000 Hz, strike azimuth 0°); (e) The inversion result of high-frequency band data(0.3~1000 Hz, strike azimuth 0°) with the initial model constructed by the inversion result of low-frequency band data (0.0001~0.3Hz, strike azimuth 45°); (f) The two-dimension inversion result of the corresponding two dimension model of L profile which is used as a comparative standard. All of the inversion models use the 2D inversion algorithm of NLCG(Rodi and Mackie, 2001), in which the grid is 91×173, and the initial model is homogeneous half space with 100Ωm resistivity at the beginning. In each figure, 'tau' means the regularization factor, 'floor Err (%)=TM 2' is the TM polarization mode data are adopted with the data floor error is 2%, 'Iter' means the number of iterations, and RMS is fitting mean square error (the same below).

(2) 分频段二维反演

我们已经通过阻抗张量成像技术得到了图 1所示立交桥模型的深浅构造电性主轴所对应的频率范围,即高频段(0.3~1000 Hz)对应于45°走向的构造,低频段(0.0001~0.3 Hz)对应于0°走向的构造.根据这些主轴方位的频率分布信息,我们很自然地想到,可以先做某一频段的二维反演,然后以其为初始模型做另一频段的二维反演.这种方案下,由于每次反演所面对的数据都是同一主轴的数据,在一定程度上可以实现各线性构造体之间的解耦,从而其二维反演结果有可能更为可靠.

对于图 1这样的较为简单的立交桥模型而言,可以有两种分解反演的方式:其一、先做高频反演,然后以其为初始模型进行低频反演;其二、先做低频段反演,然后以其为初始模型进行高频数据的反演.每一种反演过程中,网格是完全相同的,且所有模型参数都将得到修正.以这两种方式得到的二维反演结果见图 3d图 3e,其中最原始的初始模型为100 Ωm的均匀半空间.作为一个比较的标准,图 3f图 1中测线L剖面所对应的二维模型响应的二维反演结果.

根据图 3中二维反演结果的比较,我们可以得出:

(1) 45°主轴常规反演结果(图 3a),浅部高阻异常体的结构形态与图 3f几乎吻合,深部低阻异常体顶、底边界与图 3f反演结果相同,但左右边界引入了冗余构造.

(2) 0°主轴常规反演结果(图 3b),浅部高阻体结构形态扭曲变形,下方出现低阻体假构造,且低阻块体包围在高阻异常体周围,与图 3f形态差异较大.

(3) 22.5°主轴常规反演结果(图 3c),整体能较好地反映真实的模型信息,但浅部高阻体底部深度不够准确,下方出现假高导体,深部低阻体上边界左右两侧向高阻体弯曲严重.

(4) 高频反演结果作为低频数据反演的初始模型,反演结果见图 3d,浅部高阻体底边界引入了高导环假构造,深部低阻体上边界与图 3f相一致,但体积差异大.

(5) 低频反演结果作为高频数据反演的初始模型,反演结果见图 3e,深、浅异常体与实际模型吻合良好.

因此,根据上述比较结果,可以得出一个结论:先反演低频段,以其为初始模型再反演高频段数据,得到的反演结果最为可靠.这可能是因为低频数据反演不仅能反映深部电性结构,而且对浅部结构也有一定反映,以低频反演结果为初始模型进行高频数据反演,在保证低频反演结果可靠的基础上,对高频数据作出进一步的改善.

在常规反演结果中,22.5°主轴方位的反演结果也比较可靠,这主要是因为22.5°与0°和45°的差别都不是很大的原因.陈小斌等(2008)曾经对大地电磁二维反演中的数据旋转方向的选择做过讨论,认为在旋转角小于30°时,反演结果的差别不是特别大,是可以用于解释的.

1.2 模型二:横向交错模型的反演 1.2.1 模型与正演计算

图 4为本文设计的另一种交错模型.该模型的浅部有一个0°走向的低阻体、深部存在45°走向的高阻体,选择的剖面位置在两个异常体交错点的一侧,使得该剖面上存在两个构造走向不一样的测段,形成横向交错模型.对该模型进行了三维正演和L剖面所对应的二维模型进行了二维正演计算.

图 4 (a) 横向交错模型;(b)该模型结构的平面投影(图中虚线L为进行反演对比分析的目标剖面位置,其上标出了参与反演的测点名称);(c) L剖面所对应的二维模型 对(a)图三维模型利用Mackie等(1993)研发的程序进行了三维正演,计算网格为88×98×64,频段范围为0.0001~1000 Hz,共40个频点;对(c)图所示二维模型采用陈小斌等(2000)研发的程序进行了二维正演,计算网格为88×64,频率参数与三维正演完全相同. Fig. 4 Same as Fig. 1 but for the horizontal staggered model
1.2.2 阻抗张量成像分析

图 5给出了横向交错模型中测线L上所有测点和频点的电性主轴的统计成像结果.统计玫瑰图表明,该剖面上存在两个显著的主轴方位:0°和45°,其中0°特别突出,在传统的二维反演中,极可能被作为全剖面的主轴方位;由频率分布云图可知,0°主轴方位几乎占据了整个频率范围,45°主轴方位则主要分布于中低频,亦即不同的电性主轴方位分布在频率域存在重叠现象,表明不能在频率域将不同走向的线性构造区分出来;而由测点分布云图可知,大致以剖面中部的47号测点为界,0°主轴方位分布在剖面的右侧,45°主轴方位则分布在剖面的左侧.统计云图的这些特征与模型剖面的实际电性结构是相吻合的.

图 5 横向交错模型电性主轴统计成像结果(CCZ方法,蔡军涛等,2010) (a)统计玫瑰图;(b)频率分布云图;(c)测点分布云图. Fig. 5 Same as Fig. 2 but for the horizontal staggered model

因此,依据阻抗张量成像的分析结果,确定了横向交错模型中不同走向构造所对应的区段:以测点47为界,左侧(测点编号≤47)对应于45°走向的构造,而右侧(测点编号>47)对应于0°走向的构造.与分频段反演类似,下一步采用分区段的方式进行二维反演.

1.2.3 二维反演结果比较

(1) 常规二维反演

与纵向交错模型(模型一,立交桥模型)的处理方式相同,将测线L上参与反演的所有测点的所有频点数据分别在0°、45°、22.5°主轴方位角下进行二维NLGG反演,反演中初始模型均为100 Ωm的均匀半空间,反演结果见图 6abc.

图 6 横向交错模型的二维反演结果对比 (a)、(b)和(c)依次为取主轴方位为0°、45°、22.5°的反演结果;(d)先做右半剖面数据反演(测点编号>47,主轴方位0°),以其结果为初始模型再做左半剖面反演(测点编号≤47,主轴方位45°)所得到的结果;(e)为先做左半剖面反演(测点编号≤47,主轴方位45°),以其结果为初始模型再做右半剖面数据反演(测点编号>47,主轴方位0°)所得到的结果;(f)为对应的二维模型的二维反演结果,可作为比较标准.采用NLCG二维反演算法程序(Rodi and Mackie, 2001),反演网格为93×153.反演最原始的初始模型均为100 Ωm的均匀半空间模型. Fig. 6 Same as Fig. 3 but for the horizontal staggered model

(2) 分区段二维反演

根据上述阻抗张量成像分析结果,得到了左、右两个不同走向方位的区段,因此可以采用分区段的方式进行反演.与分频段反演类似,存在两种可能的分区段反演方式:方案一(先右后左方案),先对右半剖面(测点编号>47)的数据以0°主轴方位进行反演,通过将误差增大为100万倍使得左半剖面的测点数据在反演中不起作用,然后以其结果为初始模型,在45°主轴方位反演左半剖面的(测点编号≤47)数据,此时右半剖面的数据误差被增大为100万倍使其在反演中不起作用;方案二(先左后右方案),先对左半剖面(测点编号≤47)的数据以45°主轴方位进行反演,然后以其结果为初始模型,在0°主轴方位反演右半剖面(测点编号>47)的数据,误差处理方式同方案一.按照这两种方式进行NLCG二维反演的结果见图 6d图 6e.

根据图 6中二维反演结果的比较,我们可以发现:

(1) 0°主轴方位的常规反演结果(图 6a),低阻异常体与图 6f吻合较好,但高阻异常体分辨率低,体积小,且顶、底边界较实际模型差异较大.

(2) 45°主轴方位的反演结果(图 6b),高阻异常体分辨率高,较好地反映了真实的模型信息,但与图 6f相比较,体积更大,电阻率更高;浅部低阻体确实也反演出来了,但在其下方出现了显著的高阻假构造,其电阻率值甚至比剖面左边的高阻异常体还要大,该假高阻体将严重干扰反演结果的解释.

(3) 22.5°主轴方位的反演结果模型(6c)较好地反映了高、低阻体的位置特征,但与图 6f相比较,高阻体有点过于突出而低阻体得到了一定程度的压抑;此外,在低阻体下方,也出现了假高阻异常体.

(4) 分区段反演方案一(先右后左)的结果(图 6d)中,高、低阻体的电性结构都反映地较好,但相较于图 6f,高阻体的分辨率稍弱,且其上、下边界均偏浅.

(5) 分区段反演方案二(先左后右)的结果(图 6e)是所有横向交错模型的二维反演中最接近于图 6f的模型.

上述结果对比表明,采用分区段二维反演,比常规的单一主轴方位的二维反演结果要更加可靠,其出现冗余假构造的可能性更小一些.此外,我们还发现,上述分区段二维反演中,方案二(先左后右方案)的反演结果与二维理论模型的反演结果最为接近.根据图 6所示阻抗张量成像分析可知,右半剖面的0°主轴方位的显著度要远远强于45°主轴方位的显著度,意味着,我们在分区段反演时,先反演主轴方位显著度弱的所在的区段,然后在其基础上再反演主轴方位显著度强的区段.这个规律与分频段反演是一致的.分频段反演也是先反演深部的显著度弱的构造结果更好一些.

2 实测算例:郯庐断裂带南段深部电性结构特征研究

以上我们基于三维交错模型的两种典型模型:纵向交错模型(立交桥模型)和横向交错模型,详细对比分析了交错模型下不同主轴方位选择方式的二维反演结果,提出了分频段反演方法和分区段反演方法,为交错模型的二维反演提供新的技术手段.但这两种理论模型毕竟是简单的,在实际复杂的地质构造条件下,我们提出的这种反演技术,是否依然有效呢?以下我们通过实测数据来进行检验.

2.1 电性主轴的确定

实测数据来自郯庐断裂带南段的一条大地电磁宽频带探测剖面,其测点分布情况见图 7.在本研究系列(五)(陈小斌和郭春玲,2017)中,已对其进行了全面详细的阻抗张量成像分析,并得出该实测剖面所经区域存在复杂的交错式结构.本文将在该文分析结果的基础上,进一步比较常规二维反演与分频段-分区段二维反演之间的差异.

图 7 大地电磁实测剖面测点分布图(图中菱形点为大地电磁测点,蓝色测点为本文所采用,灰色测点未采用;红色圆点为历史地震分布;粗实线为活动断裂) Fig. 7 Map showing active faults (thick solid lines), historical earthquakes (red dotes) and MT sites (diamonds, blue ones are used in this work)

根据我们先前的分析(陈小斌和郭春玲,2017),得到了如图 8a所示的阻抗主轴分布结果图.依据实感应矢量分布图(如图 8b)我们可以排除阻抗张量主轴方位的90°模糊性,结果如图 9a所示.由图 9a可见,在剖面的西北段,原阻抗主轴分布图上的0.015~3 Hz频段又可分出两个频段:0.15~3 Hz频段北东向主轴分布和0.015~0.15 Hz的北西向主轴分布,二者由于相差大致为90°而在阻抗主轴分布图上不能分辨出来.

图 8 (a) 为依据阻抗张量成像确定的电性主轴分布图(据陈小斌和郭春玲,2017);(b)为实感应矢量分布图,取向上为北 Fig. 8 (a) Geo-electric strike pseudo-section determined by impedance tensor imaging (Chen and Guo, 2017); (b) Induction vector pseudo-section (arrows upwards point to north)

为了能够获得可靠的反演结果,还需要考虑主轴方位与剖面展布方位之间的关系.一般来说,在剖面设计过程中,会考虑到测区主要构造的走向而使测线尽量与之垂直布设.在本文中,郯庐带南段的测线展布方位为NE115°左右,基本与郯庐带的北北东的走向垂直,但与其邻近地区一些构造的主轴方位夹角较小(< 30°),此时单纯以电性主轴方位为反演中的数据旋转方位,结果的可靠性会降低.一般要求测线展布方位与反演主轴方位的夹角要大于45°,以保证TE、TM极化模式的正确性.而在主轴方位的选择而言,陈小斌等(2008)曾经开展过相关的研究,对于纯二维模型,当反演中数据旋转方位与电性主轴方位相差不超过15°时,反演结果的可靠性是可以保证的,而对于三维模型响应或者实际数据而言,此条件还可以放宽.因此,我们规定以下两条主轴选择的标准:

(1) 反演中主轴方位与测线的夹角不小于45°;

(2) 反演中主轴方位与阻抗张量成像确定的主轴方位的偏差不超过30°.

在上述两个标准的约束下,我们将图 9a所示的电性主轴方位分布图进一步修正为图 9b所示的反演用的电性主轴分布.至此,我们完成了二维反演所用的主轴方位的选取方案,接下来将在图 9b的主轴方位下,开展交错模型的二维反演.

图 9 (a) 为参考实感应矢量分布图后修正的阻抗张量主轴方位分布图;(b)为考虑测线展布方位和主轴偏差容量以后所获得的二维反演所用的电性主轴方位分布图 Fig. 9 (a) Modified pseudo-section of geo-electric strike with respect to induction vectors; (b) Geo-electric strike for two-dimensional inversion obtained considering the direction of the MT survey line and the tolerance of the strike deviation
2.2 二维反演及其结果

在理论交错模型的二维反演研究中已归纳出,分频段反演,宜先反演低频段数据,将其结果作为初始模型进行高频反演;分区段反演,宜先反演主轴方位显著性弱的区段,再反演主轴方位显著性强的区段.从上述电性主轴分析中已经看出,这条实测剖面的电性结构可视为一个复杂的交错模型,既需要分频段反演,又需要分区段反演.针对这一实测数据,在图 9b所示的主轴方位基础上,将分频段反演和分区段反演组合起来,其反演过程与结果如图 10所示,图中还给出了常规二维反演结果(图 10f)以为参照.显然,分频段-分区段二维反演较常规二维反演在操作上更为复杂,不过在MT-Pioneer软件中却易于实现,图 10所示的所有反演过程在一天之内即可全部完成.

图 10 实测数据二维反演结果比较 (a) 0.0001~0.015 Hz频段、所有测点NE60°主轴方位的反演结果;(b)以结果(a)为初始模型,0.015~320Hz频段、剖面东南段(测点12~32))NE25°主轴方位数据的反演结果;(c)以结果(b)为初始模型,0.015~0.15 Hz频段、剖面西北段(测点1~11)NE-20°主轴方位的反演结果;(d)以结果(c)为初始模型,0.15~3 Hz频段、剖面西北段(测点1~11)NE45°主轴方位的反演结果;(e)以结果(d)为初始模型,3~320 Hz频段、剖面西北段(测点1~11)NE0°主轴方位的反演结果,亦为分频段-分区段反演的最终结果;f为以NE35°为主轴方位的常规二维反演结果.图中C1—C3代表低阻体,R1—R3代表高阻体. Fig. 10 Two-dimensional inversion results of measured data (a) Frequency band data (0.0001~0.015 Hz) of all sites based on the corresponding geo-electric strike of NE 60°; (b) Initial model using (a), 0.015~320 Hz frequency band, sites 12 to 32, strike azimuth NE25°; (c) Initial model using (b), 0.015~0.15 Hz frequency band, sites 1 to 11, strike azimuth -20°; (d) Initial model using (c), 0.15~3 Hz frequency band, sites 1 to 11, strike azimuth 45°; (e) Initial model using (d), 3~320 Hz frequency band, sites1 to 11, strike azimuth 0°. and sub-band and sub-segment inversion method used; (f) Conventional 2D inversion with the geo-electric strike of NE 35° for all sites and all frequencies. In these figures, C1-C3 indicate low resistance bodies and R1-R3 indicate high resistance bodies.
2.3 对比分析与讨论

图 10可知,分频段-分区段反演结果(图 10e)与常规反演结果(图 10f)在浅部(< 10 km)具有较好的相似性,测区的一些标志性构造,如张八岭隆起(R1)、合肥盆地(C1)及其西侧基底高阻(R2)等都比较清楚.不过,作为该剖面探测目标、测区主要构造的郯庐断裂带,二者的反映存在较大的差异.常规的单主轴反演结果中,基本上看不出类似于断裂带特征的结构,仅隐约在张八岭隆起高阻体的西侧有所显示(图 10f中的C3),而在分频段-分区段反演结果中,类似于断裂带特征(推测为郯庐断裂带)的显示非常清晰:无论是中上地壳高阻层中的可能的分支断裂,还是深部延伸到地幔的高导通道(C3),都得到较好地分辨.

分频段-分区段反演结果虽然能提供更为丰富的信息,但这些信息是否是可靠的呢?由于无法获得岩石圈电性结构的确切信息,因而这个问题无法给出确切的答案.我们只能从以下几个方面进行一定程度的分析,力求使我们的结论能够建立在比较好的逻辑推理之上.

首先,从反演结果本身来看,图 10中分频段-分区段的反演结果更加符合我们对郯庐断裂带深部结构的预设的认识,更易于解释.常规方法中,15 km以下的深部基本没有分辨率,郯庐断裂带在深部无法分辨,这不符合我们对一条曾经左旋走滑数百公里的巨型断裂带的深部结构的经验认识.不过,我们也知道在科学研究中预设的认识经常是错误的.

此外,在常规反演的结果中,合肥盆地C1正下方出现了一个在本文中似曾相识的八字形的高阻体(R3),很可能正是这个八字形高阻体的存在使得该反演结果在深部失去了分辨能力.在理论模型反演结果(图 6)的比较中,我们已经知道,当数据旋转角与高导线性体的主轴方位存在偏差时,在高导体的下方会出现八字形假高阻异常.图 6结果表明,即使旋转角与高导体的主轴方位仅相差22.5°,其下方的八字形假高阻异常也已非常明显.而图 10f中的R3上方正好为合肥盆地的高导区,按照阻抗张量成像分析,合肥盆地的电性主轴为南北向,而在图 10f的反演中,数据旋转角为NE35°,二者相差35°,这足以引起其下方出现较大规模的八字形假高阻异常.类似的八字形高阻体我们在其他地区、其他剖面的反演中也曾多次遇到,一直被视为浅部高导区的高阻基底,现在看来其很可能是上述原因导致的假构造.因此,从这方面看,分频段-分区段的反演结果要更为可靠一些.

图 10中我们还注意到,在新方法结果中,合肥盆地西部下方的高导异常体C2位于测点07的正下方,与本研究系列(五)(陈小斌和郭春玲,2017)中阻抗张量成像分析中该处线性构造较显著的结论相一致.而反观常规结果中,C2最显著处位于测点05的正下方,而该处未见有明显线性构造分布.

郯庐断裂带是中国东部一条巨型的断裂带,曾被视为大陆裂谷,其内部存在高大地热流异常(王良书等,1995祖金华等,1996)、低速异常带和岩石圈地幔的隆升(刘福田等,1986朱光等,2002),这些现象与新方法反演结果更为吻合.朱光等(2002)认为郯庐带在白垩纪的大规模走滑运动是以莫霍面为大型滑移拆离带进行的,壳幔之间发生了显著的解耦.此结论与我们依据阻抗张量成像技术确定的郯庐带附近壳幔间可能存在不整合接触的结论是吻合的(陈小斌和郭春玲,2017).此外,地震各向异性研究结果认为,郯庐带东侧扬子地块的快波偏振方向为NEE—SWW向,其西侧华北地块的快波偏振方向为E—W向,二者与地表构造存在一定的夹角(吴萍萍等,2012),这也与本文在低频部分所采用的NEE—SWW电性主轴方位一致.

通过以上分析来看,大地电磁分频段-分区段二维反演结果较之常规固定轴二维反演结果具有更高的分辨率和可靠性.除了能提供反演结果以外,我们还可以依据不同频段、不同区段的反演结果来查看相应数据在深度和剖面距离上的影响范围,可以为大地电磁数据解释提供一些新的参考信息.

3 结语

交错模型是由不同走向的二维模型组合而成的一类特殊的三维模型,是实际情形下经常可能会遇见的一种大地构造.由于这类模型的三维特性,不适于采用传统的大地电磁二维反演技术进行反演.为了解决此一问题,本文提出了分频段-分区段二维反演方法.通过理论模型和横跨郯庐带的一条大地电磁实测剖面的反演结果对比分析,展示和讨论了该反演方法的实施过程和反演结果的可靠性和分辨率.

分频段-分区段方法的基本思想是将交错模型中不同构造走向的二维模型响应分离出来,在二维反演中对该部分数据选择对应的电性主轴方位进行反演,除了第一次反演采用传统的初始模型处理方式以外,其他反演的初始模型是其他响应的反演结果,最终自动形成所有数据的反演结果.采用阻抗张量成像方法,结合感应矢量,可以得到不同区段、不同频段的电性主轴方位,从而可以从三维交错模型的综合响应中分解出单个二维模型的响应数据.理论模型的研究表明在进行大地电磁分频段-分区段二维反演时,对于纵向交错模型(立交桥模型),宜先反演低频数据,再反演高频数据;而对于横向交错模型,宜先反演主轴方位显著度低的区段,再反演主轴方位显著度高的区段.理论交错模型、横跨郯庐带实测剖面的分频段-分区段反演与常规反演的对比研究结果表明,采用本文提出的分频段-分区段二维反演方法较常规单主轴反演方法可以获得具有更高的可靠性和分辨率的二维反演结果.

References
Anahnah F, Galindo-Zaldívar J, Chalouan A, et al. 2011. Deep resistivity cross section of the intraplate Atlas Mountains (NW Africa):New evidence of anomalous mantle and related Quaternary volcanism. Tectonics, 30: TC5014. DOI:10.1029/2010TC002859
Anahnah F, Galindo-Zaldívar J, Chalouan A, et al. 2012. Reply to the comment by A. G. Jones et al. on "Deep resistivity cross section of the intraplate Atlas Mountains (NW Africa):New evidence of anomalous mantle and related Quaternary volcanism". Tectonics, 31(5). DOI:10.1029/2012TC003116
Cai J T, Chen X B. 2010. Refined techniques for data processing and two-dimensional inversion in magnetotelluric (Ⅱ):which data polarization mode should be used in 2D inversion. Chinese J. Geophys., 53(11): 2703-2714. DOI:10.3969/j.issn.0001-5733.2010.11.018
Cai J T, Chen X B, Zhao G Z. 2010. Refined techniques for data processing and two-dimensional inversion in magnetotelluric Ⅰ:Tensor decomposition and dimensionality analysis. Chinese J. Geophys., 53(10): 2516-2526. DOI:10.3969/j.issn.0001-5733.2010.10.025
Chen X B, Cai J T, Wang L F, et al. 2014. Refined techniques for magnetotelluric data processing and two-dimensional inversion (Ⅳ):Statistical image method based on multi-site, multi-frequency tensor decomposition. Chinese J. Geophys., 57(6): 1946-1957. DOI:10.6038/cjg20140625
Chen X B, Guo C L. 2017. Refined techniques for data processing and two-dimensional inversion in magnetotelluric (Ⅴ):Detecting the linear structures of the Earth by impedance tensor imaging. Chinese J. Geophys., 60(2): 766-777. DOI:10.6038/cjg20170227
Chen X B, Zhang X, Hu W B. 2000. Application of finite-element direct iteration algorithm to MT 2-D forward computation. Oil Geophysical Prospecting, 35(4): 487-496.
Chen X B, Zhao G Z, Ma X. 2008. Preliminary discussion on selecting rotation direction in 2-D MT inversion. Oil Geophysical Prospecting, 43(1): 113-118, 128.
Chen X B, Zhao G Z, Zhan Y. 2004. A visual integrated windows system for MT data process and interpretation. Oil Geophysical Prospecting, 39(Suppl.): 11-16.
Jones A G, Kiyan D, Fullea J, et al. 2012. Comment on "Deep resistivity cross section of the intraplate Atlas Mountains (NW Africa):New evidence of anomalous mantle and related Quaternary volcanism". Tectonics, 31: TC5011. DOI:10.1029/2011TC003051
Kiyan D, Jones A G, Hogg C, et al. 2010. MT investigations over the Atlas Mountains of Morocco: Preliminary results. //20th Workshop on Electromagnetic Induction in the Earth. Giza, Egypt.
Ledo J, Jones A G, Siniscalchi A, et al. 2011. Electrical signature of modern and ancient tectonic processes in the crust of the Atlas mountains of Morocco. Phys. Earth Planet. Inter., 185(3-4): 82-88. DOI:10.1016/j.pepi.2011.01.008
Liu F T, Qu K X, Wu H, et al. 1986. Seismic tomography of north China region. Acta Geophysica Sinica, 29(5): 442-449.
Mackie R L, Madden T R, Wannamaker P E. 1993. There-dimensional magnetotelluric modeling using difference equations:Theory and comparisons to integral equation solutions. Geophysics, 58(2): 215-226. DOI:10.1190/1.1443407
Rodi W L, Mackie R L. 2001. Nonlinear conjugate gradients algorithm for 2-D magnetotelluric inversion. Geophysics, 66(1): 174-187. DOI:10.1190/1.1444893
Wang L S, Li C, Shi Y S, et al. 1995. Distributions of geotemperature and terrestrial heat flow density in lower Yangtze area. Acta Geophysica Sinica, 38(4): 469-476.
Wu P P, Wang C Y, Ding Z F, et al. 2012. Seismic anisotropy of upper mantle beneath the Dabie-Sulu and its adjacent areas. Chinese J. Geophys., 55(8): 2539-2550. DOI:10.6038/j.issn.0001-5733.2012.08.006
Zhang G W, Meng Q R, Lai S C. 1995. Structure of Qinling orogenic belt. Science in China (Series B), 25(9): 994-1003.
Zhang G W, Meng Q R, Liu S F, et al. 1997. Huge intracontinental subduction zone at south margin of North China Block and present 3-D lithospheric framework of the Qinling orogenic belt. Geological Journal of China Universities, 3(2): 129-143.
Zhu G, Niu M L, Liu G S, et al. 2002. Structural, magmatic and sedimentary events of the Tan-Lu fault belt during its early cretaceous strike-slip movement. Acta Geologica Sinica, 76(3): 325-334.
Zu J H, Wu Q F, Lian Y F. 1996. The geothermal study of the Mid-segment of the Tancheng-Lujiang fault zone and its neighboring region. Earthquake Research in China, 12(1): 43-48.
蔡军涛, 陈小斌. 2010. 大地电磁资料精细处理和二维反演解释技术研究(二)-反演数据极化模式选择. 地球物理学报, 53(11): 2703-2714. DOI:10.3969/j.issn.0001-5733.2010.11.018
蔡军涛, 陈小斌, 赵国泽. 2010. 大地电磁资料精细处理和二维反演解释技术研究(一)——阻抗张量分解与构造维性分析. 地球物理学报, 53(10): 2516-2526. DOI:10.3969/j.issn.0001-5733.2010.10.025
陈小斌, 蔡军涛, 王立凤, 等. 2014. 大地电磁资料精细处理和二维反演解释技术研究(四)——阻抗张量分解的多测点-多频点统计成像分析. 地球物理学报, 57(6): 1946-1957. DOI:10.6038/cjg20140625
陈小斌, 郭春玲. 2017. 大地电磁资料精细处理和二维反演解释技术研究(五)——利用阻抗张量成像识别大地线性构造. 地球物理学报, 60(2): 766-777. DOI:10.6038/cjg20170227
陈小斌, 张翔, 胡文宝. 2000. 有限元直接迭代算法在MT二维正演计算中的应用. 石油地球物理勘探, 35(4): 487-496.
陈小斌, 赵国泽, 詹艳. 2004. MT资料处理与解释的Windows可视化集成系统. 石油地球物理勘探, 39(Suppl.): 11-16.
陈小斌, 赵国泽, 马霄. 2008. 关于MT二维反演中数据旋转方向的选择问题初探. 石油地球物理勘探, 43(1): 113-118, 128.
刘福田, 曲克信, 吴华, 等. 1986. 华北地区的地震层面成象. 地球物理学报, 29(5): 442-449.
王良书, 李成, 施央申, 等. 1995. 下扬子区地温场和大地热流密度分布. 地球物理学报, 38(4): 469-476.
吴萍萍, 王椿镛, 丁志峰, 等. 2012. 大别-苏鲁及邻区上地幔的各向异性. 地球物理学报, 55(8): 2539-2550. DOI:10.6038/j.issn.0001-5733.2012.08.006
张国伟, 孟庆任, 赖绍聪. 1995. 秦岭造山带的结构构造. 中国科学(B辑), 25(9): 994-1003.
张国伟, 孟庆任, 刘少峰, 等. 1997. 华北地块南部巨型陆内俯冲带与秦岭造山带岩石圈现今三维结构. 高校地质学报, 3(2): 129-143.
朱光, 牛漫兰, 刘国生, 等. 2002. 郯庐断裂带早白垩世走滑运动中的构造、岩浆、沉积事件. 地质学报, 76(3): 325-334.
祖金华, 吴乾蕃, 廉雨方. 1996. 郯庐断裂带中段及邻区的地热研究. 中国地震, 12(1): 43-48.