趙義來，羅先熔，歐陽菲

(桂林理工大學 a.地球科學學院; b.廣西隱伏金屬礦產(chǎn)勘查重點實驗室,廣西 桂林 541004)

黑龍江金廠礦區(qū)熱液成礦過程計算模擬及成礦機制探討

趙義來，羅先熔，歐陽菲

(桂林理工大學 a.地球科學學院; b.廣西隱伏金屬礦產(chǎn)勘查重點實驗室,廣西 桂林 541004)

黑龍江金廠金礦是一特大型金礦床,該區(qū)中生代構(gòu)造-巖漿活動強烈,并伴隨大規(guī)模的熱液蝕變活動,形成了金廠礦床的隱爆角礫巖型、裂隙充填型、巖漿穹窿型及接觸帶型等多種礦化類型,而熱液擴容充填成礦可作為該礦區(qū)的代表性礦化類型。通過FLAC3D軟件進行金廠礦區(qū)熱液成礦過程計算模擬,結(jié)果表明:燕山晚期花崗斑巖對本區(qū)金成礦具有控制作用;本區(qū)溫度場的分布受成礦巖體規(guī)?？刂?受東西向壓應力場影響,本區(qū)出現(xiàn)了NW、NE向的剪切應變帶及NEE-SWW向展布的張應變擴容區(qū),擴容區(qū)不僅與流體匯流、實際成礦位置對應較好,還指示了新的找礦靶區(qū)。

熱液成礦;計算模擬;成礦機制;金廠金礦;黑龍江

0 引 言

礦床的形成是成礦系統(tǒng)多過程耦合演化的結(jié)果, 由于成礦系統(tǒng)演化的特殊性和復雜性, 研究一個區(qū)域的成礦動力學機制有著極大的困難, 表現(xiàn)在: 系統(tǒng)的空間尺度和時間尺度巨大, 其演化過程不可能在實驗室內(nèi)復制;成礦系統(tǒng)是一個復雜的多過程耦合系統(tǒng), 系統(tǒng)演化的結(jié)果和原因之間為非線性關(guān)系, 即系統(tǒng)初始條件的細微差異將導致系統(tǒng)演化結(jié)果的巨大差異[1-4]。因此, 在計算機技術(shù)和算法日益成熟的條件下, 計算模擬已成為研究上述復雜成礦系統(tǒng)的一種切實可行的方法。

巖漿熱液作用所成礦體的定位空間必須滿足兩個方面的動力學條件:一是這個局部空間是擴容空間,因為礦體所定位的空間原本為固液相互平衡的多孔巖石介質(zhì),如果不增加空隙的體積,運載成礦物質(zhì)的流體就不可能向這個空間匯流,礦石也難以在這個空間發(fā)生大規(guī)模堆積;二是有序的動力條件,如果沒有有序的動力作用,不可能長時間地維持流體固定的局部匯流。

黑龍江省東寧縣金廠金礦是一個特大型金礦床,位于中亞造山帶東端的吉黑成礦帶,是與中生代中酸性巖漿巖活動有關(guān)的復雜成礦系統(tǒng)[5-6]。在金廠礦區(qū)內(nèi),成礦的主要動力來源有兩個:區(qū)域應力場和巖體的熱動力條件,在上述兩個因素作用下的熱液流體流動則構(gòu)成礦質(zhì)輸運、堆積的媒介條件。由于熱液成礦作用主要發(fā)生在巖體侵位、基本固化以后,此時整個巖體-圍巖體系可以看成粘彈性的多孔介質(zhì),其力學行為可用Mohr-Colomb模型描述。因此,本文選用Itasca公司的FLAC3D(Fast Lagrangian Analysis of Continua in 3 Dimensions)軟件[7]進行金廠礦區(qū)巖體基本固化后的同構(gòu)造冷卻過程中變形-傳熱-流體流動耦合動力學的計算模擬,探討力、熱、流因素影響下的礦體定位機制。

1 礦區(qū)地質(zhì)特征

黑龍江金廠金礦是一特大型金礦床, 該礦區(qū)大地構(gòu)造位置處于太平嶺隆起與老黑山斷陷的過渡部位, 中生代構(gòu)造-巖漿活動十分強烈, 形成了一

礦區(qū)內(nèi)的金礦化類型主要有隱爆角礫巖型、裂隙充填型、巖漿穹窿型及接觸帶型4類。巖漿穹窿型礦體產(chǎn)于巖體的破碎裂隙中,因受裂隙形態(tài)控制,礦體產(chǎn)出為環(huán)狀及薄層狀;對于接觸帶型礦體,其礦石產(chǎn)出于接觸帶附近的裂隙中,且具有明顯的破碎角礫特征和類似于其他礦體的蝕變特征,并不具有典型矽卡巖礦體那樣的接觸交代及蝕變類型。因此,熱液擴容充填成礦可作為該礦區(qū)的代表性礦化類型。礦區(qū)賦礦圍巖主要有斜長花崗巖、花崗閃長巖、花崗斑巖、閃長巖、閃長玢巖等。圍巖蝕變比較發(fā)育,主要類型有硅化、絹云母化、高嶺土化、黃鐵礦化;其次有鉀長石化、碳酸鹽化、綠泥石化等,其中硅化、黃鐵礦化與金礦化關(guān)系密切。

圖1 金廠礦區(qū)地質(zhì)圖(據(jù)文獻[10]修改)

2 金廠礦區(qū)成礦動力學建模

2.1 建模過程及方法

從原始圖件到計算用網(wǎng)格化模型的建立所使用的軟件包括MapGIS、AutoCAD和ANSYS,建模工作流程及軟件操作分述如下。

(1)原始圖件為MapGIS工程文件, 系從相關(guān)研究人員及單位處收集而來的地質(zhì)圖, 由工程文件輸出為可被AutoCAD使用的*.dxf文件。

(2)在AutoCAD軟件平臺中進行的操作包括: 圖件清繪, 有用圖元的提取, 無關(guān)信息的刪減, 修改線型, 拓撲錯誤檢查與處理等; 然后通過對二維圖件的三維拉伸、旋轉(zhuǎn)及布爾運算等操作, 建立研究區(qū)的實體模型, 導出*.sat文件。

(3)ANSYS軟件為一款通用有限元分析軟件, 適用于進行結(jié)構(gòu)、流體、電場、 磁場、聲場分析, 本文使用的是該軟件的網(wǎng)格剖分功能。 將*.sat文件導入到ANSYS, 在ANSYS中對實體模型進行網(wǎng)格剖分, 以得到可用于動力學計算的網(wǎng)格模型。 剖分過程一方面要進行模型基本單元的選擇, 本文采用四面體單元對實體模型進行剖分, 該類單元對不規(guī)則的地質(zhì)體適應性強, 發(fā)生網(wǎng)格報錯的幾率也較低。 另一方面,進行模型基本單元尺寸的設(shè)置, 若尺寸設(shè)置較大, 則剖分得到的網(wǎng)格數(shù)量少, 模型粗略, 但運算速度快; 若尺寸設(shè)置較小, 剖分得到的網(wǎng)格數(shù)量多, 則模型精準, 但運算速度較慢。 一般的做法是權(quán)衡精確性和運行速度, 并結(jié)合建模經(jīng)驗進行設(shè)置;或?qū)δＰ瓦M行預運算后, 根據(jù)運算情況調(diào)整單元尺寸參數(shù), 重新進行模型剖分。 此外, 在局部地質(zhì)界線復雜部位進行網(wǎng)格加密, 以獲得較大的精確性。 在剖分過程中設(shè)置巖石單元分組信息, 以便于后續(xù)工作對不同的巖石單元賦參數(shù)及設(shè)置初始動力學條件。

(4)由ANSYS導出的網(wǎng)格化模型信息記錄于文本文檔,其格式與動力學模擬軟件FLAC3D所接受的文件格式有所不同, 通過人機交互或C++ 編程對文件格式進行轉(zhuǎn)換,即可用于FLAC3D進行模擬計算。

根據(jù)對已有資料及研究成果的判讀分析[10-17],建立了如圖2所示的FLAC3D網(wǎng)格化模型,模型大小為LEW×LSN×H=11 km×9 km×0.6 km(LEW、LSN及H分別表示模型東西、 南北向長度

圖2 金廠礦區(qū)四面體網(wǎng)格化三維模型及邊界條件

及厚度), 共剖分出29 965個四面體單元。 建模過程中,對原始圖件地質(zhì)信息的處理包括:

①在金廠礦區(qū)地表地質(zhì)圖的基礎(chǔ)上, 剝除新近系、 第四系沉積層, 將其所覆蓋的其余巖石單元的邊界進行推延, 并兼顧不同期次巖體的穿插關(guān)系；

②考慮到對模型網(wǎng)格數(shù)量的控制, 刪除大片巖體內(nèi)部的零星小侵入體, 進行此操作的前提是小侵入體與其所在的大規(guī)模巖體具有類似的巖性, 二者為巖性過渡或局部相變關(guān)系；

③將巖石單元狹長的尖滅部位進行圓滑處理, 以避免網(wǎng)格剖分過程中產(chǎn)生畸變的四面體單元, 導致在計算過程中出現(xiàn)負體積等相悖于客觀事實的情況；

④因該區(qū)大部分范圍的深部暫無鉆孔揭露, 僅基于小范圍的鉆探工程對整個區(qū)域深部巖石單元界線較難作出合理的推測, 且相對于該區(qū)較大的長寬尺度而言,模型厚度方向上的界線變化對成礦過程的影響相對較弱, 故暫不考慮該方向的界面處理。參照現(xiàn)有見礦深度, 將模型地表地質(zhì)界線作垂向順延, 取模型厚度為600m。

2.2 模型參數(shù)及條件設(shè)置

模型各巖石單元的力、熱、流參數(shù)(表1)主要根據(jù)各單元的巖石成分,在Sch?n(1998) 主編的《巖石物理性質(zhì)》[18]以及《FLAC3D應用手冊》[7]提供的數(shù)據(jù)范圍內(nèi)通過對比不同參數(shù)的模擬實驗結(jié)果而選擇。所要模擬的并不是燕山晚期巖體上侵到成礦演化的整個過程,而是巖體主體基本固結(jié)后熱液成礦的物理過程。在給定了初始和邊界條件后, 通過區(qū)域變形量的控制計算, 根據(jù)輸出力、 熱、 流等因素的計算結(jié)果進行分析。

表1 金廠模型各巖石單元物理參數(shù)

3 計算結(jié)果及成礦機制分析

3.1 計算結(jié)果

控制計算運行至該區(qū)域東西向縮短量達到10%,輸出溫度、剪切應變、體積應變、流體壓力及流體流動的結(jié)果(圖3～圖6)。這些結(jié)果表明，在熱、構(gòu)造應力和流體的三重作用下,整個區(qū)域的溫度場、應變場及流場的分布具有不均勻性,這種耦合作用下的多場不均勻分布對該區(qū)金礦成礦起主導作用。

3.2 成礦機制分析

根據(jù)上述計算結(jié)果,不但可以解釋已發(fā)現(xiàn)礦體的成礦規(guī)律,還對認識金廠礦區(qū)金礦床的成礦機制以及燕山晚期巖體在成礦中的作用提供了有力的證據(jù)。

(1)分別設(shè)置燕山晚期的花崗斑巖體及閃長質(zhì)巖體為主成礦巖體,通過對比計算發(fā)現(xiàn),花崗斑巖對本區(qū)金的成礦具有主要的控制作用,本文僅給出了最優(yōu)計算結(jié)果,未考查閃長質(zhì)巖體成礦作用的結(jié)果。

(2)本區(qū)溫度場的分布明顯受成礦巖體的控制(圖3), 經(jīng)過10Ma的熱擴散作用, 可以看到,高溫區(qū)仍環(huán)繞花崗斑巖體分布,且受巖體規(guī)模大小控制,如本區(qū)南西部位的一處花崗斑巖中心溫度為550 ℃,而區(qū)域中部呈NW向排列的幾處小巖體的溫度已降至約350 ℃。

圖3 金廠礦區(qū)溫度分布

圖4 流體流動與剪切應變關(guān)系

圖5 流體流動與體積應變關(guān)系

圖6 流體流動與孔隙流體壓力分布關(guān)系

(3)受東西向擠壓應力場控制,本區(qū)出現(xiàn)了3條大的剪切應變帶,即北西部剪應變帶、東部剪應變帶及南部剪應變帶,均遭受了較大的剪切作用(圖4)。3條剪應變帶的分布方向為NW向及NE向,是在東西向壓應力場作用下派生出的共軛剪切應變帶。剪應變的分布與規(guī)模較大的花崗斑巖體有對應關(guān)系,與小規(guī)?；◢彴邘r及其他巖石單元的對應性不強,表明剪應變的出現(xiàn)不僅受區(qū)域應力場控制,還可能受熱擴散作用影響,散熱慢的巖體對應較強的剪應變區(qū)。

流體流動的匯集區(qū)受剪切作用的控制不明顯,流體匯流部位并未出現(xiàn)在大的剪應變區(qū),說明該區(qū)剪切裂隙不是主要的控礦、容礦構(gòu)造。

(4)該區(qū)出現(xiàn)了NEE-SWW向展布的張性應變區(qū)(圖5),即正的體積應變區(qū),其具體分布為:西部花崗斑巖體的南側(cè)及北側(cè)接觸帶部位、中部巖體的北側(cè)接觸帶與金廠村南、南東部巖體與圍巖的接觸帶上。上述張應變區(qū)的展布方向與區(qū)域擠壓方向一致,是東西向壓應力場作用下派生的近東西向擴容區(qū),即等效于南北向拉伸作用的產(chǎn)物。由于擴容區(qū)均與巖體接觸帶有關(guān),推測溫度差異導致的巖體與圍巖力學性質(zhì)差異對其有一定影響。由計算結(jié)果可以看出,該區(qū)的壓應變帶主要分布在礦區(qū)南部、金廠村西,其次為礦區(qū)北東部,結(jié)合前述剪應變結(jié)果,這些區(qū)域應主要受壓扭性應力場控制。

流體匯流與擴容區(qū)的對應關(guān)系較好,主要的流體匯流部位分布在金廠村南、八號硐及其南西部、窮棒子溝中段東西兩側(cè)、礦區(qū)中部閃長玢巖南側(cè)及礦區(qū)南東角部位。這些區(qū)域不僅對應了該區(qū)NEE-SWW向的角礫巖型、裂隙充填型礦帶,還指示了金廠村南、八號硐南西部及礦區(qū)南東角部位為有利成礦部位及找礦優(yōu)選區(qū)。

上述流體匯集區(qū)的產(chǎn)生，一方面因為擴容區(qū)的存在代表張性裂隙發(fā)育的有利地段,而張性裂隙的產(chǎn)生則成為礦液運移的有利通道及匯集的場所;另一方面,孔隙流體壓力的分布也有一定貢獻(圖6),即模擬結(jié)果中流體均由高孔壓區(qū)流向低孔壓區(qū),而礦區(qū)東西兩側(cè)的具有最大流體壓力的部位均少有流體出現(xiàn),也是這一原因。另外,局部擴容影響孔壓分布,進而影響流體匯流的一般規(guī)律也為上述結(jié)果提供了佐證。

4 結(jié) 論

燕山晚期花崗斑巖對本區(qū)金的成礦具有主要的控制作用;溫度場的分布主要受巖體規(guī)模大小影響;在東西向擠壓應力場的作用下,派生出NW向及NE向的共軛剪應變帶；另外,受熱擴散作用影響,剪應變的分布與規(guī)模較大的花崗斑巖體有對應關(guān)系,流體流動的匯集區(qū)受剪切作用的控制不明顯;東西向壓應力場還派生出近東西向的擴容區(qū),流體匯流與擴容區(qū)的對應關(guān)系較好。

綜上,受控于區(qū)域東西向擠壓應力場、巖體與圍巖的溫度差異及熱力學性質(zhì)差異,巖漿期后熱液運移特征不僅揭示了該區(qū)成礦機制,印證了已發(fā)現(xiàn)礦體的產(chǎn)出規(guī)律,還指示了進一步工作的勘查方向,而該區(qū)金礦體的形成是金廠燕山晚期巖體成礦系統(tǒng)變形-熱擴散-流體運移多過程耦合作用的結(jié)果。

本文的計算模擬工作得到了中南大學計算地學研究中心的支持,在此謹表謝意!

[1]LiuLM,ZhaoYL,SunT. 3Dcomputationalshape-andcoolingprocess-modelingofmagmaticintrusionanditsimplicationforgenesisandexplorationofintrusion-relatedoredeposits:AnexamplefromtheYueshanintrusioninAnqing,China[J].Tectonophysics, 2012,526-529:110-123.

[2]LiuLM,ZhaoYL,ZhaoCB.CoupledgeodynamicsintheformationofCuskarndepositsintheTongling-Anqingdistrict,China:Computationalmodelingandimplicationsforexploration[J].JournalofGeochemicalExploration, 2010, 106(1-3):146-155.

[3]ZhaoCB,HobbsBE,MühlhausHB,etal.Computersimulationsofcoupledproblemsingeologicalandgeochemicalsystems[J].ComputerMethodsinAppliedMechanicsandEngineering, 2002,191(29-30):3137-3152.

[4]LiuLM,WanCL,ZhaoCB,etal.GeodynamicconstraintsonorebodylocalizationintheAnqingorefield,China:Computationalmodelingandfacilitatingpredictiveexplorationofdeepdeposits[J].OreGeologyReviews, 2011, 43(1):249-263.

[5]張宇, 賴勇, 卿敏, 等. 黑龍江省金廠金礦床J0 礦體流體地球化學研究[J]. 巖石學報, 2008, 24(5):1131-1144.

[6]朱成偉, 陳錦榮, 李體剛, 等. 黑龍江金廠金礦床地質(zhì)特征及成因探討[J]. 礦床地質(zhì), 2003, 22(1):56-64.

[7]GroupIC.FLAC3DUser’smanual[K].Minneapolis:ItascaConsultingGroup,Inc., 2005.

[8]賈國志, 陳錦榮, 楊兆光, 等. 金廠特大型金礦床的地質(zhì)特征與成因研究[J]. 地質(zhì)學報, 2005, 79(5): 661-670.

[9]趙玉鎖, 閆家盼, 張桂鳳, 等. 黑龍江金廠銅金礦床地質(zhì)特征及成因探討[J]. 地質(zhì)與勘探, 2013, 49(1):75-88.

[10]王艷忠, 楊言辰, 陳桂虎, 等. 黑龍江東寧金廠金礦床角礫巖型金礦體地質(zhì)特征及找礦標志[J]. 黃金, 2012, 33(11):24-28.

[11]王振忠. 黑龍江東寧金廠金礦床地質(zhì)特征及找礦前景分析[J]. 吉林地質(zhì), 2008, 27(2):31-36.

[12]趙玉鎖. 延邊-東寧成礦帶金廠斑巖金成礦系統(tǒng)[D]. 北京: 中國地質(zhì)大學(北京), 2013.

[13]李萬賀. 黑龍江東寧縣金廠金礦區(qū)成礦規(guī)律和成礦預測研究[D]. 北京: 中國地質(zhì)大學(北京), 2010.

[14]王可勇, 卿敏, 張欣娜, 等. 黑龍江金廠金礦床流體包裹體特征及成礦作用研究[J]. 巖石學報, 2011, 27(5):1275-1286.

[15]肖力, 趙玉鎖, 孔媛媛, 等. 東寧金廠礦區(qū)及外圍構(gòu)造體系厘定及其控巖控礦作用[J]. 礦床地質(zhì), 2010(S2): 63-70.

[16]陳彪，羅先熔, 邱煒, 等.[J].地電提取法在黑龍江金廠金礦區(qū)尋找隱伏金礦的應用[J]. 桂林工學院學報, 2009, 29(1):31-34.

[17]李兆誼, 羅先熔, 候冬梅, 等. 地電化學集成技術(shù)在黑龍江金廠金礦的應用研究及找礦預測[J]. 黃金科學技術(shù), 2010, 18(5):117-122.

[18]Sch?nJH.PhysicalPropertiesofRocks:FundamentalsandPrinciplesofPetrophysics[M].Oxford:Elsevier, 1998.

Numerical simulation of hydrothermal mineralization and metallogenic mechanism for Jinchang deposit in Heilongjiang

ZHAO Yi-lai, LUO Xian-rong, OUYANG Fei

(a.College of Earth Sciences; b.Guangxi Key Laboratory of Hidden Metallic Ore Deposits Exploration, Guilin University of Technology, Guilin 541004,China)

The Jinchang gold deposit is a super-large deposit with frequent Mesozoic tectonic-magmatic activities and extensive hydrothermal alteration. The orebodies are of cryptoexplosion breccia type, crack-filling type, magma dome type and contact zone type. The crack-filling type mineralization in hydrothermal dilation spaces is dominant. Computational modeling of the mineralization process is carried out on FLAC3D. The modeling results demonstrate that Late Yanshanian granite-porphyry plays an important role in the gold mineralization, and the distribution of temperature is controlled by the scale of granite-porphyry. NW and NE trending shear strain zones and NEE-SWW trending tensile strain areas are derived from the regional EW compression. The tensile dilation areas are consistent with both fluid focusing and mineralization positions. The results indicate new targets for exploration.

hydrothermal mineralization; computational modeling; metallogenic mechanism; Jinchang gold deposit;Heilongjiang

1674-9057(2015)04-0780-06

10.3969/j.issn.1674-9057.2015.04.016

2015-05-27

廣西自然科學基金項目(2014GXNSFBA118231);武警黃金指揮部項目(01-03);廣西礦冶與環(huán)境科學實驗中心項目

趙義來(1982—),男,博士,講師,地質(zhì)資源與地質(zhì)工程專業(yè),zyl@glut.edu.cn。

趙義來，羅先熔，歐陽菲.黑龍江金廠礦區(qū)熱液成礦過程計算模擬及成礦機制探討[J].桂林理工大學學報，2015,35(4):780-785.

P618.51

A