确定toc含量的方法和装置的制造方法
【技术领域】
[0001] 本发明涉及石油地球物理勘探技术，具体涉及一种确定T0C(Total Organic Carbon，总有机碳）含量的方法和装置。
【背景技术】
[0002] 目前，主要是利用地震手段来预测致密油气藏，致密油气具有两个特点：1)由于 岩性致密，孔隙度和渗透率低，油气从烃源岩往外排、运的距离非常短，因此往往储存在靠 近烃源岩的储层中，即油气具有"短排、短运"的特点；2)由于致密岩性往往成分复杂，自身 既是烃源岩又是储层，因此油气生成后有时就地存储，即油气具有"自生、自储"的特点。致 密油气的以上两个特点决定了致密油气藏预测的一项主要内容就是岩石中表征有机质丰 度指标的预测。
[0003] 总有机碳（Total Organic Carbon，T0C)含量指单位质量岩石中有机碳的质量百 分数，是以有机碳的含量表示有机质总量的一个指标。有机成因论认为石油和天然气都是 通过有机质转化而来的，因此将富含有机质的岩性称为烃源岩。由于有机质主要是碳水化 合物，因此TOC含量可以作为测定有机质含量的一个重要指标。在石油勘探中，TOC含量的 测定往往在实验室进行，是石油地质学家判断烃源岩有机质富集程度的一个重要指标。
[0004] 近年来，随着非常规油气资源的大力勘探与开发，通过地球物理技术估算TOC含 量分布成为预测有效烃源岩分布的重要手段。目前利用地球物理手段预测TOC含量的技术 主要有以下两种方式：
[0005] 方式1)实验室测定得到岩样的TOC含量数值后，与对应深度的测井曲线进行对比 分析，定性确定TOC含量对应的测井响应特征，这是为了实现由测井响应特征预测TOC含量 的目的。进一步利用地震属性与测井响应之间的关系确定测井响应对应的地震属性特征， 从而实现由地震属性间接预测TOC含量。
[0006] 这种方式与单一属性预测方法相比，可靠性更高，缺点是测井曲线统计拟合得到 的TOC曲线误差较大，与地震属性之间建立联系的方法过于简单，误差也非常大。
[0007] 方式2)结合地质分析，利用地震数据直接预测TOC含量的空间分布。这种方式与 上一种方式相比，最大的特点是以地质认识为背景，能够从宏观上把握预测结果的有效性， 但是纵向精度很低，预测结果误差更大。
[0008] 针对如何提高TOC含量的预测精度，目前尚未提出有效的解决方案。
【发明内容】
[0009] 本发明实施例提供了一种确定TOC含量的方法，以达到提高TOC含量的预测精度 的目的，该方法包括：
[0010] 获取待分析区域的测井曲线和实验室测定数据；
[0011] 根据所述测井曲线和所述实验室测定数据，估算所述待分析区域的孔隙度和含水 饱和度；
[0012] 对从所述待分析区域中获取的岩样进行X射线衍射数据分析得到XRD数据；
[0013] 对所述岩样进行分析，得到所述待分析区域岩样的TOC含量测定数据；
[0014] 获取所述岩样的体积和质量，并根据所述体积和质量计算得到无机矿物密度和 TOC密度；
[0015] 在所述待分析区域的井点处，确定所述测井曲线与所述岩样的TOC含量测定数据 之间的对应关系，并将所述对应关系作为所述井点处TOC含量数据的预测模型，根据所述 预测模型拟合得到一条连续的TOC密度曲线；
[0016] 建立由无机矿物密度和TOC密度表征的岩石骨架模型，其中，所述无机矿物密度 由所述XRD数据减去所述TOC含量测定数据得到，所述TOC密度为所述TOC密度曲线表征 的数值；
[0017] 根据所述测井曲线、所述孔隙度、所述含水饱和度、所述无机矿物密度和所述TOC 密度，建立所述岩石骨架模型表征的岩石密度；
[0018] 根据所述岩石骨架模型表征的岩石密度和所述测井曲线中的密度曲线，建立目标 函数，利用最优化分析，通过迭代校正的方式对所述TOC密度曲线进行校正，得到校正后的 TOC密度曲线；
[0019] 根据所述岩石骨架模型建立速度与弹性模量之间的定量关系；
[0020] 根据实验室对所述岩样测定得到的速度、弹性模量、TOC数据，和所述定量关系确 定微观尺度对TOC敏感的弹性参数，建立基于敏感弹性参数的TOC量版；
[0021] 根据所述TOC量版和所述测井曲线计算所述井点处对TOC敏感的弹性参数，并根 据所述井点处对TOC敏感的弹性参数和所述校正后的TOC密度曲线，从中观尺度计算得到 弹性参数与TOC密度之间的相关度的统计关系；
[0022] 获取所述待分析区域的地震数据，并对所述地震数据进行OVT域处理；
[0023] 对进行OVT域处理后的地震数据和所述测井曲线进行叠前地震反演，得到所述对 TOC敏感的弹性参数的数据体；
[0024] 根据所述统计关系，将所述对TOC敏感的弹性参数的数据体转换为TOC含量空间 分布数据体；
[0025] 根据所述TOC含量空间分布数据体得到TOC含量空间分布。
[0026] 在一个实施方式中，获取待分析区域的测井曲线，包括：
[0027] 获取所述待分析区域的一个或多个测井资料，其中，所述测井资料包括以下一种 或任意组合：纵波时差、横波时差、密度、中子、电阻率、自然伽马、自然伽马能谱铀；
[0028] 根据获得的所述一个或多个测井资料生成所述测井曲线。
[0029] 在一个实施方式中，根据所述测井曲线，估算所述待分析区域的孔隙度和含水饱 和度，包括：
[0030] 根据所述测井曲线，建立孔隙度估算模型和含水饱和度估算模型：
[0031]
[0032]
[0033] 其中，炉表示总孔隙度，Sw表示为含水饱和度，X i，X2，…\表示用于进行孔隙度估 算的实验测定数据和测井曲线，yi，y 2，…yn表示参与估算含水饱和度的测井曲线，η表示参 与估算的数据的个数，爲和fsw表示多元映射函数；
[0034] 根据所述孔隙度估算模型估算得到孔隙度，根据所述含水饱和度估算模型估算得 到含水饱和度。
[0035] 在一个实施方式中，根据所述体积和质量计算得到无机矿物密度和TOC密度，包 括：
[0036] 按照以下公式计算无机矿物密度和TOC密度：
[0037]
[0038]
[0039]
[0040] 其中，Q表示干岩样的质量，Qnk表示为无机矿物质量，Q T(K表示TOC质量，P nk表示 无机矿物密度，P 表示TOC密度，V表示干岩样的体积。
[0041] 在一个实施方式中，根据所述岩石密度和所述测井曲线中的密度曲线，建立目标 函数，利用最优化分析，通过迭代校正TOC密度校正岩石密度，以得到校正后的TOC密度曲 线，包括：
[0042] 根据所述岩石骨架模型建立响应方程：
[0043]
[0044] 其中，P表示响应方程计算的岩石体密度理论值，K表示TOC体积百分比，P nk表 示无机矿物密度，Ρτ。。表示TOC密度，@表示总孔隙度，P。表示油密度，P w表示水密度，Sw 表示含水饱和度；
[0045] 根据所述岩石密度和所述测井曲线中的密度曲线，建立如下的目标函数：
[0046]
[0047]
[0048]
[0049] 其中，F(X)表示目标函数，Pb表示实际测井密度，g(x)表示约束项，min表示通 过约束项使得目标函数达到极小值；
[0050] 以TOC密度作为所述目标函数中的X进行迭代校正，得到校正后的TOC密度曲线。
[0051] 在一个实施方式中，根据所述岩石骨架模型建立速度与弹性模量之间的定量关 系，包括：
[0052] 通过Hashin-Shtrikman界限确定所述岩石骨架模型的体积模量和剪切模量；
[0053] 通过K-T公式计算干岩石的体积模量和剪切模量；
[0054] 通过Gassmann方程计算饱含水情况下岩石的体积模量和剪切模量；
[0055] 通过饱含水情况下岩石的体积模量和剪切模量，建立速度与弹性模量之间的定量 关系。
[0056] 本发明实施例还提供了一种确定TOC含量的装置，以达到提高TOC含量的预测精 度的目的，该装置包括：
[0057] 数据获取模块，用于获取待分析区域的测井曲线和实验室测定数据；
[0058] 估算模块，用于根据所述测井曲线和所述实验室测定数据，估算所述待分析区域 的孔隙度和含水饱和度；
[0059] 分析模块，用于对从所述待分析区域中获取的岩样进行X射线衍射数据分析得到 XRD数据；
[0060] TOC含量测定模块，用于对所述岩样进行分析，得到所述待分析区域岩样的TOC含 量测定数据；
[0061] 密度确定模块，用于获取所述岩样的体积和质量，并根据所述体积和质量计算得 到无机矿物密度和TOC密度；
[0062] TOC密度曲线拟合模块，用于在所述待分析区域的井点处，确定所述测井曲线与所 述岩样的TOC含量测定数据之间的对应关系，并将所述对应关系作为所述井点处TOC含量 数据的预测模型，根据所述预测模型拟合得到一条连续的TOC密度曲线；
[0063] 岩石骨架模型建立模块，用于建立由无机矿物密度和TOC密度表征的岩石骨架模 型，其中，所述无机矿物密度由所述XRD数据减去所述TOC含量测定数据得到，所述TOC密 度为所述TOC密度曲线表征的数值；
[0064] 岩石密度确定模块，用于根据所述测井曲线、所述孔隙度、所述含水饱和度、所述 无机矿物密度和所述TOC密度，建立所述岩石骨架模型表征的岩石密度；<