一种用于执行地球物理数据的并发已编码源反演以估计物理性能模型(41)的参数的方法,尤其适于没有固定接收器的采集几何结构的勘测,例如利用移动源和接收器的海洋地震勘测。在所述源上使用编码函数(32)以生成一个或更多并发已编码源数据集合(35)并对其模拟(34),所述编码函数(32)相对于互相关是正交或伪正交的。另外,接收器也被编码,其中接收器编码被设计为使得给定接收器在勘测(38)期间对其没有监听的源较不敏感。编码函数可以是中心频率、相位或两者互不相同的时间带通滤波器。通过将若干源分组成超级源、将对应的集合分组为超级集合然后应用上述编码策略可以进一步提高所述方法的效率。
【技术实现步骤摘要】
【国外来华专利技术】【专利摘要】一种用于执行地球物理数据的并发已编码源反演以估计物理性能模型(41)的参数的方法,尤其适于没有固定接收器的采集几何结构的勘测,例如利用移动源和接收器的海洋地震勘测。在所述源上使用编码函数(32)以生成一个或更多并发已编码源数据集合(35)并对其模拟(34),所述编码函数(32)相对于互相关是正交或伪正交的。另外,接收器也被编码,其中接收器编码被设计为使得给定接收器在勘测(38)期间对其没有监听的源较不敏感。编码函数可以是中心频率、相位或两者互不相同的时间带通滤波器。通过将若干源分组成超级源、将对应的集合分组为超级集合然后应用上述编码策略可以进一步提高所述方法的效率。【专利说明】正交源和接收器编码相关申请的交叉引用本申请要求2012年3月8日提交的名称为Orthogonal Source and ReceiverEncoding(正交源和接收器编码)的美国临时专利申请61/608,435和2012年11月14日提交的名称为Orthogonal Source and Receiver Encoding(正交源和接收器编码)的美国临时专利申请61/726,319的优先权,上述两个申请的全部内容通过引用方式并入本文。
本专利技术一般涉及地球物理勘探领域,并且更具体地,涉及地球物理数据处理。具体地,本专利技术是用于对从多个地球物理源诸如地震源获得的数据进行反演的方法,包括在模拟的一次执行中计算来自许多同时活动的地球物理源的数据的地球物理模拟。
技术介绍
即使应用现代计算能力,地震全波场反演仍是计算量很大的工作。但是,用此方法获得详细的地下表示的益处预期克服该阻碍。获得更快处理时间的算法和工作流程的开发是将此技术可适用于现场大规模数据的关键步骤。地震全波形反演包括数据的正向和伴随(adjoint)模拟的若干迭代。因此,降低正向和伴随计算运行的成本的技术将允许用户在合理时间量内解决较大规模的问题。 地球物理反演试图寻找一种最佳解释观测数据并满足地质和地球物理约束的地下属性模型。存在大量已知的地球物理反演方法。这些已知方法属于迭代反演和非迭代反演两种类别中的一种。以下定义是两种类别的每种通常表示的含义。 非迭代反演一该类反演通过假设一些简单背景模型并基于输入数据更新该模型实现。这种方法不将更新的模型用作反演的下一步骤的输入。对于地震数据的情形,这些方法通常被称为成像、移植、衍射、层析或博恩(Born)反演。 迭代反演一该类反演涉及地下属性模型的反复改进,从而找到圆满解释观测数据的模型。如果这种反演收敛,则最终的模型将更好地解释观测数据并且将更接近于实际的地下属性。迭代反演相对于非迭代反演通常产生更为准确的模型,但是计算费用更为昂贵。 迭代反演通常优于非迭代反演,因为迭代反演产生更准确的地下参数模型。遗憾的是,迭代反演在计算上是如此昂贵以致于将其应用到很多目标问题是不实际的。这种高额计算费用是因所有反演技术需要很多计算密集的模拟引起。任何单个模拟的计算时间与要反演的源的数量成比例,并且地球物理数据中通常有大量的源,其中以上用到的术语源指的是源设备的有效位置。在迭代反演中,这个问题恶化,因为在反演中必须计算的模拟数量与迭代的数量成比例,并且所需的迭代数量通常成百上千。 地球物理学中最常使用的迭代反演方法是成本函数优化。成本函数优化包括相对模型M的成本函数S (M)的值的迭代最小化或最大化,成本函数S(M)是计算数据和观测数据之间的不匹配性的度量(有时也被称为目标函数),其中计算数据通过计算机使用当前地球物理属性模型以及在给定地球物理属性模型表示的媒介中的源信号的物理控制传播进行模拟。模拟计算可以通过包括但不限于有限差分、有限元或射线跟踪的若干数值方法中的任意一种来完成。模拟计算能够在频域或时域中执行。 成本函数优化方法是局部或全局的。全局方法简单地包括针对模型群体(Ml,M2,M3,...}计算成本函数S(M)并从该群体中选择近似地最小化S(M)的一个或更多模型的组。如果期望进一步改善,这个新选择的模型组接着能够被用作产生新模型群的基础,这些新模型群能够再次针对成本函数S(M)被测试。对于全局方法,测试群中的每个模型能够被认为进行迭代,或者在更高等级,每个受测群组能够被认为进行迭代。已知的全局反演方法包括蒙特卡罗(MonteCarlo)、模拟的退火、遗传和进化算法。 遗憾的是,全局优化方法通常收敛非常慢,因此大多数地球物理反演基于局部成本函数优化。算法I概括了局部成本函数优化。 1.选择初始模型2.计算关于描述模型的参数的成本函数S(M)的梯度3.搜索更新模型,所述更新模型是初始模型在负梯度方向中的扰动,其更好地解释观测数据 算法I一用于执行局部成本函数优化的算法 通过将新更新的模型用作另一次梯度搜索的初始模型来迭代该过程。该过程继续直至找到圆满解释观测数据的更新模型。常用的局部成本函数反演方法包括梯度搜索、收敛梯度和牛顿法。 声学近似中的地震数据的局部成本函数优化是常见的地球物理反演任务,并且通常示例其它类型的地球物理反演。当在声学近似中反演地震数据时,成本函数能够被写为: Ng Mr Nt S(M)=YjYjYjW(ψ,.ιιΙ(.(Μ,r,/, νν?)- Wohs(r,?, νν^))(方程 I) g=\ r=\ t=l 其中: S =成本函数, M =描述地下模型的N个参数向量Oii1, m2,...mN), g =集合指数, Wg=集合g的源函数,其是空间坐标和时间的函数,对于点源,其是空间坐标的δ函数, Ng =集合数, r=集合内的接收器指数, Nr=集合中的接收器数, t =迹线(trace)内的时间样本下标, Nt=时间样本数, W =最小化标准函数(首选为W(x) = X2,其为最小二乘法(L2)标准), Ψε3?ε =根据模型M计算的地震压力数据 ΨΛβ =测量的地震压力数据 集合能够是能够在地震正向建模程序的一次运行中被模拟的任何类型的集合。通常情况下,集合对应于地震炮点,尽管所述炮点能够比点源更普遍。关于点源,集合指数g对应于各个点源的位置。对于平面波源,g将对应于不同平面波传播方向。这种广义源数据Ψ-能够现场得到或能够利用点源得到的数据合成。另一方面,计算数据ΨΜ1。通常能够在正向建模时利用广义源函数直接计算。对于许多类型的正向建模,包括有限差分建模,广义源所需的计算时间粗略地等于点源所需的计算时间。方程式(I)能够简化为: 【权利要求】1.一种用于对测量的地球物理数据进行迭代反演以确定地下区域的物理属性模型的计算机实现方法,其包括: 使用计算机对所述测量的地球物理数据的多个已编码集合求和,每个集合与单个源或一组源相关并被不同编码函数编码,其中所述编码函数选自相对于互相关是正交或伪正交的编码函数组,从而形成表示多个源的测量的地球物理数据的并发已编码集合;然后使用来自先前迭代的假设的物理属性模型或更新的物理属性模型模拟测量的地球物理数据的所述并发已编码集合;然后计算测量所述测量的地球物理数据的并发已编码集合与所述模拟的并发已编码集合之间的不匹配性的目标函数;接着优化所述目标函数以确定模型更本文档来自技高网...
【技术保护点】
一种用于对测量的地球物理数据进行迭代反演以确定地下区域的物理属性模型的计算机实现方法,其包括:使用计算机对所述测量的地球物理数据的多个已编码集合求和,每个集合与单个源或一组源相关并被不同编码函数编码,其中所述编码函数选自相对于互相关是正交或伪正交的编码函数组,从而形成表示多个源的测量的地球物理数据的并发已编码集合;然后使用来自先前迭代的假设的物理属性模型或更新的物理属性模型模拟测量的地球物理数据的所述并发已编码集合;然后计算测量所述测量的地球物理数据的并发已编码集合与所述模拟的并发已编码集合之间的不匹配性的目标函数;接着优化所述目标函数以确定模型更新;其中接收器被编码以使所述目标函数的计算对于给定接收器的所述多个源的一个或更多个较不敏感。
【技术特征摘要】
【国外来华专利技术】...
【专利技术属性】
技术研发人员:J·R·克雷布斯,Y·H·查,S·李,P·迪米特罗夫,A·A·穆勒,N·J·唐尼,P·S·如斯,
申请(专利权)人:埃克森美孚上游研究公司,
类型:发明
国别省市:美国;US
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。