有限單元法（Finite Element Method, FEM）的優(yōu)點是適用于模擬起伏不規(guī)則的地質(zhì)體和地形起伏直流電阻率法為人工源供電,在電源附近電場衰減快,變化梯度大,模擬結(jié)果在電源位置附近出現(xiàn)誤差變大現(xiàn)象,模擬得到的網(wǎng)格節(jié)點上的電位值已不能近似表示電源附近電場實際變化特征,即點源奇異問題。針對電源奇異問題,在正演過程中利用奇異邊界校正算法對點源奇異影響進行校正,改善模擬結(jié)果在電源附近的精度,同時降低邊界條件的影響。從穩(wěn)定電流場滿足的基本微分方程出發(fā),推導了點源電場邊值問題對應的變分問題。離散積分區(qū)域時,以四邊形為最小研究單元。正演過程中求解剛度矩陣線性方程組得到的網(wǎng)格節(jié)點電位值為對應于特定波數(shù)的付氏變換電位值,需進行反付氏變換從而求取空間域電位值,從波數(shù)域電位值求取空間域電位值,利用最優(yōu)化離散波數(shù)進行反付氏變換。用最小二乘平滑約束方法進行電阻率數(shù)據(jù)二維反演,正則化因子的取值直接影響到反演結(jié)果的穩(wěn)定性及準確性,取值大反演結(jié)果模型光滑,不能反映小的構(gòu)造信息,取值大反演結(jié)果模型粗糙,可能引入一些不真實的地質(zhì)構(gòu)造信息。在反演迭代過程中,用主動約束平衡方法計算正則化因子值,根據(jù)模型參數(shù)分辨... 

【文章來源】：中南大學湖南省 211工程院校 985工程院校 教育部直屬院校

【文章頁數(shù)】：55 頁

【學位級別】：碩士

【部分圖文】：

變寬帶壓縮存儲示意圖（左幅：網(wǎng)格，右幅：系數(shù)矩陣，黑點代表非零數(shù)據(jù)）

的看出該校正算法的作用，剖分區(qū)域不外延，水平位置±20m為網(wǎng)格剖分的邊界，矩形剖分單元尺寸大小Imxlm。模擬結(jié)果如圖2-4所示。15 1—‘ ‘ ■■‘ ‘ |1 80 I I I ‘. ? - ? NoBC/singularity correction I ?? NoBC/singularitycorrection . ^i, . BC/slngularity 70" . BC/singularity correction |..10 "ft t - ； /^， ： 6。-: i-I ‘ \：.......I ： yso-i 丨 丨 I rtrr，”.；、.丨 \ IB 0 MF ； 5 40- --， i- 1 r"% 丨丨 丨.、/ 丨,..........； i—.i....-召 Currert Electrodes， V^ -5- i ；“V f. ? i ?* i! 20- V — i ；' -■ i /■10 ； ‘ 10 - Current Electrodes -.15 i i—LU 0：——：-20 -10 0 10 20 -20 -10 ^ 0 ▼ 10 20X-position(m) X-position(m)圖2- 4數(shù)值數(shù)值模擬結(jié)果曲線圖(紅色代表未利用邊界/奇異校正，綠色代表利用邊界/奇異校正。左圖：電位值曲線；右圖：相對誤差曲線）從圖2-4看出，隨著觀測點離剖分區(qū)域邊界距離越近，模擬結(jié)果精度降低，利用邊界/奇異校正處理后，數(shù)值模擬結(jié)果精度得到改善。圖2表明通過利用邊界/奇異校正可以有效改善電源奇異性問題

型的視電阻率有解析解。計算所采用的網(wǎng)格剖分參數(shù)：剖分數(shù)目個數(shù)46x52，網(wǎng)格大小不等，模型尺寸大小為1090mx3200m。圖2-7為雙頻激電對稱四級測深裝置模擬的視電阻率數(shù)值解與解析解曲線對比圖，從表格2H可以看出，與沒有經(jīng)過邊界/奇異校正處理后的數(shù)值解精度相比，經(jīng)過邊界/奇異校正處理后的數(shù)值解具有更高的精度。A3 A2 A1 M 0 N B1 B2 B3Pi=100n-m hi=12mP2=500fi ? m h2= 15mP3=10n.m圖2- 6三層水平層狀模型與對稱四級測深裝置示意圖AB/2(ffl) analytical BC NoBC rms_BC(%) rms_NoBC(%)5 101.1 102.9 116.97 ITTS 15.710 107.3 107.06 111.67 0.22 4.0715 118.86 117.29 126.54 1.32 6.4620 132.57 129.21 142.83 2.54 7.7330 155.21 149.25 164.9 3.84 6.2540 165

【參考文獻】：
期刊論文
[1]基于空間形態(tài)先驗約束的三維電阻率反演成像方法(英文)[J]. 李術才,聶利超,劉斌,宋杰,劉征宇,蘇茂鑫,徐磊.  Applied Geophysics. 2013(04)
[2]基于不等式約束的最小二乘法三維電阻率反演及其算法優(yōu)化[J]. 劉斌,李術才,李樹忱,聶利超,鐘世航,李利平,宋杰,劉征宇.  地球物理學報. 2012(01)
[3]基于非結(jié)構(gòu)化網(wǎng)格的2.5-D直流電阻率自適應有限元數(shù)值模擬（英文）[J]. 湯井田,王飛燕,任政勇.  地球物理學報. 2010(03)
[4]垂直激電測深二維自適應正則化反演[J]. 劉海飛,柳建新,阮百堯,張賽民.  同濟大學學報(自然科學版). 2009(12)
[5]變阻尼共軛梯度算法及其性能分析[J]. 劉海飛,阮百堯,柳建新.  地球物理學進展. 2008(01)
[6]三維地形直流電阻率有限元法模擬[J]. 強建科,羅延鐘.  地球物理學報. 2007(05)
[7]電阻率二維神經(jīng)網(wǎng)絡反演[J]. 徐海浪,吳小平.  地球物理學報. 2006(02)
[8]大地電磁自適應正則化反演算法[J]. 陳小斌,趙國澤,湯吉,詹艷,王繼軍.  地球物理學報. 2005(04)
[9]計算最優(yōu)化離散波數(shù)的優(yōu)化算法[J]. 柳建新,劉海飛.  物探化探計算技術. 2005(01)
[10]復雜地形條件下直流電阻率異常三維數(shù)值模擬研究[J]. 熊彬,阮百堯,羅延鐘.  地質(zhì)與勘探. 2003(04)



本文編號：3498048