Hydrocarbon generation potential and gas genesis of Cretaceous source rocks in Hesigewula Sag, western margin of Great Khingan
-
摘要:
大兴安岭西缘贺斯格乌拉凹陷为新发现凹陷，其潜在烃源岩的生烃潜力及天然气成因尚未有过系统研究。对该凹陷3口参数井及部分煤田探井白垩系160块岩心样品进行有机岩石学、有机地球化学等分析，结合构造、沉积特征，系统分析了该凹陷白垩系泥岩生烃潜力，探讨了贺地1井天然气的成因类型。贺斯格乌拉凹陷下白垩统暗色泥岩有机质丰度较高，腾格尔组一段一般泥岩TOC平均为2.20%，腾格尔组二段TOC平均为3.67%；以Ⅱ1—Ⅱ2型有机质为主；成熟度较低，目前处于未熟-低熟阶段。贺地1井天然气甲烷碳同位素较轻，含量为-70.5‰~-68.2‰，结合甲烷氢同位素特征，认为该井天然气为生物成因气，二氧化碳还原与乙酸发酵混合生成。因所采岩心样品普遍浅于1000 m，推测凹陷深部（大于2000 m）已进入到生油高峰期，具有一定的油气勘探潜力。
Abstract:The Hesigewula Sag on the western margin of the Great Khingan is a newly-discovered sag where the hydrocarbon generation potential of its potential source rocks has not been systematically studied.Based on the analysis of organic petrology and organic geochemistry of 160 core samples from three parameter wells and some coalfield exploration wells in the Hesigewula Sag, and combined with tectonic and sedimentary characteristics, the hydrocarbon generation potential of Cretaceous mudstone in the Hesigewula Sag is systematically analyzed.The genetic types of natural gas in Well Hedi 1 are also discussed.The total organic carbon of Lower Cretaceous dark mudstone in Hesigewula Sag is high.The average TOC of K1bt1mudstone is 2.20% and the K1bt2 is 3.67%.The type of the organic matter is Ⅱ1-Ⅱ2.The degree of thermal evolution is relatively low, and it is at the stage of immature-low maturity.The carbon isotope of methane in Well Hedi 1 is relatively light, with contents ranging from -70.5‰ to -68.2‰.Combined with the hydrogen isotope characteristics of methane, it is considered that the natural gas in Well Hedi 1 is biogenic gas, which is produced by the mixture of carbon dioxide reduction and acetic acid fermentation.Because the core samples are generally collected less than 1000 meters at depth, it is presumed that the deep part of the depression(more than 2000 meters) has entered the peak period of oil generation, which has a certain exploration potential.
-
IOCG型矿床，即铁氧化物-铜-金矿床，指铁氧化物含量大于20%的铜-金矿床，主要分布在澳大利亚、巴西、智利、加拿大、美国、南非等国及中国内蒙古、海南、四川、长江中下游等地。其一般具有规模大、品位高、元素多、埋藏浅、易采选等特点，是近年来颇受国内外矿业界、勘查界、矿床界等地质领域重视的一种矿床组合类型[1-7]。由于IOCG型矿床富含铁氧化物，常缺失硫化物，且蚀变范围广阔，地球物理特征明显(成矿区的磁场和重力效应明显，具有重力高、中-高幅磁异常为标志)。因此，高精度磁测往往是寻找IOCG矿床重要的有效手段之一。
对于磁测数据的处理和解译是找矿预测的重要环节。其中RGIS是中国地质调查局发展研究中心牵头研发的物探重磁电数据处理与解释软件系统。与国内外重磁专业数据处理软件相比，RGIS系统具有重磁数据整理齐备、规范，数据预处理功能丰富，重磁数据处理功能全面、实用，数据范围和格式广泛，输入输出方便、快捷、规范等特点。目前已成为中国用户最多、覆盖面最广、产学研多方应用率最高的物探专业软件，在各类地质矿产和油气勘查与评价工作中发挥着重要作用[8-22]。
为探讨高精度磁测在IOCG型铁矿勘查中的应用效果，本文以智利英格瓦塞铁矿为例，在矿区1:5000地面高精度磁测基础上，利用RGIS软件重点针对中部矿区的磁测数据，进行面积上延、面积下延、精测剖面下延及2.5D反演拟合处理，初步确定了矿体的空间分布，为后期地质找矿勘探工作奠定了很好的基础。同时，通过初步对比海南石碌铁矿，发现IOCG型铁矿皆具有富含铁氧化物的特征，磁场特征明显，因而该方法也可进一步推广应用于中国IOCG型铁矿勘查。
1. 区域地质背景
智利英格瓦塞铁矿位于美洲大陆西缘安第斯褶皱山系南美褶皱带构造区内中生代格拉岛弧弧前盆地。区域出露地层主要为中生界上三叠统—下侏罗统水之歌组、下白垩统大十二弦琴组和查尼拉尔组、上白垩统塞里约斯组、新生界古近系—新近系和第四系。区域构造主要由褶皱和断层构成。其中褶皱主要分布在区域西北部；区域内断裂构造相对发育，主要为SN向断层(区域阿塔卡玛断裂近SN向贯穿本区)，其次为NE向、NW向断层。区域内岩浆活动强烈，主要为燕山期中酸性侵入体，岩体贯穿区域南北，规模巨大，呈岩基产出，出露面积占区域面积30%以上。区域内矿产资源丰富，主要矿种为铜、金、铁。铁氧化物铜金型(IOCG)矿床分布广泛(图 1)。
![]()
图 1 智利英格瓦塞铁矿中部矿区地质图(区域位置图据参考文献[1]修改)Figure 1. Geological map of the central mining area of the Incaguasi iron deposit in Chile2. 矿区地质特征
智利英格瓦塞铁矿属于IOCG型矿床，由中部、北部、西南部3个矿区构成，中部矿区是其主要部分，位于智利阿塔卡玛大区巴耶纳尔市英格瓦塞镇北西约1.5 km处，隶属英格瓦塞镇管辖。
矿区地层出露主要为中生界下白垩统大十二弦琴组(其岩性为安山岩)和第四系(冲积物)。矿区地层总体为一单斜构造，走向NE，倾向SE，倾角约为16°，局部有扭曲现象，产状略有变化。由于本区工作程度尚低，未见明显的断裂等构造。区内岩浆岩极发育，主要发育花岗闪长岩岩体、花岗斑岩岩脉及辉石角闪石岩岩体。其中，区内广泛发育花岗闪长岩，其侵入大十二弦琴组中，使地层出露分布较凌乱，同时亦侵入辉石角闪石岩岩体中，野外可见辉石角闪石岩岩体以捕虏体形式赋存于花岗闪长岩岩体中，两者之间界线清晰。辉石角闪石岩主要分布在矿区中部，局部根据含磁铁矿异同，可分为含磁铁辉石角闪石岩、磁铁角闪石岩、角闪磁铁岩。花岗斑岩岩脉主要分布在中部矿区的东北部，总体呈NW向产出，侵入于花岗闪长岩及含磁铁角闪石岩中。此外还出露少量正长岩，呈岩株产出(图 1)。
矿区目前发现4条规模较大的矿体，主要产于含磁铁辉石角闪石岩中，围岩基本为花岗闪长岩。矿体总体呈带状，走向NE，倾向NW，倾角较陡，多70°~90°，单个矿体长约几百米至上千米，宽约几十米至上百米。矿石类型可分为含磁铁辉石角闪石岩、磁铁角闪石岩，含角闪磁铁岩。矿石主要为细粒结构，块状构造。矿石矿物成分较简单，主要为磁铁矿，其次为赤铁矿，局部见孔雀石、蓝铜矿。脉石矿物以角闪、辉石为主。地表浅部围岩矿化与蚀变主要有褐铁矿化、磁铁矿化、赤铁矿化、镜铁矿化、孔雀石化、蓝铜矿化(多在岩石裂隙面发育)、高岭土化、绿泥石化等。
3. 中部矿区磁异常特征
从区域物性参数K(磁化率，单位为4π×10-6SI)和Jr(剩磁强度，单位为10-3A/m)可知，本区花岗闪长岩磁性一般为K=37.5~3046，Jr=14.74~350；花岗斑岩磁性一般为K=23.46~3792，Jr=12.24~228；磁铁角闪石岩为K=2000~120000，Jr=2000~160000；正长斑岩K=0.41~2000，Jr=0.01~30；安山岩K=5~3000，Jr=5~400；含磁铁辉石角闪石岩一般为K=10~10000，Jr=10~3000；辉石角闪石岩一般为K=0.1~1000，Jr=0.1~400。磁性最强的为磁铁角闪石岩，其次为含磁铁辉石角闪石岩，与围岩有明显的磁性差异，因此，本区采用磁法手段勘查非常有效。
根据1:5000地面高精度磁测ΔT等值线图，结合本地区地质及岩石磁性特征，有意义的异常区圈定为大于或等于1600 nT。该工区内的ΔT最大值为5173.3 nT，等值线总体呈NE向延伸。该磁异常区内含7个大于1600 nT的规模较大的磁异常圈，编号分别为Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ、Ⅵ、Ⅶ(图 2)。
![]()
图 2 智利英格瓦塞中部矿区地面磁测ΔT等值线图Figure 2. ΔT contour map of ground magnetic survey for the central mining area of the Incaguasi iron deposit in ChileⅠ号磁异常位于中部矿区的北部，异常中心ΔT最大值为3273.1 nT。该磁异常呈不规则椭圆状，长轴呈NE向，长321 m，短轴最大宽度123 m，规模为40936 m2。该磁异常东南翼等值线较密，梯度变化较大，北西翼梯度变化较小，定性判断引起该磁异常的磁性地质体呈NE走向、向NW向倾斜。
Ⅱ号磁异常位于Ⅰ号磁异常的西南部，异常中心ΔT最大值为2229.3 nT。该磁异常呈不规则椭圆状，长轴呈EW向，长轴长340 m，短轴最大宽度为97 m，规模为39856 m2。该磁异常两翼较对称，初步认定引起该磁异常的磁性地质体呈EW走向、倾角较直立。
Ⅲ号磁异常位于矿区的中西部，异常中心ΔT最大值为2423.8 nT。该异常未完全封闭，长轴呈NE向，长轴长500 m，短轴最大宽度180 m，规模为98163 m2。该异常两翼较对称，初步认定引起该磁异常的磁性地质体呈NE走向、倾角较直立。
Ⅳ号磁异常位于Ⅲ号磁异常的南部，异常中心ΔT最大值为2814.6 nT。该磁异常呈不规则长圆状，长轴呈NE向，长370 m，短轴最大宽度120 m，规模为53993 m2。该磁异常南部为负异常且梯度变化大，初步判定引起该磁异常的磁性地质体呈NE走向、向NW倾斜。
Ⅴ号磁异常位于矿区的南部，异常中心ΔT最大值为2511.3 nT。该异常呈不规则长圆状，长轴呈NE向，长190 m，短轴最大宽度64 m，规模为14396 m2。该磁异常两翼较对称，初步认定引起该磁异常的磁性地质体呈NE走向、倾角较直立。
Ⅵ号磁异常位于矿区的东南部，异常中心ΔT最大值为4339.6 nT。该磁异常呈不规则椭圆状，长轴呈NE向，长660 m，短轴最大宽度170 m，规模为75645 m2。该磁异常南部为负异常且梯度变化大，初步判定引起该磁异常的磁性地质体呈NE走向、向NW倾斜。
Ⅶ号磁异常位于矿区的东北部，异常中心ΔT最大值为5173.3 nT。该磁异常呈椭圆状，长轴呈NE向，长300 m，短轴最大宽度130 m，规模为34740 m2。该磁异常西南翼等值线较密，梯度变化较大，东北部梯度变化较小，初步判定引起该磁异常的磁性地质体呈NE走向、向NW倾斜。
4. 中部矿区磁异常解释
磁异常是各种不同性质、不同深度、不同形态、不同规模的磁性体磁场叠加在一起反应的结果。但是在原始图件上很难识别，给地质解释及验证工作带来了难度。为了进一步提高对磁异常的分辨能力，突出更多有用信息，根据测区磁异常特征和地质解释需要，利用RGIS2010对智利英格瓦塞中部矿区磁测数据进行了面积上延、面积下延、剖面下延及2.5D反演拟合处理。
4.1 面积上延解释
面积上延计算主要目的是压制浅部干扰异常或规模较小的局部异常，突出埋藏较深和规模较大的磁异常。对英格瓦塞中部矿区面积磁测分别进行了向上延拓25 m、50 m、75 m、100 m的数据处理，经统计各矿体ΔT等值线分布情况，如表 1所示。
表 1 各地磁异常上延中心变化统计结果Table 1. Statistical table of the center variation of magnetic anomalies'upward continuation磁异常编号 延拓高度/m 中心异常/nT 与原始曲线差值/nT 每25 m变化值/nT 上延解译 Ⅰ 0 3273.1 0 0 Ⅰ号磁异常区峰值衰减速度快，说明引起该异常的磁性地质体埋藏较浅，规模较小。初步认为该磁性地质体为有限延深 +25 2024.9 -1248.2 -1248.2 +50 1548.2 -1724.9 -476.7 +75 1315.1 -1958 -233.1 +100 1175.3 -2097.8 -139.8 Ⅱ 0 2229.3 0 0 Ⅱ号磁异常区峰值衰减速度较慢，说明引起该异常的磁性地质体埋藏有一定的深度，磁性体有一定的规模。初步认为该磁性地质体为有限延深 +25 1725 -504.3 -504.3 +50 1505.9 -723.4 -219.1 +75 1337.3 -892 -168.6 +100 1277 -952.3 -60.3 Ⅲ 0 2423.8 0 0 Ⅲ号磁异常区峰值衰减速度慢，说明引起该异常的磁性地质体埋藏有深度，磁性体规模较大。初步认为该磁性地质体为有限延深 +25 1934.7 -489.1 -489.1 +50 1710.4 -713.4 -224.3 +75 1557 -866.8 -153.4 +100 1443.7 -980.1 -113.3 Ⅳ 0 2814.6 0 0 Ⅳ号磁异常区峰值衰减速度较慢，说明引起该异常的磁性地质体埋藏有一定深度，磁性体有一定的规模。初步认为该磁性地质体为有限延深 +25 2000.2 -814.4 -814.4 +50 1619.5 -1195.1 -380.7 +75 1410.9 -1403.7 -208.6 +100 1283.8 -1530.8 -127.1 Ⅴ 0 2511.3 0 0 Ⅴ号磁异常区峰值衰减速度非常快，说明引起该异常的磁性地质体埋藏较浅，规模较小。初步认为该磁性地质体为有限延深 +25 1751.2 -760.1 -760.1 +50 异常消失 +75 异常消失 +100 异常消失 Ⅵ 0 4339.6 0 0 Ⅵ号磁异常区峰值衰减速度慢，说明引起该异常的磁性地质体埋藏有深度，磁性体规模大 +25 2667 -1672.6 -1672.6 +50 1933.2 -2406.4 -733.8 +75 1641.1 -2698.5 -292.1 +100 1425.8 -2913.8 -215.3 Ⅶ 0 5173.3 0 0 Ⅶ号磁异常区峰值衰减速度快，说明引起该异常的磁性地质体埋藏浅，磁性体规模大 +25 2946.1 -2227.2 -2227.2 +50 1916.3 -3257 -1029.8 +75 1352.1 -3821.2 -564.2 +100 1013 -4160.3 -338.8 从地面磁测上延25 m、50 m、75 m、100 m的ΔT等值线统计结果(表 1)明显看出，随着向上高度的增加，Ⅰ、Ⅴ磁异常中心值迅速减小，且磁异常消失，说明引起Ⅰ、Ⅴ磁异常的磁性地质体规模小，埋藏浅且为有限延深；Ⅱ、Ⅳ、Ⅶ磁异常中心值衰减较慢，且磁异常有一定的规模，说明引起Ⅱ、Ⅳ、Ⅶ磁异常的磁性地质体有一定的规模，埋藏浅且为有限延深；Ⅲ、Ⅵ磁异常中心值衰减慢，且磁异常规模较大，说明引起Ⅲ、Ⅵ磁异常的磁性地质体有一定规模，埋藏有深度且有一定的延深。
4.2 面积下延解释
向下延拓是由实测磁场向磁源方向的延拓，利用向下延拓可以处理旁侧叠加异常，将多个单独异常从叠加异常中区分出来，同时延拓也可以用于分离上下叠加异常。对英格瓦塞中部矿区面积磁测分别进行了向下延拓25 m、50 m、75 m的数据处理，其磁异常分布情况如下(图 3、图 4、图 5)。
![]()
图 3 智利英格瓦塞中部矿区地面磁测ΔT下延25 m等值线图Figure 3. ΔT contour map with downward continuation 25 m of ground magnetic survey for the central mining area of the Incaguasi iron deposit in Chile![]()
图 4 智利英格瓦塞中部矿区地面磁测ΔT下延50 m等值线图Figure 4. ΔT contour map with downward continuation 50 m of ground magnetic survey for the central mining area of the Incaguasi iron deposit in Chile![]()
图 5 智利英格瓦塞中部矿区地面磁测ΔT下延75 m等值线图Figure 5. ΔT contour map with downward continuation 75 m of ground magnetic survey for the central mining area of the Incaguasi iron deposit in Chile从图 3可以看出，Ⅰ号、Ⅱ号、Ⅲ号、Ⅳ号、Ⅴ号、Ⅵ号磁异常中心在下延25 m时，其磁异常中心值增加速度快，说明延拓深度接近磁性地质体深度；而Ⅶ号磁异常中心在下延25 m时，其中心值增加速度快，且磁异常区及周边场源明显出现凌乱，定性分析已经接近磁性地质体场源上表。由此推断，引起Ⅶ号磁异常的磁性地质体顶板埋深在25 m左右。
从图 4可以看出，Ⅰ号、Ⅳ号磁异常中心在下延50 m时，磁异常区及周边场源明显出现凌乱，说明其已经接近磁性地质体场源上表，结合图 3推断，它们的磁性地质体顶板埋深在25~50 m之间；Ⅱ号、Ⅴ号、Ⅵ号磁异常中心在下延50 m时，磁异常区及周边场源开始出现凌乱，说明延拓深度接近磁性地质体深度；Ⅲ号磁异常中心在下延50 m时，其磁异常中心值增加速度快，说明延拓深度向磁性地质体深度靠近。
从图 5可以看出，Ⅰ号、Ⅳ号磁异常中心在下延75 m时，磁异常周围非常凌乱，证实延拓深度超过磁性地质体顶板；Ⅱ号、Ⅴ号、Ⅵ号磁异常中心在下延75 m时，磁异常区及周边场源出现明显凌乱，定性分析已经接近至磁性地质体场源上表。因此推断引起Ⅱ号、Ⅴ号、Ⅵ号磁异常的磁性地质体顶板埋深在75 m左右；Ⅲ号磁异常中心在下延75 m时，磁异常区及周边场源开始出现凌乱，说明延拓深度接近磁性地质体深度，推断引起Ⅲ号磁异常的磁性地质体顶板埋深在75 m以下。
磁异常面积下延综合表明，Ⅰ号、Ⅳ号磁性地质体顶板埋深在25~50 m之间；Ⅱ号、Ⅴ号、Ⅵ号磁性地质体顶板埋深在75 m左右；Ⅲ号磁性地质体顶板埋深在75 m以下；Ⅶ号磁性地质体顶板埋深在25 m左右。此外，地面磁测ΔT下延25 m、50 m、75 m等值线图显示，有一条NE向呈连续的串珠状的等值线，推测为一条断裂构造(TF1)。
4.3 精测剖面下延与2.5D反演拟合
通过对本区的7个异常8条(242线、246线、254线、274线、290线、292线、296线、336线)精测剖面进行曲线圆滑，高频滤波等数据处理，去掉了地表干扰和一些突跳点。利用RGIS2010对精测剖面进行向下延拓拟推断面处理和2.5D反演拟合，从下延拟推断面图和2.5D拟合结果看，8条精测剖面下延拟推与2.5D拟合结果基本一致(表 2)。
表 2 精测剖面下延拟推和2.5D反演拟合各地磁体空间分布统计结果Table 2. Statistical table of each magnetic bodies' spatial distribution by inferred profile with downward continuation and 2.5D inversion fitting线号 对应磁异常 推测磁性体顶板埋深/m 推测磁性体倾角/° 精测剖面下延拟推 2.5D反演拟合 精测剖面下延拟推 2.5D反演拟合 242 Ⅲ 118 62 75 78 246 Ⅳ 56 60 78 68 254 Ⅲ 60 55 75 75 274 Ⅱ 80 77 70 72 290 Ⅰ 37 42 75 69 292 Ⅵ 70 70 75 73 296 Ⅵ 60 80 75 70 336 Ⅶ 20 18 80 76 例如246线剖面，在对精测剖面进行向下延拓时，首先考虑到随着延拓深度的加大，在有意义的异常得到“放大”的同时，也把某些浅部的干扰或误差随着延拓深度的加大而放大，使延拓曲线发生剧烈的跳动，出现震荡现象，为克服延拓引起的震荡影响，在延拓前先对剖面异常进行“光滑”(即滤波处理)。向下延拓深度与精测剖面的取样点距(即采用数据的剖面长度)有直接关系，一般要求延拓的点距大于延拓的深度，本次取样点距为450 m。采用泊松积分公式进行向下延拓，结合向下延拓实践中出现的规律性，在向下延拓至100 m时，异常出现起伏跳动的震荡现象，就停止向下继续延拓，即到达延拓的深度。同时结合本次地面高精度磁测以大于等于1600 nT的异常确定为有意义的磁异常，推测磁性体顶部埋深56 m，倾角78°。在本剖面野外地质特征和磁异常的基础上，进行了2.5D反演拟合，剖面方位为150°，磁化强度为11271.3×10-2，模型密度为3.35 g/cm3，地磁倾角为为-26.2231°、磁偏角为-1.061°，反演结果磁铁矿体顶部埋深60 m，倾角68°。由此可见，下延拟推断面图和2.5D拟合结果基本一致(图 6、图 7)。
![]()
图 6 智利英格瓦塞中部矿区246线精测剖面曲线及下延拟推断面图Figure 6. Fine profile curve and inferred profile with downward continuation of the 246 surveying line of the central mining area of the Incaguasi iron deposit in Chile![]()
图 7 智利英格瓦塞中部矿区246线精测剖面2.5D反演拟合剖面图Figure 7. Fine profile from 2.5D inversion fitting of the 246 surveying line of central mining area of the Incaguasi iron deposit in Chile需特别说明的是，242、246、254、274、292、336号线的2.5D反演拟合磁异常剖面曲线推测为由2个到几个不等的磁性体引起，若要开展下一步地质工作，应首先对引起磁异常峰值的磁性地质体进行验证。
5. 中国IOCG型矿床勘查简论
自20世纪70年代，在探明的澳大利亚奥林匹克坝铜-铁-金-铀矿床中发现了富铁氧化物这一新特征，受到大家的关注。Hitzman等[23]通过研究澳大利亚奥林匹克坝矿床、加拿大韦尔内克山矿床、瑞典基鲁钠铁矿、美国密苏里东南铁矿及中国白云鄂博矿床，发现它们具有富含铁氧化物矿物组合的共同特征，于1992年将其统称为元古宙铁氧化物(铜-铀-金-稀土)矿床。因该类矿床巨大的潜在经济价值和理论研究意义，也逐渐受到了中国专家学者的关注，并对IOCG矿床进行了研究现状评述[2-3, 24-27]。聂凤军等[3]认为，IOCG型矿床主要产于大陆边缘俯冲带或克拉通盆地(裂谷、张裂带)，容矿围岩多为元古宙、中生代和新生代火山-沉积岩，成矿与A型、I型或磁铁矿系列侵入岩有关，矿体主要为角砾状、脉型、浸染状或块状矿石构成，矿石矿物富铁氧化物(磁铁矿和赤铁矿)而相对贫铁硫化物，围岩蚀变从下往上为钠化带、钾化带和硅化-绢云母化-赤铁矿化带，其矿床成因主要有岩浆流体说、地表/盆地流体说和变质流体说3类[3]。尽管中国对IOCG型矿床的研究刚刚起步，但是一些专家根据其地质特征，建议将长江中下游宁芜庐枞地区的玢岩铁矿、海南石碌铁矿、东准格尔地区的老山口和乔夏哈拉铁-铜-金矿床、内蒙古白云鄂博铁矿、四川拉拉矿床等矿床归为IOCG型矿床[2-3, 27]，展示了很好的资源前景。
随着传统硫化物铁型铜多金属矿床浅部资源的耗竭和深部勘查难度的增大，寻找以氧化铁为主体的铜金矿床(IOCG型矿床)则是一新的勘查方向。富铁氧化物和贫铁硫化物是IOCG型矿床的重要特征，对于勘查IOCG型铁矿而言，利用高精度磁测可以快速地圈定因富含磁铁矿引起的局部磁异常，这也是国内外勘查磁铁矿床最常见、有效的勘查方法。再利用RGIS软件对磁异常进行反演处理分析，推测和圈定矿体埋深和产出状态，能有效地减少中国此类矿床的后期勘查风险，具有很好的应用前景。
比如，中国海南石碌铁矿，其大地构造位于华南褶皱系五指山褶皱带的西段。区内地层主要为青白口系石碌群、震旦系石灰顶组，东南部分布有石炭系南好组、二叠系南龙组、峨查组-鹅顶组等。其中，石碌群中段岩性为条带状透辉透闪石化白云岩、条带状透辉透闪岩、铁质千枚岩等，是主要的赋铁矿围岩。矿区北部、西部、南部三面均被花岗岩侵入。区内褶皱、断裂发育，对赋矿地层和矿体起着重要的控制或改造作用。主矿床主要呈层状、似层状产出。铁矿石中，赤铁矿和石英分别是富铁矿石主要的矿石和脉石矿物，而贫铁矿石矿物主要为赤铁矿和磁铁矿，脉石矿物由石英、透辉石、透闪石、白云石等组成[26, 28-29]。智利英格瓦塞铁矿大地构造位于美洲大陆西缘安第斯褶皱山系南美褶皱带构造区内中生代格拉岛弧弧前盆地，矿区地层主要为中生界下白垩统大十二弦琴组(其岩性为安山岩)，铁矿体主要产于含磁铁辉石角闪石岩中，赋矿围岩基本为花岗闪长岩。矿体呈带状、层状，矿石矿物主要为磁铁矿，次为赤铁矿。虽然中国海南石碌铁矿与智利英格瓦塞铁矿在产状背景、赋矿地层、赋矿围岩、成因有很大差别，但是两者的矿石皆具有富铁氧化物的特征，与围岩的磁性参数差异较大，为高精度磁测奠定了很好的基础。例如，蒋立婷[30]通过高精度磁测在中国海南石碌铁矿矿区东南部如翁—牙加地区发现了7个磁异常，确定了2个铁多金属重点找矿靶区。
6. 结论
(1) 通过对智利英格瓦塞铁矿进行矿区地质、地球物理特征综合研究，利用RGIS软件对高精度磁测数据进行了面积上延、面积下延、精测剖面下延及2.5D反演拟合处理。从结果看，英格瓦塞中部矿区磁异常总规模为371585 m2，引起磁异常的磁性地质体整体呈NE向延伸、倾向NW、倾角较陡。除Ⅰ号和Ⅴ号磁异常规模较小外，其他几个磁异常中心找矿前景较大，有规模的磁异常区大多埋深在50~120 m之间。为本区下一步地质找矿方向提供了很好的依据。
(2) 虽然对IOCG型铁矿床的产出背景、矿床特征、成因机制还有待进一步认识，但是该类矿床富含铁氧化物和明显的地球物理异常却是其显著特征，这为高精度磁测奠定了很好的基础。利用RGIS软件对高精度磁测数据进行处理和解译，可获得磁性体相对可靠的空间产状和分布特征。因而，结合RGIS软件，利用高精度磁测对IOCG型铁矿床的找矿勘查是一种有效的手段，建议在中国大力推广应用。
致谢: 论文在撰写过程中得到中海油研究总院李威的指导与帮助，在此表示衷心的感谢。 -
![]()
图 1 二连盆地贺斯格乌拉地区构造单元[4]
Figure 1. Structural units of the Hesigewula area in Erlian Basin
![]()
图 2 贺斯格乌拉凹陷腾一段(a)和腾二段(b)泥质烃源岩厚度分布
Figure 2. Thickness of argillaceous source rocks in K1bt1(a)and K1bt2(b)in Hesigwula Sag
![]()
图 3 贺地1井、贺地2井和蒙乌参1井腾格尔组烃源岩连井对比
Figure 3. Contrast diagrams of hydrocarbon source rocks of the Tengger Formation in Hedi 1, Hedi 2 and Mengwucan 1 wells
![]()
图 4 有机碳含量(a)和生烃潜量(b)统计直方图
Figure 4. Histogram of total organic carbon(a) and hydrocarbon generation potential(b)
![]()
图 5 腾一段(a)和腾二段(b)烃源岩TOC与氢指数关系
Figure 5. Relation between TOC and HI of source rocks in the in K1bt1(a)and K1bt2(b)
![]()
图 7 干酪根有机元素划分类型范氏图解
Figure 7. Van's diagram of the classification types of organic elements in kerogen
![]()
图 10 贺地2井烃源岩甾烷异构化参数自然演化剖面
Figure 10. Natural evolution profile of sterane isomerization parameters in source rocks of well Hedi 2
表 1 腾一段与腾二段氢指数分布
Table 1 Hydrogen index distribution of K1bt1and K1bt2
地层 腾二段 腾一段 HI/(mg·g-1) 309(84) 519.5(72) 24~727 16~1014 注: 309(84), 309为均值; 84为样品数; 24~727，24为最小值; 727为最大值；腾一段类似 
下载: 导出CSV
表 2 贺地1井天然气稳定碳氢同位素特征
Table 2 Characteristics of stable carbon and hydrogen isotopes of natural gas in well Hedi 1
采样地点及深度 δ13CCH4/‰ δDCH4/‰ 贺地1井 712 m -70.5 -290 716 m -68.2 -281
下载: 导出CSV
-
李世臻, 周新桂, 王丹丹, 等. 内蒙古突泉盆地突参1井原油地球化学特征与油源分析[J]. 地质通报, 2015, 34(10): 1946-1951. doi: 10.3969/j.issn.1671-2552.2015.10.019 刘卫彬, 李世臻, 周新桂, 等. 大兴安岭西缘新发现具有油气远景的中生界凹陷[J]. 中国地质, 2018, 45(3): 634-635. https://www.cnki.com.cn/Article/CJFDTOTAL-DIZI201803019.htm 陈树旺, 公繁浩, 杨建国, 等. 松辽盆地外围油气基础地质调查工程进展与未来工作方向[J]. 中国地质调查, 2016, 3(6): 1-9. https://www.cnki.com.cn/Article/CJFDTOTAL-DZDC201606001.htm Li S Z, Liu W B, Wang D D, et al. Discovery of Hesigewula Sag on the western margin of Da Hinggan Mountains in China and its significance in petroleum geology[J]. China Geology, 2019, 2: 440-458. http://www.sciencedirect.com/science/article/pii/S2096519220300033
李思田, 李宝芳, 杨士恭, 等. 中国东北部晚中生代断陷型煤盆地的沉积作用和构造演化[J]. 地球科学, 1982, (3): 275-294. https://www.cnki.com.cn/Article/CJFDTOTAL-DQKX198203027.htm 李思田, 杨士恭, 吴冲龙, 等. 中国东北部晚中生代裂陷作用和东北亚断陷盆地系[J]. 中国科学(B辑), 1987, (2): 185-195. https://www.cnki.com.cn/Article/CJFDTOTAL-JBXK198702009.htm 云金表, 罗笃清, 李玉喜. 东北地区中生代断陷盆地群构造演化与成油关系探讨[J]. 石油勘探与开发, 1994, (6): 40-45, 115. https://www.cnki.com.cn/Article/CJFDTOTAL-SKYK406.007.htm 周荔青, 刘池阳. 中国东北油气区晚侏罗世-早白垩世断陷油气成藏特征[J]. 中国石油勘探, 2004, (2): 20-25, 1-2. doi: 10.3969/j.issn.1672-7703.2004.02.003 张文浩, 张交东, 周新桂, 等. 中国东北地区东部盆地油气资源前景展望[J]. 地质通报, 2020, 39(2/3)): 307-318. http://dzhtb.cgs.cn/gbc/ch/reader/view_abstract.aspx?file_no=2020020313&flag=1 田在艺, 韩屏. 东北地区中新生代含油气盆地构造分析与形成机制[J]. 河北地质学院学报, 1992, (5): 463-472. https://www.cnki.com.cn/Article/CJFDTOTAL-HBDX199205000.htm 卢双舫, 张敏. 油气地球化学[M]. 北京: 石油工业出版社, 2008: 198-226. 孟庆任, 胡健民, 袁选俊, 等. 中蒙边界地区晚中生代伸展盆地的结构、演化和成因[J]. 地质通报, 2002, (Z1): 224-231. http://dzhtb.cgs.cn/gbc/ch/reader/view_abstract.aspx?file_no=20020464&flag=1 陈建平, 赵长毅, 何忠华. 煤系有机质生烃潜力评价标准探讨[J]. 石油勘探与开发, 1997, (1): 1-5, 91. https://www.cnki.com.cn/Article/CJFDTOTAL-SKYK199701000.htm 李苗苗, 马素萍, 夏燕青, 等. 泌阳凹陷核桃园组湖相烃源岩微观形态特征与形成机制[J]. 岩性油气藏, 2014, 26(3): 45-50. doi: 10.3969/j.issn.1673-8926.2014.03.008 龚大兴, 林金辉, 唐云凤, 等. 上扬子地台北缘古生界海相烃源岩有机地球化学特征[J]. 岩性油气藏, 2010, 22(3): 31-37. doi: 10.3969/j.issn.1673-8926.2010.03.007 邬立言, 顾信章. 热解技术在我国生油岩研究中的应用[J]. 石油学报, 1986, 7(2): 13-19. https://www.cnki.com.cn/Article/CJFDTOTAL-SYXB198602002.htm 戴金星, 陈英. 中国生物气中烷烃组分的碳同位素特征及其鉴别标志[J]. 中国科学(B辑), 1993, (3): 303-310. https://www.cnki.com.cn/Article/CJFDTOTAL-JBXK199303011.htm 陈英, 戴金星. 关于生物气研究中几个理论及方法问题的研究[J]. 石油实验地质, 1994, 16(3): 209-219. https://www.cnki.com.cn/Article/CJFDTOTAL-SYSD403.001.htm 何家雄, 冼仲猷, 陈伟煌, 等. 莺-琼盆地生物气及生物-低成熟过渡带气特征与勘探前景[J]. 中国海上油气(地质), 2002, 16(1): 27-32. doi: 10.3969/j.issn.1673-1506.2002.01.008 Whiticar M J, Faber E, Schoell M. Biogenic methane formation in marine and freshwater environments: CO2 reduction vs. acetate fermentation-isotope evidence[J]. Geochimica et Cosmochimica Acta, 1986, 50(5): 693-709. doi: 10.1016/0016-7037(86)90346-7
Schoell M. The hydrogen and carbon isotopic composition of methane from natural gases of various origins[J]. Geochimica et Cosmochimica Acta, 1980, 44(5): 649-661. doi: 10.1016/0016-7037(80)90155-6
沈平, 徐永昌, 王先彬, 等. 气源岩和天然气地球化学特征及成气机理研究[M]. 兰州: 甘肃科学技术出皈社, 1991, 1: 66-1. 帅燕华, 张水昌, 赵文智, 等. 陆相生物气纵向分布特征及形成机理研究——以柴达木盆地涩北一号为例[J]. 中国科学(D辑), 2007, (1): 46-51. https://www.cnki.com.cn/Article/CJFDTOTAL-JDXK200701004.htm 辛碧霄, 陈一俊. 生物气成因机制与勘探前景[J]. 当代化工, 2018, 47(9): 1928-1931. doi: 10.3969/j.issn.1671-0460.2018.09.046 -
期刊类型引用(2)
1. 安国强，鲁宝亮，高新宇，朱武，李柏森. 边缘特征和深度加权约束的重力三维相关成像反演. 物探与化探. 2024(01): 113-124 . 
百度学术
2. 闫顺尚，王玲，董方营，刘晓，吕晓磊. 基于电-磁法联合勘探的火烧区范围及富水性的预测. 煤田地质与勘探. 2022(02): 132-139 . 百度学术
其他类型引用(3)
计量
- 文章访问数: 2503
- HTML全文浏览量: 453
- PDF下载量: 1550
- 被引次数: 5
