CN111751900B - 一种流体识别方法和装置 - Google Patents
一种流体识别方法和装置 Download PDFInfo
- Publication number
- CN111751900B CN111751900B CN201910247471.5A CN201910247471A CN111751900B CN 111751900 B CN111751900 B CN 111751900B CN 201910247471 A CN201910247471 A CN 201910247471A CN 111751900 B CN111751900 B CN 111751900B
- Authority
- CN
- China
- Prior art keywords
- work area
- stratum
- reservoir
- determining
- formation
- 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 223
- 238000000034 method Methods 0.000 title claims abstract description 77
- 238000005259 measurement Methods 0.000 claims abstract description 127
- 239000011435 rock Substances 0.000 claims abstract description 103
- 229930195733 hydrocarbon Natural products 0.000 claims abstract description 51
- 150000002430 hydrocarbons Chemical class 0.000 claims abstract description 49
- 239000004215 Carbon black (E152) Substances 0.000 claims abstract description 36
- 239000011148 porous material Substances 0.000 claims abstract description 26
- 230000015572 biosynthetic process Effects 0.000 claims description 174
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 99
- 238000012937 correction Methods 0.000 claims description 29
- 239000011159 matrix material Substances 0.000 claims description 21
- 229920006395 saturated elastomer Polymers 0.000 claims description 17
- 238000004590 computer program Methods 0.000 claims description 16
- 230000000704 physical effect Effects 0.000 claims description 8
- 238000005056 compaction Methods 0.000 claims description 7
- 230000010287 polarization Effects 0.000 claims description 5
- 230000035945 sensitivity Effects 0.000 abstract description 12
- 238000010586 diagram Methods 0.000 description 19
- 230000008859 change Effects 0.000 description 14
- 230000014509 gene expression Effects 0.000 description 8
- 238000012986 modification Methods 0.000 description 8
- 230000004048 modification Effects 0.000 description 8
- 230000008569 process Effects 0.000 description 6
- 238000010276 construction Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 230000006870 function Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 230000004044 response Effects 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 238000001514 detection method Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 125000001183 hydrocarbyl group Chemical group 0.000 description 2
- 238000011545 laboratory measurement Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- DMSMPAJRVJJAGA-UHFFFAOYSA-N benzo[d]isothiazol-3-one Chemical compound C1=CC=C2C(=O)NSC2=C1 DMSMPAJRVJJAGA-UHFFFAOYSA-N 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000004088 simulation Methods 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
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)
Abstract
本发明提供了一种流体识别方法和装置，该方法包括：获取工区地层的测井数据和物性参数；根据工区地层的物性参数，确定工区地层的岩石物理模板；根据工区地层的岩石物理模板和测井数据，确定工区地层中储层的实测点的孔隙度；根据工区地层中储层的实测点的孔隙度和工区地层的物性参数，确定工区地层中储层的实测点的流体因子；根据工区地层中储层的实测点的流体因子，确定工区地层中烃类的相对含量。本发明提出的基于岩石物理模板的新型流体因子有效地消除了孔隙度对流体因子敏感性的影响，能够较为准确的反映储层中流体饱和度大小，能够准确识别高饱和低孔致密储层中烃类的相对含量。
Description
技术领域
本发明涉及石油地球物理勘探领域，具体涉及一种流体识别方法和装置。
背景技术
储层流体识别是指综合利用多种地球物理资料对储层含流体(烃类和水)的特性进行识别和表征的过程。早期，储层流体识别主要依靠地震振幅属性实现，“亮点”、“暗点”、“平点”、“相位反转”等识别技术在当时广为应用。随后，根据反射系数随偏移距发生变化的原理，发展了一系列基于AVO响应特征的烃类(油和气)识别方法。目前基于叠前道集的AVO分析技术仍是储层流体识别的主要方法之一。
Smith等(1987)首次提出“流体因子”的概念，用于储层流体识别。而后诸多学者对其进行研究，提出了一系列流体因子的构建方法。这些流体因子可大致分为独立流体因子和组合流体因子两种类型。独立流体因子的物理意义更为明确，是孔隙流体变化在岩石弹性性质上产生的响应。随着研究储层的复杂化，独立流体因子无法满足储层流体识别工作的需求，理论模拟和实际数据分析表明，独立流体因子只能定性表征储层流体的性质，无法准确地描述岩石孔隙中烃类和水的相对含量。在实际生产中，饱和度过低的含气储层不具有商业价值，故预测储层流体中烃类的相对含量的工作对实际勘探开发过程有着重要意义；因此通过组合运算得到的流体因子开始发挥重要作用。较为常用的有物性参数的简单运算，以及Russell流体因子、流体因子角、泊松阻抗等更为复杂的流体因子。
岩石物理模板是连接储层物性参数与其弹性响应的桥梁，对储层流体识别有着重要的指导意义。岩石物理模板(RPT)由
和Avseth(2004)首次提出，是根据测井和地震反演数据对储层进行岩性识别和孔隙流体检测的重要手段。此后，基于岩石物理模板的流体因子构建方法不断涌现。这些方法具有简单、直观、物理意义明确等独特优势，对孔隙中的流体变化有较高的敏感性，能够有效地分辩出含烃地层。
现有的流体因子构建方法都没有考虑孔隙度的影响。同一工区同一目的层段内的岩石孔隙度并非均匀分布，而是会在横向和纵向上发生变化。若对具有不同孔隙度的岩石采用相同的流体检测标准，就会错将孔隙度变化引起的弹性响应归结为孔隙流体的变化，进而无法准确地反映岩石孔隙中烃类和水的相对含量，这一问题在高饱和低孔致密储层的流体识别过程中尤为明显，因此现有的流体因子构建方法受孔隙度的影响较大，其敏感性随孔隙度的减小而降低，无法有效识别高饱和的低孔致密储层中烃类的相对含量。
发明内容
本发明实施例提供一种流体识别方法，通过校正消除孔隙度对流体因子的影响，能够准确识别高饱和低孔致密储层中烃类的相对含量，该流体识别方法包括：
获取工区地层的测井数据和物性参数；
根据工区地层的物性参数，确定工区地层的岩石物理模板；
根据工区地层的岩石物理模板和测井数据，确定工区地层中储层的实测点的孔隙度；
根据工区地层中储层的实测点的孔隙度和工区地层的物性参数，确定工区地层中储层的实测点的流体因子；
根据工区地层中储层的实测点的流体因子，确定工区地层中烃类的相对含量；
根据工区地层中储层的实测点的孔隙度和工区地层的物性参数，确定工区地层中储层的实测点的流体因子，包括：
根据工区地层中储层的实测点的孔隙度和工区地层的物性参数，确定工区地层中储层的实测点的孔隙度校正因子；
根据工区地层中储层的实测点的孔隙度校正因子和孔隙度，确定工区地层中储层的实测点的流体因子；
按如下方式，根据工区地层中储层的实测点的孔隙度校正因子和孔隙度，确定工区地层中储层的实测点的流体因子：
其中，NFF为流体因子；
为孔隙度；
和
分别表示与工区地层中储层的实测点孔隙度相同的完全饱水岩石的纵波模量和剪切模量；Mmea为工区地层实测点的纵波模量；Gmea为工区地层实测点的剪切模量。
本发明实施例还提供一种流体识别装置，包括：
数据获取模块，用于获取工区地层的测井数据和物性参数；
岩石物理模板确定模块，用于根据工区地层的物性参数，确定工区地层的岩石物理模板；
孔隙度确定模块，用于根据工区地层的岩石物理模板和测井数据，确定工区地层中储层的实测点的孔隙度；
流体因子确定模块，用于根据工区地层中储层的实测点的孔隙度和工区地层的物性参数，确定工区地层中储层的实测点的流体因子；
烃类相对含量确定模块，用于根据工区地层中储层的实测点的流体因子，确定工区地层中烃类的相对含量；
根据工区地层中储层的实测点的孔隙度和工区地层的物性参数，确定工区地层中储层的实测点的流体因子，包括：
根据工区地层中储层的实测点的孔隙度和工区地层的物性参数，确定工区地层中储层的实测点的孔隙度校正因子；
根据工区地层中储层的实测点的孔隙度校正因子和孔隙度，确定工区地层中储层的实测点的流体因子；
按如下方式，根据工区地层中储层的实测点的孔隙度校正因子和孔隙度，确定工区地层中储层的实测点的流体因子：
其中，NFF为流体因子；
为孔隙度；
和
分别表示与工区地层中储层的实测点孔隙度相同的完全饱水岩石的纵波模量和剪切模量；Mmea为工区地层实测点的纵波模量；Gmea为工区地层实测点的剪切模量。
本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现前述流体识别方法。
本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行实现前述流体识别方法的计算机程序。
本发明提供的一种流体识别方法和装置，克服了目前常用的流体因子的敏感性大多受到孔隙度的影响、只能定性的辨别储层和非储层，而无法定量表征储层中烃类的相对含量的问题。本发明提出的基于岩石物理模板的新型流体因子有效地消除了孔隙度对流体因子敏感性的影响，能够较为准确的反映储层中流体饱和度大小，能够准确识别高饱和低孔致密储层中烃类的相对含量。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。在附图中：
图1为本发明实施例一种流体识别方法的示意图。
图2为本发明实施例一种流体识别方法基于岩石物理模板构建去除孔隙度影响的新型流体因子的流程图。
图3为本发明实施例一种流体识别方法的岩石物理模板示意图。
图4为本发明实施例一种流体识别方法的完全水饱和曲线与投影点的相对位置关系图。
图5为本发明实施例一种流体识别方法构建的新型流体因子NFF随孔隙度和流体饱和度的变化特征图。
图6为传统流体因子UFF随孔隙度和流体饱和度的变化特征图。
图7为本发明实施例一种流体识别方法构建的新型流体因子NFF在模拟砂泥岩剖面中的应用效果图。
图8为苏里格气田某井的实际测井数据在常规岩石物理模板中的投影图。
图9为本发明实施例一种流体识别方法构建的新型流体因子NFF和传统流体因子UFF应用于苏里格气田某井测井数据的结果图。
图10为以本发明实施例一种流体识别方法构建的新型流体因子NFF和反演孔隙度作为横纵坐标构建的新型岩石物理模板图。
图11为本发明实施例一种流体识别装置的示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白，下面结合附图对本发明实施例做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。
本发明实施例提供一种流体识别方法，通过校正消除孔隙度对流体因子的影响，能够准确识别高饱和低孔致密储层中烃类的相对含量，如图1本发明实施例一种流体识别方法的示意图所示，该流体识别方法包括：
步骤101，获取工区地层的测井数据和物性参数；
步骤102，根据工区地层的物性参数，确定工区地层的岩石物理模板；
步骤103，根据工区地层的岩石物理模板和测井数据，确定工区地层中储层的实测点的孔隙度；
步骤104，根据工区地层中储层的实测点的孔隙度和工区地层的物性参数，确定工区地层中储层的实测点的流体因子；
步骤105，根据工区地层中储层的实测点的流体因子，确定工区地层中烃类的相对含量。
本发明实施例提供的流体识别方法，克服了目前常用的流体因子的敏感性大多受到孔隙度的影响、只能定性的辨别储层和非储层，而无法定量表征储层中烃类的相对含量的问题。本发明提出的基于岩石物理模板的新型流体因子有效地消除了孔隙度对流体因子敏感性的影响，能够较为准确的反映储层中流体饱和度大小，能够准确识别高饱和低孔致密储层中烃类的相对含量。
高饱和低孔致密储层是当今非常规油气勘探的重要对象，但在勘探过程中现有的技术和方法无法准确地反映岩石孔隙中烃类和水的相对含量，如图2本发明实施例一种流体识别方法基于岩石物理模板构建去除孔隙度影响的新型流体因子的流程图所示，因此本发明实施例提供一种流体识别方法：首先获取工区地层的测井数据和物性参数；然后根据工区地层的物性参数，确定工区地层的岩石物理模板；接着根据工区地层的岩石物理模板和测井数据，确定工区地层中储层的实测点的孔隙度；再根据工区地层中储层的实测点的孔隙度和工区地层的物性参数，确定工区地层中储层的实测点的流体因子；最终根据工区地层中储层的实测点的流体因子，确定工区地层中烃类的相对含量。
在具体实施获取工区地层的测井数据和物性参数时，在一个实施例中，可以通过测井、实验室测量等方法，获得工区地层的测井数据和物性参数。
在获取到工区地层的测井数据和物性参数后，进一步需要确定工区地层的岩石物理模板。在一个实施例中，可以根据工区地层的物性参数，确定工区地层的岩石物理模板，可以包括：首先根据前述获取到的工区地层的物性参数，确定工区地层的纵波模量和剪切模量；具体的是，根据工区地层的物性参数和岩石物理模型，计算工区地层孔隙度和饱和度变化范围内岩石的纵波模量和剪切模量，进一步的，在岩石物理模型中输入工区地层的物性参数，选定岩石物理模板的孔隙度变化范围和饱和度变化范围；进一步的，岩石物理模型采用Xu-White模型和Gassmann方程进行正演计算，岩石物理模板的孔隙度变化范围可以是0.1％-35％，变化步长为0.5％；岩石物理模板的饱和度变化范围可以是0-100％，变化步长为0.1％；然后将工区地层的纵波模量与剪切模量进行相比，将工区地层的纵波模量与剪切模量的比，确定为工区地层的纵横波模量比；最后根据工区地层的纵波模量和纵横波模量比，确定工区地层的岩石物理模板。
在确定工区地层的岩石物理模板的过程中，一个实施例中，所需要的工区地层的物性参数，可以包括：水的体积模量、气的体积模量、岩石基质的体积模量、岩石基质的剪切模量、与孔隙形状有关的极化因子等。
获取到需要工区地层的物性参数后，在一个实施例中，按如下方式，确定工区地层的纵波模量和剪切模量，进一步的，可以根据上述获取到的工区地层的物性参数，选定该工区地层中岩石孔隙度变化范围和饱和度的变化范围，确定适用于工区地层的岩石物理模板上所有网格点的纵波模量和剪切模量：
其中，M为工区地层的岩石物理模板上网格点的纵波模量，G为工区地层的岩石物理模板上网格点的剪切模量，Ksat为岩石物理模板上网格点对应流体饱和岩石的体积模量；Gdry为岩石物理模板上网格点对应干燥岩石的剪切模量；进一步的，Ksat和Gdry可以通过下面的Gassmann方程和Xu-White模型，根据工区地层的物性参数计算得出：
Kdry＝Kma(1-φ)p,
Gdry＝Gma(1-φ)q,
其中，Kdry为岩石物理模板上网格点对应干燥岩石骨架的体积模量，Kma为岩石物理模板上网格点对应岩石基质的体积模量；Gma为岩石物理模板上网格点对应岩石基质的剪切模量，φ为岩石物理模板上网格点对应的孔隙度(一般为在岩石物理模板构建的过程中，通过对工区地层的地址特征进行分析后选定的一个变化范围)，p和q为与孔隙形状有关的极化因子；Kf为岩石孔隙流体的体积模量，Kw为水的体积模量，Kg为气的体积模量，e为经验参数。
前述提到的确定适用于工区地层岩石物理模板的纵波模量和剪切模量的表达式、推导过程和符号为举例说明，本领域技术人员可以理解，在实施时还可以根据需要对上述公式进行一定形式的变形和添加其它的参数或数据，或者提供其它的具体公式，这些变化例均应落入本发明的保护范围。
确定了岩石物理模板的纵波模量和剪切模量后，进一步的需要将岩石物理模板的纵波模量与剪切模量进行相比，将岩石物理模板的纵波模量与剪切模量的比，确定为岩石物理模板的纵横波模量比M/G；在一个实施例中，根据由岩石物理模型计算的纵波模量和纵横波模量比，确定工区地层的岩石物理模板，可以包括：根据岩石物理模板点的纵波模量M和纵横波模量比M/G建立纵波模量和纵横波模量比交会图，如图3本发明实施例一种流体识别方法的岩石物理模板示意图所示，在纵波模量和纵横波模量比交会图上建立工区地层的岩石物理模板。
建立工区地层的岩石物理模板后，进一步的需要确定工区地层中储层的实测点的孔隙度。在一个实施例中，可以根据工区地层的岩石物理模板和测井数据，确定工区地层中储层的实测点的孔隙度，具体可以包括：根据工区地层的测井数据和岩石物理模板，确定完全水饱和曲线；根据工区地层的测井数据、岩石物理模板和完全水饱和曲线，确定储层判别因子；根据储层判别因子，判定工区地层中的储层和非储层；根据工区地层的岩石物理模板，确定工区地层中的储层的实测点的孔隙度。
具体实施前述的确定完全水饱和曲线时，在一个实施例中，需要根据工区地层的测井数据和岩石物理模板，确定完全水饱和曲线，具体可以包括：将工区地层的测井数据以工区地层实测点的纵波模量为横坐标，工区地层实测点的纵横波模量比为纵坐标，投影到工区地层的岩石物理模板上，确定岩石物理模板上网格数据点；进一步的，再根据工区地层的测井数据与工区地层的岩石物理模板上网格数据点的对应关系，调节工区地层的岩石物理模板的相关参数，校准工区地层的岩石物理模板，使工区地层的岩石物理模板与工区地层的测井数据匹配；从校准的工区地层的岩石物理模板得到完全水饱和网格数据点；进一步的，再根据完全水饱和网格数据点，确定完全水饱和曲线。
具体实施前述将工区地层的测井数据以纵波模量为横坐标，纵横波模量比为纵坐标，投影到工区地层的岩石物理模板上，确定岩石物理模板上网格数据点时，在一个实施例中，工区地层的测井数据包括：工区地层实测点的纵波速度、横波速度和密度；进一步的，可以按如下方式，根据从工区地层的测井数据中的纵波速度Vp、横波速度Vs、密度ρ，确定工区地层实测点的纵波模量和剪切模量：
其中，Mmea为工区地层实测点的纵波模量；Gmea为工区地层实测点的剪切模量；Vp为工区地层实测点的纵波速度；Vs为工区地层实测点的横波速度；ρ为工区地层实测点的密度。
前述提到的确定工区地层实测点的纵波模量和剪切模量的表达式和符号为举例说明，本领域技术人员可以理解，在实施时还可以根据需要对上述公式进行一定形式的变形和添加其它的参数或数据，或者提供其它的具体公式，这些变化例均应落入本发明的保护范围。
具体实施前述的根据完全水饱和网格数据点确定完全水饱和曲线时，在一个实施例中，可以包括，拾取工区地层的岩石物理模板上的完全水饱和数据网格点，将完全水饱和数据网格点采用对数域多项式进行拟合，确定完全水饱和曲线。进一步的，可以从工区地层的岩石物理模板中提取完全水饱和的网格数据点；横坐标用Mw表示，纵坐标用Mw/Gw表示，用对数域多项式拟合从工区地层的岩石物理模板中提取的所有完全水饱和网格数据点，确定完全水饱和曲线；在具体实施时，可以采用对数域多项式进行最小二乘拟合；进一步的，在一个实施例中可以按照如下方式，确定完全水饱和曲线：
fw(x)＝exp(a·ln4x+b·ln3x+c·ln2x+d·lnx+e)，
其中，fw(x)为完全水饱和曲线；a,b,c,d,e为拟合系数。
前述提到的确定完全水饱和曲线的表达式和符号为举例说明，本领域技术人员可以理解，在实施时还可以根据需要对上述公式进行一定形式的变形和添加其它的参数或数据，或者提供其它的具体公式，这些变化例均应落入本发明的保护范围。
在确定了完全水饱和曲线后，进一步的需要确定储层判别因子。在一个实施例中，可以根据工区地层的测井数据、岩石物理模板和完全水饱和曲线，确定储层判别因子，具体可以包括：将工区地层的测井数据投影到岩石物理模板上，得到投影点；进一步的根据投影点与完全水饱和曲线的相对位置，利用工区地层实测点的纵波模量、工区地层实测点的剪切模量和岩石物理模板的完全水饱和曲线表达式，确定储层判别因子。
前述的根据投影点与完全水饱和曲线的相对位置，利用工区地层实测点的纵波模量、工区地层实测点的剪切模量和岩石物理模板的完全水饱和曲线，确定储层判别因子，在一个实施例中，可以包括：根据工区地层的测井数据，确定工区地层中的实测点的纵波模量和剪切模量；进一步的，按照如下方式，确定储层判别因子：
RF＝fw(Mmea)-Mmea/Gmea，
其中，RF为储层判别因子；fw(x)为完全水饱和曲线；Mmea为工区地层中实测点的纵波模量；Gmea为工区地层中实测点的剪切模量。
前述提到的确定储层判别因子的表达式和符号为举例说明，本领域技术人员可以理解，在实施时还可以根据需要对上述公式进行一定形式的变形和添加其它的参数或数据，或者提供其它的具体公式，这些变化例均应落入本发明的保护范围。
在确定了储层判别因子后，进一步的需要判断工区地层中的储层和非储层。在一个实施例中，可以根据储层判别因子，判定工区地层中的储层和非储层，具体可以包括：若储层判别因子的值大于等于零，则工区地层的中的实测点对应的地层是储层；若储层判别因子的值小于零，则工区地层中的实测点对应的地层是非储层。
图4是本发明实施例一种流体识别方法的完全水饱和曲线与投影点的相对位置关系图，位于完全水饱和曲线之上投影点的含水饱和度较大，可视为非储层；其下投影点的含水饱和度较小，可视为储层。
在确定了工区地层中的储层后，进一步的需要确定工区地层中储层的实测点的孔隙度。在一个实施例中，可以根据工区地层的岩石物理模板，确定工区地层中储层的实测点的孔隙度，具体可以包括：将工区地层中储层的实测点采用最小二乘方法搜索在岩石物理模板上距其最近的网格点，得到工区地层中储层的实测点的孔隙度。
在得到工区地层中储层的实测点的孔隙度后，进一步的需要确定工区地层中储层的实测点的流体因子。在一个实施例中，可以根据工区地层中储层的实测点的孔隙度和工区地层的物性参数，确定工区地层中储层的实测点的流体因子，具体可以包括：根据工区地层中储层的实测点的孔隙度和工区地层的物性参数，确定工区地层中储层的实测点的孔隙度校正因子；根据工区地层中储层的实测点的孔隙度校正因子和孔隙度，确定工区地层中储层的实测点的流体因子。
前述通过实验室测量、测井等方法获取的工区地层的物性参数，可以包括：水的体积模量、气的体积模量、岩石基质的体积模量和剪切模量等；在具体实施时，工区地层的物性参数包括：水的体积模量、气的体积模量、岩石基质的体积模量、岩石骨架的体积模量、岩石基质的剪切模量和岩石骨架的剪切模量等。
在具体实施前述的根据工区地层中储层的实测点的孔隙度和工区地层的物性参数，确定工区地层中储层的实测点的孔隙度校正因子时，一个实施例中可以包括：按如下方式确定工区地层中储层的实测点的孔隙度校正因子：
其中，
为孔隙度校正因子；Kw为水的体积模量；Kg为气的体积模量；Kma为岩石基质的体积模量；Kdry为岩石骨架的体积模量；Gma为岩石基质的剪切模量；Gdry为岩石骨架的剪切模量；
为工区地层中储层的实测点的孔隙度。
前述提到的确定工区地层中储层的实测点的孔隙度校正因子的表达式和符号为举例说明，本领域技术人员可以理解，在实施时还可以根据需要对上述公式进行一定形式的变形和添加其它的参数或数据，或者提供其它的具体公式，这些变化例均应落入本发明的保护范围。
在确定了工区地层中储层的实测点的孔隙度校正因子后，进一步的需要确定工区地层中储层的实测点的流体因子。在一个实施例中，可以按如下方式，根据工区地层中储层的实测点的孔隙度校正因子和孔隙度，确定工区地层中储层的实测点的流体因子：
其中，NFF为流体因子；
为工区地层中储层的实测点的孔隙度；
和
分别表示与工区地层中储层的实测点孔隙度相同的完全饱水岩石的纵波模量和剪切模量。
前述提到的确定工区地层中储层的实测点的流体因子的表达式和符号为举例说明，本领域技术人员可以理解，在实施时还可以根据需要对上述公式进行一定形式的变形和添加其它的参数或数据，或者提供其它的具体公式，这些变化例均应落入本发明的保护范围。
前述的
和
在一个实施例中，可以按照如下方式，将工区地层中储层的实测点的孔隙度，代入到完全水饱和情形下的Gassmann方程中计算，确定完全饱水岩石的纵波模量和剪切模量：
其中，
为完全饱水岩石的纵波模量；
为完全饱水岩石的剪切模量；Kw为水的体积模量；Kg为气的体积模量；Kma为岩石基质的体积模量；Kdry为岩石骨架的体积模量；Gma为岩石基质的剪切模量；Gdry为岩石骨架的剪切模量；
为工区地层中储层的实测点的孔隙度。
前述提到的确定完全饱水岩石的纵波模量和剪切模量的表达式和符号为举例说明，本领域技术人员可以理解，在实施时还可以根据需要对上述公式进行一定形式的变形和添加其它的参数或数据，或者提供其它的具体公式，这些变化例均应落入本发明的保护范围。
图5是本发明实施例一种流体识别方法构建的新型流体因子NFF随孔隙度和流体饱和度的变化特征图；图6是传统流体因子UFF随孔隙度和流体饱和度的变化特征图，由图6可以得出，传统流体因子UFF受孔隙度影响明显，其敏感性随孔隙度的减小而降低，无法有效识别高饱和的致密储层；而在图5中，新型流体因子NFF基本不受孔隙度影响，对各个孔隙度的岩石具有统一的检测标准，能够较为准确地表征流体饱和度的相对大小。
在确定了工区地层中储层的实测点的流体因子后，进一步的，需要确定工区地层中烃类的相对含量。在一个实施例中，可以根据工区地层中储层的实测点的流体因子，确定工区地层中烃类的相对含量，具体可以包括：
工区地层中储层的实测点的流体因子的取值介于0到1之间；
工区地层中储层的实测点的流体因子的取值为1时，表征工区地层完全含水；
工区地层中储层的实测点的流体因子的取值为0时，表征工区地层完全含烃；
工区地层中储层的实测点的流体因子的取值越小，表征工区地层中水的相对含量越小，烃类的相对含量越高；
工区地层中储层的实测点的流体因子的取值越大，表征工区地层中水的相对含量越高，烃类的相对含量越小。
图7是本发明实施例一种流体识别方法构建的新型流体因子NFF在模拟砂泥岩剖面中的应用效果图，在该剖面中包含两套具有相同流体饱和度(Sg＝0.7)的含气砂岩，上部砂岩的孔隙度为0.35，下部砂岩的孔隙度为0.1，此外还有一套高孔(Por＝0.35)的完全含水砂岩。通过对比可以发现，应用本发明实施例构建的新型流体因子NFF与实测流体饱和度更为贴合，能够解决传统流体因子UFF无法有效识别高饱和致密储层的问题。
图8是苏里格气田某井的实际测井数据在常规岩石物理模板中的投影图，可见常规岩石物理模板符合实际测井数据的规律，完全水饱和曲线之上为泥岩区，其附近为含水砂岩区，其下为含气砂岩区。说明选用的Xu-White模型和Brie流体饱和方式能够表征该工区的地质特征和岩石物理特性。
图9是本发明实施例一种流体识别方法构建的新型流体因子NFF和传统流体因子UFF应用于苏里格气田某井测井数据的结果图。从图中可以得出，第4栏中的传统流体因子在小孔情况下的数值偏低，不能有效反映岩石孔隙中的流体相对含量，若以此为标准进行流体识别，势必会忽略高饱和的致密储层；而第4栏中的新型流体因子有效地改进了传统流体因子的缺陷，在全井段内都能够较为准确地反映储层流体的相对含量，适用于高饱和致密储层流体识别工作。
图10是以本发明实施例一种流体识别方法构建的新型流体因子NFF和反演孔隙度作为横纵坐标构建的新型岩石物理模板图，相当于对原有岩石物理模板进行坐标旋转。该模板中数据点的横坐标介于0到1之间，表征岩石孔隙中的水的相对含量。横坐标值为1时储层完全含有水，横坐标为0时储层完全含烃，其间横坐标越大，表征储层含水量越高，含烃量越低。
本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现前述流体识别方法。
本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行实现前述流体识别方法的计算机程序。
本发明实施例中还提供了一种流体识别装置，如下面的实施例所述。由于该装置解决问题的原理与一种流体识别方法相似，因此该装置的实施可以参见一种流体识别方法的实施，重复之处不再赘述。
如图11所示的本发明实施例一种流体识别装置的示意图所示，本发明实施例还提供一种流体识别装置，包括：
数据获取模块1101，用于获取工区地层的测井数据和物性参数；
岩石物理模板确定模块1102，用于根据工区地层的物性参数，确定工区地层的岩石物理模板；
孔隙度确定模块1103，用于根据工区地层的岩石物理模板和测井数据，确定工区地层中储层的实测点的孔隙度；
流体因子确定模块1104，用于根据工区地层中储层的实测点的孔隙度和工区地层的物性参数，确定工区地层中储层的实测点的流体因子；
烃类相对含量确定模块1105，用于根据工区地层中储层的实测点的流体因子，确定工区地层中烃类的相对含量。
一个实施例中，岩石物理模板确定模块，用于根据工区地层的物性参数，确定工区地层的岩石物理模板，包括：
根据工区地层的物性参数，确定工区地层的纵波模量和剪切模量；
将工区地层的纵波模量与剪切模量的比，确定为工区地层的纵横波模量比；
根据工区地层的纵波模量和纵横波模量比，确定工区地层的岩石物理模板。
一个实施例中，工区地层的物性参数，包括：水的体积模量、气的体积模量、岩石基质的体积模量、岩石基质的剪切模量、与孔隙形状有关的极化因子等。
一个实施例中，岩石物理模板确定模块，用于按如下方式，根据工区地层的物性参数，确定工区地层的纵波模量和剪切模量：
G＝Gdry.
其中，M为工区地层的纵波模量；G为工区地层的剪切模量，Ksat为岩石物理模板上网格点对应流体饱和岩石的体积模量；Gdry为岩石物理模板上网格点对应干燥岩石的剪切模量，Ksat和Gdry可以可通过下面的Gassmann方程和Xu-White模型算出：
Kdry＝Kma(1-φ)p,
Gdry＝Gma(1-φ)q,
其中，Kdry为岩石物理模板上网格点对应干燥岩石骨架的体积模量和剪切模量，Kma为岩石物理模板上网格点对应岩石基质的体积模量；Gma为岩石物理模板上网格点对应岩石基质的剪切模量，φ为岩石物理模板上网格点对应的孔隙度，p和q为与孔隙形状有关的极化因子；Kf为岩石孔隙流体的体积模量，Kw为水的体积模量，Kg为气的体积模量，e为经验参数。
一个实施例中，岩石物理模板确定模块，用于根据工区地层的纵波模量和纵横波模量比，确定工区地层的岩石物理模板，包括：
在纵波模量和纵横波模量比交会图上建立工区地层的岩石物理模板。
一个实施例中，孔隙度确定模块，用于根据工区地层的岩石物理模板和测井数据，确定工区地层中储层的实测点的孔隙度，包括：
根据工区地层的测井数据和岩石物理模板，确定完全水饱和曲线；
根据工区地层的测井数据、岩石物理模板和完全水饱和曲线，确定储层判别因子；
根据储层判别因子，判定工区地层中的储层和非储层；
根据工区地层的岩石物理模板，确定工区地层中的储层的实测点的孔隙度。
一个实施例中，孔隙度确定模块，用于根据工区地层的测井数据和岩石物理模板，确定完全水饱和曲线，包括：
将工区地层的测井数据以工区地层实测点的纵波模量为横坐标，工区地层实测点的纵横波模量比为纵坐标，投影到工区地层的岩石物理模板上，确定岩石物理模板上网格数据点；
根据工区地层的测井数据与工区地层的岩石物理模板上网格数据点的对应关系，调节工区地层的岩石物理模板的相关参数，校准工区地层的岩石物理模板，使工区地层的岩石物理模板与工区地层的测井数据匹配，得到完全水饱和网格数据点；
根据完全水饱和网格数据点，确定完全水饱和曲线。
一个实施例中，孔隙度确定模块，用于将工区地层的测井数据以工区地层实测点的纵波模量为横坐标，工区地层实测点的纵横波模量比为纵坐标，投影到工区地层的岩石物理模板上，确定岩石物理模板上网格数据点，包括：
工区地层的测井数据包括：工区地层实测点的纵波速度、横波速度和密度；
按如下方式，根据工区地层的测井数据，确定工区地层实测点的纵波模量和剪切模量：
其中，Mmea为工区地层实测点的纵波模量；Gmea为工区地层实测点的剪切模量；Vp为工区地层实测点的纵波速度；Vs为工区地层实测点的横波速度；ρ为工区地层实测点的密度。
一个实施例中，孔隙度确定模块，用于根据完全水饱和网格数据点，确定完全水饱和曲线，包括：
拾取工区地层的岩石物理模板上的完全水饱和网格数据点，将完全水饱和网格数据点采用对数域多项式进行拟合，确定完全水饱和曲线。
一个实施例中，孔隙度确定模块，用于按如下方式，确定完全水饱和曲线：
fw(x)＝exp(a·ln4x+b·ln3x+c·ln2x+d·lnx+e)，
其中，fw(x)为完全水饱和曲线；a,b,c,d,e为拟合系数。
一个实施例中，孔隙度确定模块，用于根据工区地层的测井数据、岩石物理模板和完全水饱和曲线，确定储层判别因子，包括：
将工区地层的测井数据投影到岩石物理模板上，得到投影点；
根据投影点与完全水饱和曲线的相对位置，利用工区地层实测点的纵波模量、工区地层实测点的剪切模量和岩石物理模板的完全水饱和曲线，确定储层判别因子。
一个实施例中，孔隙度确定模块，用于根据投影点与完全水饱和曲线的相对位置，利用工区地层实测点的纵波模量、工区地层实测点的剪切模量和岩石物理模板的完全水饱和曲线，确定储层判别因子，包括：
根据工区地层的测井数据，确定工区地层中的实测点的纵波模量和剪切模量；
按照如下方式，确定储层判别因子：
RF＝fw(Mmea)-Mmea/Gmea，
其中，RF为储层判别因子；fw(x)为完全水饱和曲线；Mmea为工区地层中实测点的纵波模量；Gmea为工区地层中实测点的剪切模量。
一个实施例中，孔隙度确定模块，用于根据储层判别因子，判定工区地层中的储层和非储层，包括：
储层判别因子的值大于等于零，则工区地层的中的实测点对应的地层是储层；
储层判别因子的值小于零，则工区地层中的实测点对应的地层是非储层。
一个实施例中，孔隙度确定模块，用于根据工区地层的岩石物理模板，确定工区地层中储层的实测点的孔隙度，包括：
将工区地层中储层的实测点采用最小二乘方法搜索在岩石物理模板上距其最近的网格点，得到工区地层中储层的实测点的孔隙度。
一个实施例中，流体因子确定模块，用于根据工区地层中储层的实测点的孔隙度和工区地层的物性参数，确定工区地层中储层的实测点的流体因子，包括：
根据工区地层中储层的实测点的孔隙度和工区地层的物性参数，确定工区地层中储层的实测点的孔隙度校正因子；
根据工区地层中储层的实测点的孔隙度校正因子和孔隙度，确定工区地层中储层的实测点的流体因子。
一个实施例中，流体因子确定模块，用于根据工区地层中储层的实测点的孔隙度和工区地层的物性参数，确定工区地层中储层的实测点的孔隙度校正因子，包括：
按如下方式确定工区地层中储层的实测点的孔隙度校正因子：
其中，
为孔隙度校正因子；Kw为水的体积模量；Kg为气的体积模量；Kma为岩石基质的体积模量；Kdry为岩石骨架的体积模量；Gma为岩石基质的剪切模量；Gdry为岩石骨架的剪切模量；
为孔隙度。
一个实施例中，流体因子确定模块，用于按如下方式，根据工区地层中储层的实测点的孔隙度校正因子和孔隙度，确定工区地层中储层的实测点的流体因子：
其中，NFF为流体因子；
为孔隙度；
和
分别表示与工区地层中储层的实测点孔隙度相同的完全饱水岩石的纵波模量和剪切模量。
一个实施例中，烃类相对含量确定模块，用于根据工区地层中储层的实测点的流体因子，确定工区地层中烃类的相对含量，包括：
工区地层中储层的实测点的流体因子的取值介于0到1之间；
工区地层中储层的实测点的流体因子的取值为1时，表征工区地层完全含水；
工区地层中储层的实测点的流体因子的取值为0时，表征工区地层完全含烃；
工区地层中储层的实测点的流体因子的取值越小，表征工区地层中水的相对含量越小，烃类的相对含量越高；
工区地层中储层的实测点的流体因子的取值越大，表征工区地层中水的相对含量越高，烃类的相对含量越小。
综上，由于目前常用的传统流体因子的敏感性大多受到孔隙度的影响，通常只能定性的辨别储层和非储层，而无法定量表征储层中烃类的相对含量。若用传统流体因子直接进行烃类相对含量检测，常会因流体因子的敏感性在小孔隙度时较低，致使致密储层高饱和区域无法得到有效识别。因此本发明提供了一种流体识别方法和装置，克服了目前常用的流体因子的敏感性大多受到孔隙度的影响、只能定性的辨别储层和非储层，而无法定量表征储层中烃类的相对含量的问题。本发明提出的基于岩石物理模板的新型流体因子有效地消除了孔隙度对流体因子敏感性的影响，能够较为准确的反映储层中流体饱和度大小，能够准确识别高饱和低孔致密储层中烃类的相对含量，对当今较为普遍的致密储层流体识别有很好的应用效果。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (22)
1.一种流体识别方法，其特征在于，包括：
获取工区地层的测井数据和物性参数；
根据工区地层的物性参数，确定工区地层的岩石物理模板；
根据工区地层的岩石物理模板和测井数据，确定工区地层中储层的实测点的孔隙度；
根据工区地层中储层的实测点的孔隙度和工区地层的物性参数，确定工区地层中储层的实测点的流体因子；
根据工区地层中储层的实测点的流体因子，确定工区地层中烃类的相对含量；
根据工区地层中储层的实测点的孔隙度和工区地层的物性参数，确定工区地层中储层的实测点的流体因子，包括：
根据工区地层中储层的实测点的孔隙度和工区地层的物性参数，确定工区地层中储层的实测点的孔隙度校正因子；
根据工区地层中储层的实测点的孔隙度校正因子和孔隙度，确定工区地层中储层的实测点的流体因子；
按如下方式，根据工区地层中储层的实测点的孔隙度校正因子和孔隙度，确定工区地层中储层的实测点的流体因子：
其中，NFF为流体因子；
为孔隙度；
和
分别表示与工区地层中储层的实测点孔隙度相同的完全饱水岩石的纵波模量和剪切模量；Mmea为工区地层实测点的纵波模量；Gmea为工区地层实测点的剪切模量。
2.如权利要求1所述的流体识别方法，其特征在于，根据工区地层的物性参数，确定工区地层的岩石物理模板，包括：
根据工区地层的物性参数，确定工区地层的纵波模量和剪切模量；
将工区地层的纵波模量与剪切模量的比，确定为工区地层的纵横波模量比；
根据工区地层的纵波模量和纵横波模量比，确定工区地层的岩石物理模板。
3.如权利要求2所述的流体识别方法，其特征在于，所述工区地层的物性参数，包括：水的体积模量、气的体积模量、岩石基质的体积模量、岩石基质的剪切模量、与孔隙形状有关的极化因子。
4.如权利要求3所述的流体识别方法，其特征在于，按如下方式，根据工区地层的物性参数，确定工区地层的纵波模量和剪切模量：
G＝Gdry.
其中，M为工区地层的纵波模量；G为工区地层的剪切模量，Ksat为岩石物理模板上网格点对应流体饱和岩石的体积模量；Gdry为岩石物理模板上网格点对应干燥岩石的剪切模量。
5.如权利要求2所述的流体识别方法，其特征在于，根据工区地层的纵波模量和纵横波模量比，确定工区地层的岩石物理模板，包括：
在纵波模量和纵横波模量比交会图上建立工区地层的岩石物理模板。
6.如权利要求3所述的流体识别方法，其特征在于，根据工区地层的岩石物理模板和测井数据，确定工区地层中储层的实测点的孔隙度，包括：
根据工区地层的测井数据和岩石物理模板，确定完全水饱和曲线；
根据工区地层的测井数据、岩石物理模板和完全水饱和曲线，确定储层判别因子；
根据储层判别因子，判定工区地层中的储层和非储层；
根据工区地层的岩石物理模板，确定工区地层中的储层的实测点的孔隙度。
7.如权利要求6所述的流体识别方法，其特征在于，根据工区地层的测井数据和岩石物理模板，确定完全水饱和曲线，包括：
将工区地层的测井数据以工区地层实测点的纵波模量为横坐标，工区地层实测点的纵横波模量比为纵坐标，投影到工区地层的岩石物理模板上，确定岩石物理模板上网格数据点；
根据工区地层的测井数据与工区地层的岩石物理模板上网格数据点的对应关系，调节工区地层的岩石物理模板的相关参数，校准工区地层的岩石物理模板，使工区地层的岩石物理模板与工区地层的测井数据匹配，得到完全水饱和网格数据点；
根据完全水饱和网格数据点，确定完全水饱和曲线。
8.如权利要求7所述的流体识别方法，其特征在于，
将工区地层的测井数据以工区地层实测点的纵波模量为横坐标，工区地层实测点的纵横波模量比为纵坐标，投影到工区地层的岩石物理模板上，确定岩石物理模板上网格数据点，包括：
工区地层的测井数据包括：工区地层实测点的纵波速度、横波速度和密度；
按如下方式，根据工区地层的测井数据，确定工区地层实测点的纵波模量和剪切模量：
其中，Mmea为工区地层实测点的纵波模量；Gmea为工区地层实测点的剪切模量；Vp为工区地层实测点的纵波速度；Vs为工区地层实测点的横波速度；ρ为工区地层实测点的密度。
9.如权利要求8所述的流体识别方法，其特征在于，根据完全水饱和网格数据点，确定完全水饱和曲线，包括：
拾取工区地层的岩石物理模板上的完全水饱和数据网格点，将完全水饱和数据网格点采用对数域多项式进行拟合，确定完全水饱和曲线。
10.如权利要求9所述的流体识别方法，其特征在于，按如下方式，确定完全水饱和曲线：
fw(x)＝exp(a·ln4x+b·ln3x+c·ln2x+d·lnx+e)，
其中，fw(x)为完全水饱和曲线；a,b,c,d,e为拟合系数，x为横波阻抗。
11.如权利要求10所述的流体识别方法，其特征在于，根据工区地层的测井数据、岩石物理模板和完全水饱和曲线，确定储层判别因子，包括：
将工区地层的测井数据投影到岩石物理模板上，得到投影点；
根据投影点与完全水饱和曲线的相对位置，利用工区地层实测点的纵波模量、工区地层实测点的剪切模量和岩石物理模板的完全水饱和曲线，确定储层判别因子。
12.如权利要求11所述的流体识别方法，其特征在于，根据投影点与完全水饱和曲线的相对位置，利用工区地层实测点的纵波模量、工区地层实测点的剪切模量和岩石物理模板的完全水饱和曲线，确定储层判别因子，包括：
根据工区地层的测井数据，确定工区地层中的实测点的纵波模量和剪切模量；
按照如下方式，确定储层判别因子：
RF＝fw(Mmea)-Mmea/Gmea，
其中，RF为储层判别因子；fw(x)为完全水饱和曲线；Mmea为工区地层中实测点的纵波模量；Gmea为工区地层中实测点的剪切模量。
13.如权利要求6所述的流体识别方法，其特征在于，根据储层判别因子，判定工区地层中的储层和非储层，包括：
储层判别因子的值大于等于零，则工区地层的中的实测点对应的地层是储层；
储层判别因子的值小于零，则工区地层中的实测点对应的地层是非储层。
14.如权利要求6所述的流体识别方法，其特征在于，
根据工区地层的岩石物理模板，确定工区地层中储层的实测点的孔隙度，包括：
将工区地层中储层的实测点采用最小二乘方法搜索在岩石物理模板上距其最近的网格点，得到工区地层中储层的实测点的孔隙度。
15.如权利要求14所述的流体识别方法，其特征在于，
根据工区地层中储层的实测点的孔隙度和工区地层的物性参数，确定工区地层中储层的实测点的孔隙度校正因子，包括：
按如下方式确定工区地层中储层的实测点的孔隙度校正因子：
其中，
为孔隙度校正因子；Kw为水的体积模量；Kg为气的体积模量；Kma为岩石基质的体积模量；Kdry为岩石骨架的体积模量；Gma为岩石基质的剪切模量；Gdry为岩石物理模板上网格点对应干燥岩石的剪切模量；
为孔隙度。
16.如权利要求1所述的流体识别方法，其特征在于，根据工区地层中储层的实测点的流体因子，确定工区地层中烃类的相对含量，包括：
工区地层中储层的实测点的流体因子的取值介于0到1之间，包含0和1；
工区地层中储层的实测点的流体因子的取值为1时，表征工区地层完全含水；
工区地层中储层的实测点的流体因子的取值为0时，表征工区地层完全含烃；
工区地层中储层的实测点的流体因子的取值越小，表征工区地层中水的相对含量越小，烃类的相对含量越高；
工区地层中储层的实测点的流体因子的取值越大，表征工区地层中水的相对含量越高，烃类的相对含量越小。
17.一种流体识别装置，其特征在于，包括：
数据获取模块，用于获取工区地层的测井数据和物性参数；
岩石物理模板确定模块，用于根据工区地层的物性参数，确定工区地层的岩石物理模板；
孔隙度确定模块，用于根据工区地层的岩石物理模板和测井数据，确定工区地层中储层的实测点的孔隙度；
流体因子确定模块，用于根据工区地层中储层的实测点的孔隙度和工区地层的物性参数，确定工区地层中储层的实测点的流体因子；
烃类相对含量确定模块，用于根据工区地层中储层的实测点的流体因子，确定工区地层中烃类的相对含量；
根据工区地层中储层的实测点的孔隙度和工区地层的物性参数，确定工区地层中储层的实测点的流体因子，包括：
根据工区地层中储层的实测点的孔隙度和工区地层的物性参数，确定工区地层中储层的实测点的孔隙度校正因子；
根据工区地层中储层的实测点的孔隙度校正因子和孔隙度，确定工区地层中储层的实测点的流体因子；
按如下方式，根据工区地层中储层的实测点的孔隙度校正因子和孔隙度，确定工区地层中储层的实测点的流体因子：
其中，NFF为流体因子；
为孔隙度；
和
分别表示与工区地层中储层的实测点孔隙度相同的完全饱水岩石的纵波模量和剪切模量；Mmea为工区地层实测点的纵波模量；Gmea为工区地层实测点的剪切模量。
18.如权利要求17所述的流体识别装置，其特征在于，
岩石物理模板确定模块，用于根据工区地层的物性参数，确定工区地层的岩石物理模板，包括：
根据工区地层的物性参数，确定工区地层的纵波模量和剪切模量；
将工区地层的纵波模量与剪切模量的比，确定为工区地层的纵横波模量比；
根据工区地层的纵波模量和纵横波模量比，确定工区地层的岩石物理模板。
19.如权利要求18所述的流体识别装置，其特征在于，
孔隙度确定模块，用于根据工区地层的岩石物理模板和测井数据，确定工区地层中储层的实测点的孔隙度，包括：
根据工区地层的测井数据和岩石物理模板，确定完全水饱和曲线；
根据工区地层的测井数据、岩石物理模板和完全水饱和曲线，确定储层判别因子；
根据储层判别因子，判定工区地层中的储层和非储层；
根据工区地层的岩石物理模板，确定工区地层中的储层的实测点的孔隙度。
20.如权利要求17所述的流体识别装置，其特征在于，
烃类相对含量确定模块，用于根据工区地层中储层的实测点的流体因子，确定工区地层中烃类的相对含量，包括：
工区地层中储层的实测点的流体因子的取值介于0到1之间，包含0和1；
工区地层中储层的实测点的流体因子的取值为1时，表征工区地层完全含水；
工区地层中储层的实测点的流体因子的取值为0时，表征工区地层完全含烃；
工区地层中储层的实测点的流体因子的取值越小，表征工区地层中水的相对含量越小，烃类的相对含量越高；
工区地层中储层的实测点的流体因子的取值越大，表征工区地层中水的相对含量越高，烃类的相对含量越小。
21.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现权利要求1至16任一项所述流体识别方法。
22.一种计算机可读存储介质，所述计算机可读存储介质存储有执行实现权利要求1至16任一项所述流体识别方法的计算机程序。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910247471.5A CN111751900B (zh) | 2019-03-29 | 2019-03-29 | 一种流体识别方法和装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910247471.5A CN111751900B (zh) | 2019-03-29 | 2019-03-29 | 一种流体识别方法和装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111751900A CN111751900A (zh) | 2020-10-09 |
| CN111751900B true CN111751900B (zh) | 2023-04-25 |
Family
ID=72671693
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910247471.5A Active CN111751900B (zh) | 2019-03-29 | 2019-03-29 | 一种流体识别方法和装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111751900B (zh) |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107797139A (zh) * | 2016-08-31 | 2018-03-13 | 中国石油化工股份有限公司 | 页岩储层游离气含气量地震预测方法及系统 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2810097A4 (en) * | 2012-02-02 | 2015-12-23 | Headwave Inc | SYSTEM, METHOD AND COMPUTER READABLE MEDIUM FOR THE INTERACTIVE IDENTIFICATION OF UNDERGROUND REGIONS |
| CN103760081B (zh) * | 2013-12-31 | 2016-01-06 | 中国石油天然气股份有限公司 | 基于孔隙结构特征的碳酸盐岩储层的气藏预测方法及系统 |
| CN106054248B (zh) * | 2016-07-15 | 2017-07-18 | 河海大学 | 一种基于大面积致密储层地震岩石物理反演方法 |
| EP3494416B1 (en) * | 2016-08-05 | 2022-06-01 | Chevron U.S.A. Inc. | System and method for petro-elastic modeling |
| CN109115987B (zh) * | 2018-07-20 | 2021-01-29 | 中国石油天然气股份有限公司 | 一种基于岩石物理模型的流体因子的评价方法及装置 |
-
2019
- 2019-03-29 CN CN201910247471.5A patent/CN111751900B/zh active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107797139A (zh) * | 2016-08-31 | 2018-03-13 | 中国石油化工股份有限公司 | 页岩储层游离气含气量地震预测方法及系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111751900A (zh) | 2020-10-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109165375B (zh) | 基岩岩性识别和物性计算方法及装置 | |
| Nazemi et al. | The effect of carbonate reservoir heterogeneity on Archie's exponents (a and m), an example from Kangan and Dalan gas formations in the central Persian Gulf | |
| Xiao et al. | A fracture identification method for low-permeability sandstone based on R/S analysis and the finite difference method: A case study from the Chang 6 reservoir in Huaqing oilfield, Ordos Basin | |
| CN103630939B (zh) | 一种气层识别评价方法 | |
| CN105182424A (zh) | 一种基于斑块饱和模型定量预测储层孔隙度的方法和装置 | |
| CN108629459B (zh) | 储层含烃孔隙的检测方法及装置 | |
| CN111381292B (zh) | 一种预测砂岩含烃储层的测井解释方法与装置 | |
| CN110658555A (zh) | 储层流体因子构建方法及储层流体识别方法 | |
| CN110529106B (zh) | 一种利用测井资料确定煤层显微组分含量的方法 | |
| CN110017136B (zh) | 一种基于视水层电阻率的水淹层识别与产水率预测方法 | |
| CN109826623B (zh) | 一种致密砂岩储层层理缝的地球物理测井判识方法 | |
| CN104632206B (zh) | 一种石油测井标准层确定方法及装置 | |
| CN110344822B (zh) | 含钙质致密储层含水饱和度的确定方法 | |
| CN111751900B (zh) | 一种流体识别方法和装置 | |
| CN114428365A (zh) | 一种基于含气特征的致密砂岩气层测井识别方法及装置 | |
| CN108089228B (zh) | 一种确定地层岩石特性的解释数据方法和装置 | |
| CN111239809B (zh) | 一种基于属性融合原理的优势储层确定方法及系统 | |
| Qian et al. | A rock physics driven Bayesian inversion for TOC in the Fuling shale gas reservoir | |
| CN111241652B (zh) | 一种确定地层原油粘度的方法及装置 | |
| CN113267809B (zh) | I类页岩储层预测方法及装置 | |
| CN114060015B (zh) | 一种致密砂岩含气性的评价方法及装置 | |
| CN113138412B (zh) | 深层页岩气孔隙度地震预测方法及装置 | |
| CN116930023A (zh) | 致密砂岩相控分类孔隙度测井精细解释方法和装置 | |
| CN112394392B (zh) | 对烃源岩分布情况进行评价的方法和装置 | |
| CN112799135B (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 |