徐凱軍, 李 猛, 季春暉, 梁舒瑗

(中國石油大學(華東)地球科學與技術(shù)學院,山東青島 266580)

目前深部礦產(chǎn)資源勘探已成為中國勘探重點[1]。井中磁測對深部礦體的磁場變化非常敏感,在深部金屬礦勘探中發(fā)揮著重要作用[2-4]。考慮到井的數(shù)量有限,因此結(jié)合地面磁測數(shù)據(jù),進行井-地磁異常聯(lián)合反演,可以有效提高反演結(jié)果的空間分辨率[5-7]。三維磁異常反演計算量大,為了提高地球物理數(shù)據(jù)反演計算效率,一些學者開展了數(shù)據(jù)空間地球物理反演算法,極大地減少反演時間[8-11]。加入先驗信息約束可以提高反演可靠性[12],基于模糊C均值聚類(FCM)的約束反演,可以將先驗地球物理信息加入到反演中,對反演結(jié)果進行聚類,使反演的礦體邊界清晰,反演結(jié)果更符合地質(zhì)解釋[13-19]。井中磁測不僅可以觀測不同深度的三分量磁異常,同時鉆井巖心、磁化率測井曲線等信息能夠提供較精確的地下物性信息,更適合進行FCM約束反演,提高反演可靠性。筆者將數(shù)據(jù)空間反演方法和FCM約束反演方法相結(jié)合,實現(xiàn)基于模糊C均值聚類約束的井地磁異常數(shù)據(jù)空間聯(lián)合反演。

1 井地磁異常正反演理論

1.1 井中磁異常正演方法

磁法正演通常將地下剖分為一系列長方體網(wǎng)格，采用不規(guī)則網(wǎng)格，礦體區(qū)域網(wǎng)格剖分較為細密。井中觀測點接收到的磁場等于所有地下網(wǎng)格在觀測點所引起的磁場總合。觀測點處產(chǎn)生的磁場三分量異常表達式如下：

(1)

其中

式中，Bx、By和Bz分別為x、y、z三個方向的磁異常；μ0為真空磁導率；Mx、My和Mz分別為x、y、z三個方向的磁化強度。

1.2 井-地磁異常聯(lián)合反演方法

1.2.1 數(shù)據(jù)空間反演算法

三維磁法模型空間反演通常計算量大,數(shù)據(jù)空間反演是基于模型空間反演算法的一種改進,減少反演時間,提高了反演效率。

磁法數(shù)據(jù)反演的目標函數(shù)通常為

(2)

式中,φ(m)為目標函數(shù);φd(m)為數(shù)據(jù)擬合目標函數(shù);φm(m)為模型目標函數(shù);d為觀測數(shù)據(jù)向量;F(m)為正演數(shù)據(jù)向量;Wd為數(shù)據(jù)對角加權(quán)矩陣;m為模型向量;Wm為模型加權(quán)矩陣;m0為參考模型向量。反演計算就是對目標函數(shù)求解極小值的過程。因此對目標函數(shù)求偏導,并令其等于0,可得

(3)

式中,J為雅克比矩陣。基于式(3)求解,可以得到反演迭代公式,常規(guī)磁法數(shù)據(jù)反演都是基于模型空間,其反演迭代公式為

m=m0+(JTWdJ+Wm)-1JTWd(d-F(m)+J(m-m0)).

(4)

式中,JTWdJ+Wm為M×M階矩陣;M為模型空間個數(shù)。

在模型空間反演迭代公式基礎(chǔ)上進行轉(zhuǎn)換,得到數(shù)據(jù)空間反演迭代公式[8]:

(5)

1.2.2 模糊C均值聚類約束反演

模糊C均值聚類算法通過指定聚類中心的個數(shù),自動計算出每個樣本數(shù)據(jù)相對于聚類中心的“親昵程度”(隸屬度),從而確定樣本點的類別,完成數(shù)據(jù)分類的目的。FCM的數(shù)學表達式[13]為

(6)

其中

式中,M為需要聚類數(shù)據(jù)的個數(shù);C為聚類中心的個數(shù);pj為第j個數(shù)據(jù);vk為第k個聚類中心值;μjk為第j個數(shù)據(jù)pj相對第k個聚類中心vk的隸屬度;q為隸屬度μjk的模糊系數(shù),取值范圍為[1,∞)。當q=1時,聚類結(jié)果會產(chǎn)生類似“硬聚類”的效果,即每個數(shù)據(jù)對象僅分配給一個聚類中心,q值越大時聚類效果越模糊。

Sun等[13]在FCM算法的基礎(chǔ)上加入了聚類中心參考項,改進后的FCM算法為

(7)

其中

式中,tk為第k個聚類中心的參考值;ηk為評價第k個聚類中心是否有已知參考值的系數(shù)。ηk的取值范圍為0或1,當?shù)趉個聚類中心有已知參考值時,ηk=1;當?shù)趉個聚類中心沒有已知參考值時,ηk=0。

mc=vμ.

(8)

加入FCM約束的反演方法在擬合觀測數(shù)據(jù)的同時對模型的物性做出了約束限制,因此反演結(jié)果即符合觀測到的地球物理場又滿足先驗的統(tǒng)計巖石物理信息,同時FCM聚類方法使反演異常體邊界更加清晰。

1.2.3 井-地磁異常聯(lián)合反演

常規(guī)地面磁測數(shù)據(jù)具有較高的水平分辨率,但縱向分辨率較差,而井中磁測數(shù)據(jù)具有較高的縱向分辨率,但受井數(shù)量限制,控制范圍有限。充分的利用地面磁測和井中磁測的優(yōu)點,進行井-地磁異常聯(lián)合反演,能夠提高磁異常反演結(jié)果的空間分辨率。

對于井-地聯(lián)合反演,式(2)中φd可寫為

(9)

式中,γ為權(quán)重因子;γ1后的部分為地面磁測數(shù)據(jù)的殘差;γ2→γn后的部分為不同井中磁測數(shù)據(jù)的殘差。

聯(lián)合反演中要使井中磁測和地面磁測的最大振幅保持在同一個數(shù)量級內(nèi),當井中磁測的振幅遠大于地面磁測振幅時,適當降低井中磁測數(shù)據(jù)的權(quán)重,使反演結(jié)果不會因為過于擬合井中磁測數(shù)據(jù)而造成地面磁測數(shù)據(jù)的欠擬合的問題。

2 井-地磁異常特征分析

2.1 井中三分量磁異常特征

地下為單一礦體,長、寬均為200 m,高為100 m,礦體頂面埋深為400 m。布設(shè)4個井,其中井L1穿過礦體,井L2～L4未穿過礦體,模型與井位分布如圖1(a)所示。礦體的磁化強度為1 A/m,磁傾角為90°,磁偏角為0°,地磁場傾角為90°,地磁偏角為0°。井中測點從z=0起,相隔20 m設(shè)一個測點,正演計算井中測點的三分量磁異常。圖1(b)為不同井的三分量磁異常曲線,圖中灰色部分顯示了礦體的深度范圍,井位于礦體中心線上(x=500 m),由于對稱性,磁異常Bx分量為0。可以看出,當井穿過礦體時,磁異常Bz分量在礦體位置出現(xiàn)顯著變化,變化極值點可以確定磁性體的上下邊界。當井沒有穿過礦體時,磁異常垂直分量Bz的峰值基本對應礦體中心位置,磁異常水平分量By正負峰值基本對應礦體的上下邊界,隨著井離礦體越遠,磁異常幅值越小。

圖1 模型及井中三分量磁異常曲線

2.2 井-地磁異常對比分析

設(shè)置地下同時存在3個礦體的模型,模型和井分布如圖2(a)所示。礦體的磁化強度為1 A/m,磁傾角為90°,磁偏角為0°,地磁場傾角為90°,地磁偏角為0°。圖2(b)為井中三分量磁異常曲線圖?？梢钥闯?磁異常具有很好地縱向分辨率,能夠識別出離井較近的兩個縱向疊置礦體不同位置,而礦體3由于離井較遠,異常較弱,井中三分量磁異常沒有明顯響應。圖2(c)為該模型的地面磁異常等值線圖,可以看到兩個明顯的異常區(qū)域,但無法識別縱向疊置的異常體。因此井中磁測具有較高的縱向分辨率,但受限于測井位置,橫向控制范圍有限,而地面磁異常具有較好的橫向分辨率,但縱向分辨率差,因此結(jié)合兩者的優(yōu)勢開展井-地磁異常聯(lián)合反演是提高反演分辨率的有效方式。

圖2 模型及井-地磁異常

3 模型反演

3.1 數(shù)據(jù)空間反演與模型空間反演對比

建立4個頂面埋深在300 m,尺寸為500 m×500 m×480 m的磁性體,如圖3(a)所示,磁化強度為2 A/m,磁化傾角為90°,磁偏角為0°,背景區(qū)域無磁性,地磁傾角為90°,地磁偏角為0°。異常體區(qū)域剖分成50×50×40個網(wǎng)格單元,每個網(wǎng)格的尺寸為50 m×50 m×40 m。地面上從坐標原點起,在x和y方向上隔50 m設(shè)置一個測點,觀測點總數(shù)為2 500個。圖3(b)為正演計算得到的地面磁異常等值線圖?；谄胀ㄎC分別開展數(shù)據(jù)空間反演和模型空間反演,分析反演的計算效率。基于模型空間反演耗時為50 min,而基于數(shù)據(jù)空間反演耗時為15 min,可以看到基于數(shù)據(jù)空間反演明顯減少了反演所用的時間。

圖3 模型及地面磁異常

圖4分別展示了模型空間和數(shù)據(jù)空間三維反演結(jié)果示意圖,可以看出模型空間反演和數(shù)據(jù)空間反演都能反演出礦體的空間位置,但反演的物性值與真實值2 A/m都有較大偏差。

圖4 反演結(jié)果三維示意圖

3.2 地面磁測反演與井-地磁測聯(lián)合反演對比

設(shè)計兩個傾斜脈狀體共同組合成的模型,模型中設(shè)置兩口井,井位坐標分別為L1(650,1 000)、L2(1 500,1 000),如圖5(a)所示。礦體在深度剖面上呈現(xiàn)“y”形狀,背景區(qū)域的磁性為零,礦體的磁化強度為2 A/m,磁化傾角75°,磁偏角25°。圖5(b)為正演得到的地面磁異常等值線圖,地面磁異常場的總體分布特征為中間高四周低,磁異常存在著兩個高值區(qū)間?；谡莴@得井中三分量磁異常,從z=0起間隔50 m設(shè)置一個測點,圖5(c)顯示了井中磁測三分量曲線。井L1和L2的磁異常Bz分量表現(xiàn)了明顯的穿礦特征。

圖5 模型及井-地磁異常圖

針對上述模型僅使用地面磁異常數(shù)據(jù)進行反演,反演結(jié)果如圖6(a)所示,圖中僅展示磁化強度高于0.4 A/m的反演結(jié)果,以展示反演結(jié)果的整體特征。反演結(jié)果顯示深部異常體形態(tài)與實際礦體有較大差異,說明地面磁測縱向分辨率差。反演結(jié)果的磁化強度極大值為1.4 A/m,與實際礦體的磁化強度值2 A/m也有較大的差異。

圖6 模型反演結(jié)果

聯(lián)合地面磁測和井中三分量磁異常進行聯(lián)合反演,反演結(jié)果如圖6(b)所示?？梢钥闯?相比于單獨地面磁異常反演,井-地磁異常聯(lián)合反演結(jié)果更接近于真實礦體模型,“y”型礦體得到很好地恢復,說明井-地磁異常聯(lián)合反演能夠結(jié)合地面磁測橫向分辨率高和井中磁測縱向分辨率高的優(yōu)點,提高反演結(jié)果的空間分辨率,但反演結(jié)果的磁化強度極值為2.4 A/m,與實際礦體的物性值仍存在一定的差異。

3.3 FCM約束反演

針對圖5(a)模型,設(shè)置不同的聚類參數(shù),分析不同聚類參數(shù)時的反演效果。真實模型的聚類中心為2,地質(zhì)背景磁化強度為0,礦體磁化強度為2 A/m。設(shè)置二組不同聚類參數(shù),第一組聚類中心數(shù)C為2,聚類中心參考值僅輸入一個0值(即反演中背景層的磁性已知,而礦體的磁性未知);第二組聚類中心數(shù)C為4,聚類中心參考值僅輸入0和2 A/m(反演中猜測地下存在4種磁性體,但僅已知兩種巖石的物性值)。

分別使用這二組參數(shù)開展基于FCM約束的井-地磁異常聯(lián)合反演。第一組的反演結(jié)果如圖7所示。可以看出,基于FCM約束反演的物性分布直方圖更接近于真實模型。反演異常體的邊界清晰,而且反演物性值與真實模型較為一致。

圖7 第一組參數(shù)的FCM約束反演結(jié)果

第二組參數(shù)的反演結(jié)果如圖8所示,由于聚類中心數(shù)C設(shè)為4,反演結(jié)果物性分布直方圖也顯示了磁化強度主要集中分布在0、0.4、1和2 A/m 4個區(qū)間內(nèi),反演異常體邊界變得模糊,與真實模型有了明顯差異。

圖8 第二組參數(shù)的FCM約束反演結(jié)果

通過不同聚類參數(shù)的反演結(jié)果對比分析可以看出,FCM約束反演的效果與聚類中心數(shù)C的選取有關(guān),C取值越接近真實情況,反演效果越好。在使用FCM約束反演方法時,聚類中心數(shù)C的取值應盡量與地下巖石的物性種類數(shù)目一致。鉆井和測井資料能夠提供較為精確的先驗物性信息,這也為FCM約束反演提供了好的基礎(chǔ)。

3.4 冬瓜山銅礦建模及反演

冬瓜山銅礦位于安徽銅陵地區(qū),已做過不同程度的礦產(chǎn)調(diào)查工作,該區(qū)的地質(zhì)構(gòu)造情況已有較為詳細的認知。以冬瓜山銅礦的地質(zhì)模型(圖9)為基礎(chǔ),調(diào)研了相關(guān)礦物的磁化率分布區(qū)間[20-21],然后建立冬瓜山地區(qū)的三維磁化強度模型,如圖10所示。礦體的磁化角度按當?shù)氐牡卮艌鼋嵌仍O(shè)置,即磁化傾角47.4°,磁化偏角-5.54°。

圖9 冬瓜山銅礦典型地質(zhì)模型[20]

圖10 冬瓜山銅礦三維磁化強度模型

在該模型中設(shè)置3口井,井位坐標分別為L1(700,500),L2(300,600),L3(1 000,1 000),井深1 000 m。地面測點間隔100 m,井中測點間隔40 m,通過正演計算分別得到地面磁力異常和井中三分量磁力異常,然后加入5%的高斯噪音作為觀測數(shù)據(jù)。圖11顯示了地面磁異常和井中三分量磁異常。

圖11 地面和井中磁異常

分別進行地面磁異常反演、井-地磁異常聯(lián)合反演和基于FCM約束的井-地磁異常聯(lián)合反演,分析不同反演的效果。圖12顯示了3種方法的反演效果。

圖12 不同反演方法結(jié)果對比

可以看出,僅使用地面磁異常反演效果較差,銅礦脈的埋藏較深且上方存在磁性巖體遮擋,地面磁異常反演無法從中分辨出不同的磁性礦體,更無法確定出埋藏最深的銅礦脈的位置形態(tài)。井地磁異常聯(lián)合反演提高了礦體的空間分辨率,上、中、下3部分的礦體在反演結(jié)果里都得以體現(xiàn),但礦脈的連續(xù)性較差,而且反演的磁化強度數(shù)值和真實值有一定的誤差。基于FCM約束的井地磁異常反演效果得到了進一步提升,反演的礦體邊界較為清晰,磁化強度值更接近于真實的礦體物性值。

4 結(jié) 論

(1)井中磁測具有較高的縱向分辨率。磁異常垂直分量的峰值基本對應礦體中心位置,磁異常水平分量正負峰值基本對應礦體的上下邊界。

(2)與模型空間反演相比,數(shù)據(jù)空間反演能夠減少反演時間,提高反演計算效率。

(3)井-地磁異常聯(lián)合反演,綜合了地面磁測橫向分辨率高和井中磁測縱向分辨率高的優(yōu)點,可以顯著提高反演結(jié)果的空間分辨率。

(4)在反演過程中引入FCM約束使反演結(jié)果既符合觀測到的地球物理場,又滿足先驗的巖石物理信息,反演的物性值更接近真實值,而且聚類算法使反演異常體邊界清晰,反演結(jié)果更符合地質(zhì)解釋。