導語：咸淡水共存區物探資料處置方法一文來源于網友上傳，不代表本站觀點，若需要原創文章可咨詢客服老師，歡迎參考。

應用地球物理方法在廣泛的咸水分布區確定淡水資源的存在特征，是西北地區地下水資源勘查工作的一項重要內容。目前，地層電阻率值是評價地下水礦化度唯一有效的地球物理參數¨J，因而，在咸淡水共存區開展地球物理勘查工作，高分辨的音頻大地電磁法(AMT)是首選方法。在獲取高質量AMT原始數據的基礎上，如何提高縱向分辨率、去除噪聲、靜態等各種干擾，得到地層真實的電阻率值，為地下水礦化度評價工作做好鋪墊，是地球物理資料處理所要解決的問題。國內許多學者對孔隙地下水礦化度評價進行了應用研究，但從筆者近幾年來在西北地區工作的經驗來看，用地層電阻率值，尤其是通過地面物探工作獲得的地層電阻率值來評價地下水礦化度，并非簡單的套用公式就可解決，不同地區、不同水文地質環境的差異會對計算結果產生很大的影響。實際工作中，礦化度評價數學模型建立是值得探討和研究的問題。

1AMT資料處理技術

1．1增加頻點數，提高AMT縱向分辨率

雖然常規的AMT法頻點密度在對數等問隔的一個級次內的取樣數設定為7～l0個被認為是合理的，但一些學者的研究表明，對于同一地電模型，AMT響應的原始曲線頻點數多少雖未有明顯差異，但其反演結果卻存在明顯差異j，增加頻率點數確實能提高AMT反演的精度，進而提高縱向分辨率。目前，V8綜合電法儀和EH-4電導率成像系統在我國地下水勘查工作中已被廣泛使用。V8系統后續處理軟件SSMT2000和MTeditor在將時序文件轉換為阻抗文件處理的過程中能提取大量有效的頻點數據。化希瑞等人的研究表明，可通過后續資料處理技術來恢復EH一4系統軟件從功率譜文件向阻抗文件轉換過程中丟失的頻點數據，大大增加數據數量。

1．2去噪處理

資料的去噪處理方法主要針對場源噪聲和地質噪聲特點而選取，人文噪聲和觀測裝置引入的噪聲主要靠野外施工加以避免。對于由人文噪聲引起的局部視電阻率畸變，可用阻抗相位恢復的方法進行校正。

1．2．1時序文件篩選

野外測量工作完成后，即進入原始資料的預處理階段。在該階段首先是對照野外原始記錄檢查各個測點的坐標及電極距的正確性；其次，重新回放每個測點的原始時序文件，并逐點、逐屏地對時序進行挑選，剔除那些存在明顯干擾信號的時間序列段，以減少隨機干擾信號對數據的影響；最后，利用采集軟件重新處理篩選過的時序資料。這樣，最大程度地保證了估計出的大地電磁響應的質量。

1．2．2人工圓滑去噪

人工圓滑去噪實際上是人機聯作手工剔除飛點，并根據視電阻率、相位曲線的變化趨勢對資料進行人工編輯。人工圓滑對處理天然場“死頻帶”區數據以及電磁干擾下的畸變數據效果很好，但需資料處理人員對實測資料和測區地質資料有著較多的了解；同時由于人工圓滑有著很大的隨意性，編輯后的結果會因人而異，從而顯得不夠科學嚴密。

1．2．3相位資料校正畸變的視電阻率

理論研究和實測資料表明，AMT資料的相位數據的質量要優于視電阻率資料的質量。在二維地質背景下，電性主軸上阻抗相位不受電磁噪聲影響。因而，我們可以利用相位資料恢復畸變視電阻率J。利用阻抗相位計算視電阻率的遞推公式為：

1．2．4低通濾波去噪

在視電阻率斷面上，高頻信號表現為局部的振蕩和異常，致使整個斷面零亂、無規律可循。為突出全局(區域)信息，揭示全局變化趨勢，可采用低通濾波法，截斷高頻成分，保留低頻成分。圖2為某地EH-4資料高斯低通濾波前后的對比，可以看出，低通濾波后的斷面圖視電阻率等值線更平滑，更好地揭示了從山前至盆地中心地層電性的變化。由相位推算的視電阻率可能會丟失一些細節，同時由相位計算視電阻率采用的是遞推算法，會產誤差累積效應。實際應用中，建議選擇實測視電阻率畸變段兩側的正常頻點的視電阻率數據作為計算的初始值，由相位計算出畸變段的視電阻率。圖1為原始曲線與相位計算視電阻率曲線對比，可以看出，通過相位計算出的天然場“死頻帶”區的視電阻率數據與前后數據很好地銜接，Bostick反演結果也得到明顯改善。

1．3靜態校正

理論研究表明，靜態干擾如使電阻率參數改變a倍的話，那么反演后的厚度參數將改變倍。因而，靜態校正是AMT資料處理中的關鍵問題，如果校正不當，會使后續的反演解釋得出錯誤的結果，靜態校正的方法很多，但最實用的還是EMAP濾波法，EMAP處理中最重要的一點是圓滑系數的選擇，如系數小，則不能很好地壓制靜態效應；如系數太大，則又使深部信息損失，降低深部的分辨率。實際處理中選擇不同的圓滑系數進行對比，輔以相位資料對處理結果進行判斷。

1．4反演

目前AMT數據反演主要為一維和二維反演。一維反演方法有Bostick、一維連續介質、一維0c—cam；二維反演方法有RRI、OCCAM、NLCG等。實際應用中，可采用不同反演方法和初始模型獲得一組數據的多個反演結果，綜合分析研究區物探成果，并結合已知的地質資料和鉆孔資料，對反演結果進行取舍。實踐證明，對于高頻連續剖面工作模式下的AMT數據，可采用EMAP濾波后的一維Bostick反演結果作為最終成果；而對于取得的非連續的長剖面數據，二維反演結果較一維反演結果更精確。

2礦化度評價的數學模型建立

2．1礦化度評價數據模型構建方法探討

利用地層電阻率值計算松散孔隙水的礦化度其原理和評價流程可參照文獻[1]和文獻[6]。此處僅對地下水礦化度與地層電阻率值問數學關系的構建予以初步探討。文獻[6]給出了礦化度評價公式C=[(1+△)•p／5．6P]，(1)At=to+tt‘h一18；其中，C為地下水礦公度；t為大地表層溫度，即研究區地殼常溫帶溫度，此值可在區域地質報告和地方氣象資料中獲得；t為地溫梯度；h為數據點的深度；OL為溫度常數(通常取值0．025Q•m／~C)；為礦化度常數，對不同成分和溫度的溶液，等于其礦化度為1g／L時的電阻率對數值，對NaC1型水質可取值一0．95；P為地層電阻率值；P為孑L隙參數。在研究區工作程度低，未有充足的測井資料和鉆孔資料可以借鑒的情況下，可借鑒式(1)，根據地層情況，選擇適合的孔隙參數值，進行初步礦化度評價工作；礦化度評價結果與已知井水質信息進行對比，對評價公式給予適當修正，進而形成適合于全區的地下水礦化度評價數據模型。當研究區有多個鉆孔資料可以利用時，可開展孔旁頻率測深工作，對測深結果進行反演解釋，依據解釋的電性層與不同含水巖組，對地下水礦化度進行統計，利用非線性回歸法建立地層電阻率與地下水礦化度相互關系模型。

2．2礦化度評價實例

2．2．1新疆疏勒縣地下水礦化度評價

新疆疏勒縣英爾力克鄉位于蓋孜河沖洪積平原下游，第四系沉積物以細砂為主，薄層亞黏土和黏土呈不連續分布，地下水具有上咸下淡的普遍特點。鑒于工作區缺少有效的鉆孔資料，礦化度評價采用式(1)。依據工作區地層結構特征，孔隙參數P=(3一)／2，=0．42，工作區年平均氣溫為11℃，推算出地下水礦化度與地層電阻率問的數學關系為C=(0．825p／17．2)曬+0．2，(2)其中0．2為修正系數。圖3為其中一條短剖面的EH-4勘查結果，圖3a為經阻抗文件重構、去噪、靜態校正處理后的EMAP擬二維反演結果，圖3b是由式(2)計算得到的礦化度斷面。鉆孔位于剖面100m處，孔深300m，取水段位于240～300m深度，地下水礦化度為0．53g／L，地下水礦化度值與電法資料預測值相吻合。

2．2．2銀川靈沙鄉地下水礦化度評價

靈沙鄉位于銀川平原東北部，黃河西3km處，由于蒸發濃縮作用，該處潛水的礦化度普遍偏高(>5g／L)，而埋深15013''''1以深的第二承壓水由于地下水徑流滯緩，礦化度也高于寧夏回族自治區的飲用水標準，僅有埋深50～150m的第一承壓水水質相對較好。在靈沙鄉工作開展之前，在已知井旁展開了EH-4試驗工作，礦化度評價分別采用式(1)和前人利用數學統計法建立的數學關系式C：(38．37／p)一0．27，(3)結果表明，式(3)計算的礦化度與已知孔資料更接近，而利用式(1)建立起的數學模型與實際相差較大。為便于管理，井位擬定在靈沙鄉政府院內。因無法開展EH_4工作，于是沿地下水徑流方向，在距鄉政府東西各1km處進行地下水勘查的物探工作，并分別進行了礦化度評價，兩者算術平均作為擬定井位處的礦化度結果(圖4)。依據估算礦化度結果，建議鉆孔取水深度控制在50～100ITI范圍內，考慮到出水量，最后取水段選在50～1501TI范圍，礦化度為2．9g／L，符合寧夏回族自治區飲用水標準。

3結語

在西北干旱的咸淡水共存區開展地下水勘查的地球物理工作，資料處理非常重要，依據工區環境和地質、水文地質條件，選擇合理的處理方法和資料處理流程。不同物探手段和不同反演方法都會造成獲得地層電阻率值的差異性，在同一工區利用同一礦化度評價數學模型要注意工作方法和反演方法的一致性，如若不同，則需通過試驗對模型進行適當的修正。建立礦化度評價數學模型的途徑有多種，要根據實際情況進行甄別。