专利名称：一种叠前地震定量成像方法
技术领域：
本发明关于地球物理地震成像，具体的讲是一种叠前地震定量成像方法。
背景技术：
弄清地下地质体的结构是地质、地球物理工作者的首要任务，地震成像是实现这一目标的主要手段。地震成像方法有两类射线偏移和波动方程偏移。射线偏移出现在上世纪70年代初，主要用于相对简单的地质体构造形态成像。波动方程偏移最初由J. F. Claerbout提出和实现，随后发展了一系列波动方程偏移方法，如有限差分偏移、频率-波数域偏移、克希霍夫积分偏移等。克希霍夫积分偏移是目前普遍使用的方法。波动方程偏移可用于构造非常复杂的地质体界面成像，同时能通过地震波振幅信息显示界面两侧地质体差异相对大小。但是，现有的地震成像方法只限于提供地质体构造形态，即界面形态成像，不能提供地质体内部的物理参数数据，即地质体的定量成像，不能直接提供更多反映地质体内部特征的数据，如地质体的岩性、孔隙、含油气性等。
发明内容
本发明实施例提供了一种叠前地震定量成像方法，提供一种高效、高分辨率的叠前地震地质体构造形态及其内部密度和体积模量参数的成像方法，为油气藏勘探提供定量的地球物理数据。本发明的目的之一，提供一种叠前地震定量成像方法，该方法具体步骤包括(1) 采集待成像区域的地震数据，根据地震数据获得测量的多炮炮集数据；( 根据测量的炮集数据建立包含密度和体积模量参数的当前模型；C3)测量待成像区域内一井孔处的密度和纵波速度，根据井孔处的密度和纵波速度获得井孔处的体积模量，根据井孔处的密度和体积模量建立先验模型；(4)对当前模型进行地震波场正演模拟生成正演地震波场数据， 对正演地震波场数据进行反付氏变换生成模拟炮集数据；( 根据模拟炮集数据和测量的炮集数据，获得残差数据；(6)计算测量的炮集数据的协方差矩阵和模型数据的协方差矩阵，根据残差数据、当前模型、先验模型、多炮炮集数据协方差矩阵、模型协方差矩阵计算目标函数值；(7)比较目标函数值与预设值的大小判断目标函数值大于或等于所述预设值时根据获得的残差数据，获得频率域剩余地震波场数据；根据正演地震波场数据和剩余地震波场数据获得密度和体积模量的共轭修改量；根据密度和体积模量的共轭修改量对当前模型进行迭代修改，然后进行步骤(4) 步骤(7)；判断目标函数值小于所述预设值时， 根据当前模型的密度和体积模量输出定量成像结果。本发明利用频率域单程声波方程延拓和全频带炮集数据，提供地下地质体定量成像，成像精度高，计算效率高。为让本发明的上述和其他目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附图式，作详细说明如下。
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。图1本发明叠前地震数据成像流程图；图2为当前模型的生成步骤；图3为待成像区域先验模型的生成步骤；图4采集数据处理后用于成像的炮集数据图；图5炮集数据频谱图；图6初始密度模型图；图7初始体积模量模型图；图8地层密度定量成像剖面图；图9地层体积模量定量成像剖面图。
具体实施例方式下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。本发明提供了一种叠前地震定量成像方法，该方法是一种基于频率域单程声波方程延拓的地震定量成像方法，从采集的地震炮集数据中获取地下构造形态和地质体密度及体积模量参数。如图1所示，为本发明叠前地震定量成像方法的流程图步骤S101，采集待成像区域的地震数据，根据所述的地震数据获得测量的多炮炮集数据；步骤S102，根据所述的测量的炮集数据建立包含密度和体积模量参数的当前模型；步骤S103，测量待成像区域内一井孔处的密度和纵波速度，根据所述的井孔处的密度和纵波速度获得所述井孔处的体积模量，根据所述井孔处的密度和体积模量建立先验模型；步骤S104，对所述的当前模型进行地震波场正演模拟生成正演地震波场数据，对所述的正演地震波场数据进行反付氏变换生成模拟炮集数据；步骤S105，根据所述的模拟炮集数据和测量的炮集数据，获得残差数据；步骤S106，计算测量的炮集数据的协方差矩阵和模型数据的协方差矩阵，根据所述的残差数据、当前模型、先验模型、多炮炮集数据协方差矩阵、模型协方差矩阵计算目标函数值；步骤S107，判断目标函数值是否小于预设值
判断是，执行步骤S111，根据当前模型的密度和体积模量输出定量成像结果；判断否，进行下述步骤步骤S108，根据残差数据获得频率域剩余地震波场数据；步骤S109，根据正演地震波场数据和剩余地震波场数据获得的密度参数和体积模量参数的共轭修改量；步骤S110，根据密度参数和体积模量参数的共轭修改量对当前模型进行迭代修改，然后重新执行步骤S104，步骤S105，步骤S106和步骤S107，直至目标函数值小于预设值时，输出成像结果。为更好理解本发明，下面结合上述步骤对本发明实施例做进一步详细描述。本发明是一种基于频率域单程声波方程延拓的叠前地震定量成像方法，从采集的地震炮集数据中获取地下构造形态和地质体密度及体积模量参数。本发明公开的叠前地震定量成像方法在本实施例中的具体步骤如下1.采集地震数据，对地震数据进行数据处理和频谱分析，获得炮集数据和频谱数据。在该步骤中的数据处理包括解编、坏炮、坏道剔除、静校正、地表一致性振幅补偿和叠前去除噪音。频谱分析采用常规地震数据频谱分析方法。2.如图2所示为当前模型的生成步骤，从获得的炮集数据中抽取共中心点道集数据，进行纵波速度分析，层纵波速度计算和时间-深度转换，获得包含深度域体积模量和密度参数的当前模型。其中，该步骤中纵波速度分析为常规地震纵波速度分析。上述的层纵波速度计算采用纵波速度分析获得的均方根纵波速度和如公式(1)所示Dix公式进行计算
权利要求
1.一种叠前地震定量成像方法，其特征在于，所述的方法包括(1)采集待成像区域的地震数据，根据所述的地震数据获得测量的多炮炮集数据；(2)根据所述的测量的炮集数据建立包含密度和体积模量参数的当前模型；(3)测量待成像区域内一井孔处的密度和纵波速度，根据所述的井孔处的密度和纵波速度获得所述井孔处的体积模量，根据所述井孔处的密度和体积模量建立待成像区域的先验模型；(4)对所述的当前模型进行地震波场正演模拟生成正演地震波场数据，对所述的正演地震波场数据进行反付氏变换生成模拟炮集数据；(5)根据所述的模拟炮集数据和测量的炮集数据，获得残差数据；(6)计算测量的炮集数据的协方差矩阵和当前模型协方差矩阵，根据所述的残差数据、 当前模型、先验模型、测量的炮集数据协方差矩阵、当前模型协方差矩阵计算目标函数值；(7)比较所述的目标函数值与预设值的大小判断所述的目标函数值大于或等于所述预设值时根据获得的残差数据，获得频率域剩余地震波场数据；根据所述的正演地震波场数据和剩余地震波场数据获得密度和体积模量的共轭修改量；根据所述的密度和体积模量的共轭修改量对所述的当前模型进行迭代修改，然后进行步骤(4) 步骤(7)；判断所述的目标函数值小于所述预设值时，根据当前模型的密度和体积模量输出定量成像结果。
2.如权利要求1所述的方法，其特征在于，所述的步骤(2)包括从获得的测量炮集数据中抽取共中心点道集数据；对所述的共中心点道集数据进行速度分析、层速度计算和时间-深度转换，获得包含深度域体积模量和密度的当前模型。
3.如权利要求1所述的方法，其特征在于，所述的步骤(3)包括采用声波测井和密度测井获得所述井孔处的纵波速度和密度，采用常规地震解释获得的层位数据对所述的井孔处的密度和体积模量进行控制，并对所述的井孔处的密度和体积模量进行内插和外推，获得待成像区域包含密度和体积模量参数的先验模型。
4.如权利要求1所述的方法，其特征在于，所述的步骤(4)包括采用频率域单程声波方程延拓进行所述的地震波场正演模拟。
5.如权利要求1所述的方法，其特征在于，所述的步骤(5)根据所述的模拟炮集数据和测量的炮集数据获得残差数据包括采用下式计算残差数据Ad = d。bs_dcal其中，Δ d为残差数据，d。bs为测量的炮集数据，dcal为模拟炮集数据。
6.如权利要求1所述的方法，其特征在于，所述的测量的炮集数据协方差矩阵由下式获得
7.如权利要求1所述的方法，其特征在于，所述的当前模型协方差矩阵由下式获得
8.如权利要求1所述的方法，其特征在于，所述的步骤(6)中采用下式计算目标函数
9.如权利要求1所述的方法，其特征在于，所述的步骤(7)中的根据获得的残差数据， 获得频率域残差数据的剩余地震波场数据包括将获得的残差数据进行付氏变换，将付氏变换后的残差数据做为源数据，采用正演模拟的方法，在频率域进行回传延拓，获得所述的频率域残差数据的剩余地震波场数据。
10.如权利要求1所述的方法，其特征在于，所述的步骤(7)中的根据所述的正演地震波场数据和剩余地震波场数据获得密度参数和体积模量参数的共轭修改量包括在频率域分别计算所述的正演地震波场数据和剩余地震波场数据的全频带数据对时间变量的一阶导数，将所述的两个一阶导数相乘，进行反付氏变换，按下式对激发炮数进行累加，获得体积模量的共轭修改量
11.如权利要求1所述的方法，其特征在于，所述的步骤(7)中的根据所述的密度参数和体积模量参数的共轭修改量对所述的当前模型进行迭代修改包括采用常规共轭梯度法对所述的当前模型进行迭代修改。
全文摘要
本发明公开了一种叠前地震定量成像方法，该方法包括采集待成像区域的地震数据，根据地震数据获得测量的多炮炮集数据，并建立包含密度和体积模量参数的当前模型；测量待成像区域内一井孔处的密度和纵波纵波速度，根据井孔处的密度和纵波纵波速度获得井孔处的体积模量，建立先验模型；对当前模型进行地震波场正演模拟生成正演地震波场数据，并进行反付氏变换生成模拟炮集数据；根据模拟炮集数据和测量的炮集数据，获得残差数据；分别计算测量的炮集数据和模型数据的协方差矩阵，计算目标函数值，判断目标函数值小于预设值时，输出成像结果，判断目标函数值大于或等于预设值时，进行迭代修改，直至目标函数值小于所述预设值时，输出成像结果。
文档编号G01V1/30GK102385066SQ20101027424
公开日2012年3月21日 申请日期2010年9月6日 优先权日2010年9月6日
发明者石玉梅 申请人:中国石油天然气股份有限公司