CN109444958B - 致密储层孔隙度定量预测的方法 - Google Patents
致密储层孔隙度定量预测的方法 Download PDFInfo
- Publication number
- CN109444958B CN109444958B CN201811295946.XA CN201811295946A CN109444958B CN 109444958 B CN109444958 B CN 109444958B CN 201811295946 A CN201811295946 A CN 201811295946A CN 109444958 B CN109444958 B CN 109444958B
- Authority
- CN
- China
- Prior art keywords
- impedance
- porosity
- frequency component
- seismic
- frequency
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 23
- 238000005056 compaction Methods 0.000 claims abstract description 21
- 238000013213 extrapolation Methods 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 6
- 238000011160 research Methods 0.000 claims description 6
- 238000005553 drilling Methods 0.000 claims description 5
- 230000009466 transformation Effects 0.000 claims description 5
- 230000008021 deposition Effects 0.000 claims description 4
- 238000009499 grossing Methods 0.000 claims description 3
- 238000001228 spectrum Methods 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 2
- 238000011161 development Methods 0.000 abstract description 3
- 239000007787 solid Substances 0.000 abstract 1
- 238000006243 chemical reaction Methods 0.000 description 4
- 230000000704 physical effect Effects 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000005755 formation reaction Methods 0.000 description 3
- 239000011435 rock Substances 0.000 description 3
- 238000009933 burial Methods 0.000 description 2
- 238000012217 deletion Methods 0.000 description 2
- 230000037430 deletion Effects 0.000 description 2
- 235000021185 dessert Nutrition 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 239000000047 product Substances 0.000 description 2
- 239000013589 supplement Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 235000011850 desserts Nutrition 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 239000013049 sediment Substances 0.000 description 1
Images
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/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
- G01V2210/6244—Porosity
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明属于油气勘探开发技术领域，公开了一种致密储层孔隙度定量预测的方法。本发明的一种致密储层孔隙度定量预测的方法，包括步骤：1)利用研究盆地压实趋势构建超低频分量，得到超低频分量数据；2)利用压实趋势形成甚低频分量数据体；3)利用剩余阻抗曲线和低频变化趋势体得到较低频分量阻抗模型；4)将得到的超低频分量、甚低频分量和较低频分量合并，得到具有预测性的低频模型；5)进行地震约束反演，和地震中频阻抗合并得到全频的纵波阻抗，将反演阻抗转换为孔隙度体，得到最终的孔隙度定量预测成果。本发明可以很好的建立一个地区的低频模型，开展定量储层预测和孔隙度预测，为致密油气井位目标优选打下坚实基础。
Description
技术领域
本发明属于油气勘探开发技术领域，具体涉及一种致密储层孔隙度定量预测的方法。
背景技术
在致密油气勘探或复杂油气勘探时，除了和常规储层研究一样需要关注岩性的空间分布外，更需要关注孔隙度在空间内的变化规律，在致密储层中寻找相对高孔的甜点，这时储层物性的预测精度就显得尤为重要。
要想预测好三维空间里孔隙度的变化，一个好的全频波阻抗体是非常重要的，因为波阻抗是可以通过地震反演得到的数据体，同时波阻抗是岩石速度和密度的乘积，而速度和密度是直接和岩石的孔隙度相关的，大量实测统计也表明，岩石的孔隙度和波阻抗具有很好的相关性。
要想得到全频波阻抗体，需要在地震基础上开展反演，把地震波组特征去掉子波影响，得到反射系数序列，并进一步转换为地震波阻抗。
但在地震反演过程中，一个比较棘手的问题是，地震数据缺失低频信息，通常8Hz以下的信息是采集不到的，我们不得不通过其他方法来补充低频。
在油气勘探中，测井资料显示通常的频宽很宽，因此业界比较通用的补充低频的做法是利用研究区测井曲线进行插值，以此来补充低频，但很多时候，井在空间的分布不是均匀的，同时各井的测井阻抗和孔隙度差异较大，造成低频插值时绕井点画圈现象(如图1所示为常规用井插值低频模型存在问题示意图)，这显然不是地下真实的孔隙分布规律，直接影响到物性的预测精度。
在复杂致密储层地区(如断裂比较复杂、埋深差异较大、物性致密、井的分布少且不均匀的地区)如何补充地震缺失的低频，来建立合理的低频分量，进行合理的孔隙度等物性预测是业界一大难题。
发明内容
为了解决现有技术存在的上述问题，本发明目的在于提供一种致密储层孔隙度定量预测的方法。本发明探索出一套分频段补充地震频带缺失的低频的方法，进一步细分为超低频(约0-1Hz)、甚低频(约1-3Hz)、较低频(约3-8Hz)，之后分别构建的思路，较好解决了孔隙度定量预测的问题，为致密储层甜点优选奠定基础。
本发明所采用的技术方案为：
一种致密储层孔隙度定量预测的方法，包括如下具体步骤：
1)利用研究盆地所在区域压实趋势构建超低频分量，得到超低频分量数据；
2)利用去掉压实趋势后的多个井的测井阻抗曲线统计，得到一个1-3Hz频率的多井一致性的规律，然后利用断层和层位建立构造格架，之后用去掉压实趋势后的测井阻抗进行格架控制的空间插值、内插和外推，形成甚低频分量数据体；
3)利用剩余阻抗曲线和低频变化趋势体得到较低频分量阻抗模型；
4)将步骤1)中得到的超低频分量、步骤2)中得到的甚低频分量和步骤3)中得到的较低频分量合并，得到具有预测性的低频模型；
5)进行地震约束反演，和地震中频阻抗合并得到全频的纵波阻抗，将反演阻抗转换为孔隙度体，得到最终的孔隙度定量预测成果。
进一步优选的是，在步骤1)中，通过将在盆地实际钻井中得到的测井层速度和密度曲线进行相乘，得到测井阻抗曲线，并通过多个井的测井阻抗曲线的拟合，得到对应的盆地的阻抗压实趋势，然后以时间或深度为自变量，利用阻抗随时间或深度的变化规律进行函数变换，建立研究盆地所在区域的超低频三维数据体，得到超低频分量数据。
更进一步优选的是，在步骤3)中，对原始地震资料处理得到的地震速度谱进行去野值处理，并进行空间平滑，得到一个地震速度低频变化趋势体。
更进一步优选的是，将各井测井阻抗曲线减去压实趋势及纵向岩性序列变化趋势，得到剩余阻抗曲线，然后在构造和断层格架基础上将剩余阻抗曲线进行格架控制的空间内插外推，并用地震速度低频变化趋势体作为协克里金约束，得到反映各沉积时期平面沉积变化的阻抗特征，过滤到8Hz以下，得到较低频分量阻抗模型。
更进一步优选的是，在步骤5)中，利用地震资料、步骤4)中建立的低频模型以及本研究盆地所在区域提取的子波进行地震约束反演，和地震中频阻抗合并得到全频的纵波阻抗，利用井上统计的阻抗和孔隙度变化关系，将反演合并得出的纵波阻抗转换为孔隙度体，得到最终的孔隙度定量预测成果。
更进一步优选的是，利用井上统计的阻抗和孔隙度变化关系拟合得到阻抗和孔隙度的转换函数，应用在反演合并得出的纵波阻抗上，得到所述孔隙度体。
本发明的有益效果为：
常规孔隙度三维模型是直接利用测井解释成果插值，或结合某种地震属性变换。但受井的分布、测井解释大小、断层复杂程度、埋深差异等影响很大，仅靠插值算法的改进很难得到好的经得起钻井检验的孔隙度数据体，本发明利用不同频带的规律分别建立各频带的数据，再用合并综合的方法，可以很好的建立一个地区的低频模型，开展定量储层预测和孔隙度预测，为致密储层勘探开发打下坚实基础。
附图说明
图1是常规用井插值低频模型存在问题示意图；
图2是盆地的阻抗压实趋势图；
图3是合并后的低频模型示意图；
图4是反演转换后的孔隙度体图。
具体实施方式
下面结合具体实施方式对本发明作进一步的说明。本发明探索出一套细分频段补充地震频带缺失的低频的方法，将低频段进一步细分为超低频(约0-1Hz)、甚低频(约1-3Hz)、较低频(约3-8Hz)，之后对各频段分别构建阻抗模型的思路，来提高低频分量的可靠性，较好解决了孔隙度定量预测的问题，为致密储层甜点优选奠定基础。
如图2-4所示，本实施例提供一种致密储层孔隙度定量预测的方法，包括如下具体步骤：
1)通过将在研究盆地所在区域实际钻井中得到的测井层速度和密度曲线进行相乘，得到测井阻抗曲线，并通过多个井的测井阻抗曲线的拟合，得到对应的盆地的阻抗压实趋势，如图2所示；
由于阻抗压实趋势是通过多井统计得到，就可以避免由于某一口井测量误差带来的影响，得到对应的盆地的阻抗压实趋势以时间或深度为自变量，利用阻抗随时间或深度的变化规律进行函数变换，建立研究盆地所在区域的超低频三维数据体，得到超低频分量数据，构建的数据体仅和时间/深度相关，和层位及断层无关。
2)在盆地压实趋势之外，另一个比较重要的影响因素是纵向岩性序列变化趋势，即不同时期沉积地层的整体特征，在同一沉积地层，其储地比基本类似，不同沉积地层的沉积特征不同，这种不同时代沉积物的所造成的差异主要表现在甚低频频段(约1-3Hz)；
利用去掉压实趋势后的多个井的测井阻抗曲线统计，得到一个1-3Hz频率的多井一致性的规律，然后利用断层和层位建立构造格架，之后用去掉压实趋势后的测井阻抗进行格架控制的空间插值、内插和外推，形成甚低频分量数据体；
3)对原始地震资料处理得到的地震速度谱进行去野值处理，并进行空间平滑，得到一个地震速度低频变化趋势体；然后将各井测井阻抗曲线减去前面提到的压实趋势及纵向岩性序列变化趋势，得到剩余阻抗曲线，然后在构造和断层格架基础上将剩余阻抗曲线进行格架控制的空间内插外推，并用地震速度低频变化趋势体作为协克里金约束，得到反映各沉积时期平面沉积变化的阻抗特征，过滤到8Hz以下，得到较低频分量阻抗模型(大约3-8Hz)。
4)将步骤1)中得到的超低频分量、步骤2)中得到的甚低频分量和步骤3)中得到的较低频分量合并，得到更具有预测性的低频模型，如图3所示；
5)利用地震资料、步骤4)中建立的低频模型以及本研究盆地所在区域提取的子波进行地震约束反演，和地震中频阻抗合并得到全频的纵波阻抗，利用井上统计的阻抗和孔隙度变化关系，拟合得到阻抗和孔隙度的转换函数，应用在反演合并得出的纵波阻抗上，从而得到所述孔隙度体，如图4所示，得到最终的孔隙度定量预测成果，可以以此指导致密储层的油气勘探开发。
本发明中常规孔隙度三维模型是直接利用测井解释成果插值，或结合某种地震属性变换。但受井的分布、测井解释大小、断层复杂程度、埋深差异等影响很大，仅靠插值算法的改进很难得到好的经得起钻井检验的孔隙度数据体，本发明利用不同频带的规律分别建立各频带的数据，再用合并综合的方法，可以很好的建立一个地区的低频模型，开展定量储层预测和孔隙度预测，为致密储层勘探开发打下坚实基础。
本发明不局限于上述可选的实施方式，任何人在本发明的启示下都可得出其他各种形式的产品。上述具体实施方式不应理解成对本发明的保护范围的限制，本发明的保护范围应当以权利要求书中界定的为准，并且说明书可以用于解释权利要求书。
Claims (6)
1.一种致密储层孔隙度定量预测的方法，其特征在于：包括如下具体步骤：
1）利用研究盆地所在区域压实趋势构建超低频分量，得到超低频分量数据；
2）利用去掉压实趋势后的多个井的测井阻抗曲线统计，得到一个1-3Hz频率的多井一致性的规律，然后利用断层和层位建立构造格架，之后用去掉压实趋势后的测井阻抗进行格架控制的空间插值、内插和外推，形成甚低频分量数据体；
3）利用剩余阻抗曲线和低频变化趋势体得到较低频分量阻抗模型；
4）将步骤1）中得到的超低频分量、步骤2）中得到的甚低频分量和步骤3）中得到的较低频分量合并，得到具有预测性的低频模型；
5）进行地震约束反演，和地震中频阻抗合并得到全频的纵波阻抗，将反演阻抗转换为孔隙度体，得到最终的孔隙度定量预测成果。
2.根据权利要求1所述的致密储层孔隙度定量预测的方法，其特征在于：在步骤1）中，通过将在盆地实际钻井中得到的测井层速度和密度曲线进行相乘，得到测井阻抗曲线，并通过多个井的测井阻抗曲线的拟合，得到对应的盆地的阻抗压实趋势，然后以时间或深度为自变量，利用阻抗随时间或深度的变化规律进行函数变换，建立研究盆地所在区域的超低频三维数据体，得到超低频分量数据。
3.根据权利要求1所述的致密储层孔隙度定量预测的方法，其特征在于：在步骤3）中，对原始地震资料处理得到的地震速度谱进行去野值处理，并进行空间平滑，得到一个地震速度低频变化趋势体。
4.根据权利要求3所述的致密储层孔隙度定量预测的方法，其特征在于：将各井测井阻抗曲线减去压实趋势及纵向岩性序列变化趋势，得到剩余阻抗曲线，然后在构造和断层格架基础上将剩余阻抗曲线进行格架控制的空间内插外推，并用地震速度低频变化趋势体作为协克里金约束，得到反映各沉积时期平面沉积变化的阻抗特征，过滤到8Hz以下，得到较低频分量阻抗模型。
5.根据权利要求4所述的致密储层孔隙度定量预测的方法，其特征在于：在步骤5）中，利用地震资料、步骤4）中建立的低频模型以及本研究盆地所在区域提取的子波进行地震约束反演，和地震中频阻抗合并得到全频的纵波阻抗，利用井上统计的阻抗和孔隙度变化关系，将反演合并得出的纵波阻抗转换为孔隙度体，得到最终的孔隙度定量预测成果。
6.根据权利要求5所述的致密储层孔隙度定量预测的方法，其特征在于：利用井上统计的阻抗和孔隙度变化关系拟合得到阻抗和孔隙度的转换函数，应用在反演合并得出的纵波阻抗上，得到所述孔隙度体。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811295946.XA CN109444958B (zh) | 2018-11-01 | 2018-11-01 | 致密储层孔隙度定量预测的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811295946.XA CN109444958B (zh) | 2018-11-01 | 2018-11-01 | 致密储层孔隙度定量预测的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109444958A CN109444958A (zh) | 2019-03-08 |
| CN109444958B true CN109444958B (zh) | 2020-05-26 |
Family
ID=65550540
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811295946.XA Active CN109444958B (zh) | 2018-11-01 | 2018-11-01 | 致密储层孔隙度定量预测的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109444958B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113093286B (zh) * | 2021-03-15 | 2022-08-02 | 中国科学院海洋研究所 | 一种冷泉发育区储层非均质性的反演方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101206264A (zh) * | 2007-11-08 | 2008-06-25 | 符力耘 | 高分辨率非线性地震波阻抗反演方法 |
| WO2014164390A2 (en) * | 2013-03-11 | 2014-10-09 | Saudi Arabian Oil Company | Low frequency passive seismic data acquisition and processing |
| CN105425289A (zh) * | 2015-10-29 | 2016-03-23 | 中国石油天然气集团公司 | 确定低频波阻抗的方法和装置 |
| CN105866838A (zh) * | 2016-05-17 | 2016-08-17 | 中国石油天然气股份有限公司 | 一种地震数据低频信息补偿方法及装置 |
| CN106324675A (zh) * | 2016-10-09 | 2017-01-11 | 中国石油大学(华东) | 一种宽频地震波阻抗低频信息预测方法及系统 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8332154B2 (en) * | 2009-06-02 | 2012-12-11 | Exxonmobil Upstream Research Company | Estimating reservoir properties from 4D seismic data |
-
2018
- 2018-11-01 CN CN201811295946.XA patent/CN109444958B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101206264A (zh) * | 2007-11-08 | 2008-06-25 | 符力耘 | 高分辨率非线性地震波阻抗反演方法 |
| WO2014164390A2 (en) * | 2013-03-11 | 2014-10-09 | Saudi Arabian Oil Company | Low frequency passive seismic data acquisition and processing |
| CN105425289A (zh) * | 2015-10-29 | 2016-03-23 | 中国石油天然气集团公司 | 确定低频波阻抗的方法和装置 |
| CN105866838A (zh) * | 2016-05-17 | 2016-08-17 | 中国石油天然气股份有限公司 | 一种地震数据低频信息补偿方法及装置 |
| CN106324675A (zh) * | 2016-10-09 | 2017-01-11 | 中国石油大学(华东) | 一种宽频地震波阻抗低频信息预测方法及系统 |
Non-Patent Citations (1)
| Title |
|---|
| 高分辨率非线性地震波阻抗反演方法和应用;崔炯成 等;《石油物探》;20070131;第46卷(第1期);第1-12页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109444958A (zh) | 2019-03-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108931814B (zh) | 一种基于多属性融合的基岩裂缝预测的方法 | |
| CN103389513B (zh) | 应用声波测井资料约束反演提高地震资料分辨率的方法 | |
| CN106842307B (zh) | 一种基于正演约束下波形分类再检索的储层精细预测方法 | |
| CN109061764B (zh) | 一种分频融合波阻抗反演方法 | |
| CN111461386B (zh) | 基于bp神经网络的页岩气甜点预测方法 | |
| NO337423B1 (no) | Fremgangsmåte for oppdatering av en geologisk modell med seismiske data | |
| CN111522063B (zh) | 基于分频联合反演的叠前高分辨率流体因子反演方法 | |
| Khoshdel et al. | Multi attribute transform and neural network in porosity estimation of an offshore oil field—A case study | |
| CN109212627B (zh) | 一种基于细分层系的薄储层描述及增储规模确定方法 | |
| CN105116449A (zh) | 一种弱反射储层的识别方法 | |
| CN108508489A (zh) | 一种基于波形微变化匹配的地震反演方法 | |
| CN109444958B (zh) | 致密储层孔隙度定量预测的方法 | |
| CN114862609A (zh) | 一种特高含水油田的油藏开采方法 | |
| CN105388519A (zh) | 一种提高地震资料分辨率的方法 | |
| CN110568490B (zh) | 一种高速层顶薄储层的识别方法 | |
| CN111077578B (zh) | 岩层分布预测方法和装置 | |
| CN115012903B (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 | |
| CN115880455A (zh) | 基于深度学习的三维智能插值方法 | |
| CN110568149B (zh) | 沉积盆地烃源岩生排烃史精细快速定量模拟方法 | |
| Pan et al. | An approach to reserve estimation enhanced with 3-D seismic data | |
| CN111335871A (zh) | 一种基于分层产能评价的查层补孔技术方法 | |
| CN110795513A (zh) | 河流相源储异位型致密油气甜点区分布的预测方法 | |
| US11542816B2 (en) | System and method for hydrocarbon exploration based on imaging tunnel valleys | |
| CN113361111B (zh) | 一种低渗透油藏储层特征模型的划分方法 |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |