直流电阻率法与大地电磁法的二维联合反演
|
摘要
本文对电性联合反演进行了深入研究,以减少地球物理反演的多解性.将直流电阻率(DC)与大地电磁(MT)数据加入到同一反演数据集中.引入Tikhonov正则化思想建立反演目标函数,使反演过程更加高效稳定.在解决正则化反演问题过程中,分别采用了二阶最大平滑稳定因子和改进的L-curve法,提高了反演结果的稳定性和正则化因子的求取精度;最后运用非线性共轭梯度法(NLCG)对反演目标函数实现最优化求解.经研究表明:联合反演方法与单一反演方法相比,能够更加有效的约束反演模型范围;反演算法快速稳定,提高了反演精度,减少了对地下地质结构认识的模糊性.
In order to reduce the solutions of geophysical inversion,the joint inversion theory on dielectric properties were studied. The direct current resistivity( DC) and magnetotelluric( MT) data is added to the same inversion of data set. this paper introduces the idea of Tikhonov regularization establish the inversion objective function,which makes the inversion process more efficient and stable. In the process to solve the major problems of regularized inversion,We using the maximum smoothing stable second-order factor and improved L-curve method chosen regularization factor. To improve the stability of inversion result and improve the accuracy for the regularization factor; Finally,through the nonlinear conjugate gradient method( NLCG) to realize the optimization of the objective function for solving inversion. By the test shows that the joint inversion inversion is more effective constraint model range than a single inversion method; The inversion algorithm is fast and stable,improve the inversion precision and reduce the fuzziness of understanding of underground structures.
In order to reduce the solutions of geophysical inversion,the joint inversion theory on dielectric properties were studied. The direct current resistivity( DC) and magnetotelluric( MT) data is added to the same inversion of data set. this paper introduces the idea of Tikhonov regularization establish the inversion objective function,which makes the inversion process more efficient and stable. In the process to solve the major problems of regularized inversion,We using the maximum smoothing stable second-order factor and improved L-curve method chosen regularization factor. To improve the stability of inversion result and improve the accuracy for the regularization factor; Finally,through the nonlinear conjugate gradient method( NLCG) to realize the optimization of the objective function for solving inversion. By the test shows that the joint inversion inversion is more effective constraint model range than a single inversion method; The inversion algorithm is fast and stable,improve the inversion precision and reduce the fuzziness of understanding of underground structures.
引文
Bosch M,Meza R,Jiménez R,et al.2006.Joint gravity and magnetic inversion in 3D using Monte Carlo methods[J].Geophysics,71(4):G153-G156.
Chen H G,Li J X,Wu J S,et al.2012.Study on simulated-annealing MT-gravity joint inversion[J].Chinese Journal of Geophysics(in Chinese),55(2):663-670,doi:10.6038/j.issn.0001-5733.2012.02.030.
Chen J,Wen N,Chen B Y.2007.Joint inversion of gravity-magneticelectrical-seismic combination survey:Progress and prospect[J].Progress in Geophysics(in Chinese),22(5):1427-1438,doi:10.3969/j.issn.1004-2903.2007.05.013.
Chen L H,Sun J G,Wu Y G,et al.2002.Review of quasi-Linear approximation in geophysical inversion[J].Progress in Geophysics(in Chinese),17(3):464-472,doi:10.3969/j.issn.1004-2903.2002.03.016.
Chen X,Yu P,Zhang L L,et al.2010.Regularized synchronous joint inversion of MT and seismic data[J].Seismology and Geology(in Chinese),32(3):402-408.
Cheng S H.2012.Contrastive analysis of different static shift correction methods in megnetotelluric(MT)(in Chinese)[Ph.D.thesis].Xi’an:Chang’an University.
Dong H B,Wang C L.2003.Development and application of 2Dresistivity imaging surveys[J].Earth Science Frontiers(in Chinese),10(1):171-176.
Hering A,Misiek R,Gyulai A,et al.1995.A joint inversion algorithm to process geoelectric and sutface wave seismic data.Part I:Basic ideas[J].Geophysical Prospecting,43(2):135-156.
Huang Z H,Di Q Y,Hou S L.2006.CSAMT static correction and its application[J].Progress in Geophysics(in Chinese),21(4):1290-1295,doi:10.3969/j.issn.1004-2903.2006.04.038.
Jing R Z,Bao G S,Chen S Q.2003.A review of the researches for geophysical combinative inversion[J].Progress in Geophysics(in Chinese),18(3):535-540,doi:10.3969/j.issn.1004-2903.2003.03.033.
Li M,Lin W D.2014.The study of regularization inversion for 2.5dimension DC resistivity based on minimum support stabilizing factor[J].Journal of East China Institute of Technology(in Chinese),37(3):292-298.
Liu Y,Lv Q T,Meng G X,et al.2012.Joint electromagnetic and seismic inversion survey:status and prospect[J].Progress in Geophysics(in Chinese),27(6):2444-2451,doi:10.6038/j.issn.1004-2903.2012.06.019.
Moorkamp M,Heincke B,Jegen M,et al.2011.A framework for 3-Djoint inversion of MT,gravity and seismic refraction data[J].Geophysical Journal International,184(1):477-493,doi:10.1111/j.1365-246X.2010.04856.x.
Peng M,Tan H D,Jiang M,et al.2013.Three-dimensional joint inversion of magnetotelluric and seismic travel time data with crossgradient constraints[J].Chinese Journal of Geophysics(in Chinese),56(8):2728-2738,doi:10.6038/cjg20130821.
Sasaki Y.1989.Two-dimensional joint inversion of magnetotelluric and dipole-dipole resistivity data[J].Geophysics,54(2):254-262.
Sharma S P,Kaikkonen P.1999.Appraisal of equivalence and suppression problems in 1D EM and DC measurements using global optimization and joint inversion[J].Geophysical Prospecting,47(2):219-249.
Sharma S P,Verma S K.2011.Solutions of the inherent problem of the equivalence in direct current resistivity and electromagnetic methods through global optimization and joint inversion by successive refinement of model space[J].Geophysical Prospecting,59(4):760-776,doi:10.1111/j.1365-2478.2011.00952.x.
Tikhonov A N,Arsenin Y.1977.Solution of Ill-Posed Problems[M].Washington:Winston&Sons.
Vozoff K,Jupp D L B.1975.Joint inversion of geophysical data[J].Geophysical Journal of the Royal Astronomical Society,42(3):977-991.
Wan L,Lin T T,Lin J,et al.2013.Joint inversion of MRS and TEMdata based on adaptive genetic algorithm[J].Chinese Journal of Geophysics(in Chinese),56(11):3728-3740,doi:10.6038/cjg20131114.
Wei W B.2002.New advance and prospect of magnetotelluric sounding(MT)in china[J].Progress in Geophysics(in Chinese),17(2):245-254,doi:10.3969/j.issn.1004-2903.2002.02.009.
Yang H,Dai S K,Song H B,et al.2002.Overview of joint inversion of integrated geophysics[J].Progress in Geophysics(in Chinese),17(2):262-271,doi:10.3969/j.issn.1004-2903.2002.02.011.
Yang Z W,Yan J Y,Liu Y,et al.2012.Research progresses of the high-density resistivity method[J].Geology and Exploration,(in Chinese),48(5):969-978.
Yu P,Wang J L,Wu J S,et al.2006.Review and discussions on geophysical joint inversion[J].Progress in Exploration Geophysics(in Chinese),29(2):87-93,134.
Yu P,Wang J L,Wu J S,et al.2007.Comparative analysis of inversion methods of retrieving atmospheric profiles with GPS occultation measurements[J].Chinese Journal of Geophysics(in Chinese),50(2):529-538,doi:10.3321/j.issn:0001-5733.2007.02.026.
Zeyen H,Pous J.1993.3-D joint inversion of magnetic and gravimetric data with a priori information[J].Geophysical Journal International,112(2):244-256.
Zhu J Y.2003.Overview of integrated geophysics of oil&gas exploration[J].Progress in Geophysics(in Chinese),18(1):19-23,doi:10.3969/j.issn.1004-2903.2003.01.003.
陈华根,李嘉虓,吴健生,等.2012.MT-重力模拟退火联合反演研究[J].地球物理学报,55(2):663-670,doi:10.6038/j.issn.0001-5733.2012.02.030.
陈洁,温宁,陈邦彦.2007.重磁电震联合反演研究进展与展望[J].地球物理学进展,22(5):1427-1438,doi:10.3969/j.issn.1004-2903.2007.05.013.
陈丽虹,孙建国,吴燕冈,等.2002.地球物理反演的拟线性近似方法综述[J].地球物理学进展,17(3):464-472,doi:10.3969/j.issn.1004-2903.2002.03.016.
陈晓,于鹏,张罗磊,等.2010.大地电磁与地震正则化同步联合反演[J].地震地质,32(3):402-408.
程少华.2012.大地电磁静态效应校正方法对比研究[硕士论文].西安:长安大学.
董浩斌,王传雷.2003.高密度电法的发展与应用[J].地学前缘,10(1):171-176.
黄兆辉,底青云,侯胜利.2006.CSAMT的静态效应校正及应用[J].地球物理学进展,21(4):1290-1295,doi:10.3969/j.issn.1004-2903.2006.04.038.
敬荣中,鲍光淑,陈绍裘.2003.地球物理联合反演研究综述[J].地球物理学进展,18(3):535-540,doi:10.3969/j.issn.1004-2903.2003.03.033.
李曼,林文东.2014.基于最小支持的2.5维直流电阻率正则化反演研究[J].东华理工大学学报(自然科学版),37(3):292-298.
刘彦,吕庆田,孟贵祥,等.2012.大地电磁与地震联合反演研究现状与展望[J].地球物理学进展,27(6):2444-2451,doi:10.6038/j.issn.1004-2903.2012.06.019.
彭淼,谭捍东,姜枚,等.2013.基于交叉梯度耦合的大地电磁与地震走时资料三维联合反演[J].地球物理学报,56(8):2728-2738,doi:10.6038/cjg20130821.
万玲,林婷婷,林君,等.2013.基于自适应遗传算法的MRS-TEM联合反演方法研究[J].地球物理学报,56(11):3728-3740,doi:10.6038/cjg20131114.
魏文博.2002.我国大地电磁测深新进展及瞻望[J].地球物理学进展,17(2):245-254,doi:10.3969/j.issn.1004-2903.2002.02.009.
杨辉,戴世坤,宋海斌,等.2002.综合地球物理联合反演综述[J].地球物理学进展,17(2):262-271,doi:10.3969/j.issn.1004-2903.2002.02.011.
杨振威,严加永,刘彦,等.2012.高密度电阻率法研究进展[J].地质与勘探,48(5):969-978.
于鹏,王家林,吴健生,等.2006.地球物理联合反演的研究现状和分析[J].勘探地球物理进展,29(2):87-93,134.
于鹏,王家林,吴健生,等.2007.重力与地震资料的模拟退火约束联合反演[J].地球物理学报,50(2):529-538,doi:10.3321/j.issn:0001-5733.2007.02.026.
祝靓谊.2003.油气勘探综合地球物理研究方法综述[J].地球物理学进展,18(1):19-23,doi:10.3969/j.issn.1004-2903.2003.01.003.
Chen H G,Li J X,Wu J S,et al.2012.Study on simulated-annealing MT-gravity joint inversion[J].Chinese Journal of Geophysics(in Chinese),55(2):663-670,doi:10.6038/j.issn.0001-5733.2012.02.030.
Chen J,Wen N,Chen B Y.2007.Joint inversion of gravity-magneticelectrical-seismic combination survey:Progress and prospect[J].Progress in Geophysics(in Chinese),22(5):1427-1438,doi:10.3969/j.issn.1004-2903.2007.05.013.
Chen L H,Sun J G,Wu Y G,et al.2002.Review of quasi-Linear approximation in geophysical inversion[J].Progress in Geophysics(in Chinese),17(3):464-472,doi:10.3969/j.issn.1004-2903.2002.03.016.
Chen X,Yu P,Zhang L L,et al.2010.Regularized synchronous joint inversion of MT and seismic data[J].Seismology and Geology(in Chinese),32(3):402-408.
Cheng S H.2012.Contrastive analysis of different static shift correction methods in megnetotelluric(MT)(in Chinese)[Ph.D.thesis].Xi’an:Chang’an University.
Dong H B,Wang C L.2003.Development and application of 2Dresistivity imaging surveys[J].Earth Science Frontiers(in Chinese),10(1):171-176.
Hering A,Misiek R,Gyulai A,et al.1995.A joint inversion algorithm to process geoelectric and sutface wave seismic data.Part I:Basic ideas[J].Geophysical Prospecting,43(2):135-156.
Huang Z H,Di Q Y,Hou S L.2006.CSAMT static correction and its application[J].Progress in Geophysics(in Chinese),21(4):1290-1295,doi:10.3969/j.issn.1004-2903.2006.04.038.
Jing R Z,Bao G S,Chen S Q.2003.A review of the researches for geophysical combinative inversion[J].Progress in Geophysics(in Chinese),18(3):535-540,doi:10.3969/j.issn.1004-2903.2003.03.033.
Li M,Lin W D.2014.The study of regularization inversion for 2.5dimension DC resistivity based on minimum support stabilizing factor[J].Journal of East China Institute of Technology(in Chinese),37(3):292-298.
Liu Y,Lv Q T,Meng G X,et al.2012.Joint electromagnetic and seismic inversion survey:status and prospect[J].Progress in Geophysics(in Chinese),27(6):2444-2451,doi:10.6038/j.issn.1004-2903.2012.06.019.
Moorkamp M,Heincke B,Jegen M,et al.2011.A framework for 3-Djoint inversion of MT,gravity and seismic refraction data[J].Geophysical Journal International,184(1):477-493,doi:10.1111/j.1365-246X.2010.04856.x.
Peng M,Tan H D,Jiang M,et al.2013.Three-dimensional joint inversion of magnetotelluric and seismic travel time data with crossgradient constraints[J].Chinese Journal of Geophysics(in Chinese),56(8):2728-2738,doi:10.6038/cjg20130821.
Sasaki Y.1989.Two-dimensional joint inversion of magnetotelluric and dipole-dipole resistivity data[J].Geophysics,54(2):254-262.
Sharma S P,Kaikkonen P.1999.Appraisal of equivalence and suppression problems in 1D EM and DC measurements using global optimization and joint inversion[J].Geophysical Prospecting,47(2):219-249.
Sharma S P,Verma S K.2011.Solutions of the inherent problem of the equivalence in direct current resistivity and electromagnetic methods through global optimization and joint inversion by successive refinement of model space[J].Geophysical Prospecting,59(4):760-776,doi:10.1111/j.1365-2478.2011.00952.x.
Tikhonov A N,Arsenin Y.1977.Solution of Ill-Posed Problems[M].Washington:Winston&Sons.
Vozoff K,Jupp D L B.1975.Joint inversion of geophysical data[J].Geophysical Journal of the Royal Astronomical Society,42(3):977-991.
Wan L,Lin T T,Lin J,et al.2013.Joint inversion of MRS and TEMdata based on adaptive genetic algorithm[J].Chinese Journal of Geophysics(in Chinese),56(11):3728-3740,doi:10.6038/cjg20131114.
Wei W B.2002.New advance and prospect of magnetotelluric sounding(MT)in china[J].Progress in Geophysics(in Chinese),17(2):245-254,doi:10.3969/j.issn.1004-2903.2002.02.009.
Yang H,Dai S K,Song H B,et al.2002.Overview of joint inversion of integrated geophysics[J].Progress in Geophysics(in Chinese),17(2):262-271,doi:10.3969/j.issn.1004-2903.2002.02.011.
Yang Z W,Yan J Y,Liu Y,et al.2012.Research progresses of the high-density resistivity method[J].Geology and Exploration,(in Chinese),48(5):969-978.
Yu P,Wang J L,Wu J S,et al.2006.Review and discussions on geophysical joint inversion[J].Progress in Exploration Geophysics(in Chinese),29(2):87-93,134.
Yu P,Wang J L,Wu J S,et al.2007.Comparative analysis of inversion methods of retrieving atmospheric profiles with GPS occultation measurements[J].Chinese Journal of Geophysics(in Chinese),50(2):529-538,doi:10.3321/j.issn:0001-5733.2007.02.026.
Zeyen H,Pous J.1993.3-D joint inversion of magnetic and gravimetric data with a priori information[J].Geophysical Journal International,112(2):244-256.
Zhu J Y.2003.Overview of integrated geophysics of oil&gas exploration[J].Progress in Geophysics(in Chinese),18(1):19-23,doi:10.3969/j.issn.1004-2903.2003.01.003.
陈华根,李嘉虓,吴健生,等.2012.MT-重力模拟退火联合反演研究[J].地球物理学报,55(2):663-670,doi:10.6038/j.issn.0001-5733.2012.02.030.
陈洁,温宁,陈邦彦.2007.重磁电震联合反演研究进展与展望[J].地球物理学进展,22(5):1427-1438,doi:10.3969/j.issn.1004-2903.2007.05.013.
陈丽虹,孙建国,吴燕冈,等.2002.地球物理反演的拟线性近似方法综述[J].地球物理学进展,17(3):464-472,doi:10.3969/j.issn.1004-2903.2002.03.016.
陈晓,于鹏,张罗磊,等.2010.大地电磁与地震正则化同步联合反演[J].地震地质,32(3):402-408.
程少华.2012.大地电磁静态效应校正方法对比研究[硕士论文].西安:长安大学.
董浩斌,王传雷.2003.高密度电法的发展与应用[J].地学前缘,10(1):171-176.
黄兆辉,底青云,侯胜利.2006.CSAMT的静态效应校正及应用[J].地球物理学进展,21(4):1290-1295,doi:10.3969/j.issn.1004-2903.2006.04.038.
敬荣中,鲍光淑,陈绍裘.2003.地球物理联合反演研究综述[J].地球物理学进展,18(3):535-540,doi:10.3969/j.issn.1004-2903.2003.03.033.
李曼,林文东.2014.基于最小支持的2.5维直流电阻率正则化反演研究[J].东华理工大学学报(自然科学版),37(3):292-298.
刘彦,吕庆田,孟贵祥,等.2012.大地电磁与地震联合反演研究现状与展望[J].地球物理学进展,27(6):2444-2451,doi:10.6038/j.issn.1004-2903.2012.06.019.
彭淼,谭捍东,姜枚,等.2013.基于交叉梯度耦合的大地电磁与地震走时资料三维联合反演[J].地球物理学报,56(8):2728-2738,doi:10.6038/cjg20130821.
万玲,林婷婷,林君,等.2013.基于自适应遗传算法的MRS-TEM联合反演方法研究[J].地球物理学报,56(11):3728-3740,doi:10.6038/cjg20131114.
魏文博.2002.我国大地电磁测深新进展及瞻望[J].地球物理学进展,17(2):245-254,doi:10.3969/j.issn.1004-2903.2002.02.009.
杨辉,戴世坤,宋海斌,等.2002.综合地球物理联合反演综述[J].地球物理学进展,17(2):262-271,doi:10.3969/j.issn.1004-2903.2002.02.011.
杨振威,严加永,刘彦,等.2012.高密度电阻率法研究进展[J].地质与勘探,48(5):969-978.
于鹏,王家林,吴健生,等.2006.地球物理联合反演的研究现状和分析[J].勘探地球物理进展,29(2):87-93,134.
于鹏,王家林,吴健生,等.2007.重力与地震资料的模拟退火约束联合反演[J].地球物理学报,50(2):529-538,doi:10.3321/j.issn:0001-5733.2007.02.026.
祝靓谊.2003.油气勘探综合地球物理研究方法综述[J].地球物理学进展,18(1):19-23,doi:10.3969/j.issn.1004-2903.2003.01.003.