CN110703330B - 基于岩石物理的碳酸盐岩测井校正和横波速度预测方法 - Google Patents
基于岩石物理的碳酸盐岩测井校正和横波速度预测方法 Download PDFInfo
- Publication number
- CN110703330B CN110703330B CN201810754627.4A CN201810754627A CN110703330B CN 110703330 B CN110703330 B CN 110703330B CN 201810754627 A CN201810754627 A CN 201810754627A CN 110703330 B CN110703330 B CN 110703330B
- Authority
- CN
- China
- Prior art keywords
- wave velocity
- carbonate
- model
- rock
- curve
- 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
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 title claims abstract description 126
- 239000011435 rock Substances 0.000 title claims abstract description 124
- 238000000034 method Methods 0.000 title claims abstract description 41
- 238000012937 correction Methods 0.000 title claims abstract description 35
- 238000012360 testing method Methods 0.000 claims abstract description 24
- 239000011148 porous material Substances 0.000 claims description 37
- 229910052500 inorganic mineral Inorganic materials 0.000 claims description 17
- 239000011707 mineral Substances 0.000 claims description 17
- 229920006395 saturated elastomer Polymers 0.000 claims description 16
- 239000000470 constituent Substances 0.000 claims description 12
- 239000010453 quartz Substances 0.000 claims description 9
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 claims description 9
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 6
- 239000012530 fluid Substances 0.000 claims description 4
- 208000013201 Stress fracture Diseases 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 9
- 238000011161 development Methods 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 239000004576 sand Substances 0.000 description 2
- 238000005553 drilling Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 150000003839 salts Chemical class 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/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting static or dynamic corrections; Stacking
-
- 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
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
技术领域
本发明涉及一种基于岩石物理的碳酸盐岩测井校正和横波速度预测方法,属于岩石物理研究领域。
背景技术
油气勘探中测井数据的质量好坏对整个油气勘探开发至关重要。然而,由于实际测井数据会受到环境、人为、仪器等多方面影响,因此,测井数据存在着很多不确定因素,容易造成井震标定出现差错,地震解释错误。同时,由于横波在实际钻井中难获得,因此,获取高质量的测井数据至关重要。
现有常见的测井数据校正方法,多是基于经验关系的校正,如Raymer公式(Raymer等1980),或者是基于多种测井数据拟合,采用拟合公式进行其他不合理测井数据的校正。在地震岩石物理研究过程中,发现采用岩石物理模型预测出的测井数据精度较高,符合地震地质规律,应用比较成熟的是Xu-White模型,该模型是基于砂泥岩建立的岩石物理模型,能够对砂泥岩的孔隙特征、矿物成分、流体等进行模拟。Xu and Payne(2009)基于Xu-White模型建立的碳酸盐岩岩石物理模型,该模型能够描述碳酸盐岩复杂的孔隙空间,将岩石中的孔隙考虑为三种孔隙,一种是圆孔,一种是微裂隙,另一种是介于前两者之间的孔隙。Luanxiao Zhao等(2013)年在Xu and Payne模型的基础上进行改进建立了新的碳酸盐岩模型。这些岩石物理模型能够在已知矿物成分、孔隙度、孔隙结构的情况下对岩石的速度等参数进行估算,在测井中常常用来估算横波速度值。
但是,目前已有的测井数据校正方法和横波速度预测一种是基于经验关系的,一种是基于复杂岩石物理模型来估算的。由于基于经验关系的测井校正方法经验关系的区域局限性导致经验模型的不实用性,而纯粹的岩石物理模型方法则依赖于已知测井的矿物成分和孔隙度,测井的矿物成分和孔隙度又是由其他测井数据反演出来的,因此,在一定程度上保留了测井数据存在的问题。上述两种方法获得的结果都会导致预测测井曲线与实际地质信息不匹配,进一步影响井震标定和储层预测。因此,需要构建合适的基于岩石物理的碳酸盐岩测井校正方法和横波预测方法,为碳酸盐岩储层预测和开发提供支撑作用。
发明内容
针对上述问题,本发明的目的是提供一种基于岩石物理的碳酸盐岩测井校正和横波速度预测方法,其预测结果精度高,能够解决测井原始数据岩石物理规律不明确或错误的现象,以及测井横波速度缺失的问题,能够为井震标定和储层预测提供高质量的测井数据。
为实现上述目的,本发明采取以下技术方案:一种基于岩石物理的碳酸盐岩测井校正方法,其特征在于,建立碳酸盐岩密度回归关系模型,获取测井数据校正后的密度曲线;利用岩石物理模型重构模型密度曲线、纵波速度曲线和横波速度曲线,确定合格的模型纵波速度曲线和模型横波速度曲线;根据岩心测试数据密度和纵波速度的拟合关系,获取拟合的纵波速度曲线,根据岩心测试数据纵波速度和横波速度的拟合关系,获取拟合的横波速度曲线;根据拟合的纵波速度曲线和合格的模型纵波速度曲线进行分析,确定测井数据校正后的纵波速度曲线,根据拟合的横波速度曲线和合格的模型横波速度曲线进行分析,确定测井数据校正后的横波速度曲线。
在一个具体的实施例中,碳酸盐岩孔隙结构包括孔隙纵横比为0.01的微裂隙、孔隙纵横比为0.1的粒间孔和孔隙纵横比为0.9的硬孔隙。
在一个具体的实施例中,利用碳酸盐岩测井数据中非扩径储层段的深侧向电阻率、补偿声波、中子测井,建立碳酸盐岩密度回归关系模型。
在一个具体的实施例中,建立碳酸盐岩密度回归关系模型所依据的公式为:
ρ=aRd+bCNL+cAC+d (1)
式中,ρ为碳酸盐岩储层密度,Rd为碳酸盐岩测井数据中非扩径储层段的深侧向电阻率,CNL为碳酸盐岩测井数据中非扩径储层段的中子,AC为碳酸盐岩测井数据中非扩径储层段的补偿声波,a、b、c、d为拟合系数。
在一个具体的实施例中,利用Xu-Payne模型重构模型密度、纵波速度和横波速度曲线,调整Xu-Payne模型中矿物纵横比和孔隙纵横比参数使得重构的密度曲线与测井数据校正后的密度曲线差值达到极小值,确定合格的模型纵波速度曲线和模型横波速度曲线。
在一个具体的实施例中,在利用Xu-Payne模型重构模型密度、纵波速度和横波速度曲线的过程中,利用Voigt-Reuss-Hill平均模型计算碳酸盐岩饱和多种矿物骨架的岩石等效弹性模型,利用微分等效介质模型计算碳酸盐岩包含孔隙的等效岩石弹性模量,利用Gassmann方程计算碳酸盐岩孔隙饱和流体时的岩石弹性模量。
在一个具体的实施例中,利用Voigt-Reuss-Hill平均模型计算碳酸盐岩包含石英和泥质骨架的岩石等效弹性模型所依据的公式为:
式中,Mm为碳酸盐岩包含石英和泥质骨架的岩石等效弹性模量,Mi为碳酸盐岩组成岩石矿物成分第i个组成成分的模量,fi为碳酸盐岩组成岩石矿物成分第i个组成成分的体积分量,Mv为采用Voigt上限方法计算获得的岩石模量,MR为采用Reuss下限方法计算获得的岩石模量;
利用微分等效介质模型计算碳酸盐岩包含孔隙的岩石弹性模量所依据的公式为:
K*(0)=K1 (7)
μ*(0)=μ1 (8)
式中,K1为碳酸盐岩初始主相相1的体积模量,μ1为碳酸盐岩初始主相相1的剪切模量,K2为碳酸盐岩初始主相相2的体积模量,μ2为碳酸盐岩初始主相相2的剪切模量,y为碳酸盐岩初始主相相2的体积含量,P为几何因子,Q为几何因子;
利用Gassmann方程计算碳酸盐岩饱和水时的岩石弹性模量所依据的公式为:
μsat=μdry (10)
式中,Kdry为碳酸盐岩干岩石骨架的有效体积模量,Ksat为碳酸盐岩饱和岩石的有效体积模量,Kf为碳酸盐岩饱和的水的有效体积模量,φ为孔隙度,μsat为碳酸盐岩饱和岩石的有效剪切模量,μdry为碳酸盐岩干岩石骨架的有效剪切模量。
在一个具体的实施例中,确定实际岩心测试数据纵波速度和横波速度的拟合关系所依据的公式为:
VS=0.4017VP+944.68 (11)
式中,VS为横波速度,VP为纵波速度。
在一个具体的实施例中,确定实际岩心测试数据密度和纵波速度的拟合关系所依据的公式为:
VP=2005.8ρ+826.89 (12)
一种基于岩石物理的碳酸盐岩横波速度预测方法,其特征在于,包括建立碳酸盐岩密度回归关系模型,获取测井数据校正后的密度曲线;利用岩石物理模型重构密度曲线、纵波速度曲线和横波速度曲线,确定合格的模型纵波速度曲线和模型横波速度曲线;根据岩心测试数据密度和纵波速度的拟合关系,获取拟合的纵波速度曲线,根据岩心测试数据纵波速度和横波速度的拟合关系,获取拟合的横波速度曲线;根据合格的模型横波速度与拟合的横波速度值交会分析,获取校正后的碳酸盐岩模型横波速度。
本发明由于采取以上技术方案,其具有以下优点:本发明建立了碳酸盐岩密度回归关系模型,能够消除扩径、环境因素对密度曲线的影响,利用Xu-Payne模型重构密度曲线、纵波速度曲线和横波速度曲线,能够简捷、有效地确定合格的模型纵波速度曲线和横波速度曲线,结合实际岩心测试数据得到了区域拟合的纵波速度曲线和横波速度曲线,提高测井数据校正后的纵波速度曲线和横波速度曲线精度,从而能够解决测井数据岩石物理规律不明确或错误的现象,以及测井横波速度缺失的问题,碳酸盐岩模型横波速度能够较为准确地预测碳酸盐岩储层的横波速度,且更加符合岩石物理规律,本发明碳酸盐岩测井校正和横波速度预测方法能够为井震标定和储层预测提供高质量的测井数据。
本发明的其它特征和优点将在随后的说明书中阐述,并且,部分的从说明书中变得显而易见,或者通过实施本发明而了解。本发明的目的和其他优点可通过在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。
附图说明
为了更清楚的说明本发明实施例中的技术方案,下面将对实施例描述中所需要的附图做简单的介绍:
图1是本发明的一个具体实施例获取测井数据校正后的密度曲线的流程结构示意图;
图2是本发明的一个具体实施例纵波速度和横波速度校正的流程结构示意图;
图3是本发明的一个具体实施例纵波速度校正和横波速度拟合的结构示意图;
图4是本发明的一个具体实施例纵波速度校正和密度拟合的结构示意图;
图5是本发明的一个具体实施例模型重构密度曲线和测井数据校正后的密度曲线的结构示意图;
图6是本发明的一个具体实施例测井数据校正后纵波速度曲线和测井数据原始纵波速度曲线的结构示意图;
图7是本发明的一个具体实施例测井数据校正后横波速度曲线和测井数据原始横波速度曲线的结构示意图;
图8是本发明的一个具体实施例测井数据孔隙度和测井数据原始横波速度交会的结构示意图;
图9是本发明的一个具体实施例测井数据孔隙度和校正后模型横波速度交会的结构示意图。
具体实施方式
以下将结合附图及实施例来详细说明本发明的实施方式,借此对本发明如何应用技术手段来解决技术问题,并达成技术效果的实现过程能充分理解并据以实施。需要说明的是,只要不构成冲突,本发明中的各个实施例以及各实施例中的各个特征可以相互结合,所形成的技术方案均在本发明的保护范围之内。
本发明提出的基于岩石物理的碳酸盐岩测井校正,包括
1)建立碳酸盐岩密度回归关系模型,获取碳酸盐储层测井数据校正后的密度曲线
如图1所示,首先,在一套碳酸盐岩测井数据中优选出一层段的非扩径优质测井数据,接着,利用该测井数据中的深侧向电阻率测井、补偿声波测井、中子测井建立碳酸盐岩密度的回归关系模型,然后,利用该碳酸盐岩密度回归模型去对其它层段的密度曲线进行重新估算,以消除扩径、环境因素对密度曲线的影响,从而获得整个碳酸盐储层测井数据校正后的密度曲线。利用碳酸盐岩密度回归关系模型获得的测井数据校正后的密度曲线能够消除扩径和环境造成的测井数据原始密度曲线局部的失真。
建立碳酸盐岩密度回归关系模型所依据的公式为:
ρ=aRd+bCNL+cAC+d (1)
式中,ρ为碳酸盐岩储层密度,Rd为碳酸盐岩测井数据中非扩径储层段的深侧向电阻率,CNL为碳酸盐岩测井数据中非扩径储层段的中子,AC为碳酸盐岩测井数据中非扩径储层段的补偿声波,a、b、c、d为拟合系数。
在一个优选的实施例中,a=6.227×10-7,b=-0.20537,c=-0.0024,d=3.001。
在一个优选的实施例中,碳酸盐岩孔隙结构包括孔隙纵横比为0.01的微裂隙、孔隙纵横比为0.1的粒间孔和孔隙纵横比为0.9的硬孔隙。
2)利用岩石物理模型重构模型密度曲线、纵波速度曲线和横波速度曲线,通过模型密度曲线与测井数据校正后的密度曲线分析,确定合格的模型纵波速度曲线和模型横波速度曲线
如图2所示,利用Xu-Payne模型重构模型密度、纵波速度和横波速度曲线,调整Xu-Payne模型中矿物纵横比和孔隙纵横比参数,使得重构的密度曲线与测井数据校正后的密度曲线差值达到极小值(即误差最小),确定合格(合乎目标)的模型纵波速度曲线和模型横波速度曲线,此时,模型纵波速度曲线和横波速度曲线的可靠性好。
在利用Xu-Payne模型重构模型密度、纵波速度和横波速度曲线的过程中,利用Voigt-Reuss-Hill平均模型计算碳酸盐岩包含多种矿物骨架的岩石等效弹性模型,利用微分等效介质模型计算碳酸盐岩包含孔隙的等效岩石弹性模量,利用Gassmann方程计算碳酸盐岩孔隙饱和流体时的岩石弹性模量。
利用Voigt-Reuss-Hill平均模型计算碳酸盐岩包含石英和泥质骨架的岩石等效弹性模型所依据的公式为:
式中,Mm为碳酸盐岩包含石英和泥质骨架的岩石等效弹性模量,Mi为碳酸盐岩组成岩石矿物成分第i个组成成分的模量,fi为碳酸盐岩组成岩石矿物成分第i个组成成分的体积分量,Mv为采用Voigt上限方法计算获得的岩石模量,MR为采用Reuss下限方法计算获得的岩石模量。
利用微分等效介质模型计算碳酸盐岩包含孔隙的岩石弹性模量所依据的公式为:
K*(0)=K1 (7)
μ*(0)=μ1 (8)
式中,K1为碳酸盐岩初始主相相1的体积模量,μ1为碳酸盐岩初始主相相1的剪切模量,K2为碳酸盐岩初始主相相2的体积模量,μ2为碳酸盐岩初始主相相2的剪切模量,y为碳酸盐岩初始主相相2的体积含量,P为几何因子,Q为几何因子。
其中,碳酸盐岩的初始主相相1可以为石英或泥质。若碳酸盐岩的初始主相相1为石英时,则其初始主相相2为泥质。若碳酸盐岩的初始主相相1为泥质时,则其初始主相相2为石英。
利用Gassmann方程计算碳酸盐岩饱和水时的岩石弹性模量所依据的公式为:
μsat=μdry (10)
式中,Kdry为碳酸盐岩干岩石骨架的有效体积模量,Ksat为碳酸盐岩饱和岩石的有效体积模量,Kf为碳酸盐岩饱和的水的有效体积模量,φ为孔隙度,μsat为碳酸盐岩饱和岩石的有效剪切模量,μdry为碳酸盐岩干岩石骨架的有效剪切模量。
3)根据实际岩心测试数据密度和纵波速度的拟合关系,获取拟合的纵波速度曲线,根据实际岩心测试数据纵波速度和横波速度的拟合关系,获取拟合的横波速度曲线
确定实际岩心测试数据纵波速度和横波速度的拟合关系所依据的公式为:
VS=0.4017VP+944.68 (11)
式中,VS为横波速度,VP为纵波速度。
确定实际岩心测试数据密度和纵波速度的拟合关系所依据的公式为:
VP=2005.8ρ+826.89 (12)。
接着,根据实际岩心测试数据密度和纵波速度的拟合关系,获取拟合的纵波速度曲线,根据实际岩心测试数据纵波速度和横波速度的拟合关系,获取拟合的横波速度曲线。
4)根据拟合的纵波速度曲线和合格的模型纵波速度曲线进行分析,确定测井数据校正后的纵波速度曲线,根据拟合的横波速度曲线和合格的模型横波速度曲线进行分析,确定测井数据校正后的横波速度曲线
根据拟合的纵波速度曲线和合格的模型纵波速度曲线进行分析,对测井数据原始纵波速度曲线进行校正,确定测井数据校正后的纵波速度曲线。具体过程:利用重构模型建立的纵波速度曲线与拟合的纵波速度曲线和测井数据原始纵波速度曲线对比分析,以拟合的纵波速度曲线为基准,将测井数据原始纵波速度曲线中不合理数据采用模型纵波速度曲线替换,从而确定测井数据校正后的纵波速度曲线。同理,根据拟合的横波速度曲线和合格的模型横波速度曲线进行分析,对测井数据原始纵横速度曲线进行校正,确定测井数据校正后的横波速度曲线。
本发明提出的基于岩石物理的碳酸盐岩横波速度预测方法,包括
1)建立碳酸盐岩密度回归关系模型,获取测井数据校正后的密度曲线;
2)利用岩石物理模型重构密度曲线、纵波速度曲线和横波速度曲线,确定合格的模型纵波速度曲线和横波速度曲线;
3)根据实际岩心测试数据密度和纵波速度的拟合关系,获取拟合的纵波速度曲线,根据实际岩心测试数据纵波速度和横波速度的拟合关系,获取拟合的横波速度曲线;
4)根据合格的模型横波速度与拟合的横波速度值交会分析,获取校正后的碳酸盐岩模型横波速度。
在一个优选的实施例中,基于岩石物理的碳酸盐岩横波速度预测方法前三步的流程与基于岩石物理的碳酸盐岩测井校正前三步的流程相同。在此基础上,根据模型横波速度与拟合的横波速度值交会分析,将模型横波速度曲线中的不合理数据采用拟合的横波速度值替换,获取校正后的碳酸盐岩模型横波速度曲线。
下面列举一具体的实施例
如图3、图4所示,某地区实际岩石物理实验室测试交会图。对该地区的测井数据进行了校正。如图5所示,利用密度回归关系模型校正的测井数据(灰色虚线)和岩石物理模型重构的密度曲线(黑色实线),从图5中可以看出模型重构的密度曲线与校正后的密度曲线基本是一致的,除了个别层段数据的尖锐。如图6、图7所示,测井数据校正后的纵波速度曲线和测井数据校正后的横波速度曲线与测井数据原始纵波速度曲线和测井数据原始横波速度曲线的对比图,从图7中可以看出,估算的横波速度曲线与实测的横波速度曲线符合较好。如图8、图9所示,对该区域测井数据孔隙度和测井数据原始横波速度交会分析以及测井数据孔隙度和校正后模型横波速度交会分析。
虽然本发明所公开的实施方式如上,但所述的内容只是为了便于理解本发明而采用的实施方式,并非用以限定本发明。任何本发明所属技术领域内的技术人员,在不脱离本发明所公开的精神和范围的前提下,可以在实施的形式上及细节上作任何的修改与变化,但本发明的专利保护范围,仍须以所附的权利要求书所界定的范围为准。
Claims (9)
1.一种基于岩石物理的碳酸盐岩测井校正方法,其特征在于,包括
建立碳酸盐岩密度回归关系模型,获取测井数据校正后的密度曲线;
利用岩石物理模型重构模型密度曲线、纵波速度曲线和横波速度曲线,确定合格的模型纵波速度曲线和模型横波速度曲线;
根据岩心测试数据密度和纵波速度的拟合关系,获取拟合的纵波速度曲线,根据岩心测试数据纵波速度和横波速度的拟合关系,获取拟合的横波速度曲线;
根据拟合的纵波速度曲线和合格的模型纵波速度曲线进行分析,确定测井数据校正后的纵波速度曲线,根据拟合的横波速度曲线和合格的模型横波速度曲线进行分析,确定测井数据校正后的横波速度曲线;
其中,利用Xu-Payne模型重构模型密度、纵波速度和横波速度曲线,调整Xu-Payne模型中矿物纵横比和孔隙纵横比参数使得重构的密度曲线与测井数据校正后的密度曲线差值达到极小值,确定合格的模型纵波速度曲线和模型横波速度曲线。
2.根据权利要求1所述的基于岩石物理的碳酸盐岩测井校正方法,其特征在于,碳酸盐岩孔隙结构包括孔隙纵横比为0.01的微裂隙、孔隙纵横比为0.1的粒间孔和孔隙纵横比为0.9的硬孔隙。
3.根据权利要求1或2所述的基于岩石物理的碳酸盐岩测井校正方法,其特征在于,利用碳酸盐岩测井数据中非扩径储层段的深侧向电阻率、补偿声波、中子测井,建立碳酸盐岩密度回归关系模型。
4.根据权利要求3所述的基于岩石物理的碳酸盐岩测井校正方法,其特征在于,建立碳酸盐岩密度回归关系模型所依据的公式为:
ρ=aRd+bCNL+cAC+d (1)
式中,ρ为碳酸盐岩储层密度,Rd为碳酸盐岩测井数据中非扩径储层段的深侧向电阻率,CNL为碳酸盐岩测井数据中非扩径储层段的中子,AC为碳酸盐岩测井数据中非扩径储层段的补偿声波,a、b、c、d为拟合系数。
5.根据权利要求1所述的基于岩石物理的碳酸盐岩测井校正方法,其特征在于,在利用Xu-Payne模型重构模型密度、纵波速度和横波速度曲线的过程中,利用Voigt-Reuss-Hill平均模型计算碳酸盐岩包含多种矿物骨架的岩石等效弹性模型,利用微分等效介质模型计算碳酸盐岩包含孔隙的等效岩石弹性模量,利用Gassmann方程计算碳酸盐岩孔隙饱和流体时的岩石弹性模量。
6.根据权利要求5所述的基于岩石物理的碳酸盐岩测井校正方法,其特征在于,利用Voigt-Reuss-Hill平均模型计算碳酸盐岩包含石英和泥质骨架的岩石等效弹性模型所依据的公式为:
式中,Mm为碳酸盐岩包含石英和泥质骨架的岩石等效弹性模量,Mi为碳酸盐岩组成岩石矿物成分第i个组成成分的模量,fi为碳酸盐岩组成岩石矿物成分第i个组成成分的体积分量,Mv为采用Voigt上限方法计算获得的岩石模量,MR为采用Reuss下限方法计算获得的岩石模量;
利用微分等效介质模型计算碳酸盐岩包含孔隙的岩石弹性模量所依据的公式为:
K*(0)=K1 (7)
μ*(0)=μ1 (8)
式中,K1为碳酸盐岩初始主相相1的体积模量,μ1为碳酸盐岩初始主相相1的剪切模量,K2为碳酸盐岩初始主相相2的体积模量,μ2为碳酸盐岩初始主相相2的剪切模量,y为碳酸盐岩初始主相相2的体积含量,P为几何因子,Q为几何因子;
利用Gassmann方程计算碳酸盐岩饱和水时的岩石弹性模量所依据的公式为:
μsat=μdry (10)
式中,Kdry为碳酸盐岩干岩石骨架的有效体积模量,Ksat为碳酸盐岩饱和岩石的有效体积模量,Kf为碳酸盐岩饱和的水的有效体积模量,φ为孔隙度,μsat为碳酸盐岩饱和岩石的有效剪切模量,μdry为碳酸盐岩干岩石骨架的有效剪切模量。
7.根据权利要求6所述的基于岩石物理的碳酸盐岩测井校正方法,其特征在于,确定岩心测试数据纵波速度和横波速度的拟合关系所依据的公式为:
VS=0.4017VP+944.68 (11)
式中,VS为横波速度,VP为纵波速度。
8.根据权利要求7所述的基于岩石物理的碳酸盐岩测井校正方法,其特征在于,确定岩心测试数据密度和纵波速度的拟合关系所依据的公式为:
VP=2005.8ρ+826.89 (12)。
9.一种基于岩石物理的碳酸盐岩横波速度预测方法,其特征在于,包括
建立碳酸盐岩密度回归关系模型,获取测井数据校正后的密度曲线;
利用岩石物理模型重构模型密度曲线、纵波速度曲线和横波速度曲线,确定合格的模型纵波速度曲线和模型横波速度曲线;
根据岩心测试数据密度和纵波速度的拟合关系,获取拟合的纵波速度曲线,根据岩心测试数据纵波速度和横波速度的拟合关系,获取拟合的横波速度曲线;
根据合格的模型横波速度与拟合的横波速度值交会分析,获取校正后的碳酸盐岩模型横波速度;
其中,利用Xu-Payne模型重构模型密度、纵波速度和横波速度曲线,调整Xu-Payne模型中矿物纵横比和孔隙纵横比参数使得重构的密度曲线与测井数据校正后的密度曲线差值达到极小值,确定合格的模型纵波速度曲线和模型横波速度曲线。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810754627.4A CN110703330B (zh) | 2018-07-10 | 2018-07-10 | 基于岩石物理的碳酸盐岩测井校正和横波速度预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810754627.4A CN110703330B (zh) | 2018-07-10 | 2018-07-10 | 基于岩石物理的碳酸盐岩测井校正和横波速度预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110703330A CN110703330A (zh) | 2020-01-17 |
CN110703330B true CN110703330B (zh) | 2021-12-24 |
Family
ID=69192411
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810754627.4A Active CN110703330B (zh) | 2018-07-10 | 2018-07-10 | 基于岩石物理的碳酸盐岩测井校正和横波速度预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110703330B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111751878B (zh) * | 2020-05-21 | 2023-05-30 | 中国石油天然气股份有限公司 | 横波速度的预测方法和装置 |
CN112182840A (zh) * | 2020-08-28 | 2021-01-05 | 中国石油天然气集团有限公司 | 一种天然气水合物地层衰减建模方法与装置 |
CN112505755B (zh) * | 2020-10-29 | 2021-11-16 | 中国石油集团工程咨询有限责任公司 | 一种基于自适应混合骨架参数的高精度岩石物理建模方法 |
CN114578025B (zh) * | 2022-03-01 | 2023-10-17 | 江苏师范大学 | 煤矸石声速的水岩耦合模型的建立方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103278850A (zh) * | 2013-05-24 | 2013-09-04 | 西安石油大学 | 基于煤岩工业组分物理体积模型的横波时差曲线构建方法 |
CN104950331A (zh) * | 2015-06-30 | 2015-09-30 | 中国石油天然气集团公司 | 一种砂泥岩储层的孔隙度与泥质含量的地震预测方法 |
CN105370272A (zh) * | 2015-12-07 | 2016-03-02 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 页岩气储层的测井评价方法 |
CN106054248A (zh) * | 2016-07-15 | 2016-10-26 | 河海大学 | 一种基于大面积致密储层地震岩石物理反演方法 |
CN106199712A (zh) * | 2016-06-29 | 2016-12-07 | 中国石油天然气集团公司 | 一种确定压裂套管变形区域的方法及装置 |
WO2017058043A1 (en) * | 2015-09-28 | 2017-04-06 | Baker Hughes Incorporated | Correcting shale volume and measuring anisotropy in invaded zone |
CN106597544A (zh) * | 2016-11-25 | 2017-04-26 | 中国石油天然气股份有限公司 | 致密油气藏脆性的预测方法及装置 |
-
2018
- 2018-07-10 CN CN201810754627.4A patent/CN110703330B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103278850A (zh) * | 2013-05-24 | 2013-09-04 | 西安石油大学 | 基于煤岩工业组分物理体积模型的横波时差曲线构建方法 |
CN104950331A (zh) * | 2015-06-30 | 2015-09-30 | 中国石油天然气集团公司 | 一种砂泥岩储层的孔隙度与泥质含量的地震预测方法 |
WO2017058043A1 (en) * | 2015-09-28 | 2017-04-06 | Baker Hughes Incorporated | Correcting shale volume and measuring anisotropy in invaded zone |
CN105370272A (zh) * | 2015-12-07 | 2016-03-02 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 页岩气储层的测井评价方法 |
CN106199712A (zh) * | 2016-06-29 | 2016-12-07 | 中国石油天然气集团公司 | 一种确定压裂套管变形区域的方法及装置 |
CN106054248A (zh) * | 2016-07-15 | 2016-10-26 | 河海大学 | 一种基于大面积致密储层地震岩石物理反演方法 |
CN106597544A (zh) * | 2016-11-25 | 2017-04-26 | 中国石油天然气股份有限公司 | 致密油气藏脆性的预测方法及装置 |
Non-Patent Citations (2)
Title |
---|
基于岩石物理模型的测井约束横波速度计算方法研究;周中彪;《物探化探计算技术》;20100930;第32卷(第5期);第536-541页 * |
基于岩石物理测试模型的测井曲线校正方法;刘爱疆 等;《中国测试》;20130930;第39卷(第5期);第20-23页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110703330A (zh) | 2020-01-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110703330B (zh) | 基于岩石物理的碳酸盐岩测井校正和横波速度预测方法 | |
CN106405651B (zh) | 一种基于测井匹配的全波形反演初始速度模型构建方法 | |
CN103576195B (zh) | 一种随压力变化的裂隙介质横波速度预测方法 | |
CN112133377B (zh) | 一种海域天然气水合物赋存状态的判别方法及系统 | |
CN110133720B (zh) | 一种基于统计岩石物理模型的横波速度预测方法 | |
CN109490963B (zh) | 裂缝储层岩石物理建模方法及系统 | |
CN112255688B (zh) | 一种基于岩石物理理论的三维地震反演地层压力的方法 | |
CN102608664A (zh) | 深度维核磁共振反演获取横向弛豫时间谱的方法及装置 | |
CN104656162A (zh) | 一种富铀烃源岩铀含量及有机碳含量确定方法 | |
CN107387068B (zh) | 一种用于确定页岩气储层游离气含量的方法及系统 | |
WO2020215170A1 (zh) | 基于地震岩石物理实验分析的测井与地震速度匹配方法 | |
CN112505772A (zh) | 一种利用孔隙、裂隙介质弹性波理论反演岩石孔隙分布特征的方法 | |
CN109061752B (zh) | 一种含灰质地层的电阻率曲线校正方法 | |
CN105093331B (zh) | 获取岩石基质体积模量的方法 | |
CN102913240B (zh) | 一种储层流体识别方法 | |
CN110542928A (zh) | 基于vti各向异性传播矩阵的地震响应模拟方法 | |
CN111650646B (zh) | 一种用于地震处理中的声波曲线校正方法 | |
CN205297569U (zh) | 用于确定致密砂岩饱和度指数的装置 | |
CN110716235B (zh) | 一种砂泥岩测井孔隙结构反演方法 | |
CN110007348A (zh) | 一种灰质背景浊积岩储层的岩石物理建模方法 | |
CN117471540A (zh) | 基于叠后地震数据预测储层物性参数的方法、装置和设备 | |
CN113552624B (zh) | 孔隙度预测方法及装置 | |
CN106842329A (zh) | 储层含油饱和度的获取方法与装置 | |
CN110909311A (zh) | 一种计算薄煤层含气量的方法 | |
CN111594156A (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 |