本文選題：頻譜激電法 + Cole-Cole模型��； 參考：《吉林大學》2017年碩士論文

【摘要】：頻譜激電法(SIP)也稱復電阻率法(CR)具有多參數(shù)解釋異常的能力,多個Cole-Cole模型參數(shù)可以為評價異常源的性質(zhì)提供更多的信息。由于Cole-Cole模型參數(shù)較多,各參數(shù)靈敏度不同,且SIP中既包含激發(fā)極化效應又包含電磁效應,這些都制約著SIP反演問題的研究。反演是解釋地電結(jié)構(gòu)的最終手段,研究一種基于同時考慮激發(fā)極化效應和電磁感應效應正演、同時反演四個Cole-Cole模型參數(shù)的SIP反演算法對SIP的研究至關重要。本文首先基于既考慮激電效應又考慮電磁效應的一維SIP正演,實現(xiàn)了一維SIP最小二乘反演。本文采用利用SIP數(shù)據(jù)實部和虛部聯(lián)合反演的方法,將反演域從復數(shù)域變?yōu)閷崝?shù)域,降低了反演的多解性。為了更好的反演靈敏度較低的頻率相關系數(shù)和時間常數(shù),本文提出實部和虛部加權(quán)反演的思想,在定義目標函數(shù)時對實部和虛部的目標函數(shù)進行加權(quán),在反演不同參數(shù)時,通過控制實部和虛部在反演中所占的比重來提高各參數(shù)的反演效果。為此本文推導了利用SIP數(shù)據(jù)實部和虛部聯(lián)合加權(quán)反演Cole-Cole模型參數(shù)的最小二乘反演算法。模型試算結(jié)果表明:這種反演方法不僅可以很好的反演零頻電阻率和極化率,還可以顯著提高頻率相關系數(shù)和時間常數(shù)的反演效果。三維反演方面,本文基于既考慮激電效應又考慮電磁效應的三維SIP有限差分正演,利用現(xiàn)今較為成熟的NLCG技術(shù),實現(xiàn)了同時反演四個Cole-Cole模型參數(shù)的三維SIP反演算法。三維反演采用電場實部和虛部聯(lián)合反演的方案,來降低反演的多解性。由于是電場實部和虛部聯(lián)合反演,故本文推導了三維反演電場實部和虛部目標函數(shù)的靈敏度矩陣表達式,并給出利用擬正演求取實部和虛部靈敏度矩陣的方法,大大降低了計算費用,避免了直接求靈敏度矩陣耗時過長和內(nèi)存不夠等問題。為了解決一些參數(shù)靈敏度低、反演難度大的問題,在反演時采用對實部和虛部加權(quán)的方法,并設計理論模型進行了試算,證明了算法的可行性、優(yōu)越性和穩(wěn)定性。
[Abstract]:Spectral IP (complex resistivity method) has the ability to interpret anomalies with multiple parameters. Multiple Cole-Cole model parameters can provide more information for evaluating the properties of anomaly sources. Because of the large number of parameters in the Cole-Cole model, the sensitivity of each parameter is different, and the SIP includes both polarization effect and electromagnetic effect, which restrict the research of SIP inversion. Inversion is the ultimate method to explain geoelectric structure. It is very important for SIP to study a SIP inversion algorithm based on considering both induced polarization effect and electromagnetic induction effect and simultaneously retrieving the parameters of four Cole-Cole models. In this paper, the one-dimensional SIP least square inversion is realized based on the one-dimensional SIP forward modeling which considers both the IP effect and the electromagnetic effect. In this paper, by using the method of joint inversion of real and imaginary parts of SIP data, the inversion domain is changed from complex number domain to real number domain, which reduces the multiplicity of inversion solutions. In order to improve the frequency correlation coefficient and time constant with lower inversion sensitivity, the idea of weighted inversion of real and imaginary parts is presented in this paper. The objective functions of real and imaginary parts are weighted when defining objective functions. The inversion effect of each parameter is improved by controlling the proportion of real part and imaginary part in inversion. In this paper, the least square inversion algorithm for the parameters of Cole-Cole model is derived by combining real and virtual parts of SIP data. The simulation results show that this inversion method can not only invert the zero frequency resistivity and polarizability but also improve the inversion effect of frequency correlation coefficient and time constant. In the aspect of 3D inversion, based on the 3D SIP finite difference forward modeling which considers both the IP effect and the electromagnetic effect, a 3D SIP inversion algorithm for simultaneous inversion of four Cole-Cole model parameters is implemented by using the current mature NLCG technique. In order to reduce the multiplicity of inversion, the method of electric field real part and virtual part is used in 3D inversion. Because of the joint inversion of the real and imaginary parts of the electric field, this paper deduces the expression of the sensitivity matrix of the objective function of the real part and the imaginary part of the three-dimensional inversion of the electric field, and gives a method to obtain the sensitivity matrix of the real part and the imaginary part by using the pseudo forward modeling. The computational cost is greatly reduced, and the problems such as time-consuming and insufficient memory are avoided. In order to solve the problems of low parameter sensitivity and high inversion difficulty, the method of weighted real and imaginary parts is used in inversion, and a theoretical model is designed to test the feasibility, superiority and stability of the algorithm.
【學位授予單位】：吉林大學
【學位級別】：碩士
【學位授予年份】：2017
【分類號】：P631.3

【相似文獻】

相關期刊論文 前10條

1 熊翔葆;麥也爾電弧數(shù)學模型參數(shù)的數(shù)值計算及麥氏模型的校正[J];華東交通大學學報;1989年02期

2 黃夢橋,羅漢;回歸模型參數(shù)的距離估計方法[J];數(shù)學理論與應用;2001年02期

3 李抗美;預測油田產(chǎn)量的伽馬模型求解方法——與《求解Γ模型參數(shù)的新方法》作者商榷[J];石油勘探與開發(fā);2002年04期

4 朱悅新;我國人口模型參數(shù)辨識[J];西安交通大學學報;1982年02期

5 何s鉉，

本文編號：1791230