电性源瞬变电磁B场全程视电阻率计算
|
摘要
在传统电性源瞬变电磁资料解释中,常常采用Bz/T求取视电阻率值,由此带来Bz/T计算全程视电阻率多解性.而Bz随电阻率变化是单值函数,由Bz计算全程视电阻率具有减少多解性的优势.为了研究接地源瞬变电磁B场全程视电阻率,文中首先研究了接地源瞬变电磁法Bz和Bz/T响应值随电阻率的变化特性,然后设计不同的地电模型,采用二分搜索算法对B场全程视电阻率进行计算,结果表明全程视电阻率对电性目标体分辨率高,分层能力强,能较好地反映地电断面的信息;另外,全程视电阻率仅与地层的电性结构有关,与收发距无关,便于野外进行接地源瞬变电磁近源测量.在实际应用中,通过把观测的二次感应电压转换成磁场微分数据,再对磁场微分数据进行积分,求取B场数据,并进而求取视电阻率,取得了较好的应用效果.
TraditionallyBz/T is used to obtain the apparent resistivity values in the data interpretation of grounded-wire source transient electromagnetic,which lead to the problem of multiple solutions.This problem can be solved by using Bz,because the Bz with resistivity change is a single-valued function.In order to study the all-time apparent resistivity of the grounded-wire source transient electromagnetic based on B filed,the relationship between Bz orBz/T response and resistivity is analysis firstly,then all-time apparent resistivity is calculated using method of dichotomy for different geoelectric models.It is concluded that all-time apparent resistivity has high resolution reflects and only responds to the geoelectric structure well.all-time apparent resistivity has nothing to do with the offset.So,it is easy for us to observed the near soure data.Field data,the induced voltage is transformed to the B field data using the integral,then,calculated the all-time apparent resistivity is up to the real geoelectric structure.
TraditionallyBz/T is used to obtain the apparent resistivity values in the data interpretation of grounded-wire source transient electromagnetic,which lead to the problem of multiple solutions.This problem can be solved by using Bz,because the Bz with resistivity change is a single-valued function.In order to study the all-time apparent resistivity of the grounded-wire source transient electromagnetic based on B filed,the relationship between Bz orBz/T response and resistivity is analysis firstly,then all-time apparent resistivity is calculated using method of dichotomy for different geoelectric models.It is concluded that all-time apparent resistivity has high resolution reflects and only responds to the geoelectric structure well.all-time apparent resistivity has nothing to do with the offset.So,it is easy for us to observed the near soure data.Field data,the induced voltage is transformed to the B field data using the integral,then,calculated the all-time apparent resistivity is up to the real geoelectric structure.
引文
Annison C.2004.B-Field TEM data acquisition for nickel exploration[C].//ASEG 17th Geophysical Conference and Exhibition-Extended Abstracts.
Bai D H,Meju M A,Lu J,et al.2003.Numerical calculation of alltime apparent resistivity for the central loop transient electromagnetic method[J].Chinese J.Geophys.(in Chinese),46(5):697-704,doi:10.3321/j.issn:0001-5733.2003.05.018.
Chen Q Li,Yan L J,Fu Z H.2009.Algorithm for the all-Time apparent resistivity of LOTEM method[J].Chinese Journal of Engineering Geophysics(in Chinese),6(4):390-394.
Chen X D,Zhao Y,Zhang J,et al.2012.The application of HTc SQUID magnetometer to TEM[J].Chinese J.Geophys.(in Chinese),55(2):702-708,doi:10.6038/j.issn.0001-5733.2012.02.034.
Foley C P,Leslie K E.1998.Potential use of high Tc SQUIDS for airborne electromagnetics[J].Exploration Geophysics,29(2):30-34.
Huang L J.1995.A study of interpretation method for ondimensional full-region apparent resistivity of electric dipole source transient electromagnetic sounding[J].Geophysical and Geochemical Exploration(in Chinese),19(5):391-397.
Li J S,Piao H R.1993.On the calculation and application of the apparent resistivity of transient EM sounding by electric dipole[J].Computing Techniques for Geophysical and Geochemical Exploration(in Chinese),15(2):108-116.
Li W Y,Yan C W.2013.Bisection algorithm of central loop time domain transient electromagnetism all time apparent resistivity[J].Journal of Kunming University of Science and Technology(Natural Science Edition)(in Chinese),38(2):26-33.
Liu J,Yang H Y,Lin Y.2012.Binary search algorithm of all-time apparent resistivity in underground transient electromagnetic method[J].Chinese Journal of Engineering Geophysics(in Chinese),9(5):526-530.
Purss M B J,Cull J P.2003.B-field probes for down hole magnetometric resistivity surveys[J].Exploration Geophysics,34(4):233-240.
Sheng Y.1986.A single apparent resistivity expression for longoffset transient electromagnetics[J].Geophysics,51(6):1291-1297.
Tang J T,He J S.1994.A new method to define the full-zone resistivity in horizontal electric dipole frequency sounding on a layered earth[J].Acta Geophysica Sinica(in Chinese),37(4):543-552.
Wang H J,Luo Y Z.2003.Algorithm of a 2.5-Dimensional finite element method for transient electromagnetic with a central loop[J].Chinese J.Geophys.(in Chinese),46(6):855-862,doi:10.3321/j.issn:0001-5733.2003.06.020.
Wang H J.2008.Time domain transient electromagnetism all time apparent resistivity translation algorithm[J].Chinese J.Geophys.(in Chinese),51(6):1936-1942,doi:10.3321/j.issn:0001-5733.2008.06.037.
Yan L J,Hu W B,Chen Q L,et al.1999.The estimation and fast inversion of all-time apparent resistivities in long-offset transient electromagnetic sounding[J].Oil Geophysical Prospecting(in Chinese),34(5):532-538.
Yang H Y,Deng J Z,Zhang H,et al.2010.Research on full space apparent resistivity interpretation technique in mine transient electromagnetic method[J].Chinese J.Geophys.(in Chinese),53(3):651-656.
Yin C C,Piao H R.1991.A study of the definition of apparent resistivity in electromagnetic sounding[J].Geophysics and Geochemical Exploration(in Chinese),15(4):290-299.
Zhang C F,Weng A H,Sun S D,et al.2009.Computation of whole-time apparent resistivity of large rectangular loop[J].Journal of Jilin University(Earth Science Edition)(in Chinese),39(4):755-758.
陈明生,田小波.1999a.电偶源瞬变电磁测深研究(四)---瞬变电磁测深视电阻率[J].煤田地质与勘探,27(4):52-54.
陈明生,田小波.1999b.电偶源瞬变电磁测深研究(五)---实测感应电压转换成垂直磁场[J].煤田地质与勘探,27(5):63-64.
陈清礼,严良俊,付志红.2009.均匀半空间长偏移距TEM法全区视电阻率的数值计算方法[J].工程地球物理学报,6(4):390-394.
陈晓东,赵毅,张杰,等.2012.高温超导磁强计在瞬变电磁法中的应用研究[J].地球物理学报,55(2):702-708,doi:10.6038/j.issn.0001-5733.2012.02.034.
黄力军.1995.电偶源瞬变电磁测深一维全区视电阻率解释方法研究[J].物探与化探,19(5):391-397.
李文尧,晏冲为.2013.中心回线瞬变电磁法全期视电阻率的二分法求解[J].昆明理工大学学报(自然科学版),38(2):26-33.
刘建新,佟铁钢,刘春明,等.2013.E-Eφ广域视电阻率定义的改进方法及场特性识别[J].中国有色金属学报,23(9):2359-2364.
刘俊,杨海燕,林云.2012.地下瞬变电磁法全区视电阻率二分搜索算法[J].工程地球物理学报,9(5):526-530.
朴化荣.1990.电磁测深法原理[M].北京:地质出版社.
汤井田,何继善.1994.水平电偶源频率测深中全区视电阻率定义的新方法[J].地球物理学报,37(4):543-552.
提云生,刘长青,夏春妹,等.2011.激发极化法在河北隆化县多金属矿区的应用[J].西部探矿工程,23(5):153-156.
王华军,罗延钟.2003.中心回线瞬变电磁法2.5维有限单元算法[J].地球物理学报,46(6):855-862,doi:10.3321/j.issn:0001-5733.2003.06.020.
王华军.2008.时间域瞬变电磁法全区视电阻率的平移算法[J].地球物理学报,51(6):1936-1942,doi:10.3321/j.issn:0001-5733.2008.06.037.
严良俊,胡文宝,陈清礼,等.1999.长偏移距瞬变电磁测深的全区视电阻率求取及快速反演方法[J].石油地球物理勘探,34(5):532-538.
杨海燕,邓居智,张华,等.2010.矿井瞬变电磁法全空间视电阻率解释方法研究[J].地球物理学报,53(3):651-656.
白登海,Meju M A,卢健,等.2003.时间域瞬变电磁法中心方式全程视电阻率的数值计算[J].地球物理学报,46(5):698-704,doi:10.3321/j.issn:0001-5733.2003.05.018.
殷长春,朴华荣.1991.电磁测深法视电阻率定义问题的研究[J].物探与化探,15(4):290-299.
张成范,翁爱华,孙世栋,等.2009.计算矩形大定源回线瞬变电磁测深全区视电阻率[J].吉林大学学报(地球科学版),39(4):755-758.
Bai D H,Meju M A,Lu J,et al.2003.Numerical calculation of alltime apparent resistivity for the central loop transient electromagnetic method[J].Chinese J.Geophys.(in Chinese),46(5):697-704,doi:10.3321/j.issn:0001-5733.2003.05.018.
Chen Q Li,Yan L J,Fu Z H.2009.Algorithm for the all-Time apparent resistivity of LOTEM method[J].Chinese Journal of Engineering Geophysics(in Chinese),6(4):390-394.
Chen X D,Zhao Y,Zhang J,et al.2012.The application of HTc SQUID magnetometer to TEM[J].Chinese J.Geophys.(in Chinese),55(2):702-708,doi:10.6038/j.issn.0001-5733.2012.02.034.
Foley C P,Leslie K E.1998.Potential use of high Tc SQUIDS for airborne electromagnetics[J].Exploration Geophysics,29(2):30-34.
Huang L J.1995.A study of interpretation method for ondimensional full-region apparent resistivity of electric dipole source transient electromagnetic sounding[J].Geophysical and Geochemical Exploration(in Chinese),19(5):391-397.
Li J S,Piao H R.1993.On the calculation and application of the apparent resistivity of transient EM sounding by electric dipole[J].Computing Techniques for Geophysical and Geochemical Exploration(in Chinese),15(2):108-116.
Li W Y,Yan C W.2013.Bisection algorithm of central loop time domain transient electromagnetism all time apparent resistivity[J].Journal of Kunming University of Science and Technology(Natural Science Edition)(in Chinese),38(2):26-33.
Liu J,Yang H Y,Lin Y.2012.Binary search algorithm of all-time apparent resistivity in underground transient electromagnetic method[J].Chinese Journal of Engineering Geophysics(in Chinese),9(5):526-530.
Purss M B J,Cull J P.2003.B-field probes for down hole magnetometric resistivity surveys[J].Exploration Geophysics,34(4):233-240.
Sheng Y.1986.A single apparent resistivity expression for longoffset transient electromagnetics[J].Geophysics,51(6):1291-1297.
Tang J T,He J S.1994.A new method to define the full-zone resistivity in horizontal electric dipole frequency sounding on a layered earth[J].Acta Geophysica Sinica(in Chinese),37(4):543-552.
Wang H J,Luo Y Z.2003.Algorithm of a 2.5-Dimensional finite element method for transient electromagnetic with a central loop[J].Chinese J.Geophys.(in Chinese),46(6):855-862,doi:10.3321/j.issn:0001-5733.2003.06.020.
Wang H J.2008.Time domain transient electromagnetism all time apparent resistivity translation algorithm[J].Chinese J.Geophys.(in Chinese),51(6):1936-1942,doi:10.3321/j.issn:0001-5733.2008.06.037.
Yan L J,Hu W B,Chen Q L,et al.1999.The estimation and fast inversion of all-time apparent resistivities in long-offset transient electromagnetic sounding[J].Oil Geophysical Prospecting(in Chinese),34(5):532-538.
Yang H Y,Deng J Z,Zhang H,et al.2010.Research on full space apparent resistivity interpretation technique in mine transient electromagnetic method[J].Chinese J.Geophys.(in Chinese),53(3):651-656.
Yin C C,Piao H R.1991.A study of the definition of apparent resistivity in electromagnetic sounding[J].Geophysics and Geochemical Exploration(in Chinese),15(4):290-299.
Zhang C F,Weng A H,Sun S D,et al.2009.Computation of whole-time apparent resistivity of large rectangular loop[J].Journal of Jilin University(Earth Science Edition)(in Chinese),39(4):755-758.
陈明生,田小波.1999a.电偶源瞬变电磁测深研究(四)---瞬变电磁测深视电阻率[J].煤田地质与勘探,27(4):52-54.
陈明生,田小波.1999b.电偶源瞬变电磁测深研究(五)---实测感应电压转换成垂直磁场[J].煤田地质与勘探,27(5):63-64.
陈清礼,严良俊,付志红.2009.均匀半空间长偏移距TEM法全区视电阻率的数值计算方法[J].工程地球物理学报,6(4):390-394.
陈晓东,赵毅,张杰,等.2012.高温超导磁强计在瞬变电磁法中的应用研究[J].地球物理学报,55(2):702-708,doi:10.6038/j.issn.0001-5733.2012.02.034.
黄力军.1995.电偶源瞬变电磁测深一维全区视电阻率解释方法研究[J].物探与化探,19(5):391-397.
李文尧,晏冲为.2013.中心回线瞬变电磁法全期视电阻率的二分法求解[J].昆明理工大学学报(自然科学版),38(2):26-33.
刘建新,佟铁钢,刘春明,等.2013.E-Eφ广域视电阻率定义的改进方法及场特性识别[J].中国有色金属学报,23(9):2359-2364.
刘俊,杨海燕,林云.2012.地下瞬变电磁法全区视电阻率二分搜索算法[J].工程地球物理学报,9(5):526-530.
朴化荣.1990.电磁测深法原理[M].北京:地质出版社.
汤井田,何继善.1994.水平电偶源频率测深中全区视电阻率定义的新方法[J].地球物理学报,37(4):543-552.
提云生,刘长青,夏春妹,等.2011.激发极化法在河北隆化县多金属矿区的应用[J].西部探矿工程,23(5):153-156.
王华军,罗延钟.2003.中心回线瞬变电磁法2.5维有限单元算法[J].地球物理学报,46(6):855-862,doi:10.3321/j.issn:0001-5733.2003.06.020.
王华军.2008.时间域瞬变电磁法全区视电阻率的平移算法[J].地球物理学报,51(6):1936-1942,doi:10.3321/j.issn:0001-5733.2008.06.037.
严良俊,胡文宝,陈清礼,等.1999.长偏移距瞬变电磁测深的全区视电阻率求取及快速反演方法[J].石油地球物理勘探,34(5):532-538.
杨海燕,邓居智,张华,等.2010.矿井瞬变电磁法全空间视电阻率解释方法研究[J].地球物理学报,53(3):651-656.
白登海,Meju M A,卢健,等.2003.时间域瞬变电磁法中心方式全程视电阻率的数值计算[J].地球物理学报,46(5):698-704,doi:10.3321/j.issn:0001-5733.2003.05.018.
殷长春,朴华荣.1991.电磁测深法视电阻率定义问题的研究[J].物探与化探,15(4):290-299.
张成范,翁爱华,孙世栋,等.2009.计算矩形大定源回线瞬变电磁测深全区视电阻率[J].吉林大学学报(地球科学版),39(4):755-758.