CN110927787B - 基于虚拟岩性的横波速度计算方法及系统 - Google Patents
基于虚拟岩性的横波速度计算方法及系统 Download PDFInfo
- Publication number
- CN110927787B CN110927787B CN201811095695.0A CN201811095695A CN110927787B CN 110927787 B CN110927787 B CN 110927787B CN 201811095695 A CN201811095695 A CN 201811095695A CN 110927787 B CN110927787 B CN 110927787B
- Authority
- CN
- China
- Prior art keywords
- curve
- initial
- content
- sandstone
- mudstone
- 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
- 238000004364 calculation method Methods 0.000 title claims abstract description 40
- 238000000034 method Methods 0.000 claims abstract description 30
- 238000012937 correction Methods 0.000 claims abstract description 7
- 239000012530 fluid Substances 0.000 claims description 51
- 238000002156 mixing Methods 0.000 claims description 29
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 28
- 239000011435 rock Substances 0.000 claims description 22
- 239000008398 formation water Substances 0.000 claims description 8
- 239000011148 porous material Substances 0.000 description 6
- 230000000694 effects Effects 0.000 description 4
- 238000005056 compaction Methods 0.000 description 3
- 239000004576 sand Substances 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000010998 test method Methods 0.000 description 2
- 239000002023 wood Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000007596 consolidation process Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 210000003000 inclusion body Anatomy 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling Methods 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/303—Analysis for determining velocity profiles or travel times
-
- 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/622—Velocity, density or impedance
- G01V2210/6222—Velocity; travel time
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
公开了一种基于虚拟岩性的横波速度计算方法及系统。该方法可以包括:步骤1:根据砂岩含量曲线和泥岩含量曲线,获得初始岩性模型;步骤2:针对初始岩性模型,通过弹性参数计算公式,获得初始横波曲线、初始纵波曲线和初始密度曲线;步骤3:以实测纵波曲线作为目标曲线,以初始纵波曲线为校正曲线,进行迭代计算,获得虚拟砂岩含量和虚拟泥岩含量;步骤4:将虚拟砂岩含量和虚拟泥岩含量分别替代砂岩含量曲线和泥岩含量曲线,重复步骤1,获得优化岩性模型;步骤5:针对优化岩性模型,通过弹性参数计算公式,获得最终横波曲线。本发明通过虚拟岩性优化建模参数,精确简便地计算横波速度。
Description
技术领域
本发明涉及地球物理学测井岩石物理应用领域,更具体地,涉及一种基于虚拟岩性的横波速度计算方法及系统。
背景技术
横波测井资料主要是用于岩石物理分析、AVO正演、AVO属性分析以及叠前反演等工作,通过纵、横波速度信息的结合可以有效地提高储层预测以及流体检测精度。然而由于偶极子声波测井仪器与常规测井仪器不同,造成额外费用,导致很多井往往没有横波资料或是部分横波测井数据质量一般,不利于叠前反演工作的开展。横波估算主要是为了获得更多以及相对更加准确的横波资料,提高叠前储层预测及油气检测的精度。
一般来说,横波预测的方法可划分为经验公式法、多元拟合法和岩石物理建模法。其中,经验公式法通过拟合纵、横波速度关系完成横波预测,如Greenberg-Castagna纯岩性及混合岩性多项关系式法;多元拟合法实际上是对经验公式的一种拓展,如多元线性拟合法是在拟合的过程中加入密度、伽玛、电阻率等其它测井曲线;而岩石物理建模法通过构建岩石骨架模型及流体参数并将它们输入岩石物理模型计算横波速度。在碎屑岩完全固结、低孔隙、测井解释曲线比较齐全的情况下,常假定岩石由颗粒和球形或椭球型孔隙包含体组成,此类的Xu-White模型是目前估算砂泥岩横波估算最主要的方法。
在科研领域,对碎屑岩岩石物理建模理论及优化讨论的比较多,而在生产领域,利用岩石物理模型往往存在一定的问题,比如实际地层复杂的岩性、孔隙情况常不能简化成等效的纯岩性或两三种简单的岩性组合,此时利用岩石物理建模进行横波估算往往与实测值差异较大;另一方面,对于建模至关重要的总孔隙度曲线通常也很难获取,直接影响着横波估算的精度;此外,如需整段进行建模,有涉及到地层压实效应的影响,无法用一个固定的模型准确表征。
因此,有必要开发一种基于虚拟岩性的横波速度计算方法及系统。
公开于本发明背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解,而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。
发明内容
本发明提出了一种基于虚拟岩性的横波速度计算方法及系统,其能够通过虚拟岩性优化建模参数,精确简便地计算横波速度。
根据本发明的一方面,提出了一种基于虚拟岩性的横波速度计算方法。所述方法可以包括:步骤1:根据砂岩含量曲线和泥岩含量曲线,获得初始岩性模型;步骤2:针对所述初始岩性模型,通过弹性参数计算公式,获得初始横波曲线、初始纵波曲线和初始密度曲线;步骤3:以实测纵波曲线作为目标曲线,以所述初始纵波曲线为校正曲线,进行迭代计算,获得虚拟砂岩含量和虚拟泥岩含量;步骤4:将所述虚拟砂岩含量和所述虚拟泥岩含量分别替代所述砂岩含量曲线和所述泥岩含量曲线,重复步骤1,获得优化岩性模型;步骤5:针对所述优化岩性模型,通过弹性参数计算公式,获得最终横波曲线。
优选地,所述步骤1包括:根据砂岩含量曲线和泥岩含量曲线,获得初始骨架模型;根据含水饱和度曲线和含油饱和度曲线,进行流体混合,获得初始流体模型;确定所述初始骨架模型中所述砂岩的裂缝纵横比和所述泥岩的裂缝纵横比;根据所述初始骨架模型、所述初始流体模型、所述砂岩的裂缝纵横比和所述泥岩的裂缝纵横比,获得所述初始岩性模型。
优选地,根据砂岩含量曲线和泥岩含量曲线,获得初始骨架模型包括:根据岩石物理特征,确定砂岩、泥岩的体积模量、剪切模量、密度,获得砂岩含量曲线和泥岩含量曲线;根据所述砂岩含量曲线和所述泥岩含量曲线,进行砂岩和泥岩的骨架混合,获得初始骨架模型。
优选地,根据含水饱和度曲线和含油饱和度曲线,进行流体混合,获得初始流体模型包括:计算地层水、油的体积模量、剪切模量和密度,获得含水饱和度曲线和含油饱和度曲线;根据所述含水饱和度曲线和所述含油饱和度曲线,进行流体混合,获得初始流体模型。
优选地,所述步骤3包括:步骤301:以实测纵波曲线作为目标曲线,设置初始砂岩泥岩含量比和迭代次数;步骤302:针对所述初始纵波曲线进行Xu-White计算,获得校正纵波曲线;步骤303:利用所述校正纵波曲线替换所述初始纵波曲线,重复步骤302,直至重复次数达到所述迭代次数,根据获得的校正纵波曲线计算虚拟砂岩含量和虚拟泥岩含量。
根据本发明的另一方面,提出了一种基于虚拟岩性的横波速度计算系统,其特征在于,该系统包括:存储器,存储有计算机可执行指令;处理器,所述处理器运行所述存储器中的计算机可执行指令,执行以下步骤:步骤1:根据砂岩含量曲线和泥岩含量曲线,获得初始岩性模型;步骤2:针对所述初始岩性模型,通过弹性参数计算公式,获得初始横波曲线、初始纵波曲线和初始密度曲线;步骤3:以实测纵波曲线作为目标曲线,以所述初始纵波曲线为校正曲线,进行迭代计算,获得虚拟砂岩含量和虚拟泥岩含量;步骤4:将所述虚拟砂岩含量和所述虚拟泥岩含量分别替代所述砂岩含量曲线和所述泥岩含量曲线,重复步骤1,获得优化岩性模型;步骤5:针对所述优化岩性模型,通过弹性参数计算公式,获得最终横波曲线。
优选地,所述步骤1包括:根据砂岩含量曲线和泥岩含量曲线,获得初始骨架模型;根据含水饱和度曲线和含油饱和度曲线,进行流体混合,获得初始流体模型;确定所述初始骨架模型中所述砂岩的裂缝纵横比和所述泥岩的裂缝纵横比;根据所述初始骨架模型、所述初始流体模型、所述砂岩的裂缝纵横比和所述泥岩的裂缝纵横比,获得所述初始岩性模型。
优选地,根据砂岩含量曲线和泥岩含量曲线,获得初始骨架模型包括:根据岩石物理特征,确定砂岩、泥岩的体积模量、剪切模量、密度,获得砂岩含量曲线和泥岩含量曲线;根据所述砂岩含量曲线和所述泥岩含量曲线,进行砂岩和泥岩的骨架混合,获得初始骨架模型。
优选地,根据含水饱和度曲线和含油饱和度曲线,进行流体混合,获得初始流体模型包括:计算地层水、油的体积模量、剪切模量和密度,获得含水饱和度曲线和含油饱和度曲线;根据所述含水饱和度曲线和所述含油饱和度曲线,进行流体混合,获得初始流体模型。
优选地,所述步骤3包括:步骤301:以实测纵波曲线作为目标曲线,设置初始砂岩泥岩含量比和迭代次数;步骤302:针对所述初始纵波曲线进行Xu-White计算,获得校正纵波曲线;步骤303:利用所述校正纵波曲线替换所述初始纵波曲线,重复步骤302,直至重复次数达到所述迭代次数,根据获得的校正纵波曲线计算虚拟砂岩含量和虚拟泥岩含量。
其有益效果在于:
(1)得到的横波曲线更加准确:
(2)解决了复杂碎屑岩地区岩石物理建模的简单等效;
(3)解决了整段岩石物理建模横波估算不考虑压实效应的影响;
(4)与其他优化建模相比,本方法简单,不需要考虑复杂的建模过程,具备很强的实用性。
本发明具有其它的特性和优点,这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的,或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述,这些附图和具体实施方式共同用于解释本发明的特定原理。
附图说明
通过结合附图对本发明示例性实施例进行更详细的描述,本发明的上述以及其它目的、特征和优势将变得更加明显,其中,在本发明示例性实施例中,相同的参考标号通常代表相同部件。
图1示出了根据本发明的基于虚拟岩性的横波速度计算方法的步骤的流程图。
图2示出了根据本发明的一个实施例的初始横波曲线、初始纵波曲线和初始密度曲线的示意图。
图3示出了根据本发明的一个实施例的横波曲线、纵波曲线和密度曲线的示意图。
图4a、图4b分别示出了根据本发明的一个实施例的实测数据与计算的横波结果的弹性参数交会图。
具体实施方式
下面将参照附图更详细地描述本发明。虽然附图中显示了本发明的优选实施例,然而应该理解,可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反,提供这些实施例是为了使本发明更加透彻和完整,并且能够将本发明的范围完整地传达给本领域的技术人员。
图1示出了根据本发明的基于虚拟岩性的横波速度计算方法的步骤的流程图。
在该实施例中,根据本发明的基于虚拟岩性的横波速度计算方法可以包括:步骤1:根据砂岩含量曲线和泥岩含量曲线,获得初始岩性模型;步骤2:针对初始岩性模型,通过弹性参数计算公式,获得初始横波曲线、初始纵波曲线和初始密度曲线;步骤3:以实测纵波曲线作为目标曲线,以初始纵波曲线为校正曲线,进行迭代计算,获得虚拟砂岩含量和虚拟泥岩含量;步骤4:将虚拟砂岩含量和虚拟泥岩含量分别替代砂岩含量曲线和泥岩含量曲线,重复步骤1,获得优化岩性模型;步骤5:针对优化岩性模型,通过弹性参数计算公式,获得最终横波曲线。
在一个示例中,步骤1包括:根据砂岩含量曲线和泥岩含量曲线,获得初始骨架模型;根据含水饱和度曲线和含油饱和度曲线,进行流体混合,获得初始流体模型;确定初始骨架模型中砂岩的裂缝纵横比和泥岩的裂缝纵横比;根据初始骨架模型、初始流体模型、砂岩的裂缝纵横比和泥岩的裂缝纵横比,获得初始岩性模型。
在一个示例中,根据砂岩含量曲线和泥岩含量曲线,获得初始骨架模型包括:根据岩石物理特征,确定砂岩、泥岩的体积模量、剪切模量、密度,获得砂岩含量曲线和泥岩含量曲线;根据砂岩含量曲线和泥岩含量曲线,进行砂岩和泥岩的骨架混合,获得初始骨架模型。
在一个示例中,根据含水饱和度曲线和含油饱和度曲线,进行流体混合,获得初始流体模型包括:计算地层水、油的体积模量、剪切模量和密度,获得含水饱和度曲线和含油饱和度曲线;根据含水饱和度曲线和含油饱和度曲线,进行流体混合,获得初始流体模型。
在一个示例中,步骤3包括:步骤301:以实测纵波曲线作为目标曲线,设置初始砂岩泥岩含量比和迭代次数;步骤302:针对初始纵波曲线进行Xu-White计算,获得校正纵波曲线;步骤303:利用校正纵波曲线替换初始纵波曲线,重复步骤302,直至重复次数达到迭代次数,根据获得的校正纵波曲线计算虚拟砂岩含量和虚拟泥岩含量。
具体地,根据本发明的基于虚拟岩性的横波速度计算方法可以包括:
步骤1:根据岩石物理特征,确定砂岩、泥岩的体积模量、剪切模量、密度,获得砂岩含量曲线和泥岩含量曲线;根据砂岩含量曲线和泥岩含量曲线,采用Voigt-Reuss-Hill加权法进行砂岩和泥岩的骨架混合,获得初始骨架模型,计算其等效体积模量、等效剪切模量、采用等效密度公式计算等效密度;计算地层水、油的体积模量、剪切模量和密度,获得含水饱和度曲线和含油饱和度曲线;根据含水饱和度曲线和含油饱和度曲线,采用Wood公式,进行流体混合,获得初始流体模型;一般砂岩储层为粒间孔,而泥岩存在一定的裂缝,通过迭代试验的方法,通过与实测曲线的对比,确定初始骨架模型中砂岩的裂缝纵横比和泥岩的裂缝纵横比;根据初始骨架模型、初始流体模型、砂岩的裂缝纵横比和泥岩的裂缝纵横比,采用Gassmann、Wyllie和理论相结合的砂泥岩混合Xu-White模型,逐步加入泥岩孔隙和砂岩孔隙,建立初始岩性模型。
步骤2:针对初始岩性模型,通过弹性参数计算公式,获得初始横波曲线、初始纵波曲线和初始密度曲线。
步骤3包括:步骤301:以实测纵波曲线作为目标曲线,设置初始砂岩泥岩含量比和迭代次数;步骤302:针对初始纵波曲线进行Xu-White计算,获得校正纵波曲线;步骤303:利用校正纵波曲线替换初始纵波曲线,重复步骤302,直至重复次数达到迭代次数,根据获得的校正纵波曲线计算虚拟砂岩含量和虚拟泥岩含量。还可以设定方差阈值,计算校正纵波曲线与目标曲线的采样点值方差,若大于方差阈值,则以校正纵波曲线为初始纵波曲线,进行下一次迭代计算,若小于方差阈值,则根据校正纵波曲线计算虚拟砂岩含量和虚拟泥岩含量。
步骤4:将虚拟砂岩含量和虚拟泥岩含量分别替代砂岩含量曲线和泥岩含量曲线,重复步骤1,获得优化岩性模型。
步骤5:针对优化岩性模型,通过弹性参数计算公式,获得最终横波曲线。
本方法通过虚拟岩性优化建模参数,精确简便地计算横波速度。
应用示例
为便于理解本发明实施例的方案及其效果,以下给出一个具体应用示例。本领域技术人员应理解,该示例仅为了便于理解本发明,其任何具体细节并非意在以任何方式限制本发明。
根据本发明的基于虚拟岩性的横波速度计算方法包括:
步骤1:根据岩石物理特征,确定砂岩、泥岩的体积模量、剪切模量、密度,获得砂岩含量曲线和泥岩含量曲线;根据砂岩含量曲线和泥岩含量曲线,采用Voigt-Reuss-Hill加权法进行砂岩和泥岩的骨架混合,获得初始骨架模型,计算其等效体积模量、等效剪切模量、采用等效密度公式计算等效密度;计算地层水、油的体积模量、剪切模量和密度,获得含水饱和度曲线和含油饱和度曲线;根据含水饱和度曲线和含油饱和度曲线,采用Wood公式,进行流体混合,获得初始流体模型;一般砂岩储层为粒间孔,而泥岩存在一定的裂缝,通过迭代试验的方法,通过与实测曲线的对比,确定初始骨架模型中砂岩的裂缝纵横比和泥岩的裂缝纵横比;根据初始骨架模型、初始流体模型、砂岩的裂缝纵横比和泥岩的裂缝纵横比,采用Gassmann、Wyllie和理论相结合的砂泥岩混合Xu-White模型,逐步加入泥岩孔隙和砂岩孔隙,建立初始岩性模型。
图2示出了根据本发明的一个实施例的初始横波曲线、初始纵波曲线和初始密度曲线的示意图。
步骤2:针对初始岩性模型,通过弹性参数计算公式,获得初始横波曲线、初始纵波曲线和初始密度曲线,如图2所示。
步骤3包括:步骤301:以实测纵波曲线作为目标曲线,设置初始砂岩泥岩含量比和迭代次数;步骤302:针对初始纵波曲线进行Xu-White计算,获得校正纵波曲线;步骤303:利用校正纵波曲线替换初始纵波曲线,重复步骤302,直至重复次数达到迭代次数,根据获得的校正纵波曲线计算虚拟砂岩含量和虚拟泥岩含量。
步骤4:将虚拟砂岩含量和虚拟泥岩含量分别替代砂岩含量曲线和泥岩含量曲线,重复步骤1,获得优化岩性模型。
图3示出了根据本发明的一个实施例的横波曲线、纵波曲线和密度曲线的示意图。
图4a、图4b分别示出了根据本发明的一个实施例的实测数据与计算的横波结果的弹性参数交会图。
步骤5:针对优化岩性模型,通过弹性参数计算公式,获得最终横波曲线,如图3所示,实测数据与计算的横波结果的弹性参数交会图分别如图4a、图4b所示。
综上所述,本发明通过虚拟岩性优化建模参数,精确简便地计算横波速度。
本领域技术人员应理解,上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果,并不意在将本发明的实施例限制于所给出的任何示例。
根据本发明的基于虚拟岩性的横波速度计算系统,其特征在于,该系统包括:存储器,存储有计算机可执行指令;处理器,处理器运行存储器中的计算机可执行指令,执行以下步骤:步骤1:根据砂岩含量曲线和泥岩含量曲线,获得初始岩性模型;步骤2:针对初始岩性模型,通过弹性参数计算公式,获得初始横波曲线、初始纵波曲线和初始密度曲线;步骤3:以实测纵波曲线作为目标曲线,以初始纵波曲线为校正曲线,进行迭代计算,获得虚拟砂岩含量和虚拟泥岩含量;步骤4:将虚拟砂岩含量和虚拟泥岩含量分别替代砂岩含量曲线和泥岩含量曲线,重复步骤1,获得优化岩性模型;步骤5:针对优化岩性模型,通过弹性参数计算公式,获得最终横波曲线。
在一个示例中,步骤1包括:根据砂岩含量曲线和泥岩含量曲线,获得初始骨架模型;根据含水饱和度曲线和含油饱和度曲线,进行流体混合,获得初始流体模型;确定初始骨架模型中砂岩的裂缝纵横比和泥岩的裂缝纵横比;根据初始骨架模型、初始流体模型、砂岩的裂缝纵横比和泥岩的裂缝纵横比,获得初始岩性模型。
在一个示例中,根据砂岩含量曲线和泥岩含量曲线,获得初始骨架模型包括:根据岩石物理特征,确定砂岩、泥岩的体积模量、剪切模量、密度,获得砂岩含量曲线和泥岩含量曲线;根据砂岩含量曲线和泥岩含量曲线,进行砂岩和泥岩的骨架混合,获得初始骨架模型。
在一个示例中,根据含水饱和度曲线和含油饱和度曲线,进行流体混合,获得初始流体模型包括:计算地层水、油的体积模量、剪切模量和密度,获得含水饱和度曲线和含油饱和度曲线;根据含水饱和度曲线和含油饱和度曲线,进行流体混合,获得初始流体模型。
在一个示例中,步骤3包括:步骤301:以实测纵波曲线作为目标曲线,设置初始砂岩泥岩含量比和迭代次数;步骤302:针对初始纵波曲线进行Xu-White计算,获得校正纵波曲线;步骤303:利用校正纵波曲线替换初始纵波曲线,重复步骤302,直至重复次数达到迭代次数,根据获得的校正纵波曲线计算虚拟砂岩含量和虚拟泥岩含量。
本系统通过虚拟岩性优化建模参数,精确简便地计算横波速度。
以上已经描述了本发明的各实施例,上述说明是示例性的,并非穷尽性的,并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下,对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。
Claims (6)
1.一种基于虚拟岩性的横波速度计算方法,其特征在于,包括:
步骤1:根据砂岩含量曲线和泥岩含量曲线,获得初始岩性模型;
步骤2:针对所述初始岩性模型,通过弹性参数计算公式,获得初始横波曲线、初始纵波曲线和初始密度曲线;
步骤3:以实测纵波曲线作为目标曲线,以所述初始纵波曲线为校正曲线,进行迭代计算,获得虚拟砂岩含量和虚拟泥岩含量;
步骤4:将所述虚拟砂岩含量和所述虚拟泥岩含量分别替代所述砂岩含量曲线和所述泥岩含量曲线,重复步骤1,获得优化岩性模型;
步骤5:针对所述优化岩性模型,通过弹性参数计算公式,获得最终横波曲线;
其中,所述步骤1包括:
根据砂岩含量曲线和泥岩含量曲线,获得初始骨架模型;
根据含水饱和度曲线和含油饱和度曲线,进行流体混合,获得初始流体模型;
确定所述初始骨架模型中所述砂岩的裂缝纵横比和所述泥岩的裂缝纵横比;
根据所述初始骨架模型、所述初始流体模型、所述砂岩的裂缝纵横比和所述泥岩的裂缝纵横比,获得所述初始岩性模型;
其中,所述步骤3包括:
步骤301:以实测纵波曲线作为目标曲线,设置初始砂岩泥岩含量比和迭代次数;
步骤302:针对所述初始纵波曲线进行Xu-White计算,获得校正纵波曲线;
步骤303:利用所述校正纵波曲线替换所述初始纵波曲线,重复步骤302,直至重复次数达到所述迭代次数,根据获得的校正纵波曲线计算虚拟砂岩含量和虚拟泥岩含量。
2.根据权利要求1所述的基于虚拟岩性的横波速度计算方法,其中,根据砂岩含量曲线和泥岩含量曲线,获得初始骨架模型包括:
根据岩石物理特征,确定砂岩、泥岩的体积模量、剪切模量、密度,获得砂岩含量曲线和泥岩含量曲线;
根据所述砂岩含量曲线和所述泥岩含量曲线,进行砂岩和泥岩的骨架混合,获得初始骨架模型。
3.根据权利要求1所述的基于虚拟岩性的横波速度计算方法,其中,根据含水饱和度曲线和含油饱和度曲线,进行流体混合,获得初始流体模型包括:
计算地层水、油的体积模量、剪切模量和密度,获得含水饱和度曲线和含油饱和度曲线;
根据所述含水饱和度曲线和所述含油饱和度曲线,进行流体混合,获得初始流体模型。
4.一种基于虚拟岩性的横波速度计算系统,其特征在于,该系统包括:
存储器,存储有计算机可执行指令;
处理器,所述处理器运行所述存储器中的计算机可执行指令,执行以下步骤:
步骤1:根据砂岩含量曲线和泥岩含量曲线,获得初始岩性模型;
步骤2:针对所述初始岩性模型,通过弹性参数计算公式,获得初始横波曲线、初始纵波曲线和初始密度曲线;
步骤3:以实测纵波曲线作为目标曲线,以所述初始纵波曲线为校正曲线,进行迭代计算,获得虚拟砂岩含量和虚拟泥岩含量;
步骤4:将所述虚拟砂岩含量和所述虚拟泥岩含量分别替代所述砂岩含量曲线和所述泥岩含量曲线,重复步骤1,获得优化岩性模型;
步骤5:针对所述优化岩性模型,通过弹性参数计算公式,获得最终横波曲线;
其中,所述步骤1包括:
根据砂岩含量曲线和泥岩含量曲线,获得初始骨架模型;
根据含水饱和度曲线和含油饱和度曲线,进行流体混合,获得初始流体模型;
确定所述初始骨架模型中所述砂岩的裂缝纵横比和所述泥岩的裂缝纵横比;
根据所述初始骨架模型、所述初始流体模型、所述砂岩的裂缝纵横比和所述泥岩的裂缝纵横比,获得所述初始岩性模型;
其中,所述步骤3包括:
步骤301:以实测纵波曲线作为目标曲线,设置初始砂岩泥岩含量比和迭代次数;
步骤302:针对所述初始纵波曲线进行Xu-White计算,获得校正纵波曲线;
步骤303:利用所述校正纵波曲线替换所述初始纵波曲线,重复步骤302,直至重复次数达到所述迭代次数,根据获得的校正纵波曲线计算虚拟砂岩含量和虚拟泥岩含量。
5.根据权利要求4所述的基于虚拟岩性的横波速度计算系统,其中,根据砂岩含量曲线和泥岩含量曲线,获得初始骨架模型包括:
根据岩石物理特征,确定砂岩、泥岩的体积模量、剪切模量、密度,获得砂岩含量曲线和泥岩含量曲线;
根据所述砂岩含量曲线和所述泥岩含量曲线,进行砂岩和泥岩的骨架混合,获得初始骨架模型。
6.根据权利要求4所述的基于虚拟岩性的横波速度计算系统,其中,根据含水饱和度曲线和含油饱和度曲线,进行流体混合,获得初始流体模型包括:
计算地层水、油的体积模量、剪切模量和密度,获得含水饱和度曲线和含油饱和度曲线;
根据所述含水饱和度曲线和所述含油饱和度曲线,进行流体混合,获得初始流体模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811095695.0A CN110927787B (zh) | 2018-09-19 | 2018-09-19 | 基于虚拟岩性的横波速度计算方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811095695.0A CN110927787B (zh) | 2018-09-19 | 2018-09-19 | 基于虚拟岩性的横波速度计算方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110927787A CN110927787A (zh) | 2020-03-27 |
CN110927787B true CN110927787B (zh) | 2022-05-13 |
Family
ID=69856063
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811095695.0A Active CN110927787B (zh) | 2018-09-19 | 2018-09-19 | 基于虚拟岩性的横波速度计算方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110927787B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111708985A (zh) * | 2020-04-22 | 2020-09-25 | 中海石油(中国)有限公司 | 一种疏松砂岩地层测井横波速度预测方法 |
CN114114454B (zh) * | 2020-08-31 | 2023-06-23 | 中国石油化工股份有限公司 | 考虑压实作用的岩石物理建模及纵波速度估算方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4293933A (en) * | 1975-03-17 | 1981-10-06 | Schlumberger Technology Corporation | Well logging apparatus and method: synthetic logs and synthetic seismograms with extrapolated reflector dip from log measurements |
CN103954999A (zh) * | 2014-05-06 | 2014-07-30 | 成都理工大学 | 一种适用于低孔隙度砂泥岩地层的横波速度预测新方法 |
CN104181585A (zh) * | 2014-06-27 | 2014-12-03 | 中国石油化工股份有限公司 | 一种地球物理勘探横波估算方法和系统 |
CN104570064A (zh) * | 2013-10-09 | 2015-04-29 | 中国石油化工股份有限公司 | 一种砂岩地层横波速度计算方法 |
CN106814388A (zh) * | 2016-12-27 | 2017-06-09 | 中国石油大学(北京) | 砂泥岩储层地层压力的地震预测方法及装置 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8612156B2 (en) * | 2010-03-05 | 2013-12-17 | Vialogy Llc | Active noise injection computations for improved predictability in oil and gas reservoir discovery and characterization |
CN103527184B (zh) * | 2013-10-28 | 2016-04-20 | 北京大学 | 一种白云岩储层的预测方法和系统 |
-
2018
- 2018-09-19 CN CN201811095695.0A patent/CN110927787B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4293933A (en) * | 1975-03-17 | 1981-10-06 | Schlumberger Technology Corporation | Well logging apparatus and method: synthetic logs and synthetic seismograms with extrapolated reflector dip from log measurements |
CN104570064A (zh) * | 2013-10-09 | 2015-04-29 | 中国石油化工股份有限公司 | 一种砂岩地层横波速度计算方法 |
CN103954999A (zh) * | 2014-05-06 | 2014-07-30 | 成都理工大学 | 一种适用于低孔隙度砂泥岩地层的横波速度预测新方法 |
CN104181585A (zh) * | 2014-06-27 | 2014-12-03 | 中国石油化工股份有限公司 | 一种地球物理勘探横波估算方法和系统 |
CN106814388A (zh) * | 2016-12-27 | 2017-06-09 | 中国石油大学(北京) | 砂泥岩储层地层压力的地震预测方法及装置 |
Non-Patent Citations (5)
Title |
---|
Prestack forward modeling of tight reservoirs based on the Xu–White model;Zhu Chao et al.;《APPLIED GEOPHYSICS》;20150930;第12卷(第3期);第421-431页 * |
Study of the presalt dolomite reservoir prediction in the East Margin of Pre-Caspian Basin;Zhang Jing et al.;《SEG Las Vegas 2012 Annual Meeting》;20121231;第1-5页 * |
基于 Xu -White 模型的优化测井横波速度预测技术研究与应用;乔悦东等;《石油天然气学报 ( 江汉石油学院学报)》;20071031;第29卷(第5期);第100-101页 * |
松辽盆地龙凤山气田致密砂岩含气性预测研究;李瑞磊等;《石油物探》;20171130;第56卷(第6期);第874-881页 * |
虚拟孔隙度优化 Xu-White 模型法预测横波速度;王金伟等;《断块油气田》;20110731;第18卷(第4期);第446-447页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110927787A (zh) | 2020-03-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108072902B (zh) | 一种碳酸盐岩岩石物理模型的建模方法 | |
US8423337B2 (en) | Method for multi-scale geomechanical model analysis by computer simulation | |
CN107203005B (zh) | 一种定量化计算裂缝描述参数的方法 | |
US7676349B2 (en) | Integrated anisotropic rock physics model | |
US8768672B2 (en) | Method for predicting time-lapse seismic timeshifts by computer simulation | |
CN109655903B (zh) | 页岩层横波速度预测方法及系统 | |
CN103954999B (zh) | 一种适用于低孔隙度砂泥岩地层的横波速度预测方法 | |
WO2015103582A1 (en) | Multistage oilfield design optimization under uncertainty | |
CN105445791A (zh) | 一种基于多种地震属性的地层孔隙压力预测方法 | |
CN111399050B (zh) | 高角度裂缝预测方法及装置 | |
US10761231B2 (en) | Generating a high-resolution lithology model for subsurface formation evaluation | |
CN104252007A (zh) | 一种相容性岩石物理建模方法 | |
CN109799540B (zh) | 基于地质信息约束下的火山岩型铀矿床磁化率反演方法 | |
CN112965103B (zh) | 一种多重孔隙储层叠前地震概率化多道反演方法 | |
CN110927787B (zh) | 基于虚拟岩性的横波速度计算方法及系统 | |
CN110969699B (zh) | 一种高精度三维地应力建模方法及其应用 | |
CN110095811B (zh) | 膏岩层速度模型构建处理方法及装置 | |
CN110967742B (zh) | 一种孔隙度反演方法及系统 | |
CN109490988A (zh) | 建立适用于硬岩的岩石物理新模型的方法 | |
CN114428372A (zh) | 一种自适应的岩石物理建模方法 | |
CN112147700A (zh) | 速度异常区的低频模型构建方法及系统 | |
CN113917531B (zh) | 基于扰动体积模量的压力系数预测方法及电子设备 | |
NO20190217A1 (en) | Correcting biases in microseismic-event data | |
US10802177B2 (en) | Evaluating hydrocarbon reserves using tool response models | |
CN111624652B (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 |