CN106249300A - 基于井震结合确定源岩toc含量和非均质性方法和装置 - Google Patents
基于井震结合确定源岩toc含量和非均质性方法和装置 Download PDFInfo
- Publication number
- CN106249300A CN106249300A CN201610803915.5A CN201610803915A CN106249300A CN 106249300 A CN106249300 A CN 106249300A CN 201610803915 A CN201610803915 A CN 201610803915A CN 106249300 A CN106249300 A CN 106249300A
- Authority
- CN
- China
- Prior art keywords
- source rock
- toc content
- hydrocarbon source
- measured
- district
- 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
- 239000011435 rock Substances 0.000 title claims abstract description 282
- 238000000034 method Methods 0.000 title claims abstract description 59
- 239000004215 Carbon black (E152) Substances 0.000 claims abstract description 256
- 229930195733 hydrocarbon Natural products 0.000 claims abstract description 256
- 150000002430 hydrocarbons Chemical class 0.000 claims abstract description 256
- 238000009826 distribution Methods 0.000 claims abstract description 77
- 230000007774 longterm Effects 0.000 claims abstract description 19
- 238000005516 engineering process Methods 0.000 claims abstract description 14
- 238000005553 drilling Methods 0.000 claims description 58
- 230000015572 biosynthetic process Effects 0.000 claims description 31
- 238000000819 phase cycle Methods 0.000 claims description 21
- 238000005070 sampling Methods 0.000 claims description 13
- 230000018199 S phase Effects 0.000 claims description 11
- 238000010183 spectrum analysis Methods 0.000 claims description 11
- 238000002474 experimental method Methods 0.000 claims description 10
- 230000004044 response Effects 0.000 claims description 9
- 230000009897 systematic effect Effects 0.000 claims description 6
- 230000011218 segmentation Effects 0.000 claims description 4
- 230000035945 sensitivity Effects 0.000 claims description 2
- IYLGZMTXKJYONK-ACLXAEORSA-N (12s,15r)-15-hydroxy-11,16-dioxo-15,20-dihydrosenecionan-12-yl acetate Chemical compound O1C(=O)[C@](CC)(O)C[C@@H](C)[C@](C)(OC(C)=O)C(=O)OCC2=CCN3[C@H]2[C@H]1CC3 IYLGZMTXKJYONK-ACLXAEORSA-N 0.000 abstract 1
- IYLGZMTXKJYONK-UHFFFAOYSA-N ruwenine Natural products O1C(=O)C(CC)(O)CC(C)C(C)(OC(C)=O)C(=O)OCC2=CCN3C2C1CC3 IYLGZMTXKJYONK-UHFFFAOYSA-N 0.000 abstract 1
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 7
- 230000008569 process Effects 0.000 description 7
- 229910052799 carbon Inorganic materials 0.000 description 6
- 230000000694 effects Effects 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 239000005416 organic matter Substances 0.000 description 4
- 238000010276 construction Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 238000005065 mining Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 239000003208 petroleum Substances 0.000 description 2
- 238000004064 recycling Methods 0.000 description 2
- 239000004575 stone Substances 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- 230000001154 acute effect Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000006854 communication Effects 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 238000002939 conjugate gradient method Methods 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000011478 gradient descent method Methods 0.000 description 1
- 239000010439 graphite Substances 0.000 description 1
- 229910002804 graphite Inorganic materials 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000003129 oil well Substances 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 239000003079 shale oil Substances 0.000 description 1
- 238000002945 steepest descent method Methods 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
- G01V1/50—Analysing data
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明实施例提供了一种基于井震结合确定源岩TOC含量和非均质性方法和装置，其中，该方法包括：获取待测区的地质和实验数据；根据地质和实验数据正演得到待测区的烃源岩TOC含量的特征曲线；进而利用地震波形反演得到待测区的烃源岩TOC含量的反演数据体；在待测区的层序格架下，根据烃源岩TOC含量的反演数据体，利用地层切片技术，分别确定待测区长期旋回最大湖泛面位置、中期旋回最大湖泛面位置和短期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图。通过该方案可以有效解决现有方法中存在的确定烃源TOC含量精度差，准确性低，并且容易遗漏部分高TOC含量烃源岩层段的技术问题，从而提高了确定烃源岩TOC含量的精度和准确性。
Description
技术领域
本发明涉及油气勘探技术领域，特别涉及一种确定烃源岩TOC(Total OrganicCarbon，总有机碳)含量的方法和装置。
背景技术
在油气勘探领域中，高TOC含量(高丰度)烃源岩分布往往直接影响着油气的分布，尤其在页岩油气和致密油气的勘探中，这种规律非常明显。但是，由于烃源岩TOC含量具有较强的非均质性，使得利用地震资料对烃源岩TOC含量进行精细预测比较困难。因此，在实际的勘探工作中，如何利用地震资料，对待测区的烃源岩TOC含量进行精细预测是油气勘探中的重要工作。
目前，实际施工时，大都是根据波阻抗或地震属性与烃源岩TOC含量的关系，利用间接反演波阻抗或者地震属性的方法来预测烃源岩TOC含量。
然而，上述利用间接反演波阻抗或者地震属性的方法预测烃源岩TOC含量，在具体实施时往往由于方法本身的局限性，存在精度较低、准确性较差的问题，尤其在无井区进行预测时，这种方法很难在精度上满足实际的勘探需要；另外，利用上述方法在具体预测烃源岩TOC含量时，大多只是对目的层段求取算术平均，这样的处理往往还会导致遗漏部分资源潜力较好的高TOC含量的烃源岩层段，导致最终预测结果不准确。
针对上述问题，目前尚未提出有效的解决方案。
发明内容
本发明实施例提供了一种基于井震结合确定源岩TOC含量和非均质性方法和装置，以解决目前采用间接反演波阻抗或者地震属性的方法预测烃源岩TOC含量时存在的确定烃源TOC含量精度差，准确性低，并且容易遗漏部分高TOC含量烃源岩层段的技术问题，从而提高了确定烃源岩TOC含量的精度和准确性。
本发明实施例提供了一种基于井震结合确定源岩TOC含量和非均质性方法，包括：
获取待测区的地质和实验数据；
根据所述地质和实验数据正演得到所述待测区的烃源岩TOC含量的特征曲线；
以所述待测区的层序格架和所述烃源岩TOC含量的特征曲线为约束，利用地震波形反演得到所述待测区的烃源岩TOC含量的反演数据体；
在所述待测区的层序格架下，根据所述烃源岩TOC含量的反演数据体，利用地层切片技术，分别确定所述待测区长期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图，中期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图和短期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图。
在一个实施方式中，所述获取待测区的地质和实验数据，包括：
对所述待测区内的钻井的岩心和岩屑进行系统采样，得到采样品；
根据所述采样品的TOC含量实验，得到所述采样品中烃源岩TOC含量数据；
以所述采样品的X衍射及能谱分析实验数据和所述钻井的测井数据为约束，从所述烃源岩TOC含量数据中剔除砂岩的TOC含量数据；
将所述采样品中剔除砂岩的TOC含量数据后剩余的烃源岩TOC含量数据作为所述待测区的地质和实验数据。
在一个实施方式中，对所述待测区内的钻井的岩心和岩屑进行系统采样的采样间距设置为1m到2m。
在一个实施方式中，根据所述地质和实验数据正演得到所述待测区的烃源岩TOC含量的特征曲线，包括：
根据所述采样品中剔除砂岩的TOC含量数据后剩余的烃源岩TOC含量数据，参考全岩X衍射和能谱分析实验，进行岩性精细归位；
结合所述钻井的测井曲线的纵向变化特征，确定所述待测区内的钻井纵向上烃源岩TOC含量的测井响应特征；
根据所述烃源岩TOC含量的测井响应特征，做交汇图分析，确定烃源岩TOC含量的敏感测井曲线；
根据所述敏感测井曲线和所述地质和实验数据，通过分段或者分层系确定烃源岩TOC含量的测井预测模型；
根据所述测井预测模型，正演得到所述烃源岩TOC含量的特征曲线。
在一个实施方式中，以所述待测区的层序格架和所述烃源岩TOC含量的特征曲线为约束，利用地震波形反演得到所述待测区的烃源岩TOC含量的反演数据体，包括：
在所述层序格架下，以所述烃源岩TOC含量特征曲线作为约束，根据所述待测区的三维地震资料，利用所述待测区的地震波形，建立反演模型；
根据所述反演模型，求解得到所述待测区的烃源岩TOC含量的反演数据体。
在一个实施方式中，利用所述待测区的地震波形，建立反演模型，包括：
根据所述层序格架内的地震解释层位和所述待测区内的钻井数量，确定样本探井数和高频成分；
根据所述样本探井数和所述高频成分，建立初始模型；
根据所述初始模型，反演得到所述待测区的烃源岩TOC含量的反演数据体；
将所述烃源岩TOC含量的反演数据体中连井反演剖面上的烃源岩TOC含量和基于所述钻井的测井数据得到的所述钻井的烃源岩TOC含量特征曲线进行比较，得到比较误差；
如果所述比较误差小于等于预设阈值，则将所述初始模型确定为所述反演模型。
在一个实施方式中，在将所述烃源岩TOC含量的反演数据体中连井反演剖面上的烃源岩TOC含量和基于所述钻井的测井数据得到的所述钻井的烃源岩TOC含量特征曲线进行比较，得到比较误差之后，还包括：
如果所述比较误差大于所述预设阈值，则重新确定样本探井数和高频成分，并根据重新确定的样本探井数和高频成分建立初始模型，直到根据建立的初始模型，通过地震波形反演得到的所述烃源岩TOC含量的反演数据体中连井反演剖面上的烃源岩TOC含量和基于所述钻井的测井数据得到的所述钻井的烃源岩TOC含量特征曲线之间的比较误差小于等于所述预设阈值，并将当前的初始模型作为所述反演模型。
在一个实施方式中，在所述待测区的层序格架下，根据所述烃源岩TOC含量的反演数据体，利用地层切片技术，分别确定所述待测区长期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图，中期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图和短期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图，包括：
在所述待测区的层序格架下，利用地层切片技术分别获取所述待测区所述长期旋回最大湖泛面位置的沿层地层切片、所述中期旋回最大湖泛面位置的沿层地层切片和所述短期旋回最大湖泛面位置的沿层地层切片；
利用所述烃源岩TOC含量的反演数据体，根据所述长期旋回最大湖泛面位置的沿层地层切片确定所述长期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图，根据所述中期旋回最大湖泛面位置的沿层地层切片确定所述中期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图，根据所述短期旋回最大湖泛面位置的沿层地层切片确定所述短期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图；
其中，所述长期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图用于表征所述长期旋回最大湖泛面位置的TOC含量的空间分布特征，所述中期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图用于表征所述中期旋回最大湖泛面位置的TOC含量的空间分布特征，所述短期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图用于表征所述短期旋回最大湖泛面位置的TOC含量的空间分布特征。
本发明实施例还提供了一种基于井震结合确定源岩TOC含量和非均质性装置，包括：
第一获取模块，用于获取待测区的地质和实验数据；
第二获取模块，用于根据所述地质和实验数据正演得到所述待测区的烃源岩TOC含量的特征曲线；
第三获取模块，用于以所述待测区的层序格架和所述烃源岩TOC含量的特征曲线为约束，利用地震波形反演得到所述待测区的烃源岩TOC含量的反演数据体；
确定模块，用于在所述待测区的层序格架下，根据所述烃源岩TOC含量的反演数据体，利用地层切片技术，分别确定所述待测区长期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图，中期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图和短期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图。
在一个实施方式中，所述第三获取模块包括：
模型建立单元，用于在所述层序格架下，以所述烃源岩TOC含量特征曲线作为约束，根据所述待测区的三维地震资料，利用所述待测区的地震波形，建立反演模型；
数据体求解单元，用于根据所述反演模型，求解得到所述待测区的烃源岩TOC含量的反演数据体。
在本发明实施例中，由于充分利用了待测区的地震波形与待测区烃源岩TOC含量的关系，通过正演得到烃源岩的TOC含量特征曲线，进而再反演得到烃源岩的TOC含量反演数据体，考虑到TOC的非均质性，基于TOC含量反演数据体分别选择长期旋回、中期旋回和短期旋回最大湖泛面位置做沿层地层切片，进而确定这三个位置的烃源岩TOC含量等值线平面分布图，从而解决了现有方法中存在的确定烃源TOC含量精度差，准确性低和容易遗漏有效开采区的技术问题，达到了提高确定可采油气准确性的技术效果。
附图说明
此处所说明的附图用来提供对本发明的进一步理解，构成本申请的一部分，并不构成对本发明的限定。在附图中：
图1是根据本发明实施例的基于井震结合确定源岩TOC含量和非均质性方法的处理流程图；
图2是根据本发明实施例的基于井震结合确定源岩TOC含量和非均质性方法中的建立反演模型的处理流程图；
图3是根据本发明实施例的基于井震结合确定源岩TOC含量和非均质性装置的结构示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白，下面结合实施方式和附图，对本发明做进一步详细说明。在此，本发明的示意性实施方式及其说明用于解释本发明，但并不作为对本发明的限定。
考虑到目前主要采用间接反演波阻抗或者地震属性的方法预测烃源岩TOC含量时，由于现有方法本身技术上的局限性，又没有充分利用待测区的地震波形与待测区烃源岩TOC含量之间的关系，导致预测烃源岩TOC含量的结果精度低、准确性差，甚至无法满足实际的勘探需求，尤其在无井区，这种问题尤为严重。另外，又由于目前采用的预测烃源岩TOC含量的方法大都只是对目的层段求取算术平均，这种简单的求平均处理往往容易导致遗漏部分资源潜力较好的高TOC含量的烃源岩层段。针对产生上述技术问题的根本原因，本发明实施例考虑到实际施工时，虽然井与井之间及无井区的烃源岩TOC含量变化不能直接得到，但是上述区域的地震波形是已知的，可以充分利用待测区的地震波形与待测区的烃源岩TOC含量的关系，通过地震波形指示反演方法，根据反演模型确定待测区的烃源岩TOC含量反演数据体，从而，可以提高预测的精度和准确性；同时通过求取不同层序级别下，分布较为稳定的最大湖泛面位置处的烃源岩TOC含量沿层地层的切片处烃源岩TOC含量的空间变化，从而，解决现有方法中由于未充分考虑烃源岩TOC含量的非均质性，导致的容易遗漏部分资源潜力较好的高TOC含量烃源岩层段的技术问题。
基于上述思路，本发明实施例提供一种基于井震结合确定源岩TOC含量和非均质性方法，具体实施步骤可以如图1所示，包括：
步骤101：获取待测区的地质和实验数据；
上述实施方式中的TOC(Total Organic Carbon，总有机碳)是衡量岩石中有机质的丰度的主要指标之一，是指岩石中存在于有机质中的碳。它不包括碳酸盐岩、石墨中的无机碳，通常用占岩石重量的％来表示。其中，有机质丰度是指单位质量岩石中有机质的数量，在其它条件相近的前提下，一般岩石中有机质的丰度(含量)越高，其生烃能力越高。
具体实施时，一般测井数据可以通过测井记录获得，实验数据一般要先通过对待测区的测井进行采样，再对采样得到的才样品进行相关实验，从而得到上述待测区的地质和实验数据，具体可以按照以下步骤执行：对待测区内的钻井的岩心和岩屑进行系统采样，得到采样品；根据采样品的TOC含量实验，得到采样品中烃源岩TOC含量数据；以采样品的X衍射及能谱分析实验数据和该待测区内钻井的测井数据为约束，从烃源岩TOC含量数据中剔除砂岩的TOC含量数据；将采样品中剔除砂岩的TOC含量数据后剩余的烃源岩TOC含量数据作为该待测区的地质和实验数据。
在一个实施方式中，为了使得采样得到的采样品能够较好地反应待测研究区的实际情况，采样间隔的间距，一般可以设置为1m到2m。然而，值得注意的是，上述所列举的间距仅是一种示意性描述，是为了更好地说明本发明，在实际执行的过程中，还可以包括其它间距，只要基于所设置的间距可以保证采样得到采样品能较好的反应待测区的实际情况，可以根据实际需要进行选取，例如，为了得到更精确到数据，也可以将采样间距设置在1m以内，对此，本申请对此不作限定。
需要说明的是，当所使用的测井是油井时，由于含油级别较高的砂岩也会具有一定含量的TOC，会干扰本发明实施例中通过上述采样品的TOC含量实验确定的采样品中烃源岩的TOC的数据，影响最终的预测结果，因此在获取待测区的地质和实验数据时，要事先通过对采样品进行X衍射和能谱分析实验，以剔除采样品中的砂岩TOC，即：对采样品进行X衍射和能谱分析实验；根据X衍射和能谱分析实验得到的实验结果，剔除采样品中的砂岩TOC。
步骤102：根据所述地质和实验数据正演得到所述待测区的烃源岩TOC含量的特征曲线；
为了获得烃源岩TOC含量特征曲线，本发明实施例先利用实验数据和测井数据得到测井预测模型，再利用得到的测井预测模型，正演得到烃源岩TOC含量特征曲线。具体地，获得待测区的烃源岩TOC含量的特征曲线的实施步骤可以包括：根据采样品中剔除砂岩的TOC含量数据后剩余的烃源岩TOC含量数据，参考全岩X衍射和能谱分析实验，进行岩性精细归位；结合钻井的测井曲线的纵向变化特征，确定待测区内的钻井纵向上烃源岩TOC含量的测井响应特征；根据烃源岩TOC含量的测井响应特征，做交汇图分析，确定烃源岩TOC含量的敏感测井曲线；根据敏感测井曲线和地质和实验数据，通过分段或者分层系确定烃源岩TOC含量的测井预测模型；根据测井预测模型，正演得到烃源岩TOC含量的特征曲线。
其中，上述实施方式中的待测区内的钻井纵向上烃源岩TOC含量的测井响应特征，具体可以包括：高自然伽马、高声波时差、高电子、高中子和低密度。
需要说明的是，通过上述实施例得到的测井预测模型实一般会有很多个，具体实施时，本发明实施例只优选最符合要求的测井预测模型作为本发明实施例的使用的测井预测模型，这样可以使得后续步骤得到的结果尽可能的准确。具体选择时，一般将通过最小二乘法得到的测井预测模型作为符合要求的，并能够在本发明实施例中使用的测井预测模型，这种测井预测模型，相对于其他的测井预测模型往往拟合精度更高，其相关性一般可以达到0.9以上，通过这种测井预测模型正演得到的烃源岩TOC含量特征曲线往往也是符合要求，相对比较接近实际情况的曲线。
步骤103：以待测区的层序格架和烃源岩TOC含量的特征曲线为约束，利用地震波形反演得到该待测区的烃源岩TOC含量的反演数据体；
考虑到现有的采用间接反演波阻抗或者地震属性方法预测烃源岩TOC含量时由于没有利用地震波形与烃源岩TOC关系，导致的确定烃源TOC含量精度差，准确性低的技术问题，又考虑到实际施工时，尽管井与井之间及无井区的烃源岩TOC含量不能直接获得，但该地区的地震波形是已知的，通过寻找地震波形与正演得到的烃源岩TOC特征曲线的关系，可以较为准确的预测出待测区的烃源岩TOC含量空间分布。基于上述考虑，考虑可以通过地震波指示反演技术确定待测区的烃源岩TOC含量反演数据体，具体实施可以包括：在层序格架下，以烃源岩TOC含量特征曲线作为约束，根据待测区的三维地震资料，利用待测区的地震波形，建立反演模型；根据反演模型，求解得到该待测区的烃源岩TOC含量的反演数据体。
上述实施方式中的地震波形指示反演技术是指一种利用地震波形的相关特性进行地震勘探最终定位地下分布的油气资源的方法。实际进行地震勘探时一般的是在地面上人工的激发地震波，由于介质的非均匀性，当地震波在地层介质中向各方向传播时会产生反射、衍射、散射和透射，部分的地震波返回到地面就构成了所接收的地震观测数据，从这些观测数据中提取地层剖面及介质的物性参数，进而确定油气藏的位置。具体处理时，通常采用弹性动力学中的各种波动方程来描述地震波在地下介质中的传播过程；地震观测数据则可由波动方程的正演来模拟，而后可以用地震波形反演来重构地层剖面及物性参数。其中，地震波形反演方法利用叠前地震波场的运动学和动力学信息重建地层结构，具有揭示复杂地质背景下构造与储层物性的潜力，且实现过程相对比较准确、高效。目前应用较为广泛的有：Gauss-Newton法、Newton型法、梯度下降法(最速下降法)以及共轭梯度法。
需要说明的是，本发明实施例中通过地震波形反演求解是区别于现有方法，得到更加准确的待测区的烃源岩TOC含量反演数据体是一个重要步骤，而在上述具体实施方式中，利用地震波形反演确定待测区的烃源岩TOC含量反演数据体的一个重要内容又是建立反演模型，实际上也正是基于该反演模型才能得到相较于现有方法更加准确的烃源岩TOC含量反演数据体，因此建立该准确的反演模型是本发明实施例的一个重要内容，该反演模型建立的具体步骤可以如图2所示，具体包括：
步骤201：根据层序格架内的地震解释层位和待测区内的钻井数量，确定样本探井数和高频成分；
步骤202：根据样本探井数和高频成分，建立初始模型；
步骤203：根据初始模型，反演得到待测区的烃源岩TOC含量的反演数据体；
步骤204：将烃源岩TOC含量的反演数据体中连井反演剖面上的烃源岩TOC含量和基于钻井的测井数据得到的钻井的烃源岩TOC含量特征曲线进行比较，得到比较误差；
步骤205：如果该比较误差小于等于预设阈值，则将初始模型确定为反演模型。
具体实施时，如果步骤204得到的比较误差大于预设阈值时，还需要对初始模型进行反复修改调整，直到比较误差小于等于预设阈值，说明这时调整后的初始模型是符合要求接近实际情况准确的模型，可以将这时的模型作为反演模型。具体的实施步骤，可以包括：如果比较误差大于预设阈值，则重新确定样本探井数和高频成分，并根据重新确定的样本探井数和高频成分建立初始模型，直到根据建立的初始模型，通过地震波形反演得到的烃源岩TOC含量的反演数据体中连井反演剖面上的烃源岩TOC含量和基于钻井的测井数据得到的钻井的烃源岩TOC含量特征曲线之间的比较误差小于等于预设阈值，并将当前的初始模型作为反演模型。
实际实施时，为了能针对性对初始模型进行修改，使之尽快满足预设要求，可以根据比较误差，利用全局优化算法对初始模型进行修改。具体操作可以包括：根据比较误差，有针对性地修改样品数和高频成分，得到修改后的初始模型，再利用得到的修改后的初始模型重复步骤203、步骤204和步骤205，直到比较误差小于等于预设阈值。
步骤104：在待测区的层序格架下，根据烃源岩TOC含量的反演数据体，利用地层切片技术，分别确定待测区长期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图，中期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图和短期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图。
上述实施方式中的最大湖泛面(maximum flooding surface)最初是指在一个湖平面的变化周期中的最高湖平面，这时湖泊面积最广，泥岩分布面积也最广，应用到地震勘探领域，具体在井上识别湖泛面就是找广泛分布的泥岩。一般某个层序级别下的最大湖泛面位置处往往对应该层序中的分布最广、分布最稳定的位置，该位置处的烃源岩TOC含量一般是该层序中相对最高的。
考虑到现有的采用间接反演波阻抗或者地震属性方法预测烃源岩TOC含量时，没有充分地考虑到烃源岩TOC含量的非均质特性的影响，往往只是简单地对目的层段求取算术平均，因此容易导致遗漏部分资源潜力较好的高TOC含量烃源岩层段。针对这种情况，为了解决上述技术问题，本发明实施例考虑了烃源岩TOC含量的非均质性，利用地层切片技术，求取不同层序级别下分布最为稳定的最大湖泛面位置处的烃源岩TOC含量，再根据最大湖泛面位置处的烃源岩TOC含量确定待测区的烃源岩TOC含量等值线分布图。为了获得待测区长期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图，中期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图和短期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图，具体实施时，可以包括：
S1：在待测区的层序格架下，利用地层切片技术分别获取待测区长期旋回最大湖泛面位置的沿层地层切片、中期旋回最大湖泛面位置的沿层地层切片和短期旋回最大湖泛面位置的沿层地层切片；
S2：利用烃源岩TOC含量的反演数据体，根据长期旋回最大湖泛面位置的沿层地层切片确定长期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图，根据中期旋回最大湖泛面位置的沿层地层切片确定中期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图，根据短期旋回最大湖泛面位置的沿层地层切片确定短期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图；
其中，长期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图用于表征长期旋回最大湖泛面位置的TOC含量的空间分布特征，中期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图用于表征中期旋回最大湖泛面位置的TOC含量的空间分布特征，短期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图用于表征短期旋回最大湖泛面位置的TOC含量的空间分布特征。
上述实施方式中的TOC含量的非均质性可以是指在一些地质情况中，由于沉积环境和后期演化等因素的制约，导致烃源岩存在比较强烈的非均质性，进而导致烃源岩的TOC含量也相应地存在非均质性特点。这种有机质分布不均匀的特点，很容易会造成对储层中烃源岩TOC含量的预测错误，进而对烃源岩评价和资源的预测带来严重的影响。
在本发明实施例中，由于充分利用了待测区的地震波形与待测区烃源岩TOC含量的关系，通过正演得到烃源岩的TOC含量特征曲线，进而再反演得到烃源岩的TOC含量的反演数据体，又考虑到TOC的非均质性，基于TOC含量反演数据体分别选择长期旋回、中期旋回和短期旋回最大湖泛面位置做沿层地层切片，确定这三个位置的烃源岩TOC含量等值线平面分布图，从而解决了现有方法中存在的确定烃源TOC含量精度差，准确性低和容易遗漏有效开采区的技术问题，达到了提高确定可采油气准确性的技术效果。
基于同一发明构思，本发明实施例中还提供了一种基于井震结合确定源岩TOC含量和非均质性装置，如下面的实施例所述。由于基于井震结合确定源岩TOC含量和非均质性装置解决问题的原理与基于井震结合确定源岩TOC含量和非均质性方法相似，因此基于井震结合确定源岩TOC含量和非均质性装置的实施可以参见基于井震结合确定源岩TOC含量和非均质性方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。图3是本发明实施例的基于井震结合确定源岩TOC含量和非均质性装置的一种结构框图，如图3所示，确定烃源岩TOC含量的装置具体可以包括：第一获取模块301，第二获取模块302，第三获取模块303，确定模块304，下面对该结构各个模块进行具体说明。
第一获取模块301，用于获取待测区的地质和实验数据；
第二获取模块302，用于根据地质和实验数据正演得到待测区的烃源岩TOC含量的特征曲线；
第三获取模块303，用于以待测区的层序格架和烃源岩TOC含量的特征曲线为约束，利用地震波形反演得到待测区的烃源岩TOC含量的反演数据体；
确定模块304，用于在待测区的层序格架下，根据烃源岩TOC含量的反演数据体，利用地层切片技术，分别确定待测区长期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图，中期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图和短期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图。
其中，长期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图用于表征长期旋回最大湖泛面位置的TOC含量的空间分布特征，中期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图用于表征中期旋回最大湖泛面位置的TOC含量的空间分布特征，短期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图用于表征短期旋回最大湖泛面位置的TOC含量的空间分布特征
具体实施时，上述实施例中的第一获取模块301可以按照以下步骤获取待测区的实验数据：
对待测区内的钻井的岩心和岩屑进行系统采样，得到采样品；
根据采样品的TOC含量实验，得到采样品中烃源岩TOC含量数据；
以采样品的X衍射及能谱分析实验数据和钻井的测井数据为约束，从烃源岩TOC含量数据中剔除砂岩的TOC含量数据；
将采样品中剔除砂岩的TOC含量数据后剩余的烃源岩TOC含量数据作为该待测区的地质和实验数据。。
实际操作时，为了根据实验数据和测井数据，建立测井预测模型，并通过测井预测模型得到烃源岩TOC含量特征曲线，第二获取模块302具体可以按照下面步骤执行具体操作：
根据采样品中剔除砂岩的TOC含量数据后剩余的烃源岩TOC含量数据，参考全岩X衍射和能谱分析实验，进行岩性精细归位；
结合钻井的测井曲线的纵向变化特征，确定待测区内的钻井纵向上烃源岩TOC含量的测井响应特征；
根据烃源岩TOC含量的测井响应特征，做交汇图分析，确定烃源岩TOC含量的敏感测井曲线；
根据敏感测井曲线和地质和实验数据，通过分段或者分层系确定烃源岩TOC含量的测井预测模型；
根据测井预测模型，正演得到所述烃源岩TOC含量的特征曲线。
为了求解得到待测区较为准确的烃源岩TOC含量反演数据体，第三获取模块303具体可以按照下面步骤执行：
在层序格架下，以烃源岩TOC含量特征曲线作为约束，根据待测区的三维地震资料，利用待测区的地震波形，建立反演模型；
根据反演模型，求解得到待测区的烃源岩TOC含量的反演数据体。
其中，第三获取模块303建立上述反演模型过程，具体又可以包括：
根据层序格架内的地震解释层位和待测区内的钻井数量，确定样本探井数和高频成分；
根据样本探井数和高频成分，建立初始模型；
根据初始模型，反演得到待测区的烃源岩TOC含量的反演数据体；
将烃源岩TOC含量的反演数据体中连井反演剖面上的烃源岩TOC含量和基于钻井的测井数据得到的钻井的烃源岩TOC含量特征曲线进行比较，得到比较误差；
如果比较误差小于等于预设阈值，则将初始模型确定为反演模型；
如果比较误差大于所述预设阈值，则重新确定样本探井数和高频成分，并根据重新确定的样本探井数和高频成分建立初始模型，直到根据建立的初始模型，通过地震波形反演得到的烃源岩TOC含量的反演数据体中连井反演剖面上的烃源岩TOC含量和基于钻井的测井数据得到的钻井的烃源岩TOC含量特征曲线之间的比较误差小于等于预设阈值，并将当前的初始模型作为反演模型。
其中，上述步骤中，如果比较误差大于所述预设阈值，通过第三获取模块对初始模型进行一次或多次修改，具体修改可以包括：根据比较误差，针对性地修改样本探井数和高频成分，建立修改后的初始模型。
为了最终确定待测区的烃源岩TOC含量等值线分布图，确定模块304具体可以按照下面步骤处理：
在待测区的层序格架下，利用地层切片技术分别获取待测区长期旋回最大湖泛面位置的沿层地层切片、中期旋回最大湖泛面位置的沿层地层切片和短期旋回最大湖泛面位置的沿层地层切片；
利用烃源岩TOC含量的反演数据体，根据长期旋回最大湖泛面位置的沿层地层切片确定长期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图，根据中期旋回最大湖泛面位置的沿层地层切片确定中期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图，根据短期旋回最大湖泛面位置的沿层地层切片确定短期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图；
其中，长期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图用于表征长期旋回最大湖泛面位置的TOC含量的空间分布特征，中期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图用于表征中期旋回最大湖泛面位置的TOC含量的空间分布特征，短期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图用于表征短期旋回最大湖泛面位置的TOC含量的空间分布特征。
从以上描述中，可以看出，本发明实施例由于充分利用了待测区的地震波形与待测区烃源岩TOC含量的关系，采用地震波形指示反演方法确定待测区的烃源TOC含量的反演数据体，解决了现有的采用间接反演波阻抗或者地震属性方法预测烃源岩TOC含量时存在的确定烃源TOC含量精度差，准确性低的技术问题以及多解性问题，提高了确定烃源岩TOC含量的精度和准确性；此外，本发明实施例还考虑了烃源岩TOC含量的非均质性，利用地层切片技术，通过求取不同层序级别下，分布较为稳定的最大湖泛面位置处的烃源岩TOC含量沿层地层切片处TOC含量的空间变化，解决现有方法由于未考虑烃源岩TOC含量的非均质性导致容易遗漏部分资源潜力较好的高TOC含量烃源岩层段的技术问题；另外本发明实施例还通过地质和实验数据、测井数据和地震反演相结合，实现了准确确定高TOC含量烃源岩分布范围和发育程度的技术效果；最后本发明实施例还通过实验-测井正演和测井-地震反演相结合的方法，根据实际情况充分地运用正演技术和反演技术的特点，有效地预测了烃源岩TOC含量的分布特征。
显然，本领域的技术人员应该明白，上述的本发明实施例的各模块或各步骤可以用通用的计算装置来实现，它们可以集中在单个的计算装置上，或者分布在多个计算装置所组成的网络上，可选地，它们可以用计算装置可执行的程序代码来实现，从而，可以将它们存储在存储装置中由计算装置来执行，并且在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤，或者将它们分别制作成各个集成电路模块，或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样，本发明实施例不限制于任何特定的硬件和软件结合。
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (10)
1.一种基于井震结合确定源岩TOC含量和非均质性方法，其特征在于，包括：
获取待测区的地质和实验数据；
根据所述地质和实验数据正演得到所述待测区的烃源岩TOC含量的特征曲线；
以所述待测区的层序格架和所述烃源岩TOC含量的特征曲线为约束，利用地震波形反演得到所述待测区的烃源岩TOC含量的反演数据体；
在所述待测区的层序格架下，根据所述烃源岩TOC含量的反演数据体，利用地层切片技术，分别确定所述待测区长期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图，中期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图和短期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图。
2.根据权利要求1所述的方法，其特征在于，所述获取待测区的地质和实验数据，包括：
对所述待测区内的钻井的岩心和岩屑进行系统采样，得到采样品；
根据所述采样品的TOC含量实验，得到所述采样品中烃源岩TOC含量数据；
以所述采样品的X衍射及能谱分析实验数据和所述钻井的测井数据为约束，从所述烃源岩TOC含量数据中剔除砂岩的TOC含量数据；
将所述采样品中剔除砂岩的TOC含量数据后剩余的烃源岩TOC含量数据作为所述待测区的地质和实验数据。
3.根据权利要求2所述的方法，其特征在于，对所述待测区内的钻井的岩心和岩屑进行系统采样的采样间距设置为1m到2m。
4.根据权利要求2所述的方法，其特征在于，根据所述地质和实验数据正演得到所述待测区的烃源岩TOC含量的特征曲线，包括：
根据所述采样品中剔除砂岩的TOC含量数据后剩余的烃源岩TOC含量数据，参考全岩X衍射和能谱分析实验，进行岩性精细归位；
结合所述钻井的测井曲线的纵向变化特征，确定所述待测区内的钻井纵向上烃源岩TOC含量的测井响应特征；
根据所述烃源岩TOC含量的测井响应特征，做交汇图分析，确定烃源岩TOC含量的敏感测井曲线；
根据所述敏感测井曲线和所述地质和实验数据，通过分段或者分层系确定烃源岩TOC含量的测井预测模型；
根据所述测井预测模型，正演得到所述烃源岩TOC含量的特征曲线。
5.根据权利要求1所述的方法，其特征在于，以所述待测区的层序格架和所述烃源岩TOC含量的特征曲线为约束，利用地震波形反演得到所述待测区的烃源岩TOC含量的反演数据体，包括：
在所述层序格架下，以所述烃源岩TOC含量特征曲线作为约束，根据所述待测区的三维地震资料，利用所述待测区的地震波形，建立反演模型；
根据所述反演模型，求解得到所述待测区的烃源岩TOC含量的反演数据体。
6.根据权利要求5所述的方法，其特征在于，利用所述待测区的地震波形，建立反演模型，包括：
根据所述层序格架内的地震解释层位和所述待测区内的钻井数量，确定样本探井数和高频成分；
根据所述样本探井数和所述高频成分，建立初始模型；
根据所述初始模型，反演得到所述待测区的烃源岩TOC含量的反演数据体；
将所述烃源岩TOC含量的反演数据体中连井反演剖面上的烃源岩TOC含量和基于所述钻井的测井数据得到的所述钻井的烃源岩TOC含量特征曲线进行比较，得到比较误差；
如果所述比较误差小于等于预设阈值，则将所述初始模型确定为所述反演模型。
7.根据权利要求6所述的方法，其特征在于，在将所述烃源岩TOC含量的反演数据体中连井反演剖面上的烃源岩TOC含量和基于所述钻井的测井数据得到的所述钻井的烃源岩TOC含量特征曲线进行比较，得到比较误差之后，还包括：
如果所述比较误差大于所述预设阈值，则重新确定样本探井数和高频成分，并根据重新确定的样本探井数和高频成分建立初始模型，直到根据建立的初始模型，通过地震波形反演得到的所述烃源岩TOC含量的反演数据体中连井反演剖面上的烃源岩TOC含量和基于所述钻井的测井数据得到的所述钻井的烃源岩TOC含量特征曲线之间的比较误差小于等于所述预设阈值，并将当前的初始模型作为所述反演模型。
8.根据权利要求1所述的方法，其特征在于，在所述待测区的层序格架下，根据所述烃源岩TOC含量的反演数据体，利用地层切片技术，分别确定所述待测区长期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图，中期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图和短期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图，包括：
在所述待测区的层序格架下，利用地层切片技术分别获取所述待测区所述长期旋回最大湖泛面位置的沿层地层切片、所述中期旋回最大湖泛面位置的沿层地层切片和所述短期旋回最大湖泛面位置的沿层地层切片；
利用所述烃源岩TOC含量的反演数据体，根据所述长期旋回最大湖泛面位置的沿层地层切片确定所述长期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图，根据所述中期旋回最大湖泛面位置的沿层地层切片确定所述中期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图，根据所述短期旋回最大湖泛面位置的沿层地层切片确定所述短期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图；
其中，所述长期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图用于表征所述长期旋回最大湖泛面位置的TOC含量的空间分布特征，所述中期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图用于表征所述中期旋回最大湖泛面位置的TOC含量的空间分布特征，所述短期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图用于表征所述短期旋回最大湖泛面位置的TOC含量的空间分布特征。
9.一种基于井震结合确定源岩TOC含量和非均质性装置，其特征在于，包括：
第一获取模块，用于获取待测区的地质和实验数据；
第二获取模块，用于根据所述地质和实验数据正演得到所述待测区的烃源岩TOC含量的特征曲线；
第三获取模块，用于以所述待测区的层序格架和所述烃源岩TOC含量的特征曲线为约束，利用地震波形反演得到所述待测区的烃源岩TOC含量的反演数据体；
确定模块，用于在所述待测区的层序格架下，根据所述烃源岩TOC含量的反演数据体，利用地层切片技术，分别确定所述待测区长期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图，中期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图和短期旋回最大湖泛面位置的烃源岩TOC含量等值线平面分布图。
10.根据权利要求9所述的装置，其特征在于，所述第三获取模块包括：
模型建立单元，用于在所述层序格架下，以所述烃源岩TOC含量特征曲线作为约束，根据所述待测区的三维地震资料，利用所述待测区的地震波形，建立反演模型；
数据体求解单元，用于根据所述反演模型，求解得到所述待测区的烃源岩TOC含量的反演数据体。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610803915.5A CN106249300B (zh) | 2016-09-05 | 2016-09-05 | 基于井震结合确定源岩toc含量和非均质性方法和装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610803915.5A CN106249300B (zh) | 2016-09-05 | 2016-09-05 | 基于井震结合确定源岩toc含量和非均质性方法和装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106249300A true CN106249300A (zh) | 2016-12-21 |
| CN106249300B CN106249300B (zh) | 2018-09-04 |
Family
ID=57599157
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610803915.5A Active CN106249300B (zh) | 2016-09-05 | 2016-09-05 | 基于井震结合确定源岩toc含量和非均质性方法和装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106249300B (zh) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106703798A (zh) * | 2016-12-31 | 2017-05-24 | 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 | 页岩非均质性划分与对比方法 |
| CN107218029A (zh) * | 2017-06-29 | 2017-09-29 | 中国石油天然气股份有限公司 | 一种平面非均质模型、注采模拟实验装置及方法 |
| CN110501744A (zh) * | 2019-08-26 | 2019-11-26 | 中海石油(中国)有限公司上海分公司 | 烃源岩有机碳地球物理定量预测方法、装置、设备及存储介质 |
| CN112394392A (zh) * | 2019-08-13 | 2021-02-23 | 中国石油天然气股份有限公司 | 对烃源岩分布情况进行评价的方法和装置 |
| CN114791634A (zh) * | 2021-01-25 | 2022-07-26 | 中国石油化工股份有限公司 | 压扭叠合盆地大中型油气田目标区综合评价方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2011112294A1 (en) * | 2010-03-11 | 2011-09-15 | Exxonmobil Upstream Research Company | Predicting anisotropic source rock properties from well data |
| US20120095687A1 (en) * | 2010-04-21 | 2012-04-19 | Baker Hughes Incorporated | Method of predicting source rock thermal maturity from log responses |
| CN104616588A (zh) * | 2015-01-21 | 2015-05-13 | 苏州市职业大学 | 一种小型手动液压材料拉伸仪 |
| CN105158796A (zh) * | 2015-09-16 | 2015-12-16 | 中国石油天然气股份有限公司 | 确定toc含量的方法和装置 |
-
2016
- 2016-09-05 CN CN201610803915.5A patent/CN106249300B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2011112294A1 (en) * | 2010-03-11 | 2011-09-15 | Exxonmobil Upstream Research Company | Predicting anisotropic source rock properties from well data |
| US20120095687A1 (en) * | 2010-04-21 | 2012-04-19 | Baker Hughes Incorporated | Method of predicting source rock thermal maturity from log responses |
| CN104616588A (zh) * | 2015-01-21 | 2015-05-13 | 苏州市职业大学 | 一种小型手动液压材料拉伸仪 |
| CN105158796A (zh) * | 2015-09-16 | 2015-12-16 | 中国石油天然气股份有限公司 | 确定toc含量的方法和装置 |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106703798A (zh) * | 2016-12-31 | 2017-05-24 | 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 | 页岩非均质性划分与对比方法 |
| CN106703798B (zh) * | 2016-12-31 | 2020-09-01 | 中国石油化工股份有限公司 | 页岩非均质性划分与对比方法 |
| CN107218029A (zh) * | 2017-06-29 | 2017-09-29 | 中国石油天然气股份有限公司 | 一种平面非均质模型、注采模拟实验装置及方法 |
| CN112394392A (zh) * | 2019-08-13 | 2021-02-23 | 中国石油天然气股份有限公司 | 对烃源岩分布情况进行评价的方法和装置 |
| CN112394392B (zh) * | 2019-08-13 | 2023-09-26 | 中国石油天然气股份有限公司 | 对烃源岩分布情况进行评价的方法和装置 |
| CN110501744A (zh) * | 2019-08-26 | 2019-11-26 | 中海石油(中国)有限公司上海分公司 | 烃源岩有机碳地球物理定量预测方法、装置、设备及存储介质 |
| CN110501744B (zh) * | 2019-08-26 | 2021-06-25 | 中海石油(中国)有限公司上海分公司 | 烃源岩有机碳地球物理定量预测方法、装置、设备及存储介质 |
| CN114791634A (zh) * | 2021-01-25 | 2022-07-26 | 中国石油化工股份有限公司 | 压扭叠合盆地大中型油气田目标区综合评价方法 |
| CN114791634B (zh) * | 2021-01-25 | 2023-05-09 | 中国石油化工股份有限公司 | 压扭叠合盆地大中型油气田目标区综合评价方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106249300B (zh) | 2018-09-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104977618B (zh) | 一种评价页岩气储层及寻找甜点区的方法 | |
| CN103454685B (zh) | 利用测井约束波阻抗反演预测砂体厚度的方法和装置 | |
| CN106951660A (zh) | 一种海相碎屑岩水平井储层测井解释方法及装置 | |
| CN106249300B (zh) | 基于井震结合确定源岩toc含量和非均质性方法和装置 | |
| CN105675635B (zh) | 致密岩石组分相对含量及脆性指数确定方法和装置 | |
| Williams-Stroud et al. | Microseismicity-constrained discrete fracture network models for stimulated reservoir simulation | |
| Caers et al. | Stochastic integration of seismic data and geologic scenarios: A West Africa submarine channel saga | |
| CN105629325B (zh) | 前陆盆地冲积扇精细刻画与预测方法 | |
| CN104632202B (zh) | 确定干粘土三孔隙度测井参数值的方法及装置 | |
| US20150212224A1 (en) | Singularity spectrum analysis of microseismic data | |
| CN104834007B (zh) | 地震反演过程中计算碳酸盐岩缝洞型储层充填程度的方法 | |
| CN107515430B (zh) | 一种地震法探测盐湖卤水的方法 | |
| CN107229076A (zh) | 一种基于测井资料进行温度响应特征分析的方法 | |
| Duan et al. | Architectural characterization of Ordovician fault-controlled paleokarst carbonate reservoirs, Tahe oilfield, China | |
| Welch | Using geomechanical models to simulate the growth of the fracture network in the Ekofisk Formation of the Kraka structure, Danish Central Graben | |
| CN112528106A (zh) | 火山岩岩性识别方法 | |
| Zhang et al. | Architecture characteristics and characterization methods of fault-controlled karst reservoirs: A case study of the Shunbei 5 fault zone in the Tarim Basin, China | |
| Talley et al. | Dynamic reservoir characterization of Vacuum Field | |
| Gritto et al. | Seismic imaging of reservoir heterogeneity using a network with high station density at The Geysers geothermal reservoir, CA, USA | |
| Neely et al. | Seismic-based characterization of reservoir heterogeneity within the Meramec interval of the STACK play, Central Oklahoma | |
| Abbaszadeh et al. | Development of an integrated reservoir model for a naturally fractured volcanic reservoir in China | |
| Shi et al. | Identification and distribution of fractures in the Zhangjiatan shale of the Mesozoic Yanchang Formation in Ordos Basin | |
| Aydin et al. | Characterization of Geothermal Reservoirs | |
| Li et al. | Application of brittleness index to interpret microseismic event distribution in a hydraulically fractured shale formation | |
| Masoud et al. | Reservoir characterization and geostatistical model of the Cretaceous and Cambrian-Ordovician reservoir intervals, Meghil field, Sirte basin, Libya |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |