CN109143346A - 叠前混合非线性反演方法及计算机可读存储介质 - Google Patents

叠前混合非线性反演方法及计算机可读存储介质 Download PDF

Info

Publication number
CN109143346A
CN109143346A CN201710464861.9A CN201710464861A CN109143346A CN 109143346 A CN109143346 A CN 109143346A CN 201710464861 A CN201710464861 A CN 201710464861A CN 109143346 A CN109143346 A CN 109143346A
Authority
CN
China
Prior art keywords
transverse
longitudinal wave
prestack
wave reflection
reflection coefficient
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201710464861.9A
Other languages
English (en)
Other versions
CN109143346B (zh
Inventor
张远银
刘喜武
刘宇巍
霍志周
刘志远
刘炯
钱恪然
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Chemical Corp
Sinopec Exploration and Production Research Institute
Original Assignee
China Petroleum and Chemical Corp
Sinopec Exploration and Production Research Institute
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Petroleum and Chemical Corp, Sinopec Exploration and Production Research Institute filed Critical China Petroleum and Chemical Corp
Priority to CN201710464861.9A priority Critical patent/CN109143346B/zh
Publication of CN109143346A publication Critical patent/CN109143346A/zh
Application granted granted Critical
Publication of CN109143346B publication Critical patent/CN109143346B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30—Analysis
    • G01V1/303—Analysis for determining velocity profiles or travel times
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30—Analysis
    • G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/62—Physical property of subsurface
    • G01V2210/622—Velocity, density or impedance
    • G01V2210/6222—Velocity; travel time
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/63—Seismic attributes, e.g. amplitude, polarity, instant phase

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

公开了一种叠前混合非线性反演方法及计算机可读存储介质。该方法可以包括：对于叠前地震数据进行叠前线性反演，获得横波反射系数的初始值与横纵波速度比的初始值；基于横波反射系数、横纵波速度比、密度梯度，构造组合弹性参量；将叠前纵波反射系数公式表示为关于密度梯度、纵波反射系数和组合弹性参量的线性方程，对叠前纵波反射系数公式进行线性求解，获得密度梯度、纵波反射系数和组合弹性参量；构造关于横波反射系数和横纵波速度比的目标函数，进行非线性求解，进而获得最终的横波反射系数与最终的横纵波速度比。本发明将传统四参数非线性问题降解为两参数非线性问题，实现高效率、高精度地求取各种弹性参数。

Description

叠前混合非线性反演方法及计算机可读存储介质
技术领域
本发明涉及油气地球物理技术领域，更具体地，涉及一种叠前混合非线性反演方法。
背景技术
自1919年描述各向同性介质界面下纵波入射时反射和透射系数的Zoeppritz方程诞生以来，上世纪60年代以后，许多学者依据实际勘探地震需要将其纵波入射时的纵波反射(P-P)和横波反射(P-SV)公式近似成多种解析形式，目前大约已有三十多种，形成了P-P波反演、P-SV波反演和P-P与P-SV波联合反演等三类反演方法，被广泛地应用于实际储层和流体预测之中。
叠前反演通常是指AVO反演，其狭义理论基础为Zoeppritz方程P-P部分或其近似式和褶积模型。时域褶积模型理论基本观点是地下某点反射振幅是其反射系数与地震子波的褶积，而反射系数是联系地下弹性参数的桥梁，反演即据振幅等参数反演地下的弹性参数信息。依据平面简谐波传播理论，当纵波(P波)传至弹性界面时根据应力与应变的关系和连续条件，在满足Snell定律的条件下可以获得不同弹性介质的反演和透射系数。尽管叠前反演是利用多次覆盖资料求取弹性信息的必经之路，但实际资料的叠前反演却是严重的不适定问题，具有很强的多解性。不同方法或流程对于同一资料反演的结果常差别甚远，而同一方法与流程反演的不同弹性参数精度也很不相同，难于满足勘探、开发与生产的需求。
由于Zoeppritz方程完整形式过于复杂且物理意义不太明确，诸多学者从区分固液相态与突出泊松比、体现速度和密度相对变化、幂级数或射线参数、突出弹性模量等四大类进行了约十三种形式的近似，比如最为典型的Aki&Richards近似等，从而实现了弹性参数在实际资料情况下的反演。依据经典的近似关系，叠前的反演方程为非线性方程组，主要包括纵波反射系数、横波反射系数、密度反射梯度和横纵波速度比等四个未知量。忽略地震子波、地质模型、噪音干扰等因素，依据对弹性未知量的求解策略不同，现有的叠前反演方法大致可以分为叠前线性反演和叠前非线性反演两种(Zhang et al.,2013)。
叠前线性反演方法常常假设横纵波速度比为已知的常数，从而将非线性方程组线性化求解，直接求解包含纵波反射系数、横波反射系数和密度梯度的线性近似方程，这种方法操作简便，运算量小，是工业界最为常用的方法，然而由于其直接线性化的近似导致非线性量的反演精度较低，且其结果强烈依赖于初始模型和子波估算质量(Sun,1999)。非线性反演方法则直接采用非线性寻优方案对纵波反射系数、横波反射系数、密度梯度和横纵波速度比进行直接的四参量非线性求解。叠前非线性的求解算法大致可以分为从目标函数出发(比如最速下降法，牛顿法，共轭梯度法，非线性规划法，非线性最小二乘法等)和从自变量出发(比如模拟退火法，原子跃迁法，量子退火法，量子遗传法，人工神经元法，遗传算法，蚁群算法，粒子群算法，免疫算法等)两类(张远银，2015)。得益于计算机的快速进步，非线性算法可以直接被用于求解地球物理问题，然而，叠前非线性反演方法的运算量巨大，且实际地质问题常常存在多个目标函数极值从而导致运算陷入局部寻优(Zhang et al.,2013)。因此，有必要开发一种具有较高的反演精度和反演速度的叠前混合非线性反演方法。
公开于本发明背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。
发明内容
本发明提出了一种叠前混合非线性反演方法，其能够将传统四参数非线性问题降解为两参数非线性问题，实现高效率、高精度地求取各种弹性参数，提高了运算效率与精度。
根据本发明的一方面，提出了一种叠前混合非线性反演方法。所述方法可以包括：对于叠前地震数据进行叠前线性反演，获得横波反射系数Rs的初始值与横纵波速度比Rk的初始值；基于所述横波反射系数Rs、所述横纵波速度比Rk、密度梯度RD，构造组合弹性参量KS；将叠前纵波反射系数公式表示为关于所述密度梯度RD、纵波反射系数RP和组合弹性参量KS的线性方程，对所述叠前纵波反射系数公式进行线性求解，获得所述密度梯度RD、纵波反射系数RP和所述组合弹性参量KS；构造关于所述横波反射系数Rs和横纵波速度比Rk的目标函数，基于所述横波反射系数Rs的初始值与横纵波速度比Rk的初始值对所述目标函数进行非线性求解，进而获得最终的横波反射系数与最终的横纵波速度比。
优选地，所述叠前纵波反射系数公式为：
其中，Rpp(θ)表示叠前纵波反射系数，θ表示反射角与入射角的平均值。
优选地，所述组合弹性参量KS为：
优选地，所述目标函数为：
其中，F(Rs,Rk)表示目标函数。
优选地，所述目标函数为：
其中，F(Rs,Rk)表示目标函数。
优选地，通过标准粒子群算法对所述目标函数进行非线性求解。
根据本发明的另一方面，提出了一种计算机可读存储介质，其上存储有计算机程序，其中，所述程序被处理器执行时实现以下步骤：对于叠前地震数据进行叠前线性反演，获得横波反射系数Rs的初始值与横纵波速度比Rk的初始值；基于所述横波反射系数Rs、所述横纵波速度比Rk、密度梯度RD，构造组合弹性参量KS；将叠前纵波反射系数公式表示为关于所述密度梯度RD、纵波反射系数RP和组合弹性参量KS的线性方程，对所述叠前纵波反射系数公式进行线性求解，获得所述密度梯度RD、纵波反射系数RP和所述组合弹性参量KS；构造关于所述横波反射系数Rs和横纵波速度比Rk的目标函数，基于所述横波反射系数Rs的初始值与横纵波速度比Rk的初始值对所述目标函数进行非线性求解，进而获得最终的横波反射系数与最终的横纵波速度比。
优选地，所述叠前纵波反射系数公式为：
其中，Rpp(θ)表示叠前纵波反射系数，θ表示反射角与入射角的平均值。
优选地，所述组合弹性参量KS为：
优选地，所述目标函数为：
其中，F(Rs,Rk)表示目标函数。
本发明的有益效果在于：将反演方程中的线性量与非线性量分离，采用线性反演的方法快速准确反演线性量，降低方程的非线性化程度；采用粒子群非线性反演方法反演非线性量，利用线性反演的结果约束或指导非线性反演过程，提高反演的精度与速度，提高弹性参数反演精度与效率。
本发明的方法具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的，或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述，这些附图和具体实施方式共同用于解释本发明的特定原理。
附图说明
通过结合附图对本发明示例性实施例进行更详细的描述，本发明的上述以及其它目的、特征和优势将变得更加明显，其中，在本发明示例性实施例中，相同的参考标号通常代表相同部件。
图1示出了根据本发明的叠前混合非线性反演方法的步骤的流程图。
图2a、图2b分别示出了根据本发明的一个实施例的目标函数的误差值与部分放大的目标函数的误差值的示意图。
图3a、图3b分别示出了根据本发明的一个实施例的横波反射系数的误差值与横纵波速度比的误差值的示意图。
图4a、图4b、图4c、图4d分别示出了根据本发明的一个实施例的纵波反射系数、横波反射系数、密度梯度、横纵波速度比与理论值的对比示意图。
图5a、图5b、图5c分别示出了根据叠前线性反演、叠前非线性反演与本发明的一个实施例的纵波反射系数的示意图，图5d示出了三者的误差对比示意图。
图6a、图6b、图6c分别示出了根据叠前线性反演、叠前非线性反演与本发明的一个实施例的横波反射系数的示意图，图6d示出了三者的误差对比示意图。
图7a、图7b、图7c分别示出了根据叠前线性反演、叠前非线性反演与本发明的一个实施例的密度梯度的示意图，图7d示出了三者的误差对比示意图。
图8a、图8b、图8c分别示出了根据叠前线性反演、叠前非线性反演与本发明的一个实施例的横纵波速度比的示意图，图8d示出了三者的误差对比示意图。
图9示出了根据叠前线性反演、叠前非线性反演与本发明的一个实施例的计算时间的对比示意图。
具体实施方式
下面将参照附图更详细地描述本发明。虽然附图中显示了本发明的优选实施例，然而应该理解，可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了使本发明更加透彻和完整，并且能够将本发明的范围完整地传达给本领域的技术人员。
图1示出了根据本发明的叠前混合非线性反演方法的步骤的流程图。
根据本发明的叠前混合非线性反演方法可以包括：
步骤101，对于叠前地震数据进行叠前线性反演，获得横波反射系数Rs的初始值与横纵波速度比Rk的初始值。
具体地，叠前反演的基本原理如下，其基于Gildow(1987)对Aki&Richards公式(1980)的重新整理与简化：
其中，Rpp(θ)表示某一界面的叠前纵波反射系数，θ表示反射角与入射角的平均值，表示纵波反射系数，表示横波反射系数，表示密度梯度，Rk＝β/α表示横纵波速度比，α表示界面上下的平均纵波速度，β表示界面上下的平均横波速度，ρ表示界面上下的平均密度，Δα表示下界面与上界面的纵波速度差，Δβ表示下界面与上界面的横波速度差，Δρ表示下界面与上界面的密度差。
在非线性方程组的实际计算过程中，公式中的Rp，Rs，Rk和RD均为未知数据，Rpp(θ)和θ为已知数据。针对公式(1)进行线性反演则将Rk与Rs直接以常数化近似，从而可以把非线性方程组化简为线性方程组求解，即可根据公式(1)通过线性反演方法求取横波反射系数Rs的初始值与横纵波速度比Rk的初始值。
步骤102，基于横波反射系数Rs、横纵波速度比Rk、密度梯度RD，构造组合弹性参量KS；在一个示例中，组合弹性参量KS为：
步骤103，将叠前纵波反射系数公式表示为关于密度梯度RD、纵波反射系数RP和组合弹性参量KS的线性方程，对叠前纵波反射系数公式进行线性求解，获得密度梯度RD、纵波反射系数RP和组合弹性参量KS。
实际上，非线性方程组的求解过程并不一定绝对地线性或者非线性，其往往可以是多个线性量与分线性量的组合，这也意味着叠前反演可以依照未知量的属性特征分解为线性和非线性求解的组合，从而提高运算效率。将公式(1)的线性量与非线性量进行分离，可得到叠前纵波反射系数公式为：
其中，Rpp(θ)表示叠前纵波反射系数。
具体地，将横波反射系数Rs、横纵波速度比Rk、密度梯度RD代入公式(2)，构造组合弹性参量KS。将公式(1)表示为关于密度梯度RD、纵波反射系数RP和组合弹性参量KS的线性方程即为公式(3)，并进行线性求解，获得密度梯度RD、纵波反射系数RP和组合弹性参量KS。
步骤104，构造关于横波反射系数Rs和横纵波速度比Rk的目标函数，基于横波反射系数Rs的初始值与横纵波速度比Rk的初始值对目标函数进行非线性求解，进而获得最终的横波反射系数与最终的横纵波速度比。
在一个示例中，目标函数为：
其中，F(Rs,Rk)表示目标函数。
在一个示例中，目标函数为：
其中，F(Rs,Rk)表示目标函数。
具体地，构造关于横波反射系数Rs和横纵波速度比Rk的目标函数，其中目标函数可以为公式(4)或公式(5)，基于横波反射系数Rs的初始值与横纵波速度比Rk的初始值对目标函数进行非线性求解，获得最终的横波反射系数与最终的横纵波速度比。
在一个示例中，通过标准粒子群算法对目标函数进行非线性求解。
具体地，采用的非线性算法为标准粒子群算法，该算法是对鸟类群体的防御、捕食等行为仿真，将问题的搜索空间类比于鸟类的飞行空间，算法中每个优化问题的可行解都是搜索空间中的一只鸟，称之为“粒子”。每个粒子的运动可以用上述几条规则来描述，根据自己和其它粒子的“飞行经验”寻优，从而达到全空间搜索最优解的目的。假设一个n维搜索空间，每一维空间中总的粒子数为m，每个粒子所处的位置代表目标函数的一个潜在解，它们的位置量xi构成一个种群X＝(x1,x2,…xm)，每个粒子在空间中移动速度为vi。第i粒子目前为止搜索到的最优位置为pibest＝(pi1,pi2,…pim)，整个种群目前为止搜索到的最优位置为gbest＝(g1,g2,…gn)，则粒子的速度和位置可以按如下公式进行变化：
vid(t+1)＝wvid(t)+c1r1(pid(t)-xid(t))+c2r2(gd(t)-xid(t)) (6)
xid(t+1)＝xid(t)+vid(t+1)1≤i≤m 1≤d≤n (7)
其中，w表示惯性权重因子，cl为调节粒子飞向自身最好位置方向的步长，c2为调节粒子飞向全局最好位置的步长，r1,r2为[0，1]内的随机数，t为当前迭代代数。为避免由于速度过大而使粒子飞过了最优解的位置，有必要对速度的最大值加以限制，限定速度上限值为vmax，当vid>vmax时，令vid＝vmax；当vid<-vmax时，令vid＝-vmax。粒子群初始位置和初始速度随机产生，然后按照公式(6)与公式(7)进行迭代，直至找到满意的解或者达到最大迭代次数。
本方法将传统四参数非线性问题降解为两参数非线性问题，实现高效率、高精度地求取各种弹性参数，提高了运算效率与精度。
应用示例
为便于理解本发明实施例的方案及其效果，以下给出一个具体应用示例。本领域技术人员应理解，该示例仅为了便于理解本发明，其任何具体细节并非意在以任何方式限制本发明。
根据本发明的叠前混合非线性反演方法包括以下步骤：
对于叠前地震数据进行叠前线性反演，获得横波反射系数Rs的初始值与横纵波速度比Rk的初始值；
将横波反射系数Rs、横纵波速度比Rk、密度梯度RD代入公式(2)，构造组合弹性参量KS；
将公式(1)表示为关于密度梯度RD、纵波反射系数RP和组合弹性参量KS的线性方程即为公式(3)，并进行线性求解，获得密度梯度RD、纵波反射系数RP和组合弹性参量KS；
构造关于横波反射系数Rs和横纵波速度比Rk的目标函数为公式(4)，基于横波反射系数Rs的初始值与横纵波速度比Rk的初始值，通过标准粒子群算法对目标函数进行非线性求解，获得最终的横波反射系数与最终的横纵波速度比。
图2a、图2b分别示出了根据本发明的一个实施例的目标函数的误差值与部分放大的目标函数的误差值的示意图，从图中可以看出：随着迭代次数的增加，目标函数的误差值逐渐减小，基本上从第15次开始趋于稳定，稳定后的误差值很小，几乎可被忽略。
图3a、图3b分别示出了根据本发明的一个实施例的横波反射系数的误差值与横纵波速度比的误差值的示意图。虽然由于非线性量在求解过程中陷入局部寻优的概率不同，横波反射系数Rs和横纵波速度比Rk的误差变化规律有所差别，但整体处于很小的数值范围。相比Rs而言，Rk的非线性程度更高，但其获得了相对较小数值的误差程度，这是因为在本发明的相对误差的计算方式下，横纵波速度比相对反射系数较大的数值会对相对降低误差的统计额度。
图4a、图4b、图4c、图4d分别示出了根据本发明的一个实施例的纵波反射系数、横波反射系数、密度梯度、横纵波速度比与理论值的对比示意图。
进一步对比本发明的实施例所得的四项参数最终结果与理论值的差别，如图4a、图4b、图4c、图4d所示，纵波反射系数Rp和密度梯度RD的误差几乎为零，而横波反射系数Rs和横纵波速度比Rk的误差也十分小，非常准确地描述了理论值的趋势和细节。
图5a、图5b、图5c分别示出了根据叠前线性反演、叠前非线性反演与本发明的一个实施例的纵波反射系数的示意图，图5d示出了三者的误差对比示意图。
图6a、图6b、图6c分别示出了根据叠前线性反演、叠前非线性反演与本发明的一个实施例的横波反射系数的示意图，图6d示出了三者的误差对比示意图。
图7a、图7b、图7c分别示出了根据叠前线性反演、叠前非线性反演与本发明的一个实施例的密度梯度的示意图，图7d示出了三者的误差对比示意图。
图8a、图8b、图8c分别示出了根据叠前线性反演、叠前非线性反演与本发明的一个实施例的横纵波速度比的示意图，图8d示出了三者的误差对比示意图。
叠前线性反演和本发明的叠前混合非线性反演方法均取得了非常准确的Rp和RD结果，因为这两种方法在求解过程中均采用线性办法求解纵波反射系数，其结果为完全独立的线性量。另外，对于非线性量Rs和横纵波速度比Rk的各种方法求解结果均存在不同程度的误差，叠前非线性反演的误差程度最大。一方面这是因为横波参量的参数可解性本身较低，其在非线性方程组式中又是与横纵波速度比密切相关的非线性变量，因而求解难度较大；但更主要的原因还是由于四参量的非线性寻优路径或策略严重降低了单一弹性参量的求解精度。由于对横纵波速度比直接为0.5的初始运算近似，叠前线性反演的Rs和横纵波速度比结果在幅值上与理论值有所偏差；而本发明的叠前混合非线性反演方法只非线性求解两个参量，且其初始值与约束条件以相对准确的迭代反演结果为基础，因而最终结果品质最高。
图9示出了根据叠前线性反演、叠前非线性反演与本发明的一个实施例的运算时间的对比示意图。叠前非线性反演方法的运算时间最多，为233.35秒，叠前线性反演耗时最少，为0.134秒。本发明的叠前混合非线性反演方法耗时4.25秒，相比于其他两类方法的精度和时间，其工作效率最高。
综上所述，本方法将传统四参数非线性问题降解为两参数非线性问题，实现高效率、高精度地求取各种弹性参数，提高了运算效率与精度。
根据本发明提出了一种计算机可读存储介质，其上存储有计算机程序，其中，程序被处理器执行时实现以下步骤：对于叠前地震数据进行叠前线性反演，获得横波反射系数Rs的初始值与横纵波速度比Rk的初始值；基于横波反射系数Rs、横纵波速度比Rk、密度梯度RD，构造组合弹性参量KS；将叠前纵波反射系数公式表示为关于密度梯度RD、纵波反射系数RP和组合弹性参量KS的线性方程，对叠前纵波反射系数公式进行线性求解，获得密度梯度RD、纵波反射系数RP和组合弹性参量KS；构造关于横波反射系数Rs和横纵波速度比Rk的目标函数，基于横波反射系数Rs的初始值与横纵波速度比Rk的初始值对目标函数进行非线性求解，进而获得最终的横波反射系数与最终的横纵波速度比。
在一个示例中，叠前纵波反射系数公式为：
其中，Rpp(θ)表示叠前纵波反射系数，θ表示反射角与入射角的平均值。
在一个示例中，组合弹性参量KS为：
在一个示例中，目标函数为：
其中，F(Rs,Rk)表示目标函数。
在一个示例中，目标函数为：
其中，F(Rs,Rk)表示目标函数。
在一个示例中，通过标准粒子群算法对目标函数进行非线性求解。
以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择，旨在最好地解释各实施例的原理、实际应用或对市场中的技术的改进，或者使本技术领域的其它普通技术人员能理解本文披露的各实施例。

Claims (10)

1.一种叠前混合非线性反演方法，包括：
对于叠前地震数据进行叠前线性反演，获得横波反射系数Rs的初始值与横纵波速度比Rk的初始值；
基于所述横波反射系数Rs、所述横纵波速度比Rk、密度梯度RD，构造组合弹性参量KS；
将叠前纵波反射系数公式表示为关于所述密度梯度RD、纵波反射系数RP和组合弹性参量KS的线性方程，对所述叠前纵波反射系数公式进行线性求解，获得所述密度梯度RD、纵波反射系数RP和所述组合弹性参量KS；
构造关于所述横波反射系数Rs和横纵波速度比Rk的目标函数，基于所述横波反射系数Rs的初始值与横纵波速度比Rk的初始值对所述目标函数进行非线性求解，进而获得最终的横波反射系数与最终的横纵波速度比。
2.根据权利要求1所述的叠前混合非线性反演方法，其中，所述叠前纵波反射系数公式为：
其中，Rpp(θ)表示叠前纵波反射系数，θ表示反射角与入射角的平均值。
3.根据权利要求1所述的叠前混合非线性反演方法，其中，所述组合弹性参量KS为：
4.根据权利要求1所述的叠前混合非线性反演方法，其中，所述目标函数为：
其中，F(Rs,Rk)表示目标函数。
5.根据权利要求1所述的叠前混合非线性反演方法，其中，所述目标函数为：
其中，F(Rs,Rk)表示目标函数。
6.根据权利要求1所述的叠前混合非线性反演方法，其中，通过标准粒子群算法对所述目标函数进行非线性求解。
7.一种计算机可读存储介质，其上存储有计算机程序，其中，所述程序被处理器执行时实现以下步骤：
对于叠前地震数据进行叠前线性反演，获得横波反射系数Rs的初始值与横纵波速度比Rk的初始值；
基于所述横波反射系数Rs、所述横纵波速度比Rk、密度梯度RD，构造组合弹性参量KS；
将叠前纵波反射系数公式表示为关于所述密度梯度RD、纵波反射系数RP和组合弹性参量KS的线性方程，对所述叠前纵波反射系数公式进行线性求解，获得所述密度梯度RD、纵波反射系数RP和所述组合弹性参量KS；
构造关于所述横波反射系数Rs和横纵波速度比Rk的目标函数，基于所述横波反射系数Rs的初始值与横纵波速度比Rk的初始值对所述目标函数进行非线性求解，进而获得最终的横波反射系数与最终的横纵波速度比。
8.根据权利要求7所述的计算机可读存储介质，其中，所述叠前纵波反射系数公式为：
其中，Rpp(θ)表示叠前纵波反射系数，θ表示反射角与入射角的平均值。
9.根据权利要求7所述的计算机可读存储介质，其中，所述组合弹性参量KS为：
10.根据权利要求7所述的计算机可读存储介质，其中，所述目标函数为：
其中，F(Rs,Rk)表示目标函数。
CN201710464861.9A 2017-06-19 2017-06-19 叠前混合非线性反演方法及计算机可读存储介质 Active CN109143346B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710464861.9A CN109143346B (zh) 2017-06-19 2017-06-19 叠前混合非线性反演方法及计算机可读存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710464861.9A CN109143346B (zh) 2017-06-19 2017-06-19 叠前混合非线性反演方法及计算机可读存储介质

Publications (2)

Publication Number Publication Date
CN109143346A true CN109143346A (zh) 2019-01-04
CN109143346B CN109143346B (zh) 2020-10-20

Family

ID=64804592

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710464861.9A Active CN109143346B (zh) 2017-06-19 2017-06-19 叠前混合非线性反演方法及计算机可读存储介质

Country Status (1)

Country Link
CN (1) CN109143346B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111948712A (zh) * 2020-08-10 2020-11-17 中海石油(中国)有限公司 一种基于深度域地震记录的叠前线性反演方法
CN113031058A (zh) * 2021-02-26 2021-06-25 河海大学 一种基于反射系数精确式的页岩vti储层叠前混合反演方法
CN113391349A (zh) * 2020-03-11 2021-09-14 中国石油天然气股份有限公司 随机多维反射系数反演方法及系统

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012050880A3 (en) * 2010-09-29 2012-12-13 Chevron U.S.A. Inc. Device and method for imaging of non-linear and linear properties of formations surrounding a borehole
CN102854528A (zh) * 2012-07-13 2013-01-02 孙赞东 粒子群优化算法叠前非线性反演方法
CN103675911A (zh) * 2014-01-06 2014-03-26 中国石油大学(华东) 一种基于纵波和转换横波联合反演截距和梯度的方法
WO2014085614A2 (en) * 2012-11-27 2014-06-05 Chevron U.S.A. Inc. System and method for generating 3d images of non-linear properties of rock formation using surface seismic or surface to borehole seismic or both
EP2755055A2 (en) * 2013-01-11 2014-07-16 CGG Services SA Dip tomography for estimating depth velocity models by inverting pre-stack dip information present in migrated/un-migrated pre-/post-stack seismic data
CN106597544A (zh) * 2016-11-25 2017-04-26 中国石油天然气股份有限公司 致密油气藏脆性的预测方法及装置

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012050880A3 (en) * 2010-09-29 2012-12-13 Chevron U.S.A. Inc. Device and method for imaging of non-linear and linear properties of formations surrounding a borehole
CN102854528A (zh) * 2012-07-13 2013-01-02 孙赞东 粒子群优化算法叠前非线性反演方法
WO2014085614A2 (en) * 2012-11-27 2014-06-05 Chevron U.S.A. Inc. System and method for generating 3d images of non-linear properties of rock formation using surface seismic or surface to borehole seismic or both
EP2755055A2 (en) * 2013-01-11 2014-07-16 CGG Services SA Dip tomography for estimating depth velocity models by inverting pre-stack dip information present in migrated/un-migrated pre-/post-stack seismic data
CN103675911A (zh) * 2014-01-06 2014-03-26 中国石油大学(华东) 一种基于纵波和转换横波联合反演截距和梯度的方法
CN106597544A (zh) * 2016-11-25 2017-04-26 中国石油天然气股份有限公司 致密油气藏脆性的预测方法及装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
于红玉 等: "A Modified Pre-stack Inversion Method Using Chaotic Quantum Particle Swarm Optimization", 《MADRID 2015》 *
张远银 等: "P-P与P-SV波联合反演方法分类与对比", 《石油物探》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113391349A (zh) * 2020-03-11 2021-09-14 中国石油天然气股份有限公司 随机多维反射系数反演方法及系统
CN111948712A (zh) * 2020-08-10 2020-11-17 中海石油(中国)有限公司 一种基于深度域地震记录的叠前线性反演方法
CN111948712B (zh) * 2020-08-10 2021-08-10 中海石油(中国)有限公司 一种基于深度域地震记录的叠前线性反演方法
CN113031058A (zh) * 2021-02-26 2021-06-25 河海大学 一种基于反射系数精确式的页岩vti储层叠前混合反演方法
CN113031058B (zh) * 2021-02-26 2022-06-24 河海大学 一种基于反射系数精确式的页岩vti储层叠前混合反演方法

Also Published As

Publication number Publication date
CN109143346B (zh) 2020-10-20

Similar Documents

Publication Publication Date Title
CN107290782B (zh) 储层孔隙度、含水饱和度和泥质含量参数同时反演新方法
CN102854528B (zh) 粒子群优化算法叠前非线性反演方法
CN108037528B (zh) 少井区基于统计岩石物理建模的孔隙度预测方法及系统
CN109143346A (zh) 叠前混合非线性反演方法及计算机可读存储介质
CN107301282B (zh) 基于多源监测时序数据的混凝土坝力学参数反演方法
CN104570101A (zh) 一种基于粒子群算法的avo三参数反演方法
CN104237937B (zh) 叠前地震反演方法及其系统
Li et al. Long-term predictions of turbulence by implicit U-Net enhanced Fourier neural operator
CN106324674A (zh) 一种页岩气toc叠前地震反演预测方法
CN107808032A (zh) 页岩气储层含水饱和度计算方法及系统
CN103454677A (zh) 基于粒子群与线性加法器结合的地震数据反演方法
Gingery et al. Shear stress-strain curves based on the G/Gmax logic: A procedure for strength compatibility
CN104316961A (zh) 获取风化层的地质参数的方法
Zhang et al. Efficient surrogate modeling based on improved vision transformer neural network for history matching
CN114036831A (zh) 一种待检测工程场地边坡的岩土力学参数的实时检测方法
CN108763164A (zh) 煤与瓦斯突出反演相似度的评价方法
Liu et al. Prediction of ultimate axial load-carrying capacity for driven piles using machine learning methods
Stringer Implementing approximate Bayesian inference using adaptive quadrature: the aghq package
CN107643541B (zh) 基于速度模型的速度谱解释方法
CN116628886A (zh) 基于有限工程数据的盾构机穿越施工参数实时优化方法
CN110231653A (zh) 构建双向约束初始模型的方法及装置
Plante et al. A generalized stress correction scheme for the Maxwell elasto-brittle rheology: impact on the fracture angles and deformations
Scarinci et al. Robust Bayesian moment tensor inversion with optimal transport misfits: layered medium approximations to the 3-D SEG-EAGE overthrust velocity model
Wang et al. Bayesian perspective on ground property variability for geotechnical practice
Fredlund et al. Keynote Lecture: Developments in landslide analysis

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant