CN107179562A - 相控岩石物理模型指导下的储层预测方法 - Google Patents
相控岩石物理模型指导下的储层预测方法 Download PDFInfo
- Publication number
- CN107179562A CN107179562A CN201710289032.1A CN201710289032A CN107179562A CN 107179562 A CN107179562 A CN 107179562A CN 201710289032 A CN201710289032 A CN 201710289032A CN 107179562 A CN107179562 A CN 107179562A
- Authority
- CN
- China
- Prior art keywords
- reservoir
- rock
- petrophysical model
- mineral
- sedimentary facies
- 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 30
- 239000011435 rock Substances 0.000 claims abstract description 87
- 229910052500 inorganic mineral Inorganic materials 0.000 claims abstract description 58
- 239000011707 mineral Substances 0.000 claims abstract description 58
- 208000035126 Facies Diseases 0.000 claims abstract description 51
- 239000000470 constituent Substances 0.000 claims abstract description 9
- 230000008859 change Effects 0.000 claims description 20
- 238000004088 simulation Methods 0.000 claims description 18
- 239000012530 fluid Substances 0.000 claims description 16
- 239000000203 mixture Substances 0.000 claims description 15
- 238000005070 sampling Methods 0.000 claims description 9
- 238000011160 research Methods 0.000 claims description 6
- 235000021321 essential mineral Nutrition 0.000 claims description 5
- 238000010008 shearing Methods 0.000 claims description 5
- 230000008021 deposition Effects 0.000 claims description 4
- 230000005251 gamma ray Effects 0.000 claims description 4
- 238000005259 measurement Methods 0.000 claims description 4
- 241001074085 Scophthalmus aquosus Species 0.000 claims description 3
- 239000004927 clay Substances 0.000 description 6
- 239000011159 matrix material Substances 0.000 description 6
- 230000007246 mechanism Effects 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 5
- 239000002734 clay mineral Substances 0.000 description 5
- 239000004575 stone Substances 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 241000208340 Araliaceae Species 0.000 description 3
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 3
- 235000003140 Panax quinquefolius Nutrition 0.000 description 3
- 238000005553 drilling Methods 0.000 description 3
- 235000008434 ginseng Nutrition 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 239000006185 dispersion Substances 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 238000002347 injection Methods 0.000 description 2
- 239000007924 injection Substances 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- 239000011148 porous material Substances 0.000 description 2
- 239000000523 sample Substances 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 229910021532 Calcite Inorganic materials 0.000 description 1
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 229910000514 dolomite Inorganic materials 0.000 description 1
- 239000010459 dolomite Substances 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 239000010433 feldspar Substances 0.000 description 1
- 239000011799 hole material Substances 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000012067 mathematical method Methods 0.000 description 1
- 238000002156 mixing Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000001615 p wave Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 239000010453 quartz Substances 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 241000894007 species Species 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提出了一种相控岩石物理模型指导下的储层预测方法，包括：划分沉积相；得到各个沉积相内储层主要矿物成分的不同；初步确定每个沉积相内的岩石物理建模参数；建立每个沉积相内储层初步的岩石物理模型；得到模拟数据和实测数据的差异，若该差异超过预设误差允许值，则调整建模参数，重新建立岩石物理模型，直至得到每个沉积相内最终的岩石物理模型；确定对特征矿物变化最为敏感的两种弹性参数；计算得到该两种弹性参数的关系曲线；得到整个三维地震区的两种弹性参数分布；将整个三维地震区的两种弹性参数的分布划分为储层和非储层，得到整个三维地震区的储层分布。本发明能够在复杂地质条件下，提高储层预测的准确度。
Description
技术领域
本发明涉及石油勘探技术领域，特别涉及一种相控岩石物理模型指导下的储层预测方法。
背景技术
随着现代社会对油气资源需求的与日俱增，油气的勘探开发的程度不断深入，但可供勘探的大型圈闭、“亮点”型油气藏越来越少，油气勘探所面临的地质情况越来越复杂，逐渐转向岩性圈闭、复合型圈闭、非“亮点”型油气藏等，这就造成常规勘探技术难以识别油气分布，例如常规的波阻抗反演和叠前反演，在面对复杂构造和沉积情况时，往往无法识别储层，或者具有多解性，因此，需要一定的方法来提高储层预测技术的准确性，对于降低勘探风险，提高勘探成功率具有重要作用。
岩石物理分析是一种通过分析岩石的物理特征来识别储层和非储层的方法。可以通过对已钻井分析，找出储层和非储层的岩石物理参数的差异，比如纵波速度、横波速度、泊松比等的差异，如果其中一种或几种岩石物理参数能够以某一个数值为界限，很好的将储层和非储层区分开，就可以通过地震反演来预测这个参数的分布，进而预测出储层分布。但是由于地下地质条件的复杂性，往往不是一个弹性参数的数值就可以把储层界定出来的，储层和非储层的界线可能是两个弹性参数分别为纵横坐标而形成的一条曲线，这种情况下就无法通过钻井分析准确的把这条界线画出来。因此，需要建立岩石物理模型，以数字方法模拟真实的岩石特征，从正演的角度来找出储层和非储层的分界线。建立岩石物理模型首先需要确定岩石的各项组成参数及环境参数，包括矿物成分和含量，孔隙度、孔隙形状、孔隙流体类型、流体的饱和度、温度、压力等，这些参数可以从已钻井数据中获得。然后根据岩性选取合适的理论模型，理论模型是已知矿物、孔隙流体的物理性质，根据孔隙特征、环境特征来推导出整个岩石的物理性质的数学方法，包括Gassmann模型、Biot模型、Xu-White模型等，这些模型都有一定的适用条件，利用这些理论模型就可以建立出岩石的物理模型。有了岩石物理模型，就可以在计算机中，通过改变模型的矿物含量或孔隙度，分析岩石物理参数的变化规律，例如泊松比随孔隙度的变化规律。如果储层和非储层定义为以岩石中某种矿物的特定含量为界，就可以先设定岩石物理模型中此矿物为该特定含量，然后改变模型孔隙度，得到各岩石物理参数随孔隙度变化的曲线，此曲线就是用岩石物理参数值识别储层和非储层的分界线。
然而地质条件通常是复杂多变的，在一个研究区块，储层在纵向和横向上可能分为多种类型，储层岩性或矿物成分都有所不同，这种情况下应用一个同样的岩石物理模型就无法准确预测整个地区的储层分布。
发明内容
本发明的目的旨在至少解决所述技术缺陷之一。
为此，本发明的目的在于提出一种相控岩石物理模型指导下的储层预测方法，能够在复杂地质条件下，提高储层预测的准确度。
为了实现上述目的，本发明的实施例提供一种相控岩石物理模型指导下的储层预测方法，包括如下步骤：
步骤S1，利用已有的测井数据、地震数据和地质资料划分沉积相；
步骤S2，利用钻井岩心、薄片鉴定结果，得到上述划分的各个沉积相内储层主要矿物成分的不同；
步骤S3，利用所述各个沉积相内储层主要矿物成分的不同，初步确定每个沉积相内的岩石物理建模参数；
步骤S4，利用上述确定的每个所述沉积相的岩石物理建模参数，建立每个所述沉积相内储层初步的岩石物理模型；
步骤S5，利用所述岩石物理模型，模拟纵波速度、横波速度和密度，将所述模拟纵波速度、横波速度和密度与已知钻井中测量的实际纵波速度、横波速度、密度进行对比，得到模拟数据和实测数据的差异，若该差异超过预设误差允许值，则调整建模参数，重新建立岩石物理模型，并再次进行模拟和实测数据的对比，如此反复直至所述模拟和实测数据的差异位于所述预设误差允许值内，则得到每个所述沉积相内最终的岩石物理模型；
步骤S6，利用所述每个沉积相内最终的岩石物理模型，模拟当岩石中能界定储层和非储层的特征矿物在含量变化时，岩石弹性参数的变化规律，根据所述岩石弹性参数的变化规律，确定对该特征矿物变化最为敏感的两种弹性参数；
步骤S7，将所述两种弹性参数，分别为纵轴和横轴，建立坐标系，在岩石物理模型中，把特征矿物的含量设置为储层和非储层界限值，变化岩石的孔隙度，即以孔隙度为参变量，计算得到该两种弹性参数的关系曲线；
步骤S8，利用测井数据和叠前地震数据，进行叠前弹性参数反演，得到整个三维地震区的两种弹性参数分布；
步骤S9，在两种最敏感弹性参数组成的坐标系中，以两种弹性参数的关系曲线为界，将所述整个三维地震区的两种弹性参数的分布划分为储层和非储层，得到整个三维地震区的储层分布。
进一步，在所述步骤S1中，
所述测井数据包括自然伽马曲线、电阻率曲线、声波时差曲线、密度曲线；
所述地震数据包括三维叠后纯波数据、地震构造解释成果；
所述地质资料包括岩心描述、区域构造、区域沉积、以及现有研究成果。
进一步，所述步骤S2还包括如下步骤：统计多个岩心、薄片的分析鉴定结果，得到每个所述沉积相内岩石的矿物组成。
进一步，在所述步骤S3中，利用所述各个沉积相内储层主要矿物成分的不同，初步确定每个沉积相的岩石物理建模参数，包括如下步骤：初步确定模型的主要矿物组分的岩石物理参数、孔隙中所含流体类型及流体混合方式、根据岩性选择合适的理论模型。
进一步，在所述步骤S5中，所述模拟数据和实测数据随埋深变化的趋势要一致，整个地层段内平均误差不超过3％，建模参数应在矿物的实际物理性质范围内调整。
进一步，在所述步骤S6中，所述岩石弹性参数包括：纵波速度、横波速度、密度、泊松比、拉梅系数、剪切模量、体积模量。
进一步，以变化的百分比为单位，建立岩石各种弹性参数随特征矿物含量变化的关系曲线。
进一步，在所述步骤S6中，根据所述岩石弹性参数的变化规律，确定对该特征矿物变化最为敏感的两种弹性参数，包括：比较岩石各种弹性参数随特征矿物含量变化的改变幅度，幅度最大的最为敏感，以此确定对该特征矿物变化最为敏感的两种弹性参数。
进一步，在所述步骤S8中，所述测井数据包括：实测纵波时差曲线、横波时差曲线、密度曲线；所述叠前地震数据包括三维叠前共反射点CRP地震道集。
进一步，在所述步骤S9中，以两种弹性参数的值为坐标，将整个三维地震区内每个样点投到该坐标系中，提取位于储层部分的所有样点，上述样点在三维地震区内的位置即为储层位置。
根据本发明实施例的相控岩石物理模型指导下的储层预测方法，使用沉积相来控制每种储层类型的范围，在不同的沉积相内使用不同的岩石物理模型来分别进行储层预测，对多种类型的储层采用针对性的岩石物理建模方法，分别进行预测，从而降低预测的多解性，提高勘探成功率，并且能够在复杂地质条件下，提高储层预测的准确度。
本发明附加的方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解，其中：
图1为根据本发明实施例的相控岩石物理模型指导下的储层预测方法的流程图；
图2为根据本发明实施例的沉积相划分示意图；
图3为根据本发明实施例的某井碎屑岩中粘土类矿物在不同深度的含量图；
图4为根据本发明实施例的某井模拟与实测的纵波速度、横波速度、密度对比图；
图5为根据本发明实施例的模拟的某井碎屑岩储层中各弹性参数变化量与粘土类矿物含量的关系示意图；
图6为根据本发明实施例的模拟的某井中随孔隙度变化，纵波速度和密度的关系曲线。
具体实施方式
下面详细描述本发明的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，旨在用于解释本发明，而不能理解为对本发明的限制。
如图1所示，本发明实施例的相控岩石物理模型指导下的储层预测方法，包括如下步骤：
步骤S1，利用已有的测井数据、地震数据和地质资料划分沉积相。其中，所述测井数据包括自然伽马曲线、电阻率曲线、声波时差曲线、密度曲线；所述地震数据包括三维叠后纯波数据、地震构造解释成果；所述地质资料包括岩心描述、区域构造、区域沉积、以及现有研究成果。
如图2所示，利用已有的地质研究成果，对勘探区沉积相分布进行总体控制。然后，利用区内每口井的测井数据，主要包括自然伽马测井曲线、岩性解释成果、岩心照片、岩心描述记录，对每口井进行纵向的沉积相划分，利用能够反映沉积相分布的地震属性，通常为均方根振幅、分频能量属性，识别每套地层上各种沉积相的平面分布，最后结合井点处的沉积相和地震属性，综合分析调整，在平面上划分出每套地层的沉积相。
步骤S2，利用钻井岩心、薄片鉴定结果，得到上述划分的各个沉积相内储层主要矿物成分的不同。
岩石是由矿物组成的，对于碎屑岩储层，主要矿物成分为石英、长石、粘土等，对于碳酸盐岩储层，主要矿物成分为方解石、白云石等，岩心和薄片鉴定能够比较精确的分析出岩石中每种矿物的含量。如图3所示，在本步骤中，统计多个岩心、薄片的分析鉴定结果，得到每个所述沉积相内岩石的矿物组成。但缺点是样本较少，不能得出整个区域矿物的确切含量，只能定性的对比各种矿物的多少及含量范围，为岩石物理建模参数选取提供参考。
步骤S3，利用所述各个沉积相内储层主要矿物成分的不同，初步确定每个沉积相内岩石的建模参数。
由于不同沉积相的矿物成分不同，岩石物理建模参数也不同，甚至是所采用的建模方法也会不同。自然界中同一种矿物的岩石物理参数都有特定的数值或范围，根据各沉积相矿物的种类差异和自然界中每种矿物的岩石物理参数值，可以初步确定建模参数和建模方法。其中，岩石物理参数包括纵横波速度、密度、弹性模量、剪切模量等。
具体地，利用所述各个沉积相内储层主要矿物成分的不同，初步确定每个沉积相的岩石物理建模参数，包括如下步骤：初步确定模型的主要矿物组分的岩石物理参数、孔隙中所含流体类型及流体混合方式、根据岩性选择合适的理论模型。
步骤S4，利用上述确定的每个所述沉积相内的岩石物理建模参数，建立每个所述沉积相内储层初步的岩石物理模型。
岩石由骨架和孔隙流体组成，骨架由基质和孔隙组成，基质由矿物组成，因此岩石的弹性参数由四个分量决定：基质的弹性参数、干岩石骨架的弹性参数、孔隙流体的弹性参数、环境因素(压力、温度等)，岩石物理建模就是要建立这些分量间的理论关系。常用的理论模型有以下几种：
1、Gassmann模型：在低频条件下，Gassmann推导出了饱和流体状态下岩石弹性参数的理论方程。Gassmann方程用来描述从干岩石状态到饱含流体状态下的岩石弹性参数的变化，该方程只有在足够低频条件下，才是有效的，此时孔隙所受的压力在整个孔隙空间达到平衡。
2、Biot模型：Biot采用连续介质力学的方法导出了流体饱和多孔隙介质中的声波方程，建立了多孔介质中的声速、衰减与频率和多孔介质参数之间的关系。该模型反映了流体和岩石骨架中粘性和惯性相互作用机制，即包含了岩石骨架和孔隙流体对混合岩石介质弹性参数的单独作用，也包含了它们之间的耦合作用。该模型适合于任意频率条件下多孔岩石介质弹性模量的计算。
3、BISQ模型：当地震波在多孔介质中传播时，流体和固体相互作用的机制有两种：Biot流和喷射流，两种机制同时存在，Biot流描述的是宏观现象，喷射流机制反映的是局部特征，两种机制通过流体的质量守恒而统一，对地震波的衰减和频散均会产生重要影响。Dvorkin和Nur基于孔隙各向同性一维问题将这两种流体-固体相互作用的力学机制有机的结合起来，提出了同一的Biot-Squirt(BISQ)模型。BISQ模型反映了两种不同流动形式和流体特性对波速、衰减和频散的影响规律，比Biot理论更能真是的体现波在孔隙各向异性岩石介质中的传播规律。
4、Xu-White模型：基于Kuster-Toksoz模型和Gassmann理论，Xu-White提出了砂泥岩混合介质的速度模型。该模型综合考虑岩石孔隙度和粘土含量来预测声波速度，把粘土成分、压力、胶结等因素对声波的影响归于泥页岩和砂岩的孔隙几何形状和面孔率的差异。
由于多孔介质岩石物理性质的复杂性，不同的理论模型有其不同的使用条件及关键参数。在实际应用中，根据研究区的实际情况选用合适的理论模型进行岩石物理建模。
步骤S5，利用所述岩石物理模型，模拟纵波速度、横波速度和密度，将所述模拟的纵波速度、横波速度和密度与已知钻井中测量的实际纵波速度、横波速度、密度进行对比，得到模拟数据和实测数据的差异，若该差异超过预设误差允许值，则调整建模参数，重新建立岩石物理模型，并再次进行模拟和实测数据的对比，如此反复直至所述模拟和实测数据的差异位于所述预设误差允许值内，则得到每个所述沉积相内最终的岩石物理模型。
在本步骤中，利用所建立的岩石物理模型，经过理论计算，可以模拟出岩石的纵波速度、横波速度和密度。
在本发明的一个实施例中，所述模拟数据和实测数据随埋深变化的趋势要一致，整个地层段内平均误差不超过3％，建模参数应在矿物的实际物理性质范围内调整。
如图4所示，将模拟出的岩石纵波速度、横波速度和密度与测井中实测的纵波速度、横波速度、密度进行比较，分析模拟与实测的误差，如果误差超过预设误差允许值，说明建立岩石物理模型使用的参数不对，需要调整，参数应该在合理的范围内调整，比如地层压力、矿物的弹性参数等不能与测量值或经验值相差太远。调整参数后重新模拟纵波速度、横波速度、密度，并再与实测值对比，直到误差小于超过预设误差允许值，就可以认为所建立的岩石物理模型是合理的。
在本发明的一个实施例中，预设误差允许值位于2％～5％。优选的，预设误差允许值为3％。
步骤S6，利用所述每个沉积相内最终的岩石物理模型，模拟当岩石中能界定储层和非储层的特征矿物在含量变化时，岩石弹性参数的变化规律，根据所述岩石弹性参数的变化规律，确定对该特征矿物变化最为敏感的两种弹性参数。
在本发明的一个实施例中，所述岩石弹性参数包括：纵波速度、横波速度、密度、泊松比、拉梅系数、剪切模量、体积模量。
以变化的幅度(百分比)为单位，建立岩石各种弹性参数随特征矿物含量变化的关系曲线。
具体地，根据所述岩石弹性参数的变化规律，确定对该特征矿物变化最为敏感的两种弹性参数，包括：比较岩石各种弹性参数随特征矿物含量变化的改变幅度，幅度最大的最为敏感，以此确定对该特定矿物变化最为敏感的两种弹性参数。
岩石中某种矿物的含量可以界定岩石能否成为储层，例如在常规的碎屑岩地层中，砂岩是可以作为储层的，当粘土矿物含量小于40％时，岩石就可以定性为砂岩，此时特征矿物就是指粘土矿物。当在岩石物理模型中，改变特征矿物的含量时，整个岩石的弹性参数必然发生变化，通过计算每个弹性参数随矿物含量变化而变化的规律，可以对比各弹性参数对于特征矿物含量变化的敏感程度，因此可以找出两个最为敏感的弹性参数，如图5所示。图5是模拟的某井碎屑岩储层中各弹性参数变化量与粘土类矿物含量的关系示意图，此图中显示拉梅系数和纵波阻抗两个参数对粘土类矿物含量变化最敏感。
需要说明的是，对于每个沉积相，由于建模参数不同，这两个弹性参数可能是不同的。
步骤S7，将所述两种弹性参数，分别为纵轴和横轴，建立坐标系，在岩石物理模型中，把特征矿物的含量设定为储层和非储层界限值，变化岩石的孔隙度，即以孔隙度为参变量，计算得到该两种弹性参数的关系曲线。
需要说明的是，储层和非储层的界限值可以根据在先研究数据进行确定。
图6是模拟的某井中随孔隙度变化，纵波速度和密度的关系曲线。其中有粘土类矿物含量100％、40％、0％时的三条曲线，图中散点为实测的井中密度和纵波速度的交汇点，位于40％和100％线间的点为泥岩，位于40％和0％线间的点为砂岩。
如图6所示，指定特征矿物的含量为储层的界限值，例如对于碎屑岩，指定粘土类矿物含量为40％，这时改变岩石物理模型的孔隙度，可以模拟出各弹性参数随孔隙度变化而变化的规律。将两个对特定矿物含量最敏感的弹性参数分别作为纵坐标和横坐标，建立坐标系。以孔隙度为参变量，孔隙度从小到大变化时，在坐标系中就可以绘出一条这两个弹性参数的关系曲线。在这条曲线的两侧，即这条曲线所隔开的两个区域，分别为该坐标系中储层和非储层所在的区域。由于每个沉积相建模参数不同，需要对每个沉积相分别开展此步骤。
步骤S8，利用测井数据和叠前地震数据，进行叠前弹性参数反演，得到整个三维地震区的两种弹性参数分布。
在本发明的一个实施例中，所述测井数据包括：实测纵波时差曲线、横波时差曲线、密度曲线；所述叠前地震数据包括三维叠前共反射点CRP地震道集。
整个研究区的弹性参数分布只能通过叠前弹性参数反演获得。叠前弹性参数反演的基础理论方程为：
其中，θ为地震波入射角，EI(θ)为入射角θ时的弹性阻抗，α为纵波速度，β为横波速度、ρ为密度。
将叠前CRP道集按入射角大小分为3个部分，将这3部分内的道集分别叠加，得到3个入射角不同的叠后地震数据体。同时，因为已钻井上的纵波速度、横波速度和密度都已知，由式(1)可计算出已钻井上3个入射角的弹性阻抗曲线，由井上的弹性阻抗曲线制作合成地震记录，分别与对应的叠后地震数据进行标定，然后进行叠后波阻抗反演，得到全三维区的3个弹性阻抗数据体，代入式(1)中，建立三元联立方程组，解出全三维区的纵波速度、纵波速度、纵波速度，最后由纵波速度、纵波速度、纵波速度计算出体积模量、泊松比、拉梅系数、剪切模量等其它弹性参数。
步骤S9，在两种最敏感弹性参数组成的坐标系中，以两种弹性参数的关系曲线为界，将所述整个三维地震区的两种弹性参数的分布划分为储层和非储层，得到整个三维地震区的储层分布。
具体地，将反演得到的三维区的两个敏感参数数据体在上述步骤建立的坐标系中进行交汇，整个三维区的样点会被这两个敏感参数的关系曲线分隔为两部分，即储层和非储层。提取位于储层部分的所有样点，上述样点在三维地震区内的位置即为储层所在。
根据本发明实施例的相控岩石物理模型指导下的储层预测方法，使用沉积相来控制每种储层类型的范围，在不同的沉积相内使用不同的岩石物理模型来分别进行储层预测，对多种类型的储层采用针对性的岩石物理建模方法，分别进行预测，从而降低预测的多解性，提高勘探成功率，并且能够在复杂地质条件下，提高储层预测的准确度。
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
尽管上面已经示出和描述了本发明的实施例，可以理解的是，上述实施例是示例性的，不能理解为对本发明的限制，本领域的普通技术人员在不脱离本发明的原理和宗旨的情况下在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。本发明的范围由所附权利要求极其等同限定。
Claims (10)
1.一种相控岩石物理模型指导下的储层预测方法，其特征在于，包括如下步骤：
步骤S1，利用已有的测井数据、地震数据和地质资料划分沉积相；
步骤S2，利用钻井岩心、薄片鉴定结果，得到上述划分的各个沉积相内储层主要矿物成分的不同；
步骤S3，利用所述各个沉积相内储层主要矿物成分的不同，初步确定每个沉积相内的岩石物理建模参数；
步骤S4，利用上述确定的每个所述沉积相的岩石物理建模参数，建立每个所述沉积相内储层初步的岩石物理模型；
步骤S5，利用所述岩石物理模型，模拟纵波速度、横波速度和密度，将所述模拟纵波速度、横波速度和密度与已知钻井中测量的实际纵波速度、横波速度、密度进行对比，得到模拟数据和实测数据的差异，若该差异超过预设误差允许值，则调整建模参数，重新建立岩石物理模型，并再次进行模拟和实测数据的对比，如此反复直至所述模拟和实测数据的差异位于所述预设误差允许值内，则得到每个所述沉积相内最终的岩石物理模型；
步骤S6，利用所述每个沉积相内最终的岩石物理模型，模拟当岩石中能界定储层和非储层的特征矿物在含量变化时，岩石弹性参数的变化规律，根据所述岩石弹性参数的变化规律，确定对该特征矿物变化最为敏感的两种弹性参数；
步骤S7，将所述两种弹性参数，分别为纵轴和横轴，建立坐标系，在岩石物理模型中，把特征矿物含量设置为储层和非储层界限值，变化岩石的孔隙度，即以孔隙度为参变量，计算得到该两种弹性参数的关系曲线；
步骤S8，利用测井数据和叠前地震数据，进行叠前弹性参数反演，得到整个三维地震区的两种弹性参数分布；
步骤S9，在两种最敏感弹性参数组成的坐标系中，以两种弹性参数的关系曲线为界，将所述整个三维地震区的两种弹性参数的分布划分为储层和非储层，得到整个三维地震区的储层分布。
2.如权利要求1所述的相控岩石物理模型指导下的储层预测方法，其特征在于，在所述步骤S1中，
所述测井数据包括自然伽马曲线、电阻率曲线、声波时差曲线、密度曲线；
所述地震数据包括三维叠后纯波数据、地震构造解释成果；
所述地质资料包括岩心描述、区域构造、区域沉积、以及现有研究成果。
3.如权利要求1所述的相控岩石物理模型指导下的储层预测方法，其特征在于，所述步骤S2还包括如下步骤：统计多个岩心、薄片的分析鉴定结果，得到每个所述沉积相内岩石的矿物组成。
4.如权利要求1所述的相控岩石物理模型指导下的储层预测方法，其特征在于，在所述步骤S3中，利用所述各个沉积相内储层主要矿物成分的不同，初步确定每个沉积相的岩石物理建模参数，包括如下步骤：初步确定模型的主要矿物组分的岩石物理参数、孔隙中所含流体类型及流体混合方式、根据岩性选择合适的理论模型。
5.如权利要求1所述的相控岩石物理模型指导下的储层预测方法，其特征在于，在所述步骤S5中，所述模拟数据和实测数据随埋深变化的趋势要一致，整个地层段内平均误差不超过3％，建模参数应在矿物的实际物理性质范围内调整。
6.如权利要求1所述的相控岩石物理模型指导下的储层预测方法，其特征在于，在所述步骤S6中，所述岩石弹性参数包括：纵波速度、横波速度、密度、泊松比、拉梅系数、剪切模量、体积模量。
7.如权利要求6所述的相控岩石物理模型指导下的储层预测方法，其特征在于，以变化的百分比为单位，建立岩石各种弹性参数随特征矿物含量变化的关系曲线。
8.如权利要求1所述的相控岩石物理模型指导下的储层预测方法，其特征在于，在所述步骤S6中，根据所述岩石弹性参数的变化规律，确定对特征矿物变化最为敏感的两种弹性参数，包括：比较岩石各种弹性参数随特征矿物含量变化的改变幅度，幅度最大的最为敏感，以此确定对该特征矿物变化最为敏感的两种弹性参数。
9.如权利要求1所述的相控岩石物理模型指导下的储层预测方法，其特征在于，在所述步骤S8中，所述测井数据包括：实测纵波时差曲线、横波时差曲线、密度曲线；所述叠前地震数据包括三维叠前共反射点CRP地震道集。
10.如权利要求1所述的相控岩石物理模型指导下的储层预测方法，其特征在于，在所述步骤S9中，以两种弹性参数的值为坐标，将整个三维地震区内每个样点投到该坐标系中，提取位于储层部分的所有样点，上述样点在三维地震区内的位置即为储层位置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710289032.1A CN107179562B (zh) | 2017-04-27 | 2017-04-27 | 相控岩石物理模型指导下的储层预测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710289032.1A CN107179562B (zh) | 2017-04-27 | 2017-04-27 | 相控岩石物理模型指导下的储层预测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107179562A true CN107179562A (zh) | 2017-09-19 |
| CN107179562B CN107179562B (zh) | 2019-08-06 |
Family
ID=59832779
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710289032.1A Active CN107179562B (zh) | 2017-04-27 | 2017-04-27 | 相控岩石物理模型指导下的储层预测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107179562B (zh) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108491624A (zh) * | 2018-03-21 | 2018-09-04 | 中国石油化工股份有限公司 | 碳酸盐岩古暗河充填性预测方法 |
| CN110426752A (zh) * | 2019-08-20 | 2019-11-08 | 赛哲尔能源科技(北京)有限公司 | 一种基于岩石物理模型的储层参数反演方法与系统 |
| CN110489892A (zh) * | 2019-08-23 | 2019-11-22 | 中国石油大港油田勘探开发研究院 | 一种“三相耦合”储层构型刻画技术方法 |
| CN111366452A (zh) * | 2020-03-26 | 2020-07-03 | 北京科技大学 | 一种自蓄能岩体蓄能水平的测量方法 |
| CN111426616A (zh) * | 2020-02-20 | 2020-07-17 | 中国石油大学(北京) | 碳酸盐岩弹性性质与孔隙结构获取方法、装置及存储介质 |
| CN111556974A (zh) * | 2017-12-07 | 2020-08-18 | 沙特阿拉伯石油公司 | 利用岩石物理学和地震反演测绘储层的化学地层学特征 |
| CN113281825A (zh) * | 2020-02-20 | 2021-08-20 | 中国石油天然气集团有限公司 | 岩石物理模型构建方法及装置 |
| CN113960666A (zh) * | 2020-07-21 | 2022-01-21 | 中国石油天然气股份有限公司 | 一种基于叠前弹性参数的基岩气层识别方法及装置 |
| CN114910499A (zh) * | 2022-02-23 | 2022-08-16 | 河海大学 | 一种基于地层差异性成岩作用定量预测页岩气储层的方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB2321983B (en) * | 1996-12-13 | 2001-06-06 | Nigel Allister Anstey | Method of distinguishing geological sequences and their boundaries |
| CN1912655A (zh) * | 2005-08-10 | 2007-02-14 | 大庆油田有限责任公司 | 一种沉积相控制预测砂岩油藏储层属性的方法 |
| CN103412335A (zh) * | 2013-08-20 | 2013-11-27 | 成都晶石石油科技有限公司 | 一种利用地震物相体预测储层的方法 |
| CN103760081A (zh) * | 2013-12-31 | 2014-04-30 | 中国石油天然气股份有限公司 | 基于孔隙结构特征的碳酸盐岩储层的气藏预测方法及系统 |
-
2017
- 2017-04-27 CN CN201710289032.1A patent/CN107179562B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB2321983B (en) * | 1996-12-13 | 2001-06-06 | Nigel Allister Anstey | Method of distinguishing geological sequences and their boundaries |
| CN1912655A (zh) * | 2005-08-10 | 2007-02-14 | 大庆油田有限责任公司 | 一种沉积相控制预测砂岩油藏储层属性的方法 |
| CN103412335A (zh) * | 2013-08-20 | 2013-11-27 | 成都晶石石油科技有限公司 | 一种利用地震物相体预测储层的方法 |
| CN103760081A (zh) * | 2013-12-31 | 2014-04-30 | 中国石油天然气股份有限公司 | 基于孔隙结构特征的碳酸盐岩储层的气藏预测方法及系统 |
Non-Patent Citations (2)
| Title |
|---|
| 吴其林 等: "珠江口盆地H区块碳酸盐岩储层地震沉积学应用研究", 《沉积学报》 * |
| 赖锦 等: "基于岩石物理相约束的致密砂岩气储层渗透率解释建模", 《地球物理学进展》 * |
Cited By (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111556974A (zh) * | 2017-12-07 | 2020-08-18 | 沙特阿拉伯石油公司 | 利用岩石物理学和地震反演测绘储层的化学地层学特征 |
| CN111556974B (zh) * | 2017-12-07 | 2021-10-29 | 沙特阿拉伯石油公司 | 利用岩石物理学和地震反演测绘储层的化学地层学特征 |
| CN108491624A (zh) * | 2018-03-21 | 2018-09-04 | 中国石油化工股份有限公司 | 碳酸盐岩古暗河充填性预测方法 |
| CN110426752A (zh) * | 2019-08-20 | 2019-11-08 | 赛哲尔能源科技(北京)有限公司 | 一种基于岩石物理模型的储层参数反演方法与系统 |
| CN110489892B (zh) * | 2019-08-23 | 2022-04-05 | 中国石油天然气股份有限公司 | 一种“三相耦合”储层构型刻画技术方法 |
| CN110489892A (zh) * | 2019-08-23 | 2019-11-22 | 中国石油大港油田勘探开发研究院 | 一种“三相耦合”储层构型刻画技术方法 |
| CN111426616B (zh) * | 2020-02-20 | 2022-05-24 | 中国石油大学(北京) | 碳酸盐岩弹性性质与孔隙结构获取方法、装置及存储介质 |
| CN111426616A (zh) * | 2020-02-20 | 2020-07-17 | 中国石油大学(北京) | 碳酸盐岩弹性性质与孔隙结构获取方法、装置及存储介质 |
| CN113281825A (zh) * | 2020-02-20 | 2021-08-20 | 中国石油天然气集团有限公司 | 岩石物理模型构建方法及装置 |
| CN111366452B (zh) * | 2020-03-26 | 2021-01-29 | 北京科技大学 | 一种自蓄能岩体蓄能水平的测量方法 |
| CN111366452A (zh) * | 2020-03-26 | 2020-07-03 | 北京科技大学 | 一种自蓄能岩体蓄能水平的测量方法 |
| CN113960666A (zh) * | 2020-07-21 | 2022-01-21 | 中国石油天然气股份有限公司 | 一种基于叠前弹性参数的基岩气层识别方法及装置 |
| CN113960666B (zh) * | 2020-07-21 | 2024-01-30 | 中国石油天然气股份有限公司 | 一种基于叠前弹性参数的基岩气层识别方法及装置 |
| CN114910499A (zh) * | 2022-02-23 | 2022-08-16 | 河海大学 | 一种基于地层差异性成岩作用定量预测页岩气储层的方法 |
| CN114910499B (zh) * | 2022-02-23 | 2023-09-26 | 河海大学 | 一种基于地层差异性成岩作用定量预测页岩气储层的方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107179562B (zh) | 2019-08-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107179562B (zh) | 相控岩石物理模型指导下的储层预测方法 | |
| CN106054248B (zh) | 一种基于大面积致密储层地震岩石物理反演方法 | |
| CN103792581B (zh) | 基于岩石物理模型的页岩气储层地应力测井预测方法 | |
| CN102445708B (zh) | 三维等效富泥质砂岩速度预测模型 | |
| CN104502971B (zh) | 页岩储层总有机碳含量和脆性空间展布的检测方法 | |
| CN103135135B (zh) | 一种基于疏松砂岩模型进行烃类定量预测的方法和装置 | |
| CN108680951A (zh) | 一种基于地震信息判断煤层气富集沉积控制作用的方法 | |
| CN104950331A (zh) | 一种砂泥岩储层的孔隙度与泥质含量的地震预测方法 | |
| CN105182424A (zh) | 一种基于斑块饱和模型定量预测储层孔隙度的方法和装置 | |
| CN103217714A (zh) | 基于频率-尺度匹配的地震沉积学解释方法 | |
| CN108957526B (zh) | 获得裂缝融合数据体的方法 | |
| CN107045143A (zh) | 一种预测裂缝发育的方法及装置 | |
| Chen et al. | The application study on the multi-scales integrated prediction method to fractured reservoir description | |
| CN107728202A (zh) | 致密砂岩储层精细描述方法及系统 | |
| CN106597543A (zh) | 一种地层沉积相划分方法 | |
| CN105607120A (zh) | 一种基于时移测井的地震相约束构建初始模型的方法 | |
| Lu et al. | Prestack inversion identification of organic reef gas reservoirs of Permian Changxing Formation in Damaoping area, Sichuan Basin, SW China | |
| CN104820237B (zh) | 计算地层的方差体的方法 | |
| CN107219563B (zh) | 在沉积相控制下进行储层avo分析的方法 | |
| Davis | Velocity variations around Leduc reefs, Alberta | |
| CN109143397A (zh) | 碳酸盐岩储层缝洞充填物识别方法及系统 | |
| Zhonghua et al. | Application of VSP survey for shale gas exploration | |
| Shoaib et al. | Ratana field, potwar fold belt, northern pakistan: high intensity fracture zones related to major thrust faults as revealed by seismic fracture prediction | |
| Schuelke | Overview of seismic attribute analysis in shale plays | |
| Yushendri | Quantitative Seismic Interpretation of the Devonian Woodford Shale in Anadarko Basin, Oklahoma |
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 | ||
| CP01 | Change in the name or title of a patent holder | ||
| CP01 | Change in the name or title of a patent holder |
Address after: 100094 Room 401, building 4, courtyard 3, fengxiu Middle Road, Haidian District, Beijing Patentee after: New JCM Group Co.,Ltd. Address before: 100094 Room 401, building 4, courtyard 3, fengxiu Middle Road, Haidian District, Beijing Patentee before: LANDOCEAN ENERGY SERVICES CO.,LTD. |