CN108227033B - 砂岩储层的孔隙度的预测方法、装置及计算机存储介质 - Google Patents
砂岩储层的孔隙度的预测方法、装置及计算机存储介质 Download PDFInfo
- Publication number
- CN108227033B CN108227033B CN201810011059.9A CN201810011059A CN108227033B CN 108227033 B CN108227033 B CN 108227033B CN 201810011059 A CN201810011059 A CN 201810011059A CN 108227033 B CN108227033 B CN 108227033B
- Authority
- CN
- China
- Prior art keywords
- porosity
- action
- reservoir
- under
- target reservoir
- 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 84
- 238000009933 burial Methods 0.000 claims abstract description 189
- 230000009471 action Effects 0.000 claims abstract description 187
- 230000009467 reduction Effects 0.000 claims abstract description 104
- 239000011148 porous material Substances 0.000 claims abstract description 71
- 239000013049 sediment Substances 0.000 claims abstract description 62
- 238000004364 calculation method Methods 0.000 claims abstract description 28
- 230000008021 deposition Effects 0.000 claims abstract description 27
- 239000002245 particle Substances 0.000 claims description 72
- 230000015572 biosynthetic process Effects 0.000 claims description 34
- 230000000694 effects Effects 0.000 claims description 13
- 238000004422 calculation algorithm Methods 0.000 claims description 11
- 238000011038 discontinuous diafiltration by volume reduction Methods 0.000 claims description 5
- 239000000126 substance Substances 0.000 claims description 4
- 238000004590 computer program Methods 0.000 claims description 3
- 239000010802 sludge Substances 0.000 claims description 3
- 238000000151 deposition Methods 0.000 description 20
- 238000001125 extrusion Methods 0.000 description 20
- 230000008859 change Effects 0.000 description 18
- 238000005056 compaction Methods 0.000 description 18
- 230000008569 process Effects 0.000 description 17
- 239000011435 rock Substances 0.000 description 16
- 230000000875 corresponding effect Effects 0.000 description 14
- 230000006870 function Effects 0.000 description 10
- 238000004458 analytical method Methods 0.000 description 9
- 239000000463 material Substances 0.000 description 8
- 238000007405 data analysis Methods 0.000 description 7
- 239000000523 sample Substances 0.000 description 7
- 238000011161 development Methods 0.000 description 6
- 238000002474 experimental method Methods 0.000 description 6
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N Silicium dioxide Chemical compound O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 5
- 241000030538 Thecla Species 0.000 description 4
- 239000004568 cement Substances 0.000 description 4
- 238000012417 linear regression Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 230000007797 corrosion Effects 0.000 description 3
- 238000005260 corrosion Methods 0.000 description 3
- 230000007547 defect Effects 0.000 description 3
- 230000001419 dependent effect Effects 0.000 description 3
- 229910052500 inorganic mineral Inorganic materials 0.000 description 3
- 239000011707 mineral Substances 0.000 description 3
- 230000035699 permeability Effects 0.000 description 3
- 239000010453 quartz Substances 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 230000002596 correlated effect Effects 0.000 description 2
- 238000005553 drilling Methods 0.000 description 2
- 230000003628 erosive effect Effects 0.000 description 2
- 239000012634 fragment Substances 0.000 description 2
- 230000007774 longterm Effects 0.000 description 2
- 239000010445 mica Substances 0.000 description 2
- 229910052618 mica group Inorganic materials 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- 239000011044 quartzite Substances 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 239000012798 spherical particle Substances 0.000 description 2
- 208000035126 Facies Diseases 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 239000007771 core particle Substances 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000003066 decision tree Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000003921 particle size analysis Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V9/00—Prospecting or detecting by methods not provided for in groups G01V1/00 - G01V8/00
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本说明书提供砂岩储层的孔隙度的预测方法、装置及计算机存储介质，包括：根据目标储层的原始沉积体的参数，计算目标储层的原始孔隙度；其中，所述原始沉积体用于表示在沉积末期，埋藏初期所形成的沉积物；根据目标储层的埋深，计算在指定埋深作用下目标储层的孔隙度；根据目标储层的构造应力，计算在指定构造应力作用下目标储层的孔隙度减孔量；根据实测砂岩储层的孔隙度减孔量计算埋深作用下的减孔系数和构造应力作用下的减孔系数；根据所述目标储层的原始孔隙度、埋深作用下目标储层的孔隙度、构造应力作用下目标储层的孔隙度减孔量和埋深作用下的减孔系数以及构造应力作用下的减孔系数预测所述目标储层的孔隙度。提高了孔隙度计算精度。
Description
技术领域
本说明书涉及石油天然气勘探领域，特别涉及砂岩储层的孔隙度的预测方法、装置及计算机存储介质。
背景技术
孔隙度是指岩样中所有孔隙空间体积之和与该岩样体积的比值，储层的孔隙度越大，说明岩石中孔隙空间越大，孔隙不仅能够储存油气，而且可以允许油气在其中渗滤，孔隙度是储层表征的重要参数，储层孔隙度的大小直接代表了储层储集性能的好坏。
致密砂岩储层是非常重要的油气储层，属于非常规油气储层，按照国内为对超深层储层的定义，超过3500-4000m的储层为深层储层。二十世纪七十年代，美国联邦能源管理委员会将储层基质渗透率小于0.1×10-3μm2的气藏定义为致密气藏。
在我国，塔里木盆地白垩系巴什基奇克组砂岩储层埋深达6000-7000m，样品渗透率小于0.1×10-3μm2，属超深层致密砂岩储层。近年来我国非常规油气资源勘探开发取得了重大突破，新增探明储量逐年增加，占我国油气新增探明储量的60％以上。因此研究塔里木盆地白垩系巴什基奇克组超深层致密砂岩储层的储集能力，预测计算其孔隙度对指导塔里木盆地白垩系巴什基奇克组砂岩储层的油气资源勘探与开发具有十分重要的意义。
在实现本说明书的过程中，发明人发现现有技术中存在如下问题。
目前，国内外对储层孔隙度的计算预测方法大致可以分为深度指数法以及基于深度指数法的相应孔隙度预测方法；基于成岩物理实验的孔隙度预测方法；基于储层孔隙演化史孔隙度模拟方法；基于地震相的储层孔隙度预测方法以及基于局部的决策树算法的孔隙度预测方法。上述方法对储层孔隙度的预测大多基于储层成岩演化过程或是基于储层的地震属性，鉴于影响储层孔隙度发育因素的多样性，上述预测方法存在不足。
发明内容
本说明书实施方式提供砂岩储层的孔隙度的预测方法、装置及计算机存储介质。通过计算得到目标储层原始孔隙度、埋藏压实作用下目标储层的孔隙度、构造挤压作用下目标储层的孔隙度以及埋藏压实作用下的减孔系数和和构造挤压作用下的减孔系数，从而预测得到目标储层的孔隙度。
本说明书实施方式提供一种砂岩储层的孔隙度的预测方法。包括，根据目标储层的原始沉积体的参数，计算目标储层的原始孔隙度；其中，所述原始沉积体用于表示在沉积末期，埋藏初期所形成的沉积物；所述原始孔隙度用于表示所述目标储层在沉积末期，埋藏初期时的孔隙度；根据目标储层的埋深，计算在指定埋深作用下目标储层的孔隙度；根据目标储层的构造应力，计算在指定构造应力作用下目标储层的孔隙度减孔量；根据实测砂岩储层的孔隙度减孔量计算埋深作用下的减孔系数和构造应力作用下的减孔系数；根据所述目标储层的原始孔隙度、埋深作用下目标储层的孔隙度、构造应力作用下目标储层的孔隙度减孔量和埋深作用下的减孔系数以及构造应力作用下的减孔系数预测所述目标储层的孔隙度。
本说明书实施方式提供一种砂岩储层的孔隙度的预测装置。包括，原始孔隙度计算模块；根据目标储层的原始沉积体的参数，计算目标储层的原始孔隙度；其中，所述原始沉积体用于表示在沉积末期，埋藏初期所形成的沉积物；所述原始孔隙度用于表示所述目标储层在沉积末期，埋藏初期时的孔隙度；埋深作用下孔隙度计算模块；用于根据目标储层的埋深，计算在指定埋深作用下目标储层的孔隙度；构造应力作用下孔隙度减孔量计算模块；用于根据目标储层的构造应力，计算在指定构造应力作用下目标储层的孔隙度减孔量；减孔系数计算模块；根据实测砂岩储层的孔隙度减孔量计算埋深作用下的减孔系数和构造应力作用下的减孔系数；孔隙度计算模块；用于根据所述目标储层的原始孔隙度、埋深作用下目标储层的孔隙度、构造应力作用下目标储层的孔隙度减孔量和埋深作用下的减孔系数以及构造应力作用下的减孔系数预测所述目标储层的孔隙度。
本说明书实施方式提供一种计算机存储介质。所述计算机存储介质存储有计算机程序指令，在所述计算机程序指令被执行时实现：根据目标储层的原始沉积体的参数，计算目标储层的原始孔隙度；其中，所述原始沉积体用于表示在沉积末期，埋藏初期所形成的沉积物；所述原始孔隙度用于表示所述目标储层在沉积末期，埋藏初期时的孔隙度；根据目标储层的埋深，计算在指定埋深作用下目标储层的孔隙度；根据目标储层的构造应力，计算在指定构造应力作用下目标储层的孔隙度减孔量；根据实测砂岩储层的孔隙度减孔量计算埋深作用下的减孔系数和构造应力作用下的减孔系数；根据所述目标储层的原始孔隙度、埋深作用下目标储层的孔隙度、构造应力作用下目标储层的孔隙度减孔量和埋深作用下的减孔系数以及构造应力作用下的减孔系数预测所述目标储层的孔隙度。
本说明书的有益效果在于基于区域大量可靠的实验资料、砂岩动力成岩理论、数值分析和地质统计学方法，通过计算得到目标储层原始孔隙度、埋深作用下目标储层的孔隙度、构造应力作用下目标储层的孔隙度以及埋深作用下的减孔系数和构造应力作用下的减孔系数，综合计算目标储层的孔隙度。本说明书实施方式考虑了砂岩原始沉积体对储层原生孔隙及其演化的影响，还考虑了后期埋藏压实以及构造挤压作用对储层孔隙度的影响，并计算埋藏压实与构造挤压作用下的减孔系数，最终应用预测孔隙度，解决了现有技术中砂岩储层孔隙度预测过程中因为没有考虑到孔隙度发育因素多样性所带来的不足。
附图说明
为了更清楚地说明本说明书实施方式或现有技术中的技术方案，下面将对实施方式或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本说明书中记载的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本说明书实施方式中提供的砂岩储层的孔隙度的预测方法的流程图；
图2为本说明书实施方式中提供的砂岩储层的孔隙度的预测装置的结构框图；
图3为本说明书实施方式中提供的砂岩储层孔隙度与平均粒径拟合图；
图4为本说明书实施方式中提供的砂岩储层孔隙度与泥质含量拟合图；
图5为本说明书实施方式中提供的砂岩储层孔隙度与分选系数拟合图；
图6为本说明书实施方式中提供的砂岩储层孔隙度与刚性颗粒含量拟合图；
图7为本说明书实施方式中提供的实验条件下砂岩储层孔隙度随埋深变化图；
图8为本说明书实施方式中提供的砂岩储层孔隙度随埋深变化图；
图9为本说明书实施方式中提供的砂岩储层孔隙度随构造应力变化图；
图10为本说明书实施方式中提供的砂岩储层的埋藏史线、孔隙史线、地温史线。
具体实施方式
为了使本技术领域的人员更好地理解本说明书中的技术方案，下面将结合本说明书实施方式中的附图，对本说明书实施方式中的技术方案进行清楚、完整地描述，显然，所描述的实施方式仅仅是本说明书一部分实施方式，而不是全部的实施方式。基于本说明书中的实施方式，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施方式，都应当属于本说明书保护的范围。
在一个场景实施方式中，已有库车坳陷秋探1井、东秋8井白垩系巴什基奇克组储集层孔隙度，通过使用本说明书提供的实施方式对库车坳陷秋探1井、东秋8井白垩系巴什基奇克组储集层的孔隙度进行计算，并与已有的储集层孔隙度进行比较，验证孔隙度计算精度。
在本场景实施方式中，已有克拉2气田白垩系砂岩储层的地质资料，包括岩屑、岩心粒度分析结果，岩石矿物成分、埋深等资料。经过资料分析，该区域在埋深3600-4000米有残余原生孔隙型储层，胶结物含量4％-9％，构造应力为30-40MPa。在埋藏条件下，如图3和表1所示，该残余原生孔隙型储层与其原始沉积体的平均粒径、分选系数、泥质含量呈明显的负相关与刚性颗粒含量呈正相关。所述原始沉积体可以包括砂岩原始沉积时期的沉积物粒径、分选、泥质杂基、刚性颗粒(石英质)、塑性颗粒(软岩屑，云母等)等。根据原始沉积体可以确定砂岩的原始物性，可以根据原始沉积提来确定砂岩刚刚沉积之后孔隙度。
表1克拉2气田白垩系砂岩储层砂岩孔隙度与原始沉积体参数统计表
| 井号 | 深度/m | 实测岩心孔隙度/％ | S<sub>o</sub> | D/mm | V<sub>sh</sub>/％ | V<sub>Q</sub>/％ | 胶结物含量/％ |
| 克拉201 | 3 667.42 | 15.43 | 0.721 0 | 2.446 3 | 2.0 | 60 | 7.3 |
| 克拉201 | 3 675.76 | 16.86 | 0.673 0 | 3.544 3 | 5.0 | 65 | 3.0 |
| 克拉201 | 3 742.04 | 15.64 | 1.460 0 | 2.845 0 | 2.0 | 61 | 4.5 |
| 克拉201 | 3 743.04 | 16.40 | 1.314 0 | 2.665 0 | 2.0 | 63 | 6.9 |
| 克拉201 | 3 743.47 | 16.77 | 1.324 0 | 2.571 0 | 4.0 | 64 | 3.0 |
| 克拉201 | 3 743.89 | 16.04 | 1.360 0 | 2.680 0 | 4.0 | 62 | 4.0 |
| 克拉201 | 3 744.27 | 14.75 | 1.388 0 | 2.768 0 | 4.0 | 58 | 4.5 |
| 克拉201 | 3 745.15 | 14.51 | 1.597 0 | 3.353 0 | 4.0 | 59 | 9.0 |
| 克拉201 | 3 797.68 | 12.04 | 1.722 0 | 4.286 0 | 5.0 | 54 | 6.5 |
| 克拉201 | 3 800.21 | 10.40 | 1.851 0 | 4.458 0 | 7.0 | 50 | 10.5 |
| 克拉201 | 3 801.13 | 14.15 | 1.719 0 | 4.417 0 | 7.0 | 59 | 3.5 |
| 克拉201 | 3 801.79 | 10.61 | 1.797 0 | 4.105 0 | 5.0 | 50 | 6.0 |
| 克拉201 | 3 802.30 | 10.81 | 1.859 0 | 4.202 0 | 5.0 | 51 | 9.0 |
| 克拉201 | 3 802.55 | 8.65 | 4.585 0 | 4.495 0 | 11.0 | 49 | 6.0 |
| 克拉201 | 3 827.81 | 14.44 | 1.678 0 | 4.301 0 | 6.0 | 61 | 7.5 |
| 克拉201 | 3 828.29 | 15.65 | 1.523 0 | 4.250 0 | 6.0 | 60 | 5.5 |
| 克拉201 | 3 828.61 | 14.24 | 1.553 0 | 3.452 0 | 3.0 | 58 | 5.5 |
| 克拉201 | 3 830.77 | 16.21 | 1.566 0 | 2.451 0 | 3.0 | 61 | 5.0 |
| 克拉201 | 3 857.39 | 13.00 | 1.839 0 | 4.490 0 | 5.0 | 57 | 6.0 |
| 克拉201 | 3 858.18 | 12.12 | 2.282 0 | 4.593 0 | 6.5 | 56 | 7.5 |
| 克拉201 | 3 858.49 | 12.51 | 3.165 0 | 4.708 0 | 4.0 | 56 | 7.0 |
| 克拉201 | 3 858.62 | 11.30 | 2.448 0 | 4.788 0 | 3.0 | 52 | 8.5 |
| 克拉201 | 3 858.89 | 11.41 | 4.259 0 | 4.697 0 | 3.0 | 52 | 8.5 |
| 克拉201 | 3 861.62 | 12.40 | 3.726 0 | 4.792 0 | 2.0 | 58 | 7.5 |
| 克拉201 | 3 923.59 | 10.97 | 2.327 0 | 3.952 0 | 4.0 | 55 | 7.0 |
| 克拉2 | 3 575.60 | 9.68 | 2.338 7 | 1.856 8 | 13.0 | 55 | 6.0 |
| 克拉2 | 3 621.95 | 16.28 | 1.745 1 | 1.805 4 | 8.0 | 61 | 8.0 |
| 克拉2 | 3 562.49 | 17.88 | 1.851 7 | 1.318 0 | 5.0 | 69 | 4.5 |
| 克拉2 | 3 741.46 | 16.36 | 1.853 5 | 0.787 1 | 3.0 | 62 | 3.0 |
| 克拉2 | 3 746.03 | 18.81 | 1.507 1 | 2.182 4 | 3.0 | 65 | 3.5 |
| 克拉2 | 3 796.06 | 17.21 | 1.631 6 | 1.677 8 | 2.0 | 63 | 3.0 |
| 克拉2 | 3 827.76 | 13.63 | 1.676 4 | 1.233 5 | 10.0 | 57 | 7.0 |
在本场景实施方式中，接收如表1所示的原始沉积体的地质参数，可以包括原始沉积体的平均粒径D、分选系数S0、泥质含量Vsh和刚性颗粒含量VQ，所述刚性颗粒含量包括石英颗粒和硅质岩屑。基于如表1所示的地质数据，建立多元线性回归方程，可以求得回归系数，多元线性回归方程可以为：
计算原始孔隙度。其中，
表示原始孔隙度，So表示分选系数，D表示平均粒径、Vsh表示泥质含量、VQ表示刚性颗粒含量。
在本场景应用中，可以通过地质资料分析，获取库车坳陷秋探1井、东秋8井白垩系巴什基奇克组储集层原始沉积体的平均粒径、分选系数、泥质含量、和刚性颗粒含量，带入多元线性回归方程得到库车坳陷秋探1井、东秋8井白垩系巴什基奇克组储集层的原始孔隙度。
在本场景实施例中，由于克拉2气田白垩系砂岩储层经历了长期浅埋短期深埋的过程，因此进行长期浅埋短期深埋条件下砂岩孔隙度与埋深变化关系实验，通过砂岩埋藏压实实验获得，埋藏的过程实质是温压的增加过程，设定地质条件为长时间的低温(<40℃，<30Mpa)，短时间的高温(>100℃，>100Mpa)，从而得到各个阶段的砂岩储层孔隙度。如表2所示，接收克拉2气田白垩系砂岩储层砂岩孔隙度与深度的统计关系数据与实验结果数据(如图7所示)，发现储层砂岩孔隙度与埋深深度之间满足指数关系。另外，由于地质情况是复杂的，岩石类型也是多样的(包括细砂岩、中砂岩、粗砂岩，石英含量也不同，有<50％、50-65％、>65％等，根据不同类型砂岩的地质数据和储层砂岩孔隙度与埋深的数据，因此根据表2建立指数函数回归方程，按照砂岩类型的不同，如表3和图8所示，计算得到目标储层在埋深作用下的孔隙度。
表2白垩系巴什基奇克组砂岩埋深与砂岩孔隙度数据统计表
表3白垩系巴什基奇克组不同类型储层砂岩孔隙度与埋深关系
注：泥质含量1％～5％；胶结物含量2％～5％；溶蚀量0～2％
在本场景实施例中，接收在构造应力作用下白垩系巴什基奇克组砂岩储层砂岩孔隙度减孔量数据，如表4所示。进行砂岩孔隙度与构造应力变化关系实验，接收砂岩孔隙度减孔量与构造应力变化关系实验数据。根据实测数据和实验数据发现砂岩孔隙度和构造应力大小之间满足指数关系。
在本场景实施例中，根据表4数据建立指数函数回归方程，按照砂岩类型的不同，如图9所示，计算得到计算的目标储层在构造应力作用下的减孔量。
表4白垩系砂岩储集层最大古构造应力与构造挤压减孔量统计表
请参阅图10，在本场景实施例中，通过盆地模拟法获取地史数据，所述地史数据可以反映地质属性随时间的变化。通过盆地模拟软件Petromod，输入地层厚度、压实率、剥蚀量、地质时间生成储层埋藏史、热演化史、孔隙演化史曲线，；输入现今地层温度、区域地温梯度、大地热流值、地层厚度通过盆地模拟软件Petromod生成地层热史曲线；输入现今储层孔隙度、原始孔隙度、埋藏过程地层厚度、构造应力及胶结物含量、溶蚀增孔量通过盆模软件Petromod生成储层孔隙演化史曲线。
在本场景实施例中，根据所要计算的目标储层的埋深范围，根据所述埋藏史线、热演化史线和孔隙度演化史线确定在该埋深条件下的总减孔量。
在本场景实施例中，在4000米埋深之后总减孔量为16％，根据地质资料获取克拉2气田白垩系砂岩储层的最大古构造应力，计算构造应力作用下克拉2气田白垩系砂岩储层的孔隙度减孔量为11％，埋深作用下的减孔量为16％-11％＝5％，因此埋藏压实减孔系数为5％÷16％＝0.31，构造挤压减孔系数为11％÷16％＝0.69。
在本场景实施例中，根据公式，
计算目标储层的孔隙度，其中其中，
表示目标储层孔隙度，
表示原始沉积体的砂岩孔隙度，
表示构造挤压作用下目标储层的减孔量，
构造挤压作用下目标储层的减孔量，
表示埋藏压实作用下目标储层的孔隙度。
在本情景实施方式中，区域资料表明秋探1井白垩系巴什基奇克组刚性颗粒含量一般为55-60％，处于弱构造挤压带，最大古构造应力为35MPa，以分选性好的细砂岩为主，埋深6100-6300m；东秋8井白垩系巴什基奇克组刚性颗粒含量一般为65-70％，处于弱构造挤压带，最大古构造应力为30MPa，以分选性好的细砂岩为主，埋深5100-5300m，如表5所示，计算两口井孔隙度分别为11.33％，16.46％，与钻井实测/测井孔隙度误差率不足6％。
表5库车坳陷秋探1井、东秋8井白垩系巴什基奇克组储集层孔隙度综合预测表
请参阅图1，本说明书实施方式提供一种储层的孔隙度的预测方法。
在本实施方式中，所述储层的孔隙度的预测方法可以应用于电子设备。使得所述电子设备通过运行该孔隙度计算方法可以实现综合孔隙度发育多样性的特点计算孔隙度的效果。
所述储层的孔隙度的预测方法具体可以包括如下步骤。
步骤S10:根据目标储层的原始沉积体的参数，计算目标储层的原始孔隙度；其中，所述原始沉积体用于表示在沉积末期，埋藏初期所形成的沉积物；所述原始孔隙度用于表示所述目标储层在沉积末期，埋藏初期时的孔隙度。
在本实施方式中，所述原始沉积体可以包括砂岩目标储层在刚刚沉积之后，浅埋藏时期的沉积物粒径、分选、泥质杂基、刚性颗粒(石英质)、塑性颗粒(软岩屑，云母等)等。所述原始沉积体可以处于埋深3600～4000m的残余原生孔隙型储集层，胶结物含量4％～9％，构造应力为30～40MPa。可以根据原始沉积体可以确定砂岩的原始物性，可以根据原始沉积提来确定砂岩刚刚沉积之后孔隙度。
在本实施方式中，所述原始沉积体参数可以是反映原始沉积体属性的参数，可以包括岩屑、岩心粒度分析结果，岩石矿物成分、埋深等资料。
在本实施方式中，实测砂岩储层原始孔隙度是指通过区域钻井和地质资料分析，得到剩余原生孔隙型储层，测得该剩余原生孔隙型储层的孔隙度。
在本实施方式中，所述原始孔隙度可以包括，目标储层在刚刚沉积后的孔隙度。目标储层在沉积过程中，可以在目标储层所处区域有残余原生孔隙型储层，可以代表刚刚沉积期末，浅埋藏时期储层的孔隙度。
在本实施方式中，所述原始沉积体的参数可以包括所述原始沉积体的平均粒径、分选系数、泥质含量以及刚性颗粒含量等。可以根据实验分析资料和测井资料分析等方式得到所述原始沉积体的参数。可以包括岩屑、岩心粒度分析，岩石矿物成分、埋深等资料分析。
步骤S12：根据目标储层的埋深，计算在指定埋深作用下目标储层的孔隙度。
在本实施方式中，所述目标储层的埋深可以包括目标储层现今所在的位置与地表平面或者设定的基准平面之间的深度关系。具体地，例如埋深8000m的储层可以是指距离地表平面8000m的地下储层。
在本实施方式中，所述指定埋深作用下目标储层的孔隙度可以包括，在目标储层下移到现今的位置时，目标储层的孔隙度。所述目标储层在沉积完成后可以经历埋藏时期，在埋藏时期，所述目标储层上覆的沉积物可以不断加厚，所述目标储层的位置相对于地表平面不断下移，在埋藏的过程中，随着上覆沉积物的增多，在重荷压力下，目标储层孔隙中的水可以被排除，孔隙度可以减少。
在本实施方式中，计算在指定埋深作用下目标储层的孔隙度可以包括，根据已有的孔隙度与埋深数据，构建孔隙度随埋深变化的关系，将指定埋深作用代入孔隙度随埋深变化的关系中，求得目标储层在指定埋深作用下的孔隙度。
步骤S14：根据目标储层的构造应力，计算在指定构造应力作用下目标储层的孔隙度减孔量。
在本实施方式中，所述构造应力可以包括，在地壳构造运动作用力的影响下，地壳中所产生的应力，构造应力在空间上可以是垂直的也可以是水平的，在构造应力的作用下，所述目标储层的孔隙可以被挤压，可以造成孔隙度减少。所述构造应力可以通过岩心样品进行声发射法实验测试获得。
在本实施方式中，减孔量是指在构造应力的作用下，孔隙度减少的量，在本实施方式中，在构造应力的作用下，减孔因素可以有很多，可以包括构造因素、埋藏因素等，构造应力与孔隙度之间不是直接的线性关系，而与减孔量之间是直接的线性关系，因此求取减孔量与构造应力之间的关系。
在本实施方式中，计算在指定构造应力作用下目标储层的孔隙度减孔量可以包括，获取已有的孔隙度与构造应力数据，并进行构造应力与孔隙度关系的实验模拟得到孔隙度减孔量随构造应力变化的关系，将目标储层的构造应力代入孔隙度减孔量随构造应力变化的关系可以得到所述减孔量。
步骤S16：根据实测砂岩储层的孔隙度减孔量计算埋深作用下的减孔系数和构造应力作用下的减孔系数。
在本实施方式中，实测砂岩储层可以包括，已经经过实际勘探得到的其地质属性的储层。实测砂岩储层的地质属性可以包括孔隙度、埋深、构造应力、地层厚度、压实率、剥蚀率、胶结物含量等。
在本实施方式中，所述埋深作用下的减孔系数可以包括，所述目标储层在经历沉积、埋藏、挤压等作用下，孔隙度可以不断的变化，在孔隙度变化的过程中，埋深作用可以不是唯一影响孔隙度变化的因素，埋深在孔隙度变化的过程中所起到的作用可以进行量化的表示，表示埋深在孔隙度变化过程中所起到的作用的系数可以是埋深作用下的减孔系数。
在本实施方式中，所述构造应力作用下的减孔系数可以包括，所述目标储层在经理沉积、埋藏、挤压等作用下，孔隙度可以不断的变化，在孔隙度变化的过程中，构造应力作用可以不是唯一影响孔隙度变化的因素，构造应力在孔隙度变化的过程中所起到的作用可以进行量化的表示，表示构造应力在孔隙度变化过程中所起到的作用的系数可以是构造挤压作用下的减孔系数。
在本实施方式中，所述实测砂岩储层孔隙度的减孔量可以包括实测砂岩储层在埋深作用下的孔隙度减孔量、实测砂岩储层在构造应力作用下的孔隙度减孔量、实测砂岩储层在沉积末期到如今孔隙度的总减孔量。
在本实施方式中，计算所述埋深作用下的减孔系数和构造应力作用下的减孔系数可以包括，将实测砂岩储层在沉积末期到如今孔隙度的总减孔量减去在构造应力作用下的减孔量得到实测砂岩储层在埋深作用下的减孔量，分别将实测砂岩储层在埋深作用下的减孔量和实测砂岩储层在构造应力作用下的减孔量除以实测砂岩储层在沉积末期到如今孔隙度的总减孔量得到所述埋深作用下的减孔系数和构造应力作用下的减孔系数。
步骤S18：根据所述目标储层的原始孔隙度、埋深作用下目标储层的孔隙度、构造应力作用下目标储层的孔隙度减孔量和埋深作用下的减孔系数以及构造应力作用下的减孔系数预测所述目标储层的孔隙度。
在本实施方式中预测所述目标储层的孔隙度可以包括，将所述目标储层的原始孔隙度减去所述埋深作用下目标储层的孔隙度得到所述目标储层在埋深作用下的孔隙度减孔量。将所述目标储层的原始孔隙度减去埋深作用下目标储层的孔隙度减孔量和所述埋深作用下的减孔系数的乘积再减去构造应力作用下目标储层的孔隙度减孔量和所述构造应力作用下的减孔量得到所述目标储层的孔隙度。
在本说明书实施方式中，根据目标储层的原始孔隙度、在埋深作用下目标储层的孔隙度、在构造应力作用下目标储层的孔隙度减孔量以及埋深作用下的减孔系数和构造挤压作用下的减孔系数预测目标储层的孔隙度综合考虑了埋藏压实因素、构造挤压因素以及地质演化因素计算目标储层的孔隙度，使孔隙度的预测计算更加准确。
在一个实施方式中，所述地层的参数至少包括：所述原始沉积体的平均粒径、分选系数、泥质含量以及刚性颗粒含量数据。
在本实施方式中，所述平均粒径可以包括对于一个由大小和形状不相同的粒子组成的实际粒子群，与一个由均一的球形粒子组成的假想粒子群相比，如果两者的粒径全长相同，可以称此球形的粒子的直径为实际粒子群的平均粒径，平均粒径可以通过多种方法计算，可以为算数平均法、几何平均法、调和平均法、体积平均法及质量平均法等。
在本实施方式中，所述分选系数可以包括粒度累计曲线上，25％和75％处所对应的颗粒直径的比值，可以表示碎屑沉积物分选性的一种参考，碎屑大小均匀的分选性好，大小混杂的分选性差。
在本实施方式中，所述泥质含量可以包括泥质的体积占岩石总体积的比，泥质可以包括颗粒直径小于0.01mm的碎屑物质。确定泥质含量的方法可以为自然伽马法、自然电位法、中子法、交会图法等。
在本实施方式中，所述刚性颗粒含量可以包括石英颗粒、硅质岩屑等占岩石总颗粒构成的比值。
在本实施方式中，如图3至图6所示，砂岩孔隙度与其原始沉积体的平均粒径、分选系数、泥质含量呈明显的负相关，与刚性颗粒含量呈正相关，这4个参数可以是储层孔隙度的敏感因素，具有代表性且均能量化参数，因此可以优选这4个参数参与多元线性回归运算来预测计算原始沉积体的孔隙度。
在一个实施方式中，计算所述原始孔隙度包括：根据实测砂岩储层原始沉积体的孔隙度、平均粒径、分选系数、泥质含量与刚性颗粒含量通过回归算法得到所述平均粒径、分选系数、泥质含量与刚性颗粒含量对应的回归系数；根据所述与目标储层关联的指定地层的平均粒径、分选系数、泥质含量与刚性颗粒含量以及所述回归系数计算所述原始孔隙度。
在本实施方式中，所述回归算法可以包括，以所述实测砂岩储层原始沉积体的孔隙度为因变量，以所述实测砂岩储层原始沉积体的平均粒径、分选系数、泥质含量以及刚性颗粒含量为自变量，建立多元线性回归方程，已知实测砂岩储层原始沉积体的孔隙度、平均粒径、分选系数、泥质含量以及刚性颗粒含量可以得到平均粒径、分选系数、泥质含量以及刚性颗粒含量各自的回归系数，得到回归方程。
在本实施方式中，计算所述目标储层原始沉积体的原始孔隙度可以包括，获取目标储层原始沉积体的平均粒径、分选系数、泥质含量与刚性颗粒含量，代入根据回归算法所得到的回归方程，可以得到所述目标储层原始沉积体的原始孔隙度。
在本实施方式中，所述目标储层原始沉积体的原始孔隙度可以根据公式，
计算得出；其中，
表示原始沉积体的砂岩孔隙度，So表示分选系数，D表示平均粒径、Vsh表示泥质含量、VQ表示刚性颗粒含量，a、b、c、d、e分别表示回归系数。
在本实施方式中，通过已有的实测砂岩储层孔隙度的原始沉积体的孔隙度、平均粒径、分选系数、泥质含量以及刚性颗粒含量构建回归方程，并将目标储层的原始沉积体的孔隙度、平均粒径、分选系数、泥质含量代入回归方程从而得到目标储层的原始沉积体的孔隙度的方式，利用统计和数理分析的方法得到原始孔隙度与相关参数之间的关系，提高原始孔隙度的计算精度。
在一个实施方式中，计算所述指定埋深作用下目标储层的孔隙度，包括：根据实测砂岩储层孔隙度和埋深数据，通过回归算法计算回归系数；根据所述目标储层的埋深和所述回归系数计算在指定埋深作用下目标储层的孔隙度。
在本实施方式中，所述实测砂岩储层孔隙度和埋深数据可以包括，已经实际勘探的区域，储层埋深数据以及与埋深相对应的砂岩储层孔隙度。
在本实施方式中，所述回归算法可以包括，根据所述实测砂岩储层孔隙度和埋深数据可以得到所述实测砂岩储层孔隙度和埋深数据满足指数关系。可以以埋深为自变量，实测砂岩储层的孔隙度为因变量构建指数回归方程，计算回归系数，得到指数回归方程。
在本实施方式中，计算指定埋深作用下目标储层的孔隙度可以包括，将目标储层的指定深度代入已经构建的指数回归方程中，得到目标储层在指定埋深作用的埋深作用下的孔隙度。
在本实施方式中，所述指定埋深作用下目标储层的孔隙度可以根据公式：：
计算得出；其中f、g表示回归系数，H表示所述指定埋深。
在本实施方式中，通过实测砂岩储层孔隙度和埋深数据，建立回归方程，并将目标储层的埋深数据代入回归方程计算在目标储层在指定埋深作用下的孔隙度的方式。利用统计和数理分析的方法得到孔隙度与埋深之间的关系，提高所述指定埋深作用下目标储层的孔隙度计算精度。
在一个实施方式中，根据不同砂岩类型计算不同砂岩类型对应的回归系数。
在本实施方式中，所述不同砂岩类型可以包括，实测砂岩储层的性质可以包括砂岩岩性、分选性、泥质含量等，将性质相同或者相近的实测砂岩储层分为一类，具体地，例如将砂岩类型为细砂岩、分选性好、泥质含量小于3％的实测砂岩储层分为一类；将砂岩类型为中砂岩、分选性好、泥质含量为大于3％小于5％的分为一类等。实测砂岩储层的性质可以通过岩心分析、测井分析等方法得到
在本实施方式中，计算不同砂岩类型对应的回归系数可以包括，按照砂岩类型为实测砂岩储层进行分类。具体地，可以将储层为细砂岩、分选性好、刚性颗粒含量65％到75％的分为一类；将中砂岩、分选性好，刚性颗粒含量为75％到85％的分为一类。通过回归方法，根据各类实测砂岩储层的孔隙度和埋深数据计算各类砂岩储层对应的回归系数，构建回归方程。
在本实施方式中，可以通过地质资料分析的方法获得目标储层的砂岩类型，对应不同的回归方程，根据目标储层的埋深计算得到目标储层在指定埋深作用下的孔隙度。
在本实施方式中，通过根据不同砂岩类型计算不同砂岩类型对应的回归系数，目标储层可以根据自身储层的砂岩类型，计算指定埋深作用下目标储层的孔隙度进一步提高了指定埋深作用下目标储层孔隙度的计算精度。
在一个实施方式中，所述计算在指定构造应力作用下目标储层的孔隙度减孔量，包括：根据实测砂岩储层孔隙度和构造应力数据，通过回归算法计算回归系数；根据所述目标储层所受的构造应力和所述回归系数计算在指定构造应力作用下目标储层的孔隙度。
在本实施方式中，所述回归算法可以包括，根据所述实测砂岩储层孔隙度和构造应力数据以及构造应力模拟试验可以得到所述实测砂岩储层孔隙度的减孔量和构造应力之间满足指数关系。可以以构造应力为自变量，实测砂岩储层孔隙度的减孔量为因变量构建指数回归方程，计算回归系数，得到指数回归方程。
在本实施方式中，计算指定构造应力作用下目标储层的孔隙度减孔量可以包括，将目标储层的所受到的构造应力代入已经构建的指数回归方程中，得到目标储层在指定构造应力作用的埋深作用下的孔隙度减孔量。所述构造应力可以通过岩心样品进行声发射法实验测试获得。
在本实施方式中，所述指定构造应力作用下目标储层的孔隙度减孔量可以根据公式：
计算得出；其中，h、i表示回归系数，T表示所述指定构造应力。
在本实施方式中，通过实测砂岩储层孔隙度减孔量和构造应力数据，建立回归方程，并将目标储层的所受的指定构造应力代入回归方程计算目标储层在指定作用下的孔隙度减孔量的方式。利用统计和数理分析的方法得到孔隙度与构造应力之间的关系，提高所述指定构造应力作用下目标砂岩储层的孔隙度计算精度。
在一个实施方式中，根据不同砂岩类型计算不同砂岩类型对应的回归系数。
在本实施方式中，计算不同砂岩类型对应的回归系数可以包括，按照砂岩类型为实测砂岩储层进行分类。具体地，可以将储层为中砂岩、分选性中、刚性颗粒含量60％到70％的分为一类；将粗砂岩、分选性好，刚性颗粒含量为70％到80％的分为一类。通过回归方法，根据各类实测砂岩储层的孔隙度减孔量和构造应力计算各类砂岩储层对应的回归系数，构建回归方程。
在本实施方式中，可以通过地质资料分析的方法获得目标储层的砂岩类型，对应不同的回归方程，根据目标储层所受到的构造应力计算得到目标储层在指定埋深作用下的孔隙度减孔量。
在本实施方式中，通过根据不同砂岩类型计算不同砂岩类型对应的回归系数，目标储层可以根据自身储层的砂岩类型，计算指定构造应力作用下目标储层的孔隙度减孔量进一步提高了指定构造应力作用下目标储层孔隙度减孔量的计算精度。
在一个实施方式中，计算所述埋深作用下的减孔系数和构造挤压作用下的减孔系数，包括：根据实测砂岩储层的地史数据，生成所述实测砂岩储层的埋藏史线、孔隙史线、地温史线；根据所述埋藏史线、孔隙史线、地温史线得到在目标储层埋深范围内的实测砂岩储层的孔隙度减孔量；根据构造应力作用下实测砂岩储层的孔隙度减孔量和所述实测砂岩储层的孔隙度减孔量计算实测砂岩储层在埋深作用下的减孔量；将所述在构造应力作用下实测砂岩储层的孔隙度减孔量和所述实测砂岩储层在埋深作用下的减孔量分别除以所述实测砂岩储层总减孔量得到埋深作用下的减孔系数和构造应力作用下的减孔系数。
在本实施方式中，所述地史数据可以包括，实测砂岩储层所在区域反映砂岩储层地质属性随时间变化的数据。可以包括，输入地层厚度、压实率、剥蚀量、地质时间、现今地层温度、区域地温梯度、大地热流值、地层厚度、现今储层孔隙度、原始孔隙度、埋藏过程中地层厚度、构造应力、胶结物含量，溶蚀增控量等。
在本实施方式中，所述埋藏史线可以包括随着地质年代的推移，储层埋深的变化的曲线。埋藏史曲线可以通过盆地模拟软件Petromod获取，输入地层厚度、压实率、剥蚀量、地质时间可以生成；
在本实施方式中，所述地温史线可以包括，随着地质年代的推移，地层地温的变化的曲线。地层热史曲线可以通过盆地模拟软件Petromod获取，输入现今地层温度、区域地温梯度、大地热流值、地层厚度可以生成；
在本实施方式中，所述孔隙史线可以包括随着地质年代的推移，储层孔隙度的变化的曲线。储层孔隙演化史曲线可以通过盆模软件Petromod，输入现今储层孔隙度、原始孔隙度、埋藏过程地层厚度、构造应力及胶结物含量、溶蚀增孔量可以生成。
在本实施方式中，所述根据所述埋藏史线、孔隙史线、地温史线得到在目标储层埋深范围内的实测砂岩储层的孔隙度减孔量可以包括，根据目标储层的埋深范围，在埋藏史线中标定该埋深范围，由该埋深范围对应到孔隙史线，孔隙史线中表示的该埋深范围内孔隙度的差值可以是在目标储层埋深范围内的实测砂岩储层的孔隙度减孔量。
在本实施方式中，所述根据构造应力作用下实测砂岩储层的孔隙度减孔量和所述实测砂岩储层的孔隙度减孔量计算实测砂岩储层在埋深作用下的减孔量可以包括，根据地质资料得到目标储层的构造应力，计算构造应力作用下实测砂岩储层的减孔量，可以将所述在目标储层埋深范围内的实测砂岩储层的孔隙度减孔量减去所述构造应力作用下实测砂岩储层的减孔量得到实测砂岩储层在埋深作用下的减孔量。
在本说明书实施方式中，综合考虑了地质演化过程中埋藏压实和构造挤压对孔隙度改变程度的大小，可以更加精确地计算目标储层的孔隙度。另外，取构造挤压作用下孔隙度的变化数据，对于超深层致密砂岩储层孔隙度的计算，在数据处理上较为简单，求取埋藏压实和构造挤压的减孔系数也更为方便。
在一个实施方式中，所述预测所述目标储层的孔隙度，可以包括，将所述目标储层的原始孔隙度减去所述指定埋深作用下目标储层的孔隙度得到所述指定埋深作用下目标储层的孔隙度减孔量；将所述指定埋深作用下目标储层的孔隙度减孔量乘以所述埋深作用下减孔系数加上所述指定构造应力作用下目标储层的孔隙度减孔量乘以所述构造应力作用下的减孔系数得到所述预测的目标储层的孔隙度。
在本实施方式中，所述目标储层的孔隙度可以根据公式，
预测计算得出；其中，
表示目标储层孔隙度，
表示原始沉积体的砂岩孔隙度，
表示构造挤压作用下目标储层的减孔量，
构造挤压作用下目标储层的减孔量，j表示埋藏压实作用下的减孔量，k表示构造挤压作用下的减孔量，
表示埋藏压实作用下目标储层的孔隙度。
请参阅图2。本说明书实施方式还提供砂岩储层的孔隙度的预测装置，可以包括以下模块。
原始孔隙度计算模块；用于根据目标储层的原始沉积体的参数，计算目标储层的原始孔隙度；其中，所述原始沉积体用于表示所述目标储层刚刚沉积时的状态；所述原始孔隙度用于表示所述目标储层刚刚沉积时的孔隙度；埋深作用下孔隙度计算模块；用于根据目标储层的埋深，计算在指定埋深作用下目标储层的孔隙度；构造应力作用下孔隙度减孔量计算模块；用于根据目标储层的构造应力，计算在指定构造应力作用下目标储层的孔隙度减孔量；减孔系数计算模块；根据实测砂岩储层的孔隙度减孔量计算埋深作用下的减孔系数和构造应力作用下的减孔系数；孔隙度预测模块；用于根据所述目标储层的原始孔隙度、埋深作用下目标储层的孔隙度、构造应力作用下目标储层的孔隙度减孔量和埋深作用下的减孔系数以及构造应力作用下的减孔系数计算所述目标储层的孔隙度。
本说明书实施方式阐明的装置或模块等，具体可以由计算机芯片或实体实现，或者有具有某种功能的产品来实现。为了描述的方便，描述以上装置时以功能分为各种模块分别描述。当然，在申请本说明书时可以把各模块的功能在同一或多个软件和/或硬件中实现，也可以将实现同一功能的模块由多个子模块的组合实现等。以上所描述的装置实施例仅仅是示意性的，例如，所述模块的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个模块或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。
本领域技术人员也知道，除了以纯计算机可读程序代码实现控制器以外，完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件，而对其内部包括的用于实现各种功能的装置视为既可以是实现方法的软件模块，又可以是硬件部件内的结构。
本说明书实施方式可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构、类等等。也可以在分布式计算环境中实施本说明书，在这些分布式计算环境中，由于通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
本说明书实施方式还提供一种计算机存储介质，所述计算机存储介质存储有计算机程序指令，在所述计算机程序指令被执行时实现，根据目标储层的原始沉积体的参数，计算目标储层的原始孔隙度；其中，所述原始沉积体用于表示在沉积末期，埋藏初期所形成的沉积物；所述原始孔隙度用于表示所述目标储层在沉积末期，埋藏初期时的孔隙度；根据目标储层的埋深，计算在指定埋深作用下目标储层的孔隙度；根据目标储层的构造应力，计算在指定构造应力作用下目标储层的孔隙度减孔量；根据实测砂岩储层的孔隙度减孔量计算埋深作用下的减孔系数和构造应力作用下的减孔系数；根据所述目标储层的原始孔隙度、埋深作用下目标储层的孔隙度、构造应力作用下目标储层的孔隙度减孔量和埋深作用下的减孔系数以及构造应力作用下的减孔系数预测所述目标储层的孔隙度。
本说明书实施方式提供的砂岩储层的孔隙度的预测方法、装置及计算机存储介质，通过计算得到目标储层原始孔隙度、埋深作用下目标储层的孔隙度、构造应力作用下目标储层的孔隙度以及埋深作用下的减孔系数和构造应力作用下的减孔系数，综合计算目标储层的孔隙度。考虑了砂岩原始沉积体对储层原生孔隙及其演化的影响，还考虑了后期埋藏压实以及构造挤压作用对储层孔隙度的影响，并计算埋藏压实与构造挤压作用下的减孔系数，最终应用预测孔隙度，解决了现有技术中砂岩储层孔隙度预测过程中因为没有考虑到孔隙度发育因素多样性所带来的不足。
通过本说明书实施方式的描述可知，本领域的技术人员可以清楚地了解到本说明书实施方式可借助软件加必须的通用硬件平台的方式来实现。基于这样的理解，本说明书的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在存储介质中，如RAM/ROM、磁碟、光盘，包括若干指令用以使得一台计算机设备(可以是个人计算机，移动终端，服务器，或者网络设备等)执行本说明书实施方式或者实施方式的某些部分所述的方法。
本说明书的各个实施方式采用递进的方式描述，各个实施例之间相同或相似的部分互相参见即可，没各实施例重点说明的都是与其他实施例的不同之处。本说明书可用于众多通用或专用的计算机系统环境或配置中。例如：个人计算机、服务器计算机、手持设备或者便携式设备、平板型设备、多处理器系统、基于微处理器的系统、置顶盒、可编程的电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统和设备的分布式计算环境等等。
虽然本说明书提供了如实施方式或流程图所述的方法操作步骤，但基于常规或者无创造性的手段可以包括更多或更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序的一种方式，不代表唯一的执行顺序。在实际中的装置或客户端产品执行时，可以按照实施方式或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境，甚至为分布式数据处理环境)。术语“包括”、“包含”或者任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、产品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、产品或者设备所固有的要素。在没有更多限制的情况下，并不排除在包括要素的过程、方法、产品或者设备中还存在另外的相同或等同要素。
虽然通过说明书实施方式描绘了本说明书，本领域普通技术人员知道，本说明书有许多变形和变化而不脱离本说明书的精神，希望所附的权利要求包括这些变形和变化而不脱离本说明书的精神。
Claims (14)
1.一种砂岩储层的孔隙度的预测方法，其特征在于，包括：
根据目标储层的原始沉积体的参数，计算目标储层的原始孔隙度；其中，所述原始沉积体用于表示在沉积末期，埋藏初期所形成的沉积物；所述原始孔隙度用于表示所述目标储层在沉积末期，埋藏初期时的孔隙度；
根据目标储层的埋深，计算在指定埋深作用下目标储层的孔隙度；
根据目标储层的构造应力，计算在指定构造应力作用下目标储层的孔隙度减孔量；
根据实测砂岩储层的孔隙度减孔量计算埋深作用下的减孔系数和构造应力作用下的减孔系数，其中，包括：根据实测砂岩储层的地史数据，生成所述实测砂岩储层的埋藏史线、孔隙史线、地温史线；根据所述埋藏史线、孔隙史线、地温史线得到在目标储层埋深范围内的实测砂岩储层的孔隙度减孔量；根据构造应力作用下实测砂岩储层的孔隙度减孔量和所述实测砂岩储层的孔隙度减孔量计算实测砂岩储层在埋深作用下的减孔量；将在所述构造应力作用下实测砂岩储层的孔隙度减孔量和所述实测砂岩储层在埋深作用下的减孔量分别除以所述实测砂岩储层总减孔量得到埋深作用下的减孔系数和构造应力作用下的减孔系数；
根据所述目标储层的原始孔隙度、埋深作用下目标储层的孔隙度、构造应力作用下目标储层的孔隙度减孔量和埋深作用下的减孔系数以及构造应力作用下的减孔系数预测所述目标储层的孔隙度。
2.根据权利要求1所述的方法，其特征在于，所述目标储层的原始沉积体的参数至少包括：所述原始沉积体的平均粒径、分选系数、泥质含量以及刚性颗粒含量数据。
3.根据权利要求2所述的方法，其特征在于，计算所述原始孔隙度包括：
根据实测砂岩储层原始沉积体的孔隙度、平均粒径、分选系数、泥质含量与刚性颗粒含量通过回归算法得到所述平均粒径、分选系数、泥质含量与刚性颗粒含量对应的回归系数；
根据与所述目标储层关联的指定地层的平均粒径、分选系数、泥质含量与刚性颗粒含量以及所述回归系数计算所述原始孔隙度。
4.根据权利要求3所述的方法，其特征在于，所述原始孔隙度根据公式：φ0＝a×So+b×D+c×Vsh+d×VQ+e计算得出；其中，φ0表示目标储层的原始孔隙度，So表示分选系数，D表示平均粒径、Vsh表示泥质含量、VQ表示刚性颗粒含量，a、b、c、d、e分别表示回归系数。
5.根据权利要求1所述的方法，其特征在于，计算所述指定埋深作用下目标储层的孔隙度，包括：
根据实测砂岩储层孔隙度和埋深数据，通过回归算法计算回归系数；
根据所述目标储层的埋深和所述回归系数计算在指定埋深作用下目标储层的孔隙度。
6.根据权利要求5所述的方法，其特征在于，所述方法还包括，根据不同的砂岩类型计算不同砂岩类型对应的回归系数。
7.根据权利要求5所述方法，其特征在于，所述指定埋深作用下目标储层的孔隙度，根据公式：
计算得出；其中f、g表示回归系数，H表示所述指定埋深。
8.根据权利要求1所述的方法，其特征在于，所述计算在指定构造应力作用下目标储层的孔隙度减孔量，包括：
根据实测砂岩储层孔隙度和构造应力数据，通过回归算法计算回归系数；
根据所述目标储层所受的构造应力和所述回归系数计算在指定构造应力作用下目标储层的孔隙度。
9.根据权利要求8所述的方法，其特征在于，所述方法还包括，根据不同砂岩类型计算不同砂岩类型对应的回归系数。
10.根据权利要求8所述的方法，其特征在于，所述计算在指定构造应力作用下目标储层的孔隙度减孔量，根据公式
计算得出；其中，h、i表示回归系数，T表示所述指定构造应力。
11.根据权利要求1所述方法，其特征在于，所述预测所述目标储层的孔隙度，包括：
将所述目标储层的原始孔隙度减去所述指定埋深作用下目标储层的孔隙度得到所述指定埋深作用下目标储层的孔隙度减孔量；
将所述指定埋深作用下目标储层的孔隙度减孔量乘以所述埋深作用下减孔系数加上所述指定构造应力作用下目标储层的孔隙度减孔量乘以所述构造应力作用下的减孔系数得到所述预测的目标储层的孔隙度。
12.根据权利要求1所述的方法，其特征在于，所述目标储层的孔隙度，根据公式：
计算得出；其中，
表示目标储层孔隙度，
表示目标储层的原始孔隙度，
表示埋深作用下目标储层的孔隙度减孔量，
表示构造应力作用下目标储层的孔隙度减孔量，j表示埋深作用下的减孔系数，k表示构造应力作用下的减孔系数，
表示埋深作用下目标储层的孔隙度。
13.一种砂岩储层的孔隙度的预测装置，其特征在于，包括：
原始孔隙度计算模块；用于根据目标储层的原始沉积体的参数，计算目标储层的原始孔隙度；其中，所述原始沉积体用于表示在沉积末期，埋藏初期所形成的沉积物；所述原始孔隙度用于表示所述目标储层在沉积末期，埋藏初期时的孔隙度；
埋深作用下孔隙度计算模块；用于根据目标储层的埋深，计算在指定埋深作用下目标储层的孔隙度；
构造应力作用下孔隙度减孔量计算模块；用于根据目标储层的构造应力，计算在指定构造应力作用下目标储层的孔隙度减孔量；
减孔系数计算模块；根据实测砂岩储层的孔隙度减孔量计算埋深作用下的减孔系数和构造应力作用下的减孔系数，其中，包括：根据实测砂岩储层的地史数据，生成所述实测砂岩储层的埋藏史线、孔隙史线、地温史线；根据所述埋藏史线、孔隙史线、地温史线得到在目标储层埋深范围内的实测砂岩储层的孔隙度减孔量；根据构造应力作用下实测砂岩储层的孔隙度减孔量和所述实测砂岩储层的孔隙度减孔量计算实测砂岩储层在埋深作用下的减孔量；将在所述构造应力作用下实测砂岩储层的孔隙度减孔量和所述实测砂岩储层在埋深作用下的减孔量分别除以所述实测砂岩储层总减孔量得到埋深作用下的减孔系数和构造应力作用下的减孔系数；
孔隙度预测模块；用于根据所述目标储层的原始孔隙度、埋深作用下目标储层的孔隙度、构造应力作用下目标储层的孔隙度减孔量和埋深作用下的减孔系数以及构造应力作用下的减孔系数预测所述目标储层的孔隙度。
14.一种计算机存储介质，其特征在于，所述计算机存储介质存储有计算机程序指令，在所述计算机程序指令被执行时实现：
根据目标储层的原始沉积体的参数，计算目标储层的原始孔隙度；其中，所述原始沉积体用于表示在沉积末期，埋藏初期所形成的沉积物；所述原始孔隙度用于表示所述目标储层在沉积末期，埋藏初期时的孔隙度；
根据目标储层的埋深，计算在指定埋深作用下目标储层的孔隙度；
根据目标储层的构造应力，计算在指定构造应力作用下目标储层的孔隙度减孔量；
根据实测砂岩储层的孔隙度减孔量计算埋深作用下的减孔系数和构造应力作用下的减孔系数，其中，包括：根据实测砂岩储层的地史数据，生成所述实测砂岩储层的埋藏史线、孔隙史线、地温史线；根据所述埋藏史线、孔隙史线、地温史线得到在目标储层埋深范围内的实测砂岩储层的孔隙度减孔量；根据构造应力作用下实测砂岩储层的孔隙度减孔量和所述实测砂岩储层的孔隙度减孔量计算实测砂岩储层在埋深作用下的减孔量；将在所述构造应力作用下实测砂岩储层的孔隙度减孔量和所述实测砂岩储层在埋深作用下的减孔量分别除以所述实测砂岩储层总减孔量得到埋深作用下的减孔系数和构造应力作用下的减孔系数；
根据所述目标储层的原始孔隙度、埋深作用下目标储层的孔隙度、构造应力作用下目标储层的孔隙度减孔量和埋深作用下的减孔系数以及构造应力作用下的减孔系数预测所述目标储层的孔隙度。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810011059.9A CN108227033B (zh) | 2018-01-05 | 2018-01-05 | 砂岩储层的孔隙度的预测方法、装置及计算机存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810011059.9A CN108227033B (zh) | 2018-01-05 | 2018-01-05 | 砂岩储层的孔隙度的预测方法、装置及计算机存储介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108227033A CN108227033A (zh) | 2018-06-29 |
| CN108227033B true CN108227033B (zh) | 2020-09-08 |
Family
ID=62642941
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810011059.9A Active CN108227033B (zh) | 2018-01-05 | 2018-01-05 | 砂岩储层的孔隙度的预测方法、装置及计算机存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108227033B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109444957B (zh) * | 2018-10-10 | 2019-09-27 | 广州海洋地质调查局 | 基于决策树评价的水合物储层地震预测方法及处理终端 |
| CN113138414A (zh) * | 2020-07-24 | 2021-07-20 | 中国石油天然气股份有限公司 | 深层砂岩储层孔隙度的预测方法及装置 |
| CN113608260B (zh) * | 2021-07-20 | 2022-04-22 | 中国石油大学胜利学院 | 砂岩储层胶结前的多期构造压实减孔程度定量评价方法 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7933757B2 (en) * | 2004-03-31 | 2011-04-26 | Exxonmobil Upstream Research Co. | Method for simulating and estimating sandstone properties |
| CN102748016B (zh) * | 2012-07-10 | 2014-10-08 | 中国石油大学(华东) | 地质历史时期砂岩储层孔隙度演化恢复方法 |
| CN105334149B (zh) * | 2015-11-24 | 2019-06-14 | 中国石油大学(北京) | 一种致密储层微观孔隙结构评价与储层分类方法 |
-
2018
- 2018-01-05 CN CN201810011059.9A patent/CN108227033B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN108227033A (zh) | 2018-06-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107657365B (zh) | 地质资源开采价值评价方法及装置 | |
| CN103645516B (zh) | 基于岩石物理相控油气作用确定油气产能的方法 | |
| Freeze | A stochastic‐conceptual analysis of one‐dimensional groundwater flow in nonuniform homogeneous media | |
| Haldorsen | Simulator parameter assignment and the problem of scale in reservoir engineering | |
| CN109653725B (zh) | 基于沉积微相和岩石相的混积储层水淹程度测井解释方法 | |
| US10048403B2 (en) | Method and system for generation of upscaled mechanical stratigraphy from petrophysical measurements | |
| CN110318745B (zh) | 一种沉积微相约束下粒径岩性测井评价方法 | |
| CN104832166A (zh) | 一种页岩气水平井初期产能预测方法 | |
| CN107366534B (zh) | 粗化渗透率的确定方法和装置 | |
| CN108227033B (zh) | 砂岩储层的孔隙度的预测方法、装置及计算机存储介质 | |
| US9892321B2 (en) | Using maximal inscribed spheres for image-based rock property estimation | |
| CN109115987B (zh) | 一种基于岩石物理模型的流体因子的评价方法及装置 | |
| CN104239743B (zh) | 确定岩性油气藏成藏概率的方法和装置 | |
| Liu et al. | Quantitative multiparameter prediction of fractured tight sandstone reservoirs: a case study of the Yanchang Formation of the Ordos Basin, Central China | |
| Li et al. | Analysis of the key factors affecting the productivity of coalbed methane wells: A case study of a high-rank coal reservoir in the central and southern Qinshui Basin, China | |
| Wang et al. | Quantification and prediction of pore structures in tight oil reservoirs based on multifractal dimensions from integrated pressure-and rate-controlled porosimetry for the upper Triassic Yanchang formation, Ordos Basin, China | |
| CN106285652B (zh) | 确定页岩游离气体饱和度的方法 | |
| Wang et al. | Estimation of probability distribution of shear strength of slip zone soils in Middle Jurassic red beds in Wanzhou of China | |
| Aljuboori et al. | The impact of diagenesis precipitation on fracture permeability in naturally fractured carbonate reservoirs | |
| CN111323814A (zh) | 基于岩石物理模板定量确定砂岩储层含水饱和度的方法 | |
| Li et al. | Cyclic fractional plastic model for granular soils | |
| CN114460639B (zh) | 页岩油储层渗透率的预测方法及装置 | |
| CN110632665B (zh) | 一种砂体构型量化分析方法 | |
| CN116930023A (zh) | 致密砂岩相控分类孔隙度测井精细解释方法和装置 | |
| Rahimi et al. | Investigation of linear and non-linear estimation methods in highly-skewed gold distribution |
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 |