本站小編為你精心準(zhǔn)備了電法勘探反演建模的應(yīng)用參考范文，愿這些范文能點(diǎn)燃您思維的火花，激發(fā)您的寫作靈感。歡迎深入閱讀并收藏。

《桂林理工大學(xué)學(xué)報》2016年第四期

摘要:

電法勘探中各種反演解釋方法的效果與初始模型的設(shè)計有關(guān)。實(shí)測數(shù)據(jù)的方向梯度包含了與地層電性結(jié)構(gòu)有關(guān)的信息。利用電法測深數(shù)據(jù)中方向梯度和地層結(jié)構(gòu)的內(nèi)在聯(lián)系，可以設(shè)置盡可能詳細(xì)精確的反演初始模型。本文研究了電阻率測深數(shù)據(jù)的水平梯度在各種地層結(jié)構(gòu)條件下的變化規(guī)律，提出利用水平梯度信息輔助建立初始模型的思路。利用水平梯度獲得地下電性結(jié)構(gòu)橫向變化的位置，建立更為直觀和接近實(shí)際的初始模型，提高電法勘探反演解釋水平。理論研究和實(shí)例證明，水平梯度可以作為建立反演模型的輔助信息。

關(guān)鍵詞:

電法勘探;方向梯度;建模

電阻率測深法一直是資源和工程勘察的常用方法。在實(shí)際工作中，最需要解決的問題是判斷地層結(jié)構(gòu)并獲得層參數(shù)。由于反演算法的局限性和地質(zhì)條件的復(fù)雜性，在實(shí)際資料解釋工作中，會遇到以下一些問題:各種反演方法中的初始模型設(shè)置對解釋結(jié)果有很大的影響，設(shè)計過于簡單的模型，在擬合復(fù)雜模型的過程中容易向最先出現(xiàn)的區(qū)域最優(yōu)解靠近;而在缺乏先驗信息的條件下，設(shè)置初始模型常帶有隨意性或主觀性。電阻率測深或其他測深類方法的反演結(jié)果是地下介質(zhì)的電性成像，雖然擬合誤差非常理想，但地層結(jié)構(gòu)、層參數(shù)等信息都不明確，需要數(shù)據(jù)處理人員進(jìn)一步分析解釋。因此，建立一個有地層結(jié)構(gòu)、層參數(shù)等信息的初始模型將提高反演解釋的準(zhǔn)確度。為了達(dá)到這一目的，程勃等［1－3］提出基于遺傳算法和統(tǒng)計學(xué)的電阻率測深二維反演算法:識別地層結(jié)構(gòu)、設(shè)計反演初始模型部分采用統(tǒng)計學(xué)方法，修改參數(shù)部分采用遺傳算法。這種方法的優(yōu)點(diǎn)是對經(jīng)典方法的缺陷進(jìn)行了補(bǔ)充，先用統(tǒng)計判斷對地電結(jié)構(gòu)進(jìn)行定性，在定性的基礎(chǔ)上嘗試定量解釋。在統(tǒng)計學(xué)建模方法中，為確定地下介質(zhì)橫向變化的空間位置，程勃等采用了方向梯度信息。電阻率法研究點(diǎn)電源場在不同的二維或三維地層結(jié)構(gòu)條件下的分布規(guī)律。對電阻率測深來說，是研究電阻率隨深度(或水平方向)的分布規(guī)律。一條電阻率測深剖面有多個測點(diǎn)和多個極距的測量數(shù)據(jù)，同一測點(diǎn)不同極距的視電阻率測深曲線反映了地下電性結(jié)構(gòu)的分布特征，測深曲線的垂直梯度反映了相鄰地層間電阻率的數(shù)量差異。一條剖面上同一極距不同測點(diǎn)的視電阻率直觀表達(dá)了地下電性結(jié)構(gòu)的橫向變化，水平梯度的變化反映了橫向變化地層間的電性差異水平。20世紀(jì)60—80年代，電阻率測深法的反演解釋主要采用量板和對數(shù)坐標(biāo)紙為數(shù)據(jù)處理工具。

在使用量板解釋的過程中，為了利用電阻率測深曲線的細(xì)節(jié)變化，解決復(fù)雜的地質(zhì)問題，研究者利用電測深曲線的導(dǎo)數(shù)(相對極距求梯度)推測界面深度的方法，稱為反射系數(shù)(K系數(shù))方法，并總結(jié)了一些經(jīng)驗。由于方法的基礎(chǔ)是一維層狀大地，在地質(zhì)條件適合的情況下，可以獲得較好的反演結(jié)果［4－7］。以前人們利用視電阻率梯度概念主要研究一維問題，在地層結(jié)構(gòu)復(fù)雜條件下，無法解決二、三維地質(zhì)問題，反演結(jié)果的誤差較大。但隨著儀器設(shè)備和計算機(jī)技術(shù)的發(fā)展，用高密度電法獲得大量的二維測深數(shù)據(jù)，可以更精細(xì)地反映地下的地電結(jié)構(gòu)。而電法勘探數(shù)值模擬技術(shù)的發(fā)展，使各種測深類方法的正演計算易于實(shí)現(xiàn)。實(shí)際上，除了垂直梯度，電阻率測深數(shù)據(jù)中的水平梯度變化同樣反映了地電場的變化。視電阻率水平梯度概念仍然可以作為一種了解電性在不同方向變化的技術(shù)手段，用以研究視電阻率斷面沿深度以及各個水平方向上的差異，為提取層界面位置、異常體位置等信息起參考作用。這些信息在建立較真實(shí)的地層結(jié)構(gòu)模型、給出直觀的反演結(jié)果中起到主要作用［8－10］。

1電阻率測深數(shù)據(jù)的垂向梯度特征及應(yīng)用

電阻率測深數(shù)據(jù)的垂向梯度特征為人們所熟悉。例如，隨深度(極距)變化的視電阻率測深曲線可能有2層(D、G)、3層(H、K、A、Q)以及4層等多種類型，曲線的類型表示了勘探深度范圍由淺到深大地的電性變化特征及層數(shù)。電阻率測深曲線的斜率(梯度)變化能夠反映出各電性層之間的電阻率差異關(guān)系，梯度的極值點(diǎn)或符號變化點(diǎn)的位置可以反映電性界面的深度等。在水平層條件下根據(jù)電測深曲線垂直梯度的分布規(guī)律可以確定各層電阻率和層厚度。電測深曲線垂直梯度的特征點(diǎn)可以作為分層的標(biāo)志，這些特征點(diǎn)所在的極距與界面的深度也有近似的倍數(shù)關(guān)系。雖然不是非常準(zhǔn)確，但作為正演模擬的初值已經(jīng)非常接近真實(shí)值了。在地層結(jié)構(gòu)為二、三維情況下，電測深曲線受旁側(cè)地質(zhì)體影響，導(dǎo)致曲線形態(tài)和垂直梯度的數(shù)值受到影響。因此垂直梯度僅能在界面傾角較小的情況下，判斷層參數(shù)。

2電阻率測深數(shù)據(jù)的水平梯度特征及應(yīng)用

一個剖面上的電阻率測深數(shù)據(jù)，可以看成是不同極距的電剖面法數(shù)據(jù)。根據(jù)在界面附近對稱四極電阻率剖面的理論曲線可知，在界面處視電阻率會有突變，在突變點(diǎn)處水平方向梯度出現(xiàn)極值點(diǎn)。因此，利用水平梯度數(shù)據(jù)可以了解地層結(jié)構(gòu)的橫向變化的位置、兩側(cè)電性差異水平，利用水平梯度的曲線形態(tài)，可以估計界面的傾斜方向。研究不同地層結(jié)構(gòu)的水平梯度時，采用二維有限元方法對存在不同橫向變化的地層結(jié)構(gòu)模型進(jìn)行正演計算，獲得電阻率測深數(shù)據(jù)。水平梯度用式(2)計算，Δx為10m。水平梯度曲線圖中縱坐標(biāo)為水平梯度Δd，橫坐標(biāo)為測線長度。由于各極距(AB/2)的曲線型態(tài)類似，僅繪出3條有代表性的水平梯度曲線，分別為AO=20、50、100m。

2.1一個垂直接觸面模型1:覆蓋層h1=8.5m，ρ1=100Ωm，接觸面兩側(cè)ρ2=1000Ωm，ρ3=300Ωm，接觸面位于測線150m處。水平梯度曲線如圖1所示。由于ρ2大于ρ3，所以在界面附近水平梯度為負(fù)。過界面時出現(xiàn)極小值，水平梯度的極小值都集中在同一個位置上(測線160m)，水平梯度的極距越小，幅度也越小，這主要是由于覆蓋層的影響。可以看出，極距的變化影響水平梯度的幅度，但不影響極值的位置。

2.2兩個垂直接觸面模型2:覆蓋層h1=3.5m，ρ1=100Ωm，覆蓋層以下ρ21=300Ωm，ρ22=1000Ωm，ρ23=300Ωm，水平分界面位于測線100和210m處。水平梯度如圖2所示。由于ρ22大于ρ21和ρ23，所以在100m界面附近水平梯度為正，在210m界面附近水平梯度為負(fù)。過界面時出現(xiàn)極大和極小值。可以看出，極距的變化影響水平梯度的幅度，但不影響極值的位置。從圖2的計算結(jié)果可以推斷，當(dāng)?shù)貙又写嬖谟幸欢▽挾鹊母咦杌虻妥璋鍫铙w時，水平梯度會有明顯的兩個極值顯示，極值的位置可能是界面的位置。而板狀體的兩側(cè)電阻率不同時，極值的幅度也不同。根據(jù)式(2)計算出的水平梯度與界面兩側(cè)的電阻率無關(guān)，僅與兩側(cè)電阻率的比值有關(guān)，因此根據(jù)水平梯度也可以估計板狀體的電阻率。從圖1和圖2可以看出，如果有垂直接觸面存在，電測深曲線的水平梯度會出現(xiàn)明顯的異常，根據(jù)異常的位置可以確定界面的位置。這對于建立復(fù)雜地電條件下的初始模型有重要的作用。

2．3傾斜接觸面模型3:覆蓋層h1=5m，ρ1=100Ωm;接觸面兩側(cè)ρ2=1000Ωm，ρ3=300Ωm，傾斜分界面位置的頂點(diǎn)在測線175m處，界面傾角135°;模型4與模型3參數(shù)相同，頂點(diǎn)位置在160m處，界面傾角155°，水平梯度如圖3所示。由于ρ3小于ρ2，所以在界面附近水平梯度為負(fù)。在接觸面左側(cè)水平梯度幅度增大，過界面時出現(xiàn)極小值，與垂直接觸界面相同，水平梯度的極大值都位于界面頂點(diǎn)處，水平梯度在傾斜方向變化較緩，過界面后較快降至正常水平。在圖3中，水平梯度的極距越小，幅度也小，這主要是受覆蓋層的影響。可以看出極距的變化影響水平梯度的幅度，但不影響水平梯度的形態(tài)。與直立接觸面的差異在于水平梯度在傾斜方向幅度增大。因此，如果有傾斜接觸面存在，電測深曲線的水平梯度會出現(xiàn)有一定寬度的明顯的異常，與垂直界面的一個點(diǎn)的窄異常有明顯的差異，但根據(jù)極值點(diǎn)仍可以確定界面頂點(diǎn)的位置，并用水平梯度的異常寬度范圍推測界面的傾角。

2.4傾斜斷裂帶模型5:覆蓋層h1=5m，ρ1=100Ωm;斷裂帶兩側(cè)圍巖ρ2=1000Ωm;斷裂帶ρ3=100Ωm，界面傾角65°，斷裂帶位于130～160m范圍內(nèi)。水平梯度曲線見圖4，水平梯度在140和170m處出現(xiàn)極小和極大值，與界面頂點(diǎn)位置一致。水平梯度的極距越小，幅度也越小，這主要是覆蓋層的影響。從圖4中可以看出極距的變化影響水平梯度的幅度，但不影響水平梯度的形態(tài)。與直立地質(zhì)體(圖1)相比，水平梯度曲線在傾斜方向幅度增大。根據(jù)直立或傾斜接觸面(或斷裂)上的水平梯度特征可知，如果電性突變，且為直立界面，各極距的異常會集中在同一個點(diǎn);而在傾斜界面，各極距的水平梯度異常有一定的寬度，傾角越大，異常越寬，但異常的極值點(diǎn)仍然在同一個點(diǎn)上。這一特征可以作為判斷局部地質(zhì)體、接觸面或斷裂的重要標(biāo)志。

3水平方向梯度在建立初始模型中的應(yīng)用方法

在建立初始模型時，如果地層結(jié)構(gòu)被判斷為含橫向變化(如存在斷裂)，則需要確認(rèn)電性橫向變化的位置及兩側(cè)的地層結(jié)構(gòu)和層參數(shù)。為此，統(tǒng)計電阻率測深數(shù)據(jù)各極距視電阻率的水平方向梯度。水平梯度的極值位置信息可以輔助完成初始模型結(jié)構(gòu)類別判斷，并獲得斷裂或巖性接觸面的初始位置。數(shù)值模擬結(jié)果和實(shí)測資料證明，淺部的局部異常體，同樣會造成水平梯度出現(xiàn)極值，可能造成對主要異常體位置判斷失誤。因此，在水平梯度統(tǒng)計過程中根據(jù)探測目的設(shè)置統(tǒng)計權(quán)重，如在探測斷裂、巖性接觸帶或深部異常體時，對深部信息賦予較大的統(tǒng)計權(quán)重。以低阻斷裂為例，斷裂左邊界確定函數(shù)為iz=P(min［X］);(3)斷裂右邊界確定函數(shù)為iy=P(max［X］)，(4)Xi，j=Δρi，jΔiAB()2j。(5)其中:Δρi，jΔi為視電阻率在水平方向上的變化率;AB()2j為權(quán)重因子;i為測點(diǎn)號;j為極距號;P(x)的作用是求x在測線上的位置。如果是巖性接觸帶類模型，則視接觸面兩側(cè)的電阻率大小關(guān)系采用式(3)或(4)中的一個。

4應(yīng)用實(shí)例

廣西天然氣輸氣管道工程的柳州段，有斷裂或巖性接觸帶存在，是需要重點(diǎn)勘探的區(qū)域。測區(qū)大部分為第四系沖洪積(Qal+pl4)粘土、粉質(zhì)粘土，測線小號點(diǎn)約40～50m范圍內(nèi)有砂巖出露。測線大號點(diǎn)約180m處，鉆探取芯的情況為:第四系沖洪積(Qal+pl4)粘土、粉質(zhì)粘土層;含圓礫粘土、卵石層;下部為石炭系中統(tǒng)(C2)白云巖、灰?guī)r。要求利用高密度電法勘探查明斷裂(砂巖和白云巖的接觸位置)及斷裂兩側(cè)地層分布等情況。沿設(shè)計管道軸線在地表的投影以及左右兩側(cè)各10m布置了3條高密度電法測線。測量儀器為WDJD－2型多功能數(shù)字高密度電阻率儀。根據(jù)勘查要求，采用聯(lián)合三極裝置，即在一條測線上同時進(jìn)行A－MN和MN－B測量。采用的測量電極距為10m，點(diǎn)距10m，測量層數(shù)為8層，對應(yīng)供電極距AO為15～85m。處理數(shù)據(jù)之前，將聯(lián)合剖面數(shù)據(jù)轉(zhuǎn)換為電阻率測深數(shù)據(jù)。以A2線為例，圖5中右側(cè)的曲線是極距為75和85m的視電阻率曲線。可以看到，以測線150m為界，兩側(cè)有明顯的電性差異，但視電阻率有數(shù)個極值點(diǎn)。測線的水平梯度曲線如圖5所示，圖中僅繪出了極距(AB/2)為75和85m的水平梯度曲線，水平梯度的極大值點(diǎn)僅出現(xiàn)在測線的150m位置。斷層(巖性接觸面)位置用斷裂邊界確定函數(shù)式(4)確定接觸面的位置。由此，初步確定斷層在測線上的位置為150m。

5結(jié)束語

從理論分析和實(shí)測數(shù)據(jù)處理結(jié)果可以看出，利用測深數(shù)據(jù)的水平梯度方法可以提供地下電性橫向變化的信息。因此，可以利用水平梯度作為輔助信息建立接近真實(shí)地層結(jié)構(gòu)的初始模型，提高了反演的準(zhǔn)確度。激發(fā)極化測深和其他電磁法測深類方法和電阻率測深有類似之處，因此也可以應(yīng)用水平梯度作為建立反演模型的輔助方法。

參考文獻(xiàn):

［1］程勃，底青云．基于遺傳算法和統(tǒng)計學(xué)的電阻率測深反演研究［J］．地球物理學(xué)進(jìn)展，2012，27(2):788－795．

［2］程勃，底青云．基于偏導(dǎo)數(shù)和統(tǒng)計學(xué)方法的電阻率測深二維反演［J］．石油地球物理勘探，2012，47(6):1006－1013．

［3］程勃，底青云．復(fù)雜地電結(jié)構(gòu)條件下統(tǒng)計學(xué)建模法電阻率測深二維反演［J］．地球物理學(xué)報，2014，57(3):961－967．

［4］阮百堯，村上裕，徐世浙．電阻率/激發(fā)極化率數(shù)據(jù)的二維反演程序［J］．物探化探計算技術(shù)，1999，21(2):116－125．

［5］蘇朱劉，胡文寶，嚴(yán)良俊．電阻率和極化率測深法的正演修正法反演［J］．石油物探，2005，44(2):194－198．

作者：程勃 單位：桂林理工大學(xué) 廣西隱伏金屬礦產(chǎn)勘點(diǎn)實(shí)驗室