应用地球物理2017

重力勘探重力=地球引力+天体引力+离心力概念：通过观测与研究天然重力场的变化规律以查明地质构造和寻找矿产的一种 物探方法。基础：地壳中不同岩、矿石间密度差异物理学基础：牛顿万有引力定律1890年，匈牙利物理学家厄缶（18481919年）制造出了第一台测量重力变化率的扭秤。大地水准面的三级近似第一级近似为正球体，其平均半径为6376km ；第二级近似为旋转椭球体，赤道半径为6378km，极半径为6356km ；第三级近似为梨形体，其北极高出十余米， 南极凹进二十余米。地球表面的岩石密度（1）决定岩石密度大小的因素矿物成分及含量孔隙度及孔隙中的含水量埋藏深度（2）岩石密度大小的一般规律火成岩 变质岩 沉积岩火成岩的密度矿物成分及含量（主要）：随高密度铁镁暗色矿物含量的增加而增大，酸性岩V基性岩V超基性不同的成岩环境（侵入或喷发）：侵入岩 喷出岩沉积岩的密度 孔隙度（主要）：灰岩、页岩、砂岩等孔隙度高达30%40% ,干燥的岩石随 孔隙度减少密度呈线性增大充填物成分及充填孔隙所占比例：水 油 气 干燥成岩时代及埋深：上覆岩层具有压实作用，密度随着埋藏深度和时间而增加变质岩的密度原岩密度（主要）变质程度：变质作用有助于岩石孔隙的充填，并使岩石以更致密的形式再结晶， 其密度往往随变质程度的加深而增加。影响重力分布特征的地球内部因素地壳厚度的变化（莫霍面起伏）上地幔密度的不均匀性深大断裂的可能部位及延伸情况重力场强度的单位CGS 单位制 1cm/s2=1GalSI 单位制1m/s2=106g.u.质点（质量为m）的引力位：v=G*m/r重力位W等于引力位与离心力位之和重力等位面的特点：等位面上各点重力的方向与等位面的切面相垂直等位面上各点的重力等于重力位对该点等位面的内法线n的方向导数重力等位面不处处平行，既不相交也不相切，而是一个单值函数重力的方向：等位面的垂线方向 重力的大小与重力等位面之间距离 l成反比重力（重力加速度）在空间上变化的原因1）地球形状的影响：地球是近似于两极压缩的扁球体，地表起伏不平，引起约6000毫伽的重力变化。2）地球的自转：地球绕定轴旋转，惯性离心力使重力有3400毫伽的变化。3）地壳及上地幔密度分布的不均匀性：可引起几百毫伽的重力变化。正常力公式分析根据研究的需要而确定，不同学者计算出的正常重力值有所差异，不是客观上存 在的确切的。正常重力场只与计算点的纬度有关，沿经度方向无变化在赤道处最小，在两极处数值最大，相差约5伽沿经度方向的变化率与纬度有关，在45纬度处变化率最大力异常的概念、计算公式及重力异常的物理意义1）重力异常：由地下岩矿石密度分布不均匀所引起的重力变化，或地质体与围 岩密度的差异引起的重力变化。2）计算：将实测重力值减去该点正常重力值，其差值称为重力异常3）物理意义：剩余质量所产生的引力在重力方向的分量（或在铅直方向的分量）获得探测对象产生的重力异常，一般应具备的条件必须有密度不均匀体存在，但仅有密度不均匀体的分布，并不一定能产生重力异 常不仅探测对象与围岩要有一定的密度差，而且剩余质量不能太小探测对象不能埋藏过深能否取得探测对象产生的异常，还取决于该异常能否从干扰场中辨别出来重力测量原理绝对重力测量测量重力全值：自由落体运动数学摆相对重力测量测量两点间重力差值ZSM石英弹簧重力仪主要结构：弹性系统光学系统测读机构弹性系统用于感受重力变化，是仪器的心脏工作任务：预查普查详查细测零点漂移：在相对重力测量中，由于重力仪灵敏系统的弹性疲劳、温度补偿不完全等因素，仪器读数的零点值随时间而不断变化混合零点漂移：零点漂移值叠加上重力日变的影响零点漂移校正：仪器每次从基点出发，观测一段时间之后，需回到原基点（另一 基点或总基点）上观测一次，测出零点漂移量；在较短的时间内，可认为零点漂 移是随时间线性变化的k 二重力测量值的主要改正项及其计算方法i） 正常场（纬度）校正：消除测点与基点间纬度差异导致的重力变化。测点与基点间正常重力的变化量Ag2 ） 地形校正：消除测点附近地形起伏对重力观测结果数据的影响。盈余或亏损质量对测点A处引力的垂向分量3）中间层校正：消除测点基准面与基点基准面之间水平中间层的重力影响。4）高度校正：消除测点相对于基点的高程差而造成的重力数值变化。正常重力随高度变化的变化量Agh布格里力异常高度校正和中间层校正都与测点的高程有关系，通常把这两项合并起来称为布格校正观测重力差值经过正常场校正、地形校正和布格校正之后得到异常称为布格重 力异常.意义：消除了观测面以下正常密度分布的重力作用，对 异常（剩余）密度分布 的重力作用未作调布格力异常的复杂性 叠加异常 异常的组成:中、深部地质因素的重力效应（主要，分布范围较广）地壳厚度的变化（最主要）结晶基岩内部成分、构造和基底顶面的起伏浅层和局部剩余密度体的重力效应（小部分，分布范围较小）沉积岩的成分和构造（一般10毫伽以下，甚至仅十分之几毫伽金属矿及其）它矿产（一般十分之几毫伽，个别达几毫伽） 信息量巨大不同地质因素引起的异常在幅度大小、分布范围、变化快慢等特征方面都不尽 相同对实测重力值只做纬度（正常场）校正与高度校正物理意义：作正常场校正相当于从观测重力值中消除密度为正 常分布的大地椭 球体的正常重力值；只作高度改正，相当于把大地水准面以上物质都“压缩”到 大地水准面上，且未改变大地水准面以下的质量分布应用：在地形高差不太大（小于2000m）的情况下，可用于研究上地幔与地壳的物质分布情况（均衡分析）均衡力异常均衡校正值：计算“填补”的物质在测点处所产生引力的铅垂分量均衡重力异常：将均衡校正值加到布格异常中的结果大的正异常值（“补偿”不足）：物质减少或由内部物质相的变化使地表下沉（莫 霍面下移）大的负异常值（“补偿”过剩）：物质补充或是由内部物质相的变化使地表上升（莫 霍面上移）力异常的划分区域校正法：用绘图或异常取平均的办法去掉局部异常得到区域异常异常平滑曲线法、圆周平均法及网格法等异常变换法：将重力经过数学变换，换算成另外的异常形式，进而压制局部异常、突出区域异常 或压制区域异常、突出局部异常重力位高阶导数法、异常向上延拓法等力异常高阶导数的特点不同形状地质体的重力异常导数具有不同的特征，有助于对异常的解释和分类 突出浅而小地质体的异常特征，压制区域性深部地质因素的重力效应，分离不同深度和大小异常源引起的叠加异常 可分离互相靠近、埋藏深度相差不大的相邻地质体引起的叠加异常导数的次数越高分辨能力越强正演：给定地球物理模型，通过数值计算或物理模拟，得出相应的地球物理场 反演：反问题就是根据已知的地球物理场，求解可能的地球物理解释 三度体：各方向上都为有限量二度体：某一方向无限延伸，而在该方向上埋藏深度、截面形状、截面大小和 物性特点皆稳定不变密度均匀球体的重力异常平面等值线分布特征在原点处，异常取得最大值；当XT8时，异常趋近于零-分母中均为二次项，故异常相对原点对称分布异常的平面等值线为一系列的同心圆，圆心位于球心在地面的投影点处密度均匀水平IBI柱体的重力异常最大值x=0处异常取得最大值1/n极大值坐标Xi/n=D (n-1) 1/2平面等值线的分布特征：x=0时，异常取得最大值；当xts时，异常趋近于零。 x1/n为定值，异常等值线为一簇平行不等间距的直线，并以柱体中轴线在地面 的投影线为对称轴对称分布。铅垂台阶在重力异常平面图上，等值线是一系列平行于台阶走向的直线在断面两侧形成的异常变化梯度呈对称分布的等值线密集带，称为重力梯级带， 是识别断裂构造的重要标志只要台阶厚度H-h不变，不论台阶的上顶埋藏深度为何值，其gmin、Ag(0) 和gmax均不变，只是整条曲线随埋藏深度的增大而变缓根据曲线不对称的情形可以判断倾斜台阶的倾向：在倾斜方向的一侧Ag曲线较缓反演方法：直接法特征点法选择法人机交互式反演方法引起多解性的原因：场的等效性：不同的地球物理模型可以对应相同的地球物理场观测数据的有限性和离散性：实际的测量都是在有限的范围内在离散的测点上进行，有限的、离散的数据减少了地球物理场所包含的信息误差影响：实测数据都是包含误差的，使得相差不超过误差限的不同场难以区分；而且数据处理过程可能会带来误差限制多解性的方法多种方法综合应用提高观测精度提高处理精度重力异常变化总趋势从东向西逐渐降低区域性重力异常背景上存在几组巨大重力梯阶带大兴安岭一太行山一武陵山重力梯级带青藏高原周边重力梯级带钓鱼岛重力梯级带巨大重力梯阶带与我国主要褶皱山系平行或合在莫霍面深度图上表现为陡坡带大洋中脊处的重力异常特征布格重力异常：存在大约200毫伽的重力低值自由空间重力异常：比邻近海底高2030毫伽洋中脊下存在低密度的物质异常地幔海淘处的重力异常特征5公里异常的最小值位于海沟内壁上方，与自由空间重力异常相比，其最小值位置更靠近陆侧。这表明海沟内靠陆一侧下有物质亏损陆洋分界带的重力异常特征陆壳一侧为正异常，而洋壳一侧为负异常改变洋陆交界以上莫霍界面斜率，重力异常的斜率随之改变莫霍界面的斜率越大，重力异常的变化越剧烈陆壳上沉积层分布、厚度等因素也会直接影响到陆缘上的重力异常曲线的形态和幅度大小（沉积物的密度变化所引起）磁法勘探物质基础：地壳中各种岩、矿石间的磁性差异岩、矿石间的磁性差异将引起正常地磁场的变化（即磁异常），通过观测和研究磁异常来寻找有用矿产或查明地下地质构造磁法勘探的分类：地面磁测 井中磁测 航空磁测海洋磁测与重力勘探相比，磁法勘探的特点相对幅值：磁异常比重力异常大得多反映的地质因素：磁异常较为单一复杂程度：磁异常要比重力异常复杂1870年，泰朗和铁贝尔制成了寻找磁铁矿用的万能磁力仪，有人认为这是应用地球物理学开始发展的一个重要标志 磁场强度：T=1/4nu*q*Q/r2 磁位 U=1/4nu*q*Q/r磁偶极矩p ,即磁量Qm与两磁极之间距离l的乘积磁矩为:m=q/u地磁要素建立一个空间直角坐标系统：设以观测点（O）为坐标原点，x轴正向指地理北， y轴正向指东，z轴正向垂直向下。地磁场T分解为：北向分量X ；东向分量为Y ；垂直分量为Z。T在xoy面内的投影为水平分量H，H的方向即磁北方向；H与x的夹角为磁 偏角D，即磁北与地理北的夹角，东偏为正；T与H的夹角为磁倾角I，下倾为 正。X、Y、Z，H、D、I，T统称为地磁要素。磁感应强度B = uTT0（中心偶极子磁场）和Tm（大陆磁场或世界异常）之和称为地球基本磁场，编制 的世界地磁图大多为地球基本磁场的分布图T0几乎占80%85%，代表地磁场空间分布的主要特征地心偶极子磁场T0:场源是地核中物质的对流运动所形成的涡流大陆磁场Tm:地球内部构造不均匀性，在地核和地幔边缘部分存在着未知对流, 形成涡电流长期变化场基本地磁场0+随时间所发生的缓慢变化，与地球内部因素有关。第一个特征为地磁极的西向漂移。第二个特征是不同年代计算出地球磁矩的衰减变化短期变化场太阳静日变化-地磁日变变化特征-周而复始，从不间断，振幅易变，但相位不变，周期24小时-变化幅度白天大(618点)、夜间小，夏季大、冬季小-变化形态和幅度因地而异，但同纬度区变化相近30纬度区变化最大磁暴-强烈磁扰动，幅度在nnxl03nT之间-从赤道到极地均可观察到，而且几乎是全球同时发生，影响的时间为几小时 几天-对地磁场水平分量的影响特别显著，而对垂直分量影响相对较小-高纬度地区发生多、小磁暴多、春秋季多，磁暴强度从低磁纬度到高磁纬度逐渐变强监测到磁暴时，应立刻停止野外磁测工作地球基本磁场的分布特征D等偏线图曲线族分别汇聚于南、北两磁极区，这两点D在0360间变化，说明磁北的方向在磁极处无固定方向，即没有固定的磁偏角，南北两磁极为源头全球有两条零偏线(D= 0)分布，将全球分为正、负两个偏角区，即在正偏角区磁北方向为地理北偏东，而在负偏角区为地理北偏西I等倾线图等倾线分布大致和纬度平行零倾线（1=0）在地理赤道附近磁赤道，说明地磁场在赤道上为 水平方向， 垂直分量为0 （ T=H , Z=0 ），磁赤道并非直线 随纬度增大I增大，至南北磁极处I=90 ,H=0,T=Z磁赤道以北I 0 , T 下倾；磁赤道以南Iv0。，T上倾H等值线图H等值线大致沿纬度线排列随纬度增大H值减小，在赤道附近H最大，Hmax=X=40000nT ;在磁极处Hmin=0除两磁极外，全球各点的H均指向北Z等值线图Z与I等倾线分布相似，近乎与纬度线平行纬度增大Z绝对值变大，两磁极处值最大，Zmax=6000070000nT;磁赤道处Z=0，磁赤道以北Z0，方向垂直向下，反之亦然T等值线图大部分地区的T等值线也近乎与纬度线平行T值由赤道向两极逐渐增大赤道处有极小值：Tmin=H=3000040000nT，磁极处有极大值：Tmax=Z=6000070000nT正常场和磁异常一般是相对的，通常把研究对象引起的磁场部分叫做磁异常，而 周围环境和围岩引起的磁场统归为正常场磁异常是一个矢量，实际工作中需计算（测量）Ta的分量，包括垂直分量Za、水平分量Ha以及总强度磁异常AT4近似为在正常场0方向的投影 磁化强度M:单位体积内分子电流磁矩矢量和，表征在外磁场中物质可被磁化的程度磁化率K为由物质本身性质所决定的参数，表征物质被磁化的难易程度磁感应强度为场源在观测点的磁场强度与磁化物体所形成附加磁场强度之和 铁磁性物质内电子之间因自旋引起的相互作用非常强烈，在铁磁质内部形成了一 些微小磁化单元（区域）一磁畴-磁畴中，各电子的自旋磁矩排列得很整齐。 当外磁场很强时，磁畴的磁化强度方向转向外磁场方向，最后达到饱和磁化 当温度升高时铁磁性物质磁化率逐渐增加，临近某温度时达到极大值；然后急 剧下降，趋于零居里点为铁磁性物质的磁化强度陡然降低、物质由铁磁性转为顺磁性的温度 总磁化强度M=感应磁化强度+剩余磁化强度感应磁化强度：所有岩石在地 磁场作用下被磁化剩余磁化强度：在历史地磁场中获得的磁化强度火成岩具有明显的天然剩余磁性由酸性T中性T基性T超基性，磁性由弱T强同一成分的火成岩其磁性不同，喷出岩磁性 侵入岩磁性变质岩磁性与变质前的源岩和变质条件有关通常正变质岩磁性负变质岩磁性层状结构，磁性表现为各向异性沉积岩k及Mr都很小，磁性很弱，通常认为无磁性。磁化率主要取决于铁磁类副矿物的含量与成分影响岩矿石磁性的主要因素及岩矿石磁性的一般规律。1）铁磁性矿物含量含量越高，岩石磁性越强，但并不呈简单的线性关系。2）铁磁性矿物颗粒大小及结构-当铁磁性矿物含量一定时，颗粒越大，磁性越强-当磁性矿物颗粒大小、含量都相同时，颗粒相互呈胶结状者比颗粒呈分散状者磁性强。3）岩石磁性与温度、压力的关系-深度增大，地温升高，压力增大，磁性减弱。-应力作用使岩石沿应力方向磁性减小，如断裂、破碎带上磁性减弱；变质、蚀 变作用往往使岩石磁性增强。-岩石剩余磁化强度与压力关系明显，当压力为100MPa时，剩磁可减小25%。-随着深度增加，地壳内岩石的温度升高，岩石圈静压力增大，二者对岩石的磁性有很大影响，因此，在研究地壳深部岩石磁性时，需要考虑温度和压力的影响。回答剩余磁性的主要类型及其成因剩余磁性主要可分为原生剩余磁性和次生剩余磁性两大类。1）原生剩余磁性。可具体再分为热剩余磁性、碎屑剩余磁性、化学剩余磁性三 类。 热剩余磁性：岩石在冷却的过程中，受当时恒定的地磁场作用，磁化所获得 的剩磁。 碎屑剩余磁性：沉积物固结成岩后，按其碎屑的磁化方向保存下来的磁性。颗粒越细，碎屑剩磁方向与当时地磁场方向越一致 化学剩余磁性：在一定磁场中，某些磁性物质在低于居里温度的条件下，经过相变 过程（重结晶）时化学变化（氧化还原）所 获得的剩磁。2）次生剩余磁性。可具体再分为粘滞剩余磁性、等温剩余磁性两类。 粘滞剩余磁性：岩石形成后，长期处在地球磁场作用下，随着时间的推移，其 中原定向排列的磁畴逐渐地弛豫到作用磁场的方向，所形成的剩磁。 等温剩余磁性：在常温条件下，岩石因受外部磁场的作用（比如闪电作用）， 而使近地表岩矿石磁性发生大小和方向的改变而获得的剩余磁性。质子磁力仪20世纪50年代中期问世，广泛用于航空、海洋及地面等磁测具有高灵敏度、高准确度的特点测量地磁场总强度T的绝对值、梯度值由人工场方向逐渐趋向于地磁场方向时，质子围绕T旋进并趋向于T的方向，这种旋进称拉莫尔旋进决定磁性体磁化强度的因素地磁场的大小和方向（外因）磁性体本身有关因素（内因）磁性体的形状、埋深及其内部各点磁化率性体的走向，因在一个局部地区，地磁场的方向是一定的；而磁性体走向的差异 可导致其表面磁荷分布的不同顺轴磁化、无限延伸的柱体等（见pdf） 磁异常反演中存在非均匀磁性问题和反演的多解性问 特征点法经验切线法沃纳反褶积方法欧拉齐次方程方法3D可视化反演火成岩、沉积岩及变质岩的磁性特征。1）火山喷出岩 共同的特点是异常呈跳跃变化，尖锐而且梯度大，相邻测线难以对比。狭窄的 磁力高或低可能是火山喷出裂隙的反映；单个极强的峰值可能是火山口的反映。 另一个特点是异常强度随高度的增加而迅速衰减。2）火山侵入岩 侵入岩异常峰值可能很大，但形态比较圆滑相邻测线上的异常曲线可以对比。 异常形态与埋深之间有明显的依赖关系，埋藏较浅时常表现为多个孤峰，埋藏 较深时只显示异常不规则的背景。3）沉积岩 沉积岩的磁场特征一般是异常值微弱，形态平缓、光滑、梯度小，沉积厚度较 大时呈现为平静的负异常区。 主要是由于沉积岩的磁化率小，而且磁性矿物在缓慢的沉积过程中分布较为均 匀，无方向性。3）变质岩 异常的强弱与变质母岩的磁性和变质类型有关。 一般正变质岩的磁性较强，负变质岩的磁性则相对弱得多。海洋磁异常:磁异常呈条带状分布，条带的走向与洋脊平行正、负异常相间，正、负异常带宽2030km，长几百千米，异常幅值几百纳特磁异常对称于洋脊分布电法勘探地壳中不同岩层之间、岩石和矿石之间存在的电磁性质差异通过观测天然存在的或由人工建立的电场、电磁 场分布，来研究地质构造、寻 找有用矿产资源，解决工程、环境、灾害等地质问题电法勘探的特点三多：可利用的物性参数多利用的场源多方法种类多两广：应用空间广 应用范围广电阻率P是表征物质导电性的基本参数。电流垂直流过单位长度、单位截面积的体积时，该体积物质所具有的电阻值电导率。表示物质导电性的优劣，为电阻率的倒数若进行测量的地段地下岩石电性分布不均匀时，测量（或计算）的结果称之为 视电阻率。虽然不是岩石的真电阻率，但却是地下电性不均匀体和地形起伏的 种综合反映在电法勘探中，常把按电阻率划分的地质断面称为地电断面电阻率法以地壳中岩石和矿石的导电性差异为基础通过观测与研究人工建立的地下电流场（稳定场或交变场）的分布规律进行找矿 和解决地质问题的勘探方法影响岩、矿石电阻率大小的主要因素。1 ）导电矿物含量。导电矿物含量相当大（280 %）时，才对岩石电阻率有明显作。2）岩、矿石的结构、构造、孔隙度。导电矿物对岩、矿石整体电阻率影响主要决定于它们的连通情况：连通者起的作 用大，孤立者起的作用小。3）岩矿石的含水量及含水矿化度。地下水及其他天然水的电阻率均较低，并且含盐分越多，电阻率值越低。含水量 大的岩石电阻率较低，而含水量小或干燥岩石的电阻率较高。4）温度、压力。电子导电矿物（或矿石）的电阻率随温度增高而上升；离子导电岩石的电阻率随 温度增高而降低。温度的变化将引起水溶液中离子活动性的变化，所以岩石中水 溶液的电阻率也将随温度的升高而降低。冰冻条件下，地下岩石中的水溶液将由 于结冻，使岩土呈现出极高的电阻率。影响岩、矿石电阻率的因素是多方面的： 在金属矿产普查和勘探中，岩、矿石中良导电矿物的含量及结构是主要影响 因素； 在水文、工程地质调查和沉积区构造普查、勘探中，岩石的孔隙度、含水饱和 度及矿化度等成了决定性因素； 在地热研究、地震地质及深部地质构造研究中，温度及地应力变化却成了应考 虑的主要因素； 欧姆定律微分形式的推导。均匀介质中异性点电源的稳定电流场分析