CSAMT中的IP效应提取
|
摘要
可控源音频大地电磁法(CSAMT)是地球物理方法的重要组成部分。它是采用可以控制的人工场源的一种电磁法,其无论是在理论和应用方面都获得了成功。
激发极化效应(IP)是电荷聚集在金属矿物和电解液的分界面上引起的,或者是黏土或孔隙束缚离子运动而引起的。IP是一种很有意义的参数,因为它可以用来探测有经济价值的硫化矿带的赋存特点和地质构造。但是,试图从CSAMT中提取有效的IP信息,国内外学者也进行了一定的研究,至今还没有明确的定论。
本文首先回顾了从电磁法勘探中提取IP信息的发展,总结了国内外的研究现状。详细阐述了可控源音频大地电磁法的基本理论,并进一步探讨了电偶极子场的一维正演问题。接着研究了CSAMT中的IP影响,而且利用柯尔-柯尔模型模拟了可极化大地上CSAMT的IP效应。分析了电磁场各分量的特点,并且提出利用电场y分量来提取IP信息的可行性。初步尝试提取有效的IP信息。本文最后对三层介质上中间为极化层的水平电偶源视电阻率进行了理论计算。着重讨论了激电效应对于不同地电断面类型的影响,试图探求其变化规律,为CSAMT实际勘探资料中提取激电信息提供依据。通过以上研究表明,在可极化均匀大地上,IP效应对磁场y分量影响很小。而电场y分量能够很好的反应出IP效应。在层状介质中,从整体效果来讲,电场y分量定义的视电阻率对IP效应的反应好于其他定义的视电阻率。对于不同的地电断面,H型地电断面有利于发现中间极化层的存在。视电阻率对于低阻高极化层反应比较好。对于O型地电断面,中间极化层对视电阻率影响很小,很难反应出IP信息。
The Controlled Source Audio-frequency Magnetotellurics(CSAMT) is an important part of the geophysics.It is one kind of controllable source electromagnetic method which both in theory and applications have been successful.
Effect of induced polarization(IP) is the charge accumulation in metal minerals and dectrolytes caused by the boundary surface,or the clay or the hole fettering ion movement caused.IP is a significant parameter,because it can be used to detect the economic valuable sulphide ore zone characteristics and geological structure.But trying to extract from the CSAMT effective IP information,on which the domestic and foreign scholars have also conducted certain research,until now did not have the explicit conclusion.
This article first reviewed the development of extracting the IP information from the electromagnetic prospecting and summarized the domestic and foreign research present situation.Elaborated in detail the CSAMT method's basic theory, further discussed the forward problem of electromagnetic of the horizontal electrical dipole.Then has studied IP influence in the CSAMT,moreover model has simulated the CSAMT method's IP effect using Cole-Cole model on the polarizable earth.Has analyzed the electromagnetic field various components characteristic,and proposed that regarding extraction the IP information using the electric field y component the new development.Preliminary attempted to extract effective IP information.Finally, we theoretical calculated the apparent resistivity of dectromagnetic of the horizontal electrical dipole under the three-layer which the middle layer is polarized.Discussed the influence of IP on differ section of resistivity,we attempted to seek its change rule.It provides grounds for extracting IP information from CSAMT.Through the above-mentioned studies have shown that on the polarizable earth,IP effect influences H_y very small.E_y can well reflect IP effect.In layered earth,from the overall effect is concerned,the apparent resistivity of the electric field y-component definition response to IP effect better than others.As to differ section of resistivity, H section of resistivity help to find the existence of middle layer polarization.Apparent resistivity for the low resistance layer high-polarization response is better.As to Q section of resistivity,the middle layer polarization influence apparent resistivity very much little,it is difficult to reflect the IP information.
激发极化效应(IP)是电荷聚集在金属矿物和电解液的分界面上引起的,或者是黏土或孔隙束缚离子运动而引起的。IP是一种很有意义的参数,因为它可以用来探测有经济价值的硫化矿带的赋存特点和地质构造。但是,试图从CSAMT中提取有效的IP信息,国内外学者也进行了一定的研究,至今还没有明确的定论。
本文首先回顾了从电磁法勘探中提取IP信息的发展,总结了国内外的研究现状。详细阐述了可控源音频大地电磁法的基本理论,并进一步探讨了电偶极子场的一维正演问题。接着研究了CSAMT中的IP影响,而且利用柯尔-柯尔模型模拟了可极化大地上CSAMT的IP效应。分析了电磁场各分量的特点,并且提出利用电场y分量来提取IP信息的可行性。初步尝试提取有效的IP信息。本文最后对三层介质上中间为极化层的水平电偶源视电阻率进行了理论计算。着重讨论了激电效应对于不同地电断面类型的影响,试图探求其变化规律,为CSAMT实际勘探资料中提取激电信息提供依据。通过以上研究表明,在可极化均匀大地上,IP效应对磁场y分量影响很小。而电场y分量能够很好的反应出IP效应。在层状介质中,从整体效果来讲,电场y分量定义的视电阻率对IP效应的反应好于其他定义的视电阻率。对于不同的地电断面,H型地电断面有利于发现中间极化层的存在。视电阻率对于低阻高极化层反应比较好。对于O型地电断面,中间极化层对视电阻率影响很小,很难反应出IP信息。
The Controlled Source Audio-frequency Magnetotellurics(CSAMT) is an important part of the geophysics.It is one kind of controllable source electromagnetic method which both in theory and applications have been successful.
Effect of induced polarization(IP) is the charge accumulation in metal minerals and dectrolytes caused by the boundary surface,or the clay or the hole fettering ion movement caused.IP is a significant parameter,because it can be used to detect the economic valuable sulphide ore zone characteristics and geological structure.But trying to extract from the CSAMT effective IP information,on which the domestic and foreign scholars have also conducted certain research,until now did not have the explicit conclusion.
This article first reviewed the development of extracting the IP information from the electromagnetic prospecting and summarized the domestic and foreign research present situation.Elaborated in detail the CSAMT method's basic theory, further discussed the forward problem of electromagnetic of the horizontal electrical dipole.Then has studied IP influence in the CSAMT,moreover model has simulated the CSAMT method's IP effect using Cole-Cole model on the polarizable earth.Has analyzed the electromagnetic field various components characteristic,and proposed that regarding extraction the IP information using the electric field y component the new development.Preliminary attempted to extract effective IP information.Finally, we theoretical calculated the apparent resistivity of dectromagnetic of the horizontal electrical dipole under the three-layer which the middle layer is polarized.Discussed the influence of IP on differ section of resistivity,we attempted to seek its change rule.It provides grounds for extracting IP information from CSAMT.Through the above-mentioned studies have shown that on the polarizable earth,IP effect influences H_y very small.E_y can well reflect IP effect.In layered earth,from the overall effect is concerned,the apparent resistivity of the electric field y-component definition response to IP effect better than others.As to differ section of resistivity, H section of resistivity help to find the existence of middle layer polarization.Apparent resistivity for the low resistance layer high-polarization response is better.As to Q section of resistivity,the middle layer polarization influence apparent resistivity very much little,it is difficult to reflect the IP information.
引文
[1]汤井田,何继善.可控音频大地电磁法及其应用[M].湖南:中南大学出版社,2005
[2]何继善.地电磁场及其研究前沿[J].中南大学学报,1995,26(1):124-129
[3]何继善编译.可控音频大地电磁法[M].湖南:中南工业大学出版社,1990
[4]傅良魁.应用地球物理教程[M].北京:地质出版社,1991
[5]傅良魁.电法勘探教程[M].北京:地质出版社,1983
[6]张塞珍,王庆乙,罗延钟.中国电法勘探发展概况.地球物理学报,1994,37(增刊):257-260
[7]Ward S H.地球物理电磁场论[M].北京:地质出版社,1978
[8]Wait J R.大地电磁学[M].北京:地质出版社,1987
[9]A.A.考夫曼,G.V.凯轨.大地电磁测深法[M].北京:地质出版社,1987
[10]汤井田,何继善.频率测深中磁场水平分量计算的新方法[M].物化探计算技术,1993,15(1):123-129
[11]刘崧,官善友,高鹏飞.求极化椭球体真Cole-Cole参数的联合谱激电正演[J].地球物理学报,1994,37(增刊):542-551
[12]Comelis F,Du Toit.The numerical computation of Bessel function of the first and second kind for integer orders and complex argument[J].IEEE Trans on Antennas and Propagation,1990,38(9):1341-1349
[13]Anderson W L.A hybrid fast Hankel transform algorithm for electromagnefic modeling[J].Geophysics,1989,54(2):263-266.
[14]Christensen N B.Optimized fast Hankel transform filters[J].Geophysical Prospecting,1990,38(1):545-558.
[15]Snyder D D.A method for modeling the resistivity and IP response of two-dimensional bodies[J].Geophysics,1976,41:997-1015
[16]周熙襄,钟本善,等.电法勘探数值模拟技术[M].成都:四川科学技术出版社,1986
[17]Comelis F,Du Toit.The numerical computation of Bessel function of the first and second kind for integer orders and complex argument[J].IEEE Trans on Antennas and Propagation,1990,38(9):1341-1349
[18]彭仲秋.电磁波探测地下异常体的数值模拟和物理模型[J].地球物理学报,1989,32(2):354-358
[19]Lu X.Rapid relaxation inversion of CSAMT data[j].GeophysJInt,1999,138(1):381-392
[20]Zonge K L,Wynn J C.Recent advances and application in complex resistivity measurements[J].Geophysics,1975,40(5):851-864
[21]陈乐寿,等.大地电磁测深资料处理与解释[M].北京:石油工业出版社,1989
[22]王建谋,等泽.在地电磁测深译文集(第一集)[M].北京:地质出版社,1987
[23]Liu Song.A new IP decoupling scheme[J].Geophysics,1984(5),15:89-112
[24]刘崧.计算视复电阻率的新的近似公式[J].地球物理学报,1988,31(4):687-694
[25]万尼安,等.非均匀介质大地电磁测深资料解释[M].北京:地质出版社,1988
[26]日丹诺夫.电法勘探[M].武汉:中国地质大学出版社,1990
[27]马特维耶夫.电磁测深解释[M].北京:煤炭工业出版社,1978
[28]石昆法.可控源音频大地电磁法理论与应用[M].北京:科学出版社,1999
[29]吴璐萍,石昆法.可控源音频大地电磁法在地下水勘查中的应用研究[J].地球物理学报,1996,39(5):62-59
[30]王能伦.关于直流电法模型实验中的微发极化效应问题[J].地质与勘探,1965,3(4):24-31
[31]傅良魁,等.频谱激电法若干理论研[J].物探与化探,1985,9(8):45-49
[32]叶春飞,章文勋.有耗电磁场问题中复宗量整数阶圆柱函数的计算[J].微波学报,1993,32(1):54-58
[33]奚定平.贝塞尔函数IM].北京:高等教育出版社,1998
[34]岳安平,底青云,石昆法.从CSAMT信号中提取IP信息探讨[J].地球物理学进展,2007,12(6):1925-1300
[35]Guptasarma D.Effect of surface polarization on resistivity modeling[J].Geophysics,1983,48(2):98-106.
[36]蔡军涛,阮百尧,罗润林.一种快速准确求取复电阻率真频参数的反演方法[J].工程地球物理学报,2005,10(5):338-342
[37]曹中林,昌彦君,何展翔.基于演化算法的复电阻率频谱参数反演[J].工程地球物理学报,2005,2(1):33-37
[38]罗延钟,张胜业,熊彬.天然场源激电法的可行性[J].地球物理学报,2003,1(1):125-130
[39]殷长春.提高频率测深正演计算速度的方法[J].物化探计算技术,1990,12(2):45-49
[40]Brown R J.EM coupling in multi-frequency IP and a generalization of the Cole-Cole impedance model[J].Geophysical prospecting,1985,33(2):282-302
[41]罗延钟,方胜.极化椭球体上频谱激电法的异常形态[J].桂林冶金地质学院学报,1985,5(4):49-63
[42]罗延钟,张桂青.频率域激电法原理[M].北京:地质出版社,1988
[43]傅良魁.复电阻率法异常的频谱及空间分布规律[J].地质与勘探,1981,2(3):46-53
[44]傅良魁,李金铭,陈兆洪.埋藏极化体激电时间谱视参数的试验研究结果[J].地质与勘探,1985,12(5):38-43
[45]殷长春.可控源音频磁大地电流法一维正演及精度评价[J].长春地质学院学报,1994,24(4):438-453.
[46]Lee T J.Transient electromagnetic response of a polarizable ground[J].Geophysics,1981,46(8):1037-1041
[47]Raiche A P.Negative transient voltage and magnetic field responses for a half-space with Cole-Cole impedance[J].Geophysics,1983,48(6):790-791
[48]Marshall D J,Madden T R.Induced polarization.A study of its causes[J].Geophysics,1959,24(9):790-816
[49]Major J,Silic J.Restrictions on the use of Cole-Cole dispersion models in complex resistivity interpretation[J].CJeophysics,1981,46(6):916-931
[50]Luo Y,Zhang G.Theory and application of spectral induced polarization[M].Soc.Expl.Geophys,1998
[51]陈清礼,胡文宝,李金铭.由MT资料反演真谱参数的基本原理[J].石油天然气学报,2006,28(6):61-64.
[52]朴化容,薛爱民.层状介质中三度导电体的电磁响应[J].科学研究论文集,1986,1(1):25-29
[53]Lewis J G,Lee T J.The detection of induced polarization with a transient dectromagnetic system,Inst Elect Electron Eng[J].Trans Geosic Remote Sensing,1984,22(8):69-80
[54]Guptasarma D.Tree and apparent spectra of buried polarizable targets[J].Geophysics,1984,49(6):1534-1540
[55]陈久平,等.层状介质中三维大地电磁模拟[J].地球物理学报,1990,33(4):56-62
[56]朴化荣,殷长春.频率测深水平磁场的正演计算及应用[J].物化探计算技术,1988,11(3):204-213.
[57]Lee T J,Thomas L.The transient electromagnetic response of a polarizable sphere in a conducting half space[J].Geophysical Prospecting,1992,40(12):541-563.
[58]Coggon J H.New three-point formulas for inductive coupling removal in IP[J].Geophysics,1984,49(3):307-309
[59]阎述,陈明生.频率测深中磁性层的视电阻率响应.[J].地质与勘探,1997,33(3):38-42.
[60]杨进,谭捍东,傅良魁.被动源激发极化法的野外试验结果[J].现代地质,1998,9(3):437-441
[61]王若,王妙月.一维全资料CSAMT反演[J].石油地球物理勘探,2007,42(1):108-114.
[62]徐凯军.2.5维复电阻率电磁场正反演研究:[博士学位论文].长春:吉林大学,2007
[63]Yu Changming.Application of CSAMT method on gold mineral deposits exploration[J].Chinese J Geophys(Acta Geophysica Sinica),1998,41(1):133-138.
[2]何继善.地电磁场及其研究前沿[J].中南大学学报,1995,26(1):124-129
[3]何继善编译.可控音频大地电磁法[M].湖南:中南工业大学出版社,1990
[4]傅良魁.应用地球物理教程[M].北京:地质出版社,1991
[5]傅良魁.电法勘探教程[M].北京:地质出版社,1983
[6]张塞珍,王庆乙,罗延钟.中国电法勘探发展概况.地球物理学报,1994,37(增刊):257-260
[7]Ward S H.地球物理电磁场论[M].北京:地质出版社,1978
[8]Wait J R.大地电磁学[M].北京:地质出版社,1987
[9]A.A.考夫曼,G.V.凯轨.大地电磁测深法[M].北京:地质出版社,1987
[10]汤井田,何继善.频率测深中磁场水平分量计算的新方法[M].物化探计算技术,1993,15(1):123-129
[11]刘崧,官善友,高鹏飞.求极化椭球体真Cole-Cole参数的联合谱激电正演[J].地球物理学报,1994,37(增刊):542-551
[12]Comelis F,Du Toit.The numerical computation of Bessel function of the first and second kind for integer orders and complex argument[J].IEEE Trans on Antennas and Propagation,1990,38(9):1341-1349
[13]Anderson W L.A hybrid fast Hankel transform algorithm for electromagnefic modeling[J].Geophysics,1989,54(2):263-266.
[14]Christensen N B.Optimized fast Hankel transform filters[J].Geophysical Prospecting,1990,38(1):545-558.
[15]Snyder D D.A method for modeling the resistivity and IP response of two-dimensional bodies[J].Geophysics,1976,41:997-1015
[16]周熙襄,钟本善,等.电法勘探数值模拟技术[M].成都:四川科学技术出版社,1986
[17]Comelis F,Du Toit.The numerical computation of Bessel function of the first and second kind for integer orders and complex argument[J].IEEE Trans on Antennas and Propagation,1990,38(9):1341-1349
[18]彭仲秋.电磁波探测地下异常体的数值模拟和物理模型[J].地球物理学报,1989,32(2):354-358
[19]Lu X.Rapid relaxation inversion of CSAMT data[j].GeophysJInt,1999,138(1):381-392
[20]Zonge K L,Wynn J C.Recent advances and application in complex resistivity measurements[J].Geophysics,1975,40(5):851-864
[21]陈乐寿,等.大地电磁测深资料处理与解释[M].北京:石油工业出版社,1989
[22]王建谋,等泽.在地电磁测深译文集(第一集)[M].北京:地质出版社,1987
[23]Liu Song.A new IP decoupling scheme[J].Geophysics,1984(5),15:89-112
[24]刘崧.计算视复电阻率的新的近似公式[J].地球物理学报,1988,31(4):687-694
[25]万尼安,等.非均匀介质大地电磁测深资料解释[M].北京:地质出版社,1988
[26]日丹诺夫.电法勘探[M].武汉:中国地质大学出版社,1990
[27]马特维耶夫.电磁测深解释[M].北京:煤炭工业出版社,1978
[28]石昆法.可控源音频大地电磁法理论与应用[M].北京:科学出版社,1999
[29]吴璐萍,石昆法.可控源音频大地电磁法在地下水勘查中的应用研究[J].地球物理学报,1996,39(5):62-59
[30]王能伦.关于直流电法模型实验中的微发极化效应问题[J].地质与勘探,1965,3(4):24-31
[31]傅良魁,等.频谱激电法若干理论研[J].物探与化探,1985,9(8):45-49
[32]叶春飞,章文勋.有耗电磁场问题中复宗量整数阶圆柱函数的计算[J].微波学报,1993,32(1):54-58
[33]奚定平.贝塞尔函数IM].北京:高等教育出版社,1998
[34]岳安平,底青云,石昆法.从CSAMT信号中提取IP信息探讨[J].地球物理学进展,2007,12(6):1925-1300
[35]Guptasarma D.Effect of surface polarization on resistivity modeling[J].Geophysics,1983,48(2):98-106.
[36]蔡军涛,阮百尧,罗润林.一种快速准确求取复电阻率真频参数的反演方法[J].工程地球物理学报,2005,10(5):338-342
[37]曹中林,昌彦君,何展翔.基于演化算法的复电阻率频谱参数反演[J].工程地球物理学报,2005,2(1):33-37
[38]罗延钟,张胜业,熊彬.天然场源激电法的可行性[J].地球物理学报,2003,1(1):125-130
[39]殷长春.提高频率测深正演计算速度的方法[J].物化探计算技术,1990,12(2):45-49
[40]Brown R J.EM coupling in multi-frequency IP and a generalization of the Cole-Cole impedance model[J].Geophysical prospecting,1985,33(2):282-302
[41]罗延钟,方胜.极化椭球体上频谱激电法的异常形态[J].桂林冶金地质学院学报,1985,5(4):49-63
[42]罗延钟,张桂青.频率域激电法原理[M].北京:地质出版社,1988
[43]傅良魁.复电阻率法异常的频谱及空间分布规律[J].地质与勘探,1981,2(3):46-53
[44]傅良魁,李金铭,陈兆洪.埋藏极化体激电时间谱视参数的试验研究结果[J].地质与勘探,1985,12(5):38-43
[45]殷长春.可控源音频磁大地电流法一维正演及精度评价[J].长春地质学院学报,1994,24(4):438-453.
[46]Lee T J.Transient electromagnetic response of a polarizable ground[J].Geophysics,1981,46(8):1037-1041
[47]Raiche A P.Negative transient voltage and magnetic field responses for a half-space with Cole-Cole impedance[J].Geophysics,1983,48(6):790-791
[48]Marshall D J,Madden T R.Induced polarization.A study of its causes[J].Geophysics,1959,24(9):790-816
[49]Major J,Silic J.Restrictions on the use of Cole-Cole dispersion models in complex resistivity interpretation[J].CJeophysics,1981,46(6):916-931
[50]Luo Y,Zhang G.Theory and application of spectral induced polarization[M].Soc.Expl.Geophys,1998
[51]陈清礼,胡文宝,李金铭.由MT资料反演真谱参数的基本原理[J].石油天然气学报,2006,28(6):61-64.
[52]朴化容,薛爱民.层状介质中三度导电体的电磁响应[J].科学研究论文集,1986,1(1):25-29
[53]Lewis J G,Lee T J.The detection of induced polarization with a transient dectromagnetic system,Inst Elect Electron Eng[J].Trans Geosic Remote Sensing,1984,22(8):69-80
[54]Guptasarma D.Tree and apparent spectra of buried polarizable targets[J].Geophysics,1984,49(6):1534-1540
[55]陈久平,等.层状介质中三维大地电磁模拟[J].地球物理学报,1990,33(4):56-62
[56]朴化荣,殷长春.频率测深水平磁场的正演计算及应用[J].物化探计算技术,1988,11(3):204-213.
[57]Lee T J,Thomas L.The transient electromagnetic response of a polarizable sphere in a conducting half space[J].Geophysical Prospecting,1992,40(12):541-563.
[58]Coggon J H.New three-point formulas for inductive coupling removal in IP[J].Geophysics,1984,49(3):307-309
[59]阎述,陈明生.频率测深中磁性层的视电阻率响应.[J].地质与勘探,1997,33(3):38-42.
[60]杨进,谭捍东,傅良魁.被动源激发极化法的野外试验结果[J].现代地质,1998,9(3):437-441
[61]王若,王妙月.一维全资料CSAMT反演[J].石油地球物理勘探,2007,42(1):108-114.
[62]徐凯军.2.5维复电阻率电磁场正反演研究:[博士学位论文].长春:吉林大学,2007
[63]Yu Changming.Application of CSAMT method on gold mineral deposits exploration[J].Chinese J Geophys(Acta Geophysica Sinica),1998,41(1):133-138.