(1. 中南大學 地球科學與信息物理學院, 長沙 410083；
2. 中南大學 有色資源與地質(zhì)災害探查湖南省重點實驗室, 長沙 410083)
摘 要: 雙頻激電法是一種非常有效的探測極化目標體的勘探方法，可以同時觀測獲取視電阻率和極化率數(shù)據(jù)。為了提高雙頻激電觀測數(shù)據(jù)解譯的可靠性，利用粒子群優(yōu)化算法對電阻率數(shù)據(jù)和極化率進行聯(lián)合反演。在對異常目標體采用旋轉(zhuǎn)橢球體進行幾何近似模擬的基礎上，基于Core-Core散射理論實現(xiàn)雙頻激電法的快速數(shù)值模擬計算，為粒子群優(yōu)化聯(lián)合反演提供正演基礎。在對粒子群優(yōu)化算法參數(shù)進行分析設定的基礎上設計了聯(lián)合反演算法，并采用加入不同程度噪聲的模擬雙頻激電數(shù)據(jù)對算法進行了實驗測試。測試結(jié)果表明：粒子群優(yōu)化聯(lián)合反演算法能有效實現(xiàn)電阻率與極化率數(shù)據(jù)的聯(lián)合反演，具有很強的抗噪聲能力，且算法收斂快、穩(wěn)定性好。進一步的實測數(shù)據(jù)聯(lián)合反演測試還表明：該算法具有較低的模型初始信息依賴性，在給定較大搜索空間的條件下仍能反演出較為滿意的結(jié)果模型，具有較好的實用性。
關(guān)鍵字: 聯(lián)合反演；雙頻激電；粒子群優(yōu)化算法；視電阻率；極化率
（1. School of Geosciences and Info-Physics, Central South University, Changsha 410083, China;
2. Hunan Key Laboratory of Non-ferrous Resources and Geological Hazard Detection,
Central South University, Changsha 410083, China）
Abstract:Dual frequency IP method is very effective to detect polarizable targets even in mountain terrains. Both resistivity and polarizability data can be observed during dual frequency IP measurement. In order to interpret dual frequency IP data quantitatively and more accurately, a joint inversion was proposed based on particle swarm optimization (PSO). Through simulating anomalies by ellipsoids, the forward of dual frequency IP can be given based on Core-Core theory, which provides the foundation to implement joint inversion by PSO. After selecting appropriate PSO algorithm parameters, the joint inversion is designed meticulously. Then synthetic data with noise in different degrees is used to test the joint inversion algorithm. The results of testing show that the joint inversion is effective to invert dual frequency IP data with high tolerance to noise and fast convergence. Further field data tests demonstrate that the PSO inversion is not sensitive to the initial model information. It can get a satisfactory result model even in a large given search space. That is very meaningful for interpreting field data in dual frequency IP observation.
Key words: joint inversion; dual frequency IP; particle swarm optimization; resistivity; polarizability
