一种基于数据重构预测岩性体的方法
【技术领域】
[0001] 本发明涉及地球物理勘探中的地震资料解释处理技术领域，具体的说是一种利用 数据体重构来预测不同岩性体的方法。
【背景技术】
[0002] 目前，在油气勘探中岩性预测是一个比较重要的环节，不同的岩性体对油气的储 集能力基本上不一样，如高孔隙度的砂岩及溶孔发育的礁灰岩对油气具有良好的储集能 力，所以在石油勘探中是相当重要的勘探目标。如在四川盆地探区中的海相及湖相地层中 寻找这类目标对该区的油气勘探具有重要意义，钻遇此类岩性的钻井大多数经测试均获工 业气流。
[0003] 常规的岩性体预测技术方法有很多种，基本上可分为叠前、叠后反演两大类，经过 反演或属性计算可得到针对岩性体预测的属性数据体。但对于一个研究区内有多种类型的 岩性体时，往往此时的单一属性体不能很好地对目的层段中的多种岩性体进行描述，预测 结果往往具有很强的片面性;也有的方法是通过井上的岩性体、叠后波阻抗反演数据体进 行岩性体解释，得到相关岩性体空间上的分布情况。一些专利如发明《基于数学形态学的岩 性体自动追踪方法》（专利【申请号】201410154418.8)的专利公开提供一种基于数学形态学 的岩性体自动追踪方法，该方法包括选取用于岩性体自动追踪的地震属性体，并进行数据 筛选;数据筛选结束后，检测保留下来的数据点的空间位置信息，对这些数据点进行砂体分 类和编号;利用开运算对处理后得到的砂体边界做处理，断开两个粘连砂体间的窄通道，实 现砂体的分离；以及根据最小砂体大小的设置，对应用开运算处理之后得到的砂体进行重 新分类编号。可见近年来，对岩性体预测的技术方法很多，研究相对比较成熟。但是，针对岩 性体的预测技术还不够完善，主要体现在以下几个方面：
[0004] (1)采用单一数据体进行不同岩性体的预测结果往往针对性不强，准确度不高。 [0005] (2)利用数据交会分析实施不同岩性体预测往往带有人为主观因素，并且预测的 准确率不高。
[0006] 为了针对上述传统岩性体描述技术方法的不足，解决对同一工区内的多种岩性体 进行相应的区分及描述，本发明提供了一个相对简单的技术流程来实现对工区内的多种岩 性体进行预测，更方便、准确地了解不同岩性体的分布情况，从而更好的为油气勘探服务。
【发明内容】
[0007] 本发明目的在于克服现有常规岩性体预测技术中的上述和其他缺点，可以准确地 预测出不同岩性体的平面及垂向上的分布状态。为此，本发明提供一种利用数据体重构来 预测不同岩性体分布情况的方法。
[0008] 本发明的技术方案是：
[0009] -种基于数据重构预测岩性体的方法，它包括以下步骤：
[0010] S1、对待测区域的地震资料进行提取和反演，获取该区域多种岩性体的预测数据 体，并获得各岩性体所对应的预测数据体的门槛值；
[0011] S2、对步骤1所确定门槛值后的各个预测数据体进行数据重构后加权处理，得到一 个用于岩性预测的加权重构数据体，根据前述加权重构数据体获取待测区域的岩性。
[0012] 本发明的步骤S1中，预测数据体为利用叠前或叠后地震资料反演或属性提取后所 得到的数据体，或者为两个或以上的数据体之间的数学运算所得到的数据体，在用于岩性 体预测的数据体中，如果一个数据体能够实现岩性体预测时，则使用该数据体进行数据重 构、加权计算；否则应该使用两个或两个以上的数据体分别实施数据重构、加权处理后，得 到一个加权重构数据体用于岩性体预测。
[0013] 本发明的岩性体包括砂岩、泥岩、灰岩、页岩以及生物礁，前述岩性体的物理特征 各不相同，具有差异性。
[0014] 本发明的步骤S1中，不同岩性体的预测数据体门槛值设定包括以下步骤：
[0015] Sl-a、根据待测区域中各井的测井资料（由测井设备记录）、岩性资料(录井地质资 料)获取不同岩性体的深度范围数据，并根据测井资料从预测数据体中提取相关岩性的预 测数据体；（根据各井的测井资料及岩性资料统计及确定不同岩性体的深度范围数据，并根 据不同的岩性体的物理差异特征采取合适的地震资料、反演或属性计算方法提取相关的数 据体，要求相关预测数据体在数据值方面能将预测岩性体与其它非预测岩性体区分开来）
[0016] si-b、利用测井资料及岩性资料对各相关岩性的预测数据体进行井一震标定，确 定各井的时深关系；
[0017] S1-C、根据井一震标定所得到的时深关系将步骤si-a中获取的不同岩性体的深度 范围数据转换到相关数据体的时间域中，得到各岩性体所对应的预测数据体的数据值，并 综合各井相关的同一岩性体的数据构成数据集，从各岩性体的数据集中确定该岩性体预测 数据的门槛值范围，依此类推，从而实现不同岩性体的门槛值设置。
[0018] 本发明的步骤Sl-c中：对于各岩性体，从数据集中确定出最大及最小门槛值，其中 门槛值求取方式为综合各井相关的同一岩性体的数据构成数据集{Jd，门槛值求取的计算 公式如下：
[0019] Jmax=max{Ji} (1)
[0020] Jmin=min{Ji} (2)
[0021] 式(1) 一 (2)中心为井中该岩性段所对应数据体的第i个采样点的数据值，Jmax为该 岩性体的最大门槛值，J min为该岩性体的最小门槛值。
[0022]本发明的步骤S2具体为：
[0023] S2_a、采用步骤Sl-c获取的各岩性体的门槛值对步骤Sl-a获取的相关岩性的预测 数据体进行重构赋值得到重构数据体:对于不同的岩性体将各自门槛值内的数据值赋值为 某一常数值，不同的岩性体对应不同的常数值，各个门槛值外的数据均赋值为零值，以示对 岩性体的数据进行区分。依此类推，得到一系列各个岩性体对应的重构数据体。
[0024] S2_b、对各个岩性体的重构数据体进行加权处理，得到一个加权重构数据体。（首 先，针对不同的岩性体预测使用分配好的数据体进行，并设置相关岩性体的门槛值。要求一 个数据体可以用于预测一种或以上的岩性体，也可将不同的岩性体分开使用不同的数据体 进行相关的预测计算，要求不同的数据体不能预测同一岩性体;然后，不同的数据体根据分 配好预测的岩性体的门槛值进行重构赋值，门槛值外的数据设置为零值，门槛值内的数据 设置为某一个常数值，不同岩性体的常数值不能相同，并且两种或以上的代表岩性的常数 值进行加权计算时所得的结果不能与任一岩性体的常数值相等;最后，对各个数据体进行 重构赋值后针对同一道同一个采样点进行加权计算得到一个加权重构数据体，利用该数据 体进行岩性体预测。若当加权后所得到的数据值不与任何一个岩性体的常数值一致时，应 根据实际情况、地质资料及井上岩性资料进行再次重新赋值，赋值为某一个预测岩性体的 常数值。一般情况下，应根据该采样点的上、下、左及右采样点的常数值进行赋值为某一个 岩性体的常数值）
[0025] 本发明的步骤S2_a具体为:根据门槛值范围实施数据体的数据重构：
[0026] 若某一岩性的预测数据体中采样点的数据值在该岩性体的门槛值范围内，则将该 数据值设置为某一个常数K值;若采样点的数据值在门槛值范围外，则将该采样点的数据值 赋值为零即：
[0027] a、若某一个预测数据体的第i个采样点的数据值为土，假设当土<1_时，J_为该 数据体针对某一预测岩性体的最小门槛值，则& = 〇;
[0028] b、假设当Jmi"JdJmax时，Jmax为该数据体针对该预测岩性体的最大门槛值，则心 =K，要求K值一般情况下为一常数值，不同岩性体的常数K值不能相同，并且K值之间的组合 后加权结果不能与任何一个岩性体设定的Κ值一致；
[0029] c、假设当 Ji > Jmax时，则 Ji = 0。
[0030] 8、根据权利要求6所述的基于数据重构预测岩性体的方法，其特征是步骤S2_b具 体为:对各岩性体的重构数据体进行加权处理，加权处理的计算公式如下：
[0031] Μ:?=Σ