CN114578427A - 井震联合全波形反演初始速度模型建立方法与装置 - Google Patents
井震联合全波形反演初始速度模型建立方法与装置 Download PDFInfo
- Publication number
- CN114578427A CN114578427A CN202011391197.8A CN202011391197A CN114578427A CN 114578427 A CN114578427 A CN 114578427A CN 202011391197 A CN202011391197 A CN 202011391197A CN 114578427 A CN114578427 A CN 114578427A
- Authority
- CN
- China
- Prior art keywords
- seismic
- velocity
- well
- data
- wells
- 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 57
- 239000011159 matrix material Substances 0.000 claims abstract description 26
- 230000009466 transformation Effects 0.000 claims abstract description 26
- 238000011160 research Methods 0.000 claims abstract description 12
- 238000001914 filtration Methods 0.000 claims description 20
- 238000003860 storage Methods 0.000 claims description 19
- 238000004590 computer program Methods 0.000 claims description 18
- 238000000605 extraction Methods 0.000 claims description 16
- 238000012937 correction Methods 0.000 claims description 6
- 238000005516 engineering process Methods 0.000 claims description 5
- 238000010586 diagram Methods 0.000 description 19
- 238000012545 processing Methods 0.000 description 9
- 230000006870 function Effects 0.000 description 8
- 230000008569 process Effects 0.000 description 8
- 238000004891 communication Methods 0.000 description 7
- 230000003287 optical effect Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000009499 grossing Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 238000004587 chromatography analysis Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 230000002969 morbid Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
- 238000012546 transfer 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/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
-
- 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
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
- G01V1/50—Analysing data
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供一种井震联合全波形反演初始速度模型建立方法与装置，该方法包括：获取研究区的地震数据以及声波测井数据；根据所述地震数据以及所述声波测井数据获取井间的地震速度、井点的地震速度以及对应的测井速度；利用所述井点的地震速度对应的测井速度比例校正所述井间的地震速度；根据变换矩阵与校正后的井间的地震速度构建初始速度模型，其中，通过变换矩阵使井间速度插值时能够保持地震速度的横向变化趋势，融合测井与地震速度信息，建立高精度的波形反演初始速度模型。
Description
技术领域
本发明涉及地震勘探技术领域，尤其涉及一种井震联合全波形反演初始速度模型建立方法与装置。
背景技术
全波形反演是地震勘探技术领域中的研究热点之一，该技术联合利用地震波场中的运动学和动力学信息，反演出高精度的地下介质速度、密度、各向异性模型，从而预测复杂地质背景下的构造、储层及流体信息。
全波形反演利用最优化算法迭代求解观测数据与波动方程正演数据的残差来求解地下介质参数。由于地下介质参数与地震波场间存在非线性关系，因此，该方法从数学上来讲是一个病态问题，存在多解性与收敛性问题。目前解决非线性反演的常用手段是利用Born近似将非线性问题转化为线性问题，然后利用局部优化方法迭代反演速度模型。在全波形反演理论中，要满足Born近似初始模型正演数据中波的走时与实际观测数据走时之间的误差在时间域要小于半个周期，在深度域小于半个波长。如果初始模型精度不够，反演目标函数会收敛到局部极小值，出现“周波跳跃”现象。因此，初始速度建模是全波形反演的一个核心问题。通常全波形反演初始速度模型使用射线层析、偏移速度分析或者声波测井插值得到。射线层析与偏移速度分析利用地震波的走时信息及同向轴曲率，得到的速度模型精度较低，只能反映相对速度变化。利用声波测井速度插值进行速度建模，通常使用层位控制，层内采用线性或双线性内插，无法反映速度横向变化。而且，利用上述方法建立速度模型后，需要按照全波形反演多尺度框架，进行整体空间平滑，这也会引入速度误差。
发明内容
针对现有技术中的问题，本发明提供一种井震联合全波形反演初始速度模型建立方法与装置、电子设备以及计算机可读存储介质，能够至少部分地解决现有技术中存在的问题。
为了实现上述目的，本发明采用如下技术方案：
第一方面，提供一种井震联合全波形反演初始速度模型建立方法，包括：
获取研究区的地震数据以及声波测井数据；
根据所述地震数据以及所述声波测井数据获取井间的地震速度、井点的地震速度以及对应的测井速度；
利用所述井点的地震速度对应的测井速度比例校正所述井间的地震速度；
根据变换矩阵与校正后的井间的地震速度构建初始速度模型。
进一步地，所述测井数据包括：井间距离，所述变换矩阵为：
其中，d为井间距离，g1为两口井之间的测井速度差，g2为井间差值速度，v1为第一口井对应位置的地震速度，V1为第一口井对应位置的测井速度。
进一步地，所述地震数据包括：地震速度场以及层位数据；
所述根据所述地震数据以及所述声波测井数据获取井间的地震速度、井点的地震速度以及对应的测井速度，包括：
对所述地震速度场进行中值滤波；
利用井震标定技术匹配所述地震数据和所述声波测井数据；
利用所述层位数据、所述声波测井数据中的井位置提取中值滤波后的地震速度场中的井间的地震速度以及井点的地震速度；
在井震匹配后的声波测井数据中提取与各地震速度对应的测井速度。
进一步地，所述利用所述层位数据、所述声波测井数据中的井位置提取中值滤波后的地震速度场中的井间的地震速度以及井点的地震速度，包括：
利用所述层位数据建立地层框架模型；
利用所述地层框架模型及所述声波测井数据中的井位置提取中值滤波后的地震速度场中的井间的地震速度以及井点的地震速度。
进一步地，所述根据变换矩阵与校正后的井间的地震速度构建初始速度模型，包括：
将所述变换矩阵与校正后的井间的地震速度相乘得到井间差值速度；
根据各个井间差值速度构建所述初始速度模型。
第二方面，提供一种井震联合全波形反演初始速度模型建立装置，包括：
数据获取模块，获取研究区的地震数据以及声波测井数据；
速度提取模块，根据所述地震数据以及所述声波测井数据获取井间的地震速度、井点的地震速度以及对应的测井速度；
速度校正模块，利用所述井点的地震速度对应的测井速度比例校正所述井间的地震速度；
建模模块，根据变换矩阵与校正后的井间的地震速度构建初始速度模型。
进一步地，所述测井数据包括：井间距离，所述变换矩阵为：
其中，d为井间距离，g1为两口井之间的测井速度差，g2为井间差值速度，v1为第一口井对应位置的地震速度，V1为第一口井对应位置的测井速度。
进一步地，所述地震数据包括：地震速度场以及层位数据；
所述速度提取模块包括：
中值滤波单元，对所述地震速度场进行中值滤波；
井震标定单元，利用井震标定技术匹配所述地震数据和所述声波测井数据；
地震速度提取单元，利用所述层位数据、所述声波测井数据中的井位置提取中值滤波后的地震速度场中的井间的地震速度以及井点的地震速度；
测井速度提取单元，在井震匹配后的声波测井数据中提取与各地震速度对应的测井速度。
第三方面，提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现上述的井震联合全波形反演初始速度模型建立方法的步骤。
第四方面，提供一种计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现上述的井震联合全波形反演初始速度模型建立方法的步骤。
本发明提供的井震联合全波形反演初始速度模型建立方法与装置，该方法包括：获取研究区的地震数据以及声波测井数据；根据所述地震数据以及所述声波测井数据获取井间的地震速度、井点的地震速度以及对应的测井速度；利用所述井点的地震速度对应的测井速度比例校正所述井间的地震速度；根据变换矩阵与校正后的井间的地震速度构建初始速度模型，其中，通过变换矩阵使井间速度插值时能够保持地震速度的横向变化趋势，融合测井与地震速度信息，建立高精度的波形反演初始速度模型。
为让本发明的上述和其他目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附图式，作详细说明如下。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。在附图中：
图1为本发明实施例中的服务器S1与客户端设备B1之间的架构示意图；
图2为本发明实施例中的服务器S1、客户端设备B1及数据库服务器S2之间的架构示意图；
图3是本发明实施例中的井震联合全波形反演初始速度模型建立方法的流程示意图；
图4示出了本发明实施例中校正前后井间速度极值示意图；
图5示出了本发明实施例中矩阵变换插值示意图；
图6示出了本发明实施例中步骤S200的具体步骤；
图7示出了本发明实施例中步骤S230的具体步骤；
图8示出了本发明实施例中井间地层框架示意图；
图9示出了本发明实施例中步骤S400的具体步骤；
图10是本发明实施例中的井震联合全波形反演初始速度模型建立装置的结构框图；
图11示出了本发明实施例中的速度提取模块的具体结构；
图12为本发明实施例电子设备的结构图。
具体实施方式
为了使本技术领域的人员更好地理解本申请方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分的实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护的范围。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
需要说明的是，本申请的说明书和权利要求书及上述附图中的术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本申请。
利用声波测井速度插值进行速度建模，通常使用层位控制，层内采用线性或双线性内插，无法反映速度横向变化。而且，利用上述方法建立速度模型后，需要按照全波形反演多尺度框架，进行整体空间平滑，这也会引入速度误差。
有鉴于此，本申请提供了一种井震联合全波形反演初始速度模型建立装置，该装置可以为一种服务器S1，参见图1，该服务器S1可以与至少一个客户端设备B1通信连接，所述客户端设备B1可以将研究区的地震数据以及声波测井数据发送至所述服务器S1，所述服务器S1可以在线接收所述研究区的地震数据以及声波测井数据。所述服务器S1可以在线或者离线对获取的研究区的地震数据以及声波测井数据进行预处理，根据所述地震数据以及所述声波测井数据获取井间的地震速度、井点的地震速度以及对应的测井速度；利用所述井点的地震速度对应的测井速度比例校正所述井间的地震速度；根据变换矩阵与校正后的井间的地震速度构建初始速度模型。而后，所述服务器S1可以将初始速度模型在线发送至所述客户端设备B1。所述客户端设备B1可以在线接收所述初始速度模型。
另外，参见图2，所述服务器S1还可以与至少一个数据库服务器S2通信连接，所述数据库服务器S2用于存储变换矩阵。所述数据库服务器S2在线将所述历史变换矩阵发送至所述服务器S1，所述服务器S1可以在线接收变换矩阵。
可以理解的是，所述客户端设备B1可以包括智能手机、平板电子设备、网络机顶盒、便携式计算机、台式电脑、个人数字助理(PDA)、车载设备、智能穿戴设备等。其中，所述智能穿戴设备可以包括智能眼镜、智能手表、智能手环等。
在实际应用中，进行井震联合全波形反演初始速度模型建立的部分可以在如上述内容所述的服务器S1侧执行，即，如图1所示的架构，也可以所有的操作都在所述客户端设备B1中完成，且该所述客户端设备B1可以直接与数据库服务器S2进行通信连接。具体可以根据所述客户端设备B1的处理能力，以及用户使用场景的限制等进行选择。本申请对此不作限定。若所有的操作都在所述客户端设备B1中完成，所述客户端设备B1还可以包括处理器，用于进行井震联合全波形反演初始速度模型建立的具体处理。
所述服务器与所述客户端设备之间可以使用任何合适的网络协议进行通信，包括在本申请提交日尚未开发出的网络协议。所述网络协议例如可以包括TCP/IP协议、UDP/IP协议、HTTP协议、HTTPS协议等。当然，所述网络协议例如还可以包括在上述协议之上使用的RPC协议(Remote Procedure Call Protocol，远程过程调用协议)、REST协议(Representational State Transfer，表述性状态转移协议)等。
图3是本发明实施例中的井震联合全波形反演初始速度模型建立方法的流程示意图；如图3所示，该井震联合全波形反演初始速度模型建立方法可以包括一下内容：
步骤S100：获取研究区的地震数据以及声波测井数据；
其中，测井数据包括：井间距离等。地震数据包括地震速度场以及层位数据等。
步骤S200：根据所述地震数据以及所述声波测井数据获取井间的地震速度、井点的地震速度以及对应的测井速度；
步骤S300：利用所述井点的地震速度对应的测井速度比例校正所述井间的地震速度；其中，图4示出了本发明实施例中校正前后井间速度极值示意图；
步骤S400：根据变换矩阵与校正后的井间的地震速度构建初始速度模型。
其中，变换矩阵为：
其中，d为井间距离，g1为两口井之间的测井速度差，g2为井间差值速度，v1为第一口井对应位置的地震速度，V1为第一口井对应位置的测井速度。
本发明实施实例适用于全波形反演的井震联合速度建模，提供了一个变换矩阵，将地震速度的横向变化趋势加入井间速度插值中，图5示出了本发明实施例中矩阵变换插值示意图，实现井震联合速度建模，模型测试结果表明了该方法的有效性。
在一个可选的实施例中，参见图6，该步骤S200可以包括以下内容：
步骤S210：对所述地震速度场进行中值滤波；
其中，通过中值滤波去除地震速度场中的奇异值。
步骤S220：利用井震标定技术匹配所述地震数据和所述声波测井数据；
步骤S230：利用所述层位数据、所述声波测井数据中的井位置提取中值滤波后的地震速度场中的井间的地震速度以及井点的地震速度；
步骤S240：在井震匹配后的声波测井数据中提取与各地震速度对应的测井速度。
在一个可选的实施例中，参加图7，该步骤S230可以包括以下内容：
步骤S231：利用所述层位数据建立地层框架模型。
其中，图8示出了本发明实施例中井间地层框架示意图。
步骤S232：利用所述地层框架模型及所述声波测井数据中的井位置提取中值滤波后的地震速度场中的井间的地震速度以及井点的地震速度。
其中，每一个地震速度有一个对应的测井速度。
在一个可选的实施例中，参见图9，该步骤S400可以包括以下内容：
步骤S410：将所述变换矩阵与校正后的井间的地震速度相乘得到井间差值速度；
步骤S420：根据各个井间差值速度构建所述初始速度模型。
基于同一发明构思，本申请实施例还提供了一种井震联合全波形反演初始速度模型建立装置，可以用于实现上述实施例所描述的方法，如下面的实施例所述。由于井震联合全波形反演初始速度模型建立装置解决问题的原理与上述方法相似，因此井震联合全波形反演初始速度模型建立装置的实施可以参见上述方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。
图10是本发明实施例中的井震联合全波形反演初始速度模型建立装置的结构框图。如图10所示，该井震联合全波形反演初始速度模型建立装置具体包括：数据获取模块10、速度提取模块20、速度校正模块30以及建模模块40。
数据获取模块10，获取研究区的地震数据以及声波测井数据；
速度提取模块20，根据所述地震数据以及所述声波测井数据获取井间的地震速度、井点的地震速度以及对应的测井速度；
速度校正模块30，利用所述井点的地震速度对应的测井速度比例校正所述井间的地震速度；
建模模块40，根据变换矩阵与校正后的井间的地震速度构建初始速度模型。
在一个可选的实施例中，参见图11，所述地震数据包括：地震速度场以及层位数据；
所述速度提取模块包括：中值滤波单元21、井震标定单元22、地震速度提取单元23以及测井速度提取单元24。
中值滤波单元21，对所述地震速度场进行中值滤波；
井震标定单元22，利用井震标定技术匹配所述地震数据和所述声波测井数据；
地震速度提取单元23，利用所述层位数据、所述声波测井数据中的井位置提取中值滤波后的地震速度场中的井间的地震速度以及井点的地震速度；
测井速度提取单元24，在井震匹配后的声波测井数据中提取与各地震速度对应的测井速度。
上述实施例阐明的装置、模块或单元，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。一种典型的实现设备为电子设备，具体的，电子设备例如可以为个人计算机、膝上型计算机、蜂窝电话、相机电话、智能电话、个人数字助理、媒体播放器、导航设备、电子邮件设备、游戏控制台、平板计算机、可穿戴设备或者这些设备中的任何设备的组合。
在一个典型的实例中电子设备具体包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现上述的井震联合全波形反演初始速度模型建立方法的步骤。
下面参考图12，其示出了适于用来实现本申请实施例的电子设备600的结构示意图。
如图12所示，电子设备600包括中央处理单元(CPU)601，其可以根据存储在只读存储器(ROM)602中的程序或者从存储部分608加载到随机访问存储器(RAM))603中的程序而执行各种适当的工作和处理。在RAM603中，还存储有系统600操作所需的各种程序和数据。CPU601、ROM602、以及RAM603通过总线604彼此相连。输入/输出(I/O)接口605也连接至总线604。
以下部件连接至I/O接口605：包括键盘、鼠标等的输入部分606；包括诸如阴极射线管(CRT)、液晶显示器(LCD)等以及扬声器等的输出部分607；包括硬盘等的存储部分608；以及包括诸如LAN卡，调制解调器等的网络接口卡的通信部分609。通信部分609经由诸如因特网的网络执行通信处理。驱动器610也根据需要连接至I/O接口605。可拆卸介质611，诸如磁盘、光盘、磁光盘、半导体存储器等等，根据需要安装在驱动器610上，以便于从其上读出的计算机程序根据需要被安装如存储部分608。
特别地，根据本发明的实施例，上文参考流程图描述的过程可以被实现为计算机软件程序。例如，本发明的实施例包括一种计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现上述的井震联合全波形反演初始速度模型建立方法的步骤。
在这样的实施例中，该计算机程序可以通过通信部分609从网络上被下载和安装，和/或从可拆卸介质611被安装。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括暂存电脑可读媒体(transitory media)，如调制的数据信号和载波。
为了描述的方便，描述以上装置时以功能分为各种单元分别描述。当然，在实施本申请时可以把各单元的功能在同一个或多个软件和/或硬件中实现。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、商品或者设备中还存在另外的相同要素。
本领域技术人员应明白，本申请的实施例可提供为方法、系统或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
以上所述仅为本申请的实施例而已，并不用于限制本申请。对于本领域技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本申请的权利要求范围之内。
Claims (10)
1.一种井震联合全波形反演初始速度模型建立方法，其特征在于，包括：
获取研究区的地震数据以及声波测井数据；
根据所述地震数据以及所述声波测井数据获取井间的地震速度、井点的地震速度以及对应的测井速度；
利用所述井点的地震速度对应的测井速度比例校正所述井间的地震速度；
根据变换矩阵与校正后的井间的地震速度构建初始速度模型。
2.根据权利要求1所述的井震联合全波形反演初始速度模型建立方法，其特征在于，所述测井数据包括：井间距离，所述变换矩阵为：
其中，d为井间距离，g1为两口井之间的测井速度差，g2为井间差值速度，v1为第一口井对应位置的地震速度，V1为第一口井对应位置的测井速度。
3.根据权利要求1所述的井震联合全波形反演初始速度模型建立方法，其特征在于，所述地震数据包括：地震速度场以及层位数据；
所述根据所述地震数据以及所述声波测井数据获取井间的地震速度、井点的地震速度以及对应的测井速度，包括：
对所述地震速度场进行中值滤波；
利用井震标定技术匹配所述地震数据和所述声波测井数据；
利用所述层位数据、所述声波测井数据中的井位置提取中值滤波后的地震速度场中的井间的地震速度以及井点的地震速度；
在井震匹配后的声波测井数据中提取与各地震速度对应的测井速度。
4.根据权利要求3所述的井震联合全波形反演初始速度模型建立方法，其特征在于，所述利用所述层位数据、所述声波测井数据中的井位置提取中值滤波后的地震速度场中的井间的地震速度以及井点的地震速度，包括：
利用所述层位数据建立地层框架模型；
利用所述地层框架模型及所述声波测井数据中的井位置提取中值滤波后的地震速度场中的井间的地震速度以及井点的地震速度。
5.根据权利要求1所述的井震联合全波形反演初始速度模型建立方法，其特征在于，所述根据变换矩阵与校正后的井间的地震速度构建初始速度模型，包括：
将所述变换矩阵与校正后的井间的地震速度相乘得到井间差值速度；
根据各个井间差值速度构建所述初始速度模型。
6.一种井震联合全波形反演初始速度模型建立装置，其特征在于，包括：
数据获取模块，获取研究区的地震数据以及声波测井数据；
速度提取模块，根据所述地震数据以及所述声波测井数据获取井间的地震速度、井点的地震速度以及对应的测井速度；
速度校正模块，利用所述井点的地震速度对应的测井速度比例校正所述井间的地震速度；
建模模块，根据变换矩阵与校正后的井间的地震速度构建初始速度模型。
7.根据权利要求6所述的井震联合全波形反演初始速度模型建立装置，其特征在于，所述测井数据包括：井间距离，所述变换矩阵为：
其中，d为井间距离，g1为两口井之间的测井速度差，g2为井间差值速度，v1为第一口井对应位置的地震速度，V1为第一口井对应位置的测井速度。
8.根据权利要求6所述的井震联合全波形反演初始速度模型建立装置，其特征在于，所述地震数据包括：地震速度场以及层位数据；
所述速度提取模块包括：
中值滤波单元，对所述地震速度场进行中值滤波；
井震标定单元，利用井震标定技术匹配所述地震数据和所述声波测井数据；
地震速度提取单元，利用所述层位数据、所述声波测井数据中的井位置提取中值滤波后的地震速度场中的井间的地震速度以及井点的地震速度；
测井速度提取单元，在井震匹配后的声波测井数据中提取与各地震速度对应的测井速度。
9.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现权利要求1至5任一项所述的井震联合全波形反演初始速度模型建立方法的步骤。
10.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该计算机程序被处理器执行时实现权利要求1至5任一项所述的井震联合全波形反演初始速度模型建立方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011391197.8A CN114578427A (zh) | 2020-12-02 | 2020-12-02 | 井震联合全波形反演初始速度模型建立方法与装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011391197.8A CN114578427A (zh) | 2020-12-02 | 2020-12-02 | 井震联合全波形反演初始速度模型建立方法与装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114578427A true CN114578427A (zh) | 2022-06-03 |
Family
ID=81766997
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011391197.8A Pending CN114578427A (zh) | 2020-12-02 | 2020-12-02 | 井震联合全波形反演初始速度模型建立方法与装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114578427A (zh) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2002073240A1 (en) * | 2001-03-13 | 2002-09-19 | Conoco Phillips Company | Method and process for prediction of subsurface fluid and rock pressures in the earth |
| CN102129084A (zh) * | 2010-12-17 | 2011-07-20 | 中国石油天然气股份有限公司 | 一种井控获取地震薄储层速度的方法及装置 |
| CN106324678A (zh) * | 2015-06-15 | 2017-01-11 | 中国石油化工股份有限公司 | 一种基于测井数据约束的全波形反演方法及系统 |
| CN106405651A (zh) * | 2016-11-14 | 2017-02-15 | 中国石油化工股份有限公司 | 一种基于测井匹配的全波形反演初始模型构建方法 |
| CN107783187A (zh) * | 2016-08-24 | 2018-03-09 | 中国石油化工股份有限公司 | 一种将测井速度和地震速度结合建立三维速度场的方法 |
| CN108663713A (zh) * | 2017-03-27 | 2018-10-16 | 中国石油化工股份有限公司 | 一种建立深度域构造模型的方法 |
| CN111077575A (zh) * | 2018-10-18 | 2020-04-28 | 中国石油化工股份有限公司 | 一种深度域速度建模方法与装置 |
-
2020
- 2020-12-02 CN CN202011391197.8A patent/CN114578427A/zh active Pending
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2002073240A1 (en) * | 2001-03-13 | 2002-09-19 | Conoco Phillips Company | Method and process for prediction of subsurface fluid and rock pressures in the earth |
| CN102129084A (zh) * | 2010-12-17 | 2011-07-20 | 中国石油天然气股份有限公司 | 一种井控获取地震薄储层速度的方法及装置 |
| CN106324678A (zh) * | 2015-06-15 | 2017-01-11 | 中国石油化工股份有限公司 | 一种基于测井数据约束的全波形反演方法及系统 |
| CN107783187A (zh) * | 2016-08-24 | 2018-03-09 | 中国石油化工股份有限公司 | 一种将测井速度和地震速度结合建立三维速度场的方法 |
| CN106405651A (zh) * | 2016-11-14 | 2017-02-15 | 中国石油化工股份有限公司 | 一种基于测井匹配的全波形反演初始模型构建方法 |
| CN108663713A (zh) * | 2017-03-27 | 2018-10-16 | 中国石油化工股份有限公司 | 一种建立深度域构造模型的方法 |
| CN111077575A (zh) * | 2018-10-18 | 2020-04-28 | 中国石油化工股份有限公司 | 一种深度域速度建模方法与装置 |
Non-Patent Citations (5)
| Title |
|---|
| 孙月成 等: "三维数据体时深转换速度模型的精细构建方法", 大庆石油地质与开发, vol. 37, no. 03, 30 June 2018 (2018-06-30), pages 145 - 152 * |
| 张珺: "加蓬盐下复杂构造区井控高精度变速成图的方法研究", 物探与化探, vol. 41, no. 03, 30 June 2017 (2017-06-30), pages 535 - 541 * |
| 张秉铭 等: "井震结合速度建模在水平井轨迹设计中的应用", 物探化探计算技术, vol. 39, no. 04, 31 July 2017 (2017-07-31), pages 552 - 558 * |
| 李红梅 等: "复杂构造二维地震精细速度分析与应用", 中国煤炭地质, vol. 31, no. 002, 28 February 2019 (2019-02-28), pages 1674 - 1803 * |
| 王明路: "声波测井速度与地震速度的关系研究", 中国优秀硕士学位论文全文数据库基础科学辑, no. 2, 15 February 2019 (2019-02-15), pages 011 - 253 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10951500B2 (en) | Method and apparatus for warning | |
| CN103415786A (zh) | 用于利用非线性模型更新来进行地震数据反演的系统和方法 | |
| CN111766635A (zh) | 砂体连通程度分析方法和系统 | |
| CN109597964B (zh) | 理论油藏渗流参数值确定方法、系统、设备、可读介质 | |
| CN114578427A (zh) | 井震联合全波形反演初始速度模型建立方法与装置 | |
| CN111983679A (zh) | 基于沉积参数的相控型白云岩储层地震预测方法和装置 | |
| CN111967584A (zh) | 生成对抗样本的方法、装置、电子设备及计算机存储介质 | |
| US11867857B2 (en) | Method and system for updating a seismic velocity model | |
| CN112081565B (zh) | 水驱油藏体积波及系数的确定方法及装置 | |
| CN115079257A (zh) | 基于融合网络的q值估计及地震衰减补偿方法 | |
| CN111596346B (zh) | 弹性波速度反演方法和装置 | |
| Kumar Sharma et al. | New exploration on approximate controllability of fractional neutral‐type delay stochastic differential inclusions with non‐instantaneous impulse | |
| CN112241022B (zh) | 基于射线密度生成层析反演模型速度界面的方法和装置 | |
| Li et al. | A numerical computation approach for the optimal control of ASP flooding based on adaptive strategies | |
| CN114076977B (zh) | 基于区块链的地震数据处理方法和装置 | |
| Pannekoucke et al. | PDE-NetGen 1.0: from symbolic PDE representations of physical processes to trainable neural network representations | |
| CN111983674A (zh) | 地层倾斜区的地层真实厚度确定方法和装置 | |
| CN114016998A (zh) | 一种测井遇阻遇卡的预测方法及装置 | |
| CN112859171B (zh) | 利用有效裂缝进行储层连通性分析的方法和装置 | |
| CN110632650B (zh) | 全波形反演数据的分时窗提取方法及装置 | |
| Das et al. | Minuscule Corrections to Near-surface Solar Internal Rotation Using Mode Coupling | |
| CN112052631B (zh) | 完钻新井目标井段初期无阻流量参数确定方法和装置 | |
| US20240069227A1 (en) | Method of stratigraphic calibration | |
| CN111526054B (zh) | 用于获取网络的方法及装置 | |
| CN113945971A (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 |