一种基于柯西分布的叠后波阻抗反演方法技术

一种基于柯西分布的叠后波阻抗反演方法，采用常规的地震勘探方法采集地震数据对其常规处理得到叠后地震数据；对叠后地震数据进行层位拾取得到层位数据，对确定的目的层层位进行检验和校正以及内插和平滑；采用常规的测井方法得到测井声波时差曲线和密度曲线；根据叠后地震数据、层位数据和已知的钻井分层数据，把深度域的声波时差曲线和密度曲线标定为时间域的曲线，同时生成井中的时间域波阻抗曲线数据，并提取地震子波；生成初始波阻抗模型；对子波数据分别进行归一化；通过反演得到某道的反射系数序列；通过第n个采样点处的相对波阻抗的定义，通过推导得到波阻抗与反射系数的关系对所有的地震道重复上述过程，得到所有道的波阻抗反演结果。

【技术实现步骤摘要】
【国外来华专利技术】
本专利技术涉及地球物理勘探技术，属于储层预测反演技术类，是一种基于柯西分布 的叠后波阻抗反演方法。
技术介绍
地震勘探是通过人工激发地震波，在地表或地下用单分量或多分量传感器记录地 震波的地层响应，研究它们在地层中的传播规律，通过地震数据处理及反演等方法以查明 地下的地质构造岩性特征，进而寻找油气等矿产资源的地球物理勘探方法。地震勘探从认 识地下的构造形态开始，发展到直接应用地震信息判断岩性、分析岩相、定量计算岩层的物 性参数及寻找油气显示等。地震反演技术正是这一发展过程的产物。 地震反演的基本目的是利用地震波在地下介质中的传播规律，通过数据采集、处 理与解释等流程，来推测地下岩层结构和物性参数的空间分布，为勘探开发提供重要依据。 在地震反演研究中，有多种参数的反演方法，如波阻抗、速度、密度、孔隙度、渗透率、泊松比 等。由于波阻抗信息是联系地质和地球物理的一座桥梁，在叠后计算数据量相对要小，在实 际生产中应用方便而且效果明显，因此波阻抗反演在地震反演中具有特殊的地位，地震反 演通常是指波阻抗反演。 常规的地震波阻抗反演就是指利用地震贤料反演地层/岩层波阻抗的地震特殊 处理技术。与地震多参数模式识别预测储层油气、神经网络预测地层参数、振幅拟合预测储 层厚度等统计性方法相比，波阻抗反演具有明确的物理意义，它是储层岩性预测、油藏特征 描述的确定性方法，在实际应用中取得了显著的地质效果。 目前的反演方法多数是以模型为基础的方法，这些方法一般都依据测井及地质资 料建立初始模型，通过广义线性反演方法进行迭代求取岩性参数。由于该问题的非线性，所 以除了要求精细的子波外，还要求初始模型接近真实模型，才能达到可靠的结果，即反演结 果强烈依赖于初始模型的选择。除此类方法外，全局优化的反演方法（如遗传算法和模拟 退火算法等）虽然克服了基于模型方法对初始模型依赖性强的缺陷，但是由于其得到全局 最优的反演结果，所以反演速度很慢。
技术实现思路
本专利技术目的是提供。 本专利技术通过如下技术方案实现： 1)采用常规的地震勘探方法采集地震数据，对地震数据进行常规处理得到叠后地 震数据； 2)对叠后地震数据进行层位拾取得到层位数据，对确定的目的层层位进行检验和 校正以及内插和平滑； 3)采用常规的测井方法得到测井数据，得到测井声波时差曲线和密度曲线； 4)根据叠后地震数据、层位数据和已知的钻井分层数据，把深度域的声波时差曲 线和密度曲线标定为时间域的曲线，同时生成井中的时间域波阻抗曲线数据，并在井旁地 震道上提取地震子波； 所述的标定为利用测井曲线和地震子波模拟井旁地震记录，实现测井分层到地震 层位的标定和映射，由此得到时深关系曲线，由此时深关系可以将深度域的测井曲线转换 为时间域曲线。 5)利用步骤2)的层位数据和步骤4)得到的时间域波阻抗曲线，生成初始波阻抗 模型； 6)对地震数据和步骤4)得到的子波数据分别进行归一化，归一到范围之 间； 7)把步骤1)读入的一道地震数据、步骤5)生成的初始波阻抗模型的一道数据和 步骤4)提取的子波数据，代入以下公式中，通过反演得到该道的反射系数序列： r = (GtG+ μ Q+ p CtC) (GTd+ P Ct ξ ) (I) 其中d =T是叠后地震数据，N为地震数据的总采样点数；r = T是反射系数序列；G是NXN维子波褶积矩阵，上角标T代表矩阵的转置；μ 是稀疏约束因子，控制反射系数的稀疏程度；矩阵Q的对角元素为【主权项】1. ，步骤如下： 1) 采用常规的地震勘探方法采集地震数据，对地震数据进行常规处理得到叠后地震数 据； 2) 对叠后地震数据进行层位拾取得到层位数据，对确定的目的层层位进行检验和校正 以及内插和平滑； 3) 采用常规的测井方法得到测井数据，得到测井声波时差曲线和密度曲线； 4) 根据叠后地震数据、层位数据和已知的钻井分层数据，把深度域的声波时差曲线和 密度曲线标定为时间域的曲线，同时生成井中的时间域波阻抗曲线数据，并在井旁地震道 上提取地震子波； 步骤4)所述的标定为利用测井曲线和地震子波模拟井旁地震记录，实现测井分层到 地震层位的标定和映射，由此得到时深关系曲线，由此时深关系可以将深度域的测井曲线 转换为时间域曲线。 5) 利用步骤2)的层位数据和步骤4)得到的时间域波阻抗曲线，生成初始波阻抗模 型； 6) 对地震数据和步骤4)得到的子波数据分别进行归一化，归一到范围之间； 7) 把步骤1)读入的一道地震数据、步骤5)生成的初始波阻抗模型的一道数据和步骤 4)提取的子波数据，代入以下公式中，通过反演得到该道的反射系数序列： r = (GtG+ μ Q+ p CtC) (GTd+ P Ct ξ ) (I) 其中d = 是叠后地震数据，N为地震数据的总采样点数；r = T是反射系数序列；G是NXN维子波褶积矩阵，上角标T代表矩阵的转置；μ 是稀疏约束因子，控制反射系数的稀疏程度；矩阵Q的对角元素为，Qm 表示矩阵Q的第η行第η列元素的值，η为Q矩阵的行列号，矩阵Q除对角元素外其余元素 均为零，％代表噪声的标准方差，Rn是由初始波阻抗模型计算得到的第η个采样点位置处 的初始反射系数；P为模型约束因子，控制反演结果对初始模型的依赖程度；C为积分算子 矩阵，其离散形式表示为：公式（1)中上角标-1为对矩阵求逆为第η个采样点处的相对波阻 抗值，Ιο为反演时窗内第一个米样点对应的初始波阻抗值，In为反演时窗内第η个米样点 处的初始波阻抗值，In为自然对数符号，ri为反演时窗内第i个采样点的反射系数值， 表示对从第1个采样点到第η个采样点的ri进行求和计算； 8) 通过步骤7)中第η个采样点处的相对波阻抗的定义，通过推导得 到波阻抗与反射系数的关系：(2) 其中：1"为反演时窗内第η个米样点的波阻抗值，Itl为反演时窗内第一个米样点对应 的初始波阻抗值，^为步骤7)中反演得到的第i个采样点反射系数值，e表示自然对数的 底表示对从第1个采样点到第η个采样点的ri进行求和计算； 通过公式（2)将步骤7)反演得到的一道的反射系数系列转换为波阻抗序列，便可得到 该道的波阻抗反演结果； 9)对所有的地震道重复步骤7)至8)过程，得到所有道的波阻抗反演结果。【专利摘要】，采用常规的地震勘探方法采集地震数据对其常规处理得到叠后地震数据；对叠后地震数据进行层位拾取得到层位数据，对确定的目的层层位进行检验和校正以及内插和平滑；采用常规的测井方法得到测井声波时差曲线和密度曲线；根据叠后地震数据、层位数据和已知的钻井分层数据，把深度域的声波时差曲线和密度曲线标定为时间域的曲线，同时生成井中的时间域波阻抗曲线数据，并提取地震子波；生成初始波阻抗模型；对子波数据分别进行归一化；通过反演得到某道的反射系数序列；通过第n个采样点处的相对波阻抗的定义，通过推导得到波阻抗与反射系数的关系对所有的地震道重复上述过程，得到所有道的波阻抗反演结果。【IPC分类】G01V1-28, G01V1-36【公开号】CN104769458【申请号】CN201480002799【专利技术人】杨顺伟 【申请人】杨顺伟【公开日】2015年7月8日【申请日】2014年7月15日本文档来自技高网...

【技术保护点】
一种基于柯西分布的叠后波阻抗反演方法，步骤如下：1)采用常规的地震勘探方法采集地震数据，对地震数据进行常规处理得到叠后地震数据；2)对叠后地震数据进行层位拾取得到层位数据，对确定的目的层层位进行检验和校正以及内插和平滑；3)采用常规的测井方法得到测井数据，得到测井声波时差曲线和密度曲线；4)根据叠后地震数据、层位数据和已知的钻井分层数据，把深度域的声波时差曲线和密度曲线标定为时间域的曲线，同时生成井中的时间域波阻抗曲线数据，并在井旁地震道上提取地震子波；步骤4)所述的标定为利用测井曲线和地震子波模拟井旁地震记录，实现测井分层到地震层位的标定和映射，由此得到时深关系曲线，由此时深关系可以将深度域的测井曲线转换为时间域曲线。5)利用步骤2)的层位数据和步骤4)得到的时间域波阻抗曲线，生成初始波阻抗模型；6)对地震数据和步骤4)得到的子波数据分别进行归一化，归一到范围[‑1,1]之间；7)把步骤1)读入的一道地震数据、步骤5)生成的初始波阻抗模型的一道数据和步骤4)提取的子波数据，代入以下公式中，通过反演得到该道的反射系数序列：r＝(GTG+μQ+ρCTC)‑1(GTd+ρCTξ)    (1)其中d＝[d1,d2,...,dN]T是叠后地震数据，N为地震数据的总采样点数；r＝[r1,r2,...,rN]T是反射系数序列；G是N×N维子波褶积矩阵，上角标T代表矩阵的转置；μ是稀疏约束因子，控制反射系数的稀疏程度；矩阵Q的对角元素为：Qnn表示矩阵Q的第n行第n列元素的值，n为Q矩阵的行列号，矩阵Q除对角元素外其余元素均为零，σr代表噪声的标准方差，Rn是由初始波阻抗模型计算得到的第n个采样点位置处的初始反射系数；ρ为模型约束因子，控制反演结果对初始模型的依赖程度；C为积分算子矩阵，其离散形式表示为：公式(1)中上角标‑1为对矩阵求逆；为第n个采样点处的相对波阻抗值，I0为反演时窗内第一个采样点对应的初始波阻抗值，In为反演时窗内第n个采样点处的初始波阻抗值，ln为自然对数符号，ri为反演时窗内第i个采样点的反射系数值，表示对从第1个采样点到第n个采样点的ri进行求和计算；8)通过步骤7)中第n个采样点处的相对波阻抗的定义，通过推导得到波阻抗与反射系数的关系：In=I0e2Σi=1nri---(2)]]>其中：In为反演时窗内第n个采样点的波阻抗值，I0为反演时窗内第一个采样点对应的初始波阻抗值，ri为步骤7)中反演得到的第i个采样点反射系数值，e表示自然对数的底，表示对从第1个采样点到第n个采样点的ri进行求和计算；通过公式(2)将步骤7)反演得到的一道的反射系数系列转换为波阻抗序列，便可得到该道的波阻抗反演结果；9)对所有的地震道重复步骤7)至8)过程，得到所有道的波阻抗反演结果。...

【技术特征摘要】
【国外来华专利技术】

【专利技术属性】
技术研发人员：杨顺伟，
申请(专利权)人：杨顺伟，
类型：发明
国别省市：广东;44