利用重、磁异常可以研究地质体的横向不均匀性,特别是其边缘位置.文中的地质体边缘是指断裂构造线、不同地质体边界线等,实际上是具有一定密度或磁性差异的边缘位置.由于在该边缘位置存在一定的密度或磁性差异,使得地质体边缘附近的重、磁异常变化率较大,故可以利用重、磁异常的变化率来识别地质体的边缘位置.利用重、磁异常识别地质体的边缘位置有数理统计,数值计算和其它三大类[1].数值计算类边缘识别方法是研究最多、应用最广的边缘识别方法,此类方法有垂向导数[2~5]、总水平导数[6~9, 5]、解析信号振幅[10~14]三种基本方法和倾斜角[15]、θ图[16]两种基本比值方法,其他方法均是在这些方法的基础上发展而来.数值计算类边缘识别方法的理论基础是二度体铅垂台阶模型边缘处重、磁异常的变化率达到极大值或零值,故可以利用这一特点来准确识别二度体铅垂台阶的边缘位置.对于倾斜二度体、不规则二度体以及三度体边缘位置的识别均是二度体铅垂台阶模型理论的推广,但识别的边缘位置与真实位置有一定偏差(偏移).该偏差随着地质体边界形状、埋深、水平尺寸及物性差异等的变化而变化.
总水平导数最初由Cordell[6]于1979年提出,并利用重力异常总水平导数识别密度体的边缘位置;1985年,Cordell和Grauch[7]将磁力异常换算成磁源重力异常,计算磁源重力异常的总水平导数来识别磁性体的边缘位置;1987年Grauch和Cordell[8]又详细讨论了总水平导数的一些影响因素,如倾斜边界的影响、邻近地质体或异常干扰的影响、不同埋深的影响、区域异常的影响、地形起伏的影响、网格间距的影响、不规则边界的影响等;1994年余钦范和楼海[9]利用二度体倾斜台阶模型研究了总水平导数边缘识别结果的精度;2007年钟清[5]等利用大量简单规则模型的试算结果研究了重力异常总水平导数确定地质体边缘位置的误差分布规律.这些研究工作奠定了总水平导数的理论基础,但还没有系统研究总水平导数极值位置随地质体埋深的空间变化规律.针对这一问题,本文用Green函数将重、磁异常边缘识别方法统一,并研究了位场总水平导数的函数性质以及适应性.通过几种简单规则形体重力异常总水平导数和垂向导数总水平导数的极值位置的解析表达式以及剖面、平面特征,特别是断面特征,重点研究了这些简单规则形体重力异常总水平导数和垂向导数总水平导数的极值位置空间变化规律,为总水平导数边缘识别偏移量大小的研究奠定了一定的理论基础.
2 位场总水平导数边缘识别方法的性质及适应性 2.1 重力异常总水平导数和磁异常总水平导数的关系为研究重、磁异常总水平导数之间的关系,引入Green函数,并将重、磁异常及其总水平导数统一用该Green函数表示.
(1)重、磁异常的Green函数表示式
对于二度体,Green函数的形式为
(1) |
这里
对于三度体,Green函数的形式为
(2) |
这里
不论是二度体或三度体,其重力异常均可表示为
(3) |
但对于二度体磁力异常,其表达式为
(4) |
对于三度体磁力异常,其表达式为
(5) |
以上各式中,G为万有引力常数;σ为剩余密度;μ0为真空磁导率;Vz为Green函数一阶垂向导数;Vxx,Vxy,Vxz,Vyy,Vyz,Vzz为Green函数二阶混合导数;mxx=txMx,mxy=(txMy+tyMx),mxz=(txMz+tzMx),myy=tyMy,myz=(tyMz+tzMy),mzz=tzMz,{Mx,My,Mz}为磁化强度分量,{tx,ty,tz}为磁异常分量方向单位矢量.
(2)二度体重、磁异常总水平导数表达式及适应性
根据上述各式得到二度体重力异常总水平导数的表达式为
(6) |
而二度体重力异常垂向导数总水平导数的表达式为
(7) |
二度体磁力异常总水平导数的表达式为
(8) |
对于二度体化极磁力异常,{tx,tz}={0,1},{Mx,Mz}={0,M},其总水平导数的表达式为
(9) |
对于二度体化赤磁力异常,{tx,tz}={1,0},{Mx,Mz}={M,0},其总水平导数的表达式为
(10) |
由此可以看出,二度体磁力异常总水平导数受磁异常分量方向(tx,tz)和磁化方向(Mx,Mz)的影响,因而不能直接利用磁异常总水平导数的极大值位置来识别地质体的边缘位置.但二度体化极磁力异常和化赤磁力异常总水平导数不受磁异常分量方向和磁化方向的影响,因而可以直接利用其极大值位置来识别地质体的边缘位置.并且化极磁力异常总水平导数与重力异常垂向导数总水平导数的极值位置相同,也与Green函数Vxzz的极值位置相同.
(3)三度体重、磁异常总水平导数表达式及适应性
三度体重力异常总水平导数的表达式为
(11) |
而三度体重力异常垂向导数总水平导数的表达式为
(12) |
三度体磁力异常在x和y方向的一阶水平导数及总水平导数表达式分别为
(13) |
对于三度体化极磁力异常,{tx,ty,tz}={0,0,1},{Mx,My,Mz}={0,0,M},其总水平导数为
(14) |
对于三度体化赤磁力异常,{tx,ty,tz}={tx,ty,0},{Mx,My,Mz}={Mx,My,0},其水平导数和总水平导数表达式分别为
(15) |
可以看出,三度体磁力异常和化赤磁力异常的总水平导数依然受磁异常分量方向(tx,ty,tz)和磁化方向(Mx,My,Mz)的影响,只有化极磁力异常的总水平导数不受其影响.因此,不能直接使用磁异常总水平导数来识别磁性体的边缘位置,必须将其转换成化极磁力异常或磁源重力异常才可以使用.由于化极磁力异常总水平导数和重力异常垂向导数总水平导数的极值位置相同,所以本文不再研究化极磁力异常总水平导数极值位置的空间变化规律,只研究重力异常以及重力异常垂向导数的总水平导数极值位置空间变化规律.
2.2 位场总水平导数的函数性质及“奇点”问题这里通过位函数的性质来研究位场总水平导数的函数性质,并讨论“奇点”问题.对于剖面位场(重力异常或磁力异常)g(x,z),其总水平导数THDR的计算公式为
(16) |
对于平面位场g(x,y,z),其总水平导数THDR的计算公式为
(17) |
这里的gx和gy分别为位场沿x和y方向的一阶水平导数.
剖面位场g(x,z)满足Laplace方程gxx+gzz=0,其总水平导数所满足的方程为
(18) |
即当剖面位场g(x,z)是位函数时,位场总水平导数只有在THDR≠0的条件下才能保持位函数的性质.
平面位场g(x,y,z)满足Laplace方程gxx+gyy+gzz=0,其总水平导数所满足的方程为
(19) |
即当平面位场g(x,y,z)是位函数时,位场总水平导数不能保持位函数的性质.
通过以上证明得出,不论是剖面位场或平面位场,其总水平导数均不能保持位函数的性质.理论上讲不能以位场总水平导数来构造新的位场边缘识别方法,但若以位场总水平导数来构造新的边缘识别方法时,就一定要注意其“奇点”问题.如平面位场总水平导数垂向导数的表达式为
(20) |
当位场总水平导数THDR=0时总水平导数垂向导数出现了无法计算的“解析奇点”.由于该“解析奇点”的存在使得当THDR=0时位场总水平导数垂向导数的计算结果不稳定而发生畸变.这一问题在某些边缘识别方法中也存在,如以倾斜角为基础构造出的倾斜角总水平导数在总水平导数为0或解析信号振幅为0时出现了“奇点”.
3 重力异常总水平导数极值位置空间变化规律重力异常总水平导数极值位置的空间变化规律可以通过求解极值位置的解析表达式来进行研究,若难以得到其解析表达式则必须通过其极值位置的断面变化特征来进行研究.这里将选用单边界模型(铅垂台阶和倾斜台阶)、双边界模型(直立平行四边形和梯形)、多边界模型(直立平行四边形组合)、点质量模型(球体)和三度体模型(直立六面体)来研究重力异常总水平导数极值位置的空间变化规律.
3.1 直立单边界模型--铅垂台阶铅垂台阶模型只有一个直立边界,另外一个边界位于无穷远处,该模型用来模拟单边陡倾角断层或接触面.在直角坐标系中(z坐标方向向下为正),设坐标原点位于铅垂台阶边缘在水平线的投影处.又设铅垂台阶的剩余密度为σ,顶面和底面的z坐标分别为ζ1和ζ2,其引起的重力异常总水平导数(Total Horizontal Derivative)THDR为
(21) |
由(21)不难看出,THDR(-x,0)=THDR(x,0),即铅垂台阶重力异常总水平导数是一个非负偶函数.为了研究重力异常总水平导数的极值位置,对(21)沿x方向求导得
(22) |
并令其为0可得极值位置的解为xm=0,该极值为极大值,位于铅垂台阶边缘处.因此,可以利用重力异常总水平导数极大值位置来准确识别铅垂台阶的边缘位置.
3.2 倾斜单边界模型--倾斜台阶倾斜台阶只有一个倾斜边界,另外一个边界位于无穷远处,该模型用来模拟单边缓倾角断层或接触面.在直角坐标系中(z坐标方向向下为正),设坐标原点位于倾斜台阶上顶面边缘在水平线的投影点,顶面和底面的z坐标分别为ζ1和ζ2,倾斜面和底面的夹角为α,剩余密度为σ,其引起的重力异常总水平导数THDR为
(23) |
对(23)沿x方向求导,并令其为0,经推导后得到极值位置的解为
(24) |
该极值位置就是极大值位置.若倾斜台阶上顶面边缘点坐标为(0,ζ1),下底面边缘点坐标为(ξ2,ζ2),则(24)还可以表示为
(25) |
可以看出该极大值位置既不是倾斜台阶模型的上顶边缘位置,也不是下底边缘位置,是随着倾斜台阶的埋深以及倾斜角而变化的.也就是说,利用重力异常总水平导数的极大值位置能够准确识别铅垂台阶的边缘位置,但不能准确识别倾斜台阶的边缘位置,其识别结果必然有一定偏差(偏移).当ξ2=0时,xm=0;当埋深较大时,xm→ξ2/2.即重力异常总水平导数极大值位置偏离倾斜台阶上顶面边缘位置不超过倾斜面水平宽度的一半,这就是倾斜台阶重力异常总水平导数极大值位置的空间变化规律.
图 1给出了在z=0剖面上倾斜面和底面夹角分别为30°、45°和135°的重力异常总水平导数剖面图和在z=0以上半空间重力异常总水平导数极大值位置断面图(附倾斜面和底面夹角为45°时的重力异常总水平导数等值线图).
(1)从剖面图来看,重力异常总水平导数极大值位置偏向台阶倾斜方向,且在反倾向一侧的梯度变化较大,表现出不对称性;台阶倾斜角越小,重力异常总水平导数极大值位置偏离台阶上顶面边缘位置的距离越大,且倾角为α(45°)的重力异常总水平导数极大值位置与倾角为(π-α)(135°)的极大值位置以坐标原点(倾斜台阶上顶面边缘点)为对称.
(2)从断面图来看,随着倾斜台阶埋深的逐渐增加,重力异常总水平导数极大值位置首先快速向倾向一侧偏移,并很快收敛到一个稳定值--倾斜面水平宽度的一半.
3.3 直立双边界模型--直立平行四边形直立平行四边形模型的两个边界均呈直立状,上、下顶面均呈水平状,该模型用来模拟双边陡倾角断层或接触面.在直角坐标系中(z坐标方向向下为正),直立平行四边形的左、右边界x坐标分别为-ξ和ξ,顶面和底面z坐标分别为ζ1和ζ2,剩余密度为σ.其引起的重力异常总水平导数THDR可以由铅垂台阶重力异常总水平导数表达式(21)导出为
(26) |
式中不难看出,THDR(-x,0)=THDR(x,0),即重力异常总水平导数是一个非负偶函数.对该式沿x方向求导,并令其为0,得
因该方程的判别式为
所以x2有两个不等的实根.又考虑到x2必须大于等于0,故最后求得的重力异常总水平导数极值位置解析表达式为
(27) |
可以看出该极大值位置不是平行四边形的边缘位置,其值随着平行四边形的埋深以及水平尺寸的大小而变化.也就是说,利用重力异常总水平导数的极大值位置也不能准确识别平行四边形的边缘位置,其识别结果必然有一定偏差.
当ξ2远大于ζ1ζ2时,xm→±ξ,即当形体水平尺寸远大于埋深时重力异常总水平导数的极大值位置接近直立平行四边形的上顶边缘位置.
当ξ2远小于ζ1ζ2时,
对于一般情况,|xm|>ξ,表明重力异常总水平导数极大值位置位于形体上顶边缘外侧,随着形体水平尺寸的增加,重力异常总水平导数极大值位置接近形体上顶边缘位置;随着形体埋深的增加,重力异常总水平导数极大值位置接近形体平均埋深的
图 2给出了在z=0剖面上重力异常总水平导数图和在z=0以上半空间重力异常总水平导数极值位置(红色为极大值位置,黄色为极小值位置,也就是零值位置)断面图.由图 2可以看出:直立平行四边形重力异常总水平导数极大值位置位于形体上顶边缘外侧,且随埋深和水平尺寸发生变化.随着埋深的增加,直立平行四边形重力异常总水平导数极大值位置偏离上顶边缘位置的距离增大;随着水平尺寸的增大,重力异常总水平导数极大值位置偏离形体上顶边缘位置的距离变小.
梯形台阶的两个边界均呈倾斜状,上、下顶面均呈水平状,该模型用来模拟双边缓倾角断层或接触面.该模型重力异常总水平导数可以由倾斜台阶模型重力异常总水平导数表达式(23)导出.但难以得到梯形台阶模型重力异常总水平导数极值位置的解析表达式,故用其断面图来研究总水平导数极值位置的空间变化规律.
图 3为z=0剖面上的重力异常总水平导数图和z=0以上半空间重力异常总水平导数极值位置图.根据该图可以看出:重力异常总水平导数剖面图不对称,在形体中心附近达到极小值0,在形体边缘附近达到极大值,并且在反倾向一侧梯度变化较大.当形体埋深较浅时,重力异常总水平导数极大值位置接近上顶边缘位置,且向倾向一侧偏移;随着形体埋深的增加,重力异常总水平导数极大值位置先向倾向一侧偏移,到达一定深度时,开始向形体外侧偏移;随着形体埋深的进一步增大,其重力异常总水平导数极大值位置的变化规律与直立平行四边形的变化规律基本一致,即以形体中心位置为对称向形体外侧偏移.
这里设计两个相同大小的直立平行四边形组合模型来代表多边界模型.两个直立平行四边形水平宽度均为600m,垂向厚度也为600m,相距200m.根据(26)计算得到重力异常总水平导数在z=0上的剖面图和在z=0以上半空间的断面图见图 4所示.由该图可以看出:在两个互相靠近的直立平行四边形之间,埋深较浅时,总水平导数极大值位置位于直立平行四边形的上顶边缘位置;随着埋深的增加,两个形体总水平导数极大值位置首先向形体内侧偏移,之后向形体外侧偏移,直到埋深达到400 m时,两个极大值位置重合、并消失.这就是总水平导数所能分辨的最大深度.另外,随着埋深的增加,各形体外侧边缘的总水平导数极大值位置向形体外侧方向偏移,好像是一个直立平行四边形一样,其变化规律也和一个直立平行四边形总水平导数极大值位置变化规律一致.
球体可以用来研究埋深较大或等轴三度体的重力异常总水平导数极值位置空间变化规律.在直角坐标系中(z坐标方向向下为正),假设球体的剩余密度为σ、半径为R、球体中心坐标为(0,0,ζ).则重力异常总水平导数THDR的表达式为
(28) |
对(28)式沿水平方向(径向)
即球体重力异常总水平导数是关于z轴对称的函数,而且极大值位置到球心的偏移量等于球心埋深的一半.
3.7 直立边界三度体模型--直立六面体该模型用来模拟三度体陡倾角断层以及接触面,在直角坐标系中(z坐标方向向下为正),设直立六面体的坐标范围为ξ1~ξ2,η1~η2,ζ1~ζ2,剩余密度为σ.则重力异常总水平导数THDR的表达式为
(29) |
这里,
直立六面体重力异常总水平导数极值位置难以得到简明的解析表达式,只能通过重力异常总水平导数平面图和中心剖面以上半空间断面图的极大值位置来研究其空间变化规律.图 5为重力异常总水平导数平面等值线图和z=0以上半空间y=0的断面等值线图.
由该图可以看出:当直立六面体埋深较浅时,重力异常总水平导数极大值位置接近形体的上顶边缘位置;随着埋深的增加,重力异常总水平导数极大值位置向形体外侧偏离形体的上顶边缘位置.这些变化规律与二度体直立平行四边形重力异常总水平导数极大值位置的空间变化规律相似,只是偏移量的大小不同.根据球体重力异常总水平导数极值位置的空间变化规律知,当埋深较大时,其极大值位置接近形体埋深的一半.
4 重力异常垂向导数总水平导数极值位置空间变化规律重力异常垂向导数总水平导数的极值位置与化极磁力异常总水平导数的极值位置相同,其研究方法与重力异常总水平导数极值位置的研究方法也相同,选择的模型也相同.
4.1 直立单边界模型--铅垂台阶铅垂台阶模型重力异常垂向导数总水平导数(Total Horizontal Derivative of Vertical Derivative)VDR_THDR为
(30) |
由(30)不难看出,VDR_THDR(-x,0)=VDR_THDR(x,0),即重力异常垂向导数总水平导数是一个非负偶函数.对(30)沿x方向求导并令其为0,得
解得
(31) |
可以看出,铅垂台阶模型重力异常垂向导数总水平导数有三个极值点,分别位于铅垂台阶边缘处以及两侧.在铅垂台阶边缘处的极大值幅值较大,称为“主极大值”;而在两侧的极大值幅值较小,称为“次极大值”,与重力异常总水平导数峰值比较更加尖锐,其横向分辨能力更强.图 6是铅垂台阶模型重力异常垂向导数总水平导数剖面图和z=0以上半空间极大值位置断面图,由该图看出,重力异常垂向导数总水平导数“主极大值”位置不随铅垂台阶埋深的变化而变化,一直位于铅垂台阶边缘处,但“次极大值”位置随铅垂台阶埋深变化而变化.
该模型用来研究埋深较大或线质量二度体重力异常垂向导数总水平导数极值位置的空间变化规律.在直角坐标系中(z坐标方向向下为正),设无限延伸水平圆柱体的半径为R、圆柱体中心坐标为(0,ζ),剩余密度为σ.则其重力异常垂向导数总水平导数VDR_THDR为
(32) |
由(32)可看出,
即重力异常垂向导数总水平导数是一个非负偶函数,对其沿x方向求导并令其为0,得
解得
(33) |
即无限延伸水平圆柱体重力异常垂向导数总水平导数是一个非负偶函数,有四个对称的极值点,分别位于无限延伸水平圆柱体中心两侧.其中x1,2=±(
倾斜台阶模型重力异常垂向导数总水平导数VDR_THDR为
(34) |
由于难以得到重力异常垂向导数总水平导数极值位置的解析表达式,故只能通过剖面和断面特征研究其极值位置的空间变化规律.图 7给出了在z=0剖面上倾斜面和底面夹角分别为30°、45°和135°的重力异常垂向导数总水平导数图和在z=0以上半空间断面图(附倾斜面和底面夹角为45°时的重力异常垂向导数总水平导数等值线图).从剖面图来看,重力异常垂向导数总水平导数极大值位置偏向台阶倾斜方向,且在反倾向一侧的梯度变化较大,表现出不对称性;台阶倾斜角越小,重力异常垂向导数总水平导数极大值位置偏离台阶上顶面边缘位置的距离越大,且倾角为α(45°)的和倾角为(π-α)(135°)的重力异常垂向导数总水平导数极大值位置以坐标原点(倾斜台阶上顶面边缘点)为对称.从断面图来看,随着倾斜台阶埋深的逐渐增加,重力异常垂向导数总水平导数极大值位置首先快速向倾向一侧偏移,并很快收敛到一个稳定值,比倾斜台阶重力异常总水平导数极值位置的收敛速度快(即与倾斜台阶上顶面边缘位置的偏移量小).
直立平行四边形模型重力异常垂向导数总水平导数VDR_THDR为
(35) |
式中可见,VDR_THDR(-x,0)=VDR_THDR(x,0),即重力异常垂向导数总水平导数是非负偶函数.由于难以得到重力异常垂向导数总水平导数极值位置的解析表达式,故只能通过剖面和断面特征研究其极值位置的空间变化规律.图 8为z=0剖面上和z=0以上半空间断面上重力异常垂向导数总水平导数图.从剖面图来看,重力异常垂向导数总水平导数(图 8)的峰值比重力异常总水平导数(图 2)的峰值尖锐,但存在“次极大值”.从断面图来看,重力异常垂向导数总水平导数极大值位置与重力异常总水平导数极大值位置的变化规律相似,但比重力异常总水平导数的极大值位置更靠近形体上顶边缘位置.当埋深不大时,“主极大值”的位置基本保持不变,靠近地质体上顶边缘位置;当埋深逐渐增大时,“主极大值”位置逐渐自形体上顶边缘向形体外侧偏移,最后收敛到埋深的0.414倍,与无限延伸水平圆柱体重力异常垂向导数总水平导数“主极大值”位置的空间变化规律一致.
梯形台阶模型重力异常垂向导数总水平导数可以由倾斜台阶重力异常垂向导数总水平导数表达式导出.但由于难以得到梯形台阶重力异常垂向导数总水平导数极值位置的解析表达式,故用断面图来研究其极值位置的空间变化规律.
图 9为z=0剖面上的重力异常垂向导数总水平导数图和z=0以上半空间重力异常垂向导数总水平导数极值位置图.从剖面特征图来看,重力异常垂向导数总水平导数的峰值比总水平导数的峰值尖锐,但存在“次极大值”.从断面图来看,重力异常垂向导数总水平导数与重力异常总水平导数极大值位置的空间变化规律相似,但比重力异常总水平导数极大值位置更靠近形体上顶边缘位置.当形体埋深较浅时,重力异常垂向导数总水平导数“主极大值”位置向倾向一侧偏移;到达一定深度时开始向整个形体上顶边缘外侧偏移,最后收敛到形体埋深的0.414倍,与无限延伸水平圆柱体垂向导数总水平导数的极大值位置变化规律一致.
这里设计的边界模型与图 4相同.根据(35)计算得到重力异常垂向导数总水平导数在z=0上的剖面图和在z=0以上半空间的断面图见图 10所示.由该图可以看出:在两个互相靠近的直立平行四边形之间,埋深较浅时,垂向导数总水平导数极大值位置位于直立平行四边形的上顶边缘位置;随着埋深的增加,垂向导数总水平导数极大值位置首先向形体内侧偏移,之后向形体外侧偏移,直到埋深达到660m时,两个极大值位置重合、并消失,这就是垂向导数总水平导数所能分辨的最大深度.另外,随着埋深的增加,各形体外侧边缘的垂向导数总水平导数极大值位置向形体外侧方向偏移,与单一直立平行四边形垂向导数总水平导数极大值位置变化规律一致.对比图 10和图 4可以看出,重力异常垂向导数总水平导数比总水平导数极大值位置更靠近形体的边缘位置,且其分辨能力更强.
球体重力异常垂向导数总水平导数VDR_THDR的表达式为
(36) |
这里
由此解得
(37) |
此式中,r1=0和r2=2ζ为极小值位置,r3=0.389ζ为“主极大值”位置,r4=2.5687ζ为“次极大值”位置.与重力异常总水平导数极大值位置(0.5ζ)比较,球体重力异常垂向导数总水平导数“主极大值”位置更靠近球体的中心位置.
4.8 直立边界三度体模型--直立六面体直立六面体重力异常垂向导数总水平导数VDR_THDR的表达式为
(38) |
其中,
由于直立六面体重力异常垂向导数总水平导数极值位置的解析式难以得到,故通过平面图和断面图来研究其空间变化规律.图 11为重力异常垂向导数总水平导数平面等值线图和z=0以上半空间y=0的断面等值线图.从其平面变化特征来看,重力异常垂向导数总水平导数(图 11)极大值位置比重力异常总水平导数(图 5)极大值位置更靠近形体的上顶边缘位置,并且其峰值更加尖锐.从其断面图来看,随着埋深的增大,重力异常垂向导数总水平导数极大值位置自形体外侧偏离形体的上顶边缘位置,最后逐渐收敛到中心埋深的0.389倍,与球体重力异常垂向导数总水平导数的极值位置空间变化规律相同.
通过对位场总水平导数的函数性质以及对重力异常总水平导数和重力异常垂向导数总水平导数极大值位置的空间变化规律研究,得出以下结论与建议:
(1)位场总水平导数不是位函数,不能用位场总水平导数来构造新的边缘识别方法,否则会在位场总水平导数为0的点上出现“奇点”.
(2)磁力异常总水平导数受磁场分量方向和磁化方向的影响,故在利用磁力异常进行边缘识别时必须将其转换成磁源重力异常(假重力异常)或化极磁力异常才可以使用.
(3)在单一直立边界上,利用重力异常总水平导数极大值位置能够准确识别直立边界位置,而在其他任何模型上利用重力异常总水平导数极大值位置均不能准确识别其边缘位置,其偏移量随埋深、水平尺寸、倾斜角等因素发生变化.在单一倾斜边界上,重力异常总水平导数极大值位置偏向倾向一侧,并且随埋深的增大而增加,但偏移量不超过倾斜面水平宽度的一半;在双边界上,重力异常总水平导数极大值位置的空间变化规律与单一边界的空间变化规律差别很大,随埋深的增大,重力异常总水平导数极大值位置自形体上顶边缘位置向倾向一侧偏移,之后向形体外侧偏移,最终收敛于平均埋深的
(4)在单一直立边界上,利用重力异常垂向导数总水平导数极大值位置能够准确识别直立边界的位置,而在其他任何模型上的重力异常垂向导数总水平导数极大值位置与其边缘位置均有偏移,其偏移量大小随形体埋深、水平尺寸、倾斜角等因素发生变化,且小于重力异常总水平导数极大值位置的偏移量.重力异常垂向导数总水平导数比重力异常总水平导数的峰值更加尖锐、横向分辨能力更强,且出现了“次极大值”的影响.对于二度体,重力异常垂向导数总水平导数“主极大值”位置收敛于平均埋深的0.414倍;对于三度体,重力异常垂向导数总水平导数“主极大值”位置收敛于平均埋深的0.389倍.
(5)通过掌握重力异常总水平导数和重力异常垂向导数总水平导数极值位置的空间变化规律可以进一步提高位场总水平导数边缘识别结果的精度,并有可能通过不同高度极值位置的变化规律来研究边缘深度的参数反演.
致谢感谢潘作枢教授、朱光明教授以及审稿专家对本文提出的修改意见.
[1] | 王万银, 邱之云, 杨永, 史维皎. 位场边缘识别方法研究进展. 地球物理学进展 , 2009, 25(1): 196–210. Wang W Y, Qiu Z Y, Yang Y, Shi W J. Some advances in the edge recognition of the potential field. Progress in Geophysics (in Chinese) , 2009, 25(1): 196-210. |
[2] | Hood P, McClure D J. Gradient measurements in ground magnetic prospecting. Geophysics , 1965, 30(3): 403-410. DOI:10.1190/1.1439592 |
[3] | Bhattacharyya B K. Two-dimensional harmonic analysis as a tool for magnetic interpretation. Geophysics , 1965, 30(5): 829-857. DOI:10.1190/1.1439658 |
[4] | Hood P J, Teskey D J. Aeromagnetic gradiometer program of the Geological Survey of Canada. Geophysics , 1989, 54(8): 1012-1022. DOI:10.1190/1.1442726 |
[5] | 钟清, 孟小红, 刘士毅. 重力资料定位地质体边界问题的探讨. 物探化探计算技术 , 2007, 29(Suppl): 35–38. Zhong Q, Meng X H, Liu S Y. Study on the Geological body's edge detection of using gravity data. Computing Techniques for Geophysical and Geochemical Exploration (in Chinese) , 2007, 29(Suppl): 35-38. |
[6] | Cordell L. Gravimetric expression of graven faulting in Santa Fe Country and the Espanola Basin. New Mexico: New Mexico Geol. Soc. Guidebook, 30th Field Conf, 1979 : 59 -64. |
[7] | Cordell L, Grauch V J S. Mapping basement magnetization zones from aeromagnetic data in the San Juan Basin, New Mexico, in Hinze, W. J, Ed, The utility of regional gravity and magnetic anomaly maps:Soc. Explor. Geophys., 1985.181~197 |
[8] | Grauch V J S, Cordell L. Limitations of determining density or magnetic boundaries form the horizontal gradient of gravity or pseudo-gravity data. Geophysics , 1987, 52(1): 118-121. DOI:10.1190/1.1442236 |
[9] | 余钦范, 楼海. 水平梯度法提取重磁源边界位置. 物探化探计算技术 , 1994, 16(4): 363–367. Yu Q F, Lou H. Locating the boundaries of magnetic or gravity sources using horizontal gradient anomalies. Computing Techniques for Geophysical and Geochemical Exploration (in Chinese) , 1994, 16(4): 363-367. |
[10] | Nabighian M N. The analytic signal of two-dimensional magnetic bodies with polygonal cross-section:its properties and use for automated anomaly interpretation. Geophysics , 1972, 37(3): 507-517. DOI:10.1190/1.1440276 |
[11] | Nabighian M N. Toward a three dimensional automatic interpretation of potential field data via generalized Hilbert transforms:Fundamental relations. Geophysics , 1984, 49(6): 780-786. DOI:10.1190/1.1441706 |
[12] | Nelson J.B. Comparison of gradient analysis techniques for linear two-dimensional magnetic sources. Geophysics , 1988, 53(8): 1088-1095. DOI:10.1190/1.1442545 |
[13] | Roest W R, Verhoef J, Pilkington M. Magnetic interpretation using the 3-D analytic signal. Geophysics , 1992, 57(1): 116-125. DOI:10.1190/1.1443174 |
[14] | Qin S. An analytic signal approach to the interpretation of total field magnetic anomalies. Geophysical Prospecting , 1994, 42(6): 665-675. DOI:10.1111/gpr.1994.42.issue-6 |
[15] | Miller H G, Singh V. Potential field tilt-a new concept for location of potential field sources. Journal of Applied Geophysics , 1994, 32: 213-217. DOI:10.1016/0926-9851(94)90022-1 |
[16] | Wijns C, Perez C, Kowalczyk P. Theta map:Edge detection in magnetic data. Geophysics , 2005, 70(4): L39-L43. DOI:10.1190/1.1988184 |