王自國(guó)，白維燦

(1.中國(guó)煤炭地質(zhì)總局勘查研究總院，北京 100039； 2.中煤地質(zhì)工程總公司，北京 100040)

構(gòu)造變形在地殼中的分布具有一定的規(guī)律，這種規(guī)律性構(gòu)造應(yīng)力場(chǎng)受巖石力學(xué)性質(zhì)控制。構(gòu)造應(yīng)力場(chǎng)的有限元法模擬直接建立在地質(zhì)實(shí)際基礎(chǔ)之上，是實(shí)現(xiàn)構(gòu)造應(yīng)力場(chǎng)從定性描述到定量半定量計(jì)算的有效方法。在研究區(qū)域控礦構(gòu)造特征，特別是成礦期構(gòu)造特征的基礎(chǔ)上，建立了研究區(qū)地質(zhì)實(shí)體的離散化模型，模擬出研究區(qū)成礦期構(gòu)造應(yīng)力場(chǎng)，從而為成礦預(yù)測(cè)提供有價(jià)值的信息。

前孫家金礦和金華山金礦隸屬于山東省招遠(yuǎn)市蠶莊鎮(zhèn)，位于前孫家—洼孫家斷裂帶的南段，本文以這兩個(gè)金礦為例，對(duì)上述兩金礦區(qū)進(jìn)行了成礦期構(gòu)造應(yīng)力場(chǎng)數(shù)值模擬。

1 區(qū)域應(yīng)力場(chǎng)與研究區(qū)成礦期應(yīng)力場(chǎng)

對(duì)于膠東地區(qū)的構(gòu)造變形研究表明，跟金礦化有關(guān)的各種構(gòu)造形跡主要是在燕山期形成的。膠東西北部燕山期形成的主要控礦斷裂帶，早期在NW-SE向構(gòu)造應(yīng)力擠壓下，均表現(xiàn)為韌性擠壓；而燕山中晚期，區(qū)域最大主壓應(yīng)力轉(zhuǎn)為NE-SW向，在構(gòu)造抬升、巖漿冷卻以及NE向構(gòu)造擠壓作用下，三山島斷裂、招平斷裂和焦家斷裂由韌性變形轉(zhuǎn)為以張性、脆性變形為主，為后期含金熱液上升提供了通道。NE向最大主應(yīng)力作用下形成的NEE、NNE向的張性斷裂是膠東西北部主要的容礦構(gòu)造。

在大地構(gòu)造背景為NE-SW向縮短、NW-SE向伸展?fàn)顟B(tài)下，膠東西北部金礦成礦時(shí)期構(gòu)造應(yīng)力場(chǎng)以NW-SE向伸展為主。前孫家金礦床及金華山金礦床是在膠東區(qū)域成礦大背景下形成的，受焦家斷裂與靈北斷裂差異運(yùn)動(dòng)的影響而形成近SN向縮短，近EW向伸展的特殊應(yīng)力場(chǎng)。

2 礦區(qū)礦體分布

2.1 前孫家礦區(qū)

前孫家礦區(qū)主要有2號(hào)脈、4號(hào)脈、9號(hào)脈、新1號(hào)脈(圖1)?，F(xiàn)主要簡(jiǎn)述2號(hào)脈、4號(hào)脈及新1號(hào)脈特征。

圖1 前孫家礦區(qū)三中段礦體分布圖Figure 1 Orebody distribution in third middle sector, Qiansunjia gold mine

2號(hào)脈賦存在前孫家—洼孫家主破碎帶上，礦脈走向約NE30°，傾向SE，傾角65°～75°。4號(hào)脈主要分布在8-36線，總體走向NE35°左右，傾向NW，傾角55°～72°。新1號(hào)脈是在二中段以下發(fā)現(xiàn)的一個(gè)盲礦體。它賦存于2號(hào)脈和4號(hào)脈之間，走向近SN，傾向西，傾角較陡，在72°～78°。

2.2 金華山礦區(qū)

金華山礦區(qū)主要1號(hào)脈、2號(hào)脈、3號(hào)脈、4號(hào)脈(圖2)。1號(hào)脈走向北東10°?？傮w傾向北西西，局部反傾，傾角72°～78°。2號(hào)脈走向南北，傾向西，傾角70°～75°，呈透鏡狀，向南西側(cè)伏。3號(hào)脈走向北北西，傾向南西西，傾角70°左右，呈透鏡狀，向南西側(cè)伏。礦體走向南北，傾向西，傾角75°～80°。礦體呈長(zhǎng)透鏡狀。

圖2 金華山礦區(qū)七中段礦體分布圖Figure 2 Orebody distributions in seventh middle sector,Jinhuashan gold mine

3 礦區(qū)構(gòu)造應(yīng)力場(chǎng)的有限元數(shù)值模擬

3.1 邊界條件的確定

確定邊界條件一般要考慮研究對(duì)象的范圍和邊界四周的地質(zhì)情況，通常取斷層面、巖體的邊界面等天然面作為邊界，也可任取一面為邊界。前孫家金礦東西以2號(hào)脈斷裂和4號(hào)脈斷裂為界面，南北選一任意界面為邊界。使南、西兩邊固定，在北、東兩邊的邊界上相對(duì)加力，為使模型不發(fā)生剛體位移，在西南角加一固定約束。金華山金礦南、北以前孫家—洼孫家斷裂和其下盤與其平行的斷裂為界面，東、西選一任意界面為邊界。使南、西兩邊固定，在北、東兩邊的邊界上相對(duì)加力，為使模型不發(fā)生剛體位移，在西南角加一固定約束。根據(jù)鄧軍等在該區(qū)所測(cè)定的成礦期最大主應(yīng)力方向，前孫家金礦區(qū)為325°，金華山金礦區(qū)為5°。作用力的大小選一相對(duì)應(yīng)力值。

3.2 巖石力學(xué)參數(shù)的設(shè)定

研究區(qū)內(nèi)主要的巖石類型為玲瓏花崗巖和郭家?guī)X花崗閃長(zhǎng)巖及斷裂構(gòu)造巖，本文所需的巖石力學(xué)參數(shù)彈性模量(E)、泊松比(μ)引自鄧軍，巖石密度(ρ)引自山東省巖土公司(表1)。

表1 巖石力學(xué)參數(shù)表

3.3 有限單元的劃分及模擬

邊界線要盡量與構(gòu)造線、巖體與圍巖的接觸線、不同力學(xué)參數(shù)的巖層接觸線重合，重點(diǎn)地段加密，次要地段適當(dāng)放大。三角形單元盡量使其邊長(zhǎng)近等，避免出現(xiàn)鈍角 。若邊界線是曲線就用幾條直線代替曲線作為單元的邊。前孫家礦區(qū)共劃分為396個(gè)四邊形和三角形，412個(gè)結(jié)點(diǎn)(圖3)。金華山共劃分為400個(gè)四邊形和三角形，411個(gè)結(jié)點(diǎn)(圖4)。

圖3 前孫家礦區(qū)構(gòu)造應(yīng)力場(chǎng)有限元網(wǎng)格劃分圖Figure 3 Tectonic stress field finite element mesh partition in Qiansunjia gold mine

圖4 金華山礦區(qū)構(gòu)造應(yīng)力場(chǎng)有限元網(wǎng)格劃分圖Figure 4 Tectonic stress field finite element mesh partition in Jinhuashan gold mine

利用構(gòu)造應(yīng)力場(chǎng)有限元數(shù)值模擬軟件對(duì)成礦期構(gòu)造應(yīng)力場(chǎng)進(jìn)行有限元數(shù)值模擬，得出成礦期最大張應(yīng)力分布圖(圖5、圖6)。

圖5 前孫家金礦區(qū)成礦期最大張應(yīng)力分布圖Figure 5 Metallogenic epoch maximum tensile stress distribution map of Qiansunjia gold mine

圖6 金華山金礦區(qū)成礦期最大張應(yīng)力分布圖Figure 6 Metallogenic epoch maximum tensile stress distribution map of Jinhuashan gold mine

4 模擬結(jié)果的分析

決定應(yīng)力狀態(tài)的主要因素是模型的邊界條件。因此，探討和設(shè)置真實(shí)的或接近真實(shí)的邊界條件是模擬實(shí)驗(yàn)的關(guān)鍵。有時(shí)邊界條件設(shè)置后經(jīng)過(guò)運(yùn)算，還要通過(guò)反演進(jìn)行調(diào)試，使邊界條件更接近真實(shí)，使模擬的結(jié)果與自然界的實(shí)際情況相一致。

從前孫家金礦區(qū)成礦期最大張應(yīng)力分布圖(圖5)可以看出(主要研究2號(hào)脈斷裂與4號(hào)脈斷裂所夾持部分)，在2號(hào)脈斷裂與4號(hào)脈斷裂之間，應(yīng)力集中點(diǎn)一部分基本落在新1號(hào)的脈斷裂的位置，這與實(shí)際已知的新1號(hào)脈礦體賦存位置基本吻合。另一些應(yīng)力集中點(diǎn)所落的位置，也連成一條北北東走向的未知斷裂(推測(cè)斷裂)，在此位置，也可能賦存有未知的含礦斷裂。

從金華山礦區(qū)成礦期最大張應(yīng)力分布圖(圖6)可以看出(主要研究前孫家—洼孫家斷裂及下盤平行斷裂所夾持部分)，在前孫家—洼孫家斷裂與下盤平行斷裂之間，應(yīng)力集中點(diǎn)大部分落在實(shí)際已知的斷裂處，與礦區(qū)的1號(hào)脈斷裂、2號(hào)脈斷裂、3號(hào)脈斷裂、4號(hào)脈斷裂基本相對(duì)應(yīng)，與礦體賦存位置基本吻合。如1號(hào)脈斷裂的礦體較大，所對(duì)應(yīng)的應(yīng)力集中區(qū)也較大，2、3、4號(hào)脈斷裂的礦體較小，所對(duì)應(yīng)的應(yīng)力集中區(qū)也較小。另外，一少部分應(yīng)力集中點(diǎn)所落的位置， 根據(jù)礦區(qū)已知含礦斷裂的走向?yàn)楸北睎|或南北向，在這里也連成一條北北東向未知斷裂(推測(cè)斷裂1)和一條近南北向的未知斷裂(推測(cè)斷裂2)，推測(cè)在此位置也可能賦存有未知的含礦斷裂。

5 結(jié)論

由此來(lái)看，成礦期構(gòu)造應(yīng)力場(chǎng)有限元模擬的結(jié)果不僅與已知的含礦斷裂基本相吻合，而且還模擬出一些未知的斷裂，通過(guò)模擬的這些未知斷裂，可能找到新的礦體。

根據(jù)成礦期構(gòu)造應(yīng)力場(chǎng)有限元模擬的信息，結(jié)合兩礦區(qū)的構(gòu)造控礦規(guī)律，金華山礦區(qū)在4號(hào)脈的東部施工了探礦工程，在推測(cè)斷裂1的位置發(fā)現(xiàn)了一條近南北向的礦脈，在前孫家礦區(qū)模擬的推測(cè)斷裂處也發(fā)現(xiàn)了礦化蝕變帶，從而驗(yàn)證了成礦期構(gòu)造應(yīng)力場(chǎng)有限元數(shù)值模擬提供了有價(jià)值的找礦信息。

[1]張業(yè)明，孫德育，常文志.熱液金礦床成礦構(gòu)造應(yīng)力場(chǎng)有限元模擬[J].長(zhǎng)春地質(zhì)學(xué)院學(xué)報(bào)，1996，26(2)：179-190.

[2]馬德云，高振敏，楊世瑜.北衙金礦床構(gòu)造應(yīng)力場(chǎng)數(shù)值模擬[J].大地構(gòu)造與成礦學(xué)，2003，27(2)：160-166.

[3]萬(wàn)天豐.中國(guó)東部中、新生代板內(nèi)變形構(gòu)造應(yīng)力場(chǎng)及其應(yīng)用[M].北京：地震出版社.1993.

[4]朱光，王道軒，劉國(guó)生.郯廬斷裂帶的伸展活動(dòng)及其動(dòng)力學(xué)背景[J].地質(zhì)科學(xué)，2001,36(3)：269-278.

[5]鄧軍，楊立強(qiáng)，方云，等.膠東地區(qū)殼-幔作用與金礦效應(yīng)[J].地質(zhì)科學(xué)，2000,35(1)：60-70.

[6]姜喜榮，趙寅震，肖勁民.內(nèi)蒙古蓮花山金礦區(qū)成礦應(yīng)力場(chǎng)數(shù)值模擬與隱伏與隱伏礦床預(yù)測(cè)[J].地質(zhì)力學(xué)學(xué)報(bào)，1995，1(1)：82-87.

[7]王自國(guó)，錢建平，王秀娟.金華山金礦區(qū)構(gòu)造控礦規(guī)律及找礦預(yù)測(cè)[J].有色金屬(礦山部分)，2012，64(2)：58-62.

[8]鄧軍，方云，周顯強(qiáng)，等.山東膠西北金礦帶成礦構(gòu)造應(yīng)力場(chǎng)反演及其控礦作用[J]，中國(guó)區(qū)域地質(zhì)，1995(3)：252-260.