CN116413798A - 各向异性参数的求取方法及装置 - Google Patents
各向异性参数的求取方法及装置 Download PDFInfo
- Publication number
- CN116413798A CN116413798A CN202111675949.8A CN202111675949A CN116413798A CN 116413798 A CN116413798 A CN 116413798A CN 202111675949 A CN202111675949 A CN 202111675949A CN 116413798 A CN116413798 A CN 116413798A
- Authority
- CN
- China
- Prior art keywords
- anisotropic
- velocity
- final
- parameter
- isotropic
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 55
- 238000003384 imaging method Methods 0.000 claims abstract description 85
- 238000013508 migration Methods 0.000 claims abstract description 64
- 230000005012 migration Effects 0.000 claims abstract description 64
- 239000012504 chromatography matrix Substances 0.000 claims abstract description 33
- 238000004364 calculation method Methods 0.000 claims abstract description 17
- 238000012545 processing Methods 0.000 claims abstract description 13
- 238000004590 computer program Methods 0.000 claims description 26
- 230000006870 function Effects 0.000 claims description 11
- 238000003860 storage Methods 0.000 claims description 9
- 238000010276 construction Methods 0.000 claims description 8
- 239000011159 matrix material Substances 0.000 claims description 7
- 239000010410 layer Substances 0.000 description 56
- 238000010586 diagram Methods 0.000 description 16
- 230000008569 process Effects 0.000 description 12
- 238000004891 communication Methods 0.000 description 7
- 238000004587 chromatography analysis Methods 0.000 description 4
- 238000005553 drilling Methods 0.000 description 4
- 230000006872 improvement Effects 0.000 description 4
- 238000003325 tomography Methods 0.000 description 4
- 239000011229 interlayer Substances 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 238000012876 topography Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 2
- 238000007493 shaping process Methods 0.000 description 2
- 241000282414 Homo sapiens Species 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000005856 abnormality Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000003550 marker Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000001575 pathological effect Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
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/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- 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
-
- 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/34—Displaying seismic recordings or visualisation of seismic data or attributes
- G01V1/345—Visualisation of seismic data or attributes, e.g. in 3D 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/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
-
- 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/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/616—Data from specific type of measurement
- G01V2210/6161—Seismic or acoustic, e.g. land or sea measurements
-
- 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/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/616—Data from specific type of measurement
- G01V2210/6169—Data from specific type of measurement using well-logging
-
- 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/62—Physical property of subsurface
- G01V2210/626—Physical property of subsurface with anisotropy
-
- 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/64—Geostructures, e.g. in 3D data cubes
- G01V2210/643—Horizon tracking
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/70—Other details related to processing
- G01V2210/74—Visualisation of seismic data
-
- 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
技术领域
本发明涉及石油地震勘探技术领域,尤其涉及一种各向异性参数的求取方法及装置。
背景技术
本部分旨在为权利要求书中陈述的本发明实施例提供背景或上下文。此处的描述不因为包括在本部分中就承认是现有技术。
各向同性介质的常规地震成像技术和叠前时间/深度偏移地震成像技术已被广泛应用于构造解释和岩性识别,并获得了一定的应用效果,但其地震分辨率显然还难以满足薄储层识别和解决油田开发问题的需要。目前,地震勘探正在走向精细、超深层、高分辨率处理解释的方向,各向异性处理方法对于提高成像精度和消解井震误差显得尤为重要,对于提高钻井成功率也是必不可少的一步。
在各向异性层析反演方面,研究人员最初运用各向异性介质“程函方程”的哈密尔顿形式,推导出了各向异性介质走时扰动方程。后研究人员对走时扰动方程进行了修改,并且实现了三维各向异性介质走时层析成像方法,但是在实际资料中的应用效果不太好。还有研究人员在对VTI介质进行参数反演时,继而提出了对3个慢度分量(水平、垂直和动校正慢度)分别求取敏感核函数的方法。另有研究人员继而研究出了3D立体层析过程中数据空间与模型空间的线性关系,其主要原理是基于Hamiltonian算子的射线扰动理论。由层析网格剖分后建立的层析方程组,往往是病态或奇异的,属于典型的高度非线性反演问题。因此,初至波走时层析反演在本质上是不适定问题,需要通过正则化技术对反演过程施加约束,将先验信息融入反演方程组。后有研究人员通过正则化方法将先验信息引入初至波走时层析,取代了传统的外部约束模式。不过,随着地震勘探从二维到三维的发展和深入,地震层析反演方程维数成倍增加,传统基于Tikhonov正则化的反算法,存在收敛速度慢、迭代次数多、反演稳定性差等缺点。之后研究人员还提出了一种整形正则化方法,在迭代反演中通过整形算子将模型映射到平滑函数空间,显著提高了传统正则化方法的收敛速度与稳定性。
从上述各向异性层析反演研究里程来看,目前TTI各向异性参数建模与叠前深度偏移技术已逐渐成为了主流的处理手段,但传统的TTI各向异性参数的计算方法仍存在准确度低的问题,导致得到的地震偏移成像的精确度低。
发明内容
本发明实施例提供一种各向异性参数的求取方法,用以提升各向异性参数计算的准确度,得到精确度高的地震偏移成像剖面,该方法包括:
获取各向同性深度域速度体及测井分层数据;
根据各向同性深度域速度体确定解释层位体、各向同性倾角、各向同性方位角、各向异性速度体,将各向同性倾角和各向同性方位角转换为各向异性倾角和各向异性方位角;
对测井分层数据进行插值处理,得到测井层位体;
计算测井层位体与解释层位体在每一层位的深度差值,得到井震误差体;
根据井震误差体、解释层位体、各向异性速度体、各向异性倾角和各向异性方位角,构建层析矩阵;
利用层析矩阵迭代反演各向异性速度体和各向异性参数体,得到最终各向异性速度体、最终各向异性参数体,以及相应的最终地震偏移成像剖面。
本发明实施例还提供一种各向异性参数的求取装置,用以提升各向异性参数计算的准确度,得到精确度高的地震偏移成像剖面,该装置包括:
获取模块,用于获取各向同性深度域速度体及测井分层数据;
确定模块,用于根据各向同性深度域速度体确定解释层位体、各向同性倾角、各向同性方位角、各向异性速度体,将各向同性倾角和各向同性方位角转换为各向异性倾角和各向异性方位角;
插值模块,用于对测井分层数据进行插值处理,得到测井层位体;
计算模块,用于计算测井层位体与解释层位体在每一层位的层位差值,得到井震误差体;
矩阵构建模块,用于根据井震误差体、解释层位体、各向异性速度体、各向异性倾角和各向异性方位角,构建层析矩阵;
反演模块,用于利用层析矩阵迭代反演各向异性速度体和各向异性参数体,得到最终各向异性速度体、最终各向异性参数体,以及相应的最终地震偏移成像剖面。
本发明实施例还提供一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述各向异性参数的求取方法。
本发明实施例还提供一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,所述计算机程序被处理器执行时实现上述各向异性参数的求取方法。
本发明实施例还提供一种计算机程序产品,所述计算机程序产品包括计算机程序,所述计算机程序被处理器执行时实现上述各向异性参数的求取方法。
本发明实施例中,对传统的TTI各向异性参数计算方法进行了修改,将解释层位体加入各向异性层析反演当中,在没有井数据的大多区域内,以更合理的速度模型进行优化,从而能够得到更合理的各向异性参数体和速度体,进而得到更加合理的最终偏移成像剖面。此外,利用解释层位体和测井层位体计算了井震误差体,在进行层析反演时利用了该井震误差体,从而将地层信息考虑到反演过程中,在层析反演时加入了标志层的先验信息,充分考虑到标志层的重要性,获得更加合理的反演结果,再通过各向异性偏移,消除由于地层存在各向异性问题引起的井震误差,有效避免了常规通过单井点井震误差反演而产生的层位畸变问题,进一步提高了各向异性参数计算的精度。同时,利用井震误差体进行层析反演,很好的解决井震匹配中存在的误差的难题,有效提高了复杂地表地形条件的三维地震剖面的成像质量,保证了构造形态的合理性及可靠性。能够对断裂、逆掩推覆等复杂构造精确成像,为后续提高钻井成功率和储层预测提供技术保障。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。在附图中:
图1为本发明实施例中一种各向异性参数的求取方法的流程图;
图2为本发明实施例中一种各向同性深度域速度体的示意图;
图3为本发明实施例中一种基于各向同性叠前深度偏移剖面的解释层位体的示意图;
图4为本发明实施例中一种利用测井分层数据插值得到的测井层位体的示意图;
图5为本发明实施例中一种网格层析反演各向异性参数结果的平面示意图;
图6(a)为根据现有技术中各向同性叠前深度偏移剖面的示意图;
图6(b)为根据本发明实施例得到的各向异性叠前深度偏移剖面的示意图;
图7为本发明实施例中一种各向异性参数的求取装置的结构示意图;
图8为本发明实施例中一种计算机设备的结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白,下面结合附图对本发明实施例做进一步详细说明。在此,本发明的示意性实施例及其说明用于解释本发明,但并不作为对本发明的限定。
本发明实施例提供了一种各向异性参数的求取方法,如图1所示,该方法包括步骤101至步骤106:
步骤101、获取各向同性深度域速度体及测井分层数据。
各向同性深度域速度体是一个三维模型,在该各向同性深度域速度体的各个位置都具有对应的速度。示例性的,图2为得到的一种各向同性深度域速度体的示意图。
步骤102、根据各向同性深度域速度体确定解释层位体、各向同性倾角、各向同性方位角、各向异性速度体,将各向同性倾角和各向同性方位角转换为各向异性倾角和各向异性方位角。
在试验工区内,根据工区地质情况,在叠前时间偏移数据体上选取主要地质分层,并结合测井声波速度变化的界面进行地震解释层位拾取,得到解释层位体。这样获得的解释层位体通常对应地震强同相轴,信噪比较高,有利于后续走时残差的拾取。
示例性的,图3给出了一种基于各向同性深度域速度体得到的解释层位体示意图。
根据各向同性深度域速度体确定解释层位体等的确定,以及各向异性倾角和各向异性方位角的转换,均可以利用现有技术中方法进行,对于具体实现过程,在此不做赘述。
步骤103、对测井分层数据进行插值处理,得到测井层位体。
插值是一种在离散数据的基础上补插连续函数,使得连续曲线通过全部给定的离散数据点的方法,利用插值可通过函数在有限个点处的取值情况,估算出函数在其他点处的近似值。
示例性的,利用测井分层数据插值得到的测井层位体可以参照图4所示。
在具体实施时,可以先剔除插值异常值再进行插值,在必要时可以由人为对插值结果进行多次优化,以得到更符合地质认识的层位。
可以利用现有技术中方法通过测井分层数据插值得到测井层位体,具体实现过程,在此不做赘述。
步骤104、计算测井层位体与解释层位体在每一层位的深度差值,得到井震误差体。
井震误差仅存在于井位处,而层位信息是贯穿全工区的,所以通过引入层位相关的误差到后续各向异性层析反演当中,将层厚信息和井震误差融合考虑到,能得到避免反演出层厚异常的各向异性体。
以F0表示解释层位体中任一层位的速度值,F1表示测井层位体在任一层位的速度值,则按照Mistie=F1-F0计算得到各个层位上的层位误差值Mistie。将计算得到的所有层位的层位误差值确定为井震误差体。比如说,对于一个包含10个层位的解释层位体和测井层位体,先计算测井层位体第一个层位的速度与解释层位体第一个层位的深度差值,再计算测井层位体第二个层位的速度与解释层位体第二个层位的深度差值,逐个层位计算,直至得到全部10个层位的深度差值,将10个深度差值作为井震误差体。
步骤105、根据井震误差体、解释层位体、各向异性速度体、各向异性倾角和各向异性方位角,构建层析矩阵。
具体的,先将井震误差体、解释层位体、各向异性速度体、各向异性倾角和各向异性方位角作为输入,进行网格层析反演,得到层速度体模型;之后再将层速度体模型、解释层位体、井震误差体、各向异性倾角和各向异性方位角作为输入,按照给定的网格点进行射线追踪旅行时计算,建立层析矩阵。
相对于常规的各向异性参数计算方法,本发明实施例在构建层析矩阵时,将地层信息(包括解释层位体和测井层位体等)考虑到反演过程中,加入了井震误差体这一标志层的先验信息,作为各向异性参数反演的正则化条件,获得更加合理的反演结果,在之后的步骤中通过各向异性偏移,消除由于地层存在各向异性问题引起的井震误差,有效避免了常规通过单井点井震误差反演而产生的层位畸变问题,提高了各向异性参数的精度。
步骤106、利用层析矩阵迭代反演各向异性速度体和各向异性参数体,得到最终各向异性速度体、最终各向异性参数体,以及相应的最终地震偏移成像剖面。
具体的,在共成像点道集上进行自动拾取,获得深度域剩余延迟;以共成像点道集是否拉平作为误差泛函,利用深度域剩余延迟和层析矩阵,迭代反演各向异性速度体和各向异性参数体,得到最终各向异性速度体、最终各向异性参数体,以及相应的最终地震偏移成像剖面。这样一来,综合了层位控制和网格层析全局反演的优点,以地质层位为基础建立模型,在横向上沿地质层位网格和纵向上层间网格同时进行参数更新,能够在横向上沿合理的地质层位和纵向上反演优化的层间网格上双向同时向真实速度体逼近,并且将同相轴的连续性、方位角、地层倾角及测井数据等多种最终面属性信息联合应用于反演速度建模过程中,能够精确地得到速度场中高波数成分和低波数成分以及各向异性参数。
示例性的,图5给出了一种通过网格层析反演各向异性参数结果的平面示意图。
迭代反演能够得到多组对应的各向异性速度体、各向异性参数体,需要从中进行筛选得到符合精度要求的最终各向异性速度体和最终各向异性参数体,而各向异性速度体与各向异性参数体是否符合精度要求,则需要通过其得到的地震偏移成像剖面来体现。因此,以共成像点道集是否拉平作为误差泛函,利用深度域剩余延迟和层析矩阵,迭代反演各向异性速度体和各向异性参数,得到最终各向异性速度体、最终各向异性参数体,以及相应的最终地震偏移成像剖面,可以执行为如下方法:以共成像点道集是否拉平作为误差泛函,利用深度域剩余延迟求解层析矩阵,得到每个网格点的速度、各向异性参数和解释层位体的更新量,利用每个网格点的速度、各向异性参数和解释层位体的更新量反演得到中间各向异性速度体和中间各向异性参数体;利用中间各向异性速度体和中间各向异性参数体进行地震偏移成像,得到中间地震偏移成像剖面。
如果中间地震偏移成像剖面满足预设的误差要求,则将中间各向异性速度体、中间各向异性参数体和中间地震偏移成像剖面分别确定为最终各向异性速度体、最终各向异性参数体和最终地震偏移成像剖面。
如果中间地震偏移成像剖面不满足预设的误差要求,则重新求解层析矩阵并进行中间各向异性速度体和中间各向异性参数体的反演,直至得到的中间地震偏移成像剖面满足预设的误差要求,将满足预设的误差要求的中间各向异性速度体、中间各向异性参数体和中间地震偏移成像剖面分别确定为最终各向异性速度体、最终各向异性参数体和最终地震偏移成像剖面。
需要说明的是,本发明实施例中,中间各向异性速度体、中间各向异性参数体和中间地震偏移成像剖面为在反演过程中得到的一组各向异性速度体、各向异性参数体和地震偏移成像剖面的反演结果,最终各向异性速度体、最终各向异性参数体和最终地震偏移成像剖面为在反演过程中得到的中间各向异性速度体、中间各向异性参数体和中间地震偏移成像剖面在满足预设的误差要求之后,得到的最终的各项异性速度体、各项异性参数体和地震偏移成像剖面的反演结果。
图6(a)和图6(b)分别为根据现有技术中各向同性叠前深度偏移剖面的示意图以及根据本发明实施例得到的各向异性叠前深度偏移剖面的示意图,从两叠前深度偏移剖面对比来看,本发明实施例最终得到的剖面同相轴连续性和信噪比更高,断裂成像更清楚,目的层成像好,叠前深度偏移的深度不仅与老井深度一致,而且和新钻井深度一致,说明本发明求出的速度模型精度高,对构造位置的成像准确,很好地解决了该地区复杂构造成像问题。
本发明实施例中,对传统的TTI各向异性参数计算方法进行了修改,将解释层位体加入各向异性层析反演当中,在没有井数据的大多区域内,以更合理的速度模型进行优化,从而能够得到更合理的各向异性参数体和速度体,进而得到更加合理的最终偏移成像剖面。此外,利用解释层位体和测井层位体计算了井震误差体,在进行层析反演时利用了该井震误差体,从而将地层信息考虑到反演过程中,在层析反演时加入了标志层的先验信息,充分考虑到标志层的重要性,获得更加合理的反演结果,再通过各向异性偏移,消除由于地层存在各向异性问题引起的井震误差,有效避免了常规通过单井点井震误差反演而产生的层位畸变问题,进一步提高了各向异性参数计算的精度。同时,利用井震误差体进行层析反演,很好的解决井震匹配中存在的误差的难题,有效提高了复杂地表地形条件的三维地震剖面的成像质量,保证了构造形态的合理性及可靠性。能够对断裂、逆掩推覆等复杂构造精确成像,为后续提高钻井成功率和储层预测提供技术保障。
本发明实施例中还提供了一种各向异性参数的求取装置,如下面的实施例所述。由于该装置解决问题的原理与各向异性参数的求取方法相似,因此该装置的实施可以参见各向异性参数的求取方法的实施,重复之处不再赘述。
如图7所示,该装置700包括获取模块701、确定模块702、插值模块703、计算模块704、矩阵构建模块705和反演模块706。
其中,获取模块701,用于获取各向同性深度域速度体及测井分层数据;
确定模块702,用于根据各向同性深度域速度体确定解释层位体、各向同性倾角、各向同性方位角、各向异性速度体,将各向同性倾角和各向同性方位角转换为各向异性倾角和各向异性方位角;
插值模块703,用于对测井分层数据进行插值处理,得到测井层位体;
计算模块704,用于计算测井层位体与解释层位体在每一层位的深度差值,得到井震误差体;
矩阵构建模块705,用于根据井震误差体、解释层位体、各向异性速度体、各向异性倾角和各向异性方位角,构建层析矩阵;
反演模块706,用于利用层析矩阵迭代反演各向异性速度体和各向异性参数体,得到最终各向异性速度体、最终各向异性参数体,以及相应的最终地震偏移成像剖面。
在本发明实施例的一种实现方式中,矩阵构建模块705,用于:
将井震误差体、解释层位体、各向异性速度体、各向异性倾角和各向异性方位角作为输入,进行网格层析反演,得到层速度体模型;
将层速度体模型、解释层位体、井震误差体、各向异性倾角和各向异性方位角作为输入,按照给定的网格点进行射线追踪旅行时计算,建立层析矩阵。
在本发明实施例的一种实现方式中,反演模块706,用于:
在共成像点道集上进行自动拾取,获得深度域剩余延迟;
以共成像点道集是否拉平作为误差泛函,利用深度域剩余延迟和层析矩阵,迭代反演各向异性速度体和各向异性参数体,得到最终各向异性速度体、最终各向异性参数体,以及相应的最终地震偏移成像剖面。
在本发明实施例的一种实现方式中,反演模块706,用于:
以共成像点道集是否拉平作为误差泛函,利用深度域剩余延迟求解层析矩阵,得到每个网格点的速度、各向异性参数和解释层位体的更新量,利用每个网格点的速度、各向异性参数和解释层位体的更新量反演得到中间各向异性速度体和中间各向异性参数体;
利用中间各向异性速度体和中间各向异性参数体进行地震偏移成像,得到中间地震偏移成像剖面;
如果中间地震偏移成像剖面满足预设的误差要求,则将中间各向异性速度体、中间各向异性参数体和中间地震偏移成像剖面分别确定为最终各向异性速度体、最终各向异性参数体和最终地震偏移成像剖面;
如果中间地震偏移成像剖面不满足预设的误差要求,则重新求解层析矩阵并进行中间各向异性速度体和中间各向异性参数体的反演,直至得到的中间地震偏移成像剖面满足预设的误差要求,将满足预设的误差要求的中间各向异性速度体、中间各向异性参数体和中间地震偏移成像剖面分别确定为最终各向异性速度体、最终各向异性参数体和最终地震偏移成像剖面。
本发明实施例中,对传统的TTI各向异性参数计算方法进行了修改,将解释层位体加入各向异性层析反演当中,在没有井数据的大多区域内,以更合理的速度模型进行优化,从而能够得到更合理的各向异性参数体和速度体,进而得到更加合理的最终偏移成像剖面。此外,利用解释层位体和测井层位体计算了井震误差体,在进行层析反演时利用了该井震误差体,从而将地层信息考虑到反演过程中,在层析反演时加入了标志层的先验信息,充分考虑到标志层的重要性,获得更加合理的反演结果,再通过各向异性偏移,消除由于地层存在各向异性问题引起的井震误差,有效避免了常规通过单井点井震误差反演而产生的层位畸变问题,进一步提高了各向异性参数计算的精度。同时,利用井震误差体进行层析反演,很好的解决井震匹配中存在的误差的难题,有效提高了复杂地表地形条件的三维地震剖面的成像质量,保证了构造形态的合理性及可靠性。能够对断裂、逆掩推覆等复杂构造精确成像,为后续提高钻井成功率和储层预测提供技术保障。
本发明实施例还提供一种计算机设备,图8为本发明实施例中计算机设备的示意图,该计算机设备能够实现上述实施例中的各向异性参数的求取方法中全部步骤,该计算机设备具体包括如下内容:
处理器(processor)801、存储器(memory)802、通信接口(CommunicationsInterface)803和通信总线804;
其中,所述处理器801、存储器802、通信接口803通过所述通信总线804完成相互间的通信;所述通信接口803用于实现相关设备之间的信息传输;
所述处理器801用于调用所述存储器802中的计算机程序,所述处理器执行所述计算机程序时实现上述实施例中的各向异性参数的求取方法。
本发明实施例还提供一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,所述计算机程序被处理器执行时实现上述各向异性参数的求取方法。
本发明实施例还提供一种计算机程序产品,所述计算机程序产品包括计算机程序,所述计算机程序被处理器执行时实现上述各向异性参数的求取方法。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (11)
1.一种各向异性参数的求取方法,其特征在于,所述方法包括:
获取各向同性深度域速度体及测井分层数据;
根据各向同性深度域速度体确定解释层位体、各向同性倾角、各向同性方位角、各向异性速度体,将各向同性倾角和各向同性方位角转换为各向异性倾角和各向异性方位角;
对测井分层数据进行插值处理,得到测井层位体;
计算测井层位体与解释层位体在每一层位的深度差值,得到井震误差体;
根据井震误差体、解释层位体、各向异性速度体、各向异性倾角和各向异性方位角,构建层析矩阵;
利用层析矩阵迭代反演各向异性速度体和各向异性参数体,得到最终各向异性速度体和最终各向异性参数体,以及相应的最终地震偏移成像剖面。
2.根据权利要求1所述的方法,其特征在于,根据井震误差体、解释层位体、各向异性速度体、各向异性倾角和各向异性方位角,构建层析矩阵,包括:
将井震误差体、解释层位体、各向异性速度体、各向异性倾角和各向异性方位角作为输入,进行网格层析反演,得到层速度体模型;
将层速度体模型、解释层位体、井震误差体、各向异性倾角和各向异性方位角作为输入,按照给定的网格点进行射线追踪旅行时计算,建立层析矩阵。
3.根据权利要求1或2所述的方法,其特征在于,利用层析矩阵迭代反演各向异性速度体和各向异性参数体,得到最终各向异性速度体、最终各向异性参数体,以及相应的最终地震偏移成像剖面,包括:
在共成像点道集上进行自动拾取,获得深度域剩余延迟;
以共成像点道集是否拉平作为误差泛函,利用深度域剩余延迟和层析矩阵,迭代反演各向异性速度体和各向异性参数体,得到最终各向异性速度体、最终各向异性参数体,以及相应的最终地震偏移成像剖面。
4.根据权利要求3所述的方法,其特征在于,以共成像点道集是否拉平作为误差泛函,利用深度域剩余延迟和层析矩阵,迭代反演各向异性速度体和各向异性参数,得到最终各向异性速度体、最终各向异性参数体,以及相应的最终地震偏移成像剖面,包括:
以共成像点道集是否拉平作为误差泛函,利用深度域剩余延迟求解层析矩阵,得到每个网格点的速度、各向异性参数和解释层位体的更新量,利用每个网格点的速度、各向异性参数和解释层位体的更新量反演得到中间各向异性速度体和中间各向异性参数体;
利用中间各向异性速度体和中间各向异性参数体进行地震偏移成像,得到中间地震偏移成像剖面;
如果中间地震偏移成像剖面满足预设的误差要求,则将中间各向异性速度体、中间各向异性参数体和中间地震偏移成像剖面分别确定为最终各向异性速度体、最终各向异性参数体和最终地震偏移成像剖面;
如果中间地震偏移成像剖面不满足预设的误差要求,则重新求解层析矩阵并进行中间各向异性速度体和中间各向异性参数体的反演,直至得到的中间地震偏移成像剖面满足预设的误差要求,将满足预设的误差要求的中间各向异性速度体、中间各向异性参数体和中间地震偏移成像剖面分别确定为最终各向异性速度体、最终各向异性参数体和最终地震偏移成像剖面。
5.一种各向异性参数的求取装置,其特征在于,所述装置包括:
获取模块,用于获取各向同性深度域速度体及测井分层数据;
确定模块,用于根据各向同性深度域速度体确定解释层位体、各向同性倾角、各向同性方位角、各向异性速度体,将各向同性倾角和各向同性方位角转换为各向异性倾角和各向异性方位角;
插值模块,用于对测井分层数据进行插值处理,得到测井层位体;
计算模块,用于计算测井层位体与解释层位体在每一层位的深度差值,得到井震误差体;
矩阵构建模块,用于根据井震误差体、解释层位体、各向异性速度体、各向异性倾角和各向异性方位角,构建层析矩阵;
反演模块,用于利用层析矩阵迭代反演各向异性速度体和各向异性参数体,得到最终各向异性速度体、最终各向异性参数体,以及相应的最终地震偏移成像剖面。
6.根据权利要求5所述的装置,其特征在于,矩阵构建模块,用于:
将井震误差体、解释层位体、各向异性速度体、各向异性倾角和各向异性方位角作为输入,进行网格层析反演,得到层速度体模型;
将层速度体模型、解释层位体、井震误差体、各向异性倾角和各向异性方位角作为输入,按照给定的网格点进行射线追踪旅行时计算,建立层析矩阵。
7.根据权利要求5或6所述的装置,其特征在于,反演模块,用于:
在共成像点道集上进行自动拾取,获得深度域剩余延迟;
以共成像点道集是否拉平作为误差泛函,利用深度域剩余延迟和层析矩阵,迭代反演各向异性速度体和各向异性参数体,得到最终各向异性速度体、最终各向异性参数体,以及相应的最终地震偏移成像剖面。
8.根据权利要求7所述的装置,其特征在于,反演模块,用于:
以共成像点道集是否拉平作为误差泛函,利用深度域剩余延迟求解层析矩阵,得到每个网格点的速度、各向异性参数和解释层位体的更新量,利用每个网格点的速度、各向异性参数和解释层位体的更新量反演得到中间各向异性速度体和中间各向异性参数体;
利用中间各向异性速度体和中间各向异性参数体进行地震偏移成像,得到中间地震偏移成像剖面;
如果中间地震偏移成像剖面满足预设的误差要求,则将中间各向异性速度体、中间各向异性参数体和中间地震偏移成像剖面分别确定为最终各向异性速度体、最终各向异性参数体和最终地震偏移成像剖面;
如果中间地震偏移成像剖面不满足预设的误差要求,则重新求解层析矩阵并进行中间各向异性速度体和中间各向异性参数体的反演,直至得到的中间地震偏移成像剖面满足预设的误差要求,将满足预设的误差要求的中间各向异性速度体、中间各向异性参数体和中间地震偏移成像剖面分别确定为最终各向异性速度体、最终各向异性参数体和最终地震偏移成像剖面。
9.一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至4任一所述方法。
10.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有计算机程序,所述计算机程序被处理器执行时实现权利要求1至4任一所述方法。
11.一种计算机程序产品,其特征在于,所述计算机程序产品包括计算机程序,所述计算机程序被处理器执行时实现权利要求1至4任一所述方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111675949.8A CN116413798A (zh) | 2021-12-31 | 2021-12-31 | 各向异性参数的求取方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111675949.8A CN116413798A (zh) | 2021-12-31 | 2021-12-31 | 各向异性参数的求取方法及装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116413798A true CN116413798A (zh) | 2023-07-11 |
Family
ID=87055212
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111675949.8A Pending CN116413798A (zh) | 2021-12-31 | 2021-12-31 | 各向异性参数的求取方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116413798A (zh) |
-
2021
- 2021-12-31 CN CN202111675949.8A patent/CN116413798A/zh active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111077575B (zh) | 一种深度域速度建模方法与装置 | |
CN102841376A (zh) | 一种基于起伏地表的层析速度反演方法 | |
WO2011049885A2 (en) | Full-waveform inversion in the traveltime domain | |
EP3067718B1 (en) | Boundary layer tomography method and device | |
CA2876632C (en) | System and method for migration velocity modeling | |
CN113740901B (zh) | 基于复杂起伏地表的陆上地震数据全波形反演方法及装置 | |
CN107831542A (zh) | Ddw高精度深度域井震匹配方法 | |
CN107942379A (zh) | 一种提高复杂断块速度模型精度的方法 | |
CN109884700B (zh) | 多信息融合地震速度建模方法 | |
CN109655890B (zh) | 一种深度域浅中深层联合层析反演速度建模方法及系统 | |
CN104459798A (zh) | 一种基于rtm成像的速度建模方法 | |
US20200049844A1 (en) | Computer implemented method for improving a velocity model for seismic imaging | |
CN106842291B (zh) | 一种基于叠前地震射线阻抗反演的不整合圈闭储层岩性预测方法 | |
CN109188522B (zh) | 速度场构建方法及装置 | |
CN113945978B (zh) | 一种深度域地震数据优化方法及装置 | |
CN116413798A (zh) | 各向异性参数的求取方法及装置 | |
CN112649893B (zh) | 一种面向薄储层的多资料多参数融合建模方法及系统 | |
Montel et al. | Kinematics of common-image gathers—Part 2: Tomographic ray tracing and applications | |
CN111175825B (zh) | 深度域速度建模方法 | |
CN111310361A (zh) | 基于随钻地震的钻井指导方法、系统、设备及存储介质 | |
CN113433588B (zh) | 一种分偏移距扫描叠加的近地表速度分析方法 | |
CN112147700A (zh) | 速度异常区的低频模型构建方法及系统 | |
CN113267810B (zh) | 地震勘探全深度速度建模方法及装置 | |
CN113820745A (zh) | 地震速度建模方法、装置、电子设备及介质 | |
CN113960658B (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 |