激电中梯装置测量长方异常体的三维正演模拟及启示
|
摘要
激电扫面实际测量时,由于测线布设方位和地下地质异常体的走向不尽一致,从而导致测量结果不同,给物探异常的解释带来了困难,甚至会引起误导,因此这里利用数值模拟方法,计算了激电中间梯度法测量条件下二度异常体在不同方位测量条件下的三维响应规律。首先实现了有限单元法三维正演的公式推导及其数值解求解;其次通过一维模型计算了数值解和解析解的误差,验证了本文中理论计算方法的有效性;最后设计了三维模型进行了理论计算,发现了激电中间梯度装置条件下二度异常体的走向和测量剖面方向成不同夹角时的面上响应规律。研究结果表明:随着两者夹角的增加,激电异常体的规模和形态逐渐与设计异常体一致,且两者夹角在60°~90°之间探测的实测数据较为可信。
When IP(induced polarization method)scanning is applied to actual measurement,the direction of the lines' layout is not always in accordance with the orientation of geological anomaly and thus the different results in that measurement cause difficulty in the interpretation of geophysical anomaly and even misleading.Therefore the author utilizes the numerical simulation method to calculate the three-dimensional response of two-dimensional body under different azimuth's measurement when central gradient array method is used in IP.At first,the formula derivation and implementation steps of three-dimensional forward model using the finite element method are fulfilled;Secondly through the one-dimensional model to calculate the error of numerical solution and analytical solution,the validity of theoretical calculation method in this paper is verified;At last,Designing three-dimensional model for theoretical calculation,the author found response's law when the orientation of two-dimensional body and the direction of measurement's profile have different angles under the condition of central gradient array method in the IP.The results reveal that,with the increase of the angle between the orientation of two-dimensional body and the direction of measurement's profile,the size and shape of geological anomaly is consistent with the abnormal body designed in the modeling gradually and data in actual measurement that the angle of that two is between 60 and 90 is more reliable.
When IP(induced polarization method)scanning is applied to actual measurement,the direction of the lines' layout is not always in accordance with the orientation of geological anomaly and thus the different results in that measurement cause difficulty in the interpretation of geophysical anomaly and even misleading.Therefore the author utilizes the numerical simulation method to calculate the three-dimensional response of two-dimensional body under different azimuth's measurement when central gradient array method is used in IP.At first,the formula derivation and implementation steps of three-dimensional forward model using the finite element method are fulfilled;Secondly through the one-dimensional model to calculate the error of numerical solution and analytical solution,the validity of theoretical calculation method in this paper is verified;At last,Designing three-dimensional model for theoretical calculation,the author found response's law when the orientation of two-dimensional body and the direction of measurement's profile have different angles under the condition of central gradient array method in the IP.The results reveal that,with the increase of the angle between the orientation of two-dimensional body and the direction of measurement's profile,the size and shape of geological anomaly is consistent with the abnormal body designed in the modeling gradually and data in actual measurement that the angle of that two is between 60 and 90 is more reliable.
引文
[1] F.UMI.The induced polarization methed in ore investigation[J].Geophysical Prospecting,1961,9(3):459-477.
[2]余官娣.激电法寻找隐伏多金属矿床[J].地质与勘探,1989,9(2):41-43YU G D.Inquisition of inclusive polymetallic deposits by IP method[J].Geology and Prospecting,1989,9(2):41-43.(In Chinese)
[3]徐笑民.双频激电法的找矿效果[J].中国有色金属学报,1993,3(4):84-87.XU X M.The prospecting effect of dual-frequency IP method[J].The Chinese Journal of Nonferrous Metals,1993,3(4):84-87.(In Chinese)
[4]黄颖坚.应用双频激电偶极剖面法追索隐伏构造带[J].中国有色金属学报,1994,4(2):1-3.HUANG Y J.Application of dual frequency induction dipole profile method to reconnaissance of concealed structural Belt[J].The Chinese Journal of Nonferrous Metals,1994,4(2):1-3.(In Chinese)
[5]何继善.双频激电法[M].北京:高等教育出版社,2006.HE J S.Dual frequency induced polarization Method[M].Beijing:High Education Press,2006.(In Chinese)
[6]佟铁钢.双频道激电法的研究与在地质找矿中的应用[D].长沙:中南大学,2004.TONG T G.Study on dual-channel IP method and its application in geological prospecting[D].Changsha:Central South University,2004.(In Chinese)
[7]柳建新.西部特殊地貌景观区双频激电法方法及应用研究[D].长沙:中南大学,2006.LIU J X.Dual-frequency IP method and its application in special landscapes of west China[D].Changsha:Central South University,2006.(In Chinese)
[8] VAUDELET P,SCHMUTZ M,PESSEL M,et al.Mapping of contaminant plumes with geo-electrical methods.A case study in urban context[J].Journal of Applied Geophysics,2011,75(4):738-751.
[9] OLOWOFELA J A,AKINYEMI O D,OGUNGBE A S.Imaging and detecting underground contaminants in landfill sites using electrical impedance tomography(EIT):A case study of Lagos,southwestern,Nigeria[J].Research Journal of Environmental and Earth Sciences,2012,4(3):270-281.
[10]中华人民共和国国土资源部.时间域激发极化法技术规程[S].北京:全国国土资源标准化技术委员会,2016:1-27.Ministry of Land and Resources of the People's Republic of China.Technical specifications for Induced polarization in time domain[S].Beijing:National Land and ResourcesStandardizationTechnicalCommittee,2016:1-27.(In Chinese)
[11]禚冯昌.稠油油藏蒸汽驱三维物理模拟及数值模拟研究[D].东营:中国石油大学(华东),2011.ZHUO F C.Three-dimensional physical simulation and numerical simulation of steam flooding in heavy oil reservoir[D].Dongying:China University of Petroleum(East China),2011.(In Chinese)
[12]ZHANG J,MACKIE R L,MADDEN T R.3-D resistivity forward modeling and inversion using conjugate gradients[J].Geophysics,1995,60(5):1313-1325.
[13]刘小军,王家林.低电阻率差情况下三维激电法快速数值模拟[J].物探化探计算技术,2006,28(3):194-196.LIU X J,WANG J L.Fast numerical simulation of 3D induced polarization for low resistivity[J].Computing Techniques for Geophysical and Geochemical Exploration,2006,28(3):194-196.(In Chinese)
[14]刘海飞,阮百尧,柳建新,等.起伏地形电导率连续变化的三维激电数据有限元数值模拟[J].物探化探计算技术,2008,30(4):308-313LIU H F,RUAN B Y,LIU J X,et al.Finite element numerical simulation of three dimensional IP data with continuous variation of electrical conductivity in undulating terrain[J].Computing Techniques for Geophysical and Geochemical Exploration,2008,30(4):308-313.(In Chinese)
[15]乔中林.直流电法三维正演多重网格算法研究[D].北京:中国地质大学,2007.QIAO Z L.Study of 3-D forward modeling of DC resistivity method sing multi-grid method[D].Beijing:China University of Geosciences,2007.(In Chinese)
[16]王帅.直流电法三维正演并行算法研究[D].北京:中国地质大学,2011.WANG S.Research on 3-D forward modeling paratlet algorithm of DC resistivity data[D].Beijing:China University of Geosciences,2011.(In Chinese)
[17]戴前伟,江沸菠.基于混沌振荡PSO-BP算法的电阻率层析成像非线性反演[J].中国有色金属学报,2013,23(10):2897-2904.DAI Q W,JIANG F B.Nonlinear inversion for electrical resistivity tomography based on chaotic oscillation PSO-BP algorithm[J].The Chinese Journal of Nonferrous Metals,2013,23(10):2897-2904.(In Chinese)
[18]崔益安,朱肖雄,陈志学,等.污染监测视电阻率快速三维反演成像[J].中国有色金属学报,2015,25(8):2205-2211.CUI Y A,ZHU X X,CHEN Z X,et al.Fast 3Dapparent resistivity inversion for contaminant monitoring data[J].The Chinese Journal of Nonferrous Metals,2015,25(8):2205-2211.(In Chinese)
[19]黄俊革.三维电阻率/极化率有限元正演模拟与反演成像[D].长沙:中南大学,2003.HUANG J G.Three-dimensional Resistivity/Polarization Finite Element Forward Modeling and Inversion Imaging[D].Changsha:Central South university,2003.(In Chinese)
[20]黄俊革,阮百尧,鲍光淑.基于有限单元法的三维地电断面电阻率反演[J].中南大学学报(自然科学版),2004,35(2):295-299.HUANG J G,RUAN B Y,BAO G S.Resistivity inversion on 3-D section based on FEM[J].Journal of Central South University(Natural Science),2004,35(2):295-299.(In Chinese)
[21]强建科.起伏地形三维电阻率正演模拟与反演成像研究[D].武汉:中国地质大学,2006.QIANG J K.The research on 3-D resistivity forward and inversion algorithm on undulate topography[D].Wuhan:China University of Geosciences,2006.(In Chinese)
[22]强建科,何继善.椭球体上双频激电法的正演与反演算法[J].中南大学学报,2007,38(6):1199-1205.QIANG J K,He J S.Algorithm of forward and inversion of dual-frequency induced polarization method on elliptical sphere[J].Journal of Central South University,2007,38(6):1199-1205.(In Chinese)
[23]LI SHU-CAI,NIE LI-CHAO,LIU BIN,et al.3D electrical resistivity inversion using prior spatial shape constraints[J].Applied Geophysics,2013,10(4):361-372.
[24]GASPERIKOVA E,HUBBARD S S,WATSON D B,et al.Long-term electrical resistivity monitoring of recharge-induced contaminant plume behavior[J].Journal of Contaminant Hydrology,2012,142/143(11):33-49.
[2]余官娣.激电法寻找隐伏多金属矿床[J].地质与勘探,1989,9(2):41-43YU G D.Inquisition of inclusive polymetallic deposits by IP method[J].Geology and Prospecting,1989,9(2):41-43.(In Chinese)
[3]徐笑民.双频激电法的找矿效果[J].中国有色金属学报,1993,3(4):84-87.XU X M.The prospecting effect of dual-frequency IP method[J].The Chinese Journal of Nonferrous Metals,1993,3(4):84-87.(In Chinese)
[4]黄颖坚.应用双频激电偶极剖面法追索隐伏构造带[J].中国有色金属学报,1994,4(2):1-3.HUANG Y J.Application of dual frequency induction dipole profile method to reconnaissance of concealed structural Belt[J].The Chinese Journal of Nonferrous Metals,1994,4(2):1-3.(In Chinese)
[5]何继善.双频激电法[M].北京:高等教育出版社,2006.HE J S.Dual frequency induced polarization Method[M].Beijing:High Education Press,2006.(In Chinese)
[6]佟铁钢.双频道激电法的研究与在地质找矿中的应用[D].长沙:中南大学,2004.TONG T G.Study on dual-channel IP method and its application in geological prospecting[D].Changsha:Central South University,2004.(In Chinese)
[7]柳建新.西部特殊地貌景观区双频激电法方法及应用研究[D].长沙:中南大学,2006.LIU J X.Dual-frequency IP method and its application in special landscapes of west China[D].Changsha:Central South University,2006.(In Chinese)
[8] VAUDELET P,SCHMUTZ M,PESSEL M,et al.Mapping of contaminant plumes with geo-electrical methods.A case study in urban context[J].Journal of Applied Geophysics,2011,75(4):738-751.
[9] OLOWOFELA J A,AKINYEMI O D,OGUNGBE A S.Imaging and detecting underground contaminants in landfill sites using electrical impedance tomography(EIT):A case study of Lagos,southwestern,Nigeria[J].Research Journal of Environmental and Earth Sciences,2012,4(3):270-281.
[10]中华人民共和国国土资源部.时间域激发极化法技术规程[S].北京:全国国土资源标准化技术委员会,2016:1-27.Ministry of Land and Resources of the People's Republic of China.Technical specifications for Induced polarization in time domain[S].Beijing:National Land and ResourcesStandardizationTechnicalCommittee,2016:1-27.(In Chinese)
[11]禚冯昌.稠油油藏蒸汽驱三维物理模拟及数值模拟研究[D].东营:中国石油大学(华东),2011.ZHUO F C.Three-dimensional physical simulation and numerical simulation of steam flooding in heavy oil reservoir[D].Dongying:China University of Petroleum(East China),2011.(In Chinese)
[12]ZHANG J,MACKIE R L,MADDEN T R.3-D resistivity forward modeling and inversion using conjugate gradients[J].Geophysics,1995,60(5):1313-1325.
[13]刘小军,王家林.低电阻率差情况下三维激电法快速数值模拟[J].物探化探计算技术,2006,28(3):194-196.LIU X J,WANG J L.Fast numerical simulation of 3D induced polarization for low resistivity[J].Computing Techniques for Geophysical and Geochemical Exploration,2006,28(3):194-196.(In Chinese)
[14]刘海飞,阮百尧,柳建新,等.起伏地形电导率连续变化的三维激电数据有限元数值模拟[J].物探化探计算技术,2008,30(4):308-313LIU H F,RUAN B Y,LIU J X,et al.Finite element numerical simulation of three dimensional IP data with continuous variation of electrical conductivity in undulating terrain[J].Computing Techniques for Geophysical and Geochemical Exploration,2008,30(4):308-313.(In Chinese)
[15]乔中林.直流电法三维正演多重网格算法研究[D].北京:中国地质大学,2007.QIAO Z L.Study of 3-D forward modeling of DC resistivity method sing multi-grid method[D].Beijing:China University of Geosciences,2007.(In Chinese)
[16]王帅.直流电法三维正演并行算法研究[D].北京:中国地质大学,2011.WANG S.Research on 3-D forward modeling paratlet algorithm of DC resistivity data[D].Beijing:China University of Geosciences,2011.(In Chinese)
[17]戴前伟,江沸菠.基于混沌振荡PSO-BP算法的电阻率层析成像非线性反演[J].中国有色金属学报,2013,23(10):2897-2904.DAI Q W,JIANG F B.Nonlinear inversion for electrical resistivity tomography based on chaotic oscillation PSO-BP algorithm[J].The Chinese Journal of Nonferrous Metals,2013,23(10):2897-2904.(In Chinese)
[18]崔益安,朱肖雄,陈志学,等.污染监测视电阻率快速三维反演成像[J].中国有色金属学报,2015,25(8):2205-2211.CUI Y A,ZHU X X,CHEN Z X,et al.Fast 3Dapparent resistivity inversion for contaminant monitoring data[J].The Chinese Journal of Nonferrous Metals,2015,25(8):2205-2211.(In Chinese)
[19]黄俊革.三维电阻率/极化率有限元正演模拟与反演成像[D].长沙:中南大学,2003.HUANG J G.Three-dimensional Resistivity/Polarization Finite Element Forward Modeling and Inversion Imaging[D].Changsha:Central South university,2003.(In Chinese)
[20]黄俊革,阮百尧,鲍光淑.基于有限单元法的三维地电断面电阻率反演[J].中南大学学报(自然科学版),2004,35(2):295-299.HUANG J G,RUAN B Y,BAO G S.Resistivity inversion on 3-D section based on FEM[J].Journal of Central South University(Natural Science),2004,35(2):295-299.(In Chinese)
[21]强建科.起伏地形三维电阻率正演模拟与反演成像研究[D].武汉:中国地质大学,2006.QIANG J K.The research on 3-D resistivity forward and inversion algorithm on undulate topography[D].Wuhan:China University of Geosciences,2006.(In Chinese)
[22]强建科,何继善.椭球体上双频激电法的正演与反演算法[J].中南大学学报,2007,38(6):1199-1205.QIANG J K,He J S.Algorithm of forward and inversion of dual-frequency induced polarization method on elliptical sphere[J].Journal of Central South University,2007,38(6):1199-1205.(In Chinese)
[23]LI SHU-CAI,NIE LI-CHAO,LIU BIN,et al.3D electrical resistivity inversion using prior spatial shape constraints[J].Applied Geophysics,2013,10(4):361-372.
[24]GASPERIKOVA E,HUBBARD S S,WATSON D B,et al.Long-term electrical resistivity monitoring of recharge-induced contaminant plume behavior[J].Journal of Contaminant Hydrology,2012,142/143(11):33-49.