CN113530535A - 基于泥浆气同位素录井的地层热成熟度评估方法和系统 - Google Patents
基于泥浆气同位素录井的地层热成熟度评估方法和系统 Download PDFInfo
- Publication number
- CN113530535A CN113530535A CN202110831457.7A CN202110831457A CN113530535A CN 113530535 A CN113530535 A CN 113530535A CN 202110831457 A CN202110831457 A CN 202110831457A CN 113530535 A CN113530535 A CN 113530535A
- Authority
- CN
- China
- Prior art keywords
- value
- equivalent
- mud gas
- isotope
- vitrinite reflectance
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 43
- 230000015572 biosynthetic process Effects 0.000 title claims abstract description 34
- 239000004079 vitrinite Substances 0.000 claims abstract description 115
- 238000005553 drilling Methods 0.000 claims abstract description 27
- 238000001514 detection method Methods 0.000 claims abstract description 20
- 238000011156 evaluation Methods 0.000 claims abstract description 17
- 239000007789 gas Substances 0.000 claims description 131
- 239000011435 rock Substances 0.000 claims description 56
- 229930195733 hydrocarbon Natural products 0.000 claims description 33
- 150000002430 hydrocarbons Chemical class 0.000 claims description 33
- 239000004215 Carbon black (E152) Substances 0.000 claims description 26
- 239000005416 organic matter Substances 0.000 claims description 26
- 238000002310 reflectometry Methods 0.000 claims description 23
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims description 18
- 229910052799 carbon Inorganic materials 0.000 claims description 18
- 238000004458 analytical method Methods 0.000 claims description 14
- 239000002002 slurry Substances 0.000 claims description 14
- ATUOYWHBWRKTHZ-UHFFFAOYSA-N Propane Chemical compound CCC ATUOYWHBWRKTHZ-UHFFFAOYSA-N 0.000 claims description 12
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 12
- OFBQJSOFQDEBGM-UHFFFAOYSA-N n-pentane Natural products CCCCC OFBQJSOFQDEBGM-UHFFFAOYSA-N 0.000 claims description 12
- 239000000126 substance Substances 0.000 claims description 11
- 230000000155 isotopic effect Effects 0.000 claims description 10
- 229910052739 hydrogen Inorganic materials 0.000 claims description 9
- 239000001257 hydrogen Substances 0.000 claims description 9
- 239000002245 particle Substances 0.000 claims description 9
- 238000000197 pyrolysis Methods 0.000 claims description 9
- YZCKVEUIGOORGS-OUBTZVSYSA-N Deuterium Chemical compound [2H] YZCKVEUIGOORGS-OUBTZVSYSA-N 0.000 claims description 8
- 239000010426 asphalt Substances 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 8
- 229910052805 deuterium Inorganic materials 0.000 claims description 7
- OTMSDBZUPAUEDD-UHFFFAOYSA-N Ethane Chemical compound CC OTMSDBZUPAUEDD-UHFFFAOYSA-N 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 6
- 239000001294 propane Substances 0.000 claims description 6
- 239000001273 butane Substances 0.000 claims description 4
- IJDNQMDRQITEOD-UHFFFAOYSA-N n-butane Chemical compound CCCC IJDNQMDRQITEOD-UHFFFAOYSA-N 0.000 claims description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 4
- 150000002431 hydrogen Chemical class 0.000 claims description 3
- 238000012417 linear regression Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 11
- 230000006870 function Effects 0.000 description 8
- 230000008569 process Effects 0.000 description 8
- 238000004590 computer program Methods 0.000 description 7
- 238000002474 experimental method Methods 0.000 description 7
- 238000011161 development Methods 0.000 description 6
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 238000012545 processing Methods 0.000 description 4
- 238000003860 storage Methods 0.000 description 3
- 241000196324 Embryophyta Species 0.000 description 2
- 150000001335 aliphatic alkanes Chemical class 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000005336 cracking Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- OKTJSMMVPCPJKN-OUBTZVSYSA-N Carbon-13 Chemical compound [13C] OKTJSMMVPCPJKN-OUBTZVSYSA-N 0.000 description 1
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 235000021185 dessert Nutrition 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005194 fractionation Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明公开了一种基于泥浆气同位素录井的地层热成熟度评估方法和系统，方法包括在钻井时开展泥浆气同位素录井获得泥浆气同位素值，取心获得岩心样品；通过分析岩心样品得到地层的等效镜质体反射率值，通过泥浆气同位素值和等效镜质体反射率值之间的变化关系建立两者间的关系式；通过关系式推算出其他探井地层的等效镜质体反射率值，将等效镜质体反射率值与镜质体反射率序列进行比对获得地层的热成熟度。系统包括处理器、内存、泥浆气同位素录井探测装置和等效镜质体反射率值检测装置。本发明在探井时进行泥浆气同位素录井并建立泥浆气同位素值与等效镜质体反射率值之间的关系式推算地层的热成熟度，降低探测成本，提高评估的准确性。
Description
技术领域
本发明涉及资源勘探与开发技术领域，具体涉及一种基于泥浆气同位素录井的地层热成熟度评估方法和系统。
背景技术
石油、天然气和其他地下资源的勘探与开发需要消耗巨大的财力和大量的时间，给寻找例如页岩气、致密甜点区等能够形成商业开发价值的潜在区域带来很多困难。同时，现有技术无法准确区分探测区域的不同地层特征，提供给钻孔选址和后期生产的信息非常有限，也给潜在区域的开发带来不便。热成熟度是评估烃源岩生烃潜力的重要参数之一，因此在开发潜在区域的时候常常使用热成熟度来评估潜在区域的开发价值。常见的例如取芯等技术都是使用地层的热成熟度来区分不同的地层特征的，但是这些技术的实际操作大都效率低下，并且会耗费巨大的时间和财力成本。
在地层热成熟度的评估过程中，通常使用烃类来对地层的热成熟度进行评估。但是在通常情况下，烃类与其他来源的烃类混合、重力分异或者组分分馏后化学指示作用就会受到不可改变的破坏。因此在进行地层热成熟度的评估过程中，一旦烃类从源岩中排出或运移至其他地方，想要再确定地层的热成熟度就变得非常困难。
发明内容
为此，本发明所要解决的技术问题在于克服现有技术中的不足，提供一种在探井时使用泥浆气同位素录井和等效镜质体反射率值建立关系式，并由关系式对探区的地层热成熟度进行评估，降低勘探成本的同时准确获得地层热成熟度的基于泥浆气同位素录井的评估地层热成熟度的方法和系统。
为解决上述技术问题，本发明提供了一种基于泥浆气同位素录井的地层热成熟度评估方法，包括：
步骤1：对探区进行钻井的同时开展泥浆气同位素录井获得泥浆气同位素值，对探井取心获得岩心样品；
步骤2：通过分析岩心样品得到地层的等效镜质体反射率值，
步骤3：根据泥浆气的成分选择泥浆气同位素值，通过泥浆气同位素值和等效镜质体反射率值之间的变化关系建立两者之间的关系式；
步骤4：获取探区内其他探井地层的泥浆气同位素值，通过泥浆气同位素值和等效镜质体反射率值之间的关系式推算出探区内其他探井地层的等效镜质体反射率值，将推算出的等效镜质体反射率值与镜质体反射率序列进行比对获得地层的热成熟度。
进一步地，所述步骤1中对探区进行钻井的同时开展泥浆气同位素录井获得泥浆气同位素值时，对探区进行多口探井的钻探并对每口探井开展泥浆气同位素录井获得泥浆气同位素值，对探井取心获得岩心样品时获取每口探井的岩心样品；
所述步骤3中通过泥浆气同位素值和等效镜质体反射率值之间的变化关系建立两者之间的关系式时，使用的泥浆气同位素值为每口井泥浆气同位素值的平均值，使用的等效镜质体反射率值为每口井等效镜质体反射率值的平均值。
进一步地，所述步骤2中通过分析岩心样品得到地层的等效镜质体反射率值，具体为直接测定岩心样品中有机质颗粒的反射率值，作为岩心样品的等效镜质体反射率值。
进一步地，所述步骤2中通过分析岩心样品得到地层的等效镜质体反射率值，具体为对缺乏有机质颗粒的岩心样品进行岩石热解分析，获得烃类生成最大速率对应的温度值Tmax，判断岩心样品的有机质类型；
若岩心样品的有机质类型为II2-III型有机质，等效镜质体反射率值EqRo的计算公式为：EqRo＝0.016Tmax-6.173；
若岩心样品的有机质类型为I-II1型有机质，等效镜质体反射率值EqRo的计算公式为：EqRo＝0.018Tmax-7.16。
进一步地，所述步骤2中通过分析岩心样品得到地层的等效镜质体反射率值，具体为对岩心样品进行岩石热解分析，获得岩心样品中沥青或动物碎屑的反射率值，根据沥青或动物碎屑对应的等效镜质体反射率换算公式计算得到等效镜质体反射率值。
进一步地，所述步骤3中通过泥浆气同位素值和等效镜质体反射率值之间的变化关系建立两者之间的关系式，具体为通过最小二乘法建立泥浆气同位素值与等效镜质体反射率值之间的线性回归关系式。
进一步地，所述步骤3中根据泥浆气的成分选择泥浆气同位素值，具体为：如果泥浆气包含甲烷、乙烷、丙烷、丁烷和戊烷，选择碳同位素值作为泥浆气同位素值；如果泥浆气包含氢气，选择氘同位素值作为泥浆气同位素值。
进一步地，所述碳同位素值以皮迪箭石作为标准物质来计算获取，所述氘同位素值以标准平均海洋水作为标准物质来计算获取。
本发明还提供了一种基于泥浆气同位素录井的地层热成熟度评估系统，包括处理器、内存、泥浆气同位素录井探测装置和等效镜质体反射率值检测装置，所述内存中存储控制所述处理器的操作指令并通过发送操作指令控制所述处理器工作，所述泥浆气同位素录井探测装置获取钻井时的泥浆气同位素值并传送给所述处理器，所述等效镜质体反射率值检测装置获取探井处的岩心样品的等效镜质体反射率值并传送给所述处理器；
所述处理器根据泥浆气同位素值和等效镜质体反射率值建立关系式，通过关系式推算出探区内其他探井地层的等效镜质体反射率值，对比等效镜质体反射率值与镜质体反射率序列来确定地层的热成熟度。
进一步地，还包括钻机和泥浆振动筛，所述钻机用于在探区钻井，所述泥浆振动筛用于快速脱离探井泥浆中含有的地层气体并将地层气体传送到泥浆气同位素录井探测装置中。
本发明的上述技术方案相比现有技术具有以下优点：
本发明通过泥浆气同位素录井获取泥浆气同位素值并建立泥浆气同位素值与等效镜质体反射率值之间的关系式，通过关系式推算出探区内其他探井地层的等效镜质体反射率值从而得到探区内其它地层的热成熟度，这种方法在得到探区内其它地层的热成熟度时无需进行取心等耗时费钱的技术，降低了勘探成本；同时，在对探区进行钻井的同时就开展了泥浆气同位素录，有效避免了烃类从源岩中排出或运移至其它地方后化学指示作用受到破坏的情况，提高了地层热成熟度评估的准确性。
附图说明
为了使本发明的内容更容易被清楚的理解，下面根据本发明的具体实施例并结合附图，对本发明作进一步详细的说明。
图1是本发明方法的流程图。
图2是本发明系统的结构示意图。
图3是本发明实施例的实验中对三口井的Jd1和Jd3两套烃源岩在侧向上的对比示意图。
图4是本发明实施例的实验中Jd3烃源岩的等效镜质体反射率随深度的变化情况。
图5是本发明实施例的实验中Jd1烃源岩的等效镜质体反射率随深度的变化情况。
图6是本发明实施例的实验中对泥浆气碳同位素值δ13C1、δ13C2和δ13C3采用最小二乘法回归计算，建立等效镜质体反射率与同位素值之间的关系式的结果图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步说明，以使本领域的技术人员可以更好地理解本发明并能予以实施，但所举实施例不作为对本发明的限定。
在本发明的描述中，需要理解的是，术语“包括”意图在于覆盖不排他的包含，例如包含了一系列步骤或单元的过程、方法、系统、产品或设备，没有限定于已列出的步骤或单元而是可选地还包括没有列出的步骤或单元，或可选地还包括对于这些过程、方法、产品或设备固有的其他步骤或单元。
参照图1流程图所示，本发明一种基于泥浆气同位素录井的地层热成熟度评估方法的实施例，包括：
步骤1：对探区进行钻井的同时开展泥浆气同位素录井获得泥浆气同位素值，对探井取心获得岩心样品；同位素值即为同位素指纹，泥浆气同位素值包括碳13同位素值和氢2(氘)同位素值；泥浆气中的含碳和氢同位素值的同位素包括了甲烷C1、乙烷C2、丙烷C3、丁烷C4和戊烷C5。泥浆气同位素录井包括了钻井过程中循环泥浆中脱出的地层气体样品的组成和同位素测定，可以提供多种烷烃类(例如甲烷、乙烷、丙烷等)的碳、氢的同位素数据。
所述对探区进行钻井的同时开展泥浆气同位素录井获得泥浆气同位素值时，对探区进行多口探井的钻探并对每口探井开展泥浆气同位素录井获得泥浆气同位素值，对探井取心获得岩心样品时获取每口探井的岩心样品；通过最小二乘法回归建立泥浆气同位素值与等效镜质体反射率值的关系式时，使用的泥浆气同位素值为每口井泥浆气同位素值的平均值，使用的等效镜质体反射率值为每口井等效镜质体反射率值的平均值。
步骤2：通过分析岩心样品得到地层的等效镜质体反射率值。分析岩心样品的第一种方法为根据标准SY/T5124-2012《沉积岩中镜质体反射率测定方法》直接测定岩心样品中有机质颗粒的反射率值，作为岩心样品的等效镜质体反射率值，不同地质年代的岩石样品其中含有的有机质颗粒是不同的，有机质颗粒为镜质体、沥青或动物碎屑中的一种。分析岩心样品的第二种方法为对缺乏有机质颗粒的岩心样品进行岩石热解分析，获得烃类生成最大速率对应的温度值Tmax，判断岩心样品的有机质类型；若岩心样品的有机质类型为II2-III型有机质，等效镜质体反射率值EqRo的计算公式为：EqRo(％)＝0.016Tmax(℃)-6.173；若岩心样品的有机质类型为I-II1型有机质，等效镜质体反射率值EqRo的计算公式为：EqRo(％)＝0.018Tmax(℃)-7.16。第二种方法一般是在样品中缺少可测定反射率的有机质颗粒时采用。分析岩心样品的第三种方法为对岩心样品进行岩石热解分析，获得岩心样品中沥青或动物碎屑的反射率值，根据等效镜质体反射率换算公式计算得到等效镜质体反射率值。沥青或动物碎屑的等效镜质体反射率换算公式为现有技术，此处不再赘述。本发明中的温度单位为摄氏度。
步骤3：通过同位素录井技术获取泥浆气的同位素组成并根据泥浆气的成分选择泥浆气同位素值，如果泥浆气包含甲烷、乙烷、丙烷、丁烷和戊烷，选择碳同位素值；如果泥浆气包含氢气，选择氘同位素值。所述碳同位素值以皮迪箭石(PDB)作为标准物质来计算获取，皮迪箭石13C/12C值为0.0112372，所述氘同位素值以标准平均海洋水(SMOW)作为标准物质来计算获取，标准平均海洋水D/H值为155.76。以上两种标准物质是由国际原子能机构颁布的国际通用同位素标准。
通过最小二乘法回归建立泥浆气同位素值与等效镜质体反射率值之间的线性回归关系式。
步骤4：获取探区内其他探井地层的泥浆气同位素值，通过泥浆气同位素值和等效镜质体反射率值之间的关系式推算出探区内其他探井地层的等效镜质体反射率值，将推算出的等效镜质体反射率值与镜质体反射率序列进行比对获得地层的热成熟度。镜质体反射率值是热成熟度的一个指标，通过镜质体反射率值可以获得相对应的地层热成熟度。通过这种方法获得探区内其它地层的热成熟度时无需另外开展取心分析，或是其他耗时费钱的分析技术，可以降低时间和财力成本。同时，泥浆气同位素录井和热成熟度确定方法是随钻探实施的，避免了因开展其他技术而延缓钻探进程和增加钻探时间，有效避免烃类受到破坏，提高地层热成熟度评估的准确性。
参照图2的结构示意图所示，本发明一种基于泥浆气同位素录井的地层热成熟度评估系统的实施例，包括处理器、内存、泥浆气同位素录井探测装置和等效镜质体反射率值检测装置，所述内存中存储控制所述处理器的操作指令并通过发送操作指令控制所述处理器工作，所述泥浆气同位素录井探测装置获取钻井时的泥浆气同位素值并传送给所述处理器，所述等效镜质体反射率值检测装置获取探井处的岩心样品的等效镜质体反射率值并传送给所述处理器。所述处理器根据泥浆气同位素值和等效镜质体反射率值建立关系式，通过关系式推算出探区内其他探井地层的等效镜质体反射率值，对比等效镜质体反射率值与镜质体反射率序列来确定地层的热成熟度。
本实施例中的系统还包括钻机和泥浆振动筛，所述钻机用于在探区钻井，所述泥浆振动筛用于快速脱离探井泥浆中含有的地层气体并将地层气体传送到泥浆气同位素录井探测装置中。
为了进一步说明本发明的有益效果，在某探区钻三口探井并使用本发明中的方法进行评估实验。三口探井分别为井1，井2和井3，钻探过程遇到两套烃源岩层系，标记为Jd1和Jd3。Jd1和Jd3两套地层的岩性相同，为含有机质的灰黑色泥岩夹少量薄层粉砂质泥岩或泥质粉砂岩。导眼井钻探的同时，对两套烃源岩层系进行取心、岩心编录、泥浆气同位素录井工作。取心的目的是取得足够的岩心样品进行详细的烃源岩评估研究得到烃源岩评估参数，通过岩心编录判断钻进过程中钻遇地层，通过泥浆气同位素录井得到泥浆气烷烃气的碳同位素值，泥浆气同位素录井数据和烃源岩评估参数用于建立烃源岩成熟度评估方法。
对每口探井的两套烃源岩层段岩心取样进行有机碳和岩石热解分析。有机碳分析获得岩石的总有机碳含量TOC，岩石热解分析获得单位岩石游离烃量S1、干酪根裂解烃量S2、氢指数HI和干酪根裂解生烃最大速率对应的温度Tmax。通过S1，S2和HI三个参数判识干酪根类型。图3所示是三口井的Jd1和Jd3两套烃源岩在侧向上的对比，表1是Jd3和Jd1两套烃源岩岩心分析数据：
| 层位 | TOC(％) | S<sub>1</sub>(mg/grock) | S<sub>2</sub>(mg/grock) | HI(mg/gTOC) | T<sub>max</sub>(℃) |
| Jd3 | 0.49-2.98(1.08) | 0.32-4.33(1.43) | 0.51-4.78(1.88) | 47-198(108) | 434-466(451) |
| Jd1 | 0.40-4.82(1.60) | 0.22-7.62(2.06) | 0.54-9.52(2.53) | 42-174(97) | 445-483(466) |
表1Jd3和Jd1两套烃源岩岩心分析数据从图和表中可以看出，两套烃源岩具有相似的烃源岩品质，TOC含量中等，氢指数较低，属II2-III型；Jd1烃源岩Tmax值略高于Jd3烃源岩，表示热演化程度更高一些。
钻井过程中，碳同位素录井仪对三口井的Jd1和Jd3烃源岩层段的泥浆气进行了甲烷、乙烷和丙烷的碳同位素值测量，表2列出了三口井的两套烃源岩层系碳同位素录井时测量的泥浆气碳同位素值δ13C1、δ13C2和δ13C3的数据。
表2三口井的两套烃源岩层系的δ13C1、δ13C2和δ13C3的数据表
镜质体是指来源于高等植物本质组织的一种有机显微组分，对于缺乏高等植物的地层(比如深湖相沉积、前泥盆纪沉积、碳酸盐岩沉积等)一般采用其他类似镜质体的有机显微组分(如固体沥青、动物碎屑等)来预测热成熟度。这种情况下，其他类似镜质体的有机显微组分的反射率值需要换算为等效镜质体反射率。本实验中Jd1和Jd3两套烃源岩为半深湖-深湖相沉积，缺乏镜质体显微组分，镜质体反射率值无法直接测定。因此，采用岩石热解获得的Tmax值进行换算获得等效镜质体反射率值，对于本实验中的II2-III型有机质，通过公式EqRo＝0.016Tmax-6.173换算得到。经过换算之后，Jd3烃源岩的等效镜质体反射率介于0.92％(井1)至1.17％(井3)，Jd1烃源岩的等效镜质体反射率介于1.13％(井1)至1.43％(井3)。随着深度加深，等效镜质体反射率逐渐变大，图4展示了Jd3烃源岩的等效镜质体反射率随深度的变化情况，图5展示了Jd1烃源岩的等效镜质体反射率随深度的变化情况。
然后，以每口井每套烃源岩等效镜质体反射率平均值和泥浆气碳同位素平均值为X轴和Y轴，绘制图6对表2中的每一组数据(EqRo-δ13C1、EqRo-δ13C2、EqRo-δ13C3)采用最小二乘法回归计算，建立等效镜质体反射率与同位素值之间的关系式的结果图。如果有氢同位素值，同样可以建立等效镜质体反射率与氢同位素之间的关系式。得到等效镜质体反射率(x1)与同位素值δ13C1(y1)之间的关系式为y1＝35.554x1-92.965，等效镜质体反射率(x2)与同位素值δ13C2(y2)之间的关系式为y2＝31.998x2-68.551，等效镜质体反射率(x3)与同位素值δ13C3(y3)之间的关系式为y3＝22.879x3-52.074。
那么，该探区其他钻井的泥浆气同位素录井数据带入关系式可以得到烃源岩的等效镜质体反射率值，快速确定页岩层段的热成熟度和生油气潜力而无需取心开展岩心分析，降低评估成本。
本发明的上述技术方案相比现有技术具有以下优点：本发明通过泥浆气同位素录井获取泥浆气同位素值并建立泥浆气同位素值与等效镜质体反射率值之间的关系式，通过关系式推算出探区内其他探井地层的等效镜质体反射率值从而得到探区内其它地层的热成熟度，这种方法在得到探区内其它地层的热成熟度时无需进行取心等耗时费钱的技术，降低了勘探成本；同时，在对探区进行钻井的同时就开展了泥浆气同位素录，有效避免了烃类从源岩中排出或运移至其它地方后化学指示作用受到破坏的情况，提高了地层热成熟度评估的准确性。
本领域内的技术人员应明白，本申请的实施例可提供为方法、系统、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
显然，上述实施例仅仅是为清楚地说明所作的举例，并非对实施方式的限定。对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其它不同形式变化或变动。这里无需也无法对所有的实施方式予以穷举。而由此所引申出的显而易见的变化或变动仍处于本发明创造的保护范围之中。
Claims (10)
1.一种基于泥浆气同位素录井的地层热成熟度评估方法，其特征在于，包括：
步骤1：对探区进行钻井的同时开展泥浆气同位素录井获得泥浆气同位素值，对探井取心获得岩心样品；
步骤2：通过分析岩心样品得到地层的等效镜质体反射率值，
步骤3：根据泥浆气的成分选择泥浆气同位素值，通过泥浆气同位素值和等效镜质体反射率值之间的变化关系建立两者之间的关系式；
步骤4：获取探区内其他探井地层的泥浆气同位素值，通过泥浆气同位素值和等效镜质体反射率值之间的关系式推算出探区内其他探井地层的等效镜质体反射率值，将推算出的等效镜质体反射率值与镜质体反射率序列进行比对获得地层的热成熟度。
2.根据权利要求1所述的基于泥浆气同位素录井的地层热成熟度评估方法，其特征在于：所述步骤1中对探区进行钻井的同时开展泥浆气同位素录井获得泥浆气同位素值时，对探区进行多口探井的钻探并对每口探井开展泥浆气同位素录井获得泥浆气同位素值，对探井取心获得岩心样品时获取每口探井的岩心样品；
所述步骤3中通过泥浆气同位素值和等效镜质体反射率值之间的变化关系建立两者之间的关系式时，使用的泥浆气同位素值为每口井泥浆气同位素值的平均值，使用的等效镜质体反射率值为每口井等效镜质体反射率值的平均值。
3.根据权利要求1所述的基于泥浆气同位素录井的地层热成熟度评估方法，其特征在于：所述步骤2中通过分析岩心样品得到地层的等效镜质体反射率值，具体为直接测定岩心样品中有机质颗粒的反射率值，作为岩心样品的等效镜质体反射率值。
4.根据权利要求1所述的基于泥浆气同位素录井的地层热成熟度评估方法，其特征在于：所述步骤2中通过分析岩心样品得到地层的等效镜质体反射率值，具体为对缺乏有机质颗粒的岩心样品进行岩石热解分析，获得烃类生成最大速率对应的温度值Tmax，判断岩心样品的有机质类型；
若岩心样品的有机质类型为II2-III型有机质，等效镜质体反射率值EqRo的计算公式为：EqRo＝0.016Tmax-6.173；
若岩心样品的有机质类型为I-II1型有机质，等效镜质体反射率值EqRo的计算公式为：EqRo＝0.018Tmax-7.16。
5.根据权利要求1所述的基于泥浆气同位素录井的地层热成熟度评估方法，其特征在于：所述步骤2中通过分析岩心样品得到地层的等效镜质体反射率值，具体为对岩心样品进行岩石热解分析，获得岩心样品中沥青或动物碎屑的反射率值，根据沥青或动物碎屑对应的等效镜质体反射率换算公式计算得到等效镜质体反射率值。
6.根据权利要求1所述的基于泥浆气同位素录井的地层热成熟度评估方法，其特征在于：所述步骤3中通过泥浆气同位素值和等效镜质体反射率值之间的变化关系建立两者之间的关系式，具体为通过最小二乘法建立泥浆气同位素值与等效镜质体反射率值之间的线性回归关系式。
7.根据权利要求1-6任一项所述的基于泥浆气同位素录井的地层热成熟度评估方法，其特征在于：所述步骤3中根据泥浆气的成分选择泥浆气同位素值，具体为：如果泥浆气包含甲烷C1、乙烷C2、丙烷C3、丁烷C4和戊烷C5，选择碳同位素值；如果泥浆气包含氢气，选择氘同位素值。
8.根据权利要求7所述的基于泥浆气同位素录井的地层热成熟度评估方法，其特征在于：所述碳同位素值以皮迪箭石作为标准物质来计算获取，所述氘同位素值以标准平均海洋水作为标准物质来计算获取。
9.一种基于泥浆气同位素录井的地层热成熟度评估系统，其特征在于：包括处理器、内存、泥浆气同位素录井探测装置和等效镜质体反射率值检测装置，所述内存中存储控制所述处理器的操作指令并通过发送操作指令控制所述处理器工作，所述泥浆气同位素录井探测装置获取钻井时的泥浆气同位素值并传送给所述处理器，所述等效镜质体反射率值检测装置获取探井处的岩心样品的等效镜质体反射率值并传送给所述处理器；
所述处理器根据泥浆气同位素值和等效镜质体反射率值建立关系式，通过关系式推算出探区内其他探井地层的等效镜质体反射率值，对比等效镜质体反射率值与镜质体反射率序列来确定地层的热成熟度。
10.根据权利要求9所述的基于泥浆气同位素录井的地层热成熟度评估系统，其特征在于：还包括钻机和泥浆振动筛，所述钻机用于在探区钻井，所述泥浆振动筛用于快速脱离探井泥浆中含有的地层气体并将地层气体传送到泥浆气同位素录井探测装置中。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110831457.7A CN113530535A (zh) | 2021-07-22 | 2021-07-22 | 基于泥浆气同位素录井的地层热成熟度评估方法和系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110831457.7A CN113530535A (zh) | 2021-07-22 | 2021-07-22 | 基于泥浆气同位素录井的地层热成熟度评估方法和系统 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113530535A true CN113530535A (zh) | 2021-10-22 |
Family
ID=78088627
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110831457.7A Pending CN113530535A (zh) | 2021-07-22 | 2021-07-22 | 基于泥浆气同位素录井的地层热成熟度评估方法和系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113530535A (zh) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20170226851A1 (en) * | 2016-02-10 | 2017-08-10 | Saudi Arabian Oil Company | Thermal Maturity Determination of Rock Formations Using Mud Gas Isotope Logging |
| CN111638317A (zh) * | 2020-05-30 | 2020-09-08 | 中国石油天然气股份有限公司 | 一种大型淡水湖盆陆相泥页岩油形成与演化评价方法 |
| CN112327357A (zh) * | 2019-08-05 | 2021-02-05 | 中国石油天然气股份有限公司 | 基于三维地震资料的烃源岩成熟度预测方法及装置 |
-
2021
- 2021-07-22 CN CN202110831457.7A patent/CN113530535A/zh active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20170226851A1 (en) * | 2016-02-10 | 2017-08-10 | Saudi Arabian Oil Company | Thermal Maturity Determination of Rock Formations Using Mud Gas Isotope Logging |
| CN112327357A (zh) * | 2019-08-05 | 2021-02-05 | 中国石油天然气股份有限公司 | 基于三维地震资料的烃源岩成熟度预测方法及装置 |
| CN111638317A (zh) * | 2020-05-30 | 2020-09-08 | 中国石油天然气股份有限公司 | 一种大型淡水湖盆陆相泥页岩油形成与演化评价方法 |
Non-Patent Citations (4)
| Title |
|---|
| 倪春华;周小进;王果寿;王永成;: "鄂尔多斯盆地富县地区天然气地球化学特征及其成因", 西安石油大学学报(自然科学版), no. 05, 25 September 2010 (2010-09-25) * |
| 刘伟, 李景坤, 宋兰斌, 赵学栋: "汪深1井天然气成因类型及气源分析", 大庆石油地质与开发, no. 01, 1 February 2005 (2005-02-01) * |
| 刘四兵;沈忠民;罗小平;陈义才;: "保山盆地新近系生物气源岩地球化学特征", 新疆石油地质, no. 02, 15 April 2008 (2008-04-15) * |
| 杨秀然: "凯里地区下寒武统碳酸盐岩天然气特征及成因探讨", 贵州地质, no. 03, 30 September 1993 (1993-09-30) * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP3414566B1 (en) | Thermal maturity determination of rock formations using mud gas isotope logging | |
| US10083258B2 (en) | Combining downhole fluid analysis and petroleum systems modeling | |
| CN104730595B (zh) | 一种深层古油藏充注方向和途径示踪方法 | |
| RU2590921C2 (ru) | Способ определения наличия и местоположения подземной углеводородной залежи и происхождения связанных углеводородов | |
| US20090071239A1 (en) | Methods for optimizing petroleum reservoir analysis | |
| MX2014004885A (es) | Metodo para determinar en tiempo real la porosidad y la saturacion de agua de una formacion subterranea usando datos de registro de gas y perforacion. | |
| Li et al. | Prediction of CO2 leakage risk for wells in carbon sequestration fields with an optimal artificial neural network | |
| Li et al. | A new method for predicting sweet spots of shale oil using conventional well logs | |
| Miller et al. | Petrophysics in tight gas reservoirs—key challenges still remain | |
| CN104428492A (zh) | 使用矿物学确定压裂刺激点的布置的系统、方法、和计算机程序产品 | |
| CN109543915B (zh) | 基于测井信息对全井段烃源岩总有机碳含量的识别方法 | |
| Liu et al. | Temperature and pressure characteristics of Ordovician gas condensate reservoirs in the Tazhong area, Tarim Basin, northwestern China | |
| US20130056201A1 (en) | Method for evaluating hydrocarbon-containing subterrean formations penetrated by a directional wellbore | |
| Pinkett et al. | Innovative sidewall pressure coring technology improves reservoir insight in multiple applications | |
| Euzen et al. | Petroleum system analysis using unconventional gas geochemistry: examples from the Montney play of western Canada | |
| Chunhua et al. | Depression and apophysis structure type and its controls on deposition in the Jiyang Depression, Bohai Bay Basin | |
| CN113530535A (zh) | 基于泥浆气同位素录井的地层热成熟度评估方法和系统 | |
| Ungar et al. | Fluid Identification from Mud Gas in the Overburden–A Case Study for the Snorre Field | |
| Walters | Organic geochemistry at varying scales: from kilometres to ångstroms | |
| Abdi et al. | Integrated reservoir characterization and modeling of one Iranian naturally fractured reservoir using laboratory and field data | |
| US20230054254A1 (en) | Reservoir And Production Simulation Using Asphaltene Onset Pressure Map | |
| Melikov et al. | Well Placement and Operation Parameters Optimization of Horizontal Wells in the Development of the PK1 Reservoir of the Kharampurskoye Oil and Gas Condensate Field | |
| Qalandari et al. | Estimation of Rock Permeability Scores Using Machine Learning Methods | |
| Janiga et al. | The isotopic composition of gaseous hydrocarbons–tool for polish shale gas system evaluation | |
| Ahsan et al. | Advanced Gas While Drilling GWD Comparison with Pressure Volume Temperature PVT Analysis to Obtain Information About the Reservoir Fluid Composition, a Case Study from East Kuwait Jurassic Reservoir |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |