CN107367763A - 确定岩层Biot系数的方法及装置 - Google Patents
确定岩层Biot系数的方法及装置 Download PDFInfo
- Publication number
- CN107367763A CN107367763A CN201710545466.3A CN201710545466A CN107367763A CN 107367763 A CN107367763 A CN 107367763A CN 201710545466 A CN201710545466 A CN 201710545466A CN 107367763 A CN107367763 A CN 107367763A
- Authority
- CN
- China
- Prior art keywords
- formation
- rock
- modulus
- determining
- fluid
- 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
- 239000011435 rock Substances 0.000 title claims abstract description 282
- 238000000034 method Methods 0.000 title claims abstract description 59
- 239000012530 fluid Substances 0.000 claims abstract description 151
- 229920006395 saturated elastomer Polymers 0.000 claims abstract description 96
- 229910052500 inorganic mineral Inorganic materials 0.000 claims abstract description 90
- 239000011707 mineral Substances 0.000 claims abstract description 90
- 238000005755 formation reaction Methods 0.000 claims description 212
- 230000015572 biosynthetic process Effects 0.000 claims description 171
- 238000004590 computer program Methods 0.000 claims description 11
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 6
- 238000010008 shearing Methods 0.000 abstract 2
- 239000011148 porous material Substances 0.000 description 11
- 238000010586 diagram Methods 0.000 description 9
- 230000006870 function Effects 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 238000011545 laboratory measurement Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 239000004568 cement Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
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/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
- G01V1/50—Analysing data
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/616—Data from specific type of measurement
- G01V2210/6169—Data from specific type of measurement using well-logging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
- G01V2210/6242—Elastic parameters, e.g. Young, Lamé or Poisson
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本申请实施例提供了一种确定岩层Biot系数的方法及装置，该包括：根据岩层的测井曲线，确定所述岩层的饱和岩石P波模量、饱和岩石剪切模量、流体体积模量和流体密度；根据所述流体体积模量和所述流体密度确定所述岩层的矿物P波模量；根据所述饱和岩石P波模量、所述饱和岩石剪切模量、所述流体体积模量和所述矿物P波模量，确定所述岩层的Biot系数。本申请实施例的通用性更强，且得到的Biot系数的准确度更高。
Description
技术领域
本申请涉及地震勘探中的岩石物理学技术领域，尤其是涉及一种确定岩层Biot系数的方法及装置。
背景技术
Biot孔隙弹性系数(简称Biot系数)是指在相同孔隙压力下，孔隙体积变化与总体积变化之比。在地震波传播过程中，Biot系数实际反映的是孔隙空间对岩层整体性质的贡献，它是孔隙弹性特征中最重要的参数之一。因此，研究Biot系数具有重要意义。
研究表明，Biot系数与岩层组分、胶结物、岩石孔隙度、孔隙结构、连通通道几何形状、孔隙流体类型及有效压力等因素有关。目前在对岩层Biot系数的研究中，大部分研究是通过特定岩性的实验室测量和基本理论模型来实现的，获得的经验公式仅适用于某个特定区域的某种特定岩性，很难在别的区域针对不同岩性进行推广应用，从而使Biot系数的应用受到限性。
发明内容
本申请实施例的目的在于提供一种通用性更强的确定岩层Biot系数的方法及装置。
为达到上述目的，一方面，本申请实施例提供了一种确定岩层Biot系数的方法，包括：
根据岩层的测井曲线，确定所述岩层的饱和岩石P波模量、饱和岩石剪切模量、流体体积模量和流体密度；
根据所述流体体积模量和所述流体密度确定所述岩层的矿物P波模量；
根据所述饱和岩石P波模量、所述饱和岩石剪切模量、所述流体体积模量和所述矿物P波模量，确定所述岩层的Biot系数。
本申请实施例的确定岩层Biot系数的方法，所述测井曲线包括岩层的深度测井曲线、纵波波速测井曲线、横波波速测井曲线、密度测井曲线、岩石孔隙度测井曲线和含水饱和度测井曲线。
本申请实施例的确定岩层Biot系数的方法所述根据岩层的测井曲线确定所述岩层的饱和岩石P波模量，包括：
根据公式确定所述岩层的饱和岩石P波模量；
其中，Ms为岩层的饱和岩石P波模量，ρs为岩层的饱和岩石密度，Vps为岩层的饱和岩石纵波波速。
本申请实施例的确定岩层Biot系数的方法所述根据岩层的测井曲线确定所述岩层的饱和岩石剪切模量，包括：
根据公式确定所述岩层的饱和岩石剪切模量；
其中，μs为岩层的饱和岩石剪切模量，ρs为岩层的饱和岩石密度，Vss为岩层的饱和岩石横波波速。
本申请实施例的确定岩层Biot系数的方法所述根据岩层的测井曲线确定所述岩层的流体体积模量，包括：
根据公式确定所述岩层的流体体积模量；
其中，Kf为岩层的流体体积模量，Ki为岩层中第i种流体的体积模量，Pi为岩层中第i种流体的体积含量，N为岩层中所包含的流体种数。
本申请实施例的确定岩层Biot系数的方法所述根据岩层的测井曲线确定所述岩层的流体密度，包括：
根据公式确定所述岩层的流体密度；
其中，ρf为岩层的流体密度，Pi为岩层中第i种流体的体积含量，ρi为岩层中第i种流体的密度，N为岩层中所包含的流体种数。
本申请实施例的确定岩层Biot系数的方法所述根据所述流体体积模量和所述流体密度确定所述岩层的矿物P波模量，包括：
根据所述流体体积模量和所述流体密度，确定所述岩层的流体纵波波速；并根据所述流体密度确定所述岩层的矿物密度；
根据所述流体纵波波速确定所述岩层的矿物纵波波速；
根据所述矿物纵波波速和所述矿物密度，确定所述岩层的矿物P波模量。
本申请实施例的确定岩层Biot系数的方法所述根据所述流体体积模量和所述流体密度，确定所述岩层的流体纵波波速，包括：
根据公式确定所述岩层的流体纵波波速；
其中，Vf为岩层的流体纵波波速，Kf为岩层的流体体积模量，ρf为岩层的流体密度。
本申请实施例的确定岩层Biot系数的方法所述根据所述流体密度确定所述岩层的矿物密度，包括：
根据公式ρs＝(1-φ)ρm+φρf确定所述岩层的矿物密度；
其中，ρs为岩层的饱和岩石密度，φ为岩层的岩石孔隙度，ρm为岩层的矿物密度，ρf为岩层的流体密度。
本申请实施例的确定岩层Biot系数的方法所述根据所述流体纵波波速确定所述岩层的矿物纵波波速，包括：
根据公式确定所述岩层的矿物密度；
其中，Vps为岩层的饱和岩石纵波波速，φ为岩层的岩石孔隙度，Vpm为岩层的矿物纵波波速，Vf为岩层的流体纵波波速。
本申请实施例的确定岩层Biot系数的方法所述根据所述矿物纵波波速和所述矿物密度，确定所述岩层的矿物P波模量，包括：
根据公式确定所述岩层的矿物P波模量；
其中，Mm岩层的矿物P波模量，ρm为岩层的矿物密度，Vpm为岩层的矿物纵波波速。
本申请实施例的确定岩层Biot系数的方法所述根据所述饱和岩石P波模量、所述饱和岩石剪切模量、所述流体体积模量和所述矿物P波模量，确定所述岩层的Biot系数，包括：
根据公式aβ2+bβ+c＝0确定所述岩层的Biot系数；
其中，β为岩层的Biot系数，a＝A+B，b＝AC-D，c＝AD， Ms岩层的饱和岩石P波模量，Mm岩层的矿物P波模量，μd为岩层的岩石骨架剪切模量，μd由岩层的饱和岩石剪切模量替换，Kf为岩层的流体体积模量，φ为岩层的岩石孔隙度。
另一方面，本申请实施例还提供了一种确定岩层Biot系数的装置，包括：
第一参数获取模块，用于根据岩层的测井曲线，确定所述岩层的饱和岩石P波模量、饱和岩石剪切模量、流体体积模量和流体密度；
第二参数获取模块，用于根据所述流体体积模量和所述流体密度确定所述岩层的矿物P波模量；
Biot系数确定模块，用于根据所述饱和岩石P波模量、所述饱和岩石剪切模量、所述流体体积模量和所述矿物P波模量，确定所述岩层的Biot系数。
再一方面，本申请实施例还提供了一种计算机存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现以下步骤：
根据岩层的测井曲线，确定所述岩层的饱和岩石P波模量、饱和岩石剪切模量、流体体积模量和流体密度；
根据所述流体体积模量和所述流体密度确定所述岩层的矿物P波模量；
根据所述饱和岩石P波模量、所述饱和岩石剪切模量、所述流体体积模量和所述矿物P波模量，确定所述岩层的Biot系数。
由以上本申请实施例提供的技术方案可见，本申请实施例首先通过岩层的测井曲线，确定岩层的饱和岩石P波模量、饱和岩石剪切模量、流体体积模量和流体密度；然后根据流体体积模量和所述流体密度确定所述岩层的矿物P波模量；最后根据饱和岩石P波模量、饱和岩石剪切模量、流体体积模量和矿物P波模量，确定岩层的Biot系数。由于本申请实施例综合考虑了饱和岩石弹性参数、孔隙流体参数和矿物参数等多种因数对Biot系数的影响，从而提高了估算Biot系数的准确度，而且本申请实施例没有岩性限制，因此其可用于任何岩性的Biot系数估算，从而提高了Biot系数估算的通用性。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。在附图中：
图1为本申请一实施例中确定岩层Biot系数的方法流程图；
图2a～图2e为本申请一实施例中输入的测井曲线；
图3a～图3b为本申请一实施例中获得的饱和岩石弹性参数；
图4a～图4c为本申请一实施例中获得的孔隙流体参数；
图5a～图5c为本申请一实施例中获得的矿物参数；
图6为本申请一实施例中获得的Biot系数；
图7为本申请一实施例中确定岩层Biot系数的装置结构框图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护的范围。
参考图1所示，本申请实施例的确定岩层Biot系数的方法可以包括以下步骤：
S101、根据岩层的测井曲线，确定所述岩层的饱和岩石P波模量、饱和岩石剪切模量、流体体积模量和流体密度。
本申请一实施例中，所述测井曲线可以包括岩层的深度测井曲线、纵波(即P波)波速测井曲线、横波(即S波)波速测井曲线、密度测井曲线、岩石孔隙度测井曲线和含水饱和度测井曲线等。目前这些测井资料比较容易获取，同时本申请实施例无需实验室测量数据，因此具有实现方便和成本低廉的特点。
本申请一实施例中，可根据公式确定所述岩层的饱和岩石P波模量。其中，Ms为岩层的饱和岩石P波模量，ρs为岩层的饱和岩石密度(可通过密度测井曲线和含水饱和度测井曲线获得)，Vps为岩层的饱和岩石纵波波速(可从P波测井曲线中获得)。
本申请一实施例中，可根据公式确定所述岩层的饱和岩石剪切模量。其中，μs为岩层的饱和岩石剪切模量，ρs为岩层的饱和岩石密度，Vss为岩层的饱和岩石横波波速(可从S波测井曲线中获得)。
本申请一实施例中，可根据公式确定所述岩层的流体体积模量。其中，Kf为岩层的流体体积模量，Ki为岩层中第i种流体的体积模量，Pi为岩层中第i种流体的体积含量，N为岩层中所包含的流体种数。
本申请一实施例中，可根据公式确定所述岩层的流体密度。其中，ρf为岩层的流体密度，Pi为岩层中第i种流体的体积含量，ρi为岩层中第i种流体的密度，N为岩层中所包含的流体种数。
S102、根据所述流体体积模量和所述流体密度确定所述岩层的矿物P波模量。
本申请一实施例中，所述根据所述流体体积模量和所述流体密度确定所述岩层的矿物P波模量可以包括如下步骤：
首先，根据所述流体体积模量和所述流体密度，确定所述岩层的流体纵波波速；并根据所述流体密度确定所述岩层的矿物密度。
本申请一实施例中，可根据公式确定所述岩层的流体纵波波速。其中，Vf为岩层的流体纵波波速，Kf为岩层的流体体积模量，ρf为岩层的流体密度。并可根据公式ρs＝(1-φ)ρm+φρf确定所述岩层的矿物密度。其中，ρs为岩层的饱和岩石密度，φ为岩层的岩石孔隙度，ρm为岩层的矿物密度，ρf为岩层的流体密度。
然后，根据所述流体纵波波速确定所述岩层的矿物纵波波速。
本申请一实施例中，可根据公式确定所述岩层的矿物密度；
其中，Vps为岩层的饱和岩石纵波波速，φ为岩层的岩石孔隙度，Vpm为岩层的矿物纵波波速，Vf为岩层的流体纵波波速。
最后，根据所述矿物纵波波速和所述矿物密度，确定所述岩层的矿物P波模量。
本申请一实施例中，可根据公式确定所述岩层的矿物P波模量；
其中，Mm岩层的矿物P波模量，ρm为岩层的矿物密度，Vpm为岩层的矿物纵波波速。
S103、根据所述饱和岩石P波模量、所述饱和岩石剪切模量、所述流体体积模量和所述矿物P波模量，确定所述岩层的Biot系数。
本申请实施例中，可将Gassmann方程表示为Biot系数形式：
Ks＝Kd+β2M
其中，
上式中，Ks、Kd、Km和Kf分别为岩层的饱和岩石体积模量、岩石骨架体积模量、矿物体积模量和流体体积模量，φ为岩层的岩石孔隙度，β为Biot系数。
同样地，岩层的饱和岩石P波体积模量Ms、岩石骨架P波体积模量Md和Biot系数β之间具有如下关系：
Ms＝Md+β2M
将上式中的岩石骨架P波体积模量Ms用矿物P波体积模量Mm表示，即Md＝Mm(1-β)，则有：Ms＝Mm(1-β)+β2M
将带入上式，可得到关于β的一元二次方程：aβ2+bβ+c＝0，其中，β为岩层的Biot系数，a＝A+B，b＝AC-D，c＝AD， Ms岩层的饱和岩石P波模量，Mm岩层的矿物P波模量，μd为岩层的岩石骨架剪切模量，μd由岩层的饱和岩石剪切模量替换，Kf为岩层的流体体积模量，φ为岩层的岩石孔隙度。
本申请一实施例中，可根据公式aβ2+bβ+c＝0确定所述岩层的Biot系数；
其中，β为岩层的Biot系数，a＝A+B，b＝AC-D，c＝AD， Ms岩层的饱和岩石P波模量，Mm岩层的矿物P波模量，μd为岩层的岩石骨架剪切模量，μd由岩层的饱和岩石剪切模量替换，Kf为岩层的流体体积模量，φ为岩层的岩石孔隙度。
本申请实施例综合考虑了饱和岩石、岩石骨架、矿物颗粒、孔隙度及孔隙流体等多种因数对Biot系数的影响，从而提高了估算Biot系数的准确度，而且本申请实施例所涉及的岩石物理理论及模型没有岩性限制，因此其可用于任何岩性的Biot系数估算，因此，本申请实施例的方案具有较强的通用性。
为便于理解本申请，下面介绍本申请的一个示例性实施例。
在该个示例性实施例中，选择一口全波列测井资料(其测井曲线可如图2a～图2e所示)来验证本申请的有效性，该口井的测井资料包括：深度测井曲线(m)、P波波速测井曲线(m/s)、S波波速测井曲线(m/s)、密度(g/cm3)测井曲线、岩石孔隙度测井曲线及含水饱和度测井曲线。
饱和岩石弹性参数估算：以图2a～图2e所示的测井曲线作为输入，通过饱和岩石密度、P波波速及S波波速与弹性模量之间的关系估算饱和岩石P波体积模量和饱和岩石剪切模量(如图3a～图3b所示)。
饱孔隙流体参数估算：通过输入的含水饱和度等流体体积含量参数，分别估算流体体积模量、流体密度和流体P波波速(如图4a～图4c所示)。
矿物参数估算：根据上述得到的饱和岩石密度、岩石孔隙度和流体密度，估算矿物密度(如图5b所示)；根据上述得到的饱和岩石P波波速、岩石孔隙度和流体P波波速，估算矿物P波波速(如图5a所示)；根据上述估算得到的矿物密度和矿物P波波速，估算矿物P波体积模量(如图5c所示)。
Biot系数估算：将Biot-Gassmann方程表示成关于Biot系数的一元二次方程(即aβ2+bβ+c＝0)，在得到饱和岩石P波模量、饱和岩石剪切模量、流体体积模量和矿物P波模量的条件下，通过求解该一元二次方程看得到Biot系数(如图6所示)。
参考图7所示，本申请实施例确定岩层Biot系数的装置可以包括：
第一参数获取模块71，可以用于根据岩层的测井曲线，确定所述岩层的饱和岩石P波模量、饱和岩石剪切模量、流体体积模量和流体密度；
第二参数获取模块72，可以用于根据所述流体体积模量和所述流体密度确定所述岩层的矿物P波模量；
Biot系数确定模块73，可以用于根据所述饱和岩石P波模量、所述饱和岩石剪切模量、所述流体体积模量和所述矿物P波模量，确定所述岩层的Biot系数。
本申请实施例的装置与上述实施例的方法对应，因此，有关于本申请的装置细节，请参见上述实施例的方法，在此不再赘述。
为了描述的方便，描述以上装置时以功能分为各种单元分别描述。当然，在实施本申请时可以把各单元的功能在同一个或多个软件和/或硬件中实现。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中，计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。
内存可能包括计算机可读介质中的非永久性存储器，随机存取存储器(RAM)和/或非易失性内存等形式，如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括暂存电脑可读媒体(transitory media)，如调制的数据信号和载波。
还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、商品或者设备中还存在另外的相同要素。
本领域技术人员应明白，本申请的实施例可提供为方法、系统或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
以上所述仅为本申请的实施例而已，并不用于限制本申请。对于本领域技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本申请的权利要求范围之内。
Claims (14)
1.一种确定岩层Biot系数的方法，其特征在于，包括：
根据岩层的测井曲线，确定所述岩层的饱和岩石P波模量、饱和岩石剪切模量、流体体积模量和流体密度；
根据所述流体体积模量和所述流体密度确定所述岩层的矿物P波模量；
根据所述饱和岩石P波模量、所述饱和岩石剪切模量、所述流体体积模量和所述矿物P波模量，确定所述岩层的Biot系数。
2.如权利要求1所述的确定岩层Biot系数的方法，其特征在于，所述测井曲线包括岩层的深度测井曲线、纵波波速测井曲线、横波波速测井曲线、密度测井曲线、岩石孔隙度测井曲线和含水饱和度测井曲线。
3.如权利要求2所述的确定岩层Biot系数的方法，其特征在于，所述根据岩层的测井曲线确定所述岩层的饱和岩石P波模量，包括：
根据公式确定所述岩层的饱和岩石P波模量；
其中，Ms为岩层的饱和岩石P波模量，ρs为岩层的饱和岩石密度，Vps为岩层的饱和岩石纵波波速。
4.如权利要求2所述的确定岩层Biot系数的方法，其特征在于，所述根据岩层的测井曲线确定所述岩层的饱和岩石剪切模量，包括：
根据公式确定所述岩层的饱和岩石剪切模量；
其中，μs为岩层的饱和岩石剪切模量，ρs为岩层的饱和岩石密度，Vss为岩层的饱和岩石横波波速。
5.如权利要求2所述的确定岩层Biot系数的方法，其特征在于，所述根据岩层的测井曲线确定所述岩层的流体体积模量，包括：
根据公式确定所述岩层的流体体积模量；
其中，Kf为岩层的流体体积模量，Ki为岩层中第i种流体的体积模量，Pi为岩层中第i种流体的体积含量，N为岩层中所包含的流体种数。
6.如权利要求2所述的确定岩层Biot系数的方法，其特征在于，所述根据岩层的测井曲线确定所述岩层的流体密度，包括：
根据公式确定所述岩层的流体密度；
其中，ρf为岩层的流体密度，Pi为岩层中第i种流体的体积含量，ρi为岩层中第i种流体的密度，N为岩层中所包含的流体种数。
7.如权利要求6所述的确定岩层Biot系数的方法，其特征在于，所述根据所述流体体积模量和所述流体密度确定所述岩层的矿物P波模量，包括：
根据所述流体体积模量和所述流体密度，确定所述岩层的流体纵波波速；并根据所述流体密度确定所述岩层的矿物密度；
根据所述流体纵波波速确定所述岩层的矿物纵波波速；
根据所述矿物纵波波速和所述矿物密度，确定所述岩层的矿物P波模量。
8.如权利要求7所述的确定岩层Biot系数的方法，其特征在于，所述根据所述流体体积模量和所述流体密度，确定所述岩层的流体纵波波速，包括：
根据公式确定所述岩层的流体纵波波速；
其中，Vf为岩层的流体纵波波速，Kf为岩层的流体体积模量，ρf为岩层的流体密度。
9.如权利要求7所述的确定岩层Biot系数的方法，其特征在于，所述根据所述流体密度确定所述岩层的矿物密度，包括：
根据公式ρs＝(1-φ)ρm+φρf确定所述岩层的矿物密度；
其中，ρs为岩层的饱和岩石密度，φ为岩层的岩石孔隙度，ρm为岩层的矿物密度，ρf为岩层的流体密度。
10.如权利要求7所述的确定岩层Biot系数的方法，其特征在于，所述根据所述流体纵波波速确定所述岩层的矿物纵波波速，包括：
根据公式确定所述岩层的矿物密度；
其中，Vps为岩层的饱和岩石纵波波速，φ为岩层的岩石孔隙度，Vpm为岩层的矿物纵波波速，Vf为岩层的流体纵波波速。
11.如权利要求7所述的确定岩层Biot系数的方法，其特征在于，所述根据所述矿物纵波波速和所述矿物密度，确定所述岩层的矿物P波模量，包括：
根据公式确定所述岩层的矿物P波模量；
其中，Mm岩层的矿物P波模量，ρm为岩层的矿物密度，Vpm为岩层的矿物纵波波速。
12.如权利要求2所述的确定岩层Biot系数的方法，其特征在于，所述根据所述饱和岩石P波模量、所述饱和岩石剪切模量、所述流体体积模量和所述矿物P波模量，确定所述岩层的Biot系数，包括：
根据公式aβ2+bβ+c＝0确定所述岩层的Biot系数；
其中，β为岩层的Biot系数，a＝A+B，b＝AC-D，c＝AD， Ms岩层的饱和岩石P波模量，Mm岩层的矿物P波模量，μd为岩层的岩石骨架剪切模量，μd由岩层的饱和岩石剪切模量替换，Kf为岩层的流体体积模量，φ为岩层的岩石孔隙度。
13.一种确定岩层Biot系数的装置，其特征在于，包括：
第一参数获取模块，用于根据岩层的测井曲线，确定所述岩层的饱和岩石P波模量、饱和岩石剪切模量、流体体积模量和流体密度；
第二参数获取模块，用于根据所述流体体积模量和所述流体密度确定所述岩层的矿物P波模量；
Biot系数确定模块，用于根据所述饱和岩石P波模量、所述饱和岩石剪切模量、所述流体体积模量和所述矿物P波模量，确定所述岩层的Biot系数。
14.一种计算机存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现以下步骤：
根据岩层的测井曲线，确定所述岩层的饱和岩石P波模量、饱和岩石剪切模量、流体体积模量和流体密度；
根据所述流体体积模量和所述流体密度确定所述岩层的矿物P波模量；
根据所述饱和岩石P波模量、所述饱和岩石剪切模量、所述流体体积模量和所述矿物P波模量，确定所述岩层的Biot系数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710545466.3A CN107367763B (zh) | 2017-07-06 | 2017-07-06 | 确定岩层Biot系数的方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710545466.3A CN107367763B (zh) | 2017-07-06 | 2017-07-06 | 确定岩层Biot系数的方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107367763A true CN107367763A (zh) | 2017-11-21 |
| CN107367763B CN107367763B (zh) | 2019-06-11 |
Family
ID=60305969
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710545466.3A Active CN107367763B (zh) | 2017-07-06 | 2017-07-06 | 确定岩层Biot系数的方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107367763B (zh) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110219643A (zh) * | 2019-06-05 | 2019-09-10 | 中国石油天然气集团有限公司 | 油藏含水饱和度确定方法及装置 |
| CN111368245A (zh) * | 2020-03-12 | 2020-07-03 | 成都理工大学 | 一种考虑多影响因素的定量计算岩石Biot系数方法 |
| CN113281825A (zh) * | 2020-02-20 | 2021-08-20 | 中国石油天然气集团有限公司 | 岩石物理模型构建方法及装置 |
| CN114783531A (zh) * | 2022-04-02 | 2022-07-22 | 广州海洋地质调查局 | 计算砂岩型水合物中胶结型和孔隙充填型饱和度的方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5815465A (en) * | 1997-04-11 | 1998-09-29 | Turgut; Altan | Method and apparatus of classifying marine sediment |
| CN106680878A (zh) * | 2016-12-29 | 2017-05-17 | 中国石油大学(华东) | 一种基于改进Biot系数分析的横波速度估算方法 |
-
2017
- 2017-07-06 CN CN201710545466.3A patent/CN107367763B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5815465A (en) * | 1997-04-11 | 1998-09-29 | Turgut; Altan | Method and apparatus of classifying marine sediment |
| CN106680878A (zh) * | 2016-12-29 | 2017-05-17 | 中国石油大学(华东) | 一种基于改进Biot系数分析的横波速度估算方法 |
Non-Patent Citations (1)
| Title |
|---|
| 尹帅 等: ""海陆过渡相致密砂岩储层Biot系数自适应预测方法研究"", 《石油物探》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110219643A (zh) * | 2019-06-05 | 2019-09-10 | 中国石油天然气集团有限公司 | 油藏含水饱和度确定方法及装置 |
| CN113281825A (zh) * | 2020-02-20 | 2021-08-20 | 中国石油天然气集团有限公司 | 岩石物理模型构建方法及装置 |
| CN111368245A (zh) * | 2020-03-12 | 2020-07-03 | 成都理工大学 | 一种考虑多影响因素的定量计算岩石Biot系数方法 |
| CN111368245B (zh) * | 2020-03-12 | 2023-04-14 | 成都理工大学 | 一种考虑多影响因素的定量计算岩石Biot系数方法 |
| CN114783531A (zh) * | 2022-04-02 | 2022-07-22 | 广州海洋地质调查局 | 计算砂岩型水合物中胶结型和孔隙充填型饱和度的方法 |
| CN114783531B (zh) * | 2022-04-02 | 2022-11-22 | 广州海洋地质调查局 | 计算砂岩型水合物中胶结型和孔隙充填型饱和度的方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107367763B (zh) | 2019-06-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107367763B (zh) | 确定岩层Biot系数的方法及装置 | |
| CN110133725B (zh) | 地震岩石横波速度预测方法及装置 | |
| CN107942404B (zh) | 一种确定裂缝的方法及装置 | |
| CN106959469A (zh) | 地震波的速度及衰减模拟分析方法及装置 | |
| CN107367762A (zh) | 一种确定储层参数的方法及装置 | |
| CN106842304B (zh) | 一种叠前深度偏移方法及装置 | |
| CN111123354A (zh) | 基于频变反射振幅衰减预测致密气层的方法及设备 | |
| CN111208565B (zh) | 基于kt模型的孔缝参数反演方法、装置及存储介质 | |
| CN107422377B (zh) | 利用地震数据预测烃源岩总有机碳含量的方法及装置 | |
| CN110320575B (zh) | 基于岩石物理模型的页岩有机质含量确定方法及装置 | |
| CN109597964B (zh) | 理论油藏渗流参数值确定方法、系统、设备、可读介质 | |
| CN112346130B (zh) | 一种富有机质岩石横波速度预测方法、存储介质及系统 | |
| CN106960264B (zh) | 低渗透砂砾岩地层压力预测方法及装置 | |
| CN114910950B (zh) | 天然气水合物储层下方游离气饱和度的估算方法及系统 | |
| CN111208566B (zh) | 基于sca模型的孔缝参数反演方法、装置及存储介质 | |
| CN111323814B (zh) | 基于岩石物理模板定量确定砂岩储层含水饱和度的方法 | |
| CN111812716A (zh) | 页岩气储层总有机碳含量叠前定量预测方法、装置及设备 | |
| CN111435174B (zh) | 强反射地区地震资料振幅补偿方法及装置 | |
| US20090043549A1 (en) | Methods, apparatus, and products for seismic ray tracing | |
| CN108318925B (zh) | 一种岩溶型碳酸盐岩储层横波预测方法及装置 | |
| CN112817059B (zh) | 峰峰组内幕储层确定方法及装置 | |
| CN106970418B (zh) | 获取地震合成记录的方法及装置 | |
| Mandal et al. | Disturbance of SH-type waves due to moving stress discontinuity in an anisotropic soil layer overlying an inhomogeneous elastic half-space | |
| CN110219643A (zh) | 油藏含水饱和度确定方法及装置 | |
| CN111562629A (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 |