【摘要】：直流電阻率法常常被工程勘察所采用,然而工區(qū)復(fù)雜的地質(zhì)地形環(huán)境和嚴格的工程要求都對直流電法提出新的挑戰(zhàn),特別是在后期的資料解釋環(huán)節(jié)和反演成像環(huán)節(jié),高低起伏的地形給解釋和成像帶來錯誤解譯和假異常等問題。起伏地形對直流電法的影響情況和影響程度,起伏地形起伏度的變化引起的直流電法響應(yīng)特征的變化等都是工程勘察中亟待解決的問題。正演模擬的算法步驟較為繁雜,本文選取網(wǎng)格剖分和傅里葉反余弦變換做為研究重點。起伏地形的數(shù)值模擬需要非結(jié)構(gòu)化網(wǎng)格剖分技術(shù),文中討論并編程實現(xiàn)經(jīng)典的三角單元網(wǎng)格剖分算法;反余弦變換而言,因為它關(guān)聯(lián)著正演算法的精度,本文主要對反變換中的最優(yōu)化波數(shù)電極距范圍進行討論,進而改進正演模擬算法。具體來說,網(wǎng)格剖分使用限定Delaunay三角剖分算法,算法只需要提供限定條件,即地電模型的地形邊界、人工邊界和模型邊界等,就能夠自動得到地電模型的三角剖分結(jié)果,同時保證剖分結(jié)果中各條邊界一定存在,又可根據(jù)人為需要加密某部分區(qū)域,使起伏地形能夠進行正演模擬研究和反演成像研究。本文的傅里葉反余弦變換研究以最優(yōu)化波數(shù)算法為基礎(chǔ),通過對雙層水平介質(zhì)模型和垂直接觸帶模型的正反余弦變換解析公式的研究,得到計算離散波數(shù)的電極距范圍應(yīng)該根據(jù)地下空間電性變化和勘探深度要求進行設(shè)計,再采用最優(yōu)化算法算得離散波數(shù)和變換系數(shù),經(jīng)對比分析證明上述計算離散波數(shù)的方法是正確的,提高數(shù)值模擬的精度,保證正演算法的有效性,為反演成像奠定基礎(chǔ)。文中第二部分利用正演算法研究起伏地形的高密度電法視電阻率響應(yīng)特征,然后分析地形起伏度變化對響應(yīng)特征的影響區(qū)域和程度變化。最后,研討反演成像算法并設(shè)計理論模型驗證,最后將帶地形實測資料進行成像解釋。詳細工作包括采用不同的裝置對山脊、山谷和陡坎地形分別進行數(shù)值計算,通過對比溫納裝置、偶極裝置和微分裝置的視電阻率剖面圖,發(fā)現(xiàn)溫納裝置和微分裝置對應(yīng)剖面圖中的假異常形態(tài)相似,在山脊地形下方有低阻假異常,坡腳兩側(cè)是高阻假異常;而在山谷下方高阻假異常,坡腳兩側(cè)為低阻假異常;陡坎地形坡頂對應(yīng)低阻假異常,坡腳對應(yīng)高阻假異常,但微分裝置的視電阻率變化范圍更大。偶極裝置在不同地形中對應(yīng)的假異常和溫納裝置假異常位置對應(yīng)相同而電性相反,即溫納裝置的低阻假異常對應(yīng)偶極裝置的高阻假異常。起伏地形的響應(yīng)特征研究為視電阻率剖面的解釋提供參考。然后討論起伏度對正演模擬響應(yīng)的影響,設(shè)計不同縱向高度和橫向?qū)挾鹊纳郊鼓Ｐ?對比分析各山脊模型的視電阻率剖面圖,發(fā)現(xiàn)山脊的縱向高度會影響假異常的數(shù)值大小,而山脊的橫向?qū)挾葧绊懠佼惓５膮^(qū)域大小。起伏度對視電阻率響應(yīng)的影響,為評價和解釋地形對實測數(shù)據(jù)的影響提供理論支撐。最后部分是關(guān)于起伏地形的反演研究,首先敘述直流電阻率反演的基本理論,然后對偏導(dǎo)數(shù)矩陣的計算和穩(wěn)定迭代處理方法進行介紹,對理論低阻和高阻模型的模擬數(shù)據(jù)進行反演成像,對比理論模型和反演結(jié)果印證算法的正確性,最后對具有地形的實測數(shù)據(jù)反演,并對反演結(jié)果進行分析解釋,證明文中基于起伏地形的反演算法對帶地形的實測數(shù)據(jù)是可行的。
【學(xué)位授予單位】：成都理工大學(xué)
【學(xué)位級別】：碩士
【學(xué)位授予年份】：2018
【分類號】：P631.322
【圖文】：

中任意點集使用計算機繪制 Voronoi 圖，點集的 Delaunay 三角剖分諾圖的對偶圖。上述過程是計算幾何中重要組成部分，計算幾何包Delaunay 三角剖分的基礎(chǔ)問題，同時也是本文的基礎(chǔ)工作，如點線的單元的數(shù)據(jù)結(jié)構(gòu)優(yōu)化等。由于本文所涉及的問題為二維直流電阻率，因此在進行三角剖分的過程中將會涉及到限定點和限定邊兩因素元的尺度和質(zhì)量控制問題等，所以本章將會對 Delaunay 三角剖分用幾何理論進行簡要的介紹，同時詳細討論三角剖分的算法實現(xiàn)過程集的 Delaunay 三角剖分初的 Delaunay 三角剖分研究內(nèi)容主要集聚在點集和凸包，凸包可以一根橡皮筋全部罩住，橡皮筋所繪制的多邊形即是該點集的凸包。在，所以點集三角化邊界自動會形成在凸包的邊上。Delaunay 早就上點的維諾圖的對偶是一種點集的三角化，且所得到的三角形具有，或者成為 Delaunay 準則。

進行三角剖分，并設(shè)計局部變換法，由于該算法是通過判斷三角邊是否，并通過換邊達到目的的算法，所以該算法又稱為換邊算法。隨后的 Boatson 分別給出一種三角剖分的滿足 Delaunay 準則的增量算法，該算法全不同于局部換邊法的新算法，后來這兩種算法都被稱為 Delaunay 三經(jīng)典增量算法。點集的 Delaunay 三角劃分算法是本文中將限定區(qū)域劃分為單元網(wǎng)格算，而上文中又提到 Delaunay 三角化增量算法主要包括局部變換算法和yer/Watson 算法兩種，兩面將詳細敘述兩算法的具體的操作內(nèi)容。局部變換算法是一種增量算法，其主要是將點集中的點逐個添加到已經(jīng)區(qū)域，將包含點的三角形的三條邊分別以直線相連形成三個三角形，然意一條所對的兩個三角形是否滿足 Delaunay 準則，如果不滿足則將該改為兩三角形所對應(yīng)的另一條對角邊，一直迭代到所有的三角形都滿足aunay 準則。具體的圖像表示如圖 2.2 所示，當(dāng)向一個滿足 Delaunay 準化的模型中添加黑點，然后將點和所對應(yīng)的三角形的節(jié)點連接，再以添中心依次向外進行換邊操作，直至所有三角形都滿足條件為止。

【參考文獻】

相關(guān)期刊論文 前10條

1 吳小平;劉洋;王威;;基于非結(jié)構(gòu)網(wǎng)格的電阻率三維帶地形反演[J];地球物理學(xué)報;2015年08期

2 宋滔;王緒本;;點源二維電場正演的一組新的波數(shù)[J];地球物理學(xué)進展;2014年01期

3 吳曲波;張志勇;柯丹;范利飛;;點源二維直流電阻率法正演模擬[J];工程地球物理學(xué)報;2014年01期

4 潘克家;湯井田;;2.5維直流電法正演中Fourier逆變換離散波數(shù)的最優(yōu)化選取[J];中南大學(xué)學(xué)報(自然科學(xué)版);2013年07期

5 孟憲海;李吉剛;楊欽;蔡強;陳其明;;復(fù)雜限定Delaunay三角化算法[J];中國科學(xué):信息科學(xué);2010年03期

6 湯井田;公勁U

本文編號：2718884