摘要：為直觀展現(xiàn)戈塘金礦礦體空間展布規(guī)律并進(jìn)行綜合成礦分析，本文基于現(xiàn)有鉆孔與剖面等地質(zhì)資料，結(jié)合已有成礦物質(zhì)來源等研究結(jié)果，利用Petrel軟件建立了戈塘金礦三維地質(zhì)模型，在此基礎(chǔ)上討論了戈塘金礦體在地下的空間展布規(guī)律與富集特征，分析了該礦床的成礦過程，并總結(jié)了該礦床的成礦模式。研究結(jié)果表明：Au主要富集于龍?zhí)督M1段地層，大多礦體品位較低，達(dá)到工業(yè)品位的礦體多集中于斷層附近，總體上礦體品位與礦體體積呈負(fù)相關(guān)，但高品位礦體多次出現(xiàn)體積的異常增加；礦床順地層展布特征顯著，受不整合面或?qū)娱g滑脫面控制明顯，但斷層可能為含礦熱液提供了上升通道；礦床形成與巖漿活動(dòng)有較大關(guān)系，熱液沿?cái)鄬舆M(jìn)入含礦層并發(fā)生強(qiáng)烈蝕變，使Au析出預(yù)富集，在后期地下水等作用下形成金礦床，為熱液運(yùn)移-流體改造的成礦模式。上述研究結(jié)果對(duì)我國卡林型金礦成因研究及找礦工作具有一定的指示及參考意義。

關(guān)鍵詞：戈塘金礦；三維地質(zhì)建模；卡林型金礦；成礦分析

中圖分類號(hào)：P618.51文獻(xiàn)標(biāo)志碼：A卡林型金礦（Carlin-type gold deposits）最早發(fā)現(xiàn)于美國內(nèi)華達(dá)州卡林地區(qū)并因此得名［1-2］，特點(diǎn)是礦體浸染，賦存于碳酸鹽巖圍巖中［3-4］?！暗?黔-桂”金三角作為我國重要的金礦富集地之一，探明儲(chǔ)量超800 t ［5］。位于黔西南的戈塘金礦是金三角重要的組成部分［6］。前人對(duì)黔西南卡林型金礦成礦機(jī)制、成礦流體特征與流體來源、成礦年代等進(jìn)行過大量研究［7-9］。以上研究多基于樣品試驗(yàn)分析方法與二維地質(zhì)資料等分析討論，對(duì)礦區(qū)三維地質(zhì)特征缺乏深入的認(rèn)識(shí)，沒有形成立體的三維地質(zhì)格架。由于戈塘金礦開采條件極為復(fù)雜，現(xiàn)有的鉆孔及剖面等地質(zhì)資料往往因?yàn)椴粔蛑庇^而無法及時(shí)準(zhǔn)確地應(yīng)用于開采實(shí)踐中，或者導(dǎo)致開采過程中的資源浪費(fèi)等問題。因此，如何實(shí)現(xiàn)對(duì)礦床更為直觀、準(zhǔn)確地認(rèn)識(shí)成為了亟待解決的問題。

三維地質(zhì)建模的概念由加拿大地質(zhì)學(xué)家HOULDING ［10］提出，即在三維環(huán)境下，將空間信息管理、地質(zhì)解譯、空間分析和預(yù)測(cè)、地學(xué)統(tǒng)計(jì)、實(shí)體內(nèi)容分析以及圖形可視化等工具結(jié)合起來，并用于地質(zhì)分析的技術(shù)［11-13］。近年來，計(jì)算機(jī)技術(shù)及地理信息科學(xué)等的發(fā)展，開啟了地質(zhì)大數(shù)據(jù)時(shí)代［14］，數(shù)字地球［15］、數(shù)字礦山［16］、玻璃地球［17］等的提出指明了未來地質(zhì)行業(yè)的發(fā)展趨勢(shì)。三維地質(zhì)建模是實(shí)現(xiàn)以上理論重要的技術(shù)手段之一。由斯倫貝謝公司開發(fā)的Petrel軟件具有模塊化、高精度、高集成等優(yōu)點(diǎn)，被大量學(xué)者應(yīng)用于油氣儲(chǔ)層模擬［18］、第四系沉積物垂向分布規(guī)律研究［19］、儲(chǔ)層流體模擬［20］、鉆孔設(shè)計(jì)［21］等方面；在固體礦產(chǎn)方面，Petrel也被大量應(yīng)用于地質(zhì)環(huán)境表征等研究［22］。

因此，本文利用Petrel軟件，以戈塘金礦接替資源勘查區(qū)為研究對(duì)象，以現(xiàn)有鉆孔及剖面數(shù)據(jù)為基礎(chǔ)，建立戈塘金礦接替資源勘查區(qū)三維地質(zhì)模型，并進(jìn)行綜合成礦分析，結(jié)合礦區(qū)地質(zhì)資料討論礦床成礦模式，對(duì)我國卡林型金礦成因研究及找礦工作具有一定的指示及參考意義。

1地質(zhì)背景

1.1區(qū)域地質(zhì)背景

研究區(qū)大地構(gòu)造位置處于揚(yáng)子準(zhǔn)地臺(tái)西南緣的黔北臺(tái)隆六盤水?dāng)嘞葜校▓D1（a）），區(qū)內(nèi)主體構(gòu)造為戈塘穹窿背斜［1， 23］，具弧形轉(zhuǎn)折走向多變的特點(diǎn)［24］；斷裂構(gòu)造發(fā)育，主體構(gòu)造線為北東向，多為高角度正斷層［6］，且切割含礦層或破壞礦體。構(gòu)造區(qū)域內(nèi)沉積地層較為發(fā)育，出露地層以二疊系茅口組、龍?zhí)督M及長興大隆組與三疊系夜郎組及永寧鎮(zhèn)組為主［6］，巖性以灰?guī)r、砂巖及黏土巖等為主［25］?；鸪蓭r出露較少［25-26］，只在局部有峨眉山溢流玄武巖和堿性超鎂質(zhì)巖出露［6， 27］。然而重磁勘探資料表明在深部有隱伏巖漿巖侵入［25， 28-29］。區(qū)內(nèi)礦產(chǎn)以金礦與煤礦為主，金礦多集中于茅口組與龍?zhí)督M地層接觸帶，煤礦多集中于龍?zhí)督M地層中［6］ 。研究區(qū)大地構(gòu)造演化經(jīng)歷了拉張裂陷、擠壓碰撞、伸展裂陷等過程［4， 30］，此期間區(qū)內(nèi)深大斷裂形成超基性巖脈侵入體，利于卡林型金礦成礦［29-30］。盆地裂陷過程中發(fā)育的正斷層也為含礦流體上升提供了通道［4］。黔西南大地構(gòu)造位置獨(dú)特，構(gòu)造演化復(fù)雜多樣，最終形成了良好的卡林型金礦成礦條件。

1.2礦區(qū)地質(zhì)背景

戈塘金礦位于黔西南戈塘穹窿東南部（圖1（b）），主要出露地層為二疊系茅口組與龍?zhí)督M；其中，茅口組為灰色厚層至塊狀灰?guī)r，厚度大于200 m；龍?zhí)督M分5段，地層厚約250 m，平行不整合于下覆茅口組灰?guī)r之上［1， 6］。第1段含金角礫巖層為研究區(qū)金礦產(chǎn)出層位，巖性以黏土角礫巖及硅化灰?guī)r角礫巖為主，厚度不大，橫向變化明顯；第2段為黏土巖與粉砂巖，部分區(qū)域夾煤線；第3段以黏土巖為主；第4段為硅化灰?guī)r或灰?guī)r；第5段為黏土巖、粉砂巖及細(xì)砂巖互層。礦區(qū)褶皺構(gòu)造不發(fā)育；斷裂構(gòu)造發(fā)育，以北東向?yàn)橹鳎?］，對(duì)礦體破壞較小。礦石類型主要為氧化礦石（表1），圍巖蝕變以硅化為主，次為黃鐵礦化，再次為輝銻礦化、硫砷礦化、螢石化、石膏化、黏土化等。（a） 右江盆地構(gòu)造簡圖；（b） 黔西南卡林型金礦區(qū)域地質(zhì)圖

通過統(tǒng)計(jì)分析礦體空間位置、規(guī)模、礦石類型、平均品位、形態(tài)、產(chǎn)狀及賦礦巖性等，發(fā)現(xiàn)礦體平面上主要位于勘探線11-5A之間（圖2），空間上主要集中于龍?zhí)督M1段。為便于理解，本文圖2—圖9中坐標(biāo)系為CGCS2000坐標(biāo)系。為完整建立戈塘金礦礦體三維模型，水平上大致以勘探線11-5A為本次建模邊界（圖2， 3）、（表1）；由于現(xiàn)有鉆孔及剖面資料未探及茅口組底部邊界，因此，垂向上以龍?zhí)督M頂為模型頂部邊界，底部邊界設(shè)定為1 100 m。

2三維地質(zhì)建模

2.1數(shù)據(jù)準(zhǔn)備

地質(zhì)數(shù)據(jù)是三維地質(zhì)建模的基礎(chǔ)，在其它條件都相同的情況下，對(duì)各種條件數(shù)據(jù)的掌握程度（數(shù)量、精度等）就決定了模型的預(yù)測(cè)效果［31］。本研究共收集礦區(qū)相關(guān)資料包括：地質(zhì)勘查報(bào)告2份，地形地質(zhì)圖2份，勘探線剖面圖17份，鉆孔柱狀圖77份；鉆孔取樣分析結(jié)果482個(gè)；利用上述資料統(tǒng)計(jì)地層分層與斷層數(shù)據(jù)若干，將數(shù)據(jù)導(dǎo)入Petrel軟件并建立礦區(qū)地質(zhì)數(shù)據(jù)庫，以備后續(xù)建立三維模型時(shí)快速調(diào)用。

2.2斷層模型

斷層模型可以較為準(zhǔn)確地表征斷層的發(fā)育形態(tài)、斷層與斷層以及斷層與層面的接觸關(guān)系［32］。根據(jù)研究區(qū)地質(zhì)資料，建模區(qū)域內(nèi)共7條斷層，大部分為北東向高角度正斷層，斷層均切穿龍?zhí)督M地層。利用AutoCAD軟件讀取勘探線剖面圖中斷層產(chǎn)狀信息，并以斷層線方式導(dǎo)入Petrel中。將斷層線轉(zhuǎn)換為斷層柱（pillar），利用pillar控制的曲面建立斷層面，將剖面圖導(dǎo)入Petrel以約束斷層產(chǎn)狀，同時(shí)可以通過編輯pillar修改斷層面；建出斷層模型如圖4所示，為便于后續(xù)地層建模中網(wǎng)格化步驟的進(jìn)行，需將不同斷層模型相交處進(jìn)行連接，之后開始建立地層模型。

2.3地層模型

地層模型是對(duì)地層組成、結(jié)構(gòu)等的綜合表現(xiàn)。能夠清晰地展示地層在地下的展布，地層模型的空間網(wǎng)格還是建立屬性模型的前提。研究區(qū)地質(zhì)資料顯示，區(qū)內(nèi)出露地層主要為二疊系茅口組、龍?zhí)督M地層，龍?zhí)督M地層自下而上分為5段（P3l1—P3l5），因此，將模型層位數(shù)設(shè)為6。由于資料中無茅口組底部數(shù)據(jù)，為完整建立龍?zhí)督M一段底部層面，將鉆孔分層數(shù)據(jù)中茅口組底部海拔統(tǒng)一設(shè)定為1 100 m，同時(shí)為模型底部邊界。以精度為20 m×20 m網(wǎng)格化模型，邊界確定以礦體分布范圍為主。建立層面模型后利用先前導(dǎo)入勘探線剖面圖調(diào)整層面斷距，確定模型層面產(chǎn)狀，隨后建立地層空間模型（圖5）。由圖可知，龍?zhí)督M地層褶皺以單斜為主，走向?yàn)楸睎|向，可見小撓曲。為之后建立礦體品位模型，將地層模型在垂向上細(xì)分為1 m厚的小層，便于向網(wǎng)格內(nèi)部填充品位屬性。

2.4品位分布模型

Petrel軟件基于油氣藏開發(fā)，無針對(duì)礦體建模的模塊，但Petrel軟件對(duì)斷層的粗化處理相比較其他建模軟件來說更為精確，更加適合卡林型金礦三維地質(zhì)模擬，能夠更為直觀地展現(xiàn)礦體品位與斷層的關(guān)系。因此，本文選擇參考地質(zhì)模塊中屬性模型的建立方法，以測(cè)得樣品品位為測(cè)井屬性導(dǎo)入Petrel軟件，按照屬性模型建模步驟進(jìn)行測(cè)井曲線粗化與數(shù)據(jù)分析，根據(jù)傳統(tǒng)上礦體邊界品位與工業(yè)品位的定義［33］，將0.5 g/t與1 g/t分別作為礦床邊界品位與工業(yè)品位，品位高于0.5 g/t的區(qū)域定義為礦化體，高于1 g/t的區(qū)域定義為礦體并設(shè)為紅色以篩選網(wǎng)格，最終建立戈塘金礦礦體品位分布模型與礦體品位模型（圖6）。觀察發(fā)現(xiàn)龍?zhí)督M上部地層幾乎無礦體分布，礦體展布與龍?zhí)督M1段地層產(chǎn)狀高度一致，沿地層分布特征明顯；另外礦體富集區(qū)大多靠近斷層，證明斷層可能為成礦熱液上升提供了通道。

3Au富集與展布

戈塘金礦中不同地層Au富集情況差異很大。以礦區(qū)地層模型及品位模型等為基礎(chǔ)，以邊界品位為約束進(jìn)行品位提取，統(tǒng)計(jì)分析Au品位富集規(guī)律。利用軟件中屬性報(bào)告將Au品位網(wǎng)格數(shù)以0.5～13 g/t為邊界，以0.25 g/t為區(qū)間順序提取，得到Au品位分布直方圖（圖7）。經(jīng)統(tǒng)計(jì)，礦區(qū)地層平均Au品位為0.20 g/t，主要含礦層為1.02 g/t。Au品位主要集中于0.5～2 g/t之間。低于工業(yè)品位1 g/t的網(wǎng)格占比69.08%，高于5 g/t的網(wǎng)格占比13.37%。總體上來看，礦體品位越高，對(duì)應(yīng)網(wǎng)格數(shù)量逐漸減少，呈負(fù)相關(guān)關(guān)系。但在直方圖中多次出現(xiàn)較高品位礦體體積異常增加，可能是地層中存在某種作用使得Au出現(xiàn)異常富集。

為更好對(duì)比垂向上各地層Au富集情況， 將各個(gè)層位分別進(jìn)行品位統(tǒng)計(jì)，在龍?zhí)督M3段至龍?zhí)督M5段地層中未發(fā)現(xiàn)含金礦石，因此，將龍?zhí)督M1段與龍?zhí)督M2段進(jìn)行統(tǒng)計(jì)（圖8）。整體上各地層礦體仍以低品位為主，且礦體品位與體積具有負(fù)相關(guān)趨勢(shì)；礦體主要集中于龍?zhí)督M1段，龍?zhí)督M2段僅有少量礦化體。在龍?zhí)督M1段 4.5 g/t及9 g/t處與2段5～9.5 g/t處多次出現(xiàn)Au礦體體積異常增加。龍?zhí)督M1段品位分布較為連續(xù)，礦體具有較好的展布；龍?zhí)督M2段品位分布不連續(xù)，在1.5 g/t處礦體體積急劇減小，證明在該段地層中僅有少量低品位礦體與礦化體存在，且分布不連續(xù)。

整體上看，戈塘金礦床中Au的三維富集與空間展布具有較好的吻合程度，集中表現(xiàn)為低品位礦體居多，異常富集現(xiàn)象頻繁，多富集于龍?zhí)督M1段不整合面及斷層附近等。龍?zhí)督M1段中Au品位遠(yuǎn)高于其他地層，相鄰地層中礦體品位急劇降低，數(shù)量急劇減少，且多沿?cái)鄬痈浇共?，證明不整合面為戈塘金礦床重要的控礦因素，斷層可能為流體上升通道。

4礦床成因探討

4.1成礦物質(zhì)來源

卡林型金礦中Au主要以Au-S絡(luò)合物的形式賦存于成礦流體中［34-36］。在與圍巖含鐵礦物的硫化作用中富集沉淀，是卡林型金礦床金屬沉淀的重要機(jī)制之一［27， 36］。因此，反應(yīng)生成的硫化物中的S能夠指示Au來源，前人在關(guān)于黔西南卡林型金礦Au來源的研究中已經(jīng)做了大量工作［26］，對(duì)戈塘金礦的成礦物質(zhì)來源卻仍有較大爭議，LI等對(duì)礦床樣品進(jìn)行了地球化學(xué)分析，認(rèn)為沉積巖對(duì)戈塘金礦成礦具有主要貢獻(xiàn)［1］；而部分學(xué)者卻通過硫同位素及氦同位素測(cè)試結(jié)果得出戈塘金礦成礦物質(zhì)與巖漿有較大關(guān)系的結(jié)論。例如，由于黃鐵礦為戈塘金礦床中最為主要的含金礦物［23， 27， 35］，在成巖期、成礦早期、成礦晚期等均有生成［29， 36］，因此，HU等根據(jù)黃鐵礦的形態(tài)與共生特征，共鑒別出4期黃鐵礦，并通過主成礦期與成巖期黃鐵礦S同位素值差異較大排除了S來源于沉積物的可能性［6］。再由于本區(qū)內(nèi)尚未發(fā)現(xiàn)區(qū)域變質(zhì)作用蹤跡，由此可推測(cè)與礦物有關(guān)的S可能來源于深部巖漿。盡管在研究區(qū)火成巖出露較少，但證明深部可能有隱伏巖漿巖侵入的重磁勘探資料也為此提供了強(qiáng)有力的證據(jù)［25， 28-29］。對(duì)于黔西南卡林型金礦成礦流體來源前人已做過大量研究，認(rèn)為該地區(qū)卡林型金礦具有相似的成礦流體特征［3， 23］。CHEN等流體包裹體數(shù)據(jù)表明，盆地成礦流體以中低溫為主（113～294 ℃），鹽度低（0.18～13.18 wt%） 且呈弱酸性［37］，LI等與CHEN等還發(fā)現(xiàn)，由早期至晚期成礦，盆地流體溫度、鹽度與氧逸度fO2均降低，pH值升高［37-38］；LIU等則通過紫木凼礦床地球化學(xué)分析結(jié)果推斷成礦流體富含H2S[39]。對(duì)于成礦流體來源大量學(xué)者形成了幾種主要認(rèn)識(shí)包括：部分學(xué)者通過流體包裹體中氫-氧同位素分析判斷為大氣降水或變質(zhì)流體［40-41］；PENG等測(cè)試了流體包裹體中碳-氧同位素特征，認(rèn)為成礦流體主要來源于盆地［42］；然而ZHENG等通過金礦中Hg與As和Sb的強(qiáng)相關(guān)性排除了變質(zhì)流體的可能，并借助硫同位素分析判斷右江盆地卡林型金礦成礦流體可能為巖漿成因［43］；XIE等也通過測(cè)試碳酸鹽巖中Mg同位素值確定了貴州卡林型金礦的巖漿成因［44］。事實(shí)上，JIN等通過氦-氖-氬同位素測(cè)試，發(fā)現(xiàn)氖和氬同位素組成與空氣飽和水大致相當(dāng)，而氦-氖同位素?cái)?shù)據(jù)則位于地殼端元附近［25］，這證明成礦流體并非某一單獨(dú)來源，而是不同時(shí)期多種流體混合來源。綜上所述，成礦流體主要為上升巖漿流體與沉積孔隙流體的混合流體，晚期成礦流體主要為沉積孔隙流體與淺層大氣地下水的混合流體［25， 41， 45］。

4.2主要控礦因素

為更加清晰地認(rèn)識(shí)戈塘金礦床的控礦因素，將戈塘金礦各個(gè)層位礦體網(wǎng)格分別顯示并以不同顏色區(qū)分不同品位礦體，得到戈塘金礦Au品位分布三維模型（圖6（b）， （c））。另外，由于礦區(qū)斷層多為北東走向，因此，在北西-南東向建立多條近似垂直斷層走向的剖面線以進(jìn)一步研究礦體展布規(guī)律（圖9）。礦體品位模型與礦體剖面顯示戈塘金礦床中主要賦礦地層為龍?zhí)督M與茅口組間不整合面，礦體對(duì)于層位的選擇十分明顯，龍?zhí)督M3段至龍?zhí)督M5段地層僅有極少低品位礦化體零散分布。龍?zhí)督M2段地層中Au僅在礦區(qū)東北側(cè)及東南側(cè)沿層面于斷層附近少量富集，且品位較低。龍?zhí)督M1段地層Au礦體沿不整合面較為連續(xù)分布，Au富集與地層關(guān)系明顯，然而幾處礦體區(qū)域均靠近斷層，證明含礦熱液極可能沿礦區(qū)斷層向上運(yùn)移，進(jìn)入茅口組與龍?zhí)督M間地層接觸面成礦。少部分沿?cái)鄬永^續(xù)向上運(yùn)移進(jìn)入上覆地層，使得龍?zhí)督M2段中斷層附近存在少量Au富集。

從巖石能干性差異上，也能夠解釋Au在龍?zhí)督M1段的富集現(xiàn)象。能干性指巖石強(qiáng)度的大小，反映巖石變形的難易程度［46］。龍?zhí)督M1段主要由粉砂巖、泥灰?guī)r、黏土巖、硅化灰?guī)r等組成，相較下伏茅口組塊狀灰?guī)r能干性較弱［47］，地層軟硬相間，滑脫面發(fā)育［48］，含礦熱液易沿滑脫面進(jìn)入龍?zhí)督M1段；而上覆龍?zhí)督M2段黏土巖則成為封堵熱液的良好頂板［49］，使得大部分熱液不再繼續(xù)上升而成礦。綜上所述，戈塘金礦床形成與龍?zhí)督M及茅口組間地層接觸面關(guān)系密切，Au富集主要受不整合面或?qū)娱g滑脫面控制。

4.3成礦過程

龍?zhí)督M及茅口組間不整合面對(duì)戈塘金礦體的形成起到了最主要的控制作用，斷層作為成礦流體遷移通道也產(chǎn)生了一定的影響。根據(jù)戈塘金礦三維地質(zhì)模型及礦體品位展布統(tǒng)計(jì)分析結(jié)果，結(jié)合戈塘金礦床成礦流體遷移研究成果，推測(cè)Au富集過程為：

中二疊世晚期，茅口組灰?guī)r受構(gòu)造運(yùn)動(dòng)影響抬升遭受侵蝕，與上覆龍?zhí)督M地層間形成不整合面后在礦床北西部發(fā)生了大面積巖漿活動(dòng)，巖漿氣液中含有大量H2O、CO2、H2S、HF、Au、As、Sb、Hg等組份，它們以鹵化物、絡(luò)合物或揮發(fā)性氣體化合物形式，沿構(gòu)造裂隙、裂縫源源不斷進(jìn)入地塹盆地，到達(dá)不整合面后由于上部黏土巖地層封堵，大部分熱液進(jìn)入含礦層，造成圍巖強(qiáng)烈蝕變和析出大量巖石中K、Na、Ca、Mg等堿金屬，改變了水介質(zhì)成份，使其成為堿性還原環(huán)境，促使火山氣液迅速冷卻，金屬硫化物快速析出和沉淀。由圖5（c）中可知，龍?zhí)督M2段中礦體多于斷層附近富集，推測(cè)極少部分熱液沿?cái)鄬永^續(xù)向上遷移，進(jìn)入上部地層成礦。隨后含礦層中的Au、Sb等元素在溫度、壓力升高和地下水參與下發(fā)生遷移、分解、析出、沉淀，在有利部位富集成礦。因此，研究區(qū)礦床為熱液運(yùn)移-流體改造形成的金礦床（圖10）。

5結(jié)論

（1）通過建立戈塘金礦三維地質(zhì)模型并統(tǒng)計(jì)分析，發(fā)現(xiàn)Au主要富集層位為龍?zhí)督M1段，大多礦體品位較低，達(dá)到工業(yè)品位的礦體多集中于斷層附近，總體上礦體品位與礦體體積呈負(fù)相關(guān)，但多次出現(xiàn)較高品位礦體的異常體積增加現(xiàn)象。

（2）戈塘金礦床與地層關(guān)系密切，礦體集中于龍?zhí)督M底部，受不整合面或?qū)娱g滑脫面控制明顯，順層展布特征顯著，但斷層可能為含礦熱液提供了向上運(yùn)移的通道。

（3）戈塘金礦的形成與巖漿活動(dòng)有較大關(guān)系，熱液沿?cái)鄬舆M(jìn)入含礦層并發(fā)生強(qiáng)烈蝕變，蝕變類型以硅化為主，使Au析出預(yù)富集，在后期地下水等作用下形成金礦床，為熱液運(yùn)移-流體改造的成礦模式。參考文獻(xiàn)：

[1]LI X， ZHU X， GU Y， et al. Mineralogy and geochemistry characteristics and genetic implications for stratabound carlin-type gold deposits in Southwest Guizhou， China ［J］. Journal of Nanoscience and Nanotechnology， 2017， 17（9）： 6307-6317.

［2］ YAN F Z， LI Q Z. Yangshan gold deposit： the largest Carlin and Carlin-like type gold deposit in China［J］. Acta Geologica Sinica-English Edition， 2008， 82（4）： 804-810.

［3］ HU X L， DING Z J， GONG Y J， et al. Sediment-hosted disseminated gold deposits in orogenic belts： an example from the giant Jiagantan gold deposit in the West Qinling orogen， China［J］. Ore Geology Reviews， 2022， 146（1）： 104950.1-104950.18.

［4］ CHEN M H， MAO J W， BIERLEIN F P， et al. Structural features and metallogenesis of the Carlin-type Jinfeng （Lannigou） gold deposit， Guizhou Province， China［J］. Ore Geology Reviews， 2011， 43（1）： 217-234.

［5］ LI P Y， GUAN S J， HU Y Z， et al. Origin and re-Os geochronology of vein bitumen in the Nanpanjiang basin， SW China： implication for the ore-forming age of Carlin-type gold deposits［J］. Ore Geology Reviews， 2022， 149（1）： 105118.1-105118.15.

［6］ HU X， GONG Y， ZENG G， et al. Multistage pyrite in the Getang sediment-hosted disseminated gold deposit， southwestern Guizhou Province， China： insights from textures and in situ chemical and sulfur isotopic analyses［J］. Ore Geology Reviews， 2018， 99（1）：1-16.

［7］ 聶愛國. 黔西南卡林型金礦的成礦機(jī)制及成礦預(yù)測(cè)［D］. 昆明： 昆明理工大學(xué)， 2007.

［8］ 吳松洋， 侯林， 丁俊， 等. 黔西南卡林型金礦礦田控礦構(gòu)造類型及成礦流體特征［J］. 巖石學(xué)報(bào)， 2016， 32（8）： 2407-2424.

［9］ SU W C， HU R Z， XIA B， et al. Calcite Sm-Nd isochron age of the Shuiyindong Carlin-type gold deposit， Guizhou， China［J］. Chemical Geology， 2009， 258（3/4）： 269-274.

［10］HOULDING S W. 3D Geoscience Modeling［M］. Heidelberg： Springer， 1994.

［11］劉少華， 程朋根， 陳紅華. 三維地質(zhì)建模及可視化研究［J］. 桂林工學(xué)院學(xué)報(bào)， 2003， 23（2）： 154-158.

［12］方海東， 劉義懷， 施斌， 等. 三維地質(zhì)建模及其工程應(yīng)用［J］. 水文地質(zhì)工程地質(zhì)， 2002， 29（3）： 52-55.

［13］李響. 三維地質(zhì)建模技術(shù)的研究［D］. 合肥： 合肥工業(yè)大學(xué)， 2008.

［14］譚永杰， 劉榮梅， 朱月琴， 等. 論地質(zhì)大數(shù)據(jù)的特點(diǎn)與發(fā)展方向［J］. 時(shí)空信息學(xué)報(bào)， 2023， 30（3）： 313-320.

［15］李琦， 吳少巖. 數(shù)字地球：人類認(rèn)識(shí)地球的第三次飛躍［M］.北京： 北京大學(xué)出版社， 1999.

［16］吳立新， 殷作如， 鄧智毅， 等. 論21世紀(jì)的礦山： 數(shù)字礦山［J］. 煤炭學(xué)報(bào)， 2000， 25（4）： 337-342.

［17］劉樹臣. 發(fā)展新一代礦產(chǎn)勘探技術(shù)： 澳大利亞玻璃地球計(jì)劃的啟示［J］. 地質(zhì)與勘探， 2003， 39（5）： 53-56.

［18］SOLOMON O O， JOHN A A， PETER O S. Distribution of sedimentary formations in the Bornu Basin， Nigeria［J］. Petroleum Exploration and Development， 2015， 42（5）： 674-682.

［19］AI X， SUN B， CHEN X. Integrated 3D modeling of Quaternary sediments in the Beijing Plain， based on a sequential indicator simulation［J］. Geologia Croatica， 2019， 72（1）： 3-17.

[20]LI P， LU J， ZHOU D， et al. A preliminary simulation of CO2-EOR and storage in one heavy oil carbonate oilfield offshore Guangdong， China [C]//The 8th International Congress on Environmental Geotechnics （ICEG）. Singapore： Springer， 2019， 3-16.

[21]MAHMUD N， LIU X. Optimizing well structures and trajectories for maximizing oil recovery [C]//The 4th International Conference on Integrated Petroleum Engineering and Geosciences （ICIPEG）. Singapore： Springer， 2016， 783-793.

[22]DENT J， ROTTING T， WILLIAMS M. 3D geological modelling for geo-environmental characterization of mineral deposits and pragmatic management of geochemical risks [C]//The Annual Meeting of the International-Mine-Water-Association （IMWA）. Leipzig： Tech Univ Bergakademie Freiberg， 2016， 1246-1248.

［23］HUANG Y， HU R， BI X， et al. Low-temperature thermochronology of the Carlin-type gold deposits in Southwestern Guizhou， China： implications for mineralization age and geological thermal events［J］. Ore Geology Reviews， 2019， 115（1）： 103178.1-103178.14.

［24］田志偉. 貴州安龍縣戈塘金礦床地質(zhì)特征及成因探討［J］. 世界有色金屬， 2017， （04）： 141-142.

［25］JIN X Y， HOFSTRA A H， HUNT A G， et al. Noble gases fingerprint the source and evolution of ore-forming fluids of Carlin-type gold deposits in the golden triangle， south China［J］. Economic Geology， 2020， 115（2）： 455-469.

［26］WEI D T， ZHOU T F， XIA Y， et al. Ore fluid origin recorded by apatite chemistry： a case study on altered dolerite from the Badu Carlin-type gold deposit， Youjiang basin， SW China［J］. Ore Geology Reviews， 2022， 143（1）： 104745.1-104745.17.

［27］LIN S R， HU K， CAO J， et al. An in situ sulfur isotopic investigation of the origin of Carlin-type gold deposits in Youjiang Basin， southwest China［J］. Ore Geology Reviews， 2021， 134（1）： 104187.1-104187.20.

［28］WANG G， ZHU Y， ZHANG S， et al. 3D geological modeling based on gravitational and magnetic data inversion in the Luanchuan ore region， Henan Province， China［J］. Journal of Applied Geophysics， 2012， 80（1）： 1-11.

［29］JIN X Y， YANG C F， LIU J Z， et al. Source and evolution of the Ore-forming fluids of Carlin-type gold deposit in the Youjiang Basin， South China： evidences from Solute Data of Fluid Inclusion Extracts［J］. Journal of Earth Science， 2021， 32（1）： 185-194.

［30］宋威方. 南盤江-右江卡林型金礦多層次構(gòu)造滑脫成礦系統(tǒng)地球化學(xué)研究［D］. 貴陽： 貴州大學(xué)， 2022.

［31］張國強(qiáng)， 于作剛. 利用Petrel軟件建立油藏地質(zhì)模型［J］. 地下水， 2015， 37（4）： 208-210.

［32］蒼正毅. 中基性火山巖儲(chǔ)層精細(xì)描述與地質(zhì)建模［D］. 長春： 吉林大學(xué)， 2023.

［33］劉友龍， 辛培業(yè)， 朱明， 等. 甘肅某金礦礦床工業(yè)指標(biāo)論證［J］. 現(xiàn)代礦業(yè)， 2023， 39（10）： 73-76.

［34］杜放. 貴州安龍戈塘金礦角礫狀礦石流體包裹體特征研究及意義［D］. 成都： 成都理工大學(xué)， 2017.

［35］CHEN J， HUANG Z L， YANG R D， et al. Gold and antimony metallogenic relations and ore-forming process of Qinglong Sb（Au） deposit in Youjiang basin， SW China： sulfide trace elements and sulfur isotopes［J］. Geoscience Frontiers， 2021， 12（2）： 605-623.

［36］GE X， SHEN C B， HE P， et al. The roles of hydrocarbons on the mineralization of Carlin-type gold deposits， Nanpanjiang Basin， South China［J］. Ore Geology Reviews， 2022， 149（1）： 105107.1-105107.15.

［37］CHEN J， YANG R D， DU L J， et al. Mineralogy， geochemistry and fluid inclusions of the Qinglong Sb-（Au） deposit， Youjiang basin （Guizhou， SW China）［J］. Ore Geology Reviews， 2018， 92（1）： 1-18.

［38］LI J， HU R， XIAO J， et al. Genesis of gold and antimony deposits in the Youjiang metallogenic province， SW China： evidence from in situ oxygen isotopic and trace element compositions of quartz［J］. Ore Geology Reviews， 2020， 116（1）： 103257.1-103257.16.

[39]LIU Y， HU K， HAN S， et al. Carlin-type gold deposit at Zumudang， SW China [J]. Resource Geology， 2015， 65（1）： 136-159.

［40］PANG B. The characteristic and origin of ore-forming fliud from micro-disseminated gold deposits in youjiang basin［J］. Geology and Prospecting， 2005， 41（1）： 13-17.

［41］MUNTEAN J L， CLINE J S. Diversity in Carlin-style gold deposits［M］. Colorado： Society of Economic Geologists， 2018.

［42］PENG Y， GU X， ZHANG Y， et al. Ore-forming process of the Huijiabao gold district， southwestern Guizhou Province， China： evidence from fluid inclusions and stable isotopes［J］. Journal of Asian Earth Sciences， 2014， 93（15）： 89-101.

［43］ZHENG L， TAN Q， ZUO Y， et al. Two hydrothermal events associated with Au mineralization in the Youjiang Basin， southwestern China［J］. Ore Geology Reviews， 2022， 144（1）： 104816.1-104816.18.

［44］XIE Z， HUANG K J， XIA Y， et al. Heavy δ26Mg values in carbonate indicate a magmatic-hydrothermal origin of Carlin-type Au deposit［J］. Geochimica et Cosmochimica Acta， 2022， 333（1）： 166-183.

［45］HU R G， PANG B C， BAI X J， et al. Progressive crushing 40Ar/39Ar dating of a gold-bearing quartz vein from the Liaotun Carlin-type gold deposit， Guangxi， southern China［J］. Scientific Reports， 2022， 12（1）： 12793-12803.

［46］琚宜文， 王桂梁， 胡超. 海孜煤礦構(gòu)造變形及其對(duì)煤厚變化的控制作用［J］. 中國礦業(yè)大學(xué)學(xué)報(bào)， 2002， 31（4）： 374-379.

［47］王海學(xué)， 王鳳蘭， 吳桐， 等. 斷裂垂向分段生長在油氣勘探開發(fā)中的應(yīng)用及發(fā)展趨勢(shì)［J］. 地球物理學(xué)進(jìn)展， 2019， 34（6）： 2336-2345.

［48］彭美麗， 易金， 姚蓉， 等. 華南多層推滑構(gòu)造系的物理模擬實(shí)驗(yàn)［J］. 大地構(gòu)造與成礦學(xué)， 2011， 35（2）： 198-206.

［49］劉舒飛，張繼林，鄂建新.貴州戈塘金礦床地質(zhì)特征及成因淺析［C］//中國有色金屬學(xué)會(huì)金屬礦山成礦理論與深部找礦新技術(shù)研討會(huì).深圳：中南大學(xué)出版社，2012： 44-54.

（責(zé)任編輯：于慧梅）

Study on Three-Dimensional Geological Modeling and Genesis of

Getang Gold Deposit in Southwestern Guizhou

CHAI Wenju TIAN Feng FU Yong ZOU Niuniu HUANG Yongzhi

（1.College of Resources and Environmental Engineering， Guizhou University， Guiyang 550025， China;

2.State Key Laboratory of Public Big Data， Guizhou University， Guiyang 550025， China;

3.Key Laboratory of Karst Georesources and Environment， Ministry of Education， Guizhou University， Guiyang 550025， China;

4.China Petrochemical Shengli Oilfield Petroleum Development Center Co.， Ltd.， Dongying 257099， China）Abstract： In order to intuitively show the spatial distribution law of Getang gold deposit and carry out comprehensive metallogenic analysis， based on the existing geological data such as boreholes and profiles， combined with the existing research results of ore-forming material sources， a three-dimensional geological model of Getang gold deposit was established using Petrel software. On this basis， the spatial distribution and enrichment characteristics of Getang gold deposit were discussed， the metallogenic process of the deposit was analyzed， and the metallogenic model of the deposit was summarized. The results show that Au is mainly enriched in the first member of Longtan Formation， and most of the ore bodies have low grade. The ore bodies with industrial grade are mostly concentrated near the fault; in general， the grade of ore bodies is negatively correlated with the volume of ore bodies， but the abnormal increase of volume occurs many times in high grade ore bodies. The distribution characteristics of the deposit along the strata are significant， which are obviously controlled by the unconformity surface or the interlayer slip surface， but the fault may provide an upward channel for the ore-bearing hydrothermal fluid. The formation of the deposit is closely related to magmatic activity; the hydrothermal fluid enters the ore-bearing layer along the fault and undergoes strong alteration， so that Au is precipitated and pre-enriched， and the gold deposit is formed under the action of groundwater in the later stage， which is the metallogenic model of volcanic sedimentation-hydrothermal transformation. The above research results have certain indication and reference significance for the genesis research and prospecting work of Carlin-type gold deposits in China.

Key words： Getang gold mine; three-dimensional geological modeling; carlin-type gold deposit; ore-forming analysis