CN114594515A - 基于缓变各向异性的井控速度反演方法 - Google Patents

基于缓变各向异性的井控速度反演方法 Download PDF

Info

Publication number
CN114594515A
CN114594515A CN202011420614.7A CN202011420614A CN114594515A CN 114594515 A CN114594515 A CN 114594515A CN 202011420614 A CN202011420614 A CN 202011420614A CN 114594515 A CN114594515 A CN 114594515A
Authority
CN
China
Prior art keywords
well
depth
offset
ray
slowness
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
CN202011420614.7A
Other languages
English (en)
Other versions
CN114594515B (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
Geophysical Research Institute of Sinopec Shengli Oilfield Co
Original Assignee
China Petroleum and Chemical Corp
Geophysical Research Institute of Sinopec Shengli Oilfield Co
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, Geophysical Research Institute of Sinopec Shengli Oilfield Co filed Critical China Petroleum and Chemical Corp
Priority to CN202011420614.7A priority Critical patent/CN114594515B/zh
Publication of CN114594515A publication Critical patent/CN114594515A/zh
Application granted granted Critical
Publication of CN114594515B publication Critical patent/CN114594515B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/282—Application of seismic models, synthetic seismograms
    • 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/301—Analysis for determining seismic cross-sections or geostructures
    • G01V1/302—Analysis for determining seismic cross-sections or geostructures in 3D data cubes
    • 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/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/50—Corrections or adjustments related to wave propagation
    • G01V2210/51—Migration
    • G01V2210/512—Pre-stack
    • 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/66—Subsurface modeling
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30—Assessment of water resources

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

基于缓变各向异性的井控速度反演方法
技术领域
本发明涉及油气勘探地震资料处理技术领域，特别涉及一种基于缓变各向异性的井控速度反演方法。
背景技术
目前，生产中常用的速度建模方法是基于射线理论的剩余曲率速度分析，主要应用反射波的走时信息，反演的结果往往是以低波数段为主的背景速度场，难以满足高精度地震勘探的需求。如何提高速度模型的精度，实现复杂地区的精确成像是一个急需解决的关键性问题，也是地震成像领域乃至地震勘探领域的一大难点。基于波动方程理论的速度建模方法在理论上可以得到包含高频成分的速度场，分辨率高，能更好地适应速度变化剧烈的地区，但存在许多尚未解决的问题。波动方程偏移速度分析理论和实际应用尚未完善，初始模型问题和对速度模型的敏感度是该方法面临的一大问题，并且计算量庞大，处理分析不灵活。全波形反演虽然在理论上已基本完善，但其计算量太大，严重依赖低频信息和大偏移距地震数据。而对于陆上地震勘探，如何提取地震子波以及模拟复杂波场也是摆在全波形反演面前的难题。复杂地区的地震资料还存在信噪比低的问题，因此，应用全波形反演解决复杂地区的速度建模问题是短期内无法实现的。
除上述问题外，许多工区还存在各向异性现象。各向异性作用和速度相互耦合，显著影响深度成像的准确性。
发明内容
本发明的目的是针对常规偏移速度分析结果分辨率低、各向异性和速度相互耦合的问题，提供一种基于缓变各向异性的井控速度反演方法，以反演精细的地下速度场，为深度域地震成像提供技术支持。
本发明提供一种基于缓变各向异性的井控速度反演方法，该方法包括：
步骤一，整理工区内的测井速度数据，去除其中的异常值，并做平滑处理；
步骤二，基于时深转换得到的先验速度模型，对预处理的叠前地震数据进行射线类三维叠前深度偏移，获取深度偏移剖面MIG(x,y,z)和偏移距共成像点道集CIP(x,y,z,h)，其中，x、y、z、h分别为成像点的主测线方向坐标、联络测线方向坐标、深度和偏移距；
步骤三，在所述偏移剖面MIG(x,y,z)中对所述成像点进行x/y方向局部倾斜叠加扫描，以最大相似谱对应的倾角作为该点倾角，逐点扫描得x方向倾角场
Figure BDA0002822059780000022
和y方向倾角场
Figure BDA0002822059780000023
步骤四，对各所述成像点道集CIP(x,y,z,h)进行γ谱扫描，获得最大相似谱对应的γ场γ(x,y,z)；
步骤五，基于所述倾角场和γ场进行综合分析和质量监控，筛选出所述深度偏移剖面MIG(x,y,z)中的反射点；
步骤六，以各所述反射点作为出射点，进行射线追踪到地表，获得射线路径并基于射线终点之间的偏移距计算剩余深度差；
步骤七，将预处理的测井速度沿所述倾角场
Figure BDA0002822059780000024
和所述
Figure BDA0002822059780000025
外推，加权求和得到井速度插值模型；
步骤八，基于射线路径及对应的深度差、平滑约束以及从所述井速度插值模型中提取的信息，构建反演方程组；
步骤九，求解所述反演方程组，获取所述速度模型更新量，更新所述速度模型。
进一步的，在步骤三中，x/y方向的相似谱的计算公式如下：
Figure BDA0002822059780000021
其中，SEMdipx、SEMdipy分别为x、y方向的相似谱，z为深度，k为倾角的正切值，trace为地震道，x0、y0为中心道位置，xj、yj为倾斜叠加的地震道位置，L为扫描窗口道数。
进一步的，在步骤四中，γ谱扫描所采用的成像深度zmig与γ的关系
式为：
Figure BDA0002822059780000031
γ相似谱的计算公式如下：
Figure BDA0002822059780000032
其中，SEMrmo为γ相似谱，Δz为加权平均的窗口宽度。
进一步的，在步骤六中，从拾取的倾角场中获得反射点处的x/y方向倾角，计算反射点法线方向的射线以及在法线两侧以相同角度出射的射线，每次计算均以不同出射角和方位角进行射线追踪，并以到达地表边界作为该条射线的终止条件。
进一步的，在步骤六中，具体采用的三维射线追踪公式如下：
Figure BDA0002822059780000033
其中，px、py、pz为慢度分量，v是模型速度，上述方程组进一步通过Runge-Kutta法进行数值求解。
进一步的，在步骤六中，剩余深度差Δz的计算公式如下：
Figure BDA0002822059780000034
其中，z0为反射点在偏移距成像点道集中的零偏移距深度。
进一步的，在步骤八中，基于射线路径及对应的深度差构建的方程如下：
L0Δs0 cosβh-LhΔsh＝2sc·cosβh·cosα·Δz，
其中，L0为零偏移距射线路径在各网格内的长度，Δs0为零偏移距射线路径上的慢度更新量，βh为偏移距h对应的出射角，Lh为偏移距h的射线路径在各网格内的长度，Δsh为偏移距h的射线路径上的慢度更新量，sc为反射点慢度，α为地层倾角，Δz为剩余深度差。
进一步的，在步骤八中，从所述井速度插值模型中提取相邻速度的相对变化约束反演过程，具体方程如下：
Figure BDA0002822059780000041
其中，si、Δsi、
Figure BDA0002822059780000042
分别为第i点慢度、慢度更新量和井慢度，si+1、Δsi+1、
Figure BDA0002822059780000043
分别为第i+1点慢度、慢度更新量和井慢度。
进一步的，在所述步骤八中，构建的反演方程组形式如下：
Figure BDA0002822059780000044
其中，L、Δt为由权利要求7中的方程构成的系数矩阵和右端项，R为由三维拉普拉斯算子构成的平滑矩阵，W、e为由权利要求8中的方程构成的系数矩阵和右端项，Δs表示模型更新量，ε1、ε2为权重系数。
进一步的，在所述步骤九中，采用并行LSQR算法求解大型稀疏的反演方程组。
本发明有益技术效果如下：
本发明实施例提供一种基于缓变各向异性的井控速度反演方法，可以构建偏移速度场。与常规偏移速度分析相比，本发明实施例方法应用测井信息、构造信息，并考虑各向异性影响，所得的各向同性速度模型更符合实际情况，分辨率更高，为深度域的正确成像奠定基础，具有广阔的应用前景。
附图说明
图1为本发明实施例中基于缓变各向异性的井控速度反演方法的流程图；
图2为本发明实施例中工区内某口井预处理后的测井速度曲线显示图；
图3为图1中所述方法的先验速度模型及其偏移结果显示图，其中，图3(a)显示为先验速度模型，图3(b)显示为基于先验模型得到的深度偏移剖面，图3(c)显示为基于先验模型得到的某成像点道集；
图4为本发明实施例中基于深度偏移剖面倾斜叠加得到的倾角场显示图，图中黑色背景条纹为偏移剖面中的地震波同相轴；
图5为本发明实施例中基于深度偏移剖面γ扫描得到的γ场显示图，图中黑色背景条纹为偏移剖面中的地震波同相轴；
图6为本发明实施例中三维射线覆盖显示图；
图7为本发明实施例中工区井速度模型显示图，图中黑色背景条纹为偏移剖面中的地震波同相轴；
图8为本发明实施例中常规偏移速度反演的最终结果显示图，其中图8(a)为反演速度模型，其中黑圈内为对比区域；图8(b)为最终的深度偏移剖面；图8(c)为最终的成像点道集；
图9为图1所示方法的最终结果显示图；
图9(a)为反演速度模型，其中黑圈内为对比区域；
图9(b)为最终的深度偏移剖面；
图9(c)为最终的成像点道集。
具体实施方式
本部分将详细描述本发明的具体实施例，本发明之较佳实施例在附图中示出，附图的作用在于用图形补充说明书文字部分的描述，使人能够直观地、形象地理解本发明的每个技术特征和整体技术方案，但其不能理解为对本发明保护范围的限制。
如图所示，图1为本发明基于缓变各向异性的井控速度反演方法的流程图。该方法包括以下步骤：
步骤一，整理工区内的测井速度数据，去除其中的异常值，并做平滑处理，图2显示为其中一口井处理后的井速度曲线；
步骤二，基于时深转换得到先验速度模型，如图3(a)所示。对预处理的叠前地震数据进行射线类三维叠前深度偏移，获取深度偏移剖面MIG(x,y,z)如图3(b)所示，和偏移距共成像点道集CIP(x,y,z,h)，如图3(c)所示，其中，x、y、z、h分别为成像点的主测线方向坐标、联络测线方向坐标、深度和偏移距。
步骤三，在所述偏移剖面MIG(x,y,z)中对所述成像点进行x/y方向局部倾斜叠加扫描，以最大相似谱对应的倾角作为该点倾角，逐点扫描得x方向倾角场
Figure BDA0002822059780000061
和y方向倾角场
Figure BDA0002822059780000062
如图4所示。
优选的，x/y方向的相似谱的计算公式如下：
Figure BDA0002822059780000063
其中，SEMdipx、SEMdipy分别为x、y方向的相似谱，z为深度，k为倾角的正切值，trace为地震道，x0、y0为中心道位置，xj、yj为倾斜叠加的地震道位置，L为扫描窗口道数。
步骤四，对各所述成像点道集CIP(x,y,z,h)进行γ谱扫描，获得最大相似谱对应的γ场γ(x,y,z)，如图5所示。
优选的，γ谱扫描所采用的成像深度zmig与γ的关系式为：
Figure BDA0002822059780000064
γ相似谱的计算公式如下：
Figure BDA0002822059780000065
其中，SEMrmo为γ相似谱，Δz为加权平均的窗口宽度。
步骤五，基于所述倾角场和γ场进行综合分析和质量监控，筛选出所述深度偏移剖面MIG(x,y,z)中的反射点。
步骤六，以各所述反射点作为出射点，进行射线追踪到地表，获得射线路径并基于射线终点之间的偏移距计算剩余深度差，如图6为三维射线覆盖显示图。
具体地，从拾取的倾角场中获得反射点处的x/y方向倾角，计算反射点法线方向的射线以及在法线两侧以相同角度出射的射线，每次计算均以不同出射角和方位角进行射线追踪，并以到达地表边界作为该条射线的终止条件。
作为本发明的一个优选实施例，在步骤六中，具体采用的三维射线追踪公式如下：
Figure BDA0002822059780000071
其中，px、py、pz为慢度分量，v是模型速度，上述方程组进一步通过Runge-Kutta法进行数值求解。
进一步的，在步骤六中，剩余深度差Δz的计算公式如下：
Figure BDA0002822059780000072
其中，z0为反射点在偏移距成像点道集中的零偏移距深度。
步骤七，将预处理的测井速度沿所述倾角场
Figure BDA0002822059780000073
和所述
Figure BDA0002822059780000074
外推，加权求和得到井速度插值模型。
具体地，对工区内的测井速度进行处理，如去除异常值、平滑，并沿倾角场外推，加权求和得到井速度模型，如图7所示。
步骤八，基于射线路径及对应的深度差、平滑约束以及从所述井速度插值模型中提取的信息，构建反演方程组。
进一步的，在步骤八中，基于射线路径及对应的深度差构建的方程如下：
L0Δs0 cosβh-LhΔsh＝2sc·cosβh·cosα·Δz，
其中，L0为零偏移距射线路径在各网格内的长度，Δs0为零偏移距射线路径上的慢度更新量，βh为偏移距h对应的出射角，Lh为偏移距h的射线路径在各网格内的长度，Δsh为偏移距h的射线路径上的慢度更新量，sc为反射点慢度，α为地层倾角，Δz为剩余深度差。
作为本发明的一个具体实施例，在步骤八中，从井速度插值模型中提取相邻速度的相对变化约束反演过程，具体方程如下：
Figure BDA0002822059780000075
其中，si、Δsi、
Figure BDA0002822059780000076
分别为第i点慢度、慢度更新量和井慢度，si+1、Δsi+1、
Figure BDA0002822059780000077
分别为第i+1点慢度、慢度更新量和井慢度。
进一步的，在所述步骤八中，构建的反演方程组形式如下：
Figure BDA0002822059780000081
其中，L、Δt为由权利要求7中的方程构成的系数矩阵和右端项，R为由三维拉普拉斯算子构成的平滑矩阵，W、e为由权利要求8中的方程构成的系数矩阵和右端项，Δs表示模型更新量，ε1、ε2为权重系数。
步骤九，求解所述反演方程组，获取速度模型更新量，更新速度模型。
优选的，在所述步骤九中，采用并行LSQR算法求解大型稀疏的反演方程组。
经两次迭代后，可得到最终的结果。常规偏移速度反演的最终结果如图8所示，其中，图8(a)为反演速度模型，图8(b)为最终的深度偏移剖面，图8(c)为最终的成像点道集。本发明实施例方法的最终结果如图9所示，其中，图9(a)为反演速度模型，图9(b)为最终的深度偏移剖面，图9(c)为最终的成像点道集。对比图8和图9可知，两种方法的结果均可使同相轴拉平，使深度偏移进行有效成像。而对比两个速度模型图8(a)和图9(a)可知，本发明实施例方法的结果具有更高的分辨率，尤其是黑圈内的区域。
对比深度偏移剖面可知，本发明实施例方法的速度模型更符合成像剖面中的构造趋势。该具体示例证明了本发明实施例方法是一种有效的成像域层析反演方法。
以上所述仅为本发明的较佳实施例而己，并不以本发明为限制，凡在本发明的精神和原则之内所作的均等修改、等同替换和改进等，均应包含在本发明的专利涵盖范围内。

Claims (10)

1.一种基于缓变各向异性的井控速度反演方法，其特征在于，该方法包括：
步骤一，整理工区内的测井速度数据，去除其中的异常值，并做平滑处理；
步骤二，基于时深转换得到的先验速度模型，对预处理的叠前地震数据进行射线类三维叠前深度偏移，获取深度偏移剖面MIG(x,y,z)和偏移距共成像点道集CIP(x,y,z,h)，其中，x、y、z、h分别为成像点的主测线方向坐标、联络测线方向坐标、深度和偏移距；
步骤三，在所述偏移剖面MIG(x,y,z)中对所述成像点进行x/y方向局部倾斜叠加扫描，以最大相似谱对应的倾角作为该点倾角，逐点扫描得x方向倾角场
Figure FDA0002822059770000011
和y方向倾角场
Figure FDA0002822059770000012
步骤四，对各所述成像点道集CIP(x,y,z,h)进行γ谱扫描，获得最大相似谱对应的γ场γ(x,y,z)；
步骤五，基于所述倾角场和γ场进行综合分析和质量监控，筛选出所述深度偏移剖面MIG(x,y,z)中的反射点；
步骤六，以各所述反射点作为出射点，进行射线追踪到地表，获得射线路径并基于射线终点之间的偏移距计算剩余深度差；
步骤七，将预处理的测井速度沿所述倾角场
Figure FDA0002822059770000013
和所述
Figure FDA0002822059770000014
外推，加权求和得到井速度插值模型；
步骤八，基于射线路径及对应的深度差、平滑约束以及从所述井速度插值模型中提取的信息，构建反演方程组；
步骤九，求解所述反演方程组，获取所述速度模型更新量，更新所述速度模型。
2.根据权利要求1所述的基于缓变各向异性的井控速度反演方法，其特征在于，在步骤三中，x/y方向的相似谱的计算公式如下：
Figure FDA0002822059770000015
其中，SEMdipx、SEMdipy分别为x、y方向的相似谱，z为深度，k为倾角的正切值，trace为地震道，x0、y0为中心道位置，xj、yj为倾斜叠加的地震道位置，L为扫描窗口道数。
3.根据权利要求1所述的基于缓变各向异性的井控速度反演方法，其特征在于，在步骤四中，γ谱扫描所采用的成像深度zmig与γ的关系式为：
Figure FDA0002822059770000021
γ相似谱的计算公式如下：
Figure FDA0002822059770000022
其中。SEMrmo为γ相似谱，Δz为加权平均的窗口宽度。
4.根据权利要求1所述的基于缓变各向异性的井控速度反演方法，其特征在于，在步骤六中，从拾取的倾角场中获得反射点处的x/y方向倾角，计算反射点法线方向的射线以及在法线两侧以相同角度出射的射线，每次计算均以不同出射角和方位角进行射线追踪，并以到达地表边界作为该条射线的终止条件。
5.根据权利要求4所述的基于缓变各向异性的井控速度反演方法，其特征在于，在步骤六中，具体采用的三维射线追踪公式如下：
Figure FDA0002822059770000023
其中，px、py、pz为慢度分量，v是模型速度，上述方程组进一步通过Runge-Kutta法进行数值求解。
6.根据权利要求5所述的基于缓变各向异性的井控速度反演方法，其特征在于，在步骤六中，剩余深度差Δz的计算公式如下：
Figure FDA0002822059770000024
其中，z0为反射点在偏移距成像点道集中的零偏移距深度。
7.根据权利要求1所述的基于缓变各向异性的井控速度反演方法，其特征在于，在步骤八中，基于射线路径及对应的深度差构建的方程如下：L0Δs0cosβh-LhΔsh＝2sc·cosβh·cosα·Δz，其中，L0为零偏移距射线路径在各网格内的长度，Δs0为零偏移距射线路径上的慢度更新量，βh为偏移距h对应的出射角，Lh为偏移距h的射线路径在各网格内的长度，Δsh为偏移距h的射线路径上的慢度更新量，sc为反射点慢度，α为地层倾角，Δz为剩余深度差。
8.根据权利要求7所述的基于缓变各向异性的井控速度反演方法，其特征在于，在步骤八中，从所述井速度插值模型中提取相邻速度的相对变化约束反演过程，具体方程如下：
Figure FDA0002822059770000031
其中，si、Δsi、
Figure FDA0002822059770000032
分别为第i点慢度、慢度更新量和井慢度，si+1、Δsi+1、
Figure FDA0002822059770000033
分别为第i+1点慢度、慢度更新量和井慢度。
9.根据权利要求8所述的基于缓变各向异性的井控速度反演方法，其特征在于，在所述步骤八中，构建的反演方程组形式如下：
Figure FDA0002822059770000034
其中，L、Δt为由权利要求7中的方程构成的系数矩阵和右端项，R为由三维拉普拉斯算子构成的平滑矩阵，W、e为由权利要求8中的方程构成的系数矩阵和右端项，Δs表示模型更新量，ε1、ε2为权重系数。
10.根据权利要求1所述的基于缓变各向异性的井控速度反演方法，其特征在于，在所述步骤九中，采用并行LSQR算法求解大型稀疏的反演方程组。
CN202011420614.7A 2020-12-07 2020-12-07 基于缓变各向异性的井控速度反演方法 Active CN114594515B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011420614.7A CN114594515B (zh) 2020-12-07 2020-12-07 基于缓变各向异性的井控速度反演方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011420614.7A CN114594515B (zh) 2020-12-07 2020-12-07 基于缓变各向异性的井控速度反演方法

Publications (2)

Publication Number Publication Date
CN114594515A true CN114594515A (zh) 2022-06-07
CN114594515B CN114594515B (zh) 2024-03-29

Family

ID=81802848

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011420614.7A Active CN114594515B (zh) 2020-12-07 2020-12-07 基于缓变各向异性的井控速度反演方法

Country Status (1)

Country Link
CN (1) CN114594515B (zh)

Citations (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6813374B1 (en) * 2001-04-25 2004-11-02 Analogic Corporation Method and apparatus for automatic image quality assessment
US20090213692A1 (en) * 2008-02-22 2009-08-27 Martinez Ruben D Method for three dimensional seismic travel time tomography in transversely isotropic media
US20090257308A1 (en) * 2008-04-11 2009-10-15 Dimitri Bevc Migration velocity analysis methods
US20100135115A1 (en) * 2008-12-03 2010-06-03 Chevron U.S.A. Inc. Multiple anisotropic parameter inversion for a tti earth model
CN103703391A (zh) * 2011-03-30 2014-04-02 埃克森美孚上游研究公司 使用频谱整形的全波场反演的收敛速度
CN103984012A (zh) * 2014-04-16 2014-08-13 孙赞东 基于叠前高斯束深度偏移的绕射波场分离方法
US20150117144A1 (en) * 2013-10-29 2015-04-30 Hongchuan Sun Inversion-Based Reflector Dip Estimation
US20160187512A1 (en) * 2014-01-14 2016-06-30 Cgg Services Sa Detecting and estimating anisotropy errors using full waveform inversion and ray based tomography
WO2016153567A1 (en) * 2015-03-26 2016-09-29 Halliburton Energy Services, Inc. Drilling fluid property determination
CN106443777A (zh) * 2016-08-15 2017-02-22 中国石油集团川庆钻探工程有限公司地球物理勘探公司 转换波各向异性参数谱的获取方法及转换波速度分析方法
CN106569259A (zh) * 2015-10-09 2017-04-19 中国石油化工股份有限公司 基于构造倾角的正则化层析速度反演方法和装置
CN106646613A (zh) * 2016-12-28 2017-05-10 中国石油化工股份有限公司 深度域多尺度井控建模与成像联合处理方法
CN107505651A (zh) * 2017-06-26 2017-12-22 中国海洋大学 地震初至波和反射波联合斜率层析成像方法
CN109100786A (zh) * 2018-06-21 2018-12-28 中国石油天然气集团有限公司 深度域品质因子的确定方法和装置
CN109581496A (zh) * 2018-10-29 2019-04-05 中国石油天然气集团有限公司 层析反演中的双向构造倾角约束方法及装置
US20200174147A1 (en) * 2017-08-10 2020-06-04 Seismic Apparition Gmbh Method for seismic data acquisition and processing

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108072892B (zh) * 2016-11-09 2020-01-10 中国石油化工股份有限公司 一种自动化的地质构造约束层析反演方法

Patent Citations (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6813374B1 (en) * 2001-04-25 2004-11-02 Analogic Corporation Method and apparatus for automatic image quality assessment
US20090213692A1 (en) * 2008-02-22 2009-08-27 Martinez Ruben D Method for three dimensional seismic travel time tomography in transversely isotropic media
US20090257308A1 (en) * 2008-04-11 2009-10-15 Dimitri Bevc Migration velocity analysis methods
US20100135115A1 (en) * 2008-12-03 2010-06-03 Chevron U.S.A. Inc. Multiple anisotropic parameter inversion for a tti earth model
CN103703391A (zh) * 2011-03-30 2014-04-02 埃克森美孚上游研究公司 使用频谱整形的全波场反演的收敛速度
US20150117144A1 (en) * 2013-10-29 2015-04-30 Hongchuan Sun Inversion-Based Reflector Dip Estimation
US20160187512A1 (en) * 2014-01-14 2016-06-30 Cgg Services Sa Detecting and estimating anisotropy errors using full waveform inversion and ray based tomography
CN103984012A (zh) * 2014-04-16 2014-08-13 孙赞东 基于叠前高斯束深度偏移的绕射波场分离方法
WO2016153567A1 (en) * 2015-03-26 2016-09-29 Halliburton Energy Services, Inc. Drilling fluid property determination
CN106569259A (zh) * 2015-10-09 2017-04-19 中国石油化工股份有限公司 基于构造倾角的正则化层析速度反演方法和装置
CN106443777A (zh) * 2016-08-15 2017-02-22 中国石油集团川庆钻探工程有限公司地球物理勘探公司 转换波各向异性参数谱的获取方法及转换波速度分析方法
CN106646613A (zh) * 2016-12-28 2017-05-10 中国石油化工股份有限公司 深度域多尺度井控建模与成像联合处理方法
CN107505651A (zh) * 2017-06-26 2017-12-22 中国海洋大学 地震初至波和反射波联合斜率层析成像方法
US20200174147A1 (en) * 2017-08-10 2020-06-04 Seismic Apparition Gmbh Method for seismic data acquisition and processing
CN109100786A (zh) * 2018-06-21 2018-12-28 中国石油天然气集团有限公司 深度域品质因子的确定方法和装置
CN109581496A (zh) * 2018-10-29 2019-04-05 中国石油天然气集团有限公司 层析反演中的双向构造倾角约束方法及装置

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
CHANGKUN JIN ET AL.: "STEREOTOMOGRAPHY OF SEISMIC DATA ACQUIRED ON UNDULANT TOPOGRAPHY", vol. 83 *
秦宁等: "基于角道集的井约束层析速度反演", vol. 46, no. 5, pages 725 - 731 *
金昌昆等: "微测井与方位加权插值精细近地表速度建模技术", vol. 55, no. 2, pages 257 - 265 *

Also Published As

Publication number Publication date
CN114594515B (zh) 2024-03-29

Similar Documents

Publication Publication Date Title
CN108196305B (zh) 一种山地静校正方法
WO2015196779A1 (zh) 利用叠加速度求取高精度地震波速度的方法
US12007515B2 (en) Optimal survey design
CN102841376A (zh) 一种基于起伏地表的层析速度反演方法
CN109444956B (zh) 三维起伏观测面地震斜率层析成像方法
CN103713323A (zh) 一种全方位各向异性保幅成像与抽道集方法
CN103576200A (zh) 低信噪比区浅层波阻抗界面静校正方法
CN105629300B (zh) 提高复杂构造偏移数据信噪比的方法
CN106842297A (zh) 井约束非稳态相位校正方法
CN111175825B (zh) 深度域速度建模方法
CN109188522B (zh) 速度场构建方法及装置
CN109031411B (zh) 基于垂直地震剖面数据的高斯束偏移成像方法及系统
CN109521470B (zh) 分析地质构造对地震反演裂缝密度影响的方法
CN114594515A (zh) 基于缓变各向异性的井控速度反演方法
CN106526668B (zh) 原始波形提取与成像方法
CN111650645B (zh) 一种变偏移距vsp弯线校正处理方法和装置
CN109613614B (zh) 一种vsp倾角滤波器顶点的选取方法
CN114594516B (zh) 一种成像域井震联合多尺度层析反演方法
CN110208851A (zh) 一种基于加权匹配追踪的三维vsp地震资料内插方法
CN112346123B (zh) 一种地震数据处理via双参数分析方法
CN112346125B (zh) 一种地震数据处理vda双参数分析方法
CN112346124B (zh) 一种地震数据处理via双参数成像方法
CN118778114A (zh) 地震资料前处理方法、静校正处理方法及构造成像方法
CN116430450A (zh) 双复杂构造区叠前深度偏移速度建模方法及装置
CN115343781A (zh) 一种基于构造约束的井震联合复杂构造速度建模方法

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