CN106597544A - 致密油气藏脆性的预测方法及装置 - Google Patents

Abstract

本发明公开了一种致密油气藏脆性的预测方法及装置，该方法包括以下步骤：基于获取的测井曲线得到模型关系，利用模型关系再建立横波速度预测模型；基于横波速度预测模型预测得到横波速度曲线；根据横波速度曲线、密度曲线和纵波速度曲线计算弹性参数曲线；基于石英含量测井曲线和叠后地震数据体得到石英含量数据体；基于石英含量数据体得到脆性指数；基于叠前地震资料、纵波速度曲线、密度曲线得到时深曲线和地震子波；基于时深曲线、密度曲线和弹性参数曲线建立初始模型，再得到不同角度的合成角度道集，通过合成角度道集与叠前角道集建立反演方程进而反演得到杨氏模量的数据体；基于石英含量数据体和杨氏模量的数据体得到脆性指数的数据体。

Description

致密油气藏脆性的预测方法及装置

技术领域

本发明涉及石油地球物理勘探技术领域，特别涉及一种致密油气藏脆性的预测方法及装置。

背景技术

近年来非常规油气资源的勘探开发备受国内外专家的关注，致密油气储层作为其重要组成部分，成为目前油气勘探热点和难点，致密油气是指赋存于渗透率极低的致密砂岩、致密碳酸盐岩等储层的油气，单井产量未能形成工业油气流，可在一定的开采措施下形成工业油流。该类储层具有渗透率低、孔隙结构复杂、矿物组分复杂、脆性特征明显等特点。其中，脆性特征是石油地质学家衡量有效致密储层的重要指标，因此，致密储层的脆性预测引起了广大研究者的追捧。目前，储层脆性特征预测的地球物理手段主要包括以下几个方面：

1)基于岩石物理测量和岩石物理建模的致密储层的脆性预测方法。实验测量的方法模拟地下的地质环境，在高温高压的条件下，通过测量与脆性指数密切的相关的岩石特征，如岩石的破碎特征、应力应变关系、岩石的矿物组分等，然后利用这些测量的岩石信息来计算岩石的脆性特征。Kahraman(2015)选取了不同沉积环境下火成岩和变质岩的24个不同的破碎岩石样本进行测量，筛选500g样本制作成大小为19.0-9.53mm的样本在110度的高温条件下加热48小时，在实验测量岩石的破碎指数，通过测量发现岩石的脆性指数和破碎指数存在一个简单的线性关系，从而可以实现岩石脆性特征的预测。Guo(2012)、Liu(2015)充分考虑岩石的矿物组分、孔隙度和孔隙构造，建立针对特殊类型的岩石物理模型，得到岩石的弹性特征与复杂成分和特殊微裂缝属性的关系，通过岩石物理量板来定量描述岩石的脆性特征。

2)以叠前地震资料为基础，结合测井资料和地质信息的叠前反演进行致密储层预测。近年来发展了基于叠前地震反演的致密储层的脆性预测，通过结合地质信息和测井资料，利用叠前地震反演估算纵、横波速度及密度，然后利用这些弹性参数计算储层的杨氏模量和泊松比实现储层的脆性预测。郭旭光(2015)提出了利用岩石力学实验测量的结果，建立以杨氏模量和泊松比为主要参数的脆性评价系统，在此基础上利用叠前同步反演获得纵、横波阻抗和密度，然后利用这些参数计算杨氏模量和泊松比，最后，结合脆性指数剖面和平面图实现致密油储层的储层预测和井位优选。Song(2015)认为随着储层深度的增加，岩石的力学性质由脆性逐渐向塑性变化，充分利用岩石脆性和塑性转换点的特征，通过叠后和叠前反演的技术进行脆性矿物和岩石组分的预测，最终可以为致密油气的勘探和开发提供参考。

虽然上述两种方法在致密油气储层发挥了重要作用，但是对于精细的致密储层的脆性预测具有明显的不足，主要存在以下几方面的缺点：未能综合考虑岩石的力学性质和矿物组分，不能建立正确的脆性指数与储层弹性参数之间的关系，不能定量地描述岩石特征；通过简单叠前反演获取的弹性参数计算的脆性指数，缺乏定量关系的指导，使得反演精度大大降低，预测误差会明显增大；针对特殊的致密储层，震岩石物理、储层地球物理特征、反射模式和AVO特征等方面存在研究基础薄弱，不能有效地指导致密油气的勘探。

发明内容

为了克服现有技术的上述缺陷，本发明实施例中提供了一种致密油气藏脆性的预测方法及装置，其能够提高反演的精度，提高致密油气藏脆性预测的精确性。

本发明实施例的具体技术方案是：

一种致密油气藏脆性的预测方法，包括以下步骤：

基于获取的测井曲线得到横波速度与纵波速度曲线、密度曲线、GR曲线、电阻率曲线之间的模型关系，利用所述模型关系再建立横波速度预测模型；

基于所述横波速度预测模型预测得到横波速度曲线；

根据所述横波速度曲线、所述密度曲线和所述纵波速度曲线计算弹性参数曲线；

基于石英含量测井曲线和叠后地震数据体得到石英含量数据体；

基于所述石英含量数据体得到脆性指数；

基于叠前地震资料、纵波速度曲线、密度曲线得到时深曲线和地震子波；

基于所述时深曲线、所述密度曲线和所述弹性参数曲线建立初始模型，根据所述初始模型与地震子波褶积得到不同角度的合成角度道集，通过所述合成角度道集与叠前角道集建立反演方程进而反演得到杨氏模量的数据体；

基于所述石英含量数据体和所述杨氏模量的数据体得到脆性指数的数据体。

优选地，在所述基于叠前地震资料、纵波速度曲线、密度曲线得到时深曲线的步骤中具体包括：

基于叠前地震资料得到叠后地震资料；

基于所述叠后地震资料、纵波速度曲线、密度曲线得到时深曲线。

优选地，在所述基于获取的测井曲线得到横波速度与纵波速度曲线、密度曲线、GR曲线、电阻率曲线之间的模型关系，利用所述模型关系再建立横波速度预测模型的步骤中，具体为：对密度曲线、GR曲线、电阻率曲线进行环境校正和标准化处理，根据获取的测井曲线得到横波速度，构建所述横波速度和纵波速度曲线、密度曲线、GR曲线、电阻率曲线之间的模型关系，利用所述模型关系建立横波速度预测模型。

优选地，在所述根据所述横波速度曲线、所述密度曲线和所述纵波速度曲线计算弹性参数曲线的步骤中，弹性参数曲线中至少包括拉梅常数、剪切模量、体积模量、杨氏模量和泊松比其中之一。

优选地，在所述基于石英含量测井曲线和叠后地震数据体得到石英含量数据体的步骤中，具体为利用石英含量测井曲线和叠后地震数据体之间的非线性关系，运用AVF反演得到石英含量数据体。

优选地，在所述基于石英含量数据体得到脆性指数的步骤中，具体计算公式如下：

其中，B_i表示脆性系指数，E_i表示地层的杨氏模量，E_min表示最小杨氏模量，E_max表示最大杨氏模量，S_i表示岩石的石英含量。

优选地，在所述基于叠前地震资料、纵波速度曲线、密度曲线得到时深曲线和地震子波的步骤中，具体为运用纵波速度曲线和密度曲线，依托由叠前地震资料获得的叠后地震资料进行标定，使得实际井旁道和合成地震记录达到匹配，此时的时深曲线和地震子波为所求的时深曲线和地震子波。

优选地，在所述基于所述时深曲线、所述密度曲线和所述弹性参数曲线建立初始模型，根据所述初始模型与地震子波褶积得到不同角度的合成角度道集，通过所述合成角度道集与叠前角道集建立反演方程进而反演得到杨氏模量的数据体的步骤中，基于所述时深曲线、所述密度曲线和所述弹性参数曲线中的杨氏模量曲线和泊松比曲线得到的层位数据通过插值外推的方法建立杨氏模量、泊松比和密度的初始模型，根据所述初始模型计算得到不同角度的反射系数，再与地震子波褶积得到不同角度的合成角度道集，通过所述合成角度道集与叠前角道集进行比较使两者的误差达到预设条件时建立反演方程，利用最小二乘的反演算法反演得到杨氏模量的数据体。

优选地，在所述根据所述初始模型计算得到不同角度的反射系数的步骤中，具体计算公式如下：

经转化得到，

其中，θ表示入射角，k表示横纵波速度比，E、σ、ρ分别表示杨氏模量的数据体、泊松比的数据体和裂缝密度的数据体。

优选地，在所述通过所述合成角度道集与叠前角道集进行比较使两者的误差达到预设条件时建立反演方程的步骤中，反演方程如下所示：

d＝Gm

其中，E、σ、ρ分别表示杨氏模量的数据体、泊松比的数据体和裂缝密度的数据体，θ₁至θ_n表示不同的入射角。

一种致密油气藏脆性的预测装置，包括：

横波速度预测模型模块，其用于基于获取的测井曲线得到横波速度与纵波速度曲线、密度曲线、GR曲线、电阻率曲线之间的模型关系，利用所述模型关系再建立横波速度预测模型；

横波速度曲线获取模块，其用于基于所述横波速度预测模型预测得到横波速度曲线；

弹性参数曲线计算模块，其用于根据所述横波速度曲线、所述密度曲线和所述纵波速度曲线计算弹性参数曲线；

脆性指数获取模块，其用于基于石英含量测井曲线和叠后地震数据体得到石英含量数据体，再根据所述石英含量数据体得到脆性指数；

时深曲线和地震子波计算模块，其用于基于叠前地震资料、纵波速度曲线、密度曲线得到时深曲线；

杨氏模量的数据体获取模块，其用于基于所述时深曲线、所述密度曲线和所述弹性参数曲线建立初始模型，根据所述初始模型与地震子波褶积得到不同角度的合成角度道集，通过所述合成角度道集与叠前角道集建立反演方程进而反演得到杨氏模量的数据体；

脆性指数的数据体获取模块，其用于基于石英含量数据体和杨氏模量的数据体得到脆性指数的数据体。

本发明的技术方案具有以下显著有益效果：

1、本发明实施例中的致密油气藏脆性的预测方法及装置通过构建估算横波速度预测模型，利用具有明确物理意义的测井曲线，如反映孔隙信息的声波和密度曲线、反映放射特征的GR曲线、反映岩性的自然电位等来估算横波速度，克服了经典的Biot-Gassman理论、Xu-White模型和K-T模型等的输入数据通常不易直接获取的缺点，如岩石物理的测量数据和测井的解释数据，这些数据都存在一定的误差，所以本发明所使用的方法就有更好的适用性。

2、本发明实施例中的致密油气藏脆性的预测方法及装置以基于杨氏模量、泊松比和密度的PP波反射系数公式为基础，利用叠前角道集直接估算杨氏模量、泊松比等脆性指数，如此能够有效提高反演的精度，计算的脆性指数的数据体更加精确。

附图说明

在此描述的附图仅用于解释目的，而不意图以任何方式来限制本发明公开的范围。另外，图中的各部件的形状和比例尺寸等仅为示意性的，用于帮助对本发明的理解，并不是具体限定本发明各部件的形状和比例尺寸。本领域的技术人员在本发明的教导下，可以根据具体情况选择各种可能的形状和比例尺寸来实施本发明。

图1为本发明在实施例中的流程图。

图2为本发明实施例中经测井资料预处理后，利用原始测井曲线以不同的方法计算的横波速度。

图3a和图3b分别为本发明实施例中纵波阻抗和横波阻抗交会以及纵波阻抗和杨氏模量交会图。

图4a、图4b、图4c和图4d分别为本发明实施例中泊松比与横波速度、剪切模量与拉梅系数、杨氏模量与体积模量、杨氏模量与泊松比的多维交会分析图。

图5为本发明实施例中建立的岩石物理量版定量解释储层。

图6a为本发明实施例中利用测井曲线构建脆性指数曲线。

图6b为本发明实施例中利用脆性指数与横波阻抗的交会关系识别含油气砂岩和干砂岩的示意图。

图7a为本发明实施例中I类井的井旁道叠前CRP道集数据。

图7b为本发明实施例中II类井的井旁道叠前角CRP道集数据。

图7c为本发明实施例中III类井的井旁道叠前角CRP道集数据。

图8a为本发明实施例中经过叠前地震直接反演得到的杨氏模量的连井线剖面。

图8b为本发明实施例中经过叠前地震直接反演得到的杨氏模量在目的层段的平面图。

图9a为本发明实施例中经过叠前地震直接反演得到的泊松比的连井线剖面。

图9b为本发明实施例中经过叠前地震直接反演得到的泊松比在目的层段的平面图。

图10a为本发明实施例中利用反演的杨氏模量和泊松比计算的脆性指数的数据体的过井剖面。

图10b为本发明实施例中利用反演的杨氏模量和泊松比计算的脆性指数的数据体的平面图。

图11为本发明实施例中致密油气藏脆性的预测装置的结构图。

具体实施方式

结合附图和本发明具体实施方式的描述，能够更加清楚地了解本发明的细节。但是，在此描述的本发明的具体实施方式，仅用于解释本发明的目的，而不能以任何方式理解成是对本发明的限制。在本发明的教导下，技术人员可以构想基于本发明的任意可能的变形，这些都应被视为属于本发明的范围。

为了提高反演的精度，提高致密油气藏脆性预测的精确性，在本申请提供了一种致密油气藏脆性的预测方法，图1为本发明在实施例中的流程图，如图1所示，本致密油气藏脆性的预测方法包括以下步骤：

S101：基于获取的测井曲线得到横波速度与纵波速度曲线、密度曲线、GR曲线、电阻率曲线之间的模型关系，利用模型关系再建立横波速度预测模型。

在本实施方式中，对密度曲线、GR曲线、电阻率曲线、声波时差等进行环境校正和标准化处理，通过环境校正消除测井仪器、温度、压力和泥浆等因素的影响，测井资料的标准化可将整个工区中的所有测井数据校正到同一刻度下，进一步提高测井资料的可靠性。根据获取的测井曲线得到横波速度，构建横波速度和纵波速度曲线、密度曲线、GR曲线、电阻率曲线之间的模型关系，利用模型关系建立横波速度预测模型。在一个实施方式中，可以通过神经网络算法拟合一个横波速度与纵波速度曲线、密度曲线、GR曲线、电阻率曲线等参数的对应关系式，该对应关系式(也就是所谓的模型关系)即可以作为横波速度预测模型。

S102：基于横波速度预测模型预测得到横波速度曲线。

通过构建的横波速度和纵波速度曲线、密度曲线、GR曲线、电阻率曲线等之间的模型关系，然后利用该模型关系预测横波速度曲线。图2为本发明实施例中经测井资料预处理后，利用原始测井曲线以不同的方法计算的横波速度，如图2所示，图2中第1、2、3栏分别是经标准化处理后的自然电位曲线、GR曲线、电阻率曲线、声波曲线和密度曲线，第4、5、6、7是预测的横波速度与实测横波速度曲线，其分别是通过纵波速度拟合、纵波速度和密度曲线拟合、概率神经网络计算和模型法预测得到的横波速度。从图2中可以看出通过预测模型预测的横波速度较其他方法更加精确，具有更明显的优势，可以为岩石物理分析和叠前反演提供可靠的资料基础。

S103：根据横波速度曲线、密度曲线和纵波速度曲线计算弹性参数曲线。

根据横波速度曲线、密度曲线和纵波速度曲线计算弹性参数曲线，弹性参数曲线中可以至少包括拉梅常数、剪切模量、体积模量、杨氏模量和泊松比其中之一，该步骤可以从众多弹性参数中优选出对储层敏感的参数。该步骤的计算过程可以如下：

其中：ρ表示密度，Vp表示纵波速度，Vs表示横波速度，E表示杨氏模量，K表示体积模量，μ表示剪切模量，λ表示拉梅常数，σ表示泊松比。

然后，通过多维岩石物理的分析技术对每口井的储层和邻段进行交会综合分析，进而指导储层预测和油气检测。在本实施例中，对纵波阻抗和横波阻抗、纵波阻抗和杨氏模量进行交会分析，图3a和图3b分别为本发明实施例中纵波阻抗和横波阻抗交会以及纵波阻抗和杨氏模量交会图，如图3a和图3b所示，可以辨别出含气砂岩、砂岩和泥岩等，图中黑色点为含气砂岩，白色点是砂岩，灰色点为泥岩，通过该交会显示表明，横波阻抗和杨氏模量与纵波阻抗的交会图可以明显地区分储层的岩性，从而识别有利储层。

在岩性识别的基础上，利用多分岩石物理技术分析储层的含油气性，图4a、图4b、图4c和图4d分别为本发明实施例中泊松比与横波速度、剪切模量与拉梅系数、杨氏模量与体积模量、杨氏模量与泊松比的多维交会分析图，如图4a、图4b、图4c和图4d中所示，通过优选出泊松比、拉梅系数、杨氏模量和体积模量交会等敏感参数可以有效区分储层的含油气性。图5为本发明实施例中建立的岩石物理量版定量解释储层，其用于定量识别储层的含气性。

S104：基于石英含量测井曲线和叠后地震数据体得到石英含量数据体。

利用已知的泥质含量计算石英含量，进而形成石英含量测井曲线，石英含量的计算公式如下：

V_quartz＝1-V_sh

其中，V_quartz为石英含量，V_sh为泥质含量。

利用石英含量测井曲线和叠后地震数据体之间的非线性关系，运用现有技术AVF(Amplitude Versus Frequency)反演得到石英含量数据体。

分频(AVF)反演技术是近几年发展起来的一种先进的地震反演技术，依靠测井和地震数据，通过研究不同地层厚度下的振幅与频率之间的关系(AVF)，将AVF作为独立信息引入反演，合理利用地震的低、中、高频带信息，减少薄层反演的不确定性，得到一个高分辨率的反演结果。同时它也是一种无子波提取，无初始模型的高分辨率非线性反演。

S105：基于石英含量数据体得到脆性指数。

根据石英含量数据体计算得到脆性指数，具体计算公式如下：

其中，B_i表示脆性系指数，E_i表示地层的杨氏模量，E_min表示最小杨氏模量，E_max表示最大杨氏模量，S_i表示岩石的石英含量。在本实施方式中，图6a为本发明实施例中利用测井曲线构建脆性指数曲线，如图6a所示，利用杨氏模量的变化率和岩石矿物含量(石英含量数据体)计算得到岩石的脆性指数的曲线，其具体为图6a中间左边的虚线。再将脆性指数和横波阻抗进行交会分析，图6b为本发明实施例中利用脆性指数与横波阻抗的交会关系识别含油气砂岩和干砂岩的示意图，如图6b所示，可以发现横波阻抗与脆性指数的交会图可以准确地识别出含油气砂岩、干砂岩，如此说明含油气储层的脆性指数较大。

S106：基于叠前地震资料、纵波速度曲线、密度曲线得到时深曲线，其包括如下步骤：

S201：基于叠前地震资料得到叠后地震资料和叠前角道集。

通过对叠前地震资料保幅处理得到叠后地震资料和叠前角道集(即叠前CRP道集)。针对致密储层，可以对叠前角道集进行提高信噪比处理和高频信息补偿处理，最终可以得到高保真度、高分辨率和高信噪比的叠前地震资料，为致密储层的脆性特征预测提供保障。在本实施方式中，图7a为本发明实施例中I类井的井旁道叠前CRP道集数据，图7b为本发明实施例中II类井的井旁道叠前角CRP道集数据，图7c为本发明实施例中III类井的井旁道叠前角CRP道集数据，如图7a、7b、7c所示，其分别是具有AVO特征的叠前道集，可以看到：图7a中目的层段的振幅随偏移距增大而增大，AVO现象明显；图7b的目的层段的振幅也随偏移距增大而增大，AVO现象明显；图7c目的层段的振幅随偏移距发生变化但不明显，AVO异常小。所以，通过上述结构发现可以利用I类和II类AVO道集的特征识别目的层段的含气性。

S202：基于叠后地震资料、纵波速度曲线、密度曲线得到时深曲线和地震子波。

运用纵波速度和密度曲线，依托由叠前地震资料获得的叠后地震资料进行标定，使得实际井旁道和合成地震记录达到最佳匹配，最佳匹配时的时深曲线和地震子波即为所求的时深曲线和地震子波。具体而言，首先在叠后地震数据的井旁道中提取地震子波，然后利用纵波阻抗(纵波速度和密度的乘积)曲线计算反射系数，与子波褶积得到合成地震记录，然后与实际井旁道进行对比，反复重复上述，直到达到最佳匹配。

S107：基于时深曲线、密度曲线和弹性参数曲线建立初始模型，根据初始模型与地震子波褶积得到不同角度的合成角度道集，通过合成角度道集与叠前角道集建立反演方程进而反演得到杨氏模量的数据体。

基于时深曲线、密度曲线和弹性参数曲线中的杨氏模量曲线和泊松比曲线得到的层位数据通过插值外推的方法建立杨氏模量、泊松比和密度的初始模型，根据初始模型计算得到不同角度的反射系数，具体计算公式如下：

经转化得到：

其中，θ表示入射角，θ₁、…、θ_n表示不同的入射角，k表示横纵波速度比，E、σ、ρ分别表示杨氏模量的数据体、泊松比的数据体和裂缝密度的数据体，R_pp表示反射系数，其可以为纵波反射系数。

不同角度的反射系数再与地震子波褶积得到不同角度的合成角度道集，通过合成角度道集与叠前角道集进行比较使两者的误差达到预设条件时建立反演方程，预设条件可以为两者的误差达到最小值。

利用最小二乘的反演算法反演可以得到杨氏模量的数据体、泊松比的数据体、裂缝密度的数据体，反演方程如下所示：

d＝Gm

其中，E、σ、ρ分别表示杨氏模量的数据体、泊松比的数据体和裂缝密度的数据体，θ₁至θ_n表示不同的入射角。

在本实施方式中，图8a为本发明实施例中经过叠前地震直接反演得到的杨氏模量的连井线剖面，图8b为本发明实施例中经过叠前地震直接反演得到的杨氏模量在目的层段的平面图，图9a为本发明实施例中经过叠前地震直接反演得到的泊松比的连井线剖面，图9b为本发明实施例中经过叠前地震直接反演得到的泊松比在目的层段的平面图。

S108：基于石英含量数据体和杨氏模量的数据体得到脆性指数的数据体。

根据石英含量数据体和杨氏模量的数据体计算得到脆性指数的数据体，最终运用该脆性指数的数据体来评价致密储层的脆性特征。通过杨氏模量的变化率和岩石矿物含量可以表示岩石的脆性关系，图10a为本发明实施例中利用反演的杨氏模量和泊松比计算的脆性指数的数据体的过井剖面，图10b为本发明实施例中利用反演的杨氏模量和泊松比计算的脆性指数的数据体的平面图，如图10a所示，从剖面中脆性指数相对较高的部分(图10a中黑圈所示的部分)为有利的含气储层，可以计算出储层的有效厚度为19米，与已知的产气井具有很好的一致性，该井为产气井，日产气1.44×10⁴方，由此可以认为产气较高时具有较高的储层厚度，同时具有较大的脆性指数，产气较低时储层厚度相应变低，脆性指数变小，如图10b所示，平面结果与测井解释结果相符合，由此证明该方法具有一定的实用性。

通过本发明的实施例可以看出基于预测模型预测的横波速度与实测的横波速度具有很好的一致性，如图2所示。如图3、图4所示，采用岩石物理交会和多维岩石物理技术能够更好地区分储层特征和含油气性质。如图8a、8b、9a、9b、10a、10b所示，弹性参数的反演结果能够准确地识别岩性，砂岩在纵横波速度比、泊松比剖面呈中低值，表现为黑灰色(如图8a中B箭头所指)和白色(如图8a中A箭头所指)分布，在平面图中，黑色(如图8b中C箭头所指)和白色(如图8b中D箭头所指)代表砂岩储层，可以看出，弹性参数反映的储层分布规律性较好，边界清晰，且不同弹性参数之间的一致性较好；脆性指数与储层的含气性密切相关，从脆性指数剖面可以看出，脆性指数越高，目的层段的有效储层越厚，同时具有更高的含气性，这与测井解释结果完全一致，脆性指数剖面很好地反映了储层的脆性特征和含气性特征，结果证实了该发明的有效性和适用性。

本申请还提出了一种致密油气藏脆性的预测装置，图11为本发明实施例中致密油气藏脆性的预测装置的结构图，如图11所示，该装置包括：

横波速度预测模型模块，其用于基于获取的测井曲线得到横波速度与纵波速度曲线、密度曲线、GR曲线、电阻率曲线之间的模型关系，利用模型关系再建立横波速度预测模型；

横波速度曲线获取模块，其用于基于横波速度预测模型预测得到横波速度曲线；

弹性参数曲线计算模块，其用于根据横波速度曲线、密度曲线和纵波速度曲线计算弹性参数曲线；

脆性指数获取模块，其用于基于石英含量测井曲线和叠后地震数据体得到石英含量数据体，再根据石英含量数据体得到脆性指数；

时深曲线和地震子波计算模块，其用于基于叠前地震资料、纵波速度曲线、密度曲线得到时深曲线和地震子波；

杨氏模量的数据体获取模块，其用于基于时深曲线、密度曲线和弹性参数曲线建立初始模型，根据初始模型与地震子波褶积得到不同角度的合成角度道集，通过合成角度道集与叠前角道集建立反演方程进而反演得到杨氏模量的数据体；

脆性指数的数据体获取模块，其用于基于石英含量数据体和杨氏模量的数据体得到脆性指数的数据体。

控制器可以按任何适当的方式实现，例如，控制器可以采取例如微处理器或处理器以及存储可由该(微)处理器执行的计算机可读程序代码(例如软件或固件)的计算机可读介质、逻辑门、开关、专用集成电路(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)，如调制的数据信号和载波。

本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。

本申请可用于众多通用或专用的计算机系统环境或配置中。例如：个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、多处理器系统、基于微处理器的系统、置顶盒、可编程的消费电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等等。

本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。

虽然通过实施例描绘了本申请，本领域普通技术人员知道，本申请有许多变形和变化而不脱离本申请的精神，希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。

Claims (11)

1.一种致密油气藏脆性的预测方法，其特征在于，包括以下步骤：

基于获取的测井曲线得到横波速度与纵波速度曲线、密度曲线、GR曲线、电阻率曲线之间的模型关系，利用所述模型关系再建立横波速度预测模型；

基于所述横波速度预测模型预测得到横波速度曲线；

根据所述横波速度曲线、所述密度曲线和所述纵波速度曲线计算弹性参数曲线；

基于石英含量测井曲线和叠后地震数据体得到石英含量数据体；

基于所述石英含量数据体得到脆性指数；

基于叠前地震资料、纵波速度曲线、密度曲线得到时深曲线和地震子波；

基于所述时深曲线、所述密度曲线和所述弹性参数曲线建立初始模型，根据所述初始模型与地震子波褶积得到不同角度的合成角度道集，通过所述合成角度道集与叠前角道集建立反演方程进而反演得到杨氏模量的数据体；

基于所述石英含量数据体和所述杨氏模量的数据体得到脆性指数的数据体。

2.根据权利要求1所述的致密油气藏脆性的预测方法，其特征在于，在所述基于叠前地震资料、纵波速度曲线、密度曲线得到时深曲线的步骤中具体包括：

基于叠前地震资料得到叠后地震资料；

基于所述叠后地震资料、纵波速度曲线、密度曲线得到时深曲线。

3.根据权利要求1所述的致密油气藏脆性的预测方法，其特征在于，在所述基于获取的测井曲线得到横波速度与纵波速度曲线、密度曲线、GR曲线、电阻率曲线之间的模型关系，利用所述模型关系再建立横波速度预测模型的步骤中，具体为：对密度曲线、GR曲线、电阻率曲线进行环境校正和标准化处理，根据获取的测井曲线得到横波速度，构建所述横波速度和纵波速度曲线、密度曲线、GR曲线、电阻率曲线之间的模型关系，利用所述模型关系建立横波速度预测模型。

4.根据权利要求1所述的致密油气藏脆性的预测方法，其特征在于，在所述根据所述横波速度曲线、所述密度曲线和所述纵波速度曲线计算弹性参数曲线的步骤中，弹性参数曲线中至少包括拉梅常数、剪切模量、体积模量、杨氏模量和泊松比其中之一。

5.根据权利要求1所述的致密油气藏脆性的预测方法，其特征在于，在所述基于石英含量测井曲线和叠后地震数据体得到石英含量数据体的步骤中，具体为利用石英含量测井曲线和叠后地震数据体之间的非线性关系，运用AVF反演得到石英含量数据体。

6.根据权利要求1所述的致密油气藏脆性的预测方法，其特征在于，在所述基于石英含量数据体得到脆性指数的步骤中，具体计算公式如下：

$B_i=\frac{E_i-E_{\min }}{E_{max}-E_i}\·S_i$

其中，B_i表示脆性系指数，E_i表示地层的杨氏模量，E_min表示最小杨氏模量，E_max表示最大杨氏模量，S_i表示岩石的石英含量。

7.根据权利要求1所述的致密油气藏脆性的预测方法，其特征在于，在所述基于叠前地震资料、纵波速度曲线、密度曲线得到时深曲线和地震子波的步骤中，具体为运用纵波速度曲线和密度曲线，依托由叠前地震资料获得的叠后地震资料进行标定，使得实际井旁道和合成地震记录达到匹配，此时的时深曲线和地震子波为所求的时深曲线和地震子波。

8.根据权利要求1所述的致密油气藏脆性的预测方法，其特征在于，在所述基于所述时深曲线、所述密度曲线和所述弹性参数曲线建立初始模型，根据所述初始模型与地震子波褶积得到不同角度的合成角度道集，通过所述合成角度道集与叠前角道集建立反演方程进而反演得到杨氏模量的数据体的步骤中，基于所述时深曲线、所述密度曲线和所述弹性参数曲线中的杨氏模量曲线和泊松比曲线得到的层位数据通过插值外推的方法建立杨氏模量、泊松比和密度的初始模型，根据所述初始模型计算得到不同角度的反射系数，再与地震子波褶积得到不同角度的合成角度道集，通过所述合成角度道集与叠前角道集进行比较使两者的误差达到预设条件时建立反演方程，利用最小二乘的反演算法反演得到杨氏模量的数据体。

9.根据权利要求8所述的致密油气藏脆性的预测方法，其特征在于，在所述根据所述初始模型计算得到不同角度的反射系数的步骤中，具体计算公式如下：

$\begin{array}{c}{R}_{pp}\left(\mathrm{\θ}\right)=(\frac{1}{4}{\sec }^2\mathrm{\θ}-2k^2{\sin }^2\mathrm{\θ})\frac{\mathrm{\Δ}E}{\stackrel{\‾}{E}}+\[\frac{1}{4}{\sec }^2\mathrm{\θ}\frac{(2k^2-3){(2k^2-1)}^2}{k^2(4k^2-3)}\\ +2k^2{\sin }^2\mathrm{\θ}\frac{1-2k^2}{3-4k^2}\]\frac{\mathrm{\Δ}\mathrm{\σ}}{\stackrel{\‾}{\mathrm{\σ}}}+(\frac{1}{2}-\frac{1}{4}{\sec }^2\mathrm{\θ})\frac{\mathrm{\Δ}\mathrm{\ρ}}{\stackrel{\‾}{\mathrm{\ρ}}}\end{array}$

经转化得到，

其中，θ表示入射角，k表示横纵波速度比，E、σ、ρ分别表示杨氏模量的数据体、泊松比的数据体和裂缝密度的数据体，R_pp表示反射系数。

10.根据权利要求9所述的致密油气藏脆性的预测方法，其特征在于，在所述通过所述合成角度道集与叠前角道集进行比较使两者的误差达到预设条件时建立反演方程的步骤中，反演方程如下所示：

d＝Gm

$\begin{array}{ccc}d=\left[\begin{array}{c}{R}_{pp}\left({\mathrm{\θ}}_1\right)\\ R_{pp}\left({\mathrm{\θ}}_2\right)\\ \·\\ \·\\ \·\\ R_{pp}\left({\mathrm{\θ}}_n\right)\end{array}\right]& G=\left[\begin{array}{ccc}{C}_1\left({\mathrm{\θ}}_1\right)& C_2\left({\mathrm{\θ}}_1\right)& C_3\left({\mathrm{\θ}}_1\right)\\ C_1\left({\mathrm{\θ}}_2\right)& C_2\left({\mathrm{\θ}}_2\right)& C_3\left({\mathrm{\θ}}_2\right)\\ \·& \·& \·\\ \·& \·& \·\\ \·& \·& \·\\ C_1\left({\mathrm{\θ}}_n\right)& C_2\left({\mathrm{\θ}}_n\right)& C_3\left({\mathrm{\θ}}_n\right)\end{array}\right]& m=\left[\begin{array}{c}\frac{\mathrm{\Δ}E}{\stackrel{\‾}{E}}\\ \frac{\mathrm{\Δ}\mathrm{\σ}}{\stackrel{\‾}{\mathrm{\σ}}}\\ \frac{\mathrm{\Δ}\mathrm{\ρ}}{\stackrel{\‾}{\mathrm{\ρ}}}\end{array}\right]\end{array}$

其中，E、σ、ρ分别表示杨氏模量的数据体、泊松比的数据体和裂缝密度的数据体，θ₁至θ_n表示不同的入射角。

11.一种致密油气藏脆性的预测装置，其特征在于，包括：

横波速度预测模型模块，其用于基于获取的测井曲线得到横波速度与纵波速度曲线、密度曲线、GR曲线、电阻率曲线之间的模型关系，利用所述模型关系再建立横波速度预测模型；

横波速度曲线获取模块，其用于基于所述横波速度预测模型预测得到横波速度曲线；

弹性参数曲线计算模块，其用于根据所述横波速度曲线、所述密度曲线和所述纵波速度曲线计算弹性参数曲线；

脆性指数获取模块，其用于基于石英含量测井曲线和叠后地震数据体得到石英含量数据体，再根据所述石英含量数据体得到脆性指数；

时深曲线和地震子波计算模块，其用于基于叠前地震资料、纵波速度曲线、密度曲线得到时深曲线；

杨氏模量的数据体获取模块，其用于基于所述时深曲线、所述密度曲线和所述弹性参数曲线建立初始模型，根据所述初始模型与地震子波褶积得到不同角度的合成角度道集，通过所述合成角度道集与叠前角道集建立反演方程进而反演得到杨氏模量的数据体；

脆性指数的数据体获取模块，其用于基于石英含量数据体和杨氏模量的数据体得到脆性指数的数据体。