苏鲁超高压变质地体被认为是目前世界上已知的最大的超高压变质带,是世界上研究陆-陆碰撞俯冲和壳幔作用的最佳地质场所,大陆动力学研究的基地。在超高压变质作用演化的研究中,超高压变质岩石的原岩和前超高压变质事件研究一直是薄弱环节。为了探讨高压-超高压变质地体的形成演化历史,有必要对前高压-超高压变质事件的历史(包括变质基底岩石的形成及其构造环境、主要构造事件的时限等问题)以及对超高压变质作用的影响与制约进行研究。目前对苏鲁高压-超高压变质带的构造研究,主要集中在建立折返构造格架、俯冲构造及流变学的重塑、俯冲与折返事件的时限,以及俯冲与折返机制的探索等。研究超高压变质带岩石的形成、折返和相应的动力学过程,需要了解该区地壳和上地慢构造。地球物理方法是研究地球内部结构与状态、岩石圈形成及演化历史的有效手段。尽管目前已有一些有关苏鲁超高压-高压变质带地区地球物理探测的研究成果报道,但是,其解释的深度,尤其是对区域重磁场进行全面深入的处理和解释还未见。地球物理位场是一种独立存在的天然场,利用位场资料进行研究,是以地壳内各种岩石间的密度或磁性差异引起的重磁异常为物理基础,探测地壳乃至岩石圈结构,特别是在研究覆盖大范围、多尺度地壳及岩石圈结构分析中起到重要作用。因此,在前人已经有的成果基础上,充分利用已有的地质、地震、大地电磁及CCSD(?)孔等资料,将区域重磁场的特征与地壳构造构造格架关系联系起来进行合理分析和解释,具有重要的科学意义与应用价值。本文研究结合国家973计划课题“大陆板块会聚边界深部物质物理状态的高温高压实验研究”之子课题“大陆板块会聚边界深部岩石磁性结构研究”(2003CB716506)的延续研究和国家自然科学基金重大研究计划项目“不同时空背景地幔属性与克拉通破坏过程”之子课题“壳-幔边界岩石磁性研究”(90714002)的需要选本研究主题。将以苏鲁高压超高压变质带区域重磁特征研究为主线,以前人丰富的构造地质、岩石学及地球化学与地震学研究成果为基础,以研究区重要地学问题—超高压变质岩深俯冲与折返动力学为目标,以几个具体地质问题基底性质(岩性、构造与几何形态)、重要构造界线(华北和扬子地台、深俯冲及郯庐断裂等)、Moho面与磁性下界面及地壳性质等研究为桥梁,探讨超高压变质过程区域重磁响应建立多尺度重磁场结构与地壳及岩石圈结构之间的关系,为陆陆碰撞的三维结构提供证据。本文研究主要内容包括：(1)调研目前地球物理位场数据处理的主要方法(包括多尺度分离与特征分析等)并进行评述；(2)地球物理非线性反演方法研究,分析传统的非线性反演算法的面临问题并提出改进思路,利用实例进行评价分析；(3)在前人已经有的成果基础上,充分利用已有的地质、地震、大地电磁及CCSD钻孔等资料,研究区几个具体地质问题：基底性质(岩性、构造与几何形态)、重要构造界线(华北地台、扬子地台、深俯冲及郯庐断裂等)、Moho面与磁性下界面性质及下地壳性质等研究为桥梁,探讨超高压变质过程的区域重磁响应,建立多尺度重磁场结构与地壳及岩石圈结构之间的关系。在地球物理反演计算中,特别注意到如果选择适当的初始模型,对求解非线性问题的线性化反演解法可以获得较满意的结果,而非线性全局优化方法是在模型空间中随机的或系统的搜索尝试的方法,它具有对于初始模型不依赖性,因此,如果能将非线性反演方法与线性化方法两者的优势结合起来以提高计算效果,可以成为一种地球物理非线性反演算法改进思路。本文提出几种(GA-GLS、PSO-GLS等)新的地球物理非线性反演算法,例如,GA-GLS是将遗传算法的随机全局搜索特性与广义最小二乘的局部精确搜索特性优势相结合提高计算效率。通过对理论模型数据和实测数据进行反复实验,计算结果统计显示,这些新的算法与单纯非线性反演算法相比,重复实验计算结果的统计均值较接近于真实值,统计方差较小,结果较稳定。因此,这种新的改进算法主要优点有：(1)继承了单纯的非线性优化算法所具有的全局搜索优势,(2)改善了单纯的非线性算法容易陷入局部极小的特性。其在准确性、稳定性和计算效率等方面优于单纯的非线性算法。分析研究区重力场特征,获得密度分布空间结构。(1)研究区岩石密度特征。元古宙的闪长岩、榴辉岩和各时代的蛇纹岩、橄榄岩为高密度,一般在2.70gcm-3以上,其中,金红石矿化榴辉岩密度值可达4.30gcm-3以上；晚白垩纪世的砾岩砂岩、砂质页岩为低密度,一般在2.45gcm-3左右；二辉榴辉岩和橄榄岩具有高波速和高密度特征,比研究区围岩片麻岩波速高出1100m/s,密度高出0.5gcm-3,在重力场上可引起正异常,因而在沿郯庐断裂附近的超高压变质带上均分布着较强的正异常。(2)通过计算获得Moho面起伏以及布格重力异常的长波长趋势特征显示,沿郯庐断裂带附近从深部往地表有正异常突起趋势,这可能是地壳地层结构存在一定程度隆起,正异常的东西两侧可能存在比较深的断裂、甚至延伸至Moho面附近。(3)在66km截止低通滤波异常图可见,莱阳盆地和沭阳盆地附近均存在一定程度的重力值下降,但在100km截止低通滤波异常图上,莱阳盆地附近重力变化已不太明显,这可能反映出莱阳盆地沉积覆盖较浅。(4)在区域重力场中,北苏鲁超高压构造片体呈现由西向东从正异常变为负异常,而南苏鲁超高压构造片则以正异常为主,这说明在南、北苏鲁超高压构造片体在地壳深部地层密度分布特征上确实存在较大差异。这里是否与南、北带主体分别属于扬子克拉通北缘带和华北克拉通南缘带存在联系需要进一步研究。分析区域航磁异常特征,获得超高压变质带地壳磁性结构。(1)研究区磁性岩石特征。一般规律为：中酸性岩石,如花岗岩、花岗闪长岩与花岗片麻岩区表现为低磁性特征；基性超基性岩大多分布规模小,引起的磁异常不明显；基底深变质岩表现出较强磁异常。(2)根据不同尺度磁异常特征,对研究区进行了构造分区。共划分出华北克拉通、扬子克拉通及苏鲁造山带三个级构造单元郯庐断裂带、鲁西台隆、徐淮断隆、苏北胶南地体和苏北断坳五个级构造单元。这些结果与区域地质构造基本一致。(2)研究区深部存在明显北西向分布的大尺度地质体,自鲁西台隆跨过郯庐断裂带,经胶南、日照进入黄海。磁性体顶界面埋深约为20km,在郯庐断裂以西略浅,以东略深,大致反映了西侧相对隆起,东侧相对下凹的特征。与鲁西出露的泰山群变质岩磁性相当,推测该磁性体属华北克拉通太古代变质岩,揭示了下地壳华北克拉通物质连续延伸到苏北,尽管苏北地区上地壳属扬子克拉通。(3)分析航磁异常的空间变化特征,发现苏鲁地区不同波长磁异常的场源具有空间立交桥式结构特征。自上而下依次为：上部以北东向为主的中等波长磁异常场源,中间以北西向为主的中长波长磁异常场源,下面以北东向为主的长波长磁异常场源。这种“立交桥”式的空间磁性结构,可能为两克拉通碰撞和拆离的结果,反映了扬子与华北克拉通碰撞的动力学模式。利用航磁资料进行磁化强度或视磁化率填图,有效地圈定磁性岩体,划分磁性地层,突出地质构造单元界线,进而研究区域地质构造。研究区苏鲁超高压-高压变质带地区,浅层有规律地断续分布有西北-东南向、东北-西南向排列的多条磁性岩带。对照苏鲁地区岩类分布,可见有侵入岩出露的地区,一般都与磁性填图中的弱磁性到较强磁性区相对应,侵入岩的分布方向与磁性填图分布特征基本一致。(a)华北板块与莱阳盆地的西北面相接处、淄博市以东南、潍坊市以南、沂源县以东、沂水县以北地区,出现西北-东南向走向串珠状异常,可能存在一些NW-SE走向隐伏断裂,该区磁性属于中等磁性；(b)北苏鲁拆离断层以北、郯庐断裂以东的莱阳盆地地区,存在的东北-西南向和东西向线性异常,这可能反映盆地边缘存在一些断裂,该区磁性属于中等磁性；(c)北苏鲁超高压构造片体,出现西北-东南走向线性异常,这可能反映北苏鲁超高压构造片体分布结构,该区磁性属于高磁性,这里可能充填的基性或超基性岩浆岩、太古一下元古代基底和中基性火山岩及高压变质岩,区内还不均匀分布着花岗岩；(d)南苏鲁超高压构造片体地区,出现南北走向线性异常,这可能反映北苏鲁超高压构造片体分布结构,该区磁性属于中等磁性,这里可能充填的基性或超基性岩浆岩；(e)苏鲁高压构造片体地区,出现东北-西南走向线性异常,这可能反映北苏鲁高压构造片体分布结构,该区磁性属于低磁性；(f)沭阳盆地地区,出现西北-东南走向线性异常,该区磁性属于低磁性；磁性构造层底面是地球岩石圈上部磁性壳层的底界面。这一界面深度(居里温度面)以上的磁性壳层岩石具有磁性,其下岩石磁性迅速消失。针对研究区航磁异常具有分区特征,且浅源、深源磁场分布复杂的特点,应用组合磁异常的统计功率谱分析法计算磁性构造层底面深度。利用航磁数据,计算获得了研究区磁性构造层底面深度,深度大概在18.5～27km范围内变化。研究区的磁性构造层底面的起伏趋势也表现出NE走向,表明在滤除浅源磁性体影响后,磁性构造层底面计算结果突出了深部热结构特征。研究区南部磁性构造层底面深度较浅(18.5～22km),尤其是苏北盆地,其居里点深度为18.5km。而在研究区北部,特别是胶南隆起,磁性构造层底面深度相对较深(24～27km)。诸城盆地地区磁性构造层底面深度约为23.5～25km,而胶辽附近深度为23.5～24km,苏鲁超高压变质带地区深度分别为21～21.5km。这些结果显示,磁性构造层底面起伏趋势与软流圈起伏具有明显的相关性。在中朝板块与扬子板块发生陆—陆碰撞时,扬子板块深俯冲过程中,拽动相邻的中朝板块的一部分物质向下俯冲,引起上覆板块物质遭受俯冲剥蚀,导致不对称的双向俯冲。板块俯冲受阻后,板块滞留、下沉引起地幔扰动,而地幔流又在老的构造活动所造成的地壳薄弱带上涌,从而导致该区的磁性构造层底面上隆。使用场源定位分析的特殊函数和欧拉反褶积方法分别对于研究区重磁资料进行场源定位,通过对比分析获得了比较准确和直观的场源空间分布情况,寻找出研究区中可能会出现的几种特殊地质体,如接触体、板状体边缘、水平线源、点源等,分析研究区构造关系。(a)构造指数SI=-1的垂直密度接触体空间分布特征：在研究区西北部的华北板块地区深度约为5～10km处,以西北-东南向走向为主要,水平分布紧凑；北苏鲁拆离断层以北地区地下深度约为5～7km,以处东北-西南向走向主要,水平分布宽松；苏鲁超高压构造片体深度约为5～7km处,以东北-西南向、西北-东南向等多个方向走向,空间几何形态分布复杂、纵横交错；苏鲁高压构造片体深度约为5～7km处,以东北-西南向走向为主,水平分布紧凑；沭阳盆地深度约为5～7km处,以东北-西南向和西北-东南向走向为主,水平分布均匀。(b)构造指数SI=0的水平密度板状体的空间分布特征：在研究区西北部的华北板块地区深度约为10～15km处,以东北-西南向走向为主要,水平分布宽松；北苏鲁拆离断层以北地区地下深度约为5～7km,以处东北-西南向走向主要,水平分布宽松；苏鲁超高压构造片体深度约为7～13处,以东北-西南向、西北-东南向等多个方向走向,空间几何形态分布复杂、纵横交错；苏鲁高压构造片体深度约为5～7km处,以东北-西南向走向为主,水平分布紧凑；沭阳盆地深度约为5～7km处,以东北-西南向走向为主,水平分布均匀。(c)构造指数SI=1的水平密度圆柱体的空间分布特征：在研究区西北部的华北板块地区深度约为9～12km处,反映沉积盆地边缘深部地质构造；北苏鲁拆离断层以北地区地下深度约为9～15km,以西北-东南向走向为主；苏鲁超高压构造片体深度约为12～16处,围绕郯庐大断裂带和苏鲁超高压构造片体,主要以南北向走向,这可能反映了郯庐大断裂带、苏鲁超高压构造片体边缘深部地质构造；苏鲁高压构造片体深度约为12～16km处,以东西向走向为主,反映了高压构造片体边缘构造；沭阳盆地深度约为10～13km处,以西北-东南向走向为主,围绕沭阳盆地边缘。(d)构造指数SI=2的密度球体的空间分布特征：在研究区西北部的华北板块地区深度约为18～20km处,以西北-东南向走向为主,反映华北板块在该区的深部构造；北苏鲁拆离断层以北地区地下深度约为13～16km,以东西向向走向为主,水平分布集中,反映了盆地边缘深部构造；北苏鲁超高压构造片体深度约为13～17km处,主要以东西向走向,这可能反映了超高压构造片体和边缘沉积盆地边缘深部构造；南苏鲁超高压构造片体地下深度约为18～20km处,围绕郯庐大断裂带和超高压构造片体,以东西向走向为主,这可能反映了郯庐大断裂带和超高压构造片体边缘深部地质构造。苏鲁高压构造片体深度约为12～17km处,围绕郯庐大断裂带和高压构造片体,主要以南北向走向；沭阳盆地地区深度约为14～17km处,以东西向走向为主,沿沭阳盆地边缘。(e)构造指数SI=0的垂直磁性接触体的空间分布特征：在研究区西北部的华北板块地区深度约为5～10km处,以西北-东南向走向为主,水平分布紧凑；苏鲁拆离断层以北地区地下深度约为4～7km,以东北-西南向和东西向走向为主；北苏鲁超高压构造片体深度约为4～7km处,沿郯庐大断裂带和苏鲁超高压构造片体边缘,主要呈现东北-西南向走向,空间几何形态复杂、纵横交错；苏鲁高压构造片体深度约为12～17km处,以南北向偏东北-西南向为主；沭阳盆地地区深度约为5～7km处,沿郯庐大断裂带和沭阳盆地边缘,以东西向和东北-西南向等走向交错。(f)构造指数SI=1的水平磁性板状体的空间分布特征：在研究区西北部的华北板块地区深度约为7～13km处,以西北-东南向为主,水平分布集中；苏鲁拆离断层以北地区地下深度约为10～12km km,以东北-西南向和东西向走向为主；北苏鲁超高压构造片体深度约为9～12km处,以东北-西南向走向为主,反映了北苏鲁超高压构造片体水平磁性板状体分布特征；南苏鲁高压构造片体深度约为5～10km处,以东北-西南向为主,水平分布紧凑；鲁高压构造片体深度约为7～12km处,沿郯庐大断裂带和苏鲁高压构造片体边缘；沭阳盆地地区深度约为17～21km处,沿郯庐大断裂带和沭阳盆地边缘,以东北-西南向、东西向走向为主。(f)构造指数SI=2.0的磁性水平柱体的空间分布特征：在研究区西北部的华北板块地区深度约为2～6km处,沿华北板块边缘和郯庐大断裂带,磁性水平柱体分布比较集中；北苏鲁拆离断层以北地区地下深度约为2～4km,沿沉积盆地边缘和郯庐大断裂带边缘附近,磁性水平柱体分布比较松散；苏鲁超高压构造片体深度约为2～6km处,沿苏鲁超高压构造片体和郯庐大断裂带边缘,分布集中；沭阳盆地地区深度约为3～5km处,沿郯庐大断裂带和沭阳盆地边缘,分布比较集中(h)构造指数SI=3.0的磁性球体的空间分布特征：在研究区西北部的华北板块地区深度约为2～6km处,沿华北板块边缘和郯庐大断裂带,分布比较集中；北苏鲁拆离断层以北地区地下深度约为2～4km,沿沉积盆地边缘和北苏鲁超高压构造片体边缘,分布比较松散；苏鲁超高压构造片体深度约为2～6km处,沿苏鲁高压构造片体边缘,磁性球体分布比较集中；沭阳盆地地区没有明显特征；为了详细了解研究区地下结构分布,结合地表出露地震、地震反射、大地电磁、CCSD岩石物性资料,对研究区重磁资料以剖面反演结果详细描述研究区密度、磁性空间分布特征。穿过泊里到诸城(南苏鲁超高压地区)和CCSD主孔(北苏鲁超高压地区)的两条地震反射剖面分别长约44km和139.5km。参考地震反射数据的分析结果,上地壳至上中地壳结构在我们模型中厚度小于15km,反演模型曲线很好地拟合了观测曲线,模型包含有由沉淀层、花岗岩、各种片麻岩、榴辉岩和蛇纹石化橄榄岩等以透镜体、块体、地层形式组成的上地壳-上中地壳；蛇纹石化弧前地幔位于中下地壳。它们的密度从2.50到3.45gcm-3(上地幔)范围变化,磁化强度从0到4.55A/m范围变化。结合地质背景对模型进行解释分析,推测苏鲁地区深部地壳处可能存在蛇纹石化弧前地幔,为理解大别-苏鲁地区的超高压变质过程提供新的参考见解。