买专利卖专利找龙图腾，真高效！ 查专利查商标用IPTOP,全免费！专利年费监控用IP管家,真方便！

申请/专利权人：中国石油天然气集团有限公司;中国石油集团东方地球物理勘探有限责任公司

摘要：本发明公开了一种岩溶期古断裂刻画方法及装置，该方法包括：提取井旁道地震子波进行合成记录，实现岩溶储层顶界面和包含最大海泛面在内的所有上覆速度控制层的井震标定以及骨架剖面解释与加密解释；对井旁道地震子波和层位解释成果进行时深转换；对深度域地震井旁道地震子波进行沿层拉平；对拉平重构地震数据体进行成像增强处理；在处理后的地震数据上，对提取的敏感地震属性进行融合；在岩溶古断裂平面预测图上，将计算得到的断裂的平面延伸长度与最大垂直断距进行拟合；基于拟合结果对岩溶期古断裂进行分类与分级，并投影到岩溶期储层顶面古地貌图上，完成最终刻画。本发明可以实现岩溶断裂的快速识别与分级分类，提升了预测精度与效率。

主权项：1.一种岩溶期古断裂刻画方法，其特征在于，包括：提取井旁道地震子波，根据所述井旁道地震子波获得合成记录，基于合成记录进行岩溶储层顶界面和包含最大海泛面在内的所有上覆速度控制层的井震标定，以合成记录的过井地震剖面为基础，进行岩溶储层顶界面和包含最大海泛面在内的所有上覆速度控制层的骨架剖面解释与加密解释，获得层位解释成果；基于合成记录中的三维速度场，对井旁道地震子波和层位解释成果进行时深转换，获得深度域井旁道地震子波和深度域层位解释成果；以深度域层位解释成果中的最大海泛面解释层位为基准，对深度域地震井旁道地震子波进行沿层拉平，得到拉平重构地震数据体；对拉平重构地震数据体进行成像增强处理，获得处理后的地震数据；在处理后的地震数据上，提取沿岩溶储层顶界面的敏感地震属性，对敏感地震属性进行融合，获得岩溶古断裂平面预测图；在岩溶古断裂平面预测图上，计算每条岩溶期古断裂的平面延伸长度与最大垂直断距，将每条岩溶期古断裂的平面延伸长度与最大垂直断距进行拟合，获得拟合结果；基于拟合结果对岩溶期古断裂进行分类与分级，并将分类与分级后的岩溶期古断裂投影到岩溶期储层顶面古地貌图上，完成最终刻画；基于合成记录中的三维速度场，对井旁道地震子波和层位解释成果进行时深转换，获得深度域井旁道地震子波和深度域层位解释成果，包括：以合成记录标定的时深曲线为基础，以层位解释成果为约束条件进行井点速度横向插值，得到初始平均速度场；通过垂向井点标定平均速度及横向沿层平均速度计算井震绝对误差量；基于井震绝对误差量，通过反距离加权方式对初始平均速度场进行校正；利用校正后的初始平均速度场将井旁道地震子波和层位解释成果转换到深度域；对拉平重构地震数据体进行成像增强处理，获得处理后的地震数据，包括：采用平滑滤波器对拉平重构地震数据体进行边缘保持滤波，获得滤波后的地震数据；按照如下方式采用平滑滤波器对拉平重构地震数据体进行边缘保持滤波：拉平重构前的井旁道地震子波为Ds,t,l，其中，s、t和l分别表示平面属性的Inline、Crossline和Time三个方向的变量索引，平滑滤波采用沿着Inline方向进行，每次处理Inline方向的一个剖面fu,v，其中，u表示Crossline方向的变量索引下标，v表示Time方向的时间索引下标；假设剖面fu,v中的一个任意点u0,v0为中心，点u0,v0的函数值为fu0,v0，构造一个2×P+1×2×P+1的模板，其中，P为常数，并计算模板内的数据方差对模板中的每一个点ui,vi进行如下处理：若模板位置标量δui,vi＝1，否则δui,vi＝0；统计模板位置标量δ*中δui,vi＝1的位置的个数若则点u0,v0的函数值fu0,v0更新为否则点u0,v0的函数值fu0,v0更新为当点u0,v0的函数值fu0,v0更新之后，算法处理下一个点u′0,v′0的函数值fu′0,v′0；对敏感地震属性进行融合，包括：假设待融合平面地震属性记为：Ats,t，t＝1,2,K,N，那么属性At的边缘数据分析场定义为：其中s和t分别表示平面属性的Inline和Crossline两个方向的变量索引；▽A表示梯度，表示梯度的大小，表示边缘梯度的方向，表示边缘梯度的垂直方向；表示的转置矩阵；对于每一个点x＝x1,x2，令A*x1,x2＝[A1x1,x2,A2x1,x2,...,ANx1,x2]，其中Ai′x1,x2表示Ats,t在点x＝x1,x2处的函数值，其中，i′＝1,2,K,N；对于两个点x＝x1,x2和点x'＝x′1,x′2的差分为A*x-A*x′；当x与x′之间的距离趋近于无穷小时，微分算子表示为：其平方范数为：如果定义一个二阶矩阵其中i″,j″＝1,2，那么|dA*|2＝dxTGdx，其中dx＝dx1,dx2T；如果采用张量表示形式，那么表示为：矩阵的特征值λmaxmin和特征向量θmaxmin的计算方法为： 其中，对于多个地震属性融合问题，采用加权平均的方式得到张量矩阵： 其中，dicxi′表示x到原点的距离；最终的融合结果为：Vx＝λmax×θmax×signθmax×sn×▽A*，sign*表示符号函数。

全文数据：

权利要求：

百度查询： 中国石油天然气集团有限公司 中国石油集团东方地球物理勘探有限责任公司 岩溶期古断裂刻画方法及装置