CN103776981A - 一种新的岩溶期次判别方法 - Google Patents
一种新的岩溶期次判别方法 Download PDFInfo
- Publication number
- CN103776981A CN103776981A CN201310652473.5A CN201310652473A CN103776981A CN 103776981 A CN103776981 A CN 103776981A CN 201310652473 A CN201310652473 A CN 201310652473A CN 103776981 A CN103776981 A CN 103776981A
- Authority
- CN
- China
- Prior art keywords
- karst
- homogenization temperature
- time
- inclusion
- stuff
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000012850 discrimination method Methods 0.000 title abstract 2
- 238000000034 method Methods 0.000 claims abstract description 46
- CSJDCSCTVDEHRN-UHFFFAOYSA-N methane;molecular oxygen Chemical compound C.O=O CSJDCSCTVDEHRN-UHFFFAOYSA-N 0.000 claims abstract description 43
- 239000012530 fluid Substances 0.000 claims abstract description 39
- 230000000694 effects Effects 0.000 claims abstract description 35
- 238000011161 development Methods 0.000 claims abstract description 19
- 238000000265 homogenisation Methods 0.000 claims description 96
- 239000011435 rock Substances 0.000 claims description 44
- 238000005520 cutting process Methods 0.000 claims description 20
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 claims description 14
- 238000004458 analytical method Methods 0.000 claims description 12
- 238000011160 research Methods 0.000 claims description 10
- 230000004069 differentiation Effects 0.000 claims description 9
- 238000013507 mapping Methods 0.000 claims description 8
- 238000005070 sampling Methods 0.000 claims description 6
- 230000009897 systematic effect Effects 0.000 claims description 5
- 230000000052 comparative effect Effects 0.000 claims description 4
- 230000007797 corrosion Effects 0.000 claims description 4
- 238000005260 corrosion Methods 0.000 claims description 4
- 238000005259 measurement Methods 0.000 claims description 4
- 239000000463 material Substances 0.000 abstract description 9
- 239000007789 gas Substances 0.000 abstract description 8
- 230000009471 action Effects 0.000 abstract description 6
- 238000005516 engineering process Methods 0.000 abstract description 5
- 239000012266 salt solution Substances 0.000 description 11
- 239000000047 product Substances 0.000 description 9
- 229910021532 Calcite Inorganic materials 0.000 description 7
- WUKWITHWXAAZEY-UHFFFAOYSA-L calcium difluoride Chemical compound [F-].[F-].[Ca+2] WUKWITHWXAAZEY-UHFFFAOYSA-L 0.000 description 7
- 229910052799 carbon Inorganic materials 0.000 description 7
- 239000010436 fluorite Substances 0.000 description 7
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 6
- 239000003921 oil Substances 0.000 description 6
- 210000003462 vein Anatomy 0.000 description 6
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 5
- 230000015572 biosynthetic process Effects 0.000 description 5
- 239000010440 gypsum Substances 0.000 description 5
- 229910052602 gypsum Inorganic materials 0.000 description 5
- 229910052760 oxygen Inorganic materials 0.000 description 5
- 239000001301 oxygen Substances 0.000 description 5
- 229910052712 strontium Inorganic materials 0.000 description 5
- CIOAGBVUUVVLOB-UHFFFAOYSA-N strontium atom Chemical compound [Sr] CIOAGBVUUVVLOB-UHFFFAOYSA-N 0.000 description 5
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 4
- 239000013505 freshwater Substances 0.000 description 4
- 230000012010 growth Effects 0.000 description 4
- 229910052500 inorganic mineral Inorganic materials 0.000 description 4
- 239000011707 mineral Substances 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 239000000126 substance Substances 0.000 description 4
- 206010027336 Menstruation delayed Diseases 0.000 description 3
- NINIDFKCEFEMDL-UHFFFAOYSA-N Sulfur Chemical compound [S] NINIDFKCEFEMDL-UHFFFAOYSA-N 0.000 description 3
- 239000005864 Sulphur Substances 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 238000012856 packing Methods 0.000 description 3
- 239000000243 solution Substances 0.000 description 3
- 240000008042 Zea mays Species 0.000 description 2
- 235000005824 Zea mays ssp. parviglumis Nutrition 0.000 description 2
- 235000002017 Zea mays subsp mays Nutrition 0.000 description 2
- OSGAYBCDTDRGGQ-UHFFFAOYSA-L calcium sulfate Chemical compound [Ca+2].[O-]S([O-])(=O)=O OSGAYBCDTDRGGQ-UHFFFAOYSA-L 0.000 description 2
- 239000001569 carbon dioxide Substances 0.000 description 2
- 229910002092 carbon dioxide Inorganic materials 0.000 description 2
- 235000005822 corn Nutrition 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 238000005553 drilling Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 150000007524 organic acids Chemical class 0.000 description 2
- 239000011028 pyrite Substances 0.000 description 2
- NIFIFKQPDTWWGU-UHFFFAOYSA-N pyrite Chemical compound [Fe+2].[S-][S-] NIFIFKQPDTWWGU-UHFFFAOYSA-N 0.000 description 2
- 229910052683 pyrite Inorganic materials 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 102100031102 C-C motif chemokine 4 Human genes 0.000 description 1
- 101100054773 Caenorhabditis elegans act-2 gene Proteins 0.000 description 1
- 102000000584 Calmodulin Human genes 0.000 description 1
- 108010041952 Calmodulin Proteins 0.000 description 1
- 235000019738 Limestone Nutrition 0.000 description 1
- 235000008331 Pinus X rigitaeda Nutrition 0.000 description 1
- 235000011613 Pinus brutia Nutrition 0.000 description 1
- 241000018646 Pinus brutia Species 0.000 description 1
- 244000061458 Solanum melongena Species 0.000 description 1
- QAOWNCQODCNURD-UHFFFAOYSA-L Sulfate Chemical compound [O-]S([O-])(=O)=O QAOWNCQODCNURD-UHFFFAOYSA-L 0.000 description 1
- 238000006065 biodegradation reaction Methods 0.000 description 1
- 210000001124 body fluid Anatomy 0.000 description 1
- 239000010839 body fluid Substances 0.000 description 1
- 238000009933 burial Methods 0.000 description 1
- 239000010430 carbonatite Substances 0.000 description 1
- 238000005234 chemical deposition Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 239000010779 crude oil Substances 0.000 description 1
- 239000013078 crystal Substances 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 210000001163 endosome Anatomy 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 230000007773 growth pattern Effects 0.000 description 1
- 230000003760 hair shine Effects 0.000 description 1
- 238000001027 hydrothermal synthesis Methods 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 230000000155 isotopic effect Effects 0.000 description 1
- 239000006028 limestone Substances 0.000 description 1
- 238000002156 mixing Methods 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 239000002244 precipitate Substances 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 239000013535 sea water Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Landscapes
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本发明提供一种新的岩溶期次判别方法。岩溶储层是近期油气勘探开发新领域,岩溶发育期次的判别一直是国内外岩溶储层建模技术面临的难点。由于现有的岩溶期次判别方法多以暴露标志为判别依据,但岩溶作用形式多样,热液等因素发生的岩溶作用因缺少暴露标志难以判别。本发明针对岩溶期次判别难题,提出一种岩溶充填物碳氧同位素判别充填期次、包裹体均一温度判别岩溶体流体活动期次,两者综合最终实现客观判断岩溶期次的新方法。本发明建立的碳氧同位素判别图版和包裹体均一温度对比图表,可准确快速地判别岩溶的发育期次,突破了传统的岩溶期次判别以暴露标志为依据的局限,解决了地下深层缺乏暴露标志情况下岩溶期次判别的难题。
Description
技术领域
本发明涉及岩溶储层建模中一种新的岩溶期次判别技术,适用于海相碳酸盐岩盆地中岩溶期次的判定。
背景技术
岩溶储层是一种极为重要的油气储层。近期我国在塔里木、四川和鄂尔多斯三大盆地海相碳酸盐岩的岩溶储层获得了一批重要发现,成为碳酸盐岩层系最为重要的油气储层,但岩溶储层建模中岩溶发育期次的判别一直是石油地质界面临的重大难题。
前人在岩溶储层研究方面开展了大量的研究与探索,形成了以暴露标志为主要判识依据的基本方法。常用的方法包括:露头区通过不整合面、溶蚀痕迹等判断岩溶作用以及可能的岩溶期次;地下主要基于地震剖面识别截断、上超和顶超等地质现象,或者基于岩心观察识别暴露标志及岩性突变界面等方法判别岩溶作用强弱及岩溶期次。
然而,对于缺少暴露标志的岩溶作用,如顺层岩溶作用或热液岩溶作用等,就不能利用暴露标志进行岩溶作用研究。由于岩溶***中围岩和岩溶充填物属于不同期产物,岩石地球化学指标反映的是岩溶作用发生后的化学响应,因此可以采用岩石地球化学数据判别岩溶期次。
目前,利用碳氧同位素值判别岩溶充填期次、利用盐水包裹体均一温度判别流体活动期次有报道。其中:
“松辽盆地南部长岭断陷二氧化碳成藏期次”(王颖,李锡瑞,曹跃,李桂俊,段伟;《中国石油勘探》,2009年04期)公开了充分利用包裹体的岩相学特征、均一温度特征、气体组分特征以及碳同位素特征等,对松辽盆地南部长岭断陷二氧化碳气藏的成藏期次进行了分析。但该论文主要应用包裹体进行成藏期次的确定,并没有考虑到利用盐水包裹体得到成岩流体的特征。
“塔河地区奥陶系碳酸盐岩缝洞充填物的地球化学特征及其形成流体分析”(蔡春芳,李开开,李斌,等.《岩石学报》,2009,25(10):2399-2404)公开了塔河地区奥陶系方解石脉和胶结物~(87)Sr/~(86)Sr值多高于围岩与同期海水,显示具有外来的、富~(87)Sr的流体的贡献。富~(87)Sr的流体可来自奥陶系抬升时期淡水对碎屑岩的淋滤作用,但也可来自深部的热液流体。不同成因的流体所沉淀的方解石Sr含量具有明显的差异。地层水δD-δ~(18)O关系证实了存在淡水的混合作用。淡水的混合导致了原油的厌氧生物降解,产生低δ~(34)S的黄铁矿与有机酸;有机酸促进了不整合面附近的岩溶作用。而方解石的均一化温度可达150~190℃,有意义地高于地层所经历的最高温度,支持了存在深部热液流体的活动,说明Sr很可能都来自寒武系或前寒武碎屑岩地层;该热液流体促进了寒武系碳酸盐岩发生热化学硫酸盐还原作用(TSR),产生H_2S气体和黄铁矿,其δ~(34)S值为18‰~22‰。该文中利用了多种地球化学特征,但是对于成岩流体研究方面,关键参数选择不明确。
“塔河油田中下奥陶统岩溶缝洞方解石碳、氧、锶同位素地球化学特征”(刘存革,李国蓉,朱传,等.《地球科学-中国地质大学学报》,2008,33(3):377-386)中,为了探讨塔河油田中下奥陶统古岩溶产物的地球化学特征以及评价不同期次古岩溶作用的影响区域,***分析了中下奥陶统岩溶缝洞方解石的碳、氧、锶同位素。不同产状的方解石碳、氧同位素具有明显的大气水岩溶特征,早期岩溶产物以富集13C和低87Sr/86Sr比值为特征,而海西早期岩溶产物以δ13C分布较宽和高87Sr/86Sr比值为特征,两者的δ18O值分布均较宽.桑塔木组覆盖区比北部地区富集13C,来自围岩的13C比重加大,另一个原因是早期岩溶产物的δ13C值较高;桑塔木组覆盖区岩溶产物具有高87Sr/86Sr比值特征,说明有壳源锶的注入,主要是受海西早期岩溶作用的影响,并且高87Sr/86Sr比值特征比北部地区低,其原因是来源于围岩的重溶锶比重加大.δ13C和87Sr/86Sr在区域上的变化,结合区域地质背景,表明海西早期岩溶对北部地区和桑塔木组覆盖区均有重要的影响,是塔河油田中下奥陶统碳酸盐岩的主岩溶期,也是其储层的主要形成时期。该文利用碳、氧、锶同位素的地球化学特征分析,主要还是充填特征分析及古环境的判别上,并没有联系到溶蚀流体的发生期次,更没有联系到岩溶发育的期次。
“塔里木盆地北缘古岩溶充填物包裹体特征”(夏日元,唐建生,邹胜章,等.《中国岩溶》,2006,25(3):246-249)公开了塔里木盆地北缘露头区奥陶系碳酸盐岩中的古岩溶发育经历了从裸露到埋藏、从淡水作用到热液作用以及从湿润气候到干燥气候的变化,多期次古岩溶作用的叠加导致了充填物地球化学特征的差异性。古岩溶充填物主要有机械沉积充填物、化学淀积充填物和塌积充填物三大类。方解石和萤石包裹体的物理特征、化学相、盐度、均一温度和化学成分均具有多样性,其综合化学特性可划分为四种类型。该文中重点论述包裹体特征,并没有细致论述包裹体均一温度在判别成岩流体性质及期次上。
“矿物包裹体特征对古岩溶作用的指示性”(夏日元,唐建生.《地球学报》,2004,25(3):373-377)公开了碳酸盐岩溶蚀孔缝中充填矿物包裹体特征对古岩溶作用环境和期次具有较好的指示性。黄骅坳陷奥陶系充填矿物包裹体物理特征反映古岩溶作用环境为埋藏型,包裹体化学特征显示了古岩溶作用与有机质熟化程度有关。包裹体各项特征指标的明显差异性与古岩溶作用的多期性相对应:3种盐度类型反映了浅埋藏、中等埋藏和深埋藏3种环境;4个均一温度段与古岩溶发育的4个期次相吻合;该文章中充分利用了包裹体均一温度的作用,并且明确了该均一温度段与古岩溶发育期次的吻合。
《塔里木盆地奥陶系碳酸盐岩缝洞***模式及成因研究》(夏日元,邹胜章,梁彬,等.北京:地质出版社,2011)在深入分析古岩溶形成条件与控制因素基础上,初步揭示了塔里木盆地奥陶系碳酸盐岩古岩溶形成机制与演变规律,***阐明了油气储集体古岩溶缝洞***发育机理以及区域上分区差异性、垂向上分带性和时代分期性,通过古今类比建立了缝洞***发育模式。针对古岩溶发育受多因素控制的特点,建立了综合研究方法体系。该专著中提到了古岩溶研究的多个方面,针对岩溶期次方面也提出了多个方面的影像因素,但是并没有提出一套***的岩溶期次判别方法。
综上所述,在现有技术中,有单独利用包裹体均一温度进行成藏期次厘定和判别成岩流体活动的报道,或是利用碳氧同位素判别充填物性质的报道,或是分析包裹体特征及综合多种地球化学特征研究岩溶特征的报道。然而,这些仅利用单方面因素判别岩溶期次的技术在实际应用中并不可靠,而现有技术中并没有关于将多方面不同因素耦合起来判别岩溶作用期次的技术报道。
发明内容
本发明的目的在于提供一种新的岩溶期次判别方法,在对深层取心较少的情况下,快速准确地判别出岩溶期次,从而预测不同期次下岩溶储层的分布,建立岩溶储层地质模型。
为达上述目的,本发明提供了一种岩溶期次判别方法,该方法主要是从岩溶体充填物碳氧同位素和流体包裹体均一温度分布特征出发,将利用碳氧同位素值判别岩溶充填期次与利用盐水包裹体均一温度判别流体活动期次综合起来,基于充填物碳氧同位素判别充填期次与流体包裹体均一温度判别流体活动期次之间的映射对应关系原理,并综合充填物之间切割关系,客观判别岩溶形成期次。该方法特别适用于海相碳酸盐岩盆地中岩溶期次的判定。
根据本发明的具体实施方案,本发明提供的岩溶期次判别方法包括步骤:
①对岩溶***的围岩及充填物进行采集取样,对采集样品进行碳氧同位素、包裹体特征和均一温度测定;
②根据对样品进行测定得到的δ13C值、δ18O值以及包裹体均一温度,建立δ13C-δ18O交会图、包裹体均一温度频率图;
③根据建立的δ13C-δ18O交会图、包裹体均一温度频率图判别充填物充填期次和流体活动期次。
根据本发明的具体实施方案,本发明的岩溶期次判别方法中,所述步骤①中还包括对样品包裹体特征和盐度值的测定。更具体地,可以是对采集的充填物样品测量δ13C值和δ18O值,其他样品测量包裹体均一温度和盐度等数据。
根据本发明的具体实施方案,本发明的岩溶期次判别方法中,所述步骤①中还包括研究岩溶***充填物、围岩特征及其切割关系等因素。样品采集可基于岩溶***充填物、围岩特征及其切割关系研究,按照***、有代表性采样原则,采集新鲜、无污染样品。
根据本发明的具体实施方案,本发明的岩溶期次判别方法,所述步骤②中,主要在于对测量得到的δ13C值和δ18O值、包裹体均一温度和盐度等数据进行整理分析,形成δ13C-δ18O交会图、包裹体均一温度频率图等基础分析图件。
根据本发明的具体实施方案,本发明的岩溶期次判别方法,所述步骤③中可在根据建立的δ13C-δ18O交会图、包裹体均一温度频率图判别充填物充填期次和流体活动期次的基础上,进一步利用数学中的映射关系原理,结合充填物、围岩特征及其切割关系研究结果,判定岩溶发育期次。
本发明的所述步骤③,主要在于综合了实际地质现象之间的宏观关系和地球化学特征综合判别岩溶发育期次。如充填物之间的宏观切割关系体现在充填物之间的充填顺序,δ13C-δ18O交会图可以判别充填物的充填期次。
根据本发明的具体实施方案,本发明的岩溶期次判别方法还包括步骤:
④结合步骤②和步骤③,建立用于判定岩溶期次的碳氧同位素判别图版和均一温度对比图表。
本发明的所述步骤④中,主要是建立的碳氧同位素判别图版和均一温度对比图表可以作为一种岩溶期次判别图版,对岩溶期次进行判别。
在本发明的一具体实施方案中,本发明的岩溶期次判别方法包括步骤:
①研究岩溶***充填物、围岩特征及其切割关系,对岩溶***内的围岩及充填物进行***取样,对采集样品主要测试碳、氧同位素值和包裹体均一温度值,可选择性记录包裹体特征和盐度值;
②分析整理样品测定得到的δ13C值、δ18O值以及包裹体均一温度等数据,建立δ13C-δ18O交会图、包裹体均一温度频率图等;
③根据建立的δ13C-δ18O交会图、包裹体均一温度频率图等分析图件,判别充填物充填期次和流体活动期次,利用数学中的映射关系原理,结合充填物、围岩特征及其切割关系研究结果,判定岩溶发育期次;
④结合步骤②和步骤③,建立用于判定岩溶期次的碳氧同位素判别图版和均一温度对比图表。
根据本发明的具体实施方案,本发明的岩溶期次判别方法在具体实施时,应用图版进行岩溶期次判别首先通过碳氧同位素图版标定,将测量的碳氧同位素值投影到图版中,判别其所处的划定区域,从而判别出充填期次。
根据本发明的具体实施方案,本发明的岩溶期次判别方法在具体实施时,均一温度对比图表的判别主要是对比研究区与图版划定范围内的平均均一温度值的频率分布图,明确溶蚀所发生的期次,进而可以确定岩溶期次。
本发明所述的岩溶期次判别新方法,利用了大量分析化验资料,结合岩石之间的宏观特征,对投影δ13C-δ18O交会图进行区块划分,从而判别充填期次;通过分析包裹体平均均一温度值,对频率图划分不同的范围,从而判定流体活动期次。建立了用于岩溶期次判别的碳氧同位素判别图版和包裹体均一温度对比图表。方法内涵体现在将多种不同地球化学特征数据有机结合,最终实现对岩溶发育期次的客观判别,提高了岩溶期次的判别精度和准确性。
本案发明人凭借综合多种方法耦合进行储层建模的经验,将多因素耦合。其中,包裹体判别出的流体活动期次可为破坏性的或者建设性的,而碳氧同位素判别出的充填期次可以充填在同一期建设性成岩流体活动下。在本发明中充分利用了数学中的映射关系,具体指充填作用和流体活动期次之间的匹配关系进行映射和对比,结合岩石充填切割之间的宏观关系,进而判别出岩溶发育期次,在实例中具体说明。
具体实施过程中通过对单个包裹体特征分析,建立单个包裹体片下均一温度峰值分布范围,建立多种不同的峰值分布范围,最终通过统计多个包裹体平均均一温度数据建立不同平均均一温度峰值直方图。对于碳氧同位素图版的建立中主要考虑样品选择的纯净性。
综上所述,本发明提出了一种利用岩溶充填物的碳氧同位素值判别岩溶充填期次、利用盐水包裹体均一温度判别岩溶体流体活动期次,综合两者之间对应关系判别岩溶期次新方法。在本发明的一具体实施例中,利用本发明的方法,对塔里木盆地碳酸盐岩岩溶储层进行了研究,相关研究表明,利用本发明建立的碳氧同位素图版和均一温度对比图表,可准确快速地判别岩溶发育的期次,突破了传统的岩溶期次判别多以岩心宏观证据或者地球物理识别暴露标准的局限,为缺少暴露标志区岩溶期次判定提供了有效的图版支持,有效地解决了顺层和热液等岩溶作用期次判别的难题。本发明的方法特别适用于海相碳酸盐岩盆地中岩溶期次的判定。
附图说明
图1本发明中岩溶期次判别方法的流程示意图;
图2本发明中岩溶期次判别方法的实施过程示意图;
图3单个包裹体片中均一温度统计直方图;
图4均一温度范围在75~90℃之间的包裹体特征;
图5均一温度范围在105~115℃之间的包裹体特征;
图6均一温度范围在125~135℃之间的包裹体特征;
图7均一温度范围在145~155℃之间的包裹体特征;
图8均一温度范围在170℃以上的包裹体特征;
图9均一温度分布频率图;
图10碳氧同位素图版;
图11岩石之间的相互切割关系;
图12哈拉哈塘地区碳氧同位素标定图;
图13哈拉哈塘平均均一温度分布频率图;
图14塔河油田碳氧同位素标定图;
图15塔河油田平均均一温度分布频率图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合实施例及附图对本发明的实施方式作进一步地详细描述。
请参见图1和图2所示,本发明提供的岩溶期次判别方法的实施主要包括以下步骤:
①研究岩溶***充填物、围岩特征及其切割关系,对岩溶***内的围岩及充填物进行***取样,对采集样品进行碳氧同位素、古盐度以及包裹体特征和均一温度等测定;
②分析整理样品测定得到的δ13C值、δ18O值以及包裹体均一温度等数据,建立δ13C-δ18O交会图、包裹体均一温度频率图等;
③根据建立的δ13C-δ18O交会图、包裹体均一温度频率图等分析图件,判别充填物充填期次和流体活动期次,利用数学中的映射关系原理,结合充填物、围岩特征及其切割关系研究结果,判定岩溶发育期次;
④结合步骤②和步骤③,建立用于判定岩溶期次的碳氧同位素判别图版和均一温度对比图表。
本发明中以塔里木盆地西北缘大量露头资料为基础,采集了150个碳氧同位素样品,建立碳氧同位素判别图版;其中对应有55个样品进行包裹体分析,建立包裹体均一温度对比图表。具体实施过程见下文。
首先,对单个包体片进行分析主要检测结晶方解石中的流体包裹体,如西克尔地区某一包体片中检测到结晶方解石晶体中两幕盐水包裹体(表1)。其中两幕盐水包裹体中,第一幕盐水包裹体均一温度范围为85.4-94.1℃,平均均一温度为90.0℃;第二幕盐水包裹体均一温度范围为128.1-137.2℃,平均均一温度为132.9℃,从而建立单个包裹体均一温度分布直方图(图3)。
表1单个方解石片中盐水包裹体均一温度测试表
通过建立如图3所示的单个包裹体均一温度分布直方图,利用其峰值可以明确多个包裹体平均均一温度值。结合成岩期次阶段划分温度在15℃左右可划分为不同的期次,将所有的55个包裹体片中的平均均一温度进行统计后发现平均均一温度值分布集中在5个范围内,分别是75-90℃之间、105~115℃之间、125~135℃之间、145~155℃之间和170℃以上几个范围内的包裹体,都可作为不同期次范围内流体活动的证据。经过统计分析,发现不同温度范围内包裹体的形态具有不同的分布特征。如温度范围在75-90℃之间的以椭圆形和近圆形为主;温度范围在105-115℃之间的以圆形和无规则状为主;温度范围在125-135℃之间以圆形为主;温度范围在145-155℃之间以豆荚状和近圆形为主,有少量长条形,相对较少;温度范围在170℃以上无规则状为主(可参见图4~图8所示)。为此,可利用55个包裹体分析数据得到平均均一温度分布统计表(表2),进而建立不同部位不同包裹体均一温度分布频率图(图9)。经统计发现,不同均一温度峰值范围内分布的样本点数量不同。由于包裹体均一温度主要反映流体活动,故样品点数与区域内该期流体活动强度有关,样本点数量越大,可反映出流体活动越强。
表2部分岩石包裹体平均均一温度统计表
对于碳氧同位素处理相对简单,利用150个充填物和围岩的碳氧同位素数据,以充填物类型所处范围为主要依据,将图版区域划分成5个区块,其中以垮塌角砾岩为主的充填物主要分布在区块I,以黄绿色粉砂岩为主的碎屑充填物主要分布在区块II,以***砂岩为主的碎屑物和微晶方解石主要分布在区块III,以巨晶方解石为主的充填物主要分布在区块IV,以萤石和硬石膏为主的热液矿物主要分布在区块V(参见图10)。这样可以确定不同的充填物在不同充填期次下碳氧同位素值的分布范围。
充填物岩石之间的有两种相互切割关系(图11),一种是方解石、萤石、黄绿色粉砂岩3者互相接触,从三者之间的切割关系上,黄绿色砂岩被方解石和萤石切割,萤石切割方解石,这反映三者由早到晚的充填顺序为黄绿色砂岩、方解石、萤石;另一种是碎屑充填物、石膏、硬石膏、硫磺之间的相互接触,其中可看出碎屑充填物主要包括溶洞破碎带或者垮塌体中,而石膏镶嵌在碎屑物表面上,说明碎屑垮塌物充填早于石膏,硫磺主要附着在石膏之上,充填期次要晚于石膏,故可明确三者的充填顺序。这样说明有不同充填期次,再次结合包裹体均一温度峰值分布特征确定的5期流体活动,利用两者之间的映射关系,综合充填物之间的切割关系和地质背景,确定塔里木盆地加里东运动和海西运动作用下5期岩溶作用的碳氧同位素、包裹体均一温度及充填物手标本特征。最终明确了各期岩溶作用下包裹体特征、均一温度值及碳氧同位素分布值,建立了岩溶期次特征判别表(表3)。
(1)均一温度在75~90℃之间的包裹体特征主要以椭圆形和近圆形为主,对应于加里东中期I幕发育的岩溶,且与δ13CPDB值在0~1之间和δ18OPDB值在-4~-9之间所圈定的范围有较好的对应关系。
(2)均一温度在105~115℃之间的包裹体特征主要以圆形和无规则状为主,对应于加里东中期II幕发育的岩溶,且与δ13CPDB值在1~2之间和δ18OPDB值在-4~-10之间所圈定的范围有较好的对应关系。
(3)均一温度在125~135℃之间的包裹体特征主要以圆形为主,少量无规则状,对应于加里东中期III幕发育的岩溶,且与δ13CPDB值在0~-3之间和δ18OPDB值在-4~-10之间所圈定的范围较好关系。
(4)均一温度在145~155℃之间的包裹体特征主要以豆荚状和近圆形为主,有少量长条形,对应于海西早期发育的岩溶,且与δ13CPDB值在0~2.5之间和δ18OPDB值在-10~-14.5之间所圈定的范围有较好的对应关系。
(5)均一温度在170℃以上的包裹体特征主要以无规则状为主,相对较小,对应于海西晚期发育的岩溶,且与δ13CPDB值在0~-4之间和δ18OPDB值在-9.5~-16.5之间所圈定的范围有较好的对应。
表3岩溶期次判别方法指标表
岩溶期次 | δ13CPDB | δ18OPDB | 均一温度(℃) | 手标本标志 |
加里东中期I幕 | 0~1 | -4~-9 | 75-90 | 垮塌角砾岩 |
加里东中期II幕 | 1~2 | -4~-10 | 100-115 | 碎屑充填物 |
加里东中期III幕 | 0~-3 | -4~-10 | 125~135 | 微晶方解石 |
海西早期 | 0~2.5 | -10~-14.5 | 145~155 | 巨晶方解石 |
海西晚期 | 0~-4 | -9.5~-16.5 | 170以上 | 萤石、硫磺等 |
实施例1
以塔里木盆地北部哈拉哈塘地区碳酸盐岩岩溶期次判定为例,前人主要基于地震剖面上的不整合进行判定,由于研究区内钻井普遍较深,钻井较少,故很难在岩心上直接识别出暴露的标志判别岩溶期次。在研究区内的78口井中,共取得22个有效样品,对碳氧同位素值(主要是测δ13CPDB和δ18OPDB)及包裹体均一温度数据进行检测,利用所测数据(表4、表5),将碳氧同位素值标定到通过前面得到的碳氧同位素图版中,其中5个区域内均有分布(图12),说明哈拉哈塘地区也发生过5期充填作用。而利用包裹体均一温度数据建立的平均均一温度分布图(图13)与图9建立的平均均一温度值对比说明哈拉哈塘有5期流体活动。通过上面所述这样就发育5期岩溶,从而得出塔北哈拉哈塘地区同样发育5期岩溶,利用条形图中频率分布的大小进一步明确加里东中期II幕和加里东中期III幕贡献相对较大。这样就达到在数据较少的情况下快速准确地判别岩溶期次,为寻找有利油气储集体提供了基础。
表4哈拉哈塘地区岩石碳氧同位素数据
样品原号 | 岩性 | δ13CPDB | δ18OPDB | δ18OSMOW | Z值 |
H601-18-3 | 泥晶砂岩灰岩 | 0.4 | -8 | 22.6 | 124.14 |
H601-18-4 | 泥晶砂屑灰岩 | 0.4 | -8.6 | 22 | 123.84 |
H801-1 | 方解石 | -3.6 | -13.3 | 17.2 | 113.30 |
H801-2 | 泥晶砂屑灰岩 | -0.5 | -6.9 | 23.8 | 122.84 |
H12-2-3 | 亮晶砂屑灰岩 | -0.7 | -7.6 | 23.1 | 122.08 |
H7-3-2 | 深灰色亮晶砂屑灰岩 | 0.5 | -6.8 | 23.9 | 124.94 |
H11-1-1 | 碎屑物 | 1 | -6 | 24.7 | 126.36 |
H902-3 | 泥晶砂屑 | 0.6 | -5.1 | 25.6 | 125.99 |
H7-1 | 方解石脉体 | -3.9 | -12.7 | 17.8 | 112.99 |
H803-3 | 泥质充填物 | 1.2 | -6.3 | 24.4 | 126.62 |
X8H-1 | 方解石脉 | -2 | -15 | 15.5 | 115.73 |
X8H-2 | 方解石脉 | 0.1 | -14.6 | 15.8 | 120.23 |
H601-1-3 | 方解石脉 | 0.5 | -8.8 | 21.8 | 123.94 |
H6-1-4 | 方解石脉 | -1.9 | -16 | 14.4 | 115.44 |
R4 | 方解石 | -1.4 | -10.6 | 19.9 | 119.15 |
表5哈拉哈塘部分盐水包裹体均一温度值分布
实施例2
以塔河油田为例进行岩溶判定期次说明。相对哈拉哈塘而言,塔河油田的岩心相对较多,且可在部分岩心上识别出暴露标志,但仍然不足以为岩溶期次的判别提供直接的证据。依据本发明方法,利用其中32个数据碳氧同位素值(表6)标定到建立的碳氧同位素图版中,发现分布在其中的4个区域(图14),说明塔河油田发生过4期充填作用。而利用包裹体均一温度数据(表7)建立的平均均一温度分布图(图15)与图9对比说明塔河发生过4期流体活动。
表6塔河油田部分岩石碳氧同位素数据
井号 | 岩性 | δ13C(PDB) | δ18O(PDB) | 井号 | 岩性 | δ13C(PDB) | δ18O(PDB) |
AD20 | 方解石 | -1.6 | -15.4 | AD16 | 角砾岩 | -0.1 | -8.2 |
AD20 | 方解石 | -1.2 | -14.9 | AD16 | 角砾岩 | -0.5 | -3.8 |
AD20 | 方解石 | -2.7 | -18.7 | S82 | 角砾岩 | -0.5 | -6.9 |
AD3 | 方解石 | -1.7 | -14.5 | S82 | 角砾岩 | -1.3 | -3.8 |
AD3 | 方解石 | -1.5 | -14.5 | T728 | 方解石 | 1.4 | -12.3 |
AD12 | 方解石 | -2.1 | -17.9 | S82 | 方解石 | 0.9 | -11.3 |
AD12 | 方解石 | -2 | -16.8 | T728 | 方解石 | 0.7 | -12.6 |
S85 | 方解石 | -3 | -12.8 | T740 | 角砾岩 | -1.5 | -7 |
S85 | 方解石 | -2.3 | -17 | YQ3 | 角砾岩 | -1.8 | -4.5 |
AD16 | 方解石 | -3.1 | -17.5 | AD12 | 角砾岩 | -1.9 | -7.3 |
TP8 | 碎屑物质 | -2.2 | -5.9 | AD3 | 角砾岩 | -2.1 | -7.1 |
AD16 | 碎屑物质 | 1.1 | -6 | AD15 | 方解石 | -0.8 | -7.9 |
S82 | 碎屑物质 | 1.4 | -4.8 | AD20 | 方解石 | -0.5 | -7.1 |
TP8 | 碎屑物质 | 1.3 | -7.3 | AD16 | 方解石 | -0.8 | -4.3 |
AD12 | 碎屑物质 | 2 | -6.1 | YQ3 | 方解石 | -1.1 | -8.7 |
T728 | 碎屑物质 | 2.1 | -7.9 | AD16 | 方解石 | -1.3 | -4.5 |
表7塔河油田部分盐水包裹体均一温度值分布
由在本实施例得到的结论可得出塔河油田所在区块共发育4期岩溶,利用条形图中频率分布的大小可进一步明确加里东中期III幕和海西晚期是主要岩溶发育期次。该实例同样达到数据少的情况下判别出岩溶期次。
Claims (8)
1.一种岩溶期次判别方法,该方法包括步骤:
①对岩溶***的围岩及充填物进行采集取样,对采集样品进行碳氧同位素和均一温度测定;
②根据对样品进行测定得到的δ13C值、δ18O值以及包裹体均一温度数据,建立δ13C-δ18O交会图、包裹体均一温度频率图;
③根据建立的δ13C-δ18O交会图、包裹体均一温度频率图判别充填物充填期次和流体活动期次。
2.根据权利要求1所述的岩溶期次判别方法,其特征在于,所述步骤①中还包括对样品包裹体特征和盐度值的测定。
3.根据权利要求1所述的岩溶期次判别方法,其特征在于,所述步骤①中还包括研究岩溶***充填物、围岩特征及其切割关系,所述步骤③中在根据建立的δ13C-δ18O交会图、包裹体均一温度频率图判别充填物充填期次和流体活动期次的基础上,进一步利用数学中的映射关系原理,结合充填物、围岩特征及其切割关系研究结果,判定岩溶发育期次。
4.根据权利要求1所述的岩溶期次判别方法,该方法还包括步骤:
④结合步骤②和步骤③,建立用于判定岩溶期次的碳氧同位素判别图版和均一温度对比图表。
5.根据权利要求1所述的岩溶期次判别方法,该方法包括如下步骤:
①研究岩溶***充填物、围岩特征及其切割关系,对岩溶***的围岩及充填物进行***取样,对采集样品进行碳氧同位素、包裹体特征和均一温度测定;
②分析整理样品测定得到的δ13C值、δ18O值以及包裹体均一温度、盐度等数据,建立δ13C-δ18O交会图、包裹体均一温度频率图等;
③根据建立的δ13C-δ18O交会图、包裹体均一温度频率图判别充填物充填期次和流体活动期次,利用数学中的映射关系原理,结合充填物、围岩特征及其切割关系研究结果,判定岩溶发育期次;
④结合步骤②和步骤③,建立用于判定岩溶期次的碳氧同位素判别图版和均一温度对比图表。
6.根据权利要求4或5所述的岩溶期次判别方法,其特征在于,应用图版进行岩溶期次判别首先通过碳氧同位素图版标定,将测量的碳氧同位素值投影到图版中,判别其所处的划定区域,从而判别出充填期次。
7.根据权利要求4或6所述的岩溶期次判别方法,其特征在于,均一温度对比图表的判别主要是对比研究区与图版划定范围内的平均均一温度值的频率分布图,明确溶蚀所发生的期次,进而可以确定岩溶期次。
8.根据权利要求1所述的岩溶期次判别方法,该方法是用于海相碳酸盐岩盆地中岩溶期次的判定。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310652473.5A CN103776981B (zh) | 2013-12-05 | 2013-12-05 | 一种新的岩溶期次判别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310652473.5A CN103776981B (zh) | 2013-12-05 | 2013-12-05 | 一种新的岩溶期次判别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103776981A true CN103776981A (zh) | 2014-05-07 |
CN103776981B CN103776981B (zh) | 2015-08-05 |
Family
ID=50569475
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310652473.5A Active CN103776981B (zh) | 2013-12-05 | 2013-12-05 | 一种新的岩溶期次判别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103776981B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104655821A (zh) * | 2015-01-29 | 2015-05-27 | 中国石油天然气股份有限公司 | 一种碳酸盐岩储层成因判识方法 |
CN105626058A (zh) * | 2015-12-30 | 2016-06-01 | 中国石油天然气股份有限公司 | 一种确定储层岩溶发育程度的方法及装置 |
CN107966443A (zh) * | 2017-12-25 | 2018-04-27 | 中国石油大学(华东) | 一种地层古流体盐度演化史的研究方法 |
CN108267795A (zh) * | 2017-12-18 | 2018-07-10 | 中国石油天然气股份有限公司 | 岩溶垮塌角砾的确定方法和装置 |
CN108709770A (zh) * | 2018-07-05 | 2018-10-26 | 中国地质大学(武汉) | 一种单一期次油包裹体群体组分取样***及方法 |
CN109343121A (zh) * | 2018-10-25 | 2019-02-15 | 西安石油大学 | 一种确定热接触变质岩储层中岩浆热液规模的方法 |
CN110261467A (zh) * | 2019-07-01 | 2019-09-20 | 西南石油大学 | 一种识别碳酸盐岩古岩溶储层垂向分带性的方法 |
CN110320340A (zh) * | 2018-03-29 | 2019-10-11 | 中国石油化工股份有限公司 | 蒸发岩与石油硫同位素对比判断石油烃源岩时代的方法 |
CN113916880A (zh) * | 2020-06-23 | 2022-01-11 | 中国石油天然气股份有限公司 | 确定成岩流体演化历史的方法 |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109781966B (zh) * | 2019-03-14 | 2021-06-25 | 西安石油大学 | 一种海相碳酸盐岩储层成岩古流体温度校正图版的建立方法 |
-
2013
- 2013-12-05 CN CN201310652473.5A patent/CN103776981B/zh active Active
Non-Patent Citations (4)
Title |
---|
夏日元 等: "塔里木盆地北缘古岩溶充填物包裹体特征", 《中国岩溶》, vol. 25, no. 3, 30 September 2006 (2006-09-30) * |
张涛 等: "锶同位素在塔河古岩溶期次划分中的应用", 《石油实验地质》, vol. 27, no. 3, 30 June 2005 (2005-06-30) * |
汪泽成 等: "我国海相碳酸盐岩大油气田油气藏类型及分布特征", 《石油与天然气地质》, vol. 34, no. 2, 30 April 2013 (2013-04-30) * |
谢发川: "四川盆地西—北部中二叠统热液作用及对储层的影响", 《中国优秀硕士学位论文全文数据库 信息科技辑》, no. 1, 15 January 2012 (2012-01-15), pages 138 - 1136 * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104655821B (zh) * | 2015-01-29 | 2016-08-24 | 中国石油天然气股份有限公司 | 一种碳酸盐岩储层成因判识方法 |
CN104655821A (zh) * | 2015-01-29 | 2015-05-27 | 中国石油天然气股份有限公司 | 一种碳酸盐岩储层成因判识方法 |
US10641088B2 (en) | 2015-12-30 | 2020-05-05 | Petrochina Company Limited | Method and device for determining karst development degree of reservoir, computer readable storage medium and device |
CN105626058A (zh) * | 2015-12-30 | 2016-06-01 | 中国石油天然气股份有限公司 | 一种确定储层岩溶发育程度的方法及装置 |
CN108267795A (zh) * | 2017-12-18 | 2018-07-10 | 中国石油天然气股份有限公司 | 岩溶垮塌角砾的确定方法和装置 |
CN107966443A (zh) * | 2017-12-25 | 2018-04-27 | 中国石油大学(华东) | 一种地层古流体盐度演化史的研究方法 |
CN110320340A (zh) * | 2018-03-29 | 2019-10-11 | 中国石油化工股份有限公司 | 蒸发岩与石油硫同位素对比判断石油烃源岩时代的方法 |
CN110320340B (zh) * | 2018-03-29 | 2021-11-09 | 中国石油化工股份有限公司 | 蒸发岩与石油硫同位素对比判断石油烃源岩时代的方法 |
CN108709770A (zh) * | 2018-07-05 | 2018-10-26 | 中国地质大学(武汉) | 一种单一期次油包裹体群体组分取样***及方法 |
CN108709770B (zh) * | 2018-07-05 | 2023-12-01 | 中国地质大学(武汉) | 一种单一期次油包裹体群体组分取样***及方法 |
CN109343121A (zh) * | 2018-10-25 | 2019-02-15 | 西安石油大学 | 一种确定热接触变质岩储层中岩浆热液规模的方法 |
CN110261467A (zh) * | 2019-07-01 | 2019-09-20 | 西南石油大学 | 一种识别碳酸盐岩古岩溶储层垂向分带性的方法 |
CN110261467B (zh) * | 2019-07-01 | 2021-09-07 | 西南石油大学 | 一种识别碳酸盐岩古岩溶储层垂向分带性的方法 |
CN113916880A (zh) * | 2020-06-23 | 2022-01-11 | 中国石油天然气股份有限公司 | 确定成岩流体演化历史的方法 |
CN113916880B (zh) * | 2020-06-23 | 2023-12-26 | 中国石油天然气股份有限公司 | 确定成岩流体演化历史的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103776981B (zh) | 2015-08-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103776981B (zh) | 一种新的岩溶期次判别方法 | |
Blake et al. | Compositional multivariate statistical analysis of thermal groundwater provenance: A hydrogeochemical case study from Ireland | |
CN103196807B (zh) | 一种砂岩成岩过程与孔隙演化的分析方法 | |
Huang et al. | Shale gas generation and potential of the lower Cambrian Qiongzhusi formation in the southern Sichuan Basin, China | |
CN106168685B (zh) | 一种页岩气单井地质综合评价方法 | |
CN104747183B (zh) | 一种碳酸盐岩储层综合分类方法 | |
CN101943669B (zh) | 低场核磁共振测量钻井液含油率的方法 | |
CN104655821B (zh) | 一种碳酸盐岩储层成因判识方法 | |
CN103670383A (zh) | 一种识别泥页岩油藏有效储层的方法及设备 | |
CN104199121A (zh) | 一种页岩气藏建产有利区的综合判别方法 | |
CN104806232B (zh) | 一种确定碳酸盐岩储层孔隙度下限的方法 | |
Cui et al. | Oil-bearing heterogeneity and threshold of tight sandstone reservoirs: A case study on Triassic Chang7 member, Ordos Basin | |
Shuai et al. | Key oil accumulation periods of ultra-deep fault-controlled oil reservoir in northern Tarim Basin, NW China | |
CN105240005A (zh) | 一种低孔特低渗非常规储层溶蚀孔识别方法 | |
CN107102377B (zh) | 定量预测致密砂岩油气有利勘探区的方法 | |
CN106199727A (zh) | 一种应用地震层速度识别页岩气甜点的方法 | |
Li et al. | The formation environment and developmental models of argillaceous dolomite in the Xingouzui Formation, the Jianghan Basin | |
CN105134185A (zh) | 储层流体性质识别方法 | |
Hou et al. | New insights into loess formation on the southern margin of the Chinese Loess Plateau | |
Su et al. | Effect of lithofacies on shale reservoir and hydrocarbon bearing capacity in the Shahejie Formation, Zhanhua Sag, eastern China | |
CN105158802A (zh) | 坳陷湖盆重力流沉积物测井定量识别方法 | |
CN105257284B (zh) | 一种利用元素俘获能谱测井确定凝灰质含量的方法及装置 | |
Šliaupa et al. | Petrophysical and mechanical properties of the lower Silurian perspective oil/gas shales of Lithuania | |
Liu et al. | Reef-carbonate diagenesis in the Pleistocene–Holocene of the well Xike# 1, Xisha Islands, South China Sea: Implications on sea-level changes | |
CN104516025A (zh) | 碳酸盐储层物性随钻分类和评价方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |