佳木斯地块及东缘岩石圈电性结构特征

佳木斯地块及东缘是中国东北地区的重要地质构造单元. 为探测该地区地壳深部结构与构造关系，沿桦南—饶河实施了240 km的大地电磁测深的探测研究. 采用光滑模型二维反演方法对桦南—饶河大地电磁剖面的探测数据进行了二维反演和综合地球物理解释. 研究结果揭示了研究区断面域的电性结构：(1)剖面西段具有稳定的高阻特征且具有稳定的岩石圈厚度（约90km），在十几公里深度范围存在壳内高导层，该区段对应佳木斯地块. (2)剖面中部具有明显的电性梯度带，该梯度带反映了佳木斯地块的东界位置及其深部的结构形态. (3)剖面东段电性特征揭示了佳木斯地块以东区段浅部为逆冲推覆体，深部为多个高阻块体和低阻条带相间的电性结构，这些高阻块体可能为早期俯冲的岩石圈残片.     

北京-柏各庄剖面连续介质大地电磁反演解释

在华北地区大地电磁测深工作的基础上,本文对北京——柏各庄剖面上的大地电磁观测资料,使用连续介质反演方法进行反演解释,获得了电性结构的最新结果。对构制的可接受地电模型的评价表明,本剖面内地壳中部普遍存在高导电带,上地幔内高导电带的深度变化较大,电阻率为1m的等值线深度,在沧县隆起为70多km,其西侧和东侧分别在60km和50km左右。上地幔高导电带的电阻率值由西而东有减小的趋势。这些结果与表层构造的对应关系说明,上地幔软流圈及地热场的横向不均一性,对地壳内部发生的构造运动有重要意义。      

青藏高原岩石圈三维电性结构

本文报道通过综合大地电磁调查数据研究青藏高原岩石圈三维电阻率模型的初步成果.大地电磁法调查区域已经覆盖了高原大部分面积,为全区三维电阻率成像研究打下了可靠的基础.对多个测区大地电磁数据进行精细的同化处理和反演成像,取得了青藏高原可靠的岩石圈三维电阻率结构图像.成像的区域为28°N—35°N,80°E—104°E.三维反演计算时采用的网格尺寸为20km×20km,垂直方向不等间距剖分为26层.结果表明,青藏高原现今岩石圈电阻率扰动主要反映印度克拉通对亚欧大陆板块俯冲引起的热流体运动和大陆碰撞和拆离产生的构造.在岩石圈地幔,察隅地块、喜马拉雅地块和拉萨地块东部联成统一的高电阻率地块,它们反映了向北东俯冲的印度克拉通.雅鲁藏布江、班公—怒江和金沙江缝合带都有明显的低电阻率异常,表明岩石圈深处有热流体活动.雅鲁藏布江、班公—怒江和金沙江缝合带都有明显的低电阻率异常,也表明它们的岩石圈还有流体活动.青藏高原东部的低阻区沿100°E向地幔下方扩大,反映了金沙江断裂带有切穿岩石圈的趋势.地幔电阻率平面扰动的模式显示,青藏高原东西部的地体碰撞拼合形式和方向是不同的.在青藏高原西部,羌塘、拉萨和喜马拉雅等地体从北到南碰撞拼合.在青藏高原东部,羌塘—拉萨、察隅、印支、雅安和扬子等地体多方向拆离拼合,在地壳造成不正交的拆离带和压扭构造系.从高阻-低阻区的分布看,东部的地体拼合有地幔的根源,今后还会进一步发展.察隅地块岩石圈对青藏高原东部的楔入,使其北部和东部地块的岩石圈发生拆离撕裂,也造成热流体上涌的低电阻率异常.     

稀疏测线大地电磁资料三维反演研究:合成算例(英文)

受勘探成本和工区环境等因素的影响,当前大多数大地电磁实际工作采取布置稀疏测线采集数据和使用二维反演方法解释这些稀疏测线数据的方式。然而,二维反演方法在解释三维地电构造数据时存在局限性,有时甚至做出错误的地质解释。本文尝试了使用三维反演方法对大地电磁稀疏测线数据进行反演解释。使用大地电磁全信息资料三维共轭梯度反演程序对理论模型合成稀疏测线数据进行了三维反演。结果表明:这种反演方案是可行与有效的。同时,我们发现在不同数据的三维反演结果中,四个张量阻抗元素和两个倾子数据同时反演的结果相对更为准确,更接近理论模型。     

松辽盆地岩石圈减薄的深部动力学过程

松辽盆地作为东亚裂谷系的一部分,与华北克拉通一起经历了中生代岩石圈减薄的重大地质事件.对大陆岩石圈-软流圈状态和构造的整体认识,是研究大陆岩石圈减薄深部动力学过程的关键.在获得过松辽盆地的106个宽频和30个长周期大地电磁测深数据的基础上,完成测点数据二维偏离度、构造走向等计算与分析,进一步采用非线性共轭梯度算法,对TE和TM模式数据进行二维联合反演,获得了沿剖面的壳-幔电性结构,并依此构建了松辽盆地壳-幔结构模型.研究结果表明：（1）大兴安岭地区岩石圈厚度约为160 km,松辽盆地岩石圈厚度约为45 km,张广才岭岩石圈厚度在70~100 km之间,莫霍面与岩石圈底界面不呈镜像关系.软流圈整体表现为中、低阻异常,电阻率值在30 Ωm左右,其形态呈西倾约30°的蘑菇状异常,指示了软流圈物质上涌的形式,有别于软流圈垂直上涌的传统认识.（2）松辽盆地深部存在双层高导异常（电阻率小于5 Ωm）,上层为壳内高导层,呈"蛇"状分布,推断为岩浆底侵区,下层为幔内高导层,呈"哑铃"状,为软流圈上涌区.软流圈内存在两个"哑铃"状中、高阻异常,推断为拆沉的岩石圈地幔.具有冷的、高密度的下降物质流的堆积以及拆沉块体下插到两侧山岭是促使大兴安岭与张广才岭在中生代伸展环境中快速隆升重要原因;（3）松辽盆地经历了岩石圈减薄事件,与大兴安岭岩石圈厚度相比,松辽盆地岩石圈厚度减薄了近100 km,与东侧张广才岭相比减薄了70 km,而与中生代华北地台100 km的岩石圈厚度相比,减薄了近50 km,其经历了岩石圈伸展期、裂解期、拆沉期和增长期的动力学过程.     

可控源音频大地电磁三维共轭梯度反演研究

可控源音频大地电磁法在资源勘探等领域中发挥着重要的作用.我们把有限差分数值模拟方法用于可控源音频大地电磁三维正演,结合正则化反演方案和共轭梯度反演的思路,将反演中的雅可比矩阵计算问题转为求解两次"拟正演"问题,得到模型参数的更新步长,形成反演迭代,实现了可控源音频大地电磁三维共轭梯度反演算法.该反演算法可用于对有限长度电偶源激发下采集到的可控源音频大地电磁全区(近区、过渡区和远区)视电阻率和相位资料进行三维反演定量解释,获得地下三维模型的电阻率结构.理论模型合成数据的反演算例验证了所实现的可控源音频大地电磁三维共轭梯度反演算法的有效性和稳定性.     

利用大地电磁三维反演方法获得二维剖面附近三维电阻率结构的可行性

大地电磁野外实测数据月前大多为二维剖面数据.如何反演这些二维剖面数据获得较为接近实际地电情况的结果,是多数大地电磁工作者关心的问题.我们通过对理论模型的三维响应进行分析和对合成数据及实测资料的反演结果进行对比研究,讨论了利用三维反演的方法来获得大地电磁二维剖面附近三维电阻率结构的可行性.结果表明:可用三维反演的方法来解...     

滇东北地区三维电性构造研究

滇东北地区位处不同构造分区的交界地带,发育有多样的断裂构造,其中昭通—鲁甸断裂和小江断裂带的部分区段均具备积累高应力的条件,具有较强的孕震潜力.本文基于覆盖巧家、东川、鲁甸等地区采集的区域大地电磁阵列观测数据,通过大地电磁三维反演,获得了该地区的三维电性结构的水平和垂直切片,完整地揭示了滇东北地区的深部电性结构特征.反演结果表明,昭通—鲁甸断裂在10 km深度以内发育有多个局部高阻区,小江断裂带的巧家—东川段的东侧、东西两支交汇处在上地壳均表现为高阻,这些区域均展示出更强的孕震能力.根据低阻异常的分布,推测昭通—鲁甸断裂和小江断裂下方的低阻构造在中下地壳的多个位置发生连通.本文将研究区的地表热流、岩石圈有效弹性厚度、地震与电性结构结合分析,推测在滇东北深部低阻层可能存在下地壳流,并且下地壳流的塑性流动为鲁甸的地震活动提供了深部动力来源.     

华南地区岩石圈电性特征及其地球动力学意义

岩石圈-软流圈界面（The Lithosphere-Asthenosphere Boundary,LAB）是地球内部主要界面之一.大地电磁测深（Magnetotelluric,MT）是研究地球壳幔电性结构最有效的方法,利用长周期大地电磁测深数据可以较好地探测LAB.在SinoProbe-01-03课题的资助下,首次获得了华南地区4°×4°网度的高质量大地电磁测深数据.利用一维奥卡姆（Occam）算法反演了MT阻抗的反对角线元素所计算出的平均视电阻率.根据一维地电结构可以将华南地区岩石圈划分为五种类型：以湖南邵阳和贵州施秉为代表的克拉通型,以四川达州和彭州及湖北荆门为代表的构造边界型,以浙江湖州和广东云浮为代表的岩石圈中等改造型,以江西赣州、广东揭阳及福建霞浦为代表的岩石圈强烈改造型,以湖北英山为代表的造山带型.除湖南邵阳、贵州施秉及广东揭阳外,华南地区岩石圈厚度为60~145 km.本文研究表明华南地区岩石圈显示出南北两侧上抬、中部下凹、东部受不均匀改造的趋势,这一结果与之前发表的文献所揭示的华南地区岩石圈东薄西厚的典型特征是不同的.研究结果反映华南地区岩石圈稳定性较好,晚中生代以来的构造伸展作用对岩石圈的改造程度有限,可能主要以不同形式的软流圈底辟为主.     

基于递归分析和聚类的大地电磁信噪辨识及分离

为了剖析大地电磁信号和强干扰的本质特征,进一步精细分离出微弱的大地电磁有用信号,提出基于递归分析和聚类的大地电磁信噪辨识及分离方法.首先,运用递归分析法扩展大地电磁一维时间序列的维数,分析了嵌入维数、延迟时间和判别阈值对递归图的性能,并研究了不同长度的序列对递归定量分析参数的影响情况,然后,构建典型的大地电磁强干扰类型和微弱的大地电磁有用信号样本库,针对样本库讨论了强干扰和微弱大地电磁信号之间的递归定量分析参数,分析了K均值聚类和模糊C均值聚类的信噪辨识效果.最后,对实测大地电磁数据进行信噪辨识处理,并仅对辨识为强干扰的时间段采用数学形态滤波进行噪声压制.实验结果表明,递归分析能定性及定量地描述大地电磁信号时间序列的非线性特征和原动力系统的本质规律,与聚类算法相结合能对矿集区实测大地电磁信号进行信噪辨识;处理后的卡尼亚电阻率-相位曲线更为光滑、连续,其结果更为精细地保留了大地电磁信号低频段的缓变化信息,整个低频段的大地电磁数据质量得到了明显改善.     

地震资料约束下大地电磁资料反演

应用地震资料约束下大地电磁资料反演古潜山或基底内部结构的新方法. 首先利用地震及井资料,反演得到古潜山或基底之上各层的厚度及电阻率,建立了古潜山或基底以上的地电模型;并利用大地电磁一维反演资料,建立了古潜山或基底及其以下的地电模型. 在初始模型建立的基础上,应用高效快速的二维大地电磁正演算法——模式匹配方法,进行正演计算. 利用共扼梯度迭代方法,进行正反演迭代计算. 通过对胜利油田南北618剖面反演, 基底以下在地震资料没有显示的构造信息,在反演结果中得到了较好地反映.      

大地电磁资料精细处理和二维反演解释技术研究(五)——利用阻抗张量成像识别大地线性构造

大地线性构造往往与块体或物性的边界、断裂的几何延伸相关联,是地学研究的重要目标构造.在大地电磁法中,传统的方法是通过反演得到的电性结构的差异来解译线性构造,是一种间接分析方法.由于线性构造发育处走向和倾向明显,二维性较其他地方一般会更强,故而在阻抗张量成像分析中,相应的二维有效因子就会比较大,最佳主轴统计分布会更加明显.据此,本文利用大地电磁阻抗张量成像(包括最佳主轴统计成像、构造维性成像)来探测识别大地线性构造.本文首先通过三维理论模型的研究结果来验证这一技术在线性结构探测中的有效性,然后将其应用于郯庐断裂带南段实测剖面的具体分析中,展示了这一技术在线性构造探测中的实用性及其数据处理与分析过程.本文研究表明,大地电磁阻抗张量成像技术可独立于反演为地质解释尤其是线性构造的解释提供诸多参考信息.     

基于交叉梯度耦合的大地电磁与地震走时资料三维联合反演

为了有效解决目前大地电磁和地震走时资料单方法反演结果一致性不好的问题,同时克服基于岩石不同物性参数间关系耦合约束联合反演的局限性,本文研究了基于交叉梯度耦合约束的大地电磁与地震走时资料的三维联合反演算法.以较为成熟的天然地震走时资料三维正反演和大地电磁三维正反演算法为基础,实现了具有共同的反演网格,以交叉梯度结构耦合约束,并能同时获得电阻率和速度模型的三维联合反演算法.分别利用单棱柱体模型和双棱柱体模型合成数据进行了联合反演试算.结果表明:无论是单棱柱体模型还是双棱柱体模型,联合反演结果比单独反演对异常体的空间形态都有更好的恢复,其中单棱柱体模型反演的异常体电阻率更接近于真实电阻率,双棱柱体模型的联合反演结果不仅消除了围岩的部分电阻率假异常,而且增强了对异常体深部速度结构特征的恢复程度.联合反演还能同时改善电阻率和速度反向变化异常体的单独反演结果,进一步证明交叉梯度耦合不依赖于岩石物性关系,而强调地下结构的相似性,具有更普遍的适用性.     

青藏高原东南缘岩石圈有效弹性厚度及其构造意义

本文利用三维有限差分方法,基于EIGEN6C4布格重力异常和SIO V15.1地形数据,计算了青藏高原东南缘岩石圈有效弹性厚度.结果表明:青藏高原东南缘岩石圈有效弹性厚度为0~100 km,四川盆地和喜马拉雅东构造结岩石圈有效弹性厚度最大,达50~100 km;巴颜喀拉块体东部、川滇菱形块体大部、滇西等地区岩石圈强度弱,有效弹性厚度一般小于15 km;羌塘块体东部的玉树—德格附近地区岩石圈有效弹性厚度大于40 km;滇南地区岩石圈有效弹性厚度为10~30 km,大于云南北部地区.研究区域有效弹性厚度分布特征与岩石圈结构关系密切.四川盆地、喜马拉雅东构造结地区内部结构稳定,因而岩石圈强度大.川滇菱形块体等岩石圈有效弹性厚度小的地区与壳内低速、低阻/高导层分布有很好的对应关系,推测壳内岩石的部分熔融软化可能是造成高原东南缘岩石圈强度较弱的重要原因.羌塘块体东部的局部高力学强度岩石圈则可能是高原形成过程中的残留克拉通.根据本文计算的岩石圈有效弹性厚度特征,结合地震学、大地电磁等研究成果,认为青藏高原物质向东南缘挤出后受四川盆地等阻挡,造成下地壳软弱物质在理塘—稻城—丽江一带堆积,少部分物质可能穿过鲜水河断裂带的康定—道孚地区向北运动,但大部分物质向南运动,在受到滇南块体阻挡后一支流向西南的腾冲方向,另一支流向东南的攀枝花—东川方向.     

东海重磁地震综合探测剖面研究

依据东海南部一条自东海陆架直至菲律宾海盆的重磁地震综合探测剖面,采用地壳结构重磁地震综合反演解释系统开展岩石圈结构的综合研究,制作了岩石圈结构地学断面图. 结果表明基隆凹陷为一中新生代沉积凹陷,沉积厚度达到14km;冲绳海槽盆地为弧后扩张地堑型盆地,地壳厚度仅为14.5km;菲律宾板块沿北西西方向向欧亚板块俯冲,莫霍面急剧下插,导致冲绳海槽的形成与发展;岩石圈厚度约为80km,但在370km处仅为53km,在450km至540km处岩石圈厚度大于100km.     

大地电磁测深和电偶源频率电磁测深资料反演

在研究大地电磁响应函数频散关系的基础上,构制了一套滤波系数算法,以用于由一套视电阻率资料估算相应的阻抗相位.理论模型和实际大地电磁观测资料的数字试验表明,该法是行之有效的.由频散关系估算的相位值与观测的相位资料的比较,可用于检验观测资料是否满足频散关系.利用经频散关系校正的阻抗相位值,进行大地电磁阻抗的联合反演则可望获得更为可靠的结果.研究了观测频带相互衔接的电偶源频率电磁测深和大地电磁测深视电阻率的一维联合反演问题.对两个实测点两种电磁法的观测资料进行了联合反演试验,与钻井资料对比表明,所获得的电性分层参数是较为可信的.在补充了由频散关系获取的电偶源频率电磁相位资料后,对于两种电磁法的视电阻率————阻抗相位、阻抗实部视电阻率-阻抗虚部视电阻率进行了拟大地电磁反演,获得了相近的反演结果.      

中亚造山带东段软流圈分布特征——基于长周期大地电磁探测的结果

中亚造山带东段经历了古亚洲洋、蒙古—鄂霍茨克洋和古太平洋构造体制的多重作用,多期次的构造活动不仅改造了地形地貌、岩石圈结构,同时也改造了软流圈分布,因此对软流圈结构研究具有重要的科学意义.为了揭示复杂多重构造体系下软流圈的分布特征,完成了中亚造山带东段约1800km长大地电磁测深剖面,并对数据进行非线性共轭梯度反演,获得二维电性结构模型.研究发现:中亚造山带东段岩石圈内部存在多处低阻异常,地表多对应第四纪火山群或古缝合带,表明这些低阻异常与软流圈上侵有关.软流圈呈现高导特征,南部略薄,电阻率值为10～30Ωm,北部厚度较大,电阻率值为10～0.1Ωm,这种电性结构特征体现了中亚造山带东段软流圈南北两侧厚度、部分熔融程度的不均一性.而已有的层析成像结果显示中亚造山带软流圈东西向则呈现东厚西薄的特征.结合区域地质,提出软流圈与岩石圈一样在碰撞造山过程中厚度、部分熔融程度也会随之变化的动力学认识.中亚造山带东段软流圈主要经历了古亚洲洋构造体系、蒙古—鄂霍茨克构造体系和太平洋构造体系三阶段的构造事件影响,因此中亚造山带东段软流圈的南北向差异,推断为古亚洲洋闭合早于鄂霍茨克洋闭合的时限差异所致,东西向差异则主要受太平洋构造体系的影响.     

大地电磁全信息资料三维共轭梯度反演研究(英文)

在对张量阻抗数据、倾子数据和共轭梯度算法深入分析的基础上,我们实现了大地电磁全信息资料三维共轭梯度反演算法。基于全信息资料的三维共轭梯度反演研究,探讨了同时利用五个电磁场分量整理得到的大地电磁资料进行三维反演定量解释的方法以及全信息数据在三维反演中的作用。理论模型合成数据的反演结果表明,在三维反演中使用张量阻抗和倾子数据结合的全信息数据的反演结果优于只使用张量阻抗数据(或只使用倾子数据)的反演结果,提高了反演结果的分辨率和可信度。合成数据的反演算例也验证了所实现的大地电磁全信息资料三维共轭梯度反演算法的正确性和稳定性。     

花岗岩地区深部地热构造电磁法探测研究——以福建龙海地区为例

为探查福建龙海花岗岩区深部地热构造,试验音频大地电磁法(AMT)在此类地区的应用效果,在该区域布设了数条AMT测线.经过一系列的数据处理及反演,得到研究区的各条测线的电阻率深度剖面图及电阻率分布平面切片图.该测区AMT研究成果反映了地下2 km深度范围内的地质构造分布基本情况,指示了可能的深部地热构造位置,为深部构造探测提供了良好的地球物理基础资料.结果表明,音频大地电磁法在探测花岗岩地区深部地热构造方面是有效的.     

大地电磁测深在西安凹陷基底探测中的应用

为了研究西安凹陷深部地层及基底构造特征,采用大地电磁测深十字剖面法对西安凹陷的电性结构进行定性分析,并使用非线性共轭梯度(NLCG)二维反演方法结合地震勘探结果对其定量解释,划分出西安凹陷的主要地层及埋深.大地电磁测深反演结果与地震勘探解释结果相吻合,有效反映了西安凹陷深部地层结构的真实形态,为研究西安凹陷基底结构特征提供有效的依据.