地质学报
地質學報
지질학보
ACTA GEOLOGICA SINICA
2014年
4期
755-762
,共8页
范翠松%李桐林%王大勇%赵广茂%李建平
範翠鬆%李桐林%王大勇%趙廣茂%李建平
범취송%리동림%왕대용%조엄무%리건평
起伏地形%有限元法%复电阻率%2.5维反演%OCCAM反演
起伏地形%有限元法%複電阻率%2.5維反縯%OCCAM反縯
기복지형%유한원법%복전조솔%2.5유반연%OCCAM반연
topography%finite element method%complex resistivity%2 .5D Inversion%OCCAM inversion
复电阻率法以岩、矿石的频谱参数特性为勘探依据,多种参数的组合解释能为评价异常源的性质提供更为丰富的信息。但其反演问题一直没有得到较好地解决,已经严重地影响和制约了该方法的应用和发展。本文由频率域电磁场的偏微分方程出发,结合Cole-Cole数学模型进行公式推导,并利用有限单元法实现了复电阻率2.5维正演模拟,能够计算起伏地形条件下包含电磁效应的复电阻率响应。在此基础上,提出了利用多个排列的视电阻率和视相位数据的复电阻率反演方法,并在反演方程中加入了Occam法的光滑模型约束,能够对复电阻率参数直接进行大规模的全区反演,有效减弱了SIP反演中欠定性严重的问题。通过理论模型的反演算例表明,该方法能够快速、准确地反演出二维地质断面上所有单元的4种复电阻率参数,从而验证了该反演算法的正确性与稳定性。
複電阻率法以巖、礦石的頻譜參數特性為勘探依據,多種參數的組閤解釋能為評價異常源的性質提供更為豐富的信息。但其反縯問題一直沒有得到較好地解決,已經嚴重地影響和製約瞭該方法的應用和髮展。本文由頻率域電磁場的偏微分方程齣髮,結閤Cole-Cole數學模型進行公式推導,併利用有限單元法實現瞭複電阻率2.5維正縯模擬,能夠計算起伏地形條件下包含電磁效應的複電阻率響應。在此基礎上,提齣瞭利用多箇排列的視電阻率和視相位數據的複電阻率反縯方法,併在反縯方程中加入瞭Occam法的光滑模型約束,能夠對複電阻率參數直接進行大規模的全區反縯,有效減弱瞭SIP反縯中欠定性嚴重的問題。通過理論模型的反縯算例錶明,該方法能夠快速、準確地反縯齣二維地質斷麵上所有單元的4種複電阻率參數,從而驗證瞭該反縯算法的正確性與穩定性。
복전조솔법이암、광석적빈보삼수특성위감탐의거,다충삼수적조합해석능위평개이상원적성질제공경위봉부적신식。단기반연문제일직몰유득도교호지해결,이경엄중지영향화제약료해방법적응용화발전。본문유빈솔역전자장적편미분방정출발,결합Cole-Cole수학모형진행공식추도,병이용유한단원법실현료복전조솔2.5유정연모의,능구계산기복지형조건하포함전자효응적복전조솔향응。재차기출상,제출료이용다개배렬적시전조솔화시상위수거적복전조솔반연방법,병재반연방정중가입료Occam법적광활모형약속,능구대복전조솔삼수직접진행대규모적전구반연,유효감약료SIP반연중흠정성엄중적문제。통과이론모형적반연산례표명,해방법능구쾌속、준학지반연출이유지질단면상소유단원적4충복전조솔삼수,종이험증료해반연산법적정학성여은정성。
Complex resistivity method is based on the characteristics of spectral parameters of rocks and ores, and the multiple spectral parameters provide more information for evaluating anomalous bodies.However, the inverse problem of SIP has not been resolved well,so that it has restricted the application and development of SIP.Therefore,we derive the formula for the 2.5D SIP responses,which starting from maxwell’s equations combined with the Cole-Cole model,then we complete the forward modeling with finite element method.The forward method can simulate SIP response contains the electromagnetic effects with topography.On this basis,we propose a combined inversion method by using multi-arrangement of apparent resistivity and apparent phase of SIP date,and the smooth model constraint of Occam method is used to improve the stability of inversion.The algorithm can perform large-scale inversion,to deal with largely underdetermined inversion problem.The numerical test show that this method can be used to inverse all SIP parameters of geological model with synthetic data (noiseless data or noise data),which indicates the inversion algorithm is feasible.