普通物探第2-2节磁法勘探的正反问题.ppt

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,*,（华东）,2.2,磁异常正、反演,为了对磁异常做出合理的推断解释，必须要清楚异常特征与磁性体的关系，通常借助于数学物理的方法进行研究。,为便于数学手段的应用，需要对复杂的自然条件做适当的简化，以便于正反问题的开展。,一方面，了解简单是认识复杂的基础；另一方面，规则和均匀的概念都是相对的，因条件而变化。,磁异常问题的简化条件,通常所做的简化有：,磁性体形状规则；,磁性（磁化率）均匀；,不考虑剩磁或剩磁与感应磁性方向一致；,地质体的磁化磁场,-,地磁场是均匀的；,只考虑单个磁性异常体；,观测面水平。,磁异常问题的简化条件,除少数情况外，实际地质条件并不符合上述假设条件，理论上讲，只有二次曲面形体才能被均匀磁化。,由于磁法主要是研究被覆盖的地质体，形态不规则和磁化不均匀所造成的影响很小，这些假设并不失去实际意义。,主要内容,本节主要介绍计算磁性体磁场的基本公式、几种简单形体的磁异常，以及根据磁异常求取磁性体埋深的简单方法。,主要内容：,2.2.1,计算磁性体磁场的基本公式,2.2.2,简单形体的磁异常正反问题,2.2.3,磁性体埋藏深度的计算,2.2.1,计算磁性体磁场的基本公式,对一个已知地质体，磁异常的计算比重力异常要复杂得多。磁异常不但取决于物体的几何形状和磁化率，而且还与磁化方向密切相关。,目前，计算磁性题磁场的方法有多种，主要是,积分公式法,和,泊松公式法,。,磁异常计算式,根据磁场与磁位的关系：,磁异常的三分量计算式分别为：,计算磁位的面积分公式,正演计算时经常遇到由几个平面围城的形体，如各种,板状体,等，这种情况下，用磁荷面积分更方便。,根据计算磁性体磁位的基本公式：,再考虑到：,因此有：,对计算点求导数,对源点求导数,计算磁位的面积分公式,又因为：,所以：,代入前式，得到：,计算磁位的面积分公式,根据矢量积分中的高斯公式：,其中,n,为磁性体表面的外法线方向，因此：,这是场论的一个重要结论，它表明磁性体在,P,点的磁位，等于磁荷面密度 的各个表面和磁荷体密度 的磁体所引起的磁位之和。,计算磁位的面积分公式,当磁性体为均匀磁化时，,M,为常量，其散度为零，即内部磁荷体密度为零，只在其表面有磁荷分布。,此时，磁性体外部的磁位可以通过其面磁荷积分来计算，即,面积分法,，公式为：,对于由几个平面围城的磁性体，每个面上的面磁荷密度是常量，因此计算相对简单方便。,磁异常计算式,根据磁场与磁位的关系，磁场各分量的积分公式分别为：,当形体和磁化强度给定后，可以根据上述面积分计算磁性体的磁场。,2.,计算磁位的泊松公式,根据物体的重力位公式,得出,从而求出磁性体的磁位：,这就是,根据重力位计算磁位的泊松公式,。,根据场与位之间的关系，可以求出磁场矢量，进而求出各个分量。,磁异常计算式,利用重磁位场的泊松公式，磁场各分量的计算式分别为：,对于沿走向无限延伸的二度体，引力位函数与,y,无关，上式会大为简化。,2.2.2,简单形体的磁异常正反问题,对从实际地质模型简化出的规则磁性体的磁场进行数学物理的解析，从中找出异常磁场与磁性体之间的关系，作为异常数据解释推断的依据。,主要内容,点磁极的磁场,球形体的磁场,线磁极的磁场,面磁极的磁场,1.,点磁极的磁场,一个点磁极是最简单的磁体模型，用于模拟顺轴磁化，向下无限延伸的细长柱体。此时只在柱体顶端显示负磁荷，可用单点磁极的模型模拟。,根据磁库仑定律，该异常体引起的磁异常：,式中,q,m,为柱体顶端的磁荷量。,直角坐标系中，磁异常垂直分量为：,J,x,y,z,O,h,P(x,y,0),T,a,z,a,单点磁极模型,单点磁极的磁场,平面异常等值线为一系列同心圆，在中心点取极大值；剖面异常曲线呈现左右对称性质。,在过中心点的主剖面上，异常极大值：,单点磁极的磁场,单点磁极的磁场,主剖面异常半极值点横坐标,x,1/2,与磁极埋深,h,之间存在如下关系式：,这就是解反问题的特征点法。,单点磁极的磁场,双点磁极的磁场,当顺轴磁化的细长柱体长度有限时，在柱体的两端均显示磁荷聚集，磁量相等，符号相反，这可用双点磁极的模型来模拟。,此时，总的磁场是正、负磁荷磁场的叠加，即：,双点磁极模型,J,x,y,z,O,h,P(x,y,0),T,a,z,a,l,l,双点磁极的磁场,设柱体长度为,2,l,，倾角为,，其中点在地面的投影为坐标原点，正、负磁极各自产生的磁场分别是：,J,x,y,z,O,h,P(x,y,0),T,a,z,a,双点磁极模型,l,l,双点磁极的磁场,在北半球，正异常在磁性体南侧，异常幅值大，只在有限范围内分布；负异常在磁性体北侧，幅度小，分布范围要大得多。,双点磁极的磁场,双点磁极的磁场,由于磁场垂直分量向下为正，向上为负。,在北半球，正异常在磁性体南侧，异常幅值大，只在有限范围内分布；负异常在磁性体北侧，幅度小，分布范围要大得多。,双点磁极异常剖面图,在,-90,磁化情况下基本为简单的负异常；在,-45,情况下北侧负异常明显；在水平磁化时正负异常相当；在,45,时南侧正异常明显；在,90,磁化时基本为简单的正异常。,不同磁化倾角情况下的双点磁极磁场剖面图,双点磁极异常平面图,不同磁化倾角情况下的双点磁极磁场平面图,-90,-45,0,45,90,2.,球形体的磁场,对于一些尺度大小有限的地质体，当埋深远大于直径时，其地面磁场与球体磁场特征近似。,设球体的磁化强度为,M,，体积为,V,，那么，球体的磁矩,m,=,M,V,，根据球体的重力位 和重磁位场的泊松公式，相应的磁位计算式为：,是磁矩,M,与矢径,r,的夹角。,可见，,球体外的磁位与磁偶极子的磁位具有相同的形式,。,球形体的磁场,对,z,求偏导数，再令,z,=0,，得到球体在地面上引起的垂直磁异常：,式中，,I,为磁化倾角，,为方位角（测线与,x,轴的夹角）。,球形体的磁场,当测线通过球心在地面的投影，并且与磁化方向在同一铅垂面时：,当测线通过球心在地面的投影，并且与磁化方向垂直时：,南北向剖面异常,倾斜磁化球体的磁场,异常平面等值线图,东西向剖面异常,E,Z,a,球形体的磁场,球体磁场与双点磁极磁场的特征近似。,不同磁化倾角的球体磁场,倾斜磁化球体的磁场,-90,-45,0,45,90,当向下无限延伸的薄层受到顺层极化时，层的侧面不显示磁荷，只在顶端显示磁荷，这样的地质问题可用一个线磁极的模型模拟。,假设薄板顶面埋深为,h,，水平宽度为,2,b,h,，走向延伸长度为,2,L,，磁极化强度,J,方向与板的下倾方向一致，倾角为,i,，取顶面中心点在地面的投影为坐标原点，,y,轴与走向一致。,3.,线磁极的磁场,单线磁极模型,J,x,y,z,O,h,P(x,y,0),T,a,z,a,i,空间点,p,(,x,y,z,),的垂直磁异常为：,式中 表示,磁荷线密度,。,单线磁极的磁场,J,x,y,z,O,h,P(x,y,0),T,a,z,a,i,单线磁极的磁场,当线磁极无限延伸时，垂直磁异常为：,垂直于磁极线的地面主剖面上：,异常剖面曲线和平面等值线与水平圆柱体的重力异常形态特征一致。,剖面异常曲线,x,Z,a,x,y,平面异常等值线,单线磁极的磁场,双线磁极的磁场,当顺层磁化的薄层向下延伸有限宽度,2,l,时，在层的顶底面均有磁荷显示，相当于双磁极线，两极分别产生磁场，总的磁场是二者的叠加。,为简单起见，这里只给出走向无限延伸时地面上的垂直磁异常：,式中,h,为双极中心到地面的埋深，,i,为倾角。,J,x,y,z,O,h,P(x,y,0),T,a,z,a,i,双线磁极模型,双线磁极的磁场,双线磁极的磁场,Z,a,剖面曲线左右不再对称，埋藏浅的一极异常幅值大，范围小，埋藏深的一极异常幅值小，范围大。,当薄层的倾向延伸宽度,2,l,h,时，层的宽度可以忽略，此时可用偶极线模型模拟，等价于水平圆柱体模型。,双线磁极的磁异常,4.,面磁极的磁场,在顺层极化的异常体厚度较大的情况下，可用有限宽度的面磁极模型模拟。,选择坐标系,y,轴与走向一致，原点位于中心线在地面的投影，为简单起见，假设异常体沿走向无限延伸，此时主测线上一点,P,(,x,),的垂直磁异常：,式中 为磁荷面密度，是岩层顶面对计算点所张开的角。,面磁极的磁场剖面异常,半无限面磁极,岩层接触面的情况可用半无限面磁极的模型模拟。,主测线上的垂直磁异常：,无限面磁荷的情况下,x,z,Z,a,半无限面磁极的磁场剖面异常,垂直磁化岩层,对于一个厚度有限的垂直磁化岩层，顶、底面都有磁荷分布，总的异常是顶、底面磁荷分别引起的磁场的叠加。,对于半无限水平层：,其中,h,1,、,h,2,为岩层的顶、底埋深。,无限水平层的磁异常：,x,z,Z,a,半无限面磁极的磁场剖面异常,2.2.3,磁性体埋藏深度的计算,在利用磁异常资料研究含油气盆地结晶基底的起伏时，一般不能将基底表面看成一个简单的磁性分界面，像重力勘探方法一样求解，因为：,由于结晶基底的磁性通常是极不均匀的，一般不能把其表面看成一个均匀的磁性分界面。,磁异常,Z,a,和,T,通常与磁性界面的起伏之间也不存在简单的直接对应关系，而是随具体条件不同呈现复杂的关系；,磁性体埋藏深度的计算,结晶基底大多经历过强烈的褶皱和变质作用，存在大量的岩浆岩，在盖层沉积之前已遭受了长期的风化剥蚀，被夷为准平原，基底内很多孤立的、具有陡立产状的磁性岩体顶面被切平了，直接出露于基底表面，构成了基底表面的一部分，研究这些孤立磁性体顶面的深度便可了解结晶基底表面的深度变化。,求解磁性体埋深不是利用异常值的大小，而是利用其变化形态。具体方法有切线法、特征点法、选择法、积分法等。,特征点法,特征点法,是根据磁异常曲线上的某些特征点，如极大值点、极小值点、半极值点、零值点及拐点等，利用其坐标位置及其间的距离来求取地质体位置和产状的方法。,这些特征点的坐标分别用,x,max,、,x,min,、,x,1/2,和,x,1/2,、,x,0,和,x,0,、,x,g,和,x,g,等来表示。,利用特征点法计算磁性体的埋深，不同形体有不同的公式。,切线法图解,特征点法,:,球体,对于球体异常，其宽度不仅与埋深有关，还与有效磁化倾角,i,s,有关，因此，特征点的坐标和数值随有效磁化倾角,i,s,而变化。,中心埋深的近似计算式为：,切线法图解,特征点法,:,水平圆柱体,水平柱体异常曲线的特征点坐标和数值也随有效磁化倾角,i,s,而变化。,中心埋深的近似计算式为：,在,Z,a,曲线接近反对称的情况下，可用极值点间距求埋深：,切线法图解,特征点法,:,厚板状体,顺层极化、向下无限延伸的厚板状磁性体的,Z,a,曲线对称无负值，其上顶埋深的近似计算式为：,斜交极化时：,计算 在两极值间的交点与板状体中心对应；,计算辅助曲线 ；,利用,f,(,x,i,),曲线的,x,1/2,和,x,1/4,特征点计算埋深,h,。,斜交磁化厚板埋深的求法,经验切线法,经验切线法,是通过做异常曲线的切线，根据切线交点横坐标求取磁性体埋深的方法。,在,Z,a,或,T,曲线的极大值点作一条与横轴平行的切线，在两侧的极小值点分别作横轴的平行线（无极小点时作渐近线）；,在曲线两侧的拐点各作一条切线与水平切线相交；,切线法图解,经验切线法,曲线左支或右支拐点的切线与极大点和极小点切线（或渐近线）的两个交点的水平距离就大约等于磁性体顶面埋深的两倍：,实测曲线往往是不对称的，取其两侧的平均值作为磁性体的深度，即：,切线法图解,切线法校正系数,实践表明，切线法仅对顺轴磁化且走向和延伸大的地质体异常体,能取得较好的效果。对其他情况则误差较大，需要乘一个校正系数。,磁体形状,条件,校正系数,无限延伸薄板,-,i,=90,-,i,=45,-,i,=0,0.41,0.56,0.65,水平圆柱体,i,=30,i,=45,i,=90,1.03,1.10,1.30,直立柱体,顺轴磁化,0.86,球 体,i,=40-90,i,=30,i,=150,1.30,1.40,2.30,切线法求磁性体埋深的校正系数,Depth Determination,Can get very approximate depth from magnetic anomalies,Sphere or half-cylinder:,Depth to centre,of body w is roughly equal to width of anomaly peak at half its maximum value dFmax/2.,Dipping Sheet or Prism:Depth to centre of body is roughly width of linear segment of anomaly d.,半极值宽度法,Draw tangent at point of maximum slope(line 1),Find two tangents to curve with half maximum slope(lines 3,4),Depth to top of body is distance between tangent points,Peters Half-Slope Method,(theoretically-based),Map of Half-Slope Method,思考题,作业题,已知顺轴极化细长柱体长,200m,，中心点坐标（,0,0,300,），计算不同倾角情况下的垂直磁异常，绘制南北向主剖面异常曲线及平面等值线图。计算范围,x,：,-10001000,，,y,：,-10001000,。,（单号）,已知均匀磁化球体中心坐标（,0,0,300,），计算不同磁化倾角（,-90,、,-45,、,0,、,45,、,90,）情况下均匀磁化球体的垂直磁异常，绘制南北向主剖面异常曲线及平面等值线图。计算范围,x,：,-10001000,，,y,：,-10001000,。,（双号）,