基于自适应有限元的复电阻率法2.5维复杂构造电磁场模拟
|
摘要
复电阻率法在矿产、油气勘探调查中发挥着重要作用.为了认识复杂构造的复电阻率法电磁场的变化规律,本文基于自适应有限元方法,采用非结构化网格,引入Cole-Cole模型,实现了电偶源2.5D复电阻率法电磁场正演,可以模拟复杂地形和地电结构,正演结果更符合野外实际地质情况.通过将本文的计算结果与半空间模型解析解、层状介质和起伏模型结果进行对比,验证了本文算法的正确性.最后,基于复杂地电模型,通过正演模拟,系统分析了地形、激电参数、复杂构造对复电阻率法电磁场的影响特征.
The complex resistivity method plays an important role in exploration of minerals,oil and gas.In order to better understand the variation law of the electromagnetic field with complex structure for the complex resistivity method,we introduce the Cole-Cole model into the electromagnetic equation.Then we use the adaptive finite element method based on an unstructured grid to simulate the responses of electric dipole 2.5 Dcomplex resistivity.We adopt this method,because it can simulate the complex terrain and underground structure,and the forward modeling results are more consistent with actual geological conditions.We compare our results with both the analytical results and those from the literature.It is verified that the algorithm is accurate.Based on the complex geoelectric model,we analyze the influence of terrain,induced polarization parameters and geoelectric structure on the electromagnetic field in the complex resistivity method by forward modeling.
The complex resistivity method plays an important role in exploration of minerals,oil and gas.In order to better understand the variation law of the electromagnetic field with complex structure for the complex resistivity method,we introduce the Cole-Cole model into the electromagnetic equation.Then we use the adaptive finite element method based on an unstructured grid to simulate the responses of electric dipole 2.5 Dcomplex resistivity.We adopt this method,because it can simulate the complex terrain and underground structure,and the forward modeling results are more consistent with actual geological conditions.We compare our results with both the analytical results and those from the literature.It is verified that the algorithm is accurate.Based on the complex geoelectric model,we analyze the influence of terrain,induced polarization parameters and geoelectric structure on the electromagnetic field in the complex resistivity method by forward modeling.
引文
Cai J T,Ruan B,Zhao G Z.2007.Two-dimensional modeling of complex resistivity using finite element method.Chinese Journal of Geophysics(in Chinese),50(6):1869-1876.
Di Q Y,Wang R.2008.Forward and Inversion Method of Controlled Source Audio-Frequency MagnetoTelluric Data and its Application(in Chinese).Beijing:Science Press.
Farias V,Maranhao C,Rocha B,et al.2010.Induced polarization forward modelling using finite element method and the fractal model.Applied Mathematical Modelling,34(7):1849-1860.
Fan C S.2013.Research on complex resistivity forward and inversion with finite element method and itsapplication[Ph.D.thesis](in Chinese).Changchun:Jilin University.
He Z X,Wang X B.2007.Geo-electrical anomalous pattern of reservoir and oil/gas detection by electromagnetic survey.Oil Geophysical Prospecting(in Chinese),42(1):102-106.
Key K,Weiss C.2006.Adaptive finite-element modeling using unstructured grids:The 2Dmagnetotelluric example.Geophysics,71(6):G291-G299.
Key K,Ovall J.2011.A parallel goal-oriented adaptive finite element method for 2.5Delectromagnetic modelling.Geophysical Journal International,186(1):137-154.
Li Y,Key K.2007.2D marine controlled-source electromagnetic modeling:Part 1-An adaptive finite-element algorithm.Geophysics,72(2):WA51-WA62.
Li Y,Lin P R,Li T L.2011a.Finite element method for solving anomalous complex potential of 2.5-D complex resistivity.Journal of Jilin University(Earth Science Edition)(in Chinese),41(5):1596-1604.
Li Y,Lin P R,Xiao Y,et al.2011b.Induced Polarization effect on frequency-domain electromagnetic sounding with electric dipole source.Chinese Journal of Geophysics(in Chinese),54(7):1935-1944,doi:10.396/j.issn.0001-5733.2011.07.028.
Luo Y Z,Zhang G Q.1988.The Theory of Frequency-Domain Induced Polarization Method(in Chinese).Beijing:Geological Publishing House.
Mitsuhata Y.2000.2-D electromagnetic modeling by finite-element method with a dipole source and topography.Geophysics,65(2):465475.
Pelton W,Ward S,Hallof P.1978.Mineral discrimination and removal of inductive coupling with multifrequency IP.Geophysics,43(3):588-609.
Ren Z Y,Tang J T.2014.A goal-oriented adaptive finite-element approach for multi-electrode resistivity system.Geophys.J.Int.,199(1):136-145.
Ulrich C,Slater L.2004.Induced polarization measurements on unsaturated,unconsolidated sands.Geophysics,69(3):762-771.
Veeken P,Legeydo P,Davidenko Y.2009.Benefits of the induced polarization geoelectric method to hydrocarbon exploration.Geophysics,74(2):B47-B59.
Wang D Y,Li T L,Li J P,et al.2010.Forward simulation of a 3Dcomplex resistivity model.Progress in Geophysics(in Chinese),25(1):266-271,doi:10.3969/j.issn.1004-2903.2010.01.035.
Ye Y X,Li Y G,Deng J Z,et al.2014.2.5Dinduced polarization forward modeling using the adaptive finite-element method.Applied Geophysics,11(4):500-507.
Yin C C,Zhang B,Liu Y H,et al.2015.2.5Dforward modeling of the time-domain airborne EM system in areas with topographic relief.Chinese Journal of Geophysics(in Chinese),58(4):1411-1424,doi:10.6038/cjg20150427.
Zhang H,Li T L,Dong R X.2006.Modeling electromagnetic responses of complex resistivity 3-D body using volume integral equation method.Coal Geology&Exploration(in Chinese),34(1):70-73.
Zhang S Z,Li Y X,Zhang S C,et al.1984.The low frequency electrical phase spectra of mineralizedrocks(ores)and some factors which influence them in some sulfideore deposits in China.Chinese Journal of Geophysics(in Chinese),27(2):176-189.
蔡军涛,阮百尧,赵国泽等.2007.复电阻率法二维有限元数值模拟.地球物理学报,50(6):1869-1876.
底青云,王若.2008.可控源音频大地电磁数据正反演及方法应用.北京:科学出版社.
范翠松.2013.基于有限元法的复电阻率正反演研究及应用[博士论文].长春:吉林大学.
何展翔,王绪本.2007.油气藏电性异常模式与电磁法油气检测.石油地球物理勘探,42(1):102-106.
李勇,林品荣,李桐林等.2011a.基于异常复电位2.5维CR有限元数值模拟.吉林大学学报(地球科学版),41(5):1596-1604.
李勇,林品荣,肖原等.2011b.电偶源频率电磁测深激发极化效应研究.地球物理学报,54(7):1935-1944,doi:10.396/j.issn.0001-5733.2011.07.028.
罗延钟,张桂青.1988.频率域激电法原理.北京:地质出版社.
王大勇,李桐林,李建平等.2010.三维复电阻率模型电磁场正演模拟研究.地球物理学进展,25(1):266-271,doi:10.3969/j.issn.1004-2903.2010.01.035.
殷长春,张博,刘云鹤等.2015.2.5维起伏地表条件下时间域航空电磁正演模拟.地球物理学报,58(4):1411-1424,doi:10.6038/cjg20150427.
张辉,李桐林,董瑞霞.2006.体积分方程法模拟复电阻率三维体电磁响应.煤田地质与勘探,34(1):70-73.
张赛珍,李英贤,张树椿等.1984.我国几个金属矿区岩(矿)石的低频电相位频率特性及其影响因素.地球物理学报,27(2):176-189.
Di Q Y,Wang R.2008.Forward and Inversion Method of Controlled Source Audio-Frequency MagnetoTelluric Data and its Application(in Chinese).Beijing:Science Press.
Farias V,Maranhao C,Rocha B,et al.2010.Induced polarization forward modelling using finite element method and the fractal model.Applied Mathematical Modelling,34(7):1849-1860.
Fan C S.2013.Research on complex resistivity forward and inversion with finite element method and itsapplication[Ph.D.thesis](in Chinese).Changchun:Jilin University.
He Z X,Wang X B.2007.Geo-electrical anomalous pattern of reservoir and oil/gas detection by electromagnetic survey.Oil Geophysical Prospecting(in Chinese),42(1):102-106.
Key K,Weiss C.2006.Adaptive finite-element modeling using unstructured grids:The 2Dmagnetotelluric example.Geophysics,71(6):G291-G299.
Key K,Ovall J.2011.A parallel goal-oriented adaptive finite element method for 2.5Delectromagnetic modelling.Geophysical Journal International,186(1):137-154.
Li Y,Key K.2007.2D marine controlled-source electromagnetic modeling:Part 1-An adaptive finite-element algorithm.Geophysics,72(2):WA51-WA62.
Li Y,Lin P R,Li T L.2011a.Finite element method for solving anomalous complex potential of 2.5-D complex resistivity.Journal of Jilin University(Earth Science Edition)(in Chinese),41(5):1596-1604.
Li Y,Lin P R,Xiao Y,et al.2011b.Induced Polarization effect on frequency-domain electromagnetic sounding with electric dipole source.Chinese Journal of Geophysics(in Chinese),54(7):1935-1944,doi:10.396/j.issn.0001-5733.2011.07.028.
Luo Y Z,Zhang G Q.1988.The Theory of Frequency-Domain Induced Polarization Method(in Chinese).Beijing:Geological Publishing House.
Mitsuhata Y.2000.2-D electromagnetic modeling by finite-element method with a dipole source and topography.Geophysics,65(2):465475.
Pelton W,Ward S,Hallof P.1978.Mineral discrimination and removal of inductive coupling with multifrequency IP.Geophysics,43(3):588-609.
Ren Z Y,Tang J T.2014.A goal-oriented adaptive finite-element approach for multi-electrode resistivity system.Geophys.J.Int.,199(1):136-145.
Ulrich C,Slater L.2004.Induced polarization measurements on unsaturated,unconsolidated sands.Geophysics,69(3):762-771.
Veeken P,Legeydo P,Davidenko Y.2009.Benefits of the induced polarization geoelectric method to hydrocarbon exploration.Geophysics,74(2):B47-B59.
Wang D Y,Li T L,Li J P,et al.2010.Forward simulation of a 3Dcomplex resistivity model.Progress in Geophysics(in Chinese),25(1):266-271,doi:10.3969/j.issn.1004-2903.2010.01.035.
Ye Y X,Li Y G,Deng J Z,et al.2014.2.5Dinduced polarization forward modeling using the adaptive finite-element method.Applied Geophysics,11(4):500-507.
Yin C C,Zhang B,Liu Y H,et al.2015.2.5Dforward modeling of the time-domain airborne EM system in areas with topographic relief.Chinese Journal of Geophysics(in Chinese),58(4):1411-1424,doi:10.6038/cjg20150427.
Zhang H,Li T L,Dong R X.2006.Modeling electromagnetic responses of complex resistivity 3-D body using volume integral equation method.Coal Geology&Exploration(in Chinese),34(1):70-73.
Zhang S Z,Li Y X,Zhang S C,et al.1984.The low frequency electrical phase spectra of mineralizedrocks(ores)and some factors which influence them in some sulfideore deposits in China.Chinese Journal of Geophysics(in Chinese),27(2):176-189.
蔡军涛,阮百尧,赵国泽等.2007.复电阻率法二维有限元数值模拟.地球物理学报,50(6):1869-1876.
底青云,王若.2008.可控源音频大地电磁数据正反演及方法应用.北京:科学出版社.
范翠松.2013.基于有限元法的复电阻率正反演研究及应用[博士论文].长春:吉林大学.
何展翔,王绪本.2007.油气藏电性异常模式与电磁法油气检测.石油地球物理勘探,42(1):102-106.
李勇,林品荣,李桐林等.2011a.基于异常复电位2.5维CR有限元数值模拟.吉林大学学报(地球科学版),41(5):1596-1604.
李勇,林品荣,肖原等.2011b.电偶源频率电磁测深激发极化效应研究.地球物理学报,54(7):1935-1944,doi:10.396/j.issn.0001-5733.2011.07.028.
罗延钟,张桂青.1988.频率域激电法原理.北京:地质出版社.
王大勇,李桐林,李建平等.2010.三维复电阻率模型电磁场正演模拟研究.地球物理学进展,25(1):266-271,doi:10.3969/j.issn.1004-2903.2010.01.035.
殷长春,张博,刘云鹤等.2015.2.5维起伏地表条件下时间域航空电磁正演模拟.地球物理学报,58(4):1411-1424,doi:10.6038/cjg20150427.
张辉,李桐林,董瑞霞.2006.体积分方程法模拟复电阻率三维体电磁响应.煤田地质与勘探,34(1):70-73.
张赛珍,李英贤,张树椿等.1984.我国几个金属矿区岩(矿)石的低频电相位频率特性及其影响因素.地球物理学报,27(2):176-189.
© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号 地址：北京市海淀区学院路29号 邮编：100083 电话：办公室：(+86 10)66554848；文献借阅、咨询服务、科技查新：66554700 |