单元，用于通过归一化处理后的融合频率和所述三个融合频率各自对应的沉 积微相，对所述三个融合频率进行RGB融合，得到薄层沉积微相反演结果。
[0061] 在本发明实施例中，利用井点处的测井微相来标定RGB属性融合图，以得到沉积 微相的预测结果，具体地，通过测录井微相分析、测井解释、谱分解、交会图统计分析、RGB融 合等技术的应用，充分结合钻井、地质及地震信息，利用反演的思想，提出了一种薄层沉积 微相反演的方法，通过测井、录井相的约束提高了RGB融合技术在预测沉积微相时的合理 性，解决了现有技术中无法有效确定RGB融合技术中融合属性的三个频率，而导致预测得 到的沉积微相不准确的技术问题，达到了准确预测沉积微相的技术效果，进一步的，减少了 由于地质不确定性带来的地震属性预测技术的多解性。
【附图说明】
[0062] 此处所说明的附图用来提供对本发明的进一步理解，构成本申请的一部分，并不 构成对本发明的限定。在附图中：
[0063] 图1是根据本发明实施例的薄层沉积微相反演方法的方法流程图；
[0064] 图2是是根据本发明实施例的基于短时傅里叶变换的调谐频率连井剖面图；
[0065] 图3是根据本发明实施例的经过最大熵求取的调谐频率连井剖面图；
[0066] 图4是根据本发明实施例的测井相、响应频率、实测厚度交会图；
[0067] 图5是未经过本方法约束的沉积微相预测图；
[0068] 图6根据本发明实施例的方法得到的沉积微相预测图；
[0069] 图7是根据本发明实施例的经过归一化之后突出边滩有利相带的沉积微相预测 图；
[0070] 图8是根据本发明实施例的薄层沉积微相反演方法装置的结构框图。
【具体实施方式】
[0071] 为使本发明的目的、技术方案和优点更加清楚明白，下面结合实施方式和附图，对 本发明做进一步详细说明。在此，本发明的示意性实施方式及其说明用于解释本发明，但并 不作为对本发明的限定。
[0072] 如果将沉积微相（mf)作为一个模型空间，将测井数据（c)和地震数据（s)作为数 据空间，那么，对于正演而言，就有：
[0073] sc=Gmf
[0074] 其中，G表示泛函算子，G可以是一个积分算子、一个微分算子、一个矩阵或者一个 函数，也可以是一个过程。
[0075] 对于已知了测井数据（c)和地震数据（s)的情况，求取沉积微相（mf)的过程就是 反演为题，BP:
[0076] mf=G1 (s,c)
[0077] 其中，交会图可以是G1的具体表现形式，融合图是反演的具体实现方法。
[0078] 具体地，在本例中，提供了一种薄层沉积微相反演方法，如图1所示，包括以下步 骤：
[0079] 步骤101 :在工区中选择多口井，确定所述多口井中每口井的目的层的储层厚度；
[0080] 具体地，在工区中选择多口井之后，可以通过对每口井的测井解释，确定每口井的 储层厚度；
[0081] 例如，从工区中选出了 19 口井，按照对该工区的测井曲线和录井信息，对每口井 的目的层的储层厚度进行统计。例如，如果该区的目的层段以砂泥岩为主，则可以采用自然 电位（SP)曲线来确定每口井的目的层的储层厚度，具体地，首先确定SP的泥岩基线，对应 于大段泥岩层部分，此时偏离泥岩基线的异常部分可以认为是是储集层，异常幅度越大说 明泥质含量越少。在识别出渗透层之后，可以用"半幅点"的方法划分储集层的顶底，从而 确定储层的厚度，最后可以通过录井信息对储层厚度进行对比校正，从而最终确定出每口 井的目的层的储层厚度。
[0082] 步骤102:以所述多口井中每口井的目的层的储层底界作为分析界面，在地震剖 面中向上第一预定时间、向下第二预定时间作为分析时窗进行调谐分析，以确定每口井的 目的层的响应频率；
[0083] 在实际执行的过程中，可以先对多口井中每口井的目的层的储层底界进行层位标 定，标定后在全区（即全工区）此解释层位的层位追踪，并以此解释层位作为分析界面。
[0084] 具体地，可以以目的层的底界作为分析界面（即以确定的储层的下边界作为分析 界面），向上50ms，向下IOms作为分析时窗。所谓的向上向下指代的是在地震剖面中，因地 震剖面中是以时间来表示地层位置，因此在确定了分析界面后，可以将其映射到地震剖面 中的某个时间点，然后向上向下取预定的时间，就可以得到一个分析时窗。在确定了分析时 窗后，可以分别用短时傅里叶变换法和最大熵法进行调谐分析。
[0085] 在本例中，之所以用短时傅里叶变换法和最大熵法是考虑到短时傅里叶变化法相 对稳定，但是受到固定时窗的限制，频率分辨率低，而最大熵方法求取调谐频率时频率分辨 率高，但是容易带入算法误差。因此，发明人想到可以先用短时傅里叶变换法确定调谐频率 范围，然后应用最大熵法确定每口井在目的层的响应频率。如图2所示为基于短时傅里叶 变换的调谐频率连井剖面图，如图3所示为经过最大熵求取的调谐频率连井剖面图。
[0086] 具体地，短时傅里叶变换法的运算公式为：
[0088] 其中，s(t)表示地震数据，w(t_T)表示为窗函，f表示响应频率，t表示时间。
[0089] 最大熵法包括：
[0090] 先通过以下的K阶自回归差分方程模拟地震数据：
[0092] 其中，x(t)表示地震数据的地震道；(^表示自回归系数，且|ak|〈l，k= 1，2,. ..，e(t)表示预测误差，为一白噪序列，K= 2N/ln(2N)，其中，N表示地震记录的长度。
[0093] 然后，根据所述K阶自回归方程求出自回归系数app和方差后，通过以下公式计算 功率谱：
[0095]其中，P(f)表示估计功率谱，式表示e(t)的方差，表示自回归系数，f表示响 应频率，t表示时间。
[0096] 最后，从功率谱中确定出每口井的目的层的响应频率。
[0097] 因为短时傅里叶变换，算法稳定但频率分辨率低，通过此可以得到调谐频率范围， 然后采用最大熵法再对地震数据进行处理，最大熵法，算法频率分辨率高，但是容易引入算 法误差，在短时傅里叶变换得到的范围的约束下两法结合可以精确确定出每口井的目的层 的响应频率。
[0098] 步骤103 :沿所述多口井中每口井的目的层在单频体数据中提取最大振幅值，根 据所述最大振幅值在每口井的对应位置读取响应振幅；
[0099] 与响应频率不同，响应振幅的确定需要在单频体进行，可以应用小波变换法求取 单频体，对于任意函数f(t)GL2(R),其连续小波变换为：
[0101] 其中，W+表不信号时间-尺度分布，a表不尺度因子，b表不平移因子，也表不一 母小波，f表示响应频率，t表示时间，L2 (R)表示平方可积函数空间。
[0102] 在得到单频体后，沿目的层在单频体中提取最大振幅值，然后对应每口井的相应 位置读取响应振幅值。
[0103] 步骤104 :确定所述多口井中每口井在井点处的沉积微相类型，可以通过微相沉 积分析方法确定井点处的沉积微相类型；
[0104] 仍以上述的19 口井为例，在确定井点处的沉积微相的时候，一般需要考虑两个因 素，一是测井微相，二是地震剖面相。在进行沉积微相的确定的时候，可以先通过测井曲线 的钟形、箱型、锯齿型、直线型等曲线形态对应的沉积相态的不同，用测井微相标记地震数 据，然后再通过地震剖面相的识别确定井点处的沉积微相的类型。
[0105] 步骤105 :根据确定的沉积微相类型和确定的所述储层厚度，建立沉积微相与储 层厚度之间的关系；
[0106] 对于沉积微相与储层厚度的关系，可以引用沉积学中最经典的理论来证明，S卩，对 于陆相湖盆沉积而言，距离物源越近，储层越发育，厚度越大。从沉积相上来讲，辫状河沉积 环境相比曲流河沉积环境，辫状河相储层厚度更大，这一点充分证明了，储层厚度与沉积相 带之间存在密切的关系。实际应用中，在对地质背景充分认识的基础上，利用可以利用上述 步骤101和上述步骤104得到的结论，建立储层厚度与沉积微相之间的关系，在本例中，所 得到的关系为：边滩微相储层厚度最厚，天然堤相对较薄，河漫滩无储层。
[0107] 步骤106 :根据确定的储层厚度和确定的响应频率，建立储层厚度与响应频率之 间的关系；
[0108] 具体的，在本例中，整个研究区内储层厚度与响应频率之间的关系为：
[0109] y= -〇.5606x+34. 378
[0110] 对于边滩微相，储层厚度与响应频率的关系为：
[0111] y= -0.4043X+29. 913
[0112] 其中，X表示响应频率，y表示储层厚度。
[0113]