A THREE DIMENSINAL ELASTOPLASTIC FINITE ELEMENT MODEL FOR PREDICTING DYNAMICALLY RESERVOIR EFFECTIVE STRESS
-
摘要: 本文根据油藏应力与渗流相互作用的耦合机理,建立了一个用于油藏有效应力动态预测的三维弹塑性模型。该模型强调了油藏在注采交变载荷作用下应力与渗流的力学耦合效应,以及油藏岩石弹塑性变形的力学特性,并采用弹塑性增量有限元法求解,最后给出了算例.Abstract: A three dimensinal elastoplastic model for predicting dynamically reservoir effective stress is established based on the coupled interaction mechanism of reservoir stress and multiphase fluid flow.The model is characteristic of two important aspects:the coupled effect of stress and seepage,and the elastoplastic properties of reservoir rocks under alternative loads of injection and production.The model is solved using elastoplastic incremental finite element method and an example is given.
-
Key words:
- effective stress /
- elastoplastic deformation /
- seepage /
- coupled effect /
- dynamical prediction /
- finite element
-
0. 引言
波龙斑岩铜金矿床位于西藏自治区改则县北西约100 km,班公湖—怒江成矿带西段。班公湖—怒江缝合带呈狭长带状近东西向展布,东西延伸超过2000 km,是一条横贯青藏高原中部的重要缝合带[1~3]。自2002年以来,西藏地质五队在发现地表具有相似的褐铁矿化的基础上,通过持续的地质工程施工控制,陆续发现了多不杂和波龙斑岩铜金矿床[4],而波龙矿床已成为青藏高原中部发现的最大的斑岩铜矿床。
自被确定为独立的斑岩铜矿以来,波龙矿床备受关注,目前开展了波龙矿床的成矿、成岩年龄[4~6]、矿床同位素地球化学[7~9]、矿床岩石学特征[10~11]、成矿构造环境[12]及成矿作用[13~14]等方面的研究。讨论波龙矿床流体方面的文献仅有2篇[15~16],其中李丹[15]虽然已对其进行了研究,但样品数量和统计数据有限,仅仅探讨了流体包裹体的基本特征和温度、盐度、密度、压力等条件;周玉[16]进行了稀有气体同位素研究,结合包裹体气相成分及H-O同位素分析,认为成矿流体属壳幔混源。本文从流体包裹体的气液相成分入手,结合前人的研究成果,进一步分析波龙矿床成矿流体的物理化学条件。
1. 矿床地质概况
波龙斑岩铜金矿床紧邻多不杂矿床,与拿若矿床(点)等组成了班公湖—怒江成矿带西段的多龙矿集区(见图 1)。
图 1 多龙矿集区构造背景简图(a)、地质简图(b)和波龙矿床地质简图(c)[4]Figure 1. Tectonic and geological maps of Duolong metallogenic district and geological map of Bolong porphyry Cu-Au deposit矿区出露的地层主要有下侏罗统曲色组第一岩性段(J1q1)、下白垩统美日切错组第一岩性段(K1m1)和第四系(Q)。区内地层呈单斜构造,未见褶皱和断裂。矿区岩浆岩以中酸性的花岗闪长斑岩为主,其次为脉岩,少量火山岩。矿区变质作用不发育,变质程度较浅,以区域变质作用为主,可见少量变石英砂岩和粉砂质板岩。
波龙斑岩铜矿床矿体以隐伏—半隐伏状产出,平面上呈似椭圆形状,地表未见露头。矿石构造以细脉—浸染状及浸染状构造为主,具结晶结构、交代结构和固溶体分离结构。矿石矿物主要为黄铁矿、黄铜矿,其次为磁黄铁矿、斑铜矿、辉锑矿、辉钼矿、磁铁矿等;脉石矿物以石英、长石、绢云母、硬石膏为主。矿区围岩与矿化蚀变发育,蚀变类型主要为硅化、黑云母化、绢云母化、黏土化等,其次为硬石膏化、碳酸盐化、绿帘石化和绿泥石化等[7, 12]。
根据本矿区野外观察、采样和室内磨片、镜下鉴定等研究,可将波龙铜矿床的形成过程划分为3个成矿期次:岩浆晚期、热液成矿期和表生期,其中热液成矿期又分为磁铁矿-辉钼矿阶段、黄铜矿-黄铁矿阶段和硬石膏-黄铁矿阶段。
2. 包裹体一般特征
波龙斑岩铜金矿床的流体包裹体分布主要呈成群、成带或星散状。形态多样,多为负晶形、椭圆形和不规则形,大小多在5~20 μm之间,最大可达70 μm。流体包裹体类型较复杂,主要为含石盐和硫化物子晶的气-液-固三相(L+V+H)包裹体,其次为气-液两相(L+V)包裹体,还有少量的纯气相(V)和纯液相(L)包裹体,其中气液比多在10%~50%之间[15]。
3. 包裹体温度、盐度、密度与压力特征
包裹体均一温度和冰点温度的测试在成都理工大学资源勘查工程系包裹体实验室进行,所用仪器为英国Linkam THNSG600型冷热台和吉林浑江市光学仪器厂TRL-02型热台。共测试了10件样品96个流体包裹体的均一温度。成矿深度和压力依据邵洁涟[17]提出的经验公式计算。依据测得的冰点温度,采用Hall等[18]的H2O-NaCl体系公式计算了盐度,根据刘斌[19]的密度计算公式得出了密度。测试及计算结果见表 1。
表 1 波龙铜金矿床流体包裹体均一温度、盐度、密度及压力测定结果Table 1. Homogenization temperature, salinity, density and pressure of fluid inclusions in Bolong Cu-Au deposit成矿阶段 主矿物 均一温度/℃ 盐度/% 密度/(g·cm-3) 压力/105Pa 变化范围 平均 变化范围 平均 变化范围 平均 变化范围 平均 岩浆晚期 斑晶石英 334~549 467 38.24~52.16 41.86 1.0683~1.1067 1.0761 58.76~234.41 102.47 磁铁矿-辉钼矿 石英 312~533 429 34.56~50.49 42.34 1.0683~1.0975 1.0771 26.72~212.69 109.59 黄铜矿-黄铁矿 石英 262~502 388 29.47~40.79 35.01 1.0683~1.1494 1.0930 3.14~86.81 33.54 硬石膏-黄铁矿 硬石膏石英 232~471 332 28.65~36.90 32.48 1.0765~1.1598 1.1153 1.40~45.74 17.68 由表 1可以看出,包裹体均一温度在232~549 ℃之间变化。随着成矿过程从岩浆晚期→磁铁矿-辉钼矿阶段→黄铜矿-黄铁矿阶段→硬石膏-黄铁矿阶段发展,均一温度逐渐降低(467→429→388→332 ℃),包裹体盐度(NaCl的质量分数,下同)也逐渐降低(41.86%→42.34%→35.01%→32.48%)。根据波龙矿区包裹体均一温度、盐度计算出的流体密度逐渐增高(1.0761→1.0771→1.0930→1.1153 g/cm3)。该矿床各成矿阶段的成矿压力都不高,平均值在17.68×105~109.59×105 Pa之间,基本呈逐渐降低的特点。矿床主要成矿于地下0.01~0.78 km深度,因此属浅成矿床。
包裹体测试数据显示,从岩浆晚期→磁铁矿-辉钼矿阶段→黄铜矿-黄铁矿阶段→硬石膏-黄铁矿阶段,波龙铜金矿床是在高温、低压、浅成—超浅成环境下形成的,成矿流体属于高盐度、中—高密度流体。
4. 逸度
根据波龙斑岩型铜金矿床脉石矿物石英中包裹体的气相成分与均一温度测试结果,按照有关热力学逸度计算公式[20],计算了不同阶段样品的氧逸度(lgfO2)、二氧化碳逸度(lgfCO2)和硫逸度(lgfS2)(见表 2)。由表 2可得,lgfO2均值的变化为-27.280→-29.631;lgfCO2均值的变化为1.851→1.242;lgfS2均值的变化为-3.284→-4.138,表明从磁铁矿-辉钼矿阶段→黄铜矿-黄铁矿阶段,流体逸度总体呈下降趋势。
表 2 波龙铜金矿床流体逸度Table 2. Fugacity of fluid inclusion in Bolong Cu-Au deposit, Tibet样号 lgfO2 lgfCO2 lgfS2 阶段 BL002 -27.280 1.851 -3.284 磁铁矿-辉钼矿 BL019 -28.778 1.295 -3.865 黄铜矿-黄铁矿 BL045 -30.167 1.177 -4.323 BL058 -29.946 1.254 -4.227 5. pH值和Eh值
根据包裹体气液相成分测试及均一温度,运用Crerar[21]和李葆华等[22]的公式,计算出波龙矿区成矿流体的酸碱度(pH)和氧化还原电位(Eh),结果见表 3。
表 3 波龙铜金矿床流体pH和Eh值Table 3. The values of pH and Eh of the fluid inclusions in Bolong Cu-Au deposit, Tibet样号 pH Eh 阶段 BL002 5.149 0.009 磁铁矿-辉钼矿 BL019 5.427 0.013 黄铜矿-黄铁矿 BL045 5.486 0.018 BL058 5.448 0.016 表 3显示,磁铁矿-辉钼矿阶段pH值为5.149,Eh值为0.009;黄铜矿-黄铁矿阶段pH平均值为5.454,Eh平均值为0.016,表明从磁铁矿-辉钼矿阶段到黄铜矿-黄铁矿阶段,流体pH和Eh值均具有增高的趋势。
6. 总硫活度和总碳活度
根据Helgeson[23]、Crerar[24]和Haymob的lgK值(某些气体元素和盐类平衡常数,转引自李葆华等[22])及已算出的pH、fO2、fS2和fCO2,得到硫和碳的各溶解类型的活度及总硫活度和总碳活度(见表 4)。由表 4可以看出,本矿床成矿流体中总硫活度aS∑为0.536~14.067 mol/L,总碳活度aC∑为0.660~3.121 mol/L。成矿溶液中硫的溶解类型以HSO4-和H2S形式为主,碳的溶解类型以H2CO3和CO2形式为主。从磁铁矿-辉钼矿阶段→黄铜矿-黄铁矿阶段,总硫活度和总碳活度依次降低,说明随着成矿过程的进行,硫化物和碳酸盐矿物的含量依次增加。
表 4 波龙铜金矿床总硫、总碳活度计算结果Table 4. Activity of total sulphur and total carbon in Bolong Cu-Au deposit, Tibetmol/L 样品编号 lgfS2 lgfO2 lgaH2S lgaHS- lgaS2- lgaHSO4- lgaSO42- lgaΣS aΣS lgfCO2 lgaCO2 lgaH2CO3 lgaHCO3- lgaCO32- lga∑C aΣC 阶段 BL002 -3.284 -26.714 -0.785 -4.636 -11.587 1.135 -0.596 1.148 14.067 1.851 0.151 0.231 -3.530 -10.811 0.494 3.121 磁铁矿-辉钼矿 BL019 -3.865 -28.152 -0.357 -3.930 -10.603 -1.033 -2.486 -0.271 0.536 1.295 -0.405 -0.325 -3.808 -10.811 -0.062 0.867 黄铜矿-黄铁矿 BL045 -4.323 -29.691 0.184 -3.330 -9.944 -3.512 -4.906 0.184 1.529 1.177 -0.523 -0.443 -3.867 -10.811 -0.180 0.660 BL058 -4.227 -29.513 0.143 -3.409 -10.062 -3.235 -4.667 0.143 1.391 1.254 -0.446 -0.366 -3.829 -10.811 -0.103 0.788 7. Cu的迁移、沉淀机制分析
热液在流经含矿岩系时活化并萃取了大量铜和其他成矿物质,随后演化成含矿热液。
铜在热液中的迁移形式主要是氯的络合物和硫氢络合物等,如Cu(HS)3-、CuS(HS)33-、CuCl3-等[25~26]。根据公式计算得到各种铜络合物离子活度和铜的总溶解度(见表 5)。
表 5 铜络离子活度及铜的总溶解度Table 5. Activity of Cu complex ion and total solubility样品编号 lgaCu(HS)2- lgaCu(H2S)(HS)2- lgaCu+ lgaCuCl0 lgaCuCl2- lgaCuCl32- lgaCu2+ lgaCuCl+ lgaCuCl20 lgaCuCl3- lgaCuCl42- lgaΣCu 阶段 BL002 -6.674 -7.148 -6.428 -14.607 -20.786 -30.765 -13.955 -18.834 -29.214 -40.093 -51.372 -6.183 磁铁矿-辉钼矿 BL019 -5.753 -5.800 -7.065 -21.923 -34.780 -51.438 -15.230 -26.788 -43.845 -61.403 -79.361 -5.464 黄铜矿-黄铁矿 BL045 -4.883 -4.389 -7.509 -12.438 -15.366 -22.095 -16.117 -17.746 -24.875 -32.504 -40.533 -4.268 BL058 -4.983 -4.530 -7.426 -13.964 -18.502 -26.841 -15.951 -19.190 -27.928 -37.167 -46.805 -4.398 由表 5可知:从磁铁矿-辉钼矿阶段到黄铜矿-黄铁矿阶段,成矿流体中的总铜活度变化为10-6.183~10-4.268 mol/L,显示出由低到高的变化规律。
根据计算出的Cu络离子活度值,可知Cu在热液中主要以Cu(HS)2-、Cu(H2S)(HS)2-、Cu+形式存在。随着流体pH值增大,温度降低,主要成矿元素开始以硫化物的形式沉淀下来。
8. 结论
波龙斑岩铜金矿床流体包裹体类型以含石盐和硫化物子矿物的气-液-固三相包裹体为主,其次为气-液两相包裹体,还有少量的纯气相及纯液相包裹体。
从流体测试分析数据可判断出,波龙斑岩铜金矿床形成于高温(232~549 ℃)、低压(1.40×105~234.41×105 Pa)和浅成—超浅成(0.01~0.78 km)环境;成矿流体表现出高盐度(28.65%~52.16%)和中—高密度(1.0683~1.1598 g/cm3)特征。从岩浆晚期→磁铁矿-辉钼矿阶段→黄铜矿-黄铁矿阶段→硬石膏-黄铁矿阶段,成矿流体的均一温度、成矿压力和盐度逐渐降低,而成矿流体密度逐渐升高。
随着成矿阶段的发展,各阶段流体氧逸度lgfO2均值(-26.714→-29.118)、二氧化碳逸度lgfCO2均值(1.851→1.242) 和硫逸度lgfS2均值(-3.284→-4.138) 逐步降低,pH值(5.149→5.454) 和Eh值(0.009→0.016) 逐步升高;总硫活度a∑S(101.148→100.019)和总碳活度a∑C(100.494→10-0.115)逐渐降低,与矿化具有一致性,说明硫化物和碳酸盐矿物的含量在逐步增加。
在成矿热液中,铜的络离子化合物主要以Cu(H2S)(HS)2-形式存在,其次是Cu(HS)2-,少量为CuCl2-和CuCl32-的形式存在,其余活度很低。
-
计量
- 文章访问数: 195
- HTML全文浏览量: 53
- PDF下载量: 8
- 被引次数: 0