一种基于柯西分布的稀疏脉冲反演方法
【专利摘要】一种基于柯西分布的稀疏脉冲反演方法，拾取得到层位数据，校正并内插和平滑；根据叠后地震数据、层位数据和已知的钻井分层数据，把深度域的声波时差曲线和密度曲线标定为时间域的曲线，生成波阻抗曲线，并在井旁地震道上提取地震子波；利用层位数据和时间域波阻抗曲线，生成初始波阻抗模型；对地震数据和子波数据分别进行归一化，反演得到道反射系数序列，推导得到波阻抗与反射系数的关系，最终得到反演的波阻抗剖面。本发明反演结果的准确性和稳定性好，能够正确反映岩层的物理特征，不易出现假象且运算效率高。
【专利说明】一种基于柯西分布的稀疏脉冲反演方法
【技术领域】
[0001] 本发明涉及地球物理勘探技术，属于储层预测反演技术类，是一种基于柯西分布 的稀疏脉冲反演方法。
【背景技术】
[0002] 地震勘探是通过人工激发地震波，在地表或地下用单分量或多分量传感器记录地 震波的地层响应，研究它们在地层中的传播规律，通过地震数据处理及反演等方法以查明 地下的地质构造岩性特征，进而寻找油气等矿产资源的地球物理勘探方法。地震勘探从认 识地下的构造形态开始，发展到直接应用地震信息判断岩性、分析岩相、定量计算岩层的物 性参数及寻找油气显示等。地震反演技术正是这一发展过程的产物。
[0003] 地震反演的基本目的是利用地震波在地下介质中的传播规律，通过数据采集、处 理与解释等流程，来推测地下岩层结构和物性参数的空间分布，为勘探开发提供重要依据。 在地震反演研究中，有多种参数的反演方法，如波阻抗、速度、密度、孔隙度、渗透率、泊松比 等。由于波阻抗信息是联系地质和地球物理的一座桥梁，在叠后计算数据量相对要小，在实 际生产中应用方便而且效果明显，因此波阻抗反演在地震反演中具有特殊的地位，地震反 演通常是指波阻抗反演。
[0004] 常规的地震波阻抗反演就是指利用地震贤料反演地层/岩层波阻抗的地震特殊 处理技术。与地震多参数模式识别预测储层油气、神经网络预测地层参数、振幅拟合预测储 层厚度等统计性方法相比，波阻抗反演具有明确的物理意义，它是储层岩性预测、油藏特征 描述的确定性方法，在实际应用中取得了显著的地质效果。
[0005] 波阻抗反演技术在20世纪70年代开始出现，当时对地震反演的研究只是以基于 褶积模型的叠后一维波阻抗反演为主，80年代得到了蓬勃发展。1983年，Cooke介绍了地震 资料的广义线性反演方法，从而揭开了波阻抗反演技术的新篇章。周竹生等人在90年代初 期提出了综合利用地质、地震和测井资料进行约束反演，克服单一的线性反演方法的缺陷。 90年代中期，李宏兵提出了递推反演与宽带约束反演结合的方法。该方法的推出，解决了从 单道出发的反演方法不能在根本上消除噪音的困惑。在此基础之上，有人进行了无井多道 反演和有井多道反演的研究，使波阻抗反演方法更加完善。同时，非线性反演方法也得到了 极大的发展，自从Rothmans利用模拟退火法成功地解决了地震勘探中大静校正量估算问 题之后，模拟退火法受到地球物理上作者们的重视，人们开始将其应用到波阻抗反演中。遗 传算法和模拟退火算法一样，利用转移概率来指导全局搜索过程。90年代至今，围绕一维波 阻抗反演的各类算法以及应用成果层出不穷，随着研究的升温，在1997年左右开始出现了 一些反思的文章，指出了波阻抗反演中存在的一些缺陷，并提出了一些解决方案。但是随着 这类方法的引进和在生产单位的推广应用，我们逐渐发现，这些方法在实用中很难获得好 的效果，很难获得可信的反演结果。利用这些反演结果进行的储层横向预测，在己知的井点 还比较可靠，离开井点一定距离后，往往与实际钻井结果相差甚远，预测出的薄层、薄互层 及储层岩性的可信度很低。
[0006] 目前的反演方法多数是以模型为基础的方法，这些方法一般都依据测井及地质资 料建立初始模型，通过广义线性反演方法进行迭代求取岩性参数。由于该问题的非线性，所 以除了要求精细的子波外，还要求初始模型接近真实模型，才能达到可靠的结果，即反演结 果强烈依赖于初始模型的选择。除此类方法外，全局优化的反演方法(如遗传算法和模拟退 火算法等）虽然克服了基于模型方法对初始模型依赖性强的缺陷，但是由于其得到全局最 优的反演结果，所以反演速度很慢。
【发明内容】
[0007] 本发明目的是提供一种运算效率非常高的基于柯西分布的叠后波阻抗反演方法。
[0008] 本发明通过如下技术方案实现：
[0009] 1)采用常规的地震勘探方法采集地震数据，对地震数据进行常规处理得到叠后地 震数据；
[0010] 2)对叠后地震数据进行层位拾取得到层位数据，对确定的目的层层位进行检验和 校正以及内插和平滑；
[0011] 3)采用常规的测井方法得到测井数据，得到测井声波时差曲线和密度曲线；
[0012] 4)根据叠后地震数据、层位数据和已知的钻井分层数据，把深度域的声波时差曲 线和密度曲线标定为时间域的曲线，同时生成井中的时间域波阻抗曲线数据，并在井旁地 震道上提取地震子波；
[0013] 所述的标定为利用测井曲线和地震子波模拟井旁地震记录，实现测井分层到地震 层位的标定和映射，由此得到时深关系曲线，由此时深关系可以将深度域的测井曲线转换 为时间域曲线。
[0014] 5)利用步骤2)的层位数据和步骤4)得到的时间域波阻抗曲线，生成初始波阻抗 模型；
[0015] 6)对地震数据和步骤4)得到的子波数据分别进行归一化，归一到范围[-1，1]之 间；
[0016] 7)把步骤1)读入的一道地震数据、步骤5)生成的初始波阻抗模型的一道数据和 步骤4)提取的子波数据，代入以下公式中，通过反演得到该道的反射系数序列：
[0017] r= (GTG+ μ Q+ P CTC)(GTd+ P CT ξ ) (1)
[0018] 其中（ΚΙ d2, . . .，dN]T是叠后地震数据，N为地震数据的总采样点数； r=[ri，r2,. . .，rN]T是反射系数序列；G是NXN维子波褶积矩阵，上角标T代表矩阵的转置； μ是稀疏约束因子，控制反射系数的稀疏程度；矩阵Q的对角元素为：
【权利要求】
1. 一种基于柯西分布的叠后波阻抗反演方法，特点是通过以下步骤实现： 1) 采用常规的地震勘探方法采集地震数据，对地震数据进行常规处理得到叠后地震数 据； 2) 对叠后地震数据进行层位拾取得到层位数据，对确定的目的层层位进行检验和校正 以及内插和平滑； 3) 采用常规的测井方法得到测井数据，得到测井声波时差曲线和密度曲线； 4) 根据叠后地震数据、层位数据和已知的钻井分层数据，把深度域的声波时差曲线和 密度曲线标定为时间域的曲线，同时生成井中的时间域波阻抗曲线，并在井旁地震道上提 取地震子波； 5) 利用步骤2)的层位数据和步骤4)得到的时间域波阻抗曲线，生成初始波阻抗模型； 6) 对地震数据和步骤4)得到的子波数据分别进行归一化，归一到范围[-1，1]之间； 7) 把步骤1)读入的一道地震数据、步骤5)生成的初始波阻抗模型的一道数据和步骤 4)提取的子波数据，代入以下公式中，通过反演得到该道的反射系数序列：
其中cKdi，d2,. . .，dN]T是叠后地震数据，N为地震数据的总采样点数； r=[ri，r2,. . .，rN]T是反射系数序列；G是NXN维子波褶积矩阵，上角标T代表矩阵的转置； U是稀疏约束因子，控制反射系数的稀疏程度；矩阵Q的对角元素为：
Qm表示矩阵Q的第n行第n列元素的值，n为Q矩阵的行列号，矩阵Q除对角元素外其余元 素均为零，％代表噪声的标准方差，Rn是由初始波阻抗模型计算得到的第n个采样点位置 处的初始反射系数；P为模型约束因子，控制反演结果对初始模型的依赖程度；C为积分算 子矩阵，其离散形式表示为：
公式（1)中上角标-1为对矩阵求逆
^第n个采样点处的相对波阻 抗值，lei为反演时窗内第一个米样点对应的初始波阻抗值，1"为反演时窗内第n个米样点 处的初始波阻抗值，In为自然对数符号，ri为反演时窗内第i个采样点的反射系数值， 表示对从第1个采样点到第n个采样点的ri进行求和计算；
8) 通过步骤7)中第n个采样点处的相对波阻抗
的定义，通过推导得 到波阻抗与反射系数的关系：
其中：1"为反演时窗内第n个米样点的波阻抗值，]；(!为反演时窗内第一个米样点对应 的初始波阻抗值，^为步骤7)中反演得到的第i个采样点反射系数值，e表示自然对数的 底，
表示对从第1个采样点到第n个采样点的ri进行求和计算； 通过公式（2)将步骤7)反演得到的一道的反射系数系列转换为波阻抗序列，得到该道 的波阻抗反演结果； 9)对所有的地震道重复步骤7)至8)过程，得到所有道的波阻抗反演结果。
2.根据权利要求1的方法，特点是步骤4)所述的标定为利用测井曲线和地震子波模拟 井旁地震记录，实现测井分层到地震层位的标定和映射，由此得到时深关系曲线，由此时深 关系可以将深度域的测井曲线转换为时间域曲线。
【文档编号】G01V1/28GK104280767SQ201310292088
【公开日】2015年1月14日 申请日期:2013年7月12日 优先权日:2013年7月12日
【发明者】赵剑, 万忠宏, 王增波, 白雪莲, 陈亚军, 陈茂山, 贾继军, 崔京彬 申请人:中国石油天然气集团公司, 中国石油集团东方地球物理勘探有限责任公司