九江—瑞昌礦集區的3D結構及對區域找礦的啟示

鄧 震，呂慶田，嚴加永，趙金花，劉 彥

中國地質科學院礦產資源研究所，國土資源部成礦作用和資源評價重點實驗室，北京 100037

1 引 言

九江—瑞昌（下稱“九瑞”）礦集區是我國長江中下游成礦帶7個大型礦集區之一［1-3］，是贛西北及鄂東南地區重要的銅（金）多金屬基地.雖然礦集區勘探及開采工作已開展多年，但直至今天，區內平均勘查深度僅為500m左右，平均開采深度為400m左右［4］.礦集區內已探明礦產的儲量及開采能力已不能跟上日益增長的資源需求.因此尋求區內500～2000m的第二深度空間［5］的找礦突破已成為勢在必行的工作.

九瑞礦集區相繼開展過一系列的地質、物化探及礦產普查等工作，地區研究程度較高，資料豐富.然而受到當時技術所限，大部分工作僅限于地表范圍，對深部地層及巖體的三維空間形態未能給出相應的描述及解釋.隨著重磁3D物性反演成像［6-8］技術、重磁多尺度邊緣檢測［9-13］等技術的發展，使利用網格重磁數據重建地下位場場源參數的分布成為可能.這一突破將有利于研究礦集區內地殼深部結構和動力學過程，以及其對成礦系統形成、演化的制約；有利于研究礦集區深部控礦構造的性質和空間形態；有利于探測、追蹤礦集區主要含礦、容礦地層的空間分布，及與成礦有關巖體（侵入巖、次火山巖）的空間位置、規模和埋深，并大致刻畫礦集區的三維結構.

作者在前人研究的基礎上利用最新的方法技術對原有重磁資料進行了重新處理，建立了九瑞礦集區的重磁3D地質模型，重新厘定了九瑞地區線型構造的空間展布，在此基礎上對九瑞地區重、磁場特征及結構進行了研究.

2 區域地質背景

九瑞礦集區位于長江中下游成礦帶西南端轉折處的最狹窄部位.礦集區內基底地層為中元古界雙橋山群厚度巨大的深海—淺海相淺變質巖系.出露的蓋層為中上統奧陶系碳酸鹽巖、上志留統砂頁巖、上泥盆統五通組含礫石英砂巖、中石炭統黃龍組白云巖和灰巖、中下二疊統和中下三疊統碳酸鹽巖、下第三紀紫紅色砂礫巖及第四紀松散沉積物，蓋層總厚度近10km，其中基底地層為重要含礦層位和含鎢花崗巖的巖源層，震旦系—寒武系賦存含銅沉積建造，奧陶系—三疊系為該區重要的銅金賦礦層位.該礦集區經歷了揚子與華北陸塊碰撞、燕山期的陸內造山和燕山期后的陸內斷陷三個不同的發展階段.金屬成礦作用主要集中在燕山期的陸內造山演化階段.燕山期的陸內造山（板內變形）是區內蓋層構造變形的主要時期.礦集區內褶皺斷裂構造發育，六個軸向近于平行的背斜、向斜組成緊密線狀褶皺帶.褶皺相對緊閉，常見倒轉，軸向北東東，自北而南有通江嶺—鄧家山復式向斜、大橋—寶山背斜、黃橋復式向斜、丁家山—大沖背斜、賽城湖—烏石街復式向斜和城門山—長山背斜等.區內主要出露淺成—超淺成中酸性—酸性侵入體，多呈巖墻、巖株、巖枝狀產出.主要巖石類型為花崗閃長斑巖，其次為石英閃長玢巖、石英斑巖以及少量的閃長巖、石英閃長巖和其它脈巖，其中主要含礦巖體為石英閃長玢巖、花崗閃長斑巖、石英斑巖.

3 物性研究及重磁數據

物性是聯系地球物理和地質的紐帶，扎實的物性資料是可靠地質解釋的基礎.將物性資料的研究和該地區的地質巖性特征結合起來，可以建立合理的地質模型，并以此作為異常解釋的基礎.

根據所收集的九瑞地區物性資料1）地礦部第1綜合物探隊.江西省九江—瑞昌地區物探、化探綜合調查成果報告：1／5萬，1992.，2）江西省地質礦產局物化探大隊.江西省瑞昌市武山礦區外圍綜合物探勘查工作成果報告，1992.，和中國地質科學院廊坊物化探研究所2008年實測資料［14］，得到研究區部分地層、巖礦體樣本的密度及磁性資料，對九瑞礦集區的地層巖石密度和磁性進行了具體的分析.

3.1 地層及巖礦體密度特征

綜合所收集的資料，對巖石密度以地層厚度為權重進行加權平均，得到各巖層的密度變化規律（如圖1所示）.

如圖1所示，礦集區地層密度差達0.05×103kg／m3以上的密度界面共有8個.其中二疊系上統雖與上下界面之間形成密度差，但由于其厚度較小，無法識別.因此本區地層密度差達0.05×103kg／m3以上的密度界面實際共有6個.大致將研究區分為七個密度層（第四系和第三系、中三疊統—中石炭統、上泥盆統—上志留統、下志留統—中奧陶統、下奧陶統、上震旦統、下震旦統）.

（1）第四系和第三系為低密度層，密度變化范圍較大，（1.832～2.487）×103kg／m3，均值為2.32×103kg／m3.

（2）中三疊統—中石炭統為高密度層，除二疊系上統龍潭組密度較低（2.568×103kg／m3）外，整個高密度層的密度值均在2.7×103kg／m3左右.

（3）上泥盆統—上志留統、下志留統—中奧陶統為中低密度層，其密度值多出現于 （2.5～2.6）×103kg／m3之間，與上覆高密度層具有（0.1～0.2）×103kg／m3的密度差.

（4）下奧陶統為高密度層，該層密度平均值為2.765×103kg／m3，為礦集區最高密度地層，與上覆低密度層之間密度差為0.234×103kg／m3.

（5）震旦系整體為低密度層，尤其是上統，密度平均值僅為2.4×103kg／m3，與上層高密度地層的密度差達到-0.365×103kg／m3.

礦集區內巖漿巖密度變化范圍較大.其中輝綠玢巖的密度最高，平均值為2.91×103kg／m3；石英閃長玢巖次之，平均值為2.88×103kg／m3.此兩種巖漿巖密度明顯大于地層巖石密度.除石英閃長玢巖和輝綠玢巖外，巖漿巖均表現為低密度巖體，其密度值與上泥盆統—上志留統密度層的值略低，平均值為2.58×103kg／m3.

礦集區主要礦石以及含礦巖石中以含銅黃鐵礦（塊狀硫化物銅礦石）密度最高，平均值為4.17×103kg／m3；各類矽卡巖和蝕變大理巖的密度次之，平均值為（2.97～3.8）×103kg／m3.

對比各沉積巖和巖礦體的密度值我們發現：各巖漿巖與大部分沉積巖的密度值之間雖存在密度差，但由于巖體的形狀及侵入的層位較為復雜，單純的重力3D反演可能無法完美地刻畫巖體的形態，因此，在數據處理過程中需要結合精細的二維剖面反演技術；礦體及含礦巖體與各沉積巖間的密度差較大，在巖礦體規模及其埋深合適的情況下，無論其圍巖是何巖性，重力場都能有很好的反應.

3.2 地層及巖礦體磁性特征

依據對九瑞地區地層巖石磁性的分析，礦集區從奧陶系到第四系，沉積地層巖石磁化率都很小，一般在10-4CGSM以下，可視為無（弱）磁性地層.各類火成巖具有較強磁性，在本區也遵循從酸性→中性→基性磁性逐漸增強的趨勢，侵入巖體大多數有較強磁性.侵入巖中以輝石橄欖巖磁化率最高，常見值為55000×10-6SI，剩磁以閃長玢巖最高，最高達到4000×10-3A／m；火山巖中磁化率以安山巖最高，常見值29000×10-6SI，剩磁以玄武巖最高，變化范圍為2300×10-3～7300×10-3A／m.各類侵入巖體、火山巖、夕卡巖和磁性礦石（如磁黃鐵礦、磁鐵礦等）的磁化率較地層磁化率至少高出3～4個數量級.因此，認為本區的航磁異常主要為與巖漿活動相關的地質構造所引起.而在礦集區尺度上，礦體規模遠小于巖體規模，因而引起磁異常的主要因素為含磁性的巖體.

3.3 九瑞地區重、磁場特征

3.3.1 重力場特征

本文采用的布格重力數據范圍為東經115°22′—116°04′，北緯29°30′00″—29°50′，包括1∶5萬分幅楓林幅、九江幅、瑞昌幅、范家鋪幅、廬山幅和湖口縣幅共六幅經數字化后拼接而成（圖2）.布格重力數據的比例尺為1∶5萬，布格重力異常總體強度比較弱，梯度緩，形態比較規整.區內重力數據整體呈北高南低趨勢，且北面重力高、低值帶界限清晰，直觀，明確，南面重力高、低值帶界限較模糊.異常特征綜合反映了蓋層和深層構造特征和巖漿巖等因素.局部異常強度反映了九瑞地區褶皺帶從北往南埋深逐漸增大.梯度變化的場態特征顯示斷裂（如城門山—武山—東雷灣一帶，丁家山—大沖一帶）.

3.3.2 九瑞地區磁場特征

本文所用航磁數據為國土資源遙感中心提供的研究區1∶5萬航磁異常125m×125m網格數據.數據經緯度范圍與布格重力數據范圍一致.航磁測量飛機的平均飛行高度為500m.為方便研究，本文對九瑞礦集區原始航磁數據進行了化極處理，其結果如圖3所示.研究區內磁異常總體較平緩，北部出現兩個局部異常高呈V字形，中部地區異常值較低，城門山附近出現局部小異常高，南部出現磁場抬升.礦點多分布于異常變化梯度帶上（武山—丁家山—城門山—雷家灣一帶都是類似特點）.

圖1 九瑞地區巖層密度變化規律圖（陰影中數字表示平均密度）Fig.1 Variation of lithospheric density in JiuRui district（the numbers in shadow are average density）

4 重、磁反演及模擬方法

4.1 重力3D反演

布格重力異常是莫霍面起伏及地殼內各類密度不均勻體在大地水準面上的綜合響應，而本次研究深度為5km以淺范圍，因此，必須進行區域場與局部場的分離.在對位場數據進行位場分離處理時，沒有相對固定的方法，需要結合地質情況，針對具體問題選擇合適的位場分離方法［15］.本文嘗試了趨勢場法、匹配濾波法和延拓法的分離試驗，經過與實際地質情況的綜合對比，趨勢場法的二次趨勢場背景場表現為單一的北西西向梯級帶，而匹配濾波法和延拓法所求得的背景場大體上雖表現為北西西向梯級帶，但是都有不同程度的局部場未被分離出去（匹配濾波法分離效果更差）.對比其剩余異常，也是趨勢場法求取的剩余異常相對較好地保留了5km內的信息（圖4），因此，決定采用該剩余異常進行重力3D成像反演.

圖2 九瑞礦集區1∶5萬布格重力異常Fig.2 Bouguer gravity anomaly（scale 1／50000）map of JiuRui ore district

圖3 九瑞礦集區1∶5萬航磁化極異常Fig.3 Aeromagnetic anomaly（scale 1／50000）reduced to pole map of the JiuRui ore district

圖4 九瑞礦集區1∶5萬布格重力二階趨勢剩余異常Fig.4 Second－order gravity residual anomaly（scale 1／50000）map of JiuRui ore district

網格剖分：2.5km以淺深度范圍內網格單元長×寬×高＝500m×500m×250m，2.5～5km深度范圍內網格單元長×寬×高＝500m×500m×250m，網格單元總數為141×81×20個.

密度差邊界設置為0～3g／cm3，通過迭代計算完成反演.圖5為礦集區重力3D反演不同深度水平切片圖，不同深度的高密度體分布與地表剩余布格重力異常呈正相關，-1000m深度范圍內，正負密度差擾動細節豐富，反映了淺表構造的復雜性，隨著深度增加，密度異常體形態逐漸增大并有相互連通的趨勢.在-4000m以深范圍，北部和中部高密度體分別連接成一體，呈回字狀.此外，礦集區西南角上的高密度體顯示結果不真實，其原因是本區缺失西南和東北角的重力數據，但是為了網格數據的計算方便，這兩個地區的數據是利用插值算法計算得出的.

4.2 航磁3D反演

網格剖分：2.5km以淺深度范圍內網格單元長×寬×高＝500m×500m×125m，2.5～7.5km深度范圍內網格單元長×寬×高＝500m×500m×125m，網格單元總數為147×80×40個.

圖5 礦集區重力3D反演不同深度水平切片圖Fig.5 3Dhorizontal slice of JiuRui ore district at different depth

磁化率邊界設置為0～1SI，通過迭代計算完成反演，反演結果如圖7所示，圖中最上層是航磁化極異常，中間是地形地質圖，下面是航磁3D反演結果立體圖，從藍色到紅色的顏色漸變表示磁化率從低到高的變化.

根據物性測量及分析結果，判斷航磁3D反演結果中的高磁化率部分主要反映了巖體分布情況.本區侵入巖磁化率變化范圍為（410～14800）×10-6CGSM，均值在2500×10-6CGSM左右，通過對比試驗，發現該磁化率閾值基本反演了侵入巖體的邊界，因此，選擇反演結果磁化率大于2500×10-6CGSM部分代表巖體分布范圍.

圖6 九瑞礦集區1∶5萬航磁化極匹配濾波剩余淺源異常Fig.6 Matched filter and reduction to pole aeromagnetic residual anomaly（scale 1／50000）map of JiuRui ore district

圖8是根據磁化率推測的九瑞礦集區巖體分布三維圖，其形態反映了巖體的三維展布.從圖中可以看出，反演所得出巖體位置基本能與地質露頭對應，中部地區以中小型巖體成帶出現為主要特征，只有城門山巖體規模較大，北部和西北部地區兩個巖體相對區內其他巖體其規模巨大.

4.3 2.5維人機交互反演

在3D反演的基礎上，為了能對研究區地質結構做出更加細致的刻畫，又結合平面地質圖選取了8條剖面，采用人機聯合2.5維反演的方法［16］，對剖面位置的地層及結構做出更加詳盡的描述.剖面選取的位置如圖9中紅色虛線所示（剖面自西向東依次編號為L1－L8）.

剖面初始模型的構建主要是依據地震剖面、三維反演結果和地質剖面并結合所收集的地層、地質資料來完成的.雖然剖面的初始模型一部分來源于三維反演結果，但剖面擬合結果對各地層特征的刻畫比起三維反演的結果要顯得更加的具體和直觀，剖面能直接反映出各地層的厚度和構造形式，向斜、背斜、斷層等等構造形式能在剖面上一目了然地顯示出來.同時對剖面的反演與三維反演結果之間能形成相互約束的關系.即三維反演結果能說明剖面擬合的正確性，而通過剖面的擬合我們也能判斷三維反演結果的可靠性.現給出模型基本參數如下：模型厚度5000m，背景密度采用2.69g／cm3，深度范圍為地下0～2000m，為了消除邊界的影響，將模型體向南東和北西兩個方向各延長6km.在此模型的基礎上依據所收集的地質資料和其他物探資料的約束，反復修改模型所得出的部分最終地質剖面如圖10.

圖7 九瑞礦集區航磁3D反演成像結果（上層為航磁化極異常，藍色到紅色表示磁異常從低到高變化；中間層為地質圖；下層為磁化率切片，藍色到紅色表示磁化率從低到高漸變）Fig.7 3Daeromagnetic inversion imaging result of JiuRui ore district（Upper：Aeromagnetic anomaly，color change from blue to red shows change of magnetic anomaly from low to high；Middle：Geological map；Lower：Section of susceptibility，changes from blue to red shows the gradual change of susceptibility from low to high）

4.4 邊緣檢測

九瑞地區成礦受到基底斷裂、蓋層斷裂控制，同時也受到層位和巖體的控制.這四種控制成礦作用的機制中，斷裂和巖體邊界具有明顯的線性特征，而重力、航磁邊緣檢測能為我們識別并提取重磁數據中豐富的線性特征信息，因此，本文對九瑞地區重磁數據進行了多尺度邊緣檢測，得到重磁邊緣檢測結果（如圖11、12、13）.

段主任說：我不是說了嗎，像他這種情況，必須住院做個氣管鏡才能知道。你們商量一下，看看怎么辦吧？王姐又問：他們回去做行嗎？段主任點頭說：氣管鏡是小手術，在哪兒做都行。王姐說：段主任，咱們商量一下唄，在這兒做，不住院行不行？段主任笑了，擺手說：這可不行，做氣管鏡，必須住院。再說做前還有許多化驗要做呢。

圖8 0～4000m深度航磁3D成像反演推測巖體圖Fig.8 Conjectural rock mass（0～4000m）through the result of 3Dmagnetic properties inversion

圖9 人機互動反演剖面位置圖Fig.9 Human－computer interaction inversion profile location

4.4.1 布格重力異常多尺度邊緣檢測

布格重力異常多尺度邊緣檢測結果反映的是不同深度密度差異地質體的邊界，包括斷裂、板塊分界線、盆地邊界及巖體接觸線等.根據九瑞地區地質特點，判斷本文重力異常多尺度邊緣檢測結果主要表示基底斷裂、蓋層斷裂及部分巖體分界線.根據多尺度檢測結果解釋的原則［13，17］，本文在同一圖幅內，將不同延拓高度檢測結果疊加到一起，采用不同顏色、大小的線型表示不同深度尺度的信號邊界形跡，根據不同深度邊界信號形跡在分布位置上所體現的相似性來獲知該構造形跡的發育深度及傾向特征.

對九瑞礦集區1∶5萬布格重力異常向上延拓，延拓高度自500m開始，以500m為間隔依次增加，最大上延高度為6km.將每個高度的邊界檢測值疊加到一起，形成九瑞礦集區1∶5萬布格重力數據多尺度邊緣檢測結果（如圖11）.圖中藍色線到綠色線代表了不同的延拓高度：從500m～6km.依同樣方法，得到九瑞礦集區1∶20萬布格重力數據多尺度邊緣檢測結果（如圖12，最小延拓高度為5km，最大延拓高度為35km，間隔為1km，最大延拓高度結果用紅色線表示，最小高度用藍色線表示）.

4.4.2 航磁異常多尺度邊緣檢測

圖10 模型剖面圖（a）L1線剖面；（b）L3線剖面；（c）L8線剖面.橫軸為水平距離，單位m，縱軸上部為重力異常，單位mGal，下部為深度，單位m.Fig.10 The final model profile

圖11 1∶5萬布格重力異常數據多尺度邊緣檢測結果Fig.11 WORMS result of Bouguer gravity anomaly（scale 1／50000）

圖12 1∶20萬布格重力異常數據多尺度邊緣檢測結果Fig.12 WORMS result of Bouguer gravity anomaly（scale 1／200000）

航磁異常的多尺度邊緣檢測結果反映的是磁性邊界.如果重力多尺度邊緣檢測結果與航磁多尺度邊緣檢測結果吻合，則說明反映的邊界為重磁同源，如有中基性巖漿侵入的深大斷裂［18］.

依重力邊緣檢測結果的形成方法，本文得到九瑞礦集區1∶5萬航磁數據多尺度邊緣檢測結果（如圖15，最小延拓高度為1km，最大延拓高度為10km，間隔為1km），并將該結果與航磁3D反演成像結果疊合在一起（如圖13）.

5 結果及地質解釋

5.1 九瑞礦集區三維地質模型

將擬合好的剖面插入三維反演得出的高密度體中（如圖14），圖中三維反演所得出的高密度體清晰地體現了蓋層的起伏.褶皺構造是本研究區控礦構造之一，受褶皺翼部層滑斷裂帶控制的礦體在次級小背斜的軸部礦體厚度會增大.由于巖層褶皺發生層滑時，易在鞍部形成層間空隙，造成虛脫現象，從而為礦液提供了有利空間，而導致礦體厚度在該部位加大，從而形成具有開采價值的礦點、礦床.重力3D反演推測的高密度體為圈定這些構造提供了直接信息.配合一定數量的擬合剖面，我們就能很快地做到某些層位或巖體的空間定位.

圖13 航磁數據多尺度邊緣檢測結果Fig.13 WORMS result of aeromagnetic anomaly（scale 1／50000）

圖14 三維反演高密度體（黃色體）與模擬剖面的空間結構圖Fig.14 The result（yellow body）of 3Dgravity properties inversion and 8geologic profiles penetrated it

圖15 反演結果密度差＞0.015g／cm3模型體與地質圖（置于-1km）疊合圖黃色為高密度體，圖中圓點為礦點位置，其中紅、黃、藍、墨綠、綠分別表示新發現礦化點、銅金礦點、銅多金屬礦點、金礦點、鐵礦.Fig.15 The model of inversion result density difference＞0.015g／cm3Yellow for the high－density body，dot for the mine site location，and red，yellow，blue，dark green and green respectively delegating newly discovered mineralization points，Cu－Au mines，Cu polymetallic mines，Au mines and Fe mines.

圖16 -1250m深度巖體形態推測（棕紅色體）與地質圖（30%透明）疊合Fig.16 Conjectural rock mass（0～1250m）through the result of 3Dmagnetic properties inversion overlying geological map（30%transparent）

采用等值面（IsoSurface）技術，保留反演結果密度差大于0.015g／cm3的部分，將其與地質圖，礦床（點）分布集成于三維環境中（圖15）.圖中EW向和近EW向的三個高密度體正好對應了楓林大橋背斜，大浪—丁家山背斜和吳家壟北—長嶺山—城門山背斜，其間凹兜對應為上田（新莊）—黃橋復式向斜.從反演結果看，蓋層褶皺淺部地層較陡，向深部逐漸趨于平緩.背斜較緊密，向斜較寬緩，略具梳狀褶皺的特點.而這一結果在模擬剖面上也是一目了然.結合剖面擬合時的情況，向斜部位的石炭系、二疊系、三疊系灰巖密度大于背斜部位的志留系砂頁巖，僅以高密度的奧陶系地層分布難以合理解釋整個區域的重力高、低值帶.因此推測深部高密度基底與蓋層褶皺有相似的褶皺狀態.

結合礦床和礦點在礦集區的分布位置（圖中藍色和紅色點代表礦點）可以發現，在九瑞礦集區，礦床（點）多分布于高密度體上，且大多數都處于高密度體的邊緣部位.只有極個別分布于低密度區.

將地質圖與推測巖體結果在三維可視化環境中疊合在一起，截取不同深度水平切片，能更加直觀地反映巖體在不同深度的形態，同時，也便于分析礦床（點）與巖體分布之間的關系，為成礦研究和找礦分析提供豐富的信息.圖16是-1250m深度巖體與地質圖疊合到一起的三維透視圖，圖中圓點為礦點位置.從圖中可以看出，航磁3D反演推測巖體與地質填圖中巖體基本對應.

九瑞礦集區是以銅、金為主的礦集區，其成礦物質主要來源巖體深部（幔源）且淺部就位，即深部隱伏巖基控制了礦田分布，其上的高侵位小巖體對礦床的形成具有明顯的控制作用，反映出本區巖漿演化和巖體前緣成礦的基本規律.對照圖12，我們很容易看到區內已知礦床和礦點（如城門山、武山）均位于巖體上方或巖體外側.結合前人在九瑞地區找礦所總結的規律，可以認為有礦必有巖體.因此，航磁3D反演推斷的巖體（包括隱伏巖體）是找礦預測的重要依據.從巖體推斷圖（圖8，16）發現，九瑞礦集區武山，城門山，東雷灣—鄧家山3處隱伏巖體規模較大，其余巖體分布較散亂，沒有形成規模.這一方面說明九瑞地區對應的巖漿巖規模較小，巖體向下延伸有限，其離巖漿活動中心稍遠.但從另一方面看，東雷灣—鄧家山處存在的超大隱伏巖體，與九瑞地區最著名的兩個礦床（城門山、武山）所處的巖體規模相當，其行成了區內西北部磁異常的抬升區.并且其平面所在地理位置也在城門山—武山豐山洞隱伏斷裂附近，因此認為該處有比較好的深部找礦前景.

5.2 九瑞礦集區線性結構

九瑞地區密度邊界大多以NEE或近EW向出現（如圖11，12）.中淺部密度邊界線分布復雜，線束集中；中深部密度邊界線分布簡單，線束較稀疏.總體來看，區內密度邊界線以NEE向線束的展布最為常見，該組線束與褶皺軸向基本一致，主要是研究區內蓋層斷裂的顯示.從線束的位置分布看，蓋層斷裂多發育于背向斜的翼部，具體為大橋—寶山復式背斜、丁家山—大沖背斜和城門山—長山背斜的兩翼；從線束的形態看，蓋層斷裂一般切割較深，多數切割深度已經超過蓋層厚度；從線束的連續性看，蓋層斷裂的規模較大，斷裂長度最大可達30多公里.其余規模較小，連續性不強的線束，位置形態比較雜亂，且延伸深度不大，在6km以后全部消失.其顯示了九瑞地區地表復雜的構造情況.除此之外，如圖12中有一條從研究區的北東方向穿入，走向為NEE向，經城門山、瑞昌后轉為NWW向的線束.線束形態稀疏，連續性好，橫穿整個研究區，且線束的延伸非常深，直至上延30km仍有顯示，初步判斷為長江深斷裂的顯示.該斷裂切割深度大，斷裂面傾向較緩，且斷裂深、淺部位走向表現出不一致性，是九瑞礦集區重要的導礦構造，與巖漿的上侵和礦床的形成關系密切，為區內礦床的形成提供了深部礦源，控制了整個研究區的礦床形式的格局.

從圖中可以看出，磁性邊界與航磁3D反演成像結果中高磁性體邊界吻合度較好.在丁家山—大沖背斜、城門山—長山背斜和新壙向斜兩翼處磁性邊界與密度邊界也有較好的重合.而此處空間位置距長江深斷裂較近，因此推測此重磁邊界線重合與中基性巖漿侵入有較大的關系.

6 討論與結論

重磁3D反演結果及重磁異常數據邊緣檢測結果從地球物理的角度更加形象地描述了九瑞地區深部構造特征及其對成礦控礦的指示：

（1）九瑞礦集區自北向南，高密度體相間出現，從平面上看，其形態如條帶狀.結合地質資料發現，九瑞地區密度的條帶狀分布形式與區內褶皺構造的形式一一對應.自北向南高低密度帶分別對應通江嶺—鄧家山復式向斜，大橋—寶山復式背斜，黃橋復式向斜，丁家山—大沖背斜，塞城湖—烏石街復式向斜和城門山—長山背斜，新壙向斜.高密度體對應的背斜除反應了密度相對高的奧陶系、志留系地層外，同時也反應了其下部高密度基底的上隆，這說明九瑞地區基底與蓋層有相似的褶皺形式，蓋層的褶皺形式是基地褶皺形式的一種繼承.

（2）九瑞礦集區南部高磁性體的分布與長江深斷裂的走向相一致，推測為巖漿物質沿長江深斷裂上涌的反應，長江深斷裂為其形成提供了深部的物質源.研究區西北角高磁性體所對應的空間位置處于深部長江深斷裂與中淺部蓋層斷裂交叉處.因此，推測該磁性體的形成方式為：長江深斷裂為其形成提供深部的物質源，蓋層斷裂為巖漿的上涌提供上行通道.研究區北部高磁性體雖不處于長江深斷裂之上，但是其所處位置的蓋層斷裂切割較深，延伸較長，斷裂深部也可能與長江深斷裂相交.因此，推測該磁性體的形成方式為：長江深斷裂為其形成提供深部的物質源，蓋層斷裂為巖漿的上涌提供上行通道.

（3）從重磁多尺度邊緣檢測結果中找不到前人指出的城門山—武山—豐山洞隱伏深斷裂.而對于長江深斷裂，重新厘定的位置及空間展布形態也與前人研究結果有一定出入：長江深斷裂穿過本研究區段不似前人描述的為單一的NWW走向，而是自西向東由NWW走向轉為NEE走向進而穿出研究區，且斷裂深部與淺部走向也并不一致：斷裂在淺部由曹家沖以北穿出研究區，深部則在穿過瑞昌之后向豐山洞方向延伸，從研究區西北角穿出.

（4）九瑞礦集區已知礦點分布規律明顯：已知礦點分布伴隨巖體出現，空間位置相對集中，80%以上分布于NEE或近EW向的蓋層斷裂附近，斷裂走向與褶皺軸向平行，并且展布于褶皺背斜（重力3D成像反演結果中顯示為高密度體）的兩翼.礦點分布特征充分說明了巖體、斷裂和褶皺對礦點的控制關系，這也為今后的找礦提供了一個明確的方向.

（References）

［1］常印佛，劉湘培，吳言昌.長江中下游銅鐵成礦帶.北京：地質出版社，1991：1－147.Chang Y F，Liu X P，Wu Y C.Iron－Copper Mineralization Belt of Middle－Lower Yangtze Valley（in Chinese）.Beijing：Geological Publishing House，1991：1－147.

［2］Pan Y M，Dong P.The Lower Changjiang（Yangzi／Yangtze River）metallogenic belt，east central China：intrusion－and wall rock－hosted Cu－Fe－Au，Mo，Zn，Pb，Ag deposits.Ore Geology Reviews，1999，15（4）：177－242.

［3］趙文津.長江中下游金屬礦找礦前景與找礦方法.中國地質，2008，35（5）：771－802.Zhao W J.Ore prospects and ore exploration methods for metal deposits in the middle and lower Yangtze River valley.Geology in China （in Chinese），2008，35（5）：771－802.

［4］Dong S W，Xiang H S，Gao R，et al.Deep structure and ore formation within Lujiang－Zongyang volcanic ore concentrated area in Middle to Lower Reaches of Yangtze River.Acta Petrologica Sinica，2000，26（9）：2529－2542.

［5］滕吉文.強化開展地殼內部第二深度空間金屬礦產資源地球物理找礦、勘探和開發.地質通報，2006，25（7）：767－771.Teng J W.Strengthening geophysical exploration and exploitation of metallic minerals in the second deep space of the crustal interior.Geological Bulletin of China （in Chinese），2006，25（7）：767－771.

［6］Li Y G，Oldenburg D W.3－D inversion of magnetic data.Geophysics，1996，61（2）：394－408.

［7］Li Y G，Oldenburg D W.3－D inversion of gravity data.Geophysics，1998，63（1）：109－119.

［8］姚長利，鄭元滿，張聿文.重磁異常三維物性反演隨機子域法方法技術.地球物理學報，2007，50（5）：1576－1583.Yao C L，Zheng Y M，Zhang Y W.3－D gravity and magnetic inversion for physical properties using stochastic subspaces.Chinese J.Geophys.（in Chinese），2007，50（5）：1576－1583.

［9］Hornby P，Boschetti F，Horowitz F G.Analysis of potential field data in the wavelet domain.Geophys.J.Int.，1999，137（1）：175－196.

［10］Archibald N，Gow P，Boschetti F.Multiscale edge analysis of potential field data.Exploration Geophysics，1999，30（2）：38－44.

［11］Horowitz F，Hornby P，Boschetti F.Developments in the analysis of potential field data via multiscale edge representation（in Exploration beyond 2000；conference handbook）Preview 2000，84，77.

［12］Holden D J，Archibald N J，Boschetti F，et al.Inferring geological structures using wavelet－based multiscale edge analysis and forward models.Exploration Geophysics，2000，31（4）：617－621.

［13］Austin J R，Blenkinsop T G.The Cloncurry Lineament：Geophysical and geological evidence for a deep crustal structure in the Eastern Succession of the Mount Isa Inlier.Precambrian Research，2008，163（1－2）：50－68.

［14］王大勇.長江中下游礦集區綜合地質地球物理研究—以九瑞、銅陵礦集區為例［博士論文］.長春：吉林大學，2010.Wang D Y.The integrated geophysical and geological study in the ore belt of the Middle and Lower Reach of the Yangtze River—The cases study of Tongling and Jiurui ore district［Ph.D.thesis］（in Chinese）.Changchun：Jilin University，2010.

［15］Keating P，Pinet N.Use of non－linear filtering for the regional－residual separation of potential field data.Journal of Applied Geophysics，2011，73（4）：315－322.

［16］Enmark T.A versatile interactive computer program for computation and automatic optimization of gravity models.Geoexploration，1981，19（1）：47－66.

［17］管燁.云南三江地區地殼三維構造格架與礦集區關系研究［博士論文］.北京：中國地質科學院，2005.Guan Y.Relationship between 3Dtectonic frame and metallogenic district：Theory and implications in the Three Rivers Region，Yunnan ［Ph.D.thesis］（in Chinese）.Beijing：Chinese Academy of Geological Sciences，2005.

［18］嚴加永.長江中下游成礦帶深部背景綜合地球物理研究［博士論文］.北京：中國科學院地質與地球物理研究所，2010.Yan J Y.Integrated geophysical study for deep background of the Middle and Lower Reach of the Yangtze River metallogenic Belt［Ph.D.thesis］（in Chinese）.Beijing：Institute of Geology and Geophysics，CAS，2010.