CN111239811B - 一种判断储层流体趋势与压实趋势的方法及装置 - Google Patents
一种判断储层流体趋势与压实趋势的方法及装置 Download PDFInfo
- Publication number
- CN111239811B CN111239811B CN201811469482.XA CN201811469482A CN111239811B CN 111239811 B CN111239811 B CN 111239811B CN 201811469482 A CN201811469482 A CN 201811469482A CN 111239811 B CN111239811 B CN 111239811B
- Authority
- CN
- China
- Prior art keywords
- attribute
- orthogonal
- pseudo
- factor
- wave velocity
- 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
- 239000012530 fluid Substances 0.000 title claims abstract description 73
- 238000005056 compaction Methods 0.000 title claims abstract description 72
- 238000000034 method Methods 0.000 title claims abstract description 60
- 239000011435 rock Substances 0.000 claims abstract description 22
- 238000004458 analytical method Methods 0.000 claims abstract description 12
- 230000015654 memory Effects 0.000 claims description 31
- 238000004590 computer program Methods 0.000 claims description 19
- 238000003860 storage Methods 0.000 claims description 14
- 238000004364 calculation method Methods 0.000 claims description 10
- 239000011159 matrix material Substances 0.000 claims description 10
- 239000000126 substance Substances 0.000 claims description 5
- 230000002238 attenuated effect Effects 0.000 claims 1
- 230000006870 function Effects 0.000 description 17
- 238000010586 diagram Methods 0.000 description 14
- 238000009826 distribution Methods 0.000 description 12
- 230000008569 process Effects 0.000 description 11
- 230000008859 change Effects 0.000 description 9
- 230000006872 improvement Effects 0.000 description 9
- 238000012545 processing Methods 0.000 description 8
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 8
- 238000011161 development Methods 0.000 description 7
- 241001442234 Cosa Species 0.000 description 4
- 238000004891 communication Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 230000003287 optical effect Effects 0.000 description 4
- 244000089409 Erythrina poeppigiana Species 0.000 description 3
- 235000009776 Rathbunia alamosensis Nutrition 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000005553 drilling Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000001514 detection method Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000005055 memory storage Effects 0.000 description 2
- 230000000750 progressive effect Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 108010001267 Protein Subunits Proteins 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 239000007795 chemical reaction product Substances 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000005562 fading Methods 0.000 description 1
- 229910021389 graphene Inorganic materials 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 229920001296 polysiloxane Polymers 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 239000000047 product Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 229930195734 saturated hydrocarbon Natural products 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000003068 static effect Effects 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
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
本说明书公开了一种判断储层流体趋势与压实趋势的方法及装置，其中，判断储层流体趋势与压实趋势的方法包括：获取目标区域的叠前地震道集；根据目标区域的叠前地震道集获得纵波速度、横波速度和密度数据体；利用所述纵波速度、所述横波速度和所述密度数据体获得伪衰减属性体；根据所述伪衰减属性体得到基于旋转角度为参数的正交属性因子表达式，基于岩石物理分析获得旋转角度值，将所述旋转角度值代入所述正交属性因子表达式获得对应正交属性因子；基于正交属性因子判断储层流体趋势和压实趋势。
Description
技术领域
本说明书涉及石油地球物理勘探技术领域，特别涉及一种判断储层流体趋势与压实趋势的方法及装置。
背景技术
在石油勘探开发过程中，储层中流体分布规律以及储层岩石压实趋势是两项非常重要的指标，其决定了储层的开发潜力。储层压实趋势反映了储层孔隙空间的分布规律，压实越严重的区域储层储集空间越小，导致单位体积内可采资源量越少，流体趋势反映了储层中不同类型流体的分布特征，单位储层空间内含水饱和度越高，储层含烃量就越少，储层的经济性就越差。预测储层的压实趋势与流体趋势对于提高石油勘探的成功率具有重大意义，常规预测方法通常通过振幅属性分析或者叠后地震波阻抗反演建立属性与储层物性的关系，用以预测储层的压实趋势以及流体分布规律。常规方法通常是一个属性解决一个问题，每种属性代表的物理含义以及提取的方法存在较大的差异，导致在预测过程中出现相互矛盾的地方，增加了储层预测的多解性。
发明内容
为了实现对储层压实趋势与流体趋势的准确预测，本说明书提供一种判断储层流体趋势与压实趋势的方法及装置。
为实现上述目的，本说明书实施方式提供一种判断储层流体趋势与压实趋势的方法，包括：
获取目标区域的叠前地震道集；
根据目标区域的叠前地震道集获得纵波速度、横波速度和密度数据体；
利用所述纵波速度、所述横波速度和所述密度数据体获得伪衰减属性体；
根据所述伪衰减属性体得到基于旋转角度为参数的正交属性因子表达式，基于岩石物理分析获得旋转角度值，将所述旋转角度值代入所述正交属性因子表达式获得对应正交属性因子；
基于正交属性因子判断储层流体趋势和压实趋势；
其中，所述伪衰减属性体包括纵波对应的伪衰减属性和横波对应的伪衰减属性；
伪衰减属性体由以下公式计算得到：
其中，FQP、FQS分别为纵横波对应的伪衰减属性，Vp为纵波速度，Vs为横波速度，ρ为密度。
优选地，根据目标区域的叠前地震道集获得纵波速度、横波速度和密度数据体的步骤包括：
对所述目标区域的叠前地震道集进行叠加处理；
根据叠加的叠前地震道集反演得到纵波速度、横波速度和密度数据体。
优选地，所述正交属性因子包括压实趋势因子和流体趋势因子。
优选地，根据所述伪衰减属性体得到基于旋转角度为参数的正交属性因子表达式的步骤包括：
将所述伪衰减属性体进行正交旋转运算，得到基于旋转角度为参数的正交属性因子表达式。
优选地，所述正交属性因子表达式为：
其中，FQP为纵波对应的伪衰减属性；FQS为横波对应的伪衰减属性；C为压实趋势因子，C＝FQS sina-FQP cosa；F为流体趋势因子，F＝FQS cosa+FQP sina；
为正交旋转矩阵；a为旋转角；
根据岩石物理参数确定旋转角度值，旋转角度计算公式如下所示：
其中，
分别为干燥岩石骨架的纵波速度和横波速度。
为实现上述目的，本说明书实施方式提供一种判断储层流体趋势与压实趋势的装置，包括：
叠前地震道集获取单元，用于获取目标区域的叠前地震道集；
反演单元，用于根据目标区域的叠前地震道集获得纵波速度、横波速度和密度数据体；
伪衰减属性体获取单元，用于利用所述纵波速度、所述横波速度和所述密度数据体获得伪衰减属性体；
正交属性因子获取单元，用于根据所述伪衰减属性体得到基于旋转角度为参数的正交属性因子表达式，基于岩石物理分析获得旋转角度值，将所述旋转角度值代入所述正交属性因子表达式获得对应正交属性因子；
判断单元，用于基于正交属性因子判断储层流体趋势和压实趋势；
其中，所述伪衰减属性体获取单元获取的伪衰减属性体包括纵波对应的伪衰减属性和横波对应的伪衰减属性；
伪衰减属性体由以下公式计算得到：
其中，FQP、FQS分别为纵横波对应的伪衰减属性，Vp为纵波速度，Vs为横波速度，ρ为密度。
优选地，所述反演单元包括：
叠加模块，用于对所述目标区域的叠前地震道集进行叠加处理；
反演计算模块，用于根据叠加的叠前地震道集反演得到纵波速度、横波速度和密度数据体。
优选地，所述正交属性因子获取单元获取的正交属性因子包括压实趋势因子和流体趋势因子。
优选地，所述正交属性因子获取单元将所述伪衰减属性体进行正交旋转运算，得到基于旋转角度为参数的正交属性因子表达式。
优选地，所述正交属性因子获取单元获取的正交属性因子表达式为：
其中，FQP为纵波对应的伪衰减属性；FQS为横波对应的伪衰减属性；C为压实趋势因子，C＝FQS sina-FQP cosa；F为流体趋势因子，F＝FQS cosa+FQP sina；
为正交旋转矩阵；a为旋转角；
根据岩石物理参数确定旋转角度值，旋转角度计算公式如下所示：
其中，
分别为干燥岩石骨架的纵波速度和横波速度。
为实现上述目的，本说明书实施方式提供一种电子设备，包括存储器、处理器以及存储在所述存储器上并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述所述的判断储层流体趋势与压实趋势的方法。
为实现上述目的，本说明书实施方式提供一种可读存储介质，其上存储有计算机程序，所述计算机程序被执行时实现上述所述的判断储层流体趋势与压实趋势的方法的步骤。
由上可见，与现有技术相比较，本技术方案利用正交属性因子有效预测储层的流体趋势与压实趋势，克服常规储层预测的多解性，提高勘探开发的精度。
附图说明
为了更清楚地说明本说明书实施方式或现有技术中的技术方案，下面将对实施方式或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本说明书中记载的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本说明书实施方式提供一种判断储层流体趋势与压实趋势的方法流程图；
图2为伪衰减属性在岩石物理模板中的分布规律示意图；
图3为压实趋势因子和流体趋势因子在岩石物理模板中的分布规律示意图；
图4为根据压实趋势因子预测的储层压实趋势平面图；
图5为根据流体趋势因子预测的储层流体分布规律图；
图6为本说明书实施方式提供一种判断储层流体趋势与压实趋势的装置功能框图；
图7为本说明书实施例提出的一种电子设备示意图。
具体实施方式
下面将结合附图，对本公开实施例中的技术方案进行清楚、完整地描述参考在附图中示出并在以下描述中详述的非限制性示例实施例，更加全面地说明本公开的示例实施例和它们的多种特征及有利细节。应注意的是，图中示出的特征不是必须按照比例绘制。本公开省略了已知材料、组件和工艺技术的描述，从而不使本公开的示例实施例模糊。所给出的示例仅旨在有利于理解本公开示例实施例的实施，以及进一步使本领域技术人员能够实施示例实施例。因而，这些示例不应被理解为对本公开的实施例的范围的限制。
除非另外特别定义，本公开使用的技术术语或者科学术语应当为本公开所属领域内具有一般技能的人士所理解的通常意义。本公开中使用的“第一”、“第二”以及类似的词语并不表示任何顺序、数量或者重要性，而只是用来区分不同的组成部分。此外，在本公开各个实施例中，相同或类似的参考标号表示相同或类似的构件。
在实际石油勘探开发过程，需要一种基于一套岩石物理理论同时解决储层压实趋势以及流体趋势预测的技术方案，保证在理论基础一致的前提下，提高预测结果的可靠性与可对比性。为了实现对储层压实趋势与流体趋势的准确预测，本说明书提出了一种判断储层流体趋势与压实趋势的技术方案，旨在通过叠前反演数据体构建基于伪衰减特性的属性因子，并且通过对属性进行正交旋转以及岩石物理分析，推导出新的表征储层压实与流体特征的正交属性因子，利用正交属性因子有效预测储层的流体趋势与压实趋势，克服常规储层预测的多解性，提高勘探开发的精度。
基于上述描述，如图1所示，为本说明书实施方式提供一种判断储层流体趋势与压实趋势的方法流程图。包括：
步骤101)：获取目标区域的叠前地震道集。
步骤102)：根据目标区域的叠前地震道集获得纵波速度、横波速度和密度数据体。
在本步骤中，利用叠前道集数据运算得到纵波、横波速度以及密度数据体，公式如下所示：
2Rpp(θ)＝(1+tan2θ)Δln(Vp)-8Ksin2θΔln(Vs)+(1-4Ksin2θ)Δln(ρ)
其中，Rpp(θ)为叠前地震道集数据，θ为入射角，Vp为纵波速度，Vs为横波速度，ρ为密度，K为常数。
在本实施例中，选取三个角度的道集数据进行基于L2范数的叠前反演即可得到纵波速度、横波速度以及密度数据体。
步骤103)：利用所述纵波速度、所述横波速度和所述密度数据体获得伪衰减属性体。
在本步骤中，伪衰减属性体由以下公式计算得到：
其中，FQP、FQS分别为纵横波对应的伪衰减属性。Vp为纵波速度，Vs为横波速度，ρ为密度。
图2展示了伪衰减属性在岩石物理模板中的分布规律。其中，x轴为横波伪衰减属性FQS，y轴为纵波伪衰减属性FQP。在图2中，虚线箭头指示了储层压实趋势的变化方向，实线箭头指示了储层流体趋势的变化方向，可以看出伪衰减属性的变化同时受到储层压实趋势与流体趋势的影响，两种趋势变化耦合在一起同时对纵横波伪衰减属性产生影响，导致无法单独使用一种伪衰减属性来判断储层的压实趋势或者流体趋势，增大了储层预测的多解性。
步骤104)：根据所述伪衰减属性体得到基于旋转角度为参数的正交属性因子表达式，基于岩石物理分析获得旋转角度值，将所述旋转角度值代入所述正交属性因子表达式获得对应正交属性因子。
在本步骤中，将伪衰减属性体进行正交旋转运算，得到正交属性因子，包括：
定义正交旋转矩阵：
其中，a旋转角，正交旋转矩阵满足AAT＝E,AT为A矩阵的转置，E为二阶单位阵。
正交属性因子包括压实趋势因子和流体趋势因子，通过伪衰减属性旋转得到，公式如下所示：
其中，C＝FQS sina-FQP cosa为压实趋势因子，F＝FQS cosa+FQP sina为流体趋势因子。
基于岩石物理分析获得旋转角度值，并得到对应角度的正交属性因子，包括：
根据岩石物理参数确定旋转角度值，旋转角度计算公式如下所示：
其中，
分别为干燥岩石骨架的纵波速度和横波速度，两个参数通过岩石样本实验室测算得到，也可以选择研究工区默认参考值。在本实施例中选择工区默认参考值进行试算，默认参考值参照了工区内多口井的岩石物理数据，取平均后得到一个能够反映工区整体趋势的值，其中：
旋转角度a≈58°。
根据计算结果a≈58°，则正交旋转矩阵为：
正交属性因子包括压实趋势因子和流体趋势因子，其通过伪衰减属性旋转得到，公式如下所示：
其中，C＝0.85FQS-0.53FQP为压实趋势因子，F＝0.53FQS+0.85FQP为流体趋势因子，二者在岩石物理模板中的分布如图3所示，图中x轴为压实趋势因子C，y轴为流体趋势因子F，可以看到两种属性因子均为单调单向变化，即压实趋势因子C完全反映了储层的压实趋势变化，随着C值增大，储层孔隙度逐渐增大，且C不受到流体的影响。同样的，流体趋势因子F完全反映了储层的流体趋势变化，随着F增大，储层含水饱和度逐渐较小，且F不受到孔隙度的影响。两种正交趋势因子分别无干扰的反映了储层的流体趋势与压实趋势，将孔隙度与含水饱和度对储层的影响分别进行表征，大大降低储层预测的多解性，提高了流体检测与孔隙度预测的精度。
步骤105)：基于正交属性因子判断储层流体趋势和压实趋势。
储层流体趋势单向反映了储层含水饱和度的变化趋势，压实趋势单向反映了储层孔隙度的变化趋势，依据两种属性的变化规律判断选取压实程度小、含水饱和度低的区域即为有效储层发育区。
在本实施例中，随着压实趋势因子C值增大，储层孔隙度逐渐增大，随着流体趋势因子F值增大，储层含水饱和度逐渐较小。依据两种属性的变化规律判断选取压实程度小、含水饱和度低的区域即为有效储层发育区。
为了验证储层流体趋势与压实趋势预测结果的准确性，将检测结果与实际工区的钻井和测井数据进行对比分析。如图4所示，为根据压实趋势因子预测的储层压实趋势平面图。如图5所示，为根据流体趋势因子预测的储层流体分布规律图。图中虚线区域有一口井，该井为工业气井，日产气6000方，属于高产气井。预测结果显示，在图4中，在虚线高产气井区域内压实趋势因子C显示为高值异常，表明该区域为高孔隙度储层，储层物性非常好，储集空间大，这与钻井数据吻合。在图5中，在虚线高产气井区域内流体趋势因子F显示为高值异常，表明该区域为低含水高饱和烃储层，这与高产气井的结果相吻合。综合分析两种属性因子，压实趋势因子C指示了高孔隙度储层的分布范围，流体趋势因子F指示了低含水储层的分布范围，二者结合则圈定了高孔隙度、低含水饱和度优质储层的分布，排除了高含水层或者致密气层的影响，降低了储层预测的风险，大大提高了钻井成功率。
如图6所示，为本说明书实施方式提供一种判断储层流体趋势与压实趋势的装置功能框图。包括：
叠前地震道集获取单元601，用于获取目标区域的叠前地震道集；
反演单元602，用于根据目标区域的叠前地震道集获得纵波速度、横波速度和密度数据体；
伪衰减属性体获取单元603，用于利用所述纵波速度、所述横波速度和所述密度数据体获得伪衰减属性体；
正交属性因子获取单元604，用于根据所述伪衰减属性体得到基于旋转角度为参数的正交属性因子表达式，基于岩石物理分析获得旋转角度值，将所述旋转角度值代入所述正交属性因子表达式获得对应正交属性因子；
判断单元605，用于基于正交属性因子判断储层流体趋势和压实趋势。
优选地，所述反演单元包括：
叠加模块，用于对所述目标区域的叠前地震道集进行叠加处理；
反演计算模块，用于根据叠加的叠前地震道集反演得到纵波速度、横波速度和密度数据体。
优选地，所述伪衰减属性体获取单元获取的伪衰减属性体包括纵波对应的伪衰减属性和横波对应的伪衰减属性。
优选地，所述正交属性因子获取单元获取的正交属性因子包括压实趋势因子和流体趋势因子。
优选地，所述正交属性因子获取单元将所述伪衰减属性体进行正交旋转运算，得到基于旋转角度为参数的正交属性因子表达式。
优选地，所述正交属性因子获取单元获取的正交属性因子表达式为：
其中，FQP为纵波对应的伪衰减属性；FQS为横波对应的伪衰减属性；C为压实趋势因子，C＝FQS sina-FQP cosa；F为流体趋势因子，F＝FQS cosa+FQP sina；
为正交旋转矩阵；a为旋转角。
如图7所示，为本说明书实施例提出的一种电子设备示意图。包括存储器、处理器以及存储在所述存储器上并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述图1所述的判断储层流体趋势与压实趋势的方法。
本说明书实施方式提供的判断储层流体趋势与压实趋势的方法，其存储器和处理器实现的具体功能，可以与本说明书中的前述实施方式相对照解释，并能够达到前述实施方式的技术效果，这里便不再赘述。
在本实施方式中，所述存储器可以包括用于存储信息的物理装置，通常是将信息数字化后再以利用电、磁或者光学等方法的媒体加以存储。本实施方式所述的存储器又可以包括：利用电能方式存储信息的装置，如RAM、ROM等；利用磁能方式存储信息的装置，如硬盘、软盘、磁带、磁芯存储器、磁泡存储器、U盘；利用光学方式存储信息的装置，如CD或DVD。当然，还有其他方式的存储器，例如量子存储器、石墨烯存储器等等。
在本实施方式中，所述处理器可以按任何适当的方式实现。例如，所述处理器可以采取例如微处理器或处理器以及存储可由该(微)处理器执行的计算机可读程序代码(例如软件或固件)的计算机可读介质、逻辑门、开关、专用集成电路(Application SpecificIntegrated Circuit，ASIC)、可编程逻辑控制器和嵌入微控制器的形式等等。
在本实施例中，本说明书实施例还提供一种可读存储介质，其上存储有计算机程序，所述计算机程序被执行时实现上述所述的判断储层流体趋势与压实趋势的方法的步骤。
与现有技术相比较，本技术方案旨在通过对伪衰减属性进行正交旋转得到能够分别反映储层压实趋势与储层流体趋势的正交属性因子，其中伪衰减属性通过叠前反演得到的纵横波速以及密度获得，旋转角度通过岩石物理分析得到，最终得到的正交属性因子能够克服常规储层预测方法的多解性，提高储层预测的精度。
本领域技术人员也知道，除了以纯计算机可读程序代码方式实现客户端和服务器以外，完全可以通过将方法步骤进行逻辑编程来使得客户端和服务器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种客户端和服务器可以被认为是一种硬件部件，而对其内包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至，可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
通过以上的实施方式的描述可知，本领域的技术人员可以清楚地了解到本说明书可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解，本说明书的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在存储介质中，如ROM/RAM、磁碟、光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本说明书各个实施方式或者实施方式的某些部分所述的方法。
本说明书中的各个实施方式均采用递进的方式描述，各个实施方式之间相同相似的部分互相参见即可，每个实施方式重点说明的都是与其他实施方式的不同之处。尤其，针对客户端和服务器的实施方式来说，均可以参照前述方法的实施方式的介绍对照解释。
本说明书可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本说明书，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
在20世纪90年代，对于一个技术的改进可以很明显地区分是硬件上的改进(例如，对二极管、晶体管、开关等电路结构的改进)还是软件上的改进(对于方法流程的改进)。然而，随着技术的发展，当今的很多方法流程的改进已经可以视为硬件电路结构的直接改进。设计人员几乎都通过将改进的方法流程编程到硬件电路中来得到相应的硬件电路结构。因此，不能说一个方法流程的改进就不能用硬件实体模块来实现。例如，可编程逻辑器件(Programmable Logic Device,PLD)(例如现场可编程门阵列(Field Programmable GateArray，FPGA))就是这样一种集成电路，其逻辑功能由用户对器件编程来确定。由设计人员自行编程来把一个数字系统“集成”在一片PLD上，而不需要请芯片制造厂商来设计和制作专用的集成电路芯片。而且，如今，取代手工地制作集成电路芯片，这种编程也多半改用“逻辑编译器(logic compiler)”软件来实现，它与程序开发撰写时所用的软件编译器相类似，而要编译之前的原始代码也得用特定的编程语言来撰写，此称之为硬件描述语言(Hardware Description Language，HDL)，而HDL也并非仅有一种，而是有许多种，如ABEL(Advanced Boolean Expression Language)、AHDL(Altera Hardware DescriptionLanguage)、Confluence、CUPL(Cornell University Programming Language)、HDCal、JHDL(Java Hardware Description Language)、Lava、Lola、MyHDL、PALASM、RHDL(RubyHardware Description Language)等，目前最普遍使用的是VHDL(Very-High-SpeedIntegrated Circuit Hardware Description Language)与Verilog。本领域技术人员也应该清楚，只需要将方法流程用上述几种硬件描述语言稍作逻辑编程并编程到集成电路中，就可以很容易得到实现该逻辑方法流程的硬件电路。
控制器可以按任何适当的方式实现，例如，控制器可以采取例如微处理器或处理器以及存储可由该(微)处理器执行的计算机可读程序代码(例如软件或固件)的计算机可读介质、逻辑门、开关、专用集成电路(Application Specific Integrated Circuit，ASIC)、可编程逻辑控制器和嵌入微控制器的形式，控制器的例子包括但不限于以下微控制器：ARC 625D、Atmel AT91SAM、Microchip PIC18F26K20以及Silicone Labs C8051F320，存储器控制器还可以被实现为存储器的控制逻辑的一部分。本领域技术人员也知道，除了以纯计算机可读程序代码方式实现控制器以外，完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件，而对其内包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至，可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
上述实施例阐明的系统、装置、模块或单元，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。一种典型的实现设备为计算机。具体的，计算机例如可以为个人计算机、膝上型计算机、车载人机交互设备、蜂窝电话、相机电话、智能电话、个人数字助理、媒体播放器、导航设备、电子邮件设备、游戏控制台、平板计算机、可穿戴设备或者这些设备中的任何设备的组合。
虽然本说明书实施例提供了如实施例或流程图所述的方法操作步骤，但基于常规或者无创造性的手段可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的装置或终端产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境，甚至为分布式数据处理环境)。术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、产品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、产品或者设备所固有的要素。在没有更多限制的情况下，并不排除在包括所述要素的过程、方法、产品或者设备中还存在另外的相同或等同要素。
为了描述的方便，描述以上装置时以功能分为各种模块分别描述。当然，在实施本说明书实施例时可以把各模块的功能在同一个或多个软件和/或硬件中实现，也可以将实现同一功能的模块由多个子模块或子单元的组合实现等。以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
本领域技术人员也知道，除了以纯计算机可读程序代码方式实现控制器以外，完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件，而对其内部包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至，可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中，计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。
内存可能包括计算机可读介质中的非永久性存储器，随机存取存储器(RAM)和/或非易失性内存等形式，如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括暂存电脑可读媒体(transitory media)，如调制的数据信号和载波。
本领域技术人员应明白，本说明书的实施例可提供为方法、系统或计算机程序产品。因此，本说明书实施例可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本说明书实施例可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本说明书实施例可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本说明书实施例，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本说明书实施例的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
以上所述仅为本说明书实施例的实施例而已，并不用于限制本说明书实施例。对于本领域技术人员来说，本说明书实施例可以有各种更改和变化。凡在本说明书实施例的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本说明书实施例的权利要求范围之内。
虽然通过实施方式描绘了本说明书，本领域普通技术人员知道，本说明书有许多变形和变化而不脱离本说明书的精神，希望所附的权利要求包括这些变形和变化而不脱离本说明书的精神。
Claims (12)
1.一种判断储层流体趋势与压实趋势的方法，其特征在于，包括：
获取目标区域的叠前地震道集；
根据目标区域的叠前地震道集获得纵波速度、横波速度和密度数据体；
利用所述纵波速度、所述横波速度和所述密度数据体获得伪衰减属性体；
根据所述伪衰减属性体得到基于旋转角度为参数的正交属性因子表达式，基于岩石物理分析获得旋转角度值，将所述旋转角度值代入所述正交属性因子表达式获得对应正交属性因子；
基于正交属性因子判断储层流体趋势和压实趋势；
其中，所述伪衰减属性体包括纵波对应的伪衰减属性和横波对应的伪衰减属性；
伪衰减属性体由以下公式计算得到：
其中，FQP、FQS分别为纵横波对应的伪衰减属性，Vp为纵波速度，Vs为横波速度，ρ为密度。
2.如权利要求1所述的方法，其特征在于，根据目标区域的叠前地震道集获得纵波速度、横波速度和密度数据体的步骤包括：
对所述目标区域的叠前地震道集进行叠加处理；
根据叠加的叠前地震道集反演得到纵波速度、横波速度和密度数据体。
3.如权利要求1所述的方法，其特征在于，所述正交属性因子包括压实趋势因子和流体趋势因子。
4.如权利要求1所述的方法，其特征在于，根据所述伪衰减属性体得到基于旋转角度为参数的正交属性因子表达式的步骤包括：
将所述伪衰减属性体进行正交旋转运算，得到基于旋转角度为参数的正交属性因子表达式。
5.如权利要求3所述的方法，其特征在于，所述正交属性因子表达式为：
其中，FQP为纵波对应的伪衰减属性；FQS为横波对应的伪衰减属性；C为压实趋势因子，C＝FQSsin a-FQPcos a；F为流体趋势因子，F＝FQScos a+FQPsin a；
为正交旋转矩阵；a为旋转角；
根据岩石物理参数确定旋转角度值，旋转角度计算公式如下所示：
其中，
分别为干燥岩石骨架的纵波速度和横波速度。
6.一种判断储层流体趋势与压实趋势的装置，其特征在于，包括：
叠前地震道集获取单元，用于获取目标区域的叠前地震道集；
反演单元，用于根据目标区域的叠前地震道集获得纵波速度、横波速度和密度数据体；
伪衰减属性体获取单元，用于利用所述纵波速度、所述横波速度和所述密度数据体获得伪衰减属性体；
正交属性因子获取单元，用于根据所述伪衰减属性体得到基于旋转角度为参数的正交属性因子表达式，基于岩石物理分析获得旋转角度值，将所述旋转角度值代入所述正交属性因子表达式获得对应正交属性因子；
判断单元，用于基于正交属性因子判断储层流体趋势和压实趋势；
其中，所述伪衰减属性体获取单元获取的伪衰减属性体包括纵波对应的伪衰减属性和横波对应的伪衰减属性；
伪衰减属性体由以下公式计算得到：
其中，FQP、FQS分别为纵横波对应的伪衰减属性，Vp为纵波速度，Vs为横波速度，ρ为密度。
7.如权利要求6所述的装置，其特征在于，所述反演单元包括：
叠加模块，用于对所述目标区域的叠前地震道集进行叠加处理；
反演计算模块，用于根据叠加的叠前地震道集反演得到纵波速度、横波速度和密度数据体。
8.如权利要求6所述的装置，其特征在于，所述正交属性因子获取单元获取的正交属性因子包括压实趋势因子和流体趋势因子。
9.如权利要求6所述的装置，其特征在于，所述正交属性因子获取单元将所述伪衰减属性体进行正交旋转运算，得到基于旋转角度为参数的正交属性因子表达式。
10.如权利要求8所述的装置，其特征在于，所述正交属性因子获取单元获取的正交属性因子表达式为：
其中，FQP为纵波对应的伪衰减属性；FQS为横波对应的伪衰减属性；C为压实趋势因子，C＝FQSsin a-FQPcos a；F为流体趋势因子，F＝FQScos a+FQPsin a；
为正交旋转矩阵；a为旋转角；
根据岩石物理参数确定旋转角度值，旋转角度计算公式如下所示：
其中，
分别为干燥岩石骨架的纵波速度和横波速度。
11.一种电子设备，包括存储器、处理器以及存储在所述存储器上并可在所述处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1～5任意一项权利要求所述的判断储层流体趋势与压实趋势的方法。
12.一种可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被执行时实现权利要求1～5任意一项权利要求所述的判断储层流体趋势与压实趋势的方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811469482.XA CN111239811B (zh) | 2018-11-28 | 2018-11-28 | 一种判断储层流体趋势与压实趋势的方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811469482.XA CN111239811B (zh) | 2018-11-28 | 2018-11-28 | 一种判断储层流体趋势与压实趋势的方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111239811A CN111239811A (zh) | 2020-06-05 |
| CN111239811B true CN111239811B (zh) | 2022-07-05 |
Family
ID=70869356
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811469482.XA Active CN111239811B (zh) | 2018-11-28 | 2018-11-28 | 一种判断储层流体趋势与压实趋势的方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111239811B (zh) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| AU2001275708A1 (en) * | 2000-07-27 | 2002-05-09 | Schlumberger Holdings Limited | A method of processing seismic data |
| CN102749645A (zh) * | 2012-03-14 | 2012-10-24 | 中国石油天然气股份有限公司 | 一种利用角度阻抗梯度进行储层烃类检测的方法及装置 |
| CN105182424A (zh) * | 2015-08-03 | 2015-12-23 | 中国石油天然气股份有限公司 | 一种基于斑块饱和模型定量预测储层孔隙度的方法和装置 |
| CN106842310A (zh) * | 2015-12-04 | 2017-06-13 | 中国石油化工股份有限公司 | 叠前地震四参数同步反演方法 |
| CN107831543A (zh) * | 2017-09-29 | 2018-03-23 | 中国石油化工股份有限公司 | 叠前地震反演方法及系统 |
| CN108693557A (zh) * | 2018-04-27 | 2018-10-23 | 中国石油天然气集团有限公司 | 一种地震数据处理方法及装置 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB0018480D0 (en) * | 2000-07-27 | 2000-09-13 | Geco Prakla Uk Ltd | A method of processing surface seismic data |
| US8219319B2 (en) * | 2009-12-18 | 2012-07-10 | Chevron U.S.A. Inc. | Workflow for petrophysical and geophysical formation evaluation of wireline and LWD log data |
| CN103135135B (zh) * | 2013-02-26 | 2016-01-06 | 中国石油天然气股份有限公司 | 一种基于疏松砂岩模型进行烃类定量预测的方法和装置 |
| WO2018102732A1 (en) * | 2016-12-01 | 2018-06-07 | Schlumberger Technology Corporation | Coupled reservoir-geomechanical models using compaction tables |
-
2018
- 2018-11-28 CN CN201811469482.XA patent/CN111239811B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| AU2001275708A1 (en) * | 2000-07-27 | 2002-05-09 | Schlumberger Holdings Limited | A method of processing seismic data |
| CN102749645A (zh) * | 2012-03-14 | 2012-10-24 | 中国石油天然气股份有限公司 | 一种利用角度阻抗梯度进行储层烃类检测的方法及装置 |
| CN105182424A (zh) * | 2015-08-03 | 2015-12-23 | 中国石油天然气股份有限公司 | 一种基于斑块饱和模型定量预测储层孔隙度的方法和装置 |
| CN106842310A (zh) * | 2015-12-04 | 2017-06-13 | 中国石油化工股份有限公司 | 叠前地震四参数同步反演方法 |
| CN107831543A (zh) * | 2017-09-29 | 2018-03-23 | 中国石油化工股份有限公司 | 叠前地震反演方法及系统 |
| CN108693557A (zh) * | 2018-04-27 | 2018-10-23 | 中国石油天然气集团有限公司 | 一种地震数据处理方法及装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111239811A (zh) | 2020-06-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107679358B (zh) | 一种确定储层渗透率的方法及装置 | |
| CN107390268B (zh) | 低渗透砂砾岩成岩圈闭的识别方法及装置 | |
| CN108252710B (zh) | 一种确定油藏含水饱和度的方法及装置 | |
| RU2686720C1 (ru) | Способ и устройство для определения расходов компонентов многофазного флюида | |
| CN108008444B (zh) | 一种确定岩石物理模型的方法及装置 | |
| CN108828668B (zh) | 一种叠前时间偏移数据处理方法及装置 | |
| CN104391342A (zh) | 基于裂缝孔隙度反演定量预测储层渗透性的方法及装置 | |
| CN108897036B (zh) | 一种地震数据处理方法及装置 | |
| CN106202813B (zh) | 一种计算地层原油粘度的方法及装置 | |
| CN108897040B (zh) | 一种确定陆检与海底之间的耦合性的方法及装置 | |
| CN104453882A (zh) | 一种确定致密储层原油可动储量的方法及装置 | |
| CN109612900B (zh) | 一种储层岩石渗透率预测方法、装置及存储介质 | |
| CN104459774B (zh) | 一种地质岩性差异识别方法及系统 | |
| CN111323814B (zh) | 基于岩石物理模板定量确定砂岩储层含水饱和度的方法 | |
| CN106960264B (zh) | 低渗透砂砾岩地层压力预测方法及装置 | |
| CN111239811B (zh) | 一种判断储层流体趋势与压实趋势的方法及装置 | |
| CN111239832B (zh) | 3d3c-vsp成像处理方法、装置及设备 | |
| Ghedan et al. | Wettability profile of a thick carbonate reservoir by the new rise in core wettability characterization method | |
| CN112049623B (zh) | 一种中基性火山岩束缚水出水的判定方法和装置 | |
| CN108318928B (zh) | 一种大地电磁数据的处理方法及装置 | |
| CN108505991B (zh) | 一种确定岩心中油的采出程度的方法及装置 | |
| CN106908842B (zh) | 一种确定单井沉积微相曲线的方法及装置 | |
| CN108108512B (zh) | 一种储层岩相的表征方法及装置 | |
| CN107605472B (zh) | 一种确定油藏采收率的方法及装置 | |
| CN109063228B (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 |