CN105629310A - 碳酸盐岩储层的无井约束地质统计学反演的方法和装置 - Google Patents
碳酸盐岩储层的无井约束地质统计学反演的方法和装置 Download PDFInfo
- Publication number
- CN105629310A CN105629310A CN201410637165.XA CN201410637165A CN105629310A CN 105629310 A CN105629310 A CN 105629310A CN 201410637165 A CN201410637165 A CN 201410637165A CN 105629310 A CN105629310 A CN 105629310A
- Authority
- CN
- China
- Prior art keywords
- well
- inversion
- carbonate reservoir
- predicted
- seismic data
- 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
- 238000000034 method Methods 0.000 title claims abstract description 103
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 title claims abstract description 84
- 238000005553 drilling Methods 0.000 claims description 43
- 238000012937 correction Methods 0.000 claims description 8
- 238000001615 p wave Methods 0.000 claims description 8
- 239000000284 extract Substances 0.000 claims description 7
- 238000013508 migration Methods 0.000 claims description 7
- 230000005012 migration Effects 0.000 claims description 7
- 238000007689 inspection Methods 0.000 claims description 6
- 239000011435 rock Substances 0.000 abstract description 10
- 238000005457 optimization Methods 0.000 abstract description 9
- 238000012545 processing Methods 0.000 abstract description 3
- 238000005315 distribution function Methods 0.000 abstract 2
- 150000004649 carbonic acid derivatives Chemical class 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 238000003908 quality control method Methods 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000008034 disappearance Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供一种碳酸盐岩储层的无井约束地质统计学反演的方法和装置，包括：根据预先获取的待预测碳酸盐岩储层的原始测井资料统计获取待预测碳酸盐岩储层的概率密度分布函数；根据待预测碳酸盐岩储层的原始测井资料和原始地震资料进行井震精细标定处理，获取多井优化子波；根据概率密度分布函数、多井优化子波、预先获取的变差函数以及底层格架进行有井约束的地质统计学反演处理，并增加钻井作为盲井重复进行本步骤，直至所有钻井全部充当盲井时，获取第一反演结果；通过岩石物理建模正演、有井约束地质统计学反演参数优选，进行无井约束地质统计学反演获取高分辨率反演结果，提高对碳酸盐岩储层预测结果的精确性。
Description
技术领域
本发明涉及石油勘探技术，尤其涉及一种碳酸盐岩储层的无井约束地质统计学反演的方法和装置。
背景技术
对于碳酸盐岩储层的预测一直是勘探行业的一大难题，目前常用的反演预测方法包括地质统计学反演法、以波阻抗体数据体等方式的预测方式。
现有技术中，最长使用的预测方法是采用地质统计学的方式进行反演预测，地质统计学在兼顾地震横向分辨率和测井纵向分辨率的基础上，综合了地震反演与储层随机建模的优势。与常规反演手段相比，地质统计学反演具有垂向分辨率高，非均质性表征效果更好的特点。然而，对于碳酸盐岩等非均质性强的地层，测井仅反映井眼附近的地层信息，且钻井过程中容易发生钻井液漏失和放空，造成测井曲线缺失或失真。
由于碳酸盐岩储层非均质性强，储层横向变化大，增加了常规地质统计学反演进行储层预测的难度，导致对碳酸盐岩储层预测结果的精确性较低。
发明内容
本发明提供一种碳酸盐岩储层的无井约束地质统计学反演的方法和装置，用于解决现有技术中由于碳酸盐岩储层非均质性强，储层横向变化大，增加了常规地质统计学反演进行储层预测的难度，导致对碳酸盐岩储层预测结果的精确性较低的问题。
本发明一方面提供一种碳酸盐岩储层的无井约束地质统计学反演的方法，包括：
根据预先获取的待预测碳酸盐岩储层的原始测井资料统计获取所述待预测碳酸盐岩储层的概率密度分布函数；
根据所述待预测碳酸盐岩储层的所述原始测井资料和原始地震资料进行井震精细标定处理，获取多井优化子波；
根据所述概率密度分布函数、所述多井优化子波、预先获取的变差函数以及底层格架进行有井约束的地质统计学反演处理，并增加钻井作为盲井重复进行本步骤，直至所有钻井全部充当盲井时，获取第一反演结果；其中，所述地层格架为预先根据所述地震资料建立的。
可选的，所述根据所述概率密度分布函数、所述多井优化子波、预先获取的变差函数以及底层格架进行有井约束的地质统计学反演处理，并增加钻井作为盲井重复进行本步骤，直至所有钻井全部充当盲井时，获取第一反演结果，包括：
根据所述概率密度分布函数、所述多井优化子波、预先获取的变差函数以及底层格架进行有井约束的地质统计学反演处理，获取反演参数；
根据预先设置的反演结果阈值判断所述反演参数对应的反演结果是否在所述反演结果阈值范围内，若所述反演参数对应的反演结果不在所述反演参数阈值范围内，则调整所述反演参数直至所述反演结果在所述反演结果阈值范围内，则将此时的反演参数作为第一反演参数；
根据所述第一反演参数，并增加钻井作为盲井重复进行有井约束的地质统计学反演处理，直至所有钻井全部充当盲井时，实现无井约束的地质统计学反演处理，并获取第一反演结果。
可选的，所述根据预先获取的待预测碳酸盐岩储层的原始测井资料统计获取所述待预测碳酸盐岩储层的概率密度分布函数，包括：
对所述待预测碳酸盐岩储层的所述原始测井资料进行深度矫正、井眼垮塌矫正和多井一致性处理获取处理后的测井资料；
根据所述处理后的测井资料，建立岩石物理模型，统计获取所述概率密度分布函数。
可选的，所述根据所述待预测碳酸盐岩储层的所述原始测井资料和原始地震资料进行井震精细标定处理，获取多井优化子波，包括：
对所述原始地震资料进行合理性处理，获取合理的地震资料；
对所述合理的地震资料进行层位精细解释，并建立合理的地层格架；
根据所述合理的地震资料和所述地层格架，分析各层系波组特征，提取多井优化子波；
其中，所述合理性处理包括：对所述地震资料的信噪比、分辨率、能量增益和采集脚印的检验、去噪处理和偏移归位处理。
可选的，所述反演结果包括纵波阻抗和孔隙度。
本发明另一方面提供一种碳酸盐岩储层的无井约束地质统计学反演的装置，包括：
统计模块，用于根据预先获取的待预测碳酸盐岩储层的原始测井资料统计获取所述待预测碳酸盐岩储层的概率密度分布函数；
标定模块，用于根据所述待预测碳酸盐岩储层的所述原始测井资料和原始地震资料进行井震精细标定处理，获取多井优化子波；
反演模块，用于根据所述概率密度分布函数、所述多井优化子波、预先获取的变差函数以及底层格架进行有井约束的地质统计学反演处理，并增加钻井作为盲井重复进行本步骤，直至所有钻井全部充当盲井时，获取第一反演结果；其中，所述地层格架为预先根据所述地震资料建立的。
可选的，所述反演模块具体用于：
根据所述概率密度分布函数、所述多井优化子波、预先获取的变差函数以及底层格架进行有井约束的地质统计学反演处理，获取反演参数；
根据预先设置的反演结果阈值判断所述反演参数对应的反演结果是否在所述反演结果阈值范围内，若所述反演参数对应的反演结果不在所述反演参数阈值范围内，则调整所述反演参数直至所述反演结果在所述反演结果阈值范围内，则将此时的反演参数作为第一反演参数；
根据所述第一反演参数，并增加钻井作为盲井重复进行有井约束的地质统计学反演处理，直至所有钻井全部充当盲井时，实现无井约束的地质统计学反演处理，并获取第一反演结果。
可选的，所述统计模块具体用于：
对所述待预测碳酸盐岩储层的所述原始测井资料进行深度矫正、井眼垮塌矫正和多井一致性处理获取处理后的测井资料；
根据所述处理后的测井资料，建立岩石物理模型，统计获取所述概率密度分布函数。
可选的，所述标定模块具体用于：
对所述原始地震资料进行合理性处理，获取合理的地震资料；
对所述合理的地震资料进行层位精细解释，并建立合理的地层格架；
根据所述合理的地震资料和所述地层格架，分析各层系波组特征，提取多井优化子波；
其中，所述合理性处理包括：对所述地震资料的信噪比、分辨率、能量增益和采集脚印的检验、去噪处理和偏移归位处理。
可选的，所述反演模块获取的所述反演结果包括纵波阻抗和孔隙度。
本发明提供的碳酸盐岩储层的无井约束地质统计学反演的方法和装置，原始测井资料统计获取待预测碳酸盐岩储层的概率密度分布函数，并根据原始地震资料获取多井优化子波，再根据概率密度分布函数、多井优化子波、预先获取的变差函数以及底层格架进行有井约束的地质统计学反演处理，并增加钻井作为盲井重复进行直至所有钻井全部充当盲井时，获取第一反演结果，通过对岩石物理建模正演、有井约束地质统计学反演参数优选，进行无井约束地质统计学反演获取高分辨率反演结果，提高对碳酸盐岩储层预测结果的精确性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明碳酸盐岩储层的无井约束地质统计学反演的方法实施例一的流程图；
图2为本发明碳酸盐岩储层的无井约束地质统计学反演的方法实施例二的流程图；
图3为本发明碳酸盐岩储层的无井约束地质统计学反演的方法实施例三的流程图；
图4为本发明碳酸盐岩储层的无井约束地质统计学反演的方法实施例四的流程图；
图5为本发明碳酸盐岩储层的无井约束地质统计学反演的施例流程示意图；
图6为本发明碳酸盐岩储层的无井约束地质统计学反演的装置实施例一的结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
图1为本发明碳酸盐岩储层的无井约束地质统计学反演的方法实施例一的流程图，本实施例提供的无井约束地质统计学反演的方法用于对碳酸盐岩储层等非均质性强的地层进行反演，提供更可靠的反演结果，该方法的具体步骤如下：
S101：根据预先获取的待预测碳酸盐岩储层的原始测井资料统计获取所述待预测碳酸盐岩储层的概率密度分布函数。
在本实施例中，在进行地质勘探时，首先会获取待预测位置的碳酸盐岩储层的测井资料，一般情况下测井资料至少包括电阻率、岩性和孔隙度的所有特征，测井资料记录的一般都是物理资料如电阻率、自然电位、声波速度和岩石体积密度等。
S102：根据所述待预测碳酸盐岩储层的所述原始测井资料和原始地震资料进行井震精细标定处理，获取多井优化子波。
在本实施例中，原始地震资料一般通过地震勘探仪等设备在野外采集记录，或者人为设计实验获取。进行精细标定的方式主要是进行层位的标定，利用VSP资料，人工记录合成资料等建立井层关系，使钻井的地质层位和地震反射层一一对应。
S103：根据所述概率密度分布函数、所述多井优化子波、预先获取的变差函数以及底层格架进行有井约束的地质统计学反演处理，并增加钻井作为盲井重复进行本步骤，直至所有钻井全部充当盲井时，获取第一反演结果；其中，所述地层格架为预先根据所述地震资料建立的。
在本实施例中，变差函数主要是通过对地质资料的分析，结合测井统计规律来获取的。
本实施例提供的碳酸盐岩储层的无井约束地质统计学反演的方法，原始测井资料统计获取待预测碳酸盐岩储层的概率密度分布函数，并根据原始地震资料获取多井优化子波，再根据概率密度分布函数、多井优化子波、预先获取的变差函数以及底层格架进行有井约束的地质统计学反演处理，并增加钻井作为盲井重复进行直至所有钻井全部充当盲井时，获取第一反演结果，通过对岩石物理建模正演、有井约束地质统计学反演参数优选，进行无井约束地质统计学反演获取高分辨率反演结果，提高对碳酸盐岩储层预测结果的精确性。
图2为本发明碳酸盐岩储层的无井约束地质统计学反演的方法实施例二的流程图，在上述实施例的基础上，如图2所示，S103中的根据所述概率密度分布函数、所述多井优化子波、预先获取的变差函数以及底层格架进行有井约束的地质统计学反演处理，并增加钻井作为盲井重复进行本步骤，直至所有钻井全部充当盲井时，获取第一反演结果，其具体的实现步骤包括：
S201：根据所述概率密度分布函数、所述多井优化子波、预先获取的变差函数以及底层格架进行有井约束的地质统计学反演处理，获取反演参数。
S202：根据预先设置的反演结果阈值判断所述反演参数对应的反演结果是否在所述反演结果阈值范围内，若所述反演参数对应的反演结果不在所述反演参数阈值范围内，则调整所述反演参数直至所述反演结果在所述反演结果阈值范围内，则将此时的反演参数作为第一反演参数。
在本实施例中，利用上面步骤得到的概率密度函数、变差函数、多井优化子波和地层格架进行有井约束的地质统计学反演，并不断调整反演参数，得到合理的地质统计学纵波阻抗和孔隙度等反演结果。
S203：根据所述第一反演参数，并增加钻井作为盲井重复进行有井约束的地质统计学反演处理，直至所有钻井全部充当盲井时，实现无井约束的地质统计学反演处理，并获取第一反演结果。
在本实施例中，根据上一步骤中有井约束的地质统计学反演结果指导选取优化反演参数，再逐步增加盲井，最终以全部钻井充当盲井，最终实现无井约束的地质统计学反演，得到优选的纵波阻抗和孔隙度等第一反演结果。
图3为本发明碳酸盐岩储层的无井约束地质统计学反演的方法实施例三的流程图，在上述实施例一和实施例二的基础上，如图3所示，S101中的根据预先获取的待预测碳酸盐岩储层的原始测井资料统计获取所述待预测碳酸盐岩储层的概率密度分布函数，其具体的实现步骤包括：
S301：对所述待预测碳酸盐岩储层的所述原始测井资料进行深度矫正、井眼垮塌矫正和多井一致性处理获取处理后的测井资料。
S302：根据所述处理后的测井资料，建立岩石物理模型，统计获取所述概率密度分布函数。
在本实施例中，首先对原始测井资料进行深度校正、井眼垮塌校正、多井一致性等处理，在此基础上，通过建立合理的岩石物理模型正演漏失和放空的井段，统计该碳酸盐岩储层的概率密度分布函数。
图4为本发明碳酸盐岩储层的无井约束地质统计学反演的方法实施例四的流程图，在上述实施例一至实施例三的基础上，如图4所示，S102中的根据所述待预测碳酸盐岩储层的所述原始测井资料和原始地震资料进行井震精细标定处理，获取多井优化子波，其具体实现步骤包括：
S401：对所述原始地震资料进行合理性处理，获取合理的地震资料。
S402：对所述合理的地震资料进行层位精细解释，并建立合理的地层格架。
S403：根据所述合理的地震资料和所述地层格架，分析各层系波组特征，提取多井优化子波。
其中，所述合理性处理包括：对所述地震资料的信噪比、分辨率、能量增益和采集脚印的检验、去噪处理和偏移归位处理。
在本实施例中，对地震资料进行信噪比、分辨率、去噪去多次、偏移归位、能量增益及采集脚印等合理性检查，并进行层位精细解释，建立合理的地层格架；在对上述测井地震等资料严格质控的基础上，进行井震一致性研究及精细标定，分析各层系波组特征，提取多井优化子波。
特别的在上述任一实施例中，所述反演结果至少包括纵波阻抗和孔隙度，还可以包括其他常用的特征。
实施例二至实施例四提供的碳酸盐岩储层的无井约束地质统计学反演的方法，针对碳酸盐储层非均质性强等特点，通过岩石物理建模正演、有井约束地质统计学反演参数优选，进行无井约束地质统计学反演，并得到了高分辨率反演结果。与有井约束的地质统计学反演结果相比，本发明的无井约束的地质统计学反演消除了井震不匹配和碳酸盐岩非均值性强难以预测的缺点，更加接近地下实际情况，预测精度提高。
图5为本发明碳酸盐岩储层的无井约束地质统计学反演的施例流程示意图，在上述任一实施例的基础上，如图5所示，本实施例提供本方法在实际应用过程中的处理过程示意。具体的，
根据原始测井数据获取测井曲线，根据该测井曲线对多井测井进行评价并建立岩石物理模型，再根据该模型统计出概率密度函数；
对预先获取的原始地震资料进行合理性处理实现严格的质控，对合理性处理后的地震资料进行层位精细解释，并在此基础上建立底层格架；
根据待评价的碳酸盐岩储层的地址资料获取变差函数；
上述三个步骤为本发明碳酸盐岩储层的无井约束地质统计学反演的方法的基础步骤，可以并列进行也可以依次进行没有严格的顺序。
根据获取到的岩石物理建模和合理的地震资料，进行井震一致性研究及精细标定，分析各层系波组特征，提取多井优化子波。
在上述所有获取的结果基础上对该碳酸盐岩储层进行有井约束地质统计学反演，并不断的优化反演参数，获取合理的优化参数后按照该方式，并增加钻井作为盲井重复进行有井约束的地质统计学反演处理，直至所有钻井全部充当盲井时，实现无井约束的地质统计学反演处理，并获取最后的反演结果，即第一反演结果。
本实施例提供的碳酸盐岩储层的无井约束地质统计学反演的方法，与有井约束的地质统计学反演结果相比，无井约束的地质统计学反演消除了井震不匹配和碳酸盐岩非均值性强难以预测的缺点，更加接近地下实际情况，预测精度提高。
图6为本发明碳酸盐岩储层的无井约束地质统计学反演的装置实施例一的结构示意图，如图6所示，该碳酸盐岩储层的无井约束地质统计学反演的装置60，包括：统计模块61、标定模块62和反演模块63。
统计模块61，用于根据预先获取的待预测碳酸盐岩储层的原始测井资料统计获取所述待预测碳酸盐岩储层的概率密度分布函数；
标定模块62，用于根据所述待预测碳酸盐岩储层的所述原始测井资料和原始地震资料进行井震精细标定处理，获取多井优化子波；
反演模块63，用于根据所述概率密度分布函数、所述多井优化子波、预先获取的变差函数以及底层格架进行有井约束的地质统计学反演处理，并增加钻井作为盲井重复进行本步骤，直至所有钻井全部充当盲井时，获取第一反演结果；其中，所述地层格架为预先根据所述地震资料建立的。
本实施例提供的碳酸盐岩储层的无井约束地质统计学反演的装置，用于执行图1-5所示的方法实施例的技术方案，其实现原理和技术效果类似，统计模块根据原始测井资料统计获取待预测碳酸盐岩储层的概率密度分布函数，标定模块根据原始地震资料获取多井优化子波，反演模块根据概率密度分布函数、多井优化子波、预先获取的变差函数以及底层格架进行有井约束的地质统计学反演处理，并增加钻井作为盲井重复进行直至所有钻井全部充当盲井时，获取第一反演结果，通过对岩石物理建模正演、有井约束地质统计学反演参数优选，进行无井约束地质统计学反演获取高分辨率反演结果，提高对碳酸盐岩储层预测结果的精确性。
在本发明碳酸盐岩储层的无井约束地质统计学反演的装置实施例二中，在上述装置实施例一的基础上，所述反演模块63具体用于：
根据所述概率密度分布函数、所述多井优化子波、预先获取的变差函数以及底层格架进行有井约束的地质统计学反演处理，获取反演参数；
根据预先设置的反演结果阈值判断所述反演参数对应的反演结果是否在所述反演结果阈值范围内，若所述反演参数对应的反演结果不在所述反演参数阈值范围内，则调整所述反演参数直至所述反演结果在所述反演结果阈值范围内，则将此时的反演参数作为第一反演参数；
根据所述第一反演参数，并增加钻井作为盲井重复进行有井约束的地质统计学反演处理，直至所有钻井全部充当盲井时，实现无井约束的地质统计学反演处理，并获取第一反演结果。
可选的，所述统计模块61具体用于：
对所述待预测碳酸盐岩储层的所述原始测井资料进行深度矫正、井眼垮塌矫正和多井一致性处理获取处理后的测井资料；
根据所述处理后的测井资料，建立岩石物理模型，统计获取所述概率密度分布函数。
可选的，所述标定模块62具体用于：
对所述原始地震资料进行合理性处理，获取合理的地震资料；
对所述合理的地震资料进行层位精细解释，并建立合理的地层格架；
其中，所述合理性处理包括：对所述地震资料的信噪比、分辨率、能量增益和采集脚印的检验、去噪处理和偏移归位处理。
可选的，所述反演模块63获取的所述反演结果包括纵波阻抗和孔隙度。
本实施例提供的碳酸盐岩储层的无井约束地质统计学反演的装置，用于执行图1-5所示的方法实施例的技术方案，其实现原理和技术效果类似，在此不再赘述。
最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。
Claims (10)
1.一种碳酸盐岩储层的无井约束地质统计学反演的方法，其特征在于，包括：
根据预先获取的待预测碳酸盐岩储层的原始测井资料统计获取所述待预测碳酸盐岩储层的概率密度分布函数；
根据所述待预测碳酸盐岩储层的所述原始测井资料和原始地震资料进行井震精细标定处理，获取多井优化子波；
根据所述概率密度分布函数、所述多井优化子波、预先获取的变差函数以及底层格架进行有井约束的地质统计学反演处理，并增加钻井作为盲井重复进行本步骤，直至所有钻井全部充当盲井时，获取第一反演结果；其中，所述地层格架为预先根据所述地震资料建立的。
2.根据权利要求1所述的方法，其特征在于，所述根据所述概率密度分布函数、所述多井优化子波、预先获取的变差函数以及底层格架进行有井约束的地质统计学反演处理，并增加钻井作为盲井重复进行本步骤，直至所有钻井全部充当盲井时，获取第一反演结果，包括：
根据所述概率密度分布函数、所述多井优化子波、预先获取的变差函数以及底层格架进行有井约束的地质统计学反演处理，获取反演参数；
根据预先设置的反演结果阈值判断所述反演参数对应的反演结果是否在所述反演结果阈值范围内，若所述反演参数对应的反演结果不在所述反演参数阈值范围内，则调整所述反演参数直至所述反演结果在所述反演结果阈值范围内，则将此时的反演参数作为第一反演参数；
根据所述第一反演参数，并增加钻井作为盲井重复进行有井约束的地质统计学反演处理，直至所有钻井全部充当盲井时，实现无井约束的地质统计学反演处理，并获取第一反演结果。
3.根据权利要求2所述的方法，其特征在于，所述根据预先获取的待预测碳酸盐岩储层的原始测井资料统计获取所述待预测碳酸盐岩储层的概率密度分布函数，包括：
对所述待预测碳酸盐岩储层的所述原始测井资料进行深度矫正、井眼垮塌矫正和多井一致性处理获取处理后的测井资料；
根据所述处理后的测井资料，建立岩石物理模型，统计获取所述概率密度分布函数。
4.根据权利要求3所述的方法，其特征在于，所述根据所述待预测碳酸盐岩储层的所述原始测井资料和原始地震资料进行井震精细标定处理，获取多井优化子波，包括：
对所述原始地震资料进行合理性处理，获取合理的地震资料；
对所述合理的地震资料进行层位精细解释，并建立合理的地层格架；
根据所述合理的地震资料和所述地层格架，分析各层系波组特征，提取多井优化子波；
其中，所述合理性处理包括：对所述地震资料的信噪比、分辨率、能量增益和采集脚印的检验、去噪处理和偏移归位处理。
5.根据权利要求1至4任一项所述的方法，其特征在于，所述反演结果包括纵波阻抗和孔隙度。
6.一种碳酸盐岩储层的无井约束地质统计学反演的装置，其特征在于，包括：
统计模块，用于根据预先获取的待预测碳酸盐岩储层的原始测井资料统计获取所述待预测碳酸盐岩储层的概率密度分布函数；
标定模块，用于根据所述待预测碳酸盐岩储层的所述原始测井资料和原始地震资料进行井震精细标定处理，获取多井优化子波；
反演模块，用于根据所述概率密度分布函数、所述多井优化子波、预先获取的变差函数以及底层格架进行有井约束的地质统计学反演处理，并增加钻井作为盲井重复进行本步骤，直至所有钻井全部充当盲井时，获取第一反演结果；其中，所述地层格架为预先根据所述地震资料建立的。
7.根据权利要求6所述的装置，其特征在于，所述反演模块具体用于：
根据所述概率密度分布函数、所述多井优化子波、预先获取的变差函数以及底层格架进行有井约束的地质统计学反演处理，获取反演参数；
根据预先设置的反演结果阈值判断所述反演参数对应的反演结果是否在所述反演结果阈值范围内，若所述反演参数对应的反演结果不在所述反演参数阈值范围内，则调整所述反演参数直至所述反演结果在所述反演结果阈值范围内，则将此时的反演参数作为第一反演参数；
根据所述第一反演参数，并增加钻井作为盲井重复进行有井约束的地质统计学反演处理，直至所有钻井全部充当盲井时，实现无井约束的地质统计学反演处理，并获取第一反演结果。
8.根据权利要求7所述的装置，其特征在于，所述统计模块具体用于：
对所述待预测碳酸盐岩储层的所述原始测井资料进行深度矫正、井眼垮塌矫正和多井一致性处理获取处理后的测井资料；
根据所述处理后的测井资料，建立岩石物理模型，统计获取所述概率密度分布函数。
9.根据权利要求8所述的装置，其特征在于，所述标定模块具体用于：
对所述原始地震资料进行合理性处理，获取合理的地震资料；
对所述合理的地震资料进行层位精细解释，并建立合理的地层格架；
根据所述合理的地震资料和所述地层格架，分析各层系波组特征，提取多井优化子波；
其中，所述合理性处理包括：对所述地震资料的信噪比、分辨率、能量增益和采集脚印的检验、去噪处理和偏移归位处理。
10.根据权利要求6至9任一项所述的装置，其特征在于，所述反演模块获取的所述反演结果包括纵波阻抗和孔隙度。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410637165.XA CN105629310B (zh) | 2014-11-05 | 2014-11-05 | 碳酸盐岩储层的无井约束地质统计学反演的方法和装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410637165.XA CN105629310B (zh) | 2014-11-05 | 2014-11-05 | 碳酸盐岩储层的无井约束地质统计学反演的方法和装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105629310A true CN105629310A (zh) | 2016-06-01 |
| CN105629310B CN105629310B (zh) | 2018-07-13 |
Family
ID=56044421
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410637165.XA Active CN105629310B (zh) | 2014-11-05 | 2014-11-05 | 碳酸盐岩储层的无井约束地质统计学反演的方法和装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105629310B (zh) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108957529A (zh) * | 2018-04-24 | 2018-12-07 | 中国海洋石油集团有限公司 | 一种基于属性的无井区子波估计方法 |
| CN111983690A (zh) * | 2020-08-12 | 2020-11-24 | 阳泉煤业（集团）股份有限公司 | 一种煤层顶板砂岩分布地球物理预测方法 |
| CN112363244A (zh) * | 2020-10-26 | 2021-02-12 | 中国石油天然气集团有限公司 | 波阻抗反演方法与碳酸盐岩非均质储层预测方法及系统 |
| CN112904441A (zh) * | 2021-01-13 | 2021-06-04 | 中国海洋石油集团有限公司 | 一种湖相碳酸盐岩孔隙度的计算方法 |
| CN113050157A (zh) * | 2020-10-15 | 2021-06-29 | 中国石油天然气股份有限公司 | 一种基于露头资料的碳酸盐岩地震储层反演方法及系统 |
| CN113325468A (zh) * | 2020-02-28 | 2021-08-31 | 中国石油天然气集团有限公司 | 储层缝洞分布范围预测方法及装置 |
| CN113496070A (zh) * | 2020-04-01 | 2021-10-12 | 中国石油天然气股份有限公司 | 地层俘获截面曲线的处理方法、装置、设备及介质 |
| CN114488305A (zh) * | 2022-02-16 | 2022-05-13 | 重庆科技学院 | 一种无井新探区地震资料地质层位精细标定方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB2361336A (en) * | 1999-09-29 | 2001-10-17 | Nec Corp | Degree of outlier calculation device, and probability density estimation device and histogram calculation device for use therein |
| US8090555B2 (en) * | 2007-06-21 | 2012-01-03 | Schlumberger Technology Corporation | Multi-attribute seismic characterization of gas hydrates |
| CN103678899A (zh) * | 2013-12-05 | 2014-03-26 | 中国石油天然气股份有限公司 | 一种基于可变影响比的油气藏多点统计建模方法及设备 |
| CN103713319A (zh) * | 2013-12-31 | 2014-04-09 | 张远银 | 一种基于地震约束建模的叠前反演方法 |
| CN104090303A (zh) * | 2014-07-10 | 2014-10-08 | 中国海洋石油总公司 | 一种地震反演方法和装置 |
-
2014
- 2014-11-05 CN CN201410637165.XA patent/CN105629310B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB2361336A (en) * | 1999-09-29 | 2001-10-17 | Nec Corp | Degree of outlier calculation device, and probability density estimation device and histogram calculation device for use therein |
| US8090555B2 (en) * | 2007-06-21 | 2012-01-03 | Schlumberger Technology Corporation | Multi-attribute seismic characterization of gas hydrates |
| CN103678899A (zh) * | 2013-12-05 | 2014-03-26 | 中国石油天然气股份有限公司 | 一种基于可变影响比的油气藏多点统计建模方法及设备 |
| CN103713319A (zh) * | 2013-12-31 | 2014-04-09 | 张远银 | 一种基于地震约束建模的叠前反演方法 |
| CN104090303A (zh) * | 2014-07-10 | 2014-10-08 | 中国海洋石油总公司 | 一种地震反演方法和装置 |
Non-Patent Citations (2)
| Title |
|---|
| 李国会 等: "塔北隆起英买2区块缝洞型碳酸盐岩等效孔隙度模型定量描述研究", 《中国石油勘探》 * |
| 王香文 等: "地质统计学反演技术在薄储层预测中的应用", 《石油与天然气地质》 * |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108957529B (zh) * | 2018-04-24 | 2020-01-24 | 中国海洋石油集团有限公司 | 一种基于属性的无井区子波估计方法 |
| CN108957529A (zh) * | 2018-04-24 | 2018-12-07 | 中国海洋石油集团有限公司 | 一种基于属性的无井区子波估计方法 |
| CN113325468A (zh) * | 2020-02-28 | 2021-08-31 | 中国石油天然气集团有限公司 | 储层缝洞分布范围预测方法及装置 |
| CN113496070B (zh) * | 2020-04-01 | 2024-03-01 | 中国石油天然气股份有限公司 | 地层俘获截面曲线的处理方法、装置、设备及介质 |
| CN113496070A (zh) * | 2020-04-01 | 2021-10-12 | 中国石油天然气股份有限公司 | 地层俘获截面曲线的处理方法、装置、设备及介质 |
| CN111983690A (zh) * | 2020-08-12 | 2020-11-24 | 阳泉煤业（集团）股份有限公司 | 一种煤层顶板砂岩分布地球物理预测方法 |
| CN113050157B (zh) * | 2020-10-15 | 2024-03-26 | 中国石油天然气股份有限公司 | 一种基于露头资料的碳酸盐岩地震储层反演方法及系统 |
| CN113050157A (zh) * | 2020-10-15 | 2021-06-29 | 中国石油天然气股份有限公司 | 一种基于露头资料的碳酸盐岩地震储层反演方法及系统 |
| CN112363244A (zh) * | 2020-10-26 | 2021-02-12 | 中国石油天然气集团有限公司 | 波阻抗反演方法与碳酸盐岩非均质储层预测方法及系统 |
| CN112363244B (zh) * | 2020-10-26 | 2024-05-28 | 中国石油天然气集团有限公司 | 波阻抗反演方法与碳酸盐岩非均质储层预测方法及系统 |
| CN112904441B (zh) * | 2021-01-13 | 2022-04-15 | 中国海洋石油集团有限公司 | 一种湖相碳酸盐岩孔隙度的计算方法 |
| CN112904441A (zh) * | 2021-01-13 | 2021-06-04 | 中国海洋石油集团有限公司 | 一种湖相碳酸盐岩孔隙度的计算方法 |
| CN114488305A (zh) * | 2022-02-16 | 2022-05-13 | 重庆科技学院 | 一种无井新探区地震资料地质层位精细标定方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105629310B (zh) | 2018-07-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105629310A (zh) | 碳酸盐岩储层的无井约束地质统计学反演的方法和装置 | |
| CN104047597B (zh) | 含油气泥页岩地层测井曲线标准化方法 | |
| US20170205531A1 (en) | Geological modeling workflow | |
| CN104200115B (zh) | 一种基于地质统计学模拟的全地层速度建模方法 | |
| Senel et al. | Integrated reservoir modeling at the Illinois Basin–Decatur project | |
| CN105527653B (zh) | 一种基于地质信息的虚拟井构建方法 | |
| CN103867194A (zh) | 一种砂体结构的测井表征方法与钻井层段选择方法及装置 | |
| Soleimani et al. | 3D static reservoir modeling by geostatistical techniques used for reservoir characterization and data integration | |
| US20220155480A1 (en) | Secondary recovery surveillance using validated streamline-based simulation | |
| CN105093293A (zh) | 低频补偿提高洞穴型碳酸盐岩储层地震量化预测的方法 | |
| CN104914470A (zh) | 碳酸盐岩缝洞储层储量校正方法 | |
| CN106285642A (zh) | 一种基于地震资料的地层孔隙压力预测方法 | |
| Moore et al. | Tight gas sandstone reservoirs, part 2: petrophysical analysis and reservoir modeling | |
| CN105301647A (zh) | 区分灰质泥岩和砂岩的方法 | |
| CN105589098B (zh) | 碳酸盐岩去除沉积泥质横向影响的储层反演方法和系统 | |
| US11474272B2 (en) | Methods and systems for identifying and plugging subterranean conduits | |
| Jafari et al. | Estimation of acoustic-impedance model by using model-based seismic inversion on the Ghar Member of Asmari Formation in an oil field in southwestern Iran | |
| White | 3D architecture of the Aquistore reservoir: Implications for CO2 flow and storage capacity | |
| CN107765297A (zh) | 一种储层三维地质建模方法 | |
| Zerilli et al. | Resolving Complex Overthrust Features Using Magnetotellurics-The Bolivian Foothills Case Study | |
| Samantray et al. | Interpretation and application of borehole image logs in a new generation of reservoir models for a cluster of fields in southern Oman | |
| Barton et al. | Integrated Multi-Disciplinary Approach to Predict Reservoir Performance: The Diyab Reservoir Abu Dhabi | |
| CN107765333A (zh) | 一种储层三维地质建模方法 | |
| Orellana et al. | Role of Geostatistical Techniques on the Evolution of 3D Reservoir Interpretation and Modelling Outcomes-An Example | |
| Zhang et al. | Integrated geostatistical seismic inversion to predict a thin sandstone reservoir under fan delta sedimentary environment: a case study in the Junggar Basin |
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 |